平原水网圩区-圩外系统滞蓄关系优化方法及系统与流程

1.本发明涉及水文技术领域，具体涉及一种平原水网圩区-圩外系统滞蓄关系优化方法。


背景技术：

2.世界上的多数平原水网区气候湿润，降水充沛，河流纵横，河网高度发育且没有明确的流域出口，蓄泄关系复杂，同时地势低洼，因此防洪除涝问题突出。当前并没有很好的解决方案。
3.为了防洪排涝，建设了诸多圩区和城市防洪大包围，圩区(包括城市防洪大包围)建设通过筑堤、建闸、设排涝站排水等方法，防止圩外河道洪水侵袭，排除圩区内部涝水。但大量的圩区建设改变了区域洪涝水自然蓄泄特性及其在圩外河网、圩内河网的时空分配特征，当降水强度较大时圩区集中向圩外排水，圩外骨干河道排水能力有限，往往导致圩外河网水位抬升迅速，增加区域整体的防洪风险。
4.对于某个平原水网圩区-圩外水网系统来说，一场降雨下该区域洪涝水风险大小主要表现为河网水位的高低，局部时段区域河网水位偏高、洪涝风险较大，除受降雨量大等非人为可控因素的影响外，更多的是圩区-圩外水网系统滞蓄关系不当、区域洪涝水时空分配不合理的结果。而现有技术方法不能较好解决上述问题。
5.因此，在区域经济社会快速发展要求不断提升防洪除涝安全度的大背景下，寻求一种简单、可靠、可定量的平原水网圩区-圩外系统洪涝水滞蓄关系预测、优化方法，对于平原水网区洪涝风险管理显得尤为重要，也是目前区域防洪工程调度管理亟待解决的问题。


技术实现要素：

6.发明目的：针对现有技术无法模拟预测、判断并优化实时情景下水网区滞蓄关系的不足，提出一种平原水网圩区-圩外系统滞蓄关系优化方法，进一步的目的是提出一种基于上述方法的系统。
7.技术方案：
8.第一方面，提供一种平原水网圩区-圩外系统滞蓄关系优化方法，包括如下步骤：
9.s1、采集研究区域的历史水文和降雨数据，并基于所述历史水文和降雨数据识别场次降雨涨水期，建立典型历史情景集合并分类，基于河网-圩区拓扑关系，建立水文水动力模型，通过水文水动力模型计算模拟各个典型历史情景下河网水位、河网外排水量，以及各时刻下圩内与圩外河网的槽蓄量；
10.s2、构建平原水网区蓄泄特征的表征指标，所述表征指标包含典型历史情景下的区域蓄泄比sdrj、区域蓄泄比sdr、圩外河网单位面积滞蓄水量as1、圩内河网单位面积滞蓄水量as2和防洪风险指数r；
11.s3、基于典型历史情景下区域外排水量和槽蓄量的计算结果，计算典型历史情景下的区域蓄泄比sdrj，分别绘制典型历史情景下区域蓄泄比sdrj与初始水位z
(j,1)
、累计降
雨量pj的关系图，sdrj＝f(z
(j,1)
，pj)；并进行区域蓄泄比sdr与上述水文要素的相关性检验，通过相关性检验后，以各个(z
(j,1)
，ln(sdrj))为参数建立区域蓄泄比sdr～水文要素的关系；
12.s4、确定圩内、圩外单位面积滞蓄水量比的阈值c
as
；
13.计算各典型历史情景下圩外单位面积滞蓄水量as1、圩内单位面积滞蓄水量as2、圩内与圩外单位面积滞蓄水量比as2/as1，并进一步确定圩内与圩外单位面积滞蓄水量比as2/as1的阈值c
as
；
14.s5、模拟预测实时情景圩区-圩外系统滞蓄关系，并根据预测计算的实时情景防洪风险指数r，区域蓄泄比sdr～水文要素的关系，圩内与圩外单位面积滞蓄水量比as2/as1与阈值c
as
的比较，判断是否需要优化实时情景下的圩区-圩外系统滞蓄关系；
15.s6、针对防洪风险指数r》0同时圩内与圩外单位面积滞蓄水量比as2/as1＜阈值c
as
的实时情景，采用系统聚类方法对圩区进行聚类分析，再通过增加不同类型圩区调蓄水深hs的方式来优化圩区-圩外系统滞蓄关系，以降低防洪风险指数r。
16.根据本技术的一个方面，所述阈值c
as
的确定方法如下：
17.在聚类分析后生成的各类典型历史情景中，取满足防洪风险指数r＝0的若干典型历史情景中圩内与圩外单位面积滞蓄水量比as2/as1的均值作为该类情景下的c
as
，
[0018][0019]
式中，k为表示情景类别的自然数，c
as,k
为第k类情景对应的圩内与圩外单位面积滞蓄水量比阈值，j为表示情景序号的自然数，as
1,k,j
为属于k类情景下的情景j的圩外河网单位面积滞蓄水量，as
2,k,j
为属于k类情景下的情景j的圩内河网单位面积滞蓄水量，r
k,j
为k类情景中情景j的防洪风险指数，dk为k类情景中满足防洪风险指数r
k,j
＝0的情景个数。
[0020]
根据本技术的一个方面，所述s1中所述识别场次降雨涨水期，建立典型历史情景集合的过程，具体为：
[0021]
输入历史降雨、水位数据，场次降水开始的时刻记为t1，
[0022]
t1时刻水位站i的水位记为z
1,i
，区域河网的平均水位记为z1，z1取各z
1,i
的算术平均值，
[0023]
场次降水中河网水位达到峰值的时刻记为t2，
[0024]
t2时刻对应的水位为峰值水位z
max
，
[0025]
从t1到t2的降水过程定义为一个场次降水事件的涨水期，t1、t2满足以下条件：
[0026]
t1之前特定时段
△
t内平均降雨强度小于1mm/d；
[0027]
t1～t2时段内单日最大降雨量p
1d
≥10mm；
[0028]
由此建立河网初始水位z1、场次降水涨水期降雨量p的二维情景集合{z1，p}。
[0029]
根据本技术的一个方面，步骤s1所述对建立的情景集合进行分类，具体为：
[0030]
采用系统聚类方法对构建的典型历史情景集合进行分类，聚类样本为构建的若干历史情景，聚类指标为初始水位均值z1。
[0031]
根据本技术的一个方面，步骤s1所述水文水动力模型构建及计算，具体为：
[0032]
收集平原水网区地形、下垫面、水系、圩区、水文站点资料；
[0033]
构建预研究区域的圩区与圩区、圩区与河网的拓扑关系；
[0034]
基于上述拓扑关系，查找独立圩区、嵌套圩区与可分割圩区；
[0035]
将可分割圩区拆分成等效为独立圩区的若干个子圩区，并将所述子圩区插入所述独立圩区形成的集合；
[0036]
查找所述嵌套圩区的层级，并基于层级关系，标定内部圩区为二级圩区、三级圩区至n级圩区；
[0037]
构建包括产汇流子模型、河网水动力子模型的水文水动力耦合模型，并采用水文站历史水位数据对模型进行率定验证；
[0038]
采用经率定的水文水动力模型模拟计算典型历史情景下的平原水网区的河网水位、河网槽蓄量和河网外排水量。
[0039]
根据本技术的一个方面，步骤s2中所述构建平原水网区蓄泄特征表征指标中，
[0040]
所述区域蓄泄比sdrj为sdrj＝sj/w
(j,out)
[0041]
式中，sdrj为典型历史场景j下的区域蓄泄比；sj为典型历史场景j下的区域总滞蓄水量，包含圩外河网、圩内河网的滞蓄水量；w
(j,out)
为典型历史场景j下的区域河网外排水量。
[0042]
所述圩外河网单位面积滞蓄水量as
j1
为as
j1
＝s
j1
/u1；
[0043]
所述圩内河网单位面积滞蓄水量as
j2
为as j2
＝s
j2
/u2；
[0044]sj1
+s
j2
＝sj；
[0045]
式中，as
j1
为历史情景j下的圩外河网单位面积滞蓄水量，as
j2
为历史情景j下的圩内河网单位面积滞蓄水量，u1为圩外河网面积，u2为圩区面积，s
j1
为历史情景j下的圩外河网滞蓄水量，s j2
为历史情景j下的圩内河网滞蓄水量；
[0046]
为终止时刻t2时圩外的槽蓄量，为初始时刻t1时圩外的槽蓄量，为终止时刻t2时圩内的槽蓄量，为初始时刻t1时圩内的槽蓄量；其余变量含义同前；
[0047]
所述防洪风险指数r为：
[0048][0049][0050][0051]
式中：ri为水位站i的防洪风险指数，r为区域防洪风险指数，zi(t)为水位站i的水位过程，hi为水位站i的保证水位，n为站点数量；其余变量含义同前；
[0052]
当r》0时，表示区域内部分或全部水位站水位超过保证水位，区域存在一定防洪风险。
[0053]
根据本技术的一个方面，所述步骤s5进一步包括：
[0054]
s51、实时情景圩区-圩外系统滞蓄关系模拟预测：基于获取的包括河网水位、预报降雨在内的实时数据，采用所述水文水动力模型模拟预测实时情景下的河网水位、河网槽蓄量和河网外排水量，并计算得到实时情景下的区域蓄泄比sdr，圩内与圩外单位面积滞蓄水量比as2/as1、以及防洪风险指数r；
[0055]
s52、实时情景圩区-圩外系统滞蓄关系判断：首先判断实时情景防洪风险指数r值大小，若防洪风险指数r＝0，则完成实时情景圩区-圩外系统滞蓄关系计算；
[0056]
若防洪风险指数r＞0，则进一步判断圩内与圩外单位面积滞蓄水量比as2/as1与阈值c
as
大小，若圩内与圩外单位面积滞蓄水量比as2/as1＜阈值c
as
，则对实时情景圩区-圩外系统滞蓄关系进行优化，否则完成实时情景圩区-圩外系统滞蓄关系计算。
[0057]
根据本技术的一个方面，所述s6中所述圩区聚类分析方法具体为：
[0058]
构建反映圩区自然属性与社会属性的指标，圩区自然属性包括圩内水面率、堤顶高程和排涝模数，其中圩内水面率表征圩区自然状态下调蓄能力，堤顶高程差异间接反映圩区地面高程差异，排涝模数表征工程作用下圩区调蓄能力；
[0059]
社会属性指标用于区分圩区农业圩或非农业圩属性，采用圩内耕地面积占比间接表征圩区社会属性，圩内耕地面积占比越大，表示该圩区农业圩属性越大；
[0060]
基于上述自然属性和社会属性指标，采用系统聚类分析法将区域内圩区分为若干类型。
[0061]
根据本技术的一个方面，所述s6中所述实时情景滞蓄系优化，具体为：
[0062]
针对防洪风险指数r》0同时圩内与圩外单位面积滞蓄水量比as2/as1＜阈值c
as
的实时情景，在保证圩区自身安全的前提下，对于不同类型的圩区不同程度地增加圩区调蓄水深hs，以增加涨水期圩内河网滞蓄水量，优化洪涝水在圩外河网、圩内河网的时空分配，从而实现圩区-圩外系统滞蓄关系优化，降低实时情景防洪风险指数。
[0063]
根据本技术的一个方面，
[0064]
构建预研究区域的圩区与圩区、圩区与河网的拓扑关系；
[0065]
基于上述拓扑关系，查找独立圩区、嵌套圩区与可分割圩区；
[0066]
将可分割圩区拆分成等效为独立圩区的若干个子圩区，并将所述子圩区插入所述独立圩区形成的集合；
[0067]
查找所述嵌套圩区的层级，并基于层级关系，标定内部圩区为二级圩区、三级圩区至n级圩区；基于已标定的圩区层级划定圩区调控的优先级。
[0068]
在进一步的实施例中，判断是否属于可分割圩区的方法是，内部圩区与外部圩区共边(圩堰)，且内部圩区包含于外部圩区；判断是否属于嵌套圩区的方法是，外部圩区全部包含内部圩区，且内部圩区与外部圩区没有重叠的圩堰；判断两个圩区为独立圩区的方法是，两个圩区互不包含，可共用部分圩堰。
[0069]
在进一步的实施例中，在实时情景滞蓄优化时，按照圩区的层级，从内之外，逐级优化；圩内的面积为圩堰包含的所有面积。
[0070]
因此当一个大圩区包含小圩区时，在计算大圩区圩内面积时，被包含的小圩区的面积也统计在其中。第二方面，提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述实施例中任一项所述方法的步骤。
[0071]
有益效果：本发明针对某些平原水网区圩区众多，遭遇降雨时圩区-圩外系统滞蓄关系不合理导致的河网水位迅速升高影响防洪安全的问题，提出了一种平原水网圩区-圩外系统滞蓄关系优化方法，可以自动模拟预测实时情景下洪涝水的时空分配，判断并优化实时情景下水网区圩区-圩外系统的滞蓄关系，从而指导实时防洪调度，降低区域洪涝风险。
附图说明
[0072]
图1是河网-圩区拓扑关系图((a)为圩区-圩外河网拓扑关系图，(b)为单个圩区拓扑关系图)。
[0073]
图2是本发明的技术流程图。
[0074]
图3是本发明的系统架构图。
[0075]
图4是本发明实施例一圩区聚类结果。
[0076]
图5是本发明实施例一蓄泄比sdr～z1关系图。
[0077]
图6是本发明实施例一滞蓄关系优化对比图。
[0078]
图7a和图7b是本发明两种二级嵌套圩区示意图。
[0079]
图8是本发明实施例三水文水动力模型技术框架图。
具体实施方式
[0080]
在下文的描述中，给出了大量具体的细节以便提供对本发明更为彻底的理解。然而，对于本领域技术人员而言显而易见的是，本发明可以无需一个或多个这些细节而得以实施。在其他的例子中，为了避免与本发明发生混淆，对于本领域公知的一些技术特征未进行描述。
[0081]
为了解决上述问题，申请人进行了深入地研究。
[0082]
如图1中的a和b所示，对于平原河网地区，城市防洪工程和圩区建设是提高低洼地区防洪标准最直接的措施之一。此时某个区域河网的调蓄能力既包含了圩外河网的调蓄能力，也包含了圩区及城市大包围的调蓄能力。城防、圩区建设通过圈圩筑堤、建闸控制、设站排水等方法，防止外河洪水侵袭，排除圩内涝水，从而达到防洪和排涝的目的。但同时圩区建设使得各圩区自成封闭系统，诸多河道成为圩内河道，其涝水则主要排入周边的圩外河道，从而改变了自然状态下圩内、圩外洪涝水的汇流过程和时空分配。
[0083]
通过发明人的研究，认为在一定量级降雨下，降低某个河网区域洪涝风险关键在于，定量模拟和预测实时情景下洪涝水在圩区-圩外系统中的时空分配特征和河网防洪风险大小，并做出调整和优化，以求在圩区安全的前提下，圩外河网水位不至于过高，从而达到优化区域洪涝水在圩区-圩外系统中的时空分配、区域整体防洪安全最大化(防洪风险指数r最小)的目的。
[0084]
实施例一：
[0085]
本实施例以长江下游某平原河网区为例，应用本发明提出的一种平原水网圩区-圩外系统滞蓄关系优化方法来计算分析平原水网区典型历史场次降水事件涨水期圩区-圩外系统滞蓄关系，并进一步优化实时情景下的圩区-圩外系统滞蓄关系。结合附图3，其具体步骤如下：
[0086]
步骤1，根据历史水文、降雨数据识别场次降雨涨水期，建立典型历史情景集合，并对建立的情景集合进行分类，通过水文水动力模型计算模拟各历史情景下河网水位、区域外排水量，以及各时刻下圩内与圩外的槽蓄量等要素；
[0087]
步骤2，构建平原水网区蓄泄特征表征指标，包含区域蓄泄比sdrj、区域蓄泄比sdr、圩外河网单位面积滞蓄水量as1、圩内河网单位面积滞蓄水量as2和防洪风险指数r；
[0088]
步骤3，基于典型历史情景下河网水位、外排水量、槽蓄量的计算过结果，计算绘制典型历史情景j下区域蓄泄比sdrj，并分别绘制典型历史情景下sdrj与初始水位z
(j,1)
、累计降雨量pj的关系图，sdrj＝f(z
(j,1)
，pj)，并进行sdrj与上述水文要素的相关性检验，通过相关性检验后，以各个(z
(j,1)
，ln(sdrj))为参数建立区域蓄泄比sdr～水文要素的关系，为后续预测提供数据支撑。
[0089]
步骤4，确定圩内、圩外单位面积滞蓄水量比的阈值c
as
；
[0090]
计算各历史情景下圩外单位面积滞蓄水量as1、圩内单位面积滞蓄水量as2、圩内与圩外单位面积滞蓄水量比as2/as1，并进一步确定圩内与圩外单位面积滞蓄水量比as2/as1的阈值c
as
。确定圩内与圩外单位面积滞蓄水量比阈值c
as
的目的是寻求一定的水文条件下适宜的阈值c
as
值，作为某个河网区是否实现滞蓄有度的判断依据。阈值c
as
的确定方法如下：
[0091]
在聚类后生成的各类情景中，取该类情景中满足区域防洪风险指数r＝0的若干历史情景下圩内与圩外单位面积滞蓄水量比as2/as1的均值作为阈值c
as
。
[0092]
需要说明的是，尽管通过系统聚类方法将历史情景做了分类，每一个类别中水文要素相对较为接近，但由于各水位站初始水位空间分布差异以及降雨等要素的影响，即使同属于同一个类别的历史情景和未来实时情景，其河网产水量时空分布特征、圩外河网及圩区滞蓄量大小仍可能存在一定差异，并进一步导致区域防洪风险大小的差异，为此，基于偏安全的考虑，本发明采用满足防洪风险指数r＝0的若干历史情景中圩内与圩外单位面积滞蓄水量比as2/as1的均值作为该类情景下的阈值c
as
。
[0093][0094]
式中，k为表示情景类别的自然数，c
as,k
为第k类情景对应的圩内、圩外单位面积滞蓄水量比阈值，j为表示情景序号的自然数，as
1,k,j
为属于k类情景下的情景j的圩内河网单位面积滞蓄水量，as
2,k,j
为属于k类情景下的情景j的圩内河网单位面积滞蓄水量，r
k,j
为k类情景中情景j的防洪风险指数，dk为k类情景中满足防洪风险指数r
k,j
＝0的情景个数。
[0095]
步骤5，获取实时水文数据，模拟预测实时情景圩区-圩外系统滞蓄关系，并根据预测计算的实时情景防洪风险指数r，圩内与圩外单位面积滞蓄水量比as2/as1与阈值的比较，判断是否需要优化实时情景圩区-圩外系统滞蓄关系；
[0096]
步骤6，针对防洪风险指数r》0同时圩内与圩外单位面积滞蓄水量比as2/as1＜阈值c
as
的实时情景，首先采用系统聚类方法对圩区进行聚类，再通过增加不同类型圩区调蓄水深hs的方式优化圩区-圩外系统滞蓄关系，降低防洪风险指数r。
[0097]
在进一步的实施例中，根据步骤1所述方法，输入历史降雨、水位数据，场次降水开
始的时刻记为t1，t1时刻水位站i的水位记为z
1,i
，区域河网的平均水位记为z1，z1取各z
1,i
的算术平均值，场次降水中河网水位达到峰值的时刻记为t2，t2时刻对应的峰值水位记为z
max
，则t1到t2的降水过程定义为一个场次降水事件的涨水期，t1、t2满足以下条件：
[0098]
t1之前特定时段
△
t内平均降雨强度小于1mm/d；
[0099]
t1～t2时段内单日最大降雨量p
1d
≥10mm。
[0100]
由此建立河网初始水位z1、场次降水涨水期降雨量p的二维情景集合{z1，p}，本实施例中构建了满足上述条件的历史情景共计24个。
[0101]
采用系统聚类方法对构建的24个历史情景集合进行聚类分析，聚类指标为初始水位均值z1。据此，将构建的情景分为3个类别，第一类情景为z1≤3.51m，第二类情景为3.51m＜z1≤3.81m，第三类情景为z1＞3.81m，见表1。
[0102]
表1历史情景聚类结果
[0103][0104]
由于难以获取区域内完整的河网外排水量、河网槽蓄量等历史资料，在进一步的实施例中，采用构建水文水动力模型对上述二维情景进行模拟计算的方式来弥补历史实测数据的不足。收集平原水网区地形、下垫面、水系、圩区(含城市防洪大包围)、水文站点等资料，构建包括产汇流子模型、河网水动力子模型的水文水动力耦合模型，并采用水文站历史
水位数据对模型进行率定验证。采用经率定的水文水动力模型计算历史情景下的平原水网区河网水位、河网槽蓄量、河网外排水量等要素。
[0105]
在进一步的实施例中，步骤2建立平原水网区蓄泄特征表征指标，分别为典型历史情景下的区域蓄泄比sdrj、圩外河网单位面积滞蓄水量as1、圩内河网单位面积滞蓄水量as2、防洪风险指数r，各指标含义及计算方法如下：
[0106]
1)区域蓄泄比sdr
[0107]
典型情景河网总滞蓄水量与区域外排水量的比值，反映一定时段内河网(含圩内河网)发挥的调蓄作用大小，sdr越大表示河网发挥的调蓄作用越大。
[0108]
sdrj＝sj/w
(j,out)
[0109]
式中，sdrj为典型历史情景j下的区域蓄泄比；sj为典型历史情景j下的区域总滞蓄水量，包含圩外河网、圩内河网的滞蓄水量；w
(j,out)
为典型历史情景j下的区域外排水量。
[0110]
2)圩外河网单位面积滞蓄水量as1、圩内河网单位面积滞蓄水量as2[0111]
由于圩外河网、圩内河网滞蓄水量的能力与其面积有关，因此，为客观反映圩外河网、圩内河网发挥的相对调蓄作用大小，建立圩外河网、圩内河网单位面积滞蓄水量指标as
j1
、as
j2
，客观衡量圩区的调蓄作用贡献，避免圩外河网、圩区面积因素对于调蓄作用大小的影响。
[0112]
所述圩外河网单位面积滞蓄水量as
j1
为as
j1
＝s
j1
/u1；
[0113]
所述圩内河网单位面积滞蓄水量as
j2
为as j2
＝s
j2
/u2；
[0114]sj1
+s
j2
＝sj；
[0115]
式中，as1为圩外河网单位面积滞蓄水量，as2为圩内河网单位面积滞蓄水量，u1为圩外河网面积，u2为圩区面积，s
j1
为圩外河网滞蓄水量，s j2
为圩内河网滞蓄水量；为终止时刻t2时圩外的槽蓄量，为初始时刻t1时圩外的槽蓄量，为终止时刻t2时圩内的槽蓄量，为初始时刻t1时圩内的槽蓄量；其余变量含义同前。
[0116]
3)防洪风险指数r
[0117]
一场降雨下，河网水位上涨是区域河网槽蓄量增加的结果和反映。为衡量不同水文情景下圩外河网蓄水状态和防洪除涝风险差异，构建防洪风险指数r。通常认为某个水位站水位处于保证水位以下时，防洪风险基本可控，反之则认为防洪风险较大，同时防洪风险又与河网水位超保证水位历时有关，因此，防洪风险指数r按下式计算：
[0118][0119][0120][0121]
式中：ri为水位站i的防洪风险指数，r为区域防洪风险指数，zi(t)为水位站i的水位过程，hi为水位站i的保证水位，n为站点数量；其余变量含义同前。
[0122]
当r》0时，表示区域内部分或全部水位站水位超过保证水位，区域存在一定防洪风
险。
[0123]
在进一步的实施例中，步骤3具体为：绘制典型情景下区域蓄泄比sdrj与河网初始水位z
(j,1)
、pj的关系图，并进行区域蓄泄比sdrj与河网初始水位z
(j,1)
、场次降水涨水期降雨量pj的相关性检验。本实施例中，经相关性检验，区域蓄泄比sdr与河网初始水位z
(j,1)
通过显著性水平α＝0.05的相关性检验，sdrj与场次降水涨水期降雨量pj未通过显著性水平α＝0.05的相关性检验，故建立sdr～z1的关系(详见图4)如下：
[0124]
ln(sdr)＝a*z
1+b[0125]
式中，a、b为拟合后的参数，本实施例中，a、b分别为-2.570、8.796。
[0126]
在进一步的实施例中，步骤4根据前述建立的3个历史情景类别，分别计算各类历史情景中满足防洪风险指数r＝0的as2/as1的均值，作为该类情景圩内与圩外单位面积滞蓄水量比的阈值c
as
。本实施例中，第一类情景(z1≤3.51m)圩内与圩外单位面积滞蓄水量比阈值c
as,1
＝0.495，第二类情景(3.51m＜z1≤3.81m)圩内与圩外单位面积滞蓄水量比阈值c
as,2
＝0.651，第三类情景(z1＞3.81m)圩内与圩外单位面积滞蓄水量比阈值c
as,3
＝0.127。
[0127]
表2各类历史情中圩内、圩外单位面积滞蓄水量比阈值
[0128]
分类as2/as1范围阈值c
as
第一类情景0.314～0.597c
as,1
＝0.495第二类情景0.3～1.711c
as,2
＝0.651第三类情景0.049～0.246c
as,3
＝0.127
[0129]
在进一步的实施例中，针对表3中的1个实时情景开展滞蓄关系预测和优化，具体如下：
[0130]
1)通过水文遥测设备获取水位站实时水位数据，基于获取的河网水位、预报降雨等实时数据，采用水文水动力模型模拟计算实时情景下的河网水位、河网槽蓄量、河网外排水量等要素，并计算实时情景下的区域蓄泄比sdr、圩内与圩外单位面积滞蓄水量比as2/as1、防洪风险指数r等指标值。
[0131]
2)判断实时情景防洪风险指数r值大小，若r＝0，则完成实时情景圩区-圩外系统滞蓄关系计算；若r＞0，则进一步判断圩内与圩外单位面积滞蓄水量比as2/as1与阈值c
as
大小，若as2/as1＞c
as
，则对实时情景圩区-圩外系统滞蓄关系进行优化。本实施例中实时情景初始水位为3.77m，属于前述第二类情景，模拟预测的防洪风险指数r＝0.97，as2/as1＝0.113，由于r＞0且as2/as1＜c
as,2
，需要进行实时情景圩区-圩外系统滞蓄关系优化。
[0132]
3)鉴于农业型圩区、非农业型圩区的保护标准、圩区自身排泄涝水能力、调蓄空间等有差异，采用系统聚类分析法，将区域内圩区分为若干类型。构建能够反映圩区自然属性与社会属性的指标，圩区自然属性主要包括圩内水面率、堤顶高程、排涝模数等指标；社会属性指标主要用于区分圩区农业圩或非农业圩属性，采用圩内耕地面积占比间接表征圩区社会属性，圩内耕地面积占比越大，表示该圩区农业圩属性越大。基于该方法，本实施例中圩区共分为5类，分别为a类、b类、c类、d类、其他。
[0133]
4)对于不同类型的圩区不同程度地增加圩区调蓄水深hs，本实施例中，a类圩区调蓄水深由0.2m增加至0.6m，b类圩区调蓄水深由0.2m增加至0.4m，c类圩区调蓄水深由0.2m增加至0.3m，d类圩区调蓄水深由0.1m增加至0.2m，其他圩区调蓄水深不变。圩区-圩外系统滞蓄关系优化后，重新预测该实时情景下的滞蓄关系和防洪风险指数，实时情景防洪风险
指数r的预测值由初始预测的0.97降低至优化后重新预测的0.78。优化前后结果对比见表3、图6。
[0134]
表3实时情景圩区-圩外系统滞蓄关系优化前后对比
[0135][0136]
实施例二：
[0137]
本实施例以二级嵌套圩区a(某个城市防洪工程)为例，进一步说明针对二级嵌套圩区如何进行水文水动力模型概化和圩区调蓄水深优化。
[0138]
该二级嵌套圩区a内部多级嵌套圩区拓扑关系图如图7a所示。对于包含二级小圩区的城市防洪工程，在进行水文水动力模型构建时，需要将该二级嵌套圩区进行“降级”处理，即通过完整概化一级圩区内的河道，将该二级嵌套圩区处理为多个二级小圩区，再将圩区a、若干二级小圩区与其余普通圩区一起进行聚类分析。
[0139]
本实施例中，圩区a内部包含r1、r2、r3、r4、r5五条河道，并分为a-1、a-2、a-3、a-4四个独立的二级小圩区。水文水动力模型构建时，分别概化r1、r2、r3、r4四条圩外河道和r1、r2、r3、r4、r5五条圩区内部河道，以及a-1、a-2、a-3、a-4四个二级小圩区，由此完成对二级嵌套圩区的“降级”。再将圩区a和a-1、a-2、a-3、a-4二级小圩区与其余普通圩区一起进行聚类分析。在完成分类的基础上，对圩区a和a-1、a-2、a-3、a-4二级小圩区进行调蓄水深hs优化。圩区排水时，a-1、a-2、a-3、a-4四个二级小圩区内实际滞蓄水深分别超过给定的调蓄水深值时，首先通过圩区排水设施(泵站)将各自内部涝水排入圩区a内部河道r1、r2、r3、r4、r5，或直接排入r1、r2、r3、r4四条圩区外部河道；由于二级小圩区的排水，圩区a内部河道r1、r2、r3、r4、r5的实际调蓄水深将增加，当圩区a内部河道r1、r2、r3、r4、r5的实际滞蓄水深超过给定的圩区a调蓄水深值时，圩区a通过排水设施(泵站)将内部涝水排入r1、r2、r3、r4四条圩外河道。
[0140]
在部分平原河网区域，存在一些较大的圩区，例如城市或农业区域，例如图4所述的大包围。在大包围中，常常存在若干小圩区，有的小圩区如图7b所示，存在有大圩区相同的边，在这种情况下，可以近似等效为一个大圩区分割成不同区域，需要对不同区域进行独立调控，即通过对大圩区的分割，对其降级，使分割后的小圩区与其它圩区在同一调控级别。
[0141]
但是也存在下述情况。如图7b所示，当一个大圩区内具有至少一个与大圩区不存在邻边的小圩区时(或者说小圩区被大圩区完全包覆，形成于大圩区内部时)，即可形成至少二重嵌套圩区，此时，内部圩区完全被外部圩区包围，且内部圩区和外部圩区之间具有河道等区域。例如图7b中圩1和下方的圩3被大圩区全部包围。在这种情况下，如果大圩区和被其包覆的小圩区采用同样的调控策略，有可能导致调控失效。对于小圩区而言，其外部的区域是大圩区的内部，因此，两者如果采用同样的调控策略，有可能导致大圩区内部、小圩区外部的区域被两种不同的调控策略管控，导致调控优化失败。
[0142]
总之，图7a和图7b的调控优化过程如下，在上述实施例各步骤中增加相关处理方法，为了简化表述，省略了上述方法的其他内容：
[0143]
构建预研究区域的圩区与圩区、圩区与河网的拓扑关系；
[0144]
基于上述拓扑关系，查找独立圩区、嵌套圩区与可分割圩区；
[0145]
将可分割圩区拆分成等效为独立圩区的若干个子圩区，并将所述子圩区插入所述独立圩区形成的集合；
[0146]
查找所述嵌套圩区的层级，并基于层级关系，标定内部圩区为二级圩区、三级圩区至n级圩区；基于已标定的圩区层级划定圩区调控的优先级。
[0147]
在进一步的实施例中，判断是否属于可分割圩区的方法是，内部圩区与外部圩区共边(圩堰)，且内部圩区包含于外部圩区；判断是否属于嵌套圩区的方法是，外部圩区全部包含内部圩区，且内部圩区与外部圩区没有重叠的圩堰；判断两个圩区为独立圩区的方法是，两个圩区互不包含，可共用部分圩堰。
[0148]
在进一步的实施例中，在实时情景滞蓄优化时，按照圩区的层级，从内之外，逐级优化；圩内的面积为圩堰包含的所有面积。
[0149]
因此当一个大圩区(外部圩区，下同)包含小圩区(内部圩区，下同)时，在计算大圩区圩内面积时，被包含的小圩区的面积也统计在其中。
[0150]
在这一实施例中，通过对圩区内外空间关系的构建，细化对圩区与圩区、圩区与河道之间的关系的描述，并采用上述方法进行调控。
[0151]
也就是说，在步骤s1或者s2时，构建预研究区域的空间拓扑关系，包括圩区与圩区、圩区与河网(河道)的拓扑关系；
[0152]
并基于这种拓扑关系，对圩区进行分类，能够分割的圩区，进行分割，形成平行同级的圩区，对于嵌套的圩区，分成不同层级；
[0153]
在模拟预测实时情景时，按照已经分号的层级设定调控的优先级。
[0154]
最后在调控的时候，根据优先级进行调控。
[0155]
实施例三：
[0156]
提出一种平原水网圩区-圩外系统滞蓄关系优化系统，用于实现平原水网区滞蓄关系优化方法，该系统包括实时数据获取模块、水文水动力模型、滞蓄关系计算模块、滞蓄关系优化模块，如图2所示，各模块功能具体如下：
[0157]
(1)实时数据获取模块：包括水文遥测系统和在线数据库，其中，水文遥测系统用于对实际情景中各水位站的水位数据进行实时感知和接收，水位站密度不低于1个/300km2，水位数据采集频次为10min/次，在线数据库用于存储历史水文数据、水文遥测系统的实时水文数据、预报降雨数据。
[0158]
(2)水文水动力模型：基于历史水文数据进行模拟计算，或基于实时获取的水位数据、预报的降雨数据进行模拟预测。本实施例中，共包括产流计算、汇流计算两部分，其中，产流计算时区分圩外河网、圩区内部河网，并分别按水面、水田、旱地、建设用地四类下垫面进行产流计算；汇流计算时，为能更精细、准确模拟坡面汇流过程和汇流区域范围，陆域坡面汇流模拟采用分布式汇流单位线方法进行计算，圩区调度通过设置圩区内部水面滞蓄水深、泵站排涝能力等参数实现，当圩区内部某时刻实际调蓄水深超过给定的滞蓄水深hs，圩区就开启排水设施(泵站)进行排水。
[0159]
(3)滞蓄关系计算模块：基于历史水文数据完成滞蓄关系特征指标计算，以及区域蓄泄比sdr～水文要素的非线性关系建立、相关指标阈值计算等。
[0160]
(4)滞蓄关系优化模块：基于实时水文数据，模拟预测实时情景下的区域滞蓄关系；根据防洪风险指数、预设的圩内与圩外单位面积滞蓄水量比阈值，判断是否需要优化滞蓄关系，针对判断结果为“是”的情况，根据步骤5和6所述的方法完成圩区聚类、增加调蓄水深、洪涝水时空分配优化等，从而实现圩区-圩外系统滞蓄关系优化，降低实时情景防洪风险指数。
[0161]
总之，本技术针对平原水网区城市防洪工程、圩区众多的现象，通过优化区域内洪涝水在圩内、圩外两个系统的时空分布，从而达到降低区域整体防洪风险目的的一种技术。其原则是基于不同片区河网、圩区、城市防洪工程的调蓄能力和洪涝承受能力，通过提前预降水位为后期滞蓄涝水预留空间、优化圩区调蓄水深等方式，合理挖掘河网调蓄水的潜力，优化区域洪涝水在时间尺度及不同对象中的空间分配，充分发挥区域河网、圩区对于洪涝水的滞蓄作用。
[0162]
如上所述，尽管参照特定的实施例已经表示和表述了本发明，但其不得解释为对本发明自身的限制。在不脱离所附权利要求定义的本发明的精神和范围前提下，可对其在形式上和细节上做出各种变化。

技术特征：
1.一种平原水网圩区-圩外系统滞蓄关系优化方法，其特征在于，包括如下步骤：s1、采集研究区域的历史水文和降雨数据，并基于所述历史水文和降雨数据识别场次降雨涨水期，建立典型历史情景集合并分类，基于河网-圩区拓扑关系，建立水文水动力模型，通过水文水动力模型计算模拟各个典型历史情景下河网水位、河网外排水量，以及各时刻下圩内与圩外河网的槽蓄量；s2、构建平原水网区蓄泄特征的表征指标，所述表征指标包含典型历史情景下的区域蓄泄比sdr
j
、区域蓄泄比sdr、圩外河网单位面积滞蓄水量as1、圩内河网单位面积滞蓄水量as2和防洪风险指数r；s3、基于典型历史情景下区域外排水量和槽蓄量的计算结果，计算典型历史情景j下的区域蓄泄比sdr
j
，分别绘制典型历史情景下区域蓄泄比sdr
j
与初始水位z
(j,1)
、累计降雨量p
j
的关系图，sdr
j
＝f(z
(j,1)
，p
j
)；并进行区域蓄泄比sdr与上述水文要素的相关性检验，通过相关性检验后，以各个(z
(j,1)
，ln(sdr
j
))为参数建立区域蓄泄比sdr～水文要素的关系，建立蓄泄比sdr与水文要素的关系式；s4、确定圩内、圩外单位面积滞蓄水量比的阈值c
as
；计算各典型历史情景下圩外单位面积滞蓄水量as1、圩内单位面积滞蓄水量as2、圩内与圩外单位面积滞蓄水量比as2/as1，并进一步确定圩内与圩外单位面积滞蓄水量比as2/as1的阈值c
as
；s5、模拟预测实时情景圩区-圩外系统滞蓄关系，并根据预测计算的实时情景防洪风险指数r，区域蓄泄比sdr～水文要素的关系，圩内与圩外单位面积滞蓄水量比as2/as1与阈值c
as
的比较，判断是否需要优化实时情景下的圩区-圩外系统滞蓄关系；s6、针对防洪风险指数r>0同时圩内与圩外单位面积滞蓄水量比as2/as1＜阈值c
as
的实时情景，采用系统聚类方法对圩区进行聚类分析，再通过增加不同类型圩区调蓄水深h
s
的方式来优化圩区-圩外系统滞蓄关系，以降低防洪风险指数r。2.根据权利要求1所述的一种平原水网圩区-圩外系统滞蓄关系优化方法，其特征在于，所述阈值c
as
的确定方法如下：在聚类分析后生成的各类典型历史情景中，取满足防洪风险指数r＝0的若干典型历史情景中圩内与圩外单位面积滞蓄水量比as2/as1的均值作为该类情景下的c
as
，式中，k为表示情景类别的自然数，c
as,k
为第k类情景对应的圩内与圩外单位面积滞蓄水量比阈值，j为表示情景序号的自然数，as
1,k,j
为属于k类情景下的情景j的圩外河网单位面积滞蓄水量，as
2,k,j
为属于k类情景下的情景j的圩内河网单位面积滞蓄水量，r
k,j
为k类情景中情景j的防洪风险指数，d
k
为k类情景中满足防洪风险指数r
k,j
＝0的情景个数。3.根据权利要求1所述的一种平原水网圩区-圩外系统滞蓄关系优化方法，其特征在于，s1中所述识别场次降雨涨水期，建立典型历史情景集合的过程，具体为：输入历史降雨、水位数据，场次降水开始的时刻记为t1，
t1时刻水位站i的水位记为z
1,i
，区域河网的平均水位记为z1，z1取各z
1,i
的算术平均值，场次降水中河网水位达到峰值的时刻记为t2，t2时刻对应的水位为峰值水位z
max
，从t1到t2的降水过程定义为一个场次降水事件的涨水期，t1、t2满足以下条件：t1之前特定时段
△
t内平均降雨强度小于1mm/d；t1～t2时段内单日最大降雨量p
1d
≥10mm；由此建立河网初始水位z1、场次降水涨水期降雨量p的二维情景集合{z1，p}。4.根据权利要求1所述的一种平原水网圩区-圩外系统滞蓄关系优化方法，其特征在于，步骤s1所述对建立的情景集合进行分类，具体为：采用系统聚类方法对构建的典型历史情景集合进行分类，聚类样本为构建的若干历史情景，聚类指标为初始水位均值z1。5.根据权利要求1所述的一种平原水网圩区-圩外系统滞蓄关系优化方法，其特征在于，步骤s1所述水文水动力模型构建及计算，具体为：收集平原水网区地形、下垫面、水系、圩区、水文站点资料；构建预研究区域的圩区与圩区、圩区与河网的拓扑关系；基于上述拓扑关系，查找独立圩区、嵌套圩区与可分割圩区；将可分割圩区拆分成等效为独立圩区的若干个子圩区，并将所述子圩区插入所述独立圩区形成的集合；查找所述嵌套圩区的层级，并基于层级关系，标定内部圩区为二级圩区、三级圩区至n级圩区；判断是否属于可分割圩区的方法是，内部圩区与外部圩区共边(圩堰)，且内部圩区包含于外部圩区；判断是否属于嵌套圩区的方法是，外部圩区全部包含内部圩区，且内部圩区与外部圩区没有重叠的圩堰；判断两个圩区为独立圩区的方法是，两个圩区互不包含，可共用部分圩堰；构建包括产汇流子模型、河网水动力子模型的水文水动力耦合模型，并采用水文站历史水位数据对模型进行率定验证；采用经率定的水文水动力模型模拟计算典型历史情景下的平原水网区的河网水位、河网槽蓄量和河网外排水量。6.根据权利要求1所述的一种平原水网圩区-圩外系统滞蓄关系优化方法，其特征在于，步骤s2中所述构建平原水网区蓄泄特征表征指标中，所述区域蓄泄比sdr
j
为sdr
j
＝s
j
/w
(j,out)
式中，sdr
j
为典型历史场景j下的区域蓄泄比；s
j
为典型历史场景j下的区域总滞蓄水量，包含圩外河网、圩内河网的滞蓄水量；w
(j,out)
为典型历史场景j下的区域河网外排水量；所述圩外河网单位面积滞蓄水量as
j1
为as
j1
＝s
j1
/u1；所述圩内河网单位面积滞蓄水量as
j2
为as
j2
＝s
j2
/u2；式中，as
j1
为历史情景j下的圩外河网单位面积滞蓄水量，as
j2
为历史情景j下的圩内河网单位面积滞蓄水量，u1为圩外河网面积，u2为圩区面积，s
j1
为历史情景j下的圩外河网滞蓄水量，s
j2
为历史情景j下的圩内河网滞蓄水量；
为终止时刻t2时圩外的槽蓄量，为初始时刻t1时圩外的槽蓄量，为终止时刻t2时圩内的槽蓄量，为初始时刻t1时圩内的槽蓄量；其余变量含义同前；所述防洪风险指数r为：所述防洪风险指数r为：所述防洪风险指数r为：式中：r
i
为水位站i的防洪风险指数，r为区域防洪风险指数，z
i
(t)为水位站i的水位过程，h
i
为水位站i的保证水位，n为站点数量；其余变量含义同前；当r>0时，表示区域内部分或全部水位站水位超过保证水位，区域存在一定防洪风险。7.根据权利要求5所述的一种平原水网圩区-圩外系统滞蓄关系优化方法，其特征在于，所述步骤s5进一步包括：s51、实时情景圩区-圩外系统滞蓄关系模拟预测：基于获取的包括河网水位、预报降雨在内的实时数据，采用所述水文水动力模型模拟预测实时情景下的河网水位、河网槽蓄量和河网外排水量，并计算得到实时情景下的区域蓄泄比sdr，圩内与圩外单位面积滞蓄水量比as2/as1以及防洪风险指数r；s52、实时情景圩区-圩外系统滞蓄关系判断：首先判断实时情景防洪风险指数r值大小，若防洪风险指数r＝0，则完成实时情景圩区-圩外系统滞蓄关系计算；若防洪风险指数r＞0，则进一步判断圩内与圩外单位面积滞蓄水量比as2/as1与阈值c
as
大小，若圩内与圩外单位面积滞蓄水量比as2/as1＜阈值c
as
，则对实时情景圩区-圩外系统滞蓄关系进行优化，否则完成实时情景圩区-圩外系统滞蓄关系计算。8.根据权利要求1所述的一种平原水网圩区-圩外系统滞蓄关系优化方法，其特征在于，所述s6中所述圩区聚类分析方法具体为：构建反映圩区自然属性与社会属性的指标，圩区自然属性包括圩内水面率、堤顶高程和排涝模数，其中圩内水面率表征圩区自然状态下调蓄能力，堤顶高程差异间接反映圩区地面高程差异，排涝模数表征工程作用下圩区调蓄能力；社会属性指标用于区分圩区农业圩或非农业圩属性，采用圩内耕地面积占比间接表征圩区社会属性，圩内耕地面积占比越大，表示该圩区农业圩属性越大；基于上述自然属性和社会属性指标，采用系统聚类分析法将区域内圩区分为若干类型。9.根据权利要求1所述的一种平原水网圩区-圩外系统滞蓄关系优化方法，其特征在于，所述s6中所述实时情景圩区-圩外系统滞蓄关系优化，具体为：针对防洪风险指数r>0同时圩内与圩外单位面积滞蓄水量比as2/as1＜阈值c
as
的实时情景，在保证圩区自身安全的前提下，对于不同类型的圩区不同程度地增加圩区调蓄水深h
s
，以增加涨水期圩内河网滞蓄水量，优化洪涝水在圩外河网、圩内河网的时空分配，从而实现
圩区-圩外系统滞蓄关系优化，降低实时情景防洪风险指数。10.根据权利要求1所述的一种平原水网圩区-圩外系统滞蓄关系优化方法，其特征在于，所述步骤s5还包括：基于权利要求5步骤s1中已标定的圩区层级划定圩区调控的优先级；所述步骤s6还包括：在实时情景滞蓄优化时，按照圩区的层级，从内之外，逐级优化；圩内的面积为圩堰包含的所有面积。11.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至10中任一项所述方法的步骤。

技术总结
本发明公开了一种平原水网圩区-圩外系统滞蓄关系优化方法及系统，该方法包括：根据历史水文、降雨数据识别场次降雨涨水期，建立典型历史情景集合并分类，模拟典型历史情景下河网水位、河网外排水量和圩内与圩外河网的槽蓄量；构建平原水网区蓄泄特征表征指标；建立区域蓄泄特征指标~水文要素的非线性关系；确定圩内、圩外单位面积滞蓄水量比的阈值；模拟预测实时情景圩区-圩外系统滞蓄关系，并判断是否需要优化实时情景滞蓄关系；针对有优化需要的实时情景，完成实时情景圩区-圩外系统滞蓄关系优化。本发明为自动模拟预测实时情景下水网区洪涝水的时空分配，判断并优化实时情景下水网区的滞蓄关系提供了可操作的方法和智能化系统，可指导水网区实时防洪调度。可指导水网区实时防洪调度。可指导水网区实时防洪调度。


技术研发人员：龚李莉 李琛 王元元 陆志华 钱旭 白君瑞 张怡 李勇涛 潘明祥 韦婷婷 蔡梅 李蓓
受保护的技术使用者：太湖流域管理局水利发展研究中心
技术研发日：2022.02.14
技术公布日：2022/5/25