3D magnetotelluric inversions with unstructured finite-element and limited-memory quasi-Newton methods
Cao Xiao-Yue1, Yin Chang-Chun1, Zhang Bo1, Huang Xin1,2, Liu Yun-He1, and Cai Jing1
1. College of Geo-exploration Science and Technology, Jilin University, Changchun 130026, China.
2. Department of Earth Sciences, Memorial University of Newfoundland, St. John's, NL, A1B, Canada.
Abstract：
Traditional 3D Magnetotelluric (MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured-grid-based methods can model complex underground structures with high accuracy and overcome the defects of traditional methods, such as the high computational cost for improving model accuracy and the difficulty of inverting with topography. In this paper, we used the limited-memory quasi-Newton (L-BFGS) method with an unstructured finite-element grid to perform 3D MT inversions. This method avoids explicitly calculating Hessian matrices, which greatly reduces the memory requirements. After the first iteration, the approximate inverse Hessian matrix well approximates the true one, and the Newton step (set to 1) can meet the sufficient descent condition. Only one calculation of the objective function and its gradient are needed for each iteration, which greatly improves its computational efficiency. This approach is well-suited for large-scale 3D MT inversions. We have tested our algorithm on data with and without topography, and the results matched the real models well. We can recommend performing inversions based on an unstructured finite-element method and the L-BFGS method for situations with topography and complex underground structures.
. 3D magnetotelluric inversions with unstructured finite-element and limited-memory quasi-Newton methods[J]. APPLIED GEOPHYSICS, 2018, 15(3-4): 556-565.
[1]
Avdeev, D., and Avdeeva, A., 2009, 3D magnetotelluric inversion using a limited-memory quasi-Newton optimization: Geophysics, 74(3), 45−57.
[2]
Cao, X. Y., Yin, C., Zhang, B., Xin, H., Ren, X., and Qiu, C., et al., 2017, A goal-oriented adaptive finite-element algorithm for 3D anisotropic MT modelling: 79th EAGE Conference and Exhibition, H1−H4.
[3]
deGroot-Hedlin, C., and Constable, S., 1990, Occam’s inversion to generate smooth, two-dimensional models from magnetotelluric data: Geophysics,55(12), 1613−1624.
[4]
Dennis, J. E., and Schnabel, R. B., 1983, Numerical methods for unconstrained optimization and nonlinear equations, Englewood Cliffs, NJ: Prentice-Hall, New Jersey, 86−110.
[5]
Egbert, G. D., and Kelbert, A., 2012, Computational recipes for electromagnetic inverse problems:Geophysical Journal International,189(1), 251−267.
[6]
Grayver, A. V., 2015, Parallel three-dimensional magnetotelluric inversion using adaptive finite-element method, Part I: theory and synthetic study: Geophysical Journal International,202(1), 584−603.
[7]
Gürer, A., and Ìlk???k, O. M., 1997, The importance of topographic corrections on magnetotelluric response data from rugged regions of Anatolia: Geophysical prospecting, 45(1), 111−125.
[8]
Jahandari, H., and Farquharson, C. G., 2017, 3-D minimum-structure inversion of magnetotelluric data using the finite-element method and tetrahedral grids:Geophysical Journal International,211(2), 1189−1205.
[9]
Liu, Y. H., and Yin, C. C., 2013, 3D inversion for frequency-domain HEM data: Chinese J. Geophys. (in Chinese), 56(12), 4278-428.
[10]
Liu, D. C., and Nocedal, J., 1989, On the limited memory BFGS method for large scale optimization: Math. Program, 45, 503-528.
[11]
Mackie, R. L., and Madden, T. R., 1993, Three-dimensional magnetotelluric inversion using conjugate gradients: Geophysical Journal International,115(1), 215−229.
[12]
Nocedal, J., 1980, Updating quasi-Newton matrices with limited storage:Mathematics of computation,35(151), 773−782.
[13]
Newman, G. A., and Alumbaugh, D. L., 2000, Three-dimensional magnetotelluric inversion using non-linear conjugate gradients:Geophysical journal international,140(2), 410−424.
[14]
Newman, G. A., and Boggs, P. T., 2004, Solution accelerators for large-scale three-dimensional electromagnetic inverse problems:Inverse Problems,20(6), S151.
[15]
Rodi, W. L., and Mackie, R. L., 2000, Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion: Geophysics, 66(1), 174−187.
[16]
Ren, Z., Kalscheuer, T., Greenhalgh, S., and Maurer, H., 2013, A goal-oriented adaptive finite-element approach for plane wave 3-D electromagnetic modelling: Geophysical Journal International, 194(2), 700−718.
[17]
Smith, J. T., and Booker, J. R, 1991, Rapid inversion of two- and three- dimensional magnetotelluric data: Journal of Geophysical Research: Solid Earth,96(B3), 3905−3922.
[18]
Schwalenberg, K., and Edwards, R. N., 2004, The effect of seafloor topography on magnetotelluric fields: an analytical formulation confirmed with numerical results: Geophysical Journal International, 159(2), 607−621.
[19]
Siripunvaraporn, W., and Egbert, G., 2009, WSINV3DMT: vertical magnetic field transfer function inversion and parallel implementation: Physics of the Earth and Planetary Interiors, 173(3), 317−329.
[20]
Tikhonov, A. N., and Arsenin, V. Y., 1977, Solutions to Ill-posed Problems: John Wiley and Sons, New York, NY, 45−87.
[21]
Usui, Y., 2015, 3-D inversion of magnetotelluric data using unstructured tetrahedral elements: applicability to data affected by topography: Geophysical Journal International, 202(2), 828−849.
[22]
Yin, C. C., Zhang, B., and Liu, Y. H., 2017, A goal-oriented adaptive algorithm for 3D magnetotelluric forward modeling: Chinese J. Geophys (in Chinese), 60(1), 327−336.
2018, Vol. 15 
