CN113595152B - 基于火电机组调节速率包络线的电网agc指令优化分配方法及系统 - Google Patents
基于火电机组调节速率包络线的电网agc指令优化分配方法及系统 Download PDFInfo
- Publication number
- CN113595152B CN113595152B CN202110980828.8A CN202110980828A CN113595152B CN 113595152 B CN113595152 B CN 113595152B CN 202110980828 A CN202110980828 A CN 202110980828A CN 113595152 B CN113595152 B CN 113595152B
- Authority
- CN
- China
- Prior art keywords
- thermal power
- power
- unit
- power grid
- agc
- 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 57
- 238000009826 distribution Methods 0.000 claims abstract description 80
- 238000005457 optimization Methods 0.000 claims abstract description 45
- 230000004044 response Effects 0.000 claims abstract description 42
- 230000001105 regulatory effect Effects 0.000 claims description 32
- 230000006870 function Effects 0.000 claims description 21
- 238000010586 diagram Methods 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 5
- 238000004220 aggregation Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000004590 computer program Methods 0.000 description 8
- 238000003860 storage Methods 0.000 description 6
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 241000287196 Asthenes Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001360 synchronised effect Effects 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/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
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- 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)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本申请提供了一种基于火电机组调节速率包络线的电网AGC指令优化分配方法,属于电网AGC指令优化调度领域;具体方案为:根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;根据所述样本集合获得火电机组调节速率包络线;以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据火电机组调节速率包络线对电网AGC指令进行分配。本申请保证了在最短调节时间内响应完成电网AGC指令,在实际电网AGC指令优化调度方面具有良好的应用价值。本申请还提供了一种基于火电机组调节速率包络线的电网AGC指令优化分配系统。
Description
技术领域
本申请涉及电网AGC指令调度优化技术领域,特别涉及一种基于火电机组调节速率包络线的电网AGC指令优化分配方法及系统。
背景技术
本部分的陈述仅仅是提供了与本申请相关的背景技术,并不必然构成现有技术。
随着我国经济的快速发展,社会对电能的需求不断地增长,电网容量不断扩大,各大电网的峰谷差日趋增大。自动发电控制(Automatic Generation Control,AGC)能在大负荷变化或其它干扰下维持系统频率稳定性,是提高电力系统安全性的辅助服务之一。当电能的产生和使用不同步,即存在一个偏差时,控制器根据该偏差产生一个总AGC指令,电网调度中心依照某种分配方法将电网AGC指令分配给各个机组。各个机组的分散控制系统自动按分配的负荷情况,对锅炉和汽机运行工况进行调整,增大或减小出力以达到负荷指令的要求。经研究发现,机组在不同出力时具有不同的调节速率,即机组调节速率是动态变化的,所以,根据此规律充分发挥大容量火电机组的调峰作用,建立科学合理的负荷指令优化分配方法,对于电网平稳运行具有重要意义。
现有的负荷指令优化分配方法存在局限性。首先,大量现有方法没有考虑火电机组实际调节速率;其次,在考虑火电机组调节速率的分配方法中,忽略了火电机组不同出力下火电机组调节速率的差异性,以及该差异性对电网调节时间的重要影响。
发明内容
为了解决现有技术的不足,本申请提供了一种基于火电机组调节速率包络线的电网AGC指令优化分配方法及系统,基于火电机组输出功率与调节速率的样本集合,通过多项式拟合获得火电机组调节速率包络线,用于描述火电机组输出功率与调节速率的关系,根据火电机组调节速率包络线对电网AGC指令进行优化分配,使电网总功率最短时间内跟踪AGC指令变化。
为了实现上述目的,本申请采用如下技术方案:
本申请第一方面提供了一种基于火电机组调节速率包络线的电网AGC指令优化分配方法,包括以下步骤:
根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;
根据上述样本集合获得火电机组调节速率包络线;
以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据火电机组调节速率包络线对电网AGC指令进行分配。
可选地,上述根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合的步骤,包括:
将火电机组实发功率历史数据按时间序列划分成多个数据段,通过每个数据段的实发功率幅值变化量和持续时间的比值计算火电机组调节速率,即
其中,δ(m)为第m段数据的实发功率幅值变化量,τ(m)为第m段数据的持续时间;
火电机组调节速率r(m)对应的火电机组输出功率为y(nm),记P(m)=y(nm),根据分段结果得到火电机组输出功率与调节速率的样本集合
可选地,上述根据样本集合获得火电机组调节速率包络线的步骤,包括:
基于局部密度剔除上述样本集合中的边缘样本点,对筛选后的样本集合采用多项式拟合获得火电机组调节速率包络线。
可选地,上述以最短电网AGC指令响应时间为优化目标,建立电网AGC指令优化分配模型的步骤,包括:建立电网AGC指令分配模型如下:
以各台火电机组分配的AGC指令调整量之和等于电网AGC指令调整量为约束条件:
其中,T为电网AGC指令响应时间,ΔPi为分配给第i台火电机组的AGC指令调整量,ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
可选地,上述电网AGC指令响应时间T在数值上等于电网中调节时间最长的火电机组的负荷响应时间,T=max{t1,t2,…,tI},ti=ΔPi/ri。
可选地,上述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,还包括:通过拉格朗日乘子法求解电网AGC指令分配模型,获得电网AGC指令优化分配模型,得到第i台火电机组分配的AGC指令调整量为:
ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
可选地,上述根据火电机组调节速率包络线对电网AGC指令进行分配的步骤,获得第i台火电机组的AGC指令调整量为:
其中,ri(P)为第i台火电机组输出功率为P时的调节速率,ΔP为电网AGC指令调整量。
本申请第二方面提供了一种基于火电机组调节速率包络线的电网AGC指令优化分配系统。
基于火电机组调节速率包络线的电网AGC指令优化分配系统,包括:
数据处理模块,被配置为:根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;
动态关系拟合模块,被配置为:根据上述样本集合获得火电机组调节速率包络线;
指令分配模块,被配置为:以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据火电机组调节速率包络线对电网AGC指令进行分配。
可选地,上述根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合的步骤,包括:
将火电机组实发功率历史数据按时间序列划分成多个数据段,通过每个数据段的实发功率幅值变化量和持续时间的比值计算火电机组调节速率,即
其中,δ(m)为第m段数据的实发功率幅值变化量,τ(m)为第m段数据的持续时间;
火电机组调节速率r(m)对应的火电机组输出功率为y(nm),记P(m)=y(nm),根据分段结果得到火电机组输出功率与调节速率的样本集合
可选地,上述根据样本集合获得火电机组调节速率包络线的步骤,包括:
基于局部密度剔除样本集合中的边缘样本点,对筛选后的样本集合采用多项式拟合获得火电机组调节速率包络线。
可选地,上述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,包括:
建立电网AGC指令分配模型,电网AGC指令分配模型如下:
以各台火电机组分配的AGC指令调整量之和等于电网AGC指令调整量为约束条件:
其中,T为电网AGC指令响应时间,ΔPi为分配给第i台火电机组的AGC指令调整量,ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
可选地,上述电网AGC指令响应时间T在数值上等于电网中调节时间最长的火电机组的负荷响应时间,T=max{t1,t2,…,tI},ti=ΔPi/ri。
可选地,上述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,还包括:
通过拉格朗日乘子法求解电网AGC指令分配模型,获得电网AGC指令优化分配模型,得到第i台火电机组分配的AGC指令调整量为:
ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
可选地,上述根据火电机组调节速率包络线对电网AGC指令进行分配的步骤,获得第i台火电机组的AGC指令调整量为:
其中,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指令进行分配。
在一个实施例中,上述根据火电机组实发功率历史数据,获得火电机组输出功率与该输出功率对应调节速率的样本集合的步骤,包括:
将火电机组实发功率历史数据按时间序列划分成多个数据段,通过每个数据段的实发功率幅值变化量和持续时间的比值计算火电机组调节速率,即
其中,δ(m)为第m段数据的实发功率幅值变化量,τ(m)为第m段数据的持续时间;
火电机组调节速率r(m)对应的火电机组输出功率为y(nm),记P(m)=y(nm),根据分段结果得到火电机组输出功率与调节速率的样本集合
在一个实施例中,上述根据样本集合获得火电机组调节速率包络线的步骤,包括:基于局部密度剔除所述样本集合中的边缘样本点,对筛选后的样本集合采用多项式拟合获得火电机组调节速率包络线,即火电机组输出功率与调节速率的关系式。
在一个实施例中,上述以最短电网AGC指令响应时间为优化目标,建立电网AGC指令优化分配模型的步骤,包括:建立电网AGC指令分配模型,如下:
以各台火电机组分配的AGC指令调整量之和等于电网AGC指令调整量为约束条件:
其中,T为电网AGC指令响应时间,ΔPi为分配给第i台火电机组的AGC指令调整量,ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
可选地,电网AGC指令响应时间T在数值上等于电网中调节时间最长的火电机组的负荷响应时间,T=max{t1,t2,…,tI},ti=ΔPi/ri。
在一个实施例中,上述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,还包括:通过拉格朗日乘子法求解电网AGC指令分配模型,获得电网AGC指令优化分配模型,得到第i台火电机组分配的AGC指令调整量为:
ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
在一个实施例中,根据火电机组调节速率包络线对上述实施例的电网AGC指令进行分配,获得第i台火电机组的AGC指令调整量为:
其中,ri(P)为第i台机组输出功率为P时的调节速率,ΔP为电网AGC指令调整量。
实施例2
本公开实施例提供了基于火电机组调节速率包络线的电网AGC指令优化分配方法,包括以下步骤:
步骤1:对火电机组实发功率的历史数据进行分段线性表示,将时间序列分成多个短直线段,计算每个数据段的调节速率,得到火电机组输出功率与对应调节速率的样本集合。
首先,采用自底向上分段线性表示方法将一段较长的实发功率历史数据时间序列分成M个短数据段:/>则第m段可描述为/>其中nm为第m段的第一个点,nm+1为第m段的最后一个点,m∈[1,M]。那么,第m段数据/>可以用线性回归模型表示:y(n)=am+bmn+v(n),其中am和bm分别是第m段数据拟合直线的截距和斜率,v(n)是均值为零方差为/>的高斯白噪声。
使函数最小,以产生参数am和bm的估计值/>和/>即
根据最小二乘估计可以得到截距的估计值和斜率的估计值/>的表达式:
其中和/>分别是第m段数据中y和n的平均值。
因此,第m段数据段能被一条直线表示:/>而分段数M可由“L”法确定,整个序列/>的损失函数是y(n)和/>之间的拟合误差总和。
该方法将上式中由m和L(M)的“L”型线近似为两条直线,即
其中α1、β1、α2、β2为两条直线的未知参数,M0为“L”型线拐点处的分段数。
曲线的拐点作为“M”的估计值,记为则
其中上界N是时间序列/>的长度,Nmin是数据段中最小的数据样本,/>表示向下取整函数,即取不大于自身的最大整数。Nmin的引入基于一个常识,即一个数据段应持续一段时间,例如Nminh=300s时,才能分辨出该数据段的任何特征,h为实发功率时间序列的采样周期。
火电机组实发功率的PLR结果如图2所示,离散点为机组实发功率数据,连续的线段为PLR结果。那么,第m段数据的幅值变化量和持续时间可以表示为:
其中,h为实发功率时间序列的采样周期,例如1s,δ为幅值变化量,单位为MW,τ的为持续时间,单位为min。
则第m段数据估计的火电机组调节速率为:
火电机组调节速率r(m)对应的机组出力为y(nm),记P(m)=y(nm),根据分段结果可以得到火电机组对应当前出力与调节速率的集合
步骤2:基于局部密度剔除样本集合中的处于边缘的样本点,对筛选出的样本采用多项式拟合得到火电机组调节速率包络线,确定火电机组输出功率与调节速率的关系式。
从概率统计的角度来说,代表火电机组调节速率的样本点一定是多次出现的,因此通过局部密度剔除边缘样本点,从集合中选择具有足够局部密集的点,形成集合Sd。
样本点(r(m),P(m))的局部密度ρ(m)定义为
其中dc是截止距离,d(m,q)是点(r(m),P(m))和点(r(q),P(q))间的欧式距离。
χ(·)是一个截止函数,
选择具有不小于ρ0的局部密度的点组成集合Sd,其中ρ0为样本点的局部密度阈值,即
集合Sd中的样本点(r(m),P(m))在r-P平面上展示为具有一定密度的散点图,对于P轴,给定一个窗口宽度ε,确定在窗口宽度内火电机组调节速率r的最大值,通过滑动窗口ε得到K个r-P平面上轮廓散点,散点对应的样本值记为集合
对集合中的样本点(r(k),P(k))进行多项式拟合,Φ为所有次数不超过F(F≤K)的多项式构成的函数类,现求/>使得下式的W最小,其中af为次数为f的实发功率的多项式的系数,Pf为次数为f的实发功率。
而满足式(14)的称为最小二乘拟合多项式,其中,Pk为第k个样本点所对应的实发功率,/>为第k个样本点所对应的实发功率,次数为f,rk为第k个样本点所对应的速率。显然/>为a0,…,af,...,aF的多元函数,因此上述问题转化为求W=W(a0,…,af,...,aF)的极值问题。由多元函数求极值的必要条件,得
即
式(16)是关于a0,…,af,...,aF的线性方程组,其系数矩阵是一个对称正定矩阵,故有唯一解,可以解出af(f=0,...,F),进而得到火电机组调节速率包络线:
记r(P)为火电机组输出功率与变化速率关系式,并定义ri(P)为第i台机组输出功率与变化速率关系式。
步骤3:为获得最短电网AGC指令响应时间时电网AGC指令的分配结果,建立如下模型:
约束条件:
其中,T为电网AGC指令响应时间,数值上等于电网中调节时间最长的机组的负荷响应时间,即T=max{t1,t2,…,tI},ti=ΔPi/ri,其中ΔPi为分配给第i台火电机组的AGC指令调整量,ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
为处理上述约束,采用拉格朗日乘子法引入拉格朗日算子计算得到上述模型的最优解。记拉格朗日函数
对拉格朗日函数求偏导可得
由式(22)得
由式(23)和式(24)得
分配给第i台火电机组的AGC指令调整量为
根据式(17)所示的火电机组输出功率和调节速率的关系式,得到基于火电机组调节速率包络线的电网AGC指令分配结果。此时,第i台火电机组的AGC指令调整量为:
以下是本实施例的方法在具体示例中的应用。
以某区域电网为例,该区域含有4台火电机组,火电机组参数如表1所示。现电网需要增加100MW功率输出,即ΔP=100WM,计算分配到每台火电机组的调节负荷ΔPi,使电网能在最短时间内完成总负荷调整。
表1:火电机组参数信息表
/>
第一步,对火电机组实发功率的历史数据进行分段线性表示,将时间序列分成多个短直线段,计算每个数据段的调节速率,得到火电机组输出功率与该输出功率对应的调节速率的样本集合。以第一台火电机组为例,图3所示为火电机组输出功率与调节速率样本集合的散点图。
第二步,基于局部密度剔除样本集合中的处于边缘的样本点,如图3所示的五角星形数据点,筛选出具有一定局部密度的样本点,图4为具有一定局部密度的样本散点图。对筛选出的样本采用多项式拟合得到火电机组调节速率包络线,如图4所示,进而得到火电机组输出功率与调节速率的关系式。同理得到4台火电机组输出功率与调节速率的关系式,如表2所示。
表2:机组出力与最大调节速率关系式
机组编号 | 机组输出功率与调节速率的关系式 |
火电机组1 | r1=-9.21×10-4P1 2+0.44P1-44.72 |
火电机组2 | r2=-2.94×10-4P2 2+0.16P2-16.01 |
火电机组3 | r3=-0.0028P3 2+1.31P3-141.64 |
火电机组4 | r4=-0.0028P4 2+1.36P4-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 (14)
1.一种基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,包括以下步骤:
根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;
根据所述样本集合获得火电机组调节速率包络线;所述根据所述样本集合获得火电机组调节速率包络线的步骤,包括:基于局部密度剔除所述样本集合中的边缘样本点,对筛选后的样本集合采用多项式拟合获得火电机组调节速率包络线,确定火电机组输出功率与调节速率的关系式;
以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据火电机组调节速率包络线对电网AGC指令进行分配;
其中,基于局部密度剔除所述样本集合中的边缘样本点,对筛选后的样本集合采用多项式拟合获得火电机组调节速率包络线,确定火电机组输出功率与调节速率的关系式包括:
从局部密度剔除所述样本集合中选择具有足够局部密集的点,形成集合Sd;
样本点(r(m),P(m))的局部密度ρ(m)定义为:
其中,dc是截止距离,d(m,q)是点(r(m),P(m))和点(r(q),P(q))间的欧式距离χ(·)是一个截止函数;
选择具有不小于ρ0的局部密度的点组成集合Sd,其中ρ0为样本点的局部密度阈值,即:
集合Sd中的样本点(r(m),P(m))在r-P平面上展示为具有一定密度的散点图,对于P轴,给定一个窗口宽度ε,确定在窗口宽度内火电机组调节速率r的最大值,通过滑动窗口ε得到K个r-P平面上轮廓散点,散点对应的样本值记为集合
对集合中的样本点(r(k),P(k))进行多项式拟合,Φ为所有次数不超过F(F≤K)的多项式构成的函数类,求/>使得的W最小,其中af为次数为f的实发功率的多项式的系数,Pf为次数为f的实发功率;/>称为最小二乘拟合多项式,Pk为第k个样本点所对应的实发功率,/>为第k个样本点所对应的实发功率,次数为f,rk为第k个样本点所对应的速率;/>为a0,…,af,...,aF的多元函数;
将问题转化为求W=W(a0,…,af,...,aF)的极值问题;多元函数求极值的必要条件,得
即: 是关于a0,…,af,...,aF的线性方程组,其系数矩阵是一个对称正定矩阵,具有唯一解,解出af(f=0,...,F),得到火电机组调节速率包络线:
记r(P)为火电机组输出功率与变化速率关系式,并定义ri(P)为第i台机组输出功率与变化速率关系式。
2.如权利要求1所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,
所述根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合的步骤,包括:
将火电机组实发功率历史数据按时间序列划分成多个数据段,通过每个数据段的实发功率幅值变化量和持续时间的比值计算火电机组调节速率,即
其中,δ(m)为第m段数据的实发功率幅值变化量,τ(m)为第m段数据的持续时间;
火电机组调节速率r(m)对应的火电机组输出功率为y(nm),记P(m)=y(nm),根据分段结果得到火电机组输出功率与调节速率的样本集合
3.如权利要求1所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,
所述以最短电网AGC指令响应时间为优化目标,建立电网AGC指令优化分配模型的步骤,包括:建立电网AGC指令分配模型,如下:
以各台火电机组分配的AGC指令调整量之和等于电网AGC指令调整量为约束条件:
其中,T为电网AGC指令响应时间,ΔPi为分配给第i台火电机组的AGC指令调整量,ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
4.如权利要求3所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,
所述电网AGC指令响应时间T在数值上等于电网中调节时间最长的火电机组的负荷响应时间,T=max{t1,t2,…,tI},ti=ΔPi/ri。
5.如权利要求3所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,
所述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,还包括:
通过拉格朗日乘子法求解电网AGC指令分配模型,获得电网AGC指令优化分配模型,得到第i台火电机组分配的AGC指令调整量为:
ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
6.如权利要求5所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法,其特征在于,
所述根据火电机组调节速率包络线对电网AGC指令进行分配的步骤,获得第i台机组的AGC指令调整量为:
其中,ri(P)为第i台火电机组输出功率为P时的调节速率,ΔP为电网AGC指令调整量。
7.一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,包括:
数据处理模块,被配置为:根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合;
动态关系拟合模块,被配置为:根据所述样本集合获得火电机组调节速率包络线;所述根据所述样本集合获得火电机组调节速率包络线的步骤,包括:基于局部密度剔除所述样本集合中的边缘样本点,对筛选后的样本集合采用多项式拟合获得火电机组调节速率包络线,确定火电机组输出功率与调节速率的关系式;
指令分配模块,被配置为:以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型,根据火电机组调节速率包络线对电网AGC指令进行分配;
其中,基于局部密度剔除所述样本集合中的边缘样本点,对筛选后的样本集合采用多项式拟合获得火电机组调节速率包络线,确定火电机组输出功率与调节速率的关系式包括:
从局部密度剔除所述样本集合中选择具有足够局部密集的点,形成集合Sd;
样本点(r(m),P(m))的局部密度ρ(m)定义为:
其中,dc是截止距离,d(m,q)是点(r(m),P(m))和点(r(q),P(q))间的欧式距离χ(·)是一个截止函数;
选择具有不小于ρ0的局部密度的点组成集合Sd,其中ρ0为样本点的局部密度阈值,即:
集合Sd中的样本点(r(m),P(m))在r-P平面上展示为具有一定密度的散点图,对于P轴,给定一个窗口宽度ε,确定在窗口宽度内火电机组调节速率r的最大值,通过滑动窗口ε得到K个r-P平面上轮廓散点,散点对应的样本值记为集合
对集合中的样本点(r(k),P(k))进行多项式拟合,Φ为所有次数不超过F(F≤K)的多项式构成的函数类,求/>使得的W最小,其中af为次数为f的实发功率的多项式的系数,Pf为次数为f的实发功率;/>称为最小二乘拟合多项式,Pk为第k个样本点所对应的实发功率,/>为第k个样本点所对应的实发功率,次数为f,rk为第k个样本点所对应的速率;/>为a0,…,af,...,aF的多元函数;
将问题转化为求W=W(a0,…,af,...,aF)的极值问题;多元函数求极值的必要条件,得
即: 是关于a0,…,af,...,aF的线性方程组,其系数矩阵是一个对称正定矩阵,具有唯一解,解出af(f=0,...,F),得到火电机组调节速率包络线:
记r(P)为火电机组输出功率与变化速率关系式,并定义ri(P)为第i台机组输出功率与变化速率关系式。
8.如权利要求7所述的一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,
所述根据火电机组实发功率历史数据,获得火电机组输出功率与对应调节速率的样本集合的步骤,包括:
将火电机组实发功率历史数据按时间序列划分成多个数据段,通过每个数据段的实发功率幅值变化量和持续时间的比值计算火电机组调节速率,即
其中,δ(m)为第m段数据的实发功率幅值变化量,τ(m)为第m段数据的持续时间;
火电机组调节速率r(m)对应的火电机组输出功率为y(nm),记P(m)=y(nm),根据分段结果得到火电机组输出功率与调节速率的样本集合
9.如权利要求7所述的一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,
所述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,包括:
建立电网AGC指令分配模型,电网AGC指令分配模型如下:
以各台火电机组分配的AGC指令调整量之和等于电网AGC指令调整量为约束条件:
其中,T为电网AGC指令响应时间,ΔPi为分配给第i台火电机组的AGC指令调整量,ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
10.如权利要求9所述的一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,
所述电网AGC指令响应时间T在数值上等于电网中调节时间最长的火电机组的负荷响应时间,T=max{t1,t2,…,tI},ti=ΔPi/ri。
11.如权利要求9所述的一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,
所述以最短电网AGC指令响应时间为优化目标建立负荷指令优化分配模型的步骤,还包括:
通过拉格朗日乘子法求解电网AGC指令分配模型,获得电网AGC指令优化分配模型,得到第i台火电机组分配的AGC指令调整量为:
ri为第i台火电机组的调节速率,ΔP为电网AGC指令调整量。
12.如权利要求11所述的一种基于火电机组调节速率包络线的电网AGC指令优化分配系统,其特征在于,
所述根据火电机组调节速率包络线对电网AGC指令进行分配的步骤,获得第i台机组的AGC指令调整量为:
其中,ri(P)为第i台火电机组输出功率为P时的调节速率,ΔP为电网AGC指令调整量。
13.一种介质,其上存储有程序,其特征在于,该程序被处理器执行时实现如权利要求1-6任一项所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法中的步骤。
14.一种电子设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的程序,其特征在于,所述处理器执行所述程序时实现如权利要求1-6任一项所述的基于火电机组调节速率包络线的电网AGC指令优化分配方法中的步骤。
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 CN113595152A (zh) | 2021-11-02 |
CN113595152B true 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) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116454890B (zh) * | 2023-04-20 | 2024-02-06 | 中国南方电网有限责任公司 | 基于scuc模型的机组组合控制方法、装置和设备 |
Citations (3)
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 | 山东科技大学 | 一种基于机组爬坡率估计的功率指令优化分配方法 |
-
2021
- 2021-08-25 CN CN202110980828.8A patent/CN113595152B/zh active Active
Patent Citations (3)
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 | 山东科技大学 | 一种基于机组爬坡率估计的功率指令优化分配方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113595152A (zh) | 2021-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104808587B (zh) | 一种基于机加工设备运行状态的稼动率统计方法 | |
CN111092429B (zh) | 一种柔性互联配电网的优化调度方法、存储介质及处理器 | |
CN109347124B (zh) | 利用蓄热热泵群的电热微网联络线功率平抑方法与装置 | |
CN109583069B (zh) | 基于时间相关性的风速建模方法 | |
CN113595152B (zh) | 基于火电机组调节速率包络线的电网agc指令优化分配方法及系统 | |
CN105207253A (zh) | 考虑风电及频率不确定性的agc随机动态优化调度方法 | |
CN105225000B (zh) | 一种基于模糊序优化的风功率概率模型非参数核密度估计方法 | |
CN111553060B (zh) | 基于类噪声扰动的电网惯量常态化连续估计方法 | |
CN111798060A (zh) | 一种基于机组爬坡率估计的功率指令优化分配方法 | |
CN114844118A (zh) | 一种适用于微电网的多类型设备功率协调控制方法及系统 | |
CN110059938B (zh) | 一种基于关联规则驱动的配电网规划方法 | |
CN111985771B (zh) | 基于电网频率超欠调分析的电网频率调控方法及系统 | |
CN105956318A (zh) | 基于改进分裂 h-k 聚类方法的风电场机群划分方法 | |
CN115173489B (zh) | 一种基于二分法的火电集群调度方法及系统 | |
CN108988323B (zh) | 一种电力系统快速调频瓶颈边界场景的识别方法 | |
CN113536685A (zh) | 一种基于神经网络外推的风速概率模型的建模方法 | |
CN112613654A (zh) | 一种基于多类型储能的综合能源系统灵活性评估方法 | |
CN111900740B (zh) | 一种基于需求响应设备的电力系统调频方法和系统 | |
CN112560947A (zh) | 基于能源供需结构分析的聚类方法及装置 | |
CN114417918A (zh) | 一种风电场信号特征提取及去噪优化数据的方法 | |
CN106655189A (zh) | 基于虚拟发电机的最优潮流协调优化方法 | |
CN107609678B (zh) | 一种考虑矩信息不确定性的同型并行机生产调度方法 | |
CN107833149B (zh) | 一种基于多判别指标的电网发展时期划分方法及系统 | |
CN112131703A (zh) | 一种基于matlab功能模块的风机运行特性分析系统 | |
CN110707695B (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 |