发明名称 用于风电并网系统容量效益裕度计算及分配的方法
摘要 一种用于风电并网系统容量效益裕度计算及分配的方法,属于电力系统分析、运行与控制技术领域。本发明的目的是对于风电并网系统,利用一种基于系统可靠性指标的概率性方法计算系统的容量效益裕度的用于风电并网系统容量效益裕度计算及分配的方法。本发明的步骤是:采用时间序列法对风机风速进行短时预测;根据风机风速与输出功率之间的关系曲线;各风机出力求和即得到风电场综合出力;绘制系统发电裕度的概率分布曲线;求出此时系统的可靠性指标LOLE的值;求出各区域每一时刻的容量效益裕度;对容量效益裕度进行经济分配计算。本发明对于多区域风电并网系统,以经济性最佳为目标,满足可靠性为约束对容量效益裕度进行经济分配。
申请公布号 CN104269886B 申请公布日期 2017.04.05
申请号 CN201410544509.2 申请日期 2014.10.16
申请人 东北电力大学;国家电网公司;国网吉林省电力有限公司松原供电公司 发明人 李国庆;王振浩;陈洪涛;刘玢;赵键明;项福军;佟辉
分类号 H02J3/46(2006.01)I;H02J3/00(2006.01)I;G06Q50/06(2012.01)I 主分类号 H02J3/46(2006.01)I
代理机构 吉林长春新纪元专利代理有限责任公司 22100 代理人 白冬冬
主权项 一种用于风电并网系统容量效益裕度计算及分配的方法,其特征在于:(1)采用时间序列法对风机风速进行短时预测,利用Jansen模型模拟风电场的尾流效应,X是两个风电机组的距离,叶轮半径为R,自然风速和受尾流影响的风速分别是v<sub>0</sub>和v<sub>x</sub>;平坦地形尾流风速计算式为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>v</mi><mi>x</mi></msub><mo>=</mo><msub><mi>v</mi><mn>0</mn></msub><mo>&CenterDot;</mo><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>C</mi><mi>T</mi></msub></mrow><mo>)</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><mfrac><mi>R</mi><mrow><mi>R</mi><mo>+</mo><mi>k</mi><mi>X</mi></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0001121545310000011.GIF" wi="1127" he="300" /></maths>C<sub>T</sub>取0.2,衰减常数k=0.5/ln(h/z<sub>0</sub>),其中h为轮毂高度,z<sub>0</sub>为粗糙度,取0.002;利用风机功率特性曲线得到风电场出力模型W;(2)根据风机风速与输出功率之间的关系曲线,当风速小于切入风速V<sub>ci</sub>或高于切出风速V<sub>co</sub>时,机组出力为0;当风速从切入风速到额定风速V<sub>r</sub>逐渐增大时,机组出力逐步增大;当达到或超过额定风速时,功率保持为额定功率P<sub>r</sub>不变;得到其数量关系可用下式表示:<img file="FDA0001121545310000012.GIF" wi="1326" he="395" />其中,<img file="FDA0001121545310000013.GIF" wi="1293" he="254" /><maths num="0002"><math><![CDATA[<mrow><mi>B</mi><mo>=</mo><mfrac><mn>1</mn><msup><mrow><mo>(</mo><msub><mi>V</mi><mrow><mi>c</mi><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>V</mi><mi>r</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mfrac><mo>&lsqb;</mo><mn>4</mn><mrow><mo>(</mo><msub><mi>V</mi><mrow><mi>c</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>V</mi><mi>r</mi></msub><mo>)</mo></mrow><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>V</mi><mrow><mi>c</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>V</mi><mi>r</mi></msub></mrow><mrow><mn>2</mn><msub><mi>V</mi><mi>r</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>3</mn></msup><mo>-</mo><mrow><mo>(</mo><mn>3</mn><msub><mi>V</mi><mrow><mi>c</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>V</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0001121545310000014.GIF" wi="1298" he="248" /></maths><maths num="0003"><math><![CDATA[<mrow><mi>C</mi><mo>=</mo><mfrac><mn>1</mn><msup><mrow><mo>(</mo><msub><mi>V</mi><mrow><mi>c</mi><mi>i</mi></mrow></msub><mo>-</mo><msub><mi>V</mi><mi>r</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mfrac><mrow><mo>&lsqb;</mo><mrow><mn>2</mn><mo>-</mo><mn>4</mn><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>V</mi><mrow><mi>c</mi><mi>i</mi></mrow></msub><mo>+</mo><msub><mi>V</mi><mi>r</mi></msub></mrow><mrow><mn>2</mn><msub><mi>V</mi><mi>r</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>3</mn></msup></mrow><mo>&rsqb;</mo></mrow></mrow>]]></math><img file="FDA0001121545310000015.GIF" wi="1426" he="359" /></maths>(3)根据以上分析得出风机出力,将各风机出力求和即得到风电场综合出力为:<maths num="0004"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi><mi>a</mi><mi>l</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>P</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>v</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001121545310000021.GIF" wi="1094" he="354" /></maths>(4)根据一天内每小时的风电场出力模型W,负荷预测模型L以及火电机组出力模型C,则系统发电裕度模型M=W+C‑L;利用卷积公式<maths num="0005"><math><![CDATA[<mrow><msub><mi>P</mi><mi>G</mi></msub><mrow><mo>(</mo><mi>W</mi><mo>+</mo><mi>C</mi><mo>=</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mo>-</mo><mi>&infin;</mi></mrow><mi>&infin;</mi></munderover><msub><mi>P</mi><mi>W</mi></msub><mrow><mo>(</mo><mi>W</mi><mo>=</mo><mi>z</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>P</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>C</mi><mo>=</mo><mi>k</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0001121545310000022.GIF" wi="1684" he="251" /></maths><img file="FDA0001121545310000023.GIF" wi="1590" he="230" />计算得到M的概率质量函数,绘制系统发电裕度的概率分布曲线;(5)根据第(4)步的概率分布,引入系统可靠性指标缺电时间期望LOLE,对M取负值区域的概率进行累积分布处理后,求出此时系统的可靠性指标缺电时间期望LOLE的值;根据系统的可靠性指标LOLE,求出缺电区域此时时刻从其互联系统输入的CBM的大小;绘制以CBM取值为横坐标,LOLE大小为纵坐标的曲线图,得到此时互联区域间需传输的容量效益裕度的大小;(6)对于多区域电网,根据步骤(1)、(2)、(3)、(4)、(5)求出各区域每一时刻的容量效益裕度,定义容量效益裕度为负的区域为受电区域,为正的区域为送电区域,以所花费用最小为目标,在送电区域各机组间将向受电区域输送的容量效益裕度进行分配,得到送电区域各机组分担的容量效益裕度,对容量效益裕度进行经济分配的目标函数为:<maths num="0006"><math><![CDATA[<mrow><mi>min</mi><mi> </mi><mi>C</mi><mrow><mo>(</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>&lsqb;</mo><msub><mi>F</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA0001121545310000024.GIF" wi="926" he="249" /></maths>约束条件为: 0≤q<sub>i</sub>≤P<sub>i</sub>‑p<sub>i</sub>;<maths num="0007"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>q</mi><mi>i</mi></msub><mo>&GreaterEqual;</mo><msub><mi>CBM</mi><mi>A</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0001121545310000031.GIF" wi="743" he="301" /></maths><maths num="0008"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>b</mi></munderover><msub><mi>q</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>CBM</mi><mi>B</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0001121545310000032.GIF" wi="726" he="309" /></maths><maths num="0009"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>c</mi></munderover><msub><mi>q</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>CBM</mi><mi>C</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0001121545310000033.GIF" wi="726" he="284" /></maths>其中:P<sub>i</sub>为机组i最大技术出力;p<sub>i</sub>为机组i现有的有功出力;CBM<sub>A</sub>为此时刻受电区域A为满足系统可靠性要求所需要的容量效益裕度;b、c分别为此时刻B、C区域可以提供容量效益裕度的机组数量;CBM<sub>B</sub>、CBM<sub>C</sub>为此时刻送电区域B、C在满足系统可靠性约束后可以提供的容量效益裕度;(7)对容量效益裕度进行经济分配的计算过程:设定种群规模size和遗传代数G,根据各机组装机容量和已发功率确定可用功率的范围,并产生初始种群,进行迭代,判断各区域机组是否满足容量效益裕度经济分配模型的约束条件: 0≤q<sub>i</sub>≤P<sub>i</sub>‑p<sub>i</sub><maths num="0010"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>q</mi><mi>i</mi></msub><mo>&GreaterEqual;</mo><msub><mi>CBM</mi><mi>A</mi></msub></mrow>]]></math><img file="FDA0001121545310000034.GIF" wi="481" he="200" /></maths><maths num="0011"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>b</mi></munderover><msub><mi>q</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>CBM</mi><mi>B</mi></msub></mrow>]]></math><img file="FDA0001121545310000035.GIF" wi="429" he="198" /></maths><maths num="0012"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>c</mi></munderover><msub><mi>q</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>CBM</mi><mi>C</mi></msub><mo>.</mo></mrow>]]></math><img file="FDA0001121545310000036.GIF" wi="468" he="181" /></maths>
地址 132012 吉林省吉林市长春路169号
您可能感兴趣的专利