基于BTTB矩阵的快速高精度三维磁场正演

摘要：

本文改进了一种快速、高精度空间域三维正演算法,用来计算地下场源在水平观测面产生的磁异常ΔT场及其梯度场,以解决传统空间域正演计算效率低的问题.算法采用长方体对场源区域进行剖分,观测点与场源剖分单元体中心点在水平面上的投影重合.改进的算法具有以下三个特点:(1)采用无解析奇点的解析解公式计算磁异常,保证计算精度.(2)通过构造特殊的分块托普利兹(BTTB,Block-Toeplitz Toeplitz-Block)矩阵,利用其结构特性压缩核矩阵,并且用预先计算并存储中间变量,优化计算核矩阵的过程以提高计算效率.(3)基于BTTB矩阵的特殊性质,将核矩阵与磁化率向量的乘积转化为二维离散卷积的形式,因此能利用快速傅里叶变换进一步提高计算效率.模型实验显示,当剖分个数较多时,改进的快速正演算法比传统解析解方法快约5个数量级,比现有的8点高斯快速傅立叶变换(Gauss-FFT)正演算法快约两个数量级,而且绝对误差极小(最大约为10-6nT),同时将反演时核矩阵的内存占用降低约5个数量级,证明了该正演算法具有高精度、高效率、低存储量的优点.最后设计了一个合成模型实验,将改进后的快速正演算法运用到磁异常ΔT反演中,反演所得三维磁化率与真实模型特征一致,且大幅降低反演计算时间和内存占用,验证了快速正演算法的实用性.

关键词：磁异常、正演计算、BTTB矩阵、直立长方体、三维反演

所属期刊栏目：65

分类号：P631

资助基金：国家自然科学基金;国家自然科学基金;国家自然科学基金;中南大学研究生自主探索创新项目

在线出版日期：2022-03-31（万方平台首次上网日期，不代表论文的发表时间）

页数：共18页

页码：1107-1124