CN111160639B - 基于用户出行时空分布特性的电动汽车充电负荷预测方法 - Google Patents
基于用户出行时空分布特性的电动汽车充电负荷预测方法 Download PDFInfo
- Publication number
- CN111160639B CN111160639B CN201911332196.3A CN201911332196A CN111160639B CN 111160639 B CN111160639 B CN 111160639B CN 201911332196 A CN201911332196 A CN 201911332196A CN 111160639 B CN111160639 B CN 111160639B
- Authority
- CN
- China
- Prior art keywords
- charging
- road
- formula
- user
- area
- 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
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 6
- 238000000611 regression analysis Methods 0.000 claims abstract description 4
- 230000002354 daily effect Effects 0.000 claims description 15
- 150000001875 compounds Chemical class 0.000 claims description 12
- 230000003203 everyday effect Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 230000002457 bidirectional effect Effects 0.000 claims description 3
- 230000004927 fusion Effects 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 238000005457 optimization Methods 0.000 abstract 1
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 229910002092 carbon dioxide Inorganic materials 0.000 description 2
- 239000001569 carbon dioxide Substances 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000006467 substitution reaction 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/047—Optimisation of routes or paths, e.g. travelling salesman problem
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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"
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Computational Mathematics (AREA)
- Game Theory and Decision Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Probability & Statistics with Applications (AREA)
- Life Sciences & Earth Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
Abstract
本发明公开了一种基于用户出行时空分布特性的电动汽车充电负荷预测方法。首先利用Bass回归分析模型,预测了规划区电动汽车保有量。构建了融合道路类型、实时道路拥堵程度和区域连接关系的等效道路模型,并考虑不同类型日、不同出行时间、驻留时间和实时路况对用户出行的影响,采用Dijkstra路径寻优算法得到了用户出行最优路径和日出行链。此外,采用蒙特卡洛法得到多次充电负荷预测的结果,并利用非参数核密度估计和高斯拟合得到对应的充电负荷概率密度函数,概率最高的充电负荷即为对应的电动汽车充电负荷预测值。本发明能够有效的描述规划区内日充电负荷的具体分布范围,并取最大概率密度的充电负荷作为预测结果更接近实际。
Description
技术领域
本发明属于电动汽车充电负荷预测技术领域,具体涉及基于用户出行时空分布特性的电动汽车充电负荷预测方法。
背景技术
随着人民生活水平的提高,人们采用汽车出行的越来越多,而汽车保有量的不断提升将导致社会节能环保问题日益凸显。经调查研究显示,柴油汽车所排放的二氧化碳约占城市二氧化碳总排放量的40%。而电动汽车因具有节能环保的优势,能够有效缓解传统能源资源供应紧张和环境污染问题,已成为汽车发展的主流方向。
电动汽车具有较高的可移动性,在行驶过程中可在途经的任何充电站进行充电,因此电动汽车的充电负荷具有地理位置和时间尺度上的分散性和不确定性。同时,用户出行习惯也会影响电动汽车的充电负荷,因此充电负荷也会具有较强的用户主观性。此外,道路的拥堵情况和充电站所处的位置也都会影响不同区域的电动汽车充电负荷。这些因素都对规划区内电动汽车的充电负荷预测带来了很大的困难,因此如何准确描述电动汽车充电负荷是值得研究的一个问题。
发明内容
本发明针对现有技术的不足,提出了一种基于用户出行时空分布特性的电动汽车充电负荷预测方法。
由于电动汽车用户出行的时间、在停驻点驻留的时长等时间分布特性在工作日和节假日所有不同。因此本发明对工作日和节假日分开进行充电负荷预测。负荷预测时间长度为一天,时间间隔为15分钟。用户于工作日期间,一天内首次出行的时间介于早七点至八点半,出行区域多为工作区和教学区。而且用户出行次数较少,多为两至三次,而在节假日期间,首次出行时间介于上午十点至傍晚五点,相对工作日而言,出行的目的地多为商业区和其他区,并且节假日用户日出行次数较多,出行决策次数大多为三至四次。此外采用蒙特卡罗法得到多次充电负荷预测的结果,利用非参数核密度估计和高斯拟合来确定充电负荷概率密度函数,将出现概率最高的充电负荷作为对应的电动汽车充电负荷预测值。
具体按照以下步骤实施:
步骤1、建立Bass回归分析模型对未来N年电动汽车保有量进行预估;
步骤2、根据私家车出行目的地的性质不同,将规划区分为五大类,分别为住宅区、教学区,办公区,商业区和其他区域,并将住宅区视为用户默认出发地;
将所有功能区的几何中心作为功能区的节点,各个功能区域之间以节点相连,为了更好描述各个功能区之间的连接关系,将节点之间的连接线段定义为规划区的道路;
根据国家建委颁发的《城市规划定额指标暂行规定》,将国内的道路类型分为4个等级,各个道路等级rij如式(1)所示;
式中,ri,j∈R,ri,j为节点i,j之间的道路等级;
利用图论方法提取道路交通网络中的道路长度、道路连接关系、道路拥堵系数这些拓扑信息构建多源信息融合的道路模型,具体描述如式(2);
假设规划区中所有的连接道路均为双向的道路,因此ui,j如式(3)所示;
式中,i,j分别为节点;di,j为区域节点i,j之间的距离;inf为无穷大;在不融合多源信息的情况下,初始di,j由式(4)表示;
关于城市道路的拥堵程度用区域之间的平均行驶时间和自由行驶时间确定的道路拥堵系数来表达,具体函数关系式如式(5);
式中,Ci,j(t)为t时刻节点i和j之间道路的拥堵系数,系数值越大表示路段越拥挤,Ti,j(t)表示电动汽车在t时刻节点i和j之间道路的平均行驶时间;T0为每条道路的自由行驶时间,且均为定值;
由以上道路长度,路网拥堵程度和道路等级加权以后即得到的等效道路长度,等效道路长度Di,j(t)表示为:
Di,j(t)=di,j+di,jCi,j(t)+di,j/ri,j (6)
步骤3、利用Dijkstra算法寻找等效道路中的最短路径,并将此寻优的路径作为用户出行的出行链;
定义用户出行链大体分为简单链和复杂链,根据实际情况,用户出发的始发地和日出行终点均为住宅区,因此用户的出行路径构成一个闭环;用户一天内两次的出行路径为简单链,超过两次的路径记为复杂链;
用户在一天的行程中会前往一个或者多个功能区活动,因此用户日出行空间特性由一系列中间驻留点构成的出行链来表示,如式(7)所示;
Q={q1(x1,y1),q2(x2,y2),...,qs(xs,ys),...} (7)
式中,Q为出行链对应驻留功能区的集合;s为驻留功能区的序号;q1为出行链出发点;qs为出行过程中的驻留点,(xs,ys)为功能区几何中心的坐标位置;
出行链包含的路径集合由式(8)表示;
Π={ψ(q1,q2),ψ(q2,q3),...,ψ(qs,qs+1),...} (8)
设pa为两个连续驻留点间的最短路径,pa∈ψ(qs,qs+1)为一个出行链对应的一条路径,Π为出行链对应的路径集合;
步骤4、根据历史数据设定概率分布,得到每台电动汽车不同的出行特征量,此特征量包括首次出发时刻、电动汽车初始蓄电池荷电状态、电动汽车功能区驻留时长和日出行次数;
假设每个用户的出发时间tl 1近似为正态分布,其概率密度函数如式(9)所示;
式中,σD为用户出行用车的始发时间标准差;μD为期望值;
每个区域电动汽车的驻留时间也通过正态分布来描述;概率密度关系如式(10)表示;
用户每日在出行时电池的荷电状态并不全是满充的状态,假设用户每日首次用车时的电池电量服从正态分布,具体表达式由式(11)表示;
忽略出行链内途经的所有区域的驻留时长,用户日出行次数满足以下概率密度函数;
式中,I为用户日出行次数;σI和μI分别为出行次数标准差和期望值;
步骤5、电动汽车在行至目的地的过程中,途经所有功能区内的充电站对电动汽车此时的荷电状态进行评估;通过充电模式选择条件判断电动汽车是否在此充电站选择充电服务;若需要充电则进行步骤六,否则执行步骤七;
为了保证电动汽车能够顺利行驶至目的地,用户在途中会监测电池电量,而电池剩余电量通过实时荷电状态反映,从而判断在哪个功能区进行充电;EVSOCn电动汽车到达下一个功能区时蓄电池荷电状态通过式(13)表示;
式中,S100为电动汽车行驶每百公里电池损耗量;EEV为其蓄电池的额定容量;将第n台电动汽车荷电状态与到达功能区的最小荷电状态SOCMINarrive相比较,判断电动汽车在第i个区域是否需要选择充电模式,具体如下;
用户出行最终回到住宅区后,对电动汽车蓄电池荷电状态采用式(15)的判断条件确定是否充电;
步骤6、电动汽车在此功能区充电站进行充电,待充电至最大荷电状态时充电结束,并驶离此功能区;
充电时期存在慢充Pslow和快充Pfast的区别,充电期间电动汽车蓄电池荷电状态的变化如式(16)表示;
步骤7、行驶至出行链中的下一功能区,跳回循环执行步骤五,直至电动汽车顺利行驶至此次出行的目的地;
步骤8、统计规划区内各功能区在每个时刻的充电负荷量、各功能区日充电负荷总量和规划区日充电负荷总量;
步骤9、利用蒙特卡洛法对规划区充电负荷进行多次预测,采用非参数核密度估计算法对预测结果进行概率密度计算,在此基础上采用高斯拟合得到概率密度函数关系式;
步骤10:将步骤9中得到的工作日和休息日充电负荷概率密度函数关系式中对应出现概率最高的充电负荷值分布作为工作日和休息日的预测充电负荷值。
本发明方法具有的优点及有益结果为:
1)本发明为了更为准确的预测出电动汽车充电负荷,建立了基于电动汽车用户出行的时空分布特性的充电负荷预测模型,考虑用户出行的时间特性根据类型日的不同会发生改变,将工作日和休息日进行分开预测。
2)建立了考虑道路长度,路网拥堵程度和道路等级加权以后即得到的等效道路长度用于描述用户出行的空间特性。并考虑用户出行、驻留的时间特性,给出了基于用户出行时空分布特性的充电负荷预测方法,使得预测的结果更接近实际。
3)本发明为了避免单次预测中存在的随机性和不确定性,采用蒙特卡洛法得到多次充电负荷预测的结果,并利用非参数核密度估计和高斯拟合得到充电负荷概率密度函数。根据得到的概率密度,能够有效的看出规划区内日充电负荷的具体分布范围,并取最大概率密度的充电负荷作为规划区日预测充电负荷,使得充电负荷的预测更具有代表性。
附图说明
图1是本发明规划区具体分区示意图;
图2是本发明出行链的结构示意图;
图3是本发明日充电负荷的流程图。
具体实施方式
下面结合具体实施例对本发明作进一步说明,但不应该理解为本发明上述主体范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
步骤1、建立Bass回归分析模型对未来N年电动汽车保有量进行预估。
步骤2、根据私家车出行目的地的性质不同,将规划区分为五大类,分别为住宅区(home,H)、教学区(Teaching area,T),办公区(Workspace,W),商业区(Businessdistrict,B)和其他区域(Other areas,O),并将住宅区视为用户默认出发地,具体区域示意图如图1所示。
将所有功能区的几何中心作为功能区的节点,各个功能区域之间以节点相连,为了更好描述各个功能区之间的连接关系,将节点之间的连接线段定义为规划区的道路。
根据国家建委颁发的《城市规划定额指标暂行规定》,将国内的道路类型分为4个等级,各个道路等级rij如式(1)所示。
式中,ri,j∈R,ri,j为节点i,j之间的道路等级,规划区内各个功能区之间连接道路的道路等级如下。
为了更准确地描述城市道路交通对私家车出行行为的影响,本发明利用图论方法提取道路交通网络中的道路长度、道路连接关系、道路拥堵系数等拓扑信息构建多源信息融合的道路模型,具体描述如式(3)。
式中,V(G)为规划区中的节点集合;n为道路节点的个数;E(G)为规划区中的道路路段集合;两节点之间的距离即为道路距离ui,j;为节点之间的连接关系矩阵,若区域之间不存在道路连接关系则用inf表示,表达式可由式(4)表示。
假设规划区中所有的连接道路均为双向的道路,因此ui,j如式(5)所示。
式中,i,j分别为节点;di,j为区域节点i,j之间的距离;inf为无穷大。在不融合多源信息的情况下,初始di,j可以由式(6)表示。
关于城市道路的拥堵程度可以用区域之间的平均行驶时间和自由行驶时间确定的道路拥堵系数来表达,具体函数关系式如式(7)。
式中,Ci,j(t)为t时刻节点i和j之间道路的拥堵系数,系数值越大表示路段越拥挤,Ti,j(t)表示电动汽车在t时刻节点i和j之间道路的平均行驶时间;T0为每条道路的自由行驶时间,且均为定值。
由以上道路长度,路网拥堵程度和道路等级加权以后即得到的等效道路长度,等效道路长度Di,j(t)可以表示为:
Di,j(t)=di,j+di,jCi,j(t)+di,j/ri,j (8)
步骤3、利用Dijkstra算法寻找加权后道路的最短路径,并将此寻优的路径作为用户出行的出行链。
本发明定义用户出行链大体可分为简单链和复杂链,根据实际情况,用户出发的始发地和日出行终点均为住宅区,因此用户的出行路径构成一个闭环。用户一天内两次的出行路径为简单链,超过两次的路径记为复杂链。
用户在一天的行程中会前往一个或者多个功能区活动,因此用户日出行空间特性可以由一系列中间驻留点构成的出行链来表示,如式(9)所示。
Q={q1(x1,y1),q2(x2,y2),...,qs(xs,ys),...} (9)
式中,Q为出行链对应驻留功能区的集合;s为驻留功能区的序号;q1为出行链出发点;qs为出行过程中的驻留点,(xs,ys)为功能区几何中心的坐标位置。
出行链包含的路径集合可由式(10)表示。
Π={ψ(q1,q2),ψ(q2,q3),...,ψ(qs,qs+1),...} (10)
设pa为两个连续驻留点间的最短路径,pa∈ψ(qs,qs+1)为一个出行链对应的一条路径,Π为出行链对应的路径集合。出行链的路径结构示意图如图2。
图3中,和为电动汽车离开第i个驻留点的时间以及电动汽车在该驻留点停驻的时间;其中i=1,2,...,s,...,和分别为电动汽车驶离住宅区的时刻和返还住宅区的驻留时间;△Ti,i+1和di,i+1分别为第i,i+1个驻留点之间的行驶时间和行驶里程。
步骤4、根据确定的概率分布,得到每台电动汽车不同的出行特征量,此特征量包括首次出发时刻、电动汽车初始蓄电池荷电状态、电动汽车功能区驻留时长和日出行次数。
式中,σD为用户出行用车的始发时间标准差;μD为期望值。
用户出行时因为工作、学习或购物等个人原因会在出行的目的地驻留一段时间,而停驻的时间因人而异,简化起见,每个区域电动汽车的驻留时间也可通过正态分布来描述。概率密度关系如式(12)表示。
用户每日在出行时电池的荷电状态并不全是满充的状态,本发明中假设用户每日首次用车时的电池电量服从正态分布,具体表达式可由式(13)表示。
忽略出行链内途经的所有区域的驻留时长,用户日出行次数满足以下概率密度函数。
式中,I(x)为用户日出行次数;σI和μI分别为出行次数标准差和期望值。
步骤5、电动汽车在行至目的地的过程中,途经所有功能区内的充电站对电动汽车此时的荷电状态进行评估。通过充电模式选择条件判断电动汽车是否在此充电站选择充电服务。若需要充电则进行步骤六,否则执行步骤七。
为了保证电动汽车能够顺利行驶至目的地,用户在途中会监测电池电量,而电池剩余电量通过实时荷电状态反映,从而判断在哪个功能区进行充电。EVSOCn电动汽车到达下一个功能区时蓄电池荷电状态可通过式(15)表示。
式中,S100为电动汽车行驶每百公里电池损耗量;EEV为其蓄电池的额定容量。
用户出行最终回到住宅区后,对电动汽车蓄电池荷电状态采用式(17)的判断条件确定是否充电。
步骤6、电动汽车在此功能区充电站进行充电,待充电至最大荷电状态时充电结束,并驶离此功能区。
充电时期存在慢充Pslow和快充Pfast的区别,充电期间电动汽车蓄电池荷电状态的变化如式(18)表示。
步骤7、行驶至出行链中的下一功能区,跳回循环执行步骤五,直至电动汽车顺利行驶至此次出行的目的地。
步骤8、统计规划区内各功能区在每个时刻的充电负荷量、各功能区日充电负荷总量和规划区日充电负荷总量。
步骤9、利用蒙特卡洛法对规划区充电负荷进行多次预测,采用非参数核密度估计算法对预测结果进行概率密度计算,在此基础上采用高斯拟合得到概率密度函数关系式。
步骤10:将步骤9中得到的工作日和休息日充电负荷概率密度函数关系式中对应出现概率最高的充电负荷值分布作为工作日和休息日的预测充电负荷值。
Claims (1)
1.基于用户出行时空分布特性的电动汽车充电负荷预测方法,其特征在于,该方法具体包括以下步骤:
步骤1、建立Bass回归分析模型对未来N年电动汽车保有量进行预估;
步骤2、根据私家车出行目的地的性质不同,将规划区分为五大类,分别为住宅区、教学区,办公区,商业区和其他区域,并将住宅区视为用户默认出发地;
将所有功能区的几何中心作为功能区的节点,各个功能区域之间以节点相连,为了更好描述各个功能区之间的连接关系,将节点之间的连接线段定义为规划区的道路;
将国内的道路类型分为4个等级,各个道路等级rij如式(1)所示;
式中,ri,j∈R,ri,j为节点i,j之间的道路等级;
利用图论方法提取道路交通网络中的道路长度、道路连接关系、道路拥堵系数这些拓扑信息构建多源信息融合的道路模型,具体描述如式(2);
假设规划区中所有的连接道路均为双向的道路,因此ui,j如式(3)所示;
式中,i,j分别为节点;di,j为区域节点i,j之间的距离;inf为无穷大;在不融合多源信息的情况下,初始di,j由式(4)表示;
关于城市道路的拥堵程度用区域之间的平均行驶时间和自由行驶时间确定的道路拥堵系数来表达,具体函数关系式如式(5);
式中,Ci,j(t)为t时刻节点i和j之间道路的拥堵系数,系数值越大表示路段越拥挤,Ti,j(t)表示电动汽车在t时刻节点i和j之间道路的平均行驶时间;T0为每条道路的自由行驶时间,且均为定值;
由以上道路长度,路网拥堵程度和道路等级加权以后即得到的等效道路长度,等效道路长度Di,j(t)表示为:
Di,j(t)=di,j+di,jCi,j(t)+di,j/ri,j (6)
步骤3、利用Dijkstra算法寻找等效道路中的最短路径,并将此寻优的路径作为用户出行的出行链;
定义用户出行链分为简单链和复杂链,根据实际情况,用户出发的始发地和日出行终点均为住宅区,因此用户的出行路径构成一个闭环;用户一天内两次的出行路径为简单链,超过两次的路径记为复杂链;
用户在一天的行程中会前往一个或者多个功能区活动,因此用户日出行空间特性由一系列中间驻留点构成的出行链来表示,如式(7)所示;
Q={q1(x1,y1),q2(x2,y2),...,qs(xs,ys),...} (7)
式中,Q为出行链对应驻留功能区的集合;s为驻留功能区的序号;q1为出行链出发点;qs为出行过程中的驻留点,(xs,ys)为功能区几何中心的坐标位置;
出行链包含的路径集合由式(8)表示;
Π={ψ(q1,q2),ψ(q2,q3),...,ψ(qs,qs+1),...} (8)
设pa为两个连续驻留点间的最短路径,pa∈ψ(qs,qs+1)为一个出行链对应的一条路径,Π为出行链对应的路径集合;
步骤4、根据历史数据设定概率分布,得到每台电动汽车不同的出行特征量,此特征量包括首次出发时刻、电动汽车初始蓄电池荷电状态、电动汽车功能区驻留时长和日出行次数;
式中,σD为用户出行用车的始发时间标准差;μD为期望值;
每个区域电动汽车的驻留时间也通过正态分布来描述;概率密度关系如式(10)表示;
用户每日在出行时电池的荷电状态并不全是满充的状态,假设用户每日首次用车时的电池电量服从正态分布,具体表达式由式(11)表示;
忽略出行链内途经的所有区域的驻留时长,用户日出行次数满足以下概率密度函数;
式中,I为用户日出行次数;σI和μI分别为出行次数标准差和期望值;
步骤5、电动汽车在行至目的地的过程中,途经所有功能区内的充电站对电动汽车此时的荷电状态进行评估;通过充电模式选择条件判断电动汽车是否在此充电站选择充电服务;若需要充电则进行步骤6,否则执行步骤7;
为了保证电动汽车能够顺利行驶至目的地,用户在途中会监测电池电量,而电池剩余电量通过实时荷电状态反映,从而判断在哪个功能区进行充电;EVSOCn电动汽车到达下一个功能区时蓄电池荷电状态通过式(13)表示;
式中,S100为电动汽车行驶每百公里电池损耗量;EEV为其蓄电池的额定容量;
用户出行最终回到住宅区后,对电动汽车蓄电池荷电状态采用式(15)的判断条件确定是否充电;
步骤6、电动汽车在此功能区充电站进行充电,待充电至最大荷电状态时充电结束,并驶离此功能区;
充电时期存在慢充Pslow和快充Pfast的区别,充电期间电动汽车蓄电池荷电状态的变化如式(16)表示;
步骤7、行驶至出行链中的下一功能区,跳回循环执行步骤五,直至电动汽车顺利行驶至此次出行的目的地;
步骤8、统计规划区内各功能区在每个时刻的充电负荷量、各功能区日充电负荷总量和规划区日充电负荷总量;
步骤9、利用蒙特卡洛法对规划区充电负荷进行多次预测,采用非参数核密度估计算法对预测结果进行概率密度计算,在此基础上采用高斯拟合得到概率密度函数关系式;
步骤10:将步骤9中得到的工作日和休息日充电负荷概率密度函数关系式中对应出现概率最高的充电负荷值分布作为工作日和休息日的预测充电负荷值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911332196.3A CN111160639B (zh) | 2019-12-21 | 2019-12-21 | 基于用户出行时空分布特性的电动汽车充电负荷预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911332196.3A CN111160639B (zh) | 2019-12-21 | 2019-12-21 | 基于用户出行时空分布特性的电动汽车充电负荷预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111160639A CN111160639A (zh) | 2020-05-15 |
CN111160639B true CN111160639B (zh) | 2022-04-15 |
Family
ID=70557680
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911332196.3A Active CN111160639B (zh) | 2019-12-21 | 2019-12-21 | 基于用户出行时空分布特性的电动汽车充电负荷预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111160639B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111815017B (zh) * | 2020-05-29 | 2021-07-30 | 国网山东省电力公司经济技术研究院 | 一种基于出行数据的电动汽车充电负荷预测方法 |
CN112257907B (zh) * | 2020-09-30 | 2023-04-07 | 广东电网有限责任公司广州供电局 | 一种考虑电价敏感性的电动汽车负荷建模方法 |
CN113268709A (zh) * | 2021-04-21 | 2021-08-17 | 国网天津市电力公司 | 基于智能体仿真的城市电动汽车充电需求预测方法及系统 |
CN113821911B (zh) * | 2021-08-24 | 2022-05-17 | 中国能源建设集团广东省电力设计研究院有限公司 | 一种充电站负荷的预测方法、装置及系统 |
CN114049007B (zh) * | 2021-11-15 | 2024-07-12 | 东南大学 | 一种城市电网对电动汽车承载能力的计算方法 |
CN117648520B (zh) * | 2024-01-29 | 2024-05-07 | 北京理工大学 | 基于新能源汽车充电负荷分析的碳排放计算方法及系统 |
CN117833240B (zh) * | 2024-02-29 | 2024-05-31 | 江苏米特物联网科技有限公司 | 一种面向酒店场景的电动汽车充电负荷预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103903090A (zh) * | 2014-03-17 | 2014-07-02 | 东南大学 | 基于用户意愿和出行规律的电动汽车充电负荷分配方法 |
CN105719030A (zh) * | 2016-03-29 | 2016-06-29 | 武汉大学 | 一种基于效能最大化原则的电动汽车负荷预测方法 |
WO2016187902A1 (zh) * | 2015-05-28 | 2016-12-01 | 深圳市华宝新能源股份有限公司 | 电动汽车移动充电控制方法和系统 |
CN108510128A (zh) * | 2018-04-11 | 2018-09-07 | 华南理工大学广州学院 | 一种区域电动汽车充电负荷时空分布预测方法 |
-
2019
- 2019-12-21 CN CN201911332196.3A patent/CN111160639B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103903090A (zh) * | 2014-03-17 | 2014-07-02 | 东南大学 | 基于用户意愿和出行规律的电动汽车充电负荷分配方法 |
WO2016187902A1 (zh) * | 2015-05-28 | 2016-12-01 | 深圳市华宝新能源股份有限公司 | 电动汽车移动充电控制方法和系统 |
CN105719030A (zh) * | 2016-03-29 | 2016-06-29 | 武汉大学 | 一种基于效能最大化原则的电动汽车负荷预测方法 |
CN108510128A (zh) * | 2018-04-11 | 2018-09-07 | 华南理工大学广州学院 | 一种区域电动汽车充电负荷时空分布预测方法 |
Non-Patent Citations (3)
Title |
---|
The Multiobjective Based Large-Scale Electric Vehicle Charging;Yimin Zhou 等;《Hindawi Complexity》;20181016;全文 * |
电动汽车广泛接入对电网的影响及其调控策略研究;陈丽丹;《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅱ辑》;20190515;全文 * |
考虑用户出行特性的电动汽车充电负荷预测;罗平;《华中科技大学学报》;20201130;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111160639A (zh) | 2020-05-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111160639B (zh) | 基于用户出行时空分布特性的电动汽车充电负荷预测方法 | |
CN107392400B (zh) | 计及实时交通与温度的ev充电负荷时空分布预测方法 | |
CN109034648B (zh) | 一种电动汽车集群需求响应潜力评估方法 | |
CN110728396B (zh) | 一种考虑时空分布的电动汽车充电负荷综合建模方法 | |
CN106599390B (zh) | 一种计及电动出租车时空随机特性的充电负荷的计算方法 | |
CN109693576B (zh) | 一种基于模拟退火算法的电动汽车充电调度优化方法 | |
CN108596667B (zh) | 一种基于车联网的电动汽车实时充电电价计算方法 | |
CN111400662B (zh) | 一种考虑电动汽车充电需求的空间负荷预测方法 | |
CN103236179A (zh) | 一种计及交通信息与电网信息的电动汽车充电导航方法 | |
CN105160428A (zh) | 高速公路电动汽车快速充电站的规划方法 | |
CN108133329B (zh) | 考虑充电反馈效应的电动汽车出行与充电需求分析方法 | |
CN111397620A (zh) | 快充/慢充模式下的电动汽车充电导航方法及系统 | |
CN112381398B (zh) | 考虑用户有限理性出行的电动汽车充电站选址方法及系统 | |
CN111199320A (zh) | 基于出行概率矩阵的电动汽车充电负荷时空分布预测方法 | |
CN114021391A (zh) | 基于动态能耗与用户心理的电动汽车充电负荷预测方法 | |
CN116054316A (zh) | 一种基于混沌麻雀优化算法的电动汽车有序充放电方法 | |
CN111199300A (zh) | 一种“车-路-网”模式下电动汽车充电负荷时空预测方法 | |
CN115660501A (zh) | 电动汽车充电负荷可调裕度评估方法 | |
CN113112097A (zh) | 一种电动汽车负荷预测及充电设施布局优化方法 | |
Sundstroem et al. | Toward electric vehicle trip prediction for a charging service provider | |
CN112330025A (zh) | 用于城市电动汽车的时空充电负荷的预测方法 | |
CN109117972A (zh) | 一种电动汽车的充电需求确定方法 | |
CN114840984A (zh) | 一种考虑出行路径决策的电动汽车充电负荷时空预测方法 | |
CN114676885B (zh) | 一种电动汽车充放电负荷时空分布预测方法 | |
Mahmoudi et al. | Spatial-temporal prediction of electric vehicle charging demand in realistic urban transportation system of a mid-sized city in brazil |
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 |