专利名称：块追赶法用于三维地震偏移计算的地震数据处理系统的制作方法
技术领域：
本发明涉及地震勘探数据处理技术，特别是“三维地震偏移”的块追赶技术。
地震勘探工作主要由两大工作系统组成，一是地震数据采集系统，二是地震数据处理系统。其中，地震数据处理系统包括大型运算设备、数据处理程序、适用于相应程序的运算方法等。在数据处理程序中，一个重要的组成部分是“三维地震偏移”，即通过一定的技术手段来最大限度地消除在地面上测试得到的用来描述工区地质图像的三维地震数据与工区实际地质结构之间的偏差，以便得到能够更为准确地反映工区实际地质结构的勘探结果图像。到目前为止，用来实现“三维地震偏移”操作的主要途径是“波场外推”技术，即通过一定的运算方法，使得在地面上取得的地震波场数据，转化成地下震点处波场数据，再利用该数据制得反映上述工区实际地质结构情况的勘探结果图像。现有的利用“波场外推”技术来实现“三维地震偏移”操作的计算方法较多，主要是“ADI交替方向隐式法”和尚未推广应用的“降维LU分解逼近法”。前者的缺欠是计算误差明显地与地下反射层倾角、方位及速度变化有关，这不仅增加了运算难度，也使波场外推算子的圆对称性不能保证，因而计算速度与计算精度都不够理想；后者对“波场外推”方程的算子矩阵采用非线性的方法作LU分解、而后作降维处理，即把远离对角线的元素取为零，再经两次递归多项式除法后完成求解计算。这种方法的局限性首先在于引入非线性方法后，使计算过程出现混沌状态，造成计算困难。其次，到目前为止，还没有足够的理论依据来证明该方法能够保证外推算子的圆对称性，缺少定量分析，计算的可靠性尚未得到证实。
本发明的目的是为人们提供一种采用块追赶法实施“三维地震偏移”操作计算任务的地震数据处理系统，以便使“三维地震偏移”操作的计算过程更加方便、快捷，并提高运算精度，从而进一步提高地震勘探工作的总体质量水平。
本发明的目的是这样实现的，它包括地震数据处理系统中的“三维地震偏移”操作程序，实施该操作采用“波场外推”技术，其特征在于“波场外推”的计算方法为基于块两对角LU分解的块追赶法。
本发明实施后与现有技术相比具有以下优点1、相对于现有的降维LU分解逼近法而言，由于其运算中不引入非线性方法，因而运算过程筒单、明了、易于操作。其次，由于运算中不必进行降维处理，因而能较好地保持波场外推算子的圆对称性，提高波场外推算子的逼近精度。2、与现有的ADI方法相比，由于其运算过程与反射层倾角、方位、及其速度变化等复杂的附加因素无关，因而，一方面可大幅度提高运算速度，另一方面也有利于保证波场外推算子的圆对称性，进而提高运算精度。
下面结合实施例对本发明做进一步详细地说明。
利用本发明提供的地震数据处理系统，进行地震偏移操作时，也要通过“波场外推”途径，这与已有技术是相同的，而其“波场外推”的计算方法，则采用块追赶法，这与现有的降维LU分解逼近方法有着明显的区别，这种区别在于它是不降维的，不引入非线性方法。块追赶法的计算步骤及其算式如下三维波场外推方程系数矩阵为如下的块三对角阵
其中Ti，i＝1，2...，n为三对角矩阵，Ei，i＝2，3...，ni为对角阵。对于上式采用如下的块两对角形式的块LU分解，
显然，2式中的下三角阵Lk，K＝1，2，...，n是十分容易求逆的，至于Lk的计算则可以采用下面的cholesky分解LkLTk＝Tk-EkLk-1-TLk-1-1Ek(3)递进地进行。
在使用cholesky分解式3求Lk时，要对LK-1-1采用稀疏处理，尽可能地简化计算，这里对LK的计算或者说对LK-1-1中绝对值较小元素的舍弃是有要根据的，这个根据就是它们的确定是严格单调地绝对值递减的，并且这些小量的量值可以由波场外算子计算公式中的速度、道向距、cdp间距、时间采样率、外推步长等作出准确估算。
由块两对角形式的分解式，通过两次三角矩阵的求逆立刻得到了波场外推计算的块追赶法。
在计算中虽然我们作了必要的近似处理，但是这种小范围的近似处理无论在理论上，还是在实际计算中都是可行的。
权利要求
1.块追赶法用于三维地震偏移计算的地震数据处理系统，包括该系统中的“三维地震偏移”操作程序，实施该操作采用“波场外推”技术，其特征在于“波场外推”的计算方法为基于块两对角LU分解的块追赶法。
全文摘要
本发明公开了一种块追赶法用于“三维地震偏移”计算的地震数据处理系统,包它括该系统中的“三维地震偏移”操作程序,实施该操作采用“波场外推”技术,其特征在于“波场外推”的计算方法为基于块两对角LU分解的块追赶法。本发明实现后,可使其“三维地震偏移”操作的计算过程更加方便、快捷、并提高运算精度,从而进一步提高地震勘探工作的总体质量水平。
文档编号G01V1/28GK1255643SQ9912256
公开日2000年6月7日 申请日期1999年11月23日 优先权日1999年11月23日
发明者石殿祥 申请人:中油辽河油田公司勘探开发研究院