Interference noise processing method and device, electronic equipment and storage medium

Abstract

The embodiment of the invention provides an interference noise processing method and device, electronic equipment and a storage medium. The method comprises the following steps: carrying out data field conversion processing on initial seismic data in a first data field to obtain intermediate seismic data in a second data field, determining first interference noise data in the second data field according to a preset interference filtering scanning mode based on the intermediate seismic data, and performing data field inverse conversion processing on the first interference noise data to obtain second interference noise data in the first data field, and performing interference noise processing on the initial seismic data based on the second interference noise data to obtain processed effective seismic data. According to the method, effective suppression of external source interference noise is realized, and the amplitude preservation fidelity is better; meanwhile, low-frequency effective signals can be protected, the signal-to-noise ratio of seismic data is increased, the method has adaptability, and compared with a traditional exogenous interference suppression algorithm, the seismic data with huge data volume can be efficiently processed, and the suppression effect is effectively improved.

Description

Technical Field

[0001] The present invention relates to the technical field of geophysical exploration seismic data processing, and in particular to an interference noise processing method, device, electronic equipment and storage medium. Background Art

[0002] In the acquisition of offshore seismic data, there are often various noises, such as surge interference, direct waves, oblique interference, multiple waves, and external interference. With the increase of offshore facilities and the development of fisheries, the sources of interference signals for offshore data acquisition are also increasing. Therefore, for seismic profiles with strong external interference energy for offshore data, special external interference suppression work needs to be carried out.

[0003] During the process of offshore seismic data acquisition, when an external source passes by, the noise caused by its movement will appear in different ways in the seismic data records, but its characteristics are generally similar to those of normal reflection signals. From the spectrum, the frequency range overlaps with the effective signal and conforms to the law of the seismic time-distance curve. Therefore, general noise suppression methods cannot identify and suppress external interference.

[0004] Therefore, the current seismic data processing technology has the technical problems of poor suppression effect and low efficiency on external interference noise. Summary of the invention

[0005] In order to alleviate the technical problems of poor suppression effect and low efficiency of external interference noise in current seismic data processing technology, the embodiments of the present invention provide an interference noise processing method, device, electronic device and storage medium.

[0006] In a first aspect, an embodiment of the present invention provides a method for processing interference noise in seismic data, the method comprising:

[0007] Performing data domain conversion processing on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain;

[0008] Based on the intermediate seismic data, determining first interference noise data in the second data domain according to a preset interference filtering scanning method;

[0009] Performing data domain inverse conversion processing on the first interference noise data to obtain second interference noise data in the first data domain;

[0010] Based on the second interference noise data, interference noise processing is performed on the initial seismic data to obtain processed valid seismic data.

[0011] In some implementations, the first data domain includes a time domain or a frequency domain, and the second data domain includes a Toup mapping domain.

[0012] In some embodiments, performing data domain conversion processing on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain includes:

[0013] Determining seismic wave observation coordinates according to the initial seismic data;

[0014] The intermediate seismic data is generated according to the slopes and intercepts of the tangent lines corresponding to all points in the seismic wave observation coordinates.

[0015] In some embodiments, generating the intermediate seismic data according to the slopes and intercepts of the tangent lines corresponding to all points in the seismic wave observation coordinates includes:

[0016] Determine the linear dynamic correction time difference threshold according to the slope and intercept of the tangent line corresponding to all the points;

[0017] Only when the linear normalized time difference range corresponding to the seismic wave observation coordinates is less than or equal to the linear normalized time difference threshold, the intermediate seismic data is generated according to the slope and intercept of the tangent corresponding to all points in the seismic wave observation coordinates.

[0018] In some implementations, performing data domain inverse conversion processing on the first interference noise data to obtain second interference noise data in the first data domain includes:

[0019] Determining the position coordinates of the noise data in the first data domain according to the slope and intercept of the noise data in the first interference noise data;

[0020] The second interference noise data is generated according to the position coordinates of the noise data.

[0021] In some implementations, before determining the position coordinates of the noise data in the first data domain according to the slope and intercept of the noise data in the first interference noise data, the method further includes:

[0022] Acquire a channel number threshold of a common shot channel gather and an actual channel value of the first interference noise data;

[0023] Only when the actual channel value is greater than or equal to the channel value threshold, the position coordinates of the noise data in the first data domain are determined according to the slope and intercept of the noise data in the first interference noise data.

[0024] In some embodiments, the interference filtering scanning method includes a random splitting method and a median filtering method; the determining the first interference noise data in the second data domain based on the intermediate seismic data according to a preset interference filtering scanning method includes:

[0025] Determining the range and increment of common shot gathers based on the intermediate seismic data;

[0026] According to the random splitting method, the common shot gathers are split into a plurality of crossover subsets;

[0027] According to the median filtering method, determining the filter corresponding to each cross subset;

[0028] According to the filters corresponding to each cross subset, noise is suppressed on the intermediate seismic data to obtain noise suppressed seismic data;

[0029] The first interference noise data is obtained according to the intermediate seismic data and the noise-suppressed seismic data.

[0030] In a second aspect, an embodiment of the present invention provides a device for processing interference noise in seismic data, the device comprising:

[0031] The first module is used to perform data domain conversion processing on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain;

[0032] A second module is used to determine first interference noise data in the second data domain based on the intermediate seismic data according to a preset interference filtering scanning method;

[0033] A third module is used to perform data domain inverse conversion processing on the first interference noise data to obtain second interference noise data in the first data domain;

[0034] The fourth module is used to perform interference noise processing on the initial seismic data based on the second interference noise data to obtain processed valid seismic data.

[0035] In a third aspect, an embodiment of the present invention provides an electronic device, including a memory and a processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the method described in the first aspect is implemented.

[0036] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed by one or more processors, the method described in the first aspect is implemented.

[0037] Compared with the prior art, one or more embodiments of the present invention can at least bring the following beneficial effects:

[0038] An embodiment of the present invention provides an interference noise processing method, device, electronic device and storage medium; the method includes performing data domain conversion processing on initial seismic data in a first data domain to obtain intermediate seismic data in a second data domain, based on the intermediate seismic data, determining first interference noise data in the second data domain according to a preset interference filtering scanning method, performing data domain inverse conversion processing on the first interference noise data to obtain second interference noise data in the first data domain, and based on the second interference noise data, performing interference noise processing on the initial seismic data to obtain processed valid seismic data. In the method provided in the present scheme, initial seismic data in a first data domain (xt domain) are subjected to data domain conversion processing to obtain intermediate seismic data in a second data domain (τ-p domain), and then first interference noise data in the second data domain are determined according to a preset interference filtering scanning method, and then the first interference noise data are subjected to data domain inverse conversion processing to obtain second interference noise data in the first data domain, and finally the second interference noise data are removed from the initial seismic data to obtain effective seismic data, thereby achieving effective suppression of exogenous interference noise with better amplitude preservation and fidelity; at the same time, it can protect low-frequency effective signals, improve the signal-to-noise ratio of seismic data, and has adaptability. Compared with traditional exogenous interference suppression algorithms, it can efficiently process seismic data with a large amount of data and effectively improve the suppression effect. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings required for use in the embodiments are briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present invention and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without creative work.

[0040] Figure 1 A schematic diagram of a flow chart of an interference noise processing method provided in an embodiment of the present application;

[0041] Figure 2 A schematic diagram showing the τ-p transformation principle provided in an embodiment of the present application;

[0042] Figure 3 A schematic diagram showing the external interference on a gather provided in an embodiment of the present application;

[0043] Figure 4 A comparison diagram of applying the interference noise processing method provided in an embodiment of the present application to suppress external interference on the original single-shot record;

[0044] Figure 5 A schematic diagram of the spectrum comparison of different ranges of a single gun before and after the external interference suppression provided in the embodiment of the present application;

[0045] Figure 6 Another schematic diagram of the spectrum comparison of different ranges of a single gun before and after the external interference suppression provided in the embodiment of the present application;

[0046] Figure 7 A comparison diagram of superimposed cross sections before and after external interference suppression provided in an embodiment of the present application;

[0047] Figure 8 A schematic diagram of the structure of an interference noise processing device provided in an embodiment of the present application.

[0048] In the drawings, the same reference numerals are used for the same components and the drawings are not drawn to scale. DETAILED DESCRIPTION

[0049] The following will describe the implementation methods of the present application in detail with reference to the accompanying drawings and embodiments, so that the implementation process of how the present application applies technical means to solve technical problems and achieve corresponding technical effects can be fully understood and implemented accordingly. The embodiments of the present application and the various features in the embodiments can be combined with each other without conflict, and the technical solutions formed are all within the scope of protection of the present application.

[0050] Meanwhile, in the following description, many specific details are set forth for the purpose of explanation to provide a thorough understanding of the embodiments of the present invention. However, it is obvious to those skilled in the art that the present invention can be implemented without the specific details or the specific manner described herein.

[0051] Embodiment 1

[0052] Figure 1 A flow chart of the interference noise processing method provided in the embodiment of the present application is shown in FIG. Figure 1 , the interference noise processing method provided in this embodiment includes:

[0053] Step S110: Perform data domain conversion processing on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain.

[0054] In some implementations, the first data domain includes a time domain or a frequency domain, and the second data domain includes a Toup mapping domain. Specifically, the time domain or the frequency domain is an xt domain; and the Toup mapping domain is a τ-p domain.

[0055] In the present application, the linear Radon transformation method is performed on the initial seismic data through the linear Tp domain common p gather (common shot gather), and the data is converted from the xt domain to the Tp domain.

[0056] In some embodiments, the data domain conversion processing is performed on the initial seismic data in the first data domain to obtain the intermediate seismic data in the second data domain, including: determining the seismic wave observation coordinates based on the initial seismic data; and generating the intermediate seismic data based on the slope and intercept of the tangent corresponding to all points in the seismic wave observation coordinates.

[0057] Specifically, the seismic wave observation coordinates (x, t) are determined based on the initial seismic data; the ordinary seismic record xt domain is converted to the Tp domain through the linear Radon transform, that is, all points in the seismic wave observation coordinates (x, t) are represented by the slope and intercept of the tangent line, that is, the vertical wave slowness component T and the horizontal wave slowness component p of the seismic wave. The Tp positive transform can be defined as the integral along a line in the xt domain, that is,

[0058]

[0059] Where f(x, t) is the xt domain seismic data volume, x is the offset, t is the time, τ is the intercept time corresponding to the zero offset trace in the seismic data, and p is the ray parameter. After r(p, τ) conversion, the data f(x, t) is reconstructed.

[0060] The discrete form of the above transformation is:

[0061] r(p,τ)=∑f(x,t=τ+px)Δx

[0062] Figure 2 The schematic diagram showing the Tp transformation principle provided in the embodiment of the present application is as follows: Figure 2 As shown in the figure, according to the transformation principle, wave fields with straight line characteristics, such as direct waves and surface waves, will be displayed in the form of points after linear Radon transformation; wave fields with hyperbolic characteristics, such as reflected wave fields, will be displayed in the form of elliptical arcs after linear Radon transformation.

[0063] In some embodiments, the intermediate seismic data is generated according to the slopes and intercepts of the tangents corresponding to all points in the seismic wave observation coordinates, including: determining a linear dynamic correction time difference threshold according to the slopes and intercepts of the tangents corresponding to all points in the seismic wave observation coordinates; and generating the intermediate seismic data according to the slopes and intercepts of the tangents corresponding to all points in the seismic wave observation coordinates only when the linear dynamic correction time difference range corresponding to the seismic wave observation coordinates is less than or equal to the linear dynamic correction time difference threshold.

[0064] In this application, the linear dynamic correction time difference threshold can be X mzx / (Δx·f max ).

[0065] Specifically, when converting from the xt domain to the Tp domain, aliasing will occur. The aliasing occurs during the conversion from the xt domain to the Tp domain, and is affected by the maximum offset, offset increment, and maximum frequency of the input data. In the Tp domain, if the following formula is met, aliasing will occur, causing Radon transform distortion.

[0066] LMO>X mzx / (Δx·f max )

[0067] Where LMO is the linear dynamic adjustment time difference range, that is, the absolute value between the maximum time difference and the minimum time difference, X mzx is the maximum offset, Δx is the track spacing, and f max is the maximum frequency.

[0068] Since in the Tp domain, the formula LMO>X is satisfied mzx / (Δx·f max ) will produce false frequency linearity, causing Radon transform distortion. Therefore, only in the linear dynamic correction time difference range LMO≤X corresponding to the seismic wave observation coordinates mzx / (Δx·f max ), the intermediate seismic data are generated according to the slope and intercept of the tangent line corresponding to all points in the seismic wave observation coordinates.

[0069] Step S120: Based on the intermediate seismic data, determine the first interference noise data in the second data domain according to a preset interference filtering scanning method.

[0070] In some embodiments, the interference filtering scanning method includes a random splitting method and a median filtering method; the first interference noise data in the second data domain is determined based on the intermediate seismic data according to a preset interference filtering scanning method, including: based on the intermediate seismic data, determining the range and increment of the common shot track set; according to the random splitting method, splitting the common shot track set into multiple cross subsets; according to the median filtering method, determining the filter corresponding to each cross subset; according to the filter corresponding to each cross subset, performing noise suppression on the intermediate seismic data to obtain noise suppressed seismic data; according to the intermediate seismic data and the noise suppressed seismic data, obtaining the first interference noise data.

[0071] Specifically, after converting the original seismic data record from the xt domain to the Tp domain to obtain intermediate seismic data, based on the intermediate seismic data, interference filtering scanning is performed on continuous shots according to amplitude and frequency, and interference filtering and interference-free τ-p domain data are identified, and then the extracted interference noise (i.e., first interference noise data) is obtained by subtracting the interference-free τ-p domain data from the input τ-p domain (i.e., second data domain) data.

[0072] It can be understood that after Radon transformation, a common p gather (common shot gather) is generated in the Tp domain. max The time difference t min and t max The range of the common p channel sets is determined; the increment of adjacent p channel sets is determined by half of the Nyquist frequency; then the common p channel sets in the Tp domain are randomly split into multiple cross subsets, and the channel sets in each cross subset are combined with median filtering to calculate the filter used to suppress external interference, thereby performing noise suppression.

[0073] Step S130: performing data domain inverse conversion processing on the first interference noise data to obtain second interference noise data in the first data domain.

[0074] In the present application, the data domain inverse conversion process converts the Tp domain back to the xt domain.

[0075] In this embodiment, following the previous step, after the first interference noise data is determined, the first interference noise data is subjected to data domain inverse conversion processing to convert it back to the xt domain, thereby obtaining second interference noise data in the first data domain.

[0076] In some embodiments, the data domain inverse conversion processing is performed on the first interference noise data to obtain the second interference noise data in the first data domain, including: determining the position coordinates of the noise data in the first data domain based on the slope and intercept of the noise data in the first interference noise data; and generating the second interference noise data based on the position coordinates of the noise data.

[0077] Specifically, the inverse transformation from the Tp domain to the xt domain can be written as:

[0078]

[0079] Among them, f(x, t) is the xt-domain seismic data volume, x is the offset, t is the time, τ is the intercept time corresponding to the zero-offset trace in the seismic data, and p is the ray parameter.

[0080] The discrete form of the above transformation is:

[0081] f(x,t)=∑r(p,τ=t-px)Δp

[0082] In some embodiments, before determining the position coordinates of the noise data in the first data domain according to the slope and intercept of the noise data in the first interference noise data, the method further includes: obtaining a channel number threshold of a common shot channel set and an actual channel value of the first interference noise data; and determining the position coordinates of the noise data in the first data domain according to the slope and intercept of the noise data in the first interference noise data only when the actual channel value is greater than or equal to the channel number threshold.

[0083] Specifically, affected by the common p-channel gathers (i.e., common shot channel gathers) in the Tp domain, aliasing will also occur during the conversion from the Tp domain to the xt domain; the calculation formula for the number of common shot channel gathers is as follows:

[0084] N p =2(t max -t min )·f max +1

[0085] Among them, N p is the total number of p gathers, t min and t max are the maximum and minimum dynamic calibration time difference, respectively, max is the maximum frequency.

[0086] If the number of common p-gather traces is too small, the near and far traces of the seismic wave time-distance curve are truncated, and the tilt superposition range is limited, which will cause the linear Radon transform to be distorted; therefore, before determining the position coordinates of the noise data in the first data domain, the trace number threshold of the common shot traces and the actual trace value of the first interference noise data are obtained; only when the actual trace value of the first interference noise data is greater than or equal to the trace number threshold, the position coordinates of the noise data in the first data domain are determined according to the slope and intercept of the noise data in the first interference noise data.

[0087] Step S140: Based on the second interference noise data, the initial seismic data is subjected to interference noise processing to obtain processed valid seismic data.

[0088] In this embodiment, following the previous step, after obtaining the second interference noise data, it is subtracted from the initial seismic data to obtain the processed effective seismic data. The method provided by this application can suppress noise well and has good amplitude preservation and fidelity; it can lay a good foundation for subsequent amplitude preservation processing and AVO pre-stack inversion.

[0089] like Figure 3As shown in the figure, it is a schematic diagram showing the external interference on the gather. Since the external interference is strongly correlated with the source information, it is impossible to separate and suppress it in the time domain or frequency domain. When the external interference is strong, if certain measures are not taken to suppress it, the imaging quality will be greatly reduced, which will also have a great impact on the subsequent seismic data processing and fine structural characterization. Figure 4 As shown in FIG. 1 , it is a comparison diagram of suppressing external interference of the original single shot record by applying the interference noise processing method provided in the embodiment of the present application. It can be seen that after applying the method provided in the present application to process the interference noise, the external interference is well suppressed, and the effective information of the reflected wave is almost not affected at all, and the amplitude and fidelity are well preserved. Figure 5 and Figure 6 As shown, it is a comparison of the spectrum of different ranges of a single gun before and after the external interference suppression. It can be seen that after the external interference noise is processed by the method provided in this application, the spectrum shape does not change much and there is basically no loss. Figure 7 The comparison diagram of the superimposed cross-section before and after the external interference suppression provided in the embodiment of the present application is shown in FIG. Figure 7 It can be seen that the external interference has been well suppressed and the effective information has been almost not affected, which lays a good foundation for subsequent amplitude preservation processing and AVO pre-stack inversion.

[0090] In summary, this embodiment performs data domain conversion processing on the initial seismic data in the first data domain (xt domain) to obtain intermediate seismic data in the second data domain (Tp domain), and then determines the first interference noise data in the second data domain according to a preset interference filtering scanning method, and then performs data domain inverse conversion processing on the first interference noise data to obtain the second interference noise data in the first data domain, and finally removes the second interference noise data from the initial seismic data to obtain effective seismic data, thereby achieving effective suppression of exogenous interference noise with better amplitude preservation and fidelity; at the same time, it can protect low-frequency effective signals, improve the signal-to-noise ratio of seismic data, and has adaptability. Compared with traditional exogenous interference suppression algorithms, it can efficiently process seismic data with a large amount of data and effectively improve the suppression effect.

[0091] Embodiment 2

[0092] Figure 8 For a schematic diagram of the structure of the interference noise processing device provided in the embodiment of the present application, please refer to Figure 8 , the interference noise processing device provided in this embodiment includes:

[0093] The first module 810 is used to perform data domain conversion processing on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain;

[0094] The second module 820 is used to determine the first interference noise data in the second data domain based on the intermediate seismic data according to a preset interference filtering scanning method;

[0095] The third module 830 is configured to perform data domain inverse conversion processing on the first interference noise data to obtain second interference noise data in the first data domain;

[0096] The fourth module 840 is used to perform interference noise processing on the initial seismic data based on the second interference noise data to obtain processed valid seismic data.

[0097] In some embodiments, the first module 810 is further used to determine the seismic wave observation coordinates based on the initial seismic data; and generate the intermediate seismic data based on the slope and intercept of the tangent corresponding to all points in the seismic wave observation coordinates.

[0098] In some embodiments, the first module 810 is also used to determine a linear dynamic correction time difference threshold based on the slopes and intercepts of the tangents corresponding to all the points; and generate the intermediate seismic data based on the slopes and intercepts of the tangents corresponding to all the points in the seismic wave observation coordinates only when the linear dynamic correction time difference range corresponding to the seismic wave observation coordinates is less than or equal to the linear dynamic correction time difference threshold.

[0099] In some embodiments, the second module 820 is further used to determine the range and increment of the common shot gather based on the intermediate seismic data; split the common shot gather into multiple cross subsets according to the random splitting method; determine the filter corresponding to each cross subset according to the median filtering method; suppress the noise of the intermediate seismic data according to the filter corresponding to each cross subset to obtain the noise-suppressed seismic data; obtain the first interference noise data according to the intermediate seismic data and the noise-suppressed seismic data.

[0100] In some embodiments, the third module 830 is further used to determine the position coordinates of the noise data in the first data domain based on the slope and intercept of the noise data in the first interference noise data; and generate the second interference noise data based on the position coordinates of the noise data.

[0101] The specific implementation examples based on the above-mentioned modules for executing the method have been described in detail in the first embodiment and will not be repeated here.

[0102] Embodiment 3

[0103] This embodiment provides an electronic device, including a memory and a processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the interference noise processing method described in Embodiment 1 is implemented. It can be understood that the electronic device may also include an input / output (I / O) interface and a communication component.

[0104] The processor is used to execute all or part of the steps in the interference noise processing method in Embodiment 1. The memory is used to store various types of data, which may include instructions of any application or method in the terminal device, and data related to the application.

[0105] The processor can be an application specific integrated circuit (ASIC), a digital signal processor (DSP), a digital signal processing device (DSPD), a programmable logic device (PLD), a field programmable gate array (FPGA), a controller, a microcontroller, a microprocessor or other electronic components, and is used to execute the interference noise processing method in the above-mentioned embodiment 1.

[0106] The memory can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic storage, flash memory, magnetic disk or optical disk.

[0107] Embodiment 4

[0108] This embodiment also provides a computer-readable storage medium, such as a flash memory, a hard disk, a multimedia card, a card-type memory (for example, an SD or DX memory, etc.), a random access memory (RAM), a static random access memory (SRAM), a read-only memory (ROM), an electrically erasable programmable read-only memory (EEPROM), a programmable read-only memory (PROM), a magnetic memory, a disk, an optical disk, a server, an App application store, etc., wherein a computer program is stored on the storage medium, and when the computer program is executed by a processor, the following method steps can be implemented:

[0109] Performing data domain conversion processing on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain;

[0110] Based on the intermediate seismic data, determining first interference noise data in the second data domain according to a preset interference filtering scanning method;

[0111] Performing data domain inverse conversion processing on the first interference noise data to obtain second interference noise data in the first data domain;

[0112] Based on the second interference noise data, interference noise processing is performed on the initial seismic data to obtain processed valid seismic data.

[0113] The specific implementation process of the above method steps can be found in Example 1, and this embodiment will not be repeated here.

[0114] In summary, the present application provides an interference noise processing method, device, electronic device and storage medium; the method includes performing data domain conversion processing on initial seismic data in a first data domain to obtain intermediate seismic data in a second data domain, based on the intermediate seismic data, determining first interference noise data in the second data domain according to a preset interference filtering scanning method, performing data domain inverse conversion processing on the first interference noise data to obtain second interference noise data in the first data domain, and based on the second interference noise data, performing interference noise processing on the initial seismic data to obtain processed valid seismic data. In the method provided in the present scheme, initial seismic data in a first data domain (xt domain) are subjected to data domain conversion processing to obtain intermediate seismic data in a second data domain (τ-p domain), and then first interference noise data in the second data domain are determined according to a preset interference filtering scanning method, and then the first interference noise data are subjected to data domain inverse conversion processing to obtain second interference noise data in the first data domain, and finally the second interference noise data are removed from the initial seismic data to obtain effective seismic data, thereby achieving effective suppression of exogenous interference noise with better amplitude preservation and fidelity; at the same time, it can protect low-frequency effective signals, improve the signal-to-noise ratio of seismic data, and has adaptability. Compared with traditional exogenous interference suppression algorithms, it can efficiently process seismic data with a large amount of data and effectively improve the suppression effect.

[0115] It should be noted that, in this article, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "comprises a ..." does not exclude the existence of other identical elements in the process, method, article or device including the element.

[0116] Although the implementation methods disclosed in this application are as above, the contents described are only implementation methods adopted for facilitating the understanding of this application, and are not intended to limit this application. Any technician in the technical field to which this application belongs can make any modifications and changes in the form and details of implementation without departing from the spirit and scope disclosed in this application, but the scope of patent protection of this application shall still be based on the scope defined in the attached claims.

Claims

1. A method for processing interference noise in seismic data, characterized in that: include: Performing data domain conversion processing on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain; Based on the intermediate seismic data, determining first interference noise data in the second data domain according to a preset interference filtering scanning method; Performing data domain inverse conversion processing on the first interference noise data to obtain second interference noise data in the first data domain; Based on the second interference noise data, interference noise processing is performed on the initial seismic data to obtain processed valid seismic data.

2. The interference noise processing method according to claim 1, characterized in that: The first data domain includes a time domain or a frequency domain, and the second data domain includes a Toup mapping domain.

3. The interference noise processing method according to claim 2, characterized in that: The performing data domain conversion processing on the initial seismic data in the first data domain to obtain the intermediate seismic data in the second data domain includes: Determining seismic wave observation coordinates according to the initial seismic data; The intermediate seismic data is generated according to the slopes and intercepts of the tangent lines corresponding to all points in the seismic wave observation coordinates.

4. The interference noise processing method according to claim 3, characterized in that: Generating the intermediate seismic data according to the slopes and intercepts of the tangent lines corresponding to all points in the seismic wave observation coordinates comprises: Determine the linear dynamic correction time difference threshold according to the slope and intercept of the tangent line corresponding to all the points; Only when the linear normalized time difference range corresponding to the seismic wave observation coordinates is less than or equal to the linear normalized time difference threshold, the intermediate seismic data is generated according to the slope and intercept of the tangent corresponding to all points in the seismic wave observation coordinates.

5. The interference noise processing method according to claim 2, characterized in that: The performing data domain inverse conversion processing on the first interference noise data to obtain second interference noise data in the first data domain includes: Determining the position coordinates of the noise data in the first data domain according to the slope and intercept of the noise data in the first interference noise data; The second interference noise data is generated according to the position coordinates of the noise data.

6. The interference noise processing method according to claim 5, characterized in that: Before determining the position coordinates of the noise data in the first data domain according to the slope and intercept of the noise data in the first interference noise data, the method further includes: Acquire a channel number threshold of a common shot channel gather and an actual channel value of the first interference noise data; Only when the actual channel value is greater than or equal to the channel value threshold, the position coordinates of the noise data in the first data domain are determined according to the slope and intercept of the noise data in the first interference noise data.

7. The interference noise processing method according to any one of claims 1 to 6, characterized in that: The interference filtering scanning method includes a random splitting method and a median filtering method; the first interference noise data in the second data domain is determined based on the intermediate seismic data according to a preset interference filtering scanning method, including: Determining the range and increment of common shot gathers based on the intermediate seismic data; According to the random splitting method, the common shot gathers are split into a plurality of crossover subsets; According to the median filtering method, determining the filter corresponding to each cross subset; According to the filters corresponding to each cross subset, noise is suppressed on the intermediate seismic data to obtain noise suppressed seismic data; The first interference noise data is obtained according to the intermediate seismic data and the noise-suppressed seismic data.

8. A device for processing interference noise in seismic data, characterized in that: include: The first module is used to perform data domain conversion processing on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain; A second module is used to determine first interference noise data in the second data domain based on the intermediate seismic data according to a preset interference filtering scanning method; A third module is used to perform data domain inverse conversion processing on the first interference noise data to obtain second interference noise data in the first data domain; The fourth module is used to perform interference noise processing on the initial seismic data based on the second interference noise data to obtain processed valid seismic data.

9. An electronic device, characterized in that: The method comprises a memory and a processor; a computer program is stored in the memory, and when the computer program is executed by the processor, the method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by one or more processors, the method according to any one of claims 1 to 7 is implemented.

Images