CN103595062A - 基于无模型自适应控制算法的发电机广域阻尼控制方法 - Google Patents

基于无模型自适应控制算法的发电机广域阻尼控制方法 Download PDF

Info

Publication number
CN103595062A
CN103595062A CN201310316663.XA CN201310316663A CN103595062A CN 103595062 A CN103595062 A CN 103595062A CN 201310316663 A CN201310316663 A CN 201310316663A CN 103595062 A CN103595062 A CN 103595062A
Authority
CN
China
Prior art keywords
phi
wide area
theta
control
pseudo
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.)
Granted
Application number
CN201310316663.XA
Other languages
English (en)
Other versions
CN103595062B (zh
Inventor
庞小艳
赵艺
陆超
李建
段刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Sifang Automation Co Ltd
Tsinghua University
State Grid Corp of China SGCC
State Grid Sichuan Electric Power Co Ltd
Original Assignee
Beijing Sifang Automation Co Ltd
Tsinghua University
State Grid Corp of China SGCC
State Grid Sichuan Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Sifang Automation Co Ltd, Tsinghua University, State Grid Corp of China SGCC, State Grid Sichuan Electric Power Co Ltd filed Critical Beijing Sifang Automation Co Ltd
Priority to CN201310316663.XA priority Critical patent/CN103595062B/zh
Publication of CN103595062A publication Critical patent/CN103595062A/zh
Application granted granted Critical
Publication of CN103595062B publication Critical patent/CN103595062B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)

Abstract

本发明涉及基于无模型自适应控制算法的发电机广域阻尼控制方法,属于电力系统稳定分析技术领域,该方法包括:根据传统的传递函数表达式及参数计算公式得到传统发电广域阻尼控制方法的控制参数;根据得到的传统发电广域阻尼控制方法的控制参数计算发电机广域阻尼自适应控制方法的控制参数初始值,以得到的伪梯度向量初始值为初始点,利用被控系统实测的输入输出数据在线更新伪梯度向量的估计值;根据在线更新的伪梯度向量的估计值计算发电机广域阻尼自适应控制信号;重复计算直至达到所需计算的发电机广域阻尼控制信号的最大时刻。该方法能够保证发电机广域阻尼控制在系统多种运行状态下的控制效果。

Description

基于无模型自适应控制算法的发电机广域阻尼控制方法
技术领域
本发明属于电力系统稳定分析技术领域,特别涉及基于无模型自适应控制算法的广域阻尼控制方法。
背景技术
随着电力系统的快速发展以及区域电网间的不断互联,区域电网间的区间低频振荡问题日益突出。传统抑制区间低频振荡的方法是在电力系统发电机励磁侧施加电力系统稳定控制(Power Stability Stabilizer,PSS)。然而,由于PSS所采用的本地反馈信号对区间振荡模式的可观性不足,PSS对区间低频振荡的抑制效果有限。
广域测量技术的成熟,能够实现对全网广域信号的实时采集,广域信号对区间低频振荡模式的可观性明显优于本地信号,因此对于抑制区间低频振荡具有极大的潜力。这种利用广域测量系统测得的广域信号作为PSS反馈信号,抑制区间低频振荡的技术,称为发电机广域阻尼控制。传统的发电广域阻尼控制方法通常是在电力系统基础运行方式下,生成随机小幅扰动信号,采集广域阻尼控制的反馈信号;利用随机小幅扰动信号和发电机广域阻尼控制的反馈信号和扰动信号通过预报误差方法对电力系统离线辨识获得线性模型;利用线性模型获得的留数的表达式得到在电力系统基础运行方式下的移相角度和增益;利用移相角度和增益计算出传统发电机广域阻尼控制参数;利用所述控制参数对电力系统的发电机的区间低频振荡进行控制。
在实际运行中,电力系统的运行方式和结构呈现出时变的特点,传统的发电机广域阻尼控制由于其参数固定而无法保证其在电力系统的多种运行方式下都能取得较好的控制效果。目前为解决发电机广域阻尼控制的自适应问题,通常通过对电力系统在线辨识获得线性模型完成发电机广域阻尼控制。然而,电力系统本身是一个高阶复杂的非线性系统,而且在发生大的扰动或故障时,电力系统呈现出强非线性,在线辨识得到的电力系统的线性模型与实际电力系统偏差很大,基于在线辨识得到的电力系统的线性模型的控制方法往往不能保证其控制效果。
无模型自适应控制算法在20世纪90年代首次被提出,该算法适用于非线性系统,仅利用系统的输入输出数据即可实现对控制参数的自适应调整,避免了对被控系统建模的过程,为实现发电机广域阻尼自适应控制方法提供了新的解决思路。
发明内容
本发明的目的是为克服已有技术的不足之处,提出一种基于无模型自适应控制算法的发电机广域阻尼控制方法,该方法能够适应电力系统的非线性特性,利用系统的输入输出数据实现对控制参数的自适应调整,避免了对被控系统建模的过程,能够保证发电机广域阻尼控制在系统多种运行状态下的控制效果。
本发明提出的基于无模型自适应控制算法的发电机广域阻尼控制方法,其特征在于,该方法包括:
1)根据传统发电机广域阻尼控制方法的传递函数表达式及参数计算公式得到传统发电广域阻尼控制方法的控制参数;
2)定义被控电力系统的动态线性表达式、控制准则函数表达式和伪梯度向量估计准则函数表达式;
3)然后根据传统发电机广域阻尼控制方法得到的控制参数计算发电机广域阻尼自适应控制方法的控制参数初始值,包括伪梯度向量的初始值
Figure BDA00003569026700025
惩罚因子γ;
4)以步骤3)中得到伪梯度向量初始值为初始点,利用被控系统实测的输入输出数据在线更新伪梯度向量的估计值
Figure BDA00003569026700027
5)最后根据在线更新的伪梯度向量的估计值计算发电机广域阻尼自适应控制信号;
6)将步骤5)中的发电机广域阻尼自适应控制信号作用于被控电力系统的发电机励磁端,k增加1,重复步骤4)~6)直至达到所需计算的发电机广域阻尼控制信号的最大时刻。
本发明的特点及有益效果:本发明的特点在于,所述方法考虑了电力系统的非线性特点,能够有效避免对被控电力系统的建模,利用电力系统实测的输入输出数据实现了对发电机广域阻尼控制参数的自适应调整。本发明所提出的广域阻尼自适应方法,能够使得电力系统在运行方式和结构发生变化的情况下,仍然能够有效的抑制区间低频振荡,提高电力系统的稳定性。
附图说明
图1为本发明的基于无模型自适应控制算法的发电机广域阻尼控制方法实现流程框图;
图2为本发明的实施例中四机两区被控电力系统示意图;
图3为本实施例的方法与传统方法控制效果比较图。
具体实施方式
本发明提出的基于无模型自适应控制算法的发电机广域阻尼控制方法,结合附图及实施例详细说明如下:
本发明方法采用无模型自适应控制算法实现了发电机广域阻尼控制的自适应调整,使被控电力系统在多种运行方式下阻尼比均能满足要求。
本发明的基于无模型自适应控制算法的发电机广域阻尼控制方法,其特征在于,该方法包括:
1)根据传统发电机广域阻尼控制方法的传递函数表达式及参数计算公式得到传统发电广域阻尼控制方法的控制参数;
2)定义被控电力系统的动态线性表达式、控制准则函数表达式和伪梯度向量估计准则函数表达式;
3)然后根据传统发电机广域阻尼控制方法得到的控制参数计算发电机广域阻尼自适应控制方法的控制参数初始值,包括伪梯度向量的初始值
Figure BDA00003569026700036
惩罚因子γ;
4)进一步以步骤3)中得到伪梯度向量初始值为初始点,利用被控系统实测的输入输出数据在线更新伪梯度向量的估计值
Figure BDA00003569026700037
5)最后根据在线更新的伪梯度向量的估计值计算发电机广域阻尼自适应控制信号;
6)将步骤5)中的发电机广域阻尼自适应控制信号作用于被控电力系统的发电机励磁端,k增加1,重复步骤4)~6)直至达到所需计算的发电机广域阻尼控制信号的最大时刻。
本发明方法的具体实现流程如图1所示,具体包括以下步骤:
步骤1)根据传统发电机广域阻尼控制方法的传递函数表达式及参数计算公式得到传统发电广域阻尼控制方法的控制参数;具体实现方式为:
传统发电机广域阻尼控制方法的传递函数表达式如式(1)所示:
K C ( s ) = K c · ( 1 + T 1 s 1 + T 2 s ) 2 - - - ( 1 )
其中T1为移相环节的超前时间常数,T2为移相环节的滞后时间常数,Kc为发电机广域阻尼控制的增益;
参数计算公式如式(2)所示:
Figure BDA00003569026700032
其中∠K(jωd)为传统发电机广域阻尼控制方法得到的控制移相角度,|KC(jωd)|为传统发电机广域阻尼控制方法得到的控制增益,ωn为区间低频振荡模式的无阻尼自然振荡频率(已知量)。
步骤2)定义被控电力系统的动态线性表达式、控制准则函数表达式和伪梯度向量估计准则函数表达式;具体实现方式为:
所述动态线性表达式如式(3)所示:
y(k+1)=y(k)+φT(k)ψ(k)   (3)
其中φ(k)=[φ1(k)LφLu(k)φLu+1(k)LφLu+Ly(k)]T为系统的伪梯度向量,ψ(k)=[Δu(k)LΔu(k-Lu+1)Δy(k)LΔy(k-Ly+1)]T,u(k)为需要计算的发电机广域阻尼控制信号,y(k)为区间频差信号,Lu和Ly为电力系统的动态线性表达式的伪阶数,k为控制信号作用的时刻,k=1,2LTtotal1,Ttotal1是需要计算的广域阻尼控制信号的最大时刻(已知量);
所述控制准则函数表达式如式(4)所示:
J C ( u ( k ) ) = | y * ( k + 1 ) - y ( k ) - φ ^ T ( k ) ψ ( k ) | 2 + γ | u ( k ) - u ( k - 1 ) | 2 - - - ( 4 )
其中y*(k+1)为k+1时刻区间频差信号的期望值,γ为u(k)变化量的惩罚因子,其中
Figure BDA00003569026700042
为伪梯度向量在k时刻的估计值;
所述伪梯度向量估计准则函数表达式如式(5)所示:
J ( θ ^ ( k ) ) = | y ( k ) - y ( k - 1 ) - φ ^ T ( k ) ψ ( k ) | 2 + Λ μ | φ ^ ( k ) - φ ^ ( k - 1 ) | 2 - - - ( 5 )
其中Λμ变化量的对角权重矩阵。
步骤3)然后根据传统发电机广域阻尼控制方法得到的控制参数计算发电机广域阻尼自适应控制方法的控制参数初始值,包括伪梯度向量的初始值
Figure BDA00003569026700045
惩罚因子γ;具体实现方式为:
取伪阶数Lu=3和Ly=3,通过极小化控制准则函数得到发电机广域阻尼自适应控制方法的控制律表达式为:
u ( k ) = u ( k - 1 ) + ρ k φ ^ u 1 ( k ) γ + | φ ^ u 1 ( k ) | 2 [ y * ( k + 1 ) - y ( k ) - φ ^ y 1 ( k ) Δy ( k ) - φ ^ y 2 ( k ) Δy ( k - 1 ) - - - ( 6 )
- φ ^ y 3 ( k ) Δy ( k - 2 ) - φ ^ u 2 ( k ) Δu ( k - 1 ) - φ ^ u 3 ( k ) Δu ( k - 2 ) ]
其中ρk为步长序列;
为使区间低频振荡尽快平息,将k+1时刻区间频差信号的期望值y*(k+1)取为0;
Figure BDA00003569026700048
所述控制律的离散传递函数表达式为:
K M ( z ) = θ e ( k ) { [ 1 + φ ^ y 1 ( k ) ] z 3 + [ φ ^ y 2 ( k ) - φ ^ y 1 ( k ) ] z 2 + [ φ ^ y 3 ( k ) - φ ^ y 2 ( k ) ] z - φ ^ y 3 ( k ) } z 3 + ( θ e ( k ) φ ^ u 2 ( k ) - 1 ) z 2 + θ e ( k ) ( φ ^ u 3 ( k ) - φ ^ u 2 ( k ) ) z - θ e ( k ) φ ^ u 3 ( k ) - - - ( 7 )
= θ e ( k ) { [ 1 + φ ^ y 1 ( k ) ] z 3 + [ φ ^ y 2 ( k ) - φ ^ y 1 ( k ) ] z 2 + [ φ ^ y 3 ( k ) - φ ^ y 2 ( k ) ] z - φ ^ y 3 ( k ) } ( z - 1 ) [ z 2 + ( θ e ( k ) φ ^ u 2 ( k ) z + θ e ( k ) φ ^ u 3 ( k ) ]
离散传递函数KM(z)由三个零点和三个极点组成,其中一个极点为1;
(为使在k=1时刻采用发电机广域阻尼自适应控制方法的控制参数与采用传统广域阻尼控制方法的控制参数得到的控制效果相同,)将k=1时刻KM(z)的零极点取值与KC(s)的零极点相同;其中k=1时刻KM(z)的一个零点配置在极点1附近,另外两个零点和两个极点分别与KC(s)的零点和极点相同;
假设x01,x02,x03为表达式 [ 1 + φ ^ y 1 ( 1 ) ] z 3 + [ φ ^ y 2 ( 1 ) - φ ^ y 1 ( 1 ) ] z 2 + [ φ ^ y 3 ( 1 ) - φ ^ y 2 ( 1 ) ] z - φ ^ y 3 ( 1 ) = 0 的三个根;x01为在极点1附近的零点,(为减小零点x01和极点1对振荡模式的动态性能影响,)取该零点到原点的距离为区间低频振荡模式对应的特征根到虚轴距离的1/5,即
Figure BDA00003569026700054
其中Ts为采样周期,ξ0为区间低频振荡模式的阻尼比。x02=x03与KC(s)中的零点相同,即 x 02 = x 03 = e - 1 T 1 T s ;
令z=1,得到 ( 1 - x 01 ) ( 1 - x 02 ) 2 = 1 1 + φ ^ y 1 ( 1 ) ;
(由于x01,x02和x03已知,)根据上述等式得到计算
Figure BDA00003569026700057
的表达式如式(8)所示:
φ ^ y 1 ( 1 ) = 1 ( 1 - x 01 ) ( 1 - x 02 ) 2 - 1 - - - ( 8 )
(由于x01逼近1,故(1-x01)取值很小,对应的
Figure BDA00003569026700059
取值很大,)则
Figure BDA000035690267000510
KM(z)化简为:
K M ( z ) = θ e ( 1 ) { [ 1 + φ ^ y 1 ( 1 ) ] z 3 + [ φ ^ y 2 ( 1 ) - φ ^ y 1 ( 1 ) ] z 2 + [ φ ^ y 3 ( 1 ) - φ ^ y 2 ( 1 ) ] z - φ ^ y 3 ( 1 ) } z 3 + ( θ e ( 1 ) φ ^ u 2 ( 1 ) - 1 ) z 2 + θ e ( 1 ) ( φ ^ u 3 ( 1 ) - φ ^ u 2 ( 1 ) ) z - φ ^ u 3 ( 1 )
≈ θ e ( 1 ) { φ ^ y 1 ( 1 ) z 3 + [ φ ^ y 2 ( 1 ) - φ ^ y 1 ( 1 ) ] z 2 + [ φ ^ y 3 ( 1 ) - φ ^ y 2 ( 1 ) ] z - φ ^ y 3 ( 1 ) } z 3 + ( θ e ( 1 ) φ ^ u 2 ( 1 ) - 1 ) z 2 + θ e ( 1 ) ( φ ^ u 3 ( 1 ) - φ ^ u 2 ( 1 ) ) z - θ e ( 1 ) φ ^ u 3 ( 1 ) - - - ( 9 )
= θ e ( 1 ) [ φ ^ y 1 ( 1 ) z 2 + φ ^ y 2 ( 1 ) z + φ ^ y 3 ( 1 ) ] [ z 2 + θ e ( 1 ) φ ^ u 2 ( 1 ) z + θ e ( 1 ) φ ^ u 3 ( 1 ) ]
根据根与系数的关系式得到如下关系式:
- φ ^ y 2 ( 1 ) φ ^ y 1 ( 1 ) = 2 x 02 - - - ( 10 )
φ ^ y 3 ( 1 ) φ ^ y 1 ( 1 ) = ( x 02 ) 2 - - - ( 11 )
由式(10)、(11)得到计算
Figure BDA000035690267000627
的表达式如式(12)、(13)所示:
φ ^ y 2 ( 1 ) = - 2 x 02 · φ ^ y 1 ( 1 ) - - - ( 12 )
φ ^ y 3 ( 1 ) = φ ^ y 1 ( 1 ) · ( x 02 ) 2 - - - ( 13 )
假设极点xp2、xp3
Figure BDA00003569026700067
的两个根,取xp2、xp3与KC(s)中的极点相同,即 x p 2 = x p 3 = e - 1 T 2 T s ;
根据根与系数的关系式得到如下关系式:
- θ e ( 1 ) φ ^ u 2 ( 1 ) 2 = x p 2 - - - ( 14 )
θ e ( 1 ) φ ^ u 3 ( 1 ) = ( x p 2 ) 2 - - - ( 15 )
由式(14)、(15)得到计算
Figure BDA000035690267000628
的表达式如式(16)、(17)所示:
θ e ( 1 ) φ ^ u 2 ( 1 ) = - 2 x p 2 - - - ( 16 )
θ e ( 1 ) φ ^ u 3 ( 1 ) = ( x p 2 ) 2 - - - ( 17 )
| K C ( s ) | | s = j ω d = | K M ( z ) | | z = e j ω d T s , 得到如下等式:
K c = θ e ( 1 ) [ φ ^ y 1 ( 1 ) + φ ^ y 2 ( 1 ) φ ^ y 3 ( 1 ) ] [ 1 + θ e ( 1 ) φ ^ u 2 ( 1 ) + θ e ( 1 ) φ ^ u 3 ( 1 ) ] - - - ( 18 )
式(18)中
Figure BDA000035690267000629
Figure BDA000035690267000621
均已知,得到θe(1)的计算表达式如式(19)所示:
θ e ( 1 ) = K c · [ 1 + θ e ( 1 ) φ ^ u 2 ( 1 ) + θ e ( 1 ) φ ^ u 3 ( 1 ) ] [ φ ^ y 1 ( 1 ) + φ ^ y 2 ( 1 ) + φ ^ y 3 ( 1 ) ] - - - ( 19 )
根据式(19)得到的θe(1)的取值得到
Figure BDA000035690267000623
Figure BDA000035690267000624
的取值;
由于
Figure BDA000035690267000625
其中ρk和γ按照经验值选取,ρk通常在0.5~1之间取值,γ为一个大于1的值,根据θe(1)的表达式算出的取值。
步骤4)进一步以步骤3)中得到伪梯度向量初始值
Figure BDA00003569026700071
为初始点,利用被控系统实测的输入输出数据在线更新伪梯度向量的估计值
Figure BDA00003569026700072
具体实现方式为:
通过极小化伪梯度向量估计准则函数得到,伪梯度向量φ(k)在k时刻的在线更新计算表达式:
φ ^ ( k ) = φ ^ ( k - 1 ) +
[ ψ ( k - 1 ) ψ T ( k - 1 ) + μ 1 O μ L u O μ L u + L y ] - 1 · ψ ( k - 1 ) [ Δy ( k ) - ψ T ( k - 1 ) φ ^ ( k - 1 ) ] - - - ( 20 )
其中Λμ为对角权重矩阵, Λ μ = μ 1 O μ L u O μ L u + L y , 用于限制
Figure BDA00003569026700076
中每个伪梯度向量元素的变化;
Λμ根据经验值选取。
步骤5)最后根据在线更新的伪梯度向量的估计值计算发电机广域阻尼自适应控制信号;具体实现方式为:
根据步骤4)中k时刻在线更新的伪梯度向量估计值
Figure BDA00003569026700077
利用发电机广域阻尼自适应控制方法的控制律计算k时刻发电机广域阻尼自适应控制方法的控制信号:
u ( k ) = u ( k - 1 ) + ρ k φ ^ u 1 ( k ) γ + | φ ^ u 1 ( k ) | 2 [ y * ( k + 1 ) - y ( k ) - φ ^ y 1 ( k ) Δy ( k ) - φ ^ y 2 ( k ) Δy ( k - 1 ) - - - ( 21 )
- φ ^ y 3 ( k ) Δy ( k - 2 ) - φ ^ u 2 ( k ) Δu ( k - 1 ) - φ ^ u 3 ( k ) Δu ( k - 2 ) ]
步骤6)将步骤5)中的发电机广域阻尼自适应控制信号作用于被控电力系统的发电机励磁端,k增加1,重复步骤4)~6)直至达到所需计算的发电机广域阻尼控制信号的最大时刻。
实施例:四机两区系统
以四机两区系统为例来说明上述发电机广域阻尼自适应控制方法的有效性,该系统结构如图3所示。
该系统中有区域1和区域2两个区域组成,每个区域有两台发电机,两个区域通过区间联络线3、4相连,区间输送功率为413MW;系统中四台发电机完全相同,额定容量为900MW,额定电压为20kV。在发电机4的励磁端施加发电机广域自适应阻尼控制,发电机广域阻尼控制的广域反馈信号为区域2和区域1的母线频差信号。
本实施例对上述四机两区系统的发电机广域阻尼自适应控制方法,包括以下步骤:
步骤1)根据传统发电机广域阻尼控制方法的传递函数表达式及参数计算公式得到传统发电广域阻尼控制方法的控制参数;具体实现方式为:
传统发电机广域阻尼控制方法的传递函数表达式如式(1)所示:
K C ( s ) = K c · ( 1 + T 1 s 1 + T 2 s ) 2 - - - ( 1 )
其中T1为移相环节的超前时间常数,T2为移相环节的滞后时间常数,Kc为发电机广域阻尼控制的增益;
参数计算公式如式(2)所示:
Figure BDA00003569026700082
其中∠K(jωd)为25°,|KC(jωd)|为0.1,ωd为3.672rad/s;
根据式(2)计算得到:
θmax=25°/2=12.5°;
Figure BDA00003569026700083
T 2 = 1 ω n α = 1 3.672 × 1.55 = 0.22 ;
T1=αT1=0.34;
K c = 0.1 | | ( 0.34 s + 1 0.22 s + 1 ) 2 | | | s = j 3.67 = 0.064 ;
步骤2)定义被控电力系统的动态线性表达式、控制准则函数表达式和伪梯度向量估计准则函数表达式;具体实现方式为:
所述动态线性表达式如式(3)所示:
y(k+1)=y(k)+φT(k)ψ(k)   (3)
其中φ(k)=[φ1(k)LφLu(k)φLu+1(k)LφLu+Ly(k)]T为系统的伪梯度向量,ψ(k)=[Δu(k)LΔu(k-Lu+1)Δy(k)LΔy(k-Ly+1)]T,u(k)为需要计算的发电机广域阻尼控制信号,y(k)为区间频差信号,Lu和Ly为电力系统的动态线性表达式的伪阶数,k为控制信号作用的时刻,k=1,2LTtotal1,Ttotal1是需要计算的广域阻尼控制信号的最大时刻,Ttotal1=1000;
所述控制准则函数表达式如式(4)所示:
J C ( u ( k ) ) = | y * ( k + 1 ) - y ( k ) - φ ^ T ( k ) ψ ( k ) | 2 + γ | u ( k ) - u ( k - 1 ) | 2 - - - ( 4 )
其中y*(k+1)为k+1时刻区间频差信号的期望值,γ为u(k)变化量的惩罚因子,
Figure BDA00003569026700092
为伪梯度向量在k时刻的估计值;
所述伪梯度向量估计准则函数表达式如式(5)所示:
J ( θ ^ ( k ) ) = | y ( k ) - y ( k - 1 ) - φ ^ T ( k ) ψ ( k ) | 2 + Λ μ | φ ^ ( k ) - φ ^ ( k - 1 ) | 2 - - - ( 5 )
其中Λμ
Figure BDA00003569026700094
变化量的对角权重矩阵;
步骤3)然后根据传统发电机广域阻尼控制方法得到的控制参数计算发电机广域阻尼自适应控制方法的控制参数初始值,包括伪梯度向量的初始值
Figure BDA00003569026700095
惩罚因子γ;具体实现方式为:
取伪阶数Lu=3和Ly=3,通过极小化控制准则函数得到发电机广域阻尼自适应控制方法的控制律表达式为:
u ( k ) = u ( k - 1 ) + ρ k φ ^ u 1 ( k ) γ + | φ ^ u 1 ( k ) | 2 [ y * ( k + 1 ) - y ( k ) - φ ^ y 1 ( k ) Δy ( k ) - φ ^ y 2 ( k ) Δy ( k - 1 ) - - - ( 6 )
- φ ^ y 3 ( k ) Δy ( k - 2 ) - φ ^ u 2 ( k ) Δu ( k - 1 ) - φ ^ u 3 ( k ) Δu ( k - 2 ) ]
其中ρk为步长序列;
为使区间低频振荡尽快平息,将k+1时刻区间频差信号的期望值y*(k+1)取为0;
Figure BDA000035690267000911
所述控制律的离散传递函数表达式为:
K M ( z ) = θ e ( k ) { [ 1 + φ ^ y 1 ( k ) ] z 3 + [ φ ^ y 2 ( k ) - φ ^ y 1 ( k ) ] z 2 + [ φ ^ y 3 ( k ) - φ ^ y 2 ( k ) ] z - φ ^ y 3 ( k ) } z 3 + ( θ e ( k ) φ ^ u 2 ( k ) - 1 ) z 2 + θ e ( k ) ( φ ^ u 3 ( k ) - φ ^ u 2 ( k ) ) z - θ e ( k ) φ ^ u 3 ( k ) - - - ( 7 )
= θ e ( k ) { [ 1 + φ ^ y 1 ( k ) ] z 3 + [ φ ^ y 2 ( k ) - φ ^ y 1 ( k ) ] z 2 + [ φ ^ y 3 ( k ) - φ ^ y 2 ( k ) ] z - φ ^ y 3 ( k ) } ( z - 1 ) [ z 2 + ( θ e ( k ) φ ^ u 2 ( k ) z + θ e ( k ) φ ^ u 3 ( k ) ]
离散传递函数KM(z)由三个零点和三个极点组成,其中一个极点为1;
(为使在k=1时刻采用发电机广域阻尼自适应控制方法的控制参数与采用传统广域阻尼控制方法的控制参数得到的控制效果相同,)将k=1时刻KM(z)的零极点取值与KC(s)的零极点相同;其中k=1时刻KM(z)的一个零点配置在极点1附近,另外两个零点和两个极点分别与KC(s)的零点和极点相同;
假设x01,x02,x03为表达式 [ 1 + φ ^ y 1 ( 1 ) ] z 3 + [ φ ^ y 2 ( 1 ) - φ ^ y 1 ( 1 ) ] z 2 + [ φ ^ y 3 ( 1 ) - φ ^ y 2 ( 1 ) ] z - φ ^ y 3 ( 1 ) = 0 的三个根;x01为在极点1附近的零点,(为减小零点x01和极点1对振荡模式的动态性能影响,)取该零点到原点的距离为区间低频振荡模式对应的特征根到虚轴距离的1/5,即
Figure BDA00003569026700102
其中Ts=0.1s,ξ0=7.576%,ωd=3.662。x02=x03与KC(s)中的零点相同,即 x 02 = x 03 = e - 1 T 1 T s = e - 1 0.34 × 0.1 = 0.7452 ;
令z=1,得到 1 1 + φ ^ y 1 ( 1 ) = ( 1 - x 01 ) ( 1 - x 02 ) 2 = ( 1 - 0.9945 ) ( 1 - 0.7452 ) 2 ;
根据上述等式得到计算
Figure BDA00003569026700105
的表达式如式(8)所示:
φ ^ y 1 ( 1 ) = 1 ( 1 - x 01 ) ( 1 - x 02 ) 2 - 1 = 2800 - - - ( 8 )
由于KM(z)化简为:
K M ( z ) = θ e ( 1 ) { [ 1 + φ ^ y 1 ( 1 ) ] z 3 + [ φ ^ y 2 - φ ^ y 1 ( 1 ) ] z 2 + [ φ ^ y 3 ( 1 ) - φ ^ y 2 ( 1 ) ] z - φ ^ y 3 ( 1 ) } z 3 + ( θ e ( 1 ) φ ^ u 2 ( 1 ) - 1 ) z 2 + θ e ( 1 ) ( φ ^ u 3 ( 1 ) - φ ^ u 2 ( 1 ) ) z - φ ^ u 3 ( 1 )
≈ θ e ( 1 ) { φ ^ y 1 ( 1 ) z 3 + [ φ ^ y 2 ( 1 ) - φ ^ y 1 ( 1 ) ] z 2 + [ φ ^ y 3 ( 1 ) - φ ^ y 2 ( 1 ) ] z - φ ^ y 3 ( 1 ) } z 3 + ( θ e ( 1 ) φ ^ u 2 ( 1 ) - 1 ) z 2 + θ e ( 1 ) ( φ ^ u 3 ( 1 ) - φ ^ u 2 ( 1 ) ) z - θ e ( 1 ) φ ^ u 3 ( 1 ) - - - ( 9 )
= θ e ( 1 ) [ φ ^ y 1 ( 1 ) z 2 + φ ^ y 2 ( 1 ) z + φ ^ y 3 ( 1 ) ] [ z 2 + θ e ( 1 ) φ ^ u 2 ( 1 ) z + θ e ( 1 ) φ ^ u 3 ( 1 ) ]
根据根与系数的关系式得到如下关系式:
- φ ^ y 2 ( 1 ) φ ^ y 1 ( 1 ) = 2 × e - 1 0.34 × 0.1 - - - ( 10 )
φ ^ y 3 ( 1 ) φ ^ y 1 ( 1 ) = ( x 02 ) 2 = e - 1 0.34 × 0.1 × 2 = 0.555 - - - ( 11 )
由式(10)、(11)得到计算
Figure BDA000035690267001013
的表达式如式(12)、(13)所示:
φ ^ y 2 ( 1 ) = - 2 x 02 · φ ^ y 1 ( 1 ) = - 4173 - - - ( 12 )
φ ^ y 3 ( 1 ) = φ ^ y 1 ( 1 ) · ( x 02 ) 2 = 1555 - - - ( 13 )
假设极点xp2、xp3的两个根,取xp2、xp3与KC(s)中的极点相同,即 x p 2 = x p 3 = e - 1 T 2 T s = e - 1 0.22 × 0.1 = 0.635 ;
根据根与系数的关系式得到如下关系式:
- θ e ( 1 ) φ ^ u 2 ( 1 ) 2 = x p 2 - - - ( 14 )
θ e ( 1 ) φ ^ u 3 ( 1 ) = ( x p 2 ) 2 - - - ( 15 )
由式(14)、(15)得到计算
Figure BDA000035690267001126
的表达式如式(16)、(17)所示:
θ e ( 1 ) φ ^ u 2 ( 1 ) = - 2 x p 2 = 2 × e - 1 0.22 × 0.1 = - 1.27 - - - ( 16 )
θ e ( 1 ) φ ^ u 3 ( 1 ) = ( x p 2 ) 2 = e - 1 0.22 × 0.2 = 0.4029 - - - ( 17 )
| K C ( s ) | | s = j ω d = | K M ( z ) | | z = e j ω d T s , 得到如下等式:
θ e ( 1 ) [ φ ^ y 1 ( 1 ) + φ ^ y 2 ( 1 ) + φ ^ y 3 ( 1 ) ] [ 1 + θ e ( 1 ) φ ^ u 2 ( 1 ) + θ e ( 1 ) φ ^ u 3 ( 1 ) ] = K c = 0.064 - - - ( 18 )
式(18)中
Figure BDA000035690267001128
均已知,得到θe(1)的计算表达式如式(19)所示:
θ e ( 1 ) = K c · [ 1 + θ e ( 1 ) φ ^ u 2 ( 1 ) + θ e ( 1 ) φ ^ u 3 ( 1 ) ] [ φ ^ y 1 ( 1 ) + φ ^ y 2 ( 1 ) + φ ^ y 3 ( 1 ) ] = 4.6955 × 10 - 5 - - - ( 19 )
根据式(19)得到的θe(1)的取值计算出
Figure BDA000035690267001120
Figure BDA000035690267001121
由于
Figure BDA000035690267001122
按经验值选取ρk=0.7,γ=106,则
Figure BDA000035690267001129
根据θe(1)的表达式算出φu1(k)=67.3;
步骤4)进一步以步骤3)中得到伪梯度向量初始值
Figure BDA000035690267001124
为初始点,利用被控系统实测的输入输出数据在线更新伪梯度向量的估计值具体实现方式为:
通过极小化伪梯度向量估计准则函数得到,伪梯度向量φ(k)在k时刻的在线更新计算表达式:
φ ^ ( k ) = φ ^ ( k - 1 ) +
[ ψ ( k - 1 ) ψ T ( k - 1 ) + μ 1 O μ L μ O μ L u + L y ] - 1 · ψ ( k - 1 ) [ Δy ( k ) - ψ T ( k - 1 ) φ ^ ( k - 1 ) ] - - - ( 20 )
其中Λμ为对角权重矩阵,按经验值取为 Λ μ = 100 1 1 10 10 10 , 用于限制
Figure BDA00003569026700124
中每个伪梯度向量元素的变化;
步骤5)最后根据在线更新的伪梯度向量的估计值计算发电机广域阻尼自适应控制信号;具体实现方式为:
根据步骤4)中k时刻在线更新的伪梯度向量估计值
Figure BDA00003569026700125
利用发电机广域阻尼自适应控制方法的控制律计算k时刻发电机广域阻尼自适应控制方法的控制信号:
u ( k ) = u ( k - 1 ) + ρ k φ ^ u 1 ( k ) γ + | φ ^ u 1 ( k ) | 2 [ y * ( k + 1 ) - y ( k ) - φ ^ y 1 ( k ) Δy ( k ) - φ ^ y 2 ( k ) Δy ( k - 1 ) - - - ( 21 )
- φ ^ y 3 ( k ) Δy ( k - 2 ) - φ ^ u 2 ( k ) Δu ( k - 1 ) - φ ^ u 3 ( k ) Δu ( k - 2 ) ]
步骤6)将步骤5)中的发电机广域阻尼自适应控制信号作用于被控电力系统的发电机励磁端,k增加1,重复步骤4)~6)直至达到所需计算的发电机广域阻尼控制信号的最大时刻。
在该四机两区系统中的如下故障下检验上述发电机广域阻尼控制的自适应效果:联络线3中点在1s时发生三相0.1秒瞬时短路故障,30s时线路发生断线故障,采用传统发电机广域阻尼控制方法和本发明提出的发电机广域阻尼自适应控制方法的控制效果比较如图3所示。
从图3中的仿真结果可以看出,三相瞬时短路故障下,由于系统结构不变,两种控制方法的控制效果效果基本一致,当系统在30s发生断线故障后,系统结构发生变化,本发明所提出的发电机广域阻尼自适应控制能够进行在线自适应调整,控制效果优于传统发电机广域阻尼控制。

Claims (4)

1.一种基于无模型自适应控制算法的发电机广域阻尼控制方法,其特征在于,该方法包括:
1)根据传统发电机广域阻尼控制方法的传递函数表达式及参数计算公式得到传统发电广域阻尼控制方法的控制参数;
2)定义被控电力系统的动态线性表达式、控制准则函数表达式和伪梯度向量估计准则函数表达式;
3)然后根据传统发电机广域阻尼控制方法得到的控制参数计算发电机广域阻尼自适应控制方法的控制参数初始值,包括伪梯度向量的初始值
Figure FDA00003569026600011
惩罚因子γ;
4)以步骤3)中得到伪梯度向量初始值
Figure FDA00003569026600012
为初始点,利用被控系统实测的输入输出数据在线更新伪梯度向量的估计值
5)最后根据在线更新的伪梯度向量的估计值计算发电机广域阻尼自适应控制信号;
6)将步骤5)中的发电机广域阻尼自适应控制信号作用于被控电力系统的发电机励磁端,k增加1,重复步骤4)~6)直至达到所需计算的发电机广域阻尼控制信号的最大时刻。
2.根据权利要求1所述方法,其特征在于,所述步骤3)根据传统发电机广域阻尼控制方法得到的控制参数计算发电机广域阻尼自适应控制方法的控制参数初始值,包括伪梯度向量的初始值
Figure FDA00003569026600014
具体实现方式为:
取伪阶数Lu=3和Ly=3,通过极小化控制准则函数得到发电机广域阻尼自适应控制方法的控制律表达式为:
u ( k ) = u ( k - 1 ) + ρ k φ ^ u 1 ( k ) γ + | φ ^ u 1 ( k ) | 2 [ y * ( k + 1 ) - y ( k ) - φ ^ y 1 ( k ) Δy ( k ) - φ ^ y 2 ( k ) Δy ( k - 1 ) - - - ( 6 )
- φ ^ y 3 ( k ) Δy ( k - 2 ) - φ ^ u 2 ( k ) Δu ( k - 1 ) - φ ^ u 3 ( k ) Δu ( k - 2 ) ]
其中ρk为步长序列;
将k+1时刻区间频差信号的期望值y*(k+1)取为0;
Figure FDA00003569026600017
所述控制律的离散传递函数表达式为:
K M ( z ) = θ e ( k ) { [ 1 + φ ^ y 1 ( k ) ] z 3 + [ φ ^ y 2 ( k ) - φ ^ y 1 ( k ) ] z 2 + [ φ ^ y 3 ( k ) - φ ^ y 2 ( k ) ] z - φ ^ y 3 ( k ) } z 3 + ( θ e ( k ) φ ^ u 2 ( k ) - 1 ) z 2 + θ e ( k ) ( φ ^ u 3 ( k ) - φ ^ u 2 ( k ) ) z - θ e ( k ) φ ^ u 3 ( k ) - - - ( 7 )
= θ e ( k ) { [ 1 + φ ^ y 1 ( k ) ] z 3 + [ φ ^ y 2 ( k ) - φ ^ y 1 ( k ) ] z 2 + [ φ ^ y 3 ( k ) - φ ^ y 2 ( k ) ] z - φ ^ y 3 ( k ) } ( z - 1 ) [ z 2 + ( θ e ( k ) φ ^ u 2 ( k ) z + θ e ( k ) φ ^ u 3 ( k ) ]
离散传递函数KM(z)由三个零点和三个极点组成,其中一个极点为1;
将k=1时刻KM(z)的零极点取值与KC(s)的零极点相同;其中k=1时刻KM(z)的一个零点配置在极点1附近,另外两个零点和两个极点分别与KC(s)的零点和极点相同;
假设x01,x02,x03为表达式 [ 1 + φ ^ y 1 ( 1 ) ] z 3 + [ φ ^ y 2 ( 1 ) - φ ^ y 1 ( 1 ) ] z 2 + [ φ ^ y 3 ( 1 ) - φ ^ y 2 ( 1 ) ] z - φ ^ y 3 ( 1 ) = 0 的三个根;x01为在极点1附近的零点,(为减小零点x01和极点1对振荡模式的动态性能影响,)取该零点到原点的距离为区间低频振荡模式对应的特征根到虚轴距离的1/5,即其中Ts=0.1s,ξ0=7.576%,ωd=3.662;x02=x03与KC(s)中的零点相同,即 x 02 = x 03 = e - 1 T 1 T s = e - 1 0.34 × 0.1 = 0.7452 ;
令z=1,得到 1 1 + φ ^ y 1 ( 1 ) = ( 1 - x 01 ) ( 1 - x 02 ) 2 = ( 1 - 0.9945 ) ( 1 - 0.7452 ) 2 ;
根据上述等式得到计算的表达式如式(8)所示:
φ ^ y 1 ( 1 ) = 1 ( 1 - x 01 ) ( 1 - x 02 ) 2 - 1 = 2800 - - - ( 8 )
由于
Figure FDA00003569026600027
KM(z)化简为:
≈ θ e ( 1 ) { φ ^ y 1 ( 1 ) z 3 + [ φ ^ y 2 ( 1 ) - φ ^ y 1 ( 1 ) ] z 2 + [ φ ^ y 3 ( 1 ) - φ ^ y 2 ( 1 ) ] z - φ ^ y 3 ( 1 ) } z 3 + ( θ e ( 1 ) φ ^ u 2 ( 1 ) - 1 ) z 2 + θ e ( 1 ) ( φ ^ u 3 ( 1 ) - φ ^ u 2 ( 1 ) ) z - θ e ( 1 ) φ ^ u 3 ( 1 ) - - - ( 9 )
= θ e ( 1 ) [ φ ^ y 1 ( 1 ) z 2 + φ ^ y 2 ( 1 ) z + φ ^ y 3 ( 1 ) ] [ z 2 + θ e ( 1 ) φ ^ u 2 ( 1 ) z + θ e φ ^ u 3 ( 1 ) ]
根据根与系数的关系式得到如下关系式:
- φ ^ y 2 ( 1 ) φ ^ y 1 ( 1 ) = 2 × e - 1 0.34 × 0.1 - - - ( 10 )
φ ^ y 3 ( 1 ) φ ^ y 1 ( 1 ) = ( x 02 ) 2 = e - 1 0.34 × 0.1 × 2 = 0.555 - - - ( 11 )
由式(10)、(11)得到计算
Figure FDA000035690266000217
的表达式如式(12)、(13)所示:
φ ^ y 2 ( 1 ) = - 2 x 02 · φ ^ y 1 ( 1 ) = - 4173 - - - ( 12 )
φ ^ y 3 ( 1 ) = φ ^ y 1 ( 1 ) · ( x 02 ) 2 = 1555 - - - ( 13 )
假设极点xp2、xp3
Figure FDA00003569026600031
的两个根,取xp2、xp3与KC(s)中的极点相同,即 x p 2 = x p 3 = e - 1 T 2 T s = e - 1 0.22 × 0.1 = 0.635 ;
根据根与系数的关系式得到如下关系式:
- θ e ( 1 ) φ ^ u 2 ( 1 ) 2 = x p 2 - - - ( 14 )
θ e ( 1 ) φ ^ u 3 ( 1 ) = ( x p 2 ) 2 - - - ( 15 )
由式(14)、(15)得到计算
Figure FDA000035690266000323
的表达式如式(16)、(17)所示:
θ e ( 1 ) φ ^ u 2 ( 1 ) = - 2 x p 2 = 2 × e - 1 0.22 × 0.1 = - 1.27 - - - ( 16 )
θ e ( 1 ) φ ^ u 3 ( 1 ) = ( x p 2 ) 2 = e - 1 0.22 × 0.2 = 0.4029 - - - ( 17 )
| K C ( s ) | | s = j ω d = | K M ( z ) | | z = e j ω d T s , 得到如下等式:
θ e ( 1 ) [ φ ^ y 1 ( 1 ) + φ ^ y 2 ( 1 ) + φ ^ y 3 ( 1 ) ] [ 1 + θ e ( 1 ) φ ^ u 2 ( 1 ) + θ e ( 1 ) φ ^ u 3 ( 1 ) ] = K c = 0.064 - - - ( 18 )
式(18)中
Figure FDA000035690266000324
Figure FDA000035690266000315
均已知,得到θe(1)的计算表达式如式(19)所示:
θ e ( 1 ) = K c · [ 1 + θ e ( 1 ) φ ^ u 2 ( 1 ) + θ e ( 1 ) φ ^ u 3 ( 1 ) ] [ φ ^ y 1 ( 1 ) + φ ^ y 2 ( 1 ) + φ ^ y 3 ( 1 ) ] = 4.6955 × 10 - 5 - - - ( 19 )
根据式(19)得到的θe(1)的取值计算出
Figure FDA000035690266000317
Figure FDA000035690266000318
由于按经验值选取ρk=0.7,γ=106,则
Figure FDA000035690266000320
根据θe(1)的表达式算出φu1(k)=67.3。
3.根据权利要求2所述方法,其特征在于,所述步骤4)进一步以步骤3)中得到伪梯度向量初始值
Figure FDA000035690266000321
为初始点,利用被控系统实测的输入输出数据在线更新伪梯度向量的估计值
Figure FDA000035690266000322
具体实现方式为:
通过极小化伪梯度向量估计准则函数得到,伪梯度向量φ(k)在k时刻的在线更新计算表达式:
φ ( k ) = φ ^ ( k - 1 ) +
[ ψ ( k - 1 ) ψ T ( k - 1 ) + μ 1 O μ L u O μ L u + L y ] - 1 · ψ ( k - 1 ) [ Δy ( k ) - ψ T ( k - 1 ) φ ^ ( k - 1 ) ] - - - ( 20 )
其中Λμ为对角权重矩阵,按经验值取为 Λ μ = 100 1 1 10 10 10 , 用于限制中每个伪梯度向量元素的变化。
4.根据权利要求3所述方法,其特征在于,所述步骤5)最后根据在线更新的伪梯度向量的估计值计算发电机广域阻尼自适应控制信号;具体实现方式为:
根据步骤4)中k时刻在线更新的伪梯度向量估计值
Figure FDA00003569026600045
利用发电机广域阻尼自适应控制方法的控制律计算k时刻发电机广域阻尼自适应控制方法的控制信号为:
u ( k ) = u ( k - 1 ) + ρ k φ ^ u 1 ( k ) γ + | φ ^ u 1 ( k ) | 2 [ y * ( k + 1 ) - y ( k ) - φ ^ y 1 ( k ) Δy ( k ) - φ ^ y 2 ( k ) Δy ( k - 1 ) - - - ( 21 ) .
- φ ^ y 3 ( k ) Δy ( k - 2 ) - φ ^ u 2 ( k ) Δu ( k - 1 ) - φ ^ u 3 ( k ) Δu ( k - 2 ) ]
CN201310316663.XA 2013-07-25 2013-07-25 基于无模型自适应控制算法的发电机广域阻尼控制方法 Expired - Fee Related CN103595062B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310316663.XA CN103595062B (zh) 2013-07-25 2013-07-25 基于无模型自适应控制算法的发电机广域阻尼控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310316663.XA CN103595062B (zh) 2013-07-25 2013-07-25 基于无模型自适应控制算法的发电机广域阻尼控制方法

Publications (2)

Publication Number Publication Date
CN103595062A true CN103595062A (zh) 2014-02-19
CN103595062B CN103595062B (zh) 2016-05-25

Family

ID=50085097

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310316663.XA Expired - Fee Related CN103595062B (zh) 2013-07-25 2013-07-25 基于无模型自适应控制算法的发电机广域阻尼控制方法

Country Status (1)

Country Link
CN (1) CN103595062B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103997048A (zh) * 2014-05-28 2014-08-20 广东电网公司电力科学研究院 用于抑制低频振荡的电网调控方法和系统
CN104836504A (zh) * 2015-05-15 2015-08-12 合肥工业大学 凸极式永磁同步电机精确转矩输出的自适应容错控制方法
CN105388759A (zh) * 2015-12-04 2016-03-09 广西大学 煮糖结晶过程多目标优化与数据驱动控制方法
CN105487385A (zh) * 2016-02-01 2016-04-13 金陵科技学院 基于无模型自适应内模控制方法
CN107196322A (zh) * 2017-07-19 2017-09-22 国网天津市电力公司电力科学研究院 一种基于无模型自适应预测控制的低频振荡广域阻尼控制方法
CN108011380A (zh) * 2017-11-28 2018-05-08 昆明理工大学 一种基于无模型自适应的附加阻尼控制方法
CN109814388A (zh) * 2019-02-01 2019-05-28 浙江大学 参数自整定的miso异因子偏格式无模型控制方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101202451A (zh) * 2007-12-13 2008-06-18 南方电网技术研究中心 一种电力系统广域阻尼的控制系统及其方法
EP2325968A1 (en) * 2009-11-18 2011-05-25 ABB Research Ltd. Tuning a power oscillation damping unit

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101202451A (zh) * 2007-12-13 2008-06-18 南方电网技术研究中心 一种电力系统广域阻尼的控制系统及其方法
EP2325968A1 (en) * 2009-11-18 2011-05-25 ABB Research Ltd. Tuning a power oscillation damping unit

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
赵艺等: "基于辨识和留数的发电机广域附加阻尼控制器设计", 《电力系统保护与控制》, vol. 40, no. 19, 1 October 2012 (2012-10-01), pages 1 - 6 *
赵艺等: "多输入多输出无模型自适应广域阻尼控制器设计", 《清华大学学报(自然科学版)》, vol. 53, no. 4, 30 April 2013 (2013-04-30), pages 453 - 548 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103997048A (zh) * 2014-05-28 2014-08-20 广东电网公司电力科学研究院 用于抑制低频振荡的电网调控方法和系统
CN103997048B (zh) * 2014-05-28 2016-06-08 广东电网公司电力科学研究院 用于抑制低频振荡的电网调控方法和系统
CN104836504A (zh) * 2015-05-15 2015-08-12 合肥工业大学 凸极式永磁同步电机精确转矩输出的自适应容错控制方法
CN104836504B (zh) * 2015-05-15 2017-06-06 合肥工业大学 凸极式永磁同步电机精确转矩输出的自适应容错控制方法
CN105388759A (zh) * 2015-12-04 2016-03-09 广西大学 煮糖结晶过程多目标优化与数据驱动控制方法
CN105388759B (zh) * 2015-12-04 2018-04-20 广西大学 煮糖结晶过程多目标优化与数据驱动控制方法
CN105487385A (zh) * 2016-02-01 2016-04-13 金陵科技学院 基于无模型自适应内模控制方法
CN105487385B (zh) * 2016-02-01 2019-02-15 金陵科技学院 基于无模型自适应内模控制方法
CN107196322A (zh) * 2017-07-19 2017-09-22 国网天津市电力公司电力科学研究院 一种基于无模型自适应预测控制的低频振荡广域阻尼控制方法
CN108011380A (zh) * 2017-11-28 2018-05-08 昆明理工大学 一种基于无模型自适应的附加阻尼控制方法
CN108011380B (zh) * 2017-11-28 2021-05-14 昆明理工大学 一种基于无模型自适应的附加阻尼控制方法
CN109814388A (zh) * 2019-02-01 2019-05-28 浙江大学 参数自整定的miso异因子偏格式无模型控制方法

Also Published As

Publication number Publication date
CN103595062B (zh) 2016-05-25

Similar Documents

Publication Publication Date Title
CN103595062A (zh) 基于无模型自适应控制算法的发电机广域阻尼控制方法
CN103558436B (zh) 基于单相锁相环算法的检测电网电压幅值、频率和相角的方法
CN102545245B (zh) 基于端口供给能量的电力系统振荡源定位方法
CN103296688B (zh) 一种电力系统稳定器的实现方法
CN107800343A (zh) 异步电机自抗扰控制器的设计方法
CN103259472B (zh) 大型发电机组次同步振荡的反振荡自适应抑制系统
CN102510060B (zh) 一种电力系统频率特性系数的计算方法
Malla et al. Online learning control for harmonics reduction based on current controlled voltage source power inverters
CN103336882A (zh) 一种基于时域仿真的全过程动态电压稳定裕度评估方法
Qi et al. Tracking-differentiator-based dynamic virtual inertial control of offshore wind power plant for frequency regulation
CN104375876A (zh) 一种输入量突变情况下的0+误差免疫电磁暂态仿真算法
Frutos et al. Low-frequency oscillations analysis in AC railway networks using eigenmode identification
Mao et al. Improving power system dynamic performance using wide-area high-voltage direct current damping control
CN107959294A (zh) 一种基于线性化潮流的配电网灵敏度计算方法
CN102147439B (zh) 基于低频振荡阻尼灵敏度和等效交流增益的pss时间常数优化试验方法
CN101408908A (zh) 基于优化的电力系统实用时滞稳定裕度计算方法
CN111740449A (zh) 基于附加励磁信号注入法的风电机群对同步发电机组阻尼特性分析方法
CN103490435A (zh) 一种利用永磁同步风机旋转动能提高电力系统阻尼的方法
CN108808703B (zh) 基于低阶igd-irk的时滞电力系统小干扰稳定性分析方法
Wang et al. Mitigation of subsynchronous control interaction in DFIG-based wind farm using repetitive-pi
CN106026138A (zh) 一种即插即用式的电力系统稳定器设计算法
Zhu et al. Adaptive and coordinated oscillation damping control using measurement-driven approach
CN116436091A (zh) 构网型变流器暂态稳定提升控制方法和构网型变流器系统
CN110034562A (zh) 一种静止同步补偿器与发电机励磁鲁棒协调的控制方法
Farraj et al. On using distributed energy resources to reshape the dynamics of power systems during transients

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160525

Termination date: 20170725

CF01 Termination of patent right due to non-payment of annual fee