Effectiveness of time-series analysis for thermal plume propagation assessment in an open-loop groundwater heat pump plant

Thermal perturbation produced in the subsurface by open-loop groundwater heat pumps (GWHPs) represents a complex transport phenomenon that is affected by several factors, including intrinsic characteristics of the exploited aquifer, abstraction and reinjection well features, and the temporal dynamics of the accessed groundwater. Post-GWHP water may have become warmed or cooled before being reinjected into the aquifer, thereby creating a thermal plume, known as the thermal affected zone (TAZ), which can alter aquifer temperature. The TAZ is propagated mainly by advection, after which the plume tends to degrade via conductive heat transport and convection within moving water. Groundwater monitoring and multiparametric probes are used to check the dynamics of plume propagation and whether a system’s thermal plumes are generating unsuitable interference with wells, subsurface infrastructure, or land use. Analyses of time-series groundwater monitoring data can be used to monitor TAZ movement. In this paper, the thermal plume velocity was calculated by both an analytical solution and cross-correlation. Cross-correlation calculated between temperature measured in the reinjection well and control downstream piezometers can reveal plume dynamics and demonstrate the importance of advective transport in aquifer heat transfer.     

Thermal modeling of a Swiss urban aquifer and implications for geothermal heat pump systems

The progressive electrification of the building conditioning sector in recent years has greatly contributed to reducing greenhouse gas emissions by using renewable energy sources, particularly shallow geothermal energy. This energy can be exploited through open and closed shallow geothermal systems (SGS), and their performances greatly depend on the ground/groundwater temperature, which can be affected by both natural and anthropogenic phenomena. The present study proposes an approach to characterize aquifers affected by high SGS exploitation (not simulated in this work). Characterization of the potential hydro/thermogeological natural state is necessary to understand the regional flow and heat transport, and to identify local thermal anomalies. Passive microseismic and groundwater monitoring were used to assess the shape and thermal status of the aquifer; numerical modeling in both steady-state and transient conditions allowed understanding of the flow and heat transport patterns. Two significant thermal anomalies were detected in a fluvio-glacial aquifer in southern Switzerland, one created by river water exfiltration and one of anthropogenic nature. A favorable time lag of 110 days between river and groundwater temperature and an urban hot plume produced by underground structures were observed. These thermal anomalies greatly affect the local thermal status of the aquifer and consequently the design and efficiency of current and future SGS. Results show that the correct characterization of the natural thermo-hydrogeological status of an aquifer is a fundamental basis for determining the impact of boundary conditions and to provide initial conditions required to perform reliable local thermal sustainability assessments, especially where high SGS exploitation occurs.

     

Effects of regional groundwater flow on the performance of an aquifer thermal energy storage system under continuous operation

Numerical investigations and a thermohydraulic evaluation are presented for two-well models of an aquifer thermal energy storage (ATES) system operating under a continuous flow regime. A three-dimensional numerical model for groundwater flow and heat transport is used to analyze the thermal energy storage in the aquifer. This study emphasizes the influence of regional groundwater flow on the heat transfer and storage of the system under various operation scenarios. For different parameters of the system, performances were compared in terms of the temperature of recovered water and the temperature field in the aquifer. The calculated temperature at the producing well varies within a certain range throughout the year, reflecting the seasonal (quarterly) temperature variation of the injected water. The pressure gradient across the system, which determines the direction and velocity of regional groundwater flow, has a substantial influence on the convective heat transport and performance of aquifer thermal storage. Injection/production rate and geometrical size of the aquifer used in the model also impact the predicted temperature distribution at each stage and the recovery water temperature. The hydrogeological-thermal simulation is shown to play an integral part in the prediction of performance of processes as complicated as those in ATES systems.     

基于不同地质统计方法的渗透系数场对污染物运移的影响

渗透系数场的空间变异性是影响污染物运移结果的决定因素,而地质统计方法是解决渗透系数空间变异性的主要技术手段。本文利用野外场地实测数据,采用普通克里格法和指示克里格法、顺序高斯模拟法和顺序指示模拟法四种地质统计方法,插值估测和模拟再现随机渗透系数场,进而对比研究四种渗透系数场对大尺度污染物运移的影响。研究结果表明,污染羽的质心位置（一阶矩）主要由渗透系数的平均值来决定;污染羽在空间上的展布范围（二阶矩）主要受渗透系数空间变异方差的影响;条件模拟克服了估计法的平滑效果,较好地再现真实曲线的波动性,渗透系数（ lnK）估计方差与污染羽空间二阶矩随着条件模拟次数的增加而减小,并且顺序指示模拟程度更加明显。     

Impact of groundwater withdrawals on the interaction of multi-layered aquifers in the Viterbo geothermal area (central Italy)

The impact of groundwater withdrawals on the interaction between multi-layered aquifers with different water qualities in the Viterbo geothermal area (central Italy) was studied. In this area, deep thermal waters are used to supply thermal spas and public pools. A shallow overlying aquifer carries cold and fresh water, used for irrigation and the local drinking-water supply. Starting with a conceptual hydrogeological model, two simplified numerical models were implemented: a steady-state flow model of the entire groundwater system, and a steady-state flow and heat transport model of a representative area, which included complex interactions between the aquifers. The impact of increased withdrawals associated with potential future development of the thermal aquifer must be considered in terms of the water temperature of the existing thermal sources. However, withdrawals from the shallow aquifer might also influence the discharge of thermal sources and quality of the water withdrawn from the shallow wells. The exploitation of the two aquifers is dependent on the hydraulic conductivity and thickness of the intervening aquitard, which maintains the delicate hydrogeological equilibrium. Effective methods to control this equilibrium include monitoring the vertical gradient between the two aquifers and the residual discharge of natural thermal springs.     

A study of the contaminated banks of the Mahoning River, Northeastern Ohio, USA: characterization of the contaminated bank sediments and river water–groundwater interactions

The Mahoning River is one of the five most contaminated rivers in the U.S. This study characterized the contaminated sediments in the river banks and investigated the hydraulic interconnection between shallow aquifer in the banks with the river water. The study was conducted along the most polluted section of the river, which is 50-km long, using over 50 monitoring wells. The characterization part of the study investigated the sedimentology, hydraulic conductivity, and spatial distribution of the contaminated sediments. Results of the characterization revealed that the contaminated sediments consist of fine-grained sand, silt, mud, and clay. The spatial distribution of the contaminated sediment is heterogeneous and positively correlates with the hydraulic conductivity values, i.e., the greatest contamination occurs in high conductivity areas. Hydraulic conductivity was determined by the Hazen formula using 82 sediment samples. Bioremediation, which is one of the remedial options considered for the banks, is found to be hydraulically feasible because of sufficient hydraulic conductivity values (≥10^?4 cm/s) that ensure reasonable rates of nutrient delivery. Monitoring of water levels in the river and groundwater for a 10-month period shows that flow occurs from the river to groundwater and vice versa. The exchange of flow is influenced by rainfall. Flow of groundwater to the river will continually transport the dissolved contaminants in groundwater to the river. Therefore, findings of this study show that one of the remedial options that proposes dredging of channel sediments and permits no action for bank sediments cannot be chosen due to river water–groundwater interactions.     

Investigative approaches to determine exchange processes in the hyporheic zone of a low permeability riverbank

The hydraulic and hydrochemical interactions of surface water and groundwater in the riverbank of an alluvial stream in the Grand Duchy of Luxembourg were investigated. Hydraulic conductivity in the riverbank is low and the aquifer is confined. Gauged surface-water levels and groundwater levels were recorded over several months to observe potential effluent and influent aquifer conditions. Water samples were taken under different hydraulic conditions in the stream and from three observation wells in the riverbank. Samples were analysed for inorganic parameters and selected dissolved pharmaceuticals. The recorded groundwater levels in the riverbank responded almost without delay to changes in stream stage, which is very dynamic. Frequent changes from effluent to influent aquifer conditions were observed. The analyses for selected ions and pharmaceuticals support the assumption that the pharmaceutical compounds enter the riverbank via the stream only. The chosen trace organics are therefore suitable as anthropogenic tracers for groundwater and surface-water interactions at this field site. They prove that water exchanges also take place in riverbanks even where the hydraulic conductivity is low and for this reason indicate the existence of a hyporheic zone in the investigated section of the stream.     

地源热泵系统地下水热量运移阶段特性模拟研究

含水介质中地下水热量运移过程模拟是地源热泵系统这类新型浅层地温能利用工程科学论证的重要内容。通过建立饱和含水介质中热量运移对流-弥散数学模型,模拟了拟建地源热泵系统理想的对井抽-灌运行特定水流及热源条件下地下水热量运移过程。结果表明热泵系统完整运行周期(一年)内地下水热量运移过程具有显著的"阶段性"特征,"夏、冬抽-灌输运"和"春、秋停运蓄存"两阶段取决于不同的作用因素。从数学和概念模型两方面讨论了地下水流作用对热量输运过程的决定性影响,说明不同时期抽、灌井互换运行是解决拟定工作模式下热(冷)能浪费问题的理想途径。     

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.     

Synthesizing multifaceted characterization techniques to refine a conceptual model of groundwater sources to springs in valley settings (Minnesota,USA)

Springs are commonly used as low-cost monitoring locations to assess groundwater quality and long-term trends. However, spring waters in many settings are a mixture of groundwater sources that range in physical properties and water chemistry. The objective of this work was to determine water sources of springs emerging from the North American midcontinent Cambrian-Ordovician aquifer system at a fish hatchery near Lanesboro, Minnesota (USA), and compare and contrast the sources to shallower and deeper sources. The hydrology of the Lanesboro State Fish Hatchery has been studied for decades using a combination of dye tracing, thermal monitoring, geochemical sampling, and nearby borehole and outcrop observations. Previous studies are integrated with recently collected dye tracing results and geochemical data to develop a comprehensive conceptual model of groundwater flow. Dye trace findings and geochemistry indicate well-developed karst and bedrock fractures in shallowly buried unconfined carbonate formations are important transport pathways to convey anthropogenically influenced waters from the land surface to the hatchery springs. However, borehole dye traces, thermal monitoring, continuous nitrate monitoring, and mixing calculations show that a deeper confined siliciclastic aquifer is responsible for delivering relatively pristine water that accounts for about half of hatchery spring flux. Characterization of the hatchery’s groundwater systems provides fishery managers with information to protect this vital resource and improved context to interpret water-quality-monitoring data that track agricultural contaminants. The methods and results of this study may be widely applicable across a large extent of the Cambrian-Ordovician aquifer system, and to multiaquifer sedimentary bedrock systems elsewhere.

     

运用地下水对潮汐的响应识别压力传导系数

在地下水扩散方程中,压力传导系数是描述地下水运动的重要参数。传统的方法是通过抽水或注水给地下水系统一个扰动,监测地下水水位的响应,由此计算含水层的压力传导系数。文章提出用潮汐衰减率方法识别含水层压力传导系数,其适用于滨海区承压含水层的参数识别。在推导出解析解的基础上,通过数值拟合、最小二乘法、牛顿迭代法求得含水层的压力传导系数,提出潮汐衰减率的概念,建立余切函数与潮汐衰减率的线性关系,用线性关系中的斜率和截距识别压力传导系数。用潮汐衰减率方法识别出的压力传导系数与其实际值相等,说明该方法是正确有效的。潮汐信号衰减率与海水振荡的余切函数线性相关。数值仿真表明,该方法可以准确地估算出含水层的压力传导系数。潮汐衰减率方法具有少打井,经济高效等优点。潮汐衰减率方法为实际工程应用提供了可靠的理论基础,它可以用于部分实际工程中。该方法的局限性在于需要提供含水层的部分参数,如含水层的长度、海水波动振幅、频率等。     

Aquifer parameter estimation using tide-induced water-table fluctuations in the Biscayne Aquifer,Miami-Dade County,Florida (USA)

The Biscayne Aquifer (Florida, USA) is a coastal, shallow, unconfined, and heterogeneous aquifer with high water tables, composed of less-permeable sand to highly permeable karstic limestone. These properties make the Biscayne Aquifer one of the world’s most productive groundwater resources. The aquifer’s high yield and non-Darcian flow cause challenges for estimating aquifer parameters, which are essential for understanding groundwater processes and managing and protecting the groundwater resources. Water-table fluctuations in the Biscayne Aquifer are associated with astronomical tidal forces and gate operations on canal water-control structures. Analysis of observed groundwater level fluctuations can provide an understanding of the connectivity between the aquifer, Biscayne Bay, and the water level in the canals. Further, groundwater level fluctuations can be used for aquifer parameter estimation. In this research, observed ocean water levels measured at tidal stations and groundwater levels are fitted to Jacob’s analytical solution, where the amplitude of the groundwater head fluctuation decreases exponentially, and the time lag increases with distance from the shore. Observed groundwater levels were obtained from monitoring wells along the Miami-Dade shore and the barrier island of Miami Beach. Results indicate that Jacob’s solution is effective for aquifer parameter estimation in Miami Beach, where monitoring wells are closer to the shore. Estimated hydraulic conductivity appears to increase by four orders of magnitude to approximately 1 m s^–1 as the distance from shore increases. Constructing monitoring wells closer to the shore in Miami-Dade County and elsewhere would permit improved aquifer parameter estimation and support enhanced groundwater modeling efforts.

     

矿井水深层回灌过程量质耦合模拟分析

矿井水深井回灌是矿井水“转移存储”处理的主要形式,根据鄂尔多斯盆地煤矿区地质和矿井水特征,从回灌目的层地下水与矿井水的匹配性、上下岩层的隔水性、回灌层的渗透性以及封闭性角度提出了矿井水回灌目的层选取依据。并以地下水达西定律和Dupuit理论为基础,建立极坐标系完整注水井稳定流数学模型,得出在稳定注水条件下,回灌量与注水层渗透系数、厚度、回灌压力、水位埋深以及回灌井直径正相关,与影响半径负相关,与回灌层埋深无关。提出了矿井水深层回灌水动力和溶质运移耦合仿真模型构建方法,并以矿井水回灌试验案例为分析对象,模拟得出矿井水回灌过程中含水层水压形成以注水井为中心的“高位水丘”,且注水压力越大,回灌量增加较为明显,模型分析结果与现场试验结果基本一致。溶质运移范围形成以注水井为中心的“圆柱状”弥散形态,特征离子浓度沿回灌井两侧变化剧烈,回灌层特征离子浓度被迅速稀释,随着时间的延伸,弥散稀释范围增加相对较小,说明矿井水回灌对深部高浓度含水层地下水水化学影响程度不大,研究成果可为西部煤矿区矿井水高效回灌处理提供科学依据。      

Analytical solutions for transient temperature distribution in a geothermal reservoir due to cold water injection

An analytical solution to describe the transient temperature distribution in a geothermal reservoir in response to injection of cold water is presented. The reservoir is composed of a confined aquifer, sandwiched between rocks of different thermo-geological properties. The heat transport processes considered are advection, longitudinal conduction in the geothermal aquifer, and the conductive heat transfer to the underlying and overlying rocks of different geological properties. The one-dimensional heat transfer equation has been solved using the Laplace transform with the assumption of constant density and thermal properties of both rock and fluid. Two simple solutions are derived afterwards, first neglecting the longitudinal conductive heat transport and then heat transport to confining rocks. Results show that heat loss to the confining rock layers plays a vital role in slowing down the cooling of the reservoir. The influence of some parameters, e.g. the volumetric injection rate, the longitudinal thermal conductivity and the porosity of the porous media, on the transient heat transport phenomenon is judged by observing the variation of the transient temperature distribution with different values of the parameters. The effects of injection rate and thermal conductivity have been found to be profound on the results.     

类地行星热物理条件与热演化影响因素分析——以火星地幔热柱演化为例

太阳系内类地行星具有相似的岩石层包围金属核的圈层结构,在行星幔的热演化历史起源方面具有同时性和同源性,并且都在早期变形重力位能加热的基础上随放射性热能衰减而冷却。但是,由于半径、密度、粘度以及表层构造属性等物理条件的差异,其热演化历史各具特色。依据基本的热对流和热传导方程,我们计算分析了类地行星热物理条件差异对行星幔热演化历史的影响。计算表明,类地行星热演化的早期,行星幔热对流是主要的散热方式。半径较大的行星表面热流密度大,平均散热量也大。半径较小的行星内部温差小,粘滞系数高,对流能力低,提早进入传导散热状态,且传导散热的岩石层也比大行星厚。不同边界层热物理条件下,类地行星幔热演化历史会分别出现逐渐冷却的平稳式、包含热柱上涌的波动式、行星幔幕次翻转的周期式等特点不同的热演化过程。火星内部曾经存在的地幔热柱构造与火星地幔热动力学演化过程密切相关。我们从火星地幔热动力学演化模型出发,定量计算与地幔热柱构造演化相关的地幔热动力学演化特征,通过三维球壳数值模拟,研究了火星地幔热演化历史上可能存在的热柱活动造成的火星热演化历史的非单调变化,火星地幔对流环结构随时间的演变方式,以及与边界相关的地幔热柱对火星地形的影响。     

山东淄博市大武水源地裂隙岩溶水中污染物运移的数值研究

在分析研究淄博市大武水源地裂隙岩溶含水层的水力性质和污染物运移特征的基础上 ,对裂隙岩溶水的水头和污染物运移进行数值研究。目前国内外对裂隙岩溶水进行数值计算时 ,通常用等价多孔介质模型 ,但裂隙岩溶介质和多孔介质有很大不同。裂隙岩溶介质的储水和导水空间为裂隙网络 ,导水系数大 ,地下水的实际平均流速比孔隙水大得多 ,但给水度和贮水系数小。当用等价多孔介质模型进行模拟时应考虑这些特点。对于污染物运移的模拟 ,要同时求解水头方程和对流弥散方程 ,可采用MODFLOW和MT3D软件进行模拟。研究区裂隙岩溶水水头的数值计算表明 ,等效多孔介质模型水头的拟合误差能满足国标GB/T144 97- 93的要求。各时段地下水水量均衡计算的精度也满足要求。对流弥散方程的数值计算 ,由于Peclet数高达 95 .6 7,对流占绝对优势 ,可能存在数值弥散和数值振荡 ,因而采用多种方法进行了比较。对于同一问题 ,同时采用上游有限差分法 (UFDM) ,混合的欧拉拉格朗日方法 (特征线法MOC、改进特征线法MMOC和混合特征线法HMOC) ,总变异消减法(TVD)进行计算 ,并比较其结果。结果表明 ,混合特征线法 (HMOC)和总变异消减法 (TVD)比较适合于对流占优势的运移问题计算。由于渗透系数K和有效孔隙度θ对溶质运移结果的影响很大 ,?     

Characterization of the subsurface urban heat island and its sources in the Milan city area,Italy

Urban areas are major contributors to the alteration of the local atmospheric and groundwater environment. The impact of such changes on the groundwater thermal regime is documented worldwide by elevated groundwater temperature in city centers with respect to the surrounding rural areas. This study investigates the subsurface urban heat island (SUHI) in the aquifers beneath the Milan city area in northern Italy, and assesses the natural and anthropogenic controls on groundwater temperatures within the urban area by analyzing groundwater head and temperature records acquired in the 2016–2020 period. This analysis demonstrates the occurrence of a SUHI with up to 3 °C intensity and reveals a correlation between the density of building/subsurface infrastructures and the mean annual groundwater temperature. Vertical heat fluxes to the aquifer are strongly related to the depth of the groundwater and the density of surface structures and infrastructures. The heat accumulation in the subsurface is reflected by a constant groundwater warming trend between +0.1 and?+?0.4 °C/year that leads to a gain of 25 MJ/m² of thermal energy per year in the shallow aquifer inside the SUHI area. Future monitoring of groundwater temperatures, combined with numerical modeling of coupled groundwater flow and heat transport, will be essential to reveal what this trend is controlled by and to make predictions on the lateral and vertical extent of the groundwater SUHI in the study area.

     

Numerical modeling of contaminant transport in a stratified heterogeneous aquifer with dipping anisotropy

The combined influence of dip angle and adsorption heterogeneity on solute transport mechanisms in heterogeneous media can be understood by performing simulations of steady-state flow and transient transport in a heterogeneous aquifer with dipping anisotropy. Reactive and non-reactive contaminant transport in various types of heterogeneous aquifer is studied by simulations. The hydraulic conductivity (K) of the heterogeneous aquifer is generated by HYDRO_GEN with a Gaussian correlation spectrum. By considering the heterogeneity of the adsorption distribution coefficient (K _d), a perfect negative correlation between lnK and lnK _d is obtained by using the spherical grains model. The generated K and K _d are used as input to groundwater flow and transport models to investigate the effects of dipping sedimentary heterogeneity on contaminant plume evolution. Simulation results showed that the magnitude of the dip angle strongly controls the plume evolution in the studied anisotropic and heterogeneous aquifer. The retarded average pore-water velocity (v/R) of the adsorption model significantly controls the horizontal spreading of the plume. The bottom plume is intensively retarded in the zones between the dipping lenses of lower hydraulic conductivity and the no-flow bottom boundary. The implications of these findings are very important for the management of contaminated heterogeneous aquifers.     

基于人工神经网络的地下水状态方程—以山东莱州湾南岸地下水研究为例

山东莱州湾南岸地下水密度受多种因素影响,传统方法难以拟合,利用人工神经网络（简称ANN）高效的自组织学习能力和抗实验噪音能力及应用神经网络工具箱设计了三层反馈式神经网络模型,得出映射莱州湾南岸地下水密度与含盐量之间相关关系的基于ANN的地下水状态性关系的性质,基于ANN的地下水状态方程具有形式简单,易于使用等特点。     

Stream water bypass through a meander neck, laterally extending the hyporheic zone

A meander lobe neck diverts stream water into a hyporheic flow path adjacent to a low gradient stream, Little Kickapoo Creek, Illinois, USA. Hyporheic processes have been well-documented in surface water–groundwater mixing zones underlying and directly adjacent to streams. Alluvial aquifers underlying meander necks provide a further extension of the hyporheic zone. Hydraulic head and temperature data, collected from a set of wells across a meander neck, show stream water moves through the meander neck. The hydraulic gradient across the meander neck (0.006) is greater than the stream gradient (0.003) between the same points, driving the bypass. Rapid subsurface response to elevated stream stage shows a hydraulic connection between the stream and the alluvial aquifer. Temperature data and a Peclet number (Pe) of 43.1 indicate that thermal transport is dominated by advection from the upstream side to the downstream side of the meander neck. The temperature observed within the alluvial aquifer correlates with seasonal temperature variation. Together, the pressure and temperature data indicate that water moves across the meander neck. The inflow of stream water through the meander neck suggests that the meander system may host biogeochemical hyporheic zone processes.