发明名称 一种面向多电网的水电站短期调峰特征值搜索方法
摘要 本发明公开了一种面向多电网的水电站短期调峰特征值搜索方法,属于水电能源优化运行和电力系统发电优化调度技术领域。本发明以各受电电网剩余负荷的最值作为特征值,在特征值限制范围内逐步调整特征值所在时段出力,逐步修正各电网受电过程以满足水电站运行限制要求。本发明根据调度期可用水量大小将多电网短期调峰分成非弃水模式与弃水模式两种调度模式。在非弃水模式下,以考虑电网典型负荷特性及调峰要求为主,在满足下游综合用水要求的前提下,尽量发挥水电站调峰容量效益,用电高峰多发电,用电低谷少发电;在弃水模式下,水电站满发运行,在给定电网分电比要求下,以剩余负荷方差最小为目标,通过经典粒子群算法制定电网受电计划。
申请公布号 CN104636831B 申请公布日期 2016.02.17
申请号 CN201510075639.0 申请日期 2015.02.12
申请人 华中科技大学 发明人 周建中;莫莉;严冬;刘懿;张勇传;李超顺;闫宝伟;曾小凡;梁藉;孙怀卫;陈璐;赵娜;卢鹏;王超;王学敏;李纯龙;丁小玲;王华为;牛广利;谢蒙飞;朱双;吴巍
分类号 G06F17/50(2006.01)I;G06Q10/04(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06F17/50(2006.01)I
代理机构 武汉东喻专利代理事务所(普通合伙) 42224 代理人 宋业斌
主权项 一种面向多电网的水电站短期调峰特征值搜索方法,其特征在于,包括:步骤1计算调度期可用水量,并估算水电站调度期的最大出力容量,根据平均下泄流量判断水电站是否有弃水,有则执行步骤7,否则执行步骤2;步骤2各电网以所述最大出力容量为上限,按各自的预测负荷曲线单独使用逐次切负荷法得到多条出力曲线并累加,获得水电站初始出力过程,并判断所述初始出力过程是否满足时段最大最小出力约束,若不满足则将不满足约束的部分累加起来作为出力调整步长ΔE,然后执行步骤3,否则执行步骤4;步骤3当水电站发电不足时,则执行子步骤(3‑1);当水电站发电过多时,则执行子步骤(3‑2):(3‑1)按所述出力调整步长ΔE增加出力,将所述出力调整步长ΔE按各电网分电比R<sup>g</sup>分为G份:ΔE<sup>g</sup>=R<sup>g</sup>ΔE,其中,ΔE<sup>g</sup>表示需要分配到g号电网的出力调整步长,g∈[1,G],并寻找各电网余荷最大的时段,并将ΔE<sup>g</sup>分配至该时段;若受约束的限制,单一时段不能将增加的电量完全消纳,则将剩余电量分配至已调整时段之外的余荷最大时段,不断重复上述过程,即将ΔE<sup>g</sup>不断分配到已调整时段之外的余荷最大时段,直到剩余电量为零则停止继续分配,重新计算水电站总出力,执行步骤4;(3‑2)按所述出力调整步长ΔE减小出力,按所述子步骤(3‑1)的方法,在满足保证出力约束的前提下,逐个遍历各受电电网,寻找余荷最小的时段减小出力,直至ΔE分配完毕;重新计算水电站总出力,执行步骤4;步骤4处理水电站出力最小持续时间约束和出力变幅约束,其中,处理所述出力最小持续时间约束时将最小持续时间Δt内的出力平均值作为所述最小持续时间Δt内的出力值,以在满足最小持续时间约束的同时保证出力总量不变;处理所述出力变幅约束时判断相邻两个时段出力P<sub>t‑1</sub>和P<sub>t</sub>是否满足如下式所示的出力变幅约束:|P<sub>t‑1</sub>‑P<sub>t</sub>|≤ΔP其中,ΔP为水电站时段允许最大出力变幅;步骤5根据水电站出力过程计算水电站各时段水位及下泄流量过程,并对违反下泄流量及运行水位约束的时段进行处理,判断调度期末计算水位是否满足控制末水位约束,若所述计算末水位不等于所述控制末水位,则执行步骤6;否则得到水电站出力及各电网受电结果,流程结束,其中,所述运行水位约束为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>Z</mi><mi>t</mi><mi>min</mi></msubsup><mo>&le;</mo><msub><mi>Z</mi><mi>t</mi></msub><mo>&le;</mo><msubsup><mi>Z</mi><mi>t</mi><mi>max</mi></msubsup></mrow>]]></math><img file="FDA0000761070000000021.GIF" wi="334" he="79" /></maths>所述下泄流量约束为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>Q</mi><mi>t</mi><mi>min</mi></msubsup><mo>&le;</mo><msub><mi>Q</mi><mi>t</mi></msub><mo>&le;</mo><msubsup><mi>Q</mi><mi>t</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msubsup></mrow>]]></math><img file="FDA0000761070000000022.GIF" wi="341" he="79" /></maths>时段水位/下泄流量变幅约束为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open = '{' close = ''><mtable><mtr><mtd><mo>|</mo><msub><mi>Z</mi><mi>t</mi></msub><mo>-</mo><msub><mi>Z</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>|</mo><mo>&le;</mo><mi>&Delta;</mi><mi>Z</mi></mtd></mtr><mtr><mtd><mo>|</mo><msub><mi>Q</mi><mi>t</mi></msub><mo>-</mo><msub><mi>Q</mi><mrow><mi>t</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>|</mo><mo>&le;</mo><mi>&Delta;</mi><mi>Q</mi></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000761070000000023.GIF" wi="300" he="134" /></maths>其中,<img file="FDA0000761070000000024.GIF" wi="232" he="79" />分别表示水电站t时段水位上下边界;<img file="FDA0000761070000000025.GIF" wi="211" he="62" />分别表示水电站t时段下泄流量上下边界;ΔZ、ΔQ分别表示水电站时段允许最大水位变幅和下泄流量变幅;步骤6若所述计算末水位大于所述控制末水位,则水电站发电不足,按所述出力调整步长ΔE,增加出力,转至执行所述步骤3;若所述计算末水位小于所述控制末水位,则水电站发电过大,按所述出力调整步长ΔE减小出力,转至执行所述步骤3;步骤7随机生成各电网受电过程,作为粒子的初始位置,处理最小持续时间约束,其中,粒子的初始位置按下式确定:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>t</mi><mi>g</mi></msubsup><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mi>r</mi><mo>&times;</mo><msup><mi>p</mi><mi>max</mi></msup></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mrow><mo>(</mo><mi>g</mi><mo>=</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>r</mi><mo>&times;</mo><mrow><mo>(</mo><mrow><msup><mi>p</mi><mi>max</mi></msup><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>g</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>P</mi><mi>t</mi><mi>i</mi></msubsup></mrow><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mrow><mo>(</mo><mi>g</mi><mo>&gt;</mo><mn>1</mn><mi>a</mi><mi>n</mi><mi>d</mi><mi> </mi><mi>g</mi><mo>&lt;</mo><mi>G</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>p</mi><mi>max</mi></msup><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>G</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>P</mi><mi>t</mi><mi>i</mi></msubsup></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mrow><mo>(</mo><mi>g</mi><mo>=</mo><mi>G</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000761070000000031.GIF" wi="909" he="439" /></maths>调整各电网受电量的出力调整步长计算公式为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><mi>&Delta;E</mi><mi>g</mi></msup><mo>=</mo><msup><mi>R</mi><mi>g</mi></msup><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msub><mi>P</mi><mi>t</mi></msub><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><msubsup><mi>P</mi><mi>t</mi><mi>g</mi></msubsup></mrow>]]></math><img file="FDA0000761070000000032.GIF" wi="486" he="142" /></maths>其中,r为0~1的随机数;P<sub>t</sub>表示t时段出力;P<sup>max</sup>表示时段最大出力;T表示一日内调度期总时段数;t∈[1,T]表示调度时段号,并通过经典粒子群算法以剩余负荷方差最小为目标进行寻优计算,若粒子位置不满足电网分电比约束,则通过所述步骤3中的方法对受电过程进行修正,当算法达到最大迭代次数时,计算终止,以粒子群全局最优解作为各电网受电计划结果。
地址 430074 湖北省武汉市洪山区珞喻路1037号