首页 / 专利库 / 换热器 / 热电发电机 / 一种基于等微增率准则的多能虚拟电厂运行参数聚合方法

一种基于等微增率准则的多能虚拟电厂运行参数聚合方法

阅读:968发布:2020-05-12

专利汇可以提供一种基于等微增率准则的多能虚拟电厂运行参数聚合方法专利检索,专利查询,专利分析的服务。并且本 发明 涉及一种基于等微增率的多能虚拟电厂运行参数聚合方法,属于综合 能源 系统的运行控制技术领域。本发明方法将与电 力 系统相连的多能虚拟电厂中的发 电机 、 热电联产 机组、电负荷、热负荷等值为一台聚合发电机,该聚合发电机包含功率上下限及最大爬坡速率等运行参数,与电力系统通过联络线相连,便于电力系统调度,为电力系统提供额外的灵活性。本方法基于等微增率准则,尽可能减小了灵活性损失,且计算速度较快,便于多能虚拟电厂使用。,下面是一种基于等微增率准则的多能虚拟电厂运行参数聚合方法专利的具体信息内容。

1.一种基于等微增率准则的多能虚拟电厂运行参数聚合方法,其特征在于该方法包括以下步骤:
(1)设定多能虚拟电厂包含k台热电联产机组、m台发电机、电负荷和热负荷,将发电机、热电联产机组、电负荷和热负荷等值为一台聚合发电机,该聚合发电机与电系统通过联络线相连,根据等微增率准则,求解以下联立方程,分别得到多能虚拟电厂中待聚合的各热电联产机组的最低有功功率
cchp_j=αj·(pchp_j)2+βi·pchp_j+γj
其中,k表示多能虚拟电厂中待聚合的热电联产机组总台数,cchp_j表示第j台热电联产机组的成本函数,αj、βj、γj为第j台热电联产机组成本函数的二次项、一次项及常数项系数,ρj为第j台热电联产机组的热电比,上述热电联产机组的参数可从热电联产机组说明书或根据实际运行数据获取; 表示第j台热电联产机组在调度时刻t的供热功率,表示第j台热电联产机组在调度时刻t的电出力, 表示多能虚拟电厂在调度时刻t的热负荷功率,根据热负荷的历史运行数据预测得到;Υ为所有调度时刻构成的集合;
(2)将多能虚拟电厂中的k台热电联产机组、m台发电机、电负荷和热负荷聚合成一个发电机,该聚合发电机的运行参数的计算方法如下:
(2-1)计算多能虚拟电厂中的发电机在调度时刻t的有功功率上限:
(2-1-1)置发电机编号i=1,调度时刻t=1;
(2-1-2)将发电机i在调度时刻t的向上爬坡能力与该发电机的有功功率上限PG_i_max进行比较,若 则得到调度时刻t发电机i的有功功率上限
转步骤(2-1-3);若 则得到调度时刻t发电机
i的有功功率上限 转步骤(2-1-3),其中,rampG_i为第i台发电机
的爬坡速率,PG_i_max和rampG_i均可由发电机的说明书获取, 为发电机i的初始发电功率,根据发电机上一日调度计划给出;
(2-1-3)将调度时刻t与总调度时刻数T进行比较,若t(2-1-4)将发电机编号i与待聚合发电机总台数m进行比较,若i(2-2)计算多能虚拟电厂中的每个调度时刻t热电联产机组的有功功率上限:
(2-2-1)置热电联产机组编号j=1,调度时刻t=1;
(2-2-2)将热电联产机组j在调度时刻t的向上爬坡能力与该热电联产机组的有功功率上限Pchp_j_max进行比较,若 则得到调度时刻t热电联产机组j
的有功功率上限 转步骤(2-2-3);若 则得
到调度时刻t热电联产机组j的有功功率上限 转步骤(2-2-3),
其中,rampCHP_j为第j台热电联产机组的爬坡速率,Pchp_j_max和rampCHP_j由热电联产机组说明书获取, 为热电联产机组j的初始有功功率,根据热电联产机组上一日调度计划给出;
(2-2-3)将调度时刻t与总调度时刻数T进行比较,若t(2-2-4)将热电联产机组编号j与待聚合热电联产机组总台数k比较,若j(2-3)计算多能虚拟电厂中的每个调度时刻t发电机的有功功率下限:
(2-3-1)置发电机编号i=1,调度时刻t=1;
(2-3-2)将发电机i在调度时刻t的向下爬坡能力与该发电机的有功功率下限PG_i_min进行比较,若 则得到调度时刻t发电机i的有功功率下限PtG_i_min=
PG_i_min,并转步骤(2-3-3),若 则得到调度时刻t发电机i的有功
功率下限 并转步骤(2-3-3),其中,发电机的有功功率下限PG
_i_min由发电机说明书获取;
(2-3-3)将调度时刻t与总调度时刻数T进行比较,若t若t≥T,转步骤(2-3-4);
(2-3-4)将发电机编号i与待聚合发电机总台数m比较,若i(2-4)计算多能虚拟电厂中的每个调度时刻t热电联产机组的有功功率下限:
(2-4-1)置热电联产机组编号j=1,调度时刻t=1;
(2-4-2)将热电联产机组j在调度时刻t的向下爬坡能力与步骤(1)的热电联产机组j的有功功率下限 进行比较,若 则得到调度时刻t热电
联产机组j的有功功率下限 并转步骤(2-4-3);若
则得到调度时刻t热电联产机组j的有功功率下限Ptchp_j_min,
即 并转步骤(2-4-3);
(2-4-3)将调度时刻t与总调度时刻数T进行比较,若t(2-4-4)将热电联产机组编号j与待聚合热电联产机组总台数k比较,若j(2-5)利用下式,计算得到每个调度时刻t聚合发电机的有功功率上限PtTG_max、有功功率下限PtTG_min及最大爬坡速率rampTG,作为聚合发电机的聚合参数,实现多能虚拟电厂运行参数的聚合:
其中, 表示调度时刻t的电负荷功率,电负荷功率可根据历史负荷数据进行预测得到。

说明书全文

一种基于等微增率准则的多能虚拟电厂运行参数聚合方法

技术领域

[0001] 本发明涉及一种基于等微增率的多能虚拟电厂运行参数聚合方法,属于综合能源系统的运行控制技术领域。

背景技术

[0002] 为了提高多种能源的利用效率、充分利用各种能源的灵活性,减少排放量并提高新能源的渗透率,能源互联网及综合能源系统的概念应运而生。热电联产技术可充分利用燃料中的能量,在发电的同时利用余热供热,在工业生产、城市供热等领域有广泛的应用。
[0003] 近年来电系统中新能源发电的比例逐年上升,新能源发电随机性、间歇性、波动性的特点为电力系统的运行带来了新的挑战,电力系统亟需探索新的灵活性以减少新能源随机性带来的冲击。很多生产工艺同时需要电能热能的工厂均装有热电联产机组,这些机组难以由电力系统调度中心集中控制。
[0004] 聚合技术可在考虑电、热负荷平衡的同时将一定地理范围内的工厂中的自备小型发电机、分布式光伏发电装置、热电联产机组及电、热负荷聚合为一个发电机,该聚合发电机包含有功功率上下限及最大爬坡速率等运行参数,可将这些运行参数提供给电力系统调度中心,电力系统调度中心可根据该模型下发调度指令,再通过解聚合技术将电力系统调度指令分配给各机组,完成对电力系统调度计划的跟随,从而使电力系统可以调控这些机组,为电力系统提供了新的灵活性。例如专利申请号(或公开号)为201910092126.9、发明名称为“一种基于最大内嵌正方体的虚拟电厂灵活性聚合方法”的中国专利申请,将分布式光伏发电装置、电池及热负荷聚合形成虚拟电厂,并求解其最大内接正方体,以得到虚拟电厂的输出功率范围,用于电力系统调度。但该聚合方法未考虑虚拟电厂最大、最小有功功率可能随时间变化而变化的特点,因此该聚合方法过于保守,损失了一定的灵活性。

发明内容

[0005] 本发明的目的是提出一种基于等微增率准则的多能虚拟电厂运行参数聚合方法,针对已有技术中存在的最大、最小值不随时间变化或计算时间过长的问题,对聚合方法进行改进,以使虚拟电厂向电力系统提供尽可能多的灵活性。
[0006] 本发明提出的基于等微增率准则的多能虚拟电厂运行参数聚合方法,包括以下步骤:
[0007] (1)设定多能虚拟电厂包含k台热电联产机组、m台发电机、电负荷和热负荷,将发电机、热电联产机组、电负荷和热负荷等值为一台聚合发电机,该聚合发电机与电力系统通过联络线相连,根据等微增率准则,求解以下联立方程,分别得到多能虚拟电厂中待聚合的各热电联产机组的最低有功功率
[0008]
[0009]
[0010] cchp_j=αj·(pchp_j)2+βi·pchp_j+γj
[0011]
[0012] 其中,k表示多能虚拟电厂中待聚合的热电联产机组总台数,cchp_j表示第j台热电联产机组的成本函数,αj、βj、γj为第j台热电联产机组成本函数的二次项、一次项及常数项系数,ρj为第j台热电联产机组的热电比,上述热电联产机组的参数可从热电联产机组说明书或根据实际运行数据获取; 表示第j台热电联产机组在调度时刻t的供热功率,表示第j台热电联产机组在调度时刻t的电出力, 表示多能虚拟电厂在调度时刻t的热负荷功率,根据热负荷的历史运行数据预测得到;Υ为所有调度时刻构成的集合;
[0013] (2)将多能虚拟电厂中的k台热电联产机组、m台发电机、电负荷和热负荷聚合成一个发电机,该聚合发电机的运行参数的计算方法如下:
[0014] (2-1)计算多能虚拟电厂中的发电机在调度时刻t的有功功率上限:
[0015] (2-1-1)置发电机编号i=1,调度时刻t=1;
[0016] (2-1-2)将发电机i在调度时刻t的向上爬坡能力与该发电机的有功功率上限PG_i_max进行比较,若 则得到调度时刻t发电机i的有功功率上限转步骤(2-1-3);若 则得到调度时刻t发电机
i的有功功率上限 转步骤(2-1-3),其中,rampG_i为第i台发电机
的爬坡速率,PG_i_max和rampG_i均可由发电机的说明书获取, 为发电机i的初始发电功率,根据发电机上一日调度计划给出;
[0017] (2-1-3)将调度时刻t与总调度时刻数T进行比较,若t
[0018] (2-1-4)将发电机编号i与待聚合发电机总台数m进行比较,若i
[0019] (2-2)计算多能虚拟电厂中的每个调度时刻t热电联产机组的有功功率上限:
[0020] (2-2-1)置热电联产机组编号j=1,调度时刻t=1;
[0021] (2-2-2)将热电联产机组j在调度时刻t的向上爬坡能力与该热电联产机组的有功chp_j_max功率上限P 进行比较,若 则得到调度时刻t热电联产机
组j的有功功率上限 转步骤(2-2-3);若
则得到调度时刻t热电联产机组j的有功功率上限 转步骤(2-
2-3),其中,rampCHP_j为第j台热电联产机组的爬坡速率,Pchp_j_max和rampCHP_j由热电联产机组说明书获取, 为热电联产机组j的初始有功功率,根据热电联产机组上一日调度计划给出;
[0022] (2-2-3)将调度时刻t与总调度时刻数T进行比较,若t
[0023] (2-2-4)将热电联产机组编号j与待聚合热电联产机组总台数k比较,若j
[0024] (2-3)计算多能虚拟电厂中的每个调度时刻t发电机的有功功率下限:
[0025] (2-3-1)置发电机编号i=1,调度时刻t=1;
[0026] (2-3-2)将发电机i在调度时刻t的向下爬坡能力与该发电机的有功功率下限PG_i_min进行比较,若 则得到调度时刻t发电机i的有功功率下限并转步骤(2-3-3),若 则得到调度时刻t发电
机i的有功功率下限 并转步骤(2-3-3),其中,发电机的有功功率
下限PG_i_min由发电机说明书获取;
[0027] (2-3-3)将调度时刻t与总调度时刻数T进行比较,若t
[0028] (2-3-4)将发电机编号i与待聚合发电机总台数m比较,若i
[0029] (2-4)计算多能虚拟电厂中的每个调度时刻t热电联产机组的有功功率下限:
[0030] (2-4-1)置热电联产机组编号j=1,调度时刻t=1;
[0031] (2-4-2)将热电联产机组j在调度时刻t的向下爬坡能力与步骤(1)的热电联产机组j的有功功率下限 进行比较,若 则得到调度时刻t热电联产机组j的有功功率下限 并转步骤(2-4-3);若
则得到调度时刻t热电联产机组j的有功功率下限Ptchp_j_min,
即 并转步骤(2-4-3);
[0032] (2-4-3)将调度时刻t与总调度时刻数T进行比较,若t
[0033] (2-4-4)将热电联产机组编号j与待聚合热电联产机组总台数k比较,若j
[0034] (2-5)利用下式,计算得到每个调度时刻t聚合发电机的有功功率上限PtTG_max、有功功率下限PtTG_min及最大爬坡速率rampTG,作为聚合发电机的聚合参数,实现多能虚拟电厂运行参数的聚合:
[0035]
[0036]
[0037]
[0038] 其中, 表示调度时刻t的电负荷功率,电负荷功率可根据历史负荷数据进行预测得到。
[0039] 本发明提出的基于等微增率准则的多能虚拟电厂运行参数聚合方法,其特点和效果是:
[0040] 本发明的多能虚拟电厂运行参数聚合方法,考虑了电负荷平衡及热负荷平衡,将工厂中的热电联产机组、发电机、电负荷及热负荷聚合成为一台包含功率上下限及最大爬坡速率等运行参数的聚合发电机,便于电力系统调度,为电力系统提供额外的灵活性。本方法基于等微增率准则,尽可能减小了灵活性损失,且计算速度较快,便于多能虚拟电厂使用。附图说明
[0041] 图1为本发明方法涉及的电力系统与多能虚拟电厂连接关系示意图

具体实施方式

[0042] 本发明提出的基于等微增率准则的多能虚拟电厂运行参数聚合方法,包括以下步骤:
[0043] (1)设定多能虚拟电厂包含k台热电联产机组、m台发电机、电负荷和热负荷,将发电机、热电联产机组、电负荷和热负荷等值为一台聚合发电机,该聚合发电机与电力系统通过联络线相连,如图1所示;根据等微增率准则,求解以下联立方程,分别得到多能虚拟电厂中待聚合的各热电联产机组的最低有功功率
[0044]
[0045]
[0046] cchp_j=αj·(pchp_j)2+βi·pchp_j+γj
[0047]
[0048] 其中,k表示多能虚拟电厂中待聚合的热电联产机组总台数,cchp_j表示第j台热电联产机组的成本函数,αj、βj、γj为第j台热电联产机组成本函数的二次项、一次项及常数项系数,ρj为第j台热电联产机组的热电比,上述热电联产机组的参数可从热电联产机组说明书或根据实际运行数据获取; 表示第j台热电联产机组在调度时刻t的供热功率,表示第j台热电联产机组在调度时刻t的电出力, 表示多能虚拟电厂在调度时刻t的热负荷功率,根据热负荷的历史运行数据预测得到;Υ为所有调度时刻构成的集合;
[0049] 上述联立方程的求解过程如下:
[0050] 将第三式代入第一式,可得:
[0051]
[0052] 将第四式代入上式,可得:
[0053]
[0054] 将第二式与上式联立,求得 的解为则可得到考虑热负荷平衡后修正的热电联产机组的最小有功功率为:
[0055]
[0056] 对于其他形式(非二次函数)的成本函数,也可通过类似的推导方法联立第一式与第二式,求解出对应的热出力,再由热电出力之间的关系转换得到修正后热电联产机组的最小有功功率;
[0057] (2)将多能虚拟电厂中的k台热电联产机组、m台发电机、电负荷和热负荷聚合成一个发电机,该聚合发电机的运行参数(包括最大爬坡速率及有功功率上限和下限)的计算方法如下:
[0058] (2-1)计算多能虚拟电厂中的发电机在调度时刻t的有功功率上限:
[0059] (2-1-1)置发电机编号i=1,调度时刻t=1;
[0060] (2-1-2)将发电机i在调度时刻t的向上爬坡能力与该发电机的有功功率上限PG_i_max进行比较,若 则得到调度时刻t发电机i的有功功率上限转步骤(2-1-3);若 则得到调度时刻t发电机
i的有功功率上限 转步骤(2-1-3),其中,rampG_i为第i台发电机
的爬坡速率,PG_i_max和rampG_i均可由发电机的说明书获取, 为发电机i的初始发电功率,根据发电机上一日调度计划给出;
[0061] (2-1-3)将调度时刻t与总调度时刻数T进行比较,若t
[0062] (2-1-4)将发电机编号i与待聚合发电机总台数m进行比较,若i
[0063] (2-2)计算多能虚拟电厂中的每个调度时刻t热电联产机组的有功功率上限:
[0064] (2-2-1)置热电联产机组编号j=1,调度时刻t=1;
[0065] (2-2-2)将热电联产机组j在调度时刻t的向上爬坡能力与该热电联产机组的有功功率上限Pchp_j_max进行比较,若 则得到调度时刻t热电联产机组j的有功功率上限 转步骤(2-2-3);若
则得到调度时刻t热电联产机组j的有功功率上限 转步骤(2-
2-3),其中,rampCHP_j为第j台热电联产机组的爬坡速率,Pchp_j_max和rampCHP_j由热电联产机组说明书获取, 为热电联产机组j的初始有功功率,根据热电联产机组上一日调度计划给出;
[0066] (2-2-3)将调度时刻t与总调度时刻数T进行比较,若t
[0067] (2-2-4)将热电联产机组编号j与待聚合热电联产机组总台数k比较,若j
[0068] (2-3)计算多能虚拟电厂中的每个调度时刻t发电机的有功功率下限:
[0069] (2-3-1)置发电机编号i=1,调度时刻t=1;
[0070] (2-3-2)将发电机i在调度时刻t的向下爬坡能力与该发电机的有功功率下限PG_i_min进行比较,若 则得到调度时刻t发电机i的有功功率下限PtG_i_min=PG_i_min,并转步骤(2-3-3),若 则得到调度时刻t发电机i
的有功功率下限 并转步骤(2-3-3),其中,发电机的有功功率下
G_i_min
限P 由发电机说明书获取;
[0071] (2-3-3)将调度时刻t与总调度时刻数T进行比较,若t
[0072] (2-3-4)将发电机编号i与待聚合发电机总台数m比较,若i
[0073] (2-4)计算多能虚拟电厂中的每个调度时刻t热电联产机组的有功功率下限:
[0074] (2-4-1)置热电联产机组编号j=1,调度时刻t=1;
[0075] (2-4-2)将热电联产机组j在调度时刻t的向下爬坡能力与步骤(1)的热电联产机组j的有功功率下限 进行比较,若 则得到调度时刻t热电联产机组j的有功功率下限 并转步骤(2-4-3);若
则得到调度时刻t热电联产机组j的有功功率下限Ptchp_j_min,
即 并转步骤(2-4-3);
[0076] (2-4-3)将调度时刻t与总调度时刻数T进行比较,若t
[0077] (2-4-4)将热电联产机组编号j与待聚合热电联产机组总台数k比较,若j
[0078] (2-5)利用下式,计算得到每个调度时刻t聚合发电机的有功功率上限PtTG_max、有TG_min TG功功率下限Pt 及最大爬坡速率ramp ,作为聚合发电机的聚合参数,实现多能虚拟电厂运行参数的聚合:
[0079]
[0080]
[0081]
[0082] 其中, 表示调度时刻t的电负荷功率,电负荷功率可根据历史负荷数据进行预测得到。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈