毕业设计PPT660125 0708.pptVIP

各种地质体的处理 (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。其主要原因可能是区域地质构造与岩溶较发育，模型概化不能精确描述实际复杂的水文地质条件。因此，模型可以用来模拟一个时间序列浓度连续变化趋势和估计具体时间点的浓度值。 模拟及预测流线 MODFLOW计算的稳定状态的地下水位等值线分布与实际地下水位分布规律一致。计算的泉S639的水位值和95年12月的实测值相1250m吻合较好。在等水位线上明显看出，F1，F2断层的阻水作用和F3断层的导水作用。 模拟及预测流场 水头赤泥库附近山区最高，向四周递减，猫跳河下游最低，地下水从山区流向河流、百花湖，水头与流线分布基本符合实际。 模拟及预测库区总碱度污染情况 400天的总碱度渲染图 第六章 结论与建议 (1)在扎塘赤泥库污染组分中，高碱度污染影响的范围最大，建议针对每一个可能产生的渗漏通道进行连续监测。并充分做好各种污染事故情景的应急预案。 (2)影响岩溶地下水流场及污染物运移的因素众多。研究区内岩溶发育的主要展布方向及断裂构造发育方向及性质对地下水的运移有较大影响；污染物运移主要受控于水动力条件。由于地区特殊的岩溶地质水文地质条件，建立模型前要尽量收集相关资料，对其进行概化时也需综合考虑，不可以偏概全。 (3)应用地下水渗流理论及水动力弥散理论建立了数学模型。用Visual MODFLOW中的MODFLOW2000、ModPath模块模拟了地下水渗流，用MT3DMS模块模拟了污染物运移，得到污染物等值线、浓度晕并进行了预测。对岩溶地区地下水模拟具有借鉴意义。 (4)验证模型参数时，应当尽量采用近期连续多年的长系列观测资料。因为早期的资料由于源汇项、地表条件以及降落漏斗范围等均发生变化，会影响所求得的参数的精度。当水文地质条件更改时，需要继续进行反复试算，以求较高精度。 本文主要得出以下结论与建议： (1)由于受限于时间、人力和财力，水文地质调查中只有少量泉点有水位、流量、浓度观测数据。虽然曾多次尝试人工调参和用MODFLOW中自动求参的模块进行校准，以期求得复杂的参数分布，再现实际情况中水文地质参数的非均质性，然而由于观测数据过少，所以自动求参结果一直不太理想，只能通过