CN107059761A - 水库群防洪库容时空分配设计方法 - Google Patents
水库群防洪库容时空分配设计方法 Download PDFInfo
- Publication number
- CN107059761A CN107059761A CN201710465681.2A CN201710465681A CN107059761A CN 107059761 A CN107059761 A CN 107059761A CN 201710465681 A CN201710465681 A CN 201710465681A CN 107059761 A CN107059761 A CN 107059761A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- subsystem
- big system
- optimization
- flood
- 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
Classifications
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02B—HYDRAULIC ENGINEERING
- E02B1/00—Equipment or apparatus for, or methods of, general hydraulic engineering, e.g. protection of constructions against ice-strains
-
- 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
Abstract
本发明提供一种水库群防洪库容时空分配设计方法,该方法能够在空间和时间上优化设计水库群库容分配方案,同时能兼顾各个子系统所承担调度任务的差异性,其特征在于,包括以下步骤:步骤1,建立水库群系统库容时空分配的优化模型,优化模型有两级目标函数;步骤2,建立水库群大系统库容空间分配优化模型,计算水库群在空间分区上的初步库容分配方案值,实现水库群系统的分区控制;步骤3,在水库群大系统空间分区的基础上,建立各子系统内各水库分期库容分配的多目标优化模型,计算水库群在时间分期上的库容分配;步骤4,步骤2和步骤3的两级优化模型嵌套进行,计算得到水库群大系统在空间分区上和在时间分期上的库容分配非劣解。
Description
技术领域
本发明涉及水库调度技术领域,特别涉及一种水库群防洪库容时空分配设计方法。
背景技术
洪水灾害是我国最严重的自然灾害之一,水库可通过拦洪蓄水来调节水流过程,削减进入下游河道的洪峰流量,从而达到减免洪灾的目的。然而,在洪水调度实践中,过多的考虑防洪安全,较少的考虑水库兴利效益的现象,造成了一系列不合理的问题,洪水资源化这一概念被适时地提出。水库汛限水位是协调水库防洪与兴利矛盾的关键要素。我国洪水资源化的主要途径是开展水库汛限水位设计运用研究,通过调整防洪调度参数和规则来进行,其中主要是分期汛限水位和汛限水位动态控制。汛限水位设计本质上是水库防洪库容的设计,而单库汛限水位设计与运用的基础理论及试点水库的研究工作已经取得了显著的进展,形成了一系列较为成熟的理论方法。相对于单一水库汛限水位问题,水库群汛限水位的联合运用与动态控制问题更加复杂。对于水库群系统,由于水库之间存在一定的水力联系,上下游水库存在库容补偿,单纯提高某一水库的汛限水位,未必能提高梯级水库的洪水资源利用率。而且,随着水库系统中水库数量(维数)的增加,需考虑的信息越来越多,汛限水位的控制也将变得越来越复杂,目前对水库群库容分配问题的研究不多,一般为水库群系统整体构建多目标优化模型,研究侧重点多为优化求解算法,并未深入考虑水库群库容的时空分配问题。
在现有的技术中存在如下问题:(1)目前针对水库群库容分配问题仅限于空间分配角度;(2)目前的水库群联合调度运行库容分配设计方法一般为水库群大系统整体同时优化,需优化的决策变量太多,且未能考虑水库群大系统内各个水库承担调度任务的差异性。
发明内容
本发明是为了解决上述问题而进行的,目的在于提供一种水库群防洪库容时空分配设计方法,该方法能够在空间和时间上优化设计水库群库容分配方案,同时能兼顾水库群大系统内各个子系统所承担调度任务的差异性。
本发明提供了一种水库群防洪库容时空分配设计方法,其特征在于,包括以下步骤:
步骤1,将水库群划分成多个子水库群,其中,水库群视为一个大系统,多个子水库群视为大系统中的多个子系统;建立水库群大系统库容时空分配的优化模型,优化模型有两级目标函数,一级目标函数为水库群大系统整体的目标优化,二级目标函数为水库群大系统中各子系统的多目标优化,优化的决策变量为各子系统库容分配;所述一级目标函数选取为水库群大系统所承担的总调度任务;所述二级目标函数选取为各子系统所承担的调度任务;
步骤2,对水库群大系统进行空间分区,计算水库在空间分区上的各子系统库容分配初始方案值,实现水库群大系统的分区控制,具体实现如下:
步骤2-1,根据包括水库群串并联复杂分布、水库所承担的下游防洪控制点任务在内的多个因素,进行水库群大系统的空间分区,将水库群大系统分解为多个子系统;
步骤2-2,建立水库群大系统库容空间分配优化模型,以步骤1中建立的一级目标函数为优化目标,以各子系统库容分配为决策变量,进行模型优化求解,推求出相应于上述一级目标函数最优的空间分区的各子系统库容分配初始方案值;
步骤3,在水库群大系统空间分区的基础上,对各子系统各水库进行汛期分期,开展各子系统内各水库的分期库容分配设计,计算水库群大系统在时间分期上即各分期之间的库容分配,具体实现如下:
步骤3-1,对各子系统内各水库进行汛期分期,划分成前汛期、主汛期和后汛期三个汛期时段;
步骤3-2,建立各子系统内各水库分期库容分配的多目标优化模型,各子系统以步骤1中建立的二级目标函数为优化目标,以步骤2优化的各子系统库容分配初始方案值为库容约束下限值,以各子系统内各水库分期库容分配为决策变量,开展模型优化求解;
步骤4,步骤2和步骤3的两级优化模型嵌套进行,计算得到水库群大系统在空间分区上和在时间分期上的库容分配非劣解。
上述步骤1中,水库群大系统所承担的总调度任务为总发电量最大、发电效益最大或综合效益最大。
上述步骤1中,水库群大系统中各子系统所承担的分调度任务为防洪效益最大、发电效益最大、供水效益最大、灌溉效益最大、航运效益最大中的任意两个或多个。
上述步骤3中,考虑到各子系统所承担的调度任务存在差异,各子系统在满足步骤2优化的各子系统库容分配初始方案值的基础上,构建适用于各子系统本身的多目标优化模型,选取相应于各子系统本身所承担的调度任务为二级目标函数,单独开展各子系统时间分期上的库容分配的优化求解。
与现有技术相比,本发明的有益效果在于:
(1)本发明提出了一种水库群防洪库容时空分配设计方法,该方法能够同时在空间和时间上优化设计水库群库容分配方案。
(2)本发明所提出的水库群防洪库容时空分配设计方法兼顾考虑了各个子系统所承担调度任务的差异性,使得库容优化分配方案在实际调度过程中更易于操作。
附图说明
图1为本发明实施例一中水库群防洪库容时空分配设计方法的流程图;
图2为本发明实施例一中水库群大系统空间分区示意图,(a)为水库群大系统考虑串并联复杂分布因素的空间分区示意图;(b)为子系统中相邻水库的关系示意图。
具体实施方式
以下结合附图对本发明涉及的一种水库群防洪库容时空分配设计方法的具体实施方案进行详细地说明。
<实施例一>
如图1所示,本实施例一所提供的一种水库群防洪库容时空分配设计方法包括以下步骤:
步骤1,将水库群划分成多个子水库群,其中,水库群视为一个大系统,多个子水库群视为大系统中的多个子系统;建立水库群大系统库容时空分配的优化模型,优化模型有两级目标函数,一级目标函数为水库群大系统整体的目标优化,二级目标函数为水库群大系统中各子系统的多目标优化,优化的决策变量为各子系统库容分配;
一级目标函数可以选取为水库群大系统总发电量最大、发电效益最大或综合效益最大等,以防洪库容为约束条件;例如,若以水库群系统总发电量最大为目标函数,以防洪库容为约束条件,则表达式为:
目标函数:
约束条件:
式中:E总为整个水库群大系统汛期发电量,N为水库群大系统中子系统个数,i为水库大系统中第i个子系统,i=1,2,3…,N,Ei(·)为子系统i的汛期发电量,为子系统i的防洪库容分配值,T为汛期时段长,t为汛期时段长T内的第t时段,t=1,2,3…,T,ηit为子系统i在t时段的发电效率系数,qit为子系统i在t时段的出库流量,hit为子系统i在t时段的发电水头,Vf原为水库大系统原防洪库容值,为子系统i原防洪库容值;
在满足一级目标优化的初步库容分配方案值的基础上,二级目标函数可以选取为子系统所承担的调度任务,所述调度任务可以选择防洪效益最大、发电效益最大、供水效益最大及其他综合利用调度目标中的任意两个或多个,其他综合利用调度目标包括但不限于灌溉效益最大、航运效益最大。
若子系统为防洪为主,兼顾发电与供水,则二级目标可以选取为防洪库容最大、发电量最大、发电保证率最大和供水效益最大,该子系统的目标函数表达式为:
(1)防洪库容最大:
(2)供水效益最大:
(3)发电量最大:
(4)发电保证率最大:
式(1)中:k为洪水场次,Vf,i为相应于第i场次洪水条件下的防洪库容值,i=1,2,…,k,为k个洪水场次情景下的防洪库容平均值;
式(2)中,W为供水效益,T为总供水时段长,Δt为单位供水时段长,为第t时段供水量;
式(3)(4)中,n为计算时段的年数,m为年内的天数,Ni,j为第i年中的第j天的发电量,E为总发电量;Pf为保证出力,#为计数符号,若Ni,j≥Pf则加1,否则加0,Pe为发电保证率。
二级目标函数可通过进化算法、约束法、权重法等方法建立多目标优化模型,并将目标函数变形表达。
步骤2.对水库群大系统进行空间分区,计算水库群在空间分区上的各子系统库容分配初始方案值,实现水库群大系统的分区控制,具体实现如下:
步骤2-1.按水库群串并联复杂分布、水库所承担的下游防洪控制点任务等多个因素,进行水库群大系统的空间分区,将水库群大系统分解为N个子系统;
如图2(a)所示,可依据水库群大系统串并联分布情况,将水库群大系统在空间上分区为子系统A,B,C,…;各子系统分别由串联的水库群组成,例如子系统A中的单库Ai和Ai+1(1≤i≤n)的关系,为子系统A的入库流量,为水库Ai到水库Ai+1的区间入流,如图2(b);子系统A和子系统B构成并联关系,拥有共同的下游防洪控制点a,子系统A和子系统B整体和子系统C构成串联关系,下游有防洪控制点b;
步骤2-2.建立水库群大系统库容空间分配优化模型,以步骤1中建立的一级目标函数为优化目标,以各子系统库容分配为决策变量,进行模型优化求解,推求出相应于上述一级目标函数最优(最大或最小)的空间分区的各子系统库容分配初始方案值 为子系统i的防洪库容分配值,i=1,2,…,N,N为水库大系统中子系统个数。
步骤3.在水库群大系统空间分区的基础上,对各子系统各水库进行汛期分期,开展各子系统内各水库的分期库容分配(分期汛限水位)设计,计算水库群大系统在时间分期上即各分期之间的库容分配,具体实现如下:
步骤3-1.采用变点分析法等常用的分期方法对子系统内各水库进行汛期分期,划分成前汛期主汛期和后汛期三个汛期时段,i=1,2,…,N,N为水库大系统中子系统个数;
步骤3-2.建立各子系统内各水库分期库容分配的多目标优化模型,各子系统以步骤1中建立的二级目标函数为优化目标,以步骤2优化的各子系统库容分配初始方案值 为库容约束下限值,以各水库分期库容分配为决策变量,开展模型优化求解;
考虑到各子系统承担的调度任务存在差异,各子系统可在满足步骤2中库容空间分区的初始分配方案值的基础上,单独开展各子系统时间分期上的库容分配的优化求解,选取不同的二级目标函数,构建适用于各个子系统本身的多目标优化模型;例如某单库以防洪为主,兼顾发电,则二级目标可以选取为防洪库容最大、发电量最大;例如某单库以灌溉为主,兼顾发电,则二级目标可以选取为灌溉保证率最大、发电量最大。
步骤4.步骤2和步骤3的两级优化模型嵌套进行,计算得到水库群大系统在空间分区上和在时间分期上的库容分配非劣解,在子系统内部可依据决策倾向灵活选取具体的水库库容分配实施方案。
应当理解的是,本说明书未详细阐述的部分均属于现有技术。本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (4)
1.一种水库群防洪库容时空分配设计方法,其特征在于,包括以下步骤:
步骤1,将水库群划分成多个子水库群,其中,水库群视为一个大系统,多个子水库群视为大系统中的多个子系统;建立水库群大系统库容时空分配的优化模型,优化模型有两级目标函数,一级目标函数为水库群大系统整体的目标优化,二级目标函数为水库群大系统中各子系统的多目标优化,优化的决策变量为各子系统库容分配;所述一级目标函数选取为水库群大系统所承担的总调度任务;所述二级目标函数选取为各子系统所承担的调度任务;
步骤2,对水库群大系统进行空间分区,计算水库在空间分区上的各子系统库容分配初始方案值,实现水库群大系统的分区控制,具体实现如下:
步骤2-1,根据包括水库群串并联复杂分布、水库所承担的下游防洪控制点任务在内的多个因素,进行水库群大系统的空间分区,将水库群大系统分解为多个子系统;
步骤2-2,建立水库群大系统库容空间分配优化模型,以步骤1中建立的一级目标函数为优化目标,以各子系统库容分配为决策变量,进行模型优化求解,推求出相应于上述一级目标函数最优的空间分区的各子系统库容分配初始方案值;
步骤3,在水库群大系统空间分区的基础上,对各子系统各水库进行汛期分期,开展各子系统内各水库的分期库容分配设计,计算水库群大系统在时间分期上即各分期之间的库容分配,具体实现如下:
步骤3-1,对各子系统内各水库进行汛期分期,划分成前汛期、主汛期和后汛期三个汛期时段;
步骤3-2,建立各子系统内各水库分期库容分配的多目标优化模型,各子系统以步骤1中建立的二级目标函数为优化目标,以步骤2优化的各子系统库容分配初始方案值为库容约束下限值,以各子系统内各水库分期库容分配为决策变量,开展模型优化求解;
步骤4,步骤2和步骤3的两级优化模型嵌套进行,计算得到水库群大系统在空间分区上和在时间分期上的库容分配非劣解。
2.根据权利要求1所述的方法,其特征在于:
步骤1中,水库群大系统所承担的总调度任务为总发电量最大、发电效益最大或综合效益最大。
3.根据权利要求1所述的方法,其特征在于:
步骤1中,水库群大系统中各子系统所承担的分调度任务为防洪效益最大、发电效益最大、供水效益最大、灌溉效益最大、航运效益最大中的任意两个或多个。
4.根据权利要求1所述的方法,其特征在于:
步骤3中,考虑到各子系统所承担的调度任务存在差异,各子系统在满足步骤2优化的各子系统库容分配初始方案值的基础上,构建适用于各子系统本身的多目标优化模型,选取相应于各子系统本身所承担的调度任务为二级目标函数,单独开展各子系统时间分期上的库容分配的优化求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710465681.2A CN107059761B (zh) | 2017-06-19 | 2017-06-19 | 水库群防洪库容时空分配设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710465681.2A CN107059761B (zh) | 2017-06-19 | 2017-06-19 | 水库群防洪库容时空分配设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107059761A true CN107059761A (zh) | 2017-08-18 |
CN107059761B CN107059761B (zh) | 2019-06-25 |
Family
ID=59594212
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710465681.2A Active CN107059761B (zh) | 2017-06-19 | 2017-06-19 | 水库群防洪库容时空分配设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107059761B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108520323A (zh) * | 2018-04-13 | 2018-09-11 | 河海大学 | 一种基于大系统理论时间及空间分解协调的水库调度方法 |
CN109325879A (zh) * | 2018-09-14 | 2019-02-12 | 武汉大学 | 推求水电站中长期调度的变化综合出力系数的时空聚合方法 |
CN110852525A (zh) * | 2019-11-19 | 2020-02-28 | 长安大学 | 基于预报误差越级设防的水库防洪限制水位动态调节方法 |
CN116579626A (zh) * | 2023-05-15 | 2023-08-11 | 长江水利委员会水文局 | 一种基于博弈理论的梯级水库群蓄水策略计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102182159A (zh) * | 2011-03-21 | 2011-09-14 | 武汉大学 | 一种梯级水库汛限水位联合运用调度方法 |
CN103088784A (zh) * | 2013-01-21 | 2013-05-08 | 武汉大学 | 一种梯级水库汛限水位实时动态控制方法 |
CN105676890A (zh) * | 2016-01-22 | 2016-06-15 | 长江水利委员会长江科学院 | 三维及以上梯级水库汛期运行水位动态控制方法 |
CN105931131A (zh) * | 2016-04-18 | 2016-09-07 | 华电电力科学研究院 | 梯级水电站群联合调度效益分配的Shapley值法 |
CN106845712A (zh) * | 2017-01-20 | 2017-06-13 | 中国水利水电科学研究院 | 面向大系统水库群的聚合分解调度规则的提取方法 |
-
2017
- 2017-06-19 CN CN201710465681.2A patent/CN107059761B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102182159A (zh) * | 2011-03-21 | 2011-09-14 | 武汉大学 | 一种梯级水库汛限水位联合运用调度方法 |
CN103088784A (zh) * | 2013-01-21 | 2013-05-08 | 武汉大学 | 一种梯级水库汛限水位实时动态控制方法 |
CN105676890A (zh) * | 2016-01-22 | 2016-06-15 | 长江水利委员会长江科学院 | 三维及以上梯级水库汛期运行水位动态控制方法 |
CN105931131A (zh) * | 2016-04-18 | 2016-09-07 | 华电电力科学研究院 | 梯级水电站群联合调度效益分配的Shapley值法 |
CN106845712A (zh) * | 2017-01-20 | 2017-06-13 | 中国水利水电科学研究院 | 面向大系统水库群的聚合分解调度规则的提取方法 |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108520323A (zh) * | 2018-04-13 | 2018-09-11 | 河海大学 | 一种基于大系统理论时间及空间分解协调的水库调度方法 |
CN108520323B (zh) * | 2018-04-13 | 2022-04-05 | 河海大学 | 一种基于大系统理论时间及空间分解协调的水库调度方法 |
CN109325879A (zh) * | 2018-09-14 | 2019-02-12 | 武汉大学 | 推求水电站中长期调度的变化综合出力系数的时空聚合方法 |
CN109325879B (zh) * | 2018-09-14 | 2022-04-12 | 武汉大学 | 推求水电站中长期调度的变化综合出力系数的时空聚合方法 |
CN110852525A (zh) * | 2019-11-19 | 2020-02-28 | 长安大学 | 基于预报误差越级设防的水库防洪限制水位动态调节方法 |
CN116579626A (zh) * | 2023-05-15 | 2023-08-11 | 长江水利委员会水文局 | 一种基于博弈理论的梯级水库群蓄水策略计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107059761B (zh) | 2019-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103745023B (zh) | 水电站发电出力方案制作与最优负荷分配耦合建模方法 | |
CN104167730B (zh) | 一种复杂约束限制下的梯级水电站群实时优化调度方法 | |
CN109449971B (zh) | 一种新能源消纳的多目标电力系统源荷互动优化调度方法 | |
CN107059761B (zh) | 水库群防洪库容时空分配设计方法 | |
CN102298731A (zh) | 考虑顶潮供水综合要求的梯级水库短期发电优化调度方法 | |
CN104063808B (zh) | 一种跨省送电梯级水电站群调峰调度两阶段搜索方法 | |
CN104616069A (zh) | 一种考虑计划完成率和负荷率均衡的年度发电计划滚动分解优化方法 | |
CN104701888B (zh) | 一种电网水、火电工况自匹配模式调整实时负荷的方法 | |
CN104376412A (zh) | 一种新能源发电模式下高载能企业调峰控制方法 | |
CN105427017B (zh) | 一种水电富集电网特大规模电站群短期计划编制方法 | |
CN107038151B (zh) | 基于蓄能控制的梯级优化调度图绘制方法及调度方法 | |
CN102518092A (zh) | 梯级水电站水库防凌库容优化配置设计方法 | |
CN110867903B (zh) | 一种基于深度神经网络的智能电网系统实时最优功率分配方法 | |
CN103023073A (zh) | 一种水电站群电站与机组混合优化调度方法 | |
CN103199538B (zh) | 日前计划多尺度逼近的发电计划闭环控制方法 | |
CN103617455A (zh) | 基于虚拟机组子群的网厂两级负荷优化调度方法 | |
CN116231765B (zh) | 一种虚拟电厂出力控制方法 | |
CN110880761A (zh) | 一种基于p2p平台的风光抽蓄混合系统的联合调度方法 | |
CN111030101B (zh) | 一种基于多元化大数据清洁能源消纳联动调控方法及系统 | |
CN112581310A (zh) | 梯级水电站群发电指标分配方法 | |
CN109636015B (zh) | 一种梯级水电虚拟抽蓄电站调度方法 | |
CN105868841B (zh) | 一种基于风电优先上网的风水火联合调度方法 | |
CN108062022A (zh) | 一种热电协同系统的优化控制方法 | |
CN107742166B (zh) | 一种充分灌溉条件下直接补库的多库-多站系统水资源优化配置方法 | |
CN114186877B (zh) | 一种考虑水库调节库容的日内水光互补计算方法 |
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 |