床面形态驱动下潜流交换试验

为揭示河床形态特征引起的潜流交换规律,构建循环式水槽装置,通过NaCl示踪对比分析了4种河床地形驱动下的潜流交换规律,并基于扩散理论探讨了潜流交换与地表水水动力及床沙渗透特性之间的关系.结果表明,潜流交换可以发生在平坦河床地形,且交换速率随地形起伏度和雷诺数的增大而增大;在相似的地形条件下,地表流速是影响潜流交换速率的主导因素,地表水深对潜流交换速率影响较弱;此外,分析表明,有效扩散系数与河床特征粒径之间具有2次方幂率函数关系,且与河床渗透系数成正比;交换深度d正比于雷诺数Re的1/2次方.     

河床沉积物非均质性影响下的潜流交换数值模拟

为了揭示河床沉积物非均质性对潜流交换的影响,构建了沙波地形作用下的地表水-地下水耦合模型,通过生成不同的渗透系数随机场,讨论了不同地表水动力过程和河床沉积物非均质性对潜流交换通量、交换空间及平均停留时间的影响规律。结果表明,所构建的地表水-地下水耦合模型能够准确刻画水沙界面附近流场,模型具有良好的适用性;均质或非均质河床沉积物情景下,水沙界面上的平均交换通量、停留时间与雷诺数之间均呈现幂函数关系,潜流交换深度则在地表水进入完全紊流之后趋于稳定。结果还表明,较强的河床沉积物非均质性能够有效增强水沙界面上潜流交换通量和空间交换频率,但会制约潜流交换空间并缩短水流在潜流带中的停留时间。     

河流中潜流交换研究进展

河水和地下水交换——潜流交换对溶质和污染物的归宿起着重要作用。潜流交换机理主要包括泵吸交换和冲淤交换。泵吸交换是由于河床形态引起的水头梯度,这些水头梯度诱导了对流传输;冲淤交换发生是由于移动河床截留和释放孔隙水。潜流交换的主要影响因素包括:河道流量、河床水力传导性、河床形态、河道弯曲、河床不均匀和背景条件。还探讨了反应性溶质、胶体颗粒共存情况下潜流交换的规律。对潜流交换研究现存的主要问题及未来研究展望提出了看法。     

分布式水文模型中河流水深及其与地下水相互作用的耦合模拟

大尺度流域水文模型一般只模拟河道径流,对河流水深和流速并不关注。在进行河流-地下水相互作用模拟时,河流水深（水位）是一个重要因素,其时空变化及其对河流-地下水水量交换的影响应加以考虑。本文就流域分布式水文模型中河流水深的时空变化计算及其与地下水的实时耦合模拟进行了研究,提出了相应的计算方法,改进了大尺度水文模型WATLAC,并通过V-型流域考题进行了检验。模拟结果显示,模型有效地模拟了V-型流域的河流水深、地下水水位沿河道的时空变化以及河流与地下水的水交换量,揭示了河流-地下水相互作用关系在降雨过程中的变化规律及主控因子。此外就河床糙率对河流水深及河流与地下水间水交换量的影响做了模拟分析,发现河床糙率的改变将影响河流水深,从而进一步影响河流与地下水的水交换量。本文提出的算法较为真实地模拟了河流洪水演进过程及其对河流-地下水相互作用的影响,模型适用于河流-地下水相互作用明显的区域,可作为评估地表-地下水相互影响的有效计算工具。     

河床起伏条件下生物堵塞对潜流带氮素迁移转化的影响

生物作用造成河床於堵普遍存在,并受到多种环境因素的影响.以往研究较少考虑生物堵塞对潜流带氮循环的影响.采用数值模拟方法建立水流-溶质运移-生物生长耦合模型,模拟河床驱动下的潜流交换过程,探讨生物堵塞作用对氮素迁移转化的影响.研究表明:生物作用导致沙丘表面形成堵塞层和生物聚集区,并限制了NO_3~-的运移深度,潜流带中以NH_4~-为主;河水中充足的有机碳会增大堵塞程度,促进反硝化作用去除NO_3~-;强烈的河床起伏和较大的河流流速均会增大沙丘表面的压力差,增大堵塞范围和细菌生物量,从而加快硝化和反硝化作用速率,有利于NO_3~-、NH_4~-的去除.生物堵塞作用有利于氮污染物的去除,且主要受到河床表面的压力梯度和异养细菌生物量的影响.     

胶体颗粒在潜流带中沉积过程及机制

潜流带是河流地表水和地下水交混区域,是河流中重要的物质能量交换和水生生物栖息的场所,而胶体颗粒在潜流带中沉积,会改变潜流带中的水动力结构和生态环境。本文利用室内循环水槽实验和多物理场耦合的数值模拟方法,旨在研究胶体颗粒在河流上覆水与潜流带中的迁移过程和胶体颗粒在潜流带中沉积分布特征及其对不同因素的响应规律。结果表明:河流上覆水中胶体会逐渐被河床截留且截留胶体集中于河床浅层;沙波水平方向截留量呈现出迎水面较高、背水面较低的趋势;胶体在潜流带沉积的主要机制是潜流交换、颗粒沉降与河床截留作用。本文能为胶体颗粒在潜流带中的生态环境作用研究提供科学依据,并为河流生态环境修复、河流健康管理提供理论支持。     

地表水与地下水相互作用的温度示踪与模拟研究进展

刻画地表水与地下水的相互作用过程及精确计算二者的交换量始终是个挑战,而新兴的温度示踪方法对于这方面的研究具有独特优势。重点介绍了地表水与地下水相互作用的温度示踪方法原理、应用及相关模拟的研究进展。温度示踪方法的数据获取成本低且温度适宜密集与连续监测,可对地表水与地下水相互作用过程进行精细刻画;在进行地表水与地下水相互作用的水-热耦合模拟时,温度数据可用于进一步校正模型,降低模型的不确定性,以提高交换量的计算精度。监测河水温度和河床沉积物内的水温,运用温度示踪方法研究河水与地下水的相互作用过程,并将水流及热运移数值模型相耦合,利用多源数据进行模型校正,可精细刻画地表水与地下水的交换量、相互作用带内的水流途径、流速及其变化趋势和控制因素,为地表水与地下水相互作用研究提供新的模型与方法。     

河岸带潜流层温度示踪流速计算方法

为揭示低温水影响下的河岸带潜流层的温度场和流场分布特性,利用野外水温水位实时监测试验,研究河岸带潜流层温度场在不同季节、不同空间位置上的分布特性,并利用水温资料计算获得地下水流速。结果表明:河岸带潜流层温度场在夏季和冬季分别呈现出"上暖下冷"和 "上冷下暖"的温度分层现象;通过对温度示踪方法的4种计算方法进行分析比较,得到Hatch相位法计算的地下水流速具有较高的准确性,在2014年12月15—31日时段内流速大小为1.03×10-4~7.96×10-4m/s,在空间上,断面深度增加,地下水流速降低,且不同深度流速曲线接近平行。     

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

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

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

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

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.     

A new method for mapping variability in vertical seepage flux in streambeds

A two-step approach was used to measure the flux across the water-sediment interface in river channels. A hollow tube was pressed into the streambed and an in situ sediment column of the streambed was created inside the tube. The hydraulic gradient between the two ends of the sediment column was measured. The vertical hydraulic conductivity of the sediment column was determined using a falling-head permeameter test in the river. Given the availability of the hydraulic gradient and vertical hydraulic conductivity of the streambed, Darcy’s law was used to calculate the specific discharge. This approach was applied to the Elkhorn River and one tributary in northeastern Nebraska, USA. The results suggest that the magnitude of the vertical flux varied greatly within a short distance. Furthermore, the flux can change direction from downward to upward between two locations only several meters apart. This spatial pattern of variation probably represents the inflow and outflow within the hyporheic zone, not the regional ambient flow systems. In this study, a thermal infrared camera was also used to detect the discharge locations of groundwater in the streambed. After the hydraulic gradient and the vertical hydraulic conductivity were estimated from the groundwater spring, the discharge rate was calculated.     

Determination of the anisotropy of an upper streambed layer in east-central Nebraska, USA

Information on the anisotropy of streambed hydraulic conductivity (K) is a necessity for analyses of water exchange and solute transport in the hyporheic zone. An approach is proposed for the determination of K, developed from existing in-situ permeameter test methods. The approach is based on determination of vertical and horizontal hydraulic conductivity of streambed sediments on-site and eliminates the effects of vertical flow in the hyporheic zone and stream-stage fluctuation, which normally influence in situ permeameter tests. The approach was applied to seven study sites on four tributaries of the Platte River in east-central Nebraska, USA. On-site permeameter tests conducted on about 172 streambed cores for the determination of vertical hydraulic conductivity (K _v) and horizontal hydraulic conductivity (K _h) at the study sites indicate that the study sites have relatively small anisotropic ratios, ranging from 0.74 to 2.40. The ratios of K _h to K _v from individual locations within a study site show greater variation than the anisotropic ratios from the mean or median K at each of the study sites.     

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

为研究城市小型河流中污染物的物理迁移过程规律,分析基流条件下流动水体与暂态存储区之间的滞留交互作用,采用溴化锂(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个河段对比发现,城市河流渠道化、河床沉积物贫瘠等特征导致潜流交换能力弱化.     

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.     

Heterogeneity of hydraulic conductivity and Darcian flux in the submerged streambed and adjacent exposed stream bank of the Beiluo River,northwest China

Recognizing the heterogeneity of hydraulic conductivity and hyporheic flow is critical for understanding contaminant transfer and biogeochemical and hydrological processes involving streams and aquifers. In this study, the heterogeneity of hydraulic conductivity and Darcian flux in a submerged streambed and its adjacent exposed stream banks were investigated in the Beiluo River, northwest China. In the submerged streambed, Darcian flux was estimated by measurement of vertical hydraulic conductivity (K _v) and vertical head gradient (VHG) using in-situ permeameter tests. On exposed stream banks, both horizontal hydraulic conductivity (K _h) and K _v were measured by on-site permeameter tests. In the submerged streambed, K _v values gradually decreased with depth and the higher values were concentrated in the center and close to the erosional bank. Compared to the exposed stream banks, the K _v values were higher in the streambed. From stream stage to the topmost layer of tested sediment, through increasing elevation, the K _h values increased on the erosional bank, while they decreased on the depositional bank. The values of VHG along the thalweg illustrate that downwelling flux occurred in the deepest area while upwelling flux appeared in the other areas, which might result from the change of streambed elevation. The higher value of the Darcian flux in the submerged streambed existed near the erosional bank.     

Influence of streambed hydraulic conductivity on solute exchange with the hyporheic zone

A conservative solute tracer experiment was conducted in Indian Creek, a small urban stream in Philadelphia, Pennsylvania to investigate the role of subsurface properties on the exchange between streamwater and the hyporheic zone (subsurface surrounding the stream). Sodium Bromide (NaBr) was used as a conservative tracer, and it was monitored in the surface water at two stations and in the upper bed sediments (shallow hyporheic zone extending from 7.5 to 10 cm below the streambed). The hydraulic conductivity (K) of the upper bed sediments and the lower bed sediments (10–12.5 cm below the streambed) was measured in situ. High tracer concentrations were observed in the upper layer at locations where the hydraulic conductivity of the upper layer was larger than that of the lower layer. Low concentrations in the upper layer were observed in the converse case. A statistically significant relationship between the mass retained in the upper layer and the difference of K values between layers was observed.