弹性固体地震—流体井孔弹性波场有限差分数值模拟.pdfVIP

维普资讯 第 38卷 增刊 地 球 物 理 学 报 Vo]．38，Supp． 【 1995年 ACTA GEOPHYSICA SINICA 1995 一 弹性固体地层～流体井孔弹性波场的 有限差分数值模拟 ‘ s弓|舌· 王克协 (吉林大学物理系 长春 130o23) 摘 要 视 关键词妻泣擅l圭b有限差分·流体声压，一薄互层-．了2 瞧‘‘霎 ， ／／ 一 言 的解决水平分层介质的声测井 问题的方法 ．我们 曾用该方法对井外为任意水平多层 介质声场的首波传播进行了数值模拟 ]，对短源距情况 (2．5m)得到了满 意的结果． 计算 ，使半解析方法的应用受到了一定的限制．这是至今 国内外尚未解决的问题． 国家 自然科学基金委员会、中国科学院、中国石油天然气总公司和太庆石油昔理局联合资腑的谭题． I)东北师范大学理论物理研究所．长春 130024 本文 1994年 10月 6日收到．1995年 9月 27日收到修改蔷 维普资讯 地 球 物 理 学 报 38卷 中的传播规律 ．有限差分方法多用在地震勘探、工程设计等实际问题中，用该方法 研究具有圆柱面边界的声测井问题在国内还未见报道，国外也只有少数工作．StephenE 等 1985年首次将有限差分方法 用于声测井 问题 ，将位移场方程直接化为差分方程； Randall等人应用Virieux 提出的速度一应力有限差分方法．处理了各类地层的声测井 问 题 ，特别在处理非轴对称井孔的全波声测井问题取得了进展． 本文提 出声压一速度有限差分方法，研究水平多层介质包围的流体井孔中声传播的 规律，考察薄层对全波列的影响．第 1步 ，对地层为水平多层流体 (地层 的声波方程近 似， =0)情况，实现了井 内、外声场的有限差分数值模拟．第 2步将该方法发展为适于 计算井外为水平多层弹性固体情况 的声测井 问题．通过将处理复杂 问题的有限差分方法 用于处理简单的均质地层 问题 ．与仅适于处理这类 问题的Fourier变换方法进 比较，有限 差分方法的正确性得到了检验．对短源距的流体一流体薄层情况还与半解析方法进行了 比较．最后用声压一速度有限差分方法研究 了井外含有薄互层 时弹性介质地层的弹性波 的传播特性及其规律 ．得到了有意义的结果．从模拟井外有薄层介质存在时的全波曲线 看，层速度显示正确，并可 以定性判断薄层出现的位置． 二、理 论 公 式 及 算 法 考虑一各向同性脉冲点源位于半径为 的无限长 圆柱形井孔中心．选柱坐标系使坐 标原点与声源所在处重合．。轴与充流体井轴一致．设井 内为均匀各向同性流体， 和 ¨分别为流体密度和声速．井外为均匀或分层均匀的流体介质 (对应 一0)或各向同性 弹性固体介质，其密度和纵、横波速度分别为 、 和Vs(对于流体Vs=0)如图1所示． 由于所考虑的问题具有轴对称性 ，使定解 问题的解与P角无关．体系只与空间坐标 r、 有关．为减少计算量，选 r、 轴为对称轴．有限差分合成声测井体系的几何模型 图 1 井外有水平多层介质的测井声系 图2 有限差分合成声测井的几何模型 维普资讯 增刊 l