Abstract:Up until now, the key issue to the practical applications of three-dimensional magnetotelluric (MT) inversion is the insufficiency in computing resources. By further analysis and understanding on data-space inversion approach of 3D MT, we develop a massively parallel inversion scheme on the basis of frequency division and matrix decomposition, and implement its procedure by using MPI on TC5000A high-performance computing platform. The algorithm we develop includes the parallel calculation of three dimensional forward modeling, and sensitivity matrix and cross-product matrix, as well as the update of model parameters. The algorithm has the advantages of higher efficiency in computation and lower memory storage in which the storage amount of sensitivity matrix in every single computing node is 2/N times that on a PC (N is the number of nodes included in parallel computation). Furthermore, we test the implemented scheme with synthetic data from two 3D theoretical models and analyze the computational efficiency under multiple-nodes computing. The numerical experiment results show that the 3D data-space parallel inversion algorithm is feasible and efficient. Compared with the implementation on single PC, the parallel scheme is not only able to improve the computing speed and shorten the computation time, but also enlarge the calculational scale, which would advance the practicality of three dimensional magnetotelluic inversion.

Three-Dimensional Magnetotelluic Parallel Inversion Algorithm Using Data Space Method