Multi-wave data matching method based on P-wave AVO intercept minus gradient attribute

Abstract

This invention belongs to the field of seismic exploration technology for oil and gas, and discloses a multi-wave data matching method based on the P-wave AVO intercept reduction gradient attribute. The method includes: comparative tracking of P-wave marker reflectors; calculation of P-wave AVO intercept reduction gradient data from P-wave CRP gathers; comparative tracking of converted S-wave marker reflectors; low-frequency time difference correction and high-frequency time difference correction of the converted S-wave target layer. This invention is used to obtain converted S-wave data with high-precision matching of P-wave reflection times, effectively reducing multi-wave data matching errors, improving the accuracy of multi-wave data matching within formations, and laying a reliable data foundation for subsequent interpretation and application of multi-wave data.

Description

Technical Field

[0001] This invention belongs to the field of seismic exploration technology for oil and gas, and relates to a multi-wave data matching method, specifically a multi-wave data matching method based on the gradient attribute of the P-wave AVO intercept. Background Technology

[0002] In recent years, multiwave seismic exploration has developed rapidly. Compared with conventional P-wave exploration, multiwave seismic exploration can directly obtain more accurate P-wave and S-wave data of the formation. Through joint interpretation of P-wave and S-wave data, multiple elastic parameters of the formation can be obtained, thereby enabling the prediction of lithology, reservoir, and fluid-bearing properties. Multiwave seismic exploration is also generally referred to as converted S-wave exploration. With the improvement of three-component geophones and supporting acquisition equipment, large-scale digital converted wave exploration has developed rapidly.

[0003] Converted shear wave (SSW) exploration utilizes a "one-shot, two-wave" acquisition method, simultaneously acquiring formation P-wave and converted shear wave reflection data using a three-component geophone. The acquired SSW data accurately reflects the formation's shear wave information. Although years of research into acquisition and processing technologies have largely solved the dual-wave imaging problem, high-precision matching of converted shear wave and P-wave reflections still faces significant technical bottlenecks. Due to the mixed nature of converted shear wave travel time measurement (downward P-wave and upward shear wave), and the differences between the shear wave periodic signal and the P-wave periodic signal, time-domain matching between the two is quite challenging. The commonly used method for matching multi-wave data of converted shear wave and P-wave reflections is as follows: ① Based on the regional geological conditions and the P-wave and S-wave reflection profiles, select the target layer stable reflection marker layer, and use well-seismic calibration to compare and track the reflection marker layer on the P-wave and S-wave profiles; ② Calculate the time difference of the same reflection marker layer for each P-wave and S-wave, smooth the time difference appropriately, use each reflection marker layer of the S-wave as the control layer, and establish a time difference correction body by combining the smoothed time difference of each reflection marker layer; ③ Apply the time difference correction body to correct the converted shear wave data.

[0004] The conventional multi-wave data matching methods described above have the following main problems: ① The number and span of the selected marker layers with the same reflections in both P-waves and S-waves directly affect the accuracy of the target layer's internal matching; the larger the span of the top and bottom marker layers of the target layer, the lower the accuracy of the internal matching. ② Due to the difference in periodic signals between P-waves and S-waves, the difficulty of marker layer correlation varies greatly, easily leading to an increase in the overall matching time difference. ③ For strata with strong spatial heterogeneity in terrestrial clastic rocks, the accuracy requirements for multi-wave data matching within the target layer are even higher. This conventional method actually addresses the low-frequency time difference between multiple waves, but its accuracy for matching high-frequency time differences within the internal layers is insufficient. Therefore, as a crucial link in multi-wave seismic exploration, the accuracy of current techniques cannot meet production needs, and there is an urgent need to develop a new high-precision multi-wave matching technology to further improve the application effectiveness of multi-wave data. Summary of the Invention

[0005] The purpose of this invention is to provide a multi-wave data matching method based on the longitudinal wave AVO intercept subtraction gradient attribute. By using the longitudinal wave AVO intercept subtraction gradient data to perform low-frequency and high-frequency time difference correction on the converted shear wave data, the matching accuracy of the converted shear wave internals is effectively improved, laying a solid foundation for the application of converted shear waves.

[0006] To achieve the above objectives, the technical solution adopted by this invention is as follows:

[0007] A multi-wave data matching method based on the P-wave AVO intercept minus gradient property includes the following steps:

[0008] S1, Comparison and tracking of longitudinal wave marker reflector layer

[0009] S11. Using the P-wave logging curves and density logging curves of wells drilled in the study area, obtain the P-wave reflection composite record and perform well-seismic calibration on the P-wave data;

[0010] S12. Based on the top and bottom positions of the target formation during drilling, determine the top marker reflector (TOP) of the target formation on the P-wave data. P1 and the bottom marker reflector layer of the P-wave target layer (BOT) P1 Location, and mark the top of the P-wave target layer as the TOP reflector layer. P1 The scale corresponds to the top boundary (TOP) of the logging layer, and the bottom boundary (BOT) of the target layer in the P-wave is marked. P1 The scale corresponds to the bottom boundary of the well logging strata (BOT).

[0011] S13. Mark the top of the reflector layer on the P-wave target layer in the P-wave data. P1 and the bottom marker reflector layer of the P-wave target layer (BOT) P1 By conducting comparative tracking of the study area, the corresponding P-wave target layer top time horizon was obtained. P2 P-wave target layer bottom time horizon BOT P2 ;

[0012] S2. Calculate the P-wave AVO intercept subtraction gradient data of the P-wave CRP gather.

[0013] AVO attribute calculations were performed on the P-wave CRP gathers obtained after P-wave pre-stack time migration processing to obtain intercept A and gradient B. The two were then subtracted to obtain P-wave AVO intercept minus gradient data.

[0014] S3, Contrast Tracking of Transverse Wave Marker Reflector Layer

[0015] S31. Using the P-wave logging curves, S-wave logging curves, and density logging curves of wells in the study area, obtain the converted S-wave composite record and perform well-seismic calibration on the converted S-wave data.

[0016] S32. Based on the top boundary (TOP) and bottom boundary (BOT) of the logging layer in step S12, determine the top marker reflector (TOP) of the target layer in the converted shear wave data. S1 and conversion transverse wave target layer bottom marker reflector layer BOT S1 Location;

[0017] S33. The top time layer of the target layer in step S13 is the P-wave target layer TOP. P2 P-wave target layer bottom time horizon BOT P2 Projected onto the longitudinal wave AVO intercept subtraction gradient data from step S2;

[0018] S34. Compare the waveform characteristics of the projection location on the AVO intercept subtraction gradient data of the reference longitudinal wave with the top marker reflector layer (TOP) of the target layer in the converted transverse wave. S1 and conversion transverse wave target layer bottom marker reflector layer BOT S1 By comparing and tracing the waveform characteristics of the location within the study area, the corresponding converted shear wave top time layer TOP was obtained. S2 , conversion of shear wave target layer bottom time layer BOT S2 ;

[0019] S4, Low-frequency time difference correction for converted shear wave target layer

[0020] S41. Calculate the time difference ΔT at the top time level and the time difference ΔB at the bottom time level between the longitudinal wave and the converted shear wave, and smooth them respectively to obtain the top low-frequency time difference data ΔTsmooth and the bottom low-frequency time difference data ΔBsmooth.

[0021] S42. Based on the planar coordinates, analyze the top low-frequency time difference data ΔTsmooth and the converted shear wave target layer top time level (TOP). S2 As well as the bottom low-frequency time difference data △Bsmooth and the converted shear wave target layer bottom time layer BOT S2 Perform matching separately to establish a low-frequency time difference correction data volume △Tlf;

[0022] S43. Apply the low-frequency time difference correction data volume △Tlf to the converted shear wave data to perform time difference correction, and obtain the low-frequency time difference corrected converted shear wave data PS_lf.

[0023] S5, High-frequency time difference correction for converted shear wave target layer

[0024] S51. Define the correction sliding window Tcor, the sliding overlap time T, and the maximum time difference correction amount;

[0025] S52, using the top time horizon of the target layer in the P-wave phase wave P2Starting from each seismic trace, compress or stretch the TOP. P2 Within the +Tcor range, calculate the cross-correlation between the low-frequency time difference corrected converted shear wave data PS_lf and the simultaneous window P-wave AVO intercept subtraction gradient data after each compression or stretching. Record the corresponding compression or stretching amount ΔThf when the cross-correlation value is maximum. Slide the time window downwards to overlap the top of the previous time window with the bottom of the previous time window for time T, and repeat the above steps until the time window crosses the bottom time layer (BOT) of the P-wave target layer. P2 Record the compression or stretching amount △Thf for each time window;

[0026] The above calculation is iterated for each seismic trace until the calculation of the entire data volume is completed, and the high-frequency time difference correction data volume △Thf1 is obtained.

[0027] S53. Apply the high-frequency time difference correction data volume △Thf1 to the low-frequency time difference correction converted shear wave data PS_lf to perform time difference correction, and obtain the low-frequency and high-frequency time difference correction converted shear wave data PS_lf_hf.

[0028] As a limitation, the P-wave data is P-wave pre-stack time offset data, and the converted S-wave data is converted S-wave pre-stack time offset data.

[0029] As a second limitation, in step S12, the top marker reflector layer (TOP) of the target layer for the P-wave is determined on the P-wave data. P1 and the bottom marker reflector layer of the P-wave target layer (BOT) P1 When determining the location, select a marker reflection layer with clear and continuous waveform characteristics at the top and bottom of the target layer.

[0030] As a third limitation, in step S41, the time difference ΔT of the top time level is TOP. P2 -TOP S2 The time difference ΔB at the bottom time level = BOT P2 -BOT S2 .

[0031] As a fourth limitation, in step S52, the correction sliding window is limited to 60ms, the sliding overlap time is limited to 20ms, and the maximum time difference correction amount is limited to -10ms to 10ms.

[0032] The present invention, by adopting the above-described technical solution, achieves the following technical advancements compared to existing technologies:

[0033] (1) This invention uses the gradient reduction data of the longitudinal wave AVO intercept to perform low-frequency time difference correction and high-frequency time difference correction on the converted shear wave data, which effectively improves the matching accuracy of the converted shear wave internal structure, lays a solid foundation for the application of converted shear waves, and improves the application accuracy of converted shear wave exploration in the comprehensive evaluation and fine description of oil and gas reservoirs.

[0034] (2) This invention utilizes P-wave data for well seismic calibration, preferentially selects the top marker reflection layer and the bottom marker reflection layer of the P-wave target layer, conducts comparative tracking, and then projects it onto the P-wave AVO intercept sub-gradient data. Since both the P-wave AVO intercept sub-gradient data and the converted shear wave data characterize the formation shear wave information and have a high degree of similarity in periodic signals, they can effectively guide the comparison of the same marker reflection layer in the converted shear wave data, thereby improving the reliability and accuracy of the marker reflection layer comparison.

[0035] (3) Since the periodic signals of the longitudinal wave AVO intercept subtraction gradient data and the converted shear wave data are highly similar, the present invention utilizes the longitudinal wave AVO intercept subtraction gradient data to perform high-frequency time difference correction on the converted shear wave data after low-frequency time difference correction, which can effectively improve the matching accuracy of multi-wave data in the formation.

[0036] (4) This invention can effectively reduce the matching error of multi-wave data, improve the matching accuracy of multi-wave data in the formation, and lay a reliable data foundation for the subsequent interpretation and application of multi-wave data.

[0037] This invention is used for multi-wave data matching to obtain high-precision converted shear wave data with matching P-wave reflection time, which can effectively reduce multi-wave data matching error. Attached Figure Description

[0038] Figure 1 The diagram shown is a flowchart of a method according to an embodiment of the present invention. Detailed Implementation

[0039] To better explain and facilitate understanding of the present invention, the present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.

[0040] An Example: A Multi-Wave Data Matching Method Based on P-Wave AVO Intercept Subtraction Gradient Attribute

[0041] In this embodiment, the multi-wave seismic data is based on P-waves and converted S-waves that have undergone high-fidelity pre-stack time migration, belonging to the data matching and processing stage after the pre-stack time migration of P-waves and converted S-waves. Therefore, the P-wave data is P-wave pre-stack time migration data, and the converted S-wave data is converted S-wave pre-stack time migration data.

[0042] like Figure 1 As shown, this embodiment includes the following steps:

[0043] S1, Comparison and tracking of longitudinal wave marker reflector layer

[0044] S11. Using the P-wave logging curves and density logging curves of wells drilled in the study area, obtain the P-wave reflection composite record and perform well-seismic calibration on the P-wave data;

[0045] S12. Based on the top and bottom positions of the target formation during drilling, determine the top marker reflector (TOP) of the target formation on the P-wave data. P1 and the bottom marker reflector layer of the P-wave target layer (BOT) P1 Location, and mark the top of the P-wave target layer as the TOP reflector layer. P1 The scale corresponds to the top boundary (TOP) of the logging layer, and the bottom boundary (BOT) of the target layer in the P-wave is marked. P1 The scale corresponds to the bottom boundary of the well logging strata (BOT).

[0046] Among them, the top marker reflector layer (TOP) of the target layer in the P-wave data is determined. P1 and the bottom marker reflector layer of the P-wave target layer (BOT) P1 When locating, select a marker reflection layer with clear and continuous waveform characteristics at the top and bottom of the target layer;

[0047] S13. Mark the top of the reflector layer on the P-wave target layer in the P-wave data. P1 and the bottom marker reflector layer of the P-wave target layer (BOT) P1 By conducting comparative tracking of the study area, the corresponding P-wave target layer top time horizon was obtained. P2 P-wave target layer bottom time horizon BOT P2 ;

[0048] S2. Calculate the P-wave AVO intercept subtraction gradient data of the P-wave CRP gather.

[0049] AVO (Amplitude Variation with Offset) attribute calculations were performed on the P-wave CRP gathers obtained after P-wave pre-stack time migration processing to obtain the intercept A and gradient B. The two were then subtracted, i.e., AB, to obtain the P-wave AVO intercept minus gradient data.

[0050] In this step, after P-wave pre-stack time migration processing, a P-wave CRP gather with high AVO fidelity can be obtained. The P-wave CRP gather refers to the P-wave common reflection point reflection angle gather. The P-wave AVO intercept subtraction gradient data and the P-wave data have the same time measurement, which is the P-wave reflection time.

[0051] S3, Contrast Tracking of Transverse Wave Marker Reflector Layer

[0052] S31. Using the P-wave logging curves, S-wave logging curves, and density logging curves of wells in the study area, obtain the converted S-wave composite record and perform well-seismic calibration on the converted S-wave data.

[0053] S32. Based on the top boundary (TOP) and bottom boundary (BOT) of the logging layer in step S12, determine the top marker reflector (TOP) of the target layer in the converted shear wave data. S1 and conversion transverse wave target layer bottom marker reflector layer BOT S1 Location;

[0054] S33. The top time layer of the target layer in step S13 is the P-wave target layer TOP. P2 P-wave target layer bottom time horizon BOT P2 Projected onto the longitudinal wave AVO intercept subtraction gradient data from step S2;

[0055] S34. Compare the waveform characteristics of the projection location on the AVO intercept subtraction gradient data of the reference longitudinal wave with the top marker reflector layer (TOP) of the target layer in the converted transverse wave. S1 and conversion transverse wave target layer bottom marker reflector layer BOT S1 By comparing and tracing the waveform characteristics of the location within the study area, the corresponding converted shear wave top time layer TOP was obtained. S2 , conversion of shear wave target layer bottom time layer BOT S2 ;

[0056] In this step, based on the Aki_Richards equation:

[0057] R P (θ)=A+B sin 2 θ+C tan 2 θsin 2 θ,where;

[0058]

[0059] Assuming Vs / Vp = 1 / 2, then

[0060]

[0061] where: and

[0062] Among them, R P (θ) represents the longitudinal wave reflection coefficient at an incident angle of θ, A represents the intercept, B represents the gradient, C represents the curvature, and V P V represents the longitudinal wave velocity. S ΔV represents the transverse wave velocity. P The longitudinal wave velocity difference, ΔVS R represents the transverse wave velocity difference, ρ represents density, Δρ represents density difference, and R po R represents the longitudinal wave reflection coefficient at zero incident angle. so This represents the reflection coefficient of a transverse wave at zero angle of incidence.

[0063] The derivation yields:

[0064] R so =0.5(AB);

[0065] That is, intercept A minus gradient B corresponds proportionally to the zero-incident-angle shear wave reflection coefficient. Therefore, the intercept A minus gradient B data obtained by calculating the AVO attribute based on the above Aki_Richards equation in the seismic band is equivalent to shear wave reflection data.

[0066] Based on the above analysis, it can be concluded that both P-wave AVO intercept subtraction gradient data and converted S-wave data can characterize the S-wave reflection information of the formation, and the periodic signal waveforms of the two are highly similar.

[0067] S4, Low-frequency time difference correction for converted shear wave target layer

[0068] S41. Calculate the time difference ΔT at the top time level and the time difference ΔB at the bottom time level between the longitudinal wave and the converted shear wave, and smooth them respectively to obtain the top low-frequency time difference data ΔTsmooth and the bottom low-frequency time difference data ΔBsmooth.

[0069] Among them, the time difference ΔT of the top time level is TOP P2 -TOP S2 The time difference ΔB at the bottom time level = BOT P2 -BOT S2 ;

[0070] S42. Based on the planar coordinates, analyze the top low-frequency time difference data ΔTsmooth and the converted shear wave target layer top time level (TOP). S2 As well as the bottom low-frequency time difference data △Bsmooth and the converted shear wave target layer bottom time layer BOT S2 Perform matching separately to establish a low-frequency time difference correction data volume △Tlf;

[0071] S43. Apply the low-frequency time difference correction data volume △Tlf to the converted shear wave data to perform time difference correction, and obtain the low-frequency time difference corrected converted shear wave data PS_lf.

[0072] S5, High-frequency time difference correction for converted shear wave target layer

[0073] S51. Define the correction sliding window Tcor, the sliding overlap time T, and the maximum time difference correction amount;

[0074] S52, using the top time horizon of the target layer in the P-wave phase wave P2 Starting from each seismic trace, compress or stretch the TOP. P2 Within the +Tcor range, calculate the cross-correlation between the low-frequency time difference corrected converted shear wave data PS_lf and the simultaneous window P-wave AVO intercept subtraction gradient data after each compression or stretching. Record the corresponding compression or stretching amount ΔThf when the cross-correlation value is maximum. Slide the time window downwards to overlap the top of the previous time window with the bottom of the previous time window for time T, and repeat the above steps until the time window crosses the bottom time layer (BOT) of the P-wave target layer. P2 Record the compression or stretching amount △Thf for each time window; iterate the above calculation for each seismic trace until the calculation of the entire data volume is completed, and obtain the high-frequency time difference correction data volume △Thf1;

[0075] The correction sliding window is limited to 60ms, the sliding overlap time is limited to 20ms, and the compression or stretching range, i.e. the maximum time difference correction amount, is limited to -10ms to 10ms.

[0076] S53. Apply the high-frequency time difference correction data volume △Thf1 to the low-frequency time difference correction converted shear wave data PS_lf to perform time difference correction, and obtain the low-frequency and high-frequency time difference correction converted shear wave data PS_lf_hf.

[0077] The obtained converted shear wave data PS_lf_hf after low-frequency and high-frequency time difference correction is the high-precision converted shear wave data with P-wave reflection time after low-frequency and high-frequency time difference correction. This data has higher matching accuracy with the P-wave data.

[0078] It should be noted that the above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art can still modify the technical solutions described in the above embodiments or make equivalent substitutions for some of the technical features. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A multi-wave data matching method based on P-wave AVO intercept minus gradient attribute, characterized in that, The method comprises the following steps: S1, contrast tracing of the P-wave marker reflection layer S11, obtaining P-wave reflection synthetic records by using P-wave logging curves and density logging curves of the wells in the study area, and calibrating the P-wave data by well-seismic calibration; S12, according to the top and bottom positions of the drilling target layer, determining the P-wave target layer top marker reflection layer TOP on the P-wave data P1 and the P-wave target layer bottom marker reflection layer BOT P1 position, and scaling the P-wave target layer top marker reflection layer TOP P1 to the logging layer top boundary TOP, and scaling the P-wave target layer bottom marker reflection layer BOT P1 to the logging layer bottom boundary BOT; S13, on the longitudinal wave data, the top mark reflection layer TOP of the longitudinal wave target layer P1 and the bottom mark reflection layer BOT of the longitudinal wave target layer P1 Research area contrast tracking is carried out, and the corresponding top time horizon TOP of the longitudinal wave target layer is obtained P2 , the bottom time horizon BOT of the longitudinal wave target layer P2 ; S2, calculating P-wave AVO intercept minus gradient data of P-wave CRP gathers Carrying out AVO attribute calculation on the P-wave CRP gathers obtained by P-wave pre-stack time migration processing, obtaining intercept A and gradient B, and subtracting the two to obtain P-wave AVO intercept minus gradient data; S3, contrast tracing of the S-wave marker reflection layer S31, obtaining converted S-wave synthetic records by using P-wave logging curves, S-wave logging curves and density logging curves of the wells in the study area, and calibrating the converted S-wave data by well-seismic calibration; S32, determining the converted shear wave target layer top marker reflection layer TOP and the converted shear wave target layer bottom marker reflection layer BOT on the converted shear wave data according to the well logging layer top boundary TOP and the well logging layer bottom boundary BOT in step S12 S1 and the converted shear wave target layer bottom marker reflection layer BOT S1 position; S33, project the P-wave AVO intercept minus gradient data in step S2 onto the P-wave target layer top time horizon TOP P2 S34, project the P-wave AVO intercept minus gradient data in step S2 onto the P-wave target layer bottom time horizon BOT P2 S33, project the P-wave AVO intercept minus gradient data in step S2 onto the P-wave target layer top time horizon TOP S34, the waveform feature of the position where the projection on the P-wave AVO intercept minus gradient data is located, compared with the converted S-wave target layer top marker reflection layer TOP S1 and the converted S-wave target layer bottom marker reflection layer BOT S1 , the waveform feature of the position where the projection is located, the study area contrast tracking is carried out, and the corresponding converted S-wave target layer top time horizon TOP S2 , the converted S-wave target layer bottom time horizon BOT S2 ; S4, low-frequency moveout correction of the target layer of converted S-wave S41, calculating the time difference △T of the top time horizon and the time difference △B of the bottom time horizon between the P-wave and the converted S-wave, and respectively performing smoothing to obtain top low-frequency moveout data △Tsmooth and bottom low-frequency moveout data △Bsmooth; S42. Based on the planar coordinates, analyze the top low-frequency time difference data ΔTsmooth and the converted shear wave target layer top time level (TOP). S2 As well as the bottom low-frequency time difference data △Bsmooth and the converted shear wave target layer bottom time layer BOT S2 Perform matching separately to establish a low-frequency time difference correction data volume △Tlf; S43, applying the low-frequency moveout correction data body △Tlf to the converted S-wave data to perform moveout correction, and obtaining the converted S-wave data PS_lf after low-frequency moveout correction; S5, high-frequency moveout correction of the target layer of converted S-wave S51, defining a correction sliding time window Tcor, a sliding overlap time T, and a maximum moveout correction amount; S52, the top time horizon of the target layer of the longitudinal wave is compressed or stretched to obtain a time window P2 Take each seismic trace as a unit, compress or stretch the TOP P2 The converted shear wave data PS_lf after low-frequency moveout correction in the range of +Tcor, calculate the cross-correlation of each compressed or stretched converted shear wave data PS_lf after low-frequency moveout correction and the simultaneous window AVO intercept minus gradient data, and record the corresponding compression or stretching amount△Thf when the cross-correlation value is maximum; slide the time window downward with the top of the time window overlapping the bottom of the previous time window by a time T, repeat the above steps until the time window crosses the bottom time horizon BOT of the target layer of the longitudinal wave P2 , record the compression or stretching amount△Thf of each time window; iterating the above operation on each seismic trace until the operation of the entire data body is completed, and obtaining a high-frequency moveout correction data body △Thf1; S53, applying the high-frequency moveout correction data body △Thf1 to the converted S-wave data PS_lf after low-frequency moveout correction to perform moveout correction, and obtaining converted S-wave data PS_lf_hf after low-frequency and high-frequency moveout correction.

2. The multi-wave data matching method based on P-wave AVO intercept minus gradient attribute of claim 1, wherein, The P-wave data is P-wave pre-stack time migration data, and the converted S-wave data is converted S-wave pre-stack time migration data.

3. The multi-wave data matching method based on P-wave AVO intercept minus gradient attribute of claim 1, wherein, In step S12, the top marker reflection layer TOP of the P-wave target layer is determined on the P-wave data P1 and the bottom marker reflection layer BOT of the P-wave target layer P1 When the positions of the top and bottom marker reflection layers of the target layer are determined, the marker reflection layers with clear waveform features and good continuity are selected.

4. The multi-wave data matching method based on P-wave AVO intercept minus gradient attribute of claim 1, wherein, In step S41, the time difference of the top time horizon is ΔT = TOP P2 -TOP S2 ; and the time difference of the bottom time horizon is ΔB = BOT P2 -BOT S2 .

5. The method of claim 1, wherein, In step S52, the correction sliding time window is limited to 60 ms, the sliding overlap time is limited to 20 ms, and the maximum moveout correction amount is limited to -10 ms to 10 ms.

Images