CN115099702A - 基于拉格朗日松弛算法的电动公交日间运行充电优化方法 - Google Patents
基于拉格朗日松弛算法的电动公交日间运行充电优化方法 Download PDFInfo
- Publication number
- CN115099702A CN115099702A CN202210874670.0A CN202210874670A CN115099702A CN 115099702 A CN115099702 A CN 115099702A CN 202210874670 A CN202210874670 A CN 202210874670A CN 115099702 A CN115099702 A CN 115099702A
- Authority
- CN
- China
- Prior art keywords
- station
- vehicle
- charging
- line
- bus
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000009977 dual effect Effects 0.000 claims description 12
- 230000008878 coupling Effects 0.000 claims description 10
- 238000010168 coupling process Methods 0.000 claims description 10
- 238000005859 coupling reaction Methods 0.000 claims description 10
- 230000000295 complement effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 238000011478 gradient descent method Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 description 5
- 238000005265 energy consumption Methods 0.000 description 4
- 239000000446 fuel Substances 0.000 description 3
- 150000001875 compounds Chemical class 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
- B60L53/00—Methods of charging batteries, specially adapted for electric vehicles; Charging stations or on-board charging equipment therefor; Exchange of energy storage elements in electric vehicles
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Power Engineering (AREA)
- Biodiversity & Conservation Biology (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Water Supply & Treatment (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
Abstract
一种基于拉格朗日松弛算法的电动公交日间运行充电优化方法,在保证车辆预定时刻表的前提下,以最小化***总充电时间为目标,基于离散电池电量状态获得全线路网络级充电方案整数优化原问题模型,然后转换为拉格朗日松弛问题,接着分解为多个子问题,最后转化为多个紧凑子问题,通过求解多个紧凑子问题的最优解,获得还原后的原问题模型P的最优解,精准描述大规模公交网络的充电问题,一步步减少模型中的决策变量和约束条件数量,缩减问题规模,提升求解效率,进行快速高效的求解。
Description
技术领域
本发明涉及一种基于拉格朗日松弛算法的电动公交日间运行充电优化方法。
背景技术
城市化和机动化导致车辆拥有量急剧上升,造成了大量环境(如噪音和污染物排放)和社会(如能源消耗、交通拥堵)问题,交通领域产生的能源消耗已达到全球能源消耗总量的25%以上,温室气体排放量占比27%,燃油车的使用被认为是对全球变暖的严重威胁。
电动车作为燃油车的替代品,对减少城市交通对气候的负面影响和改善空气质量条件有巨大帮助。过去数十年来,全球许多城市,都在城市交通电气化方面提供大量持续投资,尤其在公共交通电气化方面取得显著成效。
而电动公交车相比于燃油公交车的最大不足在于,有限的电池容量无法满足全天的运营能耗需求,需要在日间运行过程中补电;电池充电时间与电量状态之间的非线性关系对于掌控公交行程电量带来挑战。
而现有研究中,缺乏将电动公交车的非线性充电过程与行程途中的电量状态跟踪进行融合考虑,大规模网络中公交运营时刻表与充电时序的协调尚未形成快速有效的求解办法。
发明内容
本发明要解决的技术问题是提供一种基于拉格朗日松弛算法的电动公交日间运行充电优化策略,根据车辆运营时刻表与路网信息,通过刻画电池充电的非线性过程,跟踪记录车辆到离各站点的电量状态,建立基于拉格朗日松弛算法的电动公交日间运行充电优化模型,并设计了高效的求解算法。
为了解决上述技术问题,本发明采用的技术方案是:一种基于拉格朗日松弛算法的电动公交日间运行充电优化方法,包括以下具体步骤:S1、获取电动公交数据和包含线路数量的线路运行数据,构建基于离散电池电量状态的全线路网络级充电方案整数优化模型,获得包含耦合约束的原问题模型P;S2、通过引入拉格朗日乘子将所述耦合约束转换到原问题模型P的目标函数中,从而获得具有新的约束条件的拉格朗日松弛问题P’(μ,λ);S3、根据所述线路数量,将所述拉格朗日松弛问题P’(μ,λ)分解为多个子问题Pl’,引入新的决策变量,即相应线路的车辆是否在相应站点充电的决策变量,从而将多个所述子问题Pl’转化为多个紧凑子问题Pl”;S4、通过求解多个紧凑子问题Pl”,求解出P’(μ,λ)的目标函数,即可求解出原问题模型P的目标函数,将求解原问题模型P的下界的最优解转化为求解拉格朗日对偶问题PLD,利用次梯度下降法求解拉格朗日对偶问题PLD,并在每次迭代求解后根据当前解更新拉格朗日乘子,采用k-greedy算法计算原问题模型P的上界的最优解,直至原问题模型P的上界的最优解等于下界的最优解,算法达到终止条件,获得原问题模型P的最优解,即获得各电动公交车辆日间运行充电方案,完成对电动公交日间运行充电的优化。
在某些实施方式中,所述电动公交数据包括每辆公交车在各站点的最大充电次数、车辆最低电量保护值、电动公交电量状态,所述线路运行数据包括途径站点、各站点的到发时间表、各站点的充电桩数目。
s.t.
其中,
s.t.表示约束条件为;
N和L分别表示公交站点位置和公交线路数量,l表示第l条公交线路,l∈L,n表示第n个公交站点,n∈N,Nl表示线路l途径的车站总数;
al,j表示线路l的车辆在第j个站点预定的到站时间,dl,j表示线路l的车辆在第j个站点预定的离站时间,j∈{1,…,Nl-1},将规划时段T离散为以△tmin为单位的间隔,t,t’,t”均为规划时段T内的某一时间点,充电方案由决策变量ul,j tt’或ul,j t’t”表示,当车辆在时间间隔[t,t’]或[t’,t”]时正于线路l的第j个站点充电时ul,j tt’=1或ul,j t’t”=1;
S为离散集合,S表示电动公交电量状态SOC,S={0,1,2,…,Q},集合S中的每个值s∈S代表车辆从SOC=0状态开始充电s个时间间隔△t后所达到的电量状态,即s=0表示车辆充电△t时长后电池达到的电量,s=2表示车辆充电2△t时长后电池达到的电量,令电池从SOC=0的状态充至满电所需的时间为H小时,Q表示电量状态Q,即满电状态,则Q=60H/△t;
对于线路l的所有途径车站j=1,…,Nl,变量vl,j,s表示车辆离开第j个站点时的电量状态是否为s,当变量为1时,即表示电量状态为s,变量为0时,即表示电量状态不为s,同理,变量vl,1,Q=1表示线路l的车辆离开第1个站点时的电量状态为Q,变量vl,j-1,s表示车辆离开第j-1个站点时的电量状态是否为s;
G为电动公交运输网络,G=(N,L),每条线路l的车辆最低电量保护值为Gl;Nl表示线路l途径的车站总数,n(1,j)表示线路l途径的第j个站点,n(1,j)∈Nl;
Dl,j表示线路l的车辆从其第j个站点至第j+1个站点的行程所消耗的电量,Dl,j-1表示线路l的车辆从其第j-1个站点至第j个站点的行程所消耗的电量;
线路l的车辆在每个站点允许的最大充电次数为Ul,各站点n∈N的充电桩数目为Mn。
在某些实施方式中,所述拉格朗日松弛问题P’(μ,λ)具体为:
s.t.
其中,μl,j和λn,t为步骤S2中所述的拉格朗日乘子。
在某些实施方式中,所述子问题具有|L|个,子问题Pl’是一个包含个0-1变量和条约束的0-1整数规划模型,其中,|S|表示离散集合S中元素的个数,al,n表示线路l车辆在n站点预定的到站时间,dl,n表示线路l车辆在n站点预定的离站时间,第l个子问题Pl’描述为:
s.t.
其中,充电方案由决策变量ul,n tt’或ul,n t’t”表示,当车辆在时间间隔[t,t’]或[t’,t”]时,正于线路l的第n个站点充电时ul,n tt’=1或ul,n t’t”=1。
在某些实施方式中,向所述子问题Pl’内引入新决策变量,表示[t,t+1]时,线路l的车辆正于其第j个站点充电,线路l的车辆未在第j个站点充电时,,其中j=1,…,Nl,t∈{al,j,…,dl,j},将子问题Pl’等价转换紧凑子问题Pl”,紧凑子问题Pl”具有|L|个,第l个紧凑子问题Pl”可以描述为:
s.t.
其中,变量表示线路l的车辆离开第j个站点时的电池电量状态,表示线路l的车辆离开第1个站点时的电池电量状态,表示线路l的车辆离开第j-1个站点时的电池电量状态,为新决策变量,当其为1时,表示[t,t+1]时,线路l的车辆正于其第j-1个站点充电。
在某些实施方式中,所述步骤S4具体为:步骤S41,令线路l的车辆于第j个站点的充电次数Ul,j=0,各站点t间隔内同时充电的车辆数mn,t=0;步骤S42,遍历间隔t∈[al,j+1,dl,j],若且时,令Ul,j+1;遍历间隔t∈[al,j,dl,j],;步骤S43,进而获得,,n∈N,t∈T,则所述拉格朗日松弛问题P’(μ,λ)的目标函数式为
;所述原问题模型P的目标函数式为;步骤S44,分别用μ,λ表示由所述拉格朗日松弛问题P’(μ,λ)中拉格朗日乘子组成的向量,令L(μ,λ)表示所述拉格朗日松弛问题P’(μ,λ)对应的最优解,则对原问题模型P的下界的求解可以转化为求解拉格朗日对偶问题PLD:;步骤S45,利用次梯度算法对所述拉格朗日对偶问题PLD进行求解。
在某些实施方式中,步骤S45具体包括:步骤S01,将所述拉格朗日对偶问题PLD进行初始化,设置当前迭代次数i=0,初始化拉格朗日乘子μ,λ,设置初始下界lb=-∞与上界ub=+∞;步骤S02,更新乘子μ,λ并使用商用求解器求解紧凑子问题Pl”,获得最优解,;步骤S03,利用步骤S41-步骤S44计算Ul,j (i)和mn,t (i)并获得Pl”的目标函数值objLR (i),利用k-greedy算法求解原问题并获得一个可行解对应的目标函数值obj(i),记录当前迭代结果,若松弛约束,即Ul,j (i)≤Ul和mn,t (i)≤Mn均被满足,则检验并更新上下界:若objLR (i)>lb则更新下界lb=objLR (i);若obj(i)<ub则更新上界ub=obj(i);步骤S04,利用次梯度更新拉格朗日乘子,令乘子更新步长为δ=1/2(n+1),采用如下公式更新乘子:μl,j (i+1)=max{μl,j (i)+δ·△μl,j (i),0},λn,t (i+1)=max{λn,t (i)+δ·△λn,t (i),0},其中
步骤S05,判断算法终止条件是否满足,若所有互补松弛条件均被满足时,即μl,j (i+1)(Ul,j (i)-Ul)=0且λn,t (i+1)(mn,t (i)-Mn)=0,算法终止,此时拉格朗日松弛问题对应的最优解即为原问题最优解,上述上标(i)和(i+1)分别表示第i次迭代和第i+1次迭代;当lb=ub时,算法终止,此时拉格朗日松弛问题对应的最优解即为原问题最优解;当迭代次数i超过最大迭代次数Imax时,算法终止,输出上下界对应的目标函数和最优解;若上述三个条件均未满足,循环步骤S01-步骤S05,迭代次数依次+1,直至上述三个条件中的一个满足,得出原问题最优解,即获得各电动公交车辆日间运行充电方案,完成对电动公交日间运行充电的优化。
本发明的范围,并不限于上述技术特征的特定组合而成的技术方案,同时也应涵盖由上述技术特征或其等同特征进行任意组合而形成的其它技术方案。例如上述特征与本申请中公开的(但不限于)具有类似功能的技术特征进行互相替换而形成的技术方案等。
由于上述技术方案运用,本发明与现有技术相比具有下列优点:本发明提供一种基于拉格朗日松弛算法的电动公交日间运行充电优化方法,在保证车辆预定时刻表的前提下,以最小化***总充电时间为目标,基于离散电池电量状态获得全线路网络级充电方案整数优化原问题模型,然后转换为拉格朗日松弛问题,接着分解为多个子问题,最后转化为多个紧凑子问题,通过求解多个紧凑子问题的最优解,获得还原后的原问题模型P的最优解,精准描述大规模公交网络的充电问题,一步步减少模型中的决策变量和约束条件数量,缩减问题规模,提升求解效率,进行快速高效的求解。
附图说明
附图1是本发明基于拉格朗日松弛算法的电动公交日间运行充电优化策略的流程示意图;
附图2是本发明k-greedy算法流程图。
具体实施方式
如附图1所示的一种基于拉格朗日松弛算法的电动公交日间运行充电优化方法,包括以下具体步骤:
步骤S1:获取电动公交数据和包含线路数量的线路运行数据,电动公交数据包括每辆公交车在各站点的最大充电次数、车辆最低电量保护值,线路运行数据包括途径站点、各站点的到发时间表、各站点的充电桩数目,构建基于离散电池电量状态的全线路网络级充电方案整数优化模型,获得包含耦合约束的原问题模型P,原问题模型P具体为:
s.t.
其中,s.t.表示约束条件为,
和为耦合约束,G为电动公交运输网络,G=(N,L),N和L分别表示公交站点位置和公交线路数量,l表示第l条公交线路,l∈L,Nl表示线路l途径的车站总数,n(1,j)表示该线路途径的第j个站点,n(1,j)∈Nl,al,j表示线路l的车辆在第j个站点预定的到站时间,dl,j表示线路l的车辆在第j个站点预定的离站时间,Dl,j表示线路l的车辆从其第j个站点至第j+1个站点的行程所消耗的电量,Dl,j-1表示线路l的车辆从其第j-1个站点至第j个站点的行程所消耗的电量,j∈{1,…,Nl-1},将规划时段T离散为以△tmin为单位的间隔,t,t’,t”均为规划时段T内的某一时间点,电动公交电量状态SOC表示为离散集合S={0,1,2,…,Q},电池从SOC=0的状态充至满电所需的时间为H小时,Q=60H/△t,集合中的每个值s∈S代表车辆从SOC=0状态开始充电s个时间间隔△t后所达到的电量状态,集合中的每个值s∈S代表车辆从SOC=0状态开始充电s个时间间隔△t后所达到的电量状态,即s=0表示车辆充电△t时长后电池达到的电量,s=2表示车辆充电2△t时长后电池达到的电量;由于电池的充电过程是非线性的,因此每个状态s∈S之间电池电量的增量是不一致的,即s和s+1之间的电量状态差值不一定等于s-1和s间的差值;为延长电池使用寿命,每条线路l的车辆最低电量保护值为Gl,即车辆在每次行程结束时SOC不能够低于Gl,线路l的车辆在每个站点允许的最大充电次数为Ul,考虑到电池电量在行程中的消耗,车辆需要在经过站点时充电以满足后续行程的电量需求,但由于车辆到离车站时间不同且同一车站的充电桩数目有限,先到的车辆可以中断充电并让后到且时刻表紧急的车辆充电,待后者离开后再继续充电,但每插拔一次均算入车辆充电次数,各站点n∈N的充电桩数目为Mn,充电方案由决策变量ul,j tt’或ul,j t’t”表示,当车辆在时间间隔[t,t’]或[t’,t”]时正于线路l的第j个站点充电时ul,j tt’=1或ul,j t’t”=1,对于线路l的所有途径车站j=1,…,Nl,变量vl,j,s表示车辆离开第j个站点时的电量状态是否为s,当变量为1时,即表示电量状态为s,变量为0时,即表示电量状态不为s,变量vl,1,Q=1表示线路l的车辆离开第1个站点时的电量状态为Q,变量vl,j-1,s表示车辆离开第j-1个站点时的电量状态是否为s。
原问题模型P中,具有两种类型的决策变量,分别为电动公交的充电时间表和电动公交离开车站时的最佳电池状态,原问题模型P中,具有三种类型的约束,分别为与线路时刻表相关的约束、电池电量状态约束、耦合约束。
步骤S2,通过引入拉格朗日乘子μl,j和λn,t且μl,j≥0,λn,t≥0,将所述耦合约束条件转换到原问题P的目标函数中,从而获得具有新的约束条件的拉格朗日松弛问题问题P’(μ,λ),所述拉格朗日松弛问题问题P’(μ,λ)具体为:
s.t.
步骤S03,将上述拉格朗日松弛问题P’(μ,λ)目标函数中的常数项,即移除,然后根据线路数量,将拉格朗日松弛问题P’(μ,λ)分解
为|L|个子问题,每个子问题分别是对各条线路充电方案的决策优化,其中第l个子问题Pl’
可以描述为:第l个子问题Pl’描述为:
s.t.
充电方案由决策变量ul,n tt’或ul,n t’t”表示,当车辆在时间间隔[t,t’]或[t’,t”]
时,正于线路l的第n个站点充电时ul,n tt’=1或ul,n t’t”=1,子问题Pl’是一个包含个0-1变量和条约束的0-1整数规划
模型,在大规模问题中,当每条线路的运行时刻表延续时间较长时,求解上述模型将变得极
其困难,所以继续引入新的决策变量,表示[t,t+1]时,线路l的车辆正于其第j
个站点充电,否则,其中j=1,…,Nl,t∈{al,j,…,dl,j};变量表示线路l的车辆
离开第j个站点时的电池电量状态,表示线路l的车辆离开第1个站点时的电池电量状
态,表示线路l的车辆离开第j-1个站点时的电池电量状态,为新决策变量,当其
为1时,表示[t,t+1]时,线路l的车辆正于其第j-1个站点充电;从而将多个所述子问题Pl’
等价转换紧凑子问题Pl”,紧凑子问题Pl”描述为:
s.t.
步骤S4,利用次梯度下降法求解多个所述紧凑子问题Pl”,每次迭代后根据当前解更新拉格朗日乘子,直至算法达到终止条件,获得原问题模型P的最优解,具体步骤为:
步骤S41,令线路l的车辆于第j个站点的充电次数Ul,j=0,各站点t间隔内同时充电的车辆数mn,t=0;
步骤S43,进而获得
则所述拉格朗日松弛问题P’(μ,λ)的目标函数式
等于
步骤S44,分别用μ,λ表示由所述拉格朗日松弛问题P’(μ,λ)中拉格朗日乘子组成的向量,令L(μ,λ)表示所述拉格朗日松弛问题P’(μ,λ)对应的最优解,则对原问题模型P的下界的求解可以转化为求解拉格朗日对偶问题:
步骤S45,将所述拉格朗日对偶问题PLD进行初始化,设置当前迭代次数i=0,初始化拉格朗日乘子μ,λ,设置初始下界lb=-∞与上界ub=+∞;更新乘子μ,λ并使用商用求解器求解紧凑子问题Pl”,获得最优解,;利用步骤S41-步骤S44计算Ul,j (i)和mn,t (i)并获得Pl”的目标函数值objLR (i),利用k-greedy算法求解原问题并获得一个可行解对应的目标函数值obj(i),k-greedy算法求解如图2所示,记录当前迭代结果,若松弛约束,即Ul,j (i)≤Ul和mn,t (i)≤Mn均被满足,则检验并更新上下界:若objLR (i)>lb则更新下界lb=objLR (i);若objLR (i)>lb则更新下界lb=objLR (i);利用次梯度更新拉格朗日乘子,令乘子更新步长为δ=1/2(n+1),采用如下公式更新乘子:
其中
判断算法终止条件是否满足,若所有互补松弛条件均被满足时,即μl,j (i+1)(Ul,j (i)-Ul)=0且λn,t (i+1)(mn,t (i)-Mn)=0,算法终止,此时拉格朗日松弛问题对应的最优解即为原问题最优解;当lb=ub时,算法终止,此时拉格朗日松弛问题对应的最优解即为原问题最优解;当迭代次数i超过最大迭代次数Imax时,算法终止,输出上下界对应的目标函数和最优解;若上述三个条件均未满足,循环步骤S01-步骤S05,迭代次数依次+1,直至上述三个条件中的一个满足,得出原问题最优解,即获得各电动公交车辆日间运行充电方案,完成对电动公交日间运行充电的优化。
实施例:本实施例使用了网站公开的公交网络和日间运营数据(包括车站的位置和车辆行程安排等),规划时间范围T划分为600个时间间隔△t,每个时间间隔表示1min,车辆电池容量Cap=100,电池最低电量保护值Gl=5。因为原问题P与转化后的子问题Pl’和紧凑子问题Pl”都描述为整数规划模型,均可以采用商用求解器(如Cplex、Gurobi等)直接求解。
为突显本申请方案的有效性,实施例比较了:1)使用Gurobi直接求解原问题P,和2)使用步骤S41-S45求解其中紧凑子问题Pl”使用Gurobi求解两种方法的最优解和求解效率,结果如表1所示:
表1中,算例取值N、L、Mn、Ul分别表示:网络中的站点总数、线路总数、各站点n∈N的充电桩数目Mn,各线路l的车辆在每个站点允许的最大充电次数为Ul;目标函数列表示最优解对应的目标函数值,求解耗时是以秒为单位的CPU计算时间。结果表明,本申请提出的方案能够有效解决电动公交日间运行充电优化问题,求解效率明显优于直接使用Gurobi解算器的解决方案;且当站点数和线路分别数达到300/150时,由于问题规模较大Gurobi已无法求解,而本申请方案仍表现出优良的求解能力,本申请方案比常规方法更稳定、高效。
上述实施例只为说明本发明的技术构思及特点,其目的在于让熟悉此项技术的人士能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡根据本发明精神实质所作的等效变化或修饰,都应涵盖在本发明的保护范围之内。
Claims (8)
1.一种基于拉格朗日松弛算法的电动公交日间运行充电优化方法,其特征在于:包括以下具体步骤:
S1、获取电动公交数据和包含线路数量的线路运行数据,构建基于离散电池电量状态的全线路网络级充电方案整数优化模型,获得包含耦合约束的原问题模型P;
S2、通过引入拉格朗日乘子将所述耦合约束转换到原问题模型P的目标函数中,从而获得具有新的约束条件的拉格朗日松弛问题P’(μ,λ);
S3、根据所述线路数量,将所述拉格朗日松弛问题P’(μ,λ)分解为多个子问题Pl’,引入新的决策变量,即相应线路的车辆是否在相应站点充电的决策变量,从而将多个所述子问题Pl’转化为多个紧凑子问题Pl”;
S4、通过求解多个紧凑子问题Pl”,求解出P’(μ,λ)的目标函数,即可求解出原问题模型P的目标函数,
将求解原问题模型P的下界的最优解转化为求解拉格朗日对偶问题PLD,利用次梯度下降法求解拉格朗日对偶问题PLD,并在每次迭代求解后根据当前解更新拉格朗日乘子,
采用k-greedy算法计算原问题模型P的上界的最优解,
直至原问题模型P的上界的最优解等于下界的最优解,算法达到终止条件,获得原问题模型P的最优解,即获得各电动公交车辆日间运行充电方案,完成对电动公交日间运行充电的优化。
2.根据权利要求1所述的基于拉格朗日松弛算法的电动公交日间运行充电优化方法,其特征在于:所述电动公交数据包括每辆公交车在各站点的最大充电次数、车辆最低电量保护值、电动公交电量状态,所述线路运行数据包括途径站点、各站点的到发时间表、各站点的充电桩数目。
s.t.
其中,
s.t.表示约束条件为;
N和L分别表示公交站点位置和公交线路数量,l表示第l条公交线路,l∈L,n表示第n个公交站点,n∈N,Nl表示线路l途径的车站总数;
al,j表示线路l的车辆在第j个站点预定的到站时间,dl,j表示线路l的车辆在第j个站点预定的离站时间,j∈{1,…,Nl-1},将规划时段T离散为以△tmin为单位的间隔,t,t’,t”均为规划时段T内的某一时间点,充电方案由决策变量ul,j tt’或ul,j t’t”表示,当车辆在时间间隔[t,t’]或[t’,t”]时正于线路l的第j个站点充电时ul,j tt’=1或ul,j t’t”=1;
S为离散集合,S表示电动公交电量状态SOC,S={0,1,2,…,Q},集合S中的每个值s∈S代表车辆从SOC=0状态开始充电s个时间间隔△t后所达到的电量状态,即s=0表示车辆充电△t时长后电池达到的电量,s=2表示车辆充电2△t时长后电池达到的电量,令电池从SOC=0的状态充至满电所需的时间为H小时,Q表示电量状态Q,即满电状态,则Q=60H/△t;
对于线路l的所有途径车站j=1,…,Nl,变量vl,j,s表示车辆离开第j个站点时的电量状态是否为s,当变量为1时,即表示电量状态为s,变量为0时,即表示电量状态不为s,同理,变量vl,1,Q=1表示线路l的车辆离开第1个站点时的电量状态为Q,变量vl,j-1,s表示车辆离开第j-1个站点时的电量状态是否为s;
G为电动公交运输网络,G=(N,L),每条线路l的车辆最低电量保护值为Gl;Nl表示线路l途径的车站总数,n(1,j)表示线路l途径的第j个站点,n(1,j)∈Nl;
Dl,j表示线路l的车辆从其第j个站点至第j+1个站点的行程所消耗的电量,Dl,j-1表示线路l的车辆从其第j-1个站点至第j个站点的行程所消耗的电量;
线路l的车辆在每个站点允许的最大充电次数为Ul,各站点n∈N的充电桩数目为Mn。
6.根据权利要求5所述的基于拉格朗日松弛算法的电动公交日间运行充电优化方法,其特征在于:向所述子问题Pl’内引入新决策变量,表示[t,t+1]时,线路l的车辆正于其第j个站点充电,线路l的车辆未在第j个站点充电时,,其中j=1,…,Nl,t∈{al,j,…,dl,j},将子问题Pl’等价转换紧凑子问题Pl”,紧凑子问题Pl”具有|L|个,第l个紧凑子问题Pl”可以描述为:
s.t.
7.根据权利要求6所述的基于拉格朗日松弛算法的电动公交日间运行充电优化方法,其特征在于:所述步骤S4具体为:步骤S41,令线路l的车辆于第j个站点的充电次数Ul,j=0,各站点t间隔内同时充电的车辆数mn,t=0;步骤S42,遍历间隔t∈[al,j+1,dl,j],若且时,令Ul,j+1;遍历间隔t∈[al,j,dl,j],;步骤S43,进而获得,,n∈N,t∈T,则所述拉格朗日松弛问题P’(μ,λ)的目标函数式为
8.根据权利要求7所述的基于拉格朗日松弛算法的电动公交日间运行充电优化方法,其特征在于:步骤S45具体包括:步骤S01,将所述拉格朗日对偶问题PLD进行初始化,设置当前迭代次数i=0,初始化拉格朗日乘子μ,λ,设置初始下界lb=-∞与上界ub=+∞;步骤S02,更新乘子μ,λ并使用商用求解器求解紧凑子问题Pl”,获得最优解,;步骤S03,利用步骤S41-步骤S44计算Ul,j (i)和mn,t (i)并获得Pl”的目标函数值objLR (i),利用k-greedy算法求解原问题并获得一个可行解对应的目标函数值obj(i),记录当前迭代结果,若松弛约束,即Ul,j (i)≤Ul和mn,t (i)≤Mn均被满足,则检验并更新上下界:若objLR (i)>lb则更新下界lb=objLR (i);若obj(i)<ub则更新上界ub=obj(i);步骤S04,利用次梯度更新拉格朗日乘子,令乘子更新步长为δ=1/2(n+1),采用如下公式更新乘子:μl,j (i+1)=max{μl,j (i)+δ·△μl,j (i),0},λn,t (i+1)=max{λn,t (i)+δ·△λn,t (i),0},其中
步骤S05,判断算法终止条件是否满足,若所有互补松弛条件均被满足时,即μl,j (i+1)(Ul,j (i)-Ul)=0且λn,t (i+1)(mn,t (i)-Mn)=0,算法终止,此时拉格朗日松弛问题对应的最优解即为原问题最优解,上述上标(i)和(i+1)分别表示第i次迭代和第i+1次迭代;当lb=ub时,算法终止,此时拉格朗日松弛问题对应的最优解即为原问题最优解;当迭代次数i超过最大迭代次数Imax时,算法终止,输出上下界对应的目标函数和最优解;若上述三个条件均未满足,循环步骤S01-步骤S05,迭代次数依次+1,直至上述三个条件中的一个满足,得出原问题最优解,即获得各电动公交车辆日间运行充电方案,完成对电动公交日间运行充电的优化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210874670.0A CN115099702B (zh) | 2022-07-25 | 2022-07-25 | 基于拉格朗日松弛算法的电动公交日间运行充电优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210874670.0A CN115099702B (zh) | 2022-07-25 | 2022-07-25 | 基于拉格朗日松弛算法的电动公交日间运行充电优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115099702A true CN115099702A (zh) | 2022-09-23 |
CN115099702B CN115099702B (zh) | 2022-11-15 |
Family
ID=83298144
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210874670.0A Active CN115099702B (zh) | 2022-07-25 | 2022-07-25 | 基于拉格朗日松弛算法的电动公交日间运行充电优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115099702B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115729193A (zh) * | 2022-11-24 | 2023-03-03 | 沈阳建筑大学 | 一种炼钢-连铸中间包批量计划优化方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113361073A (zh) * | 2021-05-10 | 2021-09-07 | 上海工程技术大学 | 一种基于改进拉格朗日松弛算法的库存路径联合优化方法 |
CN113887032A (zh) * | 2021-09-26 | 2022-01-04 | 云南电网有限责任公司 | 基于拉格朗日分布式算法的电动汽车有序充放电控制方法 |
-
2022
- 2022-07-25 CN CN202210874670.0A patent/CN115099702B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113361073A (zh) * | 2021-05-10 | 2021-09-07 | 上海工程技术大学 | 一种基于改进拉格朗日松弛算法的库存路径联合优化方法 |
CN113887032A (zh) * | 2021-09-26 | 2022-01-04 | 云南电网有限责任公司 | 基于拉格朗日分布式算法的电动汽车有序充放电控制方法 |
Non-Patent Citations (1)
Title |
---|
田小鹏 等: "基于拉格朗日分解的电动汽车充电路径优化", 《兰州交通大学学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115729193A (zh) * | 2022-11-24 | 2023-03-03 | 沈阳建筑大学 | 一种炼钢-连铸中间包批量计划优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115099702B (zh) | 2022-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wu et al. | A novel fast-charging stations locational planning model for electric bus transit system | |
CN108955711B (zh) | 一种应用于电动汽车智能充放电的导航方法 | |
US11087422B2 (en) | Method and device for determining a configuration for deployment of a public transportation system | |
CN103915869B (zh) | 一种基于移动设备的电动汽车智能充电***及方法 | |
CN111311122B (zh) | 面向智慧城市的充电站建设方法 | |
CN109489676B (zh) | 一种计及电网信息与充电站信息的电动汽车充电导航方法 | |
El-Taweel et al. | Integrated utility-transit model for optimal configuration of battery electric bus systems | |
Cheng et al. | Electric bus fast charging station resource planning considering load aggregation and renewable integration | |
CN107274665A (zh) | 巴士运力资源规划方法及*** | |
CN111047120B (zh) | 一种路-电耦合网络下电动汽车充电负荷预测方法 | |
Yuan et al. | Optimal scheduling for micro‐grid considering EV charging–swapping–storage integrated station | |
CN115115268B (zh) | 基于路电耦合与低碳约束的电动汽车充电桩容量规划方法 | |
CN115099702B (zh) | 基于拉格朗日松弛算法的电动公交日间运行充电优化方法 | |
CN112097783B (zh) | 基于深度强化学习的电动出租车充电导航路径规划方法 | |
Sivagnanam et al. | Minimizing energy use of mixed-fleet public transit for fixed-route service | |
CN114462864A (zh) | 一种充电设施共享策略影响下的电动公交线路车辆调度方法 | |
CN112116125A (zh) | 一种基于深度强化学习的电动汽车充电导航方法 | |
He et al. | Expansion planning of electric vehicle charging stations considering the benefits of peak‐regulation frequency modulation | |
Ruiz et al. | An optimal battery charging and schedule control strategy for electric bus rapid transit | |
Dalla Chiara et al. | Electrifying Buses for Public Transport: Boundaries with a Performance Analysis Based on Method and Experience | |
CN113222248B (zh) | 一种自动驾驶出租车充电桩选择方法 | |
Su et al. | Cost minimization control for electric vehicle car parks with vehicle to grid technology | |
CN105471043A (zh) | 充电设备的配置方法和装置 | |
CN115983568A (zh) | 一种考虑电池健康状态差异的电动公交线路车辆调度方法 | |
Liu et al. | An MHO approach for electric bus charging scheme optimization based on energy consumption estimation |
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 | ||
CP03 | Change of name, title or address |
Address after: 35 / F, block a, Suzhou City Life Plaza, 251 pinglong Road, Gusu District, Suzhou City, Jiangsu Province 215000 Patentee after: Jiangsu Boyuxin Information Technology Co.,Ltd. Country or region after: China Address before: 35 / F, block a, Suzhou City Life Plaza, 251 pinglong Road, Gusu District, Suzhou City, Jiangsu Province 215000 Patentee before: Suzhou BOYUXIN Transportation Technology Co.,Ltd. Country or region before: China |
|
CP03 | Change of name, title or address |