A study on Parallel Computation Based on Finite Element Forward Modeling of 2D Magnetotelluric

Abstract

Magnetotelluric sounding method based on the difference of the rock’s resistivity is an exploration method about doing research in earth’s resistivity and phase using the native electromagnetic field. The paper adopts 2D finite element method as the magnetotelluric forward method and calculates the total field by primary field (also named background field) plus secondary field. We can get more accurate forward result through the finite element method and we can get the result effected by the dense degree of grid slightly by the total field. But the method is not effective enough when the model is divided into relative big grid. When the frequency changes, program solves relevant equation separately. According to the feature of the algorithm, we apply MPI parallel method in the algorithm. Every process solves relevant equation. The account of frequency that a process needs to solve in parallel computation is less than the account that the process needs to solve in serial algorithm. We can see that the forward result is the same with the serial algorithm and proves the correctness of algorithm. We do statistics about the efficiency of the parallel algorithm. When the account of processes is from 2 to 8, the speedup is from 1.63 to 2.64. It proves the effectiveness of the parallel algorithm.

Keywords

Share and Cite:

Wang, M. , Tan, H. , Lin, C. , Liu, X. and Zhang, Z. (2015) A study on Parallel Computation Based on Finite Element Forward Modeling of 2D Magnetotelluric. International Journal of Geosciences, 6, 863-868. doi: 10.4236/ijg.2015.68070.

Conflicts of Interest

The authors declare no conflicts of interest.

 [1] Zhdanov, M.S. and Tolstaya, E. (2004) Minimum Support Nonlinear Parametrization in the Solution of a 3D Magnetotelluric Inverse Problem. Inverse Problems, 20, 937-952. http://dx.doi.org/10.1088/0266-5611/20/3/017 [2] Tan, H.D., Yu, Q.F., Booker, J., et al. (2003) Three-Dimensional Rapid Relaxation Inversion for the Magnetotelluric Method. Chinese Journal of Geophysics, 46, 850-855. (In Chinese) http://dx.doi.org/10.1002/cjg2.442 [3] Zeng, G. (2008) 2D Magnetotelluric Forward Modeling. Central South University, Changsha, 5-6. [4] Xu, S.Z. (1994) The Finite Element Method in Geophysics. Science Press, 15-52. (In Chinese) [5] Rodi, W.L. (1976) A Technique for Improve the Accuracy of Finite Element Solution for Magnetotelluric Data. Geophysical Journal International, 44, 483-506. http://dx.doi.org/10.1111/j.1365-246X.1976.tb03669.x [6] Xie, F. (2006) The Study of 2D Magnetotelluric Forward Modeling with Terrain. China University of Geosciences, Wuhan, 9-21. [7] Liu, X.J. (2007) The Study of Magnetotelluric Forward Modeling Based on Second Field. Journal of Tongji University, 8, 1113-1116. [8] Nabighian (1998) Electromagnetic Methods. Theory Society of Exploration Geophysicists, 1, 123-143. [9] Liu, X.J. (2007) The Magnetotelluric Focusing Inversion. Tongji University, Shanghai, 1-42. (In Chinese) [10] Shi, Y.J. (1985) Magnetotelluric Sounding Method. School of Information Technology, Chengdu University of Technology, Chengdu, 50-53. [11] Zhou, X.X. (1986) The Numerical Simulation Technology of Electrical Prospecting. Sichuan Science and Technology Press, Chengdu, 57-85. [12] Du, Z.H. (2001) Parallel Computation Technology.