CN105652325B - 基于指数拟合‑自适应卡尔曼的地空电磁数据去噪方法 - Google Patents
基于指数拟合‑自适应卡尔曼的地空电磁数据去噪方法 Download PDFInfo
- Publication number
- CN105652325B CN105652325B CN201610004452.6A CN201610004452A CN105652325B CN 105652325 B CN105652325 B CN 105652325B CN 201610004452 A CN201610004452 A CN 201610004452A CN 105652325 B CN105652325 B CN 105652325B
- Authority
- CN
- China
- Prior art keywords
- data
- air
- electromagnetic
- ground
- fitting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/12—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electromagnetic waves
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electromagnetism (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于指数拟合‑自适应卡尔曼的地空电磁数据去噪方法。本发明针对地空测量的单点电磁数据,根据地空时域电磁数据近似e指数规律衰减的特征,按照对数等间隔确定时窗将数据进行分段,提取每时段内数据的时间常数值作为拟合参数,采用e指数拟合方法对每时段内的数据进行处理,将拟合输出结果作为预测值输入到滤波器中,再应用自适应标量卡尔曼滤波方法对数据进行电磁噪声滤除,将滤波后数据进行电阻率‑深度成像。本发明的滤波方法与现有地空时域电磁数据滤波算法相比,不仅能够有效地抑制地空时域电磁数据中的电磁噪声,而且充分保留和增强了实测数据中的有效信息,提高了地空时域电磁数据的信噪比和质量。
Description
技术领域
本发明涉及一种时域电磁数据滤波方法,适用于符合e指数衰减规律的电磁数据滤波,尤其适用于地空时域电磁数据滤波。
背景技术
地空时域电磁法(Ground-airborne Time-domain electromagnetic system,简称GATEM)是利用地表铺设的长导线源发射双极性脉冲电流激励大地,在空间形成大面积的电磁波,在空中采用无人机或飞艇携带接收线圈对磁场响应进行采集。地空时域电磁法结合了地面瞬变电磁和航空时域电磁的优点。此方法不仅具有比地面瞬变电磁***高效、应对复杂地形能力强等优势,而且比航空时域电磁***空间分辨率更好、探测深度更大、信噪比更高、安全性更高。但是地空时域电磁***在飞行测量过程中,受风向气流、飞行不平稳、传感器姿态等产生运动噪声,以及工频噪声、天电干扰等电磁噪声,导致地空电磁测量信号中含有较为复杂的噪声干扰,因此严重影响后期反演结果的准确性,所以对地空电磁探测数据的去噪方法研究具有重要意义。
目前已经应用到地空电磁数据噪声消除的方法有很多,其中Wang,Y.(2013)设计了一种基于综合小波变换的基线漂移校正方法,可以去除地空电磁测量数据的基线漂移。Bouchedda,A.(2010)提出天电噪声提取和小波系数叠加方法去除天电噪声,也就是应用小波变换能够有效地抑制天电噪声,但是当局部出现较大的天电干扰时,将改变衰减曲线特性。Chen,B.(2014)应用独立主成分分析法对于特征明显的工频强干扰去噪效果较好以及滤除一部分人文噪声。
CN201110197887.4公开了长导线源瞬变电磁地空探测方法,采用长导线源瞬变电磁地空探测装置对深部地质目标体进行探测,获得瞬变电磁的观测数据,将其转换成瞬变电磁虚拟波数据,采用多点数据合成获得瞬变电磁合成孔径数据体,并对合成孔径数据体进行处理解释,获得深部地质目标体信息。
CN201510193706.9公开了基于数学形态学的航空电磁数据去噪方法,提出自适应多尺度复合形态滤波方法,根据信号局部特征和噪声特点自适应的选择结构元素类型和尺寸大小,滤除信号中的正、负脉冲噪声和随机噪声。
CN201110197887.4公开了长导线源瞬变电磁地空探测方法以及处理解释方法,但是并没有给出对地空探测数据的滤波方法;CN201510193706.9公开了一种航空去噪方法,针对航空电磁信号特点进行滤波,虽然可以采用数学形态学这种方法对地空电磁探测数据进行去噪处理,去除信号中的正、负脉冲噪声和随机噪声,但是实际地空电磁探测数据中同时含有白噪声、工频干扰、发动机、运动等电磁噪声,混合存在很难进行分离。目前针对地空电磁探测数据的去噪方法大多针对某一类电磁噪声能够进行有效地去除,但是未考虑地空电磁信号特点,导致有较大噪声干扰时衰减曲线特性被改变。因此需要设计一种考虑地空电磁信号特征并且能够有效滤除多种噪声的滤波方法。
发明内容
本发明的目的在于针对现有去噪方法的不足,根据地空电磁数据衰减特征,提供一种能够同时有效滤除多种电磁噪声,基于指数拟合-自适应卡尔曼的地空电磁数据去噪方法。
本发明的主要思想是:在提取地空电磁数据衰减特征的基础上,建立指数拟合预报方程,将拟合输出结果作为预测值,将拟合输出结果作为预测值输入到滤波器中,滤除实测数据的电磁噪声,提取有效信号,将滤波后地空电磁数据进行电阻率-深度参数计算、成像。
本发明是这样实现的,一种基于指数拟合-自适应卡尔曼的地空电磁数据去噪方法包括:
1)首先从实测连续数据中进行分割,获得单点的地空电磁测量数据;
2)利用加权均值滤波方法对电磁数据进行预处理;
3)提取地空电磁数据的衰减特征,建立指数拟合预报方程,对电磁数据进行e指数拟合;
4)建立自适应标量卡尔曼滤波器,将拟合输出结果作为预测值输入到滤波器中,滤除实测数据的电磁噪声,提取有效信号;
5)将滤波后地空电磁数据进行电阻率-深度参数计算和成像。
进一步地,步骤3中,按照对数等间隔方式将数据进行分段,提取每时段数据的时间常数τ值作为拟合参数,采用e指数拟合方法对每时段数据进行拟合。拟合方程为:
其中,vk为第k道的电磁响应预测值,a1为调整系数,a2为***增益,a3为时间常数τ,tk为第k道的延迟时间。
进一步地,步骤4中包含以下步骤:
Ⅰ、从实测连续数据中,通过计算噪声均值,提取测量***的电磁噪声统计特性;
Ⅱ、输入参数pk,bk和Rk初始值,其中pk为误差协方差,bk遗忘因子,Rk为测量噪声的方差;
Ⅲ、计算加权系数估计值其中vk为步骤3中第k取样道电磁数据的拟合结果;
Ⅳ、计算自适应标量卡尔曼滤波器的基本递归表达式参数;
Ⅴ、计算时变噪声统计估计量;
Ⅵ、对单点地空电磁测量数据的所有时间道重复Ⅲ-Ⅴ步骤,并判断是否完成全部取样道,如果未完成,则计算遗忘因子bk,返回步骤Ⅲ。
步骤4所述的的自适应标量卡尔曼滤波器,其中自适应标量卡尔曼滤波器的基本递归表达式:
kk+1=αkpk(cαkpk+Rk+1)-1 (2)
pk+1=(1-kk+1c)αkpk (3)
其中kk+1为卡尔曼增益,c为信号估计系数,yk+1为第k取样道的地空测量值,为滤波结果。
进一步地,步骤5中包含以下步骤:
a.输入测量点坐标;
b.计算收发距r,确定计算长导线源时所需剖分数N;
c.计算每个电偶极子响应的核函数;
d.计算N个电偶极子的总核函数;
e.对总核函数f(θ)进行最大值归一化处理,其中μ0为真空中磁场导率,σ为地下电阻率,t为时间自变量;
f.实测数据计算的核函数f(θ)是否为最大,确定电阻率值;
g.计算深度参数,进行电阻率-深度成像。
有益效果:本发明与现有技术相比,不仅能够有效地抑制实际地空飞行探测数据中同时含有的白噪声、工频干扰、天电噪声等电磁感应混合噪声,而且充分保留和增强了实测数据中的有效信息,提高了地空时域电磁数据的信噪比和质量,完全满足地空时域电磁要求。本方法为我国特殊地区开展地空飞行电磁探测寻找矿产资源提供新的技术保障,有利于地空时域电磁探测技术的实用化。
附图说明
图1是地空时域电磁探测***示意图;
图2是基于指数拟合-自适应卡尔曼滤波方法流程图;
图3是指数拟合方法流程图;
图4是自适应标量卡尔曼滤波方法流程图;
图5是电磁数据计算电阻率-深度参数方法示意图;
图6是本发明一个实施例的指数拟合效果图;
图7是本发明一个实施例的地空时域电磁数据滤波效果图;
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合吉林省公主岭市怀德镇地空时域电磁探测为例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例
图1是地空时域电磁探测***示意图,在吉林省公主岭市怀德镇地表铺设的长导线源发射双极性脉冲电流,在空中采用滑翔机携带接收线圈对磁场响应进行采集。基于指数拟合-自适应卡尔曼对地空电磁测量数据进行滤波处理。
参见图2结合图1所示,一种基于指数拟合-自适应卡尔曼的地空电磁数据去噪方法,包括:
1)首先从实测连续数据中进行分割,获得单点的地空电磁测量数据;
从滑翔机携带的接收***接收的连续地空电磁测量数据中,提取出接收信号中的衰减曲线,获得单点的地空电磁测量数据及其对应坐标,并对单点的地空电磁测量数据进行以下处理。
2)利用加权均值滤波方法对电磁数据进行预处理;
为了提高拟合精度,对单点的地空电磁测量数据先进行加权均值滤波,公式为:
其中,yk为第k道的地空测量值。
3)提取地空电磁数据的衰减特征,建立指数拟合预报方程,对电磁数据进行e指数拟合;
如图3所示的步骤3中,由于电磁信号在早期衰减快,晚期衰减慢,因此,在进行数据拟合时,按照对数等间隔方式将数据进行分段,提取每时段数据的时间常数τ值作为拟合参数。采用e指数拟合方法对每时段数据进行拟合,拟合方程为:
其中,vk为第k道的电磁响应预测值,a1为调整系数,a2为***增益,a3为时间常数τ,tk为第k道的延迟时间。
4)建立自适应标量卡尔曼滤波器,将拟合输出结果作为预测值输入到滤波器中,滤除实测数据的电磁噪声,提取有效信号;
如图4所示的步骤4中包含以下步骤:
Ⅰ、从实测连续数据中,通过计算噪声均值,提取测量***的电磁噪声统计特性;
地空时域电磁测量过程中,其数学模型可以由如下测量方程给出:
xk=αxk-1+wk-1 (3)
yk=cxk+ηk (4)
其中xk为k时刻的估计信号;yk为k时刻的测量信号;wk-1为测量噪声;ηk为***噪声,α、c为信号估计系数,与信号本身有关;***噪声和测量噪声的统计特性:
E[ηkηl]=Rkδkl (5)
其中Rk为测量噪声的方差,
Ⅱ、根据电性源地空数据特点,以及***噪声和测量噪声的统计特性,给出参数pk,bk和Rk初始值,其中pk为误差协方差,bk遗忘因子。
Ⅲ、计算加权系数估计值其中vk为步骤3中第k取样道电磁数据的拟合结果;
Ⅳ、计算自适应标量卡尔曼滤波器的基本递归表达式参数;
自适应标量卡尔曼滤波器的基本递归表达式:
kk+1=αkpk(cαkpk+Rk+1)-1 (6)
pk+1=(1-kk+1c)αkpk (7)
其中kk+1为卡尔曼增益,yk+1为第k道的地空测量值,为滤波结果;
Ⅴ、计算时变噪声统计估计量;
时变噪声统计估计量通过下式计算:
其中bk为遗忘因子
Ⅵ、对单点地空电磁测量数据的所有时间道重复Ⅲ-Ⅴ步骤,并判断是否完成全部取样道,如果未完成,则计算遗忘因子bk,返回步骤Ⅲ。
5)将滤波后地空电磁数据进行电阻率-深度参数计算、成像。
参见图5所示并结合图2,步骤5中感应电动势计算电阻率方法。长导线源激励响应可以采用多个偶极源响应的叠加等效,根据电偶极子的均匀半空间的解析解,设计了求解全区视电阻率算法。包含以下步骤:
a.输入测量点坐标;
b.计算收发距r,确定计算长导线源时所需剖分数N;
c.计算每个电偶极子响应的核函数;
在地面接收电磁响应时,电偶极源时域感应电动势写为:
式中S为接收线圈有效面积,I为发射电流,t为时间自变量,dl为电偶极子长度,r为收发距,x、y、z为测量点坐标,其中σ为地下电阻率,μ0为真空中磁场导率,erf为高斯误差函数。
d.计算N个电偶极子的总核函数;
长导线源感应电动势响应写为:
式中Hc为长导线源电磁响应,Hdi为第i个电偶极源的电磁响应。把(10)式代入到(11)式中,并化简可得核函数:
式中N为长导线源剖分个数,dli为第i个电偶极子长度,L为长导线源半长度,式中V(t)为测量感应电动势,
e.对总核函数f(θ)进行最大值归一化处理;
f.实测数据计算的核函数f(θ)是否为最大,确定电阻率值;
通过对核函数进行求解θ值来反求视电阻率为了保证解存在,在核函数中引入校正系数,采用核函数的最大值进行归一处理,根据核函数曲线的单调性进行划分为早期和晚期曲线,在每一段曲线上采用二分法进行搜索求解θ值计算电阻率,当求出两个电阻率值,需要和早期、晚期的近似解进行判断,最后确定视电阻率值。
g.计算深度参数,进行电阻率-深度成像。
根据步骤f计算的视电阻率值,计算第i时间道所在的深度
图6为采用图3所示的本发明一个实施例地空理论模型衰减曲线的指数拟合效果图;
图7为本发明一个实施例的地空时域电磁数据滤波后的结果,滤波后的衰减曲线信噪比明显提高,充分验证了指数拟合-自适应卡尔曼滤波方法的有效性,为地空时域电磁测量数据高精度解释提供了新的思路和方法。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于指数拟合-自适应卡尔曼的地空电磁数据去噪方法,其特征在于,包括如下步骤:
1)首先从实测连续数据中进行分割,获得单点的地空电磁测量数据;
2)利用加权均值滤波方法对电磁数据进行预处理;
3)提取地空电磁数据的衰减特征,建立指数拟合预报方程,对电磁数据进行e指数拟合;
4)建立自适应标量卡尔曼滤波器,将拟合输出结果作为预测值输入到滤波器中,滤除实测数据的电磁噪声,提取有效信号;
5)将滤波后地空电磁数据进行电阻率-深度参数计算和成像。
2.按照权利要求1所述的一种基于指数拟合-自适应卡尔曼的地空电磁数据去噪方法,其特征在于:
步骤3中,按照对数等间隔方式将数据进行分段,提取每时段数据的时间常数τ值作为拟合参数,采用e指数拟合方法对每时段数据进行拟合。
3.按照权利要求1所述的一种基于指数拟合-自适应卡尔曼的地空电磁数据去噪方法,其特征在于,步骤4中包含以下步骤:
Ⅰ、从实测连续数据中,通过计算噪声均值,提取测量***的电磁噪声统计特性;
Ⅱ、输入参数pk,bk和Rk初始值,其中pk为误差协方差,bk遗忘因子,Rk为测量噪声的方差;
Ⅲ、计算加权系数估计值其中vk为步骤3中第k取样道电磁数据的拟合结果;
Ⅳ、计算自适应标量卡尔曼滤波器的基本递归表达式参数;
Ⅴ、计算时变噪声统计估计量;
Ⅵ、对单点地空电磁测量数据的所有时间道重复Ⅲ-Ⅴ步骤,并判断是否完成全部取样道,如果未完成,则计算遗忘因子bk,返回步骤Ⅲ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610004452.6A CN105652325B (zh) | 2016-01-05 | 2016-01-05 | 基于指数拟合‑自适应卡尔曼的地空电磁数据去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610004452.6A CN105652325B (zh) | 2016-01-05 | 2016-01-05 | 基于指数拟合‑自适应卡尔曼的地空电磁数据去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105652325A CN105652325A (zh) | 2016-06-08 |
CN105652325B true CN105652325B (zh) | 2017-09-19 |
Family
ID=56491590
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610004452.6A Active CN105652325B (zh) | 2016-01-05 | 2016-01-05 | 基于指数拟合‑自适应卡尔曼的地空电磁数据去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105652325B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106814403B (zh) * | 2017-01-17 | 2019-01-04 | 中国科学院上海微***与信息技术研究所 | 一种补偿瞬变电磁信号负值的方法 |
CN107526111B (zh) * | 2017-08-15 | 2019-03-19 | 中国科学院电子学研究所 | 针对半航空电磁***的天电噪声提取与去除的方法 |
CN109343134A (zh) * | 2018-11-27 | 2019-02-15 | 中煤科工集团西安研究院有限公司 | 一种矿井瞬变电磁探测数据分析解释方法及*** |
CN113552632B (zh) * | 2020-04-23 | 2022-11-01 | 中国石油大学(北京) | 基于小波域卡尔曼滤波的地震互相关信号拾取方法和*** |
CN111780817B (zh) * | 2020-06-07 | 2022-02-11 | 承德石油高等专科学校 | 一种低频励磁电磁流量计噪声信号检测及处理的算法 |
CN111796335B (zh) * | 2020-08-28 | 2023-03-14 | 核工业航测遥感中心 | 航空瞬变电磁时间常数提取方法 |
CN113266335B (zh) * | 2021-03-19 | 2022-06-21 | 吉林大学 | 一种基于最小二乘算法的航空电磁数据***辨识的优化方法 |
CN114186584B (zh) * | 2021-12-06 | 2022-08-30 | 无锡鸣石峻致医疗科技有限公司 | 一种基于扩展指数模型对磁共振信号降噪的方法和设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104567871A (zh) * | 2015-01-12 | 2015-04-29 | 哈尔滨工程大学 | 一种基于地磁梯度张量的四元数卡尔曼滤波姿态估计方法 |
CN104793253A (zh) * | 2015-04-22 | 2015-07-22 | 吉林大学 | 基于数学形态学的航空电磁数据去噪方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
BRPI0511772A (pt) * | 2004-06-01 | 2008-01-08 | Exxonmobil Upstream Res Co | método para uso de uma abordagem de filtro de kalman para o processamento de dados eletromagnéticos |
US8121789B2 (en) * | 2006-07-25 | 2012-02-21 | Exxonmobil Upstream Research Co. | Method for correcting the phase of electromagnetic data |
-
2016
- 2016-01-05 CN CN201610004452.6A patent/CN105652325B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104567871A (zh) * | 2015-01-12 | 2015-04-29 | 哈尔滨工程大学 | 一种基于地磁梯度张量的四元数卡尔曼滤波姿态估计方法 |
CN104793253A (zh) * | 2015-04-22 | 2015-07-22 | 吉林大学 | 基于数学形态学的航空电磁数据去噪方法 |
Non-Patent Citations (2)
Title |
---|
嵌入式***测量中的指数平滑滤波新算法;金玉善 等;《吉林大学学报(工学版)》;20120930;第42卷;第314-316页 * |
电性源时域地空电磁数据小波去噪方法研究;李肃义 等;《地球物理学报》;20130930;第56卷(第9期);第3145-3152页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105652325A (zh) | 2016-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105652325B (zh) | 基于指数拟合‑自适应卡尔曼的地空电磁数据去噪方法 | |
MacLeod et al. | 3-D analytic signal in the interpretation of total magnetic field data at low magnetic latitudes | |
CN103955003B (zh) | 一种在超导瞬变电磁应用中的噪声抑制方法 | |
CN104793253B (zh) | 基于数学形态学的航空电磁数据去噪方法 | |
CN108680966B (zh) | 海洋可控源电磁勘探噪声降噪效果评估方法 | |
CN107329137B (zh) | 基于box样条小波框架的探地雷达地下目标定位方法及*** | |
CN110989002B (zh) | 地空时域电磁***浅部低阻异常体数据解释方法 | |
CN108957403B (zh) | 一种基于广义互相关的高斯拟合包络时延估计方法及*** | |
CN109324315A (zh) | 基于双层次块稀疏性的空时自适应处理雷达杂波抑制方法 | |
CN112818876A (zh) | 一种基于深度卷积神经网络的电磁信号提取与处理方法 | |
CN108169802A (zh) | 一种粗糙介质模型的时域电磁数据慢扩散成像方法 | |
CN116719088B (zh) | 一种航空瞬变电磁数据噪声压制方法 | |
CN103376443B (zh) | 探地雷达探测地面干扰快速消除方法 | |
CN106353819A (zh) | 一种井中三分量微地震初至拾取方法 | |
Fan et al. | An elaborately designed virtual frame to level aeromagnetic data | |
CN103513288B (zh) | 一种二维网格数据的补偿方向滤波方法 | |
Du et al. | Magnetic anomaly detection based on singular spectrum analysis and orthonormal basis functions | |
Wang et al. | A new data processing method for magnetic anomaly detection and localization based on 2-D orthonormal basis functions | |
Zhou et al. | Aeromagnetic anomaly detection under low SNR conditions using multiscale wavelet energy accumulation | |
CN115113163A (zh) | 一种探地雷达多分辨率低秩稀疏分解杂波抑制方法 | |
Xu et al. | Automatic source localization and attenuation of seismic interference noise using density-based clustering method | |
CN102819012B (zh) | 一种基于改进的非负ica的相控阵风廓线雷达信号处理方法 | |
CN112713907A (zh) | 一种基于字典学习的海洋csem降噪方法及*** | |
Li et al. | Research on RTP aeromagnetic gradient data and its applicability in different latitudes | |
CN109709540A (zh) | 一种mimo人体雷达信号通道选择方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |