新生界火成岩发育区地震传播机理研究.docxVIP

经典文档 下载后可复制编辑 新生界火成岩发育区地震传播机理研究 张歆，黄建东，程万里 （中国石油大学（北京），地球物理与信息工程学院） 摘要：本文主要是依据渤海湾地区的地层性质，结合渤海油田新生界火成岩发育实际靶区的地质情况，统计分析井点钻遇火成岩岩性，进行岩石物理分析，建立不同岩性、物性以及流体性质与弹性参数包括纵横波速度、密度等参数的关系，并进行对比分析，建立火成岩地球物理参数响应关系，即火成岩岩石物理模型。主要是利用岩层的物性和纵波速度对横波速度进行预测，并进行Zoeppritz方程近似式精度分析。 一、纵横波求解 根据测井测量，Castagna等人（1985）确定了水饱和状态下速度相对于孔隙度和黏土含量的经验公式。 对于泥岩（主要由黏土和粉砂大小的颗粒组成的碎屑岩硅质岩），他们发现和的关系为： （1-1） 对于砂岩，Thomsen(1986)的水饱和的砂岩和页岩的实验室—的数据。出于比较，数据点上加上了Castagna等人（1985）用这些数据拟合的最小二乘法线性方程： （1-2） 对于有较高速度的玄武岩，通过相关的资料，本文采用以下速度经验公式进行求解： （1-3） 二、Zoeppritz方程 在查阅国内外大量文献，书籍，报告和图件后，了解了 AVO 技术的发展和现状，实际操作了 AVO 正演模型的构建和方法，掌握了 AVO 正演模型的特征分析。 1919 年，Zoeppritz 通过研究在弹性介质中，平面纵波非垂直入射到不同介质的分界面时会产生反射波和透射波，由于位移连续性和应力连续性，在引入反射系数和透射系数之后，得到了入射波的入射角和振幅关系的方程，即Zoeppritz 方程。Zoeppritz 方程奠定了分析弹性波在弹性分界面上的反射振幅与入射角的关系的理论基础，同时在地震勘探中为 AVO 技术提供了理论依据。 通过对精确 Zoeppritz 方程的推导和纵波入射的分析，我们在精确 Zoeppritz 方程中的表达式进行了的化简。并对在纵波入射情况下，不同岩石物理参数对反射系数曲线和透射系数曲线的形状和变化的影响，绘出了图片，进行了对比分析，讨论了精确 Zoeppritz 方程在 AVO 正演模型分析中的适用范围。 完整的Zoeppritz方程全面考虑了平面纵波和横波入射时，在界面两侧产生的纵波和横波的反射能量和透射能量与入射角度之间的关系。 当平面纵波以入射角入射到弹性界面上时，界面的两侧产生：反射纵波，反射横波，透射纵波,透射横波。图中各参数分别为：入射角，纵波反射角，横波反射角，纵波透射角，横波透射角，上层介质的纵波速度，横波速度，密度，下层介质的纵波速度，横波速度，密度 R R = 2 \* ROMAN II VP2VS2 VP1VS1 SR PR PT ST Z X 图一： 各种波质点位移方向示意图 根据层状介质中的Snell定律，在界面上： （2-1） 为射线参数，对于任意目的层的射线参数为： （2-2） 可以得到任意地层的入射角为： （2-3） 假设入射波、反射波和透射波的位移矢量为如下射线级数解形式： （2-4） 其中表示波前类型，设介质中弹性参数、密度和位移分量分别为，，介质中弹性参数、密度和位移分量分别为，。根据作用力等于反作用力的原理，在界面上应力连续，介质质点作用于介质质点上的应力等于介质质点作用于介质质点上的应力，于是有连续边界条件： （2-5） 进行化简可以得到： （2-6） 同时介质和介质组成的弹性体是一个整体，在弹性极限内无断裂，也无滑动，因此在边界上的位移是连续的： （2-7） 已知入射波和四种转换波时，