矿渣堆积对河网水质影响的数值模拟分析.pdf

收稿日期：基金项目：国 家 自 然 科 学 基 金 资 助 项 目（，）作者简介：李文钰（），女，河南安阳人，硕士研究生，研究方向为水力学及河流动力学 通信作者：侯精明（），男，河北怀安人，教授，博士，研究方向为城市水利学：【水环境与水生态】矿渣堆积对河网水质影响的数值模拟分析李文钰，侯精明，刘占衍，张 松，栾广学，杜颖恩，韩占涛（西安理工大学 西北旱区生态水利国家重点实验室，陕西 西安；河北省邯郸水文勘测研究中心，河北 邯郸；四川水利职业技术学院，四川 成都；河南智河工程技术有限公司，河南 郑州；生态环境部 土壤与农业农村生态环境监管技术中心，北京）摘 要：掌握矿山区污染物迁移转化过程对其流域河网水环境的有效治理至关重要。以十堰市小文峪河流域为例，基于水动力水质耦合模型，充分利用各支流矿渣堆属性、流域河网本底特征及地表水、矿洞水、淋溶水的水质和流量监测数据等，构建了基于 加速技术的矿山区河网水环境高分辨率数值模型，对小文峪河流域硫酸根离子和铁离子的空间分布进行了数值模拟，并量化分析了小文峪河污染负荷来源。结果表明：所构建模型对矿山区河网水环境模拟效果较好，可以表征硫酸根离子和铁离子空间分布特征及支流与主河道交汇处污染物浓度变化过程；各支流硫酸根离子和铁离子对主河道贡献率由大到小依次为宝泉寺、郭家沟、竹园沟、清石沟、孙家沟，其中宝泉寺和郭家沟硫酸根离子贡献率分别为 和，铁离子贡献率分别为 和。关键词：水动力水质耦合模型；水质；加速技术；矿山区河网中图分类号：文献标志码：引用格式：李文钰，侯精明，刘占衍，等矿渣堆积对河网水质影响的数值模拟分析人民黄河，（）：，（，；，；，；，；，）：，：；随着人类社会的发展，对自然资源的需求不断提升，这一过程中矿产资源的开发利用造成大量矿渣堆积，威胁生态安全。如富含多种矿产资源的文峪河流域，因早年废弃矿渣未得到有效处理而致使其在河道堆积，上游来水或降水对其冲刷、淋溶和浸泡后产生了大量高浓度污染物，对下游河道水环境及居民饮用水安全造成了极大威胁。为遏止矿渣堆积对生态环境的不利影响，提出高效精准的水环境治理方案，第 卷第 期 人 民 黄 河 ，年 月 ，需要定量分析矿渣堆积后流域水系污染物的分布特征及迁移总量。近年 来 国 内 外 学 者 常 使 用 的 水 环 境 模 型 有、模型等，如：曾晓岚等使用 模型对滇池东岸花卉大棚种植区降雨径流水质、水量进行了模拟，并分析了氮磷的单位面积负荷率及输出特征；刘骞等应用 模型构建了岷江流域分布式水文和污染负荷模型，并计算了各项减排措施及气象驱动条件对流域水环境改善的贡献率；刘晨辉等采用 模型模拟了、氨氮等在长江的迁移和扩散情况。以上模型对特定区域具有较好的模拟效果，但 等一维模型无法展示污染物的对流扩散过程；等水文模型参数较多，建模困难，且根据其本底特征确定污染物输移路径后无法在模拟过程中依照水动力条件自动调整路径，适用于大尺度流域水文模拟，无法对中小尺度流域的水动力和水质变化过程进行高分辨率模拟；等水动力水质模型虽然可以进行高分辨率水环境模拟，但其模型重点在于水动力模块，且对于百万数量级以上网格，模型运行速度较慢。基于此，本文采用 加速技术构建了小文峪河流域水环境数值模型，并充分利用现有数据模拟分析小文峪河流域污染物空间分布及污染物来源，以期为小文峪河流域水环境治理提供参考。材料与方法 研究区概况文峪河位于湖北省十堰市竹山县西北部，上游主要支流有小文峪河和界岭河，下游汇入堵河。文峪河流域属于亚热带季风气候区，雨量充沛，热量充足。水量多寡取决于气候条件，呈季节性变化。小文峪河流域地表水主要污染源为矿洞和矿渣堆，矿渣堆沿自然坡面、河道沟谷及近沟坡面溜坡放置，降水时水流从矿渣堆表面和内部流过，形成矿渣堆淋溶水，严重影响小文峪河的水环境。小文峪河从上游到下游主要支流有竹园沟、清石沟、孙家沟、宝泉寺及郭家沟（见图），各支流矿渣堆面积分别为 、，矿渣方量分别为 、。本研究选取超标污染物硫酸根离子和铁离子为目标污染物，对研究区污染物时空分布进行高分辨率模拟并量化分析小文峪河污染物来源。流量及污染物数据本研究采用 年野外单场次监测数据，监测点图 研究区矿渣堆分布及监测点布设 个，见图。监测指标主要有硫酸根离子、铁离子、流量等。监测时间为旱季，流量较为稳定，因此以该场次监测数据为模型率定数据。监测点主要分布在矿渣堆上下游及支流与主河道交汇处，因此以支流矿渣堆上游监测水质和流量数据为模型边界条件，认为无矿渣堆积的支流上游来水污染物浓度为。支流上游来水对部分矿渣进行冲刷与矿洞水混合后经径流过程到达下游，发生复杂的物理变化，采用一级冲刷过程来综合概化地表水流经矿渣堆并与矿洞水混合到达下游的污染物浓度变化规律。一级支流在向主河道汇流过程中会有二级支流的汇入，由于缺乏二级支流监测数据，因此采用一级衰减过程来描述污染物浓度下降过程。小文峪河主河道上游流量和水质边界条件根据监测数据设置，即不考虑研究区域上游污染物空间分布及水量汇集特征。数值模型及求解方法模型的地表水动力模块控制方程为平面二维浅水方程（简称），主要针对具有自由液面且以平面运动为主的水流，只考虑水平方向流速，忽略垂向流动，且忽略运动黏性项、紊流黏性项、风应力和科氏力。二维浅水方程的守恒格式可用如下矢量形式来表示：（）（）（）人 民 黄 河 年第 期 （）（）式中：为时间；、分别为、方向的单宽流量；为水深；为变量矢量，包括水深及、方向的单宽流量；、分别为、方向的通量矢量；为重力加速度；为源项，包括降水源项、底坡源项及摩阻源项；、分别为、方向的流速；为入渗源项；为河床底面高程；为谢才系数，其中 为曼宁系数。水质模块控制方程为对流扩散方程：（）（）（）（）（）（）式中：为污染物垂向平均浓度；、分别为、方向的扩散系数，因对流较强，故忽略扩散作用；为点源排放的流量强度；为点源污染物垂向平均浓度；为污染物浓度对时间的全导数 ，即转化项，如污染物生化反应、生物作用等。仅考虑衰减与冲刷反应时，转化项 可表示为 （）式中：为衰减系数，为冲刷系数。采用 格式有限体积法求解二维浅水方程，选用二阶 方法进行变量重构，通过 近似黎曼求解器求解控制单元界面上的质量通量与动量通量。为提高计算效率，模型结合 架构实现 加速，同时保持精确度和高稳定性。污染物通量计算污染物通量是指水体中的污染物在单位时间内通过某一断面的总质量，是水文、地质、化学及生物等综合作用的结果。污染物通量计算公式为 （）式中：为污染物通量；为平均流量；为转化系数，值为 。模型构建地形数据的分辨率为 ，共有 个网格，模型计算采用开边界，仅在支流监测点位置设置入流流量和污染物初始浓度。研究区处于山区，且河道长期有水，可认为河道水流无下渗；参照文献，根据当地实际情况，取曼宁系数为。由于矿洞水仅有极少部分能汇入地表水，且流量较小，因此将其与各支流的发源地进行合并建模。研究区以矿渣堆为主的污染源释放的污染物会进入地表水和地下水，再随地表水和地下水向下游迁移，使下游地表水和地下水中污染物超过相应水质标准。由于矿渣堆积方量不同，因此采用 技术确定各支流矿渣堆积的范围，分别率定水流经过各支流矿渣堆及其上下游时污染物的冲刷系数和衰减系数，当地表水以一定流速通过矿渣堆时，模型将计算出各支流矿渣堆析出的污染物量。为分析矿渣堆中硫酸根离子和铁离子对流域水质的影响，采用基于 加速技术的水动力水质耦合模型对污染物的迁移转化过程进行模拟。为验证模型模拟结果的准确性，根据郭家沟、宝泉寺、孙家沟、清石沟及竹园沟的水质监测数据，对模型进行参数率定，结果见表。表 水质参数率定结果 位置硫酸根离子衰减系数冲刷系数铁离子衰减系数冲刷系数竹园沟清石沟孙家沟宝泉寺郭家沟 结果分析水质模拟结果见表，可以看出各监测点污染物模拟结果相对误差均在 以内，表明所构建模型能较好地模拟小文峪河流域硫酸根离子和铁离子在河网中的迁移转化过程。污染物空间分布特征研究区硫酸根离子和铁离子空间分布见图、图，可以看出矿渣堆上游污染物质量浓度均较低，经过矿渣堆途中污染物质量浓度逐渐上升，在矿渣堆下游污染物质量浓度达到峰值。竹园沟、清石沟、孙家沟、宝泉寺、郭家沟 条支流硫酸根离子质量浓度分别从、上升至、，铁离子质量浓度分别从、上升至、。矿渣堆被冲刷过程中，各支流污染物质量浓度总体变化趋势大致相同，在入流点人 民 黄 河 年第 期表 硫酸根离子和铁离子模拟结果和实测结果对比监测点编号硫酸根离子实测值（）模拟值（）相对误差 铁离子实测值（）模拟值（）相对误差 到矿渣堆上游污染物质量浓度呈下降趋势，在矿渣堆上游至矿渣堆下游污染物质量浓度呈上升趋势。其中竹园沟上游矿渣堆积范围较广，水流到达该支流与主河道交汇处时，污染物质量浓度达到最大，之后受上游水流汇入、污染物衰减过程及河网本底特征影响，污染物质量浓度呈减小趋势。清石沟矿渣堆面积较小，在水流经矿渣堆到达主河道之前，两种污染物质量浓度经稀释和底泥吸附等作用有所下降。支流孙家沟地表水的硫酸根离子和铁离子质量浓度均小于上游来水的，水流汇合后对主河道污染物具有一定稀释作用，使主河道污染物质量浓度下降；宝泉寺和郭家沟矿渣堆图 硫酸根离子质量浓度空间分布图 铁离子质量浓度空间分布直达下游与主河道交汇处，故在交汇前硫酸根离子和铁离子质量浓度均较大，且宝泉寺和郭家沟支流流量较大，因此水流进入主河道后，硫酸根离子和铁离子质量浓度均发生了明显变化，致使主河道下游污染物质量浓度增大。矿渣堆下游是河网污染物集中整治的关键区域。水量及污染物通量分析各支流对小文峪河干流水量贡献不一，由大到小依次为宝泉寺、郭家沟、竹园沟、清石沟和孙家沟，占比分别为、。整体水量差距相对不大，但由于各支流下游污染物质量浓度不同，因此各支流硫酸根离子和铁离子负荷对主河道的贡献率差距较大。各支流水流和日污染物通量见图。各支流输入主河道的污染物通量由大到小依次为宝泉寺、郭家沟、竹园沟、清石沟、孙家沟，其中硫酸根离子通量分别为、，铁离子通量分别为、。硫酸根离子和铁离子入河通量在不同支流变化趋势较为一致，最大值均出现在宝泉寺。受水量和矿渣堆积方量不同影响，污染物通量在各支流呈现不同变化，以硫酸根离子为例，其污染物入河通量的最大值出现在宝泉寺，最小值出现在孙家沟，最大值为最小值的 倍，而宝泉寺水量仅为孙家沟的倍，表明污染物通量的差异主要与矿渣堆积方量有关。图 污染物通量及水量 空间分布特征各支流硫酸根离子、铁离子入河通量大小顺序均为宝泉寺郭家沟竹园沟清石沟孙家沟，可见入河污染物主要来源于宝泉寺和郭家沟，二者占比合计在以上，其他支流的占比相对较小。硫酸根离子污染负荷主要来自郭家沟和宝泉寺，其对小文峪河主河道的贡献率分别为、；除郭家沟和宝泉寺外，竹园沟对河网中硫酸根离子的贡献率也较大，为人 民 黄 河 年第 期。在各支流中，宝泉寺对小文峪河主河道铁离子的贡 献 率 最 大，为 ，其 次 为 郭 家 沟，为。由各支流对小文峪河主河道污染物的贡献率可知，整体上，宝泉寺的污染负荷大于郭家沟的，其次为竹园沟、清石沟和孙家沟；流域内污染负荷主要来自下游支流宝泉寺和郭家沟，其产生的硫酸根离子、铁离子 污 染 负 荷 分 别 占 整 个 研 究 区 的 和。各支流污染物通量及贡献率不同的主要原因是，矿渣堆积方量不同，其中宝泉寺和郭家沟矿渣堆占地面积分别为、，其他 条支流矿渣堆占地面积之和为 ；宝泉寺和郭家沟流量较大，水流对矿渣堆的冲刷较为严重，导致污染物质量浓度增大，污染负荷增大。因此对小文峪河流域的水环境治理方案制定过程中，应首先考虑宝泉寺和郭家沟两条支流的整治。结论以十堰市小文峪河流域为例，基于二维水动力水质耦合模型，充分利用各支流矿渣堆属性、流域河网本底特征及地表水、矿洞水、淋溶水水质和水量监测数据等，构建了基于 加速技术的矿山区河网水环境高分辨率数值模型，对小文峪河流域硫酸根离子和铁离子的空间分布进行了可视化模拟，并量化分析了小文峪河污染负荷来源。模型模拟结果与污染物实测值相比，相对误差小于，表明模型可以对矿山区河网水环境变化过程进行模拟。小文峪河主河道的水量贡献率由大到小分别为宝泉寺、郭家沟、竹园沟、清石沟和孙家沟，占比分别为、。硫酸根离子和铁离子对小文峪河贡献率由大到小均为宝泉寺郭家沟竹园沟清石沟孙家沟，其中宝泉寺和郭家沟硫酸根离子贡献率分别为 和，铁离子贡献率分别为 和。综上，基于 加速技术的水动力水质耦合模型对矿山区河网水环境进行模拟效果较好，在小文峪河流域水环境治理方案制定过程中，应首先考虑宝泉寺和郭家沟两条支流的整治。参考文献：张浙，卢然，伍思扬，等长江经济带矿山土壤重金属污染及健康风险评价 环境科学，（）：张龙，宋波，黄凤艳，等湖南锡矿山周边土壤农作物系统锑迁移转换特征及污染评价环境科学，（）：黄锦勇，覃铭，马文英，等受酸性矿山废水影响河流悬浮物中重金属污染特征分析与生态风险评价环境化学，（）：崔飞剑，覃光雄，曾海龙，等沙河流域重污染支流表层沉积物中氮、磷和重金属的空间分布特征及污染评价环境工程，（）：王辉，赵悦铭，刘春跃，等辽河干流沉积物重金属污染特征及 潜 在 生 态 风 险 评 价 环 境 工 程，（）：，郭亚科，高燕燕，钱会，等楮河流域水体重金属时空分布特征及健康风险评价环境工程，（）：康得军，吕茳芏，脱向银，等煤矸石山覆土体系重金属污染特征及生态风险评价环境工程，（）：李桂伊，侯精明，高波，等矿渣堆积对小流域行洪能力影响的数值模拟研究水资源与水工程学报，（）：曾晓岚，王涛涛，罗万申，等设施农业生产区降雨径流和氮磷输出特征及模拟：以滇池东岸花卉大棚种植区为例湖泊科学，（）：刘骞，王维，罗彬，等 模型下“十三五”中期污染减排措施及气象条件对岷江流域水环境改善贡献环境工程，（）：刘晨辉，刘思飔，李丹，等基于 模型的长江中环排污口水质影响分析中国农村水利水电，（）：，：侯精明，李桂伊，李国栋，等高效高精度水动力模型在洪水演进中的应用研究水力发电学报，（）：栾广学，侯精明，杨露，等耦合高效高精度水动力模型的多组分污染物输移及衰减反应模型研究水资源保护，（）：，：（下转第 页）人 民 黄 河 年第 期：黄越，程静，王鹏中国北方农牧交错区生态脆弱性时空演变格局与驱动因素：以盐池县为例干旱区地理，（）：，：，：，：，：郭泽呈，魏伟，庞素菲，等基于 和遥感指数的干旱内陆河流域生态脆弱性时空演变及动因分析：以石羊河流域为例生态学报，（）：白艳芬，马海州，张宝成，等基于遥感和 技术的青海湖环湖地区生态环境脆弱性评价遥感技术与应用，（）：布和敖斯尔，马建文，韩秀珍，等中国陆地生态系统脆弱带遥感模型遥感学报，（）：安芬，李旭东，程东亚贵州省乌江流域生态脆弱性评价及其空间变化特征水土保持通报，（）：邱彭华，徐颂军，谢跟踪，等基于景观格局和生态敏感性的海南西部地区生态脆弱性分析生态学报，（）：杨雯娜，周亮，孙东琪基于分区集成的黄河流域生态脆弱性评价自然资源遥感，（）：徐辉，丁祖栋，武玲玲，等黄河下游沿黄城市生态系统健康评价人民黄河，（）：，李宝亭，吴卿，杨硕果，等水土保持弹性景观功能与生态脆弱性研究人民黄河，（）：，黄钰涵，闫浩文，李小军，等乌鲁木齐市遥感生态距离指数变化监测与评价遥感信息，（）：张娟，文广超，王恩营，等基于遥感生态指数的焦作市生态环境动态监测与评价 水土保持通报，（）：徐涵秋区域生态环境变化的遥感评价指数中国环境科学，（）：，（）：曾永年，向南平，冯兆东，等 特征空间及沙漠化遥感监测指数研究地理科学，（）：王苗，刘普幸，乔雪梅，等基于 的宁夏生态环境质量时空演变及其驱动力贡献率分析生态学杂志，（）：姚昆，周兵，李小菊，等基于 熵权模型的大渡河流域中上游地区生态环境脆弱性评价水土保持研究，（）：黄方，刘湘南，张养贞 支持下的吉林省西部生态环境脆弱态势评价研究地理科学，（）：周玲美，王世航，权玲基于遥感和地理信息系统的内蒙古呼日查干淖尔地区生态脆弱性评价生态与农村环境学报，（）：魏伟，石培基，周俊菊，等基于 和组合赋权法的石羊河流域生态环境质量评价干旱区资源与环境，（）：蔡元刚，王明田，蔡怡亨，等利用地面气象资料建立四川省日总辐射计算模型中国农业气象，（）：张笑楠，王克林，张伟，等桂西北喀斯特区域生态环境脆弱性生态学报，（）：时卉，杨兆萍，韩芳，等自然遗产地生态脆弱性分析与评价：以托木尔区域为例干旱区地理，（）：【责任编辑 吕艳梅】（上接第 页）龚佳辉，侯精明，薛阳，等城市雨洪过程模拟 加速计算 效 率 研 究 环 境 工 程，（）：，（）：田泽斌，王丽婧，李小宝，等洞庭湖出入湖污染物通量特征环境科学研究，（）：朱昕阳，蒋彩萍，马晓雁，等太湖流域浙江片区出入境水量、水质及污染物通量湖泊科学，（）：程国微，杜展鹏，严长安，等水质监测频率与极端气候对高原湖泊入湖河流氮磷通量估算的影响环境科学学报，（）：唐洪武，闫静，肖洋，等含植物河道曼宁阻力系数的研究水利学报，（）：玛哈沙提哈孜哈力，努尔夏西曼斯尔天然河道的糙率确定方法分析能源与节能，（）：【责任编辑 吕艳梅】人 民 黄 河 年第 期