随着地震勘探的不断深入,复杂近地表的地震勘探活动越来越多.在野外采集中由于受地表条件的限制,经常会需要使用不同的震源激发和检波器接收.地表激发接收条件的不一致会导致同一区块地震记录的子波存在一定的差异,表现在地震数据上就是同相轴的不一致性 (云美厚等,2006).
随着多块三维连片处理、时延地震处理、过渡带探区资料处理的普及,不同年代施工、不同采集方法、不同地表条件的小工区地震资料往往需要集中起来处理.由于采集因素不同,各块三维数据之间往往存在明显差异.提高地震资料的一致性成为这类资料处理的关键环节,甚至是决定性环节 (庄东海等,1999;云美厚等,2006;吴琼等,2008;裴江云等,2011;郭念民和吴国忱,2012).
随着地震属性反演由叠后发展到叠前,地震叠前弹性阻抗反演要求更大的入射角范围,即需要更大的偏移距范围,实现动校拉伸效应更小,这就要求在保幅条件下消除道集上的各种时差,校平同相轴,以提取正确的泊松比、AVO信息等 (王云专等,2006).
随着地震成像精度的要求越来越高,要求共中心点 (CMP) 道集或叠前偏移得到的共反射点 (CRP) 道集更加平直.另外,多方位角地震资料需要进行一致性匹配处理,以便偏移成像结果具有更好的断层阴影压制能力和更高的断层识别能力 (朱江梅等,2013).
可见,地震数据处理中波组的一致性处理是提高处理成果品质的关键.一致性处理的好坏直接关系到处理成果的信噪比、分辨率和保真度.如果不能解决好数据的一致性问题,就无法进行精细构造解释和储层反演.出于地震资料解释工作的需要,开展地震资料的一致性处理研究就显得非常重要,其现实意义和实用价值不言而喻.
1 问题的提出由于近地表的不规则性以及低降速带影响,CMP道集中反射波经常无法满足双曲线规律.在常规流程中使用动校正及地表一致性剩余静校正方法后的道集内还存在着随时间纵向变化的剩余时差.由于地震波经过地下介质的透射、反射和吸收作用,使反射波发生时间延迟、相位畸变、振幅衰减.这些因素对各个地震道是不一样的,且往往是时变的.由于实际处理中各向异性叠前深度偏移的参数难以获得准确的值,会影响到远炮检距排列的校平,难免会存在部分排列不能完全拉平的问题 (赵玲芝等,2011).种种不一致性问题,其结果造成在共深度点道集或共反射点道集叠加时,由于各道的不一致性而影响叠加剖面的品质,使信噪比和分辨率不能得到有效提高,对叠前属性提取也会造成极为不利的影响.水平叠加技术降低分辨率是影响其应用的一个主要障碍,其根本原因在于水平叠加的两个基本前提条件 (波形一致、时间对齐) 在实践中没有得到很好的满足 (夏洪瑞等,2003).
针对不一致性问题,人们研究了许多有效的技术并在生产中得到了良好的应用效果.目前,国内外地震数据一致性处理除了采用常规的地表一致性振幅补偿、地表一致性反褶积、地表一致性剩余静校正和动校正畸变切除等处理技术 (安西峰等,2010) 外,主要处理技术还有:基于相位替换为基础的时差校正技术 (林伯香和孙建国,2001),其实质是信号的相位谱中包含了信号的到达时,通过改变相位谱来改变信号到达时;基于分形理论的初至拾取算法 (刘伊克等,2001),其实质是针对观测走时存在震源子波波形的不一致性、折射波相位反转等原因引发的初至走时拾取中的相位跳跃误差,提高观测走时数据的精度;基于数学物理模型的反射振幅一致性校正技术 (李生杰等,2002),其实质是依据模型及其实现条件使用共轭梯度方法在时间域对地震记录进行谱分解处理,校正因复杂地表条件等因素对地震反射振幅而产生的影响;基于关系公式的剩余动校正量修正技术 (苑春方等,2002),其实质是通过研究带有剩余动校正量的CMP道集内反射波接收时间与偏移距之间的关系,得到了剩余时差函数;时频域有色谱校正技术 (夏洪瑞等,2002),其实质是利用模型道法、互相关法及已知波形法给出期望的有色谱,为校正提供依据,利用小波变换的正交性来克服傅氏变换引起的频泄现象;分偏移距动校正技术 (王鑫和陈新荣,2003),其实质是根据不同的偏移距对速度响应的差异,计算并消除这种差异所造成的CMP道集内的剩余时差,使同相轴的连续性和平滑性增强;匹配滤波法处理技术 (邬达理,2004),其实质是利用近炮检距数据为模型,通过匹配滤波实现远近排列数据的一致,克服大炮检距波组的低频化现象;二步法反褶积处理技术 (王西文等,2006),其实质是先采用地表一致性反褶积, 将不同震源的频带拓宽到一个标准上,再采用相位校正反褶积, 将不同震源的数据校正到相同相位上;利用地震属性求取剩余静校正技术 (魏庚雨等,2009),其实质是首先对参加相关的道集和叠加数据进行地震属性处理,然后进行剩余时差的拾取;改进型时间域剩余动校正技术 (马永军和王季,2010),其实质是通过提取同一地震道中不同时刻采样点的剩余时差来获取时移量;基于结构自适应中值滤波器的随机噪声衰减方法 (王伟等, 2012),其实质是使滤波操作窗能够最佳匹配信号的局部结构特征,解决地震剖面的随机噪声衰减和有效信号保真的问题,在增强反射同相轴的横向一致性的同时有效保持剖面内的地层边缘和细节特征,改善地震资料的品质;基于地震图像校准的共成像点道集增强技术 (李钟晓等, 2013) 其实质是实现不同偏移距道集在时间和空间上的逐点匹配对齐,增强共成像点道集中同相轴的一致性,提高叠加结果的分辨率;基于偏移成像道集的剩余静校正方法 (徐锦承和张剑锋, 2016),其实质是通过生成随炮点和检波点位置变化的偏移道集, 由偏移道集中直接拾取炮、检点的地表一致性剩余时差,避免复杂情况下同相轴归位不准确导致的剩余时差拾取误差等等.
地震资料实际处理中,各种一致性处理方法均有其适应性:有的方法是实现技术上的原因存在着一定的局限性,在经处理后的数据上仍有剩余时差存在;有的方法只校正剩余时差部分,而不考虑相位的畸变,以致使处理过的记录不能恢复到最理想的形态;有的方法的振幅保真度不够高.为了更好地解决数据的一致性问题,在大量试验处理研究的基础上,通过采用WesternGeoco公司Omega处理系统的非刚性匹配处理技术获得了较理想的处理效果.
2 非刚性匹配技术 2.1 方法原理及特点非刚性匹配方法的原理:相似于数码摄影用到的将几张照片“缝合”到一张全景照片里.它依赖于图像的相似和有局部光滑变形 (或位移) 场匹配,非刚性术语指的是该算法的局部匹配特性.该方法的基本模式是:预先使输入的地震数据体按基准的立方体的形状分解形成体积元素 (三维象素),并将为每个象素进行位移估计,属3维处理 (Nickel和Sonneland,1999).
地震数据的动态范围较大, 而且数据有正有负.为了应用数字图像处理技术, 必须首先进行数据格式的转换.采取保持原数据结构关系及其量化精度的转换方法, 其转换模型为
(1) |
式中, tm为时间采样点, 相当于图像中的行;k为测线标志, 相当于序列图像中的帧数;xn为同一测线的道标志, 每道相当于图像的每列;se(tm, xn, k) 为原地震资料在采样点 (tm, xn, k) 上的原始数据, 有正值, 也有负值;semin为原始地震数据最小值;f(m, n, k) 为转换后新的灰值图像在像素 (m, n, k) 上的数据, 无负值.
根据公式 (1) 可以将实际的叠前和叠后三维地震数据体转换成数字图像处理技术中能够使用的新的三维地震数据体 (陈凤等,2003).
非刚性匹配方法的特点和技术优势分析:
1) 由于常规的一致性处理方法,如地表一致性剩余静校正处理中只能将属于地表一致性的时差问题解决,而对于非地表一致性的时差问题则无能为力.如果采用非地表一致性剩余静校正方法,虽然可以解决一部分剩余时差问题,但由于其方法原理是基于道集中全部道记录叠加生成的模型道,只进行时差方面的校正,不能进行波形的校正.而非刚性匹配方法既能解决剩余时差问题,也能对波形进行一致性修正,且不改变振幅相对大小,有利保幅.
2) 非刚性匹配方法是基于三维像素的,故运算给出的结果更稳定可靠.
3) 非刚性匹配方法的模型道是基于近炮检距的,有利高分辨率保持,消除速度精度不够对远炮检距的时差的影响,同时也可消除动校拉伸对波形畸变的影响.
总之,该方法能使得共中心点 (CMP) 道集或叠前偏移得到的共反射点 (CRP) 道集更加平直,同一波组的波形趋于一致,且振幅保真度高.
2.2 实现过程非刚性匹配技术的应用是在预处理的基础上,首先将叠前数据分选成炮检距数据,优选部分炮检距数据形成模型,然后对每一炮检距数据逐一与模型进行非刚性匹配,实现数据的一致性处理.非刚性匹配处理技术的关键是建立高质量的模型数据,我们在实际的地震资料处理中,在大量试验的基础上,探索出了非刚性匹配处理流程 (如图 1所示),其步骤如下:
1) 叠前数据处理包括:解编、空间属性定义、近地表静校正、去噪、振幅恢复、反褶积、CMP道集分选、速度分析、剩余静校正、叠前偏移等处理.
2) 将道集数据进行规则化、分选为共炮检距域数据.
3) 选取部分近炮检距数据叠加,建立匹配的模型数据.数据选取以不降低分辨率且具有较高信噪比为原则,一般取覆盖次数的10%~20%道记录即可.
4) 应用非刚性匹配法对每一炮检距数据逐一进行处理.
5) 将地震记录再抽回到道集,进行后续的相关处理.
3 应用效果分析图 2给出了实际资料经过非刚性匹配处理前后的CRP道集记录,从中可以看出:非刚性匹配处理后的记录 (图 2b) 较非刚性匹配处理前的记录 (图 2a) 的一致性有明显的提高,消除了剩余时差对CRP道集的影响 (箭头所指处).
图 3给出了非刚性匹配处理前后的叠加剖面,可见:在没做非刚性匹配处理 (图 3a) 的情况下,剖面的品质不好;在做过非刚性匹配处理后 (图 3b) 剖面品质取得了明显改善 (箭头所指处).
图 4给出了CRP道集非刚性匹配前后波形与振幅曲线对比图,可见:在没做非刚性匹配处理的情况下,波组的波形一致性不好 (图 4a);在做过非刚性匹配处理后,波组的波形一致性取得了明显改善 (图 4c).另外,非刚性匹配处理前 (图 4b) 后 (图 4d) 的振幅曲线没有发生明显改变,振幅的强弱关系保持良好,说明非刚性匹配法的振幅保真度高.
图 5给出了非刚性匹配处理前后偏移处理成果分别求取的频谱叠合显示,可见非刚性匹配处理后的频带较未做非刚性匹配处理的频带稍宽2 Hz,说明非刚性匹配处理在保持分辨率不降低方面也是有利的.
通过对实际资料的非刚性匹配法的处理效果对比,我们可以得到如下结论:非刚性匹配技术能有效地提高数据的一致性,且振幅保真度高,有利于提高资料处理成果的信噪比,也有利于保持分辨率不降低,是改善资料处理品质的有效处理手段之一.依据该技术的特点,该技术可适用于时延地震匹配处理、叠前角道集拉平、叠前CMP/CRP道集拉平、振幅随角度变化 (AVA) 处理、PS波与PP波同相轴匹配处理等.
致谢 感谢审稿专家提出的修改意见和编辑部的大力支持![] | An X F, Liu P, Zhang X. 2010. Application of surface consistent deconvolution to seismic prospecting and its Effect[J]. Progress in Geophysics (in Chinese), 25(6): 2125–2129. DOI:10.3969/j.issn.1004-2903.2010.06.032 |
[] | Chen F, Li J Z, Huang J M, et al. 2003. Application research of digital image processing method for improving signal to noise ratio of seismic section image[J]. Progress in Geophysics (in Chinese), 18(4): 758–764. DOI:10.3969/j.issn.1004-2903.2003.04.030 |
[] | Guo N M, Wu G C. 2012. Forward-modeling and feasibility study of non-repeating acquired time-lapse seismic exploration[J]. Progress in Geophysics (in Chinese), 27(1): 232–245. DOI:10.6038/j.issn.1004-2903.2012.01.026 |
[] | Li S J, Shi X J, Zheng H M, et al. 2002. The consistent correction of seismic amplitude in complicated surface area[J]. Chinese Journal of Geophysics (in Chinese), 45(6): 862–869. DOI:10.3321/j.issn:0001-5733.2002.06.014 |
[] | 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 |
[] | 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. |
[] | Liu Y K, Chang X, Wang H, et al. 2001. Estimation of near-surface velocity and seismic tomographic static corrections[J]. Chinese Journal of Geophysics (in Chinese), 44(2): 272–278. DOI:10.3321/j.issn:0001-5733.2001.02.015 |
[] | Ma Y J, Wang J. 2010. An improved residual normal moveout correction method in time domain[J]. Geophysical Prospecting for Petroleum (in Chinese), 49(3): 245–247. |
[] | Nickel M, Sonneland L. 1999. Non-Rigid matching of migrated time-lapse seismic[C]. SEG Annual Meeting. |
[] | Pei J Y, Zhang L Y, Wang L N, et al. 2011. Multi-survey joint pre-stack time migration processing study of deep beds in Songliao basin[J]. Chinese Journal of Geophysics (in Chinese), 54(2): 294–303. DOI:10.3969/j.issn.0001-5733.2011.02.004 |
[] | Wang W, Gao J H, Chen W C, et al. 2012. Random seismic noise suppression via structure-adaptive median filter[J]. Chinese Journal of Geophysics (in Chinese), 55(5): 1732–1741. DOI:10.6038/j.issn.0001-5733.2012.05.030 |
[] | Wang X, Chen X R. 2003. Offset division normal moveout correction[J]. Geophysical Prospecting for Petroleum (in chinese), 42(2): 212–214. |
[] | Wang X W, Hu Z D, Tian Y C, et al. 2006. Study of seismic wavelet processing method in two step deconvolution[J]. Progress in Geophysics (in Chinese), 21(4): 1167–1179. DOI:10.3969/j.issn.1004-2903.2006.04.018 |
[] | 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(3): 204–208. DOI:10.3969/j.issn.1004-2903.2006.01.032 |
[] | Wei G Y, Li Z Y, Xue H. 2009. Application of seismic attributes in residual static correction[J]. Oil Geophysical Prospecting (in Chinese), 44(S1): 10–12. |
[] | Wu D L. 2004. The application of matching filter method in cannon offset data processing[J]. Geophysical Prospecting for Petroleum (in Chinese), 43(6): 599–601. |
[] | Wu Q, Yang C C, Zhang W Z, et al. 2008. Research on processing technology for matching seismic data acquired in seashore and marine prospect[J]. Progress in Geophysics (in Chinese), 23(3): 761–767. |
[] | Xia H R, Zhou K M, Chen D G, et al. 2002. Correction of colored spectrum in time frequency domain[J]. Geophysical Prospecting for Petroleum (in Chinese), 41(3): 304–308, 311. |
[] | Xia H R, Zhou K M, Chen D G, et al. 2003. Analysis of resolution-preserved processing techniques[J]. Oil Geophysical Prospecting (in Chinese), 38(S): 23–26, 30. |
[] | Xu J C, Zhang J F. 2016. Residual static correction based on migrated gathers[J]. Chinese Journal of Geophysics (in Chinese), 59(2): 655–663. DOI:10.6038/cjg20160222 |
[] | Yuan C F, Yu Q F, Wang Y C, et al. 2002. Estimation and removal of residual normal moveout[J]. Geophysical Prospecting for Petroleum (in Chinese), 41(3): 317–320. |
[] | Yun M H, Ding W, Wang K Y, et al. 2006. Study and primary application on consistency processing of seismic data[J]. Geophysical Prospecting for Petroleum (in Chinese), 45(1): 65–69. |
[] | Zhao L Z, Wang K B, Dai X Y. 2011. Application of VTI anisotropic pre-stack depth migration technique for 3-D seismic data processing in JZ area[J]. Geophysical Prospecting for Petroleum (in Chinese), 50(2): 201–205. |
[] | Zhuang D H, Xiao C Y, Xu Y, et al. 1999. 4D seismic data processing and its key[J]. Progress in Geophysics (in Chinese), 14(2): 33–43. |
[] | Zhu J M, Li L, Yang W, et al. 2013. Multi-azimuth seismic data applied analysis on the exploration and development of Wenchang sag[J]. Progress in Geophysics (in Chinese), 28(5): 2587–2596. DOI:10.6038/pg20130538 |
[] | 安西峰, 刘平, 张旭. 2010. 地表一致性反褶积在地震勘探中的应用及效果[J]. 地球物理学进展, 25(6): 2125–2129. DOI:10.3969/j.issn.1004-2903.2010.06.032 |
[] | 陈凤, 李金宗, 黄建明, 等. 2003. 利用数字图像处理技术提高地震剖面图像信噪比[J]. 地球物理学进展, 18(4): 758–764. DOI:10.3969/j.issn.1004-2903.2003.04.030 |
[] | 郭念民, 吴国忱. 2012. 非重复采集时移地震正演模拟及可行性分析[J]. 地球物理学进展, 27(1): 232–245. DOI:10.6038/j.issn.1004-2903.2012.01.026 |
[] | 李生杰, 施行觉, 郑鸿明, 等. 2002. 复杂地表条件反射振幅一致性校正[J]. 地球物理学报, 45(6): 862–869. DOI:10.3321/j.issn:0001-5733.2002.06.014 |
[] | 李钟晓, 陆文凯, 王季. 2013. 基于地震图像校准的共成像点道集增强技术[J]. 地球物理学报, 56(5): 1592–1608. DOI:10.6038/cjg20130517 |
[] | 林伯香, 孙建国. 2001. 相位替换法剩余时差校正[J]. 石油物探, 40(3): 15–22. |
[] | 刘伊克, 常旭, 王辉, 等. 2001. 三维复杂地形近地表速度估算及地震层析静校正[J]. 地球物理学报, 44(2): 272–278. DOI:10.3321/j.issn:0001-5733.2001.02.015 |
[] | 马永军, 王季. 2010. 一种改进的时间域剩余动校正方法[J]. 石油物探, 49(3): 245–247. |
[] | Nickel M, Sonneland L. 1999.非刚性匹配的延时地震偏移[C]. SEG年会. |
[] | 裴江云, 张丽艳, 王丽娜, 等. 2011. 松辽盆地深层地震资料叠前时间偏移连片处理技术研究[J]. 地球物理学报, 54(2): 294–303. DOI:10.3969/j.issn.0001-5733.2011.02.004 |
[] | 王伟, 高静怀, 陈文超, 等. 2012. 基于结构自适应中值滤波器的随机噪声衰减方法[J]. 地球物理学报, 55(5): 1732–1741. DOI:10.6038/j.issn.0001-5733.2012.05.030 |
[] | 王鑫, 陈新荣. 2003. 分偏移距动校正技术[J]. 石油物探, 42(2): 212–214. |
[] | 王西文, 胡自多, 田彦灿, 等. 2006. 地震子波处理的二步法反褶积方法研究[J]. 地球物理学进展, 21(4): 1167–1179. DOI:10.3969/j.issn.1004-2903.2006.04.018 |
[] | 王云专, 杨立伟, 李素华. 2006. 剩余时差校正及泊松比反演[J]. 地球物理学进展, 21(1): 204–208. DOI:10.3969/j.issn.1004-2903.2006.01.032 |
[] | 魏庚雨, 李振勇, 薛慧. 2009. 利用地震属性求取剩余静校正的方法[J]. 石油地球物理勘探, 44(S1): 10–12. |
[] | 邬达理. 2004. 匹配滤波法在大炮检距资料处理中的应用[J]. 石油物探, 43(6): 599–601. |
[] | 吴琼, 杨长春, 张文忠, 等. 2008. 过渡带拼接地震资料处理方法研究[J]. 地球物理学进展, 23(3): 761–767. |
[] | 夏洪瑞, 周开明, 陈德刚, 等. 2002. 时频域有色谱校正[J]. 石油物探, 41(3): 304–308, 311. |
[] | 夏洪瑞, 周开明, 陈德刚, 等. 2003. 保持分辨率处理方法分析[J]. 石油地球物理勘探, 38(S): 23–26, 30. |
[] | 徐锦承, 张剑锋. 2016. 基于偏移成像道集的剩余静校正方法[J]. 地球物理学报, 59(2): 655–663. DOI:10.6038/cjg20160222 |
[] | 苑春方, 余钦范, 王彦春, 等. 2002. 动校正剩余时差的估计与校正[J]. 石油物探, 41(3): 317–320. |
[] | 云美厚, 丁伟, 王开燕, 等. 2006. 地震资料一致性处理方法研究与初步应用[J]. 石油物探, 45(1): 65–69. |
[] | 赵玲芝, 王克斌, 戴晓云. 2011. VTI各向异性叠前深度偏移技术应用研究-以JZ工区三维地震资料处理为例[J]. 石油物探, 50(2): 201–205. |
[] | 庄东海, 肖春燕, 许云, 等. 1999. 四维地震资料处理及其关键[J]. 地球物理学进展, 14(2): 33–43. |
[] | 朱江梅, 李列, 杨薇, 等. 2013. 多方位角地震资料在文昌凹陷勘探开发中的应用分析[J]. 地球物理学进展, 28(5): 2587–2596. DOI:10.6038/pg20130538 |