【摘 要】
：
不规则地震数据的规则化是数据处理中一个非常重要的环节。本文提出了一种基于迭代加权最小二乘反演的数据规则化方法，该方法在Cauchy模意义下引入了与模型有关的加权算子对数据拟合剩余量进行加权，可以有效避免突发噪音对插值结果的影响。采用局部平面波模型作为先验信息，不仅可以保证反演过程稳定，而且能够对假频数据进行有效的插值。通过预条件共轭梯度法迭代求解，提高收敛速度。对于三维数据体，采用沿时间切片并行化
【机 构】
：
中国石油大学（华东）地球科学与技术学院
【出 处】
：
SPG/SEG2011年国际地球物理会议
论文部分内容阅读
不规则地震数据的规则化是数据处理中一个非常重要的环节。本文提出了一种基于迭代加权最小二乘反演的数据规则化方法，该方法在Cauchy模意义下引入了与模型有关的加权算子对数据拟合剩余量进行加权，可以有效避免突发噪音对插值结果的影响。采用局部平面波模型作为先验信息，不仅可以保证反演过程稳定，而且能够对假频数据进行有效的插值。通过预条件共轭梯度法迭代求解，提高收敛速度。对于三维数据体，采用沿时间切片并行化处理，可以进一步提高三维地震数据规则化的计算效率。理论模型和实际地震数据插值试验验证了本文方法速度快、效果好，具有较高的实用价值。
其他文献
三角洲沉积广泛地分布于渤海湾地区，渤海湾地区是重要的油气产区之一，因此对三角洲砂体的识别对于增储有着重要的意义。我们根据渤海湾地区实际地震资料进行模型简化，设计物理模型参数，制作了相应的物理模型，其中包括高位体系域、湖侵体系域以及低位体系域等三套砂体。在实验室条件下，采集得到物理模拟三维地震数据，对其处理得到三维的地震数据体。砂体识别主要利用谱分解技术，如短时傅里叶变换、小波变换等。本文采用局部属
PML吸收层技术已经被证明是非常有效的边界吸收技术，对体波和面波的吸收都具有非常好的效果，已经被广泛应用于弹性波的数值模拟中。但是在大角度入射即掠射的情况下传统的PML技术还是存在一定的问题，掠射情况下衰减不够，反射系数比较大，离散后会产生比较严重的假反射，降低吸收效果。而掠射是普遍存在的，比如薄片区域、震源位置接近研究区域边缘、大偏移距接收等等。本文基于一阶速度-应力方程，提出了一种卷积完全匹配
基于双程波的逆时偏移是当前精度最高的深度域成像方法，可以对全波场逆时外推，使高角度反射界面的反射波，甚至超过90°反射界面的反射波进行精确成像。通过研究，基于有限差分方法提出了精确且高效的差分格式逆时波场外推算法。针对多次反射波严重和构造十分复杂时，由互相关成像条件引入的低频噪音通过优化Laplace算子进行压制。以CPU/GPU作为计算核心，将计算量最大的波场逆时外推通过GPU实现。并利用随机速
陆相薄互层砂岩储层厚度一般明显小于调谐厚度，对于这类储层地震预测的关键在于砂体本身的识别，而非单砂体顶底界面的分辨，所以地震识别极限厚度的研究极其重要。然而，由于影响地震识别能力的因素众多，对于地震识别极限厚度及其影响因素的研究到目前为止仍然没有取得令人满意的结果。本文通过稳定互层模型的设计和参数优选，引入砂地比的概念，并通过褶积模型的解释和对比及进一步的理论推导，给出了影响地震识别储层极限厚度的
随着地震勘探的发展，越来越需要从叠前数据中获得带有方位信息的成像数据。传统方位分析考虑的地面炮检方位并不代表地下波场的真实传播方位；常用的炮检距域成像会引起入射角的混叠，在强横向变速情况下甚至还产生许多假象。适应方位各向异性介质的局部角度域叠前偏移成像方法可以很好的解决上述问题。本文基于设计的三维HTI模型，利用拟声波方程正演合成的宽方位数据，系统地测试了方位各向异性局部角度域叠前时间偏移成像算法
本文提出两种全波形反演的梯度预处理格式。第一种格式是基于时间积分残差波场建立最小二乘目标泛函，这在一定程度上解决了在传统的梯度公式中的时间导数高通滤波的影响，这常常影响反演方法的收敛效率，利用时间积分波场建立的目标泛函可以极大加快收敛率；另一种格式是受扩展成像条件的启发，通过定义广义局部地下偏移距的Hessian 矩阵，能给传统的Hessian 矩阵引入关于模型的冗余信息。本文比较了传统的炮域照明
地球实际并非完全弹性介质，地层吸收致使地震波在地层中传播存在能量衰减及速度频散。本文采用基于Futterman模型的分步傅里叶（SSF）波场延拓方法，实现了单程波动方程非零偏移距衰减介质地震波传播数值模拟。针对基于检波点波场下延原理进行炮记录计算效率低的问题，我们根据反射波场物理机制采用反射波正向传播理论，提高了计算的效率。通过数值模拟对比了单程波动方程和双程波动方程的特点，并探索衰减介质中地震波
基于构建虚同相轴来估计层间多次波的技术是克希霍夫积分表示定理的一个延伸发展。虚地震同相轴由Ikelle提出，它的存在使我们能够利用表面的散射点来估计层间多次波。它巧妙地将散射点从地下移到了表面，并且实现过程仅仅需要通过一个互相关和后续的卷积计算即可完成。在估计出层间多次波之后，由于它与实际地震记录中的多次波在振幅和频率上存在偏差，时间上也可能产生一个小的移动，因此我们采用了多道的L1范数匹配算法实
保持振幅一直受到广泛关注并取得很大进步，因此能基本满足构造勘探的要求，但是对于岩性勘探却依然存在陷阱，主要表现在：动校正叠加存在斜层反射点分散的问题，其共中心点CMP道集与一段地层的反射波相对应；叠前偏移存在混合作用及平均作用，共成像点CIP道集里包含相邻反射点振幅的影响；动校正和偏移过程中都存在拉伸作用；CRP沿层叠加是针对目标地层反射波实施的专门处理，通过沿层追踪反射波（TRALayer）技术
同态反褶积由于受到噪声及反射系数非随机的影响，计算出的子波相位稳定差，难以在实际生产中得到应用。通过结合小时窗、多道数据计算的平均复赛谱和平均自相关函数的复赛谱来计算混合相位子波，该方法增加了子波相位的稳定性，使同态反褶积的实用性大大增强。在实际资料处理过程中，又使用了指数衰减法将地震子波衰减为小相位，并对求出的子波进行能量定位，选择能量集中的部分，对计算出的子波进行窗口整形，从而减少子波旁瓣的影