CN103036498B - 一种基于pmu的同步发电机实用模型参数校核与辨识方法 - Google Patents
一种基于pmu的同步发电机实用模型参数校核与辨识方法 Download PDFInfo
- Publication number
- CN103036498B CN103036498B CN201210397510.8A CN201210397510A CN103036498B CN 103036498 B CN103036498 B CN 103036498B CN 201210397510 A CN201210397510 A CN 201210397510A CN 103036498 B CN103036498 B CN 103036498B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- mtr
- msubsup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 232
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000005259 measurement Methods 0.000 claims abstract description 15
- 238000004804 winding Methods 0.000 claims description 111
- 230000005284 excitation Effects 0.000 claims description 44
- 238000013016 damping Methods 0.000 claims description 32
- 230000001052 transient effect Effects 0.000 claims description 22
- 230000004907 flux Effects 0.000 claims description 21
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 8
- 230000002068 genetic effect Effects 0.000 claims description 4
- 238000004088 simulation Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 6
- 238000013461 design Methods 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 5
- 230000004044 response Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 238000001914 filtration Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
Landscapes
- Control Of Eletrric Generators (AREA)
Abstract
本发明提供一种基于PMU的同步发电机实用模型参数校核与辨识方法,以PMU实测数据为基础,以系统扰动后发电机高压母线电压幅值与相位、励磁电压实测值作为输入进行发电机实用参数的校核,判断实用参数的有效性;利用PMU可以实测发电机功角的功能,计算电压、电流的d、q轴分量,进而用最小二乘法分别辨识两个轴的PARK模型参数,最后求出发电机实用模型参数。本发明提高电力系统稳定分析的可信度,为电网的规划设计和调度决策提供可靠的参考。
Description
技术领域
本发明属于同步发电机参数辨识技术领域,具体涉及一种基于PMU的同步发电机实用模型参数校核与辨识方法。
背景技术
准确的发电机模型参数是电力系统稳定分析研究的关键,直接影响了电力系统的安全运行。近年来,随着系统辨识理论、相角测量技术(PMU)、全球定位技术(CPS)等新理论和技术在电力系统的广泛应用,建立在这些新技术平台上的发电机参数辨识取得了很大的发展,且各辨识算法呈现出相互交叉融合的特点。
在各种发电机参数辨识中,应用最广主要有时域辨识法、频域辨识法、抛载法等。时域辨识法基于励磁阶跃或负荷小扰动的动态响应来进行辨识,属于在线辨识方法,能够自然计及工况的影响,一旦辨识成功,那些饱和、涡流和旋转等因素就自然包含在参数估计值中,不用附加过多的假设条件。现有算法主要是最小二乘法和卡尔曼滤波法。随着数学技术的进步,时域辨识法得到了长足的发展,目前进化策略法,基因法、小波分析法、神经网络法以及粒子群优化法等均在发电机参数辨识中得到了应用,这些辨识算法在基于各种仿真数据辨识时往往可以得到比较满意的结果。然而,在针对实测数据时,常常会遇到不同试验甚至重复的试验辨识所得参数相差很大的情况。这是由于在线试验扰动较小,次暂态过程的可观测性较差,加上环境噪声和功角测量精度的影响,导致上述算法很难获得准确的辨识结果。
频域分析法可分为直流衰减法、静态频域法(SSFR)和动态频域法(OLFR)。直流衰减法由于响应持续过程短,影响频率特性的量测精度,目前较少使用。静态频域法需要大功率的变频电源作为信号源,由于需要不同频率的信号输入,试验时间很长且所得结果无法反映饱和效应。不过近年来随着多正弦输出信号电源的出现,试验耗时的问题得以解决。动态频域法得到的参数更加能够反映实际运行工况,更适用于动态稳定的研究。静态频域法和动态频域法的结合是近年来的发展趋势,即由SSFR给出参数的初值再根据一定运行条件下的OLFR法来修正参数的方法可改善算法的稳定性,且具有一定的滤波能力,不过对输人扰动信号的波形、幅值大小及其相关性要求严格。此外,频域响应分析建立在线性系统的基础上,不能反映同步电机参数的非线性特点。
抛载法试验于20世纪70年代提出后即获得了广泛的应用,其优点是易于实现,并计及工况对参数的影响,计算方法相对简单,但需要安排特定的试验。
发明内容
为了克服上述现有技术的不足,本发明提供一种基于PMU的同步发电机实用模型参数校核与辨识方法,以便提高电力系统稳定分析的可信度,为电网的规划设计和调度决策提供可靠的参考。本发明以PMU实测数据为基础,以系统扰动后发电机高压母线电压幅值与相位、励磁电压实测值作为输入进行发电机实用参数的校核,判断实用参数的有效性;利用PMU可以实测发电机功角的功能,计算电压、电流的d、q轴分量,进而用最小二乘法或遗传算法分别辨识两个轴的PARK模型参数,最后求出发电机实用模型参数。
本发明提供的基于PMU实测数据的同步发电机实用参数校核与辨识方法有两个特点:一是以PMU记录的发电机高压母线电压幅值与相角、励磁电压为输入量,仿真发电机的有功、无功输出,通过比较与实测有功、无功的误差判断发电机实用模型参数的有效性;二是通过PMU记录的扰动过程中的数据,根据发电机功角将电压、电流进行d、q分解,以电流和励磁电压为输入、机端电压为输出,采用最小二乘法分别辨识d、q轴发电机PARK模型参数及转子的转动惯量,最后求出发电机实用模型参数。
为了实现上述发明目的,本发明采取如下技术方案:
提供一种基于PMU的同步发电机实用模型参数校核与辨识方法,所述方法包括以下步骤:
步骤1:校核同步发电机实用模型参数;
步骤2:辨识同步发电机实用模型参数。
所述同步发电机实用模型参数包括同步发电机d轴暂态同步电抗x′d、同步发电机d轴次暂态同步电抗x″d、同步发电机q轴暂态同步电抗x′q、同步发电机q轴次暂态同步电抗x″q、同步发电机d轴暂态开路时间常数T′d0、同步发电机d轴次暂态开路时间常数T″d0、同步发电机q轴暂态开路时间常数T′q0和同步发电机q轴次暂态开路时间常数T″q0。
所述步骤1包括以下步骤:
步骤1-1:提取PMU记录的动态数据;所述动态数据包括有功功率P、无功功率Q、高压母线电压幅值U、高压母线电压相角和励磁电压Uf;
步骤1-2:将同步发电机组与系统解耦,形成等值系统;所述等值系统用同步发电机组微分方程、机端母线与高压母线的代数方程联立表示;
在考虑同步发电机阻尼绕组但不计饱和影响的情况下,同步发电机电压方程为
其中:
kd=(xd-x′d)/(x′d-x″d) (3)
kq=(xq-x′q)/(x′q-x″q) (4)
其中,E′q为同步发电机q轴暂态电势,E"q为同步发电机q轴次暂态电势,E′d为同步发电机d轴暂态电势,E"d为同步发电机d轴次暂态电势,id为同步发电机d轴定子电流,iq为同步发电机q轴定子电流,xd为同步发电机d轴同步电抗,xq为同步发电机q轴同步电抗,Uf为同步发电机励磁电压实测有名值,Rf为励磁绕组电阻有名值,If0为初始稳态时励磁电流有名值,
运用隐式梯形积分公式,则有:
其中:
其中,Δt为仿真步长;
忽略转速变化和定子回路电磁暂态过程影响,在t+Δt时刻的定子电压平衡方程为:
其中:ra为同步发电机定子电阻,且
将式(9)变换到XY坐标可得:
其中:
ud(t)为t时刻同步发电机d轴电压分量,uq(t)为t时刻同步发电机q轴电压分量,id(t)为t时刻同步发电机d轴定子电流,iq(t)为t时刻同步发电机q轴定子电流,δ(t)为同步发电机在t时刻的功角,为t时刻同步发电机定子电压相量,为t时刻同步发电机定子电流相量;
同步发电机转子运动方程为:
其中,Pm(t)为同步发电机t时刻机械功率,通常取0;Pe(t)为同步发电机t时刻电磁功率,Pe(t)=ud(t)*id(t)+uq(t)*iq(t),且
其中,ωs为同步发电机的同步机械角速度,Tj为同步发电机转子的转动惯量,ω(t)为同步发电机t时刻电角速度标幺值;
同步发电机端电压与高压母线电压代数方程为:
其中,RT+jXT为同步发电机出口变压器的阻抗,为高压母线电压相量;
联立(11)、(13)和(15),将带入实测值作为已知量计算出id(t+Δt)、ud(t+Δt)、uq(t+Δt)和iq(t+Δt),进而计算同步发电机输出功率:
有功功率P、无功功率Q与实测有功功率Pm实测无功功率Qm之间误差反映同步发电机参数描述同步发电机外特性的准确性,进而计算同布发电机参数误差Err,其表示为:
其中i为记录数据的长度,Δt为采样间隔;
当误差Err不满足要求时,同步发电机参数则重新进行辨识。
所述步骤2中,采用标幺值表示同步发电机PARK电压方程:
其中,ud为同步发电机d轴电压分量,uq为同步发电机q轴电压分量,uf为同步发电机励磁电压,rf为d轴励磁绕组电阻,rD为d轴阻尼绕组电阻,rg为q轴等值g绕组电阻,rQ为q轴等值Q绕组电阻,id为d轴阻尼绕组电流,if为励磁绕组电流,iD为d轴阻尼绕组电流,ig为q轴等值g绕组电流,iQ为q轴等值Q绕组电流,p为微分算子,ψd为d轴定子绕组磁链,ψf为发电机励磁绕组磁链,ψD为d轴阻尼绕组磁链,ψq为q轴定子绕组磁链,ψg为q轴等值g绕组磁链,ψQ为发电机q轴等值Q绕组磁链,ω为同步发电机角速度;
同步发电机PARK磁链方程表示为:
xd为同步发电机d轴同步电抗,xq为同步发电机q轴同步电抗,xf为发电机励磁绕组自感抗,xD为d轴阻尼绕组自感抗,xg为q轴等值g绕组自感抗,xQ为q轴等值Q绕组自感抗,xaf为发电机d轴定子绕组与励磁绕组间互感抗,xaD为发电机d轴定子绕组与阻尼绕组间互感抗,xag为发电机q轴定子绕组与q轴等值g绕组间互感抗,xaQ为发电机q轴定子绕组与q轴等值Q绕组互感抗,xfD为发电机励磁绕组与d轴阻尼绕组间互感抗,xgQ为q轴等值g绕组与等值Q绕组间的互感抗;
由于采用Xad基值系统,则有
xaf=xaD=xad(20)。
xag=xaQ=xaq
同步发电机实用模型参数分为稳态和发生扰动两种情况进行辨识;
A)在稳态情况下,同步发电机的阻尼绕组电流为0,同步发电机的初始稳态与扰动后稳态的机端电压和电流有名值方程如下:
其中,Uq0为初始稳态时q轴励磁电压有名值,Uq∞为扰动后q轴励磁电压有名值,Ud0为初始稳态时d轴励磁电压有名值,Ud∞为扰动后d轴励磁电压有名值,Ra为定子电阻有名值,Iq0为初始稳态时q轴定子电流有名值,Iq∞为扰动后q轴定子电流有名值,Id0为初始稳态时d轴定子电流有名值,Id0为扰动后d轴定子电流有名值,Xd为d轴同步电抗有名值、Xq为q轴同步电抗有名值,Xaf为d轴定子与励磁绕组互感抗有名值,If0为初始稳态时励磁电流有名值,If∞为扰动后励磁电流有名值;
通过方程(21)求出Ra、Xd、Xq和Xaf,再利用初始稳态时励磁电流与励磁电压的关系,即可计算Rf:
其中,Uf0为初始稳态时励磁电压有名值;
B)扰动后,假设励磁系统电压基值为K,作为待辨识的量,则有:
其中:Sn为同步发电机额定容量,UB为同步发电机相电压峰值,均为已知量。
发生扰动后,辨识过程包括d轴参数辨识、q轴参数辨识和同步发电机转子运动方程参数辨识,具体如下:
a)轴参数辨识
根据磁链方程推出:
进而有
在Uf实测并采用有名值的情况下:
因此建立状态方程:
量测方程为:
b)q轴参数辨识
假设同步发电机q轴有两个阻尼绕组,则阻尼绕组电流为:
因此建立状态方程:
量测方程为:
c)同步发电机转子运动方程参数辨识
其中,D为风阻系数,通常取0;Δδ(t)为同步发电机在t时刻的功角变化量;
选择Δδ、Δω作为状态量,选择Δδ(t)作为观测量,利用最小二乘法或遗传算法进行参数拟合,求出辨识参数Tj和D。
利用xf、xD、K、rD、xQ、xg、xag、rQ和rg根据所述辨识参数计算同步发电机实用模型参数:
与现有技术相比,本发明的有益效果在于:
1.本发明校核与辨识同步发电机实用参数时,不需要进行特别的试验,只是基于系统的扰动,因此简单有效,便于实现;
2.辨识参数所用的数据是发电机在线运行数据,所以辨识结果包含了饱和影响,不需要在进行额外修正;
3.本发明利用同步发电机端PMU的实测数据可以实现发电机与系统解耦,准确判断发电机参数的有效性。该算法可以直接植入PMU中,方便校核发电机参数;
4.直接发电机PARK模型参数,避免了辨识实用模型参数时多解的问题,同时也辨识出使用Xad基值系统时发电机励磁系统电压基值。
附图说明
图1是本发明中同步发电机实用参数校核原理图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
提供一种基于PMU的同步发电机实用模型参数校核与辨识方法,所述方法包括以下步骤:
步骤1:校核同步发电机实用模型参数;
步骤2:辨识同步发电机实用模型参数。
所述同步发电机实用模型参数包括同步发电机d轴暂态同步电抗x′d、同步发电机d轴次暂态同步电抗x″d、同步发电机q轴暂态同步电抗x′q、同步发电机q轴次暂态同步电抗x″q、同步发电机d轴暂态开路时间常数T′d0、同步发电机d轴次暂态开路时间常数T″d0、同步发电机q轴暂态开路时间常数T′q0和同步发电机q轴次暂态开路时间常数T″q0。
所述步骤1包括以下步骤:
步骤1-1:提取PMU记录的动态数据;所述动态数据包括有功功率P、无功功率Q、高压母线电压幅值U、高压母线电压相角和励磁电压Uf;
步骤1-2:将同步发电机组与系统解耦,形成等值系统;所述等值系统用同步发电机组微分方程、机端母线与高压母线的代数方程联立表示;
在考虑同步发电机阻尼绕组但不计饱和影响的情况下,同步发电机电压方程为
其中:
kd=(xd-x′d)/(x′d-x″d) (3)
kq=(xq-x′q)/(x′q-x″q) (4)
其中,E′q为同步发电机q轴暂态电势,E"q为同步发电机q轴次暂态电势,E′d为同步发电机d轴暂态电势,E"d为同步发电机d轴次暂态电势,id为同步发电机d轴定子电流,iq为同步发电机q轴定子电流,xd为同步发电机d轴同步电抗,xq为同步发电机q轴同步电抗,Uf为同步发电机励磁电压实测有名值,Rf为励磁绕组电阻有名值,If0为初始稳态时励磁电流有名值,
运用隐式梯形积分公式,则有:
其中:
其中,Δt为仿真步长;
忽略转速变化和定子回路电磁暂态过程影响,在t+Δt时刻的定子电压平衡方程为:
其中:ra为同步发电机定子电阻,且
将式(9)变换到XY坐标可得:
其中:
ud(t)为t时刻同步发电机d轴电压分量,uq(t)为t时刻同步发电机q轴电压分量,id(t)为t时刻同步发电机d轴定子电流,iq(t)为t时刻同步发电机q轴定子电流,δ(t)为同步发电机在t时刻的功角,为t时刻同步发电机定子电压相量,为t时刻同步发电机定子电流相量;
同步发电机转子运动方程为:
其中,Pm(t)为同步发电机t时刻机械功率,通常取0;Pe(t)为同步发电机t时刻电磁功率,Pe(t)=ud(t)*id(t)+uq(t)*iq(t),且
其中,ωs为同步发电机的同步机械角速度,Tj为同步发电机转子的转动惯量,ω(t)为同步发电机t时刻电角速度标幺值;
同步发电机端电压与高压母线电压代数方程为:
其中,RT+jXT为同步发电机出口变压器的阻抗,为高压母线电压相量;
联立(11)、(13)和(15),将带入实测值作为已知量计算出id(t+Δt)、ud(t+Δt)、uq(t+Δt)和iq(t+Δt),进而计算同步发电机输出功率:
有功功率P、无功功率Q与实测有功功率Pm实测无功功率Qm之间误差反映同步发电机参数描述同步发电机外特性的准确性,进而计算同布发电机参数误差Err,其表示为:
其中i为记录数据的长度,Δt为采样间隔;
当误差Err不满足要求时,同步发电机参数则重新进行辨识。
所述步骤2中,采用标幺值表示同步发电机PARK电压方程:
其中,ud为同步发电机d轴电压分量,uq为同步发电机q轴电压分量,uf为同步发电机励磁电压,rf为d轴励磁绕组电阻,rD为d轴阻尼绕组电阻,rg为q轴等值g绕组电阻,rQ为q轴等值Q绕组电阻,id为d轴阻尼绕组电流,if为励磁绕组电流,iD为d轴阻尼绕组电流,ig为q轴等值g绕组电流,iQ为q轴等值Q绕组电流,p为微分算子,ψd为d轴定子绕组磁链,ψf为发电机励磁绕组磁链,ψD为d轴阻尼绕组磁链,ψq为q轴定子绕组磁链,ψg为q轴等值g绕组磁链,ψQ为发电机q轴等值Q绕组磁链,ω为同步发电机角速度;
同步发电机PARK磁链方程表示为:
xd为同步发电机d轴同步电抗,xq为同步发电机q轴同步电抗,xf为发电机励磁绕组自感抗,xD为d轴阻尼绕组自感抗,xg为q轴等值g绕组自感抗,xQ为q轴等值Q绕组自感抗,xaf为发电机d轴定子绕组与励磁绕组间互感抗,xaD为发电机d轴定子绕组与阻尼绕组间互感抗,xag为发电机q轴定子绕组与q轴等值g绕组间互感抗,xaQ为发电机q轴定子绕组与q轴等值Q绕组互感抗,xfD为发电机励磁绕组与d轴阻尼绕组间互感抗,xgQ为q轴等值g绕组与等值Q绕组间的互感抗;
由于采用Xad基值系统,则有
xaf=xaD=xad(20)。
xag=xaQ=xaq
同步发电机实用模型参数分为稳态和发生扰动两种情况进行辨识;
A)在稳态情况下,同步发电机的阻尼绕组电流为0,同步发电机的初始稳态与扰动后稳态的机端电压和电流有名值方程如下:
其中,Uq0为初始稳态时q轴励磁电压有名值,Uq∞为扰动后q轴励磁电压有名值,Ud0为初始稳态时d轴励磁电压有名值,Ud∞为扰动后d轴励磁电压有名值,Ra为定子电阻有名值,Iq0为初始稳态时q轴定子电流有名值,Iq∞为扰动后q轴定子电流有名值,Id0为初始稳态时d轴定子电流有名值,Id0为扰动后d轴定子电流有名值,Xd为d轴同步电抗有名值、Xq为q轴同步电抗有名值,Xaf为d轴定子与励磁绕组互感抗有名值,If0为初始稳态时励磁电流有名值,If∞为扰动后励磁电流有名值;
通过方程(21)求出Ra、Xd、Xq和Xaf,再利用初始稳态时励磁电流与励磁电压的关系,即可计算Rf:
其中,Uf0为初始稳态时励磁电压有名值;
B)扰动后,假设励磁系统电压基值为K,作为待辨识的量,则有:
其中:Sn为同步发电机额定容量,UB为同步发电机相电压峰值,均为已知量。
发生扰动后,辨识过程包括d轴参数辨识、q轴参数辨识和同步发电机转子运动方程参数辨识,具体如下:
a)轴参数辨识
根据磁链方程推出:
进而有
在Uf实测并采用有名值的情况下:
因此建立状态方程:
量测方程为:
b)q轴参数辨识
假设同步发电机q轴有两个阻尼绕组,则阻尼绕组电流为:
因此建立状态方程:
量测方程为:
c)同步发电机转子运动方程参数辨识
其中,D为风阻系数,通常取0;Δδ(t)为同步发电机在t时刻的功角变化量;
选择Δδ、Δω作为状态量,选择Δδ(t)作为观测量,利用最小二乘法或遗传算法进行参数拟合,求出辨识参数Tj和D。
利用xf、xD、K、rD、xQ、xg、xag、rQ和rg根据所述辨识参数计算同步发电机实用模型参数:
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.一种基于PMU的同步发电机实用模型参数校核与辨识方法,其特征在于:所述方法包括以下步骤:
步骤1:校核同步发电机实用模型参数;
步骤2:辨识同步发电机实用模型参数;
所述步骤1包括以下步骤:
步骤1‐1:提取PMU记录的动态数据;所述动态数据包括有功功率P、无功功率Q、高压母线电压幅值U、高压母线电压相角和励磁电压Uf;
步骤1‐2:将同步发电机组与系统解耦,形成等值系统;所述等值系统用同步发电机组微分方程、机端母线与高压母线的代数方程联立表示;
在考虑同步发电机阻尼绕组但不计饱和影响的情况下,同步发电机电压方程为
其中:x′d为同步发电机d轴暂态同步电抗,x″d为同步发电机d轴次暂态同步电抗,x'q为同步发电机q轴暂态同步电抗,x″q为同步发电机q轴次暂态同步电抗,T′d0为同步发电机d轴暂态开路时间常数,T″d0为同步发电机d轴次暂态开路时间常数,T′q0为同步发电机q轴暂态开路时间常数,T″q0为同步发电机q轴次暂态开路时间常数,且:
kd=(xd-x'd)/(x'd-x″d) (3)
kq=(xq-x'q)/(x'q-x″q) (4)
其中,E'q为同步发电机q轴暂态电势,E"q为同步发电机q轴次暂态电势,E'd为同步发电机d轴暂态电势,E"d为同步发电机d轴次暂态电势,id为同步发电机d轴定子电流,iq为同步发电机q轴定子电流,xd为同步发电机d轴同步电抗,xq为同步发电机q轴同步电抗,Uf为同步发电机励磁电压实测有名值,Rf为励磁绕组电阻有名值,If0为初始稳态时励磁电流有名值,
运用隐式梯形积分公式,则有:
其中:
其中,Δt为仿真步长;
忽略转速变化和定子回路电磁暂态过程影响,在t+Δt时刻的定子电压平衡方程为:
其中:ra为同步发电机定子电阻,且
将式(9)变换到XY坐标可得:
其中:
ud(t)为t时刻同步发电机d轴电压分量,uq(t)为t时刻同步发电机q轴电压分量,id(t)为t时刻同步发电机d轴定子电流,iq(t)为t时刻同步发电机q轴定子电流,δ(t)为同步发电机在t时刻的功角,为t时刻同步发电机定子电压相量,为t时刻同步发电机定子电流相量;
同步发电机转子运动方程为:
其中,Pm(t)为同步发电机t时刻机械功率,通常取0;Pe(t)为同步发电机t时刻电磁功率,Pe(t)=ud(t)*id(t)+uq(t)*iq(t),且
其中,ωs为同步发电机的同步机械角速度,Tj为同步发电机转子的转动惯量,ω(t)为同步发电机t时刻电角速度标幺值;
同步发电机端电压与高压母线电压代数方程为:
其中,RT+jXT为同步发电机出口变压器的阻抗,为高压母线电压相量;
联立(11)、(13)和(15),将带入实测值作为已知量计算出id(t+Δt)、ud(t+Δt)、uq(t+Δt)和iq(t+Δt),进而计算同步发电机输出功率:
有功功率P、无功功率Q与实测有功功率Pm、实测无功功率Qm之间误差反映同步发电机参数描述同步发电机外特性的准确性,进而计算同布发电机参数误差Err,其表示为:
其中i为记录数据的长度,Δt′为采样间隔,且Δt′=Δt;
当误差Err不满足要求时,同步发电机参数则重新进行辨识。
2.根据权利要求1所述的基于PMU的同步发电机实用模型参数校核与辨识方法,其特征在于:所述步骤2中,采用标幺值表示同步发电机PARK电压方程:
其中,ud为同步发电机d轴电压分量,uq为同步发电机q轴电压分量,uf为同步发电机励磁电压,rf为d轴励磁绕组电阻,rD为d轴阻尼绕组电阻,rg为q轴等值g绕组电阻,rQ为q轴等值Q绕组电阻,if为励磁绕组电流,iD为d轴阻尼绕组电流,ig为q轴等值g绕组电流,iQ为q轴等值Q绕组电流,p为微分算子,ψd为d轴定子绕组磁链,ψf为发电机励磁绕组磁链,ψD为d轴阻尼绕组磁链,ψq为q轴定子绕组磁链,ψg为q轴等值g绕组磁链,ψQ为发电机q轴等值Q绕组磁链,ω为同步发电机角速度;
同步发电机PARK磁链方程表示为:
xd为同步发电机d轴同步电抗,xq为同步发电机q轴同步电抗,xf为发电机励磁绕组自感抗,xD为d轴阻尼绕组自感抗,xg为q轴等值g绕组自感抗,xQ为q轴等值Q绕组自感抗,xaf为发电机d轴定子绕组与励磁绕组间互感抗,xaD为发电机d轴定子绕组与阻尼绕组间互感抗,xag为发电机q轴定子绕组与q轴等值g绕组间互感抗,xaQ为发电机q轴定子绕组与q轴等值Q绕组互感抗,xfD为发电机励磁绕组与d轴阻尼绕组间互感抗,xgQ为q轴等值g绕组与等值Q绕组间的互感抗;
由于采用Xad基值系统,则有
xaf=xaD=xad
xag=xaQ=xaq (20)。
3.根据权利要求2所述的基于PMU的同步发电机实用模型参数校核与辨识方法,其特征在于:同步发电机实用模型参数分为稳态和发生扰动两种情况进行辨识;
A)在稳态情况下,同步发电机的阻尼绕组电流为0,同步发电机的初始稳态与扰动后稳态的机端电压和电流有名值方程如下:
其中,Uq0为初始稳态时q轴励磁电压有名值,Uq∞为扰动后q轴励磁电压有名值,Ud0为初始稳态时d轴励磁电压有名值,Ud∞为扰动后d轴励磁电压有名值,Ra为定子电阻有名值,Iq0为初始稳态时q轴定子电流有名值,Iq∞为扰动后q轴定子电流有名值,Id0为初始稳态时d轴定子电流有名值,Id∞为扰动后d轴定子电流有名值,Xd为d轴同步电抗有名值、Xq为q轴同步电抗有名值,Xaf为d轴定子与励磁绕组互感抗有名值,If0为初始稳态时励磁电流有名值,If∞为扰动后励磁电流有名值;
通过方程(21)求出Ra、Xd、Xq和Xaf,再利用初始稳态时励磁电流与励磁电压的关系,即可计算Rf:
其中,Rf为励磁绕组电阻有名值,Uf0为初始稳态时励磁电压有名值;
B)扰动后,假设励磁系统电压基值为K,作为待辨识的量,则有:
其中:Sn为同步发电机额定容量,UB为同步发电机相电压峰值,均为已知量。
4.根据权利要求3所述的基于PMU的同步发电机实用模型参数校核与辨识方法,其特征在于:发生扰动后,辨识过程包括d轴参数辨识、q轴参数辨识和同步发电机转子运动方程参数辨识,具体如下:
a)d轴参数辨识
根据磁链方程推出:
进而有
在Uf实测并采用有名值的情况下:
因此建立状态方程:
量测方程为:
b)q轴参数辨识
假设同步发电机q轴有两个阻尼绕组,则阻尼绕组电流为:
因此建立状态方程:
量测方程为:
c)同步发电机转子运动方程参数辨识
其中,D为风阻系数,通常取0;Δδ(t)为同步发电机在t时刻的功角变化量;
选择Δδ、Δω作为状态量,选择Δδ(t)作为观测量,利用最小二乘法或遗传算法进行参数拟合,求出辨识参数Tj和D。
5.根据权利要求4所述的基于PMU的同步发电机实用模型参数校核与辨识方法,其特征在于:利用xf、xD、K、rD、xQ、xg、xag、rQ和rg根据所述辨识参数计算同步发电机实用模型参数:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210397510.8A CN103036498B (zh) | 2012-10-18 | 2012-10-18 | 一种基于pmu的同步发电机实用模型参数校核与辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210397510.8A CN103036498B (zh) | 2012-10-18 | 2012-10-18 | 一种基于pmu的同步发电机实用模型参数校核与辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103036498A CN103036498A (zh) | 2013-04-10 |
CN103036498B true CN103036498B (zh) | 2015-07-22 |
Family
ID=48023057
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210397510.8A Active CN103036498B (zh) | 2012-10-18 | 2012-10-18 | 一种基于pmu的同步发电机实用模型参数校核与辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103036498B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103699723B (zh) * | 2013-12-09 | 2016-04-13 | 国家电网公司 | 一种发电厂机组动力系统模型校核方法 |
CN104967120B (zh) * | 2015-07-08 | 2017-07-14 | 中国电力科学研究院 | 一种基于恒功率接口的混合动态仿真方法 |
US10103666B1 (en) | 2015-11-30 | 2018-10-16 | University Of South Florida | Synchronous generator modeling and frequency control using unscented Kalman filter |
CN106936359B (zh) * | 2015-12-29 | 2021-03-02 | 中国电力科学研究院 | 一种机电暂态同步机模型及初始q轴饱和参数的确定方法 |
CN106599541B (zh) * | 2016-11-23 | 2019-04-19 | 华南理工大学 | 一种动态电力负荷模型的结构和参数在线辨识方法 |
CN107677960B (zh) * | 2017-09-21 | 2020-01-10 | 国网福建省电力有限公司 | 基于pmu量测的同步发电机参数辨识方法 |
CN107885906A (zh) * | 2017-10-16 | 2018-04-06 | 中国农业大学 | 一种基于遗传算法的电机系统能耗计算方法 |
US10985611B2 (en) | 2019-04-10 | 2021-04-20 | General Electric Company | System and method for estimating grid strength |
CN110346718A (zh) * | 2019-06-20 | 2019-10-18 | 广西电网有限责任公司电力科学研究院 | 一种同步发电机d轴参数测试及辨识方法 |
WO2020258202A1 (zh) * | 2019-06-28 | 2020-12-30 | 瑞声声学科技(深圳)有限公司 | 马达参数追踪方法及系统 |
CN110729936B (zh) * | 2019-07-25 | 2021-11-09 | 国网福建省电力有限公司电力科学研究院 | 一种计及饱和特性的发电机实用参数辨识模型及其辨识方法 |
CN113688585A (zh) * | 2021-08-03 | 2021-11-23 | 国电南瑞科技股份有限公司 | 一种基于模型参考的参数校核方法、装置及系统 |
CN117349576B (zh) * | 2023-12-05 | 2024-03-01 | 华能澜沧江水电股份有限公司 | 一种发电机投电制动转子过电压计算方法和系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101629995A (zh) * | 2009-08-10 | 2010-01-20 | 华中电网有限公司 | 相量测量单元pmu的检测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2904163B1 (fr) * | 2006-07-20 | 2008-09-12 | Schneider Toshiba Inverter | Procede d'ajustement de parametres d'un moteur synchrone et variateur de vitesse utilisant un tel procede |
-
2012
- 2012-10-18 CN CN201210397510.8A patent/CN103036498B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101629995A (zh) * | 2009-08-10 | 2010-01-20 | 华中电网有限公司 | 相量测量单元pmu的检测方法 |
Non-Patent Citations (2)
Title |
---|
同步发电机模型参数重要度等级评估;周成等;《电力科学与技术学报》;20120630;第27卷(第2期);5-9页 * |
韩超先.基于PMU数据的同步发电机参数辨识.《万方学位论文库(电子期刊)》.2012,第17- 41页. * |
Also Published As
Publication number | Publication date |
---|---|
CN103036498A (zh) | 2013-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103036498B (zh) | 一种基于pmu的同步发电机实用模型参数校核与辨识方法 | |
CN102510263B (zh) | 基于抛载试验和数值差分的同步发电机实用参数辨识方法 | |
CN107132450B (zh) | 一种海上双馈电机定子绕组匝间短路早期故障辨识方法 | |
CN107783006B (zh) | 水轮发电机转子绕组匝间短路故障检测方法 | |
CN103401238B (zh) | 一种基于总体测辨法的电力负荷建模方法 | |
CN104158189B (zh) | 一种基于参数辨识的电气化铁路牵引供电负荷建模方法 | |
CN103427751B (zh) | 永磁同步电机静态参数在线辨识的装置与方法 | |
CN102073012B (zh) | 同步发电机参数获取方法及系统 | |
CN103618492A (zh) | 一种基于时频变换的同步发电机参数辨识方法 | |
CN107590317A (zh) | 一种计及模型参数不确定性的发电机动态估计方法 | |
CN106443246A (zh) | 基于pmu量测数据的小干扰稳定参数的在线辨识方法 | |
CN103424734A (zh) | 电子式电能表在线校准方法 | |
CN102904518B (zh) | 一种同步发电机q轴参数在线辨识方法 | |
CN101566663A (zh) | 一种配电系统电压跌落源定位方法 | |
CN107565867A (zh) | 一种基于轨迹灵敏度的同步发电机参数辨识方法 | |
CN103023419B (zh) | 一种用于发电机同步电抗参数辨识的pmu数据的筛选方法 | |
CN104182644A (zh) | 一种融合分布式电源特性的配电网状态估计方法 | |
CN103995245A (zh) | 双馈风力发电机定转子电流信号检测系统的故障判断方法 | |
CN104101805B (zh) | 一种基于振荡能量注入的励磁系统负阻尼检测方法 | |
CN106610473B (zh) | 基于同步发电机静止频率响应试验的电气参数测试方法 | |
CN104393809B (zh) | 适用于scr静止变频器的抽水蓄能机组低速位置检测方法 | |
CN102280884A (zh) | 一种电网等值方法 | |
CN104267243B (zh) | 同步发电机内电势和电抗参数的测量方法及装置 | |
JP6161783B2 (ja) | コンピュータ支援により送配電網のインピーダンスを求める方法、当該方法を実施するための発電装置及びコンピュータプログラム | |
CN110601175B (zh) | 一种励磁电流改进标幺方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |