GeoEast系统反射波剩余静校正课件

GeoEast系统反射波剩余静校正GeoEast系统反射波剩余静校正主要内容概述二维地表一致性剩余静校正非地表一致性剩余静校正三维地表一致性剩余静校正主要内容概述二维地表一致性剩余静校正非地表一致性剩余静校正三由于多种因素，一个CMP道集中的各道，经过野外静校正或者初至静校后，仍然存在剩余静校正量，影响CMP叠加质量。因此在CMP叠加之前，我们必须对它进行估算和校正，实现CMP道集的同相叠加。GeoEast处理系统中，对于反射波剩余静校正主要由：剩余静校正量计算（RsStCalculat）三维剩余静校正量调整（StHarmoniz3D）模型道剩应用余时差拾取（SCRsCal3D）三维地表一致性剩余时差分解（SCRsDecom3D）叠前CMP道集时差微调（Trim3D）静校正量应用（StApply）等模块来完成。

由于多种因素，一个CMP道集中的各道，经过野外静剩余静校正模块分三类：三维地表一致性剩余时差计算

ScRscal3D

三维地表一致性剩余时差分解

ScRsDecom3D静校正应用StApply剩余静校正计算

RsStCalculat三维校正量调整

StHarmoniz3D静校正应用StApply叠前道集时差微调Trim3D剩余静校正模块分三类：三维地表一致性剩余时差计算ScRsc输入数据描述输入数据为动校正后的CMP道集数据。输出数据描述反射波剩余校正量，保存在数据库中或者数据表中供后续应用模块使用。输入数据描述主要内容概述二维地表一致性剩余静校正非地表一致性剩余静校正三维地表一致性剩余静校正主要内容概述二维地表一致性剩余静校正非地表一致性剩余静校正三剩余静校正量计算RsStCalculat在GeoEast系统中，剩余静校正计算RsStCalculat是最常用的二维地表一致性剩余静校正模块，该模块利用动校正后的CMP道集数据，计算各炮点和检波点的剩余静校正量，供静校正应用模块（StApply）调用。

剩余静校正量计算RsStCalculat应用实例如下图所示，二维反射波剩余静校正要求输入动校正之后的道集数据，注意数据输入道头字控制为CMP、trace，旗标为CMP，采用自动拾取反射层位选件处理，计算后的剩余校正量保存在数据表中供后续的模块调用。在该模块中，需要关注的主要参数有：

maximumfoldCMP道集中的最高覆盖次数；maximumtracespersource/receivercovers每个炮点或接收点所跨越的最大道数；maximumCMPspersource/receivercovers每个炮点或接收点所跨越的CMP数；statictablefilename本模块生成的剩余静校正量数据表的文件名；statictableversionnumber本模块生成的剩余静校正量数据表的版本号；numberofCMPsfordipscan倾角扫描所需要的CMP数；minimumoffset(m)输入道的最小炮检距，以米为单位；

应用实例maximumfoldCMP道集中的最高覆maximumoffset(m)输入道的最大炮检距，以米为单位；maximumstationsbetweensources在定义观测系统的数据表中，改变炮点位置时，相邻炮号间跨越的最大站数；minimumsourcenumber处理区域中的最小炮号；maximumsourcenumber处理区域中的最大炮号；minimumreceivernumber处理区域中的最小接收点号；maximumreceivernumber处理区域中的最大接收点号；

maximumoffset(m)输入道的最大炮检距，以米该模块应用的数据在完成动校正之后，需要加一个优势频带带通滤波和500ms时窗的动平衡，以便加强模型道，然后是一个基准面校正，这是针对采用CMP参考面进行静校正处理的方式，如果输入数据是在CMP参考面上，在计算剩余静校正量之前，最好加上USCM的基准面校正，也就是说剩余静校正要求在统一(浮动)基准面上进行求取，尤其是对地表起伏比较大的地区更应如此，由于剩余静校正量只是一个地震道的静态时移量，因此即使在统一面上求得的剩余校正量用在CMP面的数据上也没有问题。

该模块应用的数据在完成动校正之后，需要加一个优势频带带通滤波对于该模块二维的应用来说，最关键的参数是：dipdefinedforscanning，即定义自动拾取反射层所用的空间、时间以及扫描倾角范围。其由以下6个参数组成：fromCMP1组参数中的起始CMP号toCMP2组参数中的终止CMP号fromdip1组参数中的起始扫描倾角，以毫秒表示。todip2(ms)组参数中的终止扫描倾角，以毫秒表示。fromtime1组参数中的起始时间，以毫秒为单位。totime2(ms)组参数中的终止时间，以毫秒为单位。对于该模块二维的应用来说，最关键的参数是：dipdefindip1和dip2应在定义的时空范围内倾角最大的反射同相轴上选取。缺省时为相邻12个CMP范围内的最大时差。反射同相轴的时间沿处理方向减小时，倾角为正，反之为负。要处理的第一和最后一个CMP的倾角范围必须定义，时窗缺省时为输入道长。若中间某一个CMP的时窗或倾角范围未定义，则使用上一个CMP的相应参数。CMP的定义可以和输入模块中的CMP定义不一致，在这种情况下，以本模块的定义为准。另外需要注意的是，如果在统一面上进行剩余静校正量的计算，则时窗也要在完成统一面校正后的叠加剖面上选取。dip1和dip2应在定义的时空范围内倾角最大的反射同相轴上一是控制点方式，如下图所示，在全线选择数个控制点，并填写相应的倾角和时窗参数，控制点之间程序自动插值形成全线的计算时窗：二是区域划分法，如下图所示，将整条二维测线划分成若干段，每段使用一组时窗控制参数来定义计算时窗，段间保留一定的空间用于时窗插值。一是控制点方式，如下图所示，在全线选择数个控制点，并填写相应说明和注意事项本模块输入数据顺序为：3D：CMP-LINE；2D：CMP-TRACE在做三维处理时，为了剩余静校正效果好，须将整块三维的cmp线分开，12条cmp线（或更多一些）为一个作业计算剩余静校正量，重复4条cmp，就需要做很多个作业，本模块在处理时，将临时文件记在磁盘上，建议用户在运行本模块前先运行（SysResource）系统资源申请模块，选local参数，将临时文件存放在本节点的/u0/data/临时盘上，以免将所有运行作业的临时文件都放在共享盘data定义下的暂存目录TEMP中。定义自动拾取反射层参数中：1)起始时间应避开剖面切除区。2)倾角扫描范围为相邻12个CMP范围内的最大时差，以毫秒表示。3)反射同相轴的时间沿处理方向减小时，倾角为正，反之为负。本模块的缺省选件是最佳的参数选择值。建议用户尽量使用缺省参数值。使用三维选件时，不能使用userpick选件定义反射层时间参数。说明和注意事项主要内容概述二维地表一致性剩余静校正非地表一致性剩余静校正三维地表一致性剩余静校正主要内容概述二维地表一致性剩余静校正非地表一致性剩余静校正三1.剩余静校正量计算RsStCalculat＋静校正量调整StHarmoniz3DRsStCalculat模块不仅可以用于二维处理同样可以用于三维处理，不同的是在用于三维处理是，模块首先利用动校正后的CMP道集数据，计算并将炮点或接收点的站号、计算该站号剩余静校正量所用的道数、剩余静校正量写成静校正量二进制文件，简称为三字表文件保存在相应三维工区的datatable目录下，然后经静校正量调整模块（StHarmoniz3D）调整后，才能供静校正应用模块（StApply）使用。而且该模块在进行三维处理时，需要将一个较大的三维工区划分为若干个部分，每一部分包含若干条CMP线，而且每一个部分之间彼此略有重叠，依次对每一部分的数据进行RsStCalculat剩余静校正量计算并输出三字表文件，然后由StHarmoniz3D模块一次性完成调整，得到一套剩余校正量供后续应用。也就是说该模块在进行三维处理时，需要由两个模块配合来完成一次三维地表一致性剩余静校正量计算的。1.剩余静校正量计算RsStCalculat＋静校正量调整S1）首先使用RsStCalculat模块进行剩余校正量的求取，并输出三字表文件，如下图所示：输入应用野外静校正（以及必要的初至波剩余静校正）和动校正后的CMP道集，输入关键字控制为CMP、LINE、trace，旗标CMP，输出剩余静校正量三字表文件。1）首先使用RsStCalculat模块进行剩余校正量的求取在三维处理中，除了和二维相同的一些参数需要关注以外，以下参数也需特别关注：3Dprocessoptions三维处理选件：yesstatictablefilename生成的剩余静校正量数据表的文件名以及版本号statictableversionnumber，这两个参数是需要多个RsStCalculat作业中的一个话，不需填写，因为真正的输出是三字表文件；firstlinenumber第一条线号；totalnumberoflines线总数；以上两个参数用于划分三维工区的CMP线范围，划分的不同部分之间需要重叠几条线。startlinenumberfordipanalysis用于倾角分析的起始线号：endlinenumberfordipanalysis用于倾角分析的终止线号：以上两个参数用于定义倾角分析的CMP线范围，这两个参数于时窗参数相配合，完成整个三维工区计算时窗，倾角的定义。maximumdipincrossline(ms)被扫描的最大横向倾角绝对值，该值以横向12条线范围内反射同相轴的毫秒数表示；windowlengthfordipscanincrossline(ms)用于横向倾角扫描分析的时窗长：以毫秒为单位，它必须小于纵向拾取时窗长，缺省值为200；statictripletfilename输出剩余静校正量信息的三字表文件名。在三维处理中，除了和二维相同的一些参数需要关注以外，以下参数2）完成了全工区多个RsStCalculat作业的计算并生成了相应的三字表文件之后，即可使用StHarmoniz3D模块一次性完成剩余校正量的调整，从而得到一套完整的三维工区剩余校正量了。校正量门槛值2）完成了全工区多个RsStCalculat作业的计算并生成具体的参数如下：startsourcenumber起始炮号endsourcenumber终止炮号startreceivernumber起始接收点号endreceivernumber终止接收点号以上四个参数需要和观测系统相对应；statictablefilename剩余静校正量数据表的文件名，用于保存调整后的剩余静校正量。statictableversionnumber剩余静校正量数据表的版本号以上参数定义剩余静校正量保存的位置；statictripletfilename由RsStCalculat模块输出的剩余静校正量信息三字表文件名。最多可填60个。maximumsourcesperline每条炮线上的最大炮点个数。maximumreceiversperline每条接收线上的最大接收点个数。lowerstaticvaluelimit(ms)剩余静校正量下限值，以毫秒为单位；upperstaticvaluelimit(ms)剩余静校正量上限值，以毫秒为单位；本模块对在上、下限值之间的剩余静校正量值按规定方式加以调整，限值之外的剩余静校正量被消掉。

具体的参数如下：静校正应用StApply应用剩余静校正静校正应用StApply应用剩余静校正ＲsStCalculat前后剩余静校正前/后叠加剖面ＲsStCalculat前后剩余静校正前/后叠加剖面功能概述：

利用模型道与地震道的互相关来拾取时差和品质因子（即相似系数），并存放在指定的数据表中，供SCRsDecom3D模块求取剩余静校正量。该模块要求输入CMP道集和相应的叠加道，对于CMP道集和叠加道可以作相应的处理，例如：滤波和振幅补偿等。该模块在形成模型道过程中同时考虑inline和crossline方向的影响，因而是一种真正的三维时差拾取。

2.SCRsCal3D三维地表一致性剩余时差计算

输入数据：动校正后的按inline依次排序的三维/二维CMP道集。模型道数据（包括内部模型道和外部模型道），是三维/二维数据的叠加道或处理后的叠加道。

输出数据：拾取的时差和品质因子，存放在用户指定的数据表中，供SCRsDecom3D模块使用。功能概述：2.SCRsCal3D三维地表一致性剩余时差计参数说明：（一）

datatype 输入数据类型2D 二维数据。当数据类型为二维数据的时候，所有控制点参数的线号信息不用填写，使用缺省值03D 三维数据（缺省）。maximumnumberofCMPsperline

：inline方向测线上最大的CMP个数。地震数据中最大CMP数与最小CMP数的差值。如果此参数为缺省值0，模块将从数据库中读取，所以建议不填写本参数。如果填写必须填写正确数据。当模块运行出现此参数没有填写的错误，说明数据库中没有此参数，这时需要人工填写。本参数仅供3D数据使用。maximumnumberofCMPs：最大的CMP个数。本参数仅供2D数据使用。totalnumberofCMPlines

： 输入数据的inline方向线数。minimumlinenumber

：输入数据的最小inline线号，如果是分块拾取的情况时，所有准备一起分解的拾取数据的最小inline线号必需要相同。建议填写本工区最小线号。参数说明：参数说明（一）minimumcellnumber

：输入数据的最小面元号，如果是分块拾取的情况时，所有准备一起分解的拾取数据的最小面元号必需要相同。建议填写本工区最小面元号，在GeoEast系统中建议使用缺省值。minimumCMPnumber

： 输入数据的最小CMP号，如果是分块拾取的情况时，所有准备一起分解的拾取数据的最小CMP号必需要相同。如果此参数为缺省值0，模块将从数据库中读取，所以建议不填写本参数，如果填写必须填写正确数据。当模块运行出现此参数没有填写的错误，说明数据库中没有此参数，这时需要人工填写。inlinetraceinterval：inline方向道间距，单位为m，默认值为25m。crosslinetraceinterval：crossline方向道间距，单位为m，默认值为50m。参数说明：（二）modeltype：模型类型参数。有两种选择：external： 外部模型internal： 内部模型。当选择内部模型道的时候，会出现inline和crossline两组组合参数，这两组参数是必填的。参数说明参数说明：（二）选择内部模型道时必填的inline和crossline两组组合参数如下：inlineinlinecombinatorialinterpolationtype：建立内部模型道时inline方向组合参数的插值方式，之所以要插值是因为用户只在有限的几个面元处给定参数，而凡是没有给定参数的面元处，程序要通过插值求出其参数值。inlineglobalinterpolation： inline方向全局插值。crosslineglobalinterpolation： crossline方向全局插值。inlinecombinatorialinterpolationangle：inline方向组合的插值角度。单位为度，取值范围为（－90°～90°），缺省值为90°。inlinecombinatorialparametersformodel

：inline方向用来建立模型的参数，每一行为一个时窗的控制点，这个时窗可以有多个控制点，分多行输入。也可以有多个时窗参数。针对某一反射层为主所定义的时窗，时窗名称必须相同。因为当使用内部模型道选件时，inline和crossline方向的组合是以时窗名称相同才组合。如果数据类型为2D，则线号不用填写，使用缺省值0。参数说明：（二）参数说明：（二）选择内部模型道时必填的inline和crossline两组组合参数如下：crosslinecrosslinecombinatorialinterpolationtype

：建立内部模型道时crossline方向组合参数的插值方式。inlineglobalinterpolation

： inline方向全局插值。crosslineglobalinterpolation：crossline方向全局插值。crosslinecombinatorialinterpolationangle

：crossline方向组合的插值角度。单位为度，取值范围为（－90°～90°），缺省值为90°。crosslinecombinatorialparametersformodel：crossline方向用来建立模型的参数。如果数据类型为2D，此参数不用填写，全部使用缺省值0。参数说明：（二）参数说明：（三）maximumfold

：最大覆盖次数。modelfilename：模型数据文件名。当模型类型选择为内部模型时，该数据的各道将作为初始模型道，在程序中对初始模型道进行处理（即内部模型道）。当模型类型选择为外部模型时，该数据的各道将直接作为模型道（即外部模型道）。moveoutdatatablename：时差数据表名，计算出的时差及质量因子输出到该数据表windowinterpolationtype：时窗的插值方式。可选择：inlineglobalinterpolation：inline方向全局插值。crosslineglobalinterpolation：crossline方向全局插值。windowinterpolationangle

： 时窗的插值角度。windowparameters

：时窗参数，每一行为一个时窗的起止时间，这个时窗在不同的CMP和line上可以有不同的起止时间。也可以有多组时窗参数。如果数据类型为2D，则线号不用填写，使用缺省值0。参数说明：（三）参数说明：（三）timeshiftinterpolationtype： 相关时移参数的插值方式。可选择：inlineglobalinterpolation： inline方向全局插值。crosslineglobalinterpolation：crossline方向全局插值。timeshiftinterpolationangle：相关时移参数的插值角度。单位为度，取值范围为（－90°～90°），缺省值为90°。timeshiftlimitparameters： 相关时移限制参数，如图所示。每一行为给定的CMP和line上的相关时移参数限制值。如果数据类型为2D，则线号不用填写，使用缺省值0。

参数说明：（三）参数说明：（四）structureflag： 定义构造项参数标志。可选择：yes： 需要构造项参数no： 不需要构造项参数（缺省）structureinterpolationtype： 构造项参数的插值方式。可选择：inlineglobalinterpolation： inline方向全局插值。crosslineglobalinterpolation： crossline方向全局插值。structureinterpolationangle： 构造项参数的插值角度。单位为度，取值范围为（－90°～90°），缺省值为90°。structureparameters： 构造项的参数，如图6所示。每一行为一个时窗的构造项时间，这个时窗在不同的CMP和line上可以有不同的构造项时间。也可以有多个时窗参数。如果数据类型为2D，则线号不用填写，使用缺省值0。

参数说明：（四）参数说明：（五）startpicking： 起始拾取参数标志。可选择：yes

需要定义起始拾取参数位置（缺省）no

不需要定义起始拾取参数位置startpickingparameters： 起始拾取参数。如图所示。每一行为一个起始拾取的位置，但是同一条线上只能有一个起始CMP号。处理时应尽量选择在叠加剖面上信噪比高的部位作为时差拾取起始点。起始点最好选择靠近数据起始的位置。如果数据类型为2D，则线号不用填写，使用缺省值0。

参数说明：（五）

三维数据的外部模型道作业例子1（inline151-250）一二三四五三维数据的外部模型道作业例子1（inline151-25

三维数据的外部模型道作业例子2（inline251-350）三维数据的外部模型道作业例子2（inline251-35

生成叠后模型道数据的作业例子模型道数据和输入数据基准面必须一致。建议：对叠加数据进行修饰处理，静校正效果好。生成叠后模型道数据的作业例子模型道数据和输入数据基准面必限制和注意事项：1）本模块三维处理的数据与RsStCalculat模块所需的数据排放顺序是不同的，SCRsCal3D要求按in-line线排序，而RsStCalculat则要求按CMP线排序。2）本模块要求输入的模型道数据必须是三维叠后数据，其中内部模型道选件也需要输入模型道数据。当分块拾取的时候，模型数据最好也分块，这样可以节省内存和运行时间。如果分块拾取的时候，建议输出范围不重叠。限制和注意事项：SCRsDecom3D三维地表一致性剩余时差分解SCRsDecom3D模块根据SCRsCal3D模块拾取的时差和品质因子，利用高斯－赛德尔迭代法求取剩余静校正量，然后用求出的炮点剩余静校正量和接收点剩余静校正量建立剩余静校正数据表。最后由StApply模块完成静校正处理。SCRsDecom3D三维地表一致性剩余时差分解SCRsDe

三维地表一致性剩余时差分解的作业例子三维地表一致性剩余时差分解的作业例子作业中对于申请内存所填的参数(值小）作业宽行中统计出的参数总道数=731525+75987564作业中对于申请内存所填的参数(值小）作业宽行中统计出的参数总作业中对于申请内存所填的参数（正确）作业宽行中统计出的参数作业中对于申请内存所填的参数（正确）作业宽行中统计出的参数输出写数据表输出写数据库输出写数据表输出写数据库三、StApply静校正量应用应用数据表剩余静校正量应用数据库剩余静校正量三、StApply静校正量应用应用数据表剩余静校正量应用数据三维地表一致性剩余静校正前（左）/后（右）叠加剖面ScRscal3D三维地表一致性剩余静校正前（左）/后（右）叠加剖面ScRsc3.叠前CMP道集时差微调Trim3D叠前CMP道集时差微调Trim3D主要用于非地表一致性剩余静校正,但其还可通过叠代平均技术得到炮点和检波点的地表一致性剩余静校正量。

模块使用外部模型道stk-srd-rs-rna-pure-all计算地表一致性剩余静校正量，校正量保存在数据表trimstatic.db中，数据的输入控制为line、cmp、trace，旗标为CMP。

3.叠前CMP道集时差微调Trim3D模块使用外部模型道stmaximumstatics最大静校正量允许值。单位为毫秒，范围为0~1000；outputSCstatics是否计算并输出地表一致性静校正量，yes；inputstackfilename叠加文件名。即模型道文件名；

binarystaticsfilename用于存储静校正量的数据表文件；staticsversionnumber静校正版本号。用于向数据库中写入静校正量。缺省值为1。applySCstatics是否应用地表一致性静校正量到CMP数据，yes应用，no不应用；如果选择应用，则作业中该模块之后串接geodskout模块，保存完成静校正之后的地震数据。

maximumstatics最大静校正量允许值。单位为毫CMPvstimewindow定义时窗，每行参数包括三列，共10行。其中：CMPCMP号starttime(ms)时窗的起始时间，单位为毫秒。endtime(ms)时窗的终止时间，单位为毫秒。在计算时，对于小于第一行的CMP值的CMP点，时窗采用第一行的时窗数据，对于大于最后一有效行的CMP值的CMP点，时窗采用最后一有效行的时窗数据，对于CMP值在参数定义范围内的，采用相邻两个时窗线性内插。

CMPvstimewindow定义时窗，每行参数包括主要内容概述二维地表一致性剩余静校正非地表一致性剩余静校正三维地表一致性剩余静校正主要内容概述二维地表一致性剩余静校正非地表一致性剩余静校正三Trim3D叠前道集时差微调

3D非地表一致性剩余静校正是利用经过动校正的CMP道集内的各道与输入的模型道，在用户定义的时窗内做互相关处理，获得相对时移，即非地表一致性静校正量，并对原始