基于大地电磁二维反演的MPI并行算法研究  

Study on parallel algorithm based on inversion of 2D magnetotelluric

作　　者：汪茂[1,2] 陈霜 谭捍东 刘晓[4] WANG Mao;CHEN Shuang;TAN Han-dong;LIU Xiao(China University of Geosciences, Beijing 100083, China;Key Laboratory of Geo-detection （ China University of Geosciences, Beijing ） , Beijing 100083, China;Geological Exploration Party No. 208, CNNC, Baotou 014010, China;Nanchang Institute of Technology, Nanchang 330099, China)

机构地区：[1]中国地质大学(北京),北京100083 [2]地下信息探测技术与仪器教育部重点实验室,北京100083 [3]核工业二〇八大队,包头014010 [4]南昌工程学院,南昌330099

出　　处：《地球物理学进展》2017年第5期2085-2090,共6页Progress in Geophysics

基　　金：中国地质大学(北京)基本科研启动基金项目(2-9-2017-096)资助

摘　　要：大地电磁测深法是以岩石的电性差异为基础和前提的勘探方法.本文所采用的大地电磁二维反演方法为共轭梯度法,该方法避免了求解雅可比矩阵,效率较高,但是在将模型剖分为较细网格,多频率进行计算时效率有待提高.基于大地电磁根据各频率依次独立处理数据的特点,在本文中采用了MPI并行运算方法,用多个进程同时来计算各频率数据,最后再将数据进行收集,得到最后的计算结果.通过对正演和反演结果的比较,验证了程序的正确性.对并行算法的效率进行了统计,进程数为2~8时,加速比能达到1.63~2.64,验证了并行算法的有效性.Magnetotelluric sounding method based on the difference of the rock＇s resistivity is an exploration method which do research in earth ＇s resistivity and phase using the native electromagnetic field. The paper introduces the inversion algorithm of 2 D magnetotelluric named nonlinear conjugate gradient method,the method which avoid resolving the jacobi matrix is very effective,but the method is not effective enough when the model is divided into dense grid. This study can do more computation during the time using the parallel computation,and then gather and sort the data processed. We compare the results of the serial algorithm with the result of the parallel algorithm,it proves that the parallel algorithm is correct. When the number of the processes is 2 to 8,the speedup ratio is 1. 63 to 2. 64. It improves the effectivity of the parallel algorithm of 2 D magnetotelluric.

关 键 词：大地电磁 二维正反演 共轭梯度 并行运算 

分 类 号：P631[天文地球—地质矿产勘探]