首页> 中文期刊> 《地球物理学报》 >基于BTTB矩阵的快速高精度三维磁场正演

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

         

摘要

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

著录项

相似文献

  • 中文文献
  • 外文文献
  • 专利
获取原文

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号