---

摘要：建立了包括汽轮机、发电机、调速器和励磁调节装置在内的全系统仿真模型  ，编制了仿真软件  ，在模拟发电机组上进行了三相突然短路故障引起的失步实验和仿真分析  ，实验结果验证了仿真结果的正确性 。对某电厂300 MW汽轮发电机进行了三相短路故障引起失步的砂石生产线仿真研究  ，并对三阻抗元件失步保护和改进型三阻抗元件失步保护进行了整定计算及分析 。
关键词：暂态稳定；振荡中心；整定计算；失单电机双辊机步保护

复合式制砂机1　引言

电力系统稳定性问题是指电力系统运行中受到扰动之后能否继续保持发电机间同步运行的问题 。根据扰动大小所确定的稳定问题的性质  ，把它分为静态稳定和暂态稳定 。所谓电力系统静态稳定性  ，一般是指电力系统在运行中受到微小扰动后  ，独立地恢复到它原来的运DZSF系列直线振动筛行状态的能力 。电力系统具有暂态稳定性  ，一般是指电力系统在正常运行时  ，受到一个大的扰动后  ，能从原来的运行状态  ，过渡到新的运行状态  ，并在新的运行状态下稳定地运行  ，而暂态不稳定的后果则是系统失去同步 。

引起电力系统大扰动的螺旋传输机原因主要有下列几种：

(1)负荷的突然变化  ，如投入输送网带或切除大容量的用户等；
(2)切除或投入系统的主要元件  ，皮带输送如发电机、变压器及线路等；
(3)发生短路故障 。提升机械

其中短路故障的扰动最为严重  ，常以喂料机此作为检验系统是否具有暂态稳定的依据 。

本文针对三相突然短路故障引起的失步过程  ，利用C＋＋语言编写了仿真程序石头制砂机  ，对具有调速器和励磁环节的汽轮发电机发生三相突然短路故障后失步的动态过程进行仿真  ，并通过自编的主设备整定计算软件进行整定  ，利用仿真和整定结果对三阻抗元件失步保护和改进型三阻抗元件失步保护进行分析 。

2　板链输送机数学模型

仿真中  ，发电机组采用矿山输送机械的数学模型为单机对无穷大系统  ，如图1所示 。

正常运行时  ，开关K断开；当K闭合时  ，发生发电机机端三相突然短料仓破拱路故障；当开关K再次断开时  ，故障切除 。

2．1　发电机的数学落砂机模型

采用Pvwin德赢网官方下载石破碎机ARK方程［1］对发电机进行数学建模 。建模过程中  ，除考虑发电机机端电压和机端电流的变化以及转子励磁电流和励磁电压的变化外  ，还考虑了发电机的功角δ和转差s的变化情况 。其矩阵形式如下：

式中H为转动惯量  ，T1为发电机原动机的输入转矩  ，电磁转矩  ，ω为发电机的同步角速度 。

2．2　输电系统的数学模型

输电系统在abc坐标系下的电压方程为

其中  ，rs为线路的正序电阻；xs是线路的正序电抗；rs0、 xs0是线路的零序电阻、电抗 。

2．3　调速系统的数学模型

为了控制原动机向发电机输出的机械功率  ，以及在并列运行的发电机之间合理地分配负荷  ，每一台原动机都配置了调速器 。原动机和调速器的组合构成了发电机的调速系统 。

国内外的大型汽轮机组  ，大都采用中间再热式汽轮机  ，由于一阶模型没有计及中间再热蒸汽容积效应  ，过于粗糙；三阶模型又有些繁琐  ，在本文的精度要求下  ，对汽轮机的仿真采用二阶模型  ，其传递函数框图如图2所示 。

在动态模拟实验室进行三相短路故障引起的失步实验时  ，模拟发电机组采用的原动机是直流电动机  ，因此在对模拟发电机组的三相突然短路及故障切除进行仿真计算时  ，采用了如图3所示的直流电动机数学模型 。

2．4　励磁调节系统的数学模型

励磁系统是发电机的重要组成部分  ，对电力系统及发电机本身的安全有很大的影响 。现代的大容量汽轮发电机广泛采用半导体励磁系统  ，常用的有它励静止硅整流器励磁系统、它励旋转半导体励磁系统(即无刷励磁)、自并励半导体励磁系统及自复励半导体励磁系统等 。

本文根据励磁调节系统的原理及动作特性  ，使用了一种简化的具有可控硅励磁调节器的交流励磁系统的数学模型  ，如图4所示 。

3　仿真和实验结果对比

本文对模拟发电机组发生三相短路故障引起的失步过程进行仿真  ，并同实验结果进行了对比 。

进行实验的模拟发电机组是隐极机  ，其额定参数：SN＝30 kV·A   ，UN＝400 V  ，IN＝43．3 A  ，cosφN＝0．8 。实验记录了发电机机端三相短路引起的失步  ，振荡中心［2］分别在发变组内和发变组外的实验结果 。下面列出振荡中心在发变组内和振荡中心在发变组外两种情况下的仿真结果和实验波形 。图5中  ，线路阻抗为rs＝0．036 61  ，xs＝0．418 44  ，故障前运行状况为U＝0．89  ，I＝0．67  ，cosφ＝0．99(滞后)(本文采用常用的可逆互感系数的标幺值系统)  ，发生三相短路故障后  ，经0．225 s切除故障的实验波形和仿真波形；图6中  ，线路阻抗为rs＝0．08  ，xs＝1． 528 44  ，故障前运行状况为U＝1．03  ，I＝0．45  ，cosφ＝0．90(滞后)  ，发生三相短路故障后  ，经0．24 s切除故障的实验波形和仿真波形 。

通过对实验波形和仿真波形的比较  ，可以看出所采用的数学模型是正确的  ，数值算法也能满足计算的精度要求 。振荡中心是系统振荡时电气量变化最剧烈的一块区域(或一点)  ，它的特点是电压跌落 最为严重  ，实验与仿真结果表明线路阻抗的大小与振荡中心的位置密切相关 。

4　某电厂300 MW汽轮发电机仿真计算

本节对实际电厂300 MW汽轮发电机组进行仿真  ，其额定参数为：

发电机正常运行时的工况为：U＝1．05  ，I＝1．0  ，cosφ＝0．9(滞后) 。仿真中均令t＝0．5 s时发生机端三相突然短路故障  ，分别取如下情况对其进行仿真：

(1)励磁调节器不动作  ，经0．3 s切除故障；
(2)励磁调节器不动作  ，经0．38 s切除故障；
(3)励磁调节器动作  ，经0．38 s切除故障；
(4)励磁调节器动作  ，经0．415 s切除故障；　　

通过仿真可以发现：情况(1)、(3)经过几次振荡后  ，振幅逐渐衰减  ，最后保持稳定；情况(2)、(4)中  ，发电机振荡愈来愈强烈  ，最终导致失步 。下面列出情况(2)和(3)的仿真结果(见图7  ，图8) 。

对比以上四种情况的仿真结果  ，可以看出在相同的短路持续时间0．38 s情况下  ，由于励磁调节器的自动调节作用  ，输出的励磁电源电压随着机端电压的大小变化进行调整  ，较好地维持发电机的机端电压  ，使得发电机最终没有失步  ，而没有励磁调节器的发电机在相同的故障条件下最终失步 。可见自动 励磁调节器对提高系统的暂态稳定性有一定的作用 。但是如果短路持续时间过长  ，发电机将会与系统失去同步 。经过仿真计算  ，在有、无励磁调节器的两种情况下  ，此单机对无穷大系统在故障下的极限切除时间分别为0．413 4 s和0．375 8 s 。

5　三阻抗元件失步保护

尽管目前的失步保护有很多种  ，但是发电机的功角变化状况和发电机本身能量的平衡能力是衡量发电机运行状况与稳定性的根本  ，发电机的失步保护与预测方案也都是围绕这两点进一步发展和演变的  ，下面用仿真的结果来分析三阻抗元件失步保护 。

5．1　三阻抗元件失步保护的原理

如图9所示单机—无穷大系统发生振荡时  ，假设各阻抗的阻抗角相等  ，则发电机与系统之间的功角δ是由机端测量阻抗Zj和发电机暂态阻抗以及发电机与系统的联系电抗Zst所组成的封闭阻抗三角形来决定 。各阻抗分别乘以电流  ，可得各电压向量  ，如图10所示  ，其中Eg与Es向量之间的相

镜主轴的的交点得到 。

5．2　三阻抗元件失步保护的整定计算

通过自编的发变组保护整定计算软件对失步保护进行整定  ，首先用一次侧公用参数将发电机、变压器、线路、系统电抗等都折合到机端电压级的标幺值  ，计算过程如下：

其中xt、xl、xs分别为变压器电抗、线路电抗和系统等值电抗  ，均为标幺值  ，Zn为阻抗基值  ，l为线路长度  ，单位km  ，x1为导线单位长度的电抗 。

透镜在R－X平面上第Ⅰ象限的最远点的阻抗ZA按机端至无限大系统之间的联系阻抗Zst进行整定：

透镜在R－X平面设计第Ⅲ象限的最大点的阻抗匹配ZB  ，按同样马达的纵轴暂态电抗来整定：

可以依照滑极中间继电器开关Ⅰ段和Ⅱ段脱扣区的必须  ，电抗ZC的整定一样 先选电力变压器电抗的90％：

透镜夹头(电位差器条直线)的倾倒角φZ一般来说选购

δset的整定：

为保证机组最大负荷时的最小阻抗rL．min位于透镜之外   ，以保证保护的动作选择性  ，为此选择

一般要求检测的最大滑差频率fs．max为5 Hz  ，所以αmax＝180°×(1－0．05fs．max)＝135°  ， 通过仿真计算  ，为保证故障切除后机端阻抗点位于透镜圆外  ，取δset＝125°  ，此时fs．max＝20×(1－α／180°)≈6．1 Hz 。

5．3　三阻抗元件失步保护的应用

下面应用5．2节整定的三阻抗元件失步保护  ，分别对第4节四种情况的仿真进行讨论 。阻抗轨迹与透镜主轴的夹角(即功角)为125°、180°、235°的时刻分别记为仿真中功角是直接通过的(在实际中可以在发电机转子上安装测速装置来实现)  ，与前面用阻抗轨迹间接计算功角的结果比较  ，见表1 。

对于(1)、(3)没有失步的情况  ，三阻抗元件失步保护没有动作；对于(2)、(4)失步的情况  ，由于阻抗轨迹穿过透镜  ，而且穿越的时间大于50 ms  ，所以三阻抗元件失步保护可靠动作 。

5．4　改进型三阻抗元件失步保护

三阻抗元件失步保护尽管有鉴别短路与振荡、失步振荡与同步摇摆的能力  ，但是它不具备预测功能  ，保护动作的结果只能是切机或解列  ，没有时间采取相应的抑制失步措施保持系统稳定 。鉴于此  ，文献［3］提出在此三阻抗元件失步保护的基础上加以改进  ，计算发电机与系统之间的功角的方法与三阻抗元件失步保护中的方法相同  ，由机端测量阻抗间接算出  ，然后通过某种算法来预测功角的变化轨迹 。　　

本文采用多项式回归模型对功角进行预测：为将来某一时刻功角δ(t)的预测值  ，a0、a1…am根据历史数据δ(Δt)…δ(N·Δt)  ，用最小二乘法按Am＝B(N)·δ(N)预测  ，其中

H(N)矩阵是一个常数阵  ，所以B(N)可事先离线计算出来 。这样多项式方程就转化为差分方程的形式  ，再根据预测时刻kΔt  ，计算出预测功角值δ′(kΔt) 。　　

应用5．2节整定的ZA、ZB值  ，对上述四种情况的仿真进行失步预测分析 。取预测时间为300 ms  ，每隔10 ms观察的12个数据组成历史数据窗口  ，连续3次预测出失步后(功角大于180°)发出告警信 号  ，设其时刻为t1  ，到达功角180°的时刻为t0  ，提前预测量为Δt 。预测过程中为提高预测计算精度  ，采用滚动预测  ，不断采用新数据  ，用最新的采样值更新预测数据窗的数据 。表2分别列举了间接和直接两种方法计算功角时的失步预测结果 。

对于仿真(1)、(3)不失步的情况  ，功角预测值始终没有大于180°  ，不满足动作判据  ，保护正确不动；对于仿真(2)、(4)失步的情况  ，均提前预测出功角大于180°  ，满足动作判据  ，保护正确动作 。可见改进的三阻抗元件失步保护能提前预测出发电机的失步情况  ，其时间提前量为抑制失步、保持发电机稳定提供可能 。

由于Zst和在失步过程中不是固定不变的  ，因此通过阻抗轨迹间接计算功角并不精确  ，理论上讲基于初相角和转差计算功角的方法更准确 。通过比较表1和表2  ，两种方法计算出来的结果相差不大 。

6　结论

为研究发电机三相短路引起失步的动态过程及其保护  ，本文建立了包括励磁调节器、调速器、汽轮机和发电机在内的发电机系统模型  ，实验和仿真结果表明数学模型和计算方法是正确的 。对某电厂的300 MW汽轮发电机组分不同情况进行仿真计算  ，说明了励磁调节器对提高系统暂态稳定性的作用  ，并利用仿真结果对三阻抗元件失步保护和改进型三阻抗元件失步保护进行整定计算和分析   ，验证了改进型三阻抗元件失步保护有较好的预测功能