发明名称 一种风电场非机理等效建模方法与装置
摘要 本发明提供一种风电场非机理等效建模方法和装置,以风电场的整体外特性为建模对象,分别以风电场并网点电压、风电场输入风速为输入量,均以风电场整体有功功率和无功功率的响应为输出量,通过参数辨识方法分别建立风电场响应与并网点电压变化之间传递函数,以及风电场响应与风速变化之间的传递函数。完整的风电场非机理等效模型由这两个传递函数的乘积组成。本发明提出的方法简化了风电场等效建模的过程和所需的数据量,同时又很好地契合了风电场等效建模的最终目的,即建立能够准确描述风电场外特性的简单低阶模型。通过本发明方法建立的风电场非机理等效模型可以广泛应用于含风电电力系统的仿真分析,具有良好的工程应用前景。
申请公布号 CN104167735A 申请公布日期 2014.11.26
申请号 CN201410438859.0 申请日期 2014.08.29
申请人 河海大学;国家电网公司;江苏省电力公司;江苏省电力公司电力经济技术研究院 发明人 吴峰;龚广京;鞠平;金宇清;孔卫亚;杨晓梅;黄俊辉;王海潜;谢珍建;乔黎伟;祁万春;谈健;赵宏大
分类号 H02J3/00(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 南京理工大学专利中心 32203 代理人 朱显国
主权项 一种风电场的非机理等效建模方法,其特征在于,使用风电场并网点的电压和测风塔的风速为输入量、使用风电场整体功率响应为输出量,通过辨识方法建立基于传递函数的风电场非机理等效模型,该非机理等效建模方法包括以下步骤:步骤A:建立风电场非机理等效模型的方程,具体步骤如下:步骤A‑1:对于风电场在电网故障下整体动态响应的情况,此时忽略风速的变化,建立风电场并网点电压和输出功率的非机理等效模型,其方程为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mi>P</mi><mo>=</mo><msub><mi>P</mi><mn>0</mn></msub><msub><mi>f</mi><mi>Pu</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>f</mi><mi>Pu</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>H</mi><mi>Pu</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>g</mi><mi>Pu</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>Q</mi><mo>=</mo><msub><mi>Q</mi><mn>0</mn></msub><msub><mi>f</mi><mi>Qu</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>f</mi><mi>Qu</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>H</mi><mi>Qu</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>g</mi><mi>Qu</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000562390550000011.GIF" wi="535" he="317" /></maths>式中,P、P<sub>0</sub>分别为风电场输出的有功功率和有功功率稳态值;Q、Q<sub>0</sub>分别为此时风电场输出的无功功率和无功功率稳态值;U为风电场并网点电压;H<sub>Pu</sub>(s)和H<sub>Qu</sub>(s)分别为传递函数方程;g<sub>Pu</sub>(U)和g<sub>Qu</sub>(U)分别为以风电场并网点电压U为自变量的代数方程;传递函数方程H<sub>Pu</sub>(s)和H<sub>Qu</sub>(s)分别表达如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>H</mi><mi>Pu</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>b</mi><mrow><mi>pu</mi><mn>0</mn></mrow></msub><msup><mi>s</mi><mi>pn</mi></msup><mo>+</mo><msub><mi>b</mi><mrow><mi>pu</mi><mn>1</mn></mrow></msub><msup><mi>s</mi><mrow><mi>pn</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>b</mi><mi>pun</mi></msub></mrow><mrow><msup><mi>s</mi><mi>pm</mi></msup><mo>+</mo><msub><mi>a</mi><mrow><mi>pu</mi><mn>1</mn></mrow></msub><msup><mi>s</mi><mrow><mi>pm</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>a</mi><mi>pum</mi></msub></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>H</mi><mi>Qu</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>b</mi><mrow><mi>qu</mi><mn>0</mn></mrow></msub><msup><mi>s</mi><mi>qn</mi></msup><mo>+</mo><msub><mi>b</mi><mrow><mi>qu</mi><mn>1</mn></mrow></msub><msup><mi>s</mi><mrow><mi>qn</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>b</mi><mi>qun</mi></msub></mrow><mrow><msup><mi>s</mi><mi>qm</mi></msup><mo>+</mo><msub><mi>a</mi><mrow><mi>qu</mi><mn>1</mn></mrow></msub><msup><mi>s</mi><mrow><mi>qm</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>a</mi><mi>qum</mi></msub></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000562390550000012.GIF" wi="810" he="355" /></maths>式中,pn≤pm,qn≤qm;b<sub>pu0</sub>、b<sub>pu1</sub>、…、b<sub>pun</sub>、a<sub>pu1</sub>、…、a<sub>pum</sub>、b<sub>qu0</sub>、b<sub>qu1</sub>、…、b<sub>qun</sub>、a<sub>qu1</sub>、…、a<sub>qum</sub>为常系数;代数方程g<sub>Pu</sub>(U)和g<sub>Qu</sub>(U)的表达式为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>g</mi><mi>Pu</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mi>U</mi><msub><mi>U</mi><mn>0</mn></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mtd></mtr><mtr><mtd><msub><mi>g</mi><mi>Qu</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mi>U</mi><msub><mi>U</mi><mn>0</mn></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000562390550000013.GIF" wi="328" he="278" /></maths>式中,U、U<sub>0</sub>分别为风电场并网点电压和并网点电压稳态值;步骤A‑2:对于风电场在风速波动下电功率变化的情况,建立风电场输入风速和输出功率的非机理等效模型,其方程为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mi>P</mi><mo>=</mo><msub><mi>P</mi><mn>0</mn></msub><msub><mi>f</mi><mi>Pv</mi></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>f</mi><mi>Pv</mi></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>H</mi><mi>Pv</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>g</mi><mi>Pv</mi></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>Q</mi><mo>=</mo><msub><mi>Q</mi><mn>0</mn></msub><msub><mi>f</mi><mi>Qv</mi></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>f</mi><mi>Qv</mi></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>H</mi><mi>Qv</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>g</mi><mi>Qv</mi></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000562390550000014.GIF" wi="520" he="317" /></maths>式中,P、P<sub>0</sub>分别为风电场输出的有功功率和有功功率稳态值;Q、Q<sub>0</sub>分别为此时风电场输出的无功功率和无功功率稳态值;V为风电场输入风速;H<sub>Pv</sub>(s)和H<sub>Qv</sub>(s)分别为传递函数方程;g<sub>Pv</sub>(V)和g<sub>Qv</sub>(V)分别为以风电场输入风速V为自变量的代数方程,该传递函数方程H<sub>Pv</sub>(s)和H<sub>Qv</sub>(s)表达式为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>H</mi><mi>Pv</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>b</mi><mrow><mi>pv</mi><mn>0</mn></mrow></msub><msup><mi>s</mi><mi>pn</mi></msup><mo>+</mo><msub><mi>b</mi><mrow><mi>pv</mi><mn>1</mn></mrow></msub><msup><mi>s</mi><mrow><mi>pn</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>b</mi><mi>pvn</mi></msub></mrow><mrow><msup><mi>s</mi><mi>pm</mi></msup><mo>+</mo><msub><mi>a</mi><mrow><mi>pv</mi><mn>1</mn></mrow></msub><msup><mi>s</mi><mrow><mi>pm</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>a</mi><mi>pvm</mi></msub></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>H</mi><mi>Qv</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>b</mi><mrow><mi>qv</mi><mn>0</mn></mrow></msub><msup><mi>s</mi><mi>qn</mi></msup><mo>+</mo><msub><mi>b</mi><mrow><mi>qv</mi><mn>1</mn></mrow></msub><msup><mi>s</mi><mrow><mi>qn</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>b</mi><mi>qvn</mi></msub></mrow><mrow><msup><mi>s</mi><mi>qm</mi></msup><mo>+</mo><msub><mi>a</mi><mrow><mi>qv</mi><mn>1</mn></mrow></msub><msup><mi>s</mi><mrow><mi>qm</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>a</mi><mi>qvm</mi></msub></mrow></mfrac></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000562390550000021.GIF" wi="809" he="354" /></maths>式中,pn≤pm,qn≤qm;b<sub>pv0</sub>、b<sub>pv1</sub>、…、b<sub>pvn</sub>、a<sub>pv1</sub>、…、a<sub>pvm</sub>、b<sub>qv0</sub>、b<sub>qv1</sub>、…、b<sub>qvn</sub>、a<sub>qv1</sub>、…、a<sub>qvm</sub>为常系数;代数方程g<sub>Pv</sub>(V)和g<sub>Qv</sub>(V)的表达式为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>g</mi><mi>Pv</mi></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mi>V</mi><msub><mi>V</mi><mn>0</mn></msub></mfrac><mo>)</mo></mrow><mn>3</mn></msup></mtd></mtr><mtr><mtd><msub><mi>g</mi><mi>Qv</mi></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mi>V</mi><msub><mi>V</mi><mn>0</mn></msub></mfrac><mo>)</mo></mrow><mn>3</mn></msup></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000562390550000022.GIF" wi="310" he="278" /></maths>式中,V为风电场输入风速;V<sub>0</sub>为风电场输入风速的稳态值;步骤A‑3:对于同时考虑风速变化和电网故障下风电场整体响应的情况,建立以风电场输入风速和风电场并网点电压为输入量,以风电场整体功率为输出量的非机理等效模型,其方程为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>P</mi><mo>=</mo><msub><mi>P</mi><mn>0</mn></msub><msub><mi>f</mi><mrow><mi>P</mi><mo>_</mo><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow><msub><mi>f</mi><mrow><mi>P</mi><mo>_</mo><mi>u</mi></mrow></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>f</mi><mrow><mi>P</mi><mo>_</mo><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>H</mi><mi>Pv</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>g</mi><mi>Pv</mi></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>f</mi><mrow><mi>P</mi><mo>_</mo><mi>u</mi></mrow></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>H</mi><mi>Pu</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>g</mi><mi>Pu</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>Q</mi><mo>=</mo><msub><mi>Q</mi><mn>0</mn></msub><msub><mi>f</mi><mrow><mi>Q</mi><mo>_</mo><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow><msub><mi>f</mi><mrow><mi>Q</mi><mo>_</mo><mi>u</mi></mrow></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>f</mi><mrow><mi>Q</mi><mo>_</mo><mi>v</mi></mrow></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>H</mi><mi>Pv</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>g</mi><mi>Pv</mi></msub><mrow><mo>(</mo><mi>V</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>f</mi><mrow><mi>Q</mi><mo>_</mo><mi>u</mi></mrow></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>H</mi><mi>Qu</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><msub><mi>g</mi><mi>Qu</mi></msub><mrow><mo>(</mo><mi>U</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mo>*</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000562390550000023.GIF" wi="1274" he="542" /></maths>式中,P、P<sub>0</sub>分别为风电场输出的有功功率和有功功率稳态值;Q、Q<sub>0</sub>分别为此时风电场输出的无功功率和无功功率稳态值;U、H<sub>Pu</sub>(s)、H<sub>Qu</sub>(s)、g<sub>Pu</sub>(U)、g<sub>Qu</sub>(U)的定义和表达式同步骤A‑1,V、H<sub>Pv</sub>(s)、H<sub>Qv</sub>(s)、g<sub>Pv</sub>(V)、g<sub>Qv</sub>(V)的定义和表达式同步骤A‑2;步骤B:辨识风电场非机理等效模型的参数,具体步骤如下:采用蚁群优化算法辨识前述步骤A‑1所建立等效模型H<sub>Pu</sub>(s)和H<sub>Qu</sub>(s)以及步骤A‑2所建立等效模型H<sub>Pv</sub>(s)和H<sub>Qv</sub>(s)的参数;步骤C、根据前述步骤B的参数辨识结果,代入前述步骤A‑3中(*)式得到风电场完整的非机理等效模型。
地址 210098 江苏省南京市鼓楼区西康路1号