CN113452051A - 考虑应急电源车调度的有源配电网故障均衡供电恢复方法 - Google Patents

考虑应急电源车调度的有源配电网故障均衡供电恢复方法 Download PDF

Info

Publication number
CN113452051A
CN113452051A CN202110712635.4A CN202110712635A CN113452051A CN 113452051 A CN113452051 A CN 113452051A CN 202110712635 A CN202110712635 A CN 202110712635A CN 113452051 A CN113452051 A CN 113452051A
Authority
CN
China
Prior art keywords
power supply
load
eps
distribution network
power
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
CN202110712635.4A
Other languages
English (en)
Other versions
CN113452051B (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.)
Yanshan University
Original Assignee
Yanshan 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 Yanshan University filed Critical Yanshan University
Priority to CN202110712635.4A priority Critical patent/CN113452051B/zh
Publication of CN113452051A publication Critical patent/CN113452051A/zh
Application granted granted Critical
Publication of CN113452051B publication Critical patent/CN113452051B/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
    • 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/007Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
    • H02J3/0073Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source when the main path fails, e.g. transformers, busbars
    • 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
    • H02J3/388Islanding, i.e. disconnection of local power supply from the network
    • 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
    • H02J3/46Controlling of the sharing of output between the 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
    • H02J7/00Circuit arrangements for charging or depolarising batteries or for supplying loads from batteries
    • H02J7/0047Circuit arrangements for charging or depolarising batteries or for supplying loads from batteries with monitoring or indicating devices or circuits
    • H02J7/0048Detection of remaining charge capacity or state of charge [SOC]
    • 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]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

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

Abstract

本发明提供一种考虑应急电源车调度的有源配电网故障均衡供电恢复方法,具体步骤如下:根据有源配电网特征,分别对配电网的负荷特性、储能系统的充放电特性、应急电源车的调度特性和分布式电源的均衡供电特性进行建模;建立考虑应急电源调度的动态孤岛均衡恢复双层优化模型;基于双层优化模型和Dijkstra算法,建立交通路网模型支撑下的应急电源车优化调度方法,得到其最佳候选接入节点和最优调度路径;利用Big‑M方法对双层数学模型进行线性化处理,并利用Cplex得到全局最优故障恢复方案。本发明基于双层优化模型和dijkstra算法解决了EPS优化调度问题,减少了因调度产生的负荷中断时间,同时保证了停电周期内同等权重的负荷公平、均衡的恢复,提高了供电可靠性。

Description

考虑应急电源车调度的有源配电网故障均衡供电恢复方法
技术领域
本发明涉及配电网领域,特别涉及一种考虑应急电源车调度的有源配电网故障均衡供电恢复方法。
背景技术
近年来,化石能源的不断开发利用带来了严重的环境污染,导致全球极端天气和灾害事件频发,进而增大了配电网发生大面积停电事故的风险。因此研究极端灾害事件后配电网故障的快速恢复具有重要意义;同时,随着风电光伏等新能源接入比例的提高,我国已跃居成为风光发电装机容量榜首,大力发展新能源发电产业对保证我国电力供应和早日实现“碳达峰”和“碳中和”目标具意义重大。因此,当配电网发生大面积停电事故后,研究如何利用现存的光伏(photovoltaic,PV)、风机(Wind Turbine,WT)、储能(Energy StorageSystems,ESS)、可移动电源车(Mobile Power Vehicle,MPV)等分布式能源(DistributedEnergy Resource,DER)实现有源配电网的快速恢复是亟需解决的问题。
当前研究虽考虑了灵活性资源支撑配电网的恢复,但在制定有源配电网的恢复策略时,并未在追求负荷恢复价值最大化的同时优化恢复时间,导致抢修工作完成前同等权重的关键负荷可持续恢复时间差异显著,有失公平性和均衡性,降低了关键用户满意度;同时,传统的孤岛恢复仅利用配网内接入的风电、光伏和储能等有限分布式能源不能保证关键负荷的全部恢复,应急电源车作为一种大容量的可移动电源,可在补充和支撑配电网故障恢复和辅助电网运行方面发挥重要作用,但其优化调度问题未合理解决。
因此,有必要发明一种考虑应急电源车调度的有源配电网故障均衡供电恢复方法,使得配电网发生大面积、长时间停电事故后,可快速制定配电网内应急资源的调度策略、协调好多种发电资源的供电范围与关键负荷供电时长之间的关系,提高配电网的均衡性、供电可靠性和经济性。
发明内容
针对现有技术存在的问题,本发明提供一种考虑应急电源车调度的有源配电网故障均衡供电恢复方法,主要是通过建立动态孤岛均衡恢复双层模型和dijkstra算法解决配电网均衡恢复和应急电源车调度问题,有效提高了故障阶段配电网的供电可靠性、均衡性和经济性。
本发明提供了一种考虑应急电源车调度的有源配电网故障均衡供电恢复方法,具体实施步骤如下:
S1、根据有源配电网特征,分别对配电网的负荷特性、储能系统的充放电特性、应急电源车(Emergency Power Supply Vehicle,EPS)的调度特性和分布式电源(Distributed Energy Resource,DER)的均衡供电特性进行建模:
S11、根据负荷时变性和负荷等级建立配电网的负荷特性模型,所述负荷时变性是根据各时段的用电需求,对节点日负荷时变曲线进行积分,具体表达式为:
Figure BDA0003134296620000021
式中,Lk(a)为节点k在a到a+Δa时段的负荷功率;fk(t)为节点k的负荷预测曲线;
同时,依据负荷停供对用户经济损失、人身安全造成的影响将负荷等级分为一级、二级和三级;
S12、通过储能(Energy Storage Systems,ESS)调度策略平抑风光机组出力波动性和负荷不确定性,建立储能系统的充放电特性模型,根据充放电特性模型的公式求得各储能SOC需满足相关状态约束和充放电功率约束,具体表达式为:
Figure BDA0003134296620000022
Figure BDA0003134296620000023
Figure BDA0003134296620000024
Figure BDA0003134296620000025
式中,
Figure BDA0003134296620000026
分别为ESSk中SOC最小值和最大值;
Figure BDA0003134296620000027
分别为ESSk单位时段内最大充放电功率;
Figure BDA0003134296620000031
为二元变量,分别表示充放电状态,
Figure BDA0003134296620000032
为1表示充电,否则为0;
Figure BDA0003134296620000033
为1表示放电,否则为0;假设ESS有足够的无功补偿容量;
S13、分别根据应急电源车的充放电管理约束和应急电源车的候选接入节点配置约束建立应急电源车的调度特性模型;
S14、分别从平均恢复时间指标、第k个风光储的最大持续供电时间
Figure BDA0003134296620000034
和第k辆EPS的最大持续供电时间
Figure BDA0003134296620000035
建立分布式能源的均衡供电模型;
S2、建立考虑应急电源调度的动态孤岛均衡恢复双层优化模型:
S21、以最大化负荷恢复水平和配电网运行效能建立下层模型的目标函数,具体表达式为:
Figure BDA0003134296620000036
式中,wi为第i个负荷的权重因子;ci为单位时间内第i个负荷缺供造成的损失费用;cEPS为EPS为负荷供电时单位容量消耗费用;sEPS为EPS调度过程中单位时间产生的车程费用;
Figure BDA0003134296620000037
为EPS候选接入点变量;
Figure BDA0003134296620000038
为t时段第k辆EPS从候选接入点a调运至候选接入点b的行驶时间;Q为配电网内可接入应急电源车供电的候选节点集合;P为预放车库在内的候选节点集合,
Figure BDA0003134296620000039
b∈Q,a∈P;
Figure BDA00031342966200000310
为在t时段内EPSk为第i个负荷提供的容量;
Figure BDA00031342966200000311
为EPS应急供电决策变量;si,t=1为t时段第i个负荷由风光储或EPS恢复;
S22、根据配电网孤岛恢复均衡性最优建立上层模型的目标函数,具体表达式为:
Figure BDA00031342966200000312
式中,
Figure BDA00031342966200000313
表示第k个风光储的最大持续供电时间;
Figure BDA00031342966200000314
表示第k辆EPS的最大持续供电时间;Tnet表示配电网平均恢复时间指标;
S23、从配电网负荷的节点电压、支路功率、孤岛潮流平衡、配电网辐射运行、孤岛划分和失电负荷恢复总量一致六个因素建立动态孤岛均衡恢复双层优化模型的约束条件;
S24、根据步骤S21、S22和S23得到的目标函数和约束函数建立动态孤岛均衡恢复双层优化模型,具体表达式为:
Figure BDA0003134296620000041
式中,si,t表示t时段第i个负荷是否被恢复的二元决策变量;
Figure BDA0003134296620000042
分别为t时段第i个负荷是否被风光储和EPS恢复的二元决策变量;
Figure BDA0003134296620000043
为EPS候选接入点决策变量;
S3、基于步骤S2中双层优化模型和迪杰斯特拉算法,建立交通路网模型支撑下的应急电源车优化调度方法,以得到其最佳候选接入节点和最优调度路径;
S4、利用Big-M方法对步骤S24中双层优化模型进行线性化处理,将模型转换为混合整数线性规划,并利用Cplex求解得到全局最优故障恢复方案:
S41、利用Big-M方法等效线性化处理步骤S14的表达式分别为:
Figure BDA0003134296620000044
Figure BDA0003134296620000045
Figure BDA0003134296620000046
式中,x1~x6为6个0-1变量;M为无穷大常数;si,t表示t时段第i个负荷是否被恢复的二元决策变量;
Figure BDA0003134296620000051
分别为t时段第i个负荷是否被风光储和EPS恢复的二元决策变量;
Figure BDA0003134296620000052
为EPS候选接入点决策变量;Pi,t为第i个负荷在t时段内功率需求;
Figure BDA0003134296620000053
为DER在t时段内的状态;
Figure BDA0003134296620000054
为EPSk在时段t内的放电功率;
Figure BDA0003134296620000055
为第k个风光储系统在t时段内的最大出力;
S42、将辅助二进制变量
Figure BDA0003134296620000056
替代二进制变量
Figure BDA0003134296620000057
并增加等效逻辑约束,对步骤S21表达式中x×y类型的非线性约束实现线性等效转换,所述等效逻辑约束的表达式如下:
Figure BDA0003134296620000058
S43、对步骤S24建立的动态孤岛均衡恢复双层优化模型中的相关非线性约束进行替换,并利用CPLEX求解。
可优选的是,在步骤S12中,所述储能系统的充放电特性模型用充电状态SOC表示,其具体表达式为:
Figure BDA0003134296620000059
式中,SOCk,t为ESSk在t时段的SOC状态;
Figure BDA00031342966200000510
分别为ESSk在t-1时段内的充放电功率;ηch、ηdch分别为ESS的充放电效率;Ek为ESSk的容量;RES为DER中ESS的集合;Δt=1h。
可优选的是,在步骤S13中,所述应急电源车的充放电管理约束的具体表达式为:
Figure BDA00031342966200000511
Figure BDA00031342966200000512
Figure BDA00031342966200000513
Figure BDA00031342966200000514
Figure BDA00031342966200000515
式中,
Figure BDA00031342966200000516
分别为EPSk在时段t内的充放电功率;ηc,EPS、ηd,EPS分别为EPS的充放电效率;
Figure BDA00031342966200000517
分别为EPSk单位时段内最大充放电功率;
Figure BDA00031342966200000518
为EPSk的容量;
Figure BDA00031342966200000519
为二元变量,分别表示EPS的充放电状态;
Figure BDA00031342966200000520
Figure BDA00031342966200000521
分别表示EPS在故障时刻的初始荷电状态、故障结束时刻的荷电状态;REPS为DER中EPS的集合;
所述应急电源车的候选接入节点配置约束的具体表达式为:
Figure BDA0003134296620000061
Figure BDA0003134296620000062
Figure BDA0003134296620000063
Figure BDA0003134296620000064
式中,
Figure BDA0003134296620000065
为EPS候选接入点变量;vbk,t为节点动态孤岛划分变量;Q为配电网内可接入应急电源车供电的候选节点集合,P为包括预放车库在内的候选节点集合,
Figure BDA0003134296620000066
b∈Q,a∈P;R为分布式电源数量集合,R={1,2,....,m},k∈R。
可优选的是,在步骤S14中,所述平均恢复时间指标由Tnet表示,其具体表达式为:
Figure BDA0003134296620000067
Figure BDA0003134296620000068
式中:
Figure BDA0003134296620000069
为DER在t时段内的状态;si,t表示负荷i的连接状态;
Figure BDA00031342966200000610
为第k个风光储系统在t时段内的最大出力;Pi,t为第i个负荷在t时段内功率需求;m为配电网中DER的数量;n为配电网中失电负荷量;I表示失电负荷节点集合,I={1,2,....,n},i∈I;R表示分布式电源数量集合,R={1,2,....,m},k∈R;t0为停电起始时刻;
将第k个风光储系统作为电源,以最大化负荷恢复价值为目标函数进行孤岛划分,并计算第k个DER能够持续为该孤岛内负荷供电的时间,得到第k个风光储的最大持续供电时间
Figure BDA00031342966200000611
的具体表达式为:
Figure BDA00031342966200000612
Figure BDA00031342966200000613
式中,
Figure BDA00031342966200000614
为第k个风光储系统在t时段内的状态;
Figure BDA00031342966200000615
为节点i分配变量,若负荷i在t时段内由第k个风光储系统恢复,则
Figure BDA00031342966200000616
否则
Figure BDA00031342966200000617
根据EPS的调度得到最佳候选接入节点,并形成最优孤岛为负荷供电,得到第k辆EPS的最大持续供电时间
Figure BDA0003134296620000071
的具体表达式为:
Figure BDA0003134296620000072
Figure BDA0003134296620000073
式中,
Figure BDA0003134296620000074
为第k个EPS在t时段内的状态;
Figure BDA0003134296620000075
为应急电源车最优匹配决策变量。
可优选的是,所述步骤S23具体包括以下步骤:
S231、节点电压的约束条件表达式如下:
Ui min≤Ui≤Ui max,i=1,2,…,m
式中,Ui min和Ui max分别为节点i电压上下限;
S232、支路功率的约束条件表达式如下:
Figure BDA0003134296620000076
式中,pg max为支路g能够允许的最大功率;
S233、根据线性化的DistFlow公式,建立含DER的有源配电网潮流方程,得到孤岛潮流平衡的约束条件表达式如下:
Figure BDA0003134296620000077
Figure BDA0003134296620000078
Figure BDA0003134296620000079
Figure BDA00031342966200000710
式中,βi,j,t表示支路状态;
Figure BDA00031342966200000711
时段t内支路g的有功功率;
Figure BDA00031342966200000712
时段t内支路g的无功功率;ri支路i的电阻;xi为支路i的电抗;V1为参考电压,取V1=1;σ为电压偏差,取σ=5%;g(·,i)为功率流入母线i的支路集合;g(i,·)为功率流出母线i的支路集合;Ig为支路g直接相连的母线集合;
S234、配电网辐射运行的约束条件表达式如下:
G∈g
式中,G为重构后的配电网网络拓扑,g为辐射状网络拓扑集合;
S235、将由节点j处的风光储或EPS供电的子网称为孤岛,孤岛划分的约束条件表达式如下:
Figure BDA0003134296620000081
Figure BDA0003134296620000082
式中,zj表示节点j是否有电源,zj=1表示节点j与风光储系统或EPS连接;Q为EPS连接的候选接入节点集合;F为风光储连接的根节点集合;
S236、根据任意时段划分入孤岛恢复范围内的负荷都需由风光储和EPS完全恢复的原则,构建恢复负荷总量一致约束条件表达式如下:
Figure BDA0003134296620000083
式中,si,t表示t时段第i个负荷是否被恢复的二元决策变量;
Figure BDA0003134296620000084
分别为t时段第i个负荷是否被风光储和EPS恢复的二元决策变量。
可优选的是,在步骤S3中,建立交通路网模型支撑下的应急电源车优化调度方法具体包括以下步骤:
将应急供电平台获取的路网实时交通数据,带入通过Dijkstra算法求解EPS从预放车库到目的地候选接入节点的最短距离路径的表达式中,得到交通融合系数的具体表达式为:
Vd=VF·e-1.7r
w1=1/Vd
d=w0·w1
t=d/Vd
其中,Vd表示EPS车辆通行速度;VF为道路的零流速度;d为融合后得到的路径长度;w0为道路初始权值;w1为道路拥堵系数;r表示道路的拥塞情况,可由道路损害评估模块得到;t为EPS车辆调度时间。
本发明与现有技术相比,具有如下优点:
1.本发明采用的双层规划方法分别以关键负荷恢复均衡性最佳和配电网经济效益最大为上下两层的目标函数,通过上下层模型的关联和反馈实现了考虑孤岛均衡恢复与EPS优化调度的整体全局寻优,效果良好。
2.本发明基于双层模型和dijkstra算法解决了EPS优化调度问题,即候选接入节点的选择和调度路径的优化,减少了因调度产生的负荷中断时间,提高了供电可靠性和经济效益。
3.本发明所提均衡恢复策略可有效协调DER供电范围与关键负荷供电时长的关系,使一个停电周期内同等权重关键负荷恢复时间趋同,解决了因供电不公平带来的满意度问题,使恢复结果更合理,实用性高。
附图说明
图1为本发明考虑应急电源车调度的有源配电网故障均衡供电恢复方法的总体框架图;
图2为本发明考虑应急电源车调度的有源配电网故障均衡供电恢复方法的恢复策略流程图;
图3为本发明考虑应急电源车调度的有源配电网故障均衡供电恢复方法中EPS实时分配调度所需数据集与任务模块关系图。
具体实施方式
为详尽本发明之技术内容、结构特征、所达成目的及功效,以下将结合说明书附图进行详细说明。
本发明以实现各DER供电范围与关键负荷供电时长的均衡,使一个停电周期内具有同等权重的关键负荷持续恢复时间趋同于平均恢复时间指标;同时,制定EPS的调度策略,通过选择最佳候选接入节点和优化调度路径来减少负荷中断时间,最大化配电网的经济效益;最后对所提优化模型进行线性化处理并求解,本发明总体框架如图1所示,具体实施步骤如下:
S1、根据有源配电网特征,分别对配电网的负荷特性、储能系统的充放电特性、应急电源车(Emergency Power Supply Vehicle,EPS)的调度特性和分布式电源(Distributed Energy Resource,DER)的均衡供电特性进行建模。
S2、建立考虑应急电源调度的动态孤岛均衡恢复双层优化模型。
S3、基于步骤S2中双层优化模型和Dijkstra(迪杰斯特拉算法)算法,建立交通路网模型支撑下的应急电源车优化调度方法,以得到其最佳候选接入节点和最优调度路径。
将应急供电平台获取的路网实时交通数据,带入通过Dijkstra算法求解EPS从预放车库到目的地候选接入节点的最短距离路径的表达式中,得到交通融合系数的具体表达式为:
Vd=VF·e-1.7r
w1=1/Vd
d=w0·w1
t=d/Vd
其中,Vd表示EPS车辆通行速度;VF为道路的零流速度;d为融合后得到的路径长度;w0为道路初始权值;w1为道路拥堵系数;r表示道路的拥塞情况,可由道路损害评估模块得到;t为EPS车辆调度时间。
通过VR模块找到总权值最短的路径,并计算出EPS从预放车库到各候选接入节点的行驶时间,结合EPS的数量、容量、孤岛划分的范围以及候选接入节点的位置等影响因素开展EPS的优化调度。并基于步骤S2中双层模型求得EPS的最优候选接入节点和调度路径。完成其实时分配调度所需的数据集与任务模块之间的关系如图3所示。
S4、利用Big-M方法对步骤S24中双层优化模型进行线性化处理,将模型转换为混合整数线性规划,并利用Cplex求解得到全局最优故障恢复方案。
进一步的,步骤S1中的根据有源配电网特征,分别对配电网的负荷特性、储能系统的充放电特性、应急电源车的调度特性和分布式电源的均衡供电特性进行建模,具体过程如下:
S11、根据负荷时变性和负荷等级建立配电网的负荷特性模型,负荷时变性是根据各时段的用电需求,对节点日负荷时变曲线进行积分,具体表达式为:
Figure BDA0003134296620000101
式中,Lk(a)为节点k在a到a+Δa时段的负荷功率;fk(t)为节点k的负荷预测曲线;
同时,依据负荷停供对用户经济损失、人身安全等造成的影响将负荷等级分为一级、二级和三级。其中一级负荷需在配电网发生故障后完全恢复供电,而二级可控负荷和三级负荷,可在供电资源不足时中断或少量供电,以保证关键负荷的用电需求。
S12、通过储能(Energy Storage Systems,ESS)调度策略平抑风光机组出力波动性和负荷不确定性,建立储能系统的充放电特性模型,用充电状态(State-of-Charge,SOC)表示,其具体表达式为:
Figure BDA0003134296620000111
式中,SOCk,t为ESSk在t时段的SOC状态;
Figure BDA0003134296620000112
分别为ESSk在t-1时段内的充放电功率;ηch、ηdch分别为ESS的充放电效率;Ek为ESSk的容量;RES为DER中ESS的集合;Δt=1h。
根据充放电特性模型的公式求得各储能SOC需满足相关状态约束和充放电功率约束,具体表达式为:
Figure BDA0003134296620000113
Figure BDA0003134296620000114
Figure BDA0003134296620000115
Figure BDA0003134296620000116
式中,
Figure BDA0003134296620000117
分别为ESSk中SOC最小值和最大值;
Figure BDA0003134296620000118
分别为ESSk单位时段内最大充放电功率;
Figure BDA0003134296620000119
为二元变量,分别表示充放电状态,
Figure BDA00031342966200001110
为1表示充电,否则为0;
Figure BDA00031342966200001111
为1表示放电,否则为0。假设ESS有足够的无功补偿容量;
S13、分别根据应急电源车的充放电管理约束和应急电源车的候选接入节点配置约束建立应急电源车的调度特性模型:
在步骤S13中,应急电源车的充放电管理约束的具体表达式为:
Figure BDA00031342966200001112
Figure BDA00031342966200001113
Figure BDA00031342966200001114
Figure BDA00031342966200001115
Figure BDA0003134296620000121
式中,
Figure BDA0003134296620000122
分别为EPSk在时段t内的充放电功率;ηc,EPS、ηd,EPS分别为EPS的充放电效率;
Figure BDA0003134296620000123
分别为EPSk单位时段内最大充放电功率;
Figure BDA0003134296620000124
为EPSk的容量;
Figure BDA0003134296620000125
为二元变量,分别表示EPS的充放电状态;
Figure BDA0003134296620000126
Figure BDA0003134296620000127
分别表示EPS在故障时刻的初始荷电状态、故障结束时刻的荷电状态;REPS为DER中EPS的集合;
应急电源车的候选接入节点配置约束的具体表达式为:
Figure BDA0003134296620000128
Figure BDA0003134296620000129
Figure BDA00031342966200001210
Figure BDA00031342966200001211
式中,
Figure BDA00031342966200001212
为EPS候选接入点变量;vbk,t为节点动态孤岛划分变量;Q为配电网内可接入应急电源车供电的候选节点集合;P为包括预放车库在内的候选节点集合,
Figure BDA00031342966200001213
b∈Q,a∈P;R为分布式电源数量集合,R={1,2,....,m},k∈R。
S14、分别从平均恢复时间指标、第k个风光储的最大持续供电时间
Figure BDA00031342966200001219
和第k辆EPS的最大持续供电时间
Figure BDA00031342966200001214
建立分布式能源的均衡供电模型:
在步骤S14中,平均恢复时间指标由Tnet表示,其具体表达式为:
Figure BDA00031342966200001215
Figure BDA00031342966200001216
式中:
Figure BDA00031342966200001217
为DER在t时段内的状态;si,t表示负荷i的连接状态;
Figure BDA00031342966200001218
为第k个风光储系统在t时段内的最大出力;Pi,t为第i个负荷在t时段内功率需求;m为配电网中DER的数量;n为配电网中失电负荷量;I表示失电负荷节点集合,I={1,2,....,n},i∈I;R表示分布式电源数量集合,R={1,2,....,m},k∈R;t0为停电起始时刻;
将第k个风光储系统作为电源,以最大化负荷恢复价值为目标函数进行孤岛划分,并计算第k个DER能够持续为该孤岛内负荷供电的时间,得到第k个风光储的最大持续供电时间
Figure BDA0003134296620000131
的具体表达式为:
Figure BDA0003134296620000132
Figure BDA0003134296620000133
式中,
Figure BDA0003134296620000134
为第k个风光储系统在t时段内的状态;
Figure BDA0003134296620000135
为节点i分配变量,若负荷i在t时段内由第k个风光储系统恢复,则
Figure BDA0003134296620000136
否则
Figure BDA0003134296620000137
根据EPS的调度得到最佳候选接入节点,并形成最优孤岛为负荷供电,得到第k辆EPS的最大持续供电时间
Figure BDA0003134296620000138
的具体表达式为:
Figure BDA0003134296620000139
Figure BDA00031342966200001310
式中,
Figure BDA00031342966200001311
为第k个EPS在t时段内的状态;
Figure BDA00031342966200001312
为应急电源车最优匹配决策变量。
在步骤S14中,从优化负荷恢复时间的角度看,
Figure BDA00031342966200001313
表示在一个停电周期内各孤岛中关键负荷恢复的公平程度越来越高,当
Figure BDA00031342966200001314
表示完全公平恢复关键负荷。
进一步的,步骤S2中的建立动态孤岛均衡恢复双层优化模型的方法包括,如图2所示:
S21、以最大化负荷恢复水平和配电网运行效能建立下层模型的目标函数,具体表达式为:
Figure BDA00031342966200001315
式中,wi为第i个负荷的权重因子;ci为单位时间内第i个负荷缺供造成的损失费用;cEPS为EPS为负荷供电时单位容量消耗费用;sEPS为EPS调度过程中单位时间产生的车程费用;
Figure BDA00031342966200001316
为EPS候选接入点变量;
Figure BDA00031342966200001317
为t时段第k辆EPS从候选接入点a调运至候选接入点b的行驶时间;Q为配电网内可接入应急电源车供电的候选节点集合;P为预放车库在内的候选节点集合,
Figure BDA0003134296620000141
b∈Q,a∈P;
Figure BDA0003134296620000142
为在t时段内EPSk为第i个负荷提供的容量;
Figure BDA0003134296620000143
为EPS应急供电决策变量;si,t=1为t时段第i个负荷由风光储或EPS恢复;
上式中的目标函数是两个部分的加权和。第一项是最大限度提高失电负荷的恢复价值,尤其是保证关键负荷的供电。第二项是衡量总消耗费用的度量,目标使负荷停电损失费用、EPS供电能耗费用以及调度费用最小。总体目标是最大化配电网经济效益。故障后,将EPS调度到最佳候选接入节点,与DER协同形成各时段最佳孤岛为负荷供电。
S22、根据配电网孤岛恢复均衡性最优建立上层模型的目标函数,具体表达式为:
Figure BDA0003134296620000144
式中,
Figure BDA0003134296620000145
表示第k个风光储的最大持续供电时间;
Figure BDA0003134296620000146
表示第k辆EPS的最大持续供电时间;Tnet表示配电网平均恢复时间指标;
S23、从配电网负荷的节点电压、支路功率、孤岛潮流平衡、配电网辐射运行、孤岛划分和失电负荷恢复总量一致等六个因素建立动态孤岛均衡恢复双层优化模型的约束条件:
S231、节点电压的约束条件表达式如下:
Ui min≤Ui≤Ui max,i=1,2,…,m
式中,Ui min和Ui max分别为节点i电压上下限;
S232、支路功率的约束条件表达式如下:
Figure BDA0003134296620000147
式中,pg max为支路g能够允许的最大功率;
S233、根据线性化的DistFlow公式,建立含DER的有源配电网潮流方程,得到孤岛潮流平衡的约束条件表达式如下:
Figure BDA0003134296620000148
Figure BDA0003134296620000149
Figure BDA0003134296620000151
Figure BDA0003134296620000152
式中,βi,j,t表示支路状态;
Figure BDA0003134296620000153
时段t内支路g的有功功率;
Figure BDA0003134296620000154
时段t内支路g的无功功率;ri支路i的电阻;xi为支路i的电抗;V1为参考电压,取V1=1;σ为电压偏差,取σ=5%;g(·,i)为功率流入母线i的支路集合;g(i,·)为功率流出母线i的支路集合;Ig为支路g直接相连的母线集合;
S234、配电网辐射运行的约束条件表达式如下:
G∈g
式中,G为重构后的配电网网络拓扑,g为辐射状网络拓扑集合;
S235、将由节点j处的风光储或EPS供电的子网称为孤岛,孤岛划分的约束条件表达式如下:
Figure BDA0003134296620000155
Figure BDA0003134296620000156
式中,zj表示节点j是否有电源,zj=1表示节点j与风光储系统或EPS连接;Q为EPS连接的候选接入节点集合;F为风光储连接的根节点集合。
具体而言,在进行动态孤岛划分时,应确保各时段每个节点最多属于一个孤岛内,且每一时段各孤岛内必须有一个电源,每个电源节点作为根节点形成孤岛,具体表达式为:
Figure BDA0003134296620000157
Figure BDA0003134296620000158
式中,vik,t为节点动态孤岛划分变量,当t时段负荷节点i属于孤岛k时,vik,t=1,否则为0;上两式分别表示与电源连接的每个节点属于由其自身供电的孤岛内;没有电源不能形成孤岛;
当其父节点j属于孤岛k时,节点i才可以属于孤岛k,此时线路(i,j)也属于孤岛k,即线路处于闭合状态,否则被断开,具体表达式为:
Figure BDA0003134296620000159
Figure BDA0003134296620000161
式中,cij,t为支路的连通状态,若t时段内支路(i,j)闭合,则cij,t=1,否则cij,t=0;j为节点i的父节点;h为孤岛k中任意支路的子节点。
S236、根据任意时段划分入孤岛恢复范围内的负荷都需由风光储和EPS完全恢复的原则,构建恢复负荷总量一致约束条件表达式如下:
Figure BDA0003134296620000162
式中,si,t表示t时段第i个负荷是否被恢复的二元决策变量;
Figure BDA0003134296620000163
分别为t时段第i个负荷是否被风光储和EPS恢复的二元决策变量。
S24、根据步骤S21、S22和S23得到的目标函数和约束函数建立动态孤岛均衡恢复双层优化模型,具体表达式为:
Figure BDA0003134296620000164
式中,si,t表示t时段第i个负荷是否被恢复的二元决策变量;
Figure BDA0003134296620000165
分别为t时段第i个负荷是否被风光储和EPS恢复的二元决策变量;
Figure BDA0003134296620000166
为EPS候选接入点决策变量。
进一步的,步骤S4中的将双层优化模型进行线性化处理,转换为混合整数线性规划,利用Cplex求解得到全局最优故障恢复方案的方法包括:
S41、利用Big-M方法等效线性化处理步骤S14的表达式分别为:
Figure BDA0003134296620000167
Figure BDA0003134296620000171
Figure BDA0003134296620000172
式中,x1~x6为6个0-1变量;M为无穷大常数;si,t表示t时段第i个负荷是否被恢复的二元决策变量;
Figure BDA0003134296620000173
分别为t时段第i个负荷是否被风光储和EPS恢复的二元决策变量;
Figure BDA0003134296620000174
为EPS候选接入点决策变量;Pi,t为第i个负荷在t时段内功率需求;
Figure BDA0003134296620000175
为DER在t时段内的状态;
Figure BDA0003134296620000176
为EPSk在时段t内的放电功率;
Figure BDA0003134296620000177
为第k个风光储系统在t时段内的最大出力;
S42、将辅助二进制变量
Figure BDA0003134296620000178
替代二进制变量
Figure BDA0003134296620000179
并增加等效逻辑约束,对步骤S21表达式中x×y类型的非线性约束实现线性等效转换,等效逻辑约束的表达式如下:
Figure BDA00031342966200001710
S43、对步骤S24建立的动态孤岛均衡恢复双层优化模型中的相关非线性约束进行替换,并利用CPLEX求解。
以下结合实施例对本发明一种考虑应急电源车调度的有源配电网故障均衡供电恢复方法做进一步描述:
以某地震灾害事件为背景,以修改的PG&E 69节点标准配电系统及其灾后交通路网作为测试系统来验证本发明所提的有源配电网故障均衡恢复方法的有效性,具体实施过程如下:
S1、根据某地震灾害区的有源配电网特征,分别对配电网的负荷特性、储能系统的充放电特性、应急电源车的调度特性和分布式电源的均衡供电特性进行建模:
S11、假设接入有源配电网的DER包括PV、WT、ESS和EPS。其中ESS与PV、WT组成风光储系统,通过调度ESS平抑DER出力的波动性和负荷时变性,建立负荷特性与储能系统充放电模型,根据负荷时变性和负荷等级建立配电网的负荷特性模型,负荷时变性是根据某地震灾害区季节、作息时刻的规律性变化,因其短时波动较小,各时段的用电需求,对节点日负荷时变曲线进行积分,具体表达式为:
Figure BDA0003134296620000181
式中,Lk(a)为节点k在a到a+Δa时段的负荷功率;fk(t)为节点k的负荷预测曲线;
同时,依据负荷停供对用户经济损失、人身安全等造成的影响将负荷等级分为一级、二级和三级。其中一级负荷的缺供将严重威胁用户的生命财产安全,因此,该类关键负荷需在配电网发生故障后完全恢复供电,而针对配电网中存在的二级可控负荷和三级负荷,可在供电资源不足时中断或少量供电,以保证关键负荷的用电需求。
S12、通过储能调度策略平抑风光机组出力波动性和负荷不确定性,建立储能系统的充放电特性模型,用充电状态(SOC)表示,其具体表达式为:
Figure BDA0003134296620000182
式中,SOCk,t为ESSk在t时段的SOC状态;
Figure BDA0003134296620000183
分别为ESSk在t-1时段内的充放电功率;ηch、ηdch分别为ESS的充放电效率;Ek为ESSk的容量;RES为DER中ESS的集合;Δt=1h。
根据充放电特性模型的公式求得各储能SOC需满足相关状态约束和充放电功率约束,具体表达式为:
Figure BDA0003134296620000184
Figure BDA0003134296620000185
Figure BDA0003134296620000186
Figure BDA0003134296620000187
式中,
Figure BDA0003134296620000188
分别为ESSk中SOC最小值和最大值;
Figure BDA0003134296620000189
分别为ESSk单位时段内最大充放电功率;
Figure BDA00031342966200001810
为二元变量,分别表示充放电状态,
Figure BDA0003134296620000191
为1表示充电,否则为0;
Figure BDA0003134296620000192
为1表示放电,否则为0。假设ESS有足够的无功补偿容量。
S13、为保证某地震灾害区故障周期内关键负荷的全部恢复,并尽量减少因EPS调度时长带来的负荷损失,分别根据应急电源车的充放电管理约束和应急电源车的候选接入节点配置约束建立应急电源车的调度特性模型:
应急电源车的充放电管理约束的具体表达式为:
Figure BDA0003134296620000193
Figure BDA0003134296620000194
Figure BDA0003134296620000195
Figure BDA0003134296620000196
Figure BDA0003134296620000197
式中,
Figure BDA0003134296620000198
分别为EPSk在时段t内的充放电功率;ηc,EPS、ηd,EPS分别为EPS的充放电效率;
Figure BDA0003134296620000199
分别为EPSk单位时段内最大充放电功率;
Figure BDA00031342966200001910
为EPSk的容量;
Figure BDA00031342966200001911
为二元变量,分别表示EPS的充放电状态;
Figure BDA00031342966200001912
Figure BDA00031342966200001913
分别表示EPS在故障时刻的初始荷电状态、故障结束时刻的荷电状态;REPS为DER中EPS的集合;
应急电源车的候选接入节点配置约束的具体表达式为:
Figure BDA00031342966200001914
Figure BDA00031342966200001915
Figure BDA00031342966200001916
Figure BDA00031342966200001917
式中,
Figure BDA00031342966200001918
为EPS候选接入点变量,当t时段第k辆EPS从候选接入点a调运至候选接入点b为孤岛供电时
Figure BDA00031342966200001919
否则
Figure BDA00031342966200001920
vbk,t为节点动态孤岛划分变量,当t时段负荷节点b属于孤岛k时,vbk,t=1,否则为0;其中Q表示配电网内可接入应急电源车供电的候选节点集合,P表示包括预放车库在内的候选节点集合,
Figure BDA00031342966200001921
b∈Q,a∈P,初始时刻时a=0表示EPS从预放车库调出,R表示分布式电源数量集合,R={1,2,....,m},k∈R;以上几式分别表示当调度第k辆EPS作为孤岛内的主电源,确保其连接的候选接入节点属于该孤岛内;每个EPS供电孤岛内至少有一个EPS候选接入节点;每个候选节点最多分配一辆EPS;任意时刻一辆EPS有且只能接入一个孤岛内的候选节点。
S14、为了避免某地震灾害区恢复策略制定过程中出现因同等权重的关键负荷可持续恢复供电时间不均衡带来的用户满意度下降问题,分别从平均恢复时间指标、第k个风光储的最大持续供电时间
Figure BDA0003134296620000201
和第k辆EPS的最大持续供电时间
Figure BDA0003134296620000202
建立分布式能源的均衡供电模型:
为了使某地震灾害区DER的持续供电时间尽可能接近该指标,从而使某地震灾害区DER供电孤岛内同等权重的关键负荷恢复时间公平,提高某地震灾害区用户的满意度。令平均恢复时间指标由Tnet表示,其具体表达式为:
Figure BDA0003134296620000203
Figure BDA0003134296620000204
式中:
Figure BDA0003134296620000205
为DER在t时段内的状态,若DER可以提供的最大电量满足能够被恢复的负荷总需求量,则
Figure BDA0003134296620000206
否则
Figure BDA0003134296620000207
si,t表示负荷i的连接状态,若第i个负荷在t时段内由风光储或EPS恢复,则si,t=1,否则si,t=0;
Figure BDA0003134296620000208
为第k个风光储系统在t时段内的最大出力;Pi,t为第i个负荷在t时段内功率需求;m为某地震灾害区配电网中DER的数量;n为配电网中失电负荷量。I表示失电负荷节点集合,I={1,2,....,n},i∈I;R表示分布式电源数量集合,R={1,2,....,m},k∈R。假设停电时长不超过48小时;t0为停电起始时刻;
将第k个风光储系统作为电源,以最大化负荷恢复价值为目标函数进行孤岛划分,并计算第k个DER能够持续为该孤岛内负荷供电的时间,得到第k个风光储的最大持续供电时间
Figure BDA0003134296620000209
的具体表达式为:
Figure BDA00031342966200002010
Figure BDA00031342966200002011
式中,
Figure BDA00031342966200002012
为第k个风光储系统在t时段内的状态;
Figure BDA00031342966200002013
为节点i分配变量,若负荷i在t时段内由第k个风光储系统恢复,则
Figure BDA00031342966200002014
否则
Figure BDA00031342966200002015
根据EPS的调度得到最佳候选接入节点,并形成最优孤岛为负荷供电,得到第k辆EPS的最大持续供电时间
Figure BDA0003134296620000211
的具体表达式为:
Figure BDA0003134296620000212
Figure BDA0003134296620000213
式中,
Figure BDA0003134296620000214
为第k个EPS在t时段内的状态;
Figure BDA0003134296620000215
为应急电源车最优匹配决策变量,若负荷i在t时段内由第k辆EPS恢复,则
Figure BDA0003134296620000216
否则
Figure BDA0003134296620000217
S2、在步骤S1相关特性分析基础上建立考虑应急电源调度的动态孤岛均衡恢复双层优化模型,上层模型以各孤岛恢复均衡性最优为目标函数;下层模型以最大化负荷恢复价值和配电网经济效益为目标函数;并考虑相关约束条件:
S21、以最大化负荷恢复水平和配电网运行效能建立下层模型的目标函数,具体表达式为:
Figure BDA0003134296620000218
式中,wi为第i个负荷的权重因子,关键负荷和非关键负荷由权重因子的大小值指示;ci为单位时间内第i个负荷缺供造成的损失费用;cEPS为EPS为负荷供电时单位容量消耗费用;sEPS为EPS调度过程中单位时间产生的车程费用;
Figure BDA0003134296620000219
为EPS候选接入点变量,当t时段第k辆EPS从候选接入点a调运至候选接入点b为孤岛供电时
Figure BDA00031342966200002110
否则
Figure BDA00031342966200002111
为t时段第k辆EPS从候选接入点a调运至候选接入点b的行驶时间,也可理解为负荷缺供时间,其中Q表示配电网内可接入应急电源车供电的候选节点集合,P表示包括预放车库在内的候选节点集合,
Figure BDA00031342966200002112
b∈Q,a∈P,初始时刻时a=0表示EPS从预放车库调出;
Figure BDA00031342966200002113
为在t时段内EPSk为第i个负荷提供的容量;
Figure BDA00031342966200002114
为EPS应急供电决策变量,若在t时段内负荷i由EPSk恢复,则
Figure BDA00031342966200002115
否则
Figure BDA00031342966200002116
si,t=1时表示t时段第i个负荷由风光储或EPS恢复,当负荷i由风光储恢复时,它的中断时间在这里被认为是零,则方括号里的整体项都为零,不需要考虑总消耗费用。
S22、根据某地震灾害区配电网孤岛恢复均衡性最优建立上层模型的目标函数,具体表达式为:
Figure BDA0003134296620000221
式中,
Figure BDA0003134296620000222
表示第k个风光储的最大持续供电时间;
Figure BDA0003134296620000223
表示第k辆EPS的最大持续供电时间;Tnet表示配电网平均恢复时间指标。
S23、从某地震灾害区配电网负荷的节点电压、支路功率、孤岛潮流平衡、配电网辐射运行、孤岛划分和失电负荷恢复总量一致等六个因素建立动态孤岛均衡恢复双层优化模型的约束条件:
S231、节点电压的约束条件表达式如下:
Ui min≤Ui≤Ui max,i=1,2,…,m
式中,Ui min和Ui max分别为节点i电压上下限;
S232、支路功率的约束条件表达式如下:
Figure BDA0003134296620000224
式中,pg max为支路g能够允许的最大功率;
S233、根据线性化的DistFlow公式,建立含DER的有源配电网潮流方程,得到孤岛潮流平衡的约束条件表达式如下:
Figure BDA0003134296620000225
Figure BDA0003134296620000226
Figure BDA0003134296620000227
Figure BDA0003134296620000228
式中,βi,j,t表示支路状态,当有潮流由支路i流向支路j时βi,j,t为1,否则为0;
Figure BDA0003134296620000229
时段t内支路g的有功功率;
Figure BDA00031342966200002210
时段t内支路g的无功功率;ri支路i的电阻;xi为支路i的电抗;V1为参考电压,取V1=1,应用大M方法对两条不相连母线间电压解耦;σ为电压偏差,取σ=5%;g(·,i)为功率流入母线i的支路集合;g(i,·)为功率流出母线i的支路集合;Ig为支路g直接相连的母线集合;
S234、某地震灾害区配电网辐射运行的约束条件表达式如下:
G∈g
式中,G为重构后的配电网网络拓扑,g为辐射状网络拓扑集合;
S235、将由节点j处的风光储或EPS供电的子网称为孤岛,孤岛划分的约束条件表达式如下:
Figure BDA0003134296620000231
Figure BDA0003134296620000232
式中,zj表示节点j是否有电源,zj=1表示节点j与风光储系统或EPS连接;Q为EPS连接的候选接入节点集合;F为风光储连接的根节点集合。
S236、根据任意时段划分入孤岛恢复范围内的负荷都需由风光储和EPS完全恢复的原则,构建恢复负荷总量一致约束条件表达式如下:
Figure BDA0003134296620000233
式中,si,t表示t时段第i个负荷是否被恢复的二元决策变量;
Figure BDA0003134296620000234
分别为t时段第i个负荷是否被风光储和EPS恢复的二元决策变量。
约束条件除上述外还包括前述步骤S12和S13中的表达式。
S24、根据步骤S12、S13、S21、S22和S23得到的目标函数和约束函数建立动态孤岛均衡恢复双层优化模型,具体表达式为:
Figure BDA0003134296620000235
式中,si,t表示t时段第i个负荷是否被恢复的二元决策变量;
Figure BDA0003134296620000236
分别为t时段第i个负荷是否被风光储和EPS恢复的二元决策变量;
Figure BDA0003134296620000237
为EPS候选接入点决策变量。
S3、基于步骤S2中双层优化模型和Dijkstra(迪杰斯特拉算法)算法,建立交通路网模型支撑下的应急电源车优化调度方法,以得到其最佳候选接入节点和最优调度路径。
以应急供电平台中的VR模块为载体,获取路网实时交通数据,通过Dijkstra算法求解EPS从预放车库到目的地候选接入节点的最短距离路径,此处的“距离”实为各路段的权值,兼顾了实际道路长度与交通信息(即道路的通行速度、拥堵率)等因素的影响,并通过交通融合系数来反映,得到交通融合系数的具体表达式为:
Vd=VF·e-1.7r
w1=1/Vd
d=w0·w1
t=d/Vd
其中,Vd表示EPS车辆通行速度;VF为道路的零流速度;d为融合后得到的路径长度;w0为道路初始权值;w1为道路拥堵系数;r表示道路的拥塞情况,可由道路损害评估模块得到;t为EPS车辆调度时间。
S4、利用Big-M方法对步骤S24中双层优化模型进行线性化处理,将模型转换为混合整数线性规划,并利用Cplex求解得到全局最优故障恢复方案。
S41、利用Big-M方法等效线性化处理步骤S14的约束方程,同时,引入x1、x2、x3、x4、x5和x6变量,等效线性化处理后的表达式分别为:
Figure BDA0003134296620000241
Figure BDA0003134296620000242
Figure BDA0003134296620000243
式中,x1~x6为6个0-1变量;M为无穷大常数;si,t表示t时段第i个负荷是否被恢复的二元决策变量;
Figure BDA0003134296620000244
分别为t时段第i个负荷是否被风光储和EPS恢复的二元决策变量;
Figure BDA0003134296620000251
为EPS候选接入点决策变量;Pi,t为第i个负荷在t时段内功率需求;
Figure BDA0003134296620000252
为DER在t时段内的状态;
Figure BDA0003134296620000253
为EPSk在时段t内的放电功率;
Figure BDA0003134296620000254
为第k个风光储系统在t时段内的最大出力;
S42、针对步骤S21中表达式中x×y类型的非线性约束,将辅助二进制变量
Figure BDA0003134296620000255
替代二进制变量
Figure BDA0003134296620000256
并增加等效逻辑约束,对步骤S21表达式中x×y类型的非线性约束实现线性等效转换,等效逻辑约束的表达式如下:
Figure BDA0003134296620000257
S43、对步骤S24建立的动态孤岛均衡恢复双层优化模型中的相关非线性约束进行替换,将动态孤岛均衡恢复模型转化为混合整数线性规划问题,并利用CPLEX求解。
停电周期内考虑EPS优化调度的配电网均衡供电调节结果见下表1所示。
表1考虑EPS调度的均衡供电调节结果
Figure BDA0003134296620000258
Figure BDA0003134296620000261
由表1可知,利用本发明建立的双层模型形成了各时段最优的孤岛供电范围,保证了关键负荷的100%复供,且所有DER对各关键负荷的最大持续供电时间均为12小时左右,保证了故障周期内同等权重用户的用电公平,同时在解决均衡恢复问题时结合dijkstra算法求解得到EPS最优调度结果。
以上所述的实施例仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。

Claims (6)

1.一种考虑应急电源车调度的有源配电网故障均衡供电恢复方法,其特征在于,具体实施步骤如下:
S1、根据有源配电网特征,分别对配电网的负荷特性、储能系统的充放电特性、应急电源车的调度特性和分布式电源的均衡供电特性进行建模:
S11、根据负荷时变性和负荷等级建立配电网的负荷特性模型,所述负荷时变性是根据各时段的用电需求,对节点日负荷时变曲线进行积分,具体表达式为:
Figure FDA0003134296610000011
式中,Lk(a)为节点k在a到a+Δa时段的负荷功率;fk(t)为节点k的负荷预测曲线;
同时,依据负荷停供对用户经济损失、人身安全造成的影响将负荷等级分为一级、二级和三级;
S12、通过储能调度策略平抑风光机组出力波动性和负荷不确定性,建立储能系统的充放电特性模型,根据充放电特性模型的公式求得各储能SOC需满足相关状态约束和充放电功率约束,具体表达式为:
Figure FDA0003134296610000012
Figure FDA0003134296610000013
Figure FDA0003134296610000014
Figure FDA0003134296610000015
式中,
Figure FDA0003134296610000016
分别为ESSk中SOC最小值和最大值;
Figure FDA0003134296610000017
分别为ESSk单位时段内最大充放电功率;
Figure FDA0003134296610000018
为二元变量,分别表示充放电状态,
Figure FDA0003134296610000019
为1表示充电,否则为0;
Figure FDA00031342966100000110
为1表示放电,否则为0;假设ESS有足够的无功补偿容量;
S13、分别根据应急电源车的充放电管理约束和应急电源车的候选接入节点配置约束建立应急电源车的调度特性模型;
S14、分别从平均恢复时间指标、第k个风光储的最大持续供电时间
Figure FDA00031342966100000111
和第k辆EPS的最大持续供电时间
Figure FDA0003134296610000021
建立分布式能源的均衡供电模型;
S2、建立考虑应急电源调度的动态孤岛均衡恢复双层优化模型:
S21、以最大化负荷恢复水平和配电网运行效能建立下层模型的目标函数,具体表达式为:
Figure FDA0003134296610000022
式中,wi为第i个负荷的权重因子;ci为单位时间内第i个负荷缺供造成的损失费用;cEPS为EPS为负荷供电时单位容量消耗费用;sEPS为EPS调度过程中单位时间产生的车程费用;
Figure FDA0003134296610000023
为EPS候选接入点变量;
Figure FDA0003134296610000024
为t时段第k辆EPS从候选接入点a调运至候选接入点b的行驶时间;Q为配电网内可接入应急电源车供电的候选节点集合;P为预放车库在内的候选节点集合,
Figure FDA00031342966100000211
Figure FDA0003134296610000025
为在t时段内EPSk为第i个负荷提供的容量;
Figure FDA0003134296610000026
为EPS应急供电决策变量;si,t=1为t时段第i个负荷由风光储或EPS恢复;
S22、根据配电网孤岛恢复均衡性最优建立上层模型的目标函数,具体表达式为:
Figure FDA0003134296610000027
式中,
Figure FDA0003134296610000028
表示第k个风光储的最大持续供电时间;
Figure FDA0003134296610000029
表示第k辆EPS的最大持续供电时间;Tnet表示配电网平均恢复时间指标;
S23、从配电网负荷的节点电压、支路功率、孤岛潮流平衡、配电网辐射运行、孤岛划分和失电负荷恢复总量一致六个因素建立动态孤岛均衡恢复双层优化模型的约束条件;
S24、根据步骤S21、S22和S23得到的目标函数和约束函数建立动态孤岛均衡恢复双层优化模型,具体表达式为:
Figure FDA00031342966100000210
式中,si,t表示t时段第i个负荷是否被恢复的二元决策变量;
Figure FDA0003134296610000031
分别为t时段第i个负荷是否被风光储和EPS恢复的二元决策变量;
Figure FDA0003134296610000032
为EPS候选接入点决策变量;
S3、基于步骤S2中双层优化模型和迪杰斯特拉算法,建立交通路网模型支撑下的应急电源车优化调度方法,以得到其最佳候选接入节点和最优调度路径;
S4、利用Big-M方法对步骤S24中双层优化模型进行线性化处理,将模型转换为混合整数线性规划,并利用Cplex求解得到全局最优故障恢复方案:
S41、利用Big-M方法等效线性化处理步骤S14的表达式分别为:
Figure FDA0003134296610000033
Figure FDA0003134296610000034
Figure FDA0003134296610000035
式中,x1~x6为6个0-1变量;M为无穷大常数;si,t表示t时段第i个负荷是否被恢复的二元决策变量;
Figure FDA0003134296610000036
分别为t时段第i个负荷是否被风光储和EPS恢复的二元决策变量;
Figure FDA0003134296610000037
为EPS候选接入点决策变量;Pi,t为第i个负荷在t时段内功率需求;
Figure FDA0003134296610000038
为DER在t时段内的状态;
Figure FDA0003134296610000039
为EPSk在时段t内的放电功率;
Figure FDA00031342966100000310
为第k个风光储系统在t时段内的最大出力;
S42、将辅助二进制变量
Figure FDA00031342966100000311
替代二进制变量
Figure FDA00031342966100000312
并增加等效逻辑约束,对步骤S21表达式中x×y类型的非线性约束实现线性等效转换,所述等效逻辑约束的表达式如下:
Figure FDA0003134296610000041
S43、对步骤S24建立的动态孤岛均衡恢复双层优化模型中的相关非线性约束进行替换,并利用CPLEX求解。
2.根据权利要求1所述的考虑应急电源车调度的有源配电网故障均衡供电恢复方法,其特征在于,在步骤S12中,所述储能系统的充放电特性模型用充电状态SOC表示,其具体表达式为:
Figure FDA0003134296610000042
式中,SOCk,t为ESSk在t时段的SOC状态;
Figure FDA0003134296610000043
分别为ESSk在t-1时段内的充放电功率;ηch、ηdch分别为ESS的充放电效率;Ek为ESSk的容量;RES为DER中ESS的集合;Δt=1h。
3.根据权利要求1所述的考虑应急电源车调度的有源配电网故障均衡供电恢复方法,其特征在于,在步骤S13中,所述应急电源车的充放电管理约束的具体表达式为:
Figure FDA0003134296610000044
Figure FDA0003134296610000045
Figure FDA0003134296610000046
Figure FDA0003134296610000047
Figure FDA0003134296610000048
式中,
Figure FDA0003134296610000049
分别为EPSk在时段t内的充放电功率;ηc,EPS、ηd,EPS分别为EPS的充放电效率;
Figure FDA00031342966100000410
分别为EPSk单位时段内最大充放电功率;
Figure FDA00031342966100000411
为EPSk的容量;
Figure FDA00031342966100000412
为二元变量,分别表示EPS的充放电状态;
Figure FDA00031342966100000413
Figure FDA00031342966100000414
分别表示EPS在故障时刻的初始荷电状态、故障结束时刻的荷电状态;REPS为DER中EPS的集合;
所述应急电源车的候选接入节点配置约束的具体表达式为:
Figure FDA00031342966100000415
Figure FDA00031342966100000416
Figure FDA00031342966100000417
Figure FDA0003134296610000051
式中,
Figure FDA0003134296610000052
为EPS候选接入点变量;vbk,t为节点动态孤岛划分变量;Q为配电网内可接入应急电源车供电的候选节点集合,P为包括预放车库在内的候选节点集合,
Figure FDA00031342966100000515
R为分布式电源数量集合,R={1,2,....,m},k∈R。
4.根据权利要求1所述的考虑应急电源车调度的有源配电网故障均衡供电恢复方法,其特征在于,在步骤S14中,所述平均恢复时间指标由Tnet表示,其具体表达式为:
Figure FDA0003134296610000053
Figure FDA0003134296610000054
式中:
Figure FDA0003134296610000055
为DER在t时段内的状态;si,t表示负荷i的连接状态;
Figure FDA0003134296610000056
为第k个风光储系统在t时段内的最大出力;Pi,t为第i个负荷在t时段内功率需求;m为配电网中DER的数量;n为配电网中失电负荷量;I表示失电负荷节点集合,I={1,2,....,n},i∈I;R表示分布式电源数量集合,R={1,2,....,m},k∈R;t0为停电起始时刻;
将第k个风光储系统作为电源,以最大化负荷恢复价值为目标函数进行孤岛划分,并计算第k个DER能够持续为该孤岛内负荷供电的时间,得到第k个风光储的最大持续供电时间
Figure FDA0003134296610000057
的具体表达式为:
Figure FDA0003134296610000058
Figure FDA0003134296610000059
式中,
Figure FDA00031342966100000510
为第k个风光储系统在t时段内的状态;
Figure FDA00031342966100000511
为节点i分配变量,若负荷i在t时段内由第k个风光储系统恢复,则
Figure FDA00031342966100000512
否则
Figure FDA00031342966100000513
根据EPS的调度得到最佳候选接入节点,并形成最优孤岛为负荷供电,得到第k辆EPS的最大持续供电时间
Figure FDA00031342966100000514
的具体表达式为:
Figure FDA0003134296610000061
Figure FDA0003134296610000062
式中,
Figure FDA0003134296610000063
为第k个EPS在t时段内的状态;
Figure FDA0003134296610000064
为应急电源车最优匹配决策变量。
5.根据权利要求1所述的考虑应急电源车调度的有源配电网故障均衡供电恢复方法,其特征在于,所述步骤S23具体包括以下步骤:
S231、节点电压的约束条件表达式如下:
Uimin≤Ui≤Uimax,i=1,2,…,m
式中,Uimin和Uimax分别为节点i电压上下限;
S232、支路功率的约束条件表达式如下:
Figure FDA0003134296610000065
式中,pgmax为支路g能够允许的最大功率;
S233、根据线性化的DistFlow公式,建立含DER的有源配电网潮流方程,得到孤岛潮流平衡的约束条件表达式如下:
Figure FDA0003134296610000066
Figure FDA0003134296610000067
Figure FDA0003134296610000068
Figure FDA0003134296610000069
式中,βi,j,t表示支路状态;
Figure FDA00031342966100000610
时段t内支路g的有功功率;
Figure FDA00031342966100000611
时段t内支路g的无功功率;ri支路i的电阻;xi为支路i的电抗;V1为参考电压,取V1=1;σ为电压偏差,取σ=5%;g(·,i)为功率流入母线i的支路集合;g(i,·)为功率流出母线i的支路集合;Ig为支路g直接相连的母线集合;
S234、配电网辐射运行的约束条件表达式如下:
G∈g
式中,G为重构后的配电网网络拓扑,g为辐射状网络拓扑集合;
S235、将由节点j处的风光储或EPS供电的子网称为孤岛,孤岛划分的约束条件表达式如下:
Figure FDA0003134296610000071
Figure FDA0003134296610000072
式中,zj表示节点j是否有电源,zj=1表示节点j与风光储系统或EPS连接;Q为EPS连接的候选接入节点集合;F为风光储连接的根节点集合;
S236、根据任意时段划分入孤岛恢复范围内的负荷都需由风光储和EPS完全恢复的原则,构建恢复负荷总量一致约束条件表达式如下:
Figure FDA0003134296610000073
式中,si,t表示t时段第i个负荷是否被恢复的二元决策变量;
Figure FDA0003134296610000074
分别为t时段第i个负荷是否被风光储和EPS恢复的二元决策变量。
6.根据权利要求1所述的考虑应急电源车调度的有源配电网故障均衡供电恢复方法,其特征在于,在步骤S3中,建立交通路网模型支撑下的应急电源车优化调度方法具体包括以下步骤:
将应急供电平台获取的路网实时交通数据,带入通过Dijkstra算法求解EPS从预放车库到目的地候选接入节点的最短距离路径的表达式中,得到交通融合系数的具体表达式为:
Vd=VF·e-1.7r
w1=1/Vd
d=w0·w1
t=d/Vd
其中,Vd表示EPS车辆通行速度;VF为道路的零流速度;d为融合后得到的路径长度;w0为道路初始权值;w1为道路拥堵系数;r表示道路的拥塞情况,可由道路损害评估模块得到;t为EPS车辆调度时间。
CN202110712635.4A 2021-06-25 2021-06-25 考虑应急电源车调度的有源配电网故障均衡供电恢复方法 Active CN113452051B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110712635.4A CN113452051B (zh) 2021-06-25 2021-06-25 考虑应急电源车调度的有源配电网故障均衡供电恢复方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110712635.4A CN113452051B (zh) 2021-06-25 2021-06-25 考虑应急电源车调度的有源配电网故障均衡供电恢复方法

Publications (2)

Publication Number Publication Date
CN113452051A true CN113452051A (zh) 2021-09-28
CN113452051B CN113452051B (zh) 2022-05-27

Family

ID=77813161

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110712635.4A Active CN113452051B (zh) 2021-06-25 2021-06-25 考虑应急电源车调度的有源配电网故障均衡供电恢复方法

Country Status (1)

Country Link
CN (1) CN113452051B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114120630A (zh) * 2021-10-15 2022-03-01 国网浙江省电力有限公司杭州供电公司 一种考虑道路交通动态特征的应急电源配置方法
CN114243901A (zh) * 2021-11-12 2022-03-25 国网浙江省电力有限公司杭州供电公司 一种城市电网故障快速响应系统及方法
CN114678881A (zh) * 2022-04-06 2022-06-28 四川大学 一种v2g辅助支撑下的地震灾后电网快速恢复方法
CN114818379A (zh) * 2022-05-24 2022-07-29 四川大学 考虑多维不确定性的信息物理主动配电系统互动规划方法
CN115036954A (zh) * 2022-08-12 2022-09-09 国网天津市电力公司物资公司 基于电动物流车协同调度的负荷应急恢复方法及系统
CN115995815A (zh) * 2023-03-23 2023-04-21 国网山西省电力公司电力科学研究院 一种基于多模块嵌套迭代的负荷故障恢复方法
CN116937599A (zh) * 2023-05-16 2023-10-24 东南大学 一种考虑5g基站备用储能的配电网负荷恢复方法及系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107332234A (zh) * 2017-06-30 2017-11-07 燕山大学 一种考虑可再生能源间歇性的主动配电网多故障修复策略
CN108988322A (zh) * 2018-06-30 2018-12-11 南京理工大学 考虑系统时变性的微网运行策略优化方法
CN110263435A (zh) * 2019-06-20 2019-09-20 燕山大学 基于电-气耦合综合能源系统的双层优化故障恢复方法
CN110401229A (zh) * 2019-06-30 2019-11-01 天津大学 一种考虑微能源网支撑作用的配电网弹性提升方法
CN111478358A (zh) * 2019-11-06 2020-07-31 南瑞集团有限公司 一种计及分布式电源不确定性的配电网鲁棒恢复决策方法
CN111539566A (zh) * 2020-04-21 2020-08-14 燕山大学 一种考虑灾前预调度的配电网多故障抢修恢复方法及系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107332234A (zh) * 2017-06-30 2017-11-07 燕山大学 一种考虑可再生能源间歇性的主动配电网多故障修复策略
CN108988322A (zh) * 2018-06-30 2018-12-11 南京理工大学 考虑系统时变性的微网运行策略优化方法
CN110263435A (zh) * 2019-06-20 2019-09-20 燕山大学 基于电-气耦合综合能源系统的双层优化故障恢复方法
CN110401229A (zh) * 2019-06-30 2019-11-01 天津大学 一种考虑微能源网支撑作用的配电网弹性提升方法
CN111478358A (zh) * 2019-11-06 2020-07-31 南瑞集团有限公司 一种计及分布式电源不确定性的配电网鲁棒恢复决策方法
CN111539566A (zh) * 2020-04-21 2020-08-14 燕山大学 一种考虑灾前预调度的配电网多故障抢修恢复方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
RICARDO PÉREZ-IBACACHE: "Linear State-Feedback Primary Control for Enhanced Dynamic Response of AC Microgrids", 《IEEE TRANSACTIONS ON SMART GRID》 *
卢志刚: "基于黑板模型的配电网多故障分时段动态恢复", 《电网技术》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114120630A (zh) * 2021-10-15 2022-03-01 国网浙江省电力有限公司杭州供电公司 一种考虑道路交通动态特征的应急电源配置方法
CN114243901A (zh) * 2021-11-12 2022-03-25 国网浙江省电力有限公司杭州供电公司 一种城市电网故障快速响应系统及方法
CN114678881A (zh) * 2022-04-06 2022-06-28 四川大学 一种v2g辅助支撑下的地震灾后电网快速恢复方法
CN114678881B (zh) * 2022-04-06 2023-03-07 四川大学 一种v2g辅助支撑下的地震灾后电网快速恢复方法
CN114818379A (zh) * 2022-05-24 2022-07-29 四川大学 考虑多维不确定性的信息物理主动配电系统互动规划方法
CN114818379B (zh) * 2022-05-24 2023-10-24 四川大学 考虑多维不确定性的信息物理主动配电系统互动规划方法
CN115036954A (zh) * 2022-08-12 2022-09-09 国网天津市电力公司物资公司 基于电动物流车协同调度的负荷应急恢复方法及系统
CN115995815A (zh) * 2023-03-23 2023-04-21 国网山西省电力公司电力科学研究院 一种基于多模块嵌套迭代的负荷故障恢复方法
CN115995815B (zh) * 2023-03-23 2023-06-13 国网山西省电力公司电力科学研究院 一种基于多模块嵌套迭代的负荷故障恢复方法
CN116937599A (zh) * 2023-05-16 2023-10-24 东南大学 一种考虑5g基站备用储能的配电网负荷恢复方法及系统

Also Published As

Publication number Publication date
CN113452051B (zh) 2022-05-27

Similar Documents

Publication Publication Date Title
CN113452051B (zh) 考虑应急电源车调度的有源配电网故障均衡供电恢复方法
CN103248064B (zh) 一种复合型能源充电储能系统及其方法
CN107332234B (zh) 一种考虑可再生能源间歇性的主动配电网多故障修复方法
CN107832905B (zh) 一种适应分布式发电和储能站发展的配电网规划方法
CN110866627B (zh) 一种考虑阶梯气价的多区域电-气耦合综合能源系统优化调度方法
CN108470239A (zh) 计及需求侧管理和储能的主动配电网多目标分层规划方法
CN105375507A (zh) 雾霾环境下虚拟电厂的电力两级交互优化调度系统
CN105139147A (zh) 微电网系统的经济调度方法
CN109948849A (zh) 一种计及储能接入的配电网网架规划方法
CN110676849B (zh) 一种孤岛微电网群能量调度模型的构建方法
CN109245091A (zh) 一种面向售电公司的双层策略调度方法
Ying et al. The flexible smart traction power supply system and its hierarchical energy management strategy
Tong et al. Multi-scenario and multi-objective collaborative optimization of distribution network considering electric vehicles and mobile energy storage systems
Pal et al. Planning of EV charging station with distribution network expansion considering traffic congestion and uncertainties
CN113746123A (zh) 一种光储充一体化电站系统
CN106779250B (zh) 一种基于新型优化模型的孤立分布式电网配置方法
CN105574681A (zh) 一种多时间尺度社区能源局域网能量调度方法
CN107622332A (zh) 一种基于静态安全性约束的电网侧储能容量优化配置方法
CN112531788B (zh) 考虑多重不确定性和自趋优运行的透明微网群规划方法
CN116961057A (zh) 一种考虑电动汽车的多时段配电网故障恢复方法
CN117134409A (zh) 考虑电-氢-热互补的微网系统及其多目标优化配置方法
CN111224418B (zh) 一种基于电动汽车储能的微电网调度方法及系统
CN116307562A (zh) 一种轨道交通自洽能源系统的光储规划配置方法及装置
CN114676913A (zh) 一种结合道路网与配电网的电力应急资源优化方法及系统
Hua et al. Coordinated optimal dispatch for active distribution network under multi-time scales

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