1 Introduction
The development of high temperature geothermal energy requires a minimum temperature of 200 °C. In Europe, these conditions are realized at limited depth for exceptional high geothermal gradient areas (Iceland or the Larderello geothermal field in Italy), but generally, more or less 5 km depth is required (Gérard et al., 1991). Fluid flow through a geothermal exchanger within a deep reservoir, such as a granitic batholith (a low porosity and permeability material), depends on development of secondary structures like fault and fracture networks, with a good connectivity and permeability. These properties could be estimated by productivity tests and improved by mechanical and/or chemical stimulation processes.
Only few faults or fractures intercepting the boreholes are useable to develop the connected network (Evans et al., 2005), and fluid transfers depend on their aperture, which could control a rapid transfer pathway (Gentier et al., 2000). However, the geometry of a fault zone is more complex than a simple parallel plate model defining a homogeneous space between fault planes: it is composed of a fault core, a damage zone and the protolith (Caine et al., 1996; Sibson and Rowland, 2003). This geometry could control a large set of fluid-rock interactions as well as chemical and physical processes. The aim of this article is the quantification of different physical properties of fault zone materials that was qualitatively described by Genter et al., (2000).
2 EPS-1 borehole and sampling
Different holes were drilled through the Soultz-sous-Forêts granite to generate a circulation at great depth but only one of them (EPS-1) was cored between 1400 and 2300 m in the granite basement. Seventy-nine samples were selected along this core at different depths (see Table 1 for their locations and their petrophysical properties). Previous to the geothermal stimulation, the Hercynian and Alpine phases and especially the Tertiary extensive episode (Villemin and Bergerat, 1987) structured the granite body. This long geological history is associated with tectonic fracturing, hydrothermal fluid circulations, primary minerals dissolution and secondary minerals precipitation that created not only a network of connected open fractures but also porous and extended permeable zones around them (Genter et al., 1990; Genter and Traineau, 1992; Ledésert et al., 1993; Rosener, 2007; Sausse et al., 1998; Sausse, 2002; Surma, 2003).
Liste des échantillons et propriétés physiques.
Sample | Depth | Facies | sp. surface | Density | TC | Porosity | Permeability |
K74-2016 | 1420 | AHI | 2.31 | ||||
K75-2051 | 1428 | AHH | 2.95 E − 17 | ||||
K75 2051 | 1430.6 | 2.67 | 3.17 | 2.12 | 3.00 E − 17 | ||
K76 2059 | 1431.9 | 2.61 | 3.23 | ||||
K76 2061 | 1432.8 | 2.59 | 2.32 | 3.29 | |||
K79-2128 | 1455 | AHI | 4.33 E + 06 | 2.91 | 1.75 E − 19 | ||
K79-2128 | 1455 | AHI | 5.07 E + 06 | 3.33 | 9.18 E − 21 | ||
K81-2158 | 0.73 | ||||||
K81-2191 | 1471 | AHH à AHCH | 2.9 | 5.59 | 1.93 E − 18 | ||
K81-2191 | 1471 | AHH à AHCH | 2.96 | 5.59 | 1.93 E − 18 | ||
K81-2191 | 1471 | AHH à AHCH | 2.93 | 6.92 | 1.02 E − 18 | ||
K81-2191 | 1471 | AHH à AHCH | 2.96 | 6.92 | 1.02 E − 18 | ||
K88-2350 | 1517 | AHH à AHCH | 2.81 | ||||
K88 | 1517.5 | 2.63 | 3.68 | 3.1 | |||
K97 | 1574.3 | 2.66 | 3.71 | 0.34 | |||
K102 | 1602 | Fresh | 59100 | 0.13 | 5.87 E − 18 | ||
K108 | 1630.5 | 2.64 | 2.88 | 1.03 | 4.30 E − 17 | ||
K108-2758 | 1631 | AHI | 5.88 E + 05 | 1.14 | |||
K108-2758 | 1631 | AHI | 1.84 E + 06 | 1.21 | 4.09 E − 19 | ||
K108-2758 | 1631 | AHI | 3.47 E + 06 | 1.59 | 4.57 E − 18 | ||
K108-2758 | 1631 | AHI | 1.97 | 1.28 E − 18 | |||
K109-2782 | 1637 | AHI | 1.65 E + 06 | 0.27 | 1.82 E − 15 | ||
K109-2782 | 1637 | AHI | 0.61 | 5.29 E − 17 | |||
K109-2789 | 1639.1 | AHH | 2.69 | 3.59 | |||
K109-2793 | 1640.5 | 2.52 | 2.78 | ||||
K109-2794 | 1641 | 13.6 | |||||
K109-2795 | 1641.9 | AHH | 7.41 E + 06 | 3.93 | |||
K109-2795 | 1641.9 | AHH | 1.50 E + 06 | 1.15 | 1.58 E − 15 | ||
K109-2795 | 1641.9 | AHH | 3.52 E + 06 | 0.81 | 2.08 E − 17 | ||
K109-2795 | 1641.9 | 2.6 | 3.12 | 5.12 | 1.30 E − 17 | ||
K117 | 1670.9 | 0.56 | |||||
K117-2900 | 1679.9 | Fresh | 1.70 E + 06 | 2.84 | 0.82 | 1.35 E − 17 | |
K117-2900 | 1679.9 | Fresh | 2.84 | 0.22 | 2.36 E − 17 | ||
K117-2900 | 1679.9 | Fresh | 2.65 | 3.45 | 0.56 | ||
K122 | 1705.9 | 2.67 | 2.92 | 0.33 | |||
K134-3221 | 1776.1 | Fresh | 2.64 | 2.55 | 0.56 | ||
K136 | 1797 | Fresh | 0.44 | ||||
K136 | 1797 | Fresh | 0.38 | 2.14 E − 17 | |||
K136 | 1797 | Fresh | 0.4 | 4.63 E − 17 | |||
K136a | 1797.6 | 2.68 | 3.33 | 0.41 | 1.00 E − 17 | ||
K136b | 1797.7 | 0.41 | |||||
K138-3318 | 1807 | AHI | 2.51 E + 06 | 1.82 | 1.20 E − 18 | ||
K138-3318 | 1807 | AHI | 2.90 E + 05 | 0.5 | 2.90 E − 17 | ||
K138-3318 | 1807 | AHI | 0.8 | 1.32 E − 17 | |||
K142 | 1845.3 | 2.69 | 3.56 | 0.71 | |||
K142 | 1846 | 0.37 | 7.79 E − 18 | ||||
K142 | 1846 | 0.36 | 2.90 E − 20 | ||||
K150 | 1905 | 0.66 | |||||
K150-3645 | 1906 | AHI | 2.00 E + 06 | 0.33 | 8.17 E − 19 | ||
K150-3645 | 1906 | AHI | 0.4 | 5.07 E − 18 | |||
K150-3645 | 1906 | AHI | 0.35 | 2.44 E − 20 | |||
K160-3861 | 1973 | 2.65 | 2.75 | ||||
K165 | 2002 | AHI | 2.67 | ||||
K169-4058 | 2011 | 2.69 | 3.26 | ||||
K177-4329 | 2053 | AHCS | 1.49 E + 06 | 1.75 | 1.63 E − 20 | ||
K178-4329-4 | 2070 | 9.1 | |||||
K178 | 2071 | 10.58 | |||||
K179-4349 | 2076.4 | 0.18 | |||||
K191 | 2142 | Sain | 2.59 | ||||
K193 | 2152 | AHI | 4.26 E + 06 | 3 | 3.92 | 7.37 E − 18 | |
K193 | 2152 | AHI | 5.74 E + 06 | 3 | 5.07 | 1.70 E − 18 | |
K193-4693 | 2154.9 | ||||||
K195-4755 | 2168 | 2.67 | 3.7 | 2.3 | |||
K195-4755b | 2.66 | 3.85 | |||||
K195-4772 | 2160.5 | 2.53 | |||||
K195-4774 | 2161 | 7.07 E + 06 | 2.48 | 8.62 | 1.44 E − 18 | ||
K195-4774 | 2161 | 1.00 E + 07 | 9.87 | 7.59 E − 19 | |||
K195-4774 | 2161 | 8.49 E + 06 | 7.83 | 5.38 E − 19 | |||
K195-4777 | 2162 | AHCS | 3.49 | ||||
K195-4782 | 2163 | AHI | 2.77 | ||||
K195-4779-1 | 2174 | 2.53 | |||||
K195-4779-2 | 2174.4 | 2.53 | |||||
K195-4779-3 | 2174.4 | 2.53 | |||||
K197 | 2179.7 | 2.69 | 3.12 | ||||
K206-4984 | 2214.5 | 0.29 | |||||
K207-5009-1 | 2220.5 | ||||||
K207-5009-2 | 2220.6 | 2.67 | 3.4 | ||||
K207-5019 | 2222 | 0.56 | |||||
K207-5028 | 2224.5 | Fresh | 2.69 |
The petrographic study revealed three types of material: fresh granite, pervasive alteration granite and a fractured-altered granite, (Genter et al., 1990; Genter and Traineau, 1992; Jacquemont, 2002; Ledésert et al., 1993; Rosener, 2007; Sausse et al., 1998; Surma, 2003). Quartz, K-feldspar, plagioclase, biotite and a high density of K-feldspar phenocrysts compose the fresh granite. The first type of alteration is characterised by transformation of primary minerals like biotite, amphibole and plagioclase. Biotite is mainly altered into chlorite and sometimes in celadonite, epidote, hydrogarnet, and carbonates. Plagioclase is mainly altered into corrensite, calcite or illite. The second type of alteration is located in the matrix of the damaged zone and controlled by fluid circulation in and around fractures. The main minerals in the fractures are phyllosilicates (illite, chlorite, tosudite), carbonates and quartz, where the matrix shows mostly illite and carbonates as secondary phases. Schematically, each facies could be characterized by clay phases; the protolith is defined by chlorite formation, the damage zone by the chlorite-illite association, and the fault core by illite. Other clay phases or secondary phases are also present but we choose this simple characteristic for a quick analysis of the samples and their location in the fault model (Table 1, Fig. 1).
3 Techniques
Petrography and chemical data are available through some publications (Genter et al., 1990; Genter and Traineau, 1992; Jacquemont, 2002; Ledésert et al., 1993; Sausse et al., 1998; Surma, 2003). Petrophysical properties are measured on samples selected along the EPS1 core. Four of these are determined: porosity, permeability, thermal conductivity (TC), and specific surface. These parameters control chemical and physical fluid-rock interactions and could be later use to model these processes. All of these measurements are not systematically performed on the same samples, Table 1 summarizes all of obtained data.
Porosity measurement and threshold diameter determination are obtained from a mercury injection test. This technique is based on the Young-Laplace equation for the displacement of a non-wetting fluid in a thin capillary tube, relating the fluid pressure to the throat size. This technique gives access to the total porosity value, but also provides information about the distribution of pore size and pore access diameters. Volume of the samples varies between 3 and 5 cm3. Due to grain size issues, several samples are used to characterize each facies, in a way to analyse a representative volume.
Permeability measurements are made using a nitrogen gas permeameter working at 40 bars confining pressure. Samples are 18 mm in diameter and 25 mm length. For different constant head pressures (lower than 35 bars) and once the steady-state condition is reached, gas flow rates and pressures from the both sides of the samples (the head and foot pressures) are recorded. These values must be corrected for the Klinkenberg effect to obtain the real permeability value (Scheideffer, 1974; Debschutz et al., 1989; Klinkenberg, 1941). For several samples, these experimental conditions are not sufficient to measure flow and they are defined as impermeable, with permeability lower than 10−20 m2. It is also important to notice that these experimental conditions are much different from the deep in situ conditions.
To determine the TC, we use the optical scanning technique developed by Popov et al. (1999). With a measurement taken every 1 mm, TC normal to the analysed surface is determined as an arithmetic average of local conductivity along an entire scanning line. The relative measurement error is not more than 2% (Popov et al., 1999). Measurements are performed on samples that are 10 to 20 cm long; some of them were re-sampled afterwards for permeability measurements. Compared to other techniques, the optical scanning allows both average and pinpoint evaluations.
The BET method (Brunnauer et al., 1938) is used for specific surface area (SSA) measurements. At low temperature and low pressure, gas is adsorbed on the porosity network surface (Gregg and Sing, 1982). By measuring the quantity of gas adsorbed at different pressures, it is possible to estimate the specific surface of the sample (in m2 g−1 or m2 m−3). Nitrogen is used for our measurement.
Using SEM observations, void shapes and voids relationships with mineral components are determined to define the link between the different physical properties. Observations are performed on ships of rocks under a GEOL scanning electronic microscope.
4 Results
All of the results are summarized on the Fig. 2. Two main types of information come from the mercury injection technique, the porosity value and the distribution of the pore volume versus the threshold value. Porosity ranges from 0.5% to up to 15% and its variation is not depending on the depth of sampling. It is controlled by the location of the sample within the fault zone. Samples from the protholith have porosity as low as 0.5%. In the damage zone, porosity can increase up to 15%, and lower values can also be measured near the fault core. Porosity in the fault core varies over two orders of magnitude.
With regard to the EPS-1 well, four high porosity zones are recognized, three of them linked to fracture zones at 1620, 2070 and 2120 m depth (Fig. 2). The fourth is located at the upper part of the batholith on the Hercynian unconformity, where alteration of the granite's top and post Permian/Triassic fluid flow induced high porosity at the granite-sedimentary cover contact. Zones with lower porosity are also localized along secondary fracture zones, as at 1800 m depth. Low fracture density characterises these secondary fracture zones.
The porosity increase could be linked to an increase of the threshold values (Fig. 3). Different porosity spectra could be recognized, from relatively flat spectra for the low porosity material to a porosity network controlled by a particular threshold value, where the porosity increase is associated with the threshold increase. These porosity variations are measured on small volume samples (3 or 5 cm3) and describe matrix porosity and microcracks porosity development.
Permeability varies over five orders of magnitude from 10−20 m2 to 10−15 m2 (Fig. 2). Several high permeability zones are recognized with permeability values higher than 5.10−17 m2. In several cases, these zones could be related to high porosity zones like at 1620 or 2120 m depth, but in other cases, these high permeability zones are not linked to high porosity like at 1800 or 1900 m depth. Conversely, high porosity zones are not related to high permeability volume. This could be due to lack of permeability samples or to low cracking of the matrix around different fault zones. In the alteration zone of the upper part of the batholith, a large variation of the permeability is measured over three orders of magnitude. The permeability–porosity relationship is complex (Fig. 4) and depends on the pore structure, voids geometry and location (Fig. 5). Three types of behaviour could be distinguished depending on the porosity shape (Rosener and Géraud, 2007). For the cracks network (group 1), a sharp increase of the permeability is measured for a low porosity increase. The second group (group 2) is composed of samples where alteration and secondary phase precipitation modify the pore structure. Under these conditions the permeability decreases with porosity. Samples of the group 3 are intensively cataclased and illite and chlorite precipitate in the enlarged cracks and pores by mechanical and chemical processes. For low porosity samples, the permeability increases with porosity, and for porosity above 4% the permeability decreases. These particular relationships (decrease) between the porosity and the permeability are clearly the illustration of the tortuosity effect controlled by the infilling development in the cracks networks and the alteration of the primary phases. Especially for the group 3, the threshold value being constant for all the samples of the group, porosity parameters (volume and threshold) increase and the permeability decreases. They illustrate the influence of the tortuosity increase.
The SSA is not measured for all the samples, but a similar tendency like the ones described for permeability and porosity is observed, i.e. a large variation is measured for samples coming from the fractured zone. However, a simple relationship with permeability or porosity could not be found. The highest values are measured on samples coming from two fracture zones at 1660 and 2100 m depth. The altered zone in the upper part of the batholith gives also samples with high SSA value. These high values are directly linked to the clay content of the sample, whereas the lowest SSA is measured for the unaltered and little fractured samples (Rosener and Géraud, 2007).
TC varies between 2.3 W K−1m−1 and 3.7 W K−1m−1. TC for the fresh granite lies between 2.5 W K−1m−1 and 2.7 W K−1m−1, in accordance with the literature (Cermak and Rybach, 1982; Pribnow and Sass, 1995; Rummel, 1991; Surma and Géraud, 2003). As for the other properties, TC presents large variations within the fractured zones. Decrease of the TC value is the result of the porosity development by fracturing or by mineral alteration. This process is particularly obvious across the altered zone at the top of the batholith, where the lowest values are measured. On the other hand, a large part of the tested samples shows an increase of the TC value with alteration. This evolution is clearly associated to secondary phase precipitations, especially quartz with a TC value of 7.69 W K−1m−1 and calcite with a mean TC value 5 W K−1m−1. This TC increase could mask the effect of the porosity development.
5 Discussion: a fault zone model
From these data, physical properties of a theoretical fault zone may be proposed. A first qualitative geometry of the fault zones in the Soultz's batholith was proposed by Genter et al.(1995). Here after, our model is the synthesis of data coming from different Phd theses (Jacquemont, 2002; Rosener, 2007; Surma, 2003), their mineralogical and chemical point of view is developed in Fritz et al. (this volume). Using the model proposed by Caine et al., 1996 and Sibson and Rowland, 2003, a more refined synthesis of the fault zone could be drawn and analysed samples could be localised on this model. Essentially, a set of altered conjugate structures with lower offset is associated to the principal fault.
The central part of the main structure is described as the fault core: a highly strained material, with large ranges of porosity and permeability depending on the degree of cementation. It can also include locally open sections that could support high fluid flow after being reconnected during the stimulation. The second part is the damage zone, where previously fractured material is subject to alteration and secondary precipitations. Generally, the increase of porosity induced by fracturing and mineral dissolution is higher than the decrease of porosity associated with secondary mineral precipitation. However, due to their nature (mostly clays) secondary precipitations generate high tortuosity, and so tend to decrease the permeability of the rock (group 3 material). In these damage zones, there are also some preserved rock volumes of which some small porosity increase related to sharp permeability increase is unclear but probably induced by a crack network development (group 1). These volumes (fault core and damage zone) are generally not connected to the protolith or other preserved blocks. The last part of the system is the protholith that remains largely unaffected by the fault.
The global geometry and physical properties profile across the fault zone could be described as follows: within the damage zone, from the protolith to the fault core, permeability increases due to fracturing and primary phase alteration. Then, toward the fault plane, it could be reduced by secondary mineral precipitation (Lockner et al., 2000). Between the two surfaces of the fault plane, occasional free space could represent a high fluid flow conduit with high permeability; where in other sections, this volume could be filled more or less completely, possibly showing a lower permeability than the surrounding material. In the last case, the damage zone would then be the section showing the highest permeability. Finally, connections between the fault core and the damage zone are made possible via several secondary structures that could, at their own scale, present a similar layout as the main fault zone.
This geometry explains the high frequency variation of the different physical properties measured in the fault zones or the fractured zones and also the large variations of transport properties recognized along the borehole. It is also in accordance with the conceptual model proposed by Caine et al., 1996, where the size of the fault core and the damage zone could vary as well as their transport properties.
6 Conclusion
Petrophysical characterizations of different samples coming from the EPS1 core are used to build a fault zone model of the deep geothermal reservoir. The different elements of this fault zone are the open space between the fault planes, the fault core, the damage zone and the protolith. The characteristics of each part point out the role of the damage zone as a high porosity and permeability zone. In particular, it could be a pathway for the fluid if the fault core is filled by impermeable secondary phases. Across the fault zone, even if structure sizes are variable, petrophysical properties of the damage zone seem to be more consistent than in the fault core, where fluid-rock interactions are more intense and induce bigger variations.
Acknowledgements
F. Surma and M. Rosener were funded by “Région Alsace” and ADEME for their PhD; ADEME and Eifer supported J. Place for his PhD. This research program was supported by a European Community grant “STREP”. “Région Alsace” contributed to the founding of the Thermal Conductimeter. Many thanks to all of them for their constant support. Many thanks to the both anonymous reviewers for their constructive remarks.