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

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

Regional groundwater flow model for Abu Dhabi Emirate: scenario-based investigation

Despite the continuous increase in water supply from desalination plants in the Emirate of Abu Dhabi, groundwater remains the major source of fresh water satisfying domestic and agricultural demands. Groundwater has always been considered as a strategic water source towards groundwater security in the Emirate. Understanding the groundwater flow system, including identification of recharge and discharge areas, is a crucial step towards proper management of this precious source. One main tool to achieve such goal is a groundwater model development. As such, the main aim of this paper is to develop a regional groundwater flow model for the surficial aquifer in Abu Dhabi Emirate using MODFLOW. Up to our knowledge, this is the first regional numerical groundwater flow model for Abu Dhabi Emirate. After steady state and transient model calibration, several future scenarios of recharge and pumping are simulated. Results indicate that groundwater pumping remains several times higher than aquifer recharge from rainfall, which provides between 2 and 5% of total aquifer recharge. The largest contribution of recharge is due to subsurface inflow from the eastern Oman Mountains. While rainfall induced groundwater level fluctuation is absent in the western coastal region, it reaches a maximum of 0.5 m in the eastern part of the Emirate. In contrast, over the past decades, groundwater levels have declined annually by 0.5 m on average with local extremes spanning from 93 m of decline to 60 m of increase. Results also indicate that a further decrease in groundwater levels is expected in most of Emirate. At other few locations, upwelling of groundwater is expected due to a combination of reduced pumping and increased infiltration of water from nonconventional sources. Beyond results presented here, this regional groundwater model is expected to provide an effective tool to water resources managers in Abu Dhabi. It will help to accurately estimate sustainable extraction rates, assess groundwater availability, and identify pathways and velocity of groundwater flow as crucial information for identifying the best locations for artificial recharge.     

苏北沿海三市三维地下水流数值模拟

20世纪80年代以来,苏北沿海三市（连云港、盐城和南通）地下水过量开采导致部分地区水位持续大幅下降,引发地面沉降、水质咸化等问题。本文基于研究区大量地质、水文地质资料,以及2005—2013年地下水开采和水位观测数据,经过模型识别、验证,建立了能刻画苏北沿海三市地下水流场演化的三维非均质各向异性地下水流数值模型,模拟了研究区2005—2013年开采量由约2.03亿m³逐渐减少到1.56亿m³条件下深部地下水系统流场（包括第Ⅱ、Ⅲ、Ⅳ承压含水层）的动态演化过程。结果表明,各观测孔的水位实测值与模拟值吻合良好,水位及水均衡模拟结果显示各含水层水位下降得到有效控制,但是该区域地下水仍处于超采状态。为进一步控制水位下降,防治咸水扩散、地面沉降等地质环境灾害,模型预测并比较了现状开采和限制开采2种不同开采方案下2014—2020年苏北沿海三市深层地下水流场的变化趋势。预报结果显示,2种开采方案下,各含水层水位下降速率均减小,尤其是限制开采方案下,第Ⅱ、Ⅲ和Ⅳ承压含水层的区域水位平均下降速率依次为0.15、0.16和0.15 m/a,分别比模拟期减小了46.43%、65.21%和48.28%,在开采强度降至最低的条件下储存量变化速率仍为负值,说明需要进一步加强地下水限制开采的力度。     

Impact of the unpredictable changes in agricultural groundwater use on long-term sustainable exploitation: the case study of Bouhefna-Haffouz aquifer system,central Tunisia

Long-term exploitation schemes in many regions are often based solely on hydrodynamic factors, while the agricultural use of groundwater undergoes significant changes over time. The Bouhefna-Haffouz aquifer system in central Tunisia is one of those cases where an aquifer exploitation scheme was well designed hydrodynamically to address the political needs at the time. Fifty years later, a numerical groundwater model has been conducted to assess the sustainability of the scheme. Results show that the scheme aimed to lower the groundwater level to reduce overflow to Merguellil Wadi and maintain it at a level that benefits agricultural profitability. This caused loss of the Merguellil baseflow, forcing farmers to switch from traditional irrigation canals to deep wells and motor pumps, thereby disrupting the hydrological budget even further. The numerical model indicates that the flow to the wadi reached zero in 1978, the average flow by vertical leakage decreased from 8 hm³ in 1970 to 2 hm³ in 2020, and the horizontal percolation between the regional aquifer units increased from 1 hm³ in 1970 to 6 hm³ in 2020. Although the groundwater exploitation scheme was not previously considered a factor in local hydrological changes, the results of this study demonstrate the significant impact of societal behavior following the scheme’s implementation on the hydrological budget of Merguellil Wadi.

     

Characterization of the groundwater flow system in the hillside and coastal aquifers of the Mersin-Tarsus region (Turkey)

In the region between Mersin and Tarsus cities, located along the Mediterranean Sea coast in southern Turkey, the demand for groundwater has increased dramatically as the available surface water supplies have already been developed. Fundamental information is required to characterize the existing groundwater system in this area in order to establish a sustainable groundwater-use policy. For this purpose, hydrochemical and environmental isotopic data were collected and integrated with available geological and hydrogeological information to develop a conceptual model of the system. Results, backed up mainly by depleted stable isotope composition and infinitesimal tritium content, suggest that most of the groundwater along the coastal zone is supplied by the neighboring mountain belt while local precipitation has also contributes to aquifer recharge. The validation of the conceptual perspective by a steady-state numerical groundwater flow model reveals that about 90% of the recharge to the aquifer system is supplied by the deep flow of karstic groundwater fed from the Taurus Mountains. Monitoring of changes in the recharge regime of the mountain sector seems to be critical in establishing future groundwater use policies.     

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.     

海潮引起初始水位倾斜的滨海承压含水层地下水位 变化的数值模拟

海潮波动可以引起海岸带地下水位发生波动。建立了基于有限差分法的滨海地区一维承压含水层地下水运动数值模型。将潮汐波动概化为正弦波,模拟了滨海地区地下水位随潮汐波动的变化。通过与初始水位水平的承压含水层水位变化的比较表明,受海潮影响的滨海承压含水层地下水位与海潮有相似的波动特征,但变幅减小,受海潮的影响程度与离海岸的距离有关,随着离海岸距离的增加,地下水位的变幅及潮汐效率呈负指数函数衰减,但比前者变化程度稍缓,地下水位对海潮的滞后时间随距离呈线性增加。     

连续墙周围的地下水渗流特征及数值模拟

本文以上海市彭越浦泵站降水工程为例,分析了连续墙周围地下水的渗流特征:上部地下水不连续,底部地下水连续,呈三维流态。根据场地的工程地质与水文地质条件,提出了粘弹性越流含水层组三维修流数学模型,并进行了数值模拟。计算的结果基本上反映了连续墙周围地下水三维流特征。     

Numerical simulation of the effect of heavy groundwater abstraction on groundwater–surface water interaction in Langat Basin, Selangor, Malaysia

The paper aims at evaluating the interaction between ground and surface water along the Langat River in Malaysia through the development of a numerical simulation. Malaysia has been experiencing a rapid economic growth since the last few decades, driven by many factors such as agriculture, industry, and the like. The demand for water in these sectors has increased so tremendously that surface water has been utilized in conjunction to groundwater. Approximately 18,184 m³ of water per day is obtained from the aquifer to supply to the steel factory. There are also workshops, petroleum stations, and houses in the area thus causing the water quantity and quality to degrade. In terms of quantity, the pumping activity has altered the interaction between the groundwater and surface water. Therefore, a numerical model was proposed and two aquifer layers were simulated, with the first layer being approximately >20 m in depth and the second layer >100 m. The recharge estimated from the tank model was input into the groundwater modeling. The effects of the surface water to the aquifer were included in the simulation by defining the river conductance, river bed, and river level. The calibrated model (error about 0.9 m) was achieved and applied to predict the flow pattern in its natural state without the pumping and with the pumping states. As a result, in the first scenario, the stream was in an effluent condition influenced by the groundwater from the northeast to the west. A hyporheic flow occurred and was observed from the contour map. The flow system was changed in the second scenario when the pumping activity was included in the simulation. The groundwater lost its original function but received leakage from the stream near the pumping sites. The findings of this study will help the local authorities and other researchers to understand the aquifer system in the area and assist in the preparation of a sustainable groundwater management.     

Influence of flooding on groundwater flow in central Cambodia

Cambodia is affected by flooding from the Mekong, Tonle Sap and Bassac rivers every year, which harms human populations and damages property, as well as alters the water quality in aquifer systems. The objective of this paper is to highlight the effects of river flooding on groundwater flow using numerical simulation. A two-dimensional groundwater flow model coupled with a groundwater recharge model was applied to the research area in central Cambodia. River level variation was included in model processes, and flood areas and periods were assigned. The results showed that during flooding periods, floodwater from the three rivers played an important role in recharging groundwater. During the dry season, Tonle Sap River received groundwater supply from the northwest, and levels in the Bassac and Mekong River dropped to lower than the groundwater level. This study improves understanding of the surface water and groundwater flow system in the study area.     

Numerical simulation of sea water intrusion near Beihai, China

 A leaky aquifer system occurs in the coastal plain near Beihai, China. Seawater intrusion into the confined aquifer took place along the northern coast. Chloride concentrations at some observation wells increased steadily from 1988 and were at their peak in 1993. A quasi-three-dimensional element model has been developed to simulate the spatial and temporal evolution of hydraulic heads and chloride concentrations of the groundwater near the northern coast. The simulation model was based on the transition zone approach, which requires simultaneous solution of the governing water flow and solute transport equations. An irregular grid of a quadrangle was used to discretize the flow domain. Various aquifer parameters were verified with the numerical model in order to obtain satisfactory matches between computed values and observed data from an investigation. Three pumpage schemes were designed to use the calibrated model for prediction of future changes in water levels and chloride concentrations in groundwater in the study area. Results show that seawater intrusion would worsen in the confined aquifer if the current rates of groundwater pumpage continue. The alternative, to eliminate pumpage in the intruded area and to moderate pumpage rates from water supply wells far from the seashore, may limit seawater intrusion significantly and is considered attractive in the area. Received: 27 September 1999 · Accepted: 27 December 1999     

天津市地下水流-地面沉降耦合模型

天津市平原区地面沉降主要由地下水大量开采引起,影响范围广、危害大,已成为天津市主要的环境地质问题。分析了研究区的水文地质条件,结合地下水开发利用状况,将研究区概化为6个含水层组,地下水流考虑三维非稳定流,地面沉降选用一维固结压缩模型,运用地下水流模型Modflow 2005和地面沉降模拟模块 Sub,建立了天津市平原区地下水流-地面沉降数值耦合模型,模型面积为1.1×10⁴ km²,利用1998-2008年地下水位等值线、过程线、地面沉降过程线等资料对模型进行了识别。模拟期的地下水均衡分析表明,在多年开采条件下,越流补给、压缩释水、侧向边界流入分别占深层含水层补给量的41.84％、32.15％和24.17％。将调试后的模型应用于南水北调实施后地下水控采条件下的地面沉降趋势预测,显示出停采或减少地下水的开采,有利于减缓地面沉降下降速度,且表现出开采层位越往下,地面沉降恢复难度越大的变化趋势。     

含水岩组概化的累积导水系数法

地下水流系统理论和数值模拟技术分别是水文地质学的基本理论和技术方法,含水岩组的概化是地下水流系统分析和地下水数值模拟的重要基础,直接影响着数值模拟和水流系统分析的精度和可信度.为提高含水岩组概化的精度和可信度,提出一种含水岩组概化的新方法,即累积导水系数法.依据岩层厚度与渗透系数乘积累积值随深度的变化,以及水文地质剖面岩性分布的整体特征,概化含水介质结构.以玛纳斯河流域为例,应用该方法概化流域内的岩性剖面,结合GMS软件中TINS模块构建水文地质结构模型.结果表明,应用该方法概化后的含水层结构具有较好的合理性和仿真性,建立的三维模型很好地显示了研究区含水介质的空间展布特征,为建立地下水流模型奠定了良好的基础.      

Deterministic modelling of the cumulative impacts of underground structures on urban groundwater flow and the definition of a potential state of urban groundwater flow: example of Lyon,France

Underground structures have been shown to have a great influence on subsoil resources in urban aquifers. A methodology to assess the actual and the potential state of the groundwater flow in an urban area is proposed. The study develops a three-dimensional modeling approach to understand the cumulative impacts of underground infrastructures on urban groundwater flow, using a case in the city of Lyon (France). All known underground structures were integrated in the numerical model. Several simulations were run: the actual state of groundwater flow, the potential state of groundwater flow (without underground structures), an intermediate state (without impervious structures), and a transient simulation of the actual state of groundwater flow. The results show that underground structures fragment groundwater flow systems leading to a modification of the aquifer regime. For the case studied, the flow systems are shown to be stable over time with a transient simulation. Structures with drainage systems are shown to have a major impact on flow systems. The barrier effect of impervious structures was negligible because of the small hydraulic gradient of the area. The study demonstrates that the definition of a potential urban groundwater flow and the depiction of urban flow systems, which involves understanding the impact of underground structures, are important issues with respect to urban underground planning.     

Integrated approach of geospatial visualization and modeling for groundwater management in hard rock terrains in Nagpur Urban Area,India

Groundwater management requires complete visualization of aquifer characteristics to understand scientific aspects and hence remains a challenge, especially in hard rock terrains. In the present research paper, a comprehensive approach using 3D stratigraphic model, fence diagrams and numerical modeling has been proposed to understand the groundwater status for effective recharge. The stratigraphy, groundwater flow, and groundwater fluctuations for the period 1999–2010 were analyzed. The total volume of formations, volume of voids, storage capacity, and quantities of recharge of unconfined aquifer system in the Nagpur urban area were estimated. The steady state groundwater flow model of Basalt formation was calibrated to evaluate the subsurface system using Processing Modflow (PMWIN 5.3.2). The calibrated hydraulic head is compared with field observed head. The comparative spatial analysis presents a simple integrated approach in identifying zones with falling groundwater trends suitable for groundwater recharge in hard rock terrain in Nagpur urban area.     

Numerical simulation of the influence of management alternatives of a projected reservoir on a small alluvial aquifer affected by seawater intrusion (Almuñécar, Spain)

The present study concerns the application of a numerical approach to describe the influence of anthropogenic modifications in surface flows (operation of a projected reservoir) on the freshwater-seawater relationships in a downstream coastal aquifer which has seasonal seawater intrusion problems (River Verde alluvial aquifer, Almuñécar, southern Spain). A steady-state finite element solution to the partial differential equation governing the regional motion of a phreatic surface and the resulting sharp interface between fresh water and salt water was used to predict the regional behavior of the River Verde aquifer under actual surface flow conditions. The present model approximates, with simple triangular elements, the regional behavior of a coastal aquifer under appropriate sinks, sources, Neumann and open boundary conditions. A steady-state solution to this numerical approach has been shown to precisely calculate freshwater heads, saltwater thicknesses, and freshwater discharges along steeply sloping coasts. Hence, the adequate treatment and interpretation of the hydrogeological data which are available for the River Verde aquifer have been of main concern in satisfactorily applying the proposed numerical model. Present simulated conditions consider steady-state yearly averaged amounts of external supplies of fresh water in order to determine the influences of the projected Otívar reservoir on the further behavior of the River Verde coastal aquifer. When recharges occur at the coastline, essentially because of freshwater deficits due to groundwater overexploitation, a hypothesis of mixing for the freshwater-saltwater transition zone is made in order to still allow the model to continue calculating groundwater heads under the sea level, and, as a consequence, the resulting seawater intrusion and recharges of saltwater from the sea. Simulations show that a considerable advance in seawater intrusion would be expected in the coastal aquifer if current rates of groundwater pumping continue and a significant part of the runoff from the River Verde is channeled from the Otívar reservoir for irrigation purposes.     

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.     

基于数值模拟的区域性地下水补给过程实例研究

研究主要目的是以数值模拟的方法,再现选定的研究区域20世纪70~90年代地下水位处于较低水平的过程并推求地下水的补给过程。研究区为某流域入海口前10×12km2的矩形区域,基于研究区土层纵断面的物理条件,在平面二维的基础上,考虑含水层垂向水收支成分,利用地下水计算的基本方程式开发地下水准三维数值计算模型,以数值模拟的方法验证模型效率。对地下水位数值模拟的结果表明,模型效率可以达到80%以上。基于1966~2005年的数值计算结果,推求了研究区各含水层的补给过程。研究结果可为区域性地下水数值模拟技术及地下水资源量的评估等研究提供方法上的借鉴。     

海潮引起有越流的滨海承压含水层地下水头变化的数值模拟

海潮波动可以引起海岸带有越流的承压含水层地下水头发生波动。建立了基于有限差分法的滨海地区一维承压含水层地下水运动数值模型。通过将潮汐波动概化为正弦波,分别对初始水头水平及线性倾斜的承压含水层模拟了滨海地区有越流的承压含水层地下水头随潮汐波动的变化。通过对两种情形下的变化比较,结果表明,受海潮影响的滨海承压含水层地下水头与海潮有相似的波动特征,但变幅减小,受海潮的影响程度与离海岸的距离有关,随着离海岸距离的增加,地下水头的变幅及潮汐效率呈负指数函数衰减,水头倾斜情形下变幅更小,潮汐效率更小,滞后时间更短,地下水头对海潮的滞后时间随距离呈线性增加。     

Groundwater flow modelling in the central zone of Hanoi, Vietnam

On the basis of a review of the Quaternary sedimentary architecture in the area of Hanoi city (Vietnam), a numerical, deterministic and three-dimensional groundwater flow model has been built for a simulation between 1995 and 2004. The sedimentary architecture has been constructed on the basis of the data from 32 drillings covering the entire Quaternary sequence (but with little sedimentological detail), as well as hydrographical and hydrogeological data. Both steady- and transient-state conditions were tested. Results calculated by the model seem to indicate that the conceptual hypotheses adopted are reasonable for the modelled domain and period. The simulation allows for calculation of the regional groundwater flow trends. It is also used for assessing the relative importance of the various recharge sources of the shallow aquifer system in Hanoi, and for estimating the interactions between groundwater and the Red River.