具体实施方式
下面结合附图对本发明中基于力矩分解的积分算法识别由发电机调速系统引发电力系统功率振荡的方法作进一步详细的说明。
调速系统是控制汽门或者导叶的运动,而控制进入原动机的蒸气流量或者水流流量,蒸气或者水流产生施加在发电机组轴系上的机械转矩,而影响原动机运行状态。机械转矩影响发电机转子运动,而对电力系统动态稳定产生影响。
如图1所示,Heffron-Philips模型中机械转矩增量的表达式为ΔTM=-GGOV(s)Δx,对于某一振荡频率f,将s=jΩ=j2πf代入,得:
ΔTM=-GGOV(s)Δx=-GGOV(jΩ)Δx=-GGOV(j2πf)Δx (1)
对于给定的输入信号Δx,考察调速系统GGOV(j2πf)的相位β,即可在Δδ-Δω座标系中画出ΔTM向量的方向,就可以确定ΔTM在Δω轴上的投影方向,也就可知调速系统对电力系统动态稳定的影响。
设调速系统的传递函数为GGOV(s),当调速系统以转速偏差Δω为输入信号时,调速系统产生的机械转矩增量:
ΔTM≈ΔPM=-GGOV(s)Δω (2)
对于所关心的振荡频率f,将s=jΩ=j2πf代入,得:
ΔTM=-GGOV(jΩ)Δω=-(A+jB)Δω=-α∠β·Δω (3)
其中:A为复数GGOV(jΩ)的实部,B为复数GGOV(jΩ)的虚部,α为复数GGOV(jΩ)的模值,β为复数GGOV(jΩ)的相角。
如图2所示,调速系统产生的机械转矩ΔTM=-α∠β·Δω在Δδ-Δω座标系中在Δω轴上的投影表示了其对动态稳定的影响。投影的正负表示了调速系统对阻尼影响的性质:在Δω轴上的投影为负,提供正阻尼,在Δω轴上的投影为正,提供负阻尼;投影的长短表示了调速系统对阻尼影响的大小,在Δω轴上的投影较长,则提供较大的(正或者负的)阻尼,在Δω轴上的投影较短,则提供较小的(正或者负的)阻尼。
本例中基于力矩分解的积分算法识别由发电机调速系统引发电力系统功率振荡的方法,包括发电机机械功率PM的拟合计算和发电机速度偏差Δω的确定,具体步骤如下:
将调速系统调门信号与原动机模型进行拟合计算获得发电机机械功率PM,从发电机机械功率PM中减去发电机机械功率平均值得到发电机机械功率偏差ΔPM;
根据发电机的端电压值、电流值计算发电机内电势,进而求取发电机内电势频率;
将经过高频滤波后的发电机内电势频率替代发电机速度偏差Δω,并与发电机机械功率偏差ΔPM进行积分运算,通过积分运算得出的积分量所在象限来判断发电机调速系统提供的阻尼为正或者为负,以便对于提供负阻尼的发电机及时采取措施,快速抑制振荡。
其中,发电机机械功率平均值通过下式获得:
式中,PM为发电机机械功率;PM(i)为第i时刻发电机机械功率值,PM均由调速系统调门信号与原动机模型经过拟合计算获得。
如图3所示,调速系统调门信号可以通过测试得到,汽轮机模型的参数可以由设计值得到,调门信号经过运算可以计算出汽轮机提供的机械功率PM。
当原动机模型为汽轮机模型时,发电机机械功率PM可以通过如下A、B、C中任一方法获得:
A)当汽轮机模型为一阶模型时,设汽轮机传递函数为:
B)当汽轮机模型为二阶模型时,设汽轮机传递函数为:
C)当汽轮机模型为三阶模型时,设汽轮机传递函数为:
图3中,PM为拟合的机械功率信号,PGV为调速系统调门信号,TCH为高压缸前汽室容积时间常数,FHP为高压缸功率比例系数,TRH为再热器容积时间常数,FIP为中压缸功率比例系数,TCO为低压连通管汽室容积时间常数,FLP为低压缸功率比例系数,上述各参数为汽轮机的出厂数据;s为微分算子。
如图4所示,调速系统调门信号可以通过向量测量系统(PMU)测试得到,水轮机模型的参数可以由设计值得到,调门信号经过运算可以计算出水轮机提供的机械功率PM。
当原动机模型为水轮机模型时,发电机机械功率PM通过如下水轮机传递函数求得:
图4中,PGV为调速系统调门信号,Tw为水流惯性时间常数,上述各参数可以从水力系统设计参数中获得;s为微分算子,PM为发电机机械功率。
通常情况下,调速系统调门信号PGV可以通过向量测量系统(PMU)测量得到,原动机模型可以采用设计参数,通过计算即可获取发电机近似的机械功率PM。
先将测得的发电机的端电压值和电流值代入下式求得发电机内电势:
EQ(t)=Ut+ItRa+jItXq (4)
式中,EQ(t)为t时刻的发电机内电势;Ut为t时刻的发电机端电压值;It为t时刻的发电机电流值;Ra为发电机定子电阻;j为虚部算子;Xq为发电机交轴电抗;其中,Ra、Xq均为发电机的出厂数据,可以直接获得。
然后,设置发电机内电势中仅含有基波分量,即:
式中,A、分别表示基波电压的幅值和初相角,t表示时刻;若用f0表示额定频率,Δf表示频差,f表示真实频率,三者之间存在如下关系:
f=f0+Δf (6)
首先测取发电机端电压与电流,进而计算出发电机内电势EQ(t)。对于内电势信号EQ(t),设信号中仅含有基波分量,若用f0表示额定频率,Δf表示频差,真实频率为f。由于真实频率未知,事先只能假定系统频率为额定值f0,对时间窗[0,T0]使用傅氏算法得到向量实部ζR0与虚部ζI0:
如果令:
则:
对于下一个时间窗上式依然成立,即:
因此可以计算出真实频率f:
为计算的更为准确,可以通过多个时间窗来计算真实频率f:
计算得到的频率f就是发电机内电势真实频率,可用来替代发电机速度偏差,但需要先进行高频滤波。高频滤波器通过下述方法来实现:所述滤波器为两个串联的双二次幂滤波器,该滤波器通过下式的传递函数来滤掉发电机内电势频率(即发电机速度偏差Δω)中的高频分量,
式中,s为积分算子,s2为积分算子的二次方;ω1、ω2、δ1、δ2、δ3、δ4为双二次幂滤波器的配置参数,为配置参数ω1、ω2的平方项。
用经过上述高频滤波后的发电机内电势频率替代发电机速度偏差Δω、并与发电机机械功率偏差ΔPM通过下式进行积分运算,得到积分量
式中,t0为初始时刻,t1为结束时刻,通常取振荡波形较好的7-10个周波内的数据进行积分计算。
如图5所示,通过计算出的积分量所在象限判断发电机调速系统对振荡频率提供的阻尼为正或者为负,从而进一步识别出对于某个振荡模式提供负阻尼的发电机调速系统,以便快速采取抑制振荡。判断方法为:当积分量位于第I、II象限时,调速系统提供的为负阻尼,当积分量位于第III、IV象限时,调速系统提供的为正阻尼。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求范围当中。