发明内容
本发明基于时域仿真方法,得到故障后系统的戴维南等值参数的变化,若系统戴维南等值电势持续下降,则系统功角失稳,反之则以系统的传输功率为目标函数,将故障后系统的负荷特性作为约束条件,通过拉格朗日因子法,构造修正目标函数,求取极值点,继而判断系统为电压失稳或感应电动机失稳。通过本发明,可判断故障后在不同负荷变化特性下的电力系统的失稳本质,可运用于电力系统稳定分析以及运行控制,利于系统运行、分析人员辨别故障后系统的失稳性质,采取有效的安全控制措施,提高电力系统的稳定运行水平。
具体而言,本发明提出了一种考虑负荷变化特性的电力系统稳定性分析方法,该方法包括以下步骤:
步骤A:基于时域仿真方法,求取故障后系统中需监测的负荷母线的系统戴维南等值参数;
步骤B:若系统的戴维南等值电势持续下降,则系统功角失稳,此方法的具体内容参见中国发明专利申请《一种基于动态戴维南等值判别电压失稳与功角失稳的方法》(申请号:200810119923.3);
步骤C:若系统的戴维南等值电势无持续下降趋势,则以系统的传输功率为目标函数,由已知的负荷节点的负荷模型,以故障后的负荷特性为约束条件,利用拉格朗日因子法,构造修正目标函数;
步骤D:求修正目标函数的驻点,即在约束条件下目标函数的可能极值点,再由可能的极值点中求得目标函数的最大值点以及相应目标函数的最大值Pmax;
步骤E:由已知的负荷模型,求解故障后的负荷需求Pn;
步骤F:比较故障后系统的负荷需求与最大传输功率间的关系,若负荷需求大于系统的最大传输功率,则为电压失稳,反之则为感应电动机失稳。
其中在所述步骤A中求取故障后负荷母线的系统戴维南等值参数的方法,是利用时域仿真方法,通过每一个计算步中生成的网络代数方程求得戴维南等值母线处的综合阻抗矩阵,再通过补偿法来求解任意一个负荷母线处的系统戴维南等值参数,
t时刻,在电力系统的时域仿真过程中,必须求解如下网络方程,以获得节点电压向量
其中,
为系统导纳矩阵;
为t时刻系统各个节点的注入电流向量;
为t时刻系统内各个节点的电压向量;
在节点i处单独注入单位电流,而所有其余节点的注入电流都等于0时,求解如下方程:
可以得到等值节点i处的综合阻抗矩阵ZiT,如下所示:
采用补偿法计算开路电压
即节点i开路时,相当于流经节点i处的负荷电流为0,可以在节点i处补偿一个与注入电流量
来求取系统节点电压的变化量;
由于 此时,流经阻抗的电流相当于 则
同时,基于前面所求得的综合阻抗矩阵ZiT,可以知道,节点i处的电压变化量为
根据叠加原理,节点i处的开路电压为:
其中,为t时刻暂态稳定计算得到的节点i处的电压;
联立求解式(4)和(6),求得
节点i处短路时,相当于在原有网络的基础上,在节点i处叠加一个注入电流量
根据叠加原理,此时节点i处的电压为:
其中,
为短路后节点i处的电压。而节点i处短路时,
即可求得:
根据叠加原理,可以求得节点i处的短路电流为:
这样,基于上面计算得到的开路电压
和短路电流
通过求解两者的比值,即可得到t时刻,节点i处的系统戴维南等值阻抗Z
t,iThev,如下所示:
同理,在故障发生后任意时刻,针对不同的负荷节点,重复上述步骤,可以计算得到任意一个负荷母线处随时间变化的系统戴维南等值电势
和戴维南等值阻抗Z
Thev;
其中在所述步骤B中若系统的戴维南等值电势持续下降,则系统功角失稳,如果负荷点的电压崩溃是由功角失稳而导致的,则负荷点的电压和戴维南电势将随功角差的增大而下降;
其中在所述步骤C中的目标函数是系统传输的有功功率,即等于负荷的有功功率,它是关于负荷参数t1、t2、t3……tn的n元函数,用f(t1、t2、t3……tn)表示;
其中在所述步骤B中的约束条件为该负荷母线对应的负荷变化规律,用
表示;
其中在所述步骤D中的修正目标函数的驻点可通过求解以下方程组得到:
得到的方程组(13)的解即为目标函数的可能极值点,由其中的可能极值点中再求得最大值点,代入目标函数f(t
1、t
2、t
3……t
n)即为在负荷变化规律
下系统的最大传输功率;
其中在所述步骤E中的负荷需求Pn是指在故障后负荷需求的有功功率,它与负荷特性相关;
其中在所述步骤F中的感应电动机失稳,其是区别于电压失稳的系统失稳形式,系统的故障可能导致感应电动机的转速超过临界转速,而无法再随着故障清除,系统的恢复而恢复;
其中所述步骤C、D、E中,在感应电动机并联恒阻抗负荷等值系统的负荷模型中,故障后,恒阻抗负荷的负荷需求Pnz将随着负荷母线电压的变化而变化,如公式(14)计算,其中U为负荷母线电压,ZR为恒阻抗负荷的阻抗,RR为恒阻抗负荷的电阻:
对于恒阻抗负荷,系统的最大传输功率随系统侧的戴维南等值参数的变化而变化,其供给恒阻抗负荷的最大传输功率为故障后恒阻抗负荷消耗的有功功率,即等于Pnz,因此对于恒阻抗负荷,负荷需求等于系统传输给恒阻抗负荷的最大传输功率,不存在电压失稳机制,则无电压稳定问题,因此只需关注感应电动机的电压稳定性;
对于感应电动机负荷,当所带负荷的特性为恒定机械转矩,故障后的感应电动机的负荷需求Pnm即为感应电动机的机械功率PT,可通过机械转矩Tt求解,其中,s为感应电动机滑差,ω为同步角速度,均以标幺值计算:
PT=Tt(1-s)ω (15)
当系统恢复稳定,有1-s≈1,ω≈1,因此对于恒定机械转矩的负荷,故障后感应电动机的负荷需求为初始稳态条件下的机械功率,如公式(16)所示:
Pnm=Tt(16)
感应电动机采用一阶模型,R
1,X
1为定子电阻和电抗,R
2,X
2为转子电阻和电抗,s为滑差,X
u为励磁电抗,Z
R为恒定阻抗负荷的等效阻抗,
Z
S系统的戴维南等值电势和阻抗,进一步戴维南等值化简后,得到由戴维南等值电势U
e串联等值阻抗R
1e+jX
1e以及转子阻抗
的等值电路,其中,U
e、R
1e、x
1e的计算如下:
R1e=Re((ZS//ZL+R1+jX1)//jXμ) (18)
X1e=Im ((ZS//ZL+R1+jX1)//jXμ) (19)
感应电动机的电磁功率Pm可通过下式计算:
假设转子参数R2,X2保持不变,仅有滑差s变化。求Pm对s的极值:
求解式(21),可得临界滑差scr:
将(22)代入(20),可得系统的最大传输功率Pmax:
系统的最大传输功率将随着故障后系统戴维南等值参数的变化而变化,通过故障后系统的负荷需求与系统最大的传输功率的关系曲线,由所述步骤F中的稳定判据可以得到系统的失稳本质。
本发明的有益效果是:本发明的方法可以区分系统受扰动后的失稳过程中,负荷节点电压大幅值下降时,失稳模式属于功角失稳、电压失稳还是感应电动机失稳。通过本发明,可判断故障后在不同负荷变化特性下的电力系统的失稳本质,可应用于电力系统稳定分析以及运行控制,有利于系统运行、分析人员辨别故障后系统的失稳性质,采取有效的安全控制措施,提高电力系统的稳定运行水平。