大区域地下水流数值模拟研究现状及存在问题

在分析大区域地下水流数值模型构建缘起的前提下,系统论述了近年来地下水流数值模拟在大区域地下水资源评价、水文地质参数确定、地面沉降、溶质运移、海水入侵、盐渍化、风险评估、地下水管理及地表水与地下水的联合开发利用等方面的国内外研究应用现状;归纳、总结了目前大区域地下水流数值模型在灵敏度分析、裂隙和岩溶介质中模型建立、基于地下水流数值模拟的溶质运移模型建立、地下水流数值模型构建所需工作量等理论和方法研究及实际建模过程中存在的一些问题;展望了今后大区域地下水流数值拟在研究范围、模拟技术与方法以及与其它模型耦合等方面的发展趋势。      

某应急水源地地下水资源评价

采用地下水模拟系统软件GMS6.0,建立评价区地下水数值模拟模型。通过地下水流模拟,直观地反映并预测出开采条件下地下水降落漏斗的形态。     

某市地下水的模拟计算与分析

在详细分析研究区地质和水文地质条件的基础上,对研究区的均衡要素、流场及参数分布进行描述。建立水文地质概念模型,采用基于有限元方法的FEFLOW软件建立地下水流数值模拟模型并进行模型的识别和验证,在此基础上建立了地下水溶质运移模型,并对研究区地下水中的污染物运动进行数值模拟研究和分析。模拟结果表明：研究区西北部地区污染严重,原因主要是垃圾填埋场,建议有关部门采取相应的治理措施。     

陈崇希教授的学术思想和成就综述

在地下水资源评价理论方面, 陈崇希教授分析了"平均布井法"不符合质量守衡原理的实质, 纠正了以"地下水补给量计算可持续开采量"的错误, 提出了基于"质量守衡"的地下水资源评价原则, 强调分析"补给的增加量与排泄的减少量"在评价地下水可持续开采量时的重要意义.在地下水动力学领域, 陈崇希教授纠正了稳定井流"影响半径"模型的错误, 恢复了Dupuit"圆岛模型"的原貌, 拓展了Theis公式和Hantush公式的应用条件, 改进了地下水非稳定井流理论, 完善了其中的某些基本概念.在水文地质模拟仿真技术方面, 陈崇希教授提出确定滨海承压含水层海底边界的理论和方法; 提出地下水混合井流的模型和模拟方法, 解决了混合抽水试验确定分层水文地质参数的难题; 提出岩溶管道-裂隙-孔隙三重介质的地下水线性-非线性流动的模型; 建立了考虑井管水流雷诺数对滤管入流量分布的水平井-含水层系统的地下水流模型; 完成了"渗流-管流耦合模型"的砂槽物理模拟, 并用数值方法仿真了地下水流的规律; 最近向观测孔水位形成的传统观念———常规观测孔中的水头降深反映该孔滤水管中各点的平均降深———提出质疑.陈崇希教授建立的"渗流-管流耦合模型"使传统的基于线汇/线源的井孔-含水系统模型提高到新的水平.陈崇希教授积极倡导"防止模拟失真, 提高仿真性", 强调精细地分析水文地质条件、合理地概化模型和采用正确的仿真技术的重要性.      

长春市城区基岩裂隙水开采优化决策研究

论文阐述了长春市城区基岩裂隙水开采区地下水位、水质状态的变化特征。在研究区地下水流数值模拟的基础上建立地下水资源优化管理模型,确定了地下水优化开采方案,并提出实现地下水合理开采、优化管理的决策。     

改进的遗传算法在地下水数值模拟中的应用

地下水流数值模拟中的模型识别问题，可以转化为函数的最优化问题。鉴于遗传算法的特点，将之引入到地下水流数值法中，用以解决地下水数学模型的识别问题。在建立地下水数值模拟中模型识别问题的是优化模型后，采取将最优化模型中的目标函数嵌入到遗传算法适应度函数中的方法，实现遗传算法与地下水流数值法的耦合。基于优化模型和遗传算法的运算过程，编写计算程序，实现地下水数学模型的自动识别。根据在珲春盆地地下水资源评价实例中应用得到的结果，信纸证了改进的遗传算法在地下水数值模拟中应用的可行性与有效性。     

石家庄市东部平原FEFLOW地下水数值模拟与研究

根据研究区水文地质条件、开采条件及资料条件,将研究区内水文地质概念模型概化为二维非均质、各向同性、非稳定地下水流系统.利用先进的FEFLOW有限元地下水模拟系统建立了研究区的区域地下水流模型,预测不同开采方案下的地下水动态变化,为地下水资源开发利用和优化配置提供科学依据.     

地下水动态观测网优化设计研究

本文采用状态空间分析法,通过对地下水流系统输入变量、输出变量及地下水动态测量值的研究,提出了地下水流系统确定一随机性数值模型。运用该模型的模拟解,进行地下水动态观测网优化设计。该方法把地下水动态观测网密度、位置及观测频率优化与地下水流系统结合起来研究。为了克服该模型计算量大的缺陷,本文提出了一种改进算法,使这一观测网优化方法适于在微机上实现,减少了计算工作量。该方法曾用于陕西某地区地下水位动态观测网的优化设计,取得了显著经济效益。     

基岩岛屿地下水流场数值模拟研究

岛屿是我国领土的重要组成部分,对国家安全和国防军事意义非凡。岛屿的地下水资源尤为珍贵,地下水流场模拟是研究地下水分布规律的重要手段。水文地质条件复杂、可利用的观测井较少等原因,造成了基岩岛屿水文地质模型概化精度不高、初始条件难以获取等问题。为克服基岩岛屿地下水流场模拟的诸多困难,将珠海外伶仃岛作为研究区,利用数字高程模型数据开展地质地貌建模; 利用探地雷达法、直流电法与地质分析对岛屿进行探测,获取了地层数据; 采用地下水遥感评估法,利用实测井位数据,确定了地下水的初始水位,进而对基岩岛屿地下水流场进行建模; 最终,通过对外伶仃岛地下水流场的数值模拟得出地下水模拟流场图。岛上多个测点的探测水位值与模拟水位相关性较好,其拟合优度R²为0.872 2。由此可见,综合遥感、物探、水文地质手段等技术方法获取的数据,采用地下水模拟软件或程序实现基岩岛屿地下水流场的数值模拟,是基岩岛屿地下水资源研究的一个有效方法。     

高放废物处置库北山预选区区域水文地质概念模型

地下水流数值模拟是高放废物处置库选址中水文地质评价的重要内容.水文地质概念模型是地下水流数值模拟的基础和前提.在北山预选区水文地质分析基础上,利用DEM数据生成了研究区地表水系,据此确定了模拟范围;借助GMS软件,建立了研究区三维水文地质结构模型;再通过边界条件、流场特征、地下水均衡项和水文地质参数的分析和描述,建立了北山地区水文地质概念模型.     

Application of MODFLOW and geographic information system to groundwater flow simulation in North China Plain, China

MODFLOW is a groundwater modeling program. It can be compiled and remedied according to the practical applications. Because of its structure and fixed data format, MODFLOW can be integrated with Geographic Information Systems (GIS) technology for water resource management. The North China Plain (NCP), which is the politic, economic and cultural center of China, is facing with water resources shortage and water pollution. Groundwater is the main water resource for industrial, agricultural and domestic usage. It is necessary to evaluate the groundwater resources of the NCP as an entire aquifer system. With the development of computer and internet information technology it is also necessary to integrate the groundwater model with the GIS technology. Because the geological and hydrogeological data in the NCP was mainly in MAPGIS format, the powerful function of GIS of disposing of and analyzing spatial data and computer languages such as Visual C and Visual Basic were used to define the relationship between the original data and model data. After analyzing the geological and hydrogeological conditions of the NCP, the groundwater flow numerical simulation modeling was constructed with MODFLOW. On the basis of GIS, a dynamic evaluation system for groundwater resources under the internet circumstance was completed. During the process of constructing the groundwater model, a water budget was analyzed, which showed a negative budget in the NCP. The simulation period was from 1 January 2002 to 31 December 2003. During this period, the total recharge of the groundwater system was 49,374 × 10⁶ m³ and the total discharge was 56,530 × 10⁶ m³ the budget deficit was −7,156 × 10⁶ m³. In this integrated system, the original data including graphs and attribution data could be stored in the database. When the process of evaluating and predicting groundwater flow was started, these data were transformed into files that the core program of MODFLOW could read. The calculated water level and drawdown could be displayed and reviewed online.     

GIS assisted inverse geochemical modeling for plausible phase transfers in aquifers

A geographical information system (GIS) assisted approach that couples a groundwater flow model and an inverse geochemical model is presented to quantify the phase mole transfers between two points on the flow path within a groundwater system. It is used to investigate the plausible phase transfers in the unconfined aquifer of Mehsana district of Gujarat State, India. X-ray diffraction analysis of soil samples is carried out for mineral characterization. The groundwater flow field is simulated using MODFLOW and flow paths used for inverse geochemical modeling are traced using PMPATH, a particle tracking algorithm. The plausible phase mole transfers in the flow path are quantified using PHREEQC geochemical code. The different scenarios generated by inverse modeling routine are used as input to the geochemical model and simulation runs are taken as forward models. The obtained results are compared to the target solution chemistry by using the square of the Pearson product moment correlation coefficient. Results reveal that the groundwater is undersaturated with anhydrite, carbon dioxide (gas), fluorite, gypsum, halite, jarosite-K, and siderite. It is oversaturated with aragonite, calcite, dolomite, ferrihydrite, goethite, and hematite. Results further reveal that calcite is precipitating, while dolomite, gypsum, carbon dioxide, and fluorite are dissolving together with ion exchange in the flow path.     

虚拟地理信息技术的地下水动态模拟分析应用

针对现有地下水流模拟软件在数据预处理、三维可视化效果等方面的不足,以Visual C⁺⁺ 6.0为主要开发环境、OpenGL为三维图形显示与渲染工具,以地下水及其赋存地质体观测记录为主要数据源,实现了数据处理、建模、存储、调度、查询以及真三维可视化,在此基础上完成了地下水动态模拟分析功能,实现了虚拟地理信息技术在地下水工作中的系统应用,具有功能开放、丰富、灵活的特点。     

The impact of tunneling construction on the hydrogeological environment of “Tseng-Wen Reservoir Transbasin Diversion Project” in Taiwan

Groundwater flow and the associated surface water flow are potential negative factors on underground tunnels. Early detection of environmental impacts on water resources is of significant importance to planning, design and construction of tunnel projects, as early detection can minimize accidents and project delays during construction. The groundwater modeling software package Groundwater Modeling System (GMS), which supports the groundwater numerical codes MODFLOW and FEMWATER, was utilized to determine the impact of tunneling excavation on the hydrogeological environment in a regional area around the tunnel and a local hot springs area, at the “Tseng-Wen Reservoir Transbasin Diversion Project”, in Taiwan. A hydrogeological conceptual model was first developed to simplify structures related to the site topography, geology and geological structure. The MODFLOW code was then applied to simulate groundwater flow pattern for the hydrogeological conceptual model in the tunnel area. The automated parameter estimation method was applied to calibrate groundwater level fluctuation and hydrogeological parameters in the region. Calibration of the model demonstrated that errors between simulated and monitored results are smaller than allowable errors. The study also observed that tunneling excavation caused groundwater to flow toward the tunnel. No obvious changes in the groundwater flow field due to tunnel construction were observed far away in the surrounding regions. Furthermore, the FEMWATER code for solving 3-D groundwater flow problems, in which hydrogeological characteristics are integrated into a geographic information system (GIS), is applied to evaluate the impact of tunnel construction on an adjacent hot spring. Simulation results indicated that the groundwater drawdown rate is less than the groundwater recharge rate, and the change to the groundwater table after tunnel construction was insignificant for the hot spring area. Finally, the groundwater flow obtained via the GMS indicated that the hydrogeological conceptual model can estimate the possible quantity of tunnel inflow and the impact of tunnel construction on the regional and local groundwater resources regime of the transbasin diversion project.     

Using lithologic modeling techniques for aquifer characterization and groundwater flow modeling of the Sohag area, Egypt

Three dimensional lithologic modeling techniques have been used for detailed characterization and groundwater flow modeling of the Quaternary aquifer system of the Sohag area, Egypt. Well log data were used for building the lithologic model using RockWorks. A groundwater flow model, facilitated by MODFLOW 2000, was built using results of the lithologic model. The obtained lithologic model honored the well log data and revealed a complex sedimentary system, which is mainly composed of six lithologic categories: clay, clay and sand, fine sand, coarse sand, sand and gravel, and gravel. Inter-fingering and presence of lenses are the main characteristics of the sedimentary basin in the study area. A wide range of hydraulic conductivities is present, which vary spatially and control the groundwater flow. Heterogeneity of the aquifer system is spatially represented where different hydraulic conductivities are found in the different directions. Sandy layers tend to be connected. Hydraulic continuity is represented by inter-fingering and connection of sandy materials within the aquifer system.     

地表水文过程与地下水动力过程耦合模拟及应用

地表水和地下水是水循环系统不可缺少的部分,对相互作用与反馈的模拟计算有助于水资源精准评价。在非饱和带数值计算模型(UZF1)与模块化地下水动力模型(MODFLOW2005)耦合的基础上构建了一个地表水文过程与地下水动力过程耦合数值模拟模型,模型可以同时输出流域出口断面流量和地下水位的空间分布,有利地提高了模型的可靠性。将模型应用于资料条件较好的临涣集流域,从模拟结果分析来看,模型可以很好地模拟流域出口流量与地下水位空间动态变化过程,具有较高的模拟精度。     

基于MapObjects的C/S管线信息系统的设计与实现

基于面向对象的地理信息系统技术，针对客户机／服务器(C／S)模式设计了一个包含数据库、应用服务器和客户端三个部分的管网管线地理信息系统。该系统由测量调查、图形编辑、数据统计、信息查询、空间分析、规划、资料输出、权限管理、系统维护九大功能模块组成；并具有用户界面友好，稳定实用，安全可靠，可扩充性强等特点。并在数据存储方式、图形数据与属性数据的连接、编辑功能设计、系统兼容性、符号库建立等关键技术的实现方面有新的突破。该系统除具有一般管理系统的基本功能外，还具有地理信息系统所特有的空间数据查询和分析功能；该系统还充分发挥了网络的优势,为CIS信息社会化共享提供了一个切实可行的解决方案。     

GIS Spatial-Temporal Modeling of Water Systems in Greater Toronto Area, Canada

Modeling landscape with high-resolution digital elevation model (DEM) in a geographic information system can provide essential morphological and structural information for modeling surface processes such as geomorphologic process and water systems. This paper introduces several DEM-based spatial analysis processes applied to characterize spatial distribution and their interactions of ground and surface water systems in the Greater Toronto Area (GTA), Canada. The stream networks and drainage basin systems were derived from the DEM with 30 m resolution and the regularities of the surface stream and drainage patterns were modeled from a statistical/multifractal point of view. Together with the elevation and slope of topography, other attributes defmed from modeling the stream system, and drainage networks were used to associate geological, hydrological and topographical features to water flow in river systems and the spatial locations of artesian aquifers in the study area. Stream flow data derived from daily flow measurements recorded at river gauging stations for multi-year period were decomposed into “drainage-area dependent“ and “drainage-area independent“ flow components by two-step “frequency“ and “spatial“ analysis processes. The latter component was further demonstrated to relate most likely to the ground water discharge. An independent analysis was conducted to model the distribution of aquifers with information derived from the records of water wells. The focus was given on quantification of the likelihood of ground water discharge to river and ponds through flowing wells, springs and seepages. It has been shown that the Oak Ridges Moraine (ORM) is a unique glacial deposit that serves as a recharge layer and that the aquifers in the ORM underlain by Hilton Tills and later deposits exposed near the steep slope zone of the ridges of ORM provide significant discharge to the surface water systems (river flow and ponds) through flowing wells, springs and seepages. Various statistics (cross- and auto-correlation coefficients, fractal R/S exponent) were used in conjunction with GIS to demonstrate the influence of land types, topography and geometry of drainage basins on short- and long-term persistence of river flows as well as responding time to precipitation events. The current study has provided not only insight in understanding the interaction of water systems in the GTA, but also a base for further establishment of an on-line GIS system for predicting spatial-temporal changes of river flow and groundwater level in the GTA.     

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.     

GIS-based spatial and temporal prediction system development for regional land subsidence hazard mitigation

An integrated GIS-based approach for establishing a spatial and temporal prediction system for groundwater flow and land subsidence is proposed and applied to a subsidence-progressed Japanese coastal plain. Various kinds of fundamental data relating to groundwater flow and land subsidence are digitized and entered into a GIS database. A surface water hydrological cycle simulation is performed using a GIS spatial data operation for the entire plain, and the spatial and temporal groundwater infiltration quantity is hereby obtained. Through the data transformation from the GIS database to a groundwater flow code (MODFLOW), a 3D groundwater flow model is established and unsteady groundwater flow simulation for the past 21 years is conducted with results which compare satisfactorily with observed results. Finally, a Visual Basic code is developed for land subsidence calculations considering aquifer and aquitard deformation. Future land subsidence in the plain is predicted assuming different water pumping scenarios, and the results provide important information for land subsidence mitigation decision-making.