抑制次同步振荡的SSDC和SEDC协调控制优化方法

摘要

一种抑制次同步振荡的SSDC和SEDC协调控制优化方法，包括以下步骤：①获取电力系统元件参数；②建立SSDC和SEDC的数学模型；③建立系统元件的数学模型，并形成包含SSDC、SEDC的系统状态方程；④由式(3)求出系统的特征值

说明书

技术领域：

    本发明涉及电力系统稳定与控制技术，特别涉及一种抑制次同步振荡的SSDC和SEDC协调控制优化方法。

背景技术：

近年来，高压直流输电(HVDC)以其先进的技术和明显的经济性等优势在电力系统输电方面得到了快速发展，对我国电力的远距离输送、区域电网互联起到了重要作用。但也由于其直流换流器控制的固有反馈控制特性，HVDC容易引起与整流侧连接的汽轮发电机组的次同步振荡(SSO)问题。次同步振荡是由HVDC引起的一种电气-机械共振现象，轻微的次同步振荡可能会造成发电机轴系的疲劳损伤，严重时则可能导致发电机轴系出现裂纹甚至断裂，造成巨大的经济损失。

直流附加励磁阻尼控制器（SSDC）和附加励磁阻尼控制器（SEDC）是目前用于抑制SSO的两种较为有效的方法。SSDC属于直流系统的附加设备，通过引入机组的次同步分量，经移相放大后将输出叠加到整流侧的定电流控制器中；SEDC为附加在发电机励磁控制器上的阻尼控制器，其采取汽轮机高压缸转速作为输入信号，其输出叠加到励磁电压上，形成一个与次同步电流反向的电流分量，从而达到抑制振荡的效果。SSDC和SEDC的比例放大倍数直接关系着其向系统提供阻尼的大小，比例放大倍数越大，SSDC和SEDC提供的阻尼越大，阻尼控制器抑制次同步振荡的效果越好。但SSDC的放大倍数过大会影响直流系统的动态稳定性，SEDC则容易受到励磁控制器容量的影响，其放大倍数过大易达到励磁控制器的上限。我国目前已投入SSDC的参数大多采用制造厂商给出的参数，而SEDC的参数也仅基于单台机组的参数来整定，二者没有相互协调，亦没有从全系统的角度出来进行设计，目前次同步阻尼控制器的抑制效果尚未得到充分发挥。

原对偶内点法寻优速度快、收敛性能好，在处理大规模、非线性的寻优问题时具有较大的优势，计算时间对问题规模不敏感，不会随着问题规模的增大而增大，且整个计算过程均在可行域的内部进行求解。若能将原对偶内点法应用到附加励磁阻尼和直流附加励磁阻尼协调控制参数的优化设计上，可充分考虑到全系统的动态性能，在控制器输出不越限的情况下，使控制器达到最好的振荡抑制效果。

发明内容：

鉴于此，本发明提供一种抑制次同步振荡的SSDC和SEDC协调控制优化方法。

一种抑制次同步振荡的SSDC和SEDC协调控制优化方法，包括以下步骤：

获取电力系统元件参数，包括：发电机轴系转动惯量和弹性系数K_i, _i+1；励磁系统放大倍数和时间常数；发电机各绕组自感、互感以及各绕组等效电阻；交流线路阻抗；直流定电流控制器PI参数、逆变侧定熄弧角PI参数、换流变压器变比以及等效换向电抗、整流器和逆变器触发角、直流功率和直流电压；

建立SSDC和SEDC的数学模型；其中，SSDC为宽带通单通道结构，采用发电机质块转速差为输入信号；SSDC由信号输入、4阶带通滤波器、相位补偿模块、比例放大模块组成；SEDC为宽带通多通道结构，其通道数与振荡模态数相同，用发电机质块转速差为输入信号；SEDC输入信号经4阶带通滤波器滤除直流和低频分量后，分模态由模态带通滤波器、带阻滤波器、相位补偿模块以及比例放大模块组成；

SSDC的数学模型为

                                        (1)

其中，为SSDC的比例放大系数，为SSDC相位补偿模块传递函数的时域形式，为SSDC4阶带通滤波器传递函数的时域形式；

若发电机有三个振荡模态，SEDC的数学模型为

    (2)

其中，为SEDC各振荡模态的比例放大系数，为各振荡模态相位补偿模块传递函数的时域形式，为4阶带通滤波器传递函数的时域形式，为各振荡模态带通滤波器传递函数的时域形式，各振荡模态带阻滤波器传递函数的时域形式；

建立系统元件的数学模型，并形成包含SSDC、SEDC的系统状态方程；

                                              （3）         为状态变量，为与控制器比例放大系数相关的系数矩阵，

，其中表示发电机轴系各质量块的角速度；

由式(3)求出系统的特征值；特征值可表征系统的稳定性，将系统特征值的实部作为控制目标，若系统特征值实部的最大值离虚轴越远，则代表系统越稳定，故控制器参数的优化问题可转化为如下目标函数来求解：

                                         （4）               其中表示所有可能的控制器参数的集合，表示所有可能的运行条件的集合。

优选的，令，引入系统等式约束条件和不等式约束条件，控制器参数通过下列方程组求解

                                           （5）                                              （6）        

其中，、分别为SEDC各模态放大倍数的下限和上限，、分别取-0.03和+0.03；、分别为SSDC放大倍数的下限和上限，、分别取-0.1和+0.1；为系统各状态相量的向量形式。

优选的，式(5) 利用原对偶内点法求解，包括以下步骤：

形成原对偶内点法的修正方程；

引入松弛变量，构造罚函数，将不等式约束转化为等式约束有：

                         （7）  

式中，为罚因子，、为松弛变量；

利用拉格朗日函数法进行求解，可得到如下的拉格朗日函数

           （8）     式中，、和为拉格朗日乘子；

拉格朗日函数对公式（8）所有变量及乘子的偏导数为0，具体化后得到非线性方程组：

              （9）         式中，和分别为等式约束和不等式约束的雅克比矩阵的转置；

    由上式可得，定义对偶间隙；

将公式（9）非线性方程组线性化并 

写成矩阵形式，得到原对偶内点法的修正方程为：

    （10）     其中，

；

确定迭代次数和迭代步长，并求解方程（10）得到最优解的一个新的近似解为：

                                     （11）   

式中：k为迭代次数；和为步长，其计算公式为：

                （12）                   （13）  

上式的取值保证迭代点严格满足l>0，u>0，z>0，w>0；

判断迭代是否收敛；利用对偶间隙作为算法的收敛判据，若则算法收敛，计算结束；若不收敛，则利用重新计算惩罚因子，为设定参数，为取0.1；

通过迭代得到向量中各变量值，求得到一组最优的阻尼控制器协调控制参数。

相对于传统的SSDC、SEDC参数设计方法，本发明所提供的方法考虑了SSDC和SEDC的协调控制问题，在控制器输出不超过系统允许的范围内，以实现系统特征值最大实部最小为控制目标，得到一组最优的控制器比例放大系数值，该最优参数可使系统整体阻尼得到增强，从而达到最好的抑制效果。

附图说明：

图1为单机经直流输电接入系统模型示意图。

图2为SSDC结构图示意图。

图3为SEDC结构图示意图。

具体实施方式：

如图1所示的单机经直流输电接入系统，对本发明方法进行阐述说明。发电机、励磁的模型与参数和IEEE第一标准模型相同，轴系采用四质量块模型，直流输电采用CIGRE的直流输电标准测试模型，直流额定功率1000MW，额定电压500kV，整流侧采用定电流控制，逆变侧采用定熄弧角控制。

获取电力系统元件参数，包括：发电机轴系转动惯量和弹性系数K_i, _i+1；励磁系统放大倍数和时间常数；发电机各绕组自感、互感以及各绕组等效电阻；交流线路阻抗；直流定电流控制器PI参数、逆变侧定熄弧角PI参数、换流变压器变比以及等效换向电抗、整流器和逆变器触发角、直流功率和直流电压；

建立SSDC和SEDC的数学模型；其中，SSDC为宽带通单通道结构，如图2所示，SSDC采用发电机质块转速差为输入信号；SSDC由信号输入、4阶带通滤波器、相位补偿模块、比例放大模块组成；SEDC为宽带通多通道结构，如图3所示，SEDC通道数与振荡模态数相同，用发电机质块转速差为输入信号；SEDC输入信号经4阶带通滤波器滤除直流和低频分量后，分模态由模态带通滤波器、带阻滤波器、相位补偿模块以及比例放大模块组成；

SSDC的数学模型为

                                      (1)

其中，为SSDC的比例放大系数，为SSDC相位补偿模块传递函数的时域形式，为SSDC4阶带通滤波器传递函数的时域形式；

若发电机有三个振荡模态，SEDC的数学模型为

     (2)

其中，为SEDC各振荡模态的比例放大系数，为各振荡模态相位补偿模块传递函数的时域形式，为4阶带通滤波器传递函数的时域形式，为各振荡模态带通滤波器传递函数的时域形式，各振荡模态带阻滤波器传递函数的时域形式；

建立系统元件的数学模型，并形成包含SSDC、SEDC的系统状态方程；

直流系统模型：

SSDC的输出信号叠加到整流侧定电流控制器的电流参考环节上，故定电流控制环节的数学模型重写为：

                                           (3)

其中为整流侧触发角，、分别为整流侧PI环节的比例系数和时间常数，为直流电流偏差值。

线性化后得到直流系统数学模型的简化矩阵形式为

                                            (4)

其中状态变量；为直流电流，为逆变器熄弧角。

交流系统模型：

将整流侧交流母线电压、逆变侧交流母线电压和等值系统阻抗电流的坐标分量取为交流系统状态变量，所以交流系统为六阶模型。

其线性化方程为

                          (5)

消去中间变量，得到交流系统线性化方程的简化矩阵形式为

                            (6)

发电机系统模型：

发电机的线性化方程如下

                               (7)

式中，为发电机转速，，为定子、绕组、励磁绕组和转子等效阻尼绕组、、的磁链，选为状态变量；为对应的绕组电压，其中；为对应的绕组电阻；为对应的绕组电流；系数矩阵；。

升压变压器电压降线性化方程为

         (8)

式(8)中、为变压器的电阻和电感；、为、的初值；、为换流变母线电压的dq轴分量。

加入SEDC后，发电机磁链方程中可表示为

                 (9)

为励磁系统输出的控制状态变量，为SEDC分模态输出控制变量。

将发电机磁链方程中、消去，得到发电机磁链线性化方程。

励磁系统模型：

励磁系统的简化一阶模型的线性化方程为

                            (10)

通过前面的变压器方程式(8)消去，用磁链表示，可将励磁系统线性化方程写成简化矩阵形式

                             (11)

最后，整理并将SSDC、SEDC线性化后的数学模型代入，得到整个系统包含SSDC、SEDC的线性化状态方程：

                                                   (12)

为状态变量，为与控制器比例放大系数相关的系数矩阵，

，其中表示发电机轴系各质量块的角速度，前述。

由式(12)求出系统的特征值。特征值可表征系统的稳定性，将系统特征值的实部作为控制目标，若系统特征值实部的最大值离虚轴越远，则代表系统越稳定，则控制器参数的优化问题可转化为如下目标函数来求解：

                                             (13)

其中表示所有可能的控制器参数的集合，表示所有可能的运行条件的集合。

优选的，令，利用式(3)~(11)的系统方程可构成寻优过程的等式约束条件，以励磁系统和直流系统允许加入的模态信号幅值范围构成系统的不等式约束条件，则控制器参数通过下列方程组求解：

                                              (14)

                                         (15)

其中，、分别为SEDC各模态放大倍数的下限和上限，、分别取-0.03和+0.03；、分别为SSDC放大倍数的下限和上限，、分别取-0.1和+0.1；为系统各状态相量的向量形式。

 

优选的，利用原对偶内点法求解式（14），即可得到一组最优的控制器参数，步骤如下：

形成原对偶内点法的修正方程；

首先引入松弛变量，构造罚函数，将不等式约束转化为等式约束有：

                            (16)

式中，为罚因子，、为松弛变量。

利用拉格朗日函数法进行求解，可得到如下的拉格朗日函数

                (17)

式中，、和为拉格朗日乘子。

该问题存在极小值的必要条件是拉格朗日函数对所有变量及乘子的偏导数为0，即满足KKT条件，具体化后得到非线性方程组：

                    (18)

式中，和分别为等式约束和不等式约束的雅克比矩阵的转置。

    由上式可得，定义对偶间隙。

将上述非线性方程组线性化并 

写成矩阵形式，得到原对偶内点法的修正方程为：

         (19)

其中，

；

确定迭代次数和迭代步长，并求解方程（19）得到最优解的一个新的近似解为：

                                          (20)

式中：k为迭代次数；和为步长，其计算公式为：

                     (21)

                     (22)

上式的取值保证迭代点严格满足l>0，u>0，z>0，w>0。

判断迭代是否收敛。利用对偶间隙作为算法的收敛判据，若则算法收敛，计算结束；若不收敛，则利用重新计算惩罚因子，为设定参数，取0.1。

通过迭代得到向量中各变量值，即可得到一组最优的阻尼控制器协调控制参数的最优解。