CN113595152A - 基于火电机组调节速率包络线的电网agc指令优化分配方法及系统 - Google Patents

基于火电机组调节速率包络线的电网agc指令优化分配方法及系统 Download PDF

Info

Publication number
CN113595152A
CN113595152A CN202110980828.8A CN202110980828A CN113595152A CN 113595152 A CN113595152 A CN 113595152A CN 202110980828 A CN202110980828 A CN 202110980828A CN 113595152 A CN113595152 A CN 113595152A
Authority
CN
China
Prior art keywords
generating unit
thermal power
power generating
instruction
power grid
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
Application number
CN202110980828.8A
Other languages
English (en)
Other versions
CN113595152B (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.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Shandong Electric Power 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 State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202110980828.8A priority Critical patent/CN113595152B/zh
Publication of CN113595152A publication Critical patent/CN113595152A/zh
Application granted granted Critical
Publication of CN113595152B publication Critical patent/CN113595152B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本申请提供了一种基于火电机组调节速率包络线的电网AGC指令优化分配方法,属于电网AGC指令优化调度领域;具体方案为:根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;根据所述样本集合获得火电机组调节速率包络线;以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据火电机组调节速率包络线对电网AGC指令进行分配。本申请保证了在最短调节时间内响应完成电网AGC指令,在实际电网AGC指令优化调度方面具有良好的应用价值。本申请还提供了一种基于火电机组调节速率包络线的电网AGC指令优化分配系统。

Description

基于火电机组调节速率包络线的电网AGC指令优化分配方法 及系统
技术领域
本申请涉及电网AGC指令调度优化技术领域,特别涉及一种基于火电机组调节速率包络线的电网AGC指令优化分配方法及系统。
背景技术
本部分的陈述仅仅是提供了与本申请相关的背景技术,并不必然构成现有技术。
随着我国经济的快速发展,社会对电能的需求不断地增长,电网容量不断扩大,各大电网的峰谷差日趋增大。自动发电控制(Automatic Generation Control,AGC)能在大负荷变化或其它干扰下维持系统频率稳定性,是提高电力系统安全性的辅助服务之一。当电能的产生和使用不同步,即存在一个偏差时,控制器根据该偏差产生一个总AGC指令,电网调度中心依照某种分配方法将电网AGC指令分配给各个机组。各个机组的分散控制系统自动按分配的负荷情况,对锅炉和汽机运行工况进行调整,增大或减小出力以达到负荷指令的要求。经研究发现,机组在不同出力时具有不同的调节速率,即机组调节速率是动态变化的,所以,根据此规律充分发挥大容量火电机组的调峰作用,建立科学合理的负荷指令优化分配方法,对于电网平稳运行具有重要意义。
现有的负荷指令优化分配方法存在局限性。首先,大量现有方法没有考虑火电机组实际调节速率;其次,在考虑火电机组调节速率的分配方法中,忽略了火电机组不同出力下火电机组调节速率的差异性,以及该差异性对电网调节时间的重要影响。
发明内容
为了解决现有技术的不足,本申请提供了一种基于火电机组调节速率包络线的电网AGC指令优化分配方法及系统,基于火电机组输出功率与调节速率的样本集合,通过多项式拟合获得火电机组调节速率包络线,用于描述火电机组输出功率与调节速率的关系,根据火电机组调节速率包络线对电网AGC指令进行优化分配,使电网总功率最短时间内跟踪AGC指令变化。
为了实现上述目的,本申请采用如下技术方案:
本申请第一方面提供了一种基于火电机组调节速率包络线的电网AGC指令优化分配方法,包括以下步骤:
根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;
根据上述样本集合获得火电机组调节速率包络线;
以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据火电机组调节速率包络线对电网AGC指令进行分配。
可选地,上述根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合的步骤,包括:
将火电机组实发功率历史数据按时间序列划分成多个数据段,通过每个数据段的实发功率幅值变化量和持续时间的比值计算火电机组调节速率,即
Figure BDA0003229049420000021
其中,δ(m)为第m段数据的实发功率幅值变化量,τ(m)为第m段数据的持续时间;
火电机组调节速率r(m)对应的火电机组输出功率为y(nm),记P(m)=y(nm),根据分段结果得到火电机组输出功率与调节速率的样本集合
Figure BDA0003229049420000031
可选地,上述根据样本集合获得火电机组调节速率包络线的步骤,包括:
基于局部密度剔除上述样本集合中的边缘样本点,对筛选后的样本集合采用多项式拟合获得火电机组调节速率包络线。
可选地,上述以最短电网AGC指令响应时间为优化目标,建立电网AGC指令优化分配模型的步骤,包括:建立电网AGC指令分配模型如下:
Figure BDA0003229049420000032
以各台火电机组分配的AGC指令调整量之和等于电网AGC指令调整量为约束条件:
Figure BDA0003229049420000033
Figure BDA0003229049420000034
其中,T为电网AGC指令响应时间,ΔPi为分配给第i台火电机组的AGC指令调整量,ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
可选地,上述电网AGC指令响应时间T在数值上等于电网中调节时间最长的火电机组的负荷响应时间,T=max{t1,t2,…,tI},ti=ΔPi/ri
可选地,上述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,还包括:通过拉格朗日乘子法求解电网AGC指令分配模型,获得电网AGC指令优化分配模型,得到第i台火电机组分配的AGC指令调整量为:
Figure BDA0003229049420000035
ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
可选地,上述根据火电机组调节速率包络线对电网AGC指令进行分配的步骤,获得第i台火电机组的AGC指令调整量为:
Figure BDA0003229049420000041
其中,ri(P)为第i台火电机组输出功率为P时的调节速率,ΔP为电网AGC指令调整量。
本申请第二方面提供了一种基于火电机组调节速率包络线的电网AGC指令优化分配系统。
基于火电机组调节速率包络线的电网AGC指令优化分配系统,包括:
数据处理模块,被配置为:根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;
动态关系拟合模块,被配置为:根据上述样本集合获得火电机组调节速率包络线;
指令分配模块,被配置为:以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据火电机组调节速率包络线对电网AGC指令进行分配。
可选地,上述根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合的步骤,包括:
将火电机组实发功率历史数据按时间序列划分成多个数据段,通过每个数据段的实发功率幅值变化量和持续时间的比值计算火电机组调节速率,即
Figure BDA0003229049420000042
其中,δ(m)为第m段数据的实发功率幅值变化量,τ(m)为第m段数据的持续时间;
火电机组调节速率r(m)对应的火电机组输出功率为y(nm),记P(m)=y(nm),根据分段结果得到火电机组输出功率与调节速率的样本集合
Figure BDA0003229049420000051
可选地,上述根据样本集合获得火电机组调节速率包络线的步骤,包括:
基于局部密度剔除样本集合中的边缘样本点,对筛选后的样本集合采用多项式拟合获得火电机组调节速率包络线。
可选地,上述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,包括:
建立电网AGC指令分配模型,电网AGC指令分配模型如下:
Figure BDA0003229049420000052
以各台火电机组分配的AGC指令调整量之和等于电网AGC指令调整量为约束条件:
Figure BDA0003229049420000053
Figure BDA0003229049420000054
其中,T为电网AGC指令响应时间,ΔPi为分配给第i台火电机组的AGC指令调整量,ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
可选地,上述电网AGC指令响应时间T在数值上等于电网中调节时间最长的火电机组的负荷响应时间,T=max{t1,t2,…,tI},ti=ΔPi/ri
可选地,上述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,还包括:
通过拉格朗日乘子法求解电网AGC指令分配模型,获得电网AGC指令优化分配模型,得到第i台火电机组分配的AGC指令调整量为:
Figure BDA0003229049420000061
ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
可选地,上述根据火电机组调节速率包络线对电网AGC指令进行分配的步骤,获得第i台火电机组的AGC指令调整量为:
Figure BDA0003229049420000062
其中,ri(P)为第i台火电机组输出功率为P时的调节速率,ΔP为电网AGC指令调整量。
本申请第三方面提供了一种介质,其上存储有程序,该程序被处理器执行时实现如本申请第一方面所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法中的步骤。
本申请第四方面提供了一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的程序,所述处理器执行所述程序时实现如本申请第一方面所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法中的步骤。
与现有技术相比,本申请的有益效果是:
1、本申请实施例的方法,根据火电机组最大调节速率进行优化调度。由于不同火电机组的调节速率不同,在进行优化分配时,让调节速率强的火电机组承担更多的负荷响应任务,保证在最短时间内完成对电网总负荷指令的响应,在实际电网负荷优化调度方面具有良好的应用价值。
2、本申请实施例的方法,充分考虑了不同输出功率条件下火电机组调节速率的差异性。火电机组的调节速率不是固定不变的,即不同出力下火电机组的调节速率不同,基于火电机组历史数据信息估计火电机组出力与火电机组调节速率之间的关系。
附图说明
图1为本公开实施例1提供的基于火电机组调节速率包络线的电网AGC指令优化分配方法的流程图;
图2为本公开实施例2提供的火电机组实发功率的PLR结果示意图;
图3为本公开实施例2提供的火电机组调节速率与输出功率的集合样本分布图;
图4为本公开实施例2提供的火电机组调节速率的包络线图;
图5为本公开实施例2提供的火电机组负荷指令响应图。
具体实施方式
应该指出,以下详细说明都是示例性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
实施例1:
如图1所示,本公开实施例提供了基于火电机组调节速率包络线的电网AGC指令优化分配方法,包括以下步骤:根据火电机组实发功率历史数据,获得火电机组输出功率与该输出功率对应调节速率的样本集合;根据上述样本集合获得火电机组调节速率包络线;以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据火电机组调节速率包络线对电网AGC指令进行分配。
在一个实施例中,上述根据火电机组实发功率历史数据,获得火电机组输出功率与该输出功率对应调节速率的样本集合的步骤,包括:
将火电机组实发功率历史数据按时间序列划分成多个数据段,通过每个数据段的实发功率幅值变化量和持续时间的比值计算火电机组调节速率,即
Figure BDA0003229049420000081
其中,δ(m)为第m段数据的实发功率幅值变化量,τ(m)为第m段数据的持续时间;
火电机组调节速率r(m)对应的火电机组输出功率为y(nm),记P(m)=y(nm),根据分段结果得到火电机组输出功率与调节速率的样本集合
Figure BDA0003229049420000082
在一个实施例中,上述根据样本集合获得火电机组调节速率包络线的步骤,包括:基于局部密度剔除所述样本集合中的边缘样本点,对筛选后的样本集合采用多项式拟合获得火电机组调节速率包络线,即火电机组输出功率与调节速率的关系式。
在一个实施例中,上述以最短电网AGC指令响应时间为优化目标,建立电网AGC指令优化分配模型的步骤,包括:建立电网AGC指令分配模型,如下:
Figure BDA0003229049420000083
以各台火电机组分配的AGC指令调整量之和等于电网AGC指令调整量为约束条件:
Figure BDA0003229049420000091
Figure BDA0003229049420000092
其中,T为电网AGC指令响应时间,ΔPi为分配给第i台火电机组的AGC指令调整量,ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
可选地,电网AGC指令响应时间T在数值上等于电网中调节时间最长的火电机组的负荷响应时间,T=max{t1,t2,…,tI},ti=ΔPi/ri
在一个实施例中,上述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,还包括:通过拉格朗日乘子法求解电网AGC指令分配模型,获得电网AGC指令优化分配模型,得到第i台火电机组分配的AGC指令调整量为:
Figure BDA0003229049420000093
ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
在一个实施例中,根据火电机组调节速率包络线对上述实施例的电网AGC指令进行分配,获得第i台火电机组的AGC指令调整量为:
Figure BDA0003229049420000094
其中,ri(P)为第i台机组输出功率为P时的调节速率,ΔP为电网AGC指令调整量。
实施例2
本公开实施例提供了基于火电机组调节速率包络线的电网AGC指令优化分配方法,包括以下步骤:
步骤1:对火电机组实发功率的历史数据进行分段线性表示,将时间序列分成多个短直线段,计算每个数据段的调节速率,得到火电机组输出功率与对应调节速率的样本集合。
首先,采用自底向上分段线性表示方法将一段较长的实发功率历史数据时间序列
Figure BDA0003229049420000101
分成M个短数据段:
Figure BDA0003229049420000102
则第m段可描述为
Figure BDA0003229049420000103
其中nm为第m段的第一个点,nm+1为第m段的最后一个点,m∈[1,M]。那么,第m段数据
Figure BDA0003229049420000104
可以用线性回归模型表示:y(n)=am+bmn+v(n),其中am和bm分别是第m段数据拟合直线的截距和斜率,v(n)是均值为零方差为
Figure BDA0003229049420000105
的高斯白噪声。
使函数
Figure BDA0003229049420000106
最小,以产生参数am和bm的估计值
Figure BDA0003229049420000107
Figure BDA0003229049420000108
Figure BDA0003229049420000109
根据最小二乘估计可以得到截距的估计值
Figure BDA00032290494200001010
和斜率的估计值
Figure BDA00032290494200001011
的表达式:
Figure BDA00032290494200001012
Figure BDA00032290494200001013
其中
Figure BDA00032290494200001014
Figure BDA00032290494200001015
分别是第m段数据中y和n的平均值。
Figure BDA00032290494200001016
Figure BDA0003229049420000111
因此,第m段数据段
Figure BDA0003229049420000112
能被一条直线表示:
Figure BDA0003229049420000113
而分段数M可由“L”法确定,整个序列
Figure BDA0003229049420000114
的损失函数是y(n)和
Figure BDA0003229049420000115
之间的拟合误差总和。
Figure BDA0003229049420000116
该方法将上式中由m和L(M)的“L”型线近似为两条直线,即
Figure BDA0003229049420000117
Figure BDA0003229049420000118
其中α1、β1、α2、β2为两条直线的未知参数,M0为“L”型线拐点处的分段数。
曲线的拐点作为“M”的估计值,记为
Figure BDA0003229049420000119
Figure BDA00032290494200001110
其中上界
Figure BDA00032290494200001111
N是时间序列
Figure BDA00032290494200001112
的长度,Nmin是数据段中最小的数据样本,
Figure BDA00032290494200001113
表示向下取整函数,即取不大于自身的最大整数。Nmin的引入基于一个常识,即一个数据段应持续一段时间,例如Nminh=300s时,才能分辨出该数据段的任何特征,h为实发功率时间序列的采样周期。
火电机组实发功率的PLR结果如图2所示,离散点为机组实发功率数据,连续的线段为PLR结果。那么,第m段数据的幅值变化量和持续时间可以表示为:
Figure BDA00032290494200001114
Figure BDA0003229049420000121
其中,h为实发功率时间序列的采样周期,例如1s,δ为幅值变化量,单位为MW,τ的为持续时间,单位为min。
则第m段数据估计的火电机组调节速率为:
Figure BDA0003229049420000122
火电机组调节速率r(m)对应的机组出力为y(nm),记P(m)=y(nm),根据分段结果可以得到火电机组对应当前出力与调节速率的集合
Figure BDA0003229049420000123
步骤2:基于局部密度剔除样本集合
Figure BDA0003229049420000124
中的处于边缘的样本点,对筛选出的样本采用多项式拟合得到火电机组调节速率包络线,确定火电机组输出功率与调节速率的关系式。
从概率统计的角度来说,代表火电机组调节速率的样本点一定是多次出现的,因此通过局部密度剔除边缘样本点,从集合
Figure BDA0003229049420000125
中选择具有足够局部密集的点,形成集合Sd
样本点(r(m),P(m))的局部密度ρ(m)定义为
Figure BDA0003229049420000126
其中dc是截止距离,d(m,q)是点(r(m),P(m))和点(r(q),P(q))间的欧式距离。
Figure BDA0003229049420000127
χ(·)是一个截止函数,
Figure BDA0003229049420000128
选择具有不小于ρ0的局部密度的点组成集合Sd,其中ρ0为样本点的局部密度阈值,即
Figure BDA0003229049420000131
集合Sd中的样本点(r(m),P(m))在r-P平面上展示为具有一定密度的散点图,对于P轴,给定一个窗口宽度ε,确定在窗口宽度内火电机组调节速率r的最大值,通过滑动窗口ε得到K个r-P平面上轮廓散点,散点对应的样本值记为集合
Figure BDA0003229049420000132
对集合
Figure BDA0003229049420000133
中的样本点(r(k),P(k))进行多项式拟合,Φ为所有次数不超过F(F≤K)的多项式构成的函数类,现求
Figure BDA0003229049420000134
使得下式的W最小,其中af为次数为f的实发功率的多项式的系数,Pf为次数为f的实发功率。
Figure BDA0003229049420000135
而满足式(14)的
Figure BDA0003229049420000136
称为最小二乘拟合多项式,其中,Pk为第k个样本点所对应的实发功率,
Figure BDA0003229049420000137
为第k个样本点所对应的实发功率,次数为f,rk为第k个样本点所对应的速率。显然
Figure BDA0003229049420000138
为a0,…,af,...,aF的多元函数,因此上述问题转化为求W=W(a0,…,af,...,aF)的极值问题。由多元函数求极值的必要条件,得
Figure BDA0003229049420000139
Figure BDA00032290494200001310
式(16)是关于a0,…,af,...,aF的线性方程组,其系数矩阵是一个对称正定矩阵,故有唯一解,可以解出af(f=0,...,F),进而得到火电机组调节速率包络线:
Figure BDA0003229049420000141
记r(P)为火电机组输出功率与变化速率关系式,并定义ri(P)为第i台机组输出功率与变化速率关系式。
步骤3:为获得最短电网AGC指令响应时间时电网AGC指令的分配结果,建立如下模型:
Figure BDA0003229049420000142
约束条件:
Figure BDA0003229049420000143
Figure BDA0003229049420000144
其中,T为电网AGC指令响应时间,数值上等于电网中调节时间最长的机组的负荷响应时间,即T=max{t1,t2,…,tI},ti=ΔPi/ri,其中ΔPi为分配给第i台火电机组的AGC指令调整量,ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
为处理上述约束,采用拉格朗日乘子法引入拉格朗日算子计算得到上述模型的最优解。记拉格朗日函数
Figure BDA0003229049420000145
对拉格朗日函数求偏导可得
Figure BDA0003229049420000146
Figure BDA0003229049420000151
Figure BDA0003229049420000152
由式(22)得
Figure BDA0003229049420000153
由式(23)和式(24)得
Figure BDA0003229049420000154
分配给第i台火电机组的AGC指令调整量为
Figure BDA0003229049420000155
根据式(17)所示的火电机组输出功率和调节速率的关系式,得到基于火电机组调节速率包络线的电网AGC指令分配结果。此时,第i台火电机组的AGC指令调整量为:
Figure BDA0003229049420000156
以下是本实施例的方法在具体示例中的应用。
以某区域电网为例,该区域含有4台火电机组,火电机组参数如表1所示。现电网需要增加100MW功率输出,即ΔP=100WM,计算分配到每台火电机组的调节负荷ΔPi,使电网能在最短时间内完成总负荷调整。
表1:火电机组参数信息表
Figure BDA0003229049420000157
Figure BDA0003229049420000161
第一步,对火电机组实发功率的历史数据进行分段线性表示,将时间序列分成多个短直线段,计算每个数据段的调节速率,得到火电机组输出功率与该输出功率对应的调节速率的样本集合。以第一台火电机组为例,图3所示为火电机组输出功率与调节速率样本集合的散点图。
第二步,基于局部密度剔除样本集合中的处于边缘的样本点,如图3所示的五角星形数据点,筛选出具有一定局部密度的样本点,图4为具有一定局部密度的样本散点图。对筛选出的样本采用多项式拟合得到火电机组调节速率包络线,如图4所示,进而得到火电机组输出功率与调节速率的关系式。同理得到4台火电机组输出功率与调节速率的关系式,如表2所示。
表2:机组出力与最大调节速率关系式
机组编号 机组输出功率与调节速率的关系式
火电机组1 r<sub>1</sub>=-9.21×10<sup>-4</sup>P<sub>1</sub><sup>2</sup>+0.44P<sub>1</sub>-44.72
火电机组2 r<sub>2</sub>=-2.94×10<sup>-4</sup>P<sub>2</sub><sup>2</sup>+0.16P<sub>2</sub>-16.01
火电机组3 r<sub>3</sub>=-0.0028P<sub>3</sub><sup>2</sup>+1.31P<sub>3</sub>-141.64
火电机组4 r<sub>4</sub>=-0.0028P<sub>4</sub><sup>2</sup>+1.36P<sub>4</sub>-152.39
第三步:以最短电网AGC指令响应时间为优化目标,根据火电机组调节速率包络线对电网AGC指令进行分配,得到如图5所示的各火电机组负荷指令响应。
实施例3:
本公开实施例提供了基于火电机组调节速率包络线的电网AGC指令优化分配系统,包括:
数据处理模块,被配置为:根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;
动态关系拟合模块,被配置为:根据所述样本集合获得火电机组调节速率包络线;
指令分配模块,被配置为:以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据机组调节速率包络线对电网AGC指令进行分配。
上述系统的工作方法与上述各实施例提供的基于火电机组调节速率包络线的电网AGC指令优化分配方法相同,这里不再赘述。
实施例4:
本公开实施例提供了一种存储介质,其上存储有程序,该程序被处理器执行时实现如上述各实施例提供的基于火电机组调节速率包络线的电网AGC指令优化分配方法中的步骤,包括:
根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;
根据所述样本集合获得火电机组调节速率包络线;
以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据火电机组调节速率包络线对电网AGC指令进行分配。
上述程序实现的方法的详细步骤与上述各实施例提供的基于火电机组调节速率包络线的电网AGC指令优化分配方法相同,这里不再赘述。
实施例5:
本公开实施例提供了一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的程序,所述处理器执行所述程序时实现如上述各实施例提供的基于火电机组调节速率包络线的电网AGC指令优化分配方法中的步骤,包括:
根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;
根据所述样本集合获得火电机组调节速率包络线;
以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据火电机组调节速率包络线对电网AGC指令进行分配。
上述程序实现的方法的详细步骤与上述各实施例提供的基于火电机组调节速率包络线的电网AGC指令优化分配方法相同,这里不再赘述。
本领域内的技术人员应明白,本申请公开的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用硬件实施例、软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本公开实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(RandomAccessMemory,RAM)等。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。

Claims (16)

1.一种基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,包括以下步骤:
根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;
根据所述样本集合获得火电机组调节速率包络线;
以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据火电机组调节速率包络线对电网AGC指令进行分配。
2.如权利要求1所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,
所述根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合的步骤,包括:
将火电机组实发功率历史数据按时间序列划分成多个数据段,通过每个数据段的实发功率幅值变化量和持续时间的比值计算火电机组调节速率,即
Figure FDA0003229049410000011
其中,δ(m)为第m段数据的实发功率幅值变化量,τ(m)为第m段数据的持续时间;
火电机组调节速率r(m)对应的火电机组输出功率为y(nm),记P(m)=y(nm),根据分段结果得到火电机组输出功率与调节速率的样本集合
Figure FDA0003229049410000012
3.如权利要求2所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,
所述根据所述样本集合获得火电机组调节速率包络线的步骤,包括:
基于局部密度剔除所述样本集合中的边缘样本点,对筛选后的样本集合采用多项式拟合获得火电机组调节速率包络线。
4.如权利要求1所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,
所述以最短电网AGC指令响应时间为优化目标,建立电网AGC指令优化分配模型的步骤,包括:建立电网AGC指令分配模型,如下:
Figure FDA0003229049410000021
以各台火电机组分配的AGC指令调整量之和等于电网AGC指令调整量为约束条件:
Figure FDA0003229049410000022
Figure FDA0003229049410000023
其中,T为电网AGC指令响应时间,ΔPi为分配给第i台火电机组的AGC指令调整量,ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
5.如权利要求4所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,
所述电网AGC指令响应时间T在数值上等于电网中调节时间最长的火电机组的负荷响应时间,T=max{t1,t2,…,tI},ti=ΔPi/ri
6.如权利要求4所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,
所述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,还包括:
通过拉格朗日乘子法求解电网AGC指令分配模型,获得电网AGC指令优化分配模型,得到第i台火电机组分配的AGC指令调整量为:
Figure FDA0003229049410000031
ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
7.如权利要求6所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,
所述根据火电机组调节速率包络线对电网AGC指令进行分配的步骤,获得第i台机组的AGC指令调整量为:
Figure FDA0003229049410000032
其中,ri(P)为第i台火电机组输出功率为P时的调节速率,ΔP为电网AGC指令调整量。
8.一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,包括:
数据处理模块,被配置为:根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;
动态关系拟合模块,被配置为:根据所述样本集合获得火电机组调节速率包络线;
指令分配模块,被配置为:以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据机组调节速率包络线对电网AGC指令进行分配。
9.如权利要求8所述的一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,
所述根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合的步骤,包括:
将火电机组实发功率历史数据按时间序列划分成多个数据段,通过每个数据段的实发功率幅值变化量和持续时间的比值计算火电机组调节速率,即
Figure FDA0003229049410000041
其中,δ(m)为第m段数据的实发功率幅值变化量,τ(m)为第m段数据的持续时间;
火电机组调节速率r(m)对应的火电机组输出功率为y(nm),记P(m)=y(nm),根据分段结果得到火电机组输出功率与调节速率的样本集合
Figure FDA0003229049410000042
10.如权利要求8所述的一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,
所述根据所述样本集合获得火电机组调节速率包络线的步骤,包括:
基于局部密度剔除所述样本集合中的边缘样本点,对筛选后的样本集合采用多项式拟合获得火电机组调节速率包络线。
11.如权利要求8所述的一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,
所述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,包括:
建立电网AGC指令分配模型,电网AGC指令分配模型如下:
Figure FDA0003229049410000043
以各台火电机组分配的AGC指令调整量之和等于电网AGC指令调整量为约束条件:
Figure FDA0003229049410000051
Figure FDA0003229049410000052
其中,T为电网AGC指令响应时间,ΔPi为分配给第i台火电机组的AGC指令调整量,ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
12.如权利要求11所述的一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,
所述电网AGC指令响应时间T在数值上等于电网中调节时间最长的火电机组的负荷响应时间,T=max{t1,t2,…,tI},ti=ΔPi/ri
13.如权利要求11所述的一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,
所述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,还包括:
通过拉格朗日乘子法求解电网AGC指令分配模型,获得电网AGC指令优化分配模型,得到第i台火电机组分配的AGC指令调整量为:
Figure FDA0003229049410000053
ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
14.如权利要求13所述的一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,
所述根据火电机组调节速率包络线对电网AGC指令进行分配的步骤,获得第i台机组的AGC指令调整量为:
Figure FDA0003229049410000061
其中,ri(P)为第i台火电机组输出功率为P时的调节速率,ΔP为电网AGC指令调整量。
15.一种介质,其上存储有程序,其特征在于,该程序被处理器执行时实现如权利要求1-7任一项所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法中的步骤。
16.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的程序,其特征在于,所述处理器执行所述程序时实现如权利要求1-7任一项所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法中的步骤。
CN202110980828.8A 2021-08-25 2021-08-25 基于火电机组调节速率包络线的电网agc指令优化分配方法及系统 Active CN113595152B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110980828.8A CN113595152B (zh) 2021-08-25 2021-08-25 基于火电机组调节速率包络线的电网agc指令优化分配方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110980828.8A CN113595152B (zh) 2021-08-25 2021-08-25 基于火电机组调节速率包络线的电网agc指令优化分配方法及系统

Publications (2)

Publication Number Publication Date
CN113595152A true CN113595152A (zh) 2021-11-02
CN113595152B CN113595152B (zh) 2024-05-17

Family

ID=78239662

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110980828.8A Active CN113595152B (zh) 2021-08-25 2021-08-25 基于火电机组调节速率包络线的电网agc指令优化分配方法及系统

Country Status (1)

Country Link
CN (1) CN113595152B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116454890A (zh) * 2023-04-20 2023-07-18 中国南方电网有限责任公司 基于scuc模型的机组组合控制方法、装置和设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108985490A (zh) * 2018-06-14 2018-12-11 内蒙古大学 基于滑动窗口特征数据拟合的风力发电机功率曲线修正方法
WO2020143104A1 (zh) * 2019-01-08 2020-07-16 南京工程学院 一种考虑阻塞和储能分时电价的电网混合滚动调度方法
CN111798060A (zh) * 2020-07-08 2020-10-20 山东科技大学 一种基于机组爬坡率估计的功率指令优化分配方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108985490A (zh) * 2018-06-14 2018-12-11 内蒙古大学 基于滑动窗口特征数据拟合的风力发电机功率曲线修正方法
WO2020143104A1 (zh) * 2019-01-08 2020-07-16 南京工程学院 一种考虑阻塞和储能分时电价的电网混合滚动调度方法
CN111798060A (zh) * 2020-07-08 2020-10-20 山东科技大学 一种基于机组爬坡率估计的功率指令优化分配方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116454890A (zh) * 2023-04-20 2023-07-18 中国南方电网有限责任公司 基于scuc模型的机组组合控制方法、装置和设备
CN116454890B (zh) * 2023-04-20 2024-02-06 中国南方电网有限责任公司 基于scuc模型的机组组合控制方法、装置和设备

Also Published As

Publication number Publication date
CN113595152B (zh) 2024-05-17

Similar Documents

Publication Publication Date Title
CN105207253B (zh) 考虑风电及频率不确定性的agc随机动态优化调度方法
CN111092429B (zh) 一种柔性互联配电网的优化调度方法、存储介质及处理器
CN109508895A (zh) 烟草制丝设备的控制性能评估装置、方法以及存储介质
CN108171372B (zh) 一种存在准备时间的多品种生产分批调度方法
CN113595152A (zh) 基于火电机组调节速率包络线的电网agc指令优化分配方法及系统
CN108512258A (zh) 一种基于改进多智能体一致性算法的风电场有功调度方法
CN109507961A (zh) 一种半导体生产线动态负荷均衡投料控制方法
Malmir et al. Controlling megawatt class WECS by ANFIS network trained with modified genetic algorithm
CN114844118A (zh) 一种适用于微电网的多类型设备功率协调控制方法及系统
CN111798060A (zh) 一种基于机组爬坡率估计的功率指令优化分配方法
CN113241780B (zh) 一种基于秃鹰搜索算法的电网二次调频控制方法与系统
CN111985771B (zh) 基于电网频率超欠调分析的电网频率调控方法及系统
CN107578191A (zh) 自动进行pmc配料计算的方法
CN115173489B (zh) 一种基于二分法的火电集群调度方法及系统
CN109190212B (zh) 复杂地形陆上风电场智能化自动布机方法及其实现系统
CN111193255A (zh) 一种考虑风电不确定性的电力系统时变母线负荷模型
CN113536685B (zh) 一种基于神经网络外推的风速概率模型的建模方法
CN114138025A (zh) 控制脱硫浆液阀门开度的方法、装置及相关设备
CN112131703A (zh) 一种基于matlab功能模块的风机运行特性分析系统
CN106655189A (zh) 基于虚拟发电机的最优潮流协调优化方法
CN111900740A (zh) 一种基于需求响应设备的电力系统调频方法和系统
CN107341590A (zh) 一种基于电价时段的节能空间评估计算方法
CN117674142B (zh) 一种基于分时电量核算的电力调度方法及调度装置
JP3707103B2 (ja) 供給割当量決定装置
CN110380410B (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