地球物理学进展  2015, Vol. 30 Issue (5): 2349-2353   PDF    
动校剩余时差处理方法及应用
周鹏1, 刘志斌1, 张益明1, 马德蕻2    
1. 中海油研究总院, 北京 100028;
2. 河南地质博物馆, 郑州 450016
摘要: 常规采集和处理后的地震道集经常动校不平,存在剩余时差,对叠加成像、叠前反演和AVO属性等分析存在影响.为拉平地震道集,提出了一种与剩余时差无关的道集拉平方法-绝对值互相关方法.该方法将地震道集的每一道与参考道在选定的时窗内进行绝对值互相关计算,利用相关系数门槛值的控制剔除异常相关系数值,得到剩余校正时移量场,并对其进行平滑处理,最终利用该时移量场校平地震道集.绝对值互相关具有良好的抗噪性,能适应低信噪比道集及道集极性反转等各种异常情况.多个地区的应用表明,该方法能较好地校平叠前地震道集,改善道集的品质,并且具有较好的适应性.
关键词: 道集拉平     剩余时差     互相关     叠加成像     叠前反演     AVO属性    
The processing method and application of the residual moveout NMO
ZHOU Peng1, LIU Zhi-bin1, ZHANG Yi-ming1, MA De-hong2    
1. CNOOC Research Institute, Beijing 100028, China;
2. Henan Geological Museum, Zhengzhou 450016, China
Abstract: The seismic gathers can not be flatten after conventional processing, and there are often residual movemout. The impact of these exist on seismic imaging, pre-stack inversion and AVO attribute analysis. this paper presents a kind of method-the absolute value cross-correlation method, that has nothing to do with the remaining difference to flatten the seismic gathers. This method is that every seismic gather perform absolute value of the cross-correlation calculation with reference seismic trace within the selected window. The abnormal correlation coefficient can be excluded by using the correlation coefficient threshold.Then we can abtain the remaining amount of correction shift field and smooth it. Finally, the seismic gathers are flatten by the shift field. the absolute value cross-correlation method has good noise immunity. It can adapt to the low SNR gathers and gathers polarity reversal and other anomalies.The method has been applied in a number of areas. The good effect can be got in seismic gathers flatting and removing the remaining difference. It can improve the quality of seismic gathers and section imaging.
Key words: flatten gathers     residual movemout     cross correlation stack imaging     pre-stack inversion     AVO attribute    
 0 引 言

随着勘探逐渐从叠后走向叠前,地震道集的应用越来越广泛.但在实际的工作中,经过了常规处理后的地震道集仍然有动校正不平的现象,其中一个重要原因就是地震道集存在剩余时差.在常规采集和地震资料处理中,都有一个共中心点道集的反射同向轴的时距曲线为双曲线假设前提,通过常规速度分析与动校正,地震道集很容易被拉平.实际介质基本不满足这种假设,而且随着勘探的深入,复杂构造越来越多,地下情况也越来越复杂.地震数据受采集条件、地层介质各向异性、噪声干扰等影响,地震同向轴偏离了正常的时距曲线形态.在这样的条件下,利用常规速度分析取得的速度动校正之后,有剩余时差,道集品质偏低.

去除剩余时差,进行地震道集拉平的方法很多,如精细速度分析方法和技术(Swan,2001薛冈等,2003郭树祥等,2004程玉坤等,2008崔宝文和王维红,2010),通过各种速度分析手段或者在更密的网格上进行细致的速度分析,得到动校正速度,应用于道集,消除剩余时差;地震道集互相关技术(Hinkley et al.,2004Gulunay et al.,2007慎国强等,2010张征等,2012),通过实际道集与模型地震道的互相关来计算时移量并进行道集校正;相位匹配或替换技术(Lichman,1999林伯香和孙建国,2001王开燕等,2007),是根据相位谱的理论,首先计算模型道或者参考道的相位谱,然后对实际道集进行相位谱替换,达到校正的目的;刘洋等人基于线性连续速度模型的速度分析与动校正方法(刘洋和魏修成,2003),通过利用非双曲方法处理,在线性连续速度模型基础上,研究速度随深度的变化率来进行速度分析和动校正的方法;李钟晓等基于地震图像校正的共成像点道集增强技术(李钟晓等,2013),可以实现了不同偏移距道集在时间和空间上的逐点匹配对齐;还有其他的一些处理方法(陆文凯等,2001周青春等,2009)、道集强制拉平方法等.上述的研究方法或技术对道集拉平各有效果,但通常都有一定的应用条件,精细的速度分析方法的前提是反射同向轴满足正常的旅行时曲线方程,相位道集拉平方法及相干法要求地震资料信噪较高,邻近同相轴的剩余时差较小,而且地震同相轴没有极性反转问题.但实际地震资料情况复杂,常规速度分析手段仍然校不平的同向轴很多都不符合地震波反射双曲线规律等各种异常情况(Lichman,1999王云专等,2006刘洪林和朱秋影,2007张津海等,2012).

地震道集动校不平的这种现象普遍存在,影响了叠加成像,及后续的叠前反演、AVO属性分析等(Singleton,2009),需要进行必要的地震道集拉平处理.文中的处理方法——绝对值互相关方法,是在常规速度分析动校正的基础上,选取合适的时窗,利用地震道集的每一道与设定的参考道在选定的时窗内进行绝对值互相关计算时移量,并且根据实际资料情况,选取合适的绝对值相关系数门槛值,利用时窗、时移量、平滑参数的控制,在保真的前提下,实现各种道集拉平.该方法在多个地区的数据处理中进行了应用,取得较好效果.

1 方法原理 1.1 方法及处理流程

一个地震共反射点道集同向轴代表的是来自同一反射点不同偏移距的地震反射波,沿偏移距方向,同一反射点的地震反射波相邻道之间具有最大的相似性,因此同一段时窗内的各界面反射在不同偏移距的道之间也具有最大相似系数,绝对值互相关方法正是基于这点拉平道集.该方法通过三步拉平地震道集——首先利用绝对值互相关计算时移量,其次对时移量进行平滑校正,最后用时移量将原始道集进行拉平处理.在应用过程中根据不同地区、不同资料、不同的信噪比情况选择不同的参数加以控制,方法如下:

(1)设地震道集的每一道函数为f(t)(图 1b),选择一个近道或者利用小偏移距的一组道进行叠加形成参考道R(t)(图 1a);

(2)选择一个时窗,并将该时窗内的道集与参考道进行绝对值互相关,求得每道与参考道在该时窗内的最大绝对值互相关系数所对应的时移量,公式为

Coef为绝对值互相关系数,t0t1为时窗的起点和终点,该时窗内对全偏移距范围内都至少包含了一个波组;使用绝对值计算互相关系数,保证了相位反转同向轴与低信噪比道集的相位信息前后一致;

(3)将时窗沿时间方向滑动至道尾完成一个道集的处理工作,依次对每道处理,得到整个工区的校正时移量;并且根据地震道集的实际情况,设立最大时移量和最小绝对值互相关系数门槛值,剔除掉异常时移量,并进行适当的平滑处理;

(4)应用校正时移量场对整个地震道集进行拉平处理,去除剩余时差(图 1c).

图 1 互相关道集拉平处理示意图 Fig. 1 The cross-correlation gathers flattening processing

绝对值互相关道集拉平方法的技术流程如图 2,为了取得较好的效果,流程中通常先应用大时窗进行初步拉平,再应用小时窗进行精细处理.

图 2 绝对值互相关道集拉平处理流程 Fig. 2 Gather flattening processing flow of cross-correlation of absolute value
1.2 实际应用中的问题及解决方法

实际地震道集由于地质、采集条件及处理方法的不同,道集的信噪比、分辨率、同向轴的连续性都不相同,因此在绝对值互相关方法应用的过程中会出现相关系数的各种异常情况,应用中应注意如下几个方面:

(1)参考道的选择:参考道对结果影响较大,一般选择随机噪声小的近道或者近偏移距的几道叠加形成参考道;在选择一定偏移距范围内叠加道作为参考道时,选近道时考虑多次波残留较少的情况,选全偏移距叠加作为参考道时,不适合剩余时差较大或同相轴反转的资料;另外可以选择经长时窗初步校平的数据的全偏移距叠加道作为短时窗校平时的参考道.

(2)噪声的问题:地震数据的噪声干扰可能导致相关系数过大或者过小,出现不合适的时移量,导致地震数据出现不合理的压缩或拉伸特征.

根据数据信噪比,设置合适的相关系数门槛值,剔除低相关系数所对应时移量;设置可随偏移距变化的最大容忍时移量,从近道至远道,容忍时移量逐渐增大;多道之间的时移量通过横向一致性检查,平滑时移量的空间变化,剔除异常值.

(3)同向轴极性反转的问题:绝对值互相关的计算本身已经解决了相位反转的问题,可以校平相位反转的同向轴,但极性反转位置邻近的几道信噪比低,时移量求取难.

在极性反转处,依据对噪声问题的处理方法,设置相关系数门槛值,利用求取的临近道时移量内插得到校正时移量.

(4)校正时窗长度的选择:大时窗校正能力有限,不容易完全校正同向轴;小时窗校正效果显著,但对信噪比低的数据,会出现误校正.

采用长短时窗串联的方式,先通过大时窗初步校平整个地震道集,再应用小时窗针对目标进行精细处理;通过比较各时窗所计算出的相关系数,选出合理的时窗,并且每个时窗长度对所有的偏移距范围内的道集都至少包含目的层的一个波组.

2 数据处理 2.1 模型数据试验

为了验证方法对于复杂条件下的正确性,进行了复杂模型的模拟及处理.从图 3b的水平横线可以看到模型是一个没有校平的道集,并且该道集的相位发生了极性反转.模型随机噪声比较强,同向轴的极性反转处振幅与噪声强度基本相当.由于存在了极性反转,而且同向轴相位反转处的噪声强度高,这时采用常规相干方法或者采用相位相关方法已不适用,不能拉平道集.采用绝对值互相关方法进行校正处理,参考道的选择上,把模型前10道平均当作参考道(图 3a),互相关系数门槛值为70%以上,窗长为60 ms,包含了一个波组,最大时移量门槛值为8 ms,结果如图 3c.从图 3c可以看到,水平横线处同向轴被拉平,极性反转处没有出现畸变,道集的振幅得到了很好的保持.

图 3 道集拉平模型数据检验
(a)参考道;(b)需拉平道集;(c)拉平后道集.
Fig. 3 Data validation of gather flattening model (a)Reference trace;(b)Original gather;(c)Flatten result.
2.2 实际资料处理

作为实例的地震资料是一条二维测线,该地震资料道集采样记录为2 ms,经过了预处理、振幅补偿等常规处理,以及精细的速度分析和动校正,但在道集中仍然存在着剩余时差,同向轴没有完全被拉平.从图 4a中可以看到地震记录同向轴已经不符合动校正的双曲线时差曲线,图 4a红线标示的同向轴为了拉平近道同向轴,致使远道不平,并且常规速度分析的精度限制,致使许多同向轴动校过量(如图 4a黑线标示处),以及有些同向轴本身由于一些原因没有规律性(图 4a绿线标示处).显然常用的速度分析、动校正已经不能将该道集拉平.采用绝对值互相关拉平方法进行处理,利用前6道进行加权平均生成参考道,绝对值互相关系数门槛值为82%(大于82%的进行拉平),最大时移量为12 ms,第一次校正时窗为200 ms,进行初步校正,第二次校正时窗为70 ms(包含了一个波组),处理后的结果如图 4b.对比处理前后道集,可以看到处理后道集规律性明显增强,地震同向轴基本被拉平,连续性光滑性和分辨率得到了改善,减少了同向轴的错位现象,明显提高了道集的质量和叠加精度.并且沿着图 4中道集的同一个同向轴提取拉平前后的AVO曲线,可以看到道集拉平处理前(红色)后(蓝色)AVO曲线没有变化,地震振幅也得到了很好的保持.

图 4 道集拉平处理 Fig. 4 Gather flattening processing

将拉平前后的地震道集进行了叠加成像处理.拉平处理前的叠加剖面如图 5a,拉平处理后的叠加剖面如图 5b.从拉平前后的对比图来看,在已经拉平的道集位置处理前后叠加剖面变化不大(蓝色椭圆区域外).道集动校不平或者动校过量的区域(蓝色椭圆区域内),拉平前地震剖面振幅较弱,连续性不是很强;拉平后,剖面能量得到了加强,成像质量改善,构造变得清晰,同向轴连续性得到了提高.整个剖面的信噪比分辨率都有较大改善,给构造解释和储层预测打下了一个较好基础.

图 5 道集拉平处理后叠加剖面 Fig. 5 Stack section after gather flattening processing
3 结 论

绝对值互相关道集拉平方法,是在常规资料处理的基础上对道集进行一进步优化的方法.它可以根据实际资料情况,选取合适的参考道,将道集与参考道在选定的时窗内进行绝对制互相关计算时移量,并利用相关系数门槛值、时移量、平滑参数的控制,在保真保幅的前提下,实现道集拉平.

绝对值互相关方法能使地震道集实现同相叠加,为叠前反演和AVO分析提供合格资料,并且改善叠加效果,提高地震剖面信噪比、分辨率和连续性.通过模型和应用证明绝对值互相关方法是简单、有效、实用的地震道集拉平方法.

致 谢 感谢审稿专家提出的宝贵修改意见和编辑部的大力支持!

参考文献
[1] Cheng Y K, Ran J B, Xiao W. 2008. CRP gather based prestack data processing and its application[J]. Progress in Exploration Geophysics (in Chinese), 31(1): 38-43.
[2] Cui B W, Wang W H. 2010. Application of prestack time migration to seismic imaging of deep fault depression areas[J]. Progress in Geophysics (in Chinese), 25(5): 1703-1707, doi: 10.3969/j.issn.1004-2903.2010.05.025.
[3] Gulunay N, Magesan M, Roende H H. 2007. Gather flattening[J]. The Leading Edge, 26(12): 1538-1543.
[4] Guo S X, Han Y Z, Li J M, et al. 2004. The optimal method of velocity analysis in high resolution seismic data processing [J]. Geophysical Prospecting for Petroleum (in Chinese), 43(1): 80-82.
[5] Hinkley D, Bear G W, Dawson C. 2004. Prestack gather flattening for AVO[C].// 74th Annual International Meeting, SEG, Expanded Abstracts, 271-273.
[6] Li Z X, Lu W K, Wang J. 2013. Common image gathers enhancement based on seismic image warping[J]. Chinese Journal of Geophysics (in Chinese), 56(5): 1592-1608, doi: 10.6038/cjg20130517.
[7] Lichman E. 1999. Automated phase-based moveout correction[C].// 69th Annual International Meeting, SEG, Expanded Abstracts, 1150-1153.
[8] Lin B X, Sun J G. 2001. Residual moveout correction by using phase replacement[J]. Geophysical Prospecting for Petroleum (in Chinese), 40(3): 15-22.
[9] Liu H L, Zhu Q Y. 2007. AVO inversion and interpretation based on prestack depth migration[J]. Progress in Geophysics (in Chinese), 22(3): 905-912.
[10] Liu Y, Wei X C. 2003. Velocity analysis and normal moveout correction based on linear continuous velocity model[J]. Petroleum Exploration and Development (in Chinese), 30(5): 71-74.
[11] Lu W K, Zhang X G, Li Y D, et al. 2001. High resolution stack velocity determination based on signal estimation[J]. Chinese Journal of Geophysics (in Chinese), 44(6): 871-876.
[12] Shen G Q, Wang Y M, Meng X J, et al. 2010. Application of seismic gather flattening technique based on time-frequency analysis[J]. Journal of China University of Petroleum (Edition of Natural Science) (in Chinese), 34(1): 34-36.
[13] Singleton S. 2009. The effects of seismic data conditioning on prestack simultaneous impedance inversion [J]. The Leading Edge, 28(3): 260-267.
[14] Swan H W. 2001. Velocities from amplitude variations with offset[J]. Geophysics, 66(6): 1735-1743.
[15] Wang K Y, Li H, Guan B, et al. 2007. Application of residual moveout correction to the improvement of resolution [J]. Journal of Daqing Petroleum Institute (in Chinese), 31(6): 14-17.
[16] Wang Y Z, Yang L W, Li S H. 2006. Residual moveout correction and Poisson’s ratio inversion[J]. Progress in Geophysics (in Chinese), 21(1): 214-218.
[17] Xue G, Wang L S, Hu Z P. 2003. The comparison of NMO Method for long offset seismic data[J]. Oil Geophysical Prospecting (in Chinese), 38(2): 151-155.
[18] Zhang J H, Zhang Y Y, Sun Z D. 2012. The effects of seismic data conditioning on pre-stack AVO/AVA simultaneous inversion[J]. Oil Geophysical Prospecting (in Chinese), 47(1): 68-73.
[19] Zhang Z, Di B R, Hu Y, et al. 2012. Application of pre-stack detailed process technique on CRP traces gather[J]. Science & Technology Review (in Chinese), 30(27): 35-40.
[20] Zhou Q C, Liu H S, Kondrashkov V V, et al. 2009. Research of dual parameters evolving common reflection point stack and velocity analysis method[J]. Chinese Journal of Geophysics (in Chinese), 52(7): 1881-1890, doi: 10.3969/j.issn.0001-5733.2009.07.023.
[21] 程玉坤, 冉建斌, 肖伟. 2008. 基于CRP道集的叠前处理技术及应用[J]. 勘探地球物理进展, 31(1): 38-43.
[22] 崔宝文, 王维红. 2010. 叠前时间偏移技术在深层断陷区地震成像中的应用[J]. 地球物理学进展, 25(5): 1703-1707, doi: 10.3969/j.issn.1004-2903.2010.05.025.
[23] 郭树祥, 韩永治, 李建明,等. 2004. 高分辨率地震资料处理中的优化速度分析方法[J]. 石油物探, 43(1): 80-82.
[24] 李钟晓, 陆文凯, 王季. 2013. 基于地震图像校准的共成像点道集增强技术[J]. 地球物理学报, 56(5): 1592-1608, doi: 10.6038/cjg20130517.
[25] 林伯香, 孙建国. 2001. 相位替换法剩余时差校正[J]. 石油物探, 40(3): 15-22.
[26] 刘洪林, 朱秋影. 2007. 基于叠前深度偏移的AVO反演及解释[J]. 地球物理学进展, 22(3): 905-912.
[27] 刘洋, 魏修成. 2003. 基于线性连续速度模型的速度分析与动校正方法[J]. 石油勘探与开发, 30(5): 71-74.
[28] 陆文凯, 张学工, 李衍达,等. 2001. 基于信号估计的高分辨率叠加速度分析[J]. 地球物理学报, 44(6): 871-876.
[29] 慎国强, 王玉梅, 孟宪军,等. 2010. 基于时频分析的地震道集校平技术应用[J]. 中国石油大学学报(自然科学版), 34(1): 34-36.
[30] 王开燕, 李慧, 官波,等. 2007. 相位剩余时差校正方法在提高分辨率中的应用[J]. 大庆石油学院学报, 31(6): 14-17.
[31] 王云专, 杨立伟, 李素华. 2006. 剩余时差校正及泊松比反演[J]. 地球物理学进展, 21(1): 214-218.
[32] 薛冈, 王良书, 胡中平. 2003. 大炮检距地震资料动校正方法比较[J]. 石油地球物理勘探, 38(2): 151-155.
[33] 张津海, 张远银, 孙赞东. 2012. 道集品质对叠前AVO/AVA同时反演的影响[J]. 石油地球物理勘探, 47(1): 68-73.
[34] 张征, 狄帮让, 胡因,等. 2012. CRP道集叠前精细处理技术应用[J]. 科技导报, 30(27): 35-40.
[35] 周青春, 刘怀山, Kondrashkov V V, 等. 2009. 双参数展开CRP叠加和速度分析方法研究[J]. 地球物理学报, 52(7): 1881-1890, doi: 10.3969/j.issn.0001-5733.2009.07.023.