摘要
利用电动汽车(EV)的移动储能特性参与建筑能源系统调度,已逐渐成为推动建筑领域迈向绿色低碳发展的重要举措之一,能够有效提升建筑的能源利用效率,同时降低建筑的运行成本。针对大型光伏建筑一体化(BIPV)公共建筑,提出一种考虑EV接入的BIPV建筑能源系统双层优化配置及调度方法。首先,利用高斯混合分布模型对电动汽车的日常出行行为进行拟合,构建基于公共建筑的EV出行规律模型。其次,为同时求解出系统的配置和调度结果,建立双层优化模型,在上层中以系统年规划成本最小为优化目标得到配置结果,在下层中以系统日运行成本最小为优化目标得到电力负荷调度结果。最后,给出双层优化模型求解方法,选取某大学校园内办公建筑采暖季的典型日天气负荷数据进行算例分析。结果表明:按照求解结果对建筑能源系统进行配置与调度能够有效提高光伏利用率,延长储能设备使用寿命。并且在系统中通过对比使用不同配置形式的储能方案发现,仅利用EV储能在长短期经济性上更具优势,能够实现良好的系统运行性能。
关键词
Abstract
The integration of mobile energy storage capabilities of electric vehicles (EVs) into building energy system scheduling has gradually become a significant measure for advancing the green and low-carbon development of the construction sector, which can effectively enhance the energy efficiency of buildings while reducing operational costs. For large-scale building integrated photovoltaic (BIPV) public buildings, a two-level optimization configuration and scheduling method for BIPV building energy systems, considering EV integration, is proposed. First, a Gaussian mixture distribution model is used to fit the daily travel behavior of EVs, establishing an EV travel pattern model based on public buildings. Next, to simultaneously determine the configuration and scheduling results of the system, a two-level optimization model is established: the upper level aims to minimize the annual planning cost of the system to obtain configuration results, while the lower level aims to minimize the daily operational cost of the system to obtain power load scheduling results. Finally, a solution method for the two-level optimization model is provided, and a case study is conducted using typical daily weather load data from the heating season of an office building on a university campus. The results indicate that configuring and scheduling the building energy system according to the solution results can effectively improve photovoltaic utilization and extend the service life of energy storage equipment. Furthermore, by comparing different energy storage schemes in the system, it is found that utilizing only EV energy storage is more advantageous in both long-term and short-term economics, and achieve good system operational performance.
2020年9月,中国在第75届联合国大会上明确提出2030年“碳达峰”与2060年“碳中和”目标[1],绿色低碳已成为建筑领域发展的重要方向[2],光伏建筑一体化(building integrated photovoltaic,BIPV)能源系统将光伏发电技术与建筑相结合,即能满足自身用电需求,又可获得能源收益,是推动中国建筑行业绿色发展的重要途径[3]。由于光伏发电的波动性,在光伏富余的情况下,大多数建筑选择直接弃光,或者配置大容量储能系统,导致严重的能源浪费与较低的经济性。电动汽车(electric vehicle,EV)具有灵活的充放电能力,将其接入建筑内作为储能、供电单元,不仅能够提高光伏利用率,还可以降低固定储能容量提升系统运行经济性[4-5]。据研究表明,EV一天内大多数时段都停放在商业、办公、住宅等建筑附近且为静止状态[6]。因此对于光伏发电计划与EV充电时间重叠的公共建筑来说,调控EV资源提升光伏利用率是行之有效的方法[7]。
调控EV的技术可主要分为车连网(vehicle to grid,V2G)与车连建筑(vehicle to build,V2B)。V2G适用于调控电动汽车集群参与电网调度[8]。而为了利用建筑能源系统中的EV资源,基于EV在建筑附近长时间驻留的特点,V2G进一步发展成为V2B[9]。V2B注重EV与建筑间的能量互动,使得EV不仅可以为建筑提供电力支持,还能在电力需求高峰期缓解电网的能源压力。Luo等[10]研究了EV接入的BIPV建筑能源系统运行的可行性。Borge-diez等[11]利用V2B技术实现降低高峰需求,节约能源成本。
据统计,2030年,中国EV数量将会达到6 000万辆[12]。若建筑物盲目地接入大量EV充电,将严重影响建筑能源安全运行,因此应在满足建筑设计要求、配网设备容量、供电可靠性的条件下建立优化配置模型,以获得最佳的EV数量与其他能源设备配置参数。目前针对接入EV能源系统的优化配置的研究中,大多数通过给定数量的EV参与电网调度作为储备电源以优化固定储能容量[13-16],未考虑接入电网中EV数量的合理性,以及由此可能对电网造成的巨大冲击和V2G实现的复杂性。He等[17]通过对比不同EV数量对商业建筑电费的影响,研究在商业建筑内接入最佳的EV数量以实现建筑能源系统管理成本的降低。Barone等[18]在建筑能源系统调度中研究EV的接入对系统整体的影响,结果表明,EV的移动储能特性能够有效减少电网供电并提升系统经济性。Ahsan等[19]通过提出一个多目标经济优化框架,实现了建筑物与EV调度的经济性最优。上述文献注重V2B中调度短期经济性研究,未将短期经济性与系统投资长期经济性结合起来考虑。
针对BIPV建筑能源系统,本文将接入建筑的EV规模纳入建筑能源系统的设计参数,从接入EV数量的合理性、提升光伏利用率和优化储能容量的角度出发,提出一种基于V2B模式下的能源系统双层优化框架的最优配置及调度策略。首先,建立BIPV建筑能源系统模型,并通过真实出行数据构建基于公共建筑的EV行驶特性概率密度函数模型。其次,采用上、下两层相互传递参数多次迭代求解的思路建立双层优化模型,在上层模型中设置系统年规划成本最小为优化目标求解出建筑的配置参数,下层模型中设置系统日运行成本最小为优化目标求解出电能调度结果。最后,在某大学校园内的办公建筑中应用所提策略验证有效性,并且对EV接入建筑能源系统的长、短期经济性进行了分析。
1 BIPV建筑能源系统结构
本文所构建的BIPV建筑能源系统主要由EV集群、光伏系统、储能系统、空气源热泵和热储能系统组成,系统内能量流动见图1。图1中系统主要电能由光伏系统与电网供应,蓄电池存储富余电能,EV作为移动储能,帮助蓄电池存储超过其存储能力上限的电能,热能由空气源热泵和蓄热水罐提供。
图1系统能量流动示意
Fig.1Schematic diagram of system energy flow
2 BIPV建筑能源系统模型
2.1 EV模型
EV的行驶里程与充电行为主要受到车主的日常出行规律影响,其出行规律可近似认为与传统燃油汽车相同。因此,根据美国国家公路交通安全管理局(national highway traffic safety administration,NHTSA)[20]对全美家用车的调查统计数据结果,利用高斯混合分布模型对EV抵达办公地点时刻和离开办公地点时刻进行拟合,以获取其概率密度函数。
2.1.1 抵达办公地点时刻
抵达办公地点时刻为
(1)
式中:x为当日到达办公地点时刻,σA1、σA2为标准差,分别取值为1.39、14.89;μA1、μA2为期望值,分别取值为7.71、13.20;α1、α2为相应高斯分布的权重,分别取值为0.74、0.26。EV抵达办公地点时刻的概率分布见图2。
图2抵达办公地点时刻概率分布
Fig.2Probability distribution of arrival time at the office
2.1.2 离开办公地点时刻
离开办公地点时刻为
(2)
式中:x为当日离开办公地点时刻,σL1、σL2为标准差,分别取值为1.16、26.25;μL1、μL2为期望值,分别取值为16.86、15.51;β1、β2为相应高斯分布的权重,分别取值为0.53、0.47。EV离开办公地点时刻的概率分布见图3。
2.1.3 日行驶距离
考虑到EV日行驶距离与每百公里耗电量,根据文献[21]分析可得,EV到达办公地点的初始荷电状态这一随机变量近似服从参数为U(0.3,0.8)的均匀分布。
图3离开办公地点时刻概率分布
Fig.3Probability distribution of departure time
2.2 空气源热泵模型
空气源热泵通过消耗电能能够实现制热的效果,其输出功率受室外温度和热泵额定功率的影响,设备模型可由下式给出的数学模型表示为[22]
(3)
式中:QHP(t)为t时刻空气源热泵的输出功率,CH(t)为t时刻空气源热泵的性能指数,PHP(t)为t时刻空气源热泵的耗电量,κ为空气源热泵结霜除霜损失系数,TAIR(t)为t时刻室外空气温度。
2.3 光伏发电模型
光伏组件由串联的光伏电池组成,多个光伏组件串并联形成光伏板,整个光伏板t时刻的输出功率模型为
(4)
式中:PS(t)为t时刻光伏板的实际输出功率,NS为一个光伏组件中光伏电池串联的数量,NSS为光伏组件串联的数量,NPP为光伏组件串并联的数量,ICELL(t)为光伏电池t时刻的输出电流,VCELL(t)为光伏电池t时刻的输出电压。
上述式中的参数取决于太阳辐射照度、光伏电池材料、环境温度和光学空气质量,在文献[23]中能够获得这些参数的详细数据。
2.4 储能系统模型
在BIPV建筑能源系统中蓄电池的使用能够实现削峰填谷的作用,在负荷低谷时段充电、高峰时段放电,其充、放电时的荷电状态为
(5)
式中:Sbat(t+1)、Sbat(t)分别为t+1和t时刻蓄电池的荷电状态,v为蓄电池的自放电效率,、分别为t时刻蓄电池的充、放电功率,ηc为蓄电池充电效率,ηd为蓄电池放电效率,EB为蓄电池额定容量。
2.5 热储能系统模型
蓄热槽用于存储空气源热泵产生的热能,在用户能耗需求高峰时段释放能量,蓄热与放热的数学模型为
(6)
式中:QXR(t+1)、QXR(t)分别为t+1和t时刻的蓄热量,σXR为热损失系数,PXR,IN(t)为t时刻的蓄热功率,PXR,OUT(t)为t时刻放热功率,fXR,IN为蓄热效率,fXR,OUT为放热效率。
3 双层优化模型及其求解方法
3.1 双层优化模型
为同时求解出BIPV建筑能源系统的配置参数与调度策略,设计了一种双层优化模型。在上层模型中根据建筑自身配电容量与各设备运行的约束条件,以年规划成本经济性最佳为优化目标确定光伏组件数量、蓄电池容量和EV数量;下层模型根据系统内部能量平衡等约束条件,以日运行成本最小为优化目标得到系统内各设备的调度策略。选择典型日的天气负荷数据作为双层优化模型的输入量,上层模型求解出的配置参数作为下层模型的输入量,下层模型求解出的日运行成本作为上层模型的输入量,上、下两层模型通过传递求解变量以得到系统配置与调度策略的最优解,双层优化模型的运行框架及具体描述见图4。
图4双层优化模型运行框架及具体描述
Fig.4Two-level optimization model operation framework and specific description
3.1.1 上层目标函数
为获得系统内光伏、蓄电池和EV充电站接口数量设定的最优配置方案,以系统年规划成本最小为目标函数,表达式如下:
(7)
其中:
式中:FPV、FEV、FBA、FOP分别为系统内光伏发电年投资成本、EV充电桩年投资成本、蓄电池年投资成本和系统日运行成本,PEc(i,t)、PEd(i,t)分别为 t时刻第i辆EV的充、放电功率,cEV、dEV分别为单位时间EV充电支出和放电收益,nEV为EV数量,CRF为系统设备资本回收系数[24],i为利率,l为设备的使用寿命,CRFPV为光伏发电资本回收系数, cPV为光伏发电资本支出,nPV为光伏板发电所需块数,CRFBA为蓄电池资本回收系数,cBA为蓄电池资本支出,为蓄电池额定容量,WEV为单位EV充电接口建设资本支出。
3.1.2 上层约束条件
3.1.2.1 光伏发电所需组件数量约束
光伏发电所需组件数量约束为
(8)
式中:、和DPS分别为不同典型日天气条件下单位光伏组件发电功率的最小值、最大值和平均值,APROOF为屋顶建筑面积。
3.1.2.2 蓄电池额定容量约束
蓄电池作为系统的备用电能供应源,其额定容量应满足在外部能源供应中断等极端条件下的用电需求。为保证系统的安全运行,蓄电池额定容量约束表达式为:
(9)
(10)
式中ELOAD(t)为t时刻的负荷需求。
3.1.2.3 EV数量约束
为避免由于EV充电负荷与基本用电负荷叠加造成电力系统容量过载,EV充电接口数量应满足如下约束表达式:
(11)
式中为电力系统容量最大值。
3.1.3 下层目标函数
上层配置问题与下层调度问题间通过系统日运行成本相互耦合。下层模型在考虑到分时电价与EV充电需求的情况下,依据上层模型确定的光伏组件数量、蓄电池容量和EV充电桩数量,以日运行成本最小为目标函数制定系统内设备出力调度计划,表达式如下:
(12)
其中:
式中:FECO、FOPE,i、FENV分别为系统运行的经济成本、维护成本和环境惩罚成本,n为系统中设备数量,为t时刻购电价格,为t时刻售电价格,为t时刻电网向系统负荷输出电能,为t时刻系统向电网回馈电能,cOPE,i为第i个设备单位维护成本,PGL,i(t)为第i个设备t时刻的输出电能,i为系统中设备编号,j为污染物(CO2,SO2,NOX)种类,λj为设备污染物环境价值,μj为设备排放污染物所受惩罚费用,Pi(t)为t时刻各个设备出力功率,Rji为设备污染物排放参数,RjG为电网污染物排放参数。在本文研究中,污染物排放来源主要由电网运行消耗的化石能源产生。
3.1.4 下层约束条件
3.1.4.1 系统功率平衡约束
为维持系统内部各设备间电能和热能的实时平衡,系统电、热功率平衡约束如下。
①电功率平衡:
(13)
②热功率平衡:
(14)
3.1.4.2 系统设备运行约束
①蓄电池储能约束:
(15)
(16)
(17)
(18)
(19)
(20)
式中:是一个布尔变量,“1”表示t时刻蓄电池充电,否则为静置状态; 是一个布尔变量,“1”表示t时刻蓄电池放电,否则为静置状态;、分别为蓄电池充、放电功率的最大值,、分别为蓄电池充、放电功率的最小值,为蓄电池荷电状态的最小值,Sbat(0)为蓄电池初始荷电状态。
②空气源热泵约束条件为式(3)所示。
③蓄热槽储热约束:
(21)
(22)
(23)
(24)
(25)
式中:为一个布尔变量,“1”表示t时刻蓄热槽蓄热,否则为静置状态;为一个布尔变量,“1”表示t时刻蓄热槽放热,否则为静置状态;、分别为蓄热槽蓄热、放热功率的最大值,、分别为蓄热槽蓄热、放热功率的最小值,、分别为蓄热槽容量的最大值、最小值。
3.1.4.3 EV充、放电约束
当EV参与系统运行调度时,不仅需要引导EV有序充、放电,还需要能够在调度时段内达到车主满意的电池荷电状态。因此,对于EV应从电池的充、放电状态和车主荷电状态满意度的角度设定相应的约束条件。
①荷电状态需求约束:
(26)
(27)
(28)
(29)
(30)
式中:SEV,END为EV在调度周期结束后的电池荷电状态,SEV,ARR为EV接入系统前电池的荷电状态,xc(t)、xd(t)分别为t时刻的充、放电状态标志位,SEV,DES为EV接入系统离开后车主期望最小荷电状态,SEV,max为EV电池荷电状态的最大值,、分别为EV充、放电功率的最大值,CapEV,BAT为EV蓄电池额定电量,、分别为EV充、放电效率,TARR为EV接入系统的时刻,TDEP为EV断开连接的时刻,ΔT为EV调度期间内所用时间。
②荷电状态约束。EV在t时刻电池的荷电状态与前一时刻EV的充、放电功率和电池电量有关,为延长电池使用寿命,应在每一个时刻内对其荷电状态设置合理的范围,避免调度周期内出现过度充、放电的现象。
(31)
(32)
式中:SEV(t+1)、SEV(t)分别为t+1和t时刻EV电池的荷电状态,SEV,min为EV电池荷电状态的最小值。
③充放电状态约束。当EV通过充电桩接入系统,在一个调度时间段内其不能同时进行充、放电行为,当EV与系统断开连接,其充、放电标志位必须全部为0。
(33)
3.2 双层优化模型求解方法
BIPV建筑能源系统双层优化模型的求解是考虑约束条件的混合整数线性规划问题,模型由于上、下两层模型间求解变量相互耦合,直接求解难度较大,并且模型中存在0/1变量,难以通过KKT条件将下层模型合并到上层模型中,因此本文采用文献[25]中的双层优化算法的思路求解该模型,并且设置多次迭代求解次数以避免求解结果可能出现的随机性。在MATLAB中对上述算例进行仿真,上、下两层模型均利用Yalmip工具和CPLEX求解器进行求解。模型求解流程见图5。
图5模型求解流程图
Fig.5Flowchart of model solving
双层优化问题求解步骤具体描述如下所述。
Step1 设置系统总迭代次数m。
Step2 初始化系统参数,随机生成光伏板块数、EV数量和蓄电池容量,将参数传递至下层模型。
Step3 更新迭代次数k=k+1。
Step4 下层模型接受系统配置参数,利用CPLEX求解器求解系统设备出力和与电网交互的功率以及EV的充、放电功率,记录系统日运行成本花费并将其返回至上层模型中。
Step5 上层模型接受系统日运行成本花费,利用CPLEX求解器求解生成新的一组系统配置方案。
Step6 判断迭代次数是否达到最大值,若达到最大迭代次数则检查结果是否收敛,否则返回Step3中,重新迭代求解。
Step7 判断结果是否满足收敛条件,若满足条件则输出系统最优配置方案,否则返回Step1中,重新开始模型求解。
4 算例分析
4.1 算例参数
本文以某地大学校园内一栋办公建筑采暖季的典型日负荷数据和天气数据为例,使用图1所示的BIPV建筑能源系统结构作为算例的研究场景。在MATLAB平台上进行求解验证,设置系统调度周期为24 h、步长为1 h的仿真算例。为避免求解过程中可能会陷入局部最优,本文设置程序运行200次,采用场景消减方法将全年的负荷和天气数据缩减为30个典型日场景,选取典型日场景运行成本的平均值作为评估系统长期经济性的参数,每次求解时间大约15 min,图6为全年天气负荷数据曲线。建筑能源系统的部分设备参数见表1,设备运行维护成本[26]见表2,设备污染物排放参数、污染物环境价值标准和罚款数量级[27]见表3、4,电网分时电价[28]见表5。
图6全年天气和负荷曲线
Fig.6Annual weather and load curve
表1设备参数
Tab.1Equipment parameters
表2设备运行维护成本
Tab.2Equipment operation and maintenance cost
表3设备污染物排放参数
Tab.3Pollutant emission parameters of equipment
表4污染物环境惩罚费用
Tab.4Environmental penalty costs of pollutants
表5电网分时电价
Tab.5Power grid time-of-use price
4.2 结果对比分析
经过上、下两层模型多次迭代求解后,检验求解结果数值的稳定性。EV数量、光伏组件数量和蓄电池容量迭代结果见图7。从图7中可以得到光伏组件数量的迭代结果为933块,储能容量的迭代结果为61 kWh。EV数量由于其充电行为的不确定性导致迭代结果出现波动,但数值整体波动在14附近。因此假设EV的迭代结果数值为14,经过计算验证得到误差的平均值为0.222 7,小于1满足误差范围,系统最优配置方案见表6。
图7建筑能源系统设备参数优化结果
Fig.7Optimization results of building energy system equipment parameters
表6系统最优配置方案
Tab.6Optimal system configuration scheme
在系统最优配置方案并且其他设备参数不变的情况下,为评估EV对系统日运行调度的影响,本文首先设置两种配置场景进行对比分析。场景1为不考虑EV参与系统调度,场景2为考虑EV参与系统调度,其系统日运行调度分别见图8、9。场景1、2下蓄电池的日充电次数和光伏利用率见表7。经过对比分析得出,EV参与建筑能源系统调度,不仅提高了光伏发电的利用率,还减少了蓄电池的充放电次数。由于蓄电池的使用寿命是固定的,接入EV的配置策略有助于延长蓄电池的使用寿命。
图8储能配置为蓄电池的日运行调度结果
Fig.8Daily operation scheduling results of battery energy storage configuration
图9储能配置为EV和蓄电池的日运行调度结果
Fig.9Daily operation scheduling results of EV and battery energy storage configuration
表7EV接入对建筑能源系统运行性能的影响
Tab.7Effect of EV access on the operational performance of building energy systems
基于上述结论,为进一步分析EV接入建筑能源系统中带来的潜力,在系统中设置仅依靠EV储能的第3种场景,系统日运行调度见图10。通过对比场景1、场景2和场景3的调度结果,验证利用EV代替固定储能方案实施的可行性,表8详细展示了场景1、场景2、场景3下电能来源、年投资成本、日运行成本和年总成本的运行结果。在场景2中,由于系统中同时配置了固定储能与移动储能,导致电网电能消耗量高于场景1和场景3,增加了系统对外部电能的依赖,不利于系统的安全稳定运行。 EV接入建筑能源系统中不仅可以充分消纳富余光伏,还能够为建筑提供收益,因此随着建筑能源系统中EV承担储能能量的提升,系统的日运行成本将逐渐降低。场景2的日运行成本比场景1降低了8.6%,而场景3的日运行成本比场景1降低了16.7%。从系统长期经济性的角度分析,文中61 kWh蓄电池的年最低投资成本将达到2.2万元,14根充电桩的年投资成本为0.36万元。相比于只配置成本高昂的蓄电池的场景1,场景2的年总支出成本降低了6.3%。而场景3的年总支出成本仅为18.56万元,比场景1降低了22.8%。从场景3的运行结果可以看出,EV作为移动储能设备,与建筑之间进行能源的双向互动,不仅能够满足EV自身的充电需求,还有效地帮助建筑降低了运行成本,实现了节能减排的目标。EV接入建筑能源系统不仅在短期内具有明显的经济优势,而且在长期经济性上也表现出色,场景3的配置策略在经济效益上明显优于场景1和场景2。
图10储能配置为EV的日运行调度结果
Fig.10Daily operation scheduling results of EV energy storage configuration
表8不同储能配置下建筑能源系统运行调度结果
Tab.8Operation scheduling results of building energy system under different energy storage configurations
在建筑能源系统中应用场景3的配置方案,EV既能够缓解建筑用电高峰压力,又可以在电网负荷较低时进行充电,平衡电力供应。这种能源的双向互动,不仅提升了能源利用效率,还促进了可再生能源的消纳,有助于实现更环保的能源结构,因此考虑将EV接入建筑能源系统以取代传统蓄电池储能是一种前景广阔的策略。
5 结论
本文在“双碳”背景下针对建筑能源系统中接入不合理EV数量造成的不良影响,提出一种考虑EV接入的BIPV建筑能源系统最优配置及调度的策略,建立了系统配置和调度的双层优化模型,通过仿真试验得到的主要结论如下:
1)为确定电动汽车接入V2B充电桩后可调度的时间段,采用高斯混合分布模型构建基于公共建筑的EV用户出行概率密度函数,经过仿真验证,所建立的EV模型符合出行行为模型构建的要求,从而提高了双层优化数学模型建立的合理性。
2)本文建立的双层优化模型采用上、下层模型传递求解变量并通过多次迭代求解的策略。经过反复迭代后,对运行结果进行稳定性验证,结果表明,在不同典型日的求解中,结果表现出高度一致性,未出现较大偏差或不稳定情况。这充分证明了所提出的双层优化模型的合理性及其求解方法的有效性。
3)本文通过对接入EV的BIPV建筑能源系统应用所提双层优化模型及求解方法,经过仿真证明,该策略不仅能够提升光伏发电的利用率,还能延长蓄电池使用寿命。同时相比于在建筑能源系统中仅依靠固定储能,或混合使用固定储能和移动储能的方案,单独利用EV储能在长短期经济性上展现出更为突出的优势,使得系统整体运行更加优化,实现更长的设备使用寿命和更高的经济效益。
4)随着电动汽车的广泛应用,本文所提的最优配置及调度双层优化模型将为投资成本高昂、可再生能源利用率低的BIPV建筑能源系统提供可行的解决方案,具有良好的应用价值与现实意义,为实现低碳经济和可持续发展提供有力支持。

