The 2‐D magnetotelluric inverse problem solved with optimization

作者: Ashley E. Van Beusekom , Robert L. Parker , Randolph E. Bank , Philip E. Gill , Steven Constable

DOI: 10.1111/J.1365-246X.2010.04895.X

关键词: Optimization problemMathematical analysisNonlinear systemPartial differential equationMathematicsBoundary value problemRegularization (mathematics)Multigrid methodUpper and lower boundsInverse problem

摘要: SUMMARY The practical 2-D magnetotelluric inverse problem seeks to determine the shallow-Earth conductivity structure using finite and uncertain data collected on ground surface. We present an approach based PLTMG (Piecewise Linear Triangular MultiGrid), a special-purpose code for optimization with second-order partial differential equation (PDE) constraints. At each frequency, electromagnetic field are treated as unknowns in which misfit is minimized subject constraints that include Maxwell's equations boundary conditions. Within this framework it straightforward accommodate upper lower bounds or other conditions conductivity. In addition, underlying ill-posed, may be used apply various kinds of regularization. discuss some advantages difficulties associated PDE-constrained basis solving large-scale nonlinear geophysical problems. Combined transverse electric magnetic complex admittances from COPROD2 inverted. First, we invert penalizing size roughness giving solutions similar those found previously. second example, conventional regularization replaced by technique imposes model. both examples better than obtained previously, without any increase model complexity.

参考文章(27)
Randolph E. Bank, Philip E. Gill, Roummel F. Marcia, Interior Methods For a Class of Elliptic Variational Inequalities Springer, Berlin, Heidelberg. pp. 218- 235 ,(2003) , 10.1007/978-3-642-55508-4_13
Philip E. Gill, Walter Murray, Michael A. Saunders, SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization SIAM Review. ,vol. 47, pp. 99- 131 ,(2005) , 10.1137/S0036144504446096
P. Weidelt, P. Kaikkonen, Local 1-D interpretation of magnetotelluric B-polarization impedances Geophysical Journal International. ,vol. 117, pp. 733- 748 ,(1994) , 10.1111/J.1365-246X.1994.TB02466.X
Randolph E. Bank, Michael Holst, A New Paradigm for Parallel Adaptive Meshing Algorithms SIAM Review. ,vol. 45, pp. 291- 323 ,(2003) , 10.1137/S003614450342061
M.N. Berdichevsky, V.I. Dmitriev, E.E. Pozdnjakova, On two-dimensional interpretation of magnetotelluric soundings Geophysical Journal International. ,vol. 133, pp. 585- 606 ,(1998) , 10.1046/J.1365-246X.1998.01333.X
P. Weidelt, Inversion of two-dimensional conductivity structures Physics of the Earth and Planetary Interiors. ,vol. 10, pp. 282- 291 ,(1975) , 10.1016/0031-9201(75)90054-0
M.E. Everett, Homotopy, Polynomial Equations, Gross Boundary Data, and Small Helmholtz Systems Journal of Computational Physics. ,vol. 124, pp. 431- 441 ,(1996) , 10.1006/JCPH.1996.0070
José M. Romo, Enrique Gómez-Treviño, Francisco J. Esparza, Series and parallel transformations of the magnetotelluric impedance tensor: theory and applications Physics of the Earth and Planetary Interiors. ,vol. 150, pp. 63- 83 ,(2005) , 10.1016/J.PEPI.2004.08.021
J. Torquil Smith, John R. Booker, Rapid inversion of two- and three-dimensional magnetotelluric data Journal of Geophysical Research: Solid Earth. ,vol. 96, pp. 3905- 3922 ,(1991) , 10.1029/90JB02416