CN106777944B - 一种水电机组经直流送出系统的调速器参数整定方法 - Google Patents
一种水电机组经直流送出系统的调速器参数整定方法 Download PDFInfo
- Publication number
- CN106777944B CN106777944B CN201611116641.9A CN201611116641A CN106777944B CN 106777944 B CN106777944 B CN 106777944B CN 201611116641 A CN201611116641 A CN 201611116641A CN 106777944 B CN106777944 B CN 106777944B
- Authority
- CN
- China
- Prior art keywords
- time constant
- governor
- speed
- value range
- transmission function
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000005540 biological transmission Effects 0.000 claims abstract description 75
- 238000013016 damping Methods 0.000 claims abstract description 56
- 230000001052 transient effect Effects 0.000 claims abstract description 44
- 238000004088 simulation Methods 0.000 claims abstract description 20
- 238000012546 transfer Methods 0.000 claims abstract description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 29
- 230000001133 acceleration Effects 0.000 claims description 23
- 230000033228 biological regulation Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 239000002131 composite material Substances 0.000 abstract 1
- 230000010355 oscillation Effects 0.000 description 11
- 230000000694 effects Effects 0.000 description 7
- 230000003139 buffering effect Effects 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 241000208340 Araliaceae Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Control Of Eletrric Generators (AREA)
- Control Of Water Turbines (AREA)
Abstract
本发明公开了一种水电机组经直流送出系统的调速器参数整定方法，包括：建立水轮机调速器稳定性分析模型；得出发电机及其外部系统的传递函数，水轮机线性化模型的传递函数和调速器线性化模型的传递函数；求出水轮机调速器稳定性分析模型对应的闭环传递函数以及调速系统的闭环特征方程，然后利用劳斯稳定判据得到使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围；在使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围内进行时域仿真，得出水轮机调速器的最优参数。本发明采用了水轮机调速器稳定性分析模型，并结合劳斯稳定判据和时域仿真最终得到水轮机调速器的最优参数，鲁棒性更强。本发明可广泛应用于电力领域。
Description
技术领域
本发明涉及电力领域，尤其是一种水电机组经直流送出系统的调速器参数整定方法。
背景技术
孤岛系统中，水轮机的水锤效应可能会导致调速器不稳定，进而引发低频振荡。而水电机组呈现负阻尼特性是导致低频振荡的主要原因。一般而言，水轮机和调速器在低频段会产生负阻尼，汽轮机本身的相位滞后较小，当调速器的相位滞后也很小时，在低频段会产生正阻尼。这些都和水轮机调速器的参数整定有关。
云南电网与南方电网主网异步联网后，云南电网高水电比例和直流外送负荷高占比的特点，与大型水电厂经直流孤岛送出系统非常相似，外部阻尼的降低和水轮机的水锤效应会在低频段产生较大的相位滞后，造成超低频振荡问题。异步联网的特性实验也证明，如果云南水电机组仍采用同步联网时的参数，将会造成振荡频率为0.05Hz，波动范围在49.9Hz～50.1Hz之间的持续超低频振荡，鲁棒性较弱。
在孤岛系统中，应对水轮机调速器负阻尼特性导致的低频振荡时，一般采取的解决方案是将发电机一次调频死区设置大于直流频率限制器(Frequency LimitController，FLC)死区，但这种方法并不适应于区域电网通过直流来实现的异步互联系统。因此，亟需从源头入手，提出一种水电机组经直流送出系统的调速器参数整定方法，在相应频段进行相位补偿以避免产生超低频振荡。
发明内容
为解决上述技术问题，本发明的目的在于：提供一种鲁棒性强的水电机组经直流送出系统的调速器参数整定方法。
本发明所采取的技术方案是：
一种水电机组经直流送出系统的调速器参数整定方法，其特征在于：包括以下步骤：
建立水轮机调速器稳定性分析模型；
根据水轮机调速器稳定性分析模型得出发电机及其外部系统的传递函数，水轮机线性化模型的传递函数和调速器线性化模型的传递函数；
根据发电机及其外部系统的传递函数，水轮机线性化模型的传递函数和调速器线性化模型的传递函数求出水轮机调速器稳定性分析模型对应的闭环传递函数以及调速系统的闭环特征方程，然后根据调速系统的闭环特征方程利用劳斯稳定判据得到使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围；
在使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围内进行时域仿真，得出水轮机调速器的最优参数。
进一步，所述发电机及其外部系统的传递函数Gd(s)的表达式为：其中，ΔP(s)为机械功率变化量与电磁功率变化量的差值，Δf(s)为频率变化量，TM为发电机的惯性时间常数，D为发电机内部阻尼系数，Kf为外部阻尼系数，s为拉普拉斯算子。
进一步，所述水轮机线性化模型的传递函数Ght(s)的表达式为：其中，Tw为水启动时间，s为拉普拉斯算子。
进一步，所述水轮机调速器为加速-缓冲型调速器，所述加速-缓冲型调速器线性化模型的传递函数Ggov(s)的表达式为：
其中，bp为永态下降率；bt为暂态下降率；Td为缓冲时间常数；Tn为测量加速时间常数，T1v为测量惯性时间常数；Ty为主配压阀和主接力器时间常数，TG为导叶伺服电动机时间常数，s为拉普拉斯算子。
进一步，所述加速度-缓冲型调速器的暂态下降率bt、缓冲时间常数Td、测量加速度时间常数Tn与PID型微机调速器的比例增益KP、积分增益KI、微分增益KD的转换关系为：
进一步，所述根据发电机及其外部系统的传递函数，水轮机线性化模型的传递函数和调速器线性化模型的传递函数求出水轮机调速器稳定性分析模型对应的闭环传递函数以及调速系统的闭环特征方程，然后根据调速系统的闭环特征方程利用劳斯稳定判据得到使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围这一步骤，其包括：
根据发电机及其外部系统的传递函数Gd(s)，水轮机线性化模型的传递函数Ght(s)和调速器线性化模型的传递函数Ggov(s)求出水轮机调速器稳定性分析模型对应的闭环传递函数Gbh(s)，所述闭环传递函数Gbh(s)的表达式为：
根据闭环传递函数Gbh(s)的表达式得到调速系统的闭环特征方程，所述调速系统的闭环特征方程为：1+Ggov(s)Ght(s)Gd(s)＝0；
采用劳斯稳定判据对调速系统的闭环特征方程进行求解，得到使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围，所述劳斯稳定判据为：调速系统的特征方程各项系数符号相同且不为零，且由各项系数组成的劳斯排列阵的第一列所有项均为正。
进一步，所述采用劳斯稳定判据对调速系统的闭环特征方程进行求解，得到使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围这一步骤，其包括：
采用简化的调速器线性化模型的传递函数Ggov(s)对调速系统的闭环特征方程进行化简，得到化简后的闭环特征方程；
采用劳斯稳定判据对化简后的闭环特征方程进行求解，得到使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围。
进一步，所述简化的调速器线性化模型的传递函数Ggov(s)的表达式为：或其中，bp为永态下降率；Tn为测量加速时间常数，TG为导叶伺服电动机时间常数，s为拉普拉斯算子。
进一步，所述在使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围内进行时域仿真，得出水轮机调速器的最优参数这一步骤，其具体为：
若调速器线性化模型的传递函数Ggov(s)的表达式为则在使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围内进行时域仿真，得出水轮机调速器的最优参数范围，所述水轮机调速器的最优参数范围为：Td＝4Tw，其中，s为拉普拉斯算子，Tw为水启动时间，TM为发电机的惯性时间常数；
若调速器线性化模型的传递函数Ggov(s)的表达式为则在使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围内进行时域仿真，得出水轮机调速器的最优参数范围，所述水轮机调速器的最优参数范围为：
Td＝[5-0.5(Tw-1)]Tw，
其中，bp为永态下降率，s为拉普拉斯算子，Tw为水启动时间，TM为发电机的惯性时间常数，TG为导叶伺服电动机时间常数；
若调速器线性化模型的传递函数Ggov(s)的表达式为则在使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围内进行时域仿真，得出水轮机调速器的最优参数范围，所述水轮机调速器的最优参数范围为：
Tn＝(0.3～0.5)Tw，3Tw＜Td<5Tw，其中，s为拉普拉斯算子，Tw为水启动时间，TM为发电机的惯性时间常数，Tn为测量加速时间常数。
进一步，所述水轮机调速器的最优参数为阻尼特性最好且恢复速度最快时对应的暂态下降率取值以及缓冲时间常数取值。
本发明的有益效果是：包括建立水轮机调速器稳定性分析模型，获取水轮机调速器稳定性分析模型的闭环传递函数和调速系统的闭环特征方程，利用劳斯稳定判据得到使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围以及时域仿真得出水轮机调速器的最优参数的过程，采用了水轮机调速器稳定性分析模型，并结合劳斯稳定判据和时域仿真最终得到使调速系统稳定的水轮机调速器最优参数，能够有效避免外部阻尼的降低和水轮机的水锤效应造成的超低频振荡，鲁棒性更强。
附图说明
图1为本发明一种水电机组经直流送出系统的调速器参数整定方法的整体流程图；
图2为本发明的水轮机调速器稳定性分析模型的结构示意图；
图3为本发明的水轮机加速-缓冲型调速器的结构示意图；
图4为本发明的水轮机并联型PID型微机调速器的结构示意图。
具体实施方式
参照图1，一种水电机组经直流送出系统的调速器参数整定方法，其特征在于：包括以下步骤：
建立水轮机调速器稳定性分析模型；
根据水轮机调速器稳定性分析模型得出发电机及其外部系统的传递函数，水轮机线性化模型的传递函数和调速器线性化模型的传递函数；
根据发电机及其外部系统的传递函数，水轮机线性化模型的传递函数和调速器线性化模型的传递函数求出水轮机调速器稳定性分析模型对应的闭环传递函数以及调速系统的闭环特征方程，然后根据调速系统的闭环特征方程利用劳斯稳定判据得到使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围；
在使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围内进行时域仿真，得出水轮机调速器的最优参数。
进一步作为优选的实施方式，所述发电机及其外部系统的传递函数Gd(s)的表达式为：其中，ΔP(s)为机械功率变化量与电磁功率变化量的差值，Δf(s)为频率变化量，TM为发电机的惯性时间常数，D为发电机内部阻尼系数，Kf为外部阻尼系数，s为拉普拉斯算子。
进一步作为优选的实施方式，所述水轮机线性化模型的传递函数Ght(s)的表达式为：其中，Tw为水启动时间，s为拉普拉斯算子。
进一步作为优选的实施方式，所述水轮机调速器为加速-缓冲型调速器，所述加速-缓冲型调速器线性化模型的传递函数Ggov(s)的表达式为：
其中，bp为永态下降率；bt为暂态下降率；Td为缓冲时间常数；Tn为测量加速时间常数，T1v为测量惯性时间常数；Ty为主配压阀和主接力器时间常数，TG为导叶伺服电动机时间常数，s为拉普拉斯算子。
进一步作为优选的实施方式，所述加速度-缓冲型调速器的暂态下降率bt、缓冲时间常数Td、测量加速度时间常数Tn与PID型微机调速器的比例增益KP、积分增益KI、微分增益KD的转换关系为：
进一步作为优选的实施方式，所述根据发电机及其外部系统的传递函数，水轮机线性化模型的传递函数和调速器线性化模型的传递函数求出水轮机调速器稳定性分析模型对应的闭环传递函数以及调速系统的闭环特征方程，然后根据调速系统的闭环特征方程利用劳斯稳定判据得到使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围这一步骤，其包括：
根据发电机及其外部系统的传递函数Gd(s)，水轮机线性化模型的传递函数Ght(s)和调速器线性化模型的传递函数Ggov(s)求出水轮机调速器稳定性分析模型对应的闭环传递函数Gbh(s)，所述闭环传递函数Gbh(s)的表达式为：
根据闭环传递函数Gbh(s)的表达式得到调速系统的闭环特征方程，所述调速系统的闭环特征方程为：1+Ggov(s)Ght(s)Gd(s)＝0；
采用劳斯稳定判据对调速系统的闭环特征方程进行求解，得到使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围，所述劳斯稳定判据为：调速系统的特征方程各项系数符号相同且不为零，且由各项系数组成的劳斯排列阵的第一列所有项均为正。
进一步作为优选的实施方式，所述采用劳斯稳定判据对调速系统的闭环特征方程进行求解，得到使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围这一步骤，其包括：
采用简化的调速器线性化模型的传递函数Ggov(s)对调速系统的闭环特征方程进行化简，得到化简后的闭环特征方程；
采用劳斯稳定判据对化简后的闭环特征方程进行求解，得到使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围。
进一步作为优选的实施方式，所述简化的调速器线性化模型的传递函数Ggov(s)的表达式为：或其中，bp为永态下降率；Tn为测量加速时间常数，TG为导叶伺服电动机时间常数，s为拉普拉斯算子。
进一步作为优选的实施方式，所述在使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围内进行时域仿真，得出水轮机调速器的最优参数这一步骤，其具体为：
若调速器线性化模型的传递函数Ggov(s)的表达式为则在使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围内进行时域仿真，得出水轮机调速器的最优参数范围，所述水轮机调速器的最优参数范围为：Td＝4Tw，其中，s为拉普拉斯算子，Tw为水启动时间，TM为发电机的惯性时间常数；
若调速器线性化模型的传递函数Ggov(s)的表达式为则在使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围内进行时域仿真，得出水轮机调速器的最优参数范围，所述水轮机调速器的最优参数范围为：
Td＝[5-0.5(Tw-1)]Tw，
其中，bp为永态下降率，s为拉普拉斯算子，Tw为水启动时间，TM为发电机的惯性时间常数，TG为导叶伺服电动机时间常数；
若调速器线性化模型的传递函数Ggov(s)的表达式为则在使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围内进行时域仿真，得出水轮机调速器的最优参数范围，所述水轮机调速器的最优参数范围为：
Tn＝(0.3～0.5)Tw，3Tw＜Td<5Tw，其中，s为拉普拉斯算子，Tw为水启动时间，TM为发电机的惯性时间常数，Tn为测量加速时间常数。
进一步作为优选的实施方式，所述水轮机调速器的最优参数为阻尼特性最好且恢复速度最快时对应的暂态下降率取值以及缓冲时间常数取值。
下面结合说明书附图和具体实施例对本发明作进一步详细说明。
实施例一
针对现有技术无法避免外部阻尼的降低和水轮机的水锤效应造成的超低频振荡的问题，本发明提出了一种全新的水电机组经直流送出系统的调速器参数整定方法。该方法通过调整调速器的参数来在相应频段进行相位补偿，以避免产生超低频振荡。
如图1所示，该水电机组经直流送出系统的调速器参数整定方法包括以下步骤：
步骤1，建立如图2所示的水轮机调速器稳定性分析模型。
步骤2，根据图2所示的水轮机调速器稳定性分析模型得出发电机及其外部系统的传递函数Gd(s)，水轮机线性化模型的传递函数Ght(s)，调速器线性化模型的传递函数Ggov(s)。
1)如图2所示，发电机及其外部系统的传递函数Gd(s)为：
其中，ΔP(s)为机械功率变化量与电磁功率变化量的差值，Δf(s)为频率变化量，TM为发电机的惯性时间常数，D为发电机内部阻尼系数，Kf为外部阻尼系数(其包括负荷阻尼、发电机组一次调频作用、直流FLC(即直流频率限制)控制作用)，s为拉普拉斯算子。
在本实施例中，由于异步联网后，云南电网直流外送负荷占比很高，外部阻尼系数Kf极小，可忽略不计。同时，为简化计算，本实施例可忽略数值较小的发电机内部阻尼系数D。故本实施例在计算发电机及其外部系统的传递函数Gd(s)，将不考虑D、Kf。
2)如图2所示，水轮机线性化模型的传递函数Ght(s)为：
其中，Tw为水启动时间(即水锤效应时间常数)，Tw随发电机工况变化而变化，Tw的取值一般在0.5～4.0s之间。
3)如图3所示，忽略死区和饱和等非线性环节以及忽略辅助阀及伺服电动机环节(该环节的时间常数Tp足够小，通常Tp取0.05s)的影响，加速-缓冲型调速器线性化模型的传递函数Ggov(s)为：
其中，bp为永态下降率(即调差系数)；bt为暂态下降率；Td为缓冲时间常数；Tn为测量加速时间常数，T1v为测量惯性时间常数；Tp为辅助阀及伺服电动机时间常数；Ty为主配压阀和主接力器时间常数，Rmax为主配压阀和主接力器的过速开启系数，Rmin为主配压阀和主接力器的过速关闭系数；TG为导叶伺服电动机时间常数。
参照图4，当bp取bp＝0和忽略Ty，T1v的作用时，并联型和串联型PID结构具有相同的传递函数：
而PID型微机调速器模型与加速-缓冲型调速器模型间存在明确的参数转化关系。其中，PID型微机调速器的比例增益KP、积分增益KI、微分增益KD与加速度-缓冲型调速器的暂态下降率bt、缓冲时间常数Td、测量加速度时间常数Tn的转换关系为：
以加速缓冲性型调速器模型的传递函数为例，Ggov(s)可根据需要进行适当的简化，而采用不同化简模型其对应的最优整定参数不同。本实施例采用的简化的传递函数Ggov(s)为：
步骤3，根据发电机及其外部系统的传递函数Gd(s)，水轮机线性化模型的传递函数Ght(s)和调速器线性化模型的传递函数Ggov(s)求出水轮机调速器稳定性分析模型对应的闭环传递函数Gbh(s)，根据闭环传递函数Gbh(s)的表达式得到调速系统的闭环特征方程，然后利用劳斯稳定判据调速系统的闭环特征方程进行求解，得到使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围。
此过程可进一步细分为：
1)求出水轮机调速器稳定性分析模型对应的闭环传递函数Gbh(s)：
2)令Gbh(s)分母中的多项式为零，得到调速系统的闭环特征方程：
1+Ggov(s)Ght(s)Gd(s)＝0 (8)
在本实施例中，将简化的传递函数代入上式(8)得到的对应闭环特征方程为：
0.5btTdTMTws3+(btTdTM-TwTd)s2+(Td-Tw)s+1＝0 (9)
其中，TM和Tw均为已知量，因此为了对该闭环特征方程进行求解，还需要确定的参数为暂态下降率bt取值范围以及缓冲时间常数Td取值范围。
3)利用劳斯稳定判据(即调速系统的特征方程各项系数符号相同且不为零，且由各项系数组成的劳斯排列阵的第一列所有项均为正)对化简后的闭环特征方程式(9)进行求解，得到使调速系统稳定的bt、Td取值范围，如下式(10)和式(11)所示：
(btTdTM-TwTd)(Td-Tw)＞0.5btTdTMTw (10)
步骤4，根据时域仿真在使调速系统稳定的bt、Td取值范围内进行扫描，确定出轮机调速器的最优参数。
本发明在使调速系统稳定的bt、Td取值范围内，通过大量的时域仿真，最终确定出阻尼特性好且恢复速度快的最优参数，最终完成对水轮机组的参数整定过程。
本实施例中采用的简化的传递函数Ggov(s)为公式(6)，其对应的最优参数整定方案为：
同理，可采用本实施例的方法确定采用其他调速器模型所对应的最优参数整定方案：
a.若采用的调速器简化模型为：则通过采用本实施例的方法对水轮机调速系统的特征值进行分析和大量的时域仿真，可得出其所对应的最优参数整定方案为：
b.若采用的调速器简化模型为：通过采用本实施例的方法对水轮机调速系统的特征值进行分析和大量的时域仿真，可得出其所对应的最优参数整定方案
为：
本发明提出了一种水电机组经直流送出系统的调速器参数整定方法，该方法采用了水轮机调速器稳定性分析模型，并结合劳斯稳定判据和时域仿真最终得到使调速系统稳定的水轮机调速器最优参数，能够有效避免外部阻尼的降低和水轮机的水锤效应造成的超低频振荡，鲁棒性更强，对于合理设置大规模水电经直流送出系统的调频策略，优化异步联网系统安全防线的设计有着重要的意义。
以上是对本发明的较佳实施进行了具体说明，但本发明并不限于所述实施例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替换，这些等同的变形或替换均包含在本申请权利要求所限定的范围内。
Claims (7)
1.一种水电机组经直流送出系统的调速器参数整定方法，其特征在于：包括以下步骤：
建立水轮机调速器稳定性分析模型；
根据水轮机调速器稳定性分析模型得出发电机及其外部系统的传递函数，水轮机线性化模型的传递函数和调速器线性化模型的传递函数；
根据发电机及其外部系统的传递函数，水轮机线性化模型的传递函数和调速器线性化模型的传递函数求出水轮机调速器稳定性分析模型对应的闭环传递函数以及调速系统的闭环特征方程，然后根据调速系统的闭环特征方程利用劳斯稳定判据得到使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围；
在使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围内进行时域仿真，得出水轮机调速器的最优参数；
所述发电机及其外部系统的传递函数Gd(s)的表达式为：其中，ΔP(s)为机械功率变化量与电磁功率变化量的差值，Δf(s)为频率变化量，TM为发电机的惯性时间常数，D为发电机内部阻尼系数，Kf为外部阻尼系数，s为拉普拉斯算子；
所述水轮机线性化模型的传递函数Ght(s)的表达式为：其中，Tw为水启动时间，s为拉普拉斯算子；
所述水轮机调速器为加速度-缓冲型调速器，所述加速度-缓冲型调速器线性化模型的传递函数Ggov(s)的表达式为：
其中，bp为永态下降率；bt为暂态下降率；Td为缓冲时间常数；Tn为测量加速时间常数，T1v为测量惯性时间常数；Ty为主配压阀和主接力器时间常数，TG为导叶伺服电动机时间常数，s为拉普拉斯算子。
2.根据权利要求1所述的一种水电机组经直流送出系统的调速器参数整定方法，其特征在于：所述加速度-缓冲型调速器的暂态下降率bt、缓冲时间常数Td、测量加速度时间常数Tn与PID型微机调速器的比例增益KP、积分增益KI、微分增益KD的转换关系为：
3.根据权利要求1所述的一种水电机组经直流送出系统的调速器参数整定方法，其特征在于：所述根据发电机及其外部系统的传递函数，水轮机线性化模型的传递函数和调速器线性化模型的传递函数求出水轮机调速器稳定性分析模型对应的闭环传递函数以及调速系统的闭环特征方程，然后根据调速系统的闭环特征方程利用劳斯稳定判据得到使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围这一步骤，其包括：
根据发电机及其外部系统的传递函数Gd(s)，水轮机线性化模型的传递函数Ght(s)和调速器线性化模型的传递函数Ggov(s)求出水轮机调速器稳定性分析模型对应的闭环传递函数Gbh(s)，所述闭环传递函数Gbh(s)的表达式为：
根据闭环传递函数Gbh(s)的表达式得到调速系统的闭环特征方程，所述调速系统的闭环特征方程为：1+Ggov(s)Ght(s)Gd(s)＝0；
采用劳斯稳定判据对调速系统的闭环特征方程进行求解，得到使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围，所述劳斯稳定判据为：调速系统的特征方程各项系数符号相同且不为零，且由各项系数组成的劳斯排列阵的第一列所有项均为正。
4.根据权利要求3所述的一种水电机组经直流送出系统的调速器参数整定方法，其特征在于：所述采用劳斯稳定判据对调速系统的闭环特征方程进行求解，得到使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围这一步骤，其包括：
采用简化的调速器线性化模型的传递函数Ggov(s)对调速系统的闭环特征方程进行化简，得到化简后的闭环特征方程；
采用劳斯稳定判据对化简后的闭环特征方程进行求解，得到使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围。
5.根据权利要求4所述的一种水电机组经直流送出系统的调速器参数整定方法，其特征在于：所述简化的调速器线性化模型的传递函数Ggov(s)的表达式为： 或其中，bp为永态下降率；Tn为测量加速时间常数，TG为导叶伺服电动机时间常数，s为拉普拉斯算子。
6.根据权利要求1所述的一种水电机组经直流送出系统的调速器参数整定方法，其特征在于：所述在使调速系统稳定的暂态下降率取值范围以及缓冲时间常数取值范围内进行时域仿真，得出水轮机调速器的最优参数这一步骤，其具体为：
若调速器线性化模型的传递函数Ggov(s)的表达式为则在使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围内进行时域仿真，得出水轮机调速器的最优参数范围，所述水轮机调速器的最优参数范围为：Td＝4Tw，其中，s为拉普拉斯算子，Tw为水启动时间，TM为发电机的惯性时间常数；
若调速器线性化模型的传递函数Ggov(s)的表达式为则在使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围内进行时域仿真，得出水轮机调速器的最优参数范围，所述水轮机调速器的最优参数范围为：
其中，bp为永态下降率，s为拉普拉斯算子，Tw为水启动时间，TM为发电机的惯性时间常数，TG为导叶伺服电动机时间常数；
若调速器线性化模型的传递函数Ggov(s)的表达式为则在使调速系统稳定的暂态下降率bt取值范围以及缓冲时间常数Td取值范围内进行时域仿真，得出水轮机调速器的最优参数范围，所述水轮机调速器的最优参数范围为：
其中，s为拉普拉斯算子，Tw为水启动时间，TM为发电机的惯性时间常数，Tn为测量加速时间常数。
7.根据权利要求1-6任一项所述的一种水电机组经直流送出系统的调速器参数整定方法，其特征在于：所述水轮机调速器的最优参数为阻尼特性最好且恢复速度最快时对应的暂态下降率取值以及缓冲时间常数取值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611116641.9A CN106777944B (zh) | 2016-12-07 | 2016-12-07 | 一种水电机组经直流送出系统的调速器参数整定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611116641.9A CN106777944B (zh) | 2016-12-07 | 2016-12-07 | 一种水电机组经直流送出系统的调速器参数整定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106777944A CN106777944A (zh) | 2017-05-31 |
| CN106777944B true CN106777944B (zh) | 2019-06-18 |
Family
ID=58882357
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611116641.9A Active CN106777944B (zh) | 2016-12-07 | 2016-12-07 | 一种水电机组经直流送出系统的调速器参数整定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106777944B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107168101B (zh) * | 2017-06-07 | 2020-08-14 | 国网福建省电力有限公司 | 考虑调频及稳定约束的机组调速系统控制参数整定方法 |
| CN108107720B (zh) * | 2017-11-23 | 2021-03-16 | 中国南方电网有限责任公司 | 基于状态空间分析的水轮机调速器参数整定方法及系统 |
| CN109936319B (zh) * | 2017-12-18 | 2020-11-06 | 青岛海尔智能技术研发有限公司 | 一种整定转速控制器参数的方法及装置 |
| CN108616122B (zh) * | 2018-05-15 | 2021-06-04 | 电子科技大学 | 一种计及超低频振荡的调速器pid参数鲁棒优化方法 |
| CN109634116B (zh) * | 2018-09-04 | 2022-03-15 | 贵州大学 | 一种分数阶机械式离心调速器系统的加速自适应稳定方法 |
| CN110489824B (zh) * | 2019-07-30 | 2023-06-02 | 中国电建集团华东勘测设计研究院有限公司 | 基于时间域分析法的水轮机调节参数优化选择方法 |
| CN111193259B (zh) * | 2020-01-10 | 2023-09-22 | 燕山大学 | 一种电力电子变压器的电力系统稳定性检测方法及系统 |
| CN112985530B (zh) * | 2021-02-01 | 2022-04-22 | 南京航空航天大学 | 一种基于特征方程根轨迹的燃油计量装置设计参数调整方法 |
| CN117331305B (zh) * | 2023-11-30 | 2024-02-20 | 华中科技大学 | 水轮机调节系统控制参数与振荡特性关联方法及系统 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104808705A (zh) * | 2015-04-27 | 2015-07-29 | 贵州电力试验研究院 | 一种基于特征参数的水电机组调速系统控制参数整定方法 |
-
2016
- 2016-12-07 CN CN201611116641.9A patent/CN106777944B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104808705A (zh) * | 2015-04-27 | 2015-07-29 | 贵州电力试验研究院 | 一种基于特征参数的水电机组调速系统控制参数整定方法 |
Non-Patent Citations (1)
| Title |
|---|
| 水轮发电机及其调速系统的参数辨识方法与控制策略研究;寇攀高;《万方学位论文数据库》;20121225;第2012年卷;第1.1节，第2.1.1.1节，第2.1.3.1节，第6.2.2.1节，第6.2.3节 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106777944A (zh) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106777944B (zh) | 一种水电机组经直流送出系统的调速器参数整定方法 | |
| CN104503228B (zh) | 一种水轮机调速器功率模式下一次调频稳定域的确定方法 | |
| CN107800146B (zh) | 兼顾一次调频和超低频振荡抑制的调速器参数优化方法 | |
| CN107679769B (zh) | 含风电的电力系统频率响应模型建立方法和频率特性指标计算方法 | |
| CN107453375A (zh) | 巨型水电机组功率模式下的一次调频电量补偿精细化控制方法及装置 | |
| CN104090491B (zh) | 燃气蒸汽联合循环机组多变量受限预测函数负荷控制方法 | |
| CN105826949B (zh) | 一种提高虚拟同步控制逆变器同步稳定性的暂态控制方法 | |
| TW201316646A (zh) | 一種整合自動電壓調整器之控制系統和方法 | |
| CN106505591B (zh) | 一种多机电力系统超低频振荡分析方法 | |
| CN110518631B (zh) | 一种直驱风电机组的稳定性评估方法及系统 | |
| CN106712055A (zh) | 一种与低励限制功能相协调的电力系统稳定器配置方法 | |
| CN105512367B (zh) | 并网火电机组一次调频转速不等率的临界稳定值确定方法 | |
| CN107133749A (zh) | 一种考虑需求响应信息的电力信息物理耦合建模方法 | |
| CN101718248B (zh) | 水轮发电机组主动稳定控制方法 | |
| CN105186511B (zh) | 电池储能系统参与电网二次调频控制方法 | |
| Salhi et al. | Frequency regulation for large load variations on micro-hydro power plants with real-time implementation | |
| CN110649596B (zh) | 一种考虑系统初始状态的频率全响应解析模型 | |
| CN109103927A (zh) | 提高一次调频动态响应特性的调速系统pid控制器参数整定方法 | |
| CN109412180A (zh) | 水电机组异步联网场景下的超低频振荡控制方法 | |
| CN101915203A (zh) | 改善水轮发电机组功角振荡的阻尼注入控制方法 | |
| WO2018145498A1 (zh) | 基于强化学习算法的双馈感应风力发电机自校正控制方法 | |
| CN108809189A (zh) | 一种并网调频型飞轮储能系统充电过程中的转速控制方法 | |
| CN102629293A (zh) | 水力机组附加机械阻尼系数的模拟计算方法 | |
| CN108280245A (zh) | 一种水轮机调速系统的阻尼特性分析方法 | |
| CN106647240B (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 |