基于逆散射级数的子波提取方法.pdfVIP

·610· 中国地球物理2011 (18)储层地球物理 基于逆散射级数的子波提取方法 estimationbasedoninverse series Wavelet method scattering 刘华锋+陈小宏侯波 Liu Chen HouBo HuafengXiaohong 中国石油大学(北京)CNPc物探重点实验室北京昌平102249 在多次波衰减过程中，基于逆散射级数反演的多次波衰减方法，虽然不需要任何地下结构和速度信息， 但它需要震源子波作为已知信息。Ikelle，eta1．(1997)基于逆散射理论，以数据中含有自由表面多次波为前 提，根据一阶自由表面多次波的衰减公式，利用最小二乘准则，提出了不依赖于地下结构信息的子波估计方 法。本文针对基于逆散射级数的2D子波提取方法计算量大、耗时长的问题，对2D计算公式进行简化降维， 并考虑了震源和检波器的深度，获得了计算速度快、满足精度要求的1．5D逆散射级数子波提取方法。 1．基本原理 如果假设模型仅是深度和偏移距的函数，则(五，Xg)点处的格林函数G0(五，赡，，)仅与炮检距 Xoffs。=工。一Xs和时间t有关。从鬼波衰减后的数据Do出发，预测自由表面多次波的各阶级数可表示为， 协(七，缈)2瓦丽1 g矿‘勺”J仇一·(七，∞)侥(七，叫刀；l，2，3，… (1) 其中，k是对应于确触的波数，(13是对应于t的频率，-5-(动是震源特征，肛是参考介质的密度，￡和《 分别是震源和检波器位于水下的深度，由参考介质速度CO得出的纵向波数q的表达式为， q=sign(o．,)厄玩河 (2) 假设一阶多次波消除后的数据比消除前的数据能量低，且子波仅是时间的函数，该问题简化为一个线 性问题。同时假定通过最佳衰减一阶多次波得出的子波，也能满足高阶多次波的衰减要求。因为DI(七，缈)不 仅含有一阶多次波的能量，还包含高阶多次波的能量，而且高阶级数项中依然含有一阶自由表面多次波的 能量，所以在仅截取前两项时，会造成一定的误差。假定该误差很小，并将其定义为￡，，则有， D，(七，玫，)=Do(七，02)+爿(n，)D1(七，to)+￡r (3) 利用最小二乘准则，求取使得去除～阶自由表面多次波后数据能量最小的最佳予波彳(∞)。最佳子波 的估算是通过逐步迭代实现的，在迭代计算中，将初始爿(∞)定义为零，第…次迭代得到的彳(o’(缈)定义为 6A阳’(∞)，以简化程序设计。我们给出了两种迭代终止的判断条件：一是直接指定迭代次数；二是判断两 次迭代之问得到的校正量6彳(∞)是否已足够小。在计算过程中用到的权值函数％(七，缈)仅是波数的函数， 因此在计算每一个单频之前，首先计算出该函数，以减少计算量。 2．结论 理论模型测试结果表明，基于逆散射级数的2D予波提取方法，不需要任何地下结构信息和子波初始 模型，但要求观测系统为规则观测系统。经降维简化得到的1．5D逆散射级数子波提取方法，同样不需要 地下结构信息和子波初始模型，利用单炮频率一波数域数据，能稳定提取满足自由表面多次波预测精度要求 的子波，且大大降低了计算量。 参考文献 T Luc a1．Source estimationbasedOntheremovaloffirst-order 【11 I，et signature W a1．ABorn—WKBJinversionmethodforacousticreflection 【2】RobertC，et 1)：1559-1567．