CN111487675A - Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum - Google Patents
Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum Download PDFInfo
- Publication number
- CN111487675A CN111487675A CN202010215366.6A CN202010215366A CN111487675A CN 111487675 A CN111487675 A CN 111487675A CN 202010215366 A CN202010215366 A CN 202010215366A CN 111487675 A CN111487675 A CN 111487675A
- Authority
- CN
- China
- Prior art keywords
- time
- spectrum
- signal
- seismic
- frequency spectrum
- 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
- 238000001228 spectrum Methods 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000005311 autocorrelation function Methods 0.000 claims abstract description 18
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 18
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 12
- 238000001914 filtration Methods 0.000 claims abstract description 6
- 230000003044 adaptive effect Effects 0.000 claims abstract description 4
- 239000002131 composite material Substances 0.000 claims description 13
- 230000003595 spectral effect Effects 0.000 claims description 5
- 230000001276 controlling effect Effects 0.000 claims description 4
- 238000003384 imaging method Methods 0.000 claims description 4
- 230000001105 regulatory effect Effects 0.000 claims description 2
- 239000003245 coal Substances 0.000 abstract description 8
- 239000011159 matrix material Substances 0.000 description 11
- 239000013598 vector Substances 0.000 description 9
- 238000005457 optimization Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000001629 suppression Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000926 separation method 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
-
- 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/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a method for generating a high-signal-to-noise-ratio high-resolution time frequency spectrum of seismic data, which is realized by multi-channel weighted synthesis along a formation dip angle and a mode of constructing an adaptive mask filter; the method comprises the following steps: estimating dip angles of each time-space position on the seismic section; constructing a plurality of weighted synthesis channels aiming at each time position, and calculating an instantaneous autocorrelation function of an analytic signal of the synthesis channels; calculating standard Weiwei time-frequency distribution and corresponding flat-virtual Weiwei distribution; constructing a time-frequency spectrum self-adaptive mask filter; and realizing mask filtering of standard Weiwei time-frequency distribution to obtain an optimized seismic channel time-frequency spectrum. The method has the characteristics of high signal-to-noise ratio and high resolution, simultaneously has the capability of eliminating cross interference, and has the advantage of clearly separating the spatial distribution of the sandstone reservoir and the spatial distribution of the system containing the coal bed when the actual three-dimensional seismic data is used for reservoir prediction.
Description
Technical Field
The invention relates to the field of geophysical signal analysis, in particular to a method for generating a seismic data time frequency spectrum with high signal-to-noise ratio and high resolution.
Background
Time-frequency spectrum analysis of seismic data is a key means for reservoir geophysical research. Intuitively, the time-frequency spectrum analysis is to generate corresponding frequency spectrum point by point according to the time sampling sequence of any seismic single channel, and is to convert a one-dimensional seismic channel changing along with time into a two-dimensional spectrum changing along with time and frequency. This two-dimensional spectrum is called the time spectrum.
The Wigner-Ville time frequency distribution method is one of effective methods for calculating the time frequency spectrum of the seismic trace. The implementation process can be divided into three basic steps: calculating analytic signals of the seismic traces, calculating an instantaneous autocorrelation function of the analytic signals, and calculating Fourier transform of the autocorrelation function. Because the analytical signal and the instantaneous autocorrelation function are adopted in the calculation method of the Weiwei time-frequency distribution, the Weiwei time-frequency distribution has the characteristic of high resolution. However, it is because the computation of the instantaneous autocorrelation function leads to severe cross-interference between the components contained in the signal.
In order to suppress the cross-interference between the above-mentioned signal components, researchers have proposed a calculation method called the flat-imaginary wiry distribution. It contains two layers of meaning, one "smooth" refers to the suppression of interference in the time direction, and one "virtual" actually refers to the suppression of interference in the frequency direction. However, the suppression processing in both the time direction and the frequency direction is essentially low-pass filtering, and this has the effect of losing high-frequency components and reducing the resolution of the time spectrum.
Because the time frequency spectrum calculation method is a single-channel analysis method, the generated time frequency spectrum is influenced by the signal-to-noise ratio of seismic data, the signal-to-noise ratio of the time frequency spectrum is low, and the time frequency spectrum between channels lacks space consistency, thereby causing difficulty in the interpretation and application of the seismic time frequency spectrum.
By combining the analysis, the spectrum resolution is high but the cross interference is serious in the standard Weiwei time spectrum, the cross interference is inhibited but the resolution is reduced in the flat-virtual Weiwei time spectrum, and the calculation of the time spectrum is single-channel operation, so that the generated earthquake time spectrum has poor signal-to-noise ratio and is seriously influenced by earthquake data noise, and the application of the time spectrum and the geophysical explanation are difficult.
Disclosure of Invention
The invention aims to overcome the defects of the existing method and provides a method for generating a high-signal-to-noise-ratio high-resolution time spectrum of seismic data. The method for generating the earthquake time frequency spectrum has the advantages of high signal-to-noise ratio and high resolution.
In order to achieve the purpose, the seismic data high signal-to-noise ratio high resolution time frequency spectrum generation method provided by the invention improves the signal-to-noise ratio of the seismic time frequency spectrum by multi-channel weighted synthesis along the formation dip angle and utilizing the spatial continuity of seismic reflection homophase axes; the resolution of the frequency spectrum during earthquake is enhanced by a method of constructing an adaptive mask filter. The method has the advantages of high resolution and high signal-to-noise ratio.
As optimization, the seismic data high signal-to-noise ratio high resolution time frequency spectrum generation method specifically comprises the following steps: (1) estimating dip angles of each time-space position on the seismic section; (2) constructing a plurality of weighted synthesis channels aiming at each time position, and calculating an instantaneous autocorrelation function of an analytic signal of the synthesis channels; (3) calculating standard Weiwei (Wigner-Ville) time-frequency distribution and corresponding flat-virtual Weiwei distribution; (4) constructing a time-frequency spectrum self-adaptive mask filter; (5) and realizing mask filtering of standard Weiwei time-frequency distribution to obtain an optimized seismic channel time-frequency spectrum.
As optimization, the step (1) of estimating the dip angle of each time-space position from the actual seismic section is to construct a nonlinear equation set by adopting a discrete plane wave analytic equation mode, and estimate the time-varying and space-varying stratigraphic dip angles by a calculation method for solving a nonlinear problem.
Estimating the dip angle of each time-space position on the seismic section by adopting a plane wave analytic equation:
in the formulatIs a time sample of the seismic trace,xis the spatial location of the seismic traces,u(t,x) Is the seismic wavefield and σ is the dip of the linear event. Discretizing the plane wave analytic equation and expressing the plane wave analytic equation in a vector matrix form
f(σ) = e T [C(σ)○U]e = 0,
Where U is the data matrix, C (σ) is the matrix composed of the coefficients of the two-dimensional filter in the reverse order, [ C (σ) ○ U]Is a multiplication by elements or a Hadamard product of two matrices, e is a full 1 vector. The matrix [ C ] is realized by introducing all 1 vectors e in the formula(σ)○U]The sum of all elements in. By solving a system of non-linear equationsf(σ) ≈ 0 estimates time-varying and space-varying formation dip angles.
As an optimization, step (2) implements the following three steps at each time position: constructing a composite track by performing multi-channel weighting according to local dip anglesx(t) (ii) a Calculating an analytic signal of the composite tracez(t) (ii) a Calculating the instantaneous autocorrelation function of the analytic signal:k(t,τ) =z(t+τ/2)z*(t-τ/2) in the formulaτIs the lag time of the instantaneous auto-correlation,z*(t) Is analyzing the signalz(t) The complex number of (c) is conjugated.
As optimization, step (3) uses the instantaneous autocorrelation function obtained in step (2) to calculate the standard wiry time-frequency distribution as follows:
and the corresponding Ping-Xunwei distribution was calculated as follows:
in the formulah(τ) Is a symmetric function with time as a variable, which functions as a low-pass filter along the frequency axis,g(t) Is another function of time as a variable, which functions as a low-pass filter along the time axis, andg(0) =h(0)。
as optimization, the mask filter in step (4) is constructed as follows:
the middle ratio ofW SP (t,f)|2/ |W(t,f)|2Is the similarity of the two power spectra,is an adjustment parameter that controls the position of the mask filter boundary,ηthen it is the parameter that controls the smoothness of the mask filter boundary, settingηNot less than 2. For the sake of stability of the estimation, int,f) The power spectrum estimate at a point is the sum of the power spectra within a small two-dimensional window centered at that point.
As optimization, the constructed mask filter has the characteristics of the flattest passband and the adjusting parameters thereofEffectively controlling the boundary position of the filter; setting a boundary smoothness parameterη= 4 regulating parameterWhen = 0.05, any imaging spectral point with a power spectrum ratio value less than 20% can be filtered out from the time-frequency spectrum.
As optimization, the implementation process of the step (5) is divided into two steps: distribution of amplitude spectrum of Weiwei time-frequency distributionW(t,f) The l is filtered out and the filter is implemented,
(ii) a The wiry time frequency spectrum is reconstructed,
in the formulaƟ(t,f) The method is an original phase spectrum of Weiwei time-frequency distribution, and finally obtains a seismic channel optimized time-frequency spectrum with high signal-to-noise ratio and high resolution.
After the embodiment is adopted, the invention has the following beneficial effects: (1) the method has the advantages that the method has high resolution, the generated time frequency spectrum has the same high resolution as that of standard Weiwei distribution, and the method has the same capability of eliminating cross interference as that of flat-virtual Weiwei distribution. (2) The method has the advantage of high signal-to-noise ratio, and can clearly separate the spatial distribution of the sandstone reservoir and the spatial distribution of the coal-containing bed when the actual three-dimensional seismic data is used for reservoir prediction.
Drawings
FIG. 1 is a flow chart of a basic implementation of the seismic data high signal-to-noise ratio and high resolution time-frequency spectrum generation method of the present invention; FIG. 2 is a diagram of an actual seismic section (top) and its estimated dip (bottom) according to an embodiment of the present invention.
Fig. 3 is a mask filter proposed by the present invention, having the flattest passband characteristic. The horizontal axis in the diagram is the ratio of two power spectra +W SP (t,f)|2/|W(t,f)|2Calculating the boundary smoothness parameterηThe value is 4, and the adjusting parameters corresponding to the coefficient curves of the three mask filtersThe values are 0.1,0.05,0.01, respectively.
FIG. 4 is a diagram of an embodiment of the present invention applied to a composite trace containing 4 wavelets comparing the standard Weir distribution of the composite trace, the Ping-Vwei distribution of the composite trace, and the time-frequency spectrum generated by the method of the present invention. Fig. 4 shows that the optimized temporal spectrum generated by the method of the present invention has the same high resolution as the standard wiry distribution, while the method has the same cross-interference rejection capability as the flat-imaginary wiry distribution.
FIG. 5 is an embodiment of the present invention applied to actual three-dimensional seismic recordings for reservoir prediction. The top plot is a time slice of a three-dimensional seismic data volume. The spectrum slice of the middle graph reflects the spatial distribution of the coal measure strata of the region, and the primary frequency of the reflection earthquake of the coal measure strata of the region is 15 Hz. The spectral slice of the lower graph is the spatial distribution of the sandstone reservoir in the region, which reflects seismic dominant frequencies of 20 hertz. The figure shows the high signal-to-noise ratio characteristic of the method of the invention, enabling a clear separation of the spatial distribution of the sandstone reservoir and the coal-derived formation.
Detailed Description
The method for generating the seismic data high signal-to-noise ratio and high resolution time frequency spectrum comprises the steps of improving the signal-to-noise ratio of the seismic time frequency spectrum by means of multi-channel weighted synthesis along the dip angle of the stratum and utilizing the spatial continuity of seismic reflection in-phase axes; the resolution of the frequency spectrum during earthquake is enhanced by a method of constructing an adaptive mask filter. The method comprises the following specific steps: (1) estimating dip angles of each time-space position on the seismic section; (2) constructing a plurality of weighted synthesis channels aiming at each time position, and calculating an instantaneous autocorrelation function of an analytic signal of the synthesis channels; (3) calculating standard Weiwei (Wigner-Ville) time-frequency distribution and corresponding flat-virtual Weiwei distribution; (4) constructing a time-frequency spectrum self-adaptive mask filter; (5) and realizing mask filtering of standard Weiwei time-frequency distribution to obtain an optimized seismic channel time-frequency spectrum.
In the step (1), the estimation of the dip angle of each time-space position from the actual seismic section is to construct a nonlinear equation set by discretizing a plane wave analytic equation, and estimate the dip angles of the time-varying and space-varying stratums by a calculation method for solving a nonlinear problem. Estimating the dip angle of each time-space position on the seismic section by adopting a plane wave analytic equation:
in the formulatIs the point in time of the seismic trace,xis the spatial location of the seismic traces,u(t,x) Is the seismic wavefield and σ is the dip of the linear event. Discretizing the plane wave analytic equation and expressing the plane wave analytic equation in a vector matrix formf(σ)=e T [C(σ)○U]e =0, where U is the data matrix and C (σ) is the matrix composed of the coefficients of the two-dimensional filter in the reverse order, with dimensions identical to the data matrix U, [ C (σ) ○ U]Is element-wise multiplication or Hadamard multiplication of two matrices, e being a full "1" vector, introduction of the full "1" vector e in the above formula achieves a pair of matrices [ C (σ) ○ U]The sum of all elements in. By solving a system of non-linear equationsf(σ) ≈ 0 estimates time-varying and space-varying formation dip angles.
The step (2) realizes the following three steps at each time position: constructing a composite track by performing multi-channel weighting according to local dip anglesx(t) (ii) a Calculating an analytic signal of the composite tracez(t) (ii) a Calculating the instantaneous autocorrelation function of the analytic signal:k(t,τ)=z(t+τ/2)z*(t-τ/2) In the formulaτIs the lag time of the autocorrelation function,z*(t) Is analyzing the signalz(t) The complex number of (c) is conjugated.
And (3) calculating the standard Weiwei time-frequency distribution by using the autocorrelation function in the step (2) as follows:
and the corresponding Ping-Xunwei distribution was calculated as follows:
in the formulah(τ) Is a symmetric function with time as a variable, which functions as a low-pass filter along the frequency axis,g(t) Is another function of time as a variable, which functions as a low-pass filter along the time axis, andg(0) =h(0)。
the mask filter in the step (4) is constructed as follows:
the middle ratio ofW SP (t,f)|2/|W(t,f)|2Is the similarity of the two power spectra,is an adjustment parameter that controls the position of the mask filter boundary,ηthen it is the parameter that controls the smoothness of the mask filter boundary, settingηNot less than 2. For the sake of stability, int,f) The amplitude estimate at a point is the sum of the amplitudes within a small two-dimensional window centered at that point. The constructed mask filter has the characteristics of the flattest passband and the adjusting parameterEffectively controlling the boundary position of the filter; setting a boundary smoothness parameterηAdjustment parameter with value 4When the value is 0.05, any imaging spectrum point with the power spectrum ratio of less than 20 percent can be filtered from the time-frequency spectrum.
The implementation process of the step (5) is divided into two steps: distribution of amplitude spectrum of Weiwei time-frequency distributionW(t,f) The l is filtered out and the filter is implemented,
(ii) a The wiry time frequency spectrum is reconstructed,
in the formulaƟ(t,f) The method is an original phase spectrum of Weiwei time-frequency distribution, and finally obtains a seismic channel optimized time-frequency spectrum with high signal-to-noise ratio and high resolution.
In order to make the technical scheme and advantages of the invention more apparent, the following detailed description of the embodiments of the invention is provided in conjunction with the accompanying drawings. The invention provides a method for generating a high-signal-to-noise-ratio high-resolution time spectrum of seismic data. FIG. 1 is a flow chart of a basic implementation of the method of the present invention.
Step S101, estimating the dip angle of each time-space position on the seismic section. Using a plane wave analytic equation:
in the formulatIs the point in time of the seismic trace,xis the spatial location of the seismic traces,u(t,x) Is the seismic wavefield and σ is the dip of the linear event. Discretizing the plane wave analytic equation and expressing the plane wave analytic equation in a vector matrix form
f(σ) = e T [C(σ)○U]e = 0
Where U is the data matrix, C (σ) is the matrix composed of the coefficients of the two-dimensional filter in the reverse order, [ C (σ) ○ U]Is element-wise multiplication or Hadamard product of two matrices, e is a full 1 vector, and introducing the full 1 vector e in the above formula realizes the pair of matrices [ C (sigma) ○ U]The sum of all elements in. By solving a system of non-linear equationsf(σ) ≈ 0 estimates time-varying and space-varying formation dip angles.
FIG. 2 is an actual seismic section (top) and its estimated dip (bottom) of an embodiment of the invention. The arrows in the figure are the estimated local dip directions, coinciding with the true in-phase axes of the actual seismic event of the background.
Step S102, the following three steps are realized at each time position: constructing a composite track by performing multi-channel weighting according to local dip anglesx(t) (ii) a Calculating an analytic signal of the composite tracez(t) (ii) a Calculating the instantaneous autocorrelation function of the analytic signal:k(t,τ)=z(t+τ/2)z*(t-τ/2) In the formulaτIs the lag time of the autocorrelation function,z*(t) Is analyzing the signalz(t) The complex number of (c) is conjugated.
Step S103, calculating the standard Weiwei time-frequency distribution by using the instantaneous autocorrelation function obtained in the step S102 as follows:
and the corresponding Ping-Xunwei distribution was calculated as follows:
in the formulah(τ) Is a symmetric function with time as a variable, which functions as a low-pass filter along the frequency axis,g(t) Is another function of time as a variable, which functions as a low-pass filter along the time axis, andg(0) =h(0)。
step S104, the mask filter is constructed as follows:
the middle ratio ofW SP (t,f)|2/|W(t,f)|2Is the similarity of the two power spectra,is an adjustment parameter that controls the position of the mask filter boundary,ηthen it is the parameter that controls the smoothness of the mask filter boundary, settingηNot less than 2. For the sake of stability, int,f) The amplitude estimate at a point is the sum of the amplitudes within a small two-dimensional window centered at that point. As shown in FIG. 3, the constructed mask filter has the flattest passband characteristic, which adjusts the parametersEffectively controlling the boundaries of the filter; setting a boundary smoothness parameterηAdjustment parameter with value 4When the value is 0.05, any imaging spectrum point with the power spectrum ratio of less than 20 percent can be filtered from the time-frequency spectrum.
And step S105, realizing mask filtering of the standard Weiwei time-frequency distribution. The implementation process comprises two steps: distribution of amplitude spectrum of Weiwei time-frequency distributionW(t,f) The l is filtered out and the filter is implemented,
(ii) a The wiry time frequency spectrum is reconstructed,
in the formulaƟ(t,f) The method is an original phase spectrum of Weiwei time-frequency distribution, and finally obtains a seismic channel optimized time-frequency spectrum with high signal-to-noise ratio and high resolution.
The invention also provides two embodiments [ fig. 4, fig. 5], which illustrate that the invention has the following beneficial effects:
(1) as shown in FIG. 4, a composite trace contains 4 wavelets, and in its standard Weir distribution, the cross-interference between wavelets is very significant except for the energy blob corresponding to 4 wavelets; on the flat-virtual weiwei distribution of the composite trace, although the cross interference is completely eliminated, 4 effective energy blobs are obviously spread, so that the resolution of the time spectrum is lowered; finally, the time frequency spectrum generated by the method has the high resolution as that of the standard Weiwei distribution, and meanwhile, the method has the capability of eliminating cross interference as that of the flat-virtual Weiwei distribution.
(2) FIG. 5 is an embodiment of the present invention applied to actual three-dimensional seismic recordings for reservoir prediction. The top plot is a time slice of a three-dimensional seismic data volume. The spectral slice of the middle graph reflects the spatial distribution of the coal measure strata of the region, and the primary frequency of the reflection earthquake of the coal measure strata of the region is 15 Hz. Note that the coal measure stratigraphic space variation of the middle graph has similarities with some parts of the upper graph. The spectral slice of the lower graph is the spatial distribution of the sandstone reservoir in the region, and the reflection seismic main frequency of the sandstone reservoir in the region is 20 Hz. As shown in fig. 5, the method has the advantage of high signal-to-noise ratio, and can clearly separate the spatial distribution of the sandstone reservoir and the spatial distribution of the coal-series stratum when the actual three-dimensional seismic data is used for reservoir prediction.
Claims (8)
1. A method for generating a high-signal-to-noise-ratio high-resolution time spectrum of seismic data is characterized in that the signal-to-noise ratio of the seismic time spectrum is improved by multi-channel weighted synthesis along the dip angle of a stratum by utilizing the spatial continuity of seismic reflection event axes; the resolution of the frequency spectrum during earthquake is enhanced by a method of constructing an adaptive mask filter.
2. The method for generating the seismic data time spectrum with high signal-to-noise ratio and high resolution according to claim 1, characterized by the following steps: (1) estimating dip angles of each time-space position on the seismic section; (2) constructing a plurality of weighted synthesis channels aiming at each time position, and calculating an instantaneous autocorrelation function of an analytic signal of the synthesis channels; (3) calculating standard Weiwei (Wigner-Ville) time-frequency distribution and corresponding flat-virtual Weiwei distribution; (4) constructing a time-frequency spectrum self-adaptive mask filter; (5) and realizing mask filtering of standard Weiwei time-frequency distribution to obtain an optimized seismic channel time-frequency spectrum.
3. The method for generating high SNR and high resolution time-frequency spectrum of seismic data as claimed in claim 2, wherein the step (1) of estimating the dip angle of each time-space position from the actual seismic section is to construct a nonlinear equation system by using a discrete plane wave analytic equation, and estimate the time-varying and space-varying stratigraphic dip angles by solving a nonlinear problem.
4. The method for generating a high signal-to-noise ratio and high resolution time spectrum of seismic data as claimed in claim 2, wherein the step (2) implements the following three steps at each time position: building a composite track by multi-track weighting according to local dip anglesx(t) (ii) a Calculating an analytic signal of the composite tracez(t) (ii) a Calculating the instantaneous autocorrelation function of the analytic signal:k(t,τ)=z(t+τ/2)z*(t-τ/2) in the formulaτIs the lag time of the instantaneous auto-correlation,z*(t) Is analyzing the signalz(t) The complex number of (c) is conjugated.
5. The method for generating the seismic data time-frequency spectrum with high signal-to-noise ratio and high resolution according to claim 2, wherein the step (3) uses the instantaneous autocorrelation function obtained in the step (2) to calculate the standard wiry time-frequency distribution as follows:
and the corresponding Ping-Xunwei distribution was calculated as follows:
in the formulah(τ) Is a symmetric function with time as a variable, which functions as a low-pass filter along the frequency axis,g(t) Is another function of time as a variable, which functions as a low-pass filter along the time axis, andg(0) =h(0)。
6. the method for generating the seismic data time spectrum with high signal-to-noise ratio and high resolution as claimed in claim 2, wherein the mask filter in the step (4) is constructed as follows:
the middle ratio ofW SP (t,f)|2/ |W(t,f)|2Is the similarity of the two power spectra,is an adjustment parameter that controls the position of the mask filter boundary,ηthen it is the parameter that controls the smoothness of the mask filter boundary, settingη≥ 2。
7. The method of claim 6, wherein the mask filter is constructed to have the flattest passband characteristics and the tuning parameters thereofEffectively controlling the boundary position of the mask filter; setting a boundary smoothness parameterη= 4 regulating parameterWhen = 0.05, any imaging spectral point with a power spectrum ratio value less than 20% can be filtered out from the time-frequency spectrum.
8. The method of claim 2The method for generating the seismic data time spectrum with high signal-to-noise ratio and high resolution is characterized in that the implementation process of the step (5) is divided into two steps: distribution of amplitude spectrum of Weiwei time-frequency distributionW(t,f) The l is filtered out and the filter is implemented,
(ii) a The wiry time frequency spectrum is reconstructed,
in the formulaƟ(t,f) The method is an original phase spectrum of Weiwei time-frequency distribution, and finally obtains a seismic channel optimized time-frequency spectrum with high signal-to-noise ratio and high resolution.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010215366.6A CN111487675B (en) | 2020-03-25 | 2020-03-25 | Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010215366.6A CN111487675B (en) | 2020-03-25 | 2020-03-25 | Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111487675A true CN111487675A (en) | 2020-08-04 |
| CN111487675B CN111487675B (en) | 2021-08-27 |
Family
ID=71812363
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010215366.6A Active CN111487675B (en) | 2020-03-25 | 2020-03-25 | Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111487675B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112835100A (en) * | 2021-01-21 | 2021-05-25 | 中国石油天然气股份有限公司 | Spectrum decomposition method and device |
| CN114167499A (en) * | 2021-12-02 | 2022-03-11 | 大庆油田有限责任公司 | Method and device for automatically identifying biological reef, electronic equipment and storage medium |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102305940A (en) * | 2011-05-24 | 2012-01-04 | 中国石油集团川庆钻探工程有限公司 | Method for extracting fluid factor |
| CN102466819A (en) * | 2010-11-03 | 2012-05-23 | 中国石油天然气集团公司 | Spectrum analysis method of seismic signal and apparatus thereof |
| US20140297188A1 (en) * | 2013-03-29 | 2014-10-02 | Cgg Services Sa | Time-frequency representations of seismic traces using wigner-ville distributions |
| CN105353408A (en) * | 2015-11-20 | 2016-02-24 | 电子科技大学 | Wigner higher-order spectrum seismic signal spectral decomposition method based on matching pursuit |
| CN105954801A (en) * | 2016-06-12 | 2016-09-21 | 成都理工大学 | Frequency-conversion-attribute-based reservoir permeability evaluation method |
| CN106483555A (en) * | 2016-09-28 | 2017-03-08 | 中国地质大学（武汉） | The Green's function SPWVD Time-Frequency Analysis Method of ENPEMF data |
| CN110865410A (en) * | 2019-11-29 | 2020-03-06 | 山西大学 | Method for suppressing random noise of seismic exploration based on NAR-TFPF |
-
2020
- 2020-03-25 CN CN202010215366.6A patent/CN111487675B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102466819A (en) * | 2010-11-03 | 2012-05-23 | 中国石油天然气集团公司 | Spectrum analysis method of seismic signal and apparatus thereof |
| CN102305940A (en) * | 2011-05-24 | 2012-01-04 | 中国石油集团川庆钻探工程有限公司 | Method for extracting fluid factor |
| US20140297188A1 (en) * | 2013-03-29 | 2014-10-02 | Cgg Services Sa | Time-frequency representations of seismic traces using wigner-ville distributions |
| CN105353408A (en) * | 2015-11-20 | 2016-02-24 | 电子科技大学 | Wigner higher-order spectrum seismic signal spectral decomposition method based on matching pursuit |
| CN105954801A (en) * | 2016-06-12 | 2016-09-21 | 成都理工大学 | Frequency-conversion-attribute-based reservoir permeability evaluation method |
| CN106483555A (en) * | 2016-09-28 | 2017-03-08 | 中国地质大学（武汉） | The Green's function SPWVD Time-Frequency Analysis Method of ENPEMF data |
| CN110865410A (en) * | 2019-11-29 | 2020-03-06 | 山西大学 | Method for suppressing random noise of seismic exploration based on NAR-TFPF |
Non-Patent Citations (2)
| Title |
|---|
| 唐俊: "几种时频分析方法的对比分析", 《价值工程》 * |
| 王雨青: "地震信号分数域属性分析及储层流体识别关键技术研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112835100A (en) * | 2021-01-21 | 2021-05-25 | 中国石油天然气股份有限公司 | Spectrum decomposition method and device |
| CN114167499A (en) * | 2021-12-02 | 2022-03-11 | 大庆油田有限责任公司 | Method and device for automatically identifying biological reef, electronic equipment and storage medium |
| CN114167499B (en) * | 2021-12-02 | 2022-08-26 | 大庆油田有限责任公司 | Method and device for automatically identifying biological reef, electronic equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111487675B (en) | 2021-08-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106597532B (en) | Pre-stack seismic data frequency band expanding method combining well data and horizon data | |
| CN111487675B (en) | Method for generating seismic data high signal-to-noise ratio and high resolution time frequency spectrum | |
| CN109031422A (en) | A kind of seismic signal noise suppressing method based on CEEMDAN and Savitzky-Golay filtering | |
| CN107144879A (en) | A kind of seismic wave noise-reduction method combined based on adaptive-filtering with wavelet transformation | |
| CN109946740B (en) | Seismic resolution enhancement method based on wide flat spectrum seismic wavelet shaping | |
| CN103645507A (en) | A processing method for seismic records | |
| Liu et al. | Seismic quality factor estimation using frequency-dependent linear fitting | |
| CN111474582B (en) | Precise S transformation method for generating high-precision time frequency spectrum | |
| CN112305612B (en) | High-resolution complex spectrum decomposition time-frequency space domain amplitude variation correction method along with offset distance | |
| CN109901224B (en) | Method for protecting and suppressing noise of low-frequency signal of seismic data | |
| CN108957551B (en) | Vibroseis harmonic suppression method based on reconstructed ground force signal | |
| CN113281809B (en) | Spectrum analysis method for seismic signals | |
| CN112255690B (en) | Self-adaptive surrounding rock strong reflection separation method based on seismic phase decomposition | |
| CN116859457A (en) | Space-time-frequency wave number domain synchronous extrusion superposition tight river reservoir characterization method | |
| CN104932008B (en) | The method that the multiple time-frequency spectrum of compensation J conversion improves seismic profile resolution ratio | |
| CN114371505A (en) | Multi-wavelet inversion method and system based on seismic frequency division technology | |
| CN112213775B (en) | Fidelity frequency-boosting method for high-coverage-frequency pre-stack seismic data | |
| Porsani et al. | An adaptive local-slope SVD filtering approach to enhance events on seismic sections | |
| CN114137606A (en) | Stable spectrum simulation deconvolution method | |
| CN112925013A (en) | Seismic data high-resolution processing method based on full-band continuation fidelity | |
| CN113341462A (en) | Zero-phase processing method for marine seismic data | |
| CN113359187B (en) | Wavelet sidelobe elimination method for seismic data | |
| DONG et al. | Fast implementation technique of adaptive Kalman filtering deconvolution via dyadic wavelet transform | |
| King et al. | The role of interpolation in seismic resolution | |
| Zheng et al. | Time-Frequency Domain Deconvolution based on Synchrosqueezing Generalized S Transform |
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 | ||
| GR01 | Patent grant |