第４０卷第３期 河海大学学报（自然科学版） Ｖｏ１．４０ Ｎｏ．３ ２０１２年５月 Ｊｏｕｒｎａｌ ｏｆ Ｈｏｈａｉ Ｕｎｉｖｅｒｓｉｔｙ（Ｎａｔｕｒａｌ Ｓｃｉｅｎｃｅｓ） Ｍａｙ ２０１２ ＤＯＩ：１０．３８７６／ｊ．ｉｓｓｎ．１０００—１９８０．２０１２．０３．０１７ 碳酸盐岩生物礁地震响应特征的数值模拟 但志伟 ，肖 为 ，方中于 ，曾 驿２，汪瑞良２，张宏兵。，尚作萍。 （１．中国海洋石油总公司能源发展股份有限公司钻采工程研究院地球物理研究所，广东湛江５２４０００； ２．中海石油（中国）有限公司深圳分公司研究院，广东广州５１０２４０；３．河海大学地球科学与工程学院，江苏南京２１００９８） 摘要：针对碳酸盐岩生物礁储层具有较强非均质性的特点，利用叠前黏弹性波动方程正演数值模 拟方法，模拟了碳酸盐岩生物礁、滩的地震响应特征，为地震资料识别生物礁储层提供依据．模拟结 果显示，在生物礁两侧出现上超反射，礁顶明显见披覆现象，在礁体内部产生ｓ形前积反射、弱反射 或杂乱反射，在礁前缘斜坡相外形呈楔状且礁底为下超反射终止．此外，礁体顶底界反射产生极性 反转，底部反射同相轴可能出现上凸或下凹现象，而深度域成像剖面可以消除上凸或下凹现象． 关键词：碳酸盐岩；生物礁；叠前地震；黏弹性波动方程；地震响应 中图分类号：Ｐ６３１．４ 文献标志码：Ａ 文章编号：１０００—１９８０（２０１２）０３—０３４５—０５ Ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ ｏｆ ｓｅｉｓｍｉｃ ｒｅｓｐｏｎｓｅ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｃａｒｂｏｎａｔｅ ｒｅｅｆ ＤＡＮ Ｚｈｉ－ｗｅｉ ，ＸＩＡＯ Ｗｅｉ ，ＦＡＮＧ Ｚｈｏｎｇ．ｙｕ ，ＺＥＮＧ ｙｉｚ，ＷＡＮＧ Ｒｕｉ－ｌｉａｎ ， ＺＨＡＮＧ Ｈｏｎｇ－ｂｉｎｇ３，ＳＨＡＮＧ Ｚｕｏ－ｐｉＩｌ （１．Ｇｅｏｓｃｉｅｎｃｅｓ Ｉｎｓｔｉｔｕｔｅ ｏｆ ＣＮＯＯＣ Ｌｔｄ．，Ｚｈａｎｆｉａｎｇ ５２４０（０，Ｃｈｉｎａ； ２．Ｓｈｅｎｚｈｅｎ Ｂｒａｎｃｈ ｏｆ ＣＮＯＯＣ Ｌｔｄ．，Ｇｕａｎｇｚｈｏｕ ５１０２４０，Ｃｈｉａｎ； ３．Ｃｏｌｌｅｇｅ ｏｆ Ｓｃｉｅｎｃｅｓ ａｎｄ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｈｏｈａｉ Ｕｎｉｖｅｍｉｔｙ，Ｎａｎｆｉｎｇ ２１００９８，Ｃｈｉａｎ） Ａｂｓｔｒａｃｔ：Ｂａｓｅｄ ｏｎ ｔｈｅ ｎｏｎｈｏｍｏｇｅｎｅｏｕｓ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｃａｒｂｏｎａｔｅ ｒｅｅｆ ｒｅｓｅｒｖｏｉｒｓ，ｓｅｉｓｍｉｃ ｆｏｒｗａｒｄ ｍｏｄｅｌｉｎｇ ｗａｓ ｃｏｎｄｕｃｔｅｄ ｔｏ ｓｉｍｕｌａｔｅ ｔｈｅ ｓｅｉｓｍｉｃ ｒｅｓｐｏｎｓｅｓ ｏｆ ｃａｒｂｏｎａｔｅ ｒｅｅｆｓ ａｎｄ ｂｅａｃｈｅｓ ｕｓｉｎｇ ｈｔｅ ｐｒｅ—ｓｔａｃｋ ｖｉｓｃｏｅｌｓａｔｉｃ ｗａｖｅ ｅｑｕａｔｉｏｎ， ｉｎ ｏｒｄｅｒ ｔｏ ｐｒｏｖｉｄｅ ａ ｂａｓｉｓ ｆｏｒ ｉｄｅｎｔｉｉｆｃａｔｉｏｎ ｏｆ ｃａｒｂｏｎａｔｅ ｒｅｅｆ ｒｅｓｅｒｖｏｉｓｒ ａｃｃｏｒｄｉｎｇ ｔｏ ｓｅｉｓｍｉｃ ｄａｔａ．Ｔｈｅ ｓｉｍｕｌａｔｉｏｎ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｔｈｅ ｓｕｐｅｒｐｏｓｅｄ ｏｖｅｒｔａｋｉｎｇ ｒｅｆｌｅｃｔｉｏｎ ＯＣＣＵｌＴｅｄ ｏｎ ｔｗｏ ｓｉｄｅｓ ｏｆ ｔｈｅ ｒｅｅｆ ａｎｄ ｔｈｅ ｄｒａｐｅ ｃｏｖｅｔｉｎｇ ｐｈｅｎｏｍｅｎｏｎ ａｐｐｅａｒｅｄ ｏｖｅｒ ｈｔｅ ｒｅｅｆ．Ｓ—ｓｈａｐｅｄ ｆｏｒｗａｒｄ ｓｅｄｉｍｅｎｔａｒｙ．ｗｅａｋ．ａｎｄ ｍｕｓｓｙ ｒｅｆｌｅｃｔｉｏｎｓ ｔｏｏｋ ｐｌａｃｅ ｉｎｓｉｄｅ ｔｈｅ ｒｅｅｆ．Ｔｈｅ ｏｕｔｅｒ ｓｈａｐｅ ｏｆ ｔｈｅ ｒａｍｐ ｆａｃｉｅｓ ｉｎ ｔｈｅ ｒｅｅｆ ｆｒｏｎｔ—ｅｄｇｅ ａｐｐｅａｒｅｄ ｗｅｄｇｙ ａｎｄ ｔｈｅ ｏｖｅｒｔａｋｉｎｇ ｒｅｆｌｅｃｔｉｏｎ ｓｔｏｐｐｅｄ ａｔ ｔｈｅ ｂｏｔｔｏｍ． Ｍｏｒｅｏｖｅ￣ｒｅｖｅｒｓｅ ｐｏｌａｒｉｔｙ ｏｃｃｕｒｒｅｄ ｂｅｔｗｅｅｎ ｔｈｅ ｔｏｐ ａｎｄ ｂｏｔｔｏｍ ｒｅｆｌｅｃｔｉｏｎｓ ｉｎｓｉｄｅ ｔｈｅ ｒｅｅｆ．Ｔｈｅ ｂｏｔｔｏｍ ｒｅｆｌｅｃｔｉｏｎ ｐｒｅｓｅｎｔｅｄ ｕｐｗａｒｄ－ｃｏｎｖｅｘ ｏｒ ｄｏｗｎｗａｒｄ—ｃｏｎｃａｖｅ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ．Ｔｈｅｓｅ ｃｏｕｌｄ ｂｅ ｅｌｉｉｍｎａｔｅｄ ｉｎ ｔｈｅ ｄｅｐｔｈ ｉｍａｇｅ ｐｒｏｆｉｌｅ． Ｋｅｙ ｗｏｒｄｓ：ｃａｒｂｏｎａｔｅ ｒｅｃｋ；ｏｒｇａｎｉｃ ｒｅｅｆ；ｐｒｅ—ｓｔａｃｋ ｓｅｉｓｍｉｃ；ｖｉｓｃｏｅｌｓａｔｉｃ ｗａｖｅ ｅｑｕａｔｉｏｎ；ｓｅｉｓｍｉｃ ｒｅｓｐｏｎｓｅ 波场数值模拟技术是研究各种地震地质条件下构造、物性和岩性等各种地质因素与地震波响应特征之 间关系的一门技术，可以提高对复杂介质中地震波传播规律的认知，检验地震资料处理及解释技术的应用效 果及适用范围【１－３ｊ．所以波场数值模拟技术一直是地震勘探领域的重要研究内容． 正演数值模拟技术分为２个重要分支，一是求解偏微分方程的波动方程数值模拟，二是由积分方程以求 解波场传播旅行时为主的射线追踪数值模拟［４－９］．射线追踪数值模拟方法出现较早，在模拟相对简单的构造 模型时陕速高效，且提供的射线路径非常直观．射线追踪数值模拟方法又分为程函方程法、射线法和弯曲射 线法等，国内外大部分地球物理软件公司早期推出的数值模拟软件大多采用射线追踪数值模拟技术．此外， 出现了一些改进的射线追踪数值模拟方法，如由传统的射线追踪数值模拟方法改进的Ｇａｕｓｓ束法、射线法与 Ｋｉｍｈｈｏｆｆ近似法相结合的快速算法、采用有限差分技术求解格林函数方法等． 收稿日期：２０１１一Ｏ１—１２ 基金项目：国家科技重大专项（２００８ＺＸ０５０２５） 作者简介：但志伟（１９８ｏ一），男，湖北赤壁人，工程师，主要从事石油物探研究．Ｅ－ｍａｉｌ：ｄａｎｚｈｗ＠（３１１＂ｏｏ（３．ＣＯｉｌ３．ｃｎ 河海大学学报（自然科学版） 第４ｏ卷 射线追踪数值模拟方法很难满足高精度模拟复杂地质模型地震波场动力学特征的要求，而波动方程数 值模拟提供的波场快照图有助于细致分析地震波场在各种复杂介质内部的运动学和动力学特征．根据不同 的模拟方程，波动方程数值模拟主要分为声波模拟、弹性波模拟、黏弹性波模拟以及裂隙和孔隙弹性的Ｂｉｏｔ 模拟等ｌ１０－１４Ｊ．波动方程模拟有很多数值实现方法，如有限元、有限差分、Ｆｏｕｒｉｅｒ变换、边界元、谱元、广义屏 等，其中有限差分法是目前研究最多的波动方程数值模拟方法．为了克服频散和运算速度较慢等问题，有限 差分已从早期的低阶发展到高阶，网格划分也从早期的规则网格发展到交错网格，直至可变网格．为此，本文 基于交错网格有限差分技术的黏弹性波动方程，对几种典型的生物礁、滩地质模型进行模拟，分析不同生物 礁、滩地质模型的地震反射波组特征，为根据地震资料有效识别生物礁、滩提供理论依据． ．ｒＬ ．Ｉ一．．， １地震正演数值模拟方法 ＋ １ 一 ｒ １．１黏弹性波动方程 由地震波在介质中传播的本构方程和动量守恒方程可得三维黏弹性介质中的一阶速度一应力关系ＮＨ３ Ｊ为 ＋ √ 一 ＋ ； ￡一 ａ￡ 一 ＋ ｒ ＋ ｒ ～ 式中：ｘｉ， ——笛卡儿坐标； ——时间；仃 ——第 个应力张量； ｆ—— 方向的速度；Ｍ， ——限定参数 （限定在参考频率∞ｏ处纵、横波的相速度）；ｒ——优化参数；下角标Ｐ和ｓ分别代表纵波和横波．ｒ 是一个记 一巩 ≠ 忆变量，可以表示为 ． 一 ａ ｒｇ ａｔ一 一１［（腑 ＋ ２ ｒ ㈠ 划 Ｏ，ｒｉ＋ ８ｘｉ ｌ＋ｒｉｊ］ ｏ３Ｖ（２） （３）Ｌｊ 式中 为优化参数．此外，动量守恒方程为 ｌ０ ｌ０ ｉ一＋＾ ＋ ＝ Ｏｘｉｊ式中：１０——密度； ．—— 方向的外力． １．２数值模拟步骤 针对黏弹性波动方程式（１），采用高阶交错网格差分技术进行离散处理，本文使用四阶中心差分来近似 空间微分，二阶中心差分来近似时间微分． 此外，常规处理和速度分析基于Ｆｏｃｕｓ处理系统，采用ＧｅｏＤｅｐｔｈ叠前深度偏移等模块进行叠前偏移处 理．正演模拟步骤如下： ａ．采用人机交互的方法建立地震地质模型． ｂ．进行地震地质模型网格离散化，离散网格的大小考虑背景速度、异常体最小尺度和收敛条件． ｃ．进行常规地震资料处理，获取水平叠加时间剖面、时间／深度偏移剖面和保幅ＣＲＰ道集． 本文正演数值模拟采用与野外采集接近的观测系统和采集参数，其中道间距２５ ｍ，炮间距２５ ｍ，最小偏 移距Ｏｍ，最大偏移距３０００ｍ，接收道数１２８道，采样间隔１ ｍｓ，记录时间２．５ ｓ，雷克子波的主频分别为３０ Ｈｚ 和６０Ｈｚ． ２生物礁地质模型的地震反射响应特征 从测井资料可以看出，即使在沉积环境不同，或成岩作用也不尽相同的情况下，碳酸盐岩都具有共同的 特点：传播速度比较快，密度比较大，明显高于其上覆和下伏以砂泥岩为主的碎屑岩．碳酸盐岩Ｉ 门地震反射响 应表现为强反射特征，碳酸盐岩顶界反射以强振幅波峰出现，而底界则往往以强振幅的波谷形式出现，这为 在地震反射波记录上识别碳酸盐岩奠定了基础．碳酸盐岩因沉积环境的不同，分为生物礁相碳酸盐岩、生物 滩相碳酸盐岩和泻湖相碳酸盐岩，尤以生物礁最为有利．典型生物礁的地震反射特征主要表现为：具有强反 第３期 但志伟，等碳酸盐岩生物礁地震响应特征的数值模拟 ３４７ 射外形、礁体内部杂乱反射、礁滩翼部周围有上超反射及前积现象、礁滩上覆地层存在披覆构造、礁体底部出 现上凸或下凹现象等．生物礁在时间域成像地震剖面上的这些反射特征，可以帮助发现生物礁礁体． 因此，建立不同形态，不同沉积相碳酸盐岩生物礁的地质框架模型，通过模型正演，分析不同类型碳酸盐 岩生物礁沉积模型的地震 际地震剖面中提取的结果． ２．１塔礁概念模型及地震响应特征 特征，如振幅、频率、相位等方面的特征．本文涉及的几种类型生物礁都是从实 图１为一个塔礁地质概念模型，由２个礁体上下叠置组合而成（孔隙均不发育）．塔礁是在局部轻度隆起 的碎屑岩上生长发育而成，最终被海平面快速上升而淹死，其围岩均为碎屑岩．该生物礁具有高速（２８５０— ５２００ ｒｎ／ｓ）、高密度（２．５ ２．５７ ｃｍ３）的特征，围岩为低速（２３００ ２８６０ｍ／ｓ）、低密度（２．１８—２．３３ｇ／ｃｍ３）．模型 横向长度为１２ｋｒｎ，纵向深度为２５００ｍ，中间２层为２个不同期次重叠的礁体，设计有隆起构造．数值模拟采 用的观测方式为基于射线的零偏移距正演． 图２（ａ）为模型正演模拟时间域成像地震剖面．使用 Ｒｉｃｋｅｒ子波，其主频为６０ ＨＺ．从图２可以见到明显的地震 反射波组特征：礁体与围岩之问为强反射，礁体内部反射 非常弱；礁体顶界为正反射，礁体底界为负反射，这种极 性反转有助于识别礁体．此外，礁体中央部位呈现低频， 向边缘频率逐渐升高．在生物礁的外围上超现象明显，生 物礁顶上披覆现象明显，礁体下部及下伏界面反射呈现 明显上凸特征．利用这些振幅反射特征可以区分碳酸盐 岩和碎屑岩．图２（ｂ）为使用相同子波模拟深度域成像剖 面．对比深度域与时间域成像剖面可以发现，在深度域成 图１塔礁地质概念模型的参数模型 像剖面中，礁体下伏反射层上凸现象得到消除，与设计模 Ｆｉｇ．１ Ｐａｒａｍｅｔｅｒ ｍｏｄｅｌ ｏｆ ｇｅｏｌｏ￣ｃａｌ ｃｏｎｃｅｐｔｕａｌ ｍｏｄｅｌ 型一致．说明深度域成像效果优于时间域成像，为此对于 ｆｏｒ ａ ｔｏｗｅｒ ｒｅｅｆ 碳酸盐岩生物礁储集层要求使用深度域成像结果． ＝＝＝ ＝＝＝ 乎…… — 善 。 ＿一 ＿＿・ ＝ 越９９７ …～ ～ 一一一 ＝ 二 一＿一—～～ ６Ｏ １ ４９７ ｌ ｌ０ｌ ２０ｌ ３０１ ４０１ ５０１ ６０１ ７０１ ８０１ ９０１ １ １Ｏ１ ２０１ ３Ｏ１ ４Ｏｌ ５０ｌ ６０１ ７Ｏ１ ８０１ ９Ｏ１ ＣＤＰ ＣＤＰ （ａ）时问域成像 （ｂ）深度域成像 图２模拟主频为６０Ｈｚ的地震剖面 Ｆｉｇ．２ Ｓｉｍｕｌａｔｅｄ ｓｅｉｓｍｉｃ ｐｒｏｆｉｌｅｗｉｔｈ ｐｅａｋｆｒｅｑｕｅｎｃｙ ｏｆ６０Ｈｚ ２．２跨相带地质概念模型及地震响应特征 参照巴哈马海湾南部的水下观测和取样对现代碳酸盐前斜坡的研究成果来设计碳酸盐岩跨相带模型． 该模型以台地边缘礁为主体（以加积＋进积的成礁模型构建），从右到左分别为台坪、滩、礁、礁前斜坡相，台 地边缘礁向台地内有礁翼相滩和泻湖过渡，台地边缘礁向台地海盆方向有陡斜坡带向缓坡斜坡带过渡，以斜 坡楔的形式终结．跨不同相带的地质框架模型如图３所示． 跨不同相带地质框架模型数值模拟结果如图４所示．从图４反射波组特征可以发现，台地边缘礁相表现 出礁顶为强振幅反射、碳酸盐岩厚度明显加大、底部反射上凸、内部出现ｓ形前积现象和视频率变低等特征． 滩相的顶部为强振幅反射，碳酸盐岩厚度无明显变化，内部反射连续，滩与礁相和泻湖相间有一渐变的能量 减弱过渡现象．泻湖相的顶为强振幅反射，碳酸盐岩厚度逐渐变薄，其内部及底部有振幅减弱现象．礁前陡斜 坡相范围很窄，反射特征呈台阶式陡倾反射，且强反射突然中断．缓坡斜坡相顶面反射为较强振幅特征，内部 河海大学学报（自然科学版） 第４０卷 杂乱至断续状，上下层总的形态像喇叭状向斜坡楔收拢， 最后变成强的复合波反射． ２．３复合生物礁地质框架模型及地震响应特征 某生物礁是一个礁滩复合体，钻遇生物碳酸盐岩厚 度达５７３．５ｍ，其主体为台地边缘礁，发育多套高孔储层． 在碳酸盐岩段的中下部也有多层孔隙发育层段，局部甚 至孔隙、缝洞更为发育．经测井解释和测试，在该礁体上 部获得高产油流，下部层段为水层．从孔隙发育程度来 看，储层表现为孔隙发育、欠发育或致密层段相间出现且 厚薄不均，该特征多次出现．岩性解释结果表明，曾有多 次因海平面下降而暴露是引起上述现象的主要原因．为 此，开展精细的模型正演研究，可以帮助对地震剖面中出 现的反射特征做出合理的地质解释，预测储层发育的层 位和具体位置． 图４跨不同相带地质框架模型模拟地震剖面 Ｆｉｇ．４ Ｓｉｍｕｌａｔｅｄ ｓｅｉｓｍｉｃ ｐｒｏｆｉｌｅ ａｌｅｌｔ－（３￣ｖａｒｉｏｕｓ ｐｈａｓｅｓ Ｆｉｇ．３图３跨不同相带地质框架模型 ‘ ｏＪ０ｌ ＝ａｌ ｆｒａｍｅ ｍｏｄｅｌ籼ｓｓ ｖａｒｉｏｕｓ ｐｈａｓｅｓ 图５为该生物礁物理框架模型，其传播速度、密度值由某井 的声波和密度测井分层数据提供．模型横向长度为１００００ｍ，纵 向深度为３０００ｍ．按不同分辨率需求，模拟了２套数据（网格大 小５ ｍ×２ｍ，主频６０Ｈｚ；网格大小５ ｍ×５ ｍ，主频３Ｏ ＨＺ）．采用的 观测方式为：中间放炮，两边接收．道间距２５ｍ，炮间距２５ｍ，偏 移距一３ １７５～３ ２００ｍ，接收道数２５６道，共２００炮，６４次叠加，采 样间隔１ ＩｌｌＳ，记录时间２．５ Ｓ，第１炮的坐标位置为（１ ５００，０）． 图５某生物礁物理框架模型 模型设计时考虑了多次海平面升降的生物礁成岩环境，按 Ｆｉｇ．５ Ｐｈｙｓｉｃａｌ ｆｒａｍｅ ｍｏｄｅｌ ｆｏｒ ａ ｒｅｅｆ 退积型、加积型、进积型和加积进积型的基本模型进行组合．模 型正演结果如图６所示，图中碳酸盐岩顶界面强反射特征明显，可以清楚地区分碎屑岩与碳酸盐岩．此外，由 于横向变速造成时间成像“构造假象”（或速度陷阱），不能正确反映地下构造形态．在分辨率方面，６０Ｈｚ主频 数据明显优于３０ Ｈｚ，薄层分辨率更高． （ａ）３０Ｈｚ 图６模拟时间域成像的地震剖面 Ｆｉｇ．６ Ｓｉｍｕｌａｔｅｄ ｓｅｉｓｍｉｃ ｐｒｏｆｉｌｅ￣ａｇｍｇ ｉｎ ｔｉｍｅ ｄｏｍａｉｎ 图７为深度域成像模拟地震剖面，主频为６ｏＨｚ．从图７可以见到，礁翼附近的ｓ形前积现象比较清晰，礁的 多期次上超特点明显，礁内中间靠右的强反射同相轴表明，早期形成的礁核部位比塔礁所在位置更远离海盆方 向．对比模型正演结果与实际地震记录，两者比较接近，表明复合礁体设计合理、正演模拟方法可靠．此外，叠前 深度偏移的深度剖面外形特征与地质模型一致，说明正确的速度模型是恢复地质构造形态的关键． ３结 论 ａ．由于在生物礁、滩与上下围岩之间存在较大的波阻抗差异，其顶底分界面可以形成强地震反射以及 产生极性反转，其内部产生ｓ形前积反射、弱反射或杂乱反射，两侧出现上超反射及披覆现象．这些特征都有 助于识别碳酸盐岩生物礁． ｂ．在横向上，生物礁顶部强反射会导致其底部反射比周围反射弱．在纵向上，生物礁内不存在明显的波 第３期 但志伟，等碳酸盐岩生物礁地震响应特征的数值模拟 ３４９ 阻抗差异，其内部无法形成强反射，不利于分析生物礁内 部纵向变化特征．此外，对于大孔隙度和致密碳酸盐岩 （如礁相和台坪相），底部反射同相轴会分别出现下凹和 上凸现象，并且随着孔隙度变化幅度加大下凹和上凸现 象更加明显． Ｃ．碳酸盐岩从台坪相到滩相、台地边缘礁相及陆棚 相，其底界地震反射能量逐渐增强，而顶界反射能量变化 复杂．这种横向相带变化所引起的地震反射特征差异，有 助于识别不同相带碳酸盐岩地层． ｄ．对于碳酸盐岩生物礁储集层，高分辨率资料更能 图７模拟深度域成像的地震剖面 反映薄层的变化特征，深度域成像剖面优于时间域成像Ｆｉｇ．７ Ｓｉｍｕｌａｔｅｄ ｓｅｉｓｍｉｃ ｐｒｏｉｆｌｅ ｈｎａｇ￣ｉｎ ｄｅｐｔｈ ｄｏｍａｉｎ 剖面，可以消除生物礁下伏地层地震反射上凸或下凹现象．为此，在地震资料处理阶段要提供高分辨率的深 度域成像剖面，这也为后续的叠前或叠后波阻抗反演提供保障． 参考文献： ［１］佘德平．波场数值模拟技术［Ｊ］．勘探地球物理进展，２００４，２７（１）：１６—２１．（ＳＨＥ Ｄｅ—ｐｉｎｇ．Ｏｎ ｗａｖｅ ｆｉｅｌｄ ｎｕｍｅｒｉｃａｌ ｍｏｄｅｌｉｎｇ ｔｅｃｈｎｉｑｕｅｓｌ ＪＪ．Ｐｒｏｇｒｅｓｓ ｉｎ Ｅｘｐｌｏｒａｔｉｏｎ Ｇｅｏｐｈｙｓｉｃｓ，２１Ｘ￣４，２７（１）：１６—２１．（ｉｎ Ｃｈｉｎｅｓｅ）） ［２］张永刚．复杂介质地震波场模拟分析与应用［Ｍ］．北京：石油工业出版社，２００７． ［３］李三福，肖为，朱美娟，等．深水礁滩相储层地震模型的正演模拟及其地震响应特征分析［Ｊ］．工程地球物理学报，２０１１，８ （１）：９１—９６．（ＬＩ Ｓａｎ－ｆｕ，ＸＩＡＯ Ｗｅｉ，ＺＨＵ Ｍｅｉ—ｊｕａｎ，ｅｔ ａ１．Ｆｏｒｗａｒｄ ｍｏｄｅｌｉｎｇ ａｎｄ ａｎａｌｙｓｉｓ ｏｆ ｓｅｉｓｍｉｃ ｒｅｓｐｏｎｓｅ ｃｈａｒａｃｔｅｒｉｓｉｔｃ ｉｎ ｔｈｅ ｄｅｅｐ ｗａｔｅｒ ｒｅｅｆ－ｌｆａｔ ｒｅｓｅｒｖｏｉｒｓ ｍｏｄｅｌ［Ｊ］．Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｅｎｇｉｎｅｅｒｉｎｇ ｅＧｏｐｈｙｓｉｃｓ，２０１１，８（１）：９１—９６．（ｉｎ Ｃｈｉｎｅｓｅ）） ［４］何建军，刘家铎，鲁新便，等．基于模型正演的地震属性分析技术识别和划分碳酸盐岩储层缝洞单元［Ｊ］．石油地球物理勘 探，２００９，４４（４）：４７２－４７７．（ＨＥ Ｊｉｎａ—ｉｕｎ，ＬＩＵ Ｊｉａ—ｄｕｏ，ＬＵ Ｘｉｎ—ｂｉｎａ，ｅｔ ａ１．Ｒｅｃｏｇｎｉｔｉｏｎ ａｎｄ ｄｅｌｉｍｉｔａｔｉｏｎ ｏｆ ｆｒａｃｔｕｒｅ ａｎｄ ｃａｖｅ ｕｎｉｔｓ ｏｆ ｃａｒｂｏｎａｔｅ ｒｅｓｅｒｖｏｉｒ ｂｙ ｕｓｉｎｇ ｓｅｉｓｉｍｃ ａｔｔｉｒｂｕｔｅｓ ｎａａｌｙｔｉｃ ｔｅｃｈｎｉｕｑｅ ｂａｓｅｄ ｏｎ ｆｏｒｗａｒｄ ｍｏｄｅｌ［Ｊ］．Ｏｉｌ Ｇｅｏｐｈｙｓｉｃａｌ Ｐｒｏｓｐｅｃｔｉｎｇ，２ＯＯ９，４４（４）： ４７２－４７７．（ｉｎ Ｃｈｉｎｅｓｅ）） ［５］贺振华，黄德济．复杂油气藏地震波场特征方法理论及应用［Ｍ］．成都：四川科学技术出版社，１９９９． ［６］董良国，马在田，曹景忠，等．一阶弹性波方程交错网格高阶差分解法［Ｊ］．地球物理学报，２０００，４３（３）：８５６－８６４．（ＤＯＮＧＬｉａｎｇ —ｎｇｏ，ＭＡ Ｚａｉ—ｔｉａｎ，ＣＡＯ Ｊｉｎｇ－ｚｈｏｎｇ，ｅｔ ａ１．Ａ ｓｔａｇｇｅｒｅｄ－ｇｒｉｄ ｈｉｇｈ－ｏｒｄｅｒ ｄｉｆｅｒｅｎｃｅ ｍｅｔｈｏｄ ｏｆ ｏｎｅ—ｏｄｒｅｒ ｅｌｓａｔｉｃ ｗａｖｅ ｅｕｑａｔｉｏｎ［Ｊ］．Ｃｈｉｎｅｓｅ Ｊｏｕｍａｌ ｏｆ ｅＧｏｐｈｙｓｉｃｓ，２０００，４３（３）：８５６—８６４．（ｉｎ Ｃｈｉｎｅｓｅ）） ［７］奚先，姚姚．二维弹性随机介质中的波场特征［Ｊ］．石油地球物理勘探，２００４，３９（６）：６７９—６８５．（ＸＩ Ｘｉａｎ，ＹＡＯ Ｙａｏ．Ｗａｖｅｉｆｅｌｄ ｃｈａｒａｃｔｅｒｓ ｏｆ ２－Ｄ ｅｌｓａｔｉｃ ｒａｎｄｏｍ ｍｅｄｉｕｍ［Ｊ］．Ｏｉｌ Ｇｅｏｐｈｙｓｉｃａｌ Ｐｒｏｓｐｅｃｔｉｎｇ，２００４，３９（６）：６７９－６８５．（ｉｎ Ｃｈｉｎｅｓｅ）） ［８］白杰．射线追踪正演模拟及应用［Ｄ］．北京：中国地质大学，２００７． １ ９ ｊ ＢＵＢＥ Ｋ，ＷＡＳＨＢＯＵＲＮＥ Ｊ．Ｗａｖｅｔｒａｃｅ：ｒａｙ ｔｒａｃｉｎｇ ｆｏｒ ｔｈｅ ｐｒｏｐａｇａｔｉｏｎ ｏｆ ｂａｎｄ－ｌｉｉｍｔｅｄ ｓｉｇｎａｌｓ．Ｐａｒｔ １：Ｔｈｅｏｒｙｌ Ｃ］／／Ｅｘｐａｎｄｅｄ Ａｂｓｔｒａｃｔｓ ０ｆ ７１ Ａｎｎｕａｌ Ｉｎｔｅｍａｔ ＳＥＧ Ｍｔｇ．Ｓａｎ Ａｎｔｏｎｉｏ，ＵＳＡ：ＳＥＧ，２００ｔ：１１３１—１１３４． １１０Ｊ ＨＯＵＳＥ Ｌ，ＬＡＲＳＥＮ，ＭＡＲＴＩＮ Ｇ，ｅｔ １ａ．２Ｄ ｎａｄ ３Ｄ ｅｌｓａｔｉｃ ｍｏｄｅｈｎｇ ｗｉｈｔ ｓｈａｒｅｄ ｓｅｉｓｉｍｃ ｍｏｄｅｌｓ［Ｃ］／／Ｅｘｐａｎｄｅｄ Ａｂｓｔｒａｃｔｓ ｏｆ７２ｎａＡｎｎｕａｌ Ｉｎｔｅｍａｔ ＳＥＧ Ｍｔｇ．Ｏｂｌａｈｏｍｅ Ｃｉｔｙ：ＳＥＧ，２００２：１９７５—１９７８． ［１１］张文忠．Ｂｉｏｔ介质的交错网格差分法波场模拟研究［Ｄ］．北京：中国地质大学，２００７． ［１２］雍运动．三维黏弹性随机介质地震波场并行模拟研究［Ｄ］．吉林：吉林大学，２００７． １１３ Ｊ ＲＯＢＥＲＴＳＳＯＮ Ｊ Ｏ Ａ，ＢＬＡＮＣＨ Ｊ Ｏ，ＳＹＭＥＳ Ｗ Ｗ．Ｖｉｓｃｏｅｌａｓｔｉｃ ｉｆｎｉｔｅ—ｄｉｆｅｒｅｎｃｅ ｍｏｄｅｌｉｎｇ［Ｊ］．ｅＧｏｐｈｙｓｉｃｓ，１９９４，５９（９）：１４４４一１４５６． １１４ｊ ＣＡＲＣＩＯＮＥ Ｊ Ｍ，ＫＯＳＬＯＦＦ Ｄ，ＫＯＳＬＯＦＦ Ｒ．Ｗａｖｅ ｐｒｏｐａｇａｔｉｏｎ ｓｉｍｕｌａｔｉｏｎ ｉｎ ａ ｖｉｓｃｏｅｌｓａｔｉｃ ｍｅｉｄｍｕ［Ｊ］．ｅＧｏｐｈｙｓ Ｊ，１９８８，９５：５９７—６１１．