Empirical Analysis and Geomechanical Modelling of an Underground Water Reservoir for Hydroelectric Power Plants

: The European Union policy of encouraging renewable energy sources and a sustainable and safe low-carbon economy requires ﬂexible energy storage systems (FESSs), such as pumped-storage hydropower (PSH) systems. Energy storage systems are the key to facilitate a high penetration of the renewable energy sources in the electrical grids. Disused mining structures in closed underground coal mines in NW Spain have been selected as a case study to analyze the construction of underground pumped-storage hydropower (UPSH) plants. Mine water, depth and subsurface space in closured coal mines may be used for the construction of FESSs with reduced environmental impacts. This paper analyzes the stability of a network of tunnels used as a lower water reservoir at 450 m depth in sandstone and shale formations. Empirical methods based on rock mass classiﬁcation systems are employed to preliminarily design the support systems and to determinate the rock mass properties. In addition, 3D numerical modelling has been conducted in order to verify the stability of the underground excavations. The deformations and thickness of the excavation damage zones (EDZs) around the excavations have been evaluated in the simulations without considering a support system and considering systematic grouted rock bolts and a layer of reinforced shotcrete as support system. The results obtained show that the excavation of the network of tunnels is technically feasible with the support system that has been designed.


Introduction
In the current evolving energy context, characterized by an increasing of variable renewable energy (VRE) in the electricity mix, the development of flexible energy storage systems (FESSs), such as pumped storage hydropower (PSH) or compressed air energy storage (CAES) plants are required. PSH is the most mature technology to provide ancillary services to the electrical grid [1]. PSH systems accounted for 150 GW worldwide in 2016 (40 GW in the European Union) [2] and the capacity could be 325 GW in 2030 [3].
Underground pumped storage hydropower (UPSH) is an alternative to store large amounts of electrical energy with low environmental impacts [4,5]. Other energy storage systems such as and the power input of the turbines depends on the time at full load. The energy storage reaches 480 MWh cycle −1 considering a gross head of 450 mH2O and a useful water capacity of 450,000 m 3 .

Geology
The Asturian Central Coal Basin (ACCB) is a coal mining area located in NW Spain. It is precisely in the ACCB where the terrigenous carboniferous sediments are best represented, reaching a thickness close to 6000 m between the Namurian B and Westfalian D. The succession is divided at large scale into two sectors: one lower denominated Lena Group, characterized by limestones and thin seams of coal, and another superior called Sama Group, in which limestones are less abundant and there are levels of sandstones, shales and seams of exploitable bituminous and hard coal [8].

Underground Hydroelectric Power Plant
UPSH systems consist of two water reservoirs; the upper reservoir is located above ground, while the lower reservoir is underground. The scheme of an underground hydroelectric power plant and lower water reservoir is shown in Figure 2a. In the present work, a subsurface water reservoir conformed by a new network of tunnels with an arched roof and straight walls cross-section of 30 m 2 ( Figure 2b) is considered. The water reservoir consists of a central tunnel connected to the ventilation shaft, and 200 m length of transversal tunnels with a distance between them of 20 m. To avoid water leakage, the tunnels' surface is covered with an impermeable high-strength membrane.

Geology
The Asturian Central Coal Basin (ACCB) is a coal mining area located in NW Spain. It is precisely in the ACCB where the terrigenous carboniferous sediments are best represented, reaching a thickness close to 6000 m between the Namurian B and Westfalian D. The succession is divided at large scale into two sectors: one lower denominated Lena Group, characterized by limestones and thin seams of coal, and another superior called Sama Group, in which limestones are less abundant and there are levels of sandstones, shales and seams of exploitable bituminous and hard coal [8].

Underground Hydroelectric Power Plant
UPSH systems consist of two water reservoirs; the upper reservoir is located above ground, while the lower reservoir is underground. The scheme of an underground hydroelectric power plant and lower water reservoir is shown in Figure 2a. In the present work, a subsurface water reservoir conformed by a new network of tunnels with an arched roof and straight walls cross-section of 30 m 2 ( Figure 2b) is considered. The water reservoir consists of a central tunnel connected to the ventilation shaft, and 200 m length of transversal tunnels with a distance between them of 20 m. To avoid water leakage, the tunnels' surface is covered with an impermeable high-strength membrane.

Geology
The Asturian Central Coal Basin (ACCB) is a coal mining area located in NW Spain. It is precisely in the ACCB where the terrigenous carboniferous sediments are best represented, reaching a thickness close to 6000 m between the Namurian B and Westfalian D. The succession is divided at large scale into two sectors: one lower denominated Lena Group, characterized by limestones and thin seams of coal, and another superior called Sama Group, in which limestones are less abundant and there are levels of sandstones, shales and seams of exploitable bituminous and hard coal [8].

Underground Hydroelectric Power Plant
UPSH systems consist of two water reservoirs; the upper reservoir is located above ground, while the lower reservoir is underground. The scheme of an underground hydroelectric power plant and lower water reservoir is shown in Figure 2a. In the present work, a subsurface water reservoir conformed by a new network of tunnels with an arched roof and straight walls cross-section of 30 m 2 ( Figure 2b) is considered. The water reservoir consists of a central tunnel connected to the ventilation shaft, and 200 m length of transversal tunnels with a distance between them of 20 m. To avoid water leakage, the tunnels' surface is covered with an impermeable high-strength membrane.

Empirical Analysis
The Q-system method was developed by Barton et al. [30] in the Norwegian Geotechnical Institute (NGI) to classify rock masses. The method is based on about 200 case histories of caverns and tunnels. The Q-system has been updated and it is now based on 1260 case records [31]. Barton carried out some changes and adapted it to give support recommendations due to the increasing use of steel fibre reinforced shotcrete (SFR) in underground excavation support [32]. The quality of rock masses (Q) is obtained by applying Equation (1).
where RQD is the rock quality designation, J n is the joint set number, J r is the joint roughness, J a is the joint alteration, J w is the water reduction factor and SRF is the stress reduction factor. The Q value varies from 0.001 for exceptionally poor quality to 1000 for an exceptionally good quality rock mass. A stress-free form Q N was defined by Goel et al. [33], which is given by Equation (2).
Barton defined a new parameter Q C (Equation (3)) to improve the correlation among the engineering parameters, where σ ci is the strength of intact rock in MPa [32].
The bolt length (L b ) may be calculated in terms of the excavation width B, by the Equation (4), proposed by Barton et al. [30].

Material Properties
The physical and mechanical properties of the shale and sandstone were obtained from laboratory testing at the University of Oviedo Rock Mechanics Laboratory on intact rock samples following International Society Rock Mechanics (ISRM) methods [34]. Table 1 shows the properties of the intact rocks and the geological strength index (GSI). The rock mass properties were determined by means of empirical methods considering a blast damage factor D = 0.8. Table 2 shows the rock mass properties of the shale and sandstone formations that have been employed as input data in the numerical analysis.

Numerical Modelling
To verify the results of the empirical analysis, a 3D numerical analysis was developed. FLAC 3D commercial software was applied to obtain the deformations and the failure states considering the Appl. Sci. 2020, 10, 5853 5 of 13 excavation of the network of tunnels and the support system [9,15,35]. Figure 3 shows the geometry of the models of the lower reservoir in the form of a network of tunnels that have been selected to conduct the numerical analysis. The size of the transversal tunnels model (Figure 3a) is 82.5 m long, 12.5 m wide and 69 m high, and the central tunnel-junction zone model (Figure 3b) is 174 m long, 12.5 m wide and 80 m high. It is assumed that there are roller boundaries at the bottom and along the sides and there is an unconstrained boundary at the top of model for application of uniform vertical stress, in order to simulate the primary stress field. The mesh was refined close to the contour of the excavations and gradually was coarser at positions outwards, for increasing the accuracy of the calculations.

Numerical Modelling
To verify the results of the empirical analysis, a 3D numerical analysis was developed. FLAC 3D commercial software was applied to obtain the deformations and the failure states considering the excavation of the network of tunnels and the support system [9,15,35]. Figure 3 shows the geometry of the models of the lower reservoir in the form of a network of tunnels that have been selected to conduct the numerical analysis. The size of the transversal tunnels model (Figure 3a) is 82.5 m long, 12.5 m wide and 69 m high, and the central tunnel-junction zone model ( Figure 3b) is 174 m long, 12.5 m wide and 80 m high. It is assumed that there are roller boundaries at the bottom and along the sides and there is an unconstrained boundary at the top of model for application of uniform vertical stress, in order to simulate the primary stress field. The mesh was refined close to the contour of the excavations and gradually was coarser at positions outwards, for increasing the accuracy of the calculations.

Simulation Procedure
For the models of the network of tunnels, the solution steps included: (1) establishment of the initial stress field by applying vertical and horizontal stresses and once equilibrium was reached deformations are reset, (2) excavation of the central and transversal tunnels, and (3) installation of the support system immediately after excavation. Failure states and total displacements caused by excavations and installation of the support system were analyzed. The rock mass is considered as homogenous and linearly elastic-perfectly plastic, according to the Mohr-Coulomb (M-C) field criterion. Plastic failure occurs when the shear stress on a certain plane reaches a limit called the shear yield stress. The M-C failure criterion is expressed in Equation (5).
where σ1 and σ3 are the maximum and minimum principal stresses, respectively, in MPa, c is the cohesion, in MPa, and ϕ is the friction angle, in degrees. The initial primary stress field at the tunnels' depth of 450 m is 10.35 × 10 6 Pa in the vertical direction and 5.17 × 10 6 Pa in the horizontal direction, considering a density of 23 KN m −3 . The in-situ horizontal stresses were obtained from the empirical equations from global data and a compilation of such data within the ACCB.

Simulation Procedure
For the models of the network of tunnels, the solution steps included: (1) establishment of the initial stress field by applying vertical and horizontal stresses and once equilibrium was reached deformations are reset, (2) excavation of the central and transversal tunnels, and (3) installation of the support system immediately after excavation. Failure states and total displacements caused by excavations and installation of the support system were analyzed. The rock mass is considered as homogenous and linearly elastic-perfectly plastic, according to the Mohr-Coulomb (M-C) field criterion. Plastic failure occurs when the shear stress on a certain plane reaches a limit called the shear yield stress. The M-C failure criterion is expressed in Equation (5).
where σ 1 and σ 3 are the maximum and minimum principal stresses, respectively, in MPa, c is the cohesion, in MPa, and φ is the friction angle, in degrees. The initial primary stress field at the tunnels' depth of 450 m is 10.35 × 10 6 Pa in the vertical direction and 5.17 × 10 6 Pa in the horizontal direction, considering a density of 23 KN m −3 . The in-situ horizontal stresses were obtained from the empirical equations from global data and a compilation of such data within the ACCB. Table 3 shows the design of the support system in the transversal tunnels and the central tunnel for the existing rock masses in the study area. The bolt lengths were calculated by Equation (4), as 2.75 m;

Numerical Modelling. Support System Design
Appl. Sci. 2020, 10, 5853 6 of 13 therefore, the rock bolt length design was 3 m. In the shale formation grouted rock bolts were selected with the following specifications: 25 mm diameter, 245 KN (25 t) of load capacity, 3 m long, and spacing between bolts of 1.2 m. Furthermore, there was a layer of 150 mm thick reinforced shotcrete in the central tunnel and a layer that was 120 mm thick in transversal tunnels. In the sandstone formation, grouted rock bolts were selected with the following specifications: 25 mm in diameter, 245 kN of load capacity, 3 m long, and spacing between bolts of 1.5 m. A layer of 80 mm thick reinforced shotcrete in the central tunnel and 60 mm thick in transversal tunnels is also simulated. Table 3. Support system design.

Q-System Method
To estimate the support system it is necessary to establish an equivalent diameter for the tunnel, D e , which is the width/height of excavation divided by the ESR (excavation support ratio), which for transversal tunnels can be established at 1.3, so D e = 5.0. For the central tunnel, a value of ESR of 1.0 has been estimated, obtaining an equivalent diameter for the tunnel, D e = 6.5. The evaluations of rock masses in sandstone and shale rock masses by using the Q-system and the GSI are summarized in Table 4. By applying Equation (1), in sandstone formation a value of Q of 1.27 is obtained, while in the shale rock mass the value of Q is 0.23. The suggested support systems based on the Q-system for the central tunnel and the transversal tunnels are given in Table 5 following the recommendations proposed by Grimstad and Barton [31] using an equivalent diameter of 5.0 and 6.5 for transversal tunnel and central tunnel, respectively.

Numerical Simulations
The results of the calculations that have been conducted in sandstone and shale rock masses, in transversal tunnels, the central tunnel and the junction zone between them are shown in this section. Table 3 shows the properties of the support system used in the numerical analysis (bolt diameter, Appl. Sci. 2020, 10, 5853 7 of 13 length, spacing, load capacity, and thickness of the spayed concrete layer). Figure 4 shows the vertical and horizontal displacements in the transversal tunnels for shale rock mass according to the model indicated in Figure 3a. The maximum displacement is located at the tunnel walls (17.4 mm).

Numerical Simulations
The results of the calculations that have been conducted in sandstone and shale rock masses, in transversal tunnels, the central tunnel and the junction zone between them are shown in this section. Table 3 shows the properties of the support system used in the numerical analysis (bolt diameter, length, spacing, load capacity, and thickness of the spayed concrete layer). Figure 4 shows the vertical and horizontal displacements in the transversal tunnels for shale rock mass according to the model indicated in Figure 3a. The maximum displacement is located at the tunnel walls (17.4 mm).  Table 6 shows the summary of the numerical simulations results in transversal tunnels for shale and sandstone rock masses: vertical and horizontal displacements, thickness of the EDZ, axial load in rock bolts and axial force, bending moment and shear force in the shotcrete layer at the walls, roof and floor of tunnels for unsupported and supported cases. By installation of the support, the EDZ and total displacements (vertical and horizontal) notably decreased. Because of the lower quality of the rock mass, in shale formation the displacement values are higher than in sandstone rock mass. The maximum value of displacement reached 17.5 mm at the walls of the tunnels in shale rock mass. The maximum thickness of EDZ is located at the tunnel walls in shale formation, reaching 2.1 m. In shale formation, the vertical and horizontal displacements are reduced to 12.5 mm (27%) and 14.8 mm (15%), respectively. In sandstone formation, the vertical (roof) and horizontal (walls) displacements are reduced to 1.92 mm (34%) and 2.27 mm (22%), respectively. In addition, the thickness of the EDZ is also reduced by 38% in shale formation when the support system is applied.  Figure 5a also shows the rock bolts force for sandstone rock mass. The maximum load in the grouted rock bolts reached 30.5 kN at the tunnel walls. Figure 5b shows the failure states for  Table 6 shows the summary of the numerical simulations results in transversal tunnels for shale and sandstone rock masses: vertical and horizontal displacements, thickness of the EDZ, axial load in rock bolts and axial force, bending moment and shear force in the shotcrete layer at the walls, roof and floor of tunnels for unsupported and supported cases. By installation of the support, the EDZ and total displacements (vertical and horizontal) notably decreased. Because of the lower quality of the rock mass, in shale formation the displacement values are higher than in sandstone rock mass. The maximum value of displacement reached 17.5 mm at the walls of the tunnels in shale rock mass. The maximum thickness of EDZ is located at the tunnel walls in shale formation, reaching 2.1 m. In shale formation, the vertical and horizontal displacements are reduced to 12.5 mm (27%) and 14.8 mm (15%), respectively. In sandstone formation, the vertical (roof) and horizontal (walls) displacements are reduced to 1.92 mm (34%) and 2.27 mm (22%), respectively. In addition, the thickness of the EDZ is also reduced by 38% in shale formation when the support system is applied.  Figure 5a also shows the rock bolts force for sandstone rock mass. The maximum load in the grouted rock bolts reached 30.5 kN at the tunnel walls. Figure 5b shows the failure states for transversal tunnels in sandstone formation as the model reaches equilibrium. A combination of shear and tensile Appl. Sci. 2020, 10, 5853 8 of 13 failure initiation mechanisms are observed at the floor and walls. The failure mode changes to only shear at the roof of the tunnels. The thickness of the EDZ is 0.72 m in sandstone formation. The simulation results indicate that the designed support systems can guarantee the tunnels stability.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 13 transversal tunnels in sandstone formation as the model reaches equilibrium. A combination of shear and tensile failure initiation mechanisms are observed at the floor and walls. The failure mode changes to only shear at the roof of the tunnels. The thickness of the EDZ is 0.72 m in sandstone formation. The simulation results indicate that the designed support systems can guarantee the tunnels stability.  Figure 6a indicates the rock bolt axial force for transversal tunnels in shale formation. The maximum bolt force is located at the walls of the tunnels (130 kN). Figure 6b shows the plasticity zones for shale formation after reaching balance. Shear failure is observed at the roof, walls and floor of the transversal tunnels. The thickness of the EDZ is larger in shale rock mass, reaching a value of 2.1 m. After support installation the thickness of the EDZ is reduced to 1.3 m in shale formation and 0.37 m in sandstone formation. In addition to the stability analysis in the transversal tunnels, the state of plasticity and the total displacements have also been analyzed in the central tunnel. Figure 7 shows the axial force in the rock bolts and the state of plasticity in the central tunnel for shale formation. A combination of shear and tensile failure initiation mode is seen at the roof, walls and floor of the central tunnel. Table 7   transversal tunnels in sandstone formation as the model reaches equilibrium. A combination of shear and tensile failure initiation mechanisms are observed at the floor and walls. The failure mode changes to only shear at the roof of the tunnels. The thickness of the EDZ is 0.72 m in sandstone formation. The simulation results indicate that the designed support systems can guarantee the tunnels stability.  Figure 6a indicates the rock bolt axial force for transversal tunnels in shale formation. The maximum bolt force is located at the walls of the tunnels (130 kN). Figure 6b shows the plasticity zones for shale formation after reaching balance. Shear failure is observed at the roof, walls and floor of the transversal tunnels. The thickness of the EDZ is larger in shale rock mass, reaching a value of 2.1 m. After support installation the thickness of the EDZ is reduced to 1.3 m in shale formation and 0.37 m in sandstone formation. In addition to the stability analysis in the transversal tunnels, the state of plasticity and the total displacements have also been analyzed in the central tunnel. Figure 7 shows the axial force in the rock bolts and the state of plasticity in the central tunnel for shale formation. A combination of shear and tensile failure initiation mode is seen at the roof, walls and floor of the central tunnel. Table 7  In addition to the stability analysis in the transversal tunnels, the state of plasticity and the total displacements have also been analyzed in the central tunnel. Figure 7 shows the axial force in the rock bolts and the state of plasticity in the central tunnel for shale formation. A combination of shear and tensile failure initiation mode is seen at the roof, walls and floor of the central tunnel. Table 7 shows the summary of the numerical simulations results in the central tunnel in shale and sandstone formations for unsupported and supported cases, with the support system indicated previously in Table 3.  Finally, the stability analysis was carried out in the junction zone between the central and transversal tunnels. Figure 8a shows the horizontal displacement and Figure 8b shows the shear force in the fibre reinforced shotcrete layer.  The maximum value of displacement reached is 20.92 mm at the roof of the tunnel in shale rock mass. The thickness of the EDZ reached is 2.9 m is shale formation for the unsupported case. The vertical and horizontal displacements in sandstone formation reach 3.29 and 2.72 mm, respectively. The area of the EDZ and the maximum displacements notably decreased when the support system is applied. In shale formation, the vertical (roof) and horizontal displacements decreased down to 12.97 mm (38%) and 11.62 mm (36.5%), respectively. The rock bolt load reaches a value of 188.9 kN at the walls, while the elastic capacity of the rock bolts is 245 kN (safety factor of 1.29). The axial force, bending moment and shear force in reinforced shotcrete layer have also been analyzed, reaching 883.64 kN, 8.08 kNm, and 22.35 kN, respectively, in shale rock mass.
Finally, the stability analysis was carried out in the junction zone between the central and transversal tunnels. Figure 8a shows the horizontal displacement and Figure 8b shows the shear force in the fibre reinforced shotcrete layer.  Table 8 shows the results for the junction zone between the central and transversal tunnels. The maximum value of displacement reached 18.20 mm at the roof of the tunnels in shale rock mass. By installation of support, the thickness of the EDZ and deformations are reduced. In shale formation, the vertical and horizontal displacements decreased down to 13.84 mm at the roof and 11.51 mm in the walls, when the support system is applied. The axial force, bending moment and shear force in reinforced shotcrete layer have also been analyzed in the junction zone, reaching 1570 kN, 9.31 kNm, and 27.13 kN, respectively, in shale rock mass. The shear force indicated in Figure 8b shows peaks due to the head of the rock bolts.  Figure 9 shows the rock bolts force in the central tunnel and the transversal tunnels depending on the distance from the rock bolt head. Rock bolts located at the walls in shale and sandstone formations have been selected. The maximum axial load is reached in shale rock mass, at the central tunnel. It has been observed that in all scenarios the axial load reached is less than 245 kN (load capacity designed).  Table 8 shows the results for the junction zone between the central and transversal tunnels. The maximum value of displacement reached 18.20 mm at the roof of the tunnels in shale rock mass. By installation of support, the thickness of the EDZ and deformations are reduced. In shale formation, the vertical and horizontal displacements decreased down to 13.84 mm at the roof and 11.51 mm in the walls, when the support system is applied. The axial force, bending moment and shear force in reinforced shotcrete layer have also been analyzed in the junction zone, reaching 1570 kN, 9.31 kNm, and 27.13 kN, respectively, in shale rock mass. The shear force indicated in Figure 8b shows peaks due to the head of the rock bolts.  Figure 9 shows the rock bolts force in the central tunnel and the transversal tunnels depending on the distance from the rock bolt head. Rock bolts located at the walls in shale and sandstone formations have been selected. The maximum axial load is reached in shale rock mass, at the central tunnel. It has been observed that in all scenarios the axial load reached is less than 245 kN (load capacity designed).

Conclusions
Flexible energy storage systems allow increasing the generation of electricity by means of intermittent renewable energy sources. A closed coal mine in the Asturian Central Coal Basin (NW Spain) is proposed as a subsurface water reservoir of underground pumped storage hydropower plants. Underground pumped storage hydropower plants provide a large amount of electrical energy with rapid response and low environmental impacts. An underground reservoir conformed by a network of tunnels with an arched roof cross section of 30 m 2 has been designed. For the construction of the underground water reservoir it is necessary to analyze the stability of the excavation of the network of tunnels. Transversal tunnels, the central tunnel and junction zones have been analyzed.
Empirical analysis and three-dimensional numerical simulations have been carried out in this work. According to the empirical analysis, grouted rock bolts and a layer of fibre reinforced shotcrete were recommended as the support system. The rock mass properties used as input parameters for numerical modeling were obtained from laboratory tests and estimated from rock mass classification systems. In the numerical model, the state of plasticity, the axial load in rock bolts and the deformations around the tunnels were analyzed. In addition, the axial force, the bending moment and the shear force were also checked in the reinforced shotcrete layer. The maximum displacements and thickness of the excavation damage zone were reached in shale formation. As shown by the numerical simulations, the proposed support from empirical analysis was feasible. By applying the designed support system, the area of the excavation damage zone and the maximum displacements significantly decreased. The results of the numerical analysis show that no significant failure is expected. All this shows that a combination of empirical and numerical methods is appropriate to design a proper support system in underground infrastructures.
Based on further rock lab testing, a softening factor for the rock mass after failure should be applied in order to get more precise predictions for displacements and the extension of the excavation damage zone as well as the support loads.

Conclusions
Flexible energy storage systems allow increasing the generation of electricity by means of intermittent renewable energy sources. A closed coal mine in the Asturian Central Coal Basin (NW Spain) is proposed as a subsurface water reservoir of underground pumped storage hydropower plants. Underground pumped storage hydropower plants provide a large amount of electrical energy with rapid response and low environmental impacts. An underground reservoir conformed by a network of tunnels with an arched roof cross section of 30 m 2 has been designed. For the construction of the underground water reservoir it is necessary to analyze the stability of the excavation of the network of tunnels. Transversal tunnels, the central tunnel and junction zones have been analyzed.
Empirical analysis and three-dimensional numerical simulations have been carried out in this work. According to the empirical analysis, grouted rock bolts and a layer of fibre reinforced shotcrete were recommended as the support system. The rock mass properties used as input parameters for numerical modeling were obtained from laboratory tests and estimated from rock mass classification systems. In the numerical model, the state of plasticity, the axial load in rock bolts and the deformations around the tunnels were analyzed. In addition, the axial force, the bending moment and the shear force were also checked in the reinforced shotcrete layer. The maximum displacements and thickness of the excavation damage zone were reached in shale formation. As shown by the numerical simulations, the proposed support from empirical analysis was feasible. By applying the designed support system, the area of the excavation damage zone and the maximum displacements significantly decreased. The results of the numerical analysis show that no significant failure is expected. All this shows that a combination of empirical and numerical methods is appropriate to design a proper support system in underground infrastructures.
Based on further rock lab testing, a softening factor for the rock mass after failure should be applied in order to get more precise predictions for displacements and the extension of the excavation damage zone as well as the support loads.