CN103199546A - 动态无功补偿装置最优二次调节器及其设计方法 - Google Patents
动态无功补偿装置最优二次调节器及其设计方法 Download PDFInfo
- Publication number
- CN103199546A CN103199546A CN2013101061277A CN201310106127A CN103199546A CN 103199546 A CN103199546 A CN 103199546A CN 2013101061277 A CN2013101061277 A CN 2013101061277A CN 201310106127 A CN201310106127 A CN 201310106127A CN 103199546 A CN103199546 A CN 103199546A
- Authority
- CN
- China
- Prior art keywords
- reactive power
- power compensator
- secondary regulator
- compensation device
- matrix
- 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
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/30—Reactive power compensation
Landscapes
- Control Of Electrical Variables (AREA)
Abstract
本发明提供一种动态无功补偿装置最优二次调节器设计方法,包括以下步骤:采样无功补偿装置输出的无功功率和晶闸管放大系数,利用模型辨识算法获得无功补偿装置传递函数;按照控制系统理论能观测标准型的形式,设计无功补偿装置连续时间状态方程;建立优化指标函数;求解黎卡提方程;计算无功补偿装置最优二次调节器增益矩阵;获得无功补偿装置控制系统输入晶闸管放大系数。利用本发明方法可以按照设计者要求的动静态性能,通过求解黎卡提方程即可获得最优的控制器参数,保证系统具有最优控制性能;并且还能利用最优二次调节器的作用实现无功功率输出跟踪控制,具有非常优异的系统性能。
Description
技术领域
本发明涉及电力供电系统中动态无功补偿装置控制系统的设计方法,具体为动态无功补偿装置最优二次调节器及其设计方法。
技术背景
在钢铁、冶金等工业电力系统中,由于负荷的变化而需要对系统进行动态无功补偿,目前常用的控制方法是利用PID控制算法对系统进行闭环控制,由于负荷的变化导致控制系统在PID参数整定方面较为复杂。工程实际中应用的无功补偿装置通常采用功率因数控制、电压控制、无功功率控制等几种形式的控制策略,这些控制策略有一个共同点就是通过调整晶闸管控制角来调整无功补偿装置输出的无功功率,通过分析可以发现无功补偿装置在母线电压一定的情况下,输出的无功功率与晶闸管控制角之间对应的函数关系保持不变。为了在工程应用中减少现场的调试,可以将无功功率与晶闸管之间的对应关系设计成一个稳定的闭环控制系统,当采用电压、功率因数等不同的控制方式时,只需要调整外环控制器参数即可实现无功补偿装置的稳定运行。
发明内容
本发明要解决的技术问题是:提供一种动态无功补偿装置最优二次调节器及其设计方法,保证系统具有最优控制性能。
本发明为解决上述技术问题所采取的技术方案为:一种动态无功补偿装置最优二次调节器设计方法,其特征在于:它包括以下步骤:
S1、采样无功补偿装置输出的无功功率Q和晶闸管放大系数u,选择无功补偿装置模型阶次n,利用模型辨识算法获得无功补偿装置传递函数;
S2、按照控制系统理论能观测标准型的形式,根据步骤S1获得的无功补偿装置传递函数来设计无功补偿装置连续时间状态方程系数矩阵A、B、C,则无功补偿装置连续时间状态方程表示为Q(t)=Cx(t),其中x(t)为t时刻测量的系统状态变量值,表示x(t)的微分,u(t)为t时刻晶闸管放大系数值,Q(t)为t时刻无功补偿装置输出的无功功率;
S3、建立优化指标函数其中N为系统状态能量的偏重矩阵,R为系统控制能量的偏重矩阵,且N为n×n正定对称矩阵,R为大于零的实数,x为系统状态变量,上标“T”表示矩阵的转置运算;所述的N和R均按照经验和实际的动静态需要设定;
S4、求解黎卡提方程ATP+PA-PBR-1BTP+N=0,获得正定对称矩阵P;
S5、计算无功补偿装置最优二次调节器增益矩阵K=R-1BTP;
S6、令无功补偿装置给定无功功率为r,则无功补偿装置控制系统输入晶闸管放大系数u为:u=K(1)r-Kx(t),其中K(1)为增益矩阵K的第一个元素。
按上述方案,还包括步骤S7、对N和R进行不断的尝试取值,通过比较不同取值下获得的无功功率跟踪响应曲线,来确定N和R的优选值。
一种基于上述动态无功补偿装置最优二次调节器设计方法的动态无功补偿装置最优二次调节器,其特征在于:它包括:
无功补偿装置模型,用于采样无功补偿装置输出的无功功率Q和晶闸管放大系数u,选择无功补偿装置模型阶次n,利用模型辨识算法获得无功补偿装置传递函数;
无功补偿装置连续时间模型,用于按照控制系统理论能观测标准型的形式,根据无功补偿装置传递函数来设计无功补偿装置连续时间状态方程系数矩阵A、B、C,则无功补偿装置连续时间状态方程表示为Q(t)=Cx(t),其中x(t)为t时刻测量的系统状态变量值,表示x(t)的微分,u(t)为t时刻晶闸管放大系数值,Q(t)为t时刻无功补偿装置输出的无功功率;
优化指标装置,用于建立优化指标函数其中N为系统状态能量的偏重矩阵,R为系统控制能量的偏重矩阵,且N为n×n正定对称矩阵,R为大于零的实数,x为系统状态变量,上标“T”表示矩阵的转置运算;所述的N和R均按照经验和实际的动静态需要设定;
无功补偿装置最优二次调节器增益装置,用于计算无功补偿装置最优二次调节器增益矩阵K=R-1BTP,并对系统状态变量值进行增益计算,得到增益系统状态变量值;其中P为正定对称矩阵,由求解黎卡提方程ATP+PA-PBR-1BTP+N=0获得;
无功功率增益装置,用于对无功补偿装置给定无功功率r进行增益计算,获得增益无功功率,其增益系数为无功补偿装置最优二次调节器增益矩阵K的第一个元素;
减法器,用于将增益无功功率减去增益系统状态变量值,其结果作为无功补偿装置控制系统输入晶闸管放大系数u,即:u=K(1)r-Kx(t)。
本发明的有益效果为:利用本发明方法可以按照设计者要求的动静态性能,通过求解黎卡提方程即可获得最优的控制器参数,保证系统具有最优控制性能;并且还能利用最优二次调节器的作用实现无功功率输出跟踪控制,具有非常优异的系统性能。
附图说明
图1为本发明一实施例的原理图。
图2为调节器1获得的系统无功跟踪曲线。
图3为调节器2与调节器1获得的系统无功跟踪曲线对比图。
具体实施方式
图1为本发明一实施例的原理图,动态无功补偿装置最优二次调节器包括:无功补偿装置模型,用于采样无功补偿装置输出的无功功率Q和晶闸管放大系数u,选择无功补偿装置模型阶次n,利用模型辨识算法获得无功补偿装置传递函数;无功补偿装置连续时间模型,用于按照控制系统理论能观测标准型的形式,根据无功补偿装置传递函数来设计无功补偿装置连续时间状态方程系数矩阵A、B、C,则无功补偿装置连续时间状态方程表示为Q(t)=Cx(t),其中x(t)为t时刻测量的系统状态变量值,表示x(t)的微分,u(t)为t时刻晶闸管放大系数值,Q(t)为t时刻无功补偿装置输出的无功功率;优化指标装置,用于建立优化指标函数其中N为系统状态能量的偏重矩阵,R为系统控制能量的偏重矩阵,且N为n×n正定对称矩阵,R为大于零的实数,x为系统状态变量,上标“T”表示矩阵的转置运算;所述的N和R均按照经验和实际的动静态需要设定;无功补偿装置最优二次调节器增益装置,用于计算无功补偿装置最优二次调节器增益矩阵K=R-1BTP,并对系统状态变量值进行增益计算,得到增益系统状态变量值;其中P为正定对称矩阵,由求解黎卡提方程ATP+PA-PBR-1BTP+N=0获得;无功功率增益装置,用于对无功补偿装置给定无功功率r进行增益计算,获得增益无功功率,其增益系数为无功补偿装置最优二次调节器增益矩阵K的第一个元素;减法器,用于将增益无功功率减去增益系统状态变量值,其结果作为无功补偿装置控制系统输入晶闸管放大系数u,即:u=K(1)r-Kx(t)。
上述动态无功补偿装置最优二次调节器的设计方法包括以下步骤:S1、采样无功补偿装置输出的无功功率Q和晶闸管放大系数u,选择无功补偿装置模型阶次n,利用模型辨识算法获得无功补偿装置传递函数;S2、按照控制系统理论能观测标准型的形式,根据步骤S1获得的无功补偿装置传递函数来设计无功补偿装置连续时间状态方程系数矩阵A、B、C,则无功补偿装置连续时间状态方程表示为Q(t)=Cx(t),其中x(t)为t时刻测量的系统状态变量值,表示x(t)的微分,u(t)为t时刻晶闸管放大系数值,Q(t)为t时刻无功补偿装置输出的无功功率;S3、建立优化指标函数其中N为系统状态能量的偏重矩阵,R为系统控制能量的偏重矩阵,且N为n×n正定对称矩阵,R为大于零的实数,x为系统状态变量,上标“T”表示矩阵的转置运算;所述的N和R均按照经验和实际的动静态需要设定;S4、求解黎卡提方程ATP+PA-PBR-1BTP+N=0,获得正定对称矩阵P;S5、计算无功补偿装置最优二次调节器增益矩阵K=R-1BTP;S6、令无功补偿装置给定无功功率为r,则无功补偿装置控制系统输入晶闸管放大系数u为:u=K(1)r-Kx(t),其中K(1)为增益矩阵K的第一个元素。
本方法还可以对N和R进行不断的尝试取值,通过比较不同取值下获得的无功功率跟踪响应曲线,来确定N和R的优选值。
某钢厂6.5kV母线上接有一台TCR型无功补偿装置,按照步骤S1通过开环模型辨识方法实际采样无功功率Q和晶闸管放大系数u,选择系统阶次n=3,通过模型辨识算法获得无功补偿装置连续时间模型为
针对公式(1)中的无功补偿装置模型,按照步骤S2获得无功补偿装置状态方程如下:
获得如公式(2)和(3)所示的无功补偿装置连续时间状态方程后,按照个人经验和实际的动静态需要可以确定矩阵N和R,本实施例中N应为3×3的正定对称矩阵,R为正实数。
首先选择 R=1,按照步骤S4求解黎卡提方程可得:
按照步骤S5可得最优二次调节器增益矩阵:
此时按照步骤S6和附图1所示的控制结构图构建无功补偿装置控制系统,令此系统为调节器1,若无功功率给定值为10000乏,基于上述设计得到的无功功率跟踪响应曲线如附图2所示。从附图2中可以发现状态1跟踪给定无功功率变化,而且状态1响应曲线没有超调,动静态性能完全按照设计者给定的矩阵N和R进行最优性能设计。
为了体现本发明设计的无功补偿装置最优二次调节器能够根据设计者的要求获得最优的动静态性能,可令矩阵 R=1对系统增益矩阵进行再次设计,令此系统为调节器2,矩阵N的前后差距在于状态1的系数发生了变化,调节器2在设计中更加注重状态1的变化,获得最优二次调节器将使状态1具有更快的动态性能。按照上述步骤可得新的增益矩阵为
将调节器2的增益矩阵代入控制系统获得的无功功率跟踪曲线(即图中曲线2)与调节器1获得的无功功率跟踪曲线(即图中曲线1)如附图3所示,对比前后两次无功功率跟踪曲线可以发现,第二次系统动态响应明显快于第一次,充分体现了矩阵N对于系统性能的影响。
Claims (3)
1.一种动态无功补偿装置最优二次调节器设计方法,其特征在于:它包括以下步骤:
S1、采样无功补偿装置输出的无功功率Q和晶闸管放大系数u,选择无功补偿装置模型阶次n,利用模型辨识算法获得无功补偿装置传递函数;
S2、按照控制系统理论能观测标准型的形式,根据步骤S1获得的无功补偿装置传递函数来设计无功补偿装置连续时间状态方程系数矩阵A、B、C,则无功补偿装置连续时间状态方程表示为Q(t)=Cx(t),其中x(t)为t时刻测量的系统状态变量值,表示x(t)的微分,u(t)为t时刻晶闸管放大系数值,Q(t)为t时刻无功补偿装置输出的无功功率;
S3、建立优化指标函数其中N为系统状态能量的偏重矩阵,R为系统控制能量的偏重矩阵,且N为n×n正定对称矩阵,R为大于零的实数,x为系统状态变量,上标“T”表示矩阵的转置运算;所述的N和R均按照经验和实际的动静态需要设定;
S4、求解黎卡提方程ATP+PA-PBR-1BTP+N=0,获得正定对称矩阵P;
S5、计算无功补偿装置最优二次调节器增益矩阵K=R-1BTP;
S6、令无功补偿装置给定无功功率为r,则无功补偿装置控制系统输入晶闸管放大系数u为:u=K(1)r-Kx(t),其中K(1)为增益矩阵K的第一个元素。
2.根据权利要求1所述的动态无功补偿装置最优二次调节器设计方法,其特征在于:还包括步骤S7、对N和R进行不断的尝试取值,通过比较不同取值下获得的无功功率跟踪响应曲线,来确定N和R的优选值。
3.一种基于权利要求1或2所述的动态无功补偿装置最优二次调节器设计方法的动态无功补偿装置最优二次调节器,其特征在于:它包括:
无功补偿装置模型,用于采样无功补偿装置输出的无功功率Q和晶闸管放大系数u,选择无功补偿装置模型阶次n,利用模型辨识算法获得无功补偿装置传递函数;
无功补偿装置连续时间模型,用于按照控制系统理论能观测标准型的形式,根据无功补偿装置传递函数来设计无功补偿装置连续时间状态方程系数矩阵A、B、C,则无功补偿装置连续时间状态方程表示为Q(t)=Cx(t),其中x(t)为t时刻测量的系统状态变量值,表示x(t)的微分,u(t)为t时刻晶闸管放大系数值,Q(t)为t时刻无功补偿装置输出的无功功率;
优化指标装置,用于建立优化指标函数其中N为系统状态能量的偏重矩阵,R为系统控制能量的偏重矩阵,且N为n×n正定对称矩阵,R为大于零的实数,x为系统状态变量,上标“T”表示矩阵的转置运算;所述的N和R均按照经验和实际的动静态需要设定;
无功补偿装置最优二次调节器增益装置,用于计算无功补偿装置最优二次调节器增益矩阵K=R-1BTP,并对系统状态变量值进行增益计算,得到增益系统状态变量值;其中P为正定对称矩阵,由求解黎卡提方程ATP+PA-PBR-1BTP+N=0获得;
无功功率增益装置,用于对无功补偿装置给定无功功率r进行增益计算,获得增益无功功率,其增益系数为无功补偿装置最优二次调节器增益矩阵K的第一个元素;
减法器,用于将增益无功功率减去增益系统状态变量值,其结果作为无功补偿装置控制系统输入晶闸管放大系数u,即:u=K(1)r-Kx(t)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310106127.7A CN103199546B (zh) | 2013-03-29 | 2013-03-29 | 动态无功补偿装置最优二次调节器及其设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310106127.7A CN103199546B (zh) | 2013-03-29 | 2013-03-29 | 动态无功补偿装置最优二次调节器及其设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103199546A true CN103199546A (zh) | 2013-07-10 |
CN103199546B CN103199546B (zh) | 2015-01-28 |
Family
ID=48721950
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310106127.7A Expired - Fee Related CN103199546B (zh) | 2013-03-29 | 2013-03-29 | 动态无功补偿装置最优二次调节器及其设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103199546B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104009483A (zh) * | 2014-05-28 | 2014-08-27 | 中冶南方工程技术有限公司 | 电抗放大系数倒数与无功功率关系的获取方法及装置 |
CN105846431A (zh) * | 2016-06-06 | 2016-08-10 | 国网福建省电力有限公司 | 一种电网谐波电流信号跟踪控制方法 |
CN105958493A (zh) * | 2016-06-06 | 2016-09-21 | 国网福建省电力有限公司 | 电网谐波电流信号跟踪控制方法 |
CN108008630A (zh) * | 2016-11-01 | 2018-05-08 | 贵州火星探索科技有限公司 | 一种无人机位置环的补偿器设计方法及装置 |
CN108733091A (zh) * | 2018-05-07 | 2018-11-02 | 安徽工程大学 | 一种数字阀门流量控制系统及其最优控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030187623A1 (en) * | 2002-03-06 | 2003-10-02 | Bayard David S. | High accuracy inertial sensors from inexpensive components |
CN101630842A (zh) * | 2009-08-12 | 2010-01-20 | 江苏大学 | 有源电力滤波器的逆系统控制器的实现方法 |
US20120277888A1 (en) * | 2011-04-29 | 2012-11-01 | Georgia Tech Research Corporation | Systems and methods for parameter dependent riccati equation approaches to adaptive control |
-
2013
- 2013-03-29 CN CN201310106127.7A patent/CN103199546B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030187623A1 (en) * | 2002-03-06 | 2003-10-02 | Bayard David S. | High accuracy inertial sensors from inexpensive components |
CN101630842A (zh) * | 2009-08-12 | 2010-01-20 | 江苏大学 | 有源电力滤波器的逆系统控制器的实现方法 |
US20120277888A1 (en) * | 2011-04-29 | 2012-11-01 | Georgia Tech Research Corporation | Systems and methods for parameter dependent riccati equation approaches to adaptive control |
Non-Patent Citations (1)
Title |
---|
侯祥林等: "基于动态设计变量优化方法的代数黎卡提方程算法与应用", 《沈阳建筑大学学报(自然科学版)》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104009483A (zh) * | 2014-05-28 | 2014-08-27 | 中冶南方工程技术有限公司 | 电抗放大系数倒数与无功功率关系的获取方法及装置 |
CN105846431A (zh) * | 2016-06-06 | 2016-08-10 | 国网福建省电力有限公司 | 一种电网谐波电流信号跟踪控制方法 |
CN105958493A (zh) * | 2016-06-06 | 2016-09-21 | 国网福建省电力有限公司 | 电网谐波电流信号跟踪控制方法 |
CN108008630A (zh) * | 2016-11-01 | 2018-05-08 | 贵州火星探索科技有限公司 | 一种无人机位置环的补偿器设计方法及装置 |
CN108008630B (zh) * | 2016-11-01 | 2022-04-19 | 贵州火星探索科技有限公司 | 一种无人机位置环的补偿器设计方法及装置 |
CN108733091A (zh) * | 2018-05-07 | 2018-11-02 | 安徽工程大学 | 一种数字阀门流量控制系统及其最优控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103199546B (zh) | 2015-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103199546B (zh) | 动态无功补偿装置最优二次调节器及其设计方法 | |
CN109212974A (zh) | 区间时变时滞系统的鲁棒模糊预测容错控制方法 | |
Xu et al. | Application of multi-model switching predictive functional control on the temperature system of an electric heating furnace | |
CN103401256A (zh) | 火力发电厂大型汽轮机功频电液控制系统参数的优化方法 | |
CN103660991A (zh) | 一种常导磁浮列车非线性悬浮控制方法及其控制系统 | |
CN104950254A (zh) | 燃气蒸汽联合循环机组的一次调频试验方法、装置及系统 | |
CN103760768A (zh) | 基于正常运行数据的发电机组多变量系统辨识方法 | |
CN105470953A (zh) | 一种水轮机调节系统用电子调节器模型 | |
CN104932256B (zh) | 基于优化迭代算法的时滞广域电力系统控制器 | |
Belyakov et al. | Linear quadratic Gaussian controller design for plasma current, position and shape control system in ITER | |
CN103454921B (zh) | 飞行控制系统非线性跟踪控制器设计的正切线性化方法 | |
Zhao et al. | Design ADRC for two special kinds of systems | |
Lerm et al. | Avoiding Hopf bifurcations in power systems via set-points tuning | |
CN102820650A (zh) | 基于干扰估计的电力系统负荷频率控制器设计方法 | |
Takarics et al. | Gain scheduling for nonlinear systems via integral quadratic constraints | |
CN102866629B (zh) | 随机系统基于动静混合神经网络建模的抗干扰控制方法 | |
CN103199545B (zh) | 动态无功补偿装置最优二次高斯控制器及其设计方法 | |
CN103401006B (zh) | 一种固体氧化物燃料电池电压多模型融合控制方法 | |
CN106571624A (zh) | 考虑电网n-1约束的最优潮流方法 | |
CN102722101A (zh) | 基于闭环频域的辨识方法和系统 | |
EP2858199B1 (en) | Voltage control device, control method thereof and voltage control program | |
CN102513369B (zh) | 液压伺服模拟控制方法 | |
CN106451526B (zh) | 孤网运行时汽轮机opc保护设置方法 | |
Wan et al. | Practical nonlinear excitation control for a single-machine infinite-bus power system based on a detailed model | |
Zemtsov et al. | Using the robust PID controller to manage the population of thermostatically controlled loads |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150128 Termination date: 20160329 |