渤海湾X油田气云区地震资料关键处理技术研究

夏同星,刘　垒,明　君,周学锋

(中海石油(中国)有限公司天津分公司渤海石油研究院,天津塘沽300452)

气云区地震资料处理是地球物理界公认的世界级难题。根据地震波传播理论,由于气云区的存在,地震波能量被严重吸收衰减,地震资料信噪比较低,难以保障成像质量,同时,由于气云区粘滞气体空间分布不均匀,导致地层速度横向变化较大且无明显规律可循,给处理过程中高精度速度建模和偏移成像带来巨大挑战。

针对气云区地震资料的成像处理,通常解决思路有3种:①根据气云空间展布特征,设计合理观测系统,最大限度规避气云区对地震波传播的影响,包括采集方向、最大和最小偏移距、气枪震源容量等参数[1-2];②优化常规纵波成像处理流程,针对气云区对地震波吸收衰减等问题,开展针对性处理研究,最大限度提高气云区成像质量[3-6];③开展多波勘探及多波成像技术研究,气体的存在对体积模量影响严重,而对剪切模量影响则不大,因此,理论上可以通过横波勘探来解决气云区模糊带内的地震成像问题[7-13]。思路①和思路③都是从地震资料采集入手。本文以X油田现有地震数据为基础,按照思路②进行研究,通过多种关键技术的应用,最大限度提高气云区地震资料的成像质量。

1　油田情况及资料介绍

X油田位于郯庐断裂带的东支,渤南凸起带中—东段的结合部,是渤海湾盆地发现的特大型整装油田。

受郯庐断裂带活动影响,在油田主体区附近发育面积约20km2的气云区,从均方根振幅属性图上看到气云区呈NNE向条带状分布,气云区内地震波能量被严重吸收衰减(图1)。从过X-C井的地震剖面上可看出,气云非常发育,气云区模糊带范围较大,由浅至深地震剖面基本上无法成像(图2)。X油田X-C井钻井揭示纯油层共157.7m,但是由于该井位于气云区,难以利用地震资料对井区目的层构造和储层进行准确描述,从而影响了油田的高效开发。

图1　均方根振幅属性

图2　过X-C井地震剖面

2　关键处理技术

气云区成像难点主要表现在两个方面:一为地震波能量被严重吸收衰减,导致气云区内资料的信噪比较低,难以保障后续的成像质量;二为气体分布不均匀,导致速度横向变化且无明显规律可循,难以建立高精度速度模型。针对这两个问题,本文研究了Q值层析反演(Qtomo)能量补偿技术和潜波层析反演(DWT)浅层速度建模技术,同时辅以其它配套的处理技术应用,如匹配追踪插值(MPFI)数据规则化技术、水层多次波压制(DWD)+广义地表多次波压制(GSMP)等组合多次波压制技术等,设计形成气云区地震资料处理流程(图3),从而保障最终的成像质量。

2.1　Q值层析反演气云区能量吸收衰减补偿

气云区对纵波能量吸收衰减严重,直接影响气云区成像质量。如何对纵波能量进行有效补偿是处理的关键。

在偏移过程中采用Qtomo能量补偿技术对吸收衰减能量进行针对性补偿。补偿技术是一种比较符合地球物理实际规律的Q补偿方式,该技术主要利用层析方法对叠前地震数据进行等效Q分析,生成空变Q模型,然后在偏移过程中旅行时计算时更新Q场,在深度偏移过程中进行有效Q补偿。主要分为以下3步。

1) 利用偏移前地震数据估算衰减旅行时。目前比较典型的算法为地震数据属性分析方法,比如利用地震数据的振幅谱[14],或者利用频移特征[15]。应用过程中,可以根据地震数据的品质和方法的实用性采用合适的估算方法。

2) 利用衰减旅行时层析方法在深度偏移中生成空变Q场,这是Q值层析反演补偿过程中非常重要的一步。地震波通过粘弹性介质的吸收衰减可以表征为衰减旅行时的形式:

图3　气云区地震资料处理流程

(1)

式中:t*为衰减旅行时。对于给定的速度场,假设射线路径不受Q场影响,那么可以通过下面的线性系统来求解品质因子Q的分布。具体过程需要通过深度域的动力学射线追踪来实现。

(2)

3) 考虑吸收衰减设计深度偏移算子,在深度偏移过程中进行模型驱动的Q补偿[6,16]。

在网格层析的过程中,除了得到空变Q场,同时还生成速度模型,基于后续速度建模的考虑,该研究在第2)步中,采用DWT方法来实现衰减旅行时层析。实际补偿效果如图4所示,可看出,Qtomo较好地补偿了中深层能量,古潜山不整合面地震反射更为清晰。

2.2　潜波层析建立浅层速度模型

地震资料品质较好时,反射波层析方法可以生成较准确的速度模型。但是当地震资料品质较差时,如受气云区的影响,地震反射能量较弱、多次波干扰严重、浅层速度横向变化剧烈,反射波层析生成的速度模型存在较大不确定性。因此,在速度建模过程中,本文采用DWT技术和反射波层析技术相结合的速度建模方法。利用DWT技术建立浅层速度模型,将该模型作为反射波层析的初始速度模型来进一步反演中深层速度。

以前DWT方法主要用于处理过程中静校正量的计算[17]。在满足地层速度随深度增加逐渐增大的假设条件下,潜波在地震记录上以初至波的形式出现,表现特征相对来说清晰可靠,容易识别。因而利用潜波初至信息可以建立较好的浅部速度模型,克服了反射波层析反演在浅层速度建模中的不确定性。DWT反演是一个迭代的过程,它包括初始速度模型的建立、旅行时计算和迭代更新[18]。

初始速度模型的建立非常重要,一般而言,速度越平滑迭代越稳定,然而平滑时要综合考虑地质构造特征,速度模型太过平滑,偏离了真实速度场,那么在有限的迭代次数下无法收敛到真实的速度模型。

图4　Qtomo能量补偿前(a)、后(b)偏移叠加剖面对比

初至时间的拾取准确度直接影响着DWT计算结果的准确性。为了在背景噪声中能够有效地自动拾取初至信息,本次处理通过能量增益、随机噪声压制等方法提高气云区初至波能量。为了确保初至拾取质量,在初至拾取后,可以产生多种属性图,包括初至时间与偏移距关系、初至时间与炮点关系等,对初至信息进行有效质控。

得到可靠初至后,就可以进行DWT反演迭代。选取初始网格为200m×200m,并逐次递减,通过层析反演算法减小模型和实际数据的误差,得到可靠的近地表速度模型。由图5可看出,DWT迭代反演出的低速层范围和地震数据中低能量气云区域吻合得非常好,随着迭代次数增加,速度细节更为清晰;由图6速度切片也可看出低速信息主要集中在浅层,这因为气一般会向浅层集中,该现象在本工区尤为明显,大部分的气云都集中在浅部;图7对比了反演得到的地震波速度、声波速度和VSP层速度,可见,三者在1km以上深度范围趋势吻合较好。

在地震道集中,由于缺乏浅层反射信息,较难提取出准确的浅层速度,而利用DWT技术,得到的速度模型与测井速度趋势一致,为后续反射波层析反演速度建模提供了准确的浅层初始速度模型。

图5　DWT不同迭代次数速度模型对比a 初始速度模型; b 第一次迭代结果(200m×200m); c 第二次迭代结果(100m×100m); d 第三次迭代结果(50m×50m)

图6　DWT不同迭代次数速度模型切片对比(120m)a 初始速度模型; b 第一次迭代结果(200m×200m); c 第二次迭代结果(100m×100m); d 第三次迭代结果(50m×50m)

图7　DWT反演速度与3口井的声波速度和VSP层速度对比

2.3　其它配套处理技术研究

2.3.1MPFI数据规则化

分析X油田去噪后数据发现,该油田数据不规则,主要表现在:反射点位置不在面元中心、偏移距分布不规则、方位分布不规则[19-21],且浅层近道数据缺失严重。

常用的数据规则化方法是面元均化和加权补偿方法。面元均化技术通过向相邻面元借道补齐缺失道,破坏了地震道坐标的真实性;加权补偿法基于覆盖次数或者叠加纯波数据体能量求取比例因子,但不能弥补所有缺失的偏移距道数据。本文采用MPFI技术对数据进行规则化,该方法在数据规则化的同时,可有效避免空间假频干扰。首先通过傅里叶变换,将数据变换到f-x域。然后对每一个频率成分进行防漏频傅里叶变换[22-25]。不规则采样使得傅里叶变换的正交基函数不正交,使地震能量泄漏到其它的频率成分上,造成频谱失真。防漏频傅里叶变换方法应用非均匀傅里叶变换得到不规则数据的频谱,即:

(3)

假设造成频谱泄漏的是频谱中能量最大的系数成分,那么先从频谱中筛选出能量最大的频率成分,将其进行非均匀傅里叶反变换((4)式),再从原始的不规则数据中将分离出来的这个成分减掉((5)式),之后把余下的数据作为新输入继续分离,直到把所有的频率成分都分离出来,重构有效信号频谱。

为了防止假频,保证从频谱中筛选出的能量最大的频率成分为有效信号,可以根据f-k域主要信号特征,以低频数据的频谱为先验信息,对高频部分进行加权处理。最后将重构的频谱信息反变换到时空域,生成规则化后的地震数据。

对实际数据进行处理,MPFI规则化技术有效补偿了近偏移距缺失道(图8),同时从方差切片(图9)上可以看出,处理后的数据明显弱化了浅层采集脚印的影响,提高了信噪比。

图8　规则化前、后单偏移距剖面对比a MPFI插值前; b MPFI插值后

图9　规则化前、后方差切片对比a MPFI插值前; b MPFI插值后

2.3.2DWD+GSMP组合多次波压制

多次波压制是海上地震资料处理无法回避的问题,多次波产生类型多样,压制方法也不尽相同。X油田整体水深30m左右,由于水浅,直达波初至时间与海底反射波时间很接近,在地震资料中难以分辨海底反射信号;同时最小偏移距远大于水深,缺少近偏移距的海底反射信息。此时,直接采用预测相减法(如自由表面多次波压制方法,即SRME方法)不能有效衰减海底多次波[26-28]。因此,本文采用组合多次波压制技术,首先利用单缆鬼波去除技术(SSD)去除鬼波,在此基础上利用DWD+GSMP组合技术对多次波进行有效压制。

DWD是一项针对浅海环境开发的去除水层多次波技术,可以较好衰减与海底相关的短周期多次波。在τ-p域中,多次波周期Tm与线性动校正量P值的关系可以表示为:

(6)

式中:Tw为水底双程时间;vw为水速;Xref为参考偏移距。

与SRME方法相比,DWD方法对近道数据的要求更低,更易于实现,而且DWD方法基于模型驱动、对多次波进行非线性预测,能够较准确地得到一阶和高阶多次波振幅值,相比传统的τ-p反褶积方法它能够减少对有效波的伤害。

在实现过程中,首先通过近道自相关拾取海底反射时间,然后在τ-p域计算海底多次波周期,再进行自适应相减,最后通过τ-p反变换得到海底多次波模型。从图10可以看出,利用DWD方法处理后,多次波的能量明显得到了压制。值得注意的是,该方法只是设计用于预测和压制简单多次波和短程层间多次波,所以不能够实现其它复杂地表多次波或层间多次波的压制。

图10　DWD前、后叠加剖面对比a DWD前; b DWD后

三维GSMP基于SRME的基本思路,但不同的是GSMP是在三维数据空间通过更多的道集数据来预测多次波模型[29-30]。如图11所示,多次波在空间中具有明显的三维效应,除了来自震源和检波点的二维线方向外,还来自侧面反射,因此三维算法可以在一定孔径范围内找到了更多属于多次波的能量,建立的多次波模型更能代表复杂反射情况下的多次波,尤其是多次波模型受空间地质因素或是采集等影响而呈现复杂状态时,GSMP利用三维算法,可以模拟预测任意道数据在三维空间上的各种多次波,并进行自适应压制。

图11　多次波的三维效应

在具体计算时,首先选择目标道。这个目标道由接收点和炮点决定,然后根据工区情况决定孔径,对于孔径内的每个网格对炮点方向和检波点方向的道进行褶积,然后保存结果,所有网格计算完成后,叠加前面得到的褶积结果,然后进行下一道的计算。这样得到了多次波模型,进而在炮域或者共偏移距域进行适应性相减,得到去除多次波后的数据(图12)。

图12　GSMP前、后叠加剖面对比a GSMP前; b GSMP后; c 压制的多次波能量

3　应用效果分析

应用该流程对渤海X油田地震资料进行了处理,取得了较好效果。对比新老地震数据剖面发现,新处理数据较老资料在断层成像和气云区内成像上有较明显提高(图13):①新资料提高了花状断层的成像效果,在老地震剖面上,气云区内构造不成像,无法对断层进行合理解释,利用新地震资料可以对气云区内部的一些断层进行准确解释,如图13a所示;②新资料有效减小了气云对地震资料成像的影响范围,提高了气云内部成像质量,可以在一定程度上对气云区内的地层产状进行有效解释;③中深层的成像也得到了明显改善,尤其是对于基底界面,老资料几乎不成像,新资料成像效果有所改善,可以根据地震响应特征进行地震解释。

对比新老资料的方差切片(图14)也可以看到新资料改善了气云区周边断层的成像,虽然气云区内部仍无法准确成像,但是明显缩小了气云区的影响范围,为后续地震资料的解释提供了资料基础。

图13　新处理剖面(a)与老剖面(b)对比

图14　新处理资料(a)与老资料(b)方差体水平切片对比

4　结论

研究表明,通过优化常规纵波处理流程,对气云区地震资料开展针对性处理,可以明显改善成像质量,提高处理成果的资料品质。该处理流程在渤海X油田的应用中取得了较好的效果。

1) 能量补偿是气云区处理的关键环节,基于Qtomo的能量补偿是一种比较符合地球物理实际规律的Q补偿方法,该方法在偏移过程中,对气云区地震波能量进行有效补偿。

2) DWT反演速度建模技术利用初至波特征清晰可靠、容易识别的特点,可以在地震资料信噪比较差时,建立准确的浅层速度模型,为反射波层析反演速度建模提供准确的浅层初始速度。本文联合利用DWT和反射波层析进行速度建模,取得了较好效果。

3) 配套处理技术是处理流程设计中不可缺少的一部分,本文通过MPFI数据规则化技术,采用防漏频傅里叶变换重构有效信号,有效预测了近偏移距缺失道,弱化了采集脚印,提高了数据信噪比;利用DWD+GSMP组合多次波压制技术,较好地压制了浅水环境水体多次波和自由界面多次波。

[1]沈财余.面向地质目标的地震采集设计优化方法[J].石油地球物理勘探,2008,43(5):493-507

SHEN C Y.Optimum design method of geologic target-oriented seismic acquisition[J].Oil Geophysical Prospecting,2008,43(5):493-507

[2]李绪宣,王建花,张金淼,等.南海深水区地震资料采集设计和处理关键技术及其野外试验效果[J].中国海上油气,2013,25(6):8-14

LI X X,WANG J H,ZHANG J M,et al.Some seismic acquisition designs and key processing techniques and their testing effects in the deep water areas,South China Sea[J].China Offshore Oil and Gas,2013,25(6):8-14

[3]刘财,刘洋,王典,等.一种频域吸收衰减补偿方法[J].石油物探,2005,44(2):116-118

LIU C,LIU Y,WANG D,et al.A method to compensate strata absorption and attenuation in frequency domain[J].Geophysical Prospecting for Petroleum,2005,44(2):116-118

[4]杨学亭,刘财,刘洋,等.基于连续小波变换的时频域地震波能量衰减补偿[J].石油物探,2014,53(5):523-529

YANG X T,LIU C,LIU Y,et al.The attenuation compensation of seismic wave energy in time-frequency domain based on the continuous wavelet transform[J].Geophysical Prospecting for Petroleum,2014,53(5):523-529

[5]KIM H S,BELL M L.3-D turning-ray tomography and its application to Mississippi Delta[J].Expanded Abstracts of 70thAnnual Internat SEG Mtg,2000:593-596

[6]MITTET R.A simple design procedure for depth extrapolation operators that compensate for absorption and dispersion[J].Geophysics,2007,72(2):S105-S112

[7]何汉漪,朱宏彰,张树林.海上多波地震在莺歌海盆地的应用[J].中国海上油气(地质),1999,13(5):310-315

HE H Y,ZHU H Z,ZHANG S L.Application of offshore multi-component seismic in Yinggehai basin[J].China Offshore Oil and Gas(Geology),1999,13(5):310-315

[8]傅旦丹,刘一峰,沈丽丽,等.南海某构造二维多分量地震资料处理[J].石油物探,2004,43(2):181-184

FU D D,LIU Y F,SHEN L L,et al.Processing the 2D multicomponent seismic data of a structure in South China Sea[J].Geophysical Prospecting for Petroleum,2004,43(2):181-184

[9]杨锡锟,沈丽丽.莺歌海盆地多波地震勘探资料的采集与处理[J].中国海上油气(地质),1999,13(5):316-323

YANG X K,SHEN L L.The data acquisition and processing techniques of the multiwave seismic exploration in the Yinggehai basin[J].China Offshore Oil and Gas(Geology),1999,13(5):316-323

[10]黄中玉.多分量地震勘探的机遇和挑战[J].石油物探,2001,40(2):131-137

HUANG Z Y.The opportunity and challenge for the multi-component seismic survey[J].Geophysical Prospecting for Petroleum,2001,40(2):131-137

[11]田立新,夏庆龙,刘春成,等.多分量地震勘探在渤海海域的应用[J].石油勘探与开发,2005,32(3):78-83

TIAN L X,XIA Q L,LIU C C,et al.Application of multiple component seismic exploration in Penglai 19-3 Oilfield,Bohai Bay[J].Petroleum Exploration and Development,2005,32(3):78-83

[12]张树林,李绪宣,何汉漪,等.近海四分量地震构造勘探效果分析[J].石油地球物理勘探,2002,37(5):446-454

ZHANG S L,LI X X,HE H Y,et al.Analysis on prospecting effects of offshore 4-C seismic structure[J].Oil Geophysical Prospecting,2002,37(5):446-454

[13]李彦鹏,孙鹏远,魏庚雨,等.利用陆上三分量数据改善气云区构造成像[J].石油地球物理勘探,2009,44(4):417-424

LI Y P,SUN P Y,WEI G Y,et al.Using 3-C land seismic data to improve structural imaging in gas chimney area[J].Oil Geophysical Prospecting,2009,44(4):417-424

[14]TONN R.The determination of the seismic quality factorQfrom VSP data:a comparison of different computational methods[J].Geophysical Prospecting,1991,39(1):1-27

[15]YANG S J.GAO J H,CHEN W C,et al.Comparisons of four methods used for seismic quality factors estimation[J].Expanded Abstracts of 79thAnnual Internat SEG Mtg,2009:2472-2476

[16]TRAYNIN P,LIU J,REILLY J M.Amplitude and bandwidth recovery beneath gas zones using Kirchhoff prestack depth Q-migration[J].Expanded Abstracts of 78thAnnual Internat SEG Mtg,2008:2412-2416

[17]BELL M L,LARA R,GRAY W C.Application of turning-ray tomography to the offshore Mississippi delta[J].Expanded Abstracts of 64thAnnual Internat SEG Mtg,1994:1509-1512

[18]BISHOP T N,BUBE K P,CUTLER R T,et al.Tomographic determination of velocity and depth in laterally varying media[J].Geophysics,1985,50(6):903-923

[19]王兴芝,钟明睿,杨薇,等.海上拖缆三维地震数据规则化技术研究[J].地质科技情报,2015,34(2):207-212

WANG X Z,ZHONH M R,YANG W,et al.The regularization of 3D seismic data from marine towed streamer[J].Geological Science and Technology Information,2015,34(2):207-212

[20]王兴芝,李添才,肖二莲,等.基于反漏频傅里叶变换的数据规则化技术在海上三维拖缆地震资料处理中的应用[J].中国海上油气,2014,26(4):25-28

WANG X Z,LI T C,XIAO E L,et al.An application of data regularization based on anti-leakage Fourier transform to 3D marine seismic data processing by tow-cable[J].China Offshore Oil and Gas,2014,26(4):25-28

[21]曹国滨.频率-空间域数据规则化压制采集脚印技术研究[J].石油物探,2010,49(4):359-363

CAO G B.Study of suppressing acquisition footprint by data regularization in frequency-space domain[J].Geophysical Prospecting for Petroleum,2010,49(4):359-363

[22]贾友珠,董伟,张印堂.三维地震资料连片处理中的面元均化技术——以桩海地区为例[J].油气地质与采收率,2002,9(2):50-52

JIA Y Z,DONG W,ZHANG Y T.Bin averaging technique in merging processing of 3D seismic data-taking Zhuanghai area as an example[J].Oil & Gas Recovery Efficiency,2002,9(2):50-52

[23]XU S,ZHANG Y,PHAM D L,et al.Antileakage Fourier transform for seismic data regularization[J].Geophysics,2005,70(4):V87-V95

[24]渥·伊尔马滋.地震资料分析:地震资料处理、反演和解释[M].北京:石油工业出版社,2006:21-25

YILMAZ O.Seismic data analysis:seismic data processing,inversion and interpretation[M].Beijing:Petroleum industry press,2006:21-25

[25]孟小红,郭良辉,张致付,等.基于非均匀快速傅里叶变换的最小二乘反演地震数据重建三维地震资料连片处理中的面元均化技术[J]．地球物理学报,2008,51(1):235-241

MENG X H,GUO L H,ZHANG Z F,et al.Reconstruction of seismic data with least squares inversion based on nonuniform fast Fourier transform[J].Chinese Journal of Geophysics,2008,51(1):235-241

[26]范宝仓,谢涛.部分去除表面多次波技术在浅水海域的应用[J].石油天然气学报,2012,34(12):73-75

FAN B C,XIE T.The application of partial SRME in shallow sea[J].Journal of Oil and Gas Technology,2012,34(12):73-75

[27]刘俊,吴淑玉,高金耀,等.南黄海中部浅水区多次波衰减技术及其效果分析[J].物探与化探,2016,40(3):568-577

LIU J,WU S Y,GAO J Y,et al.An effectiveness analysis of multiple depression technique in the Shallow water of the central uplift in South Yellow Sea[J].Geophysical and Geochemical Exploration,2016,40(3):568-577

[28]张志军,魏天罡.浅水多次波的联合衰减技术在渤海海域LD地区的应用[J].中国石油勘探,2013,18(1):59-65

ZHANG Z J,WEI T G.Application of joint shallow-water multiple attenuation technique in LD zone of Bohai Sea[J].China Petroleum Exploration,2013,18(1):59-65

[29]张振波,轩义华,刘宾.基于各向异性理论的深水区地震资料叠前处理技术[J].吉林大学学报(地球科学版),2014,44(3):1031-1038

ZHANG Z B,XUAN Y H,LIU B.The pre-stack seismic data processing technology series based on anisotropic theory and its application in South China Sea deepwater Baiyun 6-1 structure[J].Journal of Jilin University(Earth Science Edition),2014,44(3):1031-1038

[30]张亚斌,施荣富,姚刚.Q-Marine技术和特色处理技术在东海海域油气区的应用[J].石油天然气学报,2013,35(6):47-52

ZHANG Y B,SHI R F,YAO G.Application of Q-Marine seismic exploration technology and specialized processing technology in Donghai oil and gas area[J].Journal of Oil and Gas Technology,2013,35(6):47-52