包气带变化及其对地下水补给的影响

【摘要】地下水是华北平原的主要供水水源,主要接受降雨入渗补给。近几十年来,由于工农业及居民用水的双重作用,长期超采地下水(尤其是山前地带)导致浅层含水层水位埋深增大,出现大面积降落漏斗,引起包气带厚度增大,进而使得垂向入渗补给过程变得更为复杂。华北平原是我国重要的棉粮基地,农业用水占总用水量的70%以上。很多地区农业用水浪费严重,灌溉水量超过作物需求,导致灌溉回水成为该地区垂向入渗补给的重要组成部分。地下水补给研究在国内外干旱-半干旱地区(包括华北平原)已开展过很多,但包气带变化条件下入渗补给的复杂过程及补给评价的不确定性依然是我们面临的难题或挑战。准确合理地评价华北平原垂向入渗补给,对制定地下水合理开发利用方案和合理的灌溉制度、寻找灌溉-开采的平衡点、防止地下水位继续下降、保护该地区的地下水资源具有重要意义；同时,获得的补给速率可以为今后华北平原区域地下水模拟提供输入数据。为此,本文以包气带水文学、溶质运移理论等为支撑,在搜集研究区前人研究数据及研究结果的基础上,综合采用野外调查、室内外实验、分析测试等方法研究因长期超采引起的包气带变化及其对补给的影响,主要包括以下三部分内容：(1)通过一组长达7个月的非饱和-饱和带场地灌溉示踪实验,利用Tecplot软件将监测数据进行三维可视化,定性分析野外包气带中水流和溶质运移规律；同时根据时间矩分析法定量分析水分及溶质在三维非均质包气带中的入渗及运移规律。(2)选取因长期超采地下水的华北山前平原(石家庄地区)和中部平原(衡水地区)作为研究区开展地下水垂向入渗补给研究,综合利用质量平衡法、人工溴示踪法、非饱和达西法和一维非饱和流数值模型评价研究区垂向入渗补给,同时收集前人的研究结果作对比分析,并将各种方法评价的结果进行相互验证,以提高补给评价的可信度。(3)综合以上补给评价结果,讨论包气带变化(厚度增大、结构变化)和外界因素变化(土地利用、极端气候)对垂向补给的影响。通过上述研究得到以下主要结论：(1)非饱和-饱和带非均质场地灌溉示踪实验在低渗透性透镜体存在的非均质介质中,相对稳定的氯离子会产生明显的绕流现象；而氚的绕流现象不明显,即水流基本保持活塞流下移；但当示踪剂(氯离子)连续输入时间较短时,示踪剂向下运移规律并不能真实反映水流的真实运移特征。分析各控制点的穿透曲线可以得出,由于氯离子的轻微吸附作用使得氯离子峰值到达时间要比氚晚；上部低渗透性透镜体的存在对下部控制点的穿透时间影响不大,即上部低渗透性透镜体的存在对潜在补给的评价影响较小,故采用人工示踪剂评价潜在补给比较可信。分析淋洗曲线知若场地存在较多低渗透性透镜体,会使得氯离子的回收率偏低,导致利用氯离子质量平衡法计算的潜在补给速率比实际情况偏低。从整体上看,若非均质场地中高渗透性孔隙介质连通性较强,水分及溶质基本以活塞流形式下移。(2)多方法对比评价地下水垂向补给采用离子质量平衡法评价石家庄地区的垂向入渗补给得：F-、Cl和S042-适用于不灌溉林地的入渗补给评价,其补给速率依次为16.9、18.8和19.4mm/a。当考虑灌溉水输入,C1-适用于灌溉农田,评价得出灌溉农田区的补给速率为87.8-155.9mm/a。通过不同土地利用类型地区开展人工溴示踪实验得到：石家庄地区冬小麦-夏玉米农田的补给速率为54.5-466.4mm/a(其中2个点位的计算时段为雨季导致补给速率偏大)、果园的补给速率为110.1-564.3mm/a、树林的补给速率为0-20.5mm/a、平均补给系数为0.173；衡水地区冬小麦-夏玉米农田的补给速率为28.9-164.8mm/a、果园的补给速率为11.9-114.6mm/a、棉花地的补给速率为87.3-186.3mm/a、树林的补给速率为0-32.2mm/a,平均补给系数为0.100。在测试河北正定试验场包气带深部增厚区土壤水力参数的基础上,用非饱和达西法计算得包气带8.0-21.0m埋深范围内入渗补给速率为25-240mm/a。采用Hydrus1D软件并选取质量平衡法和达西法所涉及的钻孔资料建立模型,根据模型运行结果计算得到的补给速率、补给系数与场地实验评价的结果非常相似,证明一维非饱和数值模拟评价补给的可靠性。其中灌溉冬小麦-夏玉米轮作农田的补给速率、补给系数分别为48-89mm/a、0.053-0.096,不灌溉杨树林的补给速率、补给系数分别为12.3mm/a、0.024,不灌溉天然杂草试验场的补给速率、补给系数分别为83.4mm/a、0.165。本文采用氯离子质量平衡法得到的山前平原灌溉农田补给速率为65.9-126.8mm/a;采用人工溴示踪法得到的山前平原灌溉农田补给速率为54.5-466.4mm/a,中部平原灌溉农田补给速率介于11.9-238.8mm/a;采用一维非饱和流数值模拟计算的山前平原灌溉农田补给速率介于48.0-89.0mm/a。同时收集得到前人补给评价结果也大多介于此范围内,故本文的研究结果总体比较可信。(3)垂向入渗补给的影响因素①包气带变化对入渗补给的影响当包气带厚度在5m与10m之间时,石家庄地区冬小麦-夏玉米灌溉农田中三个钻孔根据质量平衡法计算的补给速率均比较稳定,没有太大的波动,即包气带厚度达到一定深度后,厚度继续增加,补给速率会保持稳定状态。根据人工溴示踪实验结果知华北平原山前平原和中部平原的平均补给系数分别为0.173、0.100。而山前平原多粉土、粘土,中部平原多砂壤土、壤土,说明岩性越粗,补给速率越大,反之亦然。因浅层水位持续下降,深部包气带砂性土发生一定的压密作用,使得土壤水力特性发生改变,即土体释水能力降低,进而影响到垂向入渗补给。此外,水位波动下降,若某-深度历史水位下降速率越慢,其土壤释水能力会更差,即水位下降速率大小对垂向入渗补给也有一定的影响。②外界因素对入渗补给的影响根据离子质量平衡法得出,对比灌溉冬小麦-夏玉米农田(DH01)和不灌溉杨树林(DH02),不同离子评价的补给,都显示DH01的补给速率要大于DH02,表明灌溉在一定程度上加大了水分入渗速率,即会增大补给速率。通过野外人工溴示踪实验,得出研究区内灌溉是产生垂向补给的主要来源,灌溉量越大,地下水补给速率也越大。研究区灌溉量普遍偏大,冬小麦-夏玉米轮作的方式在一定程度上减小了输入水量对地下水补给速率的影响；制定合理的灌溉制度,使降雨和灌溉水基本满足冬小麦-夏玉米等作物生长所需,不过度灌溉和开采,才能避免水位下降过快,在一定程度上保护地下水。基于一维非饱和流数值模拟模型,改变土地利用类型、保持其他模型输入条件一致,对比其对入渗补给的影响,模型中在冬小麦、夏玉米、棉花、天然杂草、杨树林的生长期间灌溉水量分别为260、130、130、0、0mm。结果显示,当只种植冬小麦时,补给速率最大；当冬小麦-夏玉米轮作时,其补给速率略小；当土地利用为天然杂草时,补给速率最小；综合分析得出当土地利用为研究区最常见的冬小麦-夏玉米轮作时,对地下水储存最不利。此外,还研究讨论了1960年至2010年期间的历史极端气候条件(降雨、气温)对补给的影响。结果表明,日最大降雨量对补给速率的影响要比持续多天降雨量大；即使持续降雨天数达到极端大值,若在此期间日最大降雨量较小或者该极端值所在年份的年降雨量很小,都不能影响到补给速率；极端高温或者极端低温对补给速率影响不大。主要创新点：(1)通过非饱和-饱和场地示踪实验,揭示非均质包气带三维水分及溶质运移规律；(2)基于质量平衡法,综合利用Cr、F、SO42等天然示踪剂评价垂向补给,多方法相互验证,提高评价结果的可靠性；(3)揭示了包气带厚度增大和结构变化对地下水补给的影响。
【作者】林丹；
【导师】靳孟贵；
【作者基本信息】中国地质大学，水文学及水资源，2014，博士
【关键词】华北平原；包气带；地下水补给；非均质；示踪；

【参考文献】
[1]郭祥.连续刚构桥沥青混凝土桥面铺装线性及非线性受力特性分析[D].北京交通大学,2014.
[2]周耀华.掌握区域成矿规律指导全省地质找矿——浅谈江西省第二轮成矿远景区划[J].江西地质科技.1995(01)
[3]王刚.面向PaaS应用的表单定制平台的研究[D].山东师范大学,计算机软件与理论,2013,硕士.
[4]杨松令,刘亭立.基于共生理论的上市公司股东行为研究——一个研究框架及设想[J].会计研究,2009,01:81-87+97.
[5]孙明蔚,童晓阳,刘新宇,罗孝辉,甄威,王晓茹.基于时序贝叶斯知识库的电网故障诊断软件原型系统设计[J].电力系统自动化,2014,21:77-83.
[6]刘元景.江苏省赛艇公开级优秀运动员全年训练负荷研究[D].宁波大学,体育,2014,硕士.
[7]赵菲.基于网络技术的新媒体发展模式研究[D].长安大学,科学技术哲学,2014,硕士.
[8]袁勇.双工件台宏微快速交接系统研究[D].哈尔滨工业大学,仪器科学与技术,2014,硕士.
[9]刘轶.基于切向喷射的气膜冷却的机理研究[D].哈尔滨工业大学,动力工程,2013,硕士.
[10]戴俊杰.初中数学职前教师错误分析能力现状的调查研究[D].东北师范大学,学科教学,2012,硕士.
[11]董菲.从翻译适应选择论看林译《迦茵小传》[D].四川外国语大学,英语语言文学,2013,硕士.
[12]杨再朝.应用迭代法进行位场的曲化平[J].石油地球物理勘探,1989,02:200-207+216.
[13]李心鑫.慢乙肝患者原位肝移植术前后Th17细胞变化水平以及其临床意义的研究[D].河北医科大学,外科学（专业学位）,2014,硕士.
[14]沙秀娟.阿德福韦酯、替诺福韦酯治疗CHB疗效与血药浓度（谷浓度）关系的研究[D].吉林大学,内科学,2013,硕士.
[15]白海灵.交互白板环境下师生交互行为个案研究[D].华中师范大学,教育技术学,2013,硕士.
[16]李阳.汽车起重机“2C2B”涂装工艺优化、“2C1B”涂装工艺研究[D].复旦大学,材料工程,2012,硕士.
[17]徐昊.NaCl溶液中沉淀硬化不锈钢的电化学行为[D].大连海事大学,材料科学与工程,2013,硕士.
[18]叶英,穆千祥,张成平.隧道施工多元信息预警与安全管理系统研究[J].岩石力学与工程学报,2009,05:900-907.
[19]邱峰.日常生活审美化[D].华东师范大学,美术学,2013,硕士.
[20]邓鸿鹄.北京市能源消费预测方法比较研究[D].北京林业大学,统计学,2013,硕士.
[21]周成香.有效加速系数在双参数指数分布中的应用[D].华东师范大学,概率论与数理统计,2004,硕士.
[22]王全收,刘洁慈,刘春友.2925磁带机电源原理及故障处理[J].石油地球物理勘探,1997,S2:178-185+190.
[23]罗华.论哈盖尔斯特罗姆的现实主义法学思想[D].复旦大学,法学理论,2012,硕士.
[24]李杰.基于地域特色的郊野公园规划设计研究[D].安徽建筑大学,城市规划与设计,2013,硕士.
[25]张楠.S集团预算管控问题研究[D].吉林大学,会计,2013,硕士.
[26]高杰.昌台牦牛遗传多样性蛋白和基因水平的研究[D].西南民族大学,基础兽医学,2014,硕士.
[27]王晓琳.常林机械集团员工薪酬激励体系研究[D].西北大学,工商管理（专业学位）,2012,硕士.
[28]胡晓丽.垂体瘤的临床病理诊断及分类诊断的价值评估[D].郑州大学,内科学,2011,硕士.
[29]高锋,吕中秀.矿山竖井信号微机检测与控制系统[J].冶金自动化,1996,01:43-45.
[30]贾慧星,章毓晋.基于梯度方向直方图特征的多核跟踪[J].自动化学报,2009,10:1283-1289.
[31]蔡葛峰.人力资本的企业会计核算的理论与实践探讨[D].东华大学,企业管理,2004,硕士.
[32]孟娜.超高层框架—核心筒结构静动力弹塑性分析与研究[D].河北工业大学,建筑与土木工程,2013,硕士.
[33]李大连.呼叫中心在报业的成功应用[J].中国传媒科技.2006(06)
[34]全再洪.新型硬质合金用成型剂CSZ的应用研究[D].湖南大学,材料工程,2014,硕士.
[35]李卓.电子加速器脉冲电子束流测量系统研究与设计[D].吉林大学,电路与系统,2013,硕士.
[36]耿杰.环氧树脂及其复合材料微波固化工艺研究[D].南京航空航天大学,材料加工工程,2013,硕士.
[37]高莹玉.α-苯乙胺的汽液平衡和拆分研究[D].郑州大学,应用化学,2013,硕士.
[38]肖漩.爱默生超验主义自然观研究[D].湖南师范大学,外国哲学,2013,硕士.
[39]武小涵.博弈论视角下我国粮食储备主体行为研究[D].华中农业大学,产业经济学,2014,硕士.
[40]胡彦斌.GdBaCo_2O_(5+δ)基氧电极的高温电学性质和表面偏析现象的研究[D].哈尔滨工业大学,凝聚态物理,2014,硕士.
[41]倪莎.基于AVR的数字高精度磁场测量系统的设计[D].华中师范大学,电路与系统,2013,硕士.
[42]陶华.ETC电子不停车收费系统的设计与实现[D].电子科技大学,软件工程（专业学位）,2012,硕士.
[43]崔莎莎.艺体类免费师范生从教意愿及其影响因素研究[D].华中师范大学,教育经济与管理,2014,硕士.
[44]生宁.基于液质联用技术的当药多组分分析与药代动力学研究[D].河北医科大学,中药学,2014,硕士.
[45]朱圣巧.公司监事责任制度研究[D].华东师范大学,经济法学,2013,硕士.
[46]刘金杭,樊阳.二氧化钛-石墨烯纳米复合材料的合成与电化学传感性能研究[A].河南省化学会.河南省化学会2010年学术年会论文摘要集[C].河南省化学会:,2010:1.
[47]张立伟,卢震宇,张全柱,邓永红.SS_(4B)型电力机车辅助供电系统的升级改造[J].机车电传动,2014,04:96-99.
[48]夏元清,黄一,许可康,韩京清.大射电望远镜FAST馈源舱位姿控制[J].控制与决策,2004,02:195-198+203.
[49]赵世震.基于Android平台和增强现实的导航应用设计与开发[D].北京邮电大学,电子与通信工程（专业学位）,2013,硕士.
[50]吴汝正.电缆屏蔽层对地故障点的探测[J].冶金自动化,1987,03:54-56.