2014年扎塘赤泥库地下水渗流及污染运移数值模拟.pptVIP

2 溶质迁移方程的建立 * 我的中国梦 我的理想 我的未来 2014/4/19 2 建立数学模型 渗流模型： * 我的中国梦 我的理想 我的未来 2014/4/19 渗流模型： * 我的中国梦 我的理想 我的未来 2014/4/19 第五章 基于Visual Modflow的计算机模型的建立 * * 计算机模型的建模过程 计算区域的剖分图 * 计算区域的剖分图 * * 模型初始参数 (1)含水层渗透系数 根据区域特点，将含水层概化为均质各向异性潜水含水层，补给强度134.8mm/y。渗透系数Kx=0.0001m/s，Ky=0.0001 m/s， Kz=0.00001。贮水率μs设为0.00005，给水度μ设为0.2，空隙度设为0.3，有效空隙度设为0.15。 组分取碱度，线性等温吸附，分解系数为0.0000001L/mg，无反应项。弥散度为37.2m，横向/纵向弥散度比例0.1。 (2)给水度 根据本地区含水层岩性及颗粒组成，参照本区的前人研究成果和《水文地质 手册》中的经验值，取给水度μ=0.2-0.3 。 (3)降雨入渗补给强度 根据本地区包气带岩性及地貌特征，参照《水文地质手册》取降雨入渗补给 强度134.8mm/y。 (4)河床底积物渗透系数 根据野外的竖管渗水试验资料来计算，得到河床底积物的垂向渗透系数为0.0001m/s。 (5)其他系数 空隙度设为0.3，有效空隙度设为0.15。组分取碱度，线性等温吸附，分解系数为0.0000001L/mg，无反应项。弥散度为37.2m，横向/纵向弥散度比例0.1。 * 初值 最小 最大 渗透系数k Kx (m/s) 0.0001 0.00001 0.001 Ky(m/s) 0.0001 0.00001 0.001 Kz(m/s) 0.00001 0.000001 0.0001 储水率 0.0005 0.00001 0.0009 给水度 0.2 0.01 0.8 * 各种地质体的处理 (1)源汇项处理 补给主要为大气降水补给，区域相对较小，故将整个区域的补给设为统一的。 排泄项主要有向麦架河、猫跳河、及百花湖水库排泄，以泉的形式排泄。其中的泉处理为抽水井，流量根据实测资料输入。向河流、湖排泄及通用边界的水量交换模型自动计算。 浓度补给主要靠赤泥库区废水渗漏，废水碱度在14g/L左右，将其设连续浓度补给。 同时在牟老泉S010、白水泉S313、陆家湾S407、老马寨S639处设浓度观测井。 (2)断层处理 F1、F2断层位于库区西面的南北向逆掩断层，F1断层与F2断层相距50—100m，他们之间夹有一套龙潭组煤系地层，渗透系数较小。将其处理为墙，厚度0.5m，渗透系数K设为0.000001m/s。F3断层是一条北西向平移断层，倾向北东，倾角陡立，走向120°—160°，水平错距达170m，该断层东起3号库西侧，经陆家湾丫口向北西一直延伸至长冲S313泉一带。将其处理为排水沟，高程1200m。 * 运行参数的设置 (1)解算器的选择与设置 在MODFLOW2000中利用WHS解算器计算水流。 在MT3DMS中选择常规共扼梯度(GCG)解法和上游有限插分(Upstream Finite Difference)计算方法。上游有限插分方法为溶质运移模型在一个相对较短的方法来计算有限差分方程，该方法比外迭代方法更快。乘数取1.1，最大运移步长200d。输出控制中，模拟时间7300d，最大运移步长10000，指定输出时间分别为180、210、240、270、300、330、360、390、420、450、480、510、540、570、600天。 (2)时间步长 根据研究区的已知资料，选取1990年6月到1991年6月作为模型识别阶段，每月设为一个应力期，每个应力期步长设为200步。 * 泉407实测浓度与预测浓度拟合图 * 泉639实测浓度与预测浓度拟合图 * S407、S639两观测点的浓度平均误差分别为-1.3e-12mg/l，-1.6e-12 mg/l，可见正负误差相当均衡，S407相关系数平均为0.92，S639相关系数平均为0.88，系统稳定性较好，表明所建模型对水文地质条件的概化、边界条件的确定都与研究区域实际情况基本吻合，但单点实测值与模拟值有存在误差，S407的最大误差为2.699g/l，最小误差为-0.084g/l，S639的最大误差为2.008.91g/l，最小误差为0.04259g/l。其主要原因可能是区域地质构造与岩溶较发育，模型概化不能精确描述实际复杂的