本发明涉及一种辅助服务市场环境下的热电联合系统分布式调度方法,属于电力系统运行控制技术领域。
背景技术:
在我国三北地区风电弃电现象严重,在冬天受到城区供热系统的最小运行约束限制,热电联产机组无法减小其出力为风电的消纳提供空间。具体来说,冬季夜晚属于风电高发期,但夜晚同时具有很高的居民热负荷,导致了风电需要被限制来为热电联产机组提供出力空间,因此供热和风电消纳需求之间存在矛盾。由于供热管网具有良好的储热效益,综合考虑电力系统和供热系统进行调度可以实现负荷调峰,有效提高系统运行的灵活性,减小弃风现象。
在辅助服务市场环境下,常规机组的有功备用需要被考虑到备用成本中,这导致了在辅助服务市场环境下电力系统的调度与传统调度存在不同,需要针对辅助服务市场独立考虑。
电力系统和供热系统分别由电力公司与供热公司独立运行,由于调度独立性,难以实现对电力系统和供热系统完整模型进行统一调度。因此,提出辅助服务市场环境下热电联合系统分布式调度算法方法十分必要。具体来说,电力系统和供热系统可以分别对各自管辖的区域进行调度优化,通过边界条件的迭代获得辅助服务市场环境下热电联合系统调度的全局最优解。
变方向乘子算法是一种分布式求解算法,可将复杂的数学优化问题转化为若干相对简单的问题分布式进行求解,通过迭代可获得全局最优解。
技术实现要素:
本发明的目的是提出一种辅助服务市场环境下的热电联合系统分布式调度方法,针对电力系统与供热系统分别建模,考虑电力系统和供热系统的耦合约束,建立辅助服务市场环境下的热电联合调度模型。针对辅助服务市场环境下的热电联合调度模型,基于变方向乘子算法,提出辅助服务市场环境下的热电联合调度模型的分布式求解算法。
本发明提出的辅助服务市场环境下的热电联合系统分布式调度方法,包括以下步骤:
(1)建立辅助服务市场环境下热电联合系统调度模型,调度模型由目标函数和约束条件构成,具体包括:
(1-1)辅助服务市场环境下热电联合系统调度模型的目标函数:
热电联合系统调度模型的目标函数为电力系统运行成本、电力系统备用辅助服务成本与供热系统运行成本之和的最小化,表达式为:
上式中,t为调度时段集合,上标chp、tu、wd、hb和as分别表示热电联产机组、常规机组、风电机组、热锅炉和辅助服务,ichp、itu、iwd和ihb分别为热电联合系统中的热电联产机组、常规机组、风电机组和热锅炉的编号集合,下标i,t表示第i台机组在调度时段t的相应变量,
其中的热电联产机组的生产成本函数
上式中,
常规机组的生产成本函数
上式中,
风电机组的生产成本函数
上式中,
常规机组提供的备用辅助服务成本函数
上式中,上标ru和rd分别表示向上和向下旋转备用,
热锅炉的生产成本函数
上式中,
(1-2)辅助服务市场环境下热电联合系统调度模型的约束条件,包括:
(1-2-1)辅助服务市场环境下热电联合系统中电力系统的运行约束条件:
热电联产机组运行的有功功率和产热功率的约束条件为:
上式中,γ为运行极点的编号,运行极点是指热电联产机组的产热功率和发电功率所组成的平面区域的顶点,nei为第i台热电联产机组的运行极点编号集合,piγ,
热电联产机组爬坡能力约束条件为:
上式中,
常规机组运行的有功功率约束条件为:
上式中,
常规机组爬坡能力的约束条件为:
上式中,
常规机组的旋转备用约束条件为:
上式中,
风电机组运行约束条件为:
上式中,
热电联合系统中电力系统的功率平衡约束条件为:
上式中,m为电力系统中负荷的编号,上标ld代表负荷,ild为电力系统所有负荷的编号集合,dm,t为第m个负荷在第t个调度时段的负荷大小;
热电联合系统中电力系统的线路传输容量约束条件为:
上式中,j为电力系统中线路的编号,iln代表电力系统线路编号集合,l为电力系统中节点的编号,ieps代表电力系统中节点编号集合,
电力系统的旋转备用约束条件为:
上式中,srut和srdt分别代表电力系统在第t个调度时段所需的向上旋转备用和向下旋转备用,
(1-2-2)辅助服务市场环境下热电联合系统中供热系统运行约束条件,包括:
(1-2-2-1)包括有热电联产机组和热锅炉的热源供热约束条件,包括:
热源产热功率与供热系统节点的供水和回水之间温差的约束条件为:
上式中,集合
热锅炉的产热功率
上式中,
供热系统中与热源相连接的节点的供水管水温约束条件:
上式中,
(1-2-2-2)换热站运行约束条件,包括:
换热功率与供热系统节点的供水管水温、回水管水温的温差为:
上式中,n为换热站的编号,
换热站的回水管水温确保在安全范围内:
上式中,
(1-2-2-3)供热系统运行约束条件,包括:
其中,k1和k2分别为供热系统中相邻的两个节点,定义k1为流入节点,k2为流出节点,
其中,
上式中,
(2)将步骤(1)建立的辅助服务市场环境下热电联合系统调度模型转化成矩阵形式,用xe表示电力系统变量,电力系统变量包括
约束条件为:
xe∈xe
xh∈xh
dxe=exh
上式中,ce和ch分别表示电力系统和供热系统的目标函数,其中,ce为
(3)采用变方向乘子算法,对上述步骤(2)得到的矩阵形式的辅助服务市场环境下热电联合系统调度模型进行求解,步骤如下:
(3-1)初始化:将迭代次数u初始化为0,初始化
(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的相应变量,
其中的热电联产机组的生产成本函数
上式中,
常规机组的生产成本函数
上式中,
风电机组的生产成本函数
上式中,
常规机组提供的备用辅助服务成本函数
上式中,上标ru和rd分别表示向上和向下旋转备用,
热锅炉的生产成本函数
上式中,
(1-2)辅助服务市场环境下热电联合系统调度模型的约束条件,包括:
(1-2-1)辅助服务市场环境下热电联合系统中电力系统的运行约束条件:
热电联产机组运行的有功功率和产热功率的约束条件为:
上式中,γ为运行极点的编号,运行极点是指热电联产机组的产热功率和发电功率所组成的平面区域的顶点,nei为第i台热电联产机组的运行极点编号集合,piγ
热电联产机组爬坡能力约束条件为:
上式中,
常规机组运行的有功功率约束条件为:
上式中,
常规机组爬坡能力的约束条件为:
上式中,
常规机组的旋转备用约束条件为:
上式中,
风电机组运行约束条件为:
上式中,
热电联合系统中电力系统的功率平衡约束条件为:
上式中,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}),
电力系统的旋转备用约束条件为:
上式中,srut和srdt分别代表电力系统在第t个调度时段所需的向上旋转备用和向下旋转备用,
(1-2-2)辅助服务市场环境下热电联合系统中供热系统运行约束条件,包括:
(1-2-2-1)包括有热电联产机组和热锅炉的热源供热约束条件,包括:
热源产热功率与供热系统节点的供水和回水之间温差的约束条件为:
上式中,集合
热锅炉的产热功率
上式中,
供热系统中与热源相连接的节点的供水管水温约束条件:
上式中,
(1-2-2-2)换热站运行约束条件,包括:
换热功率与供热系统节点的供水管水温、回水管水温的温差为:
上式中,n为换热站的编号,
换热站的回水管水温确保在安全范围内:
上式中,
(1-2-2-3)供热系统运行约束条件,包括:
其中,k1和k2分别为供热系统中相邻的两个节点,定义k1为流入节点,k2为流出节点,
其中,
上式中,
(2)将步骤(1)建立的辅助服务市场环境下热电联合系统调度模型转化成矩阵形式,用xe表示电力系统变量,电力系统变量包括
约束条件为:
xe∈xe
xh∈xh
dxe=exh
上式中,ce和ch分别表示电力系统和供热系统的目标函数,其中,ce为
(3)采用变方向乘子算法,对上述步骤(2)得到的矩阵形式的辅助服务市场环境下热电联合系统调度模型进行求解,步骤如下:
(3-1)初始化:将迭代次数u初始化为0,初始化
(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)将得到的最优解作为辅助服务市场环境下热电联合系统调度的调度参数。