小干扰稳定计算

1 / 15 第五章 小干扰稳定计算 一、实验目的 理解电力系统分析中小干扰稳定计算的相关概念，掌握PSASP小干扰稳定计算的过程。学会根据特性值判断系统的小干扰稳定性。复习PSASP潮流计算、暂态稳定计算。 二、预习要求

复习《电力系统分析》中有关小干扰稳定计算的内容，了解有关小干扰稳定计算的功能，掌握系统小干扰稳定性的判断方法。 三、实验内容

（一）PSASP小干扰稳定计算概述 电力系统小干扰稳定是指系统受到小干扰后，不发生自发振荡或非周期性失步，自动恢复到起始运行状态的能力。系统小干扰稳定性取决于系统的固有特性，与扰动的大小无关。 从理论上来说，电力系统的小干扰稳定性相当于一般动力学系统在李亚普诺夫意义下的渐近稳定性。当前，用于研究复杂电力系统小干扰稳定的方法主要是基于李雅普诺夫一次近似法的小干扰法。该方法的基本原理如下： 系统的状态方程为：XAX＝ 其中A为n×n维系数矩阵，称为该系统的状态矩阵。对于由状态方程描述的线性系统，其小干扰稳定性由状态矩阵的所有特征值决定。如果所有的特征值实部都为负，则系统在该运行点是稳定的；只要有一个实部为正的特征值，则系统在该运行点是不稳定的；如果状态矩阵A不具有正实部特征值但具有实部为零的特征值，则系统在该运行点处于临界稳定的情况。因此，分析系统在某运行点的小干扰稳定性问题，可以归结为求解状态矩阵A的全部特征值的问题。 PSASP小干扰稳定计算程序还提供了一些相应的分析手段，使之更加实用方便。其中包括： ➢ 特征值分布及其单线图上显示的模态图； ➢ 特征值和特征向量报表； 2 / 15

➢ 线性系统频域响应曲线，包括幅频特性、相频特性、乃奎斯特(Nyquist)曲线； ➢ 线性系统时域响应曲线。 PSASP小干扰稳定的过程如下图所示：

线性化时域频域响应

用基于稀疏性

的方法求解系统特征值QR法求特性值

系统状态矩阵A

系统增广矩阵J系统元件线性化网络线性化初值计算公用数据及模型库潮流结果

（二）数据准备 以WEPRI-7节点系统为例，其系统图如下：

PSASP程序中给出了WEPRI-7节点系统的基础数据，为方便起见，就用暂态稳定计算中参数导入的方法将基础数据库（Basic、G1-CTRL）、公用参数库、单线图、地理位置接线图等数据图形导入目标数据目录（C:\XGRJS\）。 3 / 15

（三）暂态稳定作业的定义 小干扰稳定计算中要用到系统的基础参数和调节器的参数，在暂态稳定计算中，这个参数都进行了定义。故在PSASP中，小干扰稳定计算作业是基于暂态稳定计算作业的，暂态稳定计算又是基于潮流的，而在潮流计算之前又要先定义方案。故为了进行小干扰的计算分析，现要用前面学到的方法定义一个方案、一个潮流作业和一个暂态稳定计算作业。 1. 方案定义： 序号 方案名 数据组 方案解释

1 有控制 Basic、G1-CTRL 将系统的直流线改为交流线的运行方式 2. 潮流作业定义： 作业号 方案名 描述

2 有控制 有控制的潮流计算 3. 暂态稳定计算作业定义： 作业号 潮流作业 描述

1 2 常规暂态稳定分析 （四）小干扰稳定计算作业的定义和执行 小干扰稳定计算作业是基于暂态稳定计算作业的。实际上，小干扰稳定与暂稳作业中的故障、扰动及输出等信息无关，与暂稳作业基于的潮流作业所定义的初始稳态运行点有关，与暂稳作业的发电机及其调节系统、负荷、直流输电、UD等元件模型有关。小干扰稳定计算作业的定义也有文本和图形两种方式，其构成如下图： 4 / 15

暂态稳定计算作业

计算控制信息特征值搜索范围小干扰稳定计算作业

算法控制信息算法和功能选择

时域/频域输入输出信息 在文本环境窗口中点击“计算”下拉菜单中的“小干扰稳定”或在图形支持的运行模式窗口中点击“作业”下拉菜单中的“小干扰稳定”，弹出暂态稳定作业定义窗口：

PSASP暂态稳定计算中控制信息包括：算法及功能、计算信息、算法控制信息、输入/输出信息。各栏介绍如下： 作业定义栏 作业号：可选择已有的或键入新的作业号。 暂稳作业：从已定义好的暂稳计算作业中，选择其一，即确定其所基于的暂5 / 15

态稳定计算作业。 功能选择栏 “算法及功能”栏的组合框中列出了可选的特征值算法和其他计算功能，包括： ➢ QR法； ➢ 逆迭代转Rayleigh商迭代法； ➢ 同时迭代法； ➢ 线性化时域响应计算； ➢ 线性化频域响应计算。 计算信息栏 “计算信息”栏中的内容仅用于稀疏算法，包括逆迭代/Rayleigh商迭代法、同时迭代法和线性化时域/频域响应。其中：“特征值搜索范围”栏给出所要计算的特性值范围和频域响应计算的频域范围；“算法控制信息”栏中定义用于控制计算迭代的信息以及时域/频域响应的计算总时间和计算步长，各项含义如下： 特征值搜索范围栏 起始点实部（SSR）、起始点虚部（SSI）、终止点实部（ESR）、终止点虚部（ESI）：对于逆迭代/Rayleigh商迭代法和同时迭代法，为迭代初始点选取的范围；对于线性化频域响应计算，仅虚部SSI,ESI起作用，为指定的频域范围，即在频域[SSI,ESI]内进行频域响应计算。单位为弧度/秒(rad/s)。 迭代初始值个数（NSHFT）：适用于逆迭代/Rayleigh商迭代法和同时迭代法，指定迭代初始点的个数，这些迭代初始点等间距地分布于复平面中(SSR,SSI) - (ESR,ESI)线段上。 算法控制信息栏 迭代次数上限：指定迭代次数上限，若迭代次数超过该值则终止计算。 允许误差（EPS）：指定迭代算法的迭代收敛精度，实际迭代的允许误差为： ||10EPS

如若EPS=6或者EPS=-6，则 ε=10-6。 每次求解特征值个数（NEIG）：指定由每个迭代初始点开始求出的特征值个数，因此一次计算所求出的特征值个数为NSHFT * NEIG。 6 / 15

计算总时间：适用于线性化时域响应计算，指定线性化时域响应计算的总时间长度，单位为秒(s)。 计算步长：适用于线性化时域和频域响应计算。对于线性化时域响应计算，指定时域响应计算的积分步长，单位为秒(s)；对于线性化频域响应计算，指定频域响应计算的步长，单位为弧度/秒(rad/s)。 输入/输出信息栏 对于线性化时域和频域响应计算，还需指定输入和输出信息。具体方法为：先点击“线性化”按钮，在指定的初始运行点进行全系统微分方程和代数方程的线性化，之后再分别点击“输入选择”和“输出选择”按钮，指定时域或频域响应计算所需的输入和输出信息。 点击“输入选择”按钮，弹出输入选择窗口，在该窗口中依次指定输入变量的位置(母线或支路)、模型和变量，之后点击“确定”按钮，即可完成输入变量的定义。对于时域和频域响应计算，一个作业仅允许选择1个输入变量。 点击“输出选择”按钮，弹出输出选择窗口，在该窗口中依次指定输出变量的位置(母线或支路)、模型和变量，之后点击“保存”按钮，完成输出变量的定义。点击“退出”按钮，返回小干扰稳定计算信息窗口。对于时域响应计算，一个作业最多允许选择8个输出变量；对于频域响应计算，一个作业仅允许选择1个输出变量。 本实验定义5个小干扰稳定计算作业： 序号 暂稳作业号 描述 1 1 QR法计算系统的全部特性值 2 1 逆迭代转Rayleigh商迭代法计算系统的部分特征值 3 1 同时迭代法计算系统的部分特征值 4 1 线性化时域响应计算 5 1 线性化频域响应计算 在文本和图形环境下均可以定义不同的潮流计算作业，每一作业需给定一个作业号。现给出各作业的定义方法。 7 / 15

作业1: 在小干扰稳定计算信息窗口中定义作业1，如下图所示：

在该窗口中，给出小干扰稳定作业号及其所基于的暂稳作业号，在“算法及功能”下拉框中选择“QR法”即完成了采用QR法计算系统全部特征值的小干扰稳定作业定义，这时不需对计算信息栏进行定义。 作业2： 在小干扰稳定计算信息窗口中定义作业2，如下图所示： 8 / 15

在“算法及功能”下拉框中选择“逆迭代转Rayleigh商迭代法”，指定特征值的搜索范围，给出迭代次数上限，允许误差和每次求解特征值个数等信息，即完成了采用逆迭代转Rayleigh商迭代法计算系统部分特征值的小干扰稳定作业定义。 作业3： 在小干扰稳定计算信息窗口中定义作业3，如下图所示：

在该窗口中，给出小干扰稳定作业号及其所基于的暂稳作业号，在“算法及功能”下拉框中选择“同时迭代法”，指定特征值的搜索范围，给出迭代次数上限，允许误差和每次求解特征值个数等信息，即完成了采用同时迭代法计算系统部分特征值的小干扰稳定作业定义。 作业4： 在小干扰稳定计算信息窗口中定义作业4，如下图所示： 9 / 15

在“算法及功能”下拉框中选择“线性化时域响应计算”。必须首先点击“线性化”按钮，完成系统的线性化计算。 点击“输入选择”按钮，弹出“输入选择”窗口：

在该窗口中可以依次指定输入信息的位置、模型和变量。以上窗口所定义的输入信息为发电机G2的调压器输入Vs。 点击“输出选择”按钮，弹出“输出选择”窗口：