一种基于改进粒子群算法的超低频振荡抑制方法与流程

本发明涉及电力系统稳定控制技术，具体涉及一种基于改进粒子群算法的超低频振荡抑制方法。

背景技术：

随着鲁西背靠背以及渝鄂背靠背工程的相继推进，异步联网后送端电网主要特性与风险由区域电网间的暂态功角稳定问题转换为区内电网功率平衡导致的频率稳定问题。近年来，随着发电机组调速系统的技术进步，调速系统的性能得到大幅提升，响应时间越来越快，以达到电网频率出现偏差时能够迅速地调出一次调频调节量将频率拉回正常值目的。然而由于水轮机导叶频繁开启或关闭，水轮机水锤效应的逆向调节增强，将显著地对系统的动态稳定和响应特性带来不利影响，成为引发系统超低频振荡的关键因素。目前针对异步联网后送端电网出现的超低频振荡现象，相关抑制措施的研究还较少，主要研究集中在对水轮机组调速系统参数进行优化，来达到降低水轮机负阻尼的目的，但因调速系统的阻尼特性、调节性能受到液压系统、水轮机固定参数以及系统一次调频能力的限制，对调速器参数进行优化并不能理想的解决系统频率稳定性存在的问题。

技术实现要素：

本发明解决了异步联网后区内超低频振荡的问题，提供一种基于改进粒子群算法的超低频振荡抑制方法，其应用时可以通过改进粒子群的算法，提高算法的收敛速度与精度，以较少的直流功率变化实现对超低频振荡的有效抑制。

本发明通过以下技术方案实现:

一种基于改进粒子群算法的超低频振荡抑制方法，所述方法包括的步骤为：

s1：对多直流送出系统的每条直流分别施加相同的功率扰动，直至激发系统发生超低频振荡，并记录各条直流扰动下的超低频振荡模态幅值；

s2：比较各条直流扰动下的超低频振荡模态幅值，并根据超低频振荡模态幅值选择直流频率限制控制器(frequencylimitcontroller，fc)的控制敏感点；

s4：将所述fc的控制敏感点作为基础，通过tls-esprit算法辨识系统传递函数，建立系统频率偏差对时间的积分最小的粒子群算法优化模型；

s5：对所述粒子群设定最大迭代次数、初始粒子数目、加速常数、惯性权重取值范围、非线性调节因子、维度及解空间范围，并将当前迭代次数设定为1；

s6：采用粒子群算法对粒子的速度和位置进行更新；

s7：记录每个微粒经历过的最好位置，所述微粒经历过的最好位置为个体极值，并将个体极值与微粒的适应值相比较，并根据每个微粒个体极值的位置与适应值的位置，确定出当前最好的个体极值；

s8：将全局中微粒经历过的最好位置定为全局极值，并将每个微粒的适应值与全局极值作比较，并根据每个微粒全局极值的位置与适应值的位置，确定出当前最好的全局极值；

s9：使惯性权重的非线性减小，并更新惯性权重取值；

s10：判断当前迭代次数，如果当前迭代次数等于最大迭代次数，则输出fc参数优化结果，反之则将当前迭代次数加1，并返回步骤s6。目前针对异步联网后送端电网出现的超低频振荡现象，相关抑制措施的研究还较少，主要研究集中在对水轮机组调速系统参数进行优化，来达到降低水轮机负阻尼的目的，但因调速系统的阻尼特性、调节性能受到液压系统、水轮机固定参数以及系统一次调频能力的限制，对调速器参数进行优化并不能理想的解决系统频率稳定性存在的问题，所以本发明提供一种基于改进粒子群算法的超低频振荡抑制方法，通过对多直流送出系统的每条直流施加相同功率的扰动，并通过激发出的超低频震荡模态幅值选择直流频率限制控制器(fc)的控制敏感点，并以fc的控制敏感点为基础，利用tls-esprit算法辨识系统传递函数，将系统频率偏差对时间的积分最小作为优化目标，建立粒子群算法优化模型，设定粒子群的最大迭代次数、初始粒子数目、加速常数、惯性权重取值范围、非线性调节因子、维度及解空间范围，并将当前迭代次数设定为1。采用粒子群算法对粒子的速度和位置进行更新后，然后通过微粒经历过的最好位置，将每个微粒的适应值与其个体极值、全局极值作比较，确定出当前的最好个体极值与当前的最好全局极值，将惯性权重的非线性减小，更新惯性权重的取值，直到当前迭代次数等于或大于最大迭代次数的时候，输出fc参数优化结果，当前迭代次数小于最大迭代次数的时候，返回步骤s6。通过改进粒子群算法的方式来对超低频振荡进行抑制，以fc作为抑制超低频的主要手段，在对fc参数进行优化时，有效的提高了收敛精度与速度，可以更加直接、简单、有效的对多直流送出系统的超低频振荡进行抑制，保证系统的稳定运行。

进一步的，一种基于改进粒子群算法的超低频振荡抑制方法，步骤s9中的惯性权重的非线性减小采用的方法为：

式中，w代表惯性权重，wmax、wmin分别为惯性权重的最大、最小值；λ为非线性调节因子；t、tmax分别是当前进化代数和最大进化代数。在第t+1次迭代计算时，粒子i更新自己速度以及位置的规则如下：

式中，w指惯性权重；c1、c2指加速常数；r1，r2指两个在[0,1]内变化的随机数，且粒子的速度在一定的范围之内。

进一步的，一种基于改进粒子群算法的超低频振荡抑制方法，所述惯性权重的非线性减小采用的方法中引入非线性调节因子λ，并且根据不同的优化目标及条件，可改变λ的取值，来调整惯性权重w非线性减小的程度。通过引入非线性调节因子λ，可以通过改变λ的取值，来调整惯性权重w非线性减小的程度，增加算法的普遍适用性。

综上所述，本发明的以下有益效果：

1、本发明一种基于改进粒子群算法的超低频振荡抑制方法，以fc作为抑制超低频的主要手段，可以更加有效的对多直流送出系统的超低频振荡进行抑制，有效的保证系统频率的稳定。

2、本发明一种基于改进粒子群算法的超低频振荡抑制方法，通过引入非线性调节因子λ，可以有效的调整惯性权重w非线性减小的程度，使多直流送出系统的超低频振荡抑制变得可控。

3、本发明一种基于改进粒子群算法的超低频振荡抑制方法，通过将改进粒子群算法与控制敏感点相结合，可以有效的以较少的直流功率变化来实现对超低频振荡的有效抑制。

附图说明

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

图1为本发明流程示意图。

图2为某一电网实际结构示意图。

图3为标准粒子群算法的最优适用值的迭代曲线图。

图4为改进粒子群算法的最优适用值的迭代曲线图。

图5为系统的频率变化曲线。

图6为每条直流线路功率变化曲线。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。

实施例1

如图1所示，一种基于改进粒子群算法的超低频振荡抑制方法，所述方法包括的步骤为：s1：对多直流送出系统的每条直流分别施加相同的功率扰动，直至激发系统发生超低频振荡，并记录各条直流扰动下的超低频振荡模态幅值；s2：比较各条直流扰动下的超低频振荡模态幅值，并根据超低频振荡模态幅值选择直流频率限制控制器(frequencylimitcontroller，fc)的控制敏感点；s4：将所述fc的控制敏感点作为基础，通过tls-esprit算法辨识系统传递函数，建立系统频率偏差对时间的积分最小的粒子群算法优化模型；s5：对所述粒子群设定最大迭代次数、初始粒子数目、加速常数、惯性权重取值范围、非线性调节因子、维度及解空间范围，并将当前迭代次数设定为1；s6：采用粒子群算法对粒子的速度和位置进行更新；s7：记录每个微粒经历过的最好位置，所述微粒经历过的最好位置为个体极值，并将个体极值与微粒的适应值相比较，并根据每个微粒个体极值的位置与适应值的位置，确定出当前最好的个体极值；s8：将全局中微粒经历过的最好位置定为全局极值，并将每个微粒的适应值与全局极值作比较，并根据每个微粒全局极值的位置与适应值的位置，确定出当前最好的全局极值；s9：使惯性权重的非线性减小，并更新惯性权重取值；s10：判断当前迭代次数，如果当前迭代次数等于最大迭代次数，则输出fc参数优化结果，反之则将当前迭代次数加1，并返回步骤s6，步骤s9中的惯性权重的非线性减小采用的方法为：

式中，w代表惯性权重，wmax、wmin分别为惯性权重的最大、最小值；λ为非线性调节因子；t、tmax分别是当前进化代数和最大进化代数，所述惯性权重的非线性减小采用的方法中引入非线性调节因子λ，并且根据不同的优化目标及条件，可改变λ的取值，来调整惯性权重w非线性减小的程度。

实施例2

如图2～6所示，某一实际电网结构如图2所示，系统s1与系统s4通过三回±800kv的传统直流dc1、dc2、dc3相连；系统s1与系统s5通过一回±500kv的传统直流dc4相连；系统s1与系统s2通过500kv交流线路ac1、ac2、ac3相连，形成系统s12；btb为一柔直背靠背换流站，在其建成后，系统s12与外部系统s3、s4、s5的联系全部通过直流网络实现。在综合程序在综合程序psasp7.2中搭建该电网模型，考虑在丰期极小运行方式下，dc1、dc2、dc3三条特高压直流功率满送，分别为8000mw、7200mw、6400mw，而直流dc4零功率。考虑无fc的情况，分别对dc1、dc2、dc3三条直流施加200mw的功率扰动，得到系统的频率变化曲线，并利用tls-esprit算法进行分析，求得系统超低频主要振荡模态及该振荡频率下的振荡幅值，如表1所示。

表1超低频振荡模态

结果表明，当三条直流都施加200mw的功率扰动时，系统的超低频振荡模态约为0.077hz，对于该振荡模态，dc2对应的振荡幅值最大，dc1的振荡幅值次之，dc3的振荡幅值最小。因此，对于0.077hz的超低频振荡模态，dc2为控制敏感点。

进而以dc2为控制敏感点，辨识得到系统低阶线性化模型。

g(s)dc2＝(-0.3851s⁶-0.3434s⁵-0.6248s⁴-0.5813s³-0.122s²-0.2446s)

/(s⁶+0.4654s⁵+1.253s⁴+0.5422s³+0.265s²+0.1053s+0.0001466)

结合系统辨识得到的低阶线性模型，搭建优化模型，并在matlab中编写相关pso算法，以为优化的目标函数，并设置粒子群规模取30，最大迭代次数为50次，维度为2，解空间范围为[0，20]，粒子群速度范围为[-1，1]，加速常数惯性权重wmin＝0.4，wmax＝1.1，非线性调节因子λ＝0.2。得到dc2优化后的fc参数，比例系数为1.5210，时间常数为0.4201。

图3、图4分别是采用改进粒子群算法与标准粒子群算法时，最优适用值的迭代曲线。显然相比标准粒子群算法中惯性权重所采用的线性递减策略，改进粒子群算法中惯性权重所采用的非线性递减策略能够明显提升算法的收敛速度，且能保持较好的精度。

同时进一步说明dc2的fc用于抑制0.077hz超低频振荡模态最有效，对dc1、dc3两条直流也同样采用粒子群算法优化其fc参数。并设置ac3单回输电线路反生三相短路接地故障，以激发系统的超低频振荡。

由图5可以看到，当系统投入fc后，系统的超低频振荡都得到了明显抑制。dc1直流与dc2直流分别投入fc后，系统的频率变化相近，而dc3直流抑制效果相对较差。

由图6可得知，针对相同的故障三条直流的有功功率变化曲线，通过对时间进行积分，得到dc2能量变化1009.2kw·h，dc1能量变化1166.9kw·h，dc3能量变化1237.9kw·h。因此，dc2的相比于其它两条直流，其能量变化最小，也就是说dc2用了最小的能量达到了较好的抑制超低频振荡的效果。同时，若还需考虑对受端电网影响，直流功率变化也需尽可能小。

以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。