[发明专利]一种基于三维构造张量的航磁数据地质体边界识别方法

权利要求书

1.一种基于三维构造张量的航磁数据地质体边界识别方法，其特征在于，该方法包括：

获取实测航磁异常数据；

对数据进行预处理；

对预处理后的数据进行空间域计算；

对预处理后的数据进行频率域计算；

利用空间域计算的数据和频率域计算的数据得到水平导数、多阶垂向导数以及构造张量矩阵；

利用张量矩阵得到两种具有不同分辨率的边界识别滤波器NSED₁和边界识别滤波器NSED₂；

对预处理后的数据进行空间域计算包括：

采用空间域求导方式得到异常的水平导数表达式为：

其中，f为实测航磁异常数据，x、y为空间坐标的两个方向,x₁与y₁以及x₂与y₂为不同时间点航磁异常数据对应的空间坐标方向；

对预处理后的数据进行频率域计算包括：

通过频率域求导方式得到垂向导数f为实测航磁异常数据；x、y、z为空间坐标的三个方向，频率域求导先将实测航磁异常数据f做傅氏变换，用它的谱乘以波谱算子r，得到的结果再做反傅氏变换便得到垂向导数

利用空间域和频率域计算方法得到水平导数、垂向导数以及构造张量矩阵：

根据实测或计算获得的梯度值定义构造张量矩阵T，表达式为：

采用边界识别滤波器NSED₁进行航磁数据的边界识别，经验参数k选择为0.0001，NSED₁表达式为：

采用边界识别滤波器NSED₂进行航磁数据的边界识别，进一步提升边界识别结果的分辨率，通过试验确定参数k为0.0001，NSED₂表达式为：