【技术实现步骤摘要】
多尺度电磁场分量去噪方法、装置、设备及存储介质
[0001]本申请涉及地球物理勘探
,特别地涉及一种多尺度电磁场分量去噪方法、装置、设备及存储介质。
技术介绍
[0002]大地电磁测深观测的天然电磁场,具有信号弱,频带宽等特点,极易受到各种噪声的干扰。依据噪声来源,MT信号中的噪声主要包括场源噪声、人文噪声、输入端噪声和静噪声四类。人文噪声是最为普遍出现、干扰最为严重、噪声特性最为复杂的一类。含噪声的MT数据信噪比低,最终反演得到的阻抗、电阻率计算结果误差明显,从而使反演结果差。
[0003]提高电磁测量数据的信噪比,对视电阻率曲线计算及后期资料解释有着至关重要的意义。经过几十年的发展,形成了诸如Robust处理、远参考技术、小波变换、Hilbert
‑
Huang变换等行之有效的噪声分析、处理方法。目前,大地电磁各类去噪方法自有其优势,但也各有不足。互功率谱法和加权平均互功率谱法,能够得到有效压制、去除不相关噪声的信号。但是,由于接收序列中不相关噪声和相关噪声同时存在,且记录时间有限。将信号进 ...
【技术保护点】
【技术特征摘要】
1.一种多尺度电磁场分量去噪方法,其特征在于,包括:S1:利用多项式函数方式,对原始时间序列两端部分数据链进行外扩,得到外扩数据;S2:应用平稳分量的预处理方法,提取所述外扩数据中的平稳信号,得到平稳数据;S3:采用自适应小波阈值处理方式,压制所述平稳数据中的工频谐波干扰,得到去公频数据;S4:在时间域对所述去公频数据开展EMD分解,对高频分量和低频分量分别去噪,然后将去噪后的高频分量、低频分量和残余模量通过线性重构得到重构去噪的信号。2.根据权利要求1所述的方法,其特征在于,所述利用多项式函数方式,对原始时间序列两端部分数据链进行外扩的具体公式为:其中,g(k)为原始信号,g(j)表示扩边后信号,a,b为扩展系数,j为信号分解多项式;a=[a0…
a
m
];b=[b0…
b
m
],M1为信号左扩边项数,k为原始信号项数,n0为信号右扩边项数。3.根据权利要求1所述的方法,其特征在于,所述应用平稳分量的预处理方法,提取所述外扩数据中的平稳信号的具体公式为:其中Y(t)为数据量为N的时间序列,Y'(t)为提取平稳分量信号,Y(t,0:M
‑
h)、Y(t,h/2:M
‑
h/2)、Y(t,h:M)为截取不同长度的时间序列,n为扩边数据,M=N+2n,t为时间,h为滑动窗,N为时间序列长度即信号的数据量。4.根据权利要求1所述的方法,其特征在于,所述采用自适应小波阈值处理方式,压制所述平稳数据中的工频谐波干扰之前,所述方法还包括:将时间域信号变换到频率域,对分别对实部、虚部进行WT变换。5.根据权利要求1所述的方法,其特征在于,所述自适应小波阈值的选取方法包括:利用固定阈值和无偏似然估计软阈值相组合,获得最优变量阈值,适...
【专利技术属性】
技术研发人员:齐文杰,李弘,徐凤姣,肖鹏飞,
申请(专利权)人:中国石油化工股份有限公司石油物探技术研究院,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。