CN115730838A - 基于并联局域模型的电力系统暂态电压安全裕度评估方法 - Google Patents

基于并联局域模型的电力系统暂态电压安全裕度评估方法 Download PDF

Info

Publication number
CN115730838A
CN115730838A CN202211264822.1A CN202211264822A CN115730838A CN 115730838 A CN115730838 A CN 115730838A CN 202211264822 A CN202211264822 A CN 202211264822A CN 115730838 A CN115730838 A CN 115730838A
Authority
CN
China
Prior art keywords
tvsi
transient voltage
voltage safety
index
safety margin
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211264822.1A
Other languages
English (en)
Other versions
CN115730838B (zh
Inventor
王亚军
杨立波
马斌
孙广辉
贾亚飞
栗维勋
李一鹏
袁龙
张宗杰
王晓蔚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Hebei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Hebei Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Hebei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202211264822.1A priority Critical patent/CN115730838B/zh
Publication of CN115730838A publication Critical patent/CN115730838A/zh
Application granted granted Critical
Publication of CN115730838B publication Critical patent/CN115730838B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/30Reactive power compensation

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提出一种基于并联局域模型的电力系统暂态电压安全裕度评估方法,以规则形式表达的无功源控制策略,不仅可以量化暂态电压安全裕度,而且可以辅助在线预防控制的决策。该控制策略是故障相依的,而且考虑了多个时间维度。针对并联局域模型输出存在误差的问题,提出了考虑置信水平的暂态电压安全裕度评估方法,可以得出不同保守程度下的暂态电压安全裕度评估结果,进而可以识别由于并联局域模型输出误差而被忽视的暂态电压安全隐患。

Description

基于并联局域模型的电力系统暂态电压安全裕度评估方法
技术领域
本发明属于电力技术领域,具体涉及一种基于并联局域模型的电力系统暂态电压安全裕度评估方法。
背景技术
暂态电压安全评估的目标是筛选出对电网暂态电压安全有威胁的严重故障,评估严重故障下的暂态电压安全裕度,为在线控制提供有效信息,进而改善电网的暂态电压安全性。
在现有暂态电压安全裕度评估研究中,主要使用基于数值仿真的方法或者基于人工智能的方法。由于数值仿真不能直接给出电网安全裕度信息,因此基于数值仿真的方法需要通过迭代搜索来确定暂态电压安全边界,计算复杂度高。基于人工智能的方法通常基于电网样本建立一个从电网特征变量到暂态电压安全裕度的映射关系,计算速度快。不过,当在线应用时机器学习模型可能面临泛化能力差、训练速度慢等问题。目前,现有暂态电压安全裕度评估研究中主要关注对电网运行状态与暂态电压安全边界之间距离的量化,在无功源控制量对暂态电压安全裕度影响方面的研究较少。
在暂态电压安全裕度评估方面,关键是量化电网运行状态与暂态电压安全边界之间的距离,以及为在线控制提供有效信息。暂态电压安全裕度评估结果是故障相依的,而且计算复杂度与需要考虑的故障数量成正比。现有研究方法中,基于数值仿真的方法计算复杂度高;基于人工智能的方法应用于在线时可能面临泛化能力差、训练速度慢等问题。此外,现有研究中主要关注对电网运行状态与暂态电压安全边界之间距离的量化,为在线控制提供的有效信息较少。因此,如何快速、准确地评估电网暂态电压安全裕度,并且辅助在线决策,是有待进一步研究的问题。
发明内容
为此,本发明提出一种基于并联局域模型的电力系统暂态电压安全裕度评估方法,以规则形式表达的无功源控制策略,不仅可以量化暂态电压安全裕度,而且可以辅助在线预防控制的决策。该控制策略是故障相依的,而且考虑了多个时间维度。针对并联局域模型输出存在误差的问题,提出了考虑置信水平的暂态电压安全裕度评估方法,可以得出不同保守程度下的暂态电压安全裕度评估结果,进而可以识别由于并联局域模型输出误差而被忽视的暂态电压安全隐患。
暂态电压安全裕度评估是在给定的电网运行状态、预想故障、电网参数和控制系统参数条件下,量化电网运行状态与暂态电压安全边界之间的距离,为在线控制提供有效信息。现有研究中,主要关注对电网运行状态与暂态电压安全边界之间的量化,对在线控制提供的有效信息较少。
(1)以规则形式表达的无功源控制策略
采用了暂态电压安全裕度的量化评估方法和无功源控制量对暂态电压安全裕度影响的量化评估方法,进而建立了以规则形式表达的无功源控制策略。
(1.1)暂态电压安全裕度的量化评估
在电网运行层面,故障后的暂态电压过程取决于电网运行状态x0、故障flt和电网拓扑参数topo。因此,存在从x0、flt和topo到TVSI指标的映射,表达式为:
Figure BDA0003892648110000021
电网的暂态电压安全性可以基于TVSI指标进行判定。对于预设的TVSI指标阈值TVSIth,以及给定的flt和topo,如果电网运行状态x0对应的TVSI指标满足式(2),则判定为暂态电压不安全;否则,判定为暂态电压安全。
Figure BDA0003892648110000022
其中,TVSI1、TVSI2和TVSI3分别是TVSI指标的三个分量,
Figure BDA0003892648110000023
Figure BDA0003892648110000024
分别是TVSI指标阈值的三个分量,∨是布尔代数中的“或”运算。
于是,可以得出基于TVSI指标的暂态电压安全边界,表达式为:
Figure BDA0003892648110000025
其中,
Figure BDA0003892648110000026
Figure BDA0003892648110000027
是从x0、flt和topo到TVSI指标第t个分量的映射,
Figure BDA0003892648110000028
是TVSI指标阈值的第t个分量。
Figure BDA0003892648110000029
是位于暂态电压安全边界上的电网运行状态,其物理含义为:TVSI指标某一分量等于TVSI指标阈值,而且TVSI指标其他分量不大于TVSI指标阈值的电网运行状态。
暂态电压安全裕度是对电网运行状态与暂态电压安全边界之间距离的量化,可以基于x0
Figure BDA00038926481100000210
之间的距离进行评估。在flt和topo固定的前提下,TVSI指标取决于x0。因此对于一组x0
Figure BDA0003892648110000031
它们对应的TVSI指标之差TVSI-TVSIth确定的。定义暂态电压安全裕度指标TVSIres,表达式如下:
Figure BDA0003892648110000032
如果TVSIres<0,则满足式(2),说明电网运行状态是暂态电压不安全的,而且TVSIres的值越小说明暂态电压安全性越差;否则,说明电网运行状态是暂态电压安全的,而且TVSIres的值越大说明暂态电压安全性越好。因此,TVSIres可以用于量化暂态电压安全裕度,TVSIres=0对应于暂态电压安全边界。
进一步定义与TVSI指标各个时间维度对应的暂态电压安全裕度指标
Figure BDA0003892648110000033
Figure BDA0003892648110000034
其中,TVSIt是TVSI指标第t个分量,
Figure BDA0003892648110000035
是TVSI指标第t个分量对应的暂态电压安全裕度,
Figure BDA0003892648110000036
是TVSI指标第t个分量的阈值。
注意,
Figure BDA0003892648110000037
是在给定flt和topo的前提下计算的,因此暂态电压安全裕度是故障相依的。本专利基于并联局域模型计算TVSI指标,进而用于评估暂态电压安全裕度,计算速度快。
(1.2)无功源控制量对暂态电压安全裕度影响的量化评估
无功源控制量u属于电网运行状态x0中的一部分。因此如果无功源控制量u改变,则故障后电网的暂态电压过程会发生改变,进而导致暂态电压安全裕度TVSIres发生改变。如果给定故障flt和电网拓扑参数topo,并保持电网其他运行控制量不变,则可以建立从u到TVSIres的映射fres,表达式如下:
TVSIres=fres(u) (7)
由于fres中包含了电网微分代数方程,因而难以解析求解,而且不够直观、解释性差,难以用于在线决策。为此,本专利基于fres的近似表达来量化评估无功源控制量u对暂态电压安全裕度TVSIres的影响。
对fres在无功源控制量初值u0附近进行泰勒展开,则可以得到TVSIres与u之间的近似函数关系式如下:
Figure BDA0003892648110000041
其中,
Figure BDA0003892648110000042
Figure BDA0003892648110000043
分别为无功源控制量为u0和u时TVSI指标第t个分量对应的暂态电压安全裕度。
式(8)中的项数越多,对fres的近似程度越高,但是计算复杂度更高、解释性更差。为了便于在线预防控制的决策,本专利仅保留式(8)中常数项和一次项用于近似表达fres,以确保较好的解释性,表达式如下:
Figure BDA0003892648110000044
根据式(6)可得:
Figure BDA0003892648110000045
其中,TVSIt是TVSI指标第t个分量。TVSIt对无功源控制量u的一阶偏导数可以通过下式近似计算:
Figure BDA0003892648110000046
其中,
Figure BDA0003892648110000047
是TVSI指标第t个分量对第i个控制量ui的一阶灵敏度,等于TVSIt对ui的一阶偏导数。
Figure BDA0003892648110000048
是当其他控制量不变时,第i个控制量ui增加Δui时TVSI指标的第t个分量;
Figure BDA0003892648110000051
是当其他控制量不变时,第i个控制量ui减少Δui时TVSI指标的第t个分量。
(1.3)以规则形式表达的无功源控制策略
基于式(9)-式(11)可以得出在无功源控制量初值u0附近成立的表达式如下:
Figure BDA0003892648110000052
其中,
Figure BDA0003892648110000053
Figure BDA0003892648110000054
分别是无功源控制量为u0和u时TVSI指标第t个分量对应的暂态电压安全裕度。
Figure BDA0003892648110000055
是TVSI指标第t个分量的阈值。
Figure BDA00038926481100000512
是无功源控制量为u0时TVSI指标第t个分量。
Figure BDA0003892648110000057
是TVSI指标第t个分量对控制量u的一阶灵敏度向量。
如果在线运行中要求暂态电压安全裕度TVSIres≥0,则可以得出以规则形式表达的无功源控制策略如下:
Figure BDA0003892648110000058
可见,无功源控制策略的形式是线性约束,规定了控制量u的可行域。其中,灵敏度向量
Figure BDA0003892648110000059
反映了无功源控制作用的强弱,可以用于选择需要控制的无功源;
Figure BDA00038926481100000510
Figure BDA00038926481100000511
之差反映了电网运行状态与暂态电压安全边界之间的距离,可以用于判定是否需要进行预防控制。因此,无功源控制策略可以辅助在线决策。
注意,式(13)是对应于TVSI指标第t个分量的无功源控制策略。可见,本专利提出的无功源控制策略是分多个时间维度的。如果TVSI三个分量对应的线性约束都成立,则必有TVSIres≥0。
在式(13)中,控制量u可以从电网现有无功源的可控量中选取。在建立无功源控制策略中,用到了fres的近似表达式,而fres是在给定故障flt和电网拓扑参数topo的前提下得出的。因此,无功源控制策略是故障相依的,可以根据故障筛选出控制策略中需要考虑的无功源,即故障关联无功源。考虑到容抗器、变压器等离散调节设备在一天内的调节次数有限制,不可以随意调节,因此在控制量u中仅考虑可以连续调节的动态无功源的可控量。动态无功源的可控量由其接入位置的节点类型决定。例如,如果发电机位于PV型节点,则其可控量是电压;如果发电机位于PQ型节点,则其可控量是无功出力。
由于无功源控制策略的建立过程中做了很多近似和简化,导致其仅在电网运行状态初值附近精度较高,需要针对电网运行状态的改变及时进行更新。本专利基于并联局域模型计算TVSI指标,而且建立无功源控制策略时只需要考虑严重故障,可以在数秒内完成更新,能够及时响应电网运行状态的改变。
(2)考虑置信水平的暂态电压安全裕度评估方法
由于并联局域模型输出存在误差,导致基于该模型计算的暂态电压安全裕度评估结果存在误差。如果上述误差导致对暂态电压安全裕度得出乐观的估计结果,则可能忽视电网中的暂态电压安全隐患,进而可能导致严重后果。为此,需要计算暂态电压安全裕度评估结果的置信区间下界,以识别由于并联局域模型输出误差而被忽视的暂态电压安全隐患。
在置信水平β下,暂态电压安全裕度置信区间下界的表达式为:
Figure BDA0003892648110000061
其中,
Figure BDA0003892648110000062
Figure BDA0003892648110000063
分别是TVSI指标第t个分量对应的暂态电压安全裕度的真实值和置信区间下界。
式(14)的物理含义为:暂态电压安全裕度真实值
Figure BDA0003892648110000064
不小于
Figure BDA0003892648110000065
的概率为β。因此,暂态电压安全裕度置信区间下界TVSIres,β可以看作是对暂态电压安全裕度的保守估计结果。
根据式(6),暂态电压安全裕度
Figure BDA0003892648110000066
的计算表达式中TVSI指标阈值
Figure BDA0003892648110000067
是设定值,因此
Figure BDA0003892648110000068
的计算误差取决于TVSI指标TVSIt的计算误差。为了计算暂态电压安全裕度置信区间下界TVSIres,β,需要计算TVSI指标的置信区间上界,表达式为:
Figure BDA0003892648110000069
其中,
Figure BDA0003892648110000071
Figure BDA0003892648110000072
和TVSIt分别是TVSI指标第t个分量的真实值和估计值,TVSI指标的估计值是并联局域模型的输出。
Figure BDA0003892648110000073
是TVSI指标第t个分量对应的误差概率密度函数。
Figure BDA0003892648110000074
即为TVSI指标第t个分量在置信水平β下的置信区间上界。
对于训练好的并联局域模型,其输出误差的概率密度函数
Figure BDA0003892648110000075
是固定的,因此式(16)中的
Figure BDA0003892648110000076
针对某一个置信水平β只需要计算一次即可。
Figure BDA0003892648110000077
的物理含义为TVSI指标第t个分量误差的置信区间上界。
把式(15)代入式(6)可得:
Figure BDA0003892648110000078
进一步得出暂态电压安全裕度置信区间下界TVSIres,β的计算表达式为:
Figure BDA0003892648110000079
附图说明
图1为本发明暂态电压安全裕度置信区间下界的示意图;
图2(a)为实施例故障发生后的电压波形初始条件;
图2(b)为实施例故障发生后的电压波形在发电机控制前和控制后对比;
图3为实施例发生故障后的电压波形。
具体实施方式
结合实施例说明本发明的具体技术方案。
在电网运行状态x0、故障flt和电网拓扑参数topo不变的前提下,改变置信水平β会得出不同的暂态电压安全裕度置信区间下界TVSIres,β,如图1所示。其中,红色点表示暂态电压安全裕度的计算结果,黑色虚线框内的紫色线表示暂态电压安全裕度的置信区间下界。在置信水平为90%和95%时,
Figure BDA00038926481100000710
说明这两种置信水平下暂态电压安全裕度充足;在置信水平为99%时,
Figure BDA0003892648110000081
说明该置信水平下暂态电压安全裕度不足。可见,在不同的置信水平下,将会得出不同的暂态电压安全裕度评估结果。置信水平越高,则暂态电压安全裕度评估结果的保守程度越高,越不容易忽视电网中的暂态电压安全隐患。
本实施例中,展示了以规则形式表达的无功源控制策略的建立与控制效果和不同置信水平下的暂态电压安全裕度评估结果。
各算例的通用边界条件为:负荷模型是由感应电动机与恒阻抗负荷组合成的综合负荷,其中感应电动机占比60%。
(1):以规则形式表达的无功源控制策略的建立和控制效果
本算例采用IEEE39节点模型。考虑线路“Bus6-Bus7”上的三相N-1故障,故障持续时间为250ms。考虑初始的电网拓扑参数和电网运行状态。
IEEE39节点系统中包含10台发电机,无其他无功源。本算例中,把10台发电机都作为所考虑故障的关联无功源,控制量为发电机的机端电压。初始条件下,发电机机端电压如表1所示,所考虑故障发生后电网TVSI指标的三个分量分别为0.2490、0.2896和0.1513。TVSI指标与发电机机端电压之间的灵敏度系数如表2所示。
表1发电机机端电压
Figure BDA0003892648110000082
表2 TVSI指标与发电机机端电压之间的灵敏度系数
Figure BDA0003892648110000083
TVSI指标三个分量阈值分别设定为
Figure BDA0003892648110000084
Figure BDA0003892648110000085
则可以得出所考虑故障对应的无功源控制策略如下:
Figure BDA0003892648110000091
Figure BDA0003892648110000092
Figure BDA0003892648110000093
其中,V30~V39分别表示Bus30~Bus39处发电机的机端电压,
Figure BDA0003892648110000094
分别表示Bus30~Bus39处发电机机端电压的初值。式(19)-式(21)分别是TVSI指标三个分量对应的无功源控制策略。
可见,无功源控制策略包含多个时间维度。初始条件下,所考虑故障发生后的电压波形如图2(a)所示,其中展示了暂态电压安全问题最严重母线的电压波形。可见,故障后电压的恢复速度较慢,而且故障后的稳态电压低于故障前的稳态电压,说明故障引起了低电压问题。在TVSI2和TVSI3对应的无功源控制策略中,各发电机机端电压的系数都小于0,说明提高发电机机端电压可以提升电网电压水平,并且使暂态电压安全性变好。在TVSI1对应的无功源控制策略中,各发电机机端电压的系数都大于0,主要原因是当电网电压水平提升时,故障中的电压跌落幅度也会小幅度增加。因此,计算得到的无功源控制策略符合逻辑。
对比上述三个无功源控制策略中的常数项,可以发现TVSI2对应的控制策略中TVSI指标初值与TVSI指标阈值最接近,因此应该被重点关注。为了展示无功源控制策略在预防控制方面的作用,把与TVSI2之间灵敏度系数绝对值最大的三台发电机(Bus31、Bus32和Bus39)的机端电压提升到1.05p.u.,以减小TVSI2,增加电网暂态电压安全裕度。在发电机控制前和控制后,所考虑故障发生后的电压波形对比如图2(b)所示。相比之下,发电机控制后的电压波形在故障清除后恢复得更快、故障后稳态电压与故障前稳态电压更接近。发电机控制后的TVSI2为0.2306,比控制前下降了0.059,说明预防控制起到了预期效果,无功源控制策略可以辅助预防控制的决策。
此外,根据IEEE39节点模型接线图,可以发现本算例所考虑的故障线路“Bus6-Bus7”位于该电网模型的左下部分,与该故障电气距离最接近的三台发电机是位于Bus31、Bus32和Bus39的三台发电机。因此,这三台发电机机端电压与TVSI2之间的灵敏度系数绝对值最大符合逻辑,从另一方面说明了本专利提出的无功源控制策略的合理性。
(2):以规则形式表达的无功源控制策略的建立
本算例采用陕西电网模型。考虑线路“西北.7156横店线”上的三相N-1故障,故障持续时间为150ms。考虑初始的电网拓扑参数和电网运行状态。
经过分析,该故障的关联无功源中包含18台发电机,关联母线中有39个330kV电压等级以上的母线,控制量为发电机的机端电压。初始条件下,故障的18台发电机机端电压如表3所示,所考虑故障发生后电网的TVSI指标的三个分量分别为0.1487、0.5583和0.0574。TVSI指标与发电机机端电压之间的灵敏度系数如表4所示。
表3发电机机端电压
Figure BDA0003892648110000101
表4 TVSI指标与发电机机端电压之间的灵敏度系数
Figure BDA0003892648110000102
Figure BDA0003892648110000111
TVSI指标三个分量阈值分别设定为
Figure BDA0003892648110000112
Figure BDA0003892648110000113
所考虑的故障在初始运行场景下的无功源控制策略与式(19)-式(21)中的形式相似。由于本算例中故障关联无功源的数量较多,故不在此详细列出控制策略的表达式。
所考虑故障发生后的电压波形如图3所示。可见,故障后电压的恢复速度较慢,说明故障引起了低电压问题。在TVSI2和TVSI3对应的无功源控制策略中,各发电机机端电压的系数都小于0,说明提高发电机机端电压可以提升电网电压水平,并且使暂态电压安全性变好。在TVSI1对应的无功源控制策略中,各发电机机端电压的系数都大于0,主要原因是当电网电压水平提升时,故障中的电压跌落幅度也会小幅度增加。因此,计算得到的无功源控制策略符合逻辑。
此外,从表4中可以发现部分发电机对应的灵敏度系数相同。经过进一步分析,发现灵敏度系数相同的发电机都隶属于同一个高压母线,或者是在同一个电厂内,因此得出对TVSI指标灵敏度相同的结果是合理的,从另一方面说明了本专利提出的无功源控制策略的合理性。(3):不同置信水平下的暂态电压安全裕度置信区间下界
本算例采用IEEE39节点模型。考虑线路“Bus6-Bus7”上的三相N-1故障,故障持续时间为250ms。考虑初始的电网拓扑参数和电网运行状态。
初始条件下,该故障发生后电网TVSI指标的三个分量分别为0.2490、0.2896和0.1513。TVSI指标三个分量阈值分别设定为
Figure BDA0003892648110000114
Figure BDA0003892648110000115
则可以得出暂态电压安全裕度评估结果TVSIres和不同置信水平下暂态电压安全裕度置信区间下界TVSIres,β,如表5所示。
表5不同置信水平下的暂态电压安全裕度评估结果
Figure BDA0003892648110000116
Figure BDA0003892648110000121
在置信水平为90%时,TVSI指标三个分量对应的暂态电压安全裕度置信区间下界都大于0,说明该置信水平下暂态电压安全裕度充足;在置信水平为95%和99%时,TVSI2对应的暂态电压安全裕度置信区间下界都小于0,说明这两种置信水平下暂态电压安全裕度不足,电网中可能存在暂态电压安全隐患。可见,在不同的置信水平下,得出了不同的暂态电压安全裕度评估结果。

Claims (5)

1.基于并联局域模型的电力系统暂态电压安全裕度评估方法,其特征在于,包括以下步骤:
(1)建立以规则形式表达的无功源控制策略
采用暂态电压安全裕度的量化评估方法和无功源控制量对暂态电压安全裕度影响的量化评估方法,建立以规则形式表达的无功源控制策略;
(2)考虑置信水平的暂态电压安全裕度评估方法
计算暂态电压安全裕度评估结果的置信区间下界,以识别由于并联局域模型输出误差而被忽视的暂态电压安全隐患。
2.根据权利要求1所述的基于并联局域模型的电力系统暂态电压安全裕度评估方法,其特征在于,步骤(1)建立以规则形式表达的无功源控制策略中,暂态电压安全裕度的量化评估方法,包括:
在电网运行层面,故障后的暂态电压过程取决于电网运行状态x0、故障flt和电网拓扑参数topo;存在从x0、flt和topo到TVSI指标的映射,表达式为:
Figure FDA0003892648100000011
电网的暂态电压安全性基于TVSI指标进行判定;对于预设的TVSI指标阈值TVSIth,以及给定的flt和topo,如果电网运行状态x0对应的TVSI指标满足式(2),则判定为暂态电压不安全;否则,判定为暂态电压安全;
Figure FDA0003892648100000012
其中,TVSI1、TVSI2和TVSI3分别是TVSI指标的三个分量,
Figure FDA0003892648100000013
Figure FDA0003892648100000014
分别是TVSI指标阈值的三个分量,∨是布尔代数中的“或”运算;
得出基于TVSI指标的暂态电压安全边界,表达式为:
Figure FDA0003892648100000015
其中,
Figure FDA0003892648100000016
Figure FDA0003892648100000017
是从x0、flt和topo到TVSI指标第t个分量的映射,
Figure FDA0003892648100000018
是TVSI指标阈值的第t个分量;
Figure FDA0003892648100000019
是位于暂态电压安全边界上的电网运行状态,其物理含义为:TVSI指标某一分量等于TVSI指标阈值,而且TVSI指标其他分量不大于TVSI指标阈值的电网运行状态;
暂态电压安全裕度是对电网运行状态与暂态电压安全边界之间距离的量化,基于x0
Figure FDA0003892648100000021
之间的距离进行评估;在flt和topo固定的前提下,TVSI指标取决于x0;对于一组x0
Figure FDA0003892648100000022
它们对应的TVSI指标之差TVSI-TVSIth确定的;定义暂态电压安全裕度指标TVSIres,表达式如下:
Figure FDA0003892648100000023
如果TVSIres<0,则满足式(2),说明电网运行状态是暂态电压不安全的,而且TVSIres的值越小说明暂态电压安全性越差;否则,说明电网运行状态是暂态电压安全的,而且TVSIres的值越大说明暂态电压安全性越好;TVSIres用于量化暂态电压安全裕度,TVSIres=0对应于暂态电压安全边界;
定义与TVSI指标各个时间维度对应的暂态电压安全裕度指标
Figure FDA0003892648100000024
Figure FDA0003892648100000025
其中,TVSIt是TVSI指标第t个分量,
Figure FDA0003892648100000026
是TVSI指标第t个分量对应的暂态电压安全裕度,
Figure FDA0003892648100000027
是TVSI指标第t个分量的阈值。
3.根据权利要求2所述的基于并联局域模型的电力系统暂态电压安全裕度评估方法,其特征在于,步骤(1)建立以规则形式表达的无功源控制策略中,无功源控制量对暂态电压安全裕度影响的量化评估包括以下步骤:
无功源控制量u属于电网运行状态x0中的一部分;给定故障flt和电网拓扑参数topo,并保持电网其他运行控制量不变,则建立从u到TVSIres的映射fres,表达式如下:
TVSIres=fres(u) (7)
基于fres的近似表达来量化评估无功源控制量u对暂态电压安全裕度TVSIres的影响;
对fres在无功源控制量初值u0附近进行泰勒展开,得到TVSIres与u之间的近似函数关系式如下:
Figure FDA0003892648100000031
保留式(8)中常数项和一次项用于近似表达fres,以确保较好的解释性,表达式如下:
Figure FDA0003892648100000032
根据式(6)可得:
Figure FDA0003892648100000033
其中,TVSIt是TVSI指标第t个分量;TVSIt对无功源控制量u的一阶偏导数通过下式近似计算:
Figure FDA0003892648100000034
其中,
Figure FDA0003892648100000035
是TVSI指标第t个分量对第i个控制量ui的一阶灵敏度,等于TVSIt对ui的一阶偏导数;
Figure FDA0003892648100000036
是当其他控制量不变时,第i个控制量ui增加Δui时TVSI指标的第t个分量;
Figure FDA0003892648100000037
是当其他控制量不变时,第i个控制量ui减少Δui时TVSI指标的第t个分量。
4.根据权利要求3所述的基于并联局域模型的电力系统暂态电压安全裕度评估方法,其特征在于,步骤(1)建立以规则形式表达的无功源控制策略中,以规则形式表达的无功源控制策略包括以下步骤:
基于式(9)-式(11)得出在无功源控制量初值u0附近成立的表达式如下:
Figure FDA0003892648100000041
其中,
Figure FDA0003892648100000042
Figure FDA0003892648100000043
分别是无功源控制量为u0和u时TVSI指标第t个分量对应的暂态电压安全裕度;
Figure FDA0003892648100000044
是TVSI指标第t个分量的阈值;
Figure FDA00038926481000000419
是无功源控制量为u0时TVSI指标第t个分量;
Figure FDA0003892648100000046
是TVSI指标第t个分量对控制量u的一阶灵敏度向量;
如果在线运行中要求暂态电压安全裕度TVSIres≥0,则得出以规则形式表达的无功源控制策略如下:
Figure FDA0003892648100000047
灵敏度向量
Figure FDA0003892648100000048
反映了无功源控制作用的强弱,用于选择需要控制的无功源;
Figure FDA0003892648100000049
Figure FDA00038926481000000410
之差反映了电网运行状态与暂态电压安全边界之间的距离,用于判定是否需要进行预防控制。
5.根据权利要求1所述的基于并联局域模型的电力系统暂态电压安全裕度评估方法,其特征在于,步骤(2)考虑置信水平的暂态电压安全裕度评估方法,包括以下步骤:
在置信水平β下,暂态电压安全裕度置信区间下界的表达式为:
Figure FDA00038926481000000411
其中,
Figure FDA00038926481000000412
Figure FDA00038926481000000413
分别是TVSI指标第t个分量对应的暂态电压安全裕度的真实值和置信区间下界;
式(14)的物理含义为:暂态电压安全裕度真实值
Figure FDA00038926481000000414
不小于
Figure FDA00038926481000000415
的概率为β;因此,暂态电压安全裕度置信区间下界TVSIres,β看作是对暂态电压安全裕度的保守估计结果;
根据式(6),暂态电压安全裕度
Figure FDA00038926481000000416
的计算表达式中TVSI指标阈值
Figure FDA00038926481000000417
是设定值,因此
Figure FDA00038926481000000418
的计算误差取决于TVSI指标TVSIt的计算误差;为了计算暂态电压安全裕度置信区间下界TVSIres,β,需要计算TVSI指标的置信区间上界,表达式为:
Figure FDA0003892648100000051
其中,
Figure FDA0003892648100000052
Figure FDA0003892648100000053
和TVSIt分别是TVSI指标第t个分量的真实值和估计值,TVSI指标的估计值是并联局域模型的输出;
Figure FDA0003892648100000054
是TVSI指标第t个分量对应的误差概率密度函数;
Figure FDA0003892648100000055
即为TVSI指标第t个分量在置信水平β下的置信区间上界;
对于训练好的并联局域模型,其输出误差的概率密度函数
Figure FDA0003892648100000056
是固定的,因此式(16)中的
Figure FDA0003892648100000057
针对某一个置信水平β只需要计算一次即可;
Figure FDA0003892648100000058
的物理含义为TVSI指标第t个分量误差的置信区间上界;
把式(15)代入式(6)可得:
Figure FDA0003892648100000059
进一步得出暂态电压安全裕度置信区间下界TVSIres,β的计算表达式为:
Figure FDA00038926481000000510
在电网运行状态x0、故障flt和电网拓扑参数topo不变的前提下,改变置信水平β会得出不同的暂态电压安全裕度置信区间下界TVSIres,β
CN202211264822.1A 2022-10-17 2022-10-17 基于并联局域模型的电力系统暂态电压安全裕度评估方法 Active CN115730838B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211264822.1A CN115730838B (zh) 2022-10-17 2022-10-17 基于并联局域模型的电力系统暂态电压安全裕度评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211264822.1A CN115730838B (zh) 2022-10-17 2022-10-17 基于并联局域模型的电力系统暂态电压安全裕度评估方法

Publications (2)

Publication Number Publication Date
CN115730838A true CN115730838A (zh) 2023-03-03
CN115730838B CN115730838B (zh) 2023-11-07

Family

ID=85293653

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211264822.1A Active CN115730838B (zh) 2022-10-17 2022-10-17 基于并联局域模型的电力系统暂态电压安全裕度评估方法

Country Status (1)

Country Link
CN (1) CN115730838B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1835336A (zh) * 2006-03-10 2006-09-20 天津大学 基于概率的电力系统电压稳定在线监视和预防控制方法
CN102324743A (zh) * 2011-09-21 2012-01-18 国网电力科学研究院 电力系统在线暂态安全稳定评估预想故障筛选方法
CN109638815A (zh) * 2018-12-04 2019-04-16 国电南瑞科技股份有限公司 一种电力系统中长期电压安全稳定预防控制策略确定方法
CN112217207A (zh) * 2020-10-20 2021-01-12 南方电网科学研究院有限责任公司 一种暂态电压稳定裕度预测方法及装置
US11362539B1 (en) * 2021-11-29 2022-06-14 King Abdulaziz University Artificial intelligence-based power controller for low voltage ride-through control of grid connected distributed generation networks

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1835336A (zh) * 2006-03-10 2006-09-20 天津大学 基于概率的电力系统电压稳定在线监视和预防控制方法
CN102324743A (zh) * 2011-09-21 2012-01-18 国网电力科学研究院 电力系统在线暂态安全稳定评估预想故障筛选方法
CN109638815A (zh) * 2018-12-04 2019-04-16 国电南瑞科技股份有限公司 一种电力系统中长期电压安全稳定预防控制策略确定方法
CN112217207A (zh) * 2020-10-20 2021-01-12 南方电网科学研究院有限责任公司 一种暂态电压稳定裕度预测方法及装置
US11362539B1 (en) * 2021-11-29 2022-06-14 King Abdulaziz University Artificial intelligence-based power controller for low voltage ride-through control of grid connected distributed generation networks

Also Published As

Publication number Publication date
CN115730838B (zh) 2023-11-07

Similar Documents

Publication Publication Date Title
CN111628501B (zh) 一种交直流大电网暂态电压稳定评估方法及系统
CN110879917A (zh) 一种基于迁移学习的电力系统暂态稳定自适应评估方法
JP5494034B2 (ja) 信頼度評価装置、信頼度評価プログラムおよび信頼度評価方法
CN107843816B (zh) 一种考虑负载率影响的变压器绝缘缺陷状态评估方法
CN104914327A (zh) 基于实时监测信息的变压器故障检修预测方法
CN103778575A (zh) 变压器状态评估方法及系统
CN101425686A (zh) 电力系统在线安全稳定评估预想故障集自适应筛选方法
CN107203690A (zh) 一种gm(1,1)模型变压器故障预测方法
CN105488344A (zh) 一种配电设备健康指数通用评价方法
CN110619321B (zh) 一种基于深度神经网络的充电桩技术状况诊断方法及系统
Silva et al. An accurate method for fault location in HVDC systems based on pattern recognition of DC voltage signals
CN116031867A (zh) 一种电网暂态安全稳定性评估方法
CN115481726A (zh) 一种工业机器人整机健康评估方法及系统
CN113300373B (zh) 一种基于prmse评价指标的稳定裕度值预测方法及装置
CN111062569A (zh) 一种基于bp神经网络的小电流故障判别方法
CN113610119A (zh) 基于卷积神经网络识别输电线路发展性故障的方法
CN105741184A (zh) 一种变压器状态评估方法及装置
CN115730838A (zh) 基于并联局域模型的电力系统暂态电压安全裕度评估方法
CN117335426A (zh) 一种适用于短期电压稳定性评估的关键特征识别方法
Tiwary et al. Multi-dimensional ANN application for active power flow state classification on a utility system
Zhang et al. CNN‐LSTM based power grid voltage stability emergency control coordination strategy
CN111654033B (zh) 基于线性推演的大电网静态电压稳定态势评估方法及系统
CN114236332A (zh) 一种电力电缆绝缘状态判断方法及系统
CN114418287A (zh) 一种含电-气-热耦合系统的微电网频率安全评估方法
Nandagopan et al. Online prediction of DGA results for intelligent condition monitoring of power transformers

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant