MODELING GROUND WATER FLOW USING FLUX BOUNDARY CONDITIONS1

ABSTRACT: Determination of the boundary conditions for modeling ground water flow is a critical point especially in regional models. Normally the regional models require model areas that are greater than the given area of interest. This work focuses on the prediction of hydraulic heads in regional models using flux boundary conditions. The model uses flux boundary conditions that were estimated using a radial flow analog and Darcy's law. The regional model that is presented uses no parameter identification (inverse estimation) procedures. In the present work, the Houston area was used. The simulation of the hydrological conditions of the Chicot and Evangeline Aquifers that underlie the Houston area were made using the available information about the geological profile in the Houston region and the current information about the existing production wells. The regional model works as a forward problem. The system parameters such as hydraulic conductivity, specific storage, and hydrological stresses were specified, and the model predicts the hydraulic head. Actual data from piezometers operated by the U.S. Geological Survey (USGS) in many places throughout Houston were used as initial conditions. Some piezometric head data were generated using the regional variable theory called kriging to supply head estimates in areas where data were unavailable. The Modular Three Dimensional Finite Difference Groundwater Flow Model developed by the USGS was used to predict the hydraulic heads. The predicted ground water heads are compared to the actual data. The results show that the model performs well for locations where data were available.     

EVALUATION OF METHODS FOR ESTIMATING STREAM WATER QUALITY PARAMETERS IN A TRANSIENT MODEL FROM STOCHASTIC DATA1

ABSTRACT: In this study the estimation of parameters in water quality models represented by linear first order partial differential equations is investigated. Two sets of simulated input-output data, one with input noise and the other with output measurement error, were used. The parameters were estimated by a gradient technique (Bard's method) and a pattern search technique. The results indicate that the output measurement error significantly affects the values of parameter estimates as compared to the noise added to the input. Bard's method consistently gave results with a smaller sum of square value.     

GROUND WATER YIELDS AND INVENTORY FOR LAND-USE PLANNING IN CRYSTALLINE ROCK AREAS OF WAKE COUNTY,NORTH CAROLINA1

In urbanizing areas underlain by crystalline rocks an inventory of the usable ground water should form an important element in the land-use planning process. Land-use decisions are most often made upon water-well yield data alone, but these data do not address the inventory question. A method utilizing local geology, stream lowflow characteristics, and pumping test data permitted outlining portions of Wake County, North Carolina, which are more favorable and those which are less favorable for ground water supply development. Information from state-required 24-hour pumping tests on 232 wells was supplemented by information from an additional 100 wells whose initial yields were determined by shorter pumping tests. Comparison of the ground water inventory values on a per acre or per square mile basis with average water use at various residential densities provides information useful to the responsible political body as it decides about population density distribution and the need for surface water supplies. For Wake Country an average density of one residential unit per acre appears generally permissible before overdraft of the ground water supplies becomes a significant risk. By utilizing the ground water inventory and patterns of ground water yields from various rock types, county officials can maximize the effectiveness of public facilities funds. The technique appears useful for other parts of the Piedmont region of the southeastern United States.     

MEASUREMENT OF HYDRAULIC CONDUCTIVITY AND DIFFUSIVITY FOR PREDICTING THE PROCESS OF SOIL WATER INFILTRATION FROM A TRICKLE SOURCE1

ABSTRACT: Two soil water functions, hydraulic conductivity K(θ) and diffusivity D(θ), were estimated by two methods In one method D(θ) was estimated according to Bruce and Klute (1956), and K(θ) was calculated from D(θ) and the retention curve. In the second, K(θ) was obtained by field estimation, with D(θ) being calculated from K(θ) and the retention curve. The criterion of reliability for both methods was agreement between experimental and predicted distribution of soil water content. The prediction was made using the functions K(θ) and D(θ) as soil water parameters in both methods. Theoretical and experimental agreement was generally good. The first method, however, was found to be best for high soil water content and the second for low soil water content. In addition, the water content at the end of the monotonic increase of function D(θ) (estimated according to Bruce and Klute 1956) was found to be about the upper limit of field soil water content. It can be used as a boundary condition in the numerical solution of a cylindrical model of infiltration from a trickle source. It was concluded that the best agreement between theory and experiment can be found when the combined values of D(θ) and K(θ) from both methods of estimation are used.     

HISTORICAL AND ESTIMATED GROUND WATER LEVELS NEAR WINNIPEG,CANADA, AND THEIR SENSITIVITY TO CLIMATIC VARIABILITY1

ABSTRACT: Long term well hydrographs and estimated ground water levels derived from hydroclimatic and biological data were used to evaluate trends within the Upper Carbonate Aquifer (UCA) near Winnipeg, Canada, during the 20th Century. Ground water records from instruments have been kept since the early 1960s and are derived from piezometers in the overlying sediments and in open boreholes in the UCA. Some boreholes extend into an underlying Paleozoic carbonate sequence. Shallow well hydrographs show no obvious long term trends but do exhibit variations on the order of three to four years that are correlated with changes in annual temperature and precipitation at lags up to 24 months. Trends observed in deeper wells appear to be largely related to ground water usage patterns and show little correlation with climate over the past 35 years. Stepwise multiple regression modeled average annual hydraulic head in the shallow wells as a function of regional temperature, precipitation, and tree ring variables. Estimated hydraulic heads had a slightly greater range prior to the 1960s, most prominently during an interval of lowered ground water levels between 1930 and 1942. Regression results demonstrate that moisture sensitive tree ring data are viable predictors of past ground water levels and may be useful for studies of aquifers in regions that lack long, high quality precipitation records.     

Effects of calibration on L-THIA GIS runoff and pollutant estimation

Urbanization can result in alteration of a watershed's hydrologic response and water quality. To simulate hydrologic and water quality impacts of land use changes, the Long-Term Hydrologic Impact Assessment (L-THIA) system has been used. The L-THIA system estimates pollutant loading based on direct runoff quantity and land use based pollutant coefficients. The accurate estimation of direct runoff is important in assessing water quality impacts of land use changes. An automated program was developed to calibrate the L-THIA model using the millions of curve number (CN) combinations associated with land uses and hydrologic soil groups. L-THIA calibration for the Little Eagle Creek (LEC) watershed near Indianapolis, Indiana was performed using land use data for 1991 and daily rainfall data for six months of 1991 (January 1-June 30) to minimize errors associated with use of different temporal land use data and rainfall data. For the calibration period, the Nash-Sutcliffe coefficient was 0.60 for estimated and observed direct runoff. The calibrated CN values were used for validation of the model for the same year (July 1-December 31), and the Nash-Sutcliffe coefficient was 0.60 for estimated and observed direct runoff. The Nash-Sutcliffe coefficient was 0.52 for January 1, 1991 to December 31, 1991 using uncalibrated CN values. As shown in this study, the use of better input parameters for the L-THIA model can improve accuracy. The effects on direct runoff and pollutant estimation of the calibrated CN values in the L-THIA model were investigated for the LEC. Following calibration, the estimated average annual direct runoff for the LEC watershed increased by 34%, total nitrogen by 24%, total phosphorus by 22%, and total lead by 43%. This study demonstrates that the L-THIA model should be calibrated and validated prior to application in a particular watershed to more accurately assess the effects of land use changes on hydrology and water quality.     

USING SENSITIVITY ANALYSIS TO ASSIST PARAMETER ZONATION IN GROUND WATER FLOW MODEL1

ABSTRACT: For numerical modeling of ground water movement in a real aquifer system, the aquifer is usually divided into hydrogeologically defined zones, each with its own parameter values. The responses of the system, such as head or drawdown, are often available only in some of the zones. The estimated parameters of all the zones are based on the measured response in these limited zones. However, the estimates for some of the zones may be very uncertain, and these zones are therefore not justified by the data. In this paper, an approach is presented to understand which zone may produce uncertain parameter values and should be lumped with its neighbor. This approach is demonstrated using a regional numerical model for pumping test analysis in the Nottinghamshire aquifer, UK. A step-by-step process is used in identifying the aquifer zones and estimating their parameters based on the principle of using the smallest possible numbers of zones and parameters for adequate representation of the drawdown response. After the parameters of each zone are estimated, the sensitivity features of these parameters are examined. The results show that the parameters in one zone can be estimated properly by the drawdown in another zone only when there is significant sensitivity. For transmissivity, sensitivity between zones occurs when there is significant flow between them. For storativity, sufficient sensitivity can occur without large flows between the zones, provided that one zone causes significant drawdown in the other. This idea can be extended to the flow model for a large aquifer system. If the aquifer is divided in such a way that aquifer responses are not sensitive to the parameters in some of the zones, the parameters in those zones cannot be estimated properly and should be lumped into their neighboring zones. In this way, a simple but more reasonable model can be built.     

DYNAMIC MODELING OF STREAM QUALITY BY INVARIANT IMBEDDING1

ABSTRACT. The estimator equations obtained using invariant imbedding is used to estimate the parameters in river or stream pollution. By using these equations, the parameters can be estimated directly from differential equations representing the pollution model and from measured noisy data such as BOD and DO. Another advantage of this approach is that a sequential estimation scheme is obtained. By using this sequential scheme, only current data are needed to estimate current or future values of the unknown parameters. Consequently, a large amount of computer time and computer memory can be saved. Furthermore, not only the parameters but also the concentrations of pollutants can be estimated. Thus, it also forms an effective forecasting technique. The classical least squares criterion is used in the estimation. Several examples are solved to illustrate the technique. (KEY WORDS: dynamic modeling; water pollution; invariant imbedding; forecasting; least squares criterion; estimation)     

THE FLEXIBLE TOLERANCE METHOD FOR ESTIMATING HYDROLOGIC PARAMETERS IN THE ROOT ZONE1

In this study, a constrained minimization method, the flexible tolerance method, was used to solve the optimization problems for determining hydrologic parameters in the root zone: water uptake rate, spatial root distribution, infiltration rate, and evaporation. Synthetic soil moisture data were first generated using the Richards' equation and its associated initial and boundary conditions, and these data were then used for the inverse analyses. The results of inverse simulation indicate the following. If the soil moisture data contain no noise, the rate of estimated water uptake and spatial root distribution parameters are equal to the true values without using constraints. If there is noise in the observed data, constraints must be used to improve the quality of the estimate results. In the estimation of rainfall infiltration and surface evaporation, interpolation methods should be used to reduce the number of unknowns. A fewer number of variables can improve the quality of inversely estimated parameters. Simultaneous estimation of spatial root distribution and water uptake rate or estimation of evaporation and water uptake rate is possible. The method was used to estimate the water uptake rate, spatial root distribution, infiltration rate, and evaporation using long‐term soil moisture data collected from Nebraska's Sand Hills.     

Seawater Intrusion and Sustainable Yield of Basal Aquifers1

Liu, Clark C.K. and John J. Dai, 2012. Seawater Intrusion and Sustainable Yield of Basal Aquifers. Journal of the American Water Resources Association (JAWRA) 48(5): 861‐870. DOI: 10.1111/j.1752‐1688.2012.00659.x Abstract: Basal aquifers, in which freshwater floats on top of saltwater, are the major freshwater supply for the Hawaiian Islands, as well as many other coastal regions around the world. Under unexploited or natural conditions, freshwater and the underlying seawater are separated by a relatively sharp interface located below mean sea level at a depth of about 40 times the hydraulic head. With forced draft, the hydraulic head of a basal aquifer would decline and the sharp interface would move up. It is a serious problem of seawater intrusion as huge amounts of freshwater storage is replaced by saltwater. Also, with forced draft, the sharp interface is replaced by a transition zone in which the salinity increases downward from freshwater to saltwater. As pumping continues, the transition zone expands. The desirable source‐water salinity in Hawaii is about 2% of the seawater salinity. Therefore, the transition zone expansion is another serious problem of seawater intrusion. In this study, a robust analytical groundwater flow and salinity transport model (RAM2) was developed. RAM2 has a simple mathematical structure and its model parameters can be determined satisfactorily with the available field monitoring data. The usefulness of RAM2 as a viable management tool for coastal ground water management is demonstrated by applying it to determine the sustainable yield of the Pearl Harbor aquifer, a principal water supply source in Hawaii.     

ANALYSIS OF PUMPING TESTS: SIGNIFICANCE OF WELL DIAMETER,PARTIAL PENETRATION,AND NOISE1

ABSTRACT: The nonlinear least squares (NLS) method was applied to pumping and recovery aquifer test data in confined and unconfined aquifers with finite diameter and partially penetrating pumping wells, and with partially penetrating piezometers or observation wells. It was demonstrated that noiseless and moderately noisy drawdown data from observation points located less than two saturated thicknesses of the aquifer from the pumping well produced an exact or acceptable set of parameters when the diameter of the pumping well was included in the analysis. The accuracy of the estimated parameters, particularly that of specific storage, decreased with increases in the noise level in the observed drawdown data. With consideration of the well radii, the noiseless drawdown data from the pumping well in an unconfined aquifer produced good estimates of horizontal and vertical hydraulic conductivities and specific yield, but the estimated specific storage was unacceptable. When noisy data from the pumping well were used, an acceptable set of parameters was not obtained. Further experiments with noisy drawdown data in an unconfined aquifer revealed that when the well diameter was included in the analysis, hydraulic conductivity, specific yield and vertical hydraulic conductivity may be estimated rather effectively from piezometers located over a range of distances from the pumping well. Estimation of specific storage became less reliable for piezometers located at distances greater than the initial saturated thickness of the aquifer. Application of the NLS to field pumping and recovery data from a confined aquifer showed that the estimated parameters from the two tests were in good agreement only when the well diameter was included in the analysis. Without consideration of well radii, the estimated values of hydraulic conductivity from the pumping and recovery tests were off by a factor of four.     

GEOSTATISTICAL ESTIMATION OF HORIZONTAL HYDRAULIC CONDUCTIVITY FOR THE KIRKWOOD‐COHANSEY AQUIFER1

ABSTRACT: The Kirkwood‐Cohansey aquifer has been identified as a critical source for meeting existing and expected water supply needs for southern New Jersey. Several contaminated sites exist in the region; their impact on the aquifer has to be evaluated using ground water flow and transport models. Ground water modeling depends on availability of measured hydrogeologic data (e.g., hydraulic conductivity, for parameterization of the modeling runs). However, field measurements of such critical data have inadequate spatial density, and their locations are often clustered. The goal of this study was to research, compile, and geocode existing data, then use geostatistics and advanced mapping methods to develop a map of horizontal hydraulic conductivity for the Kirkwood‐Cohansey aquifer. Spatial interpolation of horizontal hydraulic conductivity measurements was performed using the Bayesian Maximum Entropy (BME) Method implemented in the BMELib code library. This involved the integration of actual measurements with soft information on likely ranges of hydraulic conductivity at a given location to obtain estimate maps. The estimation error variance maps provide an insight into the uncertainty associated with the estimates, and indicate areas where more information on hydraulic conductivity is required.     

A LUMPED APPROACH TO THE INVERSE PROBLEM IN GROUND WATER HYDROLOGY1

ABSTRACT: A solution procedure to solve the inverse problem in ground water, based on lumped approach, has been proposed. The method has the following advantages: 1) exact determination of the boundary conditions and the physical laws of flow through porous media is not required; 2) all errors of approximation in describing the boundary conditions, physical laws, and the aquifer properties are lumped into the surrogate parameters; and 3) the same mathematical model can be employed both in the identification process and in the subsequent management studies. The optimal values of the surrogate parameters are found by using a multidimensional unconstrained optimization code devised by Powell. The solution procedure and the convergence characteristics of the proposed algorithm have been illustrated by two hypothetical problems.     

Estimating Hydraulic Properties of Volcanic Aquifers Using Constant‐Rate and Variable‐Rate Aquifer Tests1

Abstract: In recent years the ground‐water demand of the population of the island of Maui, Hawaii, has significantly increased. To ensure prudent management of the ground‐water resources, an improved understanding of ground‐water flow systems is needed. At present, large‐scale estimations of aquifer properties are lacking for Maui. Seven analytical methods using constant‐rate and variable‐rate withdrawals for single wells provide an estimate of hydraulic conductivity and transmissivity for 103 wells in central Maui. Methods based on constant‐rate tests, although not widely used on Maui, offer reasonable estimates. Step‐drawdown tests, which are more abundantly used than other tests, provide similar estimates as constant‐rate tests. A numerical model validates the suitability of analytical solutions for step‐drawdown tests and additionally provides an estimate of storage parameters. The results show that hydraulic conductivity is log‐normally distributed and that for dike‐free volcanic rocks it ranges over several orders of magnitude from 1 to 2,500 m/d. The arithmetic mean, geometric mean, and median values of hydraulic conductivity are respectively 520, 280, and 370 m/d for basalt and 80, 50, and 30 m/d for sediment. A geostatistical approach using ordinary kriging yields a prediction of hydraulic conductivity on a larger scale. Overall, the results are in agreement with values published for other Hawaiian islands.     

ROOT DEPTH AS A SENSITIVE PARAMETER IN A DETERMINISTIC HYDROLOGIC MODEL1

ABSTRACT: Deterministic models of watershed hydrology require accurate a priori estimates of soil, vegetation, and watershed parameters. Physical fidelity of these values to those of the prototype natural watershed is essential. One vegetation parameter most neglected, perhaps because it is least understood, is plant root activity. Plant roots directly or indirectly affect many hydrologic processes, including evaporation, transpiration, soil moisture, and ground water. One of their more important functions is in opening surface-connected hydraulic pathways for rainfall penetration. This paper presents the results of a study in which available information on roots has been applied in hydrologic computations.     

Pesticide and Transformation Product Detections and Age‐Dating Relations from Till and Sand Deposits1

Abstract:  Pesticide and transformation product concentrations and frequencies in ground water from areas of similar crop and pesticide applications may vary substantially with differing lithologies. Pesticide analysis data for atrazine, metolachlor, alachlor, acetochlor, and cyanazine and their pesticide transformation products were collected at 69 monitoring wells in Illinois and northern Indiana to document occurrence of pesticides and their transformation products in two agricultural areas of differing lithologies, till, and sand. The till is primarily tile drained and has preferential fractured flow, whereas the sand primarily has surface water drainage and primary porosity flow. Transformation products represent most of the agricultural pesticides in ground water regardless of aquifer material – till or sand. Transformation products were detected more frequently than parent pesticides in both the till and sand, with metolachlor ethane sulfonic acid being most frequently detected. Estimated ground‐water recharge dates for the sand were based on chlorofluorocarbon analyses. These age‐dating data indicate that ground water recharged prior to 1990 is more likely to have a detection of a pesticide or pesticide transformation product. Detections were twice as frequent in ground water recharged prior to 1990 (82%) than in ground water recharged on or after 1990 (33%). The highest concentrations of atrazine, alachlor, metolachlor, and their transformation products, also were detected in samples from ground water recharged prior to 1990. These age/pesticide detection relations are opposite of what would normally be expected, and may be the result of preferential flow and/or ground‐water mixing between aquifers and aquitards as evident by the detection of acetochlor transformation products in samples with estimated ground‐water ages predating initial pesticide application.     

APPLICATION OF TABU SEARCH TO GROUND WATER PARAMETER ZONATION1

ABSTRACT: The spatial distribution of hydrogeologic parameter is an important issue in ground water simulation. One of the methods is to divide an area into several zones such that parameters are assumed to be constant within zone. The purpose of this study is to apply Tabu Search (TS) to find the best zonation of parameters that can result in the best ground water simulation. The initial zonation can be determined as the Thiessen method, and then zonation is optimized by T.S. The mean square error between simulated and observed hydraulic heads was used as the objective function. A designed confined aquifer with known zonation was used as an example to test the proposed method. Results indicated that Tabu Search can locate the optimal zonation successfully and avoid being trapped by local optimal zonations. Besides this, four other arbitrary initial zonations can be directed to the optimal zonation by TS, which proves the robustness of the proposed method. The method proposed in this study is feasible and expected to work well in the field problems with sufficient sampling of concerned parameters.     

OPTIMAL ESTIMATION OF CONTAMINANT TRANSPORT IN GROUND WATER1

ABSTRACT: Finite element and finite difference representations of the convective-dispersive equation have been widely used in determining contaminant transport in ground water. Due to inherent uncertainties of the transport process, those representations are inexact and contain errors. Errors in field measurements are unavoidable. By combining a numerical model, a measurement equation, and the Kalman filter, optimal estimates of the state variable (contaminant concentration) can be obtained. This paper describes the algorithm and gives a numerical example of contaminant transport in a two-dimensional ground water flow. The results show significant improvement in the estimated concentration distribution by using the filtering technique.     

油田井下作业清洁生产配套工艺

近几年来,江苏油田加大油田井下作业清洁生产配套工艺技术的研究力度,自行研制了地膜隔离技术以及二合一液压锚定泄油器、井口油管刮油器、抽油杆刮蜡器、抽汲环保井口装置、移动式作业废水集中处理装置等一系列新型综合防治落地油和控制作业污水排放的清洁生产配套工艺技术。这些技术在井下作业中得到了广泛应用,并在油田井下检泵作业过程中真正实现了清洁生产,取得了显著的经济和社会效益。     

GROUND WATER MONITORING NETWORK DESIGN USING MULTIPLE CRITERIA DECISION MAKING ANT) GEOSTATISTICS1

ABSTRACT: A multi-criteria approach to ground water quality monitoring network design is developed. The methodology combines multi-criteria decision making (MCDM) and modifications of geostatistical variance reduction analysis. Composite programming, a distance based optimization algorithm that employs a hierarchial structure, is used for the MCDM component of the design methodology. MCDM allows the consideration of numerous, often conflicting, design criteria. The methodology is useful for identifying the preferred combination of direct borehole and indirect geo-electric data. It also permits the use of prior information during initial stages of network development. Multi-variate kriging is employed to evaluate network performance using the combination of direct borehole data and indirect geoelectric data. Weighted measures of estimation variance are used as primary measures of performance, with the reduction in estimation variance being computed by the fictitious point method. Case study results demonstrate that the network design methodology can be used in both early and late phases of network development. It also leads to selection of the preferred combination and spatial orientation of direct and indirect data sources while considering cost-effectiveness and performance of alternative designs.