CN110486219B - 一种抑制调压室引起被控参数低频波动的pid控制方法 - Google Patents
一种抑制调压室引起被控参数低频波动的pid控制方法 Download PDFInfo
- Publication number
- CN110486219B CN110486219B CN201910657069.4A CN201910657069A CN110486219B CN 110486219 B CN110486219 B CN 110486219B CN 201910657069 A CN201910657069 A CN 201910657069A CN 110486219 B CN110486219 B CN 110486219B
- Authority
- CN
- China
- Prior art keywords
- pid
- controlled parameter
- water head
- signal
- fundamental wave
- 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
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03B—MACHINES OR ENGINES FOR LIQUIDS
- F03B15/00—Controlling
- F03B15/02—Controlling by varying liquid flow
- F03B15/04—Controlling by varying liquid flow of turbines
- F03B15/06—Regulating, i.e. acting automatically
- F03B15/08—Regulating, i.e. acting automatically by speed, e.g. by measuring electric frequency or liquid flow
- F03B15/12—Regulating, i.e. acting automatically by speed, e.g. by measuring electric frequency or liquid flow with retroactive action
-
- 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
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Feedback Control In General (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明公开了一种抑制调压室引起被控参数低频波动的PID控制方法,在传统PID调速器模型的基础上,包括如下步骤:采用傅里叶级数从被控参数测量信号i(k)中提取被控参数主波信号im(k),将其输入至PID控制器,得到PID输出值yPID(k);采用傅里叶级数从水头信号中提取水头基波信号ht(k),并将其输入至基波附加控制器中,得到水头附加控制值ya(k),然后与PID输出yPID(k)叠加后作为PID调速器的输出信号yout(k)。能使调压室引起的机组频率、功率、水头等被控参数波动迅速衰减,有效提高系统调节品质及稳定性。
Description
技术领域
本发明属于水轮机调节方法技术领域,涉及一种抑制调压室引起被控参数低频波动的PID控制方法。
背景技术
含调压室引水系统的水力发电机组在调节过程中常出现机组频率、功率、水头等参数低频率波动的现象,通常认为该波动由调压室与水库之间周期性水体输运引起,且各参数波动周期相同,通常在几十秒至几百秒之间。目前,常规PID调速器对机组频率、功率、水头等参数低频率波动的抑制作用较小,当该波动剧烈时对机组调节品质造成明显影响,甚至造成机组调节过程不稳定等问题。
发明内容
本发明的目的是提供一种抑制调压室引起被控参数低频波动的PID控制方法,解决了现有技术中存在的PID调速器无法有效抑制调压室引起的被控参数低频波动问题。
本发明所采用的技术方案是,一种抑制调压室引起被控参数低频波动的PID控制方法,包括如下步骤:
步骤1、采用傅里叶级数从被控参数测量信号i(k)中提取被控参数主波信号im(k),将其输入至PID控制器,得到PID输出值yPID(k);
步骤2、采用傅里叶级数从水头信号中提取水头基波信号ht(k),并将其输入至基波附加控制器中,得到水头附加控制值ya(k),然后与PID输出yPID(k)叠加后作为PID调速器的输出信号yout(k)。
本发明的特点还在于:
步骤1具体包括:
步骤1.1、根据傅里叶级数得到调压室引起的被控参数振荡基波过程为:
采用和式代替积分,对式(12)进行离散,得到被控参数基波信号ib(k):
其中,k为采样次数,ΔT为采样周期,N为1周期内的采样次数,Tb为调压室引起的被控参数振荡周期,t为时间,a1为余弦分量系数,b1为正弦分量系数,j为计数器;
步骤1.2、由被控参数测量信号i(k)减去被控参数基波信号ib(k)得到被控参数主波信号im(k):
im(k)=i(k)-ib(k) (14);
其中,在PID调速器工作在频率模式时,被控参数i为机组频率x,PID调速器工作在功率模式时,被控参数i为机组功率p;
步骤1.3、将被控参数主波信号im(k)输入至PID控制器进行比例、积分、微分运算,得到PID输出值yPID(k):
yPID(k)=fPID[im(k)] (15);
其中,fPID为PID计算关系;im(k)为第k个采样周期的被控参数主波信号。
步骤2具体为:
步骤2.1、采用傅里叶级数从水头信号中提取水头基波信号ht(k),并将其输入至基波附加控制器中进行比例运算得到水头附加控制值ya(k)为:
其中,ht(k)为第k个采样周期的水头基波信号,C为控制系数,Tb水头偏差的振荡周期,与被控参数振荡周期相同,a′1(k)与b′1(k)分别为第k个采样周期的水头基波系数:
其中,h(j)为第j个采样周期的水头测量信号,;
步骤2.2、将水头附加控制值ya(k)作为接力器行程控制附加量叠加至PID控制器输出值中,得到PID调速器的输出信号yout(k):
yout(k)=yPID(k)+ya(k) (18)。
水头偏差和被控参数的振荡周期Tb均为:
其中,L为隧洞长度,a为隧洞中的水击波速,F为隧洞断面面积,g为重力加速度,FT为调压室断面面积。
水头偏差和被控参数的振荡周期Tb还可以根据测量得到。
本发明的有益效果是:本发明的PID控制方法,通过傅里叶级提取被控参数的主波信号与机组水头的基波信号,将被控参数主波信号输入至PID控制器中,将机组水头基波信号进行幅值修正后附加至PID输出信号中,能使调压室引起的机组频率、功率、水头等被控参数波动迅速衰减,有效提高系统调节品质及稳定性。
附图说明
图1是本发明一种抑制调压室引起被控参数低频波动的PID控制方法的上游调压室的系统结构示意图;
图2是本发明一种抑制调压室引起被控参数低频波动的PID控制方法的下游调压室的系统结构示意图;
图3是本发明一种抑制调压室引起被控参数低频波动的PID控制方法的传统PID调速器数学模型框图;
图4是本发明一种抑制调压室引起被控参数低频波动的PID控制方法的PID调速器数学模型框图;
图5是本发明一种抑制调压室引起被控参数低频波动的PID控制方法的负荷扰动过渡过程曲线图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种抑制调压室引起被控参数低频波动的PID控制方法,建立了含调压室的系统数学模型,对于上游调压室,系统结构如图1所示,其通过隧洞与上游水库连接,通过压力管道与下游水库相接;对于下游调压室,系统结构如图2所示,其通过隧洞与下游水库连接,通过压力管道与上游水库相接。
在电站中调压室一般隧洞较长,压力管道相对较短。在保证精度的前提下为便于分析,忽略压力管道内水体造成的水击作用,同时将隧洞中的水击按照弹性水击进行考虑。
忽略隧洞中的能量损失,调压室处隧洞压力上升相对值可表示为:
其中,L为隧洞长度,m;a为隧洞中的水击波速,m/s;s为拉式算子;q为调压室处隧洞流量变化量相对值;Tw为隧洞水流惯性时间常数,s;上游调压室时取负,下游调压室时取正。
隧洞水流惯性时间常数表示为
其中,Qr为额定工况下隧洞流量,m3/s;Hr为额定工况下调压室处隧洞压力,m;F为隧洞断面面积,m2;g为重力加速度,m/s2。
不考虑调压室内水体的压缩性及调压室壁面变形,流入调压室的流量为:
其中,FT为调压室断面面积,m2;Hs为调压室水位,m。
由式(3)可得流入调压室的流量偏差相对值为:
其中,hs(s)为调压室水位偏差相对值。
根据水体连续性可得:
q(s)=q2(s)±q1(s) (5);
其中,q2为流入水力机组流量偏差相对值;上游调压室时取正,下游调压室时取负。
由式(1)、(2)、(4)、(5)得到引水系统数学模型为:
上游调压室取负,下游调压室取正。
对式(6)进行反拉氏变换可得机组流量相对值单位阶跃变化时,调压室处隧洞压力变化解析式:
上游调压室取负,下游调压室取正。
在不考虑压力管道水流惯性作用的情况下,机组水头表示为:
由式(8)可以看出,机组流量发生单位阶跃变化后,水头偏差在0附近发生周期性振荡,其振荡周期为:
对于线性时不变系统,输入正弦函数信号时输出则为同频率正弦函数信号。因此,机组频率、功率等被控参数的振荡周期必然也为Tb,被控参数和水头偏差的振荡周期Tb也可以根据测量得到。
在实际电站运行中,机组水头以及被控参数的变化过程不仅受调压室影响,同时也受运行工况、调速器、引水系统特性、负载特性等多方面因素共同影响。考虑到调压室引起的各参数变化周期长且占主要成分,而其他因素引起的各参数变化周期较短且含量较少,因此可将调压室引起的水头变化看作基波,采用傅里叶级数进行分离。而调压室引起的各参数低频波动主要是由于机组水头波动引起,因此可以将水头的基波分量作为低频波动的被控量,采用比例策略获取控制附加量。
调速器通过调整导叶开度,一方面要使被控参数到达设定值,另一方面要使调压室引起的各参数低频波动能够快速衰减。其中,使被控参数稳定在设定值是调速器基本功能,其控制方法与传统PID调速器相同,考虑到调压室引起的各参数基波分量稳态值均为0,对系统稳态无影响,仅对系统瞬态过程有影响,因此可以将被控参数中除基波分量以外的主波分量作为反馈量进行控制。
本发明一种抑制调压室引起被控参数低频波动的PID控制方法,在传统PID调速器模型的基础上,PID调速器包括PID控制器和基波附加控制器两部分,包括如下步骤:
步骤1、采用傅里叶级数从被控参数测量信号i(k)中提取被控参数主波信号im(k),将其输入至PID控制器,得到PID输出值yPID(k)。
傅里叶级数指出,对于周期为2π的函数i((x)可以展为:
其中,系数a0,aj,bj表达式如下:
步骤1.1、对于调压室引起的各参数波动过程即可以看做j为1时的分量,所以根据傅里叶级数得到调压室引起的被控参数振荡基波过程为:
采用和式代替积分,对式(12)进行离散,得到被控参数基波信号ib(k):
步骤1.2、由被控参数测量信号i(k)减去被控参数基波信号ib(k)得到被控参数主波信号im(k):
im(k)=i(k)-ib(k) (14);
其中,在PID调速器工作在频率模式时,被控参数i为机组频率x,PID调速器工作在功率模式时,被控参数i为机组功率p。
步骤1.3将被控参数主波信号im(k)输入至PID控制器进行比例、积分、微分运算,得到PID输出值yPID(k):
yPID(k)=fPID[im(k)] (15);
其中,fPID为PID计算关系;im(k)为第k个采样周期的被控参数主波信号值。
步骤2、采用傅里叶级数从水头信号中提取水头基波信号ht(k),并将其输入至基波附加控制器中进行比例运算得到水头附加控制值ya(k),然后与PID输出yPID(k)叠加作为PID调速器的输出信号yout(k)。
步骤2.1、采用傅里叶级数从水头信号中提取水头基波信号ht(k),并将其输入至基波附加控制器中进行比例运算得到水头附加控制值ya(k)为:
由上式可得在稳定工况时水头附加控制值ya为0,对机组稳态特性不造成影响。
其中,ht(k)为第k个采样周期的水头基波信号,C为控制系数,Tb水头偏差的振荡周期,与被控参数振荡周期相同,a′1(k)与b′1(k)分别为第k个采样周期的水头基波系数:
其中,h(j)为第j个采样周期的水头测量信号;
步骤2.2、将水头附加控制值ya(k)作为接力器行程控制附加量叠加至PID控制器输出值中,得到调速器的输出信号yout(k):
yout(k)=yPID(k)+ya(k) (18);
实施例
某电站为“一洞一机”引水,设有下游调压室,水轮机型号为HLD563-F13,转轮直径7.1m,机组额定出力为367MW,额定转速为100r/min,额定水头为80m,额定流量为506m3/s,机械惯性时间常数Ta=9.52s,尾水隧洞长411.06m,内径12.73m,水击波速800m/s,调压室横截面积970m2,压力管道长254.35m,内径11m,水击波速1000m/s,初始工况下机组水头为额定水头,发电机及负载自调节系数eg=1.0。
采用传统PID控制算法的控制系统框图如图3所示,调差系数bp=0.04,接力器响应时间常数Ty=0.2s,通过寻优获得10%负载扰动下调速器最优调节参数KP=0.05,KI=0.05s-1,KD=3s。其中引水系统采用特征线法求解,水轮机采用综合特性曲线进行查表获得其特性,系统10%负荷扰动过渡过程曲线如图5中虚线所示;采用本发明PID控制方法的控制系统框图如图4所示,经试算C取值为0.5,过渡过程曲线如图5中实线所示。由图5可以看出,采用本发明的PID控制方法时,调压室引起的功率波动迅速衰减,系统调节品质及稳定性得到明显改善。
通过以上方式,本发明的PID控制方法,通过傅里叶级提取被控参数的主波信号与机组水头的基波信号,将被控参数主波信号输入至PID控制器中,将机组水头基波信号进行幅值修正后附加至PID输出信号中,能使调压室引起的机组频率、功率、水头等被控参数波动迅速衰减,有效提高系统调节品质及稳定性。
Claims (3)
1.一种抑制调压室引起被控参数低频波动的PID控制方法,其特征在于,包括如下步骤:
步骤1、采用傅里叶级数从被控参数测量信号i(k)中提取被控参数主波信号im(k),将其输入至PID控制器,得到PID输出值yPID(k);
步骤1.1、根据傅里叶级数得到调压室引起的被控参数振荡基波过程为:
采用和式代替积分,对式(12)进行离散,得到被控参数基波信号ib(k):
其中,k为采样次数,ΔT为采样周期,N为1周期内的采样次数,Tb为调压室引起的被控参数振荡周期,t为时间,a1为余弦分量系数,b1为正弦分量系数,j为计数器;
步骤1.2、由所述被控参数测量信号i(k)减去被控参数基波信号ib(k)得到被控参数主波信号im(k):
im(k)=i(k)-ib(k) (14);
其中,在PID调速器工作在频率模式时,被控参数i为机组频率x,PID调速器工作在功率模式时,被控参数i为机组功率p;
步骤1.3、将所述被控参数主波信号im(k)输入至PID控制器进行比例、积分、微分运算,得到PID输出值yPID(k):
yPID(k)=fPID[im(k)] (15);
其中,fPID为PID计算关系;im(k)为第k个采样周期的被控参数主波信号;
步骤2、采用傅里叶级数从水头信号中提取水头基波信号ht(k),并将所述水头基波信号ht(k)输入至基波附加控制器中,得到水头附加控制值ya(k),然后与PID输出值yPID(k)叠加后作为PID调速器的输出信号yout(k);
步骤2.1、采用傅里叶级数从水头信号中提取水头基波信号ht(k),并将所述水头基波信号ht(k)输入至基波附加控制器中进行比例运算得到水头附加控制值ya(k):
其中,ht(k)为第k个采样周期的水头基波信号,C为控制系数,Tb水头偏差的振荡周期,与被控参数振荡周期相同,a′1(k)与b′1(k)分别为第k个采样周期的水头基波系数:
其中,h(j)为第j个采样周期的水头测量信号;
步骤2.2、将所述水头附加控制值ya(k)作为接力器行程控制附加量叠加至PID输出值yPID(k)中,得到PID调速器的输出信号yout(k):
yout(k)=yPID(k)+ya(k) (18)。
3.根据权利要求1所述的一种抑制调压室引起被控参数低频波动的PID控制方法,其特征在于,所述水头偏差和被控参数的振荡周期Tb根据测量得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910657069.4A CN110486219B (zh) | 2019-07-19 | 2019-07-19 | 一种抑制调压室引起被控参数低频波动的pid控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910657069.4A CN110486219B (zh) | 2019-07-19 | 2019-07-19 | 一种抑制调压室引起被控参数低频波动的pid控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110486219A CN110486219A (zh) | 2019-11-22 |
CN110486219B true CN110486219B (zh) | 2020-09-25 |
Family
ID=68547519
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910657069.4A Active CN110486219B (zh) | 2019-07-19 | 2019-07-19 | 一种抑制调压室引起被控参数低频波动的pid控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110486219B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112686539A (zh) * | 2020-12-29 | 2021-04-20 | 华润电力技术研究院有限公司 | 一种热力系统的性能计算方法、装置及设备 |
CN113013894B (zh) * | 2021-02-25 | 2023-05-23 | 云南电网有限责任公司电力科学研究院 | 一种水轮机调速器控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105065191A (zh) * | 2015-07-16 | 2015-11-18 | 国家电网公司 | 一种加快高水头水电机组减负荷后系统稳定的方法 |
CN107288811A (zh) * | 2016-04-05 | 2017-10-24 | 东方电机控制设备有限公司 | 利用调压井水位反馈的一管多机水轮发电系统调速控制回路及方法 |
CN107514333A (zh) * | 2017-07-19 | 2017-12-26 | 西安理工大学 | 具有抗速度饱和的微分优先水轮机pid调速器控制方法 |
CN108131238A (zh) * | 2017-12-21 | 2018-06-08 | 国网福建省电力有限公司 | 一种抑制水锤压力波动的新型pid控制方法 |
-
2019
- 2019-07-19 CN CN201910657069.4A patent/CN110486219B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105065191A (zh) * | 2015-07-16 | 2015-11-18 | 国家电网公司 | 一种加快高水头水电机组减负荷后系统稳定的方法 |
CN107288811A (zh) * | 2016-04-05 | 2017-10-24 | 东方电机控制设备有限公司 | 利用调压井水位反馈的一管多机水轮发电系统调速控制回路及方法 |
CN107514333A (zh) * | 2017-07-19 | 2017-12-26 | 西安理工大学 | 具有抗速度饱和的微分优先水轮机pid调速器控制方法 |
CN108131238A (zh) * | 2017-12-21 | 2018-06-08 | 国网福建省电力有限公司 | 一种抑制水锤压力波动的新型pid控制方法 |
Non-Patent Citations (4)
Title |
---|
"不同尾水系统对水电站系统稳定性影响";韩军颖;《工程科技Ⅱ辑》;20190131;全文 * |
"带调压井水轮机微分代数模型建模与分析";何建宇;《工程科技Ⅱ辑》;20180131;全文 * |
"机组调速系统引入调压井涌浪反馈对减负荷过程的影响";蔡龙等;《水电能源科学》;20160731;第175-177页 * |
"混流式水轮机调节系统建模及非线性动力分析";丁聪;《工程科技Ⅱ辑》;20140228;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110486219A (zh) | 2019-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110486219B (zh) | 一种抑制调压室引起被控参数低频波动的pid控制方法 | |
CN103529869A (zh) | 压力控制装置、流量控制装置、压力控制方法及流量控制方法 | |
CN112651180B (zh) | 一种一管多机水电机组调节系统微分方程计算方法 | |
CN108954286B (zh) | 一种锅炉输入动态加速指令的计算方法及系统 | |
CN103850872A (zh) | 转桨式水轮机协联关系综合优化试验方法 | |
Mao et al. | Improving internal fluid stability of pump turbine in load rejection process by co-adjusting inlet valve and guide vane | |
CN109631007A (zh) | 一种发电机组锅炉主控系统燃料反馈信号优化处理方法 | |
Zhu et al. | Dynamic modeling and eigen analysis of adjustable-speed pumped storage unit in pumping mode under power regulation | |
CN107514333A (zh) | 具有抗速度饱和的微分优先水轮机pid调速器控制方法 | |
CN105543443B (zh) | 一种用于转炉煤气加压机防喘振控制系统的方法 | |
CN105065191A (zh) | 一种加快高水头水电机组减负荷后系统稳定的方法 | |
CN109599889B (zh) | 基于模糊自抗扰的不平衡电压下的穿越控制方法、系统 | |
CN112947076B (zh) | 一种一管多机水电机组协同控制器设计方法 | |
CN111987736B (zh) | 一种补偿水轮机水锤效应的直流频率限制控制器控制方法 | |
CN110417047B (zh) | 基于复转矩系数分析双馈风机ssci阻尼特性的方法 | |
Shi et al. | Nonlinear dynamic characteristics analysis and adaptive avoid vortex-coordinated optimal control of hydropower units under grid connection | |
CN110941868A (zh) | 一种水电站调压室临界稳定断面的计算方法 | |
CN107110104B (zh) | 用于稳定具有s特性的液压机的旋转速度的方法以及用于将液压能转换成电能的设备 | |
CN112271759B (zh) | 多调频资源框架下风电调频参数的优化方法、系统及装置 | |
Ali et al. | Optimization of water level control systems using anfis and fuzzy-pid model | |
Salhi et al. | Modeling and regulation of a micro hydroelectric power plant | |
CN112952864A (zh) | 一种基于描述函数的vsc单边限幅参与的次同步振荡分析方法 | |
CN114237312A (zh) | 一种有无压管道混连的输水发电系统阀门及机组联合精准调控方法 | |
CN219676488U (zh) | 含给水动态前馈的除氧器水位三冲量控制装置 | |
CN107288811A (zh) | 利用调压井水位反馈的一管多机水轮发电系统调速控制回路及方法 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |