-
随着化石能源等传统能源的日益减少，风力发电作为一种新能源发电技术迅速崛起。然而，随着风电装机容量在电网总装机容量的比重不断提高，风力发电对电网的稳定性带来严重考验。不同类型的风电机组对电网暂态稳定性的影响不同，为了深入地了解各风力发电机组特性以及不同风机对电力系统暂态稳定性影响程度，对风力发电系统进行仿真分析是十分必要的。风电系统的暂态稳定性研究这一课题已被很多科研工作者和电网工作人员青睐，他们也做出了一定的成果。但是对这一课题的研究大都是针对单一类型风机的发电系统，本文将双馈发电系统、直驱发电系统、恒速恒频发电系统对比进行研究，得出了相关结论。另外本文所用的电力系统仿真软件PSAT是近年兴起的电力系统仿真软件，具有小信号稳定分析和时域仿真、潮流计算、连续潮流计算、最优潮流求解等模块，为用户进行电力系统仿真提供了极大的便利。
-
PSAT是一种新的电力系统分析软件，其命令行版本也是GNU Octave兼容的［1-3］。它的核心是潮流计算，会对状态变量进行初始化处理。一旦求出了系统潮流，就能执行后续的静态或动态分析。PSAT包含计算模块有：小信号稳定分析和时域仿真、潮流计算、连续潮流计算、最优潮流求解等等。
PSAT具有良好的适用于电力系统仿真的界面，和普通交互界面一样，由菜单栏、状态栏、控件栏等组件。在进行电力系统计算时，在Data file 文本框输入原始计算数据，在Perturbation file 栏录入动态更细数据，在Command line写入命令以便执行。计算边界条件实时显示在主界面右端，包括计算基准容量、频率、仿真开始及结束时间，收敛门槛等。实时计算结果包括潮流计算结果、最优潮流结果，收敛结果等显示在主界面的右边。
-
双馈风力发电系统结构图如图1所示。
图 1 双馈风力发电系统结构图
Figure 1. Double-fed wind power system
为了更好分析问题得出结论，在分析双馈风力发电系统前须作出以下假设：
（1）三相绕组在空间上以120°角对称分布，产生的磁动势以正弦波形式沿圆周分布；（2）各定子绕组转子绕组自感和互感线性对称，不考虑磁饱和；（3）转子绕组阻抗已换算到定子绕组侧，换算完成后，每相匝数相同；（4）不考虑铁芯损耗，绕组电阻不受温度和频率变化的影响。
作出以上假设后，就可以在同步旋转坐标系（d-q）下建立双馈风力发电机的数学模型，双馈风机的数学模型在很多资料中都有介绍，这里引用［6］中的模型，其中定子侧采用发电机惯例，转子侧采用电动机惯例。
电压方程如下：
(（1）) 磁链方程如下：
(（2）) 式中：rs、rr为定、转子每相电阻（Ω）；Ls、Lr、Lm为定、转子等效电感和互感（H）；ids、iqs、idr、iqr分别为定子d-q轴和转子d-q轴电流（A）；
换流器的存在使得双馈风力机组的机械部分与电气部分被解耦，本论文的换流器模型没有分为定子侧换流器和转子侧换流器来分别介绍，只是把换流器看成一个整体，与发电机电磁暂态过程相比，换流器的动态过程快很多，因此，换流器可以用一个理想电流原来代替。如文献［7］介绍。将转子电流按d-q轴分解为iqr和idr，分别用于控制转速模块与电压模块，用方程描述如下：
(（3）) (（4）) 式中：P*ω（ωm）具有功率—转速特性，当转速在某一范围内变化时，它也变化，以便获得最大的功率。同时发电机发出的无功和它的电磁转矩也由转子直轴电流idr和转子交流电流iqr控制，如下所示：
(（5）) (（6）) 由于电磁转矩和有功密切相关，由以上两式可以看出流过换流器的电流idr和iqr与发电机发出的无功和有功正相关。因为换流对流过它的电流有较高的要求，因此可以设定发电机发出的有功和无功的极值来限制idr和iqr的大小，方程如（7），即设置有功无功的最大值和最小值来限制流过换流器的电流。
(（7）) -
直驱同步发电系统见图2。
图 2 直驱风力发电系统结构图
Figure 2. Direct drive wind power system
在建立同步发电系统方程时，为了计算方便可以假定电机处于电气稳态，因为定子和转子的磁链变化比电网的变化快，因此同步机组的机械部分和电气部分完全解耦。则可以得出如下方程：
电压方程如下：
(（8）) 磁动势方程如下：
(（9）) 将换流器看作是理想电源，定子交轴电流iqs，定子直轴电流ids和定子直流电流分量idc是状态变量，将他们分别用于转子速度控制，功率控制和电压控制，由于电机速度和电压都有一个限值，因此判断机组是否失稳的标准是看控制电流是否越限。控制电流的限值方程如下所示：
(（10）) 式中：
-
恒速风力发电系统如图3所示：
图 3 恒速恒频风力发电系统结构图
Figure 3. Constant speed constant frequency wind powersystem
恒速恒频发电系统中元件的模型相对较简单，可以直接在直角坐标系中建立，以网络平衡节点为参考点。可得方程如下：
电压方程如下：
(（11）) 功率方程如下：
(（12）) 电磁转矩方程如下：
(（13）) 式中：vr，vm为定子电压的实轴和虚轴分量（V）；ir，im为电流分量（A）；bc为混合电容对应的电导（S）；
-
本论文采用的仿真模型为WSCC的三机九节点模型，见图4。
图 4 三机九节点系统图
Figure 4. 3-genreator 9-node power system
图4是以双馈风机为例，在对其它类型的风电系统进行仿真时，只需将BUS2上的双馈风机改成其它类型的风机即可。BUS8处设置的为三相短路故障。论文研究在三相短路故障和甩负荷两种情况下，各风电系统的暂态稳定情况。关于暂态稳定分析的数学原理，即功角方程解法在相关文献［8］中有介绍，本文只给出系统仿真图。
-
在BUS8处，将短路故障的开始时间设置为1秒，故障持续时间分别为120 ms、246 ms和625 ms。可得出如仿真图5。
图 5 故障持续120 ms时直驱风电系统同步电机的σ-t曲线
Figure 5. The σ-t figure of direct drive wind power system in 120 ms fault duration
图 6 故障持续120 ms时恒速机组的转速曲线
Figure 6. The rotate speed figure of constant speed wind power system in 120 ms fault duration
从图6可以看出，在故障开始后，恒速风机的转速急剧变化，但在故障清除后，转速开始回稳，大约在4 s时趋于稳定，即恒速风机本身是稳定的。
当故障持续时间为246 ms时，含双馈风机的电力系统能保持暂态稳定，但转子电流越限，如图7所示，因此机组本身不能保持稳定；含直驱同步电机的电力系统的功角变化是同步的，转子电流不越限，因此它们也都保持稳定；含恒速机组的电力系统保持稳定，但功角比故障前的功角大。
图 7 故障持续246 ms时双馈风机的转子电流曲线
Figure 7. The rotor current figure of double-fed wind power system in 246 ms fault duration
当故障持续时间为625 ms时，含双馈风机的电力系统仍然能保持稳定，但转速和功角都比原来大很多，已接近不稳定状态，转子电流继续越限，双馈机组本身不能保持稳定；含直驱风机的电力系统不稳定，机组也不稳定，仿真不能完全结束；含恒速风机的风电系统不稳定，恒速机组本身也不稳定，仿真结果不收敛。
通过以上分析，将三种发电系统以及机组本身的暂态稳定性结果填入表1。
表 1 三种风电机组及系统稳定性对比
Table 1.
Contrast of the stability between wind plants 类型 故障持续时/ms 100 246 625 双馈发电系统 稳定 稳定 不稳定 双馈机组本身 稳定 不稳定 不稳定 直驱发电系统 稳定 稳定 不稳定 直驱机组本身 稳定 稳定 不稳定 恒速发电系统 稳定 稳定 不稳定 恒速机组本身 稳定 不稳定 不稳定 同时，为了研究三种风电机组在故障持续时间为625 ms时的低电压穿越能力，得到BUS7处电压曲线如图8~图10所示。
图 8 故障持续625 ms时恒速风电系统BUS7处电压
Figure 8. The voltage of BUS7 of constant speed wind power system in 625 ms fault duration
图 9 故障持续625 ms时直驱风电系统BUS7处电压
Figure 9. The voltage of BUS7 of direct drive wind power system in 625 ms fault duration
图 10 故障持续625 ms时双馈风电系统BUS7处电压
Figure 10. The voltage of BUS7 of double-fed wind power system in 625 ms fault duration
以上图形表明，当故障持续时间为246 ms时可明显看出直驱风机的低电压穿越能力较好，当故障时间达到625 ms时，尽管三种风电机组都有一定的低电压穿越能力，但系统最终仍然失去稳定性。
-
当短路故障持续时间较长后，尽管低电压穿越能力对系统有一定积极作用，但不能使系统恢复稳定。此时，考虑在系统中加入无功补偿，容量取为风机容量的40%，为40 MVar。本算例在故障点BUS8处集中加入静止无功补偿器，仿真图形如图11~图13所示。
图 11 故障持续625 ms加入SVC后恒速风电系统BUS7处电压
Figure 11. The voltage of BUS7 of constant speed wind power system in 625 ms fault duration with SVC
图 12 故障持续625 ms加入SVC后直驱风电系统BUS7处电压
Figure 12. The voltage of BUS7 of direct drive wind power system in 625 ms fault duration with SVC
图 13 故障持续625 ms加入SVC后双馈风电系统BUS7处电压
Figure 13. The voltage of BUS7 of double-fed wind power system in 625 ms fault duration with SVC
由以上图形可看出，SVC加入对恒速风电系统的作用仍然不明显，但对双馈风电系统和直驱风电系统有较明显作用，使得BUS7处的电压能较好的恢复。但由后续分析可知，当故障持续时间继续增加时，尽管有无功补偿容量，但上述三种风电系统仍然不能保持稳定。
-
当无功负荷分别由0.45、0.5和0.45个标幺值减少到0.05、0.05和0.05个标幺值，经过仿真后可以看出各风电系统都能保持稳定性，只是系统稳定时的功角和转速与初始值差距比较大，如图14所示；双馈风机的转子电流越限，即机组本身不稳定，如图15所示。
图 14 无功负荷减少时电力系统功角曲线
Figure 14. The power angel with reduce of reactive load
图 15 无功负荷减少时双馈风机的转子电流曲线
Figure 15. The rotor current of the double-fed wind power system with reduce of reactive load
在有功负荷减少的过程中，含双馈风机的电力系统由稳定变为不稳定，机组本身也慢慢失稳，直驱风电系统和恒速风电系统也慢慢失稳。如图16描述的是在有功负荷由0.9、0.9和0.9个标幺值减少到0.6、0.6和0.6个标幺值时，双馈电力系统电压曲线。
图 16 有功负荷减少时系统电压变化
Figure 16. The change of the voltage with reduce of reactive load
-
本文在PSAT上建立了几种风机的数学模型，并进行了暂态稳定仿真，在比较的基础上得出了如下结论：
1）故障持续时间越长，系统越容易失稳；在短路情况下，双馈机组本身的稳定性较难保持，但系统稳定性保持较好。
2）双馈机组和直驱机组低电压穿越能力要好于恒速机组。
3）无功补偿容量的加入使得双馈风电系统和直驱风电系能较好的恢复稳定，但随着故障时间进一步增加，无功补偿的作用将不会很明显。
4）突然甩负荷对各风电系统的稳定性都有较大的影响。
5）同传统电力系统相比，含双馈风力发电机的电力系统受负荷影响较大，突然甩负荷使得双馈风电系统和机组本身很容易失稳，即双馈系统对负荷变化非常敏感。相对来说直驱风电系统和恒速风电系统的稳定性能较好保持，但机组本身也较容易失稳。
6）由于换流器的解耦作用，双馈风机和直驱风机暂态稳定性受发电机转速影响较小，受换流器电流大小较大的限制。因而变速风机的稳定性与换流器本身的特性有密切关系。
Research on Transient Stability for Wind Power System Based on PSAT
-
摘要:
目的 为了研究各风电系统的暂态稳定特性及各风电机组的低电压穿越能力，同时讨论无功补偿容量的加入对各风电系统稳定性的影响。 方法 在PSAT平台上建立了风速、换流器和风力发电机的数学模型，并利用三机九节点包含风力机组的电力系统网络模型，在三相短路故障情况下用时域仿真法对这个网络进行了仿真。 结果 获得了直驱风电系统、双馈风电系统、恒速恒频风电系统在不同母线短路故障持续时间下的电压曲线或功角曲线。 结论 得出了不同机组类型，不同扰动类型，以及不同持续时间下，各风电系统稳定性不同的结论。 Abstract:Introduction The paper aims to study the transient stability of the wind power system and that how the injection of reactive power to the wind power system affects the stability. Method So we established the mathematic model of the wind power generator， the frequency converter， and a three-machine-nine-node power system model was constructed based on PSAT. We realized the the time simulation in condition of the three phase fault. Result Then we get the voltage and power angel figure of direct drive wind power system，double-fed wind power system and constant speed wind power system under different bus short circuit fault time. Conclusion It demonstrates that the wind power system has different stability based on different wind power generator， disturbance and duration. -
Key words:
- PSAT /
- wind power system /
- transient stability /
- low voltage ride through
-
图 5 故障持续120 ms时直驱风电系统同步电机的σ-t曲线
Fig. 5 The σ-t figure of direct drive wind power system in 120 ms fault duration
图 6 故障持续120 ms时恒速机组的转速曲线
Fig. 6 The rotate speed figure of constant speed wind power system in 120 ms fault duration
图 7 故障持续246 ms时双馈风机的转子电流曲线
Fig. 7 The rotor current figure of double-fed wind power system in 246 ms fault duration
图 8 故障持续625 ms时恒速风电系统BUS7处电压
Fig. 8 The voltage of BUS7 of constant speed wind power system in 625 ms fault duration
图 9 故障持续625 ms时直驱风电系统BUS7处电压
Fig. 9 The voltage of BUS7 of direct drive wind power system in 625 ms fault duration
图 10 故障持续625 ms时双馈风电系统BUS7处电压
Fig. 10 The voltage of BUS7 of double-fed wind power system in 625 ms fault duration
图 11 故障持续625 ms加入SVC后恒速风电系统BUS7处电压
Fig. 11 The voltage of BUS7 of constant speed wind power system in 625 ms fault duration with SVC
图 12 故障持续625 ms加入SVC后直驱风电系统BUS7处电压
Fig. 12 The voltage of BUS7 of direct drive wind power system in 625 ms fault duration with SVC
图 13 故障持续625 ms加入SVC后双馈风电系统BUS7处电压
Fig. 13 The voltage of BUS7 of double-fed wind power system in 625 ms fault duration with SVC
图 15 无功负荷减少时双馈风机的转子电流曲线
Fig. 15 The rotor current of the double-fed wind power system with reduce of reactive load
表 1 三种风电机组及系统稳定性对比
Tab. 1.
Contrast of the stability between wind plants 类型 故障持续时/ms 100 246 625 双馈发电系统 稳定 稳定 不稳定 双馈机组本身 稳定 不稳定 不稳定 直驱发电系统 稳定 稳定 不稳定 直驱机组本身 稳定 稳定 不稳定 恒速发电系统 稳定 稳定 不稳定 恒速机组本身 稳定 不稳定 不稳定 
下载: 导出CSV
-
[1] MILANOF. Power system analysis toolbox documentation for PSAT version 2.1.6 ［EB /OL］. （2011-06-02）［2019-09-25］ https：//www.mendeley.com/catalogue/a456f9c8-56f2-3d5d-a3d8-245926956357/. [2] 王金铭，高键. 基于PSAT的电力系统时域仿真研究 ［J］. 科学技术与工程，2008，8（17）：4802-4805. WANGJ M，GAOJ. The time simulation study based on PSAT ［J］. Science Technology and Engineering，2008，8（17）：4802-4805. [3] 孙雷雷. PSAT在电力系统暂态稳定仿真中的应用 ［J］. 现代计算机（专业版），2009（10）：103-106. SUNL L. The power system transient stability simulation based on PSAT ［J］. Modern Computer（Professional），Modern Comupter，2009（10）：103-106. [4] 赵璐璐. 大型风电场并网暂态稳定性影响研究 ［D］. 北京：北京交通大学，2010. ZHAOL L. The power system transient stability study with large scale wind farm connecting to the grid ［J］. Beijing：Beijing Jiaotong University，2010. [5] 孙元章，吴俊，李国杰. 风力发电对电力系统的影响 ［J］. 电网技术，2007（20）：55-62. SUNY Z，WUJ，LIG J. The effect of the wind power generation to the power system ［J］. Power System Technology，2007（20）：55-62. [6] 陈雷，卢斯煜. 风电场出力特性与集群效应分析方法研究 ［J］. 南方能源建设，2017（1）：31-37. CHENL，LUS Y. Research on the analysis method for wind power generating output characteristic and cluster effects ［J］. Southern Energy Construction，2017（1）：31-37. [7] 郝元钊，李培强，李欣然.风电机组并网对电网暂态稳定性影响研究 ［D］. 长沙：湖南大学电气与信息工程学院，2011. HAOY Z，LIP Q，LIX R. The power system transient stability study with wind turbine connecting to the grid ［D］. Changsha：College of Electrical and Information Engineering，Hunan University，2011. [8] 李林，金先级. 数值计算方法 ［M］. 北京：高等教育出版社，2006. LIL，JINX J. Numerical computation method ［M］. Beijing：Higher Education Press，2006. -
点击查看大图
计量
- 文章访问数: 435
- HTML全文浏览量: 150
- PDF下载量: 37
- 被引次数: 0
