CN114678881A - 一种v2g辅助支撑下的地震灾后电网快速恢复方法 - Google Patents

一种v2g辅助支撑下的地震灾后电网快速恢复方法 Download PDF

Info

Publication number
CN114678881A
CN114678881A CN202210356915.0A CN202210356915A CN114678881A CN 114678881 A CN114678881 A CN 114678881A CN 202210356915 A CN202210356915 A CN 202210356915A CN 114678881 A CN114678881 A CN 114678881A
Authority
CN
China
Prior art keywords
power
disaster
distribution network
area
node
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210356915.0A
Other languages
English (en)
Other versions
CN114678881B (zh
Inventor
吴雨轩
杨祺铭
李雨濛
苏健
张过
向月
季陈林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sichuan University
Original Assignee
Sichuan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sichuan University filed Critical Sichuan University
Priority to CN202210356915.0A priority Critical patent/CN114678881B/zh
Publication of CN114678881A publication Critical patent/CN114678881A/zh
Application granted granted Critical
Publication of CN114678881B publication Critical patent/CN114678881B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/28Arrangements for balancing of the load in a network by storage of energy
    • H02J3/32Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
    • H02J3/322Arrangements for balancing of the load in a network by storage of energy using batteries with converting means the battery being on-board an electric or hybrid vehicle, e.g. vehicle to grid arrangements [V2G], power aggregation, use of the battery for network load balancing, coordinated or cooperative battery charging
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种V2G辅助支撑下的地震灾后电网快速恢复方法,包括如下步骤:建立配电网元件故障率模型,通过随机抽样生成地震灾害下的配电网故障场景;基于EV的出行特征与用户的出行需求,采用出行链理论拟合出行数据,在配电网故障场景的基础上搭建灾害前期EV调度模型;基于消费者心理学理论,将灾害前期EV调度模型纳入对EV用户灾后响应选择行为的考虑中,制定灾后反向输电激励响应机制,引导EV参与供电恢复,得到各V2G站出力;考虑配电网重构及V2G站反向充电,建立考虑二维锥约束的配网重构交流最优潮流模型;求解配网重构交流最优潮流模型,得到最优重构供电恢复方案。本发明提高地震灾后电网恢复力与供电恢复的经济性。

Description

一种V2G辅助支撑下的地震灾后电网快速恢复方法
技术领域
本发明属于电网领域,具体涉及一种V2G辅助支撑下的地震灾后电网快速恢复方法。
背景技术
现代社会的生产生活对“电”有着极高的依赖性,然而近年来全球自然灾害频发,其中地震给城市电网的安全稳定运行带来极大挑战,严重威胁到社会生产建设和人民的生命财产安全。
随着极端自然灾害逐渐增多,“弹性”概念的提出用于评价系统抵御扰动事件、适应及恢复的能力,对小概率-高损失事件具有较强恢复力稳定性的“弹性电网”受到越来越多关注。提升电网弹性意义重大,能显著减少国民经济生活和企业生产作业的停电损失。当前我国实际应用场景下的电网弹性提升方法中,一部分因系统自身电能资源有限且欠缺灵活性资源的应用造成弹性提升效果不佳,另一部分虽引入了移动电能资源,但过高的经济成本限制了其在弹性提升中的工程实际应用。
发明内容
本发明所要解决的技术问题是,提供一种V2G辅助支撑下的地震灾后电网快速恢复方法,生成灾后电网故障场景,基于该场景搭建灾害前期EV调度模型,将现有研究中受到较少关注的EV用户的自主选择行为及因此导致的时空分布特性纳入考虑范围;基于消费者心理学理论,制定灾后反向输电激励响应机制,引导V2G站内的EV参与供电恢复,考虑V2G可以提供一定量的快速响应资源,此举能够提高电网灾后恢复速度,同时充分利用了有限的电能资源,并且因减少移动电能资源的使用而降低了经济成本;建立考虑V2G反向输电和配网重构的经济性供电恢复模型并求解最优方案减少经济损失、提高响应速度和提升重要负荷恢复程度,提高了地震灾后电网的恢复力与供电恢复的经济性。
为解决上述技术问题,本发明采用的技术方案为:
一种V2G辅助支撑下的地震灾后电网快速恢复方法,包括如下步骤:
步骤1:建立配电网元件故障率模型,通过随机抽样生成地震灾害下的配电网故障场景;
步骤2:基于EV的出行特征与用户的出行需求,采用出行链理论拟合出行数据,在配电网故障场景的基础上搭建灾害前期EV调度模型;
步骤3:基于消费者心理学理论,将灾害前期EV调度模型纳入对EV用户灾后响应选择行为的考虑中,制定灾后反向输电激励响应机制,引导EV参与供电恢复,得到各V2G站出力;
步骤4:考虑配电网重构及V2G站反向充电,建立考虑二维锥约束的配网重构交流最优潮流模型;
步骤5:求解配网重构交流最优潮流模型,得到最优重构供电恢复方案。
进一步的,步骤1通过量化线路和杆塔受到地震的应力作用与故障率的关系以体现配电网元件受灾情况,采用非序贯蒙特卡洛模拟法生成配电网故障场景。
优先的,通过N-K原则筛选恶劣的故障场景,定义了系统故障场景中多重故障的阶数,K是固定值,根据决策需要进行调整,
Kmin≤∑si≤Kmax
其中,Kmin和Kmax分别为故障场景断线数的下限值和上限值,si是由非序贯蒙特卡洛模拟法采样生成的每条线路的状态,为0-1变量。
进一步的,步骤2包括如下步骤:
步骤2.1:建立基于出行需求的路网建模:
G=(V,E)
其中,G为路网集合,V为路网节点的集合,E为路网节点连通关系的集合;
步骤2.2:建立路网与电网的耦合模型;
步骤2.3:将整个路网构成的节点划分为城区U和郊区S,再将城区U分为三种功能区:居民区R、工作区W、休闲区L,则:
V={U,S}
U={R,W,L}
其中,U代表路网中属于城区的节点集合,S代表路网中属于郊区的节点集合,R代表路网城区中的居民区节点集合,W代表路网城区中的工作区节点集合,L代表路网城区中的休闲区节点集合,
路网划分的主要目的在于:筛选出V2G受灾程度低的区域作为目标区域,筛选的主要条件为区域的空旷程度,将城市划分为城区和郊区,郊区就把他作为空旷程度高的目标区域;城区的情况会复杂一些,所以细分为三个功能区,在功能区内进行以空旷程度为筛选条件进行目标区域与无关区域的划分;
步骤2.4:构建出行链中行程的状态变量集合,为:
TC={L0,Lf,P0f,T0,Tf,t0f,Ts}
其中,TC为出行链中一段行程的状态变量集合,包含起点L0、终点Lf、行驶路径P0f、出发时间T0、到达时间Tf、行驶时长T0f、停驻时长Ts
步骤2.5:对EV在目标区域和无关区域内的出行数据分别拟合,提取并判断灾害发生时每辆电动汽车的运行状态,搭建考虑用户出行行为的灾害前期EV 调度模型。
优选的,步骤2.3中出行链包括简单链和复杂链,简单链为2个节点,复杂链为3个节点。
进一步的,步骤3包括如下步骤:
步骤3.1:基于消费者心理学的激励-放电的不确定性建立响应过程的线性模型:
当经济激励水平为x时,用户m的V2G放电响应率上下界分别表示为:
Figure BDA0003583428980000031
Figure BDA0003583428980000032
其中,
Figure BDA0003583428980000033
Figure BDA0003583428980000034
分别为用户m的放电响应率上下界;
Figure BDA0003583428980000035
Figure BDA0003583428980000036
分别为不考虑和考虑不确定性情况下,放电响应率最小时的最低激励水平;
Figure BDA0003583428980000037
为用户m的最大V2G放电响应率;
步骤3.2:采用均匀分布描述某一经济激励水平下用户的V2G放电响应为:
Figure BDA0003583428980000041
其中,U表示均匀分布,
Figure BDA0003583428980000042
表示均匀分布描述下当经济激励水平为x时,用户m的V2G放电响应率;
步骤3.3:求解节点b处对应激励水平为xm的Nb个EV用户的实际V2G响应出力:
Figure BDA0003583428980000043
其中,
Figure BDA0003583428980000044
表示节点b处对应激励水平为xm的Nb个EV用户的实际V2G响应出力,
Figure BDA0003583428980000045
表示均匀分布描述下当经济激励水平为xm时,用户m的V2G放电响应率,
Figure BDA0003583428980000046
为经第一阶段到达节点b充电且剩余电量能够满足持续放电需求的EV数量,Pv2g为单辆电动汽车的反向充电功率;
步骤3.4:计算确定节点b处有意向参与配网供电恢复的EV数量:
Figure BDA0003583428980000047
其中,Nv2g表示节点b处有意向参与配网供电恢复的EV数量,
Figure BDA0003583428980000048
为与节点b耦合的V2G站额定放电功率,[]为取整符号。
进一步的,将灾后阶段划分为避难阶段和响应阶段,避难阶段是当灾害开始对电力系统造成破坏时,基于用户出行需求的EV自由调度停止,响应阶段是当灾害不再对电力系统造成破坏时,引导用户参与电力系统的灾后恢复过程,提出包含避难阶段和响应阶段的供电恢复引导补偿机制,响应阶段的供电恢复引导补偿机制参照基于消费者心理学的激励-放电的不确定性建立响应过程的线性模型;
在避难阶段的供电恢复引导机制,引入安全性驱动系数使灾前调度筛选出的目标区域与一般区域内可供避难的V2G站数量与之成正比例关系,
避难阶段实际响应并前往V2G站点避难的车辆数:
Nre=σ1Nroad+Nstation
其中,Nroad为未发生灾害时在路网行驶的EV数量;Nstation为在V2G站点停驻的EV数量;σ1为响应引导并前往V2G站点避难的EV用户比例;Nre为实际响应并前往V2G站点避难的车辆数;
其中,σ1受安全驱动系数影响:
σ1=kδ
其中,δ为安全驱动系数,k为比例因子,是大于0的常数,δ越大,σ1越大;
参与V2G反向输电的EV用户能否有足够的电量持续为电网提供支撑十分关键,因此需要对站内的EV进行筛选:
Figure BDA0003583428980000051
其中,Srem,t为能够参与供电恢复的EV剩余电量;
Figure BDA0003583428980000052
为V2G站点额定放电功率;TRescue为供电恢复抢修时间,停电到抢修恢复的时间由两部分组成。
进一步的,步骤4具体为:
将配网线路故障状态、各V2G站出力输入并更新配网系统拓扑结构,建立考虑二阶锥约束的配网重构交流最优潮流模型,制定V2G参与的灾后供电恢复方法,
其中,以最小化受灾配电网的切负荷损失为目标建立的考虑V2G反向输电和配网重构的经济性灾后恢复模型如下所示:
以考虑经济性的同时最小化受灾配电网的切负荷损失即失电量为目标函数:
Figure BDA0003583428980000053
i∈BU
bΔLS,i(Tbv-Tres)
其中,α为比例系数;ωi为节点i的重要负荷权重系数,系统中负荷按重要程度被划分为一级、二级和三级负荷;L0,i为供电恢复前的切负荷量,
Figure BDA0003583428980000054
为节点 i供电恢复后的切负荷量;bΔLS,i为节点i恢复单位负荷的综合等效收益,lLS,i为节点i单位时间内单位负荷因失电带来的经济损失,Tbv为供电恢复时长基准值,Tres为V2G参与供电恢复的响应时长,Qbasec为V2G引导机制补偿给用户的总成本;
考虑功率平衡、节点电压、支路电流和削负荷量等约束条件,并将它们松弛为二阶锥约束:
Figure BDA0003583428980000061
pij、qij分别为流过支路(i,j)的有功功率、无功功率,
Figure BDA0003583428980000062
是流过支路(i,j)的电流的平方,
Figure BDA0003583428980000063
是节点j上电压的平方。
本发明与现有技术相比,具有如下有益效果:
本发明所要解决的技术问题是,提供一种V2G辅助支撑下的地震灾后电网快速恢复方法,生成灾后电网故障场景,基于该场景搭建灾害前期EV调度模型,将现有研究中受到较少关注的EV用户的自主选择行为及因此导致的时空分布特性纳入考虑范围;基于消费者心理学理论,制定灾后反向输电激励响应机制,引导V2G站内的EV参与供电恢复,考虑V2G可以提供一定量的快速响应资源,此举能够提高电网灾后恢复速度,同时充分利用了有限的电能资源,并且因减少移动电能资源的使用而降低了经济成本;建立考虑V2G反向输电和配网重构的经济性供电恢复模型并求解最优方案减少经济损失、提高响应速度和提升重要负荷恢复程度,提高了地震灾后电网的恢复力与供电恢复的经济性。
1、本发明考虑灾害下EV出行特征及EV用户在灾害过程中的行为特征,建立了EV群体的灾害前期调度及灾后供电恢复响应模型,弥补了现有含EV的弹性电网研究中较少关注EV用户的自主选择行为及因此导致的时空分布特性的不足。
2、本发明将EV的灾前调度纳入对EV用户灾后响应选择行为的考虑中,提出包含避难阶段和响应阶段的供电恢复引导补偿机制,基于消费者心理学理论,制定灾后反向输电激励响应机制。本发明不考虑用户间竞争关系,选择剩余电能较多的EV参与供电恢复,提高电网灾后恢复速度,同时充分利用了有限的电能资源,并且因减少移动电能资源的使用而降低了经济成本。
3、本发明创造性的提出避难阶段的供电恢复引导机制,引入安全性驱动系数,使灾前调度筛选出的目标区域与一般区域内可供避难的V2G站数量与之成正比例关系。
附图说明
图1是本发明实施例中V2G辅助支撑下地震灾后电网快速恢复方法的流程图;
图2是本发明实施例中获得灾后供电恢复过程中V2G站出力的流程图;
图3是本发明实施例中各方案下失电量对比及500组配网故障场景下各方案失电量散点图;
图4是本发明实施例中各方案下等效成本对比及500组配网故障场景下各方案等效成本散点图;
图5是本发明实施例中各方案响应速度与重要负荷恢复程度对比;
图6是本发明实施例中包含重要负荷等级划分的改进的33节点配网系统示意图;
图7为本发明实施例中重构+V2G方案下一天内各小时等效成本和失电量变化曲线图;
图8为本发明实施例中确定断线场景下系统供电分布和范围示意图;
图9为本发明实施例中确定断线场景下EV群体聚集在居民/休闲区的系统供电分布和范围示意图。
具体实施方式
下面结合附图对本发明作更进一步地说明。
实际应用场景下的电网弹性提升方法的诸多灵活性资源中,电动汽车(EV,Electrical Vehicle)作为城市中广泛存在的分布式移动储能单元,能够以V2G (Vehicleto Grid,车辆到电网)的形式参与配网供电恢复的同时充分利用有限的电能资源并降低经济成本。考虑V2G可以提供一定量的快速响应资源,引入 V2G参与反向供电能够提高电网灾后恢复的速度。然而,现有研究较少关注EV 用户的自主选择行为及因此导致的时空分布特性。因此,本发明设计EV出行特征,考虑用户在灾害环境下选择行为的改变,研究灾后V2G辅助支撑的城市配电网弹性提升方法有着重要意义,未来EV应用范围必将不断扩大。
图1为本发明V2G辅助支撑下地震灾后电网快速恢复方法的流程图。如图 1所示,本发明V2G辅助支撑下地震灾后电网快速恢复方法,包括以下步骤:
S1:量化线路和杆塔受到地震的应力作用与故障率的关系以体现配电网元件受灾情况,采用非序贯蒙特卡洛模拟法生成配电网故障场景,并通过N-K原则筛选过于恶劣的故障场景。
具体的,采用导线和电线杆的故障率和受到的应力σg和MT之间的关系式来定量刻画。
Figure BDA0003583428980000081
Figure BDA0003583428980000082
其中,pfl为线路i的导线故障率,pfp为线路i的电线杆故障率;μ1和μp分别为导线和电线杆两元件强度的均值;δl和δp分别为导线和电线杆两元件强度的方差。式中导线和电线杆两元件强度的方差和均值都是由工程实际运行所得。
线路i正常工作的充要条件为线路上导线和电线杆都正常工作,因此,导线和电线杆可以理解为串联关系:
Figure BDA0003583428980000083
其中,pf,i为线路i的故障率,m1为线路i架设的电线杆数量。
接着,用非序贯蒙特卡洛模拟法采样配电网中每条线路的状态si:
Figure BDA0003583428980000084
其中,uf,i是服从(0,1)均匀分布的随机数;si为0-1变量,si=0表示线路i正常工作,si=1表示线路发生故障。
最后,根据N-K原则对生成的故障场景进行筛选,即定义了系统故障场景中多重故障的阶数,这一参数可以根据决策需要进行调整。
Kmin≤∑si≤Kmax
其中,Kmin和Kmax分别为故障场景断线数的下限值和上限值。
S2:考虑电动汽车(EV)的出行特征,对区域路网及电网进行建模,基于该模型下用户的出行需求采用出行链理论针对路网内无关区域与V2G站受灾程度较低的目标区域分别描述灾害前期EV用户的调度行为。
其中,灾害前期EV调度模型中基于出行需求的出行链建模如下所示:
G=(V,E)
将整个路网构成的区域划分为城区(U)和郊区(S),再将城区分为三种功能区:居民区(R)、工作区(W)、休闲区(L),并分别生成对应路网节点集合U,S,R,W 和L。可知,V={U,S}U={R,W,L}。
对路网与电网耦合建模,将两网耦合关系等效为路网、电网两层是否有层间节点连接。
ζ={Eαβ∈Vα×Vβ}
Figure BDA0003583428980000091
式中,ζ为路网和电网间节点的连边集合,Eαβ路网与电网构成的耦合边集合,Vα为路网节点集合,Vβ为电网节点集合,
Figure BDA0003583428980000092
为路网第i个节点,
Figure BDA0003583428980000093
为电网第 i个节点,若
Figure BDA0003583428980000094
则两网间路网节点i与电网节点j存在耦合关系。从能量传输上来看,两网耦合关系可以等效成路网上的站点从对应电网中的节点获取电能。
将行链分为简单链和复杂链两种,分别经过2和3个节点(包括终点):
TC={L0,Lf,P0f,T0,Tf,t0f,Ts}
其中TC为出行链中一段行程的状态变量集合,包含起点(不一定是出发地) L0、终点(不一定是目的地)Lf、行驶路径P0f、出发时间T0、到达时间Tf、行驶时长T0f、停驻时长Ts
对EV实时状态进行建模:
EVstatus={Srem,t,wt,Pt}
其中,EVstatus为在某时刻每辆EV实时状态信息的集合,包含该时刻每辆车的剩余电量Srem,t、行驶状态wt和离所在位置最近的V2G交互站编号Pt
S3:将EV的灾前调度纳入对EV用户灾后响应选择行为的考虑中,提出包含避难阶段和响应阶段的供电恢复引导补偿机制,基于消费者心理学理论,制定灾后反向输电激励响应机制,引导V2G站内的EV参与供电恢复,得到各V2G 站出力。
其中,基于消费者心理学的激励-放电的不确定性响应过程的线性模型如下所示:
当经济激励水平为x时,用户m的V2G放电响应率上下界分别表示为:
Figure BDA0003583428980000101
Figure BDA0003583428980000102
其中,
Figure BDA0003583428980000103
Figure BDA0003583428980000104
分别为用户m的放电响应率上下界;
Figure BDA0003583428980000105
Figure BDA0003583428980000106
分别为不考虑和考虑不确定性情况下,放电响应率最小时的最低激励水平;
Figure BDA0003583428980000111
为用户m的最大V2G放电响应率。
采用均匀分布描述某一经济激励水平下用户的V2G放电响应:
Figure BDA0003583428980000112
节点b处对应激励水平为xm的Nb个EV用户的实际V2G响应出力记为:
Figure BDA0003583428980000113
其中,
Figure BDA0003583428980000114
为经第一阶段到达节点b充电且剩余电量能够满足持续放电需求的EV数量,Pv2g为单辆电动汽车的反向充电功率。
由此可确定节点b处有意向参与配网供电恢复的EV数量:
Figure BDA0003583428980000115
其中,
Figure BDA0003583428980000116
为与节点b耦合的V2G站额定放电功率,[]为取整符号。本发明不考虑用户间竞争关系,选择剩余电能较多的EV参与供电恢复。
图2为本发明实施例中获得灾后供电恢复过程中V2G站出力的流程。如图2 所示,搭建灾害前期EV用户调度模型并制定灾后反向输电激励响应机制以得到 V2G站出力的方法为:
(1)基于图论对路网建模得出EV的最短路径行驶方案,将路网与电网的能量传输耦合关系等效成路网上的站点从对应电网中的节点获取电能;
(2)将城市划分为城区与郊区,再细分城区的功能区并分别生成对应路网节点集合,考虑对城区中三种功能区各分类讨论的结果与郊区,生成目标避难区域;
(3)采用出行链理论分别拟合该目标区域与路网中无关区域内每辆电动汽车一天内的出行数据,提取并判断灾害发生时每辆电动汽车的运行状态,搭建灾害前期EV调度模型。
(4)提出避难阶段的供电恢复引导机制,引入安全性驱动系数使灾前调度筛选出的目标区域与一般区域内可供避难的V2G站数量与之成正比例关系;
(5)在此阶段以安全性为主要导向原则进行舆论宣传引导并筛选电量足以持续为电网提供支撑的EV用户;
(6)提出响应阶段的供电恢复补偿机制,以经济利益为主要导向原则激励用户响应,基于消费者心理学的激励-放电响应采用线性模型描述放电的不确定性响应过程;
(7)采用均匀分布描述某一经济激励水平下用户的V2G放电响应,确定有意向参与配网供电恢复的EV数量和V2G服务聚合商补偿给用户的成本;
(8)计及经济激励水平在各区域间稳定与无关区域安全性驱动系数较低所导致的EV用户响应率低的情况,受灾严重的节点可引入应急发电车(MEG)参与供电恢复。
S4:将配网线路故障状态、各V2G站出力输入并更新配网系统拓扑结构,建立考虑二阶锥约束的配网重构交流最优潮流模型,制定V2G参与的灾后供电恢复方法。
其中,以最小化受灾配电网的切负荷损失为目标建立的考虑V2G反向输电和配网重构的经济性灾后恢复模型如下所示:
以考虑经济性的同时最小化受灾配电网的切负荷损失即失电量为目标函数:
Figure BDA0003583428980000121
i∈BU bΔLS,i(Tbv-Tres)
其中,α为比例系数;ωi为节点i的重要负荷权重系数,系统中负荷按重要程度被划分为一级、二级和三级负荷,系数分别为ωi=5,ωi=3和ωi=1;L0,i为供电恢复前的切负荷量,
Figure BDA0003583428980000134
为节点i供电恢复后的切负荷量;bΔLS,i为节点i恢复单位负荷的综合等效收益,lLS,i为节点i单位时间内单位负荷因失电带来的经济损失,Tbv为供电恢复时长基准值(取人工抢修从故障发生到故障修复所用平均时间),Tres为V2G参与供电恢复的响应时长,Qbasec为V2G引导机制补偿给用户的总成本。
考虑功率平衡、节点电压、支路电流和削负荷量等约束条件,并将它们松弛为二阶锥约束:
Figure BDA0003583428980000131
pij、qij分别为流过支路(i,j)的有功功率、无功功率,
Figure BDA0003583428980000132
是流过支路(i,j)的电流的平方,
Figure BDA0003583428980000133
是节点j上电压的平方。
S5:求解配电网供电恢复模型,得到地震灾后最优重构供电恢复方案。
本发明考虑地震与线路故障率关系,采用非序贯蒙特卡洛模拟法生成配电网断线场景;考虑灾害下EV出行特征及EV用户在灾害过程中的行为特征,建立了EV群体的灾害前期调度及灾后供电恢复响应模型,弥补了现有含EV的弹性电网研究中较少关注EV用户的自主选择行为及因此导致的时空分布特性的不足;考虑灾害故障场景下V2G反向出力求解最优供电恢复方法,该方法有利于提高配电网受灾后的恢复力与供电恢复的经济性。
为评估本方案的有益效果,我们选取了系统失电量、经济损失、响应速度与重要负荷恢复程度为评价指标。对所设计的三种方案进行仿真分析,仿真结果和数据分析如下:
1.系统失电量
如图3的表格所示,联络线重构对失电量的减少有一定帮助,能使灾后配网失电量降低28.26%;而引入V2G或MEG可以使失电量显著减少,配合联络线重构,可使灾后配网失电量分别降低60.18%和53.39%。
如图3所示,只重构方案的失电量分布较为分散,而重构+V2G与重构+MEG 的失电量分布大多集中在1MW以下,两者分布趋势差别不大。
由此可见,在考虑配电网联络线重构的基础上,引入V2G和引入MEG两种方案在失电量削减方面效果都较为显著,两者之间区别不大。
2.经济损失
如图4的表格所示,引入MEG的等效成本比只重构还高出43.70%,体现出 MEG参与供电恢复要付出极高的经济代价,符合如今主流弹性提升方案为尽快恢复重要负荷而不惜一切代价的情况;而引入V2G的等效成本比只重构降低了 14.68%,体现出V2G参与供电恢复有较高的经济性。
如图4所示,重构+MEG在所有场景下的等效成本都高于50000元,只重构与重构+V2G的大部分都低于50000元,但重构+V2G高于50000元的场景个数明显少于只重构的个数。
由此可见,在考虑配电网联络线重构的基础上,引入V2G的方案相比另外两方案在经济损失方面具有明显的优势。
3.响应速度与重要负荷恢复程度
如图5所示,重构和引入V2G相比引入MEG参与供电恢复响应明显更快。对于提高一级负荷恢复程度,运用联络线重构操作作用显著,只重构方案已然达到了很高的程度,引入V2G和MEG很难再有很大的提升。而对于提高二级负荷恢复程度,引入V2G相比只重构方案提高了21.98%,引入MEG相比只重构方案提高了10.70%。
由此可见,在响应速度方面,V2G与重构的方案相比引入MEG的方案有明显的优势;在重要负荷恢复程度方面,引入V2G配合联络线重构操作相比引入MEG 作用更显著。
综合上述多个评价指标的分析结果,能够验证该方法有利于提高配电网受灾后的经济性和恢复力。
采用改进的I EEE33节点配网系统和30节点路网系统为仿真场景。如图6所示,标出了V2G站接入节点的编号,负荷按照重要程度被分为一级、二级和三级负荷。
考虑EV时空特性对重构+V2G方案在灾后供电恢复过程中的影响,具体从改变灾害发生时间和选取典型的EV空间分布两方面进行探讨:
(1)灾害发生时间对重构+V2G方案在配网弹性提升中的影响
当EV群体数量恒定为300辆时,灾难发生在一天内各小时使用重构+V2G方案等效成本和失电量两指标变化如图7所示,可知重构+V2G方案受灾害发生时间影响明显,一天内不同时刻能提供的出力多少和经济性波动较大。
(2)EV空间分布对重构+V2G方案在配网弹性提升中的影响
选取确定断线场景下典型EV空间分布,研究EV聚集区域不同对供电恢复的影响。如图8所示为某确定断线场景下的系统供电分布和范围,如图9所示为在该断线场景下EV分别聚集在居民区和休闲区内时的系统供电恢复分布和范围,由此可知重构+V2G方案受EV空间分布影响明显,不同区域供电恢复效果差别较大,在EV聚集的区域系统各节点供电恢复效果较好。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (8)

1.一种V2G辅助支撑下的地震灾后电网快速恢复方法,其特征在于,包括如下步骤:
步骤1:建立配电网元件故障率模型,通过随机抽样生成地震灾害下的配电网故障场景;
步骤2:基于EV的出行特征与用户的出行需求,采用出行链理论拟合出行数据,在配电网故障场景的基础上搭建灾害前期EV调度模型;
步骤3:基于消费者心理学理论,将灾害前期EV调度模型纳入对EV用户灾后响应选择行为的考虑中,制定灾后反向输电激励响应机制,引导EV参与供电恢复,得到各V2G站出力;
步骤4:考虑配电网重构及V2G站反向充电,建立考虑二维锥约束的配网重构交流最优潮流模型;
步骤5:求解配网重构交流最优潮流模型,得到最优重构供电恢复方案。
2.根据权利要求1所述的一种V2G辅助支撑下的地震灾后电网快速恢复方法,其特征在于:步骤1通过量化线路和杆塔受到地震的应力作用与故障率的关系以体现配电网元件受灾情况,采用非序贯蒙特卡洛模拟法生成配电网故障场景。
3.根据权利要求2所述的一种V2G辅助支撑下的地震灾后电网快速恢复方法,其特征在于:通过N-K原则筛选恶劣的故障场景,定义了系统故障场景中多重故障的阶数,K是固定值,根据决策需要进行调整,
Kmin≤∑si≤Kmax
其中,Kmin和Kmax分别为故障场景断线数的下限值和上限值,si是由非序贯蒙特卡洛模拟法采样生成的每条线路的状态,为0-1变量。
4.根据权利要求1所述的一种V2G辅助支撑下的地震灾后电网快速恢复方法,其特征在于,步骤2包括如下步骤:
步骤2.1:建立基于出行需求的路网建模:
G=(V,E)
其中,G为路网集合,V为路网节点的集合,E为路网节点连通关系的集合;
步骤2.2:建立路网与电网的耦合模型;
步骤2.3:将整个路网构成的节点划分为城区U和郊区S,再将城区U分为三种功能区:居民区R、工作区W、休闲区L,则:
V={U,S}
U={R,W,L}
其中,U代表路网中属于城区的节点集合,S代表路网中属于郊区的节点集合,R代表路网城区中的居民区节点集合,W代表路网城区中的工作区节点集合,L代表路网城区中的休闲区节点集合,
路网划分的主要目的在于:筛选出V2G受灾程度低的区域作为目标区域,筛选的主要条件为区域的空旷程度,将城市划分为城区和郊区,郊区就把他作为空旷程度高的目标区域;城区的情况会复杂一些,所以细分为三个功能区,在功能区内进行以空旷程度为筛选条件进行目标区域与无关区域的划分;
步骤2.4:构建出行链中行程的状态变量集合,为:
TC={L0,Lf,P0f,T0,Tf,t0f,Ts}
其中,TC为出行链中一段行程的状态变量集合,包含起点L0、终点Lf、行驶路径P0f、出发时间T0、到达时间Tf、行驶时长T0f、停驻时长Ts
步骤2.5:对EV在目标区域和无关区域内的出行数据分别拟合,提取并判断灾害发生时每辆电动汽车的运行状态,搭建考虑用户出行行为的灾害前期EV调度模型。
5.根据权利要求4所述的一种V2G辅助支撑下的地震灾后电网快速恢复方法,其特征在于,步骤2.3中出行链包括简单链和复杂链,简单链为2个节点,复杂链为3个节点。
6.根据权利要求1所述的一种V2G辅助支撑下的地震灾后电网快速恢复方法,其特征在于,步骤3包括如下步骤:
步骤3.1:基于消费者心理学的激励-放电的不确定性建立响应过程的线性模型:
当经济激励水平为x时,用户m的V2G放电响应率上下界分别表示为:
Figure FDA0003583428970000031
Figure FDA0003583428970000032
其中,
Figure FDA0003583428970000033
Figure FDA0003583428970000034
分别为用户m的放电响应率上下界;
Figure FDA0003583428970000035
Figure FDA0003583428970000036
分别为不考虑和考虑不确定性情况下,放电响应率最小时的最低激励水平;
Figure FDA0003583428970000037
为用户m的最大V2G放电响应率;
步骤3.2:采用均匀分布描述某一经济激励水平下用户的V2G放电响应为:
Figure FDA0003583428970000038
其中,U表示均匀分布,
Figure FDA0003583428970000039
表示均匀分布描述下当经济激励水平为x时,用户m的V2G放电响应率;
步骤3.3:求解节点b处对应激励水平为xm的Nb个EV用户的实际V2G响应出力:
Figure FDA00035834289700000310
其中,
Figure FDA00035834289700000311
表示节点b处对应激励水平为xm的Nb个EV用户的实际V2G响应出力,
Figure FDA00035834289700000312
表示均匀分布描述下当经济激励水平为xm时,用户m的V2G放电响应率,
Figure FDA00035834289700000313
为经第一阶段到达节点b充电且剩余电量能够满足持续放电需求的EV数量,Pv2g为单辆电动汽车的反向充电功率;
步骤3.4:计算确定节点b处有意向参与配网供电恢复的EV数量:
Figure FDA00035834289700000314
其中,Nv2g表示节点b处有意向参与配网供电恢复的EV数量,
Figure FDA00035834289700000315
为与节点b耦合的V2G站额定放电功率,[ ]为取整符号。
7.根据权利要求6所述的一种V2G辅助支撑下的地震灾后电网快速恢复方法,其特征在于,将灾后阶段划分为避难阶段和响应阶段,避难阶段是当灾害开始对电力系统造成破坏时,基于用户出行需求的EV自由调度停止,响应阶段是当灾害不再对电力系统造成破坏时,引导用户参与电力系统的灾后恢复过程,提出包含避难阶段和响应阶段的供电恢复引导补偿机制,响应阶段的供电恢复引导补偿机制参照基于消费者心理学的激励-放电的不确定性建立响应过程的线性模型;
在避难阶段的供电恢复引导机制,引入安全性驱动系数使灾前调度筛选出的目标区域与一般区域内可供避难的V2G站数量与之成正比例关系,
避难阶段实际响应并前往V2G站点避难的车辆数:
Nre=σ1Nroad+Nstation
其中,Nroad为未发生灾害时在路网行驶的EV数量;Nstation为在V2G站点停驻的EV数量;σ1为响应引导并前往V2G站点避难的EV用户比例;Nre为实际响应并前往V2G站点避难的车辆数;
其中,σ1受安全驱动系数影响:
σ1=kδ
其中,δ为安全驱动系数,k为比例因子,是大于0的常数,δ越大,σ1越大;
参与V2G反向输电的EV用户能否有足够的电量持续为电网提供支撑十分关键,因此需要对站内的EV进行筛选:
Figure FDA0003583428970000041
其中,Srem,t为能够参与供电恢复的EV剩余电量;
Figure FDA0003583428970000042
为V2G站点额定放电功率;TRescue为供电恢复抢修时间,停电到抢修恢复的时间由两部分组成。
8.根据权利要求7所述的一种V2G辅助支撑下的地震灾后电网快速恢复方法,其特征在于,步骤4具体为:
将配网线路故障状态、各V2G站出力输入并更新配网系统拓扑结构,建立考虑二阶锥约束的配网重构交流最优潮流模型,制定V2G参与的灾后供电恢复方法,
其中,以最小化受灾配电网的切负荷损失为目标建立的考虑V2G反向输电和配网重构的经济性灾后恢复模型如下所示:
以考虑经济性的同时最小化受灾配电网的切负荷损失即失电量为目标函数:
Figure FDA0003583428970000051
i∈BU
bΔLS,i(Tbv-Tres)
其中,α为比例系数;ωi为节点i的重要负荷权重系数,系统中负荷按重要程度被划分为一级、二级和三级负荷;L0,i为供电恢复前的切负荷量,
Figure FDA0003583428970000052
为节点i供电恢复后的切负荷量;bΔLS,i为节点i恢复单位负荷的综合等效收益,lLS,i为节点i单位时间内单位负荷因失电带来的经济损失,Tbv为供电恢复时长基准值,Tres为V2G参与供电恢复的响应时长,Qbasec为V2G引导机制补偿给用户的总成本;
考虑功率平衡、节点电压、支路电流和削负荷量等约束条件,并将它们松弛为二阶锥约束:
Figure FDA0003583428970000053
pij、qij分别为流过支路(i,j)的有功功率、无功功率,
Figure FDA0003583428970000054
是流过支路(i,j)的电流的平方,
Figure FDA0003583428970000055
是节点j上电压的平方。
CN202210356915.0A 2022-04-06 2022-04-06 一种v2g辅助支撑下的地震灾后电网快速恢复方法 Active CN114678881B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210356915.0A CN114678881B (zh) 2022-04-06 2022-04-06 一种v2g辅助支撑下的地震灾后电网快速恢复方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210356915.0A CN114678881B (zh) 2022-04-06 2022-04-06 一种v2g辅助支撑下的地震灾后电网快速恢复方法

Publications (2)

Publication Number Publication Date
CN114678881A true CN114678881A (zh) 2022-06-28
CN114678881B CN114678881B (zh) 2023-03-07

Family

ID=82079026

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210356915.0A Active CN114678881B (zh) 2022-04-06 2022-04-06 一种v2g辅助支撑下的地震灾后电网快速恢复方法

Country Status (1)

Country Link
CN (1) CN114678881B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115809836A (zh) * 2023-02-09 2023-03-17 华南理工大学 考虑分布式储能应急供电能力的配电网韧性规划的方法
CN116031918A (zh) * 2023-01-09 2023-04-28 国网湖南省电力有限公司 台风灾害下的城市电力系统动态弹性供电恢复方法及系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111539566A (zh) * 2020-04-21 2020-08-14 燕山大学 一种考虑灾前预调度的配电网多故障抢修恢复方法及系统
CN111555280A (zh) * 2020-05-29 2020-08-18 山东大学 基于电-气综合能源系统的弹性配电网灾后恢复控制方法
CN111799785A (zh) * 2020-07-03 2020-10-20 国网新疆电力有限公司电力科学研究院 极端灾害后配电网电力与通信协调恢复方法及系统
CN112398122A (zh) * 2020-11-02 2021-02-23 清华大学 一种极寒灾害下热电耦合系统应急恢复方法
CN112884245A (zh) * 2021-03-22 2021-06-01 国网电力科学研究院有限公司 两阶段配电网灾后抢修调度及负荷恢复协同优化方法及系统
CN113452051A (zh) * 2021-06-25 2021-09-28 燕山大学 考虑应急电源车调度的有源配电网故障均衡供电恢复方法
CN113705964A (zh) * 2021-07-16 2021-11-26 国网上海市电力公司 一种配电网韧性恢复灾前预案制定方法及装置
CN114123167A (zh) * 2021-10-27 2022-03-01 西安交通大学 基于移动应急通信的配电网信息物理协同恢复方法及系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111539566A (zh) * 2020-04-21 2020-08-14 燕山大学 一种考虑灾前预调度的配电网多故障抢修恢复方法及系统
CN111555280A (zh) * 2020-05-29 2020-08-18 山东大学 基于电-气综合能源系统的弹性配电网灾后恢复控制方法
CN111799785A (zh) * 2020-07-03 2020-10-20 国网新疆电力有限公司电力科学研究院 极端灾害后配电网电力与通信协调恢复方法及系统
CN112398122A (zh) * 2020-11-02 2021-02-23 清华大学 一种极寒灾害下热电耦合系统应急恢复方法
CN112884245A (zh) * 2021-03-22 2021-06-01 国网电力科学研究院有限公司 两阶段配电网灾后抢修调度及负荷恢复协同优化方法及系统
CN113452051A (zh) * 2021-06-25 2021-09-28 燕山大学 考虑应急电源车调度的有源配电网故障均衡供电恢复方法
CN113705964A (zh) * 2021-07-16 2021-11-26 国网上海市电力公司 一种配电网韧性恢复灾前预案制定方法及装置
CN114123167A (zh) * 2021-10-27 2022-03-01 西安交通大学 基于移动应急通信的配电网信息物理协同恢复方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HANG SHUAI等: "Post-Storm Vehicle Routing for Distribution Grid Restoration: An OLUCT Based Learning Approach", 《2020 IEEE POWER & ENERGY SOCIETY GENERAL MEETING (PESGM)》 *
YIN XU等: "Resilience-Oriented Distribution System Restoration Considering Mobile Emergency Resource Dispatch in Transportation SystemResilience-Oriented Distribution System Restoration", 《IEEE ACCESS》 *
龙雪梅等: "考虑路网-电网交互和用户心理的电动汽车充电负荷预测", 《电力系统自动化》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116031918A (zh) * 2023-01-09 2023-04-28 国网湖南省电力有限公司 台风灾害下的城市电力系统动态弹性供电恢复方法及系统
CN116031918B (zh) * 2023-01-09 2024-04-30 国网湖南省电力有限公司 台风灾害下的城市电力系统动态弹性供电恢复方法及系统
CN115809836A (zh) * 2023-02-09 2023-03-17 华南理工大学 考虑分布式储能应急供电能力的配电网韧性规划的方法
CN115809836B (zh) * 2023-02-09 2023-05-23 华南理工大学 考虑分布式储能应急供电能力的配电网韧性规划的方法

Also Published As

Publication number Publication date
CN114678881B (zh) 2023-03-07

Similar Documents

Publication Publication Date Title
CN114678881B (zh) 一种v2g辅助支撑下的地震灾后电网快速恢复方法
AU2020294331A1 (en) Systems and methods for grid operating systems in electric power systems
CN108631306B (zh) 一种电力系统灾后恢复能力的评估方法
CN108493998B (zh) 考虑需求响应与n-1预想故障的鲁棒输电网规划方法
CN107292481B (zh) 一种基于节点重要度的电网关键节点评估方法
CN113569411B (zh) 一种面向灾害天气的电网运行风险态势感知方法
CN110212527B (zh) 一种配电网一次网架结构与配电自动化协同规划方法
CN109598377A (zh) 一种基于故障约束的交直流混合配电网鲁棒规划方法
CN107301479B (zh) 基于自然灾害风险的输电系统多场景规划方法
CN116187165A (zh) 基于改进粒子群算法的电网弹性提升方法
CN113239512A (zh) 一种考虑韧性的交直流配电网规划方案筛选方法及系统
CN110970891A (zh) 极端天气下基于多能协调的配电网弹性提升方法
CN113469750A (zh) 一种考虑极端天气的充电站与配电网协调规划方法及系统
CN113328437A (zh) 一种智能配电网cps拓扑构建方法及故障恢复方法
CN111382912B (zh) 一种交通网层能量最优分配策略确定方法及系统
CN115809836B (zh) 考虑分布式储能应急供电能力的配电网韧性规划的方法
CN109002938B (zh) 考虑n-1安全准则的交直流混合配电网双层规划方法
CN115495989A (zh) 一种考虑风特性的海上风电场集电系统拓扑结构设计方法
CN115864382A (zh) 一种高山山区城市电网渐进式优化规划设计方法
Wu et al. A Rapid Recovery Strategy for Post-Earthquake Power Grids with the Auxiliary Support of V2G
CN115459256A (zh) 基于配电网-交通网协同的配网负荷恢复方法及系统
CN111027855A (zh) 考虑输电线路气象灾害故障概率的电力系统风险控制方法
CN111092430A (zh) 一种适应于电力系统恢复的应急资源优化配置方法
CN111582585B (zh) 一种充电站与无线充电站的联合规划方法及系统
CN117220285B (zh) 配电网资源配置方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant