Please use this identifier to cite or link to this item: http://hdl.handle.net/2122/3509
DC FieldValueLanguage
dc.contributor.authorallCerv, V.; Geophysical Institute, Academy of Sciences of the Czech Republic, Prague 4, Czech Republicen
dc.contributor.authorallMenvielle, M.; Centre d’Études des Environnements Terrestre et Planétaire, Saint Maur des Fosses Cedex, Franceen
dc.contributor.authorallPek, J.; Geophysical Institute, Academy of Sciences of the Czech Republic, Prague 4, Czech Republicen
dc.date.accessioned2007-12-20T13:45:31Zen
dc.date.available2007-12-20T13:45:31Zen
dc.date.issued2007-02en
dc.identifier.urihttp://hdl.handle.net/2122/3509en
dc.description.abstractGlobal optimization and stochastic approaches to the interpretation of measured data have recently gained particular attraction as tools for directed search for and/or verification of characteristic structural details and quantitative parameters of the deep structure, which is a task often arising when interpreting geoelectrical induction data in seismoactive and volcanic areas. We present a comparison of three common global optimization and stochastic approaches to the solution of a magnetotelluric inverse problem for thick layer structures, specifically the controlled random search algorithm, the stochastic sampling by the Monte Carlo method with Markov chains and its newly suggested approximate, but largely accelerated, version, the neighbourhood algorithm. We test the algorithms on a notoriously difficult synthetic 5-layer structure with two conductors situated at different depths, as well as on the experimental COPROD1 data set standardly used to benchmark 1D magnetotelluric inversion codes. The controlled random search algorithm is a fast and reliable global minimization procedure if a relatively small number of parameters is involved and a search for a single target minimum is the main objective of the inversion. By repeated runs with different starting test model pools, a sufficiently exhaustive mapping of the parameter space can be accomplished. The Markov chain Monte Carlo gives the most complete information for the parameter estimation and their uncertainty assessment by providing samples from the posterior probability distribution of the model parameters conditioned on the experimental data. Though computationally intensive, this method shows good performance provided the model parameters are sufficiently decorrelated. For layered models with mixed resistivities and layer thicknesses, where strong correlations occur and even different model classes may conform to the target function, the method often converges poorly and even very long chains do not guarantee fair distributions of the model parameters according to their probability densities. The neighbourhood resampling procedure attempts to accelerate the Monte Carlo simulation by approximating the computationally expensive true target function by a simpler, piecewise constant interpolant on a Voronoi mesh constructed over a set of pre-generated models. The method performs relatively fast but seems to suggest systematically larger uncertainties for the model parameters. The results of the stochastic simulations are compared with the standard linearized solutions both for thick layer models and for smooth Occam solutions.en
dc.language.isoEnglishen
dc.relation.ispartofseries1/50 (2007)en
dc.subjectmagnetotelluric methoden
dc.subjectinverse problemen
dc.subjectcontrolled random searchen
dc.subjectMarkov chain Monte Carloen
dc.subjectneighbourhood algorithmen
dc.titleStochastic interpretation of magnetotelluric data, comparison of methodsen
dc.typearticleen
dc.type.QualityControlPeer-revieweden
dc.subject.INGV01. Atmosphere::01.03. Magnetosphere::01.03.06. Instruments and techniquesen
dc.relation.referencesCONSTABLE, S.C., R.L. PARKER and C.G. CONSTABLE (1987): Occam’s inversion: a practical algorithm for generating smooth models from EM sounding data, Geophysics, 52, 289-300. GELMAN, A., J.B. CARLIN, H.S. STERN and D.B. RUBIN (1995): Bayesian Data Analysis (Chapman & Hall, New York), pp. 552. GRANDIS, H., M. MENVIELLE and M. ROUSSIGNOL (1999): Bayesian inversion with Markov chains, I. The magnetotelluric one-dimensional case, Geophys. J. Int., 138, 757-768. JONES, A.G. and R. HUTTON (1979a): A multi-station magnetotelluric study in Southern Scotland, I. Fieldwork, data analysis and results, Geophys. J. R. Astron. Soc., 56, 329-349. JONES, A.G. and R. HUTTON (1979b): A multi-station magnetotelluric study in Southern Scotland, II. Monte-Carlo inversion of the data and its geophysical and tectonic implication, Geophys. J. R. Astron. Soc., 56, 351-358. MALINVERNO, A. (2002): Parsimonious Bayesian Markov chain Monte Carlo inversion in a nonlinear geophysical problem, Geophys. J. Int., 151, 675-688. MEJU, M.A. and V.R.S. HUTTON (1992): Iterative mostsquares inversion: application to magnetotelluric data, Geophys. J. Int., 108, 758-766. MENKE, W. (1989): Geophysical Data Analysis: Discrete Inverse Theory (Academic Press, London), 2nd edition, pp. 289. MONTAZ, A., A. TORN and S. VIITANEN (1997): A numerical comparison of some modified controlled random search algorithms, TUCS Technical Report No. 98. MOSEGAARD, K. and A. TARANTOLA (1995): Monte Carlo sampling of solutions to inverse problems, J. Geophys. Res., 100, 12431-12447. PORTNIAGUINE, O. and M.S. ZHDANOV (1999): Focusing geophysical inversion images, Geophysics, 64, 874-887. POUS, J., X. LANA and A.M. CORREIG (1985): Generation of Earth stratified models compatible with both ellipticity and phase velocity observations of Rayleigh waves, Pure Appl. Geophys., 123, 870-881. PRICE, W.L. (1977): A controlled random search procedure for global optimization, Computer J., 20, 367-370. SAMBRIDGE, M. (1999a): Geophysical inversion with a neighbourhood algorithm, I. Searching a parameter space, Geophys. J. Int., 138, 479-494. SAMBRIDGE, M. (1999b): Geophysical inversion with neighbourhood algorithm, II. Appraising the ensemble, Geophys. J. Int., 138, 727-746. SENN, M.K. and P.L. STOFFA (1995): Global optimization methods in geophysical inversion, Advances in Exploration Geophysics (Elsevier, Amsterdam), 4, pp. 294. TVRDÍK, J., L. MISˇÍK and I. KRˇIVY´ (2002): Competing heuristics in evolutionary algorithms, in Intelligent Technologies: from Theory to Applications, edited by P. SINCÁK, V. KVASNICKA, J. VASCÁK and J. POSPÍCHAL (IOS Press, Amsterdam), 159-165. WEAVER, J.T. and A.K. AGARWAL (1993): Automatic 1D inversion of magnetotelluric data by the method of modelling, Geophys. J. Int., 112, 115-123.en
dc.description.journalTypeJCR Journalen
dc.description.fulltextopenen
dc.contributor.authorCerv, V.en
dc.contributor.authorMenvielle, M.en
dc.contributor.authorPek, J.en
dc.contributor.departmentGeophysical Institute, Academy of Sciences of the Czech Republic, Prague 4, Czech Republicen
dc.contributor.departmentCentre d’Études des Environnements Terrestre et Planétaire, Saint Maur des Fosses Cedex, Franceen
dc.contributor.departmentGeophysical Institute, Academy of Sciences of the Czech Republic, Prague 4, Czech Republicen
item.languageiso639-1en-
item.openairecristypehttp://purl.org/coar/resource_type/c_18cf-
item.cerifentitytypePublications-
item.fulltextWith Fulltext-
item.openairetypearticle-
item.grantfulltextopen-
crisitem.author.deptGeophysical Institute, Academy of Sciences of the Czech Republic, Prague, Czech Republic-
crisitem.author.deptCentre d’Études des Environnements Terrestre et Planétaire, Saint Maur des Fosses Cedex, France-
crisitem.author.deptGeophysical Institute, Academy of Sciences of the Czech Republic, Prague, Czech Republic-
crisitem.classification.parent01. Atmosphere-
Appears in Collections:Annals of Geophysics
Files in This Item:
File Description SizeFormat
2 Cerv.pdf1.37 MBAdobe PDFView/Open
Show simple item record

Page view(s) 50

166
checked on Jul 4, 2022

Download(s) 20

320
checked on Jul 4, 2022

Google ScholarTM

Check