{"id":1690,"date":"2026-02-14T04:24:15","date_gmt":"2026-02-14T04:24:15","guid":{"rendered":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/chapter\/ch049-xhtml\/"},"modified":"2026-08-01T09:32:16","modified_gmt":"2026-08-01T09:32:16","slug":"formulating-delayed-yield-analysis","status":"publish","type":"chapter","link":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/chapter\/formulating-delayed-yield-analysis\/","title":{"raw":"10.3 Formulating Delayed Yield Analysis","rendered":"10.3 Formulating Delayed Yield Analysis"},"content":{"raw":"
\r\n
\r\n
\r\n\r\nSection 6 stated that differential equations describing the response of drawdown to pumping in unconfined settings are non\u2011linear and cannot be solved directly because the independent variables (e.g., thickness and transmissivity) are a function of the dependent variable (i.e., head). Analytical approaches describing the response of unconfined groundwater systems have been developed assuming transient drawdown behavior can be approximated by combinations of linear differential equations. Analytical equations developed to represent the delayed yield process are mathematically complex and different representations of the delayed yield analytical solutions yield slightly different values of aquifer properties (e.g., Lin et al., 2019; Nwankwor et al., 1992).\r\n\r\n<\/div>\r\n
\r\n\r\nAs described in Lohman (1972), Boulton (1954a, 1954b, 1963) along with Boulton and others (1964) developed a set of analytical equations and type curves to analyze hydraulic testing data exhibiting delayed yield. These approaches were modified by Neuman (1972, 1974, 1975). He generated a linear differential equation that assumed a line sink for the pumping well and instantaneous drainage from the capillary fringe and unsaturated zone above the water table (Moench et al., 2001).\r\n\r\n<\/div>\r\n
\r\n\r\nStrictly speaking, it is not correct to assume \u201cinstantaneous release of stored water from pores at the water table\u201d. Drainage occurs from pores within the capillary fringe and unsaturated pores above the capillary fringe but this drainage is not instantaneous. As the water table declines, a vertical hydraulic gradient occurs within the zone above the water table including the capillary fringe, especially near the pumped well (Nwankor et al, 1992). Drawdown of the capillary fringe is slower than drawdown of the water table. For this reason, it is best to use drawdown data from observation wells that are several aquifer thicknesses distant from the pumped well where vertical gradients within the saturated zone are very small.\r\n\r\n<\/div>\r\n
\r\n\r\nSince Neuman\u2019s work a number of authors have attempted to incorporate the effect of the release and movement of water in the capillary fringe and unsaturated zone within the cone of depression into delayed yield analytical solutions (e.g., Moench et al, 2001; Nwankwor et al., 1992). Lin and others (2019) provide a good summary of previous efforts and a modified method to assess unconfined hydraulic tests with delayed yield. The equations they derived are solved and linked with field test data using sophisticated numerical methods. Methods developed by Tartakovsky & Neuman (2007) and Moench (1997) are available most commercial software packages (Section 13).\r\n\r\n<\/div>\r\n
\r\n\r\nNeuman\u2019s curve matching methods are presented here to be consistent with other sections of this book. Moench and others (2001) state that for the conditions they evaluated, Neuman curve matching with late time data provided reasonable estimates of T<\/em>, Kh<\/sub><\/em>, Kz<\/sub><\/em>, and Sy<\/sub><\/em> (where Kz<\/sub><\/em> is the vertical hydraulic conductivity of the aquifer). However, their work and the work of others revealed that Neuman\u2019s approach missed the effects of storage\u2011release processes at early times. They argue that other methods to estimate storage properties should be used when relying on early and intermediate drawdowns. Such methods account for the addition of stored water to the water table by unsaturated flow processes (Moench et al., 2001).\r\n\r\n<\/div>\r\n
\r\n\r\nRecognizing these issues, this book presents the curve\u2011matching method described by Neuman (1975) which is found in many groundwater textbooks. It provides tables of variables and plotted families of curves. Readers interested in a more detailed discussion of delayed yield representations will find it in Moench and others (2001) and Lin and others (2019). Commercial software packages typically include the Neuman method as well as other methods that can be used to analyze unconfined pumping test results. Lists of software and the methods they employ are provided in Section 13 of this book.\r\n\r\n<\/div>\r\n
\r\n

Mathematical Development of a Delayed Yield Analysis Method<\/h3>\r\n
\r\n\r\nThe Neuman (1974, 1975) solution involves an integral with three variables, ua<\/sub><\/em>, uy,<\/sub><\/em> and \u03b7<\/em> as shown in Equation\u00a0(70)<\/a>. W(ua<\/sub><\/em>,uy<\/sub><\/em>,\u03b7<\/em>) is the Neuman well function. Schwartz and Zhang, (2003) note that evaluation of Neuman\u2019s analytical solution takes \u201c\u2026a large amount of time\u201d. Strengthening his original work, numerical methods have been applied to increase accuracy and efficiency of computations (Moench & Qgata, 1984; Moench 1993, 1995, 1996). The resulting general equation is described by Equations (70)<\/a>, (71)<\/a>, (72)<\/a>, and (73)<\/a>.\r\n\r\n<\/div>\r\n<\/colgroup>\r\n\r\n\r\n\r\n\r\n\r\n
<\/th>\r\n[latex]s(r,t) = \\frac{Q}{4\\pi T}\\ W\\left( u_{a},u_{y,}\\eta \\right)[\/latex]<\/th>\r\n(70)<\/span><\/th>\r\n<\/tr>\r\n
<\/td>\r\n[latex]u_{a} = \\frac{r^{2}S}{4Tt}[\/latex]<\/td>\r\n(71)<\/span><\/td>\r\n<\/tr>\r\n
<\/td>\r\n[latex]u_{y} = \\frac{r^{2}S_{y}}{4Tt}[\/latex]<\/td>\r\n(72)<\/span><\/td>\r\n<\/tr>\r\n
<\/td>\r\n[latex]\\eta = \\frac{r^{2}K_{z}}{b^{2}K_{h}}[\/latex]<\/td>\r\n(73)<\/span><\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n
\r\n\r\nwhere:\r\n\r\n<\/div>\r\n<\/colgroup>\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n
s<\/em>(r<\/em>,t<\/em>)<\/th>\r\n=<\/th>\r\ndrawdown at a radial distance r<\/em> and time t<\/em> (L)<\/th>\r\n<\/tr>\r\n<\/thead>\r\n
Q<\/em><\/td>\r\n=<\/td>\r\nconstant well discharge rate (L3<\/sup>T\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n
T<\/em><\/td>\r\n=<\/td>\r\ntransmissivity (L2<\/sup>T\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n
S<\/em><\/td>\r\n=<\/td>\r\nstorativity under early time confined conditions (bSs<\/sub><\/em>) (dimensionless)<\/td>\r\n<\/tr>\r\n
b<\/em><\/td>\r\n=<\/td>\r\npre\u2011pumping saturated thickness (L)<\/td>\r\n<\/tr>\r\n
Ss<\/sub><\/em><\/td>\r\n=<\/td>\r\nspecific storage (L\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n
Sy<\/sub><\/em><\/td>\r\n=<\/td>\r\nstorativity at late times, specific yield (dimensionless)<\/td>\r\n<\/tr>\r\n
t<\/em><\/td>\r\n=<\/td>\r\ntime (T)<\/td>\r\n<\/tr>\r\n
r<\/em><\/td>\r\n=<\/td>\r\nradial distance to the observation well (L)<\/td>\r\n<\/tr>\r\n
W (ua<\/sub><\/em>,uy<\/sub><\/em>,\u03b7<\/em>)<\/td>\r\n=<\/td>\r\nNeuman unconfined well function (dimensionless)<\/td>\r\n<\/tr>\r\n
ua<\/sub><\/em><\/td>\r\n=<\/td>\r\nearly time values of u<\/em> (dimensionless)<\/td>\r\n<\/tr>\r\n
uy<\/sub><\/em><\/td>\r\n=<\/td>\r\nlate time values of u<\/em> (dimensionless)<\/td>\r\n<\/tr>\r\n
Kz<\/sub><\/em><\/td>\r\n=<\/td>\r\nunconfined vertical hydraulic conductivity (LT\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n
Kh<\/sub><\/em><\/td>\r\n=<\/td>\r\nunconfined horizontal hydraulic conductivity (LT\u20111<\/sup>)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n
\r\n\r\nComputation of drawdown at an observation well located some distance r<\/em> from the pumping well requires values of T<\/em>, S<\/em>, and Sy<\/sub><\/em>, the initial saturated thickness of the aquifer b<\/em>, and the horizontal and vertical hydraulic conductivity values of the unconfined aquifer as indicated in Equations (70)<\/a> through (73)<\/a>. These parameters are used to generate a value of the well function, W(ua<\/sub><\/em>,uy<\/sub><\/em>,\u03b7<\/em>), from compiled tables (Figure 60<\/a> and Figure 61<\/a>). Computations to predict drawdown using Equation\u00a0(70)<\/a> can then be made for any value or r<\/em> and t<\/em>. Early time drawdown prediction uses ua <\/sub><\/em>of Equation\u00a0(71)<\/a> and late time drawdown prediction uses uy<\/sub><\/em> of Equation\u00a0(72)<\/a>.\r\n\r\n<\/div>\r\n
\r\n\r\n\"\"\r\n\r\n<\/div>\r\n
\r\n
\r\n
\r\n\r\nFigure 60 \u2011<\/strong> Table of well function values for early times presented by Kruseman and de Ridder (2000) for the Neuman well function for unconfined aquifers with delayed yield. They used a slightly different notation than presented in this book, using uA<\/sub><\/em> as equivalent to ua<\/sub><\/em> and \u03b2<\/em> in the place of \u03b7<\/em>. That is, values of \u03b2<\/em> in this table are equivalent to \u03b7<\/em> in this book (after Kruseman & de Ridder, 2000).<\/span>\r\n\r\n<\/div><\/blockquote>\r\n<\/div>\r\n
\r\n\r\n\"\"\r\n\r\n<\/div>\r\n
\r\n
\r\n
\r\n\r\nFigure 61 \u2011<\/strong> Table of well function values for late time presented by Kruseman and de Ridder (2000) for the Neuman well function for unconfined aquifers with delayed yield. They used a slightly different notation than presented in this book, using uB<\/sub><\/em> as equivalent to uy<\/sub><\/em> and \u03b2<\/em> in the place of \u03b7<\/em>. That is, values of \u03b2<\/em> in this table are equivalent to \u03b7<\/em> in this book (after Kruseman & de Ridder, 2000).<\/span>\r\n\r\n<\/div><\/blockquote>\r\n<\/div>\r\n
\r\n\r\nAs stated previously, the table of values (Figure 60<\/a> and Figure 61<\/a>) represent a subset of values that are generated by substituting combinations of ua<\/sub><\/em>, uy <\/sub><\/em>, and \u03b7<\/em> into the analytical solution to produce values of W(ua<\/sub><\/em>,uy<\/sub><\/em>,\u03b7<\/em>). The values in the tables can be plotted to create type curves. If type curves with other values of \u03b7<\/em> are required they can be generated using an available analytical solution (e.g., Neuman, 1972, 1974, 1975). Most often other values of \u03b7<\/em> are interpolated from the tables shown in Figure 60<\/a> and Figure 61<\/a>.\r\n\r\n<\/div>\r\n
\r\n\r\nExample<\/strong>\r\n\r\n<\/div>\r\n
\r\n\r\nIn this example the drawdown at an observation well located 200 m from the pumping well in an unconfined aquifer is computed. The pre\u2011pumping saturated thickness is 30 m, T<\/em> = 3,000 m2<\/sup>\/d, Sy<\/sub><\/em> = 0.15, Kz <\/sub><\/em>= 10 m\/d and Kh<\/sub><\/em> = 100 m\/d. Assume the well has been pumping at a constant rate of 2,505 m3<\/sup>\/d for 1 day and that delayed yield is likely to have ceased. The late time uy<\/sub><\/em> value is used to estimate storativity.\r\n\r\n<\/div>\r\n
\r\n
\r\n

[latex]u_{y} = \\frac{(200\\ \\text{m})^{2} 0.15}{4\\left(3000\\frac{\\text{m}^{2}}{\\text{d}}\\right)(1\\ \\text{d})} = 0.5,\\ \\text{then}\\ \\frac{1}{u_{y}} = 2 = \\frac{1}{u_{B}}\\ \\text{of Figure 61}[\/latex]<\/p>\r\n

[latex]\\eta = \\frac{(200\\ \\text{m})^{2}\\ 10\\frac{\\text{m}}{\\text{d}}}{(30\\ \\text{m})^{2}\\ 100\\frac{\\text{m}}{\\text{d}}} = 4.4 = \\beta\\ \\text{of Figure 61}[\/latex]<\/p>\r\n\r\n<\/div>\r\n<\/div>\r\n

\r\n\r\nThe following values are from Figure 61<\/a>.\r\n\r\n<\/div>\r\n<\/colgroup>\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n
1\/uB<\/sub><\/em> equivalent to 1\/uy<\/sub><\/em> of this book<\/th>\r\n\u03b2\u00a0<\/em>equivalent to \u03b7<\/em> of this book<\/th>\r\nW(uA<\/sub>,uB<\/sub><\/em>,\u03b2<\/em>)\u00a0<\/em>equivalent to W(ua<\/sub><\/em>,uy<\/sub><\/em>,\u03b7<\/em>) of this book<\/th>\r\n<\/tr>\r\n<\/thead>\r\n
1.4<\/td>\r\n4.0<\/td>\r\n0.455<\/td>\r\n<\/tr>\r\n
2.4<\/td>\r\n4.0<\/td>\r\n0.718<\/td>\r\n<\/tr>\r\n
1.4<\/td>\r\n5.0<\/td>\r\n0.430<\/td>\r\n<\/tr>\r\n
2.4<\/td>\r\n5.0<\/td>\r\n0.709<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n
\r\n\r\nUsing the values from Figure 61<\/a> and linear interpolation, first estimate the value of the well function associated with \u03b2<\/em>=4 and 1\/uB<\/sub><\/em>=2, which is W(uA<\/sub><\/em>,uB<\/sub><\/em>,\u03b2<\/em>)=0.6128. Next determine the well function value associated with \u03b2<\/em>=5 and 1\/uB<\/sub><\/em>=2, which is W(uA<\/sub><\/em>,uB<\/sub><\/em>,\u03b2<\/em>)=0.5974. Finally, interpolate between those values for the well function associated with 1\/uB<\/sub><\/em>=2 and \u03b2<\/em>=4.4, which is W(uA<\/sub><\/em>,uB<\/sub><\/em>,\u03b2<\/em>)=0.6066. Then, substituting these values into Equation\u00a0(70)<\/a>, knowing that for uy<\/sub><\/em>=0.5, 1\/uB<\/sub><\/em>=2, \u03b7<\/em>=4.4, and W(ua<\/sub><\/em>,uy<\/sub><\/em>,\u03b7<\/em>)=0.6066, s<\/em>=0.04\u00a0m as show here.\r\n\r\n<\/div>\r\n
\r\n
\r\n\r\n[latex]s = \\frac{2505\\frac{m^{3}}{d}}{4\\ (3.14)\\ 3000\\frac{m^{2}}{d}}\\ (0.6066) = 0.04\\ m[\/latex]\r\n\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>\r\n<\/div>","rendered":"
\n
\n
\n

Section 6 stated that differential equations describing the response of drawdown to pumping in unconfined settings are non\u2011linear and cannot be solved directly because the independent variables (e.g., thickness and transmissivity) are a function of the dependent variable (i.e., head). Analytical approaches describing the response of unconfined groundwater systems have been developed assuming transient drawdown behavior can be approximated by combinations of linear differential equations. Analytical equations developed to represent the delayed yield process are mathematically complex and different representations of the delayed yield analytical solutions yield slightly different values of aquifer properties (e.g., Lin et al., 2019; Nwankwor et al., 1992).<\/p>\n<\/div>\n

\n

As described in Lohman (1972), Boulton (1954a, 1954b, 1963) along with Boulton and others (1964) developed a set of analytical equations and type curves to analyze hydraulic testing data exhibiting delayed yield. These approaches were modified by Neuman (1972, 1974, 1975). He generated a linear differential equation that assumed a line sink for the pumping well and instantaneous drainage from the capillary fringe and unsaturated zone above the water table (Moench et al., 2001).<\/p>\n<\/div>\n

\n

Strictly speaking, it is not correct to assume \u201cinstantaneous release of stored water from pores at the water table\u201d. Drainage occurs from pores within the capillary fringe and unsaturated pores above the capillary fringe but this drainage is not instantaneous. As the water table declines, a vertical hydraulic gradient occurs within the zone above the water table including the capillary fringe, especially near the pumped well (Nwankor et al, 1992). Drawdown of the capillary fringe is slower than drawdown of the water table. For this reason, it is best to use drawdown data from observation wells that are several aquifer thicknesses distant from the pumped well where vertical gradients within the saturated zone are very small.<\/p>\n<\/div>\n

\n

Since Neuman\u2019s work a number of authors have attempted to incorporate the effect of the release and movement of water in the capillary fringe and unsaturated zone within the cone of depression into delayed yield analytical solutions (e.g., Moench et al, 2001; Nwankwor et al., 1992). Lin and others (2019) provide a good summary of previous efforts and a modified method to assess unconfined hydraulic tests with delayed yield. The equations they derived are solved and linked with field test data using sophisticated numerical methods. Methods developed by Tartakovsky & Neuman (2007) and Moench (1997) are available most commercial software packages (Section 13).<\/p>\n<\/div>\n

\n

Neuman\u2019s curve matching methods are presented here to be consistent with other sections of this book. Moench and others (2001) state that for the conditions they evaluated, Neuman curve matching with late time data provided reasonable estimates of T<\/em>, Kh<\/sub><\/em>, Kz<\/sub><\/em>, and Sy<\/sub><\/em> (where Kz<\/sub><\/em> is the vertical hydraulic conductivity of the aquifer). However, their work and the work of others revealed that Neuman\u2019s approach missed the effects of storage\u2011release processes at early times. They argue that other methods to estimate storage properties should be used when relying on early and intermediate drawdowns. Such methods account for the addition of stored water to the water table by unsaturated flow processes (Moench et al., 2001).<\/p>\n<\/div>\n

\n

Recognizing these issues, this book presents the curve\u2011matching method described by Neuman (1975) which is found in many groundwater textbooks. It provides tables of variables and plotted families of curves. Readers interested in a more detailed discussion of delayed yield representations will find it in Moench and others (2001) and Lin and others (2019). Commercial software packages typically include the Neuman method as well as other methods that can be used to analyze unconfined pumping test results. Lists of software and the methods they employ are provided in Section 13 of this book.<\/p>\n<\/div>\n

\n

Mathematical Development of a Delayed Yield Analysis Method<\/h3>\n
\n

The Neuman (1974, 1975) solution involves an integral with three variables, ua<\/sub><\/em>, uy,<\/sub><\/em> and \u03b7<\/em> as shown in Equation\u00a0(70)<\/a>. W(ua<\/sub><\/em>,uy<\/sub><\/em>,\u03b7<\/em>) is the Neuman well function. Schwartz and Zhang, (2003) note that evaluation of Neuman\u2019s analytical solution takes \u201c\u2026a large amount of time\u201d. Strengthening his original work, numerical methods have been applied to increase accuracy and efficiency of computations (Moench & Qgata, 1984; Moench 1993, 1995, 1996). The resulting general equation is described by Equations (70)<\/a>, (71)<\/a>, (72)<\/a>, and (73)<\/a>.<\/p>\n<\/div>\n\n\n\n\n<\/colgroup>\n\n\n\n\n\n
<\/th>\n[latex]s(r,t) = \\frac{Q}{4\\pi T}\\ W\\left( u_{a},u_{y,}\\eta \\right)[\/latex]<\/th>\n(70)<\/span><\/th>\n<\/tr>\n
<\/td>\n[latex]u_{a} = \\frac{r^{2}S}{4Tt}[\/latex]<\/td>\n(71)<\/span><\/td>\n<\/tr>\n
<\/td>\n[latex]u_{y} = \\frac{r^{2}S_{y}}{4Tt}[\/latex]<\/td>\n(72)<\/span><\/td>\n<\/tr>\n
<\/td>\n[latex]\\eta = \\frac{r^{2}K_{z}}{b^{2}K_{h}}[\/latex]<\/td>\n(73)<\/span><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n
\n

where:<\/p>\n<\/div>\n\n\n\n\n<\/colgroup>\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n
s<\/em>(r<\/em>,t<\/em>)<\/th>\n=<\/th>\ndrawdown at a radial distance r<\/em> and time t<\/em> (L)<\/th>\n<\/tr>\n<\/thead>\n
Q<\/em><\/td>\n=<\/td>\nconstant well discharge rate (L3<\/sup>T\u20111<\/sup>)<\/td>\n<\/tr>\n
T<\/em><\/td>\n=<\/td>\ntransmissivity (L2<\/sup>T\u20111<\/sup>)<\/td>\n<\/tr>\n
S<\/em><\/td>\n=<\/td>\nstorativity under early time confined conditions (bSs<\/sub><\/em>) (dimensionless)<\/td>\n<\/tr>\n
b<\/em><\/td>\n=<\/td>\npre\u2011pumping saturated thickness (L)<\/td>\n<\/tr>\n
Ss<\/sub><\/em><\/td>\n=<\/td>\nspecific storage (L\u20111<\/sup>)<\/td>\n<\/tr>\n
Sy<\/sub><\/em><\/td>\n=<\/td>\nstorativity at late times, specific yield (dimensionless)<\/td>\n<\/tr>\n
t<\/em><\/td>\n=<\/td>\ntime (T)<\/td>\n<\/tr>\n
r<\/em><\/td>\n=<\/td>\nradial distance to the observation well (L)<\/td>\n<\/tr>\n
W (ua<\/sub><\/em>,uy<\/sub><\/em>,\u03b7<\/em>)<\/td>\n=<\/td>\nNeuman unconfined well function (dimensionless)<\/td>\n<\/tr>\n
ua<\/sub><\/em><\/td>\n=<\/td>\nearly time values of u<\/em> (dimensionless)<\/td>\n<\/tr>\n
uy<\/sub><\/em><\/td>\n=<\/td>\nlate time values of u<\/em> (dimensionless)<\/td>\n<\/tr>\n
Kz<\/sub><\/em><\/td>\n=<\/td>\nunconfined vertical hydraulic conductivity (LT\u20111<\/sup>)<\/td>\n<\/tr>\n
Kh<\/sub><\/em><\/td>\n=<\/td>\nunconfined horizontal hydraulic conductivity (LT\u20111<\/sup>)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n
\n

Computation of drawdown at an observation well located some distance r<\/em> from the pumping well requires values of T<\/em>, S<\/em>, and Sy<\/sub><\/em>, the initial saturated thickness of the aquifer b<\/em>, and the horizontal and vertical hydraulic conductivity values of the unconfined aquifer as indicated in Equations (70)<\/a> through (73)<\/a>. These parameters are used to generate a value of the well function, W(ua<\/sub><\/em>,uy<\/sub><\/em>,\u03b7<\/em>), from compiled tables (Figure 60<\/a> and Figure 61<\/a>). Computations to predict drawdown using Equation\u00a0(70)<\/a> can then be made for any value or r<\/em> and t<\/em>. Early time drawdown prediction uses ua <\/sub><\/em>of Equation\u00a0(71)<\/a> and late time drawdown prediction uses uy<\/sub><\/em> of Equation\u00a0(72)<\/a>.<\/p>\n<\/div>\n

\n

\"\"<\/p>\n<\/div>\n

\n
\n
\n

Figure 60 \u2011<\/strong> Table of well function values for early times presented by Kruseman and de Ridder (2000) for the Neuman well function for unconfined aquifers with delayed yield. They used a slightly different notation than presented in this book, using uA<\/sub><\/em> as equivalent to ua<\/sub><\/em> and \u03b2<\/em> in the place of \u03b7<\/em>. That is, values of \u03b2<\/em> in this table are equivalent to \u03b7<\/em> in this book (after Kruseman & de Ridder, 2000).<\/span><\/p>\n<\/div>\n<\/blockquote>\n<\/div>\n

\n

\"\"<\/p>\n<\/div>\n

\n
\n
\n

Figure 61 \u2011<\/strong> Table of well function values for late time presented by Kruseman and de Ridder (2000) for the Neuman well function for unconfined aquifers with delayed yield. They used a slightly different notation than presented in this book, using uB<\/sub><\/em> as equivalent to uy<\/sub><\/em> and \u03b2<\/em> in the place of \u03b7<\/em>. That is, values of \u03b2<\/em> in this table are equivalent to \u03b7<\/em> in this book (after Kruseman & de Ridder, 2000).<\/span><\/p>\n<\/div>\n<\/blockquote>\n<\/div>\n

\n

As stated previously, the table of values (Figure 60<\/a> and Figure 61<\/a>) represent a subset of values that are generated by substituting combinations of ua<\/sub><\/em>, uy <\/sub><\/em>, and \u03b7<\/em> into the analytical solution to produce values of W(ua<\/sub><\/em>,uy<\/sub><\/em>,\u03b7<\/em>). The values in the tables can be plotted to create type curves. If type curves with other values of \u03b7<\/em> are required they can be generated using an available analytical solution (e.g., Neuman, 1972, 1974, 1975). Most often other values of \u03b7<\/em> are interpolated from the tables shown in Figure 60<\/a> and Figure 61<\/a>.<\/p>\n<\/div>\n

\n

Example<\/strong><\/p>\n<\/div>\n

\n

In this example the drawdown at an observation well located 200 m from the pumping well in an unconfined aquifer is computed. The pre\u2011pumping saturated thickness is 30 m, T<\/em> = 3,000 m2<\/sup>\/d, Sy<\/sub><\/em> = 0.15, Kz <\/sub><\/em>= 10 m\/d and Kh<\/sub><\/em> = 100 m\/d. Assume the well has been pumping at a constant rate of 2,505 m3<\/sup>\/d for 1 day and that delayed yield is likely to have ceased. The late time uy<\/sub><\/em> value is used to estimate storativity.<\/p>\n<\/div>\n

\n
\n

[latex]u_{y} = \\frac{(200\\ \\text{m})^{2} 0.15}{4\\left(3000\\frac{\\text{m}^{2}}{\\text{d}}\\right)(1\\ \\text{d})} = 0.5,\\ \\text{then}\\ \\frac{1}{u_{y}} = 2 = \\frac{1}{u_{B}}\\ \\text{of Figure 61}[\/latex]<\/p>\n

[latex]\\eta = \\frac{(200\\ \\text{m})^{2}\\ 10\\frac{\\text{m}}{\\text{d}}}{(30\\ \\text{m})^{2}\\ 100\\frac{\\text{m}}{\\text{d}}} = 4.4 = \\beta\\ \\text{of Figure 61}[\/latex]<\/p>\n<\/div>\n<\/div>\n

\n

The following values are from Figure 61<\/a>.<\/p>\n<\/div>\n\n\n\n\n<\/colgroup>\n\n\n\n\n\n\n\n
1\/uB<\/sub><\/em> equivalent to 1\/uy<\/sub><\/em> of this book<\/th>\n\u03b2\u00a0<\/em>equivalent to \u03b7<\/em> of this book<\/th>\nW(uA<\/sub>,uB<\/sub><\/em>,\u03b2<\/em>)\u00a0<\/em>equivalent to W(ua<\/sub><\/em>,uy<\/sub><\/em>,\u03b7<\/em>) of this book<\/th>\n<\/tr>\n<\/thead>\n
1.4<\/td>\n4.0<\/td>\n0.455<\/td>\n<\/tr>\n
2.4<\/td>\n4.0<\/td>\n0.718<\/td>\n<\/tr>\n
1.4<\/td>\n5.0<\/td>\n0.430<\/td>\n<\/tr>\n
2.4<\/td>\n5.0<\/td>\n0.709<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n
\n

Using the values from Figure 61<\/a> and linear interpolation, first estimate the value of the well function associated with \u03b2<\/em>=4 and 1\/uB<\/sub><\/em>=2, which is W(uA<\/sub><\/em>,uB<\/sub><\/em>,\u03b2<\/em>)=0.6128. Next determine the well function value associated with \u03b2<\/em>=5 and 1\/uB<\/sub><\/em>=2, which is W(uA<\/sub><\/em>,uB<\/sub><\/em>,\u03b2<\/em>)=0.5974. Finally, interpolate between those values for the well function associated with 1\/uB<\/sub><\/em>=2 and \u03b2<\/em>=4.4, which is W(uA<\/sub><\/em>,uB<\/sub><\/em>,\u03b2<\/em>)=0.6066. Then, substituting these values into Equation\u00a0(70)<\/a>, knowing that for uy<\/sub><\/em>=0.5, 1\/uB<\/sub><\/em>=2, \u03b7<\/em>=4.4, and W(ua<\/sub><\/em>,uy<\/sub><\/em>,\u03b7<\/em>)=0.6066, s<\/em>=0.04\u00a0m as show here.<\/p>\n<\/div>\n

\n
\n

[latex]s = \\frac{2505\\frac{m^{3}}{d}}{4\\ (3.14)\\ 3000\\frac{m^{2}}{d}}\\ (0.6066) = 0.04\\ m[\/latex]<\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"author":8,"menu_order":42,"template":"","meta":{"pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":[],"pb_section_license":""},"chapter-type":[],"contributor":[],"license":[],"class_list":["post-1690","chapter","type-chapter","status-publish","hentry"],"part":576,"_links":{"self":[{"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/chapters\/1690","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/wp\/v2\/users\/8"}],"version-history":[{"count":4,"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/chapters\/1690\/revisions"}],"predecessor-version":[{"id":3302,"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/chapters\/1690\/revisions\/3302"}],"part":[{"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/parts\/576"}],"metadata":[{"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/chapters\/1690\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/wp\/v2\/media?parent=1690"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/pressbooks\/v2\/chapter-type?post=1690"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/wp\/v2\/contributor?post=1690"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/books.gw-project.org\/an-introduction-to-hydraulic-testing-in-hydrogeology\/wp-json\/wp\/v2\/license?post=1690"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}