Method and apparatus for matching longitudinal and transverse waves, computing device, and storage medium
By performing P-wave and S-wave matching processing and matched filtering techniques in the depth domain, the P-wave layer is adjusted to match the converted wave layer, solving the problem of imaging gas cloud areas in P-wave exploration, achieving efficient and accurate P-wave and S-wave matching, and improving the imaging effect.
Patent Information
- Application Number
- CN202210885887.1
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2022-07-26
- Publication Date
- 2026-02-24
- Estimated Expiration
- 2042-07-26
AI Technical Summary
In existing technologies, P-wave exploration is significantly affected when imaging gas cloud areas, while time-domain calibration is difficult in converted wave processing, especially the unstable adjustment of the P-wave/S-wave velocity ratio, which affects the converted wave processing effect.
By performing P-wave and S-wave matching processing in the depth domain, depth domain result data of P-wave and converted wave are obtained. Stratification and matching filtering are then performed, adjusting the P-wave horizon to match the converted wave horizon. The S-wave velocity adjustment factor is used to adjust the P-wave horizon interpretation results, and finally, the results are converted and processed in the time domain.
It improves the efficiency and accuracy of P-wave and S-wave matching processing, solves the problem of time-domain P-wave and S-wave matching, and ensures the accuracy and consistency of imaging results.
Smart Images
Figure CN115222976B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of marine petroleum geophysical exploration technology, specifically to a method, apparatus, computing device, and storage medium for P-wave and S-wave matching processing. Background Technology
[0002] With the exploration and development of oil fields, the requirements for the imaging quality of seismic data are becoming increasingly stringent. The presence of gas cloud regions poses a significant risk to exploration and development. P-wave exploration is greatly affected by gas cloud regions, while converted waves are less affected by fluids. Using converted waves to image gas cloud regions is an excellent technique. A major challenge in converted wave processing is the time-domain calibration problem. The calibration process requires adjusting the P-wave / S-wave velocity ratio (i.e., the gamma field) and the converted wave velocity field. This method involves two variables, making calibration extremely difficult, and its instability directly affects the converted wave processing results. Summary of the Invention
[0003] In view of the above problems, the present invention is proposed to provide a method, apparatus, computing device and storage medium for longitudinal and transverse wave matching processing that overcomes or at least partially solves the above problems.
[0004] According to one aspect of the present invention, a P-wave and S-wave matching processing method is provided, the method comprising:
[0005] Acquire P-wave depth domain results data and converted wave depth domain results data;
[0006] Layer interpretation was performed on the P-wave depth domain data to obtain the P-wave layer interpretation results;
[0007] Stratigraphic interpretation was performed on the converted wave depth domain data to obtain the converted wave stratigraphic interpretation results;
[0008] The P-wave and converted wave interpretation results are superimposed according to the layers to obtain the P-wave and S-wave matching results in the depth domain.
[0009] The P-wave and S-wave matching results in the depth domain are transformed to obtain the P-wave and S-wave matching results in the time domain.
[0010] Optionally, the method further includes:
[0011] Based on the attribute information of the converted wave depth domain result data, the P-wave depth domain result data is subjected to matched filtering to obtain the filtered P-wave depth domain result data.
[0012] Further stratigraphic interpretation of the P-wave depth domain data includes:
[0013] Stratified interpretation is performed on the filtered P-wave depth domain data.
[0014] Optionally, the P-wave layer interpretation results and the converted wave layer interpretation results are superimposed according to the layers to obtain the P-wave and S-wave matching results in the depth domain, which further include:
[0015] Determine the first layer depth information of the target P-wave layer in the P-wave layer interpretation results, and the second layer depth information of the corresponding converted wave layer in the converted wave layer interpretation results;
[0016] The depth change information is calculated based on the first and second layer depth information;
[0017] The shear wave velocity adjustment factor is calculated based on the depth change information;
[0018] The P-wave layer interpretation results are adjusted based on the shear wave velocity adjustment factor.
[0019] The adjusted P-wave layer interpretation results and the converted wave layer interpretation results are superimposed according to the layer to obtain the P-wave and S-wave matching results in the depth domain.
[0020] Optionally, the P-wave depth domain data is interpreted to obtain P-wave layer interpretation results, which further include:
[0021] The main marker layers in the P-wave horizon are identified, and the data corresponding to the main marker layers in the P-wave depth domain results are interpreted to obtain the P-wave horizon interpretation results.
[0022] The converted wave depth domain data is then interpreted to obtain converted wave stratigraphic interpretation results, which further include:
[0023] The target converted wave layer data corresponding to the main marker layer in the converted wave depth domain results are determined, and the data of the target converted wave layer are interpreted to obtain the converted wave layer interpretation results.
[0024] Optionally, performing matched filtering on the P-wave depth domain result data based on the attribute information of the converted wave depth domain result data further includes:
[0025] Based on the attribute information of the converted wave depth domain results data, the matched filter operator is obtained;
[0026] A filter algorithm is constructed based on the matched filter operator, and the matched filter algorithm is used to perform matched filtering on the P-wave depth domain results data.
[0027] Optionally, the matched filtering operator includes: amplitude matched filtering operator, phase matched filtering operator and / or frequency matched filtering operator.
[0028] According to another aspect of the present invention, a P-wave and S-wave matching processing apparatus is provided, the apparatus comprising:
[0029] The data acquisition module is suitable for acquiring P-wave depth domain result data and converted wave depth domain result data;
[0030] The stratigraphic interpretation module is suitable for performing stratigraphic interpretation on P-wave depth domain data to obtain P-wave stratigraphic interpretation results; and for performing stratigraphic interpretation on converted wave depth domain data to obtain converted wave stratigraphic interpretation results.
[0031] The layer matching module is suitable for superimposing the P-wave layer interpretation results and the converted wave layer interpretation results according to the layer to obtain the P-wave and S-wave matching results in the depth domain.
[0032] The conversion module is suitable for converting the P-wave and S-wave matching results in the depth domain to obtain the P-wave and S-wave matching results in the time domain.
[0033] Optionally, the device further includes:
[0034] The matching filtering module is suitable for performing matching filtering on the P-wave depth domain result data based on the attribute information of the converted wave depth domain result data, so as to obtain the filtered P-wave depth domain result data.
[0035] The stratigraphic interpretation module is further adapted to perform stratigraphic interpretation on the filtered P-wave depth domain results data to obtain P-wave stratigraphic interpretation results.
[0036] Optionally, the layer matching module is further adapted to: determine the first layer depth information of the target P-wave layer in the P-wave layer interpretation result, and the second layer depth information of the corresponding converted wave layer in the converted wave layer interpretation result; calculate the depth variation information based on the first and second layer depth information; calculate the shear wave velocity adjustment factor based on the depth variation information; adjust the P-wave layer interpretation result based on the shear wave velocity adjustment factor; and superimpose the adjusted P-wave layer interpretation result with the converted wave layer interpretation result according to the layer to obtain the P-wave and S-wave matching result in the depth domain.
[0037] Optionally, the stratigraphic interpretation module is further adapted to: determine the main marker layer in the P-wave stratigraphic sequence, perform stratigraphic interpretation on the data corresponding to the main marker layer in the P-wave depth domain result data, and obtain the P-wave stratigraphic interpretation result; and determine the data of the target converted wave stratigraphic sequence corresponding to the main marker layer in the converted wave depth domain result data, perform stratigraphic interpretation on the data of the target converted wave stratigraphic sequence, and obtain the converted wave stratigraphic interpretation result.
[0038] Optionally, the matching filtering module is further adapted to: obtain a matching filter operator based on the attribute information of the converted wave depth domain result data; construct a filter algorithm based on the matching filter operator; and use the filter algorithm to perform matched filtering processing on the P-wave depth domain result data.
[0039] Optionally, the matched filtering operator includes: amplitude matched filtering operator, phase matched filtering operator and / or frequency matched filtering operator.
[0040] According to another aspect of the present invention, a computing device is provided, comprising: a processor, a memory, a communication interface, and a communication bus, wherein the processor, the memory, and the communication interface communicate with each other via the communication bus;
[0041] The memory is used to store at least one executable instruction, which causes the processor to perform the operation corresponding to the above-described transverse wave matching processing method.
[0042] According to another aspect of the present invention, a computer storage medium is provided, the storage medium storing at least one executable instruction that causes a processor to perform an operation corresponding to the above-described transverse wave matching processing method.
[0043] According to the present invention, the method, apparatus, computing device, and storage medium for P-wave and S-wave matching processing include: acquiring P-wave depth domain result data and converted wave depth domain result data; performing stratigraphic interpretation on the P-wave depth domain result data to obtain P-wave stratigraphic interpretation results; performing stratigraphic interpretation on the converted wave depth domain result data to obtain converted wave stratigraphic interpretation results; superimposing the P-wave stratigraphic interpretation results and the converted wave stratigraphic interpretation results according to stratigraphic layers to obtain P-wave and S-wave matching results in the depth domain; and performing conversion processing on the P-wave and S-wave matching results in the depth domain to obtain P-wave and S-wave matching results in the time domain. Through this method, P-wave and S-wave matching processing is performed directly in the depth domain, and after completing the depth domain P-wave and S-wave matching processing, a time-depth conversion process is performed to convert the depth domain P-wave and S-wave matching results into time domain P-wave and S-wave matching results, thereby solving the problem of time domain P-wave and S-wave matching. This method is simple and easy to implement, and can improve the efficiency and accuracy of P-wave and S-wave matching processing.
[0044] The above description is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it in accordance with the contents of the specification, and to make the above and other objects, features and advantages of the present invention more apparent and understandable, specific embodiments of the present invention are described below. Attached Figure Description
[0045] Various other advantages and benefits will become apparent to those skilled in the art upon reading the following detailed description of preferred embodiments. The accompanying drawings are for illustrative purposes only and are not intended to limit the invention. Furthermore, the same reference numerals denote the same parts throughout the drawings. In the drawings:
[0046] Figure 1A flowchart of the P-wave and S-wave matching processing method provided in an embodiment of the present invention is shown;
[0047] Figure 2a A schematic diagram of the time-domain imaging of P-wave results is shown in an example;
[0048] Figure 2b A schematic diagram of the converted wave results in the time domain is shown in an example;
[0049] Figure 2c A schematic diagram of the imaging of P-wave results in the depth domain is shown in an example;
[0050] Figure 2d A schematic diagram of the converted wave results in the depth domain imaging example is shown;
[0051] Figure 3 A flowchart of a P-wave and S-wave matching processing method provided in another embodiment of the present invention is shown;
[0052] Figure 4 This diagram illustrates a comparison of P-wave and converted wave data before and after matched filtering in an example.
[0053] Figure 5a A schematic diagram of the pre-matched filter P-wave imaging result is shown in one example;
[0054] Figure 5b A schematic diagram of the P-wave imaging results after matched filtering is shown in one example;
[0055] Figure 6a A schematic diagram of the interpretation results of P-wave results in the depth domain is shown in an example;
[0056] Figure 6b A schematic diagram of the joint interpretation results of P-wave and S-wave depth domain is shown in an example;
[0057] Figure 7 A schematic diagram of the longitudinal and transverse wave matching processing device provided in an embodiment of the present invention is shown.
[0058] Figure 8 A schematic diagram of the structure of a computing device provided in an embodiment of the present invention is shown. Detailed Implementation
[0059] Exemplary embodiments of the invention will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the invention are shown in the drawings, it should be understood that the invention may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this invention will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
[0060] Figure 1 A flowchart of the P-wave and S-wave matching processing method provided in an embodiment of the present invention is shown, as follows: Figure 1 As shown, the longitudinal and transverse wave matching processing method in this embodiment specifically includes the following steps:
[0061] Step S110: Obtain the P-wave depth domain result data and the converted wave depth domain result data;
[0062] Step S120: Perform stratigraphic interpretation on the P-wave depth domain data to obtain the P-wave stratigraphic interpretation results;
[0063] Step S130: Perform stratigraphic interpretation on the converted wave depth domain data to obtain the converted wave stratigraphic interpretation results;
[0064] Step S140: The P-wave layer interpretation results and the converted wave layer interpretation results are superimposed according to the layers to obtain the P-wave and S-wave matching results in the depth domain.
[0065] Step S150: The P-wave and S-wave matching results in the depth domain are transformed to obtain the P-wave and S-wave matching results in the time domain.
[0066] Because the velocities of P-waves and converted waves are inconsistent—P-waves are emitted from a P-wave source and received, both having the P-wave velocity—and the velocities of converted waves are the sum of the P-wave velocity and the S-wave velocity, while the S-wave propagation velocity is lower than the P-wave velocity, the time scales of converted wave and P-wave seismic reflections are not consistent in the time domain. This difference in propagation velocities between P and S-waves leads to inconsistencies in the imaging locations of the same strata in the time domain, making direct stratigraphic matching interpretation impossible.
[0067] Figure 2a A schematic diagram of the time-domain imaging of P-wave results is shown in an example. Figure 2b The diagram illustrates the time-domain imaging of converted wave data in an example. Comparison reveals inconsistencies in the imaging locations of the same stratum in the P-wave and S-wave seismic data within the time domain. Therefore, precise matching of P-wave and converted wave seismic data is essential before subsequent interpretation and inversion processes can be performed.
[0068] Although the propagation speeds of longitudinal waves and converted waves are different, they have the same imaging depth at the same location in the depth domain. Figure 2c A schematic diagram of the imaging of P-wave results in the depth domain is shown as an example. Figure 2d The diagram shows an example of converted wave imaging in the depth domain. The comparison shows that the imaging positions of the P-wave and the converted wave in the depth domain are at the same depth, that is, the depth scale of the converted wave seismic reflection and the P-wave seismic reflection in the depth domain is the same, and the converted wave and the P-wave can be directly matched in the depth domain.
[0069] Based on this, the method in this embodiment performs P-wave and S-wave matching processing directly in the depth domain. After completing the P-wave and S-wave matching processing in the depth domain, it performs time-depth conversion processing to convert the P-wave and S-wave matching results in the depth domain into P-wave and S-wave matching results in the time domain, thereby solving the problem of P-wave and S-wave matching in the time domain. This method is simple and easy to implement and can improve the efficiency and accuracy of P-wave and S-wave matching processing.
[0070] Figure 3 A flowchart of a P-wave and S-wave matching processing method provided in another embodiment of the present invention is shown, as follows: Figure 3 As shown, the method includes the following steps:
[0071] Step S310: Obtain the P-wave depth domain result data and the converted wave depth domain result data;
[0072] Step S320: Based on the attribute information of the converted wave depth domain result data, perform matched filtering on the P-wave depth domain result data to obtain the filtered P-wave depth domain result data.
[0073] There are many differences between the P-wave depth domain results data and the converted wave depth domain results data. These differences can affect the accuracy of the joint P-wave and S-wave layer matching interpretation. Based on this, the method in this embodiment adopts matched filtering technology to perform consistency processing on the P-wave depth domain results data based on the attributes of the converted wave depth domain results data.
[0074] In one alternative approach, step S320 specifically includes the following sub-steps: obtaining a matched filter operator based on the attribute information of the converted wave depth domain result data; constructing a filter algorithm based on the matched filter operator; and using the filter algorithm to perform matched filtering processing on the P-wave depth domain result data.
[0075] Matched filtering, in particular, uses one seismic trace as the desired output and calculates a matching filter operator to apply to another seismic trace, thereby eliminating amplitude, frequency, and / or phase differences between seismic traces. Matched filtering methods are available in single-trace and multi-trace versions, as well as pre-stack and post-stack versions. This approach can ensure good consistency between data from two seismic traces.
[0076] The attribute information of the converted wave depth domain result data includes amplitude information, phase information, and / or frequency information. The matched filtering operators include amplitude matched filtering operators, phase matched filtering operators, and / or frequency matched filtering operators. Using the matched filtering operators, a matched filtering algorithm is constructed to perform amplitude matching correction processing, phase matching correction processing, and / or frequency matching correction processing on the P-wave depth domain result data to eliminate the differences between the P-wave depth domain result data and the converted wave depth domain result data in amplitude, phase, and / or frequency.
[0077] Amplitude matching refers to spectral shaping techniques used to correct the amplitude of seismic records to a specified shape. This involves performing a Fourier transform on the seismic record, converting it to the frequency domain, calculating its amplitude spectrum, and then weighting the amplitude spectrum according to the amplitude frequency while keeping the phase spectrum unchanged. Finally, an inverse Fourier transform is performed to obtain the spectrally shaped seismic record. In this embodiment, based on the amplitude information of the converted wave depth domain data, an amplitude matching filter operator is obtained to eliminate the amplitude differences between the P-wave depth domain data and the converted wave depth domain data.
[0078] Regarding phase matching, due to significant differences in field excitation factors and inconsistent construction factors, such as varying air intensity depths and cable depths, phase differences exist in the in-phase axes of data from the same location in different work areas. Therefore, phase matching correction is necessary. In the actual acquisition process, a work area with better data quality is first selected as the operator matching direction for other work areas, i.e., the splicing direction. Before phase correction, amplitude correction processing is required. Then, the in-phase axes of the same strata can be well distinguished in the overlapping areas of the spliced work areas. Therefore, data from this overlapping area is selected, and the phase matching operator can be obtained on a single offset plane to achieve the purpose of phase correction. In this embodiment, the seismic traces corresponding to the converted wave depth domain results data are used as the operator matching direction to obtain the phase matching filter operator, which is used to eliminate the phase difference between the P-wave depth domain results data and the converted wave depth domain results data.
[0079] Regarding frequency matching, the frequency difference between P-waves and S-waves is the primary issue affecting the joint matching interpretation of P-wave and S-wave layers. P-waves decay energy slowly and acquire more signals, while S-waves decay energy quickly and acquire fewer signals. Existing techniques fail to meet the requirements for P-wave and S-wave matching interpretation, especially since the converted wave's wave group characteristics are poor and significantly different from the P-wave's, leading to inconsistencies in wave group characteristics. Therefore, matching processing is needed for both P-wave depth domain data and converted wave depth domain data to reduce their differences. In this embodiment, a frequency matching filter operator is derived based on the frequency information of the converted wave depth domain data to eliminate the frequency differences between the P-wave and converted wave depth domain data.
[0080] Figure 4 The figure shows a schematic diagram comparing the P-wave and converted wave data before and after matched filtering in an example. Curve 41 is the frequency-amplitude curve of the P-wave before matched filtering, curve 42 is the frequency-amplitude curve of the P-wave after matched filtering, and curve 43 is the frequency-amplitude curve of the converted wave. As shown in the figure, the difference between the frequencies of the P-wave and the converted wave is reduced after matched filtering. Figure 5a A schematic diagram of the P-wave imaging results before matched filtering is shown in one example. Figure 5b A schematic diagram of the P-wave imaging results after matched filtering is shown in an example. Comparison reveals that the resolution of the P-wave imaging is reduced after matched filtering. It is evident that in the method of this embodiment, the P-wave depth domain results data are consistent based on the attributes of the converted wave depth domain results data, resulting in good consistency between the two. Maintaining consistency between the two is beneficial for subsequent stratigraphic interpretation and helps improve the accuracy of P-wave and S-wave matching.
[0081] In actual production, matched filtering should be performed before pre-stack deconvolution. The reason is that after matched filtering before convolution, the wavelets and spectra of the seismic data become largely consistent, but some residual differences still exist. Deconvolution, especially multi-channel deconvolution, can further eliminate these differences. Further residual static correction can eliminate minor time differences caused by different excitation sources, ensuring consistent reflection phase axes. Finally, hybrid stacking can minimize these differences, resulting in more natural stitching points.
[0082] Step S330: Perform stratigraphic interpretation on the filtered P-wave depth domain data to obtain the P-wave stratigraphic interpretation results.
[0083] Step S340: Perform stratigraphic interpretation on the converted wave depth domain data to obtain the converted wave stratigraphic interpretation results.
[0084] In one optional approach, step S330 specifically includes the following sub-steps: determining the main marker layer in the P-wave stratigraphy, and interpreting the data corresponding to the main marker layer in the filtered P-wave depth domain data to obtain the P-wave stratigraphic interpretation result; wherein, the main marker layer refers to a layer or group of rock strata with obvious characteristics that can serve as stratigraphic correlation markers. The main marker layer should have the characteristics of obvious fossil and lithological features, stable stratigraphy, wide distribution, and easy identification. Correspondingly, step S340 specifically includes the following sub-steps: determining the data of the target converted wave stratigraphy corresponding to the main marker layer in the converted wave depth domain data, and interpreting the data of the target converted wave stratigraphy to obtain the converted wave stratigraphic interpretation result. In this approach, the main marker layer is first found in the P-wave stratigraphy and the interpretation work is completed. Then, the interpretation result is loaded onto the converted wave migration profile according to the stratigraphy to interpret the corresponding stratigraphy and obtain the converted wave stratigraphic interpretation result. This method can reduce the amount of data processed, improve the processing effect, and improve the accuracy of the matching results.
[0085] Step S350: Determine the first layer depth information of the target P-wave layer in the P-wave layer interpretation results, and the second layer depth information of the corresponding converted wave layer in the converted wave layer interpretation results; calculate the depth change information based on the first layer depth information and the second layer depth information.
[0086] When loading the P-wave layer interpretation results into the converted wave layer interpretation results, in practical applications, there may be inconsistencies in the depth position of the same layer in the two interpretation results. There are also certain differences between P-wave imaging and converted wave imaging. Therefore, adjustments are still needed to ensure that the depth position of the same layer is consistent.
[0087] The target P-wave layer can be one or a group of layers in the P-wave layer, or it can be the main marker layer mentioned above. This is equivalent to calculating the depth change based on the first and second depths of the same layer in the P-wave layer interpretation results and the converted wave layer interpretation results.
[0088] Step S360: Calculate the shear wave velocity adjustment factor based on the depth change information; adjust the P-wave layer interpretation results based on the shear wave velocity adjustment factor.
[0089] After calculating the depth change, the shear wave velocity change information, i.e. the shear wave velocity adjustment factor, is calculated according to the formula relating depth and velocity. Then, the P-wave layer interpretation result is multiplied by the shear wave velocity change information to obtain the adjusted P-wave layer interpretation result.
[0090] The formula relating depth and velocity is as follows:
[0091]
[0092] Where D represents depth, T PS Indicating the travel time of the converted wave, V P V represents the longitudinal wave velocity. S This represents the shear wave velocity. The P-wave depth D is obtained after completing the P-wave layer interpretation during the layer interpretation process. PP After being loaded with the converted wave, the local stratigraphic depth is D. PS According to D PP With D PS The difference between the two values gives the depth change Dx, which, when substituted into the above formula, yields the corresponding change in shear wave velocity.
[0093] Step S370: The adjusted P-wave layer interpretation results and the converted wave layer interpretation results are superimposed according to the layers to obtain the P-wave and S-wave matching results in the depth domain.
[0094] Finally, the adjusted P-wave layer interpretation results and the converted wave layer interpretation results are superimposed according to the layers to obtain the P-wave and S-wave matching results in the depth domain.
[0095] Figure 6a A schematic diagram of the interpretation results of P-wave results in the depth domain is shown in an example. Figure 6bThis diagram illustrates a joint depth-domain interpretation of P-waves and S-waves in an example. The P-wave depth-domain interpretation is loaded into the converted wave depth-domain interpretation, resulting in... Figure 6b The results of the joint interpretation of P-wave and S-wave depth domains are shown.
[0096] Step S380: The P-wave and S-wave matching results in the depth domain are transformed to obtain the P-wave and S-wave matching results in the time domain.
[0097] The P-wave and S-wave matching results in the depth domain are subjected to time-depth conversion processing to obtain the P-wave and S-wave matching results in the time domain. Specifically, based on the correspondence between time and depth, the P-wave and S-wave matching results in the depth domain are divided by a specified velocity to obtain the P-wave and S-wave matching results in the time domain. That is, the geological strata in the depth domain are converted to the geological strata in the time domain.
[0098] The correspondence between P-wave time or converted wave time and depth is as follows:
[0099] T = D / 2V
[0100] Where T represents time, D represents depth, and V represents velocity.
[0101] According to the P-wave and S-wave matching processing method provided in this embodiment, on the one hand, a consistent spectrum matching technique is provided. This technique processes the P-wave depth domain results data using matched filtering to eliminate differences in amplitude, frequency, and / or phase between the P-wave and converted wave results, ensuring good consistency between them. This effectively solves the problem of wave group characteristic consistency during P-wave and S-wave matching, facilitating subsequent stratigraphic interpretation. On the other hand, a depth domain P-wave and S-wave joint stratigraphic interpretation technique based on velocity updates is provided. During the joint interpretation process, S-wave variations are determined based on depth differences, and the P-wave is then adjusted according to these variations. The wave interpretation results are adjusted to achieve consistent depth of P-wave and S-wave layers, effectively carrying out joint layer interpretation of P-wave and S-wave in the depth domain, which is beneficial for layer interpretation of imaging results and S-wave velocity inversion. On the other hand, based on the consistent depth of imaging points of P-wave and S-wave in the depth domain, P-wave and S-wave matching is first performed in the depth domain, and then based on the relationship between depth, velocity and time, the P-wave and S-wave matching results in the depth domain are converted into P-wave and S-wave matching results in the time domain. This can reduce the difficulty of P-wave and S-wave matching in the time domain, improve the matching accuracy of P-wave results and converted wave results, and solve the problem of high difficulty and poor effect of P-wave and S-wave matching in the existing technology.
[0102] Figure 7 A schematic diagram of the structure of the longitudinal and transverse wave matching processing device provided in an embodiment of the present invention is shown, as follows: Figure 7 As shown, the device includes:
[0103] Data acquisition module 71 is suitable for acquiring P-wave depth domain result data and converted wave depth domain result data;
[0104] The stratigraphic interpretation module 72 is suitable for performing stratigraphic interpretation on P-wave depth domain result data to obtain P-wave stratigraphic interpretation results; and for performing stratigraphic interpretation on converted wave depth domain result data to obtain converted wave stratigraphic interpretation results.
[0105] The layer matching module 73 is suitable for superimposing the P-wave layer interpretation results and the converted wave layer interpretation results according to the layer to obtain the P-wave and S-wave matching results in the depth domain.
[0106] The conversion module 74 is suitable for converting the P-wave and S-wave matching results in the depth domain to obtain the P-wave and S-wave matching results in the time domain.
[0107] In one alternative embodiment, the apparatus further includes: a matching filtering module, adapted to perform matching filtering on the P-wave depth domain result data based on the attribute information of the converted wave depth domain result data, to obtain filtered P-wave depth domain result data;
[0108] The stratigraphic interpretation module 72 is further adapted to: perform stratigraphic interpretation on the filtered P-wave depth domain result data to obtain the P-wave stratigraphic interpretation result.
[0109] In an optional manner, the layer matching module 73 is further adapted to: determine the first layer depth information of the target P-wave layer in the P-wave layer interpretation result, and the second layer depth information of the corresponding converted wave layer in the converted wave layer interpretation result; calculate the depth variation information based on the first layer depth information and the second layer depth information; calculate the shear wave velocity adjustment factor based on the depth variation information; adjust the P-wave layer interpretation result based on the shear wave velocity adjustment factor; and superimpose the adjusted P-wave layer interpretation result and the converted wave layer interpretation result according to the layer to obtain the P-wave and S-wave matching result in the depth domain.
[0110] Optionally, the stratigraphic interpretation module 72 is further adapted to: determine the main marker layer in the P-wave stratigraphic sequence, perform stratigraphic interpretation on the data corresponding to the main marker layer in the P-wave depth domain result data, and obtain the P-wave stratigraphic interpretation result; and determine the data of the target converted wave stratigraphic sequence corresponding to the main marker layer in the converted wave depth domain result data, perform stratigraphic interpretation on the data of the target converted wave stratigraphic sequence, and obtain the converted wave stratigraphic interpretation result.
[0111] Optionally, the matching filtering module is further adapted to: obtain a matching filter operator based on the attribute information of the converted wave depth domain result data; construct a filter algorithm based on the matching filter operator; and use the filter algorithm to perform matched filtering processing on the P-wave depth domain result data.
[0112] Optionally, the matched filtering operator includes: amplitude matched filtering operator, phase matched filtering operator and / or frequency matched filtering operator.
[0113] By using the above method, P-wave and S-wave matching processing is performed directly in the depth domain. After completing the P-wave and S-wave matching processing in the depth domain, time-depth conversion processing is performed to convert the P-wave and S-wave matching results in the depth domain into P-wave and S-wave matching results in the time domain, thereby solving the problem of P-wave and S-wave matching in the time domain. This method is simple and easy to implement, and can improve the efficiency and accuracy of P-wave and S-wave matching processing.
[0114] This invention provides a non-volatile computer storage medium storing at least one executable instruction that can execute the longitudinal and transverse wave matching processing method in any of the above method embodiments.
[0115] Figure 8 The diagram shows a structural schematic of an embodiment of the computing device of the present invention. The specific embodiments of the present invention do not limit the specific implementation of the computing device.
[0116] like Figure 8 As shown, the computing device may include: a processor 802, a communications interface 804, a memory 806, and a communications bus 808.
[0117] The processor 802, communication interface 804, and memory 806 communicate with each other via communication bus 808. Communication interface 804 is used to communicate with other network elements such as clients or other servers. Processor 802 executes program 810, specifically performing the relevant steps in the above-described embodiment of the longitudinal and transverse wave matching processing method for computing devices.
[0118] Specifically, program 810 may include program code that includes computer operation instructions.
[0119] The processor 802 may be a central processing unit (CPU), an application-specific integrated circuit (ASIC), or one or more integrated circuits configured to implement embodiments of the present invention. The computing device includes one or more processors, which may be processors of the same type, such as one or more CPUs; or processors of different types, such as one or more CPUs and one or more ASICs.
[0120] Memory 806 is used to store program 810. Memory 806 may include high-speed RAM memory, and may also include non-volatile memory, such as at least one disk storage device.
[0121] The algorithms or displays provided herein are not inherently related to any particular computer, virtual system, or other device. Various general-purpose systems can also be used in conjunction with the teachings herein. The required structure for constructing such systems is apparent from the above description. Furthermore, the embodiments of the present invention are not directed to any particular programming language. It should be understood that the content of the invention described herein can be implemented using various programming languages, and the above description of specific languages is for the purpose of disclosing the best mode of implementation of the invention.
[0122] Numerous specific details are set forth in the specification provided herein. However, it will be understood that embodiments of the invention may be practiced without these specific details. In some instances, well-known methods, structures, and techniques have not been shown in detail so as not to obscure the understanding of this specification.
[0123] Similarly, it should be understood that, in order to streamline the invention and aid in understanding one or more of the various inventive aspects, features of the embodiments of the invention are sometimes grouped together in a single embodiment, figure, or description thereof in the above description of exemplary embodiments of the invention. However, this disclosure should not be construed as reflecting an intention that the claimed invention requires more features than expressly recited in each claim. Rather, as reflected in the following claims, inventive aspects lie in fewer than all features of a single foregoing disclosed embodiment. Therefore, the claims following the detailed description are hereby expressly incorporated into that detailed description, wherein each claim itself is a separate embodiment of the invention.
[0124] Those skilled in the art will understand that modules in the device of the embodiments can be adaptively changed and placed in one or more devices different from that embodiment. Modules, units, or components in the embodiments can be combined into a single module, unit, or component, and further, they can be divided into multiple sub-modules, sub-units, or sub-components. Except where at least some of such features and / or processes or units are mutually exclusive, any combination can be employed to combine all features disclosed in this specification (including the accompanying claims, abstract, and drawings) and all processes or units of any method or device so disclosed. Unless expressly stated otherwise, each feature disclosed in this specification (including the accompanying claims, abstract, and drawings) may be replaced by an alternative feature that serves the same, equivalent, or similar purpose.
[0125] Furthermore, those skilled in the art will understand that although some embodiments herein include certain features included in other embodiments but not others, combinations of features from different embodiments are intended to be within the scope of the invention and form different embodiments. For example, in the following claims, any of the claimed embodiments can be used in any combination.
[0126] The various component embodiments of the present invention can be implemented in hardware, or as software modules running on one or more processors, or a combination thereof. Those skilled in the art will understand that microprocessors or digital signal processors (DSPs) can be used in practice to implement some or all of the functions of some or all of the components according to the embodiments of the present invention. The present invention can also be implemented as a device or apparatus program (e.g., a computer program and computer program product) for performing part or all of the methods described herein. Such programs implementing the present invention can be stored on a computer-readable medium, or can be in the form of one or more signals. Such signals can be downloaded from an Internet website, provided on a carrier signal, or provided in any other form.
[0127] It should be noted that the above embodiments are illustrative of the invention and not restrictive, and that those skilled in the art can devise alternative embodiments without departing from the scope of the appended claims. In the claims, any reference signs placed between parentheses should not be construed as limiting the claims. The word "comprising" does not exclude the presence of elements or steps not listed in the claims. The word "a" or "an" preceding an element does not exclude the presence of multiple such elements. The invention can be implemented by means of hardware comprising several different elements and by means of a suitably programmed computer. In the unit claims enumerating several means, several of these means may be embodied by the same item of hardware. The use of the words first, second, and third, etc., does not indicate any order. These words can be interpreted as names. The steps in the above embodiments, unless otherwise specified, should not be construed as limiting the order of execution.
Claims
1. A method for processing P-wave and S-wave matching, characterized in that, The method includes: Acquire P-wave depth domain results data and converted wave depth domain results data; The P-wave depth domain data is interpreted to obtain the P-wave layer interpretation results; Stratigraphic interpretation is performed on the converted wave depth domain data to obtain converted wave stratigraphic interpretation results; The P-wave layer interpretation results and the converted wave layer interpretation results are superimposed according to the layer to obtain the P-wave and S-wave matching results in the depth domain. The P-wave and S-wave matching results in the depth domain are transformed to obtain the P-wave and S-wave matching results in the time domain. Specifically, based on the correspondence between time and depth, the P-wave and S-wave matching results in the depth domain are divided by a specified velocity to obtain the P-wave and S-wave matching results in the time domain.
2. The method according to claim 1, characterized in that, The method further includes: Based on the attribute information of the converted wave depth domain result data, the P-wave depth domain result data is subjected to matched filtering to obtain filtered P-wave depth domain result data. The hierarchical interpretation of the P-wave depth domain data further includes: The filtered P-wave depth domain data is then interpreted at different levels.
3. The method according to claim 1, characterized in that, The step of superimposing the P-wave layer interpretation results and the converted wave layer interpretation results according to layers to obtain the P-wave and S-wave matching results in the depth domain further includes: Determine the first layer depth information of the target P-wave layer in the P-wave layer interpretation result, and the second layer depth information of the converted wave layer corresponding to the target P-wave layer in the converted wave layer interpretation result; Depth change information is calculated based on the first layer depth information and the second layer depth information; Based on the depth change information, the shear wave velocity adjustment factor is calculated; The P-wave layer interpretation results are adjusted based on the shear wave velocity adjustment factor. The adjusted P-wave layer interpretation results are superimposed with the converted wave layer interpretation results according to the layer to obtain the P-wave and S-wave matching results in the depth domain.
4. The method according to claim 1, characterized in that, The step of interpreting the P-wave depth domain data to obtain the P-wave layer interpretation results further includes: The main marker layers in the P-wave layers are identified, and the data corresponding to the main marker layers in the P-wave depth domain results are interpreted to obtain the P-wave layer interpretation results. The step of performing stratigraphic interpretation on the converted wave depth domain data to obtain the converted wave stratigraphic interpretation results further includes: The target converted wave layer data corresponding to the main marker layer in the converted wave depth domain results are determined, and the data of the target converted wave layer are interpreted to obtain the converted wave layer interpretation results.
5. The method according to claim 2, characterized in that, The step of performing matched filtering on the P-wave depth domain result data based on the attribute information of the converted wave depth domain result data further includes: Based on the attribute information of the converted wave depth domain result data, the matched filter operator is obtained; A filter algorithm is constructed based on the matched filter operator, and the matched filter algorithm is used to perform matched filtering processing on the P-wave depth domain result data.
6. The method according to claim 5, characterized in that, The matched filtering operators include: amplitude matched filtering operator, phase matched filtering operator and / or frequency matched filtering operator.
7. A P-wave and S-wave matching processing device, characterized in that, The device includes: The data acquisition module is suitable for acquiring P-wave depth domain result data and converted wave depth domain result data; The stratigraphic interpretation module is adapted to perform stratigraphic interpretation on the P-wave depth domain result data to obtain P-wave stratigraphic interpretation results; and to perform stratigraphic interpretation on the converted wave depth domain result data to obtain converted wave stratigraphic interpretation results. The layer matching module is adapted to superimpose the P-wave layer interpretation results and the converted wave layer interpretation results according to the layers to obtain the P-wave and S-wave matching results in the depth domain. The conversion module is adapted to convert the P-wave and S-wave matching results in the depth domain to obtain the P-wave and S-wave matching results in the time domain. Specifically, based on the correspondence between time and depth, the P-wave and S-wave matching results in the depth domain are divided by a specified velocity to obtain the P-wave and S-wave matching results in the time domain.
8. The apparatus according to claim 7, characterized in that, The device further includes: The matching filtering module is adapted to perform matching filtering on the P-wave depth domain result data based on the attribute information of the converted wave depth domain result data, so as to obtain the filtered P-wave depth domain result data. The stratigraphic interpretation module is further adapted to: perform stratigraphic interpretation on the filtered P-wave depth domain result data to obtain P-wave stratigraphic interpretation results.
9. A computing device, comprising: The processor, memory, communication interface, and communication bus are provided, wherein the processor, memory, and communication interface communicate with each other via the communication bus. The memory is used to store at least one executable instruction, which causes the processor to perform the operation corresponding to the longitudinal and transverse wave matching processing method as described in any one of claims 1-6.
10. A computer storage medium storing at least one executable instruction that causes a processor to perform an operation corresponding to the P-wave and S-wave matching processing method as described in any one of claims 1-6.