CN104408330A - 水轮机调节***高阶数学模型的降阶方法 - Google Patents
水轮机调节***高阶数学模型的降阶方法 Download PDFInfo
- Publication number
- CN104408330A CN104408330A CN201410826880.8A CN201410826880A CN104408330A CN 104408330 A CN104408330 A CN 104408330A CN 201410826880 A CN201410826880 A CN 201410826880A CN 104408330 A CN104408330 A CN 104408330A
- Authority
- CN
- China
- Prior art keywords
- order
- formula
- speed response
- rotating speed
- water turbine
- 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
Landscapes
- Control Of Water Turbines (AREA)
- Control Of Turbines (AREA)
Abstract
一种水轮机调节***高阶数学模型的降阶方法,包括如下步骤:①简化***分母的最高次项,实现***的一次降阶,得到一次低阶等效***式;②从波动叠加的角度出发,将一次低阶等效***式变换为两个2阶子***相加的形式;③进行拉普拉斯反变换,得到水轮机调节***的一次低阶等效***在负荷阶跃扰动下的机组转速响应波动方程;④对完整5阶***式进行二次降阶;⑤利用二阶***进行水轮机调节***转速响应调节品质的影响因素分析。其优点是:较之其他降阶方法,本方法是从完整的5阶***出发,包含了***的所有参数和影响因数,为之后的调节品质分析提供了可靠的保障。
Description
技术领域
本发明涉及一种水轮机调节***高阶数学模型的降阶方法。
背景技术
为了进行水轮机调节***的调节品质方面的理论研究,通常需要在如下的假设:(1)忽略水体和管壁的弹性,引水隧洞和压力管道都采用刚性水击模型;(2)小波动过渡过程采用稳态工况点局部线性化的方法,用水轮机的传递系数来表述水轮机稳态特性;(3)发电机采用一阶模型,电站单独运行(在孤立电网下运行);(4)忽略调速器的非线性特性:饱和特性和转速死区。下建立完整的数学模型,求得相应的综合传递函数。但即便在这样的简化下,引水发电***的每一个子环节均用最简单的数学模型(引水隧洞动力方程为一阶方程,调压室连续性方程为一阶,压力管道动力方程为一阶方程,水轮机方程为两个线性方程,发电机采用一阶模型,调速器方程为一阶理想调速器方程),对于有调压室的电站(图1),相应的调节***框图如图2、图3所示,得到的负荷阶跃扰动下的综合传递函数为5阶:
式(1)、式(2)中:
Mg为负荷阶跃的拉普拉斯变换,输入信号;
X为转速响应的拉普拉斯变换,输出信号;
bt为暂态转差系数;
Td为缓冲装置时间常数,s;
s为拉普拉斯算子;
ai、bi为系数,是管道参数、机组参数与调速器参数的函数;
mg0为负荷阶跃相对值,即阶跃后负荷Mg与阶跃前负荷Mg0的偏差相对值,mg0=(Mg-Mg0)/Mg0。
根据伽罗瓦理论,5次及以上的方程没有公式解、5次以下的方程有公式解,故式(2)不能直接求解。为了解决这一问题。前人从基本方程的简化出发,假设调压室水位按正弦规律波动且忽略压力管道的水流惯性、整个引水道的水头损失,或者忽略压力管道的水流惯性和水头损失,得到5阶以下的频率响应方程,从而实现频率响应的解析求解。但此种做法的缺陷在于没有建立包含水轮机调节***各个子环节的完整的数学模型,简化模型无法真实的反应原始***。
发明内容
本发明的目的是针对现有技术存在的缺陷,提供一种水轮机调节***高阶数学模型的降阶方法,得到低阶等效***。此等效***不仅可以直接进行解析求解,便于开展调节品质的理论分析,更能真实准确的反映原完整5阶***,保留完整***最主要的动态特征。
一种水轮机调节***高阶数学模型的降阶方法,包括如下步骤:
1、简化***分母的最高次项a0s5,实现***的一次降阶,得到一次低阶等效***式如下
式中,E表示等效***;
2、从波动叠加的角度出发,将一次低阶等效***式(3)变换为两个2阶子***相加的形式
XE(s)=X1(s)+X2(s) (4)
式中,X1(s)、X2(s)分别表示一次低阶等效***XE(s)的两个二阶子***;C1、C2、C3、C4、C5、C6、C7、C8为系数,均可用待定系数法求得;
3、对式(4)-式(6)进行拉普拉斯反变换,得到水轮机调节***的一次低阶等效***在负荷阶跃扰动下的机组转速响应波动方程:
xE(t)=x1(t)+x2(t) (7);
式中,xE(t)、x1(t)、x2(t)分别为***XE(s)、X1(s)、X2(s)对应的转速响应波动方程;t为时间变量;
4、对完整5阶***式进行二次降阶:x1(t)、x2(t)中有一个对应***的共轭主导复极点,依据极点的分布找出主导极点对应的子波动;
以二阶***来描述***的转速响应特性方程
式中,K2为振幅,δ2为衰减度,ω为角频率,为初相位;
并计算得到调节时间Tp
式中,Δ的取值对大电网为±0.2%、小电网为±0.4%;
5、利用二阶***进行水轮机调节***转速响应调节品质的影响因素分析。
本发明的工作原理与过程如下:本发明降阶的依据为***的主导极点。闭环***的极点决定了瞬态响应的类型,而瞬态响应的形状则主要取决于闭环***的零点。可以采用各种各样的校正手段来改善***的动态品质,但从本质上看,无非都是改变其极点和零点的分布情况。主导极点是闭环***所有极点中离虚轴最近的极点,它对***瞬态过程性能的影响最大,在整个响应过程中起着主要的决定性作用;工程上往往只用主导极点估算***的动态特性;在高阶***的降阶处理中,可以略去一些实际存在的极点和零点,得到近似的结果,其中,略去的零点离虚轴越远,计算结果与实际情况的差别越小。
经过分析可知,带调压室的水轮机调节***,始终存在一对共轭主导极点,此对共轭主导极点对应的是调压室水位波动引起的机组频率低频振荡,是决定***调节品质的主要方面。当删除***分母的最高次项a0s5时,***的共轭主导极点几乎不发生变化,亦即***的动态特性几乎不发生改变,所以得到的一次低阶等效***(4阶)可以真实反映进而代替完整5阶***。又由于共轭主导极点对应着一二阶子***包含在一次低阶等效***(4阶)中,所以通过波动叠加的方式将四阶***拆成两个二阶***、依据主导极点的分布找出其对应的二阶子***,用其来代替原5阶***,进行调节品质的研究。
本发明水轮机调节***高阶数学模型的降阶方法的优点是:本发明方法理论依据充分,实现方式简单。最终得到的二阶子***,包含了***的共轭主导极点,可以反映***的动态特性,故降阶结果是准确可靠的。并且较之其他降阶方法,本方法是从完整的5阶***出发,包含了***的所有参数和影响因数,为之后的调节品质分析提供了可靠的保障。
附图说明
图1为带调压室水电站引水发电***示意图。
图2为带调压室水电站水轮机调节***结构框图。
图3为管道子***结构框图。
图4a为完整5阶***的机组转速响应对比图(bt=0.5,Td=10s,nf=1.2);
图4b为一次低阶等效***的机组转速响应对比图(bt=0.5,Td=10s,nf=1.2)。
图5a为完整5阶***的机组转速响应对比图(bt=0.8,Td=15s,nf=1.3);
图5b为一次低阶等效***的机组转速响应对比图(bt=0.8,Td=15s,nf=1.3)。
图6a为机组转速响应曲线图;
图6b为机组转速响应的主波示意图;
图6c为机组转速响应的尾波示意图。
具体实施方式
下面结合附图,对本发明进行进一步说明。
某设有上游调压室方案的水电站水轮机调节***,极点分布的计算结果如表1所示。其中:nf=F/Fth为调压室面积放大系数,Fth为调压室临界稳定断面积。
表1
从表1可以看出:***总存在1对共轭复极点(s4、s5),其实部绝对值远小于其他3个极点的实部绝对值,表明其更靠近虚轴,是***的主导极点;其他3个极点(s1、s2、s3)可以是3个实极点,也可以是1个实极点与1对共轭复极点,且均是***的非主导极点。
删去完整5阶***X(s)表达式分母的最高次项a0s5,使分母变成再进行***极点分布的计算,结果如表2所示。
表2
从表2可以看出:***仍总存在1对共轭主导复极点(s4、s5),且其取值较原***几乎不发生变化,说明删去分母的最高次项a0s5不仅不会改变主导极点的类型和数量,而且几乎不会改变主导极点的大小;另2个非主导极点(s2、s3)为实极点。
用一次低阶等效***代替完整5阶***,即XE(s)≈X(s)。对比分析这2个***的机组转速响应:额定负荷运行突减10%额定负荷(即mg0=-0.1)下的机组转速瞬态响应过程如图4、图5所示,
从图4、图5可以看出:一次低阶等效***的波动响应曲线的趋势与5阶***一致,两者吻合的很好,说明此4阶***可以真实反映完整5阶***的波动特性,用其近似代替完整5阶***的处理方法是可行的。
然后按照前文所述的方法,求得共轭主导极点对应的二阶子***。在带有调压室的水轮机调节***中,此共轭主导二阶子***对应的是***的尾波,非主导二阶子***对应的是***的主波,示意图如图6。
从图6可以看出,尾波是负荷阶跃扰动下机组转速响应波动的主体部分和和决定调节品质的主要方面,可依据尾波开展水轮机调节***转速响应调节品质的分析。
Claims (1)
1.一种水轮机调节***高阶数学模型的降阶方法,其特征在于包括如下步骤:
①简化***分母的最高次项a0s5,实现***的一次降阶,得到一次低阶等效***式如下
式中,Mg为负荷阶跃的拉普拉斯变换,输入信号;
X为转速响应的拉普拉斯变换,输出信号;
E表示等效***;
s为拉普拉斯算子;
bt为暂态转差系数;
Td为缓冲装置时间常数,s;
ai、bi为系数,是管道参数、机组参数与调速器参数的函数;
mg0为负荷阶跃相对值,即阶跃后负荷Mg与阶跃前负荷Mg0的偏差相对值,mg0=(Mg-Mg0)/Mg0;
②从波动叠加的角度出发,将一次低阶等效***式(3)变换为两个2阶子***相加的形式
XE(s)=X1(s)+X2(s) (4)
式中,X1(s)、X2(s)分别表示一次低阶等效***XE(s)的两个二阶子***;C1、C2、C3、C4、C5、C6、C7、C8为系数,均可用待定系数法求得;
③对式(4)-式(6)进行拉普拉斯反变换,得到水轮机调节***的一次低阶等效***在负荷阶跃扰动下的机组转速响应波动方程:
xE(t)=x1(t)+x2(t) (7);
式中,xE(t)、x1(t)、x2(t)分别为***XE(s)、X1(s)、X2(s)对应的转速响应波动方程;t为时间变量;
④对完整5阶***式进行二次降阶:x1(t)、x2(t)中有一个对应***的共轭主导复极点,依据极点的分布找出主导极点对应的子波动;
以二阶***来描述***的转速响应特性方程
式中,K2为振幅,δ2为衰减度,ω为角频率,为初相位;
并计算得到调节时间Tp
式中,Δ的取值对大电网为±0.2%、小电网为±0.4%;
⑤利用二阶***进行水轮机调节***转速响应调节品质的影响因素分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410826880.8A CN104408330B (zh) | 2014-12-25 | 2014-12-25 | 水轮机调节***高阶数学模型的降阶方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410826880.8A CN104408330B (zh) | 2014-12-25 | 2014-12-25 | 水轮机调节***高阶数学模型的降阶方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104408330A true CN104408330A (zh) | 2015-03-11 |
CN104408330B CN104408330B (zh) | 2017-06-09 |
Family
ID=52645961
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410826880.8A Active CN104408330B (zh) | 2014-12-25 | 2014-12-25 | 水轮机调节***高阶数学模型的降阶方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104408330B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105159068A (zh) * | 2015-07-23 | 2015-12-16 | 南京信息工程大学 | 一种改进型二阶***的参数设计方法 |
CN105956350A (zh) * | 2016-07-05 | 2016-09-21 | 华中科技大学 | 一种抽水蓄能机组过水***的建模方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN203546687U (zh) * | 2013-11-14 | 2014-04-16 | 武汉大学 | 自适应尾水位变化的水电站尾水调压室 |
CN103870709A (zh) * | 2014-03-31 | 2014-06-18 | 武汉大学 | 一种水泵水轮机全特性曲线的构造方法 |
-
2014
- 2014-12-25 CN CN201410826880.8A patent/CN104408330B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN203546687U (zh) * | 2013-11-14 | 2014-04-16 | 武汉大学 | 自适应尾水位变化的水电站尾水调压室 |
CN103870709A (zh) * | 2014-03-31 | 2014-06-18 | 武汉大学 | 一种水泵水轮机全特性曲线的构造方法 |
Non-Patent Citations (3)
Title |
---|
WENCHENG GUO ET AL.: "Effect Mechanism of Penstock on Stability and Regulation Quality of Turbine Regulating System", 《MATHEMATICAL PROBLEMS IN ENGINEERING》 * |
魏守平等: "水轮机调速器的PID调节规律", 《水利发电学报》 * |
鲍海艳: "水电站调压室设置条件及运行控制研究", 《中国博士学位全文数据库(电子期刊)工程科技II辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105159068A (zh) * | 2015-07-23 | 2015-12-16 | 南京信息工程大学 | 一种改进型二阶***的参数设计方法 |
CN105159068B (zh) * | 2015-07-23 | 2017-08-11 | 南京信息工程大学 | 一种改进型二阶***的参数设计方法 |
CN105956350A (zh) * | 2016-07-05 | 2016-09-21 | 华中科技大学 | 一种抽水蓄能机组过水***的建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104408330B (zh) | 2017-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Guo et al. | Hydropower system operation stability considering the coupling effect of water potential energy in surge tank and power grid | |
CN111027177B (zh) | 一种抽水蓄能机组调频参数优化方法及调频方法 | |
CN106532744A (zh) | 一种基于自抗扰控制的光伏电站抑制电网低频振荡的方法 | |
CN111092441B (zh) | 抽水蓄能电站区域负荷频率控制方法和*** | |
CN109861243B (zh) | 一种基于相位补偿原理的gpss抑制超低频振荡的方法 | |
Guo et al. | Nonlinear modeling and operation stability of variable speed pumped storage power station | |
Huang et al. | Improved transfer function method for flexible simulation of hydraulic-mechanical-electrical transient processes of hydro-power plants | |
CN106372296A (zh) | 一种具有共用管道的水轮机多机微分方程计算方法 | |
EP3156861A1 (en) | Controller for hydroelectric group | |
CN104408330A (zh) | 水轮机调节***高阶数学模型的降阶方法 | |
Naik et al. | IMC tuned PID governor controller for hydro power plant with water hammer effect | |
CN117458534A (zh) | 一种新型液流储能调峰调频方法及装置 | |
CN102354332B (zh) | 一种用于简化柔***直流输电***中rga计算的方法 | |
CN102496945B (zh) | 四阶混沌电力***的无源控制方法 | |
CN103956767A (zh) | 一种考虑尾流效应的风电场并网稳定性分析方法 | |
Kanchanaharuthai | Nonlinear controller design for hydraulic turbine regulating systems via immersion and invariance | |
Feng et al. | Robust Takagi-Sugeno fuzzy control for nonlinear singular time-delay hydraulic turbine governing system | |
CN116088296B (zh) | 抽蓄储能水电机组建模方法、装置及存储介质 | |
CN113437760A (zh) | 安装多直驱风力发电机并联风电场附加阻尼控制器的方法 | |
Yan et al. | Research on governor parameter optimization to suppress ultra-low frequency oscillation of power system caused by hydropower unit | |
Shi et al. | Coordinated control of generator excitation and TCSC based on Hamilton energy function | |
CN111237116A (zh) | 非线性隐极水轮发电机水门与励磁联合调节方法和*** | |
Zhu et al. | Fractional Order Fixed-Time Nonsingular Sliding Mode Control of a Fractional Hydro-Turbine Governing System | |
Nitta et al. | Output Voltage Control of Hydrogen Engine Generator Using A Hill Climbing Method by Adjusting the Hydrogen Flow Rate with A Proportional Control Valve | |
Rittenhouse et al. | Optimal sliding mode Gaussian controller for a hydropower plant |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |