论文部分内容阅读
高斯波束偏移作为射线类偏移成像方法的一种,克服了传统的射线类偏移难以处理焦散区成像的缺陷,但具有传统的射线类偏移高效和陡倾角成像的特点,且能够获得接近逆时偏移的成像结果。该方法灵活性强,能够适应多种观测装置及复杂的地表条件。到目前为止,已经提出了针对声波、弹性、各向异性和吸收介质的高斯波束偏移理论。在高斯波束偏移中,由于采用高斯波束的叠加来表示格林函数,所以要进行关于射线参数的积分。在二维情况下该积分是二重的;在三维情况下则是四重的。此外,为了实现成像,在每一个成像点处也需要完成对所有频率的积分。过去,为了提高计算的效率,文献中以损失成像精度的代价采用最速下降近似来降低关于射线参数的积分的维数。而对关于频率的积分的快速实现则鲜有报道。为了解决这个问题,本文通过改变偏移公式中积分的顺序,将两个关于频率的最内层积分看成是总复走时的实部和虚部的二维连续函数来提出一种快速算法。在确定完总复走时的实部和虚部的采样范围和采样间隔之后,在采样点处构建关于两个最内层积分的查找表。最后,根据成像点处实际的总复走时的实部和虚部在查找表中进行插值来获取两个最内层积分的值。这样就避免了在每一个成像点处都进行关于频率的积分运算。同时,由于只是针对关于频率的积分进行处理,该快速算法对于经过及未经最速下降近似的偏移公式同样适用。并行计算环境下的数值算例表明,在二维情况下采用快速算法时,经过最速下降近似的偏移公式相较未经最速下降近似的偏移公式在计算效率方面的提升并不大但却会损失成像的精度。在三维情况下,由于最速下降近似将关于射线参数的积分从四重减少到二重,因此采用快速算法时,经过最速下降近似的偏移公式相较未经最速下降近似的偏移公式在计算效率方面的提升明显。但是对于复杂构造区域,为了满足成像精度的需要,仍需要采用未经最速下降近似的偏移公式。由于上述快速算法同介质的性质无关,可同时适用于弹性及各向异性介质中的高斯波束偏移。针对单道地震数据进行处理的基于高斯波束叠加的深度域偏移也可以采用同上述快速算法类似的策略,所不同的是需要针对每一道地震数据来构建对应两个关于频率的最内层积分的查找表。对于吸收介质中的高斯波束偏移,由于需要引入衰减因子来补偿介质的吸收作用,则将两个关于频率的最内层积分看成是总复走时的实部和虚部及总衰减因子的三维连续函数。此外,本文还讨论了控制束偏移的相关问题。控制束偏移通过在局部倾斜叠加的过程中挑选出最相干的同相轴来压制实际数据中非相干噪声对于成像结果的干扰。作为CGG地球物理公司的核心技术,在公开的文献中未能查阅到其技术细节。本文利用解析道的方法来提高相干性计算的稳定性,同时优化权函数的定义方法以更好地保护相干同相轴。然后,将控制束偏移同上述快速算法相结合来提高整个偏移的计算效率。最后,利用在常倾角成像剖面中局部同相轴的相长部分的斜率同相应的倾角所对应的斜率一致,而相消部分的斜率同相应的倾角所对应的斜率有明显差别的特性,在局部倾角域中对偏移孔径进行优化来提高成像效果。同样地,在计算常倾角成像剖面中局部同相轴的相干性时采用解析道的方法来确保稳定性。