超声波测井条件下砂岩的波动方程及其压实条件探讨

本文将纯砂岩－维超声实验模型推广到三维，用声速公式得到了砂岩的超声统计模型－威利公式及其压实条件。导出欠压实砂岩超声波动方程（理论模型）。引入压实条件，得到压实砂岩的超声波动方程。最后，简单讨论了三种超声模型之间的关系，给出了超声波动方程的应用实例。     

新场气田气井产能和无阻流量的确定

利用新场气田的压力恢复试井解释成果，经理论计算求取气井的产能方程和气井的无阻流量，求得Vogel方程的系数α及新场气田计算无阻流量的一点法公式。计算结果与实测产能方程对比表明，该公式比较适合新场气田的实际，计算精度比常用的一点法公式高。     

气井产能方程系数的预测

气井在生产过程中，由于天然气粘度和气体偏差系数的变化，早期测试的产能方程系统需要根据地层压力的变化人修正，本文推导了出二项式产能方程系数Ａ，Ｂ值的预测公式，同时对“一点法”经验的公式也进行了预测。     

裂缝介质的纵波方位AVO反演研究

基于Zoeppritz方程的地震波反射系数解析式是非线性的，而利用地震波反射系数公式来反演介质的弹性参数，或是直接利用非线性反演方法，或是对复杂的反射系数公式进行线性化近似。首先，从理论上讨论了各向异性介质中的反射和透射系数方程，给出了裂缝介质的线性近似反射系数公式；然后，利用ISO／HTI模型对Vavrycuk和Ruger的近似反射系数与真实反射系数进行了对比，得出在大角度情况下Vavrycuk近似比Ruger近似更逼近真实值的结论。因此．选择Vavrycuk纵波近似反射系数公式对纵波方位AVO线性反演公式进行了推导；最后，用模型对反演公式进行了验证，反演结果能满足精度要求。     

气体流动基本方程与Colebrook-White摩擦阻力公式的综合应用

将气体流动基本方程与Colebrook White摩擦阻力公式相结合推导出新的公式 ,大大简化了气体流动基本方程原来需要多次试算、迭代的计算过程 ,并且推导过程未加任何假设。新公式的使用范围和计算精度与原Colebrook White公式相同 ,因此对气体流动基本方程的应用有着重要的改进意义     

油井脱气后脱气半径,地层压力分布及平均地层压力的求法

通过赫氏方程和裘比公式相结合的方法，推导出脱气半径、地层压力分布和平均地层压力的计算公式，并以喇嘛甸油田为例，说明了这些公式的具体应用。     

套管自振分析

本文在动力方程与相容方程的基础上，结合边界条件，采用试函数法求解，得到套管自振圆频率公式，以解决套管的自振问题，这是文献［１］特殊情况。     

辽河盆地潜山油藏高压物性计算方法研究

根据辽河盆地潜山油藏37个原油样品高压物性分析资料统计，分两个凹陷建立了潜山油藏地层油体积系数、原始气油比两项高压物性参数的经验方程。同时应用国外PVT经验公式进行计算、对比，得出结果：用本油田实际分析资料确立的经验公式应用效果要好于国外经验公式，证明该公式是可行的。     

弹性波介质中纵,横波速度反演

波动方程反演的直接目的是求得纵、横波速度。本文从弹性波动方程出发，利用分布参数控制系统的伴随理论，导出了纵、横波速度的梯度公式，然后采用优化迭代技术进行数值反演，模型试算获得了较好的结果，表明该方法是可行的。     

油气水三相流动时油井含水率预测方法研究

为研究溶解气驱油藏油井含水率的变化规律，对油气三相流动时，产油量与流压、产水量与流压的关系进行了研究，利用Ｖｏｇｅｌ方程和产水量的经验公式，推导出了含水率率预测的关系表达式。经实测数据验证，表明该公式计算简单，精度高，适应性广。     

初至波射线层析成像在复杂区静校正中的应用

在复杂地表和复杂构造地区，静校正问题突出。为此，根据目前静校正存在的问题和实际需要，进行了地震初至波射线层析成像方法的研究。给出了初至旅行时层析成像方法的基本原理，包括近地表模型的建立、射线路径和初至时间的计算、反演方程的建立以及共轭梯度反演方法等。其中提出的对复杂介质采用动态网络最短路径射线追踪方法，克服了同类方法得到的射线路径呈之字形、计算出的旅行时比实际旅行时系统偏大的缺陷。利用非线性共轭梯度优化方法来实现初至旅行时层析反演，反演过程稳定收敛。利用理论模型和实际资料对方法进行了验证，结果表明，采用上述方法求取复杂地区的静校正量，静校正效果很好，剖面的品质得到了较大的改善。     

射线参数域叠前弹性参数反演在储层预测中的应用

P道集的地震波是按照射线路径传播的,与角度道集相比,更接近地震波的实际传播路径,更能表征地下的反射特征。同时P道集可直接从地震数据计算得到,不依赖于上覆地层速度,避免了不精确的速度模型给反演带来的误差。文章在介绍P道集基本概念及射线参数域叠前弹性参数反演流程的基础上,利用分P叠加数据,通过贝叶斯稀疏脉冲反演进行了不同P值的射线弹性阻抗反演,并通过约束搜索法求取了相关的叠前弹性参数。相比于角度域的叠前弹性参数反演,其反演结果具有高分辨率、高精度且客观真实的特点。将该方法应用于川西坳陷新马滚动区以及川北坳陷元坝西区的储层预测,实际应用效果证明了该方法的可靠性和准确性,其预测结果对研究区下一步的勘探开发工作具有重要的指导意义。     

射线波阻抗反演在YA地区储层预测中的应用

本文利用射线波阻抗对苏北盆地YA地区砂泥岩薄互储层开展研究。应用结果表明射线波阻抗与声波阻抗交会能有效区分薄互储层,两者与储层物性参数的三元交会能有效判断储层物性优劣。通过YA地区的射线波阻抗反演结果准确识别了古近系始新统戴南组一段目的层两个主物源方向,精细刻画了有利砂岩展布形态并识别出砂岩上倾尖灭型圈闭,获得了良好的应用效果。     

射线层析与波动方程联合成像

射线层析与波动方程联合成像是构造和速度同时成像的新思路。由偏移共成像道集算出深度偏差,进而算出速度修正量,把速度修正量按射线路经均匀地分配,由浅入深,逐层修正速度模型;再用修正后的速度作波动方程偏移,抽取共成像道集,算出速度修正量,再将其按射线路经均匀地分配。如此反复进行波动方程偏移和射线法反演迭代,直到共成像道集深度方差满足预设的精度为止,最后得到速度和深度构造同时成像。仿真实验验证了本方法。     

综合静校正技术及其在川东高陡构造区应用效果

不同的静校正技术基于不同的模型假设,适用于不同的地表地质条件。川东高陡构造地区,地震资料处理效果证明,在老地层出露区,微测井约束的层析静校正方法取得了满意的效果;在地表平缓地段,基于EGRM的折射静校正效果最好。只有采用综合静校正技术才能优势互补,解决好一次静校正问题。然后进行初至波多域剩余静校正和反射波剩余静校正。讨论了折射静校正、层析静校正和初至波多域剩余静校正的原理、假设条件、适用范围和误差分析,结合实际资料,处理中利用野外高密度的微测井资料,获得折射静校正的风化层速度;利用折射静校正反演的速度模型和微测井建立的速度模型,建立层析成像初始速度模型;利用直达波、折射波和回折波,进行初至波射线路径追踪、走时计算和大型稀疏矩阵方程的求解。     

关于深度偏移的几个相关概念

本文仅对深度偏移和时间偏移的关系，叠前深度偏移和叠后深度偏移的关系．法向射线与成像射线、射线偏移、相干反演等相关概念进行了分析。明确了所有的偏移公式都是在深度域定义的；只有当前层以上所有各层的速度都正确时，修改速度模型使其共成像点道集上同相轴达到平直，方可认为该速度和深度接近于真实值。论述了射线偏移在建立速度一深度模型过程中的重要作用，它可以把时间域层位转化为深度域层位；进一步明确了法向射线与成像射线的概念．利用法向射线可以将叠加时间剖面上层位转换为深度域层位；利用成像射线可以把时间偏移剖面的层位转化为深度域层位；提出了利用相干反演建立初始速度深度模型的优缺点及克服其缺点的办法。只有深入地理解这些概念才能提高利用深度偏移方法解决问题的能力。     

四面体网格剖分下速度与反射界面的同时反演

针对含起伏地表、不规则反射界面及速度异常体的复杂模型的走时正、反演问题,本文采用四面体单元进行模型参数化,实现多震相地震射线的追踪计算,且推导了四面体单元下走时关于速度、走时关于反射点深度的偏导数计算公式;结合共轭梯度法求解带约束的阻尼最小二乘反演问题,实现了多震相走时联合及同时反演成像;通过在同时反演中对不同种类Jacobi偏导元素做归一化处理,有效解决了不规则界面存在散射而导致的散射点所在区域射线过密、速度和界面过度更新问题。数值模拟结果表明：采用四面体单元进行模型参数化,复杂模型正演计算结果具有较高精度;同时反演可准确刻画不规则反射界面和速度异常体。因此,本文提出的走时成像方法具有有效性和实用性。     

一种任意界面的逐段迭代射线追踪方法

本文根据同一条射线路径上满足相同射线参数的原理，首先采用一个适合于任意界面情况下计算反（折）射点的一阶近似公式，从任意给定的初始路径出发进行逐段迭代计算，求出界面上的反（折）射点。然后以求出的射线路径为基础，按上述方法逐段求取下一界面的反（折）射点，直到终点得到整条射线为止。当整条射线路径上的界面反（折）射点的修正值达到精度要求时，则认为射线追踪结束。通过模型计算表明，本方法的计算速度比“打靶法”提高了一个数量级，其精度可以根据需要来确定。本方法无射线盲区，为任意界面的地震层析成像提供了一种准确、快速的方法。     

复杂构造对下伏地层地震波动力学特征的影响

复杂构造下伏地层的CRP(共反射点)射线路径和振幅变化特征,对后续地震资料的研究如叠前反演、AVO(振幅-炮检距关系)特性分析至关重要。采用射线追踪法进行地震正演模拟,将过邻近井地震偏移剖面解释成果应用于复杂构造模型的设计,对复杂构造下伏地层地震波运动学和动力学特征的影响程度进行分析。分析表明,常规的复杂构造勘探方法所获得的CRP道集,其地表接收点关于CRP非对称性严重,陡坡带底界面上共反射点覆盖次数和反射能量横向变化不均匀性严重,仅在局部地段达到满覆盖次数,复杂构造对下伏地层反射动力学特征影响较大。