董文凯 杜文娟 王海风,
用于风电场小干扰稳定性检验的降阶模式计算方法
董文凯1杜文娟2王海风1,2
(1. 新能源电力系统国家重点实验室(华北电力大学) 北京 102206 2. 四川大学电气工程学院 成都 610065)
考虑到在大型风电场规划阶段,如对所有方案采用全阶模型进行小干扰稳定性检验,会导致工作量大、效率低,不利于寻求最优方案,该文提出一种降阶模式计算方法。首先,建立风电场全阶线性化状态空间模型。然后,考虑风电场规划阶段,可对所选型号的风电机组(WTG)采用典型模型和参数,认为各WTG动态特性近似相同。通过对风电场线性化模型中的变量做等效变换,将台WTG构成的风电场等效为个相互独立的等效子系统。各等效子系统由一台并网WTG构成。根据等效子系统线性化状态空间模型的构建方式,提出一种用于风电场小干扰稳定性检验的降阶模式计算方法,可有效降低规划阶段风电场模式分析的计算量,且所得结果准确。最后,通过仿真算例,验证了所提方法的有效性。
风电场 线性化状态空间模型 小干扰稳定性检验 模式分析 等效子系统
随着风力发电技术的发展与广泛应用,大型风电场内可具有数百台风电机组(Wind Turbine Generator, WTG),WTG之间互联后产生的动态交互作用可能会对风电场小干扰稳定性产生显著影响[1-4]。因此,有必要在规划阶段,对风电场小干扰稳定性进行检验。而由于风电场规模的扩大,如果直接对所有WTG采用详细模型,所建立的风电场模型阶数高,求解计算量大、效率低[5-14]。因此,目前常采用风电场动态等效模型对风电场进行小干扰稳定性检验。建立风电场动态等效模型的常用方法可分为单机等效法和多机等效法两类。
单机等效法以一台等效WTG代表风电场内所有WTG,并将风电场内部的集电网络等效为一个外部阻抗[5-10,12-14]。集电网络等效常采用恒功率损耗法[12-14]。等效WTG容量为风电场额定容量,参数常采用容量加权平均法(Capacity Mean Value Method)或参数辨识法(Parameter Identification)确定。容量加权平均法以风电场内各WTG的额定容量与风电场额定容量的比值作为权重,计算所有WTG参数的加权平均值作为等效WTG的参数[12-14]。采用容量加权平均法建立风电场单机等效模型,计算量低;但是将等效模型用于风电场小干扰稳定性检验,所得结果的准确性仍有待证明[5,9]。
参数辨识法将等效WTG参数的计算转换为一个优化问题,目标函数为:相同工况或外部扰动下,等效WTG与风电场在输出功率等动态响应上的误差最小;待求量为等效WTG的参数[5-10]。由于参数辨识法涉及非线性优化问题的求解[5-10],一般需采用智能优化算法,如粒子群算法[6,8]、遗传算法[5,7],以及基于它们提出的改进算法[10]等。因此,参数辨识法在提高风电场单机等效模型准确度的同时,也增加了等效模型参数求解的计算量。而在风电场规划阶段,典型工况或外部扰动下,风电场输出功率等动态响应难以实测获取,需进行非线性仿真。在多种场景下对风电场动态响应进行仿真分析,并将仿真结果用于参数辨识,可在一定程度上提高等效模型的准确度,但也会导致工作量的增加。
多机等效法首先根据运行状态对WTG进行聚类,将风电场划分为几个由运行状态相似的WTG构成的风电机群[11-14],聚类分析常采用k均值聚类、c均值聚类和支持向量机等方法,聚类指标一般选用各WTG转子转速、端电压或输入风速等能够反映WTG运行状态的数据。然后,再对各风电机群作单机等效,等效方法常采用容量加权平均法[11-14],也有少量研究中采用参数辨识法[5]。由于考虑了各WTG运行状态的差异,多机等效法提高了等效准确度;但是对WTG进行聚类时,需要WTG的运行数据作为聚类指标。而在风电场规划阶段,若要获取这些数据,需进行非线性仿真。
综上所述,容量加权平均法用于风电场小干扰稳定性检验,所得结果的准确性仍有待进一步探讨。参数辨识法或多机等效法则需要风电场或WTG的实测或仿真数据的支持,才可建立风电场动态等效模型,用于风电场规划阶段进行小干扰稳定性检验也存在一定的局限性。
本文针对风电场规划中的小干扰稳定性检验,提出了一种降阶模式计算方法。首先对风电场内所有WTG采用详细模型,建立风电场全阶线性化模型。然后,考虑在风电场规划阶段,可对所选型号的WTG采用典型模型和参数;进而通过对风电场线性化模型中的变量作等效变换,将台WTG构成的风电场等效为个由一台WTG构成的等效子系统,并提出一种风电场降阶模式计算(Reduced- Order Modal Computation, ROMC)方法。最后,通过仿真算例,验证了所提方法的有效性。
图1为一台WTG构成的风电场的示意图,L为风电外送通道的线路阻抗。此外,图1也可看作风电场内的一个风电机群,此时L为风电机群与风电场其余部分之间连接线路的阻抗。本文主要以风电场为例展开研究,所提方法也可用于将风电场内某一个风电机群看作一独立系统时,对该风电机群的小干扰稳定性检验。链式和干线式结构是风电场最常用的连接结构,本文将在这两种结构下对风电场小干扰稳定性问题展开研究。
图1 并网风电场结构
风电场内第台WTG的线性化状态空间模型可表示为
其中
在风电场规划阶段,考虑外部系统容量远大于风电场容量,可将母线C看作无穷大母线[15],对风电场作为一独立系统的小干扰稳定性进行检验。由附录中的推导,可得风电场网络方程为
由式(1)可知
由式(2)和式(3)可得风电场线性化模型为
其中
在风电场规划阶段,可对所选型号的WTG采用典型模型和参数,而且风电场内风速空间分布一般不会有太大差异,因此各WTG线性化模型相似,可假设各WTG的线性化模型相同,从而有
文献[16]中采用了一种对风电场动态系统进行等效分解的思路,将台WTG构成的风电场分解为个相互独立且由一台WTG构成的等效子系统,进而研究了风电场动态等效方法及其理论依据;该方法对于研究如何降低风电场小干扰稳定性分析的计算量具有很好的参考意义。但是文献[16]中的推导是针对所有WTG经阻抗相同的线路并联接入汇流母线这一特殊结构下的风电场展开的,其中所提等效分解方法也仅适用于上述特定场景,无法用于对本文图1所示的这样一类更具一般性的风电场的等效分解。本文借鉴这一等效分解的思想,并结合模式分析理论中对矩阵进行模式分解的思路[17],推导了WTG非对称连接所构成风电场线性化模型的等效分解。具体推导过程如下。
2.1.1 风电场网络电抗矩阵的构建
由附录式(A3),忽略线路电阻,可得
由式(6),可得到风电场网络电抗矩阵为
2.1.2 风电场线性化模型等效分解
令
为对风电场线性化模型进行等效分解,首先对式(3)中的变量做变换
然后,将式(10)与式(5)代入式(3)得
由式(8)与式(9)可得
最后,对于风电场网络方程,由式(8)的第2个方程
由式(13)和式(15),经式(10)所示变量变换,原台WTG构成的风电场的线性化模型可等效为如下个相互独立的子系统,即
其中,第个等效子系统的线性化状态空间模型为
在2.1小节将式(4)所示风电场全阶模型等效为式(17)所示的个等效子系统的推导过程中,作出了两点假设。假设1:风电场内各WTG的线性化模型相同;假设2:风电场连接线路电阻为0。对于上述两点假设作进一步分析与解释如下。
对于假设1:①一个风电场内一般采用同一型号的WTG,其初始参数设定相同。②风速差异会导致各WTG稳态功率输出具有一定差异,即各WTG的稳态运行点会有一定差异;但是在WTG型号相同的条件下,这一差异对WTG线性化模型影响不大,可认为各WTG的线性化模型“近似相同”。③一些规模较大的风电场可能由几个风电机群构成,各风电机群内WTG型号相同且风速空间分布差异不大,对此可依次将各个风电机群看作一独立系统,检验其小干扰稳定性;此时便是将图1所示系统看作风电场内的一个风电机群的情况;结合第1节的说明,风电场其余部分与外部交流系统总容量远大于该风电机群的容量,因此在规划阶段可认为母线C电压恒定,式(2)和式(4)分别为该风电机群的网络方程和全阶线性化模型。
因此,可认为:①风电场内各WTG线性化模型近似相同;②连接线路电阻的影响可忽略。
采用ROMC法检验风电场小干扰稳定性的步骤总结如下:
(1)ROMC法所得风电场模式计算结果与由全阶模型所得结果近似,因此,若风电场小干扰稳定性良好,则由ROMC法所得风电场主导模式应位于复平面左半平面且具有足够的稳定裕度;若风电场存在小干扰失稳风险,则由ROMC法所得风电场主导模式靠近虚轴或位于复平面右半平面。
(2)在风电场规划阶段,通常需要考虑多种设计方案,并从经济性、效益和稳定性等角度,对其进行综合对比和评价,以获得最优设计方案。为确保风电场实际运行中的小干扰稳定性,规划阶段应保证风电场主导模式均位于复平面左半平面,且具有一定的稳定裕度。对于各设计方案,如果均采用全阶模型作模式分析或非线性仿真分析,虽然结果准确,但也会造成小干扰稳定性检验工作量大、效率低。
(3)考虑上述分析与问题,可在风电场规划阶段引入ROMC法,一方面,可首先淘汰一部分明显具有小干扰稳定性问题或失稳风险的方案,不必对此类方案再作非线性仿真分析;另一方面,对于可初步确保小干扰稳定性的方案,也可根据ROMC法所得模式计算结果对其小干扰稳定性的优劣给出评价,并用于最优方案的选取。对于最优方案,可采用全阶模型和非线性仿真分析对其作进一步检验,以有效确保所选方案的小干扰稳定性。
综上所述,ROMC法用于风电场规划,可有效降低小干扰稳定性检验的工作量,也有助于尽可能地考虑更多方案,避免遗漏。对于大规模风电场,WTG数量大,可供考量的方案多,且风电场全阶模型阶数高,此时采用ROMC协助完成风电场小干扰稳定性检验工作,可有效提高工作效率,也有助于对比多种方案,获得最优解。
本节通过仿真算例验证ROMC法用于风电场小干扰稳定性检验的有效性。风电场采用DFIG风电场,DFIG模型见文献[22],参数见附录中附表1。仿真软件采用Matlab。
图2 算例风电场结构图(三支路)
图3 风电场模式计算结果(变化时)
Tab.1 Part of the computational results when network configuration varied when varies
由图3及附表3和表1可得:
(1)采用ROMC法所得风电场振荡模式与由全阶模型所得结果在复平面上分布情况基本一致。验证了ROMC法用于风电场小干扰稳定性检验的有效性。
图4 非线性仿真结果(变化时)
图5 其他结构下的算例风电场
对于风电场采用三支路和链式结构的情况,采用ROMC法计算风电场振荡模式,过程同上述对并联情况的处理,得到振荡模式1在复平面上的分布分别如图6b和图6c中的“×”所示。
图6 风电场模式计算结果(不同结构)
表2 部分振荡模式计算结果(不同结构)
Tab.2 Part of the computational results when network configuration varied(different configurations)
由图6和表2可得如下分析结论:
(1)采用ROMC法所得风电场振荡模式与由全阶模型所得结果在复平面上的分布情况基本一致,验证了ROMC法的有效性。
(2)当风电场采用链式结构时,有一振荡模式位于复平面右半平面,ROMC法有效检验了风电场小干扰失稳风险。非线性仿真验证结果如图7所示。仿真中的扰动设置同3.1小节。
(3)在风电场规划阶段,应注意集电网络结构变化对风电场小干扰稳定性的影响,且应尽量避免采用链式结构将多台WTG接入汇流母线的连接方式。
图7 非线性仿真结果(不同结构)
采用ROMC法检验风电场小干扰稳定性。
图8 大型算例风电场结构
图9 模式计算结果(大型风电场)
对于风电场采用图8b所示结构的情况,采用ROMC法计算风电场振荡模式,过程同上述对采用图8a所示结构的情况的处理,得到振荡模式1在复平面上的分布如图9b中“×”所示。
表3 部分振荡模式计算结果(大型风电场)
Tab.3 Part of the computational results when network configuration varied(large-scale wind farm)
由图9可得如下结论:
(1)采用ROMC法所得风电场振荡模式与由全阶模型所得结果在复平面上的分布情况基本一致,验证了ROMC法的有效性。
(2)当风电场采用图8b所示结构时,有一振荡模式位于复平面右半平面。ROMC法有效检验了风电场小干扰失稳风险。非线性仿真验证结果如图10所示。仿真中的扰动设置同3.1小节。
(3)算例结果再次表明,集电网络结构变化会影响风电场小干扰稳定性,应尽量避免采用链式结构将多台WTG接入汇流母线的连接方式。
图10 非线性仿真结果(大型风电场)
本文提出了一种用于风电场小干扰稳定性检验的ROMC法。
2)根据等效子系统线性化状态空间模型的构建方法,提出ROMC法。在规划阶段,用ROMC法检验风电场小干扰稳定性,所得结果准确,模式分析的计算量显著降低。
4)文中等效分解的推导是在风电场内各WTG线性化模型相似的条件下展开的,等效分解方法适用于由线性化模型相似的发电单元所构成的系统,如同型风电场或风电机群以及同型光伏电站,在规划阶段可认为其中设备的线性化模型相似,进而采用该方法对系统进行等效分解;对于发电单元线性化模型差异较大的系统,如多机电力系统,其中发电机容量和参数等差异较大,目前无法采用该方法分析系统动态特性。对于如何将等效分解的方法扩展至多机电力系统,则还有待进一步研究。
不考虑母线C电压波动,式(A1)中
其中
由式(A1)和式(A2)可得风电场网络方程为
式中
附表1 算例DFIG参数
App.Tab.1 Parameters of DFIG in study cases
换流器模式控制器 RSC外环有功0.14 无功0.12 内环d轴电流0.0135 q轴电流0.0110 GSC外环直流电压1.515 PLL3.189.42 内环d轴电流2200 q轴电流2200
附表2 图2风电场内部线路参数
App.Tab.2 Parameters of lines inside the wind farm of Fig.2
iii 10.004+j0.01770.004+j0.018 7130.004+j0.017 20.004+j0.015 380.004+j0.015 3140.004+j0.017 30.004+j0.018 790.003+j0.013 6150.003+j0.013 6 40.005+j0.021 3100.003+j0.012 8160.004+j0.018 7 50.004+j0.017110.004+j0.019 6170.004+j0.018 7 60.004+j0.017120.004+j0.017
附表3 3.1小节风电场网络电抗矩阵特征值计算结果
App.Tab.3 Eigenvalues of network reactance matrixes in 3.1
0.040.003 8, 0.004, 0.004 3, 0.005 3, 0.005 5, 0.006 5, 0.00 72, 0.009 9, 0.012, 0.015, 0.017, 0.028, 0.035, 0.045, 0.212, 0.254, 0.90 0.10.003 8, 0.004, 0.004 3, 0.005 3, 0.005 5, 0.006 5, 0.007 2, 0.009 9, 0.012, 0.015, 0.017, 0.028, 0.035, 0.047, 0.0212, 0.254, 1.92 0.150.003 8, 0.004, 0.004 3, 0.005 3, 0.005 5, 0.006 5, 0.007 2, 0.009 9, 0.012, 0.015, 0.017, 0.028, 0.035, 0.047, 0.212, 0.255, 2.77 0.20.003 8, 0.004, 0.004 3, 0.005 3, 0.005 5, 0.006 5, 0.007 2, 0.009 9, 0.012, 0.015, 0.017, 0.028, 0.035, 0.048, 0.212, 0.255, 3.62
附表4 3.2小节风电场网络电抗矩阵特征值计算结果
App.Tab.4 Eigenvalues of network reactance matrixes in 3.2
并联0.012 9, 0.013 6, 0.014 1, 0.015 3, 0.015 7, 0.017, 0.017, 0.017, 0.017, 0.017, 0.018, 0.019, 0.019, 0.019, 0.019 4, 0.021, 2.567 三支路同附表3 计算结果 串联0.003 8, 0.004 2, 0.004 6, 0.005, 0.005 3, 0.005 8, 0.006 7, 0.007 7, 0.01, 0.011 4, 0.015 3, 0.021, 0.033 2, 0.053 6, 0.117 6, 0.365, 4.49
[1]王一凡, 赵成勇, 郭春义. 双馈风电场孤岛经模块化多电平换流器直流输电并网系统小信号稳定性分析与振荡抑制方法[J]. 电工技术学报, 2019, 34(10): 2116-2129. Wang Yifan, Zhao Chengyong, Guo Chunyi. Small- signal stability and oscillation suppression method for islanded doubly fed induction generator-based wind farm integrated by modular multilevel converter based HVDC system[J]. Transactions of China Electrotechnical Society, 2019, 34(10): 2116-2129.
[2]苏田宇, 杜文娟, 王海风. 多直驱永磁同步发电机并联风电场次同步阻尼控制器降阶设计方法[J]. 电工技术学报, 2019, 34(1): 116-127. Su Tianyu, Du Wenjuan, Wang Haifeng. A reduced order design method for sub-synchronous damping controller of multi-PMSGs parallel wind farm[J]. Transactions of China Electrotechnical Society, 2019, 34(11): 116-127.
[3]周佩朋, 宋瑞华, 李光范, 等. 直驱风电机组次同步振荡阻尼控制方法及其适应性[J]. 电力系统自动化, 2019, 43(13): 177-188. Zhou Peipeng, Song Ruihua, Li Guangfan, et al. Damping control of subsynchronous oscillation for direct drive wind turbine generator and its adaptability[J]. Automation of Electric Power System, 2019, 43(13): 177-188.
[4]徐路遥, 兰洲, 陈飞, 等. 弱电网下无功控制对并网变流器稳定性影响分析[J]. 电力系统自动化, 2019, 43(15): 56-68. Xu Luyao, Lan Zhou, Chen Fei, et al. Analysis on impact of reactive power control on stability of grid-connected converter in weak grid[J]. Automation of Electric Power System, 2019, 43(15): 56-68.
[5]Wang Yinfeng, Lu Chao, Zhu Lipeng, et al. Comprehensive modeling and parameter identification of wind farms based on wide-area measurement systems[J]. Journal of Modern Power Systems and Clean Energy, 2016, 4(3): 383-393.
[6]El-Sharkawi M. Dynamic equivalent models for wind power plants[C]//2011 IEEE Power and Energy Society General Meeting, Detroit, 2011, DOI: 10.1109/PES.2011.6038934.
[7]Li H, Yang Chao, Zhao B, et al. Aggregated models and transient performances of a mixed wind farm with different wind turbine generator systems[J]. Electric Power Systems Research, 2012, 92: 1-10.
[8]张剑, 何怡刚. 基于轨迹灵敏度分析的永磁直驱风电场等值模型参数辨识[J]. 电工技术学报, 2020, 35(15): 3303-3313. Zhang Jian, He Yigang. Parameters identification of equivalent model of permanent magnet synchronous generator (PMSG) wind farm based on analysis of trajectory sensitivity[J]. Transactions of China Electrotechnical Society, 2020, 35(15): 3303-3313.
[9]Zhou Yuhao, Zhao Long, Lee W. Robustness analysis of dynamic equivalent model of DFIG wind farm for stability study[J]. IEEE Transactions on Industry Applications, 2018, 54(6): 5682-5690.
[10]Erlich I, Shewarega F, Feltes C, et al. Determination of dynamic wind farm equivalents using heuristic optimization [C]// Power & Energy Society General Meeting, San Diego, CA, USA, 2012, DOI: 10.1109/ PESGM.2012.6345508.
[11]Zou Jianxiao, Peng Chao, Xu Hongbin, et al. A fuzzy clustering algorithm-based dynamic equivalent modelingmethod for wind farm with DFIG[J]. IEEE Transactions on Energy Conversion, 2015, 30(4): 1329-1337.
[12]吴红斌, 何叶, 赵波, 等. 基于改进K-means聚类算法的风电场动态等值[J]. 太阳能学报, 2018, 39(11): 3232-3238. Wu Hongbin, He Ye, Zhao Bo, et al. Research on dynamic equivalent of wind farm based on improved K-means clustering algorithm[J]. Acta Energiae Solaris Sinica, 2018, 39(11): 3232-3238.
[13]蔡文超, 杨炳元. 基于RTDS风电场等值建模若干问题的解决方案[J]. 电气技术, 2017, 18(6): 13-18. Cai Wenchao, Yang Bingyuan. Solution of some problems inequivalent modeling of wind farm based on RTDS[J]. Electric Engineering, 2017, 18(6): 13-18.
[14]Ali M, Ilie I, Milanovic J, et al. Wind farm model aggregation using probabilistic clustering[J]. IEEE Transactions on Power Systems, 2013, 28(1): 309-316.
[15]Rogers G. Power System Oscillations[M]. New York: Springer Science + Business Media, 2000.
[16]Du Wenjuan, Dong Wenkai, Wang Haifeng, et al. Dynamic aggregation of same wind turbines generators in parallel connection for studying oscillation stability of a wind farm[J]. IEEE Tran- sactions on Power Systems, 2019, 34(6): 4694-4705.
[17]杜文娟, 王海风. 电力系统低频功率振荡模式分析理论与方法[M]. 北京: 科学出版社, 2017.
[18]邵冰冰, 赵书强, 高本锋, 等. 连接弱交流电网的VSC-HVDC失稳机理及判据研究[J]. 电工技术学报, 2019, 34(18): 3884-3896. Shao Bingbing, Zhao Shuqiang, Gao Benfeng, et al. Instability mechanism and criterion analysis of VSC- HVDC connected to the weak AC power grid[J]. Transactions of China Electrotechnical Society, 2019, 34(18): 3884-3896.
[19]杨苓, 陈艳东, 罗安, 等. 多机并联接入弱电网的改进型带阻滤波器高频振荡的抑制[J]. 电工技术学报, 2019, 34(10): 2079-2091. Yang Ling, Chen Yandong, Luo An, et al. Suppression method of high-frequency oscillation by improved notch filter for multi-parallel inverters connected to weak grid[J]. Transactions of China Electrotechnical Society, 2019, 34(10): 2079-2091.
[20]段贵钟, 秦文萍, 逯瑞鹏, 等. 考虑风电接入与负荷不确定性的静态电压稳定性分析[J]. 电力系统保护与控制, 2018, 46(12): 108-114. Duan Guizhong, Qin Wenping, Lu Ruipeng, et al. Static voltage stability analysis considering the wind power and uncertainty of load[J]. Power System Protection and Control, 2018, 46(12): 108-114.
[21]徐政. 交直流电力系统动态行为分析[M]. 北京: 机械工业出版社, 2004.
[22]Fan Lingling, Kavasseri R, Miao Zhixin. Modeling of DFIG-based wind farms for SSR analysis[J]. IEEE Transactions on Power Delivery, 2010, 25(4): 2073-2082.
Reduced-Order Modal Computation Method for Small-Signal Stability Examination of a Wind Farm
Dong Wenkai1Du Wenjuan2Wang Haifeng1,2
(1. State Key Laboratory of Alternate Power Systems with New Energy Resources North China Electric Power University Beijing 102206 China 2. School of Electric Engineering Sichuan University Chengdu 610065 China)
In the planning stage of a large-scale wind farm, if the small-signal stability of all the schemes is examined by using full-order model, the large workload and low working efficiency will affect the solving of the optimal scheme, thus in this paper, a reduced-order modal computation method was proposed. The full-order linearized state-space model was established at first. Secondly, it was considered that in the planning stage of a wind farm, typical model and parameters of the type of wind turbine generator (WTG) selected could be used and the dynamics of the WTGs were similar. Thus a wind farm comprised ofWTGs could be equivalently decoupled intoindependent subsystems through equivalent transformation on the variables of the linearized model. Each subsystem was comprised of a single grid-connected WTG. According to the state-space model of the equivalent subsystems, a reduced-order modal computation (ROMC) method was proposed, which could be used to exam the small-signal stability of a wind farm in planning stage accurately with computational burden obviously reduced. Finally, the effectiveness of the method proposed was validated through study cases.
Wind farm, linearized state-space model, small-signal stability examination, modal analysis, equivalent subsystem
TM712
10.19595/j.cnki.1000-6753.tces.200158
国家电网公司科技项目资助(5100-202036013A-0-0-00)。
2020-02-19
2020-04-09
董文凯 男,1992年生,博士研究生,研究方向为风电场动态等值与小干扰稳定性。E-mail:596167281@qq.com
杜文娟 女,1979年生,博士,教授,研究方向为电力系统稳定性分析与控制。E-mail:ddwenjuan@qq.com(通信作者)
(编辑 郭丽军)