CN106992555B - 一种辅助服务市场环境下的热电联合系统分布式调度方法 - Google Patents
一种辅助服务市场环境下的热电联合系统分布式调度方法 Download PDFInfo
- Publication number
- CN106992555B CN106992555B CN201710274062.5A CN201710274062A CN106992555B CN 106992555 B CN106992555 B CN 106992555B CN 201710274062 A CN201710274062 A CN 201710274062A CN 106992555 B CN106992555 B CN 106992555B
- Authority
- CN
- China
- Prior art keywords
- power
- heat
- node
- power system
- scheduling
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 230000012447 hatching Effects 0.000 title abstract 6
- 238000010438 heat treatment Methods 0.000 claims abstract description 77
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 120
- 125000004432 carbon atom Chemical group C* 0.000 claims description 39
- 238000004519 manufacturing process Methods 0.000 claims description 39
- 230000009194 climbing Effects 0.000 claims description 22
- 238000010248 power generation Methods 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000008878 coupling Effects 0.000 claims description 7
- 238000010168 coupling process Methods 0.000 claims description 7
- 238000005859 coupling reaction Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000020169 heat generation Effects 0.000 claims description 6
- 238000009987 spinning Methods 0.000 claims description 6
- 230000003190 augmentative effect Effects 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 229910052799 carbon Inorganic materials 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 230000001105 regulatory effect Effects 0.000 claims 1
- 230000033228 biological regulation Effects 0.000 abstract description 6
- 230000002452 interceptive effect Effects 0.000 abstract description 3
- 230000000694 effects Effects 0.000 abstract description 2
- 238000009825 accumulation Methods 0.000 abstract 1
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000005338 heat storage Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000026676 system process 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
-
- 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
本发明涉及一种辅助服务市场环境下的热电联合系统分布式调度方法,属于电力系统运行控制技术领域。本方法综合考虑了电力系统的调度模型和供热系统的调度模型,并考虑了辅助服务市场环境下的额外调度目标,建立了辅助服务市场环境下热电联合系统调度模型。基于变方向乘子算法,电力系统和各个供热系统分布式处理各自所管辖区域的调度问题,通过不同系统间边界数据的交互迭代获得辅助服务市场环境下的热电联合系统分布式调度最优解。本发明提出的辅助服务市场环境下的热电联合系统分布式调度方法具有良好的收敛速度,可以通过供热系统的储热效应提高系统运行的灵活性。
Description
技术领域
本发明涉及一种辅助服务市场环境下的热电联合系统分布式调度方法,属于电力系统运行控制技术领域。
背景技术
在我国三北地区风电弃电现象严重,在冬天受到城区供热系统的最小运行约束限制,热电联产机组无法减小其出力为风电的消纳提供空间。具体来说,冬季夜晚属于风电高发期,但夜晚同时具有很高的居民热负荷,导致了风电需要被限制来为热电联产机组提供出力空间,因此供热和风电消纳需求之间存在矛盾。由于供热管网具有良好的储热效益,综合考虑电力系统和供热系统进行调度可以实现负荷调峰,有效提高系统运行的灵活性,减小弃风现象。
在辅助服务市场环境下,常规机组的有功备用需要被考虑到备用成本中,这导致了在辅助服务市场环境下电力系统的调度与传统调度存在不同,需要针对辅助服务市场独立考虑。
电力系统和供热系统分别由电力公司与供热公司独立运行,由于调度独立性,难以实现对电力系统和供热系统完整模型进行统一调度。因此,提出辅助服务市场环境下热电联合系统分布式调度算法方法十分必要。具体来说,电力系统和供热系统可以分别对各自管辖的区域进行调度优化,通过边界条件的迭代获得辅助服务市场环境下热电联合系统调度的全局最优解。
变方向乘子算法是一种分布式求解算法,可将复杂的数学优化问题转化为若干相对简单的问题分布式进行求解,通过迭代可获得全局最优解。
发明内容
本发明的目的是提出一种辅助服务市场环境下的热电联合系统分布式调度方法,针对电力系统与供热系统分别建模,考虑电力系统和供热系统的耦合约束,建立辅助服务市场环境下的热电联合调度模型。针对辅助服务市场环境下的热电联合调度模型,基于变方向乘子算法,提出辅助服务市场环境下的热电联合调度模型的分布式求解算法。
本发明提出的辅助服务市场环境下的热电联合系统分布式调度方法,包括以下步骤:
(1)建立辅助服务市场环境下热电联合系统调度模型,调度模型由目标函数和约束条件构成,具体包括:
(1-1)辅助服务市场环境下热电联合系统调度模型的目标函数:
热电联合系统调度模型的目标函数为电力系统运行成本、电力系统备用辅助服务成本与供热系统运行成本之和的最小化,表达式为:
上式中,T为调度时段集合,上标CHP、TU、WD、HB和AS分别表示热电联产机组、常规机组、风电机组、热锅炉和辅助服务,ICHP、ITU、IWD和IHB分别为热电联合系统中的热电联产机组、常规机组、风电机组和热锅炉的编号集合,下标i,t表示第i台机组在调度时段t的相应变量,和分别为热电联合系统中热电联产机组、常规机组、风电机组和热锅炉的生产成本函数,为电力系统中常规机组提供的备用辅助服务成本函数,i为热电联产机组、常规机组、风电机组和热锅炉中机组或锅炉的编号,t为调度时段;
其中的热电联产机组的生产成本函数为:
上式中,和分别为生产成本系数,生产成本系数为热电联产机组固有参数,和分别为第i台热电联产机组在第t个调度时段的有功功率和产热功率;
常规机组的生产成本函数为:
上式中,和为常规机组的发电成本系数,该系数为常规机组的固有参数,为第i台常规机组在第t个调度时段的有功功率;
风电机组的生产成本函数为:
上式中,为弃风惩罚因子,取值范围一般为1~10,弃风惩罚因子的取值根据对风电的消纳需求确定,由电力系统调度中心根据调度结果反馈进行调节,为第i 台风电机组在第t个调度时段的可用有功功率,从风电的发电预测数据中获取,为第 i台风电机组在第t个调度时段的实际有功功率;
常规机组提供的备用辅助服务成本函数为:
上式中,上标RU和RD分别表示向上和向下旋转备用,和分别为电力系统在第t个调度时段的向上和向下的备用辅助服务竞价价格,该价格由辅助服务市场的参与主体中报价最高者处获取,和分别为第i台常规机组在第t个调度时段的向上旋转备用和向下旋转备用;
热锅炉的生产成本函数为:
上式中,为热锅炉的产热成本系数,为热锅炉固有参数,为第i台热锅炉在第t个调度时段的产热功率;
(1-2)辅助服务市场环境下热电联合系统调度模型的约束条件,包括:
(1-2-1)辅助服务市场环境下热电联合系统中电力系统的运行约束条件:
热电联产机组运行的有功功率和产热功率的约束条件为:
上式中,γ为运行极点的编号,运行极点是指热电联产机组的产热功率和发电功率所组成的平面区域的顶点,NEi为第i台热电联产机组的运行极点编号集合,Pi γ,分别为第i台热电联产机组第γ个运行极点的有功功率和产热功率,为第i台热电联产机组在第t个调度时段的运行点对应第γ个运行极点的凸组合系数,取值范围为0~1;
热电联产机组爬坡能力约束条件为:
上式中,和分别为第i台热电联产机组的向上爬坡能力和向下爬坡能力,从机组说明书中获取,ΔT为每个调度时段的长度;
常规机组运行的有功功率约束条件为:
上式中,和分别为第i台常规机组的有功功率上限和有功功率下限;
常规机组爬坡能力的约束条件为:
上式中,和分别为第i台常规机组的向上爬坡能力和向下爬坡能力;
常规机组的旋转备用约束条件为:
上式中,和分别为第i台常规机组在第t个调度时段的向上旋转备用和向下旋转备用;
风电机组运行约束条件为:
上式中,为第i台风电机组在第t个调度时段的可用有功功率,从风电的发电预测数据中获取,为第i台风电机组在第t个调度时段的实际有功功率;
热电联合系统中电力系统的功率平衡约束条件为:
上式中,m为电力系统中负荷的编号,上标LD代表负荷,ILD为电力系统所有负荷的编号集合,Dm,t为第m个负荷在第t个调度时段的负荷大小;
热电联合系统中电力系统的线路传输容量约束条件为:
上式中,j为电力系统中线路的编号,ILN代表电力系统线路编号集合,l为电力系统中节点的编号,IEPS代表电力系统中节点编号集合,和分别代表电力系统中与第l个节点相连接的热电联产机组、常规机组、风电机组与负荷的编号集合,Lj代表电力系统中第j个线路的最大有功功率,SFj-l为第l个节点到第j个线路的转移分布因子,转移分布因子从电力系统运行拓扑图上获取;
电力系统的旋转备用约束条件为:
上式中,SRUt和SRDt分别代表电力系统在第t个调度时段所需的向上旋转备用和向下旋转备用,和分别为第i台常规机组在第t个调度时段的向上旋转备用和向下旋转备用;
(1-2-2)辅助服务市场环境下热电联合系统中供热系统运行约束条件,包括:
(1-2-2-1)包括有热电联产机组和热锅炉的热源供热约束条件,包括:
热源产热功率与供热系统节点的供水和回水之间温差的约束条件为:
上式中,集合和分别为与供热系统节点k相连的热电联产机组编号集合和热锅炉编号集合,C为水的比热容,Mk为供热系统节点k的水流量,上标S和R分别表示供水管和回水管,下标k,t表示供热系统中第k个节点在调度时段t的相应变量,和分别表示供热系统节点k在调度时段t的供水管水温和回水管水温,集合为供热系统中与热源相连接的节点集合;
热锅炉的产热功率的约束条件,即热锅炉的产热功率需在其产热功率上限内:
上式中,表示第i个热锅炉的产热功率上限;
供热系统中与热源相连接的节点的供水管水温约束条件:
上式中,和为供热系统中任意节点k的供水管水温上限和供水管水温下限;
(1-2-2-2)换热站运行约束条件,包括:
换热功率与供热系统节点的供水管水温、回水管水温的温差为:
上式中,n为换热站的编号,为供热系统中与节点k相连接的换热站编号集合,表示第n个换热站在第t个调度时段的换热功率,表示供热系统中与换热站相连接的节点集合;
换热站的回水管水温确保在安全范围内:
上式中,和分别为供热系统中节点k的回水管水温上限和回水管水温下限;
(1-2-2-3)供热系统运行约束条件,包括:
其中,k1和k2分别为供热系统中相邻的两个节点,定义k1为流入节点,k2为流出节点,和分别为供热系统中由节点k2流向节点k1的供水流量和回水流量,和分别为供热系统中节点k1在供水管中的父节点集合和在回水管中的父节点集合,为第t个调度时段的环境温度,和分别为供热系统中由节点 k2流向节点k1的供水热传导系数和回水热传导系数,和的取值由下式测算:
其中,和分别为供热系统中由节点k2流向节点k1的供水管和回水管的单位热传导系数,和分别为供热系统中由节点k2流向节点k1的供水管和回水管的长度;
和分别为供热系统的温度变量,温度变量的计算公式为:
上式中,分别表示供热系统中由节点k2流向节点k1时供水管和回水管所需的调度时段数,符号表示向下取整;
(2)将步骤(1)建立的辅助服务市场环境下热电联合系统调度模型转化成矩阵形式,用xE表示电力系统变量,电力系统变量包括和用xH表示供热系统变量,供热系统变量包括和则辅助服务市场环境下热电联合系统调度模型转化为如下的矩阵形式:
约束条件为:
xE∈XE
xH∈XH
DxE=ExH
上式中,CE和CH分别表示电力系统和供热系统的目标函数,其中,CE为CH为约束集合XE为电力系统的约束条件集合,即步骤(1-2-1)中的全部约束条件所构成的集合;约束集合XH表示供热系统约束条件集合,即 步骤(1-2-2)中的全部约束条件所构成的集合;约束 DxE=ExH表示电力系统和供热系统耦合约束条件,既变量xE中的和变量xH中的对应相等,D、E的每一行与每一个电力系统和供热系统耦合约束条件一一对应, D的每一列与电力系统中的每一个变量一一对应,E的每一列与供热系统中的每一个变量一一对应,其中D、E的每一个元素为该元素所在行所对应的约束条件中其列所代表的变量的系数;
(3)采用变方向乘子算法,对上述步骤(2)得到的矩阵形式的辅助服务市场环境下热电联合系统调度模型进行求解,步骤如下:
(3-1)初始化:将迭代次数u初始化为0,初始化λ(u)=0,定义α=1,并定义如下式的增广拉格朗日函数A(xE,xH,λ):
(3-2)按照下式求解下面的电力系统的调度策略:
(3-3)按照下式求解下面的供热系统的调度策略:
(3-4)按照下式更新λ(u+1):
(3-5)根据上述步骤(3-2)和步骤(3-3)的计算结果进行判断,若||DxE-ExH||<10-6,则进行步骤(3-6),若||DxE-ExH||≥10-6,则将迭代次数u增加1,返回步骤(3-2)继续计算;
(3-6)将得到的最优解作为辅助服务市场环境下热电联合系统调度的调度参数。
本发明提出的辅助服务市场环境下的热电联合系统分布式调度方法,其优点是:
本方法结合了备用辅助服务市场环境,综合考虑了电力系统的调度模型和供热系统的调度模型,并考虑了辅助服务市场环境下的额外调度目标,建立了辅助服务市场环境下热电联合系统调度模型。针对所提出的辅助服务市场环境下热电联合系统调度模型,利用变方向乘子算法,实现了辅助服务市场环境下的热电联合系统分布式调度。本发明提出的辅助服务市场环境下的热电联合系统分布式调度方法中,电力系统和各个供热系统分布式处理各自所管辖区域的调度问题,通过不同系统间边界数据的交互迭代获得辅助服务市场环境下的热电联合系统分布式调度最优解。所提出的辅助服务市场环境下的热电联合系统分布式调度方法具有良好的收敛速度,可以通过供热系统的储热效应提高系统运行的灵活性。本发明方法中的迭代算法可以通过电力系统与供热系统之间边界条件的交互迭代获得全局最优解,并能显著提升供热系统的灵活性。
附图说明
图1是本发明方法涉及的热电联合系统结构示意图。
具体实施方式
本发明提出的辅助服务市场环境下的热电联合系统分布式调度方法,包括以下步骤:
(1)建立辅助服务市场环境下热电联合系统调度模型,调度模型由目标函数和约束条件构成,具体包括:
(1-1)辅助服务市场环境下热电联合系统调度模型的目标函数:
热电联合系统调度模型的目标函数为电力系统运行成本、电力系统备用辅助服务成本与供热系统运行成本之和的最小化,表达式为:
上式中,T为调度时段集合,上标CHP、TU、WD、HB和AS分别表示热电联产机组、常规机组、风电机组、热锅炉和辅助服务,ICHP、ITU、IWD和IHB分别为热电联合系统中的热电联产机组、常规机组、风电机组和热锅炉的编号集合,下标i,t表示第i台机组在调度时段t的相应变量,和分别为热电联合系统中热电联产机组、常规机组、风电机组和热锅炉的生产成本函数,为电力系统中常规机组提供的备用辅助服务成本函数,i为热电联产机组、常规机组、风电机组和热锅炉中机组或锅炉的编号,t为调度时段;
其中的热电联产机组的生产成本函数为:
上式中,和分别为生产成本系数,生产成本系数为热电联产机组固有参数,该固有参数可以从机组说明书上获取,和分别为第i台热电联产机组在第t个调度时段的有功功率和产热功率;
常规机组的生产成本函数为:
上式中,和为常规机组的发电成本系数,该系数为常规机组的固有参数,该固有参数可以从机组说明书上获取,,为第i台常规机组在第t个调度时段的有功功率;
风电机组的生产成本函数为:
上式中,为弃风惩罚因子,取值范围一般为1~10,弃风惩罚因子的取值根据对风电的消纳需求确定,由电力系统调度中心根据调度结果反馈进行调节,为第i台风电机组在第t个调度时段的可用有功功率,从风电的发电预测数据中获取,为第 i台风电机组在第t个调度时段的实际有功功率;
常规机组提供的备用辅助服务成本函数为:
上式中,上标RU和RD分别表示向上和向下旋转备用,和分别为电力系统在第t个调度时段的向上和向下的备用辅助服务竞价价格,该价格由辅助服务市场的参与主体中报价最高者处获取,和分别为第i台常规机组在第t个调度时段的向上旋转备用和向下旋转备用,其中,向上旋转备用和向下旋转备用分别指发电机组所能够提供的向上发电功率调节范围与向下发电功率调节范围;
热锅炉的生产成本函数为:
上式中,为热锅炉的产热成本系数,为热锅炉固有参数,该固有参数可以从热锅炉说明书上获取,为第i台热锅炉在第t个调度时段的产热功率;
(1-2)辅助服务市场环境下热电联合系统调度模型的约束条件,包括:
(1-2-1)辅助服务市场环境下热电联合系统中电力系统的运行约束条件:
热电联产机组运行的有功功率和产热功率的约束条件为:
上式中,γ为运行极点的编号,运行极点是指热电联产机组的产热功率和发电功率所组成的平面区域的顶点,NEi为第i台热电联产机组的运行极点编号集合,Pi γ 分别为第i台热电联产机组第γ个运行极点的有功功率和产热功率,为第i台热电联产机组在第t个调度时段的运行点对应第γ个运行极点的凸组合系数,取值范围为0~1;
热电联产机组爬坡能力约束条件为:
上式中,和分别为第i台热电联产机组的向上爬坡能力和向下爬坡能力,从机组说明书中获取,ΔT为每个调度时段的长度;
常规机组运行的有功功率约束条件为:
上式中,和分别为第i台常规机组的有功功率上限和有功功率下限,可以从机组说明书上获取;
常规机组爬坡能力的约束条件为:
上式中,和分别为第i台常规机组的向上爬坡能力和向下爬坡能力;
常规机组的旋转备用约束条件为:
上式中,和分别为第i台常规机组在第t个调度时段的向上旋转备用和向下旋转备用,其中,向上旋转备用和向下旋转备用分别指发电机组所能够提供的向上发电功率调节范围和向下发电功率调节范围;
风电机组运行约束条件为:
上式中,为第i台风电机组在第t个调度时段的可用有功功率,从风电的发电预测数据中获取,为第i台风电机组在第t个调度时段的实际有功功率;
热电联合系统中电力系统的功率平衡约束条件为:
上式中,m为电力系统中负荷的编号,上标LD代表负荷,ILD为电力系统所有负荷的编号集合,Dm,t为第m个负荷在第t个调度时段的负荷大小;
热电联合系统中电力系统的线路传输容量约束条件为:
上式中,j为电力系统中线路的编号,ILN代表电力系统线路编号集合,l为电力系统中节点的编号,IEPS代表电力系统中节点编号集合,图1是本发明方法涉及的热电联合系统结构示意图。图1所示为本发明方法涉及的热电联合系统结构示意图,以图1中的电力系统为例,IEPS={1,2,3,4,5,6},ILN={1-2,1-4,2-3,2-4,3-6,4-5,5-6}),和分别代表电力系统中与第l个节点相连接的热电联产机组、常规机组、风电机组与负荷的编号集合,Lj代表电力系统中第j个线路的最大有功功率, SFj-l为第l个节点到第j个线路的转移分布因子,转移分布因子从电力系统运行拓扑图上获取;
电力系统的旋转备用约束条件为:
上式中,SRUt和SRDt分别代表电力系统在第t个调度时段所需的向上旋转备用和向下旋转备用,和分别为第i台常规机组在第t个调度时段的向上旋转备用和向下旋转备用;
(1-2-2)辅助服务市场环境下热电联合系统中供热系统运行约束条件,包括:
(1-2-2-1)包括有热电联产机组和热锅炉的热源供热约束条件,包括:
热源产热功率与供热系统节点的供水和回水之间温差的约束条件为:
上式中,集合和分别为与供热系统节点k相连的热电联产机组编号集合和热锅炉编号集合,C为水的比热容,Mk为供热系统节点k的水流量,上标S和R分别表示供水管和回水管,下标k,t表示供热系统中第k个节点在调度时段t的相应变量,和分别表示供热系统节点k在调度时段t的供水管水温和回水管水温,供热系统的每个节点既有供水管又有回水管,集合为供热系统中与热源相连接的节点集合;
热锅炉的产热功率的约束条件,即热锅炉的产热功率需在其产热功率上限内:
上式中,表示第i个热锅炉的产热功率上限;
供热系统中与热源相连接的节点的供水管水温约束条件:
上式中,和为供热系统中任意节点k的供水管水温上限和供水管水温下限;
(1-2-2-2)换热站运行约束条件,包括:
换热功率与供热系统节点的供水管水温、回水管水温的温差为:
上式中,n为换热站的编号,为供热系统中与节点k相连接的换热站编号集合,表示第n个换热站在第t个调度时段的换热功率,表示供热系统中与换热站相连接的节点集合;
换热站的回水管水温确保在安全范围内:
上式中,和分别为供热系统中节点k的回水管水温上限和回水管水温下限;
(1-2-2-3)供热系统运行约束条件,包括:
其中,k1和k2分别为供热系统中相邻的两个节点,定义k1为流入节点,k2为流出节点,和分别为供热系统中由节点k2流向节点k1的供水流量和回水流量,和分别为供热系统中节点k1在供水管中的父节点集合和在回水管中的父节点集合,节点k1的父节点是指与节点k1直接相连的节点中水流从该节点流入至k1的节点,以图1中的供热系统为例,供水管中,节点2的父节点为节点1和节点4;回水管中,节点2的父节点为节点3和节点5,为第t个调度时段的环境温度,和分别为供热系统中由节点k2流向节点k1的供水热传导系数和回水热传导系数,和的取值由下式测算:
其中,和分别为供热系统中由节点k2流向节点k1的供水管和回水管的单位热传导系数,该系数可以从水管铭牌中获取,和分别为供热系统中由节点k2流向节点k1的供水管和回水管的长度;
和分别为供热系统的温度变量,该变量用于表示未考虑温度损失情况下由节点k2流向节点k1的供水管水温和回水管水温,温度变量的计算公式为:
上式中,分别表示供热系统中由节点k2流向节点k1时供水管和回水管所需的调度时段数,符号表示向下取整;
(2)将步骤(1)建立的辅助服务市场环境下热电联合系统调度模型转化成矩阵形式,用xE表示电力系统变量,电力系统变量包括和用xH表示供热系统变量,供热系统变量包括和则辅助服务市场环境下热电联合系统调度模型转化为如下的矩阵形式:
约束条件为:
xE∈XE
xH∈XH
DxE=ExH
上式中,CE和CH分别表示电力系统和供热系统的目标函数,其中,CE为CH为约束集合XE为电力系统的约束条件集合,即步骤(1-2-1)中的全部约束条件所构成的集合;约束集合XH表示供热系统约束条件集合,即步骤(1-2-2)中的全部约束条件所构成的集合;约束DxE=ExH表示电力系统和供热系统耦合约束条件,即 变量xE中的和变量xH中的对应相等,D、E的每一行与每一个电力系统和供热系统耦合约束条件一一对应, D的每一列与电力系统中的每一个变量一一对应,E的每一列与供热系统中的每一个变量一一对应,其中D、E的每一个元素为该元素所在行所对应的约束条件中其列所代表的变量的系数;
(3)采用变方向乘子算法,对上述步骤(2)得到的矩阵形式的辅助服务市场环境下热电联合系统调度模型进行求解,步骤如下:
(3-1)初始化:将迭代次数u初始化为0,初始化λ(u)=0,定义α=1,并定义如下式的增广拉格朗日函数A(xE,xH,λ):
(3-2)按照下式求解下面的电力系统的调度策略:
(3-3)按照下式求解下面的供热系统的调度策略:
(3-4)按照下式更新λ(u+1):
(3-5)根据上述步骤(3-2)和步骤(3-3)的计算结果进行判断,若||DxE-ExH||<10-6,则进行步骤(3-6),若||DxE-ExH||≥10-6,则将迭代次数u增加1,返回步骤(3-2)继续计算;
(3-6)将得到的最优解作为辅助服务市场环境下热电联合系统调度的调度参数。
Claims (1)
1.一种辅助服务市场环境下的热电联合系统分布式调度方法,其特征在于该方法包括以下步骤:
(1)建立辅助服务市场环境下热电联合系统调度模型,调度模型由目标函数和约束条件构成,具体包括:
(1-1)辅助服务市场环境下热电联合系统调度模型的目标函数:
热电联合系统调度模型的目标函数为电力系统运行成本、电力系统备用辅助服务成本与供热系统运行成本之和的最小化,表达式为:
上式中,T为调度时段集合,上标CHP、TU、WD、HB和AS分别表示热电联产机组、常规机组、风电机组、热锅炉和辅助服务,ICHP、ITU、IWD和IHB分别为热电联合系统中的热电联产机组、常规机组、风电机组和热锅炉的编号集合,下标i,t表示第i台机组在调度时段t的相应变量,和分别为热电联合系统中热电联产机组、常规机组、风电机组和热锅炉的生产成本函数,为电力系统中常规机组提供的备用辅助服务成本函数,i为热电联产机组、常规机组、风电机组和热锅炉中机组或锅炉的编号,t为调度时段;
其中的热电联产机组的生产成本函数为:
上式中,和分别为生产成本系数,生产成本系数为热电联产机组固有参数,和分别为第i台热电联产机组在第t个调度时段的有功功率和产热功率;
常规机组的生产成本函数为:
上式中,和为常规机组的发电成本系数,该系数为常规机组的固有参数,为第i台常规机组在第t个调度时段的有功功率;
风电机组的生产成本函数为:
上式中,为弃风惩罚因子,取值范围为1~10,弃风惩罚因子的取值根据对风电的消纳需求确定,由电力系统调度中心根据调度结果反馈进行调节,为第i台风电机组在第t个调度时段的可用有功功率,从风电的发电预测数据中获取,为第i台风电机组在第t个调度时段的实际有功功率;
常规机组提供的备用辅助服务成本函数为:
上式中,上标RU和RD分别表示向上旋转备用和向下旋转备用,和分别为电力系统在第t个调度时段的向上和向下的备用辅助服务竞价价格,该价格由辅助服务市场的参与主体中报价最高者处获取,和分别为第i台常规机组在第t个调度时段的向上旋转备用和向下旋转备用;
热锅炉的生产成本函数为:
上式中,为热锅炉的产热成本系数,为热锅炉固有参数,为第i台热锅炉在第t个调度时段的产热功率;
(1-2)辅助服务市场环境下热电联合系统调度模型的约束条件,包括:
(1-2-1)辅助服务市场环境下热电联合系统中电力系统的运行约束条件:
热电联产机组运行的有功功率和产热功率的约束条件为:
上式中,γ为运行极点的编号,运行极点是指热电联产机组的产热功率和发电功率所组成的平面区域的顶点,NEi为第i台热电联产机组的运行极点编号集合,分别为第i台热电联产机组第γ个运行极点的有功功率和产热功率,为第i台热电联产机组在第t个调度时段的运行点对应第γ个运行极点的凸组合系数,取值范围为0~1;
热电联产机组爬坡能力约束条件为:
上式中,和分别为第i台热电联产机组的向上爬坡能力和向下爬坡能力,从机组说明书中获取,ΔT为每个调度时段的长度;
常规机组运行的有功功率约束条件为:
上式中,和Pi TU 分别为第i台常规机组的有功功率上限和有功功率下限;
常规机组爬坡能力的约束条件为:
上式中,和分别为第i台常规机组的向上爬坡能力和向下爬坡能力;
常规机组的旋转备用约束条件为:
上式中,和分别为第i台常规机组在第t个调度时段的向上旋转备用和向下旋转备用;
风电机组运行约束条件为:
上式中,为第i台风电机组在第t个调度时段的可用有功功率,从风电的发电预测数据中获取,为第i台风电机组在第t个调度时段的实际有功功率;
热电联合系统中电力系统的功率平衡约束条件为:
上式中,m为电力系统中负荷的编号,上标LD代表负荷,ILD为电力系统所有负荷的编号集合,Dm,t为第m个负荷在第t个调度时段的负荷大小;
热电联合系统中电力系统的线路传输容量约束条件为:
上式中,j为电力系统中线路的编号,ILN代表电力系统线路编号集合,l为电力系统中节点的编号,IEPS代表电力系统中节点编号集合,和分别代表电力系统中与第l个节点相连接的热电联产机组、常规机组、风电机组与负荷的编号集合,Lj代表电力系统中第j个线路的最大有功功率,SFj-l为第l个节点到第j个线路的转移分布因子,转移分布因子从电力系统运行拓扑图上获取;
电力系统的旋转备用约束条件为:
上式中,SRUt和SRDt分别代表电力系统在第t个调度时段所需的向上旋转备用和向下旋转备用,和分别为第i台常规机组在第t个调度时段的向上旋转备用和向下旋转备用;
(1-2-2)辅助服务市场环境下热电联合系统中供热系统运行约束条件,包括:
(1-2-2-1)包括有热电联产机组和热锅炉的热源供热约束条件,包括:
热源产热功率与供热系统节点的供水和回水之间温差的约束条件为:
上式中,集合和分别为与供热系统节点k相连的热电联产机组编号集合和热锅炉编号集合,C为水的比热容,Mk为供热系统节点k的水流量,上标S和R分别表示供水管和回水管,下标k,t表示供热系统中第k个节点在调度时段t的相应变量,和分别表示供热系统节点k在调度时段t的供水管水温和回水管水温,集合为供热系统中与热源相连接的节点集合;
热锅炉的产热功率的约束条件,即热锅炉的产热功率需在其产热功率上限内:
上式中,表示第i个热锅炉的产热功率上限;
供热系统中与热源相连接的节点的供水管水温约束条件:
上式中,和为供热系统中任意节点k的供水管水温上限和供水管水温下限;
(1-2-2-2)换热站运行约束条件,包括:
换热功率与供热系统节点的供水管水温、回水管水温的温差为:
上式中,n为换热站的编号,为供热系统中与节点k相连接的换热站编号集合,表示第n个换热站在第t个调度时段的换热功率,表示供热系统中与换热站相连接的节点集合;
换热站的回水管水温确保在安全范围内:
上式中,和分别为供热系统中节点k的回水管水温上限和回水管水温下限;
(1-2-2-3)供热系统运行约束条件,包括:
其中,k1和k2分别为供热系统中相邻的两个节点,定义k1为流入节点,k2为流出节点,和分别为供热系统中由节点k2流向节点k1的供水流量和回水流量,和分别为供热系统中节点k1在供水管中的父节点集合和在回水管中的父节点集合,为第t个调度时段的环境温度,和分别为供热系统中由节点k2流向节点k1的供水热传导系数和回水热传导系数,和的取值由下式测算:
其中,和分别为供热系统中由节点k2流向节点k1的供水管和回水管的单位热传导系数,和分别为供热系统中由节点k2流向节点k1的供水管和回水管的长度;
和分别为供热系统的温度变量,温度变量的计算公式为:
上式中,分别表示供热系统中由节点k2流向节点k1时供水管和回水管所需的调度时段数,符号表示向下取整;
(2)将步骤(1)建立的辅助服务市场环境下热电联合系统调度模型转化成矩阵形式,用xE表示电力系统变量,电力系统变量包括和用xH表示供热系统变量,供热系统变量包括和则辅助服务市场环境下热电联合系统调度模型转化为如下的矩阵形式:
约束条件为:
xE∈XE
xH∈XH
DxE=ExH
上式中,CE和CH分别表示电力系统和供热系统的目标函数,其中,CE为CH为约束集合XE为电力系统的约束条件集合,即步骤(1-2-1)中的全部约束条件所构成的集合;约束集合XH表示供热系统约束条件集合,即步骤(1-2-2)中的全部约束条件所构成的集合;约束DxE=ExH表示电力系统和供热系统耦合约束条件,即变量xE中的和变量xH中的对应相等,D、E的每一行与每一个电力系统和供热系统耦合约束条件一一对应,D的每一列与电力系统中的每一个变量一一对应,E的每一列与供热系统中的每一个变量一一对应,其中D、E的每一个元素为该元素所在行所对应的约束条件中其列所代表的变量的系数;
(3)采用变方向乘子算法,对上述步骤(2)得到的矩阵形式的辅助服务市场环境下热电联合系统调度模型进行求解,步骤如下:
(3-1)初始化:将迭代次数u初始化为0,初始化λ(u)=0,定义α=1,并定义如下式的增广拉格朗日函数A(xE,xH,λ):
(3-2)按照下式求解下面的电力系统的调度策略:
(3-3)按照下式求解下面的供热系统的调度策略:
(3-4)按照下式更新λ(u+1):
(3-5)根据上述步骤(3-2)和步骤(3-3)的计算结果进行判断,若||DxE-ExH||<10-6,则进行步骤(3-6),若||DxE-ExH||≥10-6,则将迭代次数u增加1,返回步骤(3-2)继续计算;
(3-6)将得到的最优解作为辅助服务市场环境下热电联合系统调度的调度参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710274062.5A CN106992555B (zh) | 2017-04-25 | 2017-04-25 | 一种辅助服务市场环境下的热电联合系统分布式调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710274062.5A CN106992555B (zh) | 2017-04-25 | 2017-04-25 | 一种辅助服务市场环境下的热电联合系统分布式调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106992555A CN106992555A (zh) | 2017-07-28 |
CN106992555B true CN106992555B (zh) | 2019-06-04 |
Family
ID=59418130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710274062.5A Expired - Fee Related CN106992555B (zh) | 2017-04-25 | 2017-04-25 | 一种辅助服务市场环境下的热电联合系统分布式调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106992555B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106815661B (zh) | 2017-02-22 | 2020-10-20 | 清华大学 | 一种热电联合系统的分解协调调度方法 |
CN108155644B (zh) * | 2017-12-26 | 2020-04-24 | 清华大学 | 计及辅助服务的动态供热系统与电力系统耦合调度方法 |
CN108667078B (zh) * | 2018-04-11 | 2019-12-10 | 东南大学 | 一种区外电源引致辅助服务成本计算方法 |
CN110991845B (zh) * | 2019-11-25 | 2023-06-23 | 国网节能服务有限公司 | 一种电-热耦合系统分布式协同调度方法 |
CN111146819B (zh) * | 2019-12-26 | 2021-03-02 | 上海电力大学 | 一种考虑热网特性的电热联合系统辅助服务调度方法 |
CN113324283B (zh) * | 2021-05-12 | 2022-04-22 | 华南理工大学 | 一种计及温度半动态特性的集中供热网络等值方法 |
CN117634682B (zh) * | 2023-11-28 | 2024-06-18 | 国网吉林省电力有限公司经济技术研究院 | 一种寒冷地区电-热联供型虚拟电厂优化调控方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102506477A (zh) * | 2011-10-23 | 2012-06-20 | 西安交通大学 | 热电联产机组与风力发电联合制冷系统及其调度方法 |
CN106447122A (zh) * | 2016-10-12 | 2017-02-22 | 国网上海市电力公司 | 一种区域型能源互联网及其一体化优化规划方法 |
-
2017
- 2017-04-25 CN CN201710274062.5A patent/CN106992555B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102506477A (zh) * | 2011-10-23 | 2012-06-20 | 西安交通大学 | 热电联产机组与风力发电联合制冷系统及其调度方法 |
CN106447122A (zh) * | 2016-10-12 | 2017-02-22 | 国网上海市电力公司 | 一种区域型能源互联网及其一体化优化规划方法 |
Non-Patent Citations (2)
Title |
---|
《考虑热负荷二维可控性的微网热电协调调度》;邹云阳等;《电力系统自动化》;20170325;第41卷(第6期);第13-19页 |
Wind-Thermal Economic and Environmental Scheduling Incorporating Demand Response;E. Heydarian-Forushani等;《2013 Smart Grid Conference (SGC)》;20131218;第36-40页 |
Also Published As
Publication number | Publication date |
---|---|
CN106992555A (zh) | 2017-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106992555B (zh) | 一种辅助服务市场环境下的热电联合系统分布式调度方法 | |
CN106815661B (zh) | 一种热电联合系统的分解协调调度方法 | |
CN106339772B (zh) | 基于供热管网储热效益的热-电联合优化调度方法 | |
CN106786753B (zh) | 多用户的区域能源互联网的系统及其调节方法 | |
CN109447323A (zh) | 一种计及节点热价的综合能源系统两阶段容量配置方法 | |
CN108009683A (zh) | 一种大容量直控负荷提高风电接纳能力的调度方法 | |
CN109636037A (zh) | 一种基于pso算法的多热源环状热网水力优化调度方法 | |
CN110535128A (zh) | 基于用能舒适度的多区域综合能源系统协同调度方法 | |
CN109193628B (zh) | 一种基于一致性的综合能源系统能量管理方法 | |
CN113469430B (zh) | 综合用能园区的多能互补容量配置方法 | |
CN109389248B (zh) | 基于区域热网的综合能源系统热电功率协调调度方法 | |
CN114330827B (zh) | 多能流虚拟电厂分布式鲁棒自调度优化方法及其应用 | |
CN116341881B (zh) | 计及热网灵活性的电-热系统鲁棒超前调度方法及系统 | |
CN110688744A (zh) | 一种应用于热电耦合网络的异步分布式状态估计方法 | |
CN111062598B (zh) | 综合能源系统的分布式优化调度方法及系统 | |
CN114322044B (zh) | 一种综合能源系统及其运行控制方法 | |
CN117035202B (zh) | 考虑需求响应的电热综合能源系统双层协同扩展规划方法 | |
CN112214883B (zh) | 一种微能源网综合能源系统运行优化方法和系统 | |
CN108594658A (zh) | 一种电-气耦合系统概率最大负荷裕度多目标优化模型及其求解方法 | |
CN105279709A (zh) | 一种基于热网热惯性的电网日前优化调度方法 | |
CN114048908A (zh) | 统一电网-分布式热网系统多时间尺度加强区间优化方法 | |
CN108155644B (zh) | 计及辅助服务的动态供热系统与电力系统耦合调度方法 | |
Yuan et al. | An advanced multicarrier residential energy hub system based on mixed integer linear programming | |
CN111724026B (zh) | 一种多能源网络与配水网耦合运行的优化方法 | |
CN104766131A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190604 |
|
CF01 | Termination of patent right due to non-payment of annual fee |