需求响应调度寻优策略：是否存在使得电网侧调度策略实施最优的日最佳参与率？

针对智能电网环境下,各类不确定因素对电网需求响应调度策略的实施效果存在较大影响的问题,提出一种计及用户参与不确定性的用户响应模型,将用户参与度和电网侧激励间的关系用不确定性模型来表征。考虑到电网侧对负荷削峰填谷的要求,构造了基于负荷波动程度的需求响应目标函数。对于参与响应的用户部分,其调度策略采用有序用电调度,无响应用户部分采用无序用电。算例通过实验仿真,分析验证了在计及参与不确定性的需求响应模型下,用户侧存在着使得电网侧调度策略实施最优的日最佳参与率,也即电网侧对应存在着最优激励政策。实验仿真结果表明了计及用户参与不确定性的需求响应优化策略具有有效性和实际优越性,能为电网侧提供更全面的用电信息,电网公司从而能够针对性提高需求响应策略实施效果。

引言

智能用电是实现电网与用户双向实时交互的重要环节[1],能够满足用户的多元化用电需求,实现灵活互动的供用电新模式,是未来能源互联网背景下智能用电技术发展的趋势[2]。随着能源互联网的兴起,智能用电的规划、建设与管理运维越显重要。为加快能源互联网建设,推动互联网与电力系统各领域深度融合和创新发展,国网公司开展了全球能源互联网前瞻研究[3-8]。需求响应(demand response,DR)作为用电环节与其他各环节实现协调发展的关键支撑手段,是智能用电互动化业务中体现电能量友好交互的重要方式。但在需求响应策略实施过程中,由于用户可能存在的对激励缺乏重视、通信延迟和消费行为改变等一系列原因,用户对激励、价格信号的实际响应程度,甚至是对需求响应策略采纳的意愿往往都带有较大的不确定性[9]。考虑到采用确定性的DR模型来描述需求响应已不能满足智能用电协调运行需要[10],因此,开展计及不确定性的需求响应机理和模型,以及需求响应实施策略具有必要性。在智能电网中大规模新能源接入的环境下,调用需求响应资源参与系统调峰以应对大规模新能源接入,对电网的稳定运行具有重要意义。

当前考虑不确定性的需求响应模型及策略相关研究正逐渐成为当前的研究热点[11-16]。文献[11]构建了计及不确定性的用户响应模型,通过负荷削减能力和响应的不确定性两个特征描述需求响应随经济激励及其他因素变化的规律;文献[12]引入多智能体模拟方法,建立双向互动的改进分时电价调价模型,分析了电力公司与用户的双向互动过程,以及用户响应的不确定性对电力公司分时电价设定的影响;文献[13]针对负荷集群对电价响应的不确定性,提出一种联合调度模型,保证调度经济性的同时减少需求侧因负荷响应不确定带来的风险损失;文献[14]考虑价格型需求响应、风电出力和系统负荷的不确定性,建立了考虑风电消纳效益的电力系统源荷互动日前模糊优化调度模型;文献[15]建立了考虑风电输出不确定性的多级机组需求响应曲线模型,将电价表示成期望值和随机偏差之和的形式,最后用Benders分解进行求解;文献[16]将负荷响应量设为某区间的随机变量,并考虑权重约束,确保负荷率风险损失低于预估水平,建立随机优化模型,用抽样平均近似方法求解。在上述文献中,文献[11]和文献[12]是针对用户负荷削减的不确定性;文献[13]是针对负荷集群对电价响应的不确定性;文献[14]是针对负荷响应率的不确定性研究;文献[15]和文献[16]是针对价格弹性需求曲线的不确定性研究。

以上不确定性需求响应研究,多假设参与需求响应用户数为确定量,尚未考虑用户是否选择参与需求响应调度策略的确定性因素,因此本文针对需求响应策略中用户参与不确定性因素,首次在提出反映电网侧激励对用户用电影响的参与率(participation rate,PR)概念基础上,构建了计及用户参与不确定性的需求响应模型;同时以配电网侧波动性最小为目标函数,提出一种计及用户参与不确定性的DR调度策略优化方法,最后通过实验仿真验证了建模和优化策略求解的有效性。

1 不确定性需求响应建模

1.1 计及响应参与率的不确定性需求响应机理模型

根据消费者心理学原理,电网给出的激励政策对每个用户的刺激均存在着一个差别阈值及饱和值[17]。而针对小区用户整体而言,在不同的激励水平下,实际参与策略人数同样存在理论上限及下限(即存在一个不确定区间)。对单个用户来说,不确定性体现为是否响应DR策略,以及其响应的负荷削减比重;而以智能用电小区全体用户为对象来说,不确定性可体现为整个小区的DR参与率,即接受调度策略,参与响应的用户数占用户总数的百分比。

受文献[11]中的用户负荷削减量响应曲线模型启发,本文建立线性模型近似来表征用户参与响应的行为,如图1所示。

图1 计及不确定性的用户响应参与率曲线

Fig. 1 Curve of customer response participation rate considering uncertainty

图1中激励水平用电费折扣来表述,在某一激励水平x下,响应参与率R的上下界可分别表示为

式中：Rup(x)Rup(x)表示响应参与率的上限函数;Rdown(x)Rdown(x)表示响应参与率的下限函数;x0表示临界激励水平,激励水平在临界点x0之前时,参与率表现为死区,即激励水平达不到用户整体对策略有所响应的程度;x1表示理想饱和激励水平,当激励水平达到x1时,参与率达到理想饱和点,即参与率上限刚好可到达100%;x2表示完全饱和激励水平,当激励水平达到x2时,参与率达到完全饱和点,即参与率不确定波动区间消失,参与率成为确定量100%。

考虑到激励水平与电网侧的成本呈正相关,结合概率统计学理论,在最佳参与率RoptRopt下的策略最佳激励值xoptxopt可表述为

xopt=12[xmin+xmax] (3)

结合公式(1)与公式(2),公式(3)可化为

xopt=12[Ropt(x1+x2)]−x0⋅Ropt(4)

式中：RoptRopt表示最佳参与率;xmin与xmax分别表示在最佳参与率RoptRopt下的激励水平区间[xmin,xmax]的区间最小值和最大值,如图1中所示;xoptxopt表示在最佳参与率RoptRopt下的最佳激励水平。

本文构造的响应模型采用某一区间内的随机数来表征用户参与响应这一不确定行为,通过区间参数的选择,对策略实施情况进行多次模拟,从而尽可能在DR策略允许的参与率范围内,找到相对较优的激励水平,使得在较低激励水平条件下,调度策略还能达到较好的效果,从而使得策略的适用性提高。用户通过参与响应的不确定度来实现与电网侧的浅层互动形式,电网侧通过控制激励水平大小,进而影响用户整体对调度策略的响应程度,最终尽可能得到实现调度目标函数最小的最佳DR参与率。

1.2 目标函数

对于电网侧来说,其进行调度优化的主要目的是削峰填谷,实现整体负荷波动的最小化,从而提高经济生产效益。综上分析,目标函数的建立可考虑与DR策略调度后整体负荷波动性相关,本文按日分成T=24个时段,每个时段为1 h,则智能用电小区内n个用户第t个时段内用电总负荷Q(t)Q(t)可表示为∑k=1nQk(t)∑k=1nQk(t),每日内平均用电总负荷Q(t)¯¯¯¯¯¯¯¯¯¯Q(t)¯为Q(t)Q(t)的数学期望,即1T∑t=1TQ(t)1T∑t=1TQ(t),借鉴标准方差标示测量值的波动程度,本文选择小区日用电总负荷的标准方差作为调度后整体负荷波动性的度量,故本文所采用优化目标可表述为

min { 1T∑t=1T(Q(t)−Q(t)¯¯¯¯¯¯¯¯¯¯)2} (5)

本文考虑的响应参与不确定性实则也同时反映了用户满意度,由于是否响应策略是由用户自主选择的,因此接受或不接受激励政策的用户应均对用户日常生活的用电体验影响不大,在DR策略中,用户侧向电网侧的反馈这一环节体现的就是用户满意度因素,所以可认为目标函数从电网侧考虑是合理的。

1.3 用户用电调度模型

本文用户的负荷构成采用普通的家庭用能负荷,仅考虑基本负荷以及可中断负荷,未涉及其他新型智能用电负荷。基本负荷为居民生活必需的负荷部分,如照明、烹饪等;可中断负荷是可以根据实际需要灵活调整用电量及用电时段的部分,即可进行负荷削减或者转移,如空调、电视等家电,也有文献称为可调度负荷[18]。

以小时为单位将1 d划分为24个时段,则智能用电小区n个用户1 d内的负荷安排为

D=[Q1,Q2,⋯,Qk,⋯,Qn (6)

式中：D为24×n维的矩阵,表示总负荷安排;QkQk为24维向量,表示用户k在1d内的负荷安排;n为智能用电小区中居民用户总数。

用户k的日用电负荷QkQk可描述为

Qk=QBk+QSk (7)

式中：QBkQkB表示用户k的在t时刻的基本负荷;QSkQkS表示用户k在t时刻的可中断负荷部分。

基本负荷的约束条件为。

QBk=Qmink (8)

式中QminkQkmin表示日前预测负荷[19],可根据居民日常用电情况来获取。

可中断负荷作为可调控负荷,应满足如下约束条件：

∑t=1TQSk(t)=zSk (9)

0≤QSk(t)≤Qmaxk(t)(10)

式中QSk(t)QkS(t)表示用户k在t时刻的可中断负荷,同时QSk(t)QkS(t)还要满足1 d的总量约束zSkzkS。考虑到实际中存在用电保护约束,该负荷还须不超过最大调度量Qmaxk(t)Qkmax(t)。

2 策略寻优求解流程

2.1 算法策略流程

本文DR策略优化采用分布式求解方法,目标函数求解算法采用混合粒子群算法[20] (hybrid particle swarm optimization,HPSO),通过算法对不同用户参与率情况下的DR策略进行遍历寻优,找到最佳参与率,对应机理模型可得到电网侧最佳激励水平区间。

求解过程分为响应策略部分与非响应部分,算法流程示意图如图2所示。

算法整体流程具体包括以下几步：

1)首先,进行负荷的初始化。用户可以将其1 d的负荷使用计划上传至电网侧。

图2 算法流程示意图