CN109185036B - 基于频率响应的水轮机pid型调速器参数设计方法 - Google Patents

基于频率响应的水轮机pid型调速器参数设计方法 Download PDF

Info

Publication number
CN109185036B
CN109185036B CN201811046089.XA CN201811046089A CN109185036B CN 109185036 B CN109185036 B CN 109185036B CN 201811046089 A CN201811046089 A CN 201811046089A CN 109185036 B CN109185036 B CN 109185036B
Authority
CN
China
Prior art keywords
solution space
speed regulator
frequency
parameter
parameters
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
Application number
CN201811046089.XA
Other languages
English (en)
Other versions
CN109185036A (zh
Inventor
黄润
黄伟
吴琛
张�杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yunnan Power Grid Co Ltd
Original Assignee
Yunnan Power Grid Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yunnan Power Grid Co Ltd filed Critical Yunnan Power Grid Co Ltd
Priority to CN201811046089.XA priority Critical patent/CN109185036B/zh
Publication of CN109185036A publication Critical patent/CN109185036A/zh
Application granted granted Critical
Publication of CN109185036B publication Critical patent/CN109185036B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for AC mains or AC distribution networks
    • H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for AC mains or AC distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F05—INDEXING SCHEMES RELATING TO ENGINES OR PUMPS IN VARIOUS SUBCLASSES OF CLASSES F01-F04
    • F05B—INDEXING SCHEME RELATING TO WIND, SPRING, WEIGHT, INERTIA OR LIKE MOTORS, TO MACHINES OR ENGINES FOR LIQUIDS COVERED BY SUBCLASSES F03B, F03D AND F03G
    • F05B2270/00—Control
    • F05B2270/70—Type of control algorithm
    • F05B2270/706—Type of control algorithm proportional-integral-differential
    • 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)
  • Power Engineering (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明涉及基于频率响应的水轮机PID型调速器参数设计方法，该方法包括：建立含有系统主要发电机组的统一频率模型；选择调速器参数设计目标机组；扫频测量得到各个运行方式下统一频率系统的开环频率响应数据；求解各个运行方式下考虑系统稳定约束的参数解空间S 1；求解各个运行方式下考虑调节性能约束的参数解空间S 2；求解各个运行方式下同时考虑系统稳定与调节性能的参数解空间S；求解多运行方式下同时考虑系统稳定与调节性能的参数解空间S m ；基于参数解空间S m 设计调速器参数。本发明提出的设计方法简单易实现，可应用于多机电力系统考虑多运行方式下的水轮机调速器参数设计，使系统安全稳定运行。

Description

基于频率响应的水轮机PID型调速器参数设计方法
技术领域
本发明属于电网安全技术领域，涉及一种水轮机PID型调速器参数设计方法，具体涉及一种基于电力系统频率响应的多运行方式下多机电力系统水轮机PID型调速器参数设计方法。
背景技术
交流系统超低频振荡现象在国内、外实际电网运行中均有发现，具有持续时间长、振荡频率极低(小于0.1Hz)的特点，威胁系统安全稳定运行。自从云南电网与南方电网主网异步联网后，云南电网也出现了超低频振荡现象。已有研究对超低频振荡问题的机理做了较为深入的剖析，发现水电机组在超低频段内无法为系统提供足够的阻尼转矩，超低频振荡与系统中的水电机组一次调频关系密切。抑制超低频振荡最简单且有效的手段是切除水电机组调速器，然而切除调速器会降低机组跟踪负荷波动的性能，只能作为应急手段。
调速器参数设计需综合考虑系统稳定性与调节性能，且应能应对多运行方式下电力系统的稳定需求。如何克服现有技术的不足是目前电网安全技术领域亟需解决的问题。
发明内容
本发明的目的是为了解决现有技术的不足，提供一种抑制超低频振荡问题的水轮机调速器参数设计方法，该方法兼顾了系统的稳定性要求与调节性能要求，不需要推导系统的传递函数，而仅依靠测量得到的系统频率响应数据。
为实现上述目的，本发明采用的技术方案如下：
1、基于频率响应的水轮机PID型调速器参数设计方法，包括以下步骤：
1)建立含有系统主要发电机组在内的统一频率模型；
2)选择调速器参数设计的目标电厂机组；
3)扫频测量得到各个运行方式下统一频率模型的开环频率响应数据；
4)根据3)得到的统一频率模型的开环频率响应数据，求解各个运行方式下考虑系统稳定约束的参数解空间S1；
5)根据3)得到的统一频率模型的开环频率响应数据，求解各个运行方式下考虑调节性能约束的参数解空间S2；
6)根据4)得到参数解空间S1的和5)得到的参数解空间S2，求解各个运行方式下同时考虑系统稳定与调节性能的参数解空间S；
7)根据6)得到参数解空间S，求解多运行方式下同时考虑系统稳定与调节性能的参数解空间Sm；
8)基于参数解空间Sm设计调速器参数。
进一步，优选的是，所述的步骤1)中的建立含有系统主要发电机组的统一频率模型步骤：
(1)统计系统所含发电机组调速器和原动机模型的类型以及它们各自占系统总容量的比例；
(2)对超低频振荡发生时各电厂机组有功功率振幅从大到小排列，选择排列在第一个或前几个的电厂机组单独建模；将各类型调速器和原动机模型按照其各自占系统总容量的比例从大到小排序，选择其中占比最大的1～4种模型搭建于统一频率模型中；
(3)对于(2)中选定的调速器和原动机模型，将其参数设置为系统中出力最大的该类型机组的参数；同时按照以下方式设置统一频率模型中各类型机组的台数：统计各类型机组在系统中的总有功出力，则统一频率模型中各类型机组的台数等于其在系统中的总有功出力除以它的单台有功出力；
(4)根据系统自动发电控制的安排设置统一频率模型中参与自动发电控制的机组台数：将统一频率模型中参与自动发电控制的机组台数设置为系统中参与自动发电控制的机组总有功出力除以单台机组的有功出力的值。
进一步，优选的是，所述的步骤2)中的选择调速器参数设计的目标电厂机组步骤：
(1)计算所有单独建模电厂机组的阻尼系数MD；
对第i台机组，其阻尼系数MDi为：
其中，Re()为取实部运算，GTi(s)为电厂机组调速器-原动机传递函数，s为拉普拉斯算子，j为虚数单位，ωc为振荡频率，单位为rad/s；
(2)对于阻尼系数MD从小到大排列，选择排列在第一个或前几个的电厂机组为调速器参数设计目标电厂机组。
进一步，优选的是，所述的步骤3)中的扫频测量得到各个运行方式下统一频率模型的开环频率响应数据步骤：
在统一频率模型中，将除调速器参数设计目标电厂机组调速器以外的系统统一为开环系统P(s)；
通过扫频测量得到各个运行方式下统一频率模型的开环频率响应数据，即P(s)的频率响应数据，其中，P(s)的频率响应数据包括P(s)的幅频特性数据|P(jω)|与相频特性数据∠P(jω)。
进一步，优选的是，所述的步骤4)中的求解各个运行方式下考虑系统稳定约束的参数解空间S1步骤：
(1)目标电厂调速器与开环系统P(s)构成的反馈系统闭环稳定的等价条件为：回差矩阵行列式F(s)＝1+(KP·s+KI+KD·s2)/s·P(s)不存在右半平面的零点，其中，KP、KI、KD为PID调速器待设计参数；
(2)定义开环系统P(s)的相移数为σ＝(∠P(j∞)-∠P(-j∞))/(π/2)；
(3)定义开环系统P(s)的相对次数为
(4)将KP参数离散采样得到[KP1,KP2,…,KPn]；定义如下函数：
对于每一个KP采样点，根据上式函数曲线求解的奇数重根集为{ω1,ω2,...,ωl-1}；
(5)定义并求解如下序列：
I＝{i1,i2,...,il-1}
其中，it∈{+1,-1},t＝1,2,...,l-1；
如果rP-1是偶数，求解如下序列：
[2i1-2i2+2i3+…+(-1)l-12il-1]·Sx＝-σ-2
如果rP-1是奇数，求解如下序列：
[2i1-2i2+2i3+…+(-1)l-12il-1]·Sx＝-σ-2
(6)对于求得的第k条序列Ik＝{i1,k,i2,k,...,il-1,k}，求得如下参数解空间：
由上式可知每条序列Ik对应的解空间是由线性约束构成的解空间，记为再将所有序列的解空间求并集就得到考虑系统稳定的解空间其中U为求并集运算。
进一步，优选的是，所述的步骤5)中的求解各个运行方式下考虑调节性能约束的参数解空间S2步骤：
(1)电力系统调节性能的约束条件描述如下：
||w(s)S(s)||∞<1
其中：w(s)为性能权函数，S(s)为系统的灵敏度函数，其表达式如下：
式中：s为拉普拉斯算子，M为灵敏度函数峰值上界，ωB *为带宽频率下界，A为最大稳态跟踪误差；
该约束条件等价于目标电厂调速器与如下定义的复系数传递函数Pc(s,θr)构成的反馈系统闭环稳定：
对θ进行采样得到样本集[θ1θ2…θm]，即对于每一个样本θr，目标电厂调速器与Pc(s,θr)构成的反馈系统的回差矩阵行列式F(s)＝1+(KP·s+KI+KD·s2)/s·Pc(s,θr)不存在右半平面的零点；
(2)定义Pc(s,θr)的相移数为σ＝(∠Pc(j∞,θr)-∠Pc(-j∞,θr))/(π/2)；
(3)定义Pc(s,θr)的相对次数为
(4)将KP参数离散采样得到[KP1,KP2,…,KPn]；定义如下函数：
对于每一个KP采样点，根据上式函数曲线求解的奇数重根集为{ω1,ω2,...,ωl-1}；
(5)定义并求解如下序列：
I′＝{i′1,i′2,...,i′l-1}
其中，i′t∈{+1,-1},t＝1,2,...,l-1；
求解如下序列：
[2i′1-2i′2+2i′3+…+(-1)l-12i′l-1]]·Sx＝-σ-2
(6)对于求得的第k条序列I′k＝{i′1,k,i′2,k,...,i′l-1,k}可求得如下参数解空间：
由上式可知每条序列I′k对应的解空间是由线性约束构成的解空间，记为再将所有序列的解空间求并集就得到θr对应的解空间其中∪为求并集运算；
(7)将根据(6)求得的所有θr对应的解空间求交集就得到考虑系统调节性能约束的解空间其中∩为求交集运算。
进一步，优选的是，所述的步骤6)中的求解各个运行方式下同时考虑系统稳定与调节性能的参数解空间S＝S1∩S2；
进一步，优选的是，所述的步骤7)中的求解多运行方式下同时考虑系统稳定与调节性能的参数解空间Sm为各个运行方式下解空间S的交集。
进一步，优选的是，所述的步骤8)中的基于参数解空间Sm设计调速器参数步骤：
根据得到的解空间Sm，设置目标函数为max KP·KI，进行求解，所得到的KP、KI、KD为最终的水轮机调速器设计参数。
进一步，优选的是，利用matlab自带的函数fmincon求解max KP·KI。
本发明中M为灵敏度函数峰值上界，根据稳定要求设置；ωB *为带宽频率下界，根据跟踪性能要求设置；A为最大稳态跟踪误差，设置为很小值，如0.0001。
本发明与现有技术相比，其有益效果为：
本发明提出的方法可以应对多运行方式下电力系统水轮机PID型调速器参数设计问题，参数设计过程不需要推导系统的传递函数，而仅依靠扫频测量得到的统一频率模型的开环频率响应数据，能够给出所有满足系统的稳定与调节性能约束的线性参数解空间，再通过给定的参数目标函数max KP·KI，将水轮机调速器参数设计问题转化为一个简单易求解的非线性规划问题，设计调速器参数。该方法是PID型调速器控制参数设计的有力工具，且从未应用于解决超低频振荡问题。
附图说明
图1是本发明实施例流程图；
图2是本发明实施例提供的统一频率模型；
图3是本发明实施例提供的仿真验证中云南电网统一频率模型；
图4是本发明实施例提供的仿真验证中GS型汽轮机调速器模型结构图；
图5是本发明实施例提供的仿真验证中TB型汽轮机模型结构图；
图6是本发明实施例提供的仿真验证中GM型水轮机调速器模型结构图；
图7是本发明实施例提供的仿真验证中GA型电液伺服系统模型结构图；
图8是本发明实施例提供的仿真验证中水轮机详细模型结构图；
图9是本发明实施例提供的仿真验证中GH型调速器和TW型水轮机模型结构图；
图10是本发明实施例提供的云南电网冬大运行方式下LDL机组调速器除外的开环系统P(s)的频率响应曲线；
图11是本发明实施例提供的云南电网冬大运行方式下考虑系统稳定约束的参数解空间S1；
图12是本发明实施例提供的云南电网冬大运行方式下考虑系统调节性能约束的参数解空间S2；
图13是本发明实施例提供的云南电网冬大运行方式下同时考虑系统稳定与调节性能约束的参数解空间S；
图14是本发明实施例提供的7个方式下同时考虑系统稳定与调节性能的参数解空间Sm；
图15是本发明实施例提供的设计参数与原始参数的时域仿真结果对比。
具体实施方式
下面结合实施例对本发明作进一步的详细描述。
本领域技术人员将会理解，下列实施例仅用于说明本发明，而不应视为限定本发明的范围。实施例中未注明具体技术或条件者，按照本领域内的文献所描述的技术或条件或者按照产品说明书进行。所用材料或设备未注明生产厂商者，均为可以通过购买获得的常规产品。
在本发明实施例中，提供了基于频率响应的水轮机PID型调速器参数设计方法，如图1所示，该方法包括：
步骤101：建立含有系统主要发电机组在内的统一频率模型；
步骤102：选择调速器参数设计的目标电厂机组；
步骤103：扫频测量得到各个运行方式下统一频率系统的开环频率响应数据；
步骤104：求解各个运行方式下考虑系统稳定约束的参数解空间S1；
步骤105：求解各个运行方式下考虑调节性能约束的参数解空间S2；
步骤106：求解各个运行方式下同时考虑系统稳定与调节性能的参数解空间S＝S1∩S2；
步骤107：求解多运行方式下同时考虑系统稳定与调节性能的参数解空间Sm为各个方式下解空间S的交集；
步骤108：基于参数解空间Sm设计调速器参数。
由图1所示的流程可知，在本发明实施例中，通过建立含有系统主要发电机组调速器和原动机模型的统一频率模型，可以实现对原系统超低频振荡的模拟；通过扫频测量得到各个运行方式下统一频率系统的开环频率响应数据，作为参数设计的基础；通过求解各个运行方式下考虑系统稳定约束的参数解空间S1，获得满足稳定要求的PID调速器参数空间；通过求解各个运行方式下考虑调节性能约束的参数解空间S2，获得满足调节性能要求的PID调速器参数空间；通过将S1与S2求交集求得各个运行方式下同时考虑系统稳定与调节性能的参数解空间S；通过将各个运行方式下的解空间S求交集，获得多运行方式下同时考虑系统稳定与调节性能的参数解空间Sm；基于参数解空间Sm，设计得到调速器参数。
(1)具体实施时，统计南方电网2017冬大、冬小、夏大、夏大极限、夏小、夏小外送、汛期方式中云南电网的机组数据和潮流计算结果的基础上，建立统一频率模型。按照之前所述步骤，在Matlab/Simulink软件中建立统一频率模型如图2所示。
以冬大方式为例，上述的统一频率模型建立遵循如下步骤：
1)方式数据中火电出力为3977.2MW，水电出力为32559MW。
2)根据实际监测数据，电厂机组(NZD，LDL，XW，DCS)在振荡发生期间有功功率振幅从大到小排列排在最前，故单独建模，，采用GM型调速器模块、GA型伺服器模块以及水轮机详细模型模块。对于火电机组，典型机组类型为GS(包含GS型调速器模块与TB型汽轮机模块)。对于水电机组，典型机组类型为GH(GH模型内本身包含调速器与原动机)和GM(包含GM型调速器模块、GA型伺服器模块以及水轮机详细模型模块)。单独建模电厂机组总出力为7885MW，GS型机组总出力为3977.2MW，GH型机组总出力为7717MW，GM型机组总出力为16957MW。
3)对于2)中选定的调速器和原动机模型，将其参数设置为全系统中出力最大的该类型机组的参数，具体如下：
对于单独建模电厂机组，NZD电厂机组功率基值为722.3MW，在其自身功率基值下机组惯性为10.38。其所含GM型调速器模型如图4所示，模式选择为开度模式、硬反馈输入信号为YPID，参数如表1。
表1
TR1 0.02 TD 0.2
KW 2.9 T1 0
KP 2 T2 0
KD 1 TR2 0.02
KI 0.37 eP 0.04
其所含GA型伺服模型如图5所示，参数如表2。
表2
Tc 11.4 KD 0
To 18.5 KI 0
T2 0.02 Td 0
KP 10
其所含水轮机详细模型如图6所示，图中水系统传递函数F(s)为
其中，
其中，Tep为压力水管弹性时间常数，Ts为调压室时间常数，TWc为引水洞水击时间常数，TWp为压力水管水击时间常数，Zp为压力水管水阻抗，φc为引水洞摩擦系数，φp为压力水管摩擦系数，s为拉普拉斯算子。参数如表3。
表3
LDL电厂机组功率基值为400MW，在其自身功率基值下机组惯性为9.61。其所含GM型调速器模型如图4所示，模式选择为开度模式、硬反馈输入信号为YPID，参数如表4。
表4
其所含GA型伺服模型如图5所示，参数如表5。
表5
Tc 9.35 KD 0
To 23 KI 0
T2 0.2 Td 0
KP 15
其所含水轮机详细模型如图6所示，图中水系统传递函数F(s)为
其中，
其中，Tep为压力水管弹性时间常数，Ts为调压室时间常数，TWc为引水洞水击时间常数，TWp为压力水管水击时间常数，Zp为压力水管水阻抗，φc为引水洞摩擦系数，φp为压力水管摩擦系数，s为拉普拉斯算子。参数如如表6。
表6
XW电厂机组功率基值为778MW，在其自身功率基值下机组惯性为8.72。其所含GM型调速器模型如图4所示，模式选择为开度模式、硬反馈输入信号为YPID，参数如表7。
表7
TR1 0.02 TD 0.21
KW 1.35 T1 0
KP 2 T2 0
KD 1 TR2 0.02
KI 0.375 eP 0.04
其所含GA型伺服模型如图5所示，参数如表8。
表8
Tc 10 KD 0
To 13.1 KI 0
T2 0.02 Td 0
KP 6
其所含水轮机详细模型如图6所示，图中水系统传递函数F(s)为
其中，
其中，Tep为压力水管弹性时间常数，Ts为调压室时间常数，TWc为引水洞水击时间常数，TWp为压力水管水击时间常数，Zp为压力水管水阻抗，φc为引水洞摩擦系数，φp为压力水管摩擦系数，s为拉普拉斯算子。参数如表9。
表9
DCS电厂机组功率基值为250MW，在其自身功率基值下机组惯性为9.63。其所含GM型调速器模型如图4所示，模式选择为开度模式、硬反馈输入信号为YPID，参数如表10。
表10
TR1 0.2 TD 0.3
KW 1.8 T1 0
KP 4 T2 0
KD 0 TR2 0.02
KI 0.048 eP 0.04
其所含GA型伺服模型如图5所示，参数如表11。
表11
Tc 21.5 KD 0
To 22.5 KI 0
T2 0.2 Td 0
KP 40
其所含水轮机详细模型如图6所示，图中水系统传递函数F(s)为
其中，
其中，Tep为压力水管弹性时间常数，Ts为调压室时间常数，TWc为引水洞水击时间常数，TWp为压力水管水击时间常数，Zp为压力水管水阻抗，φc为引水洞摩擦系数，φp为压力水管摩擦系数，s为拉普拉斯算子。参数如表12。
表12
对于典型机组GS，功率基值为353MW，在其自身功率基值下机组惯性为8.5。其所含GS型调速器模型如图7所示，TB型汽轮机模型如图8所示，参数如表13。
表13
K 18.47 TRH 8
T1 0 FIP 0.4
T2 0 TCO 0.4
T3 0.5 FLP 0.3
TCH 0.3 λ 0
FHP 0.3
对于典型机组GH，功率基值为109MW，在其自身功率基值下机组惯性为8.31。其所含GH型机组模型如图9所示，参数如表14。
表14
R 0.04 TW 1
TG 0.25 Dd 0.25
TP 0.04 K’ 1
TD 5
对于典型机组GM，功率基值为667MW，在其自身功率基值下机组惯性为10.38。其所含GM型调速器模型如图4所示，模式选择为开度模式、硬反馈输入信号为YPID，参数如表15。
表15
TR1 0.2 TD 0.21
KW 1.7 T1 0
KP 2.5 T2 0
KD 1 TR2 0.02
KI 0.5 eP 0.04
其所含GA型伺服模型如图5所示，参数如表16。
表16
Tc 13 KD 0
To 13 KI 0
T2 0.2 Td 0
KP 20
其所含水轮机详细模型如图6所示，图中水系统传递函数F(s)为
其中，
其中，Tep为压力水管弹性时间常数，Ts为调压室时间常数，TWc为引水洞水击时间常数，TWp为压力水管水击时间常数，Zp为压力水管水阻抗，φc为引水洞摩擦系数，φp为压力水管摩擦系数，s为拉普拉斯算子。参数如表17。
表17
冬大方式下，单独建模电厂(NZD，LDL，XW，DCS)机组台数分别为6,3,4,5台，总出力为7885MW。由于GS型典型机组出力为240MW，总出力为3977.2MW，因此GS型典型机组台数为3977.2/240≈17台。由于GH型典型机组出力为95MW，总出力为7717MW，因此GH型典型机组台数为7717/95≈81台。由于GM型典型机组出力为600MW，总出力为24842-7885＝16957MW，因此GM型典型机组台数为16957/600≈28台。
4)AGC增益KAGC＝1。参与AGC的机组类型为GH型和GM型，其台数NGH-AGC为20台，NGM-AGC为2台。
故最终冬大方式统一频率模型中台数参数如表18。
表18
N<sub>GS</sub> 17 N<sub>AGC-GH</sub> 20
N<sub>GH</sub> 81 N<sub>AGC-GM</sub> 2
N<sub>GM</sub> 28 N<sub>NZD</sub> 6
N<sub>LDL</sub> 3 N<sub>XW</sub> 4
N<sub>DCS</sub> 5
(2)具体实施时，计算单独建模电厂机组的阻尼系数MD：NZD为-2.28，LDL为-1.08，XW为-1.40，DCS为2.45。选择阻尼系数MD从小到大排列的前三个电厂机组(NZD，LDL，XW)作为调速器参数设计的目标机组。根据云南电网的方式数据，限定各个设计参数范围为0<(KNZDP,KLDLP,KXWP)<5、0<(KNZDI,KLDLI,KXWI)<1、0<(KNZDD,KLDLD,KXWD)<3。
(3)具体实施时，根据以下方法扫频测量各个运行方式下统一频率模型的开环频率响应数据：
1)在统一频率模型中，将除调速器参数设计目标电厂机组调速器以外的系统统一为开环系统P(s)。
2)通过将s＝jω代入P(s)扫频测量得到P(s)的频率响应数据，其中包括P(s)的幅频特性数据|P(jω)|与相频特性数据∠P(jω)。以云南电网冬大运行方式为例，LDL机组调速器除外的开环系统P(s)的频率响应数据如图10所示。
(4)具体实施时，利用以下方法求解各个运行方式下考虑系统稳定约束的参数解空间S1：
以云南电网冬大运行方式为例，LDL机组PID型调速器参数为待设计参数。
1)LDL电厂PID型调速器与开环系统P(s)构成的反馈系统闭环稳定的等价条件为：回差矩阵行列式F(s)＝1+(KP·s+KI+KD·s2)/s·P(s)不存在右半平面的零点，其中，KP、KI、KD为PID调速器待设计参数；
2)以采样点Kp＝0为例，ω的奇数重根集{-0.8244,-0.3138,0,0.3138,0.8244}。
3)对于Kp＝0，求得的序列I仅有一条序列，为{1,-1,1,-1,1}。
4)对于求得的序列I＝{1,-1,1,-1,1}生成的参数解空间如下：
上式即为Kp＝0时考虑系统稳定的解空间。再求得所有KP采样点对应的解空间，即考虑系统稳定约束的参数解空间S1。求解得到的云南电网冬大运行方式下考虑系统稳定约束的参数解空间S1如图11所示；
(5)具体实施时，利用以下方法求解各个运行方式下考虑调节性能约束的参数解空间S2：
以云南电网冬大运行方式为例，LDL机组PID型调速器参数为待设计参数。
1)常见的考虑电力系统调节性能的约束条件描述如下；
||w(s)S(s)||∞<1
其中：w(s)为性能权函数，S(s)为系统的灵敏度函数，其表达式如下：
式中：s为拉普拉斯算子，M＝1.5，ωB *＝0.0001，A＝0；
2)上述约束条件等价于LDL电厂调速器与如下复系数传递函数Pc(s,θr)构成的反馈系统闭环稳定，
对θ进行采样得到样本集[θ1θ2…θm]，即对于每一个样本θr，LDL电厂调速器与Pc(s,θr)构成的反馈系统的回差矩阵行列式F(s)＝1+(KP·s+KI+KD·s2)/s·Pc(s,θr)不存在右半平面的零点。接下来以样本点θ＝0.3307为例。
3)以Kp＝0为例，ω的奇数重根集为{-0.7515-0.3118,0,0.3155,0.9055}。
4)对于Kp＝0求得的序列I仅有一条序列，为{1,-1,1,-1,1}。
5)对于求得的序列I＝{1,-1,1,-1,1}可生成额参数解空间如下：
上式即为θ＝0.3307、Kp＝0时考虑系统性能的解空间。利用上步骤将所有的θ样本点对应的解空间求交集，即得到Kp＝0时考虑系统性能的解空间。再求得所有KP采样点对应的解空间，就得到了考虑系统性能约束的解空间S2。求解得到的云南电网冬大运行方式下考虑系统性能约束的参数解空间S2如图12所示；
(6)具体实施时，通过对云南电网各个运行方式下的S1和S2求交集得到云南电网冬大方式下各个运行方式下同时考虑系统稳定与调节性能的参数解空间S如图13所示。
(7)具体实施时，通过对各个运行方式下的解空间S求交集得到云南电网多运行方式下同时考虑系统稳定与调节性能的参数解空间Sm如图14所示。
(8)具体实施时，通过以下方法基于参数解空间Sm设计LDL机组调速器参数：
1)由于上述线性约束构成的解空间Sm，选择目标函数为max KP·KI，问题转化为一个简单的非线性规划问题。
2)利用matlab自带的函数fmincon即可求解上述优化模型，得到待设计的水轮机调速器参数如表19所示。
再利用前述步骤，逐次整定得到NZD以及XW的调速器参数如表19所示。
表19
图15显示了时域仿真结果，展示了系统设计参数比原始参数增强了系统的稳定性并且提升了调速器的调节性能。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (10)

1.基于频率响应的水轮机PID型调速器参数设计方法，其特征在于，包括以下步骤：
1)建立含有系统主要发电机组在内的统一频率模型；
2)选择调速器参数设计的目标电厂机组；
3)扫频测量得到各个运行方式下统一频率模型的开环频率响应数据；
4)根据3)得到的统一频率模型的开环频率响应数据，求解各个运行方式下考虑系统稳定约束的参数解空间S1；
5)根据3)得到的统一频率模型的开环频率响应数据，求解各个运行方式下考虑调节性能约束的参数解空间S2；
6)根据4)得到参数解空间S1的和5)得到的参数解空间S2，求解各个运行方式下同时考虑系统稳定与调节性能的参数解空间S；
7)根据6)得到参数解空间S，求解多运行方式下同时考虑系统稳定与调节性能的参数解空间Sm；
8)基于参数解空间Sm设计调速器参数。
2.根据权利要求1所述的基于频率响应的水轮机PID型调速器参数设计方法，其特征在于，所述的步骤1)中的建立含有系统主要发电机组的统一频率模型步骤：
(1)统计系统所含发电机组调速器和原动机模型的类型以及它们各自占系统总容量的比例；
(2)对超低频振荡发生时各电厂机组有功功率振幅从大到小排列，选择排列在第一个或前几个的电厂机组单独建模；将各类型调速器和原动机模型按照其各自占系统总容量的比例从大到小排序，选择其中占比最大的1～4种模型搭建于统一频率模型中；
(3)对于(2)中选定的调速器和原动机模型，将其参数设置为系统中出力最大的该类型机组的参数；同时按照以下方式设置统一频率模型中各类型机组的台数：统计各类型机组在系统中的总有功出力，则统一频率模型中各类型机组的台数等于其在系统中的总有功出力除以它的单台有功出力；
(4)根据系统自动发电控制的安排设置统一频率模型中参与自动发电控制的机组台数：将统一频率模型中参与自动发电控制的机组台数设置为系统中参与自动发电控制的机组总有功出力除以单台机组的有功出力的值。
3.根据权利要求1所述的基于频率响应的水轮机PID型调速器参数设计方法，其特征在于，所述的步骤2)中的选择调速器参数设计的目标电厂机组步骤：
(1)计算所有单独建模电厂机组的阻尼系数MD；
对第i台机组，其阻尼系数MDi为：
其中，Re()为取实部运算，GTi(s)为电厂机组调速器-原动机传递函数，s为拉普拉斯算子，j为虚数单位，ωc为振荡频率，单位为rad/s；
(2)对于阻尼系数MD从小到大排列，选择排列在第一个或前几个的电厂机组为调速器参数设计目标电厂机组。
4.根据权利要求1所述的基于频率响应的水轮机PID型调速器参数设计方法，其特征在于，所述的步骤3)中的扫频测量得到各个运行方式下统一频率模型的开环频率响应数据步骤：
在统一频率模型中，将除调速器参数设计目标电厂机组调速器以外的系统统一为开环系统P(s)；
通过扫频测量得到各个运行方式下统一频率模型的开环频率响应数据，即P(s)的频率响应数据，其中，P(s)的频率响应数据包括P(s)的幅频特性数据|P(jω)|与相频特性数据∠P(jω)。
5.根据权利要求4所述的基于频率响应的水轮机PID型调速器参数设计方法，其特征在于，所述的步骤4)中的求解各个运行方式下考虑系统稳定约束的参数解空间S1步骤：
(1)目标电厂调速器与开环系统P(s)构成的反馈系统闭环稳定的等价条件为：回差矩阵行列式F(s)＝1+(KP·s+KI+KD·s2)/s·P(s)不存在右半平面的零点，其中，KP、KI、KD为PID调速器待设计参数；
(2)定义开环系统P(s)的相移数为σ＝(∠P(j∞)-∠P(-j∞))/(π/2)；
(3)定义开环系统P(s)的相对次数为
(4)将KP参数离散采样得到[KP1，KP2，...，KPn]；定义如下函数：
对于每一个KP采样点，根据上式函数曲线求解的奇数重根集为{ω1，ω2，...，ωl-1}；
(5)定义并求解如下序列：
I＝{i1，i2，...，i1-1}
其中，it∈{+1，-1}，t＝1，2，...，l-1；
如果rP-1是偶数，求解如下序列：
[2i1-2i2+2i3+…+(-1)l-12il-1]·Sx＝-σ-2
如果rP-1是奇数，求解如下序列：
[2i1-2i2+2i3+…+(-1)l-12il-1]·Sx＝-σ-2
(6)对于求得的第k条序列Ik＝{i1，k，i2，k，...，il-1，k}，求得如下参数解空间：
由上式可知每条序列Ik对应的解空间是由线性约束构成的解空间，记为再将所有序列的解空间求并集就得到考虑系统稳定的解空间其中U为求并集运算；
其中，s为拉普拉斯算子，|P(jω)|为幅频特性数据，∠P(jω)为相频特性数据，j为虚数单位，s＝jω。
6.根据权利要求4所述的基于频率响应的水轮机PID型调速器参数设计方法，其特征在于，所述的步骤5)中的求解各个运行方式下考虑调节性能约束的参数解空间S2步骤：
(1)电力系统调节性能的约束条件描述如下：
||w(s)S(s)||∞<1
其中：w(s)为性能权函数，S(s)为系统的灵敏度函数，其表达式如下：
式中：s为拉普拉斯算子，M为灵敏度函数峰值上界，ωB *为带宽频率下界，A为最大稳态跟踪误差；
该约束条件等价于目标电厂调速器与如下定义的复系数传递函数Pc(s,θr)构成的反馈系统闭环稳定：
对θ进行采样得到样本集[θ1 θ2 … θm]，即对于每一个样本θr，目标电厂调速器与Pc(s,θr)构成的反馈系统的回差矩阵行列式F(s)＝1+(KP·s+KI+KD·s2)/s·Pc(s,θr)不存在右半平面的零点；
(2)定义Pc(s,θr)的相移数为σ＝(∠Pc(j∞,θr)-∠Pc(-j∞,θr))/(π/2)；
(3)定义Pc(s,θr)的相对次数为
(4)将KP参数离散采样得到[KP1,KP2,…,KPn]；定义如下函数：
对于每一个KP采样点，根据上式函数曲线求解的奇数重根集为{ω1,ω2,...,ωl-1}；
(5)定义并求解如下序列：
I′＝{i′1,i′2,...,i′l-1}
其中，i′t∈{+1,-1},t＝1,2,...,l-1；
求解如下序列：
[2i′1-2i′2+2i′3+…+(-1)l-12i′l-1]·Sx＝-σ-2
(6)对于求得的第k条序列I′k＝{i′1,k,i′2,k,...,i′l-1,k}可求得如下参数解空间：
由上式可知每条序列Ik对应的解空间是由线性约束构成的解空间，记为再将所有序列的解空间求并集就得到θr对应的解空间其中∪为求并集运算；
(7)将根据(6)求得的所有θr对应的解空间求交集就得到考虑系统调节性能约束的解空间其中∩为求交集运算；
其中，s为拉普拉斯算子，|P(jω)|为幅频特性数据，∠P(jω)为相频特性数据，j为虚数单位，s＝jω。
7.根据权利要求1所述的基于频率响应的水轮机PID型调速器参数设计方法，其特征在于，所述的步骤6)中的求解各个运行方式下同时考虑系统稳定与调节性能的参数解空间S＝S1∩S2。
8.根据权利要求1所述的基于频率响应的水轮机PID型调速器参数设计方法，其特征在于，所述的步骤7)中的求解多运行方式下同时考虑系统稳定与调节性能的参数解空间Sm为各个运行方式下解空间S的交集。
9.根据权利要求1所述的基于频率响应的水轮机PID型调速器参数设计方法，其特征在于，所述的步骤8)中的基于参数解空间Sm设计调速器参数步骤：
根据得到的解空间Sm，设置目标函数为max KP·KI，进行求解，所得到的KP、KI、KD为最终的水轮机调速器设计参数。
10.根据权利要求9所述的基于频率响应的水轮机PID型调速器参数设计方法，其特征在于，利用matlab自带的函数fmincon求解max KP·KI。
CN201811046089.XA 2018-09-07 2018-09-07 基于频率响应的水轮机pid型调速器参数设计方法 Active CN109185036B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811046089.XA CN109185036B (zh) 2018-09-07 2018-09-07 基于频率响应的水轮机pid型调速器参数设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811046089.XA CN109185036B (zh) 2018-09-07 2018-09-07 基于频率响应的水轮机pid型调速器参数设计方法

Publications (2)

Publication Number Publication Date
CN109185036A CN109185036A (zh) 2019-01-11
CN109185036B true CN109185036B (zh) 2019-06-18

Family

ID=64915358

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811046089.XA Active CN109185036B (zh) 2018-09-07 2018-09-07 基于频率响应的水轮机pid型调速器参数设计方法

Country Status (1)

Country Link
CN (1) CN109185036B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112502883B (zh) * 2019-09-16 2023-07-14 国电南瑞南京控制系统有限公司 考虑频率稳定特性的水轮机调速器参数整定方法及装置
CN116931438B (zh) * 2023-09-12 2024-01-23 南方电网调峰调频发电有限公司检修试验分公司 一种水轮机调速器参数确定方法、装置、设备及介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107171345A (zh) * 2017-07-03 2017-09-15 云南电网有限责任公司 用于分析不确定参数对电力系统超低频振荡影响的方法
CN107476931A (zh) * 2017-08-14 2017-12-15 国网四川省电力公司电力科学研究院 一种水轮机调速器pid参数优化方法及系统
CN107658892A (zh) * 2017-11-09 2018-02-02 云南电网有限责任公司 一种多机电力系统超低频振荡稳定性评估方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107171345A (zh) * 2017-07-03 2017-09-15 云南电网有限责任公司 用于分析不确定参数对电力系统超低频振荡影响的方法
CN107476931A (zh) * 2017-08-14 2017-12-15 国网四川省电力公司电力科学研究院 一种水轮机调速器pid参数优化方法及系统
CN107658892A (zh) * 2017-11-09 2018-02-02 云南电网有限责任公司 一种多机电力系统超低频振荡稳定性评估方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于值集法对云南电网超低频振荡的稳定分析;周靖皓 等;《电网技术》;20171031;第41卷(第10期);全文
水轮机微机调速器控制策略分析与研究;陆桂明 等;《中国农村水利水电》;20101230(第7期);全文
水轮机调速器参数最佳整定—频率响应法;沈祖诒;《华东水利学院学报》;19841230(第01期);全文

Also Published As

Publication number Publication date
CN109185036A (zh) 2019-01-11

Similar Documents

Publication Publication Date Title
US12140332B2 (en) Optimized regulating and controlling method and system for integrated electricity and heat system with heat pumps
CN108964091B (zh) 一种考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法
Biegel et al. Distributed low-complexity controller for wind power plant in derated operation
Liu et al. Active power control of solar PV generation for large interconnection frequency regulation and oscillation damping
Zhang et al. Multiobjective optimization of a fractional‐order pid controller for pumped turbine governing system using an improved nsga‐Iii algorithm under multiworking conditions
Thanos et al. Dynamic data driven adaptive simulation framework for automated control in microgrids
CN107194625A (zh) 基于神经网络的风电场弃风电量评估方法
Zhang et al. Forecasting solar photovoltaic power production at the aggregated system level
CN109787251A (zh) 一种集群温控负荷聚合模型、系统参数辨识及反推控制方法
Hasan et al. Real time simulation of automatic generation control for interconnected power system
CN111509785A (zh) 一种用于电网多源最优协同控制的方法、系统及存储介质
CN107658892A (zh) 一种多机电力系统超低频振荡稳定性评估方法
CN109185036B (zh) 基于频率响应的水轮机pid型调速器参数设计方法
Neshati et al. Hardware-in-the-loop testing of wind turbine nacelles for electrical certification on a dynamometer test rig
Abdul Kadir et al. Optimal placement and sizing of photovoltaic based distributed generation considering costs of operation planning of monocrystalline and thin-film technologies
CN109345129A (zh) 一种光伏发电站健康状态评估方法和系统
Sørensen et al. Wind power variability and power system reserves in South Africa
Chen et al. Sliding-mode design with disturbance observer optimized by improved HBA algorithm for frequency regulation of multi-area power system
Karki et al. Reliability assessment of a wind power delivery system
Chen et al. Variation-cognizant probabilistic power flow analysis via multi-task learning
Zhongping et al. Integrated wind and solar power forecasting in China
Silva-Saravia et al. Effect of wind farm spatial correlation on oscillation damping in the wecc system
Pierce et al. Wind turbine control system modeling capabilities
Banta et al. Advanced control of a real smart polygeneration microgrid
Abrishambaf et al. Emulation of Energy Scheduling for an Agriculture Irrigation System Considering Field Data

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