[发明专利]一种增强重磁场边界识别方法

说明书

本发明属于重磁数据处理技术领域，具体涉及一种增强重磁场边界识别方法，包括：步骤一：输入重磁异常测量数据；步骤二：网格化重磁数据，并计算总水平导数异常；步骤三：选取滑动窗口的大小，计算各个小子域的平均值和均方差；步骤四：选取均方差较小的前三个区域，计算三个小子域的总水平导数的均方差；步骤五：输出均方差最小子域的均值；步骤六：求取滑动小子域滤波结果的垂向导数和总水平导数；步骤七：求取滑动小子域滤波结果的倾斜角，输出结果。

技术领域

本发明属于重磁数据处理技术领域，具体涉及一种增强重磁场边界识别方法。

背景技术

重磁异常是地下物质密度和磁性分布不均匀的综合反映，有横向分辨率高的优势，因此利用重磁位场在推断地质目标体的边界位置时有独特的优势。常用的方法有总水平导数法、垂向导数法、解析信号振幅、倾斜角法等，都是利用异常极大值的位置或零值的位置确定目标体的边界，但是此类方法会降低边缘识别的精度，对噪声的压制能力也较弱。国内目前还没有涉及滑动小子域滤波法与倾斜角法的有效结合来识别重磁位场的边界。

因此，需要设计一种增强重磁场边界识别方法，将滑动小子域滤波法与倾斜角法的有效结合用于有效地识别地质体的边界。

发明内容

本发明针对上述现有技术的不足，提供一种增强重磁场边界识别方法，用于解决现有技术中异常极大值位置或零值位置确定目标体的边界时边缘识别精度降低，噪声压制能力较弱的技术问题。

本发明的技术方案：

一种增强重磁场边界识别方法，包括以下步骤：

步骤一：输入重磁异常测量数据；

步骤二：网格化重磁数据，并计算总水平导数异常；

步骤三：选取滑动窗口的大小，计算各个小子域的平均值和均方差；

步骤四：选取均方差较小的前三个区域，计算三个小子域的总水平导数的均方差；

步骤五：输出均方差最小子域的均值；

步骤六：求取滑动小子域滤波结果的垂向导数和总水平导数；

步骤七：求取滑动小子域滤波结果的倾斜角，输出结果。

所述步骤二包括：利用克里金法将测量的重磁数据进行网格化处理，生成规则的网格数据，磁法数据先进行化极处理，然后根据公式(1)求取重磁异常数据的总水平导数异常：

式中，f代表重磁异常，f_TD代表重磁异常的总水平导数异常。

所述步骤三还包括：两次滑动窗口，大滑动窗口的大小为n，大窗口节点数据5n，内部小子域窗口的大小为k＝(n+1)/2，然后小子域的窗口从左到右逐点滑动，直到到达大窗口的最末端；根据公式(2)、(3)计算各个小子域的均方差作为判别准则；

式中，f代表重磁异常，表示小子域的平均值，δ代表小子域的均方差。

所述步骤四还包括：在布设的k²个小子域内选取数据均方差较小的前3个小子域作为目标区；根据公式(2)和公式(4)计算选定的3个小子域总水平导数的均值和均方差；

式中，f_TD代表重磁异常的总水平导数异常，代表小子域内总水平导数的平均值，δ_TD代表小子域内总水平导数的均方差。

所述步骤五还包括：在选定的3个小子域中选取水平导数均方差最小子域的均值作为最后的输出结果；大窗口从左到右滑动到下一点直至完成全区的计算，得到滑动小子域滤波的结果。