(文章复现)面向配电网韧性提升的移动储能预布局与动态调度策略(2)-灾后调度matlab代码

news/2024/4/19 0:41:50

参考文献:

[1]王月汉,刘文霞,姚齐,万海洋,何剑,熊雪君.面向配电网韧性提升的移动储能预布局与动态调度策略[J].电力系统自动化,2022,46(15):37-45.

1.基本原理

1. 1 目标函数

        在灾害发生后,配电网失去主网供电,设故障的持续时间可根据灾害强度和抢修资源数量做出预测。以故障持续时间内负荷削减功率加权值最小为目标,建立多源协同的灾后恢复优化模型,通过动态调度移动储能、电动汽车与柴油发电机,最大限度提升配电网韧性,目标函数如下:

1. 2 约束条件

1. 2. 1 移动储能的时空动态调度约束

        在电网-交通网融合系统中,移动储能的调度状态由其充放电状态和交通运输状态共同决定,具有时空耦合特性。考虑移动储能i在节点j与节点 k 间的交通通行时间与安装配置时间,建立移动储能的时空动态调度模型,设运输过程不消耗电能,其时空动态调度如图2所示。

        由于交通网络容易受到实时路况、突发事件影响,在不同灾害时段与路况影响下移动储能在相同道路的实际通行时间可能存在差异[21]。为了保证灾害场景下应急移动储能设备通行的可靠性,本文考虑冰雪等全覆盖型灾害对交通状况的影响,通过交通融合系数来反映移动储能的实际车速与等效通行距离之间的空间调度关系[17] 。

1. 2. 2 电动汽车充放电约束

        假设电动汽车充电桩与电动汽车均满足“车网互动”条件,即可向电网输电。在灾害发生后,接入充电桩的电动汽车可通过电池放电向负荷提供电能,其充放电功率和荷电状态需要满足以下约束:

1. 2. 3 其他约束

        在灾后恢复阶段,各时刻的负荷削减功率、分布式电源出力与配电网运行仍需要满足相应约束,约束条件形式与灾前预布局阶段相同﹐见式(8)—式(20)。

2 模型求解

        灾后恢复优化模型为混合整数二阶锥规划问题,可利用商业求解器 Gurobi 求解。

3.编程思路分析

3.1参数和变量定义

表1 相关参数

 表2 决策变量

 

3.2编程思路

        根据对文献内容的解读,可以设计下面的编程思路:

步骤1输入所需数据

        部分数据和灾前移动储能预配置阶段的相同,继续用就行其他的数据可在原文提供的附录中获取,这部分主要的问题和难点如下:

        获取各节点的距离矩阵。原文中给出的假设是交通网拓扑与电网相同,相邻电气节点间的道路距离为 2 km。但是,在文章的背景中,由于考虑了配电网重构,配电网拓扑结构可能会不断变化,交通网不可能随着电网拓扑变化。所以假设交通网和初始配电网拓扑相同。另外,需要用到的距离矩阵是所有节点两两之间的距离,目前已知的是相邻节点的距离,需要采用图论算法进行求解。

        步骤2:定义决策变量

        这一步比较简单,按照表2,初始化决策变量即可。要注意的是每个决策变量的维度以及类型(sdpvar还是binvar)不要出错。

        需要注意一点,由于文中共有两个移动储能车,33个节点,11个时段,如果采用常规变量定义方式,alpha_ME将是一个2×33×11的三维变量,这么处理会大大降低程序运行速率。为了方便起见,将定义两个33×11的二维变量,alpha_ME1alpha_ME2,就可以避免使用三维决策变量。

        步骤3:定义目标函数和约束条件,尝试求解确定性优化问题

        按照原文中给定的公式,写出相应的目标函数与约束条件即可。但文中涉及的优化问题是一个多时段的混合整数二阶锥优化,如果采用常规的思维,在写约束条件时经常会用到三层以上的循环,大大拖慢程序运行速度,可以采用一些小技巧避免多层循环的使用,提高程序效率。

        1.以约束条件(26)为例说明如何避免多重循环的使用。式(26)只写了t的取值范围,如果把其他下标的取值范围都标上,式(26)就变成了下面这样:

        按常规思维来看,需要写成5层循环的形式,matlab代码如下:

for t = 1:NTfor i = 1:NMEfor j = 1:NBfor k = 1:NBfor dt = 1:min([floor(T_ME(j,k,t)),NT-t])Constraints = [Constraints , alpha_ME(i,j,t) + alpha_ME(i,k,t+dt) <= 1];endendendend
end

        首先,在定义决策变量时,就已经把两个移动储能的决策变量分开表示,可以减小一层循环,另外,可以将j,k两层循环写成一组数对,方法如下:

j_index = zeros(1,33*33);
k_index = zeros(1,33*33);
it = 1;
for j = 1:33for k = 1:33j_index(it) = j;k_index(it) = k;it = it + 1;end
end

        其中j_indexk_index都是1089(33×33)维的数据,把这两个变量提前存在工作区,使用时读取出来,就可以把5层循环减小到2层循环,如下所示:

for t = 1:NTfor dt = 1:min([floor(T_ME(j_index(k),k_index(k),t)),NT-t])% 式26Constraints = [Constraints , alpha_ME1(j_index,t) + alpha_ME1(k_index,t+dt) <= 1 , alpha_ME2(j_index,t) + alpha_ME2(k_index,t+dt) <= 1];end
end

        经过这样处理,运行时间从原本的0.383611秒减小到0.001126秒,提升的效率达到了99.71%。我只是以约束条件(26)为例进行说明,当优化问题非常复杂时,减少循环语句的使用对程序运行效率的提升是非常惊人的。

        2.另外再强调一点,配电网最优潮流最容易踩坑的地方就是标幺值转换上。之前有朋友拿自己写的代码问我,说感觉公式模型都是按照参考文献打的,但一用求解器就是“Infeasible problem”,拿给我一看,参数中有的数值用的是实际值,有的数值用的是标幺值,非常混乱,我把参数统一修改为标幺值就可以正常运行了。建议在编程时都转换为标幺值求解。

        这一步是最难的部分,需要有耐心,反复调试。为了避免在大段约束条件中寻找错误,建议是每加上一条约束条件,便调试一下模型是否可解,如果可解,再加入下一个约束条件,否则需要反复调试,找到问题所在。

        步骤4:输出计算结果

        原文中将移动储能的预布局和优化调度拆成两阶段优化策略,第二阶段为多源协同的灾后恢复优化模型,是一个混合整数二阶锥规划模型(MISOCP)。这份代码复现了第二阶段多源协同的灾后恢复优化模型的结果,第一阶段的结果可以查看我之前的博客。

4.完整Matlab代码

        完整的matlab代码可以从这个链接获取:

🍞两阶段鲁棒韧性提升策略

5.运行结果分析

5.1表2 移动储能的动态调度结果

 

5.2 图B3 灾害发生后光伏机组的预测出力 

5.3 图4 各时段负荷功率和恢复比例 

5.4 图B4 移动储能有功功率输出与接入位置的关系

5.5 图B5 移动储能荷电状态与接入位置的关系 

5.6 图B6 灾后恢复阶段柴油发电机有功功率输出

 


http://www.ppmy.cn/news/360916.html

相关文章

总线制和多线制示意图_报警主机总线制图纸 总线制和多线制示意图

总线制报警主机和分线主机制报警主机的区别&#xff1f; 支路报警主机是前端检测设备的主机开关触发信号直接连接到主机防区实现报警。 总线报警主机将前端检测设备的开关信号连接到地址模块(或报警模块、防区扩展模块)。当探测装置被触发时&#xff0c;探测器发出报警信号&…

CAN总线转光纤/CAN光端机在海湾消防报警主机联网中应用

来可电子CAN总线转光纤/CAN光端机实现海湾主机远距离传输。也就是先把一台海湾消防主机的CAN总线的电信号通过光纤模块转化为光信号&#xff0c;通过光纤将光信号传输到下一个光纤模块&#xff0c;再把光信号转化为电信号&#xff0c;接到下一个海湾消防主机上&#xff0c;如下…

霍尼236主机说明书_你好 你有 霍尼韦尔Honeywell 236 PLUS 报警主机说明书吗  我想要谢谢你 发给我好吗...

满意答案 suf596zsz 2013.04.22 采纳率&#xff1a;58% 等级&#xff1a;12 已帮助&#xff1a;9154人 粘贴给你 2316Plus 支持十六个防区&#xff0c;内置电平触发输出、可驱动外部继电器&#xff0c;内置电话线检测功能&#xff0c;支持钥匙开关布撤防&#xff0c;采用表面…

北大青鸟消防控制器组网_北大青鸟JBF-61S20防火门监控报警主机使用说明书

一,概述: 北大青鸟环宇消防设备股份有限公司的防火门监控系统由防火门监控器主机和现场接口组成,防火门监控器主机(以下简称防火门监控器)包括三个型号,分别为 JBF-61S20 壁挂机、JBF-61S20-G 柜式机、JBF-61S20-T 台式机。现场接口主要包括 JBF6131-D、JBF6141-D、JBF6142…

CAN光端机在西门子FC18/720消防报警主机远距离联网解决方案

大家在消防火灾报警案例中&#xff0c;经常会碰到组网问题&#xff0c;例如西门子FC18/720消防火灾报警主机而言&#xff0c;主机之间组网是通过CAN总线&#xff0c;而CAN总线一般的通讯距离最长也就是1000米&#xff0c;而且需要在两端加120欧姆电阻才能保证通讯稳定。 如果组…

霍尼236主机说明书_Honeywell CK2300系列报警主机236_238_2316常用编程说明

236 主机编程速成 一、 操作&#xff1a; 布防、撤防&#xff1a;密码&#xff03; 旁路&#xff1a;密码bypass 旁区号 # (bypass 为右下角的键) 取消记忆&#xff1a;&#xff0a;&#xff11;&#xff03; 复位&#xff1a;密码&#xff0a;&#xff16;&#xff18;&#x…

北大青鸟消防控制器组网_北大青鸟消防报警主机维修与调试

一、北大青鸟消防报警主机的特点: 北大青鸟消防报警主机两总线无极性,采用地址编码技术。 整个系统只用两根总线,建筑物布线极其简单,布线路径及方式任意,且不分先后顺序,提高了布线可靠性,也便于穿线施工和线路维修,并可大大降低工程造价。 液晶中文界面,全部采用汉字…

霍尼236主机说明书_霍尼韦尔plus-236报警主机

一、功能简介 1、六防区微电脑控制、键盘编程。 2、编程简单:主叫号、被叫号、地址码、密码、延时报警时间、布防/撤防方式可由用户随意自行操作设定。 3、配置灵活:可跟据不同须求选配进口或国产探头、紧急按钮、脚跳开关、煤气探头、烟雾探测器等。 二、主要功能 1、发生警…

光纤转CAN总线在消防报警主机联网中的应用

火灾是人们不得不防范的一种灾害&#xff0c;它危及人们的生命安全&#xff0c;因此提前做好消防预警工作就十分必要了。在消防报警系统网络中&#xff0c;主要由两部分组成&#xff1a;火灾探测器和火灾报警器。在需要监视火灾的现场安装上火灾探测器&#xff0c;用于监测是否…

Cadence Allegro PCB设计88问解析(二十八) 之 Allegro中dimension environment命令使用(添加及删除尺寸标注)

一个学习信号完整性仿真的layout工程师 最近看到关于Anti Etch的设置&#xff0c;因为本人之前在layout设计是使用过这个命令。后来去到别的公司就不用了&#xff0c;从网上看到说这个命令是用来负片设计的。在这里在说下正片和负片的概念&#xff1a; 正片&#xff1a;是指在a…

JavaWeb笔记(一)

Java网络编程 在JavaSE阶段&#xff0c;我们学习了I/O流&#xff0c;既然I/O流如此强大&#xff0c;那么能否跨越不同的主机进行I/O操作呢&#xff1f;这就要提到Java的网络编程了。 **注意&#xff1a;**本章会涉及到计算机网络相关内容&#xff08;只会讲解大致内容&#x…

GaussDB云数据库SQL应用系列—索引管理

目录 一、前言 二、注意事项 三、索引创建 1、创建普通索引 2、创建唯一索引 3、创建多字段索引 4、创建部分索引 5、创建表达式索引 四、索引管理 1、查看索引信息 2、删除索引 总结 一、前言 随着互联网的快速发展&#xff0c;数据量呈现爆炸式增长。如何高效地…

猪猪猫.CN-WINXPSP2电脑城装机03CD标准版

猪猪猫.CN-WINXPSP2电脑城装机03CD标准版 相对于02标准版的不同之处&#xff1a; 升级系统补丁到2008.03.15 彻底修正无法重新安装OFFICE2003完整版的问题 MD5f8b75ed122c14f6258162a54db34e118 猪猪猫.CN-WINXPSP2电脑城装机03标准版 2008.03 特点&#xff1a;稳定安全、兼…

深雁论坛GhostXP专业装机版 V3.0

深雁论坛GhostXP专业装机版 V3.0 深雁论坛GhostXP专业装机版 V3.0 光盘介绍 ------------------------------------------------------------------------------------------ ----------------------------------- 该版采用微软官方Windows XP SP2免激活版原盘制作而成&#xf…

GhostXP_SP3会员专用版V4.0【IE6+IE8】

系统格式&#xff1a;NTFS 本版系统专为999宝藏网会员打造&#xff0c;快速稳定、极度纯净。仅供论坛技术交流、学习&#xff01; 本系统不带常用软件安装包&#xff0c;不含任何第三方浏览器、无多余插件&#xff01;常用软件请自行下载安装&#xff01; ▂▂▂▂▂▂▂▂▂▂…

【系统之家首发】Ghost_Windows7_sp1_Ultimate_x86V2011.10.10 【OEM 通用版】Windows7旗舰版好人一个出品

Windows7旗舰版Ghost_Windows7_sp1_Ultimate_x86V2011.10.10 【OEM 通用版】&#xff08;好人一个出品&#xff09; 更新说明&#xff1a; 1.系统补丁更新至2011年10月10日 2.Flash控件更新至10.3.183.10 3.更新PE的SRS驱动&#xff0c;及启动 光盘工具&#xff0c;方便装机人员…

Ghost_Windows7_sp1_Ultimate_x86V2011.10.10 【Windows7旗舰版好人一个作品】

Windows7旗舰版Ghost_Windows7_sp1_Ultimate_x86V2011.10.10 【OEM 通用版】&#xff08;好人一个出品&#xff09; 更新说明&#xff1a; 1.系统补丁更新至2011年10月10日 2.Flash控件更新至10.3.183.10 3.更新PE的SRS驱动&#xff0c;及启动 光盘工具&#xff0c;方便装机人员…

11月最新纯净版XP系统999宝藏网GhostXP_SP3会员专用版V5.0

999宝藏网GhostXP_SP3会员专用版V5.0纯净版XP系统 系统格式&#xff1a;NTFS 本版系统专为999宝藏网会员打造&#xff0c;快速稳定、极度纯净。仅供论坛技术交流、学习&#xff01; 本系统不带常用软件安装包&#xff0c;不含任何第三方浏览器、无多余插件&#xff01;常用软…

winnt32简介

U盘安装XP&#xff0c;自动驱动安装&#xff0c;18合一的OEM系统&#xff0c;智能选择 制作通用安装镜像。(在虚拟机上步骤相同&#xff09; 1、硬盘首先分区&#xff0c;再格式化。其实只要分出c区就可以了&#xff0c;其他分区可以在xp下分。 2、用U盘引导电脑进入winpe系统&…

电脑技术员联盟Ghost Xp Sp3纯净版v2.5下载 .

电脑技术员联盟ghost xp Sp3纯净版v2.5 电脑技术员联盟Ghost xp Sp3纯净版v2.5 电脑技术员联盟Ghost Xp Sp3纯净版v2.5 电脑技术员联盟Ghost Xp Sp3纯净版v2.5 技术员联盟Ghost Xp Sp3纯净版v2.5 ■本系统适合笔记本和台式机&#xff08;无分笔记本和台式机&#xff09; ★…