一种改进的独立分量分析算法在大地电磁去噪中的应用

曹小玲,严良俊

(1.长江大学油气资源与勘探技术教育部重点实验室,湖北武汉430100;2.非常规油气湖北省协同创新中心,湖北武汉430100;3.长江大学信息与数学学院,湖北荆州434023)

曹小玲,严良俊.一种改进的独立分量分析算法在大地电磁去噪中的应用[J].石油物探,2017,56(6):-897

CAO Xiaoling,YAN Liangjun.The application of an improved independent component analysis algorithm in magnetotelluric data denoising[J].Geophysical Prospecting for Petroleum,2017,56(6):-897

一种改进的独立分量分析算法在大地电磁去噪中的应用

曹小玲1,2,3,严良俊1,2

(1.长江大学油气资源与勘探技术教育部重点实验室,湖北武汉430100;2.非常规油气湖北省协同创新中心,湖北武汉430100;3.长江大学信息与数学学院,湖北荆州434023)

根据大地电磁噪声的特点和独立分量分析(ICA)中M-FastICA算法的优良性能,结合小波分析和盲源分离的相关理论,提出了一种改进的独立分量分析去噪方法。首先对观测信号进行多尺度小波分解,使信号从单道变成多道,以满足独立分量分析对观测信号的数目需求;然后采用M-FastICA算法对小波分解提取的多层高频分量进行独立分量分析以提取有效独立分量和特定独立分量;引入动态自适应因子来限制特定独立分量的权重以减小观测信号信噪比对去噪效果的影响;最后由小波低频分量和M-FastICA算法提取的两种独立分量共同构成恢复信号。模拟信号仿真实验表明,该方法的去噪性能优于传统小波阈值去噪方法。将该方法应用于实际大地电磁观测资料的去噪处理,无论是视电阻率曲线还是相位曲线,都比去噪前更加光滑和稳定,说明改进的独立分量分析算法能有效地去除大地电磁噪声。

独立分量分析;小波分析;大地电磁;去噪;M-FastICA算法

大地电磁测深(MT)是一种被动源电磁探测技术,通过在地表测量正交电、磁场分量的扰动值,得到地下介质的电性结构信息[1]。MT探测范围较宽,应用领域较广[2],但观测的天然电磁场信号较弱、频带较宽,极易受到各种噪声的干扰,给后续的数据处理、反演和地质解释结果带来误差。为此,地球物理研究人员提出了一系列去除大地电磁噪声的技术[3-7]。

传统的电磁数据去噪技术基于傅里叶变换——功率谱估计,因其无法表达信号的时域局部性质,人们先后应用小波分析方法消除大地电磁的噪声[8-13]。2004年,王书明等通过研究不同地区实测MT信号的特征,指出MT信号一般具有非高斯性、非线性和非最小相位性质[14-15],需要利用新的信号处理方法来解决大地电磁的噪声问题。近年来,独立分量分析(ICA)技术[16-19]被广泛应用于地震勘探、语音信号处理、图像处理等领域的信号分析和信号处理,在这些领域的信号去噪中表现尤其突出。该方法不需要源信号的相关信息,也不需要观测信号中源信号的混合方式,通过对观测信号的分析和处理推测甚至恢复出源信号。这个特点对于地球物理勘探来说意义重大,因为源信号的相关信息(先验知识)不太容易获得,即使获得也成本高昂。独立分量分析要求待处理信号为非依赖性、非高斯性和非白性,而大地电磁观测信号恰好符合这些要求。另外,独立分量分析一般要求处理的对象是由一组在统计意义上相互独立的源信号构成的混合信号,而大地电磁观测信号中噪声干扰和信号分别由不同的信源产生,彼此相互独立,符合独立分量分析对处理对象的这一要求。因此,独立分量分析技术可以用来去除大地电磁观测信号中的噪声干扰。该方法是在统计独立意义下对混合信号进行去噪,显然比仅仅依赖于不相关性度量的其它去噪方法效果更好。本文根据大地电磁测深数据量大的特点以及独立分量分析中M-FastICA算法的优良性能,结合小波分析和盲源分离的相关理论,研究并提出了一种改进的独立分量分析算法,利用实际大地电磁资料验证了算法的有效性。

1 独立分量分析及M-FastICA算法

独立分量分析是20世纪90年代发展起来的一种新的统计方法。它假定第i个观测信号xi是由n个相互独立的未知源信号sj混合而成:

(1)

假定每个观测信号xi和未知源信号sj都是随机变量,矢量x表示观测信号变量(x1,x2,…,xm)T,矢量s表示源信号变量(s1,s2,…,sn)T(m≥n),矩阵Am×n表示混合矩阵(aij)m×n,则公式(1)可用矢量-矩阵形式表示为:

(2)

公式(2)也称为独立分量分析模型[17]。

独立分量分析的目标是寻求分离矩阵Wm×n,得到独立分量s的估计y=WTx=WTAs,使y具有最大的非高斯性。FastICA算法[16-19]是独立分量分析的一种快速算法,该算法采用牛顿迭代法对观测变量x的大量采样点进行批处理,每次从观测信号中分离出一个独立分量,采用负熵(Negentropy)[17]作为非高斯性度量的标准,应用固定点迭代理论寻找WTx的非高斯性最大值。

(3)

实际FastICA算法运行过程中,由于采用牛顿迭代法,收敛速度至多是三阶,可能会出现收敛效果不理想或者收敛速度慢的情况,特别是对于海量数据的处理更会如此。因此,本文采用独立分量分析的M-FastICA算法[20]来对数据量庞大的大地电磁观测信号进行处理。迭代过程为:

(4)

式中,F(w)=E{xg(wTx)}-βw,其雅可比矩阵函数为J(w)=E{xxTg′(wTx)}-βI。

可以证明,(4)式的收敛阶数为m+1,且每m次迭代只需计算一次雅可比矩阵J(w),因此提高了独立分量分析的速度。一般取m=2时的(4)式为M-FastICA算法的迭代公式,即:

(5)

2 改进的独立分量分析算法

笔者在利用小波分析方法去除大地电磁噪声的过程中,参照文献[21]并进行多次试验后发现,小波变换方法对于较高信噪比的含噪声信号去噪效果较好,但对较低信噪比的含噪声信号去噪效果较差。大地电磁信号本身较弱,受到的噪声干扰相对来说可能比较强,若单独采用小波分析方法去噪可能得不到理想的效果。而属于盲源分离的独立分量分析方法对于低信噪比的含噪声信号去噪效果较好,故本文将M-FastICA算法与小波分析方法结合起来进行去噪处理。这里进行小波分析的目的有两个:一是在去噪前期消除部分噪声,二是解决将单个时间序列视为观测信号时观测信号数目不足的问题。

根据盲源分离理论[22],如果所获得的信号可以近似地认为是由各个独立的源信号经过混合之后所得的结果,那么就可以运用独立分量分析提取出源信号。就大地电磁测深而言,虽然我们无法证明源信号和噪声信号之间是严格独立的,但是因为这两种信号源的产生方式不同,因而它们之间至少是独立的。同

时,由于大地电磁源信号与噪声信号的来源不相同,导致它们的属性也不相同,因此,我们可以运用独立分量分析来对大地电磁观测信号中的源信号和噪声进行分离。

运用独立分量分析进行信噪分离时要求采集到的信号的数目大于或者等于源信号的数目[16-19],也就是说,必须同时获取至少两组数据,才能将所得数据分离为信号和噪声两组数据。而在实际野外作业中,要同时获得两组及两组以上的数据需要耗费大量的人力和物力。受文献[23—24]启发,我们对独立分量分析算法进行了改进,采用单通道模式对大地电磁时间序列信号进行去噪处理。首先对观测信号进行多尺度小波分解,使信号从单道变成多道,以满足独立分量分析对观测信号的数目需求;然后对小波分解提取的多层高频分量进行独立分量分析,以提取源信号和噪声。为了减小独立分量分析算法分离信号后幅值的不确定性,本文引入动态自适应因子α来对特定源信号的权重加以限制。α的取值根据观测信号的幅值确定,即α=max(A)/β,β=5×round{lgmax(A)}。这里A表示观测信号所有幅值的集合,round(·)表示取整数。为了尽可能地保留原始信号的基本形态，因此将小波低频分量信号加入到恢复信号中。

改进的独立分量分析算法流程如图1所示。具体实现步骤如下:

1) 将观测信号进行多层小波分解,并分别提取多层小波高频系数和某一特定层的小波低频分量,同时观察信号的幅值,确定自适应因子α的取值;

2) 运用M-FastICA算法对多层小波高频系数进行独立分量分析,提取出有效源信号;

3) 按照Sde-noised=SlW⊕0.5αSsICA⊕0.5SeICA的形式重构恢复信号,其中Sde-noised为恢复信号,SlW为特定层小波低频分量信号,⊕为求信号的合成信号,SsICA为特定独立分量信号,SeICA为有效独立分量信号(有效独立分量指M-FastICA算法分离出来的独立分量,特定独立分量指运用聚类规则从有效独立分量中挑选出来的某一个特定的独立分量)。

图1 改进的独立分量分析算法流程

3 仿真实验

通过仿真实验对改进的独立分量分析算法进行了测试。仿真实验分别采用周期和非周期两种信号,原始周期信号S1为10×sin(5t)+50×cos(3t),原始非周期信号S2为cum(t)。噪声信号为由蒙特卡洛法构造的随机信号。

在原始信号中加入不同量的噪声,分别采用如下7种方法进行去噪:①强制去噪法;②中间阈值法;③mini-max法;④自适应默认阈值去噪法;⑤史坦无偏似然估计去噪法;⑥重构信号为Sde-noised=SlW⊕SeICA的去噪法(公式中各符号含义同前);⑦本文改进的独立分量分析法(为了对比,这里采用的小波基和小波分解层数分别与方法①～⑤相同)。方法①～⑤为传统的小波阈值去噪方法[25-26]。

1) 强制去噪法。将小波分解结构中的高频系数全部置为0,即将高频部分全部滤除。

2) 中间阈值法。设小波分解得到的各层小波系数中最大的系数为F,最小的系数为f,则阈值为λ=(F+f)/2。

3) mini-max法,即最小极大方差阈值法。它是一种固定阈值的方法,设N为小波系数的长度,则阈值为:

5) 史坦无偏似然估计去噪法。采用史坦无偏似然估计原理进行自适应阈值选择。

图2和图3分别给出了信号S1和S2加入噪声后信噪比为RSN=15时的实验结果。

图2 信号S1(这里S1为含噪信号)在信噪比为RSN=15时的仿真实验去噪结果 a 信号S1的构成; b 对信号S1进行5层小波分解后的结果; c 对高频小波分量进行M-FastICA处理的结果; d 对信号S1分别运用7种方法去噪后的结果

图3 信号S2(这里S2为含噪信号)在RSN=15时的仿真实验去噪结果 a 信号S2的构成; b 对信号S2进行5层小波分解后的结果; c 对高频小波分量进行M-FastICA处理的结果; d 对信号S2分别运用7种方法去噪后的结果

为了衡量重构的恢复信号与原始信号的近似程度,这里根据系统聚类法的基本思想,采用相似系数来对两者的近似程度进行度量:

(6)

表1和表2分别是运用7种去噪方法对信号S1和S2去噪后的结果(即恢复信号)与对应的原始信号的相似系数。考虑到大地电磁信号微弱,获得的观测信号可能包含较多的噪声,故这里展示了较低信噪比情况下的去噪结果。可以看出,无论是信号S1还是信号S2,在观测信号的信噪比较低或极低的情况下,本文改进的独立分量分析方法去噪效果都比其它方法更为显著,也更为稳定。这是因为本文方法结合了小波分析和M-FastICA算法各自的优势,且所采用的恢复信号既包含低频小波系数,又包含M-FastICA算法分离出的有效源信号和特定源信号,同时通过自适应因子动态改变M-FastICA算法的分离效果,使噪声消除结果受观测信号信噪比的影响大大减小。另外,方法⑥在处理非周期信号时与方法⑦效果近似,但在处理周期信号时明显不如方法⑦,原因是方法⑥的恢复信号舍弃了特定源信号。

表1 信号S1运用7种方法去噪后的相似系数

表2 信号S2运用7种方法去噪后的相似系数

4 实际观测信号去噪处理

将本文改进的独立分量分析去噪方法应用于实际大地电磁观测信号的噪声消除,绘制视电阻率曲线和相位曲线,选取含噪声较多的两个测点进行对比分析(其中测点2中的噪声含量明显高于测点1),如图4 和图5所示。测点1应用本文方法去噪后,无论是视电阻率曲线还是相位曲线,都比去噪前更加光滑、稳定,说明本文改进的独立分量分析方法去噪效果较好。测点2因为含有的噪声多,原始信号几乎被淹没,但应用本文方法去噪后,视电阻率曲线和相位曲线中的误差棒(图中红点和黑点上的竖线)相对去噪前明显变短,甚至被消除很多,说明本文方法能比较有效地去除信噪比极低情况下的信号噪声。

图4 实际观测信号去噪前、后视电阻率曲线对比 a 去噪前的测点1; b 去噪后的测点1; c 去噪前的测点2; d 去噪后的测点2

图5 实际观测信号去噪前后相位曲线对比 a 去噪前的测点1; b 去噪后的测点1; c 去噪前的测点2; d 去噪后的测点2

5 结束语

本文根据大地电磁测深数据量大的特点以及独立分量分析M-FastICA算法的优良性能,结合小波分析和盲源分离的相关理论,提出了一种改进的独立分量分析去噪方法。为了减小观测信号信噪比对去噪效果的影响,引入动态自适应因子来限制特定独立分量的权重,恢复信号由小波低频分量和M-FastICA算法提取的两种独立分量共同构成。模拟信号仿真实验结果表明,本文方法的去噪性能优于传统的小波阈值去噪方法。将本文方法应用于实际大地电磁观测信号去噪处理,无论是噪声较少还是噪声较多的情况,都能比较有效地去除噪声,去噪后视电阻率曲线和相位曲线比去噪前更加光滑、稳定,误差棒相对去噪前明显变短、甚至消除。

[1] 石应俊,刘国栋,吴广耀,等.大地电磁测深法教程[M].北京:地震出版社,1985:1-6

SHI Y J,LIU G D,WU G Y,et al.Magnetotelluric sounding tutorial[M].Beijing:Earthquake Publishing House,1985:1-6

[2] 高利华.大地电磁测深在北天山前缘油气勘探中的应用[J].石油物探,2006,45(4):427-430

GAO L H.Application of magnetotelluric depth sounding in North Tianshan foothill belt[J].Geophysical Prospecting for Petroleum,2006,45(4):427-430

[3] 严良俊,胡文宝,陈清礼,等.远参考MT方法及其在南方强干扰地区的应用[J].石油天然气学报,1998,20(4):34-38

YAN L J,HU W B,CHEN Q L,et al.Application of remote reference MT to noisy area[J].Journal of Oil and Gas,1998,20(4):34-38

[4] 宋维琪,仝兆歧.利用泛克里格方法消除MT噪音[J].石油物探,2001,40(2):138-142

SONG W Q,TONG Z Q.Using universal Krige method to remove noise for MT data[J].Geophysical Prospecting for Petroleum,2001,40(2):138-142

[5] 张全胜,杨生.大地电磁测深资料去噪方法应用研究[J].石油物探,2002,41(4):493-499

ZHANG Q S,YANG S.An application study of noise elimination for magnetotelluric sounding data[J].Geophysical Prospecting for Petroleum,2002,41(4):493-499

[6] 蔡剑华,熊锐.基于频率切片小波变换的时频分析与MT信号去噪[J].石油物探,2016,55(6):904-912

CAI J H,XIONG R.Magnetotelluric data denosing based on time-frequency analysis of the frequency slice wavelet transform[J].Geophysical Prospecting for Petroleum,2016,55(6):904-912

[7] 陈高,金祖发,马永生,等.大地电磁测深远参考技术及应用效果[J].石油物探,2001,40(3):112-117

CHEN G,JIN Z F,MA Y S,et al.Remote reference technique in magnetotelluric sounding and its application[J].Geophysical Prospecting for Petroleum,2001,40(3):112-117

[8] TRAD D O,TRAVASSOS J M.Wavelet filtering of magnetotelluricdata[J].Geophysics,2000,65(2):482-491

[9] MANOJ C,NAGARAJAN N.The application of artificial neural networks to magnetotelluric time-series analysis[J].Geophysical Journal International,2003,153(2):409-423

[10] WECKMANN U,MAGUNIA A,RITTER O.Effective noise separation for magnetotelluric single site data processing using a frequency domain selection scheme[J].Geophysical Journal International,2005,161(3):635-652

[11] 何兰芳,王绪本,何展翔,等.MT时间序列的小波去噪分析[J].地震地质,2001,23(2):222-226

HE L F,WANG X B,HE Z X,et al.Wavelet-based de-noising of MT time series[J].Journal of Earthquake Geology,2001,23(2):222-226

[12] 凌振宝,王沛元,万云霞,等.强人文干扰环境的电磁数据小波去噪方法研究[J].地球物理学报,2016,59(9):3436-3447

LING Z B,WANG P Y,WAN Y X,et al.A combined wavelet transform algorithm used for de-noising magnetotellurics data in the strong human noise[J].Chinese Journal of Geophysics,2016,59(9):3436-3447

[13] 蔡剑华,肖晓.基于小波自适应阈值去噪的MT信号处理方法[J].地球物理学进展,2015,30(6):2433-2439

CAI J H,XIAO X.Method of processing magnetotelluric signal based on the adaptive threshold wavelet[J].Progress in Geophysics,2015,30(6):2433-2439

[14] 王书明,王家映.大地电磁信号统计特征分析[J].地震学报,2004,26(6):669-674

WANG S M,WANG J Y.Analysis on statistic characteristics of magnetotelluric signal[J].Journal of Earthquake,2004,26(6):669-674

[15] 王书明,王家映.关于大地电磁信号非最小相位性的讨论[J].地球物理学进展,2004,19(2):216-221

WANG S M,WANG J Y.Discussion on the non-minimum phase of magnetotelluric signals[J].Progress in Geophysics,2004,19(2):216-221

[16] HYVARINEN A,OJA E.A fast fixed-point algorithm for independent component analysis[J].Neural Computation,1997,9(7):1483-1492

[17] HYVARINEN A,OJA E.Independent component analysis:algorithms and applications[J].Neural Networks,2000,13(4/5):411-430

[18] 杨福生.独立分量分析的原理与应用[M].北京:清华大学出版社,2006:33-46

YANG F S.Principle and application of independent component analysis[M].Beijing:Tsinghua University Press,2006:33-46

[19] 海韦里恩.独立成分分析[M].北京:电子工业出版社,2007:6-11,138-153

HYVARINEN A.Independent component analysis[M].Beijing:Electronic Industry Press,2007:6-11,138-153

[20] 曾生根,朱宁波,包晔,等.一种改进的快速独立分量分析算法及其在图像分离中的应用[J].中国图像图形学报,2003,8(10):1159-1165

ZENG S G,ZHU N B,BAO Y,et al.A modified fast independent component analysis and its application to image separation[J].Journal of Image and Graphics,2003,8(10) 1159-1165

[21] 程娇.基于小波变换和独立分量分析的去噪方法研究[D].上海:复旦大学,2010

CHENG J.Research on de-noising method based on wavelet transform and independent component analysis[D].Shanghai:Fudan University,2010

[22] 余先川,胡丹.盲源分离理论与应用[M].北京:科学出版社,2011:55-68

YU X C,HU D.Theory and application of blind source separation[M].Beijing:Science Press,2011:55-68

[23] 陈艳,何英,朱小会.基于小波变换的独立分量分析及其在图像分离中的应用[J].现代电子技术,2007,30(24):131-134

CHEN Y,HE Y,ZHU X H.ICA based wavelet transform and its application to image separation[J].Journal of Modern Electronic Technology,2007,30(24):131-134

[24] 于刚,周以齐.单通道盲源分离算法及其在工程机械振源分析中的应用[J].机械工程学报,2016,52(10):1-8

YU G,ZHOU Y Q.Single-channel blind source separation and its application in analyzing vibration of engineering machinery[J].Journal of Mechanical Engineering,2016,52(10):1-8

[25] 余晃晶.小波降噪阈值选取的研究[J].绍兴文理学院学报,2004,24(9):34-38

YU H J.Selection of a wavelet noise-reduction threshold[J].Journal of Shaoxing University,2004,24(9):34-38

[26] DONOHO D L,JOHNSTONE I M.Ideal spatial adaption by wavelet shrinkage[J].Biometrika,1994,81(3):425-455

(编辑：戴春秋)

Theapplicationofanimprovedindependentcomponentanalysisalgorithminmagnetotelluricdatadenoising

CAO Xiaoling1,2,3,YAN Liangjun1,2

(1.KeyLaboratoryofExplorationTechnologiesforOilandGasResourcesinMinistryofEducation,Yangtzeuniversity,Wuhan430100,China;2.HubeiCooperationInnovationCenterofUnconventionalOilandGas,Wuhan430100,China;3.SchoolofInformationandMathematics,Yangtzeuniversity,Jingzhou434023,China)

Combined with the wavelet analysis theory and blind source separation,an improved independent component analysis algorithm for magnetotelluric data denoising is proposed,which is based on the M-FastICA algorithm,an improvement on the FastICA algorithm.First,the multi-scale wavelet decomposition is performed to convert the observed signal from single channel to multiple channels,to meet the quantity demand of the independent component analysis.Then,the M-FastICA algorithm is adopted to process the multi-layer high frequency components extracted by wavelet decomposition.This process extracts the effective independent components and specific independent components.Then,the dynamic adaptive factor is introduced to limit the weight of the specific independent components and reduce the SNR influence in the observed data on the denoising process.Finally,the low frequency components of the wavelet and two types of independent components extracted using the M-FastICA algorithm together constitute the recovery signal.Simulation test results on an analog signal show that the denoising performance of the proposed method is better than the traditional wavelet threshold denoising method.The application test on field magnetotelluric observation data shows that both the apparent resistivity and phase curves are smoother and more stable than that before denoising.These results illustrate that the proposed method can effectively remove magnetotelluric noise.

independent component analysis,wavelet analysis,magnetotelluric,denoising,M-FastICA algorithm

2017-04-11;改回日期2017-08-27。

曹小玲(1981—),女,讲师,博士,主要从事大地电磁数据处理方法研究工作。

国家自然科学基金(41274082,U1562109)、长江大学长江青年基金(2015cqn76)、长江大学重磁电勘探研究中心创新基金(7011201803xm)联合资助。

This research is financially supported by the National Natural Science Foundation of China (Grant Nos.41274082,U1562109),the Yangtze Youth Fund in Yangtze university (Grant No.2015cqn76),and the Innovation Fund of Gravity-magnetic-electric exploration research center in Yangtze university (Grant No.7011201803xm).

P631

A

1000-1441(2017)06-0890-08

10.3969/j.issn.1000-1441.2017.06.015