Difference between revisions of "Characterization Methods – Hydraulic Conductivity"

From Enviro Wiki
Jump to: navigation, search
m (correcting references)
 
(13 intermediate revisions by 2 users not shown)
Line 9: Line 9:
  
  
'''CONTRIBUTOR(S):''' [[Dr. Gaisheng Liu]] and [[Dr. James Butler, Jr.]]
+
'''Contributor(s):''' [[Dr. Gaisheng Liu]] and [[Dr. James Butler, Jr.]]
  
  
Line 19: Line 19:
  
 
==Introduction==
 
==Introduction==
Hydraulic conductivity is mathematically defined as the parameter ''K'' in Darcy’s Law<ref>Darcy, H. (1856). Les Fontaines Publiques de la Ville de Dijon, Dalmont, Paris. [https://doi.org/10.1029/2001wr000727 doi: 10.1029/2001WR000727]</ref> (see Fig. 1),
+
Hydraulic conductivity is mathematically defined as the parameter ''K'' in Darcy’s Law<ref>Darcy, H. (1856). Les Fontaines Publiques de la Ville de Dijon, Dalmont, Paris. [https://doi.org/10.1029/2001wr000727 doi: 10.1029/2001WR000727]</ref> (see Figure 1),
[[File:Liu1w2 Fig1.png|center|thumb|400|Figure 1. Schematic of Darcy’s Law flow experiment.]]
 
 
 
 
 
[[File:Liu1w2 Eq1.png|center|300 px]]
 
  
 +
[[File:Liu1w2 Fig1.png|thumb|500 px|left|Figure 1. Schematic of Darcy’s Law flow experiment.]]
 
{|
 
{|
 
|-
 
|-
|where:
+
| ||[[File:Liu1w2 Eq1.png|310 px]]
 +
|-
 +
|where:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 
|-
 
|-
 
|''Q''||is the flow rate across area ''A'' of a porous medium,
 
|''Q''||is the flow rate across area ''A'' of a porous medium,
Line 36: Line 35:
 
|}
 
|}
  
 
[[File:Liu1w2 Eq2.png|center|300 px]]
 
  
 
{|
 
{|
 +
|-
 +
| ||[[File:Liu1w2 Eq2.png|290 px]]
 
|-
 
|-
 
|where:
 
|where:
 
|-
 
|-
|''h<sub>1</sub>'' and ''h<sub>2</sub>''  &nbsp; &nbsp;||are the hydraulic heads at the ends of the experimental domain, and
+
|''h<sub>1</sub>'' and ''h<sub>2</sub>''  &nbsp; &nbsp;||are the [[wikipedia: Hydraulic head | hydraulic heads]] at the ends of the experimental domain, and
 
|-
 
|-
 
|''L''||is the total length.
 
|''L''||is the total length.
Line 50: Line 49:
 
Hydraulic conductivity is dependent on the properties of both water and the porous medium,
 
Hydraulic conductivity is dependent on the properties of both water and the porous medium,
  
[[File:Liu1w2 Eq3.png|center|300 px]]
 
  
 
{|
 
{|
 
|-
 
|-
|where:
+
| ||[[File:Liu1w2 Eq3.png|300 px]]
 +
|-
 +
|where:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 
|-
 
|-
 
|''k''||is the intrinsic permeability of the medium, a parameter which is solely dependent on the geometry of the interconnected pores,
 
|''k''||is the intrinsic permeability of the medium, a parameter which is solely dependent on the geometry of the interconnected pores,
 
|-
 
|-
|''g''||is the gravitational constant,
+
|''g''||is the [[wikipedia: Gravitational constant | gravitational constant]],
 
|-
 
|-
 
|''&rho;''||is the density of the pore water, and
 
|''&rho;''||is the density of the pore water, and
Line 67: Line 67:
 
Many approaches have been developed to characterize ''K''<ref name="Butler2005" />. These approaches can be grouped into two general categories based on how the ''K'' estimates are obtained:  
 
Many approaches have been developed to characterize ''K''<ref name="Butler2005" />. These approaches can be grouped into two general categories based on how the ''K'' estimates are obtained:  
  
*hydraulic methods that involve water injection or extraction and the measurement of the induced pressure response, and
+
*hydraulic methods that involve water or other fluid injection or extraction and the measurement of the induced pressure response, and
 
*indirect methods that rely on empirical correlations, often site-specific in nature, between ''K'' and other more readily evaluated formation properties (e.g., resistance to electric current).
 
*indirect methods that rely on empirical correlations, often site-specific in nature, between ''K'' and other more readily evaluated formation properties (e.g., resistance to electric current).
  
Line 74: Line 74:
 
Other classifications of ''K'' characterization approaches are possible. For example, approaches can be divided into those based on data collected in the field and those based on measurements on core samples in the laboratory. In the latter case, there can be considerable uncertainty about how representative the core sample is of field conditions.  
 
Other classifications of ''K'' characterization approaches are possible. For example, approaches can be divided into those based on data collected in the field and those based on measurements on core samples in the laboratory. In the latter case, there can be considerable uncertainty about how representative the core sample is of field conditions.  
  
Field-based approaches can be further subdivided based on whether the measurement is conducted in the subsurface or on the surface. Subsurface measurements provide the most reliable information about ''K'' variations at the scale needed for environmental site investigations. Nearly all field-based hydraulic methods are performed in the subsurface and require either boreholes or direct push (DP) installations for extracting or injecting water and monitoring the induced head changes. Compared to methods that rely on existing wells, DP approaches can be applied essentially wherever ''K'' information is needed, thus greatly expanding the spatial coverage of ''K'' measurement in the field. DP approaches, however, are generally limited to use in relatively shallow (20-30 m from land surface) unconsolidated settings.
+
Field-based approaches can be further subdivided based on whether the measurement is conducted in the subsurface or on the surface. Subsurface measurements provide the most reliable information about ''K'' variations at the scale needed for environmental site investigations. Nearly all field-based hydraulic methods are performed in the subsurface and require either boreholes or [[Direct Push (DP) Technology | direct push (DP)]] installations for extracting or injecting water and for monitoring the induced head changes. Compared to methods that rely on existing wells, DP approaches can be applied essentially wherever ''K'' information is needed, thus greatly expanding the spatial coverage of ''K'' measurement in the field. DP approaches, however, are generally limited to use in relatively shallow (20-30 m from land surface) unconsolidated settings.
  
 
Field ''K'' measurements can also be performed on the land surface or in an aircraft. Most of these are geophysical approaches that involve measuring electrical, magnetic, or seismic properties of the formation and then transforming them into ''K'' estimates using empirical relations. The biggest advantage of surface or airborne approaches is their minimal site disturbance, as no subsurface boring is needed. However, the resolution that can be provided is limited, as the measurements are typically affected by conditions over a relatively large volume of the formation. Geophysical methods can also be performed in boreholes. Borehole geophysical methods can provide a much higher resolution description of ''K'' than surface-based methods.
 
Field ''K'' measurements can also be performed on the land surface or in an aircraft. Most of these are geophysical approaches that involve measuring electrical, magnetic, or seismic properties of the formation and then transforming them into ''K'' estimates using empirical relations. The biggest advantage of surface or airborne approaches is their minimal site disturbance, as no subsurface boring is needed. However, the resolution that can be provided is limited, as the measurements are typically affected by conditions over a relatively large volume of the formation. Geophysical methods can also be performed in boreholes. Borehole geophysical methods can provide a much higher resolution description of ''K'' than surface-based methods.
Line 80: Line 80:
 
The scale of measurement is an important factor to consider when assessing various approaches. For water-supply investigations, a single estimate of ''K'' averaged over a large volume of an aquifer will usually suffice; such an estimate is commonly obtained via pumping tests in a well<ref name="Freeze1979" />. For water-quality investigations, however, ''K'' estimates over a large measurement volume are often of limited value<ref name="Butler2009">Butler Jr, J.J., 2009. Pumping tests for aquifer evaluation - Time for a change?. Groundwater, 47(5), pp.615-617. [https://doi.org/10.1111/j.1745-6584.2008.00488.x doi: 10.1111/j.1745-6584.2008.00488.x]</ref>. In that case, small-scale ''K'' measurements that provide information about local geological controls on groundwater flow and transport are usually required to obtain reliable predictions of contaminant behavior and to design effective remediation systems.  
 
The scale of measurement is an important factor to consider when assessing various approaches. For water-supply investigations, a single estimate of ''K'' averaged over a large volume of an aquifer will usually suffice; such an estimate is commonly obtained via pumping tests in a well<ref name="Freeze1979" />. For water-quality investigations, however, ''K'' estimates over a large measurement volume are often of limited value<ref name="Butler2009">Butler Jr, J.J., 2009. Pumping tests for aquifer evaluation - Time for a change?. Groundwater, 47(5), pp.615-617. [https://doi.org/10.1111/j.1745-6584.2008.00488.x doi: 10.1111/j.1745-6584.2008.00488.x]</ref>. In that case, small-scale ''K'' measurements that provide information about local geological controls on groundwater flow and transport are usually required to obtain reliable predictions of contaminant behavior and to design effective remediation systems.  
  
In this article, we focus on the approaches that can be used to obtain small-scale, localized measurements of ''K'' for environmental site investigations. Only the more commonly used approaches are discussed. Some other approaches, such as dipole flow tests<ref>Kabala, Z.J., 1993. The dipole flow test: A new single‐borehole test for aquifer characterization. Water Resources Research, 29(1), pp.99-107. [https://doi.org/10.1029/92WR01820 doi: 10.1029/92WR01820]</ref><ref>Zlotnik, V.A. and Zurbuchen, B.R., 1998. Dipole probe: Design and field applications of a single‐borehole device for measurements of vertical variations of hydraulic conductivity. Groundwater, 36(6), pp.884-893. [https://doi.org/10.1111/j.1745-6584.1998.tb02095.x doi: 10.1111/j.1745-6584.1998.tb02095.x]</ref>, tracer tests<ref>Datta‐Gupta, A., Yoon, S., Vasco, D.W. and Pope, G.A., 2002. Inverse modeling of partitioning interwell tracer tests: A streamline approach. Water Resources Research, 38(6), pp.15-1. [https://doi.org/10.1029/2001WR000597 doi: 10.1029/2001WR000597]</ref>, or hydraulic tomography<ref>Yeh, T.C.J. and Liu, S., 2000. Hydraulic tomography: Development of a new aquifer test method. Water Resources Research, 36(8), pp.2095-2105. [//www.enviro.wiki/images/a/af/2000-Yeh-Hydraulic_Tomography..Development_of_a_new_aquifer_test_method.pdf Report.pdf]</ref><ref>Bohling, G.C. and Butler Jr, J.J., 2010. Inherent limitations of hydraulic tomography. Groundwater, 48(6), pp.809-824. [https://doi.org/10.1111/j.1745-6584.2010.00757.x doi: 10.1111/j.1745-6584.2010.00757.x]</ref> will not be discussed here.
+
In this article, we focus on the approaches that can be used to obtain small-scale, localized measurements of ''K'' for environmental site investigations. Only the more commonly used approaches are discussed. Some other approaches, such as dipole flow tests<ref>Kabala, Z.J., 1993. The dipole flow test: A new single‐borehole test for aquifer characterization. Water Resources Research, 29(1), pp.99-107. [https://doi.org/10.1029/92WR01820 doi: 10.1029/92WR01820]</ref><ref>Zlotnik, V.A. and Zurbuchen, B.R., 1998. Dipole probe: Design and field applications of a single‐borehole device for measurements of vertical variations of hydraulic conductivity. Groundwater, 36(6), pp.884-893. [https://doi.org/10.1111/j.1745-6584.1998.tb02095.x doi: 10.1111/j.1745-6584.1998.tb02095.x]</ref>, tracer tests<ref>Datta‐Gupta, A., Yoon, S., Vasco, D.W. and Pope, G.A., 2002. Inverse modeling of partitioning interwell tracer tests: A streamline approach. Water Resources Research, 38(6), pp.15-1. [https://doi.org/10.1029/2001WR000597 doi: 10.1029/2001WR000597]</ref>, or hydraulic tomography<ref>Yeh, T.C.J. and Liu, S., 2000. Hydraulic tomography: Development of a new aquifer test method. Water Resources Research, 36(8), pp.2095-2105. [//www.enviro.wiki/images/a/af/2000-Yeh-Hydraulic_Tomography..Development_of_a_new_aquifer_test_method.pdf Report.pdf]</ref><ref>Bohling, G.C. and Butler Jr, J.J., 2010. Inherent limitations of hydraulic tomography. Groundwater, 48(6), pp.809-824. [https://doi.org/10.1111/j.1745-6584.2010.00757.x doi: 10.1111/j.1745-6584.2010.00757.x]</ref> will not be discussed here.
  
 
==Hydraulic Testing Approaches==
 
==Hydraulic Testing Approaches==
Line 90: Line 90:
  
 
===Slug Tests===  
 
===Slug Tests===  
 +
[[File:Liu1w2 Fig2A.png|thumb|280 px|left|Figure 2A. Schematic of a Falling-Head Slug Test (water flows from well into aquifer) Performed in a Multi-Level Slug-Test System. Head change is introduced in standpipe, which is directly connected to the interval isolated by the straddle packers (not to scale).]] 
 +
[[File:Liu1w2 Fig2B.png|thumb|400 px|right|Figure 2B. Example Data plot from a Multilevel Slug Test at GEMS. Test is initiated by sudden depressurization of a pressurized air column<ref name="Butler2003" />.]]
 +
 
The slug test is one of the most common methods for determining ''K'' at the scale of relevance for environmental site investigations<ref name="Butler1997">Butler Jr, J.J., 1997. The design, performance, and analysis of slug tests. CRC Press.</ref>. Slug tests are typically performed in existing wells at multiple horizontal and vertical locations across the site to characterize the spatial distribution of ''K'' values. However, slug tests can also be performed in temporary installations such as DP rods. In this approach, a near-instantaneous head change is applied in the well or DP rods and the subsequent head recovery is used to estimate ''K''. The initial head change can be introduced using a solid object (slug), compressed gas (pneumatic system), or by the addition/removal of a certain amount of water<ref name="Butler1997" />. The recovery data can be analyzed using different models of the well-formation configuration<ref name="Duffield2007" /><ref name="Butler1997" />. In contrast to pumping tests, slug tests provide a ''K'' estimate that is primarily a function of the materials in the vicinity of the screened interval of the test well. Slug test data analyses can be facilitated using software packages like AQTESOLV (Aquifer Test Solver), which has been developed for analyzing different test methods<ref name="Duffield2007" />.
 
The slug test is one of the most common methods for determining ''K'' at the scale of relevance for environmental site investigations<ref name="Butler1997">Butler Jr, J.J., 1997. The design, performance, and analysis of slug tests. CRC Press.</ref>. Slug tests are typically performed in existing wells at multiple horizontal and vertical locations across the site to characterize the spatial distribution of ''K'' values. However, slug tests can also be performed in temporary installations such as DP rods. In this approach, a near-instantaneous head change is applied in the well or DP rods and the subsequent head recovery is used to estimate ''K''. The initial head change can be introduced using a solid object (slug), compressed gas (pneumatic system), or by the addition/removal of a certain amount of water<ref name="Butler1997" />. The recovery data can be analyzed using different models of the well-formation configuration<ref name="Duffield2007" /><ref name="Butler1997" />. In contrast to pumping tests, slug tests provide a ''K'' estimate that is primarily a function of the materials in the vicinity of the screened interval of the test well. Slug test data analyses can be facilitated using software packages like AQTESOLV (Aquifer Test Solver), which has been developed for analyzing different test methods<ref name="Duffield2007" />.
  
 
Although a slug test is very simple in principle, considerable care must be taken in all stages of a test<ref name="Butler2005" /><ref name="Butler1997" />. Because test responses are highly sensitive to the materials immediately adjacent to the test well, that well must be appropriately developed before a test is performed. Otherwise, the ''K'' estimate can be biased by a low-''K'' zone (skin) that can form during well construction. In high-''K'' formations, the head recovery is rapid, and pneumatic methods are often used to minimize the time associated with test initiation. In addition, due to the relatively high flow velocity, the impact of pipe hydraulics on pressure readings should be considered when analyzing slug test data from high-''K'' zones<ref name="Butler2003">Butler Jr, J.J., Garnett, E.J. and Healey, J.M., 2003. Analysis of slug tests in formations of high hydraulic conductivity. Groundwater, 41(5), pp.620-631. [https://doi.org/10.1111/j.1745-6584.2003.tb02400.x doi: 10.1111/j.1745-6584.2003.tb02400.x]</ref>.  In low-''K'' formations, slug tests can take an extremely long time to complete, although test time can be significantly reduced by decreasing the effective casing radius (portion of well in which the water level is changing)<ref name="Butler1997" />.
 
Although a slug test is very simple in principle, considerable care must be taken in all stages of a test<ref name="Butler2005" /><ref name="Butler1997" />. Because test responses are highly sensitive to the materials immediately adjacent to the test well, that well must be appropriately developed before a test is performed. Otherwise, the ''K'' estimate can be biased by a low-''K'' zone (skin) that can form during well construction. In high-''K'' formations, the head recovery is rapid, and pneumatic methods are often used to minimize the time associated with test initiation. In addition, due to the relatively high flow velocity, the impact of pipe hydraulics on pressure readings should be considered when analyzing slug test data from high-''K'' zones<ref name="Butler2003">Butler Jr, J.J., Garnett, E.J. and Healey, J.M., 2003. Analysis of slug tests in formations of high hydraulic conductivity. Groundwater, 41(5), pp.620-631. [https://doi.org/10.1111/j.1745-6584.2003.tb02400.x doi: 10.1111/j.1745-6584.2003.tb02400.x]</ref>.  In low-''K'' formations, slug tests can take an extremely long time to complete, although test time can be significantly reduced by decreasing the effective casing radius (portion of well in which the water level is changing)<ref name="Butler1997" />.
 
   
 
   
Slug tests can be configured to obtain information about vertical variations in ''K'' along the screened (open) interval of a well. Test intervals can be isolated with straddle packers, and slug tests performed within that isolated interval. The straddle packers can be incrementally moved along the screened interval of a well to characterize the vertical variation of ''K'' at a relatively high resolution. Using a two-packer tool (Fig. 2A), slug tests have been performed in a number of 0.25-m intervals in a well at the Geohydrologic Experimental and Monitoring Site (GEMS) in the Kansas River valley<ref name="Butler2005" />.  At each isolated interval, multiple tests are performed, initiated with different head changes, following recommended test guidelines<ref name="Butler1997" />. An example data set is presented in Figure 2B. The multi-level slug test ''K'' estimates compare favorably with estimates obtained using other approaches (Fig. 3).  
+
Slug tests can be configured to obtain information about vertical variations in ''K'' along the screened (open) interval of a well. Test intervals can be isolated with straddle packers, and slug tests performed within that isolated interval. The straddle packers can be incrementally moved along the screened interval of a well to characterize the vertical variation of ''K'' at a relatively high resolution. Using a two-packer tool (Figure 2A), slug tests have been performed in a number of 0.25-m intervals in a well at the Geohydrologic Experimental and Monitoring Site (GEMS) in the Kansas River valley<ref name="Butler2005" />.  At each isolated interval, multiple tests are performed, initiated with different head changes, following recommended test guidelines<ref name="Butler1997" />. An example data set is presented in Figure 2B. The multi-level slug test ''K'' estimates compare favorably with estimates obtained using other approaches (Figure 3).  
  
<div><ul>
+
[[File:Liu1w2 Fig3.png|thumb|360 px|left|Figure 3. ''K'' Estimates from Different Field Methods at GEMS<ref name="Butler2005" />. Well DW is located 2 m east of GEMS4S.]]
<li style="display: inline-block;">[[File:Liu1w2 Fig2A.png|thumb|400 px|Figure 2A. Schematic of a Falling-Head Slug Test (water flows from well into aquifer) Performed in a Multi-Level Slug-Test System. Head change is introduced in standpipe, which is directly connected to the interval isolated by the straddle packers (not to scale).]] 
+
[[File:Liu1w2 Fig4.png|thumb|right|400 px|Figure 4. Schematic Illustrating Hydrostratigraphic Profiling, the combination of electrical conductivity (EC) and multilevel slug test profiling<ref name="Sellwood2005" />.]]  
<li style="display: inline-block;">[[File:Liu1w2 Fig2B.png|thumb|400 px|Figure 2B. Example Data plot from a Multilevel Slug Test at GEMS. Test is initiated by sudden depressurization of a pressurized air column<ref name="Butler2003" />.]]</li>
 
<li style="display: inline-block;">[[File:Liu1w2 Fig3.png|thumb|400 px|Figure 3. ''K'' Estimates from Different Field Methods at GEMS<ref name="Butler2005" />. Well DW is located 2 m east of GEMS4S.]]
 
</ul></div>
 
 
 
[[File:Liu1w2 Fig4.png|thumb|right|300 px|Figure 4. Schematic Illustrating Hydrostratigraphic Profiling, the combination of electrical conductivity (EC) and multilevel slug test profiling<ref name="Sellwood2005" />.]]  
 
  
 
In addition to existing wells, slug tests can also be performed in DP installations so that ''K'' estimates can be obtained virtually at any location in unconsolidated formations. Various approaches have been developed that allow slug tests to be performed at one or multiple levels in a single DP hole<ref>Hinsby, K., Bjerg, P.L., Andersen, L.J., Skov, B. and Clausen, E.V., 1992. A mini slug test method for determination of a local hydraulic conductivity of an unconfined sandy aquifer. Journal of Hydrology, 136(1-4), pp.87-106. [https://doi.org/10.1016/0022-1694(92)90006-H doi: 10.1016/0022-1694(92)90006-H]</ref><ref>Butler Jr, J.J., Healey, J.M., McCall, G.W., Garnett, E.J. and Loheide, S.P., 2002. Hydraulic tests with direct‐push equipment. Groundwater, 40(1), pp.25-36. [https://doi.org/10.1111/j.1745-6584.2002.tb02488.x doi: 10.1111/j.1745-6584.2002.tb02488.x]</ref><ref name="McCall2002">McCall, W., J.J. Butler, Jr., J.M. Healey, A.A. Lanier, S.M. Sellwood, and E.J. Garnett, 2002. A dual-tube direct-push method for vertical profiling of hydraulic conductivity in unconsolidated formations, Environ. & Eng. Geoscience, 8(2), 75-84.  [https://doi.org/10.2113/gseegeosci.8.2.75 doi: 10.2113/gseegeosci.8.2.75]</ref><ref name="Sellwood2005">Sellwood, S.M., Healey, J.M., Birk, S. and Butler, J.J., 2005. Direct‐push hydrostratigraphic profiling: coupling electrical logging and slug tests. Ground water, 43(1), pp.19-29. [http://dx.doi.org/10.1111/j.1745-6584.2005.tb02282.x doi: 10.1111/j.1745-6584.2005.tb02282.x]</ref>.  In McCall et al. (2002)<ref name="McCall2002" />, a pair of nested rod strings were driven to the test interval with a solid drive point attached to the end of inner rod string for advancement. Upon reaching the test depth, the drive point and inner rod string were retracted, and a screen was lowered to the bottom of the outer rod string. The outer rod string was pulled up while the screen was held in place, leaving the screen exposed to the surrounding formation. After the slug test was completed in the exposed screen, the screen was removed and the inner rod string with the attached solid drive point was reinserted. The nested rod strings were then driven to the next test depth. In low-''K'' formations such as silts and clays, the formation materials may not collapse completely back to the screen when the outer rod string is pulled up. In this case, the diameter of the borehole can be estimated and used in place of the screen diameter, or the problem can be avoided by using a coring tube of similar size to the screen to create a hole below the end of the outer rod string. Instead of setting the screen by pulling up the outer rod string, the screen can be directly inserted into the hole for slug testing. Regardless of how the screen is set into the formation, it is always recommended that the screen be appropriately developed before slug tests are performed. In low-''K'' formations, development may be limited to scraping the sides of the cored hole with a steel brush and removing the silty water from the screen with a low-flow pump.
 
In addition to existing wells, slug tests can also be performed in DP installations so that ''K'' estimates can be obtained virtually at any location in unconsolidated formations. Various approaches have been developed that allow slug tests to be performed at one or multiple levels in a single DP hole<ref>Hinsby, K., Bjerg, P.L., Andersen, L.J., Skov, B. and Clausen, E.V., 1992. A mini slug test method for determination of a local hydraulic conductivity of an unconfined sandy aquifer. Journal of Hydrology, 136(1-4), pp.87-106. [https://doi.org/10.1016/0022-1694(92)90006-H doi: 10.1016/0022-1694(92)90006-H]</ref><ref>Butler Jr, J.J., Healey, J.M., McCall, G.W., Garnett, E.J. and Loheide, S.P., 2002. Hydraulic tests with direct‐push equipment. Groundwater, 40(1), pp.25-36. [https://doi.org/10.1111/j.1745-6584.2002.tb02488.x doi: 10.1111/j.1745-6584.2002.tb02488.x]</ref><ref name="McCall2002">McCall, W., J.J. Butler, Jr., J.M. Healey, A.A. Lanier, S.M. Sellwood, and E.J. Garnett, 2002. A dual-tube direct-push method for vertical profiling of hydraulic conductivity in unconsolidated formations, Environ. & Eng. Geoscience, 8(2), 75-84.  [https://doi.org/10.2113/gseegeosci.8.2.75 doi: 10.2113/gseegeosci.8.2.75]</ref><ref name="Sellwood2005">Sellwood, S.M., Healey, J.M., Birk, S. and Butler, J.J., 2005. Direct‐push hydrostratigraphic profiling: coupling electrical logging and slug tests. Ground water, 43(1), pp.19-29. [http://dx.doi.org/10.1111/j.1745-6584.2005.tb02282.x doi: 10.1111/j.1745-6584.2005.tb02282.x]</ref>.  In McCall et al. (2002)<ref name="McCall2002" />, a pair of nested rod strings were driven to the test interval with a solid drive point attached to the end of inner rod string for advancement. Upon reaching the test depth, the drive point and inner rod string were retracted, and a screen was lowered to the bottom of the outer rod string. The outer rod string was pulled up while the screen was held in place, leaving the screen exposed to the surrounding formation. After the slug test was completed in the exposed screen, the screen was removed and the inner rod string with the attached solid drive point was reinserted. The nested rod strings were then driven to the next test depth. In low-''K'' formations such as silts and clays, the formation materials may not collapse completely back to the screen when the outer rod string is pulled up. In this case, the diameter of the borehole can be estimated and used in place of the screen diameter, or the problem can be avoided by using a coring tube of similar size to the screen to create a hole below the end of the outer rod string. Instead of setting the screen by pulling up the outer rod string, the screen can be directly inserted into the hole for slug testing. Regardless of how the screen is set into the formation, it is always recommended that the screen be appropriately developed before slug tests are performed. In low-''K'' formations, development may be limited to scraping the sides of the cored hole with a steel brush and removing the silty water from the screen with a low-flow pump.
  
Sellwood et al. (2005)<ref name="Sellwood2005" /> proposed a modification of the approach by McCall et al. (2002)<ref name="McCall2002" /> to reduce profiling time and gain more information about subsurface stratigraphy. Instead of performing slug tests on the way down, slug tests were performed at different depths as the outer rod string was pulled up (Fig. 4). This way, the number of changes between the outer and inner rod strings was minimized (i.e., only one change needed at the bottom of the profile). Furthermore, an electrical conductivity (EC) probe was attached to the inner rod string so that EC data could be collected as the rod strings were advanced; the EC data could then be used for selecting the intervals for slug tests as the rod strings were retracted. This approach, called hydrostratigraphic profiling, enabled the collection of information on electrical and hydraulic conductivity to be obtained at a speed and resolution that had previously not been possible.   
+
Sellwood et al. (2005)<ref name="Sellwood2005" /> proposed a modification of the approach by McCall et al. (2002)<ref name="McCall2002" /> to reduce profiling time and gain more information about subsurface stratigraphy. Instead of performing slug tests on the way down, slug tests were performed at different depths as the outer rod string was pulled up (Figure 4). This way, the number of changes between the outer and inner rod strings was minimized (i.e., only one change needed at the bottom of the profile). Furthermore, an electrical conductivity (EC) probe was attached to the inner rod string so that EC data could be collected as the rod strings were advanced; the EC data could then be used for selecting the intervals for slug tests as the rod strings were retracted. This approach, called hydrostratigraphic profiling, enabled information on electrical and hydraulic conductivity to be obtained at a speed and resolution that had previously not been possible.   
  
 
===Borehole Flowmeter Profiling===  
 
===Borehole Flowmeter Profiling===  
[[File:Liu1w2 Fig5.png|thumb|300 px|Figure 5. Schematic Diagram showing the borehole Flowmeter Profiling Procedure<ref name="Molz1989" />.]]  
+
[[File:Liu1w2 Fig5.png|thumb|350 px|left|Figure 5. Schematic Diagram showing the borehole Flowmeter Profiling Procedure<ref name="Molz1989" />.]]  
Borehole flowmeter profiling is one of the most efficient approaches for characterizing the vertical variations of ''K''<ref name="Molz1989">Molz, F.J., Morin, R.H., Hess, A.E., Melville, J.G. and Güven, O., 1989. The impeller meter for measuring aquifer permeability variations: evaluation and comparison with other tests. Water Resources Research, 25(7), pp.1677-1683. [https://doi.org/10.1029/WR025i007p01677 doi: 10.1029/WR025i007p01677]</ref>. This approach involves pumping a well at a constant rate while measuring the vertical flow rate within the screened interval (Fig. 5). A flowmeter is initially positioned at the bottom of the screen. After the pump is turned on and the pumping rate stabilizes, the flowmeter is gradually moved up, often in an incremental fashion with short stops, although continuous profiling is also done. After the entire screened interval is measured, the profile of vertical flow rates versus depths can be used for ''K'' estimation. Typically, two profiles are performed, one prior to pumping and one during pumping. The profile obtained during pumping is corrected with the profile prior to pumping to remove the impacts of ambient flow in the well. Due to the difficulty of detecting small differences in flow rates, the ''K'' estimates from flowmeter profiling are subject to more noise in less permeable settings<ref name="Bohling2012">Bohling, G.C., Liu, G., Knobbe, S.J., Reboulet, E.C., Hyndman, D.W., Dietrich, P. and Butler, J.J., 2012. Geostatistical analysis of centimeter‐scale hydraulic conductivity variations at the MADE site. Water Resources Research, 48(2). [https://doi.org/10.1029/2011wr010791 doi: 10.1029/2011WR010791]</ref>. In addition, the results of flowmeter profiling are sensitive to in-well hydraulics, such as the impedance of flowmeter to vertical flow and disturbance by the pump intake when it is too close to the measurement interval<ref name="Butler2005" />.
+
Borehole flowmeter profiling is one of the most efficient approaches for characterizing the vertical variations of ''K''<ref name="Molz1989">Molz, F.J., Morin, R.H., Hess, A.E., Melville, J.G. and Güven, O., 1989. The impeller meter for measuring aquifer permeability variations: evaluation and comparison with other tests. Water Resources Research, 25(7), pp.1677-1683. [https://doi.org/10.1029/WR025i007p01677 doi: 10.1029/WR025i007p01677]</ref>. This approach involves pumping a well at a constant rate while measuring the vertical flow rate within the screened interval (Figure 5). A flowmeter is initially positioned at the bottom of the screen. After the pump is turned on and the pumping rate stabilizes, the flowmeter is gradually moved up, often in an incremental fashion with short stops, although continuous profiling is also done. After the entire screened interval is measured, the profile of vertical flow rates versus depths can be used for ''K'' estimation. Typically, two profiles are performed, one prior to pumping and one during pumping. The profile obtained during pumping is corrected with the profile prior to pumping to remove the impacts of ambient flow in the well. Due to the difficulty of detecting small differences in flow rates, the ''K'' estimates from flowmeter profiling are subject to more noise in less permeable settings<ref name="Bohling2012">Bohling, G.C., Liu, G., Knobbe, S.J., Reboulet, E.C., Hyndman, D.W., Dietrich, P. and Butler, J.J., 2012. Geostatistical analysis of centimeter‐scale hydraulic conductivity variations at the MADE site. Water Resources Research, 48(2). [https://doi.org/10.1029/2011wr010791 doi: 10.1029/2011WR010791]</ref>. In addition, the results of flowmeter profiling are sensitive to in-well hydraulics, such as the impedance of flowmeter to vertical flow and disturbance by the pump intake when it is too close to the measurement interval<ref name="Butler2005" />.
  
 
Different types of flowmeters have been used for borehole profiling<ref>Young, S.C. and Pearson, H.S., 1995. The electromagnetic borehole flowmeter: Description and application. Groundwater Monitoring & Remediation, 15(4), pp.138-147. [https://doi.org/10.1111/j.1745-6592.1995.tb00561.x doi: 10.1111/j.1745-6592.1995.tb00561.x]</ref>.  For example, electromagnetic flowmeters are often used in environmental site investigations due to their lower flow rate detection limit and lower head loss (high head loss across the flowmeter can cause flow bypass through the filter pack of the well). However, they may be more expensive than other options and have certain practical limitations as well<ref>Bomana, G.K., Molz, F.J. and Boonec, K.D., 1997. Borehole flowmeter application in fluvial sediments: Methodology, results, and assessment. Groundwater, 35(3), pp.443-450. [https://doi.org/10.1111/j.1745-6584.1997.tb00104.x doi: 10.1111/j.1745-6584.1997.tb00104.x]</ref>.  Heat pulse flowmeters are also used for profiling in less-permeable settings.
 
Different types of flowmeters have been used for borehole profiling<ref>Young, S.C. and Pearson, H.S., 1995. The electromagnetic borehole flowmeter: Description and application. Groundwater Monitoring & Remediation, 15(4), pp.138-147. [https://doi.org/10.1111/j.1745-6592.1995.tb00561.x doi: 10.1111/j.1745-6592.1995.tb00561.x]</ref>.  For example, electromagnetic flowmeters are often used in environmental site investigations due to their lower flow rate detection limit and lower head loss (high head loss across the flowmeter can cause flow bypass through the filter pack of the well). However, they may be more expensive than other options and have certain practical limitations as well<ref>Bomana, G.K., Molz, F.J. and Boonec, K.D., 1997. Borehole flowmeter application in fluvial sediments: Methodology, results, and assessment. Groundwater, 35(3), pp.443-450. [https://doi.org/10.1111/j.1745-6584.1997.tb00104.x doi: 10.1111/j.1745-6584.1997.tb00104.x]</ref>.  Heat pulse flowmeters are also used for profiling in less-permeable settings.
  
 
The most common method of estimating ''K'' from borehole flowmeter profiles is based on the assumption that the lateral flow rate of an interval (zone between flow-rate measurements) is proportional to the ''K'' and thickness of that interval in a well fully screened in a perfectly layered aquifer,
 
The most common method of estimating ''K'' from borehole flowmeter profiles is based on the assumption that the lateral flow rate of an interval (zone between flow-rate measurements) is proportional to the ''K'' and thickness of that interval in a well fully screened in a perfectly layered aquifer,
 
[[File:Liu1w2 Eq4.png|center|400 px]]
 
  
 
{|
 
{|
 
|-
 
|-
|where:
+
| ||[[File:Liu1w2 Eq4.png|290 px]]
 +
|-
 +
|where:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 
|-
 
|-
 
|''K<sub>i</sub>''||is the ''K'' of test interval ''i'';
 
|''K<sub>i</sub>''||is the ''K'' of test interval ''i'';
Line 137: Line 135:
 
''K<sub>a</sub>'' can be obtained from pumping or slug tests at the same well, although it is generally recommended that slug tests be used, as the average ''K'' in the immediate vicinity of the well may differ from that determined from a pumping test in a laterally heterogeneous aquifer<ref name="Butler2005" />.  
 
''K<sub>a</sub>'' can be obtained from pumping or slug tests at the same well, although it is generally recommended that slug tests be used, as the average ''K'' in the immediate vicinity of the well may differ from that determined from a pumping test in a laterally heterogeneous aquifer<ref name="Butler2005" />.  
  
Figure 3 compares the results of borehole flowmeter profiling and multilevel slug tests at GEMS. Despite some local differences, the ''K'' estimates from borehole flowmeter profiling and multilevel slug tests are quite similar at that well.  
+
Figure 3 compares the results of borehole flowmeter profiling and multilevel slug tests at GEMS. Despite some local differences, the ''K'' estimates from borehole flowmeter profiling and multilevel slug tests are quite similar at that well.
  
 +
[[File:Liu1w2 Fig6A.png|thumb|200 px|right|Figure 6A. Schematic of the DPIL with a Screen Attached to the Probe Rod at its Lower End<ref name="Liu2012" />.]]
 
===DP Injection Logging===
 
===DP Injection Logging===
 
DP technology has shown great promise for characterization of ''K'' variations in shallow unconsolidated formations at the resolution, accuracy, and speed that are critically needed for practical investigations<ref>Dietrich, P. and Leven, C., 2009. Direct push-technologies. In Groundwater geophysics (pp. 347-366). Springer, Berlin, Heidelberg.</ref><ref name="McCall2005">McCall, W., Nielsen, D.M., Farrington, S.P. and Christy, T.M., 2005. Use of direct-push technologies in environmental site characterization and ground-water monitoring. In Practical Handbook of Environmental Site Characterization and Ground-Water Monitoring, Second Edition (pp. 355-482). CRC Press. [https://doi.org/10.1201/9781420032246.ch6 doi: 10.1201/9781420032246.ch6]</ref><ref>Leven, C., Weiß, H., Vienken, T. and Dietrich, P., 2011. Direct-Push-Technologien–Effiziente Untersuchungsmethoden für die Untergrunderkundung. Grundwasser, 16(4), pp.221-234. [https://doi.org/10.1007/s00767-011-0175-8 doi: 10.1007/s00767-011-0175-8]</ref>. Over the last few decades, a series of DP methods and probes have been developed for obtaining information about vertical variations in ''K''<ref name="Butler2005" /><ref name="Liu2012" /><ref name="Dietrich2008">Dietrich, P., Butler Jr, J.J. and Faiß, K., 2008. A rapid method for hydraulic profiling in unconsolidated formations. Groundwater, 46(2), pp.323-328. [https://doi.org/10.1111/j.1745-6584.2007.00377.x doi: 10.1111/j.1745-6584.2007.00377.x]</ref>.
 
DP technology has shown great promise for characterization of ''K'' variations in shallow unconsolidated formations at the resolution, accuracy, and speed that are critically needed for practical investigations<ref>Dietrich, P. and Leven, C., 2009. Direct push-technologies. In Groundwater geophysics (pp. 347-366). Springer, Berlin, Heidelberg.</ref><ref name="McCall2005">McCall, W., Nielsen, D.M., Farrington, S.P. and Christy, T.M., 2005. Use of direct-push technologies in environmental site characterization and ground-water monitoring. In Practical Handbook of Environmental Site Characterization and Ground-Water Monitoring, Second Edition (pp. 355-482). CRC Press. [https://doi.org/10.1201/9781420032246.ch6 doi: 10.1201/9781420032246.ch6]</ref><ref>Leven, C., Weiß, H., Vienken, T. and Dietrich, P., 2011. Direct-Push-Technologien–Effiziente Untersuchungsmethoden für die Untergrunderkundung. Grundwasser, 16(4), pp.221-234. [https://doi.org/10.1007/s00767-011-0175-8 doi: 10.1007/s00767-011-0175-8]</ref>. Over the last few decades, a series of DP methods and probes have been developed for obtaining information about vertical variations in ''K''<ref name="Butler2005" /><ref name="Liu2012" /><ref name="Dietrich2008">Dietrich, P., Butler Jr, J.J. and Faiß, K., 2008. A rapid method for hydraulic profiling in unconsolidated formations. Groundwater, 46(2), pp.323-328. [https://doi.org/10.1111/j.1745-6584.2007.00377.x doi: 10.1111/j.1745-6584.2007.00377.x]</ref>.
 
    
 
    
DP Injection Logging (DPIL) is one of the most powerful approaches for ''K'' characterization (in terms of resolution and speed of acquisition) and is also one of the most widely used techniques in environmental site investigations<ref name="Liu2012" /><ref>Maliva, R.G., 2016. Direct-push technology. In Aquifer Characterization Techniques (pp. 383-402). Springer, Cham. [https://doi.org/10.1007/978-3-319-32137-0_12 doi: 10.1007/978-3-319-32137-0_12]</ref>. This approach consists of advancing a probe with a single screened port (Fig. 6A). Water is injected continuously through the screen while the probe is advanced and the pressure response to injection is monitored behind the screen or at the surface. The profile of the ratio of injection rate to injection pressure is closely related to vertical variations in ''K''. However, factors other than the hydraulic conductivity can potentially affect the pressure response (e.g., injection line losses, probe advancement speed, formation alteration, and variations in specific storage), so DPIL ratios only provide relative information about ''K'' variations. Methods are needed for transforming DPIL ratios into actual ''K'' estimates; these methods typically involve correlations with nearby ''K'' data.
+
DP Injection Logging (DPIL) is one of the most powerful approaches for ''K'' characterization (in terms of resolution and speed of acquisition) and is also one of the most widely used techniques in environmental site investigations<ref name="Liu2012" /><ref>Maliva, R.G., 2016. Direct-push technology. In Aquifer Characterization Techniques (pp. 383-402). Springer, Cham. [https://doi.org/10.1007/978-3-319-32137-0_12 doi: 10.1007/978-3-319-32137-0_12]</ref>. This approach consists of advancing a probe with a single screened port (Figure 6A). Water is injected continuously through the screen while the probe is advanced and the pressure response to injection is monitored behind the screen or at the surface. The profile of the ratio of injection rate to injection pressure is closely related to vertical variations in ''K''. However, factors other than the hydraulic conductivity can potentially affect the pressure response (e.g., injection line losses, probe advancement speed, formation alteration, and variations in specific storage), so DPIL ratios only provide relative information about ''K'' variations. Methods are needed for transforming DPIL ratios into actual ''K'' estimates; these methods typically involve correlations with nearby ''K'' data.
+
[[File:Liu1w2 Fig6B.png|thumb|right|200 px|Figure 6B. Artistic Rendering of the HPT (continuous DPIL probe combined with an EC Wenner array) <ref name="Lessoff2010" />.]]
 +
 
 
There are currently two variants of DPIL: discontinuous<ref name="Dietrich2008" /><ref name="Lessoff2010">Lessoff, Steven C., Uwe Schneidewind, Carsten Leven, Philipp Blum, Peter Dietrich, and Gedeon Dagan. "Spatial characterization of the hydraulic conductivity using direct‐push injection logging." Water Resources Research 46, no. 12 (2010). [https://doi.org/10.1029/2009WR008949 doi: 10.1029/2009WR008949]</ref> and continuous<ref name="McCall2005" /><ref>McCall, W., Christy, T.M., Christopherson, T. and Issacs, H., 2009. Application of direct push methods to investigate uranium distribution in an alluvial aquifer. Groundwater Monitoring & Remediation, 29(4), pp.65-76. [https://doi.org/10.1111/j.1745-6592.2009.01258.x doi: 10.1111/j.1745-6592.2009.01258.x]</ref><ref name="Liu2009">Liu, G., Butler, J.J., Bohling, G.C., Reboulet, E., Knobbe, S. and Hyndman, D.W., 2009. A new method for high‐resolution characterization of hydraulic conductivity. Water Resources Research, 45(8). [https://doi.org/10.1029/2009WR008319 doi: 10.1029/2009WR008319]</ref> modes. In discontinuous DPIL mode, probe advancement is briefly halted at the desired depth and the injection rate is then varied in a step-wise fashion while injection rates and pressures are measured. In continuous DPIL mode, injection rate and pressure are measured as the probe is continuously advanced. Compared to continuous DPIL, discontinuous DPIL has a few advantages<ref name="Liu2012" />:  
 
There are currently two variants of DPIL: discontinuous<ref name="Dietrich2008" /><ref name="Lessoff2010">Lessoff, Steven C., Uwe Schneidewind, Carsten Leven, Philipp Blum, Peter Dietrich, and Gedeon Dagan. "Spatial characterization of the hydraulic conductivity using direct‐push injection logging." Water Resources Research 46, no. 12 (2010). [https://doi.org/10.1029/2009WR008949 doi: 10.1029/2009WR008949]</ref> and continuous<ref name="McCall2005" /><ref>McCall, W., Christy, T.M., Christopherson, T. and Issacs, H., 2009. Application of direct push methods to investigate uranium distribution in an alluvial aquifer. Groundwater Monitoring & Remediation, 29(4), pp.65-76. [https://doi.org/10.1111/j.1745-6592.2009.01258.x doi: 10.1111/j.1745-6592.2009.01258.x]</ref><ref name="Liu2009">Liu, G., Butler, J.J., Bohling, G.C., Reboulet, E., Knobbe, S. and Hyndman, D.W., 2009. A new method for high‐resolution characterization of hydraulic conductivity. Water Resources Research, 45(8). [https://doi.org/10.1029/2009WR008319 doi: 10.1029/2009WR008319]</ref> modes. In discontinuous DPIL mode, probe advancement is briefly halted at the desired depth and the injection rate is then varied in a step-wise fashion while injection rates and pressures are measured. In continuous DPIL mode, injection rate and pressure are measured as the probe is continuously advanced. Compared to continuous DPIL, discontinuous DPIL has a few advantages<ref name="Liu2012" />:  
  
Line 150: Line 150:
 
*use of multiple injection rates at each depth allows a better assessment of the formation response to injection.
 
*use of multiple injection rates at each depth allows a better assessment of the formation response to injection.
  
However, due to the need to halt advancement, discontinuous DPIL requires significantly more time than continuous DPIL, and, as a result, resolution below 10 - 20 cm is rarely possible. In common practice, continuous DPIL is generally preferred due to its speed (a 20-m profile requires about 2 hours) and resolution (one measurement per 1.5 cm vertical interval). A continuous DPIL probe, supplemented with an EC sensor, is commercially available from Geoprobe Systems under the tradename of the Hydraulic Profiling Tool (HPT; Fig. 6B).  
+
However, due to the need to halt advancement, discontinuous DPIL requires significantly more time than continuous DPIL, and, as a result, resolution below 10 - 20 cm is rarely possible. In common practice, continuous DPIL is generally preferred due to its speed (a 20-m profile requires about 2 hours) and resolution (one measurement per 1.5 cm vertical interval). A continuous DPIL probe, supplemented with an EC sensor, is commercially available from Geoprobe Systems under the tradename of the Hydraulic Profiling Tool (HPT; Figure 6B).  
  
 
The HPT is most effective in moderately permeable formations with a ''K'' range of 0.03 to 10 m/d, although the detection range can be improved by modifying equipment and operating procedures<ref name="Liu2012" />. In standard HPT practice, the probe advancement rate is 2 cm/s and the continuous water injection rate is about 300 mL/min. In high ''K'' formations (e.g., ''K'' > 10 m/d), the injection pressure response to the standard injection rate is generally too small to be reliably measured. A larger injection rate can be used to increase the pressure signal. However, the line loss between the transducer and injection screen, which increases quadratically with the injection rate, may become significant and require additional step tests to remove its impact<ref name="Liu2012" />. On the other hand, in low ''K'' formations (e.g., ''K'' < 0.03 m/d), the pressure increase from probe advancement and water injection may become so large that it will cause formation alterations and exceed the upper measurement limit of the transducer. A recent study has demonstrated that by reducing both probe advancement speed and flow injection rate, HPT can be applied to formations with ''K'' as low as 0.003 m/d<ref name="Liu2018">Liu, G., Borden, R.C. and Butler Jr, J.J., 2018. Simulation Assessment of Direct Push Injection Logging for High‐Resolution Aquifer Characterization. Groundwater. [https://doi.org/10.1111/gwat.12826 doi: 10.1111/gwat.12826]</ref>.  
 
The HPT is most effective in moderately permeable formations with a ''K'' range of 0.03 to 10 m/d, although the detection range can be improved by modifying equipment and operating procedures<ref name="Liu2012" />. In standard HPT practice, the probe advancement rate is 2 cm/s and the continuous water injection rate is about 300 mL/min. In high ''K'' formations (e.g., ''K'' > 10 m/d), the injection pressure response to the standard injection rate is generally too small to be reliably measured. A larger injection rate can be used to increase the pressure signal. However, the line loss between the transducer and injection screen, which increases quadratically with the injection rate, may become significant and require additional step tests to remove its impact<ref name="Liu2012" />. On the other hand, in low ''K'' formations (e.g., ''K'' < 0.03 m/d), the pressure increase from probe advancement and water injection may become so large that it will cause formation alterations and exceed the upper measurement limit of the transducer. A recent study has demonstrated that by reducing both probe advancement speed and flow injection rate, HPT can be applied to formations with ''K'' as low as 0.003 m/d<ref name="Liu2018">Liu, G., Borden, R.C. and Butler Jr, J.J., 2018. Simulation Assessment of Direct Push Injection Logging for High‐Resolution Aquifer Characterization. Groundwater. [https://doi.org/10.1111/gwat.12826 doi: 10.1111/gwat.12826]</ref>.  
  
 
Transforming DPIL profiles into ''K'' estimates is based on empirical relationships<ref name="McCall2005" /><ref name="Dietrich2008" /><ref name="Liu2009" /> <ref name="McCall2010">McCall, W. and Christy, T.M., 2010. Development of a hydraulic conductivity estimate for the Hydraulic Profiling Tool (HPT). In The 2010 North American Environmental Field Conference & Exposition: Conference Program with Abstracts: Session VII.</ref>. Although a general relationship has been presented<ref name="McCall2010" />, site-specific calibrations are typically used to improve the reliability of ''K'' estimates. The site-specific calibrations can be performed by comparing the DPIL profiles to nearby ''K'' estimates obtained via other means<ref name="Dietrich2008" /><ref name="Lessoff2010" /> or by a modeling approach that directly combines the DPIL data with collocated hydraulic tests such as seen in Figure 6C, which shows a series of ''K'' profiles from continuous DPIL as calibrated by collocated hydraulic tests at an alluvial aquifer in Mississippi<ref name="Liu2009" />.
 
Transforming DPIL profiles into ''K'' estimates is based on empirical relationships<ref name="McCall2005" /><ref name="Dietrich2008" /><ref name="Liu2009" /> <ref name="McCall2010">McCall, W. and Christy, T.M., 2010. Development of a hydraulic conductivity estimate for the Hydraulic Profiling Tool (HPT). In The 2010 North American Environmental Field Conference & Exposition: Conference Program with Abstracts: Session VII.</ref>. Although a general relationship has been presented<ref name="McCall2010" />, site-specific calibrations are typically used to improve the reliability of ''K'' estimates. The site-specific calibrations can be performed by comparing the DPIL profiles to nearby ''K'' estimates obtained via other means<ref name="Dietrich2008" /><ref name="Lessoff2010" /> or by a modeling approach that directly combines the DPIL data with collocated hydraulic tests such as seen in Figure 6C, which shows a series of ''K'' profiles from continuous DPIL as calibrated by collocated hydraulic tests at an alluvial aquifer in Mississippi<ref name="Liu2009" />.
 
+
[[File:Liu1w2 Fig6C.png|thumb|550 px|left|Figure 6C. Six Example DPIL (red curve) ''K'' Profiles Calibrated by Collocated Direct Push Permeameter Tests (blue bars) at the MacroDispersion Experiment (MADE) Site, Columbus, MS. The dashed curves represent the upper ''K'' limit for continuous DPIL under the given tool operating procedure<ref name="Liu2009" />.]]
<div><ul>
 
<li style="display: inline-block;">[[File:Liu1w2 Fig6A.png|thumb|Figure 6A. Schematic of the DPIL with a Screen Attached to the Probe Rod at its Lower End<ref name="Liu2012" />.]]</li>
 
<li style="display: inline-block; middle:">[[File:Liu1w2 Fig6B.png|thumb|400 px|Figure 6B. Artistic Rendering of the HPT (continuous DPIL probe combined with an EC Wenner array) <ref name="Lessoff2010" />.]]</li>
 
<li style="display: inline-block;">[[File:Liu1w2 Fig6C.png|thumb|600 px|Figure 6C. Six Example DPIL (red curve) ''K'' Profiles Calibrated by Collocated Direct Push Permeameter Tests (blue bars) at the MacroDispersion Experiment (MADE) Site, Columbus, MS. The dashed curves represent the upper ''K'' limit for continuous DPIL under the given tool operating procedure<ref name="Liu2009" />.]]</li>
 
</ul></div>
 
  
 
There are two major assumptions invoked for estimating ''K'' from DPIL profiles. First, near-rod material compaction from probe advancement has a similar impact on ''K'' at all depths of the profile. Second, the impact of formation specific storage is either similar across the profile or can be ignored. Both assumptions appear to be valid for sands and gravels. In silts and clays, particularly those with high compressibility, considerable caution must be used to assess the viability of these assumptions<ref name="Liu2012" /><ref name="Liu2018" />.
 
There are two major assumptions invoked for estimating ''K'' from DPIL profiles. First, near-rod material compaction from probe advancement has a similar impact on ''K'' at all depths of the profile. Second, the impact of formation specific storage is either similar across the profile or can be ignored. Both assumptions appear to be valid for sands and gravels. In silts and clays, particularly those with high compressibility, considerable caution must be used to assess the viability of these assumptions<ref name="Liu2012" /><ref name="Liu2018" />.
Line 168: Line 163:
 
[[File:Liu1w2 Fig7B.png|thumb|left|300 px|Figure 7B. Example DPP Test Sequence at GEMS<ref name="Liu2012" />.]]
 
[[File:Liu1w2 Fig7B.png|thumb|left|300 px|Figure 7B. Example DPP Test Sequence at GEMS<ref name="Liu2012" />.]]
 
[[File:Liu1w2 Fig8.png|thumb|right|Figure 8. Comparison of DPP and DP Multilevel Slug Test ''K'' profiles at GEMS. The slug test profile is located 2 m northwest of DPP profile 1, and 2 m southeast of DPP profile 2<ref name="Butler2005" />.]]
 
[[File:Liu1w2 Fig8.png|thumb|right|Figure 8. Comparison of DPP and DP Multilevel Slug Test ''K'' profiles at GEMS. The slug test profile is located 2 m northwest of DPP profile 1, and 2 m southeast of DPP profile 2<ref name="Butler2005" />.]]
The Direct Push Permeameter (DPP) consists of a short cylindrical screen with two pressure transducers inset into the probe at short distances above the screen (Fig. 7A)<ref name="Butler2005" /><ref name="Liu2012" /><ref name="Butler2007">Butler Jr, J.J., Dietrich, P., Wittig, V. and Christy, T., 2007. Characterizing hydraulic conductivity with the direct‐push permeameter. Groundwater, 45(4), pp.409-419. [https://doi.org/10.1111/j.1745-6584.2007.00300.x doi: 10.1111/j.1745-6584.2007.00300.x</ref>. The DPP probe is advanced to a depth at which a ''K'' estimate is needed and then a series of short-term injection tests are performed. ''K'' is then typically estimated from the spherical form of Darcy’s Law using the injection rate and the injection-induced pressure responses at the two transducers. The resulting estimate is a weighted average over the interval (approx. 0.4 m in Fig. 7A) between the screen and the farthest transducer; material outside of that interval has little influence <ref name="Liu2008">Liu, G., Bohling, G.C. and Butler, J.J., 2008. Simulation assessment of the direct‐push permeameter for characterizing vertical variations in hydraulic conductivity. Water Resources Research, 44(2). [https://doi.org/10.1029/2007wr006078 doi:10.1029/2007WR006078]</ref>. Horizontally, the DPP shows little sensitivity to conditions greater than ∼ 0.5 m away from the tool. The DPP test procedure only requires that steady-shape (constant hydraulic gradient) flow conditions be attained. True steady-state conditions are not required, which can lead to a significant reduction in the measurement time. A DPP test sequence generally requires 10 – 15 minutes per measurement interval in moderate to high-''K'' formations<ref name="Butler2007" />.  
+
The Direct Push Permeameter (DPP) consists of a short cylindrical screen with two pressure transducers inset into the probe at short distances above the screen (Figure 7A)<ref name="Butler2005" /><ref name="Liu2012" /><ref name="Butler2007">Butler Jr, J.J., Dietrich, P., Wittig, V. and Christy, T., 2007. Characterizing hydraulic conductivity with the direct‐push permeameter. Groundwater, 45(4), pp.409-419. [https://doi.org/10.1111/j.1745-6584.2007.00300.x doi: 10.1111/j.1745-6584.2007.00300.x</ref>. The DPP probe is advanced to a depth at which a ''K'' estimate is needed and then a series of short-term injection tests are performed. ''K'' is then typically estimated from the spherical form of Darcy’s Law using the injection rate and the injection-induced pressure responses at the two transducers. The resulting estimate is a weighted average over the interval between the screen and the farthest transducer (approx. 0.4 m in Figure 7A); material outside of that interval has little influence <ref name="Liu2008">Liu, G., Bohling, G.C. and Butler, J.J., 2008. Simulation assessment of the direct‐push permeameter for characterizing vertical variations in hydraulic conductivity. Water Resources Research, 44(2). [https://doi.org/10.1029/2007wr006078 doi:10.1029/2007WR006078]</ref>. Horizontally, the DPP shows little sensitivity to conditions greater than ∼ 0.5 m away from the tool. The DPP test procedure only requires that steady-shape (constant hydraulic gradient) flow conditions be attained. True steady-state conditions are not required, which can lead to a significant reduction in the measurement time. A DPP test sequence generally requires 10 – 15 minutes per measurement interval in moderate to high-''K'' formations<ref name="Butler2007" />.  
  
 
Figure 7B shows an example DPP test sequence at GEMS<ref name="Liu2012" />. Three tests with different rates are used to assess the reliability of the ''K'' estimates. The first injection rate is typically based on the expected ''K'' value in the test interval (e.g., 800 mL/min for sands). Depending on the magnitude of the pressure response in the first test, the second injection rate can be set higher or lower than the first (if pressure response is high in the first test, the flow rate in the second test is reduced, and vice versa). The third flow rate is set similar to that in the first test to check whether a similar pressure response can be obtained. The ''K'' values computed using Darcy’s Law were 66.0, 63.1 and 65.4 m/d, respectively, for the three tests in Figure 7B. A set of consistent ''K'' estimates is the goal; a lack of consistency could indicate DPP system instability or formation alteration. Figure 8 shows the comparison of DPP ''K'' estimates with those obtained from DP multilevel slug tests at GEMS<ref name="Butler2007" />. The results from DPP profiling have proven extremely useful for transforming collocated DPIL profiles into ''K'' estimates<ref name="Bohling2012" /><ref name="Liu2009" />.
 
Figure 7B shows an example DPP test sequence at GEMS<ref name="Liu2012" />. Three tests with different rates are used to assess the reliability of the ''K'' estimates. The first injection rate is typically based on the expected ''K'' value in the test interval (e.g., 800 mL/min for sands). Depending on the magnitude of the pressure response in the first test, the second injection rate can be set higher or lower than the first (if pressure response is high in the first test, the flow rate in the second test is reduced, and vice versa). The third flow rate is set similar to that in the first test to check whether a similar pressure response can be obtained. The ''K'' values computed using Darcy’s Law were 66.0, 63.1 and 65.4 m/d, respectively, for the three tests in Figure 7B. A set of consistent ''K'' estimates is the goal; a lack of consistency could indicate DPP system instability or formation alteration. Figure 8 shows the comparison of DPP ''K'' estimates with those obtained from DP multilevel slug tests at GEMS<ref name="Butler2007" />. The results from DPP profiling have proven extremely useful for transforming collocated DPIL profiles into ''K'' estimates<ref name="Bohling2012" /><ref name="Liu2009" />.
Line 181: Line 176:
 
==Indirect Approaches==
 
==Indirect Approaches==
 
===Geophysical Methods===
 
===Geophysical Methods===
Estimates for formation ''K'' can also be inferred from the measurement of more readily evaluated geophysical properties using empirical, often site-specific, relationships. The geophysical properties that may be used for ''K'' estimation include, but are not limited to, electrical conductivity (EC) (or the reverse, electrical resistivity), the hydrogen response to nuclear magnetic resonance (NMR), and natural gamma radiation.  In addition to ''K'' estimation, geophysical methods have also been used for many purposes in environmental site investigation, assessment and monitoring (see related article –[http://www.enviro.wiki/index.php?title=Geophysical_Methods Geophysical Methods].
+
Estimates for formation ''K'' can also be inferred from the measurement of more readily evaluated geophysical properties using empirical, often site-specific, relationships. The geophysical properties that may be used for ''K'' estimation include, but are not limited to, electrical conductivity (EC) (or the reverse, electrical resistivity), the hydrogen response to nuclear magnetic resonance (NMR), and natural gamma radiation.  In addition to ''K'' estimation, geophysical methods have also been used for many purposes in environmental site investigation, assessment and monitoring (see related article –[http://www.enviro.wiki/index.php?title=Geophysical_Methods Geophysical Methods]).
  
Borehole EC logging measures vertical variations in the bulk EC of the surrounding formation, which can then be used to estimate ''K'' based on empirical relationships<ref>Purvance, D.T. and Andricevic, R., 2000. On the electrical‐hydraulic conductivity correlation in aquifers. Water Resources Research, 36(10), pp.2905-2913. [https://doi.org/10.1029/2000wr900165 doi: 10.1029/2000wr900165]</ref><ref>Slater, L. and Lesmes, D.P., 2002. Electrical‐hydraulic relationships observed for unconsolidated sediments. Water Resources Research, 38(10), pp.31-1. [https://doi.org/10.1029/2001wr001075 doi:10.1029/2001WR001075]</ref>. EC has been adapted to DP equipment for more efficient assessment of subsurface conditions than the approaches that rely on existing wells<ref>Christy, C.D., Christy, T.M. and Wittig, V., 1994, May. A percussion probing tool for the direct sensing of soil conductivity. In Proceedings of the 8th National Outdoor Action Conference (Vol. 38, pp. 1-394). Dublin, Ohio: National Ground Water Association.</ref>. Variations in the EC of saturated media are primarily a function of pore-fluid chemistry, clay content, and total porosity. If variations in pore-fluid chemistry can be neglected, EC profiling can be used to delineate the major hydrostratigraphic features, and thus large-scale variations in ''K''<ref name="Schulmeister2003">Schulmeister, M.K., Butler, J.J., Healey, J.M., Zheng, L., Wysocki, D.A. and McCall, G.W., 2003. Direct‐push electrical conductivity logging for high‐resolution hydrostratigraphic characterization. Groundwater Monitoring & Remediation, 23(3), pp.52-62. [http://dx.doi.org/10.1111/j.1745-6592.2003.tb00683.x doi: 10.1111/j.1745-6592.2003.tb00683.x]</ref>. However, when electrically-conductive clay is absent, EC profiling generally provides little information about vertical variations in ''K'' in silts, sands and gravels.  
+
Borehole EC logging measures vertical variations in the bulk EC of the surrounding formation, which can then be used to estimate ''K'' based on empirical relationships<ref>Purvance, D.T. and Andricevic, R., 2000. On the electrical‐hydraulic conductivity correlation in aquifers. Water Resources Research, 36(10), pp.2905-2913. [https://doi.org/10.1029/2000wr900165 doi: 10.1029/2000wr900165]</ref><ref>Slater, L. and Lesmes, D.P., 2002. Electrical‐hydraulic relationships observed for unconsolidated sediments. Water Resources Research, 38(10), pp.31-1. [https://doi.org/10.1029/2001wr001075 doi:10.1029/2001WR001075]</ref>. EC has been adapted to DP equipment for more flexible assessment of subsurface conditions than the approaches that rely on existing wells<ref>Christy, C.D., Christy, T.M. and Wittig, V., 1994, May. A percussion probing tool for the direct sensing of soil conductivity. In Proceedings of the 8th National Outdoor Action Conference (Vol. 38, pp. 1-394). Dublin, Ohio: National Ground Water Association.</ref>. Variations in the EC of saturated media are primarily a function of pore-fluid chemistry, clay content, and total porosity. If variations in pore-fluid chemistry can be neglected, EC profiling can be used to delineate the major hydrostratigraphic features, and thus large-scale variations in ''K''<ref name="Schulmeister2003">Schulmeister, M.K., Butler, J.J., Healey, J.M., Zheng, L., Wysocki, D.A. and McCall, G.W., 2003. Direct‐push electrical conductivity logging for high‐resolution hydrostratigraphic characterization. Groundwater Monitoring & Remediation, 23(3), pp.52-62. [http://dx.doi.org/10.1111/j.1745-6592.2003.tb00683.x doi: 10.1111/j.1745-6592.2003.tb00683.x]</ref>. However, when electrically-conductive clay is absent, EC profiling generally provides little information about vertical variations in ''K'' in silts, sands and gravels.  
  
 
Borehole natural gamma logging provides a record of natural gamma radiation versus depth. This radiation is quantified by counting the gamma particles passing through a scintillation crystal in a certain time interval. A high natural gamma reading is generally associated with clay-rich intervals, while a low reading is generally associated with sands and gravels<ref>Keys, W.S., 1990. Borehole geophysics applied to ground-water investigations (No. 02-E2). USGPO; For sale by the Books and Open-File Reports Section, US Geological Survey. [https://doi.org/10.3133/twri02E2 doi: 10.3133/twri02E2]</ref>. Like EC logging, natural gamma logging provides an effective measurement of formation characteristics when clay is present, but is of less use in silts, sands and gravels.
 
Borehole natural gamma logging provides a record of natural gamma radiation versus depth. This radiation is quantified by counting the gamma particles passing through a scintillation crystal in a certain time interval. A high natural gamma reading is generally associated with clay-rich intervals, while a low reading is generally associated with sands and gravels<ref>Keys, W.S., 1990. Borehole geophysics applied to ground-water investigations (No. 02-E2). USGPO; For sale by the Books and Open-File Reports Section, US Geological Survey. [https://doi.org/10.3133/twri02E2 doi: 10.3133/twri02E2]</ref>. Like EC logging, natural gamma logging provides an effective measurement of formation characteristics when clay is present, but is of less use in silts, sands and gravels.
  
Efforts have been made to adapt nuclear magnetic resonance (NMR) profiling, a widely used borehole technique in the petroleum industry, to environmental site investigations<ref name="WalshD2013">Walsh, D., Turner, P., Grunewald, E., Zhang, H., Butler Jr, J.J., Reboulet, E., Knobbe, S., Christy, T., Lane Jr, J.W., Johnson, C.D. and Munday, T., 2013. A small‐diameter NMR logging tool for groundwater investigations. Groundwater, 51(6), pp.914-926. [https://doi.org/10.1111/gwat.12024 doi: 10.1111/gwat.12024]</ref><ref name="Walsh2011"> Walsh, D., Turner, P., Grunewald, E., Butler, J., Knight, R., Reboulet, E., Knobbe, S., Christy, T. and McCall, W., 2011, January. Field demonstration of nuclear magnetic resonance (NMR) logging tools for groundwater and environmental investigations. In Symposium on the Application of Geophysics to Engineering and Environmental Problems 2011 (pp. 69-69). Society of Exploration Geophysicists. [https://doi.org/10.4133/1.3614224 doi: 10.4133/1.3614224]</ref> <ref name="Knight2016">Knight, R., Walsh, D.O., Butler Jr, J.J., Grunewald, E., Liu, G., Parsekian, A.D., Reboulet, E.C., Knobbe, S. and Barrows, M., 2016. NMR logging to estimate hydraulic conductivity in unconsolidated aquifers. Groundwater, 54(1), pp.104-114. [https://doi.org/10.1111/gwat.12324 doi: 10.1111/gwat.12324]</ref>. This approach involves measuring the response of hydrogen atoms (protons) to a series of imposed magnetic fields at radio frequency. The response is a function of, among other things, the pore-size distribution of material in a thin cylindrical ring centered on the probe (Figs. 9A and 9B). ''K'' is then estimated from the pore-size distribution information using different empirical relationships, all of which have their origins in the Kozeny-Carman equation<ref name="Knight2016" />. Figure 9C compares the NMR results with DPP K estimates at GEMS.  
+
[[File:Liu1w2 Fig9A.png|thumb|200 px|left|Figure 9A. Measurement Domain around a NMR Borehole Probe]]
 +
[[File:Liu1w2 Fig9B.png|thumb|300 px|Figure 9B. NMR Responses as a Function of Pore Size, Where ''A'' is the Initial Magnetization, T<sub>2</sub> is NMR Relaxation Time<ref name="Walsh2011" />.]]
 +
Efforts have been made to adapt [[wikipedia: Nuclear magnetic resonance logging | nuclear magnetic resonance (NMR) profiling]], a widely used borehole technique in the petroleum industry, to environmental site investigations<ref name="WalshD2013">Walsh, D., Turner, P., Grunewald, E., Zhang, H., Butler Jr, J.J., Reboulet, E., Knobbe, S., Christy, T., Lane Jr, J.W., Johnson, C.D. and Munday, T., 2013. A small‐diameter NMR logging tool for groundwater investigations. Groundwater, 51(6), pp.914-926. [https://doi.org/10.1111/gwat.12024 doi: 10.1111/gwat.12024]</ref><ref name="Walsh2011"> Walsh, D., Turner, P., Grunewald, E., Butler, J., Knight, R., Reboulet, E., Knobbe, S., Christy, T. and McCall, W., 2011, January. Field demonstration of nuclear magnetic resonance (NMR) logging tools for groundwater and environmental investigations. In Symposium on the Application of Geophysics to Engineering and Environmental Problems 2011 (pp. 69-69). Society of Exploration Geophysicists. [https://doi.org/10.4133/1.3614224 doi: 10.4133/1.3614224]</ref> <ref name="Knight2016">Knight, R., Walsh, D.O., Butler Jr, J.J., Grunewald, E., Liu, G., Parsekian, A.D., Reboulet, E.C., Knobbe, S. and Barrows, M., 2016. NMR logging to estimate hydraulic conductivity in unconsolidated aquifers. Groundwater, 54(1), pp.104-114. [https://doi.org/10.1111/gwat.12324 doi: 10.1111/gwat.12324]</ref>. This approach involves measuring the response of hydrogen atoms (protons) to a series of imposed magnetic fields at radio frequency. The response is a function of, among other things, the pore-size distribution of material in a thin cylindrical ring centered on the probe (Figures 9A and 9B). ''K'' is then estimated from the pore-size distribution information using different empirical relationships, all of which have their origins in the [[wikipedia: Kozeny-Carman equation | Kozeny-Carman equation]]<ref name="Knight2016" />. Figure 9C compares the NMR results with DPP ''K'' estimates at GEMS.  
  
<div><ul>
+
[[File:Liu1w2 Fig9C.png|thumb|450 px|left|Figure 9C. Comparison of NMR and DPP ''K'' estimates at two GEMS locations. The dash lines show the range of the NMR ''K'' estimates due to the uncertainty in the empirical parameters used in the pore size - ''K'' relationship<ref name="Knight2016" />.]]
<li style="display: inline-block;">[[File:Liu1w2 Fig9A.png|thumb|left|Figure 9A. Measurement Domain around a NMR Borehole Probe]]</li>
 
<li style="display: inline-block;">[[File:Liu1w2 Fig9B.png|thumb|Figure 9B. NMR Responses as a Function of Pore Size, Where ''A'' is the Initial Magnetization, T<sub>2</sub> is NMR Relaxation Time<ref name="Walsh2011" />.]]</li>
 
<li style="display: inline-block;">[[File:Liu1w2 Fig9C.png|thumb|Figure 9C. Comparison of NMR and DPP ''K'' estimates at two GEMS locations. The dash lines show the range of the NMR ''K'' estimates due to the uncertainty in the empirical parameters used in the pore size - ''K'' relationship<ref name="Knight2016" />.]]</li>
 
</ul></div>
 
  
 
The NMR logging tools used in the petroleum industry are typically too large for use by the hydrological community (e.g., tools are typically well over 10 ft long and have a diameter larger than 5 inches). The equipment costs are also prohibitively high for hydrological users. After recent technology advancements, NMR measurement can be made much more effectively both at the surface and in a borehole, including with DP equipment<ref name="WalshD2013" /><ref name="Walsh2011" /><ref name="Knight2016" />. The equipment costs are significantly decreased compared to those used in the petroleum industry. The measurement time is reduced down to a few minutes per interval. The vertical sampling interval of most current logging tools is about 0.5 m. Different tool diameters are available, with the smallest being deployable in a 2-inch well. Due to these improvements, NMR has been increasingly used by environmental site investigators worldwide.
 
The NMR logging tools used in the petroleum industry are typically too large for use by the hydrological community (e.g., tools are typically well over 10 ft long and have a diameter larger than 5 inches). The equipment costs are also prohibitively high for hydrological users. After recent technology advancements, NMR measurement can be made much more effectively both at the surface and in a borehole, including with DP equipment<ref name="WalshD2013" /><ref name="Walsh2011" /><ref name="Knight2016" />. The equipment costs are significantly decreased compared to those used in the petroleum industry. The measurement time is reduced down to a few minutes per interval. The vertical sampling interval of most current logging tools is about 0.5 m. Different tool diameters are available, with the smallest being deployable in a 2-inch well. Due to these improvements, NMR has been increasingly used by environmental site investigators worldwide.
  
 
===Cone Penetrometer Profiling===
 
===Cone Penetrometer Profiling===
Cone Penetrometer (CPT) profiling is one of the most common approaches for quantifying the mechanical properties of unconsolidated sediments by geotechnical engineers. CPT profiling measures the mechanical resistance on the tip of a conical probe and on a sleeve immediately behind the tip as the probe is advanced into the subsurface. These measurements are often supplemented by pore water pressure measured behind the cone tip due to its impact on mechanical responses to advancement. These quantities are used to determine the sediment class of the material, from which ''K'' can be estimated using empirical relationships. The resulting values, however, only provide order of magnitude estimates of formation ''K'' <ref>Lunne, T., Robertson, P.K. and Powell, J.J.M., 1997. Cone penetration testing. Geotechnical Practice. Spon Press/Taylor & Francis Group, London and New York.</ref>.
+
[[wikipedia: Cone penetration test | Cone Penetrometer Test (CPT)]] profiling is one of the most common approaches for quantifying the mechanical properties of unconsolidated sediments by geotechnical engineers. CPT profiling measures the mechanical resistance on the tip of a conical probe and on a sleeve immediately behind the tip as the probe is advanced into the subsurface. These measurements are often supplemented by pore water pressure measured behind the cone tip due to its impact on mechanical responses to advancement. These quantities are used to determine the sediment class of the material, from which ''K'' can be estimated using empirical relationships. The resulting values, however, only provide order of magnitude estimates of formation ''K'' <ref>Lunne, T., Robertson, P.K. and Powell, J.J.M., 1997. Cone penetration testing. Geotechnical Practice. Spon Press/Taylor & Francis Group, London and New York.</ref>.
  
 
CPT profiling can be periodically suspended to perform pore pressure dissipation tests to obtain information about the consolidation properties of the formation at selected depths. Hydraulic conductivity can be estimated from these dissipation tests using relationships between ''K'' and the consolidation properties of the formation. Generally speaking, the ''K'' estimates from dissipation tests tend to be more reliable than those from continuous logging. Robertson (2009)<ref>Robertson, P.K., 2009. Interpretation of cone penetration tests-a unified approach. Canadian Geotechnical Journal, 46(11), pp.1337-1355. [https://doi.org/10.1139/T09-065 doi: 10.1139/T09-065]</ref> suggests that it is possible to combine continuous profiles and dissipation tests to improve the quality of CPT ''K'' estimates.  
 
CPT profiling can be periodically suspended to perform pore pressure dissipation tests to obtain information about the consolidation properties of the formation at selected depths. Hydraulic conductivity can be estimated from these dissipation tests using relationships between ''K'' and the consolidation properties of the formation. Generally speaking, the ''K'' estimates from dissipation tests tend to be more reliable than those from continuous logging. Robertson (2009)<ref>Robertson, P.K., 2009. Interpretation of cone penetration tests-a unified approach. Canadian Geotechnical Journal, 46(11), pp.1337-1355. [https://doi.org/10.1139/T09-065 doi: 10.1139/T09-065]</ref> suggests that it is possible to combine continuous profiles and dissipation tests to improve the quality of CPT ''K'' estimates.  
Line 207: Line 200:
 
Formation ''K'' can also be estimated from laboratory grain size analyses on sediment or core samples. A number of empirical and theoretical relationships have been developed for estimation of ''K'' from grain-size statistics. For example, ''K'' can be computed using the relationship developed by Hazen<ref name="Freeze1979" />,
 
Formation ''K'' can also be estimated from laboratory grain size analyses on sediment or core samples. A number of empirical and theoretical relationships have been developed for estimation of ''K'' from grain-size statistics. For example, ''K'' can be computed using the relationship developed by Hazen<ref name="Freeze1979" />,
  
[[File:Liu1w2 Eq5.png|center|400 px]]
+
 
 
{|
 
{|
 
|-
 
|-
|where:
+
| ||[[File:Liu1w2 Eq5.png|330 px]]
 +
|-
 +
|where:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 
|-
 
|-
 
|''d<sub>10</sub>''||is the particle diameter at which 10% of the sample’s mass is comprised of particles with a diameter smaller than this value, and
 
|''d<sub>10</sub>''||is the particle diameter at which 10% of the sample’s mass is comprised of particles with a diameter smaller than this value, and
Line 221: Line 216:
 
More theoretically-based relationships between grain sizes and ''K'' have been developed from the Kozeny-Carman equation that relates formation ''K'' with porosity and grain sizes assuming the porous space can be modeled as tortuous tubes of different diameters. For example, the following relationship can be used to estimate ''K''<ref>Bear, J., 1972. Dynamics of Fluids in Porous Media, Dover, 764 pp</ref>,
 
More theoretically-based relationships between grain sizes and ''K'' have been developed from the Kozeny-Carman equation that relates formation ''K'' with porosity and grain sizes assuming the porous space can be modeled as tortuous tubes of different diameters. For example, the following relationship can be used to estimate ''K''<ref>Bear, J., 1972. Dynamics of Fluids in Porous Media, Dover, 764 pp</ref>,
  
[[File:Liu1w2 Eq6.png|center|400 px]]
 
  
 
{|
 
{|
 
|-
 
|-
|where:
+
| ||[[File:Liu1w2 Eq6.png|400 px]]
 +
|-
 +
|where:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 
|-
 
|-
 
|''n''||is porosity,
 
|''n''||is porosity,
Line 234: Line 230:
 
|}
 
|}
  
The ''K'' estimates from Kozeny-Carman approaches are generally considered to be more accurate than those from the Hazen equation, although both approaches have a number of significant limitations in practice<ref>Carrier III, W.D., 2003. Goodbye, hazen; hello, kozeny-carman. Journal of Geotechnical and Geoenvironmental Engineering, 129(11), pp.1054-1056. [https://doi.org/10.1061/(ASCE)1090-0241(2003)129:11(1054) doi: 10.1061/(ASCE)1090-0241(2003)129:11(1054)]</ref>.  
+
The ''K'' estimates from Kozeny-Carman approaches are generally considered to be more accurate than those from the Hazen equation, although both approaches have a number of significant limitations in practice<ref>Carrier III, W.D., 2003. Goodbye, hazen; hello, kozeny-carman. Journal of Geotechnical and Geoenvironmental Engineering, 129(11), pp.1054-1056. [https://doi.org/10.1061/(ASCE)1090-0241(2003)129:11(1054) doi: 10.1061/(ASCE)1090-0241(2003)129:11(1054)]</ref>.
  
 
==Choosing the Right Approaches==
 
==Choosing the Right Approaches==
Line 246: Line 242:
  
 
==Summary==
 
==Summary==
A variety of methods are available for characterizing ''K'' in environmental site investigations. Historically, pumping tests and slug tests have been the most common methods for ''in situ'' measurement of ''K''.  Recently, two direct push (DP) approaches have become more common for measuring spatial variations in ''K'' in shallow unconsolidated settings: DP injection logging and the DP permeameter. DP injection logging is quick and high resolution, while the DPP provides more reliable ''K'' estimates that are not significantly impacted by a compacted zone around the probe.  
+
A variety of methods are available for characterizing ''K'' in environmental site investigations. Historically, pumping tests and slug tests have been the most common methods for ''in situ'' measurement of ''K''.  Recently, two direct push (DP) approaches have become more common for measuring spatial variations in ''K'' in shallow unconsolidated settings: DP injection logging and the DP permeameter. DP injection logging is quick and high resolution, while the DPP provides more reliable ''K'' estimates that are not significantly impacted by a compacted zone around the probe.
  
 
==References==
 
==References==

Latest revision as of 21:52, 26 April 2022

Hydraulic conductivity (K), a parameter that describes the ease with which water flows in the subsurface, is widely regarded as one of the most important hydrogeologic parameters for environmental site investigations. Mathematically, it is defined as the flowrate per unit area divided by the hydraulic gradient in the direction of flow. Traditional methods for in situ measurement of K include pumping tests and slug tests. More recent work has focused on measuring spatial variations in K and the impact of these variations on the fate and transport of contaminants and on the effectiveness of remediation systems[1][2][3][4].

Related Article(s):


Contributor(s): Dr. Gaisheng Liu and Dr. James Butler, Jr.


Key Resource(s):

Introduction

Hydraulic conductivity is mathematically defined as the parameter K in Darcy’s Law[8] (see Figure 1),

Figure 1. Schematic of Darcy’s Law flow experiment.
Liu1w2 Eq1.png
where:        
Q is the flow rate across area A of a porous medium,
K is the hydraulic conductivity, and
i is the hydraulic gradient, which can be computed as:


Liu1w2 Eq2.png
where:
h1 and h2     are the hydraulic heads at the ends of the experimental domain, and
L is the total length.

Hydraulic conductivity is dependent on the properties of both water and the porous medium,


Liu1w2 Eq3.png
where:        
k is the intrinsic permeability of the medium, a parameter which is solely dependent on the geometry of the interconnected pores,
g is the gravitational constant,
ρ is the density of the pore water, and
μ is the dynamic viscosity of the pore water.

Many approaches have been developed to characterize K[5]. These approaches can be grouped into two general categories based on how the K estimates are obtained:

  • hydraulic methods that involve water or other fluid injection or extraction and the measurement of the induced pressure response, and
  • indirect methods that rely on empirical correlations, often site-specific in nature, between K and other more readily evaluated formation properties (e.g., resistance to electric current).

Because hydraulic methods can be directly related to the mathematical definition of K through Darcy’s Law, K estimates obtained with those methods are generally considered to be more reliable than those obtained with indirect methods.

Other classifications of K characterization approaches are possible. For example, approaches can be divided into those based on data collected in the field and those based on measurements on core samples in the laboratory. In the latter case, there can be considerable uncertainty about how representative the core sample is of field conditions.

Field-based approaches can be further subdivided based on whether the measurement is conducted in the subsurface or on the surface. Subsurface measurements provide the most reliable information about K variations at the scale needed for environmental site investigations. Nearly all field-based hydraulic methods are performed in the subsurface and require either boreholes or direct push (DP) installations for extracting or injecting water and for monitoring the induced head changes. Compared to methods that rely on existing wells, DP approaches can be applied essentially wherever K information is needed, thus greatly expanding the spatial coverage of K measurement in the field. DP approaches, however, are generally limited to use in relatively shallow (20-30 m from land surface) unconsolidated settings.

Field K measurements can also be performed on the land surface or in an aircraft. Most of these are geophysical approaches that involve measuring electrical, magnetic, or seismic properties of the formation and then transforming them into K estimates using empirical relations. The biggest advantage of surface or airborne approaches is their minimal site disturbance, as no subsurface boring is needed. However, the resolution that can be provided is limited, as the measurements are typically affected by conditions over a relatively large volume of the formation. Geophysical methods can also be performed in boreholes. Borehole geophysical methods can provide a much higher resolution description of K than surface-based methods.

The scale of measurement is an important factor to consider when assessing various approaches. For water-supply investigations, a single estimate of K averaged over a large volume of an aquifer will usually suffice; such an estimate is commonly obtained via pumping tests in a well[7]. For water-quality investigations, however, K estimates over a large measurement volume are often of limited value[9]. In that case, small-scale K measurements that provide information about local geological controls on groundwater flow and transport are usually required to obtain reliable predictions of contaminant behavior and to design effective remediation systems.

In this article, we focus on the approaches that can be used to obtain small-scale, localized measurements of K for environmental site investigations. Only the more commonly used approaches are discussed. Some other approaches, such as dipole flow tests[10][11], tracer tests[12], or hydraulic tomography[13][14] will not be discussed here.

Hydraulic Testing Approaches

Pumping Tests

The pumping test is the most common method for determining K over a relatively large volume of an aquifer in water-supply investigations[9]. During a typical test, the pumping rate is kept constant, although it can be varied in time to obtain better signal-to-noise ratios in the acquired data[15]. In the constant-rate approach, a well, preferably centrally located at the site, is pumped while induced head changes are monitored at that and nearby wells. The head changes, along with the pumping rate, can then be used to estimate aquifer parameters using different models of the well-aquifer configuration[16][17]. Pumping test analyses can be facilitated using software packages like AQTESOLV (Aquifer Test Solver), which has been developed for analyzing different types of aquifer tests[18].

Pumping tests are primarily performed to obtain large-scale volumetric averages of aquifer parameters as well as information about aquifer boundaries. However, one form of the pumping test, the step-drawdown test, is specifically directed at getting information about the efficiency of the pumping well. The K estimate from a pumping test is an average over a large volume of the formation and does not provide information at the scale of most relevance for issues involving contaminant transport (e.g., meters or less). Thus, pumping tests are less commonly employed for environmental site investigations[9].

Slug Tests

Figure 2A. Schematic of a Falling-Head Slug Test (water flows from well into aquifer) Performed in a Multi-Level Slug-Test System. Head change is introduced in standpipe, which is directly connected to the interval isolated by the straddle packers (not to scale).
Figure 2B. Example Data plot from a Multilevel Slug Test at GEMS. Test is initiated by sudden depressurization of a pressurized air column[19].

The slug test is one of the most common methods for determining K at the scale of relevance for environmental site investigations[20]. Slug tests are typically performed in existing wells at multiple horizontal and vertical locations across the site to characterize the spatial distribution of K values. However, slug tests can also be performed in temporary installations such as DP rods. In this approach, a near-instantaneous head change is applied in the well or DP rods and the subsequent head recovery is used to estimate K. The initial head change can be introduced using a solid object (slug), compressed gas (pneumatic system), or by the addition/removal of a certain amount of water[20]. The recovery data can be analyzed using different models of the well-formation configuration[18][20]. In contrast to pumping tests, slug tests provide a K estimate that is primarily a function of the materials in the vicinity of the screened interval of the test well. Slug test data analyses can be facilitated using software packages like AQTESOLV (Aquifer Test Solver), which has been developed for analyzing different test methods[18].

Although a slug test is very simple in principle, considerable care must be taken in all stages of a test[5][20]. Because test responses are highly sensitive to the materials immediately adjacent to the test well, that well must be appropriately developed before a test is performed. Otherwise, the K estimate can be biased by a low-K zone (skin) that can form during well construction. In high-K formations, the head recovery is rapid, and pneumatic methods are often used to minimize the time associated with test initiation. In addition, due to the relatively high flow velocity, the impact of pipe hydraulics on pressure readings should be considered when analyzing slug test data from high-K zones[19]. In low-K formations, slug tests can take an extremely long time to complete, although test time can be significantly reduced by decreasing the effective casing radius (portion of well in which the water level is changing)[20].

Slug tests can be configured to obtain information about vertical variations in K along the screened (open) interval of a well. Test intervals can be isolated with straddle packers, and slug tests performed within that isolated interval. The straddle packers can be incrementally moved along the screened interval of a well to characterize the vertical variation of K at a relatively high resolution. Using a two-packer tool (Figure 2A), slug tests have been performed in a number of 0.25-m intervals in a well at the Geohydrologic Experimental and Monitoring Site (GEMS) in the Kansas River valley[5]. At each isolated interval, multiple tests are performed, initiated with different head changes, following recommended test guidelines[20]. An example data set is presented in Figure 2B. The multi-level slug test K estimates compare favorably with estimates obtained using other approaches (Figure 3).

Figure 3. K Estimates from Different Field Methods at GEMS[5]. Well DW is located 2 m east of GEMS4S.
Figure 4. Schematic Illustrating Hydrostratigraphic Profiling, the combination of electrical conductivity (EC) and multilevel slug test profiling[21].

In addition to existing wells, slug tests can also be performed in DP installations so that K estimates can be obtained virtually at any location in unconsolidated formations. Various approaches have been developed that allow slug tests to be performed at one or multiple levels in a single DP hole[22][23][24][21]. In McCall et al. (2002)[24], a pair of nested rod strings were driven to the test interval with a solid drive point attached to the end of inner rod string for advancement. Upon reaching the test depth, the drive point and inner rod string were retracted, and a screen was lowered to the bottom of the outer rod string. The outer rod string was pulled up while the screen was held in place, leaving the screen exposed to the surrounding formation. After the slug test was completed in the exposed screen, the screen was removed and the inner rod string with the attached solid drive point was reinserted. The nested rod strings were then driven to the next test depth. In low-K formations such as silts and clays, the formation materials may not collapse completely back to the screen when the outer rod string is pulled up. In this case, the diameter of the borehole can be estimated and used in place of the screen diameter, or the problem can be avoided by using a coring tube of similar size to the screen to create a hole below the end of the outer rod string. Instead of setting the screen by pulling up the outer rod string, the screen can be directly inserted into the hole for slug testing. Regardless of how the screen is set into the formation, it is always recommended that the screen be appropriately developed before slug tests are performed. In low-K formations, development may be limited to scraping the sides of the cored hole with a steel brush and removing the silty water from the screen with a low-flow pump.

Sellwood et al. (2005)[21] proposed a modification of the approach by McCall et al. (2002)[24] to reduce profiling time and gain more information about subsurface stratigraphy. Instead of performing slug tests on the way down, slug tests were performed at different depths as the outer rod string was pulled up (Figure 4). This way, the number of changes between the outer and inner rod strings was minimized (i.e., only one change needed at the bottom of the profile). Furthermore, an electrical conductivity (EC) probe was attached to the inner rod string so that EC data could be collected as the rod strings were advanced; the EC data could then be used for selecting the intervals for slug tests as the rod strings were retracted. This approach, called hydrostratigraphic profiling, enabled information on electrical and hydraulic conductivity to be obtained at a speed and resolution that had previously not been possible.

Borehole Flowmeter Profiling

Figure 5. Schematic Diagram showing the borehole Flowmeter Profiling Procedure[25].

Borehole flowmeter profiling is one of the most efficient approaches for characterizing the vertical variations of K[25]. This approach involves pumping a well at a constant rate while measuring the vertical flow rate within the screened interval (Figure 5). A flowmeter is initially positioned at the bottom of the screen. After the pump is turned on and the pumping rate stabilizes, the flowmeter is gradually moved up, often in an incremental fashion with short stops, although continuous profiling is also done. After the entire screened interval is measured, the profile of vertical flow rates versus depths can be used for K estimation. Typically, two profiles are performed, one prior to pumping and one during pumping. The profile obtained during pumping is corrected with the profile prior to pumping to remove the impacts of ambient flow in the well. Due to the difficulty of detecting small differences in flow rates, the K estimates from flowmeter profiling are subject to more noise in less permeable settings[26]. In addition, the results of flowmeter profiling are sensitive to in-well hydraulics, such as the impedance of flowmeter to vertical flow and disturbance by the pump intake when it is too close to the measurement interval[5].

Different types of flowmeters have been used for borehole profiling[27]. For example, electromagnetic flowmeters are often used in environmental site investigations due to their lower flow rate detection limit and lower head loss (high head loss across the flowmeter can cause flow bypass through the filter pack of the well). However, they may be more expensive than other options and have certain practical limitations as well[28]. Heat pulse flowmeters are also used for profiling in less-permeable settings.

The most common method of estimating K from borehole flowmeter profiles is based on the assumption that the lateral flow rate of an interval (zone between flow-rate measurements) is proportional to the K and thickness of that interval in a well fully screened in a perfectly layered aquifer,

Liu1w2 Eq4.png
where:        
Ki is the K of test interval i;
Ka is the average K for the entire screened interval;
ΔQi is the net lateral inflow across interval i, which can be calculated by subtracting the flow rate at the bottom of the interval from that at the top and taking ambient flow, if any, into account;
Qt is the total pumping rate;
ΔBi is the thickness of interval i; and
Bt is the total thickness of the aquifer.

Ka can be obtained from pumping or slug tests at the same well, although it is generally recommended that slug tests be used, as the average K in the immediate vicinity of the well may differ from that determined from a pumping test in a laterally heterogeneous aquifer[5].

Figure 3 compares the results of borehole flowmeter profiling and multilevel slug tests at GEMS. Despite some local differences, the K estimates from borehole flowmeter profiling and multilevel slug tests are quite similar at that well.

Figure 6A. Schematic of the DPIL with a Screen Attached to the Probe Rod at its Lower End[6].

DP Injection Logging

DP technology has shown great promise for characterization of K variations in shallow unconsolidated formations at the resolution, accuracy, and speed that are critically needed for practical investigations[29][30][31]. Over the last few decades, a series of DP methods and probes have been developed for obtaining information about vertical variations in K[5][6][32].

DP Injection Logging (DPIL) is one of the most powerful approaches for K characterization (in terms of resolution and speed of acquisition) and is also one of the most widely used techniques in environmental site investigations[6][33]. This approach consists of advancing a probe with a single screened port (Figure 6A). Water is injected continuously through the screen while the probe is advanced and the pressure response to injection is monitored behind the screen or at the surface. The profile of the ratio of injection rate to injection pressure is closely related to vertical variations in K. However, factors other than the hydraulic conductivity can potentially affect the pressure response (e.g., injection line losses, probe advancement speed, formation alteration, and variations in specific storage), so DPIL ratios only provide relative information about K variations. Methods are needed for transforming DPIL ratios into actual K estimates; these methods typically involve correlations with nearby K data.

Figure 6B. Artistic Rendering of the HPT (continuous DPIL probe combined with an EC Wenner array) [34].

There are currently two variants of DPIL: discontinuous[32][34] and continuous[30][35][36] modes. In discontinuous DPIL mode, probe advancement is briefly halted at the desired depth and the injection rate is then varied in a step-wise fashion while injection rates and pressures are measured. In continuous DPIL mode, injection rate and pressure are measured as the probe is continuously advanced. Compared to continuous DPIL, discontinuous DPIL has a few advantages[6]:

  • pore water pressure changes generated by probe advancement, which can be difficult to characterize, have little influence on the measured injection pressures;
  • background hydrostatic pressures have no influence on the measured responses; and
  • use of multiple injection rates at each depth allows a better assessment of the formation response to injection.

However, due to the need to halt advancement, discontinuous DPIL requires significantly more time than continuous DPIL, and, as a result, resolution below 10 - 20 cm is rarely possible. In common practice, continuous DPIL is generally preferred due to its speed (a 20-m profile requires about 2 hours) and resolution (one measurement per 1.5 cm vertical interval). A continuous DPIL probe, supplemented with an EC sensor, is commercially available from Geoprobe Systems under the tradename of the Hydraulic Profiling Tool (HPT; Figure 6B).

The HPT is most effective in moderately permeable formations with a K range of 0.03 to 10 m/d, although the detection range can be improved by modifying equipment and operating procedures[6]. In standard HPT practice, the probe advancement rate is 2 cm/s and the continuous water injection rate is about 300 mL/min. In high K formations (e.g., K > 10 m/d), the injection pressure response to the standard injection rate is generally too small to be reliably measured. A larger injection rate can be used to increase the pressure signal. However, the line loss between the transducer and injection screen, which increases quadratically with the injection rate, may become significant and require additional step tests to remove its impact[6]. On the other hand, in low K formations (e.g., K < 0.03 m/d), the pressure increase from probe advancement and water injection may become so large that it will cause formation alterations and exceed the upper measurement limit of the transducer. A recent study has demonstrated that by reducing both probe advancement speed and flow injection rate, HPT can be applied to formations with K as low as 0.003 m/d[37].

Transforming DPIL profiles into K estimates is based on empirical relationships[30][32][36] [38]. Although a general relationship has been presented[38], site-specific calibrations are typically used to improve the reliability of K estimates. The site-specific calibrations can be performed by comparing the DPIL profiles to nearby K estimates obtained via other means[32][34] or by a modeling approach that directly combines the DPIL data with collocated hydraulic tests such as seen in Figure 6C, which shows a series of K profiles from continuous DPIL as calibrated by collocated hydraulic tests at an alluvial aquifer in Mississippi[36].

Figure 6C. Six Example DPIL (red curve) K Profiles Calibrated by Collocated Direct Push Permeameter Tests (blue bars) at the MacroDispersion Experiment (MADE) Site, Columbus, MS. The dashed curves represent the upper K limit for continuous DPIL under the given tool operating procedure[36].

There are two major assumptions invoked for estimating K from DPIL profiles. First, near-rod material compaction from probe advancement has a similar impact on K at all depths of the profile. Second, the impact of formation specific storage is either similar across the profile or can be ignored. Both assumptions appear to be valid for sands and gravels. In silts and clays, particularly those with high compressibility, considerable caution must be used to assess the viability of these assumptions[6][37].

DP Permeameter

Figure 7A. Schematic of the DPP (not to scale).
Figure 7B. Example DPP Test Sequence at GEMS[6].
Figure 8. Comparison of DPP and DP Multilevel Slug Test K profiles at GEMS. The slug test profile is located 2 m northwest of DPP profile 1, and 2 m southeast of DPP profile 2[5].

The Direct Push Permeameter (DPP) consists of a short cylindrical screen with two pressure transducers inset into the probe at short distances above the screen (Figure 7A)[5][6][39]. The DPP probe is advanced to a depth at which a K estimate is needed and then a series of short-term injection tests are performed. K is then typically estimated from the spherical form of Darcy’s Law using the injection rate and the injection-induced pressure responses at the two transducers. The resulting estimate is a weighted average over the interval between the screen and the farthest transducer (approx. 0.4 m in Figure 7A); material outside of that interval has little influence [40]. Horizontally, the DPP shows little sensitivity to conditions greater than ∼ 0.5 m away from the tool. The DPP test procedure only requires that steady-shape (constant hydraulic gradient) flow conditions be attained. True steady-state conditions are not required, which can lead to a significant reduction in the measurement time. A DPP test sequence generally requires 10 – 15 minutes per measurement interval in moderate to high-K formations[39].

Figure 7B shows an example DPP test sequence at GEMS[6]. Three tests with different rates are used to assess the reliability of the K estimates. The first injection rate is typically based on the expected K value in the test interval (e.g., 800 mL/min for sands). Depending on the magnitude of the pressure response in the first test, the second injection rate can be set higher or lower than the first (if pressure response is high in the first test, the flow rate in the second test is reduced, and vice versa). The third flow rate is set similar to that in the first test to check whether a similar pressure response can be obtained. The K values computed using Darcy’s Law were 66.0, 63.1 and 65.4 m/d, respectively, for the three tests in Figure 7B. A set of consistent K estimates is the goal; a lack of consistency could indicate DPP system instability or formation alteration. Figure 8 shows the comparison of DPP K estimates with those obtained from DP multilevel slug tests at GEMS[39]. The results from DPP profiling have proven extremely useful for transforming collocated DPIL profiles into K estimates[26][36].

The most significant advantage of the DPP is that, unlike all other hydraulic tests discussed here, screen clogging and near-rod compaction have a very small impact on K estimates. The insensitivity to the zone of compaction is a result of head-equalizing cross-flow between the thin compacted zone and the surrounding formation[39][40]. The limited lateral sensitivity of the DPP (less than 0.5 m from the probe) is also one of its primary advantages, as it allows the tool to provide high-resolution profiles of vertical K variations in the immediate vicinity of the probe.

Laboratory Permeameter Tests on Core Samples

Laboratory permeameter tests of sediment or core samples collected during drilling are also used for acquiring information about formation K. Different hydraulic systems can be used, depending on the expected K value of the sample. In general, when the sample K is high, constant heads are used for both the influent and effluent tubes to maintain a constant hydraulic gradient across the sample; however when K is low, falling heads are used for the influent tube with either constant or rising heads in the effluent tube. By measuring both the hydraulic gradient and flow rate, Darcy’s Law can be used to calculate the sample K. In laboratory permeameter tests, considerable care must be given to all stages of the work, including sample preparation and installation into the test cell, use of appropriate ambient stress for materials with high swell potential, de-airing of the hydraulic system and sample, mobilization and resettling of fine materials, head losses in the tubes, non-Darcian flow, and biochemical reactions due to changes in temperature and pore water composition[41][42].

As discussed earlier, the question of how representative of field conditions samples are should always be considered when assessing the K estimates from laboratory analyses. Larger-scale preferential flow features, such as macropores and fractures, are usually not captured adequately by the samples, so K estimates from laboratory permeameter tests should be considered as the lower end of the range for field values. Field-based characterization methods of K are preferred over laboratory-based methods in most environmental site investigations.

Indirect Approaches

Geophysical Methods

Estimates for formation K can also be inferred from the measurement of more readily evaluated geophysical properties using empirical, often site-specific, relationships. The geophysical properties that may be used for K estimation include, but are not limited to, electrical conductivity (EC) (or the reverse, electrical resistivity), the hydrogen response to nuclear magnetic resonance (NMR), and natural gamma radiation. In addition to K estimation, geophysical methods have also been used for many purposes in environmental site investigation, assessment and monitoring (see related article –Geophysical Methods).

Borehole EC logging measures vertical variations in the bulk EC of the surrounding formation, which can then be used to estimate K based on empirical relationships[43][44]. EC has been adapted to DP equipment for more flexible assessment of subsurface conditions than the approaches that rely on existing wells[45]. Variations in the EC of saturated media are primarily a function of pore-fluid chemistry, clay content, and total porosity. If variations in pore-fluid chemistry can be neglected, EC profiling can be used to delineate the major hydrostratigraphic features, and thus large-scale variations in K[46]. However, when electrically-conductive clay is absent, EC profiling generally provides little information about vertical variations in K in silts, sands and gravels.

Borehole natural gamma logging provides a record of natural gamma radiation versus depth. This radiation is quantified by counting the gamma particles passing through a scintillation crystal in a certain time interval. A high natural gamma reading is generally associated with clay-rich intervals, while a low reading is generally associated with sands and gravels[47]. Like EC logging, natural gamma logging provides an effective measurement of formation characteristics when clay is present, but is of less use in silts, sands and gravels.

Figure 9A. Measurement Domain around a NMR Borehole Probe
Figure 9B. NMR Responses as a Function of Pore Size, Where A is the Initial Magnetization, T2 is NMR Relaxation Time[48].

Efforts have been made to adapt nuclear magnetic resonance (NMR) profiling, a widely used borehole technique in the petroleum industry, to environmental site investigations[49][48] [50]. This approach involves measuring the response of hydrogen atoms (protons) to a series of imposed magnetic fields at radio frequency. The response is a function of, among other things, the pore-size distribution of material in a thin cylindrical ring centered on the probe (Figures 9A and 9B). K is then estimated from the pore-size distribution information using different empirical relationships, all of which have their origins in the Kozeny-Carman equation[50]. Figure 9C compares the NMR results with DPP K estimates at GEMS.

Figure 9C. Comparison of NMR and DPP K estimates at two GEMS locations. The dash lines show the range of the NMR K estimates due to the uncertainty in the empirical parameters used in the pore size - K relationship[50].

The NMR logging tools used in the petroleum industry are typically too large for use by the hydrological community (e.g., tools are typically well over 10 ft long and have a diameter larger than 5 inches). The equipment costs are also prohibitively high for hydrological users. After recent technology advancements, NMR measurement can be made much more effectively both at the surface and in a borehole, including with DP equipment[49][48][50]. The equipment costs are significantly decreased compared to those used in the petroleum industry. The measurement time is reduced down to a few minutes per interval. The vertical sampling interval of most current logging tools is about 0.5 m. Different tool diameters are available, with the smallest being deployable in a 2-inch well. Due to these improvements, NMR has been increasingly used by environmental site investigators worldwide.

Cone Penetrometer Profiling

Cone Penetrometer Test (CPT) profiling is one of the most common approaches for quantifying the mechanical properties of unconsolidated sediments by geotechnical engineers. CPT profiling measures the mechanical resistance on the tip of a conical probe and on a sleeve immediately behind the tip as the probe is advanced into the subsurface. These measurements are often supplemented by pore water pressure measured behind the cone tip due to its impact on mechanical responses to advancement. These quantities are used to determine the sediment class of the material, from which K can be estimated using empirical relationships. The resulting values, however, only provide order of magnitude estimates of formation K [51].

CPT profiling can be periodically suspended to perform pore pressure dissipation tests to obtain information about the consolidation properties of the formation at selected depths. Hydraulic conductivity can be estimated from these dissipation tests using relationships between K and the consolidation properties of the formation. Generally speaking, the K estimates from dissipation tests tend to be more reliable than those from continuous logging. Robertson (2009)[52] suggests that it is possible to combine continuous profiles and dissipation tests to improve the quality of CPT K estimates.

There are a few practical challenges with CPT profiling. First, it is difficult to keep the porous element for pore pressure measurement fully saturated while the probe is driven through the vadose zone. Clogging of the porous element may also be a concern at sites with a significant amount of fine materials. Similar to DPIL, formation alteration can produce a significant impact on the pore pressure. However, unlike the HPT where the pressure port is located farther away from the tip, the porous element is immediately behind the tip so that the impact of probe advancement on the CPT pore pressure measurement will be much greater than with the HPT.

Laboratory Grain Size Analyses on Core Samples

Formation K can also be estimated from laboratory grain size analyses on sediment or core samples. A number of empirical and theoretical relationships have been developed for estimation of K from grain-size statistics. For example, K can be computed using the relationship developed by Hazen[7],


Liu1w2 Eq5.png
where:        
d10 is the particle diameter at which 10% of the sample’s mass is comprised of particles with a diameter smaller than this value, and
C is an empirically-defined coefficient dependent on grain size and sorting.

The value of C varies with the units for K, and has a large range of published values reported by different researchers.

More theoretically-based relationships between grain sizes and K have been developed from the Kozeny-Carman equation that relates formation K with porosity and grain sizes assuming the porous space can be modeled as tortuous tubes of different diameters. For example, the following relationship can be used to estimate K[53],


Liu1w2 Eq6.png
where:        
n is porosity,
K is in m/d, and
dm is the geometric mean grain size in mm, calculated as (d84 X d16)0.5.

The K estimates from Kozeny-Carman approaches are generally considered to be more accurate than those from the Hazen equation, although both approaches have a number of significant limitations in practice[54].

Choosing the Right Approaches

Due to the large uncertainty associated with core sample representativeness, field-based approaches are preferred for acquiring information about K in environmental site investigations. Under certain conditions, such as sites that are primarily composed of low-K silts and clays, hydraulic testing may be time consuming so that only a limited number of measurements can be obtained in the field. In this case, coring and subsequent laboratory analyses can be used as a complement to increase the spatial coverage of K measurements across the site.

For field-based approaches, hydraulic tests provide more reliable K estimates than indirect approaches such as geophysical methods. However, geophysical methods are generally more time efficient and can provide more information about the spatial continuity of subsurface flow and transport features. There is a growing interest in the joint use of hydraulic tests and geophysical methods, so that high-resolution characterization of K can be obtained across the entirety of a site, rather than just in the immediate vicinity of the limited number of boreholes.

Pumping tests provide K estimates averaged over a large volume of formation and do not provide information at the scale relevant to solute transport, so they are not recommended in environmental site investigations directed at getting information about spatial variations in K. Slug tests can be used to obtain high-quality K estimates at the relevant scale, but considerable care must be exercised in all stages of the work, including the proper development of the measurement interval prior to the tests. Borehole flowmeter profiling is sensitive to in-well hydraulics and generally does not provide high-quality K estimates in less permeable settings.

DP methods have great promise for measurement of spatial variations in K in shallow unconsolidated formations. Because the tools are small in diameter and no materials are removed during advancement, DP methods generate much less site disturbance than traditional borehole techniques (e.g., rotary drilling). Two of the most promising DP K methods are DP injection logging and the DP permeameter. DP injection logging is quick and high resolution, but is considered to only provide relative information on K because factors other than K can also affect the pressure response. The DPP provides reliable K estimates and is more time efficient than other formal hydraulic test methods. The most significant advantage of DPP is that a low-K zone formed by advancement-induced compaction has little impact on the estimated K due to head-equalizing cross-flow between the thin compacted zone and the surrounding formation.

Summary

A variety of methods are available for characterizing K in environmental site investigations. Historically, pumping tests and slug tests have been the most common methods for in situ measurement of K. Recently, two direct push (DP) approaches have become more common for measuring spatial variations in K in shallow unconsolidated settings: DP injection logging and the DP permeameter. DP injection logging is quick and high resolution, while the DPP provides more reliable K estimates that are not significantly impacted by a compacted zone around the probe.

References

  1. ^ Boggs, J.M., Young, S.C., Beard, L.M., Gelhar, L.W., Rehfeldt, K.R. and Adams, E.E., 1992. Field study of dispersion in a heterogeneous aquifer: 1. Overview and site description. Water Resources Research, 28(12), pp.3281-3291. doi: 10.1029/92WR01756
  2. ^ Dagan, G. and Neuman, S.P. eds., 1997. Subsurface flow and transport: a stochastic approach. Cambridge University Press.
  3. ^ Fogg, G.E., Carle, S.F. and Green, C., 2000. Connected-network paradigm for the alluvial aquifer system. Special Papers-Geological Society of America, pp.25-42.
  4. ^ Haslauer, C.P., Bárdossy, A. and Sudicky, E.A., 2017. Detecting and modelling structures on the micro and the macro scales: Assessing their effects on solute transport behaviour. Advances in water resources, 107, pp.439-450. doi: 10.1016/j.advwatres.2017.05.007
  5. ^ 5.0 5.1 5.2 5.3 5.4 5.5 5.6 5.7 5.8 5.9 Butler, J.J., 2005. Hydrogeological methods for estimation of spatial variations in hydraulic conductivity. In Hydrogeophysics (pp. 23-58). Springer, Dordrecht. doi: 10.1007/1-4020-3102-5_2
  6. ^ 6.00 6.01 6.02 6.03 6.04 6.05 6.06 6.07 6.08 6.09 6.10 Liu, G., Butler, J.J., Reboulet, E. and Knobbe, S., 2012. Hydraulic conductivity profiling with direct push methods. Grundwasser, 17(1), pp.19-29. doi: 10.1007/s00767-011-0182-9
  7. ^ 7.0 7.1 7.2 Freeze, R.A., Cherry, J.A., 1979. Groundwater. Prentice-Hall. Inc., Englewood Cliffs, NJ ISBN 0-13-365312-9.
  8. ^ Darcy, H. (1856). Les Fontaines Publiques de la Ville de Dijon, Dalmont, Paris. doi: 10.1029/2001WR000727
  9. ^ 9.0 9.1 9.2 Butler Jr, J.J., 2009. Pumping tests for aquifer evaluation - Time for a change?. Groundwater, 47(5), pp.615-617. doi: 10.1111/j.1745-6584.2008.00488.x
  10. ^ Kabala, Z.J., 1993. The dipole flow test: A new single‐borehole test for aquifer characterization. Water Resources Research, 29(1), pp.99-107. doi: 10.1029/92WR01820
  11. ^ Zlotnik, V.A. and Zurbuchen, B.R., 1998. Dipole probe: Design and field applications of a single‐borehole device for measurements of vertical variations of hydraulic conductivity. Groundwater, 36(6), pp.884-893. doi: 10.1111/j.1745-6584.1998.tb02095.x
  12. ^ Datta‐Gupta, A., Yoon, S., Vasco, D.W. and Pope, G.A., 2002. Inverse modeling of partitioning interwell tracer tests: A streamline approach. Water Resources Research, 38(6), pp.15-1. doi: 10.1029/2001WR000597
  13. ^ Yeh, T.C.J. and Liu, S., 2000. Hydraulic tomography: Development of a new aquifer test method. Water Resources Research, 36(8), pp.2095-2105. Report.pdf
  14. ^ Bohling, G.C. and Butler Jr, J.J., 2010. Inherent limitations of hydraulic tomography. Groundwater, 48(6), pp.809-824. doi: 10.1111/j.1745-6584.2010.00757.x
  15. ^ Rasmussen, T.C., Haborak, K.G. and Young, M.H., 2003. Estimating aquifer hydraulic properties using sinusoidal pumping at the Savannah River site, South Carolina, USA. Hydrogeology Journal, 11(4), pp.466-482. doi: 10.1007/s10040-003-0255-7
  16. ^ Batu, V., 1998. Aquifer hydraulics: a comprehensive guide to hydrogeologic data analysis. John Wiley & Sons.
  17. ^ Kruseman, G.P., De Ridder, N.A. and Verweij, J.M., 1990. Analysis and evaluation of pumping test data – ILRI Pub. 47. The Netherlands: International institute for land reclamation and improvement
  18. ^ 18.0 18.1 18.2 Duffield, G.M., 2007. AQTESOLV for Windows Version 4.5 User's Guide. HydroSOLVE, Reston, VA.
  19. ^ 19.0 19.1 Butler Jr, J.J., Garnett, E.J. and Healey, J.M., 2003. Analysis of slug tests in formations of high hydraulic conductivity. Groundwater, 41(5), pp.620-631. doi: 10.1111/j.1745-6584.2003.tb02400.x
  20. ^ 20.0 20.1 20.2 20.3 20.4 20.5 Butler Jr, J.J., 1997. The design, performance, and analysis of slug tests. CRC Press.
  21. ^ 21.0 21.1 21.2 Sellwood, S.M., Healey, J.M., Birk, S. and Butler, J.J., 2005. Direct‐push hydrostratigraphic profiling: coupling electrical logging and slug tests. Ground water, 43(1), pp.19-29. doi: 10.1111/j.1745-6584.2005.tb02282.x
  22. ^ Hinsby, K., Bjerg, P.L., Andersen, L.J., Skov, B. and Clausen, E.V., 1992. A mini slug test method for determination of a local hydraulic conductivity of an unconfined sandy aquifer. Journal of Hydrology, 136(1-4), pp.87-106. doi: 10.1016/0022-1694(92)90006-H
  23. ^ Butler Jr, J.J., Healey, J.M., McCall, G.W., Garnett, E.J. and Loheide, S.P., 2002. Hydraulic tests with direct‐push equipment. Groundwater, 40(1), pp.25-36. doi: 10.1111/j.1745-6584.2002.tb02488.x
  24. ^ 24.0 24.1 24.2 McCall, W., J.J. Butler, Jr., J.M. Healey, A.A. Lanier, S.M. Sellwood, and E.J. Garnett, 2002. A dual-tube direct-push method for vertical profiling of hydraulic conductivity in unconsolidated formations, Environ. & Eng. Geoscience, 8(2), 75-84. doi: 10.2113/gseegeosci.8.2.75
  25. ^ 25.0 25.1 Molz, F.J., Morin, R.H., Hess, A.E., Melville, J.G. and Güven, O., 1989. The impeller meter for measuring aquifer permeability variations: evaluation and comparison with other tests. Water Resources Research, 25(7), pp.1677-1683. doi: 10.1029/WR025i007p01677
  26. ^ 26.0 26.1 Bohling, G.C., Liu, G., Knobbe, S.J., Reboulet, E.C., Hyndman, D.W., Dietrich, P. and Butler, J.J., 2012. Geostatistical analysis of centimeter‐scale hydraulic conductivity variations at the MADE site. Water Resources Research, 48(2). doi: 10.1029/2011WR010791
  27. ^ Young, S.C. and Pearson, H.S., 1995. The electromagnetic borehole flowmeter: Description and application. Groundwater Monitoring & Remediation, 15(4), pp.138-147. doi: 10.1111/j.1745-6592.1995.tb00561.x
  28. ^ Bomana, G.K., Molz, F.J. and Boonec, K.D., 1997. Borehole flowmeter application in fluvial sediments: Methodology, results, and assessment. Groundwater, 35(3), pp.443-450. doi: 10.1111/j.1745-6584.1997.tb00104.x
  29. ^ Dietrich, P. and Leven, C., 2009. Direct push-technologies. In Groundwater geophysics (pp. 347-366). Springer, Berlin, Heidelberg.
  30. ^ 30.0 30.1 30.2 McCall, W., Nielsen, D.M., Farrington, S.P. and Christy, T.M., 2005. Use of direct-push technologies in environmental site characterization and ground-water monitoring. In Practical Handbook of Environmental Site Characterization and Ground-Water Monitoring, Second Edition (pp. 355-482). CRC Press. doi: 10.1201/9781420032246.ch6
  31. ^ Leven, C., Weiß, H., Vienken, T. and Dietrich, P., 2011. Direct-Push-Technologien–Effiziente Untersuchungsmethoden für die Untergrunderkundung. Grundwasser, 16(4), pp.221-234. doi: 10.1007/s00767-011-0175-8
  32. ^ 32.0 32.1 32.2 32.3 Dietrich, P., Butler Jr, J.J. and Faiß, K., 2008. A rapid method for hydraulic profiling in unconsolidated formations. Groundwater, 46(2), pp.323-328. doi: 10.1111/j.1745-6584.2007.00377.x
  33. ^ Maliva, R.G., 2016. Direct-push technology. In Aquifer Characterization Techniques (pp. 383-402). Springer, Cham. doi: 10.1007/978-3-319-32137-0_12
  34. ^ 34.0 34.1 34.2 Lessoff, Steven C., Uwe Schneidewind, Carsten Leven, Philipp Blum, Peter Dietrich, and Gedeon Dagan. "Spatial characterization of the hydraulic conductivity using direct‐push injection logging." Water Resources Research 46, no. 12 (2010). doi: 10.1029/2009WR008949
  35. ^ McCall, W., Christy, T.M., Christopherson, T. and Issacs, H., 2009. Application of direct push methods to investigate uranium distribution in an alluvial aquifer. Groundwater Monitoring & Remediation, 29(4), pp.65-76. doi: 10.1111/j.1745-6592.2009.01258.x
  36. ^ 36.0 36.1 36.2 36.3 36.4 Liu, G., Butler, J.J., Bohling, G.C., Reboulet, E., Knobbe, S. and Hyndman, D.W., 2009. A new method for high‐resolution characterization of hydraulic conductivity. Water Resources Research, 45(8). doi: 10.1029/2009WR008319
  37. ^ 37.0 37.1 Liu, G., Borden, R.C. and Butler Jr, J.J., 2018. Simulation Assessment of Direct Push Injection Logging for High‐Resolution Aquifer Characterization. Groundwater. doi: 10.1111/gwat.12826
  38. ^ 38.0 38.1 McCall, W. and Christy, T.M., 2010. Development of a hydraulic conductivity estimate for the Hydraulic Profiling Tool (HPT). In The 2010 North American Environmental Field Conference & Exposition: Conference Program with Abstracts: Session VII.
  39. ^ 39.0 39.1 39.2 39.3 Butler Jr, J.J., Dietrich, P., Wittig, V. and Christy, T., 2007. Characterizing hydraulic conductivity with the direct‐push permeameter. Groundwater, 45(4), pp.409-419. [https://doi.org/10.1111/j.1745-6584.2007.00300.x doi: 10.1111/j.1745-6584.2007.00300.x
  40. ^ 40.0 40.1 Liu, G., Bohling, G.C. and Butler, J.J., 2008. Simulation assessment of the direct‐push permeameter for characterizing vertical variations in hydraulic conductivity. Water Resources Research, 44(2). doi:10.1029/2007WR006078
  41. ^ Klute, A. and Dirksen, C., 1986. Hydraulic conductivity and diffusivity: Laboratory methods. Methods of Soil Analysis: Part 1 - Physical and Mineralogical Methods, (Methodsofsoilan1), pp.687-734. doi: 10.2136/sssabookser5.1.2ed.c28
  42. ^ American Society for Testing and Materials (ASTM) International, 2016. D5084. Standard test methods for measurement of hydraulic conductivity of saturated porous materials using a flexible wall permeameter. West Conshohocken, PA, USA. doi: 10.1520/d5084-16a
  43. ^ Purvance, D.T. and Andricevic, R., 2000. On the electrical‐hydraulic conductivity correlation in aquifers. Water Resources Research, 36(10), pp.2905-2913. doi: 10.1029/2000wr900165
  44. ^ Slater, L. and Lesmes, D.P., 2002. Electrical‐hydraulic relationships observed for unconsolidated sediments. Water Resources Research, 38(10), pp.31-1. doi:10.1029/2001WR001075
  45. ^ Christy, C.D., Christy, T.M. and Wittig, V., 1994, May. A percussion probing tool for the direct sensing of soil conductivity. In Proceedings of the 8th National Outdoor Action Conference (Vol. 38, pp. 1-394). Dublin, Ohio: National Ground Water Association.
  46. ^ Schulmeister, M.K., Butler, J.J., Healey, J.M., Zheng, L., Wysocki, D.A. and McCall, G.W., 2003. Direct‐push electrical conductivity logging for high‐resolution hydrostratigraphic characterization. Groundwater Monitoring & Remediation, 23(3), pp.52-62. doi: 10.1111/j.1745-6592.2003.tb00683.x
  47. ^ Keys, W.S., 1990. Borehole geophysics applied to ground-water investigations (No. 02-E2). USGPO; For sale by the Books and Open-File Reports Section, US Geological Survey. doi: 10.3133/twri02E2
  48. ^ 48.0 48.1 48.2 Walsh, D., Turner, P., Grunewald, E., Butler, J., Knight, R., Reboulet, E., Knobbe, S., Christy, T. and McCall, W., 2011, January. Field demonstration of nuclear magnetic resonance (NMR) logging tools for groundwater and environmental investigations. In Symposium on the Application of Geophysics to Engineering and Environmental Problems 2011 (pp. 69-69). Society of Exploration Geophysicists. doi: 10.4133/1.3614224
  49. ^ 49.0 49.1 Walsh, D., Turner, P., Grunewald, E., Zhang, H., Butler Jr, J.J., Reboulet, E., Knobbe, S., Christy, T., Lane Jr, J.W., Johnson, C.D. and Munday, T., 2013. A small‐diameter NMR logging tool for groundwater investigations. Groundwater, 51(6), pp.914-926. doi: 10.1111/gwat.12024
  50. ^ 50.0 50.1 50.2 50.3 Knight, R., Walsh, D.O., Butler Jr, J.J., Grunewald, E., Liu, G., Parsekian, A.D., Reboulet, E.C., Knobbe, S. and Barrows, M., 2016. NMR logging to estimate hydraulic conductivity in unconsolidated aquifers. Groundwater, 54(1), pp.104-114. doi: 10.1111/gwat.12324
  51. ^ Lunne, T., Robertson, P.K. and Powell, J.J.M., 1997. Cone penetration testing. Geotechnical Practice. Spon Press/Taylor & Francis Group, London and New York.
  52. ^ Robertson, P.K., 2009. Interpretation of cone penetration tests-a unified approach. Canadian Geotechnical Journal, 46(11), pp.1337-1355. doi: 10.1139/T09-065
  53. ^ Bear, J., 1972. Dynamics of Fluids in Porous Media, Dover, 764 pp
  54. ^ Carrier III, W.D., 2003. Goodbye, hazen; hello, kozeny-carman. Journal of Geotechnical and Geoenvironmental Engineering, 129(11), pp.1054-1056. doi: 10.1061/(ASCE)1090-0241(2003)129:11(1054)

See Also