Document Type : Original Research Paper

Authors

1 Petroleum and Geophysics Engineering, Shahrood University of Technology, Shahrood, Iran

2 Science Faculty, Arak University, Arak, Iran

3 Geological Survey of Iran, Tehran, Iran

Abstract

In this paper we used orthogonal basis functions and expansion coefficients for inverse modeling of magnetic data. The basis functions chosen are normalized eigenvectors of second derivation of the objective function (Hessian matrix) calculate for an initial model. Limited number of basis vectors obtained in this way defines a new subspace in model parameters space. A new objective function is defined in term of these new parameters and minimized in subspace of original space. As in geophysical inverse problems we need to inverse matrixes that are functions data and geometry of data and model parameters. The matrix inversion in new subspace of the original space will be better conditions due to less dimensionality in the inversion. Since the most significant eigenvectors corresponding the largest eigen values in Singular Value Decomposition ( SVD) of matrixes. Others eigenvectors have less influence in fitting data or lead inversion procedures to local minima. With apply subspace method inversion will be fast and stable against the noise. The efficiency of the method is tested with synthetic and real magnetic data (acquired from Moghan area, north-west of Iran). The results proved fast convergence and stability of inversion against the noise.

Keywords

 
Corbato, C. E., 1965- A least-square procedure for gravity interpretation. Geophysics 30,228-233
Jackson, D. D., 1972- Interpretation of inaccurate and inconsistent data, Geophys  J.R.Astr.Soc. 28, 97-109
Kennett, B. L. N. &  Sambridge, M. S., 1998- Inversion for multiple parameter classes. Geophys.J.int. 135, 304-306
Kennett, B. L. N. & Williamson, P. R., 1988- Subspace methods for large-scale nonlinear inversion, Mathematical Geophysics: a survey of recent development in seismology and geodynamics, Dordrecht. Pp. 139-154.
Kunaratnam, K., 1972- An interactive method for solution of a non linear inverse problem in magnetic interpretation. Geophysical Prospecting 20, 439-447 
Menke, W., 1989- Geophysical data analysis: discrete inverse theory. Academic press Inc.
Mickus, K. L., 1992- Inversion of gravity and magnetic data for lower surface of a 2.5 dimensional sedimentary basin. Geophysical Prospecting 40, 171-191
Minichetti, V., 1983- Simultaneous interactive magnetic and gravity inversion.  Geophysical Prospecting 31, 929-944
Mirzaei, M. & Bredewout, J. W., 1996- 3-D Microgravity data inversion for detecting cavities, European journal of environmental and engineering geophysics, 1, 249-270
Oldenburg, D. W., McGillivary, P. R. & Ellis, R. G., 1993- Generalized subspace method for large-scale inverse problems. Geophys.J.int. 114, 12-20
Oldenburg, D. W., Unsworth, M., 1995- Subspace inversion of electromagnetic data: application to mid-ocean-ridge exploration. Geophys.J.int. 123, 161-168
Pederson, L. B., 1977- Interpretation of potential field data A generalized inverse approach. Geophysical Prospecting 25, 199-230
Sambridge, M. S., 1990- Non-linear arrival time inversion: constraining velocity anomalies by seeking smooth models in 3-D. Geophys.J.int. 102, 635-677
Thurston, J. B., Smith, R. S., Guillon, J., 2002- A multimodel method for depth estimation from magnetic data. Geophysics. 67, 555-561
Wiggins, R. A., 1972- The general linear inverse problem: Implication of surface waves and free oscillation of earth structure. Rev Geophysics and space physics 10, 251-258