1 - 济南泉水流量衰减过程及复流条件的物理模拟济南泉水流量衰减过程及复流条件的物理模拟 季叶飞1 ，邹靖1 ，顾锦1 ，张维1 1.河海大学水文水资源学院，江苏南京 （210098） E-mail：jiyefei@ 摘 要： 摘 要： 为了研究济南泉流量衰减的物理机理， 本文运用物理模型对泉水的成因及其衰减过程进行模拟 由于不同时期泉流量动态的影响因素不同， 设计了三个方案对不同时期的泉流量动态变化过程进行模拟，模拟结果与实际情况相符同时在泉流量动态模拟的基础上，设计了两个使泉水复流的方案，并根据模型的模拟结果，提出了保持泉水常年喷涌的建议 关键词：关键词：济南泉水、物理模拟、流量衰减、泉水复流 山东省济南市处于我国北方碳酸盐岩岩溶发育地区[1]，水文地质条件复杂，生态环境脆弱， 在人类活动的影响下， 尤其是大量地开采地下水资源， 导致区域性地下水位持续下降、水质恶化、泉水流量衰减甚至断流等一系列环境问题[2]济南市是一个资源型缺水城市，水资源开发利用程度很不均匀，地表水开发程度很低，地下水开发程度较高，特别是市区及其周边地区，部分地区已出现超采地下水超采不仅会改变地下水本身的水动力场，形成区域性的地下水位降落漏斗，而且会袭夺地表径流，改变地下水与地表水的补排关系[3]。
济南泉水从 1972 年春季开始出现断流，以后几乎年年有停喷现象，而且停喷时间逐年延长[4]，由于过量开采地下水，漏斗区面积已扩大到 289 km2泉流量大小与附近地下水位高低成正比，也与附近含水层渗透性能有关[5] 研究济南泉水的衰减原因对济南市的保泉供水及其它地区类似问题的研究有着极其重要的参考价值， 本文通过物理模拟实验对济南泉域泉水流量衰减的影响因素进行分析，并提出了相应的解决方案 1.研究区简介 1.研究区简介 济南市位于鲁中山地西北，地理坐标为东经 116°11'-117°44'，北纬 36°01'-37°32'，面积 8154 km2济南市东临环渤海经济圈，西临山西能源基地，北靠京津唐城市带，南接沪宁杭经济区，是我国东部沿海经济开发区重要的中心城市之一，也是华北与华东、沿海与内地的陆路交通枢纽，地理位置十分优越 济南泉域为一相对闭合流域， 地处中纬度内陆地带， 属温带大陆性气候， 春季干燥多风，夏季炎热多雨，秋季晴爽，冬季寒冷干燥多年平均气温 14.20℃，多年平均降水量为675.31mm[6]泉域总面积 818.5 km2，其地域范围为：东到埠村向斜构造(章丘地界)，南到莱芜盆地、 锦绣川以北的区域分水岭， 西到济南火成岩体(隔水)西南部及与之相连的刘长山——万灵山岩溶弱发育带， 北到济南阻水的火成岩体。
济南市区的泉水依其地理分布大体可基金资助：河海大学大学生创新基金项目 - 2 - 分为四大泉群：趵突泉泉群、五龙潭泉群、黑虎泉泉群和珍珠泉泉群它们集中分布在东起青龙桥，西至匡市街，南起正觉寺街，北至大明湖这一面积约 2.6km2的范围内（图 1）[7] 图 1 济南岩溶泉水分布图 图 1 济南岩溶泉水分布图 济南市南部为泰山山地，北部为黄河冲积平原，地势南高北低，地质构造总体上是一个以古生代地层为主体的北倾单斜构造，大气降水在南部山区以直接和间接的方式补给地下水 南部山区的大气降水沿裂隙入渗， 地下水沿地层倾向流到市区北部受到岩浆岩体的阻挡，沿相对封闭的顶板裂隙上升喷出地表而形成泉水 图 2 济南泉水成因示意图 图 2 济南泉水成因示意图 注：引自山东省地矿局第一水文地质队资料 由济南泉域 1959-1987 年的降水量、趵突泉流量和地下水开采量的变化过程线（图 3）可知：1968 年是一个明显的转折点，1968 年以前地下水开采量较小，小于泉流量，泉流量只与降水量存在较高的相关关系；1968 年以后，地下水开采量迅速增大，大于泉流量，泉水流量急剧下降且与地下水开采量、降水量变化密切相关。
- 3 - 020040060080010001200140019591963196719711975197919831987年份降水量（毫米/年）0102030405060地下水开采量（万立方米/日）降水量泉水流量地下水开采量图 3 1959～1987 年济南市降水量、趵突泉泉流量、地下水开采量变化过程线 2.物理模拟 2.1 模型简介 图 3 1959～1987 年济南市降水量、趵突泉泉流量、地下水开采量变化过程线 2.物理模拟 2.1 模型简介 本实验器材由英国 Armfield 公司生产的 BHS 系统（Basic Hydrology System） ，其结构示意图见图 4该系统主要由以下几个部件构成：降雨装置、测压管、含水介质模拟区、出水口和开采井等降雨系统主要为实验提供不同的降雨条件；测压管用于测量地下水头；含水介质模拟区用于填放实验的含水介质（如岩块、粗砂、中砂等）以模拟岩石中的空隙；出水口用于模拟泉水出流；开采井则用于控制地下水的开采量 该实验器材的主要功能包括： （1）模拟降雨与径流之间的关系； （2）模拟地下水流及其影响因素之间的关系等第一个功能分为六个子功能，第二个功能分为 8 个子功能，在此不一一列举。
本文主要运用该系统的第二个功能， 对泉流量的动态变化及其与影响因素之间关系的模拟 图 4 BHS 结构示意图 2.2 裂隙介质的选定及其排列方式 图 4 BHS 结构示意图 2.2 裂隙介质的选定及其排列方式 由于裂隙岩溶水介质的复杂性， 要建立实验室模型较为困难， 尤其是设计出近似于实际含水层的实验室模型更加困难[8]为了使物理模型与实际的地质条件相似，选择石灰岩作为含水介质的模拟材料，其为规则的长方体，长 30cm、宽 15cm、厚 2cm - 4 - 实验槽的规格为长 2m、宽 1m，高 0.2m介质的局部排列方式如图 5 所示在实验仪器中共铺设六层石灰岩块， 用石料间的缝隙来模拟实际的裂隙 经实验测得裂隙模拟区的裂隙率为 0.25 图 5 模拟的裂隙介质局部排列方式 图 5 模拟的裂隙介质局部排列方式 2.3 方案确定 2.3 方案确定 根据上述分析， 为了更好地了解泉流量的动态变化， 拟定了以下的三个模拟方案 （表 1） ： 表 1 泉流量衰减过程模拟方案 表 1 泉流量衰减过程模拟方案 方案 具体的实施方法 1-1 将降雨量从 0.14L/s 降到 0.10L/s，无地下水开采 1-2 将降雨量从 0.14L/s 降到 0，无地下水开采 1-3 将地下水开采量从 0 增大至 0.035L/s，降水量维持 0.10L/s 不变 注：方案 1-1 中的降雨量从 0.14L/s 降到 0.10L/s 之前降雨量一直维持在 0.14L/s，且泉的出流处于一个相对稳定的状态，方案 1-2，1-3 中均是如此。
2.4 模拟结果 方案 1-1：2.4 模拟结果 方案 1-1： 将实验结果点绘成图 （图 6） ， 由图 6 可明显看出： 在无地下水开采的情况下，泉流量的大小受降雨量控制， 且呈正相关关系 泉流量在 215s 内从 0.13L/s 衰减为 0.09L/s，之后便一直维持在 0.09L/s，即泉流量 1s 内衰减 0.00019L/s - 5 - 图 6 泉流量与降雨量的关系曲线 方案 1-2：图 6 泉流量与降雨量的关系曲线 方案 1-2： 将实验结果点绘成图 （图 7） ， 由图 7 可明显看出： 在无地下水开采的情况下，泉流量受降雨量控制，且呈正相关关系泉流量在 250s 之内，由 0.13L/s 降为 0 即断流，泉水流量 1s 内衰减 0.00052L/s降雨停止后，泉水流量迅速衰减，衰减速率约为方案 1-1的 3 倍 图 7 泉流量与降雨量的关系曲线 方案 1-3:图 7 泉流量与降雨量的关系曲线 方案 1-3:实验结果如图 8 所示，由图 8 可看出：在降雨量恒定的情况下，地下水开采后，泉流量明显减小，泉流量受地下水开采控制，且呈现负相关关系 - 6 - 图 8 泉流量与地下水开采量、降雨量的关系曲线 图 8 泉流量与地下水开采量、降雨量的关系曲线 泉流量衰减过程的模拟有利于保泉工作的开展， 泉水复流模拟则能更好的指导泉水复流工作的开展。
3.泉水复流方案的模拟分析 3.1 方案确定 3.泉水复流方案的模拟分析 3.1 方案确定 根据上述的泉流量衰减的物理模拟分析， 本文拟定了以下的方案来实现泉水的复流 具体的方案见表 2 表 2 泉水复流条件的模拟方案 表 2 泉水复流条件的模拟方案 方案 具体的实施方法 2-1 将地下水开采量从 0.09/s 降到 0.03L/s，同时保持降水量为 0.06L/s 不变 2-2 将降雨量从 0.06L/s 增大到 0.08L/s，维持地下水开采量为 0.06L/s 不变 注：方案 2-1 中，在降雨量增大到 0.08L/s 之前，降雨量一直维持 0.06L/s，地下水开采量也维持 0.06L/s不变，同时，泉流量处于一个相对稳定的状态方案 2 也是如此 3.2 模拟结果及其分析 方案 2-1：3.2 模拟结果及其分析 方案 2-1：将实验数据点绘成图（图 9） ，由图 9 可明显看出：当泉断流后，如果及时地减小地下水开采量，泉水就有复流可能 - 7 - 图 9 泉流量与地下水开采量、降雨量的关系曲线 方案 2-2：图 9 泉流量与地下水开采量、降雨量的关系曲线 方案 2-2：将实验数据点绘成图（图 10） ，由图 10 可很明显看出：当泉断流后，如果及时地增大一定数量的降水量，如丰水年、或人工补给，那么泉水就会有复流的可能。
图 10 泉流量与降雨量、地下水开采量的关系曲线 4.物理模拟的结论 图 10 泉流量与降雨量、地下水开采量的关系曲线 4.物理模拟的结论 根据上述的模拟结果可得出以下的结论： 泉流量受降水量与地下水开采量的控制， 且在不同的情况下，这两个因素中哪个占主导位置是不同的即在降雨丰富、地下水开采水平较低的情况下，降雨量是主导控制因素；枯季、地下水开采量很大的情况下，地下水开采量为主要控制因素 因此， 可通过人工补给或限制地下水开采这两种积极的措施使泉水尽快复流将上述实验方案及结果运用到实际可以得到以下结果： 在降雨量充沛的季节 （雨季） ， 泉基本能够一直保持喷涌状态； 在降雨量缺乏的季节 （枯季） ，要使泉保持喷涌状态应开展地下水人工补给或减小地下水开采量，使地下水位抬升， - 8 - 从而使泉水复流 济南市区土地硬化面积的增加以及南部山区的工程建设， 对水土流失影响很大， 减少了泉城的地下水入渗补给[9]在补给区内实施人工回灌是恢复泉涌的重要措施之一[10] 5.致谢 5.致谢 感谢河海大学水文水资源学院束龙仓教授在整个实验过程和论文完成过程中的悉心指导，感谢硕士研究生王茂枚在整个实验过程中提供的帮助。
参考文献 参考文献 [1] 汪家权，吴义锋，钱家忠，李福林． 济南泉域岩溶地下水三维等参有限元数值模拟[J]， 煤田地质与勘探，2005（3） ：39-41 [2] 李砚阁，杨昌兵，耿雷华等．北方岩溶大泉流量动态模拟及其管理[J]，水科学进展，1998（3） ：275-281 [3] 张效平，宋儒．洪山泉流量动态及影响因素分析[J]，中国煤田地质，2002（3） ：31-32 [4] 章亦兵．济南市人工回灌补源保护泉水的研究[D]，南京：河海大学，2005 [5] 卢文喜．地下水系统的模拟预测和优化管理[M]，北京：科学出版社，1999 [6] 。