基于灰色系统的含异常值航磁数据处理方法及电子设备

技术特征：
1.一种基于灰色系统的含异常值航磁数据处理方法，其特征在于，包括以下步骤：s1：获取原始磁测数据的序列，表示为b0＝(b
01
,b
02
,
…
,b
0i
,
…
,b
0n
)，n表示序列中磁测数据的个数，i表示磁测数据的序号，i＝1,2,...,n；s2：将序列b0从小到大进行排序得到新的序列b
asc
；s3：对序列b
asc
进行逐位累加处理得到数据累加曲线b
acc
；s4：计算磁测数据的平均值其平均值累加曲线为s5：构建磁测数据的灰色区域图：将平均值累加曲线b
ava
作为灰色区域图的上界函数，将平均值累加曲线与数据累加曲线的差值曲线b
inf
作为灰色区域图的下界函数，并且磁测数据的中间点i
p
为差值曲线b
inf
的转折点其中h为灰色识别系数，δb
max
为数据平均值累加序列与数据累加序列的差值δb的最大值；s6：基于所述灰色区域图进行磁测数据异常值剔除，具体流程为：s6.1：当b
inf
(i)≤b
acc
(i)≤b
ava
(i)时，磁测数据全部包含在灰色区域图的内部，即不存在异常值，其中b
inf
(i)、b
acc
(i)、b
ava
(i)分别表示差值序列、数据累加序列和平均值累加序列；s6.2：当b
acc
(i)≤b
inf
(i)或b
ava
(i)≤b
acc
(i)时，存在灰色区域图以外的数据，即存在异常值，灰色区域图两端会产生数据累加序列与差值序列的交叉点，通过分析交叉点两端的数据信息，在原始磁测数据中将两端数据进行剔除处理；s7：将异常值的数据剔除后，返回重复步骤s2-s6，直到所有磁测数据均满足b
inf
(i)≤b
acc
(i)≤b
ava
(i)，即实现了磁测数据中的全部异常值剔除。2.根据权利要求1所述的含异常值航磁数据处理方法，其特征在于，步骤s6.2包括：s6.2.1：若b
acc
(i)≤b
inf
(i)，i＜i
p
，说明b
acc
前i个数据小于正常磁测值的最小界限，且b
asc
(1)数值最小，将b
asc
(1)判定为异常值；s6.2.2：若b
ava
(i)≤b
acc
(i)，i≥i
p
，说明b
acc
后n-i个数据大于正常磁测值的最大界限，因下界函数b
inf
与数据累加曲线b
acc
汇于最后一个测量点，b
ava
(n)＝b
acc
(n)，无法判断是否存在异常值，所以除b
acc
(n)外，b
ava
(n-1)值最大，将b
asc
(n-1)判定为异常值。3.根据权利要求1所述的含异常值航磁数据处理方法，其特征在于，当磁测数据中不存在异常值时，认为序列b
asc
是等差序列，并且表示为b
asc
(i)＝b
asc
(1) r0·
i；其中，b
asc
(1)为序列中的第一个数，则b
asc
(i)写为：此时，则b
ava
(i)表示为：
数据平均值累加序列与数据累加序列的差值δb为：所以，当中间点时，δb达到最大δb
max
。4.根据权利要求3所述的含异常值航磁数据处理方法，其特征在于，所述下界函数b
inf
表示为：式中，代入上式得5.根据权利要求1所述的含异常值航磁数据处理方法，其特征在于，所述灰色识别系数h的取值为5。6.根据权利要求1所述的含异常值航磁数据处理方法，其特征在于，在步骤s7之后，还包括以下步骤：基于异常值剔除后的磁测数据，利用最小二乘法获取磁干扰系数，进行航磁补偿。7.一种电子设备，其特征在于，所述电子设备包括存储器和处理器，所述存储器上存储有计算机程序，该计算机程序被所述处理器执行时，执行如权利要求1-6任一项所述的含异常值航磁数据处理方法。

技术总结
本发明公开了一种基于灰色系统的含异常值航磁数据处理方法及电子设备，该方法包括：获取原始磁测数据的序列B0；将序列B0从小到大排序得到新序列B


技术研发人员：葛健 吴涛 胡祥云 董浩斌 王科 杨文君 王光旭 徐为 肖敏 唐悦能
受保护的技术使用者：中国地质大学（武汉）
技术研发日：2022.08.15
技术公布日：2022/12/12