2. 5 dimension electromagnetic forward parallel computing based on MPI + OpenMP
Bai Hongtao1,3,4, Fan Cuisong2, Li Ang4,5, He Lili4,5, Li Tonglin1
1. College of Earth Survey Science and Technology, Jilin University, Changchun, Jilin 130026, China;
2. Tianjin Center, China Geological Survey, Tianjin 300170, China;
3. Center for Computer Fundamental Education, Jilin University, Changchun, Jilin 130012, China;
4. Key Laboratory of Symbolic Computation and Knowledge Engineering of Ministry of Education, Jilin University, Changchun, Jilin 130012, China;
5. College of Computer Science and Technology, Jilin University, Changchun, Jilin 130012, China
Abstract:We propose in this paper a complex resistivity 2. 5-Dimension electromagnetic forward parallel computing algorithm based on MPI + OpenMP to improve low efficiency of large-scale numerical computation in geophysical exploration. This hybrid algorithm puts different wave domains to nodes of computer clusters using MPI. In the algorithm interior, the computation of sparse linear equations with the incomplete Cholesky biconjugate gradient (ICBCG) method is parallized on OpenMP. All waves are assigned to the nodes with an equal interval to have load balance. And to reduce communication overhead, each node reads its own model parameters. Different experiments demonstrate that our parallel algorithm can not only acquire the same data accuracy as that of serial computation methods, but also gain a near-linear speed up.
Wang Dayong, Li Tonglin, Gao Yuan et al. The application of CASMT and TEM to exploration buried deposits in Longhu mountain area at Tongling, Anhui Province. Journal of Jilin University (Earth Science Edition), 2009, 39(6):1134-1140.
[2]
An Zhiguo, Di Qingyun. Application of the CSAMT method for exploring deep coal mines in Fujian province, Southeastern China. Journal of Environmental and Engineering Geophysics, 2010, 15(4):243-249.
Huang Yi, Shi Xueming, Fan Jianke et al. Review on parallel computing and its application in exploration geophysics. Progress in Geophysics, 2010, 25(2):642-649.
Gu Guanwen, Liang Meng, Wu Wenli. Pseudo-2D inversion interpretation of CSAMT data based on parallel computation. Geophysical & Geochemical Exploration, 2010, 34(3):399-402.
[6]
Coggon J H. Electromagnetic and electrical modeling by the finite element method. Geophysics, 1971,36(1): 132-155.
[7]
Unsworth J M, Bryan J T, Alan D C. Electromagnetic induction by a finite electric dipole source over a 2-D earth. Geophysics, 1993, 58(2):198-214.
[8]
Mitsuhata Y. 2-D electromagnetic modeling by finite-element method with a dipole source and topography. Geophysics, 2000, 65(2):465-475.
[9]
Yuguo L, Kerry K. 2D marine controlled-source electromagnetic modeling: Part 1-An adaptive finite-element algorithm. Geophysics, 2007,72(2): 51-62.
[10]
底青云, Martyn U, 王妙月. 复杂介质有限元法2.5维可控源音频大地电磁法数值模拟.地球物理学报, 2004,47(4):723-730.
Di Qingyun, Martyn U, Wang Miaoyue. 2.5-D CSAMT modeling with the finite element method over 2-D complex earth media. Chinese Journal of Geophysics, 2004,47(4):723-730.
Zhang Jifeng, Tang Jingtian Yu Yan et al. Finite element numerical simulation on line controlled source based on quadratic interpolation. Journal of Jilin University (Earth Science Edition), 2009,39(5):929-935.
Li Yan, Hu Xiangyun, Yang Wencai et al. A study on parallel computation for 3D magnetotelluric modeling using the staggered-grid finite difference method. Chinese Journal of Geophysics, 2012,55(12):4036-4043.
[16]
Tan Handong, Tong Tuo, Lin Changhong. The parallel 3D magnetotelluric forward modeling algorithm. Applied Geophysics, 2006,3(4):197-202.