声反射成像测井在地层中的三维波场模拟方法研究.pdfVIP

中国科学 D 辑 :地球科学 2008 年 第 38 卷 增刊 Ⅰ: 166 ~ 173 SCIENCE IN CHINA PRESS 声反射成像测井在地层中的三维波场模拟 方法研究 * 陶果, 何峰江, 王兵, 王华, 陈萍 油气资源与探测国家重点实验室( 中国石油大学(北京)), 北京 102249 * Email: taoguo@ 收稿日期: 2007-04-20; 接受日期: 2008-01-30 国家自然科学基金项目(批准号:和国家重点基础研究发展计划项目(编号: 2007CB209601)资助 摘要 声反射成像测井方法和仪器近来成为勘探地球物理界的新研究热点. 但有利于 关键词 成像的反射波信号与井中的模式波叠加而难以分离出来. 为了改进仪器设计, 研究更有 反射声波成像 效的数据处理和解释方法, 必须用以波动方程为基础的数值模拟方法来精确描述声反射 三维波场模拟 成像测井的波场. 为此建立了非均匀交错网格、加PML 吸收边界的、能够用于正交各向 非均匀交错网格 PML 吸收边界 异性介质的 3 维有限差分数值模拟方法和并行算法程序. 将其应用于典型各向同性和各 并行算法 向异性地层中反射声波成像测井的三维波场模拟, 并与已公开发表的应用离散波数方法 的相应结果对比, 证明了此方法正确、算法高效. 模拟的波场动态范围能够大于 60 dB, 计算速度比常规算法提高了 3.5 倍. 这些技术能够正确、有效地模拟大尺度(在 10 m 左 右) 、较高频率(10 kHz 左右) 的反射声波成像测井波场. 为我们研究声反射成像测井仪器 和资料处理解释方法提供了有力的定量分析工具. 近年来, 声反射成像测井方法和仪器成为勘探地 位于分界面上倾(或下倾)一侧时, 会接收到P-P 波反 球物理界的新研究热点之一. 国际上, Schlumberger公 射信号, 如图 2 所示; 当发射器和接收器分别置于分 司继数年前研制出的BARS(Borehole Acoustic 界面的上、下倾一侧时, 会接收到P-S 转换波信号和 Reflective Survey)[1]仪器之后, 2006 年又推出了含有 S-P转换波信号, 如图 3 所示. 实际上在图2 所示情形 [ 2] 改进的声反射成像功能的新仪器SonicScanner , 下, 还可能接收到P-S , S-P 以及S-S转换波. 用这些井 Baker Hughse 公司开发出了基于XMACII仪器的声 附近的反射波信号成像能够获得有关构造位置, 几 反射成像处理软件 [ 3] . 我国大港油田测井公司于 何形态, 倾角, 方位等重要信息 [4,5]. 然而, 目前这项 2005 年开发了声反射成像功能的新仪器 1). 这些仪器 重要的应用往往遇到很大困难, 因为有利于成像的 在井内测量时其声源辐射到地层中的能量以P 波和S 反射波信号被淹没在沿井壁传播的直达波信息中从 波的形式向外传播. 当遇到声阻抗不连续界面时会 而变得模糊不清, 这些反射波与井中的模式波叠加, 发生反射和折射, 如图 1 所示, 并伴随着P波和S波之 难以分离出来. 以往采用几何声学的方法来分析反 间的相互转换. 由Snell定律, 当发射和接收换能器都 射信号的信