CN103590969B - 基于多工况时域响应的pid水轮机调速器参数优化方法 - Google Patents
基于多工况时域响应的pid水轮机调速器参数优化方法 Download PDFInfo
- Publication number
- CN103590969B CN103590969B CN201310588803.9A CN201310588803A CN103590969B CN 103590969 B CN103590969 B CN 103590969B CN 201310588803 A CN201310588803 A CN 201310588803A CN 103590969 B CN103590969 B CN 103590969B
- Authority
- CN
- China
- Prior art keywords
- operating mode
- model
- parameter
- pid
- water
- 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.)
- Expired - Fee Related
Links
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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/20—Hydro energy
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于多工况时域响应的PID水轮机调速器参数优化方法,包括以下步骤:(1)建立水轮机调节***模型,具体包括:PID水轮机调速器模型,水轮机-引水管道模型和发电机模型;(2)通过现场试验和参数辨识的方式求得水轮机调节***模型参数;(3)设计基于多工况时域响应的综合适应度函数,确定优化目标为:找到最优PID控制参数,使得综合适应度函数值达到最小;(4)利用智能优化算法求取最优控制参数。该方法求取的控制参数使得***能够在不同的工况下都保持满意的动态特性,增强了***鲁棒性。
Description
技术领域
本发明涉及水轮机调节技术,具体涉及一种基于多工况时域响应的PID水轮机调速器参数优化方法。
背景技术
水轮机调节的基本任务是根据电网负荷的变化,不断调节发电机组的有功功率输出,并维持机组转速在规定的范围内,其调节对象包括:压力引水***和泄水***、水轮机、发电机和机组并入运行的电网。由于该对象是一个将水力过程、机械过程和电气过程综合于一体、彼此密切联系的结构复杂和参数时变的被控制***,在水电站控制中通常采用比例积分微分(PID)调速器与其构成闭环水轮机调节***,以确保该***安全稳定运行,进而为电网提供高质量的供电。
对水轮机调节***而言,水轮机调速器是保证其拥有满意的动态特性并可靠运行的重要控制设备,而调速器的控制参数选择是否得当则决定了调速器控制品质的优劣。不合理的参数选择会导致水电机组性能不佳,限制其调节能力的发挥,并存在严重的安全问题。
许多先进的理论和方法已广泛用于PID水轮机调速器参数优化,但大多数方法都是在***处于某一特定工况的基础上进行控制参数优化,这样得到的最优参数不能保证***在其它工况下也拥有满意的调节品质,因而导致水轮机调节***对工况变化的鲁棒性较差。
发明内容
针对现有技术的不足,本发明的目的在于提出一种基于多工况时域响应的PID水轮机调速器参数优化方法,该方法求取的控制参数使得***能够在不同的工况下都保持满意的动态特性,***鲁棒性得以增强。
为实现以上发明目的,本发明采用以下技术方案:
一种基于多工况时域响应的PID水轮机调速器参数优化方法,具体包括以下步骤:
(1)在某一具体工况Cm下,根据最大超调量M,稳定时间T以及振荡次数N这三个时域响应指标,对一组给定的PID控制参数K,定义单一工况适应度函数为:
其中,系数ai用于各个指标之间的加权,i=1,2,3,4;函数e(·)表示频率偏差;K为PID控制器比例系数Kp、积分系数Ki和微分系数Kd构成的向量, t表示时间。
(2)在多个工况Cm(m=1,2,...,M)的情况下,定义PID控制参数K的综合适应度函数为:
(3)确定优化目标为:找到最优PID控制参数K*,使得综合适应度函数值达到最小,即:
与现有技术相比,本发明具有以下有益效果:使水轮机调节***在多个工况条件下能够保持满意的动态特性,因而显著提高了***对工况变化的鲁棒性。
附图说明
参照下面的说明,结合附图,可以对本发明有最佳的理解。在附图中,相同的部分可由相同的标号表示。
图1为PID控制器模型框图
图2为液压***模型框图
图3为水轮机-引水管道模型框图
图4为发电机模型框图
图5为水轮机调节***模型框图
图6为参数优化策略图
图7(a)为C1工况时各单一工况最优PID控制参数的***动态响应过程比较图
图7(b)为C2工况时各单一工况最优PID控制参数的***动态响应过程比较图
图7(c)为C3工况时各单一工况最优PID控制参数的***动态响应过程比较图
图8(a)为C1工况时多工况综合最优PID控制参数的***动态响应过程比较图
图8(b)为C2工况时多工况综合最优PID控制参数的***动态响应过程比较图
图8(c)为C3工况时多工况综合最优PID控制参数的***动态响应过程比较图
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及示例性实施例,对本发明进行进一步详细说明。应当理解,此处所描述的示例性实施例仅用以解释本发明,并不用于限定本发明的适用范围。
本发明提出的基于多工况时域响应的PID水轮机调速器参数优化方法,其基本思想为:以非线性水轮机调节***模型为基础,根据***在多个工况下的时域响应指标定义综合适应度函数,并利用智能优化算法对PID水轮机调速器控制参数进行综合优化。具体包括以下步骤:
(1)建立水轮机调节***非线性模型,包括PID水轮机调速器模型、水轮机-引水管道模型以及发电机模型。
(1.1)PID水轮机调速器模型
水轮机调速器由PID控制器和液压***构成。
PID控制器模型框图如图1所示,其控制规律描述如下:
其中,s为拉格朗日算子;Kp、Ki和Kd分别为调速器比例、积分和微分系数;e为频率偏差;Fg和Ft分别为给定频率和机组频率;Ff为人工频率死区;T1v为微分衰减常数;bp为永态转差系数;Yg、Yt和Ypid分别为给定开度值、反馈开度值和PID控制器输出开度值;PIDMAX和PIDMIN分别为PID输出饱和最大值与最小值。
液压***模型框图如图2所示,其传递函数如下:
其中,Y为液压***输出,即接力器行程;Ty为接力器反应时间;um为主配压阀死区;ρmax和ρmin分别为开度速率饱和最大值与最小值;Ymax和Ymin分别为开度位置饱和最大值与最小值。
(1.2)水轮机-引水管道模型
水轮机-引水管道模型框图如图3所示,其中主要关系描述如下。
水轮机流量q的表达式为:
其中,At为水轮机增益;h为水轮机水头。
引水管道内水头h和流量q的动态关系为:
其中,h0为初始稳态水头;Tw为水流惯性时间常数;Te为引水管道弹性时间常数;f为引水管道摩擦系数。
水轮机输出功率p为:
p=(q-qNL)h-DAtYΔx(5)
其中,qNL为空载流量;D为水轮机阻尼系数。
(1.3)发电机模型
发电机模型框图如图4所示,其传递函数为:
其中,Ta为机组惯性时间常数;eg为发电机自调节系数。
方程(1)到(6)完整描述了水轮机调节***的动态特性,***模型框图如图5所示,其中,xr为给定转速;x为发电机转速输出。
(2)通过现场试验求得水轮机调节***模型参数,包括死区、接力器反应时间、最短开启、关闭时间,并利用参数辨识方法对水轮机-引水管道模型中的水轮机增益、水流惯性时间常数、引水管道弹性时间常数、水轮机阻尼系数等模型参数进行辨识。
具体而言,死区、接力器反应时间、最短开启、关闭时间等非线性环节参数可以利用接力器关闭规律试验、调速器静特性试验等现场试验方法测得。另外,在不同工况Cm(m=1,2,...,M)下,对***施加阶跃频率变化以激励***,并实时记录水轮机导叶开度(或接力器行程)、水轮机输出转矩(或发电机电磁功率)、机组转速数据,直到动态过程重新恢复稳定。根据所记录的数据,采用参数辨识的方法可实现对不同工况下的水轮机增益、水流惯性时间常数、引水管道弹性时间常数、水轮机阻尼系数等水轮机-引水管道模型参数进行辨识。
在本实施例中,所述参数辨识方法可以采用以下文献中的方法:
用于水轮机-引水管道参数辨识的改进型人工鱼群算法[J].刘昌玉,何雪松,李崇威等.电力自动化设备,2013,33(11):61-65。
由于参数辨识方法属于现有技术,本发明在此不再赘述。
(3)设计基于多工况时域响应的综合适应度函数,确定优化目标为:找到最优PID控制参数,使得综合适应度函数值达到最小。
此为本发明的核心部分。其中,将PID调速器的比例系数Kp、积分系数Ki和微分系数Kd作为优化整定的控制参数,令水轮机调节***在不同的控制参数组合下会产生不同的调节品质,为了量化***在一组控制参数Kp、Ki和Kd下调节品质的优劣,需要进行适应度函数的计算。
考虑到机组在实际运行中,***调节品质主要体现在以下三个时域响应指标:最大超调量M,稳定时间T以及振荡次数N,对于一组给定的PID控制参数K,在某一具体工况Cm下,单一工况适应度函数定义为:
其中,系数ai用于各个指标之间的加权,i=1,2,3,4;函数e(·)表示频率偏差;t表示时间。
在此基础之上,PID控制参数K在多个工况Cm(m=1,2,...,M)情况下的综合适应度函数定义如下:
该综合适应度函数量化了***取同一组PID控制参数时,在多个工况下的调节品质。因此,确定优化目标为:找到最优PID控制参数K*,使得综合适应度函数值达到最小,即:
通过最小化Fit(K)函数,***将同时在多个工况下获得满意的动态特性。
(4)利用智能优化算法求取最优控制参数
图6展示了在多工况条件下进行PID调速器参数优化的基本策略。具体包括:
(4.1)任意取定一组控制参数K,并依次在不同工况Cm(m=1,2,...,M)下为步骤(1)和(2)建立的水轮机调节***模型施加阶跃扰动,记录***的时域响应过程;
(4.2)根据***时域响应过程,利用公式(7)计算各自单一工况下的适应度函数值并由公式(8)求得多工况下的综合适应度函数值Fit(K);
(4.3)采用智能优化算法对参数K进行更新,并不断重复步骤(4.1)与步骤(4.2),直至最小化综合适应度函数值Fit(K),求得最优控制参数K*。
在本实施例中,上述智能优化算法可以采用以下文献中的算法:
基于改进粒子群算法的水轮机调速***建模[J].刘昌玉,李崇威,洪旭钢等.水电能源科学,2011,29(12):124-127。
本领域技术人员应当理解,不失一般性,步骤(4.3)中还可以采用其他智能优化算法,由于改进粒子群算法等相关智能优化算法为现有技术,在此不再赘述。
接下来,以国内某水电站在三个工况(C1=18.1%开度、C2=19.8%开度、C3=16.3%开度)下的水轮机调节***模型为基础,通过与传统单一工况参数优化方法相比较,可以展示本发明提出方法的优越性。
表1第二列总结了在三个工况下,分别使用传统单一工况参数优化方法得到的最优控制参数,而图7展示了这三个最优控制参数在不同工况条件下的时域响应。很明显,在具体某一工况下(如:C1),***使用其对应的最优PID控制参数产生的动态响应过程相当优秀,但使用其他工况最优参数 所产生的动态响应过程则出现严重恶化。这表明,单一工况优化的PID控制参数只能适用于该工况,当工况发生变化时,***调节品质会急剧下降。
图8展示了由本发明方法获得的多工况综合最优PID控制参数在各个工况下的动态响应过程。可以观察到,虽然在每个工况下,***的调节品质不能到达最优,但同一组控制参数却使***在不同的工况下都保持满意的动态特性,***对工况变化的鲁棒性得以增强。
表1最优PID控制参数优化结果
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于多工况时域响应的PID水轮机调速器参数优化方法,具体包括以下步骤:
(1)在某一具体工况Cm下,根据最大超调量M、稳定时间T以及振荡次数N这三个时域响应指标,对一组给定的PID控制参数K,定义单一工况适应度函数为:
其中,系数ai用于各个指标之间的加权,i=1,2,3,4;函数e(·)表示频率偏差;K为PID控制器比例系数Kp、积分系数Ki和微分系数Kd构成的向量, t表示时间;
(2)在多个工况Cm(m=1,2,…,M)的情况下,定义PID控制参数K的综合适应度函数为:
(3)确定优化目标为:找到最优PID控制参数K*,使得综合适应度函数值达到最小,即:
2.根据权利要求1所述的优化方法,其中,所述步骤(1)之前还包括水轮机调节***模型建立步骤和求取模型参数的步骤,其中,
所述水轮机调节***模型建立步骤包括:建立水轮机调节***非线性模型,包括PID水轮机调速器模型、水轮机-引水管道模型以及发电机模型;
所述求取模型参数的步骤包括:通过现场试验求得水轮机调节***模型参数,包括死区、接力器反应时间、最短开启、关闭时间,并利用参数辨识方法对水轮机-引水管道模型中的水轮机增益、水流惯性时间常数、引水管道弹性时间常数、水轮机阻尼系数进行辨识。
3.根据权利要求1所述的优化方法,其中,所述步骤(3)之后还包括步骤(4):利用智能优化算法求取最优控制参数。
4.根据权利要求3所述的优化方法,其中,所述步骤(4)包括:
(4.1)任意取定一组控制参数K,并依次在不同工况Cm(m=1,2,…,M)下为水轮机调节***模型施加阶跃扰动,记录***的时域响应过程;
(4.2)根据***时域响应过程,利用公式(1)计算各自单一工况下的适应度函数值并由公式(2)求得多工况下的综合适应度函数值Fit(K);
(4.3)采用智能优化算法对参数K进行更新,并不断重复步骤(4.1)与步骤(4.2),直至最小化综合适应度函数值Fit(K),求得最优控制参数K*。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310588803.9A CN103590969B (zh) | 2013-11-20 | 2013-11-20 | 基于多工况时域响应的pid水轮机调速器参数优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310588803.9A CN103590969B (zh) | 2013-11-20 | 2013-11-20 | 基于多工况时域响应的pid水轮机调速器参数优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103590969A CN103590969A (zh) | 2014-02-19 |
CN103590969B true CN103590969B (zh) | 2016-04-20 |
Family
ID=50081248
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310588803.9A Expired - Fee Related CN103590969B (zh) | 2013-11-20 | 2013-11-20 | 基于多工况时域响应的pid水轮机调速器参数优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103590969B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104503260B (zh) * | 2014-12-05 | 2018-09-28 | 广东电网有限责任公司电力科学研究院 | 调速器参数设定方法及装置 |
CN104533701B (zh) * | 2014-12-23 | 2017-03-08 | 华中科技大学 | 一种水轮机调速***控制参数的自动整定方法 |
CN105114242B (zh) * | 2015-07-22 | 2017-07-28 | 重庆邮电大学 | 基于模糊自适应的depso算法的水轮机调速器参数优化方法 |
CN105138739B (zh) * | 2015-07-30 | 2019-02-12 | 南方电网科学研究院有限责任公司 | 一种计及死区效应的电力***最低频率快速计算方法 |
CN105425612B (zh) * | 2015-11-10 | 2018-02-23 | 华中科技大学 | 一种水轮机调节***控制参数的优选方法 |
CN105844044A (zh) * | 2016-04-06 | 2016-08-10 | 国家电网公司 | 水轮机调节***被控对象参数辨识满意度函数设计方法 |
CN106681133B (zh) * | 2016-12-29 | 2019-06-28 | 中国南方电网有限责任公司电网技术研究中心 | 一种水电机组模型改进型子空间闭环辨识方法 |
CN107514333B (zh) * | 2017-07-19 | 2019-02-01 | 西安理工大学 | 具有抗速度饱和的微分优先水轮机pid调速器控制方法 |
CN108223241B (zh) * | 2018-01-12 | 2020-11-10 | 国家电网公司 | 一种导叶反馈故障判断逻辑方法 |
CN110262223B (zh) * | 2019-07-16 | 2022-10-18 | 华北水利水电大学 | 一种基于分数阶pid调速***的水轮机综合性模型建模方法 |
CN112502883B (zh) * | 2019-09-16 | 2023-07-14 | 国电南瑞南京控制***有限公司 | 考虑频率稳定特性的水轮机调速器参数整定方法及装置 |
CN110778447A (zh) * | 2019-11-04 | 2020-02-11 | 云南电网有限责任公司电力科学研究院 | 一种混流式水轮发电机调速器空载pid参数优化方法 |
CN112012875B (zh) * | 2020-07-23 | 2022-03-11 | 国网江西省电力有限公司电力科学研究院 | 一种水轮机调节***pid控制参数的优化方法 |
CN112682392B (zh) * | 2020-12-15 | 2023-05-05 | 苏州伟创电气科技股份有限公司 | 液压控制方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101344065A (zh) * | 2008-08-26 | 2009-01-14 | 昆明理工大学 | 基于满意控制的水轮机调速***参数优化方法 |
CN102052233A (zh) * | 2010-12-23 | 2011-05-11 | 广东电网公司电力科学研究院 | 一种电力***稳定分析用水轮机调节***模组 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20040100613A (ko) * | 2003-05-23 | 2004-12-02 | 두산중공업 주식회사 | 발전소에서 터빈 제어 시스템 |
US7069171B2 (en) * | 2003-07-09 | 2006-06-27 | Cummins, Inc. | EMI noise filter for eddy current turbo speed sensor |
JP2010229962A (ja) * | 2009-03-30 | 2010-10-14 | Mitsubishi Electric Corp | 水車またはポンプ水車用調速制御装置 |
JP5868100B2 (ja) * | 2011-09-28 | 2016-02-24 | 日立三菱水力株式会社 | 水車またはポンプ水車の調速機 |
-
2013
- 2013-11-20 CN CN201310588803.9A patent/CN103590969B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101344065A (zh) * | 2008-08-26 | 2009-01-14 | 昆明理工大学 | 基于满意控制的水轮机调速***参数优化方法 |
CN102052233A (zh) * | 2010-12-23 | 2011-05-11 | 广东电网公司电力科学研究院 | 一种电力***稳定分析用水轮机调节***模组 |
Also Published As
Publication number | Publication date |
---|---|
CN103590969A (zh) | 2014-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103590969B (zh) | 基于多工况时域响应的pid水轮机调速器参数优化方法 | |
CN107476931B (zh) | 一种水轮机调速器pid参数优化方法及*** | |
CN104503228B (zh) | 一种水轮机调速器功率模式下一次调频稳定域的确定方法 | |
CN107800146A (zh) | 兼顾一次调频和超低频振荡抑制的调速器参数优化方法 | |
CN104389733B (zh) | 基于不确定性模型的水轮机pid调速器控制参数整定方法 | |
CN109103927B (zh) | 调速***pid控制器参数整定方法 | |
CN109445273B (zh) | 一种水轮机调速器控制方法 | |
CN104865830A (zh) | 一种机组负荷双重智能优化控制方法 | |
CN106837678A (zh) | 基于改进tlbo算法的水轮发电机组pid调速器参数优化 | |
Gil-González et al. | Passivity-based control and stability analysis for hydro-turbine governing systems | |
CN108490790A (zh) | 一种基于多目标优化的过热汽温自抗扰串级控制方法 | |
CN103606919B (zh) | 一种含多负荷扰动的电力***的频率控制器设计方法 | |
Jain et al. | Three area power system load frequency control using fuzzy logic controller | |
Ding et al. | Sliding mode/H∞ control of a hydro-power plant | |
CN105449698A (zh) | 一种新型水电机组负荷与频率控制器 | |
CN107942664A (zh) | 一种基于灵敏度分析的水轮机调速器参数整定方法及*** | |
CN107989665B (zh) | 供热发电机组的三侧一次调频控制方法及*** | |
CN108879721A (zh) | 一种基于风电消纳的频率控制方法 | |
CN112947076B (zh) | 一种一管多机水电机组协同控制器设计方法 | |
Zou et al. | Eigen-Structure Assignment-Based Differential Evolution Algorithm for TS Fuzzy Control Tuning Applied to Water-Turbine Governing System | |
CN104466979A (zh) | 一种输电网无功优化调控***及方法 | |
CN116799839A (zh) | 一种多资源协同功率控制方法、装置及电子设备 | |
CN105576712B (zh) | 一种基于自适应极点配置的风电场有功功率在线控制方法 | |
CN115031225A (zh) | 一种蒸汽热网储能辅助调峰控制方法、***、设备及介质 | |
CN115347622A (zh) | 水风光互补下水电站协调控制仿真方法、设备及存储介质 |
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: 20160420 Termination date: 20171120 |
|
CF01 | Termination of patent right due to non-payment of annual fee |