南海浅水区域三维τ -p域鬼波压制分析❋

2020-06-23谢玉洪张迎朝欧阳敏方中于李三福张健男

中国海洋大学学报（自然科学版） 2020年7期

谢玉洪，张迎朝,李列,欧阳敏,谭军，方中于,李三福,张健男

(1.中国海洋石油集团有限公司,北京 100010;2.中国海洋石油湛江分公司,广东湛江 524057;3.中国海洋大学海洋地球科学学院,山东青岛 266100;4.中国海洋大学海底科学与探测技术教育部重点实验室,山东青岛 266100;5.中海油服物探事业部特普公司,广东湛江 524057)

三维地震勘探工区位于珠三拗陷中部文昌A凹陷，主要发育文昌组、恩平组两套四种类型的烃源岩[1]，并具有中新统珠江组海相砂岩储集层与海侵泥岩盖层[2]，从而形成良好的生、储、盖成藏组合。该区域还发育一系列背斜、断背斜及断块构造等圈闭[1]，构成了“下生上储，陆生海储”的古近-新近系含油气系统及油气聚集区。文昌A凹陷低压构造脊带呈北东向展布，形成了多个局部构造，其中存在低幅度小背斜和异常体。该区三维地震资料分别为三个年度采集，部分年度的数据质量差、信噪比低，导致常规地震处理结果为圈闭落实、有利目标评价带来极大困难。

在海洋地震勘探中，需将震源和检波器放置于海平面以下一定深度。气枪震源激发的能量其中一部分直接向上传播到海面，然后发生反射并向地下深部传播，遇到地下界面再反射回来被电缆检波器所接收，形成尾随于一次反射波之后的震源鬼波；与此类似，震源激发后能量向下传播，遇到地下界面后反射回海面，再发生下行反射后被电缆检波器所接收，同样会形成位于一次反射波之后的电缆鬼波[3]。由于鬼波的波形、频率、视速度等均与一次反射波接近，其不但导致地震资料出现陷频，还造成地震资料出现多轴、多相位现象，从而严重影响局部构造解释——尤其会误导针对低幅度小背斜、异常体的识别分析。

为了实现鬼波的压制处理，进而获得具有较宽频带的地震资料，国内外的地球物理学家提出了联合反褶积鬼波压制[4-7]、基于格林定理的鬼波压制[8-13]以及基于逆散射级数法的鬼波压制[14-17]等多类方法。上述各类方法中，联合反褶积具有明显的效率优势，当虚反射与一次波时差较小(譬如海上常规拖缆地震)时，若能够求出较为准确的虚反射时间延迟量，可获得理想的鬼波压制效果。在频率——空间域联合反褶积鬼波压制方法[18-19]的基础上，Wang等[20-21]提出了在F-K域、τ -p域进行的联合反褶积鬼波压制算法，将入射角引入到鬼波延迟时估算而减小了误差，从而在实际地震处理中获得了更优的效果。

在浅水区域文昌A凹陷的三维地震资料鬼波压制中，由于边缆的纵向偏移距较大，相应单缆记录的二维τ -p变换结果出现了严重的空间假频，致使基于二维τ -p变换的鬼波剔除遇到了困难。为此，本文实现了基于三维τ -p变换的鬼波压制技术，然后对工区三个年度的三维地震资料进行了鬼波压制处理，其最终偏移数据体中紧随一次波信号的鬼波同相轴旁瓣得到明显压制，地层反射更加清晰、信噪比显著提高，从而验证了三维τ -p域鬼波压制技术的有效性。

1 三维τ -p域鬼波压制的基本原理

τ -p域鬼波压制方法首先将偏移距域地震数据变换到τ -p域，即实现地震数据的平面波分解，然后对每个平面波分量进行鬼波压制。若已知震源(或检波器)沉放深度与海水的地震波速，鬼波延迟时可视为波场入射(或出射)角度的函数。因此，在τ -p域可根据射线参数(水平慢度)p估计出具有不同入射(或出射)角度平面波分量的鬼波延迟时，据此实现震源(或检波器)鬼波的精准压制。

1.1 三维τ -p变换的基本原理

三维离散空间——时间域的τ -p正、反变换的计算公式[22]分别为

(1)

式中：d(x,y,t)为偏移距域数据；u(px,py,τ)为τ-p变换的结果；d′(x,y,t)为τ-p反变换的地震记录；t为双程旅行时，τ为双程零偏移距旅行时；x、y分别为x、y方向的偏移距；px、py分别为x、y方向的射线参数(水平慢度)。

为了保证τ-p变换过程的稳定性，令输入记录d(x,y,t)与反变换结果d′(x,y,t)的误差平方和最小，从而导出频率域的最小平方三维τ-p正、反变换的表达式[23-24]：

(2)

式中:d、u分别表示偏移距域、(2)τ-p域的单频数据矩阵；d′为τ-p反变换的结果；β为常数，I表示单位矩阵；上标T表示矩阵转置；*表示复共轭。

根据偏移距信息(x,y)与射线参数(px,py)确定复矩阵L，即

(3)

在公式(2)中，复矩阵LT*为L的共轭转置矩阵，因此有

LT*=

(4)

对于三维τ -p变换而言，由于射线参数P(px,py)的规模极为庞大，而输入记录与待求解的τ -p域数据中均存在大量零值，因此利用公式中第1个方程求取u的过程属于大规模稀疏矩阵的线性方程组求解问题，通常采用共轭梯度法进行迭代计算。

1.2 基于三维τ -p变换的鬼波压制

在假设地震波的传播满足平面波的前提条件下，鬼波延迟时Δt(鬼波与一次波的旅行时差)仅与震源(或拖缆检波器)沉放深度、地震波入射(或出射)角度有关。在图1所示一次波与鬼波的传播过程示意图中，设震源的沉放深度为h，平面波的入射(或出射)角度为θ，则鬼波的延迟时Δt为

(5)

式中：h为震源(或拖缆检波器)深度；v为地震波在水中传播的速度。

图1 平面波假设条件下的震源鬼波传播过程示意图Fig.1 Schematic diagram of the propagation of source ghost under the assumption of plane wave

在τ -p域中，射线参数p与地震波入射(或出射)角度θ存在以下关系

(6)

将式(6)代入到式(5)中，可得

(7)

在三维情况下，射线参数p可表示为

(8)

因此，三维情况下的鬼波延迟时Δt为

(9)

根据波场延拓理论，鬼波可视为一次波波场在海面发生下行反射所形成的下行波场，即

g(px,py,ω)=R(px,py,ω)u′(px,py,ω)e-iωΔt。

(10)

式中，g表示下行的鬼波波场；u′表示一次波波场；R为鬼波入射到平界面的反射系数，可基于Zoeppritz方程组予以求解，由于海面为自由界面，R通常取-1。因此，总波场u为

u(px,py,ω)=u′(px,py,ω)+g(px,py,ω)=

u′(px,py,ω)+

R(px,py,ω)u′(px,py,ω)e-iωΔt。

(11)

利用式可求出一次波波场，即消除鬼波后的结果

(12)

为简便起见，通常利用鬼波算子G来描述形成鬼波的波场延拓过程，即

G(px,py,ω)=1+R(px,py,ω)e-iωΔt。

(13)

因此，式(12)所描述的鬼波压制过程可简单地表示为

(14)

其中ε为保证分母不为0的阻尼因子。

图2展示了鬼波压制的详细流程，主要包括以下步骤：

①将包含多缆数据的三维炮集记录由时间域变换到频率域，针对每个频率成分ω解方程u=(LT*L+βI)·LT*d(见公式(2))，即可获得频率域的τ -p域记录。

图2 τ -p域鬼波压制流程图Fig.2 The processing flow chart of τ -p domain deghosting

2 文昌A凹陷区域三维鬼波压制分析

目标工区具有三个时期采集的地震数据，其采集时间分别为1997、2006和2014年，且野外观测参数具有明显差异(见表1)。为了获得较优的成像效果，需要基于工区内的三批数据进行地震联片处理。

表1 与鬼波压制有关的野外采集参数Table 1 Acquisition parameters are related of deghosting

Note:①Year；②Cable number；③Trace number/Cable；④Cable interval；⑤Vertical minimum offset；⑥Vertical trace interval；⑦Source depth；⑧Receiver depth

输入衰减多次波后的炮记录，基于图2所示的处理流程进行鬼波压制。为保证三维τ -p变换的精度，需采用较大的px、py值范围与较小的px、py值间隔，通常根据主要反射同相轴在最大偏移距道中的时移量γ来描述τ -p变换参数：

其中，γx、γy分别为反射同相轴在X、Y方向最大偏移距道中的最大时移量，据此可求出射线参数p。

(15)

因此，根据表1中的野外观测参数可计算出px、py的范围与间隔。此外，在基于公式(14)进行鬼波压制时，为保证公式中分母不为0，设阻尼因子ε为0.01。

采用上述参数进行三维τ -p域的鬼波压制处理，图2、3分别展示了鬼波压制前、后的一组炮集记录对比。为清晰起见，对图3、4中所示数据的局部记录(道号范围：1～200，时间范围：0～2 s)进行放大显示，分别见图5、6。通过对比可知，输入炮记录中紧随一次波信号的鬼波旁瓣得到明显压制，地震数据的信噪比显著提高。抽取同震源炮集的第11道与第101道组成共检波器道集，压制鬼波前、后的共检波器道集分别如图7、8所示，因震源鬼波与检波器鬼波叠加所形成的同相轴旁瓣已被去除，两剖面中的强反射同相轴处于“反相”状态，并且由原来的多相位变成单相位。创建图7、8所示道集的自相关谱，所得结果分别如图9、10所示，其中箭头指向的一次波与鬼波旁瓣相关所形成的能量已被完全消除。

(2006年采集，1～3缆。Acquisition in 2006,1～3 Cable.)图3 鬼波压制前的炮集记录示例Fig.3 Shot record before deghosting

(2006年采集，1～3缆。Acquisition in 2006,1～3 Cable.)图4 鬼波压制后的炮集记录示例Fig.4 Shot record after deghosting

(道号范围：1～200，时窗范围：0～2 s。Trace number:1～200,Time window:0～2 s.)
图5 图3所示记录的局部放大显示
Fig.5 Enlarged part of the record as shown in Fig.3

(道号范围：1～200，时窗范围：0～2 s。Trace number:1～200,Time window:0～2 s.)
图6 图4所示记录的局部放大显示
Fig.6 Enlarged part of the record as shown in Fig.4

图7 压制鬼波前的共检波器道集Fig.7 Common receiver gathers before deghosting

图8 压制鬼波后的共检波器道集Fig.8 Common receiver gathers after deghosting

图9 压制鬼波前共检波器道集的自相关谱Fig.9 Autocorrelation spectrum of common receiver gathers before deghosting

图10 压制鬼波后共检波器道集的自相关谱Fig.10 Autocorrelation spectrum of common receiver gathers after deghosting

图11、12展示了压制鬼波前、后的叠前时间偏移剖面。在原始成像剖面中，时差较小的鬼波信号尾随于一次波之后，且叠加在一起形成了同相轴旁瓣，导致剖面中的反射层位处于多相位状态，这对后续的解释分析是非常不利的。在图12所示去除鬼波的偏移剖面中，紧随一次波信号的鬼波旁瓣得到明显压制，地层反射更加清晰，中深层低频能量增强，信噪比显著提高。