犞犜犐介质迭代正则化偏移反演算法.pdfVIP

年 月 第 卷 第 期 ２０１４ １２ ４９ ６ 　 　 　 ·偏移成像 · 文章编号： （ ） １０００７２１０２０１４０６１０８３０８ 犞犜犐介质迭代正则化偏移反演算法 张洪宙①② 王彦飞① 　 （ 中国科学院地质与地球物理研究所，中国科学院油气资源研究重点实验室，北京 １０００２９； 中国科学院大学，北京 １０００４９） ① ② 张洪宙，王彦飞 介质迭代正则化偏移反演算法 石油地球物理勘探， ， （）： ．犞犜犐 ． ２０１４４９６ １０８３１０９０． 摘要 ＶＴＩ介质地震数据成像的关键问题之一是解决反问题的不适定性，特别是当正演算子的伴随算子无法对 　 其逆算子充分逼近时，可以采用迭代偏移反演方法求取更精确的逆算子。本文研究了ＶＴＩ介质地震偏移反演成 像的计算方法，导出了偏移成像算子方程，并提出了固定点迭代正则化算法。该方法是一种改进的迭代正则化算 法，与一般的迭代正则化方法相比，无需矩阵分解，具有更高的计算效率。本文对典型模型数据进行了数值模拟， 并对实际数据进行了试验性处理，结果表明，该偏移反演成像方法在地震成像领域具有一定的应用前景。 关键词 ＶＴＩ介质 地震偏移 正则化反演 　 　 　 中图分类号： 文献标识码： Ｐ６３１ Ａ 　　 方向传播的波并适应速度变化剧烈的介质，但由于 １ 引言 其在描述大角度传播的波时存在相位改变和振幅削 　 弱的问题，导致其难以对陡倾角界面成像。为了弥 随着地震勘探技术的不断发展，各向异性因素 补单程波偏移的倾角限制，实际地震资料处理中通 逐渐受到重视。已经证实地震各向异性普遍存在于 常采用 Ｋｉｒｃｈｈｏｆｆ偏移方法处理陡倾角地层的成像 地下介质中，特别是在沉积岩地区，它主要表现为地 问题，但 Ｋｉｒｃｈｈｏｆｆ偏移方法也存在难以适应介质 震波传播速度是传播方向的函数、体波间的相互耦 速度横向变化剧烈和不能处理单个成像点的多值旅 ［］ 合以及横波发生分裂等现象。Ｔｓｖａｎｋｉｎ１ 基于弹性 行时问题。因此这些偏移方法在对复杂构造进行处 模量张量中独立元素的个数，将具有对称特性的各 理时都得不到最佳的成像结果。随着对复杂勘探对 向异性介质分为：三斜晶系（Ｔｒｉｃｌｉｎｉｃ）、单斜晶系 象成像精度要求的不断提高，传统偏移方法的不足 （ ）、正交晶系（ ）、横向各向 日益明显。 Ｍｏｎｏｃｌｉｎｉｃ Ｏｒｔｈｏｒｈｏｍｂｉｃ 同性（ＴｒａｎｓｖｅｒｓｅｌＩｓｏｔｒｏｉｃ）。如今，地震波各向 以上分析了目前较流行的常规偏移成像方法存