《omega剩余静校正相关分析拾取反射时间模块解析》由会员分享,可在线阅读,更多相关《omega剩余静校正相关分析拾取反射时间模块解析(29页珍藏版)》请在金锄头文库上搜索。
1、.Omega地震资料处理系统XPERT和REFLECTION_MISER模块剖析摘要:对XPERT和REFLECTION_MISER模块进行剖析,结合理论知识和实际工作中所积累的一些经验对如何正确完成其中的每一个步骤进行详细的指导。关键词:CMP,XCMP,in-line,crossline,REFLECTION_MISER,XPERT,模型,模型片段,模型扩展,模型segment,时移量限定值(shiftlimit),相对模型加权,修改过模型片段加权,时窗标记。前言静校正问题在地震资料数据处理中是一个十分重要的问题,对于地形起伏剧烈的四川和有类似情况的其它工区来说尤其的关键。静校问题解决的好
2、坏直接影响到处理出的地震资料成像效果和是否能真实反映地下构造的形态。为了能在生产中更好的解决静校正问题,在此特地对Omega地震资料处理系统中静校正处理模块XPERT和REFLECTION_MISER,特别是XPERT进行了深入的剖析。一、模块概述相关分析拾取反射时间模块(XEPERT)为MISER这些剩余静值分析模块求取反射时间和质量因子,其也可用于非地表一致性或平滑静校正(trim)。XPERT两个输入项一为按CMP点排序并经动校后的地震道文件,另一为相同道的叠加文件或外部模型文件。XEPERT的用途是为后面所用的模块(通常是MISER)提供拾取时间,该模块确定地表一致性的剩余反射静值。剩
3、余静校正量的应用对所有陆地和一些浅海资料的处理是必不可少的部分。二维和三维数据在XEPERT中用相同方法处理,二维被看作是三维的特例。除一些特别说明的例子外,文章里所说的都适用于二维和三维资料。文章中的CMP指的是二维CMP和三维面元CMP。数据模型由用户控制、倾角调向平均值窗口、围绕被拾取的CMP道集的叠加数据组成,数据对从以下两个叠加道中选出的模型道有效:第一个由已被XEPERT处理过的动校后的叠加道构成;另一个数据为输入的叠加数据。对信噪比够的数据,可对含有相关分析的倾角在自动倾角调向上使用自动追踪来制造模型和追踪表面构造。用户提供的倾角区域允许很大的增加在低覆盖或变覆盖地区、低信噪比数
4、据区、三维工区边界和二维测线末端的拾取成果,在REFLECTIONMISER中也能增加预测构造的稳定性以用于控制不耦合。用户提供的参数,如时移量限定值(SHIFT-LIMITS)、模型制造窗口、空间扩展和道加权都是空变的,其他的可选的参数在三维处理中对三维资料使用。XEPERT输出的是拾取静校正量文件、模型片段文件,例如向上提交的高级叠加道和最终模型文件。二、输入数据和数据准备XEPERT在内部或外部模型模式中运行,模式由被地球物理语言确定了的输入模型(在INIT_MODEL_FRAC或EXTERNAL_MODEL中)由用户选择决定。对内部模型模式,输入NMO和模型片段文件;对外部模型模式,输
5、入NMO和外部模型数据。1、NMO数据准备用于拾取的输入数据包含有常规处理流程中的误差值,处理员可以作一些工作来改善拾取处理和静校正。对NMO数据的一些必须的常规指导如下:1. 数据必须经过预处理和按CMP点排序。2. 数据应该作过合适的球面扩散补偿和衰减以及反褶积。注意在噪音大的数据中一道接一道(tracebytrace)的反褶积能造成明显的道对道(tracetotrace)的时移,这种并不是单纯的静校正问题,因此会影响时间拾取。对这类数据,在剩余静校正前作地表一致性反褶积是合适的反褶积方法,虽然其在压制震源产生的噪声方面的效果没有一道接一道(tracebytrace)的反褶积好,并且也要影
6、响时间拾取。滤波、增益或外科切除处理可用来解决这个问题,而这种方法是否有助于拾取没有经过反褶积的数据,现在仍没有定论。3.数据应该进行过动校正。NMO速度应该是有实质地质意义,并且是在速度分析中从很少或没有静校正量影响的区域中选取的最佳速度。在第一次MISER运行时,应用的速度在横向上应该是缓慢变化的(在整个构造的趋势上对测线或整个工区仅有大范围的一致性),如果数据很平缓或测线不长、工区较小时,一个速度函数对第一次运行MISER就足够了(无论甲方要求进行多少次分析都没问题)。另外应该对数据进行切除,就像数据必须被叠加一样。MISER和NMO速度应用的关系很重要,MISER试图找到与先前所应用的
7、NMO速度场相符合的静校正问题的解决方法。换种说法,即MISER试图用一个像你开始所用的速度场的最终速度场来完成,因为velcity_analysis_implied速度中速度与速度(velantovelan)间的误差要影响静校正,应用这种速度来进行拾取可能会拾取和MISER中数据真实的静校条件,所以将造成叠加质量变差。另外,我们在工作发现,对于能量较差、静校正问题严重的资料,作第一次速度分析时你会发现速度谱中能量团不收敛,进行准确的速度分析是很困难的。对这种资料,处理员可以不作第一次速度分析,直接用中心速度来作动校正,对初叠数据作了第一次REFLECTION_MISER后再进行第一次速度分析
8、。这时你会发现速度谱中的能量团收敛性变好,将有利于你进行准确的速度分析。如果在一个地质构造上较平缓、有严重的静校问题的区域,你进行了很多次速度分析,不加选择的拾取波峰,并应用易变化的速度场对数据进行拾取,这样就会因使用错误的速度而造成对解决静校正问题无法弥补的损害。如果你在平滑动校速度场时采用了two-cable-lengths滤波器就会降低出现以上损害的风险。三或更多cable-lengths的滤波器的效果可能会更好,但过于相信这种平滑滤波器的作用可能会造成速度的不稳定。对类似情况,最好在测线品质较好的一段选择一个速度函数并应用于整条测线,如果测线显示有区域性的倾角,则在测线的两端选择两个速
9、度函数并进行平滑的内插。应该使用区域NMO速度直到解决了静校正问题为止,因为应用的速度函数如变化太快的话要影响MISER的解决。在应用速度函数时,MISER对速度平滑变化的情况可以不考虑,能得到相对较好的结果,但在速度变化迅速的情况下,则不要期望MISER对静校正问题的解决有什么帮助。4.滤波因改善了相关分析时窗里数据的信噪比从而改善了道对模型的相关分析,使拾取效果更好并且改善了MISER得到的结果。各种各样的滤波和增益都可用来提高数据的信噪比,例如对严重的地滚波就可进行外科切除、低截滤波或F-K滤波。滤波的目标就是对影响了拾取结果的噪声进行衰减。通常,如果要作一个简单的频率滤波可以选择以下滤
10、波器:时间/空间/变偏移距信号道滤波器、multi-channel或其他(如F-K)滤波器等,在数据中有噪声的区域划出滤波操作时窗对去噪也有很好的效果。减弱数据中信号控制的部分与频率或时空中的任一个因素都将影响MISER的解决。multi-channel操作,如F-K滤波,将会影响由静校问题所引起的道对道(tracetotrace)的数据的时移,因此可能降低静校正的质量。不经过试验很难决定用哪一种滤波器才能有好的效果,所以只有当你确定选用的滤波器能改善MISER的质量时才对数据进行滤波。如果你选用了一个简单的频率滤波器,它对模型道进行滤波是很省机时的而对NMO数据则不然(一个滤波后的道和没滤波
11、的模型将被当作没滤波的道和滤波后的模型产生一个相同的相关分析函数)。既然模型道的滤波不能在XPERT中实现,那么就必须在XPERT外进行滤波。对NMO数据的预拾取进行滤波的相对于对数据进行滤波两倍的效果。从滤波后的动校数据中选出的模型道同没滤波的数据中选出的滤波后模型是一样的。当这个模型同滤波后的NMO数据进行相关时,相关分析函数获得了应用两次滤波的效果和两倍的陡度。为了在外部模型模式下避免这个问题,用户应对外部模型道而不是NMO数据进行滤波。在内部模型模式中,用陡度只有期望效果滤波器一半的频率滤波器对NMO数据预滤波。对滤波后的NMO数据进行叠加或对没滤波的NMO数据叠加后再对其用半陡度滤波
12、器进行滤波,就能生成模型片段文件。最后将滤波后的文件提交给XPERT,这样作的效果就跟内部生成的模型道被所期望的全陡度的滤波器进行滤波一样。5.增益为了均衡相关分析窗口里的振幅、对球面扩散与衰减进行补偿,时间/偏移距/地表一致性振幅调节、大时窗(线形变化)增益和(或)道均衡等可以应用。这些处理使那些被采样来用于相关分析和制造模型的每一道数据看起来都很均衡。小时窗增益(包括使数据的振幅能够很快的变得相同的自动增益)可用来去噪,当在相关时窗内噪声的振幅比信号的振幅大时,增益的作用同滤波是一样的,即相对于信号来说减弱了噪声的振幅从而改善了相关。但是要注意的是,在信噪比高的数据上运用小时窗增益可能因提
13、高数据中小振幅的振幅而使有效波同相轴信噪比降低、减弱大振幅而使同相轴信噪比提高,这样可能使相关分析变差,因此在类似的情况下最好不使用小时窗增益。通常情况下,对那些持续时间较短的噪声,较好的解决方法是外科切除、区域切除或减少噪声数据区域。6.基准面NMO数据应参考和校正到一个水平基准面。对陆地资料,应该在动校之前将数据校正到CMP基准面,因此在HyperbolicandLinealMoveout模块中道头字应选用STATCOR_DETECT和STATCOR_SOURCE。动校后在基准面校正模块(DatumCorrection)中选用道头字CMP_DATCOR将数据校正到水平基准面。包含了过渡区域
14、数据和应用了折射静值的海洋资料同陆地资料一样,采用同样的方法校正到水平基准面。对只有剩余静值的海洋资料,数据只影响到震源和检波点高程间一半位置的一个平面上,对该类数据不计算CMP基准面,只应用前面所作的MISER成功运行后所产生的剩余静值。2、初始模型片段数据对一个叠加数据进行拾取。对于内部模型模式来说,需要模型片段数据的目的是靠输入叠加数据来进行两次覆盖,在道集被相关来用于生成构成模型的模型片段文件之前,XPERT不用运行,在XPERT前生成的数据(尤其是三维数据)是很庞大的。对于外部模型模式,输入各种叠加数据来构成模型都是被允许的。输入没有进行增益的NMO数据的叠加数据的效果较好,另外一些
15、叠加数据也是可以选择但不是必须的,如作了随机噪声衰减(RNA)处理的叠加数据。无论输入哪一种叠加数据,它都应该同NMO数据一样校正到同一个水平基准面上。对每一个在NMO数据中被拾取的CMP点应该有一个活跃的叠加道存在,如果没有活跃的叠加道,则CMP点的模型道不能制造,CMP点不能被拾取。3、外部模型数据用户由外部模型输入直接向XPERT提交模型。应该使用确定始末时间的加添窗口(INRTERPOL_WNDWS)来画出外部模型时窗。在这个版本的XPERT中,当选用外部模型模式时,不要期望输出模型片段数据。另外,外部模型同样也要被校正到同动校数据一样的水平基准面上。三、输出数据共输出了三套数据:静值
16、拾取数据、最终模型数据和模型片段数据。1、静值拾取数据静值拾取数据是XPERT的主要输出并作为MISER的输入。在确认输入的NMO数据中的道头字TRACE_BALANCE_FRCTOR和CRL_STACK_WEIGHT中包含了拾取静值、道对模型的相关分析确定了相关系数后,这些道头字被一对一的复制到静值拾取数据中。相关函数和(或)属性被输出,这个输出项是可选的,其被相关拾取参数设置所控制,以允许的格式输出到道数据缓冲区上。如果选择输出相关分析函数,相关函数里的时间代表着动校道相对于其模型的反射时间。例如,如果相关函数中拾取的波峰显示的时间是5ms,则指NMO道比其模型反射时间大5ms。相关函数的零lag值在输出上产生0ms反射时间,相关函数将下一个100ms的MULTIPLE的两侧的零lag值