[发明专利]一种地磁日变数据补缺方法及处理终端

说明书

本发明涉及一种地磁日变数据补缺方法及处理终端，所述方法包括如下步骤：步骤1：获取待补缺地磁日变数据，以及至少一台参考地磁台站的参考地磁台站数据；步骤2：对待补缺地磁日变数据中为磁平静日的所有地磁日变数据进行数据分离，得到不同日期的静日数据，若为磁扰日，则对参考地磁台站数据进行数据分离，得到磁扰数据；通过相关分析法和聚类分析法处理，得到各个日期的静日数据之间的相关性和样本静日数据与基准静日数据之间的聚类程度，根据聚类程度和相关性强度得到第一静日数据，从而得到最终静日数据，完成地磁日变数据补缺。本发明具有补缺效果非常良好；磁平静日情况下，能够对单独的地磁台站数据自身不齐全进行地补缺。

技术领域

本发明涉及地磁数据处理技术领域，具体是一种地磁日变数据补缺方法及处理终端。

背景技术

日变改正是地磁测量数据处理中的一个重要环节，其处理后的数据质量直接影响日变改正精度，进而影响最终的地磁测量成果。对于日变改正，通常选择工区(地磁测量作业区域)附近的固定地磁台站或野外测量期间而在工区内临时投放的观测地磁台站来获得地磁日变数据，但受制于测量环境等因素，经常会出现采集的地磁日变数据部分缺失，导致部分测量数据无法进行日变改正，而根据目前的行业规范要求，对于无法进行日变改正的测线需要作废处理，直接影响地磁测量的工作效率，甚至导致某次地磁测量失败，无法提交成果。

目前，对于日变数据补缺大多数是利用多个地磁台站进行推算得到，例如，单汝俭等在1990年发表的文章《局部地区地磁日变及拟合方法研究》中提出了三种局部地区地磁日变的拟合方法，分别为二维多项式最小二乘拟合法、时空拟合法和线性内插法；边刚等在2009年发表的文章《海洋磁力测量中多站地磁日变改正值计算方法》中采用多站同步实例分析了加权平均法和函数拟合法在海洋磁测中的应用，其忽略经度效应的影响，采用纬度方向的距离来计算某点的日变改正值；卞光浪等在2010年发表的文章《基于纬差加权法的海洋磁力测量多站地磁日变改正值计算》中提出了基于纬差加权法的海洋磁力测量多站地磁日变改正值计算，其摆脱了距离上的相对关系，而将日变数据的差异只建立在纬度差关系上，文中通过对沿岸几个日变站的数据进行计算验证，证实该方法较距离加权法效果更好；彭飞等在2015年发表的文章《调和分析方法在海洋磁力测量日变改正中的应用》中分析了调和分析法，该方法是基于傅里叶级数建立了日变数据处理谐波分析模型，将磁静和磁扰数据进行了合理分离，运用该方法可辅助地磁日变数据的推算。

以上的地磁日变数据补缺都是利用某一个或某几个地磁台站数据作为样本地磁台站数据进行推算得到某一位置处的地磁日变数据，上述方法能够实现的前提是样本地磁台站数据齐全，而对于单一的地磁台站缺失部分数据而需要进行补缺的情况却没有相应的解决方法。实际工作中，经常会遇到获取某一单独的地磁台站数据自身不齐全，导致在缺少部分数据的情况下，无法对单一地磁台站缺失数据进行补缺，直接影响日变改正工作的开展，进而影响地磁测量数据的成果交付。

发明内容

针对现有技术的不足，本发明的目的之一提供一种地磁日变数据补缺方法，其能够解决地磁日变数据补缺的问题；

本发明的目的之二提供一种处理终端，其能够解决地磁日变数据补缺的问题。

实现本发明的目的之一的技术方案为：一种地磁日变数据补缺方法，包括如下步骤：

步骤1：获取除缺失数据时间段之外的待补缺地磁台站的待补缺地磁日变数据，以及若所述待补缺地磁台站的缺失数据时间段为磁扰日，则获取至少一台参考地磁台站的参考地磁台站数据；

步骤2：对待补缺地磁日变数据中为磁平静日的所有地磁日变数据进行数据分离，得到不同日期的静日数据，若所述缺失数据时间段为磁扰日，则还包括对参考地磁台站数据进行数据分离，得到磁扰数据；

步骤3：对任意两日之间的静日数据通过相关分析法进行处理，得到任意两日之间的静日数据之间的相关性，根据各个日期的静日数据之间的相关性强度，筛选出相关性排序靠前的若干日期的静日数据，作为样本静日数据，