[发明专利]航磁数据向下延拓方法

权利要求书

1.一种航磁数据向下延拓方法，其特征在于，包括：

对观测面的航磁垂向导数进行向上延拓，得到分布在观测面以上，七个不同高度的航磁垂向导数；

利用观测面的航磁数据、航磁垂向导数以及其向上延拓七个高度的航磁垂向导数，构成向下延拓的八阶Adams-Bashforth格式；

利用八阶Adams-Bashforth格式计算航磁数据的向下延拓初步结果；

利用积分迭代方式，对向下延拓初步结果进行迭代，计算向下延拓最终结果；

利用积分迭代方式，对向下延拓初步结果进行迭代，得到向下延拓最终结果，包括：

将所得到的八阶Adams-Bashforth格式向下延拓的结果作为高度-h处的迭代初始值u(x,y,z₀-h)₀，利用向上延拓方法计算高度为z₀的观测面的航磁数据u(x,y,z₀)₀，比较迭代前后航磁数据u(x,y,z₀)和u(x,y,z₀)₀的差异Δu₀，将差异Δu₀代入向下延拓值u(x,y,z₀-h)₀中，反复此过程，直到迭代前后航磁数据的差异Δu₀满足设计要求为止。

2.根据权利要求1所述的航磁数据向下延拓方法，其特征在于，八阶Adams-Bashforth格式由如下公式表示：

其中，u_z(x,y,z₀+h)、u_z(x,y,z₀+2h)、u_z(x,y,z₀+3h)、u_z(x,y,z₀+4h)、u_z(x,y,z₀+5h)、u_z(x,y,z₀+6h)、u_z(x,y,z₀+7h)为z₀高度以上h、2h、3h、4h、5h、6h、7h七个不同高度的航磁垂向导数，h表示延拓高度，z₀表示观测面的z轴坐标。

3.根据权利要求1或2所述的航磁数据向下延拓方法，其特征在于，还包括：

在对观测面的航磁垂向导数进行向上延拓，得到分布在观测面以上，七个不同高度的航磁垂向导数之前，获取观测面的航磁垂向导数。

4.根据权利要求3所述的航磁数据向下延拓方法，其特征在于，获取观测面的航磁垂向导数，包括：

判断是否具有观测面的航磁垂向导数的实测数据；

若有，获取观测面的航磁垂向导数；

若无，利用ISVD方法计算观测面的航磁垂向导数。

5.根据权利要求4所述的航磁数据向下延拓方法，其特征在于，利用ISVD方法计算观测面的航磁垂向导数，包括：

利用波数域的航磁数据换算方法计算标量位；

利用有限差分法计算标量位水平方向的二阶导数；

利用拉普拉斯方程，得到标量位的二阶垂直方向导数，即位场的一阶垂直方向导数。