Interference noise processing method and device, electronic equipment and storage medium
By converting seismic data from the time or frequency domain to the Taupe mapping domain, and using Radon transform and interference filtering to identify and suppress external interference noise, the problem of poor suppression of external interference noise in marine seismic data is solved, achieving efficient noise removal and signal protection.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- CHINA PETROLEUM & CHEMICAL CORP
- Filing Date
- 2023-09-12
- Publication Date
- 2026-05-12
AI Technical Summary
Existing technologies are ineffective and inefficient at suppressing external interference noise in marine seismic data processing, and cannot effectively identify and remove external interference signals.
By transforming seismic data from the time or frequency domain to the Taupe mapping domain, external interference noise is identified and suppressed using linear Radon transform and interference filtering scanning, and then the data is inversely transformed back to the initial data domain for removal.
It effectively suppresses external interference noise, maintains good amplitude and fidelity, improves the signal-to-noise ratio of seismic data, and can efficiently process large-scale seismic data, thus enhancing the suppression effect.
Smart Images
Figure CN119620184B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of geophysical exploration seismic data processing technology, specifically to a method, apparatus, electronic device, and storage medium for processing interference noise. Background Technology
[0002] In marine seismic data acquisition, various types of noise are often encountered, such as swell interference, direct waves, oblique waves, multiple waves, and external interference. With the increase in marine facilities and the development of fisheries, the sources of interference signals in marine data acquisition are also increasing. Therefore, for seismic profiles with strong external interference energy in marine data, specialized external interference suppression work is required.
[0003] During the acquisition of marine seismic data, when an external source passes by, the noise generated during its movement will manifest in different ways in the seismic data record. However, its characteristics are generally similar to those of normal reflected signals. From the spectrum, the frequency range overlaps with the effective signal and conforms to the seismic time-distance curve. Therefore, general noise suppression methods cannot identify and suppress external interference.
[0004] Therefore, current seismic data processing technology suffers from poor suppression of external interference noise and low efficiency. Summary of the Invention
[0005] To alleviate the technical problems of poor suppression effect and low efficiency of current seismic data processing technology for external interference noise, embodiments of the present invention provide an interference noise processing method, device, electronic device and storage medium.
[0006] In a first aspect, embodiments of the present invention provide a method for processing interference noise in seismic data, the method comprising:
[0007] The initial seismic data in the first data domain is transformed to obtain the intermediate seismic data in the second data domain.
[0008] Based on the intermediate seismic data, the first interference noise data in the second data domain is determined according to the preset interference filtering scanning method;
[0009] The first interference noise data is subjected to data domain inverse transformation to obtain the second interference noise data in the first data domain.
[0010] Based on the second interference noise data, the initial seismic data is processed to obtain processed effective 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 Taupe mapping domain.
[0012] In some implementations, the step of performing data domain transformation on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain includes:
[0013] Based on the initial seismic data, determine the observation coordinates of the seismic waves;
[0014] The intermediate seismic data is generated based on the slope and intercept of the tangents corresponding to all points in the seismic wave observation coordinates.
[0015] In some implementations, generating the intermediate seismic data based on the slope and intercept of the tangents corresponding to all points in the seismic wave observation coordinates includes:
[0016] The linear dynamic calibration time difference threshold is determined based on the slope and intercept of the tangents corresponding to all the points.
[0017] The intermediate seismic data are generated only when the linear dynamic time difference range corresponding to the seismic wave observation coordinates is less than or equal to the linear dynamic time difference threshold, based on the slope and intercept of the tangent corresponding to all points in the seismic wave observation coordinates.
[0018] In some embodiments, the step of performing a data domain inverse transformation on the first interference noise data to obtain second interference noise data within the first data domain includes:
[0019] Based on the slope and intercept of the noise data in the first interference noise data, determine the position coordinates of the noise data in the first data domain;
[0020] The second interference noise data is generated based on the location coordinates of the noise data.
[0021] In some embodiments, before determining the position coordinates of the noise data within the first data domain based on the slope and intercept of the noise data in the first interference noise data, the method further includes:
[0022] Obtain the threshold number of the common shot gather and the actual number of the first interference noise data;
[0023] Only when the actual channel number is greater than or equal to the channel number threshold, the position coordinates of the noise data in the first data domain are determined based on the slope and intercept of the noise data in the first interference noise data.
[0024] In some implementations, the interference filtering scanning method includes a random splitting method and a median filtering method; the step of determining the first interference noise data in the second data domain based on the intermediate seismic data and according to the preset interference filtering scanning method includes:
[0025] Based on the intermediate seismic data, the scope and increment of the common shot gather were determined;
[0026] According to the random splitting method, the common artillery trajectory set is split into multiple intersecting subsets;
[0027] Based on the median filtering method, determine the filter corresponding to each cross subset;
[0028] Based on the filters corresponding to each cross subset, noise suppression is applied to the intermediate seismic data to obtain noise-suppressed seismic data.
[0029] The first interference noise data is obtained based on the intermediate seismic data and the noise-suppressed seismic data.
[0030] In a second aspect, embodiments of the present invention provide an interference noise processing device for seismic data, the device comprising:
[0031] The first module is used to perform data domain transformation processing on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain.
[0032] The second module is used to determine the first interference noise data in the second data domain based on the intermediate seismic data and according to a preset interference filtering scanning method.
[0033] The third module is used to perform data domain inverse transformation on the first interference noise data to obtain the second interference noise data in the first data domain.
[0034] The fourth module is used to process the initial seismic data based on the second interference noise data to obtain processed effective seismic data.
[0035] Thirdly, embodiments of the present invention provide an electronic device, including a memory and a processor, wherein the memory stores a computer program, and the computer program, when executed by the processor, implements the method described in the first aspect.
[0036] Fourthly, embodiments of the present invention provide a computer-readable storage medium storing a computer program, which, when executed by one or more processors, implements the method described in the first aspect.
[0037] Compared with the prior art, one or more embodiments of the present invention can bring at least the following beneficial effects:
[0038] This invention provides an interference noise processing method, apparatus, 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; determining first interference noise data in the second data domain based on the intermediate seismic data 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 performing interference noise processing on the initial seismic data based on the second interference noise data to obtain processed effective seismic data. The method provided in this scheme transforms the initial seismic data in the first data domain (xt domain) to obtain intermediate seismic data in the second data domain (τ-p domain). Then, based on a preset interference filtering scanning method, the first interference noise data in the second data domain is determined. Next, the first interference noise data undergoes inverse data domain transformation to obtain the second interference noise data in the first data domain. Finally, the second interference noise data is removed from the initial seismic data to obtain the effective seismic data. This method effectively suppresses external interference noise with better amplitude and fidelity preservation. It also protects low-frequency effective signals, improves the signal-to-noise ratio of seismic data, and has adaptive properties. Compared with traditional external interference suppression algorithms, it can efficiently process large amounts of seismic data and effectively improve the suppression effect. Attached Figure Description
[0039] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of the present invention and should not be regarded as a limitation on the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.
[0040] Figure 1 This is a schematic flowchart of an interference noise processing method provided in an embodiment of this application;
[0041] Figure 2 This is a schematic diagram illustrating the principle of τ-p transformation provided in the embodiments of this application;
[0042] Figure 3 A schematic diagram illustrating the external interference on the gather provided in the embodiments of this application;
[0043] Figure 4 A comparison diagram showing the application of the interference noise processing method provided in this application embodiment to suppress external interference in the original single-shot record;
[0044] Figure 5 This is a schematic diagram showing a comparison of the spectrum of a single gun before and after suppression of external interference, provided in an embodiment of this application.
[0045] Figure 6 This is another schematic diagram showing the comparison of the spectrum of a single gun at different ranges before and after the suppression of external interference, as provided in the embodiments of this application.
[0046] Figure 7 A comparison diagram of superimposed cross-sections before and after external interference suppression provided in the embodiments of this application;
[0047] Figure 8 This is a schematic diagram of the interference noise processing device provided in the embodiments of this application.
[0048] In the accompanying drawings, the same parts are referred to by the same reference numerals, and the drawings are not drawn to scale. Detailed Implementation
[0049] The following detailed description of the embodiments of this application, in conjunction with the accompanying drawings, will provide a thorough understanding of how this application uses technical means to solve technical problems and achieve corresponding technical effects, enabling its implementation. The embodiments of this application and the various features within them can be combined with each other without conflict, and all resulting technical solutions are within the protection scope of this application.
[0050] Furthermore, numerous specific details are set forth in the following description for purposes of explanation, in order to provide a thorough understanding of the embodiments of the invention. However, it will be apparent to those skilled in the art that the invention may be practiced without the specific details herein or the particular methods described.
[0051] Example 1
[0052] Figure 1 Please refer to the flowchart of an interference noise processing method provided in the embodiments of this application. Figure 1 The interference noise processing method provided in this embodiment includes:
[0053] Step S110: Perform data domain transformation 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 Taupe-Phase domain. Specifically, the time domain or frequency domain is the xt domain; the Taupe-Phase domain is the τ-p domain.
[0055] In this application, the initial seismic data are transformed from the xt domain to the τ-p domain by using the linear τ-p domain common p gather (common shot gather).
[0056] In some implementations, the step of performing data domain transformation processing on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain includes: 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 tangents corresponding to all points in the seismic wave observation coordinates.
[0057] Specifically, based on the initial earthquake data, the observed coordinates (x, t) of the seismic waves are determined. A linear Radon transform is then used to convert the ordinary seismic record from the xt domain to the τ-p domain. This means representing all points in the observed coordinates (x, t) of the seismic waves using the slope and intercept of the tangent line, i.e., the vertical wave slowness component τ and the horizontal wave slowness component p. The τ-p positive transform can be defined as the integral along a line in the xt domain, i.e.
[0058]
[0059] in It is the seismic data volume in the XT domain. This is the offset distance. It is time. It is the intercept time corresponding to the zero offset trace in the seismic data. These are the ray parameters. After... After conversion, the data It was reconstructed.
[0060] The discrete form of the above transformation is:
[0061]
[0062] Figure 2 This is a schematic diagram illustrating the principle of τ-p transformation provided in the embodiments of this application, such as... Figure 2 As shown, according to the transformation principle, wave fields with linear characteristics, such as direct waves and surface waves, will be represented as points after the linear Radon transform; wave fields with hyperbolic characteristics, such as reflected wave fields, will be represented as elliptical arcs after the linear Radon transform.
[0063] In some implementations, generating the intermediate seismic data based on the slope and intercept of the tangents corresponding to all points in the seismic wave observation coordinates includes: determining a linear dynamic correction time difference threshold based on the slope and intercept of the tangents corresponding to all points; and generating the intermediate seismic data 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 calibration time difference threshold can be .
[0065] Specifically, spurious frequency distortion occurs during the transformation from the xt domain to the τ-p domain. This distortion arises during the transformation and is influenced by the maximum offset, the offset increment, and the maximum frequency of the input data. In the τ-p domain, spurious frequency distortion occurs if the following formula is satisfied, causing distortion in the Radon transform.
[0066]
[0067] in This represents the linear dynamic time difference range, i.e., the absolute value between the maximum and minimum time difference. It is the maximum offset distance. It is a distance, That is the maximum frequency.
[0068] Since in the τ-p domain, the formula is satisfied This can produce spurious frequency linearity, causing distortion in the Radon transform. Therefore, linear dynamic time difference correction is only applied within the range corresponding to the seismic wave observation coordinates. Intermediate seismic data are generated based on the slope and intercept of the tangents 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 the preset interference filtering scanning method.
[0070] In some implementations, the interference filtering scanning method includes a random splitting method and a median filtering method; determining the first interference noise data in the second data domain based on the intermediate seismic data and according to the preset interference filtering scanning method includes: determining the range and increment of the common shot set based on the intermediate seismic data; splitting the common shot set into multiple cross subsets according to the random splitting method; determining the filter corresponding to each cross subset according to the median filtering method; performing noise suppression on the intermediate seismic data according to the filter corresponding to each cross subset to obtain noise-suppressed seismic data; and obtaining the first interference noise data based on the intermediate seismic data and the noise-suppressed seismic data.
[0071] Specifically, after converting the original seismic data records from the xt domain to the τ-p domain to obtain intermediate seismic data, based on this intermediate seismic data, interference filtering is performed on the continuous shot according to the amplitude and frequency, and the interference-filtered and interference-free τ-p domain data are identified. Then, by subtracting the interference-free τ-p domain data from the input τ-p domain data (i.e., the second data domain), the extracted interference noise (i.e., the first interference noise data) is obtained.
[0072] It is understandable that the Radon transform produces a common p-gathering (common gun gather) in the τ-p domain, at the maximum offset. time difference and The range of the common p-channels is determined; the increment of adjacent p-channels is determined by half of the Nyquist frequency; then the common p-channels in the τ-p domain are randomly split into multiple cross subsets, and the channels in each cross subset are combined with median filtering to calculate a filter for suppressing external interference, thereby performing noise suppression.
[0073] Step S130: Perform data domain inverse transformation on the first interference noise data to obtain the second interference noise data in the first data domain.
[0074] In this application, the data domain inverse transformation process is to transform the τ-p domain back to the xt domain.
[0075] In this embodiment, following the previous step, after determining the first interference noise data, the first interference noise data is subjected to data domain inverse transformation to convert it back to the xt domain, thereby obtaining the second interference noise data in the first data domain.
[0076] In some implementations, the step of performing a data domain inverse transformation on the first interference noise data to obtain second interference noise data within the first data domain includes: determining the position coordinates of the noise data within 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 τ-p domain back to the xt domain can be written as:
[0078]
[0079] in, It is the seismic data volume in the XT domain. This is the offset distance. It is time. It is the intercept time corresponding to the zero offset trace in the seismic data. These are the ray parameters.
[0080] The discrete form of the above transformation is:
[0081]
[0082] In some implementations, before 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, the method further includes: obtaining a common shot gather trace number threshold and the actual trace value of the first interference noise data; and determining the position coordinates of the noise data in the first data domain only when the actual trace value is greater than or equal to the trace number threshold, based on the slope and intercept of the noise data in the first interference noise data.
[0083] Specifically, due to the influence of the common p-gathering (i.e., common gun-gathering) in the τ-p domain, spurious frequency phenomena will also occur during the transformation from the τ-p domain to the xt domain; the formula for calculating the number of common gun-gathering is as follows:
[0084]
[0085] in, For a total of p channels, and These are the maximum and minimum time differences for dynamic calibration, respectively. That is the maximum frequency.
[0086] If the number of traces in the common p-gather is too small, the near and far traces of the seismic wave travel curve will be truncated, and the tilt superposition range will be limited, which will cause distortion of the linear Radon transform. Therefore, before determining the position coordinates of the noise data in the first data domain, the trace number threshold of the common shot gather 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 interference noise data 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 effective 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 in this application can effectively suppress noise and has good amplitude and fidelity preservation; it can lay a good foundation for subsequent amplitude preservation processing and AVO pre-stack inversion.
[0089] like Figure 3 The diagram illustrates the effect of external interference on the gather. Because external interference is strongly correlated with source information, it cannot be separated and suppressed in the time or frequency domains. When external interference is strong, without suppression measures, it will significantly reduce image quality and greatly impact subsequent seismic data processing and detailed structural characterization. Figure 4The image shown is a comparison of the application of the interference noise processing method provided in this application to suppress external interference in the original single-shot record. It can be seen that after applying the method provided in this application to process the interference noise, the external interference is well suppressed, and the effective information of the reflected wave is almost unaffected, demonstrating good amplitude and fidelity preservation. Figure 5 and Figure 6 As shown, the spectrum of a single gun is compared before and after the suppression of external interference. It can be seen that after applying the method provided in this application to process external interference noise, the spectrum shape does not change much and there is basically no loss. Figure 7 This is a comparison diagram of the superimposed cross-sections before and after external interference suppression provided in the embodiments of this application. Figure 7 It can be seen that external interference was well suppressed and the effective information was almost unaffected, laying a good foundation for subsequent amplitude preservation processing and AVO pre-stack inversion.
[0090] In summary, this embodiment performs data domain transformation on the initial seismic data in the first data domain (xt domain) to obtain intermediate seismic data in the second data domain (τ-p domain). Then, it determines the first interference noise data in the second data domain according to a preset interference filtering scanning method. Next, it performs data domain inverse transformation on the first interference noise data to obtain the second interference noise data in the first data domain. Finally, it removes the second interference noise data from the initial seismic data to obtain the effective seismic data. This achieves effective suppression of external interference noise with better amplitude and fidelity preservation. At the same time, it can protect the effective low-frequency signal, improve the signal-to-noise ratio of seismic data, and has adaptiveness. Compared with traditional external interference suppression algorithms, it can efficiently process large amounts of seismic data and effectively improve the suppression effect.
[0091] Example 2
[0092] Figure 8 Please refer to the schematic diagram of the interference noise processing device provided in the embodiments of this application. Figure 8 The interference noise processing device provided in this embodiment includes:
[0093] The first module 810 is used to perform data domain transformation 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 and according to a preset interference filtering scanning method.
[0095] The third module 830 is used to perform data domain inverse conversion processing on the first interference noise data to obtain the 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 effective seismic data.
[0097] In some implementations, the first module 810 is further configured to determine the seismic wave observation coordinates based on the initial seismic data; and to generate the intermediate seismic data based on the slope and intercept of the tangents corresponding to all points in the seismic wave observation coordinates.
[0098] In some implementations, the first module 810 is further configured to determine a linear dynamic time difference threshold based on the slope and intercept of the tangents corresponding to all points; and to generate the intermediate seismic data only when the linear dynamic time difference range corresponding to the seismic wave observation coordinates is less than or equal to the linear dynamic time difference threshold, based on the slope and intercept of the tangents corresponding to all points in the seismic wave observation coordinates.
[0099] In some embodiments, the second module 820 is further configured to: determine the range and increment of the common shot set based on the intermediate seismic data; divide the common shot set into multiple cross subsets according to the random splitting method; determine the filter corresponding to each cross subset according to the median filtering method; perform noise suppression on the intermediate seismic data according to the filter corresponding to each cross subset to obtain noise-suppressed seismic data; and obtain the first interference noise data based on the intermediate seismic data and the noise-suppressed seismic data.
[0100] In some embodiments, the third module 830 is further configured 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] Specific implementation methods based on the above modules have been detailed in Implementation 1 and will not be repeated here.
[0102] Example 3
[0103] This embodiment provides an electronic device, including a memory and a processor. The memory stores a calculator program, which, when executed by the processor, implements the interference and noise processing method described in Embodiment 1. It is understood that the electronic device may further include an input / output (I / O) interface and communication components.
[0104] The processor is used to execute all or part of the steps in the interference noise processing method as described in Embodiment 1. The memory is used to store various types of data, which may include, for example, instructions for any application or method in the terminal device, as well as application-related data.
[0105] The processor may be implemented as an Application Specific Integrated Circuit (ASIC), Digital Signal Processor (DSP), Digital Signal Processing Device (DSPD), Programmable Logic Device (PLD), Field Programmable Gate Array (FPGA), controller, microcontroller, microprocessor, or other electronic components, and is used to execute the interference noise processing method in Embodiment 1 above.
[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] Example 4
[0108] This embodiment also provides a computer-readable storage medium, such as flash memory, hard disk, multimedia card, card-type memory (e.g., SD or DX memory), random access memory (RAM), static random access memory (SRAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), programmable read-only memory (PROM), magnetic memory, disk, optical disk, server, app store, etc., on which a computer program is stored. When the computer program is executed by a processor, it can implement the following method steps:
[0109] The initial seismic data in the first data domain is transformed to obtain the intermediate seismic data in the second data domain.
[0110] Based on the intermediate seismic data, the first interference noise data in the second data domain is determined according to the preset interference filtering scanning method;
[0111] The first interference noise data is subjected to data domain inverse transformation to obtain the second interference noise data in the first data domain.
[0112] Based on the second interference noise data, the initial seismic data is processed to obtain processed effective seismic data.
[0113] For a detailed description of the above method steps, please refer to Example 1. This example will not be repeated here.
[0114] In summary, this application provides an interference noise processing method, apparatus, 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; determining first interference noise data in the second data domain based on the intermediate seismic data 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 performing interference noise processing on the initial seismic data based on the second interference noise data to obtain processed effective seismic data. The method provided in this scheme transforms the initial seismic data in the first data domain (xt domain) to obtain intermediate seismic data in the second data domain (τ-p domain). Then, based on a preset interference filtering scanning method, the first interference noise data in the second data domain is determined. Next, the first interference noise data undergoes inverse data domain transformation to obtain the second interference noise data in the first data domain. Finally, the second interference noise data is removed from the initial seismic data to obtain the effective seismic data. This method effectively suppresses external interference noise with better amplitude and fidelity preservation. It also protects low-frequency effective signals, improves the signal-to-noise ratio of seismic data, and has adaptive properties. Compared with traditional external interference suppression algorithms, it can efficiently process large amounts of seismic data and effectively improve the suppression effect.
[0115] It should be noted that, in this document, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.
[0116] Although the embodiments disclosed in this application are as described above, the content is merely for the purpose of facilitating understanding of this application and is not intended to limit this application. Any person skilled in the art to which this application pertains may make any modifications and changes in the form and details of the implementation without departing from the spirit and scope disclosed in this application; however, the scope of patent protection of this application shall still be determined by the scope defined in the appended claims.
Claims
1. A method for processing interference noise in seismic data, characterized in that, include: The initial seismic data in the first data domain is transformed to obtain the intermediate seismic data in the second data domain. The initial seismic data is seismic data obtained from marine seismic data acquisition; Based on the intermediate seismic data, the first interference noise data in the second data domain is determined according to the preset interference filtering scanning method; The first interference noise data is subjected to data domain inverse transformation to obtain the second interference noise data in the first data domain. Based on the second interference noise data, the initial seismic data is subjected to interference noise processing to obtain processed effective seismic data; The interference filtering scanning method includes a random splitting method and a median filtering method; the step of determining the first interference noise data in the second data domain based on the intermediate seismic data and according to the preset interference filtering scanning method includes: determining the range and increment of the common shot gather based on the intermediate seismic data; According to the random splitting method, the common shot set is split into multiple cross subsets; according to the median filtering method, the filter corresponding to each cross subset is determined; according to the filter corresponding to each cross subset, the intermediate seismic data is noise-suppressed to obtain noise-suppressed seismic data; according to the intermediate seismic data and the noise-suppressed seismic data, the first interference noise data is obtained.
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 Taupe mapping domain.
3. The interference noise processing method according to claim 2, characterized in that, The process of performing data domain transformation on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain includes: Based on the initial seismic data, determine the observation coordinates of the seismic waves; The intermediate seismic data is generated based on the slope and intercept of the tangents corresponding to all points in the seismic wave observation coordinates.
4. The interference noise processing method according to claim 3, characterized in that, The process of generating the intermediate seismic data based on the slope and intercept of the tangents corresponding to all points in the seismic wave observation coordinates includes: The linear dynamic calibration time difference threshold is determined based on the slope and intercept of the tangents corresponding to all the points. The intermediate seismic data is generated only when the linear dynamic time difference range corresponding to the seismic wave observation coordinates is less than or equal to the linear dynamic time difference threshold, based on 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 step of performing a data domain inverse transformation on the first interference noise data to obtain the second interference noise data within the first data domain includes: Based on the slope and intercept of the noise data in the first interference noise data, determine the position coordinates of the noise data in the first data domain; The second interference noise data is generated based on the location 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 within the first data domain based on the slope and intercept of the noise data in the first interference noise data, the method further includes: Obtain the threshold number of the common shot gather and the actual number of the first interference noise data; Only when the actual channel number is greater than or equal to the channel number threshold, the position coordinates of the noise data in the first data domain are determined based on the slope and intercept of the noise data in the first interference noise data.
7. A device for processing interference and noise in seismic data, characterized in that, include: The first module is used to perform data domain transformation processing on the initial seismic data in the first data domain to obtain intermediate seismic data in the second data domain. The initial seismic data is seismic data obtained from marine seismic data acquisition; The second module is used to determine the first interference noise data in the second data domain based on the intermediate seismic data and according to a preset interference filtering scanning method. The third module is used to perform data domain inverse transformation on the first interference noise data to obtain the second interference noise data in the first data domain. The fourth module is used to process the initial seismic data based on the second interference noise data to obtain processed effective seismic data. The interference filtering scanning method includes a random splitting method and a median filtering method; the second module is further configured to determine the range and increment of the common shot set based on the intermediate seismic data; split the common shot set into multiple cross subsets according to the random splitting method; determine the filter corresponding to each cross subset according to the median filtering method; perform noise suppression on the intermediate seismic data according to the filter corresponding to each cross subset to obtain noise-suppressed seismic data; and obtain the first interference noise data based on the intermediate seismic data and the noise-suppressed seismic data.
8. An electronic device, characterized in that, It includes a memory and a processor; the memory stores a computer program that, when executed by the processor, implements the method as described in any one of claims 1 to 6.
9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by one or more processors, implements the method as described in any one of claims 1 to 6.