CN108427397A - 一种水轮机调速系统控制参数的自动整定方法 - Google Patents
一种水轮机调速系统控制参数的自动整定方法 Download PDFInfo
- Publication number
- CN108427397A CN108427397A CN201810211577.5A CN201810211577A CN108427397A CN 108427397 A CN108427397 A CN 108427397A CN 201810211577 A CN201810211577 A CN 201810211577A CN 108427397 A CN108427397 A CN 108427397A
- Authority
- CN
- China
- Prior art keywords
- control parameter
- turbine governor
- governor system
- optimal
- fuzzy
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000005457 optimization Methods 0.000 claims abstract description 17
- 239000002245 particle Substances 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 230000001133 acceleration Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 230000003362 replicative effect Effects 0.000 claims description 3
- 241000406668 Loxodonta cyclotis Species 0.000 claims 1
- 239000000686 essence Substances 0.000 claims 1
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0243—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/20—Pc systems
- G05B2219/24—Pc safety
- G05B2219/24065—Real time diagnostics
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种水轮机调速系统控制参数的自动整定方法,包含以下步骤:(1)采集实际机组水轮机调速系统动态过程数据;(2)建立控制对象的T‑S模糊模型;(3)建立基于T‑S模糊模型的仿真系统;(4)建立PID控制参数优化目标函数;(5)运用GSA算法求解PID控制参数优化目标函数,得到最优控制参数。本发明采用T‑S模糊模型辨识得到实际机组水轮机调速系统控制对象的仿真系统,在该仿真系统基础上进行控制参数优化,优化后的控制参数能直接应用于水轮机调速系统,可以得到最优控制参数向量。
Description
技术领域
本发明涉及一种水轮机控制技术,具体是一种水轮机调速系统控制参数的自动整定方法。
背景技术
水轮机调节系统是水电机组的核心控制系统,承担着稳定机组频率和调节机组功率的重任,本发明采用T-S模糊模型辨识得到实际机组水轮机调速系统控制对象的仿真系统,在该仿真系统基础上进行控制参数优化,优化后的控制参数能直接应用于水轮机调速系统。
发明内容
本发明的目的在于提供一种水轮机调速系统控制参数的自动整定方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:
一种水轮机调速系统控制参数的自动整定方法,包含以下步骤:
(1)采集实际机组水轮机调速系统动态过程数据;
(2)建立控制对象的T-S模糊模型;
(3)建立基于T-S模糊模型的仿真系统;
(4)建立PID控制参数优化目标函数;
(5)运用GSA算法求解PID控制参数优化目标函数,得到最优控制参数,fITAE(Kp,Ki,Kd)是PID控制参数的函数,当PID控制参数发生变化时,利用步骤(3)中建立的仿真系统得到一组机组频率响应曲线,计算对应的fITAE值,极小化目标函数,得到最优控制参数。
作为本发明再进一步的方案:所述步骤(1)具体是:搜集不同PID控制参数、不同频率扰动值所对应的动态过程数据样本(ck,uk,xk)j组成数据样本集,ck为k时刻的频 率扰动值,uk为k时刻调速器输出的控制量序列,xk为k时刻的机组频率,k=1,..,Ns,Ns为动态过程最大采样点数;j=1,…,Nc,表示不同的动态过程,总共采集Ns*Nc个样本组成样本集,将样本集中若干动态过程样本作为训练样本;其中所述水轮机调速系统包括PID控制器、接力器、引水系统、水轮机,其中,接力器、引水系统、水轮机组成控制对象,由PID控制器进行调节控制。
作为本发明再进一步的方案:所述步骤(2)具体是:提取训练样本,以(ck,uk)为输入数据,xk为输出数据,通过训练样本集建立控制对象的T-S模糊模型;T-S模糊模型训练完成后,通过测试样本测试T-S模糊模型精度,将测试样本中的(ck,uk)作为模糊模型输入数据,计算T-S模糊模型输出tk与测试样本xk的均方根误差,均方根误差越小,说明T-S模糊模型与实际机组水轮机调速系统中的控制对象的误差越小,从而选择最优T-S模糊模型。
作为本发明再进一步的方案:所述步骤(3)具体是:用得到的控制对象的T-S模糊模型替换控制对象,得到基于T-S模糊模型的水轮机调速系统的仿真系统,其中PID控制器按照实际机组水轮机调速系统PID控制器进行设计,Kp、Ki和Kd分别为比例、积分和微分增益,是需要整定的PID控制参数。
作为本发明再进一步的方案:所述步骤(4)具体是:采用离散形式ITAE指标作为控制参数优化的目标函数,目标函数定义为:其中c(k)为频率扰动值,为一常数,x为机组频率,Ns为采样点数,T为时间序列。
作为本发明再进一步的方案:所述步骤(5)具体包含以下步骤:
A、算法初始化:设置算法参数,包括群体规模Np、最大迭代次数Nmax_it、引力常数初值G0、衰减指数α;确定优化控制参数范围,在此区间初始化群体中粒子的位置向量分别代表控 制参数Kp、Ki和Kd;初始化历史最优目标函数值fbest和最优解Xbest,设置迭代次数t=0;
B、计算目标函数值:从位置向量中取出控制参数,将控制参数输入步骤(3)中建立的基于T-S模糊模型的水轮机调速系统的仿真系统,得到仿真机组频率,按照步骤(4)计算粒子目标函数值fi,i=1,…,Np;
C、求出群体最优目标函数值fbest_current=min({fi,i=1,…,Np})和对应的粒子位置Xbest_current,如果fbest_current<fbest,则fbest<fbest_current,Xbest=Xbest_current;
D、计算粒子质量Mi,i=1,…,Np:
E、其中,fworst_current=max({fi,i=1,…,Np});
计算引力Fi d:
其中,randj为[0,1]范围随机数,d=1,…,3;
F、计算加速度:
G、更新粒子速度和粒子位置
其中,randi为[0,1]范围随机数;
H、t=t+1;如果t<Nmax_it,转至Step 2,否则结束循环,得到最优目标函数值fbest和最优位置向量Xbest,所述最优位置向量即为最优控制参数向量
与现有技术相比,本发明的有益效果是:本发明采用T-S模糊模型辨识得到实际机组水轮机调速系统控制对象的仿真系统,在该仿真系统基础上进行控制参数优化,优化后的控制参数能直接应用于水轮机调速系统,可以得到最优控制参数向量。
具体实施方式
下面对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例中,一种水轮机调速系统控制参数的自动整定方法,包含以下步骤:
步骤(1):采集实际机组水轮机调速系统动态过程数据:搜集不同PID控制参数、不同频率扰动值所对应的动态过程数据样本(ck,uk,xk)j组成数据样本集,ck为k时刻的频率扰动值,uk为k时刻调速器输出的控制量序列,xk为k时刻的机组频率,k=1,..,Ns,Ns为动态过程最大采样点数;j=1,…,Nc,表示不同的动态过程,总共采集Ns*Nc个样本组成样本集,将样本集中若干动态过程样本作为训练样本;其中所述水轮机调速系统包括PID控制器、接力器、引水系统、水轮机,其中,接力器、引水系统、水轮机组成控制对象,由PID控制器进行调节控制;
步骤(2):建立控制对象的T-S模糊模型:提取训练样本,以(ck,uk)为输入数据,xk为输出数据,通过训练样本集建立控制对象的T-S模糊模型;T-S模糊模型训练完成后,通过测试样本测试T-S模糊模型精度,将测试样本中的(ck,uk)作为模糊模型输入数据,计算T-S模糊模型输出tk与测试样本xk的均方根误差,均方根误差越小,说明T-S模糊模型与实际机组水轮机调速系统中的控制对象的误差越小,从而选择最优T-S模糊模型。
步骤(3):建立基于T-S模糊模型的仿真系统:用得到的控制对象的T-S模糊模型替换控制对象,得到基于T-S模糊模型的水轮机调速系统的仿真系统,其中PID控制器按照实际机组水轮机调速系统PID控制器进行设计,Kp、Ki和Kd分别为比例、积分和微分增益,是需要整定的PID控制参数;
步骤(4):建立PID控制参数优化目标函数:采用离散形式ITAE指标作为控制参数优化的目标函数,目标函数定义为:
其中c(k)为频率扰动值,为一常数,x为机组频率,Ns为采样点数,T为时间序列。
步骤(5):运用GSA算法求解PID控制参数优化目标函数,得到最优控制参数。fITAE(Kp,Ki,Kd)是PID控制参数的函数,当PID控制参数发生变化时,利用步骤(3)中建立的仿真系统得到一组机组频率响应曲线,计算对应的fITAE值,极小化目标函数,得到最优控制参数。
A:算法初始化:设置算法参数,包括群体规模Np、最大迭代次数Nmax_it、引力常数初值G0、衰减指数α;确定优化控制参数范围,在此区间初始化群体中粒子的位置向量分别代表控制参数Kp、Ki和Kd;初始化历史最优目标函数值fbest和最优解Xbest,设置迭代次数t=0;
B:计算目标函数值:从位置向量中取出控制参数,将控制参数输入步骤(3)中建立的基于T-S模糊模型的水轮机调速系统的仿真系统,得到仿真机组频率,按照步骤(4)计算粒子目标函数值fi,i=1,…,Np;
C:求出群体最优目标函数值fbest_current=min({fi,i=1,…,Np})和对应的粒子位置Xbest_current,如果fbest_current<fbest,则fbest<fbest_current,Xbest=Xbest_current;
D:计算粒子质量Mi,i=1,…,Np:
其中,fworst_current=max({fi,i=1,…,Np});
E:计算引力Fi d:
其中,randj为[0,1]范围随机数,d=1,…,3;
F:计算加速度:
G:更新粒子速度和粒子位置
其中,randi为[0,1]范围随机数;
H:t=t+1;如果t<Nmax_it,转至Step 2,否则结束循环,得到最优目标函数值fbest和最优位置向量Xbest,所述最优位置向量即为最优控制参数向量。
对于本领域技术人员而言,显然本发明不限于上述示范性实施例的细节,而且在不背离本发明的精神或基本特征的情况下,能够以其他的具体形式实现本发明。因此,无论从哪一点来看,均应将实施例看作是示范性的,而且是非限制性的,本发明的范围由所附权利要求而不是上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (6)
1.一种水轮机调速系统控制参数的自动整定方法,其特征在于,包含以下步骤:
(1)采集实际机组水轮机调速系统动态过程数据;
(2)建立控制对象的T-S模糊模型;
(3)建立基于T-S模糊模型的仿真系统;
(4)建立PID控制参数优化目标函数;
(5)运用GSA算法求解PID控制参数优化目标函数,得到最优控制参数,fITAE(Kp,Ki,Kd)是PID控制参数的函数,当PID控制参数发生变化时,利用步骤(3)中建立的仿真系统得到一组机组频率响应曲线,计算对应的fITAE值,极小化目标函数,得到最优控制参数。
2.根据权利要求1所述的一种水轮机调速系统控制参数的自动整定方法,其特征在于,所述步骤(1)具体是:搜集不同PID控制参数、不同频率扰动值所对应的动态过程数据样本(ck,uk,xk)j组成数据样本集,ck为k时刻的频率扰动值,uk为k时刻调速器输出的控制量序列,xk为k时刻的机组频率,k=1,..,Ns,Ns为动态过程最大采样点数;j=1,…,Nc,表示不同的动态过程,总共采集Ns*Nc个样本组成样本集,将样本集中若干动态过程样本作为训练样本;其中所述水轮机调速系统包括PID控制器、接力器、引水系统、水轮机,其中,接力器、引水系统、水轮机组成控制对象,由PID控制器进行调节控制。
3.根据权利要求1所述的一种水轮机调速系统控制参数的自动整定方法,其特征在于,所述步骤(2)具体是:提取训练样本,以(ck,uk)为输入数据,xk为输出数据,通过训练样本集建立控制对象的T-S模糊模型;T-S模糊模型训练完成后,通过测试样本测试T-S模糊模型精度,将测试样本中的(ck,uk)作为模糊模型输入数据,计算T-S模糊模型输出tk与测试样本xk的均方根误差,均方根误差越小,说明T-S模糊模型与实际机组水轮机调速系统中的控制对象的误差越小,从而选择最优T-S模糊模型。
4.根据权利要求1所述的一种水轮机调速系统控制参数的自动整定方法,其特征在于,所述步骤(3)具体是:用得到的控制对象的T-S模糊模型替换控制对象,得到基于T-S模糊模型的水轮机调速系统的仿真系统,其中PID控制器按照实际机组水轮机调速系统PID控制器进行设计,Kp、Ki和Kd分别为比例、积分和微分增益,是需要整定的PID控制参数。
5.根据权利要求1所述的一种水轮机调速系统控制参数的自动整定方法,其特征在于,所述步骤(4)具体是:采用离散形式ITAE指标作为控制参数优化的目标函数,目标函数定义为:其中c(k)为频率扰动值,为一常数,x为机组频率,Ns为采样点数,T为时间序列。
6.根据权利要求1所述的一种水轮机调速系统控制参数的自动整定方法,其特征在于,所述步骤(5)具体包含以下步骤:
A、算法初始化:设置算法参数,包括群体规模Np、最大迭代次数Nmax_it、引力常数初值G0、衰减指数α;确定优化控制参数范围,在此区间初始化群体中粒子的位置向量分别代表控制参数Kp、Ki和Kd;初始化历史最优目标函数值fbest和最优解Xbest,设置迭代次数t=0;
B、计算目标函数值:从位置向量中取出控制参数,将控制参数输入步骤(3)中建立的基于T-S模糊模型的水轮机调速系统的仿真系统,得到仿真机组频率,按照步骤(4)计算粒子目标函数值fi,i=1,…,Np;
C、求出群体最优目标函数值fbest_current=min({fi,i=1,…,Np})和对应的粒子位置Xbest_current,如果fbest_current<fbest,则fbest<fbest_current,Xbest=Xbest_current;
D、计算粒子质量Mi,i=1,…,Np:
E、其中,fworst_current=max({fi,i=1,…,Np});
计算引力Fi d:
其中,randj为[0,1]范围随机数,d=1,…,3;
F、计算加速度:
G、更新粒子速度和粒子位置 其中,randi为[0,1]范围随机数;
H、t=t+1;如果t<Nmax_it,转至Step2,否则结束循环,得到最优目标函数值fbest和最优位置向量Xbest,所述最优位置向量即为最优控制参数向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810211577.5A CN108427397A (zh) | 2018-03-15 | 2018-03-15 | 一种水轮机调速系统控制参数的自动整定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810211577.5A CN108427397A (zh) | 2018-03-15 | 2018-03-15 | 一种水轮机调速系统控制参数的自动整定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108427397A true CN108427397A (zh) | 2018-08-21 |
Family
ID=63158540
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810211577.5A Withdrawn CN108427397A (zh) | 2018-03-15 | 2018-03-15 | 一种水轮机调速系统控制参数的自动整定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108427397A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111027178A (zh) * | 2019-11-08 | 2020-04-17 | 华中科技大学 | 一种抽水蓄能机组调频参数优化方法 |
CN114065501A (zh) * | 2021-11-11 | 2022-02-18 | 中兵勘察设计研究院有限公司 | 生产线工艺仿真优化方法及系统 |
-
2018
- 2018-03-15 CN CN201810211577.5A patent/CN108427397A/zh not_active Withdrawn
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111027178A (zh) * | 2019-11-08 | 2020-04-17 | 华中科技大学 | 一种抽水蓄能机组调频参数优化方法 |
CN114065501A (zh) * | 2021-11-11 | 2022-02-18 | 中兵勘察设计研究院有限公司 | 生产线工艺仿真优化方法及系统 |
CN114065501B (zh) * | 2021-11-11 | 2024-07-09 | 中兵勘察设计研究院有限公司 | 生产线工艺仿真优化方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104533701B (zh) | 一种水轮机调速系统控制参数的自动整定方法 | |
CN106681150A (zh) | 基于闭环参考模型的分数阶pid控制器参数优化整定方法 | |
CN105425612B (zh) | 一种水轮机调节系统控制参数的优选方法 | |
CN108427397A (zh) | 一种水轮机调速系统控制参数的自动整定方法 | |
CN110262223A (zh) | 一种基于分数阶pid调速系统的水轮机综合性模型建模方法 | |
CN104734588B (zh) | 一种生物质气内燃发电机组转速控制方法 | |
CN114168796B (zh) | 一种建立飞行器高空气动力数据库的方法 | |
CN104616072B (zh) | 一种基于区间优化的提高谷氨酸发酵产物浓度的方法 | |
CN105807632A (zh) | 一种基于改进大林算法的加热炉温度控制器设计方法 | |
CN102509152A (zh) | 一种基于rbf神经网络的开关磁阻电机在线建模方法 | |
CN111502914B (zh) | 基于线性变参数系统的风机变桨控制器设计方法 | |
CN108398878A (zh) | 基于血红细胞繁殖模型的分数阶pd控制器设计方法 | |
CN105389614B (zh) | 一种神经网络自更新过程的实现方法 | |
CN103809446B (zh) | 飞行器多回路模型簇颤振抑制复合频率鲁棒控制器设计方法 | |
CN106406100A (zh) | 基于模糊自整定单神经元pid控制的转子动平衡控制系统及其方法 | |
CN108549207A (zh) | 一种水轮机调节系统控制参数的方法 | |
CN113328435A (zh) | 一种基于强化学习的主动配电网有功无功联合控制方法 | |
CN113239643A (zh) | 一种海上风电场自动发电控制流场动态建模方法 | |
CN103809434B (zh) | 纵向飞行模型簇复合根轨迹多级pid控制器设计方法 | |
CN109653947A (zh) | 一种侧偏型小型风力发电机的动态特性评估方法 | |
CN103809442B (zh) | 飞行器多回路模型簇复合频率鲁棒控制器设计方法 | |
CN106940526A (zh) | 一种去伪控制的碳纤维凝固浴工艺 | |
CN110244852A (zh) | 一种大脑情绪控制器控制参数整定方法 | |
CN103823376B (zh) | 纵向飞行模型簇复合pid控制器设计方法 | |
CN103809452B (zh) | 纵向飞行模型簇颤振抑制复合根轨迹多级pid鲁棒控制器设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20180821 |
|
WW01 | Invention patent application withdrawn after publication |