An Empirical Approach for Tunnel Support Design through Q and RMi Systems in Fractured Rock Mass

: Empirical systems for the classiﬁcation of rock mass are used primarily for preliminary support design in tunneling. When applying the existing acceptable international systems for tunnel preliminary supports in high-stress environments, the tunneling quality index ( Q ) and the rock mass index ( RMi ) systems that are preferred over geomechanical classiﬁcation due to the stress characterization parameters that are incorporated into the two systems. However, these two systems are not appropriate when applied in a location where the rock is jointed and experiencing high stresses. This paper empirically extends the application of the two systems to tunnel support design in excavations in such locations. Here, the rock mass characterizations and installed support data of six tunnel projects are used. The back-calculation approach is used to determine the Q value using the Q -system support chart, and these values are then used to develop the equations and charts to characterize the stress reduction factor ( SRF ), which is also numerically evaluated. These equations and charts reveal that the SRF is a function of relative block size, strength–stress ratio, and intact rock compressive strength. Furthermore, the RMi -suggested supports were heavier than the actual installed ones. If the approximate inverse relation between stress level (SL) and SRF is used, the difference between the actual and the recommended supports increases when using the RMi -recommended rock support chart for blocky ground. An alternate system is made for support recommendation using a Q -system support chart. In this system, the ground condition factor is modiﬁed from the available parameters, and a correlation is developed with a modiﬁed Q system. analysis through thrust-bending moment interaction diagrams.


Introduction
For the purposes of rock engineering design, eight modeling methods have been categorized, including pre-existing standard methods, analytical methods, basic and extended numerical methods, precedent type analysis, empirical classification systems, the basic system approach, and the combined system approach [1]. The empirical classification systems that have been developed for rock mass are the backbone for empirical support design in tunnels, and are currently widely used. These systems have been continuously updated for this particular purpose. In 1946, Terzaghi was the first to develop a rock mass classification system, which was known as rock load factor classification, and a number of systems for rock mass quality determination have been developed since. These systems consider the most essential parameters affecting the rock mass behavior to rate its quality. In designing tunnel supports, the geomechanical classifications, which are also known as the Rock Mass Rating (RMR) system, the tunneling quality index (Q) system, and the Rock Mass index (RMi), are currently internationally accepted, and are widely used in the field of tunneling for preliminary support designs [2]. Originally, these empirical classification systems were developed, refined, and updated for the tunnel support [3]. However, they also have applications in other engineering projects such as mining [4][5][6][7], slopes [8][9][10][11][12], and tunnel-boring machine (TBM) tunnels [13,14]. Comparing the application of these systems for tunneling in high-stress environments, the Q and RMi systems are preferred over the RMR, because their parameters are more relevant for the purpose [15]. The Q-system classification was developed in 1974 [16], and major changes in its characterizations and classifications were made in 1993 and 2002 [17,18]. In characterization, the foremost changes include the modification of the stress reduction factor (SRF) and intact rock uniaxial compressive strength (σ c ), respectively. The system was initially developed for tunneling cases in hard and jointed rock mass [19], and thus, this system is more often applied in jointed rock mass [20]. However, its application in tunnel support design is limited to an equivalent dimension (D e ) ranging from 2.5 to 30, and a Q value ranging from 0.1 to 40 [21].
The Rock Mass index (RMi) system was proposed by Palmström in 1995 primarily for the reduced rock mass strength caused by jointing [22] and for rock support estimation in tunneling. The system was revised, and major developments and simplifications were made after five years of successful use and experience based on practical applications [23]. These simplifications are in terms of reduction in the input parameters and correlation between block volume (V b ) and volumetric joint count (J v ). This leads to easier and quicker calculations of the tunnel support parameters; i.e., the ground condition factor G c and size ratio factor S r , especially cases where limited field investigations have been performed.
In tunneling, in situ stresses are the key factors-besides rock mass and project-related features-that control the ground behavior, but empirical classification systems have limitations in application in high stress-environments [15,24]. The limitations of empirical classification systems, which are based on the ground behavior for underground excavation in rock engineering, are shown in Table 1. In comparing the Q system and RMi for preliminary tunnel support design in the case of high in-situ stresses, the Q system has a stress reduction factor (SRF) for the purpose, but it is limited to massive rock (Table 2) [17,25]. The RMi system uses stress level (SL) for the first three conditions of Table 2 for jointed rock mass, and has an inverse relation with SRF; it can also reveal that a high stress has a positive impact on the ground condition [22,23]. In high-stress rock mass environments, the presence of joints has a positive impact on tunneling [26], and for this reason de-stress blasting is successfully used as a rockburst preventive measure [27,28]. Neither system gives guidelines for SRF and SL in high-stressed jointed rock mass (the last three categories of Table 2), which is actually a limitation of the data on which these systems are based.  Fitness rating of different systems: 1, suitable; 2, fair; 3, poor; and 4, not applicable [15].
In this study, for the subclasses 4 and 5 of Table 2, application of the Q and RMi systems is extended from already-supported tunnel sections, which are excavated in high-stressed jointed rock mass. Equations and charts are recommended for the determination of SRF in order to extend the application of the Q system. An alternate approach of using the Q-system chart to recommend supports from a modified RMi system is taken.

Materials and Methods
From the tunneling viewpoint, the Himalayas pose the most challenging ground conditions in the world [29]. This mountain range is located in the north of the subcontinent and has great potential for hydropower generation; to harness this energy, deep tunneling here is the most tough job [30]. Along with these hydropower tunnels (completed, under construction, and planned), a number of road tunnels have also been completed recently in the northern part of Pakistan.

Project Descriptions
This section provides a brief overview of the six projects shown in Figure 1, whose data are used in this work. The first four projects are completed, and are presented in the manuscript as Project 1, 2, 3, and 4, which are based on the number of sections used in this work. The last two, Project 5 and 6, are the planned one, and their geotechnical data is used in the numerical modeling for validation. Plastic behavior (initial)  4  3-4  3  Squeezing ground  4  3  3  Raveling from slaking or friability  4  4  4  Swelling ground  4  3  3  Flowing ground  4  4  4  Water ingress  4  4  4 Fitness rating of different systems: 1, suitable; 2, fair; 3, poor; and 4, not applicable [15].
In this study, for the subclasses 4 and 5 of Table 2, application of the Q and RMi systems is extended from already-supported tunnel sections, which are excavated in high-stressed jointed rock mass. Equations and charts are recommended for the determination of SRF in order to extend the application of the Q system. An alternate approach of using the Q-system chart to recommend supports from a modified RMi system is taken.

Materials and Methods
From the tunneling viewpoint, the Himalayas pose the most challenging ground conditions in the world [29]. This mountain range is located in the north of the subcontinent and has great potential for hydropower generation; to harness this energy, deep tunneling here is the most tough job [30]. Along with these hydropower tunnels (completed, under construction, and planned), a number of road tunnels have also been completed recently in the northern part of Pakistan.

Project Descriptions
This section provides a brief overview of the six projects shown in Figure 1, whose data are used in this work. The first four projects are completed, and are presented in the manuscript as Project 1, 2, 3, and 4, which are based on the number of sections used in this work. The last two, Project 5 and 6, are the planned one, and their geotechnical data is used in the numerical modeling for validation.

Project 1
The Lowari Rail Tunnel (LRT) is 8.5 km, and is the longest transportation tunnel in Pakistan, connecting the north district, Chitral, with the remainder of the country for vehicle transportation on an all-weather link road. The cross-sectional design of the tunnel was modified by the National Highway Authority (NHA) of Pakistan after the excavation of the LRT (7.12 m in span), for a road tunnel and named the Modified Road Tunnel (MRT). This MRT (11.17 m in span) is a part of the extension of road network for connecting Central Asian countries with Pakistan through the Chitral district, traversing Vakhan (Afghanistan). The geological formation of the project area is located between the active Eurasian continental plate and the Indian plate, and due to the ongoing subduction of the Indian plate, this area belongs to the active seismic zone. Based on the geological investigations for the LRT, five major geological units were defined along the tunnel route: the granite unit, the meta sediment unit, the meta volcanic unit, the meta igneous unit, and the biotite granite unit [31]. The LRT and MRT were completed in 2009 and 2016, respectively, and became operational for traffic in 2017.

Project 2
The Neelum-Jhelum Hydroelectric Project (NJHEP) of the Water and Power Development Authority (WAPDA) of Pakistan is located in Muzaffarabad, within the Kashmir region in northeast Pakistan. A number of underground structures are part of this 969 MW run-of-the-river project, including a diversion tunnel, headrace tunnels, a tailrace tunnel, access tunnels, a powerhouse cavern, and pressure shafts. The 28.5-km headrace tunnels convey the water of the Neelum River and pass beneath the Jhelum River. These headrace tunnels are a combination of single and twin tunnels. The twin tunnels (each with a span of 9.68 m) were excavated by the conventional method and a tunnel-boring machine (TBM). These tunnels pass through the Murree Formation of alternating beds of sandstone (medium to fine-grained) and siltstone with layers of occasionally thin mudstone [30,32,33]. The sandstone (SS) is further divided into two subcategories, SS-1 and SS-2, based on their mechanical properties. SS-1 is made up of comparatively strong, well-cemented sandstone, of fine to medium grain. It is also thinly to thickly bedded, massive at places, blocky, moderately to closely jointed, and at places fractured [34].

Project 3
The Kohat Tunnel (KT)-also known the Pak-Japan Friendship Tunnel or simply the Friendship Tunnel-is an alternate route for traversing Kohat Pass in the Khyber Pakhtunkhwa province of Pakistan along the major national highway of NHA (known as the Indus highway), and follows the west bank of the Indus River. The 1.9-km, two-lane road tunnel was the first long road tunnel (9.4-m wide) in the country, and was open to traffic in 2003. Difficult geological sites were encountered during excavation. Along the tunnel route, the major rock types were limestone and shale [35].

Project 4
The Golen Gol Hydropower Project (GGHP) is the most recently completed 108-MW run-of-the-river project of the WAPDA. The project is located in the Chitral district, Khyber Pakhtunkhwa province, Pakistan. A number of underground structures, including a headrace tunnel, a pressure shaft, and a pressure tunnel, are the major components of the project. The total length of the 3.7-m wide single headrace tunnel for this project is 3.8 km. The maximum longitudinal section of the headrace tunnel passes through granite rock (2.8 m), followed by quartz mica schist, marble, and calcareous quartzite [36].

Project 5
Located at the boundary of the Khyber Pakhtunkhwa and Gilgit Baltistan provinces of Pakistan, the multipurpose planned Diamer Basha Dam project has a capacity of 4500 megawatts of hydroelectric power. The project is within a location where the subduction of the Indian plate beneath the Eurasian plate occurs and major rock types in this location is GN (gabbronorite) and UMA (ultramafic association) [37][38][39]. The planned 15.4-m span diversion tunnels on the right bank of the Indus River will be the first underground water passageway during the project construction.

Project 6
The planned Kohala hydropower project, a run-of-the-river project, is located on the Jhelum River in the Muzaffarabad district, Kashmir, Pakistan. The project includes a 17.5-km-long, 8.5-m span headrace tunnel for the generation of 1100 MW of electric power. Since the project is located in the same geological formation as Project 2, and also from local exploration, the expected major rock units will be SS-1 and SS-2 along the tunnel route [37]. The in situ stress measurement results show that high in situ stresses will be the major factor controlling the rock mass behavior [40].

Characterization of Rock Mass Based on the Q System
The tunneling data-based Q system categorizes rock mass into nine classes on a logarithmic scale, ranging in value from 0.001 to 1000 (exceptionally poor to exceptionally good, respectively). It is calculated through Equation (1) [16]: where RQD is the rock quality designation, J n denotes the rating of the number of joint sets, J r is the rating for the joint surface roughness, J a denotes the rating for the degree of alteration or clay-filling joint set, J w denotes the ratings for groundwater inflow and pressure effects, and SRF is the stress reduction factor. The ratio RQD/J n is the relative block size in a selected domain, J r /J a is a measure of the inter-block frictional strength in that domain, and J w /SRF is the effect of water, the strength-stress ratio, faulting, swelling or squeezing, and defined active stress. The role of σ c is substantial in the properties of rock mass. Therefore, a normalization factor was applied to Equation (1) and modified to Q c (Equation (2)) [18]: The New Austrian Tunneling Method (NATM) approach, which considers the ground behavior, was used for the excavation and support of the tunnels in the first four projects. This approach also included the geological documentation in the construction of conventional tunneling [41]. Therefore, the available data from the excavation face mapping is appropriate for the characterization of the RQD, J n , J r , J a , J w , and σ c .
Rockburst is a dynamic phenomenon that is linked with highly stressed environments in underground excavations. It is dominant in massive rock, and excavation in jointed rock mass has a positive impact under the same stress state as compared to massive rock. The increase in SRF that is shown in Table 2 is for massive rock based on RQD/J n [17]. In jointed rock mass, the extreme value of the old SRF (1974 version) is proposed [42]. Equation (3) [43] and Equation (4) [42] are based on mining cases in South Africa and Australia, respectively, for the characterization of SRF (1974 version) for competent rock having rock stress problems.
The average unit weight of the intact rock for maximum sections of the first four projects in this study is about 0.027 MN/m 3 . Equation (3) can be rewritten as Equation (5): Equation (6) was developed for the calculation of SRF in the Q system using the data of the Nathpa Jhakri project in India [44]: Equation (6) can be rewritten, with 0.027 MN/m 3 as the unit weight of the rock, as Equation (7): Barton [18] divided the rock mass into massive and jointed, based on its RQD/J n ratio, and revealed that a massive rock can have a ratio from 25 to 200, while a jointed rock mass has a typical ratio of 10. According to the rock mass fabric index approach, the maximum limit of the RQD/J n ratio for a jointed rock mass is 16.5 [45]. The relationship between the RQD/J n ratio and SRF in hard rock under high stress directs that, for jointed rock mass, the SRF value that was suggested for massive rock in Table 2 is too high [17]. The Q-system parameters were characterized for those selected sections where the RQD/J n ratio <16.5, and the intact rock strength to the major principal stress ratio (σ c /σ 1 ) is in the range of two to five.
The characterization parameters of Equation (2) for the first four projects are tabulated in Table 3, and details of RQD/J n and σ c /σ 1 are displayed in Figure 2. In Figure 2a, the strength-stress ratio is in the range of two to five for Project 1 and 2. The data statistics show that all of the sections of Project 3 experienced a strength-stress ratio higher than 3.5. In the case of Project 4, this ratio is less than 3.5. In Figure 2b, the relative block sizes in the majority of the sections are in the range of six to eight.
Equation (6) was developed for the calculation of SRF in the Q system using the data of the Nathpa Jhakri project in India [44]: Equation (6) can be rewritten, with 0.027 MN/m 3 as the unit weight of the rock, as Equation (7): Barton [18] divided the rock mass into massive and jointed, based on its RQD/Jn ratio, and revealed that a massive rock can have a ratio from 25 to 200, while a jointed rock mass has a typical ratio of 10. According to the rock mass fabric index approach, the maximum limit of the RQD/Jn ratio for a jointed rock mass is 16.5 [45]. The relationship between the RQD/Jn ratio and SRF in hard rock under high stress directs that, for jointed rock mass, the SRF value that was suggested for massive rock in Table 2 is too high [17]. The Q-system parameters were characterized for those selected sections where the RQD/Jn ratio <16.5, and the intact rock strength to the major principal stress ratio (σc/σ1) is in the range of two to five.
The characterization parameters of Equation (2) for the first four projects are tabulated in Table  3, and details of RQD/Jn and σc/σ1 are displayed in Figure 2. In Figure 2a, the strength-stress ratio is in the range of two to five for Project 1 and 2. The data statistics show that all of the sections of Project 3 experienced a strength-stress ratio higher than 3.5. In the case of Project 4, this ratio is less than 3.5. In Figure 2b, the relative block sizes in the majority of the sections are in the range of six to eight.

SRF Calculation Through Back Analysis
A back-analysis procedure for the Q-value calculation of rock mass, as shown in Figure 3, was suggested using the installed rock bolts and fiber-reinforced, sprayed shotcrete [46]. The back-calculation is explained through Example 1 for the calculation of the Q value. The latest support chart of the Q system [19] is used for the back-calculation for all of the sections of the four projects. The support chart indicated that the shotcrete thickness is a function of D e , in addition to the rock quality (Q or Qc); however, rock bolt spacing (S b ) is a function of rock quality only. This spacing of rock bolts was used for the back-calculation of the tunneling quality index (Q 1 ). The thickness of the fiber-reinforced shotcrete is a function of the tunnel span and the tunneling quality index. The thickness of the fiber-reinforced shotcrete, along with the tunnel span, was also used for the back-calculation of the tunneling quality index (Q 2 ). The D e is the ratio of tunnel width/span and the excavation support ratio (ESR). The ESR value depends on the country safety standards [20] as well as on the quality of the rock mass [19]. The relation of rock bolt length (L b ) to D e for hydropower tunnels, i.e., Project 2 and Project 4, reveals that a value of one for ESR is an appropriate rating for headrace tunnels, and was also used in this study.

SRF Calculation Through Back Analysis
A back-analysis procedure for the Q-value calculation of rock mass, as shown in Figure 3, was suggested using the installed rock bolts and fiber-reinforced, sprayed shotcrete [46]. The backcalculation is explained through Example 1 for the calculation of the Q value. The latest support chart of the Q system [19] is used for the back-calculation for all of the sections of the four projects. The support chart indicated that the shotcrete thickness is a function of De, in addition to the rock quality (Q or Qc); however, rock bolt spacing (Sb) is a function of rock quality only. This spacing of rock bolts was used for the back-calculation of the tunneling quality index (Q1). The thickness of the fiberreinforced shotcrete is a function of the tunnel span and the tunneling quality index. The thickness of the fiber-reinforced shotcrete, along with the tunnel span, was also used for the back-calculation of the tunneling quality index (Q2). The De is the ratio of tunnel width/span and the excavation support ratio (ESR). The ESR value depends on the country safety standards [20] as well as on the quality of the rock mass [19]. The relation of rock bolt length (Lb) to De for hydropower tunnels, i.e., Project 2 and Project 4, reveals that a value of one for ESR is an appropriate rating for headrace tunnels, and was also used in this study.

Example 1
Tunnel span or De (m) = 11.17 Rock bolt spacing (m) = 1.6 Shotcrete thickness (cm) = 9 Q1 = 0.7 Q2 = 1.7 Q or Qc = 1.03 Although the modified Q value (Qc) was due to the additional normalized factor in Equation (1), the support chart remained unchanged. The resultant tunneling quality index from back analysis was defined as Q or Qc. Equations (1) and (2) were rearranged to equations (8) and (9), respectively, and two different SRFs (SRFQ or SRFQC) were calculated, taking the back-calculated rock quality value of either Q or Qc, respectively.

Example 1
Tunnel span or D e (m) = 11.17 Rock bolt spacing (m) = 1.6 Shotcrete thickness (cm) = 9 Although the modified Q value (Q c ) was due to the additional normalized factor in Equation (1), the support chart remained unchanged. The resultant tunneling quality index from back analysis was defined as Q or Q c . Equations (1) and (2) were rearranged to Equations (8) and (9), respectively, and two different SRFs (SRF Q or SRF QC ) were calculated, taking the back-calculated rock quality value of either Q or Q c , respectively.

Rock Mass Characterization by RMi System
The RMi system is a relative strength index for rock masses that numerically expresses the quality of the rock mass. The principle of this system is based on the reduction of σ c due to the penetrated joints. The RMi system indicates the approximate uniaxial compressive strength of a rock mass by combining σ c and a jointing parameter (JP) using Equation (10): where: σ c = Intact rock uniaxial compressive strength JP = Joint parameter (a parameter that combines the main joint features in the rock mass) jC = Joint condition factor (Equation (12)) The superscript D in Equation (10) is defined through Equation (11): In the RMi system, JP is a function of V b and jC. The jC can be estimated by joint roughness (jR), joint alteration (jA), and joint size/length (jL) using Equation (12):

Tunnel Support Design by RMi System in Discontinuous Ground
By definition, RMi is the uniaxial compressive strength (UCS) of rock mass, which only includes the inherent features affecting its strength. RMi has various applications in rock engineering, including tunnel support design. After the effective use of RMi in the tunneling field, some significant changes and adjustments for estimating preliminary rock supports by using the block volume and tunnel diameter were made [23]. The stability in blocky (jointed) ground is mainly influenced by the block size and shape, the shear strength of the joints delineating the blocks, and the orientation of the same joints relative to the tunnel axis. For the support application, RMi and tunnel span were adjusted based on the ground condition factor (G c ) and the size ratio factor (S r ), respectively, using Equations (13) and (14), as shown in Figure 4. The support parameters, which include all of these features, are used in the support chart. Here: SL = Adjustment parameter for stress level C = Gravity adjustment parameter for support in the roof or wall CF = Continuity factor, which is the ratio of the tunnel diameter (D t ) and the block diameter (D b ), differentiating between continuous and discontinuous ground and discontinuous ground having CF in the range from five to 100 C 0 = Adjustment for joint orientation N j = Adjustment for number of joint sets  . Flowchart for the tunnel support estimation using the RMi system in discontinuous ground.
The RMi system uses more parameters than the Q system, and it is relatively more complicated to use than the Q system. According to Palmström and Stille [15], the groundwater factor is not properly included in the RMi system.
In this study, the rating for jL is equal to one, because the persistence of joints is from three to 10 meters. According to the Q system [19], support at the wall and roof is independent of its location in a high-stress environment; therefore, the rating for the gravity adjustment factor (C) is also one. The ratings for joint roughness and joint alteration show that there are four different cases based on the Jr/Ja and jR/jA values, as shown in Table 4. The RMi system uses block volume, which is formed by joints in a three-dimensional network in an underground environment. This parameter is usually determined on a two-dimensional surface or by a one-dimensional borehole or scan line. Only a limited part of a joint can be correctly measured in a given location, and the direct characterization of each block size is not practical at the excavation face. Since it can be difficult to recognize the actual size and shape of individual blocks, the block size and shape must sometimes be determined from reasonable simplifications. Indirect approaches are used for the determination of block size. The generalized correlation between the block volume (Vb) and the volumetric joint count (Jv) is given in Equation (15), where block shape influences the correlation [47]: Rock support chart for blocky/discontinuous ground Figure 4. Flowchart for the tunnel support estimation using the RMi system in discontinuous ground.
The RMi system uses more parameters than the Q system, and it is relatively more complicated to use than the Q system. According to Palmström and Stille [15], the groundwater factor is not properly included in the RMi system.
In this study, the rating for jL is equal to one, because the persistence of joints is from three to 10 m. According to the Q system [19], support at the wall and roof is independent of its location in a high-stress environment; therefore, the rating for the gravity adjustment factor (C) is also one. The ratings for joint roughness and joint alteration show that there are four different cases based on the J r /J a and jR/jA values, as shown in Table 4. The RMi system uses block volume, which is formed by joints in a three-dimensional network in an underground environment. This parameter is usually determined on a two-dimensional surface or by a one-dimensional borehole or scan line. Only a limited part of a joint can be correctly measured in a given location, and the direct characterization of each block size is not practical at the excavation face.
Since it can be difficult to recognize the actual size and shape of individual blocks, the block size and shape must sometimes be determined from reasonable simplifications. Indirect approaches are used for the determination of block size. The generalized correlation between the block volume (V b ) and the volumetric joint count (J v ) is given in Equation (15), where block shape influences the correlation [47]: where β is the block shape factor, which depends on the ratio of the longest and shortest dimension of the block; the common value is 36. J v is the number of joints intersecting a volume of one m 3 . Based on J v , the degree of jointing is classified into six classes, from a very low degree of jointing (J v < 1) to crushed rock (J v > 60). Random joints may represent a significant part of the total number of discontinuities, and neglecting them would lead to the erroneous quantification of the discontinuity nature of the rock mass [48]. A rule of thumb for the calculation of J v is given in Equation (16) [49]: where S 1 , S 2 , and S 3 are the average spacings for the joint sets, N r is the number of random joints in the actual location, and A is the area in m 2 . It is not possible to measure all of the joints in a unit volume. An alternate method to find J v is from RQD or the joint frequency, as given in Equations (17) and (18) [22,47]: where K 1 = 2 for the average condition, and N 1 is the number of joints per unit length (joint frequency). Joint frequency can be calculated directly or determined indirectly. The first indirect methodology used for fracture frequency is the mathematical equation between RQD and fracture frequency [50], which is given as follows: where λ is the joint frequency. Another method to determine the joint frequency is to use RQD and the joint spacing rating, by using Equation (20) [51] as follows: where N f is the fracture frequency rating. The joint frequency can be calculated from N f using a graph of rating for the discontinuity density [52]. In this case, when the data are available both in terms of the number of joint sets and the average spacing between the joints of each joint set, the joint frequency is determined using Equation (21) as follows [46]: λ = Number of joint sets Average spacing between joints (m) .
The details of the ground conditions based on the RMi system for all of the sections of the four projects are summarized in Table 5. Table 5. Ground conditions (G c ) and size ratios (S r ) for all of the tunnel sections of the four tunnels.

Back-Calculated SRF Comparison with Previous Work
Equations (3) through (7) are for the SRF characterization of Equation (1); therefore, the back-calculated SRF Q values are compared with them. For the given values, Equation (4) gives a higher value of SRF than Equation (5). This exponential Equation (4) agreed with Barton for the calculation of SRF in highly stressed jointed rock mass [53]. The second part of Equation (5) has a negligible effect on the SRF value for σ c /σ 1 < 5. SRF was calculated from Equation (4) for different field stress ratios, and the details of the differences compared with back-calculated SRF (SRF Q ) values are shown in Figure 5.

Back-Calculated SRF Comparison with Previous Work
Equations (3) through (7) are for the SRF characterization of Equation (1); therefore, the backcalculated SRFQ values are compared with them. For the given values, Equation (4) gives a higher value of SRF than Equation (5). This exponential Equation (4) agreed with Barton for the calculation of SRF in highly stressed jointed rock mass [53]. The second part of Equation (5) has a negligible effect on the SRF value for σc/σ1 <5. SRF was calculated from Equation (4) for different field stress ratios, and the details of the differences compared with back-calculated SRF (SRFQ) values are shown in Figure 5. Due to the superscript value (0.001) in Equation (7), the variation of SRF is negligible with σc/σ1.

Suggested SRF for Highly Stressed Jointed Rock Mass
The back-analysis approach is used for the tunneling quality index calculation, using the procedure shown in Figure 3, and taking this value as Q. The SRF calculated from Equation (8)  Due to the superscript value (0.001) in Equation (7), the variation of SRF is negligible with σ c /σ 1 .

Suggested SRF for Highly Stressed Jointed Rock Mass
The back-analysis approach is used for the tunneling quality index calculation, using the procedure shown in Figure 3, and taking this value as Q. The SRF calculated from Equation (8) is defined as SRF Q for all of the sections. The SRF Q values are plotted against RQD/J n for different σ c /σ 1 values, as shown in Figure 6a. These values are also plotted against σ c /σ 1 for different RQD/J n values, as shown in Figure 6b. A strong relationship was found between SRF Q and RQD/J n and σ c /σ 1 . The relations between SRF Q , RQD/J n , and σ c /σ 1 were used, and Equation (22)  The tunneling quality index obtained from the back analysis is taken as Qc, and the SRF value calculated from Equation (9) is taken as SRFQC. This SRFQC value can also be obtained from Equation (23), using the normalized factor: The resultant SRFQC values calculated from Equation (9) or Equation (23) were plotted against RQD/Jn for different σc values and ranges of σc/σ1. The variation in SRFQC with RQD/Jn is shown in Figure 7. The tunneling quality index obtained from the back analysis is taken as Qc, and the SRF value calculated from Equation (9) is taken as SRF QC . This SRF QC value can also be obtained from Equation (23), using the normalized factor: The resultant SRF QC values calculated from Equation (9) or Equation (23) were plotted against RQD/J n for different σ c values and ranges of σ c /σ 1 . The variation in SRF QC with RQD/J n is shown in Figure 7.  The tunneling quality index obtained from the back analysis is taken as Qc, and the SRF value calculated from Equation (9) is taken as SRFQC. This SRFQC value can also be obtained from Equation (23), using the normalized factor: The resultant SRFQC values calculated from Equation (9) or Equation (23) were plotted against RQD/Jn for different σc values and ranges of σc/σ1. The variation in SRFQC with RQD/Jn is shown in Figure 7.

Comparison of the Actual Installed Supports against the RMi-Suggested Supports
In the RMi system, the suggested rock support chart that was developed for blocky ground is based on CF (5< CF <100); therefore, only the sections that satisfy these limitations are considered here for rock support determination using the RMi system. Since the minimum value of CF is 7.84 in Table 5, cases with CF >100 are omitted. After applying these limitations, details of the Jv classification are given in Figure 8. Comparing the actual installed supports with those suggested by RMi, out of the 49 sections of Project 3, only a single section showed that the actual rock bolt spacing was closer (5.88%) than the RMi-suggested spacing. In Project 4, the RMi system suggested rock bolt spacing of 2.5 m for five sections, but the actual spacing was 1.7 m. In Project 2, five sections were supported with rock bolts

Comparison of the Actual Installed Supports against the RMi-Suggested Supports
In the RMi system, the suggested rock support chart that was developed for blocky ground is based on CF (5 < CF < 100); therefore, only the sections that satisfy these limitations are considered here for rock support determination using the RMi system. Since the minimum value of CF is 7.84 in Table 5, cases with CF >100 are omitted. After applying these limitations, details of the J v classification are given in Figure 8.

Comparison of the Actual Installed Supports against the RMi-Suggested Supports
In the RMi system, the suggested rock support chart that was developed for blocky ground is based on CF (5< CF <100); therefore, only the sections that satisfy these limitations are considered here for rock support determination using the RMi system. Since the minimum value of CF is 7.84 in Table 5, cases with CF >100 are omitted. After applying these limitations, details of the Jv classification are given in Figure 8. Comparing the actual installed supports with those suggested by RMi, out of the 49 sections of Project 3, only a single section showed that the actual rock bolt spacing was closer (5.88%) than the RMi-suggested spacing. In Project 4, the RMi system suggested rock bolt spacing of 2.5 m for five sections, but the actual spacing was 1.7 m. In Project 2, five sections were supported with rock bolts Comparing the actual installed supports with those suggested by RMi, out of the 49 sections of Project 3, only a single section showed that the actual rock bolt spacing was closer (5.88%) than the RMi-suggested spacing. In Project 4, the RMi system suggested rock bolt spacing of 2.5 m for five sections, but the actual spacing was 1.7 m. In Project 2, five sections were supported with rock bolts spaced at 1.4 m and shotcrete that was 12 cm thick, but the RMi system suggested a special design for these sections. The details of all of the sections of the four projects in terms of the difference (%) in rock bolt spacing and shotcrete thickness are shown in Figure 9. spaced at 1.4 m and shotcrete that was 12 cm thick, but the RMi system suggested a special design for these sections. The details of all of the sections of the four projects in terms of the difference (%) in rock bolt spacing and shotcrete thickness are shown in Figure 9.

Tunnel Support Design by RMi System on Discontinuous Ground Using the Q-System Support Chart
As shown in Table 6, SL ≈ 1/SRF [22]. If this relation for SL characterization in highly stressed, jointed rock mass is used, the already heavier supports suggested by the RMi system would be further increased. Since the SRF calculation was based on the Q-system support chart, an attempted has been made to use it for the RMi system, again in this particular case of a highly stressed, jointed rock mass. The modified Q system from Equation (2) can be rewritten as Equation (24) using c = Jr/Ja. Similarly, the ground condition factors in Equation (13) can be rewritten as Equation (25) by introducing Nj and C0 into it from Equation (14), and replacing Vb in terms of RQD using Equation (17). The constants a and b depend on jR/jA: where SRF = f (σc, σc/σ1 and RQD/Jn), and the constants a, b, and c are given in Table 7, depending on the case number as in Table 4. Table 7. Constants a, b, and c for equations (24) and (25)   . Actual and RMi-suggested support differences in term of (a) rock bolt spacing (b) shotcrete thickness (positive value on the x-axis indicates that the recommended support is heavier than the actual).

Tunnel Support Design by RMi System on Discontinuous Ground Using the Q-System Support Chart
As shown in Table 6, SL ≈ 1/SRF [22]. If this relation for SL characterization in highly stressed, jointed rock mass is used, the already heavier supports suggested by the RMi system would be further increased. Since the SRF calculation was based on the Q-system support chart, an attempted has been made to use it for the RMi system, again in this particular case of a highly stressed, jointed rock mass. The modified Q system from Equation (2) can be rewritten as Equation (24) using c = J r /J a . Similarly, the ground condition factors in Equation (13) can be rewritten as Equation (25) by introducing N j and C 0 into it from Equation (14), and replacing V b in terms of RQD using Equation (17). The constants a and b depend on jR/jA: where SRF = f (σ c , σ c /σ 1 and RQD/J n ), and the constants a, b, and c are given in Table 7, depending on the case number as in Table 4.
As J w and C 0 are not properly included in the RMi and Q c systems, respectively [15,18], taking J w = one (dry condition) and C 0 = one (favorable joint orientation), the above two equations are based on the same parameters as the relation of N j , and the number of joint sets is given in Equation (26) [23]: where n j is the number of joint sets. The ratio between G * c and Q c can be developed from the available data using Equations (24) and (25), keeping the limitations of the data of the four projects in this study.
Knowing the ratio, Table 8 can be used to adjust G * c for joint orientation and corresponding Q c . This Q c value is then adjusted for J w and used for the recommended support from the Q-system support chart. Ultimately, the groundwater factor is included in the RMi system and C 0 is included in the Q system.

Numerical Modeling
The results from seven boreholes (BDR-08, 10, 21, 22, 24, 25, and 26) along the diversion tunnels of Project 5 show that 75% of the RQD/J n values were less than 16.5, and only 7.6% of the RQD/J n values were greater than 25. Similarly, for Project 6, of 13 boreholes (BH-01, 02, 03, 04, 05, 08, 09, 10,11,12,13,15,26), in most of them, the degree of jointing was high to very high, resulting in RQD < 50%. Only BH-09 and BH-26 in Project 6 had predominantly moderate jointing, which corresponds to a range of RQD from 50% to 75% [37]. For simplicity and to reduce the number of cases that were numerically evaluated, J r = 3, J a = 1, and J w = 1 were used for empirical support design.
The average properties of the intact rock for different rock types of projects 5 and 6 are given in Table 9. For the numerical validation of Equations (22) and (23), different overburdens were used for different scenarios of σ c /σ 1 , within the limitations of this study (5 > σ c /σ 1 > 2). For different RQD/J n and σ c /σ 1 values, SRF values were determined using Equations (22) and (23). Based on the calculated tunneling quality index (Q or Qc) from Equations (1) and (2), preliminary supports were determined from the support chart of the Q system; the details are shown in Table 10.  Table 11 shows the rock mass properties that were extrapolated from the intact rock strength properties of Table 9, and the geological strength index (GSI) values through the generalized Hoek-Brown failure criterion [54]. The GSI values that are shown in Table 11 are for the particular cases of Table 10, and are based on the rock mass fabric index approach, which was calculated from RQD/J n and J r /J a values [45]. The empirically determined support system was used in numerical modeling to analyze the results before and after the support in terms of the plastic zones in the vicinity of the tunnels and displacement in rock mass. The total displacement at different points along the tunnel perimeter ( Figure 10) are shown in Table 12. For the analysis, FLAC (Fast Lagrangian Analysis of Continua) version 7.0 is used. FLAC is an explicit two-dimensional (2D) finite difference program that is suited for sequential excavation modeling. A horseshoe-shaped tunnel was selected because it is a better shape for the drilling and blasting methods that are used in such circumstances, and is the most commonly used tunnel shape in Pakistan.
The two excavation stages (top heading and bench) and three construction steps (excavation, applying soft shotcrete and rock bolt, and hard shotcrete) in each stage are simulated for analysis. The model dimensions are 90 m × 80 m for both projects 5 and 6. The modified Hoek-Brown model was used for the analysis. This model is based on the nonlinear relation between major and minor principal stresses, σ 1 and σ 3 , and the criterion is used for plastic yielding when σ 3 is compressive. The fine mesh was simulated around the tunnel boundary for better results. The model was fixed at the sides and bottom, and vertical stresses (σ yy ) were applied at the top of the model. The in situ stress environment was created using gravity, σ yy , and the FISH (a programming language embedded within FLAC that enables the user to define new variables and functions) function. During the three construction steps for each excavation stage, 40% 30%, and 30% relaxation were used. To calculate the total displacement at points along the tunnel perimeter, the x and y displacement components were obtained from the grid points data. The plastic zones before and after support were analyzed for all four cases of the two projects. The plastic state zone was retrieved from the plot command using the zone condition in the model. For the SS-1 rock case of Project 6, the yield in shear was still the comparative maximum of all the cases, as shown in Figure 10; however, with the application of support, this decreased enormously.
Capacity diagrams, which graphically represent the relationships between thrust-bending moments and thrust-shear forces, are used to evaluate the stability of liners [55]. For this, the structure's liner information was obtained from the utility function command. Based on the proposed shotcrete thickness from Table 10, thrust-bending diagrams were plotted ( Figure 11) for shotcrete thicknesses of 12 cm and 9 cm for factor of safety (FOS) values equal to 1.5, 2, and 2.5. the plot command using the zone condition in the model. For the SS-1 rock case of Project 6, the yield in shear was still the comparative maximum of all the cases, as shown in Figure 10; however, with the application of support, this decreased enormously. Capacity diagrams, which graphically represent the relationships between thrust-bending moments and thrust-shear forces, are used to evaluate the stability of liners [55]. For this, the structure's liner information was obtained from the utility function command. Based on the proposed shotcrete thickness from Table 10, thrust-bending diagrams were plotted ( Figure 11) for shotcrete thicknesses of 12 cm and 9 cm for factor of safety (FOS) values equal to 1.5, 2, and 2.5. (a) 9-cm-thick shotcrete (b) 12-cm-thick shotcrete Figure 11. Shotcrete stability analysis through thrust-bending moment interaction diagrams.

Discussion
The application of rock mass classification in tunnel support design is now acceptable worldwide, and a number of empirical systems have been developed since the first classification system was developed by Terzaghi in 1946, which was known as rock load classification. The application of such systems is not universal, due to the case histories on which each system is based. However, a few systems, including RMR, Q, and RMi, have gotten international attention. Compared to the other classification systems, these systems have still been widely used in tunnel design since the plot command using the zone condition in the model. For the SS-1 rock case of Project 6, the yield in shear was still the comparative maximum of all the cases, as shown in Figure 10; however, with the application of support, this decreased enormously. Capacity diagrams, which graphically represent the relationships between thrust-bending moments and thrust-shear forces, are used to evaluate the stability of liners [55]. For this, the structure's liner information was obtained from the utility function command. Based on the proposed shotcrete thickness from Table 10, thrust-bending diagrams were plotted ( Figure 11) for shotcrete thicknesses of 12 cm and 9 cm for factor of safety (FOS) values equal to 1.5, 2, and 2.5. (a) 9-cm-thick shotcrete (b) 12-cm-thick shotcrete Figure 11. Shotcrete stability analysis through thrust-bending moment interaction diagrams.

Discussion
The application of rock mass classification in tunnel support design is now acceptable worldwide, and a number of empirical systems have been developed since the first classification system was developed by Terzaghi in 1946, which was known as rock load classification. The application of such systems is not universal, due to the case histories on which each system is based. However, a few systems, including RMR, Q, and RMi, have gotten international attention. Compared to the other classification systems, these systems have still been widely used in tunnel design since Figure 11. Shotcrete stability analysis through thrust-bending moment interaction diagrams.

Discussion
The application of rock mass classification in tunnel support design is now acceptable worldwide, and a number of empirical systems have been developed since the first classification system was developed by Terzaghi in 1946, which was known as rock load classification. The application of such systems is not universal, due to the case histories on which each system is based. However, a few systems, including RMR, Q, and RMi, have gotten international attention. Compared to the other classification systems, these systems have still been widely used in tunnel design since their development, because they are continuously upgraded in terms of characterizations or support, or both [3]. The evaluations of the three systems indicate that the Q and RMi systems incorporate stress-related factors, making their application more suitable for tunnel support design, where stress is the major factor governing the ground behavior. However, these systems do not cover tunnels designed in jointed rock mass and high in situ stress fields. The two systems are limited by the case histories on which they are based. In this study, the application of these two systems are empirically extended from the relevant tunnel case records.
In the literature, the SRF factor is either a function of RQD/J n or σ c /σ 1 . In this study, the empirical equation that was developed for the characterization of SRF is a function of both RQD/J n and σ c /σ 1 , along with σ c . Figure 6a shows that SRF Q (the SRF value that was developed for the original Q-system equation, Equation (1)) significantly increases with relative block size. A maximum change in RQD/J n of less than seven results in a variation of 5.4 for SRF Q , while for a maximum change in RQD/J n greater than seven, SRF Q varies by approximately 15.05, which is larger by a factor of about 2.75. A total SRF Q change of 20.82 can be observed with RQD/J n for a given value of σ c /σ 1 . With an increase in RDQ/J n , the difference of SRF Q for different values of σ c /σ 1 is comparatively low for a given value of RQD/J n .
An opposite trend can be seen in Figure 6b. For RQD/J n , the rate of decrease of SRF Q with respect to low values of σ c /σ 1 is large compared to high values of σ c /σ 1 . For a σ c /σ 1 value below 3.5, the maximum fluctuation in σ c /σ 1 for a given value of RQD/J n results in a variation of the SRF Q value of about 2.49. Fluctuation in σ c /σ 1 by an equal amount in its upper limit (above 3.5) can change SRF Q by 1.35 approximately, which is about 1.84 times smaller. The maximum change in SRF Q is 4.09 with σ c /σ 1 for any given value of RQD/J n . The available data indicate that the maximum value of the relative block size is 12. Figure 7 shows that the variation of SRF QC with RQD/J n also depends on σ c as well as σ c /σ 1 . For the same value of σ c /σ 1 , high principal stress (σ 1 ) requires rock with a high σ c value. The high level of stress means that more support is required, and hence, a high value of SRF QC . The maximum variation of SRF QC with RQD/J n is 18.5 for σ c = 100 MPa, with σ c /σ 1 ranging from three to four (Figure 7b), and the minimum variation is 7.38 for σ c = 37.5 MPa (Figure 7c). During the design stage of a tunnel, when inadequate information exists, then, with knowledge of σ c (which is the most basic parameter to be determined) and the location of the tunnel, Figure 7 can be used to approximate values of SRF QC for different values of RQD/J n . In the construction phase of a tunnel, when most of the information is in hand at the excavation face, Equations (22) and (23) can be used to calculate SRF QC . The maximum rating of σ c is 100, as shown in Table 3; therefore, the calculated SRF Q and SRF QC have the same value for tunnel sections where σ c = 100 MPa. In the remaining sections, SRF QC < SRF Q , and the difference depends on σ c .
The RMi system uses SL in its structure for the stress factor (Equation (13)), and its ratings are shown in Table 6. Comparing Table 6 with Table 2, SL is only applicable for the first three conditions of Table 2 where σ c /σ 1 > 5. Equation (13) shows that a high value of SL is favorable, but it is not applicable when σ c /σ 1 < 5, as shown in Table 2. The comparison of actual installed supports with those supports that are suggested by RMi showed that the suggested support was heavier than the actual. The equations for the calculation of Q c and G c are rewritten in such a way that each equation is the function of the same parameter, and a relation is developed between the two, keeping in view the limitation of data.
The extension of the two systems is based on Equations (22) and (23); therefore, these equations are used for determining the empirical supports of four rock units, using the Q-system support chart. The support application is numerically evaluated in terms of a decrease in the total displacement at different points along the tunnel periphery, and in the plastic zone. The maximum decrease in total displacement occurred at the tunnel crowns, with a value of 1.3 mm for the SS-2 rock unit of Project 6 ( Table 12). According to plasticity theory, a plastic zone occurs around a tunnel after excavation when the induced stresses exceed the rock mass strength. The application of empirically recommended support not only reduced the number of yield elements, but also substantially decreased the extent of the plastic zone in all of the cases. With support, yield in tension was not observed, while yield in shear was either zero or substantially reduced.
As can be seen from the capacity diagram in Figure 11, none of the cases cross the limit of FOS = 1.5. The capacity diagram for the liner shows that for the SS-2 rock unit of Project 6, the dominant failure mode was in compression; however, it too does not cross FOS = 1.5. In the case of the GN rock unit of Project 5, the dominant failure mode was in tension, but FOS > 1.5.

Conclusions
RMR, Q, and RMi classifications are based on jointed rock mass data, and are extensively used in the support design of tunnels, but are not used in highly stressed, jointed environments due to the limitations of the relevant data on which these systems are based. As a result, applications of these systems are limited under such conditions.
The Q system does not provide SRF values for highly stressed, jointed rock mass. To extend the Q-system application for such a situation, empirical equations are suggested for SRF characterization in competent rock that has rock stress problems. For jointed rock mass, the relevant data from the four existing tunnel projects that already have supports installed are used. In the proposed equations, SRF depends upon the intact rock strength, the ratio of intact rock strength to the major principal stress, and the relative block size.
The stress level (SL) in the RMi system has no allowance for highly stressed, jointed rock mass. Accordingly, the RMi system recommends heavier support than what is already installed. Using SL ≈ 1/SRF, the ground quality factor decreases, which results in a heavier support.
Therefore, an alternate approach using a tunneling quality index support chart to extend application of the RMi system in highly stressed, jointed rock mass was used. The equations that were proposed in this study for SRF characterization were used for the rock quality index value determination, as applied to two tunnels that were in the design stage. The empirical support suggested by the rock quality index system support chart was numerically evaluated, and found that with application of that support, the total displacement and the extent of plastic zone decreased. The various thicknesses of shotcrete were also evaluated with a capacity diagram to confirm their stability.