Next Article in Journal
Hydrological Drought in Dongting Lake Area (China) after the Running of Three Gorges Dam and a Possible Solution
Next Article in Special Issue
Evaluation of Ground Temperature Changes by the Operation of the Geothermal Heat Pump System and Climate Change in Korea
Previous Article in Journal
The Soil Water Evaporation Process from Mountains Based on the Stable Isotope Composition in a Headwater Basin and Northwest China
Previous Article in Special Issue
Changes in Geochemical Composition of Groundwater Due to CO2 Leakage in Various Geological Media
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distortion of the Estimated Hydraulic Conductivity from a Hydraulic Test in Fractured Rock Due to Excessive Injection or Extraction

Radioactive Waste Disposal Research Division, Korea Atomic Energy Research Institute, Daejeon 34057, Korea
*
Author to whom correspondence should be addressed.
Water 2020, 12(10), 2712; https://doi.org/10.3390/w12102712
Submission received: 26 August 2020 / Revised: 19 September 2020 / Accepted: 24 September 2020 / Published: 28 September 2020

Abstract

:
In this study, we discussed distortion of the estimated hydraulic conductivity from a hydraulic test due to excessive injection or extraction of groundwater by evaluating the influence of nonlinear flow. Pulse, slug, and constant head withdrawal tests with various head displacements were conducted in fractured granite rock, and the changes of representative Reynolds numbers (Re) during the tests were calculated. The Forchheimer equation and cubic law were used to evaluate the influence of nonlinear flow on the hydraulic tests, and thus the possibility of distortion of the estimated hydraulic conductivity. Our results showed that there was little possibility that nonlinear flow occurred during the pulse tests in the test zones. In the slug tests at several test zones, however, the estimated hydraulic conductivities were likely to be distorted due to nonlinear flow. Except for the test zones with low permeability, the scale effects of the estimated hydraulic conductivities from different types of tests were observed. These results indicated that the scale effect and distortion of the hydraulic parameters can be evaluated by conducting various types of hydraulic tests.

1. Introduction

The safety assessment of a subsurface repository for high-level radioactive waste (HLW) is to estimate the dose rates of the nuclides reached to the biosphere from the repository. One of the major pathways of the nuclides from the repository to the biosphere is groundwater flow, and the reliability of the safety assessment is controlled by the uncertainty in hydrogeological characterization of the disposal site. Crystalline fractured rock is one of the preferred host rocks for HLW disposal repositories. In a fractured rock, groundwater flows predominantly through fractures rather than matrix, which makes crystalline rock very heterogeneous and discrete hydraulically. Hydraulic tests for hydrogeological characterization are generally analyzed with the semi-analytical solutions induced from the assumption of a homogeneous and continuous medium, thus the hydrogeological characterization of a fractured rock is more uncertain than a porous medium.
The semi-analytical solutions for the hydraulic test analyses also assume the linear relation between the pressure gradient and the flux. However, it was reported that fluid flow in a fracture deviates from the linear relation at sufficiently large Reynolds numbers (Re), e.g., [1,2,3], and it was suggested that the nonlinear flow occurs due to the geometrical characteristics of a fracture and significant inertia of flowing fluid, e.g., [4,5,6]. In the hydraulic tests, such as pumping and slug tests in fractured rocks, nonlinear flow was also observed. Using the results of constant rate injection tests with various rates in a fractured sedimentary rock, Quinn et al. [7] analyzed the relation between the imposed injection rates and the resulting stabilized head changes and found that that relation deviated from linearity as the injection rate increased. From this result, they identified the nonlinear flow regime at the test zone. Quinn et al. [8] determined the critical flow rate, where the groundwater flow began to depart from the linear flow regime, from the constant-head step tests in a fractured sedimentary rock. By comparing the determined critical flow rate to the measured flow rates during slug tests, they suggested the appropriate initial displacement, below which the linear flow regime was guaranteed during a slug test at the test zone. Ji and Koh [9] conducted slug tests in fractured granite rock with various initial displacements. They found that a nonlinear flow arose during the slug test and influenced the test results when the initial displacement was large. They identified the nonlinear flow regimes during the slug tests by calculating a representative Re with the slug test and borehole logging results. Then, they suggested the way to estimate the hydraulic conductivity in a linear flow regime by using the Forchheimer equation and data from the slug tests in a nonlinear flow regime. Chen et al. [10] introduced the relationship curves between flow rate and injection pressure from the constant head injection tests in fractured granite rock, and divided the flow regimes during the tests into a laminar, non-Darcian, and hydraulic fracturing phases based on the relation between the flow rate and pressure. Then, they estimated the hydraulic conductivities of the test zones in different flow regimes using an Izbash’s law based empirical model.
As reviewed in the previous researches, excessive injection or extraction of water may induce a nonlinear flow during a hydraulic test in fractured rock. In this study, we evaluated the possibility of distortion of the estimated hydraulic properties from a hydraulic test due to nonlinear flow. We conducted pulse, slug, and constant head withdrawal tests with various head displacements using a double packer system in a borehole at Korea Atomic Energy Research Institute (KAERI) Underground Research Tunnel (URT), hereafter KURT. With the geophysical logging data for the test borehole and the monitored hydraulic heads during the tests, changes of the representative Re were calculated for identifying the nonlinear flow regimes during the tests. Considering the influence of nonlinear flow, distortion of the estimated hydraulic conductivity from pulse, slug, and constant head withdrawal tests were discussed.

2. Materials and Methods

2.1. Study Site and Hydraulic Tests

The study site was KURT, which is a small-scale underground research facility for radioactive waste disposal, located in Daejeon, South Korea (Figure 1). The host rock of KURT is Mesozoic two-mica granite. It was tunneled from the side of mountain with a dipping angle of about 10%, and the maximum depth of KURT is about 120 m below the surface. It has 6 research galleries (RG), where various in situ tests have been conducted.
The test borehole, BDZ-E-1, is located at RG-6 (Figure 1). It is 20 m long and vertical borehole with a diameter of 10.2 cm (4 inches), and it is not screened. The intersecting fractures were logged through the acoustic televiewing. In 12 packed-off sections with a length of 1.45 m, pulse, slug, and constant head withdrawal tests were conducted with various head displacements. The static groundwater levels at the packed-off sections were as high as ~5.3 m to ~23.7 m from the top of the casing (TOC) (Table 1), and under artesian conditions. Pulse, slug, and constant head withdrawal tests all require a sudden change in hydraulic head at the start of the test, and Figure 2 shows the special packer system allowing isolation of the test zone and abrupt change in water pressure for the tests under artesian condition. In pulse tests, the shut-in valve was initially closed to sever the hydraulic connection between the test zone and the packer access pipe, and the water level in the packer access pipe was adjusted considering the magnitude of initial head displacement ( H 0 ). After stabilizing the hydraulic head in the packed-off test zone, the shut-in valve was opened for several seconds to apply the initial head displacement to the test zone, and then it was closed. The change in the volume of water occurred while the shut-in valve was opened was estimated with the measured water level change in the packer access pipe, and the recovery of the hydraulic head in the test zone was monitored using a transducer located in the chamber between the shut-in valve and upper packer. Slug tests proceeded the same procedures as pulse tests except that the shut-in valve remained open after applying the initial head displacement. In constant head withdrawal tests, the shut-in valve kept open. After the hydraulic head in the test zone became static, water was bled off using the valve for discharging water to maintain the specified head displacement ( s w ), and the leakage rates were measured. Note that the hydraulic heads in the test zone were measured at 1 s intervals during all tests. In each test zones, the pulse and slug tests were conducted 3–5 times each with various initial head displacements, while the constant head withdrawal test was performed only once due to the long test period.

2.2. Estimation of Hydraulic Conductivities

To evaluate distortion of the estimated hydraulic conductivity by nonlinear groundwater flow, the flow regimes at fractures during the hydraulic tests were identified by calculating representative Re following Ji and Koh [9]. The Re at a fracture is defined as
Re = ρ ve μ = ρ Q w μ   ,
where ρ (M/L3) is the fluid density, v (L/T) is the flow velocity, e (L) is the fracture aperture, μ (M/L T) is the fluid viscosity, Q (L3/T) is the flow rate, and w (L) is the fracture width perpendicular to flow. Following Ji and Koh [9], we assumed that the open fractures can be hydraulically active, and the number and characteristics of the open fractures from the acoustic televiewer logging data were used to calculate the representative Re for the fractures in the test zone. Table 1 includes the characteristics of open fractures used for calculation of the representative Re. w was assumed as the circumferences of the traces of open fractures at the borehole wall, and the mean fracture width ( w m ) for a fracture in the test zone was estimated as
w m = i N f π ( r w cos θ i + r w ) N f   ,
where N f is the number of open fractures in the test zone, r w (L) is the radius of the borehole, and θ i is the dipping angle of the i -th fracture in the test zone. Then, the mean flow rate ( Q m ) for calculating the representative Re was calculated with
Q m = i N f Q i N f   ,
where Q i (L3/T) is the flow rate of the i -th fracture in the test zone. The total flow rate in Equation (3) was estimated differently depending on the type of test while it was directly measured in the constant head withdrawal tests. It was given by
i N f Q i = π r e 2 Δ h Δ t   ,
i N f Q i = π r p 2 Δ h Δ t   ,
for pulse and slug tests, respectively, where r e (L) is the effective radius for pulse test analysis given by r e = Δ V w π H 0 [11], r p (L) is the radius of the packer access pipe, Δ h (L) is the change in the hydraulic head in the test zone during time Δ t (T), and Δ V w (L3) is the change in the volume of water in the packer access pipe occurred while the shut-in valve was opened for pulse extraction during a pulse test.
The results of pulse and slug tests were analyzed with the analytical methods using the Bouwer–Rice model [12] and Cooper–Bredehoeft–Papadopulos model for a slug test in a confined aquifer [13,14], and those of constant head withdrawal tests were the Jacob–Lohman model [15] and straight line model [16]. The hydraulic conductivity of the test zone was estimated by minimizing the residuals between the observed data and the type curves from the semi-analytical solution of each model with the Gauss–Newton optimization algorithm [17]. Note that the residual is defined as
Residual = i n ( y i - y i ^ ) 2 n - p   ,
where y i and y i ^ are the i -th observed data and estimate of y i from the type curve, n is the number of observations, and p is the number of estimated hydraulic parameters. Then, the hydraulic conductivity of the test zone was determined as the one of the models showing small residual among the Bouwer–Rice and Cooper–Bredehoeft–Papadopulos models for the pulse and slug tests, and the Jacob–Lohman and straight line models for the constant head withdrawal tests.
The Forchheimer equation has been used to describe the nonlinear flow in an aquifer including a single fracture, e.g., [18,19]. To clarify the nonlinear flow regime, Zimmerman et al. [20] coupled the cubic law to the Forchheimer equation, and the coupled equation is given by
K = K 0 1 + α Re   ,
where K (L/T) is the hydraulic conductivity, K 0 (L/T) is the hydraulic conductivity at a low Re in a linear flow regime, α (dimensionless) is a constant whose definition is α = 4 π 2 r w 2 eb μ K 0 / ρ , and b (T2/L6) is a constant. Using Equation (7), the hydraulic conductivity of the test zone in a linear flow regime was suggested following Ji and Koh [9], and distortion of the estimated K from various hydraulic tests was evaluated.

3. Results and Discussion

Pulse, slug, and constant head withdrawal tests were conducted in each test zone except several sections due to technical issues. Figure 3 shows the results of pulse, slug, and constant head withdrawal tests conducted at section 3 as an example. Note that s w and Q in Figure 3c are the specified head displacement for a constant head withdrawal test and measured outflow rate during the test, respectively. Using the logged data on open fractures, the monitored hydraulic heads and flow rates, and then the representative Re at a hydraulically conductive fracture in the test zone during the tests, were estimated using Equations (1)–(5) as shown in Figure 3. Regardless of the type of test, the representative Re were initially maximum and then gradually decreased, as were the head displacements and flow rates. This pattern of change in the representative Re was the same for all test sections. This is because the representative Re is proportional to the flow rate as in Equation (1).
To estimate K of the test zones, the results of the conducted hydraulic tests were analyzed using the type curves from the semi-analytical solutions for pulse, slug, and constant head withdrawal tests. Figure 3 indicates the optimal matched type curves for the results of pulse, slug, and constant head withdrawal tests at section 3 as an example. Most of pulse and slug tests were optimally matched with the Cooper–Bredehoeft–Papadopulos model except sections 4, 6, and 11. The Bouwer–Rice model for pulse and slug tests was developed from the Thiem equation which assumes the steady state groundwater flow [12], but the Cooper–Bredehoeft–Papadopulos model was induced from the partial differential equation describing the transient groundwater flow in a confined aquifer [14]. Thus, the optimization results of pulse and slug tests with the type curves show that the effect of storage was relatively large in the most test zones. Most of constant head withdrawal tests were matched well with the straight line model, which means that the storage coefficients of the test zones were enough small to approximate the well function, W ( u ) , to 2.30 log 10 2.25 Tt / r w 2 S , where u (dimensionless) is a variable defined as u = r w 2 S / 4 Tt , T (L2/T) is the transmissivity, t (T) is the elapsed time since the test began, and S (dimensionless) is the storage coefficient [16].
Figure 4 shows the residuals between the observed data and optimized type curves. Generally, the residuals of pulse tests were 2 to 9 times larger than those of slug tests, depending on the test zones. This result indicates that groundwater flow during a pulse test was more likely to deviate from the type curve than a slug test and the uncertainty in interpretation of the pulse test was greater than the slug test. Since each test zone has different hydraulic characteristics, it is meaningless to evaluate the relation between the residual and maximum Re during hydraulic tests without considering the test zone. In pulse tests, there was no clear correlation when the relation between the residual and maximum Re for each test zone was examined (Figure 4a). In slug tests, however, the residual generally increased at a larger maximum Re for several test zones (Figure 4b), which means that groundwater flow during a slug test deviated from the theoretical behavior expressed by the type curve with a linear groundwater flow assumption. From the results of the constant head withdrawal tests, the relation between the residual and the maximum Re could not be evaluated because the constant head withdrawal test was conducted only once for each test zone (Figure 4c).
The relations between the maximum Re and estimated K from pulse, slug, and constant head withdrawal tests for each test zone were shown in Figure 5. For pulse tests, the estimated K for each test zone showed similar values, and there was no clear correlation between the maximum Re and estimated K . In section 10, the estimated K seems to decrease as the maximum Re increases, but it is not clear because the range of maximum K is small (Figure 5i). The maximum Re during pulse tests were less than ~7, except section 10. Accordingly, it seems that the nonlinear flow did not occur during the pulse tests. This is because a pulse test is a hydraulic test that utilize pressure changes rather than actual water movement and the flow rates were very small despite the large H 0 .
For slug tests, it was difficult to see the correlation between the estimated K and maximum Re in sections 3–5, 7, 9, and 12 (Figure 5b–d,f,h,k). The maximum Re of those test zones were below 5 except section 12. Considering the reported critical Re and the graphs showing the estimated K and maximum Re, the groundwater flow regimes during slug tests at those test zones were likely to be linear ones, and the probability that the estimated K were distorted by nonlinear flow was small. In sections 1, 6, 8, 10, and 11, however, the estimated   K decreased with an increase in the maximum Re, and the relations between the estimated K and maximum Re were fitted to Equation (7) (Figure 5a,e,g,i,j). This shows that nonlinear flows were likely to have effects on estimating K of sections 1, 6, 8, 10, and 11 from slug tests. From the fitted curves, K 0 of sections 1, 6, 8, 10, and 11 were optimized to 2.1 × 10 8 m/s, 5.5 × 10 8 m/s, 6.0 × 10 8 m/s, 2.5 × 10 7 m/s, and 9.0 × 10 8 m/s, respectively. The critical Re for slug tests in those test zones were ranged from ~1.1 to ~5.9. Note that the reported critical Re varied in a range from 0.001 to 24.8, which was suggested through numerical, laboratory and field experimental analyses: 0.001–25 [3], 1 [5], 5 [21], 3–12 [9], and 3.5–24.8 [22].
Because the water leakage valve for draining water is located on the ground as in Figure 2, small s w were applied to the constant head withdrawal tests. Thus, the maximum Re during the constant head withdrawal tests at all test zones except section 10 were smaller than the critical Re for slug tests, but it was difficult to discuss the effect of nonlinear flow on the constant head withdrawal tests in our test borehole because it was conducted only once at each test zone.
When the permeability of a test zone was low like sections 3 and 4 (Figure 5b,c), the estimated K from pulse, slug, and constant head withdrawal tests were similar to each other. This may be because the radii of influence for pulse, slug, and constant head withdrawal tests were similar due to the low permeability. In other cases, the estimated K from constant head withdrawal tests generally tended to be the largest, those from slug tests were the next largest, and those from pulse tests were the smallest. This trend was evident in the test zones with high permeability such as sections 10 and 11 (Figure 5i,j). For pulse, slug, and constant head withdrawal tests, the test durations were different as shown in Figure 3. This difference in test duration may cause a difference in the radii of influence for pulse, slug, and constant head withdrawal tests, and thus the estimated K . These results were consistent with the reports that the hydraulic conductivity tends to increase with the involved rock volume, e.g., [23,24,25].

4. Conclusions

Distortion of the estimated K by excessive injection or extraction of groundwater was discussed by evaluating the influence of nonlinear flow on hydraulic tests. Pulse, slug, and constant head withdrawal tests were conducted at an installed borehole in KURT, and the changes of representative Re at a hydraulically conductive fracture during the hydraulic tests were calculated. Then, the relations between the maximum Re during the tests and estimated K from the tests were matched to an equation introduced from the Forchheimer equation and cubic law, and the matching results were used to identify the effects of nonlinear flow on the tests and to evaluate the possibility of distortion of the estimated K . Our results showed that nonlinear flow was unlikely to occur during the pulse tests in our test zones and there was little possibility that the estimated K from pulse tests distorted from nonlinear flow. In slug tests, however, the relations between the maximum Re and estimated K indicated that estimation of hydraulic parameters for several test zones was likely to be distorted due to nonlinear flow during the tests, and the estimated K for those zones were corrected following Ji and Koh [9]. Then, the estimated K from pulse, slug, and constant head withdrawal tests were compared with each other. The comparison results showed that the scale effect was generally invoked in estimation of K for our test zones and the estimated K increased as the test scale grew from pulse tests to constant head withdrawal tests. When the test zones had low permeability, however, the scale effect was reduced and the estimated K from the different types of hydraulic tests became similar.
The hydraulic properties obtained through site characterization are reflected in the safety assessment for a subsurface HLW repository, and the safety can be overestimated if the hydraulic conductivity of the disposal site is underestimated. Nonlinear flow from excessive injection or extraction can induce underestimation of the hydraulic conductivity. Our study indicated the way to evaluate the effects of nonlinear flow on a hydraulic test and distortion of the hydraulic conductivity estimated from the test. Considering the nonlinear flow and scale effect, our results showed that a pulse test is recommended for the zones with low permeability in our study site. In other cases, however, a constant head withdrawal test is proposed for the conservative safety assessment.

Author Contributions

Conceptualization, S.-H.J.; methodology, S.-H.J.; formal analysis, S.-H.J. and B.-H.P.; investigation, B.-H.P. and K.-W.P.; writing—original draft preparation, S.-H.J.; writing—review and editing, S.-H.J., B.-H.P., and K.-W.P.; and visualization, S.-H.J. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Nuclear Research and Development Program of the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT (No. 2017M2A8A5014858).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ReReynolds number
HLWHigh-Level radioactive Waste
KURTKorea Atomic Energy Research Institute Underground Research Tunnel
RGResearch Gallery
TOCTop Of the Casing

References

  1. Ranjith, P.G.; Darlington, W. Nonlinear single-phase flow in real rock joints. Water Resour. Res. 2007, 43, W09502. [Google Scholar] [CrossRef] [Green Version]
  2. Ji, S.-H.; Lee, H.-B.; Yeo, I.W.; Lee, K.-K. Effect of nonlinear flow on DNAPL migration in a rough-walled fracture. Water Resour. Res. 2008, 44, W11431. [Google Scholar] [CrossRef]
  3. Javadi, M.; Sharifzadeh, M.; Shahriar, K.; Mitani, Y. Critical Reynolds number for nonlinear flow through rough-walled fractures: The role of shear processes. Water Resour. Res. 2014, 50, 1789–1804. [Google Scholar] [CrossRef] [Green Version]
  4. Zimmerman, R.W.; Bodvarsson, G.S. Hydraulic conductivity of rock fractures. Transp. Porous Media 1996, 23, 1–30. [Google Scholar] [CrossRef] [Green Version]
  5. Brush, D.J.; Thomson, N.R. Fluid flow in synthetic rough-walled fractures: Navier-Stokes, Stokes, and local cubic law simulations. Water Resour. Res. 2003, 39, 1085. [Google Scholar] [CrossRef]
  6. Zou, L.; Jing, L.; Cvetkovic, V. Roughness decomposition and nonlinear fluid flow in a single rock fracture. Int. J. Rock Mech. Min. Sci. 2015, 75, 102–118. [Google Scholar] [CrossRef]
  7. Quinn, P.M.; Cherry, J.A.; Parker, B.L. Quantification of non-Darcian flow observed during packer testing in fractured sedimentary rock. Water Resour. Res. 2011, 47, W09533. [Google Scholar] [CrossRef]
  8. Quinn, P.M.; Parker, B.L.; Cherry, J.A. Validation of non-Darcian flow effects in slug tests conducted in fractured rock boreholes. J. Hydrol. 2013, 486, 505–518. [Google Scholar] [CrossRef]
  9. Ji, S.-H.; Koh, Y.-K. Nonlinear groundwater flow during a slug test in fractured rock. J. Hydrol. 2015, 520, 30–36. [Google Scholar] [CrossRef]
  10. Chen, Y.-F.; Hu, S.-H.; Hu, R.; Zhou, C.-B. Estimating hydraulic conductivity of fractured rocks from high-pressure packer tests with an Izbash’s law-based empirical model. Water Resour. Res. 2015, 51, 2096–2118. [Google Scholar] [CrossRef]
  11. Steiner, W.; Thut, A.; Fisch, H.; Gysi, H.-J. Geohydraulische Versuche in Fels; Eidgenössisches Departement für Umwelt, Verkehr, Energie und Kommunikation: Bern, Switzerland, 2006; pp. 37–39. [Google Scholar]
  12. Bouwer, H. The Bouwer and Rice slug test—An update. Groundwater 1989, 27, 304–309. [Google Scholar] [CrossRef]
  13. Cooper, H.H.; Bredehoeft, J.D.; Papadopulos, I.S. Response of a finite-diameter well to an instantaneous change of water. Water Resour. Res. 1967, 3, 263–269. [Google Scholar] [CrossRef]
  14. Bredehoeft, J.D.; Papadopulos, S.S. A method for determining the hydraulic properties of tight formations. Water Resour. Res. 1980, 16, 233–238. [Google Scholar] [CrossRef]
  15. Jacob, C.E.; Lohman, S.W. Nonsteady flow to a well of constant drawdown in an extensive aquifer. Transactions 1952, 33, 559–569. [Google Scholar] [CrossRef]
  16. Lohman, S.W. Ground-Water Hydraulic, Geological Survey Professional Paper 708; United States Government Printing Office: Washington, DC, USA, 1972; pp. 23–27.
  17. Marquardt, D.W. An algorithm for least-squares estimation of nonlinear parameters. J. Soc. Ind. Appl. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
  18. Jacob, C.E. Effective radius of drawdown test to determine artesian well. Proc. Am. Soc. Civ. Eng. 1946, 72, 629–646. [Google Scholar]
  19. McElwee, C.D.; Zenner, M.A. A nonlinear model for analysis of slug-test data. Water Resour. Res. 1998, 34, 55–66. [Google Scholar] [CrossRef]
  20. Zimmerman, R.W.; Al-Yaarubi, A.; Pain, C.C.; Grattoni, C.A. Non-linear regimes of fluid flow in rock fractures. Int. J. Rock Mech. Min. Sci. 2004, 41, 384. [Google Scholar] [CrossRef]
  21. Konzuk, J.S.; Kueper, B.H. Evaluation of cubic law based models describing single-phase flow through a rough-walled fracture. Water Resour. Res. 2004, 40, W02402. [Google Scholar] [CrossRef]
  22. Zhang, Z.; Nemcik, J. Fluid flow regimes and nonlinear flow characteristics in deformable rock fractures. J. Hydrol. 2013, 477, 139–151. [Google Scholar] [CrossRef]
  23. Illman, W.A.; Neuman, S.P. Steady-state analysis of cross-hole pneumatic injection tests in unsaturated fractured tuff. J. Hydrol. 2003, 281, 54–72. [Google Scholar] [CrossRef]
  24. Illman, W.A. Analysis of permeability scaling within single boreholes. Geophys. Res. Lett. 2004, 31, L06503. [Google Scholar] [CrossRef]
  25. Martinez-Landa, L.; Carrera, J. An analysis of hydraulic conductivity scale effects in granite (Full-scale Engineered Barrier Experiment (FEBEX), Grimsel, Switzerland). Water Resour. Res. 2005, 41, W03006. [Google Scholar]
Figure 1. (a) Location; and (b) layout of the Korea Atomic Energy Research Institute (KAERI) Underground Research Tunnel (KURT). The test borehole, BDZ-E-1 is located at RG-6 in KURT.
Figure 1. (a) Location; and (b) layout of the Korea Atomic Energy Research Institute (KAERI) Underground Research Tunnel (KURT). The test borehole, BDZ-E-1 is located at RG-6 in KURT.
Water 12 02712 g001
Figure 2. Double packer system for pulse, slug, and constant head withdrawal tests at the test zone under artesian condition.
Figure 2. Double packer system for pulse, slug, and constant head withdrawal tests at the test zone under artesian condition.
Water 12 02712 g002
Figure 3. An example of monitored data and calculated representative Re for the (a) pulse; (b) slug; and (c) constant head withdrawal tests at section 3. The red solid lines are the optimized type curves.
Figure 3. An example of monitored data and calculated representative Re for the (a) pulse; (b) slug; and (c) constant head withdrawal tests at section 3. The red solid lines are the optimized type curves.
Water 12 02712 g003
Figure 4. Relation between the residual and maximum Re for the (a) pulse; (b) slug; and (c) constant head withdrawal tests.
Figure 4. Relation between the residual and maximum Re for the (a) pulse; (b) slug; and (c) constant head withdrawal tests.
Water 12 02712 g004
Figure 5. Relation between the estimated K and maximum Re at (a) section 1; (b) section 3; (c) section 4; (d) section 5; (e) section 6; (f) section 7; (g) section 8; (h) section 9; (i) section 10; (j) section 11; and (k) section 12. The solid and dashed lines are the fitted curves using Equation (7) to evaluate the influence of nonlinear flow on the tests and discuss distortion of the estimated K from the tests.
Figure 5. Relation between the estimated K and maximum Re at (a) section 1; (b) section 3; (c) section 4; (d) section 5; (e) section 6; (f) section 7; (g) section 8; (h) section 9; (i) section 10; (j) section 11; and (k) section 12. The solid and dashed lines are the fitted curves using Equation (7) to evaluate the influence of nonlinear flow on the tests and discuss distortion of the estimated K from the tests.
Water 12 02712 g005
Table 1. Locations of the test zones and the static groundwater levels and intersecting open fractures in the test zones.
Table 1. Locations of the test zones and the static groundwater levels and intersecting open fractures in the test zones.
Test ZoneDepth From TOC (m)Static Water Level Above TOC (m)Number of Intersecting Open FracturesOrientations of Intersecting Open Fractures
Section 118.00–19.4523.73(205.1, 70.0) 1
(261.6, 55.4)
(266.0, 73.7)
Section 216.55–18.008.70-
Section 315.10–16.5523.31(88.7, 85.9)
Section 413.65–15.1016.03(63.0, 33.3)
(274.0, 86.0)
(0, 0)
Section 512.20–13.656.01(102.9, 87.6)
Section 610.75–12.205.53(28.4, 31.2)
(254.5, 68.2)
(251.8, 58.9)
Section 79.30–10.7513.32(262.5, 8.2)
(87.8, 17.5)
Section 87.85–9.3010.84(95.8, 63.4)
(94.9, 83.1)
(93.1, 84.2)
(157.8, 17.5)
Section 96.40–7.858.03(16.0, 35.4)
(262.5, 49.8)
(32.8, 16.1)
Section 104.95–6.406.61(180.0, 63.3)
Section 113.50–4.9510.12(206.6, 51.0)
(182.7, 67.6)
Section 122.05–3.505.31(117.0, 68.6)
1 (Dip direction, Dipping angle).

Share and Cite

MDPI and ACS Style

Ji, S.-H.; Park, B.-H.; Park, K.-W. Distortion of the Estimated Hydraulic Conductivity from a Hydraulic Test in Fractured Rock Due to Excessive Injection or Extraction. Water 2020, 12, 2712. https://doi.org/10.3390/w12102712

AMA Style

Ji S-H, Park B-H, Park K-W. Distortion of the Estimated Hydraulic Conductivity from a Hydraulic Test in Fractured Rock Due to Excessive Injection or Extraction. Water. 2020; 12(10):2712. https://doi.org/10.3390/w12102712

Chicago/Turabian Style

Ji, Sung-Hoon, Byeong-Hak Park, and Kyung-Woo Park. 2020. "Distortion of the Estimated Hydraulic Conductivity from a Hydraulic Test in Fractured Rock Due to Excessive Injection or Extraction" Water 12, no. 10: 2712. https://doi.org/10.3390/w12102712

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

Article Metrics

Back to TopTop