CN108964091B - 一种考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法 - Google Patents
一种考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法 Download PDFInfo
- Publication number
- CN108964091B CN108964091B CN201811045389.6A CN201811045389A CN108964091B CN 108964091 B CN108964091 B CN 108964091B CN 201811045389 A CN201811045389 A CN 201811045389A CN 108964091 B CN108964091 B CN 108964091B
- Authority
- CN
- China
- Prior art keywords
- frequency
- power plant
- model
- unit
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000013461 design Methods 0.000 title claims abstract description 19
- 230000010355 oscillation Effects 0.000 claims abstract description 21
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 50
- 238000012546 transfer Methods 0.000 claims description 19
- 238000013016 damping Methods 0.000 claims description 13
- 238000010408 sweeping Methods 0.000 claims description 10
- 238000010248 power generation Methods 0.000 claims description 8
- 230000035945 sensitivity Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 3
- 238000003491 array Methods 0.000 claims description 2
- 230000001172 regenerating effect Effects 0.000 claims description 2
- ZZUFCTLCJUWOSV-UHFFFAOYSA-N furosemide Chemical compound C1=C(Cl)C(S(=O)(=O)N)=CC(C(O)=O)=C1NCC1=CC=CO1 ZZUFCTLCJUWOSV-UHFFFAOYSA-N 0.000 claims 3
- 238000004088 simulation Methods 0.000 abstract description 12
- 230000002401 inhibitory effect Effects 0.000 abstract description 3
- 238000012795 verification Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
Classifications
-
- 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
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明涉及一种考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法，用于同时设计多个水电厂机组的调速器参数，以在抑制系统超低频振荡现象的同时保证跟踪性能满足要求，属于电网安全技术领域。该评估方法包括：建立含有系统主要发电机组的统一频率模型；选择调速器参数设计目标机组；确定扫频频率；将统一频率模型重构为模型；寻找使得模型在扫频频率内结构奇异值上界小于1的调速器参数为所设计的参数。本发明提出的设计方法简单易实现，对不确定性的模拟精确，且只要运行方式变化带来摄动不超过预设范围就能保证所设计参数的有效性，可应用于多机电力系统考虑运行方式不确定性时水轮机调速器参数设计，使系统安全稳定运行。
Description
技术领域
本发明属于电网安全技术领域，涉及一种水轮机调速器参数设计方法，具体涉及一种考虑方式不确定性的多机电力系统水轮机调速器参数设计方法。
背景技术
交流系统超低频振荡现象在国内外实际电网运行中均有发现，具有持续时间长、振荡频率极低(小于0.1Hz)的特点，威胁系统安全稳定运行。自从云南电网与南方电网主网异步联网后，云南电网也出现了超低频振荡现象。已有研究对超低频振荡问题的机理做了较为深入的剖析，发现水电机组在超低频段内无法为系统提供足够的阻尼转矩，超低频振荡与系统中的水电机组一次调频关系密切。抑制超低频振荡最简单且有效的手段是切除水电机组调速器，然而切除调速器会降低机组跟踪负荷波动的性能，只能作为应急手段。
调速器参数设计需综合考虑系统稳定性与跟踪性能，且应具有足够的鲁棒性应对运行方式变化导致的系统摄动。结构奇异值理论能够兼顾鲁棒稳定与鲁棒性能，是调速器控制参数设计的有力工具，且从未应用于解决超低频振荡问题。
发明内容
本发明的目的是为了解决现有技术的不足，提供一种抑制超低频振荡问题的水轮机调速器参数设计方法，该方法兼顾了系统的稳定性要求与跟踪性能要求，且只要系统实际运行中由运行方式变化带来的摄动不超过预定值就能保证所设计的调速器参数有效性。
为实现上述目的，本发明采用的技术方案如下：
一种考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法，包括以下步骤：
1)建立含有全系统主要发电机组在内的统一频率模型；
2)选择调速器参数设计的目标电厂机组；
3)确定扫频频率；
4)将统一频率模型重构为M-Δ模型；
5)基于M-Δ模型扫频设计调速器参数。
进一步，优选的是，所述的步骤1)中的建立含有全系统主要发电机组在内的统一频率模型步骤：
(1)统计全系统所含发电机组调速器和原动机模型的类型以及它们各自占全系统总容量的比例；
(2)对超低频振荡发生时各电厂机组有功功率振幅从大到小排列，选择其中振幅最大的一个或几个电厂机组单独建模；对除单独建模的电厂机组外的剩余机组，将各类型调速器和原动机模型按照其各自占剩余机组有功出力总容量的比例从大到小排序，选择其中占比最大的1～4种模型搭建于统一频率模型中；
(3)对于(2)中选定的调速器和原动机模型，将其参数设置为全系统中出力最大的该类型机组的参数；同时按照以下方式设置统一频率模型中各类型机组的台数：统计各类型机组在全系统中的总有功出力，则统一频率模型中各类型机组的台数等于其在全系统中的总有功出力除以它的单台有功出力；
(4)根据全系统自动发电控制的安排设置统一频率模型中参与自动发电控制的机组台数：将统一频率模型中参与自动发电控制的机组台数设置为全系统中参与自动发电控制的机组总有功出力除以单台机组的有功出力的值。
进一步，优选的是，所述的步骤2)中的选择调速器参数设计的目标电厂机组步骤：
(1)计算所有单独建模电厂机组的阻尼系数MD；
对第i台机组，其阻尼系数MDi为：其中，Re()为取实部运算，GTi(s)为电厂机组调速器-原动机传递函数，s为拉普拉斯算子，j为虚数单位，ωc为振荡频率，单位为rad/s；
(2)对于阻尼系数MD从小到大排列，选择排列在第一个或前几个的电厂机组为调速器参数设计目标电厂机组。
进一步，优选的是，所述的步骤3)中的确定扫频频率中，扫频频率[ω1,ω2,...,ωN]为一组数组，且需包含全系统发生超低频振荡时所有可能的振荡频率区间；其中ω1,ω2,...,ωN为扫频频率点，单位为rad/s。
进一步，优选的是，所述的步骤4)中的将统一频率模型重构为M-Δ模型步骤：假设编号为1～n的电厂机组为目标电厂机组；
(1)在统一频率模型中，将除调速器参数设计的目标电厂机组外的系统统一为外部系统P(s)：
其中：s为拉普拉斯算子，MS为外部系统等效惯性常数，D为负荷阻尼系数，GTi(s)为第i台电厂机组调速器-原动机传递函数；
(2)统一频率模型中由运行方式不确定性带来的摄动如下：Δ1,…,Δn为电厂机组数目N1,…,Nn的摄动，ΔP为外部系统的摄动；
(3)在目标电厂机组总功率偏差∑ΔPm-∑ΔPe与目标电厂机组总负荷波动∑ΔPe之间插入用于表示性能的虚拟摄动Δf，且Δf的无穷范数|Δf|∞<1；其中，∑ΔPm为目标电厂机组总有功功率波动；
(4)设置摄动Δ1,…,Δn,ΔP和Δf的权函数分别为W1,…,Wn,WP(s)和WS(s)；其中：W1,…,Wn和WP(s)应能够包络N1,…,Nn和P(s)所有可能偏离标称情况的最大范围；WS(s)为性能权函数，表达式为：
式中：s为拉普拉斯算子，M为灵敏度函数峰值上界，ωB *为带宽频率下界，A为最大稳态跟踪误差；
(5)通过推导[u1,...,un,uP,uf]T＝M·[v1,...,vn,vP,vf]T的关系得到M矩阵，从而将统一频率模型重构为具有对角结构型摄动Δ＝diag{Δ1,...,Δn,ΔP,Δf}的M-Δ模型：
式中
其中：u1,…,un,uP,uf为Δ1,…,Δn,ΔP和Δf的输入信号，v1,…,vn,vP,vf为Δ1,…,Δn,ΔP和Δf的输出信号，GTi为目标电厂机组i的调速器-原动机传递函数，Ni为目标电厂机组i数目的标称值，GTj为目标电厂机组j的调速器-原动机传递函数，Nj为目标电厂机组j数目的标称值，P外部系统传递函数，Wi,WP,Wf分别为Δi,ΔP,Δf的权函数。
进一步，优选的是，所述的步骤5)基于M-Δ模型扫频设计调速器参数步骤：
(1)对外部系统Pk(s)扫频，设置每个扫频点处标称系统P(jωm)＝avg(Pk(jωm))，权函数WP(jωm)＝max(|Pk(jωm)-P(jωm)|)，其中：Pk(jωm)为方式k的外部系统在扫频频率点ωm处的频率特性，j为虚数单位，ωm为扫频频率点，k为运行方式编号；
(2)设置目标电厂机组数目的标称值和权函数；对目标电厂机组i，电厂机组数目的标称值Ni＝(max(Nik)+min(Nik))/2，权函数Wi＝(max(Nik)-min(Nik))/2，其中：Nik为方式k中目标电厂机组i数目，k为运行方式编号；
(3)设置性能权函数WS(s)参数，带入并扫频得到每个扫频点处的WS(jωm)，其中：j为虚数单位，ωm为扫频频率点；
s为拉普拉斯算子，M为灵敏度函数峰值上界，ωB *为带宽频率下界，A为最大稳态跟踪误差；
(4)在各目标电厂机组调速器参数KP、KI、KD容许变化范围内，随机生成一组调速器参数KP、KI、KD，其中：KP为比例参数，KI为积分参数，KD为微分参数；
(5)将上述调速器参数带入各电厂机组调速器-原动机传递函数GTi(s)，扫频得到每个扫频点处目标电厂机组调速器-原动机传递函数频率特性GTi(jωm)，m＝1,2,…,N，其中：j为虚数单位，ωm为扫频频率点，m为扫频频率点序号；
(6)在每个扫频点处将(1)～(5)中得到结果带入M-Δ模型，计算系统结构奇异值上界其中，j为虚数单位，ωi为扫频频率点，m为扫频频率序号；
(7)若以下条件满足
其中，j为虚数单位，ωm为扫频频率点，m为扫频频率点序号；
则(4)中生成的KP、KI、KD满足要求，为最终的设计参数；否则返回(4)，重新生成一组调速器参数，并重复(5)～(6)，直至满足(7)所示的条件。
本发明中，M为灵敏度函数峰值上界，是根据稳定要求设置；ωB *为带宽频率下界，是根据跟踪性能要求设置；A为最大稳态跟踪误差，设置为很小值，如0.0001。
本发明与现有技术相比，其有益效果为：
本发明基于结构奇异值理论，将运行方式不确定性处理为统一频率模型相应位置的摄动，且通过在功率偏差与负荷波动间插入虚拟摄动的方法将调速器设计需要考虑的鲁棒跟踪性能问题转化为鲁棒稳定问题。将摄动与系统分离，原系统被重构为包含对角型摄动的M-Δ模型，从而设计过程为寻找使M-Δ模型结构奇异值上界小于1的水轮机调速器参数。
现有技术大部分只能针对单一工况设计单个调速器参数，无法应对多工况、多调速器参数整定问题。本发明可以同时设计多个调速器参数，且所设计的参数能够保证在所考虑的工况范围内系统性能均满足要求，解决了现有技术的不足。
附图说明
图1为本发明方法的流程示意图。
图2为本发明方法中统一频率模型示意图。
图3为本发明实施例仿真验证中云南电网统一频率模型示意图。
图4为本发明实施例仿真验证中GS型调速器模型。
图5为本发明实施例仿真验证中TB型汽轮机模型。
图6为本发明实施例仿真验证中GH型机组模型。
图7为本发明实施例仿真验证中GM型调速器模型。
图8为本发明实施例仿真验证中GA型伺服器模型。
图9为本发明实施例仿真验证中水轮机详细模型。
图10为本发明实施例仿真验证中重构后的统一频率模型示意图。
图11为本发明实施例仿真验证中权函数WP(s)的幅频特性。
图12为本发明实施例仿真验证中系统结构奇异值上界。
具体实施方式
下面结合实施例对本发明作进一步的详细描述。
本领域技术人员将会理解，下列实施例仅用于说明本发明，而不应视为限定本发明的范围。实施例中未注明具体技术或条件者，按照本领域内的文献所描述的技术或条件或者按照产品说明书进行。所用材料或设备未注明生产厂商者，均为可以通过购买获得的常规产品。
如图1所示，采用本发明方法进行处理时，通过建立含有全系统主要发电机组调速器和原动机模型的统一频率模型，可以实现对原全系统超低频振荡的模拟，避免了使用复杂的详细模型，减小了模型复杂度；通过确定包括振荡频率在内的扫频频率[ω1,ω2,...,ωN]，可以确定需要计算结构奇异值上界的频率点；基于建立的统一频率模型，可以将全系统重构为参数设计使用的M-Δ模型；通过将统一频率模型重构为M-Δ模型，并且随机生成一组目标机组调速器参数KP、KI、KD，扫频计算各扫频频率[ω1,ω2,...,ωN]处的结构奇异值上界若在各扫频频率[ω1,ω2,...,ωN]处结构奇异值上界则随机生成的目标机组调速器参数KP、KI、KD为最终的设计参数。
本发明的具体实施例如下：
(1)统计南方电网2017冬大、冬小、夏大、夏大极限、夏小、夏小外送、汛期方式中云南电网的机组数据和潮流计算结果的基础上，建立统一频率模型。按照之前所述步骤，在Matlab/Simulink软件中建立统一频率模型如图3所示。
以冬大方式为例，上述的统一频率模型建立遵循如下步骤：
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型调速器模型如图7所示，模式选择为开度模式、硬反馈输入信号为YPID，参数如表1。
表1
其所含GA型伺服模型如图8所示，参数如表2。
表2
| Tc | 11.4 | KD | 0 |
| To | 18.5 | KI | 0 |
| T2 | 0.02 | Td | 0 |
| KP | 10 |
其所含水轮机详细模型如图9所示，图中水系统传递函数F(s)为
其中，
其中，Tep为压力水管弹性时间常数，Ts为调压室时间常数，TWc为引水洞水击时间常数，TWp为压力水管水击时间常数，Zp为压力水管水阻抗，φc为引水洞摩擦系数，φp为压力水管摩擦系数，s为拉普拉斯算子。参数如表3。
表3
LDL电厂机组功率基值为400MW，在其自身功率基值下机组惯性为9.61。
其所含GM型调速器模型如图7所示，模式选择为开度模式、硬反馈输入信号为YPID，参数如表4。
表4
| TR1 | 0.02 | TD | 0.198 |
| KW | 1.14 | T1 | 0 |
| KP | 2.5 | T2 | 0 |
| KD | 1 | TR2 | 0.02 |
| KI | 0.5 | eP | 0.04 |
其所含GA型伺服模型如图8所示，参数如表5。
表5
| Tc | 9.35 | KD | 0 |
| To | 23 | KI | 0 |
| T2 | 0.2 | Td | 0 |
| KP | 15 |
其所含水轮机详细模型如图9所示，图中水系统传递函数F(s)为
其中，
其中，Tep为压力水管弹性时间常数，Ts为调压室时间常数，TWc为引水洞水击时间常数，TWp为压力水管水击时间常数，Zp为压力水管水阻抗，φc为引水洞摩擦系数，φp为压力水管摩擦系数，s为拉普拉斯算子。参数如如表6。
表6
XW电厂机组功率基值为778MW，在其自身功率基值下机组惯性为8.72。其所含GM型调速器模型如图7所示，模式选择为开度模式、硬反馈输入信号为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型伺服模型如图8所示，参数如表8。
表8
| Tc | 10 | KD | 0 |
| To | 13.1 | KI | 0 |
| T2 | 0.02 | Td | 0 |
| KP | 6 |
其所含水轮机详细模型如图9所示，图中水系统传递函数F(s)为
其中，
其中，Tep为压力水管弹性时间常数，Ts为调压室时间常数，TWc为引水洞水击时间常数，TWp为压力水管水击时间常数，Zp为压力水管水阻抗，φc为引水洞摩擦系数，φp为压力水管摩擦系数，s为拉普拉斯算子。参数如表9。
表9
DCS电厂机组功率基值为250MW，在其自身功率基值下机组惯性为9.63。其所含GM型调速器模型如图7所示，模式选择为开度模式、硬反馈输入信号为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型伺服模型如图8所示，参数如表11。
表11
| Tc | 21.5 | KD | 0 |
| To | 22.5 | KI | 0 |
| T2 | 0.2 | Td | 0 |
| KP | 40 |
其所含水轮机详细模型如图9所示，图中水系统传递函数F(s)为
其中，
其中，Tep为压力水管弹性时间常数，Ts为调压室时间常数，TWc为引水洞水击时间常数，TWp为压力水管水击时间常数，Zp为压力水管水阻抗，φc为引水洞摩擦系数，φp为压力水管摩擦系数，s为拉普拉斯算子。参数如表12。
表12
对于典型机组GS，功率基值为353MW，在其自身功率基值下机组惯性为8.5。其所含GS型调速器模型如图4所示，TB型汽轮机模型如图5所示，参数如表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型机组模型如图6所示，参数如表14。
表14
对于典型机组GM，功率基值为667MW，在其自身功率基值下机组惯性为10.38。其所含GM型调速器模型如图7所示，模式选择为开度模式、硬反馈输入信号为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型伺服模型如图8所示，参数如表16。
表16
| Tc | 13 | KD | 0 |
| To | 13 | KI | 0 |
| T2 | 0.2 | Td | 0 |
| KP | 20 |
其所含水轮机详细模型如图9所示，图中水系统传递函数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)作为调速器参数设计的目标机组。
(3)由于对所有方式，全系统超低频振荡频率的范围为0.25～0.38rad/s，因此设定扫频频率范围为0.19～0.63rad/s。
(4)将统一频率模型重构为如图10所示的M-Δ模型。
上述的重构遵循如下步骤，
1)将除NZD、LDL、XW电厂机组外的系统统一为外部系统P(s)：
其中：s为拉普拉斯算子，MS为外部系统等效惯性常数，D为负荷阻尼系数。
2)统一频率模型中由运行方式不确定性带来的摄动如下：ΔNZD、ΔLDL、ΔXW为电厂机组数目NNZD、NLDL、NXW的摄动，ΔP为外部系统的摄动；
3)在目标电厂机组总功率偏差∑ΔPm-∑ΔPe与目标电厂机组总负荷波动∑ΔPe之间插入用于表示性能的虚拟摄动Δf，且Δf的无穷范数|Δf|∞<1，以把鲁棒性能问题转化为鲁棒稳定问题；其中，∑ΔPm为目标电厂机组总有功功率波动；
4)设置摄动ΔNZD、ΔLDL、ΔXW、ΔP和Δf的权函数WNZD、WLDL、WXW、WP和WS。其中：WNZD、WLDL、WXW、WP应能够包络NNZD、NLDL、NXW和P可能偏离标称情况的最大范围；WS为性能权函数，表达式为
式中：s为拉普拉斯算子，M为灵敏度函数峰值上界(根据稳定要求设置)，ωB *为带宽频率下界(根据跟踪性能要求设置)，A为最大稳态跟踪误差(设置为很小值，如0.0001)；
5)通过推导[uNZD,uLDL,uXW,uP,uf]T＝M·[vNZD,vLDL,vXW,vP,vf]T的关系得到M矩阵，从而将全系统重构为具有对角结构型摄动Δ＝diag{ΔNZD,ΔLDL,ΔXW,ΔP,Δf}的M-Δ模型，其中：uNZD,uLDL,uXW,uP,uf为ΔNZD、ΔLDL、ΔXW、ΔP和Δf的输入信号，vNZD,vLDL,vXW,vP,vf为ΔNZD、ΔLDL、ΔXW、ΔP和Δf的输出信号。由图10易知以下关系成立
推导得到如下M-Δ模型
其中
(5)随机生成NZD、LDL、XW电厂机组水轮机调速器参数KNZDP、KNZDI、KNZDD，KLDLP、KLDLI、KLDLD，KXWP、KXWI、KXWD。其中各参数随机生成的范围为0<(KNZDP,KLDLP,KXWP)<10、0<(KNZDI,KLDLI,KXWI)<0.5、0<(KNZDD,KLDLD,KXWD)<3。
(6)利用Matlab的robust control toolbox中的mussv()函数，扫频计算各扫频点处结构奇异值上界检验(5)中随机生成的参数是否满足要求。
详细步骤如下：
1)对所有七个方式外部系统Pk(s)扫频，设置每个扫频点处标称系统P(jωm)＝avg(Pk(jωm))，权函数WP(jωm)＝max(|Pk(jωm)-P(jωm)|)，其中：j为虚数单位，ωm为扫频频率点，k为运行方式(冬大、冬小、夏大、夏大极限、夏小、夏小外送、汛期)。WP幅频特性如图11所示；
2)设置电厂机组的标称数目和权函数。由于在所有七个运行方式中，NZD电厂机组数目变化范围为2～5，LDL电厂机组数目变化范围为1～3，XW电厂机组数目变化范围为2～4，因此设置标称数目NNZD＝3.5，NLDL＝2，NXW＝3，WNZD＝1.5，WLDL＝1，WXW＝1；
3)设置性能权函数WS(s)参数，带入并扫频得到每个扫频点处的WS(jωm)，其中：j为虚数单位，ωm为扫频频率点，
s为拉普拉斯算子，M＝10，ωB *＝0.13，A＝0.0001；
4)将随机生成的KNZDP、KNZDI、KNZDD，KLDLP、KLDLI、KLDLD，KXWP、KXWI、KXWD带入GTNZD,GTLDL,GTXW,
s为拉普拉斯算子。扫频得到每个频率点处的GTNZD(jωm),GTLDL(jωm),GTXW(jωm)；
5)在每个扫频点处，带入前述步骤计算得到的结果，得到每个扫频点处的M矩阵为一个数值矩阵，利用matlab robust control toolbox中的mussv()函数计算每个扫频点处的结构奇异值上界若满足如下条件
则随机生成的KNZDP、KNZDI、KNZDD，KLDLP、KLDLI、KLDLD，KXWP、KXWI、KXWD为最终设计得到的调速器参数，否则回到步骤(5)重新随机生成一组调速器参数后重复上述步骤。最终设计得到的调速器参数如表19。
表19
| K<sub>NZDP</sub> | 0.87 | K<sub>LDLP</sub> | 7.51 | K<sub>XWP</sub> | 0.95 |
| K<sub>NZDI</sub> | 0.06 | K<sub>LDLI</sub> | 0.15 | K<sub>XWI</sub> | 0.14 |
| K<sub>NZDD</sub> | 1 | K<sub>LDLD</sub> | 1 | K<sub>XWD</sub> | 1 |
其结构奇异值上界如图12所示，可见原始参数不满足鲁棒性能，设计参数满足鲁棒性能。
为了验证分析结果，计算了各方式全系统在原始参数和设计参数下的主导特征值与闭环带宽频率如表20。
表20
可见对于所有运行方式，设计参数在保证全系统闭环带宽(跟踪性能)的基础上大幅提高了主导模式的阻尼(稳定性)，从而说明了本发明方法的有效性。
因此，本发明方法能够用于设计满足稳定性与跟踪性能要求的水轮机调速器参数，且所设计的参数有足够的鲁棒性应对方式不确定性，本发明方法具有突出显著的技术效果。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (5)
1.一种考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法，其特征在于，包括以下步骤：
1)建立含有全系统主要发电机组在内的统一频率模型；
2)选择调速器参数设计的目标电厂机组；
3)确定扫频频率；
4)将统一频率模型重构为M-Δ模型；
5)基于M-Δ模型扫频设计调速器参数；所述的步骤4)中的将统一频率模型重构为M-Δ模型步骤：假设编号为1～n的电厂机组为目标电厂机组；
(1)在统一频率模型中，将除调速器参数设计的目标电厂机组外的系统统一为外部系统P(s)：
其中：s为拉普拉斯算子，MS为外部系统等效惯性常数，D为负荷阻尼系数，GTi(s)为第i台电厂机组调速器-原动机传递函数；
(2)统一频率模型中由运行方式不确定性带来的摄动如下：Δ1,…,Δn为电厂机组数目N1,…,Nn的摄动，ΔP为外部系统的摄动；
(3)在目标电厂机组总功率偏差∑ΔPm-∑ΔPe与目标电厂机组总负荷波动∑ΔPe之间插入用于表示性能的虚拟摄动Δf，且Δf的无穷范数|Δf|∞<1；其中，∑ΔPm为目标电厂机组总有功功率波动；
(4)设置摄动Δ1,…,Δn,ΔP和Δf的权函数分别为W1,…,Wn,WP(s)和WS(s)；其中：W1,…,Wn和WP(s)应能够包络N1,…,Nn和P(s)所有可能偏离标称情况的最大范围；WS(s)为性能权函数，表达式为：
式中：s为拉普拉斯算子，M为灵敏度函数峰值上界，ωB *为带宽频率下界，A为最大稳态跟踪误差；
(5)通过推导[u1,...,un,uP,uf]T＝M·[v1,...,vn,vP,vf]T的关系得到M矩阵，从而将统一频率模型重构为具有对角结构型摄动Δ＝diag{Δ1,...,Δn,ΔP,Δf}的M-Δ模型：
式中
其中：u1,…,un,uP,uf为Δ1,…,Δn,ΔP和Δf的输入信号，v1,…,vn,vP,vf为Δ1,…,Δn,ΔP和Δf的输出信号，GTi为目标电厂机组i的调速器-原动机传递函数，Ni为目标电厂机组i数目的标称值，GTj为目标电厂机组j的调速器-原动机传递函数，Nj为目标电厂机组j数目的标称值，P为外部系统传递函数，Wi,WP,Wf分别为Δi,ΔP,Δf的权函数。
2.根据权利要求1所述的考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法，其特征在于，所述的步骤1)中的建立含有全系统主要发电机组在内的统一频率模型步骤：
(1)统计全系统所含发电机组调速器和原动机模型的类型以及它们各自占全系统总容量的比例；
(2)对超低频振荡发生时各电厂机组有功功率振幅从大到小排列，选择排列在第一个或前几个的电厂机组单独建模；对除单独建模的电厂机组外的剩余机组，将各类型调速器和原动机模型按照其各自占剩余机组有功出力总容量的比例从大到小排序，选择其中占比最大的1～4种模型搭建于统一频率模型中；
(3)对于(2)中选择的调速器和原动机模型，将其参数设置为全系统中出力最大的该类型机组的参数；同时按照以下方式设置统一频率模型中各类型机组的台数：统计各类型机组在全系统中的总有功出力，则统一频率模型中各类型机组的台数等于其在全系统中的总有功出力除以它的单台有功出力；
(4)根据全系统自动发电控制的安排设置统一频率模型中参与自动发电控制的机组台数：将统一频率模型中参与自动发电控制的机组台数设置为全系统中参与自动发电控制的机组总有功出力除以单台机组的有功出力的值。
3.根据权利要求1所述的考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法，其特征在于，所述的步骤2)中的选择调速器参数设计的目标电厂机组步骤：
(1)计算所有单独建模电厂机组的阻尼系数MD；
对第i台机组，其阻尼系数MDi为：
其中，Re()为取实部运算，GTi(s)为电厂机组调速器-原动机传递函数，s为拉普拉斯算子，j为虚数单位，ωc为振荡频率，单位为rad/s；
(2)对于阻尼系数MD从小到大排列，选择排列在第一个或前几个的电厂机组为调速器参数设计目标电厂机组。
4.根据权利要求1所述的考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法，其特征在于，所述的步骤3)中的确定扫频频率中，扫频频率[ω1,ω2,...,ωN]为一组数组，且需包含全系统发生超低频振荡时所有可能的振荡频率区间；其中ω1,ω2,...,ωN为扫频频率点，单位为rad/s。
5.根据权利要求1所述的考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法，其特征在于，所述的步骤5)基于M-Δ模型扫频设计调速器参数步骤：
(1)对外部系统Pk(s)扫频，设置每个扫频点处标称系统P(jωm)＝avg(Pk(jωm))，权函数WP(jωm)＝max(|Pk(jωm)-P(jωm)|)，其中：Pk(jωm)为方式k的外部系统在扫频频率点ωm处的频率特性，j为虚数单位，ωm为扫频频率点，k为运行方式编号；
(2)设置目标电厂机组数目的标称值和权函数；对目标电厂i，电厂机组数目的标称值Ni＝(max(Nik)+min(Nik))/2，权函数Wi＝(max(Nik)-min(Nik))/2，其中：Nik为方式k中目标电厂机组i数目，k为运行方式编号；
(3)设置性能权函数WS(s)参数，带入并扫频得到每个扫频点处的WS(jωm)，
s为拉普拉斯算子，M为灵敏度函数峰值上界，ωB *为带宽频率下界，A为最大稳态跟踪误差；
(4)在各目标电厂机组调速器参数KP、KI、KD容许变化范围内，随机生成一组调速器参数KP、KI、KD；
5)将上述调速器参数带入各电厂机组调速器-原动机传递函数GTi(s)，扫频得到每个扫频点处目标电厂机组调速器-原动机传递函数频率特性GTi(jωm)，m＝1,2,…,N；N为扫频频率点的个数；
(6)在每个扫频点处将(1)～(5)中得到结果带入M-Δ模型，计算系统结构奇异值上界其中，j为虚数单位，m为扫频频率序号；
(7)若以下条件满足N为扫频频率点的个数；
则(4)中生成的KP、KI、KD满足要求，为最终的设计参数；否则返回(4)，重新生成一组调速器参数，并重复(5)～(6)，直至满足(7)所示的条件。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811045389.6A CN108964091B (zh) | 2018-09-07 | 2018-09-07 | 一种考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811045389.6A CN108964091B (zh) | 2018-09-07 | 2018-09-07 | 一种考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108964091A CN108964091A (zh) | 2018-12-07 |
| CN108964091B true CN108964091B (zh) | 2019-09-20 |
Family
ID=64476110
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811045389.6A Active CN108964091B (zh) | 2018-09-07 | 2018-09-07 | 一种考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108964091B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109499696B (zh) * | 2018-12-17 | 2019-11-08 | 东北大学 | 一种多机驱动高频振动磨机的参数确定方法 |
| CN110739728B (zh) * | 2019-11-30 | 2024-10-01 | 中国南方电网有限责任公司 | 一种基于机组惯性功率比的电网机组运行方式的确定方法 |
| CN112003317B (zh) * | 2020-06-05 | 2022-01-25 | 国电南瑞科技股份有限公司 | 一种水电机组功率调节输入输出优化方法及装置 |
| CN113255485B (zh) * | 2021-05-13 | 2022-05-17 | 云南电网有限责任公司 | 一种水力发电机组并网模式的识别方法及装置 |
| CN117134345B (zh) * | 2023-09-05 | 2024-06-18 | 北京交通大学 | 基于结构奇异值法的抑制超低频振荡的控制方法及系统 |
| CN116931438B (zh) * | 2023-09-12 | 2024-01-23 | 南方电网调峰调频发电有限公司检修试验分公司 | 一种水轮机调速器参数确定方法、装置、设备及介质 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107171345B (zh) * | 2017-07-03 | 2019-11-05 | 云南电网有限责任公司 | 用于分析不确定参数对电力系统超低频振荡影响的方法 |
| CN107658892B (zh) * | 2017-11-09 | 2018-08-28 | 云南电网有限责任公司 | 一种多机电力系统超低频振荡稳定性评估方法 |
-
2018
- 2018-09-07 CN CN201811045389.6A patent/CN108964091B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108964091A (zh) | 2018-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108964091B (zh) | 一种考虑运行方式不确定性的多机电力系统水轮机调速器参数设计方法 | |
| Xu et al. | Modeling, tuning, and validating system dynamics in synthetic electric grids | |
| Yao et al. | An optimal power-flow approach to improve power system voltage stability using demand response | |
| CN103138256A (zh) | 一种新能源电力消纳全景分析系统及方法 | |
| CN108695857B (zh) | 风电场自动电压控制方法、装置及系统 | |
| CN106130004B (zh) | 一种综合考虑稳定特性的并网点新能源接纳能力的评估方法 | |
| CN110429648B (zh) | 考虑风速随机波动的小干扰稳定裕度概率评估方法 | |
| CN107658892A (zh) | 一种多机电力系统超低频振荡稳定性评估方法 | |
| CN109193750A (zh) | 基于pidd2控制器的风电集群参与agc方法 | |
| Neshati et al. | Hardware-in-the-loop testing of wind turbine nacelles for electrical certification on a dynamometer test rig | |
| Wang et al. | Cooperative output regulation of large-scale wind turbines for power reserve control | |
| CN109185036B (zh) | 基于频率响应的水轮机pid型调速器参数设计方法 | |
| Kharrazi | Tuning of power system stabilizer in Single Machine Infinite Bus (SMIB) using genetic algorithm and Power Factory Modal Analysis | |
| Singh et al. | Optimal tuning of temporary droop governor of hydro power plant using genetic algorithm | |
| CN116339134B (zh) | 一种大扰动火电机组调频优化控制系统 | |
| CN104167735B (zh) | 一种风电场非机理等效建模方法与装置 | |
| Zhang et al. | Model-Free Fast Frequency Support of Wind Farms for Tracking Optimal Frequency Trajectory | |
| CN118094965B (zh) | 一种虚拟电厂频率响应聚合等值方法、装置及存储介质 | |
| Rosen et al. | Optimal tuning of governors on synchronous generators in a multi-area power system with a large photovoltaic plant | |
| Xu et al. | Frequency response model of offshore wind turbines participating in frequency modulation | |
| CN114513005B (zh) | 考虑电流限制和出力水平的分布式光伏机电暂态建模方法 | |
| CN114529208B (zh) | 考虑chp机组快速爬坡能力约束的电热耦合系统动态优化调度方法 | |
| Curtis Saunders et al. | The effect of timescales on wind farm power variability with nonlinear model predictive control | |
| Mugo et al. | Optimization of Turbine Blade Pitch Angle of a Home-built Wind Turbine for Maximum Power Output | |
| CN105226704A (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 |