Smoothest Model and Sharp Boundary Based Two‐Dimensional Magnetotelluric Inversion

How to get rapid and stable inversion results and clear geological interfaces is a focus problem in current Magnetotelluric (MT) inversion. With the purpose to represent clear boundary with resistivity discontinuity, the paper integrates the ideas of OCCAM inversion and sharp boundary inversion (SBI) to construct a new objective function for inversion by using the minimum gradient support functional and conjugate gradient inversion method to solve the optimal minimum. So the 2-D MT inversion method based on OCCAM and SBI is fulfilled. Through the model tests and comparisons with the related published inversion results, this method is proved accurate in determining sharp electrical interfaces based on the smoothest and stable models. It is also applied to process the field observed MT data in the regions of Guangdong Xuwen, the results indicate this method's effect and feasibility.

[1]  G. Newman,et al.  Three-dimensional magnetotelluric inversion using non-linear conjugate gradients , 2000 .

[2]  Steven Constable,et al.  Inversion of magnetotelluric data for 2D structure with sharp resistivity contrasts , 2004 .

[3]  M. Zhdanov,et al.  3‐D magnetic inversion with data compression and image focusing , 2002 .

[4]  S. Constable,et al.  Occam's inversion to generate smooth, two-dimensional models from magnetotelluric data , 1990 .

[5]  Michael S. Zhdanov,et al.  Focusing geophysical inversion images , 1999 .

[6]  Xiao‐Ping Wu,et al.  Study on 3‐D Resistivity Inversion Using Conjugate Gradient Method , 2000 .

[7]  Michael S. Zhdanov,et al.  Three-dimensional regularized focusing inversion of gravity gradient tensor component data , 2004 .

[8]  Liang‐Jun Yan,et al.  Nonlinear Inversion with Quadratic Function Approaching Method for Magnetotelluric Data , 2004 .

[9]  Philip E. Wannamaker,et al.  A stable finite element solution for two-dimensional magnetotelluric modelling , 1987 .

[10]  Gary D. Egbert,et al.  An efficient data-subspace inversion method for 2-D magnetotelluric data , 2000 .

[11]  R. Parker,et al.  Occam's inversion; a practical algorithm for generating smooth models from electromagnetic sounding data , 1987 .

[12]  Douglas W. Oldenburg,et al.  Calculation of sensitivities for the frequency-domain electromagnetic problem , 1994 .

[13]  Francis T. Wu,et al.  THE INVERSE PROBLEM OF MAGNETOTELLURIC SOUNDING , 1968 .

[14]  Hui Yang,et al.  Constrained Joint Inversion of Magneto-telluric and Seismic Data Using Simulated Annealing Algorithm , 2002 .

[15]  J. Wang.  NEW DEVELOPMENT OF MAGNETOTELLURIC SOUNDING IN CHINA , 1997 .

[16]  Philip E. Wannamaker,et al.  Calculating the two‐dimensional magnetotelluric Jacobian in finite elements using reciprocity , 1996 .

[17]  P. Xu.  IMPROVEMENT OF OCCAM'S INVERSION FOR MT DATA , 1998 .

[18]  Gerald W. Hohmann,et al.  A rapid inversion technique for transient electromagnetic soundings , 1989 .

[19]  J. T. Smith,et al.  Rapid inversion of two‐ and three‐dimensional magnetotelluric data , 1991 .

[20]  R. Wu,et al.  Wave Equation Based Pre‐Stack Depth Migration Using Generalized Screen Propagator , 2002 .

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

[22]  D. A. H. Jacobs,et al.  A Generalization of the Conjugate-Gradient Method to Solve Complex Systems , 1986 .

[23]  J. Coggon.  Electromagnetic and electrical modeling by the finite element method , 1971 .

[24]  D. Oldenburg,et al.  Approximate sensitivities for the electromagnetic inverse problem , 1996 .