基于比值均方根的杨氏模量反演方法

杨氏模量和泊松比是表征岩石脆性、评价储层可压性以及描述储层流体特征的重要参数,叠前地震反演是从地震数据中得到此类参数的有效手段.常规的叠前地震反演方法多使用Zeopppritz方程的近似式计算反射系数,近似式的误差影响了反演结果的精度.针对这个问题,推导了比值均方根形式的杨氏模量、泊松比和密度的Zoeppritz方程,并基于广义线性反演构建了杨氏模量和泊松比直接反演方法,有效地避免了使用近似式的局限性,提高了反演的精度.模型测试和实际资料试算结果表明：采用比值均方根形式构建的反演流程稳定,能够从叠前地震数据中获得可信的杨氏模量、泊松比和密度,提供了一种可靠的杨氏模量和泊松比直接反演方法.     

非均匀介质地震AVO响应模拟及分析

地震AVO （amplitude versus offset）技术是一项利用振幅信息研究岩性、检测油气的地震勘探技术。常规方法基于Zoeppritz方程计算模型界面处的反射系数,而实际地下非均匀介质中地震反射特征不仅与入射角度、物性差异有关,还与入射波频率、地层厚度、薄互层结构等因素有关。为此,应用传播矩阵理论充分考虑与这些因素有关的调谐干涉等传播效应,针对实际数据计算非均匀地下介质的高精度合成地震记录,对比Zoeppritz方程、Shuey二项近似方程、Shuey三项近似方程以及传播矩阵算法的模拟效果。研究发现：在小角度入射时Zoeppritz方程、Shuey二项近似方程和Shuey三项近似方程的反射振幅和波形基本一致,大角度入射时Zoeppritz方程与Shuey三项近似方程接近;Zoeppritz算法的模拟结果在小角度入射和浅层情况下与传播矩阵算法差别较小,而在大角度入射和深层情况下与传播矩阵差别较大,说明层间多次波的调谐干涉等传播效应不可忽略。     

介质密度反演偏导矩阵的精确计算方法

实现反演偏导矩阵的计算是基于导数最优化反演方法的关键,然而目前的地震反演几乎都是基于Zoeppritz方程近似实现的,使计算精度和适应范围受到限制.本文利用Zoeppritz方程建立了反射系数对地层介质密度比偏导方程,导出了Zoeppritz方程矩阵元对介质密度比的导数.通过求解偏导方程获得了反射系数对介质密度比偏导数的精确计算(考虑了速度中含介质密度的问题).利用数值算例分析了反射系数对介质密度比偏导数的变化特点.本文采用直接解法求解偏导矩阵方程组,获得了快的计算速度和高的计算精度,为实现地层介质密度反演(包括大角度反演)提供了偏导矩阵的计算方法.     

转换波反射系数近似公式的比较及分析

多分量地震资料解释的首要问题就是多分量的标定和对比,其中最关键的是转换波反射系数的求取及子波的选取,然而Zoeppntz方程在数学表达上极其复杂,在物理意义上也不直观,因此很难得到直接的应用.为此,在详细介绍五种近似公式的基础上,利用六种模型分析了随入射角的变化各公式的适应性,最后得出结论,Aki和Richard的近似公式适应性普遍比较好.     

关于VTI介质中多分量转换波地震数据采集的最佳接收时窗

以P_P波时窗和P_SV波时窗的交叉部分作为VTI介质中P_P波、P_SV波联合采集的最佳时窗,在有效时窗内通过基于VTI介质的射线追踪法确定每个反射点的P_P波、P_SV波射线路径和相角,求解Zoeppritz方程得出P_SV波与P_P波反射系数(或透射系数)的比值,并将该比值的绝对值定义为控制因子G. 以G值作为控制VTI介质中P_P波、P_SV波联合采集观测系统设计的主要参数,通过分析G值确定检波器位置的最佳接收时窗,得出优化VTI介质中多分量转换波地震数据采集的接收时窗的方法.     

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

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

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

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

长偏移距P波AVO反射系数的一种近似

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

节理面透射模型及其隔振性能研究

振动波入射到节理界面上,与软弱断层一样会产生波场的分解,能够起到隔振作用。根据线性变形节理的透反射规律研究发现,隔振后各点的波场可由透射P波与透射S波叠加而成,其能量分配服从Zoeppritz方程,不连续位移造成了透射波的衰减;法向刚度增加透射纵波系数会增加,而切向刚度与透射纵波无关,法向与切向刚度对透射S波均有影响,刚度增加透射S波振幅逐渐降低,逐步逼近于0;由于透射后波场相互干涉,不同波场延时不同,透射合成波与入射波相比,振幅降低,持续时间增加,波形尾段有畸变;总能量系数与节理面的刚度无关,而随入射角增加不断降低。     

联合PP、PS波叠前AVA反演方法

从单一pp波数据中提取出来的弹性参数与岩性以及孔隙内流体之间的联系是不充分的,这就造成了在油藏的解释中包含有大量的不确定性.通过加入PS波数据信息,可以有效降低反演的多解性,提高反演的稳定性.笔者所研究的反演方法是基于精确的Zoeppritz方程,利用广义线性反演方法,通过迭代的思想来反演出纵波速度、横波速度和密度.该方法相对单一pp反演可以有效提高横波速度的估计精度.