CN116720358A - 面向配电-交通系统韧性多阶段提升的资源优化配置方法 - Google Patents
面向配电-交通系统韧性多阶段提升的资源优化配置方法 Download PDFInfo
- Publication number
- CN116720358A CN116720358A CN202310680065.4A CN202310680065A CN116720358A CN 116720358 A CN116720358 A CN 116720358A CN 202310680065 A CN202310680065 A CN 202310680065A CN 116720358 A CN116720358 A CN 116720358A
- Authority
- CN
- China
- Prior art keywords
- road
- scene
- stage
- typhoon
- traffic
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000005457 optimization Methods 0.000 title claims abstract description 23
- 238000009826 distribution Methods 0.000 claims abstract description 27
- 238000012423 maintenance Methods 0.000 claims description 25
- 230000008439 repair process Effects 0.000 claims description 23
- 230000002787 reinforcement Effects 0.000 claims description 19
- 230000008569 process Effects 0.000 claims description 17
- 230000009977 dual effect Effects 0.000 claims description 15
- 230000009467 reduction Effects 0.000 claims description 15
- 230000008878 coupling Effects 0.000 claims description 13
- 238000010168 coupling process Methods 0.000 claims description 13
- 238000005859 coupling reaction Methods 0.000 claims description 13
- 238000005516 engineering process Methods 0.000 claims description 13
- 230000004044 response Effects 0.000 claims description 12
- 230000003993 interaction Effects 0.000 claims description 11
- 230000003014 reinforcing effect Effects 0.000 claims description 11
- 238000013468 resource allocation Methods 0.000 claims description 11
- 230000006872 improvement Effects 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 9
- 230000009471 action Effects 0.000 claims description 7
- 230000002441 reversible effect Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 230000003190 augmentative effect Effects 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000003416 augmentation Effects 0.000 claims description 4
- 239000000446 fuel Substances 0.000 claims description 4
- 230000000903 blocking effect Effects 0.000 claims description 3
- 230000009194 climbing Effects 0.000 claims description 3
- 239000000295 fuel oil Substances 0.000 claims description 3
- 230000004083 survival effect Effects 0.000 claims description 3
- 238000005728 strengthening Methods 0.000 claims 1
- 102100028862 Delta-aminolevulinic acid dehydratase Human genes 0.000 description 11
- 238000004364 calculation method Methods 0.000 description 10
- 230000007123 defense Effects 0.000 description 5
- 125000004122 cyclic group Chemical group 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000750 progressive effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了面向配电‑交通系统韧性多阶段提升的资源优化配置方法，涉及配电网领域，本发明考虑了台风场景下配电网和交通网的多种交互影响，构建了包含防御资源、应急响应资源和修复资源在内的两阶段随机规划模型，解决了不确定故障背景下耦合系统韧性提升资源的协同优化配置问题，并结合基于罚函数的高斯赛德尔算法和场景求解跳过技术，实现了所提出模型的高效、精确求解。
Description
技术领域
本发明涉及配电网领域，具体是面向配电-交通系统韧性多阶段提升的资源优化配置方法。
背景技术
近年来，台风等频发的自然灾害对以配电网、交通网为代表的城市基础设施系统带来巨大威胁。此外，随着电动汽车和充电站规模的快速扩张，配电网和交通网间的耦合日益密切，进一步使得极端天气下耦合系统级联故障风险急剧增加。因此，研究耦合系统的韧性提升方法对于维持两网安全稳定运行，保障人民生产生活具有重要意义。
配置相应资源以提升系统抵御灾害、应急响应、快速恢复的能力是提高系统韧性的一种有效途径。现有研究多分别从配电网或交通网角度出发，研究一种或某几种韧性提升资源的优化配置方法，忽视了两网之间的交互影响和防御、应急响应、修复多阶段韧性提升资源之间的配合，一定程度上影响了配置效果。且韧性提升资源的优化配置需考虑系统故障的不确定性，相关研究多采用鲁棒优化(Robust Optimization，RO)或随机规划(Stochastic Programming，SP)模型。RO模型考虑最差场景下的最优配置策略，结果相对保守；SP模型以多个场景的运行成本期望值最小为目标制定配置策略，计算负担随场景的增多和系统规模的增大而加重。为减少SP模型的计算量，常将原优化问题基于场景分解为多个子问题并行求解，如采用渐进对冲(Progressive Hedging，PH)算法，但该算法适用于凸优化问题，对于包含整数变量的SP问题存在不收敛的情况。
发明内容
本发明的目的在于提供面向配电-交通系统韧性多阶段提升的资源优化配置方法，通过考虑台风场景下配电网和交通网的多种交互影响和防御、应急响应以及修复资源间的相互配合，实现耦合系统韧性提升资源的协同优化配置，并结合基于罚函数的高斯赛德尔算法和场景求解跳过技术，实现含混合整数的两阶段随机规划问题的高效、精确求解，解决了上述背景技术中提出的问题。
为实现上述目的，本发明提供如下技术方案：
面向配电-交通系统韧性多阶段提升的资源优化配置方法，包括如下步骤：
①构建台风场景下考虑配电-交通耦合系统多阶段韧性提升资源的两阶段随机优化配置模型。多阶段韧性提升资源包括台风过程中的防御性资源(加固电力设备、加固行道树)和台风后的应急响应(搭建联络线、配备交通反向设施)及修复资源(建造仓库和雇佣维修队)。两网交互影响包括：配电网杆塔倒塌会造成附近道路阻塞；配电网故障会影响充电站供电，进而影响电动汽车的出行选择；交通网道路故障会改变电动汽车行驶路线，进而改变充电负荷的时空分布；交通网道路故障会影响维修队的出行，进而影响配电网故障修复情况。其中，第一阶段以总投资成本最小为目标制定资源配置决策，目标函数具体为：
其中，CINV为总投资成本，c1,c2,c3,c4,c5和c6分别为加固电力设备、加固行道树、搭建联络线、配备交通反向设施、建造仓库和雇佣维修队的单位成本，和/> 分别为代表线路加固情况、道路加固情况、联络线建设情况、道路反向资源配置情况、仓库建造情况和维修队雇佣情况的0，1变量。第一阶段中资源配置约束具体为：
其中，NDH,NTH,NTL和NCON为相应资源的上限，CCREW为维修队集合，DDEP为仓库集合。第二阶段以特定台风场景下耦合系统损失最小为目标，目标函数具体为：
其中，第一项是负荷削减成本，第二项是交通拥堵惩罚，第三项是故障持续时间惩罚。cL,cEVC分别为常规负荷削减和电动汽车充电负荷削减的成本系数，cTN,cREP分别为交通拥堵和故障持续时间的惩罚费用系数。pLC,pEVC分别为常规负荷削减功率和电动汽车充电负荷削减功率，xa为道路车流量，ca为道路容量，μij,t为表示线路运行状态的0，1变量，εa为避免分母为零的小常数。第二阶段中，配电网运行约束具体为：
其中，EB,EL分别为配电网节点和线路集合，A(j),B(j)分别为以节点j为起点和终点的线路的集合，P,Q分别为线路传输的有功和无功功率，pG,qG分别发电机输出的有功和无功功率，pL,qL分别为常规有功和无功负荷，qLC为无功负荷削减功率，pEVC为电动汽车充电负荷，U为节点电压幅值的平方，r,x分别为线路电阻和电抗。Pmax,Qmax,Umax,/> Umin分别为相应变量的限值，pramp为发电机最大爬坡限值，Mij为大常数。交通网运行约束具体为：
其中，TA为道路集合，分别为电动汽车和燃油车行驶路径集合，C(j)为节点j连接的充电站集合。δa,p为表示路径与道路的关系的0，1变量，当路径p经过道路a时，该变量为1，否则为0。/>分别为选择路径p的电动汽车和燃油车流量，/>别为电动汽车和燃油车出行需求。PEV为电动汽车充电功率。台风中防御资源动作约束具体为：
其中，Na为建设在道路附近的电力线路集合，其倒塌有可能影响道路通行。μij,0,μa,0分别为电力线路和交通道路在台风中的生存情况。分别为采取和不采取电力设备加固措施下某一台风场景s中线路的故障情况，其值在场景生成过程中由采样得到。考虑两网交互，电网杆塔倒塌有可能造成交通网阻塞，/>分别为同时采取电力设备和行道树加固，只有行道树加固，只有电力设备加固和未采取加固措施下道路故障情况，其值在场景生成过程中由采样得到。台风后应急响应资源动作约束具体为：
其中，ETL为联络线集合。λa为表示道路反向情况的0，1变量，a'为与道路a方向相反的道路。台风后修复资源动作约束具体为：
其中，yc,ij为表示维修队c是否前去维修线路(i,j)的0，1变量，pathc,ij,a为表示前往维修线路(i,j)的维修队c选择的行驶路径是否经过道路a的0，1变量，pathc,ij为前往维修线路(i,j)的维修队c选择的行驶路径。为维修队到达故障地点的时间，Ta为道路通行时间，/> 为维修队修复故障所用时间，△t为每个时间间隔的长度。/>为故障被修复的时刻，zij,t为故障的修复状态。Anode×road为交通网节点与道路关系矩阵，ADEP×node,Aline×node分别为仓库和线路地理位置矩阵，Anode×node为单位阵。两阶段随机优化配置模型目标函数具体为：
其中，ps为每个台风场景出现的概率，S为台风场景集合。两阶段随机优化配置模型的约束包含第一阶段与第二阶段的所有约束，台风场景集由场景生成方法获得。
②所述场景生成方法叙述如下：确定台风相关参数、登陆坐标和移动路径；根据台风风速经验模型，计算台风过程中电网、交通网各设备所经历的最大风速。所述风速经验模型具体为：
其中，v为风速，xe为设备与风眼之间的距离，vm为最大风速，Rm为最大风速半径，R为台风影响范围的最大半径,此处风速为vm/κ，K和κ均为可调参数。根据设备经受的最大风速和脆弱性曲线，确定设备加固和不加固情况下的故障概率，采用均匀采样法，确定各台风场景下线路和道路故障状态
③结合基于罚函数的高斯赛德尔算法和场景求解跳过技术，求解两阶段随机优化配置模型。所述基于罚函数的高斯赛德尔(Penalty-based Gauss Seidel，PBGS)算法具体为:
将第一阶段变量x按场景复制为xs，并施加约束：
其中，z为辅助变量。将该约束进行拉格朗日松弛，获得原问题的对偶问题，并引入正则惩罚项保证两者间的对偶间隙为零。引入惩罚项后的增广拉格朗日对偶问题具体为：
其中，ys为第二阶段变量，为第一阶段变量可行域，/>为第二阶段变量可行域。ρs为惩罚系数，[f]-＝-min{0,f}，β为大于1的参数，u为求解过程中的外层循环迭代次数。增广拉格朗日对偶问题的解即为原问题的解，设置外层循环更新惩罚系数，更新规则具体为：
其中，γs,为给定参数。外层循环直至所有场景第一阶段变量/>满足小于等于设定误差或循环达到最大设定次数时中断。设置内层循环，在惩罚系数固定的情况下，将增广拉格朗日对偶问题按场景分解为子问题进行并行求解，即并行求解内层循环采用高斯赛德尔算法，交替更新(xs,ys)和z，直至/> 小于等于设定误差或循环达到最大设定次数时中断，l为内层循环迭代次数。PBGS算法结合场景求解跳过技术，可在循环过程中避免部分已求得最优解的场景子问题的重复求解，明显减少了计算时间，场景求解跳过技术具体为：当满足/>时，子问题在第u次循环中的求解可以跳过，其中，lu-1为第(u-1)次外循环中的内循环次数。
与现有技术相比，本发明的有益效果是：本发明通过考虑配电网和交通网的多种交互影响和防御、应急响应以及修复资源间的相互配合，实现耦合系统韧性提升资源的协同优化配置，并结合基于罚函数的高斯赛德尔算法和场景求解跳过技术，实现含混合整数的两阶段随机规划问题的高效、精确求解。
附图说明
图1为本发明方法流程图。
图2为本发明实施例中33节点配电网系统示意图。
图3为本发明实施例中12节点交通网系统示意图。
图4为本发明实施例中所提出的算法与PH算法在不同场景集规模下的收敛曲线图。
图5为本发明实施例中所提出的算法与PBGS算法实际求解场景数示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参阅图1，本发明实施例中，①构建台风场景下考虑配电-交通耦合系统多阶段韧性提升资源的两阶段随机优化配置模型。多阶段韧性提升资源包括台风过程中的防御性资源(加固电力设备、加固行道树)和台风后的应急响应(搭建联络线、配备交通反向设施)及修复资源(建造仓库和雇佣维修队)。两网交互影响包括：配电网杆塔倒塌会造成附近道路阻塞；配电网故障会影响充电站供电，进而影响电动汽车的出行选择；交通网道路故障会改变电动汽车行驶路线，进而改变充电负荷的时空分布；交通网道路故障会影响维修队的出行，进而影响配电网故障修复情况。其中，第一阶段以总投资成本最小为目标制定资源配置决策，目标函数具体为：
其中，CINV为总投资成本，c1,c2,c3,c4,c5和c6分别为加固电力设备、加固行道树、搭建联络线、配备交通反向设施、建造仓库和雇佣维修队的单位成本，和/> 分别为代表线路加固情况、道路加固情况、联络线建设情况、道路反向资源配置情况、仓库建造情况和维修队雇佣情况的0，1变量。第一阶段中资源配置约束具体为：
其中，NDH,NTH,NTL和NCON为相应资源的上限，CCREW为维修队集合，DDEP为仓库集合。第二阶段以特定台风场景下耦合系统损失最小为目标，目标函数具体为：
其中，第一项是负荷削减成本，第二项是交通拥堵惩罚，第三项是故障持续时间惩罚。cL,cEVC分别为常规负荷削减和电动汽车充电负荷削减的成本系数，cTN,cREP分别为交通拥堵和故障持续时间的惩罚费用系数。pLC,pEVC分别为常规负荷削减功率和电动汽车充电负荷削减功率，xa为道路车流量，ca为道路容量，μij,t为表示线路运行状态的0，1变量，εa为避免分母为零的小常数。第二阶段中，配电网运行约束具体为：
其中，EB,EL分别为配电网节点和线路集合，A(j),B(j)分别为以节点j为起点和终点的线路的集合，P,Q分别为线路传输的有功和无功功率，pG,qG分别发电机输出的有功和无功功率，pL,qL分别为常规有功和无功负荷，qLC为无功负荷削减功率，pEVC为电动汽车充电负荷，U为节点电压幅值的平方，r,x分别为线路电阻和电抗。Pmax,Qmax,Umax,/> Umin分别为相应变量的限值，pramp为发电机最大爬坡限值，Mij为大常数。交通网运行约束具体为：
其中，TA为道路集合，分别为电动汽车和燃油车行驶路径集合，C(j)为节点j连接的充电站集合。δa,p为表示路径与道路的关系的0，1变量，当路径p经过道路a时，该变量为1，否则为0。/>分别为选择路径p的电动汽车和燃油车流量，/>别为电动汽车和燃油车出行需求。PEV为电动汽车充电功率。台风中防御资源动作约束具体为：
其中，Na为建设在道路附近的电力线路集合，其倒塌有可能影响道路通行。μij,0,μa,0分别为电力线路和交通道路在台风中的生存情况。分别为采取和不采取电力设备加固措施下某一台风场景s中线路的故障情况，其值在场景生成过程中由采样得到。考虑两网交互，电网杆塔倒塌有可能造成交通网阻塞，/>分别为同时采取电力设备和行道树加固，只有行道树加固，只有电力设备加固和未采取加固措施下道路故障情况，其值在场景生成过程中由采样得到。台风后应急响应资源动作约束具体为：
其中，ETL为联络线集合。λa为表示道路反向情况的0，1变量，a'为与道路a方向相反的道路。台风后修复资源动作约束具体为：
其中，yc,ij为表示维修队c是否前去维修线路(i,j)的0，1变量，pathc,ij,a为表示前往维修线路(i,j)的维修队c选择的行驶路径是否经过道路a的0，1变量，pathc,ij为前往维修线路(i,j)的维修队c选择的行驶路径。为维修队到达故障地点的时间，Ta为道路通行时间，/> 为维修队修复故障所用时间，△t为每个时间间隔的长度。/>为故障被修复的时刻，zij,t为故障的修复状态。Anode×road为交通网节点与道路关系矩阵，ADEP×node,Aline×node分别为仓库和线路地理位置矩阵，Anode×node为单位阵。两阶段随机优化配置模型目标函数具体为：
其中，ps为每个台风场景出现的概率，S为台风场景集合。两阶段随机优化配置模型的约束包含第一阶段与第二阶段的所有约束，台风场景集由场景生成方法获得。
②所述场景生成方法叙述如下：确定台风相关参数、登陆坐标和移动路径；根据台风风速经验模型，计算台风过程中电网、交通网各设备所经历的最大风速。所述风速经验模型具体为：
其中，v为风速，xe为设备与风眼之间的距离，vm为最大风速，Rm为最大风速半径，R为台风影响范围的最大半径,此处风速为vm/κ，K和κ均为可调参数。根据设备经受的最大风速和脆弱性曲线，确定设备加固和不加固情况下的故障概率，采用均匀采样法，确定各台风场景下线路和道路故障状态
③结合基于罚函数的高斯赛德尔算法和场景求解跳过技术，求解两阶段随机优化配置模型。所述基于罚函数的高斯赛德尔(Penalty-based Gauss Seidel，PBGS)算法具体为:
将第一阶段变量x按场景复制为xs，并施加约束：
其中，z为辅助变量。将该约束进行拉格朗日松弛，获得原问题的对偶问题，并引入正则惩罚项保证两者间的对偶间隙为零。引入惩罚项后的增广拉格朗日对偶问题具体为：
其中，ys为第二阶段变量，为第一阶段变量可行域，/>为第二阶段变量可行域。ρs为惩罚系数，[f]-＝-min{0,f}，β为大于1的参数，u为求解过程中的外层循环迭代次数。增广拉格朗日对偶问题的解即为原问题的解，设置外层循环更新惩罚系数，更新规则具体为：
其中，γs,为给定参数。外层循环直至所有场景第一阶段变量/>满足小于等于设定误差或循环达到最大设定次数时中断。设置内层循环，在惩罚系数固定的情况下，将增广拉格朗日对偶问题按场景分解为子问题进行并行求解，即并行求解内层循环采用高斯赛德尔算法，交替更新(xs,ys)和z，直至/> 小于等于设定误差或循环达到最大设定次数时中断，l为内层循环迭代次数。PBGS算法结合场景求解跳过技术，可在循环过程中避免部分已求得最优解的场景子问题的重复求解，明显减少了计算时间，场景求解跳过技术具体为：当满足/>时，子问题在第u次循环中的求解可以跳过，其中，lu-1为第(u-1)次外循环中的内循环次数。
作为本发明进一步的实施例，请参阅图2和图3，本实施例涉及IEEE 33节点的配电系统和12节点环状交通网系。耦合系统中共有8个快充电站，各韧性资源候选安装位置及数量限制如表1所示：
表1各韧性资源配置限制
模型中各成本系数设置如表2所示：
表2成本系数设置
33节点配电网系统的基准电压为12.66kV，基准容量为10MVA，节点电压幅值的安全范围为0.89p.u.～1.05p.u.。电动汽车充电功率为50kW。台风登陆坐标为(-150km,-120km)，最大风速半径，最大台风影响半径和最大风速分别为25km,300km和70m/s。求解过程中，内外层循环收敛判据设置为0.01，初始惩罚系数ρ0设置为1500，γs和与ρ0相等，β设置为1.7。
表3给出了不同场景集规模下，直接求解两阶段随机规划模型、采用PH算法、采用PBGS算法和采用本实施例提出的结合场景跳过技术的PBGS算法的计算结果和计算时间。当场景集规模为20时，直接求解法无法在60000s计算时间限制内获得最优结果；其他情况下，本实施例提出的算法求解结果与直接求解结果间误差低于1.1％，验证了本实施例提出的求解方法的精度。
图4给出了不同场景集规模下，本实施例提出的结合场景跳过技术的PBGS算法与常用的PH算法的收敛曲线，最大循环次数设置为30。所有场景集规模下，本实施例提出的算法均能收敛，PH算法均不能收敛。结合表3，相比于PH算法，本实施提出的算法计算时间平均减少57.2％，且计算结果更接近于直接求解方法，说明本实施提出的算法在计算精度和效率上均优于PH算法。
图5给出不同场景规模下，采用PBGS算法和本实施例提出的算法的实际求解场景数。随着场景集规模的增大，本实施例提出的求解算法求解过程中跳过的场景数随之增加，结合表3，其求解速度相较于不采用场景跳过技术的PBGS算法平均提高了7.6％，与PBGS算法结果间的误差均小于0.3％，说明本实施例提出的算法中结合的场景跳过技术在保证计算精度的同时进一步提升了计算效率。
表3不同求解方法效果对比
接下来，设置6个对比案例以验证本实施例所提优化配置策略的有效性。案例1：不配置韧性提升资源；案例2：只配置交通网侧的韧性提升资源，即加固行道树和配备道路反向资源；案例3：只配置配电网侧的韧性提升资源，即加固电力设备、安装联络线、建设仓库和雇佣维修队；案例4：同时配置两网韧性提升资源，但不考虑两网之间的交互影响；案例5：采用本实施例提出的耦合系统多阶段韧性提升资源的优化配置方法；案例6：不考虑资源间的相互配合，防御资源、应急响应资源和修复资源独立优化配置。案例1-5第一阶段总投资成本，第二阶段负荷削减成本、交通拥堵惩罚和故障持续时间惩罚，投资收益率如表4所示。可以看出，本实施例所提优化配置方法的投资收益率最高，说明了在所提方法中同时考虑配电网和交通网侧资源及两网交互的有效性。表5给出了案例5、6中各类资源的投资成本和所有场景下第二阶段目标函数的期望值，可以看出本实施例所提优化配置方法中考虑了多阶段韧性提升资源的配合，避免了过度投资和低效投资。
表4案例1-5结果对比
表5案例5、6结果对比
以上算法中，为本发明独创、从未被公开且其工作方式与任何现有文献记载均不相同的是：本发明通过考虑配电网和交通网的多种交互影响和防御、应急响应以及修复资源间的相互配合，实现耦合系统韧性提升资源的协同优化配置，并结合基于罚函数的高斯赛德尔算法和场景求解跳过技术，实现含混合整数的两阶段随机规划问题的高效、精确求解。
对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。
Claims (1)
1.面向配电-交通系统韧性多阶段提升的资源优化配置方法，其特征在于，包括如下步骤：
①构建台风场景下考虑配电-交通耦合系统多阶段韧性提升资源的两阶段随机优化配置模型；
多阶段韧性提升资源包括台风过程中的防御性资源和台风后的应急响应及修复资源，两网交互影响包括：配电网杆塔倒塌造成附近道路阻塞；配电网故障影响充电站供电，进而影响电动汽车的出行选择；交通网道路故障改变电动汽车行驶路线，进而改变充电负荷的时空分布；交通网道路故障影响维修队的出行，进而影响配电网故障修复情况，其中，第一阶段以总投资成本最小为目标制定资源配置决策，目标函数具体为：
其中，CINV为总投资成本，c1,c2,c3,c4,c5和c6分别为加固电力设备、加固行道树、搭建联络线、配备交通反向设施、建造仓库和雇佣维修队的单位成本，和/> 分别为代表线路加固情况、道路加固情况、联络线建设情况、道路反向资源配置情况、仓库建造情况和维修队雇佣情况的0，1变量，第一阶段中资源配置约束具体为：
其中，NDH,NTH,NTL和NCON为相应资源的上限，CCREW为维修队集合，DDEP为仓库集合，第二阶段以特定台风场景下耦合系统损失最小为目标，目标函数具体为：
其中，第一项是负荷削减成本，第二项是交通拥堵惩罚，第三项是故障持续时间惩罚，cL,cEVC分别为常规负荷削减和电动汽车充电负荷削减的成本系数，cTN,cREP分别为交通拥堵和故障持续时间的惩罚费用系数，pLC,pEVC分别为常规负荷削减功率和电动汽车充电负荷削减功率，xa为道路车流量，ca为道路容量，μij,t为表示线路运行状态的0，1变量，εa为避免分母为零的小常数，第二阶段中，配电网运行约束具体为：
其中，EB,EL分别为配电网节点和线路集合，A(j),B(j)分别为以节点j为起点和终点的线路的集合，P,Q分别为线路传输的有功和无功功率，pG,qG分别发电机输出的有功和无功功率，pL,qL分别为常规有功和无功负荷，qLC为无功负荷削减功率，pEVC为电动汽车充电负荷，U为节点电压幅值的平方，r,x分别为线路电阻和电抗，Pmax,Qmax,Umax,/> Umin分别为相应变量的限值，pramp为发电机最大爬坡限值，Mij为大常数，交通网运行约束具体为：
其中，TA为道路集合，分别为电动汽车和燃油车行驶路径集合，C(j)为节点j连接的充电站集合，δa,p为表示路径与道路的关系的0，1变量，当路径p经过道路a时，该变量为1，否则为0，/>分别为选择路径p的电动汽车和燃油车流量，/>别为电动汽车和燃油车出行需求，PEV为电动汽车充电功率，台风中防御资源动作约束具体为：
其中，Na为建设在道路附近的电力线路集合，μij,0,μa,0分别为电力线路和交通道路在台风中的生存情况，分别为采取和不采取电力设备加固措施下某一台风场景s中线路的故障情况，其值在场景生成过程中由采样得到，/>分别为同时采取电力设备和行道树加固，只有行道树加固，只有电力设备加固和未采取加固措施下道路故障情况，其值在场景生成过程中由采样得到，台风后应急响应资源动作约束具体为：
其中，ETL为联络线集合，λa为表示道路反向情况的0，1变量，a'为与道路a方向相反的道路，台风后修复资源动作约束具体为：
其中，yc,ij为表示维修队c是否前去维修线路(i,j)的0，1变量，pathc,ij,a为表示前往维修线路(i,j)的维修队c选择的行驶路径是否经过道路a的0，1变量，pathc,ij为前往维修线路(i,j)的维修队c选择的行驶路径，为维修队到达故障地点的时间，Ta为道路通行时间，/> 为维修队修复故障所用时间，△t为每个时间间隔的长度，/>为故障被修复的时刻，zij,t为故障的修复状态，Anode×road为交通网节点与道路关系矩阵，ADEP×node,Aline×node分别为仓库和线路地理位置矩阵，Anode×node为单位阵，两阶段随机优化配置模型目标函数具体为：
其中，ps为每个台风场景出现的概率，S为台风场景集合，两阶段随机优化配置模型的约束包含第一阶段与第二阶段的所有约束，台风场景集由场景生成方法获得；
②所述场景生成方法叙述如下：确定台风相关参数、登陆坐标和移动路径；根据台风风速经验模型，计算台风过程中电网、交通网各设备所经历的最大风速，所述风速经验模型具体为：
其中，v为风速，xe为设备与风眼之间的距离，vm为最大风速，Rm为最大风速半径，R为台风影响范围的最大半径,此处风速为vm/κ，K和κ均为可调参数，根据设备经受的最大风速和脆弱性曲线，确定设备加固和不加固情况下的故障概率，采用均匀采样法，确定各台风场景下线路和道路故障状态
③结合基于罚函数的高斯赛德尔算法和场景求解跳过技术，求解两阶段随机优化配置模型，所述基于罚函数的高斯赛德尔算法具体为：
将第一阶段变量x按场景复制为xs，并施加约束：
其中，z为辅助变量，将该约束进行拉格朗日松弛，获得原问题的对偶问题，并引入正则惩罚项保证两者间的对偶间隙为零，引入惩罚项后的增广拉格朗日对偶问题具体为：
其中，ys为第二阶段变量，为第一阶段变量可行域，/>为第二阶段变量可行域，ρs为惩罚系数，[f]-＝-min{0,f}，β为大于1的参数，u为求解过程中的外层循环迭代次数，增广拉格朗日对偶问题的解即为原问题的解，设置外层循环更新惩罚系数，更新规则具体为：
其中，γs,为给定参数，外层循环直至所有场景第一阶段变量/>满足/>小于等于设定误差或循环达到最大设定次数时中断，设置内层循环，在惩罚系数固定的情况下，将增广拉格朗日对偶问题按场景分解为子问题进行并行求解，即并行求解/>内层循环采用高斯赛德尔算法，交替更新(xs,ys)和z，直至/> 小于等于设定误差或循环达到最大设定次数时中断，l为内层循环迭代次数，场景求解跳过技术具体为：当满足/>时，子问题在第u次循环中的求解可以跳过，其中，lu-1为第(u-1)次外循环中的内循环次数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310680065.4A CN116720358B (zh) | 2023-06-09 | 2023-06-09 | 面向配电-交通系统韧性多阶段提升的资源优化配置方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310680065.4A CN116720358B (zh) | 2023-06-09 | 2023-06-09 | 面向配电-交通系统韧性多阶段提升的资源优化配置方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116720358A true CN116720358A (zh) | 2023-09-08 |
| CN116720358B CN116720358B (zh) | 2024-02-02 |
Family
ID=87874595
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310680065.4A Active CN116720358B (zh) | 2023-06-09 | 2023-06-09 | 面向配电-交通系统韧性多阶段提升的资源优化配置方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116720358B (zh) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160335223A1 (en) * | 2014-06-27 | 2016-11-17 | University Of South Florida | Methods and systems for computation of bilevel mixed integer programming problems |
| CN112132427A (zh) * | 2020-09-10 | 2020-12-25 | 国家电网有限公司 | 一种考虑用户侧多种资源接入的电网多层规划方法 |
| CN113742917A (zh) * | 2021-09-06 | 2021-12-03 | 东南大学 | 一种考虑多阶段恢复过程的综合能源系统韧性提升方法 |
| CN114218753A (zh) * | 2021-11-22 | 2022-03-22 | 华南理工大学 | 台风灾害下配电网弹性评估系统及其二阶段弹性提升方法 |
| CN115130745A (zh) * | 2022-06-24 | 2022-09-30 | 国网河南省电力公司开封供电公司 | 变频空调负荷参与输配电网备用分布式协调优化配置方法 |
| CN115663910A (zh) * | 2022-07-20 | 2023-01-31 | 东南大学 | 配电网规划层-运行层两阶段分布鲁棒优化治理方法 |
| CN115688963A (zh) * | 2022-08-01 | 2023-02-03 | 杭州电子科技大学 | 考虑源网荷不确定性的新型电力系统输配协同规划方法 |
| CN115765015A (zh) * | 2022-11-30 | 2023-03-07 | 国网湖南省电力有限公司 | 面向电网实际应用场景的源网荷储协同互动方案制定方法 |
| CN115912330A (zh) * | 2022-09-27 | 2023-04-04 | 国网湖北省电力有限公司电力科学研究院 | 一种有源配电网拓扑演化模型的两阶段机会约束优化方法和系统 |
| CN115995790A (zh) * | 2022-11-25 | 2023-04-21 | 华中科技大学 | 一种配电网故障恢复方法、系统、设备及介质 |
| CN116090840A (zh) * | 2022-10-11 | 2023-05-09 | 云南电网有限责任公司电力科学研究院 | 基于储能规划的配电网韧性提升方法、装置、设备和介质 |
| CN116090677A (zh) * | 2023-04-10 | 2023-05-09 | 湖南大学 | 考虑电力-通信-交通网络耦合的空地应急资源规划方法 |
-
2023
- 2023-06-09 CN CN202310680065.4A patent/CN116720358B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20160335223A1 (en) * | 2014-06-27 | 2016-11-17 | University Of South Florida | Methods and systems for computation of bilevel mixed integer programming problems |
| CN112132427A (zh) * | 2020-09-10 | 2020-12-25 | 国家电网有限公司 | 一种考虑用户侧多种资源接入的电网多层规划方法 |
| CN113742917A (zh) * | 2021-09-06 | 2021-12-03 | 东南大学 | 一种考虑多阶段恢复过程的综合能源系统韧性提升方法 |
| CN114218753A (zh) * | 2021-11-22 | 2022-03-22 | 华南理工大学 | 台风灾害下配电网弹性评估系统及其二阶段弹性提升方法 |
| CN115130745A (zh) * | 2022-06-24 | 2022-09-30 | 国网河南省电力公司开封供电公司 | 变频空调负荷参与输配电网备用分布式协调优化配置方法 |
| CN115663910A (zh) * | 2022-07-20 | 2023-01-31 | 东南大学 | 配电网规划层-运行层两阶段分布鲁棒优化治理方法 |
| CN115688963A (zh) * | 2022-08-01 | 2023-02-03 | 杭州电子科技大学 | 考虑源网荷不确定性的新型电力系统输配协同规划方法 |
| CN115912330A (zh) * | 2022-09-27 | 2023-04-04 | 国网湖北省电力有限公司电力科学研究院 | 一种有源配电网拓扑演化模型的两阶段机会约束优化方法和系统 |
| CN116090840A (zh) * | 2022-10-11 | 2023-05-09 | 云南电网有限责任公司电力科学研究院 | 基于储能规划的配电网韧性提升方法、装置、设备和介质 |
| CN115995790A (zh) * | 2022-11-25 | 2023-04-21 | 华中科技大学 | 一种配电网故障恢复方法、系统、设备及介质 |
| CN115765015A (zh) * | 2022-11-30 | 2023-03-07 | 国网湖南省电力有限公司 | 面向电网实际应用场景的源网荷储协同互动方案制定方法 |
| CN116090677A (zh) * | 2023-04-10 | 2023-05-09 | 湖南大学 | 考虑电力-通信-交通网络耦合的空地应急资源规划方法 |
Non-Patent Citations (3)
| Title |
|---|
| 孔顺飞;胡志坚;谢仕炜;杨黎;郑云飞;: "含电动汽车充电站的主动配电网二阶段鲁棒规划模型及其求解方法", 电工技术学报, no. 05 * |
| 杨火明;徐潇源;严正;: "考虑配电网韧性的储能系统选址定容优化方法", 电力建设, no. 01 * |
| 赵磊;杨莉;: "计及复杂场景运行的主动配电网多阶段规划", 能源工程, no. 04 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116720358B (zh) | 2024-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112751350B (zh) | 一种移动储能时空联合优化调度策略的制定方法 | |
| CN110460036B (zh) | 一种考虑风电不确定性的交直流配电网分布式优化方法 | |
| CN111555280B (zh) | 基于电-气综合能源系统的弹性配电网灾后恢复控制方法 | |
| CN103699950A (zh) | 一种考虑交通网络流量的电动汽车充电站规划方法 | |
| CN105740556A (zh) | 基于客流需求的列车运行图自动编制方法 | |
| CN104600713A (zh) | 含风力/光伏发电配电网日前无功调度的生成装置及方法 | |
| CN108649605A (zh) | 一种基于双层场景区间潮流的der并网准入容量规划方法 | |
| CN109740975B (zh) | 考虑潮汐资源评估与流速方向的潮汐流能发电场规划方法 | |
| CN104156788A (zh) | 一种基于禁忌搜索算法的配网资源抢修优化调度方法 | |
| CN115995790A (zh) | 一种配电网故障恢复方法、系统、设备及介质 | |
| CN115310378A (zh) | 一种极端台风灾害下电网韧性评估及差异化规划方法 | |
| CN115459256A (zh) | 基于配电网-交通网协同的配网负荷恢复方法及系统 | |
| CN115439000A (zh) | 考虑风光荷功率不确定性和相关性的配电网区块划分方法 | |
| CN114928110A (zh) | 基于n-k故障的规划调度协同优化的三层鲁棒优化方法 | |
| Li et al. | An MILP method for optimal planning of electric vehicle charging stations in coordinated urban power and transportation networks | |
| CN116720358B (zh) | 面向配电-交通系统韧性多阶段提升的资源优化配置方法 | |
| CN116961057A (zh) | 一种考虑电动汽车的多时段配电网故障恢复方法 | |
| CN113987729B (zh) | 一种基于土地利用的大规模城市轨道网络自动生成方法 | |
| CN115102192A (zh) | 一种应对极端天气的配电网弹性控制方法 | |
| CN116415708A (zh) | 考虑置信水平的电网鲁棒规划方法 | |
| CN113762588A (zh) | 以城市路网为主体的充电基础设施配置方法 | |
| CN114239206A (zh) | 基于强化学习的灾后电-路协同修复方法、系统与设备 | |
| CN104572438A (zh) | 多变电站设备运行维护优化方法 | |
| Du et al. | Network optimization of conventional public transit based on urban rail transit | |
| Luo et al. | Path and transport mode selection in multimodal transportation with time window |
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 |