CN114740530A - Medium-high frequency quasi-linear noise suppression method and device based on hyperbolic time window constraint - Google Patents
Medium-high frequency quasi-linear noise suppression method and device based on hyperbolic time window constraint Download PDFInfo
- Publication number
- CN114740530A CN114740530A CN202110020312.9A CN202110020312A CN114740530A CN 114740530 A CN114740530 A CN 114740530A CN 202110020312 A CN202110020312 A CN 202110020312A CN 114740530 A CN114740530 A CN 114740530A
- Authority
- CN
- China
- Prior art keywords
- medium
- high frequency
- linear noise
- single shot
- frequency quasi
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000001629 suppression Effects 0.000 title claims description 56
- 238000001914 filtration Methods 0.000 claims abstract description 21
- 238000004590 computer program Methods 0.000 claims description 16
- 238000009826 distribution Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 8
- 238000003825 pressing Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000000926 separation method Methods 0.000 claims description 4
- 230000006835 compression Effects 0.000 abstract description 4
- 238000007906 compression Methods 0.000 abstract description 4
- 230000000007 visual effect Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 13
- 230000003044 adaptive effect Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
The invention discloses a method and a device for suppressing medium-high frequency quasi-linear noise based on hyperbolic time window constraint, wherein the method comprises the following steps: fitting the medium-frequency and high-frequency quasi-linear noise by utilizing hyperbolic characteristics of a seismic reflection longitudinal wave travel time equation, and then separating single shot records; extracting medium-high frequency quasi-linear noise by utilizing frequency division scanning and mixing processing, and further processing by utilizing an FK domain multi-channel dip filtering method according to the medium-high frequency quasi-linear noise and the medium-high frequency quasi-linear noise to obtain a single shot record after the medium-high frequency quasi-linear noise is finally suppressed; and obtaining the complete single shot record after the final compression of the medium-high frequency quasi-linear noise according to the single shot record without the medium-high frequency quasi-linear noise and the single shot record after the final compression. The invention can comprehensively consider the difference between the medium-high frequency quasi-linear noise and the effective signal in frequency, amplitude and visual speed to achieve the purpose of suppressing the medium-high frequency quasi-linear noise and improve the signal-to-noise ratio of the seismic data.
Description
Technical Field
The invention relates to the technical field of exploration geophysics, in particular to a method and a device for suppressing medium-high frequency quasi-linear noise based on hyperbolic time window constraint.
Background
This section is intended to provide a background or context to the embodiments of the invention that are recited in the claims. The description herein is not admitted to be prior art by inclusion in this section.
When land seismic data are collected, when strong absorption attenuation media such as sand dunes, thick loess and the like exist on the ground surface, and a stable reflecting layer with wave impedance difference with a low-speed zone medium exists in a certain depth range of the shallow layer of the ground surface, the collected seismic data often generate shallow layer noise (hereinafter referred to as medium-high frequency quasi-linear noise) which has energy weaker than that of surface wave, but is stronger than normal effective reflected wave energy, has a speed range of 1000m/s to 2000m/s, has a frequency distribution range of more than 20Hz and is similar to linear interference characteristics. The existence of the medium-high frequency quasi-linear noise seriously affects the signal-to-noise ratio of shallow data, increases the difficulty of subsequent data analysis, and particularly covers the effective signal in-phase axis on a superposed section. The existing linear noise suppression technology generally performs linear noise suppression according to the difference of linear noise and effective signal in visual velocity, position and energy, and according to the difference of linear noise and effective signal between different domains. Because the frequency range of the medium-high frequency quasi-linear noise is highly overlapped with the effective signal, and the apparent velocity is equivalent to the shallow effective wave velocity, the conventional single linear noise suppression method is difficult to obtain a good noise suppression effect.
Disclosure of Invention
The embodiment of the invention provides a hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression method, which is used for solving the technical problem that a conventional single linear noise suppression method is difficult to obtain a good noise suppression effect and comprises the following steps:
acquiring a single shot record;
fitting medium-high frequency quasi-linear noise in the single shot record by utilizing hyperbolic characteristics of an earthquake reflection longitudinal wave travel time equation, determining a root-mean-square speed interval, and separating the single shot record by utilizing a sector area formed by the root-mean-square speed interval to obtain the single shot record containing medium-high frequency quasi-linear noise and the single shot record not containing medium-high frequency quasi-linear noise;
performing frequency division scanning on the single shot record containing the medium-high frequency quasi-linear noise, determining the dominant frequency range of the medium-high frequency quasi-linear noise, extracting the medium-high frequency quasi-linear noise from the single shot record containing the medium-high frequency quasi-linear noise according to the dominant frequency range and mixing processing, and performing self-adaptive subtraction on the single shot record containing the medium-high frequency quasi-linear noise and the extracted medium-high frequency quasi-linear noise to obtain a single shot record for preliminarily suppressing the medium-high frequency quasi-linear noise;
processing the single shot record of preliminarily suppressed medium-high frequency quasi-linear noise by using an FK domain multi-channel dip filtering method to obtain a single shot record of finally suppressed medium-high frequency quasi-linear noise;
and carrying out self-adaptive addition on the single-shot record without the medium-high frequency quasi-linear noise and the single-shot record after the medium-high frequency quasi-linear noise is finally pressed to obtain the complete single-shot record after the medium-high frequency quasi-linear noise is finally pressed.
The embodiment of the invention also provides a device for suppressing medium-high frequency quasi-linear noise based on hyperbolic time window constraint, which is used for solving the technical problem that the conventional single linear noise suppression method is difficult to obtain good noise suppression effect, and comprises the following steps:
the single shot record acquisition module is used for acquiring a single shot record;
the fitting and separating module is used for fitting the medium-high frequency quasi-linear noise in the single shot record by utilizing the hyperbolic characteristic of the seismic reflection longitudinal wave travel time equation, determining a root-mean-square velocity interval, and separating the single shot record by utilizing a sector area formed by the root-mean-square velocity interval to obtain the single shot record containing the medium-high frequency quasi-linear noise and the single shot record not containing the medium-high frequency quasi-linear noise;
the preliminary suppression module is used for performing frequency division scanning on the single shot record containing the medium-high frequency quasi-linear noise, determining the dominant frequency range of the medium-high frequency quasi-linear noise, extracting the medium-high frequency quasi-linear noise from the single shot record containing the medium-high frequency quasi-linear noise according to the dominant frequency range and the mixing processing, and performing self-adaptive subtraction on the single shot record containing the medium-high frequency quasi-linear noise and the extracted medium-high frequency quasi-linear noise to obtain the single shot record for preliminarily suppressing the medium-high frequency quasi-linear noise;
the first final suppression module is used for processing the single shot record for preliminarily suppressing the medium-high frequency quasi-linear noise by using an FK domain multi-channel inclination filtering method to obtain the single shot record for finally suppressing the medium-high frequency quasi-linear noise;
and the second final pressing module is used for performing self-adaptive addition on the single-shot record without the medium-high frequency quasi-linear noise and the single-shot record after the medium-high frequency quasi-linear noise is finally pressed to obtain the complete single-shot record after the medium-high frequency quasi-linear noise is finally pressed.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, where the processor implements the method for suppressing middle-high frequency quasi-linear noise based on hyperbolic time window constraint when executing the computer program.
An embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium stores a computer program for executing the hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression method.
In the embodiment of the invention, compared with the technical scheme that a good noise suppression effect is difficult to obtain by a conventional single linear noise suppression method in the prior art, the method determines the distribution range of noise signals by fitting the hyperbolic characteristic of a seismic reflection longitudinal wave travel time formula to the medium-high frequency quasi-linear noise, and performs the primary separation of the noise and effective signals according to the range, thereby reducing the damage of a subsequent denoising process to other signals outside the range; according to the difference between the noise and the effective signal in the spectrum characteristic and the amplitude characteristic, the middle-high frequency quasi-linear noise is extracted by using spectrum scanning and mixing processing, and a predicted noise model is established, so that the aim of preliminarily suppressing the middle-high frequency quasi-linear noise is fulfilled; and further constructing a noise prediction model by using an FK domain multi-channel dip filtering technology to finally suppress the medium-frequency and high-frequency quasi-linear noise. The method for suppressing the medium-high frequency quasi-linear noise based on the hyperbolic time window constraint can comprehensively consider the differences of the medium-high frequency quasi-linear noise and effective signals in frequency, amplitude and apparent velocity to achieve the purpose of suppressing the medium-high frequency quasi-linear noise and improve the signal-to-noise ratio of seismic data.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments or the prior art descriptions will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts. In the drawings:
FIG. 1 is a flow chart of a method for suppressing high-frequency quasi-linear noise based on hyperbolic time window constraint according to an embodiment of the present invention;
FIG. 2 is a diagram of single shot records of different frequency bands after frequency division scanning is performed on the single shot records by using a zero-phase band-pass filter in the embodiment of the present invention;
FIG. 3 is a diagram of a single shot record before the medium-high frequency quasi-linear noise is suppressed by applying the present invention in an embodiment of the present invention;
FIG. 4 is a diagram of a single shot record after the suppression of medium-high frequency quasi-linear noise according to the present invention;
FIG. 5 is a graph of medium-high frequency quasi-linear noise records suppressed by the present invention;
FIG. 6 is a cross-sectional view of a stack of the embodiment of the present invention before suppressing the high frequency quasi-linear noise;
FIG. 7 is a cross-sectional view of a stack of signals after suppressing high frequency quasi-linear noise according to an embodiment of the present invention;
fig. 8 is a block diagram of a structure of a medium-high frequency quasi-linear noise suppression device based on hyperbolic time window constraint in the embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
In order to effectively suppress medium-high frequency quasi-linear noise, a medium-high frequency quasi-linear noise suppression method based on hyperbolic time window constraint is invented, and the purposes of effectively suppressing medium-high frequency quasi-linear noise and improving the signal-to-noise ratio of seismic data are achieved by comprehensively utilizing the differences of effective waves and medium-high frequency quasi-linear noise in frequency, amplitude and apparent speed.
Fig. 1 is a flowchart of a method for suppressing medium-high frequency quasi-linear noise based on hyperbolic time window constraint in an embodiment of the present invention, as shown in fig. 1, the method includes:
at step 101, a conventional single shot seismic record (FIG. 3) is entered and flow proceeds to step 102.
In step 102, fitting the medium-high frequency quasi-linear noise in the single shot record by using the hyperbolic characteristic of the seismic reflection longitudinal wave travel time equation, determining a root-mean-square velocity interval, and separating the single shot record by using a sector area formed by the root-mean-square velocity interval to obtain the single shot record containing the medium-high frequency quasi-linear noise and the single shot record without the medium-high frequency quasi-linear noise.
Specifically, the rms speed interval is determined from the rms speed scan, wherein the rms speed interval may be used to describe the noise distribution area.
The formula of the travel time of the seismic reflection longitudinal wave can be expressed as:
wherein t is the two-way travel time of the seismic reflection longitudinal wave0When the double-journey travel is carried out by self-excitation and self-collection of the seismic reflected longitudinal wave, x isDistance, V, of shot point to demodulator probenmoThe velocity is dynamically corrected.
Specifically, a sector interval including middle-high frequency quasi-linear noise can be formed according to the root-mean-square velocity interval and the reflected longitudinal wave travel time equation. According to the sector interval, the original single shot record can be divided into a single shot record containing medium-high frequency quasi-linear noise and a single shot record containing no medium-high frequency quasi-linear noise. The flow proceeds to step 103.
In step 103, the separated single shot record containing the medium-high frequency quasi-linear noise is used as input data, and frequency division scanning is performed on the input data by using a zero-phase band-pass filter, so that the dominant frequency range of the medium-high frequency quasi-linear noise is determined, and fig. 2 is a frequency division single shot record obtained by performing frequency division scanning on the single shot seismic record by using the zero-phase band-pass filter. Then, the high-frequency quasi-linear noise is preliminarily extracted from the single shot record containing the medium-high-frequency quasi-linear noise by using a band-pass filter according to the dominant frequency range. And then, analyzing the amplitude characteristics of the high-frequency quasi-linear noise and the effective signal on the basis of the extracted noise data, and further extracting the noise data by utilizing mixing processing according to the difference of the noise and the effective signal in the amplitude characteristics, thereby preliminarily establishing a prediction model of the medium-high frequency quasi-linear noise. The flow proceeds to step 104. The effective signal is the single shot record with the medium-high frequency quasi-linear noise removed from the single shot record with the medium-high frequency quasi-linear noise.
Wave mixing treatment: and performing band-pass filtering in different frequency band ranges on the single shot record containing the medium-high frequency quasi-linear noise to judge whether the signal in a certain frequency range is mainly the medium-high frequency quasi-linear noise signal or other signals except the medium-high frequency quasi-linear noise, and counting the amplitude difference of the two signals. According to the amplitude difference in different frequency band ranges, medium-high frequency quasi-linear noise is extracted.
In step 104, the single shot record containing the medium-high frequency quasi-linear noise obtained in step 102 and the prediction model of the medium-high frequency quasi-linear noise obtained in step 103 are used as input, and the mode matching adaptive subtraction technology is applied to perform adaptive subtraction on the input single shot record to obtain a single shot record after the medium-high frequency quasi-linear noise is preliminarily suppressed. In the pattern-matched adaptive subtraction technique, the data in each time window is further subdivided in the temporal and spatial directions, and a non-stationary Prediction Error Filter (PEF) in the spatial-temporal domain is calculated for each subdivided data. First, the PEF (hereinafter referred to as PEFm) of the noise preliminary prediction model is calculated. PEFm is then convolved with the input noisy single shot record to obtain a rough estimate of the valid signal. This rough estimate is then used to compute the PEF (called PEFp) of the valid signal. The final output valid signal is then estimated from PEFm, PEFp and the input single shot record. The pattern matching adaptive subtraction technique can be expressed as:
wherein, p is the final output effective signal; n is PEF of initial input model data; ' is conjugate transpose; p is PEF of initial input single shot record; eta is a noise attenuation coefficient, and the larger the value is, the heavier the noise suppression is; d is the initially entered single shot record. The flow proceeds to step 105.
In step 105, according to the difference between the apparent velocity of the medium-high frequency quasi-linear noise and the effective signal and the defect of the conventional FK domain filter is avoided, the method adopts an FK domain multi-channel inclination filtering technology to test the dominant parameter and further extract a medium-high frequency quasi-linear noise model. FK domain high-pass dip filter H1(kxY, ω) is expressed as:
wherein k isxWave number of X axis, alpha is 24/3πfnvnFor undetermined coefficients related to apparent velocity, inclination, dominant frequency, fnVideo frequency, v, being noisenApparent velocity as noise; i is an imaginary unit, ω is an angular frequency, y is an offset distance along the crossline direction, and n is a sampling point number.
And carrying out Fourier transformation on the single shot record subjected to the preliminary medium-high frequency quasi-linear noise suppression, and carrying out filtering processing on the single shot record subjected to the Fourier transformation and subjected to the preliminary medium-high frequency quasi-linear noise suppression by using an FK domain high-pass filter, so as to extract medium-high frequency quasi-linear noise signals and construct a medium-high frequency quasi-linear noise model again. The flow proceeds to step 106.
In step 106, the single shot record after the preliminary suppression of the medium-high frequency quasi-linear noise obtained in step 104 and the medium-high frequency quasi-linear noise prediction model obtained in step 105 are used as input, and the two are subjected to self-adaptive subtraction to obtain the single shot record after the final suppression of the medium-high frequency quasi-linear noise. The flow proceeds to step 107.
In step 107, the single shot record obtained in step 102 and the single shot record obtained in step 106 after the final suppression of the medium-high frequency quasi-linear noise are taken as input, and the two are added in a self-adaptive manner to obtain a complete single shot record (fig. 4) after the final suppression of the medium-high frequency quasi-linear noise.
In order to illustrate the processing effect of the invention, the original single-shot seismic record in step 101 is subtracted from the single-shot seismic record after the final suppression of the medium-high frequency quasi-linear noise in step 107, so as to obtain the medium-high frequency quasi-linear noise suppressed by the invention (fig. 5). Comparing fig. 3, 4 and 5, it can be seen that the medium-high frequency quasi-linear noise is well suppressed and no effective signal loss is seen. In order to further verify the processing effect of the invention, the seismic data before and after the application of the invention are stacked under the condition of applying the same stacking flow and parameters, and the stacking results before and after the application of the invention are respectively shown in fig. 6 and fig. 7. From the view of the superimposed profile, the arc noise caused by the existence of the medium-high frequency quasi-linear noise is well suppressed, the effective signal is highlighted, and the signal-to-noise ratio of the superimposed profile is greatly improved.
The embodiment of the invention also provides a medium-high frequency quasi-linear noise suppression device based on hyperbolic time window constraint, and the device is described in the following embodiment. The problem solving principle of the device is similar to that of the hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression method, so the implementation of the device can refer to the implementation of the hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression method, and repeated parts are not repeated.
Fig. 8 is a block diagram of a structure of a medium-high frequency quasi-linear noise suppression device based on hyperbolic time window constraint in an embodiment of the present invention, and as shown in fig. 8, the device includes:
a single shot record obtaining module 02 for obtaining a single shot record;
the fitting and separating module 04 is used for fitting the medium-high frequency quasi-linear noise in the single shot record by utilizing the hyperbolic characteristic of the seismic reflection longitudinal wave travel time equation, determining a root-mean-square velocity interval, and separating the single shot record by utilizing a sector area formed by the root-mean-square velocity interval to obtain the single shot record containing the medium-high frequency quasi-linear noise and the single shot record not containing the medium-high frequency quasi-linear noise;
the preliminary pressing module 06 is used for performing frequency division scanning on the single shot records containing the medium-high frequency quasi-linear noise, determining the dominant frequency range of the medium-high frequency quasi-linear noise, extracting the medium-high frequency quasi-linear noise from the single shot records containing the medium-high frequency quasi-linear noise according to the dominant frequency range and mixing processing, and performing self-adaptive subtraction on the single shot records containing the medium-high frequency quasi-linear noise and the extracted medium-high frequency quasi-linear noise to obtain the single shot records for preliminarily pressing the medium-high frequency quasi-linear noise;
the first final suppression module 08 is used for processing the single shot record for preliminarily suppressing the medium-high frequency quasi-linear noise by using an FK domain multi-channel inclination filtering method to obtain the single shot record for finally suppressing the medium-high frequency quasi-linear noise;
and the second final suppression module 10 is configured to perform adaptive addition on the single shot record without the medium-high frequency quasi-linear noise and the single shot record after the medium-high frequency quasi-linear noise is finally suppressed to obtain a complete single shot record after the medium-high frequency quasi-linear noise is finally suppressed.
In an embodiment of the present invention, the fitting separation module is specifically configured to:
determining a sector area by combining the root-mean-square velocity interval with an earthquake reflected longitudinal wave travel time equation, wherein the sector area completely contains the distribution range of medium-high frequency quasi-linear noise on the single shot record, and the single shot record without the medium-high frequency quasi-linear noise is outside the sector area;
and separating the single shot records according to the sector area.
In an embodiment of the invention, the preliminary press module is specifically configured to:
and performing frequency division scanning on the single shot record containing the medium-high frequency quasi-linear noise by using a zero-phase band-pass filter.
In an embodiment of the invention, the preliminary press module is specifically configured to:
according to the dominant frequency range, preliminarily extracting high-frequency quasi-linear noise from the single shot record containing the medium-high frequency quasi-linear noise by using band-pass filtering;
and analyzing the amplitude characteristics of the high-frequency quasi-linear noise and the effective signal, combining the frequency distribution characteristics and the amplitude characteristics of the medium-high frequency quasi-linear noise, and further extracting the noise by utilizing mixing processing, wherein the effective signal is the single shot record from which the medium-high frequency quasi-linear noise is removed in the single shot record containing the medium-high frequency quasi-linear noise.
In an embodiment of the invention, the first final press module is specifically configured to:
carrying out Fourier transformation on the single shot record subjected to the preliminary compression of the medium-high frequency quasi-linear noise, and carrying out filtering processing on the single shot record subjected to the Fourier transformation and subjected to the preliminary compression of the medium-high frequency quasi-linear noise by using an FK domain high-pass filter to extract a medium-high frequency quasi-linear noise signal;
and carrying out self-adaptive subtraction on the single shot record of the preliminarily pressed medium-high frequency quasi-linear noise and the medium-high frequency quasi-linear noise extracted by utilizing the FK domain multi-channel inclination filtering to obtain the single shot record of the finally pressed medium-high frequency quasi-linear noise.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor executes the computer program to realize the hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression method.
An embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium stores a computer program for executing the hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression method.
In the embodiment of the invention, compared with the technical scheme that a good noise suppression effect is difficult to obtain by a conventional single linear noise suppression method in the prior art, the method determines the distribution range of noise signals by fitting the hyperbolic characteristic of a seismic reflection longitudinal wave travel time formula to the medium-high frequency quasi-linear noise, and performs the primary separation of the noise and effective signals according to the range, thereby reducing the damage of a subsequent denoising process to other signals outside the range; according to the difference between the noise and the effective signal in the spectrum characteristic and the amplitude characteristic, the middle-high frequency quasi-linear noise is extracted by using spectrum scanning and mixing processing, and a predicted noise model is established, so that the aim of preliminarily suppressing the middle-high frequency quasi-linear noise is fulfilled; and further constructing a noise prediction model by using an FK domain multi-channel dip filtering technology to finally suppress the medium-frequency and high-frequency quasi-linear noise. The hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression method can comprehensively consider the differences of medium-high frequency quasi-linear noise and effective signals in frequency, amplitude and visual speed to achieve the purpose of suppressing the medium-high frequency quasi-linear noise and improve the signal-to-noise ratio of seismic data.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (12)
1. A method for suppressing medium-high frequency quasi-linear noise based on hyperbolic time window constraint is characterized by comprising the following steps:
acquiring a single shot record;
fitting medium-high frequency quasi-linear noise in the single shot record by utilizing hyperbolic characteristics of an earthquake reflection longitudinal wave travel time equation, determining a root-mean-square speed interval, and separating the single shot record by utilizing a sector area formed by the root-mean-square speed interval to obtain the single shot record containing medium-high frequency quasi-linear noise and the single shot record not containing medium-high frequency quasi-linear noise;
performing frequency division scanning on the single shot record containing the medium-high frequency quasi-linear noise, determining the dominant frequency range of the medium-high frequency quasi-linear noise, extracting the medium-high frequency quasi-linear noise from the single shot record containing the medium-high frequency quasi-linear noise according to the dominant frequency range and mixing processing, and performing self-adaptive subtraction on the single shot record containing the medium-high frequency quasi-linear noise and the extracted medium-high frequency quasi-linear noise to obtain a single shot record for preliminarily suppressing the medium-high frequency quasi-linear noise;
processing the single shot record of the preliminarily pressed medium-high frequency quasi-linear noise by using an FK domain multi-channel dip filtering method to obtain a single shot record of the finally pressed medium-high frequency quasi-linear noise;
and carrying out self-adaptive addition on the single shot record without the medium-high frequency quasi-linear noise and the single shot record after the medium-high frequency quasi-linear noise is finally suppressed to obtain the complete single shot record after the medium-high frequency quasi-linear noise is finally suppressed.
2. The hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression method according to claim 1, wherein the separating of the single shot records by forming a sector area with the root-mean-square velocity interval comprises:
determining a sector area by utilizing the root-mean-square velocity interval and combining with an earthquake reflection longitudinal wave travel time equation, wherein the sector area completely comprises the distribution range of medium-high frequency quasi-linear noise on the single shot record, and the single shot record without the medium-high frequency quasi-linear noise is outside the sector area;
and separating the single shot records according to the sector area.
3. The hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression method according to claim 1, wherein performing frequency division scanning on a single shot record containing medium-high frequency quasi-linear noise comprises:
and performing frequency division scanning on the single shot record containing the medium-high frequency quasi-linear noise by using a zero-phase band-pass filter.
4. The hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression method according to claim 1, wherein extracting medium-high frequency quasi-linear noise from a single shot record containing medium-high frequency quasi-linear noise according to the dominant frequency range and the mixing process comprises:
according to the dominant frequency range, preliminarily extracting high-frequency quasi-linear noise from the single shot record containing the medium-high frequency quasi-linear noise by using band-pass filtering;
and analyzing the amplitude characteristics of the high-frequency quasi-linear noise and the effective signal, combining the frequency distribution characteristics and the amplitude characteristics of the medium-high frequency quasi-linear noise, and further extracting the noise by utilizing mixing processing, wherein the effective signal is the single shot record from which the medium-high frequency quasi-linear noise is removed in the single shot record containing the medium-high frequency quasi-linear noise.
5. The hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression method according to claim 1, wherein a single shot record of preliminarily suppressed medium-high frequency quasi-linear noise is processed by using an FK domain multi-channel dip filtering method to obtain a single shot record of finally suppressed medium-high frequency quasi-linear noise, and the method comprises:
carrying out Fourier transformation on the single shot record subjected to preliminary medium-high frequency quasi-linear noise suppression, and carrying out filtering processing on the single shot record subjected to Fourier transformation and subjected to preliminary medium-high frequency quasi-linear noise suppression by using an FK domain high-pass filter to extract medium-high frequency quasi-linear noise signals;
and carrying out self-adaptive subtraction on the single shot record of the preliminarily pressed medium-high frequency quasi-linear noise and the medium-high frequency quasi-linear noise extracted by utilizing the FK domain multi-channel inclination filtering to obtain the single shot record of the finally pressed medium-high frequency quasi-linear noise.
6. A kind of medium-high frequency quasi-linear noise suppression device based on hyperbolic time window constraint, characterized by that, including:
the single shot record acquisition module is used for acquiring a single shot record;
the fitting and separating module is used for fitting the medium-high frequency quasi-linear noise in the single shot record by utilizing the hyperbolic characteristic of the seismic reflection longitudinal wave travel time equation, determining a root-mean-square velocity interval, and separating the single shot record by utilizing a sector area formed by the root-mean-square velocity interval to obtain the single shot record containing the medium-high frequency quasi-linear noise and the single shot record not containing the medium-high frequency quasi-linear noise;
the preliminary suppression module is used for performing frequency division scanning on the single shot record containing the medium-high frequency quasi-linear noise, determining the dominant frequency range of the medium-high frequency quasi-linear noise, extracting the medium-high frequency quasi-linear noise from the single shot record containing the medium-high frequency quasi-linear noise according to the dominant frequency range and the mixing processing, and performing self-adaptive subtraction on the single shot record containing the medium-high frequency quasi-linear noise and the extracted medium-high frequency quasi-linear noise to obtain the single shot record for preliminarily suppressing the medium-high frequency quasi-linear noise;
the first final suppression module is used for processing the single shot record for preliminarily suppressing the medium-high frequency quasi-linear noise by using an FK domain multi-channel inclination filtering method to obtain the single shot record for finally suppressing the medium-high frequency quasi-linear noise;
and the second final pressing module is used for performing self-adaptive addition on the single shot record without the medium-high frequency quasi-linear noise and the single shot record after the medium-high frequency quasi-linear noise is finally pressed to obtain the complete single shot record after the medium-high frequency quasi-linear noise is finally pressed.
7. The hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression device according to claim 6, wherein the fitting separation module is specifically configured to:
determining a sector area by utilizing the root-mean-square velocity interval and combining with an earthquake reflection longitudinal wave travel time equation, wherein the sector area completely comprises the distribution range of medium-high frequency quasi-linear noise on the single shot record, and the single shot record without the medium-high frequency quasi-linear noise is outside the sector area;
and separating the single shot records according to the sector area.
8. The hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression device according to claim 6, wherein the preliminary suppression module is specifically configured to:
and performing frequency division scanning on the single shot record containing the medium-high frequency quasi-linear noise by using a zero-phase band-pass filter.
9. The hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression device according to claim 6, wherein the preliminary suppression module is specifically configured to:
according to the dominant frequency range, preliminarily extracting high-frequency quasi-linear noise from single shot records containing medium-high frequency quasi-linear noise by using band-pass filtering;
and analyzing the amplitude characteristics of the high-frequency quasi-linear noise and the effective signal, combining the frequency distribution characteristics and the amplitude characteristics of the medium-high frequency quasi-linear noise, and further extracting the noise by utilizing mixing processing, wherein the effective signal is the single shot record from which the medium-high frequency quasi-linear noise is removed in the single shot record containing the medium-high frequency quasi-linear noise.
10. The hyperbolic time window constraint-based medium-high frequency quasi-linear noise suppression device of claim 6, wherein the first final suppression module is specifically configured to:
carrying out Fourier transformation on the single shot record subjected to preliminary medium-high frequency quasi-linear noise suppression, and carrying out filtering processing on the single shot record subjected to Fourier transformation and subjected to preliminary medium-high frequency quasi-linear noise suppression by using an FK domain high-pass filter to extract medium-high frequency quasi-linear noise signals;
and carrying out self-adaptive subtraction on the single shot record of the preliminarily pressed medium-high frequency quasi-linear noise and the medium-high frequency quasi-linear noise extracted by utilizing the FK domain multi-channel inclination filtering to obtain the single shot record of the finally pressed medium-high frequency quasi-linear noise.
11. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 5 when executing the computer program.
12. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110020312.9A CN114740530B (en) | 2021-01-07 | Medium-high frequency quasi-linear noise suppression method and device based on hyperbolic time window constraint |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110020312.9A CN114740530B (en) | 2021-01-07 | Medium-high frequency quasi-linear noise suppression method and device based on hyperbolic time window constraint |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114740530A true CN114740530A (en) | 2022-07-12 |
| CN114740530B CN114740530B (en) | 2024-06-25 |
Family
ID=
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102681014A (en) * | 2012-05-23 | 2012-09-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Regular linear interference suppressing method based on polynomial fitting |
| CN103869370A (en) * | 2014-03-31 | 2014-06-18 | 中国石油大学(北京) | Method and system for filtering out linear interference of seismic data |
| WO2014195508A2 (en) * | 2013-06-07 | 2014-12-11 | Cgg Services Sa | Systems and methods for de-noising seismic data |
| CN105242309A (en) * | 2015-10-19 | 2016-01-13 | 中国石油天然气集团公司 | Earthquake data regular interference suppressing method and device |
| CN106125134A (en) * | 2016-07-04 | 2016-11-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Based on the geological data signal-noise ratio computation method of window during hyperbolic |
| CN106842323A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing |
| CN111562616A (en) * | 2019-02-14 | 2020-08-21 | 中国石油天然气股份有限公司 | Seismic data scattering noise suppression method and device |
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102681014A (en) * | 2012-05-23 | 2012-09-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Regular linear interference suppressing method based on polynomial fitting |
| WO2014195508A2 (en) * | 2013-06-07 | 2014-12-11 | Cgg Services Sa | Systems and methods for de-noising seismic data |
| CN103869370A (en) * | 2014-03-31 | 2014-06-18 | 中国石油大学(北京) | Method and system for filtering out linear interference of seismic data |
| CN105242309A (en) * | 2015-10-19 | 2016-01-13 | 中国石油天然气集团公司 | Earthquake data regular interference suppressing method and device |
| CN106842323A (en) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | A kind of slip scan harmonic wave disturbance suppression method based on scaling down processing |
| CN106125134A (en) * | 2016-07-04 | 2016-11-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Based on the geological data signal-noise ratio computation method of window during hyperbolic |
| CN111562616A (en) * | 2019-02-14 | 2020-08-21 | 中国石油天然气股份有限公司 | Seismic data scattering noise suppression method and device |
Non-Patent Citations (4)
| Title |
|---|
| JIE ZHANG 等: "Noise Attenuation for Seismic Data by Hyperbolic-Trace Time-Frequency Peak Filtering", IEEE GEOSCIENCE AND REMOTE SENSING LETTERS, vol. 12, no. 3, 16 September 2014 (2014-09-16), pages 601 - 605, XP011560941, DOI: 10.1109/LGRS.2014.2352671 * |
| 孙万军 等: "拟线性空间变换叠前F-X地震去噪新方法", 大庆石油地质与开发, vol. 24, no. 003, 31 December 2005 (2005-12-31), pages 97 - 99 * |
| 杨文成: "提高深层地震资料分辨率方法研究", 中国博士学位论文全文数据库 基础科学辑, no. 2, 15 February 2018 (2018-02-15), pages 011 - 77 * |
| 段中钰: "拟线性噪声的形态滤波方法研究", 北京信息科技大学学报(自然科学版), vol. 31, no. 03, 30 June 2016 (2016-06-30), pages 19 - 22 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108983284B (en) | F-p domain ghost wave compression method suitable for offshore inclined cable data | |
| Chen et al. | Random noise reduction using a hybrid method based on ensemble empirical mode decomposition | |
| US20150168573A1 (en) | Geologic quality factor inversion method | |
| CN112835103B (en) | Adaptive ghost wave removing and broadband quasi-zero phase deconvolution combined processing method and system | |
| CN105242309B (en) | The drawing method and device of rule interference in a kind of geological data | |
| CN104614769B (en) | A kind of Beamforming for suppressing seismic surface wave | |
| Zhou et al. | Robust noise attenuation based on nuclear norm minimization and a trace prediction strategy | |
| CN109633752B (en) | Offshore towing cable data self-adaptive ghost wave compression method based on three-dimensional fast Radon transformation | |
| CN105445801B (en) | A kind of processing method for eliminating 2-d seismic data random noise | |
| CN106707342B (en) | Common-shot-gather multistage surface wave pressing method and device | |
| CN106680876A (en) | Joint denoising method for seismic data | |
| CN111708087A (en) | Method for suppressing seismic data noise based on DnCNN neural network | |
| CN107255833B (en) | The determination method and apparatus of Seismic Stacked Section | |
| CN110618449A (en) | Method and system for processing seismic data | |
| CN114740530B (en) | Medium-high frequency quasi-linear noise suppression method and device based on hyperbolic time window constraint | |
| CN114740530A (en) | Medium-high frequency quasi-linear noise suppression method and device based on hyperbolic time window constraint | |
| CN105954799B (en) | A kind of time-frequency domain Processing Seismic Data based on weighted superposition | |
| CN112764108B (en) | Novel seismic data noise suppression algorithm based on improved empirical wavelet transformation | |
| CN110596756B (en) | Desert seismic exploration noise suppression method based on self-adaptive mixed complex diffusion model | |
| CN110837119B (en) | Seismic data processing method and system for enhancing inverse Q value compensation stability | |
| CN113093282A (en) | Desert data denoising method based on geometric modal characteristic parallel network | |
| CN112363218A (en) | 3D-Shearlet transformation and scale self-adaption based seismic random noise suppression method | |
| CN110737019B (en) | Processing method and system for improving velocity spectrum precision | |
| Zhou et al. | Seismic noise attenuation using an improved variational mode decomposition method | |
| CN112198553B (en) | Man-machine interactive deep sea water body three-dimensional velocity analysis method integrating reflection amplitude compensation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |