一种计及死区效应的电力系统最低频率快速计算方法
技术领域
本发明应用于电力系统仿真领域,具体涉及一种计及死区效应的电力系统最低频率快速计算方法。
背景技术
随着我国电网规模的逐渐扩大,频率稳定越来越成为电力规划部门关心的重要问题之一。我国电力系统存在远距离、大规模送电的情况,一旦大容量输电通道被破坏、大容量机组被切除等,很容易导致系统频率的大幅下降,甚至引发频率崩溃。因此,快速计算扰动后系统的最低频率,对于有效实施紧急控制措施防止系统频率跌落有重要意义。
调速器对频率变化有至关重要的影响,调速器对发电机出力的控制决定了系统的有功功率平衡,其响应速度决定了系统最低频率的水平和出现时间,而调速器死区是客观存在的,不适当的死区设置都对系统有不良影响。现有的频率计算方法均未考虑调速器死区的影响,造成了与实际情况误差较大的情况,不利于系统频率安全。另一方面,调速器死区的存在为系统带来强烈的非线性问题,如何处理含死区的调速系统、定量考虑死区对最低频率的影响是亟待解决的问题。
综上,含死区的调速系统在实际电力系统中广泛存在,死区的大小对最低频率的影响较大,建立数学模型定量描述死区特性是提高最低频率计算准确性的有效方法。
发明内容
针对现有技术的缺点,本发明的目的是提供一种计及死区效应的电力 系统最低频率快速计算方法。
为了实现上述目的,本发明提供了一种计及死区效应的电力系统最低频率快速计算方法,其包括如下步骤:
A、获取相关模型参数
输入电力系统模型,包括所有机组的发电机、调速器与原动机模型以及模型的参数,扰动功率ΔpL;
B、简化调速器-原动机模型
将每台机组的调速器-原动机模型简化为传递函数为的一阶模型,其中Ni(bi,A,ki)为包含调速器单位调节功率ki、调速器死区bi、以及扰动后系统惯性中心最大频率偏差A的函数;时间常数Ti通过最小二乘法拟合在确定阶跃输入下的原始调速器-原动机模型的响应曲线确定;
C、建立用于频率偏差计算的系统等值模型
基于用于负荷频率控制分析的典型系统等值模型进行修改;将系统频率对调速器的反馈环节断开,用随时间线性变化的频率偏差Δω(t)作为新的调速器输入;每台机组的调速器-原动机模型用步骤B中的一阶简化模型替代;忽略负荷的频变效应,所有发电机的惯性等值为惯性中心的值;
D、计算扰动后系统最低频率
基于步骤C中建立的系统等值模型,根据步骤A中获取的参数和扰动功率ΔpL迭代求解系统惯性中心的最低频率。
快速准确地计算扰动后电力系统的最低频率,对于有效实施紧急控制措施防止发生频率崩溃有重要意义。最低频率发生在一次调频的过程中,调速器对于最低频率起了重要的决定作用,传统的最低频率预测方法都忽略了调速器的死区或没有计及死区的作用,这使得计算出的最低频率存在误差。本发明提出了一种计及调速器死区效应的电力系统扰动后最低频率快速计算方法,该方法能够定量考虑死区对最低频率的影响,提高了电力系统中切机故障后的最低频率计算的准确性。
本发明推导出的考虑调速器死区的扰动后最低频率计算方法,能够通过迭代计算快速得出结果,可以方便地用于在线应用。
本发明的方法将死区特性线性化,能定量描述含死区调速系统的输入输出关系,解决了其他计算方法无法考虑调速器死区的问题,减小了最低频率计算的误差。
本发明提供了一种定量描述调速器死区特性的频率偏差计算模型,基于该模型的最低频率计算方法能够处理死区带来的非线性问题,使本发明的方法能更有效地应用于实际系统中,减小因忽略死区而带来的计算误差。
根据本发明另一具体实施方式,步骤A具体包括如下步骤:
A1、获取发电机组i的参数与配置调速器-原动机模型
根据电力系统时域仿真动态数据,确定发电机组i的容量Si、惯性时间常数Mi,以及该机组的调速器-原动机传递函数模型,特别地,获取该调速器的死区大小bi,单位调节功率ki;
A2、获取扰动功率ΔpL:最低频率的产生通常是由于切机故障产生的,因此稳态时被切除发电机的出力PGi就是扰动功率大小;
上述获取模型和参数的方法通过电力系统时域仿真的动态模型数据即可获得,数据量较小,任何投运的实际电力系统均有相应的动态数据。
根据本发明另一具体实施方式,步骤B具体包括如下步骤:
B1、将发电机组i的调速器-原动机模型简化为一阶模型,传递函数为
其中Ni(bi,A,ki)为包含调速器单位调节功率ki、调速器死区bi、以及扰动后系统惯性中心最大频率偏差A的函数,前两个量都通过A步获取的参数得到;
B2、死区线性化
对于含死区的非线性调速器-原动机系统,输入与输出之间的关系可以 利用间隙特性进行描述,如图1所示,可以利用描述函数法进行线性化。
输出与输入的关系为Y=F(X),非线性函数F(X)可以用傅里叶级数展开,并仅取前两项得,
F(X)=N0+N·XL
式中
利用死区特性的典型函数描述函数可以得到上式满足以下关系:
本发明在处理死区特性时,采用描述函数法将死区线性化,可以做到定量计算死区对最低频率的影响,这是其他最低频率计算方法所不具有的。
B3、机组i的调速器-原动机系统的频率阶跃响应
根据A步中获得的模型传递函数和参数建立机组i的调速器-原动机的数学模型,给该模型值为0.01pu的频率阶跃输入,记录机械功率输出随时间变化的曲线;
B4、拟合响应曲线确定时间常数Ti
令A=0.01,根据式(1)计算一阶模型传递函数的分子项Ni(bi,A,ki),并设置0.01pu的频率阶跃输入,即
可得到一阶模型的输出
经过拉普拉斯反变换可得时域中含参数的输出
通过最小二乘法利用上式拟合B3步中获得的响应曲线,最终确定该机组的惯性时间常数Ti;
B5、利用步骤B2、步骤B3,获得所有机组的一阶简化模型惯性时间常数;
本发明按以上简化方法可以处理系统中存在不同调速系统的情况,并且能够定量考虑含死区的调速系统的传递函数,简化后的调速器-原动机模型不但能够方便地应用于多机系统计算中,还简化了计算提高了计算速度。
根据本发明另一具体实施方式,步骤C具体包括如下步骤:
C1、解耦系统对调速器的频率反馈
传统的用于负荷频率控制的系统等值模型如图2所示,
Meq是系统惯性中心的惯性时间常数
D为负荷的频变系统,断开系统对调速器频率反馈环节,并令调速器输入为
修改后的等值模型如图3所示
C2、简化调速器-原动机模型
将原调速器-原动机模型用B步中的一阶模型替代,并且忽略负荷的频变效应,系统等值模型变为图4所示。
此时调速器-原动机系统的输出为
进一步令ΔPGi(t)简化为随时间线性变化的函数,即
当到达最低频率时,上面两式应相等,即
tmin为到达最低频率的时间,由上式可得到Ci的表达式
由此得到用于频率偏差计算的等值模型如图5所示,m为切机后系统内的发电机数量;
该等值模型通过将系统频率对调速系统的反馈开环处理,将频率偏差输入线性化,简化了发电机转动方程;并且该模型不用考虑系统结构和网络方程,避免了复杂的矩阵运算。
根据本发明另一具体实施方式,步骤D具体包括如下步骤:
D1、迭代求解Ci
根据图5所示的模型,频率偏差的计算式为
求解式(7)可以得到
由于在最大频率偏差处因此可以得到
式(6)中包含最大频率偏差A,将式(9)代入式(8)可以得到Δω(tmin)的值,这就是最大频率偏差,即
通过上式推导,可以得到以下方程组
方程组(11)可以通过迭代求解,最终计算出Ci值;
D2、最低频率计算
由D1步求解出的Ci值,可以计算最低频率时的频率偏差值为
最终最低频率由fmin=f0(1-Δωmin)得到。
与现有技术相比,本发明具备如下有益效果:
本发明的方法将死区特性线性化,能定量描述含死区调速系统的输入输出关系,解决了其他计算方法无法考虑调速器死区的问题,减小了最低频率计算的误差。
下面结合附图对本发明作进一步的详细说明。
附图说明
图1是含死区的调速器-汽轮机间隙特性图;
图2是频率分析的系统等值模型图;
图3是开环的系统等值模型图;
图4是一阶模型的开环系统等值模型图;
图5是一阶模型的开环系统等值模型图;
图6是实施例1中,IEEE9节点的算例系统结构图,系统包括3台发电机,3个负荷节点,3个联络节点(不含负荷);
图7是实施例1中,调速器-原动机系统GS-TB模型的传递函数图;
图8是实施例1中,G1机组调速器-原动机模型的0.01pu频率阶跃响应曲线和一阶简化模型拟合曲线。
具体实施方式
实施例1
本实施例以IEEE9节点标准算例系统为算例演示具体实施方式,测试系统结构图如图6所示,其包括如下步骤:
A、获取相关模型参数
输入电力系统模型,包括所有机组的发电机、调速器与原动机模型以及模型的参数,扰动功率ΔpL;
获取相关参数的具体作法为:
A1、获取发电机组i的参数与配置调速器-原动机模型
根据电力系统时域仿真动态数据,确定发电机组i的容量Si、惯性时间常数Mi,以及该机组的调速器-原动机传递函数模型,特别地,获取该调速器的死区大小bi,单位调节功率ki;
A2、获取扰动功率ΔpL:最低频率的产生通常是由于切机故障产生的,因此稳态时被切除发电机的出力PGi就是扰动功率大小。
B、简化调速器-原动机模型
将每台机组的调速器-原动机模型简化为传递函数为 的一阶模型,其中Ni(bi,A,ki)为包含调速器单位调节功率ki、调速器死区bi、以及扰动后系统惯性中心最大频率偏差A的函数;时间常数Ti通过最小二乘法拟合在确定阶跃输入下的原始调速器-原动机模型的响应曲线确定。
上述简化机组i的调速器-原动机模型,具体作法如下:
B1、将发电机组i的调速器-原动机模型简化为一阶模型,传递函数为
其中Ni(bi,A,ki)为包含调速器单位调节功率ki、调速器死区bi、以及扰动后系统惯性中心最大频率偏差A的函数,前两个量都通过A步获取的参数得到;
B2、死区线性化
对于含死区的非线性调速器-原动机系统,输入与输出之间的关系可以利用间隙特性进行描述,如图1所示,可以利用描述函数法进行线性化。
输出与输入的关系为Y=F(X),非线性函数F(X)可以用傅里叶级数展开,并仅取前两项得,
F(X)=N0+N·XL
式中
利用死区特性的典型函数描述函数可以得到上式满足以下关系:
B3、机组i的调速器-原动机系统的频率阶跃响应
根据A步中获得的模型传递函数和参数建立机组i的调速器-原动机的数学模型,给该模型值为0.01pu的频率阶跃输入,记录机械功率输出随时 间变化的曲线;
B4、拟合响应曲线确定时间常数Ti
为了拟合B3步中的0.01pu的频率阶跃输入下的原调速器-原动机响应曲线,令A=0.01,根据式(1)计算一阶模型传递函数的分子项Ni(bi,A,ki),并设置0.01pu的频率阶跃输入,即
可得到一阶模型的输出
经过拉普拉斯反变换可得时域中含参数的输出
通过最小二乘法利用上式拟合B3步中获得的响应曲线,最终确定该机组的惯性时间常数Ti;
B5、利用B2、B3步中的步骤,获得所有机组的一阶简化模型惯性时间常数。
C、建立用于频率偏差计算的系统等值模型
基于用于负荷频率控制分析的典型系统等值模型进行修改;将系统频率对调速器的反馈环节断开,用随时间线性变化的频率偏差Δω(t)作为新的调速器输入;每台机组的调速器-原动机模型用B步中的一阶简化模型替代;忽略负荷的频变效应,所有发电机的惯性等值为惯性中心的值。
建立用于频率偏差计算的系统等值模型,具体作法如下:
C1、解耦系统对调速器的频率反馈
传统的用于负荷频率控制的系统等值模型如图2所示,
Meq是系统惯性中心的惯性时间常数
D为负荷的频变系统,断开系统对调速器频率反馈环节,并令调速器输入为
修改后的等值模型如图3所示。
C2、简化调速器-原动机模型
将原调速器-原动机模型用B步中的一阶模型替代,并且忽略负荷的频变效应,系统等值模型变为图4所示。
由上图的传递函数可得调速器-原动机系统的输出为
式(5)描述的调速器输出的变化率在初始时为0,稳态时达到 在系统达到最低频率过程中,调速器输出变化近似为线性变化,因此进一步令ΔPGi(t)简化为随时间线性变化的函数,即
Ci应是一个介于0与的值,当到达最低频率时,上面两式应相等,即
tmin为到达最低频率的时间,由上式可得到Ci的表达式
由此得到用于频率偏差计算的等值模型如图5所示,m为切机后系统 内的发电机数量;
D、计算扰动后系统最低频率
基于C步中建立的系统等值模型,根据A步中获取的参数和扰动功率ΔpL迭代求解系统惯性中心的最低频率。
计算扰动后系统最低频率,具体作法如下:
D1、迭代求解Ci
根据图5所示的模型,频率偏差的计算式为
求解式(7)可以得到
由于在最大频率偏差处因此可以得到
式(6)中包含最大频率偏差A,将式(9)代入式(8)可以得到Δω(tmin)的值,这就是最大频率偏差,即
通过上式推导,可以得到以下方程组
方程组(11)可以通过迭代求解,最终计算出Ci值;
D2、最低频率计算
由D1步求解出的Ci值,可以计算最低频率时的频率偏差值为
最终最低频率由fmin=f0(1-Δωmax)得到。
仿真实验:
采用IEEE9节点标准系统对本发明的方法进行测试,系统包含3台发电机,3个联络节点,3个负荷节点,系统结构图如图6所示。G1-G3发电机的调速器和原动机都采用汽轮机调速器模型GS和串联组合、单再热器汽轮机模型TB,模型传递函数图如图7所示,模型参数如表1和表2所示。
T<sub>1</sub> |
T<sub>2</sub> |
T<sub>3</sub> |
VEL<sub>open</sub> |
VEL<sub>close</sub> |
T<sub>CH</sub> |
0 |
0 |
0.5 |
2.0 |
2.0 |
0.2 |
F<sub>HP</sub> |
T<sub>RH</sub> |
F<sub>IP</sub> |
T<sub>CO</sub> |
F<sub>LP</sub> |
λ |
0.333 |
10.0 |
0.667 |
0 |
0 |
0 |
表1:G1-G3机组的调速器-原动机系统统一固定参数
ID |
M |
S |
k |
G1 |
19.1 |
247.5 |
20 |
G2 |
13.3 |
192 |
25 |
表2:G1-G3发电机惯性、额定容量、单位调节功率
利用MATLABR2012软件中的simulink工具箱建立GS-TB的模型,通过设置输入0.01,得到各机组的调速器-原动机响应曲线,并利用一阶简化模型去拟合该响应曲线最终得到各机组的时间常数T1~T3。即使是同一个调速器-原动机模型,对于不同的调速器死区,都需要进行0.01pu的频率阶跃响应实验以确定一阶简化模型的时间常数。图8的实线是G1机组的频率阶跃响应曲线,虚线为一阶简化模型的拟合曲线,结合表3可以看出死区大小对拟合的时间常数存在影响。
死区大小(pu) |
0.003 |
0.004 |
0.005 |
0.006 |
0.008 |
T<sub>1</sub>(s) |
6.3409 |
6.4381 |
6.5354 |
6.6283 |
6.8040 |
表3:G1一阶简化模型在不同死区的拟合时间常数
负荷采用恒功率负荷模型,扰动设置为在0时刻在节点BUSA上增加50MW有功功率,则扰动功率以系统发电机总容量为基准。为了检验算法的有效性,分别设置不同的调速器死区,比较时域仿真与计算结果的差异。
仿真结果与本实施例方法的对比结果如表4所示。前3个情境下的各机组死区相同并且逐渐增大,最低频率也逐渐降低,本文方法计算出的频率偏差与时域仿真结果误差最小能够达到0.003Hz,十分接近。后3个情境的各机组死区均不相同,计算误差保持在0.05Hz左右,证明本发明方法对机组间不同死区大小的情况同样适用。
表4:G1-G3机组在不同死区下的系统扰动后最低频率
虽然本发明以较佳实施例揭露如上,但并非用以限定本发明实施的范围。任何本领域的普通技术人员,在不脱离本发明的发明范围内,当可作些许的改进,即凡是依照本发明所做的同等改进,应为本发明的范围所涵盖。