发明内容
针对上述存在的问题,如何使价格型需求响应参与综合能源***运行优化是实现需求侧和供给侧之间良好互动机制的基础,同时如何利用精细化的调控策略对需求响应资源进行建模分析,是需要解决的主要技术问题。
利用价格型需求响应模型进行综合能源***运行优化的基本思想是:优化前,实现能源结构的确定、设备的性能参数及功率约束参数的确定、各类能源价格的导入、用户响应负荷的计算以及引导用户实施需求响应的实时电价的确定;优化中,按照设定的优化场景设定优化目标,利用NSGA-II方法和模糊隶属度方法进行寻优,最终确定最佳的目标结果;优化后,实施优化目标结果以及负荷削减情况的对比分析,研究***负荷和运行成本在需求响应前后的变化。
为了解决所述技术问题,本发明采用的技术方案是:一种基于价格型需求响应模型的综合能源运行优化方法,包括以下步骤:
S1:获取基础数据及模型参数
获取综合能源***运行的基础数据以及模型对应的参数,包括各类设备的性能参数及功率约束参数、各类能源价格信息、用户的能源负荷需求、用户实施需求响应的实时电价等数据。
S2:构建基于实时电价的需求响应模型
以分时电价为基础建立实时电价需求响应模型,实时电价计算公式如下:
式中,Pspot(t)表示园区综合能源***对外的实时售电价格;Q(t)表示用户t时刻的电负荷;N表示仿真的时段数,仿真时间取24小时,仿真间隔为1小时;Ptou(t)表示基准电价,取国网公布的分时电价,即园区综合能源***购买电力的分时电价;Qmin和Qmax表示对外售电的价格的约束下限和上限。
根据用户需求价格弹性原理,价格型需求响应前后负荷量之间的关系式,如下所示:
式中,Q0(t)表示响应前t时刻的用户负荷需求;σ(t,t0)表示用户在(t,t0)时刻的需求价格弹性;P0(t0)表示响应前t0时刻的电价,P(t0)表示响应后t0时刻的电价。
S3:构建模型目标函数与约束条件
园区综合能源***经济运行的最大目标是在既定的投资水平下既满足负荷需求又最大限度降低能源成本,同时实现碳排放量降低的目标。因此本发明采取双目标优化的方式对***的运行成本和碳排放量两个目标进行优化,***根据负荷情况产生实时电价,用户根据该价格参与需求响应。
S31):***总运行成本的计算
综合能源***的运行成本主要包含燃气轮机消耗天然气的成本、燃气锅炉消耗天然气的成本、向电网购电成本、各设备单元的运行维护成本。以一个调度周期,以综合需求响应以及合理的安排各单位出力,使园区综合能源***总运行成本最小,收益最大。计算公式如下:
Min F1=Fope-Fin,
式中,Fope表示运行成本;Fin表示运行收入;Pgrid(t)是综合能源***和电网之间的交换功率,kW;Cgrid(t)为交换功率的单位价格,元/kW;△t为调度周期,小时;Cgas(t)为购买天然气的单位价格,元/m3;Qgb(t)为燃气锅炉的产热功率,kW;ηgb为燃气锅炉的产热效率,%;Pgt(t)为燃气轮机的发电功率,kW;a、b和c为燃气轮机的燃料成本系数;HVgas为天然气的低热值,kWh/m3;Pbess,c/bess.f(t)为储能电池的充放电功率,kW;Capbess为储能电池的总容量,kW;Ccir为储能电池单次完全充放电的成本;Mcircle为储能电池的充放电次数;Xbess表示储能电池的数量;其中假设储能电池每次充放电的使用成本保持相同。Mj equ表示***中第j个设备的单位功率维护成本,元/kW;Pj equ(t)表示***中第j个设备的出力功率,kW;Qcool(t)表示***内冷负荷;Pcool(t)表示满足用户冷负荷的实时价格,Qhot(t)表示***内热负荷,Phot(t)表示满足用户热负荷的实时价格,Pele(t)表示***内的电负荷,Pspot(t)表示园区综合能源***对外的实时售电价格。
S32):***碳排放量的计算
有效的降低碳排放量将不仅对环境保护发挥重大作用,而且也有利于园区的长远运行,计算公式如下:
式中,Ccarb是碳排放量,kg;α是发电过程中污染物排放系数;Pgrid(t)是综合能源***和电网之间的交换功率;Pin,NG(t)表示天然气输入功率;β是天然气燃烧中污染物排放系数,HVgas为天然气的低热值。
S33):约束条件的确定
在园区综合能源***中,考虑到***和外部的交互关系以及***与能源网的安全,包括***内部设备自身以及***与外部网络之间能量交换的功率必须遵循约束条件。约束条件分为能源平衡约束、能源供应约束以及能源***运行约束三种。
能源平衡约束包含电功率平衡约束、热功率平衡约束和冷功率平衡约束。
电功率平衡约束:
Pgrid(t)+Pgt(t)+Ppv(t)+Pbess,d(t)=Pele(t)+Per(t)+Pbess,c(t),
式中,Pgrid(t)表示综合能源***和电网之间t时刻的交互功率;Pgt(t)表示燃气轮机t时刻的出力;Ppv(t)表示光伏在t时刻的出力;Pbess,d(t)表示储能电池在t时刻的放电功率;Pele(t)、Per(t)和Pbess,c(t)分别表示***内t时刻的电负荷、电制冷机的电负荷以及储能电池的充电功率。
热功率平衡约束:
式中,ηgt表示燃气轮机的产热效率;Pgt(t)表示燃气轮机在t时刻的发电功率;ηwhb表示余热锅炉的产热效率;ω表示热分配系数,用以分配供给制冷机和制热水的热量;Qgb(t)表示燃机锅炉在t时刻的出力;Qhot(t)表示***内t时刻的热负荷;θ表示热负荷的调节系数。
冷功率平衡约束:
ηlbuQlbu(t)+ηerQer(t)=Qcool(t),
式中,ηlbu表示溴化锂制冷机的制冷效率;Qlbu(t)表示溴化锂制冷机在t时刻的制冷功率;ηer表示电制冷机的制冷效率;Qer(t)表示电制冷机在t时刻的制冷功率;Qcool(t)表示***内t时刻的冷负荷。
为保证***与能源网的安全,综合能源***需遵循能源供应约束:
Pe_min≤|Pe_grid|≤Pe_max
Pg_min≤|Pg_network|≤Pg_max,
式中,Pe_min(t)表示电网与综合能源***之间的最小电量交换功率;Pe_max(t)表示电网与综合能源***之间的最大电量交换频率;Pg_max(t)表示天然气网向***供应天然气的最大功率;Pg_min(t)表示天然气网向***供应天然气的最小功率。能源***运行约束包含设备运行约束、储能电池约束、天然气连接节点压力约束和分配系数约束。
设备运行约束:
式中,Pgrid,max(t)为综合能源***向电网购电的最大功率;Pi,min(t)为第i种产电设备在t时刻的电功率下限;Pi(t)表示***中第i种产电设备在t时刻的产电功率;Qi,min(t)为第i种产热设备在t时刻的热功率下限;Qi(t)表示***中第i种产热设备在t时刻的产热功率;Ei(t)表示***中第i种产电设备在t时刻的出力;Smin为***中设备向上爬坡速率。
储能电池约束:
式中,soct和soct+1分别为储能电池在t和t+1时刻的荷电状态;δbess为储能电池的自放电率;ηbess,c为储能电池的充电效率;ηbess,d为储能电池的放电效率;αbess,c和αbess,d分别为储能电池的充电状态和放电状态,取值为0-1变量,αbess,c取1表示充电,αbess,d取1表示放电,不可同时充放电;pbess,c和pbess,d分别表示储能电池的充电和放电功率;socmin和socmax分别表示储能电池荷电状态的上下限;socL和socT分别表示储能电池的起始和终止荷电状态;pbess,c max和pbess,d max分别表示储能电池的最大充放电功率;Imax ch和Imax dis分别表示电池允许的最大充放电电流;Pinv表示***中逆变器的容量;Vbat表示电池的额定电压;Qmax表示电池容量。
天然气连接节点压力约束:
式中,ρmin j和ρmax j表示气网允许的最小压力和最大压力。
分配系数约束:
式中,vc为某种形式能源分配到不同能源转换装置的分配系数;c为能源转换装置的编号。
S4:求解综合能源***运行成本与碳排放量最低运行计划
综合能源***的运行需要从经济成本、环境成本等多方面的因素综合考虑,而为寻找最佳规划解而设置的***模拟运行过程又是典型的多目标约束下的优化问题,涉及需要线性、非线性、随机性的因素。相比多约束、单目标优化模型获得唯一解,本发明建立的多指标约束、多目标协同优化模型获得由多组解组成的帕累托解。综合考虑每个能源供应单元的输出功率和园区综合能源***中的***能量需求以及模型的特征,本发明采用带精英策略的非支配排序的遗传算法(NSGA-Ⅱ)求解***运行优化问题,求解得出的帕累托解采用模糊隶属度函数进行计算评价。
S41):优化问题描述及变量识别
实现经济成本最低、环境成本可承受是一个园区健康持续发展的必然要求,在此基础上通过建立多指标约束、多目标协同优化模型是满足该要求的关键步骤。多目标优化问题可综合概括为多目标混合整数非线性规划的求解,其一般形式如下:
式中,x、y为优化变量;f(x,y)为优化目标函数;X为变量集;Y为所有可行解的集合;J(x,y)为不等式约束;H(x,y)为等式约束。
在综合能源***运行优化中,通常优化所得的是***内各个设备的运行状态和每个小时的设备出力值。本发明采用二进制的方式来表示设备的运行状态,设X_s=1表示设备处于运行状态,X_s=0表示设备处于待机状态,采用X_p表示设备的出力水平,当设备的自动调解功率输出范围没有限制时,设备的实际出力为Preal=X_p;当设备的功率调节范围具有上下限约束时,设备的实际出力为Preal=X_p×X_s。综合本发明设定的优化目标函数、优化变量的类型等情况,最后模型的优化变量可识别如下:
式中,X_p表示***内各个设备的输出功率集合;X_s表示***内各个设备的运行状态;Sgt(t)表示燃气轮机运行状态;Spv(t)表示光伏发电设备运行状态;Swhb(t)表示余热锅炉运行状态;Slbs(t)表示溴化锂制冷机运行状态;Sgfb(t)表示燃气锅炉运行状态;Ser(t)表示电制冷机运行状态;Sbees(t)表示储能电池运行状态;Pgt(t)表示燃气轮机输出功率;Ppv(t)表示光伏发电输出功率;Pwhb(t)表示余热锅炉输出功率;Plbs(t)表示溴化锂制冷机输出功率;Pgfb(t)表示燃气锅炉输出功率;Per(t)表示电制冷机输出功率;Pbees(t)表示储能电池输出功率。
综合能源***的运行优化问题涉及较多因素,建立的优化模型中变量较多且目标函数复杂,包含大量的非线性方程和变量,常规数学模型求解难度较大。在优化模型中存在多个变量、多种约束、目标函数较为复杂,传统的数学方法不适合求解多目标问题。综合考虑每个能源供应单元输出功率和园区综合能源***中的***能量需求以及模型的特征,采用带精英策略的非支配排序的遗传算法(NSGA-II)求解***运行优化问题。
S42):采用NSGA-II算法对模型求解
本发明采用NSGA-II算法对模型进行求解,通过计算获得优化后的由多组解组成的非劣解集。首先随机产生规模为N的初始种群,非支配排序后通过遗传算法的选择、交叉、变异三个基本操作得到第一代子代种群;从第二代开始将父代种群与子代种群合并,进行快速非支配排序,同时对每个非支配层中的个体进行拥挤度计算,根据非支配关系以及个体的拥挤度选取合适的个体组成新的父代种群。以此类推,当满足设定的程序结束条件时输出帕累托解集。
S43):选择最优解
模糊隶属度理论通过合成每个目标函数的模糊隶属度来从帕累托解集中选择最满意解。首先通过公式计算帕累托集中目标函数的隶属度;根据决策者的偏好设置权重值,并计算多目标函数的最优解的隶属权重值,最大值对应的帕累托解集是最优解。本发明中,权重设置为相等权重,隶属度值计算公式为:
式中,Fi表示目标函数的第i个目标解;Fmin i表示目标函数的最小解;Fmax i表示目标函数的最大解。N表示设定的种群数目;Di表示设定的目标函数的权重;Nobj表示设定的目标函数的数量;uk表示目标函数的隶属度值。
S5:算例优化结果分析
将得到的基础数据与边界参数输入价格型需求响应模型,按“自发自用,余电上网”的模式,对指定算例进行仿真分析,仿真周期为一天,仿真步长为一小时。通过优化目标结果的对比分析负荷削减与成本节约情况,评估模型对综合能源***的优化程度。
本发明的有益效果:
1.已有的综合能源***需求侧响应评估多集中考虑电负荷响应的单一评估,本发明考虑了综合能源***的电、热、冷三种负荷同时响应的效果评估。
2.已有的综合能源***需求侧响应优化方法多对***运行成本单一目标进行优化,本发明采取双目标优化的方式对***的运行成本和碳排放量两个目标进行优化。
3.NSGA-II是一种多目标遗传算法,它在第一代NSGA基础上提出了快速非支配排序算法,并引进了精英策略和拥挤距离的概念,具有运行速度快、优化结果精度高、解集收敛性好的优点。
4.本发明结合国内综合能源***需求侧响应技术发展现状,以园区综合能源***为研究对象,统筹考虑电力市场价格数据、园区设备特性数据、园区负荷数据等因素,建立综合能源***价格型需求响应模型,首次采用NSGA-II和模糊隶属度函数求出基于价格型需求响应模型的综合能源***运行优化最优解。相关计算过程原理简单,实用性强,可有效降低运行成本以及***负荷,促进用户侧参与综合能源***运行。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的说明。
本实施例以某工业园区为例,公开一种基于价格型需求响应模型的园区综合能源***运行优化方法,具体描述如下:
S1)获取基础数据及模型参数
获取综合能源***运行的基础数据以及模型对应的参数,包括各类设备的性能参数及功率约束参数、各类能源价格信息、用户的能源负荷需求、用户实施需求响应的实时电价等数据,具体:综合能源***和电网之间交换功率Pgrid、交换功率单位价格Cgrid、调度周期△t、天然气单位价格Cgas、燃气锅炉产热功率Qgb、燃气锅炉产热效率ηgb,、燃气轮机发电功率Pgt、燃气轮机燃料成本系数a、b和c、天然气低热值HVgas、储能电池充放电功率Pbess,c/bess.f、储能电池总容量Capbess、储能电池单次完全充放电成本Ccir、储能电池充放电次数Mcircle、储能电池数量Xbess、***中第j个设备的单位功率维护成本Mj equ、***中第j个设备出力功率Pj equ、碳排放量Ccarb、发电过程中污染物排放系数α、电网输入功率PE、天然气输入功率Pin,NG、天然气燃烧中污染物排放系数β、光伏出力Ppv、储能电池放电功率Pbess,d、***内电负荷Pele、电制冷机的电负荷Per、储能电池充电功率Pbess,c、燃气轮机产热效率ηgt、余热锅炉产热效率ηwhb、热分配系数ω、***内热负荷Qhot、热负荷调节系数θ、溴化锂制冷机制冷效率ηlbu、溴化锂制冷机在t时刻制冷功率Qlbu、电制冷机制冷效率ηer、电制冷机制冷功率Qer、、***内冷负荷Qcool、综合能源***向电网购电最大功率Pgrid,max、第i种产电设备电功率下限Pi,min、***中第i种产电设备产电功率Pi、第i种产热设备热功率下限Qi,min、***中第i种产热设备产热功率Qi、***中第i种产电设备出力Ei、***中设备向上爬坡速率Smin、储能电池在t时刻荷电状态soct、储能电池自放电率δbess、储能电池充电效率ηbess,c;储能电池放电效率ηbess,d、储能电池充电状态αbess,c、储能电池放电状态αbess,d、储能电池荷电状态的上下限socmin和socmax、储能电池的起始和终止荷电状态socL和socT、电池允许的最大充放电电流Imax ch和Imax dis、***中逆变器容量Pinv、电池额定电压Vbat、电池容量Qmax、气网允许的最小压力和最大压力ρmin j和ρmax j、某种形式能源分配到不同能源转换装置的分配系数vc、园区综合能源***对外的实时售电价格Pspot、电网公布的分时电价Ptuo、对外售电价格的约束下限Qmin和上限Qmax、需求价格弹性σ;其中***内的能源价格情况如表1-1所示。
表1-1***能源价格分布表
本发明所研究的园区综合能源***内的光伏、燃气轮机、燃气锅炉、余热锅炉、电制冷机、吸收式制冷机、电储能的性能参数如表1-2所示,各个设备的功率上下限约束参数如表1-3所示。
表1-2设备性能参数情况表
表1-3设备功率参数表
S2)构建基于实时电价的需求响应模型
以分时电价为基础建立实时电价需求响应模型,实时电价计算公式如下:
根据用户需求价格弹性原理,价格型需求响应前后负荷量之间需满足关系式,如下所示:
S3)构建模型目标函数与约束条件
根据S1中的数据建立综合能源***运行优化模型,综合能源***运行优化模型以最小运行成本和碳排放量为目标,将综合能源***内部设备运行特性和负荷需求作为约束条件。
目标函数:
(1)***总成本
Min F1=Fope-Fin
(2)***碳排放量
约束条件:
(1)电功率平衡约束
Pgrid(t)+Pgt(t)+Ppv(t)+Pbess,d(t)=Pele(t)+Per(t)+Pbess,c(t),
(2)热功率平衡约束:
(3)冷功率平衡约束:
ηlbuQlbu(t)+ηerQer(t)=Qcool(t),
(4)能源供应约束:
Pe_min≤|Pe_grid|≤Pe_max
Pg_min≤|Pg_network|≤Pg_max,
(5)设备运行约束:
(6)储能电池约束:
(7)天然气连接节点压力约束:
(8)分配系数约束:
S3)求解综合能源***运行成本与碳排放量最低运行计划
综合考虑每个能源供应单元的输出功率和园区综合能源***中的***能量需求以及模型的特征,本发明采用带精英策略的非支配排序的遗传算法(NSGA-Ⅱ)求解***运行优化问题,求解得出的帕累托解Fi采用隶属度函数进行计算最终得到目标函数的隶属度值uk,根据隶属度值的大小得到优化的最佳解。
S31):优化问题描述及变量识别
实现经济成本最低、环境成本可承受是一个园区健康持续发展的必然要求,在此基础上通过建立多指标约束、多目标协同优化模型是满足该要求的关键步骤。本发明采用二进制的方式来表示设备的运行状态,设X_s=1表示设备处于运行状态,X_s=0表示设备处于待机状态,采用X_p表示设备的出力水平。综合设定的优化目标函数、优化变量的类型等情况,最后模型的优化变量可识别如下:
式中,X_p表示***内各个设备的输出功率集合;X_s表示***内各个设备的运行状态。
S32):采用NSGA-II算法对模型求解
本发明采用NSGA-II算法对模型进行求解,通过计算获得优化后的由多组解组成的非劣解集。首先随机产生规模为N的初始种群,非支配排序后通过遗传算法的选择、交叉、变异三个基本操作得到第一代子代种群;从第二代开始将父代种群与子代种群合并,进行快速非支配排序,同时对每个非支配层中的个体进行拥挤度计算,根据非支配关系以及个体的拥挤度选取合适的个体组成新的父代种群。以此类推,当满足设定的程序结束条件时输出帕累托解集。
S33):采用隶属度函数选择最优解
模糊隶属度理论通过合成每个目标函数的模糊隶属度来从非劣解集中选择最满意解。首先通过公式计算帕累托集中目标函数的隶属度;根据决策者的偏好设置权重值,并计算多目标函数的最优解的隶属权重值,最大值对应的非劣解集解是最优解。本发明中,权重设置为相等权重,隶属度值计算公式为:
式中,Fi表示目标函数的第i个目标解;Fmin i表示目标函数的最小解;Fmax i表示目标函数的最大解。N表示设定的种群数目;Di表示设定的目标函数的权重;Nobj表示设定的目标函数的数量;uk表示目标函数的隶属度值。
S4)算例优化结果分析
本发明以园区综合能源***运行成本最小化和碳排放量最低化为优化目标,将得到的基础数据与边界参数输入价格型需求响应模型,按“自发自用,余电上网”的模式,对指定算例进行仿真分析,仿真周期为一天,仿真步长为一小时。通过需求响应参与优化与不参与优化两种场景进行优化,对比分析优化目标结果以及负荷削减情况,评估模型对综合能源***的优化程度。
以下分别为需求响应模式与非需求响应模式下优化目标值对比和电力负荷削减水平对比:
表4-1需求响应前后目标值
表4-2需求响应前后电力负荷特性值
由上表和图2可看出,需求响应模式下的购电电费比非需求响应模式下的购电电费降低18.45%,而购买天然气的费用需求响应模式下仅比非需求响应模式下上涨0.24%。在实时电价的影响下,用户负荷需求在低谷时间段明显增加,在高峰时间段负荷需求明显降低,电力负荷削减程度达到47.23%。因此,本发明所研究的基于价格型需求响应模型能够很好的削峰填谷并为***运行方带来较好的收益。
以上描述的仅是本发明的基本原理和优选实施例,本领域计算人员根据本发明做出的改进和替换,属于本发明的保护范围。