CN113341210B - 母线暂态电压稳定性评估方法、系统、介质及计算设备 - Google Patents
母线暂态电压稳定性评估方法、系统、介质及计算设备 Download PDFInfo
- Publication number
- CN113341210B CN113341210B CN202110549792.8A CN202110549792A CN113341210B CN 113341210 B CN113341210 B CN 113341210B CN 202110549792 A CN202110549792 A CN 202110549792A CN 113341210 B CN113341210 B CN 113341210B
- Authority
- CN
- China
- Prior art keywords
- distribution network
- bus
- moment
- network bus
- 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
- 230000001052 transient effect Effects 0.000 title claims abstract description 94
- 238000011156 evaluation Methods 0.000 title claims description 8
- 238000000034 method Methods 0.000 claims abstract description 54
- 238000013097 stability assessment Methods 0.000 claims abstract description 18
- 238000009826 distribution Methods 0.000 claims description 147
- 238000005070 sampling Methods 0.000 claims description 102
- 230000010355 oscillation Effects 0.000 claims description 71
- 230000008569 process Effects 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 8
- 230000015654 memory Effects 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000368 destabilizing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/165—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values
- G01R19/16533—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values characterised by the application
- G01R19/16538—Indicating that current or voltage is either above or below a predetermined value or within or outside a predetermined range of values characterised by the application in AC or DC supplies
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种母线暂态电压稳定性评估方法、系统、介质及计算设备，本发明在两区域发生功角振荡期间，将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，计算各采样时刻配网母线的暂态电压稳定裕度，根据功角振荡期间配网母线暂态电压稳定裕度正负变化情况，判断出母线是否暂态电压失稳，通过实时跟踪电网运行状态，及时识别出由于功角振荡而发生暂态电压失稳的母线，评估结果能够为紧急控制提供辅助决策，从而能够对发生暂态电压失稳的母线采取有针对性的切负荷措施，既能防止事故进一步扩大，又能最大限度地减少负荷损失，有助于提高电网供电的可靠性和运行的稳定性。
Description
技术领域
本发明涉及一种母线暂态电压稳定性评估方法、系统、介质及计算设备，属于电力系统及其自动化技术领域。
背景技术
故障或扰动会导致互联的两区域电网之间发生功角振荡。在功角振荡过程中，两区域联络线各点电压也会随之波动。联络线上同一点的电压幅值随着功角摆开幅度的增大而降低。在实际运行中，当系统发生故障或扰动后，经常发生两互联区域功角振荡但并没有发生失步，然而两区域之间联络线所接负荷母线却发生暂态电压失稳的问题。由于功角振荡几个周期之后很快恢复稳定，所以负荷母线发生暂态电压失稳具有突发性和隐蔽性。虽然刚开始只是局部负荷母线发生暂态电压失稳，但是如果不及时针对失稳母线采取相应的紧急控制措施，往往会演化为多种失稳形态交织的大范围失稳问题，大大增加了分析和控制的难度，从而造成更大的负荷损失。
随着大区互联电网规模的日益扩大，我国电网结构日趋复杂，更容易发生故障或扰动引起的区域间功角振荡。当互联的两个区域相对较远时，需要通过中间的枢纽变电站向两区域之间的负荷供电。因此，经枢纽变电站互联的两个区域电网发生功角相对振荡时，通过枢纽变电站供电的配网负荷母线容易发生暂态电压失稳。目前，我国判断负荷母线暂态电压失稳的方法依然沿用工程经验判据，即通过电压跌落门槛值和持续时间的二元表来判断暂态电压稳定性。我国《电力系统安全稳定导则》(GB 38755-2019)将暂态电压稳定的判据阐述为：在电力系统受到扰动后的暂态或者动态过程中，负荷母线电压能够恢复到规定的运行电压水平以上。该判据是在没有准确的暂态电压稳定判别方法时采用的一种权宜之计，很大程度上取决于以往的计算及运行经验。另外，运用现有方法评估暂态电压失稳往往需要考察从故障开始后的10s以内的电压响应情况，对于实际电网来说，这种方法只适用于离线计算分析，不适用于紧急控制。
因此传统的母线暂态电压稳定性评估已不适应电网的发展，现在急需一种新的母线暂态电压稳定性评估方法。
发明内容
本发明提供了一种母线暂态电压稳定性评估方法、系统、介质及计算设备，解决了背景技术中披露的问题。
为了解决上述技术问题，本发明所采用的技术方案是：
母线暂态电压稳定性评估方法，包括：
在两区域发生功角振荡期间，将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，计算各采样时刻配网母线的暂态电压稳定裕度；其中，配网母线连接枢纽变电站母线；
若T1到T2期间，配网母线的所有暂态电压稳定裕度均小于0，则判断配网母线在该功角振荡过程中发生暂态电压失稳；其中，T1为配网母线的暂态电压稳定裕度第一次小于0的时刻，T2为功角振荡结束时刻。
确定功角振荡的起始时刻和结束时刻的过程为，
对两区域母线的电压相角进行实时采样，计算采样时刻两区域母线电压相角差的变化量；
若前一采样时刻变化量的绝对值等于0、且当前采样时刻变化量的绝对值大于0，则当前采样时刻为功角振荡的起始时刻；
若前一采样时刻变化量的绝对值大于0、且当前采样时刻变化量的绝对值等于0，则当前采样时刻为功角振荡的结束时刻。
计算采样时刻两区域母线电压相角差的变化量，具体公式为，
其中，ai、bi分别为两区域母线；t为采样时刻；t-1采样时刻的前一刻；为t时刻两区域母线电压相角差的变化量；/>为t时刻两区域母线电压相角差；/>为t-1时刻两区域母线电压相角差。
将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，计算采样时刻配网母线的暂态电压稳定裕度，具体过程为，
将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，计算采样时刻配网母线的负荷等值阻抗、以及采样时刻从配网母线看进去的戴维南等值阻抗；
根据负荷等值阻抗、戴维南等值阻抗和配网母线电压相量，计算采样时刻从配网母线看进去的戴维南等值电势；
根据戴维南等值阻抗，计算采样时刻配网母线的戴维南等值阻抗阻抗角；
根据采样时刻配网母线的负荷有功功率和负荷无功功率，计算采样时刻配网母线的负荷功率因数角；
根据负荷母线电压幅值、戴维南等值电势幅值、戴维南等值阻抗阻抗角和负荷功率因数角，计算采样时刻配网母线的暂态电压稳定裕度。
将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，计算采样时刻配网母线的负荷等值阻抗、以及采样时刻从配网母线看进去的戴维南等值阻抗，具体过程为，
将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，简化电力系统；
根据采样时刻配网母线的负荷有功功率、负荷无功功率和电压幅值，计算采样时刻配网母线的负荷等值阻抗；
根据简化电力系统的网络节点导纳矩阵，计算采样时刻从配网母线看进去的戴维南等值阻抗。
计算采样时刻配网母线的负荷等值阻抗，具体公式为，
其中，t为采样时刻；为t时刻配网母线j的负荷等值阻抗；/>分别为t时刻配网母线j的负荷有功功率和负荷无功功率；Uj,t为t时刻配网母线j的电压幅值。
计算采样时刻从配网母线看进去的戴维南等值电势，具体公式为，
其中，t为采样时刻；为t时刻从配网母线j看进去的戴维南等值电势；/>为t时刻配网母线j的负荷等值阻抗；Zjj,t为t时刻从配网母线j看进去的戴维南等值阻抗；为t时刻配网母线j的电压相量。
计算采样时刻配网母线的戴维南等值阻抗阻抗角，具体公式为，
其中，t为采样时刻；ξj,t为t时刻配网母线j的戴维南等值阻抗阻抗角；分别为戴维南等值阻抗的虚部和实部。
计算采样时刻配网母线的负荷功率因数角，具体公式为，
其中，t为采样时刻；为t时刻配网母线j的负荷功率因数角；/>分别为t时刻配网母线j的负荷有功功率和负荷无功功率。
计算采样时刻配网母线的暂态电压稳定裕度，具体公式为，
其中，t为采样时刻；ηj,t为t时刻配网母线j的暂态电压稳定裕度；Uj,t为t时刻配网母线j的电压幅值；
为t时刻配网母线j的临界电压，/>为t时刻从配网母线j看进去的戴维南等值电势幅值，ξj,t为t时刻配网母线j的戴维南等值阻抗阻抗角；/>为t时刻配网母线j的负荷功率因数角。
母线暂态电压稳定性评估系统，包括：
暂态电压稳定裕度计算模块：在两区域发生功角振荡期间，将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，计算各采样时刻配网母线的暂态电压稳定裕度；其中，配网母线连接枢纽变电站母线；
判断模块：若T1到T2期间，配网母线的所有暂态电压稳定裕度均小于0，则判断配网母线在该功角振荡过程中发生暂态电压失稳；其中，T1为配网母线的暂态电压稳定裕度第一次小于0的时刻，T2为功角振荡结束时刻。
一种存储一个或多个程序的计算机可读存储介质，所述一个或多个程序包括指令，所述指令当由计算设备执行时，使得所述计算设备执行母线暂态电压稳定性评估方法。
一种计算设备，包括一个或多个处理器、一个或多个存储器以及一个或多个程序，其中一个或多个程序存储在所述一个或多个存储器中并被配置为由所述一个或多个处理器执行，所述一个或多个程序包括用于执行母线暂态电压稳定性评估方法的指令。
本发明所达到的有益效果：本发明在两区域发生功角振荡期间，将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，计算各采样时刻配网母线的暂态电压稳定裕度，根据功角振荡期间配网母线暂态电压稳定裕度正负变化情况，判断出母线是否暂态电压失稳，通过实时跟踪电网运行状态，及时识别出由于功角振荡而发生暂态电压失稳的母线，评估结果能够为紧急控制提供辅助决策，从而能够对发生暂态电压失稳的母线采取有针对性的切负荷措施，既能防止事故进一步扩大，又能最大限度地减少负荷损失，有助于提高电网供电的可靠性和运行的稳定性。
附图说明
图1为拓扑关系图；
图2为简化后的拓扑关系图；
图3为本发明方法的流程图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
母线暂态电压稳定性评估方法，包括以下步骤：
步骤1，在两区域发生功角振荡期间，将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，计算各采样时刻配网母线的暂态电压稳定裕度；其中，配网母线连接枢纽变电站母线；
步骤2，若T1到T2期间，配网母线j的所有暂态电压稳定裕度均小于0，则判断配网母线j在该功角振荡过程中发生暂态电压失稳；其中，T1为配网母线的暂态电压稳定裕度第一次小于0的时刻，T2为功角振荡结束时刻。
上述方法主要运用在存在若干如图1所示拓扑关系的电力系统，图中区域A和B通过枢纽变电站连接配网，具体为：配网母线j连接枢纽变电站母线i，枢纽变电站母线i连接两区域母线ai、bi，配网母线j是枢纽变电站母线i所接配网区域中的一条母线。当两区域发生功角振荡时，电力系统采用上述方法对所有拓扑关系(图1)中的配网母线和枢纽变电站母线进行评估，输出所有暂态电压失稳母线，直到检测到功角振荡停止。
这里通过两区域母线电压相角差的变化量判断是否发生功角振荡，即通过两区域母线电压相角差的变化量确定功角振荡的起始时刻和结束时刻，具体过程如下：
A)对两区域母线的电压相角进行实时采样，计算采样时刻两区域母线电压相角差的变化量；
对两区域母线的电压相角进行实时采样，分别记为同一采样时刻，两个电压相角的相角差为/>通过相邻时刻的相角差可计算出采样时刻两区域母线电压相角差的变化量，具体公式如下：
其中，ai、bi分别为两区域母线；t为采样时刻；t-1采样时刻的前一刻；为t时刻两区域母线电压相角差的变化量；/>为t时刻两区域母线电压相角差；/>为t-1时刻两区域母线电压相角差；
B)若前一采样时刻变化量的绝对值等于0、且当前采样时刻变化量的绝对值大于0，则当前采样时刻为功角振荡的起始时刻；若前一采样时刻变化量的绝对值大于0、且当前采样时刻变化量的绝对值等于0，则当前采样时刻为功角振荡的结束时刻；
如果则判定此刻t产生功角振荡，即t为功角振荡的起始时刻；如果/>则判定此刻t功角振荡结束，即t为功角振荡的结束时刻。
在两区域发生功角振荡期间，将与两区域母线连接的枢纽变电站母线i视为电动势幅值随功角振荡变化的电压源，计算采样时刻配网母线j的暂态电压稳定裕度，具体过程为：
11)将与两区域母线连接的枢纽变电站母线i视为电动势幅值随功角振荡变化的电压源，简化电力系统，具体简化图1的拓扑关系，简化后的拓扑关系如图2所示。
12)根据采样时刻配网母线j的负荷有功功率、负荷无功功率和电压幅值，计算采样时刻配网母线j的负荷等值阻抗；
其中，t为采样时刻；为t时刻配网母线j的负荷等值阻抗；Uj,t为t时刻配网母线j的电压幅值，该值可从变电站的就地测量装置获取；/>分别为t时刻配网母线j的负荷有功功率和负荷无功功率。
13)根据简化电力系统的网络节点导纳矩阵，计算采样时刻从配网母线j看进去的戴维南等值阻抗；
网络节点导纳矩阵Yt包括各节点的自导纳和互导纳，具体可表示为：
其中，Ynn为节点n的自导纳；yn0为节点n的对地导纳；ynk为节点n和k之间的支路导纳；Ynk为节点n和节点k的互导纳；
利用三角分解法求解YtUt＝It，其中，Ut为t时刻配网各独立节点电压列向量，It为t时刻配网各独立节点注入电流列向量，在求解过程中，需要求哪个节点的自阻抗，就将It中对应该节点的注入电流设为1，其余节点注入电流设为0即可，Yt的矩阵形式如下：
将YtUt＝It左右两边同时左乘Yt -1即可得Ut＝Yt -1It，由于Yt -1＝Zt，因此，Ut＝ZtIt，将该式写为详细的矩阵表达式为：
计算上式可得，从而可得，/>
利用三角分解法求得配网母线j处的电压即为配网母线j的自阻抗Zjj,t，也就是从配网母线j处看进去的戴维南等值阻抗/>其中，/>
15)根据负荷等值阻抗、戴维南等值阻抗和配网母线电压相量，计算采样时刻从配网母线j看进去的戴维南等值电势；具体采用的公式如下：
其中，为t时刻从配网母线j看进去的戴维南等值电势；Zjj,t为t时刻从配网母线j看进去的戴维南等值阻抗；/>为t时刻配网母线j的电压相量，可通过本地量测装置得到。
16)根据戴维南等值阻抗，计算采样时刻配网母线j的戴维南等值阻抗阻抗角；具体采用的公式如下：
其中，ξj,t为t时刻配网母线j的戴维南等值阻抗阻抗角；分别为戴维南等值阻抗的虚部和实部。
17)根据采样时刻配网母线j的负荷有功功率和负荷无功功率，计算采样时刻配网母线j的负荷功率因数角；具体采用的公式如下：
其中，t为采样时刻；为t时刻配网母线j的负荷功率因数角。
18)根据配网母线电压幅值、戴维南等值电势幅值、戴维南等值阻抗阻抗角和负荷功率因数角，计算采样时刻配网母线j的暂态电压稳定裕度；具体采用的公式如下：
其中，t为采样时刻；ηj,t为t时刻配网母线j的暂态电压稳定裕度；Uj,t为t时刻配网母线j的电压幅值；为t时刻配网母线j的临界电压；
根据戴维南等值电势幅值、戴维南等值阻抗阻抗角和负荷功率因数角计算获得，具体公式如下：
其中，为t时刻从配网母线j看进去的戴维南等值电势幅值；ξj,t为t时刻配网母线j的戴维南等值阻抗阻抗角；/>为t时刻配网母线j的负荷功率因数角。
在整个功角振荡期间内，获得各采样时刻配网母线的暂态电压稳定裕度，查看暂态电压稳定裕度的变化，定义暂态电压稳定裕度第一次小于0的时刻为T1，功角振荡结束时刻为T2，若T1到T2期间，配网母线的所有暂态电压稳定裕度均小于0，即则判断配网母线在该功角振荡过程中发生暂态电压失稳，将其放入暂态电压失稳负荷母线集，否则配网母线j不为暂态电压失稳母线。
以配网母线j为例，上述方法的具体过程如图3所示，包括：
1)确定目标对象并赋予初值：确定电力系统中所有如图1所示的拓扑关系，设配网母线失稳母线集设功角振荡的起始时刻和结束时刻分别为T0＝0和T2＝0，设功角振荡过程中，配网母线暂态电压稳定裕度第一次由正变负的时刻T1＝0；
2)实时计算ai、bi电压相角差的变化量；
3)若当前时刻(即t时刻)将与两区域母线连接的枢纽变电站母线i视为电动势幅值随功角振荡变化的电压源，简化电力系统，见图2，判断T0是否大于0，若是，则直接转至4)；若不是，则将t赋值给T0并转至步骤4)；
若先判断上一时刻是否发生功角振荡，若没有发生功角振荡则直接结束本方法；若发生功角振荡，则进一步判断T1到T2期间，配网母线j的所有暂态电压稳定裕度均小于0，若是，则判断配网母线j在该功角振荡过程中发生暂态电压失稳，将其加入失稳母线集Hi并输出失稳母线集合，结束本方法；若不是，则判断配网母线j未发生暂态电压失稳，不将其加入失稳母线集，直接结束本方法；
4)计算
5)计算简化电力系统的网络节点导纳矩阵Yt；
6)计算t时刻从配网母线j看进去的戴维南等值阻抗和戴维南等值电势/>
7)计算ξj,t和
8)计算t时刻配网母线j的暂态电压稳定裕度ηj,t；
9)若ηj,t＜0，则判断表征配网母线j暂态电压稳定裕度第一次变为负值的时间T1是否大于0，若是，则直接进入下一时刻的计算；若不是，将当前时间赋值给T1并进入下一时刻的计算；若ηj,t≥0，则直接进入下一时刻的计算。
上述方法在两区域发生功角振荡期间，将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，计算各采样时刻配网母线的暂态电压稳定裕度，根据功角振荡期间配网母线暂态电压稳定裕度正负变化情况，判断出母线是否暂态电压失稳，通过实时跟踪电网运行状态，及时识别出由于功角振荡而发生暂态电压失稳的母线，评估结果能够为紧急控制提供辅助决策，从而能够对发生暂态电压失稳的母线采取有针对性的切负荷措施，既能防止事故进一步扩大，又能最大限度地减少负荷损失，有助于提高电网供电的可靠性和运行的稳定性。
母线暂态电压稳定性评估系统，包括：
暂态电压稳定裕度计算模块：在两区域发生功角振荡期间，将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，计算各采样时刻配网母线的暂态电压稳定裕度；其中，配网母线连接枢纽变电站母线；
判断模块：若T1到T2期间，配网母线的所有暂态电压稳定裕度均小于0，则判断配网母线在该功角振荡过程中发生暂态电压失稳；其中，T1为配网母线的暂态电压稳定裕度第一次小于0的时刻，T2为功角振荡结束时刻。
一种存储一个或多个程序的计算机可读存储介质，所述一个或多个程序包括指令，所述指令当由计算设备执行时，使得所述计算设备执行母线暂态电压稳定性评估方法。
一种计算设备，包括一个或多个处理器、一个或多个存储器以及一个或多个程序，其中一个或多个程序存储在所述一个或多个存储器中并被配置为由所述一个或多个处理器执行，所述一个或多个程序包括用于执行母线暂态电压稳定性评估方法的指令。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上仅为本发明的实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均包含在申请待批的本发明的权利要求范围之内。
Claims (11)
1.母线暂态电压稳定性评估方法，其特征在于，包括：
在两区域发生功角振荡期间，将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，简化电力系统；根据采样时刻配网母线的负荷有功功率、负荷无功功率和电压幅值，计算采样时刻配网母线的负荷等值阻抗；根据简化电力系统的网络节点导纳矩阵，计算采样时刻从配网母线看进去的戴维南等值阻抗；根据负荷等值阻抗、戴维南等值阻抗和配网母线电压相量，计算采样时刻从配网母线看进去的戴维南等值电势；根据戴维南等值阻抗，计算采样时刻配网母线的戴维南等值阻抗阻抗角；根据采样时刻配网母线的负荷有功功率和负荷无功功率，计算采样时刻配网母线的负荷功率因数角；根据配网母线电压幅值、戴维南等值电势幅值、戴维南等值阻抗阻抗角和负荷功率因数角，计算采样时刻配网母线的暂态电压稳定裕度；其中，配网母线连接枢纽变电站母线；
若T1到T2期间，配网母线的所有暂态电压稳定裕度均小于0，则判断配网母线在该功角振荡过程中发生暂态电压失稳；其中，T1为配网母线的暂态电压稳定裕度第一次小于0的时刻，T2为功角振荡结束时刻。
2.根据权利要求1所述的母线暂态电压稳定性评估方法，其特征在于：确定功角振荡的起始时刻和结束时刻的过程为，
对两区域母线的电压相角进行实时采样，计算采样时刻两区域母线电压相角差的变化量；
若前一采样时刻变化量的绝对值等于0、且当前采样时刻变化量的绝对值大于0，则当前采样时刻为功角振荡的起始时刻；
若前一采样时刻变化量的绝对值大于0、且当前采样时刻变化量的绝对值等于0，则当前采样时刻为功角振荡的结束时刻。
3.根据权利要求2所述的母线暂态电压稳定性评估方法，其特征在于：计算采样时刻两区域母线电压相角差的变化量，具体公式为，
其中，ai、bi分别为两区域母线；t为采样时刻；t-1采样时刻的前一刻；为t时刻两区域母线电压相角差的变化量；/>为t时刻两区域母线电压相角差；/>为t-1时刻两区域母线电压相角差。
4.根据权利要求1所述的母线暂态电压稳定性评估方法，其特征在于：计算采样时刻配网母线的负荷等值阻抗，具体公式为，
其中，t为采样时刻；为t时刻配网母线j的负荷等值阻抗；/>分别为t时刻配网母线j的负荷有功功率和负荷无功功率；Uj,t为t时刻配网母线j的电压幅值。
5.根据权利要求1所述的母线暂态电压稳定性评估方法，其特征在于：计算采样时刻从配网母线看进去的戴维南等值电势，具体公式为，
其中，t为采样时刻；为t时刻从配网母线j看进去的戴维南等值电势；/>为t时刻配网母线j的负荷等值阻抗；Zjj,t为t时刻从配网母线j看进去的戴维南等值阻抗；/>为t时刻配网母线j的电压相量。
6.根据权利要求1所述的母线暂态电压稳定性评估方法，其特征在于：计算采样时刻配网母线的戴维南等值阻抗阻抗角，具体公式为，
其中，t为采样时刻；ξj,t为t时刻配网母线j的戴维南等值阻抗阻抗角；分别为戴维南等值阻抗的虚部和实部。
7.根据权利要求1所述的母线暂态电压稳定性评估方法，其特征在于：计算采样时刻配网母线的负荷功率因数角，具体公式为，
其中，t为采样时刻；为t时刻配网母线j的负荷功率因数角；/>分别为t时刻配网母线j的负荷有功功率和负荷无功功率。
8.根据权利要求1所述的母线暂态电压稳定性评估方法，其特征在于：计算采样时刻配网母线的暂态电压稳定裕度，具体公式为，
其中，t为采样时刻；ηj,t为t时刻配网母线j的暂态电压稳定裕度；Uj,t为t时刻配网母线j的电压幅值；
为t时刻配网母线j的临界电压，/>为t时刻从配网母线j看进去的戴维南等值电势幅值，ξj,t为t时刻配网母线j的戴维南等值阻抗阻抗角；为t时刻配网母线j的负荷功率因数角。
9.母线暂态电压稳定性评估系统，其特征在于，包括：
暂态电压稳定裕度计算模块：在两区域发生功角振荡期间，将与两区域母线连接的枢纽变电站母线视为电动势幅值随功角振荡变化的电压源，简化电力系统；根据采样时刻配网母线的负荷有功功率、负荷无功功率和电压幅值，计算采样时刻配网母线的负荷等值阻抗；根据简化电力系统的网络节点导纳矩阵，计算采样时刻从配网母线看进去的戴维南等值阻抗；根据负荷等值阻抗、戴维南等值阻抗和配网母线电压相量，计算采样时刻从配网母线看进去的戴维南等值电势；根据戴维南等值阻抗，计算采样时刻配网母线的戴维南等值阻抗阻抗角；根据采样时刻配网母线的负荷有功功率和负荷无功功率，计算采样时刻配网母线的负荷功率因数角；根据配网母线电压幅值、戴维南等值电势幅值、戴维南等值阻抗阻抗角和负荷功率因数角，计算采样时刻配网母线的暂态电压稳定裕度；其中，配网母线连接枢纽变电站母线；
判断模块：若T1到T2期间，配网母线的所有暂态电压稳定裕度均小于0，则判断配网母线在该功角振荡过程中发生暂态电压失稳；其中，T1为配网母线的暂态电压稳定裕度第一次小于0的时刻，T2为功角振荡结束时刻。
10.一种存储一个或多个程序的计算机可读存储介质，其特征在于：所述一个或多个程序包括指令，所述指令当由计算设备执行时，使得所述计算设备执行根据权利要求1至8所述的方法中的任一方法。
11.一种计算设备，其特征在于，包括：
一个或多个处理器、一个或多个存储器以及一个或多个程序，其中一个或多个程序存储在所述一个或多个存储器中并被配置为由所述一个或多个处理器执行，所述一个或多个程序包括用于执行根据权利要求1至8所述的方法中的任一方法的指令。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110549792.8A CN113341210B (zh) | 2021-05-20 | 2021-05-20 | 母线暂态电压稳定性评估方法、系统、介质及计算设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110549792.8A CN113341210B (zh) | 2021-05-20 | 2021-05-20 | 母线暂态电压稳定性评估方法、系统、介质及计算设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113341210A CN113341210A (zh) | 2021-09-03 |
| CN113341210B true CN113341210B (zh) | 2024-02-23 |
Family
ID=77469955
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110549792.8A Active CN113341210B (zh) | 2021-05-20 | 2021-05-20 | 母线暂态电压稳定性评估方法、系统、介质及计算设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113341210B (zh) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101572410A (zh) * | 2009-06-01 | 2009-11-04 | 国网电力科学研究院 | 电力系统暂态安全稳定模式中元件参与因子识别方法 |
| CN102983573A (zh) * | 2012-11-09 | 2013-03-20 | 天津大学 | 基于安全域的安全约束经济调度方法 |
| CN103236719A (zh) * | 2013-04-17 | 2013-08-07 | 国家电网公司 | 风电和火电接入电网后的风电打捆控制方法 |
| CN103762579A (zh) * | 2013-12-22 | 2014-04-30 | 南京南瑞集团公司 | 一种提高电网暂态功角稳定性的直流功率紧急控制性能指标计算方法 |
| CN109038538A (zh) * | 2017-06-09 | 2018-12-18 | 刘光晔 | 一种用广义电阻裕度评价静态电压稳定与功角稳定的方法 |
| JP2019030065A (ja) * | 2017-07-26 | 2019-02-21 | 株式会社東芝 | 電力系統信頼度評価システム |
| CN109638815A (zh) * | 2018-12-04 | 2019-04-16 | 国电南瑞科技股份有限公司 | 一种电力系统中长期电压安全稳定预防控制策略确定方法 |
| CN110165660A (zh) * | 2019-04-26 | 2019-08-23 | 国电南瑞科技股份有限公司 | 一种受端电网暂态电压预防控制优化方法及其系统 |
| CN110569485A (zh) * | 2019-08-28 | 2019-12-13 | 云南电网有限责任公司 | 一种计及新能源电源出力波动性的新能源电源节点状态脆弱性监测方法 |
| CN110571790A (zh) * | 2019-03-21 | 2019-12-13 | 中国电力科学研究院有限公司 | 基于戴维南等值的在线电压稳定预警的方法、系统 |
| CN111327058A (zh) * | 2020-03-02 | 2020-06-23 | 国网宁夏电力有限公司 | 基于暂稳态电压约束的直流无功控制死区整定方法和装置 |
| CN111371089A (zh) * | 2020-03-10 | 2020-07-03 | 南瑞集团有限公司 | 一种电网动态等值的量化评估方法及系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9876356B2 (en) * | 2014-10-02 | 2018-01-23 | Mitsubishi Electric Research Laboratories, Inc. | Dynamic and adaptive configurable power distribution system |
-
2021
- 2021-05-20 CN CN202110549792.8A patent/CN113341210B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101572410A (zh) * | 2009-06-01 | 2009-11-04 | 国网电力科学研究院 | 电力系统暂态安全稳定模式中元件参与因子识别方法 |
| CN102983573A (zh) * | 2012-11-09 | 2013-03-20 | 天津大学 | 基于安全域的安全约束经济调度方法 |
| CN103236719A (zh) * | 2013-04-17 | 2013-08-07 | 国家电网公司 | 风电和火电接入电网后的风电打捆控制方法 |
| CN103762579A (zh) * | 2013-12-22 | 2014-04-30 | 南京南瑞集团公司 | 一种提高电网暂态功角稳定性的直流功率紧急控制性能指标计算方法 |
| CN109038538A (zh) * | 2017-06-09 | 2018-12-18 | 刘光晔 | 一种用广义电阻裕度评价静态电压稳定与功角稳定的方法 |
| JP2019030065A (ja) * | 2017-07-26 | 2019-02-21 | 株式会社東芝 | 電力系統信頼度評価システム |
| CN109638815A (zh) * | 2018-12-04 | 2019-04-16 | 国电南瑞科技股份有限公司 | 一种电力系统中长期电压安全稳定预防控制策略确定方法 |
| CN110571790A (zh) * | 2019-03-21 | 2019-12-13 | 中国电力科学研究院有限公司 | 基于戴维南等值的在线电压稳定预警的方法、系统 |
| CN110165660A (zh) * | 2019-04-26 | 2019-08-23 | 国电南瑞科技股份有限公司 | 一种受端电网暂态电压预防控制优化方法及其系统 |
| CN110569485A (zh) * | 2019-08-28 | 2019-12-13 | 云南电网有限责任公司 | 一种计及新能源电源出力波动性的新能源电源节点状态脆弱性监测方法 |
| CN111327058A (zh) * | 2020-03-02 | 2020-06-23 | 国网宁夏电力有限公司 | 基于暂稳态电压约束的直流无功控制死区整定方法和装置 |
| CN111371089A (zh) * | 2020-03-10 | 2020-07-03 | 南瑞集团有限公司 | 一种电网动态等值的量化评估方法及系统 |
Non-Patent Citations (5)
| Title |
|---|
| Application_of_FACTS_Devices_in_a_Multi-Machine_Power_System_Network_for_Transient_Stability_Enhancement_A_Case_Study_of_the_Nigerian_330kV_48-Bus_System;Imo E. Nkan;2019 IEEE AFRICON;1-9 * |
| 在线暂态安全稳定评估的分类滚动故障筛选方法;徐泰山;《电力系统自动化》;第42卷(第13期);182-188 * |
| 基于失稳机理分析的暂态电压稳定性量化指标;李宣;《电力自动化设备》;第43卷(第11期);195-202 * |
| 电力系统安全稳定预防控制在线计算方法的评述;方勇杰;《电力系统自动化》;第39卷(第11期);199-207 * |
| 电力系统戴维南等值分析应用及DG灵活接入后的影响;王宾;《华北电力技术》;29-34 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113341210A (zh) | 2021-09-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zeng et al. | Online estimation of power system inertia constant under normal operating conditions | |
| Yuan et al. | A comparative study of measurement-based Thevenin equivalents identification methods | |
| JP7219542B2 (ja) | 電力系統制御装置、電力系統制御システム、電力系統制御方法及び電力系統制御プログラム | |
| CN103762581A (zh) | 一种同塔双回直流输电系统周期性换相失败的防御方法 | |
| Wang et al. | Resilience enhancement strategy using microgrids in distribution network | |
| Qi et al. | Adaptive optimal PMU placement based on empirical observability gramian | |
| Weckesser et al. | Real-time remedial action against aperiodic small signal rotor angle instability | |
| Su et al. | Novel voltage-to-power sensitivity estimation for phasor measurement unit-unobservable distribution networks based on network equivalent | |
| Eidiani et al. | Determination of online DATC with uncertainty and state estimation | |
| Chandra et al. | Model-free angle stability assessment using wide area measurements | |
| CN113341210B (zh) | 母线暂态电压稳定性评估方法、系统、介质及计算设备 | |
| Leelaruji et al. | Voltage stability monitoring using sensitivities computed from synchronized phasor measurement data | |
| CN112086970B (zh) | 一种储能电站的配置方法、装置、电子设备和存储介质 | |
| Li et al. | Inclusion of classic HVDC links in a PMU-based state estimator | |
| Al-Omari et al. | Power system dynamic state estimation using smooth variable structure filter | |
| Chen et al. | Generalized injection shift factors and application to estimation of power flow transients | |
| CN105186440A (zh) | 基于机电暂态方法的继电保护定值整定方法 | |
| Peng et al. | A real-time enhanced thevenin equivalent parameter estimation method for PLL synchronization stability control in VSC | |
| Ghasemkhani et al. | Decentralized dynamic state estimation with missing and delayed PMU measurements | |
| Seo et al. | Determination of reactive power compensation considering large disturbances for power flow solvability in the korean power system | |
| Song et al. | Determination of load shedding for power-flow solvability using outage-continuation power flow (OCPF) | |
| KR102348138B1 (ko) | 온라인 발전기 고장 파급 시스템의 안정도 마진을 산정하기 위한 장치 및 그 방법 | |
| Duong et al. | Architecture and laboratory implementation of a testing platform for Wide Area Monitoring Systems | |
| Ali et al. | Online voltage instability detection of distribution systems for smart-grid applications | |
| Gubitz et al. | Online Detection and Quantification of Transient Instability using Lyapunov Exponents from PMU Data |
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 |