山东淄博市大武水源地裂隙岩溶水中污染物运移的数值研究.pdfVIP

第 卷第 期 A 地学前缘（中国地质大学，北京） M6 9 A TM 9 年 月 !## % （ ， ） GH 9 !## FGH7I J5K4L54 )HML7K4HN IKLG ;LKO4HNK7P MQ -4MN5K4L54N .4KRKLS !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 山东淄博市大武水源地裂隙岩溶水中 污染物运移的数值研究 朱学愚，刘建立 （南京大学 地球科学系，江苏 南京!##$% ） 摘 要：在分析研究淄博市大武水源地裂隙岩溶含水层的水力性质和污染物运移特征的基础上，对裂 隙岩溶水的水头和污染物运移进行数值研究。目前国内外对裂隙岩溶水进行数值计算时，通常用等 价多孔介质模型，但裂隙岩溶介质和多孔介质有很大不同。裂隙岩溶介质的储水和导水空间为裂隙 网络，导水系数大，地下水的实际平均流速比孔隙水大得多，但给水度和贮水系数小。当用等价多孔 介质模型进行模拟时应考虑这些特点。对于污染物运移的模拟，要同时求解水头方程和对流弥散方 程，可采用’()*’+ 和,%( 软件进行模拟。研究区裂隙岩溶水水头的数值计算表明，等效多孔介 质模型水头的拟合误差能满足国标-./ ,00$1 2 $% 的要求。各时段地下水水量均衡计算的精度也满 足要求。对流弥散方程的数值计算，由于345647 数高达$8 9:1 ，对流占绝对优势，可能存在数值弥散和 数值振荡，因而采用多种方法进行了比较。对于同一问题，同时采用上游有限差分法（ ），混合的 ;)( 欧拉 拉格朗日方法（特征线法 ’、改进特征线法 ’ 和混合特征线法 =’ ），总变异消减法 （ ）进行计算，并比较其结果。结果表明，混合特征线法（ ）和总变异消减法（ ）比较适合于 ,( =’ ,( 对流占优势的运移问题计算。由于渗透系数 ! 和有效孔隙度 对溶质运移结果的影响很大，而对于 ! 裂隙岩溶含水层又难以准确确定，故用等价多孔介质模型进行裂隙岩溶含水层中污染物运移计算，其 精度要比水头计算的精度差些。 关键词： 裂隙岩溶水；污染物运移；对流弥散方程；数值计算；欧拉 拉格朗日方法 中图分类号： ； 文献标识码： 文章编号： （ ） ?8!% 3:0 @ ##8 !%! !## # #1 #A 山东淄博市大武水源地是一个特大型的水源 中，造成地下水的石油污染，在氯碱厂附近还有氯离 地，多年来平均开采量为 % ，供给 子污染。污染物向下游运移，威胁到大武水源地的 8!# ###80# ### B / C 该市张店和临淄两区8## ### 人口的城镇生活供水 安全，因此有必要对地下水中污染物的运移进行数 和若干大型国有企业的工业生产用水。含水地层为 值研究。 奥陶系中统马家沟组碳酸盐岩，裂隙岩溶发育，形成 一个统一的含水层。 世纪 年代中期，在水源 !# A# 研究区概况及水文地质条件 地西南部的补给区建设了%## ### 7 乙烯工程（包括 烯烃厂、氯碱厂、塑料厂等）。许多管道直接座落在 ! !