CN104076396A - Variable Q value calculation method based on unsteady state deconvolution - Google Patents
Variable Q value calculation method based on unsteady state deconvolution Download PDFInfo
- Publication number
- CN104076396A CN104076396A CN201410151541.4A CN201410151541A CN104076396A CN 104076396 A CN104076396 A CN 104076396A CN 201410151541 A CN201410151541 A CN 201410151541A CN 104076396 A CN104076396 A CN 104076396A
- Authority
- CN
- China
- Prior art keywords
- value
- time
- unstable state
- window
- die
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 13
- 238000000034 method Methods 0.000 claims abstract description 23
- 238000010521 absorption reaction Methods 0.000 claims description 22
- 230000008859 change Effects 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000011002 quantification Methods 0.000 claims description 5
- 238000007619 statistical method Methods 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims description 2
- 238000003379 elimination reaction Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 6
- 230000002411 adverse Effects 0.000 abstract description 4
- 230000007246 mechanism Effects 0.000 abstract description 2
- 238000001228 spectrum Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 description 5
- 238000012937 correction Methods 0.000 description 5
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 238000013016 damping Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000011551 log transformation method Methods 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005728 strengthening Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000000280 densification Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a variable Q value calculation method based on unsteady state deconvolution. According to the method, energy relation in hyperbolic stripes is fully utilized, counted and analyzed in a logarithm domain, a basic technique idea of Q calculation is established by utilizing linear fitting slope information, adverse effects of seismic wavelet are suppressed constantly with the help of an iteration mechanism, and the longitudinal and transverse stability and practicability of estimation of Q value are improved. According to the method, the Q value is calculated by utilizing an energy statistic relation for the first time, compared with a traditional method, data of a whole time-frequency domain can be utilized fully, information amount is sufficient, estimation results are more reliable, and matching situation with geological background is better. The idea of utilizing iteration linear fitting for suppressing the influences of the seismic wavelet constantly is emphasized, and approximation to precision of the underground Q value is further achieved.
Description
Technical field
The present invention relates to belong to geophysical survey field, relate in particular to a kind of change Q value computation based on unstable state deconvolution.
Background technology
Unstable state feature is an important research aspect in seismic data process field.Existing most of deconvolution are all based upon on the basis of classical Robinson convolution model, but this convolution model is thought not temporal evolution of seismic wavelet approx, the unstable state feature understanding intrinsic to seismic signal is not enough, therefore this class deconvolution is called again static deconvolution, and its practical value is restricted to a great extent.Be different from static deconvolution technology, Gabor deconvolution is taking unstable state convolution model as basis, emphasis has been considered attenuation by absorption impact, the unstable state feature that removes seismic signal by remove the mode of attenuation by absorption and source wavelet (propagating wavelet) simultaneously, is called again unstable state deconvolution with this Gabor deconvolution.In the past decade, a large amount of scholars and researchist improve many technical elements of Gabor deconvolution, strengthen its stability and reliability, window function when content comprises, and time-frequency spectrum is estimated, phase correction and more rational smoothing method.These technological layers have also been consolidated status and the field of Gabor deconvolution in seismic data process field.
Different from traditional static deconvolution, Gabor deconvolution can solve the unstable state feature being caused by attenuation by absorption.According to the most widely used Q theory, attenuation by absorption shows as in time to be increased with frequency, and amplitude is exponential damping
α(t，f)=e
-πtf/Q-iH(πtf/Q) (I)
Wherein,
T, f represents respectively Time And Frequency;
Q represents stratum quality factor;
H is Hilbert conversion.
It is generally acknowledged, attenuation by absorption has minimum phase property, so the imaginary part in above formula provides minimum phase information, H represents Hilbert conversion.In fact, non-stationary seismic trace s (t) is at attenuation by absorption α (t, f) and reflection coefficient r (t) dynamically on the basis of convolution, then with source wavelet
carry out the result of static convolution,
Be equivalent to
By means of Gabor conversion, above-mentioned unstable state convolution model can be decomposed into three amassing, as follows
Wherein,
with
represent the Gabor time-frequency spectrum of seismic trace and reflection coefficient.
The Gabor that the Gabor conversion of formula (4) expression unstable state seismic trace is approximately equal to Fourier transform, attenuation by absorption and the reflection coefficient of source wavelet converts this three products.If only consider magnitude relation, have
Wherein, front two ripples of varitron while making to describe become possibility, are collectively referred to as the amplitude spectrum of making to propagate wavelet
in brief, formula (5) has been described the unstable state convolution model in Gabor transform domain, has formed the basis of Gabor deconvolution (being unstable state deconvolution).Different from classical inverse convolution, unstable state deconvolution is by analyzing and the energy relationship of balanced each hyperbolic band is realized attenuation by absorption and compensated, and its outstanding feature is not need to estimate that Q has just realized the correction to unstable state feature.In fact, this feature has also limited academia and has utilized the thinking of unstable state deconvolution estimation Q value to a certain extent.
Summary of the invention
The technical problem to be solved in the present invention is that current unstable state deconvolution research is trapped in raising resolution aspect, and asks calculation aspect never to relate to excavating its profound intension-Q value.In addition, current Q values is asked the information wretched insufficiency of calculation technology to time-frequency domain, causes and asks calculation result can depart from reality.
In order to solve the problems of the technologies described above, a kind of change Q value computation based on unstable state deconvolution of the present invention is, in log-domain, make full use of and statistical study hyperbolic band in energy relationship, ask the basic fundamental thinking of calculation by utilizing the slope information of linear fit to set up Q, thereby and strengthen by the adverse effect of the continuous compacting of iteration mechanism seismic wavelet vertical lateral stability and the practicality that Q value is estimated.
The concrete steps of this method are as follows:
Step 1: take the logarithm on the basis of the seismic trace that does not carry out attenuation compensation being carried out to Gabor conversion, enter logarithm time-frequency domain;
Step 2: at logarithm time-frequency domain along time orientation by the true timing window of actual conditions, time window divide hyperbolic band and add up its energy relationship, form die-away curve;
Step 3: the major part that die-away curve is carried out is carried out linear fit, determines the Q value of this analysis window center according to fit slope;
Step 4: successively estimate Q size by seismic trace timesharing window;
Step 5: obtaining on the basis of Q achievement, rejecting wild value, interpolation, level and smooth and carry out histogram adjusting according to Given information, forming Q value and estimate result.
In step 2, in the time of logarithm Time-Frequency Analysis, in window, the die-away curve being formed by each hyperbolic band is realized being caused that by attenuation by absorption the quantification of unstable state feature represents, the comprehensive embodiment of Q value information in window when in fact its linear relationship has embodied.
In step 3, utilize best linear fit mode to realize the estimation to Q value, the method itself is a kind of statistical analysis means, has utilized more Time-Frequency Information, and the Q value estimation result of acquisition is more stable.
Change Q value computation based on unstable state deconvolution of the present invention compared with prior art has following beneficial effect.
1, realized first and utilized energy statistics relation to ask calculation Q value, compared with classic method, can utilize more fully the data of whole time-frequency domain, quantity of information abundance, estimation result is more reliable, with the identical better off of geologic background.
2, emphasize to utilize the thinking of iteration linear fit constantly to suppress seismic wavelet impact, and then realize the precision of underground Q value is approached.
The ultimate principle of the change Q value computation based on unstable state deconvolution provided by the invention is as follows:
At log-domain, unstable state convolution model is converted to linear equation
Wherein,
with
represent the log-transformation of the Gabor time-frequency spectrum of seismic trace, source wavelet and reflection coefficient.
Slope and the quality factor q of this linear equation are closely related, and slope reacts the combined influence of source wavelet and reflection coefficient.At logarithm time-frequency domain, attenuation by absorption isoline still present hyperbolic form, difference is to transfer linear relationship to by exponential relationship between isoline.Similarly, the logarithm time-frequency spectrum of seismic trace carried out to the division of hyperbolic band and calculate each band average to form die-away curve.Say in essence, this die-away curve has been realized being represented by a kind of quantification of attenuation by absorption unstable state feature.Slope characteristics on this die-away curve is directly associated with Q value, according to the first term of formula (6), does not exist in the stratum of attenuation by absorption in perfect elasticity, and this die-away curve should present the straight line of basic horizontal, represents that stratum Q value is infinitely great; Being formed by compacted rock in weak attenuation by absorption stratum, this die-away curve there is less slope characteristics, represent that stratum Q value is relatively high; And at the near Faults being caused by unconsolidated formation or tectonic movement, the slope of this die-away curve is larger, represent that stratum Q value is relatively little.In sum, the slope of die-away curve has direct indicative function to the size of stratum quality factor q in essence, is the important foundation that the present invention solves Q value.
But be subject to the impact of seismic wavelet, and there is relatively large deviation in initial decay slope of a curve and actual value, and the present invention uses the iterative analysis in logarithm time-frequency domain to reduce wavelet impact, the true form of constantly approaching die-away curve for this reason.On the basis that obtains die-away curve, its major part is carried out to linear fit, the process of this linear fit is close to a least square inverting.Mathematical notation is
A'AX=A'b (7)
Wherein,
A represents parameter matrix,
X represents parameter vector to be solved.
For strengthening stability, the adverse effect of the disturbance that can produce in factors such as equation left end increase damping term reduction noises to estimation slope,
(A'A+μΙ)X=A'b (8)
In formula,
μ represents the prewhitening factor,
Ι represents unit matrix (matrix size is 1) herein.
Utilize the slope information X producing to ask calculation to Q value, obtain
Realize the Q value estimation in analysis window with the method.According to above-mentioned thinking, to seismic section press time-frequency figure along time orientation according to debugging and determining a series of analysis windows, successively add up die-away curve, linear fit, estimation Q value, and to its carry out smoothly, the whole Q value section of processing acquisition such as interpolation and correction.
Brief description of the drawings
Below in conjunction with the drawings and specific embodiments, the change Q value computation based on unstable state deconvolution of the present invention is described in further detail.
The hyperbolic of Fig. 1 (a)-Fig. 1 (b) logarithm time-frequency domain is divided and to slope iterative approach process;
In Fig. 2 (a)-Fig. 2 (e) actual seismic road, utilize time frequency distribution map and the die-away curve form of unstable state Deconvolution Correction front and back;
Fig. 3 (a)-Fig. 3 (b) does not make the actual seismic section of attenuation by absorption compensation and the Q value Profile Correlation that utilizes the present invention to solve;
Western certain work area of Fig. 4 (a)-Fig. 4 (b) China is contrasted by the coherence slice of conventional pre-stack time migration imaging results and the coherence slice based on the present invention calculates the prestack Q compensating offset imaging results of Q value.
Fig. 5 is the change Q value computation process flow diagram based on unstable state deconvolution of the present invention.
Embodiment
As shown in Figure 1.Obtain the process of true slope feature with the division of numerical model test specification hyperbolic band and the impact of iteration compacting wavelet at this.In this model, have three layers of reflectivity model, this model is divided into three layers, and the reflection coefficient amplitude level in middle layer is upper and lower two-layer 1/3.In model measurement, source wavelet and this reflection coefficient sequence synthetic seismogram of minimum phase that to use dominant frequency be 40hz, given Q=32 decays to this seismic trace at time-frequency domain also subsequently, obtains thus final theogram.In Fig. 1 a, signal has shown 20 hyperbolic bands divisions (being the isoline of time and frequency product).If the data mean value feature of adding up for the first time under each hyperbolic band is set up die-away curve as shown in blue curve in Fig. 1 b, and affected by source wavelet, first linear fit result and true slope (black line in figure) certainly exist obvious errors.For this reason, adopt iteration thinking constantly estimate source wavelet and remove, progressively suppress its impact and realize constantly approaching true slope with this.After three iteration, can find out that fit slope is in close proximity to true slope.This example has confirmed that taking numerical model as basis iteration proposes the feasibility of slop estimation.
As shown in Figure 2, the die-away curve in actual seismic road also has linear feature conventionally.Fig. 2 a has shown the seismic section in certain marine site, and extracting wherein its corresponding time frequency distribution map is together 2b, and its die-away curve (see that Fig. 2 c) has obvious linear feature, caused by filtering in early stage by smooth top and bottom.Utilize unstable state deconvolution technology, the frequency difference (as shown in Figure 2 d) of this seismic trace top of balance and bottom gradually, smooths the die-away curve feature of this inclination that (Fig. 2 e) substantially.In fact this technical indicator that, unstable state deconvolution technology contains in realizing proposes initial inspiration of the present invention and power just.(to the explanation of Fig. 2 c?)
As shown in Figure 3, Fig. 3 a is a seismic section of the western Carbonate Reservoir of China, and the upper left quarter of this seismic section is fault development district, affected by the factors such as rock failure, and this district should have stronger attenuation by absorption effect, analyzes in theory this district and should present compared with low reactance-resistance ratio; The bottom of this seismic section is the Cambrian system rock of compacting densification, and its attenuation by absorption effect is very little, should have high Q value.Fig. 3 b becomes Q value section by utilizing the present invention to process to be obtained to this section.Compare with Basic Geological background, Q value Changing Pattern and the above-mentioned geologic background of Fig. 3 b match, the validity of proved invention in actual Q value estimation.
As shown in Figure 4, Fig. 4 a be based on traditional pre-stack time migration imaging results target reservoir position, deep along layer coherence slice.For improving deep reservoir resolution, first utilize the present invention to ask for and become Q value information, and be applied to the imaging of Q compensating offset, corresponding target reservoir position along layer coherence slice as shown in Figure 4 b.Contrast two kinds of results, be not difficult to find that the latter has significant resolution, can effectively identify pinching line (shown in red arrow) and the river feature (shown in green arrow) in deep.
As shown in Figure 5, a kind of change Q value computation based on unstable state deconvolution of the present invention is, in log-domain, analyze the energy statistics relation in hyperbolic band, set up linear fit equation, realize the calculation of asking to Q value from the slope information obtaining, and utilize the means such as elimination of burst noise, interpolation, level and smooth and histogram adjusting to strengthen vertical lateral stability and the practicality that Q value is estimated.
The concrete steps of this method are as follows:
Step 1: take the logarithm on the basis of the seismic trace that does not carry out attenuation compensation being carried out to Gabor conversion, enter logarithm time-frequency domain;
Step 2: at logarithm time-frequency domain along time orientation by the true timing window of actual conditions, time window divide hyperbolic band and add up its energy relationship, form die-away curve;
Step 3: the major part that die-away curve is carried out is carried out linear fit, determines the Q value of this analysis window center according to fit slope;
Step 4: successively estimate Q size by seismic trace timesharing window;
Step 5: obtaining on the basis of Q achievement, rejecting wild value, interpolation, level and smooth and carry out histogram adjusting according to Given information, forming Q value and estimate result.
In step 2, in the time of logarithm Time-Frequency Analysis, in window, the die-away curve being formed by each hyperbolic band is realized being caused that by attenuation by absorption the quantification of unstable state feature represents, the comprehensive embodiment of Q value information in window when in fact its linear relationship has embodied.
In step 3, utilize best linear fit mode to realize the estimation to Q value, the method itself is a kind of statistical analysis means, has utilized more Time-Frequency Information, and the Q value estimation result of acquisition is more stable.
Change Q value computation based on unstable state deconvolution of the present invention compared with prior art has following beneficial effect.
1, realized first and utilized energy statistics to be related to ball calculation Q value, compared with classic method, can utilize more fully the data of whole time-frequency domain, quantity of information abundance, estimation result is more reliable, with the identical better off of geologic background;
2, emphasize to utilize the thinking of iteration linear fit constantly to suppress seismic wavelet impact, and then realize the precision of underground Q value is approached.
The ultimate principle of the change Q value computation based on unstable state deconvolution provided by the invention is as follows:
At log-domain, unstable state convolution model is converted to linear equation
Wherein,
with
represent the log-transformation of the Gabor time-frequency spectrum of seismic trace, source wavelet and reflection coefficient.
Slope and the quality factor q of this linear equation are closely related, and slope reacts the combined influence of source wavelet and reflection coefficient.At logarithm time-frequency domain, attenuation by absorption isoline still present hyperbolic form, difference is to transfer linear relationship to by exponential relationship between isoline.Similarly, the logarithm time-frequency spectrum of seismic trace carried out to the division of hyperbolic band and calculate each band average to form die-away curve.Say in essence, this die-away curve has been realized being represented by a kind of quantification of attenuation by absorption unstable state feature.Slope characteristics on this die-away curve is directly associated with Q value, according to the first term of formula (6), does not exist in the stratum of attenuation by absorption in perfect elasticity, and this die-away curve should present the straight line of basic horizontal, represents that stratum Q value is infinitely great; Being formed by compacted rock in weak attenuation by absorption stratum, this die-away curve there is less slope characteristics, represent that stratum Q value is relatively high; And at the near Faults being caused by unconsolidated formation or tectonic movement, the slope of this die-away curve is larger, represent that stratum Q value is relatively little.In sum, the slope of die-away curve has direct indicative function to the size of stratum quality factor q in essence, is the important foundation that the present invention solves Q value.
But be subject to the impact of seismic wavelet, and there is relatively large deviation in initial decay slope of a curve and actual value, and the present invention uses the iterative analysis in logarithm time-frequency domain to reduce wavelet impact, the true form of constantly approaching die-away curve for this reason.On the basis that obtains die-away curve, its major part is carried out to linear fit, the process of this linear fit is close to a least square inverting.Mathematical notation is
A'AX=A'b (7)
Wherein,
A represents parameter matrix,
X represents parameter vector to be solved.
For strengthening stability, the adverse effect of the disturbance that can produce in factors such as equation left end increase damping term reduction noises to estimation slope,
(A'A+μΙ)X=A'b (8)
In formula,
μ represents the prewhitening factor,
Ι represents unit matrix (matrix size is 1) herein.
Utilize the slope information X producing to ask calculation to Q value, obtain
Realize the Q value estimation in analysis window with the method.According to above-mentioned thinking, to seismic section press time-frequency figure along time orientation according to debugging and determining a series of analysis windows, successively add up die-away curve, linear fit, estimation Q value, and to its carry out smoothly, the whole Q value section of processing acquisition such as interpolation and correction.
Above embodiment is to illustrate the invention and not to limit the present invention.
List of references
Margrave?G?F.Theory?of?nonstationary?linear?filtering?in?Fourier?domain?with?application?to?time-variant?filtering:Geophysics,1998,63,244-259.
Margrave,G.F.,M.P.Lamoureux,and?D.C.Henley,2011,Gabor?deconvolution:Estimating?reflectivity?by?nonstationary?deconvolution?of?seismic?data:Geophysics,76,W15-W30.
Margrave,G.F.,M.P.Lamoureux,D.C.Henley,2011,Gabor?deconvolution:Estimating?reflectivity?by?nonstationary?deconvolution?of?seismic?data:Geophysics,76,W15-W30.
Sacchi,M.D.,1997,Reweighting?strategies?in?seismic?deconvolution:Geophys.J.Int.,129,651-656.
Sun X.K., Sun Z.D., and Zhou X.Y., etc., 2013, Gabor deconvolution based on hyperbolic smoothing in log spectra:75th EAGE Conference & Exhibition Extended Abstract.
Sun X.K., Sun Z.D., Peng T., waits 2012, Gabor deconvolution:hyperbolic smoothing in logarithmic magnitude spectrum:82nd SEG annual meeting Expanded Abstract.
Sun,S.Z.,Y.Bai,S.Wu,and?D.Wang,2012,Two?promising?approaches?for?amplitude-preserved?resolution?enhancement,The?Leading?Edge,31,206-210.
Liu likes the pretty auspicious Zhang Gailan in Wuning, and 2009, Cauchy sparse constraint Bayesian estimates that earthquake blind deconvolution framework and algorithm improve: petroleum exploration, 48 (5), 459-464
Yuan Hemin, 2011, multiple dimensioned joint inversion and application thereof, China University Of Petroleum Beijing's Master's thesis.
Zhang Fanchang, Liu Jie, prints Xing Yao etc., and 2008, revise Cauchy and retrain earthquake blind deconvolution method: geophysical prospecting for oil 2008,43 (4), 391-396.
Claims (4)
1. the change Q value computation based on unstable state deconvolution, it is characterized in that: in log-domain, analyze the energy statistics relation in hyperbolic band, set up linear fit equation, realize the calculation of asking to Q value from the slope information obtaining, and utilize the means such as elimination of burst noise, interpolation, level and smooth and histogram adjusting to strengthen vertical lateral stability and the practicality that Q value is estimated.
2. the change Q value computation based on unstable state deconvolution according to claim 1, is characterized in that; The method step is as follows:
Step 1: take the logarithm on the basis of the seismic trace that does not carry out attenuation compensation being carried out to Gabor conversion, enter logarithm time-frequency domain;
Step 2: at logarithm time-frequency domain along time orientation by the true timing window of actual conditions, time window divide hyperbolic band and add up its energy relationship, form die-away curve;
Step 3: the major part that die-away curve is carried out is carried out linear fit, determines the Q value of this analysis window center according to fit slope;
Step 4: successively estimate Q size by seismic trace timesharing window;
Step 5: obtaining on the basis of Q achievement, rejecting wild value, interpolation, level and smooth and carry out histogram adjusting according to Given information, forming Q value and estimate result.
3. the change Q value computation based on unstable state deconvolution according to claim 2, it is characterized in that: in step 2, in the time of logarithm Time-Frequency Analysis in window, the die-away curve being formed by each hyperbolic band is realized being caused that by attenuation by absorption the quantification of unstable state feature represents, the comprehensive embodiment of Q value information in window when in fact its linear relationship has embodied.
4. the change Q value computation based on unstable state deconvolution according to claim 2, it is characterized in that, in step 3, utilize best linear fit mode to realize the estimation to Q value, the method itself is a kind of statistical analysis means, utilized more Time-Frequency Information, the Q value estimation result of acquisition is more stable.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410151541.4A CN104076396B (en) | 2014-04-16 | 2014-04-16 | Based on the change Q value computation of unstable state deconvolution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410151541.4A CN104076396B (en) | 2014-04-16 | 2014-04-16 | Based on the change Q value computation of unstable state deconvolution |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104076396A true CN104076396A (en) | 2014-10-01 |
| CN104076396B CN104076396B (en) | 2016-01-27 |
Family
ID=51597778
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410151541.4A Expired - Fee Related CN104076396B (en) | 2014-04-16 | 2014-04-16 | Based on the change Q value computation of unstable state deconvolution |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104076396B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110988990A (en) * | 2019-12-10 | 2020-04-10 | 北京化工大学 | High-precision seismic attribute inversion method |
| CN111474586A (en) * | 2020-05-06 | 2020-07-31 | 中南大学 | Frequency domain multi-scale crack weakness inversion method |
| CN112415587A (en) * | 2019-08-21 | 2021-02-26 | 中国石油化工股份有限公司 | Reservoir seismic wave attenuation characteristic analysis method and reservoir reflection coefficient inversion method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050086005A1 (en) * | 2003-10-16 | 2005-04-21 | Taner M. T. | Method for determining formation quality factor from seismic data |
| CN102305941A (en) * | 2011-05-25 | 2012-01-04 | 东北石油大学 | Method for determining stratum stack quality factor by direct scanning of prestack time migration |
| CN102854530A (en) * | 2012-07-13 | 2013-01-02 | 孙赞东 | Hyperbolic smooth dynamic deconvolution method based on logarithm time-frequency domain |
-
2014
- 2014-04-16 CN CN201410151541.4A patent/CN104076396B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050086005A1 (en) * | 2003-10-16 | 2005-04-21 | Taner M. T. | Method for determining formation quality factor from seismic data |
| CN102305941A (en) * | 2011-05-25 | 2012-01-04 | 东北石油大学 | Method for determining stratum stack quality factor by direct scanning of prestack time migration |
| CN102854530A (en) * | 2012-07-13 | 2013-01-02 | 孙赞东 | Hyperbolic smooth dynamic deconvolution method based on logarithm time-frequency domain |
Non-Patent Citations (1)
| Title |
|---|
| 刘局新 等: "基于Gabor变换的地震资料Q值计算方法及应用", 《中国地球物理2010——中国地球物理学会第二十六届年会、中国地震学会第十三次学术大会论文集》, 31 December 2010 (2010-12-31), pages 557 - 558 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112415587A (en) * | 2019-08-21 | 2021-02-26 | 中国石油化工股份有限公司 | Reservoir seismic wave attenuation characteristic analysis method and reservoir reflection coefficient inversion method |
| CN110988990A (en) * | 2019-12-10 | 2020-04-10 | 北京化工大学 | High-precision seismic attribute inversion method |
| CN111474586A (en) * | 2020-05-06 | 2020-07-31 | 中南大学 | Frequency domain multi-scale crack weakness inversion method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104076396B (en) | 2016-01-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Enhancing resolution of nonstationary seismic data by molecular-Gabor transform | |
| CN104849756B (en) | A kind of seismic data resolution that improves strengthens the method for effective weak signal energy | |
| US20140372043A1 (en) | Full Waveform Inversion Using Perfectly Reflectionless Subgridding | |
| CN103645507B (en) | The disposal route of seismologic record | |
| CN107272062B (en) | A kind of Q estimation methods of underground medium of data-driven | |
| CN108614295B (en) | Stratum Q value calculation method based on generalized seismic wavelets | |
| CN104122581B (en) | A kind of poststack sound impedance inversion method | |
| CN102692647B (en) | Stratum oil-gas possibility prediction method with high time resolution | |
| Wang et al. | Seismic data denoising for complex structure using BM3D and local similarity | |
| Li et al. | Seismic acoustic impedance inversion with multi-parameter regularization | |
| CN105089652A (en) | Pseudo-acoustic curve rebuilding and sparse pulse joint inversion method | |
| Zoukaneri et al. | A combined Wigner-Ville and maximum entropy method for high-resolution time-frequency analysis of seismic data | |
| CN103412333A (en) | Static correction base level determining method | |
| Xin et al. | Robust Q tomographic inversion through adaptive extraction of spectral features | |
| CN104391324A (en) | Seismic trace set dynamic correction stretching correction pre-processing technology before AVO inversion depending on frequency | |
| CN104076396A (en) | Variable Q value calculation method based on unsteady state deconvolution | |
| CN104597485B (en) | A kind of miniature fault detection method and fault detection device | |
| CN102854530B (en) | Hyperbolic smooth dynamic deconvolution method based on logarithm time-frequency domain | |
| Zhang et al. | Combining CEEMD and recursive least square for the extraction of time-varying seismic wavelets | |
| Li et al. | A generalized seismic attenuation compensation operator optimized by 2-D mathematical morphology filtering | |
| Liu et al. | Improving the resolution of seismic data based on S-transform and modified variational mode decomposition, an application to Songliao Basin, Northeast China | |
| Cheng et al. | Q estimation based on the logarithmic spectral area double difference | |
| Zhao et al. | Automatic events extraction in pre-stack seismic data based on edge detection in slant-stacked peak amplitude profiles | |
| Jiang et al. | Source-Independent Full-Waveform Inversion Based on Convolutional Wasserstein Distance Objective Function | |
| Zou et al. | Low-frequency impedance modeling based on well logs and seismic attributes |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| ASS | Succession or assignment of patent right |
Owner name: SUN XUEKAI Free format text: FORMER OWNER: SUN ZANDONG Effective date: 20150520 Owner name: SUN ZANDONG Effective date: 20150520 |
|
| C41 | Transfer of patent application or patent right or utility model | ||
| C53 | Correction of patent of invention or patent application | ||
| CB03 | Change of inventor or designer information |
Inventor after: Sun Xuekai Inventor after: Sun Zandong Inventor after: Han Jianfa Inventor after: Peng Tao Inventor before: Sun Xuekai Inventor before: Sun Zandong Inventor before: Han Jianfa |
|
| COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: SUN XUEKAI SUN ZANDONG HAN JIANFA TO: SUN XUEKAI SUN ZANDONG HAN JIANFA PENG TAO |
|
| TA01 | Transfer of patent application right |
Effective date of registration: 20150520 Address after: 102249, China University of Petroleum, 18, Xuefu Road, Beijing, Changping District (Beijing) Applicant after: Sun Xuekai Applicant after: Sun Zandong Address before: 100000 Beijing city Changping District Road No. 18, China University of Petroleum (Beijing) Applicant before: Sun Zandong |
|
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190315 Address after: 102249 18 Fu Xue Road, Changping District, Beijing Patentee after: China University of Petroleum (Beijing) Address before: 102249 Xuefu Road, Changping District, Beijing, China University of Petroleum (Beijing) Co-patentee before: Sun Zandong Patentee before: Sun Xuekai |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160127 Termination date: 20210416 |