[发明专利]一种基于位场数据向下延拓的场源深度获得方法

权利要求书

1.一种基于位场数据向下延拓的场源深度获得方法，其特征在于，包括以下步骤：

S1：获取实测位场数据，所述实测位场数据包括目标体重力场或磁力场数据；

S2:依据所获得实测位场数据计算其解析信号数据；

S3：根据实测位场数据及解析信号数据，调用预设的深度计算模型，获得目标体的深度位置，其中所述深度计算模型根据改进的Chebyshev-Padé逼近向下延拓计算方法，对所述解析信号数据进行向下延拓计算，并对计算结果进行简单归一化计算，基于归一化结果寻找局部极大值的位置，该局部极大值的位置即为场源的深度位置；

S4：对步骤S3得到的归一化结果三维可视化，直观的显示场源的深度位置；

所述步骤S2中解析信号数据的计算方法为：

其中，f为位场函数，和分别为x-、y-和z-方向的一阶导数；所述步骤S3中基于改进的Chebyshev-Padé逼近的向下延拓的计算方法包括以下步骤：

a、建立频率域Chebyshev-Padé逼近的向下延拓因子；

b、依据频率域向上延拓的低通滤波特征，建立适合深度计算的低通滤波因子；

c、对频率域Chebyshev-Padé逼近的向下延拓因子乘以低通滤波因子，获得改进的向下延拓算子，对解析信号数据进行向下延拓计算；

d、对向下延拓计算结果进行归一化计算；

所述步骤a中Chebyshev-Padé向下延拓因子为：

其中Δh为延拓距离，当Δh为正数，则为向下延拓，Δh为负数，则为向上延拓，(ω_x，ω_y)是x-和y-方向的波数，其中为圆波数；

所述步骤b中向上延拓计算因子为：

其中(ω_x，ω_y)是x-和y-方向的波数，为圆波数，Δh为延拓距离；

根据其滤波特征，适合深度计算的低通滤波因子为：

其中(ω_x，ω_y)是x-和y-方向的波数，为圆波数，Δh为延拓距离；α为滤波系数，当场源为三度体时取α＝0.8，当场源为二度体时取α＝0.95；所述步骤c中视向下延拓计算因子为：

其中(ω_x，ω_y)是x-和y-方向的波数，为圆波数，Δh为延拓距离；可得视向下延拓计算结果为：

其中为解析信号频率域值，ASM_h(ω_x，ω_y，Δh)为向下延拓后值，F^-1[]为傅里叶逆变换；

所述步骤d中归一化计算方法为：

其中，ln(||)为求取解析信号延拓结果的绝对值并计算自然对数。