划分地下水源地保护区的数值模拟方法

本文提出了一种划分地下水水源地保护区的方法,首先用地下水模拟系统中的MODFLOW建立研究区的数值模拟模型,获得地下水流场;然后再将水源地划分为3级保护区,用MODPATH确定各保护区的范围;最后根据国家饮用水标准,采用MT3D确定主要污染质在保护区中的可释放阈值,并给出了应用实例。     

Groundwater flow modelling of Hindon-Yamuna interfluve region, Western Uttar Pradesh

The study area Hindon -Yamuna interfluve region is underlain by a thick pile of unconsolidated Quaternary alluvial deposits and host multiple aquifer system. Excessive pumping in the last few decades, mainly for irrigation, has resulted in a significant depletion of the aquifer. Therefore, proper groundwater management of Hindon-Yamuna interfluve region is necessary. For effective groundwater management of a basin it is essential that careful zone budget study should be carried out. Keeping this in view, groundwater flow modelling was attempted to simulate the behavior of flow system and evaluate zone budget. Visual MODFLOW, pro 4.1 is used in this study to simulate groundwater flow. The model simulates groundwater flow over an area of about 1345 km² with a uniform grid size of 1000 m by 1000 m and contains three layers, 58 rows and 37 columns. The horizontal flows, seepage losses from unlined canals, recharge from rainfall and irrigation return flows were applied using different boundary packages available in Visual MODFLOW, pro 4.1. The river — aquifer interaction was simulated using the river boundary package. Simulated pumping rates of 500 m³/day, 1000 m³/day and 1500 m³/day were used in the pumping well package.The zone budget for the steady state condition of study area indicated that the total annual direct recharge is 416.10 MCM and the total annual groundwater draft through pumping is of the order of 416.63 MCM. Two scenarios were considered to predict aquifer system response under different conditions. Sensitivity analysis on model parameters was conducted to quantitatively evaluate the impact of varying model inputs. Based on the results obtained from the sensitivity analysis, it was found that the model is more sensitive to hydraulic conductivity and recharge parameter. Present study deals with importance of groundwater modelling for planning, design, implementation and management of groundwater resources.     

宁夏甜水河煤田地下水模拟及涌水量预测

地下水模拟技术是进行煤田地下水环境影响评价和矿区疏干排水设计的主要工具。该文以宁夏甜水河煤田为研究对象,根据矿区地质条件和水文地质条件的分析,确定模型边界条件并建立了水文地质概念模型,应用基于有限差分法的MODFLOW程序,建立地下水流三维模型,进而对煤田各开采阶段的地下水流场和涌水量进行预测,从而为煤田地下水有效防治和安全生产提供依据。     

Modelling nitrate pollution pressure using a multivariate statistical approach: the case of Kinshasa groundwater body,Democratic Republic of Congo

A multivariate statistical modelling approach was applied to explain the anthropogenic pressure of nitrate pollution on the Kinshasa groundwater body (Democratic Republic of Congo). Multiple regression and regression tree models were compared and used to identify major environmental factors that control the groundwater nitrate concentration in this region. The analyses were made in terms of physical attributes related to the topography, land use, geology and hydrogeology in the capture zone of different groundwater sampling stations. For the nitrate data, groundwater datasets from two different surveys were used. The statistical models identified the topography, the residential area, the service land (cemetery), and the surface-water land-use classes as major factors explaining nitrate occurrence in the groundwater. Also, groundwater nitrate pollution depends not on one single factor but on the combined influence of factors representing nitrogen loading sources and aquifer susceptibility characteristics. The groundwater nitrate pressure was better predicted with the regression tree model than with the multiple regression model. Furthermore, the results elucidated the sensitivity of the model performance towards the method of delineation of the capture zones. For pollution modelling at the monitoring points, therefore, it is better to identify capture-zone shapes based on a conceptual hydrogeological model rather than to adopt arbitrary circular capture zones.     

Groundwater flow modelling of Yamuna-Krishni interstream,a part of central Ganga Plain Uttar Pradesh

Groundwater is a major source of water for agricultural and domestic requirements in western Uttar Pradesh. Due to increasing agricultural requirements the abstraction of groundwater has increased manifold in the last two-to-three decades. The quaternary alluvium hosts the aquifer in the region. The study area forms a part of Yamuna-Krishni interfluve. Although the area hosts potential aquifers these have been adversely affected by poor management. For effective groundwater management of a basin it is essential that a careful water balance study should be carried out. Keeping this in mind groundwater flow modelling was attempted to simulate the behaviour of the flow system and evaluate the water balance. The groundwater flow modelling was carried out. The horizontal flows, seepage losses from unlined canals, recharge from rainfall and irrigation return flows were applied using different boundary packages available in Visual MODFLOW, Pro 4.1. The river-aquifer interaction was simulated using the river boundary package. Hydraulic conductivity values were applied to specific zones and these ranged from 9.8 to 26.6m/day. Recharge due to rainfall and irrigation returns were assigned to respective zones. Pumping rates of 500m³/day, 1000m³/day, 1500m³/day, 2000m³/day and 2500m³/day were applied to appropriate areas of the model to simulate areas of stress. The zone budget shows a water balance deficit for the period June 2006 to June 2007. The total recharge to the study area is 160.21 million m³ (Mcum). The groundwater draft through pumping is of the order of 233.56 Mcum, thus leaving a deficit balance of −73.35 Mcum. The sensitivity of the model to input parameters was tested by varying the parameters of interest over a range of values, monitoring the response of the model and determining the root mean square error of the simulated groundwater heads to the measured heads. These analyses showed that the model is most sensitive to hydraulic conductivity and recharge parameters. Three scenarios were considered to predict aquifer responses under varied conditions of groundwater bstraction.     

基于GMS的三维TOUGH2模型及模拟

GMS和TOUGH2均是应用广泛的地下流动系统数值模拟软件,其中,GMS界面友好,功能强大;而TOUGH2虽具有强大的数值计算能力,但缺乏友好的可视化前后处理界面。本文借助GMS强大的前处理能力,基于概念模型建立了三维复杂模型,把GMS/MODFLOW三维数值模型(包括网格数据、岩性数据、初始和边界数据等)转化为TOUGH2数值模型进行数值计算。通过两个计算实例(含水层水平没有起伏和有起伏)对比分析了GMS/MODFLOW和TOUGH2计算结果的差异。结果显示,该方法可以快速地建立刻画复杂地质条件的TOUGH2模型,计算结果与GMS/MODFLOW差异很小,说明两个软件均有很高的可信度;同时,该方法发挥了两个软件各自的优势,为进行更为复杂的多相流动数值模拟提供了可行性。     

Climate and irrigation scenario analyses using three-dimensional numerical modelling: a case study of the Nasia sub-basin in the White Volta Basin,Ghana

A groundwater resource characterisation and assessment model was developed for Nasia river sub-basin in the White Volta Basin, Ghana. The model is useful to policymakers for planning and sustainable management of groundwater resources in the basin for domestic and irrigation purposes. A conceptual model was constructed that characterized boundary conditions and hydrostratigraphy, and estimated recharge rates and hydraulic and storage parameters. From current understanding of the hydrogeological dynamics, three hydrostratigraphic layers were delineated. The conceptual model was converted to a three-dimensional steady-state groundwater flow model using MODFLOW. Recharge rates estimated from the base model indicate a minimum of 1.1% and maximum of 6.2% of the total rainfall. The hydraulic conductivity ranged between 0.20 and 15 m/day. Four possible scenarios were simulated: (1) increased population, (2) climate variations (reduced recharge), (3) increased abstraction for irrigation, and (4) worst-case scenario which is a combination of the first three scenarios. Results from scenarios 1 and 2 indicated that, under such conditions, the groundwater resources could be sustained and no significant effect on any of the water budget indicators was observed. For scenario 3, there was significant drop in hydraulic head in the central portions of the study area. The scenario 4 simulation indicated that there was significant reduction in groundwater levels and groundwater discharge into streams under these stressors. Such reduction can affect stream levels in the basin and, subsequently, the ecosystem. These findings are valid within the limits of uncertainty in the hydrogeological data that were used in this study.

     

基于MODFLOW的地下水流模型前处理优化

MODFLOW中现有的降水补给数据和井流数据是基于每个剖分网格输入的,导致基于MODFLOW建立大区域水流数值模型前处理的降水补给文件和井流文件存储开销过大、读取效率低。为此,通过改进现有的降水补给子程序包RCH及开发新的子程序包RAW,给出了一种基于面状补排项的数据输入新方法。改进后的MODFLOW程序通过读取RCH文件或RAW文件中每个分区的补排量数据,以及每个网格对应的分区编号,在程序内部实现了补排量向每个模型网格的分配。RAW子程序包实现了多层面状补排量的表达,可用于面状的地下水开采、农业回灌等源汇项的处理。相对于原始的源汇项数据存储方式,基于新方法建立的华北平原地下水流模型,RCH及RAW文件大小分别减小为原来的1/145和1/255,整个模型数据的读取时间的加速比为7.46。     

Geographical information system approach for environmental management in coastal area (Hardelot-Plage,France)

The use of geographic information system (GIS) minimizes the effort and improves the efficiency of numerical models. The GIS provides a platform for high capacity collection, management, manipulation, analysis, modeling and display of spatial data. The conceptual model is created using GIS objects including points, arcs and polygons so that it can accurately represent real world condition. According to the research problem, the geographical model is based on Hypergraph Based Data Structure method, and a conceptual data model has been created from which a physical data model was elaborated in ArcGIS9.3 platform. The groundwater modeling system (GMS) provides a powerful tool for hydrodynamics modeling and it is able to solve complex problems such as the groundwater flow and seawater intrusion. The sand-dune system of Hardelot-Plage (North of France) suffers from a lack of well-developed foredune. This problem is linked to the almost constant saturation of beach sand which is the potential source of dune nourishment. In the south of Hardelot, the coastline is slowly, but constantly retreating. To remedy this situation, a coupling between a GIS and GMS was adopted, in order to find the possible scenarios which could lower the piezometric surface in the concerned area and allow dune nourishment again. The GMS used supports the Modflow-2000 code. A direct approach to designing Modflow finite difference model is tedious and less intuitive, specifically for complex boundary and initial conditions. Therefore, a Modflow model can be developed either using a grid or conceptual model approach. The preparation of input data modeling is tedious and takes a long time. The model created in GMS was calibrated against the historical and observed water level data for 1995–2006. Then a hydrodispersive model (MT3d code in GMS) was launched for evaluating sea-water intrusion. The model was run to generate groundwater and salt concentration scenario during pumping tests.     

Simulation of a groundwater fall caused by geological discontinuities

In some areas close to a reverse fault with significant vertical displacement, or the boundary of a mining area in multi-layer rock and the boundary of perched water, the phreatic water-level changes abruptly from one side of the fault or boundary to the other. There is no direct hydraulic connection between the two sides and a groundwater fall occurs across the fault or boundary. This groundwater fall cannot be easily handled by the modules available in MODFLOW. Semi-analytical solutions were derived for water discharge from the upstream to the downstream side of the fault or boundary, and verified by air-water two-phase numerical simulation. The numerical results show that the derived semi-analytical solutions can be used to accurately estimate the groundwater flux at the fault even though the unsaturated zone is not considered. A module package, GWF, was developed to simulate the groundwater fall, which can be embedded directly into MODFLOW. A theoretical example is presented to show how the package GWF is used to simulate perched water. This package was also applied successfully to build a regional groundwater model of the Urumqi River Basin, Xinjiang, China, and the simulation results showed good agreement with the local hydrogeologic conditions.     

区域含水层系统结构的体视化表达及地下水位预测——以江苏省宿迁市为例

本文结合宿市区域水文地质条件的实际状况,利用GMS对宿迁市区域含水层系统结构进行了体视化表达,使建立的区域水文地质概念模型更加直观,在概念模型的基础上建立了相应的数学模型,并利用Visual MODFLOW对宿迁市不同开采条件下的地下水位变化进行了预测。结果表明：2009年末将出现两个较大的水位降落漏斗区。根据模型预测结果,提出了该区地下水合理开采的建议。     

Visual MODFLOW在大庆龙西地区地下水数值模拟中的应用

Visual MODFLOW是综合MODFLOW、MODPATH、MT3D、RT3D和WinPEST等地下水模型而开发的可视化地下水模拟软件，具有水的质点向前、向后示踪流线模拟、任意区域的水均衡、简化数值模拟数据前处理和后处理，将模拟的复杂性降到最小等特点。龙西地区是大庆油田新增水源地，使用Visual MODFLOW在该地区进行了地下水资源的模拟与预报，建立准三维模型，发现PCG2求解方法精确度最高。另外，把随机模型和MODFLOW模型相结合，解决了边界水位在计算中不断变化的问题。按照设计开采量对地下水水位进行预报，根据模型预报结果，水：原地设计开采量是可行的。Visual MODFLOW成功地运用于大庆龙西地区地下水资源的评价，操作简便，预报准确、可靠，说明具有很好的应用和推广前景。     

山东淄博市地下水资源评价及其合理开发利用研究

本文采用有限单元法进行地下水资源评价以解决评价区供水面积大、精度要求高的矛盾。首先根据实际的开采量和长期观测孔的水位资料,校正数学模型,反求有关参数;然后给定开采方案进行未来水头的预测;最后研究地下水资源的合理开发利用问题。利用有限单元法进行水源地或单一含水层的地下水资源评价已被证明是先进的、有效的。本文是对水文地质条件非常复杂的大面积的区域,用有限单元法进行地下水资源评价的一种尝试,可供类似地区参考。     

鄂尔多斯白垩系地下水盆地水文地质概念模型

水文地质概念模型是地下水模拟与水资源评价的基础。为了构建地下水数值模拟模型,以鄂尔多斯白垩系地下水盆地为研究对象,以丰富翔实的勘查资料为依据,从白垩系地下水系统边界条件、水文地质结构、水流系统、水文地质参数、源汇项等方面进行概化,联合采用ArcGIS、RS、GMS同位素、基流分割、地质统计学、随机模拟等技术方法,建立了一组水文地质概念模型及其相应的空间离散表达,并讨论了地下水模拟面临的问题与发展趋势。     

第二松花江干流傍河开采潜力评价

在第二松花江傍河开采适宜性评价的基础上,根据研究区的地质和水文地质条件,建立了水文地质概念模型及相应的数学模型,建立了基于Visual MODFlow的地下水流数值模型;以地下水水位降深不大于含水砂层厚度的1/3为准则,通过使预测流场与允许降深条件下的流场尽量接近,给出了第二松花江干流傍河开采条件下的地下水允许开采量。结果表明:现状年条件下区域地下水资源开采强度较小,而傍河开采条件下各适宜区地下水开采潜力巨大,在多年平均降水保证率条件下,地下水允许开采量可达140×104 m3/d;适度规模傍河开采可引起地表水与地下水转化关系的变化,但不会对河流生态功能产生不利影响。     

Numerical simulation of groundwater flow and contamination transport in Nahavand plain aquifer,west of Iran

Numerical simulation of groundwater flow used for the estimation of hydraulic and hydrologic parameters which is an important tool for the management of aquifers. This study presents the results of a mathematical model developed for the simulation of groundwater flow in Nahavand plain aquifer in the southwest Hamadan province. For this purpose Groundwater Modeling Software (GMS) was used which supports the MODFLOW-2000 code. After gathering required data such as the hydrological, hydrogeological and topography maps, a 3D hydrogeological model of plain was constructed with borehole and surface elevation data. Then MODFLOW was used for simulation of flow. After initial simulation of the flow, the model was calibrated in steady state with trial-and-error and parameter estimation methods the observed head of groundwater table monitoring data of 1997. Results of calibration show that error between observed head and computed head is in allowable range. Also results of computed head with model show that groundwater flow is in the direction of the dominate slope (southeast to northwest). Finally MODPATH code which simulates advective transport of particles was used for estimation of flow path and source of contaminants.     

Groundwater resources management through flow modeling in lower part of Bhagirathi — Jalangi interfluve,Nadia, West Bengal

Groundwater utilization in Nadia district, West Bengal has been subjected to rapid exploitation in the wake of increasing urbanization and production of agricultural commodities. It is, therefore, necessary to evaluate the existing trend and availability of groundwater in time and space and its movement for proper planning in future. In the present study, an attempt has been made develop a groundwater management model using Visual MODFLOW software. The groundwater flow model for the study area was formulated by using input hydrogeological data and appropriate boundary conditions. The groundwater flow pattern of the study area indicates the occurrence of base flow which feed both the Rivers Bhagirathi and Jalangi throughout the year. The computed hydraulic heads were calibrated by comparing with observed groundwater level data for years 2004 to 2006 and were verified with the data of 2007. The outcome of modeling shows that this model can be used for prediction purpose in the future by updating input boundary conditions and hydrologic stresses during the preceding year. The model optimized unit draft for deep tube well as 556.5 m³/day and the same for shallow tube well as 41 m³/day keeping the existing tube well structures in running condition and maintaining the present and recent past trends of groundwater level. The model can be further improved if more spatial and temporal input parameters are available and can be incorporated into the model for more realistic characterization of groundwater flow.     

Simulation of horizontal well performance using Visual MODFLOW

A proposed horizontal well or radial collector well installation in shallow aquifers to enhance water withdrawal rates in Pintu Geng well field in Kelantan, Malaysia was simulated using the Drainage Package of MODFLOW groundwater model. The modelling exercise aimed at identifying an optimum pumping rate that would safely achieve the desired drawdown of less than 2 m in an area of 300 m radius surrounding the Pintu Geng horizontal collector well. The model also would serve as a basis for the design of the horizontal well components. High degree of grid refinement for the well location is needed to simulate the real field installation. However, for the purpose of designing water withdrawal systems, it is important to obtain the correct production rate of these wells for a given drawdown. A transient groundwater flow model was calibrated and validated with few assumptions of the horizontal well hydraulic properties. The model demonstrates that under natural flow condition at ?3 m depth, the six collectors (drains) tap a volume of 19,200–43,700 m³/day. A steady-state model was also developed to predict the capture zone delineation. Attention is also given to the impact of the well installation to the surrounding 300 m radius by inspecting the degree of the drawdown.     

Numerical groundwater-flow modeling to evaluate potential effects of pumping and recharge: implications for sustainable groundwater management in the Mahanadi delta region,India

Process-based groundwater models are useful to understand complex aquifer systems and make predictions about their response to hydrological changes. A conceptual model for evaluating responses to environmental changes is presented, considering the hydrogeologic framework, flow processes, aquifer hydraulic properties, boundary conditions, and sources and sinks of the groundwater system. Based on this conceptual model, a quasi-three-dimensional transient groundwater flow model was designed using MODFLOW to simulate the groundwater system of Mahanadi River delta, eastern India. The model was constructed in the context of an upper unconfined aquifer and lower confined aquifer, separated by an aquitard. Hydraulic heads of 13 shallow wells and 11 deep wells were used to calibrate transient groundwater conditions during 1997–2006, followed by validation (2007–2011). The aquifer and aquitard hydraulic properties were obtained by pumping tests and were calibrated along with the rainfall recharge. The statistical and graphical performance indicators suggested a reasonably good simulation of groundwater flow over the study area. Sensitivity analysis revealed that groundwater level is most sensitive to the hydraulic conductivities of both the aquifers, followed by vertical hydraulic conductivity of the confining layer. The calibrated model was then employed to explore groundwater-flow dynamics in response to changes in pumping and recharge conditions. The simulation results indicate that pumping has a substantial effect on the confined aquifer flow regime as compared to the unconfined aquifer. The results and insights from this study have important implications for other regional groundwater modeling studies, especially in multi-layered aquifer systems.     

Application of MODFLOW for simulating groundwater flow in the Trifilia karst aquifer, Greece

Karst aquifers can have a complex flow as a result of the formation of large conduits from dissolution features. As a result, a three-dimensional finite-difference groundwater flow model (equivalent porous media) may not apply as the dual porosity nature of karst features and the effects of turbulent flow cannot be directly simulated. Statistical analysis of karst hydrographs of the Trifilia aquifer in Greece showed the existence of a slightly karstified mass with high primary porosity that regulates the flow. An equivalent porous media model was developed to simulate the Trifilia karst aquifer using MODFLOW. Steady state and transient state calibration gave encouraging results for the equivalent porous media approach, which does not consider pipe flow or turbulence. Detailed hydrogeological research conducted in the area helped define the aquifer hydraulic conductivity zones and extent; and flux to/from the aquifer. Only hydraulic conductivity and specific yield were adjusted during calibration, as the flux to/from the system was considered known and applied as boundary conditions. Small mean absolute and RMS piezometric head error of the model under both steady and transient state conditions were achieved.