CN112766565B - 梯级水库群防洪优化调度分数阶逐步优化方法 - Google Patents
梯级水库群防洪优化调度分数阶逐步优化方法 Download PDFInfo
- Publication number
- CN112766565B CN112766565B CN202110054616.7A CN202110054616A CN112766565B CN 112766565 B CN112766565 B CN 112766565B CN 202110054616 A CN202110054616 A CN 202110054616A CN 112766565 B CN112766565 B CN 112766565B
- Authority
- CN
- China
- Prior art keywords
- reservoir
- scheduling
- flood control
- fractional order
- reservoir group
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 73
- 238000005457 optimization Methods 0.000 title claims abstract description 36
- 238000004364 calculation method Methods 0.000 claims abstract description 29
- 230000008569 process Effects 0.000 claims abstract description 26
- 239000000126 substance Substances 0.000 claims description 3
- 230000007787 long-term memory Effects 0.000 abstract description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 19
- 230000008901 benefit Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 230000004069 differentiation Effects 0.000 description 4
- 230000015654 memory Effects 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000007123 defense Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 241000239290 Araneae Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000011144 upstream manufacturing 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Economics (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Data Mining & Analysis (AREA)
- General Business, Economics & Management (AREA)
- Development Economics (AREA)
- Primary Health Care (AREA)
- Entrepreneurship & Innovation (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Algebra (AREA)
- Quality & Reliability (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及梯级水库群防洪调度技术领域,公开了一种梯级水库群防洪优化调度分数阶逐步优化方法,以逐步优化方法为基础框架,在给定初始调度过程后,首先将多阶段水库群防洪调度问题分解成多个两阶段子问题;然后对每个子问题采用分数阶微分进行转换,利用具有长时记忆能力的分数阶更新方程进行求解;最后通过迭代寻优逐次逼近全局最优解。本发明梯级水库群防洪优化调度分数阶逐步优化方法,操作简单、方便易行,避免传统方法的维数灾问题,计算速度快、搜索能力强、计算精度高,结果有效合理,是大规模水库群优化调度问题求解的一种新方法。
Description
技术领域
本发明涉及梯级水库群防洪调度技术领域,具体涉及一种梯级水库群防洪优化调度分数阶逐步优化方法。
背景技术
洪涝灾害历来是中华民族心腹之患,加快经济发展方式转变,必须营造安全、和谐、秀美的环境,这对洪水防御提出了更高要求。水库是流域治理开发中的重要工程措施,具有不可替代的防洪作用。伴随着大江大河建设进程的逐步推进,防洪调度是国内外研究的热点。
水库群防洪优化调度旨在以确保大坝枢纽安全为基础,通过有效拦蓄、削峰、调峰,最大程度地避免或减少洪灾损失。但水库群防洪优化调度具有强约束、高维度、多阶段和非线性,其高效求解备受关注。传统动态规划能较好地解决单库防洪优化调度,但在求解梯级水库防洪优化调度问题存在维数灾问题。同时,粒子群算法、差分进化算法、群居蜘蛛算法等智能优化算法在求解该问题时,也易早熟收敛,存在全局搜索能力不足的缺陷。
逐步优化方法(POA)是一种求解多阶段多状态优化问题的改进动态规划方法,其基本原理是将多阶段问题分解成多个两阶段子问题,每次计算均固定其他阶段变量,只针对当前阶段及相邻两阶段的目标函数值进行迭代寻优,如此逐阶段进行计算,直至输出最终解。该方法减少了可行解的离散状态,相比动态规划方法可大幅提高计算效率,但当水库群数目较多时,仍会面临严重的维数灾问题。因此,亟待对POA进行有效改进,提高水库群联合调度问题的计算效率。
分数阶微分作为数学分析的一个分支,由整数阶微积分发展而来,其阶次不再仅限于整数而推广到任意实数阶次,已成功应用于流体动力学、生物工程、系统工程控制、拟合回归模型等领域。由于分数阶导数项是无限的,且隐含的表示过去所有事件信息,故分数阶具有记忆性,使得分数阶微分非常适合用来描述具有记忆性质的状态变化情况,抗干扰能力强。
考虑到分数阶微分固有的长期记忆特性十分适合描述搜索寻优过程,本次发明将POA的历史轨迹状态进行记忆和学习,建立考虑长时记忆特性的分数阶微分方程,通过调整控制参数进行迭代寻优,来有效克服POA的维数灾问题,提高方法的求解效率和计算精度,对梯级水库群优化调度问题具有良好的支撑应用价值。
发明内容
本发明的目的就是针对上述技术的不足,提供一种梯级水库群防洪优化调度分数阶逐步优化方法,解决POA在水库群防洪优化调度问题求解时存在的维数灾缺陷。
为实现上述目的,本发明所设计的梯级水库群防洪优化调度分数阶逐步优化方法,包括如下步骤:
1)确定初始计算条件,包括水库群防洪优化调度的目标函数、约束条件和决策变量,设定计算参数,包括最大迭代次数M、各水库初始离散步长、初始离散数目K、收敛精度ε、水库数量N、时段总数T;
3)设定迭代次数m=1;
5)从m=M′+1次开始,执行分数阶逐步优化方法,该方法包括如下步骤:
①首先,将调度问题分解成T-1个两阶段问题;
②然后,对各两阶段问题构建分数阶微分方程,执行分数阶寻优,在第m次迭代的第s个两阶段问题中,为水库i在两阶段s的调度过程,Δi,s为水库i在两阶段s的离散步长,为水库i在两阶段s时第m次迭代选择的离散数目,有根据Grünwald-Letnikov定义,POA中变化关系的优化过程转换为:
7)令m=m+1,若m>M,则转至步骤8),否则转至步骤5);
8)停止计算,输出最终的调度过程。
优选地,所述初始离散数目K为奇数。
优选地,所述分数阶次0.8≤α<1。
本发明与现有技术相比,具有以下优点:充分利用了历史轨迹知识信息,采用分数阶微分方程替代两阶段子问题中的全面枚举操作,操作简单、方便易行,计算速度快、搜索能力强,有效克服了传统算法的维数灾问题。
附图说明
图1为本发明梯级水库群防洪优化调度分数阶逐步优化方法的流程图;
图2是本发明方法在1965年1%来水条件下两河口水库的防洪调度计算结果图;
图3是本发明方法在1965年1%来水条件下锦屏一级水库的防 洪调度计算结果图;
图4是本发明方法在1965年1%来水条件下二滩水库的防洪调度计算结果图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的详细说明。
如图1所示,一种梯级水库群防洪优化调度分数阶逐步优化方法,包括如下步骤:
1)确定初始计算条件,包括水库群防洪优化调度的目标函数、约束条件和决策变量,设定计算参数,包括最大迭代次数M、各水库初始离散步长、初始离散数目K、收敛精度ε、水库数量N、时段总数T,初始离散数目为奇数,其中,水库群防洪优化调度模型可以描述为:已知各水库调度期初始水位、末水位和入库洪水过程以及区间洪水过程,在满足各水库及防洪控制站复杂约束的情况下,以水库群各级水库出库流量过程平方和最小为优化目标,水库群防洪优化调度的目标函数如下:
式中:N为水库数量;i为水库序号,i=1,2,…,N,T是调度期的时段总数,j为时段序号,且j=1,2,…,T,Qi,j为水库i-1在j时段的出库流量(m3/s),Ri,j为水库i+1在j时段的区间流量(m3/s);
需要满足的约束条件,主要包括:
(1)水量平衡约束:
式中,Vi,j为水库i在j时段末的库容值(m3),Ii,j为水库i在j时段的入库流量(m3/s),Δj为时段长度(h);
(2)水力联系约束:
式中,Qi-1,j为上游水库i-1在j时段的出库流量(m3/s),Ri,j为水库i在j时段的区间流量(m3/s);
(3)时段水位约束:
(4)出库流量约束:
(5)初末水位约束:
Zi,O=Zi,start,Zi,T=Zi,end
式中,Zi,start和Zi,end分别为水库i的调度期初始水库水位(m)和调度期末控制水位(m);
(6)水库泄流能力约束
式中,Qi,max(Zi,j)为水库i在j时段相应于水位Zi,j的最大下泄能力;
(7)非负约束:各种变量均为非负值;
决策变量为水位;
3)设定迭代次数m=1;
5)从m=M′+1次开始,执行分数阶逐步优化方法,该方法包括如下步骤:
①首先,将调度问题分解成T-1个两阶段问题;
②然后,对各两阶段问题构建分数阶微分方程,执行分数阶寻优,在第m次迭代的第s个两阶段问题中,为水库i在两阶段s的调度过程,Δi,s为水库i在两阶段s的离散步长,为水库i在两阶段s时第m次迭代选择的离散数目,有本 次对这一变化关系引入分数阶微分思想,认为任意调度过程变化与之前的历史状态有关,具有长时记忆特性,根据Grünwald-Letnikov定义,分数阶微分表达式为:
从其定义来看,分数阶微分与之前所有的历史状态有关,具有长期连续记忆特性,其离散形势的近似实现为:
式中,Dα[b(t)]为序列b(t)的分数阶微分,Γ()为Gamma函数,α为分数阶次,U为采样周期,本次迭代次数更新一次,调度过程改变一次,故U=1,r为截断阶次,根据需要选择记忆中前r次的历史状态,在本专利中取r=4,得到:
式中,α为分数阶次,通过改变参数值来调节状态过程对当前状态的影响程度,达到控制求解方法的搜索速度和计算精度的目的,为区间上的随机整数,0.8≤α<1,通过在相应取值范围内随机生成α和进行调整,不断优化本次状态,来优化更新调度过程;
7)令m=m+1,若m>M,则转至步骤8),否则转至步骤5);
8)停止计算,输出最终的调度过程。
以我国雅砻江为例进行研究。
雅砻江是金沙江第一大支流,干流全程约1571km,水能资源丰富且非常集中。按照长江防御洪水方案要求,雅砻江梯级水库中两河口、锦屏一级和二滩水库分别为长江干流预留防洪库容20亿m3、16亿m3和9亿m3。现以两河口、锦屏一级和二滩组成的梯级水库群防洪优化调度问题为例进行研究。
选取1965年四种来水频率(1%、2%、3.3%、5%)作为来水系列,分别采用本发明方法与POA开展梯级水库群防洪优化调度,结合图2、图3及图4,表1列出了本发明方法与DDDP的防洪优化调度计算结果。
表1 本发明方法和POA计算结果对比情况
结果表明,本发明方法相比于POA,其优势体现在:①在所有来水情况下,本发明方法的计算结果与POA基本一致,均小于POA计算结果,且出库过程较为平稳,最末梯级二滩水库的出库流量也更小;②本发明方法计算时间均小于POA,耗时仅为POA的12%,随着水库计算规模增加,计算性能优势将更为凸显。
可见,本发明方法计算结果更优、计算效率越快,表明了本发明方法应用于水库群防洪优化调度问题高效求解的可行性和有效性
显然可见,本发明与POA相比,通过分数阶微分来向学习来优化搜索过程,搜索能力强,计算效率快,为求解高维度、非线性、强约束的梯级水库群防洪优化调度提供了一种实用有效的工具,对于流域防洪具有明显的错峰削峰、补偿调节和防灾减灾作用。
Claims (3)
1.一种梯级水库群防洪优化调度分数阶逐步优化方法,其特征在于:包括如下步骤:
1)确定初始计算条件,包括水库群防洪优化调度的目标函数、约束条件和决策变量,设定计算参数,包括最大迭代次数M、各水库初始离散步长、初始离散数目K、收敛精度ε、水库数量N、时段总数T;
3)设定迭代次数m=1;
5)从m=M′+1次开始,执行分数阶逐步优化方法,该方法包括如下步骤:
①首先,将调度问题分解成T-1个两阶段问题;
②然后,对各两阶段问题构建分数阶微分方程,执行分数阶寻优,在第m次迭代的第s个两阶段问题中, 为水库i在两阶段s的调度过程,Δi,s为水库i在两阶段s的离散步长,为水库i在两阶段s时第m次迭代选择的离散数目,有根据Grünwald-Letnikov定义,POA中变化关系的优化过程转换为:
7)令m=m+1,若m>M,则转至步骤8),否则转至步骤5);
8)停止计算,输出最终的调度过程。
2.根据权利要求1所述梯级水库群防洪优化调度分数阶逐步优化方法,其特征在于:所述初始离散数目K为奇数。
3.根据权利要求1所述梯级水库群防洪优化调度分数阶逐步优化方法,其特征在于:所述分数阶次0.8≤α<1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110054616.7A CN112766565B (zh) | 2021-01-15 | 2021-01-15 | 梯级水库群防洪优化调度分数阶逐步优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110054616.7A CN112766565B (zh) | 2021-01-15 | 2021-01-15 | 梯级水库群防洪优化调度分数阶逐步优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112766565A CN112766565A (zh) | 2021-05-07 |
CN112766565B true CN112766565B (zh) | 2022-05-03 |
Family
ID=75701810
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110054616.7A Active CN112766565B (zh) | 2021-01-15 | 2021-01-15 | 梯级水库群防洪优化调度分数阶逐步优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112766565B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113971362B (zh) * | 2021-10-26 | 2022-05-10 | 中国水利水电科学研究院 | 基于能量准则的混联水库群防洪优化调度方案生成方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106937295A (zh) * | 2017-02-22 | 2017-07-07 | 沈阳航空航天大学 | 基于博弈论的异构网络高能效功率分配方法 |
CN110286595A (zh) * | 2019-08-12 | 2019-09-27 | 金陵科技学院 | 一类受饱和非线性输入影响的分数阶系统自适应控制方法 |
CN111092441A (zh) * | 2019-11-29 | 2020-05-01 | 江苏方天电力技术有限公司 | 抽水蓄能电站区域负荷频率控制方法和系统 |
-
2021
- 2021-01-15 CN CN202110054616.7A patent/CN112766565B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106937295A (zh) * | 2017-02-22 | 2017-07-07 | 沈阳航空航天大学 | 基于博弈论的异构网络高能效功率分配方法 |
CN110286595A (zh) * | 2019-08-12 | 2019-09-27 | 金陵科技学院 | 一类受饱和非线性输入影响的分数阶系统自适应控制方法 |
CN111092441A (zh) * | 2019-11-29 | 2020-05-01 | 江苏方天电力技术有限公司 | 抽水蓄能电站区域负荷频率控制方法和系统 |
Non-Patent Citations (4)
Title |
---|
Multi-objective Decision Making Method for Reservoir Flood Operation Based on Cumulative Prospect Theory and Maximum Entropy Theory;Zou Qiang etal.;《Water Resources and Power》;20181231;第57-60页 * |
Research for optimal flood dispatch model for reservoir based on POA;Qin Xubao etal.;《Water Resources and Power》;20081231;第60-62页 * |
基于分数信号处理的变步长自适应振动主动控制方法研究;杨晓京 等;《重庆邮电大学学报( 自然科学版)》;20191031;第31卷(第5期);第729-736页 * |
长江上游水库群多目标优化调度模型及应用研究I:模型原理及求解;黄草 等;《水利学报》;20140930;第1009-1018页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112766565A (zh) | 2021-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106951985B (zh) | 一种基于改进人工蜂群算法的梯级水库多目标优化调度方法 | |
CN107491635B (zh) | 一种基于嵌套降维算法的梯级水库水沙联合优化调度方法 | |
CN103971174B (zh) | 基于改进量子粒子群算法的水电站群优化调度方法 | |
CN109447405B (zh) | 一种承担调峰任务的一库多级式梯级库群短期计划制定方法 | |
CN109345010B (zh) | 一种梯级泵站的多目标优化调度方法 | |
CN109858173B (zh) | 一种风光水互补发电系统联合调度的多目标正切方法 | |
CN110598919B (zh) | 一种梯级水电站群生态调控方法与系统 | |
CN109636043A (zh) | 一种梯级水电系统发电调度自适应优化方法与系统 | |
CN102043905A (zh) | 基于自适应算法的小水电群智能优化错峰调度方法 | |
CN110598983A (zh) | 一种自适应改进粒子群算法的梯级水库优化调度方法 | |
CN104166887A (zh) | 一种梯级水电站群联合优化调度的正交离散微分动态规划方法 | |
CN112766565B (zh) | 梯级水库群防洪优化调度分数阶逐步优化方法 | |
CN107704965A (zh) | 一种兼顾生态保护需求的水库兴利调度图优化方法 | |
CN111476475A (zh) | 多约束条件下梯级水电站短期优化调度方法 | |
CN110766210B (zh) | 一种梯级水库群短期优化调度方法与系统 | |
CN108537370A (zh) | 基于混合智能降维算法的特大流域水电站群优化调度方法 | |
CN111476477A (zh) | 基于发电效益目标的梯级水电站中长期优化调度方法 | |
CN105160443A (zh) | 基于扩展线性二次高斯方法的复杂水库群优化调度方法 | |
CN111861137A (zh) | 梯级水库群并行多目标调度方法 | |
CN105976058B (zh) | 一种基于灰色系统预测的离散微分动态规划方法 | |
CN109635999A (zh) | 一种基于粒子群-细菌觅食的水电站调度方法及系统 | |
CN112766564B (zh) | 梯级水电站群优化调度逐阶段反向学习降维优化方法 | |
CN106779182B (zh) | 一种基于动态规划的水库调度图优化方法 | |
CN112633578B (zh) | 一种引调水工程影响下梯级水库群优化调度方法 | |
CN116960988B (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 |