法律状态公告日
法律状态信息
法律状态
2017-04-12
授权
授权
2015-04-22
实质审查的生效 IPC(主分类):G01R31/00 申请日:20141212
实质审查的生效
2015-03-25
公开
公开
技术领域
本发明涉及电力系统运行与控制领域,具体涉及电力系统阻尼比的评估,提出了一种 基于功率振荡次数评估系统阻尼比的方法。
背景技术
随着区域电网的互联,有功功率振荡逐渐成为影响大型互联电网安全稳定运行的主要 威胁之一。电网网架的日益庞大,快速励磁、电液伺服调速装置等的广泛应用,使得互联 电网出现有功功率振荡的风险大大增加,一方面抑制了电网的输电能力,另一方面也增加 了电网出现停电事故的概率。因此,大电网有功功率振荡问题是国内外电网发展必须面临 的一个难题。对于运行调度人员而言,根据广域测量系统(WAMS)的数据监控并预测电网 阻尼比的变化,在互联电网阻尼减弱时采取相应措施是有效避免有功功率振荡的方法之 一。
广域测量系统(WAMS)源自互联电网对于时间上同步和空间上广域的实际需求,它是 以相量测量单位(PMU)为基层单元采集信息,经过通信系统上传至调度中心,实现对互联 电网的监测,构成一个实时监控系3统。PMU利用全球定位系统(GPS)时钟同步的特点, 测量各节点以及线路在同一时间坐标下的各种状态量。WAMS信息目前主要被用于在线 分析互联电网振荡阻尼,扰动定位等方面,选取的是监测的实时数据,而其历史数据在电 网中应用很少。广域测量的历史数据中有电网每日发生的有效扰动次数,这里有效扰动定 义为:
当广域测量系统中有联络线有功功率振荡曲线同时满足:
a、振荡频率在0.1~2.5Hz;
b、联络线功率振荡连续r个周波的幅值大于给定阈值M;
时,判定该联络线发生了有效扰动;否则,判定该联络线没有发生有效扰动;M即为有效 扰动阈值。
电力系统分析综合程序(Power System Analysis Software Package)简称PSASP。它是 一套历史长久、功能强大、使用方便的商用电力系统分析程序,由中国电力科学研究院技 术分公司开发,具有我国自主知识产权,是资源共享,使用方便,高度集成和开放的大型 软件包。电力系统分析综合程序基于电网基础数据库、固定模型库以及用户自定义模型库 的支持,可进行电力系统(输电、供电和配电系统)的各种计算分析。能实现如下功能:潮 流计算,暂态稳定,短路电流,网损分析,电压稳定,静态安全分析,静态和动态等值, 直接法暂态稳定,小干扰稳定,最优潮流和无功优化,参数优化协调,继电保护整定与仿 真等。
互联电网小干扰线性化分析互联电网静态稳定的步骤(何仰赞,温增银.电力系统分 析(下册)[M].华中科技大学出版社,2002:226-229.)如下:
(1)列写电力系统各元件的微分方程以及联系各元件间关系的代数方程;
(2)分别对微分方程和代数方程线性化;
(3)消去方程中的非状态变量,求出线性化小扰动状态方程及矩阵A;
(4)进行给定运行情况的初态计算,确定A矩阵各元素的值;
(5)确定或判断A矩阵特征值实部的符号,判断系统在给定的运行条件下是否具有 静态稳定性。
杨东俊等的“基于WAMS量测数据的低频振荡机理分析”(电力系统自动化, 2009,33(23):24-28)针对华中电网WAMS)基于记录到的一次区域间低频振荡事件,采用 Prony方法分析振荡各阶段特征确定其振荡类型,并结合电网历史小干扰事件频率分布统 计情况验证结论。但是文中没有对小干扰事件频率分布统计情况进行分析解释。杨东俊等 的“基于WAMS的负阻尼低频振荡与强迫振荡的特征判别”(电力系统自动化,2013,37(13): 57-62)在WAMS的历史实测数据基础上,结合负阻尼振荡和强迫功率振荡事件发生前振 荡模式分布统计特征,对电力系统低频振荡性质进行了分析和判断。但是文中并未给出联 络线有效扰动次数与系统阻尼比对应关系的理论分析,难以定量确定有效扰动次数与系统 阻尼比的对应关系。
发明内容
本发明的目的在于提供一种基于功率振荡次数评估系统阻尼比的方法,该方法首先通 过广域测量系统获取所有节点的电压幅值和相角以及有效扰动次数的信息,通过电力系统 分析综合程序软件得到互联电网中各发电机的参数和所有线路的参数,将上述信息和参数 生成互联电网的小干扰线性化方程;然后求解小干扰线性化方程,得到各支路功率振荡的 时域表达式;结合各支路功率振荡的时域表达式和有效扰动的判定方法求得区域As和区 域At间在[f1,f2]的频率范围内联络线lh的有效扰动次数和互联电网总有效扰动次数的理论 比值与区域As和区域At间阻尼比的关系式;通过广域测量系统获取历史数据中区域As和区域At间在[f1,f2]的频率范围内各联络线每日有效扰动次数和互联电网总每日有效扰动 次数,计算得到区域As和区域At间各联络线每日有效扰动次数和互联电网总每日有效扰 动次数的实测比值;最后将理论比值与区域As和区域At间阻尼比的关系式绘制成曲线图, 根据实测比值查图得到对应的阻尼比。该方法具有机理明确、计算简单的特点,为运行调 度人员提供互联电网运行状态的参考。
为了达到上述目的,本发明采用如下技术方案:
一种基于功率振荡次数评估系统阻尼比的方法,该方法包括以下步骤:
(a)在有H个区域、L条区域间联络线、N条支路、n个节点、其中发电机节点m 个的互联电网中装设广域测量系统,该电网在所有支路上均装设同步相量测量单元;
(b)在广域测量系统获取所有节点的电压幅值和相角以及有效扰动次数的信息,通 过电力系统分析综合程序软件得到互联电网中各发电机的参数和所有线路的参数,将上述 信息和参数生成互联电网的小干扰线性化方程;
(c)求解步骤(b)的小干扰线性化方程,得到各支路功率振荡的时域表达式
其中:Δθi和Δθj分别表示第i节点和j个节点的电压相角变化量,ΔVi和ΔVj分别表示第i 节点和j个节点的电压幅值的变化量,θi0和θj0分别表示第i节点和j个节点的稳态电压相 角,Vi0和Vj0分别表示第i节点和j个节点的稳态电压幅值,Xij表示第i节点和j个节点之 间的支路电抗值,上述i=1,2,…,n;j=1,2,…,n;
(d)在区域As和区域At间有g条联络线,这两个区域间的振荡频率在[f1,f2]之间,s、 t=1,2,…,H,s≠t,结合步骤(c)中各支路功率振荡的时域表达式和有效扰动的判定方法 求得在[f1,f2]的频率范围内联络线lh的有效扰动次数qp和互联电网总有效扰动次数Q的理 论比值与区域As和区域At间阻尼比的关系式
其中:p=1,2,…,L,M为有效扰动阈值,Kl(p)(ξ)为第 p条联络线功率振荡幅值在t=0时刻的初始值,μ1、σ1分别为冲击扰动面积F服从的正态 分布的均值和方差,h=1,2,…,g;
(e)通过广域测量系统获取历史数据中区域As和区域At间在[f1,f2]的频率范围内各 联络线每日有效扰动次数qaw和互联电网总每日有效扰动次数Q',计算得到区域As和区域 At间各联络线每日有效扰动次数和互联电网总每日有效扰动次数的实测比值 QAw=qaw/Q',w=1,2,…,g;
(f)将步骤(d)的关系式绘制成区域As和区域At间各联络线对应的曲线图,共g 个,并根据步骤(e)求得的g个联络线每日有效扰动次数qaw和互联电网总每日有效扰动 次数Q'的实测比值QAw查图得到g个比值对应的阻尼比ξp,p=1,2,…,g;取作 为区域As和区域At间阻尼比的评估值。
本发明包括以下优点和技术效果:
1、新颖性:电网中广域测量系统的历史数据使用较少,广域测量系统的历史数据 中统计有电网每日发生的有效扰动次数,该方法确立了历史数据中有效扰动次 数与系统阻尼之间存在的关系式,从而实现系统阻尼比的评估,为低频振荡的 在线监测和预警提供经验判据;
2、计算简单:理论计算方法容易且准确性强,避免了大量重复的蒙特卡洛抽样;
3、实用性强:提供了在[f1,f2]的频率范围内两区域间联络线的有效扰动次数和互联 电网总有效扰动次数的理论比值与区域间阻尼比的关系式。并通过广域测量系 统获取历史数据中区域间在[f1,f2]的频率范围内各联络线每日有效扰动次数和 互联电网总每日有效扰动次数的实测比值对系统阻尼比进行评估,为实际运行 人员预估系统阻尼比提供参考依据。
附图说明
图1是一种基于功率振荡次数评估系统阻尼比的方法的流程图。
图2为4机11节点互联电网结构示意图。
图3为4机11节点互联电网联络线7-8第一回线上有效扰动次数q1和互联电网总有 效扰动次数Q的理论比值Q1与区域一和区域二间阻尼比ξ的关系曲线。
图4为4机11节点互联电网联络线7-8第二回线上有效扰动次数q2和互联电网总有 效扰动次数Q的理论比值Q2与区域一和区域二间阻尼比ξ的关系曲线。
具体实施方式
以下结合附图对本发明进行进一步的说明。
本发明提供了一种基于功率振荡次数评估系统阻尼比的方法,该方法首先通过广域测 量系统获取所有节点的电压幅值和相角以及有效扰动次数的信息,通过电力系统分析综合 程序软件得到互联电网中各发电机的参数和所有线路的参数,将上述信息和参数生成互联 电网的小干扰线性化方程;然后求解小干扰线性化方程,得到各支路功率振荡的时域表达 式;结合各支路功率振荡的时域表达式和有效扰动的判定方法求得区域As和区域At间在 [f1,f2]的频率范围内联络线lh的有效扰动次数和互联电网总有效扰动次数的理论比值与区 域As和区域At间阻尼比的关系式;通过广域测量系统获取历史数据中区域As和区域At间在[f1,f2]的频率范围内各联络线每日有效扰动次数和互联电网总每日有效扰动次数,计 算得到区域As和区域At间各联络线每日有效扰动次数和互联电网总每日有效扰动次数的 实测比值;最后将理论比值与区域As和区域At间阻尼比的关系式绘制成曲线图,根据实 测比值查图得到对应的阻尼比。
实施例1
由图1可知,一种基于功率振荡次数评估系统阻尼比的方法,该方法包括以下步骤:
(a)在有H个区域、L条区域间联络线、N条支路、n个节点、其中发电机节点m 个的互联电网中装设广域测量系统,该电网在所有支路上均装设同步相量测量单元;
(b)在广域测量系统获取所有节点的电压幅值和相角以及有效扰动次数的信息,通 过电力系统分析综合程序软件得到互联电网中各发电机的参数和所有线路的参数,将上述 信息和参数生成互联电网的小干扰线性化方程;
首先说明的是,在本发明中,将互联电网中的发电机母线、负荷母线、联络母线称为 节点。通过广域测量系统获取互联电网节点的电压幅值Vi和相角θi、各联络线每日有效扰 动次数qaw和互联电网总每日有效扰动次数Q';i=1,2,…,n;w=1,2,…,g。n为互联电网总节 点数;g为区域As和区域At间的联络线条数,s、t=1,2,…,H,s≠t;H为互联电网区域数。 通过电力系统分析综合程序软件获取的参数包括:同步发电机本体的模型结构及对应的电 气参数,例如对同步发电机四阶模型,需要获取的电气参数有发电机额定容量SGN,额定 有功功率输出PGN,额定功率因数发电机惯性时间常数TJ,发电机定子电阻Ra、 定子漏抗Xl,发电机d轴同步电抗Xd、暂态电抗X'd、开路瞬变时间常数T'd0、q轴同步 电抗Xq、暂态电抗X'q、开路瞬变时间常数T'q0;励磁系统模型结构及对应的结构参数; 调速系统模型结构及对应的结构参数;电力系统稳定器(PSS)模型结构及对应的结构参数; 线路电气参数,包括线路电阻Rline、线路电抗Xline、线路电容Bline。
使用潮流方程作为网络方程,互联电网线性化方程(Sauer PW and Pai MA,“Power system dynamics and stability,”Upper Saddle River,NJ:Prentice Hall,1998:222-230)表示 为:
C2ΔX+D3ΔIg+D4ΔVg+D5ΔVl=0 \*MERGEFORMAT(1)
D6ΔVg+D7ΔVl=0
其中,四个式子分别为线性化后的发电机微分方程、发电机定子电压方程、发电机节点的 潮流方程和非发电机节点的潮流方程。式中,ΔX表示互联电网的状态变量变化量;表 示互联电网的状态变量变化量对时间的导数;ΔIg=[ΔId1ΔIq1…ΔIdiΔIqi…ΔIdmΔIqm]T,ΔIdi、 ΔIqi分别为第i个发电机节点的d轴、q轴电流分量的变化量,i=1,2,…,m;ΔVg=[Δθ1ΔV1…ΔθiΔVi…ΔθmΔVm]T,Δθi、ΔVi分别为第i个发电机节点的电压相角和幅值的变化量, i=1,2,…,m;ΔVl=[Δθm+1ΔVm+1…ΔθiΔVi…ΔθnΔVn]T,Δθi、ΔVi分别为第i个非发电机节点 的电压相角和幅值的变化量,i=m+1,2,…,n;m为发电机节点数,n为互联电网总节点数。 A1到D7的表达式与发电机、励磁系统、调速系统、电力系统稳定器(PSS)选取的模型结构 均密切相关,以不考虑励磁系统、调速系统和电力系统稳定器(PSS),发电机采用四阶模 型为例,各矩阵表示为:
式中,
式中,
通过矩阵运算,将式(1)变换为:
式中,A为互联电网状态空间矩阵;ΔX表示互联电网的状态变量变化量;表示互联 电网的状态变量变化量对时间的导数。
互联电网出现功率扰动表达式为ΔP1。以负荷节点i(i=m+1,2,…,n)为例,将ΔP1添加 在线性化后的网络方程上,得到节点i的网络方程为
ΔP1+[D6(i,1)...D6(i,j)...D6(i,m)]ΔVg+[D7(i,m+1)...D7(i,j)...D7(i,n)]ΔVl=0 \*MERGEFORMAT(15)
式中,
D6(i,j)=-Vi0[Vj0Yijsin(θi0-θj0-αij)Yijcos(θi0-θj0-αij)](j=1,2,...,m) \*MERGEFORMAT(16)
D7(i,j)=-Vi0[Vj0Yijsin(θi0-θj0-αij)Yijcos(θi0-θj0-αij)] \*MERGEFORMAT(17)
(j=m+1,...,n且j≠i)
(j=m+1,...,n且j=i)
式中,Vi0、Vk0、θi0、θk0分别为i、k节点的线性化稳态运行点的电压幅值和相角,Yik、 αik分别为i、k节点间导纳的幅值和相角。
将式(1)中的非发电机网络方程改写为
C4ΔP1+D6ΔVg+D7ΔVn=0 \*MERGEFORMAT(19)
式中,C4=[0 … 1 0 … 0]T,为1处是第i号节点的有功功率线性化方程,i=m+1,2,…,n;
D6=[D6(i,1)...D6(i,j)...D6(i,m)];D7=[D7(i,m+1)...D7(i,j)...D7(i,n)]。
对改写后的式(1)进行矩阵变换,原小干扰线性化方程转化为:
式中,A为原小干扰线性化方程状态空间矩阵。KF为分配因子,表达式写为:
式中,
同理,如果功率扰动发生在发电机节点i(i=1,2,…,m)上,式(1)中的发电机节点潮流 方程变为:
C3ΔP1+C2ΔX+D3ΔIg+D4ΔVg+D5ΔVn=0 \*MERGEFORMAT(22)
式中,C3=[0…10…0]T,为1处是该发电节点的节点潮流有功线性化方程;互联电网 动态方程仍用式(8)表示,KF表达式为:
至此,互联电网小干扰线性化方程表达式及所有参数参数含义已全部给出。
(c)求解步骤(b)的小干扰线性化方程,得到各支路功率振荡的时域表达式
对步骤(b)计算得到的状态空间矩阵A进行特征值计算,能够得到特征根对角阵 Λ=diag(λ1,λ2,...,λN),右特征向量矩阵U=[u1,u2,...uN],左特征向量矩阵的转置VT=U-1。 根据左右特征向量的定义有:
U-1AU=Λ \*MERGEFORMAT(24)
令ΔX=UΔZ,则在互联电网某一节点加入功率扰动后,有
式中,ΔZ是ΔX经过线性变换后解耦的状态变量的变化量,是ΔZ对时间的导数。 互联电网中观测到的主要是阻尼较弱的区域间振荡模式,设这一主振荡模式对应的特征值 为λi=-α+jωd,λj=-α-jωd(i、j表示该振荡模式的特征根在对角阵A中的行数,α为特征根 实部的相反数,ωd为特征根虚部),振荡频率为fd=ωd/(2π),功率扰动是面积为F的冲击扰 动Fδ(t),则主振荡模式对应的方程为:
式中,ΔZi(t)、ΔZj(t)分别是第i、j个解耦状态变量的变化量,ViT、VjT分别表示VT的第i、 j行。
再由ΔX=UΔZ,得到
式中,Δx(t)是互联电网状态变量ΔX的时域表达形式,(-α)为特征根实部,ωd为特征根虚 部,将式(28)代入式(1)中,求解出ΔIg、ΔVg、ΔVl,计算式如下:
式中,
联络线i-j上的功率振荡ΔPij用节点i、j的电压幅值和相角表示为
式中,Δθi、Δθj、ΔVi、ΔVj分别为第i、j个节点的电压相角和幅值的变化量;θi0、θj0、 Vi0、Vj0分别为第i、j个节点的稳态电压相角和幅值;Xij为第i、j个节点间的联络线电 抗值;i=1,2,…,n;j=1,2,…,n。在ΔV中提取i、j节点对应的电压相角和幅值的变化量, 代入式(32)中得到ΔPij的表达式为
ΔPij=TΔV=T'ΔX \*MERGEFORMAT(33)
式中,ΔV=[ΔVgΔVl]T,是2n×1的矩阵;T是1×2n的矩阵,
联络线i-j上的功率振荡ΔPij的时域表达式写为:
其中,Kl=T'KX,其余参数含义同式(28)。
至此,联络线功率振荡的时域表达式已全部给出。
(d)在区域As和区域At间有联络线g条,这两个区域间的振荡频率在[f1,f2]之间,s、 t=1,2,…,H,s≠t,结合步骤(c)中各支路功率振荡的时域表达式和有效扰动的判定方法 求得在[f1,f2]的频率范围内联络线lp的有效扰动次数qp和互联电网总有效扰动次数Q的理 论比值Qp与区域As和区域At间阻尼比ξ的关系式
其中:p=1,2,…,g,M为有效扰动阈值,Kl(p)(ξ)为 第p条联络线功率振荡幅值在t=0时刻的初始值,μ1、σ1分别为冲击扰动面积F服从的正 态分布的均值和方差;
当广域测量系统中有联络线有功功率振荡曲线同时满足:a、振荡频率在0.1~2.5Hz, b、联络线功率振荡连续r个周波,幅值大于M时,判定该联络线发生了一次有效扰动; 否则,判定该联络线没有发生有效扰动。
以扰动类型为冲击扰动为例,其表达式为Fδ(t),冲击面积F服从正态分布这里,μ1为正态分布的均值,σ1为正态分布的方差。取r=5,则联络线有功功率振荡被 计作一次有效扰动的条件为:
f(ξ,t)|t=5T=e-α(ξ)tKl(ξ)F|t=5T=h(ξ)>M \*MERGEFORMAT(35)
式中,ξ是主振荡模式对应的阻尼比,M为联络线功率振荡阈值,Kl定义同式(34),其余参数含义同式(28)。
那么,联络线有功功率振荡被记录成一次有效扰动的概率为
式中,
式中各参数含义同式(35)。
在区域As和区域At间有联络线g条,这两个区域间的振荡频率fd在[f1,f2]之间,s、 t=1,2,…,H,s≠t,联络线lp的有效扰动次数qp和互联电网总有效扰动次数Q的理论比值Qp与区域As和区域At间阻尼比ξ的关系式为
式中,p=1,2,…,g,M为有效扰动阈值,Kl(p)(ξ)为第p条联络线 功率振荡幅值在t=0时刻的初始值,μ1、σ1分别为冲击扰动面积F服从的正态分布的均值 和方差;
(e)通过广域测量系统获取历史数据中区域As和区域At间在[f1,f2]的频率范围内各 联络线每日有效扰动次数qaw和互联电网总每日有效扰动次数Q',计算得到区域As和区域 At间各联络线每日有效扰动次数和互联电网总每日有效扰动次数的实测比值 QAw=qaw/Q',w=1,2,…,g;
(f)将步骤(d)的关系式绘制成区域As和区域At间各联络线对应的曲线图,共g 个,并根据步骤(e)求得的g个联络线每日有效扰动次数qaw和互联电网总每日有效扰动 次数Q'的实测比值QAw查图得到g个比值对应的阻尼比ξp,p=1,2,…,g;取作 为区域As和区域At间阻尼比的评估值。
实施例2
在4机11节点标准测试算例上进行本发明的测试,其网络结构图如图2所示。研究中 针对互联电网中某一阻尼较弱模式,即主要关注4机11节点互联电网中阻尼最弱的区域 间振荡模式,调整PSS使其余振荡模态阻尼比较高,得到该互联电网的振荡模式如表1 所示。
表1 4机11节点互联电网振荡模式
外加冲击扰动面积F服从N(μ,σ2),取μ=0.5,σ2=0.08,施加冲击扰动的节点位置分 布服从平均分布,计算得到在4机11节点互联电网上在[0.50Hz,0.55Hz]的频率范围内区 域一和区域二间2条联络线上有效扰动次数和互联电网总有效扰动次数的理论比值Qp与 区域一和区域二间阻尼比ξ的关系曲线,如图3-图4所示。在仿真模型系统上施加1000次 冲击扰动,记录到在2条联络线上的有效扰动次数分别为q1=795,q2=795;互联电网总有 效扰动次数为Q'=3878,实测比值分别为QA1=0.205,QA2=0.205;通过查询对应曲线得 到阻尼比分别为ξ1=0.04,ξ2=0.04,故与表1中模式1对应的阻尼比0.039 基本相同。由此,实现了对系统阻尼比的评估。
机译: 一种利用优化技术计算功率振荡阻尼比的方法
机译: 一种利用优化技术计算功率振荡阻尼比的方法
机译: 储罐阻尼比估计方法,以及基于储罐阻尼比估计方法的减振罐设计方法