Comparative Study on the Application of Different Slug Test Models for Determining the Permeability Coefﬁcients of Rock Mass in Long-Distance Deep Buried Tunnel Projects

: In large-scale water diversion projects, especially in the central and western regions of China, long-distance deep buried tunnels are generally involved. Therefore, it is essential to carry out ﬁeld tests to obtain the permeability of the rock mass through which tunnels pass. However, the test holes of large-scale water diversion projects are basically located in mountain areas with complex hydrogeological conditions. Meanwhile, the test holes are far apart and large in depth. As a result, traditional pumping tests cannot meet the requirements. Therefore, the slug test was chosen as the main test method, and the calculation results of the water injection test, the water pressure test and the slug test are analyzed and compared. The calculation results of the three test methods are basically consistent. However, the water injection test and the water pressure test are difﬁcult to implement at a large scale due to many environmental constraints, complex test equipment, long test periods and other factors. Furthermore, the Kipp model, the CBP model and the proposed HWS model, considering the effect of the ﬁnite thickness well-skin layer for the ﬁrst time, were used to analyze and process the slug test data, respectively. The curve ﬁtting effect of the Kipp model was the best, but the calculations were generally larger. The difference between the CBP model and the proposed HWS model is smaller in the calculation results; however, the curve ﬁtting effect of the CBP model is the worst, and the CBP model needs to be further improved. The curve ﬁtting effect of the proposed HWS model was between that of the Kipp model and the CBP model, and the proposed HWS model can be applied to the parameter calculations of the slug test with well-skin. In general, with reference to the criteria for the damping coefﬁcient of the aquifer in the Kipp model, the Kipp model was applicable to the slug test for test holes without well-skin and an aquifer damping coefﬁcient between 0.1 and 5.0. The CBP model was applicable to the slug test under the conditions of no well-skin and an aquifer damping coefﬁcient greater than 2.0. The novel proposed HWS model was applicable to the slug test when the aquifer damping coefﬁcient was greater than 1.0 under the conditions of no well-skin, positive well-skin and negative well-skin.


Introduction
The permeability coefficient of rock mass is a very important hydrogeological parameter in groundwater movement, especially for large-scale water diversion projects with deep buried tunnels. The permeability of the rock mass through which tunnels pass directly affects the construction safety and project cost. Meanwhile, the tunnel construction will disturb the surrounding groundwater environment and affect the drinking groundwater resources for residents in the surrounding mountain areas. Therefore, it is necessary to determine the permeability coefficient of rock mass in the study area during the preliminary investigation. Traditional test methods for determining the permeability coefficients include the pumping test, water pressure test, water injection test, etc. The traditional test

Fundamentals of the Kipp Model
Kipp [5] derived differential equations characterizing water level variations in the test hole based on the unsteady flow control equation in the aquifer, the mass conservation equation of water between the test hole and aquifer and the momentum balance equation of water motion in the test hole. The solution of the equation was obtained by Laplace transform and its inverse, and then the hydrogeological parameters were determined by fitting the standard curve.
The groundwater continuity equation in the aquifer is as follows: S T According to the flow rate balance, the change in the flow in the well caused by oscillation is equal to the radial flow of the aquifer around the well wall: The differential equation of the oscillation of water flow in the test hole based on the mass conservation and momentum balance is as follows: where L e is the effective water column length, L e = L + r 2 c r 2 w B 2 ; With initial conditions: when t = 0, h = h 0 , w = w 0 , dw dt = w 0 , h w = L = h 0 ; With boundary conditions: when r → ∞ , h → h 0 . These equations can be converted to a dimensionless form by the dimensionless variables and parameters and solved by the application of Laplace transform and inverse transform. Then, a set of w D ∼ τ standard curves can be obtained, as shown in Figure 1. transform and its inverse, and then the hydrogeological parameters were determined by fitting the standard curve.
The groundwater continuity equation in the aquifer is as follows: According to the flow rate balance, the change in the flow in the well caused by oscillation is equal to the radial flow of the aquifer around the well wall: The differential equation of the oscillation of water flow in the test hole based on the mass conservation and momentum balance is as follows: where is the effective water column length, = + ; With initial conditions: when = 0, ℎ = ℎ , = , = , ℎ = = ℎ ; With boundary conditions: when → ∞, ℎ → ℎ . These equations can be converted to a dimensionless form by the dimensionless variables and parameters and solved by the application of Laplace transform and inverse transform. Then, a set of ~ standard curves can be obtained, as shown in Figure 1. The aquifer permeability coefficient is calculated as follows: The − measured curve can be plotted according to the slug test data. Keep the vertical coordinate axis of the measured curve and the standard curve aligned at the same horizontal height. Then, the measured curve is shifted along the horizontal axis to make the measured curve fit the standard curve, and the corresponding and values of the fitted standard curve are recorded. At the same time, select any matching point and record The aquifer permeability coefficient is calculated as follows: The w − t measured curve can be plotted according to the slug test data. Keep the vertical coordinate axis of the measured curve and the standard curve aligned at the same horizontal height. Then, the measured curve is shifted along the horizontal axis to make the measured curve fit the standard curve, and the corresponding ζ and β values of the fitted standard curve are recorded. At the same time, select any matching point and record the corresponding coordinate values ([τ] and [t]) of the standard curve and the measured curve. Thus, the aquifer permeability coefficient can be obtained by Equation (4).
Among them, h is the aquifer water level variation (m); t is the time (s); τ is the dimensionless time; S is the storage coefficient; T is the hydraulic conductivity coefficient m 2 /s ; K is the permeability coefficient (m/s); r c is the radius of the well casing (m); r w is the radius of the well screen (m); w is the well water level variation (m); w 0 is the initial well water level variation (m); L is the static water column length above the top of the aquifer (m); B is the aquifer thickness (m); h w is the water head of the screen (m); ζ is the dimensionless damping parameter; β is the dimensionless inertial parameter, obtained by iterative calculation; w D is the dimensionless water level variation of the standard curve; w is the dimensionless water level variation of the measured curve.

Fundamentals of the CBP Model
The CBP model is a typical conductive slug test theory whose purpose is to obtain the permeability coefficient by fitting a standard curve with a semi-log plot of the water level variation after being dimensionless in the test hole with time [4]. The groundwater continuity equation in the aquifer is the same as Equation (1).
Cooper derived the water head variation relationship as follows: Among them, w D = w w 0 ; F(ϕ, τ) = 8ϕ A set of w D − τ standard curves can be obtained, as shown in Figure 2.
the corresponding coordinate values ([ ] and [ ]) of the standard curve and the measured curve. Thus, the aquifer permeability coefficient can be obtained by Equation (4). Among them, ℎ is the aquifer water level variation (m); is the time (s); is the dimensionless time; is the storage coefficient; is the hydraulic conductivity coefficient (m /s); is the permeability coefficient (m/s); is the radius of the well casing (m); is the radius of the well screen (m); is the well water level variation (m); is the initial well water level variation (m); is the static water column length above the top of the aquifer (m); is the aquifer thickness (m); ℎ is the water head of the screen (m); is the dimensionless damping parameter; is the dimensionless inertial parameter, obtained by iterative calculation; is the dimensionless water level variation of the standard curve; ′ is the dimensionless water level variation of the measured curve.

Fundamentals of the CBP Model
The CBP model is a typical conductive slug test theory whose purpose is to obtain the permeability coefficient by fitting a standard curve with a semi-log plot of the water level variation after being dimensionless in the test hole with time [4]. The groundwater continuity equation in the aquifer is the same as Equation (1).
Cooper derived the water head variation relationship as follows: A set of − standard curves can be obtained, as shown in Figure 2. The aquifer permeability coefficient is calculated as follows:  The aquifer permeability coefficient is calculated as follows: The w − t measured curve can be plotted according to the slug test data. Keep the vertical coordinate axis of the measured curve and the standard curve aligned at the same horizontal height. Then, the measured curve is shifted along the horizontal axis to fit the standard curve. At the same time, select any matching point and record the corresponding coordinate values ([τ] and [t]). Thus, the aquifer permeability coefficient can be obtained by Equation (6).
Among them, ϕ is the dimensionless storage coefficient; I 0 (u) and K 0 (u) are the firstand second-type modified Bessel functions of the zero order; I 1 (u) and K 1 (u) are the firstand second-type modified Bessel functions of the order 1; the rest of the symbols have the same meaning as they do above.

Fundamentals of the HWS Model
A novel HWS model (HHU Well-Skin model, abbreviated as the HWS model) is proposed by us for the first time. Meanwhile, a parameter calculation method of the slug test considering the influence of the finite thickness well-skin layer based on the HWS model is proposed.
The groundwater continuity equation is established between the confined aquifer and the finite thickness well-skin layer. The water balance equation is established in the contact surface between the confined aquifer and the finite thickness well-skin layer and between the finite thickness well-skin layer and the well. The theoretical model of the slug test under the influence of the well-skin is shown in Figure 3.
Among them, is the dimensionless storage coefficient; ( ) and ( ) are the first-and second-type modified Bessel functions of the zero order; ( ) and ( ) are the first-and second-type modified Bessel functions of the order 1; the rest of the symbols have the same meaning as they do above.

Fundamentals of the HWS Model
A novel HWS model (HHU Well-Skin model, abbreviated as the HWS model) is proposed by us for the first time. Meanwhile, a parameter calculation method of the slug test considering the influence of the finite thickness well-skin layer based on the HWS model is proposed.
The groundwater continuity equation is established between the confined aquifer and the finite thickness well-skin layer. The water balance equation is established in the contact surface between the confined aquifer and the finite thickness well-skin layer and between the finite thickness well-skin layer and the well. The theoretical model of the slug test under the influence of the well-skin is shown in Figure 3. (1) Establishment of the theoretical model The governing equation of the groundwater radial flow in a well-skin layer and a confined aquifer is as follows: A water balance equation is established in the contact surface between the confined aquifer and the well-skin layer of finite thickness as follows: A water balance equation is established in the contact surface between the well-skin layer of finite thickness and the well as follows: (1) Establishment of the theoretical model The governing equation of the groundwater radial flow in a well-skin layer and a confined aquifer is as follows: A water balance equation is established in the contact surface between the confined aquifer and the well-skin layer of finite thickness as follows: A water balance equation is established in the contact surface between the well-skin layer of finite thickness and the well as follows: Among them, subscripts 1 and 2 represent the well-skin layer and the confined aquifer parameters, respectively, r s is the distance from the center of the well to the outer boundary of the well-skin layer (m) and the rest of the symbols have the same meaning as they do above. (2) Solution of the theoretical model First, these equations can be converted to a dimensionless form by the following dimensionless variables and parameters: Then, Laplace transform is used to solve the set of Equations (7)-(10) in turn, as follows: According to the general solution of the modified Bessel equation of the order 0, the solutions of Equations (12) and (13) are obtained as follows: Among them, C 1 and C 2 , C 3 and C 4 are constants determined by boundary conditions in the well-skin layer and the aquifer, respectively.
By substituting C 1 , C 2 , C 3 and C 4 into Equations (16) and (17), respectively, the solutions of Equations (12) and (13) can be obtained as follows: The final solution in the Laplace transform domain is obtained as follows: The analytical solution of the model in the Laplace transform domain is obtained by using the AWG (Azari-Wooden-Gaver) algorithm [18,19] in the inverse Laplace transform in Equation (20), as follows: Among them, w D (τ) is the dimensionless water level variation; w D (p) is the dimensionless water level variation in the Laplace transform domain; N is the number of summation terms; the rest of the symbols have the same meaning as they do above.
According to Equations (20) and (21), the standard curve data corresponding to different γ values when the well-skin factor is α = 1, β = 1, ρ s = 1 are calculated. Then, a set of w D ∼ τ standard curves without the influence of the well-skin effect can be obtained by plotting the curve data, as shown in Figure 4. In the same way, the standard curve data corresponding to different α values when β = 50, γ = 1000 and ρ s = 3 are calculated. Then, a set of w D ∼ τ standard curves under the influence of the well-skin effect can be obtained by plotting the curve data, as shown in Figure 5. Among them, ( ) is the dimensionless water level variation; ( ) is the dimensionless water level variation in the Laplace transform domain; is the number of summation terms; the rest of the symbols have the same meaning as they do above.
According to Equations (20) and (21), the standard curve data corresponding to different values when the well-skin factor is = 1, = 1, = 1 are calculated. Then, a set of ~ standard curves without the influence of the well-skin effect can be obtained by plotting the curve data, as shown in Figure 4. In the same way, the standard curve data corresponding to different values when = 50, = 1000 and = 3 are calculated. Then, a set of ~ standard curves under the influence of the well-skin effect can be obtained by plotting the curve data, as shown in Figure 5.  The aquifer permeability coefficient is calculated as follows: Among them, ( ) is the dimensionless water level variation; ( ) is the dimensionless water level variation in the Laplace transform domain; is the number of summation terms; the rest of the symbols have the same meaning as they do above.
According to Equations (20) and (21), the standard curve data corresponding to different values when the well-skin factor is = 1, = 1, = 1 are calculated. Then, a set of ~ standard curves without the influence of the well-skin effect can be obtained by plotting the curve data, as shown in Figure 4. In the same way, the standard curve data corresponding to different values when = 50, = 1000 and = 3 are calculated. Then, a set of ~ standard curves under the influence of the well-skin effect can be obtained by plotting the curve data, as shown in Figure 5.  The aquifer permeability coefficient is calculated as follows: The aquifer permeability coefficient is calculated as follows: Similarly, the w − t measured curve can be plotted according to the slug test data. Keep the vertical coordinate axis of the measured curve and the standard curve aligned at the same horizontal height. Then, the measured curve is shifted along the horizontal axis to fit the standard curve, and the corresponding α, β and γ values of the fitted standard curve are recorded. At the same time, select any matching point, and record the corresponding coordinate values ([τ] and [t]). Thus, the aquifer permeability coefficient can be obtained by Equation (23).

Fundamentals of the Water Injection Test
According to the Chinese "Code of water injection test for water resources and hydropower engineering" (SL 345-2007) [20], when carrying out the variable water head water injection test in the test hole, it is necessary to first seal the water at the top and bottom of the test section. Then, inject a certain amount of water into the casing in order to make the water level in the casing higher than the groundwater level by a certain height, and, finally, record the change in water level in the casing with time. Draw the relation curve between the water head ratio and the time on semi-logarithmic coordinate paper. When the relationship between the water head ratio and the time is not linear, it should be checked and retested. The test can be ended when the water head of the test hole drops to 0.3 times the initial water head. The variable water head water injection test is applicable to the silt, cohesive soil or rock mass with a low permeability coefficient below the groundwater level.
The aquifer permeability coefficient is calculated as follows: Among them, T 0 is the characteristic time of the water injection test (s), the value of t corresponding to H t /H 0 = 0.37. It can be determined on the ln(H t /H 0 ) − t curve. H t is the water head value when the time is t (m). H 0 is the initial water head value of the water injection test (m). r c is the inner radius of the well casing (m). A is the shape factor (m).

Fundamentals of the Conventional Water Pressure Test
According to the Chinese "Specification for the Water Pressure Test in the Borehole of Hydropower Projects" (NB/T 35113-2018) [21], generally, a conventional water pressure test [22,23] is always carried out with five phases (P 1 − P 2 − P 3 − P 4 (= P 2 ) − P 5 (= P 1 ), P 1 < P 2 < P 3 ) and three pressure stages, P 1 , P 2 and P 3 , which are 0.3 MPa, 0.6 MPa and 1 MPa, respectively. The water permeability of the test section is calculated by the pressure value (P 3 ) and flow rate (Q 3 ) of the third pressure stage in Equation (25) [24].
Among them, q is the water permeability of the test section (Lu or L/min/m); L is the length of the test section (m); Q 3 is the flow rate of the third pressure stage (L/min); P 3 is the test section pressure of the third pressure stage (MPa).

Central Yunnan Water Diversion Project Profile
The central Yunnan water diversion project is the largest and most invested water resources allocation project in Southwest China. The total length of the main water diversion channel is about 661 km, and there are 63 main tunnels with a total length of about 610 km. The underground tunnels account for more than 90% of the total length of the total route length [25]. The water diversion route comprises the following sections: Dali I, Dali II, Chuxiong, Kunming, Yuxi and Honghe. The location and water delivery route of the central Yunnan water diversion project are shown in Figure 6. It passes through some active faults of Holocene and Late-Pleistocene, and the impact of these active faults on the diversion route during the construction-especially the groundwater environment-cannot be ignored [26]. Among them, the Xianglushan tunnel is the control project of the central Yunnan water diversion project, located at the head of Dali I. It is the longest deep-buried tunnel of the project, with a total length of about 63 km and a maximum buried depth of 1512 m. Meanwhile, it is a large-slope inclined shaft tunnel with the most difficult construction and complex hydrogeological conditions along the whole line. The tunnel construction has problems such as water inrush and mud inrush, a high external water pressure, crossing active faults, a large buried depth, a difficult construction and groundwater treatment [27,28]. After the completion of the project, water can be diverted from the main stream of Jinsha River, with relatively abundant water, to central Yunnan, which will solve severe water shortages in the central Yunnan region, improve the ecological and environmental conditions of rivers and lakes in the region and effectively promote the sustainable economic and social development of Yunnan [29,30].
Dali II, Chuxiong, Kunming, Yuxi and Honghe. The location and water delivery route of the central Yunnan water diversion project are shown in Figure 6. It passes through some active faults of Holocene and Late-Pleistocene, and the impact of these active faults on the diversion route during the construction-especially the groundwater environment-cannot be ignored [26]. Among them, the Xianglushan tunnel is the control project of the central Yunnan water diversion project, located at the head of Dali I. It is the longest deepburied tunnel of the project, with a total length of about 63 km and a maximum buried depth of 1512 m. Meanwhile, it is a large-slope inclined shaft tunnel with the most difficult construction and complex hydrogeological conditions along the whole line. The tunnel construction has problems such as water inrush and mud inrush, a high external water pressure, crossing active faults, a large buried depth, a difficult construction and groundwater treatment [27,28]. After the completion of the project, water can be diverted from the main stream of Jinsha River, with relatively abundant water, to central Yunnan, which will solve severe water shortages in the central Yunnan region, improve the ecological and environmental conditions of rivers and lakes in the region and effectively promote the sustainable economic and social development of Yunnan [29,30].

Slug Test of the Central Yunnan Water Diversion Project
As the central Yunnan water diversion project is located in a mountainous area with complex geological conditions, there are many difficulties in conducting hydrogeological tests in this area. For example, it is difficult to adopt the traditional pumping test for investigation, and the water pressure test is also limited by the site and equipment. A lot of field test work has been carried out in order to seek the theory and method applicable to the field test of the rock mass hydrogeological parameters in deep test hole. Among them, the borehole CYWD-1 is specially selected to conduct sectional comparative tests. With top-pressure plug water sealing, seven groups of sectional slug tests and variable water head water injection tests are carried out every 5-25 m according to the drilling progress and are numbered CYWD-1-1, CYWD-1-2, CYWD-1-3, CYWD-1-4, CYWD-1 CYWD-1-5, CYWD-1-6 and CYWD-1-7. Among them, the first test section is located in the borehole depth of 130-135 m; the second test section is located in the borehole depth of 135-143 m; the third test section is located in the borehole depth of 143-161 m; the fourth test section is located in the borehole depth of 161-173 m; the fifth test section is located in the borehole depth of 173-182 m; the sixth test section is located in the borehole depth of 182-196 m; and the seventh test section is located in the borehole depth of 182-196 m. The test data are respectively calculated and analyzed with the variable water head injection test and three theoretical methods of the slug test. The slug test data are fitted with the standard curves of three slug test models according to the Kipp model, CBP model and HWS model. The fitted standard curves are shown in Figure 7.

Slug Test of the Central Yunnan Water Diversion Project
As the central Yunnan water diversion project is located in a mountainous area with complex geological conditions, there are many difficulties in conducting hydrogeological tests in this area. For example, it is difficult to adopt the traditional pumping test for investigation, and the water pressure test is also limited by the site and equipment. A lot of field test work has been carried out in order to seek the theory and method applicable to the field test of the rock mass hydrogeological parameters in deep test hole. Among them, the borehole CYWD-1 is specially selected to conduct sectional comparative tests. With top-pressure plug water sealing, seven groups of sectional slug tests and variable water head water injection tests are carried out every 5-25 m according to the drilling progress and are numbered CYWD-1-1, CYWD-1-2, CYWD-1-3, CYWD-1-4, CYWD-1 CYWD-1-5, CYWD-1-6 and CYWD-1-7. Among them, the first test section is located in the borehole depth of 130-135 m; the second test section is located in the borehole depth of 135-143 m; the third test section is located in the borehole depth of 143-161 m; the fourth test section is located in the borehole depth of 161-173 m; the fifth test section is located in the borehole depth of 173-182 m; the sixth test section is located in the borehole depth of 182-196 m; and the seventh test section is located in the borehole depth of 182-196 m. The test data are respectively calculated and analyzed with the variable water head injection test and three theoretical methods of the slug test. The slug test data are fitted with the standard curves of three slug test models according to the Kipp model, CBP model and HWS model. The fitted standard curves are shown in Figure 7.  As shown in Figure 7, the standard curve shape of the three models is basically consistent. When the water level recovery rate gets faster, the CBP model will at first be poorly fitted. As the rate is further accelerated, especially when the aquifer damping parameter is less than 1, the HWS model will also appear to be unable to completely fit in the part of the beginning and end. Since the Kipp model has a large variation in the shape of the standard curve, it can match the measured curves at different water level recovery rates to the greatest extent, so the Kipp model has the best fit in the matching process, and the HWS model has the second best fit. The overall shape of the standard curve cluster of the CBP model is relatively gradual, and the fitting effect of the CBP model is relatively poor when the water level recovers rapidly. In general, the fitting effect of all three models can meet the requirements of the calculation.
In addition, an on-site investigation of the boreholes along the Shigu water source area and the Xianglushan tunnel of the central Yunnan water diversion project is adopted. The whole hole section slug test studies are carried out in the boreholes of CYWD-2, CYWD-3, CYWD-4 and CYWD-5. Three slug test models are used for the fitting curve and parameter calculation. The fitted standard curves are shown in Figure 8. As shown in Figure 7, the standard curve shape of the three models is basically consistent. When the water level recovery rate gets faster, the CBP model will at first be poorly fitted. As the rate is further accelerated, especially when the aquifer damping parameter is less than 1, the HWS model will also appear to be unable to completely fit in the part of the beginning and end. Since the Kipp model has a large variation in the shape of the standard curve, it can match the measured curves at different water level recovery rates to the greatest extent, so the Kipp model has the best fit in the matching process, and the HWS model has the second best fit. The overall shape of the standard curve cluster of the CBP model is relatively gradual, and the fitting effect of the CBP model is relatively poor when the water level recovers rapidly. In general, the fitting effect of all three models can meet the requirements of the calculation.
In addition, an on-site investigation of the boreholes along the Shigu water source area and the Xianglushan tunnel of the central Yunnan water diversion project is adopted. The whole hole section slug test studies are carried out in the boreholes of CYWD-2, CYWD-3, CYWD-4 and CYWD-5. Three slug test models are used for the fitting curve and parameter calculation. The fitted standard curves are shown in Figure 8. As shown in Figure 8, the standard curve shapes fitted by the three models are basically consistent. Among them, the data of the water level recovery of the boreholes CYWD-2, CYWD-3 and CYWD-4 are relatively stable, and the curve shape is comparatively smooth during the slug test. Meanwhile, the water level recovery time of these three boreholes is relatively long. The qualitative analysis shows that the rock mass permeability around these three boreholes is relatively low. The data of the water level recovery of As shown in Figure 8, the standard curve shapes fitted by the three models are basically consistent. Among them, the data of the water level recovery of the boreholes CYWD-2, CYWD-3 and CYWD-4 are relatively stable, and the curve shape is comparatively smooth during the slug test. Meanwhile, the water level recovery time of these three boreholes is relatively long. The qualitative analysis shows that the rock mass permeability around these three boreholes is relatively low. The data of the water level recovery of borehole CYWD-5 have some fluctuations due to some certain disturbances during the test. The speed of the water level recovery of borehole CYWD-5 is fast, and the rock mass permeability around this borehole should be high. According to the analysis of the matching curve, the water level recovery curve of the three boreholes CYWD-2, CYWD-3 and CYWD-4 is basically the same, and all of them fit the standard curve with a damping coefficient of 1 in the Kipp model. The CBP model has a poor fitting to these three groups of data, as the overall shape of the standard curve is relatively gradual, and there is a comparatively large deviation at the beginning and the end. The HWS model can basically fit the three groups of data well. Although the recovery time of the water level of borehole CYWD-5 is relatively short during the slug test, the water level is disturbed to a certain extent, leading to a certain fluctuation effect in the water level recovery data. In general, the Kipp model fits best, the HWS model fits the second best and the CBP model fits relatively poorly.

Water Diversion Project from the Three Gorges Reservoir to Hanjiang River Profile
The Water Diversion Project from the Three Gorges Reservoir to Hanjiang River was started in Danjiangkou City, Hubei Province on 7 July 2022. The project of the water diversion project from the Three Gorges Reservoir to Hanjiang River is the first project in the follow-up project of the South to North Water Diversion Project. The project takes water from Longtan Creek on the left bank of the Three Gorges Reservoir of the Yangtze River, passes through Yichang, Xiangyang and Danjiangkou of Hubei Province and ends at the Anle River mouth on the right bank of Hanjiang River, which is downstream of Danjiangkou Reservoir. The total length is 194.8 km, and the design diversion flow is 170 m 3 /s. The location and water delivery route of the project is shown in Figure 9. The whole route of the project adopts tunnels for water conveyance, which is the longest pressurized water conveyance tunnel under construction in China. Approximately 2.2 × 10 10 m 3 of tunnel muck would be taken out from the underground [31].

Slug Test of the Water Diversion Project from the Three Gorges Reservoir to Hanjiang River
According to the route comparison scheme of the river diversion project and the boreholes situation, the permeability of six boreholes located in the main fault fracture zone was studied, and two sections of fractured strata were selected for testing in each borehole. In total, 12 groups of slug tests were carried out, namely, WDTGH-1-1, WDTGH-1-2, WDTGH-2-1, WDTGH-2-2, WDTGH-3-1, WDTGH-3-2, WDTGH-4-1, WDTGH-4-2, WDTGH-5-1, WDTGH-5-2, WDTGH-6-1 and WDTGH-6-2. Among them, After the implementation of the project, it will connect the South-to-North Water Diversion Project and the Three Gorges project, as well as the Three Gorges Reservoir and the Danjiangkou Reservoir, and, moreover, link the Yangtze River, the Han River Basin and the Beijing-Tianjin-Hebei-Henan area and improve the national important water network layout. In the meantime, it will also supplement water to the middle and downstream areas of the Hanjiang River, which play an important role in improving the water resources allocation capacity of the Hanjiang River Basin and the water ecological environment of the middle and downstream area of the Hanjiang River [32,33]. The project adopts the pressure single tunnel artesian water conveyance, which is the most difficult long-distance water diversion tunnel project under construction in China. The maximum buried depth is 1182 m, and the tunnel with a buried depth of over 600 m accounts for 45%. The route is long, and the buried depth is large. With high mountains, deep valleys, developed faults and folds, the complex terrain and geological conditions, the project faces multiple challenges, such as water inrush, mud inrush, large faults, etc.

Slug Test of the Water Diversion Project from the Three Gorges Reservoir to Hanjiang River
According to the route comparison scheme of the river diversion project and the boreholes situation, the permeability of six boreholes located in the main fault fracture zone was studied, and two sections of fractured strata were selected for testing in each borehole. In total, 12 groups of slug tests were carried out, namely, WDTGH-1-1, WDTGH-1-2, WDTGH-2-1, WDTGH-2-2, WDTGH-3-1, WDTGH-3-2, WDTGH-4-1, WDTGH-4-2, WDTGH-5-1, WDTGH-5-2, WDTGH-6-1 and WDTGH-6-2. Among them, two groups of water pressure tests were carried out synchronously in borehole WDTGH-1 for comparative study. Similarly, according to the theories of the Kipp model, CBP model and HWS model, the slug test data were fitted with the standard curves of the three slug test models. The fitted standard curves are shown in Figure 10.

Slug Test of the Water Diversion Project from the Three Gorges Reservoir to Hanjiang River
According to the route comparison scheme of the river diversion project and the boreholes situation, the permeability of six boreholes located in the main fault fracture zone was studied, and two sections of fractured strata were selected for testing in each borehole. In total, 12 groups of slug tests were carried out, namely, WDTGH-1-1, WDTGH-1-2, WDTGH-2-1, WDTGH-2-2, WDTGH-3-1, WDTGH-3-2, WDTGH-4-1, WDTGH-4-2, WDTGH-5-1, WDTGH-5-2, WDTGH-6-1 and WDTGH-6-2. Among them, two groups of water pressure tests were carried out synchronously in borehole WDTGH-1 for comparative study. Similarly, according to the theories of the Kipp model, CBP model and HWS model, the slug test data were fitted with the standard curves of the three slug test models. The fitted standard curves are shown in Figure 10.   As shown in Figure 10, the water level recovery rate of these 12 groups of tests is basically gradual, which all fit best with the Kipp model standard curve of the damping coefficient greater than 1. For the CBP model, the smaller the water level recovery rate is, the larger the damping coefficient is, leading to a better fitting effect of the CBP model. When the damping coefficient is less than 3, there will still be poor fitting in some areas, and the damping coefficient is greater than or equal to 3. The fitting effect of the CBP model can basically reach the best. The fitting effect of the HWS model is better than that of the CBP model in general. When the damping coefficient is greater than 1, the HWS model can basically achieve a better fitting effect. The shape difference between different standard curves of the Kipp model is relatively large. The standard curve of different shapes can be selected for fitting according to the actual measured data, so the Kipp model can have more choices in the matching curve, and the fitting effect of the matching curve is also the best. The shape of the standard curve of the CBP model has little difference, and the overall gradient is small. Thus, when conducting the slug test in the high permeability stratum, the measured data can only be matched with the standard curve of the dimensionless water storage coefficient ϕ = 10 −5 with the largest gradient. When the water level recovers quickly, the fitting effect using the CBP model is not ideal. The standard curve shape of the HWS model is between the other two models. When using the HWS model for fitting, there is a large choice, and its fitting effect is also between the other two models. In addition, according to the qualitative analysis of the water level recovery time in Figure 10, it can be qualitatively understood that the rock mass permeability in the tested sections of WDTGH-1-1, WDTGH-4-2 and WDTGH-6-1 with a slow water level recovery is low.

Comparative Analysis of the Sectional Variable Water Head Water Injection Test and Slug Test of Borehole CYWD-1
A sectional comparison test was carried out in the borehole CYWD-1. The sectional variable water head water injection test and three different slug test theories were used for analysis and calculation. The statistics of the calculation results are shown in Table 1, and a comprehensive comparison of the calculation results is shown in Figure 11.  Figure 11. Comparison of the calculation results water injection test and different slug test theory models in borehole CYWD-1.
It is shown in Table 1 and Figure 11 that the calculated results using the Kipp model are generally large, those using the CBP model are generally small and the calculated results using the water injection test method fluctuate relatively largely. Through a comprehensive comparison, the results calculated using the HWS model are most consistent with the average value, and the results obtained are more representative. This indicates that the slug test is more stable and reliable compared to the water injection test. In general, the water injection test results of each test section are basically consistent with the corresponding slug test results, and the ratio of the maximum and minimum permeability coefficients calculated by different methods is within 2.5 times. Therefore, the calculated results of the sectional water injection test and slug test carried out in the field are reliable, and it also shows that the HWS model derived by considering the finite thickness wellskin effect is also very practical and can be popularized when the well-skin factor is 1.

Comparative Analysis of the Whole Hole Slug Test of the Central Yunnan Water Diversion Project
The whole hole slug test is carried out in the boreholes CYWD-2, CYWD-3, CYWD-4 and CYWD-5. Three different slug test models are used for analysis and calculation. The statistics of the calculation results are shown in Table 2, and a comprehensive comparison of the calculation results is shown in Figure 12. It is shown in Table 1 and Figure 11 that the calculated results using the Kipp model are generally large, those using the CBP model are generally small and the calculated results using the water injection test method fluctuate relatively largely. Through a comprehensive comparison, the results calculated using the HWS model are most consistent with the average value, and the results obtained are more representative. This indicates that the slug test is more stable and reliable compared to the water injection test. In general, the water injection test results of each test section are basically consistent with the corresponding slug test results, and the ratio of the maximum and minimum permeability coefficients calculated by different methods is within 2.5 times. Therefore, the calculated results of the sectional water injection test and slug test carried out in the field are reliable, and it also shows that the HWS model derived by considering the finite thickness well-skin effect is also very practical and can be popularized when the well-skin factor is 1.

Comparative Analysis of the Whole Hole Slug Test of the Central Yunnan Water Diversion Project
The whole hole slug test is carried out in the boreholes CYWD-2, CYWD-3, CYWD-4 and CYWD-5. Three different slug test models are used for analysis and calculation. The statistics of the calculation results are shown in Table 2, and a comprehensive comparison of the calculation results is shown in Figure 12.  According to Table 2 and Figure 12, during the full borehole test, the calculated results of the Kipp model are still the largest among the three models. The results calculated by the CBP model and HWS model are close to and consistent with the average value, and those of the HWS model are relatively smaller. Nevertheless, the calculated results of the three models are basically consistent, and the ratio of the maximum and minimum permeability coefficients calculated by different methods is about 1.5 times. Meanwhile, combined with the qualitative analysis results in Figure 8, the quantitative calculation results in Table 2 also show that the permeability of the rock mass in the boreholes CYWD-2, CYWD-3 and CYWD-4 is low, and the permeability of the rock mass in borehole CYWD-5 is high. The qualitative and quantitative analysis results are completely consistent, which shows that the slug test has a good field operation stability, and the results calculated by the three models are reliable. The sectional comparison test is carried out in borehole WDTGH-1. The analysis and calculations were carried out by using the water pressure test and three different slug test According to Table 2 and Figure 12, during the full borehole test, the calculated results of the Kipp model are still the largest among the three models. The results calculated by the CBP model and HWS model are close to and consistent with the average value, and those of the HWS model are relatively smaller. Nevertheless, the calculated results of the three models are basically consistent, and the ratio of the maximum and minimum permeability coefficients calculated by different methods is about 1.5 times. Meanwhile, combined with the qualitative analysis results in Figure 8, the quantitative calculation results in Table 2 also show that the permeability of the rock mass in the boreholes CYWD-2, CYWD-3 and CYWD-4 is low, and the permeability of the rock mass in borehole CYWD-5 is high. The qualitative and quantitative analysis results are completely consistent, which shows that the slug test has a good field operation stability, and the results calculated by the three models are reliable. The sectional comparison test is carried out in borehole WDTGH-1. The analysis and calculations were carried out by using the water pressure test and three different slug test models. The statistics of the calculation results are shown in Table 3. According to the relationship between water permeability and the permeability coefficient in the Chinese "Code for engineering geological investigation of water resources and hydropower" (GB50487-2008) [34], the comprehensive comparison of the calculation results is shown in Figure 13. relationship between water permeability and the permeability coefficient in the Chinese "Code for engineering geological investigation of water resources and hydropower" (GB50487-2008) [34], the comprehensive comparison of the calculation results is shown in Figure 13. It can be seen from Table 3 and Figure 13 that, during the water pressure test and slug test in borehole WDTGH-1, the Kipp model still has a large calculated result among the three slug test models, and the HWS model has a small one at the 91-101 m test section. The calculated results of the CBP model and the HWS model are close at the 217-227 m test section, and the CBP model is relatively small. The calculated results of the three models are basically consistent. The ratio of the maximum and the minimum of the permeability coefficient calculated by the different methods is about two times, and the statistical relationship with the results of the water pressure test is also within the confidence intervals. This shows that both the water pressure test and slug test can obtain the rock mass permeability in the fault fracture zone, but the equipment and test requirements for the water pressure test are higher, and high-pressure water easily causes borehole collapse in the test section of the fractured rock, so only the slug test is carried out in the other boreholes. It can be seen from Table 3 and Figure 13 that, during the water pressure test and slug test in borehole WDTGH-1, the Kipp model still has a large calculated result among the three slug test models, and the HWS model has a small one at the 91-101 m test section. The calculated results of the CBP model and the HWS model are close at the 217-227 m test section, and the CBP model is relatively small. The calculated results of the three models are basically consistent. The ratio of the maximum and the minimum of the permeability coefficient calculated by the different methods is about two times, and the statistical relationship with the results of the water pressure test is also within the confidence intervals. This shows that both the water pressure test and slug test can obtain the rock mass permeability in the fault fracture zone, but the equipment and test requirements for the water pressure test are higher, and high-pressure water easily causes borehole collapse in the test section of the fractured rock, so only the slug test is carried out in the other boreholes.

Comparison and Analysis of the Sectional Slug Test of the Water Diversion Project from the Three Gorges Reservoir to Hanjiang River
Sectional slug tests were carried out in the boreholes WDTGH-2, WDTGH-3, WDTGH-4, WDTGH-5 and WDTGH-6. Two sections of broken zone strata were selected from each borehole for testing. Three different slug tests models were used for the analysis and calculation in each test. The statistics of the calculation results are shown in Table 4, and a comprehensive comparison of the calculation results is shown in Figure 14. Sectional slug tests were carried out in the boreholes WDTGH-2, WDTGH-3, WDTGH-4, WDTGH-5 and WDTGH-6. Two sections of broken zone strata were selected from each borehole for testing. Three different slug tests models were used for the analysis and calculation in each test. The statistics of the calculation results are shown in Table 4, and a comprehensive comparison of the calculation results is shown in Figure 14.  It can be seen in Table 4 and Figure 14 that when the sectional slug test is carried out in five boreholes, the calculation results of the Kipp model are still larger, those of the CBP model are in the middle and those of the HWS model are the smallest. Therefore, the calculation results of the CBP model are closest to the average value. However, when the CBP model is used for parameter calculation, due to the similarity of the standard curve shape, it is of great uncertainty to determine the relevant parameters. The principle of the CBP model can be further improved. Hence, this fitting method based on the CBP model It can be seen in Table 4 and Figure 14 that when the sectional slug test is carried out in five boreholes, the calculation results of the Kipp model are still larger, those of the CBP model are in the middle and those of the HWS model are the smallest. Therefore, the calculation results of the CBP model are closest to the average value. However, when the CBP model is used for parameter calculation, due to the similarity of the standard curve shape, it is of great uncertainty to determine the relevant parameters. The principle of the CBP model can be further improved. Hence, this fitting method based on the CBP model may lead to the inaccurate estimation of the aquifer permeability coefficient. The calculated results based on the three slug test models are relatively close, and the ratio between the maximum and the minimum of the permeability coefficient calculated by different slug test models is about 2.5 times. Similarly, combined with the qualitative analysis results in Figure 10, the quantitative calculation results in Tables 3 and 4 also show that the rock mass permeability in the test sections WDTGH-1-1, WDTGH-4-2 and WDTGH-6-1 is low, and the qualitative and quantitative analysis results are completely consistent, which further indicates the stability, reliability and accuracy of the slug test in determining the rock mass permeability coefficient.

1.
The water injection test, water pressure test and slug test can be applied to obtain the rock mass permeability coefficients. The permeability coefficients of the rock mass calculated by the three test methods are basically consistent. However, the equipment and operation of the water pressure test are complicated. It is easy to cause the borehole to collapse when encountering broken strata, such as fault zones. The applicability of this water pressure test is limited. In addition, the water pressure test and water injection test generally need a large amount of water injection when carrying out tests; therefore, it is not convenient to put into practice in some test sites where transportation is inconvenient and the water source is distant. Through comprehensive comparison, when determining the rock mass permeability coefficient in large-scale water diversion projects and long-distance deep buried tunnels, the slug test can meet the test requirements to the greatest extent, with a simple operation and strong applicability.

2.
Three theoretical models of the slug test are used to synchronously analyze and compare the data of the slug test. Among them, the calculated results of the Kipp model are generally large, and those of the CBP model and the HWS model are relatively close. However, the principle of the CBP model can be further improved.
The difference of the calculated results of the three models is about two times, and all of them can satisfy the requirements of determining the permeability of the rock mass. Meanwhile, it is verified that the Kipp model, the CBP model and the novel proposed HWS model have good applicability in the parameter calculation of the deep borehole slug test.

3.
The standard curves of the Kipp model are based on the standard solution form of the oscillatory equation, and the standard curves cover the situation from overdamping to underdamping. Therefore, the Kipp model has more choices in fitting curves and the best effect of fitting curves among the three models. The CBP model has a single standard curve cluster shape and a relatively gradual curve. When the slug test is conducted in high-permeability fractured rock mass, the measured data can only basically match the standard curves with the largest gradient, where the dimensionless water storage coefficient ϕ is equal to 10 − 5. When the aquifer damping coefficient is less than 1, the fitting effect of the CBP model will be relatively poor among the three models. The fitting effect of the HWS model is between the other two models. Through the research, it can be found that the calculated results of the Kipp model are larger, and the calculated results of the HWS model are smaller, so we should pay attention to the application scope and characteristics of different models. Therefore, when selecting different slug test models for determining the permeability coefficients, the permeability coefficients should be reasonably selected in combination with the drilling core data to avoid the errors of different slug test models.

4.
The novel proposed HWS model is a theoretical model of a slug test considering the influence of the finite thickness well-skin layer. In contrast, the other two models do not consider the influence of well-skin. Test holes are selected in the bedrock area with good integrity, and the disturbance around the test holes is small. Therefore, the standard curves with a well-skin factor of 1 are selected in the calculation. In the process of drilling, if the test hole is excessively washed or the wall protection material is used, the Kipp model and the CBP model will produce relatively large errors in calculation, so the HWS model is very practical and can be better popularized.