Instructions for running the analytical code PAT (Purge Analyzer Tool) for computation of in-well time of travel of groundwater under pumping conditions

Abstract

Introduction Understanding the optimal time needed to purge a well while pumping to collect a representative groundwater sample requires an understanding of groundwater flow in wells (in-well flow). Parameters that affect in-well flow include the hydraulic properties of the aquifer, well construction, drawdown from pumping, and pump rate. The time of travel relative to in-well flow is affected by the pump’s intake location. The Purge Analyzer Tool (PAT) incorporates hydraulic calculations to help assess the optimal purge times required to vertically transport groundwater in the well to the pump intake (Harte, 2017). Harte (2017) includes a discussion on the rationale for determining in-well groundwater flow and time of travel and also discusses the limitations inherent in the PAT; an understanding of the limitations is important to ensure proper use. The PAT calculates flow by use of the Dupuit-Theim equation (Lohman, 1979) that assumes steady-state radial flow and a total inflow from the well opening or screen equal to the pumping rate (eq. 1). A bulk average hydraulic conductivity (K avg ) is derived from this relationship. Once K avg is calculated, the program calculates incremental (layered) horizontal radial inflow into the well over user defined increments (layers). These defined increments represent the screen or well opening as a fraction of the total inflow. The amount of inflow per layer is proportional to the user-defined layered distribution of hydraulic conductivity (K layer ) because drawdown is assumed to be uniformly distributed in the well. The water budget equation that guides the solution of the PAT (eq. 1) is specified as: Q p = Q v + Q H + Q w (1) where Q P is pumping rate, Q v is vertical flow entering the boundary of the mixing zone (M z ) from the summation of layered radial flow (∑ Q hl-n ) where l-n denotes number of layers, Q H is horizontal radial flow into the mixing zone (M z ), and Q w is flow from wellbore storage effects. The in-well flow is computed from the convergence of incremental (layered) radial inflows (Q hl-n ) summed to the total vertical flow (Q V ) entering the adjacent zone to the pump intake (called mixing zone [M z ]) as shown in figure 1. The Q v is transported as one-dimensional piston flow. Within the M z , it's assumed that flow to the pump is dominated by horizontal radial flow (Q H ) when the pump is in the open interval of the well. Flow from the wellbore storage (Q w ) is computed from the volume of water pumped from the well at the time of the drawdown (s) measurement(s). Aquifer storage effects are unaccounted for but are likely to be problematic when (1) dewatering within the well opening occurs or (2) when the water table is close to the top of the well screen or open interval where additional flow into the upper portion of the well opening may occur. For fully saturated wells tens of feet below the water table, storage effects are likely to be more uniformly distributed across the well screen or open interval (regardless of confined or unconfined conditions). Therefore, radial inflow from storage will be less prominent under pump rates commonly used in groundwater sampling either for volumetric sampling ( < 3 gallons per minute) or low-flow sampling ( < 0.5 liters per minute). A major benefit of the use of the PAT is the understanding of time-varying, vertical integration of captured pump water. The analytical model computes aquifer (formation) capture intervals relative to the open interval of the well. This information is displayed graphically (called aquifer fraction graphs) and can be used to assess the likely formation intervals contributing water to the sample at any time.