海上粘弹介质FCT有限差分正演模拟.pdfVIP

石油天然气学报 (江汉石油学院学报) 2011年12月 第33卷 第12期 JournalofOilandGasTechnology(J．JPI) Dec．2011Vo1．33 No．12 海上粘弹介质 FCT有限差分正演模拟 ‘成景旺，顾汉明，刘 琳 (中国地质大学 (武汉)地球物理与空间信息学院，湖北武汉430074) 刘春成 ，刘志斌 (中海油研究总院，北京100027) [摘要]高阶有限差分正演模拟将时间和空间离散，以差分代替微分，其求解精度受子波主频、网格大小 以及有限差分精度 的影响，不当的话会 引起数值频散。海上勘探采用 的气枪震源子波频带较宽，主频高， 波长短，数值频散现象严重，而减小空间网格大小和提高有限差分精度需要较大的计算量和较长的计算 时间。针对这一特点，采用通量校正传输 (FCT)技术与有限差分相结合的方法来减少数值频散。在分 析气枪震源子波的特征及该子波在粘弹介质中传播频谱特征的基础上，呆用FCT有限差分模拟对粘弹介 质盐丘模型进行 了正演模拟 。结果表 明，FCT有限差分能够有效地压制数值频散 ，可 以提高模拟 的精度 和运算效率。 [关键词]粘弹介质；有限差分；数值频散；通量校正传输；海洋地震勘探 [中图分类号]P631．44 [文献标识码]A [文章编号]1000—9752(2Ol1)12—0083—05 海上石油天然气等矿产资源的巨大潜力使得海上地震勘探开发和研究变得极其重要 。而正演数值模 拟技术在地震勘探的各个阶段都有重要的作用。交错网格高阶有限差分法内存 占用少 、模拟精度高以及 计算快使得其成为 目前最为流行 的方法之一 1『]。有限差分将时间和空间离散，以差分代替微分 ，在求 解过程中会产生不期望的数值频散现象 。数值频散的严重程度与子波的主频、空间网格大小以及有限差 分格式精度有着很大的关系。陆地勘探采用的爆炸震源 ，由于近地表品质因子较小 ，吸收 比较严重 ，穿 过近地表后子波主频迅速减小 。而海上气枪激发是在海水里进行 ，海水 的品质因子较大，对子波 的吸 收几乎没有 。故气枪震源子波频率较高 ，最高甚至达到 200Hz左右 ，每一个波长 内的采样点很少，数 值频散严重，虽然减小空间网格大小或提高有限差分精度可 以得到改善，但是将会大大增加计算量 。 Boris等 在流体力学连续方程的求解过程中，发展 了一种通量校正传输 (FCT)算法 ，将其应用于声 波方程的求解中，有效地压制了在粗 网格下差分计算产生的数值频散 、侣]。而且在实际的海上地震勘探 中，当地震波在粘弹介质中传播时，介质粘滞性引起 的吸收作用是导致地震波衰减的主要因素之一 ，所 以进行粘弹介质的正演模拟是必不可少的，粘弹介质中的数值模拟更接近于实际 【]，基于小波域提取 地层品质因子的方法为粘弹介质的正演模拟提供 了前提条件 。笔者从粘弹性理论人手研究波 的吸收 和衰减，推导了波尔茨曼介质的一阶应力一速度粘弹性波动方程 ，将有限差分与FCT相结合，实现了适 合海上勘探的粘弹介质波动方程的数值模拟。 1 海上气枪震源子波的特征 在当前的海洋地震勘探 中，为了获得较好的气枪子波特性 ，提高震源能量，都采用气枪阵列组合 。 其中，空气枪采用调谐组合或相干组合 ，都可以实现主脉冲的加强及对气泡脉冲和以后的气泡脉冲的干 扰的压制，以满足高能量、高气泡比的海上震源要求 。常常通过不同参数的试验 ，得到最佳激发参数下 的最佳气枪子波及其频谱 (图 1)。 由图1可以看到，气枪子波是一个近似的最小相位子波，其延续时间较短，频谱成分较宽，有一个 [收稿 日期]2011—10—10 [基金项 目]国家科技重大专项 (2011ZX05025一oo1)。 [作者简介]成景旺 (1987一)，男，2009年大学毕业，博士生，现主要从事地震正演模拟与数据处理的研究工作。 石油天然气学报 (江汉石油学院学报) 2011年 12月 薹 管 运 恕 罂 鞲