用非零偏P-P波、P-SV波计算零偏反射纵波和反射横波

在P波入射时，地下同一点产生的P－P波和P－SV波是入射角、界面两侧纵波阻抗和横波阻抗的函数，并且两者是耦合的，利用两者的耦合关系，可在入射角确定的情况下提取与入射角无关的零偏纵波反射和横波反射。本文对此进行了理论推导，结论是：入射角为零时的纵波与横波反射系数是入射角非零时的P－P波和P－SV波反射系数的线性组合，进而可以用纵横波反射计算泊松比反向间面，把多波叠前反演转化到叠后进行，通过合成弹性参数反射剖面，直接反演弹性参数等，都很有意义。文中给出了有关的公式，并对具体实现方法和有关参数进行了讨论。     

基于等效Thomsen参数的P-SV波AVO属性研究

在等效Thomsen各向异性参数的P-SV波反射系数近似公式基础上研究了反射系数的多种AVO(振幅随炮检距变化)属性特征,针对不同的属性特征构建了多属性AVO交绘图;并利用反射系数公式对三类含气砂岩AVO的特征进行分析.结果表明P-SV波反射系数公式可以有效的区分第三类含气砂岩;岩石的孔隙度、流体饱和度等信息是影响地震波AVO的重要因素.利用Gassmann方程进行了对上层为HTI介质,下层为孔隙储层的介质模型进行了流体替换计算,分析了孔隙度、含气饱和度和各向异性参数变化对地层AVO的影响.     

转换波AVO近似及其在PP/PS联合反演中的应用(英文)

多分量勘探相比常规纵波勘探而言有许多的优势,PP/PS联合AVO分析和反演是强而有效的储层识别方法。在本文中,我们推导了更精确的转换波AVO反射系数公式;并在入射角小于30度时进行了进一步简化,其对于转换波AVO分析和反演而言更加简单有效。基于该近似,我们进行了PP/PS联合AVO反演。实际资料的实例表明,反演得到的纵横波速度比可很好地识别岩性及油气。反演得到的流体因子和泊松比等其他属性在储层处也显示出明显异常,效果显著。     

基于贝叶斯理论的叠前多波联合反演弹性模量方法

AVO反演可以获得地层岩性和流体信息,而叠前反演问题都是高维的和非适定的,因此获得可靠稳定的解对叠前反演至关重要. 本文给出了一种基于贝叶斯理论的纵波和转换波联合反演密度比和模量比的方法. 鉴于剪切模量比、体积模量比可以更好地指示油气,基于岩石物理中速度比与模量比之间的关系,将此关系式代入Zoeppritz方程的近似形式Aki-Richards公式中,得到与模量比有关的反射系数近似公式. 联合纵波和转换波,利用最小二乘准则构建目标函数,最终反演出密度比、剪切模量比、体积模量比三个参数. 在反演过程中引入贝叶斯理论,假定先验信息服从高斯分布,待求参数服从改进的Cauchy分布,并去除待求参数之间的相关性. 利用模型数据和实际数据对本文方法进行测试,并与常规的单独利用纵波数据来反演方法进行比较,结果表明联合反演稳定性更好、精度更高、抗噪音能力更强,验证了本文方法的可行性和有效性.     

基于精确Zoeppritz方程三变量柯西分布先验约束的广义线性AVO反演

常规AVO三参数反演是通过Zoeppritz方程的近似公式来建立AVO正演模拟的过程,然而在P波入射角过临界角和弹性参数在纵向上变化剧烈的情况下,Zoeppritz方程近似公式精度有限.针对这种情况,可以使用精确的Zoeppritz方程来构建反演目标函数,由于精确Zoeppritz方程中P波反射系数和弹性参数之间是一种复杂的非线性关系,通常解决途径是利用非线性的优化算法来进行数值计算,但是非线性优化算法的缺点是计算量过大;另外一种途径是利用广义线性反演的方法,通过泰勒一阶展开式将P波反射振幅展开后,用线性关系近似表达非线性关系,经过几次迭代后,在理论上可以达到很高的精度,但是广义线性反演算法的核心部分--Jacobian矩阵由于矩阵条件数过大,往往会造成反演算法的不稳定,其应用范围得到了限制.贝叶斯反演方法是通过引入模型参数的先验分布结合噪声的似然函数,生成模型参数的后验分布,通过求取模型参数的最大后验概率分布来得到模型参数的反演解,由于引入模型参数的先验分布信息,可以有效的降低反演的不适定问题.本文将两种反演算法的思想相结合,利用广义线性反演算法的思想,构建AVO正演模拟的过程来提高大角度地震数据反演的精度,同时结合贝叶斯理论,通过引入模型参数的先验分布信息构建反演目标函数的正则化项,可以有效降低由于Jacob矩阵条件数过大带来的反演不适定问题,该算法假设模型参数服从三变量柯西分布.     

P-P波和P-SV波联合AVO反演研究

目前转换波AVO(Amplitude Versus Offset)反演或纵波、转换波联合AVO反演是基于Aki等的反射系数近似公式,通常反演结果为上、下层速度比值和速度平均比值.本文在Shuey和郑晓东近似公式的基础上,运用最小平方算法,提出一种p-p波和P-SV波联合AVO反演下层介质纵、横波速度及泊松比的方案,并应用于双层含气砂岩模型上,误差在10%以内.此方案要求上、下层物性差别不大且入射角小于30°,但不需要速度比值为常数2的假设.     

TTI介质各向异性参数多波反演与PS波AVO分析

把遗传算法引入到了TTI介质AVO信息反演各向异性参数的过程中,依据TTI介质PP波、PS波反射系数公式,建立Thomsen参数和TTI介质对称轴倾角、方位角的目标函数,分别通过PP波和PS波的反射系数反演出了各向异性参数和对称轴倾角、方位角等信息.文中对反演结果的精确度和稳定性进行了分析,发现PS波的反演结果优于PP波反演结果;对称轴倾角的反演准确性明显优于对称轴方位角.本文通过模型正演合理解释了这一现象的原因.最后,本文通过对PS波AVO梯度的研究,提出了利用PS波振幅定性分析TTI介质对称轴倾角的方法.     

基于精确Zoeppritz方程的非线性AVO三参数反演

常规AVO三参数反演通常存在密度反演不准确的问题,而密度参数对常规油气藏中的流体识别、流体饱和度计算、孔隙度计算以及非常规油气藏中TOC含量计算、裂缝预测等都至关重要,因此对于研究如何利用大偏移距振幅信息和富含密度信息的PS波地震资料来提高密度反演结果的稳定性和精度显得尤为重要.研究基于贝叶斯反演理论框架,引入三变量Cauchy分布先验约束,利用精确Zoeppritz方程构建了AVO三参数联合反演的目标函数,对目标函数进行Taylor二阶非线性简化,得到模型参数的迭代更新公式,实现了大偏移距地震振幅信息的利用和PP波、PS波联合反演.合成数据和实际地震数据的方法测试结果表明,新方法不仅可以直接反演纵波速度、横波速度和密度,而且还具有很高的精度,尤其是密度反演结果.基于合成数据的PP波、PS波单独反演结果与PP波和PS波联合反演结果对比显示,联合反演稳定性更好,精度更高,抗噪能力更强,验证了该方法的可行性和有效性.与基于Aki-Richards近似公式的反演结果对比表明,该反演方法具有更高的反演精度和更好的抗噪性.     

AVF反演Q值方法研究

品质因子Q的准确估计是叠前、叠后地震数据处理及AVO分析的关键环节.理论分析和正演模拟表明界面处单纯Q值的变化也会引起反射,该反射系数随频率变化而变化(Amplitude Versus Frequency,AVF).利用其线性近似公式,通过AVF反演可以准确地获得品质因子Q.AVF反演方法为Q值的求取提供了新的思路,与传统Q值求取方法不同,该方法主要考虑衰减对反射系数的影响,通过反射系数获得目的层的Q值.将上覆层为弹性介质的假设扩展到非弹性介质,构建新的近似公式,该近似公式应用更具普遍性与适应性,并利用时频分析方法获得多层模型介质的Q值;测试结果显示AVF反演的Q值准确、稳定.     

An approximation to the P-wave reflection coefficient for amplitude variation with long-offset seismic data

近年来,长偏移距排列逐渐被大量应用到海上勘探。传统AVO分析都是小入射角,对于长偏移距数据来说有很大不足。为了探讨拓宽AVO分析的研究范围,本文推导了长偏移距时P波反射系数的一种近似公式,并讨论了其精确性。本文的结果和广泛应用的Shuey公式相类似,Shuey公式可以看作本文结果在小入射角情况的简化。     

弱各向异性介质中的P-SV波近似反射系数计算

We introduce the Thomsen anisotropic parameters into the approximate linear reflection coefficient equation for P-SV wave in weakly anisotropic HTI media. From this we get a new, more effective, and practical reflection coefficient equation. We performed forward modeling to AVO attributes, obtaining excellent results. The combined AVO attribute analysis of PP and PS reflection data can greatly reduce ambiguity, obtain better petrophysical parameters, and improve parameter accuracy.     

Data acquisition and pre-processing required for simultaneous P-SV inversion

The goal of seismic reflection surveys is the derivation of petrophysical subsurface parameters from surface measurements. Today's well established technique in data acquisition, as well as processing terms, is based on the acoustic approximation to the real world's wave propagation. In recent years a lot of work has been done to extend the technique to the elastic approximation. There was especially an important trend towards elastic inversion techniques operating on plane-wave seismograms, called simultaneous P-SV inversion (or short P-SV inversion) within this paper. Being still under investigation, some important aspects of P-SV inversion concerning data acquisition as well as pre-processing, should be pointed out. To fit the assumptions of P-SV inversion schemes, at least a two-dimensional picture of the reflected wavefield with vertical and in-line horizontal receivers has to be recorded. Moreover, the theoretical work done suggests that in addition to a survey with a compressional wave source, a second survey should be done using sources radiating vertically polarized shear waves, is needed. Finally, proper slant stacking must be performed to get plane-wave seismograms. The P/S separated plane-wave seismograms are then well prepared for feeding into the inversion algorithms. P/S separated planewave seismograms are then well prepared for feeding into the inversion algorithm.s In this paper, a tutorial overview of the data acquisition and pre-processing in accordance with the P-SV inversion philosophy is given and illustrated using synthetic seismograms. A judgement on the feasibility of the P-SV inversion philosophy must be left to ongoing research.     

弱各向异性介质中P-SV波弹性阻抗与流体识别因子(英文)

基于Petr Jilek正交各向异性介质中的P-SV波反射系数精确公式,通过近似得到VTI和HTI介质中的P-SV波近似公式,同时结合Kenneth横波弹性阻抗原理,推导建立了VTI和HTI介质中P-SV波的弹性阻抗近似公式;在此基础上构建了弱各向异性介质中的流体识别因子公式,通过理论模型对Castagna和Smith的岩性组合进行了识别;最后详细分析了各向异性参数的取值对P-SV波弹性阻抗及识别因子产生的影响及变化趋势。     

基于Russell近似的纵横波联合反演方法研究

PP波和PS波联合反演方法作为有效的地震技术,比单纯纵波反演精度要高,能够提高地震储层识别的精度.以Russell近似理论为基础,推导了新的转换波AVO近似公式,双层模型界面的反射特征数值模拟显示,新公式具有较高的近似精度,且具备直接反演流体因子f、剪切模量μ和密度ρ等参数的优势,有效避免间接反演带来的误差.结合纵横波联合反演理论,提出了基于贝叶斯理论的新型联合反演算法.在实际应用中,对纵波和转换波角道集进行同相轴匹配处理,综合利用纵波和转换波资料携带的信息,实现基于Russell近似的多波联合反演.模型数据和实际资料测试结果表明,反演结果与真实值或测井结果匹配度较高,证实该方法真实有效.     

Analytical wavefront curvature correction to plane‐wave reflection coefficients for a weak‐contrast interface

Most amplitude versus offset (AVO) analysis and inversion techniques are based on the Zoeppritz equations for plane‐wave reflection coefficients or their approximations. Real seismic surveys use localized sources that produce spherical waves, rather than plane waves. In the far‐field, the AVO response for a spherical wave reflected from a plane interface can be well approximated by a plane‐wave response. However this approximation breaks down in the vicinity of the critical angle. Conventional AVO analysis ignores this problem and always utilizes the plane‐wave response. This approach is sufficiently accurate as long as the angles of incidence are much smaller than the critical angle. Such moderate angles are more than sufficient for the standard estimation of the AVO intercept and gradient. However, when independent estimation of the formation density is required, it may be important to use large incidence angles close to the critical angle, where spherical wave effects become important. For the amplitude of a spherical wave reflected from a plane fluid‐fluid interface, an analytical approximation is known, which provides a correction to the plane‐wave reflection coefficients for all angles. For the amplitude of a spherical wave reflected from a solid/solid interface, we propose a formula that combines this analytical approximation with the linearized plane‐wave AVO equation. The proposed approximation shows reasonable agreement with numerical simulations for a range of frequencies. Using this solution, we constructed a two‐layer three‐parameter least‐squares inversion algorithm. Application of this algorithm to synthetic data for a single plane interface shows an improvement compared to the use of plane‐wave reflection coefficients.     

PP and PS joint AVO inversion and fluid prediction

The technique of amplitude variation with offset or angle (AVO or AVA) can be used to extract fluid and lithology information from prestack seismic data. Based on three-term AVO equations, three elastic parameters can be inverted for by linear AVO inversion. However, many theoretical and numerical studies have demonstrated that by using offset limited data, a three-term AVO inversion may have problems of instability and inaccuracy while inverting for the density term. We have searched for an elastic parameter that contains density information and inverted this parameter in a more stable manner using offset limited data. First, we test the sensitivity of elastic parameters to hydrocarbon reservoirs and select the optimal fluid factor (ρf) that contains density information and has an excellent performance as an inversion parameter used to detect hydrocarbons. Then, we derive approximate PP and PS reflection coefficient equations in terms of the fluid factor. The derived equations allow us to directly estimate the fluid factor of the reservoir. Finally, we apply these equations to synthetic data by employing a joint AVO inversion technique. The results show that the method is stable and unambiguous.     

双相介质的AVO资料中蕴含的流体信息的定量分析

振幅随偏移距变化(AV0)弹性参数的反演是储层预测的一个重要方法.利用正演技术分析双相介质的地震资料中含有的流体信息,有助于提高AVO弹性参数反演的可靠性.本文应用Russell推出的计算纵波反射系数的AVO线性近似公式,对双相介质的纵波资料中所含有的流体信息进行定量分析,研究从AV0资料的反演中可以确定哪些弹性参数或它们的组合关系.文中考虑的弹性参数有流体因子、剪切模量和密度三者的相对扰动量.在时间城内对反射系数的正演计算公式进行线性化,得到灵敏度矩阵,并用奇异值分解(SVD)计算其特征值和特征向量采分析弹性参数的反演.由计算结果可知,在入射角较小时由AV0资料反演只能得到流体因子和密度两者的相对扰动量;若想同时获得三个参数的相对扰动量,就必须增大入射角.     

Nonlinear simultaneous inversion of pore structure and physical parameters based on elastic impedance

Carbonate reservoirs have complex pore structures, which not only significantly affect the elastic properties and seismic responses of the reservoirs but also affect the accuracy of the prediction of the physical parameters. The existing rockphysics inversion methods are mainly designed for clastic rocks, and the inversion objects are generally porosity and water saturation. The data used are primarily based on the elastic parameters, and the inversion methods are mainly linear approximations. To date, there has been a lack of a simultaneous pore structure and physical parameter inversion method for carbonate reservoirs. To solve these problems, a new Bayesian nonlinear simultaneous inversion method based on elastic impedance is proposed. This method integrates the differential effective medium model of multiple-porosity rocks, Gassmann equation,Amplitude Versus Offset(AVO) theory, Bayesian theory, and a nonlinear inversion algorithm to achieve the simultaneous quantitative prediction of the pore structure and physical parameters of complex porous reservoirs. The forward modeling indicates that the contribution of the pore structure, i.e., the pore aspect ratio, to the AVO response and elastic impedance is second only to that of porosity and is far greater than that of water saturation. The application to real data shows that the new inversion method for determining the pore structure and physical parameters directly from pre-stack data can accurately predict a reservoir's porosity and water saturation and can evaluate the pore structure of the effective reservoir.     

页岩气储层纵横波叠前联合反演方法

杨氏模量与密度乘积(Eρ)能够突显页岩气储层的异常特征,泊松比能够指示储层的含流体性.与常规叠前弹性参数反演相比,基于Eρ、泊松比和密度的叠前纵横波联合反演可以获得更加精确的弹性参数,为页岩气储层识别和流体预测提供可靠的依据.首先,推导了基于Eρ、泊松比和密度的纵波和转换波反射系数近似方程,利用典型模型对新推导的反射系数方程做精度分析,当入射角小于30°时,新推导的反射系数公式具有较高的精度;其次,充分利用纵波和转换波的信息,在贝叶斯的框架下,建立叠前纵横波联合反演流程,进行Eρ、泊松比和密度的直接反演,避免了间接反演带来的累积误差;最后,利用实际工区井模型数据进行算法测试,结果表明,该反演方法所获得的Eρ、泊松比和密度的估计值与真实值吻合较好,满足精细地震反演的精度要求.