【电力系统暂态分析】-第二章-同步发电机-前三节

第二章 同步发电机

一、同步发电机在空载情况下定子突然三相短路后的电流波形分析

1、同步发电机的组成及原理

组成

• 定子

  装有a、b、c三相电枢绕组

• 转子

  装有直流励磁绕组

原理

• 发电机励磁绕组中通入直流电流$i_f$，在转子周围产生磁场，定子a、b、c三相绕组切割磁力线，在电枢绕组感应出三相电动势

2、同步发电机发生三相短路

当同步发电机突然三相短路后，电枢绕组的电流将会对转子绕组产生影响，由于定、转子之间的电磁耦合从而造成同步发电机的暂态过程变得复杂

3、同步发电机空载情况下突然三相短路后的电流波形

同步发电机空载情况下，**三相短路后电流波形（图1）**如下图所示，（a）为三相定子电流波形（短路电流），（b）为励磁回路电流波形

由图1可知

• 在突然短路瞬间，定子短路电流和励磁回路电流均不突变，即三相定子电流均为0，励磁回路电流为$i_{f|0|}$（感性回路的电流不发生突变）
• 定子侧短路后三相绕组
  • 基频交流分量、直流分量和倍频分量（幅值较小）
• 转子侧短路后励磁绕组
  • 直流分量、基频交流分量
• 暂态过程中
  • 定子绕组的基频交流分量与转子中直流分量衰减时间常数相同
  • 定子侧直流分量与转子中基频交流分量衰减时间常数相同

上图的**短路电流波形的分解（图2）**如下图所示，（a）为三相直流分量，（b）为交流分量包络线的衰减

三相短路电流中的直流分量如图2（a）所示

• 三相直流分量大小不等
• 按照相同的指数规律衰减，最后衰减至0，直流分量的衰减时间常数：$T_a$

交流分量如图2（b）所示

• 逐渐衰减，最终衰减为稳态值$I_{m\infty }$，由图可知，短路初始值为$I_m^{\prime \prime }$，即$I_m^{\prime \prime }-I_{m\infty }$衰减为0
• 该衰减过程按两个时间常数衰减，分别为小的时间常数$T_d^{\prime \prime }$和大的时间常数$T_d^{\prime }$，一开始按照$T_d^{\prime \prime }$衰减，随后按照$T_d^{\prime }$衰减
• 交流分量幅值的表达式：$I_{\mathrm{m}}(t)=\left(I_{\mathrm{m}}^{\prime \prime }-I_{\mathrm{m}}^{\prime }\right){\mathrm{e}}^{-t/T_{\mathrm{d}}^{\prime \prime }}+\left(I_{\mathrm{m}}^{\prime }-I_{m\infty }\right){\mathrm{e}}^{-t/T_{\mathrm{d}}^{\prime }}+I_{m\infty }$
  • 次暂态分量：$I_{\mathrm{m}}(t)=\left(I_{\mathrm{m}}^{\prime \prime }-I_{\mathrm{m}}^{\prime }\right){\mathrm{e}}^{-t/T_{\mathrm{d}}^{\prime \prime }}$
    • 幅值：$I_{\mathrm{m}}^{\prime \prime }-I_{\mathrm{m}}^{\prime }$
    • 时间常数：$T_d^{\prime \prime }$（次暂态衰减）
  • 暂态分类：$\left(I_{\mathrm{m}}^{\prime }-I_{m\infty }\right){\mathrm{e}}^{-t/T_{\mathrm{d}}^{\prime }}$
    • 幅值：$I_{\mathrm{m}}^{\prime }-I_{m\infty }$
    • 时间常数：$T_d^{\prime }$（暂态衰减）
  • 稳态分量：$I_{m\infty }$
• 短路后过程：次暂态——暂态——稳态

二、同步发电机稳态运行情况及暂态参数

在进行分析过程中，认为同步发电机为理想电机（假设条件：电力系统暂态分析P25）

在同步发电机暂态过程分析中，主要研究凸极同步发电机，隐极机可作为凸极机的特例。对于凸极机的阻尼绕组，一般在转子的直轴（d轴）和交轴（q轴）上各设置一个等值阻尼绕组，分别为D绕组和Q绕组

1、定子、转子各绕组磁轴、电流的正方向

同步发电机各绕组轴线和电流正方向如下图所示：

如图可知：

• 定子三绕组分别为a、b、c，各相差$120^{\circ }$
  • 定子绕组各相绕组轴线的正方向为各相绕组磁链的正方向
• 转子极中心线用d轴（直轴）表示，极间轴线用q轴（交轴）表示，q轴超前d轴$90^{\circ }$
  • 励磁绕组和直轴阻尼绕组磁链的正方向与d轴正方向相同
  • 交轴阻尼绕组磁链的正方向与q轴正方向相同
• 励磁绕组f与d轴重合
• 阻尼绕组用两个互相正交的短接绕组等效，与d轴重合的称为D阻尼绕组，与q轴重合的称为Q阻尼绕组

定子各相绕组电流产生的磁通方向与各该相绕组轴线的正方向相反时电流为正值，转子各绕组电流产生的磁通方向与d轴或q轴正方向相同时电流为正值。

• 即定子绕组中正电流产生负磁通，转子绕组中正电流产生正磁通

2、磁链守恒原理及同步发电机双反应原理

磁链守恒原理

• 任何一个闭合的线圈，它的磁链在同一瞬间不能从一个数值跳变到另外一个数值，如果存在外界因素迫使线圈的磁链发生突变，该线圈会感应出一个自由电流分量，从而产生一个反方向磁链抵制外来磁链，以维持原线圈所匝链的磁链不变
• 是对发电机突然短路暂态过程进行物理分析的理论基础

双反应原理

• 电枢反应

  同步发电机空载时，定子电流为0，气隙中只有励磁电流建立的磁场。当同步发电机带负荷时，气隙中除了转子磁场外，还存在着由定子三相电流产生的同步旋转的电枢磁动势，从而气隙中的磁动势变成为合成磁动势，从而使气隙中原有的磁动势大小和位置都发生改变

• 前提

  对于凸极机，由于结构不对称，同一电枢磁动势作用在不同位置时电枢反应不一样

• 双反应理论

  当电枢磁动势$F_a$作用于交、直轴间的任意位置时，可将其分解为直轴电枢反应分量$F_{ad}$和交轴电枢反应分量$F_{aq}$，分别求出后再进行叠加

3、同步发电机的参数

a、同步发电机的稳态参数

• 空载电动势$E_q$
  • 励磁电流$i_f$产生的主磁通${\phi }_0$切割定子绕组，感应出空载电动势
  • 当磁路不饱和时，$E_q$与$i_f$是线性关系，$E_q$正比于$i_f$，即励磁电流的大小决定了空载电动势的大小
  • 方向沿q轴时，相量表达式：$\dot{E_q}=j{E}_q$
• 同步电抗
  • 直轴同步电抗（$X_d$）
    • 组成：$X_d=X_{ad}+X_{\sigma }$ （直轴电枢反应电抗$X_{ad}$、定子漏抗$X_{\sigma }$）
    • 反应d轴方向磁路的磁导，路径包括转子铁心、气隙和定子铁心
    • 等效电路：$X_{ad}$与$X_{\sigma }$串联
  • 交轴同步电抗
    • 组成：$X_q=X_{aq}+X_{\sigma }$（交轴电枢反应电抗$X_{aq}$、定子漏抗$X_{\sigma }$）
    • 反应q轴方向磁路的磁导，路径无转子绕组（忽略阻尼绕组）
    • 等效电路：$X_{aq}$与$X_{\sigma }$串联
  • $X_q<X_d$
• 电枢反应电抗
  • 直轴电枢反应电抗（$X_{ad}$）
    • 由d轴电枢反应磁通${\phi }_{ad}$产生，等效为去磁效应
    • 满足$E_{aq}=I_d{X}_{ad}$（$E_{ad}$为电枢反应引起的压降）
  • 交轴电枢反应电抗（$E_{a}q$）
    • 由q轴电枢反应磁通${\phi }_{aq}$产生，无励磁绕组影响
    • 满足$E_{aq}=I_q{X}_{aq}$
• 定子漏抗（$X_{\sigma }$）
  • 定子绕组自身电流产生的漏磁通${\phi }_{\sigma }$对应的电抗
  • 与电枢反应无关，仅与绕组结构有关
  • $X_d和X_q$中均包含$X_{\sigma }$
• 端电压方程
  • 纯感性负载（$I_d\ne 0,I_q=0$）
    • $\dot{U_q}=\dot{E_q}-j{I}_d{X}_d$
  • 纯阻性负载（$I_d=0,I_q\ne 0$）
    • $\dot{U_d}=-j{I}_q{X}_q$
  • 一般负载
    • 需叠加d、q轴分量
• 励磁绕组参数（$X_f$）
  • $X_f=X_{f\sigma }+X_{ad}$（励磁绕组漏抗$X_{f\sigma }$、d轴电枢反应电抗$X_{ad}$）
  • 反应励磁电流$i_f$产生的总磁链（主磁通+漏磁通）

b、同步发电机正常稳态运行时的电压平衡关系及相量图

同步发电机正常稳态运行时的相量图如下图所示:

图（a）为凸极机，图（b）为隐极机

图中参数：

$\dot{E_q}$：空载电动势；$\dot{U}$：发电机端电压；$\dot{I}$：电流；$\delta$：$\dot{E_q}、\dot{U}$间的夹角，即功角；$\phi$：$\dot{U}、\dot{I}$间的夹角，即功率因数角

$\dot{U_d}、\dot{U_q}、\dot{I_d}、\dot{I_q}$：电压和电流在d、q轴上的分量

$x_d、x_q$：定子直轴和交轴同步电抗

$r$：定子绕组电阻

计及电阻后，q、d轴电压方程如下所示：

$\{\begin{array}{l}{\dot{E}}_q={\dot{U}}_q+r{\dot{I}}_q+j{x}_d{\dot{I}}_d\\ 0={\dot{U}}_d+r{\dot{I}}_d-j{x}_q{\dot{I}}_q\end{array}$

且：$\dot{U}={\dot{U}}_d+{\dot{U}}_q$、$\dot{I}={\dot{I}}_d+{\dot{I}}_q$，将q、d轴电压方程相加，可得图（a）和图（b）相量图所对应的方程，分别为：

${\dot{E}}_q=\dot{U}+r\dot{I}+j{x}_d{\dot{I}}_d+j{x}_q{\dot{I}}_q$（凸极机$x_d\ne x_q$）

${\dot{E}}_q=\dot{U}+r\dot{I}+j{x}_d\dot{I}$（隐极机$x_d=x_q$）

当忽略$r$后，方程分别为：

${\dot{E}}_q=\dot{U}+j{x}_d{\dot{I}}_d+j{x}_q{\dot{I}}_q$（凸极机$x_d\ne x_q$）

${\dot{E}}_q=\dot{U}+j{x}_d\dot{I}$（隐极机$x_d=x_q$）

对于以上两式，$\dot{U}、\dot{I}$可知，分别为$\dot{U}=U\angle 0^{\circ }、\dot{I}=I\angle -\phi$，但其他数据未知

• 对于隐极机

  可根据$\dot{U}、\dot{I}$直接求出$\dot{E_q}$，进而确定$\dot{U_d}、\dot{U_q}、\dot{I_d}、\dot{I_q}$

• 对于凸极机

  不能够直接求出，需借助虚构电动势$\dot{E_Q}$

  可将凸极机电压公式改写成：

  $\begin{array}{rl}{\dot{E}}_q& =\dot{U}+j{x}_d{\dot{I}}_d+j{x}_q{\dot{I}}_q\\ & =\dot{U}+j{x}_d{\dot{I}}_d+j{x}_q{\dot{I}}_q+j{x}_q{\dot{I}}_d-j{x}_q{\dot{I}}_d\\ & =\dot{U}+j{x}_q({\dot{I}}_q+{\dot{I}}_d)+j(x_d-x_q){\dot{I}}_d\\ & =\dot{U}+j{x}_q\dot{I}+j(x_d-x_q){\dot{I}}_d\\ & ={\dot{E}}_Q+j(x_d-x_q){\dot{I}}_d\end{array}$

  其中：

  ${\dot{E}}_Q=\dot{U}+j{x}_q\dot{I}$

  由此即可得出虚构电动势的相量图，相量图如下图所示：

  

  由图可知，虚构电动势$\dot{E_Q}$的方向在q轴上，故可通过求出$\dot{E_Q}$来确定q轴的方向，进而求出其他待求量。

c、忽略阻尼绕组时的发电机暂态电动势和暂态电抗

• 短路瞬间过程分析

  三相短路后，定子电流交流分量增大 ——> 定子绕组旋转磁动势增大 ——> 定子增大的磁通试图穿过励磁绕组（匝链励磁绕组） ——> 由于磁链守恒，励磁绕组感应出反向直流电流$\Delta i_{f\alpha }$ ——> 总的励磁电流$i_{f|0|}+\Delta i_{f\alpha }$ ——> 对应的空载电动势$E_{q0}$突变，突变量$\Delta E_{q0}=x_{ad}\Delta i_{f\alpha }$，抵制磁通的穿入 ——> 迫使磁路路径发生改变，改走漏磁路径
  $|0|$：短路前瞬间，$0$：短路后瞬间

• 参数变化

  • 励磁电流突变 ——> 暂态电动势$E_q^{{}^{\prime }}$
  • 电枢反应磁通路径的变化 ——> 暂态电抗$X_d^{\prime }$

• 暂态电动势$E_q^{’}$

  • 短路前后保持不变
  • 用于计算短路电流初始量
  • $E_q^{’}=E_{q0}+\Delta E_{q0}，\Delta E_{q0}=x_{ad}\Delta i_{f\alpha }$

• 暂态电抗$X_d^{’}$

  • $X_d^{’}=x_{\sigma }+\frac{x_{f\sigma }{x}_{ad}}{x_{f\sigma }+x_{ad}}=x_{\sigma }+x_{ad}^{{}^{\prime }}$
    • $x_{\sigma }$：定子漏抗
    • $x_{ad}^{’}$：暂态直轴电枢反应电抗
  • 交轴暂态电抗 —— $x_q^{’}=x_q$

d、计及阻尼绕组时发电机的次暂态电动势和次暂态电抗

因为计及阻尼绕组，所以转子d轴上有两个绕组，一个为励磁绕组，一个为阻尼绕组

• 物理过程分析

  三相短路后，定子电流突增产生的磁链会同时匝链励磁绕组和阻尼绕组 ——> 由于磁链守恒，励磁绕组和阻尼绕组共同维持磁链不变 ——> 从而定子电枢反应磁通被迫绕开转子绕组，改走漏磁路径

• 参数变化

  • 转子绕组磁链守恒 ——> 次暂态电动势$E_q^{”}、E_d^{"}$
  • 电枢反应路径受限 ——> 次暂态电抗$x_d^{”}、x_q^{”}$

• 次暂态电动势

  • 直轴次暂态电动势${\dot{E}}_{d0}^{”}$

    ${\dot{E}}_{d0}^{”}={\dot{E}}_{d|0|}^{”}$

    ${\dot{E}}_{d0}^{”}={\dot{U}}_{d|0|}+j{\dot{I}}_{q|0|}{x}_q^{”}$

  • 交轴次暂态电动势${\dot{E}}_{q0}^{”}$

    ${\dot{E}}_{q0}^{”}={\dot{E}}_{q|0|}^{”}$

    ${\dot{E}}_{q0}^{”}={\dot{U}}_{q|0|}+j{\dot{I}}_{d|0|}{x}_d^{”}$

• 次暂态电抗

  • 直轴次暂态电抗$E_d^{”}$

    $E_d^{”}=x_{\sigma }+\frac{1}{\frac{1}{x_{ad}}+\frac{1}{x_{f\sigma }}+\frac{1}{x_{D\sigma }}}$

    定、转子等效为三绕组变压器，如下图所示：

    

    图（a）为等值三绕组变压器，图（b）为等效电路

  • 交轴次暂态电抗$E_q^{”}$

    $X_q^{\prime \prime }=X_{\sigma }+\frac{X_{aq}{X}_{Q\sigma }}{X_{aq}+X_{Q\sigma }}$

    定子和阻尼绕组可等效为双绕组变压器，如下图所示：

    

    图（a）为等值三绕组变压器，图（b）为等效电路

三、同步发电机空载突然三相短路电流分析

1、同步发电机突然短路后定子电流的分析

短路瞬间物理过程（不考虑发电机的阻尼绕组）

• 短路前状态

  发电机空载运行（定子电流为0），励磁电流产生主磁通，定子绕组感应出空载电动势。

• 短路后状态

  • 定子绕组产生三相交流分量（$i_{a\omega }、i_{b\omega }、i_{c\omega }$），用于抵消主磁通的变化，阻止磁链突变

  • 定子绕组产生三相直流分量（$\Delta i_{a\alpha }、\Delta i_{b\alpha }、\Delta i_{c\alpha }$），三相幅值不等但相加为零，用于维持短路瞬间定子绕组磁链不变

  • 交流分量与主磁链反向（去磁效应），直流分量与初始磁链反向（维持守恒）

• 二倍频波动

  定子三相绕组中直流分量所合成的静止磁动势与转子的直轴与交轴的磁阻不相同，频率两倍于基频，即直流电流不是恒定的，幅值以两倍频波动

• 短路瞬间电流为零

2、同步发电机突然短路后励磁绕组电流的分析

• 短路前状态

  • 定子绕组：空载，电流为零
  • 励磁绕组：存在恒定的励磁电流$i_{f|0|}$，产生空载电动势$E_{q0}$

• 短路后状态

  • 定子绕组
    • 交流分量$i_{i\omega }$：三相对称，产生与转子同步旋转的电枢反应磁动势，
    • 直流分量$\Delta i_{i\alpha }$：三相不对称，合成静止磁动势，被励磁绕组以${\omega }_0$转速旋转切割
    • 倍频分量$\Delta i_{i2\omega }$：由转子侧交流分量$\Delta i_{f\omega }$的脉动磁场感应产生
  • 励磁绕组（转子）
    • 直流分量$\Delta i_{f\alpha }$：阻止d轴电枢反应磁链${\psi }_{ad}$穿过励磁绕组，维持励磁绕组磁链守恒
    • 交流分量$\Delta i_{f\omega }$：单相基频，因磁链绕组以${\omega }_0$转速切割磁动势
  • 磁链守恒方程
    • ${\psi }_{fi}=-({\psi }_{ad}+{\psi }_{f\omega })$
    • ${\psi }_{ad}$：d轴电枢反应磁链；${\psi }_{f\omega }$：定子电流分量引起的转子交流分量磁链

• 电流分量关系

  定子电流与转子电流各分量间的关系如下图所示：

  

  各对应分量相互依存、相互影响

3、短路过程中短路电流各分量的衰减规律

• 短路电流分量的组成
  • 定子侧
    • 交流分量$i_{i\omega }$
    • 直流分量$\Delta i_{i\alpha }$
    • 倍频分量$\Delta i_{i2\omega }$
  • 转子侧
    • 直流分量$\Delta i_{f\alpha }$
    • 交流分量$\Delta i_{f\omega }$
• 衰减规律与时间常数
  • 定子侧直流分量$\Delta i_{i\alpha }$衰减
    • 无源直流分量，流过电阻消耗能量，最终衰减至0
    • 时间常数$T_a$：取决于定子回路的电阻和从定子侧看向转子侧的等值电感，$T_a=\frac{L_a}{R_a}$
    • $\Delta i_{i\alpha }$衰减 ——> $\Delta i_{f\omega }$和$\Delta i_{i2\omega }$同步衰减
  • 转子直流分量$\Delta i_{f\alpha }$衰减
    • 励磁绕组电阻导致能量耗散
    • 时间常数$T_d^{’}$：取决于转子回路的电阻和从转子侧看向定子侧的等值电感，$T_d^{’}=\frac{L_f^{’}}{R_f}$
    • $\Delta i_{f\alpha }$衰减 ——> $i_{i\omega }$同步衰减

4、短路电流基频交流分量的初始值及稳态有效值

• 短路电流稳态值$I_{\infty }$

  • 当短路达到稳态后，各自由分量衰减完毕，仅剩恒定的励磁电流$i_{|f0|}$产生强制基频分量
  • $I_{\infty }=I_d=\frac{E_q}{x_d}$，其中：$E_q$：空载电动势，$x_d$：直轴同步电抗

• 短路电流的初始值

  • 忽略阻尼绕组时的初始值$I^{’}$（暂态电流）

    • $\dot{I^{’}}=\dot{I_d^{’}}=\frac{{\dot{E}}_q^{’}}{j{x}_d^{’}}$

      $I^{’}=I_d^{’}=\frac{E_q^{’}}{x_d^{’}}$，其中：$E_q^{{}^{\prime }}$：暂态电动势，$x_d^{’}$：直轴暂态电抗

    • $x_d^{’}<x_d$ ——> $I^{’}>I_{\infty }$

  • 计及阻尼绕组时的初始值$I^{”}$（次暂态电流）

    • $\dot{I^{”}}=\dot{I_d^{”}}=\frac{{\dot{E}}_q^{”}}{j{x}_d^{”}}$

      $I^{"}=I_d^{"}=\frac{E_q^{"}}{x_d^{"}}$，其中：$E_q^{”}$：次暂态电动势，$x_d^{”}$：直轴次暂态电抗

    • $X_d’’ < X_d’ < X_d $ ——> $I^{\prime \prime }>I^{\prime }>I_{\infty }$

5、空载时短路电流的近似表达式

基频交流分量的近似表达式

• 短路电流交流分量幅值随时间的变化过程
  $I_m(t)=\sqrt{2}{E}_{q|0|}\left[\left(\frac{1}{X_d^{\prime \prime }}-\frac{1}{X_d^{\prime }}\right)e^{-t/T_d^{\prime \prime }}+\left(\frac{1}{X_d^{\prime }}-\frac{1}{X_d}\right)e^{-t/T_d^{\prime }}+\frac{1}{X_d}\right]$

  • 次暂态项：$\left(\frac{1}{X_d^{\prime \prime }}-\frac{1}{X_d^{\prime }}\right)e^{-t/T_d^{\prime \prime }}$
  • 暂态项：$\left(\frac{1}{X_d^{\prime }}-\frac{1}{X_d}\right)e^{-t/T_d^{\prime }}$
  • 稳态项：$\frac{1}{X_d}$

• 三相交流电流瞬时值

  $\{\begin{array}{l}{i}_{a\omega }=I_m(t)\cos ({\theta }_0+{\omega }_{0}t)\\ i_{b\omega }=I_m(t)\cos ({\theta }_0+{\omega }_{0}t-120^{\circ })\\ i_{c\omega }=I_m(t)\cos ({\theta }_0+{\omega }_{0}t+120^{\circ })\end{array}$

全电流的近似表达式

$\{\begin{array}{l}{i}_a=\sqrt{2}{E}_{q|0|}\left[\left(\frac{1}{x_d^{\prime }}-\frac{1}{x_d^{\prime }}\right)e^{-t/T_d^{\prime }}+\left(\frac{1}{x_d^{\prime }}-\frac{1}{x_d}\right)e^{-t/T_d^{\prime }}+\frac{1}{x_d}\right]\cos ({\theta }_0+{\omega }_{0}t)-\frac{\sqrt{2}{E}_{q|0|}}{\frac{n}{x_d}}\cos {\theta }_0{e}^{-t/T_a}\\ i_b=\sqrt{2}{E}_{q|0|}\left[\left(\frac{1}{x_d^{\prime }}-\frac{1}{x_d^{\prime }}\right)e^{-t/T_d^{\prime }}+\left(\frac{1}{x_d^{\prime }}-\frac{1}{x_d}\right)e^{-t/T_d^{\prime }}+\frac{1}{x_d}\right]\cos ({\theta }_0+{\omega }_{0}t-120^{\circ })-\frac{\sqrt{2}{E}_{q|0|}}{\frac{n}{x_d}}\cos ({\theta }_0-120^{\circ })e^{-t/T_a}\\ i_c=\sqrt{2}{E}_{q|0|}\left[\left(\frac{1}{x_d^{\prime }}-\frac{1}{x_d^{\prime }}\right)e^{-t/T_d^{\prime }}+\left(\frac{1}{x_d^{\prime }}-\frac{1}{x_d}\right)e^{-t/T_d^{\prime }}+\frac{1}{x_d}\right]\cos ({\theta }_0+{\omega }_{0}t+120^{\circ })-\frac{\sqrt{2}{E}_{q|0|}}{\frac{n}{x_d}}\cos ({\theta }_0+120^{\circ })e^{-t/T_a}\end{array}$