A method and system for determining shear wave splitting of converted waves, device, storage medium

Abstract

The application relates to a converted wave shear wave splitting calculation method and system, equipment and a storage medium, and belongs to the technical field of oil and natural gas geophysical exploration. The method comprises the following steps: respectively pre-processing, sorting and dynamic correction of converted wave radial component data and converted wave tangential component data; synchronously dynamic correcting radial component gather data and tangential component gather data; respectively performing azimuthal division on the obtained radial component data and tangential component data after synchronization; based on the azimuthal division obtained common receiver point radial component and tangential component azimuthal division gather, a phase inversion method and a phase-constrained least square method, the shear wave splitting is calculated according to a vertical layer. The application combines the advantages of the phase inversion method and the least square method, has lower requirements on input data, is suitable for the situation of low signal-to-noise ratio data and irregular acquisition data, and improves the shear wave splitting calculation precision.

Description

Technical Field

[0001] This invention belongs to the field of petroleum and natural gas geophysical exploration technology, and specifically relates to a method, system, equipment, and storage medium for obtaining converted wave shear wave splitting. Background Technology

[0002] Existing methods for shear wave splitting analysis include the phase inversion method and the least squares method, which are calculated independently. The phase inversion method requires regular azimuth angle data input, necessitating high-density data acquisition. The least squares method requires a high signal-to-noise ratio, and the quality of the original data acquisition directly affects its calculation accuracy. Therefore, existing technologies mainly suffer from the following problems and drawbacks: First, they place very high demands on the density of the original data acquisition. Second, when the data signal-to-noise ratio is low, the calculation accuracy is low, and it is easily affected by environmental noise. Summary of the Invention

[0003] To address the aforementioned problems, this invention provides a method, system, equipment, and storage medium for obtaining converted wave shear wave splitting. The method of this invention can be applied to the field of geophysical exploration, particularly for studying rock anisotropy. The method uses converted waves (P-to-S or S-to-P waves) in seismic wave propagation for analysis to obtain information about the physical properties of subsurface rocks. Specifically, it can be applied to the following sub-fields: a. Anisotropy Study: The converted wave shear wave splitting method is mainly used to study the anisotropic characteristics of subsurface rocks. This anisotropy can be caused by fractures, mineral textures, rock composition distribution, etc., and by analyzing the splitting phenomenon of shear waves, the degree and direction of anisotropy in subsurface rocks can be inferred. b. Oil and Gas Exploration: Studying the anisotropy of subsurface media such as fractures and pores helps to better understand the characteristics of oil and gas reservoirs and improve exploration accuracy. c. Converted Wave Imaging: This method can improve the accuracy and resolution of converted wave imaging, which is of great significance for subsequent joint interpretation of P-waves and S-waves.

[0004] The first objective of this invention is to provide a method for determining the splitting of converted wave shear waves, comprising:

[0005] Preprocessing, sorting, and dynamic correction are performed on the radial component data and the tangential component data of the converted wave, respectively.

[0006] Radial component gather data and tangential component gather data obtained by synchronous dynamic correction;

[0007] The radial and tangential component data obtained after synchronization are superimposed separately by azimuth angle;

[0008] Based on the azimuth superposition gather of radial and tangential components of the common receiver point obtained by azimuth superposition, and the phase reversal method and the phase-constrained least squares method, shear wave splitting is completed according to the vertical layer.

[0009] In a specific embodiment of the present invention, the preprocessing includes noise attenuation processing, static correction processing, and deconvolution processing.

[0010] In a specific embodiment of the present invention, the sorting includes:

[0011] The radial component data of the converted wave or the tangential component data of the converted wave are sorted into common receiver point gathers according to the trace header definition.

[0012] In a specific embodiment of the present invention, the formula for dynamic correction is as follows:

[0013] t = t p +t s

[0014] Where t is the conversion fluctuation correction time, t p For the longitudinal wave time, t s Transverse wave time.

[0015] In a specific embodiment of the present invention, the radial component gather data and tangential component gather data obtained by synchronous dynamic correction include:

[0016] In the common detector point domain, the corrected radial component data of the trace word and the corrected tangential component data of the trace word are marked with different labels.

[0017] The radial component data and tangential component data of the track header words with different labels are merged and processed.

[0018] Based on the shot point number and receiver point number information in the merged data, shot points and receiver points with the same radial and tangential components are selected, thus obtaining the synchronously processed radial and tangential component data.

[0019] In a specific embodiment of the present invention, the step of performing azimuth-based superposition of the radial component data and tangential component data obtained after synchronization includes:

[0020] The radial component gather data obtained after synchronization is divided into several parts according to the azimuth angle, and each part is superimposed according to the bisector angle to obtain the radial component azimuth superimposed gather of the common detector point.

[0021] The tangential component gather data obtained after synchronization is divided into several parts according to the azimuth angle. Each part is superimposed according to the bisector angle to obtain the tangential component azimuth superimposed gather of the common detector point.

[0022] The bisected angle is 360 degrees divided by the total number of azimuth angles.

[0023] In a specific embodiment of the present invention, the method of obtaining the shear wave splitting by vertical layer using the azimuth-based superposition of radial and tangential components of the common receiver point, along with the phase reversal method and the phase-constrained least squares method, includes:

[0024] Estimating the direction of fast wave: On the common detector point superposition gather of the tangential component, the phase reversal method is used to scan along layer H1 to obtain the angle at which the phase reversal occurs, which is the estimated direction of fast wave.

[0025] Calculate the fast wave direction: Using the estimated fast wave direction as the center, perform least squares fitting on the radial component gather data to calculate the fast wave direction;

[0026] Calculate the time delay attributes of fast and slow waves: Based on the calculated direction of the fast wave, the cross-correlation method is used to estimate the time delay of the fast and slow waves, and obtain the time delay attributes of layer H1, that is, the time delay attributes of fast and slow waves.

[0027] One rotation: Based on the calculated fast wave direction, rotate the radial component gather to the fast wave S1 direction and the slow wave S2 direction, and at the same time rotate the tangential component gather to the fast wave S1 direction and the slow wave S2 direction, thus obtaining the fast wave S1 gather and the slow wave S2 gather.

[0028] Correction: Based on the time delay properties of the fast and slow waves, and in order to align with the direction of the fast wave S1, the slow wave S2 gather is time-corrected to obtain the corrected slow wave S2′ gather.

[0029] Secondary rotation: Rotate the fast wave S1 gather and the slow wave S2′ gather to their initial radial and tangential directions respectively, thus obtaining the gather data and gather data after shear wave splitting compensation.

[0030] Repeat: Repeat the steps of "estimate fast wave direction", "calculate fast wave direction", "calculate fast and slow wave delay attributes", "first rotation", "correction" and "second rotation" to perform shear wave splitting calculation and compensation for the next layer after layer H1 until all in-phase axes on the tangential component data have completed shear wave splitting calculation and compensation, that is, the shear wave splitting is obtained.

[0031] A second objective of this invention is to provide a phase-reversal constrained system for determining the splitting of converted wave transverse waves, comprising:

[0032] Preprocessing module: used to preprocess, sort, and dynamically correct the radial component data and tangential component data of the converted wave, respectively;

[0033] Synchronization module: used to synchronize the radial component gather data and tangential component gather data obtained from dynamic correction;

[0034] The superposition module is used to superimpose the radial component data and tangential component data obtained after synchronization by azimuth angle.

[0035] The acquisition module is used to obtain shear wave splitting based on the azimuth superposition gather of radial and tangential components of the common receiver point obtained by azimuth superposition, as well as the phase reversal method and the phase-constrained least squares method, according to the vertical layer.

[0036] In a specific embodiment of the present invention, the obtaining module includes an estimation submodule, a calculation submodule, a rotation submodule, a correction submodule, and a repetition submodule:

[0037] Estimation submodule: Used to estimate the fast wave direction: On the common detector point superposition gather of the tangential component, the phase reversal method is used to scan along layer H1 to obtain the angle at which the phase reversal occurs, which is the estimated fast wave direction.

[0038] The calculation submodule is used to calculate the fast wave direction and the fast and slow wave time delay attributes. Specifically, the fast wave direction is calculated by performing least squares fitting on the radial component gather data with the estimated fast wave direction as the center. The fast and slow wave time delay attributes are calculated by estimating the time delay of the fast and slow waves using the cross-correlation method based on the calculated fast wave direction, thereby obtaining the time delay attributes of layer H1, i.e., the fast and slow wave time delay attributes.

[0039] The rotation submodule is used to perform a first rotation and a second rotation. The first rotation involves rotating the radial component gather to the fast wave S1 direction and the slow wave S2 direction according to the calculated fast wave direction, and simultaneously rotating the tangential component gather to the fast wave S1 direction and the slow wave S2 direction, thus obtaining the fast wave S1 gather and the slow wave S2 gather. The second rotation involves rotating the fast wave S1 gather and the slow wave S2' gather to the initial radial and tangential directions, respectively, thus obtaining the gather data after shear wave splitting compensation.

[0040] The correction submodule is used to perform correction: based on the time delay properties of the fast and slow waves, and in order to align with the direction of the fast wave S1, the slow wave S2 gather is time-corrected to obtain the corrected slow wave S2' gather.

[0041] The repeat module is used to repeat the steps of "estimating the fast wave direction", "calculating the fast wave direction", "calculating the fast and slow wave delay attributes", "first rotation", "correction" and "second rotation" to perform shear wave splitting calculation and compensation for the next layer after layer H1, until all in-phase axes on the tangential component data have completed shear wave splitting calculation and compensation, that is, the shear wave splitting is obtained.

[0042] A third object of the present invention is to provide an electronic device comprising: a processor coupled to a memory;

[0043] The memory is used to store computer programs;

[0044] The processor is configured to execute the computer program stored in the memory, so that the electronic device performs the method described above.

[0045] A fourth object of the present invention is to provide a computer-readable storage medium storing a program or instructions that, when executed on a computer, cause the computer to perform the method described above.

[0046] The beneficial effects of this invention are:

[0047] This invention discloses a method, system, device, and storage medium for determining transverse wave splitting. The method combines the advantages of the phase reversal method and the least squares method, has lower requirements for input data, and is suitable for low signal-to-noise ratio data and irregularly acquired data. It can solve the problem that currently acquired three-component seismic data is limited by insufficient sampling density and low data signal-to-noise ratio, making it difficult to accurately determine the azimuth and magnitude of transverse wave splitting in anisotropic media, thus improving the accuracy of transverse wave splitting determination.

[0048] The method for obtaining converted wave shear wave splitting in this invention has broad application prospects in the fields of earth science and exploration, and is of great significance for deepening the understanding of underground structures and properties and solving resource exploration and environmental problems.

[0049] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures pointed out in the description, claims and drawings. Attached Figure Description

[0050] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0051] Figure 1 A flowchart of a method for determining the splitting of a converted wave according to an embodiment of the present invention is shown;

[0052] Figure 2 A schematic diagram illustrating a first rotation, correction, and second rotation according to an embodiment of the present invention is shown;

[0053] Figure 3The comparison results of the azimuth angle superimposed trace synthesis record and the actual data in Example 1 of the present invention are shown;

[0054] Figure 4 The R-component and T-component gathers are shown using the method described in Example 1 of the present invention.

[0055] Figure 5 The R-component and T-component gathers are shown after using the acquisition method according to Example 1 of the present invention;

[0056] Figure 6 A schematic diagram of the preconversion wave on the common conversion point (ACP) superposition profile using the method of obtaining the wave according to Example 1 of the present invention is shown.

[0057] Figure 7 A schematic diagram of the converted wave on the common conversion point (ACP) superposition profile after using the method of obtaining it according to Example 1 of the present invention is shown;

[0058] Figure 8 The R and T components before and after using the acquisition method according to Example 2 of the present invention are shown;

[0059] Figure 9 The superimposed profile of the R components before using the determination method according to Example 2 of the present invention is shown;

[0060] Figure 10 The superimposed profile of the R components after using the method of determination according to Example 2 of the present invention is shown;

[0061] Figure 11 This illustrates the fast and slow shear wave gathers of the well transition wave used in Example 3 of the present invention for gas content detection;

[0062] Figure 12 A framework diagram of a converted wave shear wave splitting determination system according to an embodiment of the present invention is shown;

[0063] Figure 13 A frame diagram of an electronic device according to an embodiment of the present invention is shown;

[0064] In the figure: preprocessing module 10; synchronization module 20; superposition module 30; acquisition module 40; electronic device 300; processor 301; memory 302. Detailed Implementation

[0065] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0066] like Figure 1 As shown, a method for determining the splitting of a converted wave according to an embodiment of the present invention includes:

[0067] X1. Preprocess, sort, and dynamically correct the radial component data and tangential component data of the converted wave, respectively.

[0068] X2, radial component gather data and tangential component gather data obtained by synchronous dynamic correction;

[0069] X3. Perform azimuth-based superposition on the radial and tangential component data obtained after synchronization.

[0070] X4. Based on the azimuth superposition gather of radial and tangential components of the common receiver point obtained by azimuth superposition, and the phase reversal method and the least squares method with phase constraint, the shear wave splitting is completed according to the vertical layer.

[0071] In this embodiment of the invention, the vertical layer refers to a position deeper underground from the ground (vertical coordinate is zero), where the converted wave is divided into some distinct interfaces, which are identified based on the phase axis, i.e., the layer position of the vertical layer.

[0072] In some embodiments of the present invention, the preprocessing in step X1 includes noise attenuation processing, static correction processing, and deconvolution processing. Noise attenuation processing, static correction processing, and deconvolution processing are all common operations in the technical field, and will not be described in detail here.

[0073] In some embodiments of the present invention, step X1, the sorting includes: sorting the converted wave radial component data or the converted wave tangential component data into common detector point gathers according to the trace head definition.

[0074] In some embodiments of the present invention, in step X1, the formula for the dynamic correction is as shown in equation (1):

[0075] t = t p +t s (1)

[0076] In equation (1), t is the conversion fluctuation correction time. p For the longitudinal wave time, t s Transverse wave time.

[0077] In some embodiments of the present invention, step X2 includes:

[0078] X2-1. In the common detector point domain, the corrected radial component data and the corrected tangential component data of the trace word are marked with different labels, that is, the corrected radial component data of the trace word is marked as a, and the corrected tangential component data of the trace word is marked as b.

[0079] X2-2. Merge the radial component data and tangential component data of the track header with different labels, that is, combine the two data into one data;

[0080] X2-3. Based on the shot point number and receiver point number information in the merged data, shot points and receiver points with the same radial and tangential components are selected. That is, based on the shot point number and receiver point number information in the data merged in step X2-2, the shot point number and receiver point number of the aforementioned data are selected. Due to the different labels marked in step X2-1, the radial component gather and tangential component gather are obtained from the shot point number and receiver point number of the aforementioned data, respectively. Thus, the synchronously processed radial component data and the synchronously processed tangential component data are obtained.

[0081] In some embodiments of the present invention, step X3 includes:

[0082] X3-1. Divide the radial component gather data obtained after synchronization into several parts according to the azimuth angle, and superimpose each part according to the bisector angle to obtain the radial component azimuth superimposed gather of the common detector point.

[0083] X3-2. Divide the tangential component gather data obtained after synchronization into several parts according to the azimuth angle, and superimpose each part according to the bisector angle to obtain the tangential component azimuth superimposed gather of the common detector point.

[0084] Wherein, the bisected angle is 360 degrees divided by the total number of azimuth angles;

[0085] In this embodiment of the invention, the division of the azimuth angle into several parts is not specifically limited, and can be flexibly set according to the amount of calculation and calculation efficiency in specific implementation.

[0086] In some embodiments of the present invention, step X4 includes:

[0087] X4-1. Estimating the direction of fast wave: On the common detector point superposition gather of the tangential component, the phase reversal method is used to scan along layer H1 to obtain the angle at which the phase reversal occurs, which is the estimated direction of fast wave.

[0088] X4-2. Calculate the fast wave direction: Using the estimated fast wave direction as the center, perform least squares fitting on the radial component gather data to calculate the fast wave direction.

[0089] X4-3. Calculate the time delay attributes of fast and slow waves: Based on the calculated fast wave direction, the cross-correlation method is used to estimate the time delay of fast and slow waves, and obtain the time delay attributes of layer H1, that is, the time delay attributes of fast and slow waves.

[0090] X4-3. First rotation: Based on the calculated fast wave direction, rotate the radial component gather to the fast wave S1 direction and the slow wave S2 direction, and at the same time rotate the tangential component gather to the fast wave S1 direction and the slow wave S2 direction, thus obtaining the fast wave S1 gather and the slow wave S2 gather.

[0091] X4-4. Correction: Based on the time delay properties of the fast and slow waves, and in order to align with the direction of the fast wave S1, the slow wave S2 gather is time-corrected to obtain the corrected slow wave S2' gather.

[0092] X4-5, Secondary rotation: Rotate the fast wave S1 gather and the slow wave S2′ gather to the initial radial and tangential directions respectively to obtain the gather data and gather data after shear wave splitting compensation.

[0093] X4-6. Repeat: Repeat the steps of "estimate fast wave direction", "calculate fast wave direction", "calculate fast and slow wave delay attributes", "first rotation", "correction" and "second rotation" to perform shear wave splitting calculation and compensation for the next layer below layer H1 (that is, in the vertical layer direction, perform shear wave splitting calculation and compensation in the order from top to bottom) until all in-phase axes on the tangential component data have completed shear wave splitting calculation and compensation (that is, there are no obvious converted wave in-phase axes on the tangential data), and the shear wave splitting is completed.

[0094] For detailed diagrams of the first rotation, correction, and second rotation in steps X4-3, X4-4, and X4-5 above, please refer to [link / reference]. Figure 2 .

[0095] To facilitate the explanation of the converted wave shear wave splitting method provided in the embodiments of the present invention, the specific process of the method is illustrated by the following example.

[0096] Example 1

[0097] Phase 1

[0098] Step 1:

[0099] Select the radial component (R component) and tangential component (T component) data of the converted wave, extract the data into the common detector point domain, and perform converted wave correction processing respectively.

[0100] In this step, the specific method is as follows: after completing the noise attenuation, static correction and deconvolution processing of the converted wave radial component data or tangential component data, the data is sorted into common detector point gathers according to the trace head word definition, and then converted wave correction is performed. The dynamic correction satisfies equation (1).

[0101] Step Two:

[0102] Synchronize radial component gather data and tangential component gather data.

[0103] The procedure for this step is as follows:

[0104] Ⅰ. In the common detector point domain, set the radial component data of the trace word to 1 and the tangential component data of the trace word to 2;

[0105] II. Merge these two sets of data;

[0106] III. Based on the shot point number and receiver point number information in the data, shot points and receiver points with the same radial and tangential components are selected. Then, according to the information of trace word 1 and trace word 2, the radial component data and tangential component data after synchronization processing are obtained, thus realizing the synchronization processing of data.

[0107] Step 3:

[0108] The radial and tangential component data after synchronization are superimposed separately by azimuth angle.

[0109] The specific method is as follows: The radial and tangential components of the common receiver point gather data are divided into 36 parts according to azimuth angle, and each part is superimposed at 10° intervals to obtain the superimposed azimuth gather of the radial and tangential components of the common receiver point. (Refer to...) Figure 3 , Figure 3 This is a comparison chart of the results of the overlay of azimuth angles and the actual data.

[0110] Step Four:

[0111] On the azimuth stacked gather, the fast wave direction is calculated using the phase reversal method by using the T-component azimuth stacked data, and the fast wave direction information θ is obtained.

[0112] The operation method of this step is as follows: on the common detector point superposition gather of the T component, perform a phase reversal scan along layer H1 to determine the angle at which the phase reversal occurs, i.e., the fast wave direction θ.

[0113] Step 5: Based on the fast wave direction information θ, use the least squares method to calculate the fast wave direction θ' and the time attributes of fast and slow wave delays.

[0114] The operation method of this step is as follows: with the fast wave direction θ as the center, perform least squares fitting on the R component gather data to obtain a more accurate fast wave direction θ'. After obtaining the fast wave direction θ', use the cross-correlation method to estimate the time delay of the fast wave and the slow wave, and obtain the time delay attribute delta of H1 at this position.

[0115] Phase Two

[0116] Step 6: Rotate the R component gather and T component gather to the fast wave S1 direction and the slow wave S2 direction respectively to obtain the fast wave S1 gather and the slow wave S2 gather.

[0117] The operation method for this step is as follows: rotate the R component gather data and the T component gather data to the fast wave direction θ' to obtain the fast wave S1 gather and the slow wave S2 gather.

[0118] Step 7: Perform time correction on the slow wave S2 gather data according to the delay time attribute to obtain the corrected slow wave S2' gather data.

[0119] The operation method of this step is as follows: Based on the slow wave delay time data obtained in step five, the slow wave S2 gather data is adjusted for delay time to align it with the fast wave S1 direction, that is, to make it consistent with the fast wave data, thus obtaining the slow wave S2' gather data.

[0120] Step 8: Rotate the fast wave S1 gather and the slow wave S2' gather to their original radial direction R and tangential direction T to obtain the R' and T' gather data after shear wave splitting compensation.

[0121] The procedure for this step is as follows: rotate the fast wave S1 gather and slow wave S2' gather data back to their original R and T directions to obtain the R' and T' component data after shear wave splitting correction.

[0122] Step 9: Repeat steps 4 to 8 to perform the second layer of shear wave splitting and compensation calculations to obtain the final shear wave splitting and compensation results.

[0123] The operation method for this step is as follows: repeat steps four to eight to perform the second layer of shear wave splitting calculation and compensation, and so on, until there is no obvious converted wave phase axis in the tangential data, and the shear wave splitting compensation result is obtained.

[0124] refer to Figure 4-7 , Figure 4 , 5 This is a diagram showing the effect of applying it to a common detector point gather. Figure 6 , 7 This is an image showing the effect of applying it to an ACP overlay profile. Figure 4 This is the common receiver gather of the R and T components before applying this method. Figure 5This is a common receiver gather for the R and T components after applying this method. Figure 6 This is a superimposed profile of the common conversion points (ACP) of the R components before applying this method; Figure 7 This is a superimposed profile of the common conversion points (ACP) of the R components after applying this method.

[0125] Through the above steps, higher-precision calculation and application of converted wave shear wave splitting properties were achieved.

[0126] Example 2

[0127] The converted wave shear wave was obtained by performing steps X1-X4, and the results are as follows: Figures 8-10 .

[0128] Figure 8 The data shows that applying the converted wave shear wave splitting method of this invention improves the quality of the converted wave gather:

[0129] refer to Figure 8 The present invention compares the R-component and T-component gathers of the converted wave before and after converting wave shear wave splitting compensation. After the R-component fast and slow wave time difference correction, the gathers are straighter, which is beneficial to improve in-phase superposition and enhance the quality of profile imaging; the T-component energy is weakened.

[0130] Figures 9-10 The data shows that applying the phase-converted wave shear wave splitting method of this invention improves the quality of converted wave profile imaging:

[0131] refer to Figure 9 and 10 , Figure 9 This is the superimposed profile of the converted wave R components before converted wave shear wave splitting compensation using this invention. Figure 10 The converted wave R component superposition profile after applying the present invention for converted wave shear wave splitting compensation has better phase axis focusing and higher longitudinal resolution.

[0132] Example 3

[0133] The converted wave shear wave was obtained by performing steps X1-X4, and the results are as follows: Figure 11 .

[0134] Figure 11 The data in the figure demonstrates a converted wave shear wave splitting method according to an embodiment of the present invention, used for gas content detection in target strata:

[0135] refer to Figure 11The converted wave fast and slow shear wave gathers obtained using this invention are used to collect the wellbore pass-point converted wave data. In the wellbore pass-point data, wells 1, 2, and 3 contain gas, and the converted wave exhibits shear wave splitting, with a significant difference in the conversion wave fast and slow shear wave time. Well 4 does not contain gas, and the converted wave shear wave splitting effect is very weak, with the conversion wave fast and slow shear wave time difference almost zero. By using the converted wave fast and slow shear wave gathers obtained using this invention and comparing the time differences at target formations, gas-bearing detection results can be obtained as a reference.

[0136] like Figure 12 As shown, a phase-reversal constrained converted wave shear wave splitting determination system according to the present invention includes:

[0137] Preprocessing module 10: used to preprocess, sort, and dynamically correct the radial component data and tangential component data of the converted wave, respectively;

[0138] Synchronization module 20: used to synchronize the radial component gather data and tangential component gather data obtained from dynamic correction;

[0139] Overlay module 30: used to overlay the radial component data and tangential component data obtained after synchronization by azimuth angle respectively;

[0140] Module 40: Used to obtain shear wave splitting based on the radial and tangential components of the common receiver point azimuth superposition gather obtained by azimuth superposition, as well as the phase reversal method and the phase-constrained least squares method, according to the vertical layer.

[0141] In some embodiments of the present invention, the obtaining module 40 includes an estimation submodule, a calculation submodule, a rotation submodule, a correction submodule, and a repetition submodule:

[0142] Estimation submodule: Used to estimate the fast wave direction: On the common detector point superposition gather of the tangential component, the phase reversal method is used to scan along layer H1 to obtain the angle at which the phase reversal occurs, which is the estimated fast wave direction.

[0143] The calculation submodule is used to calculate the fast wave direction and the fast and slow wave time delay attributes. Specifically, the fast wave direction is calculated by performing least squares fitting on the radial component gather data with the estimated fast wave direction as the center. The fast and slow wave time delay attributes are calculated by estimating the time delay of the fast and slow waves using the cross-correlation method based on the calculated fast wave direction, thereby obtaining the time delay attributes of layer H1, i.e., the fast and slow wave time delay attributes.

[0144] The rotation submodule is used to perform a first rotation and a second rotation. The first rotation involves rotating the radial component gather to the fast wave S1 direction and the slow wave S2 direction according to the calculated fast wave direction, and simultaneously rotating the tangential component gather to the fast wave S1 direction and the slow wave S2 direction, thus obtaining the fast wave S1 gather and the slow wave S2 gather. The second rotation involves rotating the fast wave S1 gather and the slow wave S2' gather to the initial radial and tangential directions, respectively, thus obtaining the gather data after shear wave splitting compensation.

[0145] The correction submodule is used to perform correction: based on the time delay properties of the fast and slow waves, and in order to align with the direction of the fast wave S1, the slow wave S2 gather is time-corrected to obtain the corrected slow wave S2' gather.

[0146] The repeat module is used to repeat the steps of "estimating the fast wave direction", "calculating the fast wave direction", "calculating the fast and slow wave delay attributes", "first rotation", "correction" and "second rotation" to perform shear wave splitting calculation and compensation for the next layer after layer H1, until all in-phase axes on the tangential component data have completed shear wave splitting calculation and compensation, that is, the shear wave splitting is obtained.

[0147] like Figure 13 As shown, in some embodiments of the present invention, an electronic device is provided, the electronic device 300 including: a processor 301 coupled to a memory 302;

[0148] The memory 302 is used to store computer programs;

[0149] The processor 301 is configured to execute the computer program stored in the memory 302, so that the electronic device performs the method described in the above embodiments.

[0150] In some embodiments of the present invention, a computer-readable storage medium is provided that stores a program or instructions that, when executed on a computer, cause the computer to perform the methods described in the above embodiments.

[0151] According to embodiments of the present invention, the computer-readable storage medium may be a non-volatile computer-readable storage medium, such as including, but not limited to: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof. In the present invention, the computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, electronic device, or apparatus.

[0152] Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for determining the splitting of a converted transverse wave, characterized in that, include: Preprocessing, sorting, and dynamic correction are performed on the radial component data and the tangential component data of the converted wave, respectively. Radial component gather data and tangential component gather data obtained by synchronous dynamic correction; The radial and tangential component data obtained after synchronization are superimposed separately by azimuth angle; Based on the azimuth superposition gather of radial and tangential components of the common receiver point obtained by azimuth superposition, and the phase reversal method and the phase-constrained least squares method, shear wave splitting is completed according to the vertical layer.

2. The method for determining the splitting of a converted wave shear wave according to claim 1, characterized in that, The preprocessing includes noise attenuation processing, static correction processing, and deconvolution processing.

3. The method for determining the splitting of a converted wave shear wave according to claim 1, characterized in that, The sorting includes: The radial component data of the converted wave or the tangential component data of the converted wave are sorted into common receiver point gathers according to the trace header definition.

4. The method for determining the splitting of a converted wave shear wave according to claim 1, characterized in that, The formula for dynamic correction is as follows: t=t p +t s Where t is the conversion fluctuation correction time, t p For the longitudinal wave time, t s Transverse wave time.

5. The method for determining the splitting of a converted wave shear wave according to claim 1, characterized in that, The radial component gather data and tangential component gather data obtained by synchronous dynamic correction include: In the common detector point domain, the corrected radial component data of the trace word and the corrected tangential component data of the trace word are marked with different labels. The radial component data and tangential component data of the track header words with different labels are merged and processed. Based on the shot point number and receiver point number information in the merged data, shot points and receiver points with the same radial and tangential components are selected, thus obtaining the synchronously processed radial and tangential component data.

6. The method for determining the splitting of a converted wave shear wave according to claim 1, characterized in that, The step of separately superimposing the radial and tangential component data obtained after synchronization by azimuth angle includes: The radial component gather data obtained after synchronization is divided into several parts according to the azimuth angle, and each part is superimposed according to the bisector angle to obtain the radial component azimuth superimposed gather of the common detector point. The tangential component gather data obtained after synchronization is divided into several parts according to the azimuth angle. Each part is superimposed according to the bisector angle to obtain the tangential component azimuth superimposed gather of the common detector point. The bisected angle is 360 degrees divided by the total number of azimuth angles.

7. A method for determining the splitting of a converted wave shear wave according to any one of claims 1-6, characterized in that, The method of obtaining the shear wave splitting by using the azimuth-based stacked gather of radial and tangential components of the common receiver point, along with the phase reversal method and the phase-constrained least squares method, and then performing the shear wave splitting by vertical layer includes: Estimating the direction of fast wave: On the common detector point superposition gather of the tangential component, the phase reversal method is used to scan along layer H1 to obtain the angle at which the phase reversal occurs, which is the estimated direction of fast wave. Calculate the fast wave direction: Using the estimated fast wave direction as the center, perform least squares fitting on the radial component gather data to calculate the fast wave direction; Calculate the time delay attributes of fast and slow waves: Based on the calculated direction of the fast wave, the cross-correlation method is used to estimate the time delay of the fast and slow waves, and obtain the time delay attributes of layer H1, that is, the time delay attributes of fast and slow waves. One rotation: Based on the calculated fast wave direction, rotate the radial component gather to the fast wave S1 direction and the slow wave S2 direction, and at the same time rotate the tangential component gather to the fast wave S1 direction and the slow wave S2 direction, thus obtaining the fast wave S1 gather and the slow wave S2 gather. Correction: Based on the time delay properties of the fast and slow waves, and in order to align with the direction of the fast wave S1, the slow wave S2 gather is time-corrected to obtain the corrected slow wave S2′ gather. Secondary rotation: Rotate the fast wave S1 gather and the slow wave S2′ gather to their initial radial and tangential directions respectively, thus obtaining the gather data and gather data after shear wave splitting compensation. Repeat: Repeat the steps of "estimate fast wave direction", "calculate fast wave direction", "calculate fast and slow wave delay attributes", "first rotation", "correction" and "second rotation" to perform shear wave splitting calculation and compensation for the next layer after layer H1 until all in-phase axes on the tangential component data have completed shear wave splitting calculation and compensation, thus completing the shear wave splitting acquisition.

8. A system for obtaining the splitting of a converted wave with phase reversal constraint, characterized in that, include: Preprocessing module: used to preprocess, sort, and dynamically correct the radial component data and tangential component data of the converted wave, respectively; Synchronization module: used to synchronize the radial component gather data and tangential component gather data obtained from dynamic correction; Overlay module: used to overlay the radial component data and tangential component data obtained after synchronization by azimuth angle respectively; The acquisition module is used to obtain shear wave splitting based on the azimuth superposition gather of radial and tangential components of the common receiver point obtained by azimuth superposition, as well as the phase reversal method and the phase-constrained least squares method, according to the vertical layer.

9. A phase-reversal constrained converted wave transverse wave splitting determination system according to claim 8, characterized in that, The calculation module includes an estimation submodule, a calculation submodule, a rotation submodule, a correction submodule, and a repetition submodule: Estimation submodule: Used to estimate the fast wave direction: On the common detector point superposition gather of the tangential component, the phase reversal method is used to scan along layer H1 to obtain the angle at which the phase reversal occurs, which is the estimated fast wave direction. The calculation submodule is used to calculate the fast wave direction and the fast and slow wave time delay attributes. Specifically, the fast wave direction is calculated by performing least squares fitting on the radial component gather data with the estimated fast wave direction as the center. The fast and slow wave time delay attributes are calculated by estimating the time delay of the fast and slow waves using the cross-correlation method based on the calculated fast wave direction, thereby obtaining the time delay attributes of layer H1, i.e., the fast and slow wave time delay attributes. The rotation submodule is used to perform a first rotation and a second rotation. The first rotation involves rotating the radial component gather to the fast wave S1 direction and the slow wave S2 direction according to the calculated fast wave direction, and simultaneously rotating the tangential component gather to the fast wave S1 direction and the slow wave S2 direction, thus obtaining the fast wave S1 gather and the slow wave S2 gather. The second rotation involves rotating the fast wave S1 gather and the slow wave S2' gather to the initial radial and tangential directions, respectively, thus obtaining the gather data after shear wave splitting compensation. The correction submodule is used to perform correction: based on the time delay properties of the fast and slow waves, and in order to align with the direction of the fast wave S1, the slow wave S2 gather is time-corrected to obtain the corrected slow wave S2′ gather. The repeat module is used to repeat the steps of "estimating the fast wave direction", "calculating the fast wave direction", "calculating the fast and slow wave delay attributes", "first rotation", "correction" and "second rotation" to perform shear wave splitting calculation and compensation for the next layer after layer H1, until all in-phase axes on the tangential component data have completed shear wave splitting calculation and compensation, that is, the shear wave splitting is obtained.

10. An electronic device, characterized in that, include: Processor, the processor being coupled to memory; The memory is used to store computer programs; The processor is configured to execute the computer program stored in the memory to cause the electronic device to perform the method as described in any one of claims 1 to 7.

11. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a program or instructions that, when executed on a computer, cause the computer to perform the method as described in any one of claims 1 to 7.

Images