地震波形反演技术的应用.pdfVIP

地质与勘探 　　　　　　　　　　　　　　天 　然 　气 　工 　业 　　　　　　　　　　　　　　　　　2009 年 2 月 地震波形反演技术的应用 1 2 张荣忠 　袁 诚 1. 中国石化胜利油田分公司物探研究院　2 . 中国石油川庆钻探工程公司 　　张荣忠等. 地震波形反演技术的应用. 天然气工业 ,2009 ,29 (2) :5254 . 　　摘 　要 　地震波形反演利用的信息不只波至时间和振幅信息 ,还包括波形细节变化信息 。在总结了波形反演 方法的分类 、求解方法 、实现波形反演的优化解法 , 以及波形反演当前的应用现状等基础上 ,指出旅行时和波形反 演相结合 ,叠前反演与叠后反演相结合是当前发展波形反演技术的两种可行思路 。 　　关键词 　地震数据 　反演 　时间域 　频率域 　应用 　　DO I :10 . 3787/ j . i ssn .. 2009 . 02 . 0 12 个问题 。因而 ,找到在波形有噪音时也稳定的目标 1 　波形反演的求解方法 函数是很有价值的。通常在存在地震噪音时 ,L 1 模 1 . 1 　梯度法非线性波形反演法 比L 2 模更为稳健 。然而 ,当剩余值变为零时 ,L 1 模 　　非线性最优化反演方法可分为随机全局最优化 存在一个奇点问题 。即使剩余值不为零 ,而是很小 方法和迭代最优化方法 。1984 年 ,Albert Tarantola 时 ,反演过程也会失败 。 [ 1 ] 提出的梯度法 非线性波形反演就是迭代最优化方 2 　实现波形反演的优化方法 法的典型代表 。地震波形反演通常用梯度类的方法 来解 ,但该算法存在局限性 : 由浅到深“逐面”反演降 2 . 1 　时间域平面波全波形反演 　　Dene s Vigh 和 E. William St ar r [3 ] 研究了时间 低了运算速度 ,浅部的反演结果影响深部的反演结 域平面波全波形反演求取速度场的方法 。共检波器 果 ;且容易陷入局部极小点 , 从而影响反演结果的 τ 道集通过 —p 变换后 , 把反射数据变换到平面波 精度 。 域 。平面波全波形反演方法就是在变换到平面波域 1. 2 　Gauss - Ne wton 法[ 2] 的反射数据上进行的。这种变换后的平面波域提供 　　Tarantola ( 1987) 定义 N ewton 法为 : 了几个选择 ,如改变平面数与在两个方向上增加面 n +1 n α( t ) - 1 t Δ n α - 1 t Δ m = m - J J J d = m - [ H ] J d a 元尺寸等 。这些特点使 3D 数据的反演速度较快 。 式中: Ha 是近似 He ssian 矩阵 。由偏导数波场的零 2 . 2 　大规模