论文部分内容阅读
时间域航空电磁法是以飞机为载体,通过测量地下介质感应产生的随时间变化的二次场,来探测地下介质的电性分布的一种地球物理勘探方法。具有勘探面积广、成本相对较低、受地形地貌影响较小、勘探效率较高以及具有良好的水平分辨率等优点,因此,可以被广泛的应用在众多领域中,比如:地下水分布调查和矿物资源勘探、未爆炸物的检测、地质填图以及地质灾害监测等。时间域航空电磁法通常采集的是信号较弱的二次场响应,其有效信号频带较宽,从而使得接收到的信号容易受到多种噪声的影响,制约了航空电磁探测系统对地下介质的勘探深度以及后期的数据处理和反演成像精度。所以,对时间域航空电磁数据进行预处理,提高数据的信噪比,对提高数据处理质量和反演解释精度有着重要的意义。常规的抑制噪声手段大多是从单一噪声入手,从时域、频域、信号分解的角度进行滤波,但航空电磁信号与噪声在时频域存在重叠,常规处理手段往往达不到理想效果。由于许多现象或信号可以看作是随机信号,因此,从统计学角度对信号进行处理研究是目前重要手段之一。本文在“十三五”国家重点研发计划项目课题“直升机航空电磁探测数据解释软件系统研发”支持下,利用高斯过程回归方法对时间域航空电磁数据进行去噪研究,并对其有效性和可行性进行分析。具体来说,主要研究内容如下:(1)时间域航空电磁模拟数据的生成。首先,使用一维正演模拟方法获得基于中心回线装置的一维航空瞬变电磁响应数据。然后,建立层状二维地电模型,并利用2.5维正演模拟方法生成模拟的二维电磁剖面数据。最后,根据实测数据中噪声的高斯特性和天电噪声的特点,在一、二维电磁数据的基础上添加高斯白噪声和天电噪声得到模拟后的含噪数据,从而为后面的去噪方法研究及去噪结果分析提供了基础测试数据。(2)基于高斯过程回归模型的去噪方法。首先,针对高斯过程回归中核函数的选择问题对常见的平方指数协方差函数和有理二次协方差函数定义的高斯先验做采集实验,分析哪种核函数优越性高为后续去噪实验做模型支撑;然后针对高斯过程回归中超参数选择问题,从单一参数(长度规模尺度)入手,做正态线性拟合实验,分析不同参数值对模型拟合效果有何影响,为后续去噪实验自适应获取超参数值提供理论依据。最后,利用高斯过程回归方法对时间域航空电磁一维模拟数据进行去噪,并对其有效性和可行性进行分析。(3)基于稀疏高斯过程回归的去噪方法。由于高斯过程回归方法中需要对协方差函数进行求逆,导致了传统高斯过程回归方法需要消耗大量的计算成本,因此本文采用稀疏高斯过程回归方法降低运算成本,提高预测效率。首先,对正弦模型添加一定噪声做预测实验,证明稀疏伪输入法的去噪可行性。然后,在固定伪输入集个数的情况下基于不同训练点数做曲线拟合实验,分析训练点数对预测结果和预测速度的影响。最后,在固定训练点数的情况下针对不同伪输入集个数将稀疏伪输入法和传统高斯过程回归方法应用到时间域航空电磁一维模拟数据去噪中,并通过信噪比和均方根误差评价两者的去噪效果,通过训练时间和预测时间评价两方法的去噪效率。(4)高斯过程回归和稀疏高斯过程回归去噪方法性能比较。首先利用高斯过程回归和稀疏高斯过程回归两种方法对模拟二维剖面电磁响应数据进行去噪;然后对去噪后的数据进行反演成像,分析上述两种方法的去噪效果;最后使用这两种方法对新疆某探区实测数据进行处理,进一步验证所推荐去噪方法的可行性。