专利名称：一种基于波场分离原理提高地震资料信噪比的去噪方法
技术领域：
本发明属于石油地球物理勘探领域，具体涉及一种基于波场分离原理提高地震资料信噪比的去噪方法。
背景技术：
随机噪声是地震资料处理中遇到的一类主要干扰波，在地震记录上主要表现为波形杂乱无章，频带很宽，视速度不确定，无一定的传播方向，在一定程度上能够满足统计规律。随机噪声的存在会降低地震资料的信噪比，影响资料处理的质量。压制随机噪声的方法主要分为3类I，基于预测技术的方法(f-x反褶积和投影反褶积等)；2，基于统计原理的方法(中值滤波和多项式拟合等)；3，基于数学变换的方法(小波包[I]、奇异值分解技术(SVD/K-L) [2] ,Radon变换和模式分解等(此处的[I]、[2]、[3]指参考文献，其中，[I]王振国，汪恩华.小波包相关阀值去噪[J].石油物探，2002，41(4) =400-405 ； [2]吕景贵，刘振彪，管叶君.压制叠前相干噪音的速度变换域滤波方法[J].石油物探，2001，40 (4) 94-99 ； [3]张军华，吕宁，雷凌，等.抛物线拉冬变换消除多次波的应用要素分析[J].石油地球物理勘探，2004，39 (4)398-405))。其中f-x(y)域反褶积是地震资料处理中压制随机噪声的主流技术。从20世纪70年代开始，SVD技术被广泛应用于信号处理，包括SVD (基于奇异值分解)技术压制噪声及有损压缩数字图像；噪声压制、倾角滤波、垂直地震剖面数据波场分离和剩余静校正。这些技术计算量大，需要复杂的修正才能适用于倾角数据。而且基于SVD的噪声压制技术是在t-x域中实现的，一般要求信号为水平，并不适合同相轴交错的情况。当有效信号同相轴呈倾斜和弯曲时，SVD会损害有效信号
发明内容
本发明的目的在于解决上述现有技术中存在的难题，提供一种基于波场分离原理提高地震资料信噪比的去噪方法，将数据分块技术与波场分离技术相结合来压制随机噪声，无需对同相轴进行任何相对时差校正，在假设线性同相轴条件下能适应同相轴交错的复杂情况。本发明是通过以下技术方案实现的一种基于波场分离原理提高地震资料信噪比的去噪方法，所述方法利用数据分块技术与波场分离技术相结合来压制随机噪声进而提高地震资料的信噪比；所述方法首先对地震数据进行分块处理，然后对每个数据块进行复数域SVD滤波，利用波场分离技术进行特征值分解，利用优选特征值重构数据，并替换原始频率切片进而实现压制随机噪声。所述方法包括以下步骤(I)按时间-空间方向对地震数据进行分块处理，获得一组数据块；
(2)复数域SVD滤波步骤将步骤⑴中分解得到的每一个数据块变换到频率域中，构建Hankel矩阵，然后利用波场分离技术对Hankel矩阵进行特征值分解，并从大到小选取I或2个奇异值，最后利用优选特征值重构原始数据，并替换原始频率切片进而实现压制随机噪声。
所述步骤(I)包括以下步骤(11)计算时空窗的个数N= (T-S-D-Off)/(1-0) W(1-1)(1-1)式中T为纵、横向最大长度向量；S为纵、横向起始向量；D为纵、横向延迟向量；0为纵、横向重叠百分比；w为时空窗纵、横向跨度；(12)计算起始位置T0 = S+D+(i-l)ff(l-0) I ^ i ^ N(1-2)(13)采用镶边函数对边界进行平滑过渡，所述镶边函数采用Hanning函数，即
(、—l + cos〔^) Ixl <αW(x) = \2{a J(1-3)
\x\>a 0 1 1(1-3)式中a为镶边带长度。(14)设分时窗操作算子为P，Cadzow滤波算子为F，则整个数据分块处理过程表示为d = [W1PW2FPjd(1-4)(1-4)式中=WJPW2表示分块正、反过程的加权系数算子，由(1-3)式计算得到的，d指的是原始地震数据。步骤⑵包括以下步骤(21)对每一个数据块进行傅里叶正变换；(22)对选定频率范围的，构建nXn的复数Hankel矩阵Α，其元素为每道的时间频率变换值；选择合适的奇异值个数k，利用SVD计算获得选取有限个奇异值后构建的矩阵Fk(A);沿反对角线平均，恢复矩阵Hankel结构；利用Fk(A)中的值替换每道Fourier变换值；(23)进行傅里叶反变换，获得计算结果。与现有技术相比，本发明的有益效果是I)与其它传统的SVD特征值技术不同，本发明的方法对同相轴倾斜或交错情况没有任何特殊要求，并且不需要对同相轴进行相对时差校正；2)假定特定时空范围内剖面含有有限个倾角，模型数据和实际数据处理结果表明，只取第一个特征值会对有效信号造成显著损害，而本发明取前3个显著特征值进行去噪效果比较稳健。3)在小时空窗内地震数据的倾角变换可以看作是稳态的，实际资料处理结果表明，本发明采用分块处理的方式可以克服数据倾角的非稳态变化。4) SVD和重构需要的时间随着数据体的增大而增大，本发明采用分块处理的方法可以大大缩短SVD和重构占用的时间。
图I是本发明方法实施例中的分块处理的示意图。图2是本发明方法实施例中的复数域波场分离的效果图。图3是采用现有技术中的常规f_x域反褶积滤波结果图。·
具体实施例方式下面结合附图对本发明作进一步详细描述一种提高地震资料信噪比的去噪方法，所述方法利用数据分块技术与波场分离技术相结合来压制随机噪声进而提高地震资料的信噪比；所述方法首先对地震数据进行分块处理，然后对每个数据块进行复数域SVD滤波，利用波场分离技术(又称为复数域奇异值分解技术、奇异值滤波、奇异值分解及重构，简称SVD)进行特征值分解，利用优选特征值重构数据，并替换原始频率切片进而实现压制随机噪声。所述方法包括以下步骤(I)数据分块处理步骤将地震数据按时间-空间方向对地震数据进行分块处理(图I给出的是分块的示意图)地震信号同相轴在不同的时空位置上倾角不同，从统计规律上来说，同相轴的倾角变化为非稳态的。应用奇异值滤波一般选取较小的秩数字(1-3)，即假定倾角个数不超过3，如果将其应用在整个剖面上，难免会损害有效信号。在数据处理过程中，本发明采用分块处理的方法来克服上述非稳态的矛盾。另外奇异值分解及重构(即SVD技术)所占用的时间随数据体阶数的增大呈级数增加，而本发明采用数据分块处理处理(又称为数据分时空窗处理、时窗分割处理等)可以大大缩短SVD所占用的时间。分块个数的计算公式为N= (T-S-D-Off)/(1-0) W (1-1)式中T为纵、横向最大长度向量；S为纵、横向起始向量；D为纵、横向延迟向量；0为纵、横向重叠百分比；W为时空窗纵、横向跨度。起始位置计算公式为T0 = S+D+(i-l)ff(l-0) I ^ i ^ N (1-2)为了减弱截断效应，在数据分块处理时必须对边界进行平滑过渡，采用的镶边函数为Hanning函数，即
,、—l + cos〔^) Ixl <αW(x) = \2)(1-3)
\x\>a 0 1 1式中a为镶边带长度。设分时窗操作算子为P，Cadzow滤波算子为F，则整个数据分块处理过程可以表示为d = [W1PW2FPjd(1-4)式中=W1和W2表示分块正、反过程的加权系数算子，由(1-3)式计算得到的，d指的是原始地震数据。(1-4)是用数学向量的形式描述“分块->处理->合并”的过程，描述的是整个过程。综上所述，利用(1-1)式首先计算每个空间方向的块数，然后利用(1-2)式计算每个分块的起始位置，同时利用(1-3)式计算每个分块的加权系数及反加权系数，这样就完成了数据分块处理。(2)复数域SVD滤波步骤将步骤⑴中分解得到的每一个数据块变换到频率域中，构建Hankel矩阵，然后利用复数域奇异值分解(SVD)技术对Hankel矩阵进行特征值分解，从大到小选取I或2个奇异值(又称为优选特征值)，最后重构数据。对于时间频率域中的数据，数据的有效秩数不超过数据的倾角个数。设有nXn空间网格，其空间坐标为(xp，yq), P = I…n，q= 1···η。设B为nXn维矩阵，其元素为每一道频率为ω的离散傅里叶变化值。如果只有k个倾角，矩阵B具有如下形式[Bi = V= V(、
L Xp,q11(2-1)
/=1 /=1式中，Ci和Cli为与倾角和频率有关的实数；bi为复数值，其幅度为该频率处倾斜同相轴的振幅。定义η个向量A和gi (i = Ρ··10 [f ]p =ρ = 1".η(2-2)[g/]g = Qiq q — X-'Yl(2-3)则有
k
_2] B = ^fgT(2-4)
/=1矩阵B的各列为fi(i = 1-k)的线性组合，因此矩阵B的秩最大为k，前k个特征值不为零，可以用前k个特征图像加权和重构数据。对于二维频率域地震数据，利用Hankel矩阵估算二维数据的有效秩数。设有N个地震道(此N与式(1-1)中的时空窗个数N不是同样的意思，此处的为标量，而上面的为矢量)(这些空间位置不必等间隔分布)，其空间坐标为xp，P = I…N。对于一个时间频率ω，其空间采样值分别为t1； t2, . . . . , tN假设地震信号具有k个倾角，根据f_x域地震信号具有预测性的特点可以将地震
信号表不为k个正弦信号之和k
_7] tp=Kbl(2-5)
7 = 1式中ai为复数屯为实数。数据Hankel矩阵形式为
权利要求
1.一种基于波场分离原理提高地震资料信噪比的去噪方法，其特征在于所述方法利用数据分块技术与波场分离技术相结合来压制随机噪声进而提高地震资料的信噪比；所述方法首先对地震数据进行分块处理，然后对每个数据块进行复数域SVD滤波，利用波场分离技术进行特征值分解，利用优选特征值重构数据，并替换原始频率切片进而实现压制随机噪声。
2.根据权利要求I所述的基于波场分离原理提高地震资料信噪比的去噪方法，其特征在于所述方法包括以下步骤 (1)按时间-空间方向对地震数据进行分块处理，获得一组数据块； (2)复数域SVD滤波步骤将步骤(I)中分解得到的每一个数据块变换到频率域中，构建Hankel矩阵，然后利用波场分离技术对Hankel矩阵进行特征值分解，并从大到小选取I或2个奇异值，最后利用优选特征值重构原始数据，并替换原始频率切片进而实现压制随机噪声。
3.根据权利要求I所述的基于波场分离原理提高地震资料信噪比的去噪方法，其特征在于所述步骤(I)包括以下步骤 (11)计算时空窗的个数N= (T-S-D-Off)/(1-0) W (1-1) (1-1)式中T为纵、横向最大长度向量；S为纵、横向起始向量；D为纵、横向延迟向量；O为纵、横向重叠百分比；W为时空窗纵、横向跨度； (12)计算起始位置T0 = S+D+(i-l)ff(l-0) I ^ i ^ N (1-2) (13)采用镶边函数对边界进行平滑过渡，所述镶边函数采用Hanning函数，即
4.根据权利要求I所述的基于波场分离原理提高地震资料信噪比的去噪方法，其特征在于所述步骤(2)包括以下步骤 (21)对每一个数据块进行傅里叶正变换； (22)对选定频率范围的，构建nXn的复数Hankel矩阵Α，其元素为每道的时间频率变换值；选择合适的奇异值个数k，利用SVD计算获得选取有限个奇异值后构建的矩阵Fk(A)；沿反对角线平均，恢复矩阵Hankel结构；利用Fk(A)中的值替换每道Fourier变换值； (23)进行傅里叶反变换，获得计算结果。
全文摘要
本发明提供了一种提高地震资料信噪比的去噪方法，属于石油地球物理勘探领域。本发明利用数据分块技术与波场分离技术相结合来压制随机噪声进而提高地震资料的信噪比；所述方法首先对地震数据进行分块处理，然后对每个数据块进行复数域SVD滤波，利用波场分离技术进行特征值分解，利用优选特征值重构数据，并替换原始频率切片进而实现压制随机噪声。与传统的假设同相轴水平的SVD技术相比，该方法无需对同相轴进行任何相对时差校正，且利用本方法的去噪效果比较稳健，本发明采用的分块处理技术能有效地克服实际地震资料中倾角变化的非稳态，且可以大大缩短SVD和重构占用的时间。
文档编号G01V1/36GK102854533SQ20111018427
公开日2013年1月2日 申请日期2011年7月1日 优先权日2011年7月1日
发明者崔树果, 郭全仕, 吕秋玲, 王跃, 崔连军 申请人:中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院