Hydraulic fracturing at the European HDR/HFR test site Soultz-sous-Forêts (France) a conceptual model

Similar documents
New Concepts for Extracting Geothermal Energy from One Well: The GeneSys-Project

Review of the Hydraulic Development in the Multi-Reservoir / Multi-Well EGS Project of Soultz-sous-Forêts

Synergetic CO2 Storage and gas/oil production from tight reservoirs

ERIC ROP RESERVOIR ENGINEER KENGEN, KENYA

NUMERICAL STUDY ON HEAT EXTRACTION FROM SUPERCRITICAL GEOTHERMAL RESERVOIR

Review of the Ogachi HDR Project in Japan

Micro-earthquake Observation Results in the Ogachi HDR Project

Flow Testing Results from Habanero EGS Project

A Simple Engineering Analysis Method to Evaluate the Geothermal Potential from a Hot Water Aquifer

Quantification of Exploration Risks for Hydrogeothermal Wells

Engineering Geothermal Systems in Oil and Gas Reservoirs. Susan Petty Bill Livesay Black Mountain Technology

UTILITY OF THE DATA GATHERED FROM THE FENTON HILL PROJECT FOR DEVELOPMENT OF ENHANCED GEOTHERMAL SYSTEMS

New Large Geothermal Storage Project in Zurich/Switzerland for Heating and Cooling

Multiple Zone Geothermal Stimulation Case Study: Fluid Diversion Using Thermo-Degradable Zonal Isolation Materials

REMOVAL OF WATER FOR CARBON DIOXIDE-BASED EGS OPERATION

Simulation study of hydraulic fracturing in carbonate oil reservoir by using KGD model

EFFECT OF NATURAL CONVECTION PATTERNS ON OPTIMAL LOCATION AND SIZE OF A HEAT SINK IN A GEOTHERMAL RESERVOIR

Petroleum Production Systems

Geothermal reservoir simulation of hot sedimentary aquifer system using FEFLOW

NET POWER CAPACITY OF GEOTHERMAL WELLS VERSUS RESERVOIR TEMPERATURE A PRACTICAL PERSPECTIVE

ECOGI : an EGS project for the industry in the Upper Rhine Graben. Martino LACIRIGNOLA, M.Eng. ADEME French Environment and Energy Management Agency

PRESSURE TRANSIENT BEHAVIOUR OF A HORIZONTAL STEAM INJECTION WELL IN A NATURALLY FRACTURED RESERVOIR

Geothermal reservoir candidates in deep crystalline and sedimentary formations: tracer-assisted evaluation of hydraulic stimulation tests

3D simulations of an injection test done into an unsaturated porous and fractured limestone

Well Productivity Enhancement by Drilling Multi-Legged Wells A Quantitative Assessment

- Case study Groß Schönebeck -

Numerical study of the influence of injection/production well perforation location on CO 2 -EGS system

Metal sealing packer technology in deep boreholes

SMU Geothermal Conference March 13, 2013

Dr. J. Wolters. FZJ-ZAT-379 January Forschungszentrum Jülich GmbH, FZJ

Analysis of Sedimentary Geothermal Systems Using an Analytical Reservoir Model

An experimental study of permeability determination in the lab

Erik Lindeberg and Per Bergmo. SINTEF Petroleum Research, NO-7465 Trondheim, Norway

Temperature Plume Migration in Aquifers: The necessary first step to geochemical evaluation of thermally-mobilized constituents

Simulation of Fracture Design Generation, Production Characteristics and Temperature Development of a Hot Dry Rock Geothermal Reservoir.

On Minimizing the Levelized Cost of Electric Power from Enhanced Geothermal Systems

THE DEVELOPMENT OF A COMPUTERIZED PRESSURE TRANSIENT TEST SYSTEM (PART 1)

Unconventional Reservoir Stimulation. R. D. Barree Barree & Associates LLC

Geothermal Energy Storage

Geothermal Well Productivity: Why Hotter is Not Always Better

A NOVEL TECHNIQUE FOR EXTRACTION OF GEOTHERMAL ENERGY FROM ABANDONED OIL WELLS

NOTICE CONCERNING COPYRIGHT RESTRICTIONS

Optimization of Experimental Parameters for Comparing Enhanced Geothermal Working Fluids in the Lab

Flexural behaviour of a polymeric foam/glass-fibre composite: FE modelling and experimental testing

Performance of the Subsurface Hydraulics in a Doublet System Using the ThermoGIS Calculator

Real Time Performance of Injection Wells Andrey Ryazanov. FORCE seminar on Water Injection Safety 4 December 2013

CORE ANALYSIS AS A KEY TO UNDERSTANDING FORMATION DAMAGE AFTER HYDRAULIC FRACTURING TREATMENT

Modeling of Usage of Air Injection Well in a Geothermal System

Laboratory Testing of Safety Relief Valves

A NOVEL CONCEPT TO ENGINEERED GEOTHERMAL SYSTEMS

Permeability, Flow Rate, and Hydraulic Conductivity Determination for Variant Pressures and Grain Size Distributions

Chapter 2: Mechanical Behavior of Materials

Discrete Fracture Embedding to Match Injection Returns in a Single Porosity Model

Simulation of a Novel Groundwater Lowering Technique using. Arbitrary Lagrangian-Eulerian (ale) Method

Influence of Vertical Acceleration on Seismic Response of Endbearing

4 Image Analysis of plastic deformation in the fracture of paper

An Efficient and Optimized Approach to Multiple Fracturing of Horizontal Wells

New techniques of coal seam panels' pre-drainage - based on project experience

Prediction of Water Production in CBM wells

Diagnostic plots applied to well-tests in karst systems

Interpretation of the Injection and Heat Up tests at Sabalan geothermal field, Iran

The model was developed on the TETRAD reservoir simulator.

EFFECT OF LOCAL PLASTIC STRETCH OM TOTAL FATIGUE LIFE EVALUATION

Production prediction model of coalbed methane wells at stable production stage

IMPROVING GEOTHERMAL PROJECT ECONOMICS WITH MULTI-ZONE STIMULATION: RESULTS FROM THE NEWBERRY VOLCANO EGS DEMONSTRATION

Simulation of Pumping Induced Groundwater Flow in Unconfined Aquifer Using Arbitrary Lagrangian-Eulerian Method

Horizontal Well Spacing and Hydraulic Fracturing Design Optimization: A Case Study on Utica-Point Pleasant Shale Play

REVIEW ON SHEAR SLIP OF SHEAR KEYS IN BRIDGES

POROSITY, SPECIFIC YIELD & SPECIFIC RETENTION. Physical properties of

PHASES OF GEOTHERMAL DEVELOPMENT

Using underground experiments to improve the understanding of induced seismicity. Domenico Giardini, Florian Amann & the DUG-Lab group

INJEÇÃO DE CO 2 PARA PRODUÇÃO ACRESCIDA DE METANO DE CARVÃO EM CAMADA CO 2 INJECTION FOR ENHANCED COALBED METHANE (ECBM)

Geothermal Power Capacity of Wells in Non-Convective Sedimentary Formations

THE UK GEOTHERMAL HOT DRY ROCK R&D PROGRAMME

Proposed Enhanced geothermal system demonstration at the Nevada Test Site

Numerical Modeling of Slab-On-Grade Foundations

Research Article Fuzzy Comprehensive Evaluation on the Effect of Measures Operation for Oil-Water Well

Standard Test Method for Determining Transmissivity and Storativity of Low Permeability Rocks by In Situ Measurements Using Pressure Pulse Technique 1

CFD Simulation of Pore Pressure Oscillation Method for the Measurement of Permeability in Tight Porous-Media

An Alternative and Modular Approach to Enhanced Geothermal Systems

8.0 FEASIBILITY OF PASSIVE SITE REMEDIATION. 8.1 Introduction

THERMAL EXTRACTION ANALYSIS OF FIVE LOS AZUFRES PRODUCTION WELLS. Luis Quijano GPG Com. Fed. de Elec. Morelia, Mich.

ACARP Project C Hydraulic Fracturing for Gas Recovery and Its Impact on Subsequent Coal Mining

EVALUATION OF MICROGRAVITY BACKGROUND AT THE UNDISTURBED OGUNI GEOTHERMAL FIELD, JAPAN

Comparison of Borehole Heat Exchangers (BHEs): State of the Art vs. Novel Design Approaches

Simulation of Contaminant Transportation in Saturated Transparent Rock Fractures

Baigcng Hu. Thermal Engineering Depaarnent, Tsinghua University, Beijing, , P. R. China

Analysis of Fracture Design in a Vertical Well with a Horizontal Fracture

Considering the Local Anisotropy of Short Fibre Reinforced Plastics: Validation on Specimen and Component

EGS Update Annex III September 21, 2103

Three-dimensional computer simulation of soil nailing support in deep foundation pit

Perforating for Hydraulic Fracturing Next Generation Technologies for Unconventional Completions

Numerical Analysis of MCY01 Production Well Mechanism: Understanding of Well Coughing Phenomenon and Irregular Productivity

EGS in the Permian of the North German Basin, Europe: a borehole doublet utilizing a former exploration well

INCORPORATION OF GROUNDWATER FLOW INTO NUMERICAL MODELS AND DESIGN MODELS

Reservoir Simulation & Forecasting Course Outline

Proxy Model Based on Butler s SAGD Theory

Geomechanical Effects of CO 2 Injection under Thermal and Fracturing Conditions

APPROPRIATE USE OF USGS VOLUMETRIC HEAT IN PLACE METHOD AND MONTE CARLO CALCULATIONS

Physical Simulation for Water Injection of Horizontal Wells in Naturally Fractured Buried- Hill Reservoirs

Transcription:

Hydraulic fracturing at the European HDR/HFR test site Soult-sous-Forêts (France) a conceptual model G. Grecksch 1, R. Jung 1, T. Tischner and R. Weidler 3 1 Leibni Institute for Applied Geosciences GGA, Stilleweg, 30655 Hannover, Germany, grecksch@gga-hannover.de, r.jung@gga-hannover.de Federal Institute for Geosciences and Natural Resources BGR, Hannover, Germany 3 Q-con GmbH, Bad Bergabern, Germany Abstract The HDR/HFR (Hot-Dry-Rock, Hot-Fractured-Rock) technique allows the extraction of heat from deep crystalline rock by circulating fluids through an underground heat exchanger which is formed by a hydraulically stimulated fracture system. During the last 10 years massive hydraulic fracturing tests have been performed successfully at the European test site Soult. By injection of water at high rates and pressure, large fracture systems could be generated in the crystalline basement between 3 and 5 km depth. A clear understanding of the hydro-frac process is necessary for the development of more efficient stimulation techniques in order to achieve higher flow rates and larger fracture surfaces. Thus, a conceptual model for hydraulic fracturing is proposed which has been developed based on hydraulic and seismic observations. It is assumed that the hydro-frac process creates tensile artificial fractures within or between natural faults and major fracture ones. The opening of artificial fractures during hydraulic stimulation leads to irreversible shearing along the natural faults. The stimulation experiments at Soult revealed that a self-propping mechanism kept the fractures open. This observation is attributed to the shearing mechanism. Preliminary results from analytical and numerical model calculations confirm the validity of the conceptual model. Keywords: HDR/HFR/EGS, hydraulic fracturing, stimulation, bilinear flow, modeling Introduction The European HDR/HFR test site is located in the Rhine graben near Soult-sous-Forêts, northern Alsace, France. The overall aim of the European HDR programme is the realiation of power production by extracting geothermal energy from hot deep crystalline rocks. Fluid will be circulated through a 3-borehole system (1 injection, extraction boreholes) connected to each other by a large underground heat exchanger in 5 km depth. A scientific pilot plant for geothermal power production of about 6 MWe will be realied. The heat exchanger has to be created artificially by injection of fluid at high flow rates and pressure. These hydraulic fracturing tests (hydro-fracs) are a common method to increase permeability of a reservoir. The conventional hydro-frac technique includes the injection of proppants to keep the artificially created fractures open [1]. For geothermal purposes the use of proppants is not recommended because of their limited penetration depth and their relatively high costs. Hydraulic fracturing tests using water without proppants require an irreversible shearing

process for persistent fracture opening (self propping effect). Since 1993, several water-frac tests have been carried out successfully in the Soult boreholes GPK1 and GPK. Two hydraulically separated artificial fracture systems could be created in dephts of about 3500 m and 5000 m. In 1997, a four-month circulation test in the shallow reservoir was carried out sucessfully []. After deepening GPK down to a depth of 5 km in 1999, a deep reservoir was stimulated which is hydraulically separated from the shallow reservoir. A temperature of 00 C was reached. In 00, a second 5 km deep well, GPK3, was completed and will be stimulated in 003. The injectivity and productivity of the boreholes GPK1 and GPK could be enhanced by more than one order of magnitude [3] showing that the water-frac technology is suitable for the creation of large geothermal underground heat exchangers. In order to achieve higher flow rates and larger fracture surfaces for future commercial geothermal power plants, an enhancement of the efficiency of hydraulic stimulation methods is essential. A clear understanding of the hydraulic fracturing process is necessary for the development of more efficient stimulation techniques. As contribution to the improvement of understanding, a new conceptual model for the hydro-frac process was proposed by Jung and Weidler [4]. The model has been developed based on hydraulic and seismic observations in Soult (shallow reservoir). Preliminary results from analytical and numerical modeling addressing the validity of the conceptual model will be presented. Hydraulic and seismic observations The hydraulic test programm of the deep reservoir is running at present and revealed indications for a relatively tight and fault dominated reservoir system [3]. These observations are preliminary and we focussed this analysis on the shallow reservoir. The main results of hydraulic tests in the boreholes GPK1 and GPK between 1993 and 1996 are as follows: The injection overpressure during water-frac experiments was found to reach maximum values of 1.5 MPa. This is in agreement with the minimum horiontal stress at Soult [5] and indicates a pressure limitation due to normal opening of existing fractures or tensile fracture propagation. Pressure build-up and decline curves recorded during frac and post-frac injection tests exhibit a transient pressure behaviour changing linearly with the fourth root of time (Fig. 1a) which is characteristic for bilinear flow [6]. The bilinear flow period comprises two simultaneous linear flows, one occurs within a fracture of finite conductivity, the other permeates into the formation, perpendicular to the fracture plane. Fig. 1a shows a bilinear flow period persisting for approximately one day. All post-frac hydraulic tests with step-like flowrate history displayed a quadratic relationship between downhole overpressure and corresponding flowrate (Fig. 1b). This behaviour indicates the presence of turbulent flow [7]. The post-frac injectivity or productivity shows a linear relationship with increasing flowrate applied during water-frac stimulation (Fig. a). Massive hydro-frac stimulation produces micro-seismic events. Their spatial distribution is assumed to represent the shape and sie of the stimulated fracture system. Fig. (b) displays the top view of the recorded seismic cloud (after [8]). Some linear structures are obvious.

7 6 (a) 14 1 (b) differential pressure (MPa) 5 4 3 1 0 differential pressure (MPa) 10 8 6 4 0 95JUN16 (frac) 95JUN10 (pre-frac) 95JUL01 (post-frac) 96SEP9 (post-frac) 0 5 10 15 0 0 10 0 30 40 50 60 time^1/4 (s^1/4) flowrate (l/s) Fig. 1: (a) Pressure decline during a hydraulic post-frac injection test in the Soult borehole GPK in 1996 (test 96SEP9, 0 s 1/4 =shut-in). The decline curve shows a fourth root of time behaviour which indicates the presence of a distinct bilinear flow period. (b) Downhole overpressure during hydraulic tests vs. corresponding flowrates. Post-frac curves show a quadratic behaviour with increasing flowrate. 14 injectivity or productivity (l/(mpa s)) 1 10 8 6 4 0 (a) GPK1, 1993 GPK, 1996 GPK, 1995 (b) 0 10 0 30 40 50 60 70 80 flowrate during stimulation (l/s) Fig. : (a) Post-frac injectivity or productivity is linearly related to the flowrate applied during hydro-frac stimulation. (b) Locations of micro-seismic events during hydro-frac experiments in the Soult boreholes GPK1 and GPK, orientation of minimum horiontal stress S h. Linear structures are obvious in the seismic cloud [8].

A conceptual model for hydraulic fracturing A clear understanding of the physical processes during hydraulic fracturing is essential for the development of more efficient stimulation techniques. A conceptual hydro-frac model was developed by Jung and Weidler [4] which is based on the hydraulic and seismic observations in the shallow reservoir. Fault or fracture ones of a scale of at least several 100 m to some km are believed to represent the important natural discontinuities in a crystalline basement. The conceptual model assumes that a tensile artificial fracture is created within or between these natural fault ones (Fig. 3). The observed increase of the borehole productivity and injectivity is attributed to this mechanism. The vertical growth of artificial fractures is stopped as soon as they intersect the natural discontinuities. The fracture heights may amount to some 10 m. The fracture length may reach several 100 m or even more than 1 km and is limited due to fluid losses into the fracture one. The normal opening of artificial fractures leads to irreversible shearing along the natural discontinuities. The stimulation experiments at Soult revealed that a self-propping mechanism kept the fractures open. This observation and also linear structures identified in the seismic cloud are attributed to this shearing mechanism. The hydraulic observations can be explained as follows: The overpressure during the waterfrac experiment is limited due to normal opening of existing fractures or tensile fracture propagation. The natural and artificial fractures together form a bilinear flow system as indicated by the arrows in Fig. 3. Due to the limited fracture height the flow velocity within the artificial fracture is high which leads to turbulent flow conditions. The linear increase of post-frac injectivity or productivity with increasing flowrate can be understood as result of the fracture propagation that is controlled by the fluid losses into the matrix. The propagation stops when the amount of fluid losses equals the injection flowrate. Thus, the fracture lengh L and the post-frac injectivity or productivity is proportional to the injection flowrate during the water-frac test. Fig. 3: Proposed conceptual model of hydraulic water-frac stimulation in a crystalline rock mass: During the hydro-frac process, a tensile artificial fracture is created within or between natural fracture ones. L is the fracture lengh, h its height, and d is the penetration depth of the pressure front into the fault one (from [4]).

Estimation of bilinear flow parameters After an analytical solution for bilinear flow developed by Cinco-Ley and Samaniego-V. [6], the differential pressure curve is given by 1/ 4 P = mt with 34. 97qµ B m = (eq. 1) 1 ( T h) / ( T S µ ) 1/ 4 f and with differential pressure P [kpa], time t [hr], flowrate q [m³/d], dynamic fluid viscosity µ [Pa s], formation volume factor B [m³/m³], fracture transmissibility T f [Dm], fracture height h [m], total transmissibility of the fault one T [Dm], and total storage coefficient of the fault one S [m/pa]. Eq. 1 was developed for laminar flow behaviour, but the bilinear flow period seems to be preserved even under turbulent fracture flow conditions [7]. From eq. 1, the product T f h can be expressed by (B = 1) 34.97 q µ T f h =. (eq.) m ( T S µ ) 1/ The Reynolds number Re indicates the beginning of turbulent flow. For rough fracture surfaces, deviations from laminar Darcian flow may be found at Reynolds numbers as low as 100 [7]. From Re, a maximum fracture height h can be estimated as follows qc ρ q ρ Re = > 100 h < c (eq. 3) hµ 100µ with a flowrate of q c at the onset of turbulence. Combining eqs. and 3, a minimum transmissibility T f can be calculated (rf. Table 1). The minimum effective fracture width is given by the cubic law [9]: ( T ) 1/ 3 w eff = 1 f. (eq. 4) The half lengh of the artificial fracture L/ can be determined from the duration of the bilinear flow period which ends at the time t b [6]: L ( ) 1/ 4 4 1/ 4 ( ) 1/ 10 T h ( T S ) 1/ t = f µ b β (eq. 5) with unit conversion constant β = 3.6e-09. The penetration distance of the pressure front into the surrounding fault one, d/, can be estimated from the hydraulic diffusivity D and may be expressed by ( Dt ) 1/ d = with b T D =. (eq. 6) S µ The parameters for the bilinear flow system as described above (rf. Fig. 3) have been estimated based on data from the post-frac step injection test 96SEP9 (rf. Fig. 1) in the Soult borehole GPK in 1996. Table 1 summaries the adopted parameters and the results.

Input parameters from: Flowrate q 0.0375 m³/s = 340 m³/d test data Dynamic fluid viscosity µ 3.e-04 Pa s assumed Total transmissibility of fault one T < 1.e-1 m³ = 1 Dm assumed Total storage coefficient of fault one S < 1.e-09 m/pa assumed Slope of the fourth root of time curve m 3.e05 Pa/s 1/4 =.3e03 kpa/hr 1/4 Figure 1a Duration of bilinear flow period t b 1/4 t b = 17 s 1/4 =. hr 1/4 Figure 1a Flowrate at onset of turbulence q c 0.005 m³/s [7] Fluid density ρ 1000 kg/m³ assumed Parameters of bilinear flow system Product T f h > 3.9e-10 m 4 eq. Fracture height h < 170 m eq. 3 Fracture transmissibility T f >.3e-1 m³ eq. and 3 Effective fracture width w eff > 0.3 mm eq. 4 Fracture half length L/ > 800 m eq. 5 Penetration depth d/ 530 m eq. 6 Table 1: Estimated parameters of the bilinear flow system. The adopted values of transmissibility T and storage coefficient S represent reasonable maxima. Smaller values of T and S lead to a higher estimate of fracture transmissibility T f and subsequently to higher values of fracture half length L/ and smaller values of penetration depth d/. The estimated values of L/ and d/ are in reasonable agreement with the sie of the stimulated fracture system which is given by the spatial distribution of micro-seismic events recorded during water-frac stimulation. Modeling results Analytical and numerical model calculations addressing the validity of the above described conceptual model have been performed. The post-frac step injection test 96SEP9 was analysed using KAPPA s pressure transient interpretation software Saphir Advanced. The selected analytical model comprises a well intercepted by a vertical finite-conductivity fracture in an isotropic, homogeneous reservoir with planar constant pressure boundary at a distance of some 100 m. The modeled transient behaviour of a vertically fractured well includes linear fracture flow followed by a bilinear flow period, formation linear flow, and finally a pseudoradial flow period. Bilinear flow is observed in finite-conductivity fractures and before fracture tip effects affect the well behaviour [6]. The observed downhole pressure during the step injection test and the fitted curve are displayed in Fig. 4. The model calculation yields a quite well fit of the observed pressure history. This indicates that the observed pressure curve which exhibits a distinct bilinear flow period can be reasonably explained by a well-fracture system as proposed in the conceptual

model of hydraulic fracturing. Nevertheless, further calibrations of the model are necessary in order to determine the hydraulic fracture parameters. Additionally, numerical finite-element modeling based on the above described conceptual model have been started using the FE-code RockFlow. In order to focus on the underlying physical processes, it is necessary to reduce the uncertainties by using simple model geometries. The finite-element model consists of an artificial fracture of reasonable dimensions (rf. Table 1) within a fault one. The fault one is given by a porous permeable matrix representing altered crystalline rock in the fault. The natural discontinuities are replaced by this matrix. Turbulent fracture flow is implemented by Forchheimer s approach. Preliminary results based on the flowrate history of the injection test 96SEP9 also suggest that the FE model is suitable to explain the observed pressure curves. differential pressure (MPa) 7 6 5 4 3 1 0 6. 15 5 0 40 80 10 160 00 40 time (hr) Q (l/s) 37.5 pressure (MPa) model fit Fig. 4: Analysis of hydraulic post-frac test 96SEP9 in the Soult well GPK; downhole differential pressure in MPa (gauge depth 300 m), and fitted pressure (vertically fractured well, KAPPA s Saphir Advanced) vs. time; injection rate Q in liters/sec vs. time. Conclusions The water-frac technology is a suitable method to create large scale underground heat exchangers for geothermal energy use. In order to develop more efficient stimulation techniques, a clear understanding of the water-frac process is necessary. Main hydraulic and seismic observations during hydro-fracs at the HDR/HFR test site Soult and a conceptual model for the hydro-frac process as proposed by Jung and Weidler [4] have been reviewed. The conceptual model assumes the creation of a tensile artificial fracture

within a natural fault one or between major natural discontinuities. The fracture opening leads to irreversible shearing along the natural faults and a self-propping mechanism kept the fractures open. The conceptual model is in good qualitative agreement with the hydraulic observations at Soult. However, the spatial distribution of micro-seismic events suggest a stimulated fracture system which is much more complicated. But the conceptual model seems to reflect the underlying mechanisms of hydraulic fracturing at Soult. Preliminary results from analytical and numerical model calculations indicate that the observed post-frac transient pressure behaviour can be reasonably explained by the proposed conceptual model. Acknowledgement The present work at the Leibni Institute for Applied Geosciences GGA in frame of the European deep geothermal HDR project Soult is performed under the financial support of the EU and the Projektträger Jülich (PTJ) on behalf of the Federal Ministry for Economy and Labour (BMWA). References [1] M.J. Economides and K.G. Nolte: Reservoir Stimulation, John Wiley & Sons Ltd, Chichester, 3rd Edition, 000, 657 p. [] J. Baumgärtner, A. Gérard, R. Baria, R. Jung, T. Tran-Viet, T. Gandy, L. Aquilina and J. Garnish: Circulating the HDR Reservoir at Soult: Maintaining Production and Injection Flow in Complete Balance. Proceedings, 3. Workshop on Geothermal Reservoir Engineering, Stanford CA (1998) [3] R. Weidler, A. Gerard, R. Baria, J. Baumgärtner and R. Jung: Hydraulic and Microseismic Results of a Massive Stimulation Test at 5 km Depth at the European Hot-Dry-Rock Test Site Soult, France, Proceedings, 7. Workshop on Geothermal Reservoir Engineering, Stanford CA (00) [4] R. Jung and R. Weidler: A Conceptual Model for the Stimulation Process of the HDR- System at Soult, Geothermal Resources Council Transactions, 4, 143-147 (000) [5] T. Hettkamp, G. Klee and F. Rummel: Stress Regime and Permeability at Soult Derived From Laboratory and In-Situ Tests, Geologisches Jahrbuch, R. Baria, J. Baumgärtner, A. Gérard, R. Jung (Eds.): Int. Conf. 4th HDR Forum, Strasbourg 1998, 165-17 (00) [6] H. Cinco-Ley and F. Samaniego-V.: Transient Pressure Analysis for Fractured Wells, Journal of Petroleum Technology, 1749-1766 (1981) [7] T. Kohl, K.F. Evans, R.J. Hopkirk, R. Jung and L. Rybach: Observation and Simulation of Non-Darcian Flow Transients in Fractured Rock, Water Resources Research, 33, No. 3, 407-418 (1997) [8] J. Baumgärtner, R. Jung, A. Gérard, R. Baria and J. Garnish: The European HDR project at Soult-sous-Forêts: stimulation of the second deep well and first circulation experiments. Proceedings, 1. Workshop on Geothermal Reservoir Engineering, Stanford CA (1996) [9] P.A. Witherspoon, J.S.Y. Wang, K. Iwai and J. E. Gale: Validity of cubic law for fluid flow in a deformable rock fracture, Water Resources Research, 16, No. 6, 1016-104 (1980)