电力系统阻尼比在线趋优方法、存储介质和计算设备

摘要

本发明公开了一种电力系统阻尼比在线趋优方法、存储介质和计算设备，方法包括先获取电力系统在t时刻下的运行状态，得到实时运行数据；对电力系统的状态空间进行辨识，获取电力系统中的主导模式，并计算出在t时刻下主导模式的阻尼比；然后判断当前时刻主导模式的阻尼比是否高于设定的安全阈值，在是的情况下对当前时刻是否达到设定时间进行判断；在否的情况下根据实时运行数据构建阻尼比在线趋优的优化模型；接着求解优化模型，得到优化后的阻尼比；最后判断当前时刻是否达到设定时间，若是，则结束趋优；若否，则令t＝t+1，并继续获取实时运行数据进行阻尼比的趋优。本发明利用运行数据即可提升主导模式的阻尼比，可广泛应用且提高了小干扰稳定。

权利要求

1.一种电力系统阻尼比在线趋优方法，其特征在于，步骤如下：S1、通过实际电力系统中已有的广域监测系统和数据采集与监视控制系统获取电力系统在t时刻下的运行状态，得到实时运行数据；

S2、对电力系统的状态空间进行辨识，以获取电力系统中的主导模式，并计算出在t时刻下主导模式的阻尼比ξ[t]；

S3、判断当前时刻主导模式的阻尼比是否高于设定的安全阈值，若是，则跳至S6；若否，则跳至S4；

S4、根据实时运行数据构建阻尼比在线趋优的优化模型，该优化模型具体如下：主导模式的阻尼比ξ[t]与状态参量x的函数关系表示为：在函数关系式中，状态参量x是指除了平衡机以外的所有发电机的有功出力，x是列向量，维度为p×1，p是发电机的总台数；

是阻尼比灵敏度，是行向量，维度为1×p，行向量的求解过程如下：在[t‑1,t)时间段内获取N种运行方式的阻尼比值和机组出力，并分别组成对应的矩阵ξ和X，矩阵ξ的维度为N×1；矩阵X的维度为N×p；

利用局部加权回归分析方法获得阻尼比灵敏度其中，W为权重矩阵，T表示矩阵转置；

目标函数为：

Δx为电力系统中机组的调节量；

根据机组出力平衡约束和机组出力限制约束以及机组爬坡速率约束确定机组出力调整的上下限，其中，机组出力平衡约束为： i是发电机的序号；Δxi是第i台发电机的有功出力调节量；

机组出力限制约束为：Δxl≤Δx≤Δxu；Δxl为各个机组容量的20％与当前出力的差值；Δxu为各个机组的容量与当前出力的差值；

机组爬坡速率约束为：Δx≤|Δxramp|；其中，Δxramp为各个机组的爬坡速率，| |为绝对值；

S5、对优化模型进行求解，得到优化后的阻尼比ξ；

S6、判断当前时刻是否达到设定时间，若是，则结束趋优；

若否，则令t＝t+1，并返回步骤S1继续进行阻尼比的趋优。

2.根据权利要求1所述的电力系统阻尼比在线趋优方法，其特征在于，实时运行数据包括系统的潮流分布、各母线的电压、电力系统初始状态的阻尼比。

3.根据权利要求1所述的电力系统阻尼比在线趋优方法，其特征在于，在步骤S2中，具体是采用随机子空间辨识方法对电力系统的状态空间进行辨识；

所获取到的电力系统中的主导模式λ为：λ＝α+jβ；

根据阻尼比计算公式计算出在t时刻下主导模式的阻尼比ξ[t]，计算公式如下：其中，j是虚数单位；α是λ的实部；β是λ的虚部。

4.根据权利要求1所述的电力系统阻尼比在线趋优方法，其特征在于，安全阈值位于

0.03～0.04这一范围。

5.根据权利要求1所述的电力系统阻尼比在线趋优方法，其特征在于，在步骤S5中，采用单纯形法对优化模型进行求解，得到优化后状态参量x的值，将该值代入函数关系的表达式中，得到优化后的阻尼比ξ。

6.根据权利要求1所述的电力系统阻尼比在线趋优方法，其特征在于，设定时间设置为若干天或者选取一天中的时间段。

7.一种存储介质，存储有程序，其特征在于，所述程序被处理器执行时，实现权利要求1至6中任一项所述的电力系统阻尼比在线趋优方法。

8.一种计算设备，包括处理器以及用于存储处理器可执行程序的存储器，其特征在于，所述处理器执行存储器存储的程序时，实现权利要求1至6中任一项所述的电力系统阻尼比在线趋优方法。

说明书

电力系统阻尼比在线趋优方法、存储介质和计算设备

技术领域

[0001] 本发明涉及电力系统辨识和控制处理交叉技术领域，特别涉及一种电力系统阻尼比在线趋优方法、存储介质和计算设备。

背景技术

[0002] 现代电力系统正发生一系列重大变革：在电源侧，以风、光电为主的可再生能源高比例渗入，在降低电力碳排放率的同时，却加强了发电侧的随机波动性；在负荷侧，除了传
统负荷外，电动汽车作为新型负荷，其充放电行为给负荷侧带来了极大的随机性。未来的电
力系统将呈现出强随机性和不确定性，给电网运行调度及安全防御控制带来极大的挑战。

[0003] 在此背景下，电力系统振荡问题尤为突出，系统主导模式的最小阻尼比值在某些场景下可能会低于3％‑4％的安全范围。传统方法通常通过设置阻尼控制器或者调整系统
运行方式来提高系统阻尼。目前的阻尼控制器多设计于某一确定场景，难以在时变系统中
发挥出较好的效果。而现有的调整系统运行方式以离线灵敏度做支撑，对系统模型的要求
较高。这两者都需要准确的电力系统模型作为支撑，限制了其在新背景下的应用。

发明内容

[0004] 本发明的第一目的在于克服现有技术的缺点与不足，提供一种电力系统阻尼比在线趋优方法，该方法利用运行数据即可提升主导模式的阻尼比，可广泛应用且提高了小干
扰稳定。

[0005] 本发明的第二目的在于提供一种存储介质。

[0006] 本发明的第三目的在于提供一种计算设备。

[0007] 本发明的第一目的通过下述技术方案实现：一种电力系统阻尼比在线趋优方法，步骤如下：

[0008] S1、通过实际电力系统中已有的广域监测系统和数据采集与监视控制系统获取电力系统在t时刻下的运行状态，得到实时运行数据；

[0009] S2、对电力系统的状态空间进行辨识，以获取电力系统中的主导模式，并计算出在t时刻下主导模式的阻尼比；

[0010] S3、判断当前时刻主导模式的阻尼比是否高于设定的安全阈值，若是，则跳至S6；若否，则跳至S4；

[0011] S4、根据实时运行数据构建阻尼比在线趋优的优化模型；

[0012] S5、对优化模型进行求解，得到优化后的阻尼比ξ；

[0013] S6、判断当前时刻是否达到设定时间，

[0014] 若是，则结束趋优；

[0015] 若否，则令t＝t+1，并返回步骤S1继续进行阻尼比的趋优。

[0016] 优选的，实时运行数据包括系统的潮流分布、各母线的电压、电力系统初始状态的阻尼比。

[0017] 优选的，在步骤S2中，具体是采用随机子空间辨识方法对电力系统的状态空间进行辨识；

[0018] 所获取到的电力系统中的主导模式λ为：λ＝α+jβ；

[0019] 根据阻尼比计算公式计算出在t时刻下主导模式的阻尼比ξ[t]，计算公式如下：

[0020]

[0021] 其中，j是虚数单位；α是λ的实部；β是λ的虚部。

[0022] 优选的，安全阈值位于0.03～0.04这一范围。

[0023] 优选的，在步骤S4中，阻尼比值在线趋优的优化模型具体如下：

[0024] 主导模式的阻尼比ξ[t]与状态参量x的函数关系表示为：

[0025]

[0026] 在函数关系式中，状态参量x是指除了平衡机以外的所有发电机的有功出力； 是阻尼比灵敏度；

[0027] 目标函数为：

[0028]

[0029] Δx为电力系统中机组的调节量；

[0030] 根据机组出力平衡约束和机组出力限制约束以及机组爬坡速率约束确定机组出力调整的上下限，其中，机组出力平衡约束为： i是发电机的序号；p是发电机的
总台数；Δxi是第i台发电机的有功出力调节量；

[0031] 机组出力限制约束为：Δxl≤Δx≤Δxu；Δxl为各个机组容量的20％与当前出力的差值；Δxu为各个机组的容量与当前出力的差值；

[0032] 机组爬坡速率约束为：Δx≤|Δxramp|；其中，Δxramp为各个机组的爬坡速率，||为绝对值。

[0033] 更进一步的，状态参量x是列向量，维度为p×1；阻尼比灵敏度 是行向量，维度为1×p，行向量的求解过程如下：

[0034] 在[t‑1,t)时间段内获取N种运行方式的阻尼比值和机组出力，并分别组成对应的矩阵ξ和X，矩阵ξ的维度为N×1；矩阵X的维度为N×p；

[0035] 利用局部加权回归分析方法获得阻尼比灵敏度

[0036]

[0037] 其中，W为权重矩阵，T表示矩阵转置。

[0038] 更进一步的，在步骤S5中，采用单纯形法对优化模型进行求解，得到优化后状态参量x的值，将该值代入函数关系的表达式中，得到优化后的阻尼比ξ。

[0039] 优选的，设定时间设置为若干天或者选取一天中的时间段。

[0040] 本发明的第二目的通过下述技术方案实现：一种存储介质，存储有程序，所述程序被处理器执行时，实现本发明第一目的所述的电力系统阻尼比在线趋优方法。

[0041] 本发明的第三目的通过下述技术方案实现：一种计算设备，包括处理器以及用于存储处理器可执行程序的存储器，所述处理器执行存储器存储的程序时，实现本发明第一
目的所述的电力系统阻尼比在线趋优方法。

[0042] 本发明相对于现有技术具有如下的优点及效果：

[0043] 本发明电力系统阻尼比在线趋优方法，首先获取电力系统在t时刻下的运行状态，得到实时运行数据；对电力系统的状态空间进行辨识，获取电力系统中的主导模式，并计算
出在t时刻下主导模式的阻尼比；然后判断当前时刻主导模式的阻尼比是否高于设定的安
全阈值，在是的情况下对当前时刻是否达到设定时间进行判断；在否的情况下根据实时运
行数据构建阻尼比在线趋优的优化模型；接着求解优化模型，得到优化后的阻尼比；最后判
断当前时刻是否达到设定时间，若是，则结束趋优；若否，则令t＝t+1，并继续获取实时运行
数据进行阻尼比的趋优。本发明方法针对电力系统由于新能源功率波动而出现的某些弱阻
尼场景，兼顾机组自身及调整过程中存在的约束，科学量化了机组出力的调整量，在经过机
组出力调整，实现了电力系统主导模式的最小阻尼比值的提升，进而提高电力系统的小干
扰稳定，适合广泛使用。并且，本发明方法由运行数据驱动，通过根据运行数据构建的优化
模型进行阻尼比的趋优，相较于依赖离线模型的传统方法，只利用了系统在某段时间内的
运行数据，并不需要测试系统的完整模型，因此能够克服传统方法离线模型与在线实际不
匹配的问题。

附图说明

[0044] 图1是本发明电力系统阻尼比在线趋优方法的流程图。

[0045] 图2是二区四机测试系统的结构图。

[0046] 图3是图2系统的主导模式阻尼比值的变化图。

[0047] 图4是图2系统的机组出力的调整图。

具体实施方式

[0048] 下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。

[0049] 实施例1

[0050] 本实施例公开了一种电力系统阻尼比在线趋优方法，如图1所示，步骤如下：

[0051] S1、通过实际电力系统中已有的广域监测系统和数据采集与监视控制系统获取电力系统在t时刻下的运行状态，得到实时运行数据。

[0052] 实时运行数据具体包括系统的潮流分布、各母线的电压、电力系统初始状态的阻尼比。

[0053] S2、对电力系统的状态空间进行辨识，以获取电力系统中的主导模式，并计算出在t时刻下主导模式的阻尼比。

[0054] 其中，状态空间的辨识采用已知的随机子空间辨识方法(SSI，Stochastic Subspace Identification)，从电力系统的输出数据中辨识出系统的随机状态方程，求出
状态方程中的状态矩阵，根据状态矩阵得到主导模式。

[0055] 电力系统中的主导模式λ具体为：λ＝α+jβ；

[0056] 根据阻尼比计算公式计算出在t时刻下主导模式的阻尼比ξ[t]，该模式下的阻尼比值最小，计算公式如下：

[0057]

[0058] 其中，j是虚数单位；α是λ的实部；β是λ的虚部。

[0059] S3、判断当前时刻主导模式的阻尼比是否高于设定的安全阈值，若是，则跳至S6；若否，则跳至S4。本实施例的安全阈值位于0.03～0.04这一范围。

[0060] S4、根据实时运行数据构建阻尼比在线趋优的优化模型，该优化模型具体如下：

[0061] (1)主导模式的阻尼比ξ[t]与状态参量x的函数关系可用下式表示：

[0062]

[0063] 在函数关系式中，状态参量x是指除了平衡机以外的所有发电机的有功出力；x是列向量，维度为p×1，p是发电机的总台数； 是阻尼比灵敏度，是行向量，维度为1×p，行
向量的求解过程如下：

[0064] 在[t‑1,t)时间段内获取N种运行方式的阻尼比值和机组出力，并分别组成对应的矩阵ξ和X，矩阵ξ的维度为N×1；矩阵X的维度为N×p；本实施例N取200。

[0065] 利用局部加权回归分析方法获得阻尼比灵敏度

[0066]

[0067] 其中，W为权重矩阵，T表示矩阵转置。

[0068] (2)目标函数为：

[0069]

[0070] Δx为电力系统中机组的调节量；

[0071] (3)根据机组出力平衡约束和机组出力限制约束以及机组爬坡速率约束确定机组出力调整的上下限，其中，机组出力平衡约束为： i是发电机的序号；Δxi是第i
台发电机的有功出力调节量；

[0072] 机组出力限制约束为：Δxl≤Δx≤Δxu；Δxl为各个机组容量的20％与当前出力的差值；Δxu为各个机组的容量与当前出力的差值；

[0073] 机组爬坡速率约束为：Δx≤|Δxramp|；其中，Δxramp为各个机组的爬坡速率，||为绝对值。

[0074] S5、对优化模型进行求解，得到优化后的阻尼比ξ。

[0075] 在本实施例中，由于优化模型为线性模型，采用单纯形法可对优化模型进行求解，得到优化后状态参量x的值，将该值代入步骤S4函数关系的表达式中，即可得到优化后的阻
尼比ξ。

[0076] S6、判断当前时刻是否达到设定时间，设定时间可以根据实际情况设置为若干天例如一天，或者选取一天中的时间段，

[0077] 若是，则结束趋优；

[0078] 若否，则令t＝t+1，并返回步骤S1继续进行阻尼比的趋优。

[0079] 如图2所示，本实施例的电力系统为二区四机测试系统，包括4台发电机、2处负荷和2个并联电容，二区即是指区域1和区域2，该系统的额定频率为50Hz，主要电压等级为230
和20kV，其中3号机为平衡机。状态参量x即是指除了3号平衡机G3以外的3台发电机(G1机组
1、G2机组2、G4机组4)的有功出力。

[0080] 图3是初始运行方式和优化运行方式下的主导模式阻尼比的变化图，其中较为平缓的曲线是初始运行方式下的主导模式阻尼比的变化曲线，可以看出由于功率波动，原始
制定的运行方式出现了弱阻尼状态，通过本实施例的阻尼比趋优方法，改变机组的出力，得
到优化运行方式下的主导模式阻尼比的变化曲线，参见图3中剩余的那条斜率较大的曲线。
结合两条曲线，可以得出该二区四机测试系统经过调整，主导模式的阻尼比得到了明显的
上升。

[0081] 图4是机组出力调整过程，通过图4展示了各个机组每一次的出力调整量。

[0082] 在调整过程中，本实施例方法只利用了系统在某段时间内的运行数据，并不需要测试系统的完整模型，克服传统方法依赖于离线模型，离线模型与在线实际不匹配的问题。
同时，经过本实施例方法的调整，可以提升电力系统的小干扰稳定，因此本实施例方法适合
广泛应用。

[0083] 实施例2

[0084] 本实施例公开了一种存储介质，存储有程序，所述程序被处理器执行时，实现实施例1中所述的电力系统阻尼比在线趋优方法，具体如下：

[0085] S1、通过实际电力系统中已有的广域监测系统和数据采集与监视控制系统获取电力系统在t时刻下的运行状态，得到实时运行数据；

[0086] S2、对电力系统的状态空间进行辨识，以获取电力系统中的主导模式，并计算出在t时刻下主导模式的阻尼比；

[0087] S3、判断当前时刻主导模式的阻尼比是否高于设定的安全阈值，若是，则跳至S6；若否，则跳至S4；

[0088] S4、根据实时运行数据构建阻尼比在线趋优的优化模型；

[0089] S5、对优化模型进行求解，得到优化后的阻尼比ξ；

[0090] S6、判断当前时刻是否达到设定时间，

[0091] 若是，则结束趋优；

[0092] 若否，则令t＝t+1，并返回步骤S1继续进行阻尼比的趋优。

[0093] 本实施例中的存储介质可以是磁盘、光盘、计算机存储器、只读存储器(ROM，Read‑Only Memory)、随机存取存储器(RAM，Random Access Memory)、U盘、移动硬盘等介质。

[0094] 实施例3

[0095] 本实施例公开了一种计算设备，包括处理器以及用于存储处理器可执行程序的存储器，所述处理器执行存储器存储的程序时，实现实施例1中所述的电力系统阻尼比在线趋
优方法，具体如下：

[0096] S1、通过实际电力系统中已有的广域监测系统和数据采集与监视控制系统获取电力系统在t时刻下的运行状态，得到实时运行数据；

[0097] S2、对电力系统的状态空间进行辨识，以获取电力系统中的主导模式，并计算出在t时刻下主导模式的阻尼比；

[0098] S3、判断当前时刻主导模式的阻尼比是否高于设定的安全阈值，若是，则跳至S6；若否，则跳至S4；

[0099] S4、根据实时运行数据构建阻尼比在线趋优的优化模型；

[0100] S5、对优化模型进行求解，得到优化后的阻尼比ξ；

[0101] S6、判断当前时刻是否达到设定时间，

[0102] 若是，则结束趋优；

[0103] 若否，则令t＝t+1，并返回步骤S1继续进行阻尼比的趋优。

[0104] 本实施例中所述的计算设备可以是台式电脑、笔记本电脑、智能手机、PDA手持终端、平板电脑或其他具有处理器功能的终端设备。

[0105] 上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，
均应为等效的置换方式，都包含在本发明的保护范围之内。