肖晓^1,2,3,,
黄保尚¹,
任政勇^1,2,3,,,
汤井田^1,2,3
1. 中南大学地球科学与信息物理学院, 长沙 410083
2. 中南大学有色金属成矿预测与地质环境监测教育部重点实验室, 长沙 410083
3. 有色资源与地质灾害探查湖南省重点实验室, 长沙 410083

基金项目: 国家自然科学基金（41574120，4171101400，41830107），国家高技术研究发展计划（2014AA06A602），湖南省自然科学基金（2016JJ2139），中南大学创新驱动计划（2016CX005），青年973（2015CB060200）资助


详细信息

作者简介: 肖晓, 男, 1981年生, 副教授, 主要从事电磁法正反演及数据处理研究.E-mail:csuxiaox@csu.edu.cn

通讯作者: 任政勇, 男, 1983年生, 教授, 主要从事电磁场数值模拟及位场正反演理论的研究.E-mail:renzhengyong@csu.edu.cn

中图分类号: P631

收稿日期:2017-07-10
修回日期:2018-06-20
上线日期:2019-03-05

Large-scale forward modeling of magnetic data using an adaptive multi-level fast multipole method

XIAO Xiao^1,2,3,,
HUANG BaoShang¹,
REN ZhengYong^1,2,3,,,
TANG JingTian^1,2,3
1. School of Geosciences and Info-physics of Central South University, Changsha 410083, China
2. The Key Laboratory of Metallogenic Prediction of Nonferrous Metals of Ministry of Education, Central South University, Changsha 410083, China
3. Key Laboratory of Non-Ferrous Resources and Geological Hazard Detection, Changsha 410083, China


More Information

Corresponding author: REN ZhengYong,E-mail:renzhengyong@csu.edu.cn

MSC: P631

--> Received Date: 10 July 2017
Revised Date: 20 June 2018
Available Online: 05 March 2019


摘要
摘要:提出了一种基于非结构化四面体以及带地形模型的自适应多层快速多极大规模磁法快速正演算法.该算法弥补了传统积分方法采用FFT加速计算时不能采用非结构化网格的缺陷；同时采用自适应快速多极算法突破积分求和方法求解大规模磁法问题耗时长的突出问题.首先，采用非结构化的四面体网格剖分技术能够更好的模拟复杂模型以及带地形模型，实现磁法模型的高精度模拟；其次，采用一种自适应多层快速多极（AMFM）算法实现大规模磁法正演求解.通过将计算区域划分为近区和远区，对近区采用解析计算高精度求解，对远区采用自适应多层快速多极算法进行加速计算，假设有M个观测点，N个四面体源单元，可将计算复杂度由传统积分求和法的O（MN）减少到O（MlogN）.本文设计了组合体模型以及安徽怀宁地区的实际地形模型，模型计算结果体现了采用该方法进行大规模复杂模型三维磁法正演模拟的高效性和准确性.
关键词: 磁法正演/
自适应多层快速多极算法/
四面体网格/
地形

Abstract:In this paper, we propose a fast and high-precision approach for large-scale magnetic forward modeling using the adaptive multi-level fast multipole (AMFM) method and an unstructured grid. The algorithm can overcome the shortcomings of the Fast fourier transfom (FFT) method which cannot adopt unstructured grids. In addition, the adaptive multi-level fast multipole (AMFM) method is chosen to solve the problem of high computational cost in traditional integral summation methods. Firstly, the unstructured tetrahedral meshing technique is used to properly approximate the complex model, for instance the terrain model, and realize the high-precision simulation of the magnetic model. Secondly, an adaptive multi-level fast multipole (AMFM) method is employed to accelerate large-scale magnetic forward modeling. Then, we divide the integral region into two parts. One is the near part in which we use analytical solution to get high-precision results, and the other is far part in which we choose the the adaptive multi-level fast multipole method to accelerate the calculation. The computational complexity can be reduced from in the traditional integral summation method to, where M and N are the numbers of observation sites and source elements, respectively. A combinational model and a complex model based on the digital elevation model (DEM) of Huaining County in Anhui province are used to test the proposed method. The results indicate that the proposed method is effective and accurate in three-dimensional magnetic forward modeling.
Key words:Magnetic forward modeling/
Adaptive multi-level fast multipole (AMFM) method/
Tetrahedral mesh/
Terrain



PDF全文下载地址:

http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_14909