Simulated annealing for controlled-source audio-frequency magnetotelluric data inversion

Simulated annealing (SA) is used to invert 1D controlled source audio-frequency magnetotelluric (CSAMT) data. In the annealing process, the system energy is taken as the root-mean-square fitting error between model responses and real data. The model parameters are the natural logarithms of the resistivity and the thickness in each layer of the earth. The annealing temperature decreases exponentially, while the model is refreshed randomly according to the temperature and is accepted according to a Boltzmann probability. We first tested the SA on synthetic data and developed a cooling schedule of model updates specifically for CSAMT data inversion. The redesigned cooling schedule reduces the magnitude of the model updating, and makes the solution converge rapidly and stably. For a three-layer model whose resistivity increases with depth, SA has difficulty in obtaining the global solution for the middle layer. However, the solution for such a layer can be significantly improved by using the mean value of the...

[1]  Mrinal K. Sen,et al.  2-D resistivity inversion using spline parameterization and simulated annealing , 1996 .

[2]  Rapid relaxation inversion of CSAMT data , 1999 .

[3]  Jiansheng Wu,et al.  Constrained Joint Inversion of Gravity and Seismic Data Using the Simulated Annealing Algorithm , 2007 .

[4]  Diquan Li,et al.  The Application of Genetic Algorithm to CSAMT Inversion for Minimum Structure , 2008 .

[5]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[6]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[7]  Walter L. Anderson,et al.  A hybrid fast Hankel transform algorithm for electromagnetic modeling , 1989 .

[8]  Greg Hodges,et al.  Simulated annealing for airborne EM inversion , 2007 .

[9]  Kris Vasudevan,et al.  A comparison of several cooling schedules for simulated annealing implemented on a residual statics problem , 1993 .

[10]  Nils Ryden,et al.  Fast simulated annealing inversion of surface waves on pavement using phase-velocity spectra , 2006 .

[11]  S. Sharma,et al.  Global nonlinear inversion of transient EM data from conducting surroundings using a free‐space plate model , 2000 .

[12]  Comparison of induced polarization and controlled-source audio-magnetotellurics methods for massive chalcopyrite exploration in a volcanic area , 1997 .

[13]  Mrinal K. Sen,et al.  Global Optimization Methods in Geophysical Inversion , 1995 .

[14]  P. Wannamaker.  Tensor CSAMT survey over the Sulphur Springs thermal area, Valles Caldera, New Mexico, U.S.A., Part I: Implications for structure of the western caldera , 1997 .

[15]  J. J. Zach,et al.  Inversion of marine CSEM data using up‐down wavefield separation and simulated annealing , 2007 .

[16]  Douglas W. Oldenburg,et al.  Inversion of controlled source audio-frequency magnetotellurics data for a horizontally layered earth , 1999 .

[17]  L. Cagniard.  Basic theory of the magneto-telluric method of geophysical prospecting , 1953 .

[18]  J. Wright,et al.  Tensor CSAMT studies at the Buchans Mine in central Newfoundland , 1993 .

[19]  L. Ingber.  Very fast simulated re-annealing , 1989 .