Surface wave suppression method and device based on energy replacement and electronic equipment

Abstract

The invention relates to the technical field of oil and gas geophysics, and discloses a surface wave suppression method and device based on energy replacement and electronic equipment. The surface wave suppression method based on energy replacement comprises the following steps: determining a bandwidth range of surface wave interference; separating the original data to obtain data in a surface wave bandwidth range and data outside the surface wave bandwidth range; performing multi-time window gain analysis, and calculating a gain operator of each seismic channel in each time window; performing gain operator replacement; the large gain operator after replacement is used for carrying out anti-gain calculation on the single-shot signal in the frequency width range of the surface wave to realize suppression of the surface wave; and adding the data outside the surface wave bandwidth range separated from the original data and the data in the surface wave bandwidth range after energy replacement to obtain single-shot seismic data after surface wave suppression. According to the technical scheme, surface waves can be effectively suppressed, low-frequency effective components of signals are basically not affected, the amplitude preservation performance is good, the signal-to-noise ratio of data is improved, and effective reflection information does not exist in residual errors.

Description

Technical Field

[0001] The invention relates to the technical field of oil and gas geophysics, and in particular to a surface wave suppression method, device and electronic equipment based on energy replacement. Background Art

[0002] Surface waves are interference waves that propagate near the surface and are excited by the earthquake source in onshore seismic exploration. The energy of surface waves in seismic records is generally several times or even dozens of times that of effective waves, and is intertwined with the phase axis of the reflection wave, which seriously reduces the signal-to-noise ratio of seismic records, has a great impact on the effective use of pre-stack amplitude information, and reduces the quality of seismic data processing. Effective suppression of surface waves is conducive to improving the quality of seismic data and improving the signal-to-noise ratio of seismic data. It is crucial for subsequent pre-stack attribute analysis, full waveform velocity inversion, and improving the imaging quality of deep geological targets. Surface waves are a kind of regular interference that is widely present in seismic exploration data acquisition. They are linearly distributed on the shot gather and have the characteristics of high intensity, low frequency, low velocity, and long vibration duration. They seriously affect the signal-to-noise ratio and resolution of data, especially for fracture-cavity reservoirs. Surface waves are mixed with low-frequency diffraction waves. Poor suppression of surface waves will damage effective waves, thereby affecting the imaging of fractures and caves.

[0003] In the prior art, there are many methods for suppressing surface waves:

[0004] The regional abnormal noise attenuation commonly used in production is to set a window in the surface wave area and suppress the noise within the window. This method filters out the surface wave according to the surface wave frequency, which will lose the low-frequency effective information.

[0005] Another technology that is widely used in production is the noise suppression technology of the analog surface roll adaptive phase subtraction detection line domain (SWAMI). This technology can suppress some noise in the surface roll, but some low-frequency noise will remain. If the surface roll noise is not removed cleanly, it will cause arcing in the later stage.

[0006] The coherent noise subtraction method (NUCNS) can remove the linear noise in the surface roll by suppressing the surface roll, but the nonlinear noise will remain. After the removal, empty channels with zero amplitude value are likely to appear in the CMP domain gather, which makes it difficult to remove the empty channels.

[0007] Abnormal noise attenuation technology (AAA) can only suppress abnormal noise, environmental noise, and wild values ​​with large amplitudes in the artillery field, and it is difficult to suppress surface waves.

[0008] In addition, Wang (2013) proposed a joint surface roll suppression method based on time transform and shearlet transform. Ishiyama et al. (2014) proposed a method to estimate and suppress surface waves by using characteristics such as dispersion.

[0009] However, these methods cannot be applied to the actual production of massive data processing. In addition, with the deepening of ultra-deep carbonate fracture reservoir exploration, the protection of low-frequency information is becoming more and more important. It is very important to develop a technology that effectively suppresses surface wave noise and protects low-frequency information. Summary of the invention

[0010] The present application provides a surface roll suppression method, device and electronic device based on energy replacement to solve the above-mentioned technical problem that the prior art suppresses surface rolls poorly in protecting low-frequency information.

[0011] According to one aspect of the present application, an embodiment provides a surface roll suppression method based on energy replacement, comprising:

[0012] Analyze single shot scans containing surface wave seismic data to determine the bandwidth of surface wave interference;

[0013] Based on the frequency bandwidth of the surface roll interference, the original data is separated to obtain data within the frequency bandwidth of the surface roll and data outside the frequency bandwidth of the surface roll;

[0014] Based on the energy difference between the surface roll and the effective signal, a multi-time window gain analysis is performed on the data within the surface roll bandwidth range and the data outside the surface roll bandwidth range, and a gain operator of each seismic trace in each time window is calculated;

[0015] Performing gain operator substitution on the relatively larger gain operator in the effective wave outside the surface wave bandwidth range with the data within the surface wave bandwidth range; and using the larger gain operator after substitution to perform inverse gain calculation on the single shot signal within the surface wave bandwidth range to achieve surface wave suppression;

[0016] The data outside the surface roll frequency bandwidth separated from the original data are added to the data within the surface roll frequency bandwidth after energy replacement to obtain single-shot seismic data after surface roll suppression.

[0017] In one embodiment, the scanning analysis is set to double frequency scanning and spectrum analysis.

[0018] In one embodiment, the multi-time window gain analysis includes:

[0019] A multi-time window root mean square amplitude gain analysis is performed on the data within the surface roll frequency bandwidth and the data outside the surface roll frequency bandwidth channel by channel along the time axis.

[0020] In one embodiment, performing gain operator replacement of a relatively large gain operator in a valid wave outside the surface roll frequency bandwidth with data within the surface roll frequency bandwidth includes:

[0021] By comparing the gain operators, the relatively large gain operators in each time window of the effective wave outside the surface roll frequency bandwidth are replaced with the gain operators of the data within the surface roll frequency bandwidth.

[0022] In one embodiment, in the inverse gain calculation of the single shot signal within the frequency bandwidth of the surface wave, the larger gain operator in each time window after the substitution is used.

[0023] In one embodiment, in the surface wave suppression, the strong energy of the surface wave is greatly suppressed, while the amplitude of the effective signal outside the surface wave bandwidth range remains unchanged.

[0024] According to one aspect of the present application, an embodiment provides a surface wave suppression device based on energy replacement, characterized in that it includes:

[0025] A determination module, used for analyzing a single shot scan containing surface wave seismic data to determine the frequency bandwidth range of surface wave interference;

[0026] A separation module, used for separating the original data based on the frequency bandwidth of the surface roll interference to obtain data within the surface roll frequency bandwidth and data outside the surface roll frequency bandwidth;

[0027] A first processing module is used to perform multi-time window gain analysis on the data within the surface roll frequency bandwidth and the data outside the surface roll frequency bandwidth based on the energy difference between the surface roll and the effective signal, and calculate the gain operator of each seismic trace in each time window;

[0028] A suppression module is used to replace the relatively large gain operator in the effective wave outside the surface wave bandwidth with the data within the surface wave bandwidth; and use the larger gain operator after replacement to perform inverse gain calculation on the single shot signal within the surface wave bandwidth to suppress the surface wave; and

[0029] The second processing module is used to add the data outside the surface roll frequency bandwidth separated from the original data and the data within the surface roll frequency bandwidth after energy replacement to obtain the single-shot seismic data after surface roll suppression.

[0030] In one embodiment, the first processing module includes:

[0031] The first gain analysis module is used to perform multi-time window root mean square amplitude gain analysis on the data within the surface roll frequency bandwidth channel by channel along the time axis;

[0032] A second gain analysis module is used to perform multi-time window root mean square amplitude gain analysis on the data outside the surface roll frequency bandwidth along the time axis one by one; and

[0033] The calculation module is used to calculate the gain operator of each seismic trace in each time window.

[0034] According to one aspect of the present application, an embodiment provides an electronic device, including:

[0035] Memory; and

[0036] processor;

[0037] Wherein, the memory is used to store one or more computer instructions; the one or more computer instructions are executed by the processor to implement the surface roll suppression method based on energy replacement as described in any one of the above.

[0038] According to one aspect of the present application, an embodiment provides a readable storage medium having computer instructions stored thereon; wherein, when the computer instructions are executed by a processor, the surface roll suppression method based on energy replacement as described in any one of the above items is implemented.

[0039] The technical solution of the above embodiment of the present application adopts the frequency division energy exchange surface wave suppression technology for the surface waves existing in the seismic data, which can effectively suppress the surface waves, and the low-frequency effective components of the signal are basically not affected, the amplitude preservation is good, the data signal-to-noise ratio is improved, and there is no effective reflection information in the residual. Improve the structural imaging effect of the target layer such as pinch-out bodies, beads, caves, small faults and small-scale cracks, weak reflections and micro-structural geological bodies, and provide more accurate and reliable imaging results data for geological research. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 is a flow chart of a surface roll suppression method based on energy replacement in an embodiment;

[0041] Figure 2 is a flowchart of a surface roll suppression method based on energy replacement in an embodiment;

[0042] Figure 3 is a schematic diagram of single-shot frequency-doubled scanning in an embodiment;

[0043] Figure 4 is a schematic diagram of an original single shot and spectrum analysis in one embodiment;

[0044] Figure 5 is a schematic diagram of single shot frequency separation in an embodiment;

[0045] Figure 6 is a schematic diagram of single shot RMS gain analysis in one embodiment;

[0046] Figure 7 is a schematic diagram of an energy displacement single shot in one embodiment;

[0047] Figure 8 is a schematic diagram of noise comparison before and after noise suppression in an embodiment;

[0048] Fig. 9 is a schematic diagram of comparative analysis of an original single shot and a regional surface wave suppression and energy replacement single shot in an embodiment;

[0049] Fig.10 is a schematic diagram of comparative analysis of single-shot spectrum of regional surface wave suppression and energy displacement in an embodiment;

[0050] Fig.11 is a schematic diagram of a comparison of superimposed cross sections before and after noise suppression in an embodiment;

[0051] Fig.12 is a schematic diagram of a cross section of an offset comparison of new and old materials in an embodiment;

[0052] Fig.13 The diagram is a schematic structural diagram of a surface wave suppression device based on energy replacement in an embodiment. DETAILED DESCRIPTION

[0053] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0054] In order to enable those skilled in the art to better understand the present application, the technical solutions in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only embodiments of a part of the present application, not all embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in the field without creative work should fall within the scope of protection of the present application.

[0055] It should be noted that the terms "first", "second" etc. in the present application are used to distinguish similar objects, and need not be used to describe a specific order or sequential order. It should be understood that the data used in this way can be interchanged in appropriate circumstances, so that the embodiments of the present application described here. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, the process, method, system, product or equipment comprising a series of steps or units need not be limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or equipment.

[0056] It should be understood that when an element (such as a layer, film, region, or substrate) is described as being "on" another element, the element may be directly on the other element, or there may be an intermediate element. Moreover, in this application, when an element is described as being "connected" to another element, the element may be "directly connected" to the other element, or "connected" to the other element through a third element.

[0057] Embodiment 1

[0058] Please refer to Figure 1 , Figure 2 An embodiment provides a surface roll suppression method based on energy replacement, comprising the following steps:

[0059] 101. Perform single shot scanning analysis on surface wave seismic data to determine the frequency bandwidth of surface wave interference, wherein the scanning analysis is configured as frequency multiple scanning and spectrum analysis.

[0060] 102. Based on the frequency bandwidth of the surface roll interference, the original data is separated to obtain data within the surface roll frequency bandwidth and data outside the surface roll frequency bandwidth.

[0061] 103. Based on the energy difference between the surface roll and the effective signal, a multi-time window gain analysis is performed on the data within the surface roll bandwidth and the data outside the surface roll bandwidth, and a gain operator of each seismic trace in each time window is calculated. The multi-time window gain analysis includes: performing a multi-time window root mean square amplitude gain analysis on the data within the surface roll bandwidth and the data outside the surface roll bandwidth along the time axis.

[0062] 104. Perform gain operator replacement on the data within the surface roll frequency bandwidth with the relatively large gain operator in the effective wave outside the surface roll frequency bandwidth. Specifically, compare the gain operators and perform gain operator replacement on the data within the surface roll frequency bandwidth with the relatively large gain operator in each time window in the effective wave outside the surface roll frequency bandwidth.

[0063] By using the larger gain operator after substitution, preferably, using the larger gain operator in each time window after substitution, the single shot signal within the surface wave bandwidth is subjected to inverse gain calculation to achieve surface wave suppression. In the implementation of surface wave suppression, the strong energy of the surface wave is greatly suppressed, while the amplitude of the effective signal outside the surface wave bandwidth remains unchanged.

[0064] 105. The data outside the surface roll frequency bandwidth separated from the original data are added to the data within the surface roll frequency bandwidth after energy replacement to obtain single-shot seismic data after surface roll suppression.

[0065] This application proposes a surface wave suppression method based on energy replacement of frequency domain constraints, which is aimed at the problem that the surface waves existing in the original data of oil and gas exploration seismic acquisition cause low seismic signal-to-noise ratio and affect the structural imaging of seismic profiles. In general, this method first performs frequency analysis and spectrum scanning on the original single shot, and separates the data into two parts of data within the surface wave bandwidth and outside the surface wave bandwidth according to the surface wave frequency range. Secondly, according to the difference in energy between the surface wave and the effective signal, multi-time window root mean square gain analysis is performed on the data inside and outside the surface wave bandwidth along the time axis, and the gain operator of the seismic channel in each time window is calculated, and then the gain operator is replaced. Finally, the separated surface wave and effective signal are superimposed to obtain the single shot data after surface wave suppression, so as to achieve lossless suppression of the effective signal and improve the signal-to-noise ratio.

[0066] Embodiment 2

[0067] Please refer to Figure 1 , Figure 2 An embodiment provides a surface roll suppression method based on energy replacement. The following will explain the technical solution of the present application with the ideas of the inventor, which will help those skilled in the art to be familiar with the ideas of the present application. Among them, some technical points are explained in detail, so that those skilled in the art can fully understand the technical solution of the present application.

[0068] The surface wave suppression method based on energy replacement includes the following:

[0069] Step 1: Perform frequency scanning and analysis on a single shot to determine the frequency bandwidth of the surface wave.

[0070] Step 2: Separate the surface wave signals within the bandwidth and the signals outside the bandwidth from the original single shot.

[0071] Step 3: According to the energy difference between the surface wave and the effective signal, perform multi-time window root mean square amplitude gain analysis on the data inside and outside the surface wave bandwidth along the time axis, and calculate the gain operator of each seismic trace in each time window.

[0072] Step 4, compare each gain factor, and replace the relatively large gain operator in each time window of the effective wave outside the surface wave bandwidth with the data within the surface wave bandwidth; use the larger gain operator in each time window after replacement to perform inverse gain calculation on the single shot signal within the surface wave bandwidth, the strong energy of the surface wave will be greatly suppressed, while the amplitude of the effective signal outside the surface wave bandwidth remains unchanged, thereby achieving surface wave suppression.

[0073] Step 5: Add the previously separated effective signal outside the surface wave bandwidth to the signal within the surface wave bandwidth after energy replacement to obtain the single-shot seismic data after surface wave suppression.

[0074] The surface wave suppression method based on energy replacement in the present application is a noise suppression method based on energy exchange in a frequency-constrained domain, which mainly suppresses noise based on the difference in amplitude energy between surface waves and effective waves. Among them, within the surface wave frequency range, by suppressing the difference between the surface wave amplitude energy and the effective wave energy, the surface wave suppression effect can be guaranteed without losing the effective low-frequency signal in the surface wave area, thereby achieving the purpose of suppressing surface waves to protect low-frequency information and improving the imaging of ultra-deep fractures and holes.

[0075] Embodiment 3

[0076] Please refer to Figure 3-Figure 12 An embodiment provides a surface roll suppression method based on energy replacement, which is explained below through a case to facilitate further understanding of its technical solution.

[0077] The diffraction waves of the main target layer in a certain western exploration area in China are developed. The strong surface wave energy is mixed with the effective wave, which interferes with the imaging of the diffraction wave crack. It is difficult to effectively remove it by conventional methods, which is easy to damage the effective low-frequency information and it is difficult to achieve satisfactory results. According to the method of this application, and in conjunction with the accompanying drawings, a detailed description is given as follows:

[0078] Step 101, perform frequency scanning and spectrum analysis on a single shot containing typical surface waves to determine the frequency bandwidth of the surface wave interference, such as Figure 3 , Figure 4 .

[0079] Step 102, the original data is separated according to the frequency bandwidth of the surface wave interference into data within the surface wave frequency bandwidth and data outside the range, such as Figure 5 .

[0080] Step 103, according to the energy difference between the surface roll and the effective signal, the surface roll and the effective signal data are distributed / separated along the time axis to perform multi-time window root mean square amplitude gain analysis, and the gain operator of each seismic trace in each time window is calculated, such as Figure 6 .

[0081] Step 104, compare the gain factors, and replace the relatively large gain operator in each time window of the effective wave outside the surface wave bandwidth with the data within the surface wave bandwidth; use the larger gain operator in each time window after replacement to perform inverse gain calculation on the single shot signal within the surface wave bandwidth, the strong energy of the surface wave will be greatly suppressed, while the amplitude of the effective signal outside the surface wave bandwidth remains unchanged, so as to suppress the surface wave. Figure 7 ;

[0082] Step 105, add the previously separated effective signal outside the surface wave bandwidth range to the signal within the surface wave bandwidth range after energy replacement to obtain the single shot seismic data after surface wave suppression, such as Figure 8 .

[0083] like Figure 9-12 Based on the energy replacement noise suppression technology, good results have been achieved. From the comparison of single shots, the energy replacement suppression of surface wave noise has been well suppressed, and the noise suppression is relatively clean, which is significantly better than the regional surface wave suppression effect. The surface wave spectrum has been well broadened. From the perspective of the migration imaging effect, the signal-to-noise ratio is significantly better than the old data, the fault fracture is clear, and the beading anomaly is obvious, achieving ideal results.

[0084] Embodiment 4

[0085] Please refer to Fig.13 , an embodiment provides a surface wave suppression device based on energy replacement, which takes the following form:

[0086] The determination module 10 is used to analyze the single shot scanning containing surface wave seismic data to determine the frequency bandwidth range of the surface wave interference.

[0087] The separation module 20 is used to separate the original data based on the frequency bandwidth of the surface roll interference to obtain data within the surface roll frequency bandwidth and data outside the surface roll frequency bandwidth.

[0088] The first processing module 30 is used to perform multi-time window gain analysis on the data within the surface roll bandwidth and the data outside the surface roll bandwidth based on the energy difference between the surface roll and the effective signal, and calculate the gain operator of each seismic trace in each time window.

[0089] The suppression module 40 is used to replace the relatively larger gain operator in the effective wave outside the surface wave bandwidth range with the data within the surface wave bandwidth range; and use the larger gain operator after replacement to perform inverse gain calculation on the single shot signal within the surface wave bandwidth range to achieve surface wave suppression.

[0090] The second processing module 50 is used to add the data outside the surface roll frequency bandwidth separated from the original data to the data within the surface roll frequency bandwidth after energy replacement to obtain single-shot seismic data after surface roll suppression.

[0091] In one embodiment, the first processing module 30 includes: a first gain analysis module, a second gain analysis module and a calculation module. The first gain analysis module is used to perform multi-time window root mean square amplitude gain analysis on each channel of the data within the surface roll bandwidth along the time axis; the second gain analysis module is used to perform multi-time window root mean square amplitude gain analysis on each channel of the data outside the surface roll bandwidth along the time axis; and the calculation module is used to calculate the gain operator of each seismic channel in each time window.

[0092] It should be noted that the above-mentioned surface roll suppression device based on energy replacement is used to implement the surface roll suppression method based on energy replacement in the above-mentioned embodiments, and each module in the device corresponds to each step in the method.

[0093] Embodiment 5

[0094] Based on the same inventive concept, an embodiment of the present application provides an electronic device, comprising: a memory and a processor; wherein the memory is used to store one or more computer instructions; the one or more computer instructions are executed by the processor to implement the surface roll suppression method based on energy replacement as described in any one of the above embodiments.

[0095] Embodiment 6

[0096] Based on the same inventive concept, an embodiment of the present application provides a readable storage medium, on which computer instructions are stored; wherein, when the computer instructions are executed by a processor, the surface roll suppression method based on energy replacement described in any one of the above embodiments is implemented.

[0097] Among them, the above one or more computer instructions can be formed into a program.

[0098] The above program can be run in the processor, or it can be stored in the memory (or computer-readable medium), which includes permanent and non-permanent, removable and non-removable media. Information storage can be achieved by any method or technology. Information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, read-only compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, tape disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined in this article, computer-readable media does not include temporary computer-readable media (transitory media), such as modulated data signals and carrier waves.

[0099] These computer programs can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, so that the instructions executed on the computer or other programmable device provide instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps of the functions specified in one or more blocks can be implemented by different modules corresponding to different steps.

[0100] The above description is only the preferred embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A surface roll suppression method based on energy replacement, characterized in that: include: Analyze single shot scans containing surface wave seismic data to determine the bandwidth of surface wave interference; Based on the frequency bandwidth of the surface roll interference, the original data is separated to obtain data within the frequency bandwidth of the surface roll and data outside the frequency bandwidth of the surface roll; Based on the energy difference between the surface roll and the effective signal, a multi-time window gain analysis is performed on the data within the surface roll bandwidth range and the data outside the surface roll bandwidth range, and a gain operator of each seismic trace in each time window is calculated; Performing gain operator substitution on the relatively larger gain operator in the effective wave outside the surface wave bandwidth range with the data within the surface wave bandwidth range; and using the larger gain operator after substitution to perform inverse gain calculation on the single shot signal within the surface wave bandwidth range to achieve surface wave suppression; The data outside the surface roll frequency bandwidth separated from the original data are added to the data within the surface roll frequency bandwidth after energy replacement to obtain single-shot seismic data after surface roll suppression.

2. The surface roll suppression method based on energy replacement according to claim 1, characterized in that: The scanning analysis is set as frequency doubling scanning and spectrum analysis.

3. The surface wave suppression method based on energy replacement according to claim 1, characterized in that: The multi-time window gain analysis includes: A multi-time window root mean square amplitude gain analysis is performed on the data within the surface roll frequency bandwidth and the data outside the surface roll frequency bandwidth channel by channel along the time axis.

4. The surface wave suppression method based on energy replacement according to claim 1, characterized in that: The step of performing gain operator replacement of a relatively large gain operator in a valid wave outside the surface wave frequency bandwidth with data within the surface wave frequency bandwidth comprises: By comparing the gain operators, the relatively large gain operators in each time window of the effective wave outside the surface roll frequency bandwidth are replaced with the gain operators of the data within the surface roll frequency bandwidth.

5. The surface roll suppression method based on energy replacement according to claim 1, characterized in that: In the inverse gain calculation of the single shot signal within the frequency bandwidth of the surface wave, the larger gain operator in each time window after the substitution is used.

6. The surface wave suppression method based on energy replacement according to claim 1, characterized in that: In the implementation of suppressing the surface wave, the strong energy of the surface wave is greatly suppressed, while the amplitude of the effective signal outside the frequency bandwidth of the surface wave remains unchanged.

7. A surface wave suppression device based on energy replacement, characterized in that: include: A determination module, used for analyzing a single shot scan containing surface wave seismic data to determine the frequency bandwidth range of surface wave interference; A separation module, used for separating the original data based on the frequency bandwidth of the surface roll interference to obtain data within the surface roll frequency bandwidth and data outside the surface roll frequency bandwidth; A first processing module is used to perform multi-time window gain analysis on the data within the surface roll frequency bandwidth and the data outside the surface roll frequency bandwidth based on the energy difference between the surface roll and the effective signal, and calculate the gain operator of each seismic trace in each time window; A suppression module is used to replace the relatively large gain operator in the effective wave outside the surface wave bandwidth with the data within the surface wave bandwidth; and use the larger gain operator after replacement to perform inverse gain calculation on the single shot signal within the surface wave bandwidth to suppress the surface wave; and The second processing module is used to add the data outside the surface roll frequency bandwidth separated from the original data and the data within the surface roll frequency bandwidth after energy replacement to obtain the single-shot seismic data after surface roll suppression.

8. The surface wave suppression device based on energy replacement according to claim 7, characterized in that: The first processing module comprises: The first gain analysis module is used to perform multi-time window root mean square amplitude gain analysis on the data within the surface roll frequency bandwidth channel by channel along the time axis; A second gain analysis module is used to perform multi-time window root mean square amplitude gain analysis on the data outside the surface roll frequency bandwidth along the time axis one by one; and The calculation module is used to calculate the gain operator of each seismic trace in each time window.

9. An electronic device, characterized in that: include: Memory; and processor; Wherein, the memory is used to store one or more computer instructions; the one or more computer instructions are executed by the processor to implement the surface roll suppression method based on energy replacement as described in any one of claims 1 to 6.

10. A readable storage medium, characterized in that: The readable storage medium stores computer instructions; wherein, when the computer instructions are executed by a processor, the surface roll suppression method based on energy replacement as described in any one of claims 1 to 6 is implemented.

Images