第8章CRS叠加.docVIP

第八章 共反射面叠加 §8.1 二维共反射面叠加 在常规处理中，时间域成像和叠后偏移的效果不是很理想，对复杂构造的成像更是如此。于是一些新的时间域成像技术和叠前偏移方法应运而生，共反射面（CRS）叠加便是其中的一种。CRS叠加是一种不依赖于速度信息，可以直接由多次覆盖反射数据得到零炮检距（ZO）剖面的叠加成像技术。对确定CRS叠加面的三个参数（ZO射线的出射角，以及与ZO射线有关的两个波前曲率半径和）进行优化可使CRS走时面最佳地拟合反射同相轴。CRS叠加不仅能够改进模拟ZO剖面，提高深层的信噪比（S/N），而且给出了可用于反演速度场的三参数剖面。模型数据的试算和实际资料的处理验证了该方法的有效性和实用性。 一．概述 CRS叠加的思想最早来源于Hubral教授。首先给出了与ZO射线有关的法向入射点波波前曲率半径的概念（Hubral, 1983），以后相继提出了基于包括ZO射线在地面的出射角和在内的两参数优化的共反射元（CRE）叠加法（Cruz et al, 2000）和基于包括和与ZO射线有关的两个波前曲率半径-和法向波波前曲率半径在内的三参数优化的CRS叠加法（Muller er al, 1998; Hubral et al, 1999; Garabito, 2001; J?ger et al, 2001; Bergler, 2002; zhang et al, 2002）。CRS叠加是一仅依赖于近地表速度而与宏观速度模型无关的地震成像方法。它的理论基础是几何地震学，考虑了反射层的局部特征和第一菲涅耳带内的全部反射。理论分析表明这一方法在提高信噪比和成像精度方面优于其它叠加技术、偏移到零炮检距（MZO）、NMO/DMO/叠加和克希霍夫叠前深度偏移（PreSDM）技术，从而为地震叠加和成像技术提供了一种新方法。CRS叠加解决深层复杂构造和岩性地震反射成像的能力使其成为提高深层地震资料质量的一种有效手段。高质量的三参数剖面可用于深层速度反演。 二．方法原理 图7-1a（Hubral et al, 1999）是MZO原理图，它是基于ZO等时线叠加原理，只在反射点处准确，一离开反射点误差明显增大，同相性显著变差；图7-1b（Hubral et al, 1999）是Kirchhoff PreSDM原理图，它是基于绕射曲线叠加原理，在零炮检距处最准确，在近炮检距处较准确，在远炮检距处误差变大；而CRS叠加原理可见图7-1c（Hubral et al, 1999）。与MZO和Kirchhoff PreSDM相比较，CRS叠加考虑了反射层的局部特征和第一菲涅耳带内的全部反射，反射波的同相性好，有效利用了多次复盖反射数据，信噪比明显提高。 利用射线理论和二阶泰勒展式，可得出以地震三参数和表示的双曲近似和抛物近似的时距方程。在以中心点和半炮检距h建立的坐标系中，CRS面的双曲走时近似公式为 （7-1） CRS面的抛物走时近似公式为 （7-2） 图7-1a 基于零炮检距叠加的MZO原理图 图7-1b 基于绕射曲线叠加的Kirchhoff PreSDM原理图 图1c CRS叠加原理图 其中，是近地表速度；地震三参数和的含义如图7-2所示。是ZO射线在地表的出射角；是法向入射点波波前曲率半径，法向入射点波波前对应于反射界面上点源产生的波前；为法向波波前曲率半径，法向波波前对应于爆炸反射面产生的波前。由于地震资料的时距关系更符合双曲规律，在此我们使用了CRS面的双曲走时近似公式。 基于方程（7-1），对ZO剖面上的每一点，都须确定最佳三参数和。由这些最佳参数，公式（7-1）给出的走时面能最好地拟合反射同相轴。求取这些参数的一种方法就是测试所有的三参数，从中选取具有最大相干值的三参数，当然，测试无穷个三参数是不实际的，它必须局限到一有限三维网格，然而，这仍需昂贵的计算成本。因此，必须寻找一种更有效的方法来求取最优三参数，即以最短计算时间求取相干值的最大绝对值，该相干值是三个独立变量和的函数。 图7-2 法向入射点波和法向波形成图 其中， 和分别表示ZO射线在地面的出射角、法向入射点波波前、法向波波前、入射波波前、透射波波前和反射层的曲率半径，表示层速度 这是一典型的全局优化问题，一般说来，求取全局极值是一很困难的事情，而且所有这些最优化算法的成功应用都大大取决于所考虑的函数性质。只有给出的初始参数非常接近于最佳三参数，最优化算法才是成功的。否则，最优化算法很可能给出一局部极大值。因此，为更好地使用最优化算法，必须确定一初始三参数()作为最优化算法的起点，这就是我们提出参数多级优化的原因。 易于