非均匀地震数据重建方法和其在三维叠前深度偏移处理中的应用.pdfVIP

非均匀地震数据重建方法及其在三维叠前 深度偏移处理中的应用 刘喜武 刘洪 李幼铭 1引言 研究地震数据空间插值的学者大致可以分为两类，—类研究等间隔反假频道内插，如f-x域预 测滤波插僧ⅥJ，f墩域反假频地震道插值!圳；—类是研究非均匀道插值，但不能反假频，如基于Fourier 变换信号重建的方法％ 本文基于Fourier变换带限信号重构，既能进行非均匀采样数据重建，又可以去除假频。插值或 外推时，由不完全数据(均匀或非均匀脓复出完全数据的Fourier系划51。最小二乘反演，引人DFr- 加权范数规则化项16】，用预条件共轭梯度法(PCG)求解，解稳定，收敛速度快。根据线性同相轴 在，一工域可预测理论，处理线性同相轴假频问题时，依据采样定理，采用Yule-Walker方程线性预 测方法，由带限信号的低频功率谱预测高频功率谱，去除假频。同样Yule-Walker方程亦采用预条 件共轭梯度法求解，保证解的稳定性和收敛速度。 在炮集三维叠前深度偏移处理中。要求炮集数据均匀化目尽可能覆盖速度体范围，应用本文方 法实现不规则三维炮集数据内插、外推，避免以往植入空道的简单化处理，有效地扩大了偏移孔径 信息，经实际三维地震资松}试算，偏移效果得到改善。表明本文方法的有效性和实用性。 2自适应DFT加权约束最zb-乘非均匀地震数据重建 带限地震数据重建问题栽丕为 Y=TBx=T二FH芭Fx— M (1) 其中，y为已知不规则观测值，T为取样算子，B为带限算子，X为欲重建的未知数据，F为离散 变换。带限数据重建问题可以表述为从不完全数据或不规则数据恢复未知数据的Fourier系数的反演 问题。目标函数中，选取规则化项为D删口权范数，其解受谱支持和功率谱形状约束，重建数据 的频率成分与加权函数形状相似。 ．，～(x)=ITBx-y幢+Ⅳ2X； (2) 建数据功率谱相似的加权函数。解为： 囊～-(rT+p2吉rP2F)。1Try (3) 由于X未知，加权函数要用修正的周期图迭代得到。采用预条件共轭梯度法(PCG)i匾常不直接求解 (3)式，而是求解等价的原始问题超定线性系统， 187 f T 、I ／．-、 ④ J l者阡卜。10 范数约束，故可以自适应完成非均匀地震数据重建。 3线性预测去除空间假频 对导致空间假频现象的空间采样不足、坏道、空道等情况，在插值、补道等数据重建时，很重 要的问题是去除空间假频。去假频就是要恢复出位于Nyquist波数以外信号的波数成分，需借助于 时间或频率域的有关信息，建立不同频率的波数谱之间的关系。对于t-x域有限个线性同相轴或称 平面波(在小范围应该成立)，频率谱可以用AR(自回归)模型描述，高频的AR系数可以由低 频AR系数预测得到。由于低频频率谱无假频，初始高频频率谱由低频AR系数通过线性预测得到， 再求取功率谱，则去除假频，并将其作为高频数据估计的加权函数，结合到自适应DFT加权约束 最dx--乘非均匀地震数据重建中。 去假频策略为：(1)由自适应DFF加权约束最小二乘方法估计的低频频率谱；(2)对每个低 条件共轭梯度法求解； (3)用低频得到滤波系数或AR系数预测高频频率谱。 4在炮集三维叠前深度偏移中的应用 首先对大庆某实际叠加数据进行试算，测试。图1为某实际二维叠加剖面(有缺失道)与F— K谱，图2为未进行周期图迭代的重建结果与H(谱阁3为进行周期图迭代的重建剖面及其FL K谱。对比结果表明周期图迭代重要性，D叫Ⅱ权约束的有效性，以及反假频策略的效果。说明 反假频自适应D叫Ⅱ权约束重建方法既能适于非均匀采样数据，又能反假频，目．具有自适应性。