2. 国家海洋局第二海洋研究所 海底科学重点实验室, 杭州 310012
2. Second Institute of Oceanography, SOA, Hangzhou 310012, China
短排列无定位拖缆多道地震勘探在近海浅水工程地质勘察中应用广泛,该方法不仅可以获取海底的地层展布、沉积特征及构造形态等地质情况,还可以揭示海底下可能存在的灾害地质(如活动断层、浅层气、麻坑、泥火山与海底滑坡等).而实际应用中,无定位拖缆多道地震勘探在资料处理方面还存在多个急需解决的问题.
首先,海上多道地震勘探面临着严重的多次波干扰,特别是在浅水硬质海底区,自由表面相关多次波干扰尤为明显,严重影响多道地震资料的处理与解释.针对深水长排列油气勘探中的自由表面相关多次波,荷兰Delft大学Berkhout与Verschuur提出了SRME有效压制方法(Verschuur et al., 1992;Berkhout and Verschuur, 1997;Berkhout, 1999),该方法基于实际数据,通过波动方程多次波预测与匹配相减两过程(李鹏等,2007;李学聪等,2010;石颖等,2010;宋家文等,2014),可以有效压制深水长排列多道地震记录中的多次波.但该方法理论上不适用于浅水长排列多次波压制,潘军等(2015)详细分析了其原因,认为勘探水深与排列长度的比值大小,决定了SRME压制多次波的效果,因浅水长排列其比值太大(超过100),SRME方法压制效果就差.然而,针对小偏移距短接收排列,其长度较短,相应比值也小,SRME方法可以有效应用(徐强等,2015).Moore和Bisley (2006)提出水深与最小偏移距的比值及其反射数据能量的一致性(无相互干涉同相轴)也决定了SRME方法在浅水短排列中的应用效果.
浅水短排列地震勘探满足了水深与排列长度及最小偏移距比值大小的要求,SRME方法能否有效应用还受其他因素影响,如无定位拖缆野外观测系统的准确定义、浅水直达波与折射波干涉、近偏移距数据缺失等影响.浅水直达波与折射波对海底反射同相轴的干涉,可在τ-p域内通过谱图有效切除来解决(Spitz,1999;李列等,2015).近偏移距反射数据缺失,可通过小偏移距高精度抛物Radon变换插值来解决(Kabir and Verschuur, 1995).然而,无定位拖缆野外观测系统的准确定义及CMP位置有效归算问题,目前还没有可靠方法,现有文献只介绍定位拖缆浅水勘探的野外观测系统定义及多次波压制情况(徐锦玺等,2005;周健等,2012;孟庆生等,2007),少有介绍无定位拖缆野外观测系统定义的问题.针对短排列无定位拖缆多道地震勘探,实际工作中常采用完全直线观测假设,简单定义野外勘探的炮点和接收点空间关系.基于直线假设的野外观测系统平面位置偏差较大,距实际勘探航迹偏离较远,严重影响CMP抽道集及位置归算精度,制约了SRME方法对浅水多次波的压制效果及反射层位的成像精度.
2 无定位拖缆野外观测系统定义本文以舟山东极岛海域30 m水深范围内的小道距短排列多道地震勘探数据为例,分析常规方法定义的无定位拖缆观测系统缺陷,由此提出一种新的计算方法,并以SRME方法压制浅水多次波的效果,说明新方法的可靠性.
2.1 基于直线假设定义海上长排列油气多道地震勘探一般采用大道距(12.5 m以上)长排列(120道以上)接收,在接收缆的首尾以及中间段按固定间隔安装有磁罗经、声学鸟、深度传感器与尾标等定位设备,通过这些设备,野外勘探时以勘探船上的GPS天线为参考点,可以实时计算出各接收道实际经纬度的平面投影坐标,资料处理时利用这些准确坐标可以精确定义野外观测系统参数.
与深水长排列油气勘探不同,浅水短排列(24道12.5 m以内道间距)地震勘探一般不在接收缆上安装声学定深、罗经与尾标等定位设备,各接收道的实际坐标现场无法获取,常规资料处理常采用人为假设方法来定义野外观测系统参数.具体方法是人工假设第一炮最大偏移距接收道为坐标0点,再根据炮间距、最小偏移距和道间距等参数,按照完全直线激发接收的假设,沿直线航向依次计算各炮点与接收点的相对位置,最后利用这些相对位置简单生成野外观测系统图.各接收点与炮点的具体计算公式为:
(1) |
式中,Lri、Lsn为接收点与炮点相对直线长度,Imax为排列总道数,图 1中取值24,ΔX为道间距,Sn为炮序号,ΔS为炮间距,δmin为最小偏移距.通过(1) 式可以简单计算出所有炮点与检波点相对于假设0点的直线距离,而不是实际坐标值,用相对距离定义的观测系统无法计算各炮点与接收点的实际X、Y坐标,会影响后续CMP抽道集运算,并给共反射点叠加次数与位置归算带来偏差.直线假设简单定义观测系统示意如图 1.
利用上面方法及公式(1),对舟山东极岛海域多道地震勘探所采用的短排列无定位拖缆作假定观测系统定义,其结果如图 2所示,图 2中圆点为导航软件控制激发并记录的实际炮点位置轨迹,而灰色方格直线为采用上述假设方法计算的CMP位置与叠加次数图.从图 2可以看出,完全直线假设生成观测系统方法虽然简单易行,但基于直线假设的CMP位置与实际炮点轨迹偏差很大.由于实际勘探船是沿设计测线左右偏航跑线(即羽角的存在),且船速不稳定,若等时间间隔放炮,炮间距就不是假设中的直线等间距规则炮点,即使等距离放炮,完全直线假设的勘探排列与实际情况也相差甚远,造成CMP抽道集处理及位置归算不准,严重影响后续的数据处理精度与成像效果,特别是多次波压制处理.另外,由于无法准确计算实际大地坐标值,图 1假设方法计算的CMP点无法获得图 2中的实际X、Y坐标值(为方便对比,图 2将其放置于同一图中),这给剖面的高精度解释带来不利.
浅水无定位拖缆野外勘探时,震源炮点的激发是由导航软件通过等距或等时设定来控制的,激发后,炮点的实际DGPS (Differential Global Positioning System,差分全球定位系统)位置由导航软件实时记录存储,利用存储的炮点位置文件,我们可以生成图 2中所示的精确炮点轨迹激发图.但浅水短排列拖缆因未装配声学定位设备,现场无法获得各接收道的具体位置,也很难计算其拖行状态,唯一可利用的就是炮点轨迹坐标.由于浅水短排列道间距小,接收道数少,在海况较好的情况下,调查船正常走航勘探时,接收拖缆基本沿炮点轨迹拖行.因此,利用上面实际的炮点轨迹,我们可以模拟并计算出各接收道在炮点轨迹上的空间位置及具体坐标,利用这些坐标可以实现对野外观测系统的准确定义.
具体计算分两步实现,首先计算无定位拖缆各接收道相对于炮点轨迹的空间位置关系,其次计算各接收道的实际X、Y坐标值.
2.2.1 接收道空间关系计算震源炮点实际轨迹坐标是计算空间位置关系的基础,为了获得可靠的高精度位置,我们可以将GPS天线安装在船尾并靠近震源的安全位置,量测好震源枪距GPS天线的距离,同时在导航软件中,将震源所在位置设为导航定位的中心,由导航软件实时获取并记录震源点的实时坐标,有了这些坐标,采用如下的运算公式,可以获得无定位拖缆相对于震源炮点的空间位置.
(2) |
式中,Dri为待求接收道与当前排列震源炮点间距离,i为排列上道序号,n表示待求接收道同震源间的距离Dri是炮间距ΔS的多少倍,由n值可得炮点轨迹上离待求接收道最近的震源炮点位置,Q为待求接收道在震源排列上的方位角,Yn为待求接收道附近最近的震源炮点坐标,Ys为当前排列震源炮点坐标,Dsn为两炮点间距.利用震源炮点实际轨迹坐标,再由上面公式中的变量n与Q,可以计算出各接收道在炮点轨迹上的合理空间位置.
2.2.2 接收道实际坐标计算根据上面相对空间关系计算公式(2),利用实际炮点轨迹,我们可以模拟接收排列在勘探船后拖行的延伸形态,并得到各接收道在炮点轨迹中的相对空间位置,然后通过震源炮点实际坐标值,结合反距离比线性加权插值算法,可以计算出各接收道的实际坐标值.
反距离比反映了距离接收道越远的炮点坐标,对当前接收道坐标值的贡献越小,实际意义表明越远的炮点位置,勘探船对接收缆的拖拉力越小,对拖缆的拖行形态影响就小,综合接收道附近多个不同位置的炮点,可以模拟勘探船沿测线走航时接收缆的拖行状态,并计算其具体位置.拖缆中各接收道具体坐标计算公式为:
(3) |
式中,xri、yri为待求接收道实际坐标,Xsn、Ysn为接收道附近参与计算的炮点坐标,Drisn为待求接收道与炮点间距离,Ds1sN为参与计算的首尾两炮点间距离,N为参与计算的总炮数,且N≥2.从公式(3) 中可以看出,距离待求接收点越近的炮点,其距离比系数越大,对接收道实际坐标的计算贡献越大,再综合利用附近多个炮点的坐标值,通过反距离比线性插值,可以较准确地解算出拖缆中各接收道的实际坐标值.
综合公式(2) 及公式(3),通过编制相应的运算程序,可以计算出观测系统定义中所需要的接收点实际坐标值,结合各炮点实际坐标,将其应用于上述舟山东极岛实际地震勘探数据中.野外勘探时采用的是24道12.5 m道间距固态拖缆,震源为210 in3 GI气枪(Generator and Injector,生成与注射两次充气式气枪),炮间距25 m,最小偏移距40 m.按照理论计算,该反射地震CDP覆盖次数平均为6次.利用公式(2) 与(3),依据炮点轨迹坐标,计算出某炮接收排列的空间位置如图 3所示.图 3中933号圆点为当前待计算排列中炮点,灰色三角形为拖缆中待求解24个接收道位置,求解时,N值取3,即待求解接收道前取2个炮点,后取1个炮点参与计算,通过反距离比线性插值算法,计算出图 3所示的各接收道相对于炮点轨迹的空间关系及其实际坐标图.最后利用专业处理软件导入这些实际坐标,定义对应的观测系统图,专业软件会自动计算出共反射点的实际位置与叠加次数,并绘制如图 4所示的共中心点覆盖次数与排列关系图.从图 4可以看出,采用新方法定义的观测系统,获得的共中心点位置符合弯曲测线实际勘探航迹,共中心点覆盖次数也比较均匀,未出现图 2中出现的较大偏差.
SRME方法是压制海水面等自由表面相关多次波的首选方法,它能够有效压制震源相关类多次波(Lokshtanov, 2000),压制过程中不需要任何的层位速度与深度等信息,技术发展已很成熟,国内外发表了大量相关文献介绍其在深水长排列地震数据中的应用,然而在浅水无定位拖缆中的应用介绍非常少.浅水勘探鸣震类自由表面相关多次波干扰尤为明显,特别是在硬质海底(如基岩面、砂土层等)区,强能量多次波会严重干扰海底下地层有效反射,模糊有效层位成像,影响浅水勘探效果.有效压制该类多次波是目前浅水多道地震勘探资料处理中的重点和难点问题.
由于无定位拖缆在勘探现场无法直接获取各接收道的实际位置坐标,采用上述常规的基于完全直线假设简单定义观测系统后,再应用SRME方法作浅水自由表面相关多次波压制,其效果往往比较差.为提高SRME方法的应用效果,改进浅水区无定位拖缆地震数据的成像质量与位置精度,我们采用文中设计的新方法进行准确的野外观测系统定义,再进行同样方法的多次波压制.
首先计算无定位拖缆中各接收道相对于震源轨迹的空间方位和附近小距离震源点位置,再通过反距离比线性插值算法,利用上面公式计算各接收道的实际坐标,综合所有震源点和接收点的实际坐标,定义野外观测系统相关参数,最后利用专业处理软件,采用SRME方法压制自由表面相关多次波,在同样的参数设置下,新方法定义的观测系统数据获得了比较满意的压制效果.下面利用上述两类观测系统定义的实际多道反射地震数据,根据叠加剖面有效反射层位的成像效果,说明采用两类不同方法定义野外观测系统后,SRME方法压制浅水鸣震类自由表面相关多次波的剖面效果.
利用图 2与图 4两套观测系统计算的CMP叠加次数及对应位置,文中采用相同参数的SRME方法作多次波压制处理,最后经共反射点叠加成像,获得了对应图 5和图 6两套叠加剖面成果图.从两图可以明显看出,因直线假设造成CMP抽道计算误差及位置偏差,图 5定义的观测系统数据其多次波压制效果欠佳,叠加剖面上鸣震类海底多次波反射干扰依然明显(B1处),海底下留存的虚震源类多次波依然严重干扰有效倾斜地层反射(M2、M4处),重点探测的基底反射同相轴能量模糊不连续,强能量海底鸣震处的基底反射中断(B1、B2处),基底虚震源类多次波能量依然比较明显(M1、M3处).相对应的图 6中,海底鸣震类多次波能量衰减较好(B1处),海底下倾斜地层反射同相轴能量突出,重点探测的基底反射同相轴连续且能量清晰,虚震源类多次波压制较好,叠加剖面反射地层成像明显优于图 5.
另外,因直线假设定义的观测系统无法获得实际坐标,经抽道集叠加处理后,图 5剖面中显示的X、Y坐标都为0;新方法根据实际炮点坐标,计算出了各接收道的实际X、Y位置,通过定义准确的观测系统,再利用商业软件作常规处理后,图 6中显示了准确的X、Y坐标值.
因受海底浅层气对纵波能量的屏蔽影响,两图中都有三处反射紊乱区,对应正下方的基底反射能量也受到干扰,使基底反射同相轴中断(G1、G2、G3处).对于浅层气的屏蔽问题,理论上采用数字处理方法无法解决,油气勘探中一般采用纵横波联合勘探的方法来避开气层屏蔽,清晰揭示气层及其下方的有效反射同相轴.
4 结论野外观测系统定义是多道地震勘探工作的重点,也是后续所有处理方法成功应用的基础.浅水无定位拖缆多道地震勘探因野外无法直接获取各接收点的实际坐标,影响后续观测系统的准确定义.采用直线假设简单定义观测系统是目前比较通用的手段,但该方法定义的位置与实际弯曲测线不符,计算出的共中心点叠加次数也有偏差,严重影响后续资料处理.为避免常规方法缺陷,本文以计算无定位拖缆实际坐标为出发点,提出一种利用炮点轨迹计算接收点实际坐标的新方法.该方法分为两步,首先计算各接收道相对于炮点轨迹的方位及相对位置,其次利用炮点坐标,采用反距离比线性加权插值算法计算各接收道的实际坐标,综合这些实际坐标,经野外观测系统准确定义后,可以获得符合野外弯曲测线勘探实际的观测系统图.最后利用两类不同观测系统定义数据,对浅水无定位短排列多道地震数据作SRME海底鸣震与虚震源类多次波压制处理,结果表明采用文中新方法所定义的观测系统数据,其处理效果明显优于基于传统直线假设的观测系统数据.经实际浅水地震勘探数据应用证明,新方法不仅避免了传统假设方法的缺陷,还提高了叠加剖面成像效果及位置精度,为后续地质解释提供了基础.
致谢感谢野外数据采集团队提供的实际勘探数据,感谢编辑对本文修改给予的帮助,特别感谢匿名评审专家提出的宝贵详细改进意见!
Berkhout A J, Verschuur D J. 1997. Estimation of multiple scattering by iterative inversion, Part I:Theoretical considerations. Geophysics, 62(5): 1586-1595. DOI:10.1190/1.1444261 | |
Berkhout A J. 1999. Multiple removal based on the feedback model. The Leading Edge, 18(1): 127-131. DOI:10.1190/1.1438140 | |
Kabir M M N, Verschuur D J. 1995. Restoration of missing offsets by parabolic radon transform. Geophysical Prospecting, 43(3): 347-368. DOI:10.1111/gpr.1995.43.issue-3 | |
Li L, Xie Y H, Li Z N, et al. 2015. Research progress on offshore multiple attenuation and imaging. Progress in Geophysics, 30(1): 446-453. DOI:10.6038/pg20150166 | |
Li P, Liu Y K, Chang X, et al. 2007. Application of the equipoise pseudomultichannel matching filter in multiple elimination using wave-equation method. Chinese Journal of Geophysics, 50(6): 1844-1853. DOI:10.3321/j.issn:0001-5733.2007.06.027 | |
Li X C, Liu Y K, Chang X, et al. 2010. The adaptive subtraction of multiple using the equipoise multichannel L1-norm matching. Chinese Journal of Geophysics, 53(4): 963-973. DOI:10.3969/j.issn.0001-5733.2010.04.021 | |
Lokshtanov D. 2000. Suppression of water-layer multiples-from deconvolution to wave-equation approach.//70th Meeting Society of Exploration Geophysicists. Expanded Abstract, 1981-1984. | |
Meng Q S, Chu X F, Guo X J, et al. 2007. The application of high resolution seismic data processing technique in multi-channel shallow offshore engineering seismic surveys. Progress in Geophysics, 22(3): 1006-1010. DOI:10.3969/j.issn.1004-2903.2007.03.053 | |
Moore I, Bisley R. 2006. Multiple attenuation in shallow-water situations.//68th Meeting, EAGE. Expanded Abstracts, F018. | |
Pan J, Luan X W, Sun Y B, et al. 2015. The application of SRME in high resolution seismic exploration of shallow sea environment. Progress in Geophysics, 30(1): 429-434. DOI:10.6038/pg20150163 | |
Shi Y, Liu H, Zou Z. 2010. Surface-related multiples prediction based on wave equation and adaptive subtraction investigation. Chinese Journal of Geophysics, 53(7): 1716-1724. DOI:10.3969/j.issn.0001-5733.2010.07.023 | |
Shi Y, Wang W H. 2012. Surface-related multiple suppression approach by combining wave equation prediction and hyperbolic Radon transform. Chinese Journal of Geophysics, 55(9): 3115-3125. DOI:10.6038/j.issn.0001-5733.2012.09.029 | |
Song J W, Verschuur D J, Chen X H. 2014. Research status and progress in multiple elimination. Progress in Geophysics, 29(1): 240-247. DOI:10.6038/pg20140134 | |
Spitz S. 1999. Pattern recognition, spatial predictability, and subtraction of multiple events. The Leading Edge, 18(1): 55-58. DOI:10.1190/1.1438154 | |
Verschuur D J, Berkhout A J, Wapenaar C P A. 1992. Adaptive surface-related multiple elimination. Geophysics, 57(9): 1166-1177. DOI:10.1190/1.1443330 | |
Verschuur D J, Berkhout A J. 1997. Estimation of multiple scattering by iterative inversion, Part Ⅱ:Practical aspects and examples. Geophysics, 62(5): 1596-1611. DOI:10.1190/1.1444262 | |
Xu J X. 2005. Seismic exploration acquisition technique application of Beach shallow sea. Progress in Geophysics, 20(1): 66-70. DOI:10.3969/j.issn.1004-2903.2005.01.013 | |
Xu Q, Wang Z, Shi Z Y, et al. 2015. Shallow water multiple attenuation. Oil Geophysical Prospecting, 50(2): 238-242. | |
Zhou J, Wang H Y, Zhu L J. 2012. Research of seismic acquisition technique in beach shallow sea.//Oil Geophysical Prospecting, Symposium on Geophysical Prospecting Technology (in Chinese). Harbin:Chinese Petroleum Society, Chinese Geophysical Society, 398-401. | |
李列, 谢玉洪, 李志娜, 等. 2015. 海上多次波压制与成像方法研究进展. 地球物理学进展, 30(1): 446–453. DOI:10.6038/pg20150166 | |
李鹏, 刘伊克, 常旭, 等. 2007. 均衡拟多道匹配滤波法在波动方程法压制多次波中的应用. 地球物理学报, 50(6): 1844–1853. DOI:10.3321/j.issn:0001-5733.2007.06.027 | |
李学聪, 刘伊克, 常旭, 等. 2010. 均衡多道1范数匹配多次波衰减的方法与应用研究. 地球物理学报, 53(4): 963–973. DOI:10.3969/j.issn.0001-5733.2010.04.021 | |
孟庆生, 楚贤峰, 郭秀军. 2007. 高分辨率数据处理技术在近海工程地震勘探中的应用. 地球物理学进展, 22(3): 1006–1010. DOI:10.3969/j.issn.1004-2903.2007.03.053 | |
潘军, 栾锡武, 孙运宝, 等. 2015. SRME技术在海洋浅水高分辨率地震勘探中的应用. 地球物理学进展, 30(1): 429–434. DOI:10.6038/pg20150163 | |
石颖, 刘洪, 邹振. 2010. 基于波动方程表面多次波预测与自适应相减方法研究. 地球物理学报, 53(7): 1716–1724. DOI:10.3969/j.issn.0001-5733.2010.07.023 | |
石颖, 王维红. 2012. 基于波动方程预测和双曲Radon变换联合压制表面多次波. 地球物理学报, 55(9): 3115–3125. DOI:10.6038/j.issn.0001-5733.2012.09.029 | |
宋家文, VerschuurD J, 陈小宏. 2014. 多次波压制的研究现状与进展. 地球物理学进展, 29(1): 240–247. DOI:10.6038/pg20140134 | |
徐锦玺. 2005. 滩浅海地震勘探采集技术应用. 地球物理学进展, 20(1): 66–70. DOI:10.3969/j.issn.1004-2903.2005.01.013 | |
徐强, 王征, 史增园, 等. 2015. 浅水多次波衰减. 石油地球物理勘探, 50(2): 238–242. | |
周健, 王沪钰, 朱鲁京. 2012. 滩浅海地震勘探采集技术方法研究. //中国石油学会2012年物探技术研讨会论文集. 哈尔滨: 中国石油学会, 中国地球物理学会, 398-401. | |