**Radiogenic Heat Generation in Western Australia — Implications for Geothermal Energy**

Mike F. Middleton

[7] Nasruddin, M., Alhamid, I., Daud, Y., Surachman, A., Sugiyono, A., Aditya, H.B., Mahila, T.M.I. 2016. Potential of geothermal energy for electricity generation in Indo‐

[8] Resosudarmo, B.P., Alisjahbana, A., Nurdianto, D.A. 2010., "Energy Security in Indo‐ nesia, ANU Working Papers in Trade and Development No. 2010/08, Canberra. [9] Fiscal Policy Office of the Ministry of Finance) & World Bank. 2011. "Phase 1 assess‐ ment of emissions – key findings," viewed 16 March 2011, at (http://www.esmap.org/ esmap/sites/esmap.org/files/factsheet04.pdf). IGA-International geothermal associa‐

[10] Fadillah, R.D. 2011. "More German companies eye renewable energy sector," *The Ja‐ karta Post*, 25 October, viewed 16 March 2011, at (http://www.thejakartapost.com/

news/2011/10/25/moregerman-companies-eye-renewable-energy-sector.html).

Emissions), The Secretary of the Cabinet of the Republic of Indonesia, Jakarta.

[12] The Secretary of the Cabinet of the Republic of Indonesia 2011. "Cabinet Secretary: President Yudhoyono signed presidential decree on national action plan to reduce GHG emissions," Press Release, 26 September 2011, the Secretary of Cabinet of the

[13] WWF., 2012. Igniting the ring of fire: A vision for developing Indonesia's geothermal

[14] Tanoto, Y., Wijaya, M.E. 2011. Economic and Environmental Emissions Analysis in Indonesian Electricity Expansion Planning: Low-rank Coal and Geothermal Energy Utilization Scenarios, *International Journal of Renewable Energy Research*, Vol. 1, 2011.

[16] The World Bank Energy Sector Management Assistance Program. 2010. Geothermal handbook: planning and financing geothermal power generation. Washington DC:

[15] de Wilde, A. 2010. Accelerating Geothermal Development in Indonesia. Jakarta.

[11] President of the Republic of Indonesia. 2011. "Peraturan Presiden no. 61 tahun 2011 tentang Rencana Aksi Nasional Penurunan Emisi Gas Rumah Kaca" (Presidential Decree No. 61 of 2011 on National Action Plan for the Reduction of Greenhouse Gas

nesia: A review. *Renew Sustain Energy Rev*, 53, 733-740.

tion, http://www.geothermal-energy.org

Republic of Indonesia, Jakarta.

The World Bank. 2010.

power.

48 Advances in Geothermal Energy

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/61963

#### **Abstract**

The chapter reviews heat generation in crystalline rocks and influences on overlying sedi‐ mentary basins in Western Australia (WA). Regions of elevated thorium and uranium will cause elevated heat generation, which in turn can cause elevated heat flow. Western Australia hosts several large sedimentary basins with the potential for hot sedimentary aquifers (HSAs). These include the Perth, Carnarvon, and Canning basins. Parts of these basins are underlain by crystalline rocks that contain high levels of heat-generating ele‐ ments, such as uranium, thorium, and potassium. Also, the Pilbara Craton, which con‐ tains both sedimentary and crystalline rocks, that entertains a number of active mines, which may benefit from geothermal energy, is investigated. Further, the southern part of the Perth Basin (Vasse Shelf), which is underlain by crystalline rocks with elevated con‐ centrations of thorium and uranium, is shown to possess higher than usual temperatures. From observations, and geothermal modeling, it is concluded that the Perth Basin has a high potential for medium- to low-temperature geothermal energy developments. In oth‐ er parts of Western Australia, the Carnarvon Basin has elevated temperatures in artesian groundwater. Heat flow in the Canning Basin is briefly reviewed; this basin has some ge‐ othermal potential, but it is far from the major population centers.

**Keywords:** Radiogenic heat, Uranium, Thorium, Potassium, heat flow

## **1. Introduction**

The information presented in this study reflects the view in Western Australia to move toward mechanisms for sustainable energy into the future. The Western Australian legislative framework supports this view, and continuous studies have been carried out since 2008 to implement this vision into the future. The initial formal Western Australian governmental view was to use geothermal energy for electricity generation, and legislation was formulated

© 2016 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

to accommodate this. By late 2014, it became apparent that geothermal energy in Western Australia was taking a different path. The chapter focuses, both formally published works and less formal studies, on radiogenic heat generation within rocks in Western Australia and their contribution to heat flow, elevated temperatures, and geothermal energy potential. It adopts the approach that heat flow from basement (essentially crystalline) rocks beneath sedimentary basins will elevate temperatures within the sedimentary basin, and thus create a natural reservoir of hot fluids in the sediments that can be used for geothermal energy purposes. These are referred to as hot sedimentary aquifers.

The connection between radioactive elements in the earth and heat generation has been known for many years. Among many, the following workers have developed this science [1–3], and in specific reference to Australia [4–7]. About this period, an increasing awareness was growing for the potential need for geothermal energy in the Australian energy mix, and a landmark book was published by Beardsmore and Cull [7]. Studies were undertaken by Geoscience Australia, the national geoscience survey organization in Canberra [8–10]. In Western Australia, one of the significant early workers to recognize that geothermal energy may be viable in the state was by Bestow [11]. In his study, Bestow [11] recognized parallels of Western Australian basins to the Paris Basin and the need to achieve temperatures in the vicinity of 65–85 °C. At this early stage, Bestow also clearly recognized that the geothermal potential in Western Australia was best based around the concept of hot sedimentary aquifers, rather than either the engineered geothermal systems (EGSs) being advocated elsewhere in Australia [12,13] or the more commonly exploited volcanogenic heat sources.

Early work on heat flow and radiogenic heat production was carried out by Jaeger [4], Sass et al. [5], Middleton [6], Cull and Denham [14], and Cull [15]. These studies identified the background to understanding heat flow and heat generation in Western Australia. Perhaps rather more by serendipity than by intention, the study by Jaeger [4] identified a rather unique site in the "Wheat Belt" of Western Australia (Figure 1). This site was in a granite quarry near the small town of Doodlakine, and returned a heat generation value of about 21.9 heat generation units (hgu; 1 hgu = 0.418 µW m−3) at the surface and 21.2 hgu as an average of 30 samples from an associated bore hole. This translates to 9.15 µW m−3 in SI units. At the time, this was the largest heat generation value published for a Western Australian granitoid (a term used herein to refer to a felsic igneous rock with largely granitic mineralogy). More recent work on heat generation in Western Australia was carried out by Middleton [16], Middleton and Stevens [17], and Middleton et al. [18]. During these studies, the Doodlakine site was revisited and now initially reported; surface measurements of heat generation of the granite near the original measurements published by Jaeger [4] were made, using an RS-125 gamma-ray scintillometer (see below for description of the technique). Twelve measurements were made on the granite outcrop, near where Jaeger (1970) carried out his investigation, and these yielded a mean heat generation of 8.5 µW m−3 (20.3 hgu used by Jaeger [4]). This is a good agreement, given slightly different assumptions made for conversion of uranium (U), thorium (Th), and potassium (K) concentrations to heat generation, made by Jaeger [4] as compared with Middleton [16].

Radiogenic Heat Generation in Western Australia — Implications for Geothermal Energy http://dx.doi.org/10.5772/61963 51

to accommodate this. By late 2014, it became apparent that geothermal energy in Western Australia was taking a different path. The chapter focuses, both formally published works and less formal studies, on radiogenic heat generation within rocks in Western Australia and their contribution to heat flow, elevated temperatures, and geothermal energy potential. It adopts the approach that heat flow from basement (essentially crystalline) rocks beneath sedimentary basins will elevate temperatures within the sedimentary basin, and thus create a natural reservoir of hot fluids in the sediments that can be used for geothermal energy purposes. These

The connection between radioactive elements in the earth and heat generation has been known for many years. Among many, the following workers have developed this science [1–3], and in specific reference to Australia [4–7]. About this period, an increasing awareness was growing for the potential need for geothermal energy in the Australian energy mix, and a landmark book was published by Beardsmore and Cull [7]. Studies were undertaken by Geoscience Australia, the national geoscience survey organization in Canberra [8–10]. In Western Australia, one of the significant early workers to recognize that geothermal energy may be viable in the state was by Bestow [11]. In his study, Bestow [11] recognized parallels of Western Australian basins to the Paris Basin and the need to achieve temperatures in the vicinity of 65–85 °C. At this early stage, Bestow also clearly recognized that the geothermal potential in Western Australia was best based around the concept of hot sedimentary aquifers, rather than either the engineered geothermal systems (EGSs) being advocated elsewhere in

Early work on heat flow and radiogenic heat production was carried out by Jaeger [4], Sass et al. [5], Middleton [6], Cull and Denham [14], and Cull [15]. These studies identified the background to understanding heat flow and heat generation in Western Australia. Perhaps rather more by serendipity than by intention, the study by Jaeger [4] identified a rather unique site in the "Wheat Belt" of Western Australia (Figure 1). This site was in a granite quarry near the small town of Doodlakine, and returned a heat generation value of about 21.9 heat generation units (hgu; 1 hgu = 0.418 µW m−3) at the surface and 21.2 hgu as an average of 30 samples from an associated bore hole. This translates to 9.15 µW m−3 in SI units. At the time, this was the largest heat generation value published for a Western Australian granitoid (a term used herein to refer to a felsic igneous rock with largely granitic mineralogy). More recent work on heat generation in Western Australia was carried out by Middleton [16], Middleton and Stevens [17], and Middleton et al. [18]. During these studies, the Doodlakine site was revisited and now initially reported; surface measurements of heat generation of the granite near the original measurements published by Jaeger [4] were made, using an RS-125 gamma-ray scintillometer (see below for description of the technique). Twelve measurements were made on the granite outcrop, near where Jaeger (1970) carried out his investigation, and these yielded a mean heat generation of 8.5 µW m−3 (20.3 hgu used by Jaeger [4]). This is a good agreement, given slightly different assumptions made for conversion of uranium (U), thorium (Th), and potassium (K) concentrations to heat generation, made by Jaeger [4] as compared with

Australia [12,13] or the more commonly exploited volcanogenic heat sources.

are referred to as hot sedimentary aquifers.

50 Advances in Geothermal Energy

Middleton [16].

**Figure 1.** WA map showing tectonic provinces. Note that the Doodlakine location is halfway between Perth and Kal‐ goorlie.

Uranium distribution, as with Th distribution, in Western Australia is poorly known at depth. However, a reasonable knowledge of surface U distribution is known from airborne radio‐ metric data compiled in the "uranium merge" and "thorium merge" maps of the state produced by the Geological Survey of Western Australia [19]. The U and Th merge maps are shown in Figure 2. It must be recognized that these maps have limitations because of the attenuation of radiation, due to many effects, for example, the distance from surface radioactive source rocks to the sensor in the aircraft, time of day of acquisition, the limitations (because of different instrumentation and calibration) in the merge process for various survey datasets, and the difference in surface geology across the state [20,21]. From the maps, it is observed that Th concentration (parts per million, ppm) is generally higher than the U concentration. Thorium to uranium ratio (Th/U) has been reviewed by Middleton et al. [18] for various surface locations (Table 1). It should be noted that the surface U and Th concentrations from airborne measurements are commonly lower than surface measurements of concentration, because of the environmental effects commonly observed in deriving U and Th from the surface "merge" maps [20,21]. From the data in Table 1, the Th/U ratio showed a variation between 2.0 and 13.0 for outcrop rocks. Uranium and thorium are often concentrated in deeply emplaced felsic igneous (granitoid) rocks. However, these elements can also be concentrated in sedimentary rocks. In recognizing this, some surface measurements of U and Th in sedimentary basins may be higher than normal due to sedimentary depositional and enrichment processes.

**Figure 2.** Map with surface U and Th surface distribution in WA. These are based on "merge maps" produced by Brett [19] for the Geological Survey of Western Australia.

Radiogenic Heat Generation in Western Australia — Implications for Geothermal Energy http://dx.doi.org/10.5772/61963 53


Uranium distribution, as with Th distribution, in Western Australia is poorly known at depth. However, a reasonable knowledge of surface U distribution is known from airborne radio‐ metric data compiled in the "uranium merge" and "thorium merge" maps of the state produced by the Geological Survey of Western Australia [19]. The U and Th merge maps are shown in Figure 2. It must be recognized that these maps have limitations because of the attenuation of radiation, due to many effects, for example, the distance from surface radioactive source rocks to the sensor in the aircraft, time of day of acquisition, the limitations (because of different instrumentation and calibration) in the merge process for various survey datasets, and the difference in surface geology across the state [20,21]. From the maps, it is observed that Th concentration (parts per million, ppm) is generally higher than the U concentration. Thorium to uranium ratio (Th/U) has been reviewed by Middleton et al. [18] for various surface locations (Table 1). It should be noted that the surface U and Th concentrations from airborne measurements are commonly lower than surface measurements of concentration, because of the environmental effects commonly observed in deriving U and Th from the surface "merge" maps [20,21]. From the data in Table 1, the Th/U ratio showed a variation between 2.0 and 13.0 for outcrop rocks. Uranium and thorium are often concentrated in deeply emplaced felsic igneous (granitoid) rocks. However, these elements can also be concentrated in sedimentary rocks. In recognizing this, some surface measurements of U and Th in sedimentary basins may

52 Advances in Geothermal Energy

be higher than normal due to sedimentary depositional and enrichment processes.

**Figure 2.** Map with surface U and Th surface distribution in WA. These are based on "merge maps" produced by Brett

[19] for the Geological Survey of Western Australia.

**Table 1.** Western Australian RS-125 measurement campaign for U, Th, K & Ao for the years 2011 and 2015.

Potassium usually has a different genesis to U and Th. It is a very commonly occurring element in the earth's crust. The abundance (concentration) of K in most igneous rocks is commonly observed to be about 2–6%, and it also occurs abundantly in many sedimentary rocks. Beardsmore and Cull [7] have discussed the specific abundance of the radioactive isotope (K-40) within normal potassium distributions in rocks. The surface K "merge" distribution map [19] can be misleading, as the K concentrations tend to follow recent sedimentary depositional features, such as seasonally dry riverbeds. In such deposits, the time of year of that an airborne radiometric survey was acquired has implications for interpretation of the observed K concentration [20–22]. Further, as will be seen in the following sections, the contribution of K to radiogenic heat generation is usually less than the contribution of U and Th [7,16,23].

In this review, surface-based geophysical or laboratory measurements of U, Th, and K are emphasized rather than airborne measurements. This is because surface-based (on-theground) measurements of outcropping rocks tend to be more directly indicative of the particular rock type being investigated, rather than a "bulk" estimate of the chemistry of the rocks being sampled by an airborne survey. This will be discussed further below.

## **2. Theory**

This section reviews the basic, but well-known, concepts to understand heat generation, heat flow, and simple temperature distributions with the earth's crust.

## **2.1. Heat flow and temperature due to heat generation in a radiogenic layer**

Heat flow at the surface on the earth (Qs) over a series of *n* layers with internal heat generation and stable tectonic environment can be determined by the equation

$$Q\_s = Q\_b + A\_1 H\_1 + A\_2 H\_2 + \dots \ + A\_n H\_{n'} \tag{1}$$

where *Qb* is the heat flow from the base of the *n* layers being considered, typically the upper mantle, *A*1 is the heat generation in layer 1, *H*1 is the thickness of layer 1, *A*2 is the heat generation in layer 2, *H*2 is the thickness of layer 2, *An* is the heat generation in layer *n*, and *Hn* is the thickness of layer *n*.

Temperature within a surface layer of thickness *L*, and uniform heat generation of *A*0, is shown by Carslaw and Jaeger [24], described by the equation:

$$T = \left[\frac{A\_0 z}{2K}\right] \left(2L - z\right) + \frac{Q\_b z}{K} + T\_{s'} \tag{2}$$

where *T* is the temperature, *z* is the depth of observation, *K* is the thermal conductivity, *Ts* is the mean annual surface temperature, and *Qb* is the heat flow at the base of the layer. The geometry applicable to this equation is shown in Figure 3, and this equation can be applied to estimate the temperature versus depth in a granite batholith with uniform radiogenic heat generation.

Potassium usually has a different genesis to U and Th. It is a very commonly occurring element in the earth's crust. The abundance (concentration) of K in most igneous rocks is commonly observed to be about 2–6%, and it also occurs abundantly in many sedimentary rocks. Beardsmore and Cull [7] have discussed the specific abundance of the radioactive isotope (K-40) within normal potassium distributions in rocks. The surface K "merge" distribution map [19] can be misleading, as the K concentrations tend to follow recent sedimentary depositional features, such as seasonally dry riverbeds. In such deposits, the time of year of that an airborne radiometric survey was acquired has implications for interpretation of the observed K concentration [20–22]. Further, as will be seen in the following sections, the contribution of K to radiogenic heat generation is usually less than the contribution of U and

In this review, surface-based geophysical or laboratory measurements of U, Th, and K are emphasized rather than airborne measurements. This is because surface-based (on-theground) measurements of outcropping rocks tend to be more directly indicative of the particular rock type being investigated, rather than a "bulk" estimate of the chemistry of the

This section reviews the basic, but well-known, concepts to understand heat generation, heat

Heat flow at the surface on the earth (Qs) over a series of *n* layers with internal heat generation

where *Qb* is the heat flow from the base of the *n* layers being considered, typically the upper mantle, *A*1 is the heat generation in layer 1, *H*1 is the thickness of layer 1, *A*2 is the heat generation in layer 2, *H*2 is the thickness of layer 2, *An* is the heat generation in layer *n*, and *Hn* is the

Temperature within a surface layer of thickness *L*, and uniform heat generation of *A*0, is shown

*b*

*s*

(2)

( ) <sup>0</sup> 2 , <sup>2</sup>

*A z Q z T Lz T K K*

= -+ + ê ú

é ù

ë û

11 22 ... , *Q Q AH AH AH s b n n* =+ + + + (1)

rocks being sampled by an airborne survey. This will be discussed further below.

flow, and simple temperature distributions with the earth's crust.

and stable tectonic environment can be determined by the equation

by Carslaw and Jaeger [24], described by the equation:

**2.1. Heat flow and temperature due to heat generation in a radiogenic layer**

Th [7,16,23].

54 Advances in Geothermal Energy

**2. Theory**

thickness of layer *n*.

**Figure 3.** theory for temperature model. Note that the same modeling geometry is used for relatively thin enrichment layer models in section 7. In the latter case, the radiogenic granite layer is a layer where U- and Th-rich fluids have been rapidly emplaced due to chemical alteration processes, largely due to weathering and erosion.

The more general time-variable equation describing the gradual rise in temperature due to thermal adjustment of the emplacement of a layer of thickness *L* at an initial time of *t*=0 is given by Carslaw and Jaeger [24] for 0 < *z* < *L*:

$$T = \left[\frac{\alpha A\_0 t}{K}\right] \left[1 - 4i^2 \text{erfc}\left(\frac{z}{\tau}\right) + 2i^2 \text{erfc}\left(\frac{L+z}{\tau}\right) - 2i^2 \text{erfc}\left(\frac{L-z}{\tau}\right)\right] + T\_0 \text{erf}\left(\frac{z}{\tau}\right) + \frac{Q\_b z}{K} \tag{3}$$

where *T* is the temperature, *z* is the depth, *t* is the time, *α* is the thermal diffusivity, *K* is the thermal conductivity, *A*<sup>0</sup> is the heat generation, *T*<sup>0</sup> is the initial temperature of layer at *t* = 0, *Qb* is the heat flow at the base of layer, erf(*x*) is the error function, erfc(*x*) is the complimentary error function [1 − erf(*x*)], *i* 2 erfc(*x*) is the second integration of the error function, and *τ* = (4*αt*) ½. The temperature modeling below is largely based on these three equations. A theoretical case, based on Eq. (3), where a radiogenically enriched layer of variable thickness occurs, by chemical alteration processes, at the earth's surface is also investigated in Section 7.

The case of sedimentary layers overlying a radiogenic granite is also important. This situation has the potential to cause a "thermal blanketing effect," where the overlying sediments possess significantly lower thermal conductivity than the underlying rocks. In many cases, this has been shown to create anomalously high temperatures within the upper sedimentary sequence. Such cases have been treated theoretically by Carslaw and Jaeger [24]. In the Western Aus‐ tralian perspective, the "thermal blanketing" effect was investigated in the case of a geothermal (temperature) anomaly in the northern Perth Basin (Gillingarra anomaly) and is discussed below. This investigation showed that the wavelength of the observed heat flow anomaly was too low (and magnitude too high) to accommodate a radiogenic heat source from an estimated 10-km basement depth.

## **2.2. Heat generation from U, Th, and K**

Studies of radiogenic heat generation from the decay of unstable isotopes of U, Th, and K have been summarized by Kappelmeyer and Haenel [23], Jessop [25], Beardsmore and Cull [7], and Middleton [16]. The generalized relation [16] used for the present study is

$$A\_0 = \begin{pmatrix} 0.26 \end{pmatrix} \text{(U ppm)} + \begin{pmatrix} 0.07 \end{pmatrix} \text{(Th ppm)} + \begin{pmatrix} 0.1 \end{pmatrix} \text{(K } \% \text{)}.\tag{4}$$

where *A*<sup>0</sup> is the heat generation, (U ppm) is the concentration of uranium in ppm, (Th ppm) is the concentration of thorium in ppm, and (K %) is the concentration of potassium in percent. The relationship used in Eq. (4) is consistent with that of Beardsmore and Cull [7]. It is recognized that by using these constants in Eq. (4), the density of all rocks is assumed to be a constant value of 2700 kg m−3 (as is similarly assumed in the example of Beardsmore and Cull [7]). Essentially, this assumption represents a generalized conversion formula.

### **2.3. Field measurements with an RS-125 gamma-ray spectrometer**

For studies reviewed in this chapter, many of the sites measured for U, Th, and K assay were obtained using an RS-125 gamma-ray spectrometer. The RS-125 instrument is calibrated, via "assay mode", to provide an approximation of U, Th, and K concentrations in a half-space being sampled. The manufacturer indicates that an instrumental error of about 10% is incurred in the RS-125 in the assay mode, due to the programmed calibration matrix. The volume sampled is assumed to be a perfect half-space, and variations from this perfect half-space can also cause some errors in the assay determination. However, observation locations were chosen so as to minimize this type of error. The manufacturer's specifications of the RS-125 spectrometer indicate that the instrument reflects a sample area of a radius of approximately 1 m in the half-space (Figure 3).

Calibration procedures were carried out prior to, and during, the field acquisition with the RS-125. An important factor for calibration is the sample time of each measurement. In determining this, the manufacturer asserts that the calibration matrix, which is used to optimize the spectral definition of the three isotopes, is designed to provide the minimum error (< 5%) for sample times of approximately 4 minutes. However, this has to be balanced against the requirement for long sample times (often over many hours) for rocks with low U, Th, and K concentrations, as expected in many of the granites being investigated. It was finally decided to base all measurements for studies reported in this review to be 5 minutes [16]. Further, the RS-125 field assays were compared with U, Th, and K assays obtained by laboratory meas‐ urements for sample sites in the Leeuwin Complex in the southwest of Western Australia, and good correspondence was observed (Figure 4).

The case of sedimentary layers overlying a radiogenic granite is also important. This situation has the potential to cause a "thermal blanketing effect," where the overlying sediments possess significantly lower thermal conductivity than the underlying rocks. In many cases, this has been shown to create anomalously high temperatures within the upper sedimentary sequence. Such cases have been treated theoretically by Carslaw and Jaeger [24]. In the Western Aus‐ tralian perspective, the "thermal blanketing" effect was investigated in the case of a geothermal (temperature) anomaly in the northern Perth Basin (Gillingarra anomaly) and is discussed below. This investigation showed that the wavelength of the observed heat flow anomaly was too low (and magnitude too high) to accommodate a radiogenic heat source from an estimated

Studies of radiogenic heat generation from the decay of unstable isotopes of U, Th, and K have been summarized by Kappelmeyer and Haenel [23], Jessop [25], Beardsmore and Cull [7], and

where *A*<sup>0</sup> is the heat generation, (U ppm) is the concentration of uranium in ppm, (Th ppm) is the concentration of thorium in ppm, and (K %) is the concentration of potassium in percent. The relationship used in Eq. (4) is consistent with that of Beardsmore and Cull [7]. It is recognized that by using these constants in Eq. (4), the density of all rocks is assumed to be a constant value of 2700 kg m−3 (as is similarly assumed in the example of Beardsmore and Cull

For studies reviewed in this chapter, many of the sites measured for U, Th, and K assay were obtained using an RS-125 gamma-ray spectrometer. The RS-125 instrument is calibrated, via "assay mode", to provide an approximation of U, Th, and K concentrations in a half-space being sampled. The manufacturer indicates that an instrumental error of about 10% is incurred in the RS-125 in the assay mode, due to the programmed calibration matrix. The volume sampled is assumed to be a perfect half-space, and variations from this perfect half-space can also cause some errors in the assay determination. However, observation locations were chosen so as to minimize this type of error. The manufacturer's specifications of the RS-125 spectrometer indicate that the instrument reflects a sample area of a radius of approximately

Calibration procedures were carried out prior to, and during, the field acquisition with the RS-125. An important factor for calibration is the sample time of each measurement. In determining this, the manufacturer asserts that the calibration matrix, which is used to optimize the spectral definition of the three isotopes, is designed to provide the minimum error (< 5%) for sample times of approximately 4 minutes. However, this has to be balanced against

*A*<sup>0</sup> = 0.26 U ppm + 0.07 Th ppm + 0.1 K % , ( ) ( ) ( ) ( ) ( ) ( ) (4)

Middleton [16]. The generalized relation [16] used for the present study is

[7]). Essentially, this assumption represents a generalized conversion formula.

**2.3. Field measurements with an RS-125 gamma-ray spectrometer**

10-km basement depth.

56 Advances in Geothermal Energy

**2.2. Heat generation from U, Th, and K**

1 m in the half-space (Figure 3).

**Figure 4.** A diagram showing the RS-125 sampling zone, which is essentially a 1-m-radius hemisphere below the earth's surface.

**Figure 5.** Correlation of RS-125 to laboratory assays for Leeuwin Complex samples.

#### **2.4. Heat generation derived from gamma logs**

Ryback [26] and Bücker and Ryback [27] proposed empirical equations to determine heat generation (*A*0, µWm−3) from wire-line gamma-ray (GR) logs run in (petroleum and sometimes mineral) exploration wells or boreholes. The Bücker and Ryback [27] correlation appears to be superior to the earlier correlation by Ryback [26], and is expressed as

$$A\_0 = 0.0158 \text{ (GR)} \text{ - 0.0126 }, \tag{5}$$

where *A*<sup>0</sup> is the heat generation (in units of µWm−3) and GR is the logged gamma-ray value in API (American Petroleum Institute) units. A gamma-ray API is defined, and further discussed in relation to other radioactivity measurement units, in Society of Petroleum Engineers (SPE) [28]. When looking at gamma-ray measurements greater than 20 API units, this equation simplifies to *A*0 ≈ 0.0158 (GR). For this study, it is assumed that *A*0 = 0.0158 (GR).

Beardsmore and Cull [7] investigated this relationship with regard to the East Yeeda 1 well from the Canning Basin in Western Australia. They derived a relationship between heat generation and the gamma-ray and density logs from this well, which had total count and spectral gamma-ray logs. These workers derived a relationship between the heat generation and the two wire-line-log values to be

$$A\_0 \approx 0.005 \text{ (GR) (RHOB)},\tag{6}$$

where RHOB (g cm−3) is the bulk density and GR (API units) are the readings from the wireline log. The authors suggest that an uncertainty in *A*0 is about ± 0.5 µW m−3. These workers place the caveat that this relationship may be local correlation, rather than a general relation‐ ship. If we assume a general rock density of 2700 kg m−3, then the Beardsmore and Cull relationship becomes

$$A\_0 \approx 0.014 \text{ (GR)}.\tag{7}$$

On the basis of these considerations, Eqs (5) and (7), heat generation may be in a very general sense correlated to the wire-line log gamma-ray value by a relationship represented by the range 0.014 (GR) < *A*0 < 0.016 (GR), with GR in API units and *A*0 in µW m−3.

From 2011 to 2015, a series of studies was carried out with a commercial gamma-ray spectral scintillation detector, RS-125 instrument, as described below. This instrument was used to measure many of the heat generation estimates reported in this study [16,18]. It is of importance to establish how this instrument performs, and provides results, consistent with other radioactive measuring instruments. Since 2011, the results shown in Figure 6 have been accumulated, which can provide an approximation of dose rate (DR; as reported by the instrument) with respect to uranium (ppm), thorium (ppm), and potassium (K). This correla‐ tion can provide a correlation of dose rate to heat generation via Eq. 4. Correlations, according to the data in Figure 6, can be established between DR and *A*0, with heat generation values (*A*0) being determined according to Eq. (4). Regression curves have been established for (1) DR total and (2) DR <300 µR/hr. It is believed that DR (µR/hr) < 300 may be more appropriate for correlation studies, because of the nonlinearity of the DR versus *A*<sup>0</sup> relationship above the value of 300 µR/hr. In the case of all the data being incorporated, a second-degree polynomial curve DR = 0.6292 (*A*0) 2 + 47.186 (*A*0), with regression coefficient of *R*<sup>2</sup> = 0.9254, is found to provide the best match (Figure 6).

**2.4. Heat generation derived from gamma logs**

58 Advances in Geothermal Energy

and the two wire-line-log values to be

relationship becomes

superior to the earlier correlation by Ryback [26], and is expressed as

Ryback [26] and Bücker and Ryback [27] proposed empirical equations to determine heat generation (*A*0, µWm−3) from wire-line gamma-ray (GR) logs run in (petroleum and sometimes mineral) exploration wells or boreholes. The Bücker and Ryback [27] correlation appears to be

where *A*<sup>0</sup> is the heat generation (in units of µWm−3) and GR is the logged gamma-ray value in API (American Petroleum Institute) units. A gamma-ray API is defined, and further discussed in relation to other radioactivity measurement units, in Society of Petroleum Engineers (SPE) [28]. When looking at gamma-ray measurements greater than 20 API units, this equation

Beardsmore and Cull [7] investigated this relationship with regard to the East Yeeda 1 well from the Canning Basin in Western Australia. They derived a relationship between heat generation and the gamma-ray and density logs from this well, which had total count and spectral gamma-ray logs. These workers derived a relationship between the heat generation

where RHOB (g cm−3) is the bulk density and GR (API units) are the readings from the wireline log. The authors suggest that an uncertainty in *A*0 is about ± 0.5 µW m−3. These workers place the caveat that this relationship may be local correlation, rather than a general relation‐ ship. If we assume a general rock density of 2700 kg m−3, then the Beardsmore and Cull

On the basis of these considerations, Eqs (5) and (7), heat generation may be in a very general sense correlated to the wire-line log gamma-ray value by a relationship represented by the

From 2011 to 2015, a series of studies was carried out with a commercial gamma-ray spectral scintillation detector, RS-125 instrument, as described below. This instrument was used to measure many of the heat generation estimates reported in this study [16,18]. It is of importance to establish how this instrument performs, and provides results, consistent with other radioactive measuring instruments. Since 2011, the results shown in Figure 6 have been accumulated, which can provide an approximation of dose rate (DR; as reported by the instrument) with respect to uranium (ppm), thorium (ppm), and potassium (K). This correla‐

range 0.014 (GR) < *A*0 < 0.016 (GR), with GR in API units and *A*0 in µW m−3.

simplifies to *A*0 ≈ 0.0158 (GR). For this study, it is assumed that *A*0 = 0.0158 (GR).

*A*<sup>0</sup> = 0.0158 GR - 0.0126 , ( ) (5)

*A*<sup>0</sup> » 0.005 GR RHOB , ( ) ( ) (6)

*A*<sup>0</sup> » 0.014 GR . ( ) (7)

**Figure 6.** Correlation between RS-125 dose rate in (üR/hr) and heat generation (µW/m3 ) from the data derived by field measurements [16,17].

In the special case of DR being < 300 µR/hr, a regression match is found to be DR = 43.21 (*A*0), with a regression coefficient of *R*<sup>2</sup> = 0.8866. In this case, an intercept at zero is also assumed. The observed data for *A*0 less than about 1.3 µWm−3 (DR < about 15 µR/hr) appear to be showing a limitation in the calibration of the scintillation crystal to accurately measure the uranium, thorium, and potassium elemental concentration spectra, and more research is required to look into why this is happening. Despite the concerns of the calibration of the RS-125 crystal at low U, Th, and K concentrations, the author believes, that sufficient confidence can be placed in the data for DR >15 µR/hr, that an acceptable correlation between DR and *A*<sup>0</sup> can be sustained to be DR = 0.43.21 (*A*0). Assuming the SPE correlation between DR and API gamma-ray measurements is DR = 0.67 API units [28], the data found by Middleton [16] indicate a correlation of DR ≈ 0.0155 (GR). This correlation is very close to the Bücker and Ryback [27] approximation.

It is proposed that an approximate conversion factor for GR (API units) to *A*<sup>0</sup> (µWm−3 units) for this investigation to be

$$A\_0 = \text{0.015 (GR)}\,\text{.}\tag{8}$$

This equation is the mean of the Bücker and Ryback [27] and Beardsmore and Cull [7] approximations. Some estimates of heat generation, derived from this equation, of the gammaray response of basement rocks observed in wire-line logs from the Perth Basins are shown below.

## **3. Perth basin**

The Perth Basin lies in the southwest corner of Western Australia (Figure 1). It contains sediments of Phanerozoic age, which are surrounded by crystalline granitoid and metamor‐ phic rocks of the Yilgarn Craton, Leeuwin Complex, and Northampton Complex (Figure 1). The Perth Basin will be treated in three parts: the southern Perth Basin [29], the central Perth Basin [30–32], and the northern Perth Basin [33].

## **3.1. Southern perth basin**

Heat generation in the granites of the Leeuwin Block, which immediately crops out to the west of the Vasse Shelf, may have relevance for the exploration for geothermal energy "hot spots" in the Vasse Shelf sediments (Figure 1). The significance of measuring heat generation in the Leeuwin Complex is that it is expected that these rocks underlie the elements of the Perth Basin immediately to the east of these outcropping sediments. The elements of the Perth Basin in question are the Vasse Shelf sediments (largely of Permian age) and the deeper and thicker Bunbury Trough sediments (largely of Cretaceous to Permian age) [29]. It is uncertain if the Leeuwin Complex crystalline rocks extend under the Bunbury Trough, which has been described geologically by Crostella and Backhouse [29].

Western Australia has several regions where hot granitoids are known to occur in outcrop. The locations of high radioactivity rocks, including hot granitoids, can initially be detected by airborne radiometric surveys. These data have been acquired for Western Australia through a number of cooperative projects between the Geological Survey of Western Australia and Geoscience Australia, and maps showing these data are publically available. The Vasse region was flown in 2011. The U concentration map for the northern part of the Vasse Region, derived from this survey, is shown in Figure 2a, and the Th map for the same region is shown in Figure 2b. A ground-truthing survey was carried out by the members of the Petroleum Division of the Department of Mines and Petroleum in May 2013. This survey used a RS-125 gamma-ray spectrometer to obtain surface values of U, Th, and K to compare with the airborne data. Middleton [16,34] carried out a similar study to investigate geothermal heat generation for the Darling Range granitoids. Surface sampling of granitoids of the Leeuwin Complex was previously carried out by the Geological Survey of Western Australia and Geoscience Aus‐ tralia, and these sites were also revisited during 2011 and 2012 to compare the spectrometer values to those obtained by the previous laboratory-based measurements of U, Th, and K (main locations shown in Figure 8). The U, Th, and K concentrations in the granitoids can be converted into approximate radiogenic heat generation using Eq. (4); see Figure 7.

It is proposed that an approximate conversion factor for GR (API units) to *A*<sup>0</sup> (µWm−3 units)

This equation is the mean of the Bücker and Ryback [27] and Beardsmore and Cull [7] approximations. Some estimates of heat generation, derived from this equation, of the gammaray response of basement rocks observed in wire-line logs from the Perth Basins are shown

The Perth Basin lies in the southwest corner of Western Australia (Figure 1). It contains sediments of Phanerozoic age, which are surrounded by crystalline granitoid and metamor‐ phic rocks of the Yilgarn Craton, Leeuwin Complex, and Northampton Complex (Figure 1). The Perth Basin will be treated in three parts: the southern Perth Basin [29], the central Perth

Heat generation in the granites of the Leeuwin Block, which immediately crops out to the west of the Vasse Shelf, may have relevance for the exploration for geothermal energy "hot spots" in the Vasse Shelf sediments (Figure 1). The significance of measuring heat generation in the Leeuwin Complex is that it is expected that these rocks underlie the elements of the Perth Basin immediately to the east of these outcropping sediments. The elements of the Perth Basin in question are the Vasse Shelf sediments (largely of Permian age) and the deeper and thicker Bunbury Trough sediments (largely of Cretaceous to Permian age) [29]. It is uncertain if the Leeuwin Complex crystalline rocks extend under the Bunbury Trough, which has been

Western Australia has several regions where hot granitoids are known to occur in outcrop. The locations of high radioactivity rocks, including hot granitoids, can initially be detected by airborne radiometric surveys. These data have been acquired for Western Australia through a number of cooperative projects between the Geological Survey of Western Australia and Geoscience Australia, and maps showing these data are publically available. The Vasse region was flown in 2011. The U concentration map for the northern part of the Vasse Region, derived from this survey, is shown in Figure 2a, and the Th map for the same region is shown in Figure 2b. A ground-truthing survey was carried out by the members of the Petroleum Division of the Department of Mines and Petroleum in May 2013. This survey used a RS-125 gamma-ray spectrometer to obtain surface values of U, Th, and K to compare with the airborne data. Middleton [16,34] carried out a similar study to investigate geothermal heat generation for the Darling Range granitoids. Surface sampling of granitoids of the Leeuwin Complex was

*A*<sup>0</sup> = 0.015 GR . ( ) (8)

for this investigation to be

60 Advances in Geothermal Energy

below.

**3. Perth basin**

**3.1. Southern perth basin**

Basin [30–32], and the northern Perth Basin [33].

described geologically by Crostella and Backhouse [29].

**Figure 7.** Correlation of U or Th concentration to heat generation, using Eq. (4).

Previous works [4,5,16,34,35] have been carried out on heat generation in granitoids of the western part of the Yilgarn Craton. The heat generation of outcropping hot granites in the Yilgarn Craton has been measured in the range of 1–10 µWm−3 and is similar to published values for the Cooper Basin Hot Dry Rock (HDR) granites (ca. 10 µWm−3) [6]. Some previous heat generation values were published for the Leeuwin Complex by Middleton and Stevens [17], although U, Th, and K analyses have been measured by laboratory techniques from outcrop samples for selected locations (see Figure 7). More recent heat generation values have been acquired in the northern part of the Leeuwin Complex and are herein reported.

Table 1 shows the mean U, Th, and K concentrations (measured with at least three repeats) at the 16 localities in this Leeuwin and south Perth Basin reviewed in this study. The elemental abundances observed at these sites were converted to heat generation in units of µWm−3 using the factors in Eq. (4). Table 1 shows the relationship between U, Th, K and heat generation (designated as *A*0). The table indicates that the observed mean heat generation at the locations investigated ranged between 3.8 µWm−3 and 40.1 µWm−3 throughout the Leeuwin Complex region. The heat generation data exhibit a good linear relationship between U and Th concen‐ trations. The very high U and Th concentrations, and calculated heat generation of 40.1 µWm −3, at the Meelup Park anomaly location (MP in Figure 8) near Dunsborough are indeed anomalous, which was investigated because of the high airborne radiometric signals (see Figures 8–10). The reason may be due to the concentration of these elements in a thin laterite or gravel layer residing on the top of a hot underlying granite body. No obvious outcrop of the underlying granitic rocks was observed during the survey at this locality.

**Figure 8.** South Perth Basin structural and other tectonic elements.

The airborne signature suggests that the underlying hot granite may have an areal extent of about 5 km2 . Figure 11 shows the thorium anomalies superimposed on the aeromagnetic data at the Meelup Park location. The anomalous thorium appears to be located in a low magnetic region in the range of −50 to −200 nT (green to light blue in Figure11). The low magnitude of the magnetic signature of the thorium-rich regions implies that the thorium may be located in low magnetic crystalline rocks. This tends to suggest that the thorium (and uranium) may be associated with the more felsic rocks in the Leeuwin Complex, which is also supported by evidence from the Sugarloaf Rock locality. The Sugarloaf Rock locality is located on the coast about 3 km south of Cape Naturaliste (see Figure 8). Laboratory measurements by the Geological Survey of Western Australia indicate at least two igneous rock suites at the locality: (a) with U ~3.9 ppm and Th ~92.6 ppm and (b) U ~1.7 ppm and Th ~25.4 ppm. Surface measurements with the spectrometer yielded that (a) U ranged between 3–8.5 ppm and Th ~108 ppm for an observed felsic suite of rocks and (b) U ~1.3 ppm and Th ~21 ppm for a mafic suite. It should be noted that the felsic suite of igneous rocks are expected to have a low magnetic signature and the mafic suite a high magnetic signature. Measurements at this site support the observations at Meelup Park ("MP" in Figure 8) that the low-magnetic rocks seem to possess a higher concentration of radioactive elements than do the high-magnetic rocks.

A simple 1D model [24] is used to determine the temperature at depth in sediments, which are underlain by a granitic body with uniform heat generation and limited depth extent. The main unknown parameters in this modeling exercise are heat generation within the layer of hot granite (*A*0) and its thickness (*L*). The other parameters, such as thermal conductivity (*K*), basal heat flow (*Qb*), and surface temperature (*Ts*), are relatively well known [4,35,36], and commonly observed values have been assumed (see Table 2). Several cases with different heat generation and depth extents of hot granites are modeled to investigate temperatures at depth within the western part of the Vasse Shelf. Granitoid bodies with depth extents of 2000, 3500, and 6000 m are considered for this modeling. Heat generation is assumed to vary between 4 and 20 µWm−3. Figure 12 shows temperature versus depth for each of the four cases with correspond‐ ing parameters shown in Table 2. The results suggest that temperatures at 3000-m depth in the Vasse Shelf may have temperatures as high as 180 °C, if underlain by a very radioactive granitoid *A*0 (> 20 µWm−3) of sufficient thickness (> 6000 m).


**Table 2.** Parameters used in south Perth Basin temperature modeling.

investigated ranged between 3.8 µWm−3 and 40.1 µWm−3 throughout the Leeuwin Complex region. The heat generation data exhibit a good linear relationship between U and Th concen‐ trations. The very high U and Th concentrations, and calculated heat generation of 40.1 µWm −3, at the Meelup Park anomaly location (MP in Figure 8) near Dunsborough are indeed anomalous, which was investigated because of the high airborne radiometric signals (see Figures 8–10). The reason may be due to the concentration of these elements in a thin laterite or gravel layer residing on the top of a hot underlying granite body. No obvious outcrop of

The airborne signature suggests that the underlying hot granite may have an areal extent of

at the Meelup Park location. The anomalous thorium appears to be located in a low magnetic region in the range of −50 to −200 nT (green to light blue in Figure11). The low magnitude of the magnetic signature of the thorium-rich regions implies that the thorium may be located in low magnetic crystalline rocks. This tends to suggest that the thorium (and uranium) may be associated with the more felsic rocks in the Leeuwin Complex, which is also supported by evidence from the Sugarloaf Rock locality. The Sugarloaf Rock locality is located on the coast about 3 km south of Cape Naturaliste (see Figure 8). Laboratory measurements by the Geological Survey of Western Australia indicate at least two igneous rock suites at the locality: (a) with U ~3.9 ppm and Th ~92.6 ppm and (b) U ~1.7 ppm and Th ~25.4 ppm. Surface

. Figure 11 shows the thorium anomalies superimposed on the aeromagnetic data

the underlying granitic rocks was observed during the survey at this locality.

**Figure 8.** South Perth Basin structural and other tectonic elements.

about 5 km2

62 Advances in Geothermal Energy

It is worth noting that the hydrogeological study by Wharton [37] reported temperatures within the sediments in the range of 38–40 °C at approximately 1000-m depth. The modeled temperatures in Figure 12 are consistent with the temperatures reported by Wharton [37] for cases 2, 3, 4, and 7 in Table 2. For example, the average case of a radiogenic-rich granitoid with a thickness of 3500 m and heat generation of 8 µWm−3 can satisfactorily yield the temperatures observed in the study by Wharton [37], as can a granitoid of thickness of 2000 m and heat generation of 20 µWm−3. Interestingly, in the former case, the temperature extrapolated to 5000 m will be 111°C and 133 °C in the latter case. However, it is unlikely that sediment thicknesses of 5000 m exist in the part of the Vasse Shelf.

### **3.2. Central Perth basin**

#### *3.2.1. Darling range*

The Darling Range is a physiographic feature that typifies part of the western margin of the Yilgarn Craton in the southwest of Western Australia. It is dominated by the Darling Fault (Figure 1), which separates the Yilgarn Craton from the Perth Basin [30,32,38,39], and is one of the longest fault zones in the world, extending for about 1000 km. In the vicinity of Perth, the Darling Fault is expressed as a surface scarp of about 300 m, but the subsurface expression of the fault zone on seismic reflection data [40] is observed to extend to over 20 seconds twoway time (approximately 50 km). Most importantly, it is not known if the same rocks that comprise the Yilgarn Craton underlie the Perth Basin, as it is interpreted as a zone of major continental collision and orogenesis. Commonly within the Craton, there is an upper laterite (Regolith) profile that has been formed on the Darling Range granitoids and associated igneous lithologies. In many places, the upper laterite profile has been removed by erosion, generally in incised valleys, and in those places relatively fresh granite is exposed [40,41]. Tertiary channels also occur within the regolith in the Darling Range, and these have also been investigated for surface U, Th, and K (Table 1). Field measurements of abundance of naturally occurring radiogenic isotopes were made in order to estimate radiogenic heat generation of granitoids in the Darling Range region of the Yilgarn Craton. Fifteen locations are reported from this region (Table 1). The field measurements on outcropping granitoids, laterite, and tertiary channel sediments indicate a range of heat generation between 3.3 and 10.2 µW m−3.

An alternative approach to evaluating heat generation is to use Eq. (4) to estimate U, Th, and K concentrations from calibrated airborne radiometric data, which was investigated by Middleton [16]. This is now reviewed. A heat generation map (Figure 6), based on airborne radiometric data, was originally published by Middleton [34]. This type of map is based on U, Th, and K airborne-radiometric maps released by Geoscience Australia [10,22]. Essentially, these data represent an approximation of the surface heat generation for all the rocks in each map pixel, which has dimensions of approximately 1.3 × 1.3 km of variable rock outcrop. However, such an average may not be representative of any particular outcrop that may fall in the approximate pixel area of 1.3 × 1.3 km. Accordingly, while such maps are a good general guide to high heat generation regions, the actual heat generation values may be several times less than specific geological units. As shown in Figure 6, the Darling Range has a surface radiometric heat generation in the range of 2.5–5.0 µW m−3, which is lower than the field measurement of all the outcropping granites in the region. This, however, might be expected, as the data in Figure 6 represent an average of outcropping rocks in the pixel area of 1.3 × 1.3 km. Nevertheless, the map may provide a good general guide to regions of high surface radiogenic heat generation.

Very little is known about spatial variation of heat generation (and indeed U and Th concen‐ tration) in radiogenic granitoids in the Darling Range of Western Australia. Table 2 of Middleton [16], and included in Table 1 of this work, cites the variation between four sites at the Glen Forrest (J Forrest Pk 1–4 in Table 1) location. These four sites were within a total distance of about 100 m, and are within and around an abandoned granite quarry. The purpose of these closely spaced measurements was to investigate the spatial variability of heat generation over small distances within a batholith. The table indicates that local heat genera‐ tion within this batholith can vary between 6 and 8 µW m−3 over the 100 m. These observations suggest that heat generation within the batholith is relatively constant over such distances. The study of Jaeger [4] also found fairly constant heat generation values (mean of 8.9 µW m−3 with a standard error of 0.6 µW m−3) over a borehole depth interval of 300 m within a radiogenic granite near the town of Doodlakine, located also in the Yilgarn Craton (see Figure 1) Jaeger's measurements lie about 1 km to the east of where the present measurements were made (see Table 1, 8.2 µW m−3). Also, a study of several surface sites within a hot granitoid body in the Dunsborough region of the Southern Perth Basin (Table 1) exhibits a range of heat generation (12.3–15.7 µW m−3, based on U, Th, and K concentrations) within an approximately 500-m stretch of land, which is assumed to comprise the same granitoid body. Spatial variability of U, Th, and *A*0 within single granitoid bodies is a subject that requires considerably further work.

**3.2. Central Perth basin**

64 Advances in Geothermal Energy

radiogenic heat generation.

The Darling Range is a physiographic feature that typifies part of the western margin of the Yilgarn Craton in the southwest of Western Australia. It is dominated by the Darling Fault (Figure 1), which separates the Yilgarn Craton from the Perth Basin [30,32,38,39], and is one of the longest fault zones in the world, extending for about 1000 km. In the vicinity of Perth, the Darling Fault is expressed as a surface scarp of about 300 m, but the subsurface expression of the fault zone on seismic reflection data [40] is observed to extend to over 20 seconds twoway time (approximately 50 km). Most importantly, it is not known if the same rocks that comprise the Yilgarn Craton underlie the Perth Basin, as it is interpreted as a zone of major continental collision and orogenesis. Commonly within the Craton, there is an upper laterite (Regolith) profile that has been formed on the Darling Range granitoids and associated igneous lithologies. In many places, the upper laterite profile has been removed by erosion, generally in incised valleys, and in those places relatively fresh granite is exposed [40,41]. Tertiary channels also occur within the regolith in the Darling Range, and these have also been investigated for surface U, Th, and K (Table 1). Field measurements of abundance of naturally occurring radiogenic isotopes were made in order to estimate radiogenic heat generation of granitoids in the Darling Range region of the Yilgarn Craton. Fifteen locations are reported from this region (Table 1). The field measurements on outcropping granitoids, laterite, and tertiary channel sediments indicate a range of heat generation between 3.3 and 10.2 µW m−3.

An alternative approach to evaluating heat generation is to use Eq. (4) to estimate U, Th, and K concentrations from calibrated airborne radiometric data, which was investigated by Middleton [16]. This is now reviewed. A heat generation map (Figure 6), based on airborne radiometric data, was originally published by Middleton [34]. This type of map is based on U, Th, and K airborne-radiometric maps released by Geoscience Australia [10,22]. Essentially, these data represent an approximation of the surface heat generation for all the rocks in each map pixel, which has dimensions of approximately 1.3 × 1.3 km of variable rock outcrop. However, such an average may not be representative of any particular outcrop that may fall in the approximate pixel area of 1.3 × 1.3 km. Accordingly, while such maps are a good general guide to high heat generation regions, the actual heat generation values may be several times less than specific geological units. As shown in Figure 6, the Darling Range has a surface radiometric heat generation in the range of 2.5–5.0 µW m−3, which is lower than the field measurement of all the outcropping granites in the region. This, however, might be expected, as the data in Figure 6 represent an average of outcropping rocks in the pixel area of 1.3 × 1.3 km. Nevertheless, the map may provide a good general guide to regions of high surface

Very little is known about spatial variation of heat generation (and indeed U and Th concen‐ tration) in radiogenic granitoids in the Darling Range of Western Australia. Table 2 of Middleton [16], and included in Table 1 of this work, cites the variation between four sites at the Glen Forrest (J Forrest Pk 1–4 in Table 1) location. These four sites were within a total distance of about 100 m, and are within and around an abandoned granite quarry. The purpose

*3.2.1. Darling range*

Heat flow at the base of granite batholiths (*Qb*) in the Yilgarn Craton and Perth Basin is poorly known. Conclusions that can be drawn from the work of Jaeger (1970) are that (1) heat flow at the base of his presumed Yilgarn upper-crustal heat-generation layer is about 0.63 heat flow units (hfu, which convert to 26.4 mW m−2) and (2) the radiogenic upper-crustal layer is about 4.5 km thick. Later works [34,35] in the central Perth Basin indicate that it commonly possesses a heat flow in the vicinity of 86 mW m−2. However, in a locality at the Gillingarra 6 water bore (G6 in Figure 16), surface heat flow was found to be approximately 130 mW m−2 [34], which is anomalously high for the central Perth Basin. One could conclude that this thermal anomaly is due to anomalously elevated heat generation in the granitoids underlying the G 6 locality. On the basis of the results from the anomalous G 6 data, a range for the radiogenic heat generation from granitoids under the Perth Basin may be proposed from a simplification of Eq. (1), which is *Qs* = *Qb* + *A*1*H*1. Thus, in the non-anomalous situation, if *Qs* = 86 mW m−2, *Qb* = 26 mW m−2, and *H*1 = 5 km, then the radiogenic heat production in the sub-basin granitoid layer is estimated to be 17.2 µW m−3. However, in the anomalous case at G 6, if *Qs* = 130 mW m−2, *Qb* = 26 mW m−2, and *H*1 = 5 km, then the radiogenic heat production in the sub-basin granitoid layer is estimated to be 26 µW m−3. Both of these values (17.2 and 26 µW m−3) for heat generation appear to be unreasonably high. Perhaps, as suggested by Middleton [34,41], the G 6 temper‐ ature anomaly may be caused by either anomalous thermal conductivity or convective groundwater flow. This currently remains unresolved. These observations are consistent with other regional studies [41].

For temperature modeling, it is important to establish mean annual surface temperature data that are necessary to estimate geothermal gradients [7,16,35]. These data can be derived from Australian Bureau of Metrology data covering the past five or six decades in the Perth Basin and other Western Australian basins. The difficulty is to how to process these data to arrive at a meaningful mean surface temperature that is relevant to geothermal processes. This is beyond the scope of this paper. For the purposes of modeling in the current study, a value for *Ts* of 21 °C is assumed, which is consistent with values proposed by Jaeger [4].

Temperature modeling for the Darling Range granitoids [16], and Table 3, was carried out using the geometry shown in Figure 3 and based upon Eq. (3). The purpose of the modeling is to determine if temperatures in the metropolitan Darling Range are high enough to support geothermal applications. In this modeling, the assumed parameters are the thicknesses of the radiogenic granite batholith (*L*) varying between 3 and 6 km. The complete temperature profile from surface (*z*=0) to a maximum depth of 4 km (N.B. to only 3 km in the 3-km-thick batholith case) has been modeled. The temperature versus depth for uniform heat generation of 4, 6, 8, 10, and 12 µW m−3 within the granitoids is calculated for batholith thicknesses of 3, 4.5, and 6 km. The results are shown in Figure 14 with *L*= 3 km (Figure 14a), *L*= 4.5 km (Figure 14b) and *L*= 6 km (Figure 14c). The thermal conductivity is assumed to be uniform at 3 W m−1 °C−1 [36]. The red curve also shown in Figure 14 represents the temperature versus depth in the absence of the radiogenic granite layer. The temperature in this case is estimated from a constant background heat flow of 24 mW m−2 (geothermal gradient of 8 ºC km−1). In all cases, the temperature gradient is less in the absence of the radiogenic granite. The modeling suggests that a maximum temperature of approximately 100 °C may occur at 3000 m within a granitoid body with depth extent of 6 km that possesses a uniform heat generation of 12 µW m−3. Certainly, such bodies have been observed from surface measurements, but the complete geometry and petrophysical makeup of these bodies are poorly understood.


**Table 3.** Parameters assumed for the temperature modeling for the Darling Range region for the seven cases investigated, wherein the temperature profile versus depth is shown in Figure 13.

Perhaps a more realistic situation may entail a granitoid of depth extent of 4.5 km (which is consistent with Jaeger [4] and uniform heat generation of 8 mW m−2). In this case, one is looking at almost 50 ºC at about 2000 m. Although this scenario cannot be used for most air conditioning and electricity uses, it is viable for recreational activities. It should be finally noted that this modeling is applicable only for the Darling Range vicinity of the Perth metropolitan area, which essentially comprises granitoid or latertic–granitoid-derived rocks. The Perth metro‐ politan area to the west of the Darling Range comprises over 15 km of Paleozoic sediments. It is presumed that the underlying rocks are similar to the Darling Range. This has yet to be conclusively confirmed.

## *3.2.2. Perth metropolitan area*

Temperature modeling for the Darling Range granitoids [16], and Table 3, was carried out using the geometry shown in Figure 3 and based upon Eq. (3). The purpose of the modeling is to determine if temperatures in the metropolitan Darling Range are high enough to support geothermal applications. In this modeling, the assumed parameters are the thicknesses of the radiogenic granite batholith (*L*) varying between 3 and 6 km. The complete temperature profile from surface (*z*=0) to a maximum depth of 4 km (N.B. to only 3 km in the 3-km-thick batholith case) has been modeled. The temperature versus depth for uniform heat generation of 4, 6, 8, 10, and 12 µW m−3 within the granitoids is calculated for batholith thicknesses of 3, 4.5, and 6 km. The results are shown in Figure 14 with *L*= 3 km (Figure 14a), *L*= 4.5 km (Figure 14b) and *L*= 6 km (Figure 14c). The thermal conductivity is assumed to be uniform at 3 W m−1 °C−1 [36]. The red curve also shown in Figure 14 represents the temperature versus depth in the absence of the radiogenic granite layer. The temperature in this case is estimated from a constant background heat flow of 24 mW m−2 (geothermal gradient of 8 ºC km−1). In all cases, the temperature gradient is less in the absence of the radiogenic granite. The modeling suggests that a maximum temperature of approximately 100 °C may occur at 3000 m within a granitoid body with depth extent of 6 km that possesses a uniform heat generation of 12 µW m−3. Certainly, such bodies have been observed from surface measurements, but the complete

geometry and petrophysical makeup of these bodies are poorly understood.

CASE 1 CASE 2 CASE 3 CASE 4 CASE 5 CASE 6 CASE 7

Ao (µWm-3) 4 10 20 4 10 20 8 K (Wm-1K-1) 3 3 3 3 3 3 3

L (km) 2 2 2 6 6 6 3.5 Qb (mWm-2) 26.3 26.3 26.3 26.3 26.3 26.3 26.3 Ts (°C) 20 20 20 20 20 20 20

Perhaps a more realistic situation may entail a granitoid of depth extent of 4.5 km (which is consistent with Jaeger [4] and uniform heat generation of 8 mW m−2). In this case, one is looking at almost 50 ºC at about 2000 m. Although this scenario cannot be used for most air conditioning and electricity uses, it is viable for recreational activities. It should be finally noted that this modeling is applicable only for the Darling Range vicinity of the Perth metropolitan area, which essentially comprises granitoid or latertic–granitoid-derived rocks. The Perth metro‐ politan area to the west of the Darling Range comprises over 15 km of Paleozoic sediments. It is presumed that the underlying rocks are similar to the Darling Range. This has yet to be

**Table 3.** Parameters assumed for the temperature modeling for the Darling Range region for the seven cases

**PARAMETER VALUE**

investigated, wherein the temperature profile versus depth is shown in Figure 13.

conclusively confirmed.

66 Advances in Geothermal Energy

The concept of recovering geothermal energy from the Yilgarn Craton is significantly enhanced by the presence of the city of Perth, the capital city of Western Australia, located near the western margin of the craton. The Darling Range is a physiographic expression of the Yilgarn Craton on the eastern edge of the Perth metropolitan area. Before 2010, only a few values of heat generation were available in the literature for the Yilgarn Craton [4,16]. However, studies have been carried out by the Western Australian Geothermal Centre of Excellence (WAGCoE) from 2009 to 2013 [42]. This centre carried out extensive studies on the potential of geothermal energy in Western Australia, and perhaps culminated in the concept of the "geothermal city" by Regenauer-Lieb et al. [43] and geothermal modeling of the Perth Basin by Schilling et al. [44].

The interest in geothermal energy was initially directed toward electricity generation. How‐ ever, Ballesteros [45] noted that there had been steady growth in the direct use of geothermal heat in Western Australia, especially in the Perth metropolitan area. He reported that 12 sites were using, or were in development, geothermal heat for community or industrial projects. Pujol et al. [46] and Ricard and Pujol [47] have reviewed the previous 20 years of exploitation of the Yarragadee aquifer in the Perth Basin for direct-use geothermal heat. These authors have indicated how the face of geothermal energy in Western Australia has changed from a perception of electricity creation to direct heat energy extraction.

## **3.3. Northern Perth basin**

The northern Perth Basin is bounded by the Northampton Complex to the north and the Yilgarn Craton to the east. Reasonable U and Th data are known from the Northampton Complex and airborne data are available for the Yilgarn Craton to the east. The Northampton Complex (see Figure 1 for location), which is part of the Pinjarra Orogen, has been mapped on the Geological Survey of Western Australia's Geraldton and Ajana 1:250,000 geological sheets [48]. The Northampton Complex forms the pre-sedimentary "basement" of the northern part of the Perth Basin, which is considered to possess appreciable geothermal energy potential, due to high temperatures observed within the sediments of the adjacent basin [12,50]. The high temperatures are presumed to be caused by elevated heat generation in the underlying Northampton Complex rocks.

Little recent geological work has been carried on the Northampton Complex. The Northamp‐ ton Complex (see Figures 1 and 2 for location) has been mapped on the Geraldton and Ajana 1:250,000 geological map sheets [48], and is part of the larger Pinjarra Orogen, which lies along the western margin of onshore Western Australia. On the Ajana 1:250,000 geological sheet, the Northampton Complex is dominated by rocks described as granulite and gneiss. These rocks have been metamorphosed under pressures in the range of 600–900 MPa (6–9 kBar) and temperatures in the range of 600– 800 °C (upper amphibolite-to-granulite facies metamorphic conditions; Sanders and McGuiness [50]). These rocks are not dissimilar to those of the Leeuwin Complex, which is also considered to be part of the Pinjarra Orogen. The merged thorium aero-radiometric map (Figure 2) suggests much of the Northampton Complex region, which consists of rocks with greater than 19.5 ppm thorium, noting that the proposed upper cutoff for the aero-radiometric maps for thorium is 19.5 ppm, as discussed above. However, also mentioned above and evident from ground-based studies in the Leeuwin Complex and Darling Range, one should be aware that 19.5 ppm appears to be the limit on the calibration technique of aero-radiometric observations for thorium, rather than a peak value that might be observed at the earth's surface. Laboratory-based values of concentration reported by Sanders and McGuiness [50] suggest that some localities may be significantly higher (i.e., > 39 ppm) than the aero-radiometric threshold maximum of 19.5 ppm.

Sanders and McGuiness [50] reviewed the geochemistry of the Ajana 1:250,000 geological sheet, and therein presented the valuable data concerning the surface distribution of thorium in the regolith. For the purposes of this study, the regolith is defined as unconsolidated or indurated weathered rock, and includes residual and transported material that can cover and obscure the underlying bedrock. The thorium concentrations reported in this review are from the Geological Survey of Western Australia map of Sanders and McGuiness [50], and the values were measured by the recognized laboratory methods referenced in the Sanders and McGui‐ ness' study. Their map is reproduced as Figure 15 in this chapter. The study presented a systematic grid of surface thorium concentrations in outcropping rocks at a spacing of approximately 5–10 km, nominally one sample per 16 km2 . The map of Sanders and McGuiness displays a range of thorium concentration between approximately 5 and > 39 ppm. Unfortu‐ nately, the upper limit of thorium concentration in the study, although acknowledged to be greater than 39 ppm, is not reported in their report.

Rocks from the Northampton Complex are very likely to extend under the north Perth Basin, and they probably contribute to the high temperatures observed in petroleum wells in that region [12,35]. Little definitive work has been carried out to date on the influence of radiogenic elements in the basement rocks upon elevated temperatures in the sediments of the north Perth Basin. Basement rocks have also been intersected in a number of petroleum wells. Table 4 shows an analysis of the gamma-ray log response found in the basement rocks in this part of the basin. The table also shows the gamma-ray log response for high-organic shale formations and some sandstones in the observation well for comparison with the basement rocks. Equation 8 is used to convert gamma-ray (GR) data from the wells to heat generation (*A*0) data. Heat generation in basement rocks determined by this method is seen to be in the range of 2.0– 5.6 µWm−3.

Thorium values mapped on the Ajana 1:250,000 geological sheet from Sanders and McGuiness [50] reflect relative heat generation as expressed in Eq. (4). Surprisingly, high Th is also found in the sedimentary rocks on the margin of the northern Perth Basin, but it does not appear to be especially dominant in the surrounding igneous and metamorphic rocks. The provenance of high radiogenic elements in sediments of the northern Perth Basin (perhaps also the southern Carnarvon Basin) is currently an imponderable. Further studies are needed to fully understand the geology, geochemistry, and geophysics of the northern Perth Basin and Northampton districts.


which consists of rocks with greater than 19.5 ppm thorium, noting that the proposed upper cutoff for the aero-radiometric maps for thorium is 19.5 ppm, as discussed above. However, also mentioned above and evident from ground-based studies in the Leeuwin Complex and Darling Range, one should be aware that 19.5 ppm appears to be the limit on the calibration technique of aero-radiometric observations for thorium, rather than a peak value that might be observed at the earth's surface. Laboratory-based values of concentration reported by Sanders and McGuiness [50] suggest that some localities may be significantly higher (i.e., > 39

Sanders and McGuiness [50] reviewed the geochemistry of the Ajana 1:250,000 geological sheet, and therein presented the valuable data concerning the surface distribution of thorium in the regolith. For the purposes of this study, the regolith is defined as unconsolidated or indurated weathered rock, and includes residual and transported material that can cover and obscure the underlying bedrock. The thorium concentrations reported in this review are from the Geological Survey of Western Australia map of Sanders and McGuiness [50], and the values were measured by the recognized laboratory methods referenced in the Sanders and McGui‐ ness' study. Their map is reproduced as Figure 15 in this chapter. The study presented a systematic grid of surface thorium concentrations in outcropping rocks at a spacing of

displays a range of thorium concentration between approximately 5 and > 39 ppm. Unfortu‐ nately, the upper limit of thorium concentration in the study, although acknowledged to be

Rocks from the Northampton Complex are very likely to extend under the north Perth Basin, and they probably contribute to the high temperatures observed in petroleum wells in that region [12,35]. Little definitive work has been carried out to date on the influence of radiogenic elements in the basement rocks upon elevated temperatures in the sediments of the north Perth Basin. Basement rocks have also been intersected in a number of petroleum wells. Table 4 shows an analysis of the gamma-ray log response found in the basement rocks in this part of the basin. The table also shows the gamma-ray log response for high-organic shale formations and some sandstones in the observation well for comparison with the basement rocks. Equation 8 is used to convert gamma-ray (GR) data from the wells to heat generation (*A*0) data. Heat generation in basement rocks determined by this method is seen to be in the range of 2.0–

Thorium values mapped on the Ajana 1:250,000 geological sheet from Sanders and McGuiness [50] reflect relative heat generation as expressed in Eq. (4). Surprisingly, high Th is also found in the sedimentary rocks on the margin of the northern Perth Basin, but it does not appear to be especially dominant in the surrounding igneous and metamorphic rocks. The provenance of high radiogenic elements in sediments of the northern Perth Basin (perhaps also the southern Carnarvon Basin) is currently an imponderable. Further studies are needed to fully understand the geology, geochemistry, and geophysics of the northern Perth Basin and

. The map of Sanders and McGuiness

ppm) than the aero-radiometric threshold maximum of 19.5 ppm.

approximately 5–10 km, nominally one sample per 16 km2

greater than 39 ppm, is not reported in their report.

5.6 µWm−3.

68 Advances in Geothermal Energy

Northampton districts.

**Reference lithology code for Table 4**. HS: Holmwood Shale (Permian); PCS: Permian coals and shale; IRCM: Irwin River Coal Measures; V: Volcanic intrusives (? Cretaceous); K: Kockatea Shale (Triassic); TS: Tumblagooda Sandstone; NCR: No Petrophysical Data (for gamma-ray). See Mory and Iasky (1996) for reference to stratigraphic units.

**Table 4.** Northern Perth Basin with U, Th, K & Ao derived from petroleum well data. The figure shows showing gamma-ray log (GR) responses from basement rocks with established crystalline basement intersections and derived sediments. Reference gamma-ray units (API) within the overlying formations for various wells in the Perth Basin are shown for confidence in gamma-ray logging in API units. The heat generation calculation, derived from observed logderived GR data is Ao = 0.015 (GR) (see text for explanation).

## **4. Pilbara region**

The Pilbara region of Western Australia contains some of the state's most significant mineral deposits. However, difficulties can arise for the development of these mineral resources due to the remoteness and lack of economic energy supply. The energy supply for mines at remote locations in this region is commonly sourced by diesel-powered generators, although some mines are supplied with gas from the North West Shelf. Geothermal energy may provide an economic alternative, or additional energy supply, if access to deep radiogenic granitoids is available. Hot Dry Rock geothermal systems are one of three ways that are considered viable for using geothermal energy in Australia. The three systems are (1) direct heat extraction, (2) Hot Sedimentary Aquifers (HSA) energy conversion, and (3) the Hot Dry Rock geothermal regime [12,51]. The HSA system, as well as direct heat extraction, is currently being explored by a number of companies as a possible system for both power generation and air conditioning in Western Australia. An HDR system is being exploited in the Cooper Basin by the geothermal company Geodynamics Limited. The HDR system entails injecting cool water into hot fractured rocks at depths of about 4 km, allowing the water to percolate through the hot rocks (which heats the water), and then extracting the hot water for geothermal energy applications. Geothermal explorers commonly seek temperatures in excess of 200 ºC for such systems. The hot water is fed into various types of turbines for electricity generation; the efficiency of electricity generation is dependent on the temperature of the fluid being introduced into the turbine.

Two regions of high radiogenic heat generation can be seen in Figure 2. The two hot regions are identified to be within the Numbana and Coolegong mesogranite intrusions. These intrusions are late-stage, highly fractionated granites [52,53] known to possess high U and Th concentrations, which crop out at the surface. Table 1 shows U and Th versus radiogenic heat generation (*A*0) of samples of the Numbana and Coolegong granites (data courtesy of the Geological Survey of Western Australia). An excellent linear trend has been found between U and heat generation (*A*0). The linear regression curve through U and *A*<sup>0</sup> data from Pilbara and Yilgarn Craton granitoids exhibit a simple relationship [16]:

$$\mathcal{U} = \mathcal{Z}\mathcal{Z}A\_0 - \text{5.4},\tag{9}$$

where *U* is the uranium concentration (ppm) and *A*<sup>0</sup> is the radiogenic heat generation (µWm −3); the coefficient of linear regression (*R*<sup>2</sup> ) is 0.92. As mentioned previously, the U and Th (essentially reflecting heat generation) of the rock samples are generally higher than the values observed from the airborne data [16]. These data suggest that the radiogenic heat generation in the mesogranites of the Pilbara Craton may be as high as 13 µWm−3, which is in the vicinity of those observed in the Leeuwin Complex in south west part of Western Australia.

Figure 18 shows temperature models displaying expected temperature versus depth within the two granitoids in the Pilbara region. The numbers on each curve represent the thickness (in km) of a granitoid with a uniform heat generation of 12 µW m−3. Temperatures at depths within the granitoids of 3–4 km may be over 100 ºC, but these require drilling with largediameter bore holes to those depths through crystalline rocks. Currently, this does not present an economic development.

New generations of turbines, for example, Kalina cycle and variable phase turbines [54,55], are allowing the use of lower-temperature geofluids for electricity generation. A stand-alone variable phase turbine has been designed to generate about 1 MW of electricity with input fluid temperatures of about 110 ºC, although greater generation efficiencies are obtained at higher temperatures. These new generation of turbines may hold significant implications for generating power at Pilbara mining sites. The key ingredient is the proximity of a hot granite to a mine location. Figure 2 shows the approximate radiogenic heat generation of surface rocks on the Marble Bar 1:250,000 map sheet. The map is generated from calibrated uranium, thorium, and potassium elemental abundances determined from airborne radiometric surveys [9,34]. The U, Th, and K abundances are converted to approximate radiogenic heat generation by a method described by Kappelmeyer and Haenel [23].

## **5. Carnarvon and Canning basins**

Apart from Perth Basin, HAS geothermal resources have been shown potentially to exist in the Carnarvon and Canning basins [12,51]. Some further work has been carried out for this review on the Carnarvon Basin, which explores, in a preliminary sense, the anomalous elevated temperature regime in the southern Carnarvon Basin. For the Canning Basin (Figure 1), little new geothermal insight has been gained since the study of Hot Dry Rocks [57], which was commissioned by the Geological Survey of Western Australia. Beardsmore and Cull [7] have reported some important heat generation data derived from gamma-ray logs in that basin. Their approach has been used earlier in this chapter to attempt to understand heat generation and the presence of radiogenic elements. Given the lack of new data from the Canning Basin, only the Carnarvon Basin will be discussed further.

## **5.1. Carnarvon basin**

**4. Pilbara region**

70 Advances in Geothermal Energy

turbine.

The Pilbara region of Western Australia contains some of the state's most significant mineral deposits. However, difficulties can arise for the development of these mineral resources due to the remoteness and lack of economic energy supply. The energy supply for mines at remote locations in this region is commonly sourced by diesel-powered generators, although some mines are supplied with gas from the North West Shelf. Geothermal energy may provide an economic alternative, or additional energy supply, if access to deep radiogenic granitoids is available. Hot Dry Rock geothermal systems are one of three ways that are considered viable for using geothermal energy in Australia. The three systems are (1) direct heat extraction, (2) Hot Sedimentary Aquifers (HSA) energy conversion, and (3) the Hot Dry Rock geothermal regime [12,51]. The HSA system, as well as direct heat extraction, is currently being explored by a number of companies as a possible system for both power generation and air conditioning in Western Australia. An HDR system is being exploited in the Cooper Basin by the geothermal company Geodynamics Limited. The HDR system entails injecting cool water into hot fractured rocks at depths of about 4 km, allowing the water to percolate through the hot rocks (which heats the water), and then extracting the hot water for geothermal energy applications. Geothermal explorers commonly seek temperatures in excess of 200 ºC for such systems. The hot water is fed into various types of turbines for electricity generation; the efficiency of electricity generation is dependent on the temperature of the fluid being introduced into the

Two regions of high radiogenic heat generation can be seen in Figure 2. The two hot regions are identified to be within the Numbana and Coolegong mesogranite intrusions. These intrusions are late-stage, highly fractionated granites [52,53] known to possess high U and Th concentrations, which crop out at the surface. Table 1 shows U and Th versus radiogenic heat generation (*A*0) of samples of the Numbana and Coolegong granites (data courtesy of the Geological Survey of Western Australia). An excellent linear trend has been found between U and heat generation (*A*0). The linear regression curve through U and *A*<sup>0</sup> data from Pilbara and

where *U* is the uranium concentration (ppm) and *A*<sup>0</sup> is the radiogenic heat generation (µWm

(essentially reflecting heat generation) of the rock samples are generally higher than the values observed from the airborne data [16]. These data suggest that the radiogenic heat generation in the mesogranites of the Pilbara Craton may be as high as 13 µWm−3, which is in the vicinity

Figure 18 shows temperature models displaying expected temperature versus depth within the two granitoids in the Pilbara region. The numbers on each curve represent the thickness (in km) of a granitoid with a uniform heat generation of 12 µW m−3. Temperatures at depths

of those observed in the Leeuwin Complex in south west part of Western Australia.

<sup>0</sup> *U A* = - 2.2 5.4, (9)

) is 0.92. As mentioned previously, the U and Th

Yilgarn Craton granitoids exhibit a simple relationship [16]:

−3); the coefficient of linear regression (*R*<sup>2</sup>

High temperatures have been recognized in the Southern Carnarvon Basin (Figure 1). These elevated temperatures have been reported in various studies [12,51,56]. This region is probably the most prospective region in Western Australia for geothermal energy, behind the Vasse Shelf region in the southern Perth Basin. However, very little geothermal exploration has been carried out in the region. However, substantial evidence supports the premise that the Carnarvon Basin may be able to supply the Pilbara mining communities with not only (1) longterm geothermal-sourced electricity (looking beyond North West Shelf-sourced gas supply, i.e., beyond perhaps 2050), but also (2) intermediate-term nearby resource projects and tourism enterprises. It is of note that Davidson [58,59] described the Peninsular Hot Spring tourist enterprise in southern Victoria, which entertains an estimated \$75 million per annum regional economic benefit. Such enterprises are potentially available to environmentally sensitive, but tourist intensive, coastal resort sites from Carnarvon north to Exmouth.

## *5.1.1. Geology and geophysics*

The geology of the southern Carnarvon Basin is well understood, having been described by Hocking et al. [60], Hocking [61], and the geothermal potential by Ghori [12]. The geology of the North West Cape, which is especially relevant to the current article, was studied by Malcolm et al. [62]. It has been shown from petroleum drilling data that an unusually hot region exists along the coastal region from Shark Bay to the Exmouth Peninsula, and maybe for some distance inland toward, and perhaps including, the outcropping Pilbara Craton. It is reason‐ able to expect that the Pilbara Craton underlies a reasonable amount of the onshore Carnarvon Basin between Shark Bay and Exmouth, including the North West Cape. Further, it is expected that the Pilbara Craton may contain some radiogenic granitoids in this region. There appears to be no detailed scientific study of why the elevated heat flow occurs in this region, despite it being covered since 2010 by seven GEPs. According to the studies by Ghori [12,51], temper‐ atures in excess of 150 ºC may exist between drillable depths of 3000–4000 m, and these wells have a geothermal gradient in the range of 60–100 ºC/km (figure 2 from Ghori [12]).

## *5.1.2. A geothermal puzzle*

The unanswered geothermal puzzle of the southern Carnarvon Basin is centered on the cause of the elevated geothermal gradients in the Exmouth and Gascoyne subbasins. Essentially, the puzzle may condense into whether the elevated temperatures reported by Ghori [12,51] are caused by underlying radiogenic granitoids or by the flow of deep hot groundwater into surface sediments. No extensive studies have been carried out to understand this phenomenon. The source of the observed high heat flow in petroleum wells appears to be generated from high radiogenic granitoids that underlie the deep sedimentary blanket in the region. However, this is only a circumstantial evidence, and without sufficient further investigation, it cannot be stated with certainty. Convective heat transport should not be ruled out as the cause of these elevated temperatures. Comparison to the surface geological units to elsewhere in Western Australia suggests that the sedimentary cover may have a surface heat generation in the vicinity of 1 µWm−3 (see Table 4, and Beardsmore and Cull [7]). On the adjacent Pilbara Craton, outcropping granitoids possess heat generation values of up to 5 µWm−3 in the adjacent Pilbara Craton; these may translate up to 12 µWm−3 in specific geologic units on the ground. Therefore, the granitoids that are interpreted to underlie the Carnarvon Basin in the cross section in Figure 19 may possess a heat generation of up to 12 µWm−3 or more. This heat generation value is typical of radiogenic granitoids underlying the Perth Basin, and also the central Australian Cooper Basin [6], where geothermal energy is currently being developed. The WAPET Cape Range petroleum exploration wells are known to encounter bottom-hole temperatures in excess of 140–160 ºC at depths between 4300 and 4600 m; these temperatures are probably greater after fluid stabilization in the drill hole. Figure 19 shows that the total sedimentary blanket over interpreted high heat-generating granites may be up to 9 km [61,62].

#### *5.1.3. Temperature modeling*

Based on the cross section [56] shown in Figure 19, several numerical models of the temperature versus depth for various heat-generation scenarios are shown in Table 5a. These are based on the previously described 1D models for such situations as described in Eqs (1) and (2). The assumptions used to generate Table 5a are shown in Table 5b, and they are based on the surface radiogenic heat data derived from Figure 2 and Eq. (4), and also studies reported elsewhere in the literature [4,5,16,61]. Ghori [12] has shown from petroleum wells that the temperature gradient in the southern Carnarvon Basin can range between 30 and 100 ºC/km, and these are consistent with the temperatures generated in the models in Table 1. If convection is not involved, then the governing factor appears to be the heat generation of underlying granites. However, a disturbing observation from hydrogeological studies is that there are significant indications of hot artesian water movement in this basin. The static models in Table 5a demonstrate that temperatures at depths greater than 3000 m may have temperatures signif‐ icantly greater than 200 ºC, depending on the heat generation of underlying granite bodies. In the upper side of the models, one is looking at full steam-type turbines for the generation of electricity, which are the current optimal geothermal energy-producing technology, and is common place in New Zealand, Indonesia, and California. A caveat needs to be placed on the geothermal activity in this basin centered around the balance of radiogenic heat input versus convective heat input. More investigation is required.


**Table 5.** Carnarvon Basin modeling results (a) temperatures and (b) assumptions. See also the geological cross-section in Figure 19 for geological consideration discussed in the text.

## **6. Time-variant considerations**

*5.1.1. Geology and geophysics*

72 Advances in Geothermal Energy

*5.1.2. A geothermal puzzle*

*5.1.3. Temperature modeling*

The geology of the southern Carnarvon Basin is well understood, having been described by Hocking et al. [60], Hocking [61], and the geothermal potential by Ghori [12]. The geology of the North West Cape, which is especially relevant to the current article, was studied by Malcolm et al. [62]. It has been shown from petroleum drilling data that an unusually hot region exists along the coastal region from Shark Bay to the Exmouth Peninsula, and maybe for some distance inland toward, and perhaps including, the outcropping Pilbara Craton. It is reason‐ able to expect that the Pilbara Craton underlies a reasonable amount of the onshore Carnarvon Basin between Shark Bay and Exmouth, including the North West Cape. Further, it is expected that the Pilbara Craton may contain some radiogenic granitoids in this region. There appears to be no detailed scientific study of why the elevated heat flow occurs in this region, despite it being covered since 2010 by seven GEPs. According to the studies by Ghori [12,51], temper‐ atures in excess of 150 ºC may exist between drillable depths of 3000–4000 m, and these wells

have a geothermal gradient in the range of 60–100 ºC/km (figure 2 from Ghori [12]).

blanket over interpreted high heat-generating granites may be up to 9 km [61,62].

Based on the cross section [56] shown in Figure 19, several numerical models of the temperature versus depth for various heat-generation scenarios are shown in Table 5a. These are based on

The unanswered geothermal puzzle of the southern Carnarvon Basin is centered on the cause of the elevated geothermal gradients in the Exmouth and Gascoyne subbasins. Essentially, the puzzle may condense into whether the elevated temperatures reported by Ghori [12,51] are caused by underlying radiogenic granitoids or by the flow of deep hot groundwater into surface sediments. No extensive studies have been carried out to understand this phenomenon. The source of the observed high heat flow in petroleum wells appears to be generated from high radiogenic granitoids that underlie the deep sedimentary blanket in the region. However, this is only a circumstantial evidence, and without sufficient further investigation, it cannot be stated with certainty. Convective heat transport should not be ruled out as the cause of these elevated temperatures. Comparison to the surface geological units to elsewhere in Western Australia suggests that the sedimentary cover may have a surface heat generation in the vicinity of 1 µWm−3 (see Table 4, and Beardsmore and Cull [7]). On the adjacent Pilbara Craton, outcropping granitoids possess heat generation values of up to 5 µWm−3 in the adjacent Pilbara Craton; these may translate up to 12 µWm−3 in specific geologic units on the ground. Therefore, the granitoids that are interpreted to underlie the Carnarvon Basin in the cross section in Figure 19 may possess a heat generation of up to 12 µWm−3 or more. This heat generation value is typical of radiogenic granitoids underlying the Perth Basin, and also the central Australian Cooper Basin [6], where geothermal energy is currently being developed. The WAPET Cape Range petroleum exploration wells are known to encounter bottom-hole temperatures in excess of 140–160 ºC at depths between 4300 and 4600 m; these temperatures are probably greater after fluid stabilization in the drill hole. Figure 19 shows that the total sedimentary

> It may be worth examining some time-variant models for heat flow and temperature variation for Western Australian geothermal scenarios. These are largely based on applications of Eq.

(3) [24], which considers the case where a layer of rock with high heat generation is rapidly placed at the earth's surface at a time *t* = 0. A typical situation for Western Australia may occur where lateritic enrichment occurs in the Tertiary Period. The question being posed is: does this enrichment in U and Th alter the temperature profile in the vicinity of the chemical alteration process, and, if so, how great is the temperature disturbance and how long does it take to stabilize? An initial intuitive guess would be that stabilization occurs fairly rapidly and there is not much disturbance to the temperature or heat flow regime. A lateritic profile in the Darling Range may be up to several hundred meters thick [40], but generally less than 100 m. A weathering profile on the Leeuwin Complex may be of similar thickness [63]. The question posed revolves around the difference that occurs for a heat-generating profile of thickness between 50 and 300 m at the surface. We can use the Meelup Park enrichment value (Table 1) of 40 µW m−3 as an initial test.

Carslaw and Jaeger [24] propose that the surface geothermal gradient for a geometry shown in Figure 3, where a heat-generating layer (injection of fluids or erosional enrichment) occurs rapidly, can be described by an equation of the type:

$$
\Delta Q\_{\aleph} = \left. 2A\_0 \left( at \right)^{\natural\_{\natural}} \left| \left< \left( \pi \right)^{\natural\_{\natural}} - \text{ierfc} \left[ L \left< \left( \mathbf{4} at \right)^{\natural\_{\natural}} \right] \right> \right| \right. \tag{10}
$$

where Δ*QR* is the heat flow contribution of the radiogenic layer in units of mW m−2, *A*0 is the uniform heat generation in the layer in µW/m3 , *α* is the thermal diffusivity of the rocks in m2 s−1, *t* is the time in seconds, *L* is the thickness of the radiogenic enrichment layer in m, and ierfc[*x*] is the integrated error function. This equation has been modified from the Carslaw and Jaeger [24] equation 6 by multiplication by thermal conductivity (*K)* to be expressed as heat flow rather than thermal gradient. Carslaw and Jaeger also indicate that this equation will reduce to the simple form of Δ*QR* = *A*<sup>0</sup> *L*, where time becomes very large (> 10 million years) and *L* is less than 50 km. This was referred to earlier with regard to Eq. (1).

Some simple calculations with Eq. (9) indicate that stabilization of an enrichment layer, such as observed at Meelup Park (see Table 1 and Figure 9 for magnitude of Th anomaly) in the Leeuwin Complex and the Vasse Shelf, southern Perth Basin, will occur well within 1 million years of enrichment, where a 100 m of enrichment zone is assumed. A zone of 100 m may be expected, based on geological and environmental mapping [63]. Similar durations for thermal stabilization may also be expected for the lateritic zones in the Darling Range in the central Perth Basin, although such high levels of enrichment are not observed from current observa‐ tions (see Table 1). These conclusions, based on numerical modeling, are consistent with expectations. In the case where Δ*QR* = *A*<sup>0</sup> *L*, which is the stabilized case above, it is interesting to understand the influence that an enriched zone may have on heat flow. Table 6 shows a range of incremental heat flow (Δ*QR*) values that will occur for a range of heat generation (*A*0) and enrichment layer thicknesses (*L*). The table indicates that enrichment zones may contribute an appreciable input to localized heat flow, especially if the thickness is greater than 100 m and heat generation is greater than 80 µW m−3.

(3) [24], which considers the case where a layer of rock with high heat generation is rapidly placed at the earth's surface at a time *t* = 0. A typical situation for Western Australia may occur where lateritic enrichment occurs in the Tertiary Period. The question being posed is: does this enrichment in U and Th alter the temperature profile in the vicinity of the chemical alteration process, and, if so, how great is the temperature disturbance and how long does it take to stabilize? An initial intuitive guess would be that stabilization occurs fairly rapidly and there is not much disturbance to the temperature or heat flow regime. A lateritic profile in the Darling Range may be up to several hundred meters thick [40], but generally less than 100 m. A weathering profile on the Leeuwin Complex may be of similar thickness [63]. The question posed revolves around the difference that occurs for a heat-generating profile of thickness between 50 and 300 m at the surface. We can use the Meelup Park enrichment value (Table 1)

Carslaw and Jaeger [24] propose that the surface geothermal gradient for a geometry shown in Figure 3, where a heat-generating layer (injection of fluids or erosional enrichment) occurs

> ( ) { ( ) ( ) } ½ ½ <sup>½</sup> 0 2 1 ierfc 4 , *Q At <sup>R</sup>*

where Δ*QR* is the heat flow contribution of the radiogenic layer in units of mW m−2, *A*0 is

Some simple calculations with Eq. (9) indicate that stabilization of an enrichment layer, such as observed at Meelup Park (see Table 1 and Figure 9 for magnitude of Th anomaly) in the Leeuwin Complex and the Vasse Shelf, southern Perth Basin, will occur well within 1 million years of enrichment, where a 100 m of enrichment zone is assumed. A zone of 100 m may be expected, based on geological and environmental mapping [63]. Similar durations for thermal stabilization may also be expected for the lateritic zones in the Darling Range in the central Perth Basin, although such high levels of enrichment are not observed from current observa‐ tions (see Table 1). These conclusions, based on numerical modeling, are consistent with expectations. In the case where Δ*QR* = *A*<sup>0</sup> *L*, which is the stabilized case above, it is interesting to understand the influence that an enriched zone may have on heat flow. Table 6 shows a range of incremental heat flow (Δ*QR*) values that will occur for a range of heat generation (*A*0) and enrichment layer thicknesses (*L*). The table indicates that enrichment zones may contribute an appreciable input to localized heat flow, especially if the thickness is greater than 100 m

s−1, *t* is the time in seconds, *L* is the thickness of the radiogenic enrichment layer in m, and ierfc[*x*] is the integrated error function. This equation has been modified from the Carslaw and Jaeger [24] equation 6 by multiplication by thermal conductivity (*K)* to be expressed as heat flow rather than thermal gradient. Carslaw and Jaeger also indicate that this equation will reduce to the simple form of Δ*QR* = *A*<sup>0</sup> *L*, where time becomes very large (> 10 million years) and *L* is less than 50 km. This was referred to earlier with regard to

 a*L t* é ù D = - ê ú ë û (10)

, *α* is the thermal diffusivity of the rocks

ap

of 40 µW m−3 as an initial test.

74 Advances in Geothermal Energy

in m2

Eq. (1).

rapidly, can be described by an equation of the type:

the uniform heat generation in the layer in µW/m3

and heat generation is greater than 80 µW m−3.

**Figure 9.** Thorium distribution in Dunsborough region (GSWA database).

**Figure 10.** Uranium distribution in Dunsborough region (GSWA database).

Radiogenic Heat Generation in Western Australia — Implications for Geothermal Energy http://dx.doi.org/10.5772/61963 77


**Table 6.** Modeling results for shallow enrichment zones.

**Figure 10.** Uranium distribution in Dunsborough region (GSWA database).

76 Advances in Geothermal Energy

**Figure 11.** Magnetic anomalies and Th anomaly locations for the Leeuwin Complex and vase Shelf (south Perth Basin). The Th and associated heat generation anomalies appear to lie in magnetic medium to low areas. The map shows the aeromagnetics in the vicinity of Dunsborough with the location of the Th anomalies superimposed on the left-handside (LHS) image. The right-hand-side (RHS) image shows the aeromagnetics without the Th anomalies superimposed, and the color scale of magnetic intensity in units of nanotesla (nT). The magnitude of the Th anomalies (outline in black on LHS) is represented as green indicating 50–100 ppm, orange indicating 100–150 ppm, and red indicating 150– 200 ppm. Ground measurements of Th indicate concentrations are in the order of 280 ppm in the orange regions.

**Figure 12.** Temperature versus depth for the seven cases proposed in Table 2 for the south Perth Basin (Vasse Shelf). The number on each curve identifies the appropriate case in Table 2. The Wharton [37] study suggests that the temper‐ ature at 1000-m depth is approximately 40 °C. Cases 2, 3, 4, and 7 most closely satisfy this observation.

**Figure 13.** Image of the central Perth Basin (in the vicinity of the Perth metropolitan area and adjacent Darling Range): with some measurement localities (modified from Middleton [16]).

Radiogenic Heat Generation in Western Australia — Implications for Geothermal Energy http://dx.doi.org/10.5772/61963 79

**Figure 12.** Temperature versus depth for the seven cases proposed in Table 2 for the south Perth Basin (Vasse Shelf). The number on each curve identifies the appropriate case in Table 2. The Wharton [37] study suggests that the temper‐

**Figure 13.** Image of the central Perth Basin (in the vicinity of the Perth metropolitan area and adjacent Darling Range):

with some measurement localities (modified from Middleton [16]).

ature at 1000-m depth is approximately 40 °C. Cases 2, 3, 4, and 7 most closely satisfy this observation.

78 Advances in Geothermal Energy

**Figure 14.** Temperature versus depth for the Darling Range [16]. The values shown in the circles are the assumed uni‐ form heat generation in the radiogenic (granitoid) layer in units of µW m−3. Figure 14(a) shows temperature versus depth with the radiogenic layer of depth extent of 3 km. Figure 14(b) shows temperature versus depth with the radio‐ genic layer of depth extent of 4.5 km. Figure 14(c) shows temperature versus depth with the radiogenic layer of depth extent of 6 km. The curves in red on 14a, 14b, and 14c represent the temperature versus depth profile, where there is no heat generation in the granitoid layer.

**Figure 15.** Thorium values Ajana 1:250,000 geological sheet from Sanders and McGuiness [50] are shown in this figure. The map also reflects the main contributing sources to heat generation in the vicinity of the northern Perth Basin and the rocks expected to underlie it. High Th is found in the sedimentary rocks on the margin of the northern Perth Basin, and it does not appear to be significantly less than the surrounding igneous and metamorphic rocks. The provenance of high radiogenic elements in sediments of the northern Perth Basin (perhaps also the southern Carnarvon Basin) is currently an imponderable. More work is definitely needed to understand the geology, geochemistry, and geophysics of the northern Perth Basin and Northampton district

Radiogenic Heat Generation in Western Australia — Implications for Geothermal Energy http://dx.doi.org/10.5772/61963 81

**Figure 16.** North Perth Basin well locations with basement intersections for potential heat generation estimations (com‐ bined with Table 4).

**Figure 15.** Thorium values Ajana 1:250,000 geological sheet from Sanders and McGuiness [50] are shown in this figure. The map also reflects the main contributing sources to heat generation in the vicinity of the northern Perth Basin and the rocks expected to underlie it. High Th is found in the sedimentary rocks on the margin of the northern Perth Basin, and it does not appear to be significantly less than the surrounding igneous and metamorphic rocks. The provenance of high radiogenic elements in sediments of the northern Perth Basin (perhaps also the southern Carnarvon Basin) is currently an imponderable. More work is definitely needed to understand the geology, geochemistry, and geophysics

of the northern Perth Basin and Northampton district

80 Advances in Geothermal Energy

**Figure 17.** Pilbara Region hot granites. A map showing the integration of airborne U, Th, and K concentrations to re‐ flect surface heat generation in the Pilbara region) [34]. The "hot" Numbana and Coolegong granitoid bodies are seen as the areas as shown in the image.

**Figure 18.** Pilbara temperatures for hot granites. The curves represent temperature versus depth for a radiogenic gran‐ itoid of uniform heat generation of 12 µW m−3. The numbers on the curves represent the thickness of the radiogenic (granitoid) layer.

**Figure 19.** Carnarvon Basin summary (cross section). The red represents interpreted granitoid rocks that may possess radiogenic heat potential.

**Figure 17.** Pilbara Region hot granites. A map showing the integration of airborne U, Th, and K concentrations to re‐ flect surface heat generation in the Pilbara region) [34]. The "hot" Numbana and Coolegong granitoid bodies are seen

**Figure 18.** Pilbara temperatures for hot granites. The curves represent temperature versus depth for a radiogenic gran‐ itoid of uniform heat generation of 12 µW m−3. The numbers on the curves represent the thickness of the radiogenic

as the areas as shown in the image.

82 Advances in Geothermal Energy

(granitoid) layer.

**Figure 20.** Figure showing known temperature regimes and possible geothermal energy development scenarios for Western Australia as of 2015.

## **7. Conclusion**

In conclusion of this study, there are a number of generally low-temperature geothermal energy applications effectively operating in Western Australia [45–47]. The principal challenge for geothermal energy is to compete in the gas-rich economy of Western Australia. Figure 20 shows the current range of scenarios for temperatures versus depth in Western Australia. Case 1 is the current best case, where at about 2000 m, approximately 105 °C is encountered. It is anticipated that it can happen, but is not expected to be common on current geological and geophysical knowledge. Developments in the cases 3–5 are well known and are currently happening, and are common in the Perth metropolitan area. The green area in Figure 20 (between curves 2 and 3) is the space of current geothermal development, and regions in Western Australia that may take advantage of this temperature/depth range may be in the tourist areas of the southwest of Western Australia and the Carnarvon Basin-Cape Range regions. The transition to electricity (curve 5 and beyond) is challenging in the current economic climate, but may be anticipated to occur within the next 10–15 years.

The review recognizes that temperatures between 3000 and 4000 m in the Vasse Shelf may fall in the range of 90–150 °C. This temperature range is sufficient to permit electricity generation with organic Rankine cycle (ORC) or variable phase cycle (VPC) turbines [54,55], albeit with low efficiency. This temperature range is ideal for direct heat use for the tourist industry, for which the region is well known. Heat generation in radiogenic granitoids in the Darling Range is observed to fall in the range of 4–10 µW m−3. The heat generation appears to be quite variable within the Darling Range region; however, there is some evidence that the heat generation may be fairly uniform within any particular batholith or location. In addition, there is a movement toward direct use for tourism and small-scale industrial usage. Significant geo‐ thermal projects are in the development stage for district air-conditioning, because Perth, which builds housing for a mild Mediterranean climate, experiences a wide temperature range from as low as 1 °C in winter to over 40 °C in summer. This reflects the direction of the current change in paradigm for geothermal in Western Australia.

Finally, in consideration of commercial geothermal energy production in the Western Aus‐ tralian region, the modeled temperatures cannot support current concepts of large-scale electricity generation with geothermal energy based on the *Hot-Dry-Rock-type* scenario, where it is generally accepted that temperatures are well in excess of 220 ºC.However, the new generation of low-temperature, organic-Rankine cycle turbines [54,55], and adsorption chillers, which operate in the vicinity of 100–150 ºC, could place the Vasse Shelf, Darling Range, and mining regions in the Pilbara in an advantageous economic position for alternative energy supply close to an existing market.

## **Acknowledgements**

The technical help and advice of Mark K. Stevens, Ameed Ghori, Klaus Gessner, Lynn Reid, and Kerem Kanadarikrik from the Western Australian Department of Mines and Petroleum is gratefully acknowledged for input into my geothermal studies over the past six years. The editorial skills of Karina Jonasson from the same organization have made a significant impact on the quality of the department's geothermal and other publications over these years. Encouragement from and discussions with, his friend and colleague, Mark Ballesteros, set the author on a unique direction within the geothermal path in 2009. This resulted in two extremely valuable symposia: Western Australian Geothermal Energy Symposium (WAGES) 2011 and (WAGES) 2012. Finally, the author wishes to acknowledge the Executive Director of the Petroleum Division in the Western Australian Department of Mines and Petroleum, Jeff Haworth, for permission to publish this work.

## **Author details**

**7. Conclusion**

84 Advances in Geothermal Energy

In conclusion of this study, there are a number of generally low-temperature geothermal energy applications effectively operating in Western Australia [45–47]. The principal challenge for geothermal energy is to compete in the gas-rich economy of Western Australia. Figure 20 shows the current range of scenarios for temperatures versus depth in Western Australia. Case 1 is the current best case, where at about 2000 m, approximately 105 °C is encountered. It is anticipated that it can happen, but is not expected to be common on current geological and geophysical knowledge. Developments in the cases 3–5 are well known and are currently happening, and are common in the Perth metropolitan area. The green area in Figure 20 (between curves 2 and 3) is the space of current geothermal development, and regions in Western Australia that may take advantage of this temperature/depth range may be in the tourist areas of the southwest of Western Australia and the Carnarvon Basin-Cape Range regions. The transition to electricity (curve 5 and beyond) is challenging in the current

economic climate, but may be anticipated to occur within the next 10–15 years.

change in paradigm for geothermal in Western Australia.

supply close to an existing market.

**Acknowledgements**

The review recognizes that temperatures between 3000 and 4000 m in the Vasse Shelf may fall in the range of 90–150 °C. This temperature range is sufficient to permit electricity generation with organic Rankine cycle (ORC) or variable phase cycle (VPC) turbines [54,55], albeit with low efficiency. This temperature range is ideal for direct heat use for the tourist industry, for which the region is well known. Heat generation in radiogenic granitoids in the Darling Range is observed to fall in the range of 4–10 µW m−3. The heat generation appears to be quite variable within the Darling Range region; however, there is some evidence that the heat generation may be fairly uniform within any particular batholith or location. In addition, there is a movement toward direct use for tourism and small-scale industrial usage. Significant geo‐ thermal projects are in the development stage for district air-conditioning, because Perth, which builds housing for a mild Mediterranean climate, experiences a wide temperature range from as low as 1 °C in winter to over 40 °C in summer. This reflects the direction of the current

Finally, in consideration of commercial geothermal energy production in the Western Aus‐ tralian region, the modeled temperatures cannot support current concepts of large-scale electricity generation with geothermal energy based on the *Hot-Dry-Rock-type* scenario, where it is generally accepted that temperatures are well in excess of 220 ºC.However, the new generation of low-temperature, organic-Rankine cycle turbines [54,55], and adsorption chillers, which operate in the vicinity of 100–150 ºC, could place the Vasse Shelf, Darling Range, and mining regions in the Pilbara in an advantageous economic position for alternative energy

The technical help and advice of Mark K. Stevens, Ameed Ghori, Klaus Gessner, Lynn Reid, and Kerem Kanadarikrik from the Western Australian Department of Mines and Petroleum

Mike F. Middleton\*

Address all correspondence to: mike.middleton@dmp.wa.gov.au

Department of Mines and Petroleum, East Perth, WA, Australia

## **References**


[8] Chopra, P.N, and Holgate, F., 2005, A GIS analysis of temperature in the Australian crust, *Proceedings of the World Geothermal Congress 2005*, Antalya, Turkey, 24–29 April

[9] Budd, A.R., 2011, Geothermal Energy Project, Mapping Heat Across the Australian

[10] Budd, A.R., Gerner, E., Meixner, A.J., Kirkby, A.I., Weber, R. and Spring, J., 2011, Ge‐ oscience Australia's onshore energy security program: geothermal energy, in: *Pro‐ ceedings of the 2011 Australian Geothermal Energy Conference*, 16–18 November,

[11] Bestow, T.T., 1982, The potential for geothermal energy in Western Australia, *Geologi‐*

[12] Ghori, K.A.R., 2009, Petroleum data: leading the search for geothermal resources in

[13] Reid, P. W., Messeiller M., Llanos E. M. and Hasting, M., 2011, Paralana 2 – Well Testing and Stimulation Australian Geothermal Energy Conference 2011, Geoscience

[14] Cull, J.P. and Denham, D., 1979, Regional variations in Australian heat flow, *BMR*

[15] Cull, J.P., 1982, An appraisal of Australian heat flow data, *BMR Journal of Australian*

[16] Middleton, M.F., 2013, Radiogenic heat generation in the Darling Range, Western

[17] Middleton, M. F. and Stevens, M.K.S., 2014, Heat Generation in the Vasse Shelf Re‐ gion: Implications for Geothermal Exploration, Petroleum in Western Australia,

[18] Middleton, M.F., 2014, Brett, J. and Flint, D., 2014, Thorium occurrence: geological and geophysical implications for Western Australia, Chapter 6. In: Thorium: Chemi‐ cal Properties, Uses and Environmental Effects, (R. Jamison, Ed.) Nova Science Pub‐

[19] Brett, J, 2013, *Thorium merge map of WA, Geological Survey of Western Australia,* elec‐

[20] Minty, B.R.S., 1997, Fundamentals of airborne gamma-ray spectrometry, AGSO *Jour‐*

[21] Minty, B., 2014, The Gamma-ray Spectrometric Method for Mineral Exploration and

Melbourne, Geoscience Australia, Record 2011/43, (Budd, A.R. Ed.), 31–36.

Continent, 3rd Hot Rock Energy Conference 2011.

*cal Survey of Western Australia*, Record 82/6. 67p.

Australia Record 2011/43, (A. Budd, Ed.), 193–196.

*Journal of Australian Geology and Geophysics*, 4, 1–13.

Australia, *Exploration Geophysics*, 44, 206–214.

lishers Inc., ISBN: 978-1-63321-309-8, pp. 133–178.

*nal of Australian Geology and Geophysics*, v. 17, No. 2, 39–50.

Environmental Mapping, ASEG OzSTEP Workshop2014.

*Geology and Geophysics*, 7, 11–21.

March issue 2014.

tronic publication.

Western Australia, *APPEA Journal*, 365–382.

2005, 1–7.

86 Advances in Geothermal Energy


18 November, Melbourne, Geoscience Australia, Record 2011/43, (A. Budd, Ed.), 77– 80.

[50] Sanders, A. J. and McGuiness, S. A., 2000, Geochemical mapping of the Ajana 1:250 000 map sheet: *Western Australia Geological Survey*, 1:250 000 Regolith Geochemistry Series Explanatory Notes, 55p.

[36] Schön, J.H., 1996, Physical Properties of Rocks: Fundamentals and Principles of Pet‐

[37] Wharton, P.H., 1982, The geology and hydrogeology of the Quindalup boreline, southern Perth Basin, Western Australia, *Geological Survey of Western Australia*, Re‐

[38] Playford, P. E., Cockbain, A. E. and Low, G. H., 1976, Geology of the Perth Basin,

[39] Wilde, S.A., Middleton, M.F. and Evans, B.J., 2006, Terrane accretion in the south‐ western Yilgarn Craton: evidence from a deep seismic crustal profile, *Precambrian Re‐*

[40] Annand, R.R. and Butt, C.R.M., 2003, Distribution and evolution of 'laterites' and lat‐ eritic weathering profiles, Darling Range, Western Australia, *Australian Geomechanics*,

[41] Middleton, M.F., 2010, The Gillingarra temperature anomaly, northern Perth Basin, Australia, *Expanded Abstracts, Proceedings of SEG International Exposition and 80th An‐*

[42] Regenauer-Lieb, K. and the Team of the Western Australian Centre of Excellence, 2011, Geothermal Cities, Australian Geothermal Energy Conference 2011, Geoscience

[43] Regenauer-Lieb, K. and the Team of the Western Australian Centre of Excellence, 2012, The Western Australian Geothermal Research Centre and its future, in: Middle‐ ton, M. and Gessner, K. (Eds) *West Australian Geothermal Energy Symposium Abstracts*,

[44] Schilling, O, Sheldon, H.A., Reid, L.B. and Corbel, S., 2013, Hydrothermal models of the Perth metropolitan area, Western Australia: implications for geothermal energy,

[45] Ballesteros, M., 2012, personal communication, concerning geothermal developments

[46] Pujol, M., Ricard, L., and Bolton, G., 2015, 20 years of exploitation of the Yarragadee aquifer in the Perth Basin of Western Australia for direct-use geothermal heat, *Geo‐*

[47] Ricard, L and Pujol, M., 2015, Geothermal production in the Perth Metropolitan Area, 20 years of success. Petroleum in Western Australia, September issue, 2015. [48] Hocking, R. M., van de Graaff., W. J. E., Blockley, J. G. and Butcher, B.P., 1982. Ajana, W.A.: *Western Australia Geological Survey*, 1:250 000 Geological Series Explanatory

[49] Ghori, K.A.R. and Gibson, H., 2011, Geothermal simulation, Perth Basin, Western Australia, in: 2011, *Proceedings of the 2011 Australian Geothermal Energy Conference*, 16–

Western Australia: *Western Australia Geological Survey*, *Bulletin* 124.

rophysics, Pergamon.

cord 1982/2.

88 Advances in Geothermal Energy

*search*, 78, 179–196.

v. 36, No. 4, 41–58.

*nual Meeting, Denver, Colorado*, October 2010.

Australia Record 2011/43, (A. Budd, Ed.), 185–188.

2 & 3 April 2012, v. 2, 23–24 (ISBN 978-1-74052-250-2.

*Hydrogeology Journal* 21 (3), 605–621.

in the Perth metropolitan area.

*thermics*, v. 57, 39–55.

Notes.


[62] Lennard, E.L., 1991, YALLINGUP SHEET 1930 IV and part Sheet 18301, Geological Survey of Western Australia. 1:50 000 Environmental Geology Series.
