1、连云港市水利局 颜建1 刘二军1 刘赟 1 宋彦生1 江苏连云港 222003
2、中国水利水电科学研究院 李岩桃2 张迪2 北京 100038
摘要：随着社会发展，城市水环境是城市建设和可持续发展的重要保障。本文以连云港市为例，通过构建一维水动力-水质耦合模型，对连云港市区河网联合调控典型情景进行了模拟，分析了市区主要污染物氨氮在不同初始浓度和引水流量条件下，在河道中输移扩散规律，由此提出了面向水环境改善的连云港市区河网联合调控优化方案建议。为连云港市实现“高质发展、后发先至”提供水环境保障和技术支持。
关键词：连云港市；河网联合调控；水动力-水质耦合模型；污染物输移扩散
1引言
城市水环境是城市建设和可持续发展的重要保障。如何改善城市水环境、如何应用新一代信息技术实现城市河道的长效久治和智慧化管理，是自然科学领域的热点问题和前沿问题，也是城市发展的关键科学问题之一。
本文以连云港市为例，通过摸底调查，梳理连云港市区河网水环境现状，识别主要污染源和水环境胁迫因素，系统构建一维水动力-水质耦合模型，综合模拟和预测不同情景、不同调控方案下的污染物输移扩散过程，提出面向水环境改善的河网多闸联动调控方案，探索长效久治的智慧化河道管理模式。为城市河网的联合调度、智慧管理、科学决策和快速响应奠定坚实基础，为连云港市实现“高质发展、后发先至”提供水环境保障和技术支持。
2研究内容
2.1研究区域
连云港市地处江苏省北端，位于北纬33°58′55″～35°08′30″，东经118°24′03″～119°54′51″之间。连云港市水系基本属于淮河流域沂沭泗水系，位于淮河流域沂沭泗水系的最下游，沂沭地区的主要排洪河道新沂河、新沭河等均从市内入海，处在江苏省供水网络的末梢，既是水资源较为紧缺的地区，又是著名的洪水走廊[1]。研究区域涉及玉带河、大浦河、龙尾河、东盐河和排淡河5条河流，电厂闸、玉带河闸等9座闸门，2处省考断面及3处国考断面。本文以连云港市区4条典型调水路线做为研究目标，进行河网联合调控研究，各调水路线如图1所示。
.png)
2.2主要污染源和水环境胁迫因素
2016年至今，连云港市相继实施了“清水进城”工程建设项目、黑臭水体整治项目，开展了生态调水工作。目前，市区河道已初步实现了从“清水进城”到“碧水绕城”的蜕变，连云港市区水环境问题得到了有效改善，但市区河道仍存在突发污染的风险。相关资料显示，连云港市区河道主要污染物为氨氮，主要污染源为生活污水，2017～2018年国考断面监测数据显示，个别月份氨氮浓度超出Ⅳ类，甚至Ⅴ类水标准，水环境容量为负[1][2]。本文以氨氮做为连云港市区水环境胁迫因素，构建连云港市区河道一维水动力-水质耦合模型，模拟不同情景下连云港市区河道内氨氮的输移扩散规律，为连云港市实现“高质发展、后发先至”提供水环境保障和技术支持。
3模型建立
3.1模型控制方程
3.1.1一维水动力模型控制方程
采用一维圣维南方程描述河流水动力特征，由反映质量守恒定律的连续性方程和反映动量守恒定律的动量方程共同刻画。根据流域上下游拓扑关系，将模型计算区域分成不同单元和节点进行分区计算，控制方程如下：
.png)
式中，Q为流量；A河道过水断面面积；Z为水位； u为流速； cn为河段的谢才系数；R水力半径； n为河床糙率系数；q为区间入流。
3.1.2一维水质模型控制方程
描述河流水系污染物输移运动和浓度变化规律的控制方程为对流扩散方程，其方程形式如下：
.png)
式中，D为弥散系数；C为河段中某种污染物的浓度；t为时间；x为河水流动的距离；u为河段水流的平均流速。
河流的弥散系数是分子扩散系数、紊动扩散系数与纵向离散系数之和，但纵向离散系数一般要比分子扩散系数、紊动扩散系数大得多，当河流的弥散系数很小时，一般可以忽略，即D=0，入式（4）。对于均匀河段，流量和排污稳定时，各断面的污染浓度不随时间变化，如式（5）。
3.2河网概化
对研究区域河网进行概化，玉带河4.2km、大浦河10.2km、龙尾河5.1km、东盐河27.2km和排淡河17.1km河道断面根据设计资料确定；研究区域内涉及的闸门、泵站按照模型中可控建筑物进行概化，其中玉带河闸与玉龙泵站为分流节点，大浦闸、新城闸和大板跳闸为水位控制节点。
3.3边界与参数
边界：模型上游边界以流量做为控制条件，位于电厂闸处；模型下游边界以水位进行控制，由于缺乏基础资料，水位根据明渠均匀流公式进行推求，位于大浦闸、新城闸和大板跳闸处。
糙率：根据经验，治理后混凝土河道糙率取0.02[3]。
分流比：玉带河闸处分流比为1:2，上游来水的1/3流向大浦河，2/3流向东盐河；玉龙泵站处分流比为1:2，上游来水的1/3流向龙尾河，同时满足龙尾河过流量≤8m3/s。
氨氮目标浓度：以研究区域考核断面水质标准作为目标，氨氮浓度设定为2.0mg/L（Ⅴ类水）。
3.3情景工况设计
本文连云港市区河道水质指标中最容易超标的氨氮为例进行模拟，由于调水路线较短，氨氮在河道中的降解有限，因此在模拟过程中不考虑其降解。构建的一维水动力-水质耦合模型中氨氮的运动包括输移和扩散两种，其中扩散主要发生在两个不同浓度水体的交界面，从浓度低向浓度高方向扩散，扩散速度有限，同时本次模拟涉及的距离较短，在模拟过程中氨氮的扩散作用相较输移作用基本可以忽略不计[3][4][5]。
在此基础上，本文结合现场调查结果及相关资料，设计了4个电厂闸引水流量梯度，分别为10m3/s、20m3/s、30m3/s和40m3/s，引水水体的水质按V类设定，氨氮浓度为2.0mg/L；河段内氨氮浓度设计3个梯度，分别在Ⅴ类水的基础上，分布按照超标10%、30%和50%设定，共组合设计12种典型情景如表1所示。
.png)
4模拟结果分析
将4条调水路线与12种典型情景进行组合，形成48种调水方案进行模拟，对模拟结果进行统计分析，模拟结果如表2所示，发现48种调水方案整体上存在相同趋势。以调水路线1猴嘴闸上游500m处省控断面为例，在相同的流量下，随着河段内氨氮初始浓度的增加，同一考核断面氨氮浓度达标的时间基本一致（图2）；在氨氮浓度相同下，随着引水流量的增加，同一考核断面水质达标的时间有所减少，但减少的幅度随着流量的增加而呈现下降趋势（图3）。
.png)
.png)
.png)
小时
图3 调水路线1同浓度下猴嘴闸上游500m省控断面模拟结果
5河网联合调控建议
对48种调水方案模拟结果进行统计，分析各考核断面水质达标时间、氨氮浓度和电厂闸引水流量关系。结果显示，以水质达标所需的水量最少为考量时，推荐在连云港市区河道省控或者国控断面水质超标时选择情景1进行调水，此情景下电厂闸引水流量为10m3/s，调水路线需根据考核断面位置选择；综合考量水质达标所需时间和水量时，推荐在连云港市区河道省控或者国控断面水质超标时选择情景8进行调水，此情景下电厂闸引水流量为30m3/s，调水路线需根据考核断面位置选择。
图4以猴嘴闸上游500m处省控断面为例，展示了不同调水路线下，实现水质达标所需的时间，结果显示，调水路线1、2水质达标所需的时间基本相同，调水路线4在引水流量较低时，水质达标所需时间较其余两条调水路线长，引水流量较高时，水质达标所需的时间和其余两条调水路线接近。因此，在实际调水过程中，可根据此规律合理选择调水路线。
.png)
6结论
本文以连云港市为例，通过调查分析识别其主要污染源和水环境胁迫因素，构建了一维水动力-水质耦合模型，模拟了48种河网联合调控典型情景下，市区主要污染物氨氮在河道中输移扩散规律，提出了面向水环境改善的连云港市区河网联合调控建议，主要结论如下：
1、相同引水流量下，随河段内氨氮初始浓度的增加，同一考核断面氨氮浓度达标的时间基本一致；
2、相同氨氮初始浓度下，随着引水流量的增加，同一考核断面水质达标的时间有所减少，但减少的幅度随着流量的增加而呈现下降趋势；
3、基于氨氮在河道中输移扩散规律，建议连云港市在水质调控过程中，以水质达标所需的水量最少为考量时，推荐选择10m3/s的小流量方案调水；综合考量水质达标所需时间和水量时，推荐选择30 m3/s的中等流量方案调水，调水路线均需根据考核断面位置选择。
参考文献
[1]潘志富, 刘二军, 孙佑祥, 潘琼, 张勇军.连云港市区水环境调查与治理规划[J]. 水科学与工程技术, 2019(4): 1-3.
[2]张言, 宗国营, 王玉甫, 占飞. 连云港市通榆河水环境综合治理浅析[J]. 西北水电, 2021(1): 22-27.
[3]刘洪燕. 北方小基流河流氨氮综合降解系数特征研究[J]. 能源与节能, 2016(9): 109-128.
[4]王蓉, 黄天寅, 吴玮. 典型城市河道氮、磷自净能力影响因素[J]. 湖泊科学, 2016, 28(1): 105-113.
[5]孙远军, 卢士强, 邵一平, 杨漪帆, 李柱. 城市河流水体污染物降解规律及降解系数研究[J]. 上海环境科学, 2019(24): 99-104.
基金项目：江苏省第五期“333工程”科研项目，编号：BRA2019245。
作者简介：颜建（1978—），男，江苏省宿迁市人，研究生，高级工程师，研究方向：水利工程管理与水环境等。