[发明专利]一种重磁位场边界弱信息提取方法

权利要求书

1.一种重磁位场边界弱信息提取方法，其特征在于，包括：

获取重磁异常测量数据；

对所述重磁异常测量数据进行网格化和化极处理得到网格化的重磁异常数据；

使用滑动小子域滤波法对所述网格化的重磁异常数据进行处理得到滑动小子域滤波结果；

根据所述滑动小子域滤波结果确定地质体边界的位置。

2.根据权利要求1所述的一种重磁位场边界弱信息提取方法，其特征在于，所述使用滑动小子域滤波法对所述网格化的重磁异常数据进行处理得到滑动小子域滤波结果，包括：

在所述网格化重磁异常数据上构建滑动窗口，并计算各个小子域窗口内的小子域重磁异常数据的平均值和均方差；

按照从小到大的顺序排列所述均方差，并选取前三个均方差所对应的小子域窗口作为目标区；

计算各个目标区的总水平导数的均方差，并选取最小的总水平导数均方差所对应的小子域窗口输出的小子域重磁异常数据的平均值作为所述滑动窗口的输出结果；

不断滑动所述滑动窗口直到完成整个网格化的重磁异常数据的去噪得到滑动小子域滤波结果。

3.根据权利要求2所述的一种重磁位场边界弱信息提取方法，其特征在于，所述在所述网格化重磁异常数据上构建滑动窗口，并计算各个小子域窗口内的小子域重磁异常数据的平均值和均方差，包括：

采用公式：

计算各个小子域窗口内的小子域重磁异常数据的平均值和均方差；其中，表示小子域的平均值，k表示小子域窗口的大小，δ表示小子域的均方差。

4.根据权利要求3所述的一种重磁位场边界弱信息提取方法，其特征在于，目标区的总水平导数的均方差计算公式为：

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

5.根据权利要求4所述的一种重磁位场边界弱信息提取方法，其特征在于，所述根据所述滑动小子域滤波结果确定地质体边界的位置，包括：

求取所述滑动小子域滤波结果的总水平导数异常与总梯度模量异常；

根据所述总水平导数异常和所述总梯度模量异常得到滑动小子域滤波结果的Theta图；

根据所述Theta图确定地质体边界的位置。

6.根据权利要求5所述的一种重磁位场边界弱信息提取方法，其特征在于，根据所述总水平导数异常和所述总梯度模量异常得到滑动小子域滤波结果的Theta图，包括：

采用公式：

得到滑动小子域滤波结果的Theta图；其中，MSSf_TD为滑动小子域滤波结果的总水平导数异常，MSSf_ASM为滑动小子域滤波结果的总梯度模量异常，MSSf为滑动小子域滤波结果，MSSTheta为Theta图计算公式，f为重磁异常数据。