一种考虑大量预想故障的暂态电压安全预防控制优化方法转让专利
申请号 : CN202110295002.8
文献号 : CN113097997B
文献日 : 2022-05-13
发明人 : 郭庆来 , 张振安 , 王彬 , 崔惟 , 孙宏斌 , 单瑞卿 , 葛怀畅 , 林银鸿
申请人 : 清华大学 , 国网河南省电力公司电力科学研究院 , 国家电网有限公司
摘要 :
权利要求 :
1.一种考虑大量预想故障的暂态电压安全预防控制优化方法,其特征在于,该方法基于衡量节点的暂态电压跌落程度评估指标,根据初始的预想故障集,建立由目标函数和约束条件构成的考虑大量预想故障的暂态电压安全预防控制优化模型;然后通过在该模型的暂态电压安全约束中引入松弛变量,同时在目标函数中加入相应的惩罚项,对原模型进行松弛;最后通过逐次线性化的方法求解该松弛化后的模型,得到各无功设备的优化设定值;
该方法包括以下步骤:
(1)定义一种衡量节点的暂态电压跌落程度评估指标TVSI,如下式所示:式中,V0表示节点电压的初始值,v(t)表示故障发生后t时刻的节点电压,Tflt、Tend分别表示故障发生的时刻以及故障后节点电压第一次恢复到初始值V0的时刻;DPena表示惩罚函数;ΔT表示电压标幺值跌落至0.8以下时可以容忍的持续时间;v(Tflt+ΔT)表示故障发生经过ΔT后的节点电压;
(2)基于步骤(1)的暂态电压跌落程度评估指标,获取预想故障集SFlt并作为初始的预想故障集,建立考虑大量预想故障的暂态电压安全预防控制优化模型,该模型由目标函数和约束条件构成;具体步骤如下:
(2‑1)确定模型的目标函数如下:其中, 表示系统中发电机j的无功出力调整量,SG表示当前系统中的发电机集合;
表示容抗器s无功补偿容量的调整量,SC表示当前系统中的容抗器集合;ΔTh表示变压器h的分接头调整量,SH表示系统中可调节变比的变压器集合;ω1、ω2、ω3为权重系数;
(2‑2)确定模型的约束条件,具体如下:(2‑2‑1)发生预想故障前,电力系统的节点电压增量与无功设备的调整量近似线性化约束,表达式如下:
其中,ΔVi、Vi,0、Vi分别表示节点i的电压增量、电压初始值以及调整无功设备后的电压幅值; 分别表示发电机j调整前的无功出力和调整后的无功出力; 分别表示容抗器s调整前的无功补偿容量和调整后的无功补偿容量;Th,0、Th分别表示变压器h抽头调整前后的档位; 表示节点i的电压对于发电机j无功出力的偏导数, 表示节点i的电压对于容抗器s无功补偿容量的偏导数, 表示节点电压i对于变压器h分接头档位的偏导数;
(2‑2‑2)电力系统稳态运行时各节点电压幅值、发电机无功出力、容抗器无功补偿容量、变压器可调分接头满足的不等式约束如下:Vmin<Vi<Vmax,i=1,2…NbusTh,min<Th<Th,max,h=1,2…NH其中,Vmin、Vmax分别表示节点电压的下限和上限; 分别表示发电机j出力的下限和上限; 分别表示容抗器s无功补偿容量的下限和上限;Th,min、Th,max分别表示变压器h可调分接头的下限和上限;Nbus、NG、NC、NH分别表示电力系统中的节点、发电机、容抗器、分接头可调的变压器的数量;
(2‑2‑3)暂态电压安全约束;
电力系统发生故障k后,节点i在t时刻的电压 满足如下的暂态电压安全约束:其中, 表示在初始状态下当系统发生故障k之后,节点i的暂态电压跌落程度评估指标; 表示系统发生故障k之后,节点i的暂态电压跌落程度评估指标对发电机j无功出力的偏导数;TVSIth表示暂态电压跌落程度评估指标的安全门槛值;NFlt表示预想故障集SFlt的故障数量;
(3)将步骤(2)中的暂态电压安全约束进行松弛,如下式:ζk≥0
k=1,2…NFlt
其中,ζk表示为了满足故障k对应的暂态电压安全约束所引入的松弛变量;
同时,在目标函数中加入一个惩罚项,更新后的目标函数如下式:其中,M为正数; 表示惩罚项;
(4)求解经过步骤(3)松弛后的考虑大量预想故障的暂态电压安全预防控制优化模型,得到各无功设备的当前最优设定值,包括:各发电机的无功出力 各容抗器的无功补偿容量 各变压器分接头档位
调整系统中的各无功设备的运行点为当前最优设定值,然后进行时域仿真,扫描初始的预想故障集SFlt;
(5)根据初始的预想故障集SFlt的故障扫描结果,得到每种故障发生后各节点在各个时刻的电压 计算得到各故障下各节点的暂态电压跌落程度评估指标 通过对各进行判定,生成优化后不满足暂态电压安全约束的故障组成的集合为 如下式:其中,Fltk表示第k个故障;Sbus表示系统的节点集合;
(6)判断当前最优设定值是否满足如下式所示的收敛条件:其中ε表示收敛阈值;
如果满足,则求解结束,输出当前最优设定值 为最终的最优设定值,优化结束;否则,重新返回步骤(4),将步骤(5)得到的 作为更新后的SFlt,继续下一轮求解。
说明书 :
一种考虑大量预想故障的暂态电压安全预防控制优化方法
技术领域
背景技术
满足用电负荷的需求,目前已经基本形成了西电东送通道,在东部和南部负荷中心则形成
了若干个大型的多直流落点受端电网。这一系列新技术的引入使得电网中出现了大量的电
力电子设备,而这些设备对电压变化很敏感。与此同时,受端电网中动态负荷尤其是感应电
动机的比例也在日益增长。一旦受端电网发生短路故障,这些动态负荷以及直流换流站将
从交流电网中吸收大量无功,导致系统的无功电压支撑不足,容易出现电压失稳甚至崩溃,
造成大面积停电。因此,亟需针对受端电网采取考虑暂态电压安全约束的预防控制措施,实
现“未雨绸缪”。
开机不足,而发电机又是电力系统中重要的动态无功支撑设备,因此这将导致系统中的可
调无功源明显不足。此时如果直接针对所有的预想故障建立一个考虑暂态电压安全约束的
预防控制优化模型,很可能导致优化问题不可行,即在一些预想故障下的电压安全约束无
法同时得到满足,需要通过新增无功补偿设备来解决问题。由于新增设备的成本一般较高,
因此在这种情况下,如何通过预防控制尽可能满足更多故障的暂态电压安全约束就显得十
十分重要。然而,现有的研究中一般仅在优化模型中考虑少量指定的故障约束,或者是同时
考虑在电网中新增无功补偿设备,而对于如何通过预防控制来满足实际电网大量预想故障
安全约束的研究仍然较少。
发明内容
故障的安全约束,通过本发明得到的优化结果可降低直流换流站发生换流失败的风险。
标函数和约束条件构成的考虑大量预想故障的暂态电压安全预防控制优化模型;然后通过
在该模型的暂态电压安全约束中引入松弛变量,同时在目标函数中加入相应的惩罚项,对
原模型进行松弛;最后通过逐次线性化的方法求解该松弛化后的模型,得到各无功设备的
优化设定值。该方法包括以下步骤:
罚函数;ΔT表示电压标幺值跌落至0.8以下时可以容忍的持续时间;v(Tflt+ΔT)表示故障
发生经过ΔT后的节点电压;
函数和约束条件构成;具体步骤如下:
压器h的分接头调整量,SH表示系统中可调节变比的变压器集合;ω1、ω2、ω3为权重系数;
别表示容抗器s调整前的无功补偿容量和调整后的无功补偿容量;Th,0、Th分别表示变压器h
抽头调整前后的档位; 表示节点i的电压对于发电机j无功出力的偏导数, 表示节
点i的电压对于容抗器s无功补偿容量的偏导数, 表示节点电压i对于变压器h分接头档
位的偏导数;
别表示变压器h可调分接头的下限和上限;Nbus、NG、NC、NH分别表示电力系统中的节点、发电
机、容抗器、分接头可调的变压器的数量;
无功出力的偏导数;TVSIth表示暂态电压跌落程度评估指标的安全门槛值;NFlt表示预想故
障集SFlt的故障数量;
补偿容量 各变压器分接头档位
对各 进行判定,生成优化后不满足暂态电压安全约束的故障组成的集合为 如下
式:
轮求解。
系统的暂态电压安全水平,降低直流换流站发生换相失败的风险。
具体实施方式
的暂态电压跌落程度评估指标TVSI,如下式所示:
罚函数:当节点电压不满足WECC的电压稳定判据时,该惩罚函数的取值为10,否则取值为0。
ΔT表示电压标幺值跌落至0.8以下时可以容忍的持续时间,一般取值为0.4s。v(Tflt+ΔT)
表示故障发生经过ΔT后的节点电压。
该模型由目标函数和约束条件构成;具体步骤如下:
表示当前系统中的容抗器集合。ΔTh表示变压器h的分接头调整量,SH表示系统中可调节变
比的变压器集合。ω1、ω2、ω3为权重系数,可以根据用户需求设定。
别表示容抗器s调整前的无功补偿容量和调整后的无功补偿容量。Th,0、Th分别表示变压器h
抽头调整前后的档位。 表示节点i的电压对于发电机j无功出力的偏导数, 表示节
点i的电压对于容抗器s无功补偿容量的偏导数, 表示节点电压i对于变压器h分接头档
位的偏导数。
容量的下限和上限,Th,min、Th,max分别表示变压器h可调分接头的下限和上限,这些限值由设
备本身的运行参数来决定。Nbus、NG、NC、NH分别表示电力系统中的节点、发电机、容抗器、分接
头可调的变压器的数量。
方法见步骤(1)。 表示系统发生故障k之后,节点i的暂态电压跌落程度评估指标对发
电机j无功出力的偏导数。TVSIth表示暂态电压跌落程度评估指标的安全门槛值,取值为10。
NFlt表示预想故障集SFlt的故障数量。
行适当松弛,如下式:
补偿容量 各变压器分接头档位 其中,若是第一次求解该模型,则该模型的预想
故障集采用初始的预想故障集SFlt;若不是第一次求解该模型,则该模型的故障集采用前一
次模型求解后由步骤(5)生成的不满足暂态电压安全约束的故障组成的集合 作为更新
后的SFlt。
到各故障下各节点的暂态电压跌落程度评估指标 通过对各 进行判定,生成优化
后不满足暂态电压安全约束的故障组成的集合为 如下式:
(4),将步骤(5)得到的 作为更新后的SFlt,继续下一轮求解。