Geosciences

p-ISSN: 2163-1697    e-ISSN: 2163-1719

2019;  9(1): 1-7

doi:10.5923/j.geo.20190901.01

 

Least-Square Geothermal Models in Crustal Rocks with Non-Radiogenic and Radiogenic Heat Generation Sources in the Niger Delta Basin

G. I. Alaminiokuma1, E. S. Osegbowa2

1Department of Earth Sciences, Federal University of Petroleum Resources, Effurun, Effurun, Nigeria

2Department of Physics, Rivers State University of Science & Technology, Nkpolu-Oroworukwo, Nigeria

Correspondence to: G. I. Alaminiokuma, Department of Earth Sciences, Federal University of Petroleum Resources, Effurun, Effurun, Nigeria.

Email:

Copyright © 2019 The Author(s). Published by Scientific & Academic Publishing.

This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/

Abstract

Geothermal models in crustal rocks with non-radiogenic and radiogenic heat generation sources were developed using temperature data from Continuous Temperature Logs in 26 wells randomly distributed across the Niger Delta Basin. Data analyses and interpretation by the method of Least-squares generated results which show that Ground Surface Temperature (GST) ranged between 5.4 and 53.5°C with a mean of 27.98°C; Geothermal Gradient (Gg) ranged between 0.37 to 3.84°C/100m with a mean of 2.14°C/100m. The mean geothermal distribution functions obtained for both the non-radiogenic and radiogenic heat generation models in the Niger Delta are T=27.983+0.0214y and T=27.983+0.0214y-2.64x10-10y2 respectively. This is very vital in effectively and accurately providing geotemperature information essential in estimating the earth’s heat flow; developing the subsidence history including the formation, maturity and migration of hydrocarbon in sedimentary basins, and searching for deep-seated geologic structures favourable to hydrocarbon accumulation.

Keywords: Least-square, Geothermal, Non-Radiogenic, Radiogenic, Crustal Rocks and Niger Delta

Cite this paper: G. I. Alaminiokuma, E. S. Osegbowa, Least-Square Geothermal Models in Crustal Rocks with Non-Radiogenic and Radiogenic Heat Generation Sources in the Niger Delta Basin, Geosciences, Vol. 9 No. 1, 2019, pp. 1-7. doi: 10.5923/j.geo.20190901.01.

1. Introduction

The present efforts of exploring deeper prospects towards recovering more untapped hydrocarbon in the Niger Delta Basin require accurate data analyses and in-depth knowledge of geotemperature distribution. Temperature distribution in the earth’s subsurface is either of radiogenic source or of formation source or both [1]. Geotemperature distribution in the Continental Crust is governed both by either the conductive heat loss to the surface or heat generated internally by the decay of radioactive isotopes in rocks and heat that flows upward from sub-continental mantle [2].
Several studies have been conducted using the simple gradient method, statistical and geostatistical analyses, thermal resistance method, simple linear increase of temperature with depth among others to determine the geothermal characteristics in different parts of the Niger Delta Basin. A brief review of some of these studies related to geotemperature and its applications and variations are given in Table 1.
Table 1. A Review of Studies Related to Geotemperatures
This research is focused on generating geothermal models considering the non-radiogenic and radiogenic heat generation sources which contributes a significant amount of heat to the surface heat flow in the Niger Delta Basin and employing the Least-Squares method to fit the best curves, compute sets of geotemperature data using Normal equations to estimate the Ground Surface Temperature, determine the Geothermal Distribution Functions, predict the vertical trend of geothermal gradient and examine the effects of radiogenic heat generation on the thermal parameters in the Niger Delta Basin. This has the advantage of effectively and accurately providing geotemperature information essential in estimating the earth’s heat flow; developing sedimentary basins, its subsidence history including the formation, maturity and migration of hydrocarbon; searching for deep-seated geologic structures (faults, anticlines) favourable to hydrocarbon accumulation in oil exploration.

2. Location and Geology of Niger Delta

The area of study is the Niger Delta sedimentary Basin (Figure 1). It is situated on the Continental margin of the Gulf of Guinea in Equatorial West Africa between Latitudes 3 °N and 6 °N and Longitude 5 °E and 8 °E.
Figure 1. Map of the Niger Delta Sedimentary Basin
Early studies of the Niger Delta Basin [13] subdivided the basin into three formations namely the Akata (base), the Agbada (middle), and the Benin (top) Formations (Figure 2). These Formations represent the prograding depositional facies that are distinguished mostly on the basis of their sand-shale ratios.
Figure 2. Map showing the stratigraphy of Niger Delta
The Akata Formation comprises mostly marine shale sequences and lowstand turbidite sand. The Agbada Formation comprises alternating sequences of sandstone and shale and is interpreted as a cyclic parallic sequence comprising marine and fluvial deposits. The Benin Formation comprises continental, massive sands with clay intercalations ([4]; [13]; [14]; [15]).

3. Methodology

3.1. Field Data

Geotemperature data from Continuous Temperature Logs in 26 wells randomly distributed across the Niger Delta Basin were used in this study. These data are reliable since the wells have attained thermal equilibrium after some period of drilling [7].

3.2. Geothermal Model Assumptions

A homogenous earth model [16] is assumed and the dependence of temperature variation on lithology is neglected so as to observe the variations of temperature with depth ([7]; [17]). The earth is modelled as an isotropic slab to the “Moho” depth of about 36 km to allow for analytic solutions to the Normal equations governing heat conduction problems. The earth’s crust approximates a plane due to its negligible thickness compared with the thickness of the entire earth. This is a 1-layer problem and 1-dimensional heat flow by conduction is assumed due to the solid nature of the earth’s crust.

3.3. Geothermal Model Equations

3.3.1. Non-Radiogenic Heat Generation Case
For a 1-D steady-state heat conduction through an earth crust that approximates a slab without radiogenic heat generation, external heat penetration and shielded from topographic as well as seasonal changes, the diffusion equations is:
(1)
Given to the great age of the crust, time-dependent effects are generally neglected [2].
Applying the boundary conditions (at y = 0):
(2a)
(2b)
Then qs = -qs because heat is assumed to be flowing out of the earth. Solving equation (1) and applying the boundary conditions, we obtain:
(3)
Where:
(4)
The geotemperature distribution function now becomes:
(5)
Where: Ts = Ground Surface Temperature (GST) (°C), qs = Heat flow (W/m2), k = Thermal Conductivity (W/m°C) and y = Vertical Depth (m).
Equation (5) defines temperature at any depth, y in the earth for a homogenous earth model without radiogenic heat generation. The Least Square equation for this case is given as:
(6)
With the following Normal equations:
(7a)
(7b)
Comparing equations (5) and (6), we have: and .
3.3.2. Radiogenic Heat Generation Case
If the radiogenic heat (H) is generated per unit volume per unit time, then the heat conduction equation becomes:
(8)
Applying the boundary conditions (equations 2a and 2b) and solving equation (8), we obtain:
The geotemperature distribution function:
(9)
The geotemperature gradient function:
(10)
Equation (9) defines temperature at any depth, y in an earth model with radiogenic heat generation. The Least Square equation is given as:
(11)
With the following Normal equations:
(12a)
(12b)
(12c)
Comparing equations (9) and (11), we have: , and .
Solutions to the Normal equations above are obtained by Matrix equations:
(13)
Where B, A and x are matrices then:
(14)
Matrix x is that of required solutions a0, a1, and a2 of the sets of Normal equations. Matrices A and B are developed using in the Normal equations.
Solutions to these equations are obtained by solving the Matrix Algebra using MS Excel package.

4. Results and Discussion

4.1. Ground Surface Temperature (GST)

As stated by [18], the ground surface is an important layer which temperature forms an important boundary condition for the determination of geothermal gradient, a parameter which is useful in evaluating the thermal energy resources of any region. Table 2 shows the ground surface temperature (GST) for the 26 wells in this research.
Table 2. Ground Surface Temperatures (GST) and Least-Square Geothermal Distribution Functions for Crustal Rocks with Non-Radiogenic and Radiogenic Heat Generation for the 26 Wells in the Study Area
The GST ranges between 5.369°C for well 7 and 53.350°C for well 6 for both the non-radiogenic and radiogenic models with a mean value of 27.983°C computed for the Niger Delta. This is in agreement with an earlier value of 27°C predicted for the southern Nigeria sedimentary Basin by [3].

4.2. Geothermal Distribution Functions for Non-Radiogenic and Radiogenic Sources

Table 2 also shows the geothermal distribution functions for crustal rocks with non-radiogenic and radiogenic heat generation sources for the 26 wells studied. The Least-Square equations indicate linear functions for the non-radiogenic heat generation model and quadratic functions for the radiogenic heat generation model. The mean geotemperature distribution functions for wells 1-26 are computed respectively as:
For Non-radiogenic Model:
(15)
For Radiogenic Model:
(16)

4.3. Geotemperature Gradient of the Niger Delta with Radiogenic Heat Generation

Table 3 shows the geotemperatures and geothermal gradients for depths extrapolated to the “Moho” in the Niger Delta Basin using the mean geotemperature functions (equations 15 and 16) with Non- Radiogenic and Radiogenic heat generations respectively for wells 1-26.
Table 3. Mean Geotemperatures and Geothermal Gradient Functions to the “Moho” in the Niger Delta Basin
The ‘Moho’ depth (36000 m) mean temperature for wells 1-26 are 813.503°C for the non-radiogenic model and 813.161°C for the radiogenic model. The trend suggests a distributive pattern of geotemperature for the Crust in the Basin. The effect of radiogenic heat generation is pronounced at great depth and negligible at shallow depth. Geothermal gradient is observed to be constant for the non-radiogenic model but decreases with depth for the radiogenic heat model for all the wells studied. The mean geothermal gradient is constant for wells 1-26 at a value of 2.182 (°C/100m) for the non-radiogenic case but decreases with depth from 2.182 to 2.181 (°C/100m) for the radiogenic model.
Figures 3 (a and b) are plots of mean geotemperature (°C) versus depth (m) for wells 1-26 for the non-radiogenic and radiogenic models respectively.
Figure 3a. Geotemperature (°C) versus Depth (m) for the Non-Radiogenic Heat Generation
Figure 3b. Geotemperature (°C) versus Depth (m) for the Radiogenic Heat Generation
The curves approximate straight lines due to the insignificant effects of the radiogenic heat sources on the geotemperature gradient. As a result of this, the Moho surface temperatures for both models are similar.
Geotemperature gradient is observed to decrease linearly with depth for the radiogenic heat generation model since the third term of the radiogenic heat generation model equation is infinitesimally small and negligible. This is observed in the linear nature of the graph of radiogenic heat generation model irrespective of the quadratic nature (Figure 4).
Figure 4. Geothermal Gradient for the Radiogenic Heat Generation
The geotemperature gradient obtained for the non-radiogenic heat generation model is comparable to the results by [4]; [5]; [6]; [8]; [9]; and [12]. The agreement of the results in this study with those of other researchers, irrespective of the different methods used, further supports the low rates of radiogenic heat generation in the Niger Delta Basin obtained in this research.

5. Conclusions

The application of the Least Square method for the computations of ground surface temperature, geothermal gradients and generation of geotemperature distribution functions for the non-radiogenic and radiogenic models is suitable for accurate geotemperature analyses since it gives very close results to those of other researchers using different methods in the Niger Delta.
The radiogenic heat generation source contributes two significant effects on the subsurface temperature distribution in the region. The low value of radiogenic heat generation per unit thermal conductivity accounts for the small difference between the linear and non-linear Least Square geotemperature model values. The radiogenic model also shows that geotemperature gradient decreases infinitesimally with depth. This can be because of the concentration of radiogenic heat sources with depth based on the stratigraphy in the Niger Delta [13]. The increase in shaliness in the Basin can be viewed as decreasing with depth. This is substantiated by the general trend of decrease in radiogenic heat sources from Crust to the Mantle [1].
It is obvious from the foregoing that the geothermal models of the Niger Delta crust is not significantly affected by the radiogenic heat production of the region.

6. Recommendations

The geothermal models developed for the Niger Delta Basin crustal rocks considering radiogenic heat generation and employing the Least Squares method provide new frontiers for research in hydrocarbon and geothermal energy resources appraisal in the Basin. A detailed experimental work using Temperature Logs, Gamma Ray Logs (to obtain radiogenic heat generation) and core samples (to determine thermal conductivity) should be conducted in the study area to obtain results for comparison with the approach in this study.
Further research should be devoted to the study of the effects of other thermodynamic conditions apart from radiogenic heat generation on geotemperature distribution in the Niger Delta and other sedimentary Basins.

ACKNOWLEDGEMENTS

The authors wish to thank Shell Petroleum Development Company Nigeria for making available the data used in this research and to the Department of Petroleum Resources for the permission to use the data.

References

[1]  W. Lowrie, Fundamentals of Geophysics, 1st Edition, Cambridge University Press, Cambridge, UK, Pp 178-202, 1997.
[2]  L. T. Donald, and S. Gerald, Geodynamics, John Wiley and Sons Inc. U.S.A. 135-194, 1982.
[3]  Nwachukwu, S. O., 1976, Approximate Geothermal Gradients in Niger Delta Sedimentary Basin. AAPG Bull, 60: 1073-1077.
[4]  Evamy, B. D., Haremboure, J., Kamerling, P., Knaap, W. A. and Molloy F. A., 1978, Hydrocarbon Habitat of Tertiary Niger Delta. Am. Associat. Petrol. Geol. Bull., 62: 277-298.
[5]  Chukwueke, C., Thomas, G. and Delfraud, J., 1992, Sedimentary Processes, Eustatism, Subsidence and Heat Flow in the Distal Part of the Niger Delta. Bull. Centres Rech. Explor. Prod. Elf-Acquitaine, 16: 137-186.
[6]  Uko, E. D. 1996, Thermal Modelling of the Northern Niger Delta unpublished Ph.D. Thesis, University of Science and Technology, Port Harcourt.
[7]  Akpabio, I. O. and Ejedawe, J. E., 2001, Temperature Variation in the Niger Delta Subsurface from Continuous Temperature Logs. Global Journal of Pure and Applied Sciences. Vol. 7(1), 137-141.
[8]  Akpabio, I. O., Ejedawe, J. E., Ebeniro, J. O. and Uko, E. D., 2003, Geothermal Gradients in the Niger Delta Basin from Continuous Temperature Logs. Global J. Pure Applied Sci., 9: 265-271.
[9]  Uko, E. D. and Eze, C. L., 2003, Estimation of Geothermal Gradients in the North-Western Niger Delta, Nigeria. Journal of Applied Sciences. 6(3), 3814-3823.
[10]  Emujakporue, G. O. and Ekine, A. S., 2014, Determination of geothermal gradient in the Eastern Niger Delta Sedimentary Basin from bottom hole temperatures. J. Earth Sci. Geotechnical Eng., 4: 109-114.
[11]  Odumodu, C. F. R. and Mode, A. W., 2016, Geothermal Gradients and Heat Flow Variations in Parts of the Eastern Niger Delta, Nigeria. Journal Geological Society of India. Vol. 88, Pp.107-118.
[12]  Emujakporue, G. O. and Nwosu, L. I., 2017, Spatial Variation Modeling of Geothermal Gradient and Heat Flow in Eastern Parts of Niger Delta Sedimentary Basin, Nigeria. Physical Science International Journal 14(1): 1-13.
[13]  Short, K. C. and Stauble, A. J., 1967, Outline of Geology of Niger Delta. AAPG Bulletin volume 51, pp. 761- 779.
[14]  Weber, K. J., 1987, Hydrocarbon Distribution Patterns in Nigerian Growth Fault Structures Controlled by Structural Style and Stratigraphy: Journal of Petroleum Science and Engineering, v. 1, p. 91-104.
[15]  Reijers, T. J. A., Petters, S. W., and Nwajide, C. S., 1997, The Niger Delta Basin, in Selley, R.C., ed., African Basins-Sedimentary Basin of the World 3: Amsterdam, Elsevier Science, pp. 151-172.
[16]  C. M. R. Fowler, The Solid Earth, Cambridge University Press, Cambridge, Pp 505, 1990.
[17]  Uko, E. D., Amakiri, A. R. C. and Alagoa, K. O., 2002, Effects of Lithology on Geothermal Gradient on the Southeast Nigeria Delta, Nigeria. Global J. Pure Applied Sci., 8: 325-338.
[18]  S. K. M. Ali and D. M. Orazulike, Estimation of Ground Surface Temperature from Meteorological Records, from Bauchi, North East, Nigeria. Atil Press Nigeria, 276-299, 2003.