Next Article in Journal
Differences in Geochemical Signatures and Petrogenesis between the Van Canh and Ben Giang-Que Son Granitic Rocks in the Southern Kontum Massif, Vietnam
Previous Article in Journal
The Petrology and Geochemistry of REE-Enriched, Alkaline Volcanic Rocks of Ambitle Island, Feni Island Group, Papua New Guinea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Slip Activation Potential of Fractures in the Crystalline Basement Rocks of Kuujjuaq (Nunavik, Canada) to Assess Enhanced Geothermal Systems Development

1
INRS—Institut National de la Recherche Scientifique, 490 rue de la Couronne, Quebec, QC G1K 9A9, Canada
2
Department of Earth & Environmental Sciences, Centre for Environmental and Information Technology (EIT), 200 University Ave. W, Waterloo, ON N2L 3G1, Canada
3
Natural Resources Canada, 1 Haanel Drive, Ottawa, ON K1A 1M1, Canada
*
Author to whom correspondence should be addressed.
Geosciences 2023, 13(11), 340; https://doi.org/10.3390/geosciences13110340
Submission received: 13 August 2023 / Revised: 29 September 2023 / Accepted: 31 October 2023 / Published: 7 November 2023

Abstract

:
This work presents an estimate of the slip activation potential of existing fractures in a remote northern community located on Canadian Shield rocks for geothermal purposes. To accomplish this objective, we analyzed outcrop analogues and recorded geometrical properties of fractures, namely the strike and dip. Then, we estimated the stress regime in the study area through an empirical approach and performed a probabilistic slip tendency analysis. This allowed us to determine the slip probability of the pre-existing fractures at the current state of stress, the orientation of fractures that are most likely to be activated and the fluid pressures needed for the slip activation of pre-existing fractures, which are key aspects for developing Enhanced Geothermal Systems. The results of this simple, yet effective, analysis suggest that at the current state of stress, the pre-existing natural fractures are relatively stable, and an injection pressure of about 12.5 MPa/km could be required to activate the most optimally oriented fractures to slip. An injection of water at this pressure gradient could open the optimally oriented pre-existing fractures and enhance the permeability of the reservoir for geothermal fluid extraction. The information described in this paper provides a significant contribution to the geothermal research underway in remote northern communities.

1. Introduction

Canada is an energy-intensive and energetically contrasting country (Figure 1; [1]). Canadians in the densely populated southern region are connected to a provincially interconnected electricity grid, powered by a variety of energy sources. In the Northwest Territories and the Yukon, most communities are connected to independent territorial grids dominated by hydro power. The situation in Nunavut and Nunavik (an administrative region of Quebec) is however quite different as there are no territorial grids, and communities generate their own electricity using generators connected to micro-grids [1,2]. Despite some renewable power projects and several initiatives to diversify the Nunavut and Nunavik communities’ energy portfolio, diesel is still their main source of energy [1,2,3].
Beyond being not connected to the provincial power grid, the communities in Nunavut and Nunavik are also not connected by road. Airplanes and boats are the only means of transportation for goods, materials, and people to reach the communities. The residents of Nunavut and Nunavik are thus compelled to provide for themselves, since the communities cannot share infrastructural assets and resources.
Another particularity of the communities in Nunavut and Nunavik is that they are located on the crystalline basement rocks of the Canadian Shield. Metamorphic and igneous rocks, Archean to Paleoproterozoic in age, dominate the geology of these regions [4,5].
An assessment of the Canadian geothermal potential carried out by Grasby et al. [6] suggested that engineered/enhanced geothermal systems (EGS) could be an option to extract the geothermal resource in the communities located on the Canadian Shield (Figure 1).
Figure 1. Distribution of geothermal potential in Canada by end use and location of power grid and remote communities (adapted from [6,7]). Provinces: BC—British Columbia; AB—Alberta; SK—Saskatchewan; MB—Manitoba; ON—Ontario; QC—Quebec; NL—Newfoundland and Labrador; NB—New Brunswick; NS—Nova Scotia; PE—Prince Edward. Territories: YT—Yukon Territory; NT—Northwest Territories; NU—Nunavut.
Figure 1. Distribution of geothermal potential in Canada by end use and location of power grid and remote communities (adapted from [6,7]). Provinces: BC—British Columbia; AB—Alberta; SK—Saskatchewan; MB—Manitoba; ON—Ontario; QC—Quebec; NL—Newfoundland and Labrador; NB—New Brunswick; NS—Nova Scotia; PE—Prince Edward. Territories: YT—Yukon Territory; NT—Northwest Territories; NU—Nunavut.
Geosciences 13 00340 g001
EGS is an emerging technology that arose from a concept initiated in Los Alamos (USA), Cornwall (UK), and Soultz-sous-Forêts (France) for exploiting geothermal resources in low-permeability rocks [8,9,10,11,12,13,14]. The central concept is to engineer fluid flow pathways between two or more wells to permit circulation of fluid through the reservoir rock mass to extract heat at rates of commercial interest.
EGS technology is needed in low-porosity fractured crystalline rocks whose natural permeability is very low due to poor hydraulic connectivity within the natural fracture network, requiring stimulation to increase the permeability of the rock mass [12]. Generally, hydraulic stimulation techniques are used for this objective.
During hydraulic stimulation, a high-pressure fluid is injected through the wellbore into the crystalline rock mass, leading to the shearing and opening of natural fractures, enhancing the hydraulic connectivity of these pre-existing fractures in the reservoir far field [11,12,13,15,16], and perhaps developing new fractures and extending existing fractures. In a successful hydraulic stimulation treatment, geothermal wells become better connected to the reservoir far field, and the geothermal resource can be exploited by circulating fluid between wells through the newly stimulated reservoir where the newly conductive fractures act as natural heat exchangers.
Permeability enhancement via hydraulic stimulation is commonly associated with two main mechanisms [15]: (1) hydraulic fracturing (opening of existing fractures and initiation and propagation of new tensile fractures), and (2) hydraulic shearing (activation of existing discontinuities in shear, leading to irreversible dilation). These two mechanisms are often referred to as Mode I and Mode II fracturing and correspond to the “end members” of the stimulation mechanisms while mixed-mode stimulation mechanisms can also occur [17,18]. This is because the pressure required for natural fractures to fail in shear may be lower than the pressure required for pure tensile fracturing (Figure 2; [19]).
Fluid injection changes the in situ pore pressure, and the effective normal stresses on the relatively weaker slip planes are reduced, diminishing the frictional strength. This leads to the activation of slip along faults and joints within the fracture network that are favorably aligned with the principal stress orientations (i.e., at a 30° ± 5° from the maximum horizontal stress [21]). Shearing and dilation of natural pre-existing fractures result in a permanent increase in the rock mass permeability to levels suitable for the extraction of geothermal resources [12,13,15,16]. A review of hydraulic stimulation experiments highlights that hydroshearing will be more effective in areas with high differential stresses and hence high shear stress (Figure 2; [22]), a condition referred to as “critically stressed” if the joint is close to the onset of the shear slip.
However, the activation of natural fractures can generate seismic events. In fact, a major issue associated with hydraulic stimulation is induced seismicity, especially in critically stressed crystalline basement rocks [23,24].
Monitoring of injection-induced seismicity is commonly used to assess how successful hydraulic stimulation treatments were and to evaluate the development of flow pathways within the reservoir. The microseismic event cloud collected as part of the monitoring program helps in designing further drilling into the reservoir [13,25].
Although induced seismicity is necessary for reservoir development, if not properly assessed, predicted, and controlled, it can become a major barrier to the development of geothermal energy [26]. Most induced events are of low magnitude (below local magnitude ML = 2) and produce ground motions far below the threshold to be felt [27]. However, large magnitude (ML > 2.5) events may be felt at the surface and lead to the suspension of an EGS project [27], such as occurred with the Basel Deep Heat Mining project in Switzerland [22,26].
A primary induced seismicity assessment tool is to test the likelihood of fracture activation (i.e., shear slip) related to stress field perturbations [28]. This can be carried out by applying the technique of slip tendency analysis [29] and Mohr diagrams to answer fundamental questions that arise at early stages of geothermal exploration [20,30]:
  • What is the probability of shear slip on pre-existing fractures at the current state of stress?
  • Which orientations of fractures are most likely to be activated?
  • What in situ fluid pressure is required to overcome the shear stress and activate pre-existing fractures?
The propensity of a surface to undergo shear slip depends on its frictional characteristics and the ratio of shear to normal effective stress acting on the surface—the slip tendency [29]. Slip tendency analysis allows us to investigate slip potential along any fracture orientation with respect to the ambient stress field, and consequently, assess the slip probability subject to a pore pressure perturbation.
Usually, a slip tendency analysis is carried out in a deterministic manner, considering a single analysis as definitive. However, subsurface conditions are uncertain and each geomechanical parameter that influences the rate and magnitude of injection-induced seismicity has inherent levels of uncertainty. Understanding this uncertainty is helpful to make well-informed decisions regarding user-controlled parameters (e.g., injection pressure, flow rate, fluid volume, etc.) [31]. Thus, a probabilistic slip tendency analysis is a relevant risk management and mitigation tool [31].
A probabilistic slip tendency analysis considers the inherent uncertainties associated with each input variable, including stress magnitudes and orientations, fault dip directions, angles, and frictional strengths. In a probabilistic slip tendency analysis, the uncertainties of the input variables (e.g., pressure, distance) are considered, allowing for a more comprehensive evaluation of slip probability in various scenarios. The slip tendency is directly influenced by the magnitude of the principal stresses, meaning that any changes in these stresses will directly impact the outcome. Therefore, a probabilistic approach provides a more thorough and suitable method for assessing the slip probability across multiple scenarios.
Minimum requirements to carry out a slip tendency analysis, either deterministic or probabilistic, are an accurate knowledge of the stress field and of the geometrical properties of fractures, namely the strike and dip.
In the absence of subsurface data, outcrop analogues are useful tools to acquire information regarding the geometrical properties of fractures [32,33,34,35]. Care is necessary when analyzing fractures from exposed outcrops, since these may be a result of exhumation and stress relief and may not persist in the deep subsurface [32,35]. Nevertheless, outcrop analogues are particularly useful in remote environments where subsurface data are unavailable, and it can be assumed that for a first-order analysis in crystalline rock masses, there are no significant differences between rocks at the surface and at depth.
The International Society for Rock Mechanics (ISRM) suggests overcoring and hydraulic fracturing and/or hydraulic testing of pre-existing fractures as methods to estimate rock stress [36,37,38]. However, in remote areas, these methods often cannot be applied, and a first-order estimate of the stress field relies on a compilation of available structural and rock mass fabric data and the application of empirical relationships. This is the first step to establish an in situ stress model [39].
The objective of this work is to estimate the slip activation potential of existing fractures in subsurface crystalline basement rocks from outcrop analogues and empirical stress predictions. The interactive computer tool developed by Yaghoubi et al. [31,40] for a probabilistic assessment of the slip tendency of faults at the Alberta No. 1 geothermal project site in Alberta, Canada, and in the Montney Formation, is used to make the necessary calculations and to plot Mohr diagrams to analyze the relationships between stresses, fluid pressure, and fractures, and provide a graphic representation of the effective stress states and slip potential.
First, we present the results of the fracture network characterization based on outcrop analogues, and the estimates of the stress field based on data compilation and empirical relationships. This information is then used in the slip tendency analysis. The remote community of Kuujjuaq in Nunavik (study area; Figure 1) is used as an example. This community was selected since geothermal research has been previously conducted at this location [41,42,43,44,45,46,47,48,49,50], and this work represents a further contribution.

2. Fracture Network Characterization

Geometrical properties of fractures (i.e., strike and dip) were sampled on exposed outcrop surfaces using the scanline sampling method (Figure 3) and a transit compass corrected for the −21°1′ magnetic declination in Kuujjuaq [51]. Since no distinction between joints and faults was made while collecting information, the term fracture is preferred and used in this work to group joints and faults.
In this method, a tape is laid on the outcrop surface perpendicular to the fracture sets observed [52,53]. Information about the strike and dip of each fracture intersecting the tape is then recorded.
Six fracture sampling areas around the community of Kuujjuaq were selected (Figure 4) considering their lithology, quality and extension of the rock exposure, and proximity to the community.
Sampling areas A1 and A2 are paragneiss outcrops, sampling areas A3 and A4 are outcrops of diorite, sampling area A5 is an outcrop in the paragneiss–tonalite contact, and sampling area A6 is an outcrop showing the contact of tonalite with gabbro.
Sampling area A1 is characterized by two main fracture sets: E-W and NW-SE. Sampling area A2 is also characterized by two main sets: E-W and NNE-SSW. Sampling area A3 is characterized by one main fracture set oriented WNW-ESE and a minor set oriented NW-SE. Sampling area A4 is characterized by one main fracture set oriented E-W and two minor sets oriented NNW-SSE and NE-SW. Sampling area A5 is characterized by two main sets, one oriented N-S and the other oriented NE-SW. In sampling area A6, three main sets were identified: NW-SE, E-W, and NE-SW. Grouping fracture information from these six sampling areas suggests four main fracture sets: F1—E-W; F2—N-S; F3—NNW-SSE; and F4—NW-SE (Figure 5).
A major issue in Kuujjuaq is that the majority of analyzed outcrops have only horizontal surfaces. An exception is outcrop A6 (Figure 6). Here, a vertical surface was available for the analysis, and the dip and dip direction of each fracture could be accurately characterized. The fractures observed suggest a high-angle inclination, with the dip ranging between 60 and 80°.

3. Empirical Stress Regime Estimates

Stress data are almost nonexistent in remote northern regions in Canada (Figure 7a). This, associated with the lack of recorded seismic events (Figure 7b), makes empirical stress predictions a useful tool for a first-order assessment of the activation potential of fractures.
The first step in assessing the stress regime through an empirical approach is to compile the available information on the orientation and magnitude of the principal stress components (Table 1). This compilation suggests a regional compression trend of NE–SW for the maximum horizontal stress in the Canadian Shield [58,59,60,61,62,63,64], except at the Raglan mining site where the maximum horizontal stress is oriented N–S [64]. The available stress data suggest a reverse faulting regime [58,60,63,64,65]. The gradient of the vertical stress is found to be about 26.0 MPa km−1 [59,60,61,62,63,64,65,66]. The gradient of the minimum horizontal stress ranges between 23.0 and 37.8 MPa km−1, while the gradient for the maximum horizontal stress is between 34.4 and 51.3 MPa km−1 [59,60,63,64,65,66].
The second step of the empirical approach is applying empirical correlations for a first-order approximation of the magnitude of the principal stresses and in situ fluid pressure.
The vertical stress component can be simply estimated based on the weight of the overlying rock at depth as [67]
S V = ρ g z
where SV (Pa) is the vertical stress, ρ (kg m−3) is the density of the geological materials, g (m s−2) is the gravitational acceleration, and z (m) is depth.
The horizontal stress components were estimated based on the horizontal to vertical stress ratio as [68]
S H m a x = k m a x S V S h m i n = k m i n S V
where k (–) is the stress ratio coefficient and SHmax and Shmin (Pa) are the maximum and minimum horizontal stresses, respectively. Stress ratio expressions obtained for the Canadian Shield (Table 2) were used in this work.
The in situ fluid pressure can be estimated based on the pore-fluid factor [69,70]:
P p = 0.4 × S V
where Pp (Pa) is the in situ fluid pressure and 0.4 is the pore to fluid factor assuming a hydrostatic regime.
The Monte Carlo-based sensitivity analysis carried out by Miranda et al. [48] to estimate the magnitude of the principal stresses suggests the following gradients: 27 ± 1.3 MPa km−1 for the vertical stress, 42 ± 5.7 MPa km−1 for the maximum horizontal stress, and 30 ± 3 MPa km−1 for the minimum horizontal stress (Table 3). The gradient for the in situ fluid pressure was estimated as 11 ± 0.5 MPa km−1 (Table 3). The value on the left side of the ± sign corresponds to the mean, while the value on the right side of the ± sign corresponds to one standard deviation. The choice of using gradient values, i.e., MPa km−1, instead of absolute values in the analysis was to make the results independent of depth.
Since each stress was estimated individually, the distributions overlap (Figure 8), and some scenarios may be unrealistic (e.g., a normal fault regime). The compilation of stress data measured in the Canadian Shield suggests a reverse faulting regime, with SHmax = S1 and SV = S3 (Table 1). This allows us to eliminate unrealistic scenarios: for example, assuming that the SV distribution is realistic, then normal faulting regimes can be eliminated. Similarly, scenarios where Shmin > SHmax arises because of the random sampling are rejected.

4. Slip Tendency Analysis and Reactivation Potential

A well-established method for assessing the fracture slip tendency under different stress conditions is using Mohr diagrams and the Mohr–Coulomb yield criterion [29,71]. The method involves calculating the shear strength of a fracture and comparing it with the maximum shear stress that the fracture can withstand before it slips [71]. The basic equation for the Mohr–Coulomb yield criterion is
τ = c + ( σ n P p ) × tan ( ϕ )
where τ (Pa) is the shear stress acting on the fracture, c (Pa) is the cohesion of the fracture, σn (Pa) is the normal stress acting on the fracture, Pp (Pa) is the pore pressure in the fracture and ϕ (°) is the angle of internal friction of the fracture. ϕ is equal to the friction coefficient (µ) as
μ = t a n ( ϕ )
For a cohesionless fracture, i.e., c = 0, the equation for the Mohr–Coulomb yield criterion becomes
τ = μ ( σ n P p )
The fracture is unstable and will slip when the resolved shear stress becomes greater than the effective normal stress multiplied by the friction coefficient acting on the fracture plane:
τ > μ ( σ n P p )
Meeting this condition generally requires the fracture to be hydraulically connected to a source of pore pressure perturbation but also mechanically unstable, potentially leading to the generation of seismic events as the slip takes place. The shear and normal stress acting on a fracture plane can be estimated using mathematical equations that take into account the stress tensor and the orientation of the fracture [71].
Figure 9 shows the 3D Mohr diagram produced considering the average stress magnitude gradients based on the empirical stress regime estimates described in the previous section (SHmax 42.1 MPa km−1, Shmin 30 MPa km−1, Sv 27 MPa km−1).
A point in the 3D Mohr diagram can represent a fracture plane with shear stress (τ) and normal stress (σn) acting on it. The stability of the plane can be assessed by comparing the shear and normal stresses to the Mohr–Coulomb yield criterion. For the case in our study, with respect to the average stress magnitude gradient, the minimum required coefficient of friction is 0.35 (Figure 9). This indicates that the ratio of shear stress (τ) to normal stress (σn) acting on the most likely fracture to slip should be less than 0.35 for the rock mass to remain stable.
This coefficient of friction is smaller than the typical coefficient of friction for joint and fault surfaces in crystalline rock; usually it ranges between 0.6 and 1.0 [72]. Considering the minimum required coefficient of friction is 0.35, most fractures are considered stable, meaning they are not prone to slip under the existing stress conditions. However, the stability of fractures is contingent upon the state of stress remaining constant. If the stress conditions change, for example through thermoelastic processes as well as pore pressure changes, the stability of fractures can be altered [73].
The pore pressure parameter in the Mohr–Coulomb yield criterion can be changed by human intervention and thus alter the state of stress surrounding the fractures. This occurs during hydraulic stimulation techniques used to increase the productivity of geothermal reservoirs. Fluid injection increases the pore pressure in the surrounding interconnected fractures, resulting in a reduction of the effective normal stress, thereby allowing stressed fractures in the right orientation to slip more easily (Equation (9)). Furthermore, hydraulic stimulation of the rock mass is intended to create new pathways for fluid flow [22], so that a pore pressure perturbation can affect more distant fractures and potentially a critically stressed fault. As the flow efficiency of the geothermal reservoir is improved by creating better fracture conductivity in the natural fracture system, more structures are involved in the pressure perturbation.
The critical pore pressure perturbation to trigger slip is expressed as
P = P i n j e c t i o n P p
where Pinjection (Pa) is the fluid pressure of the injected fluid.
Figure 10, displayed as a lower hemisphere stereonet, provides a visual representation of the changes in pore pressure required for fractures to initiate slip considering the most likely stress regime. Each point on the stereonet corresponds to a fracture in a particular orientation; the color indicates the magnitude of the necessary pore pressure perturbation for slip initiation. Fractures in the red regions require the smallest ΔP perturbation to cause slip, while those in blue tones require the largest perturbation.
The stereonet suggests that fractures oriented perpendicular to the maximum horizontal stress with a low dip angle, i.e., 30 ± 10°, are under the most critical stress conditions. Specifically, for this scenario, the injection pressure should be approximately ΔP = 10 MPa km−1 higher than the current estimated in situ fluid pressure of 10.9 MPa km−1.
However, the fractures observed in the field have high dip angles (60–80°), making them less critically stressed. Optimally oriented sets that require less injection fluid pressure to be activated are WNW-ESE and N-S (ΔP about 12.5 MPa km−1). The fracture sets striking E-W and NE-SW would require a ΔP of about 17.5 MPa km−1 to induce slip. The least slip-sensitive set is NW-SE, which is oriented parallel to the minimum horizontal stress.
The slip tendency of a fracture is expressed as
τ ( σ n P p ) μ
In a Mohr–Coulomb shear slip assessment, different input parameters (e.g., stress, friction angle) can be treated as random variables with specific statistical parameters. The probabilities of slip can then be described as
P f a i l u r e = P ( τ ( σ n P p ) × μ 0 )
This probabilistic slip tendency approach becomes particularly useful when dealing with large uncertainty in each input parameter. As previously mentioned, and discussed by Miranda et al. [48], there are uncertainties associated with the magnitude and orientation of each principal stress as well as the value of pore pressure. These uncertainties associated with the magnitude of three principal stresses are shown by the green error bars in Figure 9. Uncertainties in the dip direction and dip angle of fractures also exist. A joint/fault is seldom a straight discontinuity plane: it can have different shapes and orientations, ranging from straight or planar fractures to irregular or curved cracks. The latter have even lower slip tendencies because of their macroscopic roughness (deviation from planarity).
To account for these uncertainties in stress, pore pressure and fracture properties, we utilized Monte Carlo probabilistic methods to generate a range of possible outcomes. The Mohr–Coulomb yield criterion was then applied to each of these outcomes to determine the probability of fracture slip. This approach allowed for a more comprehensive analysis of the potential for fracture slip, considering the different uncertainties in the system. Figure 9 and Table 3 present distributions of three principal stress magnitudes and pore pressure as inputs to the analysis. Recall that a constraint was imposed to neglect unrealistic stress scenarios. Our analysis assumes that the friction coefficient ranges between 0.5 and 0.7 and that fracture dip direction and dip angle deviate by up to 5° compared to the field measurement. This analysis was conducted using 10,000 random combinations of parameters for each mapped fracture segment to assess the conditional probability of slip as a function of pore pressure perturbation (ΔP).
Using the uncertainty distributions and fracture properties as randomly selected inputs, the distribution of critical injection pressure calculated from 10,000 scenarios shows a wide range (Figure 11). The curve shows the cumulative probability of exceeding the critical injection pressure, with higher probabilities indicating a greater likelihood of fracture slip. The distribution can be used to estimate the range of injection pressures that are likely to cause fracture slip, as well as the probability of such events occurring under different injection scenarios. Results of the analysis indicate a very low probability of fracture slip that increases with increasing pore pressure. However, when the injection pressure is raised to 26 MPa/km, the probability of slip increases to 68%, which further escalates with a higher injection pressure.
Tornado plots were used to visualize the sensitivity of injection pressure to slip for each input parameter to evaluate the impact of uncertainties in specific parameters on the slip tendency of a particular fracture plane. The parameters considered are ranked based on the extent of variation in the resulting pore pressure to induce slip. The magnitude of the variation indicates which parameters exert the most significant influence on the calculated pore pressure required to induce slip. Figure 12 is an example of a tornado plot, focusing on the fracture that has the highest slip tendency. For this specific case, the fracture is particularly sensitive to variations in the Shmin gradient, as indicated by the error bar associated with Shmin in Figure 9. It is important to note that the sensitivity analysis may vary for each fracture, depending on its unique strike and dip angle.

5. Discussion

The EGS is an emerging technology that could revolutionize energy supply in remote northern communities located on Canadian Shield rocks. However, the EGS is still a challenging technology with many development issues that need to be addressed. Pollack et al. [74] compiled challenges experienced in world-wide EGS projects and found that:
  • 6 sites ceased operations temporarily or shut down due to seismicity or seismicity concerns;
  • 24 EGS projects were delayed or terminated due to drilling and plant operation issues, such as holes and cracks in wellbore casing, stuck drill strings, and well collapses;
  • 18 sites faced challenges in reservoir creation and circulation, such as insufficient connectivity between the injection and production wells or water loss.
Despite all these potential roadblocks, Pollack et al. [74] also found that 29 sites had successful stimulation experiences, and these sites are active and improving today’s knowledge about EGS technology.
The main objective of this technology is to develop a low-impedance hydraulic connection between injector and producer wells that will allow reasonable flow rates, minimal thermal drawdown, minimal water losses, and induced seismicity at levels imperceptible to the population. For these reasons, it is important to forecast how fractures will behave during hydraulic stimulation treatments. This is tightly linked to the natural fracture network properties at the target site, the in situ stress field and the friction properties of the fractures.
Slip tendency analysis is a simple yet effective tool to answer questions that arise at early stages of geothermal exploration:
  • What is the probability of shear slip of pre-existing fractures at the current state of stress?
  • Which orientations of fractures are most likely to be activated?
  • What in situ fluid pressure is required to overcome the shear stress and activate pre-existing fractures?
Considering the results of our analysis, field fracture observations, and empirical stress regime estimates, none of the fractures at the target site are near a critical state of stress, meaning a very low (<1%) probability of slippage under the current state of stress. Nevertheless, if pore pressures change (i.e., hydraulic stimulation treatments or higher circulation pressures), the fracture sets making a 30° angle with the orientation of the maximum horizontal stress (i.e., WNW-ESE and N-S) may be reactivated at an injection pressure (ΔP) of about 12.5 MPa km−1. A ΔP of about 17.5 MPa km−1 may be necessary to reactivate the fracture sets parallel to the maximum horizontal stress (i.e., NE-SW) and the E-W fracture set. The least optimally oriented set is the NW-SE fracture set, which is oriented parallel to the minimum horizontal stress (depicted in blue in Figure 10). Again, it is important to emphasize that these estimates are empirical, and several factors influence how fractures behave during hydraulic stimulation. A more robust analysis including field experiments would be necessary to validate our estimates.
This observation is somewhat in line with what was found with the first hydraulic stimulation experiments at the Carnmenellis Granite, Cornwall, UK [25]. In this EGS site, micro-seismicity patterns indicated anisotropic fracture-controlled flow closely related to the orientation of the maximum horizontal stress [25]. It was deduced that, due to the impact of the anisotropic stress field on fracture apertures, the preferential flow paths were fractures favorably oriented to the maximum horizontal stress. These tend to be activated in shear or in transtension. At Kuujjuaq, our study site, a detailed structural analysis would be the next logical step to verify if a relationship between the fracture aperture and orientation exists.
Additionally, it is important to emphasize that the injection pressure values to activate the fractures are only valid if the most likely scenario in terms of stress prevails in the study area. In instances where the differential stresses are smaller or larger than the average value, the required activation injection pressure is different.
In a high differential stress scenario (σ1σ3 = 29 MPa km−1), the optimally oriented fracture is near their critical stress state, and an injection pressure of ΔP ~ 2 MPa km−1 would likely lead to slip, and induced seismicity would be far more likely. However, in a low differential stress scenario (σ1σ3 = 9 MPa km−1), this fracture would be activated at injection pressures of ΔP ~ 24 MPa km−1 (Figure 11).
The influence of differential stress on hydraulic stimulation effectiveness has been highlighted by Xie et al. [22] who compiled data from seven EGS sites and noticed that where the differential stresses were large, less additional fluid pressure (ΔP) was required to activate shear slip of natural fractures. Their observations showed that the Hijiori EGS project (Japan) is the least stressed and has the smallest differential stress. For this reason, it required the most effort to trigger shear slip compared to, for instance, Rosemanowes (UK), Soultz (France), and Basel (Switzerland) EGS projects. In other words, when the differential stress is small, more fluid pressure is needed to move the Mohr stress circle to meet the slip yield envelope.
It is also important to mention the role played by the intermediate principal stress. The activation of pre-existing fracture planes is dependent on not only the change in the maximum differential stress, but also the value of the intermediate stress. For example, if the intermediate stress was equal to the minimum stress instead of equal to the maximum stress in the worst-case stress scenario, then the required injection pressure to reactivate the fractures would be less than 24 MPa km−1. Similarly, if the intermediate stress was larger than assumed in the best-case stress scenario, then the fractures would not be at a critical state of stress and an injection pressure of about 11 MPa km−1 would be necessary to reactivate them.
The influence of external factors, such as an injection of cold water into a warm medium, can also alter the state of stress, increasing the deviatoric stress. If this increases the shear stress along favorably oriented discontinuities while reducing the normal stress, the fractures move toward a more critical condition. This study only takes into account pressure changes, but further research can be undertaken to consider temperature changes as well.
Another aspect worth discussing is the lithological heterogeneity. The study area is composed mainly by paragneiss rocks. However, diorite, gabbro, tonalite, and granites can be found interlayering the paragneiss. Some rock types may exhibit higher or lower slip activation potential, and this could cause a heterogeneity in the system that would change the analysis undertaken. To avoid adding more uncertainty sources to the analysis, we considered the system to be homogeneous in terms of lithology. However, further research is needed and should be undertaken to deepen the simple analysis carried out in this work. Additionally, anisotropy was also not considered, but it could be included in a more complex model.
The behavior of the Lac Pingiajjulik fault (see Figure 3 for location) during hydraulic stimulation is another question that deserves further investigation. This corresponds to a large-scale shear zone, and one of two scenarios can happen when it is stimulated: either it acts as a preferential flow pathway, or it acts as a barrier to flow circulation. In either case, the EGS must be properly designed to meet the goal of high flow rates, minimal thermal drawdown, minimal water losses, and induced seismicity at levels imperceptible to the population.
Hydraulic shearing experiments were conducted at the Grimsel Test Site in Switzerland under the framework of the In-situ Stimulation and Circulation experiment [15,16]. The rock volume tested is intercepted by two sets of shear zones of ductile (S1.0–S1.3) and brittle-ductile (S3.1–S3.2) origin [75]. The first set of shear zones (S1.0 to S1.3) includes three ductile shear zones, which show strong foliation and mylonitization. The second set (S3.1 and S3.2) includes two shear zones that bind a densely fractured zone in between.
The hydraulic stimulation experiments reveal that fractures associated with one out of two shear zone types were hydraulically reactivated [76]. The authors also noted that the two shear zone types not only differ in terms of tectonic genesis and architecture, but transmissivity change, jacking pressure, and seismic activity were also different for the two shear zone types. These observations led them to suggest that shear zone architectures govern the seismo-hydromechanical response. Based on the in situ experiments, Krietsch et al. [76] observed that elevated fracture fluid pressures associated with the stimulations propagated mostly along the stimulated shear zones, and that flow is channelized within the shear zones. Furthermore, Villiger et al. [77] observed three seismic clusters associated with the stimulation of fractures in the damage zone of the shear zone.
The hydraulic stimulation experiments carried out at the Grimsel Test Site highlight the need to properly characterize the tectonic genesis and architecture of the shear zone crossing Kuujjuaq. Further structural analysis of the area is the next logical step to gain insights into its behavior.
It is also worth discussing how useful hydraulic experiments could be to decrease the uncertainty found in this research work. Hydraulic stimulation experiments and microseismic monitoring are nowadays relatively feasible and could provide important insights into the reservoir development. However, carrying out these field experiments is quite challenging in northern and remote environments where there is no road linking the southern and northern areas and the equipment needs to be shipped by plain or boat. A first-order assessment in this situation is necessary to trigger interest for further investments in geothermal research.

6. Conclusions

Crystalline rocks can host important geothermal resources, but the permeability of these reservoirs is usually too low for the system to be commercially viable. The EGS is an emerging technology that can help increase the productivity of reservoirs by enhancing hydraulic connections within the fracture network. However, the EGS is still a challenging technology: issues associated with reservoir development and circulation and induced seismicity are the main roadblocks to the successful deployment of this technology.
Therefore, understanding how the fracture network will behave when hydraulically stimulated is a key step toward the development of the EGS. The objective of this work was to estimate the slip activation potential of existing fractures in subsurface crystalline basement rocks from outcrop analogues in the Kuujjuaq area and empirical stress predictions.
With this work, we addressed questions that arise at early stages of geothermal exploration, namely “What is the probability of shear slip of the pre-existing fractures at the current state of stress?”, “Which orientations of fractures are most likely to be activated?”, and “What in situ fluid pressure is required to overcome the shear stress and activate pre-existing fractures?”.
Our analysis using a probabilistic slip tendency approach suggests that, at the current state of stress, the probability of the pre-existing fractures slipping is very low (<1%). The analysis also suggests that the optimally oriented fracture sets that can be activated at lower injection pressures are WNW-ESE and N-S. An injection pressure of 12.5 MPa km−1 could be necessary to activate shear slip; however, this is only valid considering the most likely stress scenario. If the differential stress is larger or smaller than the average scenario, then larger or smaller injection pressures are required to reactivate the optimally oriented fractures.
The next logical and needed steps should include more accurate stress measurements than the empirical stress estimates used in this work and a comprehensive structural characterization of the joints, faults, and shear zones in the study area. Hydraulic stimulation experiments would also be needed to not only gain insights into the reservoir and its behavior when pressurized, but also to properly design the EGS to meet the targets of high flow rates, minimal thermal drawdown, minimal water losses, and induced seismicity at levels imperceptible to the population.

Author Contributions

Conceptualization, M.M.M. and A.Y.; methodology, M.M.M. and A.Y.; software, A.Y.; validation, M.M.M., A.Y., J.R., A.W. and M.B.D.; formal analysis, M.M.M. and A.Y.; investigation, M.M.M. and A.Y.; resources, M.M.M. and A.Y.; data curation, M.M.M. and A.Y.; writing—original draft preparation, M.M.M. and A.Y.; writing—review and editing, M.M.M., A.Y., J.R., A.W. and M.B.D.; visualization, M.M.M. and A.Y.; supervision, J.R. and M.B.D.; project administration, J.R.; funding acquisition, J.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Institut Nordique du Québec (INQ) through a research chair awarded to Jasmin Raymond to evaluate the geothermal potential of northern Québec.

Data Availability Statement

All data are presented in the manuscript.

Acknowledgments

The Centre d’études nordiques (CEN), supported by the Fonds de recherche du Québec—nature et technologies (FRQNT), and the Observatoire Homme Milieu Nunavik (OHMI) are acknowledged for helping with field campaigns and logistics. The authors are also thankful to Cynthia Brind’Amour-Coté for her support during the field campaign. The authors would also like to acknowledge the two anonymous reviewers whose comments helped to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Canada’s Renewable Power. Available online: https://www.cer-rec.gc.ca/en/data-analysis/energy-commodities/electricity/report/canadas-renewable-power/index.html (accessed on 14 July 2023).
  2. Gabrys, M. Energy Transition in Canada’s Northern Territories. In Energy Policy Advancement, 1st ed.; Kurochkin, D., Crawford, M.J., Shabliy, E.V., Eds.; Springer: Cham, Switzerland, 2022; pp. 123–143. [Google Scholar]
  3. Pinto, H.; Gates, I.D. Why is it so difficult to replace diesel in Nunavut, Canada? Renew. Sustain. Energy Rev. 2022, 157, 112030. [Google Scholar] [CrossRef]
  4. Geology of Nunavut. Available online: https://geoscan.nrcan.gc.ca/starweb/geoscan/servlet.starweb?path=geoscan/fulle.web&search1=R=223684 (accessed on 14 July 2023).
  5. Services Web Cartographiques. Available online: https://sigeom.mines.gouv.qc.ca/signet/classes/I1108_afchCarteIntr (accessed on 14 July 2023).
  6. Grasby, S.E.; Allen, D.M.; Bell, S.; Chen, Z.; Ferguson, G.; Jessop, A.; Kelman, M.; Ko, M.; Majorowicz, J.; Moore, M.; et al. Geothermal Energy Resource Potential of Canada; Report No.: Open File 6914; Geological Survey of Canada: Ottawa, ON, Canada, 2012.
  7. Arriaga, M.; Brooks, M.; Moore, N. Energy Access—The Canadian Context; Waterloo Global Initiative’s OpenAccess Energy Blueprint; WGSI: Waterloo, ON, Canada, 2017. [Google Scholar]
  8. MacDonald, P.; Stedman, A.; Symons, G. The UK Geothermal Hot Dry Rock R&D Programme. In Proceedings of the Seventeenth Workshop on Geothermal Reservoir Engineering, Stanford, CA, USA, 29–31 January 1982. [Google Scholar]
  9. Richards, H.G.; Parker, R.H.; Green, A.S.P.; Jones, R.H.; Nicholls, J.D.M.; Nicol, D.A.C.; Randall, M.M.; Richards, S.; Stewart, R.C.; Willis-Richards, J. The performance and characteristics of the experimental hot dry rock geothermal reservoir at Rosemanowes, Cornwall (1985–1988). Geothermics 1994, 23, 73–109. [Google Scholar] [CrossRef]
  10. Brown, D.W.; Duchane, D.V. Scientific progress on the Fenton Hill HDR project since 1983. Geothermics 1999, 28, 591–601. [Google Scholar] [CrossRef]
  11. Parker, R. The Rosemanowes HDR project 1983–1991. Geothermics 1999, 28, 603–615. [Google Scholar] [CrossRef]
  12. Genter, A.; Evans, K.; Cuenot, N.; Fritsch, D.; Sanjuan, B. Contribution of the exploration of deep crystalline fractured reservoir of Soultz to the knowledge of enhanced geothermal systems (EGS). Comptes Rendus Geosci. 2010, 342, 502–516. [Google Scholar] [CrossRef]
  13. Brown, D.W.; Duchane, D.V.; Heiken, G.; Hriscu, V.T. Mining the Earth’s Heat: Hot Dry Rock Geothermal Energy, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  14. Kelkar, S.; WoldeGabriel, G.; Rehfeldt, K. Lessons learned from the pioneering hot dry rock project at Fenton Hill, USA. Geothermics 2016, 63, 5–14. [Google Scholar] [CrossRef]
  15. Amann, F.; Gischig, V.; Evans, K.; Doetsch, J.; Jalali, R.; Valley, B.; Krietsch, H.; Dutler, N.; Villiger, L.; Brixel, B.; et al. The seismo-hydromechanical behavior during deep geothermal reservoir stimulations: Open questions tackled in a decameter-scale in situ stimulation experiment. Solid Earth 2018, 9, 115–137. [Google Scholar] [CrossRef]
  16. Gischig, V.S.; Giardini, D.; Amann, F.; Hertrich, M.; Krietsch, H.; Doetsch, S.; Doonechaly, N.G.; Driesner, T.; Dutler, N.; Evans, K.F.; et al. Hydraulic stimulation and fluid circulation experiments in underground laboratories: Stepping up the scale towards engineered geothermal systems. Geomech. Energy Environ. 2020, 24, 100175. [Google Scholar] [CrossRef]
  17. McClure, M.W.; Horne, R.N. An investigation of stimulation mechanisms in Enhanced Geothermal Systems. Int. J. Rock Mech. Min. Sci. 2014, 72, 242–260. [Google Scholar] [CrossRef]
  18. Norbeck, J.H.; McClure, M.W.; Horne, R.N. Field observations at the Fenton Hill enhanced geothermal system test site support mixed-mechanism stimulation. Geothermics 2018, 74, 135–149. [Google Scholar] [CrossRef]
  19. Li, S.; Wang, S.; Tang, H. Stimulation mechanism and design of enhanced geothermal systems: A comprehensive review. Renew. Sustain. Energy Rev. 2022, 155, 111914. [Google Scholar] [CrossRef]
  20. Peacock, D.C.P.; Sanderson, D.J.; Leiss, B. Use of Mohr Diagrams to Predict Fracturing in a Potential Geothermal Reservoir. Geosciences 2021, 11, 501. [Google Scholar] [CrossRef]
  21. Sibson, R.H. Reverse fault rupturing: Competition between non-optimal and optimal fault orientations. Geol. Soc. Lond. Spec. Publ. 2012, 367, 39–50. [Google Scholar] [CrossRef]
  22. Xie, L.; Min, K.-B.; Song, Y. Observations of hydraulic stimulations in seven enhanced geothermal system projects. Renew. Energy 2015, 79, 56–65. [Google Scholar] [CrossRef]
  23. Majer, E.L.; Baria, R.; Stark, M.; Oates, S.; Bommer, J.; Smith, B.; Asanuma, H. Induced seismicity associated with Enhanced Geothermal Systems. Geothermics 2007, 36, 185–222. [Google Scholar] [CrossRef]
  24. Buijze, L.; van Bijsterveldt, L.; Cremer, H.; Paap, B.; Veldkamp, H.; Wassing, B.B.T.; van Wees, J.-D.; van Yperen, G.C.N.; ter Heege, J.H.; Jaarsma, B. Review of induced seismicity in geothermal systems worldwide and implications for geothermal systems in the Netherlands. Neth. J. Geosci. 2019, 98, e13. [Google Scholar] [CrossRef]
  25. Pine, R.J.; Batchelor, A.S. Downward migration of shearing in jointed rock during hydraulic injections. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1984, 21, 249–263. [Google Scholar] [CrossRef]
  26. Deichmann, N.; Giardini, D. Earthquakes Induced by the Stimulation of an Enhanced Geothermal System below Basel (Switzerland). Seismol. Res. Lett. 2009, 80, 784–798. [Google Scholar] [CrossRef]
  27. Zang, A.; Oye, V.; Jousset, P.; Deichmann, N.; Gritto, R.; McGarr, A.; Majer, E.; Bruhn, D. Analysis of induced seismicity in geothermal reservoirs—An overview. Geothermics 2014, 52, 6–21. [Google Scholar] [CrossRef]
  28. Moeck, I.; Kwiatek, G.; Zimmermann, G. Slip tendency analysis, fault reactivation potential and induced seismicity in a deep geothermal reservoir. J. Struct. Geol. 2009, 31, 1174–1182. [Google Scholar] [CrossRef]
  29. Morris, A.; Ferrill, D.A.; Henderson, D.B. Slip-tendency analysis and fault reactivation. Geology 1996, 24, 275–278. [Google Scholar] [CrossRef]
  30. Ito, T.; Hayashi, K. Role of Stress-controlled Flow Pathways in HDR Geothermal Reservoirs. Pure Appl. Geophys. 2003, 160, 1103–1124. [Google Scholar] [CrossRef]
  31. Yaghoubi, A.; Hickson, C.J.; Leonenko, Y.; Dusseault, M.B. Probabilistic assessment of induced seismicity at the Alberta No. 1 geothermal project site. Can. J. Earth Sci. 2023, 60, 294–306. [Google Scholar] [CrossRef]
  32. Bauer, J.F.; Krumbholz, M.; Meier, S.; Tanner, D.C. Predictability of properties of a fractured geothermal reservoir: The opportunities and limitations of an outcrop analogue study. Geotherm. Energy 2017, 5, 24. [Google Scholar] [CrossRef]
  33. Wydte, L.M.; Bär, K.; Colombero, C.; Comina, C.; Deb, P.; Lepillier, B.; Mandrone, G.; Milsch, H.; Rochelle, C.A.; Vagnon, F.; et al. Outcrop analogue study to determine reservoir properties of the Los Humeros and Acoculco geothermal fields, Mexico. Adv. Geosci. 2018, 45, 281–287. [Google Scholar] [CrossRef]
  34. Lepillier, B.; Bruna, P.-O.; Bruhn, D.; Bastesen, E.; Daniilidis, A.; Garcia, O.; Torabi, A.; Wheeler, W. From outcrop scanlines to discrete fracture networks, an integrative workflow. J. Struct. Geol. 2020, 133, 103992. [Google Scholar] [CrossRef]
  35. Peacock, D.C.P.; Sanderson, D.J.; Leiss, B. Use of Analogue Exposures of Fractured Rock for Enhanced Geothermal Systems. Geosciences 2022, 12, 318. [Google Scholar] [CrossRef]
  36. Hudson, J.A.; Cornet, F.H.; Christiansson, R. ISRM Suggested Methods for rock stress estimation—Part 1: Strategy for rock stress estimation. Int. J. Rock Mech. Min. Sci. 2003, 40, 991–998. [Google Scholar] [CrossRef]
  37. Sjöberg, J.; Christiansson, R.; Hudson, J.A. ISRM Suggested Methods for rock stress estimation—Part 2: Overcoring methods. Int. J. Rock Mech. Min. Sci. 2003, 40, 999–1010. [Google Scholar] [CrossRef]
  38. Haimson, B.C.; Cornet, F.H. ISRM Suggested Methods for rock stress estimation—Part 3: Hydraulic fracturing (HF) and/or hydraulic testing of pre-existing fractures (HTPF). Int. J. Rock Mech. Min. Sci. 2003, 40, 1011–1020. [Google Scholar] [CrossRef]
  39. Stephansson, O.; Zang, A. ISRM Suggested Methods for Rock Stress Estimation—Part 5: Establishing a Model for the In Situ Stress at a Given Site. In The ISRM Suggested Methods for Rock Characterization, Testing and Monitoring: 2007–2014, 1st ed.; Ulusay, R., Ed.; Springer: Cham, Switzerland, 2014; pp. 187–201. [Google Scholar]
  40. Yaghoubi, A.; Dusseault, M.B.; Leonenko, Y. Injection-induced fault slip assessment in Montney Formation in Western Canada. Sci. Rep. 2022, 12, 11551. [Google Scholar] [CrossRef]
  41. Giordano, N.; Raymond, J. Alternative and sustainable heat production for drinking water needs in a subarctic climate (Nunavik, Canada): Borehole thermal energy storage to reduce fossil fuel dependency in off-grid communities. Appl. Energy 2019, 252, 113463. [Google Scholar] [CrossRef]
  42. Kanzari, I. Évaluation du Potentiel des Pompes à Chaleur Géothermique Pour La Communauté Nordique de Kuujjuaq. Master’s Thesis, Université du Québec, Institut national de la Recherche Scientifique, Québec, QC, Canada, 2019. [Google Scholar]
  43. Gunawan, E.; Giordano, N.; Jensson, P.; Newson, J.; Raymond, J. Alternative heating systems for northern remote communities: Techno-economic analysis of ground-coupled heat pumps in Kuujjuaq, Nunavik, Canada. Renew. Energy 2020, 147, 1540–1553. [Google Scholar] [CrossRef]
  44. Miranda, M.M.; Giordano, N.; Raymond, J.; Pereira, A.J.S.C.; Dezayes, C. Thermophysical properties of surficial rocks: A tool to characterize geothermal resources of remote northern regions. Geotherm. Energy 2020, 8, 4. [Google Scholar] [CrossRef]
  45. Miranda, M.M.; Raymond, J.; Dezayes, C. Uncertainty and Risk Evaluation of Deep Geothermal Energy Source for Heat Production and Electricity Generation in Remote Northern Regions. Energies 2020, 13, 4221. [Google Scholar] [CrossRef]
  46. Miranda, M.M.; Velez Marquez, M.I.; Raymond, J.; Dezayes, C. A numerical approach to infer terrestrial heat flux from shallow temperature profiles in remote northern regions. Geothermics 2021, 93, 102064. [Google Scholar] [CrossRef]
  47. Miranda, M.M.; Raymond, J.; Willis-Richards, J.; Dezayes, C. Are Engineered Geothermal Energy Systems a Viable Solution for Arctic Off-Grid Communities? A Techno-Economic Study. Water 2021, 13, 3526. [Google Scholar] [CrossRef]
  48. Miranda, M.M.; Raymond, J.; Dezayes, C. Estimating theoretical stress regime for engineered geothermal energy systems in an arctic community (Kuujjuaq, Canada). Comptes Rendus Geosci. 2023, 355, 85–108. [Google Scholar] [CrossRef]
  49. Miranda, M.M.; Raymond, J.; Dezayes, C.; Wigston, A.; Perreault, S. Multiscale fracture networks and their impact on hydroshearing response in the Canadian Shield (Kuujjuaq, Canada). Geomech. Geophys. Geo-Energy Geo-Resour. 2023, 9, 85. [Google Scholar] [CrossRef]
  50. Miranda, M.M.; Raymond, J. Assessing Kuujjuaq’s (Nunavik, Canada) Deep Geothermal Energy Potential: Core Analysis, Thermal Properties Characterization and Surface Heat Flux Estimation of a 234 m Deep Geothermal Exploration Well; Research Report R2109; INRS: Québec, QC, Canada, 2023. [Google Scholar]
  51. What is Magnetic Declination in Kuujjuaq, Canada? Available online: https://www.magnetic-declination.com/Canada/Kuujjuaq/335810.html (accessed on 16 July 2023).
  52. Zeeb, C.; Gomez-Rivas, E.; Bons, P.D.; Blum, P. Evaluation of sampling methods for fracture network characterization using outcrops. AAPG Bull. 2013, 97, 1545–1566. [Google Scholar] [CrossRef]
  53. Watkins, H.; Bond, C.E.; Healy, D.; Butler, R.W.H. Appraisal of fracture sampling methods and a new workflow to characterise heterogeneous fracture networks at outcrop. J. Struct. Geol. 2015, 72, 67–82. [Google Scholar] [CrossRef]
  54. Grapher. Available online: https://www.goldensoftware.com/products/grapher (accessed on 16 July 2023).
  55. Stereonet. Available online: https://www.rickallmendinger.net/stereonet (accessed on 16 July 2023).
  56. World Stress Map Project. Available online: https://www.world-stress-map.org/ (accessed on 16 July 2023).
  57. Simplified Seismic Hazard Map for Canada, the Provinces and Territories. Available online: https://www.seismescanada.rncan.gc.ca/hazard-alea/simphaz-en.php#QC (accessed on 16 July 2023).
  58. Adams, J. Crustal stresses in Eastern Canada. In Earthquakes at North-Atlantic Passive Margins: Neotectonics and Postglacial Rebound; Gregersen, S., Basham, P.W., Eds.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1989; pp. 289–298. [Google Scholar]
  59. Arjang, B. Pre-mining stresses at some hard rock mines in the Canadian Shield. CIM Bull. 1991, 84, 80–86. [Google Scholar]
  60. Arjang, B. Canadian crustal stresses and their application in mine design. In Mine Planning and Equipment Selection 1998; Singhal, R.K., Ed.; A. A. Balkema: Rotterdam, The Netherlands, 1998; pp. 269–274. [Google Scholar]
  61. Herget, G. High stress occurrences in the Canadian Shield. In Proceedings of the 23rd US Symposium on Rock Mechanics (USRMS), Berkeley, CA, USA, 25 August 1982. [Google Scholar]
  62. Herget, G. Stress assumptions for underground excavations in the Canadian shield. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1987, 24, 95–97. [Google Scholar] [CrossRef]
  63. Herget, G. Rock stresses and rock stress monitoring in Canada. In Rock Testing and Site Characterization: Principles, Practice and Projects; Hudson, J.A., Ed.; Elsevier: Berlin, Germany, 1993; pp. 473–496. [Google Scholar]
  64. Young, S.; Maloney, S. An Update to the Canadian Shield Stress Database; NWMO-TR-2015-18; NWMO: Toronto, ON, Canada, 2015. [Google Scholar]
  65. Andrieux, P.P.; Li, H.; Brummer, R.K.; Caumartin, R.; Rispoli, A. Rock mechanics approach for the recovery of the zone G crown pillar at the Raglan Katinniq Mine. In Proceedings of the 3rd CANUS Rock Mechanics Symposium, Toronto, ON, Canada, 9–15 May 2009. [Google Scholar]
  66. Kaiser, P.K.; Maloney, S.M.; Yong, S. Role of large scale heterogeneities on in-situ stress and induced stress fields. In Proceedings of the 50th US Rock Mechanics/Geomechanics Symposium, Houston, TX, USA, 26–29 June 2016. [Google Scholar]
  67. Hoek, E.; Brown, E.T. Underground Excavation in Rock; Institute of Mining & Metallurgy: London, UK, 1980. [Google Scholar]
  68. Zhang, L. Engineering Properties of Rocks; Elsevier: Oxford, UK, 2017. [Google Scholar]
  69. Sibson, R.H. Controls on maximum fluid overpressure defining conditions for mesozonal mineralization. J. Struct. Geol. 2004, 26, 1127–1136. [Google Scholar] [CrossRef]
  70. Zoback, M.D.; Townend, J. Implications of hydrostatic pore pressures and high crustal strength for the deformation of intraplate lithosphere. Tectonophysics 2001, 336, 19–30. [Google Scholar] [CrossRef]
  71. Jaeger, J.C.; Cook, N.G.W.; Zimmerman, R.W. Fundamentals of Rock Mechanics; Blackwell Publishing: Oxford, UK, 2007. [Google Scholar]
  72. Byerlee, J. Friction of rocks. Pure Appl. Geophys. 1978, 116, 615–626. [Google Scholar] [CrossRef]
  73. Mahbaz, S.B.; Yaghoubi, A.; Dehghani-Sanij, A.; Sarvaramini, E.; Leonenko, Y.; Dusseault, M.B. Well-Doublets: A First-Order Assessment of Geothermal SedHeat Systems. Appl. Sci. 2021, 11, 697. [Google Scholar] [CrossRef]
  74. Pollack, A.; Horne, R.; Mukerji, T. What Are the Challenges in Developing Enhanced Geothermal Systems (EGS)? Observations from 64 EGS Sites. In Proceedings of the World Geothermal Congress 2020+1, Reykjavik, Iceland, 24–27 October 2021. [Google Scholar]
  75. Krietsch, H.; Doetsch, J.; Dutler, N.; Jalali, M.; Gischig, V.; Loew, S.; Amann, F. Comprehensive geological dataset describing a crystalline rock mass for hydraulic stimulation experiments. Sci. Data 2018, 5, 180269. [Google Scholar] [CrossRef]
  76. Krietsch, H.; Gischig, V.S.; Doetsch, J.; Evans, K.F.; Villiger, L.; Jalali, M.; Valley, B.; Löw, S.; Amann, F. Hydromechanical processes and their influence on the stimulation effected volume: Observations from a decameter-scale hydraulic stimulation project. Solid Earth 2020, 11, 1699–1729. [Google Scholar] [CrossRef]
  77. Villiger, L.; Gischig, V.S.; Kwiatek, G.; Krietsch, H.; Doetsch, J.; Jalali, M.; Amann, F.; Giardini, D.; Wiemer, S. Metre-scale stress heterogeneities and stress redistribution drive complex fracture slip and fracture growth during a hydraulic stimulation experiment. Geophys. J. Int. 2021, 225, 1689–1703. [Google Scholar] [CrossRef]
Figure 2. Example of Mohr diagram showing the fields in which extension (1), hybrid (2) and shear (3) fractures occur [20]. Reprinted/adapted with permission from Ref. [20]. Copyright 2021, The Authors.
Figure 2. Example of Mohr diagram showing the fields in which extension (1), hybrid (2) and shear (3) fractures occur [20]. Reprinted/adapted with permission from Ref. [20]. Copyright 2021, The Authors.
Geosciences 13 00340 g002
Figure 3. Example of scanline sampling method on an exposed outcrop surface in the study area.
Figure 3. Example of scanline sampling method on an exposed outcrop surface in the study area.
Geosciences 13 00340 g003
Figure 4. Geological map of Kuujjuaq (adapted from [5]) with indication of each fracture sampling area (A1–A5) and rose diagrams produced from data collection. LP—Lac Pingiajjulik fault. The rose diagrams were produced using the software GrapherTM version 19.3.323 26 March 2022, from Golden Software [54].
Figure 4. Geological map of Kuujjuaq (adapted from [5]) with indication of each fracture sampling area (A1–A5) and rose diagrams produced from data collection. LP—Lac Pingiajjulik fault. The rose diagrams were produced using the software GrapherTM version 19.3.323 26 March 2022, from Golden Software [54].
Geosciences 13 00340 g004
Figure 5. (a) Stereographic projection of the poles to planes and 1% area contour, (b) rose diagram and (c) relative frequency of the fracture strike planes. F—Lac Pingiajjulik fault plane. Stereographic projection of the poles to planes and 1% area contour were produced using the software Stereonet version 11.3.7 of Allmendinger et al. [55]. Rose diagram was produced using the software GrapherTM version 19.3.323 26 March 2022, from Golden Software [54].
Figure 5. (a) Stereographic projection of the poles to planes and 1% area contour, (b) rose diagram and (c) relative frequency of the fracture strike planes. F—Lac Pingiajjulik fault plane. Stereographic projection of the poles to planes and 1% area contour were produced using the software Stereonet version 11.3.7 of Allmendinger et al. [55]. Rose diagram was produced using the software GrapherTM version 19.3.323 26 March 2022, from Golden Software [54].
Geosciences 13 00340 g005
Figure 6. Outcrop on sampling area A6 fractures identified and stereonet canvas produced from fracture interpretation [49]. Reprinted/adapted with permission from Ref. [49]. Copyright 2023, The Authors. Red line in the stereonet indicates the outcrop orientation. Stereonet canvas produced using the software Stereonet version 11.3.7 of Allmendinger et al. [55].
Figure 6. Outcrop on sampling area A6 fractures identified and stereonet canvas produced from fracture interpretation [49]. Reprinted/adapted with permission from Ref. [49]. Copyright 2023, The Authors. Red line in the stereonet indicates the outcrop orientation. Stereonet canvas produced using the software Stereonet version 11.3.7 of Allmendinger et al. [55].
Geosciences 13 00340 g006
Figure 7. (a) Stress map [56] and (b) seismic hazard map [57] of Québec province [48]. Reprinted/adapted with permission from Ref. [48]. Copyright 2023, The Authors. The arrows in (a) indicate the regional trend of the contemporary stress field and the orientation of the maximum horizontal compression based on Adams [58].
Figure 7. (a) Stress map [56] and (b) seismic hazard map [57] of Québec province [48]. Reprinted/adapted with permission from Ref. [48]. Copyright 2023, The Authors. The arrows in (a) indicate the regional trend of the contemporary stress field and the orientation of the maximum horizontal compression based on Adams [58].
Geosciences 13 00340 g007
Figure 8. Distribution of each principal stress and in situ fluid pressure. SV—vertical stress; Shmin—minimum horizontal stress; SHmax—maximum horizontal stress; Pp—in situ fluid pressure.
Figure 8. Distribution of each principal stress and in situ fluid pressure. SV—vertical stress; Shmin—minimum horizontal stress; SHmax—maximum horizontal stress; Pp—in situ fluid pressure.
Geosciences 13 00340 g008
Figure 9. A 3D Mohr diagram for the average stress magnitude gradient and uncertainty associated with each principal stress component.
Figure 9. A 3D Mohr diagram for the average stress magnitude gradient and uncertainty associated with each principal stress component.
Geosciences 13 00340 g009
Figure 10. Lower hemisphere stereonet showing critical injection pressure perturbation leading to fracture slip. Each point in the graph represents a fracture pole identified in the field.
Figure 10. Lower hemisphere stereonet showing critical injection pressure perturbation leading to fracture slip. Each point in the graph represents a fracture pole identified in the field.
Geosciences 13 00340 g010
Figure 11. Each curve represents the cumulative probability function of slip on each fracture at current pressure. The hot color curve represents fractures with a higher slip tendency. The higher the injection pressure, the more likely the fractures are to slip. An injection pressure of 30 MPa/km is expected to result in a 95% fracture slip.
Figure 11. Each curve represents the cumulative probability function of slip on each fracture at current pressure. The hot color curve represents fractures with a higher slip tendency. The higher the injection pressure, the more likely the fractures are to slip. An injection pressure of 30 MPa/km is expected to result in a 95% fracture slip.
Geosciences 13 00340 g011
Figure 12. Tornado plot displaying a sensitivity analysis of uncertain parameters affecting the slip tendency of the optimally oriented fracture.
Figure 12. Tornado plot displaying a sensitivity analysis of uncertain parameters affecting the slip tendency of the optimally oriented fracture.
Geosciences 13 00340 g012
Table 1. Principal stresses in the Canadian Shield based on the literature review [59,60,61,62,63,64,65,66].
Table 1. Principal stresses in the Canadian Shield based on the literature review [59,60,61,62,63,64,65,66].
Principal StressOrientationMagnitudeObservationsRef.
SV---(0.0260–0.0324)z0 < z < 2200 mSV < SH,average[61,62]
SH,average---9.9 + 0.0371z0 < z < 900 m
33.4 + 0.0111z900 < z 2200 m
SV---(0.0266 ± 0.008)z60 < z < 1890 mSV < Shmin < SHmax[59]
SH,average---5.9 + 0.0349z
SHmax---8.2 + 0.0422z
Shmin---3.6 + 0.0276z
SV---0.0285z0 < z < 2200 mSV < Shmin < SHmax[63]
S1N248°/10°12.1 + (0.0403 ± 0.0020)z
S2N300–340°/0°6.4 + (0.0293 ± 0.0019)z
S3Vertical1.4 + (0.0225 ± 0.0015)z
SV---0.0260z0 < z < 6000 mSV < Shmin < SHmax[60]
S1NE/horizontal13.5 + 0.0344z
S2NW/sub-horizontal8.0 + 0.0233z
S3Vertical3.0 + 0.0180z
S1N-S/horizontal0.0513z---SV < Shmin < SHmax[65]
S2E-W/horizontal0.0378z
S3Vertical0.0270z
SV---(0.0258–0.0263)z12 < z < 2552 mSV < Shmin < SHmax[64]
S1N227°/02°(0.040 ± 0.001)z − (9.2 ± 1.5)
S2N310°/08°(0.029 ± 0.001)z + (4.6 ± 1.159)
S3N270°/88°(0.021 ± 0.001)z − (0.8 ± 0.872)
SV---0.021z0 < z < 1300 mSV < Shmin < SHmax[66]
S1---0.012z + 42.4660 < z < 1300 m
S2---0.013z + 24.1
S3---0.007z + 9.7
SV—vertical stress; SHmax—maximum horizontal stress; Shmin—minimum horizontal stress; SH,average—average horizontal stress ( S H m a x + S h m i n 2 ); S1—maximum principal stress; S2—intermediate principal stress; S3—minimum principal stress.
Table 2. Horizontal to vertical stress ratios inferred for the Canadian Shield [60,61,62].
Table 2. Horizontal to vertical stress ratios inferred for the Canadian Shield [60,61,62].
Stress Ratio CoefficientExpressionReference
kmax 357 z + 1.46 [62]
kmin 167 z + 1.10 [62]
kmax 272 ± 8 z + 1.72 [63]
kmin 30 ± 4 z + 0.86 [63]
kmax 7.44 z 0.198 [60]
kmin 2.81 z 0.120 [60]
kmax—maximum horizontal to vertical stress ratio; kmin—minimum horizontal to vertical stress ratio; z—depth in meters.
Table 3. Estimated stress components [47].
Table 3. Estimated stress components [47].
Depth
(km)
Statistical ParametersPrincipal Stresses
SV
(MPa)
Shmin
(MPa)
SHmax
(MPa)
Pp
(MPa)
1Mean27304211
St dev1.33.05.70.5
Median26.830.342.910.9
[Min–Max][24–30][22–38][28–56][10–12]
2Mean54577522
St dev2.65.010.20.9
Median53.757.075.721.7
[Min–Max][48–61][44–71][49–98][19–24]
3Mean818310733
St dev3.96.814.41.3
Median80.583.6108.132.6
[Min–Max][73–91][65–103][71–139][29–36]
4Mean10811013844
St dev5.28.718.71.8
Median107.4109.9140.043.5
[Min–Max][97–122][85–136][93–185][39–48]
5Mean13513617054
St dev6.610.722.72.2
Median134.2137.1171.254.3
[Min–Max][121–152][107–168][113–226][49–61]
SV—vertical stress; SHmax—maximum horizontal stress; Shmin—minimum horizontal stress; Pp—in situ fluid pressure.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Miranda, M.M.; Yaghoubi, A.; Raymond, J.; Wigston, A.; Dusseault, M.B. Slip Activation Potential of Fractures in the Crystalline Basement Rocks of Kuujjuaq (Nunavik, Canada) to Assess Enhanced Geothermal Systems Development. Geosciences 2023, 13, 340. https://doi.org/10.3390/geosciences13110340

AMA Style

Miranda MM, Yaghoubi A, Raymond J, Wigston A, Dusseault MB. Slip Activation Potential of Fractures in the Crystalline Basement Rocks of Kuujjuaq (Nunavik, Canada) to Assess Enhanced Geothermal Systems Development. Geosciences. 2023; 13(11):340. https://doi.org/10.3390/geosciences13110340

Chicago/Turabian Style

Miranda, Mafalda M., Ali Yaghoubi, Jasmin Raymond, Andrew Wigston, and Maurice B. Dusseault. 2023. "Slip Activation Potential of Fractures in the Crystalline Basement Rocks of Kuujjuaq (Nunavik, Canada) to Assess Enhanced Geothermal Systems Development" Geosciences 13, no. 11: 340. https://doi.org/10.3390/geosciences13110340

APA Style

Miranda, M. M., Yaghoubi, A., Raymond, J., Wigston, A., & Dusseault, M. B. (2023). Slip Activation Potential of Fractures in the Crystalline Basement Rocks of Kuujjuaq (Nunavik, Canada) to Assess Enhanced Geothermal Systems Development. Geosciences, 13(11), 340. https://doi.org/10.3390/geosciences13110340

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop