邹城唐村—西龙河水源地岩溶水资源数值模拟.PDF

第５４卷　 第２期　　 地 　 质　 论 　 评 　　　 ＧＥＯＬＯＧＩＣＡＬＲＥＶＩＥＷ　　 Ｖｏｌ．５４　Ｎｏ．２ ２００８年 ３月 Ｍａｒ．　２００８ 邹城市唐村—西龙河水源地岩溶水资源数值模拟 曹丁涛 兖州煤业股份有限公司，山东邹城，２７３５００ 内容提要：邹城市唐村—西龙河水源地位于郭里集单斜东缘。该单斜区已建有多个水源地，目前处于一种供稍 大于求的弱平衡状态，正确评价唐村—西龙河新建水源地开采量十分重要。在深入分析研究区及周边水文地质条 件基础上，分别按照平水年和枯水年设置开采方案，采用ＶｉｓｕａｌＭＯＤＦＬＯＷ地下水流数值模拟模型软件，对唐村—西 龙河水源地岩溶水资源进行了精细的数值模拟，经模型识别和检验之后，对水源地岩溶地下水可开采量作出了评 价，通过分析比较，选定了方案三。在水位约束评价时，考虑了实际存在的“井损”问题，因而，评价结果更趋于客观、 实际。模拟结果可以作为水源地岩溶地下水可开采量的依据。 关键词：岩溶水资源；数值模拟；“井损”；邹城唐村—西龙河水源地 　　２０世纪６０年代以来，随着计算机技术的迅猛 地下水数值模拟也存在一些问题有待改进和提 发展，数值方法在地下水研究领域得到了快速发展。 高。如野外试验获取的第一手资料，也是至关重要 数值模拟已成为地下水资源定量计算与评价的重要 的资料，不足或代表性不强；各有关学科之间缺乏沟 手段。目前，国外主要针对数值模拟在某些方面的 通（郝治福等，２００６）；随机模拟研究滞后；基础理论 薄弱环节，提出一些新的思维方法，采用新的数学手 的实验研究重视不够；轻视具体地质条件的研究，过 段，分析不同尺度下的变化，从而合理地描述地下水 多地依赖模拟技术；模型建了很多，但投入运行后的 系统中某些不确定性或模糊因素。Ｗｏｏｄ（１９９６）提 实际效果尤其长期效果如何，却很少有人问津（薛 出了二维地下水运动有限元计算的时间步长条件。 禹群等，１９９９）。此外，水资源计算与评价时，忽视 Ｓｃｈｅｉｂｅ等（１９９８）模拟分析了不同尺度下的地下水 “井损”因素也是较为普遍的现象。尽管“井损”现 流及其运移行为。Ｍａｚｚｉａ等（２００２）提出一种特别 象早就为水文地质工作者发现并进行了研究 的数值方法求解重盐地下水运移模拟的二维非线性 （Ｐｒｉｃｋｅｔｔ，１９６７；Ｐｅａｃｅｍａｎ，１９７８；李祖贻等，１９８４；曹 动力学控制方程。Ｍｅｈｌ等（２００２）提出了二维局部 洪等，２００３），但在实践中，尤其在大型水源地水资 网格细分法的有限差分地下水模型，提供了新的插 源评价时，却往往又把井损抛开了。由于井损的存 值和误差分析方法。ＬｉＳｈｕｇｕａｎｇ等（２００３）开创性 在，抽水井的实际水位要比模型计算的理论水位低 地提出了一种随机地下水模型，解决均值分部和小 许多，常常导致预测、评价结果与实际差距较大，甚 尺度过程的不同尺度问题。我国的水文地质学家也 至达一倍以上，给政府或企业造成很大被动或损失。 做了卓有成效的探索、研究工作。卞锦宇等（２００２） 故在唐村—西龙河水源地资源评价中，特别注意考 较好地解决了相对隔水层缺失区越流系数无法调试 虑了井损因素，评价结果更趋客观、真实。 的难题。张祥伟等（２００４）根据地质统计、逆问题理 １　研究区基本地质概况 论和地下水运动理论提出了大区域地下水系统数值 模拟的理论和方法。薛禹群等（２００４）将 ＭｓＦＥＭ 研究区位于邹城市（原邹县）南部郭里集单斜 东缘，面积约 ２ （多尺度有限元法）应用于非均质多孔介质中的流