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

说明书

本发明提供了一种航磁数据向下延拓方法。该技术包括：对观测面的航磁垂向导数进行向上延拓，得到分布在观测面以上，七个不同高度的航磁垂向导数；利用观测面的航磁数据、航磁垂向导数以及其向上延拓七个高度的航磁垂向导数，构成向下延拓的八阶Adams‑Bashforth格式；利用八阶Adams‑Bashforth格式计算航磁数据的向下延拓初步结果；利用积分迭代方式，对向下延拓初步结果进行迭代，计算向下延拓最终结果。本发明提供的航磁数据向下延拓方法在保证准确和稳定的前提下，提高了向下延拓的深度，对航空磁测数据具有增强的效果，进而实现航磁数据解释的精细化改善。

技术领域

本发明涉及勘查技术与工程技术领域，特别是涉及一种航磁数据向下延拓方法。

背景技术

传统地面磁法测量方式受到地表条件等制约因素影响，实施生产作业困难且效率不高。航空磁法测量是以有人驾驶飞机、无人驾驶飞机等为搭载平台，利用磁力仪等设备进行空中飞行采集地磁场数据。航磁测量可以解决地面磁测的沙漠、沼泽、密林、高山等复杂地表问题，但是地磁场随高度的增加呈指数衰减，导致航磁数据的信息被消减，将航磁数据归算到地面或者水准面可以增强信息，减少衰减效应。这种归算技术称为向下延拓。但是一般方法不能实现大深度(超过15倍点距)稳定的航磁数据向下延拓。徐世浙于2006年提出积分迭代向下延拓方法，对无噪声数据能够实现稳定向下延拓，且向下延拓点距大(甚至可达20倍点距)，但因其所用的起始迭代公式逼近效果差，迭代次数较多，导致计算效率降低和噪声累加(陈生昌和肖鹏飞，2007)；张冲等于2017年提出的三阶Adams-Bashforth向下延拓方法，利用实测的垂向导数，通过微分方程的数值解法Adams-Bashforth格式，建立了直接的向下延拓公式。但是三阶Adams-Bashforth向下延拓方法精度不高、且在实践过程中需要利用实测的航磁垂向导数。

发明内容

本发明要解决的技术问题是提供一种航磁数据向下延拓方法，在保证准确和稳定的前提下，提高了向下延拓的深度，对航空磁测数据具有增强的效果，进而实现航磁数据解释的精细化改善。

为解决上述技术问题，本发明提供了一种航磁数据向下延拓方法，所述方法包括：对观测面的航磁垂向导数进行向上延拓，得到分布在观测面以上，七个不同高度的航磁垂向导数；利用观测面的航磁数据、航磁垂向导数以及其向上延拓七个高度的航磁垂向导数，构成向下延拓的八阶Adams-Bashforth格式；利用八阶Adams-Bashforth格式计算航磁数据的向下延拓初步结果；利用积分迭代方式，对向下延拓初步结果进行迭代，计算向下延拓最终结果。

在一些实施方式中，八阶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轴坐标。