论文信息 - A NLCG 3‐D Inversion Method of Magnetotellurics with Parallel Structure

A NLCG 3‐D Inversion Method of Magnetotellurics with Parallel Structure

In this paper, we improve and verify the inversion algorithm and parallel computing program of three-dimensional magnetotelluric sounding, and the program does not require much physical memory and CPU quantity. It makes possible to do three-dimensional inversion on ordinary PC computers. Throughout nonlinear conjugate gradient algorithm proposed by Newman and Alumbaugh (2000) and preconditioning method proposed by Rodi and Mackie (2001), we improve the preconditioning method add parallel structure, and verify the program with theoretical models. According to the inversion results of measured data in Japanese KAYABE area, we demonstrate the practicality of the algorithm.

[1] Gerald W. Hohmann,et al. Magnetotelluric responses of three-dimensional bodies in layered earths , 1982 .

[2] C. M. Reeves,et al. Function minimization by conjugate gradients , 1964, Comput. J..

[3] J. T. Smith,et al. Rapid inversion of two‐ and three‐dimensional magnetotelluric data , 1991 .

[4] R. Mackie,et al. Three-dimensional magnetotelluric inversion using conjugate gradients , 1993 .

[5] G. Newman,et al. Three-dimensional magnetotelluric inversion using non-linear conjugate gradients , 2000 .

[6] S. Constable,et al. Occam's inversion to generate smooth, two-dimensional models from magnetotelluric data , 1990 .

[7] J. T. Smith. Conservative modeling of 3-D electromagnetic fields, Part II: Biconjugate gradient solution and an accelerator , 1996 .

[8] J. T. Smith. Conservative modeling of 3-D electromagnetic fields, Part I: Properties and error analysis , 1996 .

[9] Yongwimon Lenbury,et al. Three-dimensional magnetotelluric inversion : data-space method , 2005 .

[10] William Rodi,et al. Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion , 2001 .