基于入射角的两项流体阻抗反演方法

本文在孔隙弹性介质理论的指导下,基于入射角AVO近似方程推导了包含Russell流体项的两项AVO近似方程和相应的弹性阻抗方程,通过分析可知其精度符合反演要求.在贝叶斯理论框架下,建立了包含正则化约束的弹性阻抗反演方法,在此基础上直接提取Russell流体项.该方法可在缺少大角度叠前地震资料的情况下进行叠前直接反演得到流体因子,减少传统方法带来的累积误差.模型试算表明,该方法具有较好的准确度和稳定性.实际工区应用取得了良好的效果,表明该方法有实用性.     

基于对数变换的扩展弹性阻抗技术在岩性识别中的应用

叠后地震数据,是共反射点道集叠加的结果,虽然提高了信噪比,但是也损失了丰富的叠前地震道集信息,建立在叠后数据基础上的叠后波阻抗反演不可避免的掩盖了很多AVO类信息.实际地震数据的采集都是非零炮检距的炮集记录,这种叠前地震数据提供了多种波型地震信息,对流体和岩性更加敏感.因此,利用叠前道集或者部分道集信息的叠前地震反演能够大大提高预测的精度,扩展弹性阻抗技术不仅使用了丰富的叠前地震道集信息,同时还解决了常规弹性阻抗技术不稳定和实用性的问题,使叠前反演更具有实际应用价值.通过弹性阻抗公式和扩展弹性阻抗公式的对比分析,验证了扩展弹性阻抗反演更加稳定.在利用扩展弹性阻抗技术进行反演前先对扩展弹性阻抗公式进行对数变换,建立梯度阻抗和声波阻抗的坐标系统,通过计算最佳旋转角度下的弹性阻抗,进行岩性预测.同时指出了计算声阻抗和梯度阻抗的方法,并分析了岩性预测基本原理.应用扩展弹性阻抗技术在我国西部某气田河道砂体识别中取得较好的效果.     

AVO反演的不确定性分析

叠前地震数据反演可以得到比常规叠后波阻抗反演更丰富、更有效的岩性信息,但叠前数据体的噪声及其它因素严重影响了AVO反演的稳定性,如何评估AVO反演结果的可靠性显得尤为重要.本文从贝叶斯理论出发,假定均匀先验分布、高斯噪音分布,推出不确定性分析方程,利用协方差矩阵中的方差描述反演问题的不确定性,模型研究显示反演不确定性与叠前信噪比、纵横波速度比、覆盖次数及反演采用的角度范围相关,方法预测的反演误差可定量解释反演结果的可靠性,为解释人员提供有效的质量监控手段.     

非线性AVO反演方法研究

与叠后地震数据相比,叠前地震数据包含有更多的反映地下地层特征的信息,利用AVO( Amplitude Versus Offset,振幅随偏移距的变化)信息通过求解Zoeppritz方程的近似公式,叠前反演可直接得到反映地下岩石特征的弹性参数——密度、纵波速度和横波速度.从本质上讲,叠前地震反演是非线性的,但目前多采用线...     

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

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

基于改进反射系数近似方程的纵横波阻抗同步反演

基于地震波反射系数近似公式的叠前反演是油气勘探的重要工具.本文在已有研究的基础上,推导了一个改进的射线参数域地震纵波反射系数近似方程.该方程建立了地震纵波反射系数与纵波阻抗和横波阻抗的非线性关系,在中、小角度的范围内较现有的反射系数线性近似公式精度更高.另外,由于该方程仅包含纵波和横波阻抗反射系数项,因此基于新方程的反演能够有效地降低同步反演纵波速度、横波速度、密度三个参数的不适定性.在此基础上,结合广义线性反演法(GLI)理论和贝叶斯理论,相应地发展了一种叠前地震同步反演方法.模型测试和实际资料的应用表明,基于新方程的反演方法能够利用有限角度(偏移距)的数据稳定地反演纵波和横波阻抗,由于在反演过程中,不需要假设纵横波速度为常数,因此该方法还能有效地提高反演结果的精度.     

LD地区致密气叠前关键技术及应用

鄂尔多斯盆地LD地区为典型的黄土塬地貌,其砂体规模小、横向变化大、埋藏深、厚度薄、储层致密.以高保真地震分偏移距叠加资料为基础,在岩石物理分析的指导下,创新探索引入了几项关键叠前技术:球面波反射AVO分析技术研究山1段大偏移距道集AVO特征;叠前角度域中分析地震资料吸收规律,精度相对较高,能分辨5m储层;引入果蝇优化算法构建虚拟井,解决叠前同时反演过程中缺少横波测井问题.这几项黄土塬叠前关键技术为LD地区深层致密气勘探取得新进展发挥了重要作用.     

多波时移地震AVO反演研究

数值模拟了油藏含油饱和度与有效压力变化时移地震AVO的响应,确定利用时移地震AVO区分油藏参数的变化、实现油藏定量解释的可行性.从Aki等 AVO近似方程出发,详细推导了P_P波和P_S转换波时移地震AVO计算公式.结合岩石物理近似关系和本文推导的时移地震AVO计算公式,推导了利用多波时移地震AVO反演油藏含油饱和度和压力变化的方程.数据试验表明,文中推导的多波时移地震AVO方程能较好地反演油藏含油饱和度变化和有效压力变化,实现油藏定量解释.     

P-P、P-SV波叠前联合反演在致密气储层预测中的应用

常规叠前反演大都是利用P-P纵波单分量资料基于各向同性理论的近似公式来反演,仅利用P-P纵波资料反演地层岩性等弹性参数,无论在理论基础还是反演精度上都存在不足,且反演结果具有多解性.而同时用P-P纵波和P-SV转换波资料联合反演,可以增加约束条件,提高反演结果的稳定性、可靠性.本文以各向异性理论为基础,利用实际3D3C资料,进行了P-P纵波、P-SV转换波联合叠前反演计算,P-P纵波、P-SV转换波叠前联合反演利用了纵波和转换波的叠前资料,尽可能保留了地层的AVO信息,输入数据增加了转换横波数据,同时利用了纵波和转换波两个独立变量信息,采用完全Knott-Zoeppritz方程求解,反演结果比P-P纵波单分量资料叠前反演更稳定、可靠、反演精度明显提高.研究区实际资料计算结果表明,反演结果可以更真实的反映储层段的物性和含气性,P-P纵波、P-SV转换波联合叠前反演比P-P纵波单分量资料叠前反演预测储层符合率提高了十个百分点,该反演方法是一种有效的反演方法.     

基于叠前弹性反演的含气储层预测技术

叠前弹性反演是目前岩性气藏储层预测的主要技术,仅仅使用纵波阻抗反演预测储层在理论和准确性方面都存在着一些不足,反演的多解性更是其致命的缺陷.本文从横波测井资料不同岩性、流体的岩石物理参数分析入手,优选能够区分岩性和流体的最佳敏感弹性参数,通过叠前反演获得反映储层岩性和含气性的弹性参数体,最后进行地震弹性参数交汇解释预测储层和识别流体.将纵向分辨率较高的测井岩石物理参数分析和横向分辨率较高的地震叠前反演结合起来,应用多个弹性参数,明显提高了储层预测的精度.地震数据中的AVO信息得到了充分的应用,采用纵波独立变量求解Zoeppritz方程.为有效储层预测和流体识别探索出了新途径.在靶区利用叠前弹性纵横波速比的属性,刻画出五级层序单砂体,在此基础上确定建议的水平井井位及其轨迹,通过叠前弹性反演比P波数据独立反演钻井成功率提高了10%以上.     

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.     

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

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

Multi-parameter nonlinear inversion with exact reflection coefficient equation

Amplitude variation with amplitude or angle (AVO/AVA) inversion has been widely utilized in exploration geophysics to estimate the formation of elastic parameters underground. However, conventional AVO/AVA inversion approaches are based on different approximate equations of Zoeppritz equations under various hypotheses, such as limited incident angles or weak property contrast, which reduces their prediction precision theoretically. This study combines the exact P-wave Zoeppritz equation with a nonlinear direct inversion algorithm to estimate the six parameters imbedded in the exact equation simultaneously. A more direct and explicit expression of the Zoeppritz equation is discussed in the case of P-wave exploration, under which condition the incident longitudinal wave produces the reflected longitudinal (P–P) wave and upgoing converted shear (P–SV) wave. Utilizing this equation as the forward solver, a nonlinear direct inversion method is introduced to implement the direct inversion of the six parameters including P-wave velocities, S-wave velocities, and densities in the upper and lower media around an interface, respectively. This nonlinear algorithm is able to estimate the inverse of the nonlinear function in terms of model parameters directly rather than in a conventional optimization way. Model tests illustrate that the nonlinear direct inversion method shows great potential to estimate multiple parameters with the exact Zoeppritz equation.     

An Accurate Method of Computing the Gradient of Seismic Wave Reflection Coefficients (SWRCs) for the Inversion of Stratum Parameters

The optimization inversion method based on derivatives is an important inversion technique in seismic data processing, where the key problem is how to compute the Jacobian matrix. The computational precision of the Jacobian matrix directly influences the success of the optimization inversion method. Currently, most of the AVO (amplitude versus offset) inversions are based on approximate expressions for the Zoeppritz equations to obtain the derivatives of the seismic wave reflection coefficients (SWRCs) with respect to the stratum parameters. As a result, the computational precision and range of applications of these AVO inversions are restricted undesirably. In order to improve the computational precision and to extend the range of applications of AVO inversions, the partial derivative equations of the Zoeppritz equations are established, with respect to the ratios of wave velocities and medium densities. By solving the partial derivative equations of the Zoeppritz equations accurately, we obtained the partial derivative of SWRCs with respect to the ratios of seismic wave velocities and medium densities. With the help of the chain rule for derivatives, the gradient of the SWRCs can be accurately computed. To better understand the behavior of the gradient of the SWRCs, we plotted the partial derivative curves of the SWRCs, analyzed the characteristics of these curves, and gained some new insight into the derivatives. Because only a linear system of equations is solved in our method without adding any new restrictions, the new computational method has both high precision and a quick running speed; it is not only suitable for small incident angles and weak reflection seismic waves but also for large incident angles and strong reflection seismic waves. With the theoretical foundations established in the article, we can further study inversion problems for layered stratum structures and we can further improve the computational speed and precision of the inversions.     

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

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

HTI介质中的反射纵波方位属性

利用横波响应进行裂缝性各向异性介质的检测在实际应用中取得了很好的效果,但技术复杂、成本较高使该方法的广泛使用受到限制,而纵波资料采集和处理技术的精细有效保持了纵波的各种属性,这为直接利用纵波资料进行裂缝检测创造了条件.地下垂直定向裂缝通常用HTI介质模型来描述,为此,本文利用射线追踪和反射率法计算了层状各向同性介质背景下的HTI介质顶、底界面的反射纵波旅行时和反射系数,并分析了这些属性随观测方位的变化规律.研究表明,HTI介质底界面反射纵波旅行时和HTI介质顶界面反射系数表现出了明显的方位各向异性;旅行时、振幅和AVO梯度属性均在0°观测方位和90°观测方位上存在最大差异,可以用多种属性联合来精确判定裂缝的发育方向.     

P-SV波AVO方法研究进展

为有效地利用AVO(amplitude versus offset)信息来反演岩性参数和预测油气储层，详细介绍了P—SV转换波AVO方法的研究现状和主要进展，并针对几种具有代表性的方法的基本思路、方法特点及参数反演等方面进行了对比和评述。同时，给出了两个不同形式的近似公式，并讨论了利用这两种近似公式进行AVO分析的思路和参数反演的方法，通过对比和讨论可以得出：不同近似的主要目的是为了提高近似精度并体现不同的岩性参数对反射系数的影响、敏感程度及其在参数反演中的意义；根据弹性模量和波速之间的关系，目前所有的PP波和P—SV波反射系数近似公式都可以统一表示成射线参数幂级数的形式，这些结论对于AVO理论研究和参数反演都具有重要意义。     

Joint PP and PS AVO inversion based on Zoeppritz equations

Considering Zoeppritz equations, reflections of PP and PS are only the function of ratios of density and velocity. So the inversion results will be the same if the ratios are the same but values of density, velocities of P-wave and S-wave are different without strict constraint. This paper makes efforts to explore nonlinear simultaneous PP and PS inversion with expectation to reduce the ambiguity of AVO analysis by utilizing the redundancy of multi-component AVO measurements. Accurate estimation of ratio parameters depends on independence of input data. There are only two independent AVO attributes for PP reflectivity (i.e. intercept and gradient) and two for PS reflectivity (i.e. pseudo-intercept and pseudo-gradient or extreme amplitude), respectively. For individual PP and PS inversion, the values of least-squares objective function do not converge around a large neighborhood of chosen true model parameters. Fortunately for joint PP and PS inversion the values of the least-squares objective function show closed contours with single minima. Finally the power function fitting is used to provide a higher precision AVO attributes than traditional polynomial fitting. By using the four independent fitting attributes (two independent attributes for PP and PS respectively), the inversion of four ratio parameters (velocities and densities) would be estimated with less errors than that in traditional method.     

带粒子滤波约束的PP-PS联合反演的稀疏解算法

随着地震勘探目标从构造型油气藏向岩性油气藏的转变,地震勘探难度日益增大,这就要求从地震数据中获得更多可靠且具有明确地质含义的属性信息,并充分利用这些属性信息来对储层的岩性、岩相进行分析.AVO三参数反演能够从振幅随炮检距的变化信息中直接提取纵波速度、横波速度以及密度来估计岩石和流体的性质,进而对储层进行预测.然而,AVO反演本身是一个不适定的问题,加上地震纵波反射系数对横波速度和密度的不敏感,会造成单纯利用纵波地震数据进行反演的结果误差大.随着地震接收和数据处理技术的发展,越来越多的学者对PP-PS联合反演方法进行了研究并在实际资料中得以运用.融合转换横波地震数据的联合反演在一定程度上提高了反演的精度,降低了解的不稳定性.但是在信噪比较低的情况下,联合反演的效果受到了限制.本文从优化理论出发,提出了基于粒子滤波提供先验知识的l₁范数约束极小化问题的稀疏解算法.并将上述方法运用到了不同的模型中,通过比较分析,证实了该方法在不同信噪比资料中的有效性和在信噪比较低情况下的优势.     

Cauchy prior distribution-based AVO elastic parameter estimation via weakly nonlinear waveform inversion

Cauchy priori distribution-based Bayesian AVO reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. For the inversion, the computation of the covariance matrix and regularized terms requires prior estimation of model parameters, which makes the iterative inversion weakly nonlinear. At the same time, the relations among the model parameters are assumed linear. Furthermore, the reflectivities, the results of the inversion, or the elastic parameters with cumulative error recovered by integrating reflectivities are not well suited for detecting hydrocarbons and fuids. In contrast, in Bayesian linear AVO inversion, the elastic parameters can be directly extracted from prestack seismic data without linear assumptions for the model parameters. Considering the advantages of the abovementioned methods, the Bayesian AVO reflectivity inversion process is modified and Cauchy distribution is explored as a prior probability distribution and the time-variant covariance is also considered. Finally, we propose a new method for the weakly nonlinear AVO waveform inversion. Furthermore, the linear assumptions are abandoned and elastic parameters, such as P-wave velocity, S-wave velocity, and density, can be directly recovered from seismic data especially for interfaces with large reflectivities. Numerical analysis demonstrates that all the elastic parameters can be estimated from prestack seismic data even when the signal-to-noise ratio of the seismic data is low.