CN113765075B - 一种用于电力系统功角失稳的主动解列控制方法及系统 - Google Patents
一种用于电力系统功角失稳的主动解列控制方法及系统 Download PDFInfo
- Publication number
- CN113765075B CN113765075B CN202110877923.5A CN202110877923A CN113765075B CN 113765075 B CN113765075 B CN 113765075B CN 202110877923 A CN202110877923 A CN 202110877923A CN 113765075 B CN113765075 B CN 113765075B
- Authority
- CN
- China
- Prior art keywords
- branch
- sbttc
- power system
- index
- power
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000005520 cutting process Methods 0.000 claims abstract description 25
- 238000012216 screening Methods 0.000 claims abstract description 6
- 238000000926 separation method Methods 0.000 claims abstract description 3
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 230000001687 destabilization Effects 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 4
- 230000009467 reduction Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 3
- 230000004075 alteration Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
- H02H7/00—Emergency 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/26—Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The 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指数；
根据如下公式：
其中，Δθi＝|Δθmi-Δθni|，m、n为示一条线路上取的两个量测点的位置，i为第i条支路；
对每条支路的sBTTC指数从小到大进行排序，定义具有最小sBTTC指数值的支路为关键支路，并将关键支路的支路编号记为k。
可选的，根据sBTTC指数定位关键支路后，定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
可选的，控制参数还包括：电力系统中待考察支路总个数Nb、用于判定系统失稳的严重程度的门槛值Uvth。
可选的，关键支路k满足预设判据，具体包括：
当电力系统发生功角失稳时，电力系统内支路两端节点的电压相位差会不断增大，关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小，即有功功率满足如下：
且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、n的电压相位差Δθ增大而不断增大，即频率偏差Δf满足如下：
定义支路垂足电压Uv的计算公式如下：
其中，umy＝|Im(Um)|、unx＝|Re(Un)|、umx＝|Re(Um)|，即，umx、umy分别是Um实部、虚部的绝对值，unx是Un实部的绝对值，Um为关键支路k端点m的电压，Un为关键支路k端点n的电压；
为判断支路垂足电压Uv是否位于支路之上，则ξv满足：
其中，ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、Um节点m与Uv的相位差、Un节点n与Uv的相位差，ΔUmn、ΔUmv、ΔUnv根据下式计算：
关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据，如下：
Uvk<Uvth。
可选的，若将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询，不能匹配到可解列割集，判定是否满足kq<Nb；
其中，kq为提取的具有次小sBTTC指数的支路累计个数，其初值设置为1，Nb为待考察支路总个数，最小值为1，最大值为系统内全部支路的总个数；
若满足，令kq＝kq+1再次进行匹配查询；
若不满足，判定全部待考察支路中是否能够有构成割集的支路组合，若有则再次进行匹配查询，若无，则将故障清除后的延时时间加上Δt时间后，再次判断故障清除后的延时时间是否大于TFS，直到故障清除后的延时时间是大于TFS。
可选的，从多个割集中筛选出唯一的待解列割集，计算公式如下：
根据下式：
其中，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指数；
根据如下公式：
其中，Δθi＝|Δθmi-Δθni|，m、n为示一条线路上取的两个量测点的位置，i为第i条支路；
对每条支路的sBTTC指数从小到大进行排序，定义具有最小sBTTC指数值的支路为关键支路，并将关键支路的支路编号记为k。
可选的，根据sBTTC指数定位关键支路后，定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
可选的，控制参数还包括：电力系统中待考察支路总个数Nb、用于判定系统失稳的严重程度的门槛值Uvth。
可选的，关键支路k满足预设判据，具体包括：
当电力系统发生功角失稳时，电力系统内支路两端节点的电压相位差会不断增大，关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小，即有功功率满足如下：
且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、n的电压相位差Δθ增大而不断增大，即频率偏差Δf满足如下：
定义支路垂足电压Uv的计算公式如下：
其中，umy＝|Im(Um)|、unx＝|Re(Un)|、umx＝|Re(Um)|，即，umx、umy分别是Um实部、虚部的绝对值，unx是Un实部的绝对值，Um为关键支路k端点m的电压，Un为关键支路k端点n的电压；
为判断支路垂足电压Uv是否位于支路之上，则ξv满足：
其中，ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、Um节点m与Uv的相位差、Un节点n与Uv的相位差，ΔUmn、ΔUmv、ΔUnv根据下式计算：
关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据，如下：
Uvk<Uvth。
可选的，若将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询，不能匹配到可解列割集，判定是否满足kq<Nb；
其中，kq为提取的具有次小sBTTC指数的支路累计个数，其初值设置为1，Nb为待考察支路总个数，最小值为1，最大值为系统内全部支路的总个数；
若满足，令kq＝kq+1再次进行匹配查询；
若不满足，判定全部待考察支路中是否能够有构成割集的支路组合，若有则再次进行匹配查询，若无，则将故障清除后的延时时间加上Δt时间后，再次判断故障清除后的延时时间是否大于TFS，直到故障清除后的延时时间是大于TFS。
可选的，从多个割集中筛选出唯一的待解列割集，计算公式如下：
根据下式：
其中，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指数；
根据如下公式：
其中，Δθi＝|Δθmi-Δθni|，m、n为示一条线路上取的两个量测点的位置，i为第i条支路；
对每条支路的sBTTC指数从小到大进行排序，定义具有最小sBTTC指数值的支路为关键支路，并将关键支路的支路编号记为k。
其中，根据sBTTC指数定位关键支路后，定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
其中，控制参数还包括：电力系统中待考察支路总个数Nb、用于判定系统失稳的严重程度的门槛值Uvth。
其中，关键支路k满足预设判据，具体包括：
当电力系统发生功角失稳时，电力系统内支路两端节点的电压相位差会不断增大，关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小，即有功功率满足如下：
且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、n的电压相位差Δθ增大而不断增大，即频率偏差Δf满足如下：
定义支路垂足电压Uv的计算公式如下：
其中，umy＝|Im(Um)|、unx＝|Re(Un)|、umx＝|Re(Um)|，即，umx、umy分别是Um实部、虚部的绝对值，unx是Un实部的绝对值，Um为关键支路k端点m的电压，Un为关键支路k端点n的电压；
为判断支路垂足电压Uv是否位于支路之上，则ξv满足：
其中，ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、Um节点m与Uv的相位差、Un节点n与Uv的相位差，ΔUmn、ΔUmv、ΔUnv根据下式计算：
关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据，如下：
Uvk<Uvth。
其中，若将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询，不能匹配到可解列割集，判定是否满足kq<Nb；
其中，kq为提取的具有次小sBTTC指数的支路累计个数，其初值设置为1，Nb为待考察支路总个数，最小值为1，最大值为系统内全部支路的总个数；
若满足，令kq＝kq+1再次进行匹配查询；
若不满足，判定全部待考察支路中是否能够有构成割集的支路组合，若有则再次进行匹配查询，若无，则将故障清除后的延时时间加上Δt时间后，再次判断故障清除后的延时时间是否大于TFS，直到故障清除后的延时时间是大于TFS。
其中，从多个割集中筛选出唯一的待解列割集，计算公式如下：
根据下式：
其中，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的电压相位。
根据式(2)对每条支路sBTTC指数从小到大进行排序，定义具有最小sBTTC指数值的支路为关键支路，其支路编号记为k，定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
步骤四：判定关键支路k是否满足式(3)。
由于sBTTC指数具有随功角稳定性恶化而单调递减的特性，当关键支路k的sBTTC指数值小于设定门槛值εasth时，表明系统的功角稳定性已显著恶化，系统已经发生功角失稳。
sBTTCk<εasth (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)特征作为功角失稳的判据之一，以提高本流程判断的准确性。
定义支路垂足电压Uv的计算公式如下：
其中，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。
其中，ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、Um节点m与Uv的相位差、Un节点n与Uv的相位差，ΔUmn、ΔUmv、ΔUnv根据式(7)和(8)计算：
对于复杂的包含多台发电机的电力系统，还可增加关键支路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为可解列割集的个数。
步骤十一：对待解列割集内的所有支路同时实施解列控制，并结束本流程。
本发明还提出了一种用于电力系统功角失稳的主动解列控制系统200，如图4所示，包括：
参数确定模块201，针对电力系统的功角失稳情况，设置电力系统的控制参数，所述控制参数包括：为避免在电力系统功角失稳进行故障清除后，初始阶段电气量波动影响设置的延时时间TFS、用于判定系统的功角稳定性是否已经恶化的门槛值εasth；
定位模块202，在电力系统发生功角失稳后，判断故障清除后的延时时间是否大于TFS，若大于，量测电力系统中支路参数，根据支路参数确定电力系统中每一条支路的sBTTC指数，根据sBTTC指数定位关键支路；
判定模块203，在关键支路k的sBTTC指数值小于设定门槛值εLth时，且关键支路k满足预设判据，将关键支路k放入待考察支路集合，将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询，判定是否能够匹配到割集；
控制模块204，在匹配到多个可解列割集后，从多个割集中筛选出唯一的待解列割集，对待解列割集内的所有支路同时实施解列控制。
其中，量测电力系统中支路参数，根据支路参数确定电力系统中每一条支路的sBTTC指数，根据sBTTC指数定位关键支路，具体包括：
量测电力系统每一条支路两端节点的电压幅值Umi和Uni和电压相位Δθi，i为电力系统内支路总数，并根据下式计算每一条支路的sBTTC指数；
根据如下公式：
其中，Δθi＝|Δθmi-Δθni|，m、n为示一条线路上取的两个量测点的位置，i为第i条支路；
对每条支路的sBTTC指数从小到大进行排序，定义具有最小sBTTC指数值的支路为关键支路，并将关键支路的支路编号记为k。
其中，根据sBTTC指数定位关键支路后，定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
其中，控制参数还包括：电力系统中待考察支路总个数Nb、用于判定系统失稳的严重程度的门槛值Uvth。
其中，关键支路k满足预设判据，具体包括：
当电力系统发生功角失稳时，电力系统内支路两端节点的电压相位差会不断增大，关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小，即有功功率满足如下：
且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、n的电压相位差Δθ增大而不断增大，即频率偏差Δf满足如下：
定义支路垂足电压Uv的计算公式如下：
其中，umy＝|Im(Um)|、unx＝|Re(Un)|、umx＝|Re(Um)|，即，umx、umy分别是Um实部、虚部的绝对值，unx是Un实部的绝对值，Um为关键支路k端点m的电压，Un为关键支路k端点n的电压；
为判断支路垂足电压Uv是否位于支路之上，则ξv满足：
其中，ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、Um节点m与Uv的相位差、Un节点n与Uv的相位差，ΔUmn、ΔUmv、ΔUnv根据下式计算：
关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据，如下：
Uvk<Uvth。
其中，若将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询，不能匹配到可解列割集，判定是否满足kq<Nb；
其中，kq为提取的具有次小sBTTC指数的支路累计个数，其初值设置为1，Nb为待考察支路总个数，最小值为1，最大值为系统内全部支路的总个数；
若满足，令kq＝kq+1再次进行匹配查询；
若不满足，判定全部待考察支路中是否能够有构成割集的支路组合，若有则再次进行匹配查询，若无，则将故障清除后的延时时间加上Δt时间后，再次判断故障清除后的延时时间是否大于TFS，直到故障清除后的延时时间是大于TFS。
其中，从多个割集中筛选出唯一的待解列割集，计算公式如下：
根据下式：
其中，asTSC功角失稳目标割集，Pscl0为支路l的稳态有功功率，Nsc为可选解列割集个数。
其中，若关键支路k不满足预设判据，关键支路k的sBTTC指数值大于设定门槛值εLth，或故障清除后的延时时间小于TFS，将故障清除后的延时时间加上Δt时间后，再次判断故障清除后的延时时间是否大于TFS，直到故障清除后的延时时间是大于TFS。
本发明在系统发生功角失稳时，可以通过对搜索得到的解列割集支路实施主动解列控制，快速将失步电网解列，从而避免机组功角摇摆引起网络节点电压降低所造成的连锁故障，缩小失稳影响范围。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本发明实施例中的方案可以采用各种计算机语言实现，例如，面向对象的程序设计语言Java和直译式脚本语言JavaScript等。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。
Claims (12)
1.一种用于电力系统功角失稳的主动解列控制方法，所述方法包括：
针对电力系统的功角失稳情况，设置电力系统的控制参数，所述控制参数包括：为避免在电力系统功角失稳进行故障清除后，初始阶段电气量波动影响设置的延时时间TFS、用于判定系统的功角稳定性是否已经恶化的门槛值εasth；
在电力系统发生功角失稳后，判断故障清除后的延时时间是否大于TFS，若大于，量测电力系统中支路参数，根据支路参数确定电力系统中每一条支路的sBTTC指数，根据sBTTC指数定位关键支路；
若关键支路k的sBTTC指数值小于设定门槛值εLth，且关键支路k满足预设判据，将关键支路k放入待考察支路集合，将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询，判定是否能够匹配到割集；
若匹配到多个可解列割集，从多个割集中筛选出唯一的待解列割集，对待解列割集内的所有支路同时实施解列控制；
所述量测电力系统中支路参数，根据支路参数确定电力系统中每一条支路的sBTTC指数，根据sBTTC指数定位关键支路，具体包括：
量测电力系统每一条支路两端节点的电压幅值Umi和Uni和电压相位Δθi，i为电力系统内支路总数，并根据下式计算每一条支路的sBTTC指数；
根据如下公式：
其中，Δθi＝|Δθmi-Δθni|，m、n为示一条线路上取的两个量测点的位置，i为第i条支路；
对每条支路的sBTTC指数从小到大进行排序，定义具有最小sBTTC指数值的支路为关键支路，并将关键支路的支路编号记为k；
所述关键支路k满足预设判据，具体包括：
当电力系统发生功角失稳时，电力系统内支路两端节点的电压相位差会不断增大，关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小，即有功功率满足如下：
且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、n的电压相位差Δθ增大而不断增大，即频率偏差Δf满足如下：
定义支路垂足电压Uv的计算公式如下：
其中，umy＝|Im(Um)|、unx＝|Re(Un)|、umx＝|Re(Um)|，即，umx、umy分别是Um实部、虚部的绝对值，unx是Un实部的绝对值，Um为关键支路k端点m的电压，Un为关键支路k端点n的电压；
为判断支路垂足电压Uv是否位于支路之上，则ξv满足：
其中，ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、Um节点m与Uv的相位差、Un节点n与Uv的相位差，ΔUmn、ΔUmv、ΔUnv根据下式计算：
关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据，如下：
Uvk<Uvth。
2.根据权利要求1所述的方法，所述根据sBTTC指数定位关键支路后，定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
3.根据权利要求1所述的方法，所述控制参数还包括：电力系统中待考察支路总个数Nb、用于判定系统失稳的严重程度的门槛值Uvth。
4.根据权利要求1所述的方法，所述将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询，不能匹配到可解列割集，判定是否满足kq<Nb；
其中，kq为提取的具有次小sBTTC指数的支路累计个数，其初值设置为1，Nb为待考察支路总个数，最小值为1，最大值为系统内全部支路的总个数；
若满足，令kq＝kq+1再次进行匹配查询；
若不满足，判定全部待考察支路中是否能够有构成割集的支路组合，若有则再次进行匹配查询，若无，则将故障清除后的延时时间加上Δt时间后，再次判断故障清除后的延时时间是否大于TFS，直到故障清除后的延时时间是大于TFS。
5.根据权利要求1所述的方法，所述从多个割集中筛选出唯一的待解列割集，计算公式如下：
根据下式：
其中，asTSC功角失稳目标割集，Pscl0为支路l的稳态有功功率，Nsc为可选解列割集个数。
6.根据权利要求1所述的方法，所述若关键支路k不满足预设判据，关键支路k的sBTTC指数值大于设定门槛值εLth，或故障清除后的延时时间小于TFS，将故障清除后的延时时间加上Δt时间后，再次判断故障清除后的延时时间是否大于TFS，直到故障清除后的延时时间是大于TFS。
7.一种用于电力系统功角失稳的主动解列控制系统，所述系统包括：
参数确定模块，针对电力系统的功角失稳情况，设置电力系统的控制参数，所述控制参数包括：为避免在电力系统功角失稳进行故障清除后，初始阶段电气量波动影响设置的延时时间TFS、用于判定系统的功角稳定性是否已经恶化的门槛值εasth；
定位模块，在电力系统发生功角失稳后，判断故障清除后的延时时间是否大于TFS，若大于，量测电力系统中支路参数，根据支路参数确定电力系统中每一条支路的sBTTC指数，根据sBTTC指数定位关键支路；
判定模块，在关键支路k的sBTTC指数值小于设定门槛值εLth时，且关键支路k满足预设判据，将关键支路k放入待考察支路集合，将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询，判定是否能够匹配到割集；
控制模块，在匹配到多个可解列割集后，从多个割集中筛选出唯一的待解列割集，对待解列割集内的所有支路同时实施解列控制；
所述量测电力系统中支路参数，根据支路参数确定电力系统中每一条支路的sBTTC指数，根据sBTTC指数定位关键支路，具体包括：
量测电力系统每一条支路两端节点的电压幅值Umi和Uni和电压相位Δθi，i为电力系统内支路总数，并根据下式计算每一条支路的sBTTC指数；
根据如下公式：
其中，Δθi＝|Δθmi-Δθni|，m、n为示一条线路上取的两个量测点的位置，i为第i条支路；
对每条支路的sBTTC指数从小到大进行排序，定义具有最小sBTTC指数值的支路为关键支路，并将关键支路的支路编号记为k；
所述关键支路k满足预设判据，具体包括：
当电力系统发生功角失稳时，电力系统内支路两端节点的电压相位差会不断增大，关键支路k上有功功率P会随着支路k两端节点m、n的电压相位差Δθ增大而不断减小，即有功功率满足如下：
且关键支路k两端节点m、n的频率偏差Δf会随着支路k两端节点m、n的电压相位差Δθ增大而不断增大，即频率偏差Δf满足如下：
定义支路垂足电压Uv的计算公式如下：
其中，umy＝|Im(Um)|、unx＝|Re(Un)|、umx＝|Re(Um)|，即，umx、umy分别是Um实部、虚部的绝对值，unx是Un实部的绝对值，Um为关键支路k端点m的电压，Un为关键支路k端点n的电压；
为判断支路垂足电压Uv是否位于支路之上，则ξv满足：
其中，ΔUmn、ΔUmv、ΔUnv分别为支路两端节点m、n的电压相位差、Um节点m与Uv的相位差、Un节点n与Uv的相位差，ΔUmn、ΔUmv、ΔUnv根据下式计算：
关键支路k的垂足电压Uvk低于设定门槛值Uvth的判据，如下：
Uvk<Uvth。
8.根据权利要求7所述的系统，所述根据sBTTC指数定位关键支路后，定义具有次小sBTTC指数值的若干条支路为潜在的能与关键支路形成割集的准关键支路簇。
9.根据权利要求7所述的系统，所述控制参数还包括：电力系统中待考察支路总个数Nb、用于判定系统失稳的严重程度的门槛值Uvth。
10.根据权利要求7所述的系统，所述将待考察支路集合中的每一条支路与离线设定的可解列割集进行匹配查询，不能匹配到可解列割集，判定是否满足kq<Nb；
其中，kq为提取的具有次小sBTTC指数的支路累计个数，其初值设置为1，Nb为待考察支路总个数，最小值为1，最大值为系统内全部支路的总个数；
若满足，令kq＝kq+1再次进行匹配查询；
若不满足，判定全部待考察支路中是否能够有构成割集的支路组合，若有则再次进行匹配查询，若无，则将故障清除后的延时时间加上Δt时间后，再次判断故障清除后的延时时间是否大于TFS，直到故障清除后的延时时间是大于TFS。
11.根据权利要求7所述的系统，所述从多个割集中筛选出唯一的待解列割集，计算公式如下：
根据下式：
其中，asTSC功角失稳目标割集，Pscl0为支路l的稳态有功功率，Nsc为可选解列割集个数。
12.根据权利要求7所述的系统，所述关键支路k不满足预设判据，关键支路k的sBTTC指数值大于设定门槛值εLth，或故障清除后的延时时间小于TFS，将故障清除后的延时时间加上Δt时间后，再次判断故障清除后的延时时间是否大于TFS，直到故障清除后的延时时间是大于TFS。
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 CN113765075A (zh) | 2021-12-07 |
| CN113765075B true 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)
| 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 | 中国电力科学研究院有限公司 | 一种基于低频低压的主动解列控制方法及系统 |
-
2021
- 2021-07-30 CN CN202110877923.5A patent/CN113765075B/zh active Active
Patent Citations (4)
| 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 |
|---|---|
| CN113765075A (zh) | 2021-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105893740B (zh) | 一种用于挖掘电力网络中脆弱节点和线路的方法 | |
| EP2720332B1 (en) | Online protection coordination for distribution management system | |
| Mahmood et al. | Real-time reduced steady-state model synthesis of active distribution networks using PMU measurements | |
| CN114254818B (zh) | 一种基于风险评估模型的低压配电网连锁故障预警方法 | |
| CN106548265B (zh) | 一种基于连锁故障事故链搜索的输电网可靠性评估方法 | |
| CN109936113B (zh) | 一种基于随机森林算法的保护动作智能诊断方法及系统 | |
| CN110166483B (zh) | 一种电网故障以及网络攻击的辨识方法、装置及设备 | |
| CN109375050B (zh) | 一种基于初始故障线路割集的电力系统关键断面辨识方法 | |
| CN110147827B (zh) | 一种基于iaalo-svm与相似性度量的故障预测方法 | |
| CN108921725A (zh) | 一种复杂配电网可靠性指标快速解析计算方法 | |
| CN116152016A (zh) | 一种配电网故障定位方法、装置、设备及存储介质 | |
| CN113765075B (zh) | 一种用于电力系统功角失稳的主动解列控制方法及系统 | |
| CN104852383A (zh) | 一种海上油田群电网的紧急控制方法 | |
| CN109873418B (zh) | 一种暂态功角稳定模式的强相关输电断面识别方法及系统 | |
| CN104462404A (zh) | 电力系统潮流转移危险线路快速搜索方法 | |
| CN113824094B (zh) | 一种用于电力系统电压失稳的主动解列控制方法及系统 | |
| CN113922343B (zh) | 配电网保护定值整定方法、终端设备和存储介质 | |
| CN109167353B (zh) | 暂态功角稳定视角下的电网强关联断面识别方法 | |
| CN115566678A (zh) | 一种基于综合评价指标的电力系统网络重构方法 | |
| CN115879271A (zh) | 一种基于割集断面树的电网割集型关键输电断面搜索方法 | |
| CN113918416A (zh) | 一种基于事件图谱的软件系统异常定位方法 | |
| CN114511194A (zh) | 一种电力物联网的运行风险预测方法、系统和电子设备 | |
| CN110994564A (zh) | 配电网保护投入开关自动选取和定值整定方法 | |
| CN112713584B (zh) | 适用于稳控系统预想故障集外故障的控制方法和装置 | |
| CN117741345B (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 |