CN112068202A - High-precision time-varying wavelet extraction method and system - Google Patents
High-precision time-varying wavelet extraction method and system Download PDFInfo
- Publication number
- CN112068202A CN112068202A CN202010943833.7A CN202010943833A CN112068202A CN 112068202 A CN112068202 A CN 112068202A CN 202010943833 A CN202010943833 A CN 202010943833A CN 112068202 A CN112068202 A CN 112068202A
- Authority
- CN
- China
- Prior art keywords
- time
- frequency
- seismic data
- frequency spectrum
- target
- 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
- 238000000605 extraction Methods 0.000 title claims abstract description 28
- 238000001228 spectrum Methods 0.000 claims abstract description 94
- 230000009466 transformation Effects 0.000 claims abstract description 53
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 230000006870 function Effects 0.000 claims description 18
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 6
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 6
- 238000011551 log transformation method Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000009933 burial Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000011426 transformation method Methods 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
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 provides a high-precision time-varying wavelet extraction method and a system, comprising the following steps: acquiring target seismic data of a target area; the target seismic data carries geological information of an underground space of a target area; performing time-frequency transformation on the target seismic data to obtain time-frequency seismic data; decomposing time-frequency seismic data by utilizing logarithmic transformation to obtain a time-varying sub-wave time frequency spectrum and a reflection coefficient time frequency spectrum; performing Fourier series fitting on the time-varying wavelet time-frequency spectrum to obtain a fitted time-varying wavelet time-frequency spectrum; and obtaining a reflection coefficient based on the time-frequency seismic data and the fitted time-varying wavelet time-frequency spectrum. The method and the device solve the technical problem that in the prior art, the calculation of the reflection coefficient is inaccurate due to low extraction precision of the time-varying wavelet.
Description
Technical Field
The invention relates to the technical field of seismic data processing, in particular to a high-precision time-varying wavelet extraction method and system.
Background
The extraction of the seismic wavelets is one of the important steps in the seismic data processing, and has important significance on the data processing result. In the process of transmitting seismic waves in an underground medium, due to the existence of internal friction and heterogeneity, frequency components of actual seismic waves can change along with the transmission, so that the seismic wavelets show attenuation characteristics in the process of wave field transmission, and seismic signals are non-stationary signals. The estimation of the seismic wavelet has important influence on seismic data deconvolution, directly determines the longitudinal resolution of the seismic data, and has important significance for identifying buried geology with complex structure, small scale and deep burial depth.
The traditional resolution improvement method is mostly based on a seismic trace convolution model, and the seismic wavelet is assumed to be a stable signal and does not change along with the propagation time. In fact, seismic waves propagate in a geometrically dispersive manner in a heterogeneous, viscoelastic subsurface space as non-stationary signals. The existing methods mostly assume that the reflection coefficient of the underground medium is white-noise distribution and does not accord with the actual data characteristics. Therefore, the technical problem of inaccurate calculation of the reflection coefficient caused by low extraction precision of the time-varying wavelet exists in the prior art.
Disclosure of Invention
In view of the above, the present invention provides a method and a system for extracting a time-varying wavelet with high precision, so as to alleviate the technical problem in the prior art that the reflection coefficient is not calculated accurately due to low accuracy of extracting the time-varying wavelet.
In a first aspect, an embodiment of the present invention provides a high-precision time-varying wavelet extraction method, including: acquiring target seismic data of a target area; the target seismic data carries geological information of an underground space of the target area; performing time-frequency transformation on the target seismic data to obtain time-frequency seismic data; decomposing the time-frequency seismic data by utilizing logarithmic transformation to obtain a time-varying sub-wave time frequency spectrum and a reflection coefficient time frequency spectrum; performing Fourier series fitting on the time-varying wavelet time-frequency spectrum to obtain a fitted time-varying wavelet time-frequency spectrum; and obtaining a reflection coefficient based on the time-frequency seismic data and the fitted time-varying wavelet time-frequency spectrum.
Further, the time-frequency transformation is preset generalized S transformation; performing time-frequency transformation on the target seismic data to obtain time-frequency seismic data, wherein the time-frequency seismic data comprises the following steps: by the formula:to the abovePerforming preset generalized S transformation on the target seismic data to obtain time-frequency seismic data; wherein τ is time, f is frequency, x (t) is the target seismic data, λ, b and p are three preset constants, VsAnd x (tau, f) is the time-frequency seismic data.
Further, decomposing the time-frequency seismic data by using logarithmic transformation to obtain a time-varying wavelet time-frequency spectrum and a reflection coefficient time-frequency spectrum, including: decomposing the time-frequency seismic data into the product of time-varying sub-wave time frequency spectrum and reflection coefficient time frequency spectrum; and decomposing the log-form time-frequency seismic data into the sum of log-form time-varying wavelet time-frequency spectrum and log-form reflection coefficient time-frequency spectrum by using log transformation.
Further, performing fourier series fitting on the time-varying wavelet time-frequency spectrum to obtain a fitted time-varying wavelet time-frequency spectrum, including: establishing a target function when Fourier series fitting is carried out on the time-varying sub-wave time frequency spectrum through a least square method; wherein the mathematical expression of the objective function is:f(τ,fm) In order to perform Fourier series expansion,is the square of the norm of L2,for fitting error, m is the frequency point, fmThe frequency of the mth frequency point; and carrying out iterative solution on the objective function through a confidence domain algorithm to obtain a fitted time-varying wavelet time-frequency spectrum.
Further, obtaining a reflection coefficient based on the time-frequency seismic data and the fitted time-varying wavelet time-frequency spectrum, including: obtaining a target reflection coefficient time frequency spectrum based on the time-frequency seismic data, the fitted time-varying wavelet time frequency spectrum and a target relational expression between the time-frequency seismic data and the time-varying wavelet time frequency spectrum; the target relation is: ln | Vsx(τ,f)|≈ln|σ(τ,f)|+ln|Vsr(τ,f)|,Vsx (tau, f) is the time-frequency seismic data, sigma (tau, f) is the time-varying sub-wave time-frequency spectrum, Vsr (τ, f) is the reflection coefficient time spectrum; and performing time-frequency transformation inverse transformation on the target reflection coefficient time-frequency spectrum to obtain a reflection coefficient.
In a second aspect, an embodiment of the present invention further provides a high-precision time-varying wavelet extraction system, including: the system comprises an acquisition module, a transformation module, a decomposition module, a fitting module and a calculation module, wherein the acquisition module is used for acquiring target seismic data of a target area; the target seismic data carries geological information of an underground space of the target area; the transformation module is used for carrying out time-frequency transformation on the target seismic data to obtain time-frequency seismic data; the decomposition module is used for decomposing the time-frequency seismic data by utilizing logarithmic transformation to obtain a time-varying sub-wave time frequency spectrum and a reflection coefficient time frequency spectrum; the fitting module is used for performing Fourier series fitting on the time-varying wavelet time-frequency spectrum to obtain a fitted time-varying wavelet time-frequency spectrum; and the computing module is used for obtaining a reflection coefficient based on the time-frequency seismic data and the fitted time-varying wavelet time-frequency spectrum.
Further, the time-frequency transformation is preset generalized S transformation; the transformation module is further configured to: by the formula:performing preset generalized S transformation on the target seismic data to obtain time-frequency seismic data; wherein τ is time, f is frequency, x (t) is the target seismic data, λ, b and p are three preset constants, VsAnd x (tau, f) is the time-frequency seismic data.
Further, the fitting module is further configured to: establishing a target function when Fourier series fitting is carried out on the time-varying sub-wave time frequency spectrum through a least square method; wherein the mathematical expression of the objective function is:f(τ,fm) In order to perform Fourier series expansion,is the square of the norm of L2,for fitting error, m is the frequency point, fmThe frequency of the mth frequency point; and carrying out iterative solution on the objective function through a confidence domain algorithm to obtain a fitted time-varying wavelet time-frequency spectrum.
In a third aspect, an embodiment of the present invention further provides an electronic 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 steps of the method according to the first aspect when executing the computer program.
In a fourth aspect, the present invention also provides a computer-readable medium having non-volatile program code executable by a processor, where the program code causes the processor to execute the method of the first aspect.
According to the high-precision time-varying wavelet extraction method and system, Fourier series expansion fitting is adopted in the wavelet extraction process, the accuracy of wavelet fitting is improved, high-precision wavelet extraction is achieved, the calculated reflection coefficient is more accurate, and the technical problem that the calculation of the reflection coefficient is inaccurate due to the fact that the extraction precision of the time-varying wavelet is low in the prior art is solved.
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 description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flowchart of a method for extracting high-precision time-varying wavelets according to an embodiment of the present invention;
fig. 2 is a schematic diagram of a high-precision time-varying wavelet extraction system according to an embodiment of the present invention.
Detailed Description
The technical solutions of the present invention will be described clearly and completely with reference to the accompanying drawings, and it should be understood that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The first embodiment is as follows:
FIG. 1 is a diagram illustrating a method for extracting a high-precision time-varying wavelet according to an embodiment of the present invention. As shown in fig. 1, the method specifically includes the following steps:
step S102, target seismic data of a target area are obtained; the target seismic data carries the geological information of the underground space of a target area, and the target area is an area to be processed.
Preferably, in an embodiment of the present invention, the target seismic data is shot gather data.
And step S104, performing time-frequency transformation on the target seismic data to obtain time-frequency seismic data.
Optionally, the time-frequency transform is a preset generalized S-transform.
Specifically, by the formula:performing preset generalized S transformation on target seismic data to obtain time-frequency seismic data; wherein tau is time, f is frequency, x (t) is target seismic data, lambda, b and p are three preset constants, VsAnd x (tau, f) is time-frequency seismic data.
In the embodiment of the invention, the preset generalized S transformation is an improved transformation of the traditional generalized S transformation, the problem that the energy of the seismic signal in an S domain leaks along a time axis is considered, and the preset constant b is an adjustable parameter related to the frequency f.
And S106, decomposing the time-frequency seismic data by utilizing logarithmic transformation to obtain a time-varying sub-wave time frequency spectrum and a reflection coefficient time frequency spectrum.
Optionally, decomposing the time-frequency seismic data into a product of a time-varying wavelet time-frequency spectrum and a reflection coefficient time-frequency spectrum; for example, time-frequency seismic data may be represented as: i Vsx(τ,f)|≈|σ(τ,f)||Vsr(τ,f)|,|Vsx (tau, f) | is time-frequency seismic data obtained by non-stationary seismic channel generalized S transformation, | sigma (tau, f) | is time-varying wavelet time-frequency spectrum, | Vsr (τ, f) | is the reflection coefficient time spectrum obtained by the generalized S transform of the reflection coefficient sequence.
And decomposing the log-form time-frequency seismic data into the sum of log-form time-varying wavelet time-frequency spectrum and log-form reflection coefficient time-frequency spectrum by using log transformation. For example, taking the logarithm of time-frequency seismic data may result in: ln | Vsx(τ,f)|≈ln|σ(τ,f)|+ln|Vsr(τ,f)|。
And S108, performing Fourier series fitting on the time-varying wavelet time-frequency spectrum to obtain the fitted time-varying wavelet time-frequency spectrum.
And step S110, obtaining a reflection coefficient based on the time-frequency seismic data and the fitted time-varying wavelet time-frequency spectrum.
According to the high-precision time-varying wavelet extraction method provided by the invention, Fourier series expansion fitting is adopted in the wavelet extraction process, the accuracy of wavelet fitting is improved, high-precision wavelet extraction is realized, the calculated reflection coefficient is more accurate, and the technical problem of inaccurate calculation of the reflection coefficient caused by low extraction precision of the time-varying wavelet in the prior art is solved.
Optionally, the fitting process in step S108 specifically includes:
step S1081, establishing a target function when Fourier series fitting is carried out on a time-varying sub-wave time frequency spectrum through a least square method; wherein, the mathematical expression of the objective function is:
f(τ,fm) Is a Fourier seriesUnfolding:
is the square of the norm of L2,for fitting error, m is the frequency point, fmIs the frequency of the m-th frequency point, a0As fitting coefficient, anAnd bnIs the nth fitting coefficient, kτIs a positive number and N is the fitting order.
And step S1082, performing iterative solution on the objective function through a confidence domain algorithm to obtain a fitted time-varying wavelet time-frequency spectrum.
Optionally, step S110 further includes the steps of:
step S1101, fitting a time-varying wavelet time-frequency spectrum after fitting based on time-frequency seismic data, and a target relational expression between the time-varying wavelet time-frequency spectrum and the time-varying wavelet seismic data to obtain a target reflection coefficient time-frequency spectrum; the target relationship is: ln | Vsx(τ,f)|≈ln|σ(τ,f)|+ln|Vsr(τ,f)|,Vsx (tau, f) is time-frequency seismic data, sigma (tau, f) is time-varying sub-wave time-frequency spectrum, Vsr (τ, f) is the reflection coefficient time spectrum.
Step S1102, time-frequency transformation inverse transformation is carried out on the target reflection coefficient time-frequency spectrum to obtain a reflection coefficient. The time-frequency transformation is preset generalized S transformation, and the mathematical expression of the time-frequency transformation inverse transformation is as follows:
where s (t) is the reflection coefficient.
The invention provides a high-precision time-varying wavelet extraction method, which comprises the steps of firstly, obtaining target seismic data of a target area; then, performing time-frequency transformation on the target seismic data to obtain time-frequency seismic data; decomposing time-frequency seismic data by utilizing logarithmic transformation to obtain a time-varying sub-wave time frequency spectrum and a reflection coefficient time frequency spectrum; then, Fourier series fitting is carried out on the time-varying wavelet time-frequency spectrum to obtain a fitted time-varying wavelet time-frequency spectrum; and finally, obtaining a reflection coefficient based on the time-frequency seismic data and the fitted time-varying wavelet time-frequency spectrum. The time-frequency transformation method has the advantages that the problem of energy leakage of seismic signals in a time-frequency domain is considered, and the assumption that the reflection coefficient is white is overcome; meanwhile, Fourier series expansion fitting is adopted in the wavelet extraction process, the accuracy of wavelet fitting is improved, high-precision wavelet extraction is achieved, and further, the reflection coefficient is calculated according to the extracted high-precision wavelet, so that the obtained reflection coefficient is more accurate.
Example two:
FIG. 2 is a diagram illustrating a high-precision time-varying wavelet extraction system according to an embodiment of the present invention. As shown in fig. 2, the system includes: an acquisition module 10, a transformation module 20, a decomposition module 30, a fitting module 40 and a calculation module 50.
Specifically, the acquiring module 10 is configured to acquire target seismic data of a target area; the target seismic data carries subsurface space geological information for the target region.
Preferably, in an embodiment of the present invention, the target seismic data is shot gather data.
And the transformation module 20 is configured to perform time-frequency transformation on the target seismic data to obtain time-frequency seismic data.
Optionally, the time-frequency transform is a preset generalized S-transform.
Specifically, the transformation module 20 is further configured to: by the formula:performing preset generalized S transformation on target seismic data to obtain time-frequency seismic data; wherein tau is time, f is frequency, x (t) is target seismic data, lambda, b and p are three preset constants, VsAnd x (tau, f) is time-frequency seismic data.
And the decomposition module 30 is configured to decompose the time-frequency seismic data by using logarithmic transformation to obtain a time-varying sub-wave time-frequency spectrum and a reflection coefficient time-frequency spectrum.
And the fitting module 40 is configured to perform fourier series fitting on the time-varying wavelet time-frequency spectrum to obtain a fitted time-varying wavelet time-frequency spectrum.
And the calculating module 50 is used for obtaining the reflection coefficient based on the time-frequency seismic data and the fitted time-varying wavelet time-frequency spectrum.
According to the high-precision time-varying wavelet extraction system, Fourier series expansion fitting is adopted in the wavelet extraction process, the accuracy of wavelet fitting is improved, high-precision wavelet extraction is achieved, the calculated reflection coefficient is more accurate, and the technical problem that the calculation of the reflection coefficient is inaccurate due to the fact that the extraction accuracy of the time-varying wavelet is low in the prior art is solved.
Optionally, the fitting module 40 is further configured to: establishing a target function when Fourier series fitting is carried out on a time-varying sub-wave time spectrum by a least square method; and (4) carrying out iterative solution on the objective function through a confidence domain algorithm to obtain a fitted time-varying wavelet time-frequency spectrum.
Wherein, the mathematical expression of the objective function is:
f(τ,fm) In order to perform Fourier series expansion,is the square of the norm of L2,for fitting error, m is the frequency point, fmIs the frequency of the mth frequency point.
The embodiment of the present invention further provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and capable of running on the processor, and when the processor executes the computer program, the steps of the method in the first embodiment are implemented.
The embodiment of the invention also provides a computer readable medium with a non-volatile program code executable by a processor, wherein the program code causes the processor to execute the method in the first embodiment.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (10)
1. A high-precision time-varying wavelet extraction method is characterized by comprising the following steps:
acquiring target seismic data of a target area; the target seismic data carries geological information of an underground space of the target area;
performing time-frequency transformation on the target seismic data to obtain time-frequency seismic data;
decomposing the time-frequency seismic data by utilizing logarithmic transformation to obtain a time-varying sub-wave time frequency spectrum and a reflection coefficient time frequency spectrum;
performing Fourier series fitting on the time-varying wavelet time-frequency spectrum to obtain a fitted time-varying wavelet time-frequency spectrum;
and obtaining a reflection coefficient based on the time-frequency seismic data and the fitted time-varying wavelet time-frequency spectrum.
2. The method of claim 1, wherein the time-frequency transform is a preset generalized S-transform; performing time-frequency transformation on the target seismic data to obtain time-frequency seismic data, wherein the time-frequency seismic data comprises the following steps:
by the formula:performing preset generalized S transformation on the target seismic data to obtain time-frequency seismic data;
wherein τ is time, f is frequency, x (t) is the target seismic data, λ, b and p are three preset constants, VsAnd x (tau, f) is the time-frequency seismic data.
3. The method of claim 1, wherein decomposing the time-frequency seismic data using a logarithmic transformation to obtain a time-varying wavelet time-frequency spectrum and a reflection coefficient time-frequency spectrum comprises:
decomposing the time-frequency seismic data into the product of time-varying sub-wave time frequency spectrum and reflection coefficient time frequency spectrum;
and decomposing the log-form time-frequency seismic data into the sum of log-form time-varying wavelet time-frequency spectrum and log-form reflection coefficient time-frequency spectrum by using log transformation.
4. The method of claim 2, wherein fitting a fourier series to the time-varying wavelet time-frequency spectrum to obtain a fitted time-varying wavelet time-frequency spectrum comprises:
establishing a target function when Fourier series fitting is carried out on the time-varying sub-wave time frequency spectrum through a least square method; wherein the mathematical expression of the objective function is: f(τ,fm) In order to perform Fourier series expansion,is the square of the norm of L2,for fitting error, m is the frequency point, fmIs the frequency of the m-th frequency point;
And carrying out iterative solution on the objective function through a confidence domain algorithm to obtain a fitted time-varying wavelet time-frequency spectrum.
5. The method of claim 1, wherein deriving reflection coefficients based on the time-frequency seismic data and the fitted time-varying wavelet time-frequency spectrum comprises:
obtaining a target reflection coefficient time frequency spectrum based on the time-frequency seismic data, the fitted time-varying wavelet time frequency spectrum and a target relational expression between the time-frequency seismic data and the time-varying wavelet time frequency spectrum; the target relation is: ln | Vsx(τ,f)|≈ln|σ(τ,f)|+ln|Vsr(τ,f)|,Vsx (tau, f) is the time-frequency seismic data, sigma (tau, f) is the time-varying sub-wave time-frequency spectrum, Vsr (τ, f) is the reflection coefficient time spectrum;
and performing time-frequency transformation inverse transformation on the target reflection coefficient time-frequency spectrum to obtain a reflection coefficient.
6. A high accuracy time-varying wavelet extraction system, comprising: an acquisition module, a transformation module, a decomposition module, a fitting module and a calculation module, wherein,
the acquisition module is used for acquiring target seismic data of a target area; the target seismic data carries geological information of an underground space of the target area;
the transformation module is used for carrying out time-frequency transformation on the target seismic data to obtain time-frequency seismic data;
the decomposition module is used for decomposing the time-frequency seismic data by utilizing logarithmic transformation to obtain a time-varying sub-wave time frequency spectrum and a reflection coefficient time frequency spectrum;
the fitting module is used for performing Fourier series fitting on the time-varying wavelet time-frequency spectrum to obtain a fitted time-varying wavelet time-frequency spectrum;
and the computing module is used for obtaining a reflection coefficient based on the time-frequency seismic data and the fitted time-varying wavelet time-frequency spectrum.
7. The system according to claim 6, wherein the time-frequency transform is a preset generalized S-transform; the transformation module is further configured to:
by the formula:performing preset generalized S transformation on the target seismic data to obtain time-frequency seismic data;
wherein τ is time, f is frequency, x (t) is the target seismic data, λ, b and p are three preset constants, VsAnd x (tau, f) is the time-frequency seismic data.
8. The system of claim 7, wherein the fitting module is further configured to:
establishing a target function when Fourier series fitting is carried out on the time-varying sub-wave time frequency spectrum through a least square method; wherein the mathematical expression of the objective function is: f(τ,fm) In order to perform Fourier series expansion,is the square of the norm of L2,for fitting error, m is the frequency point, fmThe frequency of the mth frequency point;
and carrying out iterative solution on the objective function through a confidence domain algorithm to obtain a fitted time-varying wavelet time-frequency spectrum.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the steps of the method of any of the preceding claims 1 to 5 are implemented when the computer program is executed by the processor.
10. A computer-readable medium having non-volatile program code executable by a processor, wherein the program code causes the processor to perform the method of any of claims 1-5.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010943833.7A CN112068202B (en) | 2020-09-09 | 2020-09-09 | High-precision time-varying wavelet extraction method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010943833.7A CN112068202B (en) | 2020-09-09 | 2020-09-09 | High-precision time-varying wavelet extraction method and system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112068202A true CN112068202A (en) | 2020-12-11 |
CN112068202B CN112068202B (en) | 2021-08-31 |
Family
ID=73663232
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010943833.7A Active CN112068202B (en) | 2020-09-09 | 2020-09-09 | High-precision time-varying wavelet extraction method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112068202B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102096101A (en) * | 2010-11-24 | 2011-06-15 | 中国石油天然气集团公司 | Method and device for extracting hybrid-phase seismic wavelets |
CN102176054A (en) * | 2011-02-18 | 2011-09-07 | 中国石油化工股份有限公司 | Near-surface comprehensive information processing explanation method |
US20120287755A1 (en) * | 2011-05-13 | 2012-11-15 | Conocophillips Company | Seismic true estimated wavelet |
CN106125139A (en) * | 2016-06-20 | 2016-11-16 | 中国神华能源股份有限公司 | A kind of D seismic modeling method and system |
CN111208561A (en) * | 2020-01-07 | 2020-05-29 | 自然资源部第一海洋研究所 | Seismic acoustic wave impedance inversion method based on time-varying wavelet and curvelet transformation constraint |
-
2020
- 2020-09-09 CN CN202010943833.7A patent/CN112068202B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102096101A (en) * | 2010-11-24 | 2011-06-15 | 中国石油天然气集团公司 | Method and device for extracting hybrid-phase seismic wavelets |
CN102176054A (en) * | 2011-02-18 | 2011-09-07 | 中国石油化工股份有限公司 | Near-surface comprehensive information processing explanation method |
US20120287755A1 (en) * | 2011-05-13 | 2012-11-15 | Conocophillips Company | Seismic true estimated wavelet |
CN106125139A (en) * | 2016-06-20 | 2016-11-16 | 中国神华能源股份有限公司 | A kind of D seismic modeling method and system |
CN111208561A (en) * | 2020-01-07 | 2020-05-29 | 自然资源部第一海洋研究所 | Seismic acoustic wave impedance inversion method based on time-varying wavelet and curvelet transformation constraint |
Non-Patent Citations (1)
Title |
---|
吴昊骏,等: "基于雷管实际延时范围的逐孔爆破振动合成计算与应用", 《爆炸与冲击》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112068202B (en) | 2021-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sun et al. | Cross-correlation analysis and time delay estimation of a homologous micro-seismic signal based on the Hilbert–Huang transform | |
WO2022057305A1 (en) | Signal processing method and apparatus, terminal device and storage medium | |
CN113777650B (en) | Sparse time-frequency spectrum decomposition method, device and equipment based on mixed norm and wavelet transform and storage medium | |
CN113093272A (en) | Time domain full waveform inversion method based on convolutional coding | |
CN109557578A (en) | A kind of reservoir gas-bearing property detection method and device | |
Aghamiry et al. | Interval-Q estimation and compensation: An adaptive dictionary-learning approach | |
CN110794458B (en) | Gas content detection method and device based on time-frequency analysis and storage medium | |
CN111060961B (en) | Quality factor determination method, device and system based on multi-information constraint inversion | |
Qi et al. | Seismic time–frequency spectrum analysis based on local polynomial Fourier transform | |
Fu et al. | Generalized orthogonal matching pursuit with singular value decomposition | |
Yang et al. | A seismic interpolation and denoising method with curvelet transform matching filter | |
CN112068202B (en) | High-precision time-varying wavelet extraction method and system | |
CN111505709B (en) | Attenuation qualitative analysis method based on sparse spectral decomposition | |
CN112904412A (en) | Mine microseismic signal P wave first arrival moment extraction method and system | |
CN115079257A (en) | Q value estimation and seismic attenuation compensation method based on fusion network | |
Chai et al. | Q-compensated acoustic impedance inversion of attenuated seismic data: Numerical and field-data experiments | |
CN111352157B (en) | Shear wave static correction method and system | |
Geetha et al. | Microseismic signal denoising based on variational mode decomposition with adaptive non-local means filtering | |
US20230021093A1 (en) | Electromagnetic wave field data processing method and apparatus, and medium | |
Xu et al. | Ground-penetrating radar time-frequency analysis method based on synchrosqueezing wavelet transformation | |
CN109633744B (en) | Method, device and equipment for extracting seismic wavelets and storage medium | |
CN109459788B (en) | Stratum quality factor calculation method and system | |
Huang et al. | Shannon entropy-based seismic local correlation measure and enhancement | |
Tan et al. | Shape dynamic time warping for seismic waveform inversion | |
CN116736378B (en) | Lq norm constraint-based velocity analysis method for seismic exploration |
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 |