CN110880048B - 一种梯级水库生态随机优化调度模型及求解方法 - Google Patents

一种梯级水库生态随机优化调度模型及求解方法 Download PDF

Info

Publication number
CN110880048B
CN110880048B CN201911077800.2A CN201911077800A CN110880048B CN 110880048 B CN110880048 B CN 110880048B CN 201911077800 A CN201911077800 A CN 201911077800A CN 110880048 B CN110880048 B CN 110880048B
Authority
CN
China
Prior art keywords
reservoir
period
nth
ecological
time period
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
Application number
CN201911077800.2A
Other languages
English (en)
Other versions
CN110880048A (zh
Inventor
李黄强
童华敏
许阳
刘江鹏
王炜
徐康
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yichang Power Supply Co of State Grid Hubei Electric Power Co Ltd
Original Assignee
Yichang Power Supply Co of State Grid Hubei Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yichang Power Supply Co of State Grid Hubei Electric Power Co Ltd filed Critical Yichang Power Supply Co of State Grid Hubei Electric Power Co Ltd
Priority to CN201911077800.2A priority Critical patent/CN110880048B/zh
Publication of CN110880048A publication Critical patent/CN110880048A/zh
Application granted granted Critical
Publication of CN110880048B publication Critical patent/CN110880048B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • 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
    • 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply

Abstract

一种梯级水库生态随机优化调度模型及求解方法，优化调度模型的建立包括：首先建立修正全年流量偏差AAPFD指标最小的梯级水库生态随机调度模型目标函数；接着描述梯级水库生态随机调度模型的约束条件，包括：水位约束、下泄流量约束、库容约束以及保证出力约束。利用强化学习中的SARSA算法对该生态随机调度模型进行求解将克服随机动态规划算法所带来的维数灾问题。本发明能在保证梯级水库水电站出力的前提下提高生态需求，对于水库水资源的生态利用具有较高指导意义。同时提出的SARSA算法能够很好地在较短时间内获得梯级水库随机生态优化调度模型的结果。

Description

一种梯级水库生态随机优化调度模型及求解方法
技术领域
本发明一种梯级水库生态随机优化调度模型及求解方法，属于水库优化调度领域。
背景技术
水电作为目前开发技术最为成熟、开发规模最大的可再生能源，以其低廉运行成本、良好调节性能和快速负荷响应能力，在全球电力能源供给中发挥着重要作用。水库具有发电、防洪、灌溉等兴利作用。但是水库在实际运行过程中，为提高发电等效益，常会改变水库所在河流的天然径流，破坏了河流水文变化的天然规律，给河流生态环境带来极大影响，造成河流生态环境恶化。随着我国加快推进绿色发展，实行最严格的环境保护制度，要求水库在生态保护中发挥作用。水库生态调度在传统调度基础上增加生态因子的考虑，以改善流域生态环境为目标，优化水资源的可持续利用性，实现河流生态系统的良性循环。
发明内容
本发明目的一方面在于提供一种梯级水库生态随机优化调度模型，通过引入修正全年流量偏差AAPFD，以修正全年流量偏差AAPFD平方的最小期望值为目标函数，同时考虑水库保证出力约束及其他约束，确定对应的最优放水策略。本发明目的另一方面在于针对所述梯级水库生态随机优化调度模型而提供一种求解方法，该方法针对水库中长期来水的随机性，通过建立梯级水库入库流量的随机过程，利用强化学习中的SARSA算法对模型进行求解，使得求解结果计算时间更短，且求解结果更符合实际。
本发明采取的技术方案为：
一种梯级水库生态随机优化调度模型，该模型的目标函数为：
Figure BDA0002263031370000011
式中：f*为梯级水库最小AAPFD期望值的平方，N为梯级水电站的个数，T为调度周期总的时段数，Ns为入库径流离散状态数，Rn,t为第n个水库第t时段的指标函数，fn,t为第n个水库第t时段的AAPFD值。
Rn,t、fn,t表达式如下所示：
Figure BDA0002263031370000021
式中：Vn,i为第n级电站i状态时对应的水位，
Figure BDA0002263031370000022
为第n级电站第t时段的出库流量，
Figure BDA0002263031370000023
为第n个水库第t时段离散的j个入库流量值，j、k为入库径流相互之间的状态转移，
Figure BDA0002263031370000024
为第n个水库第t时段的出库流量，
Figure BDA0002263031370000025
为第n个水库第t时段的入库径流，
Figure BDA0002263031370000026
为第n个水库第t时段处于j状态转移到k状态的转移概率，
Figure BDA0002263031370000027
为调度周期天然入库流量的平均值。
所述目标函数约束条件为：
①、龙头水库水量平衡约束：
Figure BDA0002263031370000028
上式中：V1 t+1为龙头水库第t+1时段的水位，V1 t为龙头水库第t时段的水位，
Figure BDA0002263031370000029
为龙头水库第t+1时段的入库流量，
Figure BDA00022630313700000210
为龙头水库第t+1时段的出库流量，△t为面临时段的时长；
②、其他水库水量平衡约束：
Figure BDA00022630313700000211
上式中：
Figure BDA00022630313700000212
为第n+1级水库第t+1时段的水位，
Figure BDA00022630313700000213
为第n+1级水库第t时段的水位，
Figure BDA00022630313700000214
为第n级水库第t时段的出库流量，
Figure BDA00022630313700000215
为第t时段第n+1级水库上游区间来水流量，
Figure BDA00022630313700000216
为第n+1级水库第t时段的出库流量，△t为面临时段的时长，N为梯级水库的总数。
③、水库蓄水量约束：
Figure BDA00022630313700000217
上式中：Vn,min为第n级水库的最小蓄水量，
Figure BDA00022630313700000218
为第n级水库第t时段的水位，Vn,max为第n级水库的最大蓄水量，N为梯级水库的总数。
④、水库下泄流量约束：
Figure BDA00022630313700000219
上式中：Qn,min,ck为第n级水库最小下泄流量，
Figure BDA00022630313700000220
为第n级水库第t时段的下泄流量，Qn,max,ck为第n级水库最大下泄流量，N为梯级水库的总数。
⑤、保证出力约束：Wn,bz≤Wn≤Wn,zj n＝1,...,N；
上式中：Wn,bz为第n级水库的保证出力，Wn为第n级水库的出力，Wn,jz为第n级水库的机组最大出力。
所述优化调度模型中，对梯级电站按照上下游的顺序从上到下依次编号为1,...,n，并根据上下游水库之间的水量关系，建立梯级各水库的水量平衡方程：
Figure BDA0002263031370000031
上式中：Ln+1,t为流域内t时段两级水库之间的区间来水。
所述优化调度模型中，由于水库中长期径流具有随机性，因此各个水库的入库径流不再以传统的方式，采用某一特定年份历史观测数据作为水库的入库径流。而是采用随机入库过程作为梯级水库入库径流的表示方式，即根据皮尔逊Ⅲ型分布法计算所得频率离散为从小到大的各水库各时段随机来水的拟合值。
所述优化调度模型中，生态指标利用修正全年流量偏差AAPFD指标，该指标反映一个调度周期内流量变化对河流生态环境的影响，其值越小，表示水库调度后流量变化对河流生态系统的影响越小，河流生态环境越好，
对于单一水库，修正全年流量偏差AAPFD表示为：
Figure BDA0002263031370000032
式中：T为调度周期总时段数，
Figure BDA0002263031370000033
为第t时段水库出库流量值；
Figure BDA0002263031370000034
为第t时段水库天然入库流量；
Figure BDA0002263031370000035
为调度周期天然入库流量的平均值。
采用SARSA算法，对梯级水库生态调度随机优化模型进行优化计算，得到修正全年流量偏差AAPFD期望值的平方最大，以及各水库各时段的最优出库流量
Figure BDA0002263031370000036
包括以下步骤：
步骤1：根据梯级水库历年入库流量，求解各水库各时段的皮尔逊Ⅲ型概率密度曲线的统计参数，并根据统计样本参数绘制各水库各时段的径流频率曲线，最终通过皮尔逊Ⅲ型分布模拟得出的径流随机值计算状态转移概率矩阵。
步骤2：对每一个水库任一时段始末库容分别为Vn,i,Vn,j时，进行阶段收益求解，通过状态转移矩阵求解和Bellman方程求解在每一个水库在下一时刻水位随机情况下的余留收益，从而得到强化学习算法中的收益矩阵。
步骤3：根据SARSA算法试错的原理，从第一阶段梯级水库起调水位V0开始进行尝试，对每一组[Vi,Vj]均对出库寻优计算，直到进行到最后一个阶段，完成一次寻优，然后根据时序差分法TD通过下一时刻的价值更新前一时刻的价值的思想对值函数方程，进行如下计算：
Figure BDA0002263031370000041
式中：
Figure BDA0002263031370000042
为t+1时段采用策略π得到状态动作对(s,a)时的Q值函数，
Figure BDA0002263031370000043
为t时段采用策略π得到状态动作对(s,a)时的Q值函数，α为SARSA算法的学习率，
Figure BDA0002263031370000044
为采取动作a后状态由s转移到s′的回报值。γ为SARSA算法的奖赏折扣率，
Figure BDA0002263031370000045
为t时段采用策略π得到状态动作对(s′,a′)时的Q值函数。
通过对Q值函数进行不断的试错更新，最终形成一个以s，a分别为横纵坐标，以离散的水位数为维数的Q矩阵。
步骤4：重复对当前状态Vi进行观察，并根据贪婪决策选择下一阶段的执行动作，然后对当前Q矩阵进行更新，直到Q矩阵收敛。然后从Q矩阵中选取最大值作为策略，即可得到修正全年流量偏差(AAPFD)期望值的平方最大，以及各水库各时段的最优出库流量
Figure BDA0002263031370000046
本发明一种梯级水库生态随机优化调度模型及求解方法，技术效果如下：
1)所述模型能在保证梯级水库水电站出力的前提下提高生态需求，对于水库水资源的生态利用具有较高指导意义。同时提出的SARSA算法能够很好地在较短时间内获得梯级水库随机生态优化调度模型的结果。
2)利用强化学习中的SARSA算法，对该生态随机调度模型进行求解，将克服随机动态规划算法所带来的维数灾问题。
3)采用随机模型，考虑不同时段之间入库径流的马尔可夫性，能够更好的反映天然径流的随机，避免了采用特定年份历史数据作为入库径流时的偶然性。
4)考虑生态的调度模型相较于单纯的发电模型，其得出的的调度规则能够更好地适应流域内的生态环境，最终达到发电生态多目标最优化的效果。
附图说明
图1为本发明的模型求解方法流程图。
图2为水位库容曲线图。
图3为水库皮尔逊Ⅲ型概率密度曲线图。
具体实施方式
一种梯级水库生态随机优化调度模型，该模型的目标函数为：
Figure BDA0002263031370000051
式中：f*为梯级水库最小AAPFD期望值的平方，N为梯级水电站的个数，T为调度周期总的时段数，Ns为入库径流离散状态数，Rn,t为第n个水库第t时段的指标函数，fn,t为第n个水库第t时段的AAPFD值；
Rn,t、fn,t表达式如下所示：
Figure BDA0002263031370000052
式中：Vn,i为第n级电站i状态时对应的水位，
Figure BDA0002263031370000053
为第n级电站第t时段的出库流量，
Figure BDA0002263031370000054
为第n个水库第t时段离散的j个入库流量值，j、k为入库径流相互之间的状态转移，
Figure BDA0002263031370000055
为第n个水库第t时段的出库流量，
Figure BDA0002263031370000056
为第n个水库第t时段的入库径流，
Figure BDA0002263031370000057
为第n个水库第t时段处于j状态转移到k状态的转移概率，
Figure BDA0002263031370000058
为调度周期天然入库流量的平均值。
所述目标函数约束条件为：
①、龙头水库水量平衡约束：
Figure BDA0002263031370000059
上式中：V1 t+1为龙头水库第t+1时段的水位，V1 t为龙头水库第t时段的水位，
Figure BDA00022630313700000510
为龙头水库第t+1时段的入库流量，
Figure BDA00022630313700000511
为龙头水库第t+1时段的出库流量，△t为面临时段的时长；
②、其他水库水量平衡约束：
Figure BDA00022630313700000512
上式中：
Figure BDA0002263031370000061
为第n+1级水库第t+1时段的水位，
Figure BDA0002263031370000062
为第n+1级水库第t时段的水位，
Figure BDA0002263031370000063
为第n级水库第t时段的出库流量，
Figure BDA0002263031370000064
为第t时段第n+1级水库上游区间来水流量，
Figure BDA0002263031370000065
为第n+1级水库第t时段的出库流量，△t为面临时段的时长，N为梯级水库的总数。
③、水库蓄水量约束：
Figure BDA0002263031370000066
上式中：Vn,min为第n级水库的最小蓄水量，
Figure BDA0002263031370000067
为第n级水库第t时段的水位，Vn,max为第n级水库的最大蓄水量，N为梯级水库的总数。
④、水库下泄流量约束：
Figure BDA0002263031370000068
上式中：Qn,min,ck为第n级水库最小下泄流量，
Figure BDA0002263031370000069
为第n级水库第t时段的下泄流量，Qn,max,ck为第n级水库最大下泄流量，N为梯级水库的总数。
⑤、保证出力约束：Wn,bz≤Wn≤Wn,zj n＝1,...,N；
上式中：Wn,bz为第n级水库的保证出力，Wn为第n级水库的出力，Wn,jz为第n级水库的机组最大出力。
所述优化调度模型中，对梯级电站按照上下游的顺序从上到下依次编号为1,...,n，并根据上下游水库之间的水量关系，建立梯级各水库的水量平衡方程：
Figure BDA00022630313700000610
上式中：Ln+1,t为流域内t时段两级水库之间的区间来水。
所述优化调度模型中，由于水库中长期径流具有随机性，因此各个水库的入库径流不再以传统的方式，采用某一特定年份历史观测数据作为水库的入库径流。而是采用随机入库过程作为梯级水库入库径流的表示方式，即根据皮尔逊Ⅲ型分布法，计算所得频率离散为从小到大的各水库各时段随机来水的拟合值。
所述优化调度模型中，生态指标利用修正全年流量偏差AAPFD指标，该指标反映一个调度周期内流量变化对河流生态环境的影响，其值越小，表示水库调度后流量变化对河流生态系统的影响越小，河流生态环境越好。
对于单一水库，修正全年流量偏差AAPFD表示为：
Figure BDA0002263031370000071
式中：T为调度周期总时段数，
Figure BDA0002263031370000072
为第t时段水库出库流量值；
Figure BDA0002263031370000073
为第t时段水库天然入库流量；
Figure BDA0002263031370000074
为调度周期天然入库流量的平均值。
采用SARSA算法，对梯级水库生态调度随机优化模型进行优化计算，得到修正全年流量偏差AAPFD期望值的平方最大，以及各水库各时段的最优出库流量
Figure BDA0002263031370000075
具体包括以下步骤：
S1.1：根据梯级水库中库容与水位的长序列数据进行拟合，得到梯级各水库的水位库容曲线，如图2所示，利用梯级水库中库容与水位的长序列数据，导入matlab中绘制出如图2的V-Z二次拟合曲线，并求解出二次曲线所代表的二次轨迹方程。由此可知在任意水位下的库容值。
S1.2：根据梯级水库历年入库流量，求出各水库入库流量的皮尔逊Ⅲ型概率密度曲线的统计参数：入库径流平均值，均方差、变差系数和偏差系数。并根据所求参数绘制出皮尔逊Ⅲ型概率密度曲线，如图3所示，由求得的入库径流平均值，均方差、变差系数和偏差系数等相关统计参数，在海参几率格纸上绘制各个水库的皮尔逊Ⅲ型概率密度曲线，并通过曲线筛选出水库各时段对应来水频率的模拟径流随机值。
S1.3：求出各水库来水在相邻两个时段之间的相关性以及转移概率：
马尔可夫相关性检验公式：
Figure BDA0002263031370000076
式中，rt,t+1为水库相邻时段不同径流值之间转移的相关系数，Ns为样本的总年份数，
Figure BDA0002263031370000077
为n水库t时段和t+1时段第i年的实际入库径流值，
Figure BDA0002263031370000078
为n水库t时段和t+1时段第i年的实际入库径流均值。σt为t时段和t+1时段实际入库径流的均方差。
求解出t时段与t+1时段之间的相关系数，并判定其马尔科夫性，并求出第t时段与t+1时段之间的状态转移矩阵。按照频率将入库径流划分为N个区段，有pjk＝P(St+1|St)，其中，j,k分别表示第t和t+1时段入库径流序列所在的区段，第t个时段状态转移概率矩阵形式为：
Figure BDA0002263031370000081
式中，Pt为状态转移概率矩阵，pi,j为模拟径流随机值由i状态转移到j状态的概率值
S2：根据各水库约束条件中设置的水位的上下限值，按照指定的步长s1从下限值到上限值进行等距离散化处理，得出离散化的水位集合S作为SARSA算法中的状态集合。同时对拟合出来的随机入库流量值也按照一定的步长k1进行离散化处理，得出的离散化的入库径流集合K作为SARSA算法中的入库流量状态集合。
S3：利用SARSA算法进行求解，包括以下步骤：
S3.1：(1)、确定最后一个时段的AAPFD期望值。最后一个阶段的余留效益取0。该时段N种代表入库流量下所对应的AAPFD期望值即为本阶段期望值。
(2)、确定倒数第二时段的AAPFD期望值。假设本时段初末水位均离散为M种状态，且该时段入库流量同样取N种。先求出每一种代表流量下本时段初末状态水位的阶段AAPFD期望值，然后再根据状态转移概率矩阵可得倒数第二时段在N种代表流量下对应于最后一个时段的余留效益值，最后根据求倒数第二时段的累计AAPFD期望值＝阶段AAPFD期望值+余留效益值。
(3)、后续时段直至初始时段的累计AAPFD期望值，均采用上述方法进行求解。
S3.2：设各时段初离散水位状态集为S，时段出库流量的动作集为A，学习率α，衰减因子γ，以及贪婪决策的探索率ε。设置每组状态动作集合[Vn,i,Vn,j]所对应的阶段AAPFD期望值为奖励矩阵R。设置初始Q值矩阵Q(s,a)和迭代轮数n。
S3.3：设置梯级各水库的起调水位，随机产生一个0-1之间的均匀分布小数ξ。使用贪婪决策，若ξ≤ε则采取探索策略，在动作集A中随机选择一个动作a为当前状态s选择的动作。反之则在Q值矩阵Q(s,a)中选择状态s对应的最大Q值，然后通过Q值确定下一动作a。
S3.4：通过当前状态动作集合(s,a)，得出在始末水位状态下本时段采取动作时，所获得的奖励值R(s,a)。将当前状态s采取动作a获得的新状态s′赋值给s。再次使用贪婪决策随机产生一个0-1之间的均匀分布小数ξ，若ξ≤ε则采取探索策略，在动作集A中随机选择一个动作a′作为当前状态s选择的动作。根据时序差分法(Temporal Difference,TD)通过下一时刻的价值更新前一时刻的价值的思想对值函数方程进行如下计算：
Figure BDA0002263031370000091
式中：
Figure BDA0002263031370000092
为t+1时段采用策略π得到状态动作对(s,a)时的Q值函数，
Figure BDA0002263031370000093
为t时段采用策略π得到状态动作对(s,a)时的Q值函数，α为SARSA算法的学习率，
Figure BDA0002263031370000094
为采取动作a后状态由s转移到s′的回报值。γ为SARSA算法的奖赏折扣率，
Figure BDA0002263031370000095
为t时段采用策略π得到状态动作对(s′,a′)时的Q值函数。通过对Q值函数进行不断的试错更新，最终形成一个以s，a分别为横纵坐标，以离散的水位数为维数的Q矩阵。
S3.5：将下一阶段状态s′赋值给s，下一阶段的a′赋值给a，若此时的s为目标状态，则此轮迭代完成，否则转移到步骤S3.3,直到迭代次数完成，或Q(s,a)收敛。
S3.6：根据已经收敛的Q(s,a)，通过对Q矩阵进行最大值逐步寻优，得出最优的策略π。进而得出从初始水位开始到最末水位的梯级水库修正全年流量偏差(AAPFD)期望值，以及各水库各时段的最优出库流量
Figure BDA0002263031370000096

Claims (6)

1.一种梯级水库生态随机优化调度模型，其特征在于，该模型的目标函数为：
Figure FDA0002263031360000011
式中：f*为梯级水库最小AAPFD期望值的平方，N为梯级水电站的个数，T为调度周期总的时段数，Ns为入库径流离散状态数，Rn,t为第n个水库第t时段的指标函数，fn,t为第n个水库第t时段的AAPFD值；
Rn,t、fn,t表达式如下所示：
Figure FDA0002263031360000012
式中：Vn,i为第n级电站i状态时对应的水位，
Figure FDA0002263031360000013
为第n级电站第t时段的出库流量，
Figure FDA0002263031360000014
为第n个水库第t时段离散的j个入库流量值，j、k为入库径流相互之间的状态转移，
Figure FDA0002263031360000015
为第n个水库第t时段的出库流量，
Figure FDA0002263031360000016
为第n个水库第t时段的入库径流，
Figure FDA0002263031360000017
为第n个水库第t时段处于j状态转移到k状态的转移概率，
Figure FDA0002263031360000018
为调度周期天然入库流量的平均值。
2.根据权利要求1所述一种梯级水库生态随机优化调度模型，其特征在于：所述目标函数约束条件为：
①、龙头水库水量平衡约束：
Figure FDA0002263031360000019
上式中：V1 t+1为龙头水库第t+1时段的水位，V1 t为龙头水库第t时段的水位，
Figure FDA00022630313600000110
为龙头水库第t+1时段的入库流量，
Figure FDA00022630313600000111
为龙头水库第t+1时段的出库流量，△t为面临时段的时长；
②、其他水库水量平衡约束：
Figure FDA00022630313600000112
上式中：
Figure FDA00022630313600000113
为第n+1级水库第t+1时段的水位，
Figure FDA00022630313600000114
为第n+1级水库第t时段的水位，
Figure FDA0002263031360000021
为第n级水库第t时段的出库流量，
Figure FDA0002263031360000022
为第t时段第n+1级水库上游区间来水流量，
Figure FDA0002263031360000023
为第n+1级水库第t时段的出库流量，△t为面临时段的时长，N为梯级水库的总数；
③、水库蓄水量约束：
Figure FDA0002263031360000024
上式中：Vn,min为第n级水库的最小蓄水量，
Figure FDA0002263031360000025
为第n级水库第t时段的水位，Vn,max为第n级水库的最大蓄水量，N为梯级水库的总数；
④、水库下泄流量约束：
Figure FDA0002263031360000026
上式中：Qn,min,ck为第n级水库最小下泄流量，
Figure FDA0002263031360000027
为第n级水库第t时段的下泄流量，Qn,max,ck为第n级水库最大下泄流量，N为梯级水库的总数；
⑤、保证出力约束：Wn,bz≤Wn≤Wn,zj n＝1,...,N；
上式中：Wn,bz为第n级水库的保证出力，Wn为第n级水库的出力，Wn,jz为第n级水库的机组最大出力；N为梯级水库的总数。
3.根据权利要求1所述一种梯级水库生态随机优化调度模型，其特征在于：所述优化调度模型中，对梯级电站按照上下游的顺序从上到下依次编号为1,...,n，并根据上下游水库之间的水量关系，建立梯级各水库之间的水量平衡方程：
Figure FDA0002263031360000028
上式中：Ln+1,t为流域内t时段两级水库之间的区间来水。
4.根据权利要求1所述一种梯级水库生态随机优化调度模型，其特征在于：所述优化调度模型中，由于水库中长期径流具有随机性，因此各个水库的入库径流采用随机入库过程作为梯级水库入库径流的表示方式，即根据皮尔逊Ⅲ型分布法计算所得频率离散为从小到大的各水库各时段随机来水的拟合值。
5.根据权利要求1所述一种梯级水库生态随机优化调度模型，其特征在于：所述优化调度模型中，生态指标利用修正全年流量偏差AAPFD指标，该指标反映一个调度周期内流量变化对河流生态环境的影响，其值越小，表示水库调度后流量变化对河流生态系统的影响越小，河流生态环境越好；
对于单一水库，修正全年流量偏差AAPFD表示为：
Figure FDA0002263031360000031
式中：T为调度周期总时段数，
Figure FDA0002263031360000032
为t时段水库出库流量值；
Figure FDA0002263031360000033
为t时段水库天然入库流量；
Figure FDA0002263031360000034
为调度周期天然入库流量的平均值。
6.采用SARSA算法对权利要求1-5中所述任意一种梯级水库生态调度随机优化模型进行优化计算，得到修正全年流量偏差AAPFD期望值的平方最大，以及各水库各时段的最优出库流量
Figure FDA0002263031360000035
包括以下步骤：
步骤1：根据梯级水库历年入库流量，求解各水库各时段的皮尔逊Ⅲ型概率密度曲线的统计参数，并根据统计样本参数绘制各水库各时段的径流频率曲线，最终通过皮尔逊Ⅲ型分布模拟得出的径流随机值计算状态转移概率矩阵；
步骤2：对每一个水库任一时段始末库容分别为Vn,i,Vn,j时，进行阶段收益求解，通过状态转移矩阵求解和Bellman方程求解在每一个水库在下一时刻水位随机情况下的余留收益，从而得到强化学习算法中的收益矩阵；
步骤3：根据SARSA算法试错的原理，从第一阶段梯级水库起调水位V0开始进行尝试，对每一组[Vi,Vj]均对出库寻优计算，直到进行到最后一个阶段，完成一次寻优，然后根据时序差分法TD通过下一时刻的价值更新前一时刻的价值的思想对值函数方程，进行如下计算：
Figure FDA0002263031360000036
式中：
Figure FDA0002263031360000037
为t+1时段采用策略π得到状态动作对(s,a)时的Q值函数，
Figure FDA0002263031360000038
为t时段采用策略π得到状态动作对(s,a)时的Q值函数，α为SARSA算法的学习率，
Figure FDA0002263031360000039
为采取动作a后状态由s转移到s′的回报值；γ为SARSA算法的奖赏折扣率，
Figure FDA00022630313600000310
为t时段采用策略π得到状态动作对(s′,a′)时的Q值函数；
通过对Q值函数进行不断的试错更新，最终形成一个以s，a分别为横纵坐标，以离散的水位数为维数的Q矩阵；
步骤4：重复对当前状态Vi进行观察，并根据贪婪决策选择下一阶段的执行动作，然后对当前Q矩阵进行更新，直到Q矩阵收敛；然后从Q矩阵中选取最大值作为策略，即可得到修正全年流量偏差(AAPFD)期望值的平方最大，以及各水库各时段的最优出库流量
Figure FDA0002263031360000041
CN201911077800.2A 2019-11-06 2019-11-06 一种梯级水库生态随机优化调度模型及求解方法 Active CN110880048B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911077800.2A CN110880048B (zh) 2019-11-06 2019-11-06 一种梯级水库生态随机优化调度模型及求解方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911077800.2A CN110880048B (zh) 2019-11-06 2019-11-06 一种梯级水库生态随机优化调度模型及求解方法

Publications (2)

Publication Number Publication Date
CN110880048A CN110880048A (zh) 2020-03-13
CN110880048B true CN110880048B (zh) 2022-06-21

Family

ID=69729092

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911077800.2A Active CN110880048B (zh) 2019-11-06 2019-11-06 一种梯级水库生态随机优化调度模型及求解方法

Country Status (1)

Country Link
CN (1) CN110880048B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111461419B (zh) * 2020-03-25 2023-05-30 长江水利委员会长江科学院 一种通航河流水库水沙联合调度方案生成方法
CN111626512B (zh) * 2020-05-28 2023-11-14 国家电网公司西南分部 一种考虑电网断面约束的梯级水电日前调度方法
CN111637435B (zh) * 2020-06-05 2021-09-28 沈阳航空航天大学 基于sarsa的核动力系统蒸汽发生器水位控制方法
CN111950982B (zh) * 2020-08-10 2024-04-19 中国电力科学研究院有限公司 梯级水电水平衡模型构建方法、系统、设备及存储介质
CN112036633B (zh) * 2020-08-26 2021-06-04 三峡大学 一种基于水库生态发电多目标中长期随机调度模型的优化调度方法
CN112036632B (zh) * 2020-08-26 2021-06-04 三峡大学 一种基于梯级水库生态发电多目标中长期随机调度模型的优化调度方法
CN112054561B (zh) * 2020-09-16 2022-06-14 三峡大学 基于SARSA(λ)算法的风电-抽蓄联合系统日随机动态调度方法
CN113191543B (zh) * 2021-04-27 2022-07-29 西安理工大学 考虑水资源供需过程关联性与随机性的水库优化调度方法
CN113110067A (zh) * 2021-05-20 2021-07-13 国网湖南省电力有限公司 一种小水电站的水头优化方法、装置、设备及介质
CN113487249B (zh) * 2021-09-07 2021-12-07 长江水利委员会水文局 一种自适应的水电站智能生态调控方法
CN115630800B (zh) * 2022-09-22 2023-07-04 广东省水利水电科学研究院 水利枢纽防洪发电联合优化调度方法、系统、装置及存储介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO20022204D0 (no) * 2001-05-09 2002-05-08 Inst Francais Du Petrole Fremgangsmåte for å bestemme ved numerisk simulering av restitueringsbetingelser, av fluidene i reservoaret, av enkompleks brönn skadet under boreoperasjoner
CA2435626A1 (fr) * 2003-07-28 2005-01-28 Benoit Godin Methode d`echantillonnage et de gestion du risque lors de l`inspection metallurgique de conduits ou reservoirs
CN1672499A (zh) * 2005-03-29 2005-09-28 西北农林科技大学 新低压滴灌系统滴头工作压力与毛管优化设计方法
CN101714186A (zh) * 2009-12-29 2010-05-26 北京师范大学 兼顾人类和生态需求的供水型水库调度图优化确定方法
CN102080366A (zh) * 2011-01-01 2011-06-01 国网电力科学研究院 一种梯级水库联合调度图制作方法
CN102817335A (zh) * 2012-07-12 2012-12-12 贵州乌江水电开发有限责任公司 一种梯级水库群联合防洪优化调度的方法及系统
CN104123589A (zh) * 2014-06-24 2014-10-29 华中科技大学 一种梯级水电站短期优化调度方法
EP2819086A1 (en) * 2013-06-28 2014-12-31 Repsol, S.A. Method to generate scenarios of hydrocarbon reservoirs based on limited amount of information on a target hydrocarbon reservoir
CN108734415A (zh) * 2018-06-06 2018-11-02 广西桂冠电力股份有限公司 一种梯级水库的入库流量自动校正系统
CN109886473A (zh) * 2019-01-24 2019-06-14 河海大学 一种考虑下游生态的流域风光水系统多目标优化调度方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6823297B2 (en) * 2003-03-06 2004-11-23 Chevron U.S.A. Inc. Multi-scale finite-volume method for use in subsurface flow simulation
US20100286917A1 (en) * 2009-05-07 2010-11-11 Randy Doyle Hazlett Method and system for representing wells in modeling a physical fluid reservoir

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NO20022204D0 (no) * 2001-05-09 2002-05-08 Inst Francais Du Petrole Fremgangsmåte for å bestemme ved numerisk simulering av restitueringsbetingelser, av fluidene i reservoaret, av enkompleks brönn skadet under boreoperasjoner
CA2435626A1 (fr) * 2003-07-28 2005-01-28 Benoit Godin Methode d`echantillonnage et de gestion du risque lors de l`inspection metallurgique de conduits ou reservoirs
CN1672499A (zh) * 2005-03-29 2005-09-28 西北农林科技大学 新低压滴灌系统滴头工作压力与毛管优化设计方法
CN101714186A (zh) * 2009-12-29 2010-05-26 北京师范大学 兼顾人类和生态需求的供水型水库调度图优化确定方法
CN102080366A (zh) * 2011-01-01 2011-06-01 国网电力科学研究院 一种梯级水库联合调度图制作方法
CN102817335A (zh) * 2012-07-12 2012-12-12 贵州乌江水电开发有限责任公司 一种梯级水库群联合防洪优化调度的方法及系统
EP2819086A1 (en) * 2013-06-28 2014-12-31 Repsol, S.A. Method to generate scenarios of hydrocarbon reservoirs based on limited amount of information on a target hydrocarbon reservoir
CN104123589A (zh) * 2014-06-24 2014-10-29 华中科技大学 一种梯级水电站短期优化调度方法
CN108734415A (zh) * 2018-06-06 2018-11-02 广西桂冠电力股份有限公司 一种梯级水库的入库流量自动校正系统
CN109886473A (zh) * 2019-01-24 2019-06-14 河海大学 一种考虑下游生态的流域风光水系统多目标优化调度方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"三峡及金沙江下游梯级水库群蓄水期联合调度策略";刘强 等;《南水北调与水利科技》;20161031;全文 *
"基于SARSA算法的水库长期随机优化调度研究";李文武 等;《水电能源科学》;20180930;全文 *
"基于随机动态规划的混合式抽水蓄能电站水库中长期优化调度研究";李文武 等;《电力系统保护与控制》;20130501;全文 *

Also Published As

Publication number Publication date
CN110880048A (zh) 2020-03-13

Similar Documents

Publication Publication Date Title
CN110880048B (zh) 一种梯级水库生态随机优化调度模型及求解方法
CN110851977B (zh) 基于生态流量的供水-发电-生态多目标调度图优化方法
CN112036633B (zh) 一种基于水库生态发电多目标中长期随机调度模型的优化调度方法
Liu et al. Deriving optimal refill rules for multi-purpose reservoir operation
CN103049671B (zh) 一种自适应气候变化的水库多目标优化调度图编制方法
CN109670650B (zh) 基于多目标优化算法的梯级水库群调度模型的求解方法
CN103088784B (zh) 一种梯级水库汛限水位实时动态控制方法
CN112036632B (zh) 一种基于梯级水库生态发电多目标中长期随机调度模型的优化调度方法
CN104715292A (zh) 基于最小二乘支持向量机模型的城市短期用水量预测方法
CN109345068B (zh) 一种基于余留期效益函数近似的水电站水库两阶段随机优化调度方法
CN107563538B (zh) 大电网平台下的面向关键水位控制的多目标水库群调度优化方法
CN105243502A (zh) 一种基于径流区间预测的水电站调度风险评估方法及系统
CN111080157B (zh) 一种梯级水电站排磷量的调度方法和系统
CN109636054A (zh) 基于分类和误差组合预测的太阳能光伏发电量预测方法
CN104504455B (zh) 一种梯级蓄能控制下水电站群长期优化调度方法
CN112700080A (zh) 一种梯级水电多阶段优化调度方法
CN115271304A (zh) 基于气候变化的跨流域调水工程优化调度方法
CN108537370A (zh) 基于混合智能降维算法的特大流域水电站群优化调度方法
CN113363976B (zh) 一种基于场景图的风光水互补发电系统中期优化调度方法
CN110633849A (zh) 一种基于启发式大系统分解协调思想的水资源配置方法
CN108330896B (zh) 一种梯级水库调度方法及系统
CN112132379B (zh) 考虑经济性的新能源跨区域消纳评估方法和存储介质
CN109214610A (zh) 一种基于长短期记忆神经网络的饱和电力负荷预测方法
CN111476477A (zh) 基于发电效益目标的梯级水电站中长期优化调度方法
CN115081681B (zh) 一种基于prophet算法的风电功率预测方法

Legal Events

Date Code Title Description
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