CN110570014A - A method for forecasting electric vehicle charging load based on Monte Carlo and deep learning - Google Patents
A method for forecasting electric vehicle charging load based on Monte Carlo and deep learning Download PDFInfo
- Publication number
- CN110570014A CN110570014A CN201910725799.3A CN201910725799A CN110570014A CN 110570014 A CN110570014 A CN 110570014A CN 201910725799 A CN201910725799 A CN 201910725799A CN 110570014 A CN110570014 A CN 110570014A
- Authority
- CN
- China
- Prior art keywords
- charging
- electric
- load
- time
- initial
- 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
- 238000007600 charging Methods 0.000 title claims abstract description 313
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000013135 deep learning Methods 0.000 title claims abstract description 24
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 230000029305 taxis Effects 0.000 claims abstract description 16
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 8
- 238000012614 Monte-Carlo sampling Methods 0.000 claims abstract description 6
- 230000006870 function Effects 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000011161 development Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 5
- 239000000446 fuel Substances 0.000 claims description 4
- 230000007423 decrease Effects 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- QBPFLULOKWLNNW-UHFFFAOYSA-N chrysazin Chemical compound O=C1C2=CC=CC(O)=C2C(=O)C2=C1C=CC=C2O QBPFLULOKWLNNW-UHFFFAOYSA-N 0.000 claims 1
- 238000010606 normalization Methods 0.000 claims 1
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 description 4
- 229910052744 lithium Inorganic materials 0.000 description 4
- 238000013528 artificial neural network Methods 0.000 description 3
- 238000013277 forecasting method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000003930 cognitive ability Effects 0.000 description 1
- 238000010280 constant potential charging Methods 0.000 description 1
- 238000010277 constant-current charging Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013278 delphi method Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- 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)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Evolutionary Computation (AREA)
- Marketing (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Secondary Cells (AREA)
Abstract
Description
技术领域technical field
本发明属于电力系统领域,具体地说是一种基于蒙特卡洛和深度学习的电动汽车充电负荷预测方法。The invention belongs to the field of power systems, in particular to a method for predicting charging loads of electric vehicles based on Monte Carlo and deep learning.
背景技术Background technique
由于能源安全和环境污染问题的日益突出,近年来新能源获得大力的发展,其中电动汽车更是发展迅速,将逐步替代传统的燃油汽车,成为未来主要的交通工具,具有广阔的发展前景。与其他低碳负荷相比,电动汽车具有规模化应用的基础,将成为电网负荷的重要组成部分。当电动汽车规模化应用时,其对电网的影响不可忽视。因此,充分考虑电动汽车的发展规律和充电规律对城市配电网的长远发展具有重要意义。对其保有量及比例进行预测,充电负荷进行建模,给出未来的日负荷预测结果是开展电动汽车接入对电网的影响分析、配电网规划与控制运行、电动汽车与电网双向互动及电动汽车与其他能源、交通等系统协调研究的基础。Due to the increasingly prominent problems of energy security and environmental pollution, new energy has been vigorously developed in recent years, among which electric vehicles are developing rapidly, and will gradually replace traditional fuel vehicles and become the main means of transportation in the future, with broad development prospects. Compared with other low-carbon loads, electric vehicles have the basis for large-scale application and will become an important part of grid load. When electric vehicles are applied on a large scale, their impact on the power grid cannot be ignored. Therefore, fully considering the development law and charging law of electric vehicles is of great significance to the long-term development of urban distribution network. Predict its inventory and proportion, model the charging load, and give the future daily load forecast results. The basis for coordinated research between electric vehicles and other energy, transportation and other systems.
目前,在电动汽车负荷预测方面,已提出了一些指标和方法,例如常见的负荷预测方法有单耗法、趋势分析法、弹性系数法、回归分析法、时间序列法、灰色模型法、神经网络法、德尔菲法、专家系统法以及优选组合分析法等方法。At present, some indicators and methods have been proposed for electric vehicle load forecasting. For example, common load forecasting methods include unit consumption method, trend analysis method, elastic coefficient method, regression analysis method, time series method, gray model method, and neural network method. method, Delphi method, expert system method, optimal combination analysis method and other methods.
现有电动汽车充电负荷预测研究仍存在不足之处:普适性不强。大部分研究仅考虑了充电行驶里程、充电起始时间、起始荷电状态等建模条件的简单因素,尚无涵盖各种关键因素和多种车辆类型的理论负荷模型。且传统法以及浅层学习法自适应能力不够,对非线性特征负荷认知能力不够,预测精度不高。There are still deficiencies in the existing research on electric vehicle charging load forecasting: the universality is not strong. Most studies only consider simple factors such as charging mileage, charging start time, initial state of charge and other modeling conditions, and there is no theoretical load model covering various key factors and multiple vehicle types. Moreover, the traditional method and the shallow learning method have insufficient adaptive ability, insufficient cognitive ability for nonlinear feature load, and low prediction accuracy.
可见,电动汽车充电负荷预测方法还有待改进。It can be seen that the electric vehicle charging load forecasting method needs to be improved.
发明内容Contents of the invention
为解决上述技术问题,本发明提出一种基于蒙特卡洛和深度学习的电动汽车充电负荷预测方法。该方法建立了负荷因素的影响概率模型,搭建了不同类型电动汽车充电负荷计算模型,采用蒙特卡洛模拟方法得到各时刻充电负荷,然后采用深度学习算法求取电动汽车充电负荷预测曲线。In order to solve the above-mentioned technical problems, the present invention proposes a method for predicting electric vehicle charging load based on Monte Carlo and deep learning. This method establishes the probability model of the influence of load factors, builds the charging load calculation model of different types of electric vehicles, uses the Monte Carlo simulation method to obtain the charging load at each time, and then uses the deep learning algorithm to obtain the electric vehicle charging load prediction curve.
本发明采用以下技术方案实现:The present invention adopts following technical scheme to realize:
一种基于蒙特卡洛和深度学习的电动汽车充电负荷预测方法,包括以下步骤:A method for forecasting electric vehicle charging load based on Monte Carlo and deep learning, comprising the following steps:
首先,根据电动汽车的特征将电动汽车划分为电动公交车、电动出租车、电动私家车和电动公务车4种类型,建立负荷影响因素的概率模型,进而得到不同类型电动汽车充电功率的计算模型。First of all, according to the characteristics of electric vehicles, electric vehicles are divided into four types: electric buses, electric taxis, electric private cars and electric official vehicles, and the probability model of load influencing factors is established, and then the calculation model of charging power of different types of electric vehicles is obtained .
其次,根据电动汽车保有量预测结果,采用蒙特卡洛模拟方法抽取电动汽车的起始荷电状态、起始充电时间等来计算电动汽车的充电负荷,得到各时刻的充电负荷。Secondly, according to the prediction results of electric vehicle ownership, Monte Carlo simulation method is used to extract the initial state of charge and initial charging time of electric vehicles to calculate the charging load of electric vehicles, and obtain the charging load at each moment.
最后,根据蒙特卡洛抽样得到的各时刻电动汽车充电负荷,采用LSTM深度学习算法对电动汽车充电负荷进行预测,从而得到电动汽车充电负荷曲线。Finally, according to the charging load of electric vehicles at each time obtained by Monte Carlo sampling, the LSTM deep learning algorithm is used to predict the charging load of electric vehicles, so as to obtain the charging load curve of electric vehicles.
上述技术方案中,进一步地,步骤1)所述建立负荷影响因素概率模型,得到电动汽车充电功率的计算模型,包括:In the above technical solution, further, step 1) establishes a probability model of load influencing factors to obtain a calculation model of electric vehicle charging power, including:
(1)电动汽车的充电负荷受多方面因素的影响,主要影响因素有起始充电时刻、日行驶里程、充电时长、电动汽车保有量、起始荷电状态(SOC)、电池容量等。(1) The charging load of electric vehicles is affected by many factors, the main influencing factors are initial charging time, daily mileage, charging time, electric vehicle ownership, initial state of charge (SOC), battery capacity, etc.
根据电动汽车的特征将电动汽车划分为电动公交车、电动出租车、电动私家车和电动公务车4种类型,构建起始充电时间模型,通过数据拟合处理得到电动汽车的起始充电时间满足下式所示的正态分布:According to the characteristics of electric vehicles, electric vehicles are divided into four types: electric buses, electric taxis, electric private cars and electric official vehicles, and the initial charging time model is constructed, and the initial charging time of electric vehicles is obtained through data fitting processing The normal distribution given by:
式中:t为初始充电时间,即最后一次出行的结束时刻;μa和σa分别为起始充电时间的期望和标准差,由于四种类型的汽车拟合出来的起始充电时刻的正态分布结果不同,因此不同类型电动汽车的期望和标准差不同;In the formula: t is the initial charging time, that is, the end time of the last trip; μ a and σ a are the expectation and standard deviation of the initial charging time, respectively, since the positive The results of the state distribution are different, so the expectations and standard deviations of different types of electric vehicles are different;
(2)构建日行驶里程模型,日行驶里程是电动汽车驾驶特性的重要指标,反映了汽车在一天内的耗电量,进而影响电动汽车的充电时间。同理,用传统燃油车的出行特性代替电动汽车的出行特性进行分析,电动汽车的日行驶里程服从对数正态分布,其概率密度函数为:(2) Construct a daily mileage model. The daily mileage is an important indicator of the driving characteristics of electric vehicles, reflecting the power consumption of the vehicle in a day, which in turn affects the charging time of electric vehicles. Similarly, using the travel characteristics of traditional fuel vehicles instead of electric vehicles for analysis, the daily mileage of electric vehicles follows a lognormal distribution, and its probability density function is:
式中:s为日行驶里程,单位为km;μb和σb分别为行驶里程s的对数lns的期望和方差,随着不同类型电动汽车行驶特性的不同而变化。In the formula: s is the daily mileage in km; μ b and σ b are the expectation and variance of the logarithm lns of the mileage s, respectively, which vary with the driving characteristics of different types of electric vehicles.
(3)构建起始荷电状态模型,在一个充电周期内,电动汽车充电所需的电力需求是随时间变化的。要确定电动汽车的充电负荷必须要获得电池充电开始时刻的荷电状态,即起始SOC。起始SOC是上次充电后的电动汽车行驶距离的随机函数,取值范围为0到100%。假设电动汽车的荷电状态随着行驶里程是线性下降的,则可通过车辆的行驶里程估算出起始充电时间的荷电状态,用概率密度函数来表示:(3) Construct the initial state-of-charge model. In a charging cycle, the electric power demand for electric vehicle charging changes with time. To determine the charging load of electric vehicles, it is necessary to obtain the state of charge at the beginning of battery charging, that is, the initial SOC. The starting SOC is a random function of the distance traveled by the EV since the last charge, with values ranging from 0 to 100%. Assuming that the state of charge of an electric vehicle decreases linearly with the mileage, the state of charge at the initial charging time can be estimated from the mileage of the vehicle, expressed by a probability density function:
SOC=(SOC0-s/smax)×100%S OC =(S OC0 -s/s max )×100%
式中:SOC表示电池充电的起始SOC;SOC0表示最近一次充电后电池的荷电状态值。由于电动汽车的充电时间和地点更为分散,SOC0往往不为1;smax表示电池充满后可行驶的最大里程数,单位km。In the formula: S OC represents the initial SOC of battery charging; S OC0 represents the state of charge value of the battery after the latest charge. Since the charging time and location of electric vehicles are more scattered, S OC0 is often not 1; s max indicates the maximum mileage that can be driven after the battery is fully charged, in km.
(4)构建充电时长模型,目前,电动汽车大多使用锂电池,锂电池的充电过程是恒压-恒流两阶段充电过程。假设整个充电过程是恒定功率的,则充电时长有以下两种计算方式:(4) Build a charging time model. At present, most electric vehicles use lithium batteries, and the charging process of lithium batteries is a two-stage charging process of constant voltage and constant current. Assuming that the entire charging process is a constant power, the charging time can be calculated in the following two ways:
基于荷电状态计算充电时长,则有如下公式:To calculate the charging time based on the state of charge, there is the following formula:
式中:Te为充电时长,单位为h;U为电池容量,单位为kW·h;P为充电功率,单位为kW;η为充电效率;In the formula: Te is the charging time, the unit is h; U is the battery capacity, the unit is kW h; P is the charging power, the unit is kW; η is the charging efficiency;
基于日行驶里程的方式计算充电时长,则如下式所示:Calculate the charging time based on the daily mileage, as shown in the following formula:
式中:W100为汽车每行驶100km的耗电量,单位为(kW·h)/百公里。In the formula: W 100 is the power consumption of the car every 100km, and the unit is (kW·h)/100km.
更进一步地,步骤2)计算各时刻充电负荷,包括:Furthermore, step 2) calculates the charging load at each moment, including:
(1)划分不同类型电动汽车出行特征:将电动汽车按照不同的用途分为公交车、出租车、私家车和公务车4类,具体分析不同类型电动汽车的出行特征,以得到负荷预测模型的参数,方法如下:(1) Divide the travel characteristics of different types of electric vehicles: divide electric vehicles into four categories according to different purposes: buses, taxis, private cars, and official vehicles, and analyze the travel characteristics of different types of electric vehicles in order to obtain the load forecasting model. parameters, as follows:
电动公交车electric bus
公交车的行驶特性相对非常的固定,主要采取轮班运营制。根据文献《基于蒙特卡洛模拟的电动汽车充电负荷预测》中的数据,公交车的日行驶里程大约为70km。考虑到安全运行,一天一充难以满足电动公交车的运营需求,需要一天两充。公交车的运营时间、路线相对集中,可以进行集中充电。在中午时段进行快速充电,晚上下班后进行常规充电。一般而言,公交车的充电时段为9.30至16.00、23:00至次日05:00,分别服从正态分布,具体分布参数可根据城市调研数据获取。The driving characteristics of the bus are relatively fixed, and the operation system is mainly adopted in shifts. According to the data in the literature "Electric Vehicle Charging Load Prediction Based on Monte Carlo Simulation", the daily mileage of the bus is about 70km. Considering safe operation, charging once a day is difficult to meet the operating needs of electric buses, and charging twice a day is required. The operating hours and routes of buses are relatively concentrated, and centralized charging can be carried out. Fast charging at noon and regular charging at night after get off work. Generally speaking, the charging period of the bus is from 9.30 to 16.00, and from 23:00 to 05:00 the next day, which respectively obey the normal distribution. The specific distribution parameters can be obtained according to the urban research data.
电动出租车electric taxi
出租车的运营时间大致为06:00--24:00。根据文献《北京市私人机动车交通出行特征及发展对策》中的数据,出租车的日行驶里程大约为400km。同电动公交车一样,电动出租车一般采取一天两充的模式,充电时间选择在中午换班及晚上。由于出租车的休息时间有限,但需要及时补充电量,因此电动出租车选择快速充电模式。根据以上分析,出租车的充电时段为02:00—05:00、11:30—14:30,分别服从正态分布。The operating hours of taxis are roughly from 06:00 to 24:00. According to the data in the document "Travel Characteristics and Development Strategies of Private Motor Vehicles in Beijing", the daily mileage of taxis is about 400km. Like electric buses, electric taxis generally adopt the mode of charging twice a day, and the charging time is selected at noon and night. Since the rest time of the taxi is limited, but the battery needs to be replenished in time, the electric taxi chooses the fast charging mode. According to the above analysis, the charging time of taxis is 02:00-05:00 and 11:30-14:30, respectively obeying the normal distribution.
电动私家车electric private car
相较于公交车和出租车,私家车的行驶特性更具有随机性和任意性。电动私家车每日充电一次,充电时段分为上午09:00—12:00、下午14:00--17:00和晚上19:00至次日07:00。上午和下午在工作单位的停车场充电,晚上下班后在居民区的停车场充电。在单位停车场充电选择快速充电模式,在居民区停车场选择常规充电模式。根据文献中的数据,不考虑长途出行,私家车的日行驶里程为40km。Compared with buses and taxis, the driving characteristics of private cars are more random and arbitrary. Electric private cars are charged once a day, and the charging time is divided into 09:00-12:00 in the morning, 14:00-17:00 in the afternoon, and 19:00 in the evening to 07:00 the next day. Charge in the parking lot of the workplace in the morning and afternoon, and charge in the parking lot of the residential area after get off work in the evening. Choose the fast charging mode for charging in the unit parking lot, and choose the regular charging mode for the parking lot in the residential area. According to the data in the literature, regardless of long-distance travel, the daily mileage of private cars is 40km.
电动商务车Electric commercial vehicle
公务车主要用作日常公务出行,若不考虑长途出行,其行驶特性和私家车相似。公务车的充电时间一般在晚上下班后于单位的停车场充电。公务车一天一充即可,且采取常规充电模式,充电时段为19:00至次日07:00。Official vehicles are mainly used for daily business trips. If long-distance travel is not considered, their driving characteristics are similar to private cars. The charging time of official vehicles is generally charged in the parking lot of the unit after get off work at night. Official vehicles can be charged once a day, and the regular charging mode is adopted, and the charging period is from 19:00 to 07:00 the next day.
(2)预测各种类型电动汽车保有量:以《中国汽车产业发展报告(2012)》中指出的2020年的数据为基数,根据不同类型电动汽车的比例,对各种类型的汽车数量进行预测,(2) Predict the number of various types of electric vehicles: Based on the data in 2020 pointed out in the "China Automobile Industry Development Report (2012)", the number of various types of vehicles is predicted according to the proportion of different types of electric vehicles ,
基于蒙特卡洛仿真对电动汽车充电负荷进行计算:Calculate the charging load of electric vehicles based on Monte Carlo simulation:
做合理假设:Make reasonable assumptions:
a.各类型电动汽车的起始充电时间、日行驶里程和充电功率为相互独立的随机变量;a. The initial charging time, daily driving mileage and charging power of various types of electric vehicles are independent random variables;
b.各类型电动汽车的充电功率视为恒功率模型;b. The charging power of various types of electric vehicles is regarded as a constant power model;
c.所有车辆每次都充满电量;c. All vehicles are fully charged every time;
d.所有车辆的最后一次出行结束时刻即为该车辆的起始充电时间。d. The end time of the last trip of all vehicles is the initial charging time of the vehicle.
(3)选择蒙特卡洛方法模拟法对不同类型的电动汽车进行采样,通过抽取起始充电时间和日行驶里程,得到单辆电动汽车的充电时长和充电功率。然后累加所有电动汽车充电功率,得到充电负荷;以时间为横坐标,各时刻的充电负荷为纵坐标,可得到充电负荷曲线。则j时刻的充电负荷为:(3) Choose the Monte Carlo simulation method to sample different types of electric vehicles, and obtain the charging time and charging power of a single electric vehicle by extracting the initial charging time and daily mileage. Then add up the charging power of all electric vehicles to obtain the charging load; take time as the abscissa and the charging load at each moment as the ordinate to obtain the charging load curve. Then the charging load at time j is:
式中:Ptotal,j为在j时刻的充电负荷,单位为kw;Nb,Nt,Ns,Nw分别为电动公交车、出租车、私家车和商务车的总数,单位为辆;将一天划分为T个时间段;Pibj,Pibj,Pibj,Pibj为不同类型的电动汽车在j时刻的充电功率,单位为kW。In the formula: P total, j is the charging load at time j, the unit is kw; N b , N t , N s , N w are the total number of electric buses, taxis, private cars and commercial vehicles respectively, and the unit is vehicle ; Divide a day into T time periods; P ibj , P ibj , P ibj , P ibj are the charging power of different types of electric vehicles at time j, and the unit is kW.
(4)由于充电功率恒定,则充电负荷仅与充电时长相关,将电动汽车的充电时长作为起始充电时间的抽样约束条件,即可选取起始充电时间的抽样范围,然后抽取起始充电时间,从而求得各个抽样点对应时刻的充电负荷;充电时长的计算方法如下:(4) Since the charging power is constant, the charging load is only related to the charging time. Taking the charging time of the electric vehicle as the sampling constraint of the initial charging time, the sampling range of the initial charging time can be selected, and then the initial charging time can be extracted , so as to obtain the charging load at the corresponding time of each sampling point; the calculation method of charging time is as follows:
电动汽车有常规充电和快速充电两种模式,具体充电情况比较复杂。为方便建模,假设汽车优先选择常规充电,常规充电时基于日行驶里程计算充电时长,先利用蒙特卡洛方法抽取日行驶里程并计算所需充电时长;以满足充电所需要时长的约束条件为前提来选取起始充电时间的抽样范围,再抽取起始充电时间,进行负荷计算。而选择快速充电模式时基于荷电状态计算充电时长,同样基于蒙特卡洛方法抽取起始充电时间,计算充电时间段内的剩余充电时长和满足到下一次充电时的行驶需求的充电时长,取两者时间短的一个为实际充电时长。Electric vehicles have two modes of regular charging and fast charging, and the specific charging situation is more complicated. For the convenience of modeling, it is assumed that the car chooses conventional charging first. During conventional charging, the charging time is calculated based on the daily driving mileage. First, the Monte Carlo method is used to extract the daily driving mileage and calculate the required charging time; the constraints to meet the required charging time are The premise is to select the sampling range of the initial charging time, and then extract the initial charging time for load calculation. When the fast charging mode is selected, the charging time is calculated based on the state of charge, and the initial charging time is also extracted based on the Monte Carlo method, and the remaining charging time in the charging time period and the charging time to meet the driving demand for the next charging are calculated. The shorter of the two times is the actual charging time.
更进一步地,步骤3)用LSTM深度学习算法预测电动汽车充电负荷,包括以下步骤:Furthermore, step 3) predicts the electric vehicle charging load with the LSTM deep learning algorithm, including the following steps:
(1)在进行电动汽车充电负荷预测前需要对输入数据进行归一化处理,以保证网络训练和应用能够有较好的结果。采用min-max标准化,公式如下所示:(1) It is necessary to normalize the input data before electric vehicle charging load prediction to ensure better results in network training and application. Using min-max standardization, the formula is as follows:
式中:Xo表示经过标准化后的电动汽车负荷矩阵;X表示原电动汽车负荷矩阵,该矩阵只有一行,对应的是各个时刻的负荷值;xmin、xmax分别表示原电动汽车负荷矩阵中负荷的最小值和最大值。In the formula: X o represents the standardized electric vehicle load matrix; X represents the original electric vehicle load matrix, which has only one row, corresponding to the load value at each moment; x min and x max respectively represent the Load minimum and maximum.
(2)确定历史充电负荷信息的存放。即遗忘门的输出:(2) Determine the storage of historical charging load information. That is, the output of the forget gate:
ft=sigmoid(Wf·[ht-1,xt]+bf)f t =sigmoid(W f ·[h t-1 ,x t ]+b f )
式中:ft表示遗忘门输出;Wf表示从输入到遗忘门的权重系数;ht-1表示上一时刻的输出;xt表示网络当前时刻输入;bf表示遗忘门偏置。In the formula: f t represents the output of the forget gate; W f represents the weight coefficient from the input to the forget gate; h t-1 represents the output at the previous moment; x t represents the input of the network at the current moment; b f represents the bias of the forget gate.
(3)确定新充电负荷信息的存放。第一部分为通过sigmoid函数得到的输入门输出;另一部分为通过tanh建立的新的候选向量:(3) Determining the storage of new charging load information. The first part is the input gate output obtained by the sigmoid function; the other part is the new candidate vector established by tanh:
it=sigmoid(Wi·[ht-1,xt]+bi)i t = sigmoid(W i ·[h t-1 ,x t ]+b i )
式中:Wi表示输入门权重参数;bi表示输入门偏置;表示细胞充电负荷状态候选值,WC表示细胞充电负荷状态候选值计算权重参数;bC表示细胞充电负荷状态候选值计算偏置。所述的细胞表示这个时刻的单元,可根据之前蒙特卡洛抽样点间隔确定;In the formula: W i represents the weight parameter of the input gate; b i represents the bias of the input gate; Indicates the candidate value of the cell charging load state, W C represents the weight parameter for calculating the candidate value of the cell charging load state; b C represents the calculation bias for the candidate value of the cell charging load state. The cell represents the unit at this moment, which can be determined according to the previous Monte Carlo sampling point interval;
(4)更新细胞充电负荷状态:(4) Update the cell charging load status:
(5)确定当前时刻充电负荷输出。LSTM中,细胞充电负荷状态经过tanh函数处理和输出门过滤后,得到最终输出:(5) Determine the current charging load output. In LSTM, the cell charging load state is processed by the tanh function and filtered by the output gate to obtain the final output:
ot=sigmoid(Wo·[ht-1,xt]+bo)o t =sigmoid(W o ·[h t-1 ,x t ]+b o )
ht=ot*tanh(Ct)h t =o t *tanh(C t )
式中:ot表示输出门输出;Wo表示输入到输出门的权重参数;bo表示输出门偏置;ht表示当前时刻的输出。In the formula: o t represents the output of the output gate; W o represents the weight parameter input to the output gate; b o represents the bias of the output gate; h t represents the output at the current moment.
(6)将蒙特卡洛仿真获得的电动汽车充电负荷样本数据采用上述(1)-(5)的步骤进行处理,最终得到LSTM网络的电动汽车充电负荷预测结果。(6) The sample data of electric vehicle charging load obtained by Monte Carlo simulation is processed by the above steps (1)-(5), and finally the prediction result of electric vehicle charging load of LSTM network is obtained.
(7)对电动汽车充电负荷预测结束后,采用均方根误差预测精度评价指标:(7) After the electric vehicle charging load prediction is completed, the root mean square error prediction accuracy evaluation index is used:
本发明的的有益效果在于:The beneficial effects of the present invention are:
本发明建立了一种基于蒙特卡洛和深度学习的电动汽车充电负荷预测方法,根据该方法可得到不同类型电动汽车充电功率的计算模型和各时刻电动汽车的充电负荷。本发明基于蒙特卡洛仿真和深度学习算法对未来一天的电动汽车充电负荷进行预测,得出未来电动汽车充电负荷的整体趋势。仿真结果表明,本发明所提方法能够精确预测各时刻的充电负荷,有效地证明了本发明充电负荷预测方法的科学性、客观性。The present invention establishes a method for predicting the charging load of electric vehicles based on Monte Carlo and deep learning. According to the method, the calculation models of charging power of different types of electric vehicles and the charging load of electric vehicles at each moment can be obtained. The present invention predicts the charging load of the electric vehicle in the coming day based on the Monte Carlo simulation and the deep learning algorithm, and obtains the overall trend of the charging load of the electric vehicle in the future. The simulation results show that the method proposed by the invention can accurately predict the charging load at each moment, which effectively proves the scientificity and objectivity of the charging load prediction method of the invention.
附图说明Description of drawings
图1不同类型电动汽车充电方式及充电时间分布;Figure 1 Charging methods and charging time distribution of different types of electric vehicles;
图2蒙特卡洛抽样得到的10天电动汽车充电负荷数据;Figure 2 The 10-day electric vehicle charging load data obtained by Monte Carlo sampling;
图3LSTM结构图;Figure 3 LSTM structure diagram;
图4利用前9天数据学习后预测得到第10天数据;Figure 4 uses the data of the first 9 days to predict the data of the 10th day after learning;
图5隐藏层节点数为400LSTM负荷预测结果与实际值;Figure 5. The number of nodes in the hidden layer is the 400LSTM load prediction result and actual value;
图6本发明的一种基于蒙特卡洛和深度学习的充电负荷预测方法流程图。Fig. 6 is a flowchart of a charging load forecasting method based on Monte Carlo and deep learning of the present invention.
具体实施方式Detailed ways
下面结合附图和具体实施例对本发明作进一步的说明。The present invention will be further described below in conjunction with the accompanying drawings and specific embodiments.
本发明的一种基于蒙特卡洛和深度学习的电动汽车充电负荷预测方法,包括如下步骤:A kind of electric vehicle charging load prediction method based on Monte Carlo and deep learning of the present invention comprises the following steps:
S01.构建起始充电时刻模型,通过数据拟合处理得到电动汽车的起始充电时刻满足下式所示的正态分布:S01. Build the initial charging time model, and obtain the initial charging time of the electric vehicle through data fitting processing to meet the normal distribution shown in the following formula:
式中:t为初始充电时刻,即最后一次出行的结束时刻;μa和σa分别为起始充电时刻的期望和方差,不同类型电动汽车的期望和方差不同。In the formula: t is the initial charging time, that is, the end time of the last trip; μ a and σ a are the expectation and variance of the initial charging time, respectively, and the expectations and variances of different types of electric vehicles are different.
S02.构建日行驶里程模型,日行驶里程是电动汽车驾驶特性的重要指标,反映了汽车在一天内的耗电量,进而影响电动汽车的充电时间。同理,用传统燃油车的出行特性代替电动汽车的出行特性进行分析,电动汽车的日行驶里程服从对数正态分布,其概率密度函数为:S02. Construct a daily mileage model. The daily mileage is an important indicator of the driving characteristics of an electric vehicle, reflecting the power consumption of the vehicle in a day, which in turn affects the charging time of the electric vehicle. Similarly, using the travel characteristics of traditional fuel vehicles instead of electric vehicles for analysis, the daily mileage of electric vehicles follows a lognormal distribution, and its probability density function is:
式中:s为日行驶里程,单位为km;μb和σb分别为行驶里程s的对数lns的期望和方差,随着不同类型电动汽车行驶特性的不同而变化。In the formula: s is the daily mileage in km; μ b and σ b are the expectation and variance of the logarithm lns of the mileage s, respectively, which vary with the driving characteristics of different types of electric vehicles.
S03.构建起始荷电状态模型,在一个充电周期内,电动汽车充电所需的电力需求是随时间变化的。要确定电动汽车的充电负荷必须要获得电池充电开始时刻的荷电状态,即起始SOC。起始SOC是一个上次充电后的电动汽车行驶距离的随机函数,取值范围为0到100%。假设电动汽车的荷电状态随着行驶里程是线性下降的,则可通过车辆的行驶里程估算出起始充电时刻的荷电状态,用概率密度函数来表示:S03. Construct the initial state-of-charge model. In a charging cycle, the power demand required for charging an electric vehicle changes with time. To determine the charging load of electric vehicles, it is necessary to obtain the state of charge at the beginning of battery charging, that is, the initial SOC. The starting SOC is a random function of the distance traveled by the electric vehicle since the last charge, with values ranging from 0 to 100%. Assuming that the state of charge of an electric vehicle decreases linearly with the mileage, the state of charge at the initial charging time can be estimated from the mileage of the vehicle, expressed by a probability density function:
SOC=(SOC0-s/smax)×100%S OC =(S OC0 -s/s max )×100%
式中:SOC表示电池充电的起始SOC;SOC0表示最近一次充电后电池的荷电状态值,s为日行驶里程。由于电动汽车的充电时间和地点更为分散,SOC0往往不为1;smax表示电池充满后可行驶的最大里程数,单位km。In the formula: S OC represents the initial SOC of battery charging; S OC0 represents the state of charge value of the battery after the latest charge, and s is the daily mileage. Since the charging time and location of electric vehicles are more scattered, S OC0 is often not 1; s max indicates the maximum mileage that can be driven after the battery is fully charged, in km.
S04.构建充电时长模型,目前,电动汽车大多使用锂电池,锂电池的充电过程是恒压-恒流两阶段充电过程。假设整个充电过程是恒定功率的。充电时长受较多因素的制约,基于荷电状态计算充电时长,则有如下公式:S04. Build a charging time model. At present, most electric vehicles use lithium batteries, and the charging process of lithium batteries is a two-stage charging process of constant voltage and constant current. Assume that the entire charging process is constant power. The charging time is restricted by many factors. To calculate the charging time based on the state of charge, there is the following formula:
式中:Te为充电时长,单位为h;U为电池容量,单位为kW·h;P为充电功率,单位为kW;η为充电效率。In the formula: T e is the charging time, the unit is h; U is the battery capacity, the unit is kW h; P is the charging power, the unit is kW; η is the charging efficiency.
S05.基于日行驶里程的方式计算充电时长,则如下式所示:S05. Calculate the charging time based on the daily mileage, as shown in the following formula:
式中:W100为汽车每行驶100km的耗电量,单位为(kW·h)/百公里。In the formula: W 100 is the power consumption of the car every 100km, and the unit is (kW·h)/100km.
S06.划分不同类型电动汽车出行特征:将电动汽车按照不同的用途分为公交车、出租车、私家车和公务车4类,具体分析不同类型电动汽车的出行特征,以得到负荷预测模型的参数。4类电动汽车充电方式及充电时间分布如附图1所示。S06. Divide the travel characteristics of different types of electric vehicles: divide electric vehicles into four categories according to different purposes: buses, taxis, private cars, and official vehicles, and analyze the travel characteristics of different types of electric vehicles in detail to obtain the parameters of the load forecasting model . The charging methods and charging time distribution of the four types of electric vehicles are shown in Figure 1.
S07.预测各种类型电动汽车保有量:以《中国汽车产业发展报告(2012)》中指出的2020年的数据为基数,根据不同类型电动汽车的比例,对各种类型的汽车数量进行预测。不同类型电动汽车的保有量预测情况如表1所示。S07. Predict the number of various types of electric vehicles: Based on the data of 2020 pointed out in the "China Automobile Industry Development Report (2012)", the number of various types of vehicles is predicted according to the proportion of different types of electric vehicles. The forecast of ownership of different types of electric vehicles is shown in Table 1.
表1 不同类型电动汽车的保有量预测情况Table 1 Prediction of ownership of different types of electric vehicles
S08.基于蒙特卡洛仿真对电动汽车充电负荷进行计算:S08. Calculate the charging load of electric vehicles based on Monte Carlo simulation:
做合理假设:Make reasonable assumptions:
a.各类型电动汽车的起始充电时刻、日行驶里程和充电功率为相互独立的随机变量;a. The initial charging time, daily mileage and charging power of various types of electric vehicles are independent random variables;
b.各类型电动汽车的充电功率视为恒功率模型,没有恒压和恒流充电这2个阶段;b. The charging power of various types of electric vehicles is regarded as a constant power model, and there are no two stages of constant voltage and constant current charging;
c.所有车辆每次都充满电量;c. All vehicles are fully charged every time;
d.所有车辆的最后一次出行结束时刻即为该车辆的起始充电时刻。d. The end time of the last trip of all vehicles is the starting charging time of the vehicle.
选择蒙特卡洛方法模拟法对不同类型的电动汽车进行采样,通过抽取起始充电时间和日行驶里程,得到单辆电动汽车的充电时长和充电功率。然后累加所有电动汽车充电功率,得到各时刻的充电负荷。以时间为横坐标,各时刻的充电负荷为纵坐标,可得到充电负荷曲线。则j时刻的充电负荷为:The Monte Carlo simulation method is selected to sample different types of electric vehicles, and the charging time and charging power of a single electric vehicle are obtained by extracting the initial charging time and daily mileage. Then add up the charging power of all electric vehicles to get the charging load at each moment. Taking time as the abscissa and the charging load at each moment as the ordinate, the charging load curve can be obtained. Then the charging load at time j is:
式中:其中,Ptotal,j为在j时刻的充电负荷,单位为kw;Nb,Nt,Ns,Nw分别为电动公交车、出租车、私家车和商务车的总数,单位为辆;将一天划分为T个时间段;Pibj,Pibj,Pibj,Pibj为不同类型的电动汽车在j时刻的充电功率,单位为kW。In the formula: among them, P total, j is the charging load at time j, the unit is kw; N b , N t , N s , N w are the total number of electric buses, taxis, private cars and commercial vehicles, respectively, and the unit is is a vehicle; divide a day into T time periods; P ibj , P ibj , P ibj , P ibj are the charging power of different types of electric vehicles at time j, and the unit is kW.
S09.由于充电功率恒定,则充电负荷仅与充电时长相关,将电动汽车的充电时长作为起始充电时间的抽样约束条件,即可选取起始充电时间的抽样范围,然后抽取起始充电时间,从而求得各个抽样点对应时刻的充电负荷;充电时长的计算方法如下:S09. Since the charging power is constant, the charging load is only related to the charging time. Taking the charging time of the electric vehicle as the sampling constraint of the initial charging time, the sampling range of the initial charging time can be selected, and then the initial charging time can be extracted. In this way, the charging load at the corresponding time of each sampling point is obtained; the calculation method of the charging time is as follows:
假设汽车优先选择常规充电,常规充电时基于日行驶里程计算充电时长,先利用蒙特卡洛方法抽取日行驶里程并计算所需充电时长;以满足充电所需要时长的约束条件为前提来选取起始充电时间的抽样范围,再抽取起始充电时间,进行负荷计算。而选择快速充电模式时基于荷电状态计算充电时长,同样基于蒙特卡洛方法抽取起始充电时间,计算充电时间段内的剩余充电时长和满足到下一次充电时的行驶需求的充电时长,取两者时间短的一个为实际充电时长。蒙特卡洛仿真所得10天的电动汽车充电负荷数据如附图2所示。Assume that the car chooses conventional charging first. During conventional charging, the charging time is calculated based on the daily driving mileage. First, the Monte Carlo method is used to extract the daily driving mileage and calculate the required charging time; the starting point is selected on the premise of satisfying the constraints of the required charging time. The sampling range of the charging time, and then extract the initial charging time for load calculation. When the fast charging mode is selected, the charging time is calculated based on the state of charge, and the initial charging time is also extracted based on the Monte Carlo method, and the remaining charging time in the charging time period and the charging time to meet the driving demand for the next charging are calculated. The shorter of the two times is the actual charging time. The 10-day electric vehicle charging load data obtained by Monte Carlo simulation is shown in Figure 2.
S10.在进行电动汽车充电负荷预测前需要对输入数据进行归一化处理,以保证网络训练和应用能够有较好的结果。采用min-max标准化,公式如下所示:S10. It is necessary to normalize the input data before electric vehicle charging load prediction, so as to ensure better results in network training and application. Using min-max standardization, the formula is as follows:
式中:Xo表示经过标准化后的电动汽车负荷矩阵;X表示原电动汽车负荷矩阵,该矩阵只有一行,对应的是各个时刻的负荷值;xmin、xmax分别表示原电动汽车负荷矩阵中负荷的最小值和最大值。In the formula: X o represents the standardized electric vehicle load matrix; X represents the original electric vehicle load matrix, which has only one row, corresponding to the load value at each moment; x min and x max respectively represent the Load minimum and maximum.
S11.LSTM结构图如附图3所示。确定历史充电负荷信息的存放。即遗忘门的输出:S11. The structure diagram of LSTM is shown in Figure 3. Determine the storage of historical charging load information. That is, the output of the forget gate:
ft=sigmoid(Wf·[ht-1,xt]+bf)f t =sigmoid(W f ·[h t-1 ,x t ]+b f )
式中:ft表示遗忘门输出;Wf表示从输入到遗忘门的权重系数;xt表示网络当前时刻输入;bf表示遗忘门偏置。In the formula: f t represents the output of the forget gate; W f represents the weight coefficient from the input to the forget gate; x t represents the input of the network at the current moment; b f represents the bias of the forget gate.
S12.确定新充电负荷信息的存放。第一个部分为通过sigmoid函数得到的输入门输出;另一部分为通过tanh建立的新的候选向量:S12. Determining storage of new charging load information. The first part is the input gate output obtained by the sigmoid function; the other part is the new candidate vector established by tanh:
it=sigmoid(Wi·[ht-1,xt]+bi)i t = sigmoid(W i ·[h t-1 ,x t ]+b i )
式中:Wi表示输入门权重参数;bi表示输入门偏置;表示细胞充电负荷状态候选值;WC表示细胞充电负荷状态候选值计算权重参数;bC表示细胞充电负荷状态候选值计算偏置。In the formula: W i represents the weight parameter of the input gate; b i represents the bias of the input gate; Indicates the candidate value of the cell charging load state; W C represents the calculation weight parameter of the cell charging load state candidate value; b C represents the calculation bias of the cell charging load state candidate value.
S13.更新细胞充电负荷状态:S13. Update the cell charging load state:
S14.确定当前时刻充电负荷输出。LSTM中,细胞充电负荷状态经过tanh函数处理和输出门过滤后,得到最终输出:S14. Determine the charging load output at the current moment. In LSTM, the cell charging load state is processed by the tanh function and filtered by the output gate to obtain the final output:
ot=sigmoid(Wo·[ht-1,xt]+bo)o t =sigmoid(W o ·[h t-1 ,x t ]+b o )
ht=ot*tanh(Ct)h t =o t *tanh(C t )
式中:ot表示输出门输出;Wo表示输入到输出门的权重参数;bo表示输出门偏置;ht表示当前时刻的输出。In the formula: o t represents the output of the output gate; W o represents the weight parameter input to the output gate; b o represents the bias of the output gate; h t represents the output at the current moment.
S15.将蒙特卡洛仿真获得的电动汽车充电负荷样本数据经过上述步骤后,最终得到LSTM网络的电动汽车充电负荷预测结果。S15. After the sample data of the electric vehicle charging load obtained by the Monte Carlo simulation undergoes the above steps, the prediction result of the electric vehicle charging load of the LSTM network is finally obtained.
S16.对电动汽车充电负荷预测结束后,采用均方根误差预测精度评价指标:S16. After the electric vehicle charging load prediction is completed, the root mean square error prediction accuracy evaluation index is used:
应用例:Application example:
某市有电动汽车共12.5万辆,不同类型的电动汽车比例如表1所示。其中,电动公交车一天两充,充电时段分别服从正态分布N(13,12)和N(23,12),快速充电和常规充电的功率分别为108kW、60kW。日行驶里程服从对数正态分布N(4.3,0.342)。公交车以“比亚迪K9”为例,电动公交车电池容量为324kW·h,续航里程为250km,行驶100km的耗电量为140(kW·h)/百公里。电动出租车一天两充,充电时段分别服从正态分布N(3.5,12)和N(13,12),快速充电的功率为60kW。日行驶里程服从对数正态分布N(5.2,0.342)。以“比亚迪E6”型电动汽车为电动出租车、私家车和商务车的分析对象。该型号的电动汽车电池容量为82kW·h,续航里程为400km,行驶100km的耗电量为20.5(kW·h)/百公里,其快速充电模式的充电功率为60kW,常规充电模式的充电功率为7kW。电动私家车在上午、下午和晚上的充电概率分别设为0.2,0.1,0.7具体充电时间分别服从正态分布N(9.5,12)、N(14,12)、N(19,22),日行驶里程服从对数正态分布N(3.1,O.862)。公务车的充电时段服从正态分布N(19,22),日行驶里程服从对数正态分布N(3.1,O.862)。There are a total of 125,000 electric vehicles in a city, and the proportions of different types of electric vehicles are shown in Table 1. Among them, electric buses are charged twice a day, and the charging periods follow the normal distribution N(13, 1 2 ) and N(23, 1 2 ), respectively, and the powers of fast charging and conventional charging are 108kW and 60kW, respectively. The daily driving mileage obeys lognormal distribution N(4.3, 0.34 2 ). Taking "BYD K9" as an example, the battery capacity of the electric bus is 324kW h, the cruising range is 250km, and the power consumption for driving 100km is 140(kWh)/100km. Electric taxis are charged twice a day, and the charging periods follow the normal distribution N(3.5, 1 2 ) and N(13, 1 2 ), respectively, and the fast charging power is 60kW. The daily driving mileage obeys the logarithmic normal distribution N(5.2, 0.34 2 ). Take the "BYD E6" electric vehicle as the analysis object of electric taxis, private cars and commercial vehicles. The battery capacity of this type of electric vehicle is 82kW h, the cruising range is 400km, the power consumption for driving 100km is 20.5(kW h)/100km, the charging power of the fast charging mode is 60kW, and the charging power of the conventional charging mode It is 7kW. The charging probabilities of electric private cars in the morning, afternoon and evening are set to 0.2 , 0.1 and 0.7 respectively . ), the daily mileage obeys lognormal distribution N(3.1, O.86 2 ). The charging period of official vehicles obeys the normal distribution N(19, 2 2 ), and the daily mileage obeys the logarithmic normal distribution N(3.1, O.86 2 ).
以15分钟为一个采样点,采用蒙特卡洛模拟一天抽取共96个点来模拟全部电动汽车的充电负荷,共采取9天的负荷数据作为样本数据。采用提出的LSTM深度学习算法预测得到第10天的负荷数据如附图4所示,与通过蒙特卡洛仿真得到的第10天采样数据比较如附图5所示。Taking 15 minutes as a sampling point, Monte Carlo simulation is used to extract a total of 96 points a day to simulate the charging load of all electric vehicles, and a total of 9 days of load data are taken as sample data. The load data on the 10th day predicted by the proposed LSTM deep learning algorithm is shown in Figure 4, and the comparison with the sampling data on the 10th day obtained through Monte Carlo simulation is shown in Figure 5.
可见,采用LSTM深度学习预测得到的电动汽车充电负荷与通过蒙特卡洛采样得到的充电负荷情况基本吻合,误差比较小。其中LSTM预测结果误差与BP神经网络预测结果误差如表2所示,可见,深度学习所得结果相比BP神经网络这样的浅层学习结果更佳。It can be seen that the electric vehicle charging load predicted by LSTM deep learning is basically consistent with the charging load obtained by Monte Carlo sampling, and the error is relatively small. The LSTM prediction result error and the BP neural network prediction result error are shown in Table 2. It can be seen that the result of deep learning is better than that of shallow learning such as BP neural network.
表2 LSTM与BP预测误差比较Table 2 Comparison of LSTM and BP prediction errors
4种类型电动汽车的总充电负荷的峰值处于14:00—15:00,约为3900MW。这是因为选择电动汽车在该时段充电的电动汽车比较多,且多采用快速充电模式进行充电。充电负荷曲线表明,充电负荷在一天内的波动甚是剧烈,对电网的稳定具有潜在的威胁。因此,该发明结果为电网方面考虑相应指导性的充电策略以引导规范用户的充电行为提供支持。The peak of the total charging load of the four types of electric vehicles is between 14:00 and 15:00, which is about 3900MW. This is because there are more electric vehicles that choose to charge electric vehicles during this period, and most of them use fast charging mode for charging. The charging load curve shows that the charging load fluctuates violently within a day, which poses a potential threat to the stability of the power grid. Therefore, the result of the invention provides support for the power grid to consider corresponding guiding charging strategies to guide and regulate the charging behavior of users.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910725799.3A CN110570014B (en) | 2019-08-07 | 2019-08-07 | Electric vehicle charging load prediction method based on Monte Carlo and deep learning |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910725799.3A CN110570014B (en) | 2019-08-07 | 2019-08-07 | Electric vehicle charging load prediction method based on Monte Carlo and deep learning |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110570014A true CN110570014A (en) | 2019-12-13 |
CN110570014B CN110570014B (en) | 2022-04-01 |
Family
ID=68774762
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910725799.3A Active CN110570014B (en) | 2019-08-07 | 2019-08-07 | Electric vehicle charging load prediction method based on Monte Carlo and deep learning |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110570014B (en) |
Cited By (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111291999A (en) * | 2020-02-21 | 2020-06-16 | 东南大学 | A Load Modeling Method for Dynamic Wireless Charging Sections Based on Hierarchical Agglomerative Clustering |
CN111626494A (en) * | 2020-05-22 | 2020-09-04 | 广东电网有限责任公司 | Power load prediction method, power load prediction device, computer equipment and storage medium |
CN111626514A (en) * | 2020-05-29 | 2020-09-04 | 深圳供电局有限公司 | A kind of electric vehicle charging load prediction method and device |
CN111626527A (en) * | 2020-06-10 | 2020-09-04 | 太原理工大学 | Intelligent power grid deep learning scheduling method considering fast/slow charging/discharging form of schedulable electric vehicle |
CN112001564A (en) * | 2020-09-04 | 2020-11-27 | 国网江西省电力有限公司电力科学研究院 | User leaving time prediction method based on long-term and short-term memory network |
CN112215415A (en) * | 2020-09-29 | 2021-01-12 | 长安大学 | Prediction method of vehicle charging load scene based on optimal quantile of probability model |
CN112241596A (en) * | 2020-10-26 | 2021-01-19 | 国网上海市电力公司 | A mid- and long-term production simulation method for urban power grids considering electric vehicles |
CN112329984A (en) * | 2020-10-08 | 2021-02-05 | 国网河南省电力公司封丘县供电公司 | Electric vehicle optimized scheduling method based on electric vehicle cluster system division |
CN112330025A (en) * | 2020-11-06 | 2021-02-05 | 国网冀北电力有限公司张家口供电公司 | Prediction method of space-time charging load for urban electric vehicle |
CN112613731A (en) * | 2020-12-21 | 2021-04-06 | 国网河南省电力公司经济技术研究院 | Energy efficiency evaluation system and method for electric vehicle charging station based on complex network theory |
CN112613682A (en) * | 2020-12-29 | 2021-04-06 | 国网江苏省电力有限公司宜兴市供电分公司 | Electric vehicle charging load prediction method |
CN112785050A (en) * | 2021-01-12 | 2021-05-11 | 国网浙江省电力有限公司湖州供电公司 | Ordered charging scheduling method based on electric vehicle charging load prediction |
CN113268815A (en) * | 2021-06-24 | 2021-08-17 | 南方电网科学研究院有限责任公司 | Power distribution network double-layer planning method, device, equipment and storage medium |
CN113298298A (en) * | 2021-05-10 | 2021-08-24 | 国核电力规划设计研究院有限公司 | Charging pile short-term load prediction method and system |
CN113420953A (en) * | 2021-05-24 | 2021-09-21 | 国网上海市电力公司电力科学研究院 | Flexible load adjustability analysis method based on capacity credibility |
CN113435663A (en) * | 2021-07-15 | 2021-09-24 | 国网冀北电力有限公司唐山供电公司 | CNN-LSTM combined load prediction method considering electric vehicle charging load influence |
CN113554205A (en) * | 2020-04-26 | 2021-10-26 | 国网上海市电力公司 | Monte Carlo algorithm-based electric vehicle charging demand load research method |
CN113595113A (en) * | 2021-07-23 | 2021-11-02 | 上海电机学院 | Power distribution network regional electric vehicle charging load prediction method and storage medium |
CN113629709A (en) * | 2021-09-10 | 2021-11-09 | 广东电网有限责任公司 | Electric vehicle charging and discharging model obtaining method, system, equipment and storage medium |
CN114048894A (en) * | 2021-10-26 | 2022-02-15 | 国网福建省电力有限公司经济技术研究院 | Load prediction system and method for electric vehicle charging time point based on spatiotemporal characteristics |
CN114239967A (en) * | 2021-12-20 | 2022-03-25 | 深圳供电局有限公司 | Electric vehicle load prediction method, system and storage medium |
CN114254817A (en) * | 2021-12-09 | 2022-03-29 | 深圳供电局有限公司 | Method and system for predicting total charging load of electric vehicle |
CN114498634A (en) * | 2022-02-17 | 2022-05-13 | 四川大学 | Electric vehicle charging load prediction method based on ammeter data |
CN114580789A (en) * | 2022-04-08 | 2022-06-03 | 东南大学 | Multi-type electric vehicle load prediction method and system based on charging behavior |
CN114742118A (en) * | 2020-12-23 | 2022-07-12 | 中国科学院广州能源研究所 | Electric vehicle cluster charging and discharging load combination prediction method |
CN114781733A (en) * | 2022-04-28 | 2022-07-22 | 国网上海市电力公司 | A charging demand prediction method based on the spatiotemporal distribution of electric vehicles |
CN114881387A (en) * | 2021-12-16 | 2022-08-09 | 西安交通大学 | An orderly scheduling method and system for a vehicle-network cooperative cluster |
CN114914899A (en) * | 2022-06-15 | 2022-08-16 | 国网智慧能源交通技术创新中心(苏州)有限公司 | Electric vehicle power demand prediction method and system, storage medium and charging station |
CN114971067A (en) * | 2022-06-17 | 2022-08-30 | 长沙理工大学 | Building-new energy automobile load prediction method considering personnel characteristics |
WO2022252426A1 (en) * | 2021-05-31 | 2022-12-08 | 国网江苏省电力有限公司电力科学研究院 | Method for determining controllability of electric vehicle cluster, scheduling method and system |
CN116203448A (en) * | 2023-02-07 | 2023-06-02 | 江南大学 | Prediction method of power battery remaining life based on Monte Carlo and deep learning |
CN116373666A (en) * | 2023-03-28 | 2023-07-04 | 无锡市政公用新能源科技有限公司 | Charging dispatching system for active power distribution of charging pile |
CN117493878A (en) * | 2023-10-09 | 2024-02-02 | 华中科技大学 | Power prediction data enhancement method and power prediction method of electric vehicle charging station |
CN118539420A (en) * | 2024-05-17 | 2024-08-23 | 南京工程学院 | EV charging load prediction method based on hierarchical quantum clustering and user portrayal |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015022469A (en) * | 2013-07-18 | 2015-02-02 | 三菱重工業株式会社 | Power management device of electric automobile and power demand prediction method of electric automobile |
CN107067110A (en) * | 2017-04-14 | 2017-08-18 | 天津大学 | Charging electric vehicle load spatio-temporal prediction method under car Road network pattern |
CN107392400A (en) * | 2017-09-04 | 2017-11-24 | 重庆大学 | Meter and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodology |
CN109034498A (en) * | 2018-08-31 | 2018-12-18 | 国网上海市电力公司 | Consider the electric car charging load forecasting method of user's charge frequency and charge power variation |
-
2019
- 2019-08-07 CN CN201910725799.3A patent/CN110570014B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015022469A (en) * | 2013-07-18 | 2015-02-02 | 三菱重工業株式会社 | Power management device of electric automobile and power demand prediction method of electric automobile |
CN107067110A (en) * | 2017-04-14 | 2017-08-18 | 天津大学 | Charging electric vehicle load spatio-temporal prediction method under car Road network pattern |
CN107392400A (en) * | 2017-09-04 | 2017-11-24 | 重庆大学 | Meter and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodology |
CN109034498A (en) * | 2018-08-31 | 2018-12-18 | 国网上海市电力公司 | Consider the electric car charging load forecasting method of user's charge frequency and charge power variation |
Non-Patent Citations (2)
Title |
---|
乔文娟: "基于蒙特卡洛模拟的电动汽车充电负荷预测", 《华电技术》 * |
廖斌杰: "电动出租车充电负荷的时空随机分布特性", 《电力建设》 * |
Cited By (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111291999A (en) * | 2020-02-21 | 2020-06-16 | 东南大学 | A Load Modeling Method for Dynamic Wireless Charging Sections Based on Hierarchical Agglomerative Clustering |
CN113554205A (en) * | 2020-04-26 | 2021-10-26 | 国网上海市电力公司 | Monte Carlo algorithm-based electric vehicle charging demand load research method |
CN111626494A (en) * | 2020-05-22 | 2020-09-04 | 广东电网有限责任公司 | Power load prediction method, power load prediction device, computer equipment and storage medium |
CN111626514A (en) * | 2020-05-29 | 2020-09-04 | 深圳供电局有限公司 | A kind of electric vehicle charging load prediction method and device |
CN111626527A (en) * | 2020-06-10 | 2020-09-04 | 太原理工大学 | Intelligent power grid deep learning scheduling method considering fast/slow charging/discharging form of schedulable electric vehicle |
CN111626527B (en) * | 2020-06-10 | 2023-02-03 | 太原理工大学 | A deep learning scheduling method for smart grid considering the fast/slow charging and discharging forms of dispatchable electric vehicles |
CN112001564A (en) * | 2020-09-04 | 2020-11-27 | 国网江西省电力有限公司电力科学研究院 | User leaving time prediction method based on long-term and short-term memory network |
CN112001564B (en) * | 2020-09-04 | 2022-05-10 | 国网江西省电力有限公司电力科学研究院 | User leaving time prediction method based on long-short term memory network |
CN112215415A (en) * | 2020-09-29 | 2021-01-12 | 长安大学 | Prediction method of vehicle charging load scene based on optimal quantile of probability model |
CN112329984A (en) * | 2020-10-08 | 2021-02-05 | 国网河南省电力公司封丘县供电公司 | Electric vehicle optimized scheduling method based on electric vehicle cluster system division |
CN112241596A (en) * | 2020-10-26 | 2021-01-19 | 国网上海市电力公司 | A mid- and long-term production simulation method for urban power grids considering electric vehicles |
CN112330025A (en) * | 2020-11-06 | 2021-02-05 | 国网冀北电力有限公司张家口供电公司 | Prediction method of space-time charging load for urban electric vehicle |
CN112613731B (en) * | 2020-12-21 | 2023-02-24 | 国网河南省电力公司经济技术研究院 | Energy Efficiency Evaluation Method for Electric Vehicle Charging Stations Based on Complex Network Theory |
CN112613731A (en) * | 2020-12-21 | 2021-04-06 | 国网河南省电力公司经济技术研究院 | Energy efficiency evaluation system and method for electric vehicle charging station based on complex network theory |
CN114742118B (en) * | 2020-12-23 | 2023-10-27 | 中国科学院广州能源研究所 | Electric automobile cluster charge-discharge load combination prediction method |
CN114742118A (en) * | 2020-12-23 | 2022-07-12 | 中国科学院广州能源研究所 | Electric vehicle cluster charging and discharging load combination prediction method |
CN112613682A (en) * | 2020-12-29 | 2021-04-06 | 国网江苏省电力有限公司宜兴市供电分公司 | Electric vehicle charging load prediction method |
CN112785050A (en) * | 2021-01-12 | 2021-05-11 | 国网浙江省电力有限公司湖州供电公司 | Ordered charging scheduling method based on electric vehicle charging load prediction |
CN113298298B (en) * | 2021-05-10 | 2023-12-29 | 国核电力规划设计研究院有限公司 | Short-term load prediction method and system for charging pile |
CN113298298A (en) * | 2021-05-10 | 2021-08-24 | 国核电力规划设计研究院有限公司 | Charging pile short-term load prediction method and system |
CN113420953A (en) * | 2021-05-24 | 2021-09-21 | 国网上海市电力公司电力科学研究院 | Flexible load adjustability analysis method based on capacity credibility |
CN113420953B (en) * | 2021-05-24 | 2022-08-02 | 国网上海市电力公司电力科学研究院 | Flexible load adjustability analysis method based on capacity credibility |
WO2022252426A1 (en) * | 2021-05-31 | 2022-12-08 | 国网江苏省电力有限公司电力科学研究院 | Method for determining controllability of electric vehicle cluster, scheduling method and system |
CN113268815A (en) * | 2021-06-24 | 2021-08-17 | 南方电网科学研究院有限责任公司 | Power distribution network double-layer planning method, device, equipment and storage medium |
CN113435663A (en) * | 2021-07-15 | 2021-09-24 | 国网冀北电力有限公司唐山供电公司 | CNN-LSTM combined load prediction method considering electric vehicle charging load influence |
CN113595113A (en) * | 2021-07-23 | 2021-11-02 | 上海电机学院 | Power distribution network regional electric vehicle charging load prediction method and storage medium |
CN113629709B (en) * | 2021-09-10 | 2024-05-17 | 广东电网有限责任公司 | Method, system, equipment and storage medium for acquiring charge and discharge model of electric vehicle |
CN113629709A (en) * | 2021-09-10 | 2021-11-09 | 广东电网有限责任公司 | Electric vehicle charging and discharging model obtaining method, system, equipment and storage medium |
CN114048894A (en) * | 2021-10-26 | 2022-02-15 | 国网福建省电力有限公司经济技术研究院 | Load prediction system and method for electric vehicle charging time point based on spatiotemporal characteristics |
CN114254817A (en) * | 2021-12-09 | 2022-03-29 | 深圳供电局有限公司 | Method and system for predicting total charging load of electric vehicle |
CN114881387A (en) * | 2021-12-16 | 2022-08-09 | 西安交通大学 | An orderly scheduling method and system for a vehicle-network cooperative cluster |
CN114239967A (en) * | 2021-12-20 | 2022-03-25 | 深圳供电局有限公司 | Electric vehicle load prediction method, system and storage medium |
CN114498634A (en) * | 2022-02-17 | 2022-05-13 | 四川大学 | Electric vehicle charging load prediction method based on ammeter data |
CN114498634B (en) * | 2022-02-17 | 2023-08-29 | 四川大学 | Electric automobile charging load prediction method based on ammeter data |
CN114580789B (en) * | 2022-04-08 | 2024-05-03 | 东南大学 | A method and system for predicting load of multiple types of electric vehicles based on charging behavior |
CN114580789A (en) * | 2022-04-08 | 2022-06-03 | 东南大学 | Multi-type electric vehicle load prediction method and system based on charging behavior |
CN114781733A (en) * | 2022-04-28 | 2022-07-22 | 国网上海市电力公司 | A charging demand prediction method based on the spatiotemporal distribution of electric vehicles |
CN114914899A (en) * | 2022-06-15 | 2022-08-16 | 国网智慧能源交通技术创新中心(苏州)有限公司 | Electric vehicle power demand prediction method and system, storage medium and charging station |
CN114914899B (en) * | 2022-06-15 | 2024-09-06 | 国网智慧能源交通技术创新中心(苏州)有限公司 | Electric demand prediction method and system for electric automobile, storage medium and charging station |
CN114971067A (en) * | 2022-06-17 | 2022-08-30 | 长沙理工大学 | Building-new energy automobile load prediction method considering personnel characteristics |
CN116203448B (en) * | 2023-02-07 | 2023-10-03 | 江南大学 | Power battery residual life prediction method based on Monte Carlo and deep learning |
CN116203448A (en) * | 2023-02-07 | 2023-06-02 | 江南大学 | Prediction method of power battery remaining life based on Monte Carlo and deep learning |
CN116373666A (en) * | 2023-03-28 | 2023-07-04 | 无锡市政公用新能源科技有限公司 | Charging dispatching system for active power distribution of charging pile |
CN117493878A (en) * | 2023-10-09 | 2024-02-02 | 华中科技大学 | Power prediction data enhancement method and power prediction method of electric vehicle charging station |
CN118539420A (en) * | 2024-05-17 | 2024-08-23 | 南京工程学院 | EV charging load prediction method based on hierarchical quantum clustering and user portrayal |
Also Published As
Publication number | Publication date |
---|---|
CN110570014B (en) | 2022-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110570014B (en) | Electric vehicle charging load prediction method based on Monte Carlo and deep learning | |
Chen et al. | An on-line predictive energy management strategy for plug-in hybrid electric vehicles to counter the uncertain prediction of the driving cycle | |
Tao et al. | Data-driven optimized layout of battery electric vehicle charging infrastructure | |
Xydas et al. | A data-driven approach for characterising the charging demand of electric vehicles: A UK case study | |
CN111400662B (en) | Space load prediction method considering charging requirements of electric automobile | |
CN107169273A (en) | The charging electric vehicle power forecasting method of meter and delay and V2G charge modes | |
CN115327422A (en) | Evaluation method of electric bus power battery health based on charging and discharging behavior | |
CN112215434A (en) | LSTM model generation method, charging duration prediction method and medium | |
CN109934403A (en) | Analysis and prediction method of electric vehicle charging load in residential area based on mathematical model | |
CN106295860A (en) | A kind of electric automobile scale charge requirement Forecasting Methodology based on Monte Carlo Analogue Method | |
CN111047120B (en) | Electric automobile charging load prediction method under road-electricity coupling network | |
CN110232219B (en) | A Data Mining-Based Approval Method for Dispatchable Capacity of Electric Vehicles | |
CN109840640B (en) | Method and system for site selection of electric vehicle charging pile | |
CN112766677A (en) | Modeling method and charging system for residential electric vehicle charging demand based on fuzzy comprehensive evaluation method | |
CN107609670A (en) | One kind is based on copula algorithm electric automobile charging station load forecasting methods | |
Liu et al. | Load characteristics forecasting of Hubei power grid up to year 2030 with the development of electric vehicles | |
Chen et al. | On the relationship between energy consumption and driving behavior of electric vehicles based on statistical features | |
CN110826813A (en) | Grid optimization method based on differential charging demands of household electric vehicle users | |
CN112765726B (en) | Life prediction method and device | |
CN103701176B (en) | A kind of electric automobile soon, the computational methods of electrically-charging equipment allocation ratio at a slow speed | |
CN106203719A (en) | A kind of electric automobile accesses the load forecasting method of electrical network | |
He et al. | Multi-time simulation of electric taxicabs' charging demand based on residents' travel characteristics | |
CN116029140A (en) | Electric automobile short-term charging load prediction model considering temperature influence | |
Cheng et al. | Electric vehicles charging load prediction based on Monte Carlo method | |
CN115907128A (en) | Highway annual scale electricity load prediction method and device |
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 |