Inversion of Magnetotelluric Measurements Using Multigoal Oriented hp-adaptivity
Autor: | Hélène Barucq, Julen Álvarez-Aramberri, David Pardo |
---|---|
Přispěvatelé: | University of the Basque Country/Euskal Herriko Unibertsitatea (UPV/EHU), Laboratoire de Mathématiques et de leurs Applications [Pau] (LMAP), Université de Pau et des Pays de l'Adour (UPPA)-Centre National de la Recherche Scientifique (CNRS), Advanced 3D Numerical Modeling in Geophysics (Magique 3D), Université de Pau et des Pays de l'Adour (UPPA)-Centre National de la Recherche Scientifique (CNRS)-Université de Pau et des Pays de l'Adour (UPPA)-Centre National de la Recherche Scientifique (CNRS)-Inria Bordeaux - Sud-Ouest, Institut National de Recherche en Informatique et en Automatique (Inria)-Institut National de Recherche en Informatique et en Automatique (Inria), Vassil Alexandrov, Michael Lees, Valeria Krzhizhanovskaya and Jack Dongarra and Peter M.A. Sloot, University of the Basque Country [Bizkaia] (UPV/EHU) |
Jazyk: | angličtina |
Rok vydání: | 2013 |
Předmět: |
Computer science
Computation Inverse Problems Constrained optimization Inversion (meteorology) 010103 numerical & computational mathematics hp-adaptivity Inverse problem Conductivity 01 natural sciences 010101 applied mathematics Finite Element Method (FEM) Magnetotellurics Electric field Goal-Oriented adaptivity General Earth and Planetary Sciences 0101 mathematics Electrical impedance Algorithm [MATH.MATH-NA]Mathematics [math]/Numerical Analysis [math.NA] General Environmental Science |
Zdroj: | Procedia Computer Science ICCS 2013-International Conference on Computational Science ICCS 2013-International Conference on Computational Science, Jun 2013, Barcelona, Spain. pp.1564-1573, ⟨10.1016/j.procs.2013.05.324⟩ ICCS |
DOI: | 10.1016/j.procs.2013.05.324⟩ |
Popis: | International audience; The inversion of magnetotelluric measurements require the computation of the quantity of interest (solution at the receivers) at different positions. Using a multigoal-oriented adaptive strategy for an hp Finite Element algorithm, we ensure accurate solutions in all receivers. With them, and using a limited memory algorithm for bound constraint optimization (L-BFGS-U method) to solve the Inverse Problem, we show that we are able to invert properly the simulated measurements in order to recover the subsurface conductivity distribution. We emphasize that the choice of the unknown of the Inverse Problem affects significantly the convergence of the inversion. We also exhibit the limitations of the electric field as the quantity of interest. When the distribution of the subsurface contains big jumps, the solution of the problem is not unique. The fact that the quantity of interest is not sensible enough to the variation of the conductivity suggests the use of the impedance as quantity of interest as a candidate to bring better results in a future research. |
Databáze: | OpenAIRE |
Externí odkaz: |