Multiple identification and attenuation method and device, electronic equipment and medium
By performing dynamic and reaction correction of multiple waves in the hierarchical geological model, combined with Ladong transformation and multi-channel convolution, the problem of multiple wave identification and removal in physical simulated seismic data is solved, and imaging quality and computing efficiency are improved.
Patent Information
- Application Number
- CN202311648817.5
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2023-12-04
- Publication Date
- 2025-06-06
AI Technical Summary
In oil and gas geophysical exploration, it is difficult for the prior art to effectively identify and remove multiple waves in physically simulated seismic data, resulting in low imaging quality and signal leakage problems.
By selecting multiple coordinate points in the hierarchical geological model, time-depth correction of the hierarchical well in the time domain, pressing down multiple waves underwater, and using the multiple wave velocity to perform dynamic correction and reaction correction, combining Ladong transformation and multi-channel convolution, multiple waves are gradually removed.
It improves the imaging quality of physical simulation data, reduces signal leakage, enhances the reliability of imaging mechanisms, and improves computing efficiency.
Smart Images
Figure CN120103423A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the field of oil and gas geophysical exploration, and more specifically, to a multiple wave identification and attenuation method, device, electronic equipment and medium. Background Art
[0002] In actual seismic data processing, identifying and removing multiple waves has always been a difficult problem. In physical simulation, the multiple waves of seismic data are relatively developed, with strong and concentrated energy. The main types are water bottom-related long-range multiple waves, interlayer short-range multiple waves, and single multiple waves and multiple multiple waves of these two types.
[0003] Taking advantage of the differences in properties between multiple waves and primary waves such as apparent velocity, inclination, frequency and periodicity, the current methods for removing multiple waves mainly include filtering methods (predictive deconvolution, τp transform, Radon transform, etc.) and predictive subtraction methods (SRME free surface multiple wave prediction suppression, etc.). However, each method has certain defects. The predictive deconvolution method is effective for short-period multiple waves in shallow layers. The classical Radon transform requires that the time difference between the primary wave and the multiple waves is greater than 30ms, and the primary wave velocity needs to be accurately picked up, otherwise the removed multiple waves contain valid information, resulting in effective signal leakage. The predictive subtraction method based on wave equation theory can attenuate the effect, but it is necessary to give the layer or range that generates the multiple waves, otherwise it is difficult to obtain the ideal effect, and this type of method has a relatively large amount of calculation and efficiency problems.
[0004] There is a need to develop a multiple wave identification and attenuation method for physically simulated seismic data.
[0005] The information disclosed in the background technology section of the present invention is only intended to deepen the understanding of the general background technology of the present invention, and should not be regarded as acknowledging or suggesting in any form that the information constitutes the prior art already known to those skilled in the art. Summary of the invention
[0006] The present invention proposes a multiple wave identification and attenuation method, device, electronic equipment and medium, which can effectively identify multiple waves, improve the imaging quality of physical simulation data and promote the reliability of revealing the imaging mechanism.
[0007] In a first aspect, an embodiment of the present disclosure provides a multiple wave identification and attenuation method, comprising:
[0008] Select multiple coordinate points in the stratigraphic geological model to obtain the corresponding stratigraphic wells in the time domain, and suppress the multiple waves related to the water bottom in the shot gather domain;
[0009] According to the cmp gathers, the attenuated gathers are obtained;
[0010] Performing a counter-correction on the attenuated gathers by using the multiple wave velocities, and performing a dynamic correction by using the primary wave velocity, to obtain a gather after dynamic correction and flattening;
[0011] Attenuating the multiple waves based on the dynamic correction and flattened gathers, subtracting the result from the data of the cmp gathers to obtain the multiple wave gathers;
[0012] Obtaining multiple wave attenuation gathers according to the cmp gather data and the attenuated gather data;
[0013] The near-track multiple waves are attenuated for the gathers after dynamic correction and flattening to obtain the near-track multiple wave attenuation gathers;
[0014] Multi-channel convolution is applied to the near-channel multiple wave attenuation gather to obtain a cmp gather with multiple waves removed.
[0015] As a specific implementation of the embodiment of the present disclosure, selecting multiple coordinate points in the stratigraphic geological model to obtain the corresponding stratigraphic wells in the time domain includes:
[0016] Select multiple coordinate points in the stratigraphic geological model, use the coordinate point position-stratigraphic depth to make multiple wells, perform time-depth correction on the velocity of each layer corresponding to the velocity model for the well, and obtain the stratigraphic well in the time domain.
[0017] As a specific implementation of the embodiment of the present disclosure, obtaining the attenuated gathers according to the cmp gathers includes:
[0018] The data after removing the bottom multiple waves are sorted into the cmp domain, and the multiple wave velocities are picked up on the cmp gathers;
[0019] Performing dynamic correction on the cmp gather by using the multiple wave velocity to obtain an attenuated gather.
[0020] As a specific implementation manner of the embodiment of the present disclosure, a dynamic correction is performed based on the multiple wave velocity through Radon transform.
[0021] As a specific implementation method of the embodiment of the present disclosure, the multiple waves are attenuated by Radon transformation on the flattened gathers.
[0022] As a specific implementation manner of the embodiment of the present disclosure, the attenuated gather data is subtracted from the cmp gather data to obtain a multiple wave attenuation gather.
[0023] As a specific implementation method of the embodiment of the present disclosure, FK dip filtering or incision is performed on the track gather after dynamic correction and flattening to achieve attenuation of the residual near-track multiple waves.
[0024] In a second aspect, the embodiment of the present disclosure further provides a multiple wave identification and attenuation device, including:
[0025] The suppression module selects multiple coordinate points in the stratigraphic geological model, obtains the corresponding stratigraphic wells in the time domain, and suppresses the multiple waves related to the water bottom in the shot gather domain;
[0026] The sorting module obtains the attenuated gathers according to the cmp gathers;
[0027] A correction module, for performing a reaction correction on the attenuated gathers by using the multiple wave velocities, and performing a dynamic correction by using the primary wave velocity, to obtain a gather after dynamic correction and flattening;
[0028] A first attenuation module attenuates the multiple waves based on the dynamic correction and flattened gathers, and subtracts the result from the data of the cmp gathers to obtain multiple wave gathers;
[0029] A calculation module, which obtains a multiple wave attenuation gather according to the cmp gather data and the attenuated gather data;
[0030] The second attenuation module attenuates the near-channel multiple waves for the gathers after dynamic correction and flattening to obtain the near-channel multiple wave attenuation gathers;
[0031] The convolution calculation module applies multi-channel convolution to the near-channel multiple wave attenuation gather to obtain a cmp gather with multiple waves removed.
[0032] As a specific implementation of the embodiment of the present disclosure, selecting multiple coordinate points in the stratigraphic geological model to obtain the corresponding stratigraphic wells in the time domain includes:
[0033] Select multiple coordinate points in the stratigraphic geological model, use the coordinate point position-stratigraphic depth to make multiple wells, perform time-depth correction on the velocity of each layer corresponding to the velocity model for the well, and obtain the stratigraphic well in the time domain.
[0034] As a specific implementation of the embodiment of the present disclosure, obtaining the attenuated gathers according to the cmp gathers includes:
[0035] The data after removing the bottom multiple waves are sorted into the cmp domain, and the multiple wave velocities are picked up on the cmp gathers;
[0036] Performing dynamic correction on the cmp gather by using the multiple wave velocity to obtain an attenuated gather.
[0037] As a specific implementation manner of the embodiment of the present disclosure, a dynamic correction is performed based on the multiple wave velocity through Radon transform.
[0038] As a specific implementation method of the embodiment of the present disclosure, the multiple waves are attenuated by Radon transformation on the flattened gathers.
[0039] As a specific implementation manner of the embodiment of the present disclosure, the attenuated gather data is subtracted from the cmp gather data to obtain a multiple wave attenuation gather.
[0040] As a specific implementation method of the embodiment of the present disclosure, FK dip filtering or incision is performed on the track gather after dynamic correction and flattening to achieve attenuation of the residual near-track multiple waves.
[0041] In a third aspect, an embodiment of the present disclosure further provides an electronic device, the electronic device comprising:
[0042] A memory storing executable instructions;
[0043] A processor is used to execute the executable instructions in the memory to implement the multiple wave identification and attenuation method.
[0044] In a fourth aspect, an embodiment of the present disclosure further provides a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, the multiple wave identification and attenuation method is implemented.
[0045] Its beneficial effects are:
[0046] The present invention improves the classical high-precision Radon transform processing method to avoid the effective signal leakage problem caused by inaccurate speed, improves the operation efficiency and ensures the fidelity of the result. Since the model of physical simulation is basically known and the large layer is accurately described, the layer of the physical model is extracted while combining the above method, similar to the combination with well seismic, and the time-depth conversion is performed to the time domain using the marked layers at different coordinates, which is helpful to help identify which are multiple problems on the time gather.
[0047] The methods and apparatus of the present invention have other features and advantages that will be apparent from, or will be described in detail in, the accompanying drawings and subsequent detailed descriptions incorporated herein, which together serve to explain the specific principles of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS
[0048] The above and other objects, features and advantages of the present invention will become more apparent through a more detailed description of exemplary embodiments of the present invention in conjunction with the accompanying drawings, wherein like reference numerals generally represent like components throughout the exemplary embodiments of the present invention.
[0049] Figure 1 A flow chart showing the steps of a multiple identification and attenuation method according to one embodiment of the present invention.
[0050] Figure 2A block diagram of a multiple wave identification and attenuation device according to an embodiment of the present invention is shown.
[0051] Description of reference numerals:
[0052] 201, suppression module; 202, sorting module; 203, correction module; 204, first attenuation module; 205, calculation module; 206, second attenuation module; 207, convolution calculation module. DETAILED DESCRIPTION
[0053] The preferred embodiments of the present invention will be described in more detail below. Although the preferred embodiments of the present invention are described below, it should be understood that the present invention can be implemented in various forms and should not be limited to the embodiments set forth herein.
[0054] To facilitate understanding of the solutions and effects of the embodiments of the present invention, six specific application examples are given below. Those skilled in the art should understand that the examples are only for facilitating understanding of the present invention, and any specific details thereof are not intended to limit the present invention in any way.
[0055] Example 1
[0056] Figure 1 A flow chart showing the steps of a multiple identification and attenuation method according to one embodiment of the present invention.
[0057] like Figure 1 As shown, the multiple wave identification and attenuation method includes: step 101, selecting multiple coordinate points in the stratum geological model, obtaining the corresponding stratum wells in the time domain, and suppressing the multiple waves related to the water bottom in the shot gather domain; step 102, obtaining the attenuated gather according to the cmp gather; step 103, performing anti-correction correction by using the multiple wave velocity for the attenuated gather, and applying the primary wave velocity for dynamic correction to obtain the gather after dynamic correction and leveling; step 104, attenuating the multiple waves based on the gather after dynamic correction and leveling, and subtracting the result from the data of the cmp gather to obtain the multiple wave gather; step 105, obtaining the multiple wave attenuation gather according to the data of the cmp gather and the attenuated gather data; step 106, performing attenuation of near-track multiple waves for the gather after dynamic correction and leveling to obtain the near-track multiple wave attenuation gather; step 107, applying multi-channel convolution to the near-track multiple wave attenuation gather to obtain the cmp gather with the multiple waves removed.
[0058] In one example, multiple coordinate points are selected in the stratigraphic geological model to obtain corresponding stratigraphic wells in the time domain, including:
[0059] Select multiple coordinate points in the stratigraphic geological model, use the coordinate point position-stratigraphic depth to make multiple wells, perform time-depth correction on the velocity of each layer corresponding to the velocity model for the well, and obtain the stratigraphic well in the time domain.
[0060] In one example, obtaining the attenuated gathers according to the cmp gathers includes:
[0061] The data after removing the bottom multiple waves are sorted into the cmp domain, and the multiple wave velocities are picked up on the cmp gathers;
[0062] The cmp gathers are corrected by multiple wave velocities to obtain the attenuated gathers.
[0063] In one example, the NMO correction is performed based on the multiple wave velocity by Radon transform.
[0064] In one example, the multiple waves are attenuated by Radon transform on the flattened gathers.
[0065] In one example, the attenuated gather data is subtracted from the cmp gather data to obtain the multiple wave attenuation gather.
[0066] In one example, FK dip filtering or incision is performed on the track gather after NMO leveling to achieve attenuation of the residual near-track multiple waves.
[0067] Specifically, multiple coordinate points are selected in the stratigraphic geological model, and multiple wells are made using the coordinate point position-stratigraphic depth. The velocity model is used to perform time-depth correction on the velocity of each layer of the well to obtain the stratigraphic well in the time domain. The stratigraphic well is used to calibrate which axes in the time domain are valid. Multiple waves related to the water bottom are suppressed in the shot gather domain to remove multiple waves generated by the water layer.
[0068] According to the cmp data set, the attenuated data set is obtained; the data without the bottom multiple waves are sorted into the cmp domain, and the multiple wave velocities are picked up on the cmp data set; the cmp data set is subjected to dynamic correction based on the multiple wave velocities through Radon transform to obtain the attenuated data set, and only the multiple waves are left after the effective waves are removed to protect the primary wave.
[0069] For the attenuated data set, the multiple wave velocity is used for counter-correction, and the primary wave velocity is used for dynamic correction to obtain the data set after dynamic correction and flattening. The multiple waves are flattened by dynamic correction using the multiple wave velocity, and then the data set is obtained by correction using the primary wave velocity.
[0070] The multiple waves are attenuated by Radon transform on the flattened channel gathers, and the result is subtracted from the data of cmp channel gathers to obtain the multiple wave channel gathers.
[0071] Subtract the attenuated gather data from the cmp gather data to obtain the multiple wave attenuation gather, thus protecting the primary wave.
[0072] Perform FK dip filtering or incision on the track gather after dynamic correction and leveling to attenuate the residual near-track multiple waves and obtain near-track multiple wave attenuation track gathers. Since the near-track multiple waves and primary wave corrections can be nearly leveled, the near-track multiple waves will enter the effective waves, and these multiple waves will be removed in this step.
[0073] Multi-channel convolution is applied to the near-channel multiple wave attenuation gathers to obtain the cmp gathers with multiple waves removed. The prediction step size and operator are determined according to the effective fruit wave length, which has a compression effect on some inter-layer multiple wave side lobes.
[0074] Example 2
[0075] The present invention also provides a multiple wave identification and attenuation device, comprising:
[0076] The suppression module selects multiple coordinate points in the stratigraphic geological model, obtains the corresponding stratigraphic wells in the time domain, and suppresses the multiple waves related to the water bottom in the shot gather domain;
[0077] The sorting module obtains the attenuated gathers according to the cmp gathers;
[0078] The correction module performs a reaction correction on the attenuated gathers by using multiple wave velocities, and applies a primary wave velocity to perform a dynamic correction to obtain the gathers after dynamic correction and flattening;
[0079] The first attenuation module attenuates the multiple waves based on the gathers after NMO and subtracts the result from the data of cmp gathers to obtain the multiple wave gathers;
[0080] The calculation module obtains the multiple wave attenuation gathers according to the cmp gather data and the attenuated gather data;
[0081] The second attenuation module attenuates the near-channel multiple waves for the gathers after dynamic correction and flattening to obtain the near-channel multiple wave attenuation gathers;
[0082] The convolution calculation module applies multi-channel convolution to the near-channel multiple wave attenuation gathers to obtain the cmp gathers with multiple waves removed.
[0083] In one example, multiple coordinate points are selected in the stratigraphic geological model to obtain corresponding stratigraphic wells in the time domain, including:
[0084] Select multiple coordinate points in the stratigraphic geological model, use the coordinate point position-stratigraphic depth to make multiple wells, perform time-depth correction on the velocity of each layer corresponding to the velocity model for the well, and obtain the stratigraphic well in the time domain.
[0085] In one example, obtaining the attenuated gathers according to the cmp gathers includes:
[0086] The data after removing the bottom multiple waves are sorted into the cmp domain, and the multiple wave velocities are picked up on the cmp gathers;
[0087] The cmp gathers are corrected by multiple wave velocities to obtain the attenuated gathers.
[0088] In one example, the NMO correction is performed based on the multiple wave velocity by Radon transform.
[0089] In one example, the multiple waves are attenuated by Radon transform on the flattened gathers.
[0090] In one example, the attenuated gather data is subtracted from the cmp gather data to obtain the multiple wave attenuation gather.
[0091] In one example, FK dip filtering or incision is performed on the track gather after NMO leveling to achieve attenuation of the residual near-track multiple waves.
[0092] Specifically, multiple coordinate points are selected in the stratigraphic geological model, and multiple wells are made using the coordinate point position-stratigraphic depth. The velocity model is used to perform time-depth correction on the velocity of each layer of the well to obtain the stratigraphic well in the time domain. The stratigraphic well is used to calibrate which axes in the time domain are valid. Multiple waves related to the water bottom are suppressed in the shot gather domain to remove multiple waves generated by the water layer.
[0093] According to the cmp data set, the attenuated data set is obtained; the data without the bottom multiple waves are sorted into the cmp domain, and the multiple wave velocities are picked up on the cmp data set; the cmp data set is subjected to dynamic correction based on the multiple wave velocities through Radon transform to obtain the attenuated data set, and only the multiple waves are left after the effective waves are removed to protect the primary wave.
[0094] For the attenuated data set, the multiple wave velocity is used for counter-correction, and the primary wave velocity is used for dynamic correction to obtain the data set after dynamic correction and flattening. The multiple waves are flattened by dynamic correction using the multiple wave velocity, and then the data set is obtained by correction using the primary wave velocity.
[0095] The multiple waves are attenuated by Radon transform on the flattened channel gathers, and the result is subtracted from the data of cmp channel gathers to obtain the multiple wave channel gathers.
[0096] Subtract the attenuated gather data from the cmp gather data to obtain the multiple wave attenuation gather, thus protecting the primary wave.
[0097] Perform FK dip filtering or incision on the track gather after dynamic correction and leveling to attenuate the residual near-track multiple waves and obtain near-track multiple wave attenuation track gathers. Since the near-track multiple waves and primary wave corrections can be nearly leveled, the near-track multiple waves will enter the effective waves, and these multiple waves will be removed in this step.
[0098] Multi-channel convolution is applied to the near-channel multiple wave attenuation gathers to obtain the cmp gathers with multiple waves removed. The prediction step size and operator are determined according to the effective fruit wave length, which has a compression effect on some inter-layer multiple wave side lobes.
[0099] Example 3
[0100] (1) Five coordinate points are selected in the stratigraphic geological model for designing the physical model, which can represent the stratigraphic layers and their changes in the entire model. Five "wells" are made using the coordinate point position-stratigraphic depth. The velocity model is used to perform time-depth correction on the wells to obtain stratigraphic wells in the time domain. The wells can be projected onto the time data set and are an accurate reference for identifying multiple waves in the time domain data set.
[0101] (2) Using the conventional wave equation extension method to suppress the multiple waves related to the water bottom in the shot gather domain;
[0102] (3) Data preparation: sort the data into the cmp domain, implement the method in the common center point or common reflection point domain, and pick up the multiple wave velocity on the cmp gather;
[0103] (4) Apply multiple wave velocity to the gathers for dynamic correction, select appropriate dynamic correction time difference parameters, and use high-precision Radon transform to attenuate the effective waves;
[0104] (5) Apply the multiple wave velocity to the gathers after the attenuation of the "multiple waves" for the anti-motion correction calculation, and apply the primary wave velocity to the dynamic correction calculation;
[0105] (6) Applying high-precision Radon transform to the corrected channel gathers to attenuate the multiple waves, and subtracting the result from the data before Radon transform to obtain the multiple wave channel gathers;
[0106] (7) Subtract the multiple wave gathers in step (4) from the initial input data to obtain the multiple wave attenuation gathers;
[0107] (8) For the gathers after multiple wave removal in (5), apply FK dip filtering (or incision) to attenuate the residual near-track multiple waves;
[0108] (9) Finally, multi-channel prediction deconvolution is used to suppress the inter-layer multiple times to obtain the cmp data set with multiple waves removed.
[0109] Example 4
[0110] Figure 2 A block diagram of a multiple wave identification and attenuation device according to an embodiment of the present invention is shown.
[0111] like Figure 2 As shown, the multiple wave identification and attenuation device comprises:
[0112] The suppression module 201 selects a plurality of coordinate points in the stratigraphic geological model, obtains the corresponding stratigraphic wells in the time domain, and suppresses the multiple waves related to the water bottom in the shot gather domain;
[0113] A sorting module 202 obtains attenuated gathers according to the cmp gathers;
[0114] The correction module 203 performs a counter-correction on the attenuated gathers by using multiple wave velocities, and applies a primary wave velocity to perform a dynamic correction to obtain a gather after dynamic correction and flattening;
[0115] The first attenuation module 204 attenuates the multiple waves based on the NMO-flattened gathers, and subtracts the result from the cmp gather data to obtain the multiple wave gathers;
[0116] The calculation module 205 obtains the multiple wave attenuation gathers according to the cmp gather data and the attenuated gather data;
[0117] The second attenuation module 206 performs attenuation of near-channel multiple waves on the gathers after the dynamic correction and flattening to obtain near-channel multiple wave attenuation gathers;
[0118] The convolution calculation module 207 applies multi-channel convolution to the near-channel multiple wave attenuation gathers to obtain cmp gathers with multiple waves removed.
[0119] As an optional solution, multiple coordinate points are selected in the stratigraphic geological model to obtain the corresponding stratigraphic wells in the time domain, including:
[0120] Select multiple coordinate points in the stratigraphic geological model, use the coordinate point position-stratigraphic depth to make multiple wells, perform time-depth correction on the velocity of each layer corresponding to the velocity model for the well, and obtain the stratigraphic well in the time domain.
[0121] As an optional solution, based on the cmp gathers, the attenuated gathers are obtained including:
[0122] The data after removing the bottom multiple waves are sorted into the cmp domain, and the multiple wave velocities are picked up on the cmp gathers;
[0123] The cmp gathers are corrected by multiple wave velocities to obtain the attenuated gathers.
[0124] As an alternative, a normal moveout correction is performed based on the multiple wave velocities by Radon transform.
[0125] As an optional solution, the multiple waves can be attenuated by Radon transform on the flattened gathers.
[0126] As an optional solution, the attenuated gather data is subtracted from the cmp gather data to obtain the multiple wave attenuation gather data.
[0127] As an optional solution, FK dip filtering or incision can be performed on the track gathers after dynamic correction and leveling to achieve attenuation of the residual near-track multiple waves.
[0128] Example 5
[0129] This embodiment provides an electronic device, which includes: a memory storing executable instructions; and a processor, which runs the executable instructions in the memory to implement the above-mentioned multiple wave identification and attenuation method.
[0130] An electronic device according to an embodiment of the present disclosure includes a memory and a processor.
[0131] The memory is used to store non-temporary computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory (cache), etc. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.
[0132] The processor may be a central processing unit (CPU) or other forms of processing units having data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of the present disclosure, the processor is used to run the computer-readable instructions stored in the memory.
[0133] Those skilled in the art should be able to understand that in order to solve the technical problem of how to obtain a good user experience, the present embodiment may also include well-known structures such as a communication bus and an interface, and these well-known structures should also be included in the protection scope of the present disclosure.
[0134] For detailed description of this embodiment, reference may be made to the corresponding descriptions in the aforementioned embodiments, which will not be repeated here.
[0135] Example 6
[0136] This embodiment provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the multiple wave identification and attenuation method is implemented.
[0137] According to the computer-readable storage medium of the embodiment of the present disclosure, non-transitory computer-readable instructions are stored thereon. When the non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the above-mentioned methods of each embodiment of the present disclosure are executed.
[0138] The above-mentioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tapes or mobile hard disks), media with built-in rewritable non-volatile memory (e.g., memory cards), and media with built-in ROM (e.g., ROM boxes).
[0139] Those skilled in the art should understand that the purpose of the above description of the embodiments of the present invention is only to exemplarily illustrate the beneficial effects of the embodiments of the present invention, and is not intended to limit the embodiments of the present invention to any given examples.
[0140] The embodiments of the present invention have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.
Claims
1. A multiple wave identification and attenuation method, It is characterized in that include: Select multiple coordinate points in the stratigraphic geological model to obtain the corresponding stratigraphic wells in the time domain, and suppress the multiple waves related to the water bottom in the shot gather domain; According to the cmp gathers, the attenuated gathers are obtained; Performing a counter-correction on the attenuated gathers by using the multiple wave velocities, and performing a dynamic correction by using the primary wave velocity, to obtain a gather after dynamic correction and flattening; Attenuating the multiple waves based on the dynamic correction and flattened gathers, subtracting the result from the data of the cmp gathers to obtain the multiple wave gathers; Obtaining multiple wave attenuation gathers according to the cmp gather data and the attenuated gather data; The near-track multiple waves are attenuated for the gathers after dynamic correction and flattening to obtain the near-track multiple wave attenuation gathers; Multi-channel convolution is applied to the near-channel multiple wave attenuation gather to obtain a cmp gather with multiple waves removed.
2. The multiple wave identification and attenuation method according to claim 1, in, Select multiple coordinate points in the stratigraphic geological model to obtain the corresponding stratigraphic wells in the time domain, including: Select multiple coordinate points in the stratigraphic geological model, use the coordinate point position-stratigraphic depth to make multiple wells, perform time-depth correction on the velocity of each layer corresponding to the velocity model for the well, and obtain the stratigraphic well in the time domain.
3. The multiple wave identification and attenuation method according to claim 1, in, According to the cmp gathers, the attenuated gathers include: The data after removing the bottom multiple waves are sorted into the cmp domain, and the multiple wave velocities are picked up on the cmp gathers; Performing dynamic correction on the cmp gather by using the multiple wave velocity to obtain an attenuated gather.
4. The multiple wave identification and attenuation method according to claim 3, in, A dynamic correction is performed based on the multiple wave velocity by Radon transformation.
5. The multiple wave identification and attenuation method according to claim 1, in, The multiple waves are attenuated by using Radon transform on the flattened gathers.
6. The multiple wave identification and attenuation method according to claim 1, in, The attenuated gather data is subtracted from the cmp gather data to obtain the multiple wave attenuation gather.
7. The multiple wave identification and attenuation method according to claim 1, in, Perform FK dip filtering or incision on the gathers after dynamic correction and leveling to achieve attenuation of the residual near-track multiple waves.
8. A multiple wave identification and attenuation device, It is characterized in that include: The suppression module selects multiple coordinate points in the stratigraphic geological model, obtains the corresponding stratigraphic wells in the time domain, and suppresses the multiple waves related to the water bottom in the shot gather domain; The sorting module obtains the attenuated gathers according to the cmp gathers; A correction module, for performing a reaction correction on the attenuated gathers by using the multiple wave velocities, and performing a dynamic correction by using the primary wave velocity, to obtain a gather after dynamic correction and flattening; A first attenuation module attenuates the multiple waves based on the dynamic correction and flattened gathers, and subtracts the result from the data of the cmp gathers to obtain multiple wave gathers; A calculation module, which obtains a multiple wave attenuation gather according to the cmp gather data and the attenuated gather data; The second attenuation module attenuates the near-channel multiple waves for the gathers after dynamic correction and flattening to obtain the near-channel multiple wave attenuation gathers; The convolution calculation module applies multi-channel convolution to the near-channel multiple wave attenuation gather to obtain a cmp gather with multiple waves removed.
9. An electronic device, It is characterized in that The electronic device comprises: A memory storing executable instructions; A processor, wherein the processor runs the executable instructions in the memory to implement the multiple wave identification and attenuation method according to any one of claims 1 to 7.
10. A computer-readable storage medium, It is characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the multiple wave identification and attenuation method according to any one of claims 1 to 7 is implemented.