Euler Deconvolution of Gravity Data for Structural Mapping and Heat Source Targeting In the Silali Geothermal Field, Kenya

Daniel N. Nyakundi^1, , Erick R. Nyakundi¹ and Solomon Namaswa²

¹Department of Physics, Kisii University, Nairobi, Kenya

²Department of Physics, Multimedia university of Kenya, Nairobi, Kenya

Cite this paper:
Daniel N. Nyakundi, Erick R. Nyakundi and Solomon Namaswa. Euler Deconvolution of Gravity Data for Structural Mapping and Heat Source Targeting In the Silali Geothermal Field, Kenya. Journal of Geosciences and Geomatics. 2025; 13(2):43-51. doi: 10.12691/jgg-13-2-2

Abstract

Gravity data analysis is a widely applied technique for imaging the Earth’s subsurface. This study interprets satellite-derived gravity data over Silali geothermal field to delineate subsurface structures and identify zones favourable for geothermal resource development. Standard reduction procedures were applied to raw gravity data to obtain a complete Bouguer Anomaly (CBA) map. Residual gravity anomalies were extracted using Butterworth and bandpass filters in Oasis Montaj. The anomalies predominantly trend in North-South, Northwest –Southeast directions, forming grabens and horsts that localize magmatic intrusions and geothermal up flow zones. Further enhancement using upward continuation, analytic signal, Power spectrum analysis, and located Euler deconvolution revealed intermediate and deep seated structures. Spectral analysis indicates sediment thickness between 1500 and 5300 metres sufficient to sustain geothermal energy systems. The analytic signal delineates potential intrusions and faults, particularly in the central, South-East, and Western zones. Euler deconvolution solutions (structural index = 0.5, 1.0, 2.0) cluster mainly North and West of the caldera, outlining dense magmatic bodies and fault-controlled features influencing geothermal fluid flow. These findings provide valuable insights for targeting drilling locations in the Silali geothermal system.

Keywords:
Gravity anomalies geothermal delineate power spectrum analysis Euler Deconvolution

This work is licensed under a Creative Commons Attribution 4.0 International License. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

References:

[1]	Pham, L.T., et al., Mapping the structural configuration of the western Gulf of Guinea using advanced gravity interpretation methods. Physics and Chemistry of the Earth, Parts A/B/C, 2023. 129: p. 103341.
[2]	Pham, L.T., et al., Subsurface structural mapping from high-resolution gravity data using advanced processing methods. Journal of King Saud University-Science, 2021. 33(5): p. 101488.
[3]	Nzeuga, A.R., et al., Contribution of advanced edge-detection methods of potential field data in the tectono-structural study of the southwestern part of Cameroon. Frontiers in Earth Science, 2022. 10: p. 970614.
[4]	Ekwok, S.E., et al., Implications of tectonic anomalies from potential field data in some parts of Southeast Nigeria. Environmental Earth Sciences, 2022. 81: p. 1-15.
[5]	Prasad, K., L. Pham, and A. Singh, A novel filter “ImpTAHG” for edge detection and a case study from Cambay Rift Basin, India. Pure and Applied Geophysics, 2022. 179(6): p. 2351-2364.
[6]	Ekwok, S.E., et al., Application of the enhanced horizontal gradient amplitude (EHGA) filter in mapping of geological structures involving magnetic data in southeast Nigeria. Journal of King Saud University-Science, 2022. 34(8): p. 102288.
[7]	Pham, L.T. and K. Prasad, Analysis of gravity data for extracting structural features of the northern region of the Central Indian Ridge. Vietnam Journal of Earth Sciences, 2023. 45(2): p. 147-163.
[8]	Hinze, W.J., et al., New standards for reducing gravity data: The North American gravity database. Geophysics, 2005. 70(4): p. J25-J32.
[9]	Keller, G., Exploration for geothermal energy, in Developments in geophysical exploration methods. 1981, Springer. p. 107-149.
[10]	Simiyu, S.M. and G.R. Keller, An integrated geophysical analysis of the upper crust of the southern Kenya rift. Geophysical journal international, 2001. 147(3): p. 543-561.
[11]	Kearey, P., M. Brooks, and I. Hill, An Introduction to Geophysical Exploration. 2013: Wiley.
[12]	Cooper, G. and D. Cowan, Filtering using variable order vertical derivatives. Computers & Geosciences, 2004. 30(5): p. 455-459.
[13]	Mariita, N.O., An integrated geophysical study of the northern Kenya rift crustal structure: implications for geothermal energy prospecting for Menengai area. 2003: The University of Texas at El Paso.
[14]	Kana, J.D., et al., A review of geophysical methods for geothermal exploration. Renewable and Sustainable Energy Reviews, 2015. 44: p. 87-95.
[15]	Forte, E., et al., Imaging and characterization of a carbonate hydrocarbon reservoir analogue using GPR attributes. Journal of Applied Geophysics, 2012. 81: p. 76-87.
[16]	LaFehr, T., Standardization in gravity reduction. Geophysics, 1991. 56(8): p. 1170-1178.
[17]	Kearey, P., M. Brooks, and I. Hill, An introduction to geophysical exploration. Vol. 4. 2002: John Wiley & Sons.
[18]	Mulwa, J.K., Integrated geophysical study of Lake Bogoria basin, Kenya: Implications for geothermal energy prospecting. 2011, University of Nairobi, Kenya.
[19]	Telford, W.M., L.P. Geldart, and R.E. Sheriff, Applied geophysics. 1990: Cambridge university press.
[20]	Thompson, D., EULDPH: A new technique for making computer-assisted depth estimates from magnetic data. Geophysics, 1982. 47(1): p. 31-37.
[21]	Reid, A.B., et al., Magnetic interpretation in three dimensions using Euler deconvolution. Geophysics, 1990. 55(1): p. 80-91.
[22]	Reid, A.B., J. Ebbing, and S.J. Webb, Avoidable Euler errors–the use and abuse of Euler deconvolution applied to potential fields. Geophysical Prospecting, 2014. 62(5): p. 1162-1168.
[23]	Dunkley, P., et al., The geothermal activity and geology of the northern sector of the Kenya Rift Valley. 1993.
[24]	Omenda, P., et al., Ranking of geothermal prospects in the Kenya rift. Kenya Electricity Generating Company Limited. 2000, internal report.
[25]	Achieng, J., et al., Mapping of hydrothermal minerals related to geothermal activities using remote sensing and GIS: case study of Paka Volcano in Kenyan Rift Valley. International Journal of Geosciences, 2017. 8(5): p. 711-725.
[26]	Omenda, P., Geothermal exploration in Kenya. Short Course VI on Exploration for Geothermal Resources”, organized by UNU-GTP, GDC and KenGen, Lake Naivasha, Kenya, 2010.
[27]	McCall, G., The five caldera volcanoes of the Central Rift Valley, Kenya. Proc. Geol. Soc. London, 1968. 1647: p. 54-59.
[28]	Macdonald, R., et al., Petrogenesis of Silali volcano, Gregory Rift, Kenya. Journal of the Geological Society, 1995. 152(4): p. 703-720.
[29]	Arnason, K., et al. The resistivity structure of high-temperature geothermal systems in Iceland. in Proceedings of the World Geothermal Congress 2000, Kyushu-Tohoku, Japan. 2000.
[30]	Smith, M., et al., Geochronology, stratigraphy and structural evolution of Silali volcano, Gregory Rift, Kenya. Journal of the Geological Society, 1995. 152(2): p. 297-310.
[31]	Mutonga, M., Geothermal geology of Silali volcano in Kenya. Geothermal Res Counc Trans, 2012. 36: p. 943-950.
[32]	Lichoro, C.M., K. Arnason, and W. Cumming, Joint interpretation of gravity and resistivity data from the Northern Kenya volcanic rift zone: Structural and geothermal significance. Geothermics, 2019. 77: p. 139-150.
[33]	Zucchi, M., Faults controlling geothermal fluid flow in low permeability rock volumes: An example from the exhumed geothermal system of eastern Elba Island (northern Tyrrhenian Sea, Italy). Geothermics, 2020. 85: p. 101765.
[34]	Woollard, G.P., The new gravity system—Changes in international gravity base values and anomaly values. Geophysics, 1979. 44(8): p. 1352-1366.
[35]	Blakely, R.J., Potential theory in gravity and magnetic applications. 1996: Cambridge university press.
[36]	Hinze, W.J., Bouguer reduction density, why 2.67? Geophysics, 2003. 68(5): p. 1559-1560.
[37]	Harrison, J. and M. Dickinson, Fourier transform methods in local gravity modeling. Bulletin géodésique, 1989. 63: p. 149-166.
[38]	Nyakundi, E., J. Githiri, and W. Ambuso, Geophysical Investigation of Geothermal Potential of the Gilgil Area. Nakuru County, Kenya Using Gravity Method, 2017.
[39]	Mohammadzadeh Moghaddam, M., et al., Integrated magnetic and gravity surveys for geothermal exploration in Central Iran. Arabian Journal of Geosciences, 2016. 9: p. 1-12.
[40]	Reid, A.B. and J.B. Thurston, The structural index in gravity and magnetic interpretation: Errors, uses, and abuses. Geophysics, 2014. 79(4): p. J61-J66.
[41]	Urrutia-Fucugauchi, J. and J.H. Flores-Ruiz, Bouguer gravity anomalies and regional crustal structure in central Mexico. International Geology Review, 1996. 38(2): p. 176-194.
[42]	Zheng, T., et al., Crust–mantle structure difference across the gravity gradient zone in North China Craton: seismic image of the thinned continental crust. Physics of the Earth and Planetary Interiors, 2006. 159(1-2): p. 43-58.
[43]	Senosy, M., M. Youssef, and M.A. Zaher, Sedimentary cover in the South Western Desert of Egypt as deduced from Bouguer gravity and drill-hole data. Journal of African Earth Sciences, 2013. 82: p. 1-14.
[44]	Archibald, N., P. Gow, and F. Boschetti, Multiscale edge analysis of potential field data. Exploration Geophysics, 1999. 30(2): p. 38-44.
[45]	Roest, W.R., J. Verhoef, and M. Pilkington, Magnetic interpretation using the 3-D analytic signal. Geophysics, 1992. 57(1): p. 116-125.
[46]	Li, J., et al., Crustal tectonic framework of northern Xinjiang and adjacent regions and its formation. Acta Geologica Sinica (English Edition), 2006. 80(1): p. 143-168.