基于鲁棒自适应最小方差信号无畸变响应波束形成的高密度数据室内组合方法研究

蓝 阳,谢俊法,杨志鹏,崔保生

(1.中国石油天然气股份有限公司新疆油田分公司勘探开发研究院,新疆克拉玛依834000;2.中国石油天然气股份有限公司勘探开发研究院西北分院,甘肃兰州730020)

随着“两宽一高”地震勘探的推进,单点高密度地震采集技术越来越受到重视。刘欣欣等[1]详细论述了国内外单点高密度地震勘探技术的发展。钱荣钧[2]指出,高密度、均匀空间采样技术既是解决当前地震勘探所面对的复杂问题的技术对策,也是今后地震勘探的发展趋势。李燕燕等[3]指出,全方位高密度、超多道数接收仍是不变的地震采集技术发展方向,预期2020年后形成全方位高密度技术、超多道数固定排列采集技术。曲寿利[4]结合国内外应用实例分析了高密度采集技术的关键问题,指出高密度三维地震技术是老油区二次勘探的关键技术之一。正是由于单点高密度采集技术的重要作用,多年来国内学者对单点高密度数据特点进行了详细分析。曹务祥[5]、杨贵祥等[6]、杨照海等[7]、CAI等[8]、冯刚等[9]研究了数字检波器单点采集地震资料的特点,并结合实际数据分析了不同组合间距对地震信号的影响,同时指出接收道组合对资料前期的预处理十分必要。于世焕等[10]对比分析了数字检波器单点采集与组合接收野外试验资料,发现单点采集有利于保护有效波频率成分,但会降低地震数据的信噪比。张永刚等[11]、胡莲莲等[12]结合数值模拟和实测数据,从分辨率、信噪比和数据处理的角度剖析了单点高密度地震数据与常规地震数据的区别。近年来单点高密度采集技术在川东碳酸盐岩[13]、塔里木盆地缝洞型碳酸盐岩[14]、辽河青龙台[15]等地区发挥了重要作用。

目前国内针对单点高密度地震数据的室内组合,主要利用数据一致性处理技术消除组合数据间的差异,然后采用算术平均法进行组合叠加。曹务祥[5]提出对地震资料进行严格的去相干噪声、道间振幅和相位的差异调整以及道间时差的校正等处理后,再进行接收道组合的方法。孙成禹等[16]利用自适应匹配

滤波一致性校正方法对组合数据进行一致性校正,改善了单点高密度地震数据室内组合的效果。刘志鹏等[17]利用海上高密度地震数据对比分析了相位谱校正组合与自适应匹配组合方法。在对地震数据进行去噪、静校正、一致性校正等预处理后,虽然尽可能地消除了组合内波形间的差异,但仍然无法达到完全一致,此时采用算术平均法进行组合,会导致高频信号的损失。

本文通过引入波束形成理论,分析了基于窗函数的波束形成算法和基于PCA的鲁棒自适应MVDR波束形成算法,并将其应用于点高密度地震数据的室内组合中,拓宽了室内组合方法的研究思路,并且利用Marmousi2模型正演模拟数据对本文组合方法进行了应用测试。

1 方法原理

对单点高密度采集的地震数据,组合地震道经过一致性处理后,采用算术平均法合成数据,如图1所示。地震道的组合可以视为地震信号经过波束形成器的过程,其输出信号为[18-20]:

(1)

由公式(1)可以发现,地震道的组合受权重矢量的影响。在数据处理中,各种窗函数的应用十分广泛,将窗函数引入权重矢量可以拓宽理论研究组合的思路。目前常用的窗函数有均匀窗、余弦窗、汉明窗和布莱克曼窗[21]等,分别如下:

图1 常规波束形成器示意

1) 均匀加权。

(2)

其中,n=1,2,…,N。

2) 余弦加权。

(3)

3) 汉明加权。

(4)

4) 布莱克曼加权。

(5)

分别将公式(2)～公式(5)代入公式(1)中,可得到不同窗函数加权下的输出信号。

自适应MVDR波束形成算法能在有效保障信号保真度的前提下提高地震数据的信噪比,建立二者之间的平衡关系,但需已知较为准确的标准道。考虑到单点高密度地震数据一致性校正处理后,已尽可能地消除了组合内波形间的差异,而PCA可以得到组合数据的主成分,其第一主成分反映了有效信号的信息,因此,可以利用第一主成分重构地震数据,以此作为组合数据的标准道。针对单点高密度地震数据的特点,本文引入自适应MVDR波束形成算法,进一步利用PCA和基于线性组合的对角载入算法对其进行改进,构建鲁棒自适应MVDR加权函数。

1.1 MVDR波束形成器

利用公式(1)构建自适应MVDR滤波器,即:在保证信号不失真的情况下,使得信噪比最大[22-23]。

(6)

其中,Rs=E{S(t)S(t)T}和Ri=E{I(t)I(t)T}分别是有效信号的协方差矩阵和噪声的协方差矩阵。

利用拉格朗日乘数法构建目标函数:

(7)

其中,λ为拉格朗日乘数。利用(7)式对WT求偏导并令其为零,可得到:

(8)

对公式(8)进行简单变换,即:

(9)

(10)

在实际应用中,噪声的协方差矩阵Ri难以获取,当利用实际数据的协方差矩阵Rx=E{X(t)X(t)T}=Ri+Rs代替Ri时,(10)式的最优解不会发生改变,此时:

(11)

由公式(11)可以发现,要求得权重矢量Wmax需已知有效信号的协方差矩阵Rs,因此,需要利用PCA提取组合信号的有效信息[24-25]。具体实现步骤如下:

1) 设XN×M表示N道地震数据,每道地震数据的点数为M,对XN×M进行归一化(零均值或者去均值处理),得到δXN×M;

2) 计算δXN×M的协方差矩阵DN×M;

3) 对协方差矩阵进行特征值分解,即D·V=Λ·V。其中,Λ为特征值矩阵,Λ=diag(λ1,λ2,…,λN),且λ1≥λ2≥…≥λN≥0;V是N×N的特征向量矩阵,且VVT=1;

4) 计算δXN×M在V中的投影,即Y=VTδX,得到δXN×M矩阵的主成分y1,y2,…,yN;

1.2 基于对角载入的鲁棒MVDR自适应波束形成算法

在实际应用中,准确的协方差矩阵Rx和Rs无法获得,但可以估计得到,因此有效信号协方差矩阵和信号协方差矩阵可以近似为[22]:

式中:M表示每道数据的采样点数。

此时,公式(11)变为:

(14)

(15)

(16)

式中:‖·‖表征矢量范数。

将公式(15)代入公式(16)中并分别对α,β求偏导,经过推导得到最优解α0,β0[26]:

(17)

(18)

(19)

此时,公式(14)变为:

(20)

将公式(20)代入公式(1)中,可以得到经过鲁棒自适应MVDR波束形成器后的输出信号。

2 数值分析

为了进一步对比分析不同加权组合情况下地震记录的特点,分别从组合后输出数据的信噪比、幅频特性、同相轴连续性等几个方面进行分析比较。

根据地震波传播的特征,选取雷克子波作为有效信号。雷克子波的表达式如下[27]:

(21)

式中:t表示时间;fp为子波主频。

为了分析不同加权组合对地震信号信噪比的影响,利用公式(21)设计10道主频fp=50Hz的地震信号S(t),每道地震信号加入随机噪声I(t),则第n道输入信号的信噪比为[21]:

(22)

当采用加权组合合成地震信号时,组合后输出信号的信噪比满足[22]:

(23)

图2为10道含随机噪声且不含时差的组合数据,表1为不同加权接收对该组合数据信噪比的改善情况。图3为10道含随机噪声和时差的组合数据,其中第2道和第8道有较小的时差,表2为不同加权接收对该组合数据信噪比的改善情况。可以发现加权组合可以提高地震数据的信噪比:①在信号不存在时差的情况下,各加权算法对信噪比的改善差异较小;②当信号中存在时差时,相比于其它几种加权组合,基于PCA的鲁棒自适应MVDR加权算法能进一步提高信号的信噪比。

图2 不含时差的组合数据

图3 含时差的组合数据

表1 不同加权组合对不含时差的组合数据信噪比改善情况

表2 不同加权组合对含时差的组合数据信噪比改善情况

为了分析不同加权对地震数据幅频特性及地震记录同相轴连续性的影响,在不考虑噪声的情况下,构建含气的Marmousi2模型,如图4所示。模型大小为1001×451,网格大小为3m×3m,时间采样间隔为1ms,记录时间为2s,激发主频为30Hz,炮点在第470道距地面4个采样点(12m)处。对模型进行正演模拟,得到炮记录如图5所示。分别对炮记录进行均匀加权组合、余弦加权组合、汉明加权组合、布莱克曼加权组合和基于PCA的鲁棒自适应MVDR加权组合并显示在同一坐标系下,其正演记录如图6所示。

图4 Marmousi2含气模型

由图6可知,相比于常规的均匀加权组合,基于其它窗函数的加权组合和基于PCA的鲁棒自适应MVDR加权组合均能使炮记录上同相轴的能量更加均衡,更易识别追踪,其中基于PCA的鲁棒自适应MVDR加权组合效果更为明显。同时,抽取第90道地震记录进行频谱分析(图7)发现,相比于均匀加权,其它几种加权算法均能更好地保护中高频成分,取得较好的保真效果。

图5 炮记录

图6 不同加权组合的正演记录a 均匀加权; b 余弦加权; c 汉明加权; d 布莱克曼加权; e 鲁棒自适应MVDR加权

图7 第90道的幅频曲线

为了进一步对比分析各加权组合对信号的影响和它们之间的差异,对图6加权组合正演的地震记录进行F-K变换并计算各加权地震记录F-K变换结果与均匀加权变换结果的差值(图8～图12)。由图8 ～图12可以发现:与常规均匀加权相比,基于窗函数的加权组合能保护部分高频信号;基于PCA的鲁棒自适应MVDR加权组合对地震记录各频率均有一定的保护作用,图12a中20～40Hz频段范围内信号的保真效果更为明显,而进行正演模拟时采用的子波主频为30Hz,体现了MVDR波束形成器的优势,公式(6)中WTRsW=1保证了有效信号在不失真的情况下进行组合,同时,对比差值剖面可以发现,基于PCA的鲁棒自适应MVDR加权组合对低波数信号的保护更加明显。由于地层速度随着深度的增加呈增大趋势,因此在频率一定的情况下,深层信号的波数越小。可见,相比于其它组合方式,基于PCA的鲁棒自适应MVDR加权组合对有效信号乃至深层信号能起到更为显著的保护作用。

图8 均匀加权地震记录F-K变换结果

图9 余弦加权地震记录F-K变换结果(a)及其与图8的差值剖面(b)

图10 汉明加权地震记录F-K变换结果(a)及其与图8的差值剖面(b)

图11 布莱克曼加权地震记录F-K变换结果(a)及其与图8的差值剖面(b)

3 结论与认识

本文重点研究了基于窗函数的波束形成算法和基于PCA的鲁棒自适应MVDR波束形成算法,将其与目前常用的均匀加权(算术平均)算法对比,得到以下认识:①组合能提高地震数据的信噪比,余弦加权、汉明加权、布莱克曼加权能进一步保护地震数据的高频成分;②基于PCA的鲁棒自适应MVDR加权组合实现了信号保真与信噪比的有效统一,对深层信号的保护作用较为明显。

针对单点高密度地震数据的特点,对地震数据的去噪、静校正和一致性校正等预处理很有必要,基于波束形成理论的数字组合方法能够对当前单点高密度地震数据的室内组合处理技术进行适当的补充,下一步将研究该方法在单点高密度实际资料处理中的应用。

[1] 刘欣欣,吴国忱,梁锴.单点高密度地震勘探技术研究综述[J].地球物理学进展,2009,24(4):1354-1366

LIU X X,WU G C,LIANG K.The review of point-source/point-receiver high density seismic exploration technology[J].Progress in Geophysics,2009,24(4):1354-1366

[2] 钱荣钧.关于地震采集空间采样密度和均匀性分析[J].石油地球物理勘探,2007,42(2):235-243

QIAN R J.Analysis on spatial sampling density and uniformity of seismic acquisition[J].Oil Geophysical Prospecting,2007,42(2):235-243

[3] 李燕燕,赵殿栋,于世焕,等.中石化陆上地震采集技术现状与发展趋势[J].石油物探,2013,52(4):363-371

LI Y Y,ZHAO D D,YU S H,et al.Status and trend on land seismic acquisition technique of SINOPEC[J].Geophysical Prospecting for Petroleum,2013,52(4):363-371

[4] 曲寿利.高密度三维地震技术——老油区二次勘探的关键技术之一[J].石油物探,2006,45(6):557-562

QU S L.High-density 3-D seismic technique:one of the key technique of the secondary exploration in old oil region[J].Geophysical Prospecting for Petroleum,2006,45(6):557-562

[5] 曹务祥.单道接收地震资料的室内组合方法[J].石油地球物理勘探,2006,41(6):615-618

CAO W X.Indoor array for unite single trace receiving seismic data[J].Oil Geophysical Prospecting,2006,41(6):615-618

[6] 杨贵祥,杨振升,仲伯军.单点单分量高密度地震采集技术及应用[J].油气藏评价与开发,2011,1(3):12-18

YANG G X,YANG Z S,ZHONG B J.Single-point and single-component high density seismic acquisition technology and its application[J].Reservoir Evaluation and Development,2011,1(3):12-18

[7] 杨照海,凡正才,赵前华.单点数字检波器道组合技术[J].石油地球物理勘探,2008,43(增刊2):88-93

YANG Z H,FAN Z C,ZHAO Q H.Overlapping technique of single-point digital geophone[J].Oil Geophysical Prospecting,2008,43(S2):88-93

[8] CAI X L,LIU X W,DENG C Y,et al.Characteristics analysis on high density spatial sampling seismic data[J].Applied Geophysics,2006,3(1):48-54

[9] 冯刚,毕丽飞,李建明,等.单点数字检波器地震资料特点及处理对策[J].石油地球物理勘探,2008,43(增刊2):115-120

FENG G,BI L F,LI J M,et al.Seismic data character of single-point digital geophone and processing strategy[J].Oil Geophysical Prospecting,2008,43(S2):115-120

[10] 于世焕,赵殿栋,于晨.数字检波器单点地震采集与组合接收对比试验[J].石油物探,2012,51(3):280-284

YU S H,ZHAO D D,YU C.Comparative experiment on digital geophone single point and array receiving in the seismic acquisition[J].Geophysical Prospecting for Petroleum,2012,51(3):280-284

[11] 张永刚,王赟,尹军杰.单点高密度地震数据处理分析与初步评价[J].石油地球物理勘探,2010,45(2):201-207

ZHANG Y G,WANG Y,YIN J J.Single point high density seismic data processing analysis and initial evaluation[J].Oil Geophysical Prospecting,2010,45(2):201-207

[12] 胡莲莲,叶树刚,芦俊,等.关于单点高密度地震数据的室内组合分析[J].地球物理学进展,2010,25(6):2115-2124

HU L L,YE S G,LU J,et al.Analysis of indoor combination of Q-Land seismic data[J].Progress in Geophysics,2010,25(6):2115-2124

[13] 黄锐.川东北碳酸盐岩地区地震勘探技术难点与对策[J].石油物探,2008,47(5):476-482

HUANG R.Technical difficulties in and countermeasures for seismic exploration in the carbonate area of northeast Sichuan basin[J].Geophysical Prospecting for Petroleum,2008,47(5):476-482

[14] 郭念民,陈猛,崔永福,等.碳酸盐岩储层单点高密度采集三维地震勘探实例[J].石油物探,2016,55(6):771-780

GUO N M,CHEN M,CUI Y F,et al.The application of single-point high-density seismic acquisition for carbonate reservoir 3D seismic exploration[J].Geophysical Prospecting for Petroleum,2016,55(6):771-780

[15] 李明生,李雅楠,董文波.辽河青龙台地区高密度全数字三维地震采集技术及效果[J].中国石油勘探,2017,22(1):106-112

LI M S,LI Y N,DONG W B.High-density digital 3D seismic acquisition technology and its application results in Qinglongtai area,Liaohe Oilfield[J].China Petroleum Exploration,2017,22(1):106-112

[16] 孙成禹,孙甜甜,彭洪超,等.单点地震数据的一致性组合叠加方法[J].石油地球物理勘探,2009,44(S1):23-28

SUN C Y,SUN T T,PENG H C,et al.Consistent array stacking method for single point seismic data[J].Oil Geophysical Prospecting,2009,44(S1):23-28

[17] 刘志鹏,赵伟,张云鹏.海上高密度地震数据数字组合方法研究与应用[J].中国海上油气,2013,25(5):21-25

LIU Z P,ZHAO W,ZHANG Y P.A research and application of the digital grouping methods for offshore high-density seismic data[J].China Offshore Oil and Gas,2013,25(5):21-25

[18] 杨益新.声呐波束形成与波束域高分辨方位估计技术研究[D].陕西西安:西北工业大学,2002

YANG Y X.Studies on beamforming and beamspace high resolution bearing estimation techniques in sonar systems[D].Xi’an:Northwestern Polytechnical University,2002

[19] 蓝阳,孙成禹,邵婕.基于波束形成理论的组合检波技术研究[J].石油物探,2015,54(5):493-500

LAN Y,SUN C Y,SHAO J.The geophone array design based on beam-forming[J].Geophysical Prospecting for Petroleum,2015,54(5):493-500

[20] 郭伟.水下监测系统中目标探测若干关键技术研究[D].湖南长沙:国防科学技术大学,2011

GUO W.Research on key methods of target detecting for underwater monitoring system[D].Changsha:National University of Defense Technology,2011

[21] VAN TREES H L.Optimum array processing,part Ⅳ of detection estimation,and modulation[M].New York:John Wiley & Sons,2002:1-1443

[22] PANEA I,DRIJKONINGEN G.The spatial data-adaptive minimum-variance distortionless-response beamformer on seismic single-sensor data[J].Geophysics,2008,73(5):Q29-Q42

[23] SHAHBAZPANAHI S,GERSHMAN A B,LUO Z Q,et al.Robust adaptive beamforming for general-rank signal models[J].IEEE Transactions on Signal Processing,2003,51(9):2257-2269

[24] 武莹,陆从德,杜兴忠,等.主成分分析在航空瞬变电磁去噪中的应用[J].物探化探计算技术,2014,36(2):170-176

WU Y,LU C D,DU X Z,et al.A denoising method based on principal component analysis for airborne transient electromagnetic data[J].Computing Techniques for Geophysical and Geochemical Exploration,2014,36(2):170-176

[25] 陈斌,陆从德,刘光鼎.基于核主成分分析的时间域航空电磁去噪方法[J].地球物理学报,2014,57(1):295-302

CHEN B,LU C D,LIU G D.A denoising method based on kernel principal component analysis for airborne time domain electromagnetic data[J].Chinese Journal of Geophysics,2014,57(1):295-302

[26] STOICA P,LI J,ZHU X,et al.On using a priori knowledge in space-time adaptive processing[J].IEEE Transactions on Signal Processing,2008,56(6):2598-2602

[27] 孙成禹,李振春.地震波动力学基础[M].北京:石油工业出版社,2011:1-197

SUN C Y,LI Z C.The base of seismic wave dynamics[M].Beijing:Petroleum Industry Press:1-197