《GeoEast系统反射波剩余静校正(2).ppt》由会员分享，可在线阅读，更多相关《GeoEast系统反射波剩余静校正(2).ppt（56页珍藏版）》请在三一办公上搜索。
1、GeoEast系统反射波剩余静校正,主要内容,二维地表一致性剩余静校正,由于多种因素，一个CMP道集中的各道，经过野外静校正或者初至静校后，仍然存在剩余静校正量，影响CMP叠加质量。因此在CMP叠加之前，我们必须对它进行估算和校正，实现CMP道集的同相叠加。GeoEast 处理系统中，对于反射波剩余静校正主要由：剩余静校正量计算（RsStCalculat）三维剩余静校正量调整（StHarmoniz3D）模型道剩应用余时差拾取（SCRsCal3D）三维地表一致性剩余时差分解（SCRsDecom3D）叠前CMP道集时差微调（Trim3D）静校正量应用（StApply）等模块来完成。,剩余静校正模块
2、分三类：,三维地表一致性剩余时差计算 ScRscal3D 三维地表一致性剩余时差分解 ScRsDecom3D静校正应用 StApply,剩余静校正计算 RsStCalculat三维校正量调整 StHarmoniz3D静校正应用 StApply,叠前道集时差微调Trim3D,输入数据描述 输入数据为动校正后的CMP道集数据。输出数据描述 反射波剩余校正量，保存在数据库中或者数据表中供后续应用模块使用。,主要内容,概述,二维地表一致性剩余静校正,剩余静校正量计算RsStCalculat在GeoEast系统中，剩余静校正计算RsStCalculat是最常用的二维地表一致性剩余静校正模块，该模块利用动
3、校正后的CMP道集数据，计算各炮点和检波点的剩余静校正量，供静校正应用模块（StApply）调用。,应用实例如下图所示，二维反射波剩余静校正要求输入动校正之后的道集数据，注意数据输入道头字控制为CMP、trace，旗标为CMP，采用自动拾取反射层位选件处理，计算后的剩余校正量保存在数据表中供后续的模块调用。在该模块中，需要关注的主要参数有：,maximum fold CMP 道集中的最高覆盖次数；maximum traces per source/receiver covers 每个炮点或接收点所跨越的最大道数；maximum CMPs per source/receiver covers 每
4、个炮点或接收点所跨越的CMP数；static table filename 本模块生成的剩余静校正量数据表的文件名；static table version number 本模块生成的剩余静校正量数据表的版本号；number of CMPs for dip scan 倾角扫描所需要的CMP数；minimum offset(m)输入道的最小炮检距，以米为单位；,maximum offset(m)输入道的最大炮检距，以米为单位；maximum stations between sources 在定义观测系统的数据表中，改变炮点位置时，相邻炮号间跨越的最大站数；minimum source numb
5、er 处理区域中的最小炮号；maximum source number 处理区域中的最大炮号；minimum receiver number 处理区域中的最小接收点号；maximum receiver number 处理区域中的最大接收点号；,该模块应用的数据在完成动校正之后，需要加一个优势频带带通滤波和500ms时窗的动平衡，以便加强模型道，然后是一个基准面校正，这是针对采用CMP参考面进行静校正处理的方式，如果输入数据是在CMP参考面上，在计算剩余静校正量之前，最好加上USCM的基准面校正，也就是说剩余静校正要求在统一(浮动)基准面上进行求取，尤其是对地表起伏比较大的地区更应如此，由于剩余
6、静校正量只是一个地震道的静态时移量，因此即使在统一面上求得的剩余校正量用在CMP面的数据上也没有问题。,对于该模块二维的应用来说，最关键的参数是：dip defined for scanning，即定义自动拾取反射层所用的空间、时间以及扫描倾角范围。其由以下6个参数组成：from CMP1 组参数中的起始CMP号 to CMP2 组参数中的终止CMP号 from dip1 组参数中的起始扫描倾角，以毫秒表示。to dip2(ms)组参数中的终止扫描倾角，以毫秒表示。from time1 组参数中的起始时间，以毫秒为单位。to time2(ms)组参数中的终止时间，以毫秒为单位。,dip1和di
7、p2应在定义的时空范围内倾角最大的反射同相轴上选取。缺省时为相邻12个CMP范围内的最大时差。反射同相轴的时间沿处理方向减小时，倾角为正，反之为负。要处理的第一和最后一个CMP的倾角范围必须定义，时窗缺省时为输入道长。若中间某一个CMP的时窗或倾角范围未定义，则使用上一个CMP的相应参数。CMP的定义可以和输入模块中的CMP定义不一致，在这种情况下，以本模块的定义为准。另外需要注意的是，如果在统一面上进行剩余静校正量的计算，则时窗也要在完成统一面校正后的叠加剖面上选取。,一是控制点方式，如下图所示，在全线选择数个控制点，并填写相应的倾角和时窗参数，控制点之间程序自动插值形成全线的计算时窗：,二
8、是区域划分法，如下图所示，将整条二维测线划分成若干段，每段使用一组时窗控制参数来定义计算时窗，段间保留一定的空间用于时窗插值。,说明和注意事项本模块输入数据顺序为：3D：CMP-LINE；2D：CMP-TRACE在做三维处理时，为了剩余静校正效果好，须将整块三维的cmp线分开，12条cmp线（或更多一些）为一个作业计算剩余静校正量，重复4条cmp，就需要做很多个作业，本模块在处理时，将临时文件记在磁盘上，建议用户在运行本模块前先运行（SysResource）系统资源申请模块，选local参数，将临时文件存放在本节点的/u0/data/临时盘上，以免将所有运行作业的临时文件都放在共享盘data定
9、义下的暂存目录TEMP中。定义自动拾取反射层参数中：1)起始时间应避开剖面切除区。2)倾角扫描范围为相邻12个CMP范围内的最大时差，以毫秒表示。3)反射同相轴的时间沿处理方向减小时，倾角为正，反之为负。本模块的缺省选件是最佳的参数选择值。建议用户尽量使用缺省参数值。使用三维选件时，不能使用user pick选件定义反射层时间参数。,主要内容,概述,二维地表一致性剩余静校正,三维地表一致性剩余静校正,1.剩余静校正量计算RsStCalculat静校正量调整StHarmoniz3DRsStCalculat模块不仅可以用于二维处理同样可以用于三维处理，不同的是在用于三维处理是，模块首先利用动校正后
10、的CMP道集数据，计算并将炮点或接收点的站号、计算该站号剩余静校正量所用的道数、剩余静校正量写成静校正量二进制文件，简称为三字表文件保存在相应三维工区的datatable目录下，然后经静校正量调整模块（StHarmoniz3D）调整后，才能供静校正应用模块（StApply）使用。而且该模块在进行三维处理时，需要将一个较大的三维工区划分为若干个部分，每一部分包含若干条CMP线，而且每一个部分之间彼此略有重叠，依次对每一部分的数据进行RsStCalculat剩余静校正量计算并输出三字表文件，然后由StHarmoniz3D模块一次性完成调整，得到一套剩余校正量供后续应用。也就是说该模块在进行三维处理
11、时，需要由两个模块配合来完成一次三维地表一致性剩余静校正量计算的。,1）首先使用RsStCalculat模块进行剩余校正量的求取，并输出三字表文件，如下图所示：输入应用野外静校正（以及必要的初至波剩余静校正）和动校正后的CMP道集，输入关键字控制为CMP、LINE、trace，旗标CMP，输出剩余静校正量三字表文件。,在三维处理中，除了和二维相同的一些参数需要关注以外，以下参数也需特别关注：3D process options 三维处理选件：yesstatic table filename生成的剩余静校正量数据表的文件名以及版本号static table version number，这两个参
12、数是需要多个RsStCalculat作业中的一个话，不需填写，因为真正的输出是三字表文件；first line number第一条线号；total number of lines线总数；以上两个参数用于划分三维工区的CMP线范围，划分的不同部分之间需要重叠几条线。start line number for dip analysis用于倾角分析的起始线号：end line number for dip analysis用于倾角分析的终止线号：以上两个参数用于定义倾角分析的CMP线范围，这两个参数于时窗参数相配合，完成整个三维工区计算时窗，倾角的定义。maximum dip in crosslin
13、e(ms)被扫描的最大横向倾角绝对值，该值以横向12条线范围内反射同相轴的毫秒数表示；window length for dip scan in crossline(ms)用于横向倾角扫描分析的时窗长：以毫秒为单位，它必须小于纵向拾取时窗长，缺省值为200；static triplet filename输出剩余静校正量信息的三字表文件名。,2）完成了全工区多个RsStCalculat作业的计算并生成了相应的三字表文件之后，即可使用StHarmoniz3D模块一次性完成剩余校正量的调整，从而得到一套完整的三维工区剩余校正量了。,校正量门槛值,具体的参数如下：start source number
14、 起始炮号 end source number 终止炮号 start receiver number 起始接收点号 end receiver number 终止接收点号 以上四个参数需要和观测系统相对应；static table filename 剩余静校正量数据表的文件名，用于保存调整后的剩余静校正量。static table version number 剩余静校正量数据表的版本号 以上参数定义剩余静校正量保存的位置；static triplet filename 由RsStCalculat模块输出的剩余静校正量信息三字表文件名。最多可填60个。maximum sources per li
15、ne 每条炮线上的最大炮点个数。maximum receivers per line 每条接收线上的最大接收点个数。lower static value limit(ms)剩余静校正量下限值，以毫秒为单位；upper static value limit(ms)剩余静校正量上限值，以毫秒为单位；本模块对在上、下限值之间的剩余静校正量值按规定方式加以调整，限值之外的剩余静校正量被消掉。,静校正应用StApply,sStCalculat,前,后,剩余静校正前/后叠加剖面,功能概述：利用模型道与地震道的互相关来拾取时差和品质因子（即相似系数），并存放在指定的数据表中，供SCRsDecom3D模块求取
16、剩余静校正量。该模块要求输入CMP道集和相应的叠加道，对于CMP道集和叠加道可以作相应的处理，例如：滤波和振幅补偿等。该模块在形成模型道过程中同时考虑inline和crossline方向的影响，因而是一种真正的三维时差拾取。,2.SCRsCal3D 三维地表一致性剩余时差计算,输入数据：动校正后的按inline依次排序的三维/二维CMP道集。模型道数据（包括内部模型道和外部模型道），是三维/二维数据的叠加道或处理后的叠加道。输出数据：拾取的时差和品质因子，存放在用户指定的数据表中，供SCRsDecom3D模块使用。,参数说明：（一）data type 输入数据类型2D 二维数据。当数据类型为二
17、维数据的时候，所有控制点参数的线号信息不用填写，使用缺省值0 3D 三维数据（缺省）。maximum number of CMPs per line：inline方向测线上最大的CMP个数。地 震数据中最大CMP数与最小CMP数的差值。如果此参数为缺省值0，模块将从数据库中读取，所以建议不填写本参数。如果填写必须填写正确数据。当模块运行出现此参数没有填写的错误，说明数据库中没有此参数，这时需要人工填写。本参数仅供3D数据使用。maximum number of CMPs：最大的CMP个数。本参数仅供2D数据使用。total number of CMP lines：输入数据的inline方向线数
18、。minimum line number：输入数据的最小inline线号，如果是分块拾取的情况时，所有准备一起分解的拾取数据的最小inline线号必需要相同。建议填写本工区最小线号。,参数说明（一）minimum cell number：输入数据的最小面元号，如果是分块拾取的情况时，所有准备一起分解的拾取数据的最小面元号必需要相同。建议填写本工区最小面元号，在GeoEast系统中建议使用缺省值。minimum CMP number：输入数据的最小CMP号，如果是分块拾取的情况时，所有准备一起分解的拾取数据的最小CMP号必需要相同。如果此参数为缺省值0，模块将从数据库中读取，所以建议不填写本参数
19、，如果填写必须填写正确数据。当模块运行出现此参数没有填写的错误，说明数据库中没有此参数，这时需要人工填写。inline trace interval：inline方向道间距，单位为m，默认值为25m。crossline trace interval：crossline方向道间距，单位为m，默认值为50m。,参数说明：（二）model type：模型类型参数。有两种选择：external：外部模型internal：内部模型。当选择内部模型道的时候，会出现inline和crossline两组组合参数，这两组参数是必填的。,参数说明：（二）选择内部模型道时必填的inline和crossline两组组
20、合参数如下：inlineinline combinatorial interpolation type：建立内部模型道时inline方向组合参数的插值方式，之所以要插值是因为用户只在有限的几个面元处给定参数，而凡是没有给定参数的面元处，程序要通过插值求出其参数值。inline global interpolation：inline方向全局插值。crossline global interpolation：crossline方向全局插值。inline combinatorial interpolation angle：inline方向组合的插值角度。单位为度，取值范围为（9090），缺省值为90
21、。inline combinatorial parameters for model：inline方向用来建立模型的参数，每一行为一个时窗的控制点，这个时窗可以有多个控制点，分多行输入。也可以有多个时窗参数。针对某一反射层为主所定义的时窗，时窗名称必须相同。因为当使用内部模型道选件时，inline和crossline方向的组合是以时窗名称相同才组合。如果数据类型为2D，则线号不用填写，使用缺省值0。,参数说明：（二）选择内部模型道时必填的inline和crossline两组组合参数如下：crosslinecrossline combinatorial interpolation type：建立
22、内部模型道时crossline方向组合参数的插值方式。inline global interpolation：inline方向全局插值。crossline global interpolation：crossline方向全局插值。crossline combinatorial interpolation angle：crossline方向组合的插值角度。单位为度，取值范围为（9090），缺省值为90。crossline combinatorial parameters for model：crossline方向用来建立模型的参数。如果数据类型为2D，此参数不用填写，全部使用缺省值0。,参数说明
23、：（三）maximum fold：最大覆盖次数。model filename：模型数据文件名。当模型类型选择为内部模型时，该数据的各道将作为初始模型道，在程序中对初始模型道进行处理（即内部模型道）。当模型类型选择为外部模型时，该数据的各道将直接作为模型道（即外部模型道）。moveout data table name：时差数据表名，计算出的时差及质量因子输出到该数据表window interpolation type：时窗的插值方式。可选择：inline global interpolation：inline方向全局插值。crossline global interpolation：cross
24、line方向全局插值。window interpolation angle：时窗的插值角度。window parameters：时窗参数，每一行为一个时窗的起止时间，这个时窗在不同的CMP和line上可以有不同的起止时间。也可以有多组时窗参数。如果数据类型为2D，则线号不用填写，使用缺省值0。,参数说明：（三）time shift interpolation type：相关时移参数的插值方式。可选择：inline global interpolation：inline方向全局插值。crossline global interpolation：crossline方向全局插值。time shift
25、 interpolation angle：相关时移参数的插值角度。单位为度，取值范围为（9090），缺省值为90。time shift limit parameters：相关时移限制参数，如图所示。每一行为给定的CMP和line上的相关时移参数限制值。如果数据类型为2D，则线号不用填写，使用缺省值0。,参数说明：（四）structure flag：定义构造项参数标志。可选择：yes：需要构造项参数no：不需要构造项参数（缺省）structure interpolation type：构造项参数的插值方式。可选择：inline global interpolation：inline方向全局插值。
26、crossline global interpolation：crossline方向全局插值。structure interpolation angle：构造项参数的插值角度。单位为度，取值范围为（9090），缺省值为90。structure parameters：构造项的参数，如图6所示。每一行为一个时窗的构造项时间，这个时窗在不同的CMP和line上可以有不同的构造项时间。也可以有多个时窗参数。如果数据类型为2D，则线号不用填写，使用缺省值0。,参数说明：（五）start picking：起始拾取参数标志。可选择：yes 需要定义起始拾取参数位置（缺省）no 不需要定义起始拾取参数位置st
27、art picking parameters：起始拾取参数。如图所示。每一行为一个起始拾取的位置，但是同一条线上只能有一个起始CMP号。处理时应尽量选择在叠加剖面上信噪比高的部位作为时差拾取起始点。起始点最好选择靠近数据起始的位置。如果数据类型为2D，则线号不用填写，使用缺省值0。,三维数据的外部模型道作业例子1（inline151-250）,一,二,三,四,五,三维数据的外部模型道作业例子2（inline251-350）,生成叠后模型道数据的作业例子,模型道数据和输入数据基准面必须一致。建议：对叠加数据进行修饰处理，静校正效果好。,限制和注意事项：1）本模块三维处理的数据与RsStCalcu
28、lat模块所需的数据排放顺序是不同的，SCRsCal3D要求按in-line线排序，而RsStCalculat则要求按CMP线排序。2）本模块要求输入的模型道数据必须是三维叠后数据，其中内部模型道选件也需要输入模型道数据。当分块拾取的时候，模型数据最好也分块，这样可以节省内存和运行时间。如果分块拾取的时候，建议输出范围不重叠。,SCRsDecom3D三维地表一致性剩余时差分解,SCRsDecom3D模块根据SCRsCal3D模块拾取的时差和品质因子，利用高斯赛德尔迭代法求取剩余静校正量，然后用求出的炮点剩余静校正量和接收点剩余静校正量建立剩余静校正数据表。最后由StApply模块完成静校正处理
29、。,三维地表一致性剩余时差分解的作业例子,作业中对于申请内存所填的参数(值小）,作业宽行中统计出的参数,总道数=731525+75987564,作业中对于申请内存所填的参数（正确）,作业宽行中统计出的参数,三、StApply静校正量应用,3.叠前CMP道集时差微调Trim3D叠前CMP道集时差微调Trim3D主要用于非地表一致性剩余静校正,但其还可通过叠代平均技术得到炮点和检波点的地表一致性剩余静校正量。,模块使用外部模型道stk-srd-rs-rna-pure-all计算地表一致性剩余静校正量，校正量保存在数据表trimstatic.db中，数据的输入控制为line、cmp、trace，旗标
30、为CMP。,maximum statics 最大静校正量允许值。单位为毫秒，范围为01000；output SC statics 是否计算并输出地表一致性静校正量，yes；input stack filename 叠加文件名。即模型道文件名；,binary statics filename 用于存储静校正量的数据表文件；statics version number静校正版本号。用于向数据库中写入静校正量。缺省值为1。apply SC statics 是否应用地表一致性静校正量到CMP数据，yes 应用，no 不应用；如果选择应用，则作业中该模块之后串接geodskout模块，保存完成静校正之后
31、的地震数据。,CMP vs time window 定义时窗，每行参数包括三列，共10行。其中：CMP CMP号 start time(ms)时窗的起始时间，单位为毫秒。end time(ms)时窗的终止时间，单位为毫秒。在计算时，对于小于第一行的CMP值的CMP点，时窗采用第一行的时窗数据，对于大于最后一有效行的CMP值的CMP点，时窗采用最后一有效行的时窗数据，对于CMP值在参数定义范围内的，采用相邻两个时窗线性内插。,主要内容,概述,二维地表一致性剩余静校正,非地表一致性剩余静校正,Trim3D 叠前道集时差微调,3D非地表一致性剩余静校正是利用经过动校正的CMP道集内的各道与输入的模型
32、道，在用户定义的时窗内做互相关处理，获得相对时移，即非地表一致性静校正量，并对原始CMP道集在频率域通过相位变换实现时移。模块还可以通过叠代平均技术得到炮点和检波点的地表一致性剩余静校正量。该模块既可用于二维处理，也可用于三维处理。,非地表一致性静校正要点,只能在静校正基本解决的情况下使用,采用较大时窗和较小时差的非地表一致性静校正,模块输入动校正后的道集数据，输入控制为line、CMP、trace，旗标cmp（二维为CMP、trace，旗标CMP），使用模型道文件stk-srd-rs-rna-pure-all计算非地表一致性剩余静校正并应用道数据体上（非地表一致性剩余静校正无法输出校正量文件
33、），输出地震数据供后续处理使用。同该模块的时窗参数的定义与地表一致性静校正计算一样。该模块处理输出的地震数据可直接用于后续处理，而在计算时由于动校正已经完成了数据体切除且不可恢复，所以在进行该处理时，切除库的定义需要谨慎。,叠加数据要和输入的道集数据基准面一致,与该模块匹配的输入模块关键字，3D：LINE-CMP；2D：CMP-TRACE,叠前道集时差微调前/后叠加剖面,Trim3D,前,后,剩余静校正和速度分析的迭代过程：,野外静校正,速度分析,折射波静校正,剩余静校正,CMP参考面校正,统一面校正,迭代,1、优势频带里求取剩余静校正量,2、剩余静校正与速度分析迭代,3、剩余静校正时窗选取,
34、有连续的强反射层追踪,构造有起伏应空变，角度要算准,无标准层，需开更大时窗注意边界处时窗的选取,小结,4、质量控制 检查剩余静校正数据有无异常 剩余静校正前后大道集对比 剩余静校正前后叠加剖面对比 剩余静校正前后速度谱的对比,5、注意 剩余静校正前要加AMPEQU L500 强化模型道 填足覆盖次数、跨炮道数、跨CMP数 标准层平时可适当加大倾角扫描道数,小结,静校正是地震资料处理的基础，静校正的好坏直接关系到整个地震资料处理的成败，尤其是对复杂地区的地震资料处理而言更是如此。因此，在实际地震资料处理过程中，如何选择合适的处理思路和相应静校正技术，把静校正工作做好是我们一直探索的目标。GeoEast系统为我们提供了多种室内静校正处理的方法，并形成了一套比较完整的静校正配套技术，就目前而言，能够较好的解决静校正、尤其是复杂地区的静校正问题。,Thank You!,
