含水层物理堵塞的回灌试验与数值模拟（Recharge test and numerical simulation of physical plugging of aquifer）.pdfVIP

谨以此论文献给我的母校——中国海洋大学，以及所 有支持我、鼓励我、关心我的师长和朋友们！ 单蓓蓓 含水层物理堵塞的回灌试验与数值模拟 摘 要 以雨洪水作为回灌水源，人工补给地下含水层是调节和控制地下水储量的有 效措施。然而，回灌过程中产生的堵塞问题（特别是物理堵塞），将严重影响回 灌效果，限制回灌工程的使用和推广。 本文借助室内模拟试验系统研究了定水头回灌条件下雨洪水中所含悬浮颗 粒对含水层堵塞的影响机理，着重研究了回灌液悬浮物浓度和水动力条件对回灌 堵塞发生及发展过程的影响。此外，本文就悬浮颗粒在介质中的运移-沉积过程 建立数学模型，以此来预测物理堵塞的发生和发展过程。该模型综合考虑了悬浮 颗粒在多孔介质中的吸附-脱附作用，并结合其对介质孔隙度的影响，得到介质 渗透系数的变化。同时，利用此数学模型，结合青岛市大沽河流域地下含水层的 实际回灌条件，对该区域地表漫灌和井灌注水过程中含水层的物理堵塞开展预测。 研究得出以下几点结论: (1) 在垂直入渗回灌条件下，含水层的物理堵塞现象均呈现“在回灌初期， 砂柱相对渗透系数 K’下降幅度较大；随着回灌的进行，K’下降幅度减小并趋于 稳定”的规律。 (2) 回灌液悬浮物浓度是影响含水介质堵塞程度的重要因素。悬浊液浓度增 大，则堵塞越严重。悬浮颗粒首先沉积在砂柱的表层，这在一定程度上抑制了悬 浮颗粒的向下迁移，所以回灌对深部砂层渗透性影响较弱。通过不同悬浮物浓度 回灌试验，当悬浮物浓度为25 、100、200mg/L 时，回灌100h 后，砂柱表层OB 段（0~11cm ）的相对渗透K’分别降为0.75、0.32 和0.12，整个砂柱K’分别变为 0.85、0.66 和0.56 。 (3) 回灌过程中的渗流速度影响堵塞程度。流速越大，流体对悬浮颗粒的剪 切力越大，颗粒越容易从孔隙壁上被“剥离”，悬浮颗粒运移深度增加，造成柱体 内部的渗透性减弱。在不同的流速条件下，表层相对渗透系数 K’变化趋势基本 相同，最终都降到0.10 左右；当水头差为40.00cm 和105.00cm 时，回灌结束后 CF 段（19~34cm）K’分别稳定在0.44 和0.54, 整个砂柱K’降为初始的61%和55% 。 (4) 悬浮颗粒在介质中的吸附-脱附是动态变化的，受渗流速度的影响。当流 速较小时，颗粒只发生吸附作用；随着流速的增大，颗粒开始逐渐从介质表面脱 I 落，此刻吸附和脱附现象并存。颗粒的吸附和脱附造成了颗粒的运移和沉积，使 得介质孔隙体积减小导致堵塞。可以通过实验测定颗粒吸附系数、脱附系数与渗 流速度之间的函数关系，以此来表征其动态变化过程。 (5) 基于悬浮颗粒的吸附-脱附现象，建立颗粒在介质中的运移-沉积模型， 可以对回灌过程中含水介质的物理堵塞现象进行定量描述，预测在不同回灌条件 下的堵塞程度。通过模型敏感度分析可以得知回灌液中悬浮颗粒的浓度是影响堵 塞程度的最重要因素。 (6) 结合研究区实际场地条件和回灌水源水质特征，利用物理堵塞模型，并 耦合地表漫灌和潜水井灌过程中的地下水渗流模型，可以预测含水层的物理堵塞 情况。预测结果表明无论是丰水期还是枯水期，物理堵塞主要发生在靠近回灌池 底的表层含水层（漫灌）以及注水井壁10cm 范围内含水层 （井灌）。 (7) 对回灌水源进行预处理（如通过过滤滤除部分悬浮颗粒，以降低回灌水 中悬浮物浓度）是预防含水层物理堵塞的主要措施。针对不同的回灌方式，可以 采取不同的堵塞治理手段，如对回灌池底进行周期性的清淤（漫灌），对回灌井 进行周期性的反冲刷 （井灌）。 关键词：人工回灌 多孔介质 悬浮颗粒 物理堵塞 数值模拟 II Recharge Experiment and Numerical Simulation on Aquifer Physical Clogging Abstract Artificial Recharge, whi