２０１４年第２９卷第１期 褂掉虫点缸 地球物理学进展Ｐｒｏｇｒｅｓｓ ｉｎ Ｇｅｏｐｈｙｓｉｃｓ ＩＳＳＮ １００４—２９０３ ２０１４，２９（１）：０４００—０４０５ ｈｔｔｐ／／ｗｗｗ．ｐｒｏｇｅｏｐｈｙｓ．ｃｎ ＣＮ １ １—２９８２／Ｐ 王伟国，舒虎，邢涛，等．２０１４．南海神狐海域天然气水合物叠后逆时偏移处理及效果．地球物理学进展，２９（１）：０４００—０４０５，ｄｏｉ：１０． ６Ｏ３８／ｐｇ２Ｏ１４Ｏ１５７． ＷＡＮＧ Ｗｅｉ—ｇｕｏ，ＳＨＵ Ｈｕ，ＸＩＮＧ Ｔａｏ，ｅｔ ａ１．２０１４．Ｎａｔｕｒａｌ ｇａｓ ｈｙｄｒａｔｅｓ ｐｏｓｔ—ｓｔａｃｋ ｒｅｖｅｒｓｅ－ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｐｒｏｃｅｓｓｉｎｇ ａｎｄ ｅｆｆｅｃｔｓ ｉｎ Ｓｈｅｎｈｕ ａｒｅａ ｏｆ Ｓｏｕｔｈ Ｃｈｉｎａ Ｓｅａ．Ｐｒｏｇｒｅｓｓ ｉｎ Ｇｅｏｐｈｙｓｉｃｓ（ｉｎ Ｃｈｉｎｅｓｅ），２９（１）：０４００—０４０５，ｄｏｉ：１０．６Ｏ３８／ｐｇ２Ｏ１４Ｏ１５７． 南海神狐海域天然气水合物叠后逆时偏移处理及效果 Ｎａｔｕｒａｌ ｇａｓ ｈｙｄｒａｔｅｓ ｐｏｓｔ—ｓｔａｃｋ ｒｅｖｅｒｓｅ－ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｐｒｏｃｅｓｓｉｎｇ ａｎｄ ｅｆｆｅｃｔｓ ｉｎ Ｓｈｅｎｈｕ ａｒｅａ ｏｆ Ｓｏｕｔｈ Ｃｈｉｎａ Ｓｅａ 王伟国，舒 虎，邢 涛，徐华宁，李丽青，张宝金 ＷＡＮＧ Ｗｅｉ—ｇｕｏ，ＳＨＵ Ｈｕ，ＸＩＮＧ Ｔａｏ，ＸＵ Ｈｕａ－ｎｉｎｇ，ＬＩ Ｌｉ—ｑｉｎｇ，ＺＨＡＮＧ Ｂａｏ—ｊｉｎ 广州海洋地质调查局，广州５１０７６０ ＧｕａｎｇｚｈｏｕＭａｒｉｎｅ Ｇｅｏｌｏｇｉｃａｌ Ｓｕｒｖｅｙ，Ｇｕａｎｇｚｈｏｕ ５１０７６０，Ｃｈｉｎａ 摘要应用不同的偏移方法对Ｍａｒｍｏｕｓｉ速度模型进行了偏 Ａｂｓｔｒａｃｔ Ｕｓｉｎｇ ａ ｖａｒｉｅｔｙ ｏｆ ｍｅｔｈｏｄｓ ｆｏｒ Ｍａｒｍｏｕｓｉ ｍｉｇｒａｔｉｏｎ ｅｘｐｅｒｉｍｅｎｔ．Ｔｈｅ ｒｅｓｕ［ｔｓ ｓｈｏｗ ｔｈａｔ ｐｏｓｔ－ｓｔａｃｋ ｒｅｖｅｒｓｅ－ｔｉｍｅ ｍｉｇｒａｔｉｏｎ，ｔｏ ｃｏｍｐｌｅｘ ｓｔｒｕｃｔｕｒｅ，ｅｓｐｅｃｉａｌｌｙ ｔｈｅ ｉｎｔｅｒｎａｌ ｓｔｒｕｃｔｕｒｅ ｏｆ ｓａｌｔ，ｉｍａｇｅｓ ｃｌａｒｉｔｙ，ａｃｃｕｒａｔｅ，ａｎｄ ｄｅｅｐ ｉｍａｇｉｎｇ ｏｆ 移实验，实验结果表明，叠后逆时深度偏移成像清晰，对于复杂 构造，尤其是丘体内部构造刻画准确，成像效果优于所比较的 其他方法．与传统的偏移方法相比，叠后逆时偏移精度相对较 高、无倾角，能很好地适应强横向速度变化．对南海神狐海 域天然气水合物调查的一条二维测线应用了叠后逆时偏移，偏 ｓｔｒｕｃｔｕｒｅ ａｒｅ ｍｏｒｅ ｅｆｆｅｃｔｉｖｅ ｔｈａｎ ｏｔｈｅｒ ｍｅｔｈｏｄｓ ｅｍｐｌｏｙｅｄ ｈｅｒｅ． Ｃｏｍｐａｒｅｄ ｗｉｔｈ ｃｏｎｖｅｎｔｉｏｎａｌ ｍｉｇｒａｔｉｏｎ ｍｅｔｈｏｄ，ｔｈｅ ｐｏｓｔ—ｓｔａｃｋ ｒｅｖｅｒｓｅ－ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｉｓ ｍｏｒｅ ａｃｃｕｒａｔｅ，ｎｏ ｄｉｐ ｌｉｍｉｔａｔｉｏｎ， ｍｏｒｅ ａｄａｐｔｉｖｅ ｔｏ ｌａｔｅｒａｌ ｖｅｌｏｃｉｔｙ ｖａｒｉａｔｉｏｎ．Ａｐｐｌｙｉｎｇ ｔｈｅ ｐｏｓｔ— ｓｔａｃｋ ｒｅｖｅｒｓｅ－ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｔｏ ａ ｌｉｎｅ ｏｎ Ｓｈｅｎｈｕ Ａｒｅａ ｏｆ Ｓｏｕｔｈ Ｃｈｉｎａ Ｓｅａ，ｔｈｅ ｉｍａｇｅ ｉｓ ｂｅｔｔｅｒ ｔｈａｎ ｔｈａｔ ｏｂｔａｉｎｅｄ ｆｒｏｍ ｏｎｅ－ｗａｙ ｗａｖｅ ｅｑｕａｔｉｏｎ ｐｒｅ－ｓｔａｃｋ ｄｅｐｔｈ ｍｉｇｒａｔｉｏｎ．Ｔｈｕｓ，ｐｏｓｔ—ｓｔａｃｋ ｒｅｖｅｒｓｅ—ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｓｈｏｕｌｄ ｂｅ ｇｉｖｅｎ ｍｏｒｅ ｅｙｅｓ ｉｎ ｇａｓ ｈｙｄｒａｔｅ 移效果甚至优于单程波波动方程的叠前深度偏移．因此，天然 气水合物地震资料处理中，叠后逆时偏移应是一种值得重视的 偏移方法． 关键词叠后逆时偏移；Ｍａｒｍｏｕｓｉ模型；南海；神狐海域；天然 气水合物 中图分类号文献标识码Ｐ６３１，Ｐ７３８ Ａ ｄａｔａ ｐｒｏｃｅｓｓｉｎｇ． Ｋｅｙｗｏｒｄｓ ｐｏｓｔ—ｓｔａｃｋ ｒｅｖｅｒｓｅ－ｔｉｍｅ ｍｉｇｒａｔｉｏｎ；Ｍａｒｍｏｕｓｉ ｄｏｉ：１０．６Ｏ３８／ｐｇ２Ｏ１４Ｏ１５７ ｍｏｄｅｌ；Ｓｏｕｔｈ Ｃｈｉｎａ Ｓｅａ；Ｓｈｅｎｈｕ ａｒｅａ；ｇａｓ ｈｙｄｒａｔｅｓ ０引言 震资料，详细对比了经典Ｋｉｒｃｈｈｏｆｆ偏移方法、单程波动方程 偏移方法与逆时偏移方法．美国ＧＸＴ公司采用调谐算法和 提高计算机计算能力等手段打破了逆时偏移计算效率的瓶 颈，先后在墨西哥湾、西非、北海等地完成了近３Ｏ个逆时偏 逆时偏移理论成熟于２０世纪８Ｏ年代，随后国内外很多 学者对逆时偏移理论进行了研究（Ｗｕ，１９９６；尧德中等， １９９３；何兵寿等，２００８；Ｚｈａｎｇ ｅｔ ａ１．，２００７；Ｆｌｅｔｃｈｅｒ ｅｔ ａ１．， ２００９；张美根等，２００１；李博等，２０１０；杜启振等，２００９）．１９８２ 年ＳＥＧ年会上，Ｗｈｉｔｅｍｏｒｅ（Ｗｈｉｔｅｍｏｒｅ，１９８２）（Ｗｈｉｔｅｍｏｒｅ， 移项目（杨琴勇等，２０１０）．值得欣喜的是，目前国内在逆时偏 移的应用方面取得了很大的突破，中国科学院自主开发的所 谓ＣＰＰＣ技术（刘洪等，２０１０），即ＧＰＵ／ＣＰＵ协同并行计算 技术已经应用于国内一部分的油田企业． Ｎ．Ｄ，１９８２）首次提出了逆时偏移的思想．此后，Ｌｏｅｗｅｎｔｈａｌ （Ｌｏｅｗｅｎｔｈａｌ ｅｔ ａ１．，１９８３）等探讨了空间一频率域的逆时偏移 原理，并对纵波资料应用了叠后逆时偏移，取得了较好的应 目前，叠前处理成为地震资料处理的主流，叠前时间偏 移更是已经常态化，对于复杂构造成像，叠前深度偏移也有 了较为广泛的应用，叠后偏移处理已经越来越不为处理员所 应用．然而，对于我国南海海域天然气水合物勘探来说，由于 天然气水合物独有的一些特点，如振幅空白带、极性反转及 与地层斜交（Ｋｖｅｎｖｏｌｄｅｎ，１９９３；宋海滨等，２００１；徐华宁等， 用效果．此时的逆时偏移成像技术主要还在理论研究阶段， 鲜有实际地震资料的应用． 进入２１世纪，计算机硬件有了飞速的发展，同时烃类资 源勘探难度的不断增加，使得逆时偏移技术应用于实际成为 可能．Ｆａｒｍｅｒ（Ｆａｒｍｅｒ ｅｔ ａ１．，２００９）等应用墨西哥湾实际地 收稿日期 ２０１３—０７—１０；修回日期２０１３—１１—１８． 投稿网址２０１０），从而针对天然气水合物的处理就需以保幅、保真为前 ｈｔｔｐ／／ｗｗｗ．ｐｒｏｇｅｏｐｈｙｓ．ｃｎ 基金项目 国家自然科学基金面上项目（４１１７６０５６）资助． 作者简介 王伟国，男，１９８２年生，江西萍乡人，硕士研究生，工程师，主要从事海洋地震资料成像研究．（Ｅ－ｍａｉｌ：ｈｕａｙｉｎｇ１００９＠１６３．ｃｏｍ） ４０４ 地球物理学进展Ⅵ ｐｒｏｇｅｏｐｈｙｓ．ｃｎ ＡｔＶ（ｚ， ）Ｍ（ｚ， ，ｎａｔ）． （６） 假设Ｔ为时间剖面最大的记录时间，ｕ（ｘ，ｚ，Ｔ＋Ａｔ）和ｕ（ｘ， ，Ｔ＋２Ａｔ）为零，由方程（５）（６）就可以完成从最大记录时间 Ｔ到ｔ一０的逆时延拓计算，由爆炸反射界面成像原理可知， 当ｔ＝Ｏ时刻时，深度剖面ｕ（ｘ， ，Ｏ）即为逆时深度偏移的最 终偏移剖面（贺振华，１９８９）． １．２稳定性条件 关于叠后逆时深度偏移的稳定性条件，张文生等于 １９９８年有过详细的探讨，此处不再赘述，具体请参考文献 （张文生等，１９９８）． ２模型数据实验 以Ｍａｒｍｏｕｓｉ速度模型为例，如图１所示，模型横向 ９２００ ｍ，纵向３０００ ｍ，横向采样间隔１２．５ ｍ，纵向采样间隔 ４ ｒｒＬ利用爆炸反射理论进行叠后模拟，炸药震源激发，震源 子波为１５ Ｈｚ的Ｒｉｃｋｅｒ子波，时间采样４ ｍｓ，利用四阶有限 差分模拟． 在本例中，Ｆ－Ｘ域的单程波波动方程叠前深度偏移（以 下简称”Ｆ－Ｘ偏移”）是基于ＣＤＰ道集实验的，观测系统为左 端接收，每炮９６道接收，最小偏移距２００ ｍ，炮间隔２５ ｍ，道 间隔２５ ｍ，首炮位于ｚ一２５７５ ｍ处，末炮位与ｘ＝８５５０ ｍ，总 共２４０炮，震源和采样率和叠后模拟相同．其对应的速度模 型段为图（ａ）红色方框所示，偏移结果如图２（ｃ）所示．对比三 种偏移方法所得到的剖面，如图２所示，单从成像质量上看， Ｋｉｒｃｈｈｏｆｆ叠后深度偏移显然对岩丘边界和岩丘下部成像很 不理想，畸变严重；对比Ｆ－Ｘ偏移和叠后逆时偏移剖面，二者 均能对复杂构造进行有效的偏移，但笔者认为叠后逆时偏移 对岩丘体内部以及２５００ ｍ处油性圈闭的成像更为清晰， １５００ ｍ以上的层位细节表现更为突出，而且从偏移的耗时 来说，叠后逆时偏移更为经济一些． ３实际资料应用 以南海神狐海域天然气水合物调查的一条二维测线为 例．基本采集参数如表１所示： 表１测线采集参数 Ｔａｂｌｅ １ Ａｃｑｕｉｓｉｔｉｏｎ ｐａｒａｍｅｔｅｒ 测线水深超过３１００ ｍ，通过一系列的叠前噪音消除、海 底多次波衰减、反褶积等处理，得到较高质量的共中心点道 集（以下称“ＣＤＰ道集”），对ＣＤＰ道集进行叠加，得到叠后逆 时偏移所需的叠加剖面，如图３（ｂ）所示；另外，基于ＣＤＰ道 集分别做了叠后Ｋｉｒｃｈｈｏｆｆ深度偏移和Ｆ－Ｘ域叠前深度偏 移，三种方法所用的速度一深度模型如图３（ａ）所示，目的层位 为海底以及凹陷的底界．速度模型是横向变化的，海底起伏 较大，沉积层位集中，沉积底界坡度陡，崎岖不平． 首先从整体效果来看，图４（ｂ）和（ｃ）层位突出，波组特征 明显，与图４（ａ）相比，好像图４（ａ）有着更好的分辨率，同相 轴更细，其实不然．通过放大对比可发现，图４（ａ）中绕射波有 很大一部分没有收敛，同相轴上有很多往上翘的”尾巴（残余 绕射波）”，目的层不够突出，整体波组特征不强，这就导致解 释人员无法很好地追踪层位．再看图４（ｂ）和（ｃ）的放大细节， 层位归位比较准确，”画弧”少，尤其是图４（ｃ），目的层较图４ （ｂ）更为连续，分辨率更高．当然，对于叠后逆时偏移剖面而 言，如图４（ｃ），剖面上还是存在一些偏移噪音，但是对于天然 气水合物地震勘探而言，其主要目的还是寻找似海底反射 （Ｂｏｔｔｏｍ Ｓｉｍｕｌａｔｉｎｇ Ｒｅｆｌｅｃｔｏｒ，简称ＢＳＲ）（姚伯初，１９９８）以 及游离气等特征，天然气水合物出现的位置通常在海底以下 ４００ ｍ范围之内，属于浅层，而叠前深度偏移算法对于中深 部复杂构造成像有着较好的效果，从图４（ｂ）的放大效果图 可以看出，对于ＣＤＰ号从７４００到８４００，深度从４６００ ｍ到 ５０００ ｍ的这一段剖面，同相轴比较模糊，分辨率低，当然造 成这样的结果不一定是算法本身的原因，可能跟速度模型准 确度、叠前道集处理等因素有关．但同时也可以看到，在图４ （ｃ）的放大效果剖面，ＣＤＰ道集、速度模型和叠前深度偏移所 用到的是一样的，而叠后逆时偏移（深度域）剖面中的同相轴 就更为清晰，视分辨率也较图４（ｂ）要高． ４结论与认识 实际地震资料处理过程中，不同的偏移方法所应用的范 围不同．模型数据的对比以及实际资料的应用效果表明，叠 后逆时偏移在起伏海底地形以及构造比较复杂地区成像效 果理想，表现在： （１）无倾角．在实际二维地震资料中，地层断面以及 沉积层边界角度较陡，甚至接近垂直方向，如图３所示，但叠 后逆时偏移仍然可以很好地对构造进行成像． （２）各种波场能量得到很好地聚焦．对目的层位的成像 可以看到，残余绕射波能量很少，层位连续，界面清晰，波组 特征明显，分辨率提高． （３）浅部构造刻画突出．对于天然气水合物地震资料处 理而言，浅中层的成像更为关键，和叠前深度偏移方法相比， 叠后逆时偏移的浅层成像能力更为突出． 当然，叠后逆时偏移方法由于是基于爆炸反射界面理论 的假设，基于共中心点叠加得到的叠加剖面也存在着诸如水 平地层、速度横向变化小等等假设条件，所以叠后逆时偏移 也存在着”先天性”的不足，但是从偏移成像对比的结果来 看，叠后深度偏移与叠前深度偏移还是存在可比性，浅部成 像甚至优于叠前深度偏移．对于天然气水合物地震资料处理 来说，叠后深度偏移方法应该是一种较好的选择，有一定的 推广价值，值得重视． 致谢对审稿专家提出的众多有益的修改意见表示感谢 Ｄｕ Ｑ Ｚ，Ｑｉｎ ２００９．Ｍｕｈｉｃｏｍｐｏｎｅｎｔ ｐｒｅｓｔａｃｋ ｒｅｖｅｒｓｅ￣ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｏｆ ｅｌａｓｔｉｃ ｗａｖｅｓ ｉｎ ｔｒａｎｓｖｅｒｓｅ ｉｓｏｔｒｏｐｉｃ ｍｅｄｉｕｍ［Ｊ］． Ｃｈｉｎｅｓｅ Ｊ．Ｇｅｏｐｈｙｓ．（ｉｎ Ｃｈｉｎｅｓｅ），５２（３）：８０１—８０７． Ｆｌｅｔｃｈｅｒ Ｒ Ｐ，Ｄｕ Ｘ，Ｆｏｗｌｅｒ Ｐ Ｊ．２００９．Ｒｅｖｅｒｓｅ ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｉｎ ２０１４，２９（１） 王伟国，等：南海神狐海域天然气水合物叠后逆时偏移处理及效果（Ⅵｎ ｐｒｏｇｅｏｐｈｙｓ．ｃｎ）４０５ ｔｉｌｔｅｄ ｔｒａｎｓｖｅｒｓｅｌｙ ｉｓｏｔｒｏｐｉｃ（ＴＴＩ）ｍｅｄｉａ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，７４ （６）：Ｗ＿ＣＡ１７９一ＷＣＡ１８７． Ｆａｒｍｅｒ Ｐ，Ｚｈｏｕ Ｚ Ｚ，Ｊｏｎｅｓ ｎ ２００９．Ｔｈｅ ｒｏｌｅ ｏｆ ｒｅｖｅｒｓｅ ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｉｎ ｉｍａｇｉｎｇ ａｎｄ ｍｏｄｅｌ ｅｓｔｉｍａｔｉｏｎ［Ｊ］．Ｔｈｅ Ｌｅａｄｉｎｇ Ｅｄｇｅ，２８（４）：４３６—４４１． Ｆｏｒｎｂｅｒｇ １３．１９８７．Ｔｈｅ ｐｓｅｕｄｏｓｐｅｃｔｒａｌ ｍｅｔｈｏｄ；Ｃｏｍｐａｒｉｓｏｎｓ ｗｉｔｈ ｆｉｎｉｔｅ ｄｉｆｆｅｒｅｎｃｅｓ ｆｏｒ ｔｈｅ ｅｌａｓｔｉｃ ｗａｖｅ ｅｑｕａｔｉｏｎ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ， ５２（４）：４８３－５０１． Ｇａｚｄａｇ Ｊ．１９８１．Ｍｏｄｅｌｉｎｇ ｏｆ ｔｈｅ ａｃｏｕｓｔｉｃ ｗａｖｅ ｅｑｕａｔｉｏｎ ｗｉｔｈ ｔｒａｎｓｆｏｒｍ ｍｅｔｈｏｄｓ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，４６（６）：８５４－８５９． Ｈｅ Ｂ Ｓ，Ｚｈａｎｇ Ｈ Ｘ，Ｚｈａｎｇ Ｊ．２００８．Ｐｒｅｓｔａｃｋ ｒｅｖｅｒｓｅ－ｔｉｍｅ ｄｅｐｔｈ ｍｉｇｒａｔｉｏｎ ｏｆ ａｒｂｉｔｒａｒｉｌｙ ｗｉｄｅ－ａｎｇｌｅ ｗａｖｅ ｅｑｕａｔｉｏｎｓ［Ｊ］．Ａｃｔａ Ｓｅｉｓｍｏｌｏｇｉｃａ Ｓｉｎｉｃａ（ｉｎ Ｃｈｉｎｅｓｅ），３０（５）：４９１—４９９． Ｈｅ Ｚ Ｈ．１９８９．Ｒｅｆｌｅｃｔｉｏｎ ｓｅｉｓｍｉｃ ｄａｔａ ｍｉｇｒａｔｉｏｎ ｐｒｏｃｅｓｓｉｎｇ ａｎｄ ｉｎｖｅｒｓｉｏｎ ｍｅｔｈｏｄｓ［Ｍ］．Ｃｈｏｎｇｑｉｎｇ：Ｃｈｏｎｇｑｉｎｇ Ｕｎｉｖｅｒｓｉｔｙ ｐｒｅｓｓ． Ｋｖｅｎｖｏｌｄｅｎ Ｋ八１９９３．Ｇａｓ ｈｙｄｒａｔｅｓ—ｇｅｏｌｏｇｉｃａｌ ｐｅｒｓｐｅｃｔｉｖｅ ａｎｄ ｇｌｏｂａｌ ｃｈａｇｎｅ［Ｊ］．Ｒｅｖｉｅｗｓ ｏｆ Ｇｅｏｐｈｙｓｉｃｓ，３１（２）：１７３—１８７． Ｋｒｅｉｓｓ Ｈ ０，Ｏｌｉｇｅｒ Ｊ．１９７２．Ｃｏｍｐａｒｉｓｏｎ ｏｆ ａｃｃｕｒａｔｅ ｍｅｔｈｏｄｓ ｆｏｒ ｔｈｅ ｉｎｔｅｇｒａｔｉｏｎ ｏｆ ｈｙｐｅｒｂｏｌｉｃ ｅｑｕａｔｉｏｎｓ［Ｊ］．Ｔｅｌｌｕｓ，２４（３）： １９９—２１５． Ｌｏｅｗｅｎｔｈａ１ Ｄ，Ｍｕｌｆｔｉ Ｉ Ｒ １９８３．Ｒｅｖｅｒｓｅｄ ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｉｎ ｓｐａｔｉａｌ ｆｒｅｑｕｅｎｃｙ ｄｏｍａｉｎ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，４８（５）：６２７—６３５． Ｌｏｅｗｅｎｔｈａｌ Ｄ，Ｌｕ Ｌ，Ｒｏｂｅｒｓｏｎ Ｒ，ｅｔ ａ１．１９７６．Ｔｈｅ ｗａｖｅ ｅｑｕａｔｉｏｎ ａｐｐｌｉｅｄ ｔｏ ｍｉｇｒａｔｉｏｎ［Ｊ］．Ｇｅｏｐｈｙｓｉｃａｌ Ｐｒｏｓｐｅｃｔｉｎｇ，２４（２）： ３８０—３９９． Ｕ Ｂ，Ｌｉｕ Ｈ Ｗ，Ｌｉｕ Ｇ Ｆ，ｅｔ ａ１．２０１０．Ｃｏｍｐｕｔａｔｉｏｎａ１ ｓｔｒａｔｅｇｙ ｏｆ ｓｅｉｓｉｍｃ ｐｒｅ－ｓｔａｃｋ ｒｅｖｅｒｓｅ ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｏｎ ＣＰＵ／ＧＰＵ［Ｊ］． Ｃｈｉｎｅｓｅ Ｊ．Ｇｅｏｐｈｙｓ．（ｉｎ Ｃｈｉｎｅｓｅ），５３（１２）：２９３８—２９４３． Ｌｉｕ Ｈ，Ｔｏｎｇ Ｘ Ｌ，Ｌｉｕ Ｑ，ｅｔ ａ１．２０１０．Ｃｏｌｌａｂｏｒａｔｉｖｅ ｐａｒａｌｌｅ１ ｃｏｍｐｕｔｉｇｎ ｔｅｃｈｎｏｌｏｇｙ ｏｎ ｏｉｌ ａｎｄ ｇａｓ ｓｅｉｓｍｉｃ ｄａｔａ ｐｒｏｃｅｓｓｉｎｇ ｗｉｔｈ ＧＰＵ／ＣＰＵ［Ｍ］．Ｂｅｂｉｎｇ：Ｐｅｔｒｏｌｅｕｍ Ｉｎｄｕｓｔｒｙ Ｐｒｅｓｓ． Ｓ０ｎｇ Ｈ Ｂ，Ｓｏｎｇ Ｌ Ｘ，Ｙａｎｇ Ｓ Ｘ，ｅｔ ａ１．２００１．Ｇｅｏｐｈｙｓｉｃａｌ ｒｅｓｅａｒｃｈｅｓ ｏｎ ｍａｒｉｎｅ ｇａｓｈｙｄｒａｔｅｓ（ＩＩ）：ｓｅｉｓｍｉｃ ｍｅｔｈｏｄｓ［Ｊ］． Ｐｒｏｇｒｅｓｓ ｉｎ Ｇｅｏｐｈｙｓｉｃｓ（ｉｎ Ｃｈｉｎｅｓｅ），１６（３）：１１０—１１８． Ｗｈｉｔｅｍｏｒｅ Ｎ ｎ １９８２．Ｍｉｇｒａｔｉｏｎ：Ｆｕｎｄａｍｅｎｔａｌｉｓｓｕｅｓ ａｎｄ ｆｕｔｕｒｅ ｄｅｖｅｌｏｐｍｅｎｔｓ［Ａ］／／５２州ＡｎｒＬ Ｉｎｔｅｒｎｅｔ．Ｍｅｅｔｉｎｇ［Ｃ］．Ｄａｌｌａｓ： Ｓｏｃ．Ｅｘｐｌｏｒ．Ｇｅｏｐｈｙｓ．Ｐｕｂ１．Ｃｏ．，５２０－５２０． Ｗｕ Ｗ Ｊ．１９９６．Ａｎａｌｙｓｉｓ ｏｆ ｈｉｇｈｅｒ－ｏｒｄｅｒ，ｆｉｎｉｔｅ－ｄｉｆｆｅｒｅｎｃｅ ｓｃｈｅｍｅｓ ｉｎ ３一Ｄ ｒｅｖｅｒｓｅ－ｔｉｍｅ ｍｉｇｒａｔｉｏｎ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，６１（３）：８４５—８５６． Ｗａｎｇ Ｗ Ｇ，Ｘｉｏｎｇ Ｓ Ｊ，Ｘｕ Ｈ Ｎ。ｅｔ ａ１．２０１２．Ｒｅｖｅｒｓｅ－ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｕｓｉｇｎ ｐｓｅｕｄｏ－ｓｐｅｃｔｒａ１ ｍｅｔｈｏｄ ｆｏｒ ｔｉｌｔｅｄ ＴＴＩ ｍｅｄｉａ ［Ｊ］．ＯＧＰ（ｉｎ Ｃｈｉｎｅｓｅ），４７（４）：５６６—５７２． Ｘｕ Ｈ Ｎ，Ｓｈｕ Ｈ，Ｌｉ Ｌ Ｑ，ｅｔ ａ１．２００９．３－Ｄ ｓｅｉｓｍｉｃ ｄａｔａ ｐｒｏｃｅｓｓｉｎｇ ｔｅｃｈｎｉｑｕｅｓ ｆｏｒ ｇａｓ ｈｙｄｒａｔｅ ｂｙ ｔｈｅ ｓｉｎｇｌｅ－ｓｏｕｒｃｅ ａｎｄ ｓｉｎｇｌｅ－ｃａｂｌｅ ａｃｑｕｉｓｉｔｉｏｎ ｍｅｔｈｏｄ［Ｊ］．Ｐｒｏｇｒｅｓｓ ｉｎ Ｇｅｏｐｈｙｓｉｃｓ（ｉｎ Ｃｈｉｎｅｓｅ），２４ （５）：１８０１—１８０６． Ｘｕ Ｈ Ｎ，ｒａｎｇ Ｓ Ｘ，Ｚｈｅｎｇ Ｘ Ｄ，ｅｔ ａ１．２０１０．Ｓｅｉｓｍｉｃ ｉｄｅｎｔｉｆｉｃａｔｉｏｎ ｏｆ ｇａｓ ｈｙｄｒａｔｅ ａｎｄ ｉｔｓ ｄｉｓｔｒｉｂｕｔｉｏｎ ｉｎ Ｓｈｅｎｈｕ Ａｒｅａ．Ｓｏｕｔｈ Ｃｈｉｎａ Ｓｅａ［Ｊ］．Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｇｅｏｐｈｙｓｉｃｓ（ｉｎ Ｃｈｉｎｅｓｅ），５３（７）： １６９１—１６９８． Ｙａｇｎ Ｑ Ｙ，Ｄｕａｎ Ｘ Ｂ［Ｊ］．２０１０．Ｄｅｖｅｌｏｐｍｅｎｔ ｓｔａｔｕｓ ａｎｄ ｔｒｅｎｄ ｏｆ ｒｅｖｅｒｓｅ ｔｉｍｅ ｍｉｇｒａｔｉｏｎ［Ｊ］．ＧＰＰ（ｉｎ Ｃｈｉｎｅｓｅ），４９（１）：９２—９８． Ｙａｏ Ｄ Ｚ，Ｚｈｏｕ Ｘ Ｘ．１９９３．Ｅｌａｓｔｉｃ ｗａｖｅ ｒｅｖｅｒｓｅ ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｏｆ ｌｏｎｇｉｔｕｄｉｎａｌ ａｎｄ ｔｒａｎｓｖｅｒｓｅ ｗａｖｅ ｓｅｐａｒａｔｅ ｉｍａｇｅｒｙ［Ｊ］．Ｃｈｉｎｅｓｅ Ｊ．Ｇｅｏｐｈｙｓ．（ｉｎ Ｃｈｉｎｅｓｅ），３６（５）：６６５—６７３． Ｙａｏ Ｂ Ｃ．１９９８．Ｐｒｅｌｉｍｉｎａｒｙ ｅｘｐｌｏｒａｔｉｏｎ ｏｆ ｇａｓ ｈｙｄｒａｔｅ ｉｎ ｔｈｅ ｎｏｒｔｈｅｒｎ ｍａｒｇｉｎ ｏｆ ｔｈｅ Ｓｏｕｔｈ Ｃｈｉｎａ Ｓｅａ［Ｊ］．Ｍａｒｉｎｅ Ｇｅｏｌｏｇｙ＆ Ｑｕａｔｅｒｎａｒｙ Ｇｅｏｌｏｇｙ（ｉｎ Ｃｈｉｎｅｓｅ），１８（４）：１１－１８． Ｚｈａｎｇ Ｙ，Ｓｕｎ Ｊ，Ｇｒａｙ Ｓ ２００７．Ｒｅｖｅｒｓｅ－ｔｉｍｅ ｍｉｇｒａｔｉｏｎ．．ａｍｐｌｉｔｕｄｅ ａｎｄ ｉｍｐｌｅｍｅｎｔａｔｉｏｎ ｉｓｓｕｅｓ ｒＣ］／／Ｅｘｐａｎｄｅｄ Ａｂｓｔｒａｃｔｓ ｏｆ ７７ｔｈ ＳＥＧ Ｍｔｇ．，２１４５—２１４８． Ｚｈａｎｇ Ｍ Ｇ，Ｗａｎｇ Ｍ ２００１．Ｐｒｅｓｔａｃｋ ｆｉｎｉｔｅ ｅｌｅｍｅｎｔ ｒｅｖｅｒｓｅ－ｔｉｍｅ ｍｉｇｒａｔｉｏｎ ｆｏｒ ａｎｉｓｏｔｒｏｐｉｃ ｅｌａｓｔｉｃ ｗａｖｅｓ［Ｊ］．Ｃｈｉｎｅｓｅ Ｊ．Ｇｅｏｐｈｙｓ． （ｉｎ Ｃｈｉｎｅｓｅ），４４（５）：７１１—７１９． Ｚｈａｎｇ Ｗ Ｓ，Ｚｈａｎｇ Ｇ Ｑ，Ｈａｏ Ｘ Ｊ．１９９８．Ｐｏｓｔｓｔａｃｋ ｒｅｖｅｒｓｅ－ｔｉｍｅ ｄｅｐｔｈ ｉｍｇｒａｔｉｏｎ ｆｏｒ ｃｏｍｐｌｅｘ ｓｔｒｕｃｔｕｒｅｓ［Ｊ］．ＧＰＰ（ｉｎ Ｃｈｉｎｅｓｅ）， ３７（４）：５５—６３． 附中文参考文献 杜启振，秦童．２００９．横向各向同性介质弹性波多分量叠前逆时偏 移［Ｊ］．地球物理学报，５２（３）：８０１—８０７． 何兵寿，张会星，张晶．２００８．任意广角波动方程叠前逆时深度偏 移［Ｊ］．地震学报，３０（５）：４９１—４９９． 贺振华．１９８９．反射地震资料偏移处理与反演方法［Ｍ］．重庆：重庆 大学出版社． 李博，刘红伟，刘国锋，等．２０１０．地震叠前逆时偏移算法的ＣＰＵ／ ＧＰＵ实施对策［Ｊ］．地球物理学报，５３（１２）：２９３８—２９４３． 刘洪，佟小龙，刘钦，等．２０１０．油气地震勘探数据处理ＧＰＵ／ＣＰＵ 协同并行计算技术［Ｍ］．北京：石油工业出版社． 宋海斌，松林修，杨胜雄，等．２００１．海洋天然气水合物的地球物理 研究（ＩＩ）：地震方法［Ｊ］．地球物理学进展，１６（３）：１１０－１１８． 王伟国，熊水金，徐华宁，等．２０１２．ＴＴＩ介质各向异性伪谱法逆时 偏移［Ｊ］．石油地球物理勘探，４７（４）：５６６—５７２． 徐华宁，舒虎，李丽青，等．２００９．单源单缆方式采集的天然气水合 物三维地震数据处理技术［Ｊ］．地球物理学进展，２４（５）： １８０卜１８０６． 徐华宁，杨胜雄，郑晓东，等．２０１０．南中国海神狐海域天然气水合 物地震识别及分布特征［Ｊ］．地球物理学报，５３（７）： １６９卜１６９８． 杨勤勇，段心标．２０１０．逆时偏移技术发展现状与趋势［Ｊ］．石油物 探，４９（１）：９２—９８． 尧德中，周熙襄．１９９３．纵横波成像的弹性波逆时偏移方法 ［Ｊ］．地球物理学报，３６（５）：６６５－６７３． 姚伯初．１９９８．南海北部陆缘天然气水合物初探口］．海洋地质与第 四纪地质，１８（４）：１１－１８． 张美根，王妙月．２００１．各向异性弹性波有限元叠前逆时偏移［Ｊ］． 地球物理学报，４４（５）：７１１－７１９． 张文生，张关泉，郝现军．１９９８．复杂构造叠后逆时深度偏移［Ｊ］． 石油物探，３７（４）：５５—６３．
