[发明专利]一种AMT去噪方法、系统及计算机可读存储介质

说明书

本发明涉及一种AMT去噪方法、系统及计算机可读存储介质，属于地球物理电磁法勘探领域。通过设定可根据实际数据情况做变更的结构元素参数并选择三角型、抛物线型、圆盘型和余弦型四种结构元素构建广义形态滤波器，分别对高频、中频、低频文件中包含的Ex、Ey、Hx、Hy等电场分量做组合广义形态滤波处理，提取到各序列的噪声形态，再将之从原序列中删除，最后重构去噪后的高频、中频、低频文件。上述方法能够去除时间序列中的多种形态类噪声。通过计算机可读存储介质上存储的计算机程序，能直接生成去噪后的TS2、TS3、TS4文件，操作方便，可灵活使用。

技术领域

本发明涉及地球物理电磁法勘探领域，特别是涉及一种AMT去噪方法、系统及计算机可读存储介质。

背景技术

音频大地电磁测深(AMT)法已在油气远景区推断、地热勘查、寻找水源和良导性矿体、探测沉积盆地构造等领域取得了广泛的应用。野外利用V8(主机)结合AMTC30(磁探头)的方式开展数据采集，获取到TS2(高频)、TS3(中频)、TS4(低频)三个时间序列文件，对应的采样率分别为24000、2400、150Hz。每种TS文件都包含了Ex、Ey、Hx、Hy、Hz(E表示电场，H表示磁场)这5种电磁场分量时间序列。

但由于AMT法以微弱的天然大地电磁场为场源来探测地下电性结构，并且随着国家工业的发展，采集数据的工作区受多种噪声干扰的情况愈发严重，导致采集到的信号被淹没在噪声中，数据质量差。

因此，亟需开发一种可针对不同噪声进行去噪的方法、系统及计算机可读存储介质。通过该方法可去除AMT数据中的不同形态类噪声，从而得到准确的地下电性结构信息。

发明内容

本发明的目的是提供一种AMT去噪方法、系统及计算机可读存储介质，从而得到去噪后的高频、中频、低频时间序列文件。

为实现上述目的，本发明提供了如下方案：

一种AMT去噪方法，所述方法包括：

采用AMT法获取高频、中频、低频时间序列文件，所述高频、中频、低频时间序列文件为地下电性结构信息的文件；

选择三角型、抛物线型、圆盘型和余弦型四种结构元素构建广义形态滤波器；

根据所述高频时间序列文件的数据设定结构元素参数K2、L2，所述结构元素参数K2、L2根据所述高频时间序列文件的数据不同而不同，其中，K2为极值参数，L2为零点参数；

根据所述广义形态滤波器及结构元素参数K2、L2，分别对所述高频时间序列文件中包含的Ex、Ey、Hx、Hy电磁场分量时间序列做组合广义形态滤波处理，获得所述高频时间序列文件中包含的所述Ex、Ey、Hx、Hy各自的噪声序列，所述Ex和Ey分别为电场在x轴和y轴的分量，所述Hx和Hy分别为磁场在x轴和y轴的分量，所述组合广义滤波处理包含负结构元素的广义形态滤波处理和正结构元素的广义形态滤波处理；

分别从所述高频时间序列文件中包含的所述Ex、Ey、Hx、Hy中删除各自的噪声序列，获得去噪后所述高频时间序列文件的Ex序列、Ey序列、Hx序列和Hy序列；

利用所述去噪后的所述高频时间序列文件的Ex序列、Ey序列、Hx序列和Hy序列重构去噪后的高频文件；

根据所述中频时间序列文件的数据设定结构元素参数K3、L3，所述结构元素参数K3、L3根据所述中频时间序列文件的数据不同而不同；

根据所述广义形态滤波器及结构元素参数K3、L3，分别对所述中频时间序列文件中包含的Ex、Ey、Hx、Hy电磁场分量时间序列做组合广义形态滤波处理，获得所述中频时间序列文件中包含的所述Ex、Ey、Hx、Hy各自的噪声序列；