发明名称 水电站发电出力方案制作与最优负荷分配耦合建模方法
摘要 本发明提供一种水电站发电出力方案制作与最优负荷分配耦合建模方法,包括:以调峰容量效益最大为目标建立水电站发电出力方案模型,以耗水量最小为目标建立厂内机组负荷最优分配模型;求解水电站发电出力方案模型,得到水电站初始出力方案;将水电站初始出力方案作为所述厂内机组负荷最优分配模型的输入,对厂内机组负荷最优分配模型进行模拟计算,通过一体化自适应调整模式对水电站初始出力方案进行算法循环迭代和模型轨迹嵌套修正,耦合后的输出即为水电站最终出力方案。实现具有自适应调整特性的水电站发电出力方案制作与厂内机组负荷最优分配一体化无缝耦合运作,获得满足厂内机组负荷最优分配条件下的精细的水电站发电出力方案。
申请公布号 CN103745023B 申请公布日期 2015.08.12
申请号 CN201310601268.6 申请日期 2013.11.22
申请人 华中科技大学;中国长江三峡集团公司 发明人 周建中;肖舸;莫莉;陈建国;严冬;李晖;张睿;陈忠贤;闫宝伟;刘志武;欧阳硕;张慧;王学敏;李鹏;王超;梁君;李超顺;阮燕云;叶磊;吴江;李纯龙;卢鹏;袁柳;朱双;张华杰;廖想;孟长青;谢蒙飞;吉鹏;丁小玲;张海荣;王婷婷;卢韦伟;江兴稳;王华为;牛广利;张德发
分类号 G06F17/50(2006.01)I;G06Q50/06(2012.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京市盛峰律师事务所 11337 代理人 赵建刚
主权项 一种水电站发电出力方案制作与最优负荷分配耦合建模方法,其特征在于,包括以下步骤:S1,以调峰容量效益最大为目标建立水电站发电出力方案模型,以耗水量最小为目标建立水电站厂内机组负荷最优分配模型;其中,所述水电站发电出力方案模型通过以下方式建立:根据水电站预报来水、机组检修安排情况、电站调度规程以及其它水能综合利用要求,结合电网给定的典型负荷曲线形式,以调峰容量效益最大为目标,建立所述水电站发电出力方案模型;该水电站发电出力方案模型的目标函数为水电站调峰容量效益最大,即<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>F</mi><mo>=</mo><mi>Max</mi><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mrow><mo>(</mo><msub><mi>Q</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>,</mo><msub><mi>H</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>&Delta;T</mi><mo>&CenterDot;</mo><msub><mi>&beta;</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA0000711988250000011.GIF" wi="680" he="149" /></maths>式中,F为控制期内梯级电站调峰容量效益,万kW;N为电站个数;T为调度期时段数;ΔT为时段长;P<sub>i,t</sub>为电站i在t时段的出力;Q<sub>i,t</sub>为电站i在t时段的发电流量;H<sub>i,t</sub>为电站i在t时段的平均水头;β<sub>i,t</sub>为电站i在t时段的调峰权重参数;其中,所述厂内机组负荷最优分配模型通过以下方式建立:根据水电站水位控制及来水情况,建立以耗水量最小为目标的厂内机组负荷最优分配模型;该厂内机组负荷最优分配模型的目标函数是耗水量最小,即<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>W</mi><mi>i</mi></msub><mo>=</mo><mi>min</mi><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mo>{</mo><msub><mi>q</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>[</mo><msub><mi>h</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>,</mo><msub><mi>N</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>]</mo><mo>&CenterDot;</mo><mi>&Delta;T</mi><mo>+</mo><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>,</mo><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&CenterDot;</mo><msub><mi>q</mi><mrow><mi>i</mi><mo>,</mo><mi>sk</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>,</mo><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>q</mi><mrow><mi>i</mi><mo>,</mo><mi>ck</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>}</mo><mo>;</mo></mrow>]]></math><img file="FDA0000711988250000012.GIF" wi="1515" he="155" /></maths>式中,W<sub>i</sub>为电站i在调度期内给定负荷任务时的总耗水量;K为电站机组台数;T为调度期时段数;ΔT为时段长;N<sub>i,k,t</sub>为机组k在t时段的出力;h<sub>i,k,t</sub>为机组k在t时段的净水头;q<sub>i,k,t</sub>为t时段机组k在净水头h<sub>i,k,t</sub>下出力为N<sub>i,k,t</sub>时的发电引用流量;u<sub>i,k,t</sub>为t时段机组k的开停机状态变量;S2,以水电站空间最优流量分配表为基础,在第一约束函数集的约束下求解所述水电站发电出力方案模型,得到水电站初始出力方案;弃水模式下的水电站初始出力方案制作:弃水模式下电站各时段均满发,出力过程不考虑电网特性要求:①设初始时段t=1,计算相应水头下电站满发流量,若出库小于满发流量,降低时段末水位,加大出库流量;反之,则升高时段末水位,减小出库流量,迭代求解至出库流量等于满发流量或时段水位到达极限为止;②令t=t+1,重复使用与①同样的方法进行满发情况下时段水位计算直至t=T;③判断调度期末水位是否满足给定控制蓄水位要求,若不满足则相应的增加或减小各时段出库,将差值水量平均分配到各个时段,重新进行水量平衡计算,直至满足调度期末水位约束为止;非弃水模式下水电站初始出力方案制作:非弃水模式下,水电站以充分发挥电站容量调峰效益为目标,出力应与电网负荷曲线趋势保持一致:①由上一级电站下泄、区间入流预报以及中期分配可用水量估算电站调度期发电能力,并转换为电站平均出力N;②以给定电网负荷曲线形式为基础,由P<sub>i,t</sub>=N·C<sub>t</sub>经电力电量平衡计算电站初始出力过程,并由时段出力计算调峰权重参数β<sub>i,t</sub>,其中C<sub>t</sub>为电网负荷曲线峰段系数,N为不同电网供电比例;③设初始时段t=1,由时段电站负荷、入库流量计算时段末水位、出库流量和机组组合;然后令t=t+1,重复由时段电站负荷、入库流量计算时段末水位、出库流量和机组组合直至t=T;④根据出力过程进行调度期水量优化计算,若末水位不满足给定蓄水位要求,根据β<sub>i,t</sub>调整出力至调度期末水位满足要求,具体为:如果水位偏高,则优先加大β<sub>i,t</sub>大的时段出力,反之,如果水位偏低,则优先减小β<sub>i,t</sub>小的时段出力,获得改进的电站运行方案;或者,查找权重系数最小的时段与最大的时段,并按一定步长分别减小与增大出力,获得改进的电站运行方案,计算β<sub>i,t</sub>并循环执行本步骤,直到末水位满足给定蓄水位要求,然后执行5;5,若搜索次数到达上限或目标值不再改变则停止,输出出力制作结果;否则,查找最小的时段与最大的时段,并按一定步长分别减小出力与增大出力,根据给定目标获取改进的电站运行方案,然后转到步骤③;其中,所述水电站发电出力方案模型的第一约束函数集为(1.4.1)—(1.4.3)、(1.4.5)—(1.4.8);另外,所述厂内机组负荷最优分配模型具有第二约束函数集为(1.4.1)—(1.4.2)、(1.4.4)—(1.4.8);(1.4.1)水库水力联系:I<sub>i,t</sub>=Q<sub>i‑1,t‑τ</sub>+S<sub>i‑1,t‑τ</sub>+R<sub>i,t</sub>式中,I<sub>i,t</sub>为电站i在t时段的入库流量;S<sub>i‑1,t‑τ</sub>为第i‑1个水电站在t‑τ时段的弃水流量;τ为i‑1与i电站间水流时滞;R<sub>i,t</sub>为i‑1与i电站之间的区间入流;Q<sub>i‑1,t‑τ</sub>为电站i‑1在t‑τ时段的发电流量;(1.4.2)水量平衡约束:V<sub>i,t</sub>=V<sub>i,t‑1</sub>+(I<sub>i,t</sub>‑Q<sub>i,t</sub>‑S<sub>i,t</sub>)·Δt式中,V<sub>i,t</sub>为t时段末水库i的蓄水量;(1.4.3)电站库容/流量/出力约束:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>VL</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&le;</mo><msub><mi>V</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&le;</mo><msub><mi>VU</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>QL</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&le;</mo><mrow><mo>(</mo><msub><mi>Q</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>S</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>)</mo></mrow><mo>&le;</mo><msub><mi>QU</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>PL</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&le;</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&le;</mo><msub><mi>PU</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000711988250000031.GIF" wi="490" he="219" /></maths>式中,VU<sub>i,t</sub>与VL<sub>i,t</sub>分别为t时段水库i蓄水容量上下限;QU<sub>i,t</sub>与QL<sub>i,t</sub>分别为t时段电站i下泄流量上下限;PU<sub>i,t</sub>与PL<sub>i,t</sub>分别为t时段电站i出力上下限;(1.4.4)末水位控制约束:Z<sub>i,t</sub>=Z<sub>i,end</sub>式中,Z<sub>i,t</sub>为电站i调度期末水位;Z<sub>i,end</sub>为电站i调度期末水位控制值;(1.4.5)单站出力/水位/流量变幅约束:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&le;</mo><msub><mi>PCH</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Z</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><msub><mi>Z</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&le;</mo><msub><mi>ZCH</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><msub><mi>Q</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>-</mo><msub><mi>Q</mi><mrow><mi>i</mi><mo>,</mo><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>&le;</mo><msub><mi>QCH</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000711988250000032.GIF" wi="337" he="201" /></maths>式中,PCH<sub>i</sub>、ZCH<sub>i</sub>、QCH<sub>i</sub>分别为电站i在t时段允许的最大出力变幅、水位变幅以及流量变幅;(1.4.6)单站负荷平衡约束:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>Load</mi><mi>i</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>N</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>u</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>,</mo><mi>t</mi></mrow></msub></mrow>]]></math><img file="FDA0000711988250000033.GIF" wi="361" he="133" /></maths>式中,Load<sub>i</sub>表示电网调度中心下达给电站i的负荷任务值;(1.4.7)机组稳定运行约束:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>N</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi><mo>,</mo><mi>t</mi></mrow></msub><mo>&NotElement;</mo><msub><mi>MV</mi><mrow><mi>i</mi><mo>,</mo><mi>k</mi></mrow></msub></mrow>]]></math><img file="FDA0000711988250000036.GIF" wi="238" he="77" /></maths>式中,NU<sub>i,k</sub>与NL<sub>i,k</sub>分别为i电站k号机组出力上下限;(1.4.8)机组最短开停机时间约束:<maths num="0007" id="cmaths0007"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>off</mi></mrow><mi>t</mi></msubsup><mo>&GreaterEqual;</mo><msub><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>down</mi></mrow></msub></mtd></mtr><mtr><mtd><msubsup><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>on</mi></mrow><mi>t</mi></msubsup><mo>&GreaterEqual;</mo><msub><mi>T</mi><mrow><mi>i</mi><mo>,</mo><mi>up</mi></mrow></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000711988250000034.GIF" wi="254" he="164" /></maths>式中,T<sub>i,up</sub>、T<sub>i,down</sub>分别为机组k允许的最短开、停机时间限制;<img file="FDA0000711988250000035.GIF" wi="222" he="86" />分别为机组k在t‑1时段以前的持续开、停机历时;S3,判断所研究的水电站在调度期是否处于弃水模式,如果为弃水模式,则所述水电站初始出力方案即为水电站最终出力方案,直接输出所述水电站最终出力方案,并结束流程;如果为非弃水模式,则执行S4;S4,将所述水电站初始出力方案作为所述厂内机组负荷最优分配模型的输入,对所述厂内机组负荷最优分配模型进行模拟计算,然后通过一体化自适应调整模式对所述水电站初始出力方案进行算法循环迭代和模型轨迹嵌套修正,耦合后的输出即为兼顾调峰容量收益最大和耗水量最小为目标的水电站最终出力方案;其中,S2中,所述电站空间最优流量分配表通过以下方式制作:S2.1,根据电站运行水头范围,以一定步长对电站毛水头H及电站发电流量Q进行离散;选择一组电站毛水头H及电站发电流量Q离散组合作为当前输入条件;S2.2,根据电站毛水头H查找各类型机组稳定运行区间以及最大最小出力,确定机组各出力限制点,并由机组NHQ曲线查找各出力限制点对应的机组可运行流量范围,由动态规划最优化原理进行递推寻优,在可运行流量范围内进行机组间流量分配,得到在当前给定电站毛水头H及发电流量Q组合下的机组最优流量分配方案;S2.3,遍历求解所有电站毛水头H及发电流量Q组合下的最优流量分配方案,将优化结果集存入数据库,所述优化结果集即为所述电站空间最优流量分配表;其中,S4具体包括以下步骤:S4.1,将所述水电站初始出力方案作为所述厂内机组负荷最优分配模型的输入,对所述厂内机组负荷最优分配模型进行模拟计算,得到厂内机组负荷最优分配模拟结果,所述厂内机组负荷最优分配模拟结果具体包括:水电站相应出库流量过程、水库水位过程、最优机组组合及机组间负荷分配方案;S4.2,对于所述厂内机组负荷最优分配模拟结果,首先判断水电站时段出库流量是否满足最小出库流量限制;如果满足,则转至S4.3;如果不满足,则将该时段电站出力值加大,其他时段出力值保持不变,并根据调整后出力过程重新制作水电站发电出力方案,并转至S4.1,用重新制作的水电站发电出力方案替换所述水电站初始出力方案,循环进行;S4.3,进一步判断水电站调度期末水位是否满足给定调度期末水位约束;如果水电站调度期末水位大于给定调度期末水位,则加大电站出力,将可运用水量最大限度安排至用电高峰时段,对于按峰荷比例调峰方式,各峰段优先级按负荷序列排序,按峰荷比例增加峰段出力;若峰段均满发,则加大平段出力,然后根据调整后出力过程重新制作水电站发电出力方案,转至S4.1,重新进行厂内机组负荷最优分配模拟计算,算法循环迭代和模型轨迹嵌套调整出力直至满足水位约束为止,最终得到兼顾调峰容量收益最大和耗水量最小为目标的水电站最终出力方案。
地址 430074 湖北省武汉市洪山区珞喻路1037号