发明名称 一种协调水电站汛期防洪风险与发电效益的调度方法
摘要 本发明公开了一种协调水电站汛期防洪与发电效益的调度方法,包括:步骤1,将水电站汛期调度概化为防洪与发电双目标两阶段协调控制问题;步骤2,分析预报信息不确定性及其可能产生的附加防洪风险;步骤3,以两阶段的发电效益最大和防洪风险最小为优化目标建立优化调度模型;步骤4,分析目标函数的经济学特性,给出防洪与发电二者的最佳权衡条件,进而求解模型。本发明考虑水文预报信息及其不确定性影响建立水电站汛期双目标两阶段优化调度模型,充分利用了水文预报信息,依据模型求解结果进行汛期调度,可以最大化水电站发电和防洪综合效益。此外本发明依据KKT条件解释与分析模型最优解,大大降低了模型的求解难度。
申请公布号 CN105869065A 申请公布日期 2016.08.17
申请号 CN201610173639.9 申请日期 2016.03.24
申请人 大连理工大学 发明人 彭勇;丁伟;张弛;周惠成;唐榕;黄朔
分类号 G06Q50/06(2012.01)I 主分类号 G06Q50/06(2012.01)I
代理机构 大连理工大学专利中心 21200 代理人 梅洪玉;潘迅
主权项 一种协调水电站汛期防洪风险与发电效益的调度方法,其特征在于,包括以下步骤:第一步,将水电站汛期调度概化为防洪与发电双目标两阶段的协调控制问题;所述的两阶段指阶段1决策阶段和阶段2余留阶段;第二步,分析决策阶段预报信息和余留阶段预报信息的不确定性,阶段1按确定性预报处理,若下游组合流量不超过其安全泄水量,则不考虑防洪风险,阶段2按不确定性预报处理,考虑防洪风险;量化预报信息不确定性产生的防洪风险,得到预报信息不确定性产生的防洪风险率计算式,具体步骤为:2.1)由水量平衡方程建立水电站泄水与来水之间的关系式;r<sub>t+1</sub>=s<sub>t+1</sub>+q<sub>t+1</sub>‑s<sub>t+2</sub>        (1)其中,r<sub>t+1</sub>为余留阶段的水库泄流,s<sub>t+1</sub>为阶段1末库容,q<sub>t+1</sub>为余留阶段的来水,s<sub>t+2</sub>为余留阶段的末库容;2.2)基于防洪风险率概念,得到防洪风险率计算式;<maths num="0001"><math><![CDATA[<mrow><mi>R</mi><mi>i</mi><mi>s</mi><mi>k</mi><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>P</mi><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&gt;</mo><msub><mi>r</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><msub><mi>r</mi><mi>max</mi></msub><mrow><mo>+</mo><mi>&infin;</mi></mrow></msubsup><mi>f</mi><mrow><mo>(</mo><msub><mi>r</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><msub><mi>dr</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000949303920000011.GIF" wi="1342" he="103" /></maths>其中,Risk(r<sub>t+1</sub>)为利用预报信息后的下游防洪风险率,r<sub>max</sub>为下游安全泄量,P(r<sub>t+1</sub>>r<sub>max</sub>)为泄流r<sub>t+1</sub>超过下游安全泄量r<sub>max</sub>的概率,f(r<sub>t+1</sub>)为泄流r<sub>t+1</sub>的概率密度函数;2.3)当q<sub>t+1</sub>≤s<sub>t+2</sub>‑s<sub>t+1</sub>‑r<sub>max</sub>时,r<sub>t+1</sub>>r<sub>max</sub>,得到预报信息不确定性产生的防洪风险率计算式;<maths num="0002"><math><![CDATA[<mrow><mi>R</mi><mi>i</mi><mi>s</mi><mi>k</mi><mrow><mo>(</mo><msub><mi>s</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>P</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&le;</mo><msub><mi>s</mi><mrow><mi>t</mi><mo>+</mo><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>s</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>r</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mi>&infin;</mi></mrow><mrow><msub><mi>s</mi><mrow><mi>t</mi><mo>+</mo><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>s</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>r</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow></msubsup><mrow><mi>f</mi><mrow><mo>(</mo><msub><mi>q</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><msub><mi>dq</mi><mrow><mi>t</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000949303920000012.GIF" wi="1517" he="95" /></maths>其中,Risk(s<sub>t+1</sub>)为利用预报信息后的下游防洪风险率,f(q<sub>t+1</sub>)为来水q<sub>t+1</sub>的概率密度函数,P(q<sub>t+1</sub>≤s<sub>t+2</sub>‑s<sub>t+1</sub>‑r<sub>max</sub>)为泄流超过下游安全泄量r<sub>max</sub>的概率;第三步,以两阶段的发电效益最大和防洪风险最小为优化目标建立优化调度模型;所述的优化调度模型的目标函数为:<img file="FDA0000949303920000021.GIF" wi="1716" he="166" />其中,G<sub>1</sub>为发电目标,G<sub>2</sub>为风险率目标,ω为发电目标权重,N<sub>装</sub>为装机容量,T为预报信息可利用预见期,Δt为一个时段,g<sub>t</sub>(s<sub>t</sub>,s<sub>t+1</sub>)为阶段1的发电量,Eg<sub>t+1</sub>(s<sub>t+1</sub>,s<sub>t+2</sub>)为阶段2的发电量,Risk(s<sub>t+1</sub>)为利用预报信息后的下游防洪风险率;第四步,分析目标函数的经济学特性,给出防洪与发电二者的最佳权衡条件,运用Karush–Kuhn–Tucker条件,求解优化调度模型的目标函数的最优解;当所有约束条件均能满足时,蓄水量带来的发电量增加值和风险率增加值相等,则该蓄水量为最优蓄水量,符合边际效益原则;否则针对不满足的约束进行反推,得到最优蓄水量。
地址 116024 辽宁省大连市甘井子区凌工路2号