222Rn solution by the circulating fluids in a “hot dry rock” geothermal reservoir

The mechanism of²²²Rn release into fracture fluids by direct alpha-recoil, lattice and grain boundary/micro-crack diffusion is discussed. Experimental measurements of²²²Rn release into surrounding air and water phases have been made for crystalline rock specimens with well defined surface areas. The²²²Rn flux from an infinite plane surface and hence the effective diffusion length of²²²Rn in the rock matrix has been estimated.The²²²Rn flux from plane crystalline rock surfaces has been used in conjunction with a simple hydrological model of the reservoir to calculate the²²²Rn content of the return fluids of a geothermal doublet circulation system. For given production rate and piezometric difference between the injection and production wells, the²²²Rn content of the return fluid is dependent upon the distribution of flow path lengths and fracture apertures in the reservoir. Matching of the calculated and experimental²²²Rn contents of the return fluids has been used to select appropriate parameters for the reservoir model and hence to estimate the extent of the heat-transfer surface. The model estimates the fracture width of the flow paths, total swept surface area and fracture volume within the reservoir.     

Rn solution by the circulating fluids in a “hot dry rock” geothermal reservoir

The mechanism of²²²Rn release into fracture fluids by direct alpha-recoil, lattice and grain boundary/micro-crack diffusion is discussed. Experimental measurements of²²²Rn release into surrounding air and water phases have been made for crystalline rock specimens with well defined surface areas. The²²²Rn flux from an infinite plane surface and hence the effective diffusion length of²²²Rn in the rock matrix has been estimated.The²²²Rn flux from plane crystalline rock surfaces has been used in conjunction with a simple hydrological model of the reservoir to calculate the²²²Rn content of the return fluids of a geothermal doublet circulation system. For given production rate and piezometric difference between the injection and production wells, the²²²Rn content of the return fluid is dependent upon the distribution of flow path lengths and fracture apertures in the reservoir. Matching of the calculated and experimental²²²Rn contents of the return fluids has been used to select appropriate parameters for the reservoir model and hence to estimate the extent of the heat-transfer surface. The model estimates the fracture width of the flow paths, total swept surface area and fracture volume within the reservoir.     

Controls on pore-fluid concentration of 4He and 222Rn and the calculation of 4He/222Rn ages

Controls governing the production of ⁴He and ²²²Rn in the solid phase as well as parameters and processes contributing to their transport into the pore fluid are discussed. ²²²Rn activity is used to quantify the uranium sources for ⁴He and the result is a simplified pore-fluid age equation which is virtually independent of the porosity, the water/rock ratio, the rock density and the uranium concentration, but does require a Th/U estimate. The crucial parameter is the ratio of the release factors for the two species, /GL_Rn//GL_He, which is discussed in terms of three possible release mechanisms: (1) recoil; (2) recoil followed by diffusion; and (3) weathering release of accumulated ⁴He. It was found that /GL_Rn//GL_He can vary over several orders of magnitude, but can be expressed in terms of the effective grain size r, and the diffusive half-length for ²²²Rn decay, r_e.⁴He measurements are used to “date” gases from known gas fields and the agreement with the assumed source-rock age is good. Application of ⁴He/ ²²²Rn measurements to continental freshwater springs indicates that the weathering release of accumulated ⁴He dominates the input and results in a large overestimate of groundwater age. Measurement in the Lardarello geothermal field indicates that the ⁴He/²²²Rn method can indicate relative transport direction. Other possible applications in various geochemical fields are suggested.     

Controls on pore-fluid concentration of He and Rn and the calculation of He/Rn ages

Controls governing the production of ⁴He and ²²²Rn in the solid phase as well as parameters and processes contributing to their transport into the pore fluid are discussed. ²²²Rn activity is used to quantify the uranium sources for ⁴He and the result is a simplified pore-fluid age equation which is virtually independent of the porosity, the water/rock ratio, the rock density and the uranium concentration, but does require a Th/U estimate. The crucial parameter is the ratio of the release factors for the two species, /GL_Rn//GL_He, which is discussed in terms of three possible release mechanisms: (1) recoil; (2) recoil followed by diffusion; and (3) weathering release of accumulated ⁴He. It was found that /GL_Rn//GL_He can vary over several orders of magnitude, but can be expressed in terms of the effective grain size r, and the diffusive half-length for ²²²Rn decay, r_e.⁴He measurements are used to “date” gases from known gas fields and the agreement with the assumed source-rock age is good. Application of ⁴He/ ²²²Rn measurements to continental freshwater springs indicates that the weathering release of accumulated ⁴He dominates the input and results in a large overestimate of groundwater age. Measurement in the Lardarello geothermal field indicates that the ⁴He/²²²Rn method can indicate relative transport direction. Other possible applications in various geochemical fields are suggested.     

Radon modelling and the heat transfer surface area of the British hot dry rock geothermal reservoir

A simple thermal model is developed to evaluate the heat exchange surface area of a hot dry rock geothermal reservoir. This model, in conjunction with the Rn model of Andrews and coworkers, is applied to RH12/RH15 system of the British HDR reservoir. Results suggest that although the estimated Rn transfer surface area represents the actual swept surface area, it is between 25 and 45 times larger than the heat transfer surface area. The difference is explained as due to clustering of several fractures within a range of thermal interaction, over the duration of the circulation. It is also shown that the decline in the measured temperatures of circulation fluids since 1985 to the present is consistent with the heat exchange surface area going through a maximum in 1987.     

Helium, radon and radiocarbon studies on a regional aquifer system of the North Gujarat-Cambay region, India

The study reports the age evolution of groundwater as it flows from the recharge area through a regional alluvial aquifer system in North Gujarat-Cambay region in western India. Radiocarbon (¹⁴C), ⁴He and ⁴He / ²²²Rn dating methods have been employed. Sediments from a drill core in the Cambay Basin were also analysed for uranium (U) and thorium (Th) concentrations and the measured values have been used to estimate the ⁴He and ²²²Rn production rate for groundwater age calculations. Additionally, factors controlling the distribution of ²²²Rn, ⁴He and temperature anomalies in groundwater, vis-à-vis their relation to the tectonic framework and lithology of the study area, have also been examined.The multi-isotope study indicated a reasonable correspondence in groundwater age estimates by the three methods employed. The groundwater ¹⁴C ages increased, progressively, in the groundwater flow direction: from the foothills of Aravalli Mountains in the east, and reached a value of ∼35 ka towards the region of lowest elevation, linking Little Rann of Kachchh (LRK)-Nalsarovar (NS)-Gulf of Khambhat (GK) in the western part of the study area. In this region, groundwater ages obtained for free flowing thermal wells and springs employing ⁴He and ⁴He / ²²²Rn systematics are in the order of million years. Such anomalous ages are possibly due to enhanced mobilisation and migration of ‘excess helium’ from hydrothermal circulation vents along deep-seated faults. Excluding such anomalous cases and considering all uncertainties, presently estimated ⁴He and ⁴He / ²²²Rn groundwater ages are in reasonable agreement with ¹⁴C age estimates in the Cambay Basin for helium release factor (Λ_He) value of 0.4 ± 0.3. The ⁴He method also indicated west-southwards progression of groundwater ages up to ∼100 ka beyond the Cambay Basin.Large ‘excess helium’ concentrations are also seen to be generally associated with anomalous groundwater temperatures (> 35 °C) and found to overlie some of the basement faults in the study area, particularly along the east and the west flanks of the Cambay Basin. Groundwater ²²²Rn activities in most of the study area are 800 ± 400 dpm/l. But, a thermal spring at Tuwa on the east flank of the Cambay Basin, having granitic basement at shallow depth, recorded the highest ²²²Rn activity (∼63,000 dpm/l).     

Hydrogeological and geochemical control of the variations of Rn concentrations in a hard rock aquifer: Insights into the possible role of fracture-matrix exchanges

To investigate the possible variations of Rn concentration in crystalline rocks as a function of flow conditions, a field study was carried out of a fractured aquifer in granite. The method is based on the in situ measurement of Rn in groundwater, aquifer tests for the determination of hydraulic characteristics of the aquifer and laboratory measurement of Rn exhalation rate from rocks. A simple crack model that simulates the Rn concentration in waters circulating in a fracture intersecting a borehole was also tested. The Rn concentrations in groundwaters from boreholes of the study site ranged from 192 to 1597 Bq L⁻¹. The Rn exhalation rates of selected samples of granite and micaschist were determined from laboratory experiments. The results yielded fluxes varying from 0.5 to 1.3 mBq m⁻² s⁻¹ in granite and from 0.5 to 0.9 mBq m⁻² s⁻¹ in micaschists. Pumping tests were performed in the studied boreholes to estimate the transmissivity and calculate the equivalent hydraulic aperture of the fractures. Transmissivities ranged from 10⁻⁵ to 10⁻³ m² s⁻¹. Using the cubic law, hydraulic equivalent fracture apertures were calculated to be in the range of 0.5–2.3 mm.     

Radioisotope dynamics - the origin and fate of nuclides in groundwater

The use of radionuclides as clocks for groundwater dating and as probes to investigate the geometry and spatial extent of the contact area between rocks and water is reviewed. Subsurface production rates for²²²Rn³⁷Ar,⁸⁵Kr,³⁹Ar,³⁶CI,³He,⁴He and⁴⁰Ar in various rock types are listed. Measured Rn fluxes from the surface of sandstone grains and from pieces of granite point to scale-dependent diffusion coefficients. The temporal evolution of subsurface-produced²²²Rn-,³⁷Ar-,⁸⁵Kr- and³⁹Ar-activities in groundwaters yields radionuclide escape factors between 0.1 % and 9% for the Stripa granite (Sweden) and between 1% and 4% for the Milk River sandstone (Canada). The combination of³H,⁸⁵Kr,³⁹Ar,¹⁴C,³⁶CI,⁴He in the UK Triassic sandstone aquifer allows groundwater dating up to 40 000 a. Very old groundwaters can be studied using Cl,³⁶Cl and⁴He evolution as demonstrated in the Milk River aquifer in Canada.     

Granite is an Effective Helium Source Rock: Insights from the Helium Generation and Release Characteristics in Granites from the North Qinling Orogen,China

Global helium(He) shortage is a challenging problem; however, the types of helium source rock and the mechanisms of He generation and release therein remain still poorly understood. In this study, in order to evaluate the potential of granite as an effective helium source rock, we collected granitic samples from the North Qinling Orogen, Central China, in the south of the helium-rich Weihe Basin. The helium generation and release behaviors in granite were studied through analysis of U and Th concentrations, EMPA images, and He and Ar concentrations and isotopic ratios extracted by crushing and stepwise heating. The results indicate that Ar has a better retention and a lower mobility than He. 3 He/4 He ratios released by crushing and stepwise heating are 0.016–0.056 RA and 0.003–0.572 RA, respectively, where RA is the atmospheric 3 He/4 He of 1.4×10-6, reflecting a crustal and radiogenic source. Helium concentrations extracted by the two ways are 0.13–0.95 ucm3 STP/g and 7.82–115.62 ucm3 STP/g, respectively, suggesting that matrix-sited He accounts for more than 98% of total helium preserved in granite. In addition, the total generated He amounts in granites are calculated based on the measured U and Th concentrations in granitic samples. Dividing the preserved He quantities by the generated He amounts, it turned out that less than 10% of He produced since the formation of the granite is preserved in the rock over geological time, suggesting that more than 90% generated He can be transferred to the Weihe Basin. Temperature and fracture are the two critical factors controlling He release. Based on the relationship between He diffusivity of granites and temperature and the He closure temperatures of a variety of U-and Th-rich minerals(27–250°C), we estimate that He can be partially released out of granite at the depths 400 m and totally released at the depths 7800 m. Fractures provide effective transfer of free He from deep source rocks to shallow reservoirs. Finally, a model on granite as an effective helium source rock is established. We suggest exploring He resources in hydrocarbon basins with granitic basement(or adjacent to granite bodies), high geothermal field, and young active fractures.     

The implications of laboratoryRn flux measurements to the radioactivity in groundwaters: the case of a karstic limestone aquifer

Laboratory time-scale experiments were conducted on Carboniferous Limestone gravels from the Mendip Hills area, England, with the purpose of evaluating the release of²²²Rn to the water phase. The specific surface areas of the samples were 4.14 and 1.69 cm² g⁻¹ , which provided, respectively, values of 50.6 and 12.7 pCi for the released Rn. These results allowed the calculation of the emanation coefficient of this rock matrix with respect to the release of Rn, where completely different values corresponding to 23% and 6% were found, suggesting that the extent to which grain boundaries or imperfections in aggregates of micro-crystals of calcite intersect the particle surface certainly affects the Rn release. They also permitted the evaluation of models for the generation of Rn in rocks and transfer to water, in order to interpret the radioactivity due to this gas in groundwaters from the karstic aquifer of the Mendip Hills area, where the calculated activities in groundwater based on the values of 23% and 6% for the emanation coefficient were about 51 and 15 times higher than actually measured in groundwater. Therefore, the emanation coefficient in nature is considerably smaller than in the lab experiment, and another factork (0 < k < 1) may be introduced into the equations related to the modelling, with the aim of adjusting the theoretical-practical results.     

Time lapse tracer and SP measurements to characterize the hydrodynamics of fractured granite aquifer: A case study

An attempt is made to identify and delineate the groundwater flow direction and rate in a fractured hard rock aquifer in Maheshwaram granite watershed in Andhra Pradesh using multiple tracers such as bromide, iodide and rhodamine-B under both natural and induced conditions. A main well in the center and three observation wells at 25 m in north-northwest, southeast and southwest directions respectively were constructed based on ²²²Rn anomalies and ⁴He measurements. Streaming potential survey was also carried out in the area before, during and after pumping test. The tracer tests indicated preferential groundwater flow in N-S direction. The groundwater flow velocity estimated based on arrival of tracer under natural and induced conditions are 0.52 m/d and 375 m/d respectively.     

Exploration of geothermal activity using time series analysis of subsurface gases data from Bakreswar hot springs area,eastern India

Geothermal power seems to be a potential source of green energy in India. But these renewable energy resources are still ignored in India even after having a lot of potential sources as seen in more than 300 hot springs scattered throughout different geothermal areas of the country. Many of them could be utilized for power generation using the earth’s internal heat. More hours are needed to explore these geothermal areas using geochemical, geophysical techniques, and statistical analysis to qualitatively estimate power harnessing capabilities and sustainability of the areas for generation of geothermal power. In the present paper, attention has been focused to investigate Bakreswar geothermal field of India by continuous (24*7) and online monitoring of terrestrial gases such as He and radioactive gases (²²²Rn) in hot spring emanations of the geothermal area for more than 5 years along with discrete measurement of some other geochemical and geophysical parameters. The discrete measures imply that the seven hot spring vents at the study area are linked to the aquifer through different ways covering non-uniform rock assemblies and ascertain the presence of high amount of radioactive minerals at the underneath terrains. Stable activities with a high-out flux of ²²²Rn and He for a prolonged time period is expected within the reservoir present at the study area as consequence of analysis of the temporal variations and statistical measures of the continuous data sets. Seasonal variation of time series data also recommends that the high amount of radioactive sources present at the crust of the reservoir is able to produce enough quantity of heat irrespective of meteorological effects. The investigation on the power spectra interferes that the geothermal system is still in quite active phase on the influence of tectonic activities. Therefore, the geothermal reservoir present at the Bakreswar geothermal area may be utilized as a constant and continuous heat source for a long time period to run a geothermal power plant.     

An experimental study of lake water-sediment interaction rates

Dissolution rates of sediments obtained from the Oued Cherf reservoir were measured in closed-system batch reactors at 25 °C in fluids sampled concurrently from the same locations as the sediments. The BET surface areas of the sediments ranged from 16 to 45 m²/g and consisted primarily of quartz, calcite, and clay minerals. After a brief initial period, release rates of Si, Mg, Ca, Cl, SO₄, and NO₃ from these sediments are approximately linear with time over the course of the experiments, which lasted from 3 to 5 months. BET surface area normalized Si release rates ranged from 10^–17.4 to 10^–18.4 mol/cm²/s. These release rates match closely Si release rates from quartz and clay minerals determined from laboratory dissolution rates reported in the literature. This coherence suggests that laboratory measured silicate dissolution rates can be used with confidence to predict the dissolution behavior of sediments in natural surface waters.     

Fluid–rock interaction in the Qitianling granite and associated tin deposits, South China: Evidence from boron and oxygen isotopes

The behavior of boron and its isotopes in fluid-dominant processes and hydrothermal alteration of granites is examined using fresh and altered granite samples from the Qitianling granite and associated hydrothermal tin deposits in South China. Boron concentrations are highest in the fresh granite (37 ppm) and depleted as a result of two stages of fluid mobility and fluid–rock interaction within the granite. Constraints provided both by δ¹¹B and δ¹⁸O data suggest that the first stage was related to exsolution of aqueous fluids from the granite magma at a temperature of > 450 °C. This was followed by further boron depletion in the granite by hydrothermal circulation of meteoric water at lower temperatures (~ 350 °C) and low water/rock ratios. The sensitivity of coupled boron and oxygen isotope systematics to these processes suggests that they can provide valuable constraints on of fluid mobility in granite and associated mineralization.     

Radioelements,radiogenic helium and age relationships for groundwaters from the granites at Stripa,Sweden

The solution of radioelements and radiogenic ⁴He by groundwaters in fractured rocks is dependent upon the radioelement distribution in the rock matrix and the extent of the rock-water interface. The ${}^{234}\text{U}{}^{238}\text{U}$ activity ratio and the dissolved U, Rn and He contents of such groundwaters respond to changes in the flow regime with time. Although ${}^{234}\text{U}{}^{238}\text{U}$ activity ratios change with groundwater residence time as a consequence of ²³⁴Th-recoil induced solution of ²³⁴U, the activity ratio is strongly influenced by the U distribution within fractures, by the extent of the rock-water interface and by the amount of ²³⁸U in solution. A model for the quantitative evaluation of these effects is presented.Groundwaters from depths up to 880 m in the Stripa granite have variable dissolved uranium contents and ${}^{234}\text{U}{}^{238}\text{U}$ activity ratios. The uranium geochemistry is primarily determined by variations in flow path rather than by groundwater age.Dissolved radiogenic ⁴He in the groundwaters increases with their depth of origin, and is dependent upon the U content of the granite and upon its fracture porosity. It increases with groundwater residence time but movement of ⁴He by diffusion and transport processes make the actual groundwater age indeterminate.     

Uranium, radon, helium and other trace elements and gases in well waters of parts of the St. Lawrence Lowlands, (Ottawa region) Canada

To test the usefulness of groundwater for U prospecting, 130 domestic wells were sampled from a 22,000 km² rectangular area extending south and west from Ottawa, Ontario. The waters were analyzed for twenty variables including ten gases and five trace elements.The dissolved gases give information on subsurface chemical and geological environments. Carbon dioxide highs occur at the contact of granites and limestones; He highs are along major faults; CH₄ and H₂S reveal strongly reducing environments, and O₂ in wells indicates shallow wells or active water turnover. Rn reveals U mineralization and granitic rocks.The South March U-Cu surface occurrence gives coincident U, Cu, Zn, Pb and Rn highs. A weak U-Rn-Cu high in the Richmond-North Gower area suggests more radioactive subcropping March formation there.A strong regional He anomaly coincident with the outcropping Rockcliffe Formation and controlled by the Hazeldean Fault is believed to be due to deep-seated He, but its coincidence with a weak regional Cu anomaly, a Rn anomaly in the Bells Corners area, its linear northeast trend, the presence of Helikian U-rich source rocks to the southwest, and northeast-dipping Paleozoic sandstones provide both indications of, and the right environment for epigenetic type U-Cu occurrences.     

Water–granite interaction: Clues from strontium,neodymium and rare earth elements in soil and waters

Strontium-, Nd-, and rare-earth-element-isotope data are presented from rock, weathered rock (arene) and saprolite, sediment and soil, shallow and deep groundwater (e.g. mineral-water springs), and surface waters in the Margeride massif, located in the French Massif Central. Granitoid rock and gneiss are the main lithologies encountered in the Margeride, which corresponds to a large and 5-km-deep laccolith. Compared to bedrock, the Sr isotopes in arene, regolith, sediment and soil strongly diverge with a linear increase in the ⁸⁷Sr/⁸⁶Sr and Rb/Sr ratios. Neodymium isotopes fluctuate least between bedrock and the weathering products. In order to characterise the theoretical Sr isotopic signature IRf(Sr) of water interacting with granite, a dissolution model was applied, based on the hypothesis that most of the Sr comes from the dissolution of plagioclase, K-feldspar and biotite. Similar to the Sr model, an approach was developed for modelling the theoretical Nd isotopic signature IRf(Nd) of water interacting with a granite, assuming that most Nd originates from dissolution of the same minerals as those that yield Sr, plus apatite. The IRf(Sr) ratio of water after equilibration with the Sr derived from minerals was calculated for the Margeride granite and compared to values measured in surface- and groundwaters. Comparison of the results shows agreement between the calculated IRf(Sr) and the observed ⁸⁷Sr/⁸⁶Sr ratios. When calculating the IRf(Nd) ratio of water after equilibration with the Nd derived from minerals of the Margeride granite, the results indicated good agreement with surface-water values, whereas mineralised waters analysed within the Margeride hydrosystem could not be directly linked to weathering of the granite alone. Because the recharge area of deep groundwater is located on the Margeride massif, very deep circulation involving interaction with other rocks (e.g. shales) at depths of >5 km must be considered.     

Soil gas prospection of He,222Rn and CO2: Vulcano Porto area,Aeolian Islands,Italy

In March 1994, soil gases were sampled in the area of Vulcano Porto, on the island of Vulcano, using a grid of about 200 points/km². Analysed gases were CO₂, He and²²²Rn and, over a smaller area, H₂S. Some of the samples were also analysed for the isotope composition of CO₂C. Three anomalous CO₂ degassing areas were identified: Grotta dei Palizzi, the area near the Telephone Exchange, and the area near the beach fumaroles. The behaviour of He and²²²Rn is different in these 3 areas. The concentration of He is much lower than that of atmospheric He (down to −3950 ppb) in the isthmus, and only in the area near Grotta dei Palizzi does it have values significantly higher than atmospheric ones (up to 1900 ppb). The activity of²²²Rn, always significantly positively related to CO₂ concentrations, peaked in the 2 fumarole areas (isthmus, Telephone Exchange).     

Main characteristics of the deep geothermal brine (5 km) at Soultz-sous-Forêts (France) determined using geochemical and tracer test data

Three deep wells (5000 m) have been drilled into a fractured granite basement at Soultz-sous-Forêts, within the Tertiary Rhine Graben, in order to develop a heat exchanger and produce electricity after the creation of an EGS reservoir. Very few analyses representative of the deep geothermal fluids are available because of frequent contamination by drilling fluids or injected waters. These indicate similar chemical and isotopic compositions (NaCl fluids) and high salinities (about 100 g/l) suggesting a common sedimentary origin and identical water–rock interaction processes at equilibrium temperatures close to 230 °C in a sedimentary rather than a granite reservoir. The latter would be situated closer to the Graben centre where the Triassic Buntsandstein formation is deepest and hottest. Tracer tests conducted after 2000 show that the deep native geothermal brine is omnipresent in the fluids discharged during the production and circulation tests. Its natural convective flux was estimated at 1–1.2 m³/h.     

Mechanism of transport of U-Th series radioisotopes from solids into ground water

A set of experiments conducted to understand the mechanism responsible for release of large quantities of ²²²Rn from solids into ground water reveals that a major part of Rn moving into the intergranular water comes from within the grains (solids), not by recoil from the outer surface. We have deduced that the solids we have studied are permeated, to a varying degree, with pores having very large wall area but very small volume, i.e., having width of opening of the order of 10–20 nm. Recoil from the walls of these openings introduces Rn and other isotopes into the water contained in these openings (nanopores). Radon, an inert gas, is able to diffuse from nanopore water into the intergranular water in substantial quantities. In contrast, nongaseous isotopes of the U-Th series recoiled into the nanopore water in comparable quantities are adsorbed within the nanopores on the large internal area available there. Further, it appears that a large part of the inventory of long-lived isotopes produced by decay of radon in intergranular water is slowly transported to and adsorbed on the large internal area available in the nanopores. Thus, most of the inventory of all long-lived non-gaseous emanation isotopes exist in an adsorbed state in the nanopores, while radon is located in intergranular water.