汾河兰村—柴村段潜流带水动力特征研究.docx

本文格式为Word版，下载可任意编辑汾河兰村—柴村段潜流带水动力特征研究 季晨龙 李海明 肖瀚 李梦雨 王瑞芬 杨丽霞 张雪丽 摘 要：采用数值模拟方法研究汾河兰村—柴村段的潜流带水动力特征，量化潜水含水层对河流水位瞬时抬升1.5 m的响应幅度与影响距离，计算该条件下河流的侧渗补给强度结果说明：随着与河岸线垂直距离的增大，潜水响应幅度以及河流的侧渗补给流量、侧渗速度逐步减小，一致时刻起始点x=0处的侧渗速度为上游midstream>downstream. Within 5 to 90 d， the response amplitude and lateral seepage flow rate gradually decrease from the peak value to zero after 1 000 m vertical distance from the river riparian line. The response range of phreatic aquifer water level increases with time and the response range of phreatic water reaches 0.18 -0.25 m at 90 days at 500 m away from the river riparian line. The lateral seepage velocity and the flow rate reach the peak value within 5-30 d at the distance of 20-1 000 m from each section to the riparian line of the river and then decrease with time. The influence at the same time is the farthest from section 1-1′ and 2-2′ and 3-3′ are roughly equal. Key words： subsurface flow; hydrodynamic characteristics; lateral seepage; Lancun-Chaicun section of Fenhe River 隨着汾河流域经济社会进展速度加快，地下水位急剧下降，生态环境遭到严重破坏。

为解决汾河生态问题，山西省通过调引客水、污染防治、源头及干流两侧污染企业综合整治等方式举行大规模治理，流域地下水位得到持续上升[1]，但流域自身的产流功能仍未得到有效恢复[2]汾河流域早期曾拥有丰富的湿地资源，在教养水源、降解污染物、防洪蓄水等方面起到重要作用[3]，但后来大量滩涂、湿地被开发利用，湿地资源不断裁减[4]为恢复汾河流域湿地资源，山西省实施了汾河流域生态修复中游核心区干流蓄水工程，该工程主要通过河道蓄水、向堤外湿地补水，改善水生态环境该工程蓄水期间使流域湿地水位抬升1～2 m，在湿地补给地下水的过程中存在诸多问题有待深入研究 潜流带是指河流河床内水分饱和的沉积物层，是河流地表水与地下水相互作用的交汇区域，也是河床与河水举行物质和能量交换的区域[5]影响潜流交换的因素有好多种，其中河水流速与流量是主要因素[6-7]HUNT Bruce[8]研究说明，对于假设饱和模式下的河流与含水层系统，河水与地下水之间的水量交换强度同河水位与地下水位差成线性关系;PACKMAN Aaron I.等[9]通过水槽试验分析得出了潜流交换率与河水流速等因素之间的关系;HARVEY Judson W.等[10]得出了在河床邻近潜流交换随河流径流量变化的规律;LIERMANN Catherine Reidy[11]通过对Fauntleroy溪调查说明，不同流速对潜流交换势的影响不同。

对河床渗透性能时空变异性及其影响因素的识别不仅是研究河水与地下水水量交换的关键与难点[12-15]，也是河流与含水层相互作用研究的热点[16-17] 目前潜流交换的研究多是在小尺度上举行的，例如ZHANG Guotao等[18]研究了潜流带水交换与河岸线曲率的关系如何量化大尺度潜流交换已成为一个重要课题，侯宏冰等[19]对黄河温县段不同水位变化处境下的侧渗量举行了模拟计算，得出黄河水位的变化会引起地下水位的相应变化;XI Haiyang等[20]在黑河流域下游举行的田间试验说明，河流侧渗主要受流量及持续时间的操纵对河流与含水层系统的数值模拟是区域尺度上河水与地下水交换研究的重要手段[19，21-24]，其中对河流与含水层系统模型的概化以及对模型参数的估算是影响研究结果的关键[17，25] 鉴于目前对大尺度潜流交换的研究较少，笔者采用数值模拟的方法，从侧渗速度和侧渗补给流量两方面，研究位于汾河流域生态修复中游核心区干流蓄水工程受水区的上兰村—柴村桥段水位瞬时抬升1.5 m时的潜流带水动力特征，以期为汾河流域湿地生态修复供给参考 1 研究区水文地质条件 研究区位于太原市西张盆地，最高海拔为850 m，地势总体呈北高南低、东西两侧高中间低的特点，地质年头为第四系全新统及上更新统，地貌类型为冲洪积倾斜平原地貌、冲积平原地貌以及冲洪积交接凹地地貌。

研究区内汾河长度约10 km，属平原性河流研究区地下水主要接收汾河的渗漏补给以及大气降水入渗补给，地下水位低于河水位，河水同时向河两岸含水层侧向渗透补给地下水降水冬春少而夏秋多，降水量集中在7—9月，多年（1956—2022年）平均降水量为507 mm研究区潜水径流主要受地形、含水层岩性、补给来源等操纵，总体上由北向南、由东西两侧向中部滚动研究区潜水的排泄方式主要有蒸发、越流排泄及含水层的侧向径流3种形式研究区地下水为松散岩类孔隙水，含水层多为全新统砂卵砾石层，在0～160 m埋深范围内分为潜水含水层和弱透水层，包气带厚度为10～22 m，岩性多为粉土、亚砂土潜水含水层顶板高程620～824 m，厚度30～150 m，岩性多为砂卵石弱透水层厚度10～80 m，底板埋深40 ～160 m，岩性多为黏土 2 地下水流数值模型构建 2.1 水文地质概念模型的建立 2.1.1 边界条件 根据地形地貌、地层布局及区域流场特征，把潜水含水层南边界与东边界中下部概化为定流量边界，东边界上部、西边界以及北边界概化为隔水边界，弱透水层均概化为隔水边界（见图1） 2.1.2 含水层空间概化 本文只关注与河流水力联系紧密的潜水含水层，根据已有的地质图、剖面图、地质报告等资料，将含水层概化为非均质各向同性潜水含水层，地下水流以水平运动为主。

研究区含水层三维布局与剖面见图2 2.2 地下水流数值模型的建立 2.2.1 网格剖分 对模拟研究区举行网格剖分，共剖分为20 000个网格（100行×100列×2层），其中活动单元格8 368个，垂向上划分为两层（见图3） 2.2.2 模型识别应力期选取 根据现有资料，模型以2022年1月—2022年12月为校正期，2022年1月—2022年12月為验证期，以月为应力期，每个应力期为一个时间步长 2.2.3 源汇项处理 本文模型中潜水含水层地下水补给项主要包括降水入渗补给、地表水渗漏补给，排泄项主要包括蒸发、径流排泄 2.3 模型的校正与验证 将研究区校正期和验证期潜水流场水位模拟值与实测值举行对比，察觉两者的拟合程度较高，说明所建立的水文地质概念模型和数值模型根本合理，能较好地模拟研究区地下水动态变化过程以2022年1月地下水位为初始条件，模拟预料当汾河水位瞬时抬升1.5 m时河水侧渗补给地下水的速度和补给量以及地下水位的时空变化规律 3 结果与议论 3.1 不同断面地下水位变化特征 为了更好地熟悉河流补给地下水的水动力特征，沿汾河水流方向选取不同断面，研究其地下水位变化特征、河水补给地下水侧渗速度和补给量的变化特征。

在汾河上游（上兰村）、中游（汾河大桥）、下游（柴村桥）各设置一个断面，分别记为1—1′、2—2′、3—3′，以河流断面与河岸线的交错点为原点，以垂直河流中心线且远离河流的方向为x轴将各断面的潜水水位与初始水位的差值记为Δh （1）不同时刻Δh与河岸线距离的关系曲线见图4，由图4可知：Δh随与河岸线距离的增大总体呈下降趋势，在x=0处，Δh=1.5 m，随后Δh值随x增大在0～20 m急速下降，降幅占河流水位瞬时抬升高度的78.7%～97.3%，在20 m以后缓慢下降直至在距河岸线足够远处为0①1—1′上兰村断面河水位抬升5 d时Δh从河岸线处的1.5 m降为20 m处的0.04 m，降幅为1.46 m，占河流水位瞬时抬升高度的97.3%，之后从20 m到3 000 m的距离内，从0.04 m降至0 m;90 d时Δh在距河岸线20 m处从1.5 m降至0.32 m，降幅为1.18 m，占河流水位瞬时抬升高度的78.7%，之后从20 m到3 000 m的距离内，从0.32 m降至0.03 m，降幅为0.29 m②2—2′汾河大桥断面河水位抬升5 d时Δh在距河岸线20 m处降幅占河流水位瞬时抬升高度的96.7%，之后从20 m到3 000 m的距离内，降至0 m;90 d时Δh在距河岸线20 m处降幅占河流水位瞬时抬升高度的81.3%，之后从20 m到3 000 m的距离内降至0 m。

③3—3′柴村桥断面河水位抬升5 d时Δh在距河岸线20 m处降幅占河流水位瞬时抬升高度的96.7%，之后从20 m到3 000 m的距离内降至0 m;90 d时Δh在距河岸线20 m处降幅占河流水位瞬时抬升高度的80%，之后从20 m到3 000 m的距离内降至0 m （2）为了更直观地查看不同断面Δh随距离的变化规律，选取各断面x=500 m处绘制Δh随时间的变化曲线（见图5），察觉各断面Δh随时间推移均呈增大趋势，并且各断面在同一时刻的Δh值1—1′上兰村断面最大，2—2′汾河大桥断面与3—3′柴村桥断面大致接近，说明受河流水位抬升影响后的地下水位在上游涨幅较大1—1′上兰村断面距河岸线500 m处经过90 d后Δh增幅最大，为0.25 m2—2′汾河大桥断面与3—3′柴村桥断面Δh增幅大致接近，2—2′汾河大桥断面Δh在90 d内由0增加到0.18 m，3—3′柴村桥断面在90 d内由0增加到0.19 m故距河岸线500 m处90 d潜水响应幅度范围为0.18～0.25 m （3）定义Δh=0.01 m所对应的垂直于河岸线的距离为河流水位瞬时抬升1.5 m时对潜水水位的影响距离。

根据不同断面的影响距离（见表1）可得：一致时间下，影响距离1—1′上兰村断面最远，2—2′汾河大桥断面与3—3′柴村桥断面相近1—1′上兰村断面5、10、30 d影响距离分别为458、851、2 150 m，60 d时超出研究区范围;2—2′汾河大桥断面5、10、30、60、90 d影响距离分别为349、657、1 438、2 041、2 538 m;3—3′柴村桥断面5、10、30、60、90 d影响距离分别为364、642、1 360、2 030、2 517 m 3.2 侧渗速度的变化规律 根据Darcy定律，河流的侧渗速度v可用如下公式来表示[26]： 式中：K为渗透系数;Δh0，t为在距河岸线0 m的位置经过t时间水位上升的幅度;a为压力传导系数，a=Khm/μ（hm为含水层厚度，μ为给水度）;x为距离河岸线距离 渗透系数K与给水度μ根据数值模拟结果输入，含水层厚度hm。