Te Kopia geothermal system (New Zealand) — the relationship between its structure and extent

The Te Kopia geothermal system is one of several high-temperature systems in the Taupo Volcanic Zone (TVZ) of New Zealand. It is located along the Paeroa Fault Zone, a major active fault system trending NE in the central TVZ. Three independent studies, i.e. resistivity survey, magnetic interpretation, and detailed topographic analysis of faults and fractures, indicate the existence of another fault system, trending NW, that also significantly influences the Te Kopia geothermal system. Results from these studies also show that, at Te Kopia, a resistivity low and hydrothermally demagnetised rocks (both are indicators of a geothermal reservoir in volcanic rocks) clearly coincide with a zone of high fault and fracture density. Hence, the Te Kopia field is a good example of the significant influence that geological structures (major fault systems) have on the extent of a geothermal reservoir, by creating zones of fractured rocks that provide permeable paths for thermal fluids.     

Spatial data analysis and integration for regional-scale geothermal potential mapping,West Java,Indonesia

Conceptual modeling and predictive mapping of potential for geothermal resources at the regional-scale in West Java are supported by analysis of the spatial distribution of geothermal prospects and thermal springs, and their spatial associations with geologic features derived from publicly available regional-scale spatial data sets. Fry analysis shows that geothermal occurrences have regional-scale spatial distributions that are related to Quaternary volcanic centers and shallow earthquake epicenters. Spatial frequency distribution analysis shows that geothermal occurrences have strong positive spatial associations with Quaternary volcanic centers, Quaternary volcanic rocks, quasi-gravity lows, and NE-, NNW-, WNW-trending faults. These geological features, with their strong positive spatial associations with geothermal occurrences, constitute spatial recognition criteria of regional-scale geothermal potential in a study area. Application of data-driven evidential belief functions in GIS-based predictive mapping of regional-scale geothermal potential resulted in delineation of high potential zones occupying 25% of West Java, which is a substantial reduction of the search area for further exploration of geothermal resources. The predicted high potential zones delineate about 53–58% of the training geothermal areas and 94% of the validated geothermal occurrences. The results of this study demonstrate the value of regional-scale geothermal potential mapping in: (a) data-poor situations, such as West Java, and (b) regions with geotectonic environments similar to the study area.     

A helium isotope perspective on the Dixie Valley,Nevada, hydrothermal system

Fluids from springs, fumaroles, and wells throughout Dixie Valley, NV were analyzed for noble gas abundances and isotopic compositions. The helium isotopic compositions of fluids produced from the Dixie Valley geothermal field range from 0.70 to 0.76 Ra, are among the highest values in the valley, and indicate that ∼7.5% of the total helium is derived from the mantle. A lack of recent volcanics or other potential sources requires flow of mantle-derived helium up along the valley bounding Stillwater Range Front Fault, from which the geothermal fluids are produced. Using a one-dimensional flow model, a lower limit fluid flow rate up through the fault of 7 mm/yr is estimated, corresponding to a mantle ³He flux of ∼10⁴ atoms m⁻² s⁻¹.A comparison between the fluids from Dixie Valley springs, fumaroles, and wells and the fluids produced from the geothermal field reveals a mixing trend between the geothermal fluid and younger, cooler groundwaters. The exceptions are those features that either emanate directly from the Stillwater fault or wells that penetrate and extract fluids from the fault zone, all of which have helium isotopic compositions that are indistinguishable from the geothermal production fluids. The results of our study indicate that the Stillwater Range Front Fault system must act as a permeable conduit that can sustain high vertical fluid flow rates from deep within the crust and crust-mantle boundary and that high permeability may exist along most of its length. This suggests that the geothermal potential of the Stillwater fault may be significantly greater than the 6–8 km long system presently under production. Since all the numerous springs, wells, and fumaroles in the valley also contain a fluid component that is indistinguishable from the geothermal/Stillwater fault fluid, the potential for an additional deeper and more pervasive geothermal system also exists and should be further evaluated. Furthermore, we suggest that elevated helium isotope compositions in regions with little or no recent magmatism are an indicator of the deep crustal permeability that is required to drive and sustain extensional geothermal systems.     

Basement structure, lithology and permeability at Kawerau and Ohaaki geothermal fields, New Zealand

Poorly permeable basement rocks commonly occur in geothermal regions around the world, and the Quaternary Taupo Volcanic Zone (TVZ) of New Zealand is no exception. Production from basement terrane requires detailed knowledge of its geological and geophysical parameters, as shown by the history of Kawerau and Ohaaki, the only geothermal fields in the TVZ where Mesozoic Torlesse terrane greywacke (litharenite) basement is commonly penetrated at drilled depths of 1–2.5 km. In both fields the basement is step-faulted down into the TVZ. Although hot and hydrothermally altered, the greywackes have little permeability. Some production wells feed from elusive basement faults at Kawerau, but rarely at Ohaaki. Greywackes at Ohaaki are of “granite-rhyolite” provenance, and have more interbedded argillite than the “andesite-dacite” derived Kawerau greywackes. In consequence, the Kawerau basement may sustain brittle fracture at higher temperatures and depths than the more ductile Ohaaki basement, allowing convective circulation of higher enthalpy fluids into permeable Quaternary aquifers.     

Geophysical investigations of the Seferihisar geothermal area, Western Anatolia, Turkey

Self-potential (SP), magnetic and very low frequency electromagnetic (EM-VLF) surveys were carried out in the Seferihisar geothermal area to identify major and minor fault zones and characterize the geothermal system. The SP study provided useful information on the local faults and subsurface fluid flow. The main SP anomalies appear mostly along and near active fault zones in the area of the Cumalı, Tuzla and Doğanbey hot springs. Two of these anomalies near the Tuzla hot springs were further evaluated by SP modelling. Total magnetic field values increase from the Doğanbey to the Cumalı hot springs. Modelling performed on the magnetic data indicates that between these two spring areas are four different regions or units that can be distinguished on the basis of their magnetic susceptibility values. Fraser filtering of EM-VLF data also indicates that there are three significant conductive zones in the regions around the Cumalı, Tuzla and Doğanbey hot springs, and that they lie between important fault systems. The EM-VLF and total (stacked) SP data show that the conductive tilt anomalies obtained by Fraser filtering generally coincide with negative SP areas.According to our geophysical investigations, new exploratory wells should be drilled into the conductive zones located between the Cumalı and Tuzla hot springs. We further recommend that resistivity and magnetotelluric methods be carried out in the area to obtain additional information on the Seferihisar geothermal system.     

Geologic structure and volcanic history of the Yanaizu-Nishiyama (Okuaizu) geothermal field, Northeast Japan

The Yanaizu-Nishiyama geothermal field, also known as Okuaizu, supports a 65 MWe geothermal power station. It is located in the western part of Fukushima Prefecture, northeast Japan. This field is characterised by rhyolitic volcanism of about 0.3–0.2 Ma that formed Sunagohara volcano. Drillcore geology indicates that volcanism began with a caldera-forming eruption in the center of this field, creating a 2-km-diameter funnel-shaped caldera. Subsequently, a fault-bounded block including this caldera subsided to form a 5-km-wide lake that accumulated lake sediments. Post-caldera volcanism formed lava domes and intrusions within the lake, and deposited ash-flow tuffs in and around the lake. The hydrothermal system of this field is strongly controlled by subvertical faults that have no relation to the volcanism. The principal production zone occurs at a depth of 1.0–2.6 km within fractured Neogene formations along two northwest-trending faults to the southeast of the caldera. These faults also formed fracture zones in the lake sediments, but there was no apparent offset of the sediments. Stratigraphic studies suggest that post-caldera activities of Sunagohara volcano have migrated southeastward to the present high-temperature zone. The source magma of Sunagohara volcano may contribute to the thermal potential of this field.     

Intersecting faults and sandstone stratigraphy at the cerro prieto geothermal field

The NW - SE trending Cerro Prieto fault zone is part of a major regional lineament that extends into Sonora, Mexico, and has characteristics of both a wrench fault and an oceanic transform fault. The zone includes a number of separate identified faults and apparently penetrates deep into the basement and crustal rocks in the area. The zone serves as a conduit for both large and rapid heat flow. Near well M-103, where the Michoacán fault obliquely intersects a shorter NE - SW trending fault (i.e., the Pátzcuaro fault), large circulation losses during drilling indicate greater permeability and hence increased natural convective fluid flow. Temperature contour maps for the southern portion of the field suggest that a shear fault zone also exists in the vicinity of wells M-48, M-91 and M-101. This shear zone aids in rapidly distributing geothermal fluid away from the Cerro Prieto fault zone, thus enhancing recharge to the western part of the reservoir.We have studied the distribution of lithologies and temperature within the field by comparing data from well cuttings, cores, well logs and geochemical analyses. Across the earliest developed portion of the field, in particular along a 1.25 km NE - SW section from well M-9 to M-10, interesting correlations emerge that indicate a relationship among lithology, microfracturing and temperature distribution. In the upper portion of the reservoir of this section, between 1200 and 1400 m, the percentage of sandstones ranges from 20 to 55. Well logs, calcite isotope maxima, and the Na - K - Ca geothermometer indicate temperatures of 225–275°C. The isothermal high in this vicinity corresponds to the lowest total percentage of sandstones. Scanning electron microphotographs of well cores and cuttings from sandstone and shale units reveal open microfractures, mineral dissolution and mineral precipitation along microfractures and in pores between sand grains. Our working hypothesis is that these sandy shale and siltstone facies are most amenable to increased microfracturing and, in turn, such microfracturing allows for higher temperature fluid to rise to shallower depths in this part of the reservoir.Our ongoing research is aimed at achieving a coherent geological model that provides a basis for estimating reservoir capacity, and that illustrates our understanding of fluid flow along major faults, laterally through fault shear zones, and within predominantly silty and shaley deltaic clastics that have been microfractured.     

Numerical modeling of Basin and Range geothermal systems

Basic qualitative relationships for extensional geothermal systems that include structure, heat input, and permeability distribution have been established using numerical models. Extensional geothermal systems, as described in this paper, rely on deep circulation of groundwater rather than on cooling igneous bodies for heat, and rely on extensional fracture systems to provide permeable upflow paths. A series of steady-state, two-dimensional simulation models is used to evaluate the effect of permeability and structural variations on an idealized, generic Basin and Range geothermal system of the western U.S.Extensional geothermal systems can only exist in a relatively narrow range of basement (bulk) permeability (10⁻¹⁵ m² to 10⁻¹⁶ m²). Outside of this window, shallow subsurface fault zone temperatures decrease rapidly. Mineral self-sealing does not significantly affect the flow system until the flow path is almost completely sealed off. While topography gives an extra “kick” to convective circulation, it is not a requirement for geothermal system development. Flow from the ranges to the fault dominates the circulation, while secondary flow systems exist on the range front slopes. A permeable fault in one valley can also induce cross-range flow if there are no equally good upflow paths in the adjacent valleys. When bulk permeability is high enough, additional deep circulation cells develop in adjacent valleys, diverting heat and fluid from the fault and consequently reducing temperatures in the fault itself. Qualitative comparison between temperature–depth logs from actual geothermal systems and from the generic models is a significant aid to understanding real-world geothermal fluid flow, and suggests new or better interpretations of existing systems.     

Application of chemical geothermometers to thermal springs of the Maghreb,North Africa

In order to assess the geothermal potential in the Maghrebian region, several studies have been undertaken in the three countries concerned, Morocco, Algeria and Tunisia, during the past decade. Research programmes have considered the surface evidence (thermal springs) and underground thermal information from deep and shallow wells. The main chemical characteristics of the sampled thermal springs and the results of the application of geothermometers as result from these studies are presented. Of the 238 inventoried thermal springs, 169 have been selected, on the basis of complete water analyses and acceptable ionic balances. Measured temperatures range from 22.5 to 98°C, thermal indexes from 0.5 to 78°C and salinities from 0.13 to 52.5 g/L. Most studied springs are sodium-chloride type waters. These basic data allow identification of the main thermal anomalies in the Maghrebian zone, which are located in regions of the Libyan-Tunisian, Algerian-Moroccan and Algerian-Tunisian frontiers, of northern Tunisia, the Eastern Rif and the northern part of the Saharan Atlas.Several chemical geothermometers have been applied to selected springs: NaK, NaKCa, NaKCaMg, Na/Li, Mg/Li, K²/Mg, quarts, chalcedony (Fournier) and chalcedony (Arnorsson). The NaK, NaKCa, NaKCaMg, Na/Li and Mg/Li geothermometers seem to give unreliable results, while K²/Mg and silica temperatures are apparently reasonable. However, dissolved silica seems to be governed by quartz solubility for some thermal springs and by chalcedony solubility for others. The results are tentatively compared with known geothermal gradients and geological features.     

The nestos delta geothermal field and its relation to the associated self potential (sp) field

Prospecting for geothermal reservoirs by geophysical methods has proved to be a challenge in recent years. In the case of Nestos geothermal field, considered to be a blind field (no surface manifestation), the geological and geophysical data were studied and intercorrelated. The geophysical results from SP, gravity, and VES data were compared with basement morphology and tectonics, as well as with the high temperatures measured in the area of main geothermal interest. As a result, the highly conductive subsurface zone and high temperatures observed on the ground surface were closely associated with a specific fracturing system. The latter was successfully mapped by the SP method. The SP method also defined certain fracture zones which, being highly electrically polarized, are hypothesized as future geothermal targets.All existing geophysical data have been re-evaluated, along with the SP data acquired over the geothermal field, and compared with the thermal contour maps and regional tectonics of the area.A deep borehole, drilled earlier in the area of the geothermal field, detected high temperatures in the basement (115°C), but very low flowrate of the geothermal fluids. According to the results of this study, this is due to the fact that the borehole was located outside the area of the main fracture zones of geothermal interest.     

Geology of the tengchong geothermal field and surrounding area, West Yunnan, China

This paper describes the regional tectonics, active structures and thermal springs of the Tengchong geothermal field and surrounding area. Regional tectonic evidence shows that Tengchong and the surrounding area was a microcontinent between the Gondwanaland and Eurasian plates. The distribution of the thermal springs is controlled by active faults in the “arched” and circular structures. The majority of the thermal areas have a deep circulation within the fault system. Only a few high-temperature hydrothermal systems issuing sodium chloride may be heated by a local magmatic source.     

The use of soil Hg to delineate zones of upwelling in low-to-moderate temperature geothermal systems

A soil mercury survey was conducted near the town of Calistoga, California to identify and delineate a buried fault system that is thought to control the upwelling of low-to-moderate temperature geothermal fluids in the upper Napa Valley. Soil samples were collected at 100 m intervals along traverses that crossed hot springs and existing geothermal well sites. Strong mercury anomalies occur along a broadly-defined zone and in close proximity to surface thermal activity including active hot springs and silica sinter formations. In contrast, background mercury concentrations are present in locations with little or no indication of subsurface thermal activity, such as along the margins of the valley or near groundwater wells producing non-thermal water. Analysis of smoothed Hg values reveals a N65W-trending lineament of high Hg concentrations. These results suggest that soil mercury surveys can be a useful and cost-effective method for the identification and mapping of structures controlling subsurface fluid flow in low-to-moderate temperature geothermal systems.     

Geothermal resources in the shallow, unsaturated zone of the Wiesbaden spa district, Germany

The geothermal energy potential of the Wiesbaden spa district continues to be largely untapped. Although the thermal water is being utilized, any other activity liable to disturb the hydraulic regime of the thermal springs is prohibited. The geothermal potential of the unsaturated zone in the spa district has been mapped, although the zone must be regarded as a thermal insulator rather than a good conductor. The results of the mapping revealed the presence of a number of heat anomalies that could be exploited in the future; further economic benefits could be gained from installing the heat transfer units during road works. The modelling studies considered two possible scenarios: direct heating and heat pump usage. The results indicate that heat pumps are the more efficient option, yielding a thermal capacity of approximately 100 W/m².     

Identification of structures inside basement at Soultz-sous-Foret (France) by triaxial drill-bit vertical seismic profiling

Seismic signals produced while drilling a Hot Dry Rock (HDR) geothermal well in Soultz-sous-Forêt (France) were detected using a 4-component seismic detector and two single-component detectors in the granite basement. The recorded signals had a signal-to-noise ratio of around 8 dB with respect to the background noise and show relatively high spectral whiteness. Results from three-dimensional hodogram analysis were successfully interpreted using a model in which the drill bit generates predominantly continuous shear waves.Reflected waves were detected using both the spatial and temporal cross-correlation of hodograms. The distribution of reflectivity from the tri-axial drill-bit vertical seismic profiling (TAD–VSP) suggests that the bottom of the artificial geothermal reservoir, created by hydraulic stimulation, has a high reflection coefficient. The pre-existing permeable zone inside the reservoir was also detected by the TAD–VSP. The reflection image of the TAD–VSP agrees well with that found by the AE (Acoustic Emission) reflection method. The TAD–VSP was successfully used to detect the internal structure of the deep basement.     

Geothermal resources in Algeria

The geothermal resources in Algeria are of low-enthalpy type. Most of these geothermal resources are located in the northeastern of the country. There are more than 240 thermal springs in Algeria. Three geothermal zones have been delineated according to some geological and thermal considerations: (1) The Tlemcenian dolomites in the northwestern part of Algeria, (2) carbonate formations in the northeastern part of Algeria and (3) the sandstone Albian reservoir in the Sahara (south of Algeria). The northeastern part of Algeria is geothermally very interesting. Two conceptual geothermal models are presented, concerning the northern and southern part of Algeria. Application of gas geothermometry to northeastern Algerian gases suggests that the reservoir temperature is around 198 °C. The quartz geothermometer when applied to thermal springs gave reservoir temperature estimates of about 120 °C. The thermal waters are currently used in balneology and in a few experimental direct uses (greenhouses and space heating). The total heat discharge from the main springs and existing wells is approximately 642 MW. The total installed capacity from producing wells and thermal springs is around 900 MW.     

High-temperature measurements in well WD-1A and the thermal structure of the kakkonda geothermal system, Japan

The New Energy and Industrial Technology Development Organization (NEDO) drilled well WD-1a between 1994 and 1995 in the Kakkonda geothermal field as part of their Deep Seated Geothermal Resources Survey project. High-temperature measurements were carried out in WD-1a. Logging temperatures above 414°C were confirmed at 3600 m and 3690 m depth after 82 h standing time. Simple Horner extrapolations based on observed temperatures up to 82 h after shut-in suggested a temperature of about 500°C at 3500 m depth. Temperatures between 500°C and 510°C were also confirmed at 3720 m depth after 129–159 h standing time, using calibrated melting .tablets. These are the highest temperatures measured in a geothermal well. These results suggest a thermal structure consisting of three layers. Layer one is a shallow permeable zone of the reservoir, at less than 1500 m depth, at 230°C to 260°C. The second layer is a deep zone of the reservoir, which is less permeable and has a temperature of 350°C to 360°C from 1500 m to about 3100 m depth. The third layer is a zone of heat conduction. The transition between the hydrothermal-convection zone and the deeper heat-conduction zone is at 3100 m depth in well WD-1a.     

Monitoring production using surface deformation: the Hijiori test site and the Okuaizu geothermal field,Japan

Production in geothermal reservoirs often leads to observable surface displacement. As shown in this paper, there is a direct relationship between such displacement and reservoir dynamics. This relationship is exploited in order to image fluid flow at two geothermal field sites. At the first locality, the Hijiori Hot Dry Rock (HDR) test site, 17 tilt meters record deformation associated with a 2.2 km deep injection experiment. Images of fluid migration along a ring fracture system of the collapsed Hijiori caldera are obtained. At the Okuaizu geothermal field, leveling and tilt meter data provide constraints on long- and short-term fluid movement within the reservoir. A set of 119 leveling data suggest that the north-to-northeast trending Takiyagawa fault acts as a barrier to flow. The northwesterly oriented Chinoikezawa and Sarukurazawa faults appear to channel fluid from the southeast. The tilt data from Okuaizu indicate that a fault paralleling the Takiyagawa fault zone acts as a conduit to transient flow, on a time scale of several weeks. The volume strain in a region adjacent to the injection wells reaches a maximum and then decreases with time. The transient propagation of fluid along the fault may be due to pressure build-up, resulting from the re-initiation of injection.     

Geochemistry of thermal springs around Lake Abhe,Western Djibouti

The Republic of Djibouti, occupying an area of 23,180 km², falls within the arid zone of East Africa and is located above the ‘Horn of Africa’, adjacent to the Red Sea. This country has several thermal springs and fumaroles distributed over three regions – Lake Assal, Lake Hanle and Lake Abhe. The most characteristic feature of Lake Abhe is the presence of several linear chains of travertine chimneys. The thermal waters are typical of the Na-Cl type near neutral waters rich in CO₂. These waters show an oxygen shift, indicating reservoir temperatures>200°C. The chemical signature of the thermal springs and the geology of the Lake Abhe region are very similar to the Tendaho geothermal area of Ethiopia. The geology, temperature gradient and its proximity to Damah Ale volcano make the Lake Abhe region a potential site for geothermal power development.     

Geothermal areas in Pakistan

In this paper an attempt has been made to correlate the tectonic and geologic features with surface manifestations of geothermal activity in Pakistan to delineate prospective areas for exploration and development of geothermal energy. Underthrusting of the Arabian plate beneath the Eurasian plate has resulted in the formation of Chagai volcanic arc which extends into Iran. Quaternary volcanics in this environment, along with the presence of thermal springs, is an important geotectonic feature revealing the possible existence of geothermal fields. Geothermal activity in the northern areas of Pakistan, as evidenced by thermal springs, is the likely result of collision and underthrusting of the Indian plate beneath the Eurasian plate. Numerous hot springs are found along the Main Mantle thrust and the Main Karakorum thrust in Chilas and Hunza areas respectively. The concentration of hot springs in Sind Province is also indicative of geothermal activity. A string of thermal seepages and springs following the alignment of the Syntaxial Bend in Punjab Province is also noteworthy from the geothermal viewpoint.In Baluchistan Province (southwest Pakistan), Hamun-e-Mushkhel, a graben structure, also shows geothermal prospects on the basis of aeromagnetic studies.     

Hydrogeochemical exploration of geothermal prospects in the Tecuamburro Volcano region, Guatemala

Chemical and isotopic analyses of thermal and nonthermal waters and of gases from springs and fumaroles are used to evaluate the geothermal potential of the Tecuamburro Volcano region, Guatemala. Chemically distinct geothermal surface manifestations generally occur in separate hydrogeologic areas within this 400 km² region: low-pressure fumaroles with temperatures near local boiling occur at 1470 m elevation in a sulfur mine near the summit of Tecuamburro Volcano; non-boiling acid-sulfate hot springs and mud pots are restricted to the Laguna Ixpaco area, about 5 km NNW of the sulfur mine and 350–400 m lower in elevation; steam-heated and thermal-meteoric waters are found on the flanks of Tecuamburro Volcano and several kilometers to the north in the andesitic highland, where the Infernitos fumarole (97°C at 1180 m) is the primary feature; neutral-chloride hot springs discharge along Rio Los Esclavos, principally near Colmenares at 490 m elevation, about 8–10 km SE of Infernitos. Maximum geothermometer temperatures calculated from Colmenares neutral-chloride spring compositions are 180°C, whereas maximum subsurface temperatures based on Laguna Ixpaco gas compositions are 310°C. An exploration core hole drilled to a depth of 808 m about 0.3 km south of Laguna Ixpaco had a bottom-hole temperature of 238°C but did not produce sufficient fluids to confirm or chemically characterize a geothermal reservoir. Hydrogeochemical data combined with regional geologic interpretations indicate that there are probably two hydrothermal-convection systems, which are separated by a major NW-trending structural boundary, the Ixpaco fault. One system with reservoir temperatures near 300°C lies beneath Tecuamburro Volcano and consists of a large vapor zone that feeds steam to the Laguna Ixpaco area, with underlying hot water that flows laterally to feed a small group of warm, chloriderich springs SE of Tecuamburro Volcano. The other system is located beneath the Infernitos area in the andesitic highland and consists of a lower-temperature (150–190°C) reservoir with a large natural discharge that feeds the Colmenares hot springs.