CN111582581B - 考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法 - Google Patents
考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法 Download PDFInfo
- Publication number
- CN111582581B CN111582581B CN202010383239.7A CN202010383239A CN111582581B CN 111582581 B CN111582581 B CN 111582581B CN 202010383239 A CN202010383239 A CN 202010383239A CN 111582581 B CN111582581 B CN 111582581B
- Authority
- CN
- China
- Prior art keywords
- path
- decision
- model
- dfrl
- robust
- 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.)
- Active
Links
- 238000010187 selection method Methods 0.000 title claims abstract description 12
- 238000000034 method Methods 0.000 claims abstract description 19
- 238000005315 distribution function Methods 0.000 claims abstract description 11
- 230000008569 process Effects 0.000 claims abstract description 6
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000009977 dual effect Effects 0.000 claims description 4
- 239000013598 vector Substances 0.000 claims description 4
- 230000008859 change Effects 0.000 description 8
- 238000005457 optimization Methods 0.000 description 7
- 230000006399 behavior Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 101150096839 Fcmr gene Proteins 0.000 description 1
- 241000764238 Isis Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000010248 power generation 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
- G06Q10/043—Optimisation of two dimensional placement, e.g. cutting of clothes or wood
-
- 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
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0201—Market modelling; Market analysis; Collecting market data
-
- 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
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0201—Market modelling; Market analysis; Collecting market data
- G06Q30/0202—Market predictions or forecasting for commercial activities
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S50/00—Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
- Y04S50/14—Marketing, i.e. market research and analysis, surveying, promotions, advertising, buyer profiling, customer management or rewards
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Economics (AREA)
- Accounting & Taxation (AREA)
- Finance (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Game Theory and Decision Science (AREA)
- Tourism & Hospitality (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法,包括:获取运输网络所对应的OD对集合和选址点集合;根据OD对集合生成路径,并根据线性衰减函数表示路径需求截获比例;利用不确定分布函数集合建立DR‑DFRL模型;通过增加虚拟节点对运输网络进行扩展;基于扩展网络,以选址决策、路径选择决策和服务单个用户收益决策为决策变量,重构DR‑DFRL模型;利用外部估计算法求解DR‑DFRL模型,求解过程中迭代加入约束条件,得到符合约束条件的决策变量最优解。通过本发明的技术方案,同时考虑用户路径选择行为和需求不确定性,使决策更加符合实际运营情况,实用性更强,降低了决策风险,且可实现高效求解。
Description
技术领域
本发明涉及运营管理及运筹学技术领域,尤其涉及一种考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法。
背景技术
随着全球能源危机与气候变化的加剧,世界各国都在寻求解决方案。在运输领域,电动汽车的推广受到了大力支持。当前,电动汽车能源补充有两种模式:充电模式和换电模式。但是充电模式有以下不足之处:充电时间长、投资成本高和车辆行驶里程有限。为了解决电动汽车充电模式的不足,电动汽车行业开始探索换电模式。虽然换电模式在2013年由于Better Place公司的破产而宣布失败。但是近年来,多家企业开始重新探索换电模式。例如,北汽新能源推出“擎天柱计划”,旨在将新能源汽车、动力电池、换电站、光伏发电进行深度融合,计划到2022年在全国超过100个城市,建成3000座光储换电站,投放换电车辆50万台。
电动汽车换电站的选址对于换电服务运营至关重要。决策者需要在短期投资成本和长期收益之间权衡,这主要取决于电动汽车换电需求量。关于选址问题,最早出现的是确定性的选址模型。确定性选址模型可以分为三类:基于节点的选址模型、基于弧的选址模型和基于路径的选址模型。目前,大多数选址模型都是基于路径的模型。Hodgson最早开始研究基于路径的选址模型,他提出了截流选址(Flow Capturing Location,FCL)模型,在给定建站个数的约束下,该模型的目标是最大化给定起始终点(Origin-Destination,OD)对间服务的总需求。该模型假设一条路径上只要有一个站,这条路径上的需求即可被服务,但是该模型没有考虑到车辆的有限续航里程。在Hodgson提出FCL模型之后,很多学者提出了改进的模型:考虑车辆续航里程的流量续航选址(Flow Refueling Location,FRL)模型,考虑车站容量约束的有约束的流量续航选址(Capacitated Flow Capturing Location,CFRL)模型等。但上述模型都需要预生成大量的可行选址点组合,因此有学者对其进行改进提出新的模型:新的基于节点和基于弧的FRL模型,扩展网络模型。但这些模型都没有考虑到需求的不确定性和用户路径选择行为。
从宏观层面上讲,给定的OD对之间的换电需求取决于人口分布、电动汽车普及率以及其他的经济技术因素,需求通常都是动态变化的,因此在基础设施规划阶段需求是不确定的。确定性的流量续航选址模型认为每一个OD对间的需求都是确定的,而不确定性的存在使得确定性问题的最优解并非最优甚至不可行。故而出现了针对换电站选址问题中不确定性的研究,有随机选址模型和鲁棒选址模型。
从微观层面上讲,换电站服务的换电需求还取决于用户的路径选择。由于有限的车辆行驶里程和换电站数量,最短路径并不总是可行的。用户可能必须要走非最短路才可行。目前的路径偏移模型认为旅行需求随偏移距离的增加而减少,但是大都认为需求是确定的。
发明内容
针对上述问题,本发明提供了一种考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法,同时考虑用户路径选择行为和需求不确定性,采用分布式鲁棒优化方法,提出了分布鲁棒-偏移流量续航选址(Distributionally Robust-Deviation Flow RefuelingLocation,DR-DFRL)模型。该模型为实际的换电站选址提供决策,更加符合实际运营情况,实用性更强。该模型利用了不确定参数的均值和协方差信息,相比于传统的鲁棒优化模型,保守性更低。此外,采用线性衰减函数表示用户路径选择行为需求量,使得DR-DFRL模型能够构建为混合整数非线性规划问题。该非线性规划模型可采用外部估计(OuterApproximation,OA)算法高效求解。相比已有的选址建模方法得到更加符合实际情况的选址结果,从而降低了决策风险。
为实现上述目的,本发明提供了一种考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法,包括:获取运输网络所对应的OD对集合和选址点集合;根据所述OD对集合生成路径,并根据线性衰减函数表示路径需求截获比例;利用不确定分布函数集合建立DR-DFRL模型;通过增加虚拟节点对所述运输网络进行扩展,得到扩展网络;基于所述扩展网络,以选址决策、路径选择决策和服务单个用户收益决策为决策变量,重构所述DR-DFRL模型;利用OA算法求解所述DR-DFRL模型,求解过程中迭代加入约束条件,得到符合约束条件的决策变量最优解。
在上述技术方案中,优选地,利用线性衰减函数表达所述OD对间截获的路径需求截获比例,根据所述路径需求截获比例与所述OD对的随机需求得到所述OD对间截获的换电需求量,并根据换电服务费用、车辆行驶里程和路径选择决策变量,得到所述运输网络上的服务单个用户收益决策变量。
在上述技术方案中,优选地,所述约束条件包括流量平衡约束、节点允许流量通过约束、路径选择数量约束、服务单个用户收益约束和自变量定义域约束。
在上述技术方案中,优选地,采用K-最短路算法获得符合所述用户绕行偏移距离的路径,并将不符合所述用户绕行偏移距离的路径及相应的节点删除以获得简化后的所述运输网络。
在上述技术方案中,优选地,通过在所述运输网络中增加虚拟起点和虚拟终点,基于换电站建站规则进行路径扩展,构成所述扩展网络。
在上述技术方案中,优选地,利用OA算法对所述DR-DFRL模型进行求解过程中,通过求解原始非线性规划问题的线性近似主问题得到原始问题的上界可行解,将所述上界可行解代入所述DR-DFRL模型的目标函数,得到原始问题的最优值的一个下界值,逐一迭代加入所述约束条件,得到符合所有约束条件的最优解。
在上述技术方案中,优选地,在所述OA算法中,在迭代次数达到预设的最大迭代次数或相对间隙小于等于预设的相对间隙允许值时,停止迭代。
在上述技术方案中,优选地,所述OD对上的换电需求构成向量,根据所述换电需求的均值和协方差矩阵得到所述换电需求可能的概率密度函数所属的分布函数集合;利用所述分布函数集合构建出目标函数,以寻求最坏情况下平均利润最大化的最优选址集合;根据拉格朗日对偶理论,将所述目标函数转换为所述DR-DFRL模型的目标函数。
与现有技术相比,本发明的有益效果为:通过同时考虑用户路径选择行为和需求不确定性,采用分布式鲁棒优化方法,提出了DR-DFRL模型,为实际的换电站选址提供决策,更加符合实际运营情况,实用性更强,利用了不确定参数的均值和协方差信息,相比于传统的鲁棒优化模型,保守性更低。此外,采用线性衰减函数表示用户路径选择行为需求量,使得DR-DFRL模型能够构建为混合整数非线性规划问题,采用OA算法可高效求解,相比已有的选址建模方法得到更加符合实际情况的选址结果,从而降低了决策风险。
附图说明
图1为本发明一种实施例公开的考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法的流程示意图;
图2为本发明一种实施例公开的运输网络的路径构成示意图;
图3为本发明一种实施例公开的路径的扩展网络示意图;
图4为本发明一种实施例公开的OA算法的流程示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
下面结合附图对本发明做进一步的详细描述:
如图1示,根据本发明提供的一种考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法,包括:获取运输网络所对应的OD对集合和选址点集合;根据OD对集合生成路径,并根据线性衰减函数表示路径需求截获比例;利用不确定分布函数集合建立DR-DFRL模型;通过增加虚拟节点对运输网络进行扩展,得到扩展网络;基于扩展网络,以选址决策、路径选择决策和服务单个用户收益决策为决策变量,重构DR-DFRL模型;利用OA算法求解DR-DFRL模型,求解过程中迭代加入约束条件,得到符合约束条件的决策变量最优解。
但由于车辆行驶里程的限制以及有限的换电站数量,用户选择OD对间的最短路并不一定能够完成往返里程。因此,电动汽车用户会根据现有的换电站建站信息选择成本最低的偏移路径。但由于偏移距离可能为用户带来不便,并增加旅行成本,因此,并不是所有的用户都会选择偏移路径,而是只有一部分用户选择偏移路径。
在该实施例中,同时考虑用户路径选择行为和需求不确定性,采用较新的分布式鲁棒优化方法,提出了DR-DFRL模型,为实际的换电站选址提供决策。该DR-DFRL模型是两阶段优化模型。在第一阶段,决策者从一系列的候选节点中选址建站,此时,我们只知道换电需求的均值-协方差信息。在第二阶段,换电站建设完成,则用户换电需求量确定,用户可选择的最短出行路径确定。
在该实施例中,运输网络由G=(N,A)表示,其中N是网络上的节点集合,A是网络上有向弧的集合。N中所有的节点都可以当做起点、终点和候选建站节点,因此,给定N可以得到个起始-终点(OD)对,从所有的OD对集合中选择有往返需求的部分组成集合Q,fq表示OD对q(q∈Q)间的随机需求。假设均值为E(f)=μ,协方差矩阵为Cov(f)=Σ,令为换电站候选点集合。
假设用户的出行成本只与距离有关,对用户而言,成本越高,用户选择该方式出行的意愿就越低。对任意OD对q,用表示OD对q间最短路径长度。当前换电站建设处于起步阶段,由于地形条件以及投资成本的限制,换电站建设数量较少,用户选择最短路径不一定能够完成行驶里程。由于顾客愿意行驶的绕远距离有限,对任意OD对q我们考虑在长度不超过(1+β)Distq的路径中选择一条服务该OD对间的需求。
优选地,采用K-最短路算法获得符合用户绕行偏移距离的路径,并将不符合用户绕行偏移距离的路径及相应的节点删除以获得简化后的运输网络。具体地,在本发明中,采用K-最短路算法获得长度不大于(1+β)Distq的路径的集合Pq。因此,能够服务Q中所有OD对的路径集合为P=∪p∈QPq,P中所有路径上所有节点构成选址点集合基于P与F,原运输网络上不在集合P中路径上的节点与弧均可以删掉,得到简化网络。不失一般性地,以下仍然用G=(N,A)表示简化后的网络。
在上述实施例中,优选地,利用线性衰减函数表达OD对间截获的路径需求截获比例,根据路径需求截获比例与OD对的随机需求得到OD对间截获的换电需求量,并根据换电服务费用、车辆行驶里程和路径选择决策变量,得到运输网络上的服务单个用户收益决策变量。
具体地,当选址组合给定时,可以得到能够服务某个OD对q间需求的所有可行路径集合,假设是第一阶段中决策出的选址点集合,令P(S,q)为可服务OD对q间需求的路径集合,其中从而可得到整个运输网络上选址给定时可选路径集合P=∪q∈QP(S,q)。由于用户出行考虑的因素是成本,因此每个OD对间的用户都会在P(S,q)中选择成本最低的路p(S,q),即
p(S,q)∈argm inp∈P(S,q)Dist(p). (1)其中,Dist(p)是路径p的路径长度。
令Distq为OD对q间最短路径的距离,那么对于OD对q间选定的路径p,用户选择的路径p偏移距离为Δ(S,q)=Dist(p(S,q))-Distq,OD对q(q∈Q)间截获的路径需求截获比例可用如下线性衰减函数表示:
其中,β>0是给定的参数。线性衰减函数说明OD对q(q∈Q)间没有用户会选择超过(1+β)Distq长度的路径。当路径长度不超过(1+β)Distq时,截获的需求随偏移距离的增加而减少。因此,对于给定的需求fq,OD对q间截获的需求为fq·gq(S)。
令fq表示OD对q(q∈Q)间的随机需求,Q中所有OD对上的需求构成向量f。假设均值为E(f)=μ,协方差矩阵为Cov(f)=Σ。
令s为单个换电站为单个用户换电一次的服务费用,R为单块电池充满电的车辆行驶里程。给定选址集合S,令bq(S)表示路径p(S,q)上的单个用户为完成旅程所需支付的换电费用,OD对q间单个用户选择路径p(S,q)完成旅程所需换电的次数近似为所以路径p(S,q)上的单个用户为完成旅程所需支付的换电费用为给定OD对q间的换电需求fq,运营商服务OD对q间需求获得的换电服务收益为hq(fq,S)=fqgq(S)bq(S),其中αS,q=bq(S)gq(S),则服务运输网络上的换电服务收益为
由于bq(S),gq(S)和αq(S)只与选择的路径有关,因此,可将其用bq,gq和αq表示,并且αp=bpgp。
因此,当给定选址集合S,在需求f的概率密度函数给定的情况下,运输网络上服务需求的平均利润为
其中,ci是换电站i(i∈F)的建站成本。
在上述实施例中,优选地,OD对上的换电需求构成向量,根据换电需求的均值和协方差矩阵得到换电需求可能的概率密度函数所属的分布函数集合;利用分布函数集合构建出目标函数,以寻求最坏情况下平均利润最大化的最优选址集合;根据拉格朗日对偶理论,将目标函数转换为DR-DFRL模型的目标函数。
具体地,实际生活中,在选址阶段,需求f的真实概率密度函数通常是未知的,但已知其均值与协方差信息,假设f的概率密度函数属于以下分布函数集合:
D={f:(E[f]-μ)TΣ-1(E[f]-μ)≤?2} (4)
其中,ò>0是一个正数,决定估计精度。
由于DR-DFRL模型的目标是找到一个最优选址集合S*使得最坏情况下的平均利润最大化。即,目标函数具有以下形式:
令t=E[f],λ∈R+为(6)式中约束的拉格朗日乘子,根据朗格朗日对偶理论,DR-DFRL模型可等价于
由于选址集合S的组合较多,从可能的S集合中计算求解得到最优的S,会使得计算复杂度增加,故而在此采用MirHassani提出的基于扩展网络的方法得到有效的选址组合,在运输网络中增加虚拟起点和虚拟终点,基于换电站建站规则进行路径扩展,构成扩展网络,从而降低计算复杂度。
对每一条路径p∈P,生成与该路径对应的扩展网络Gp=(Np,Ap)。Np是路径p扩展后的节点集合,Ap是路径p扩展后的弧集合。基于2个换电站建站规则(两个换电站之间距离不能超过R;起点处没有换电站时,用户电池只有50%的电量,需能够到达第一个换电站),在此进行路径扩展。以图2为例,原始路径p为A→B→C→D,原始节点集合为Np={A,B,C,D},原始路径上弧的集合Ap={(A,B),(B,C),(C,D)},假设R=100km,每条弧上面的数字表示弧的长度。进行网络扩展的具体步骤如下:
第一步:增加两个虚拟节点到路径p上,即虚拟起点sp和虚拟终点tp,添加(sp,A)、(D,tp)以及(tp,sp)到Ap中。此时Np={A,B,C,D,sp,tp},此时用户若可以从虚拟起点sp出发经过中间节点顺利到达虚拟tp,那么也可以从tp顺利返回sp;
第二步:对于p上的任一节点i(除sp和tp),若d(sp,i)≤R/2,将sp与i连接,添加(sp,i)到Ap中;
第三步:对于p上的任一节点i(除sp和tp),若d(i,tp)≤R/2,将i与tp连接。(i,tp)添加到Ap中;
第四步:对于p上的任意两个节点i,j,若d(i,j)≤R,并且i的索引<j的索引,将i与j连接。(i,j)添加到Ap中。
最终得到如图3所示的路径扩展图。令G=∪p∈P Gp,路径轨迹与选择组合一一对应,如sp-B-D-tp对应有效选址组合{B,D}。
基于上述路径扩展,为了找到最优的路径以建站,我们定义如下三类决策变量:
选址决策:
模型中的候选选址点集合为F,集合F中的任意一个节点均需要决定是否选址,令yi=1表示在i∈F处建设换电站,否则,yi=0。
路径选择决策:
该变量可将路径选择与选址变量联系起来,当该路径上点i,j处建站时,用户才有被该弧截获的可能。
服务单个用户收益决策:
基于扩展网络,目标函数(P1)可表示:
本模型的约束条件包括流量平衡约束、节点允许流量通过约束、路径选择数量约束、服务单个用户收益约束和自变量定义域约束,具体如下:
流量平衡约束:
公式(9)表示扩展路径p上从虚拟起点发出的流量等于从虚拟终点流入虚拟起点的流量;
公式(10)表示扩展路径p上流入虚拟终点的流量等于从虚拟终点流入虚拟起点的流量;
公式(11)表示对于扩展路径p上的中间节点而言,该节点的流入量等于流出量。
节点允许流量通过约束:
只有在某一点处建站,该点才能有流量通过:
路径选择数量约束:
每个OD对间可服务需求的路径,最多只能选择一条,即若某个OD对q上的需求能够被服务,那么所有用户都选择同一条路径:
服务单个用户收益约束:
自变量定义域约束:
选址变量是0-1变量:
弧变量是0-1变量:
OD对间服务单个用户收益非负:
此时,将P1重构为一个混合整数非线性规划模型。
OA算法是基于切平面算法的经典方法,用于求解凸混合整数非线性规划问题。OA算法的思想是不断用一系列线性函数近似该非线性部分。OA算法构建了一个主问题,主问题是原始非线性规划问题的线性近似,通过求解主问题可得到原问题的一个上界可行解,将该可行解代入公式(8)中可得到原问题最优值的一个下界。每次求解主问题都会添加一个约束,通过数次迭代之后即可得到最优解。
(MPm)maxx,y,z,tμTz-òt-cTy
s.t.(9)-(17)
t≥0 (19)
该算法需要一个初始可行解(x0,y0,z0),可通过启发式算法或者求解ò=0时的(P2)得到。
在该算法中,预设TRG是相对间隙允许值,MaxIt是最大迭代次数,令vm为MPm的最优值。每一次迭代首先求主问题(MPm),最优值vm用来更新上界UB,第m次迭代得到的最优解(xm,ym,zm,tm)用来更新下界LB。由于(xm,ym,zm)是(P2)的一个可行解,因此,μTzm-g(zm)-cTym是(P2)最优值的一个下界。当达到迭代次数或者RG≤TRG时,停止迭代;否则,根据zm添加约束令m=m+1,继续下一次迭代。如图4所示,具体的算法流程如下:
1.初始化:m=1,上界UB=+∞,下界LB=μTz0-g(z0)-cTy0,相对间隙RG=+∞。
2.求解主问题(MPm)得到最优解(xm,ym,zm,tm)和最优值vm。
3.若vm<UB,更新上界UB=vm;
更新下界LB=max{LB,μTzm-g(zm)-cTy0};
4.若m达到最大迭代次数MaxIt或相对间隙RG小于等于相对间隙允许值TRG,输出最优解(x*,y*,z*,t*)=(xm,ym,zm,tm),上界UB、下界LB和相对间隙RG;否则,添加线性约束令m=m+1,返回步骤2。
以上仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法,其特征在于,包括:
获取运输网络所对应的OD对集合和选址点集合;
根据所述OD对集合生成路径,并根据线性衰减函数表示路径需求截获比例;
利用不确定分布函数集合建立分布鲁棒-偏移流量续航选址DR-DFRL模型,具体包括:
所述OD对上的换电需求构成向量,根据所述换电需求的均值和协方差矩阵得到所述换电需求对应的概率密度函数所属的分布函数集合;
利用所述分布函数集合构建出目标函数,以寻求最坏情况下平均利润最大化的最优选址集合;
根据拉格朗日对偶理论,将所述目标函数转换为所述DR-DFRL模型的目标函数;
通过增加虚拟节点对所述运输网络进行扩展,得到扩展网络;
基于所述扩展网络,以选址决策、路径选择决策和服务单个用户收益决策为决策变量,重构所述DR-DFRL模型,其中,所述选址决策为所述选址点集合中的任意一个节点决定是否选址的决策,所述路径选择决策为用户选择某一路径出行的决策,当该路径上点i,j处建站时,用户才有被该路径上的弧(i,j)截获的可能;
其中,利用线性衰减函数表达所述OD对间截获的路径需求截获比例,根据所述路径需求截获比例与所述OD对的随机需求得到所述OD对间截获的换电需求量,并根据换电服务费用、车辆行驶里程和路径选择决策变量,得到所述运输网络上的服务单个用户收益决策变量;
利用外部估计算法求解所述DR-DFRL模型,求解过程中迭代加入约束条件,得到符合约束条件的决策变量最优解。
2.根据权利要求1所述的考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法,其特征在于,所述约束条件包括流量平衡约束、节点允许流量通过约束、路径选择数量约束、服务单个用户收益约束和自变量定义域约束。
3.根据权利要求1所述的考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法,其特征在于,采用K-最短路算法获得符合所述用户绕行偏移距离的路径,并将不符合所述用户绕行偏移距离的路径及相应的节点删除以获得简化后的所述运输网络。
4.根据权利要求1所述的考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法,其特征在于,通过在所述运输网络中增加虚拟起点和虚拟终点,基于换电站建站规则进行路径扩展,构成所述扩展网络。
5.根据权利要求2所述的考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法,其特征在于,利用外部估计算法对所述DR-DFRL模型进行求解过程中,通过求解原始非线性规划问题的线性近似主问题得到原始问题的上界可行解,将所述上界可行解代入所述DR-DFRL模型的目标函数,得到原始问题的最优值的一个下界值,逐一迭代加入所述约束条件,得到符合所有约束条件的最优解。
6.根据权利要求5所述的考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法,其特征在于,在所述外部估计算法中,在迭代次数达到预设的最大迭代次数或相对间隙小于等于预设的相对间隙允许值时,停止迭代。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010383239.7A CN111582581B (zh) | 2020-05-08 | 2020-05-08 | 考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010383239.7A CN111582581B (zh) | 2020-05-08 | 2020-05-08 | 考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111582581A CN111582581A (zh) | 2020-08-25 |
CN111582581B true CN111582581B (zh) | 2022-06-21 |
Family
ID=72125457
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010383239.7A Active CN111582581B (zh) | 2020-05-08 | 2020-05-08 | 考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111582581B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112581313B (zh) * | 2020-12-23 | 2022-02-22 | 北京理工大学 | 一种光伏充电站资源的分配调节方法和系统 |
CN112613659B (zh) * | 2020-12-23 | 2022-11-15 | 东南大学 | 基于三维时空网络的新能源货运充电站选址-路径优化方法 |
CN112784204B (zh) * | 2020-12-28 | 2024-06-25 | 北京交通大学 | 面向不确定需求的列车时刻表与客流控制鲁棒优化方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105760960A (zh) * | 2016-02-29 | 2016-07-13 | 东南大学 | 基于轨道交通的停车换乘设施最优选址和容量确定方法 |
WO2017045294A1 (zh) * | 2015-09-17 | 2017-03-23 | 华南理工大学 | 一种城市常规公交线网设计方法 |
CN108665089A (zh) * | 2018-04-04 | 2018-10-16 | 清华大学 | 一种用于选址问题的鲁棒优化模型求解方法 |
CN109034465A (zh) * | 2018-07-16 | 2018-12-18 | 河海大学 | 考虑充电站选址与出行路径耦合的充电站两层规划方法 |
-
2020
- 2020-05-08 CN CN202010383239.7A patent/CN111582581B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017045294A1 (zh) * | 2015-09-17 | 2017-03-23 | 华南理工大学 | 一种城市常规公交线网设计方法 |
CN105760960A (zh) * | 2016-02-29 | 2016-07-13 | 东南大学 | 基于轨道交通的停车换乘设施最优选址和容量确定方法 |
CN108665089A (zh) * | 2018-04-04 | 2018-10-16 | 清华大学 | 一种用于选址问题的鲁棒优化模型求解方法 |
CN109034465A (zh) * | 2018-07-16 | 2018-12-18 | 河海大学 | 考虑充电站选址与出行路径耦合的充电站两层规划方法 |
Non-Patent Citations (2)
Title |
---|
"Distributionally Robust Deviation Flow Refueling Location for Electric Vehicle Swapping Stations";Yuli Zhang,Yue Guo,Ningwei Zhang;《2020 Chinese Control And Decision Conference(CCDC)》;20200811;全文 * |
"考虑行程时间可靠性的城市冷链末端配送站选址研究";尹小庆 等;《交通运输系统工程与信息》;20191231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111582581A (zh) | 2020-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111582581B (zh) | 考虑随机交通流量偏移的电动汽车换电站鲁棒选址方法 | |
CN109034465B (zh) | 考虑充电站选址与出行路径耦合的充电站两层规划方法 | |
Pal et al. | Placement of public fast-charging station and solar distributed generation with battery energy storage in distribution network considering uncertainties and traffic congestion | |
CN107180274B (zh) | 一种电动汽车充电设施规划典型场景选取和优化方法 | |
CN115547052B (zh) | 改进自适应大邻域算法的动态需求响应电动公交调度方法 | |
CN112686441B (zh) | 基于交通均衡的电动汽车充电导航及路径选择方法 | |
Bayram et al. | Location analysis of electric vehicle charging stations for maximum capacity and coverage | |
CN110826841A (zh) | 考虑用户充电体验和配电网运行风险的充电站规划方法 | |
CN110189025A (zh) | 考虑不同负荷增长的电动汽车充电站规划方案获取方法 | |
CN112418514A (zh) | 一种运用蚁群系统优化校园公交线路规划的方法 | |
Zhang et al. | Deploying public charging stations for battery electric vehicles on the expressway network based on dynamic charging demand | |
Zhang et al. | An immune algorithm for solving the optimization problem of locating the battery swapping stations | |
CN115879657A (zh) | 考虑多站点容量设计的电动车换电站选址路径优化方法 | |
Bayram et al. | Could petrol stations play a key role in transportation electrification? a gis-based coverage maximization of fast ev chargers in urban environment | |
Shen et al. | Integrated optimization of electric vehicles charging location and allocation for valet charging service | |
Liu et al. | Cooperative Game-based Charging-Discharging Efficiency Optimization of Electric Vehicles in 6G-enabled V2G | |
Cipriani et al. | A multimodal transit network design procedure for urban areas | |
Pourvaziri et al. | Planning of electric vehicle charging stations: An integrated deep learning and queueing theory approach | |
CN111651899A (zh) | 考虑用户选择行为的换电站鲁棒选址定容方法和系统 | |
CN117091616A (zh) | 一种电动汽车充电导航与调度方法 | |
Zhang et al. | A second order cone programming model for pev fast-charging station planning | |
Al-Alwash et al. | Optimal charging scheme for electric vehicles (EVs) in SDN-based vehicular network, scheduling, and routing | |
Shi et al. | The economic impact of electric vehicle routing and charging strategy on traffic-power integrated networks | |
Li et al. | Cross-entropy-based approach to multi-objective electric vehicle charging infrastructure planning | |
Huang et al. | Electric vehicle charging strategy for intercity travel: Impact of user perception and battery degradation |
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 |