一种三维地震偏移成像的块追赶方法技术

本发明专利技术特别涉及一种通过提出的简便算法将得到的地震反射波数据进行计算，然后嵌入到地震数据处理软件平台中形成图像，从而准确地确定油井或其它地下资源最佳位置的方法。本发明专利技术把二维计算中的正定三对角系数矩阵的追赶法推广到三维分块三对角系数矩阵的分裂逼近中，实现了对三维地震数据的快速准确的波场外推处理，把三维有限差分法偏移成像处理技术真正用于勘探工业生产。本发明专利技术的目的在于解决目前所存在的计算量大、精度差、无法避免混沌现象及无法克服方向畸变等方面存在的问题。本发明专利技术在数据处理过程中完全保持波场传播的圆对称性，使最终成像精度得到可靠保证，又由于计算效率比传统方法提高许多，计算成本是工业中可以接受的。

【技术实现步骤摘要】

本专利技术属于勘探地球物理领域中地震偏移成像方面所研究的课题，地震勘探数据是通过现场仪器设备记录到磁性介质并且在地震数据处理软件平台中形成图像，本专利技术提出一种的简便计算方法对地震反射波数据进行偏移成像处理计算，为了解地下介质结构从而准确地确定勘探井位或其它地下资源所处位置提供了一种的切实可行的可靠方法。
技术介绍
利用地震反射波数据资料对地层结构进行成像是地震勘探的核心技术之一。这种称为偏移成像(migration imaging)的技术，很像医学中的CT成像，但是其复杂程度却是难以类比。地震偏移成像的困难在于资料数据量特别大；信号源不确定；介质结构复杂；地震波传播中干扰噪声强大，绕射、屏蔽、吸收等无法掌握的不利因素都对最终成像精度产生坏的影响。与此同时，就成本而言，一口勘探井的工业投入一般达到数千万圆人民币。因此工业上总是要求快速准确的偏移成像成果，作为部署勘探方案、确定井位的依据。因此三维地震偏移成像技术的进步是期待已久的，在地球物理学领域，从理论到工业应用是长期以来受到普遍关注、吸引了大批优秀科学家和工程技术人员来从事这方面的研究。近二、三十年来，借助于计算机技术的迅猛发展，在专门用于地震资料处理大型计算机系统或网络机群上已经建立了不少各具特色的软件平台，这是能源工业勘探中必不可少的基本的常规信号处理工具。但是，任何一种平台都迫切需要提高偏移成像处理的速度和精度，目前还没有一种平台是令人满意的，大批的技术力量在致力于研究开发新的方法。当前，受到普遍重视的是三维隐格式有限差分偏移方法，这种方法的优点是适应地下速度的横向变化，精度高；缺点是计算量太大。1、数据量问题对于10km×10km的叠后数据资料，如果线距和CDP间距为25m，采样间隔为4.00ms，那么在每个时间点上的数据达到400×400＝160000个，在这个时间点上需要求解的方程阶数为160000×160000，对于记录长度为6秒的资料需要求解上述阶数的方程3000-6000次。2、计算效率问题由于数据量大，计算繁杂，通常已知的求解方法无法胜任。现有的计算理论一般只研究一次或几次求解的算法。但是对于一块很小面积的三维地震数据却需要数千次地求解阶数如此巨大的方程。这其实已经不再是效率和成本的问题了，而是有没有工业应用的前景问题了。3、方向畸变问题由于波场外推过程中数千次，使用已算出的数据，如果处理中不能控制好畸变，那么误差的传递积累，足以使得最终成果毫无实际意义。长期以来，人们期待的是快速准确的计算技术，计算技术方面的问题解决了，成像方面的问题也就变的非常简单。现有已知的计算技术中分裂求逆是比较常用的，大致有交替方向隐式(ADI)法也称为沿X-Y方向分裂(Brown，1983等)、基于螺旋边界的降维LU分解(Rickett and Claerbout，1998)。后者是一种非线性方法，计算中无法避免混沌现象，边界处理也有问题。方向畸变问题主要是由于方向分裂算法产生的，由于方向分裂算法十分快捷所以一直被沿用着，但是方向畸变问题却几十年来都没有地解决，在外推过程中由于圆对称性被破坏，无法克服的方向畸变使最终成果被歪曲。在后面的图例和应用例，将给出具体说明。
技术实现思路
1、专利技术目的本专利技术把二维计算中的正定三对角系数矩阵的追赶法推广到三维分块三对角系数矩阵的分裂逼近计算中，其主要目的在于解决现有的地震偏移成像领域中所存在的计算量大、精度差、无法避免混沌现象及无法克服方向畸变等方面存在的问题。2、技术方案本专利技术给出，具体的处理技术是按下列步骤进行a、首先将有限差分法偏移的波场外推方程的系数矩阵做定点的块LU分解分解前矩阵为 其中Tk，Ek为子块矩阵，Tk为三对角矩阵，Ek为对角矩阵，作块两对角分解，就有A＝LLT，其中L为块下三角矩阵，具有如下块两对角形式 块下三角矩阵L的主对角线上为对称正定三对角矩阵Lk，k＝1，2，…，n，主对角线上非零矩阵的左侧除相邻列有非零矩阵之外，其余位置为零矩阵；Ll为L的转置矩阵；b、对称正定三对角矩阵的计算为了计算分块矩阵L，本方法递推地使用Cholesky分解LkLkT=Tk-EkLk-1-TLk-1-1Ek]]>k＝2，3…，n，对上式进行分块逐行递推是一种可以节省存储的方法，因为计算只涉及k-1，k两个参数，因此可以循环使用存储单元；c、稀疏化处理利用上面的Cholesky分解计算Lk时，本方法对Lk-1-1采用了稀疏化处理，由下三角矩阵求逆的递推算法表明逆矩阵中远离主对角线的元素在主对角元占优时，是绝对值迅速变小的，同时Ek的元素的绝对值也是远小于主对角元绝对值的，利用了这种严格单调地绝对值递减规律，大幅度地简化了计算；d、在嵌入到地震数据处理软件平台中的软件模块中，使用上述的块追赶方法，可以通过选定合适的步长等一系列参数得到精度较高的图像。3、优点及效果本专利技术为解决长期困扰勘探地球物理工业界的三维地震偏移成像精度问题，把二维计算中的正定三对角系数矩阵的追赶法计算技术推广到三维分块三对角系数矩阵的分裂逼近中，这是传统的分裂逼近技术中没法实现的一种新方法。本方法实现了对三维地震数据的快速准确的波场外推处理。把三维有限差分法偏移成像处理技术真正用于勘探工业生产。此前由于要么是计算量过于巨大已至于没法承受，要么是波场外推中产生畸变，使有限差分法三维偏移长期停留在理论研究的状态。本专利技术在数据处理过程中完全保持了波场传播的圆对称性，使最终成像精度得到可靠保证，又由于计算效率比传统方法提高许多，计算成本是工业中可以接受的。本专利技术为勘探地球物理工业解决了一个长期的技术难题。实际应用中已经能够证实，所采用的计算方法在辽河油田计算所实施中具有明显的技术效果和显著的经济效益。四附图说明附图1为三维地震偏移方法的分类图；附图2为针对波场外推算子圆对称性畸变所作的各种校正图；附图3为CG偏移的脉冲响应图；附图4为波场外推圆对称性效果图；附图5为另一波场外推圆对称性效果图；附图6为本专利技术的波场外推圆对称性效果图；附图7为0351线偏移前叠加剖面图；附图8为0351线偏移效果图；附图9为另一0351线偏移效果图；附图10为本专利技术的0351线偏移效果图；附图11为三种0351线偏移效果对比图；附图12为0381线偏移效果图；附图13为另一0381线偏移效果图；附图14为本专利技术的0381线偏移效果图；附图15为三种0251线偏移效果对比图。五具体实施例方式本专利技术主要是通过快速准确的计算技术，将地震反射波数据资料进行计算，然后对得到的数据结果通过地震数据处理软件平台中形成图像，从而准确地完成了对地层结构的勘探工作。现有的三维偏移成像方法有很多种，见附图1所示。但对于地震勘探形成地层结构的图像却很难做出来，主要是由于计算量过大、计算中无法避免混沌现象，边界处理也有问题等方面的原因。目前较常采用的计算方法有交替方向隐式(ADI)法和螺旋边界非线性降维LU分解法，螺旋边界非线性降维LU分解法在计算中总出现混沌现象，附图2、3、4、5就采用了交替方向隐式(ADI)法形成的图像。附图2为针对波场外推算子圆对称性畸变所作的各种校正，资料中采用9点、13点、17点、25点和六角形McCllalen变换。附图3为外国软件CG软件偏移的脉冲本文档来自技高网...

【技术保护点】
一种三维地震偏移成像的块追赶方法，其特征在于：该方法主要按下列步骤进行：ａ、首先将有限差分法偏移的波场外推方程的系数矩阵做定点的块ＬＵ分解：分解前矩阵为Ａ＝＊＊＊，其中Ｔ↓［ｋ］，Ｅ↓［Ｋ］为子块矩阵，Ｔ↓［ｋ］为三 对角矩阵，Ｅ↓［ｋ］为对角矩阵，作块两对角分解，就有Ａ＝ＬＬ↑［Ｔ］，其中Ｌ为块下三角矩阵，具有如下块两对角形式：Ｌ＝＊＊＊，块下三角矩阵Ｌ的主对角线上为对称正定三对角矩阵Ｌ↓［ｋ］，ｋ＝１，２，…，ｎ，主对角线上非零矩阵的 左侧除相邻列有非零块矩阵之外，其余位置为零矩阵；Ｌ↑［Ｔ］为Ｌ的转置矩阵；ｂ、对称正定块三对角矩阵的计算：为了计算分块矩阵Ｌ，本方法递推地使用Ｃｈｏｌｅｓｋｙ分解：Ｌ↓［ｋ］Ｌ↓［ｋ］↑［Ｔ］＝Ｔ↓［ｋ］－Ｅ↓［ｋ］ Ｌ↓［ｋ－１］↑［－Ｔ］Ｌ↓［ｋ－１］↑［－１］Ｅ↓［ｋ］ｋ＝２，３…，ｎ，对上式进行分块逐行递推计算。这是一种可以节省存储的方法，因为计算只涉及ｋ－１，ｋ两个参数，因此可以循环使用存储单元；ｃ、稀疏化处理：利用 上面的Ｃｈｏｌｅｓｋｙ分解计算Ｌ↓［ｋ］时，本方法对Ｌ↓［ｋ－１］↑［－１］采用了稀疏化处理，由下三角矩阵求逆的递推算法表明逆矩阵中远离主对角线的元素在主对角元占优时，是绝对值迅速变小的，同时Ｅ↓［ｋ］的元素的绝对值也是远小于主对角元绝对值的，利用了这种严格单调地绝对值递减规律，大幅度地简化了计算；ｄ、分块计算时，实际上是把三维的问题降为二维，可以在其中使用非常快捷常规的追赶法计算；ｅ．在嵌入到地震数据处理软件平台中的软件模块中，使用上述步骤中所描述的块追赶法 计算方法，通过选定合适的步长等一些参数，可得到精度较高的计算结果和较为精确的地下介质结构的图像，当然也可以在精度要求不很严格时获得较高的计算速度从而满足计算效率方面的需要。...

【技术特征摘要】

【专利技术属性】
技术研发人员：石殿祥，
申请(专利权)人：石殿祥，
类型：发明
国别省市：21[中国|辽宁]