信息物理融合的智慧能源系统多级对等协同优化

司方远 汪晋宽 韩英华 赵强

随着分布式可再生能源技术的推广及售电侧市场逐步开放,电力能源在生产、传输、存储、消费和交易等各个环节呈现出信息与物理深度融合的特点,电力网络基础设施作为物理系统便于与集通信、计算和控制于一身的信息系统无缝集成和互动.因此,信息物理融合系统(Cyber-physical systems,CPS)[1-2]正是上述特定集成和交互问题的解决方案,而信息物理融合能源系统(Cyber-physical energy systems,CPES)[3]已成为当前工业与能源领域研究的热点问题,并推动着能源电力系统从数字化、信息化向智能化、金融化的方向发展.

目前,微电网、智能电网以及能源互联网都是典型的CPES系统.其中,微电网和智能电网取得成功的关键在于利用CPES实现信息系统与电力能源物理系统的高效互动,以有效提供经济、安全和可持续的能源供应[4-6].而能源互联网的技术需求源于能源供需矛盾和新型大规模可再生能源的接入,同时强调通过能量枢纽和能源路由器实现能源流、信息流和业务流的深度融合[7-9].此外,能源的本质是服务于人类社会的.反过来,人类活动直接参与调控能源生产、消费与交易过程.为此,Zhang等[10]提出并诠释了社会能源概念,结合能源的生产、消费过程和环境保护等方面深入剖析了人类活动与能源系统的交互作用.考虑这种交互,人类将有机会向社会挖掘第五燃料—节能.文献[11]从复杂系统理论角度结合知识自动化理念进一步提出了基于信息—物理—社会融合系统的能源5.0.因此,能源系统迫切需要与新的CPS技术相融合,以构建满足人类社会需求的智慧能源系统(Intelligent energy systems,IES)[12-14].IES系统的精髓是智慧,是集信息网络和能源物理网络于一体,综合考虑能源的社会和技术双重属性,拥有自组织、自协调、自治理等功能,具有开放、共享和即插即用等特征,并满足安全、经济和可持续等社会需求的未来能源形式,是智慧时代背景下能源革命的迫切需要.

在此背景下,分布式能源管理方法近年来受到了广泛关注和深入研究.Roche等[15]对用于电网能源管理的几种多智能体系统进行了对比分析.Duan等[16]采用协同合作多智能体系统,针对电力基础设施中利益相关者的互操作性,提出了一种未来电力零售市场运行机制.Vergados等[17]将大量电力消费者聚类分别作为单一虚拟集群参与市场,并通过减少负荷预测的不确定性来降低总能源成本.文献[18]基于Nikaido-Isoda松弛算法进一步实现了零售电力市场的一般框架.分布式能源通过考虑纳什均衡的价格竞标策略参与电力市场,从而最大化售电利润.

上述文献虽然以双边或多边等形式考虑未来自由化的电力市场运行,但涉及的能源实体均为单纯的发电侧或单纯的负荷侧,且交互过程受双方调度中心的制约:1)大量发电或负荷实体的接入将导致调度中心的计算复杂度增加、运行效率降低,难以满足电力资源优化配置的高效性和实时性;2)调度中心拥有发、配电设备的完全控制权,导致能源的社会属性受到制约,难以满足用户的个性化需求;3)调度中心的信任机制依赖于第三方背书,交易实体与中心机构间容易产生信任问题,难以确保能源交易的有效性与公平性;4)调度中心存在信息安全隐患,中心机构一旦受到攻击,其树形结构将导致大量信息泄露及整个调度系统的瘫痪,难以维护用户的隐私和能源交易的安全性.

针对上述问题,新兴的弱中心化策略为分布式智慧能源系统能量优化调度提供了可行的技术路径.作为虚拟货币核心技术的区块链技术[19-20],其特点可以归纳为:分布式、自治、可编程、可追溯、去信任、按照智能合约执行等[21-22],所强调的分布式特性与分布式电力交易系统的构造不谋而和,可编程的智能合约可解决电力能源灵活交易执行等问题.另外,区块链的加密机制和算法可维护产消者的隐私和电力交易的安全性,进一步实现用户友好的电力能源交易模式.

此外,IES系统的构建还需兼顾可再生能源和智能负荷的整合配置以及更加灵活的分布式能源优化调度和交易机制.因此,构建集能源生产者、能源消费者、能源买家以及能源卖家多重角色于一体的产消者(Prosumer)[23-25]成为实现IES系统分布式能源管理的可行解决方案.然而,如何构建新型产消者市场分布式交易模型,以及采用何种求解策略目前还是一个开放性的问题.

针对这一问题,Zhang等[26]利用双边非合作博弈竞价搭建了产消者对等能源交易平台,结果表明基于产消者的对等能源交易形式有助于能源的合理利用以及改善本地能源的供需平衡.文献[27]基于对偶内点法求解实时电力市场平衡交易问题,结果显示不同的控制参数对算法收敛性影响较大.针对含风电场的电力系统的优化调度问题,文献[28]提出了含风电场的实时调度模型,重点研究了混沌量子粒子群优化方法(Chaos quantum-behaved particle swarm optimization,CQPSO),结果表明引入混沌初始化和混沌扰动可有效避免陷入局部最优.针对含多种分布式能源的中压微电网优化调度问题,文献[29]提出了考虑经济和环境等因素的协同优化模型,重点研究了带有自适应权重和加速度系数的粒子群算法优化方法(Particle swarm optimization with adaptive weight and acceleration coefficients,AWCPSO),有效地提高了算法的收敛精度.

基于上述分析,如何构建产消者市场的分布式交易模型,以及如何有效应用模型信息指导优化算法的搜索过程以提高计算效率,是IES系统优化运行亟待解决的问题.本文提出了一种信息物理融合的智慧能源系统多级对等协同优化方法.基于经济学Stone-Geary效用函数和双向拍卖机制构建了分布式产消者群体电力资源优化分配的效用函数.采用邻域再调度粒子群(Neighborhood re-dispatch particle swarm optimization,NR-PSO)算法[30],通过收敛判定域的引导并结合全局随机寻优与区域定向寻优的策略,进一步提高了算法的收敛速度,在最大化社会福利的同时确保实时电力供求平衡.此外,为解决分布式电力交易的信任问题并维护产消者的隐私和电力交易的安全性,引入可编程的智能合约技术进一步实现用户友好的电力交易模式.最后,仿真实例验证了本文方法的有效性和可行性.

1 问题描述

1.1 问题分析

IES具有以下组合特征:1)混合.如前所述,IES是能源电力系统与CPS的结合.因此,在建模、设计、优化和分析的过程中需要有效的理论来推理动态的、离散的或连续的混合系统.2)异构.由于IES的组件的种类和结构不一,系统需要考虑多个网络平台和不同计算模型接口间的互操作性.3)分布式.在IES架构下,物理层包含了大量的分布式产消者或产消者联盟;信息层则需要提供相应的分布式计算,例如:并行计算以及云计算等社会计算[31].而金融层主要由分布式智能合约及分布式账本构成,以实现分布式对等电力交易服务.4)大规模.由于产消者不断向IES渗透,系统的规模不断扩张,终将导致传感器、执行器、计算和通信设备的集群互联并产生海量相关数据.5)动态.IES所处的能源政策等社会环境在不断发展,因此系统的设计和运行必须考虑在这种动态环境下的自适应能力.6)人在环路(Human-in-the-loop).在IES系统电力市场交易中,产消者作为理性人可根据自身需求产生理性的购电意愿,同时对交易的电能理性报价.因此,IES系统有必要将人在环路的角色和接口设计视为一个核心方面.

1.2 IES多级对等协同优化架构

基于上述特征,本文提出了IES多级对等协同优化架构,如图1所示.由下至上将系统分为三层网络[32-33]:能源物理系统网络层、能源信息系统网络层以及能源金融系统网络层.

在能源物理系统网络层,每一个电力用户被视为一个产消者单元,由智能电表、能量枢纽、可再生能源发电装置(如家用小型风力发电机、光伏面板)、储能设施(如蓄电池)以及智能负荷(如智能家用电器、电动汽车)等组成.其中,每个用户与其他用户相比,各自的负荷特征存在较大差别,且呈现出明显的波动性,而在小范围内一天中可再生能源出力被认为是等量的,这种供需差额将导致用户产生购/售电意愿,且二者呈显著的正相关,即差额越大,用户的意愿越强烈.因此,有必要将若干个用户结成产消者联盟(图1),通过局域协调,充分满足用户意愿.然而局域协调很难实现联盟内能源供需平衡,导致整个联盟产生向其他联盟或主电网购/售电的意愿,因此需要通过广域协调进一步最小化能源供需的整体偏差.其中,能源路由器可作为能源聚合器,完成局域协调后电能的导出或导入.

在能源信息系统网络层,基于现有通信技术或下一代通信技术,智能电表和能源路由器又可作为信息聚合器实时采集其他节点的能源信息,并以数字货币(Digital currency)的形式完成电力交易.此外,能源路由器参与分布式电力优化调度.一方面,在产消者联盟内,能源路由器根据产消者节点的能源信息协调所有节点的对等电力交易.另一方面,在产消者联盟间,底层产消者节点可通过能源路由器与其他联盟内产消者节点达成对等电力交易.

能源金融系统网络层基于区块链技术向用户提供可编程智能合约和分布式交易记录等服务,其信任机制不是来源于第三方的背书,而是所有用户对于共识机制的认同.因此,用户可根据自身资源情况制定合约实现对等电力交易,并记录交易结果.

需要强调的是,虽然系统为分层架构,但三层网络间并不独立,而是通过能量流、信息流与业务流优化互动,且每一层的运行过程均受到人类意愿的制约.

1.3 IES系统的运行

结合平行系统理论[34],构建IES系统电力调度体系,如图2所示.纵向来看,系统可分为实际能源系统层、平行交互层和人工能源系统层.其中,人工能源系统通过平行交互的方式利用社会计算成果引导实际系统运行,从而激励产消者参与对等电力交易,实现电力资源的有效配置.横向来看,实际能源系统与人工能源系统均以社会福利最大为系统的控制目标输入,系统控制器由面向产消者节点的局域协调控制器和面向联盟路由器的广域协调控制器构成,二者组合互动,以提高用户满意度.

IES系统的电力调度依赖于事件驱动的信息交互体系.图3为[t,t+1]时段两个联盟构成的信息交互体系，图中编号代表的事件如下：

图3(a):t时刻,联盟内所有产消者向本地能源路由器发送供需差额数据;

图3(b):局域范围内,本地能源路由器根据差额数据以局域社会福利最大为目标协调对等电力交易,各产消者制定交易电量,并根据双向拍卖机制理性报价,同时向金融系统网络发送交易请求信息;

图3(c):金融系统认证交易信息后向能源路由器发送交易合约,并记录交易结果;

图3(d):本地能源路由器通过局域协调控制器执行合约规定的局域电力调度任务;

图3(e):局域协调后,部分仍存在供需差额的产消者再次向本地能源路由器反馈差额数据;

图1 IES系统多级对等协同优化架构Fig.1 The architecture of multilevel peer-to-peer co-optimization in IES

图3(f):广域范围内,各能源路由器根据反馈数据以广域社会福利最大为目标协调不同联盟间产消者对等电力交易,各产消者制定交易电量,并根据双向拍卖机制理性报价,同时向金融系统网络发送交易请求信息;

图3(g):金融系统认证交易信息后向各能源路由器发送交易合约,并记录交易结果;

图3(h):各能源路由器通过广域协调控制器执行合约规定的广域电力调度任务,并准备下一时段的信息交互.

其中,基于区块链技术的能源金融系统网络层通过数据加密、时间戳、共识机制等技术手段,在不依赖第三方可信背书的情况下实现用户友好的分布式电力交易[35-36].

2 协同优化策略

针对上述IES系统,综合考虑多级产消者对等互动过程中的社会福利、供求平衡和交易信息等因素,重构优化问题的目标函数,并采用智能算法进一步提高求解速度.同时,基于智能合约技术实现用户友好的电力交易模式.

2.1 社会福利问题

在IES系统中,产消者节点与联盟路由器节点实时执行由人工系统引导的电能交易以及电价定价.为量化用户对商品的需求效用,爱尔兰统计学家

图2 IES系统电力调度体系Fig.2 The power dispatching architecture of IES

图3 IES系统信息交互体系Fig.3 The information interaction architecture of IES

Geary在文献[37]中通过整合生活成本恒定效用指数,首次提出了Stone-Geary效用函数.在此基础上,文献[38]利用该效用函数模型描述电动汽车用户对车用电能的需求与其收入及电价的关系,更准确地反映了用户的消费偏好.由于电能作为商品可被用于电力交易,同时产消者节点作为理性人可根据自身需求产生理性购电意愿,产消者需求电量可以分为实际需求电量和基本需求电量,因此采用Stone-Geary效用函数来描述产消者对电能的需求效用.

当买方节点i(i=1,2,···,I)的发电量小于所需负荷,该节点需要向发电量丰富或负荷量较小的卖方节点j(j=1,2,···,J)请求电力交易,设实际需求量为(每个买方节点i的需求量由J个卖方节点提供),基本需求量为,每个买方节点的购电意愿为wi,则IES系统的Stone-Geary效用函数为

对式(1)的最后一项进行单调变换—取对数,可得到表征每一个买方节点购电意愿的效用函数.

对于卖方节点的效用函数,这里考虑用二次函数来表征.其中,交易双方的效用函数均为正向递增的可微严格凹函数,这里将社会福利问题定义为调度时间范围内最大化双方节点效用函数之差.另外,根据上述框架又将社会福利问题具体分为面向产消者节点的局域协调社会福利问题和面向联盟路由器节点的广域协调社会福利问题,以充分挖掘可再生能源的利用效率,提高系统全局经济效益.

2.1.1 局域协调(节点级)

在每一个产消者联盟内,定义产消者买方节点x=1,2,···,X,卖方节点y=1,2,···,Y,双方交易的电量分别为,则买方节点的效用函数表达如下:

对于产消者卖方节点,其成本函数表达如下:

其中,l1,l2为卖方节点的售电成本系数,且l1＞0.

由于联盟内产消者买方的意愿是使其效用函数最大,而产消者卖方的意愿是使其售电成本最小,从整个联盟的角度看,为充分满足双方意愿并实现联盟内有效的市场均衡[39],联盟需处理社会福利最大化问题以合理分配产消者的交易电量.因此,面向产消者节点的局域协调社会福利问题(Local coordination of social welfare problem,LCSWP)可以表示为

其中,αx,βx,γx,λxy,µxy为式 (4)的拉格朗日乘子,对应的集合分别为α,β,γ,λ,µ. 根据 Karush-Kuhn-Tucker(KKT)条件可知,LCSWP的最优解应满足以下条件:

对于LCSWP,人工系统有必要获取联盟内所有产消者节点效用函数和成本函数的完整信息,并利用式(8)和式(9)求解局域社会福利最大时的最优交易电量.

2.1.2 广域协调(路由级)

考虑到用户不同的负荷特征以及可再生能源有限的容量配置,仅通过局域协调仍很难实现联盟内供需平衡,因此产消者需要通过能源路由器与联盟外产消者进行广域协调,以进一步消除供需差额.定义参与广域协调的产消者买方节点m=1,2,···,M,卖方节点n=1,2,···,N,在广域协调的基础上,双方交易的电量分别为,则各节点的效用函数表达如下:

由此,面向联盟路由器的广域协调社会福利问题(Wide-area coordination of social welfare problem,WCSWP)可以表示为

由此,可利用式(13)和式(14)求解广域社会福利最大时的最优交易电量.

上述局域协调和广域协调的实质是实现电力资源的最优分配,通过人工系统量化产消者节点间和联盟路由器节点间的电力交易量,并引导实际系统完成交易电量的传输.然而,无论是买方节点还是卖方节点考虑的都是自身经济效益最大化的问题,仅根据人工系统提供的电量分配不足以满足节点的经济效益需求.为解决这一问题,还需要制定合理的定价机制以激励各节点理性报价并参与系统资源优化配置过程.

2.2 双向拍卖机制

在IES系统中,考虑产消者作为理性人,可根据实际交易电量理性报价,同时为提高各节点的经济效益,采用双向拍卖机制[40]制定电力成交价格.对于系统内所有节点,定义买方节点i=1,2,···,I,卖方节点j=1,2,···,J,此时买方节点实际购入电量为chaij,买方购价为buyij,卖方节点实际出售电量为disji,卖方售价为sellji.由于买方节点购电意愿越大,其购电消费预算越大,这里考虑用买方节点购价替代购电意愿,为获取更加符合实际消费情况的买方节点报价,买方节点的需求电量以实际购入电量为基准,则表征所有买方节点购电消费偏好的Stone-Geary效用函数为

对式(15)进行单调变换— 取对数,可得到表征同样偏好的效用函数.

由此,最优交易价格即成交价格可通过求解下面的资源优化配置问题(Optimal allocation problems,OAP)得到.

由于在优化定价的同时,需要确保局域协调和广域协调的社会福利最大,则式(17)的优化问题需要服从式(6)或式(12)的约束条件.构建拉格朗日乘子函数式L3,并根据Karush-Kuhn-Tucker(KKT)条件可知,OAP的最优解应满足以下条件:

其中,ξ=η或ξ=ρ表示产消者节点间或联盟路由器间电力传输效率.对比式(8),(9)和式(13),(14)可知

2.3 NR-PSO算法实现

由于传统基于梯度的求解方法难以确定合适的搜索步长,在最优值附近容易产生锯齿现象,此外,传统的智能优化方法在新的迭代周期通常采取随机初始化搜索或遍历搜索策略,导致算法的求解效率降低,故本文采用NR-PSO算法求解上述局域协调和广域协调组合互动模型.利用优化模型信息构造算法搜索过程的收敛判定域,通过收敛判定域的引导并结合全局随机寻优与区域定向寻优的策略,进一步提高算法局部搜索能力及收敛速度.

PSO算法作为一种并行进化算法,具有收敛速度快、全局搜索能力强等优势,在各领域中被广泛应用.粒子通过不断调整位置和速度更新公式来寻找全局最优解.然而,由于受迭代后期随机振荡现象的影响,在全局最优值附近的粒子搜索过程需较长时间才能最终收敛,故传统PSO算法存在局部搜索能力弱、易陷入局部最优等缺点.其进化方程组如下:

其中,ω为惯性权重,c1和c2为学习因子,r1和r2为区间(0,1)上的随机数,为第k代第u个粒子第d维个体历史最优位置,为第k代第d维粒子群体全局最优位置,为第k代第u个粒子第d维的速度和位置.

针对本文优化问题,定义空间维度为D,对于局域协调,D=X+Y.对于广域协调,D=M+N.解空间的宽度表示为,其中分别表示搜素空间的上界和下界粒子位置.对于买方节点,分别表示局域协调或广域协调下买方节点的可行购电量粒子范围.对于卖方节点,分别表示局域协调或广域协调下卖方节点的可行售电量粒子范围.同时,将粒子周围区域定义为收敛判定域,定义s为收敛因子,由此可获得收敛判定域宽度Wconv=sW.

此外,定义e为邻域因子,且e＞s,用Wneigh=eW表示邻域再调度的解空间宽度.第k代第u个粒子位置和全局最优粒子位置之间的距离可以表示为.迭代过程中,若distu＜Wconv,可按区域定向寻优的位置更新公式搜索.

在迭代前期,按全局随机寻优的位置更新公式搜索可加快收敛速度.而在迭代中后期,转为由收敛判定域引导的区域定向寻优过程以进一步提高算法局部搜索能力.由此,求解局域协调和广域协调组合互动模型的NR-PSO算法的具体步骤如下:

步骤1.初始化粒子种群参数,生成初始位置和初始速度;

步骤2.根据LCSWP的目标函数和约束求解,并判断是否满足终止条件,若满足,停止迭代,否则转到步骤3;

步骤3.当迭代k≤εkmax时,按全局随机寻优的式(22)和式(23)更新所有粒子的位置和速度,ε为迭代控制参数;

步骤4.当迭代k＞εkmax时,若粒子位置与全局最优粒子位置的距离小于收敛判定域宽度,即distu＜Wconv,则按式(24)更新粒子位置;

步骤5.生成新种群,返回步骤2;

步骤6.当k=kmax时,输出局域协调优化结果;

步骤7.结合局域协调优化结果判断是否需要进行广域协调,若不需要转到步骤8,否则根据WCSWP的目标函数和约束求解,重复步骤3～6,输出广域协调优化结果;

步骤8.根据两次寻优结果,代入式(20)和式(21),输出各节点最优交易电量及最优定价.

2.4 基于智能合约的对等电力交易

上述流程经两次求解局域协调和广域协调优化问题,分别得到各节点最优交易电量及最优定价信息.在IES多级对等协同优化架构下,对于系统内所有买方节点i和卖方节点j,基于智能合约的对等电力交易具体流程如下:

节点i预测其分布式可再生能源将在短时间[t0,t1]内波动,或电力负荷的突然增加需要更多的电力供应.经过上一节所述的局域协调和广域协调等社会计算后,将自动通过智能合约在能源金融系统网络层发布需求电量、需求时段以及交易报价等信息,且系统中所有节点都可以接收到需求信息.信息发布要求节点i处理自己的公钥iPublickey和私钥iPrivatekey,生成用于交易验证和信息传输的地址Ti-add=hash(iPublickey),其中,iPublickey=hash(iPrivatekey),Ti-add代表节点i的加密地址,是经过SHA256和Base58哈希计算之后的结果.发布的需求信息如下:

其中,chaij和buyij分别代表[t0,t1]时段节点i的需求电量和最优报价.

节点j接收到节点i的需求信息并计算可以出售的电量,即在[t0,t1]时段为节点i提供相应的电量disji,并通过生成相应的地址发送响应信息

式中,LTID是节点j在区块链中上一个交易信息的记录地址,Sign.j(chaij,[t0,t1],LT)是节点j在本次交易和上次交易内容的基础上通过jPublickey生成的加密签名.节点j接收到节点i的需求信息后,用自己的私钥iPrivatekey对j发送的消息摘要进行解密,得到Qi-add,j的全部内容.同时,通过jPublickey获取签名Sign.j中的详细信息.

由于存在多个电力供应节点,节点i确认交易的原则是最大限度地满足自身电力需求,因此,如果i接受j提供的电量,则i把智能合约脚本发送给j.如果不接受,则终止交易.

此外,为确保资源避免出现“二次支付”现象,记录交易内容的脚本信息必须包括区域电力网络运营商的签名Sign.DEO.因此,交易记录的脚本应包含交易双方的签名,区域电力网络运营商的签名以及交易内容摘要表示如下:

一旦Script由三方完成签名,则i需要按智能合约进行支付,j需要按智能合约向i传送约定量的电能.交易完成后,本次交易信息被广播到整个能源金融系统网络.在所有节点验证交易信息之后,本次交易信息将被写入当前区块.

3 实验与分析

本文以城市低压配电网中某小型住宅区域的IES系统为例验证所提优化模型和方法的有效性.该区域IES系统的光伏和风机装机容量分别为60kW 和40kW,参与IES系统服务的产消者规模为60户住宅.该区域的基本负荷水平、光伏、风机出力曲线如图4所示.为简化模型,考虑将产消者群体分为4组,每组中的可再生能源装机容量均匀分配.根据上述数据构建IES系统多级对等协同优化模型并进行仿真实验.T=24h,以5min为周期进行IES系统实时多级对等电力交易.实验基于MATLAB R2016a编程环境,计算主机为Intel Core i7 CPU,主频为3.4GHz,运行内存为16GB.

图4 光伏、风机出力及基本负荷水平曲线Fig.4 Renewable energy output and basic load curves

3.1 系统参数

在NR-PSO算法中,种群规模u=200,最大迭代次数kmax=20,收敛因子s=0.01,邻域因子e=0.03.基本购电需求,.其他相关参数如表1所示.

表1 优化模型相关参数Table 1 The values of related parameters

3.2 结果与讨论

3.2.1 算法性能分析

基于IES系统多级对等协同优化模型,为验证本文方法的有效性,分别采取NR-PSO算法、CQPSO[28]算法、AWCPSO[29]算法以及迭代下降算法[40]求解系统一个周期内的优化问题,算法的收敛特性如图5所示.

图5 算法性能的对比分析Fig.5 Performance comparison and analysis of different algorithms

结果表明,CQPSO算法和AWCPSO算法需要更多的迭代次数才能最终收敛.其中,CQPSO算法针对量子粒子群易陷入局部最优问题,利用混沌的遍历性,使粒子的搜索范围覆盖整个解空间,同时在迭代过程中多次引入扰动机制使粒子位置跳出局部最优.然而,该方法的遍历搜索和扰动过程大幅增加了迭代过程的计算负担,从而降低了求解效率.AWCPSO算法通过自适应减小权重和加速度系数,指导迭代过程中粒子飞行的速度和方向,可在一定程度上减小迭代后期随机振荡现象的影响.但在新的迭代周期随机初始化搜索,使得算法寻优过程依然十分缓慢,迭代下降算法平均在第12次迭代收敛.然而在全局最优值附近,每次迭代后目标函数值行进的距离越来越小,产生锯齿现象,从而降低算法在迭代后期的收敛速度.而通过收敛判定域引导并结合全局随机寻优与区域定向寻优策略的NR-PSO算法,平均在第8次迭代(ε=0.2)收敛至相应的社会福利最优值,相比迭代下降算法减少了33%的迭代次数,有效地提高了算法局部搜索能力及收敛速度.

此外,仿真还对比了NR-PSO算法迭代控制参数ε不同取值下的迭代次数和总计算时间,如表2所示.结果显示,迭代次数和总计算时间随ε取值的削减而逐渐减小.在ε取0.2时达到最佳的迭代次数(8次)和总计算时间(162.64s),进一步验证了NR-PSO算法的有效性.

表2 不同控制参数下的迭代次数和计算时间Table 2 The number of iterations and calculation time with different control parameters

考虑IES系统实时电力交易以及一段时间内各产消者节点的实时信息交互,由计算结果进一步可知,当算法迭代控制参数ε取值小于0.4时,算法总求解时间均不超过本文的优化调度周期5min.因此,该算法能够满足系统的实时性要求.

3.2.2 供求平衡分析

本文以4个不同联盟产消者的协同优化结果为例,具体分析电力供求平衡情况.通过供求平衡分析,不仅可以了解产消者间交易电量情况,而且可以了解IES系统多级对等协同优化模型的实时运行状况.

如图6～9所示,实际能量状态表示当前时段内负荷与可再生能源出力之间的偏差电量,正值表示需要购买的电量,负值表示需要出售的电量.局域协调交易电量表示该产消者与联盟内其他产消者交易电量情况.广域协调交易电量表示该产消者通过能源路由器与其他联盟产消者交易电量情况.

IES系统通过局域协调与广域协调的组合互动来保障实时电力供求平衡.在每一个调度周期开始时刻,产消者向本地能源路由器发送其供需差额数据.本地能源路由器根据该数据进行局域双向拍卖,并基于NR-PSO算法求解局域协调优化问题,以最大限度满足各产消者基本需求电量约束及局域社会福利最大化目标.然而部分产消者的实际需求电量未得到满足,则需再次向本地能源路由器反馈差额数据.此时局域范围内已无可调度电量,则本地能源路由器需请求与其他路由器进行广域双向拍卖,并再次调用优化算法求解广域协调优化问题,以保障本地产消者电力供求平衡及广域社会福利最大化.

然而,在某些时段局域协调和广域协调的组合互动仍然不能实现严格的供求平衡,例如图6中的21:00～23:00时段,图7中的13:00～15:00和23:00～24:00时段,图8中的1:00～3:00和23:00～24:00时段,以及图9中的17:00～19:00时段等.一方面,受可再生能源出力不确定性和用户负荷随机性的影响,产消者的偏差电量呈现出相应的波动性.另一方面,在上述时段系统内参与交易的一方需要处理多方请求信息,然而这些供应方或需求方给出的供应量或需求量的累加不能完全消除偏差电量,从而导致该时段的供求失衡.此时,无法消除的少量偏差电量可通过与主电网交易或通过调度储能设施的充放等方式来弥补.

图6 联盟1产消者3实时电力供求平衡Fig.6 Real-time supply and demand balance of 3rd prosumer in the first union

图7 联盟2产消者10实时电力供求平衡Fig.7 Real-time supply and demand balance of 10th prosumer in the second union

3.2.3 多级对等电力交易结果

以联盟1中产消者3节点在3:15～3:40时段内的交易结果为例,进一步分析IES系统中产消者多级对等电力交易的具体过程,如表3所示.

图8 联盟3产消者6实时电力供求平衡Fig.8 Real-time supply and demand balance of 6th prosumer in the third union

图9 联盟4产消者14实时电力供求平衡Fig.9 Real-time supply and demand balance of 14th prosumer in the fourth union

根据IES系统的信息交互流程,在产消者向金融系统网络发布交易信息时,买方根据当前的偏差电量制定报价策略,密封交易信息摘要并发布.卖方接收需求信息,并根据买方报价制定售电策略.需要指出的是,在3:15～3:20时段的偏差电量为0.457kWh,而从卖方购买的总电量为0.453kWh,这主要是由于当前时段没有适合的卖方或卖方不愿提供额外的0.004kWh电量.其他时段的交易信息及结果如表3所示.多级对等电力交易结果说明,基于智能合约的辅助服务可实现用户友好的电力交易模式,其分布式对等交易结构有效地解决了用户与中心机构之间的信任问题,同时维护了用户的隐私和电力交易的安全性.

表3 联盟1产消者3某时段多级对等电力交易结果Table 3 Multilevel peer-to-peer power trading results of the third prosumer node in the first union at a certain time

4 结论

针对能源电力系统的优化管理与控制问题,提出了信息物理融合的智慧能源系统多级对等协同优化方法.该方法基于Stone-Geary函数和双向拍卖机制重构了电力资源优化分配问题的目标函数,更好地反映了分布式电力交易中的产消者偏好.同时,利用收敛判定域引导的全局随机寻优与区域定向寻优策略,进一步提高了算法的收敛性.在保障系统社会福利最大化的前提下,有效地实现了产消者电力资源的优化分配,促进了电力资源的经济合理利用.此外,通过智能合约的辅助服务,进一步实现了用户友好的电力交易模式.