基于部分左特征结构配置的电力系统优化阻尼控制方法与流程

本发明涉及电力系统阻尼控制领域，特别是一种基于部分左特征结构配置的电力系统优化阻尼控制方法。

背景技术：

随着西电东输工程的建设和全国电力联网的实施，我国将进入跨区域和远距离传输巨大功率的超高压、交直流混合输电的时代。以此带来的电力系统阻尼作用的减小，从而在一些扰动情况下容易引发电力系统低频振荡的问题，已成为限制电力系统功率传输能力和影响大型互联电网安全稳定运行的重要因素。

目前，在电力系统阻尼控制中广泛采用的是基于留数和相位补偿方法配置的传统pss，针对单个振荡模态设计控制器参数，配置闭环系统的特征值向左平面移动，增加系统阻尼抑制低频振荡。在实际电力系统中往往同时存在多个需要抑制的振荡模态，这也使得多pss的协调控制成为一个难题。另外，从多输入多输出系统的控制角度来看，基于特征值配置的输出反馈控制中仍然存在着多余的设计自由度，可用于电力系统的优化阻尼控制。为解决以上问题，有效地抑制现代大型互联电网的低频振荡，设计一种有效可靠的电力系统优化阻尼控制器十分重要。

技术实现要素：

本发明的目的在于提供一种基于部分左特征结构配置的电力系统优化阻尼控制方法，以克服现有技术中存在的缺陷。

为实现上述目的，本发明的技术方案是：一种基于部分左特征结构配置的电力系统优化阻尼控制方法，包括如下步骤：

步骤s1：给定主导机电振荡模式；

步骤s2：采用参与因子法和留数法选择控制器安装位置和广域反馈信号；

步骤s3：引入动态补偿型pss；

步骤s4：建立最小化系统输出振荡幅度的子目标函数f1(φ)；

步骤s5：将主导特征值移动到正确位置，并建立子目标函数f2(φ)；

步骤s6：通过限制闭环系统未配置特征分布，建立子目标函数f3(φ)和f4(φ)；

步骤s7：限制闭环系统控制能量，建立最小化控制代价的子目标函数f5(φ)；

步骤s8：通过合并，获取单目标函数f(φ)；

步骤s9：随机初始种群并根据f(φ)计算适应度函数；

步骤s10：采用群搜索优化算法求解控制器状态空间参数；

步骤s11：判断是否满足终止条件，若满足则输出阻尼控制器状态空间参数，否则转到所述步骤s9继续进行参数计算。

在本发明一实施例，在所述步骤s2中，还包括如下步骤：

步骤s21：分别以机组功角δ和转子角频率ω作为状态变量，计算各机组对所述主导机电振荡模式的参与因子，综合考虑，筛选出在所述主导机电振荡模式下，有最大参与因子的发电机组作为阻尼控制器的最佳安装位置；

步骤s22：引入机组转子角频率ω作为反馈信号构成广域控制回路，计算主导振荡模式中ω对励磁参考电压的留数矩阵rⁱ＝cφiψi^tb；其中，ψ和φ分别为电力系统全状态空间模型特征矩阵的左右特征向量；所述留数矩阵中各机组对应所述步骤s21最佳安装位置的最大元素，即为各机组转子角频率对主导振荡模式的留数；

步骤s23：在步骤s22计算的主导振荡模式留数中，筛选出对应于所给振荡模式留数最大的机组作为本地信号和优选的广域反馈信号；

步骤s24：将电力系统模型在工作点附近通过线性时不变状态方程来描述：

其中，x(t)∈r^n×1、y(t)∈r^r×1、u(t)∈r^m×1分别为系统状态向量、输出向量和控制向量，a∈r^n×n、b∈r^n×m、c∈r^r×n分别为状态矩阵、控制矩阵和输出矩阵，n、r、m为对应矩阵的维度。

在本发明一实施例，在所述步骤s3中，还包括如下步骤：

步骤s31：引入动态补偿器型pss：其中，z(t)∈r^a为控制状态向量，且该动态补偿器型pss以所述步骤s23中本地反馈信号和优选的广域量测信号作为输入，所述参数矩阵d、e、f和g通过所述步骤s10中状态空间参数获取；

步骤s32：引入动态补偿器后，将闭环电力系统表示成等值增广系统：其中，所述等值增广系统的输出与多变量输入呈线性比例关系；所述等值增广系统为一个比例输出反馈控制系统，其状态变量输入和输出变量数分别为：其中，a为动态补偿器的状态向量维度；所述

在本发明一实施例，在所述步骤s4中，通过开环振荡模态参与因子p近似系统初始状态，建立所述最小化系统输出振荡幅度的子目标函数：其中，ω为系统左特征向量，m为左特征向量个数。

在本发明一实施例，在所述步骤s5中，利用待配置的特征根，对控制器参数的灵敏度估计待配置特征值的迁移量和并计算待配置闭环特征值的估计值将λi和差值的归一化作为优化目标子函数：保证待配置特征根移动到正确位置。

在本发明一实施例，在所述步骤s6中，通过限制闭环系统未配置特征分布，建立子目标函数：

其中，为阻尼比，为实部值，ζmin为未配置部分闭环系统特征值的最小阻尼比，αmin为实部最小值。

在本发明一实施例，在所述步骤s7中，限制闭环系统控制能量，将比例反馈控制增益k中最大元素kmax限制在区间[-20,20]范围内，建立最小化控制代价的子目标函数：

在本发明一实施例，在所述步骤s8中，引入权值系数βi，将所述步骤s4、步骤s5、步骤s6和步骤s7的控制器优化多个子目标函数合并为单目标函数，也即为群搜索优化算法中的适应度函数，：

在本发明一实施例，在所述步骤s9中，随机初始种群并根据f(φ)计算适应度函数；在搜索空间rⁿ中，第i个成员在第k次迭代中的坐标位置为搜索角度为

相较于现有技术，本发明具有以下有益效果：本发明提出的一种基于部分左特征结构配置的电力系统优化阻尼控制方法，该方法将多余的设计自由度用于配置特征向量，引入动态补偿器型pss增加系统输入，满足部分特征结构配置方法中控制自由度要求，并实现控制器的参数化。优选远端广域量测反馈信号，提高阻尼控制效果。最后，建立加权多目标函数，将控制器参数优化归纳为一个无约束的非线性规划问题，并利用群搜索优化算法对控制器参数求解。依据该方法设计的阻尼控制器能有效抑制低频振荡，较基于留数法设计的传统pss有更好阻尼控制效果和鲁棒性。通过将多余的设计自由度配置到特征向量上，能够同时对电力系统的多个主导振荡模态进行有效的抑制。通过引入动态补偿器型pss，满足部分左特征结构配置方法对系统输出反馈变量的要求，同时增加闭环电力系统的阶次，可以使电力系统特征结构在更高的维度中优化配置。

附图说明

图1为本发明中基于部分左特征结构配置的电力系统优化阻尼控制方法的流程图。

图2为本发明一实施例中的ieee-39节点电力系统模型单线图。

图3为本发明一实施例中的三种联络线停运条件下机组对主导振荡模态的留数计算结果图。

图4为本发明一实施例中的四种运行条件下闭环控制效果图。

图5为本发明一实施例中的短路故障下发电机g1转子频率响应曲线。

图6为本发明一实施例中的短路故障下发电机g9转子频率响应曲线。

图7为本发明一实施例中的短路故障下发电机组g1和g9的相对功角动态响应曲线。

具体实施方式

下面结合附图，对本发明的技术方案进行具体说明。

本发明一种基于部分左特征结构配置的电力系统优化阻尼控制方法，如图1所示，

包括以下步骤：

步骤s1：给定主导机电振荡模式，在如图2所示的ieee-39节点仿真系统上设计阻尼控制器，分析本发明对电力系统低频振荡的抑制效果。

步骤s2：分别以机组功角δ和转子角频率ω作为状态变量，计算各机组对主导振荡模式的参与因子，筛选出对应于所给主导振荡模式中较大参与因子的发电机组，作为阻尼控制器的最佳安装位置。

步骤s3：引入机组转子角频率ω作为反馈信号构成广域控制回路，计算主导振荡模式中ω对励磁参考电压的留数矩阵rⁱ＝cφiψi^tb；其中，ψ和φ分别为电力系统全状态空间模型特征矩阵的左右特征向量；留数矩阵中各机组对应步骤s2最佳安装位置的最大元素，即为各机组转子角频率对主导振荡模式的留数，选择运行方式1为基本工作模式，联络线4-14、17-18、25-26停运分别为运行方式2、3、4，在三种联络线停运条件下机组对主导振荡模态的留数计算结果示于附图3中。

步骤s4：在步骤s3计算的主导振荡模式留数中，筛选出对应于所给振荡模式留数最大的机组作为本地信号和优选的广域反馈信号；采用该输入信号，得到的闭环系统特征根分布示于附图4中。

步骤s5：电力系统模型在工作点附近用线性时不变状态方程来描述：其中，x(t)∈r^n×1、y(t)∈r^r×1、u(t)∈r^m×1分别为系统状态向量、输出向量和控制向量，a∈r^n×n、b∈r^n×m、c∈r^r×n分别为状态矩阵、控制矩阵和输出矩阵，n、r、m为对应矩阵的维度。

步骤s6：引入动态补偿器型pss：其中，z(t)∈r^a为控制状态向量。

步骤s7：引入动态补偿器后，将闭环电力系统表示成等值增广系统：其中，等值增广系统的输出与多变量输入呈线性比例关系；所述等值增广系统为一个比例输出反馈控制系统，其状态变量输入和输出变量数分别为：其中，a为动态补偿器的状态向量维度。

步骤s8：用开环振荡模态参与因子p近似系统初始状态，建立最小化系统输出振荡幅度的子目标函数：

步骤s9：利用待配置的特征根对控制器参数的灵敏度估计待配置特征值的迁移量和并计算待配置闭环特征值的估计值将λi和差值的归一化作为优化目标子函数：保证待配置特征根移动到正确位置。

步骤s10：通过限制闭环系统未配置特征分布，建立子目标函数：

其中，

步骤s11：限制闭环系统控制能量，将比例反馈控制增益k中最大元素kmax限制在区间[-20,20]范围内，建立最小化控制代价的子目标函数：

步骤s12：引入权值系数βi，将步骤s8、步骤s9、步骤s10和步骤s11的控制器优化多个子目标函数合并为单目标函数，即为群搜索优化算法中的适应度函数：

步骤s13：随机初始种群并根据计算适应度函数；在搜索空间rⁿ中，第i个成员在第k次迭代中的坐标位置为搜索角度为

步骤s14：将步骤s12中的权值系数βi设为相同的0.2，并设置适应度阈值为0.001，采用群搜索优化算法求解控制器状态空间参数。

步骤s15：判断是否满足终止条件，若满足则输出阻尼控制器状态空间参数，否则转到步骤s13继续进行参数计算；

进一步的，在本实施例中，步骤s6中的动态补偿器型pss采用步骤s4中的本地反馈信号和优选的广域量测信号作为输入。

进一步的，在本实施例中，步骤14中的状态空间参数构成等值增广矩阵的反馈控制器k，即为步骤s6中的参数矩阵d、e、f和g；

进一步的，在本实施例中，如图5、图6和图7所示，投入依据本实施例提供的基于部分左特征结构配置方法设计的动态补偿型控制器后，机组g1、g9母线频率和其相对功角的振荡幅度明显减小，并且振荡信号能够更快衰减和到达稳态。

以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作用未超出本发明技术方案的范围时，均属于本发明的保护范围。