CN109856689A - 一种超导航磁梯度张量数据抑噪处理方法和*** - Google Patents
一种超导航磁梯度张量数据抑噪处理方法和*** Download PDFInfo
- Publication number
- CN109856689A CN109856689A CN201910150000.2A CN201910150000A CN109856689A CN 109856689 A CN109856689 A CN 109856689A CN 201910150000 A CN201910150000 A CN 201910150000A CN 109856689 A CN109856689 A CN 109856689A
- Authority
- CN
- China
- Prior art keywords
- data
- magnetic field
- field value
- magnetic gradient
- reference background
- 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.)
- Granted
Links
Landscapes
- Measuring Magnetic Variables (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开一种超导航磁梯度张量数据抑噪处理方法及***。方法包括:获取在飞行测量平台坐标系下的大地背景参考磁场值、大地背景参考磁梯度场值和飞行测量线上的背景参考磁场值的时间变化率;根据获得的数据对三分量磁强计数据进行补偿校正处理,得到磁场校正补偿处理数据;根据磁场校正补偿处理数据、大地背景参考磁场值和飞行测量线上的背景参考磁场值的时间变化率对磁梯度计单一通道测量数据进行补偿,得到磁梯度数据补偿模型;根据磁梯度数据补偿模型,得到补偿后的磁梯度数据;对补偿后的磁梯度数据进行低通滤波处理,得到处理后的磁梯度数据;根据处理后的磁梯度数据对补偿效果进行评估。本发明能够提高实际测量数据的补偿有效性。
Description
技术领域
本发明涉及地球物理航空磁法勘探领域,特别是涉及一种超导航磁梯度张量数据抑噪处理方法和***。
背景技术
航空磁法勘探是一种重要并广泛应用于矿产和环境勘测的地球物理勘探技术手段。近年来,得益于基于超导量子干涉仪(SQUID)的磁梯度测量***的大力发展,已使得直接测量高精度的磁梯度数据成为可能,而超导全张量磁梯度仪因其高信息量、高磁场灵敏度、体积小等诸多优点,被认为是第三代航磁探测的发展方向。目前,我国已建设有基于SQUID的全张量磁梯度测量平台原型机,不同于传统的采用机载方式进行测量的航磁梯度测量飞行平台,该原型机将磁梯度测量仪及相关设备(GPS、惯导***(INS)、读入读出电路等)装载于一个吊舱内,并通过直升机拖曳的方式进行飞行测量,该方式将大大降低直升机对磁测量干扰的影响,甚至在拖曳缆绳足够长的条件下直升机影响可以忽略不计,但吊舱内仍然存在大量会产生磁干扰的设备,因此有效地通过磁补偿处理提高飞行平台的测量精度,以实现高分辨率磁梯度测量具有重要的意义。
目前针对该问题的方法还十分有限,除了通过硬件设备压制干扰的方法外,已有的软补偿方法多数仅考虑传统干扰源对测量的影响,但在实际飞行测量过程中往往噪声***更为复杂,仅通过一种模型及相关的单一流程很难将复杂的各类干扰充分滤除。
发明内容
本发明的目的是提供一种超导航磁梯度张量数据抑噪处理方法和***,能够提高实际测量数据的补偿有效性。
为实现上述目的,本发明提供了如下方案:
一种超导航磁梯度张量数据抑噪处理方法,包括:
获取在飞行测量平台坐标系下的大地背景参考磁场值、大地背景参考磁梯度场值和飞行测量线上的背景参考磁场值的时间变化率;
根据所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对三分量磁强计数据进行补偿校正处理,得到磁场校正补偿处理数据;
根据所述磁场校正补偿处理数据、所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对磁梯度计单一通道测量数据进行补偿,得到磁梯度数据补偿模型;
根据所述磁梯度数据补偿模型,得到补偿后的磁梯度数据;
对所述补偿后的磁梯度数据进行低通滤波处理,得到处理后的磁梯度数据;
根据所述处理后的磁梯度数据对补偿效果进行评估。
可选的,所述获取在飞行测量平台坐标系下的大地背景参考磁场值、大地背景参考磁梯度场值和飞行测量线上的背景参考磁场值的时间变化率,具体包括:
获取在大地坐标系下的大地背景参考磁场值和大地背景参考磁梯度场值;
将所述大地坐标系下的大地背景参考磁场值和大地背景参考磁梯度场值转换到飞行测量平台坐标系中,得到在飞行测量平台坐标系下的大地背景参考磁场值和大地背景参考磁梯度场值;
根据所述飞行测量平台坐标系下的大地背景参考磁场值,通过数值差分或傅里叶变换的方式,得到飞行测量线上的背景参考磁场值的时间变化率。
可选的,所述根据所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对三分量磁强计数据进行补偿校正处理,得到磁场校正补偿处理数据,具体包括:
通过正则化最小二乘拟合或Huber范数拟合的方式求解方程得到第一校正补偿系数c1;其中,A为直接获取的观测数据;
对测量数据Bm进行第一阶段校正补偿计算,得到初始处理数据
通过正则化最小二乘拟合或Huber范数拟合的方式求解方程得到第二校正补偿系数c2;
对数据B1m进行第二阶段补偿计算,通过公式得到磁场校正补偿处理数据。
可选的,所述根据所述磁场校正补偿处理数据、所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对磁梯度计单一通道测量数据进行补偿,得到磁梯度数据补偿模型,具体包括:
根据所述磁场校正补偿处理数据、所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对磁梯度计单一通道测量数据进行补偿,得到磁梯度数据补偿模型
其中,l为通道编号;为通道测量的磁梯度值;为通道l的补偿磁梯度值;He=(Hex,Hey,Hez)为参考背景磁场,dHe=(dHex,dHey,dHez)为参考背景磁场值He的时间变化率,Br=(Brx,Bry,Brz)为测量磁场校正补偿处理后数据,t为惯导数据采样时间,fs为对应采样频率,p,q为干扰源信号滤波基长度,相应的,A为由各类观测数据及相关计算量构成的N×M阶矩阵,x=(s1,o1,a1,x,a2,x,K,a2p+1,z,b1,x,K,b2p+1,z,c1,x,K,c2q+1,z)T为实际计算中使用的补偿参数,其中s1=1/sl,o1=ol/sl,sl为测量通道的尺度偏差;ol为测量通道的偏差系数。
可选的,所述根据所述磁梯度数据补偿模型,得到补偿后的磁梯度数据,具体包括:
根据所述磁梯度数据补偿模型求解补偿系数;
根据所述补偿系数,确定补偿后的磁梯度数据。
可选的,所述对所述补偿后的磁梯度数据进行低通滤波处理,得到处理后的磁梯度数据,具体包括:
对所述补偿后的磁梯度数据采用巴特沃斯低通滤波方法进行低通滤波处理,得到处理后的磁梯度数据。
可选的,所述根据所述处理后的磁梯度数据对补偿效果进行评估,具体包括:
根据所述处理后的磁梯度数据,采用标准化参数对补偿效果进行评估,所述标准化参数包括均方根、标准差和改善比。
一种超导航磁梯度张量数据抑噪处理***,包括:
获取模块,用于获取在飞行测量平台坐标系下的大地背景参考磁场值、大地背景参考磁梯度场值和飞行测量线上的背景参考磁场值的时间变化率;
磁场校正补偿处理数据确定模块,用于根据所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对三分量磁强计数据进行补偿校正处理,得到磁场校正补偿处理数据;
磁梯度数据补偿模型确定模块,用于根据所述磁场校正补偿处理数据、所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对磁梯度计单一通道测量数据进行补偿,得到磁梯度数据补偿模型;
补偿磁梯度数据确定模块,用于根据所述磁梯度数据补偿模型,得到补偿后的磁梯度数据;
磁梯度数据处理模块,用于对所述补偿后的磁梯度数据进行低通滤波处理,得到处理后的磁梯度数据;
补偿效果评估模块,用于根据所述处理后的磁梯度数据对补偿效果进行评估。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供一种超导航磁梯度张量数据抑噪处理方法,获取在飞行测量平台坐标系下的大地背景参考磁场值、大地背景参考磁梯度场值和飞行测量线上的背景参考磁场值的时间变化率;根据所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对三分量磁强计数据进行补偿校正处理,得到磁场校正补偿处理数据;根据所述磁场校正补偿处理数据、所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对磁梯度计单一通道测量数据进行补偿,得到磁梯度数据补偿模型;根据所述磁梯度数据补偿模型,得到补偿后的磁梯度数据;对所述补偿后的磁梯度数据进行低通滤波处理,得到处理后的磁梯度数据;根据所述处理后的磁梯度数据对补偿效果进行评估。本发明能够有效的解决多种测量平台***噪声和环境噪声的干扰,提高了SQUID航磁梯度仪的磁梯度数据补偿的精度和有效性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例超导航磁梯度张量数据抑噪处理方法流程图;
图2为在飞行测量坐标系下对三分量磁强计测量值的校正补偿结果示意图,其中,(a)表示Bx分量校正前、校正后与相应参考场数据的对比结果;(b)表示By分量校正前、校正后与相应参考场数据的对比结果;(c)表示Bz分量校正前、校正后与相应参考场数据的对比结果;(d)表示磁场模量|B|校正前、校正后与相应参考场数据的对比结果;
图3为磁场三分量数据补偿校正前后,各磁场分量Fourier频率幅值谱分析的对比图,其中,(a)表示Bx分量校正前、校正后与相应参考场数据的对比结果;(b)表示By分量校正前、校正后与相应参考场数据的对比结果;(c)表示Bz分量校正前、校正后与相应参考场数据的对比结果;
图4为在飞行测量坐标系下对G1通道的磁梯度计测量值的补偿及滤波处理结果示意图,其中,(a)表示补偿前与补偿后数据的对比结果;(b)表示进行补偿处理后,滤波前与滤波后数据的对比结果;
图5为对G1通道的磁梯度计测量值进行补偿及滤波处理后,测量数据、补偿后数据、滤波后数据的Fourier频率幅值谱分析对比图;
图6为本发明实施例超导航磁梯度张量数据抑噪处理***结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种超导航磁梯度张量数据抑噪处理方法和***,能够提高实际测量数据的补偿有效性。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例超导航磁梯度张量数据抑噪处理方法流程图。如图1所示,一种超导航磁梯度张量数据抑噪处理方法,包括:
获取在飞行测量平台坐标系下的大地背景参考磁场值、大地背景参考磁梯度场值和飞行测量线上的背景参考磁场值的时间变化率;
步骤101:根据所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对三分量磁强计数据进行补偿校正处理,得到磁场校正补偿处理数据;
具体包括:
获取在大地坐标系下的大地背景参考磁场值和大地背景参考磁梯度场值;
将所述大地坐标系下的大地背景参考磁场值和大地背景参考磁梯度场值转换到飞行测量平台坐标系中,得到在飞行测量平台坐标系下的大地背景参考磁场值和大地背景参考磁梯度场值;
根据所述飞行测量平台坐标系下的大地背景参考磁场值,通过数值差分或傅里叶变换的方式,得到飞行测量线上的背景参考磁场值的时间变化率。
首先通过飞行测量平台搭载的GPS***获取各测量点处的经纬高数据(经度,纬度,高度)、通过全球地磁场参考模型IGRF-12以及区域地磁台网观测数据,可以计算出大地坐标系下(NEU),大地背景参考磁场值He0=(He0x,He0y,He0z)T、大地背景参考磁梯度场值Hge0=(H0xx,H0xy,H0xz,H0yy,H0yz,H0zz)。
然后,通过INS***获取各测量点处的测量平台的飞行姿态数据(偏转角η、俯仰角λ、翻滚角ψ),可将大地坐标系下的He0及Hge0转换到飞行测量平台坐标系下。其中,转换后的背景参考磁场值He及参考磁梯度场值Hge可以分别表示为:
He=(Hex,Hey,Hez)T=T2He (1)
其中,T2为大地坐标系到飞行测量坐标系的坐标转换矩阵,即:
由于在实际中受飞行测量设备的限制,使得GPS、INS的数据采样率与磁强计Bm、磁梯度计Gi的采样率不一致(本实施例中,前者采样率要远低于后者),因此造成相应计算得到的参考磁场数据He及Hge的分辨率要低于磁传感器测量信号Bm及Gi的分辨率。对此,在插值精度允许的条件下,利用三次样条插值对He及Hge进行数据加密,同时对Bm及Gi进行下采样处理,使得原来两个采样率不同的数据集经过处理后获得同样长度的数据样本,以便用于下面磁校正及补偿处理。为了简便表达式,经过插值或下采样处理的数据用原有同样的符号表示。
最后,将He看作时间序列He(t),通过数值差分或FFT变换的方式,计算在飞行测量线上的背景磁场值He的时间变化率dHe=(dHx,dHy,dHz)=(He(t+1/fs)-He(t))·fs,其中fs为数据采样频率。
步骤102:根据所述磁场校正补偿处理数据、所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对磁梯度计单一通道测量数据进行补偿,得到磁梯度数据补偿模型;
磁梯度张量数据测量仪是由6个不共面的一阶SQUID平面磁梯度计和3轴正交参考SQUID磁强计组成的9通道信号输出设备,由于超导芯片制备工艺的限制使得平面梯度计的两个线圈具有不一致性,因此每一个梯度计信号中还会包含寄生的三分量磁场信号。虽然多梯度计模块的几何构建可以降低自身固有的不平衡性,但该精度仍然无法满足在地磁背景场下条件的实际测量需求,因此需要在梯度计附近安装三轴磁强计,同时测量磁场三分量场,用于测量的磁梯度场数据的平衡度补偿处理。由于三轴磁强计与磁梯度计的测量环境相同,同样会受到多种干扰噪声的影响,因此在用于磁梯度计补偿前需要对其本身采集的磁场数据进行校正补偿处理。
首先建立基于飞行测量平台坐标系下的磁场三分量数据校正补偿综合模型。考虑磁测传感器的低频噪声干扰主要来源于两个部分:内部噪声(多个传感器敏感度差异性、传感器偏差、安装误差及软磁性干扰)及外部干扰(测量平台的剩磁、感磁、涡流磁场干扰),结合三轴磁强计校正模型及Tolles-Lawson磁补偿模型,建立如下的三分量磁场数据校正补偿综合模型:
Bm=SDB0+O1+O2+M1He+M2dHe+n (4)
其中,B0=(Bx,By,Bz)为在飞行测量坐标系下的真实磁场三分量值,在高空测量条件下应与He近似相等;
S=diag(sx,sy,sz)为尺度偏差因子;
为安装误差矩阵,偏差角axy为实际测量x’轴在直角坐标系xoy平面上的投影与x轴的小夹角,其他偏差角定义类同;
O1=diag(o1x,o1y,o1z)为传感器偏差系数矩阵;
O2=diag(o2x,o2y,o2z)为测量平台的剩磁干扰系数矩阵
M1=(m1ij)3×3为测量平台的感磁干扰系数矩阵;
M2=(m2ij)3×3为测量平台的涡流干扰系数矩阵;
n为传感器高斯白噪声。
因此,上述模型可以整理得到如下磁场校正补偿系数求解方程:
其中,为实际计算使用的校正补偿系数,其各个部分与原始具有物理含义的噪声干扰系数的关系为:F(Bm,1,dHe)为由数据Bm,dHe,1形成的方程系数矩阵,当采样数据量为N时,该矩阵维度为3N×21。
该校正补偿系数求解方程可以通过正则化最小二乘拟合或Huber范数拟合的方式进行求解,其求解精度由系数矩阵F的条件数决定。由理论模型测试分析可知,由Bm和dHe构成的向量组成线性相关,因此在存在数据噪声的情况下直接求解该方程会造成严重的计算不稳定,因此实际计算中将校正补偿系数c分为两组及分别进行求解,并对磁场数据进行分阶段校正补偿处理,具体方法如下所述:
通过正则化最小二乘拟合或Huber范数拟合的方式求解方程得到第一校正补偿系数c1;其中,A为直接获取的观测数据;
对测量数据Bm进行第一阶段校正补偿计算,得到初始处理数据
通过正则化最小二乘拟合或Huber范数拟合的方式求解方程得到第二校正补偿系数c2;
对数据B1m进行第二阶段补偿计算,通过公式得到磁场校正补偿处理数据。
步骤103:根据所述磁梯度数据补偿模型,得到补偿后的磁梯度数据;
首先,由于是对磁梯度计单一通道测量数据进行补偿,因此需要将位于飞行测量坐标系下的参考背景磁梯度场Hge坐标旋转到与待补偿通道平行的方向,设坐标旋转矩阵为T3,在实施例中,六个平面磁梯度计成六棱台分布,每一片在安装支架平面圆的投影成夹角为60度的均匀分布,具体的转化矩阵形式如下所示:
其中,分别为与6个平面磁梯度计通道方向平行的参考背景磁梯度场值,α为平面梯度计倾斜面与安装支架平面的夹角(该采样设备倾角值α=63°)。
然后,建立针对单一磁梯度计通道测量数据的综合补偿模型:
根据平面磁梯度计的工艺说明及建立三分量磁场校正补偿综合模型的原理可知,磁梯度计中单一通道的低频噪声干扰主要来源于三个部分:内部偏差(多个传感器敏感度差异性、传感器偏差、安装误差及软磁性干扰)、外部干扰(测量平台的剩磁、感磁、涡流磁场)以及寄生磁场。
综合考虑,外部干扰源等价的磁矩矢量M与数据He,dHe的关系可以表示为:
其中,Mr,Mind,Mv分别表示外部干扰源等价的剩磁、感磁、涡流磁矩矢量;k=(kx,ky,kz)为感磁系数,e=(ex,ey,ez)为涡流系数,其分别与干扰源的磁性结构及几何形态分布有关。
而磁矩矢量M与磁梯度张量分量Bij的关系可以表示为如下形式,也同样为一种线性关系:
其中,Gij=(Gijx,Gijy,Gijz)对应于磁梯度分量的核函数向量,其具体公式为:
(1)单位球形源:
(2)正六面体源:
rp(xp,yp,zp)为测量点位置,rQ(xQ,yQ,zQ)为干扰源位置;ΩQ为干扰源的几何形状,如将该区域干扰源离散近似为一系列等价球形偶极子源,则外部干扰源产生的磁梯度场可以近似表示为如下形式:
其中,N表示离散的磁球个数,V表示离散的磁球体积。
结合上述的飞行测量坐标系到测量通道平面的转换公式T3(公式(6)),以及外部干扰源的磁矩模型(公式(7)),可以将外部干扰源沿单一通道方向的磁梯度测量干扰场表示为如下形式:
其中,以及g0,g1是由坐标转换变换T3 l、磁梯度核函数Gij、离散近似磁球体积V、以及感磁系数k,涡流系数e等常值参数构成的线性算子。经过一系列等价转换后,初始的复杂干扰场方程被转化为一个关于向量1,He,dHe的简单线性组合形式,这种表示可以理解为外部干扰场空间的局部可以表示为如下一个线性空间:
Vout={Gout|Gout∈span{1,He,dHe}}。 (13)
因此,基于类似信号线性空间理论,类比建立了磁梯度测量的寄生磁场空间:
Vbalance={Gbalance|Gbalance∈span{Br}}。 (14)
进一步,考虑数据采集中的***误差及数据预处理过程造成的数值误差,比如不同采集设备、同一采集设备不同采集通道之间的信号延迟,以及步骤101中对于采样频率不一致的数据采用插值或下采样的处理方式,均会造成不同采样数据间无法做到完全同步。因此,如若干扰场的基信号仅选取当前时间点的一个值进行线性估计则很可能造成不可忽略的偏差量。为了更好的处理这一问题,本发明借鉴了线性滤波的思想对基本模型进行了优化:选取基信号会首先利用临近点信息对当前点的基信号进行线性估计,然后再利用估计量作为基向量对干扰场进行线性重构,此时需要将采集信号看作时间序列。据此,修正后的外部干扰场空间和寄生磁场空间可以表示为:
其中,2p+1,2q+1为两组基向量的估计域长度。综上所述,结合采集通道信号的内部偏差校正量,具体实施中采用如下模型对单一通道采集信号进行表示:
其中,l为通道编号;为通道测量的磁梯度值;Gl为通道测量点处的真实磁梯度值;sl为测量通道的尺度偏差;为测量通道的***偏差,为外部干扰源的剩磁干扰,两者由于在一定时间内对应单一通道测量值均为恒定不变的常值,可综合为一个偏差系数ol;为外部干扰源感磁干扰的线性算子,由(2p+1)×3个系数组成;为外部干扰源涡流干扰的线性算子,由(2p+1)×3个系数组成;h0为寄生磁场干扰的线性算子,由(2q+1)×3个系数组成;n为高频Gauss白噪声。通过对上述采集信号模型进行整理,可以得到用于磁梯度数据补偿的计算模型,如下所示:
其中,A为由各类观测数据及相关计算量构成的N×M阶矩阵,x=(s1,o1,a1,x,a2,x,K,a2p+1,z,b1,x,K,b2p+1,z,c1,x,K,c2q+1,z)T为实际计算中使用的补偿参数,共M=12p+6q+11个;需要注意的是,在实际计算中,由于补偿系数具有一定的物理含义,如sl,因此需要根据实际传感器性能定义求解值域约束条件,如实施例中设约束区间sl∈[0.5,1]。
根据上述磁梯度数据综合补偿模型,建立如下的约束优化问题求解补偿系数:
(1)对于测量区域内梯度值很小的情况下,采用带约束的线性最小二乘数据拟合的方法求解补偿系数x,求解算法可采用信赖域反射算法(Trust region reflective),有效集算法(Active set)或内点法(Interiorpoint)等。实施例中采用的是常用的有效集算法求解如下问题:
(2)对于测量区域内存在局部较大梯度值的情况下,为了保证求解的稳定性,采用Huber数据拟合的方式求解补偿系数x,求解算法采用交替方向乘子方法(ADMM):
通过如上计算方案求得补偿系数x之后,使用如下公式计算信号通道l采集的磁梯度补偿值:
其中,l为通道编号;为通道测量的磁梯度值;为通道l的补偿磁梯度值;He=(Hex,Hey,Hez)为参考背景磁场,dHe=(dHex,dHey,dHez)为参考背景磁场值He的时间变化率,Br=(Brx,Bry,Brz)为测量磁场校正补偿处理后数据,t为惯导数据采样时间,fs为对应采样频率,p,q为干扰源信号滤波基长度,相应的,A为由各类观测数据及相关计算量构成的N×M阶矩阵,x=(s1,o1,a1,x,a2,x,K,a2p+1,z,b1,x,K,b2p+1,z,c1,x,K,c2q+1,z)T为实际计算中使用的补偿参数,其中s1=1/sl,o1=ol/sl,sl为测量通道的尺度偏差;ol为测量通道的偏差系数。
综上所述,新建立的磁梯度数据补偿综合模型不仅可以有效的滤除飞行测量平台在实际航磁测量时可能存在的多种类型的干扰噪声,而且基于信号线性空间的思想,该模型也有利于进一步的扩展和修正,只要能分析得到与干扰信号频带相同的特征信号即可用于扩展磁测干扰场空间域,通过求解相关干扰场的滤波补偿系数,消除干扰噪声而对测量数据进行有效补偿。
步骤104:对所述补偿后的磁梯度数据进行低通滤波处理,得到处理后的磁梯度数据;
通过上述一系列的数据补偿处理后,补偿后的数据还存在一部分在仪器输出条件下无法滤除的高频噪声,因此需要根据实际应用问题的信号频带特征,以及已知的传感器芯片的本底噪声频带特征,对补偿后的数据进一步进行低通滤波处理。在本发明的实施例中,已知航空测量的地磁场的有效信号频率应小于10Hz,因此最后通过巴特沃斯低通滤波方法对补偿数据进行滤波处理,获得最终的磁梯度数据处理结果。
步骤105:根据所述处理后的磁梯度数据对补偿效果进行评估。
完成上述磁补偿及数据滤波处理后,还需要对补偿质量进行评价。实施例中对补偿后数据使用如下标准化参数对补偿结果进行评价:均方根RMS、标准差σc、改善比IR。参数值的具体定义如下所示:
其中,dc为补偿及滤波处理后的数据,d0为原始的采集数据,N为同一类型的数据个数。
图6为本发明实施例超导航磁梯度张量数据抑噪处理***结构图。如图6所示,一种超导航磁梯度张量数据抑噪处理***,包括:
获取模块201,用于获取在飞行测量平台坐标系下的大地背景参考磁场值、大地背景参考磁梯度场值和飞行测量线上的背景参考磁场值的时间变化率;
磁场校正补偿处理数据确定模块202,用于根据所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对三分量磁强计数据进行补偿校正处理,得到磁场校正补偿处理数据;
磁梯度数据补偿模型确定模块203,用于根据所述磁场校正补偿处理数据、所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对磁梯度计单一通道测量数据进行补偿,得到磁梯度数据补偿模型;
补偿磁梯度数据确定模块204,用于根据所述磁梯度数据补偿模型,得到补偿后的磁梯度数据;
磁梯度数据处理模块205,用于对所述补偿后的磁梯度数据进行低通滤波处理,得到处理后的磁梯度数据;
补偿效果评估模块206,用于根据所述处理后的磁梯度数据对补偿效果进行评估。
具体实施例1:
下面通过实测数据对本发明所提出的基于SQUID磁梯度计单通道测量数据补偿方法的效果进行验证和说明。本发明中涉及的SQUID磁梯度测量仪共有9个磁传感器输出通道,ch0-ch2为三分量磁强计输出通道,ch3-ch8为磁梯度计输出通道,本实施例使用ch0-ch3通道测量的数据对磁梯度数据(ch3)补偿效果进行说明。
图2为在飞行测量坐标系下对三分量磁强计测量值的校正补偿结果示意图,其中,(a)表示Bx分量校正前、校正后与相应参考场数据的对比结果;(b)表示By分量校正前、校正后与相应参考场数据的对比结果;(c)表示Bz分量校正前、校正后与相应参考场数据的对比结果;(d)表示磁场模量|B|校正前、校正后与相应参考场数据的对比结果;图3为磁场三分量数据补偿校正前后,各磁场分量Fourier频率幅值谱分析的对比图,其中,(a)表示Bx分量校正前、校正后与相应参考场数据的对比结果;(b)表示By分量校正前、校正后与相应参考场数据的对比结果;(c)表示Bz分量校正前、校正后与相应参考场数据的对比结果;图4为在飞行测量坐标系下对G1通道的磁梯度计测量值的补偿及滤波处理结果示意图,其中,(a)表示补偿前与补偿后数据的对比结果;(b)表示进行补偿处理后,滤波前与滤波后数据的对比结果;图5为对G1通道的磁梯度计测量值进行补偿及滤波处理后,测量数据、补偿后数据、滤波后数据的Fourier频率幅值谱分析对比图。
相应的,补偿和滤波处理结果对应的质量评价指标由表1给出,结合图4、图5和表1可知,经过补偿和滤波处理后的磁梯度数据均方根小于20pT,改善比IR可达2.36e3,因此基于本发明提出的方法能够获得良好的补偿结果,可达到仪器分辨率的要求。
表1基于综合磁梯度补偿方法的处理数据评价参数统计
至此,已经结合附图对本发明的基本方法和实施例进行了详细描述。依据以上描述,本领域技术人员应当对本发明的基于SQUID航磁梯度仪单一采集通道信号的磁梯度数据补偿方法有了清楚的认识。需要说明的是,在附图或说明书正文中,未绘示或描述的实现方式,均为所属技术领域中普通技术人员所知的形式,并未进行详细说明。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的***而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种超导航磁梯度张量数据抑噪处理方法,其特征在于,包括:
获取在飞行测量平台坐标系下的大地背景参考磁场值、大地背景参考磁梯度场值和飞行测量线上的背景参考磁场值的时间变化率;
根据所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对三分量磁强计数据进行补偿校正处理,得到磁场校正补偿处理数据;
根据所述磁场校正补偿处理数据、所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对磁梯度计单一通道测量数据进行补偿,得到磁梯度数据补偿模型;
根据所述磁梯度数据补偿模型,得到补偿后的磁梯度数据;
对所述补偿后的磁梯度数据进行低通滤波处理,得到处理后的磁梯度数据;
根据所述处理后的磁梯度数据对补偿效果进行评估。
2.根据权利要求1所述的超导航磁梯度张量数据抑噪处理方法,其特征在于,所述获取在飞行测量平台坐标系下的大地背景参考磁场值、大地背景参考磁梯度场值和飞行测量线上的背景参考磁场值的时间变化率,具体包括:
获取在大地坐标系下的大地背景参考磁场值和大地背景参考磁梯度场值;
将所述大地坐标系下的大地背景参考磁场值和大地背景参考磁梯度场值转换到飞行测量平台坐标系中,得到在飞行测量平台坐标系下的大地背景参考磁场值和大地背景参考磁梯度场值;
根据所述飞行测量平台坐标系下的大地背景参考磁场值,通过数值差分或傅里叶变换的方式,得到飞行测量线上的背景参考磁场值的时间变化率。
3.根据权利要求1所述的超导航磁梯度张量数据抑噪处理方法,其特征在于,所述根据所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对三分量磁强计数据进行补偿校正处理,得到磁场校正补偿处理数据,具体包括:
通过正则化最小二乘拟合或Huber范数拟合的方式求解方程得到第一校正补偿系数c1;其中,A为直接获取的观测数据;
对测量数据Bm进行第一阶段校正补偿计算,得到初始处理数据
通过正则化最小二乘拟合或Huber范数拟合的方式求解方程得到第二校正补偿系数c2;
对数据B1m进行第二阶段补偿计算,通过公式得到磁场校正补偿处理数据。
4.根据权利要求1所述的超导航磁梯度张量数据抑噪处理方法,其特征在于,所述根据所述磁场校正补偿处理数据、所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对磁梯度计单一通道测量数据进行补偿,得到磁梯度数据补偿模型,具体包括:
根据所述磁场校正补偿处理数据、所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对磁梯度计单一通道测量数据进行补偿,得到磁梯度数据补偿模型
其中,l为通道编号;为通道测量的磁梯度值;为通道l的补偿磁梯度值;He=(Hex,Hey,Hez)为参考背景磁场,dHe=(dHex,dHey,dHez)为参考背景磁场值He的时间变化率,Br=(Brx,Bry,Brz)为测量磁场校正补偿处理后数据,t为惯导数据采样时间,fs为对应采样频率,p,q为干扰源信号滤波基长度,相应的,A为由各类观测数据及相关计算量构成的N×M阶矩阵,x=(s1,o1,a1,x,a2,x,K,a2p+1,z,b1,x,K,b2p+1,z,c1,x,K,c2q+1,z)T为实际计算中使用的补偿参数,其中s1=1/sl,o1=ol/sl,sl为测量通道的尺度偏差;ol为测量通道的偏差系数。
5.根据权利要求1所述的超导航磁梯度张量数据抑噪处理方法,其特征在于,所述根据所述磁梯度数据补偿模型,得到补偿后的磁梯度数据,具体包括:
根据所述磁梯度数据补偿模型求解补偿系数;
根据所述补偿系数,确定补偿后的磁梯度数据。
6.根据权利要求1所述的超导航磁梯度张量数据抑噪处理方法,其特征在于,所述对所述补偿后的磁梯度数据进行低通滤波处理,得到处理后的磁梯度数据,具体包括:
对所述补偿后的磁梯度数据采用巴特沃斯低通滤波方法进行低通滤波处理,得到处理后的磁梯度数据。
7.根据权利要求1所述的超导航磁梯度张量数据抑噪处理方法,其特征在于,所述根据所述处理后的磁梯度数据对补偿效果进行评估,具体包括:
根据所述处理后的磁梯度数据,采用标准化参数对补偿效果进行评估,所述标准化参数包括均方根、标准差和改善比。
8.一种超导航磁梯度张量数据抑噪处理***,其特征在于,包括:
获取模块,用于获取在飞行测量平台坐标系下的大地背景参考磁场值、大地背景参考磁梯度场值和飞行测量线上的背景参考磁场值的时间变化率;
磁场校正补偿处理数据确定模块,用于根据所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对三分量磁强计数据进行补偿校正处理,得到磁场校正补偿处理数据;
磁梯度数据补偿模型确定模块,用于根据所述磁场校正补偿处理数据、所述大地背景参考磁场值和所述飞行测量线上的背景参考磁场值的时间变化率对磁梯度计单一通道测量数据进行补偿,得到磁梯度数据补偿模型;
补偿磁梯度数据确定模块,用于根据所述磁梯度数据补偿模型,得到补偿后的磁梯度数据;
磁梯度数据处理模块,用于对所述补偿后的磁梯度数据进行低通滤波处理,得到处理后的磁梯度数据;
补偿效果评估模块,用于根据所述处理后的磁梯度数据对补偿效果进行评估。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910150000.2A CN109856689B (zh) | 2019-02-28 | 2019-02-28 | 一种超导航磁梯度张量数据抑噪处理方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910150000.2A CN109856689B (zh) | 2019-02-28 | 2019-02-28 | 一种超导航磁梯度张量数据抑噪处理方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109856689A true CN109856689A (zh) | 2019-06-07 |
CN109856689B CN109856689B (zh) | 2020-12-29 |
Family
ID=66899362
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910150000.2A Expired - Fee Related CN109856689B (zh) | 2019-02-28 | 2019-02-28 | 一种超导航磁梯度张量数据抑噪处理方法和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109856689B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110146839A (zh) * | 2019-05-30 | 2019-08-20 | 中国海洋大学 | 一种移动平台磁梯度张量***校正方法 |
CN110596619A (zh) * | 2019-09-16 | 2019-12-20 | 中国科学院上海微***与信息技术研究所 | 一种全张量磁梯度测量组件及其优化方法 |
CN111079285A (zh) * | 2019-12-16 | 2020-04-28 | 山东大学 | 一种全张量磁梯度数据补偿优化方法及*** |
CN112946760A (zh) * | 2021-02-02 | 2021-06-11 | 清华大学 | 基于正则化方法的未爆弹三维立体成像方法、装置及*** |
CN113671582A (zh) * | 2021-08-26 | 2021-11-19 | 吉林大学 | 基于三分量squid的电性源感应-极化效应探测方法 |
CN113739821A (zh) * | 2021-08-31 | 2021-12-03 | 北京航空航天大学 | 一种基于pid算法的原子自旋陀螺仪全自动磁补偿方法 |
CN117607974A (zh) * | 2023-11-15 | 2024-02-27 | 中国科学院空天信息创新研究院 | 一种数据补偿方法、装置、电子设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110013481A1 (en) * | 2007-11-12 | 2011-01-20 | David Alan Clark | Method and apparatus for detecting marine deposits |
CN105509737A (zh) * | 2015-11-26 | 2016-04-20 | 哈尔滨工业大学 | 一种不受地磁变化影响的航空运动平台磁干扰补偿方法 |
CN106997035A (zh) * | 2017-03-29 | 2017-08-01 | 吉林大学 | 一种基于磁梯度不变量的磁梯度计校正方法 |
-
2019
- 2019-02-28 CN CN201910150000.2A patent/CN109856689B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110013481A1 (en) * | 2007-11-12 | 2011-01-20 | David Alan Clark | Method and apparatus for detecting marine deposits |
CN105509737A (zh) * | 2015-11-26 | 2016-04-20 | 哈尔滨工业大学 | 一种不受地磁变化影响的航空运动平台磁干扰补偿方法 |
CN106997035A (zh) * | 2017-03-29 | 2017-08-01 | 吉林大学 | 一种基于磁梯度不变量的磁梯度计校正方法 |
Non-Patent Citations (1)
Title |
---|
SHUANGXI JI 等: "Compensation of SQUID gradient data in magnetic tensor gradiometer system for airborne exploration", 《THE INTERNATIONAL SYMPOSIUM ON DEEP EARTH EXPLORATION AND PRACTICES》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110146839B (zh) * | 2019-05-30 | 2021-11-16 | 中国海洋大学 | 一种移动平台磁梯度张量***校正方法 |
CN110146839A (zh) * | 2019-05-30 | 2019-08-20 | 中国海洋大学 | 一种移动平台磁梯度张量***校正方法 |
CN110596619A (zh) * | 2019-09-16 | 2019-12-20 | 中国科学院上海微***与信息技术研究所 | 一种全张量磁梯度测量组件及其优化方法 |
CN110596619B (zh) * | 2019-09-16 | 2021-07-09 | 中国科学院上海微***与信息技术研究所 | 一种全张量磁梯度测量组件及其优化方法 |
CN111079285B (zh) * | 2019-12-16 | 2022-02-08 | 山东大学 | 一种全张量磁梯度数据补偿优化方法及*** |
CN111079285A (zh) * | 2019-12-16 | 2020-04-28 | 山东大学 | 一种全张量磁梯度数据补偿优化方法及*** |
CN112946760A (zh) * | 2021-02-02 | 2021-06-11 | 清华大学 | 基于正则化方法的未爆弹三维立体成像方法、装置及*** |
CN112946760B (zh) * | 2021-02-02 | 2024-02-06 | 清华大学 | 基于正则化方法的未爆弹三维立体成像方法、装置及*** |
CN113671582A (zh) * | 2021-08-26 | 2021-11-19 | 吉林大学 | 基于三分量squid的电性源感应-极化效应探测方法 |
CN113671582B (zh) * | 2021-08-26 | 2023-08-22 | 吉林大学 | 基于三分量squid的电性源感应-极化效应探测方法 |
CN113739821B (zh) * | 2021-08-31 | 2022-06-17 | 北京航空航天大学 | 一种基于pid算法的原子自旋陀螺仪全自动磁补偿方法 |
CN113739821A (zh) * | 2021-08-31 | 2021-12-03 | 北京航空航天大学 | 一种基于pid算法的原子自旋陀螺仪全自动磁补偿方法 |
CN117607974A (zh) * | 2023-11-15 | 2024-02-27 | 中国科学院空天信息创新研究院 | 一种数据补偿方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109856689B (zh) | 2020-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109856689A (zh) | 一种超导航磁梯度张量数据抑噪处理方法和*** | |
CN105891755B (zh) | 航空器挂载式磁通门磁梯度张量仪的校正方法 | |
CN107544042B (zh) | 一种磁力计阵列校正方法 | |
CN104345348B (zh) | 航空超导全张量磁梯度测量***相关参数获取装置及方法 | |
CN111433634A (zh) | 基于航磁补偿误差模型的磁补偿方法 | |
Munschy et al. | Scalar, vector, tensor magnetic anomalies: Measurement or computation? | |
Schiffler et al. | Calibration of SQUID vector magnetometers in full tensor gradiometry systems | |
AU2008231589B2 (en) | Terrain correction systems | |
CN106997035A (zh) | 一种基于磁梯度不变量的磁梯度计校正方法 | |
CN110133544B (zh) | 航空超导全张量磁补偿系数的获取方法、终端及存储介质 | |
CN109814163A (zh) | 一种基于扩展补偿模型的航磁张量数据抑噪方法及*** | |
CN107132587B (zh) | 航空超导全张量磁梯度测量***安装误差标定方法及装置 | |
CN109856690A (zh) | 基于混合范数拟合的航磁梯度张量数据抑噪方法及*** | |
Shockley | Ground vehicle navigation using magnetic field variation | |
CN113156355B (zh) | 一种超导全张量磁梯度测量装置的磁干扰补偿方法 | |
CN103562753A (zh) | 重力梯度计勘测技术 | |
CN110274586A (zh) | 包含多光系原子磁力仪方向误差补偿的航空磁补偿方法 | |
Getscher et al. | Magnetic gradient tensor framework for attitude-free position estimation | |
Chen et al. | An improved geomagnetic navigation method based on two-component gradient weighting | |
Yin et al. | Calculation of the magnetic gradient tensor from total magnetic anomaly field based on regularized method in frequency domain | |
Shen et al. | Geometry structure optimization of hexagonal pyramidal full tensor magnetic gradient probe | |
Wooldridge | Review of modern airborne gravity focusing on results from GT-1A surveys | |
CN110596619B (zh) | 一种全张量磁梯度测量组件及其优化方法 | |
Zhang et al. | A component compensation method for magnetic interferential field | |
Li et al. | Compensation method for the carrier magnetic interference of underwater magnetic vector measurement system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201229 |
|
CF01 | Termination of patent right due to non-payment of annual fee |