Stochastic modeling of gravity data for 2D basement reliefs without regularization coefficient

نوع مقاله : مقاله تحقیقی‌ (پژوهشی‌)

نویسندگان

1 M.Sc Graduate Institute of Geophysics, University of Tehran, Tehran, Iranدانشگاه تهران

2 Assistant Professor, Institute of Geophysics, University of Tehran, Tehran, Iran

چکیده

In this paper, a modified version of strength Pareto evolutionary algorithm SPEA (II) is used as a multi-objective optimization method in gravity data modelling. In this method, a two-dimensional gravity inversion problem is solved by iteratively random creation of forward models. It is shown that it can be used as a fast and effective inversion tool in the depth modelling of two-dimensional layer problems with applications in depth-to-basements, geometry of bedrocks and sedimentary basins modelling cases. Owing to the direct use of the regularization term as a separate objective function, smooth models have a high chance of being selected as final solutions, which makes the results more acceptable and easier to interpret. The most important advantages of this method are that it works independently of the regularization coefficient; thus, there is no need to run the algorithm so many times to find a proper regularization parameter. Furthermore, there is no need to directly deal with inverse formulations, and last but not least, by using a multi-objective algorithm as a global optimization method, convergence to a stable solution does not depend on the initial model, the way classical inversion methods do. For testing the algorithm, a synthetic model is used for layer boundary modelling and to assess the stability of this algorithm, white Gaussian noise is added to the synthetic model. To evaluate the validity of this method, real data from the Recôncavo basin in Brazil is considered for processing and inversion, and the results are compared to the ones from previous studies. All computations have been done in the GNU Octave 5.1.0 environment.

کلیدواژه‌ها

موضوعات


عنوان مقاله [English]

Stochastic modeling of gravity data for 2D basement reliefs without regularization coefficient

نویسندگان [English]

  • Mahmoud Reshadati 1
  • Seyed-Hani Motavalli-Anbaran 2
1 M.Sc Graduate Institute of Geophysics, University of Tehran, Tehran, Iran
2 Assistant Professor, Institute of Geophysics, University of Tehran, Tehran, Iran
چکیده [English]

In this paper, a modified version of strength Pareto evolutionary algorithm SPEA (II) is used as a multi-objective optimization method in gravity data modelling. In this method, a two-dimensional gravity inversion problem is solved by iteratively random creation of forward models. It is shown that it can be used as a fast and effective inversion tool in the depth modelling of two-dimensional layer problems with applications in depth-to-basements, geometry of bedrocks and sedimentary basins modelling cases. Owing to the direct use of the regularization term as a separate objective function, smooth models have a high chance of being selected as final solutions, which makes the results more acceptable and easier to interpret. The most important advantages of this method are that it works independently of the regularization coefficient; thus, there is no need to run the algorithm so many times to find a proper regularization parameter. Furthermore, there is no need to directly deal with inverse formulations, and last but not least, by using a multi-objective algorithm as a global optimization method, convergence to a stable solution does not depend on the initial model, the way classical inversion methods do. For testing the algorithm, a synthetic model is used for layer boundary modelling and to assess the stability of this algorithm, white Gaussian noise is added to the synthetic model. To evaluate the validity of this method, real data from the Recôncavo basin in Brazil is considered for processing and inversion, and the results are compared to the ones from previous studies. All computations have been done in the GNU Octave 5.1.0 environment.

کلیدواژه‌ها [English]

  • Gravity
  • Inversion
  • regularization coefficient
  • multi-objective
  • Genetic
  • basement relief
Akca., I., Gunther, T., Muller-Petke, M., Basokur, A. T., Yaramanci, U., 2014, Joint parameter estimation from magnetic resonance and vertical electric soundings using a multi-objective genetic algorithm: Geophysical Prospecting, 62, 364–376.
Barbosa, V. C. F., Silva, J. B. C., and Medeiros, W. E., 1997, Gravity inversion of basement relief using approximate equality constraints on depths: Geophysics, 62, 1745–1757.
Barbosa, V. C. F., Silva, J. B. C., and Medeiros, W. E., 1999, Gravity inversion of a discontinuous relief stabilized by weighted smoothness constraints on depth: Geophysics, 64(5), 1429-1437.
Blickle, T., and Thiele, L., 1996, A comparison of selection schemes used in evolutionary algorithms: Evolutionary Computation, 4(4), 361–394.
Coello, C. A. C., 1999, A comprehensive survey of evolutionary-based multi-objective optimization techniques: Knowledge and Information Systems, 1(3), 269-308.
Golub, G. H., Heath, M., and Wahba, G., 1979, Generalized cross validation as a method for choosing a good ridge parameter: Technometrics, 21, 215–23.
Hansen, P. C., 1992, Analysis of discrete ill-posed problems by means of the L-curve: SIAM Review.
Hubert, M. K. A., 1948, A line-integral method for computing the gravimetric effect of two-dimensional masses: Geophysics, 13, 215–222.
Lawrence, K. P., and Phillips, R. J., 2003, Gravity/topography admittance inversion on Venus using niching genetic algorithms: Geophys. Res. Lett., 30(19), 1994.
Le˜ao, J. W. D., Menezes, P. T. L., Beltr˜ao, J. F., and Silva, J. B. C., 1996, Gravity inversion of basement relief constrained by the knowledge of depth at isolated points: Geophysics, 61, 1702–1714.
Li, Y., and Oldenburg, D. W., 1999, 3D Inversion of DC resistivity data using an L-curve criterion: 69th Annual International Meeting., Society of Exploration Geophysicists, Expanded Abstracts, 251–254.
Miernik, K., Bogacz, A., Kozubal, A., Danek, T., and Wojdyla, M., 2016, Pareto joint inversion of 2D magnetotelluric and gravity data – towards practical applications: Acta Geophysica, 64(5), 1655–1672.
Miller, B. L., and Golberg, D. E, 1995, Genetic algorithms tournament selection and the effects of noise: Complex Systems, 9(3), 193-212.
Montesinos, F. G., Arnoso, J., and Vieira, R., 2005. Using a genetic algorithm for 3-D inversion of gravity data in Fuerteventura (Canary Islands): International Journal of Earth Sciences: 94, 301–316.
Morozov, V. A., 1966, On the solution of functional equations by the method of regularization: Soviet Mathematics Doklady, 7, 414–417.
Silva, J. B. C., Costa, D. C. L., and Barbosa, V. C. F., 2006, Gravity inversion of basement relief and estimation of density contrast variation with depth: Geophysics, 71, J51–J58.
Silverman, B. W., 1986, Density estimation for statistics and data analysis: Chapman and Hall.
Telford, W. M., Geldart, L. P., and Sheriff, R. E., 1990, Applied Geophysics: Cambridge University Press.
Talwani, M., Worzel, J. L., and Landisman, M., 1959, Rapid gravity computations for two-dimensional bodies with application to the Mendocino submarine fracture zone: J. . Geophys. Res., 64, 49-59.
Zhang, J., Wang, C., Shi, Y., Cai, Y., Chi, W., Dreger, D., Cheng, W., and Yuan, Y., 2004, Three-dimensional crustal structure in central Taiwan from gravity inversion with a parallel genetic algorithm: Geophysics, 69, 917–924, DOI:10.1190/1.1778235.
Zitzler, E., Laumanns, M., and Thiele, L., 2001, SPEA2: Improving the Strength Pareto Evolutionary Algorithm: TIK-Report, 103.
Zitzler, E., and Thiele, L., 1999, Multi-objective evolutionary algorithms: A comparative case study and the strength Pareto approach: IEEE Transactions on Evolutionary Computation, 3(4), 257–271.