1887
Volume 69, Issue 6
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

A Matlab‐based optimization algorithm is introduced for inverting fault structures from observed gravity anomalies. A convenient graphical user interface is also presented for incorporating the input parameters without any technical complexity to any users. The inversion code uses particle swarm optimization, and all control parameters are tuned initially for faster convergence. There is no requirement of prior choice of an initial model, that is the advantage of using global optimization. The optimization technique is versatile enough to handle any depth‐varying density distributions. The maximum number of iterations and stopping criterion is fixed initially for getting the best optimized solution. The inverted model's output in terms of fault structure, observed and inverted gravity anomalies and dip, and vertex location of fault plane can be viewed in the graphical user interface at the end of the optimization process. The optimization algorithm is applied to different synthetic models with fixed and depth‐varying density contrasts. All synthetic models are further contaminated with white Gaussian noise for sensitivity analysis, and detailed uncertainty appraisal was also performed for the reliability estimation. Finally, the optimization is implemented for fault structure inversion of the Aswaraopet boundary fault, India, and found that the optimized solution provides a good agreement with the previously published literature. Optimized results indicate that this novel optimization approach demonstrates a robust implementation of fault inversion for any depth‐varying density distributions.

Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.13094
2021-06-14
2021-07-30
Loading full text...

Full text loading...

References

  1. Abdelrahman, E. and Essa, K. (2015) Three least‐squares minimization approaches to interpret gravity data due to dipping faults. Pure and Applied Geophysics, 172(2), 427–438.
    [Google Scholar]
  2. Abdelrahman, E., Essa, K. and Abo‐Ezz, E. (2013) A least‐squares window curves method to interpret gravity data due to dipping faults. Journal of Geophysics and Engineering, 10(2), 025003.
    [Google Scholar]
  3. Amjadi, A. and Naji, J. (2013) Application of genetic algorithm optimization and least square method for depth determination from residual gravity anomalies. Journal of Science, Engineering and Technology, 11, 114–123.
    [Google Scholar]
  4. Anderson, N., Essa, K.S. and Elhussein, M. (2020) A comparison study using particle swarm optimization inversion algorithm for gravity anomaly interpretation due to a 2D vertical fault structure. Journal of Applied Geophysics, page 104120.
    [Google Scholar]
  5. Assaad, F.A. (2009) Surface geophysical petroleum exploration methods. In Field Methods for Petroleum Geologists. Springer, pp. 21–23.
    [Google Scholar]
  6. Barbosa, V.C., Menezes, P.T. and Silva, J.B. (2007) Gravity data as a tool for detecting faults: in‐depth enhancement of subtle Almada's basement faults, Brazil. Geophysics, 72(3), B59–B68.
    [Google Scholar]
  7. 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.
    [Google Scholar]
  8. Biswas, A. (2015) Interpretation of residual gravity anomaly caused by simple shaped bodies using very fast simulated annealing global optimization. Geoscience Frontiers, 6(6), 875–893.
    [Google Scholar]
  9. Biswas, A. and Sharma, S.P. (2017) Interpretation of self‐potential anomaly over 2‐D inclined thick sheet structures and analysis of uncertainty using very fast simulated annealing global optimization. Acta Geodaetica et Geophysica, 52(4), 439–455.
    [Google Scholar]
  10. Chai, Y. and Hinze, W.J. (1988) Gravity inversion of an interface above which the density contrast varies exponentially with depth. Geophysics, 53(6), 837–845.
    [Google Scholar]
  11. Chakravarthi, V. and Sundararajan, N. (2004a) Automatic 3‐D gravity modeling of sedimentary basins with density contrast varying parabolically with depth. Computers & Geosciences, 30(6), 601–607.
    [Google Scholar]
  12. Chakravarthi, V. and Sundararajan, N. (2004b) Ridge‐regression algorithm for gravity inversion of fault structures with variable density. Geophysics, 69(6), 1394–1404.
    [Google Scholar]
  13. Chakravarthi, V. and Sundararajan, N. (2005) Invgrafalt: a Fortran code for Marquardt inversion of gravity anomalies of faulted beds with varying density. Computers & Geosciences, 31(10), 1234–1240.
    [Google Scholar]
  14. Eberhart, R. and Kennedy, J. (1995) A new optimizer using particle swarm theory. In MHS'95. Proceedings of the Sixth International Symposium on Micro Machine and Human Science. IEEE, pp. 39–43.
    [Google Scholar]
  15. Ekinci, Y.L., Balkaya, Ç. and Göktürkler, G. (2019) Parameter estimations from gravity and magnetic anomalies due to deep‐seated faults: differential evolution versus particle swarm optimization. Turkish Journal of Earth Sciences, 28(6), 860–881.
    [Google Scholar]
  16. Ekinci, Y.L., Balkaya, Ç., Göktürkler, G. and Turan, S. (2016) Model parameter estimations from residual gravity anomalies due to simple‐shaped sources using differential evolution algorithm. Journal of Applied Geophysics, 129, 133–147.
    [Google Scholar]
  17. Essa, K.S. (2013) Gravity interpretation of dipping faults using the variance analysis method. Journal of Geophysics and Engineering, 10(1), 015003.
    [Google Scholar]
  18. Essa, K.S. and Géraud, Y. (2020) Parameters estimation from the gravity anomaly caused by the two‐dimensional horizontal thin sheet applying the global particle swarm algorithm. Journal of Petroleum Science and Engineering, 107421.
    [Google Scholar]
  19. Essa, K.S., Mehanee, S.A. and Elhussein, M. (2021) Gravity data interpretation by a two‐sided fault‐like geologic structure using the global particle swarm technique. Physics of the Earth and Planetary Interiors, 311, 106631.
    [Google Scholar]
  20. Essa, K.S., Mehanee, S.A., Soliman, K.S. and Diab, Z.E. (2020) Gravity profile interpretation using the R‐parameter imaging technique with application to ore exploration. Ore Geology Reviews, 126, 103695.
    [Google Scholar]
  21. Essa, K.S. and Munschy, M. (2019) Gravity data interpretation using the particle swarm optimisation method with application to mineral exploration. Journal of Earth System Science, 128(5), 123.
    [Google Scholar]
  22. Fernández‐Martínez, J.L., Fernández‐Muñiz, Z., Pallero, J. and Pedruelo‐González, L.M. (2013) From Bayes to Tarantola: new insights to understand uncertainty in inverse problems. Journal of Applied Geophysics, 98, 62–72.
    [Google Scholar]
  23. Fernández Martínez, J.L., Mukerji, T., Garcia Gonzalo, E. and Suman, A. (2012) Reservoir characterization and inversion uncertainty via a family of particle swarm optimizers. Geophysics, 77(1), M1–M16.
    [Google Scholar]
  24. Florio, G. (2020) The Estimation of depth to basement under sedimentary basins from gravity data: review of approaches and the ITRESC method, with an application to the Yucca Flat Basin (Nevada). Surveys in Geophysics, 41(5), 935–961.
    [Google Scholar]
  25. Gao, Y.‐l., An, X.‐h. and Liu, J.‐m. (2008) A particle swarm optimization algorithm with logarithm decreasing inertia weight and chaos mutation. In 2008 International Conference on Computational Intelligence and Security, Vol. 1. IEEE, pp. 61–65.
    [Google Scholar]
  26. Kaftan, İ. (2017) Interpretation of magnetic anomalies using a genetic algorithm. Acta Geophysica, 65(4), 627–634.
    [Google Scholar]
  27. Lelièvre, P.G., Farquharson, C.G. and Hurich, C.A. (2012) Joint inversion of seismic traveltimes and gravity data on unstructured grids with application to mineral exploration. Geophysics, 77(1), K1–K15.
    [Google Scholar]
  28. Li, H.‐R. and Gao, Y.‐L. (2009) Particle swarm optimization algorithm with exponent decreasing inertia weight and stochastic mutation. In 2009 Second International Conference on Information and Computing Science, Vol. 1. IEEE, pp. 66–69.
    [Google Scholar]
  29. Li, Y. and Oldenburg, D.W. (1996) 3‐D inversion of magnetic data. Geophysics, 61(2), 394–408.
    [Google Scholar]
  30. Litinsky, V.A. (1989) Concept of effective density: key to gravity depth determinations for sedimentary basins. Geophysics, 54(11), 1474–1482.
    [Google Scholar]
  31. Liu, S., Hu, X. and Liu, T. (2014) A stochastic inversion method for potential field data: ant colony optimization. Pure and Applied Geophysics, 171(7), 1531–1555.
    [Google Scholar]
  32. Liu, S., Liang, M. and Hu, X. (2018) Particle swarm optimization inversion of magnetic data: field examples from iron ore deposits in China. Geophysics, 83(4), J43–J59.
    [Google Scholar]
  33. Malik, R.F., Rahman, T.A., Hashim, S.Z.M. and Ngah, R. (2007) New particle swarm optimizer with sigmoid increasing inertia weight. International Journal of Computer Science and Security, 1(2), 35–44.
    [Google Scholar]
  34. Murty, B. and Raghavan, V. (2002) The gravity method in groundwater exploration in crystalline rocks: a study in the peninsular granitic region of Hyderabad, India. Hydrogeology Journal, 10(2), 307–321.
    [Google Scholar]
  35. Nagy, D. (1966) The gravitational attraction of a right rectangular prism. Geophysics, 31(2), 362–371.
    [Google Scholar]
  36. Nickabadi, A., Ebadzadeh, M.M. and Safabakhsh, R. (2011) A novel particle swarm optimization algorithm with adaptive inertia weight. Applied Soft Computing, 11(4), 3658–3670.
    [Google Scholar]
  37. Pallero, J., Fernandez‐Martinez, 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.
    [Google Scholar]
  38. Paul, M. and Goodacre, A. (1984) The gravity profile and its role in positioning the edge of a two‐dimensional faulted structure having an arbitrary vertical variation of density. Geophysics, 49(7), 1097–1104.
    [Google Scholar]
  39. Perez, R.l. and Behdinan, K. (2007) Particle swarm approach for structural design optimization. Computers & Structures, 85(19‐20), 1579–1588.
    [Google Scholar]
  40. Pham, L.T., Oksum, E. and Do, T.D. (2018) Gch_gravinv: a matlab‐based program for inverting gravity anomalies over sedimentary basins. Computers & Geosciences, 120, 40–47.
    [Google Scholar]
  41. Portniaguine, O. and Zhdanov, M.S. (2002) 3‐D magnetic inversion with data compression and image focusing. Geophysics, 67(5), 1532–1541.
    [Google Scholar]
  42. Roshan, R. and Singh, U.K. (2017) Inversion of residual gravity anomalies using tuned PSO. Geoscientific Instrumentation, Methods and Data Systems, 6(1), 71.
    [Google Scholar]
  43. Roy, A., Dubey, C.P. and Prasad, M. (2021) Gravity inversion for heterogeneous sedimentary basin with b‐spline polynomial approximation using differential evolution algorithm. Geophysics, 86(3), 1–63.
    [Google Scholar]
  44. Roy, L., Sen, M.K., Blankenship, D.D., Stoffa, P.L. and Richter, T.G. (2005) Inversion and uncertainty estimation of gravity data using simulated annealing: an application over Lake Vostok, East Antarctica. Geophysics, 70(1), J1–J12.
    [Google Scholar]
  45. Shi, Y. and Eberhart, R. (1998) A modified particle swarm optimizer. In 1998 IEEE International Conference on Evolutionary Computation Proceedings. IEEE World Congress on Computational Intelligence . IEEE, pp. 69–73.
    [Google Scholar]
  46. Silva, J.B., Costa, D.C. and Barbosa, V.C. (2006) Gravity inversion of basement relief and estimation of density contrast variation with depth. Geophysics, 71(5), J51–J58.
    [Google Scholar]
  47. Singh, A. (2020) Triangular grid‐based fuzzy cross‐update inversion of gravity data: case studies from mineral exploration. Natural Resources Research, 29(1), 459–471.
    [Google Scholar]
  48. Singh, A. and Biswas, A. (2016) Application of global particle swarm optimization for inversion of residual gravity anomalies over geological bodies with idealized geometries. Natural Resources Research, 25(3), 297–314.
    [Google Scholar]
  49. Stern, T.A. (1978) Gravity Survey of the Taylor Glacier, Victoria Land, Antarctica. Antarctic Research Centre, Victoria University of Wellington.
    [Google Scholar]
  50. Sundararajan, N. and Srinivas, Y. (1996) A modified Hilbert transform and its application to self potential interpretation. Journal of Applied Geophysics, 36(2‐3), 137–143.
    [Google Scholar]
  51. Talwani, M., Worzel, J.L. and Landisman, M. (1959) Rapid gravity computations for two‐dimensional bodies with application to the Mendocino submarine fracture zone. Journal of Geophysical Research, 64(1), 49–59.
    [Google Scholar]
  52. Tikhonov, Andrey N. and Arsenin, V.Y. (1977) Solutions of Ill‐posed Problems. New York: Halsted Press, pp. 1–30.
    [Google Scholar]
  53. Toushmalani, R. (2013) Gravity inversion of a fault by particle swarm optimization (PSO). SpringerPlus, 2(1), 315.
    [Google Scholar]
  54. Wan, L. and Zhang, J. (2019) Analytical solutions of gravity vector and gravity gradient tensor caused by a 2D polygonal body with a 2D polynomial density contrast. Surveys in Geophysics, 1–33.
    [Google Scholar]
  55. Winckel, G.V. (2004) Legendre‐gauss quadrature weights and nodes. MATLAB Central File Exchange, available at https://www.mathworks.com/matlabcentral/fileexchange/4540‐legendre‐gauss‐quadrature‐weights‐and‐nodes. [Accessed 11 May, 2004].
  56. Xin, J., Chen, G. and Hai, Y. (2009) A particle swarm optimizer with multi‐stage linearly‐decreasing inertia weight. In 2009 International Joint Conference on Computational Sciences and Optimization, Vol. 1. IEEE, pp. 505–508.
    [Google Scholar]
  57. Zhou, X. (2008) 2D vector gravity potential and line integrals for the gravity anomaly caused by a 2D mass of depth‐dependent density contrast. Geophysics, 73(6), I43–I50.
    [Google Scholar]
  58. Zhou, X. (2009) General line integrals for gravity anomalies of irregular 2D masses with horizontally and vertically dependent density contrast. Geophysics, 74(2), I1–I7.
    [Google Scholar]
  59. Zhou, X. (2013) Gravity inversion of 2D bedrock topography for heterogeneous sedimentary basins based on line integral and maximum difference reduction methods. Geophysical Prospecting, 61(1), 220–234.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/1365-2478.13094
Loading
/content/journals/10.1111/1365-2478.13094
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): Gravity , Inversion , Modelling and Parameter estimation
This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error