CN110705171A - 一种基于mike模型的感潮河网水环境治理方法 - Google Patents

Abstract

本发明提供了一种基于MIKE模型的感潮河网水环境治理方法，包括如下步骤：步骤1，基于MIKE构建河网模型，并通过空间步长控制、地形复核、以时间序列赋值群闸操作的方法，增加模型稳定性和精度；步骤2，确定补水活水方案；步骤3，确定补水过程和补水量。针对感潮河网复杂的水流特性，利用MIKE11模型的HD，AD模块构建了一维水动力与水质模型，模型可直接确定水闸和泵站的规划流量、污染源汇入点流量和浓度，输出河涌各断面的水位过程线及污染物浓度的过程线，为河网区域的补水活水、河道整治项目实施安排提供决策依据，充分利用潮汐动力和现有水闸调度补水，项目运行成本较低，有利于实施，河涌互连成片，改善河涌水环境同时有利于洪水调节。

Description

一种基于MIKE模型的感潮河网水环境治理方法

技术领域

本发明涉及水环境治理方法，特别是一种基于MIKE模型的感潮河网水环境治理方法。

背景技术

为了深入推进水污染治理工作，提升水环境质量，需对内河涌采取工程与非工程措施进行水环境整治，很多项目所在水系为感潮河网地区，水流受潮汐的影响形成双向往复流，水系纵横交错、水情极其复杂，且水工建筑众多，这些因素都大大增加了水动力模拟的难度，因此采用的数学模型需要既能够模拟其复杂的河道水动力条件，还能够进行水质计算。目前模型应用较多的WASP、QUAL2K、EFDC等，主要应用于水量模拟、污染物的模拟预测等领域，模型均用于模拟河流流量、水位、泥沙输送，并没有进行水质模拟的模块，从而无法科学分析评价利用潮汐动力及泵闸调度的补水量以及补水活水后水环境质量的改善效果。

因此需要建立新的模型进行水环境治理前后的方案策划和效果评估。

发明内容

为克服现有技术的缺陷，本发明的目的在于提供一种基于MIKE模型的感潮河网水环境治理方法，包括如下步骤：

步骤1，基于MIKE构建河网模型，并通过空间步长控制、地形复核、以时间序列赋值群闸操作的方法，增加模型的稳定性和精度；

步骤2，确定补水活水方案；

步骤3，确定补水过程和补水量。

优选地，步骤1所述河网模型建模的基本原理包括：

(1)河网模型的控制方程

研究水系属于感潮河网水系，模型控制方程为描述一维非恒定流运动的圣维南方程组：

式中，z为水位；Q为流量；B为过水断面水面宽；q为旁侧入流量；A为过水断面面积；c为谢才系数；R为水力半径；x，t为位置和时间坐标；方程组利用Abbott-Ionescu六点隐式有限差分格式离散控制方程组，该格式在每一个网格点不同时计算水位和流量，按顺序交替计算水位或流量，分别称为h点和Q点；

(2)水工结构物

研究区域建有多个闸泵等水工结构物，MIKE11将水闸作为特殊的汊点进行考虑，描述各种不同类型的水闸的调度，当水闸全开时，该汊点满足简单衔接条件，当水闸关闭或部分关闭时，水闸上下游河段作为边界河段考虑，其边界条件为流量边界，当水闸不完全闭合时根据上下游水位差和水闸的堰流或孔流流量公式计算出过闸流量。

优选地，所述步骤1包括利用水动力HD模块进行项目片区河网水动力计算，利用可控水工建筑物SO模块对河网内的水闸、泵站进行模拟，利用对流扩散AD模块进行水质计算，从而分析评价利用潮汐动力及闸泵调度的补水量及补水活水后水环境质量的改善效果。

优选地，所述步骤1所述构建数学模型包括：相应区域的大河网模型以及相应区域的河网水动力模型，在率定和验证的基础上，利用大河网模型提供边界驱动条件给河网水动力模型，并且将河网水动力模型与大河网模型的对应站点水位和流量进行对比，从而为所述河网水动力模型进行间接验证。

优选地，所述水动力HD模块进行项目片区河网水动力计算，利用可控水工建筑物SO模块对河网内的水闸、泵站进行模拟包括：

步骤11，确定模型的主要参数，包括河床糙率以及时间步长和空间步长，河网汊点处地形若有较大落差，容易导致模型发散，此时在不影响计算结果的前提下，可以人为的在两个突变的汊点断面之间以两个断面底高程均值插值出一个断面作为过渡；

步骤12，确定边界控制条件；

步骤13，根据设计资料及工程实际情况，建立的河网模型对多个水闸和多个泵站的位置、结构、调度方式进行设置，采用时间序列的方式准确赋值闸门开度和开泵流量等数据，时间精确到分钟，闸门开度精确到mm，泵站流量精确到0.01立方米/秒；

步骤14，模型验证，通过实测数据和大河网水动力模型所计算的数据，对河网水动力模型进行率定和验证，从而确保河网水动力模型的水位和流量模拟效果较好，计算值和实测值基本一致。

优选地，还包括：

步骤15，建立水质模型，所述水质模型为一维稳态水质模型，计算方程为：

C(x)＝C₀exp(-k·x/u)

式中，C(x)为流经x距离后的污染物浓度，mg/L；x为沿河段的纵向距离，m；u为设计流量下河道断面的平均流速，m/s；K为污染物综合衰减系数，l/s。

优选的，所述水质模型由河道断面及源汇项组成，其中源汇项包括沿程接纳的污染物，汇项包括相关的取水设施及河道型污水处理设施的取水口。

优选的，所述水质模型的计算参数包括：

水质指标：根据内河涌主要污染物，计算选取COD，NH₃-N2个水质指标；

水质边界：上游水质边界采用V类水标准，下游水质边界即引水水质采用实测LC入海口水质；

污染源强：补水方案按照整治后河涌条件计算，水质模型污染源不考虑底泥污染交换项等内源污染，考虑污水排放等外源污染，污染物排放量采用水污染普查工作的分析结果，浓度按一般生活污水浓度，即COD350mg/L，NH₃-N35mg/L；

污染物降解系数：参照相关地区境内河流调查分析结果，COD的讲解系数取0.1/d，NH₃-N的降解系数取0.07/d；

扩散参数：根据相近工程研究成果，扩散参数取10m²/s；

优选地，所述步骤2包括：确定水流线路，根据各河涌现行条件，结合项目实施的阶段性，计算分析多种方案不同的工况条件，确定水闸和泵站的调度运行原则。

本发明的有益效果：针对感潮河网复杂的水流特性，利用MIKE11模型的HD，AD模块构建了一维水动力与水质模型，模型可直接确定水闸和泵站的规划流量、污染源汇入点流量和浓度，输出河涌各断面的水位过程线及污染物浓度的过程线，为河网区域的补水活水、河道整治等治理项目实施安排提供决策依据，充分利用潮汐动力和现有水闸调度进行补水，项目运行成本较低，有利于实施，河涌互连成片，改善河涌水环境的同时，有利于洪水调节。

根据下文结合附图对本发明具体实施例的详细描述，本领域技术人员将会更加明了本发明的上述以及其他目的、优点和特征。

附图说明

图1表示根据本发明实施例的基于mike模型的感潮河网水环境治理方法流程图；

图2表示根据本发明实施例的河流一维水质模型模拟示意图；

图3表示根据本发明实施例的补水过程示意图；

图4表示根据本发明实施例的不同工况补水量及泵站运行时间；

图5表示根据本发明实施例的RMC水质改善效果分析示意图；

图6表示根据本发明实施例的LC水质改善效果分析示意图。

具体实施方式

参见图1，一种基于MIKE模型的感潮河网水环境治理方法，包括如下步骤：

步骤1，基于MIKE构建数学模型，并通过空间步长控制、地形复核、以时间序列赋值群闸操作的方法，增加模型的稳定性和精度；

步骤2，确定补水活水方案；

步骤3，确定补水过程和补水量。

MIKE11是由丹麦水力研究所DHI开发的一维水环境模拟程序包。水动力模型HD模块式MIKE11模型的核心，利用Abbott六点隐式格式求解河流一维非恒定流控制方程组。Abbott六点隐式格式无条件稳定，在Courant数较大的取值范围内仍可以保持计算稳定，较长的时间步长取值也可以保证模型收敛。河网模型建模的基本原理包括：

(1)河网模型的控制方程

研究水系属于感潮河网水系，模型控制方程为描述一维非恒定流运动的圣维南方程组：

式中，z为水位；Q为流量；B为过水断面水面宽；q为旁侧入流量；A为过水断面面积；c为谢才系数；R为水力半径；x，t为位置和时间坐标。计算的上游流量边界采用多年的平均径流量，下游潮位边界采用2018年2月某闸站实测潮位过程。

方程组利用Abbott-Ionescu六点隐式有限差分格式离散控制方程组，该格式在每一个网格点不同时计算水位和流量，按顺序交替计算水位或流量，分别称为h点和Q点。Abbott-Ionescu格式稳定性好、计算精度高，离散后用追赶法求解线性方程组。

(2)水工结构物

研究区域建有多个闸泵等水工结构物，MIKE11将水闸作为特殊的汊点进行考虑，描述各种不同类型的水闸的调度。当水闸全开时，该汊点满足简单衔接条件。当水闸关闭或部分关闭时，水闸上下游河段作为边界河段考虑，其边界条件为流量边界。当水闸不完全闭合时根据上下游水位差和水闸的堰流或孔流流量公式计算出过闸流量。

步骤1包括利用水动力HD模块进行项目片区河网水动力计算，利用可控水工建筑物SO模块对河网内的水闸、泵站进行模拟，利用对流扩散AD模块进行水质计算，从而分析评价利用潮汐动力及闸泵调度的补水量及补水活水后水环境质量的改善效果。

由于缺乏匹配的边界驱动水文数据，为获得相应区域河网水动力模型的边界驱动条件，构建相应区域的大河网模型。步骤1所述构建数学模型包括：相应区域的大河网模型以及相应区域的河网水动力模型。在率定和验证的基础上，利用大河网模型提供边界驱动条件给河网水动力模型，并且将河网水动力模型与大河网模型的对应站点水位和流量进行对比，从而为所述河网水动力模型进行间接验证。

水动力HD模块进行项目片区河网水动力计算，利用可控水工建筑物SO模块对河网内的水闸、泵站进行模拟包括：

步骤11，确定模型的主要参数，包括河床糙率以及时间步长和空间步长，本实施例中，结合各河涌地形图，同时选取目标区域河网水动力数学模型同步的水文测验资料和实测数据进行模型参数率定，目标区域河网水动力模型的糙率值范围为0.025-0.030。对于时间步长和空间步长，MIKE11模型采用Abbott六点隐式格式，理论上对时间步长没有限制，但建模过程表明，过大的时间补偿，容易使计算结果过于坦化而失真，而过小的时间步长则容易使得某些非线性的小扰动得到响应而导致计算失稳。同样，过大的空间步长无法真实反映河道现状，但过密的河道断面容易导致模型发散。考虑到模型计算稳定性要求以及计算效率，河网水动里模型的时间步长设定为5秒，河道地形断面间距根据河涌不同控制在50-500m之间。此外，河网汊点处地形若有较大落差，容易导致模型发散，此时在不影响计算结果的前提下，可以人为的在两个突变的汊点断面之间以两个断面底高程均值插值出一个断面作为过渡。

步骤12，确定边界控制条件，本实施例河网模型采用的组合式上游边界黄沙，沙洛围和大石设定流量过程线，下边界黄埔左和黄埔右设定水位过程线，流量和水位过程均从大河网水动力模型提取。

步骤13，根据设计资料及工程实际情况，建立的河网模型对40个水闸和6个泵站的位置、结构、调度方式进行设置，为了能准确还原闸泵调度操作，采用时间序列的方式准确赋值闸门开度和开泵流量等数据，时间精确到分钟，闸门开度精确到mm，泵站流量精确到0.01立方米/秒。计算结果表明，通过时间序列赋值闸泵操作的方式，比通过调度判断条件的水工模拟方式，模拟精度更高。

步骤14，模型验证，通过实测数据和大河网水动力模型所计算的数据，对河网水动力模型进行率定和验证，从而确保河网水动力模型的水位和流量模拟效果较好，计算值和实测值基本一致。

步骤15，建立水质模型，所述水质模型为一维稳态水质模型，计算方程为：

C(x)＝C₀exp(-k·x/u)

式中，C(x)为流经x距离后的污染物浓度，mg/L；x为沿河段的纵向距离，m；u为设计流量下河道断面的平均流速，m/s；K为污染物综合衰减系数，l/s。

模型主要由河道断面及源汇项组成，其中源汇项主要包括沿程接纳的污染物，汇项主要包括相关的取水设施及河道型污水处理设施的取水口等，河流一维水质模型模拟概化图如图2所示。

水质模型的计算参数包括：

水质指标：根据内河涌主要污染物，计算选取COD，NH₃-N2个水质指标。

水质边界：上游水质边界采用V类水标准，下游水质边界即引水水质采用实测LC入海口水质。

污染源强：本实施例中目前河涌进行边坡整治和清淤疏浚工作，补水方案按照整治后河涌条件计算，水质模型污染源不考虑底泥污染交换项等内源污染，主要考虑污水排放等外源污染。污染物排放量采用水污染普查工作的分析结果，浓度按一般生活污水浓度，即COD350mg/L，NH₃-N 35mg/L。

污染物降解系数：污染物在河道中的降解系数与很多因素有关。本实施例水质模拟各污染物降解系统参照相关地区境内河流调查分析结果，COD的讲解系数取0.1/d，NH₃-N的降解系数取0.07/d。

扩散参数：根据相近工程研究成果，扩散参数取10m²/s。

步骤2包括：确定水流线路，根据各河涌现行条件，结合项目实施的阶段性，计算分析多种方案不同的工况条件，确定水闸和泵站的调度运行原则。

本实施例中，海水由LC入内河涌，经连通箱涵进入RMC，由RMC、KZC流入MZH，最后返回海湾。对于水闸的调度运行原则，外江潮位盖于内河涌水位时，打开LC水闸，关闭RMC水闸，利用外海与内河涌的水位差，将海湾潮水引入LC。当外江潮位下降且低于内河涌水位时，关闭LC水闸，待MZH水位低于内河涌水位时，打开RMC水闸排水。对于泵站的调度运行，当RMC内水位高于0m且低于外江水位，且LC水闸进水时，开启RMC泵站将内河涌的水抽向MZH。参见表1，获得4种不同的工况条件。

表1 4种方案不同工况条件

对于步骤3，本实施例中：

(1)补水过程：此海域的潮汐属于不规则半日潮，即在一个太阴日内出现两次高潮和两次低潮，其潮高、潮差和潮历时各不相等。随着潮位的变化，每天进行2次补水，如图3所示。

(2)补水量：4个方案不同工况条件的补水量结果参见表2和图4

表2 4种方案补水倍数

仅利用潮汐动力和水闸的调度，日均可引潮量达到20多万立方米，约相当于内河涌每天换水1.5次，启用RMC泵站后，引潮量翻倍，但会增加泵站运行费用。对比方案1和2，方案3和4，泵站的使用能显著增加补水量，这是因为泵站能将内河涌的水及时抽排至外江，增大外海与内河涌的水位差，从而增大补水量。对比方案1和3，方案2和4，打通KZC断头涌能增大补水量，但由于KZC水闸排水能力有限，补水量增大不明显，但KZC目前为断头涌，将其与RMC连通后，能使死水变活水，同时改善KZC的水环境。

本实施例通过MIKE11水质模型对河涌补水的水质状况进行分析，结果参见图5和图6。通过补水，RMC水质得到明显改善，由于补水量是随着潮汐变化而变化，不同时段的水质有所差别。COD最小值能达到V类水标准。COD平均值从超V类水3倍改善到超1倍之内，NH₃-N从超V类水10倍改善到超2-3倍。NH₃-N最小值仍超V类，因为补水来源LC入海口的海水NH₃-N超V类水标准。LC现状水质优于RMC。通过补水，LC的COD平均值从超V类水2倍改善到接近V类水标准，NH₃-N从超V类水7倍改善到超2倍以内。COD改善效果优于NH₃-N改善效果。通过截污治污工程减少排入内河涌的污染物，结合流域综合整治的实施所带来的外海水质的改善，补水效果更加明显，可以达到V类水标准。

虽然本发明已经参考特定的说明性实施例进行了描述，但是不会受到这些实施例的限定而仅仅受到附加权利要求的限定。本领域技术人员应当理解可以在不偏离本发明的保护范围和精神的情况下对本发明的实施例能够进行改动和修改。

Claims (9)

1.一种基于MIKE模型的感潮河网水环境治理方法，其特征在于包括如下步骤：

步骤1，基于MIKE构建河网模型，并通过空间步长控制、地形复核、以时间序列赋值群闸操作的方法，增加模型的稳定性和精度；

步骤2，确定补水活水方案；

步骤3，确定补水过程和补水量。

2.根据权利要求1所述的一种基于MIKE模型的感潮河网水环境治理方法，其特征在于步骤1所述河网模型建模的基本原理包括：

(1)河网模型的控制方程

研究水系属于感潮河网水系，模型控制方程为描述一维非恒定流运动的圣维南方程组：

式中，z为水位；Q为流量；B为过水断面水面宽；q为旁侧入流量；A为过水断面面积；c为谢才系数；R为水力半径；x，t为位置和时间坐标；方程组利用Abbott-Ionescu六点隐式有限差分格式离散控制方程组，该格式在每一个网格点不同时计算水位和流量，按顺序交替计算水位或流量，分别称为h点和Q点；

(2)水工结构物

研究区域建有多个闸泵等水工结构物，MIKE11将水闸作为特殊的汊点进行考虑，描述各种不同类型的水闸的调度，当水闸全开时，该汊点满足简单衔接条件，当水闸关闭或部分关闭时，水闸上下游河段作为边界河段考虑，其边界条件为流量边界，当水闸不完全闭合时根据上下游水位差和水闸的堰流或孔流流量公式计算出过闸流量。

3.根据权利要求1所述的一种基于MIKE模型的感潮河网水环境治理方法，其特征在于所述步骤1包括利用水动力HD模块进行项目片区河网水动力计算，利用可控水工建筑物SO模块对河网内的水闸、泵站进行模拟，利用对流扩散AD模块进行水质计算，从而分析评价利用潮汐动力及闸泵调度的补水量及补水活水后水环境质量的改善效果。

4.根据权利要求1所述的一种基于MIKE模型的感潮河网水环境治理方法，其特征在于所述步骤1所述构建数学模型包括：相应区域的大河网模型以及相应区域的河网水动力模型，在率定和验证的基础上，利用大河网模型提供边界驱动条件给河网水动力模型，并且将河网水动力模型与大河网模型的对应站点水位和流量进行对比，从而为所述河网水动力模型进行间接验证。

5.根据权利要求4所述的一种基于MIKE模型的感潮河网水环境治理方法，其特征在于所述水动力HD模块进行项目片区河网水动力计算，利用可控水工建筑物SO模块对河网内的水闸、泵站进行模拟包括：

步骤11，确定模型的主要参数，包括河床糙率以及时间步长和空间步长，河网汊点处地形若有较大落差，容易导致模型发散，此时在不影响计算结果的前提下，可以人为的在两个突变的汊点断面之间以两个断面底高程均值插值出一个断面作为过渡；

步骤12，确定边界控制条件；

步骤13，根据设计资料及工程实际情况，建立的河网模型对多个水闸和多个泵站的位置、结构、调度方式进行设置，采用时间序列的方式准确赋值闸门开度和开泵流量等数据，时间精确到分钟，闸门开度精确到mm，泵站流量精确到0.01立方米/秒；

步骤14，模型验证，通过实测数据和大河网水动力模型所计算的数据，对河网水动力模型进行率定和验证，从而确保河网水动力模型的水位和流量模拟效果较好，计算值和实测值基本一致。

6.根据权利要求5所述的一种基于MIKE模型的感潮河网水环境治理方法，其特征在于还包括：

步骤15，建立水质模型，所述水质模型为一维稳态水质模型，计算方程为：

C(x)＝C₀exp(-k·x/u)

式中，C(x)为流经x距离后的污染物浓度，mg/L；x为沿河段的纵向距离，m；u为设计流量下河道断面的平均流速，m/s；K为污染物综合衰减系数，l/s。

7.根据权利要求6所述的一种基于MIKE模型的感潮河网水环境治理方法，其特征在于所述水质模型由河道断面及源汇项组成，其中源汇项包括沿程接纳的污染物，汇项包括相关的取水设施及河道型污水处理设施的取水口。

8.根据权利要求6所述的一种基于MIKE模型的感潮河网水环境治理方法，其特征在于所述水质模型的计算参数包括：

水质指标：根据内河涌主要污染物，计算选取COD，NH₃-N2个水质指标；

水质边界：上游水质边界采用V类水标准，下游水质边界即引水水质采用实测LC入海口水质；

污染源强：补水方案按照整治后河涌条件计算，水质模型污染源不考虑底泥污染交换项等内源污染，考虑污水排放等外源污染，污染物排放量采用水污染普查工作的分析结果，浓度按一般生活污水浓度，即COD350mg/L，NH₃-N 35mg/L；

污染物降解系数：参照相关地区境内河流调查分析结果，COD的讲解系数取0.1/d，NH₃-N的降解系数取0.07/d；

扩散参数：根据相近工程研究成果，扩散参数取10m²/s。

9.根据权利要求1所述的一种基于MIKE模型的感潮河网水环境治理方法，其特征在于所述步骤2包括：确定水流线路，根据各河涌现行条件，结合项目实施的阶段性，计算分析多种方案不同的工况条件，确定水闸和泵站的调度运行原则。

Images