换相序技术应用于稳定控制的条件及最优控制策略

黄少锋1 李 慧1 李轶凡1 张月品2

(1. 新能源电力系统国家重点实验室(华北电力大学) 北京 102206 2. 北京四方继保自动化股份有限公司 北京 100085)

摘要 换相序技术是电力系统稳定控制的一种方案。在电力系统即将失稳的情况下,利用电力电子开关迅速进行相序切换,发电机侧ABC相换相序连接至系统侧cab相,进而使系统功角减小width=21.3,height=13.15,达到控制稳定的目的。基于能量函数法,分析换相序前后总能量的突变机理,得到换相序瞬间的能量变化量;在此基础上推导换相序技术的应用条件,通过换相序前后系统稳定性判断,给出换相序后系统稳定裕度增量,量化分析换相序技术提高系统暂态稳定性的程度,并以此为依据提出换相序技术的最优控制策略。基于Matlab的仿真结果表明,换相序最优控制策略能够实现稳定控制的最优效果,获得最大稳定裕度,并验证了应用条件的正确性。

关键词:换相序技术 稳定控制 能量函数分析 暂态稳定性 应用条件 最优控制策略

0 引言

稳定控制是交流系统最关注的课题之一,伴随着电力规模的逐渐扩大,系统结构和模型日益复杂,电力系统的运行也更加灵活多变,对电力系统稳定控制提出了更大的挑战[1-4]

电力系统中许多破坏性事故都源于大扰动下发电机的功角失稳,通常以功角在width=25.05,height=11.9范围内变化来表征。大扰动后电磁功率和机械功率的不平衡是导致系统暂态稳定破坏的主要原因。长期以来,电力系统采用了许多提高系统暂态稳定性的措施,例如故障的快速切除和自动重合闸装置的应用,采用强行励磁、快关汽门控制、紧急切机等[5-10]。这些措施在一定条件下起到了改善系统稳定性的作用,但也牺牲了电网完整性。

文献[11]提出了一种利用断路器实现width=20.05,height=13.75相位旋转的稳定控制方法,但传统断路器分合闸过程较慢,在断路器分闸、合闸过程中功角会发生很大改变。近年来,随着电力电子技术的广泛应用,基于电力电子元件的交流固态断路器表现出优异性能[12-14],分合闸过程实现在20ms以内,对于振荡第一摆,20ms对应的功角偏差不大。同时,量测信息的实时性和全局性[15-16]为电力系统在线稳定控制提供了技术支撑[17]。文献[18]提出了一种稳定控制措施——换相序技术(Phase Sequence Exchange Technology, PSET),利用电力电子开关迅速进行相序切换,从而使系统功角减小width=21.3,height=13.15,并基于等面积定则分析了换相序技术的基本原理。该措施有利于电网结构的完整性,提高了电力系统稳定性。文献[19]计算了换相序过程产生的冲击电流及冲击转矩,研究了系统能够承受该冲击的条件,并提出了固态断路器分相分合闸的换相序方法,以减小系统冲击。文献[20]基于拓展等面积法,在多机系统中提出了一种换相序控制方法,并验证了换相序技术在多机系统中的有效性。

本文基于文献[18],分析了换相序前后系统暂态稳定特性,在此基础上进行暂态能量函数分析,得到了换相序瞬间的能量变化量及换相序技术使系统趋于稳定的必要条件。基于此,推导了换相序技术的应用条件,以稳定裕度增量为依据,量化分析了换相序技术提高系统暂态稳定性的程度,并提出了换相序的最优控制策略。

1 换相序前后的暂态稳定特性分析

1.1 换相序前后功角转换

根据参考文献[18],基于单机无穷大系统进行换相序分析,其相量图如图1所示。无穷大系统侧三相电压为width=55.1,height=15.65,发电机三相电动势为width=54.45,height=15.65。以width=11.25,height=11.25相为例,系统功角width=10,height=11.9width=15.05,height=15.65width=13.75,height=15.65夹角。当系统即将失稳,功角为width=50.1,height=16.9时进行换相序操作,利用电力电子开关对相序进行快速切换[19],即发电机侧width=24.4,height=11.9相与系统侧width=18.15,height=13.15相断开后对应连接至系统侧width=18.15,height=13.15相,其示意图如图1d、图1e所示。换相序瞬间,原有电动势相量并未发生突变,仅将原发电机电动势width=54.45,height=15.65标号改为系统侧标号顺序width=58.25,height=15.65,于是,换相序结束瞬间功角变为width=11.9,height=15.05,如图1c所示。显然,换相序后系统功角减小了width=21.3,height=13.15,以width=11.9,height=13.15表示,从而有

width=87.65,height=28.15 (1)

式中,width=10,height=11.9width=11.9,height=13.15分别为换相序前、后功角;下标width=9.4,height=11.9代表换相序瞬间。

width=222.75,height=156

图1 换相序技术相量图及示意图

Fig.1 PSET phasor diagrams and schematic diagrams

下面从width=18.25,height=15.05坐标系统阐述换相序前后功角转换,首先分析换相序前后坐标系统的转换。换相序前,发电机width=24.7,height=11.8坐标系统经Park变换转换为width=18.25,height=15.05坐标系统,width=10.75,height=10.75轴与width=9.65,height=12.9轴夹角为width=9.65,height=11.8;换相序后,发电机切换为width=24.7,height=11.8相序,经Park变换转换为width=25.8,height=15.05坐标系统,width=10.75,height=10.75轴与width=10.75,height=12.9轴夹角仍为width=9.65,height=11.8。由于width=24.7,height=11.8坐标系统超前width=24.7,height=11.8坐标系统width=21.5,height=12.9,因此,width=25.8,height=15.05坐标系统超前原width=18.25,height=15.05坐标系统width=21.5,height=12.9,如图2所示。

进一步通过Park变换推导width=25.8,height=15.05width=18.25,height=15.05坐标系统下电动势相量的关系。换相序前瞬间,发电机电动势width=15.05,height=16.1width=33.3,height=16.1经Park变换转换为width=50.5,height=17.2,从而有

width=63.4,height=54.8(2)

width=212.25,height=111.75

图2 Park变换示意图

Fig.2 dq0 transform diagrams

式中,width=50.5,height=17.2width=18.25,height=15.05坐标系统下的发电机电动势;width=10.75,height=10.75为Park变换矩阵,即

width=195.6,height=91.35(3)

式中,width=9.65,height=11.8width=10.75,height=10.75轴与width=9.65,height=12.9轴夹角。

换相序结束瞬间,发电机电动势width=58.05,height=16.1经Park变换转换为width=54.8,height=17.2,从而有

width=66.65,height=48.35(4)

式中,width=54.8,height=17.2width=25.8,height=15.05坐标系统下的发电机电动势;P矩阵中,width=9.65,height=11.8width=10.75,height=10.75轴与width=10.75,height=12.9轴夹角。

比较图1b、图1c可以发现,width=54.8,height=16.1width=16.1,height=16.1width=36.55,height=16.1关系为

width=74.15,height=48.35(5)

将式(5)代入式(4)得

width=81.65,height=48.35(6)

再将式(2)代入式(6)得

width=72,height=48.35(7)

将图2与式(7)结合分析,图2中,换相序结束瞬间,转子位置由于惯性并未发生突变,但由于width=25.8,height=15.05坐标系统超前原width=18.25,height=15.05坐标系统width=21.5,height=12.9,使得式(7)中width=54.8,height=17.2滞后width=50.5,height=17.2width=21.5,height=12.9,也就是说,换相序结束瞬间,Park变换后的发电机电动势相角减小了width=21.5,height=12.9,因此,功角减小了width=21.5,height=12.9

以上从width=24.7,height=11.8width=18.25,height=15.05两种坐标系统阐明了换相序前后功角转换,如式(1)所示。下面分析换相序前后发电机功角特性。

1.2 换相序前后功角特性分析

假设发电机采用经典二阶模型,此时可忽略原动机、调速器、励磁系统的影响,于是,机械功率width=13.95,height=15.05保持不变,故障切除后电磁功率width=18.25,height=17.2[21]

width=115,height=30.1 (8)

式中,width=9.65,height=11.8为功角,即发电机等效电动势width=26.85,height=12.9与无穷大母线width=25.8,height=13.95之间的角度;width=17.2,height=16.1为故障切除后系统总电抗;width=65.55,height=16.1为故障切除后电磁功率最大值。

换相序后,width=13.95,height=15.05保持不变。此外,电动势幅值width=12.9,height=11.8width=11.8,height=11.8width=17.2,height=16.1不变,功角变为width=11.8,height=12.9,于是,换相序后电磁功率width=20.4,height=17.2

width=122.5,height=30.1 (9)

将式(1)代入式(9)可得

width=97.8,height=29 (10)

将式(8)~式(10)绘制于图3,以便阐述换相序过程。图中,上标不带“′”的变量对应width=20.4,height=11.8坐标系,上标带“′”的变量对应width=24.7,height=12.9坐标系。在width=11.8,height=15.05处换相序时,系统功角减小width=21.5,height=12.9,突变为width=11.8,height=15.05,且电磁功率从width=18.25,height=17.2曲线上的width=9.65,height=12.9点突变为width=20.4,height=17.2曲线上width=10.75,height=12.9点。此后,沿width=20.4,height=17.2曲线运行,功角为width=54.8,height=12.9

width=185.25,height=99

图3 换相序前后功角特性曲线

Fig.3 Power-angle curves before and after phase sequence exchanging

1.3 换相序前后平衡点的关系

换相序前,width=11.8,height=15.05为故障切除后稳定平衡点,width=11.8,height=15.05为不稳定平衡点,引用文献[22]对平衡点的描述,从而有

width=84.9,height=63.4(11)

设换相序后平衡点为width=11.8,height=15.05width=11.8,height=15.05,由于平衡点处满足width=40.85,height=17.2,将式(9)代入其中,解得

width=83.8,height=63.4(12)

由于width=115,height=17.2,则width=11.8,height=15.05为换相序后稳定平衡点,记为width=11.8,height=15.05width=11.8,height=15.05为换相序后不稳定平衡点,记为width=11.8,height=15.05。再结合式(11)、式(12)可得

width=38.7,height=35.45 (13)

式(13)表明,换相序后系统平衡点保持不变,如图3所示。

因此,换相序后功角减小width=21.5,height=12.9,电磁功率变为width=20.4,height=17.2,但是平衡点保持不变,如式(1)、式(9)、式(13)所示。此结论是下文能量函数分析的基础。

2 换相序技术的能量函数分析

基于暂态能量函数法,给出换相序前、后系统总能量,进而得到换相序瞬间的能量变化量。

2.1 换相序前能量函数分析

依据经典能量函数法[23],定义系统动能width=11.8,height=15.05

width=66.65,height=26.85 (14)

式中,width=10.75,height=9.65为转速差;width=13.95,height=10.75为发电机惯性时间常数。

以稳定平稳点width=11.8,height=15.05为势能参考点,定义系统势能width=11.8,height=16.1

width=234.25,height=22.55(15)

式中,右侧第一项为系统电抗中存储的磁能;第二项为转子的位置能量。

故障切除后至换相序前,系统总能量width=10.75,height=11.8

width=214.95,height=45.15(16)

当系统功角为width=11.8,height=15.05,转速差为width=15.05,height=16.1时进行换相序操作,则换相序前瞬间动能为width=33.3,height=18.25、势能为width=31.15,height=17.2,总能量width=11.8,height=15.05

width=226.75,height=47.3

2.2 换相序后能量函数分析

在刚完成换相序操作的瞬间,功角及转速差为width=11.8,height=15.05width=15.05,height=16.1。假设换相序过程迅速,对发电机转速影响较小,则

width=36.55,height=16.1 (18)

此外,width=11.8,height=15.05width=11.8,height=15.05关系如式(1)所示。此时,利用换相序前后瞬间功角及转速差的关系,即可分析换相序前后各能量变化情况。

考虑式(18)的关系,则换相序结束瞬间动能width=33.3,height=18.25

width=161.2,height=26.85 (19)

因此,换相序瞬间动能保持不变,则能量变化只取决于势能变化。

换相序后,以稳定平衡点width=11.8,height=15.05为势能参考点,则系统势能width=30.1,height=17.2

width=203.1,height=40.85

将式(1)、式(13)代入式(20)得

width=224.6,height=61.25

式(21)表明,在width=22.55,height=16.1坐标系下,width=12.9,height=16.1相当于width=11.8,height=16.1右向平移width=21.5,height=12.9,由此绘制width=11.8,height=16.1width=12.9,height=16.1曲线于图4b,以便后续对比与分析。

width=31.15,height=15.05width=29,height=15.05代入式(21),即可得到换相序结束瞬间势能width=31.15,height=17.2

width=200.25,height=216

图4 width=44.9,height=14.95时换相序前后能量变化图

Fig.4 Energy change diagrams before and after phase sequence exchanging whenwidth=45.15,height=15.05

width=214.95,height=61.25

于是,换相序结束瞬间系统总能量width=12.9,height=15.05

width=212.8,height=78.45

2.3 换相序瞬间能量变化量

根据上述分析,可得换相序瞬间能量的变化量DV

width=209.55,height=38.7

将式(19)的结论代入式(24),可得

width=205.25,height=77.35

式(25)表明,总能量的变化量DV亦为势能的变化量。此外,width=36.55,height=16.1width=11.8,height=15.05的取值决定了能量变化量DV的正负。当width=31.15,height=11.8时,换相序后系统总能量增大;当width=31.15,height=11.8时,换相序后系统总能量减小。

从物理意义角度进一步分析DV中各项含义,式(25)中,width=40.85,height=15.05为转子位置势能的变化。虽然换相序瞬间转子位置不变,但参考坐标系统的改变使转子角滞后width=21.5,height=12.9,因而转子位置势能增加了width=32.25,height=15.05width=83.8,height=18.25为电抗中存储磁能的变化量,当该项大于0时,其计算值即为换相序瞬间磁性势能的减小量。当width=120.35,height=18.25,即width=31.15,height=11.8时,势能减小量大于增加量,因而消减了系统能量。因此,不采取其他措施,仅通过一次或多次换相序控制使系统趋于稳定的必要条件为

width=31.15,height=11.8 (26)

将换相序前后能量变化过程绘制出来,如图4所示。图4a中,轨迹width=13.95,height=15.05width=13.95,height=15.05分别为换相序前、后相轨迹,轨迹width=13.95,height=15.05width=9.65,height=15.05点对应换相序前瞬间功角width=11.8,height=15.05、转速差width=15.05,height=16.1,对应能量为换相序前瞬间总能量width=11.8,height=15.05,如式(17)所示。由于相轨迹为等能量曲线,width=9.65,height=15.05点总能量可直观表示为等能量曲线width=13.95,height=15.05。轨迹width=13.95,height=15.05width=10.75,height=15.05点对应换相序结束瞬间功角width=60.2,height=15.05、转速差width=36.55,height=16.1,对应能量为换相序结束瞬间总能量width=12.9,height=15.05,如式(23)所示,则width=10.75,height=15.05点总能量可直观表示为等能量曲线width=13.95,height=15.05,这样,width=18.25,height=15.05间总能量差值即为等能量曲线width=13.95,height=15.05width=13.95,height=15.05间的总能量差值。图4b中,根据式(15)、式(21)绘制换相序前后势能变化曲线,如曲线③、④所示。width=27.95,height=15.05点分别为换相序前、后瞬间势能width=31.15,height=17.2width=31.15,height=17.2。由图4可以看出,在width=11.8,height=15.05处的换相序瞬间,系统总能量由width=9.65,height=15.05突变至width=10.75,height=15.05点,势能由width=10.75,height=15.05突变至width=10.75,height=15.05点。根据式(24)、式(25),width=18.25,height=15.05间总能量差值(亦为等能量曲线width=13.95,height=15.05width=13.95,height=15.05间的能量差值)与width=20.4,height=15.05间势能差值相等,均为能量变化量DV。此时,width=31.15,height=11.8,换相序减小了系统能量。

因此,换相序技术的关键在于控制瞬间能量减小,其能量变化量DV如式(25)所示。

3 换相序技术应用条件及最优控制策略分析

基于上述结论,推导了换相序技术的应用条件,并以稳定裕度增量为依据,提出了换相序技术的最优控制策略。

3.1 换相序技术应用条件分析

将式(25)代入式(26)的稳定必要条件,可得

width=107.45,height=31.15 (27)

再结合式(11)的width=11.8,height=15.05表达式,则式(27)整理为

width=122.5,height=33.3(28)

式中,width=11.8,height=15.05为换相序前瞬间功角。

因此,式(28)即为基于能量函数法的换相序技术应用条件。换句话说,width=11.8,height=15.05满足式(28)条件后,则DV满足式(26),在width=29,height=15.05处换相序才能够提高系统暂态稳定性。

3.2 换相序技术提高稳定裕度分析

为了进一步量化换相序技术对系统稳定性的提升程度,从稳定裕度增量角度进行对比分析。通过换相序前后系统临界能量分析,判断系统稳定性。

3.2.1 换相序前后系统稳定性判断

1)判断换相序前系统是否稳定

将故障切除时刻的功角width=11.8,height=15.05及转速差width=12.9,height=15.05代入式(16)即可得到对应的能量为width=61.25,height=17.2。此外,不采取换相序措施时,以故障切除后不稳定平衡点width=11.8,height=15.05处的势能作为临界能量width=12.9,height=15.05,即

width=165.5,height=17.2(29)

于是,比较width=10.75,height=15.05width=12.9,height=15.05,即可判定系统是否稳定。也就是说,若width=35.45,height=15.05,则系统稳定;反之若width=31.15,height=15.05,则系统失稳,需采取换相序控制措施,如图5中轨迹width=13.95,height=15.05所示。

width=189.75,height=78

图5 width=40.85,height=13.95时换相序前后相轨迹

Fig.5 Phase trajectory diagram before and after phase sequence exchanging whenwidth=40.85,height=13.95

2)判断换相序后系统是否稳定

为了以故障切除时刻能量width=10.75,height=15.05判断换相序后系统是否稳定,需分析width=10.75,height=15.05对应的临界能量。由于直接求解较为困难,考虑到width=10.75,height=15.05与换相序前瞬间总能量width=11.8,height=15.05相等,故将式(24)转换为

width=54.8,height=15.05 (30)

因此,上述问题转换为求解换相序后总能量width=12.9,height=15.05对应的临界能量,在此基础上加width=18.25,height=11.8,即可得到width=10.75,height=15.05对应的临界能量。

由图4b可知,在不稳定平衡点width=11.8,height=15.05处,换相序后的势能取得最大值。于是,以width=11.8,height=15.05处的势能作为width=12.9,height=15.05对应的临界能量width=12.9,height=15.05,即

width=168.7,height=17.2(31)

将式(13)代入式(31)可得

width=127.9,height=17.2 (32)

因此,width=10.75,height=15.05对应的临界能量width=36.55,height=15.05

width=214.95,height=61.25

于是,比较故障切除时刻能量width=10.75,height=15.05与其对应的临界能量width=36.55,height=15.05,即可判定该次换相序后系统是否稳定。也就是说,若width=59.1,height=15.05,则该次换相序能够使系统稳定,如图5中width=10.75,height=15.05width=36.55,height=15.05、轨迹width=13.95,height=15.05所示;反之若width=54.8,height=15.05,则该次换相序后系统仍失稳,可进行多次换相序控制,逐次减小系统能量直至系统稳定。

3.2.2 稳定裕度增量分析

根据上述分析,以临界能量与故障切除时刻能量的差值作为系统稳定裕度的定量描述,并进行规格化处理。于是,不采取换相序措施时,稳定裕度width=18.25,height=15.05[23]

width=60.2,height=30.1 (34)

式中,Vkc为故障切除时的动能,width=61.25,height=16.1

同理,采取换相序措施后,稳定裕度width=20.4,height=15.05

width=81.65,height=30.1 (35)

稳定裕度width=18.25,height=15.05width=21.5,height=15.05小于0时,系统失稳;等于0时,为临界稳定;大于0时,系统稳定,且数值越大,系统稳定程度越高。

比较式(34)、式(35)可以发现,换相序后系统稳定裕度增量width=24.7,height=16.1

width=141.85,height=30.1 (36)

再将式(29)、式(32)代入式(36),可得

width=50.5,height=29 (37)

式(37)表明,增加的稳定裕度width=35.45,height=15.05由能量变化量DV决定,DV满足式(26)的必要条件,则width=38.7,height=16.1。因此,换相序提高系统稳定裕度的本质为控制瞬间能量减小。

3.3 换相序技术最优控制策略分析

为了实现换相序最优效果,获得最大的稳定裕度,需分析换相序最优控制策略。

由式(37)可知,width=18.25,height=11.8越大,换相序消减的能量越大,稳定裕度增量width=24.7,height=16.1越大,越有利于系统稳定。因此,若能推导出DV最大时对应的换相序功角值,即可得到换相序最优控制策略。为此,令式(25)的DV表达式中width=36.55,height=16.1保持不变,求得DV关于换相序前瞬间功角width=11.8,height=15.05的导数为

width=110.7,height=30.1 (38)

为求得DV驻点,令式(38)等于0且width=50.5,height=17.2,得到

width=35.45,height=27.95 (39)

将式(39)与图4b结合分析可以确定:当width=40.85,height=15.05时换相序,DV为正最大值,从而达到换相序最优效果。此时,稳定裕度增量width=24.7,height=16.1最大,其数值为

width=99.95,height=40.85 (40)

因此,最优控制策略为width=40.85,height=15.05时换相序。将width=40.85,height=15.05代入式(28),即可得到最优控制策略对应的应用条件为

width=43,height=15.05 (41)

也就是说,当width=11.8,height=15.05满足式(41)时,最优控制策略下稳定裕度增量最大值width=30.1,height=16.1必定大于0。

综上所述,应用于稳定控制的换相序技术能够极大地提高系统暂态稳定性,但是需满足式(28)的应用条件。此外,最优控制策略为功角等于width=21.5,height=12.9时换相序,从而获得最大稳定裕度。

4 仿真验证

算例采用文献[24]中的单机无穷大系统。初始功角width=65.55,height=15.05width=25.8,height=11.8width=21.5,height=11.8时刻发生故障,故障期间电磁功率最大值width=51.6,height=17.2width=38.7,height=11.8故障切除,故障切除后width=46.2,height=16.1width=66.65,height=15.05。其他参数为width=47.3,height=11.8width=27.95,height=15.05,阻尼width=26.85,height=11.8。为验证上述理想二阶模型下理论的准确性,本文采用数值积分法进行仿真验证,积分步长为0.001s,并根据上述能量函数公式计算各状态下的能量。

1)是否采取换相序最优控制策略的相轨迹对比如图6a所示。故障切除时刻能量width=45.15,height=15.05、临界能量width=47.3,height=15.05,由于width=31.15,height=15.05,系统沿轨迹width=13.95,height=15.05失稳;若在功角width=40.85,height=15.05处换相序,由式(23)、式(30)计算得,换相序后能量突变为width=45.15,height=15.05,能量减小了width=87.05,height=15.05。此外,由式(33)计算得换相序的临界能量width=45.15,height=15.050.585,由于width=54.8,height=15.05,则系统沿轨迹width=13.95,height=15.05趋于稳定。

width=177.75,height=264.75

图6 是否采取换相序最优控制策略对比

Fig.6 The comparison diagrams of whether optimal control strategy of PSET is implemented

上述情况对应的功角摇摆曲线如图6b所示,不采取换相序措施时功角为width=9.65,height=11.8,换相序后功角为width=54.8,height=12.9。由式(34)~式(37)可得,换相序控制下,暂态稳定裕度由-0.352提高到0.244,换相序增加的稳定裕度为0.596。由此可以知道,换相序技术使系统由失稳趋于稳定,极大地提高了系统稳定裕度。

2)当采取最优控制策略width=40.85,height=15.05时,式(28)的应用条件可整理为式(41)的形式。在上述仿真中width=66.65,height=15.05,满足式(41)的应用条件。为了验证应用条件的正确性,进一步验证式(41)对应的临界情况,即width=98.85,height=15.05。此时,仿真参数修改为width=94.55,height=16.1,对应的临界情况相轨迹如图7所示。经计算,式(17)的换相序前能量width=11.8,height=15.05与式(23)换相序后能量width=12.9,height=15.05相等,均为0.437,于是,换相序瞬间能量变化量width=69.85,height=15.05,因而换相序后系统仍然失稳。对比图6a、图7可知,只有满足换相序技术应用条件,一次或多次换相序控制才能够减小系统能量,使系统趋于稳定。

width=171.75,height=120

图7 临界情况仿真曲线

Fig.7 Simulation curves in the critical case

图8为width=65.55,height=15.05width=40.85,height=15.05时不满足式(41)应用条件的仿真曲线。仿真结果表明,虽然换相序后系统仍失稳,但换相序技术起延缓系统失稳作用。

width=228,height=120.75

图8 不满足应用条件时仿真曲线

Fig.8 The simulation curves when the application condition is not satisfied

3)图9、表1为不同功角处换相序的仿真结果。由式(25)、式(34)~式(37)计算得不同功角处换相序的能量变化量DV、不换相序的稳定裕度width=18.25,height=15.05、换相序后稳定裕度width=21.5,height=15.05、稳定裕度增量width=24.7,height=16.1,见表1所示。分析图9、表1可知,功角为width=40.85,height=15.05时换相序,系统剩余能量最小,减小的能量DV最大,因而稳定裕度增量width=24.7,height=16.1最大,换相序效果最优。

width=174.75,height=125.25

图9 不同功角下换相序的仿真曲线

Fig.9 Simulation curves of phase sequence exchanging at different power angles

表1 不同功角处换相序的仿真结果

Tab.1 Simulation results of phase sequence exchanging at different power angles

/radDV 0.2088-0.3520.1580.51 0.23530.2230.575 0.24390.2440.596 0.23490.2220.574 0.20660.1530.505

此外,还验证了多种情况,均证明了上述理论的正确性,不再赘述。

5 结论

针对换相序稳定控制技术,本文分析了换相序前后系统的暂态稳定特性,采用能量函数法分析换相序前后能量变化,得到了换相序瞬间的能量变化量,并得到了如下结论:

1)从能量函数角度出发,推导了换相序技术的应用条件。

2)功角为width=21.5,height=12.9时进行换相序控制,能够获得最大稳定裕度,实现稳定控制的最优效果。

3)通过换相序前后系统临界能量分析,推导了换相序技术对暂态稳定性的提高程度(即稳定裕度增量)。

理论分析和仿真结果均表明,本文的最优控制策略能够实现稳定控制的最优效果,并验证了应用条件的正确性。

参考文献

[1]唐飞, 刘扬, 施浩波, 等. 一种考虑风电场并网的大电网快速主动解列策略[J]. 电工技术学报, 2019, 34(10): 2092-2101. Tang Fei, Liu Yang, Shi Haobo, et al. A fast active islanding strategy for large power grid considering wind farm integration[J]. Transactions of China Electrotechnical Society, 2019, 34(10): 2092-2101.

[2]陈鹏远, 黎灿兵, 周斌, 等. 异步互联电网柔性直流输电紧急功率支援与动态区域控制偏差协调控制策略[J]. 电工技术学报, 2019, 34(14): 3025-3034. Chen Pengyuan, Li Canbing, Zhou Bin, et al. VSC-HVDC emergency power support and dynamic area control error coordinated control strategy for improving the stability of asynchronous interconnected power grids[J]. Transactions of China Electrotechnical Society, 2019, 34(14): 3025-3034.

[3]林旻威, 温步瀛. 大规模风电接入对电力系统暂态稳定性影响研究综述[J]. 电气技术, 2017, 18(4): 1-8. Lin Minwei, Wen Buying. The overview of influence of large scale wind generation on transient stability of power system[J]. Electrical Engineering, 2017, 18(4): 1-8.

[4]辛建波, 王玉麟, 舒展, 等. 特高压交直流接入对江西电网暂态稳定的影响分析[J]. 电力系统保护与控制, 2019, 47(8): 71-79. Xin Jianbo, Wang Yulin, Shu Zhan, et al. Transient stability impact analysis of UHV AC/DC access to Jiangxi power grid[J]. Power System Protection and Control, 2019, 47(8): 71-79.

[5]顾卓远, 汤涌, 孙华东, 等. 基于响应的电力系统安全稳定综合防御技术[J]. 中国电机工程学报, 2019, 39(4): 943-953. Gu Zhuoyuan, Tang Yong, Sun Huadong, et al. Study on framework of comprehensive defense architecture for power system security and stability[J]. Proceedings of the CSEE, 2019, 39(4): 943-953.

[6]石访, 张林林, 胡熊伟, 等. 基于多属性决策树的电网暂态稳定规则提取方法[J]. 电工技术学报, 2019, 34(11): 2364-2374. Shi Fang, Zhang Linlin, Hu Xiongwei, et al. Power system transient stability rules extraction based on multi-attribute decision tree[J]. Transactions of China Electrotechnical Society, 2019, 34(11): 2364-2374.

[7]刘九良, 王彤, 朱劭璇, 等. 计及保护信息的电力系统暂态稳定裕度解析算法[J]. 电工技术学报, 2020, 35(3): 542-552. Liu Jiuliang, Wang Tong, Zhu shaoxuan, et al. An analytic method for power system transient stability margin considering protection information[J]. Transactions of China Electrotechnical Society, 2020, 35(3): 542-552.

[8]孙黎霞, 白景涛, 周照宇, 等. 基于双向长短期记忆网络的电力系统暂态稳定评估[J]. 电力系统自动化, 2020, 44(13): 64-72. Sun Lixia, Bai Jingtao, Zhou Zhaoyu, et al. Transient stability assessment of power system based on bi-directional long-short-term memory network[J]. Automation of Electric Power Systems, 2020, 44(13): 64-72.

[9]韦肖燕, 李欣然, 钱军, 等. 采用储能电源辅助的暂态稳定紧急控制方法[J]. 电工技术学报, 2017, 32(18): 286-300. Wei Xiaoyan, Li Xinran, Qian Jun, et al. Power system transient stability emergency control method assisted by energy storage[J]. Transactions of China Electrotechnical Society, 2017, 32(18): 286-300.

[10]陈学通, 凌超, 薛峰, 等. 一种基于贪心算法的紧急控制策略优化搜索方法[J]. 电力系统保护与控制, 2017, 45(23): 74-81. Chen Xuetong, Ling Chao, Xue Feng, et al. A fast method of searching optimal emergency control strategy based on greedy algorithm[J]. Power System Protection and Control, 2017, 45(23): 74-81.

[11]Cresap R L, Taylor C W, Kreipe M J. Transient stability enhancement by 120-degree phase rotation[J]. IEEE Transactions on Power Apparatus and Systems, 1981, 100(2): 745-753.

[12]Song S M, Kim J Y, Choi S S, et al. New simple-structured ac solid-state circuit breaker[J]. IEEE Transactions on Industrial Electronics, 2018, 65(11): 8455-8463.

[13]卢其威, 高志宣, 滕尚甫, 等. 基于MOSFET的限流式固态断路器及其过电压抑制[J]. 电工技术学报, 2017, 32(24): 42-52. Lu Qiwei, Gao Zhixuan, Teng Shangfu, et al. Current limiting solid state circuit breaker based on MOSEFT and its over voltage suppression[J]. Transactions of China Electrotechnical Society, 2017, 32(24): 42-52.

[14]Meyer C, Schroder S, De Doncker R W. Solid-state circuit breakers and current limiters for medium-voltage systems having distributed power systems[J]. IEEE Transactions on Power Electronics, 2004, 19(5): 1333-1340.

[15]许苏迪, 刘灏, 毕天姝, 等. 一种适用于同步相量测量装置校准器的高精度相量测量方法[J]. 电工技术学报, 2020, 35(2): 372-382. Xu Sudi, Liu Hao, Bi Tianshu, et al. A high accuracy phasor estimation algorithm for phasor measurement units calibrator[J]. Transactions of China Electrotechnical Society, 2020, 35(2): 372-382.

[16]罗深增, 李银红, 石东源. 广域测量系统可观性概率评估及其在PMU优化配置中的应用[J]. 电工技术学报, 2018, 33(8): 1844-1853. Luo Shenzeng, Li Yinhong, Shi Dongyuan. Wide area monitoring system observability probabilistic evaluation and it’s application in optimal PMU placement[J]. Transactions of China Electrotechnical Society, 2018, 33(8): 1844-1853.

[17]赵晋泉, 邓晖, 吴小辰, 等. 基于广域响应的电力系统暂态稳定控制技术评述[J]. 电力系统保护与控制, 2016, 44(5): 1-9. Zhao Jinquan, Deng Hui, Wu Xiaochen, et al. Review on power system transient stability control technologies based on PMU/WAMS[J]. Power System Protection and Control, 2016, 44(5): 1-9.

[18]Huang Shaofeng, Li Yifan, Li Hui, et al. A new technology applied to power system stability control: phase sequence exchange technology[J]. IEEE Access, 2019, 7: 93002-93009.

[19]Li Yifan, Huang Shaofeng, Li Hui, et al. A new phase sequence exchanging control method for reducing impulse current and voltage[J]. IEEE Access, 2019, 7: 164734-164745.

[20]Li Yifan, Huang Shaofeng, Li Hui, et al. Application of phase sequence exchange in emergency control of a multi-machine system[J]. International Journal of Electrical Power & Energy Systems, 2020, DOI: 10.1016/j.ijepes.2020.106136.

[21]李光琦. 电力系统暂态分析[M]. 北京: 中国电力出版社, 2007.

[22]闵勇, 陈磊, 姜齐荣. 电力系统稳定分析[M]. 北京: 清华大学出版社, 2016.

[23]倪以信, 陈寿孙, 张宝霖. 动态电力系统的理论和分析[M]. 北京: 清华大学出版社, 2002.

[24]刘笙, 陈陈. 电力系统暂态稳定的能量函数分析[M]. 北京: 科学出版社, 2014.

The Condition of Phase Sequence Exchange Technology Applied to Stability Control and Optimal Control Strategy

Huang Shaofeng1 Li Hui1 Li Yifan1 Zhang Yuepin2

(1. State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources North China Electric Power University Beijing 102206 China 2. Beijing Sifang Automation Co. Ltd Beijing 100085 China)

Abstract Phase sequence exchange technology (PSET) is a scheme of power system stability control. When the system is about to lose stability, the phase sequence is quickly exchanged by the power electronic switch. The ABC phase of generator side is disconnected and then connected to cab phase of system side correspondingly, thereby reducing the power angle bywidth=21.3,height=13.15and achieving the purpose of stability control. Based on the energy function method, the mutation mechanism of total energy before and after phase sequence exchanging was analyzed and the mutational energy at the instant of phase sequence exchanging was obtained. On this basis, the application condition of PSET was derived. By judging the system stability before and after phase sequence exchanging, the stability margin increment after phase sequence exchanging was given, thus analyzing the degree of improving the power system transient stability by PSET quantitatively. Based on this, the optimal control strategy of PSET was proposed. Finally, simulation results based on Matlab show that the optimal control strategy of PSET can achieve the optimal effect of stable control and obtain the maximum stability margin. Besides, the correctness of application condition is verified.

keywords:Phase sequence exchange technology, stability control, energy function analysis, transient stability, application condition, optimal control strategy

中图分类号:TM712

DOI:10.19595/j.cnki.1000-6753.tces.200457

中央高校基本科研业务费资助项目(2019QN107)。

收稿日期 2020-05-07

改稿日期 2020-08-01

作者简介

黄少锋 男,1958年生,教授,博士生导师,研究方向为电力系统继电保护。E-mail:huangsf@sf-auto.com

李 慧 女,1994年生,博士研究生,研究方向为电力系统稳定与控制。E-mail:student089@163.com(通信作者)

(编辑 赫蕾)