CN113765075A - 一种用于电力系统功角失稳的主动解列控制方法及系统 - Google Patents

一种用于电力系统功角失稳的主动解列控制方法及系统 Download PDF

Info

Publication number
CN113765075A
CN113765075A CN202110877923.5A CN202110877923A CN113765075A CN 113765075 A CN113765075 A CN 113765075A CN 202110877923 A CN202110877923 A CN 202110877923A CN 113765075 A CN113765075 A CN 113765075A
Authority
CN
China
Prior art keywords
branch
power system
sbttc
power
index
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
CN202110877923.5A
Other languages
English (en)
Other versions
CN113765075B (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.)
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
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 China Electric Power Research Institute Co Ltd CEPRI filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN202110877923.5A priority Critical patent/CN113765075B/zh
Publication of CN113765075A publication Critical patent/CN113765075A/zh
Application granted granted Critical
Publication of CN113765075B publication Critical patent/CN113765075B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02HEMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
    • H02H7/00Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
    • H02H7/26Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/241The oscillation concerning frequency

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种用于电力系统功角失稳的主动解列控制方法及系统,属于电力系统安全稳定分析技术领域。本发明方法,包括:针对电力系统的功角失稳情况,设置电力系统的控制参数;根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路;将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,判定是否能够匹配到割集;从多个割集中筛选出唯一的待解列割集,对待解列割集内的所有支路同时实施解列控制。本发明在系统发生功角失稳时,可以通过对搜索得到的解列割集支路实施主动解列控制,快速将失步电网解列,从而避免机组功角摇摆引起网络节点电压降低所造成的连锁故障,缩小失稳影响范围。

Description

一种用于电力系统功角失稳的主动解列控制方法及系统
技术领域
本发明涉及电力系统安全稳定分析技术领域,并且更具体地,涉及一种用于电力系统功角失稳的主动解列控制方法及系统。
背景技术
受大容量特高压直流投运工程、大型能源基地集中开发与并网外送工程、偏远地区长距离弱互联工程的影响,我国电网功角稳定问题十分突出。随着电网特性日趋复杂,预先设定的电网运行方式和功角失稳情况下的控制措施容易出现与实际电网运行方式、实际电网发生功角失稳情况下所需控制措施不匹配的问题。在这种情况下,切机、切负荷等控制措施将无法使系统恢复功角稳定,需要及时采取解列控制,断开部分线路以降低功角失稳影响范围、减小功角失稳造成损失。
发明内容
针对上述问题,本发明提出了一种用于电力系统功角失稳的主动解列控制方法,包括:
针对电力系统的功角失稳情况,设置电力系统的控制参数,所述控制参数包括:为避免在电力系统功角失稳进行故障清除后,初始阶段电气量波动影响设置的延时时间TFS、用于判定系统的功角稳定性是否已经恶化的门槛值εasth
在电力系统发生功角失稳后,判断故障清除后的延时时间是否大于TFS,若大于,量测电力系统中支路参数,根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路;
若关键支路k的sBTTC指数值小于设定门槛值εLth,且关键支路k满足预设判据,将关键支路k放入待考察支路集合,将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,判定是否能够匹配到割集;
若匹配到多个可解列割集,从多个割集中筛选出唯一的待解列割集,对待解列割集内的所有支路同时实施解列控制。
可选的,量测电力系统中支路参数,根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路,具体包括:
量测电力系统每一条支路两端节点的电压幅值Umi和Uni和电压相位Δθi,i为电力系统内支路总数,并根据下式计算每一条支路的sBTTC指数;
Figure BDA0003189201230000021
根据如下公式:
Figure BDA0003189201230000022
其中,Δθi=|Δθmi-Δθni|,m、n为示一条线路上取的两个量测点的位置,i为第i条支路;
对每条支路的sBTTC指数从小到大进行排序,定义具有最小sBTTC 指数值的支路为关键支路,并将关键支路的支路编号记为k。
可选的,根据sBTTC指数定位关键支路后,定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
可选的,控制参数还包括:电力系统中待考察支路总个数Nb、用于判定系统失稳的严重程度的门槛值Uvth
可选的,关键支路k满足预设判据,具体包括:
当电力系统发生功角失稳时,电力系统内支路两端节点的电压相位差会不断增大,关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小,即有功功率满足如下:
Figure BDA0003189201230000031
且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、 n的电压相位差Δθ增大而不断增大,即频率偏差Δf满足如下:
Figure BDA0003189201230000032
定义支路垂足电压Uv的计算公式如下:
Figure BDA0003189201230000033
其中,umy=|Im(Um)|、unx=|Re(Un)|、umx=|Re(Um)|,即,umx、umy分别是Um实部、虚部的绝对值,unx是Un实部的绝对值,Um为关键支路k端点 m的电压,Un为关键支路k端点n的电压;
为判断支路垂足电压Uv是否位于支路之上,则ξv满足:
Figure BDA0003189201230000034
其中,ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、 Um节点m与Uv的相位差、Un节点n与Uv的相位差,ΔUmn、ΔUmv、ΔUnv根据下式计算:
Figure DEST_PATH_1
Figure BDA0003189201230000036
关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据,如下:
Uvk<Uvth
可选的,若将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,不能匹配到可解列割集,判定是否满足kq<Nb
其中,kq为提取的具有次小sBTTC指数的支路累计个数,其初值设置为1,Nb为待考察支路总个数,最小值为1,最大值为系统内全部支路的总个数;
若满足,令kq=kq+1再次进行匹配查询;
若不满足,判定全部待考察支路中是否能够有构成割集的支路组合,若有则再次进行匹配查询,若无,则将故障清除后的延时时间加上Δt时间后,再次判断故障清除后的延时时间是否大于TFS,直到故障清除后的延时时间是大于TFS
可选的,从多个割集中筛选出唯一的待解列割集,计算公式如下:
根据下式:
Figure BDA0003189201230000041
其中,asTSC功角失稳目标割集,Pscl0为支路l的稳态有功功率,Nsc为可选解列割集个数。
可选的,若关键支路k不满足预设判据,关键支路k的sBTTC指数值大于设定门槛值εLth,或故障清除后的延时时间小于TFS,将故障清除后的延时时间加上Δt时间后,再次判断故障清除后的延时时间是否大于TFS,直到故障清除后的延时时间是大于TFS
本发明还提出了一种用于电力系统功角失稳的主动解列控制系统,包括:
参数确定模块,针对电力系统的功角失稳情况,设置电力系统的控制参数,所述控制参数包括:为避免在电力系统功角失稳进行故障清除后,初始阶段电气量波动影响设置的延时时间TFS、用于判定系统的功角稳定性是否已经恶化的门槛值εasth
定位模块,在电力系统发生功角失稳后,判断故障清除后的延时时间是否大于TFS,若大于,量测电力系统中支路参数,根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路;
判定模块,在关键支路k的sBTTC指数值小于设定门槛值εLth时,且关键支路k满足预设判据,将关键支路k放入待考察支路集合,将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,判定是否能够匹配到割集;
控制模块,在匹配到多个可解列割集后,从多个割集中筛选出唯一的待解列割集,对待解列割集内的所有支路同时实施解列控制。
可选的,量测电力系统中支路参数,根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路,具体包括:
量测电力系统每一条支路两端节点的电压幅值Umi和Uni和电压相位Δθi,i为电力系统内支路总数,并根据下式计算每一条支路的sBTTC指数;
Figure BDA0003189201230000051
根据如下公式:
Figure BDA0003189201230000052
其中,Δθi=|Δθmi-Δθni|,m、n为示一条线路上取的两个量测点的位置,i为第i条支路;
对每条支路的sBTTC指数从小到大进行排序,定义具有最小sBTTC 指数值的支路为关键支路,并将关键支路的支路编号记为k。
可选的,根据sBTTC指数定位关键支路后,定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
可选的,控制参数还包括:电力系统中待考察支路总个数Nb、用于判定系统失稳的严重程度的门槛值vth
可选的,关键支路k满足预设判据,具体包括:
当电力系统发生功角失稳时,电力系统内支路两端节点的电压相位差会不断增大,关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小,即有功功率满足如下:
Figure BDA0003189201230000053
且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、 n的电压相位差Δθ增大而不断增大,即频率偏差Δf满足如下:
Figure BDA0003189201230000054
定义支路垂足电压Uv的计算公式如下:
Figure BDA0003189201230000055
其中,umy=|Im(Um)|、unx=|Re(Un)|、umx=|Re(Um)|,即,umx、umy分别是Um实部、虚部的绝对值,unx是Un实部的绝对值,Um为关键支路k端点 m的电压,Un为关键支路k端点n的电压;
为判断支路垂足电压Uv是否位于支路之上,则ξv满足:
Figure BDA0003189201230000056
其中,ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、Um节点m与Uv的相位差、Un节点n与Uv的相位差,ΔUmn、ΔUmv、ΔUnv根据下式计算:
Figure 241146DEST_PATH_1
Figure BDA0003189201230000062
关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据,如下:
Uvk<Uvth
可选的,若将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,不能匹配到可解列割集,判定是否满足kq<Nb
其中,kq为提取的具有次小sBTTC指数的支路累计个数,其初值设置为1,Nb为待考察支路总个数,最小值为1,最大值为系统内全部支路的总个数;
若满足,令kq=kq+1再次进行匹配查询;
若不满足,判定全部待考察支路中是否能够有构成割集的支路组合,若有则再次进行匹配查询,若无,则将故障清除后的延时时间加上Δt时间后,再次判断故障清除后的延时时间是否大于TFS,直到故障清除后的延时时间是大于TFS
可选的,从多个割集中筛选出唯一的待解列割集,计算公式如下:
根据下式:
Figure BDA0003189201230000063
其中,asTSC功角失稳目标割集,Pscl0为支路l的稳态有功功率,Nsc为可选解列割集个数。
可选的,若关键支路k不满足预设判据,关键支路k的sBTTC指数值大于设定门槛值εLth,或故障清除后的延时时间小于TFS,将故障清除后的延时时间加上Δt时间后,再次判断故障清除后的延时时间是否大于TFS,直到故障清除后的延时时间是大于TFS
本发明在系统发生功角失稳时,可以通过对搜索得到的解列割集支路实施主动解列控制,快速将失步电网解列,从而避免机组功角摇摆引起网络节点电压降低所造成的连锁故障,缩小失稳影响范围。
附图说明
图1为本发明方法的流程图;
图2为本发明稳定分析的典型系统模型;
图3为本发明主动解列控制流程原理图;
图4为本发明系统的结构图。
具体实施方式
现在参考附图介绍本发明的示例性实施方式,然而,本发明可以用许多不同的形式来实施,并且不局限于此处描述的实施例,提供这些实施例是为了详尽地且完全地公开本发明,并且向所属技术领域的技术人员充分传达本发明的范围。对于表示在附图中的示例性实施方式中的术语并不是对本发明的限定。在附图中,相同的单元/元件使用相同的附图标记。
除非另有说明,此处使用的术语(包括科技术语)对所属技术领域的技术人员具有通常的理解含义。另外,可以理解的是,以通常使用的词典限定的术语,应当被理解为与其相关领域的语境具有一致的含义,而不应该被理解为理想化的或过于正式的意义。
本发明还提出了一种用于电力系统功角失稳的主动解列控制方法,如图1所示,包括:
针对电力系统的功角失稳情况,设置电力系统的控制参数,所述控制参数包括:为避免在电力系统功角失稳进行故障清除后,初始阶段电气量波动影响设置的延时时间TFS、用于判定系统的功角稳定性是否已经恶化的门槛值εasth
在电力系统发生功角失稳后,判断故障清除后的延时时间是否大于TFS,若大于,量测电力系统中支路参数,根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路;
若关键支路k的sBTTC指数值小于设定门槛值εLth,且关键支路k满足预设判据,将关键支路k放入待考察支路集合,将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,判定是否能够匹配到割集;
若匹配到多个可解列割集,从多个割集中筛选出唯一的待解列割集,对待解列割集内的所有支路同时实施解列控制。
其中,量测电力系统中支路参数,根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路,具体包括:
量测电力系统每一条支路两端节点的电压幅值Umi和Uni和电压相位Δθi,i为电力系统内支路总数,并根据下式计算每一条支路的sBTTC指数;
Figure BDA0003189201230000081
根据如下公式:
Figure BDA0003189201230000082
其中,Δθi=|Δθmi-Δθni|,m、n为示一条线路上取的两个量测点的位置,i为第i条支路;
对每条支路的sBTTC指数从小到大进行排序,定义具有最小sBTTC 指数值的支路为关键支路,并将关键支路的支路编号记为k。
其中,根据sBTTC指数定位关键支路后,定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
其中,控制参数还包括:电力系统中待考察支路总个数Nb、用于判定系统失稳的严重程度的门槛值Uvth
其中,关键支路k满足预设判据,具体包括:
当电力系统发生功角失稳时,电力系统内支路两端节点的电压相位差会不断增大,关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小,即有功功率满足如下:
Figure BDA0003189201230000083
且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、 n的电压相位差Δθ增大而不断增大,即频率偏差Δf满足如下:
Figure BDA0003189201230000091
定义支路垂足电压Uv的计算公式如下:
Figure BDA0003189201230000092
其中,umy=|Im(Um)|、unx=|Re(Un)|、umx=|Re(Um)|,即,umx、umy分别是Um实部、虚部的绝对值,unx是Un实部的绝对值,Um为关键支路k端点 m的电压,Un为关键支路k端点n的电压;
为判断支路垂足电压Uv是否位于支路之上,则ξv满足:
Figure BDA0003189201230000093
其中,ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、 Um节点m与Uv的相位差、Un节点n与Uv的相位差,ΔUmn、ΔUmv、ΔUnv根据下式计算:
Figure 33653DEST_PATH_1
Figure BDA0003189201230000095
关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据,如下:
Uvk<Uvth
其中,若将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,不能匹配到可解列割集,判定是否满足kq<Nb
其中,kq为提取的具有次小sBTTC指数的支路累计个数,其初值设置为1,Nb为待考察支路总个数,最小值为1,最大值为系统内全部支路的总个数;
若满足,令kq=kq+1再次进行匹配查询;
若不满足,判定全部待考察支路中是否能够有构成割集的支路组合,若有则再次进行匹配查询,若无,则将故障清除后的延时时间加上Δt时间后,再次判断故障清除后的延时时间是否大于TFS,直到故障清除后的延时时间是大于TFS
其中,从多个割集中筛选出唯一的待解列割集,计算公式如下:
根据下式:
Figure BDA0003189201230000096
其中,asTSC功角失稳目标割集,Pscl0为支路l的稳态有功功率,Nsc为可选解列割集个数。
其中,若关键支路k不满足预设判据,关键支路k的sBTTC指数值大于设定门槛值εLth,或故障清除后的延时时间小于TFS,将故障清除后的延时时间加上Δt时间后,再次判断故障清除后的延时时间是否大于TFS,直到故障清除后的延时时间是大于TFS
下面结合具体实施例进行说明:
本发明根据功角稳定分析的典型系统模型进行分析,模型如图2所示,流程原理如图3所示,包括如下:
步骤一:设置控制相关参数TFS、εasth、Nb、Uvth
TFS为避免故障清除后的初始阶段电气量波动影响而人为设置的延时时间;
εasth为人为设定的一个门槛值,用于判定测试系统的功角稳定性是否已经恶化;
Nb为待考察支路总个数,对于不同的测试系统,通过对测试系统内全部支路进行计数得到该参数;
Uvth为人为设定的一个门槛值,用于判定测试系统失稳的严重程度。
步骤二:判定故障清除后的持续时间是否大于TFS
TFS为避免故障清除后的初始阶段电气量波动影响而设置的延时时间。
若故障清除后的持续时间大于TFS,则进入步骤三;
若故障清除后的持续时间小于TFS,则使t=t+Δt,并返回步骤二。
步骤三:量测每一条支路两端节点的电压幅值Umi、Uni和电压相位Δθi (i=1、2…N,N为系统内支路的总数),并根据式(1)计算每一条支路的 sBTTC指数。按照式(2)排序并定位其中的关键支路。
Δθi=|Δθmi-Δθni|,表示支路i两端节点的电压相位差。Umi、Uni、Δθmi、Δθni在图2中均已标出,下标m、n表示一条线路上取的两个量测点的位置,通常取线路两端。Umi是支路i上节点m的电压幅值,Uni是支路i上节点n的电压幅值,Δθmi是支路i上节点m的电压相位,Δθni是支路i上节点n的电压相位。
Figure BDA0003189201230000111
根据式(2)对每条支路sBTTC指数从小到大进行排序,定义具有最小sBTTC指数值的支路为关键支路,其支路编号记为k,定义具有次小 sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
Figure BDA0003189201230000112
步骤四:判定关键支路k是否满足式(3)。
由于sBTTC指数具有随功角稳定性恶化而单调递减的特性,当关键支路k的sBTTC指数值小于设定门槛值εasth时,表明系统的功角稳定性已显著恶化,系统已经发生功角失稳。
sBTTCkasth (3)
若满足,则进入步骤五;
若不满足,使t=t+Δt,并返回步骤二。
步骤五:判定关键支路k是否同时满足式(4)、(5)、(7)和(10)。
当系统发生功角失稳时,支路两端节点的电压相位差会不断增大。关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小,即满足式(4)的描述;且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、n的电压相位差Δθ增大而不断增大,即满足式(5)的描述。因此,可以将支路所呈现出的式(4)和式(5)特征作为功角失稳的判据之一,以提高本流程判断的准确性。
Figure BDA0003189201230000113
Figure BDA0003189201230000114
定义支路垂足电压Uv的计算公式如下:
Figure BDA0003189201230000115
其中,umy=|Im(Um)|、unx=|Re(Un)|、umx=|Re(Um)|,也即, umx、umy分别是Um实部、虚部的绝对值,unx是Un实部的绝对值。
为判断支路垂足电压Uv是否位于支路之上,定义位置系数ξv如式(7) 所示。若Uv位于支路之上,则必有ΔUnv<ΔUmn且ΔUmv<ΔUmn,即ξv大于1且小于2;若Uv位于支路之外,则必有ΔUnv>ΔUmn或ΔUmv>ΔUmn,对应ξv小于1。
Figure BDA0003189201230000121
其中,ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、 Um节点m与Uv的相位差、Un节点n与Uv的相位差,ΔUmn、ΔUmv、ΔUnv根据式(7)和(8)计算:
Figure BDA0003189201230000122
Figure BDA0003189201230000123
对于复杂的包含多台发电机的电力系统,还可增加关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据,即以式(9)判断系统扰动失稳的严重程度。满足式(9)即意味着系统已临近功角失稳或已经发生功角失稳,可以进行实施主动解列控制的决策。
Uvk<Uvth (10)
若满足,则进入步骤六;
若不满足,使t=t+Δt,并返回步骤二。
步骤六:根据式(2),提取一条具有次小sBTTC指数且满足式(4)、 (5)、(7)和(10)的支路,放入待考察支路集合,将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,判定是否能够匹配到割集。
若能够匹配到割集,则进入步骤七;
若不能匹配到割集,则进入步骤八。
步骤七:判定匹配到的割集是否单一。
若匹配到多个割集,则进入步骤十;
若匹配到单一割集,则进入步骤十一。
步骤八:判定是否满足kq<Nb
定义kq为提取的具有次小sBTTC指数的支路累计个数,其初值设置为1。定义Nb为待考察支路总个数,该值由人为设定,最小值为1,最大值为系统内全部支路的总个数。
若满足,令kq=kq+1并返回步骤六;
若不满足,说明提取的具有次小sBTTC指数的支路累计个数已达到最大值,也即已经遍历了全部待考察支路,进入步骤九。
步骤九:判定全部待考察支路中是否能够有构成割集的支路组合。
若有能够有构成割集的支路组合,进入步骤七;
若没有能够有构成割集的支路组合,使t=t+Δt,并返回步骤二;
步骤十:根据式(11),从搜索匹配到的多个割集中筛选出唯一的待解列割集。
当搜索匹配到多个可解列割集,则以割集中全部支路稳态有功功率之和最小作为优化选择原则并确定功角失稳目标解列割集(asTSC,angle stability target splittingcutset)。
式中,Lsc为割集SC中所包含的支路总数、Pscl0为割集SC中支路l的稳态有功功率、Pscl0为可解列割集的个数。
Figure BDA0003189201230000131
步骤十一:对待解列割集内的所有支路同时实施解列控制,并结束本流程。
本发明还提出了一种用于电力系统功角失稳的主动解列控制系统200,如图4所示,包括:
参数确定模块201,针对电力系统的功角失稳情况,设置电力系统的控制参数,所述控制参数包括:为避免在电力系统功角失稳进行故障清除后,初始阶段电气量波动影响设置的延时时间TFS、用于判定系统的功角稳定性是否已经恶化的门槛值εasth
定位模块202,在电力系统发生功角失稳后,判断故障清除后的延时时间是否大于TFS,若大于,量测电力系统中支路参数,根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路;
判定模块203,在关键支路k的sBTTC指数值小于设定门槛值εLth时,且关键支路k满足预设判据,将关键支路k放入待考察支路集合,将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,判定是否能够匹配到割集;
控制模块204,在匹配到多个可解列割集后,从多个割集中筛选出唯一的待解列割集,对待解列割集内的所有支路同时实施解列控制。
其中,量测电力系统中支路参数,根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路,具体包括:
量测电力系统每一条支路两端节点的电压幅值Umi和Uni和电压相位Δθi,i为电力系统内支路总数,并根据下式计算每一条支路的sBTTC指数;
Figure BDA0003189201230000141
根据如下公式:
Figure BDA0003189201230000142
其中,Δθi=|Δθmi-Δθmi|,m、n为示一条线路上取的两个量测点的位置,i为第i条支路;
对每条支路的sBTTC指数从小到大进行排序,定义具有最小sBTTC 指数值的支路为关键支路,并将关键支路的支路编号记为k。
其中,根据sBTTC指数定位关键支路后,定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
其中,控制参数还包括:电力系统中待考察支路总个数Nb、用于判定系统失稳的严重程度的门槛值Uvth
其中,关键支路k满足预设判据,具体包括:
当电力系统发生功角失稳时,电力系统内支路两端节点的电压相位差会不断增大,关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小,即有功功率满足如下:
Figure BDA0003189201230000143
且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、 n的电压相位差Δθ增大而不断增大,即频率偏差Δf满足如下:
Figure BDA0003189201230000144
定义支路垂足电压Uv的计算公式如下:
Figure BDA0003189201230000145
其中,umy=|Im(Um)|、unx=|Re(Un)|、umx=|Re(Um)|,即,umx、umy分别是Um实部、虚部的绝对值,unx是Un实部的绝对值,Um为关键支路k端点 m的电压,Un为关键支路k端点n的电压;
为判断支路垂足电压Uv是否位于支路之上,则ξv满足:
Figure BDA0003189201230000151
其中,ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、 Um节点m与Uv的相位差、Un节点n与Uv的相位差,ΔUmn、ΔUmv、ΔUnv根据下式计算:
Figure 875707DEST_PATH_1
Figure BDA0003189201230000153
关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据,如下:
Uvk<Uvth
其中,若将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,不能匹配到可解列割集,判定是否满足kq<Nb
其中,kq为提取的具有次小sBTTC指数的支路累计个数,其初值设置为1,Nb为待考察支路总个数,最小值为1,最大值为系统内全部支路的总个数;
若满足,令kq=kq+1再次进行匹配查询;
若不满足,判定全部待考察支路中是否能够有构成割集的支路组合,若有则再次进行匹配查询,若无,则将故障清除后的延时时间加上Δt时间后,再次判断故障清除后的延时时间是否大于TFS,直到故障清除后的延时时间是大于TFS
其中,从多个割集中筛选出唯一的待解列割集,计算公式如下:
根据下式:
Figure BDA0003189201230000154
其中,asTSC功角失稳目标割集,Pscl0为支路l的稳态有功功率,Nsc为可选解列割集个数。
其中,若关键支路k不满足预设判据,关键支路k的sBTTC指数值大于设定门槛值εLth,或故障清除后的延时时间小于TFS,将故障清除后的延时时间加上Δt时间后,再次判断故障清除后的延时时间是否大于TFS,直到故障清除后的延时时间是大于TFS
本发明在系统发生功角失稳时,可以通过对搜索得到的解列割集支路实施主动解列控制,快速将失步电网解列,从而避免机组功角摇摆引起网络节点电压降低所造成的连锁故障,缩小失稳影响范围。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本发明实施例中的方案可以采用各种计算机语言实现,例如,面向对象的程序设计语言Java和直译式脚本语言JavaScript等。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (16)

1.一种用于电力系统功角失稳的主动解列控制方法,所述方法包括:
针对电力系统的功角失稳情况,设置电力系统的控制参数,所述控制参数包括:为避免在电力系统功角失稳进行故障清除后,初始阶段电气量波动影响设置的延时时间TFS、用于判定系统的功角稳定性是否已经恶化的门槛值εasth
在电力系统发生功角失稳后,判断故障清除后的延时时间是否大于TFS,若大于,量测电力系统中支路参数,根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路;
若关键支路k的sBTTC指数值小于设定门槛值εlth,且关键支路k满足预设判据,将关键支路k放入待考察支路集合,将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,判定是否能够匹配到割集;
若匹配到多个可解列割集,从多个割集中筛选出唯一的待解列割集,对待解列割集内的所有支路同时实施解列控制。
2.根据权利要求1所述的方法,所述量测电力系统中支路参数,根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路,具体包括:
量测电力系统每一条支路两端节点的电压幅值Umi和Uni和电压相位Δθi,i为电力系统内支路总数,并根据下式计算每一条支路的sBTTC指数;
Figure FDA0003189201220000011
根据如下公式:
Figure FDA0003189201220000012
其中,Δθi=|Δθmi-Δθni|,m、n为示一条线路上取的两个量测点的位置,i为第i条支路;
对每条支路的sBTTC指数从小到大进行排序,定义具有最小sBTTC 指数值的支路为关键支路,并将关键支路的支路编号记为k。
3.根据权利要求2所述的方法,所述根据sBTTC指数定位关键支路后,定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
4.根据权利要求1所述的方法,所述控制参数还包括:电力系统中待考察支路总个数Nb、用于判定系统失稳的严重程度的门槛值Uvth
5.根据权利要求1所述的方法,所述关键支路k满足预设判据,具体包括:
当电力系统发生功角失稳时,电力系统内支路两端节点的电压相位差会不断增大,关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小,即有功功率满足如下:
Figure FDA0003189201220000021
且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、n的电压相位差Δθ增大而不断增大,即频率偏差Δf满足如下:
Figure FDA0003189201220000022
定义支路垂足电压Uv的计算公式如下:
Figure FDA0003189201220000023
其中,umy=|Im(Um)|、unx=|Re(Un)|、umx=|Re(Um)|,即,umx、umy分别是Um实部、虚部的绝对值,unx是Un实部的绝对值,Um为关键支路k端点m的电压,Un为关键支路k端点n的电压;
为判断支路垂足电压Uv是否位于支路之上,则ζv满足:
Figure FDA0003189201220000024
其中,ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、Um节点m与Uv的相位差、Un节点n与Uv的相位差,ΔUmn、ΔUmv、ΔUnv根据下式计算:
Figure 1
Figure FDA0003189201220000026
关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据,如下:
Uvk<Uvth
6.根据权利要求1所述的方法,所述若将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,不能匹配到可解列割集,判定是否满足kq<Nb
其中,kq为提取的具有次小sBTTC指数的支路累计个数,其初值设置为1,Nb为待考察支路总个数,最小值为1,最大值为系统内全部支路的总个数;
若满足,令kq=kq+1再次进行匹配查询;
若不满足,判定全部待考察支路中是否能够有构成割集的支路组合,若有则再次进行匹配查询,若无,则将故障清除后的延时时间加上Δt时间后,再次判断故障清除后的延时时间是否大于TFS,直到故障清除后的延时时间是大于TFS
7.根据权利要求1所述的方法,所述从多个割集中筛选出唯一的待解列割集,计算公式如下:
根据下式:
Figure FDA0003189201220000031
其中,asTSC功角失稳目标割集,Pscl0为支路1的稳态有功功率,Nsc为可选解列割集个数。
8.根据权利要求1所述的方法,所述若关键支路k不满足预设判据,关键支路k的sBTTC指数值大于设定门槛值εLth,或故障清除后的延时时间小于TFS,将故障清除后的延时时间加上Δt时间后,再次判断故障清除后的延时时间是否大于TFS,直到故障清除后的延时时间是大于TFS
9.一种用于电力系统功角失稳的主动解列控制系统,所述系统包括:
参数确定模块,针对电力系统的功角失稳情况,设置电力系统的控制参数,所述控制参数包括:为避免在电力系统功角失稳进行故障清除后,初始阶段电气量波动影响设置的延时时间TFS、用于判定系统的功角稳定性是否已经恶化的门槛值εasth
定位模块,在电力系统发生功角失稳后,判断故障清除后的延时时间是否大于TFS,若大于,量测电力系统中支路参数,根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路;
判定模块,在关键支路k的sBTTC指数值小于设定门槛值εLth时,且关键支路k满足预设判据,将关键支路k放入待考察支路集合,将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,判定是否能够匹配到割集;
控制模块,在匹配到多个可解列割集后,从多个割集中筛选出唯一的待解列割集,对待解列割集内的所有支路同时实施解列控制。
10.根据权利要求9所述的系统,所述量测电力系统中支路参数,根据支路参数确定电力系统中每一条支路的sBTTC指数,根据sBTTC指数定位关键支路,具体包括:
量测电力系统每一条支路两端节点的电压幅值Umi和Uni和电压相位Δθi,i为电力系统内支路总数,并根据下式计算每一条支路的sBTTC指数;
Figure FDA0003189201220000041
根据如下公式:
Figure FDA0003189201220000042
其中,Δθi=|Δθmi-Δθni|,m、n为示一条线路上取的两个量测点的位置,i为第i条支路;
对每条支路的sBTTC指数从小到大进行排序,定义具有最小sBTTC指数值的支路为关键支路,并将关键支路的支路编号记为k。
11.根据权利要求10所述的系统,所述根据sBTTC指数定位关键支路后,定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
12.根据权利要求9所述的系统,所述控制参数还包括:电力系统中待考察支路总个数Nb、用于判定系统失稳的严重程度的门槛值Uvth
13.根据权利要求9所述的系统,所述关键支路k满足预设判据,具体包括:
当电力系统发生功角失稳时,电力系统内支路两端节点的电压相位差会不断增大,关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小,即有功功率满足如下:
Figure FDA0003189201220000051
且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、n的电压相位差Δθ增大而不断增大,即频率偏差Δf满足如下:
Figure FDA0003189201220000052
定义支路垂足电压Uv的计算公式如下:
Figure FDA0003189201220000053
其中,umy=|Im(Um)|、unx=|Re(Un)|、umx=|Re(Um)|,即,umx、umy分别是Um实部、虚部的绝对值,unx是Un实部的绝对值,Um为关键支路k端点m的电压,Un为关键支路k端点n的电压;
为判断支路垂足电压Uv是否位于支路之上,则ξv满足:
Figure FDA0003189201220000054
其中,ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、Um节点m与Uv的相位差、Un节点n与Uv的相位差,ΔUmn、ΔUmv、ΔUnv根据下式计算:
Figure 1
Figure FDA0003189201220000056
关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据,如下:
Uvk<Uvth
14.根据权利要求9所述的系统,所述若将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询,不能匹配到可解列割集,判定是否满足kq<Nb
其中,kq为提取的具有次小sBTTC指数的支路累计个数,其初值设置为1,Nb为待考察支路总个数,最小值为1,最大值为系统内全部支路的总个数;
若满足,令kq=kq+1再次进行匹配查询;
若不满足,判定全部待考察支路中是否能够有构成割集的支路组合,若有则再次进行匹配查询,若无,则将故障清除后的延时时间加上Δt时间后,再次判断故障清除后的延时时间是否大于TFS,直到故障清除后的延时时间是大于TFS
15.根据权利要求9所述的系统,所述从多个割集中筛选出唯一的待解列割集,计算公式如下:
根据下式:
Figure FDA0003189201220000061
其中,asTSC功角失稳目标割集,Pscl0为支路1的稳态有功功率,Nsc为可选解列割集个数。
16.根据权利要求9所述的系统,所述若关键支路k不满足预设判据,关键支路k的sBTTC指数值大于设定门槛值εLth,或故障清除后的延时时间小于TFS,将故障清除后的延时时间加上Δt时间后,再次判断故障清除后的延时时间是否大于TFS,直到故障清除后的延时时间是大于TFS
CN202110877923.5A 2021-07-30 2021-07-30 一种用于电力系统功角失稳的主动解列控制方法及系统 Active CN113765075B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110877923.5A CN113765075B (zh) 2021-07-30 2021-07-30 一种用于电力系统功角失稳的主动解列控制方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110877923.5A CN113765075B (zh) 2021-07-30 2021-07-30 一种用于电力系统功角失稳的主动解列控制方法及系统

Publications (2)

Publication Number Publication Date
CN113765075A true CN113765075A (zh) 2021-12-07
CN113765075B CN113765075B (zh) 2023-11-03

Family

ID=78788271

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110877923.5A Active CN113765075B (zh) 2021-07-30 2021-07-30 一种用于电力系统功角失稳的主动解列控制方法及系统

Country Status (1)

Country Link
CN (1) CN113765075B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108649559A (zh) * 2018-04-25 2018-10-12 中国电力科学研究院有限公司 一种评估电力系统紧急控制策略的效果的方法和系统
CN109428327A (zh) * 2017-09-01 2019-03-05 中国电力科学研究院 基于响应的电网关键支路与主导稳定模式识别方法及系统
CN112332417A (zh) * 2021-01-07 2021-02-05 中国电力科学研究院有限公司 基于校正振荡中心电压的广域协同失步解列方法及装置
CN112421579A (zh) * 2020-09-21 2021-02-26 中国电力科学研究院有限公司 一种基于低频低压的主动解列控制方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109428327A (zh) * 2017-09-01 2019-03-05 中国电力科学研究院 基于响应的电网关键支路与主导稳定模式识别方法及系统
CN108649559A (zh) * 2018-04-25 2018-10-12 中国电力科学研究院有限公司 一种评估电力系统紧急控制策略的效果的方法和系统
CN112421579A (zh) * 2020-09-21 2021-02-26 中国电力科学研究院有限公司 一种基于低频低压的主动解列控制方法及系统
CN112332417A (zh) * 2021-01-07 2021-02-05 中国电力科学研究院有限公司 基于校正振荡中心电压的广域协同失步解列方法及装置

Also Published As

Publication number Publication date
CN113765075B (zh) 2023-11-03

Similar Documents

Publication Publication Date Title
CN109118098B (zh) 高比例风电并网的连锁故障风险评估方法及系统
US20140309952A1 (en) Method for Analyzing Faults in Ungrounded Power Distribution Systems
CN103928914B (zh) 继电保护设备整定方法及装置
CN106548265B (zh) 一种基于连锁故障事故链搜索的输电网可靠性评估方法
CN109375050B (zh) 一种基于初始故障线路割集的电力系统关键断面辨识方法
CN115469259B (zh) 基于rbf神经网络的ct误差状态在线定量评估方法及装置
CN108921725A (zh) 一种复杂配电网可靠性指标快速解析计算方法
Wang et al. A DRL-aided multi-layer stability model calibration platform considering multiple events
CN112684295A (zh) 一种基于相似离度的高渗透率下配电网故障选线方法及系统
CN109873418B (zh) 一种暂态功角稳定模式的强相关输电断面识别方法及系统
CN112649695A (zh) 一种基于节点全覆盖的配电网电能质量评估方案
CN109066651A (zh) 风电-负荷场景的极限传输功率的计算方法
CN112989587B (zh) 一种电容式电压互感器劣化成因在线分析方法及系统
Wang et al. Neural networks based lyapunov functions for transient stability analysis and assessment of power systems
CN113765075A (zh) 一种用于电力系统功角失稳的主动解列控制方法及系统
JP7077250B2 (ja) 電力系統安定化システム
CN104112240A (zh) 一种基于拓扑岛分析的停电分析方法
CN109167353B (zh) 暂态功角稳定视角下的电网强关联断面识别方法
Bernardes et al. Topological processing of mutually coupled circuits for directional overcurrent protection
CN108448569B (zh) 一种大电网机电暂态仿真数值收敛性的识别和提高方法
CN113824094B (zh) 一种用于电力系统电压失稳的主动解列控制方法及系统
CN109038544B (zh) 一种基于电压轨迹特征的暂态稳定识别方法
CN110994564A (zh) 配电网保护投入开关自动选取和定值整定方法
CN112713584B (zh) 适用于稳控系统预想故障集外故障的控制方法和装置
CN113363987B (zh) 一种配电网主从式故障自愈控制方法及系统

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