CN104570071A - Bayesian linear AVA and AVF retrieval method of sticky sound medium - Google Patents
Bayesian linear AVA and AVF retrieval method of sticky sound medium Download PDFInfo
- Publication number
- CN104570071A CN104570071A CN201310476278.1A CN201310476278A CN104570071A CN 104570071 A CN104570071 A CN 104570071A CN 201310476278 A CN201310476278 A CN 201310476278A CN 104570071 A CN104570071 A CN 104570071A
- Authority
- CN
- China
- Prior art keywords
- obs
- sigma
- omega
- fix
- theta
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a Bayesian linear AVA and AVF retrieval method of a sticky sound medium and belongs to the field of geophysical prospecting for petroleum. The method comprises the following steps: S1, inputting angle-domain seismic data, frequency-domain seismic data and log data; S2, establishing prior information and a likelihood function of the log data; S3, carrying out Bayesian linear AVA and AVF retrieval by using the prior information and the likelihood function so as to obtain a velocity and a quality factor.
Description
Technical field
The invention belongs to field of petroleum geophysical exploration, be specifically related to a kind of sticky linear AVA and AVF inversion method of acoustic medium Bayes, utilize bayesian theory to carry out AVA and AVF inverting based on sticky acoustic medium and obtain acoustic velocity and decay factor.
Background technology
Traditional AVO (amplitude offset distance) or AVA (Amplitudeversusangle) is always based on the perfectly elastic hypothesis of underground medium.And VSP record, Petrophysical measurement in log data and laboratory all explicitly seismic wave in communication process, there will be decay and dispersion phenomenon, especially for the stratum containing hydrocarbon, the decay of seismic event is more obvious.In fact, the impact ignoring seismic wave attenuation brings huge risk by analyzing AVO.
Base area earthquake wave propagation theory, the reflection coefficient that it absorbs by force owing to underground medium by external some scholars.White (1975), De Hoop (1991), Ursin (2002) etc. are about how calculating the discussion absorbing and to give in reflection coefficient in series of theories, general absorption inverse Problem is reduced to the reflection coefficient of the single interface of inverting by Innanen and Weglein (2007), and Innanen (2011) has derived the expression formula of seismic event reflection and transmission coefficient when inciding absorbing medium by the non-attenuation medium of elasticity in detail.In fact, the reflection coefficient absorbed in complete description or inverting underground medium, had both needed the change (AVA) considering amplitude angle, also will consider the variation relation (AVF) of amplitude frequency simultaneously.Although established some mechanism to explain the variation relation of reflection coefficient with frequency at present, how to solve indirect problem and remained an international difficult problem.
Domestic scholars notices the amplitude anomaly in gas-bearing reservoir regional earthquake low frequency part very early, and has carried out correlative study based on the method such as generalized S-transform or wavelet transformation, but not yet someone carries out amplitude frequency change always) systematic study.Appoint (2009) to carry out the mutation analysis research of pore media amplitude frequency in University of Houston, and derived in detail vertical incidence time reflection coefficient with the change of frequency.With reference to AVO classification, appoint (2009) to be divided three classes equally by AVF and distinguish corresponding different gas-bearing reservoirs.This research work is subject to the attention of more domestic scholars very soon, derived further in the basis that Liu (2011) etc. are in office inclined seismic wave time amplitude frequency variation relation, and achieve some useful understanding.Nonetheless, the research about sticky acoustic medium AVA/AVF inverting work still still belongs to blank at home.
Summary of the invention
The object of the invention is to solve the difficult problem existed in above-mentioned prior art, a kind of sticky linear AVA and AVF inversion method of acoustic medium Bayes is provided, geophysics fine description for reservoir is target, the sticky linear AVA/AVF inverting of acoustic medium Bayes is utilized to obtain speed and the dampening information of underground medium, improve reliability and the accuracy of reservoir oil and gas prediction, there is face the practice, bright characteristics that application is strong.
The present invention is achieved by the following technical solutions:
A kind of sticky linear AVA and AVF inversion method of acoustic medium Bayes, comprising:
The first step, input angle territory geological data, frequency field geological data and log data
Second step, sets up prior imformation and the likelihood function of log data;
3rd step, utilizes described prior imformation and likelihood function, carries out linear AVA and the AVF inverting of Bayes, obtains speed and quality factor.
In the described first step, the geological data of angle domain geological data and frequency field is all expressed as d
obs, these data are the 3-D data volume of time, angle and frequency; Log data is expressed as: m=[lnc, α ,]
t, wherein, m is the vector of elastic parameter composition, l
ncrepresent the logarithm of speed respectively, α represents decay factor, and the two all obtains from well-log information.
Described second step is achieved in that
The prior imformation setting up log data is as follows:
Wherein,
represent n the sampling point Normal Distribution about m, μ
m, ∑
mrepresent the expectation and variance of m respectively, by what obtain the statistics of well logging sampling point;
Set up likelihood function as follows:
Wherein,
represent that observation geological data meets Gaussian distribution,
represent the expectation of geological data,
represent the variance of geological data;
Expectation and variance meets respectively:
∑ "
mrepresenting asks second order to lead to the statistical variance of elastic parameter, ∑
erepresent the variance of noise;
Wherein, S is:
Can be write as each S:
For AVA inverting, adopt:
For AVF inverting, adopt:
Wherein, s
1(θ
i, ω
i) expression incident angle is θ
i, dominant frequency is ω
iwavelet; Subscript 1,2 until nwave represents the time delay of wavelet; A
c(θ
n θ, ω
fix) represent frequencies omega
fixfixing, depend on incidence angle θ
n θvelocity coefficient; A
α(θ
n θ, ω
fix) represent frequencies omega
fixfixing, depend on incidence angle θ
n θdecay factor coefficient; A
c(θ
fix, ω
n ω) represent incidence angle θ
fixfixing, depend on frequencies omega
n ωvelocity coefficient; A
α(θ
fix, ω
n ω) represent incidence angle θ
fixfixing, depend on frequencies omega
n ωdecay factor coefficient.
Described 3rd step is achieved in that
Seismologic record and elastic parameter meet joint distribution:
Wherein
At given d
obswhen m Posterior probability distribution be:
Posterior error and variance are respectively:
The maximum a posteriori solution (MAP) of m just equals posterior error value
covariance matrix is evaluated in order to the uncertainty final to inversion result.
Formula (23) is expectation and the variance of m, and namely that m expression is velocities Vp and quality factor q p, and namely the expectation therefore obtaining m obtains the expectation of velocities Vp and quality factor q p, namely maximum probability solution; Variance is then in order to assess the uncertainty of these two parameters.
Compared with prior art, the invention has the beneficial effects as follows: relative to existing technology, invention introduces geological data frequency-domain information, under the guidance of bayesian theory, build the likelihood function of prior imformation needed for inverting and frequency dependence, inversion result not only can obtain the speed of underground medium but also can obtain the decay factor change of medium; The result provided with the form of probability is that follow-up uncertainty analysis provides possibility, also evaluates more comprehensively for reservoir prediction provides simultaneously.
Accompanying drawing explanation
Fig. 1 gas-bearing sandstone reservoir model schematic.
Interface, Fig. 2 reservoir top accurate reflection coefficient is with the change of incident angle and frequency.
Interface, Fig. 3 reservoir top is similar to the change of reflection coefficient with incident angle and frequency.
The interface reflection coefficients error analysis of Fig. 4 reservoir top.
Fig. 5-1 simulated seismogram is with the change (angle-frequency field reflection record data volume) of incident angle and frequency.
Fig. 5-2 simulated seismogram is with the change (being 10 ° of angles incidence sections and 20hz frequency slice) of incident angle and frequency.
The sticky acoustic medium Bayes AVA inversion result of Fig. 6-1.
The sticky acoustic medium Bayes AVF inversion result of Fig. 6-2.
Fig. 7 is the step block diagram of the inventive method.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail:
The approximate expression of sticky acoustic medium AVA and AVF and bayesian theory combine by the present invention, not only utilize the information that amplitude incident angle converts, and take full advantage of the information converting of amplitude frequency, carry out the linear AVA/AVF inverting of Bayes by the expectation under structure probability meaning and covariance matrix, obtain speed and the decay factor of zone of interest.
The inventive method specifically comprises:
(1) derivation and the application of acoustic medium AVA/AVF approximate expression is glued
For sticky acoustic medium, phase velocity formula can be labeled as:
H(ω)=1+Q
-1F(ω),
Wherein c (ω) is phase velocity, and c is the speed of reference frequency.For normal density medium, by the ess-strain condition of continuity of reflecting interface, reflection coefficient can be written as:
Linearization is carried out to acoustic medium reflection coefficient:
?
place carries out Taylor expansion, and retains single order item (weak velocity contrast is similar to):
Phase velocity formula is substituted into, obtains:
?
place carries out Talor expansion, and retains single order item:
Wherein, c is the speed of reference frequency, and θ, ω are respectively angle and frequency, and H is the intermediate variable with frequency dependence, specifically sees formula (1).
Formula (6) is sticky acoustic medium Bayes AVA/AVF inverting approximate expression used.
(2) prior imformation of Bayes's inverting and likelihood function are set up;
1. prior imformation is set up:
Bayes in Bayesian frame, prior model definition be the statistical models of elastic parameter prior imformation, expressed by the form of probability density function.Meanwhile, prior imformation independent of geological data, can only must be obtained by other channel (as well logging information and geologic knowledge).Acoustic medium speed c and attenuation factor=1/Q is glued to describe underground medium for Chang Midu, and supposes that parameter { ln c (t), α (t) } meets Gaussian distribution.Adopt this hypothesis to be on the one hand because it compares to meet the fact, a large amount of well-log informations shows that the logarithm of speed and decay factor are all close to Gaussian distribution; On the other hand, the hypothesis of Gaussian distribution can bring great convenience to calculating.
First a continuous Gaussian vector field be made up of elastic parameter is defined:
m(t)=[ln c(t),α(t)]
T(7)
Being contemplated to be of it:
E[m(t)]=μ(t)=[μ
c,μ
α]
T(8)
Lnc (t), α (t) in the variance of time sampling point t and s are:
Cov[m(t),m(s)]=∑(t,s)=∑
0v
t(τ) (9)
Wherein v
t(τ) be a time correlation function, that τ represents is time delay τ=t-s, and ∑
0represent be one time constant covariance matrix:
In this matrix, diagonal entry represents the variance of wave velocity and decay factor respectively, v
c αthe related coefficient of representation speed and density.Time correlation function v
t(τ) must be a positive definite function, value be in [-1,1], and v
t(0)=1.Such as bi-exponential function can be expressed as:
Wherein d is the parameter of portraying temporal correlation.Thus, discrete elastic parameter prior imformation can be write as:
Wherein, m is the vector of elastic parameter composition, and lnc represents respectively and takes the logarithm to speed, and α represents decay factor,
represent n the sampling point Normal Distribution about m, μ
m, ∑
mrepresent the expectation and variance of m respectively.
2. likelihood function is set up:
Based on sticky acoustic medium AVA/AVF approximate expression, the equation both sides of formula (6) are obtained to time t differentiate the situation that this equation is extended to interface continuous time:
In Bayes AVA/AVF inverting, likelihood function is defined as and obtains seismic observation data d when setting models parameter m
obsprobability, can be write as p (d
obs| m).It is described that to what extent these model parameters can obtain seismic observation data.Be written as with a matrix type:
R=Am′ (14)
Wherein m ' is for m is to the derivative of time, and A is the sparse matrix be made up of coefficient in equation,
For AVA inverting:
For AVF inverting:
Just drilling the seismologic record obtained to be expressed as:
d
obs=SR+e=SAm′+e (16)
Wherein e is the noise in seismologic record, and S is:
Can be write as each S:
Assuming that it is 0 that noise e meets the expectation, variance is ∑
egaussian distribution:
such observation data also meets Gaussian distribution, and likelihood function is written as:
Formula (18) is exactly the likelihood function used by AVA and AVF, wherein, and d
obsrepresent observation geological data, m represents elastic parameter,
represent that observation geological data meets Gaussian distribution,
represent the expectation of geological data,
represent the variance of geological data.So-called likelihood function is when known underground parameter m, obtains observation geological data d
obsprobability have much.Owing to meeting Gaussian distribution, therefore likelihood function here can be expressed with expectation and variance, and expectation and variance meets respectively:
(3) for the linear AVA/AVF inverting of Bayes of speed and decay;
Seismologic record and elastic parameter meet joint distribution:
Wherein
At given d
obswhen m Posterior probability distribution be:
Posterior error and variance are respectively:
The maximum a posteriori solution (MAP) of m just should equal posterior error value
covariance matrix is evaluated in order to the uncertainty final to inversion result.
As shown in Figure 7, the inventive method comprises:
The first step, input angle territory geological data, frequency field geological data and log data; The geological data of angle domain geological data and frequency field is all expressed as d
obs, these data are the 3-D data volume of time, angle and frequency, and user can extract corresponding data out when carrying out AVA and AVF inverting, d
obsspecifically angle domain geological data or frequency field geological data, depends on user and uses what inversion method.Log data is expressed as: m=[lnc, α ,]
t, m is the vector of elastic parameter composition, and lnc represents the logarithm of speed respectively, and α represents decay factor, and the two all obtains from well-log information;
Second step, calculates the prior imformation of log data;
3rd step, utilizes prior imformation and likelihood function, carries out linear AVA and the AVF inverting of Bayes, obtains speed and quality factor.
Below by a synthesis road collection, effect of the present invention is described, algorithm C language is write, and picture Matlab makes.Model is accompany sandstone reservoir (as shown in Figure 1) in the middle of shale layer, and the parameter of three bed interfaces is set to respectively: Vp1=300m/s, Qp1=100; Vp2=4000, Qp2=10; Vp3=3500, Qp3=80.Wavelet is ricker wavelet, dominant frequency 35hz; Time sampling interval is: 2ms; Angular range is from 0 °-45 ° changes, and frequency range is: 2hz-80hz.
First the comparative analysis error (as Suo Shi Fig. 2, Fig. 3 and Fig. 4) of sticky acoustic medium AVA/AVF approximate expression and accurate reflection coefficient, result is presented at that the upper and lower velocity contrast of reflecting interface is little, quality factor poor little when, degree of approximation is very high.And these 2 approximate conditions are all easier to meet, the sticky acoustic medium AVA/AVF approximate expression of therefore deriving has very high precision and using value.Then, the record (as shown in Fig. 5-1 and Fig. 5-2) of seismic amplitude with incident angle and frequency change is obtained by convolution model.Next amplitude is utilized to carry out the linear AVA/AVF inverting of Bayes at incident angle domain information and amplitude in frequency-domain information respectively, obtain speed and the decay factor change (as shown in Fig. 6-1 and Fig. 6-2) of underground medium, the value of initial model is the average of medium velocity.Be not difficult to find out, the sticky linear AVA inverting of acoustic medium Bayes can be finally inversed by the conversion of medium velocity accurately, but limited for the inversion accuracy of decay factor; Bayes linear AVF inverting then point-devicely can be finally inversed by speed and the decay factor of underground medium.
Along with the aggravation of oil-gas exploration difficulty, people have to sight to turn to more hidden lithologic deposit from structural deposit.Therefore, traditional structure imaging has been difficult to the demand of current oil-gas exploration.And the wide-azimuth acquisition technique developed on this basis and relative amplitude preserved processing technology all achieve significant progress, also played a great role in reservoir prediction by the prestack inversion comprising AVO technology simultaneously and provided guarantee.But traditional AVO analyzing total is hypothesis underground medium is perfectly elastic.Decay and velocity dispersion phenomenon can be there is in the Petrophysical measurement in VSP data, well-log information and laboratory all explicitly seismic wave in actual propagation process.Especially for the region containing hydrocarbon, decay clearly.Ignore seismic wave attenuation and can bring huge risk to AVO analysis and reservoir prediction.
The attenuation by absorption of seismic event is brought into whole inverting framework and is got off by the present invention, make use of the information of seismic reflection data angle domain and frequency field simultaneously, it is also combined with Bayes's inversion method and builds new prior imformation and likelihood function by the Reflection Efficient Approximation of frequency of having derived angle domain, inversion result not only can obtain the change of sticky acoustic medium medium velocity, can obtain the change of underground medium decay factor simultaneously.This invention is not only outstanding to seism processing work meaning, more finds subterranean oil gas reservoir and provides important instruction.
Technique scheme is one embodiment of the present invention, for those skilled in the art, on the basis that the invention discloses application process and principle, be easy to make various types of improvement or distortion, and the method be not limited only to described by the above-mentioned embodiment of the present invention, therefore previously described mode is just preferred, and does not have restrictive meaning.
Claims (4)
1. linear AVA and the AVF inversion method of sticky acoustic medium Bayes, is characterized in that: described method comprises:
The first step, input angle territory geological data, frequency field geological data and log data
Second step, sets up prior imformation and the likelihood function of log data;
3rd step, utilizes described prior imformation and likelihood function, carries out linear AVA and the AVF inverting of Bayes, obtains speed and quality factor.
2. linear AVA and the AVF inversion method of sticky acoustic medium Bayes according to claim 1, is characterized in that: in the described first step, the geological data of angle domain geological data and frequency field is all expressed as d
obs, these data are the 3-D data volume of time, angle and frequency; Log data is expressed as: m=[lnc, α ,]
t, wherein, m is the vector of elastic parameter composition, and lnc represents the logarithm of speed respectively, and α represents decay factor, and the two all obtains from well-log information.
3. linear AVA and the AVF inversion method of sticky acoustic medium Bayes according to claim 2, is characterized in that: described second step is achieved in that
The prior imformation setting up log data is as follows:
Wherein,
represent n the sampling point Normal Distribution about m, μ
m, ∑
mrepresent the expectation and variance of m respectively, by what obtain the statistics of well logging sampling point;
Set up likelihood function as follows:
Wherein,
represent that observation geological data meets Gaussian distribution,
represent the expectation of geological data,
represent the variance of geological data;
Expectation and variance meets respectively:
∑ "
mrepresenting asks second order to lead to the statistical variance of elastic parameter, ∑
erepresent the variance of noise;
Wherein, S is:
Can be write as each S:
For AVA inverting, adopt:
For AVF inverting, adopt:
Wherein, s
1(θ
i, ω
i) expression incident angle is θ
i, dominant frequency is ω
iwavelet; Subscript 1,2 until nwave represents the time delay of wavelet; A
c(θ
n θ, ω
fix) represent frequencies omega
fixfixing, depend on incidence angle θ
n θvelocity coefficient; A
α(θ
n θ, ω
fix) represent frequencies omega
fixfixing, depend on incidence angle θ
n θdecay factor coefficient; A
c(θ
fix, ω
n ω) represent incidence angle θ
fixfixing, depend on frequencies omega
n ωvelocity coefficient; A
α(θ
fix, ω
n ω) represent incidence angle θ
fixfixing, depend on frequencies omega
n ωdecay factor coefficient.
4. linear AVA and the AVF inversion method of sticky acoustic medium Bayes according to claim 3, is characterized in that: described 3rd step is achieved in that
Seismologic record and elastic parameter meet joint distribution:
Wherein
At given d
obswhen m Posterior probability distribution be:
Posterior error and variance are respectively:
The maximum a posteriori solution (MAP) of m just equals posterior error value
covariance matrix is evaluated in order to the uncertainty final to inversion result.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310476278.1A CN104570071B (en) | 2013-10-12 | 2013-10-12 | A kind of linear AVA and AVF inversion methods of viscous acoustic medium Bayes |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310476278.1A CN104570071B (en) | 2013-10-12 | 2013-10-12 | A kind of linear AVA and AVF inversion methods of viscous acoustic medium Bayes |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104570071A true CN104570071A (en) | 2015-04-29 |
CN104570071B CN104570071B (en) | 2017-08-18 |
Family
ID=53086585
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310476278.1A Active CN104570071B (en) | 2013-10-12 | 2013-10-12 | A kind of linear AVA and AVF inversion methods of viscous acoustic medium Bayes |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104570071B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112558151A (en) * | 2020-12-22 | 2021-03-26 | 同济大学 | AVA waveform simulation method and system in attenuation medium |
CN113219536A (en) * | 2021-06-25 | 2021-08-06 | 成都理工大学 | Pre-stack seismic inversion method of longitudinal and transverse wave attenuation parameters depending on frequency |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005066660A2 (en) * | 2004-01-09 | 2005-07-21 | Statoil Asa | Processing seismic data representing a physical system |
CN102004263A (en) * | 2010-09-17 | 2011-04-06 | 中国石油大学(华东) | Angular part superposition seismic data-based fluid identification method |
CN103293552A (en) * | 2013-05-24 | 2013-09-11 | 中国石油天然气集团公司 | Pre-stack seismic data retrieval method and system |
-
2013
- 2013-10-12 CN CN201310476278.1A patent/CN104570071B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005066660A2 (en) * | 2004-01-09 | 2005-07-21 | Statoil Asa | Processing seismic data representing a physical system |
CN102004263A (en) * | 2010-09-17 | 2011-04-06 | 中国石油大学(华东) | Angular part superposition seismic data-based fluid identification method |
CN103293552A (en) * | 2013-05-24 | 2013-09-11 | 中国石油天然气集团公司 | Pre-stack seismic data retrieval method and system |
Non-Patent Citations (4)
Title |
---|
ARILD BULAND 等: ""Bayesian Linearized AVO Inversion"", 《GEOPHYSICS》 * |
KRISTOPHER A. INNANEN: ""Inversion of the seismic AVF/AVA signatures of highly attenuative targets"", 《GEOPHYSICS》 * |
LONG TENG 等: ""Probabilistic petrophysical properties estimation integrating viscoelastic AVA inversion with statistical rock physics"", 《SEG TECHNICAL PROGRAM EXPANDED ABSTRACTS 2012》 * |
黄捍东 等: ""基于贝叶斯理论的薄层反演方法",,", 《石油地球物理勘探》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112558151A (en) * | 2020-12-22 | 2021-03-26 | 同济大学 | AVA waveform simulation method and system in attenuation medium |
CN113219536A (en) * | 2021-06-25 | 2021-08-06 | 成都理工大学 | Pre-stack seismic inversion method of longitudinal and transverse wave attenuation parameters depending on frequency |
CN113219536B (en) * | 2021-06-25 | 2022-03-01 | 成都理工大学 | Pre-stack seismic inversion method of longitudinal and transverse wave attenuation parameters depending on frequency |
Also Published As
Publication number | Publication date |
---|---|
CN104570071B (en) | 2017-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yuan et al. | Multiscale adjoint waveform tomography for surface and body waves | |
Dasgupta et al. | Estimation of Q from surface seismic reflection data | |
CN104090302B (en) | The method of work area underground medium frequency domain anomaly analysis | |
Larsen | AVO inversion by simultaneous PP and PS inversion | |
CN105388518A (en) | Centroid frequency and spectral ratio integrated borehole seismic quality factor inversion method | |
CN103293552A (en) | Pre-stack seismic data retrieval method and system | |
CN104237938B (en) | Method for matching longitudinal wave and converted wave data through dynamic time adjustment | |
Chen et al. | Numerical simulation of frequency-dependent seismic response and gas reservoir delineation in turbidites: A case study from China | |
CN105044777A (en) | Method for detecting strong reflection amplitude elimination of seismic marker layer based on empirical mode decomposition | |
CN110618450B (en) | Intelligent gas-bearing property prediction method for tight reservoir based on rock physical modeling | |
Pan et al. | Matrix-fluid-fracture decoupled-based elastic impedance variation with angle and azimuth inversion for fluid modulus and fracture weaknesses | |
CN110261897A (en) | Based on four parameter inversion method of prestack that group is sparse | |
Pafeng et al. | Prestack waveform inversion of three-dimensional seismic data—An example from the Rock Springs Uplift, Wyoming, USA | |
Leiceaga et al. | Enhanced density estimation from prestack inversion of multicomponent seismic data | |
CN104391324A (en) | Seismic trace set dynamic correction stretching correction pre-processing technology before AVO inversion depending on frequency | |
Rubino et al. | Thin-bed prestack spectral inversion | |
CN110187389A (en) | A kind of AVA inversion method based on thin bed reflection theory | |
CN109581500B (en) | Reflection seismic record frequency-variable velocity analysis method | |
CN104570071A (en) | Bayesian linear AVA and AVF retrieval method of sticky sound medium | |
Sun et al. | The stability problem of reverse time migration for viscoacoustic VTI media | |
CN106226814B (en) | Utilize converted shear wave seismic data inversion reservoir S-wave impedance and the method for density | |
Yang et al. | Mitigating velocity errors in least-squares imaging using angle-dependent forward and adjoint Gaussian beam operators | |
CN106249294A (en) | A kind of reservoir detecting method of hydrocarbon | |
CN105929453B (en) | State estimation method for infinite distribution time lag neural network system with channel fading | |
Wu et al. | An unsupervised inversion method for seismic brittleness parameters driven by the physical equation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |