Access the full text.
Sign up today, get DeepDyve free for 14 days.
Karlsruher Institute of Technology (KIT), Adenauerring 20b, 76131 Background: The characterization of fault zones in the Bavarian Molasse Basin plays Karlsruhe, Germany a major role for further geothermal reservoir development. Hence, their identification, Full list of author information is available at the end of the article geological origin, and hydraulic characterization are discussed extensively. Methods: Stress indicators and fractures are interpreted from image and caliper logs of three highly deviated wells at the Sauerlach site. We transform the identified stress field into the borehole coordinate system and compare the observed orientation to the modeled stress field which assumes a homogeneous borehole surrounding. Results: High breakout occurrence, cross-cutting fractures, and a fracture orientation from N-S to NNE-SSW are observed in Sauerlach Th1. In Sauerlach Th2 and Th3, fractures strike primarily ENE-WSW and N-S to NNE-SSW. Drilling-enhanced natural fractures and drilling-induced tensile fractures are observed in all three wells and indicate the orientation of tensile stress at the borehole wall. In Sauerlach Th2 and Th3, stress transformation indicates a S -dir. ~ N 10°E in a strike-slip stress regime. The modeled stress orientations match the observed orientations within the well Sauerlach Th1 if either S -dir. is N 320°E in a strike-slip regime or S -dir. is N 10°E in a H H normal faulting regime. Conclusion: This approach improves the detection of the local stress field especially for non-vertical wells, which has, in combination with the facture pattern, a major impact on the hydraulic system of the geothermal reservoir. Background The use of geothermal energy in the Bavarian Molasse Basin started in the late 1990s when first district heating projects emerged. To encourage further development of geo- thermal energy use, the Bavarian Geothermal Atlas was published in 2004. Indeed, this led to an increased exploration activity for heat and power generation. For geothermal district heating projects, an incentive program followed on the national level, resulting in a total of 12 successfully developed district heating projects with an installed ther- mal capacity of 141.1 MW so far. After the first power generation project in 2004, th three more projects followed after 2008/2009 in the area south of Munich, the installed capacity being 22.36 MW and 42 MW (GTV 2014; Dorsch and Pletl 2012). el th The Bavarian Molasse Basin offers favorable conditions for the utilization of geother- mal energy. Among a total number of 46 wells, 44 were operated successfully. The ini- tial exploration strategy focused on faults and fault zones. They were identified mainly based on old seismic industry data showing an offset of significant reflectors especially © 2015 Seithel et al.; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. Seithel et al. Geothermal Energy (2015) 3:4 Page 2 of 22 at thebaseof the Tertiary strata. With advancing reservoir knowledge resulting from the evaluation of the drilled geothermal wells and the acquisition of new seismic data, the lithofacial aspects of the Malm reservoir were studied. The findings were incorporated in the exploration strategy (Boehm et al. 2012; Lueschen et al. 2014). Nevertheless, the two failed drillings still give rise to questions with respect to diagenetic processes, along with dolomitization as well as hydraulically active fractures and faulting. A key aspect towards an improved evaluation of reservoir rock properties of the Malm is the integration of stress field parameters, since critically stressed fractures are known to be more often hydraulically active than uncritically stressed fractures (Barton et al. 1995). It is well known that the major geothermal provinces are located at structurally favorable settings dominated by step-over/transferring, overlapping, or terminating normal faults. They are preferably striking parallel to the maximum horizontal stress (Faulds et al. 2011). The majority of the Molasse basin fault zones identified are ENE-WSW-trending normal faults (Lemcke 1988). In relation to the roughly N-S-oriented maximum horizontal stress (Reinecker et al. 2010), the N-S or NNE-SSW (approximately N 30°E)-oriented fracture systems show a high tendency to develop tensile or shearing character. Recent seismic in- vestigations revealed inhomogeneities near normal faults, which may have a significant im- pact on the hydraulic properties of the reservoir (Lueschen et al. 2014). In order to study fracture systems, stress indicators, and facies structures within the res- ervoir, image logs can be used (Shamir and Zoback 1992; Hickman and Zoback 2004). This study covers three wells of the Sauerlach site with highly deviated reservoir sections. Their fracture sets as well as stress indicators (drilling-induced tensile fractures (DITF), drilling-enhanced natural fractures (DENF), and breakouts (BO)) are interpreted in the image logs. For highly deviated wells, in situ stress interpretation based on stress indicators is carried out by transforming the regional stress field into the wellbore coordinate system. Our approach is to compare the observed orientations from the different stress indicators with the calculated (modeled) stress orientations. For the regional stress field, we generally assumed a strike-slip stress regime (Reinecker et al. 2010). The difference between the orientation of the interpreted stress indicators and the calculated stress distribution in the well can be determined. By varying the stress regime, the difference between the observed and calculated stress direction is minimized to study the local stress field near the wells. The geothermal Sauerlach project We analyzed data from the geothermal project at Sauerlach which is situated SE of Mun- ich. It is one of the biggest heat- and power-supplying geothermal project in Germany (GTV 2014) and is operated by the Munich utility company Stadtwerke München. Three wells, two for injection (Sauerlach Th2 and Th3) and one for production (Sauerlach Th1), were drilled in 2008 and 2009. The production yield is 110 L/s with a temperature of 140°C from the 4 km deep reservoir (Pletl et al. 2010). Under present conditions, an about 4MW of electrical power and a maximum heat production of 4 MW (4 GWh/a) are e th supplied. The plant has been connected to the grid since spring 2014. All three wells are drilled from one drill pad with a nearly vertical orientation down to 1,000 m depth. Starting with the third well section, they are inclined in the direction of the reservoir targets: Sauerlach Th1 is oriented to NW-SE (N 300°E) with an inclin- ation between 50° and 60° (from the vertical) and reaches the reservoir at a depth of 3,739 m TVD (Figure 1a/red point). Sauerlach Th2 reaches the reservoir at 3,571 m Seithel et al. Geothermal Energy (2015) 3:4 Page 3 of 22 Figure 1 Study area and tectonic background. a. Map with well trajectories of Th1, Th2, and Th3 and the top of the reservoir section marked by red points as well as fault zones of the Geothermie Atlas (STWIVT 2010). b. Map of the region south of Munich with the tectonic sections of the unfolded and folded Molasse Basin, geothermal wells, and the fault zones found in the tertiary structures of the Geothermie Atlas. c. Overview map from Geotis (Schulz et al. 2009) and the investigation area marked by a black asterisk. TVD (Figure 1a/red point) and is directed towards N (N 0°E) with an inclination be- tween 40° and 60°. Sauerlach Th3 reaches the reservoir at 3,872 m TVD and is oriented to SSE (N 160°E) (Figure 1a/red point) with an inclination between 35° and 60°. Directly after drilling, 6-arm caliper (Schlumberger EMS-GR, Houston, TX) and image logs (Schlumberger FMI, Houston, TX) were measured to image the borehole wall (Table 1). In addition, drilling safety tests, such as formation integrity tests (FIT), and the reports of cementation jobs revealed borehole pressure data indicating the lower bound of the minimum principle stress magnitude. Well test analysis shows a considerably higher productivity of Sauerlach Th1 com- pared to Sauerlach Th2 and Th3, although permeability values of the reservoir in all three wells were similar. A dominant radial flow regime underlines the significance of matrix porosity to reservoir characterization in all three wells during expansion of the depression cone. However, additional to its better productivity, Sauerlach Th1 stands out with its very high negative skin effect at the very beginning of pumping, which is an indicator of a good hydraulic connection between the borehole and the reservoir (personal communication Savvatis 2014). Seithel et al. Geothermal Energy (2015) 3:4 Page 4 of 22 Table 1 Database for structural and stress interpretation of well logs Well Tool Hole Hole azimuth/ Depth interval Depth interval diameter inclination (MD) (TVD) Th1 EMS-GR (6-arm caliper) 23″ Vertical 49 to 987 m 49 to 987 m EMS-GR (6-arm caliper) 16″ Vertical 982 to 2,351 m 982 to 2,351 m EMS-GR (6-arm caliper) 12″ Vertical to 230° to 290°/45° 2,346 to 3,990 m 2,351 to 3,739 m FMI 8,5″ 300°/45° to 60° 3,981 to 4,757 m 3,739 to 4,177 m Th2 FMI 8,5″ 0°/45° to 60° 4,025 to 4,850 m 3,571 to 4,086 m Th3 FMI 6,125″ 160°/45° to 60° 4,343 to 5,490 m 3,872 to 4,438 m Reservoir geology The Upper Jurassic (Malm) forms the major hydrothermal reservoir in Bavarian Molasse Basin. In the Munich area, it reaches depths between 2,000 and 4,000 m, while the thick- ness is about 550 m. With temperatures between 80° and 140°C and a low salinity, favor- able conditions exist for heat and power generation (STWIVT 2010). The Malm is located below the Bavarian Molasse Basin as part of the North Alpine Foreland Basin. The Bavarian Molasse Basin evolved as a typically wedge-shaped foreland basin from the Upper Eocene to the Upper Miocene in response to Alpine tectonics accompanied by ero- sion and uplift (Kuhlemann and Kempf 2002). Hydrogeologically speaking, the Malm exhibits the behavior typical of carbonatic aquifers, with a complex interaction of karstification, fracture, and matrix porosity, where structural and lithofacial properties determine the productivity of the wells (Lueschen et al. 2014; Steiner et al. 2014). The Malm lithology is characterized by two major sedimentation cycles with marls and marly limestones in Malm alpha and gamma and predominantly limestones in Malm beta and delta/epsilon. Within the younger Malm zeta, a diversification of the sedimentary environment took place (Meyer and Schmidt-Kaler 1989). During the lat- est Jurassic and earliest Cretaceous, sedimentation changed to a brackish or hypersaline ‘Purbeck’ facies. Since these units are both carbonatic and cannot be separated in the seismic data, the latter is considered as a part of the Malm reservoir (Lemcke 1988). Regarding the lithofacies of the Malm, two types of hyper-facies are distinguished: The so-called ‘bedded facies’ with typically thin-bedded marly or micrits limestones and the so-called ‘mass facies’ with thick-bedded limestone or dolomitic units and reef structures. Boehm et al. (2012) pointed out that in the Munich region, good reservoir properties can be expected from dolomitized limestones in the mass facies of the Malm zeta and in the Malm delta/epsilon due to parts having a favorable matrix porosity. In addition, the massy dolomitized limestones tend to have karstification potential (Stier and Prestel 1991) and also show a higher fracture density under structural impact due to more brittle properties. Structurally speaking, the Bavarian Molasse Basin has undergone different tectonic stages of compressive and extensive regimes (Ziegler 1987). The tectonic map at the base Tertiary published in the geothermal atlas (2010) shows a dominant WSW-ENE system in the Munich area (Figure 1b). It consists mainly of antithetic normal fault systems with displacements of up to 150 m, which are oriented parallel to the Alpine Orogen due to basin subsidence as well as subduction in the south (Bachmann et al. 1987). During the Upper Miocene, the compressive regime started with thrust faulting Seithel et al. Geothermal Energy (2015) 3:4 Page 5 of 22 forming the folded Molasse with still a high pore pressure declining into the unfolded northern foreland (Mueller et al. 1988). It is still being discussed where and to what an extent the younger alpine tectonics af- fects the unfolded foreland, since fracture porosity can be influenced significantly by the activity of faults. Previous publications of Betz and Wendt (1983) and Illies et al. (1981) presented conjugated lineament sets from satellite imagery in the western part and Unger (1999) described a set of strike-slip faults in the eastern part of the Bavarian Molasse Basin. Kraemer (2009) identified strike-slip faults in the SW of Munich close to the Alpine forefront. Megies and Wassermann (2014) analyzed recent seismic activ- ity and concluded a NW-oriented strike-slip movement who identified a N 20° to 45°E fault system from 3D seismic measurements (Lueschen et al. 2014). The present‐day maximum horizontal stress orientation in the Bavarian Molasse Basin is roughly perpendicular to the strike of the Alpine front, indicating that the first-order, far field tectonics in Western Europe is less dominant in the area than the second‐ and third‐order stress pattern (Tingay et al. 2006). A regional study of the stress field in the Bavarian Molasse Basin based on breakout analysis by Reinecker et al. (2010) revealed a rather homogeneously distributed N-S-oriented stress with a dominating strike-slip or thrust faulting stress regime. For the WSW-ENE-oriented fault systems, a limited reactivation potential can be as- sumed. The hydraulic potential of these faults probably is also limited. Consequently, faults related to a strike-slip regime as identified by the analysis of focal mechanisms of seismic events (Megies and Wassermann 2014) would also be of importance to geo- thermal exploration. Methods Determination of the local stress field The stress field can be divided into first-order stresses as a result of forces generated at plate boundaries by e.g. global lithospheric motion and second-order stresses resulting from lithospheric flexure e.g. due to glacial loading and unloading or lateral density contrasts (Zoback 1992). On a local scale, third-order stresses in sedimentary basins are controlled by geological structures (Bell 1996a). On the reservoir scale, active faults, lateral or vertical contrasts in material parameters (for example, salt structures, de- coupling horizons) can lead to mechanically modified parameters and deviations of the stress orientation and tectonic regime from the regional or plate-wide stress pattern (Tingay et al. 2006). According to Anderson (1951), three tectonic regimes can be distinguished based on the magnitudes of the vertical stress (S ), the maximum horizontal stress (S ), and the v H minimum horizontal stress (S ). In a normal faulting regime (NF), S is the maximum h v principle stress (S ). It is the intermediate principle stress (S ) in a strike-slip regime 1 2 (SS) and the minimal principal stress (S ) in a thrust faulting regime (TF). In regions with little topography, it is common to assume S to be a principal stress. In areas where no data of the stress magnitude are available, a method known as fric- tional equilibrium theory can be applied to reduce the range of possible stress ratios (Sibson 1974). Differential stress (S -S ) magnitudes in the brittle crust are limited by 1 3 the frictional strength of optimally oriented faults (Jaeger et al. 2007). Byerlee (1978) has shown in laboratory tests that the coefficient of static friction (μ) is independent of Seithel et al. Geothermal Energy (2015) 3:4 Page 6 of 22 the rock type but depends on the magnitude of normal stress. Accordingly, for normal stress bigger than 200 MPa, the coefficient of static friction of μ = 0.6 fits best and for lower normal stress μ = 0.85 can be applied. On the basis of the frictional equilibrium theory, stress polygons for frictionally stable areas in different stress regimes and pos- sible horizontal stresses (S and S ) can be defined (Zoback et al. 2003). In Figure 2a, H h the stress polygon for both static friction coefficients (μ = 0.6 and μ = 0.85) is shown. The determination of the stress components requires a number of theoretical consid- erations and a reliable data assessment. The S component can be estimated from the density of the overburden. Pore pressure (P ) can be estimated by hydraulic pumping tests or pressure measurements in the wells, and S can be estimated from the vertical stress by making certain assumptions (e.g. no lateral strain), if measurements, such as mini-fractures or leak-off tests, are not available (Hubbert and Willis 1972; Zoback 2007). If FIT or cementation reports are available, it is also possible to calculate a min- imal value for the least principal stress, which helps to define a lower limit of the stress gradient (Zoback et al. 2003). In this study, stress-induced phenomena, such as DITF, DENF, and BO recorded at the borehole wall by both FMI and 6-arm caliper logs, were used. The stress state in a borehole with its circumferential/tangential (σ ) and radial stress (σ ) components are θθ rr described by Kirsch (1898) for a vertical well. For inclined wells, further computations with an analytical solution for the stress distribution are needed (Hiramatu and Oka 1968). According to these computations, the orientation of BO and DITFs does not only vary with the stress magnitudes but also with the orientation of the borehole tra- jectory relative to the stress tensor (Mastin 1988). Peska and Zoback (1995) described a method to determine the direction of the maximal (σ ) and minimal tangential tmax Figure 2 Stress field investigations. a). Stress polygon for a reservoir depth of 4,000 m TVD. Different stress regimes are illustrated (NF normal faulting (green), SS strike-slip faulting (yellow), TF thrust faulting (red)) with differential stress ratios (Anderson 1951) for μ = 0.85 (solid line) and (Bachmann et al. 1987) for μ = 0.6 (dotted line). The red dotted line describes the lower threshold for S as obtained from FIT and cementation pressure data. The blue asterisk shows the stress tensor combining both stress polygons for μ = 0.6 as well as μ = 0.85 and the S values resulting from FIT and cementation pressures. b). Stress profiles: pore pressure (P ) for a groundwater level of 225 m and a pressure measurement at 3,450 m TVD result in a gradient of 9.814 MPa /km to 2.25 MPa; vertical stress (S ) for a mean density of the overburden of 2.3 g/cm ; minimal horizontal stress (S ) gradients for the FIT data in the greater Munich area (S /FIT) and h h minimal horizontal stress gradient from cementation pressures at the Sauerlach wells (S /Cem.); reservoir sections for Th1, Th2, and Th3 (black dotted and solid lines). Seithel et al. Geothermal Energy (2015) 3:4 Page 7 of 22 stresses (σ ) in deviated boreholes for effective stresses (S =S − P ). It consists in a tmin s p transformation of the principal stress tensor (S ) into the stress tensor in the geographic coordinate system (S ) using the transformation matrix R . Then, S is transformed into g s g the cylindrical borehole coordinate system (S )by R . The complete transformation b b from principal (S ) to geographic (S ) and to the borehole coordinates (S )is: s g b T T S ¼ R R S R R b b s s s b The stress state in the cylindrical borehole coordinate system (S ) is referred to an angle θ, which is rotated in clockwise direction from the bottom of hole (BOH) or top of hole (TOH) along the borehole axis, respectively. The following equations describe the cylindrical stress defined by Peska and Zoback (1995). Note the correction required for the Poisson's ratio (ν) and the pressure difference (ΔP) between mud pressure and P (Zoback 2007): σ ¼ σ −2νσðÞ −σ cos2θ−4νσ sin2θ zz 33 11 22 12 σ ¼ σ þ σ −2ðÞ σ −σ cos2θ−4σ sin2θ−ΔP θθ 11 22 11 22 12 σ ¼ 2ðÞ σ cosθ−σ sinθ θz 23 13 σ ¼ ΔP rr where σ is the stress in axial direction, σ the stress tangential to the borehole, and zz θθ σ represents the shear stress. The maximum (σ ) and minimum tangential stresses θz tmax (σ ) are defined by: tmin qﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃ ðÞ σ ¼ σ þ σ þ σ −σ þ 4σ tmax zz θθ zz θθ θz qﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃ σ ¼ σ þ σ − ðÞ σ −σ þ 4σ tmin zz θθ zz θθ θz In order to determine the angles of σ and σ , θ is varied in 0.2° steps from 0° to tmax tmin 360°. This enables us to determine the direction of the maximum of σ (θ ) tmax maxmodel and the minimum of σ (θ ). From this ‘forward modeling’, the compressive tmin minmodel σ as well as tensile stresses σ and their orientations (θ /θ ) are de- tmax tmin maxmodel minmodel termined for a given well trajectory in an arbitrary stress field (see Chapter 4.2). For the Sauerlach well trajectories, the stress distributions are calculated for the well surface and compared to stress indicators like DITFs, DENFs, and BOs. This ‘stress inversion’ is applied in Chapter 4.3. Stress interpretation approach Information on the stress distribution in the Sauerlach project area was derived from st nd rd oriented caliper measurements with EMS-GR-tools available for the 1 ,2 , and 3 sections of Sauerlach Th1 and from image measurements from FMI tools for all three wells of the reservoir section (see also Table 1). Oriented caliper logs measure mechan- ically the diameter of the borehole and link it with geographical information. The FMI tool measures the resistivity of the borehole wall in high resolution and display an ori- ented borehole resistivity wall image on this basis. In image log data, natural fractures can be detected by their fracture resistivity con- trasting to the surroundings and sinusoidal trace. Completely sinusoidal features with Seithel et al. Geothermal Energy (2015) 3:4 Page 8 of 22 low resistivity (open fractures filled with mud) are distinguished from sinusoidal fea- tures, which only show a low resistivity in some areas (partly open fracture filled with mud), and sinusoidal features of higher resistivity than the formation (closed/healed fracture) (Trice 1999). Since the distribution of fracture orientation identified from the image log depends on the angle (d°) between the fracture plane and the well trajectory, a bias correction for the fracture density, called ‘Terzaghi correction’, is applied (Terza- ghi 1965). As a threshold for bias correction, it is used a value of d° <15°. BO occurs when σ exceeds the rock mass strength. They occur in the direction of tmax maximum tangential stress (θ ) (Bell and Gough 1979). These structures can be de- max tected by caliper or image logs or borehole televiewer tools. For the present study, the ‘SAC’ 6-arm caliper interpretation software was used (Wagner et al. 2004) and the cri- teria given by Reinecker et al. (2003) for 4-arm caliper tools were adapted to the 6-arm caliper tools. Averaging was performed using the circular statistics of Mardia (1972). BO analyses in image log data were described as blurry, conductive, symmetrical fea- tures (Figure 3a) by Bell (1996a). During the process of drilling in highly fractured reservoirs, stress rearrangement at the borehole wall can enhance existing fractures in the direction of σ (θ ). Barton tmin min and Moos (2010) call these fractures drilling-enhanced natural fractures (DENFs). These structures are relatively short segments of fine-aperture natural fractures, which are not visible as complete sinusoidal traces (Figure 3c). DENFs are aligned to existing sinusoidal traces and can be fitted by flexible sinusoidal traces (Barton and Moos 2010). Due to the interaction of fracture planes and tensile stress, it is sometimes diffi- cult to precisely determine the orientation of the tensile areas with DENFs. In compari- son to DENFs, DITFs occur as axisymmetric pairs parallel to the borehole axis. They are formed when σ reaches the tensile strength of the rock mass (Barton and Moos tmin Figure 3 Typical stress indicators found in FMI logs. a). Axially oriented breakouts with blurry appearance in the directions of 90° and 270° TOH. b). Drilling-induced tensile fractures (DITFs) occurring as axial fractures in the directions of 180° and 355° TOH parallel to one principal stress. c). Drilling-enhanced natural fractures (DENF) appear as inclined non-continuous fractures (green line) next to continuous natural fractures (blue sinusoidal curve) in the directions of 190° and 15° TOH. Anderson (1951) shows the alignment of these features to existing natural fractures. Bedding planes (green dotted sinusoidal curve) crosscutting the fracture planes are unaffected by stress-induced borehole wall features. Seithel et al. Geothermal Energy (2015) 3:4 Page 9 of 22 2010) and are promoted by cooling down the well or increasing the mud weight (Davatzes and Hickman 2010). For wells that are inclined to one principal stress axis, DITF occurs as en-echelon sets of fractures (E-DITFs) at a small angle to the borehole axis (Peska and Zoback 1995). Thus, axial fractures or E-DITF in combination with knowledge of the well trajectory indicate the principal stress direction. Petal-centerline fractures are features which develop during the drilling process in front of the drill bit in the plane of S due to enhanced tensile stress. They can propagate into DITFs. The difference is that they develop as complex, non-planar features that typically are not co- axially oriented to the borehole (Davatzes and Hickman 2010). These features might be but not used in the dataset and can be responsible for data scattering of drilling- induced features. In this study, we focus on DENFs and DITFs identified in the Sauerlach wells to analyze the orientation of σ (θ ) and BO in order to determine the orientation of tmin min σ (θ ) in the borehole coordinate system. These data observed are finally com- tmax max pared to the modeled stresses as described above (Chapter 3.1) and discussed. Stress regime The orientation of the identified BO as derived from the interpretation of the 6-arm caliper measurements in Sauerlach Th1 is presented in Figure 4. From the top to 500 m TVD, numerous washout zones can be seen, which limit the interpretation and cannot be included in the results. From 500 m TVD to 2,600 m TVD, clearly identifi- able BO zones exist. An inclination built-up beginning from 2,400 m TVD to more than 20° from vertical, however, makes further BO interpretation due to severe tool decentralization impossible (Figure 4). The 220 BOs identified with a total length of 625 m indicate an S - direction of approximately N-S (N 7.5°E ± 9°). The section where image logs are acquired ranges from a depth of 3,700 to 4,100 m TVD. This open-hole section is highly deviated and covered by mainly Cenozoic com- pacted sand and clay stones. For this coverage, no density measurements were per- formed. Thus, for these sedimentary layers, we assume a mean density of 2,300 kg/m (Fjaer et al. 2008). At reservoir level (4,000 m TVD), S will be 92 MPa with a stress gradient of 23 MPa/km (cf. Figure 2b). Due to low topography and investigations of focal mechanisms (Megies and Wassermann 2014), the vertical stress can be considered as a principal stress at the Sauerlach site. Pore pressure measurements were conducted after drilling at 3,450 m TVD and re- veal a static P of 31.6 MPa 550 m above the reservoir and a groundwater level of 225 m below surface. To extrapolate P to reservoir depth (4,000 m TVD), a low min- eralized model water with ρ = 1.000227 g/cm (10°C, 1 bar) is used. So the P can be extrapolated to 37 MPa (0.4 S ) at a depth of 4,000 m TVD, which results in a gradient of 9.814 MPa/km - 2.25 MPa (Figure 2b). Technical drilling operations like FIT and cementation pressures from the Sauerlach drillings help to estimate the lower boundary of S within the Molasse sediments. In Figure 2b, test pressures are illustrated as S /FIT and S /Cem. They show a mean S h h h gradient of about 15.5 MPa/km (Figure 2b). S is determined by applying the stress polygon method (Zoback et al. 2003). For a given depth (here, we use the mean reservoir depth of 4,000 m TVD), the gradients of S ,P (introduced in the previous paragraphs), and the static friction coefficient (μ) give v p Seithel et al. Geothermal Energy (2015) 3:4 Page 10 of 22 Figure 4 Breakout data measured by oriented 6-arm caliper tools in the 1st, 2nd, and 3rd sections of well Th1. Column a). Well trajectory data: the 1st and 2nd sections are nearly vertical; in the 3rd section, well inclination is increased to 45° with an azimuth of 230° to 290°. Column b). Caliper diameter (black lines) and bitsize diameter (red line): the bitsize decreases from 23″ in the 1st section to 16″ in the 2nd section and 12″ in the 3rd section. Column c). Breakout azimuth (0° to 180°) determined with the ‘SAC’ software (Wagner et al. 2004): increasing well inclination at 2,700 m MD makes breakout interpretation impossible because of tool decentralization. Column d). Stratigraphic units: Upper Freshwater Molasse (OSM) at the middle Miocene to upper Lower Miocene (50 to 500 m) is characterized by numbers of washouts; Upper Marine Molasse (UMM) at the middle Lower Miocene (500 to 800 m), good well stability and less breakouts are observed; Lower Freshwater Molasse (USM) at the middle Lower Miocene to the Lower Oligocene (800 to 2,800 m), clear breakout orientation, high number of breakouts. the frictional limits for S as well as S for each stress regime (NF, SS, TF) (Figure 2a). H h According to Jaeger et al. (2007), the stress limits for optimally oriented faults can be described by the static friction coefficient, as it is exemplarily done for μ = 0.6 (lower static friction coefficient according to Byerlee). Seithel et al. Geothermal Energy (2015) 3:4 Page 11 of 22 S −P 1 p 1= 2 2 ¼ μ þ 1 þ μ P ¼ 0:4S ; μ ¼ 0:6 p v S −P 3 p S þ 0:85 S ¼ 3:12 S 1 v 3 For normal faulting (S >S >S ), strike-slip faulting (S >S >S ), and thrust faulting v H h H v h (S >S >S ), the limits can hence be defined as: H h v S ¼ 0:58 S ðÞ S ¼ S ; S ¼ S ðÞ NF h v 1 v 3 h S ¼ 3:12 S −0:85 SðÞ S ¼ S ; S ¼ S ðÞ SS H h v 1 H 3 h S ¼ 2:27 S ðÞ S ¼ S ; S ¼ S ðÞ TF H v 1 H 3 v In Figure 2a, the stress polygon for a static friction coefficient of 0.6 (black dotted line) and for a static friction coefficient of 0.85 (black solid line) is illustrated. This method defines values for S (0.52 S for μ = 0.85; 0.58 S for μ = 0.6) which are smaller h v v than the values measured with FIT and cementation pressures (0.65 S ) (red dotted line). Using μ = 0.85, the stress polygon shows higher possible magnitudes of S ,if μ = 0.6 lower values of S are obtained for the same S /S . For the S magnitude H h V h analyzed based on FIT and cementation reports (S = 0.65 S ), S magnitudes between h v H 1.25 S (for μ = 0.6) and 1.5 S (for μ = 0.85) are obtained. Several authors give insights v v into the stress regime in the Bavarian Molasse Basin and describe it as being dominated by strike slip or thrust faulting (Illies et al. 1981; Reinecker et al. 2010), whereas local stress anomalies/perturbations (see Riem Th2 in Reinecker et al. (2010)) indicate simi- lar magnitudes of S and S . H h For the following analyses, the stress regime is defined at the intersection point of the lower boundary of S and the stress polygon for μ = 0.6 (blue asterisk in Figure 2a). The stress gradient for S is 15.5 MPa/km, for S 23 MPa/km, and for S 28 MPa/km. With h v H regard to the stress orientation, BO data confirm the pervasive N-S orientation of S (Reinecker et al. 2010). Identified natural fractures and drilling-induced stress indicators For the natural fractures analyzed according to the methods described in Chapter 3.2, a Terzaghi correction was applied. The Rose diagrams of the fracture orientations are il- lustrated in Figure 5 and statistical data presented in Table 2. The results of Sauerlach Th2 and Th3 show a bi-directional distribution in a main fracture set directed towards the ENE-WSW and a secondary orientation towards the N-S to NNE-SSW. In contrast to this, Sauerlach Th1 shows a singular natural fracture set directed N-S to NNE-SSW. Both fracture sets are inclined by 70° to 90°. In Sauerlach Th1 oriented to N 300°E (see Figure 5a black arrow), the N-S fracture set exhibits an angle of d° = 60° and the ENE- WSW fracture set has an angle of d° = 30°. Therefore, the bias correction could be ap- plied. In Sauerlach Th2 and Th3 oriented to N 0°E and to N 160°E (see Figure 5b, and c black arrow), by contrast, the angle of the NNE-SSW fracture set is d° = 0° to 15°, while that of the ENE-WSW fracture set is d° = 80°. In these wells, bias correction can be applied for the ENE-WSW fracture set but due to d° values below the critical limit of 15°, this correction cannot be applied to the NNE-SSW fracture set. Seithel et al. Geothermal Energy (2015) 3:4 Page 12 of 22 Figure 5 Data obtained from interpreting the FMI logs of Sauerlach Th1 (a), Th2 (b), and Th3 (c). First column: Rose display of the Terzaghi-corrected fractures related to the geographic coordinate system. Light blue colors indicate partly conductive/partly open fractures, dark blue colors conductive/open fractures, and red colors resistive/healed fractures (Trice 1999). The black arrow indicates the mean well orientation in the logged interval and N the total number of fractures. Second column: Drilling-induced tensile fractures (DITF) in the borehole coordinate system relative to the top of hole (TOH). N indicates the number and θ the mean orientation with standard deviation. Third column: Drilling-enhanced natural fractures (DENF). Fourth column: Breakout (BO) data. All stress indicator data are measured in the borehole coordinate system against TOH due to high well inclination in the reservoir sections. It has to be mentioned here that fracture interpretation in Sauerlach Th1 is of limited quality between 4,000 m and 4,390 m MD due to extensive intersection of different fracture orientations, which is called crossfracturing (Figure 6a). This implies a partly low fracture density, although a highly fractured or deformed reservoir can be observed in the image log. In the transition zone of the alternating strata of dolomite and lime- stone to a compact limestone in the mass facies (4,350 m MD), the fracture density in- creases and cross-fracturing disappears. Such a change in fracture density is frequently observed in the Malm section of geothermal wells. This is also observed in Sauerlach Th2 and Th3 in the transition area to the compact deep dolomite units in the mass fa- cies (4,600 m MD - Sauerlach Th2; 5,000 m MD - Sauerlach Th3). Seithel et al. Geothermal Energy (2015) 3:4 Page 13 of 22 Table 2 FMI interpretation data from the Sauerlach Th1, Th2, and Th3 Sauerlach Th1 Sauerlach Th2 Sauerlach Th3 Logged length (m) 776 825 1,147 Fractures Number 1,660 961 2,460 Orientation (north) N-S, NNE-SSW ENE-WSW; NNE-SSW ENE-WSW; NNE-SSW DITF Number 70 823 158 Total length (m) 21 207 61 Orientation (TOH) 38.5° ± 16.5° 176.3° ± 9.6° 1.52° ± 15.7° RENF Number 378 492 621 Total length (m) 66 105.5 130 Orientation (TOH) 27.2° ± 11.2° 173.9° ± 10.4° 179.9° ± 15.3° BO Number 212 80 7 Total length (m) 95 24 2.6 Orientation (TOH) 113.9° ± 14.4° 92.0° ± 14.4° 104,8° ± 12.9° Logged length describes the reservoir depth where the image log is measured. The number and orientation for picked fractures, as well as the number, total length, and orientation for drilling-induced fractures (DITF), reopen natural fractures (RENF), and breakouts (BO), are measured. BO indicating drilling-induced stress occurs in Sauerlach Th1 with a mean azimuth of 115° ± 14° to TOH in the upper part of the reservoir between 4,050 m and 4,390 m MD depth, where crossfracturing dominates (Figure 6). Deeper than 4,390 m MD, the drilling process leads to a reactivation of intense fracturing and DENF occur with a dir- ection of 27° ± 11° TOH (Figure 5a). In Sauerlach Th2 the direction of DITFs and DENFs is 176° ± 9.6° TOH/173.9° ± 10.4°. Additionally, BOs can be found in a short area between 4,020 m and 4,140 m MD with an orientation of 95° ± 14 TOH (Figure 6). In Sauerlach Th3 mainly DENFs occur over the whole reservoir interval with an orien- tation of 179.9° ± 15.3° TOH (Figure 6). Interpretation and discussion of results Fracture system and tectonic environment Fractures connected to normal faults form networks which are often subparallel to fault zones, resulting in a WSW-ENE-oriented fracture set in the Molasse Basin. Indeed, the frac- tures in the Sauerlach wells Th2 and Th3 show this orientation (Figure 5b/c). Sauerlach Th1, by contrast, is oriented between two terminating normal faults and shows a unimodal fracture set distribution between N 350° and 10°E (Figure 5a), although fractures oriented towards the WSW-ENE would be partly affected by observation bias. Near to fault zones fracture density frequently increases and rock mass strength is re- duced. In general, fault zones possess a complex structure so that it is hard to identify them solely from borehole data. Additional indication of a fault may be found in the upper part of Sauerlach Th1 (4,000 to 4,390 m MD). It is dominated by cross- fracturing which is difficult to identify, but shows a deformed rock mass. The lower part (4,390 m to 4,760 m MD) does not exhibit any deformation and only tension frac- tures occur. Typically, the Malm carbonates show a high rock mass strength reflected by a good well stability in Sauerlach Th2 and Th3 (Figure 6b/c), but cross-fracturing and the high breakout frequency between 4,000 m to 4,390 m MD in Sauerlach Th1 in- dicate cataclastic rocks and a reduced rock strength (Figure 6a). Seithel et al. Geothermal Energy (2015) 3:4 Page 14 of 22 Figure 6 Comparison of the stress indicators and the modeled tangential stress. First column: Well inclination (0° to 90°) relative to the vertical (black line), well azimuth (0° to 360°)(red line). Second column: Corrected fracture density calculated in 5-m steps from 0 to 15/m and interpreted cross-fracturing in well Th1. Third column: Stress indicators (blue - drilling-induced tensile fractures, green - drilling-enhanced natural fractures, red - breakouts) and the direction of θ (red line) and θ (blue line) for an N-S-oriented maxmodel minmodel strike-slip stress regime. All angles are measured in the borehole coordinate system with the θ angle (0° to 360°) relative to TOH. Fourth column: Facies interpretation with mass and bedded facies. The boxes show the mean direction of the observed (θ and θ ) and modeled stress indicators (θ and θ ) and the max min maxmodel minmodel difference between the observed and modeled orientations (Δθ). These data indicate that an N-S-oriented strike-slip stress regime predominates in the wells Th2 and Th3 but not in well Th1. Since no further direct information on the structural setting around Sauerlach Th1, e.g. from 3D seismic data, is available, the fracture system as a single source of information could indicate an N-S-oriented fault zone. However, an N-S-oriented fault system has not yet been discussed for the Bavarian Molasse Basin. When considering strike-slip envi- ronments (Lueschen et al. 2014; Megies and Wassermann 2014) with faults oriented to Seithel et al. Geothermal Energy (2015) 3:4 Page 15 of 22 N 20° to 45°E, however, such a fault system may exist, since fracture in strike-slip fault zones are not always oriented parallel to the fault zone (Sylvester 1988). Stress indicators and modeled stress directions The wellbore trajectory and the stress field control the state of stress at the wellbore wall and the occurrence of BOs, DENFs, and DITFs (see Chapter 3.2). The input stress magnitudes of our model approach are shown in Figure 2b (see Chapter 3.3). Figure 6 presents the results of θ (blue line) and θ (red line) with the picked minmodel maxmodel BOs, DENFs, and DITFs. Moreover, Figure 6a shows a mismatch between the orientation of observed features and those assuming a strike-slip regime. The mismatch differs in two distinct sections: (i) between 4,000 m and 4,380 m MD, there is a significant offset between θ maxmodel and θ of 35° to 45° (Figure 6/Sauerlach Th1) and (ii) between 4,380 m and 4,760 m max MD the offset is shifted by a mean of 15° to 30° with respect to TOH. In contrast to this, the observed θ from BOs and θ from DENFs and DITF in max min Sauerlach Th2 correlate well with the modeled stress orientations and confirm the stress conditions assumed for modeling (Figure 6/Sauerlach Th2). The data set of Sauerlach Th2 is considered to be representative due to the low standard deviation and high number of stress indicators in the well interval (BO, DENF, and DITF). Modeled stress orientations in Sauerlach Th3 generally confirm the assumed stress regime with DITF and DENF (Figure 6/Sauerlach Th3); however, few BOs to localize θ can be max interpreted. The observations made, in combination with modeling, suggest local stress changes in the vicinity of the wells. The wells oriented north (Sauerlach Th2) and south (Sauerlach Th3) exhibit a strike-slip stress regime with S oriented N-S. The modeled stress distribu- tion of Sauerlach Th1 oriented NW-SE differs. The mismatch between observed and mod- eled data requires a critical review of potentially influencing factors and will be discussed in the following sections. Sensitivity study of Sauerlach Th1 and Th2 Heterogeneous structures, such as highly fractured zones or locally changed lithological composition in the reservoir, can have an influence on the local stress regime and stress trajectories (Bell 1996b; Homberg et al. 1997). Depending on the distance to the fault or its properties, the direction and the magnitude of the maximum horizontal stress can sig- nificantly vary in the vicinity of these structures (Shamir and Zoback 1992; Hickman and Zoback 2004). The direction and the stress regime based on the stress indicators are investigated to identify the local stress field in the vicinity of the wells. The depth (TVD) and the well trajectory at the position of every stress indicator are included, and the angle θ maxmodel as well as θ (see Chapter 3.1) are calculated. For every stress regime, the average minmodel values of the square roots of the directional differences of the observed (θ ) vs. mod- obs eled stresses (θ ) are calculated. For every stress regime, the mean directional differ- model ence (Δθ), including all stress indicators, is defined. sﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃﬃ jj θ −θ obs model Δθ ¼ N Seithel et al. Geothermal Energy (2015) 3:4 Page 16 of 22 This value can be used to better describe the validity of a stress model applied to stress indicators. First, the influence of the S magnitude is studied by varying its magnitude. In a second step, the magnitude and the orientation of S are varied and S is kept constant. To study the influence of S , the magnitudes of S and S are kept constant in a h H v strike-slip regime (S = 1.25 S ) (Figure 7b) and the magnitude of S is varied from low H v h values (0.65 S ) to higher values (S ). In the same way, the influence of S on stress dis- v v h tribution in a normal faulting regime is verified (S = 0.9 S ;S < 0.9 S ) (Figure 7a). H v h v Figure 7 illustrates the Δθ values for S /S ratios at θ and θ . In a strike-slip re- h v max min gime, the standard deviations for different S values are between 0.8° to 6° (Figure 7b); in a normal faulting regime (Figure 7a), these values are between 0.5° and 4° and emphasize that the S magnitude does not have a significant impact on stress distribu- tion at the borehole wall. Thus, in the following consideration, S will be kept constant to simplify the analyses. In the next step, the influence of the S direction as well as of the S magnitude on H H Δθ values for compressive stress indicators is analyzed. Here, the focus lies on the wells Th1 and Th2 because of the large numbers of stress indicators. As illustrated in Figure 7, S has a small influence on Δθ. It is concluded that it remains unaffected by stress perturbations and is kept constant at 0.65 S . The direction of S is rotated v H clockwise from N 0° to 180°E (north to south), and the stress regime is varied from a normal faulting regime with nearly isotropic horizontal stress (S =0.7 S ) H v to a strike-slip faulting regime (S =1.25 S ). This variation can be traced in the H v stress polygon from the blue asterisk to the lower limit of the red dotted line (Figure 2a). Figure 7 Mean difference between the modeled and observed stress orientations (Δθ) for varied minimal horizontal stresses. θ from DITF and DENF as well as θ from BO are compared to the min max modeled directions (θ , θ )in (a). a normal faulting (S = 0.9 S ) and (b). a strike-slip stress maxmodel minmodel H v regime (S = 1.25 S ) for a variable S magnitude. The value of Δθ is calculated to evaluate the consistency H v h of the model (see Chapter 4.3). The S /S ratio's independence of the stress distribution results in constant h v Δθ values and low standard deviations. In the normal faulting regime (a)., the Δθ values for Th1, Th2, and Th3 are particularly constant for different S /S ratios. Such constant Δθ values are also seen in the h v strike-slip regime (b). for Th2 and Th3. The data of Th1 in b.) illustrate a change for increasing S /S ratios, h v but the standard deviation is even lower than 6°. These results show that the S /S ratio does not have a h v significant impact in stress distribution in the well. So, further analyses will be based on fixed S /S ratios h v with modified SH magnitude and SH orientation. Seithel et al. Geothermal Energy (2015) 3:4 Page 17 of 22 Figure 8 shows the Δθ values for compressive stress indicators (BO) in Th1 and Th2 as a function of S orientation and S magnitude. White areas define stress regimes H H with Δθ values lower than 20°. Gray to black areas indicate high Δθ values of up to 70°. Consequently, the white areas represent stress regimes with the best agreement of ob- served and modeled data. In Sauerlach Th1, the observations do not correlate with the assumed strike-slip stress regime which is marked by the blue asterisk for low values of Δθ.The best fit is marked by the black line (Figure 8/Sauerlach Th1). A good fit between model and data is obtained for either a rotation of 40° to a direction of N 150°E in a strike-slip regime (S >S ) (Figure 8/orange asterisk) or for unperturbed stress orientations H V (S = N 10°E) by a reduction of S The latter indicates a stress regime change to Hdir H. an intermediate strike-slip/normal faulting regime (S =S ) or a normal faulting re- H v gime with S =0.85 S (Figure 8/yellow asterisk). H v Figure 8 Differences of the modeled and observed stress indicators for a varied direction and stress regime. The vertical (S ) and the minimum horizontal stresses (S ) are kept constant (0.65 S ), and the stress v h v direction and the magnitude of the maximum horizontal stress (S ) are changed (S -dir: 0° to 180°, S :0.7 to H H H 1.25 S ). The location of the strike-slip stress tensor is shown by the blue asterisk. The colors illustrate the mean difference (Δθ) between the observed and modeled maximum tangential stresses. Th1 shows a rotation of the stress regime (orange asterisk), or a change in the stress regime to normal faulting (yellow asterisk) fits our observation best. In Th2, the N-S-oriented strike-slip regime is confirmed, but this analysis also suggests that stress direction is independent in a normal faulting regime. Seithel et al. Geothermal Energy (2015) 3:4 Page 18 of 22 The best fit for well Sauerlach Th2 shows that without stress rotation (S =N 10°E), Hdir the observed stress indicators can be modeled in a strike-slip regime (1.25 S >S >S ), v H v intermediate strike-slip/normal faulting regime (S =S ) or a normal faulting regime H v (S >S ) (Figure 8/Sauerlach Th2). The line of best fit correlates well with the N-S- v H oriented strike-slip stress regime marked by the blue asterisk in Figure 2. It is worth mentioning that in Sauerlach Th2 for S /S < 0.8 (normal faulting regime), the results H v are independent of the S orientation. Therefore, it is not possible to deduce a prefer- ential stress orientation. Stress modeling indicates that stress indicators observed in Sauerlach Th1 show the influence of local stress perturbations. Stress magnitude determination is based on the existence and the orientation of breakouts, but we did not use e.g. breakout depths and width in combination with strength to estimate the stress magnitude. For additional analyses, geomechanical data and logging information, such as acoustic logs, cored sec- tions, and extended pressure tests, are required. Conceptual model at the Sauerlach site According to Illies et al. (1981) and Reinecker et al. (2010), the stress regime in the Bavarian Molasse Basin is primarily influenced by the alpine orogeny and oriented accordingly per- pendicular to it (N-S to NW-SE). Breakout interpretation in the Sauerlach wells based on the caliper data in the Cenozoic sediments from 500 to 2,600 m TVD reveals a mean S direction of N 7.5°E ±9°. The total breakout length (625 m) and total logged interval with breakout zones (2,100 m) give a ratio of 0.30 and indicate a mean differential stress for breakout formation. This phenomenon is interpreted by Illies et al. (1981) and Reinecker et al. (2010) to be an indicator of a strike-slip or thrust faulting regime and also the focal mechanisms analysis (Megies and Wassermann 2014) support this interpretation. The number and total length of breakouts within the Cenozoic sedimentary layers in Sauerlach Th1 (Figure 4) is much higher than in the Malm reservoir in Sauerlach Th1, Th2, and Th3 (Table 2). This can be the result of lower tangential compressive stress due to the well trajectory, a changed stress regime or higher rock strength in the Malm reser- voir. Nevertheless, the analysis of stress indicator orientation indicates an N-S-oriented stress field with a strike-slip character for the wells Sauerlach Th2 and Th3. According to our stress inversion, a transitional strike-slip/normal faulting (S =S ) or a normal faulting v H regime (S <S ) might exist when judging from the stress data. We therefore conclude a H V strike-slip to normal faulting regime where the compressional character decreases with depth due to higher vertical stress gradients in comparison to the horizontal stress gradi- ents. The absence of stress rotation with depth points to significantly different horizontal stress magnitudes as also observed by Reinecker et al. (2010). Our stress inversion in Sauerlach Th1 indicates either a normal faulting regime (S >S ) v H with S oriented to N 10°E or a rotation of S oriented to N 150°E in a strike-slip stress H H regime (S >S ). Compared to the observations in Sauerlach Th2 and Th3, this H v would indicate a regionally perturbed stress field of 2nd- or 3rd-order stress pat- terns. Such a perturbed stress field at the Sauerlach site is not a unique phenomenon in the Molasse Basin; according to Reinecker et al. (2010), at least ten wells indicate a locally rotated stress field. Within the geothermal reservoir in the Bavarian Molasse Basin, units with a perturbed stress field can either be linked to lithofacial homogeneous bodies Seithel et al. Geothermal Energy (2015) 3:4 Page 19 of 22 (e.g. reef complexes) with a significant size or tectonic structures, as e.g. observed in a set of combined strike-slip/extensional fault system with several en-echelon normal faults building several small relay ramps (Lueschen et al. 2014). Consequently, wells drilled into the normal faults, such as Sauerlach Th2 and Th3, may exhibit a fracture set that differs significantly from that of Sauerlach Th1 which is drilled in-between two terminating nor- mal faults. The identified mainly N-S- to NNE-SSW-oriented fractures in Sauerlach Th1 are ei- ther influenced from a dilatational tendency in a normal faulting regime with an N-S- oriented SH (Figure 9/white arrows in Case II) or a shearing character in the rotated strike-slip stress regime (Figure 9/shear movement in Case I). In both scenarios, the fractures would have an opening character and provide a good hydraulic activity. Therefore, the above-presented structural concept goes along with the observations during the early stage of pumping tests corresponding with a high negative skin effect indicating a good connectivity of Sauerlach Th1 to the reservoir (see Chapter 2.1). From this perspective, it is presumed that these open fractures contribute to the high productivity of the well by a low turbulent inflow into the borehole. Compared to the fractures in Sauerlach Th2 and Th3 which are mainly oriented ENE-WSW, they are identified to be closing within the observed stress field. Therefore, it seems that these are unfavorable for geothermal targeting though a Figure 9 Stress field interpretation near the well Th1 and the stress pattern in Th2 and Th3. Gray arrows present the orientation of the maximum horizontal stress. Here, local variations of a strike-slip regime (SS) and a normal faulting regime (NF) can be observed. For well Th1, a local perturbation of the stress regime is observed and determined to be strike slip and S oriented to 320° (Case I) or normal faulting regime with S oriented towards the N-S (Case II). In Th2 and Th3, the strike-slip regime is observed with S having an N-S orientation. The Rose diagram illustrates the fracture sets. The principal fracture set in Th1 for stress state Case I indicates a high shear tendency. For Case II, these fracture sets exhibit a dilatational tendency (white arrows). Fracture sets in Th2 and Th3 show a dilatational tendency (white arrows) (NNE-SSW) as well as a closing tendency (black arrows) (ENE-WSW) in the strike-slip stress regime. Seithel et al. Geothermal Energy (2015) 3:4 Page 20 of 22 significant displacement along with faulting and was initially encountered in the seismic data. Conclusion The stress analysis made within the framework of this study has shown that the strike-slip stress regime with a N-S-oriented S can be found in the Cenozoic layers in one well (Sauerlach Th1) and even in the Malm reservoir in two other wells (Sauerlach Th2 and Th3). Stress inversion of the breakout data measured in the Malm reservoir of Sauerlach Th1 and Th2 indicates a perturbed stress regime for Sauerlach Th1. The fracture systems within the Malm reservoir of Sauerlach Th2 and Th3, both drilled towards the WSW- ENE faults, have a predominant WSW-ENE orientation. In contrast, Sauerlach Th1, drilled between the northern and southern WSW-ENE fault zones, exhibits an N-S- to NNE-SSW-oriented fracture system. This tectonic situation can lead to consequences for the hydraulic system for both inferred stress regimes: A) In a stress regime with nearly N-S-oriented S , the N-S- to NNE-SSW-oriented fracture system, which is dominant in the Sauerlach Th1 but very minor seen in Sauerlach Th2 and Th3, can have a dilatational character (Figure 9/white arrows in Case I). B) In a N 150°E-oriented strike-slip stress regime, the N-S- to NNE-SSW-oriented fracture planes will reach highest shear stress (Figure 9/shear movement in Case I) and most likely experience shear movement. In both cases, the fracture aperture might be increased, leading to better hydraulic con- nectivity. This is supported by the fact that Th1 has a higher productivity. Overall, the iden- tified stress effects help us to understand the hydraulic properties of the well Sauerlach Th1 compared to the wells Sauerlach Th2 and Th3. In this study, we could assess the stress distribution in highly deviated wells at the Sauerlach site. The approach chosen, however, is not site specific and can be transferred to other wells. The ability to identify and precisely locate stress perturbations can improve the predictability of hydraulic performance and thus contribute to a better assessment strategy. Competing interests The authors declare that they have no competing interests. Authors' contributions RS carried out the data analysis and modeling and wrote the manuscript. US and BM participated in analyzing the results and in drafting the manuscript. TK initiated the study and participated in the design, coordination, and quality control of the project and in drafting the manuscript. CH completed the data set and was involved in drafting the manuscript. All authors read and approved the final manuscript. Acknowledgements Grateful thank goes to the reviewers for their helpful expertice in improving the manuscript. We would like to thank the Stadtwerke München for the provision of the data set and the permission to publish the results. Special thanks go to Erdwerk GmbH for the possibility to carry out such a project and the support during data collection. Additionally, we thank the EnBW AG and especially the Department of Research and Development for the support of the Chair for Geothermal Research of KIT. Author details 1 2 Karlsruher Institute of Technology (KIT), Adenauerring 20b, 76131 Karlsruhe, Germany. Erdwerk GmbH, Bonner Platz 1, 80803 München, Germany. Stadtwerke München GmbH (SWM), Emmy-Noether-Straße 2, 80992 München, Germany. Received: 15 June 2014 Accepted: 11 December 2014 Seithel et al. Geothermal Energy (2015) 3:4 Page 21 of 22 References Anderson EM. The dynamic of faulting and dyke formation with application to Britain. Edinburgh: Oliver and Boyd; 1951. Bachmann GH, Mueller M, Weggen K. Evolution of the Molasse Basin (Germany, Switzerland). Tectonophysics. 1987;137(1–4):77–92. doi:10.1016/0040-1951(87)90315-5. Barton C, Moos M. Geomechanical wellbore imaging: key to managing the asset life cycle. In: Poppelreiter MG-CC, Kraaijveld M, editors. Dipmeter and Borehole Image Log Technology. AAPG Memoir, vol 92. 2010. p 81–112. doi:10.1306/13181279M922689. Barton CA, Zoback MD, Moos D. Fluid-flow along potentially active faults in crystalline rock. Geology. 1995;23(8):683–6. doi:10.1130/0091-7613(1995)023 < 0683:Ffapaf > 2.3.Co;2. Bell JS. Petro Geoscience 1. In situ stresses in sedimentary rocks (part 1): measurement techniques. Geosci Can. 1996a;23(2):85–100. doi:10.12789/gs.v23i2.3902. Bell JS. Petro Geoscience 2. In situ stresses in sedimentary rocks (part 2): applications of stress measurements. Geosci Can. 1996b;23(3):135–53. doi:10.12789/gs.v23i3.3910. Bell JS, Gough DI. Northeast-Southwest compressive stress in Alberta - evidence from oil wells. Earth Planet Sc Lett. 1979;45(2):475–82. doi:10.1016/0012-821x(79)90146-8. Betz D, Wendt A. Neuere Ergebnisse der Aufschluss- und Gewinnungstätigkeiten auf Erdöl und Erdgas in Sueddeutschland. Bulletin Verein Schweizer Petroleum-Geologen und -Ingenieure. 1983;49(117):9–36. Boehm F, Savvatis A, Steiner U, Schneider M, Koch R. Lithofacies and characterization of the geothermal Malm reservoir in the greater area of Munich. Grundwasser. 2012;18(1):3–13. doi:10.1007/s00767-012-0202-4. Byerlee J. Friction of rocks. Pure Appl Geophys. 1978;116(4–5):615–26. doi:10.1007/Bf00876528. Davatzes NC, Hickman SH. Stress, fracture, and fluid-flow analysis using acoustic and electrical image logs in hot fractured granites of the Coso Geothermal Field, California, U.S.A. In: Pöppelreiter M, García-Carballido C, Kraaijveld M, editors. Dipmeter and borehole image log technology. AAPG Memoir, vol 92. 2010. p 259–293. doi:10.1306/13181288M923134. Dorsch K, Pletl C. Bayerisches Molassebecken - Erfolgsregion der Tiefengeothermie in Mitteleuropa. Geothermische Energie. 2012;73:14–8. Faulds JE, Hinz NH, Coolbaugh MF, Cashman PH, Kratt C, Derling G, et al. Assessment of favorable structural settings of geothermal systems in the Great Basin. San Diego, California, USA: Western USA. Paper presented at the Geothermal Resources Council Transactions; 2011. Fjaer E, Holt RM, Horsrud P, Raaen AM, Risnes R. Petroleum related rock mechanics, vol. 2. Amsterdam, Netherlands: Elsevier; 2008. GTV (2014) Liste der tiefen Geothermieprojekte in Deutschland. www.geothermie.de/fileadmin/useruploads/wissenswelt/ Projekte/Projektliste_Tiefe_Geothermie_2014_alphabetisch.pdf. Accessed 12.06 2014. Hickman SH, Zoback MD. Stress orientation and magnitudes in the SAFOD pilot hole. Geophys Res Lett. 2004;31. doi:10.1029/2004GL020043. Hiramatu Y, Oka Y. Determination of stress in rock unaffected by boreholes or drifts from measured strains or deformations. Int J Rock Mech Min. 1968;5(4):337–53. doi:10.1016/0148-9062(68)90005-3. Homberg C, Hu JC, Angelier J, Bergerat F, Lacombe O. Characterization of stress perturbations near major fault zones: insights from 2-D distinct-element numerical modelling and field studies (Jura Mountains). J Struct Geol. 1997;19(5):703–18. doi:10.1016/S0191-8141(96)00104-6. Hubbert MK, Willis DG. Mechanics of hydraulic fracturing. AIME: Petroleum Transaction; 1972. p. 210. Illies JH, Baumann H, Hoffers B. Stress pattern and strain release in the Alpine Foreland. Tectonophysics. 1981;71(1–4):157–72. doi:10.1016/0040-1951(81)90059-7. Jaeger J, Cook NG, Zimmerman R. Fundamentals of rock mechanics. Malden: Blackwell Pub.; 2007. Kirsch G. Die Theorie der Elastizität und die Bedürfnisse der Festigkeitslehre. Zeitschrift des Vereins deutscher Ingenieure. 1898;42:797–807. Kraemer C. Identifikation horizontaler Blattverschiebungen im Untergrund des Molassetrogs - neue Interpretationsergebnisse auf Grundlage von 2D-Seismik. Paper presented at the Geothermiekongress. Germany: Bochum; 2009. Kuhlemann J, Kempf O. Post-Eocene evolution of the North Alpine Foreland Basin and its response to Alpine tectonics. Sediment Geol. 2002;152(1–2):45–78. doi:10.1016/S0037-0738(01)00285-8. Lemcke K. Geologie von Bayern - Das bayerische Alpenvorland vor der Eiszeit. Stuttgart: E. Schweizerbart' sche Verlagsbuchhandlung; 1988. Lueschen E, Wolfgramm M, Fritzer T, Dussel M, Thomas R, Schulz R. 3D seismic survey explores geothermal targets for reservoir characterization at Unterhaching, Munich, Germany. Geothermics. 2014;50:167–79. doi:10.1016/j.geothermics.2013.09.007. Mardia KV. Statistics of directional data: probability and mathematical statistics. London: London Academic Press; 1972. Mastin L. Effect of borehole deviation on breakout orientation. J Geophys Res. 1988;93(B3):9187–95. doi:10.1029/ JB093iB08p09187. Megies T, Wassermann J. Microseismicity observed at a non-pressure-stimulated geothermal power plant. Geothermics. 2014;52:36–49. doi:10.1016/j.geothermics.2014.01.002. Meyer RKF, Schmidt-Kaler H. Paläogeographischer Atlas des sueddeutschen Oberjura (Malm), vol. 115. Hannover: Geologisches Jahrbuch; 1989. Mueller M, Nieberding F, Wanninger A. Tectonic style and pressure distribution at the Northern Margin of the Alps between Lake Constance and the River Inn. Geol Rundsch. 1988;77(3):787–96. doi:10.1007/Bf01830185. Peska P, Zoback MD. Compressive and tensile failure of inclined well bores and determination of in situ stress and rock strength. J Geophys Res-Sol Ea. 1995;100(B7):12791–811. doi:10.1029/95jb00319. Pletl C, Stoyke R, Toll H. Bohrerfahrungen bei Sueddeutschlands größtem. Geothermieprojekt bbr. 2010;3:38–47. Reinecker J, Tingay M, Mueller B. Borehole breakout analysis from four-arm caliper logs. World stress analysis Guidelines. www.world-stress-map.org Reinecker J, Tingay M, Mueller B, Heidbach O. Present-day stress orientation in the Molasse Basin. Tectonophysics. 2010;482(1–4):129–38. doi:10.1016/j.tecto.2009.07.021. Savvatis A. Personal communication of well analyses. 2014. Seithel et al. Geothermal Energy (2015) 3:4 Page 22 of 22 Schulz R, Agemar T, Alten J-A, Brunken J, Heber M, Kuder J, et al. Aufbau eines geothermischen Informationssystems für Deutschland. Hannover: Leibniz-Institut für Angewandte Geophysik; 2009. Shamir G, Zoback MD. Stress orientation profile to 3.5 km depth near the San Andreas Fault at Cajon Pass, California. J Geophys Res. 1992;97(B4):5059–80. doi:10.1029/91JB02959. Sibson RH. Frictional constraints on thrust, wrench and normal faults. Nature. 1974;249:542–4. doi:10.1038/249542a0. Steiner U, Savvatis A, Boehm F, Schubert A. Explorationsstrategie tiefer geothermischer Ressourcen am Beispiel des sueddeutschen Oberjuras (Malm). In: Bauer M, Freeden W, Jacobi H, editors. Handbuch Tiefe Geothermie. Berlin Heidelberg: Springer Spektrum Verlag; 2014. p. 429–57. doi:101007/978-3-642-54511-5. Stier P, Prestel R. Der Malmkarst im süddeutschen Molassebecken - Ein hydrogeologischer Überblick. In: LfW B, LGR B, editors. Hydrogeologische Energiebilanz und Grundwasserhaushalt des Malmkarsts im süddeutschen Molassebeckens, vol. 03E-6240 A/B. München: Schlussbericht zum Forschungsvorhaben; 1991. p. 12. STWIVT. Bayerischer Geothermieatlas - Hydrothermale Ernergiegewinnung. München: Bayerisches Staatsministerium für Wirtschaft Infrastruktur Verkehr und Technologie (STWIVT); 2010. Sylvester AG. Strike-slip faults. Geol Soc Am Bull. 1988;100(11):1666–703. doi:10.1130/0016-7606(1988)100 < 1666:Ssf > 2.3.Co;2. Terzaghi RD. Sources of error in joint surveys. Geotechnique. 1965;15(3):287–304. doi:10.1680/geot.19184.108.40.2067. Tingay M, Mueller B, Reinecker J, Heidbach O. State and origin of the present-day stress field in sedimentary basins: new results from the World Stress Map Project, Paper presented at the Golden Rocks 2006. Colorado, USA: The 41st U.S. Symposium on Rock Mechanics (USRMS), 17–21 June, Golden; 2006. Trice R. A methodology for applying a non unique, morphological classification to sine wave events picked from borehole image log data. In: Lovell MAWG, Harvey PK, editors. Borehole imaging: applications and case histories, vol special publications. London: Geological Society; 1999. p. 77–90. doi:10.1144/GSL.SP.1999.159.01.04. Unger H-J. Die tektonischen Strukturen der bayerischen Ostmolasse. Documenta Naturae. 1999;124:1–16. Wagner D, Mueller B, Tingay M. Correction for tool decentralisation of oriented six-arm caliper logs for determination of contemporary tectonic stress orientation. Petrophysics. 2004;42:1049–76. Ziegler PA. Compressional intra-plate deformation in the Alpine Foreland - an introduction. Tectonophysics. 1987;137:1–5. doi:10.1016/0040-1951(87)90223-X. Zoback MD. Reservoir geomechanics. Cambridge: Cambridge University Press. 2007. Zoback MD, Barton CA, Brudy M, Castillo DA, Finkbeiner T, Grollimund BR, et al. Determination of stress orientation and magnitude in deep wells. Int J Rock Mech Min. 2003;40(7–8):1049–76. doi:10.1016/j.ijrmms.2003.07.001. Zoback ML. 1st-order and 2nd-order patterns of stress in the lithosphere - the World Stress Map Project. J Geophys Res-Sol Ea. 1992;97(B8):11703–28. doi:10.1029/92jb00132. Submit your manuscript to a journal and beneﬁ t from: 7 Convenient online submission 7 Rigorous peer review 7 Immediate publication on acceptance 7 Open access: articles freely available online 7 High visibility within the ﬁ eld 7 Retaining the copyright to your article Submit your next manuscript at 7 springeropen.com
Geothermal Energy – Springer Journals
Published: Feb 12, 2015
Access the full text.
Sign up today, get DeepDyve free for 14 days.