三维地震数据离散光滑插值的共轭梯度法-北京玛格泰克科技发展有限.PDFVIP

第 45 卷 第 5 期 地 球 物 理 学 报 Vol . 45 , No . 5 　　　　　　 　　　　　　 　2002 年 9 月 CHINESE JOURNAL OF GEOPHYSICS 　 Sep . , 2 0 0 2 ( ) [文章编号] 　0001 - 5733 2002 - 05 - 0691 - 09 [ 中图分类号] 　P 631 三维地震数据离散光滑插值的共轭梯度法 李 　冰 　刘 　洪 　李幼铭 ( 中国科学院地质与地球物理研究所 ,北京 　100101) [摘 　要] 　针对三维地震数据插值 ,提出采用 Laplacian 算子进行光滑约束的插值方法 ,并借 鉴 Mallet 研究的离散光滑插值思路 ,采用预条件共轭梯度法 ,直接生成网格节点上的值 ,从而 回避寻求满足插值方程的函数. 为了实现其中Laplacian 算子的快速求逆 ,文中引入 Claerbout 螺旋坐标系谱因式分解理论. 在螺旋坐标系下 ,Laplacian 算子的表示矩阵具有 Toeplitz 结构 , 其快速求逆可由谱法LU 分解实现. 基于二维离散光滑插值 ,文中还给出共轭梯度法与 NMO 相结合的沿时间切片逐层处理的离散光滑插值流程. 最后 ,应用该方法对模型数据和实际三 维地震数据进行了处理. [ 关键词] 　Laplacian 算子 ,离散光滑道插值 ,预条件共轭梯度 ,螺旋坐标系谱因式分解 ,NMO. 1 　引　　言 众所周知 ,无论是陆上还是海上 ,地震勘探数据采集通常都会出现空间采样密度不 足 、坏道和数据分布不规则等问题[ 1 —5 ] . 这些问题常常影响到地震处理效果 ,并最终妨碍 正确的地质解释 。解决办法有两种 :加大采样密度和在处理时进行空间内插. 加大采样 密度会增加成本 ,降低效益. 而由于不规则地形 、禁区或仪器故障等因素也会引发类似问 题 ,此时即使加大采样密度 ,也无济于事. 近年来 ,利用三维地震资料进行的三维波动方程叠前深度偏移在油气勘探中的作用 已越发显著. 但是三维地震勘探野外施工难 、成本高 ,许多油田退而采取大线距和道距的 经济三维数据采集的方法. 而用波动方程叠前深度偏移法进行实际资料处理时 ,首先就 必须将经济三维地震数据插值成处理流程所需的线距和道距尺度. 虽然 自1981 年Larnar [6 ] 等 首次实现二维道内插之后 , 国内外的地球物理学家都十分重视道插值的研究工作 ,发 展了诸如多项式拟合法[7 ] 、最佳内插算子法[8 ,9 ] 、f x 域内插法[ 10 —13 ] 、人工智能内插法[ 14 ] 、 [ 15 ] [ 16 ] [ 17 —20 ] τ p 域内插法 、非线性相关法 等多种道插值技术 . 但是 , 当上述插值方法由二 维插值推广应用于三维插值时 ,却都存在诸如运算时间长 、计算精度不够 、参数选择过多、 结果不稳定等问题. 因此 ,寻求快速而稳定的插值方法一直是三维地震资料预处理的关 [ 收稿 日期] 　2001 - 02 - 05 收到 ,2002 - 01 - 28 收到修定稿. ( ) ( ) [基金项 目] 　中国科学院知识创新工程重大项 目资助课题 KZCXIY01