Pre-stack gather energy equalization processing method, device and equipment and storage medium

Abstract

The application provides a pre-stack gather energy balancing processing method, device and equipment and a storage medium, and comprises the following steps: collecting post-stack seismic data, pre-stack CRP gathers and a marked geological horizon in a known complex structure area; calculating a pre-stack CRP gather structure energy correction factor by using the post-stack seismic data, the pre-stack CRP gathers and the marked geological horizon; and obtaining corrected pre-stack CRP gathers by using the pre-stack CRP gather structure energy correction factor and the pre-stack CRP gathers. The application can correct the energy of the pre-stack CRP gathers under the geological conditions of complex structures and high-steep structures, and can improve the accuracy of pre-stack seismic inversion and further reflect the real elasticity, physical properties and fluid of the reservoir more accurately without reprocessing.

Description

Technical Field

[0001] This application relates to the field of geophysical exploration, and in particular to a method, apparatus, equipment and storage medium for pre-stack gather energy equalization processing. Background Technology

[0002] The quality of pre-stack CRP gathers directly affects the signal-to-noise ratio and accuracy of pre-stack elastic inversion. If the pre-stack CRP gathers do not meet the requirements of pre-stack elastic inversion, even advanced inversion algorithms will affect the final inversion results. Therefore, before conducting pre-stack elastic inversion, it is essential to perform data quality evaluation and amplitude-preserving denoising processing on the pre-stack CRP gathers to provide a high-quality seismic data foundation for subsequent pre-stack elastic inversion.

[0003] Typical pre-stack gather optimization steps include the following: ① Denoising, primarily to improve the signal-to-noise ratio (SNR) of the pre-stack gather. ② Angle-shifting gathers, mainly converting the offset-domain pre-stack CRP gathers to the angle domain based on the velocity field for pre-stack inversion. ③ Residual time difference correction, mainly to flatten the gather's phase axis and eliminate spurious AVO phenomena during the inversion process. ④ Super-gather processing, with a similar purpose to denoising, improving the SNR through pixel stacking. ⑤ Angular partial stacking, stacking approximately 40 channels from a single CDP angle gather into 3-5 angle partial stacked data as input for pre-stack inversion; this stacking improves the stability of the pre-stack inversion. Sometimes, specific compensation methods are employed to address irregularities in the acquisition and observation system.

[0004] However, few people notice the systematic differences in energy among pre-stack gathers in steep or other complex structural regions. These energy differences typically arise during seismic data processing. Seismic data processing for steep structures usually involves ray compensation after data stacking to achieve energy equalization of post-stack seismic data. However, energy differences due to structural morphology still exist in pre-stack CRP gathers. If these structurally-induced energy differences are not equalized and are carried over into pre-stack inversion, they will ultimately lead to biases in our understanding of the reservoir. Summary of the Invention

[0005] To address the aforementioned issues, this application provides a method, apparatus, device, and storage medium for pre-stack gather energy equalization processing.

[0006] This application provides a pre-stack gather energy equalization processing method, including:

[0007] S1: Data preparation: Collect post-stack seismic data, pre-stack CRP gathers, and key geological horizons for known complex tectonic zones;

[0008] S2: Calculation of pre-stack CRP gather tectonic energy correction factor: The pre-stack CRP gather tectonic energy correction factor is calculated using the post-stack seismic data, pre-stack CRP gathers and key geological horizons.

[0009] S3: Calculation of the corrected pre-stack CRP gather: The energy correction factor and the pre-stack CRP gather are constructed by applying the pre-stack CRP gather to obtain the corrected pre-stack CRP gather.

[0010] In some embodiments, the specific method for calculating the tectonic energy correction factor of the pre-stack CRP gather using the post-stack seismic data, pre-stack CRP gathers, and key geological horizons includes:

[0011] (1) Extract the root mean square amplitude attribute of post-stack seismic data with a large time window;

[0012] (2) Extract the large time window root mean square amplitude attribute of the pre-stack CRP gather;

[0013] (3) The large time window root mean square amplitude attribute of the post-stack seismic data and the large time window root mean square amplitude attribute of the pre-stack CRP gather are used to obtain the structural energy correction factor of the pre-stack CRP gather.

[0014] In some embodiments, the specific method for extracting the large time window root mean square amplitude attribute of post-stack seismic data includes:

[0015] Based on key geological strata, large time window root mean square amplitude attributes of post-stack seismic data are extracted.

[0016] In some embodiments, the specific method for extracting the large time window root mean square amplitude attribute of the pre-stack CRP gather includes:

[0017] (1) Perform full stacking of pre-stack CRP gathers to obtain post-stack data based on direct stacking of pre-stack CRP gathers;

[0018] (2) Based on the landmark geological strata, the large time window root mean square amplitude attribute is extracted from the post-stack data directly superimposed from the pre-stack CRP gathers to obtain the large time window root mean square amplitude attribute of the pre-stack CRP gathers.

[0019] In some embodiments, the selection principle for the large time window is that the time window is not less than 200ms, and all data within the time window is valid throughout the entire work area.

[0020] In some embodiments, the specific method for obtaining the pre-stack CRP gather tectonic energy correction factor by applying the large-window root-mean-square amplitude attribute of the post-stack seismic data and the large-window root-mean-square amplitude attribute of the pre-stack CRP gather includes:

[0021] The structural energy correction factor of the pre-stack CRP gather = the root mean square amplitude attribute of the large time window of the post-stack seismic data / the root mean square amplitude attribute of the large time window of the pre-stack CRP gather.

[0022] In some embodiments, the specific method for constructing an energy correction factor and a pre-stack CRP gather using the pre-stack CRP gather to obtain the corrected pre-stack CRP gather includes:

[0023] Corrected pre-stack CRP gather = pre-stack CRP gather construction energy correction factor × pre-stack CRP gather.

[0024] This application provides a pre-stack gather energy equalization processing device, including:

[0025] Data preparation module, pre-stack CRP gather construction energy correction factor calculation module, and corrected pre-stack CRP gather calculation module;

[0026] The data preparation module collects post-stack seismic data, pre-stack CRP gathers, and key geological strata in known complex tectonic zones.

[0027] The pre-stack CRP gather tectonic energy correction factor calculation module: calculates the pre-stack CRP gather tectonic energy correction factor using the post-stack seismic data, pre-stack CRP gathers and key geological strata;

[0028] The corrected pre-stack CRP gather calculation module: applies the pre-stack CRP gather to construct the energy correction factor and the pre-stack CRP gather, and obtains the corrected pre-stack CRP gather.

[0029] This application provides a pre-stack gather energy equalization processing device, including a memory and a processor. The memory stores a computer program, which, when executed by the processor, performs any of the pre-stack gather energy equalization processing methods described above.

[0030] This application provides a storage medium storing a computer program that can be executed by one or more processors and can be used to implement the pre-stack gather energy equalization processing method described in any of the above claims.

[0031] This application provides a method, apparatus, equipment, and storage medium for pre-stack gather energy equalization processing.

[0032] To address the systematic energy differences in pre-stack CRP gathers under complex and steep geological conditions, a method for energy correction of pre-stack gathers is proposed without requiring reprocessing. This improves the accuracy of pre-stack seismic inversion and more accurately reflects the true elastic, physical, and fluid properties of the reservoir. Attached Figure Description

[0033] The present application will be described in more detail below based on embodiments and with reference to the accompanying drawings.

[0034] Figure 1 A schematic diagram illustrating the implementation process of a pre-stack gather energy equalization processing method provided in an embodiment of this application;

[0035] Figure 2 The marker layer Wufeng-Longmaxi Formation layer diagram provided for embodiments of this application;

[0036] Figure 3 A cross-sectional view of the short axis direction based on direct stacking of pre-stack CRP gathers, provided for an embodiment of this application;

[0037] Figure 4 The pre-stack CRP gather stacking data provided in this application are based on a large time window root mean square amplitude attribute map of a landmark geological stratum;

[0038] Figure 5 The post-stack result data provided for the embodiments of this application are based on the large time window root mean square amplitude attribute map of the landmark geological strata;

[0039] Figure 6 This is a planar property map of the pre-stack gather energy correction factor provided in an embodiment of this application;

[0040] Figure 7 The energy-corrected pre-stack gather stacking data provided in the embodiments of this application are based on a large time window root mean square amplitude attribute map of a landmark geological stratum;

[0041] In the accompanying drawings, the same parts are referred to by the same reference numerals, and the drawings are not drawn to scale. Detailed Implementation

[0042] To make the objectives, technical solutions, and advantages of this application clearer, the application will be further described in detail below with reference to the accompanying drawings. The described embodiments should not be regarded as limitations on this application. All other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0043] In the following description, references are made to “some embodiments,” which describe a subset of all possible embodiments. However, it is understood that “some embodiments” may be the same subset or different subsets of all possible embodiments and may be combined with each other without conflict.

[0044] If the application documents contain similar descriptions such as "first, second, third", the following explanation shall be added: In the following description, the terms "first, second, third" are used only to distinguish similar objects and do not represent a specific order of objects. It is understood that "first, second, third" may be interchanged in a specific order or sequence where permitted, so that the embodiments of this application described herein can be implemented in an order other than that illustrated or described herein.

[0045] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used herein is for the purpose of describing embodiments of this application only and is not intended to limit this application.

[0046] Before introducing the pre-stack gather energy equalization processing method provided in the embodiments of this application, a brief introduction is given to the problems existing in the related technologies:

[0047] To address the problems existing in related technologies, this application provides a pre-stack gather energy equalization processing method. This method is applied to a pre-stack gather energy equalization processing device, which can be an electronic device, such as a computer or mobile terminal. The functions implemented by the pre-stack gather energy equalization processing method provided in this application can be achieved by the processor of the electronic device calling program code, wherein the program code can be stored in a computer storage medium.

[0048] Example 1

[0049] This application provides a pre-stack gather energy equalization processing method. Figure 1 This is a schematic diagram illustrating the implementation process of a pre-stack gather energy equalization processing method provided in an embodiment of this application, as shown below. Figure 1 As shown, it includes:

[0050] S1: Data Preparation: Collecting post-stack seismic data from known complex tectonic zones Pre-stack CRP collection and iconic geological strata ;

[0051] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0052] S2: Prestack CRP gather structure energy correction factor Calculation: Applying the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0053] S3: Corrected pre-stack CRP gather Calculation: Energy correction factor constructed using pre-stack CRP gathers and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. .

[0054] This application provides a pre-stack gather energy equalization processing method. Under complex and steep geological conditions, pre-stack CRP gathers exhibit systematic energy differences related to the structure. Without requiring reprocessing, this method proposes a way to correct the energy of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and more accurately reflecting the true elastic, physical, and fluid properties of the reservoir.

[0055] Example 2

[0056] Based on the foregoing embodiments, this application further provides a pre-stack gather energy equalization processing method, including:

[0057] S1: Data Preparation: Collecting post-stack seismic data from known complex tectonic zones Pre-stack CRP collection and iconic geological strata ;

[0058] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0059] S2: Prestack CRP gather structure energy correction factor Calculation: Applying the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0060] In some embodiments, the application of the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor Specific methods include:

[0061] (1) Extracting post-stack seismic data Large time window root mean square amplitude property ;

[0062] (2) Extracting pre-stack CRP gathers Large time window root mean square amplitude property ;

[0063] (3) Apply the large time window root mean square amplitude attribute of the post-stack seismic data Large time window root mean square amplitude property of pre-stack CRP gather The pre-stack CRP gather construction energy correction factor was obtained. ;

[0064] S3: Corrected pre-stack CRP gather Calculation: Energy correction factor constructed using pre-stack CRP gathers and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. .

[0065] This application provides a pre-stack gather energy equalization processing method. Under complex and steep geological conditions, pre-stack CRP gathers exhibit systematic energy differences related to the structure. Without requiring reprocessing, this method proposes a way to correct the energy of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and more accurately reflecting the true elasticity, physical properties, and fluidity of the reservoir.

[0066] Example 3

[0067] Based on the foregoing embodiments, this application further provides a pre-stack gather energy equalization processing method, including:

[0068] S1: Data Preparation: Collecting post-stack seismic data from known complex tectonic zones Pre-stack CRP collection and iconic geological strata ;

[0069] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0070] S2: Prestack CRP gather structure energy correction factor Calculation: Applying the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0071] In some embodiments, the application of the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor Specific methods include:

[0072] (1) Extracting post-stack seismic data Large time window root mean square amplitude property ;

[0073] (2) Extracting pre-stack CRP gathers Large time window root mean square amplitude property ;

[0074] (3) Apply the large time window root mean square amplitude attribute of the post-stack seismic data Large time window root mean square amplitude property of pre-stack CRP gather The pre-stack CRP gather construction energy correction factor was obtained. ;

[0075] In some embodiments, the specific method for extracting the large time window root mean square amplitude attribute of post-stack seismic data includes:

[0076] Based on key geological strata Extracting the large time window root mean square amplitude attribute of post-stack seismic data ;

[0077] S3: Corrected pre-stack CRP gather Calculation: Energy correction factor constructed using pre-stack CRP gathers and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. .

[0078] This application provides a pre-stack gather energy equalization processing method. Under complex and steep geological conditions, pre-stack CRP gathers exhibit systematic energy differences related to the structure. Without requiring reprocessing, this method proposes a way to correct the energy of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and more accurately reflecting the true elasticity, physical properties, and fluidity of the reservoir.

[0079] Example 4

[0080] Based on the foregoing embodiments, this application further provides a pre-stack gather energy equalization processing method, including:

[0081] S1: Data Preparation: Collecting post-stack seismic data from known complex tectonic zones Pre-stack CRP collection and iconic geological strata ;

[0082] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0083] S2: Prestack CRP gather structure energy correction factor Calculation: Applying the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0084] In some embodiments, the application of the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor Specific methods include:

[0085] (1) Extracting post-stack seismic data Large time window root mean square amplitude property ;

[0086] (2) Extracting pre-stack CRP gathers Large time window root mean square amplitude property ;

[0087] (3) Apply the large time window root mean square amplitude attribute of the post-stack seismic data Large time window root mean square amplitude property of pre-stack CRP gather The pre-stack CRP gather construction energy correction factor was obtained. ;

[0088] In some embodiments, the specific method for extracting the large time window root mean square amplitude attribute of post-stack seismic data includes:

[0089] Based on key geological strata Extracting the large time window root mean square amplitude attribute of post-stack seismic data ;

[0090] In some embodiments, the extraction of the large time window root mean square amplitude attribute of the pre-stack CRP gather Specific methods include:

[0091] (1) Pre-stack CRP gather Perform full stacking to obtain post-stack data based on direct stacking of pre-stack CRP gathers. ;

[0092] (2) Based on landmark geological strata For post-stack data based on direct stacking of pre-stack CRP gathers Large time window root mean square amplitude attributes are extracted to obtain the large time window root mean square amplitude attributes of the pre-stack CRP gather. ;

[0093] S3: Corrected pre-stack CRP gather Calculation: Energy correction factor constructed using pre-stack CRP gathers and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. .

[0094] This application provides a pre-stack gather energy equalization processing method. Under complex and steep geological conditions, pre-stack CRP gathers exhibit systematic energy differences related to the structure. Without requiring reprocessing, this method proposes a way to correct the energy of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and more accurately reflecting the true elasticity, physical properties, and fluidity of the reservoir.

[0095] Example 5

[0096] Based on the foregoing embodiments, this application further provides a pre-stack gather energy equalization processing method, including:

[0097] S1: Data Preparation: Collecting post-stack seismic data from known complex tectonic zones Pre-stack CRP collection and iconic geological strata ;

[0098] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0099] S2: Prestack CRP gather structure energy correction factor Calculation: Applying the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0100] In some embodiments, the application of the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor Specific methods include:

[0101] (1) Extracting post-stack seismic data Large time window root mean square amplitude property ;

[0102] (2) Extracting pre-stack CRP gathers Large time window root mean square amplitude property ;

[0103] (3) Apply the large time window root mean square amplitude attribute of the post-stack seismic data Large time window root mean square amplitude property of pre-stack CRP gather The pre-stack CRP gather construction energy correction factor was obtained. ;

[0104] In some embodiments, the specific method for extracting the large time window root mean square amplitude attribute of post-stack seismic data includes:

[0105] Based on key geological strata Extracting the large time window root mean square amplitude attribute of post-stack seismic data ;

[0106] In some embodiments, the extraction of the large time window root mean square amplitude attribute of the pre-stack CRP gather Specific methods include:

[0107] (1) Pre-stack CRP gather Perform full stacking to obtain post-stack data based on direct stacking of pre-stack CRP gathers. ;

[0108] (2) Based on landmark geological strata For post-stack data based on direct stacking of pre-stack CRP gathers Large time window root mean square amplitude attributes are extracted to obtain the large time window root mean square amplitude attributes of the pre-stack CRP gather. ;

[0109] In some embodiments, the selection principle for the large time window is that the time window is not less than 200ms, and all data within the time window is valid across the entire work area;

[0110] S3: Corrected pre-stack CRP gather Calculation: Energy correction factor constructed using pre-stack CRP gathers and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. .

[0111] This application provides a pre-stack gather energy equalization processing method. Under complex and steep geological conditions, pre-stack CRP gathers exhibit systematic energy differences related to the structure. Without requiring reprocessing, this method proposes a way to correct the energy of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and more accurately reflecting the true elasticity, physical properties, and fluidity of the reservoir.

[0112] Example 6

[0113] Based on the foregoing embodiments, this application further provides a pre-stack gather energy equalization processing method, including:

[0114] S1: Data Preparation: Collecting post-stack seismic data from known complex tectonic zones Pre-stack CRP collection and iconic geological strata ;

[0115] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0116] S2: Prestack CRP gather structure energy correction factor Calculation: Applying the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0117] In some embodiments, the application of the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor Specific methods include:

[0118] (1) Extracting post-stack seismic data Large time window root mean square amplitude property ;

[0119] (2) Extracting pre-stack CRP gathers Large time window root mean square amplitude property ;

[0120] (3) Apply the large time window root mean square amplitude attribute of the post-stack seismic data Large time window root mean square amplitude property of pre-stack CRP gather The pre-stack CRP gather construction energy correction factor was obtained. ;

[0121] In some embodiments, the specific method for extracting the large time window root mean square amplitude attribute of post-stack seismic data includes:

[0122] Based on key geological strata Extracting the large time window root mean square amplitude attribute of post-stack seismic data ;

[0123] In some embodiments, the extraction of the large time window root mean square amplitude attribute of the pre-stack CRP gather Specific methods include:

[0124] (1) Pre-stack CRP gather Perform full stacking to obtain post-stack data based on direct stacking of pre-stack CRP gathers. ;

[0125] (2) Based on landmark geological strata For post-stack data based on direct stacking of pre-stack CRP gathers Large time window root mean square amplitude attributes are extracted to obtain the large time window root mean square amplitude attributes of the pre-stack CRP gather. ;

[0126] In some embodiments, the selection principle for the large time window is that the time window is not less than 200ms, and all data within the time window is valid across the entire work area;

[0127] In some embodiments, the specific method for obtaining the pre-stack CRP gather tectonic energy correction factor by applying the large-window root-mean-square amplitude attribute of the post-stack seismic data and the large-window root-mean-square amplitude attribute of the pre-stack CRP gather includes:

[0128] The pre-stack CRP gather structural energy correction factor = large-window root mean square amplitude attribute of post-stack seismic data / large-window root mean square amplitude attribute of pre-stack CRP gather, i.e.

[0129] ;

[0130] S3: Corrected pre-stack CRP gather Calculation: Energy correction factor constructed using pre-stack CRP gathers and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. .

[0131] This application provides a pre-stack gather energy equalization processing method. Under complex and steep geological conditions, pre-stack CRP gathers exhibit systematic energy differences related to the structure. Without requiring reprocessing, this method proposes a way to correct the energy of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and more accurately reflecting the true elasticity, physical properties, and fluidity of the reservoir.

[0132] Example 7

[0133] Based on the foregoing embodiments, this application further provides a pre-stack gather energy equalization processing method, including:

[0134] S1: Data Preparation: Collecting post-stack seismic data from known complex tectonic zones Pre-stack CRP collection and iconic geological strata ;

[0135] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0136] S2: Prestack CRP gather structure energy correction factor Calculation: Applying the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0137] In some embodiments, the application of the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor Specific methods include:

[0138] (1) Extracting post-stack seismic data Large time window root mean square amplitude property ;

[0139] (2) Extracting pre-stack CRP gathers Large time window root mean square amplitude property ;

[0140] (3) Apply the large time window root mean square amplitude attribute of the post-stack seismic data Large time window root mean square amplitude property of pre-stack CRP gather The pre-stack CRP gather construction energy correction factor was obtained. ;

[0141] In some embodiments, the specific method for extracting the large time window root mean square amplitude attribute of post-stack seismic data includes:

[0142] Based on key geological strata Extracting the large time window root mean square amplitude attribute of post-stack seismic data ;

[0143] In some embodiments, the extraction of the large time window root mean square amplitude attribute of the pre-stack CRP gather Specific methods include:

[0144] (1) Pre-stack CRP gather Perform full stacking to obtain post-stack data based on direct stacking of pre-stack CRP gathers. ;

[0145] (2) Based on landmark geological strata For post-stack data based on direct stacking of pre-stack CRP gathers Large time window root mean square amplitude attributes are extracted to obtain the large time window root mean square amplitude attributes of the pre-stack CRP gather. ;

[0146] In some embodiments, the selection principle for the large time window is that the time window is not less than 200ms, and all data within the time window is valid across the entire work area;

[0147] In some embodiments, the specific method for obtaining the pre-stack CRP gather tectonic energy correction factor by applying the large-window root-mean-square amplitude attribute of the post-stack seismic data and the large-window root-mean-square amplitude attribute of the pre-stack CRP gather includes:

[0148] The pre-stack CRP gather structural energy correction factor = large-window root mean square amplitude attribute of post-stack seismic data / large-window root mean square amplitude attribute of pre-stack CRP gather, i.e.

[0149] ;

[0150] S3: Corrected pre-stack CRP gather Calculation: Energy correction factor constructed using pre-stack CRP gathers and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. ;

[0151] In some embodiments, the specific method for constructing an energy correction factor and a pre-stack CRP gather using the pre-stack CRP gather to obtain the corrected pre-stack CRP gather includes:

[0152] Corrected pre-stack CRP gather = pre-stack CRP gather construction energy correction factor × pre-stack CRP gather, that is, .

[0153] For each CDP point Any path, any time Add gain to all CRP gather with energy correction .

[0154] This application provides a pre-stack gather energy equalization processing method. Under complex and steep geological conditions, pre-stack CRP gathers exhibit systematic energy differences related to the structure. Without requiring reprocessing, this method proposes a way to correct the energy of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and more accurately reflecting the true elasticity, physical properties, and fluidity of the reservoir.

[0155] Example 8

[0156] Based on the method in Example 7, this application provides an example of pre-stack gather energy consistency correction for an unconventional shale gas reservoir under complex geological conditions in southern China, using real data:

[0157] The target stratum in this work area, the Longmaxi Formation, is a marine shale with dramatic structural undulations, varying in depth from 500 meters to 4000 meters. Figure 2 This is a time-domain burial depth map of the Longmaxi Formation. Figure 3 The post-stack seismic profile obtained by stacking CRP gathers along the short axis clearly shows that the energy weakens in the steep reconstructed areas, and the overall energy of the profile is uneven. Figure 4 To perform full stacking based on pre-stack CRP gathers, and for post-stack data directly stacked based on pre-stack CRP gathers, a marker layer is used to extract the large-window root-mean-square amplitude attribute of the same time window. From Figure 4 It can be clearly seen that there is a good correlation between the planar amplitude properties and the burial depth of the marker layer Longmaxi Formation, indicating that the energy difference is affected by the complex structure.

[0158] Figure 5 To extract large-window root-mean-square amplitude attributes of post-stack data based on landmark geological strata, the same time window is used. As can be seen, the energy of the post-stack results is more balanced and reasonable, unaffected by the structure. Based on the extracted... Attributes and The energy correction factor is calculated using the properties and formula. :

[0159] ;

[0160] Figure 6 The obtained energy correction factor is a plane property.

[0161] Based on formula The corrected CRP gather can be calculated as follows: The gathers are then overlaid, and the root mean square amplitude attribute of the large time window is extracted based on the marker layer. ,like Figure 7 As shown, the energy-corrected gather has a more balanced and reasonable energy distribution, unaffected by structural features. Preprocessing this gather before performing pre-stack elasticity and property inversion yields more reasonable and reliable reservoir and fluid information.

[0162] Example 9

[0163] Based on the foregoing embodiments, this application provides a pre-stack gather energy equalization processing device, comprising:

[0164] Data preparation module, pre-stack CRP gather construction energy correction factor calculation module, and corrected pre-stack CRP gather calculation module;

[0165] The data preparation module collects post-stack seismic data from known complex tectonic zones. Pre-stack CRP collection and iconic geological strata ;

[0166] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0167] The pre-stack CRP gather tectonic energy correction factor calculation module applies the post-stack seismic data. Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0168] The corrected pre-stack CRP gather calculation module: constructs an energy correction factor using the pre-stack CRP gather. and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. .

[0169] This application provides a pre-stack gather energy equalization processing device. In complex and steep geological conditions, there are systematic energy differences in pre-stack CRP gathers related to the structure. Without the need for reprocessing, a method is proposed to achieve energy correction of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and further reflecting the true elasticity, physical properties and fluid properties of the reservoir more accurately.

[0170] Example 10

[0171] Based on the foregoing embodiments, this application provides a pre-stack gather energy equalization processing device, comprising:

[0172] Data preparation module, pre-stack CRP gather construction energy correction factor calculation module, and corrected pre-stack CRP gather calculation module;

[0173] The data preparation module collects post-stack seismic data from known complex tectonic zones. Pre-stack CRP collection and iconic geological strata ;

[0174] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0175] The pre-stack CRP gather tectonic energy correction factor calculation module applies the post-stack seismic data. Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0176] In some embodiments, the application of the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor Specific methods include:

[0177] (1) Extracting post-stack seismic data Large time window root mean square amplitude property ;

[0178] (2) Extracting pre-stack CRP gathers Large time window root mean square amplitude property ;

[0179] (3) Apply the large time window root mean square amplitude attribute of the post-stack seismic data Large time window root mean square amplitude property of pre-stack CRP gather The pre-stack CRP gather construction energy correction factor was obtained. ;

[0180] The corrected pre-stack CRP gather calculation module: constructs an energy correction factor using the pre-stack CRP gather. and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. .

[0181] This application provides a pre-stack gather energy equalization processing device. In complex and steep geological conditions, there are systematic energy differences in pre-stack CRP gathers related to the structure. Without the need for reprocessing, a method is proposed to achieve energy correction of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and further reflecting the true elasticity, physical properties and fluid properties of the reservoir more accurately.

[0182] Example 11

[0183] Based on the foregoing embodiments, this application provides a pre-stack gather energy equalization processing device, comprising:

[0184] Data preparation module, pre-stack CRP gather construction energy correction factor calculation module, and corrected pre-stack CRP gather calculation module;

[0185] The data preparation module collects post-stack seismic data from known complex tectonic zones. Pre-stack CRP collection and iconic geological strata ;

[0186] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0187] The pre-stack CRP gather tectonic energy correction factor calculation module applies the post-stack seismic data. Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0188] In some embodiments, the application of the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor Specific methods include:

[0189] (1) Extracting post-stack seismic data Large time window root mean square amplitude property ;

[0190] (2) Extracting pre-stack CRP gathers Large time window root mean square amplitude property ;

[0191] (3) Apply the large time window root mean square amplitude attribute of the post-stack seismic data Large time window root mean square amplitude property of pre-stack CRP gather The pre-stack CRP gather construction energy correction factor was obtained. ;

[0192] In some embodiments, the specific method for extracting the large time window root mean square amplitude attribute of post-stack seismic data includes:

[0193] Based on key geological strata Extracting the large time window root mean square amplitude attribute of post-stack seismic data ;

[0194] The corrected pre-stack CRP gather calculation module: constructs an energy correction factor using the pre-stack CRP gather. and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. .

[0195] This application provides a pre-stack gather energy equalization processing device. In complex and steep geological conditions, there are systematic energy differences in pre-stack CRP gathers related to the structure. Without the need for reprocessing, a method is proposed to achieve energy correction of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and further reflecting the true elasticity, physical properties and fluid properties of the reservoir more accurately.

[0196] Example 12

[0197] Based on the foregoing embodiments, this application provides a pre-stack gather energy equalization processing device, comprising:

[0198] Data preparation module, pre-stack CRP gather construction energy correction factor calculation module, and corrected pre-stack CRP gather calculation module;

[0199] The data preparation module collects post-stack seismic data from known complex tectonic zones. Pre-stack CRP collection and iconic geological strata ;

[0200] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0201] The pre-stack CRP gather tectonic energy correction factor calculation module applies the post-stack seismic data. Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0202] In some embodiments, the application of the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor Specific methods include:

[0203] (1) Extracting post-stack seismic data Large time window root mean square amplitude property ;

[0204] (2) Extracting pre-stack CRP gathers Large time window root mean square amplitude property ;

[0205] (3) Apply the large time window root mean square amplitude attribute of the post-stack seismic data Large time window root mean square amplitude property of pre-stack CRP gather The pre-stack CRP gather construction energy correction factor was obtained. ;

[0206] In some embodiments, the specific method for extracting the large time window root mean square amplitude attribute of post-stack seismic data includes:

[0207] Based on key geological strata Extracting the large time window root mean square amplitude attribute of post-stack seismic data ;

[0208] In some embodiments, the extraction of the large time window root mean square amplitude attribute of the pre-stack CRP gather Specific methods include:

[0209] (1) Pre-stack CRP gather Perform full stacking to obtain post-stack data based on direct stacking of pre-stack CRP gathers. ;

[0210] (2) Based on landmark geological strata For post-stack data based on direct stacking of pre-stack CRP gathers Large time window root mean square amplitude attributes are extracted to obtain the large time window root mean square amplitude attributes of the pre-stack CRP gather. ;

[0211] The corrected pre-stack CRP gather calculation module: constructs an energy correction factor using the pre-stack CRP gather. and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. .

[0212] This application provides a pre-stack gather energy equalization processing device. In complex and steep geological conditions, there are systematic energy differences in pre-stack CRP gathers related to the structure. Without the need for reprocessing, a method is proposed to achieve energy correction of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and further reflecting the true elasticity, physical properties and fluid properties of the reservoir more accurately.

[0213] Example 13

[0214] Based on the foregoing embodiments, this application provides a pre-stack gather energy equalization processing device, comprising:

[0215] Data preparation module, pre-stack CRP gather construction energy correction factor calculation module, and corrected pre-stack CRP gather calculation module;

[0216] The data preparation module collects post-stack seismic data from known complex tectonic zones. Pre-stack CRP collection and iconic geological strata ;

[0217] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0218] The pre-stack CRP gather tectonic energy correction factor calculation module applies the post-stack seismic data. Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0219] In some embodiments, the application of the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor Specific methods include:

[0220] (1) Extracting post-stack seismic data Large time window root mean square amplitude property ;

[0221] (2) Extracting pre-stack CRP gathers Large time window root mean square amplitude property ;

[0222] (3) Apply the large time window root mean square amplitude attribute of the post-stack seismic data Large time window root mean square amplitude property of pre-stack CRP gather The pre-stack CRP gather construction energy correction factor was obtained. ;

[0223] In some embodiments, the specific method for extracting the large time window root mean square amplitude attribute of post-stack seismic data includes:

[0224] Based on key geological strata Extracting the large time window root mean square amplitude attribute of post-stack seismic data ;

[0225] In some embodiments, the extraction of the large time window root mean square amplitude attribute of the pre-stack CRP gather Specific methods include:

[0226] (1) Pre-stack CRP gather Perform full stacking to obtain post-stack data based on direct stacking of pre-stack CRP gathers. ;

[0227] (2) Based on landmark geological strata For post-stack data based on direct stacking of pre-stack CRP gathers Large time window root mean square amplitude attributes are extracted to obtain the large time window root mean square amplitude attributes of the pre-stack CRP gather. ;

[0228] In some embodiments, the selection principle for the large time window is that the time window is not less than 200ms, and all data within the time window is valid across the entire work area;

[0229] The corrected pre-stack CRP gather calculation module: constructs an energy correction factor using the pre-stack CRP gather. and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. .

[0230] This application provides a pre-stack gather energy equalization processing device. In complex and steep geological conditions, there are systematic energy differences in pre-stack CRP gathers related to the structure. Without the need for reprocessing, a method is proposed to achieve energy correction of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and further reflecting the true elasticity, physical properties and fluid properties of the reservoir more accurately.

[0231] Example 14

[0232] Based on the foregoing embodiments, this application provides a pre-stack gather energy equalization processing device, comprising:

[0233] Data preparation module, pre-stack CRP gather construction energy correction factor calculation module, and corrected pre-stack CRP gather calculation module;

[0234] The data preparation module collects post-stack seismic data from known complex tectonic zones. Pre-stack CRP collection and iconic geological strata ;

[0235] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0236] The pre-stack CRP gather tectonic energy correction factor calculation module applies the post-stack seismic data. Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0237] In some embodiments, the application of the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor Specific methods include:

[0238] (1) Extracting post-stack seismic data Large time window root mean square amplitude property ;

[0239] (2) Extracting pre-stack CRP gathers Large time window root mean square amplitude property ;

[0240] (3) Apply the large time window root mean square amplitude attribute of the post-stack seismic data Large time window root mean square amplitude property of pre-stack CRP gather The pre-stack CRP gather construction energy correction factor was obtained. ;

[0241] In some embodiments, the specific method for extracting the large time window root mean square amplitude attribute of post-stack seismic data includes:

[0242] Based on key geological strata Extracting the large time window root mean square amplitude attribute of post-stack seismic data ;

[0243] In some embodiments, the extraction of the large time window root mean square amplitude attribute of the pre-stack CRP gather Specific methods include:

[0244] (1) Pre-stack CRP gather Perform full stacking to obtain post-stack data based on direct stacking of pre-stack CRP gathers. ;

[0245] (2) Based on landmark geological strata For post-stack data based on direct stacking of pre-stack CRP gathers Large time window root mean square amplitude attributes are extracted to obtain the large time window root mean square amplitude attributes of the pre-stack CRP gather. ;

[0246] In some embodiments, the selection principle for the large time window is that the time window is not less than 200ms, and all data within the time window is valid across the entire work area;

[0247] In some embodiments, the specific method for obtaining the pre-stack CRP gather tectonic energy correction factor by applying the large-window root-mean-square amplitude attribute of the post-stack seismic data and the large-window root-mean-square amplitude attribute of the pre-stack CRP gather includes:

[0248] The pre-stack CRP gather structural energy correction factor = large-window root mean square amplitude attribute of post-stack seismic data / large-window root mean square amplitude attribute of pre-stack CRP gather, i.e.

[0249] ;

[0250] The corrected pre-stack CRP gather calculation module: constructs an energy correction factor using the pre-stack CRP gather. and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. .

[0251] This application provides a pre-stack gather energy equalization processing device. In complex and steep geological conditions, there are systematic energy differences in pre-stack CRP gathers related to the structure. Without the need for reprocessing, a method is proposed to achieve energy correction of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and further reflecting the true elasticity, physical properties and fluid properties of the reservoir more accurately.

[0252] Example 15

[0253] Based on the foregoing embodiments, this application provides a pre-stack gather energy equalization processing device, comprising:

[0254] Data preparation module, pre-stack CRP gather construction energy correction factor calculation module, and corrected pre-stack CRP gather calculation module;

[0255] The data preparation module collects post-stack seismic data from known complex tectonic zones. Pre-stack CRP collection and iconic geological strata ;

[0256] in, The coordinates of the earthquake observation system are the x and y coordinates. For longitudinal time domain sampling of seismic data, It records the offset or incident angle information of the pre-stack CRP gather;

[0257] The pre-stack CRP gather tectonic energy correction factor calculation module applies the post-stack seismic data. Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor ;

[0258] In some embodiments, the application of the post-stack seismic data Pre-stack CRP collection and iconic geological strata Calculate the pre-stack CRP gather construction energy correction factor Specific methods include:

[0259] (1) Extracting post-stack seismic data Large time window root mean square amplitude property ;

[0260] (2) Extracting pre-stack CRP gathers Large time window root mean square amplitude property ;

[0261] (3) Apply the large time window root mean square amplitude attribute of the post-stack seismic data Large time window root mean square amplitude property of pre-stack CRP gather The pre-stack CRP gather construction energy correction factor was obtained. ;

[0262] In some embodiments, the specific method for extracting the large time window root mean square amplitude attribute of post-stack seismic data includes:

[0263] Based on key geological strata Extracting the large time window root mean square amplitude attribute of post-stack seismic data ;

[0264] In some embodiments, the extraction of the large time window root mean square amplitude attribute of the pre-stack CRP gather Specific methods include:

[0265] (1) Pre-stack CRP gather Perform full stacking to obtain post-stack data based on direct stacking of pre-stack CRP gathers. ;

[0266] (2) Based on landmark geological strata For post-stack data based on direct stacking of pre-stack CRP gathers Large time window root mean square amplitude attributes are extracted to obtain the large time window root mean square amplitude attributes of the pre-stack CRP gather. ;

[0267] In some embodiments, the selection principle for the large time window is that the time window is not less than 200ms, and all data within the time window is valid across the entire work area;

[0268] In some embodiments, the specific method for obtaining the pre-stack CRP gather tectonic energy correction factor by applying the large-window root-mean-square amplitude attribute of the post-stack seismic data and the large-window root-mean-square amplitude attribute of the pre-stack CRP gather includes:

[0269] The pre-stack CRP gather structural energy correction factor = large-window root mean square amplitude attribute of post-stack seismic data / large-window root mean square amplitude attribute of pre-stack CRP gather, i.e.

[0270] ;

[0271] The corrected pre-stack CRP gather calculation module: constructs an energy correction factor using the pre-stack CRP gather. and pre-stack CRP collection The corrected pre-stack CRP gather was obtained. ;

[0272] In some embodiments, the specific method for constructing an energy correction factor and a pre-stack CRP gather using the pre-stack CRP gather to obtain the corrected pre-stack CRP gather includes:

[0273] Corrected pre-stack CRP gather = pre-stack CRP gather construction energy correction factor × pre-stack CRP gather, that is, .

[0274] For each CDP point Any path, any time Add gain to all CRP gather with energy correction .

[0275] This application provides a pre-stack gather energy equalization processing device. In complex and steep geological conditions, there are systematic energy differences in pre-stack CRP gathers related to the structure. Without the need for reprocessing, a method is proposed to achieve energy correction of pre-stack gathers, thereby improving the accuracy of pre-stack seismic inversion and further reflecting the true elasticity, physical properties and fluid properties of the reservoir more accurately.

[0276] It should be noted that, in the embodiments of this application, if the above-described pre-stack gather energy equalization processing method is implemented as a software functional module and sold or used as an independent product, it can also be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the embodiments of this application, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, mobile hard drives, read-only memory (ROM), magnetic disks, or optical disks. Thus, the embodiments of this application are not limited to any specific hardware and software combination.

[0277] Accordingly, this application provides a storage medium storing a computer program thereon, characterized in that the computer program, when executed by a processor, implements the steps in the pre-stack gather energy equalization processing method provided in the above embodiments.

[0278] Example 16

[0279] This application provides a pre-stack gather energy equalization processing device including a memory and a processor. The memory stores a computer program. When the computer program is executed by the processor, the processor is configured to execute a program of a pre-stack gather energy equalization processing method stored in the memory, so as to implement the steps in the pre-stack gather energy equalization processing method provided in the above embodiment.

[0280] The descriptions of the display device and storage medium embodiments above are similar to those of the method embodiments above, and have similar beneficial effects. For technical details not disclosed in the computer device and storage medium embodiments of this application, please refer to the descriptions of the method embodiments of this application for understanding.

[0281] It should be noted that the descriptions of the storage medium and device embodiments above are similar to the descriptions of the method embodiments above, and have similar beneficial effects. For technical details not disclosed in the storage medium and device embodiments of this application, please refer to the descriptions of the method embodiments of this application for understanding.

[0282] It should be understood that the phrase "one embodiment" or "an embodiment" throughout the specification means that a specific feature, structure, or characteristic related to the embodiment is included in at least one embodiment of this application. Therefore, "in one embodiment" or "in an embodiment" appearing throughout the specification does not necessarily refer to the same embodiment. Furthermore, these specific features, structures, or characteristics can be combined in any suitable manner in one or more embodiments. It should be understood that in the various embodiments of this application, the sequence numbers of the above-described processes do not imply a sequential order of execution; the execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of this application. The sequence numbers of the above-described embodiments are merely descriptive and do not represent the superiority or inferiority of the embodiments.

[0283] It should be noted that, in this document, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes that element.

[0284] In the several embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. The device embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and in actual implementation, there may be other division methods, such as: multiple units or components can be combined, or integrated into another system, or some features can be ignored or not executed. In addition, the coupling, direct coupling, or communication connection between the various components shown or discussed can be through some interfaces, and the indirect coupling or communication connection between devices or units can be electrical, mechanical, or other forms.

[0285] The units described above as separate components may or may not be physically separate. The components shown as units may or may not be physical units. They may be located in one place or distributed across multiple network units. Some or all of the units may be selected to achieve the purpose of this embodiment according to actual needs.

[0286] In addition, each functional unit in the various embodiments of this application can be integrated into one processing unit, or each unit can be a separate unit, or two or more units can be integrated into one unit; the integrated unit can be implemented in hardware or in the form of hardware plus software functional units.

[0287] Those skilled in the art will understand that all or part of the steps of the above method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When the program is executed, it performs the steps of the above method embodiments. The aforementioned storage medium includes various media that can store program code, such as mobile storage devices, read-only memory (ROM), magnetic disks, or optical disks.

[0288] Alternatively, if the integrated units described above are implemented as software functional modules and sold or used as independent products, they can also be stored in a computer-readable storage medium. Based on this understanding, the technical solutions of the embodiments of this application, or the parts that contribute to the prior art, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a controller to execute all or part of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as mobile storage devices, ROMs, magnetic disks, or optical disks.

[0289] The above description is merely an embodiment of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for pre-stack gather energy equalization, characterized in that, include: S1: Data preparation: Collect post-stack seismic data, pre-stack CRP gathers, and key geological horizons for known complex tectonic zones; S2: Calculation of pre-stack CRP gather tectonic energy correction factor: The pre-stack CRP gather tectonic energy correction factor is calculated using the post-stack seismic data, pre-stack CRP gathers and key geological horizons. S3: Calculation of the corrected pre-stack CRP gather: The energy correction factor and the pre-stack CRP gather are constructed by applying the pre-stack CRP gather to obtain the corrected pre-stack CRP gather.

2. The method according to claim 1, characterized in that, The specific methods for calculating the tectonic energy correction factor of the pre-stack CRP gather using the post-stack seismic data, pre-stack CRP gathers, and key geological horizons include: (1) Extract the root mean square amplitude attribute of post-stack seismic data with a large time window; (2) Extract the large time window root mean square amplitude attribute of the pre-stack CRP gather; (3) The large time window root mean square amplitude attribute of the post-stack seismic data and the large time window root mean square amplitude attribute of the pre-stack CRP gather are used to obtain the structural energy correction factor of the pre-stack CRP gather.

3. The method according to claim 2, characterized in that, The specific method for extracting the large time-window root mean square amplitude attribute of post-stack seismic data includes: Based on key geological strata, large time window root mean square amplitude attributes of post-stack seismic data are extracted.

4. The method according to claim 3, characterized in that, The specific method for extracting the large time-window root mean square amplitude attribute of the pre-stack CRP gather includes: (1) Perform full stacking of pre-stack CRP gathers to obtain post-stack data based on direct stacking of pre-stack CRP gathers; (2) Based on the landmark geological strata, the large time window root mean square amplitude attribute is extracted from the post-stack data directly superimposed from the pre-stack CRP gathers to obtain the large time window root mean square amplitude attribute of the pre-stack CRP gathers.

5. The method according to claim 4, characterized in that, The selection principle for the large time window is that the time window is not less than 200ms, and all data within the time window is valid within the entire work area.

6. The method according to claim 5, characterized in that, The specific method for obtaining the structural energy correction factor of the pre-stack CRP gather by applying the large-window root-mean-square amplitude attribute of the post-stack seismic data and the large-window root-mean-square amplitude attribute of the pre-stack CRP gather includes: The structural energy correction factor of the pre-stack CRP gather = the root mean square amplitude attribute of the large time window of the post-stack seismic data / the root mean square amplitude attribute of the large time window of the pre-stack CRP gather.

7. The method according to claim 6, characterized in that, The specific method for constructing the energy correction factor and the pre-stack CRP gather using the pre-stack CRP gather, and obtaining the corrected pre-stack CRP gather, includes: Corrected pre-stack CRP gather = pre-stack CRP gather construction energy correction factor × pre-stack CRP gather.

8. A pre-stack gather energy equalization processing device, characterized in that, include: Data preparation module, pre-stack CRP gather construction energy correction factor calculation module, and corrected pre-stack CRP gather calculation module; The data preparation module collects post-stack seismic data, pre-stack CRP gathers, and key geological strata in known complex tectonic zones. The pre-stack CRP gather tectonic energy correction factor calculation module: calculates the pre-stack CRP gather tectonic energy correction factor using the post-stack seismic data, pre-stack CRP gathers and key geological strata; The corrected pre-stack CRP gather calculation module: applies the pre-stack CRP gather to construct the energy correction factor and the pre-stack CRP gather, and obtains the corrected pre-stack CRP gather.

9. A pre-stack gatherer energy equalization processing device, characterized in that, It includes a memory and a processor, wherein the memory stores a computer program that, when executed by the processor, performs the pre-stack gather energy equalization processing method as described in any one of claims 1 to 7.

10. A storage medium, characterized in that, The computer program stored in the storage medium can be executed by one or more processors and can be used to implement the pre-stack gather energy equalization processing method as described in any one of claims 1 to 7.

Images