航空电磁时间域数据工频干扰研究与去除.pdf
《航空电磁时间域数据工频干扰研究与去除.pdf》由会员分享,可在线阅读,更多相关《航空电磁时间域数据工频干扰研究与去除.pdf(7页珍藏版)》请在咨信网上搜索。
1、 第 47 卷第 5 期物 探 与 化 探Vol.47,No.5 2023 年 10 月GEOPHYSICAL&GEOCHEMICAL EXPLORATION Oct.,2023doi:10.11720/wtyht.2023.1458黄威,肖都,贲放,等.航空电磁时间域数据工频干扰研究与去除J.物探与化探,2023,47(5):1281-1287.http:/doi.org/10.11720/wtyht.2023.1458Huang W,Xiao D,Ben F,et al.Study and removal of power frequency interference in time-do
2、main airborne electromagnetic dataJ.Geophysical and Geo-chemical Exploration,2023,47(5):1281-1287.http:/doi.org/10.11720/wtyht.2023.1458航空电磁时间域数据工频干扰研究与去除黄威1,2,3,肖都1,2,3,贲放1,2,3,李军峰1,2,3,李飞1,2,3,胥值礼1,2,3(1.自然资源部 地球物理电磁法探测技术重点实验室,河北 廊坊 065000;2.中国地质科学院 地球物理地球化学勘查研究所,河北 廊坊 065000;3.国家现代地质勘查工程技术研究中心,河北
3、 廊坊 065000)摘 要:时间域航空电磁系统在飞行作业过程中,信号受人文噪声中的高压输电线影响,导致测得的响应中含有工频干扰信息而不能真实地反映地下介质分布情况,因此本文开展了工频干扰噪声去除研究,提出一种可直接在时间域中对数据进行处理的梳状滤波器,该滤波器的提出不仅节省了时频和频时转换的时间,也避免了繁琐的转换步骤,降低了转换过程中可能出现误差的风险。利用梳状滤波器的关键点是通过理想阻带来抑制工频频率信号及其 k 次谐波,达到对工频干扰压制的目的。本文不仅研究了 IIR 型和 FIR 型梳状滤波器,同时对比了单一陷波滤波器对工频干扰的去除效果,期望通过该研究可以为时间域航空电磁数据处理提
4、供一种高效的工频干扰去除手段。关键词:时间域航空电磁法;梳状滤波器;工频干扰中图分类号:P631 文献标识码:A 文章编号:1000-8918(2023)05-1281-07收稿日期:2022-09-21;修回日期:2023-02-08基金项目:国家自然科学基金(42204147);基本科研业务费专项经费(AS2022J01);中国地质调查局项目(DD20230245)第一作者:黄威(1989-),女,高级工程师,2016 年博士毕业于吉林大学,主要从事航空电磁正反演理论及数据处理与解释研究工作。Email:huangwei2012511 通讯作者:肖都(1976-),男,教授级高级工程师,主
5、要从事地球物理勘查技术研究工作。Email:xdu 0 引言在数据采集过程中,时间域航空电磁信号易受到来自各方面的干扰,导致有用信号淹没在噪声中无法正确解释出地下介质的电性信息,主要噪声有运动噪声、工频干扰和天电噪声等。工频干扰的频率根据不同国家和地区交流电频率而定,我国电网标准频率是 50 Hz,习惯上称为工频。工频干扰是指电力系统运行过程中对地球物理探测信号引起的一种干扰,是以利用电磁感应电流和接地漏电电流等方式,存在于交流高压输电、变电和用电设施及大型建筑物周边的一种干扰1。根据目标任务的不同,工频干扰的作用也不同,如若目的是得到二次场的响应曲线,工频干扰则被视为电磁信号的噪声予以去除;
6、如若判断高压线接地情况时,可利用该工频干扰信号来分析高压线塔基接地情况2。在时间域航空电磁对地进行勘探时,工频干扰视为噪声被去除。工频干扰的传播方式分为空间辐射和传导 2 种。在电磁勘探中,工频干扰以电磁波的辐射形式对测量信号产生正弦波或其他信号与正弦波的叠加干扰,在响应中会出现 50 Hz 及其谐波振荡,导致数据曲线形态发生畸变。正确、精准地对工频干扰进行去除,可以提高数据质量,是正确反映地下电性构造的必要手段。根据工频干扰对不同研究领域信号影响程度的不同,学者们纷纷提出了解决工频干扰的方法。在医学生物信号领域中,去除心电信号中工频干扰较常用的方法是单一陷波滤波器、频率域梳状滤波器,后出现了
7、 LMS 和 RLS 自适应滤波方法3。在机械振动研究中,学者利用奇异向量和奇异值分解的方法对工频干扰进行压制4。在地学领域中,地震勘探处理主要在频率域里进行去除,常规的方法为单物 探 与 化 探47 卷 一陷波滤波和梳状滤波器,但由于地震工作频率5与工频的频率相近,不适宜使用陷波器对工频干扰进行去除,在研究了陷波滤波器对信号损伤基础上,陆续出现了利用 s 变换6、频谱剪切7等技术方法;在时域中出现了正余弦函数逼近法、余弦函数逼近法、自相关算法8和时域压制等方法。在大地电磁工作中,传统方法主要通过远程参考处理、编辑时间序列和基于统计参数的数据段选择来压制和消除工频干扰,后面又出现了频率域稀疏分
8、解9、神经网络10等方法。在瞬变电磁法中,通过双极性同步采样来抑制噪声的前提下11,一般采用改变叠加次数或者更换采样基频的方式来压制工频干扰,或者将时间域数据转换到频率域进行滤波12,也有使用自适应滤波器的方法13对工频干扰进行去除。时间域航空电磁采样密集会导致数据量庞大,在对原始数据流进行处理时保证滤波方法准确的前提下,需加快处理速度。自适应滤波处理方法需要经过大量数据对滤波系数进行调整,因此减缓了数据处理的速度,显然该方法不适合应用于时间域航空电磁的野外数据实时处理中。综上所述,本文提出一种可去除时间域航空电磁数据工频干扰的方法 梳状滤波器。首先,通过单一陷波器的研究了解陷波原理,通过分析
9、其对50 Hz 干扰的去除效果,进而扩展到梳状滤波器对50 Hz 及谐波干扰的去除。研究发现,梳状滤波器可设计成 FIR 型和 IIR 型,并通过实测数据处理来观察两者的滤波效果。1 工频干扰去除方法研究时间域航空电磁采集信号中含有多种噪声,其中工频干扰、运动噪声对信号影响较大。运动噪声对信号的影响主要是低频的非平稳振荡,本文为减小其干扰影响,选取了被运动噪声影响微弱的信号进行研究14。如图 1 所示,图中分别给出了总场数据无工频干扰和有工频干扰噪声的曲线变化情况,从中可见曲线较平稳,运动噪声影响非常小。图中信号为总场信息,包含了背景场信号和二次场信号。一般在数据预处理过程中最后进行背景场去除
10、得到二次场信号,工频干扰去除属于数据预处理前期工作,所以本文中研究工频干扰去除的数据示例为总场数据。a无工频干扰的信号特征;b工频干扰噪声特征asignal characteristics without power frequency interference;bpower frequency interference noise characteristics图 1 总场响应无干扰和有工频干扰噪声形态Fig.1 Power frequency interference noise patterns 由于工业用电对电磁信号产生了噪声干扰,使得信号产生畸变。工频干扰在信号中主要体现为 50 H
11、z 及其谐波的振荡,在频谱中主要影响 50 Hz 及 k 次频率周围幅值的变化,导致瞬变电磁有用信号被淹没,出现了一定幅度的正弦或余弦形态噪声,幅度的大小取决于信号接收器与干扰的距离。图 1b 中实测数据的基频为 12.5 Hz,一个周期内含有 4 个周期的工频干扰,即频率为 50 Hz。陷波滤波器、梳状滤波器可对工频干扰产生抑制作用,因此,本文对陷波滤波器和梳状滤波器进行详细的分析。1.1 陷波滤波器陷波滤波器是通过对输入信号在某一频点快速衰减,从而达到对该信号频点进行抑制的一种滤波器。陷波滤波器是一种带阻滤波器,其阻带非常狭窄,理想陷波滤波器的频率响应为:H ej()=0,=01,其他(1
12、)其中:=2f,f 是陷波频率。如图 2 所示,陷波频率点的信号响应为 0,其他频点的响应值为 1。滤波器中3 dB 频带的取值影响陷波的阻带宽度,3 dB 频带2821 5 期黄威等:航空电磁时间域数据工频干扰研究与去除取值越小,陷波滤波器越接近理想结果,其滤波效果越好。为此,利用陷波滤波器去除 50 Hz 工频干扰,给出 z 域中陷波滤波器的二阶传递函数:H z()=K1-2cos0z-1+z-21-2rcos0z-1+r2z-2,0 r 1,(2)式中:K 为滤波器的增益,由幅度响应来确定;r 为陷波滤波器的极半径;0=2f0,f0是陷波频点。将系统转换到时域后,再进行逆变换可得到对时间
13、域数据进行处理的方法:y n()=a0 x n()+a1x n-1()+a2x n-2()+b1y n-1()+b2y n-2(),(3)a0=K,a1=-Kcos0,a2=K,b1=2rcos0,b2=-r2,(4)式中:x(n)可视为时间序列输入信号;y(n)为输出信号,根据式(3)将待滤波信号进行输入,即可得到去噪后结果。图 2 50 Hz 频点陷波滤波器的特性曲线Fig.2 The characteristic curves of notch filter at 50 Hz1.2 梳状滤波器作为单一陷波滤波器的特殊形式,梳状滤波器是只允许某些特定频率范围信号通过的,由许多按一定频率间隔
14、相同排列的通带和阻带组成的一种滤波器,其特性曲线像梳子一样,故称为梳状滤波器。梳状滤波器不仅可以去除某一固定频率干扰,同时可去除该频率产生的 k 次谐波,即用来抑制周期性的噪声。梳状滤波器根据冲激响应的不同,分为有限冲激响应(FIR)滤波器和无限冲激响应(IIR)滤波器。结构形式上 FIR 型滤波器的输出取决于当前和过去的输入信号值,而 IIR 型滤波器的输出不仅取决于当前和过去的输入信号,也取决于过去的输出信号值;相位特性上,前者相位是线性的,后者相位是非线性的。根据 2 种类型滤波器的特点,本文分别对FIR 型和 IIR 型梳状滤波器进行研究。首先,给出 FIR 型梳状滤波器的系统函数:H
15、 z()=1-z-N()/2,(5)将其转换为时间域的表达式:y n()=x n()-x n-N()2,(6)式中:N 为滤波器的零点个数;N 为信号的半周期采样点数。从式(6)中可以看出,输出信号只取决于当前和过去的输入信号。通过对该系统函数幅频响应的分析(采样率为100 kHz,3 dB 带宽为 2,N=2 000),如图 3a 中的蓝色曲线所示,虽然 FIR 型梳状滤波器幅频特性的过渡带较为平缓,但可以看出对于 50 Hz 及其谐频可以起到抑制作用。同时,设计了IIR型梳状滤波器,aFIR 型梳状滤波器幅频特征;bIIR 型梳状滤波器幅频特征athe amplitude-frequenc
16、y characteristics of FIR comb filter;bthe amplitude-frequency characteristics of IIR comb filter图 3 FIR 和 IIR 梳状滤波器特性曲线Fig.3 The characteristic curves of FIR and IIR comb filters3821物 探 与 化 探47 卷 通过 3 dB 带宽来限制阻带的宽窄度,其系统函数为:H z()=1+rN21-z-N1-rNz-N,其中令 a=1+rN2,b=rN,(7)将其进行逆变换,并转到时间域为:y n()=ax n()-x n-
17、N()+by n-N(),(8)式(8)中明确了各参量之间的关系:N 为滤波器的零点个数;N为信号的半周期采样点数。IIR 型梳状滤波器的输出信号不仅取决于输入信号 x(n)和x(n-N),同时与输出信号 y(n-N)也相关。由此可以给出该滤波器的信号流(图 4)。图 4 IIR 型梳状滤波器的信号流Fig.4 The signal flow diagram of IIR comb filter 该信号流较清晰地从另一角度展示了各参量间的输入、输出关系。通过对时间域函数特性曲线进行分析,如图 3b 中红色曲线所示,可见 IIR 型梳状滤波器的阻带也可以很好地限制 50 Hz 及其 k 次谐波的
18、通过。2 工频干扰去除示例分析实测数据预处理中包含:工频干扰、运动噪声、天电噪声、人文干扰、背景场去除、叠加抽道等技术11-13,第一步去噪工作即工频干扰去除。为减小其他噪声对本文研究内容的影响,选取其他噪声较小的数据段作为本文的主要研究对象。实测数据中响应为总场值,包含了背景场值和二次场值,因二次场和背景场幅值差异较大,因此供电阶段曲线形态基本反映的是背景场值变化;断电早期受背景场影响不能直观地看出二次场的涡流反冲形态,需进行背景场去除后,才能进行直观显示;除早期外断电响外应基本为二次场值。本文选取由中国地质科学院地球物理地球化学勘查研究所自主研发的固定翼时间域航空电磁系统(iFTEM-II
19、)对在吉林省大安市龙沼镇地区的实测飞行数据进行研究。系统数据采样率为 100 kHz,基频为 12.5 Hz,发射波形采用双极性半正弦波,供电时间为 4.15 ms。首先,利用陷波滤波器对工频干扰进行去除。图 5 和图 6 分别给出了陷波器滤波前后的曲线形态和幅频特性,黑色曲线是原始数据,绿色曲线是滤波后数据。从图 5 中可以看出,工频干扰对数据产生较大的影响,去除工频干扰前的数据(黑色曲线)形态已被工频干扰所污染。利用陷波滤波器进行去除后(绿色曲线),响应中存在某些频率的工频干扰噪声,曲线形态仍处于被噪声破坏的形态。从图 6 幅频特征中可以看出,滤波后仅压制了 50 Hz 的噪声干扰,响应中
20、仍存在 50 Hz 及谐波噪声(红色圆圈区域)。因此说明,陷波滤波器在工频干扰去除方面存在着不足,需要寻找其他更好的方法对工频干扰进行去除。图 5 陷波滤波器去除工频干扰前后效果对比Fig.5 The comparison of the effect of notch filter removing power-line interference图 6 陷波滤波器去噪前后幅值特征Fig.6 The amplitude-frequency characteristics of notch filters4821 5 期黄威等:航空电磁时间域数据工频干扰研究与去除 其次,利用 FIR 型和 IIR
21、 型梳状滤波器对工频干扰进行去除。图 7 给出了 FIR 型和 IIR 型梳状滤波器数据处理前后的效果,黑色曲线是原始数据,天蓝色曲线是 FIR 型梳状滤波器滤波后数据,红色虚线是 IIR 型梳状滤波器滤波后数据(极点半径 r 为0.95)。图中可以看出,2 种滤波器较好地去除了原数据中存在的工频干扰,两者滤波效果相同,曲线重合。这里 FIR 型采用了 N=2 000 的零点个数。供电阶段滤波效果较小,因总场响应值远大于工频干扰,滤波后仔细观察可看到峰值存在一定变化,即去除了工频干扰的影响。断电后总场响应较工频干扰来说幅值小很多,因此滤波后,较明显地看出工频干扰被去除。图 7 FIR 型和 I
22、IR 型梳状滤波器滤波效果对比Fig.7 The comparison of effects for FIR and IIR comb filter 为更好地分析 2 种滤波器对工频干扰的去除效果,对滤波前后数据进行傅里叶变换,观察其频谱特征。图 8 给出了滤波前后数据的频谱特征,黑色曲线是原始数据频谱,红色曲线是 FIR 型和 IIR 型梳状滤波器滤波后数据频谱,两者滤波效果整体一致,频谱曲线呈现吻合状态,所以统一用红色曲线来代替。从滤波效果可以看出,梳状滤波器可以较好地滤除 50 Hz 及其谐波产生的工频干扰,如图中 50、100、150、n50 Hz。由于在工频干扰去除后,数据曲线形态发
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 航空 电磁 时间 数据 干扰 研究 去除
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。