发明名称 基于摄动法的风电场无功功率对电压的灵敏度计算方法
摘要 本发明涉及一种基于摄动法的风电场无功功率对电压的灵敏度计算方法,属于电力系统自动电压控制技术领域。本方法依据摄动法得到的风电场并网点处的无功功率对其电压的灵敏度辨识出风电场所接入的外部电网等值阻抗,从而和风电场内的电力网络共同组成完整的电力网络,并结合风电场内实时量测形成雅克比矩阵,再计算得到风电场内各节点间的无功功率对电压的灵敏度矩阵,最后再用摄动法得到的无功功率对电压的灵敏度对计算结果进行校核,若校核不通过,则对外部电网等值阻抗重新进行辨识,并重新计算无功功率对电压的灵敏度矩阵,直至灵敏度矩阵通过校核。依此方法,最终得到准确度较高的无功功率对电压的灵敏度矩阵,供风电场内电压控制所用。
申请公布号 CN102769300A 申请公布日期 2012.11.07
申请号 CN201210228752.4 申请日期 2012.07.02
申请人 清华大学;华北电网有限公司 发明人 郭庆来;孙宏斌;刘翊枫;许晓菲;江长明;贾琳;张伯明;吴文传;王彬;陈锐
分类号 H02J3/38(2006.01)I 主分类号 H02J3/38(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 罗文群
主权项 1.一种基于摄动法的风电场无功功率对电压的灵敏度计算方法,其特征在于该方法包括以下步骤:(1)在时间段T内,使风电场并入电网的并网点高压母线处的无功功率值变化量为:<img file="FDA00001842745400011.GIF" wi="130" he="55" />测定高压母线的电压变化量为<img file="FDA00001842745400012.GIF" wi="123" he="55" />则风电场并网点高压母线处的无功功率对该处电压的灵敏度摄动值<img file="FDA00001842745400013.GIF" wi="77" he="57" />为<img file="FDA00001842745400014.GIF" wi="272" he="132" />其中0秒&lt;T&lt;5秒;(2)设风电场并网点高压母线的电压基准值为<img file="FDA00001842745400015.GIF" wi="95" he="63" />根据上述灵敏度摄动值<img file="FDA00001842745400016.GIF" wi="104" he="57" />得到风电场外部电网等值阻抗Z<sub>eq</sub>为<img file="FDA00001842745400017.GIF" wi="332" he="63" />(3)根据风电场内部的电力网络模型和风电场外部的电网等值阻抗Z<sub>eq</sub>,构成风电场电力网络模型,并形成极坐标下的雅克比矩阵:<maths num="0001"><![CDATA[<math><mrow><mi>J</mi><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>V</mi><mi>P</mi></msub></mtd><mtd></mtd></mtr><mtr><mtd></mtd><mtd><msub><mi>V</mi><mi>Q</mi></msub></mtd></mtr></mtable></mfenced><mfenced open='(' close=')'><mtable><mtr><mtd><msup><mi>H</mi><mo>&prime;</mo></msup></mtd><mtd><msup><mi>N</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><msup><mi>M</mi><mo>&prime;</mo></msup></mtd><mtd><msup><mi>L</mi><mo>&prime;</mo></msup></mtd></mtr></mtable></mfenced><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>V</mi><mi>P</mi></msub></mtd><mtd></mtd></mtr><mtr><mtd></mtd><mtd><msub><mi>V</mi><mi>Q</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>上式中V<sub>P</sub>,V<sub>Q</sub>是风电场内电力网络节点的电压幅值的对角线矩阵,<maths num="0002"><![CDATA[<math><mrow><msup><mi>J</mi><mo>&prime;</mo></msup><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msup><mi>H</mi><mo>&prime;</mo></msup></mtd><mtd><msup><mi>N</mi><mo>&prime;</mo></msup></mtd></mtr><mtr><mtd><msup><mi>M</mi><mo>&prime;</mo></msup></mtd><mtd><msup><mi>L</mi><mo>&prime;</mo></msup></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mi>B</mi><mi>cos</mi><mi>&theta;</mi><mo>-</mo><mi>G</mi><mi>sin</mi><mi>&theta;</mi><mo>+</mo><mi>Q</mi></mtd><mtd><mo>-</mo><mi>G</mi><mi>cos</mi><mi>&theta;</mi><mo>-</mo><mi>B</mi><mi>sin</mi><mi>&theta;</mi><mo>-</mo><mi>P</mi></mtd></mtr><mtr><mtd><mi>G</mi><mi>cos</mi><mi>&theta;</mi><mo>+</mo><mi>B</mi><mi>sin</mi><mi>&theta;</mi><mo>-</mo><mi>P</mi></mtd><mtd><mi>B</mi><mi>cos</mi><mi>&theta;</mi><mo>-</mo><mi>G</mi><mi>sin</mi><mi>&theta;</mi><mo>-</mo><mi>Q</mi></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><img file="FDA000018427454000110.GIF" wi="355" he="56" /><img file="FDA000018427454000111.GIF" wi="371" he="56" />P<sub>i</sub>、Q<sub>i</sub>和V<sub>i</sub>分别为风电场内电力网络节点i的有功功率、无功功率和电压;B cosθ是B<sub>ij</sub> cosθ<sub>ij</sub>构成的矩阵,B sinθ是B<sub>ij</sub> sinθ<sub>ij</sub>构成的矩阵,G cosθ是G<sub>ij</sub> cosθ<sub>ij</sub>构成的矩阵,G sinθ是G<sub>ij</sub> sinθ<sub>ij</sub>构成的矩阵,G<sub>ij</sub>、B<sub>ij</sub>分别是风电场内电力网络节点导纳矩阵中导纳的实部和虚部,θ<sub>ij</sub>是节点i和节点j间的电压相角差;(4)根据上述极坐标下的雅克比矩阵,得到风电场内电力网络节点i的无功功率对风电场内电力网络节点j的电压的灵敏度<img file="FDA000018427454000112.GIF" wi="517" he="63" />依照上式计算风电场内电力网络所有节点间的无功功率对电压的灵敏度,进而得到风电场无功功率对电压的灵敏度矩阵S;(5)在时间段T内,使风电场内任意风机k的无功功率值变化量为<img file="FDA000018427454000113.GIF" wi="130" he="55" />保持风电场内其它设备的无功功率不变,测定风机k的机端电压值的变化量为<img file="FDA000018427454000114.GIF" wi="120" he="55" />则风机k的无功功率对风机k机端电压的灵敏度摄动值<img file="FDA000018427454000115.GIF" wi="74" he="61" />为:<img file="FDA000018427454000116.GIF" wi="373" he="120" />(6)从上述步骤(4)的风电场无功功率对电压的灵敏度矩阵S中,得到风机k的无功功率对风机k机端电压的灵敏度计算值<img file="FDA00001842745400021.GIF" wi="141" he="61" />设β为计算精度阈值,β的取值范围为0~20%,根据步骤(5)的灵敏度摄动值<img file="FDA00001842745400022.GIF" wi="99" he="61" />对灵敏度计算值<img file="FDA00001842745400023.GIF" wi="119" he="61" />进行以下判断:(6-1)若<img file="FDA00001842745400024.GIF" wi="384" he="152" />则回到步骤(1)重新开始计算;(6-2)若<img file="FDA00001842745400025.GIF" wi="384" he="151" />完成风电场无功功率对电压的灵敏度矩阵计算。
地址 100084 北京市海淀区清华园1号