CN104331620A - 近区密集火电多直流外送系统次同步振荡的分析方法 - Google Patents

Abstract

一种近区密集火电多直流外送系统次同步振荡的分析方法，包含有以下步骤：步骤1、获取近区密集火电多直流外送系统的静态参数；步骤2、利用所述静态参数，筛选出存在次同步振荡风险的机组，作为待研究机组；步骤3、构建网络部分的复频域方程和每个所述待研究机组的发电机电磁部分的复频域模型；步骤4、根据所述复频域方程和所述复频域模型计算出每个所述待研究机组的复转矩系数；步骤5、根据每个所述待研究机组的复转矩系数，分析出该待研究机组的存在次同步振荡风险的频率。无需进行复杂的拓扑分析生成网络常态树来确定网络的状态变量，无需按照严格的规则和正方向规定对支路和节点进行编号，可方便的生成网络部分状态空间，计算复杂度小。

Description

近区密集火电多直流外送系统次同步振荡的分析方法

技术领域

本发明电力系统风险检测技术领域，尤其是近区密集火电多直流外送系统次同步振荡的分析方法。

背景技术

随着“西电东送”的推进，大型火电基地经串联电容补偿线路、高压直流输电(HVDC)的点对网输电已成为我国重要的输电模式。上述输电模式在提高系统输电能力的同时，也可能会造成电力系统的次同步振荡。近年来，我国一些火电基地均不同程度发生了次同步振荡问题，系统包括串联电容补偿或直流远距离输电工程且可能涉及临近的不同机组和电厂，次同步振荡特性复杂，不能简单的用单机简化系统对其进行特性模拟和研究。因此，需要对近区密集火电多直流外送系统的次同步振荡特性进行深入的分析和研究。特征值分析法是研究次同步振荡问题的精确的数学方法。特征值分析法能够明确系统结构、参数等因素对次同步振荡阻尼的影响；同时，结合控制理论，特征值分析法也可为次同步振荡抑制措施控制器的设计提供理论基础和方法指导。但目前特征值分析的研究对象主要为单机系统，复杂系统的SSO特征值分析还存在一定的困难。

电力系统次同步振荡分析需计及网络元件的动态特性，，而实际电网中元件的接线较为复杂，且网络独立状态变量的选取不易，造成复杂电网的网络部分状态空间建立十分复杂。目前的研究主要建立在对网络进行良好的拓扑分析获得网络的常态树的基础上，然后选取树枝电容电压和连枝电感电流为网络部分的状态变量，并按照严格的规则和正方向规定对节点和支路进行编号排序形成网络基本割集和基本回路矩阵，最后，利用全网的KCL(Kirchhoff's Circuit Laws)、KVL(Kirchhoff's Voltage Laws)消去除状态变量和输入、输出变量外的其他变量，得到网络部分的状态空间；由于复杂系统网络常态树的生成十分复杂，且全网的KCL和KVL运算量大，因此，该方法应用于复杂系统的次同步振荡分析具有一定的实现难度。此外，在系统网络结构变化时，上述方法需重新对网络变量进行排序，重新建立模型网络元件的动态特性矩阵、KCL/KVL系数矩阵，计算复杂度高。

发明内容

本发明的目的之一是提供一种近区密集火电多直流外送系统次同步振荡的分析方法，以解决现有技术中机组的次同步振荡风险检测复杂的问题。

在一些说明性实施例中，所述近区密集火电多直流外送系统次同步振荡的分析方法，包括：

步骤1、获取近区密集火电多直流外送系统的静态参数；

步骤2、利用所述静态参数，筛选出存在次同步振荡风险的机组，作为待研究机组；

步骤3、构建网络部分的复频域方程和每个所述待研究机组的发电机电磁部分的复频域模型；

步骤4、根据所述复频域方程和所述复频域模型计算出每个所述待研究机组的复转矩系数；

步骤5、根据每个所述待研究机组的复转矩系数，分析出该待研究机组的存在次同步振荡风险的频率。

与现有技术相比，本发明的说明性实施例包括以下优点：

利用基于复频域端口等效导纳矩阵(CPCM)的复转矩系数法计算机组复转矩系数，分析其次同步振荡阻尼特性。该方法无需进行复杂的拓扑分析生成网络常态树来确定网络的状态变量，无需按照严格的规则和正方向规定对支路和节点进行编号，可方便的生成网络部分状态空间，计算复杂度小；在网络结构发生变化时可快速的对网络部分状态空间进行修改。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1是按照本发明的说明性实施例的流程图；

图2是按照本发明的说明性实施例的三相对称的R-L串联支路和C支路示意图；

图3是本发明实施例中#1机组电气弹性系数和电气阻尼系数频域图；

图4是本发明实施例中#1机组机械弹性系数和机械阻尼系数频域图；

图5是本发明实施例中#1机组综合弹性系数和综合阻尼系数频域图；

图6是本发明实施例中#2机组电气弹性系数和电气阻尼系数频域图；

图7是本发明实施例中#2机组机械弹性系数和机械阻尼系数频域图；

图8是本发明实施例中#2机组综合弹性系数和综合阻尼系数频域图。

具体实施方式

在以下详细描述中，提出大量特定细节，以便于提供对本发明的透彻理解。但是，本领域的技术人员会理解，即使没有这些特定细节也可实施本发明。在其它情况下，没有详细描述众所周知的方法、过程、组件和电路，以免影响对本发明的理解。

为了更好的理解本发明中的说明性实施例，下面对本发明说明性实施例中的主要思想进行简单说明。

本发明通过对现有各种次同步振荡分析方法进行深入分析，在已有的调研基础上，改进了复转矩系数法，配合机组作用系数法提出了一种新的近区密集火电多直流外送系统次同步振荡分析方法。该方法无需进行复杂的拓扑分析生成网络常态树来确定网络的状态变量，无需按照严格的规则和正方向规定对支路和节点进行编号，可方便的生成网络部分状态空间，计算复杂度小；在网络结构发生变化时可快速的对网络部分状态空间进行修改。该方法采用分散消元的方式来克服的“维数灾”问题，且节点分析的方法可方便的建立全系统的复频域模型，因此较适于进行多机电力系统的次同步振荡分析方法采用分散消元的方式来克服的“维数灾”问题，且节点分析的方法可方便的建立全系统的复频域模型，因此较适于进行多机电力系统的次同步振荡分析，并可以分析研究系统中SVC等具有发电机转速反馈控制的SSO阻尼控制装置对机组轴系扭振特性的影响。该方法可为火电厂的系统接入方案规划以及SSO抑制措施的设计提供指导和理论支撑，对我国火电基地与直流输电系统的规划建设具有重要意义。

在一些说明性实施例中，公式中所使用的符号A、B、E、F、G、H、M均为推导过程中用到的系数矩阵参数，并没有具体意义。

如图1，公开了一种近区密集火电多直流外送系统次同步振荡的分析方法，其特征在于，包含有以下步骤：

步骤1、获取近区密集火电多直流外送系统的静态参数；

步骤2、利用所述静态参数，筛选出存在次同步振荡风险的机组，作为待研究机组；

步骤3、构建网络部分的复频域方程和每个所述待研究机组的发电机电磁部分的复频域模型；

步骤4、根据所述复频域方程和所述复频域模型计算出每个所述待研究机组的复转矩系数；

步骤5、根据每个所述待研究机组的复转矩系数，分析出该待研究机组的存在次同步振荡风险的频率。

在一些说明性实施例中，所述步骤1中获取到的所述静态参数至少包括：近区密集火电多直流外送系统的线路参数、机组参数以及控制参数。

在一些说明性实施例中，所述步骤2中采用机组作用系数法筛选出存在次同步振荡风险的机组。

在一些说明性实施例中，所述采用机组作用系数法筛选出存在次同步振荡风险的机组，具体包括：

1)、利用如下公式计算出每个机组的机组作用系数；

${\mathrm{UIF}}_i=\frac{S_{\mathrm{HVDC}}}{S_i}{(1-\frac{{\mathrm{SC}}_i}{{\mathrm{SC}}_{\mathrm{\Σ}}})}^2$

其中，S_HVDC为高压直流输电的额定容量，S_i为发电机组的额定容量，SC_i为不计机组i影响时整流站交流母线的短路容量，SC_∑为计算机组i影响时整流站交流母线的短路容量；

2)、当UIF_i>0.1时，判定第i个机组为存在次同步振荡风险的机组。

当汽轮发电机组i与交流系统的电气联系很强或距离高压直流输电整流站的电气距离很远时，其对整流站交流母线短路容量的贡献很小，即有SC_i≈SC_∑，进一步可以得到UIF_i≈0。若对第i台机组有UIF_i<0.1，则一般认为该机组与高压直流输电的耦合程度较弱，不存在发生由高压直流输电引起的次同步振荡问题。在计算连接于同一母线的多台相同型号机组的机组作用系数UIF时，须将其按容量合并为一台等值机组进行。

在一些说明性实施例中，所述步骤3中构建网络部分的复频域方程的过程，具体包括：

利用三相对称的R-L串联支路和C支路，如图2所示；构建系统等值网络及所述系统等值网络中的每个网络元件的复频域端口等效导纳矩阵；

利用网络元件的复频域端口等效导纳矩阵，采用节点分析法构建网络部分的复频域方程：

I＝Y_RLC(s)U

其中，I为支路的线性化电流矩阵，U为支路的线性化电压矩阵，Y_RLC(s)为支路的等效矩阵。

在一些说明性实施例中，构建每个网络元件的复频域端口等效导纳矩阵，具体包括：

建立在xy同步坐标系下的动态方程：

$\left[\begin{array}{c}\mathrm{\&Delta;}{u}_{\mathrm{ijx}}\\ \mathrm{\&Delta;}{u}_{\mathrm{ijy}}\end{array}\right]=\left[\begin{array}{cc}{x}_l& 0\\ 0& x_l\end{array}\right]s\left[\begin{array}{c}\mathrm{\&Delta;}{i}_{\mathrm{ijx}}\\ \mathrm{\&Delta;}{i}_{\mathrm{ijy}}\end{array}\right]+\left[\begin{array}{cc}r& -x_l\\ x_l& r\end{array}\right]\left[\begin{array}{c}\mathrm{\&Delta;}{i}_{\mathrm{ijx}}\\ \mathrm{\&Delta;}{i}_{\mathrm{ijy}}\end{array}\right]$

$\left[\begin{array}{c}\mathrm{\&Delta;}{i}_{\mathrm{Cx}}\\ \mathrm{\&Delta;}{i}_{\mathrm{Cy}}\end{array}\right]=\frac{1}{x_C}\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]s\left[\begin{array}{c}\mathrm{\&Delta;}{u}_{\mathrm{Cx}}\\ \mathrm{\&Delta;}{u}_{\mathrm{Cy}}\end{array}\right]+\frac{1}{x_C}\left[\begin{array}{cc}0& -1\\ 1& 0\end{array}\right]\left[\begin{array}{c}\mathrm{\&Delta;}{u}_{\mathrm{Cx}}\\ \mathrm{\&Delta;}{u}_{\mathrm{Cy}}\end{array}\right]$

其中，I＝[Δi_x Δi_y]^T，U＝[Δu_x Δu_y]^T，Δ表示变量相对于其稳态值的微变量，x_l为网络元件电抗，x_c为网络元件电容，r为网络元件电阻，s为拉普拉斯算子。

在一些说明性实施例中，构建每个所述待研究机组的发电机电磁部分的复频域模型，具体包括：

1)、将所述建立在xy同步坐标系下的动态方程线性转换为在同步发电机dq坐标系下的动态方程；

其中，所述转换的关系如下：

$\left[\begin{array}{c}\mathrm{\&Delta;}{X}_x\\ \mathrm{\&Delta;}{X}_y\end{array}\right]=\left[\begin{array}{cc}\sin {\mathrm{\&delta;}}_{G0}& \cos {\mathrm{\&delta;}}_{G0}\\ -\cos {\mathrm{\&delta;}}_{G0}& \sin {\mathrm{\&delta;}}_{G0}\end{array}\right]\left[\begin{array}{c}\mathrm{\&Delta;}{X}_d\\ \mathrm{\&Delta;}{X}_q\end{array}\right]+\left[\begin{array}{cc}\cos {\mathrm{\&delta;}}_{G0}& -\sin {\mathrm{\&delta;}}_{G0}\\ \sin {\mathrm{\&delta;}}_{G0}& \cos {\mathrm{\&delta;}}_{G0}\end{array}\right]\left[\begin{array}{c}{X}_{d0}\\ X_{q0}\end{array}\right]\mathrm{\&Delta;}{\mathrm{\&delta;}}_G$

$\left[\begin{array}{c}\mathrm{\&Delta;}{X}_d\\ \mathrm{\&Delta;}{X}_q\end{array}\right]=\left[\begin{array}{cc}\sin {\mathrm{\&delta;}}_{G0}& -\cos {\mathrm{\&delta;}}_{G0}\\ \cos {\mathrm{\&delta;}}_{G0}& \sin {\mathrm{\&delta;}}_{G0}\end{array}\right]\left[\begin{array}{c}\mathrm{\&Delta;}{X}_x\\ \mathrm{\&Delta;}{X}_y\end{array}\right]+\left[\begin{array}{cc}\cos {\mathrm{\&delta;}}_{G0}& \sin {\mathrm{\&delta;}}_{G0}\\ -\sin {\mathrm{\&delta;}}_{G0}& \cos {\mathrm{\&delta;}}_{G0}\end{array}\right]\left[\begin{array}{c}{X}_{x0}\\ X_{y0}\end{array}\right]\mathrm{\&Delta;}{\mathrm{\&delta;}}_G$

其中，变量X可为发电机机端电压u或电流i，下标“0”表示变量的稳态值；

2)、根据在同步发电机dq坐标系下的动态方程构建每个所述待研究机组的发电机电磁部分的线性化状态空间：

$s\left[\begin{array}{c}{X}_G\\ Z_G\end{array}\right]=A_G\left[\begin{array}{c}{X}_G\\ Z_G\end{array}\right]+B_G{U}_G+B_G^{\&prime;}\mathrm{\&Delta;}{\mathrm{\&omega;}}_G$

其中，X_G＝[Δi_d Δi_q Δi_f Δi_D Δi_Q Δi_g]^T为发电机定子/转子绕组电流矩阵，U_G＝[Δu_d Δu_q]^T为发电机的机端电压矩阵；Z_G为除X_G外发电机电磁部分的状态变量矩阵；

3)、根据每个所述待研究机组的发电机电磁部分的线性化状态空间获得每个待研究机组的发电机电磁部分的复频域模型：

G_G(s)X_G＝H_G(s)U_G+M_G(s)Δω_G

优选地，上述实施例的具体过程如下：

发电机本体仍采用转子d轴有1个阻尼绕组、q轴有2个阻尼绕组的6绕组模型。建立在自身dq坐标系下的发电机电磁部分建立在自身dq坐标系下的发电机电磁部分线性化状态空间为：

$s\left[\begin{array}{c}{X}_G\\ Z_G\end{array}\right]=A_G\left[\begin{array}{c}{X}_G\\ Z_G\end{array}\right]+B_G{U}_G+B_G^{\&prime;}\mathrm{\&Delta;}{\mathrm{\&omega;}}_G$

式中X_G＝[Δi_d Δi_q Δi_f Δi_D Δi_Q Δi_g]^T为发电机定/转子绕组电流矩阵，U_G＝[Δu_d Δu_q]^T为发电机的机端电压矩阵；Z_G为除X_G外发电机电磁部分的状态变量矩阵，含励磁控制器、PSS(Power System Stabilizer)等装置模型的状态量。消去Z_G，公式可化简为：

sX_G＝E_G(s)X_G+F_G(s)U_G+F′_G(s)Δω_G

从而可以得到发电机电磁部分的复频域模型为：

G_G(s)X_G＝H_G(s)U_G+M_G(s)Δω_G

其中，式中有G_G(s)＝sI-E_G(s)，H_G(s)＝F_G(s)，M_G(s)＝F′_G(s)。

另外，获得的网络部分复频域方程和发电机电磁部分复频域模型，即得到除发电机轴系机械部分外全系统的复频域模型。

在一些说明性实施例中，所述步骤4中计算出每个所述待研究机组的复转矩系数的过程，包括：

1)、利用每个待研究机组的复频域模型，得到该待研究机组的第一电磁转矩：

-ΔT_e＝K_M(s)Δδ_G

其中，K_M(s)为该机组的机械复转矩系数，s为拉普拉斯算子；

2)、利用所述网络部分的复频域方程，得到第二电磁转矩：

ΔT_e＝K_E(s)Δδ_G

其中，K_E(s)为该机组的电气复转矩系数；

3)、令s＝jω，则得到该待研究机组在ω频率下的复转矩系数：

K_M(ω)＝K_m(ω)+jωD_m(ω)

K_E(ω)＝K_e(ω)+jωD_e(ω)

其中，K_m(ω)和D_m(ω)分别为该待研究机组在ω频率下的机械弹性系数和机械阻尼系数，K_e(ω)和D_e(ω)分别为该待研究机组在ω频率下的电气弹性系数和电气阻尼系数。

优选地，上述实施例具体分析如下：

对于待研究发电机组而言，系统的状态空间可以分为机械子系统、电气子系统两个部分。机械子系统为待研究发电机组的轴系，其一般采用多质量块模型；电气子系统则包含待研究发电机组的电磁部分、外部网络和系统中的其他发电机组。上述两个子系统通过待研究发电机组的电磁转矩和发电机质量块位置角两个状态量连接为整个系统的状态空间。利用发电机组轴系的线性化状态方程，可以得到待研究机组线的性化的电磁转矩表达式为：

-ΔT_e＝K_M(s)Δδ_G

式中K_M(s)为该机组的机械复转矩系数，s为Laplace算子。

利用系统剩余部分的线性化状态方程，该机组线性化的电磁转矩又可表示为：

ΔT_e＝K_E(s)Δδ_G

式中K_E(s)为该机组的电气复转矩系数。

令s＝jω，即得到该机组在ω频率下的复转矩系数为：

K_M(ω)＝K_m(ω)+jωD_m(ω)

K_E(ω)＝K_e(ω)+jωD_e(ω)

式中，K_m(ω)和D_m(ω)分别为该机组在ω频率下的机械弹性系数和机械阻尼系数，K_e(ω)和D_e(ω)分别为该机组在ω频率下的电气弹性系数和电气阻尼系数。在整个次同步频率范围内对公式进行频率扫描，即可得到待研究机组的机械复转矩系数和电气复转矩系数。

在一些说明性实施例中，所述步骤5中所述分析出该待研究机组的存在次同步振荡风险的频率的过程，包括：

在轴系的自然扭振频率ω处，机械弹性系数和电气弹性系数满足如下关系：

K_m(ω)+K_e(ω)＝0

通过上述公式找到每个待研究机组的所有满足该公式的频率ω；

另外，在满足上述公式的基础上，再有机械阻尼系数和电气阻尼系数满足下式关系时，频率为ω的扭振扭振模态存在次同步振荡不稳定风险：

$D_m\left(\mathrm{\&omega;}\right)+D_e\left(\mathrm{\&omega;}\right)<0_{[K_m\left(\mathrm{\&omega;}\right)+K_e\left(\mathrm{\&omega;}\right)\&cong;0]}$

以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种近区密集火电多直流外送系统次同步振荡的分析方法，其特征在于，包含有以下步骤：

步骤1、获取近区密集火电多直流外送系统的静态参数；

步骤2、利用所述静态参数，筛选出存在次同步振荡风险的机组，作为待研究机组；

步骤3、构建网络部分的复频域方程和每个所述待研究机组的发电机电磁部分的复频域模型；

步骤4、根据所述复频域方程和所述复频域模型计算出每个所述待研究机组的复转矩系数；

步骤5、根据每个所述待研究机组的复转矩系数，分析出该待研究机组的存在次同步振荡风险的频率。

2.根据权利要求1所述的分析方法，其特征在于，所述步骤1中获取到的所述静态参数至少包括：近区密集火电多直流外送系统的线路参数、机组参数以及控制参数。

3.根据权利要求1所述的分析方法，其特征在于，所述步骤2中采用机组作用系数法筛选出存在次同步振荡风险的机组，具体包括：

1)、利用如下公式计算出每个机组的机组作用系数；

${\mathrm{UIF}}_i=\frac{S_{\mathrm{HVDC}}}{S_i}{(1-\frac{{\mathrm{SC}}_i}{{\mathrm{SC}}_{\mathrm{\&Sigma;}}})}^2$

其中，S_HVDC为高压直流输电的额定容量，S_i为发电机组的额定容量，SC_i为不计机组i影响时整流站交流母线的短路容量，SC_∑为计算机组i影响时整流站交流母线的短路容量；

2)、当UIF_i>0.1时，判定第i个机组为存在次同步振荡风险的机组。

4.根据权利要求1所述的分析方法，其特征在于，所述步骤3中构建网络部分的复频域方程的过程，具体包括：

构建系统等值网络及所述系统等值网络中的每个网络元件的复频域端口等效导纳矩阵；

利用网络元件的复频域端口等效导纳矩阵，采用节点分析法构建网络部分的复频域方程：

I＝Y_RLC(s)U

其中，I为支路的线性化电流矩阵，U为支路的线性化电压矩阵，Y_RLC(s)为支路的等效矩阵。

5.根据权利要求4所述的分析方法，其特征在于，构建每个网络元件的复频域端口等效导纳矩阵，具体包括：

建立在xy同步坐标系下的动态方程：

$\left[\begin{array}{c}\mathrm{\&Delta;}{u}_{\mathrm{ijx}}\\ \mathrm{\&Delta;}{u}_{\mathrm{ijy}}\end{array}\right]=\left[\begin{array}{cc}{x}_l& 0\\ 0& x_l\end{array}\right]s\left[\begin{array}{c}\mathrm{\&Delta;}{i}_{\mathrm{ijx}}\\ \mathrm{\&Delta;}{i}_{\mathrm{ijy}}\end{array}\right]+\left[\begin{array}{cc}r& -x_l\\ x_l& r\end{array}\right]\left[\begin{array}{c}\mathrm{\&Delta;}{i}_{\mathrm{ijx}}\\ \mathrm{\&Delta;}{i}_{\mathrm{ijy}}\end{array}\right]$

$\left[\begin{array}{c}\mathrm{\&Delta;}{i}_{\mathrm{Cx}}\\ \mathrm{\&Delta;}{i}_{\mathrm{Cy}}\end{array}\right]=\frac{1}{x_C}\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]s\left[\begin{array}{c}\mathrm{\&Delta;}{u}_{\mathrm{cx}}\\ \mathrm{\&Delta;}{u}_{\mathrm{Cy}}\end{array}\right]+\frac{1}{x_C}\left[\begin{array}{cc}0& -1\\ 1& 0\end{array}\right]\left[\begin{array}{c}\mathrm{\&Delta;}{u}_{\mathrm{Cx}}\\ \mathrm{\&Delta;}{u}_{\mathrm{Cy}}\end{array}\right]$

其中，I＝[Δi_x Δi_y]^T，U＝[Δu_x Δu_y]^T，Δ表示变量相对于其稳态值的微变量，x_l为网络元件电抗，x_c为网络元件电容，r为网络元件电阻，s为拉普拉斯算子。

6.根据权利要求5所述的分析方法，其特征在于，构建每个所述待研究机组的发电机电磁部分的复频域模型，具体包括：

1)、将所述建立在xy同步坐标系下的动态方程线性转换为在同步发电机dq坐标系下的动态方程；

其中，所述转换的关系如下：

$\left[\begin{array}{c}\mathrm{\&Delta;}{X}_x\\ \mathrm{\&Delta;}{X}_y\end{array}\right]=\left[\begin{array}{cc}\sin {\mathrm{\&delta;}}_{G0}& \cos {\mathrm{\&delta;}}_{G0}\\ -\cos {\mathrm{\&delta;}}_{G0}& \sin {\mathrm{\&delta;}}_{G0}\end{array}\right]\left[\begin{array}{c}\mathrm{\&Delta;}{X}_d\\ \mathrm{\&Delta;}{X}_q\end{array}\right]+\left[\begin{array}{cc}\cos {\mathrm{\&delta;}}_{G0}& -\sin {\mathrm{\&delta;}}_{G0}\\ \sin {\mathrm{\&delta;}}_{G0}& \cos {\mathrm{\&delta;}}_{G0}\end{array}\right]\left[\begin{array}{c}{X}_{d0}\\ X_{q0}\end{array}\right]\mathrm{\&Delta;}{\mathrm{\&delta;}}_G$

$\left[\begin{array}{c}\mathrm{\&Delta;}{X}_d\\ \mathrm{\&Delta;}{X}_q\end{array}\right]=\left[\begin{array}{cc}\sin {\mathrm{\&delta;}}_{G0}& -\cos {\mathrm{\&delta;}}_{G0}\\ \cos {\mathrm{\&delta;}}_{G0}& \sin {\mathrm{\&delta;}}_{G0}\end{array}\right]\left[\begin{array}{c}\mathrm{\&Delta;}{X}_x\\ \mathrm{\&Delta;}{X}_y\end{array}\right]+\left[\begin{array}{cc}\cos {\mathrm{\&delta;}}_{G0}& \sin {\mathrm{\&delta;}}_{G0}\\ -\sin {\mathrm{\&delta;}}_{G0}& \cos {\mathrm{\&delta;}}_{G0}\end{array}\right]\left[\begin{array}{c}{X}_{x0}\\ X_{y0}\end{array}\right]\mathrm{\&Delta;}{\mathrm{\&delta;}}_G$

其中，变量X可为发电机机端电压u或电流i，下标“0”表示变量的稳态值；

2)、根据在同步发电机dq坐标系下的动态方程构建每个所述待研究机组的发电机电磁部分的线性化状态空间：

$s\left[\begin{array}{c}{X}_G\\ Z_G\end{array}\right]=A_G\left[\begin{array}{c}{X}_G\\ Z_G\end{array}\right]+B_G{U}_G+B_G^{\&prime;}\mathrm{\&Delta;}{\mathrm{\&omega;}}_G$

其中，X_G＝[Δi_d Δi_q Δi_f Δi_D Δi_Q Δi_g]^T为发电机定子/转子绕组电流矩阵，U_G＝[Δu_d Δu_q]^T为发电机的机端电压矩阵；Z_G为除X_G外发电机电磁部分的状态变量矩阵；

3)、根据每个所述待研究机组的发电机电磁部分的线性化状态空间获得每个待研究机组的发电机电磁部分的复频域模型：

G_G(s)X_G＝H_G(s)U_G+M_G(s)Δω_G。

7.根据权利要求1所述的分析方法，其特征在于，所述步骤4中计算出每个所述待研究机组的复转矩系数的过程，包括：

1)、利用每个待研究机组的复频域模型，得到该待研究机组的第一电磁转矩：

-ΔT_e＝K_M(s)Δδ_G

其中，K_M(s)为该机组的机械复转矩系数，s为拉普拉斯算子；

2)、利用所述网络部分的复频域方程，得到第二电磁转矩：

ΔT_e＝K_E(s)Δδ_G

其中，K_E(s)为该机组的电气复转矩系数；

3)、令s＝jω，则得到该待研究机组在ω频率下的复转矩系数：

K_M(ω)＝K_m(ω)+jωD_m(ω)

K_E(ω)＝K_e(ω)+jωD_e(ω)

其中，K_m(ω)和D_m(ω)分别为该待研究机组在ω频率下的机械弹性系数和机械阻尼系数，K_e(ω)和D_e(ω)分别为该待研究机组在ω频率下的电气弹性系数和电气阻尼系数。

8.根据权利要求1所述的分析方法，其特征在于，所述步骤5中所述分析出该待研究机组的存在次同步振荡风险的频率的过程，包括：

找到满足公式K_m(ω)+K_e(ω)＝0的频率ω；

当存在时，判定该待研究机组在ω频率下存在次同步振荡风险。