米乐地下水环境影响评价中数值模拟的关键问题讨论

　　米乐建立三维地质模型后，将边界条件、水文地质参数和初始水头输入模型后，即可进行水流模型 的计算。数值法求解地下水非稳定流动问题需要 给出初始条件，即各节点初始时刻的水头。对于 网格剖分后形成的庞大数目的节点，实际的水位 观测数据显然无法满足需求。因此通常的处理方 法是利用稳定流计算得到模拟区的天然流场，依 据实际观测水位进行调参、反演和拟合工作以获 得能够反映地下水流场特征的初始水头。此过程 中的关键点在于边界条件及参数的输入、调参技 巧和参数反演与模型验证。

　　（1）边界条件及参数的输入。对于边界条件， 在实际操作中关键在于提高边界赋值的精度和质 量，以河流做定水头边界为例，在河道坡降比较小 且边界较短时可赋予同一定值，否则应沿河流进 行线性插值，ＦＥＦＬＯＷ 提供指定两点水头自动线 性插值的功能，在使用时应注意对于水头差较大 的情况宜分段进行插值，并对照节点的高程进行 检查，避免出现水头高于地表的情况。在参数输 入时，渗透系数、给水度和孔隙度等含水层特征参 数应使用野外试验和室内土工试验参数，在资料 不足时可先用经验参数代替。而降雨量、蒸发量 等源汇项参数可参考区域气象资料且精度应达到 多年月平均级别，避免使用常量。此外应根据水 文地质资料和实际野外调查情况合理的各类参数 进行分区，这是表达非均质性的主要途径，一般以 水文地质图为基础进行分区即可，但在图件资料比例尺较小或某评价区非均质性较强时，需特别 注重野外调查结果的综合分析。图５为湖南某工 程渗透系数分区情况，分区依据来源于野外调查。

　　（2）模型的校准米乐。由于地质体的非连续、非 均匀特点，模型中的各类参数远远无法表达一套 含水层的本身特性，在此情况下需要对模型进行 校准，通过参数的调整来寻找某个等效值，从而使 使用等效值的含水层整体特性逼近真实的含水层 特性，即模型校准是一个调整模型输入参数，直到 模型输出变量（或因变量）与野外观测值达到适当 匹配程度的过程，一般分为正演和反演两大类。在实际工作中正演试错法是最常用的方法，它并 无通用的执行标准，对调参经验和水文地质认识 要求高，需要统观模型的各个方面，当参数组合较 多时整个过程十分冗长复杂，本文仅根据经验就 水流模型的校准总结如下：① 调参准备。明确待 校准的参数数量及调整范围，主要对象是一般入 渗补给量和渗透系数，其中渗透系数是关键。此 外还应设置好合理的水位观测点，详细记录调整 过程、各观测点的水位变化及模型水均衡情况。② 调参过程。整体过程分为前中后三期。前期 可使用单因子调参原则和大跳跃的参数间隔，即 每次只改变某一参数且调整幅度可成倍甚至成数 量级变化，主要目的是熟悉模型基本情况，宏观掌 握各参数的敏感性，确定重点调参对象；此后即可 进入调参中期，可同时调整两个甚至三个参数，确 定它们彼此间的相对大小关系，调整幅度一般不 出现数量级变化，主要任务是在明确各参数相对 关系并缩小参数调整值域；最后进入调参后期，在 维持参数间相对关系的前提下进行精细化微调米乐， 调整幅度以参数值的１０％甚至更低为基准线，以 求得可靠的参数组合。需要注意的是，整体趋势 的拟合比单点的拟合精度更为重要。③ 分析思 路。在调参过程中确定调整方向是关键，应认真 分析每个观测点的变化情况，从补给量和渗透性 入手确定调整方案，如在第四系孔隙水仅接受大 气降雨补给时，若孔隙水模拟水位偏低可考虑增 大降雨入渗补给、调低第四系垂向入渗或调低下 层含水层渗透性；在第四系孔隙水与河流边界有 水力联系时，若孔隙水模拟水位偏低则可考虑增大转化系数。此外可依据实际情况对地质模型、 网格剖分、边界条件和参数分区等做合理的调整。④模型设置。模型中模式设置对计算结果影响极 大，其中最关键的是潜水面模式的设置。ＦＥＦＬＯＷ 对潜水面的处理方式有 Ｆｒｅｅ和 Ｐｈｒｅａｔｉｃ 两种。其中Ｆｒｅｅ模式只模拟饱和部分，适用于地 下水水位只在某一层内变化的情况，对潜水面陡 变的情景不太适用；而Ｐｈｒｅａｔｉｃ模式允许单元变 干或部分饱和，适用于地下水位梯度较陡且跨越 多层的情景。值得注意的是，完全适用于Ｆｒｅｅ模 式的情景较为少见，Ｐｈｒｅａｔｉｃ模式同样可用于 Ｆｒｅｅ模式下问题的求解，虽然计算量大、收敛相 对困难，但其计算精度高且收敛后的水均衡情况 好，可优先考虑使用。此外关于是否约束水位在 顶底面之内的设置，应根据实际情况进行动态调 整，如模拟抽水情景时，某一层的水位低于其底板 属正常情况，则不应进行约束。关于参数的反演， 实质上是多目标函数求最优解问题，目前国际上 已有若干水流模型自动参数识别程序，如 ＦＥ－ ＦＬＯＷ 自带的 ＰＥＳＴ 程序。需要注意的是自动 调参是参数数值上的调整并不涉及参数结构，一 开始就使用很难获得满意的结果，而在正演的中 后期，进行参数微调时使用效果较好。

　　（3）模型验证。当完成校准工作后需对结果 进行验证，主要原因是模型校准是一个非唯一性 的过程，即很多参数组合可能显著不同，但均能提 供与观测值同等合理匹配的模拟结果。一般应从 水位、水动态和水均衡三方面进行验证。然而在 实际项目中，受数据限制基本上仅开展了基于水 位的验证工作，但对非稳定流模拟来说，仅用稳定 流数据校准是远不够的，当未充分验证的模型加 载各类水量、水质工况并进行非稳定流计算时，模 型可能出现水动态异常、水均衡失稳、溶质浓度出 现较大负值、溶质迁移规律明显错误和模型求解 超出误差范围与迭代次数等各类问题。故在资料 偏少的情况下，根据专业经验针对水动态和水均 衡开展半定量甚至定性的验证工作仍十分必要， 可虚拟一定数量的水位观测点，记录各类边界条 件、各含水层、重要区域和模型整体的水均衡情 况，调整模型为非稳定流并输入动态降雨数据，运行 若干年并对水动态和水均衡数据进行专业判读。

　　弥散度是研究污染物在土壤及地下水中迁移 转化规律的重要参数之一，弥散度随着溶质运移 距离和研究问题尺度增大而增大的现象称为多孔介质水动力弥散的尺度效应。对于造成水动力 弥散尺度效应的原因，目前趋于一致的看法是野 外条件下介质的不均匀性造成了室内试验结果与 野外试验结果之间的巨大差别。目前一般不推荐开展弥散试验工作，弥散度 参数可参照ＧｅｌｈａｒＬＷ 等对世界范围内所收 集的５９个大区域弥散资料进行的整理分析结果 确定模型所对应的参数级别，在风险最大化原则 下进行取值。

　　关于工况情景的设计，考虑防渗和工况两种 因素将其分为四种情景，见表２。其中情景１、３ 模拟意义相似，但情景１可在一定程度上验证长 期渗漏的情况下防渗层的可靠性，故一般选择情 景１进行模拟，此外情景２、４分别代表长期低量 和短期高量污染下的环境风险，应作为重点模拟 情景。

　　将各类工况情景中的污染泄漏地点、污染物 类型、浓度和总量输入溶质模型后即可执行计算米乐， 而污染物总量的控制和污染源强的刻画是其中的 关键问题。污染物总量的控制有两种实现方式：

　　①方法 １。可采用修改入渗补给量并设置定浓度边界的 方法，即将泄漏区每天的入渗补给量修改为污染 物每天泄漏的体积，并将泄漏区设置为定浓度边 界，浓度取污染物在污水中的浓度。

　　②方法２。采用设置注水井并施加注入浓度的方法，即将污 染物每天泄漏的体积换算到泄漏区的每个节点 上，将所有节点设置为注水井并指定注入浓度为 污染物浓度。

　　上述两种方法通过控制总体积和浓 度来刻画污染物总量，值得注意的是，两种方法均 有一定的缺陷和适用范围。在不考虑地面工程场 地的固化情况下，方法１改变了泄漏区的入渗补 给量，当泄漏区面积相对于模拟面积不可忽略时， 需慎用方法１。而对于方法２来说，面状污染源 的节点数较庞大，全部设置为注水井后将对水流 模型的求解产生较大影响，此方法对网格的剖分、 注水井个数及注水量等问题有着较高的设计要 求。在实际工作中，应根据具体情况进行选择，必 要时可进行多次试算来辅助判断。

　　对污染源强来说，将污染区设置为定浓度边 界后即为连续源强，而间断源强的控制方法是利 用时间序列加上边界限制条件米乐。需注意的是时间 序列中的值代表的是赋值操作，并非为进入地下 水中的污染物浓度值，故需要施加限制条件以解 决此问题，即限制泄漏区范围内地下水中污染物 浓度最小值为零以使得赋值操作失效。此外在间 断时间点前后插入间隔较小的时间点，有助于模 型的求解。

　　特别声明：北极星转载其他网站内容，出于传递更多信息而非盈利之目的，同时并不代表赞成其观点或证实其描述，内容仅供参考。版权归原作者所有，若有侵权，请联系我们删除。

　　摘要：当前阶段我国地下水已经污染的较为严重，很多城市的地下水都普遍遭受到了污染。因此，如何解决地下水污染问题已经逐渐成为当前研究的主要课题。本篇文章全面分析了当前阶段地下水污染的现状，并阐述了地下水污染的概念及其评价情况，最后详细给出了关于地下水污染问题的解决对策。关键词：地下水

　　地下水系统是什么？本篇就来讲解一下地下水系统和地下含水系统以及地下水流动系统，详情如下：

　　摘要：岩溶地下水资源是西南地区社会和经济发展的重要保障，随着城市化进程的不断加快，西南地区岩溶地下水污染问题也日益凸显，直接威胁当地群众的健康安全，已成为社会热点。以西南地区典型岩溶地下水系统污染实例为研究对象，通过分析每个岩溶地下水系统污染实例的水文地质条件、污染源特征、污染物

　　一、微纳米气泡“云”透析池塘水污染治理技术1、核心技术简介池塘水污染治理需靠人工创造较好的水生态基础环境。其治理方法由三个装置、一个系统组成。三个装置为：一是：控制水的流动装置，即推流装置；二是：超微米气泡透析装置；三是：浮渣污染物收集装置。一个系统,即水生态立体系统。其技术路线为

　　12月8日，在河北涿州市刁窝镇，涿州市供电公司发展建设部员工崔光凯等人参考“电网水文地质信息一张图”（以下简称“一张图”），对10千伏刁豆线防洪防涝改造工程方案开展研究。借助“一张图”，电网规划设计人员能够准确掌握重要水文地质信息，更加科学高效地做好受损电网恢复重建方案编制。在以往电

　　在水文地质中，根据地下埋藏条件的不同，地下水可分为上层滞水、潜水和承压水三大类。下面就给大家通俗易懂的讲解一下他们的区别和联系。01上层滞水上层滞水(perchedwater)存在于潜水面以上包气带中的局部不透水层或弱透水层之上的重力水。上层滞水的形成是在大面积透水的水平或缓倾斜岩层中，有相对隔

　　摘要：一直以来，污染场地环境治理都是一项难度较高的工作，归根究底在于所用技术缺乏一定的合理性，难以保证勘察结果的准确性，无法被应用到环境治理工作中，给相关工作发展造成了不小的阻碍。新时期的来临诞生了大量的科技含量较高的新型技术，倘若能够将这些技术运用到污染场地环境治理工作中，必然

　　水文地质、工程地质和环境地质条件是矿山开采的主要基础技术条件，不但对矿山建设、开采和生产起重要控制作用，同样对矿山生态修复的工程设计起着重要的限制作用。一、微地貌及水文调查（一）微地貌调查矿区所处地貌单元多为二级或三级，要从地貌形态入手，调查矿区所处的原始微地貌类型，以及人工地貌

　　作者：丁爱中北师大水科学研究院教授。主持/完成国家、部(省)级科研项目30多项，包括国家“十一五”科技支撑重大、重点项目子课题、“863”重大项目子课题、北京科委重点项目、“十五”重点项目、国家科委攻关项目米乐、水利部重大科技项目、国家自然科学基金项目、广东省自然基金和攻关项目等。01地下水基

　　摘要：伴随着国家矿产资源开发的黄金十年，铁矿的资源开采被提上了日程，但开采的先决条件即解决矿山开采与地下水资源保护和安全建设之间的矛盾。关键词：大水矿山；帷幕注浆；地面预注浆；矿坑水回灌；治理；在基建及生产过程中，坑下开拓工程揭露灰岩顶板时均出现大量涌水，并造成地面塌陷、地表水倒

　　天津地区在沉积类型上有显著的横向分带特点，其中滨海沉积地层占据了天津市总面积的四分之三，在水文地质方面有不同于非滨海沉积区的显著特征。中科环境修复（天津）股份有限公司周胜教授针对这些特征如何影响污染场地调查及修复，在首届天津土壤污染防治高峰论坛就“天津地区污染场地环境调查中水文地

　　在系统调查场地是否被污染及污染程度和范围的过程中，水文地质调查越来越多地得到使用。本文主要介绍场地污染调查的现状及存在问题，探讨场地水文地质调查的必要性，进而为污染场地修复提供依据。延伸阅读：【PPT】扒一扒污染场地环境修复那些事关于我国污染场地风险控制和修复技术决策的几个观点[$New

　　在系统调查场地是否被污染及污染程度和范围的过程中，水文地质调查越来越多地得到使用。技术人员通过借鉴地质行业中已然规范成型的水文地质调查技术方法，并将其与环境科学的专业有机结合，不仅相当程度上丰富了当今国内的污染场地调查技术手段，同时也体现了行业内部重视水文与环保相结合的一致共识。

　　因群众反对意见较大，14日惠州市政府新闻发言人表示，广东“惠州市生态环境园项目”选址尚在论证阶段，又一个垃圾焚烧厂陷入“宣布-反对-搁置”窘境。一些公共环境项目“一闹就停”的背后，折射的是垃圾焚烧技术与成本、知情权、选址等多重矛盾。根据《“十二五”全国城镇生活垃圾无害化处理设施建设规划》要求，“十二五”期间全国需要新增每日20万吨的垃圾焚烧处理能力，垃圾焚烧厂“落地”难题如何解?“上马—