高寒山区河水与地下水相互作用的温度示踪:以黑河上游葫芦沟流域为例

高寒山区的地表水与地下水相互作用的定量研究对水资源的评价及管理等具有重要意义,而目前在高寒山区开展的地表水与地下水相互作用的定量研究相对较少.以黑河上游葫芦沟流域为研究区域,采用温度示踪方法对高寒山区河水与地下水的相互作用进行了研究,并对温度示踪方法在高寒山区的适用性进行了讨论.监测了研究区两个时段的地温、河水水位、地下水水位以及河床沉积物底部不同深度处的温度,并对温度系列数据进行定量分析,计算了不同位置处河水入渗流速.结果表明:研究区河水水位普遍高于地下水水位;河床底部温度在9月份整体低于7月;流速计算结果表明监测时段内主要为河水入渗补给地下水,入渗速率整体介于2×10-6~5×10-5 m/s.温度示踪法在高寒山区的适用性分析表明:在地下水受多途径补给时,温度示踪法仅指示河水对地下水的补给,而其他水源对地下水的补给还要通过同位素方法和数值模拟等其他手段进行计算.影响高寒山区河水对地下水补给的因素主要有:河水与地下水水位、河床沉积物的水力传导系数与热容.      

利用稳定同位素识别黑河流域地上水的补给来源

本文利用稳定同位素(2H和18O)及水化学方法识别黑河流域地下水的补给来源,估算黑河水与地下水的转化数量。研究结果表明,黑河流域地下水的主要补给来自山区出山河流,山前戈壁带是地下水快速补给区,中下游盆地地下水补给来源为引河灌溉和河流侧渗。黑河干流出山河水在张掖以上河段约4.4×108m3/a渗漏补给地下水,约占出山迳流量的27%。张掖—正义峡河段道地下水向河道平均排泄量为11.4×108m3/a,占该段河流迳流量的69%。研究成果不仅对黑河流域地下水的开发管理有着重要意义,对我国西北类似的内陆盆地地下水的开发管理有着借鉴意义。     

黑河干流浅层地下水与地表水相互转化的水化学特征

通过分析黑河干流地表水与地下水的水化学特征,识别沿黑河干流不同地带地下水与地表水的相互转化关系.研究结果表明:(1)在祁连山区,地下水与地表水的转化以地下水向河流排泄为主.(2)南部盆地,在山前戈壁带,出山河水入渗转化为地下水;溢出带地下水以泉的形式转化为地表水;进入细土平原后,汛期河水补给地下水,非汛期地下水补给河水;在农灌区引河水通过田间入渗补给地下水.(3)北部盆地,在金塔灌区,地下水主要接受引水灌溉入渗补给;在金塔灌区到额济纳旗,河流入渗转化为地下水.     

基于分布式光纤温度示踪探测裂隙岩体地下水渗流特征

认识深部裂隙岩体中的地下水渗流特征(流速、渗流路径等),是深部地质工程开发建设的重要前提。近年来,分布式光纤测温技术作为识别深部裂隙岩体地下水渗流特征的有效方法,在国外开展了大量的研究,但在国内鲜少见在实际场地开展的相关工作。本研究以我国首个地下实验室场址甘肃北山新场花岗岩岩体中的两个钻孔(BSQ02及BSQ03)为试验对象,开展基于分布式光纤测温(Fiber-Optic Distributed Temperature Sensing,FO-DTS)的现场温度-水力试验,实现了对钻孔地下水温度的高精度、连续性观测。通过分析现场试验获取的钻孔温度-深度剖面随时间的变化,推断BSQ02在试验过程中存在外源地下水的流入,然后结合钻孔柱状图对钻孔中的入流导水裂隙进行了定位;基于现场观测数据建立了钻孔的渗流-传热耦合数值模型,反演估算出钻孔中地下水平均流速为0.01 m·s-1,通过裂隙流入地下水温度小于钻孔中原地下水温度,两者之间的温度差为0.7 ℃,通过裂隙流入的地下水流速为1×10-5 m·s-1,获取了地下水的渗流特征。该项工作可为基于分布式光纤测温技术的裂隙介质地下水渗流规律研究提供借鉴与指导。     

格尔木河222Rn同位素变化及其对地表水-地下水交互关系的指示意义

我国西北内陆干旱区水资源匮乏,生态环境脆弱,在全球气候变化和人类活动干扰背景下,采用同位素方法进行精细尺度地表水-地下水交互作用研究是探求当地水循环变化和水资源管理的基本要求。通过测量格尔木河流域河水、地下水样品2019年5月和8月的²²²Rn浓度和典型断面流量,结果发现：山区河段河水²²²Rn浓度最高,平均值为948.72 Bq·m^-3,指示基岩裂隙水是山区河段重要补给来源;山前冲洪积扇河水²²²Rn浓度最低,平均值为76.71 Bq·m^-3,地下水补给较少;溢出带地区河水²²²Rn浓度上升至平均676 Bq·m^-3,地下水溢出补给河水,向下至细土平原,河水²²²Rn浓度呈下降趋势。时间变化上,8月与5月相比,河水²²²Rn浓度下降,表明地下水补给减少。溢出带S1~S2河段河水与地下水交互关系以双向转化为主,基于质量守恒原理计算河水与地下水交互通量,5月和8月累积河水渗漏通量分别为3.87 m³?s^-1和0.9 m³?s^-1,地下水补给通量分别为0.51 m³?s^-1和0.47 m³?s^-1,河水渗漏强度大于地下水补给,二者交互通量存在时空差异。     

利用222Rn、EC及断面测流分析黑河中游地表水与地下水转换关系

本文以黑河中游为研究对象,采用222Rn、EC及断面测流法对莺落峡至正义峡河段地表水、地下水转化关系进行了分析。结果表明,莺落峡至张掖黑河大桥段为地表水补给地下水,其单长渗漏率为0.248～0.433m3/s.km;张掖黑河大桥至正义峡段为地下水补给地表水,其单长渗出率为0.042～1.057m3/s.km。在相互补给率计算基础上对不同河段不同时段可分别计算相互的补给量,为水资源评价及规划提供科学依据。     

黑河流域中游盆地地表水与地下水转化机制研究

地表水与地下水相互转化是中国西北干旱内流盆地水循环的显著特征,转化机制研究是盆地水循环规律认知和水资源可持续管理的重要基础。以我国西北干旱内流河黑河流域中游的张掖盆地和盐池盆地为研究区,建立了黑河主干河道时变水平衡模型和地表水地下水耦合数值模型,研究了长周期水文变化和人类活动双重影响下地表水与地下水转化机制,得到如下认识:（1）补给条件由以天然条件下河流渗漏为主的线状补给演变为以河流与引水渠道渗漏的线状补给和灌区田间入渗面状补给,排泄条件由以泉水溢出和天然湿地排泄演变为以泉水溢出与地下水开采为主的排泄。（2）张掖盆地黑河干流河道入渗段和溢出段大致以G312 大桥为界,亦称为地表水与地下水转化的转折点。莺落峡—G312 大桥段为悬河渗漏段,河道入渗补给主要受控于进入河道的实际过水量。其中,莺落峡—草滩庄段河道入渗补给率为28.20 %;草滩庄—G312 大桥段河道入渗补给量与河道过水量的关系可用分段函数表达,河道过水量大于或等于0.37×108 m3/mon时呈幂函数关系,小于则呈线性函数关系。G312 大桥—正义峡段为地下水溢出段,其中G312大桥—平川大桥段地下水溢出量约占全部溢出量的70%,溢出峰值出现在高崖水文站下游约6 km处,其单长溢出量可达0.46 m3/(s·km)。（3）研究区是一个相对完整的河流—含水层系统,近31年来经历了连枯和连丰的水文变化,地下水补给排泄条件及与地表水转化机制均发生了相应的变化。地表水与地下水转化最强烈的地区为张掖盆地中部的黑河—梨园河倾斜平原。1990—2001 年连枯期,灌区引水量总体逐年减少,以河道入渗和渠系渗漏为主的补给量平均以0.06×108 m3/a速率减少,农田灌溉面积增加导致灌溉用水增加,地下水开采量显著增加,地下水水位逐年下降,储存量累计减少5.77×108 m3,地下水溢出量平均减少0.16×108 m3/a;而2002—2020 年连丰期,灌区引水量总体逐年减少,河道入渗量呈增加趋势,地下水总补给量平均增加0.15×108 m3/a,灌溉面积继续扩大,农灌开采量随之增加,以河道入渗量增加为主导,地下水水位持续上升,储存量累计增加5.45×108 m3,地下水溢出量平均增加0.08×108 m3/a。总之,补给和排泄条件变化较大,地下水储存量先减后增,地下水溢出总量变化较为平缓,反映了该区巨厚含水层系统的巨大调蓄功能。（4）位于张掖盆地东部的诸河倾斜平原地下水水位长期处于持续下降状态,这是由于地表水开发过度,补给量锐减。黑河侵蚀堆积平原地下水水位基本稳定。30 多年来盐池盆地倾斜平原地下水水位长期处于持续下降状态,这是由于移民开垦导致地下水过量开采。（5）内流盆地天然悬河入渗段是珍贵的地下水补给通道,无论连枯期还是连丰期,河道实际过水量是河道渗漏补给量的关键,保护上游天然河道和一定的河道实际过水量是内流盆地水资源可持续管理的关键。     

河道潜流交换水动力过程现场监测研究进展

河道潜流带是河水与地下水交互作用的关键过渡区,具有高度的动态性和敏感性,频繁的水力交换易引发其内部的生物地球化学反应,并对河流生态安全产生重要影响。为提高河道潜流交换水动力监测结果的真实性和有效性,介绍渗流、水压、溶质、温度等监测方法及其适用性,归纳河床和河岸的监测点布设参数和不同监测手段的数据采集频率,系统总结河段、河床及河岸的潜流交换理论与计算方法,探讨导致潜流监测结果不确定的因素。针对现有监测方法的不足和未来研究的发展方向,提出今后需加强多方法联合监测和多技术集成的应用,开展河岸侧向潜流监测及不确定性分析等工作,为潜流带水文-生物地球化学过程耦合机制研究提供可靠数据。     

祁连山区黑河上游多年冻土分布考察

高山多年冻土的分布及土壤季节冻融过程对地表水文过程、 生态系统、 碳循环及寒区工程建设等都有很大的影响. 黑河上游地处祁连山中东部, 属于高原亚寒带半干旱气候, 研究黑河流域多年冻土分布对于系统理解该区域的生态-水文过程、 气候与环境变化以及水资源评价、 工程建设等非常重要. 2011年6-8月对黑河干流源头西支开展了多年冻土调查, 沿二尕公路(S204)在热水大坂垭口至石棉矿岔口之间区域, 完成测温孔7眼, 并布设测温管进行地温监测. 根据勘察、 钻探及测温资料, 确定了黑河源头地区山地多年冻土下界为海拔3 650~3 700 m之间. 受高度地带性的控制, 随着海拔的降低, 活动层厚度由在海拔4 132 m时的1.6 m增加至在多年冻土下界处的约4.0 m, 多年冻土年平均地温也相应的由-1.7℃增加到0.0℃左右, 而多年冻土厚度由100 m以上减小到多年冻土下界处的0.0 m. 同时, 坡度和坡向、 岩性、 含水(冰)量、 地下水、 河水等局地因素对多年冻土温度和厚度也有重要的作用.     

黑河中上游段河道渗漏量计算方法的试验研究

通过河道断面法对黑河中上游河段河道渗漏量进行了比测试验,结果表明:河床岩性组成、河流平面几何形态变化以及人类活动对输水方式的改变等对河道渗漏量影响显著.黑河实施全流域调水前后河道渗漏量发生显著的跳跃性变化,莺落峡至312桥区间河道渗漏量调水前为1.2543×108 m3,调水后为4.3747×108 m3,调水后较调水前河道渗漏量增加3.1203×108 m3;河道流量越小其相应损失率愈大,河道流量越大,其相应损失率愈小,由此得出各河段不同流量单位河长的损失率服从倒数函数分布.研究结果是在黑河中上游实体河道上进行上下断面流量比测试验的基础上得出的,符合该河段的实际情况,可作为黑河中游地表水、地下水转化分析、地下水流系统模拟、流域水资源重复利用、水量调度计算各种模拟模型验证河道渗漏量的基础资料.     

Heat tracing to examine seasonal groundwater flow beneath a low-gradient stream in rural central Illinois,USA

The thermal profile of a streambed is affected by a number of factors including: temperatures of stream water and groundwater, hydraulic conductivity, thermal conductivity, heat capacity of the streambed, and the geometry of hyporheic flow paths. Changes in these parameters over time cause changes in thermal profiles. In this study, temperature data were collected at depths of 30, 60, 90 and 150 cm at six streambed wells 5 m apart along the thalweg of Little Kickapoo Creek, in rural central Illinois, USA. This is a third-order low-gradient baseflow-fed stream. A positive temperature gradient with inflection at 90-cm depth was observed during the summer period. A negative temperature gradient with inflection at 30 cm was observed during the winter period, which suggests greater influence of stream-water temperatures in the substrate during the summer. Thermal models of the streambed were built using VS2DHI to simulate the thermal profiles observed in the field. Comparison of the parameters along with analysis of temperature envelopes and Peclet numbers suggested greater upwelling and stability in temperatures during the winter than during the summer. Upwelling was more pronounced in the downstream reach of the pool in the riffle and pool sequence.     

Investigating spatial variability of vertical water fluxes through the streambed in distinctive stream morphologies using temperature and head data

Investigating the interaction of groundwater and surface water is key to understanding the hyporheic processes. The vertical water fluxes through a streambed were determined using Darcian flux calculations and vertical sediment temperature profiles to assess the pattern and magnitude of groundwater/surface-water interaction in Beiluo River, China. Field measurements were taken in January 2015 at three different stream morphologies including a meander bend, an anabranching channel and a straight stream channel. Despite the differences of flux direction and magnitude, flux directions based on vertical temperature profiles are in good agreement with results from Darcian flux calculations at the anabranching channel, and the Kruskal-Wallis tests show no significant differences between the estimated upward fluxes based on the two methods at each site. Also, the upward fluxes based on the two methods show similar spatial distributions on the streambed, indicating (1) that higher water fluxes at the meander bend occur from the center of the channel towards the erosional bank, (2) that water fluxes at the anabranching channel are higher near the erosional bank and in the center of the channel, and (3) that in the straight channel, higher water fluxes appear from the center of the channel towards the depositional bank. It is noted that higher fluxes generally occur at certain locations with higher streambed vertical hydraulic conductivity (K _v) or where a higher vertical hydraulic gradient is observed. Moreover, differences of grain size, induced by stream morphology and contrasting erosional and depositional conditions, have significant effects on streambed K _v and water fluxes.     

Heat tracing to determine spatial patterns of hyporheic exchange across a river transect

Significant spatial variability of water fluxes may exist at the water-sediment interface in river channels and has great influence on a variety of water issues. Understanding the complicated flow systems controlling the flux exchanges along an entire river is often limited due to averaging of parameters or the small number of discrete point measurements usually used. This study investigated the spatial pattern of the hyporheic flux exchange across a river transect in China, using the heat tracing approach. This was done with measurements of temperature at high spatial resolution during a 64-h monitoring period and using the data to identify the spatial pattern of the hyporheic exchange flux with the aid of a one-dimensional conduction-advection-dispersion model (VFLUX). The threshold of neutral exchange was considered as 126 L m^?2 d^?1 in this study and the heat tracing results showed that the change patterns of vertical hyporheic flux varied with buried depth along the river transect; however, the hyporheic flux was not simply controlled by the streambed hydraulic conductivity and water depth in the river transect. Also, lateral flow dominated the hyporheic process within the shallow high-permeability streambed, while the vertical flow was dominant in the deep low-permeability streambed. The spatial pattern of hyporheic exchange across the river transect was naturally controlled by the heterogeneity of the streambed and the bedform of the stream cross-section. Consequently, a two-dimensional conceptual illustration of the hyporheic process across the river transect is proposed, which could be applicable to river transects of similar conditions.     

Statistical and geostatistical features of streambed hydraulic conductivities in the Platte River, Nebraska

This paper presents streambed hydraulic conductivities of the Platte River from south-central to eastern Nebraska. The hydraulic conductivities were determined from river channels using permeameter tests. The vertical hydraulic conductivities (K _v ) from seven test sites along this river in south-central Nebraska belong to one statistical population. Its mean value is 40.2 m/d. However, the vertical hydraulic conductivities along four transects of the Ashland test site in eastern Nebraska have lower mean values, are statistically different from the K _v values in south-central Nebraska, and belong to two different populations with mean values of 20.7 and 9.1 m/d, respectively. Finer sediments carried from the Loup River and Elkhorn River watersheds to the eastern reach of the Platte River lowers the vertical hydraulic conductivity of the streambed. Correlation coefficients between water depth and K _v values along a test transect indicates a positive correlation – a larger K _v usually occurs in the part of channel with deeper water. Experimental variograms derived from the vertical hydraulic conductivities for several transects across the channels of the Platte River show periodicity of spatial correlation, which likely result from periodic variation of water depth across the channels. The sandy to gravelly streambed contains very local silt and clay layers; spatially continuous low-permeability streambed was not observed in the river channels. The horizontal hydraulic conductivities were larger than the vertical hydraulic conductivities for the same test locations.     

A multi-method approach to quantify groundwater/surface water-interactions in the semi-arid Hailiutu River basin,northwest China

Identification and quantification of groundwater and surface-water interactions provide important scientific insights for managing groundwater and surface-water conjunctively. This is especially relevant in semi-arid areas where groundwater is often the main source to feed river discharge and to maintain groundwater dependent ecosystems. Multiple field measurements were taken in the semi-arid Bulang sub-catchment, part of the Hailiutu River basin in northwest China, to identify and quantify groundwater and surface-water interactions. Measurements of groundwater levels and stream stages for a 1-year investigation period indicate continuous groundwater discharge to the river. Temperature measurements of stream water, streambed deposits at different depths, and groundwater confirm the upward flow of groundwater to the stream during all seasons. Results of a tracer-based hydrograph separation exercise reveal that, even during heavy rainfall events, groundwater contributes much more to the increased stream discharge than direct surface runoff. Spatially distributed groundwater seepage along the stream was estimated using mass balance equations with electrical conductivity measurements during a constant salt injection experiment. Calculated groundwater seepage rates showed surprisingly large spatial variations for a relatively homogeneous sandy aquifer.     

Analytical methods that use natural heat as a tracer to quantify surface water–groundwater exchange, evaluated using field temperature records

Two methods applying natural heat as a tracer to quantify surface water–groundwater exchange were evaluated using field data. Arrays capable of monitoring and recording the streambed response to diurnal temperature variations in the surface water were deployed for a 2-month period in three locations in perennial pools at Maules Creek, New South Wales, Australia. Multi-level array design, field deployment and parameter estimation are discussed. The applicability of analytical solutions derived from the heat transport equation to the streambed environments was analysed using the recorded temperature time series. The stream was found to lose water to the aquifer, which was supported by simultaneously recorded hydraulic gradients. However, the one-dimensional (1D) analytical solutions did not adequately describe the observed streambed thermal response at two locations. The resulting artefacts in the estimated flow velocities are discussed. It was hypothesised that the artefacts originate from model limitation due to streambed heterogeneity and application of 1D solutions to multi-dimensional and dynamic streambed flow. This consequently imposes limitations on the field applicability of the methods. Nevertheless, in combination with time series of surface water and streambed water levels, the use of heat as a tracer provided a powerful tool for better understanding the shallow hydrogeological system.     

Surface water–groundwater interaction and chemistry in a mineral-armored hydrothermal outflow channel, Yellowstone National Park, USA

Small quantities of groundwater interact with hydrothermal surface water to drive in-stream geochemical processes in a silica-armored hot-spring outflow channel in Yellowstone National Park, USA. The objective of this study was to characterize the hydrology and geochemistry of this unique system in order to (1) learn more about the Yellowstone Plateau’s subsurface water mixing between meteoric and hydrothermal waters and (2) learn more about the chemical and physical processes that lead to accumulation of streambed cements, i.e., streambed armor. A combination of hydrological, geochemical, mineralogical, microscopic, and petrographic techniques were used to identify groundwater and surface-water exchange. Interaction could be identified in winter because of differences in surface water and groundwater composition but interaction at other times of the year cannot be ruled out. Dissolved constituents originating from groundwater (e.g., Fe(II) and Mg) were traced downstream until oxidation and/or subsequent precipitation with silica removed them, particularly where high affinity substrates like cyanobacterial surfaces were present. Because the stream lies in a relatively flat drainage basin and is fed mainly by a seasonally relatively stable hot spring, this system allowed study of the chemical processes along a stream without the obscuring effects of sedimentation.     

Methods to identify and reduce potential surface stream water losses into abandoned underground mines

Methods to identify and subsequently seal surface water loss zones in stream channels were tested by the United States Bureau of Mines at Staub Run, a first-order stream near Frostburg, Maryland, that partially overlies abandoned coal mine workings. Conventional stream gauging was conducted to establish discharge patterns before and after stream sealing. Electromagnetic terrain conductivity surveys were performed within the stream channel to identify zones of increased relative water saturation to depths less than 15 m. Zones of increased conductivity were generally found to be associated with areas exhibiting statistically significant (P 0.05) gauged flow losses. Conversely, zones that exhibited declining conductivity delineated areas where between-station flows were not significantly different. Using this information on potential loss zones, an experimental grouting procedure was applied by injecting an expandable polyurethane grout to a depth less than one meter into the alluvial streambed over a 180-m section of the stream channel. Before grouting, the study section exhibited a 24 I/sec flow loss; first-phase grouting reduced this to a 14 I/sec flow loss; with a second-phase grouting the losses were only 3 I/sec.     

城市小型河流水动力弥散和潜流交换过程

为研究城市小型河流中污染物的物理迁移过程规律,分析基流条件下流动水体与暂态存储区之间的滞留交互作用,采用溴化锂(LiBr)作为保守性示踪剂进行野外现场示踪试验,结合一维溶质运移存储模型(One-dimensional Transport with Inflow and Storage model, OTIS)定量解析潜流交换特性,估算纵向弥散系数(D)、潜流交换面积(A_s)、主河道断面面积(A)和潜流交换系数(α).模型度量指标D_aI值和均方根误差值结果表征参数模拟结果可靠性高,拟合效果理想.由泵入点O至下游1 300 m设置的A、B、C、D 4处监测点的模拟结果表明,水文参数D、A_s、A和α均随水文条件而变,OB河段(0~600 m)潜流交换能力较弱,主要以对流弥散过程为主;BD河段(600~1 300 m)具有较强的暂态存储能力,对溶质的滞留时间长;BC(600~1 000 m)和CD(1 000~1 300 m)河段交换系数分别为(3.42×10-6±0.65×10-6)s-1和(2.87×10-6±0.81×10-6 )s-1;河段BC存在2.2×10-5m3/(s·m)的侧向补给流量.4个河段对比发现,城市河流渠道化、河床沉积物贫瘠等特征导致潜流交换能力弱化.     

河床水力传导度及其各向异性的测定

采用直接测定法观测了黑河中游段河床水力传导度及其各向异性,结果表明:河床水力传导度不仅存在较强的各向异性,而且存在空间尺度上的变异性.河床中心位置在垂直、水平和θ=30°方向的平均水力传导度分别为0.45、22.49和1.71 m/d,河床边分别为5.95、29.69、16.80 m/d.在同一测点,水力传导度随着与水平方向的夹角增大呈幂函数曲线下降.试验结果表明河床边是河水的主要渗漏区,并且以侧渗为主.