Non-linear stochastic inversion of 2D gravity data using evolution strategy (ES)

نوع مقاله : مقاله پژوهشی‌

نویسندگان

1 موسسه ژئوفیزیک دانشگاه تهران

2 دانشکده علوم زیستی - دانشگاه خوارزمی

چکیده

In the current work, a 2D non-linear inverse problem of gravity data is solved using the evolution strategies (ES) to find the thickness of a sedimentary layer in a deep-water situation where a thick sedimentary layer usually exists. Such problems are widely encountered in the early stages of petroleum explorations where potential field data are used to find an initial estimate of the basin geometry. However, the gravity data are the non-unique problem, and classical deterministic inversions can only offer a single approximation of the solution. Conversely, ES and evolutionary algorithms in general due to their random nature can offer a range of solutions that all fit the data within the acceptable threshold. The inverse problem is formulated as a single objective unconstrained numerical optimization. In evolutionary algorithms, the random nature of the search follows the same rules of the Darwinian biological evolution. Hence, the search is not exhaustive and requires less computational resources compared to Monte Carlo methods. Herein, first, the algorithm is introduced, and a classical synthetic problem is formed and successfully solved in the presence of white Gaussian noise. Then, a two-layered synthetic oceanic crust is formed. ES is successfully tested for solving the formulated under-determined inverse problem of estimating both the base of the sedimentary layer and the crust (i.e., Moho boundary). Finally, using the proposed method, the thickness of the sedimentary layer of the Caspian basin is found along an E-W profile crossing the Caspian Sea. The results have a good agreement with the previous estimates by deep seismic sounding method. The proposed method could be of particular interest because, in deep-water situations, the high water content of sediments, and the expected large thickness of the sediments among other factors make the use of reflective seismic methods unfeasible.

کلیدواژه‌ها


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

Non-linear stochastic inversion of 2D gravity data using evolution strategy (ES)

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

  • Khadij Ghasemi 1
  • Seyed-Hani Motavalli-Anbaran 1
  • Gilda Karimi 2
1 Institute of Geophysics, University of Tehran
2 , Department of Biological Sciences, Kharazmi University
چکیده [English]

In the current work, a 2D non-linear inverse problem of gravity data is solved using the evolution strategies (ES) to find the thickness of a sedimentary layer in a deep-water situation where a thick sedimentary layer usually exists. Such problems are widely encountered in the early stages of petroleum explorations where potential field data are used to find an initial estimate of the basin geometry. However, the gravity data are the non-unique problem, and classical deterministic inversions can only offer a single approximation of the solution. Conversely, ES and evolutionary algorithms in general due to their random nature can offer a range of solutions that all fit the data within the acceptable threshold. The inverse problem is formulated as a single objective unconstrained numerical optimization. In evolutionary algorithms, the random nature of the search follows the same rules of the Darwinian biological evolution. Hence, the search is not exhaustive and requires less computational resources compared to Monte Carlo methods. Herein, first, the algorithm is introduced, and a classical synthetic problem is formed and successfully solved in the presence of white Gaussian noise. Then, a two-layered synthetic oceanic crust is formed. ES is successfully tested for solving the formulated under-determined inverse problem of estimating both the base of the sedimentary layer and the crust (i.e., Moho boundary). Finally, using the proposed method, the thickness of the sedimentary layer of the Caspian basin is found along an E-W profile crossing the Caspian Sea. The results have a good agreement with the previous estimates by deep seismic sounding method. The proposed method could be of particular interest because, in deep-water situations, the high water content of sediments, and the expected large thickness of the sediments among other factors make the use of reflective seismic methods unfeasible.

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

  • non-linear inversion of gravity
  • evolution strategy
  • Caspian Sea
  • deep-water sediments
Back, T., 1996, Evolutionary Algorithms in Theory and Practice: Evolution Strategies, Evolutionary Programming, Genetic Algorithms: Oxford University Press.
Bäck, T., and Schwefel, H. P., 1993, An overview of evolutionary algorithms for parameter optimization: Evolutionary computation, 1(1), 1-23.
Barbosa, V. C. F., Silva, J. B., and Medeiros, W. E., 1997, Gravity inversion of basement relief using approximate equality constraints on depths: Geophysics, 62(6), 1745-1757.
Beyer, H. G., 2013,The Theory of Evolution Strategies: Springer Science and Business Media.
Blum, C., Chiong, R., Clerc, M., De Jong, K., Michalewicz, Z., Neri, F., and Weise, T. , 2012, Evolutionary optimization, in Variants of Evolutionary Algorithms for Real-World Applications, edited, 1-29: Springer.
Bott, M., 1965, The upper mantle beneath Iceland: Geophysical Journal International, 9(2-3), 275-277.
Čuma, M., Wilson, G. A., and Zhdanov, M. S. , 2012, Large-scale 3D inversion of potential field data: Geophysical Prospecting, 60(6), 1186-1199.
Entezar-Saadat, V., Motavalli-Anbaran, S.H., and Zeyen, H., 2017, Lithospheric structure of the Eastern Iranian plateau from integrated geophysical modeling, A transect from Makran to the Turan platform: Journal of Asian Earth Sciences, 138, 357-366.
Farquharson, C. G., and Oldenburg, D. W., 2004, A comparison of automatic techniques for estimating the regularization parameter in non-linear inverse problems: Geophysical Journal International, 156(3), 411-425.
Fogel, D. B., Fogel, L. J., and Atmar, J. W., 1991, Meta-evolutionary programming, paper presented at Signals, systems and computers: Conference record of the twenty-fifth asilomar conference on, IEEE.
Goldberg, D. E., 1989, Genetic algorithms in search, optimization, and machine learning. Reading: Addison-Wesley.
Hinze, W. J., Von Frese, R. R., and Saad, A. H., 2013, Gravity and magnetic exploration: Principles, practices, and applications, Cambridge University Press.
Jackson, J., Priestley, K., Allen, M., and Berberian, M. , 2002, Active tectonics of the south Caspian basin: Geophysical Journal International, 148(2), 214-245.
Jamasb, A., Motavalli-Anbaran, S. H., and Zeyen, H., 2017, Non-linear stochastic inversion of gravity data via quantum-behaved particle swarm optimisation: application to Eurasia–Arabia collision zone (Zagros, Iran): Geophysica Prospecting, 65(S1), 274-294.
Kramer, O., 2014, A Brief Introduction to Continuous Evolutionary Optimization: Springer.
Krasnopolsky, V. M., and Schiller, H., 2003, Some neural network applications in environmental sciences. Part I: forward and inverse problems in geophysical remote measurements: Neural Networks, 16(3), 321-334.
Martins, C. M., Lima, W. A., Barbosa, V. C., and Silva, J. B., 2011, Total variation regularization for depth-to-basement estimate, Part 1—Mathematical details and applications: Geophysics, 76(1), I1-I12.
Mead, J., and Hammerquist C., 2013, χ^2 Tests for the Choice of the Regularization Parameter in Nonlinear Inverse Problems: SIAM Journal on Matrix Analysis and Applications, 34(3), 1213-1230.
Menke, W., 1989, Discrete Inverse Theory.
Mojica, O., and Bassrei, A., 2015, Application of the Generalized Simulated Annealing Algorithm to the Solution of 2D Gravity Inversion of Basement Relief: 3rd Latin American Geosciences Student Conference.
Nagy, D., Papp, G., and Benedek, J., 2000, The gravitational potential and its derivatives for the prism: Journal of Geodesy, 74(7-8), 552-560.
Pallero, J., Fernández-Martínez, J. L., Bonvalot, S., and Fudym, O., 2015, Gravity inversion and uncertainty assessment of basement relief via Particle Swarm Optimization: Journal of Applied Geophysics, 116, 180-191.
Rechenberg, I., 1989, Evolution strategy: Nature’s way of optimization, in Optimization: Methods and applications, possibilities and limitations, edited, 106-126: Springer.
Sambridge, M., and Mosegaard, K., 2002, Monte Carlo methods in geophysical inverse problems: Reviews of Geophysics, 40(3).
Sen, M. K., and Stoffa P. L., 2013, Global Optimization Methods in Geophysical Inversion: Cambridge University Press.
Silva, J. B., Santos, D. F., and Gomes, K. P., 2014, Fast gravity inversion of basement relief: Geophysics, 79(5), G79-G91.
Silva, J. B., Medeiros, W. E., and Barbosa, V. C., 2001, Potential-field inversion: Choosing the appropriate technique to solve a geologic problem: Geophysics, 66(2), 511-520.
Silva, J. S., and Tenreyro, S., 2006, The log of gravity: The Review of Economics and statistics, 88(4), 641-658.
Snopek, K., 2005, Inversion of gravity data with application to density modeling of the Hellenic subduction zone: Ph. D. Thesis, University of Bochum.
Vatankhah, S., Ardestani, V. E., and Renaut, R. A., 2015, Application of the χ2 principle and unbiased predictive risk estimator for determining the regularization parameter in 3D focusing gravity inversion: Geophysical Journal International, 200(1), 265-277.
Zeyen, H. and Pous, J., 1993, 3D joint inversion of magnetic and gravimetric data with a priori information: Geophysical Journal International, 112(2), 244-256.
Zhdanov, M. S., 2002, Geophysical Inverse Theory and Regularization Problems: Elsevier.