一种三维地震偏移成像的定点分解方法

摘要

一种通过提出的简便算法将得到的地震反射波数据进行计算，然后嵌入到地震数据处理软件平台中形成图像，从而准确地确定油井或其它地下资源最佳位置的方法。本发明采用一种新的分裂逼近方法，实现了一种不完全Cholesky分解(因为完全Cholesky分解的系数矩阵是稠密的不能用于计算)，规定极少数位置上的元素非零，采用误差最小的原则计算出这些非零元，在保证系数矩阵高度稀疏条件下完成逼近。本发明主要解决了现有技术中计算量大、计算中无法避免混沌现象，边界处理不好等方面存在的问题，由于不采用沿方向分裂，避免了方向畸变，保证偏移成像的精度。

主权项

1.一种三维地震偏移成像的定点分解方法，为一种提高地震数据处理精度的方法，其包括首先采用常规手段采集地震勘探数据，然后对采集到的地震勘探数据使用三维地震偏移成像的定点分解方法处理；所述三维地震偏移成像的定点分解方法按下列步骤进行：a、首先将有限差分法偏移的波场外推方程的系数矩阵做定点的块LU分解：分解前矩阵为其中T_k，E_K为子块矩阵，T_k为三对角矩阵，E_k为对角矩阵，作块两对角分解为A＝LL^T+R；其中L为块下三角矩阵，并且在主对角线上有非零元，主对角线上非零矩阵的左侧除相邻列有非零元之外，其余位置为零矩阵；L^T为L的转置矩阵；R是块不完全Cholesky分解的剩余矩阵；零矩阵不需要进行数值计算，计算仅涉及主对角线上的非零元和主对角线两侧的两个非元，得到系数矩阵A的定点LU分解或称定点BIC分解；b、然后对定点BIC分解的残差进行计算：剩余矩阵R＝R_L+R_U为一个在上、下三角部分各有沿对角线方向的一列非零元的矩阵，其中R_L为上三角部分的非零元矩阵，R_T为下三角部分的非零元矩阵，用x＝L^-TL^-1b-Rx得到一种迭代算法x^(k+1)＝L^-TL^-1b-Rx^(k)；c、根据地震勘探数据确定主元优势，进行外推步长计算，选择合适的步长；d、使用上述步骤中得到的A的定点分解式以及所选定的步长，得到精度较高的地下介质结构的图像；e.通过得到的地下介质结构的图像，为部署勘探方案确定钻井井位或勘