基于缓蓄快放系统的西安市护城河水质改善方法研究.pdf

第 卷 第 期 人民珠江 年 月 ：基金项目：陕西省水利科技计划项目（）；陕西省自然科学基金项目（）；西安市科技计划资助项目（）收稿日期：作者简介：孙超（），男，硕士，高级工程师，主要从事水生态修复和水质综合治理方向研究。：孙超，韩伟刚，农晓英，等 基于缓蓄快放系统的西安市护城河水质改善方法研究 人民珠江，（）：基于缓蓄快放系统的西安市护城河水质改善方法研究孙超，韩伟刚，农晓英，王慧，王文玉，高泽海（陕西省水利电力勘测设计研究院，陕西西安 ；西安理工大学省部共建西北旱区生态水利国家重点实验室，陕西西安 ）摘要：研究缓蓄快放系统，探索水资源短缺地区的河流水动力提升方法及其对水质影响规律。以西安市护城河南西线为研究区，根据护城河现状水质与水动力数据，利用 构建了西安市护城河水动力 水质耦合模型。依护城河现有坝体设计蓄水闸门，模拟河道上游快速泄水对下游河道的冲淤能力，探究不同蓄水高度和不同水质条件快速泄水对护城河水动力和水质的影响规律。基于 的数值模拟结果表明，不同蓄水高度对水体的水动力特性影响显著，蓄水水质对污染物浓度具有显著影响。基于 的西安市护城河水动力 水质耦合模型模拟结果表明采用 蓄水高度与类水补给相结合的工况可以最大程度实现护城河水体的水质、水动力改善。西安市护城河采用符合地表水类的补水水源，通过提升护城河水动力的方式可以有效改善西安市护城河水质状态。关键词：水动力提升；水质改善；数值模拟；西安护城河中图分类号：文献标识码：文章编号：（）牞 牞 牞 牞 牞 牗 牞 牞 牷 牞 牞 牞 牘 牶 牞 牞 牞 牶 牷 牷 牷 牷 人民珠江 年第 期随着社会经济的迅速发展和城市化建设的不断加快，水体污染和水资源短缺问题愈发突出，已成为制约社会可持续发展和人民生活水平提高的重要因素。城市景观水体作为现代城市建设的基础设施之一，起着防洪调蓄、维持区域生物多样性、气候微调和景观文娱等功能，是城市经济与生态环境平衡发展的重要标志 。因此，探索城市河流治理方法、保障城市河流水环境是城市水系高效管理重要研究方向之一。西安市护城河作为一条人工开挖典型的城市河流，是古代城垣建筑体系的一个重要组成部分。在西安市城市化建设过程中，护城河受到过不同程度的污染，影响到了城市生态环境和居民健康，同时也制约了城市经济和社会的发展 。西安市护城河极小的坡降和河道闸坝工程使河流流速极低，死水区与缓流区频繁出现，水体置换时间较长，污染物迁移扩散缓慢，水体复氧能力降低，水体自净能力减弱，造成水质恶化、水体富营养化的问题 。由于中国西北干旱半干旱地区水资源短缺、时空差异显著，能够用于提升护城河水体的水资源极为有限，制约着西安市护城河的水环境发展 。针对西安市护城河水环境存在的问题，利用有限的水资源提升水动力，是改善护城河水质的有效途径之一 。水环境模型是分析水环境问题的有效技术手段，其中 软件和 （）是 模 拟 河 流 水 质 的 代 表 性 软件 。目前的研究主要集中在水质指标空间分布和典型控制断面的指标变化特征方面 。水动力模型是通过数值模型模拟分析复杂流场的演进过程，针对不同的需求对研究区进行不同维度的模拟，其中 与 软件，被国内外众多学者用于模拟水流在复杂边界或实际地形上的流动过程 。鉴于 既能模拟水动力又能模拟水质，可用于构建水动力 水质耦合模型，且模拟结果可视化效果清晰，被广泛用于上游调水对下游库区水质影响、水库联合调度水质状态问题，分析污染物运移扩散规律、水环境容量计算、水质改善方案优化和水生态修复方案优选等方面 ，特别是在中小型城市河湖低流速小流量补水模拟方面有着显著的优势 。本研 究 以 西 安 市 护 城 河 为 研 究 区，利 用 构建西安市护城河水动力 水质耦合模型，并通过护城河现状水质、水动力监测结果验证所建模型的准确性。设计基于溃坝理论的“缓蓄快放”系统，即在河道某断面处修建闸门，通过闸门将上游小流量缓慢输入的来水积蓄起来，达到一定水量和水头时全开闸门，快速泄放，增大水流流速，置换缓流水体，冲刷底泥、污染物，达到保持水质的效果。通过缓慢蓄水快速泄流的方式解决水资源短缺条件下的水动力提升能力。利用所构建的水动力 水质耦合模型，研究护城河不同来水水质和不同蓄水高度条件下“缓蓄快放”系统对水体流速、水质浓度分布及水体交换特性的影响，以期为水资源短缺地区的城市景观水体水质改善提供参考依据。研究区概况西安市（，）是中国西北地区最大的城市，坐落于秦岭以南的关中平原上，属于温带半湿润大陆性季风气候。西安市夏季炎热多雨，冬季寒冷干燥，年平均降雨量 。西安市护城河全长 ，总库容 万 。西安市护城河补水水源分别来自大峪水库和北石桥再生水厂。西安护城河以东南城角作为分水岭分为东北线与南西线，河水分别途经东北线与南西线自西北角退水系统进入汉城湖截流涵，再经过漕运明渠流入渭河（图 ）。西安护城是古城堡防御体系的一个重要组成部分，为了保持其高墙深壕的历史原貌，护城河全线综合改造规划建设 座拦河坝以壅高水位，护城河改造段各库区上下游最大河道高程差仅为 ，最大比降低于 。西安市护城河全线综合改造建设的拦河坝在壅高水位的同时，造成水体断面面积较大，致使现状供水条件难以维持护城河水体的自然流态，形成缓流区和死水区。较差的水动力状况致使水体置换时间较长、污染物的迁移扩散缓慢，将加剧水质的恶化。人民珠江 年第 期图 西安市护城河地理信息 水动力 水质耦合模型构建 水动力数学模型方程 水动力模型是基于不可压缩的沿水深积分的流体雷诺平均 方程，该模型包括连续性、动量和密度方程，且服从布辛涅斯克（）和静水压力假设 。平面二维水流连续方程见式（）：（）平面二维水流的动量方程见式（）、（）：()()()（）()()()（）且有，。式中、笛卡尔坐标系；时间；河底高程；静水深，即基准面到底床的距离；总水深，；、方向的流速分量；、基于水深的平均流速；重力加速度；水的密度；水的参考密度；大气压力；、各方 向 的 应 力 张 量；、辐射应力张量各方向的分量；源汇项产生的流量，由电源进入周围水体的流速表示为（，）；侧向应力项，包括黏滞摩擦、湍流摩擦、差异平流，据沿水深平均流速的梯度计算，分为流体黏性应力、紊流应力 和水平对流 ，表达式如下。（）人民珠江 年第 期 ()（）（）水质数学模型方程 对流扩散模型假定污染物垂向混合均匀，污染物的浓度只在水平面上发生变化 。因此污染物的输移扩散可采用如下二维对流扩散方程。()()()()()（）()（）式中、方向的扩散系数；污染物浓度；线性衰减系数；水深；、方向上的流速；源汇项流量；源汇项污染物质的浓度。水动力 水质耦合模型的数值求解采用基于单元中心的有限体积法对二维浅水控制方程进行空间离散求解 。该方法将空间划分为不重复的控制单元，能够保证水动力模型中水量和动量守恒，模型在时间上的离散方法为显性欧拉法。西安市护城河水动力 水质耦合模型 网络剖分西安市护城河南西线，从护城河东南角至西北角共建有 座堤坝，各堤坝的距离及主要参数见表。本研究设计“缓蓄快放”系统，模拟在 座堤坝上修建不同高度的闸门，以增加水头，并通过快速泄放提升水动力，置换缓流水体，实现水质提升。构建水动力 水质耦合模型时，假定每一段坝的坝前常水位与该段坝顶高程齐平，各段水质以该段采样点监测的水质数据平均值表征。表 护城河南西线堤坝参数单位：堤坝坝顶高程河底高程坝宽与进水口距离号坝 号坝 号坝 号坝 根据收集到的护城河地形数据图，利用 前处理工具网格生成器，将长约为 的护城河南西线划分为 个三角形网格，包含 个节点，导入高程数据文件，并基于剖分后的网格进行插值，得到研究区的地形见图 ，局部网格剖分示意见图 。图 研究区域地形高程图 局部网格的剖分示意 模型边界设置将位于东南城角的补水点的流量定义为上边界，位于护城河西北角的退水口定义为下边界。水位边界为边界条件，其余岸线则定义为陆地边界。研究不考虑降雨、蒸发及风场的作用，在水动力模型的构建中主要涉及干湿水深、涡黏系数、床底糙率、源汇项以及水工建筑物等因素。模型参数校准 模拟周期与初始条件模拟时间从 年 月 日至 月 日，总时间步数为 ，主时间步长为 。本文所建水动力 水质耦合模型，模拟初始水位与流速设置为护城河现场实际测点多次监测的水位与流速数据平均值。水质对流扩散模拟的初始条件为研究区的初始污染物浓度分布情况，首次运行模型时选取历次取样点各污染物指标的平均值作为模型的初始条人民珠江 年第 期件。初始值为 年 月 日实测结果平均值。为了降低模拟结果对初始条件的依赖，提高模拟结果的合理性与准确性，本次模拟采用之前模拟结果为初始条件的热启动方法。参数设置对于干湿判别技术的设计采用模型的预设值，即干水深为 ，淹没水深为 ，湿水深为 。模拟设定 系数为 ，曼宁系数为 。模拟选取总氮、总磷、氨氮及高锰酸钾指数作为对流扩散模拟的污染物指标 。对流扩散模块中，研究区扩散系数和降解系数的率定结果表明，扩散系数的取值在某一范围内时，率定精度基本相同，超过这一范围，模拟值与实测值的均方根误差显著增加。降解系数对率定结果较为敏感，降解系数的微小变动都会造成模拟值的波动。在查阅相关文献及对模型多次率定的基础上，最终确定各污染物的扩散系数为 ，总氮、总磷、氨氮、高锰酸钾的降解系数分别为 、。模拟结果分析选取模拟期间建国门及朱雀门的实测水位及流速数据作为验证点，测点水位和流速的模拟值与实测值见图 、。由图 、可知，模拟结果的水位及流速与实测值相吻合。为定量分析模型准确性，采用均方根误差（，）、平均绝对误差（，）和决定系数作为模型评价指标 ，结果见表 。）建国门）朱雀门图 水位模拟值与实测值对比）建国门）朱雀门图 流速模拟值与实测值对比表 水位、流速均方根误差及平均绝对误差位置 水位流速 水位流速水位流速建国门 朱雀门 人民珠江 年第 期由表可知，模拟误差较小，所构建的护城河二维水动力模型能很好地反映研究区的水动力特性，可用于后续的水动力提升模拟。由护城河长期实测水质指标结果可知，总氮、总磷、氨氮及高锰酸钾指数是西安市护城河水质主要超标指标，因此，本文为通过分析这 个指标改善效果来验证所提“缓蓄快放”系统的有效性。图 、分别为建国门和朱雀门总氮、总磷、氨氮及高锰酸钾指数模拟数值与实测数据的对比结果。）总氮）总磷）氨氮）高锰酸钾指数图 建国门水质指标模拟值与实测值对比）总氮）总磷）氨氮图 朱雀门水质指标模拟值与实测值对比人民珠江 年第 期）高锰酸钾指数续图 朱雀门水质指标模拟值与实测值对比个水质采样点各水质指标模拟值与实测值的 、和 见表 。表 采样点水质指标均统计结果指标 建国门 朱雀门 建国门 朱雀门建国门 朱雀门总氮 总磷 氨氮 高锰酸钾指数 由模拟分析结果可知，本研究所构建模型能够良好地模拟水体物质的对流扩散，具有较好的模型精度，可用于后续的水质变化过程模拟。模型应用 基于水动力的水质改善设计西安市护城河现状补水水源分别为大峪水库来水和西安市北石桥污水处理厂再生水。大峪水库执行 地表水环境质量标准 类标准。西安市北石桥污水处理厂于 年 月前，执行 城镇污水处理厂污染物排放标准 一级 标准的工艺要求。年 月执行 陕西省黄河流域污水综合排放标准 标准的工艺要求。年 月至今，主控指标执行 地表水环境质量标准 类标准对护城河进行补水。因此，根据西安市护城河补水水质历史数据，设计不同的补水水质和坝前水头，研究不同的工况下，通过“缓蓄快放”系统快速泄流对护城河水质与水动力的影响规律，工况设计见表 。表 不同模拟工况的设置工况坝前水头 补水水质类别类 类 类 类类 类 类 类工况 为现状蓄水高度，未执行 地表水环境质量标准 类标准时的补水情况。工况 为现状蓄水高度，执行 地表水环境质量标准 类标准时的补水情况。本文考虑水资源的稀缺性和水质改善的有效性，旨在利用现有资源进行补水的条件下实现西安市护城河水质的改善。在同样补水情况下，对于不同的蓄水高度，所需蓄水时间不同。因此，本文以较短的蓄水周期和良好的水质提升能力作为不同工况下“缓蓄快放”系统效果评价标准。水质水动力改善效果分析 流速特征比较不同模拟工况下，快速泄水后，待上下游水位相等，水流流速基本稳定时刻，各取样点的流速值见表。显然，不同蓄水高度对各点的流速值影响显著，不同蓄水高度下各点的流速值存在差异。由表 可知，蓄水 、蓄水 及蓄水 的工况对流速均有所改善，随着蓄水高度的增加，下游各点的流速值均有所增加。蓄水 时和平门处流速的增长率最大，增长率为 。蓄水 时和平门处流速的增长率最大为 ，蓄水 时朱雀门流速增长最大，为蓄水 时的 倍。不同蓄水工况下研究区的流速分布见图。对比各蓄水工况下的流速分布可知，随着蓄水高度的增加，坝后流速场的分布有明显差异。人民珠江 年第 期表 不同工况下各取样点流速值单位：工况东南角建国门和平门文昌门朱雀门含光门西北角 ）蓄水 ）蓄水 ）蓄水 图 不同蓄水工况下的流速分布蓄水 时，号坝处流速为 ，为现状工况下流速的 倍。此时 号坝下游 、处的流速分别为 、，对应位置较现状工况下的流速分别增加 、。蓄水时，号坝处流速为 ，为现状工况下流速的 倍。此时 号坝下游约 、处的流速分别为 、，对应位置较现状工况下的流速分别增加 、。蓄水 时，坝后流速值大于 的范围增加，大量水体的突然下泄导致局部流速较大，虽然对局部冲刷明显，但对于沿程其他断面流速的改善程度与蓄水工况差异不大。从蓄泄周期以及流速改善程度两方面综合分析认为蓄水 ，水动力改善效果最佳。水质特征比较不同模拟工况下，快速泄水后，待上下游水位相等，水流流速基本稳定时刻，各取样点在模拟时段内的总氮平均浓度见表 。从平均浓度而言，采用类水补水条件下，“缓蓄快放”系统对水体总氮浓度的改善效果微弱，随着蓄水高度的增加，总氮浓度呈现微弱减小趋势。图 为类水补水条件下，不同蓄水高度对文昌门和朱雀门附近水体总氮浓度变化曲线。由表 可知，相比于采用类水补水，利用类水补水时，总氮浓度显著降低。同样采用类水补水时，总氮浓度随蓄水高度增高而降低。“缓蓄快放”系统对下游水质的改善，随蓄水高度增高而增强。从总氮含量下降速率而言，蓄水高度由 人民珠江 年第 期变化为 时，总氮含量下降效果最为显著。蓄水高度由 增加到 时，通过快速泄流尽管能进一步降低总氮含量，但是相比于蓄水高度由 变化为 时，对总氮含量降低比率减弱。图 为类水补水条件下，不同蓄水高度对文昌门和朱雀门附近水体总氮浓度变化曲线。表 不同工况下各取样点的总氮平均浓度单位：工况东南角 建国门 和平门 文昌门 朱雀门 含光门 ）文昌门）朱雀门图 工况 总氮浓度）文昌门）朱雀门图 工况 总氮浓度由图 、可知，“缓蓄快放”系统蓄水高度越高，水体置换耗时越短、置换率越高。蓄水高度从 增高为 时，水质改善效果最为明显，蓄水高度从 增高为 时，通过蓄泄水体来改善水质效果与蓄水高度 差异不大。针对水资源短缺的西北地区，西安市护城河采用蓄水 的方案具有较短的蓄泄周期，水龄较小，水体交换能力强，水质提升效果显著。综合水动力和水质实验结果可知，通过将上游水位蓄高，快速泄水的方式是提高护城河水动力的有效途径。考虑蓄泄周期和流速改善程度，上游水位蓄水 对水动力改善综合效果最佳。西安市护城河通过使用地表类水作为补水水源时，综合考虑蓄水成本和水质提升效果，明确在蓄水高度 时对下游水体水质改善效果最佳。因此，西安市护城河采用“缓蓄快放”系统通过将符合地表类水的护城河补水水体蓄高 ，快速泄流能够有效提升护城河水动力，改善护城河水质。结论本文设计缓蓄快放系统，通过缓慢蓄水快速泄人民珠江 年第 期放的方法提升西安市护城河水动力、改善护城河水质。构建了基于 护城河水动力 水质耦合模型，并通过实测水位、流速及水质数据对所建耦合模型进行率定及验证。通过设置不同的蓄水高度及上游补水水质，模拟了 种方案下的水质、水动力改善效果，探究了缓蓄快放系统不同蓄水高度和不同补水水质共同作用对水体流速和水质关键指标浓度分布影响特性。结果表明：不同蓄水高度对流速影响显著，随着蓄水高度的增加，流速增加明显，不同补水水质对下游水体的水动力无显著影响；补水水质对护城河水质影响显著，随着蓄水高度的增加，“缓蓄快放”系统对水质改善程度先增加后减弱；西安市护城河采用地表类水作为补水水源、蓄水高度 工况下，通过“缓蓄快放”系统对护城河水质改善最优。参考文献：陈军，吴程，李爽，等 灞河流域水环境数值模拟研究 西北农林科技大学学报（自然科学版），（）：谢启姣，刘进华 年武汉城市湖泊时空演变及其生态服务价值响应 生态学报，（）：牞 牞 牞 犤 犦 牞 牞 牶 王浩，孟现勇，林晨 黑臭河道治理关键技术及其应用 人民长江，（）：牞 牞 牞 牞 犤 犦 牞 牞 牞 牶 黄子奇，张建丰，李涛，等 库区几何要素对溃坝洪水特性影响的模拟研究 西安理工大学学报，（）：胡晖，张建丰，刘武刚，等 极端流态溃坝过程中水体流速与流量变化的模拟研究 环境污染与防治，（）：杨晨，徐明德，郭媛 基于 的汾河水库突发环境事件数值模拟 灌溉排水学报，（）：牞 牞 犤 犦 牞 牞 牶 牞 牞 牞 犤 犦 牞 牞 牶 牞 牞 牞 犤 犦 牞 牞 牗 牘 牶 牞 牞 牞 牶 犤 犦 牞 牞 牶 薛树红，李明，吴巍，等 基于多耦合模型的盛桥河流域水环境模拟研究 西北农林科技大学学报（自然科学版），（）：牞 牞 牞 犤 犦 牞 牞 牶 牞 牞 牞 犤 犦 牞 牞 牶 牞 牞 牞 犤 犦 牞 牞 牶 尚福强，蔡一，唐莉华，等 郓城南湖公园水体环境改善与提升方案研究 水利水电技术（中英文），（）：蔡一，尚福强，唐莉华，等 基于数值模拟的南湖流场优化与水质改善研究 水资源保护，（）：闫军波，冯健，张守平 基于 的小南海湖补水方案效果分析 人民珠江，（）：周瞡，张俊，徐凌云 城市森林公园水系流场模拟及水环境综合整治：以上海共青森林公园为例 人民珠江，（）：牞 牞 牞 犤 犦 牞 牞 牶 牞 牞 犤 犦 牞 牞 牗 牘 牶 牞 牞 牞 牶 牞 牞 犤 犦 牞 牞 牶 夏妍，窦鹏，楼春华，等 不同植被覆盖情景下景观水体水质改善效果模拟研究 环境科学学报，（）：人民珠江 年第 期 牞 牞 牶 牞 牞 犤 犦 牞 牞牗 牘 牶 牞 牞 牞 犤 犦 牞 牞 牶 牞 牞 犤 犦 牞 牞 牗 牘 牶 杨程，吕岁菊，李春光，等 基于非结构网格的三维水流数值模拟研究 西北农林科技大学学报（自然科学版），（）：陈虹，李家科，李亚娇，等 暴雨洪水管理模型 的研究及应用进展 西北农林科技大学学报（自然科学版），（）：牞 牞 牞 犤 犦 牞 牞 牞 牶 牞 牞 犤 犦 牞 牞 牶 牞 牞 牞 犤 犦 牞 牞 牶 牞 牞 牞 犤 犦 牞 牞 牗 牘 牶 （责任编辑：高天扬）（上接第 页）王艺霖，苏程佳，陈晓宏，等 磨刀门水道枯水期咸潮变化趋势及成因 水资源保护，（）：牞 牞 牞 犤 犦 牞 牞 牶 牞 牞 犤 犦 牞 牞 牶 犤 犦 牞 牞 牶 牞 牞 牞 牞 牞 牞 牞 犤 犦 牞 牞 牶 牞 牞 牞 犤 犦 牞 牞 牶 牶 犤 犦 牞 牞 牗 牘 牶 （责任编辑：程茜）