CN112101818A - 一种适用于复杂水力联系的海绵城市洪涝优化调度方法 - Google Patents
一种适用于复杂水力联系的海绵城市洪涝优化调度方法 Download PDFInfo
- Publication number
- CN112101818A CN112101818A CN202011092833.7A CN202011092833A CN112101818A CN 112101818 A CN112101818 A CN 112101818A CN 202011092833 A CN202011092833 A CN 202011092833A CN 112101818 A CN112101818 A CN 112101818A
- Authority
- CN
- China
- Prior art keywords
- flood
- scheduling
- flow
- water
- urban
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 69
- 238000005457 optimization Methods 0.000 title claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims abstract description 22
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 171
- 238000003860 storage Methods 0.000 claims description 39
- 230000008569 process Effects 0.000 claims description 21
- 238000010586 diagram Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 238000005086 pumping Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 201000004569 Blindness Diseases 0.000 abstract description 3
- 244000207740 Lemna minor Species 0.000 description 15
- 235000006439 Lemna minor Nutrition 0.000 description 15
- 235000001855 Portulaca oleracea Nutrition 0.000 description 15
- 241000272814 Anser sp. Species 0.000 description 12
- 239000010977 jade Substances 0.000 description 8
- 238000011144 upstream manufacturing Methods 0.000 description 6
- 230000009467 reduction Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 4
- 230000001276 controlling effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 230000001105 regulatory effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 235000018996 Nuphar luteum Nutrition 0.000 description 1
- 244000209105 Nuphar luteum Species 0.000 description 1
- 241001536526 Nuphar pumila Species 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000008239 natural water Substances 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Geometry (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Operations Research (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Flow Control (AREA)
Abstract
本发明公开了一种适用于复杂水力联系的海绵城市洪涝优化调度方法,包括:构建城区水力联系拓扑关系图;确定洪涝调度的主要节点、优化计算时段间隔;根据洪涝调度的主要控制目标,确定洪涝调度的目标函数;确定洪涝调度的约束条件;根据洪涝调度的目标函数和约束条件,建立洪涝调度模型;按照洪涝调度的目标函数优先级,采用分层优化法,结合洪涝调度的优化计算时段间隔,优化求解洪涝调度模型,集成洪涝调度系统,并根据洪涝调度模型的优化调度结果,指导城区洪涝调度。本发明针对具有复杂网状水力联系的海绵城市洪涝调度问题,可避免现有经验调度方式的盲目性,实现城市洪涝问题的科学和精细化调度,消除或最大程度减轻城市洪涝灾害损失。
Description
技术领域
本发明涉及城市防洪减灾技术领域,更具体的涉及一种适用于复杂网状水力联系的海绵城市洪涝优化调度方法。
背景技术
城市防洪减灾是海绵城市建设的主要目的之一,海绵城市建设旨在通过渗、滞、蓄、净、用、排等多种技术,维持或恢复城市的“海绵”功能,实现城市良性水文循环。从非工程措施角度分析,为系统而又高效地发挥工程对洪水的削减、转输、调蓄等功能,达到有效消除或减轻城市雨洪灾害目的,开展城市洪涝优化调度工作,形成科学的调度决策和方案,这些非工程措施是海绵城市建设不可回避的重要工作。
现有的城市洪涝灾害研究和防洪减灾方法主要集中在城市洪涝过程模拟和洪涝预报,城市洪涝调度方面则存在欠缺,以经验调度为主,优化调度鲜有涉及。海绵城市建设重构了城市下垫面、水系空间,输调配水工程以及河湖连通工程使得城市水系统拓扑结构从单纯的树状结构变为复杂的网状结构。与其他区域水资源优化调度相比,海绵城市雨水系统优化调度具有较多特殊性,比如水资源系统环状网络结构、洪水顶托、天然水道和人工渠化水道并存的水量转输和调蓄等问题。经验调度方式和传统的水资源调度方法不能考虑复杂网状水力联系,湖泊蓄泄、控制闸操作和泵排流量决策的时空组合极多,各约束耦合关系复杂,经验调度方式无法做到全域、全时空的统一调度,无法发挥海绵城市工程的最佳效用,因此不能适用于复杂的海绵城市洪涝调度问题。基于此,本发明提出一种适用于复杂网状水力联系的海绵城市洪涝优化调度方法。
发明内容
本发明实施例提供一种适用于复杂水力联系的海绵城市洪涝优化调度方法,用以解决上述背景技术中提出的问题。
本发明实施例提供一种适用于复杂水力联系的海绵城市洪涝优化调度方法,包括:
分析海绵城市河湖管网系统,构建城区水力联系拓扑关系图;
根据城区水力联系拓扑关系图中的节点,确定洪涝调度的主要节点;并根据主要节点之间的水流传播时间,确定洪涝调度的优化计算时段间隔;
根据洪涝调度的主要控制目标,确定洪涝调度的目标函数;
根据河段流量演算、湖泊水库蓄水限制和水量平衡、关键节点水位和流量控制、河道水位-流量-断面面积关系、控制闸过流能力和相邻时段流量变幅,确定洪涝调度的约束条件;
根据洪涝调度的目标函数和洪涝调度的约束条件,建立洪涝调度模型;
按照洪涝调度的目标函数优先级,采用分层优化法,结合洪涝调度的优化计算时段间隔,优化求解洪涝调度模型;
集成洪涝调度系统,并根据洪涝调度模型的优化调度结果,指导城区洪涝调度。
进一步地,所述根据海绵城市河湖管网系统,构建城区水力联系拓扑关系图,具体包括:
实地考察海绵城市河湖管网系统,结合海绵城市规划设计资料和卫星地图,依据海绵城区防洪排涝主次关系,构建城区水力联系拓扑关系图;
其中,所述海绵城市河湖管网系统包括:河、湖、水库、连通渠、控制闸、排涝泵站、排水管网。
进一步地,所述城区水力联系拓扑关系图中的节点,包括:
关键水位控制点、水文/水位站点、河渠控制闸、排涝泵站。
进一步地,所述洪涝调度的目标函数,包括:
洪涝灾害损失最小,用调度洪峰流量最小、关键水位控制点淹没时长表示;
调度成本最低,用泵排总流量、控制闸操作次数表示。
进一步地,所述洪涝调度的约束条件,具体包括:
考虑水流滞时的城区湖泊(水库)水量平衡,表示为:
考虑水流滞时的城区河道水量平衡,表示为:
其中,i表示湖泊,j表示河道,P表示区间汇流量,Q表示入流量,O表示出流量,A表示流入控制闸集合,B表示流出控制闸集合,m表示控制闸,τm表示控制闸相应渠道的水流传播时间,V表示蓄水量,t表示时段;
湖泊(水库)蓄水限制,表示为:
其中,Vmin、Vmax分别是湖泊允许的最小和最大蓄水量,此处用蓄水量表示与水量平衡一致,避免水位-蓄水量关系的换算,优化计算完成之后再将蓄水量换算为水位;
控制闸过流能力限制,表示为:
控制闸相邻时段流量变幅应尽可能平稳,表示为:
其中,Δqk是控制闸k在相邻时段的流量变幅;
划分的河段节点处水位-流量-断面面积关系由水文站点数据拟合,缺少基本数据时根据邻近水文站点和水力学公式推算;水位-流量-断面面积关系采用二次表达式拟合,不宜采用高次多项式,以免规划求解中出现数值病态;若曲线关系的线性相关度高,采用线性表达式拟合,同时降低计算复杂度;
河道断面面积-流量关系,表示为:
河道水位-流量关系,表示为:
其中,Q表示流量,S表示断面面积,Z表示水位;
洪涝调度模型中各变量均取非负值。
进一步地,所述洪涝调度的约束条件,具体还包括:
关键节点水位和流量控制、水位顶托关系、调度期湖泊水库初末蓄水量或水位。
本发明实施例提供的适用于复杂水力联系的海绵城市洪涝优化调度方法,还包括:对洪涝调度模型中特殊形式的约束进行处理,使之满足数学规划的要求;
对流量变幅控制涉及的绝对值约束式(5),进行等价处理如下:
进一步地,所述优化求解洪涝调度模型,具体包括:
求解优先级最高的目标,得到优化调度的最大洪峰流量;
以最大洪峰流量为约束,即在不超过最大洪峰流量的前提下,对泵排流量和控制闸流量变幅进行优化;
在分层优化求解过程中,两级优化均是非线性规划,采用基于Python的优化建模工具Pyomo编程实现洪涝调度模型,并且调用Ipopt采用内点法进行求解。
进一步地,所述集成洪涝调度系统,并根据洪涝调度模型的优化调度结果,指导城区洪涝调度,具体包括:
将优化后的洪涝调度模型集成到洪涝调度系统中,通过txt文本交换数据,输入是预测的城区洪水过程,输出是经过优化计算得到的湖泊水库蓄泄过程和控制闸流量过程;洪涝调度系统接受调度计算命令后,在后台执行调度模型程序,然后在洪涝调度系统界面上显示最优的调度决策,根据最优的调度决策指导城区洪涝调度。
本发明实施例提供一种适用于复杂水力联系的海绵城市洪涝优化调度方法,与现有技术相比,其有益效果如下:
本发明针对具有复杂网状水力联系的海绵城市洪涝调度问题,通过水资源系统分析方法,运用数学规划理论,构建科学有效的优化调度模型求解方法,为城市洪涝调度和防洪减灾提供科学的方案和决策;该洪涝优化调度方法可避免现有经验调度方式的盲目性,实现城市洪涝问题的科学和精细化调度,消除或最大程度减轻城市洪涝灾害损失;系统而又高效地发挥海绵城市工程对洪水的的削减、转输、调蓄等功能,形成科学的调度决策和方案,达到有效消除或减轻城市雨洪灾害的目的。
附图说明
图1为本发明实施例提供的海绵城市洪涝调度模型和方法的流程图;
图2为本发明实施例提供的萍乡城区河湖管网系统拓扑结构;
图3为本发明实施例提供的萍乡水文站历史实测和模拟调度的水位-流量对比关系图;
图4为本发明实施例提供的五丰水文站(五丰河口)历史实测和模拟调度的水位-流量对比关系图;
图5为本发明实施例提供的各湖泊蓄水过程;
图6a为本发明实施例提供的控制闸k=1流量过程;
图6b为本发明实施例提供的控制闸k=2流量过程;
图6c为本发明实施例提供的控制闸k=3流量过程;
图6d为本发明实施例提供的控制闸k=5流量过程;
图6e为本发明实施例提供的控制闸k=6流量过程;
图6f为本发明实施例提供的控制闸k=7流量过程。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参见图1,本发明实施例提供一种适用于复杂水力联系的海绵城市洪涝优化调度方法,该方法包括:
步骤1,分析海绵城市河湖管网系统水力关系和结构,构建水力联系拓扑关系图。
步骤2,根据研究区域的河湖水力联系、关键水位控制点、水文/水位站点、河渠控制闸和排涝泵站等,确定洪涝调度的主要节点,考虑主要节点之间的水流传播时间,确定调度计算的时段间隔;
步骤3,根据洪涝调度的主要控制目标,确定洪涝调度的目标函数;
步骤4,根据河段流量演算、湖泊(水库)蓄水限制和水量平衡、关键节点水位和流量控制、河道水位-流量-断面面积关系、控制闸过流能力和相邻时段流量变幅,确定洪涝调度的约束条件;
步骤5,根据洪涝调度的目标函数和洪涝调度的约束条件,建立洪涝优化调度模型,对模型中特殊形式的约束进行处理,使之满足数学规划的要求;
步骤6,求解并验证所述洪涝调度模型;
步骤7,基于所述洪涝调度模型的优化调度结果指导城区洪涝调度。
对上述步骤1~7的具体说明如下:
对于步骤1:海绵城市输调配水工程以及河湖连通工程使得城市水系统拓扑结构从单纯的树状结构变为复杂的网状结构,系统主要由河、湖(水库)、连通渠、控制闸、排涝泵站、排水管网等组成。通过实地考察,结合海绵城市规划设计资料和卫星地图,依据城区防洪排涝主次关系,构建城区水力联系拓扑关系图。
对于步骤2:海绵城市防涝调度主要针对场次暴雨,即以场次暴雨历时为调度期。优化调度计算采用的时间间隔应考虑模型的精度需求,以及时段间隔与关键节点之间水流传播时间的耦合关系。
对于步骤3:海绵城市洪涝调度常用的控制目标是洪涝灾害损失最小和调度成本最低,前者可用调度洪峰流量最小、关键水位控制点淹没时长等表示,后者可用泵排总流量、控制闸操作次数等表示。
对于步骤4:具体包括以下步骤:
步骤4.1,一般地,考虑水流滞时的城区湖泊(水库)水量平衡可表示为
考虑水流滞时的城区河道水量平衡可表示为
其中,i表示湖泊,j表示河道,P表示区间汇流量,Q表示入流量,O表示出流量,A表示流入控制闸集合,B表示流出控制闸集合,m表示控制闸,τm表示控制闸相应渠道的水流传播时间,V表示蓄水量,t表示时段。
步骤4.2,湖泊(水库)蓄水限制一般可表示为
其中Vmin、Vmax分别是湖泊允许的最小和最大蓄水量。此处用蓄水量表示与水量平衡一致,避免水位-蓄水量关系的换算,优化计算完成之后再将蓄水量换算为水位即可。
步骤4.3,控制闸过流能力限制一般可表示为
其中Δqk是控制闸k在相邻时段的流量变幅。
步骤4.4,划分的河段节点处水位-流量-断面面积关系由水文站点数据拟合,缺少基本数据时可根据邻近水文站点和水力学公式推算。水位-流量-断面面积关系采用二次表达式拟合,不宜采用高次多项式,以免规划求解中出现数值病态;若曲线关系的线性相关度很高,也可采用线性表达式拟合,同时降低计算复杂度。河道断面面积-流量关系一般可表示为
河道水位-流量关系一般可表示为
其中Q表示流量,S表示断面面积,Z表示水位。
步骤4.5,模型中各变量均取非负值,涉及的其他约束还包括:关键节点水位和流量控制,水位顶托关系,调度期湖泊(水库)初末蓄水量(或水位)等。
对于步骤5:对模型中的非标准约束进行处理,使其满足数学规划的要求。这里主要是流量变幅控制涉及的绝对值约束,式(5)可等价处理如下:
对于步骤6:具体包括以下步骤:
步骤6.1,调度模型一般涉及多个目标函数,按目标函数的优先级采用分层优化法求解。其中防洪效益最为重要,即先求解优先级最高的目标,得到优化调度的最大洪峰流量;然后以此为约束,即在不超过此洪峰流量的前提下,对泵排流量和控制闸流量变幅进行优化。
步骤6.2,在分层优化求解过程中,两级优化均是非线性规划,采用基于Python的优化建模工具Pyomo编程实施上述调度模型,调用Ipopt采用内点法进行求解。
步骤6.3,以城区典型洪水过程为输入,采用上述优化调度模型和求解方法进行模拟调度,分析计算结果和调度决策的合理性、可操作性,对调度模型和求解方法进行修正。
对于步骤7:以预测的城区洪水过程为输入,运用上述调度模型和方法,得到城区河湖各控制闸和排涝泵站的调度方案。在实际运用中,可根据城区水雨情监测和预报,每隔若干时段滚动更新,不断生成新的调度决策,进一步提高调度的科学性。
实施例
江西省萍乡市是全国第一批海绵城市建设试点城市之一,城区河湖连通工程情况如下。萍水河田中水文站以上原弯曲河段经裁弯取直之后形成萍水湖,河道行洪能力增加,并且可以配合运用萍水湖的调蓄能力。五丰河上游修建赤山隧洞与黄土开水库下游河道连通,隧洞两侧布设多向闸门,通过闸门的组合运用控制水流方向。枯水期可引黄土开水库水流至玉湖,增加景观效益,洪水期可拦截五丰河上游流量。玉湖位于五丰河中游,出口设闸门控制。五丰河中下游于万龙湾建设路桥下通过箱涵引水至鹅湖,减轻万龙湾等内涝区的防洪压力,鹅湖通过泵站排水至萍水河。五丰河口为防止萍水顶托,设闸门和泵站,必要时关闸排水。
萍乡市洪涝灾害频发,2019年7月萍乡市出现大面积洪涝灾害,其固有缺乏大型调蓄工程的缘故,但也未充分利用现有条件进行合理调度,以减轻其灾害程度。海绵城市建设后,城区河湖关系更为密切,系统更为复杂,亟需科学的调度模型和方法。实施例建立萍乡城区洪涝调度模型,针对2019年7月洪水进行调度,对比分析调度的防洪效益。
(1)河湖系统分析
萍乡市河湖管网雨水系统拓扑关系如图2所示,其中圆圈表示湖泊或水库,环形圆圈表示赤山隧洞,线条表示河流或渠道,箭头表示水流方向。
控制闸及其渠道视为整体,可调控流量,主要属性是最大过流能力。控制闸在图2中用符号k表示,详见表1。
表1控制闸及其渠道
控制闸k | 位置及渠道 | 最大过流能力(m/s) |
1 | 萍水湖进水口 | 100 |
2 | 萍水湖出水口 | 100 |
3 | 经赤山隧洞引水至萍水 | 17.5 |
4 | 经赤山隧洞供水至玉湖 | 17.5 |
5 | 玉湖出水口 | 20 |
6 | 鹅湖进水口(箱涵) | 20 |
7 | 鹅湖出水口 | 51.4(泵) |
8 | 五丰河注入萍水 | 15(泵) |
五丰河和萍水河是主要研究对象,依据水力关系、湖泊和控制闸划分节点。五丰河和萍水河的节点位置及沿河距离分别见表2和表3。水流传播时间如下:赤山隧洞经上游支流调水至萍水1h;五丰河1-3河段1h,3-5河段1h;萍水河1-3河段1h,3-7河段1h。
表2五丰河主要节点
序号m | 节点位置 | 沿河距离(km) |
1 | 赤山隧洞 | 0 |
2 | 玉湖入口 | 4.52 |
3 | 玉湖出口 | 5.98 |
4 | 鹅湖入口(箱涵、万龙湾) | 7.85 |
5 | 河口(注入萍水) | 10.01 |
表3萍水河节点
序号n | 节点位置 | 沿河距离(km) |
1 | 上游支流汇入 | 0 |
2 | 萍水湖入口 | 8.34 |
3 | 萍水湖出口(福田河汇入) | 9.2 |
4 | 鹅湖出口 | 13.26 |
5 | 五丰河汇入 | 15.06 |
6 | 白源河汇入 | 15.43 |
7 | 萍乡水文站 | 20 |
(2)变量定义
Vi min-湖泊最小蓄水量限制
Vi max-湖泊最大蓄水量限制
②河道(五丰河、萍水河)涉及的主要变量如下:
③控制闸及其渠道涉及的主要变量如下:
(3)目标函数
调度期(T):场次暴雨
时段间隔(t):小时
②抽排水费用最小,直接用抽水量(鹅湖泵站)表示
③控制闸流量应尽可能平稳,即变幅最小
(4)约束条件
①湖泊蓄水量限制
②五丰河水量关系
玉湖(i=2)水量平衡
五丰河3-5河段水量平衡(河段蓄水量=河段上下断面平均面积×河段长度)
经箱涵引水至鹅湖
鹅湖(i=3)水量平衡
③萍水河水量关系
萍水河1-3河段水量平衡
萍水湖(i=1)水量平衡
萍水河3-7河段水量平衡
萍水河下断面(萍乡水文站)最大流量限制
④河湖节点水位关系
萍水湖下泄时,萍水湖出口水位应高于萍水河水位
鹅湖泵站工作外水位限制
⑤河道水位-流量-断面面积关系
五丰水文站流量-断面面积关系(多项式系数为a)
萍水河上断面流量-断面面积关系(多项式系数为b)
田中水文站流量-断面面积关系(多项式系数为c)
萍乡水文站流量-断面面积关系(多项式系数为d)
田中水文站水位-流量关系(多项式系数为e)
萍乡水文站水位-流量关系(多项式系数为f)
⑥控制闸约束
控制闸最大过流能力限制
控制闸相邻时段流量变幅控制
⑦决策变量约束
模型中的各变量均取非负值。
(5)求解方法
上述洪水调度模型涉及大量的等式和不等式约束,其中河道的水位-流量-断面面积曲线是非线性关系,这里采用二次多项式表示这些曲线关系(高次多项式的微小系数可能导致数值病态)。因此,该模型是一个大规模的非线性规划问题。
模型涉及三个目标函数,按目标函数的优先级采用分层优化法求解。其中防洪效益最为重要,即先求解优先级最高的目标1,得到优化调度的最大洪峰流量;然后以此为约束,即在不超过此洪峰流量的前提下,对调度费用和控制闸流量及其变幅进行优化,这里直接把目标2和3相加作为目标函数进行求解。
(6)模拟调度分析
2019年7月萍乡城区洪水从8日2:00至10日20:00,共67个时段。萍乡水文站实测和调度的水位-流量对比如图3所示,图中实线表示实测、虚线表示调度。2019年7月8日2:00,水文站实测水位84.18m,相应流量101.83m3/s;至7月8日22:00,水文站实测水位87.04m,相应流量488.78m3/s,已超过淹没高程87m,开始产生内涝;淹没时间23h,最大淹没深度出现在7月9日10:00—11:00(第33个时段),水位88.32m,最大内涝深度1.32m,相应流量725.42m3/s。优化调度后,9日0:00至22:00出现内涝,最高水位87.70m,最大内涝深度0.70m,内涝深度消减了0.62m,相应流量610m3/s,消减流量115.42m3/s。
五丰水文站实测和调度的水位-流量对比如图4所示,图中实线表示实测、虚线表示调度,。优化调度后,五丰水文站最高水位89.3m,最大流量控制在15m3/s,五丰河主城区不出现内涝。第23个时段(9日0:00—1:00),五丰河出口处萍水河水位(89.37m)高于五丰河口水位(89.30m),开始关闸启动泵排,一直维持到第46个时段(9日23:00—24:00),共24h,泵排流量15m3/s,泵排时长23h。
由此可见,上述调度模型和方法效果明显,有效解决了五丰河主城区的洪涝问题,同时大幅度消减城区出口断面(萍乡水文站)洪峰流量。城区各湖泊调度蓄水过程见图5,各控制闸调度流量过程见图6a~6f。
本发明的洪涝调度模型和方法针对具有复杂网状水力联系的海绵城市洪涝问题,能够充分运用河湖的调蓄能力,河渠控制闸和排涝泵站协调运作,实现了全域、全时空的统一优化调度,得到了洪涝期间各时段各湖泊的蓄泄决策和各控制闸的过流量决策。该优化调度模型和方法可避免现有经验调度方式的盲目性,实现城市洪涝问题的科学和精细化调度,消除或最大程度减轻城市洪涝灾害损失。
以上公开的仅为本发明的几个具体实施例,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围,但是,本发明实施例并非局限于此,任何本领域的技术人员能思之的变化都应落入本发明的保护范围内。
Claims (9)
1.一种适用于复杂水力联系的海绵城市洪涝优化调度方法,其特征在于,包括:
分析海绵城市河湖管网系统,构建城区水力联系拓扑关系图;
根据城区水力联系拓扑关系图中的节点,确定洪涝调度的主要节点;并根据主要节点之间的水流传播时间,确定洪涝调度的优化计算时段间隔;
根据洪涝调度的主要控制目标,确定洪涝调度的目标函数;
根据河段流量演算、湖泊水库蓄水限制和水量平衡、关键节点水位和流量控制、河道水位-流量-断面面积关系、控制闸过流能力和相邻时段流量变幅,确定洪涝调度的约束条件;
根据洪涝调度的目标函数和洪涝调度的约束条件,建立洪涝调度模型;
按照洪涝调度的目标函数优先级,采用分层优化法,结合洪涝调度的优化计算时段间隔,优化求解洪涝调度模型;
集成洪涝调度系统,并根据洪涝调度模型的优化调度结果,指导城区洪涝调度。
2.如权利要求1所述的适用于复杂水力联系的海绵城市洪涝优化调度方法,其特征在于,所述根据海绵城市河湖管网系统,构建城区水力联系拓扑关系图,具体包括:
实地考察海绵城市河湖管网系统,结合海绵城市规划设计资料和卫星地图,依据海绵城区防洪排涝主次关系,构建城区水力联系拓扑关系图;
其中,所述海绵城市河湖管网系统包括:河、湖、水库、连通渠、控制闸、排涝泵站、排水管网。
3.如权利要求1所述的适用于复杂水力联系的海绵城市洪涝优化调度方法,其特征在于,所述城区水力联系拓扑关系图中的节点,包括:
关键水位控制点、水文/水位站点、河渠控制闸、排涝泵站。
4.如权利要求1所述的适用于复杂水力联系的海绵城市洪涝优化调度方法,其特征在于,所述洪涝调度的目标函数,包括:
洪涝灾害损失最小,用调度洪峰流量最小、关键水位控制点淹没时长表示;
调度成本最低,用泵排总流量、控制闸操作次数表示。
5.如权利要求1所述的适用于复杂水力联系的海绵城市洪涝优化调度方法,其特征在于,所述洪涝调度的约束条件,具体包括:
考虑水流滞时的城区湖泊(水库)水量平衡,表示为:
考虑水流滞时的城区河道水量平衡,表示为:
其中,i表示湖泊,j表示河道,P表示区间汇流量,Q表示入流量,O表示出流量,A表示流入控制闸集合,B表示流出控制闸集合,m表示控制闸,τm表示控制闸相应渠道的水流传播时间,V表示蓄水量,t表示时段;
湖泊(水库)蓄水限制,表示为:
Vi min≤Vi t≤Vi max (3)
其中,Vmin、Vmax分别是湖泊允许的最小和最大蓄水量,此处用蓄水量表示与水量平衡一致,避免水位-蓄水量关系的换算,优化计算完成之后再将蓄水量换算为水位;
控制闸过流能力限制,表示为:
控制闸相邻时段流量变幅应尽可能平稳,表示为:
其中,Δqk是控制闸k在相邻时段的流量变幅;
划分的河段节点处水位-流量-断面面积关系由水文站点数据拟合,缺少基本数据时根据邻近水文站点和水力学公式推算;水位-流量-断面面积关系采用二次表达式拟合,不宜采用高次多项式,以免规划求解中出现数值病态;若曲线关系的线性相关度高,采用线性表达式拟合,同时降低计算复杂度;
河道断面面积-流量关系,表示为:
河道水位-流量关系,表示为:
其中,Q表示流量,S表示断面面积,Z表示水位;
洪涝调度模型中各变量均取非负值。
6.如权利要求5所述的适用于复杂水力联系的海绵城市洪涝优化调度方法,其特征在于,所述洪涝调度的约束条件,具体还包括:
关键节点水位和流量控制、水位顶托关系、调度期湖泊水库初末蓄水量或水位。
8.如权利要求1所述的适用于复杂水力联系的海绵城市洪涝优化调度方法,其特征在于,所述优化求解洪涝调度模型,具体包括:
求解优先级最高的目标,得到优化调度的最大洪峰流量;
以最大洪峰流量为约束,即在不超过最大洪峰流量的前提下,对泵排流量和控制闸流量变幅进行优化;
在分层优化求解过程中,两级优化均是非线性规划,采用基于Python的优化建模工具Pyomo编程实现洪涝调度模型,并且调用Ipopt采用内点法进行求解。
9.如权利要求1所述的适用于复杂水力联系的海绵城市洪涝优化调度方法,其特征在于,所述集成洪涝调度系统,并根据洪涝调度模型的优化调度结果,指导城区洪涝调度,具体包括:
将优化后的洪涝调度模型集成到洪涝调度系统中,通过txt文本交换数据,输入是预测的城区洪水过程,输出是经过优化计算得到的湖泊水库蓄泄过程和控制闸流量过程;洪涝调度系统接受调度计算命令后,在后台执行调度模型程序,然后在洪涝调度系统界面上显示最优的调度决策,根据最优的调度决策指导城区洪涝调度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011092833.7A CN112101818B (zh) | 2020-10-13 | 2020-10-13 | 一种适用于复杂水力联系的海绵城市洪涝优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011092833.7A CN112101818B (zh) | 2020-10-13 | 2020-10-13 | 一种适用于复杂水力联系的海绵城市洪涝优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112101818A true CN112101818A (zh) | 2020-12-18 |
CN112101818B CN112101818B (zh) | 2024-04-23 |
Family
ID=73783610
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011092833.7A Active CN112101818B (zh) | 2020-10-13 | 2020-10-13 | 一种适用于复杂水力联系的海绵城市洪涝优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112101818B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112950125A (zh) * | 2021-04-02 | 2021-06-11 | 拉扎斯网络科技(上海)有限公司 | 一种配送资源安全状态数据的获得方法、装置、电子设备 |
CN112989253A (zh) * | 2021-04-02 | 2021-06-18 | 南昌工程学院 | 一种基于JupyterHub的水库优化调度实践教学平台的搭建方法 |
CN115017665A (zh) * | 2022-06-29 | 2022-09-06 | 珠江水利委员会珠江水利科学研究院 | 用于城市暴雨洪涝防御的模拟方法、装置、设备及存储介质 |
CN117236673A (zh) * | 2023-11-16 | 2023-12-15 | 水利部交通运输部国家能源局南京水利科学研究院 | 城市水网多尺度防洪排涝联合优化调度方法和系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102296562A (zh) * | 2010-06-25 | 2011-12-28 | 华东电网有限公司 | 防洪与发电相耦合的梯级水库联合洪水调度优化方法 |
CN109284531A (zh) * | 2018-08-03 | 2019-01-29 | 中国水利水电科学研究院 | 一种基于空间拓扑的一二维水动力学耦合方法 |
CN109872063A (zh) * | 2019-02-11 | 2019-06-11 | 南昌工程学院 | 平原城市河网防洪、排涝、活水联合优化调度方法及系统 |
CN109919397A (zh) * | 2019-04-12 | 2019-06-21 | 黄河勘测规划设计研究院有限公司 | 平原河网地区排涝容量校验优化方法 |
CN110838079A (zh) * | 2019-10-18 | 2020-02-25 | 南京中科智慧生态科技有限公司 | 一种智慧城市防洪排水远程监控系统及方法 |
CN111046574A (zh) * | 2019-12-20 | 2020-04-21 | 武汉大学 | 一种平原湖区湖河闸泵系统防洪排涝演算方法 |
CN111274690A (zh) * | 2020-01-16 | 2020-06-12 | 无锡北大建筑工程有限公司 | 一种海绵城市的雨水源头调蓄系统的构建方法 |
CN111476676A (zh) * | 2020-04-02 | 2020-07-31 | 中国水利水电科学研究院 | 一种城镇小流域分类分级智能调度方法 |
AU2020101587A4 (en) * | 2020-07-31 | 2020-09-10 | China Institute Of Water Resources And Hydropower Research | Coupling method of 1D and 2D hydrodynamic models based on spatial topology |
CN111651885A (zh) * | 2020-06-03 | 2020-09-11 | 南昌工程学院 | 一种智慧型海绵城市洪涝预报方法 |
-
2020
- 2020-10-13 CN CN202011092833.7A patent/CN112101818B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102296562A (zh) * | 2010-06-25 | 2011-12-28 | 华东电网有限公司 | 防洪与发电相耦合的梯级水库联合洪水调度优化方法 |
CN109284531A (zh) * | 2018-08-03 | 2019-01-29 | 中国水利水电科学研究院 | 一种基于空间拓扑的一二维水动力学耦合方法 |
CN109872063A (zh) * | 2019-02-11 | 2019-06-11 | 南昌工程学院 | 平原城市河网防洪、排涝、活水联合优化调度方法及系统 |
CN109919397A (zh) * | 2019-04-12 | 2019-06-21 | 黄河勘测规划设计研究院有限公司 | 平原河网地区排涝容量校验优化方法 |
CN110838079A (zh) * | 2019-10-18 | 2020-02-25 | 南京中科智慧生态科技有限公司 | 一种智慧城市防洪排水远程监控系统及方法 |
CN111046574A (zh) * | 2019-12-20 | 2020-04-21 | 武汉大学 | 一种平原湖区湖河闸泵系统防洪排涝演算方法 |
CN111274690A (zh) * | 2020-01-16 | 2020-06-12 | 无锡北大建筑工程有限公司 | 一种海绵城市的雨水源头调蓄系统的构建方法 |
CN111476676A (zh) * | 2020-04-02 | 2020-07-31 | 中国水利水电科学研究院 | 一种城镇小流域分类分级智能调度方法 |
CN111651885A (zh) * | 2020-06-03 | 2020-09-11 | 南昌工程学院 | 一种智慧型海绵城市洪涝预报方法 |
AU2020101587A4 (en) * | 2020-07-31 | 2020-09-10 | China Institute Of Water Resources And Hydropower Research | Coupling method of 1D and 2D hydrodynamic models based on spatial topology |
Non-Patent Citations (2)
Title |
---|
SHUHAN ZHANG等: "Storm Water Management and Flood Control in Sponge City Construction of Beijing", 《WATER》, vol. 10, pages 1 - 11 * |
丑天驰: "基于文化框架改进粒子群算法的河湖联控防洪排涝调度模型研究", 《中国优秀硕士学位论文全文数据库(电子期刊)》, no. 3, pages 14 - 30 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112950125A (zh) * | 2021-04-02 | 2021-06-11 | 拉扎斯网络科技(上海)有限公司 | 一种配送资源安全状态数据的获得方法、装置、电子设备 |
CN112989253A (zh) * | 2021-04-02 | 2021-06-18 | 南昌工程学院 | 一种基于JupyterHub的水库优化调度实践教学平台的搭建方法 |
CN115017665A (zh) * | 2022-06-29 | 2022-09-06 | 珠江水利委员会珠江水利科学研究院 | 用于城市暴雨洪涝防御的模拟方法、装置、设备及存储介质 |
CN117236673A (zh) * | 2023-11-16 | 2023-12-15 | 水利部交通运输部国家能源局南京水利科学研究院 | 城市水网多尺度防洪排涝联合优化调度方法和系统 |
CN117236673B (zh) * | 2023-11-16 | 2024-01-26 | 水利部交通运输部国家能源局南京水利科学研究院 | 城市水网多尺度防洪排涝联合优化调度方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112101818B (zh) | 2024-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112101818A (zh) | 一种适用于复杂水力联系的海绵城市洪涝优化调度方法 | |
CN109685685B (zh) | 一种基于宏观配置方案的水资源多目标均衡调度方法 | |
CN108108838B (zh) | 一种高水量利用率的季调节水库优化调度方法 | |
CN109272245B (zh) | 梯级水电站泄洪闸门开度实时自动分配方法及系统 | |
CN108764573B (zh) | 一种面向内陆干旱区的水资源多维均衡配置系统 | |
CN105678067B (zh) | 城市河涌防洪排涝动态预报控制方法及系统 | |
CN113407897B (zh) | 基于多源互济供水模式的分布式水循环模型的设计方法 | |
CN111046574B (zh) | 一种平原湖区湖河闸泵系统防洪排涝演算方法 | |
CN109740808B (zh) | 一种风光水互补发电计划计算方法及系统 | |
CN107239657A (zh) | 一种面向对象的水动力学建模要素管理方法 | |
CN107992989A (zh) | 一种水库群洪水资源化预蓄预泄风险决策方法 | |
CN115271304A (zh) | 基于气候变化的跨流域调水工程优化调度方法 | |
Thomas Jr et al. | On the efficient use of High Aswan Dam for hydropower and irrigation | |
CN109919397B (zh) | 平原河网地区排涝容量校验优化方法 | |
CN117674293A (zh) | 一种梯级水电站的长期发电优化调度方法及装置 | |
Martin | Optimal daily operation of surface-water systems | |
CN108090671A (zh) | 一种多阶段域间-域内反馈迭代的跨流域负荷分配方法 | |
CN109537520A (zh) | 一种建设项目防排洪设小水库留空库容削峰调洪系统及方法 | |
Wu et al. | Multi-objective model of irrigation water distribution based on particle swarm optimization | |
Zhou | Integrated allocation model of water quantity and quality and its application in the Yellow River | |
Puig et al. | Model predictive control for combined water supply and navigability/sustainability in river systems | |
Ghumman et al. | Study of feasibility of night-closure of irrigation canals for water saving | |
CN117650581B (zh) | 一种梯级多电站的联合优化调度方法及系统 | |
CN116562133B (zh) | 一种自流灌溉灌区的渠系多目标配水实时优化方法 | |
Li et al. | Research on the key technology of scheduling and operation of the irrigation and drainage canal |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |