SEISMOLOGY AND GEOLOGY ›› 2010, Vol. 32 ›› Issue (3): 392-401.DOI: 10.3969/j.issn.0253-4967.2010.03.005

• Brief Report • Previous Articles     Next Articles

PARALLEL COMPUTATION OF 2-D MAGNETOTELLURIC FORWARD MODELING BASED ON MPI

LI Yan1, HU Xiang-yun1, WU Gui-ju1, YE Yi-xin2, LIAO Guo-zhong1   

  1. 1. Institute of Geophysics and Geomatics, China University of Geoscience, Wuhan 430074, China;
    2. School of Nuclear Engineering and Technology, East China Institute of Technology, Fuzhou 344000, China
  • Received:2010-02-23 Revised:2010-08-05 Online:2010-09-30 Published:2011-01-14

基于MPI的二维大地电磁正演的并行计算

李焱1, 胡祥云1, 吴桂桔1, 叶益信2, 廖国忠1   

  1. 1. 中国地质大学地球物理与空间信息学院, 武汉, 430074;
    2. 东华理工大学核工程技术学院地球物理系, 抚州, 344000
  • 作者简介:李焱,男,1984年生,中国地质大学(武汉)2008级地球探测与信息技术专业在读硕士研究生,主要研究方向为大地电磁正反演的并行计算,电话:13986283144,E-mail:liyan063042@gmail.com.
  • 基金资助:
    国家深部专项(Sinoprobe010302);国家自然科学基金(40974040)共同资助

Abstract: Nowadays,the 2-D inversion of magnetotelluric data has been put into practice,however,the 3-D inversion is still in its infancy,not because of the backward inversion methods or theories from the real yield but the lag of the hardware development.It is the huge amount of time to calculate in 2-D and 3-D inversion processes that traditional single-computer calculation speed is unacceptable.The most time-consuming part of the inversion is forward calculation,so it is of great significance to take parallel processing technology to improve the efficiency of parallel calculation.The program for 2-D magnetotelluric forward problems presented in this paper is based on the classic finite element method put forward by Wannamaker.It combines FORTRAN and MPICH 2 on Linux system in compilation.The forward computing for two-dimensional magnetotellurics is done separately in each frequency,and the value of electromagnetic field in the corresponding frequency is independent,therefore it is fit for frequency-decomposing parallel computing.This paper presents the two-dimensional magnetotelluric forward parallel algorithm,based on the finite element.Using the development tools of Fortran and MPICH 2,we write a parallel program that combines the advantages of MPI,and adopt the master-slave parallel mode and the parallel program of frequency-decomposing parallel computing.In order to verify the program,we designed two models,and compared with the serial program to validate the feasibility and accuracy of this algorithm,the result provides a research base for the parallel computation for two-dimensional inversion and three-dimensional forward and inversion.

Key words: magnetotellurics, MPI, two-dimensional, finite element, forward computation

摘要: 大地电磁二维正演对每一个频率分别进行计算,各频率对应的电磁场值间相互独立,根据这一特点可以将程序按频率划分粒度,将每个频点计算分配到各个进程同时进行计算,并行执行。文中给出了大地电磁二维有限元正演的并行算法,它结合了MPI的优点,采用主从并行模式、分频并行计算的并行方案来执行。为了检验编写的并行程序,对设计的2个模型进行试算,并与串行程序对比验证了该算法的可行性、正确性,为二维反演、三维正反演的并行计算提供了研究基础。

关键词: 大地电磁, MPI, 二维, 有限元, 正演

CLC Number: