青岛李哥庄地下咸水体恢复方案的数值优化

海水入侵导致地下水变咸对青岛李哥庄地区经济发展产生了不利影响,急需开展治理该地区咸水体的工作。根据研究区已有的降雨、蒸发及水文地质等资料,结合现场调查和监测,利用地下水模拟软件Visual-MODFLOW中的SEAWAT模块建立了该地区的数值模型,并利用实测资料识别和验证了该模型参数,表明所建模型能反映李哥庄地区的实际水文地质条件。为了模拟研究区地下咸水体的恢复治理情况,提出了连续抽水、间歇式抽水和抽注水结合三种咸水恢复方案,并利用该模型对三种恢复方案进行了模拟优化。结果表明,抽取地下咸水只是咸水恢复的一个因素,周边淡水的驱替作用则是咸水恢复另一个的因素,抽注水结合恢复方案为三种方案中最优的方案。     

Modeling approaches and strategies for data-scarce aquifers: example of the Dar es Salaam aquifer in Tanzania

Management of groundwater resources can be improved by using groundwater models to perform risk analyses and to improve development strategies, but a lack of extensive basic data often limits the implementation of sophisticated models. Dar es Salaam in Tanzania is an example of a city where increasing groundwater use in a Pleistocene aquifer is causing groundwater-related problems such as saline intrusion along the coastline, lowering of water-table levels, and contamination of pumping wells. The lack of a water-level monitoring network introduces a problem for basic data collection and model calibration and validation. As a replacement, local water-supply wells were used for measuring groundwater depth, and well-top heights were estimated from a regional digital elevation model to recalculate water depths to hydraulic heads. These were used to draw a regional piezometric map. Hydraulic parameters were estimated from short-time pumping tests in the local wells, but variation in hydraulic conductivity was attributed to uncertainty in well characteristics (information often unavailable) and not to aquifer heterogeneity. A MODFLOW model was calibrated with a homogeneous hydraulic conductivity field and a sensitivity analysis between the conductivity and aquifer recharge showed that average annual recharge will likely be in the range 80–100 mm/year.     

Groundwater resource sustainability in the Wadi Watir delta, Gulf of Aqaba, Sinai, Egypt

The Wadi Watir delta, in the arid Sinai Peninsula, Egypt, contains an alluvial aquifer underlain by impermeable Precambrian basement rock. The scarcity of rainfall during the last decade, combined with high pumping rates, resulted in degradation of water quality in the main supply wells along the mountain front, which has resulted in reduced groundwater pumping. Additionally, seawater intrusion along the coast has increased salinity in some wells. A three-dimensional (3D) groundwater flow model (MODFLOW) was calibrated using groundwater-level changes and pumping rates from 1982 to 2009; the groundwater recharge rate was estimated to be 1.58?×?10⁶ m³/year. A variable-density flow model (SEAWAT) was used to evaluate seawater intrusion for different pumping rates and well-field locations. Water chemistry and stable isotope data were used to calculate seawater mixing with groundwater along the coast. Geochemical modeling (NETPATH) determined the sources and mixing of different groundwaters from the mountainous recharge areas and within the delta aquifers; results showed that the groundwater salinity is controlled by dissolution of minerals and salts in the aquifers along flow paths and mixing of chemically different waters, including upwelling of saline groundwater and seawater intrusion. Future groundwater pumping must be closely monitored to limit these effects.     

0-1整数规划在水源地开采井最优布局中的应用研究

新建水源地如何最优布井,对水源地开采运行具有重大影响。本文介绍了0-1整数规划及其求解方法,并结合一个实际问题进行说明。作者应用0-1整数规划方法建立了北京密怀顺平原区拟建水源地最优布井的地下水管理模型,采用响应矩阵法将地下水模拟模型和管理模型相耦合,运用分支-定界法求出最优井位,并将优化井位与原拟定井位的计算结果进行了对比。结果表明,0-1整数规划能很好地解决最优布井问题,对水源地规划、基坑降水、地下水人工回灌等方面都具有一定的指导意义。     

应用解析法确定傍河水源地地下水开采方案

合理设计开采方案对傍河水源地地下水资源开发具有重要意义。以德惠典型傍河地下水水源地为研究区,在查明研究区水文地质条件的基础上,以管井为例,运用解析法确定了最佳开采方案和水源地地下水允许开采量,并对最佳开采方案的合理性和可靠性进行了分析。结果表明,当开采井沿河流呈直线形单排布井方式,井数为12眼,井河距离为200 m,井间距离为150 m,单井开采量为2 960 m3/d,总开采量为35 520 m3/d时,取水保证程度较高,面对河流突发性污染具有一定的抵御能力,不会对区域水环境和其他用户造成较大影响。     

Groundwater resources assessment for irrigable agricultural lands in the Wadi Araba area, southern Jordan

The main target of this research paper was to the hydrogeological assessment of the groundwater resources to irrigate 600 ha of irrigable agricultural lands, distributed along the Dead Sea–Aqaba Highway in Umm, Methla, Wadi Musa, Qa’ Saideen and Rahma, southern Jordan. Therefore, a comprehensive groundwater study was commenced by drilling eight new wells which can be used to supply irrigable areas with the existing groundwater that would be enriched by the yield of three proposed recharge dams on Wadi Musa, Wadi Abu-Burqa, and Wadi Rahma. The evaluation of the pumping test data of the drilled was carried out using the standard methods of pumping test interpretation. This was based on the available water table measurements at well locations and knowledge of water flow in the general. The sustainable yield of each well was calculated based on the pumping test parameters. The obtained results indicate that pumping out of Beer Mathkor wells should not exceed 1,100 m³/day in the case of continuous pumping and 8,700 m³/day in the case of intermittent pumping. Since the water table did not significantly change with small changes in pumping (it took eightfolds of magnitude increase in pumping from approximately 1,100 to 8,700 m³/day to show a significant drop in the water table equivalent to about 5.5 MCM per year from the aquifer.     

Use of stream response functions to determine impacts of replacing surface-water use with groundwater withdrawals

A regional-scale numerical groundwater model is used to study the impacts of replacing surface-water use with groundwater wells to improve low-flow stream conditions for endangered species within the Bertrand and Fishtrap watersheds, southern British Columbia, Canada and Washington, USA. Stream response functions ranging from 0 to 1.0 were calculated for individual wells placed within a steady-state groundwater flow model at varying distances from the streams to determine the impact that these replacement wells, operating under sustained pumping rates, would have on summer instream flows. Lower response ratios indicate groundwater pumping will have less of an impact on streamflow than taking an equivalent amount of water directly from a surface-water source. Results show that replacing surface-water use with groundwater withdrawals may be a viable alternative for increasing summer streamflows. Assuming combined response factors should be ≤0.5 for irrigators to undergo the expense of installing new wells, ~57% of the land area within 0.8 km of Bertrand Creek would be suitable for replacement wells. Similarly, 70% of the land area within 0.8 km of Fishtrap Creek was found to be appropriate. A visual analysis tool was developed using STELLA to allow stakeholders to quickly evaluate the impact associated with moving their water right.     

Spatial variability estimation and risk assessment of the aquifer level at sparsely gauged basins using geostatistical methodologies

The spatial variability evaluation of the water table level of an aquifer provides useful information in water resources management plans. Three different approaches are applied to estimate the spatial variability of the water table in the study basin. All of them are based on the Kriging methodology. The first is the classical Ordinary Kriging approach, while the second involves information from a secondary variable (surface elevation) and the application of Residual Kriging. The third calculates the probability to lie below a certain groundwater level limit that could cause significant problems in groundwater resources availability. The latter is achieved by means of Indicator Kriging. A recently developed non-linear normalization method is used to transform both data and residuals closer to normal distribution for improved prediction results. In addition, the recently developed Spartan variogram model is applied to determine the spatial dependence of the measurements. The latter proves to be the optimal model, compared to a series of models tested, which provides in combination with the Kriging methodologies the most accurate cross validation estimations. The variogram form is explained with respect to the radius of influence of the pumping wells representing the spatial impact of the pumping activity. Groundwater level and probability maps are developed providing the ability to assess the spatial variability of the groundwater level in the basin and the risk that certain locations have in terms of a safe groundwater level limit that has been set for the sustainability of the groundwater resources of the basin.     

Change of groundwater chemistry from 1896 to present in the Mid-Levels area, Hong Kong

In this study, groundwater quality information collected in 1896 (well waters), 1980/1981 (piezometric and seepage samples) and 2002/2003 (seepage samples) in the regions centered by the Mid-Levels area, Hong Kong Island, was compared to illustrate how groundwater quality has changed over a century and the processes controlling it. As shown by saline ammonia and nitrate levels in the late nineteenth century, groundwater was severely polluted by widespread and obvious leakage from poorly designed wastewater collection systems, although groundwater was still a drinking water source for local residents. The extremely high residual chlorines in groundwater demonstrated that large doses of disinfection agents were added to wells at that time. In view of the decline in saline ammonia and nitrate levels, groundwater became less organically polluted in the 1980s probably due to significant improvement of the design of underground sewers. However, more leakage from sources such as salty flushing water and fresh water pipes emerged in the past few decades which added complexity to groundwater chemical systems. Some chemicals were used to identify possible locations of leakages. The temporal variations of the distribution of these chemicals over the area may shed light on the rate of leakage. Leakage from service pipes seems to have improved from the early 1980s to 2002/2003. However, the area is still suffering from widespread and small-scale leakage from service pipes. More efforts should be paid to control small leakages in the future. The findings will be instructive to various government organizations such as the Water Supplies Department and Drainage Services Department to identify possible locations of unobvious leakages in the area.     

In-well time-of-travel approach to evaluate optimal purge duration during low-flow sampling of monitoring wells

A common assumption with groundwater sampling is that low (<0.5 L/min) pumping rates during well purging and sampling captures primarily lateral flow from the formation through the well-screened interval at a depth coincident with the pump intake. However, if the intake is adjacent to a low hydraulic conductivity part of the screened formation, this scenario will induce vertical groundwater flow to the pump intake from parts of the screened interval with high hydraulic conductivity. Because less formation water will initially be captured during pumping, a substantial volume of water already in the well (preexisting screen water or screen storage) will be captured during this initial time until inflow from the high hydraulic conductivity part of the screened formation can travel vertically in the well to the pump intake. Therefore, the length of the time needed for adequate purging prior to sample collection (called optimal purge duration) is controlled by the in-well, vertical travel times. A preliminary, simple analytical model was used to provide information on the relation between purge duration and capture of formation water for different gross levels of heterogeneity (contrast between low and high hydraulic conductivity layers). The model was then used to compare these time–volume relations to purge data (pumping rates and drawdown) collected at several representative monitoring wells from multiple sites. Results showed that computation of time-dependent capture of formation water (as opposed to capture of preexisting screen water), which were based on vertical travel times in the well, compares favorably with the time required to achieve field parameter stabilization. If field parameter stabilization is an indicator of arrival time of formation water, which has been postulated, then in-well, vertical flow may be an important factor at wells where low-flow sampling is the sample method of choice.     

Optimal Groundwater Development in Coastal Aquifers Near Beihai, China

INTRODUCTIONThe city of Beihai,located on the south coast ofGuangxi,China,relies heavily on groundwater for its potablewater supply and agricultural irrigation.With rapid increasein population (for instance,from 134 0 0 0 in 1987to 47930 0in1995 ) and in developm ent program s,the demand for freshwater has been growing. Approxim ately 170 0 0 0 m3/ d ofgroundwater has been pumped from the productive coastalaquifers in recent years.Contamination of the fresh water inthe coastal aquifers b…     

Case history of one of the few successful Superfund remediation sites: A site at Salinas, California, USA

In 1985, a former tire manufacturing plant surrounded by agricultural fields in the Salinas Valley was designated a Superfund Site by the US Environmental Protection Agency. The plant had been operating for seventeen years, from 1963 until 1980. When dismantling of the plant was started, it was determined that toxic hydrocarbon solvents and oils from the plant had contaminated soil and groundwater in alluvial deposits alongside the plant. It was determined later that the groundwater contamination also lay beneath the agricultural fields in a narrow groundwater plume that extends about 4.3 km downgradient from the plant. Because of the complex architecture of the aquifer system, the gradient, and extensive pumping of agricultural wells, the contaminants migrated northwestward and downward to deeper levels away from the plant.

The agricultural fields are underlain by an unconfined shallow aquifer and by a system of confined aquifers that extend to more than 180 m below surface. Aquifers are discontinuous beds of channel sand and gravel; confining beds are overbank clay and silt, and estuarine clay.

Geophysical data, logs of existing agricultural and other wells, and careful consideration of the stratigraphic architecture of the depositional environment provided the basis for a conceptual hydrogeologic model and for locating characterization wells for detailed visual and geophysical logging and hydrologic testing. Successive refinements of the characterization by sequential installation of wells indicated optimal locations for installation of extraction and monitoring wells. Validity of the concept of the hydrogeologic regime was verified by close match of predictions made by modeling with the later results of pumping from the extraction wells in a pump-and-treat system.

Successful remediation was accomplished by analyzing data from 110 agricultural wells, the few domestic water wells, nearly 200 sequentially installed stratigraphic-characterization and monitoring wells, 25 extraction wells, and by close cooperation among federal, state and local agencies, and the ranchers and growers.

Total contaminants recovered from activated-charcoal strippers of the treatment system totalled < 230 kg. Large quantities were harmlessly volatilized and dispersed into the atmosphere by air strippers and by agricultural sprinkling systems spraying water onto the fields. Crop testing showed no contamination of food crops. The activity has taken seven years and has cost more than US$22 million.     

An integrated approach to investigate saline water intrusion and to identify the salinity sources in the Central Godavari delta, Andhra Pradesh, India

The Central Godavari delta is located along the Bay of Bengal Coast, Andhra Pradesh, India, and is drained by Pikaleru, Kunavaram and Vasalatippa drains. There is no groundwater pumping for agriculture as wells as for domestic purpose due to the brackish nature of the groundwater at shallow depths. The groundwater table depths vary from 0.8 to 3.4 m and in the Ravva Onshore wells, 4.5 to 13.3 m. Electrical Resistivity Tomography (ERT) surveys were carried out at several locations in the delta to delineate the aquifer geometry and to identify saline water aquifer zones. Groundwater samples collected and analyzed for major ions for assessing the saline water intrusion and to identify the salinity origin in the delta region. The results derived from ERT indicated low resistivity values in the area, which can be attributed to the existence of thick marine clays from ground surface to 12–15 m below ground level near the coast and high resistivity values are due to the presence of coarse sand with freshwater away from the coast. The resistivity values similar to saline water <0.01 Ω m is attributed to the mixing of the saline water along surface water drains. In the Ravva Onshore Terminal low resistivity values indicated up coning of saline water and mixing of saline water from Pikaleru drain. The SO ₄ ^?2 /Cl^?and Na⁺²/Cl^?ratios did not indicate saline water intrusion and the salinity is due to marine palaeosalinity, dilution of marine clays and dissolution of evaporites.     

Application of Numerical Modeling for Groundwater Flow System Analysis in the Akaki Catchment,Central Ethiopia

A three dimensional steady-state finite difference groundwater flow model is used to quantify the groundwater fluxes and analyze the subsurface hydrodynamics in the Akaki catchment by giving particular emphasis to the well field that supplies water to the city of Addis Ababa. The area is characterized by Tertiary volcanics covered with thick residual and alluvial soils. The model is calibrated using head observations from 131 wells. The simulation is made in a two layer unconfined aquifer with spatially variable recharge and hydraulic conductivities under well-defined boundary conditions. The calibrated model is used to forecast groundwater flow pattern, the interaction of groundwater and surface water, and the effect of pumping on the well field under different scenarios. The result indicates that the groundwater flows regionally to the south converging to the major well field. Reservoirs and rivers play an important role in recharging the aquifer. Simulations made under different pumping rate indicate that an increase in pumping rate results in substantial regional groundwater level decline, which will lead to the drying of springs and shallow hand dug wells. Also, it has implications of reversal of flow from contaminated rivers into productive aquifers close to main river courses. The scenario analysis shows that the groundwater potential is not enough to sustain the ever-growing water demand of the city of Addis Ababa. The sensitivity and scenario analysis provided important information on the data gaps and the specific sites to be selected for monitoring, and may be of great help for transient model development. This study has laid the foundation for developing detailed predictive groundwater model, which can be readily used for groundwater management practices.     

辽河油田地下水资源优化利用方案研究

辽河油田以开采上第三系地下水为主，已出现降落漏斗、水质恶化及咸水入侵等问题。文章应用Visual MODFLOW模型，建立了地下水流数值模拟模型，对不同开采方案下的地下水动态变化进行了预测和分析，确定了最佳开采量和最优开采布局。建议调减开采量和调整地下水开采布局以实现辽河油田地区水资源可持续利用。     

Determining the optimal pumping duration of transient pumping tests for estimating hydraulic properties of leaky aquifers using global curve-fitting method: a simulation approach

Global curve-fitting method (GCFM) is regarded as an effective approach in hydrogeological parameter estimation, as it integrates and uses pumping data and water recovery data of a transient pumping test for parameter estimation. The impacts of pumping duration on hydrogeological parameter estimation by GCFM were investigated in the present study using 2 in situ pumping tests and 24 simulated transient pumping tests. Empirical formulas for determining the optimal pumping duration were derived. The study results suggest that pumping duration will have impacts on the accuracy of hydrogeological parameter estimation. When pumping duration is longer than a certain period, relative errors of hydrogeological parameter estimation keep relatively stable within an acceptable limit. Therefore, it is unnecessary to continue the pumping for a very long time after the groundwater level has become stable. When the change rate of drawdown over time (γ) in an observation well located within a distance of 100 m to the pumping well reaches 0.134, the pumping can be stopped. If there are more than one observation wells in a pumping test, the smallest γ value should be selected to determine the optimal pumping duration. This research is meaningful in the instruction of pumping tests, and will reduce test costs greatly.     

Hydrothermal contamination of public supply wells in Napa and Sonoma Valleys,California

Groundwater chemistry and isotope data from 44 public supply wells in the Napa and Sonoma Valleys, California were determined to investigate mixing of relatively shallow groundwater with deeper hydrothermal fluids. Multivariate analyses including Cluster Analyses, Multidimensional Scaling (MDS), Principal Components Analyses (PCA), Analysis of Similarities (ANOSIM), and Similarity Percentage Analyses (SIMPER) were used to elucidate constituent distribution patterns, determine which constituents are significantly associated with these hydrothermal systems, and investigate hydrothermal contamination of local groundwater used for drinking water. Multivariate statistical analyses were essential to this study because traditional methods, such as mixing tests involving single species (e.g. Cl or SiO₂) were incapable of quantifying component proportions due to mixing of multiple water types. Based on these analyses, water samples collected from the wells were broadly classified as fresh groundwater, saline waters, hydrothermal fluids, or mixed hydrothermal fluids/meteoric water wells. The Multivariate Mixing and Mass-balance (M3) model was applied in order to determine the proportion of hydrothermal fluids, saline water, and fresh groundwater in each sample. Major ions, isotopes, and physical parameters of the waters were used to characterize the hydrothermal fluids as Na–Cl type, with significant enrichment in the trace elements As, B, F and Li. Five of the wells from this study were classified as hydrothermal, 28 as fresh groundwater, two as saline water, and nine as mixed hydrothermal fluids/meteoric water wells. The M3 mixing-model results indicated that the nine mixed wells contained between 14% and 30% hydrothermal fluids. Further, the chemical analyses show that several of these mixed-water wells have concentrations of As, F and B that exceed drinking-water standards or notification levels due to contamination by hydrothermal fluids.     

Modeling of seawater intrusion in a coastal aquifer of Karaburun Peninsula,western Turkey

Seawater intrusion is a major problem to freshwater resources especially in coastal areas where fresh groundwater is surrounded and could be easily influenced by seawater. This study presents the development of a conceptual and numerical model for the coastal aquifer of Karareis region (Karaburun Peninsula) in the western part of Turkey. The study also presents the interpretation and the analysis of the time series data of groundwater levels recorded by data loggers. The SEAWAT model is used in this study to solve the density-dependent flow field and seawater intrusion in the coastal aquifer that is under excessive pumping particularly during summer months. The model was calibrated using the average values of a 1-year dataset and further verified by the average values of another year. Five potential scenarios were analyzed to understand the effects of pumping and climate change on groundwater levels and the extent of seawater intrusion in the next 10 years. The result of the analysis demonstrated high levels of electrical conductivity and chloride along the coastal part of the study area. As a result of the numerical model, seawater intrusion is simulated to move about 420 m toward the land in the next 10 years under “increased pumping” scenario, while a slight change in water level and TDS concentrations was observed in “climate change” scenario. Results also revealed that a reduction in the pumping rate from Karareis wells will be necessary to protect fresh groundwater from contamination by seawater.     

Managing groundwater rise: Experimental results and modelling of water pumping from a quarry lake in Milan urban area (Italy)

An innovative approach to solve the problem of lowering water table was carried out in a quarry lake south of the city of Milan (northern Italy): the project, based upon pumping out water at a rate of 1,000 L/s can be considered a strategic medium to long-term solution to hinder the rise of groundwater level interfering with underground structures (foundation, construction, subway) in urban areas. The basic idea is to pump a high groundwater rate as close as possible to the stagnation point of the piezometric depression located in the city. After a pilot-test was carried out in November 1998, experimental activities started in July 2001 and lasted one year; water withdrawal was discharged into artificial channels used in agricultural practice. Maximum drawdowns measured in the quarry lake by the monitoring network resulted in more than 5 m, and a significant drawdown was registered up to 1.5 km of distance from the quarry in the important historical site of Chiaravalle Abbey, threatened by groundwater rise. The results of this pumping activity confirm the importance of the project, its lower cost compared with traditional solutions (such as drainage by wells) and remarkable effects on the improvement of surface water quality. A groundwater model was implemented to evaluate further scenarios of discharge rate and pumping location, too.     

Migration of arsenic in multi-aquifer system of southern Bengal Basin: analysis via numerical modeling

A heterogeneous anisotropic steady-state groundwater flow model for the multi-aquifer system of a part of southern Bengal Basin shows that human intervention has changed the natural groundwater flow system. At present, the shallow groundwater flow is restricted within the aquifer, with very short travel time of tens of years and vertical path length. The deep aquifer is fed by surface water or rainwater from distant locations with travel time of thousands of years and has no hydraulic connection with the arsenic-rich shallow aquifer. Numerical simulations indicate that the future pumping of deep groundwater is not likely to drive in arsenic from the shallow aquifer. Therefore, new wells may be installed in the deep aquifer. High pumping of shallow unpolluted aquifer consisting of brown sand will drive in groundwater containing organic matter from the post-Last Glacial Maximum aquifer-aquitard system. The organic matter drives reduction of manganese oxides at strip interfaces between palaeo-channel and palaeo-interfluve. After the completion of manganese reduction, FeOOH reduction may take place in the marginal palaeo-interfluvial aquifer and release sorbed arsenic. Arsenic then moves into the interior of palaeo-interfluvial aquifer polluting its fresh groundwater. Arsenic migration rates ranges between 0.21 and 6.3 and 1.3 × 10^?2 and 0.4 m/year in horizontal and vertical directions, respectively. Therefore, palaeo-interfluvial aquifer will remain arsenic-free for hundreds to thousands of years to supply safe drinking water.