CN107392400A - Meter and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodology - Google Patents
Meter and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodology Download PDFInfo
- Publication number
- CN107392400A CN107392400A CN201710787528.1A CN201710787528A CN107392400A CN 107392400 A CN107392400 A CN 107392400A CN 201710787528 A CN201710787528 A CN 201710787528A CN 107392400 A CN107392400 A CN 107392400A
- Authority
- CN
- China
- Prior art keywords
- mrow
- charging
- time
- car
- traffic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 64
- 230000002123 temporal effect Effects 0.000 title claims abstract description 41
- 238000005265 energy consumption Methods 0.000 claims abstract description 27
- 238000004088 simulation Methods 0.000 claims abstract description 22
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 11
- 238000012546 transfer Methods 0.000 claims abstract description 11
- 238000012360 testing method Methods 0.000 claims abstract description 8
- 238000013459 approach Methods 0.000 claims abstract description 5
- 230000009471 action Effects 0.000 claims description 44
- 238000004378 air conditioning Methods 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 13
- 150000001875 compounds Chemical class 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims description 6
- 239000004744 fabric Substances 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 230000005611 electricity Effects 0.000 claims description 4
- 238000011068 loading method Methods 0.000 claims description 4
- 230000000977 initiatory effect Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 7
- 230000008859 change Effects 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 4
- 230000007613 environmental effect Effects 0.000 description 3
- 230000032258 transport Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000005381 potential energy Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005295 random walk Methods 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000009452 underexpressoin Effects 0.000 description 1
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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
- G06Q10/06375—Prediction of business process outcome or impact based on a proposed change
-
- 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
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Traffic Control Systems (AREA)
Abstract
The present invention relates to a kind of meter and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodologies, belong to intelligent grid field.The method comprising the steps of:S1:First all kinds of electric automobiles are classified, analyze its space transfer and charge characteristic;S2:Make abstract simulation to urban traffic network with Graph-theoretical Approach using Trip chain is theoretical;S3:The Time-spatial diversion model of all kinds of electric automobiles is set up using Monte Carlo simulation and the random shortest path method based on Markovian decision;S4:Temperature and traffic energy consumption model are established according to electric automobile actual test data.The present invention has versatility, analyzed suitable for regional distribution, spatial and temporal distributions situation in urban area can conveniently and effectively be calculated, and the change for the load that charged when traffic changes with temperature regime can be predicted, its prediction result provides foundation for researchs such as charging station planning, load schedulings.
Description
Technical field
The invention belongs to intelligent grid field, is related to meter and real-time traffic and the EV of temperature charging load spatial and temporal distributions are predicted
Method.
Background technology
As world energy sources structural adjustment and auto industry technology develop, electric automobile (EV) is the main of new-energy automobile
Developing direction.Extensive electric automobile access, very important influence, including load will be brought to Power System Planning and operation
Growth, operation of power networks optimal control difficulty increase, influence the quality of power supply, new requirement etc. proposed to distribution network planning.Solution
Certainly the basis of problem above is that charging load spatial and temporal distributions are effectively predicted, therefore accurately and effectively load is established in active demand
Forecast model.
The load prediction of foreign scholar's early start charging electric vehicle is studied.Early stage research is in charging set in advance
Between with charging place under, directly using trip survey data complete electric automobile load modeling.The load that existing literature is derived
Demand shows that electric automobile load can increase the peak-valley difference of power network.Because trip survey data are based on conventional truck, mostly
Assume that charging interval and place are fixed value in number existing literature, be generally set to fill in house after last time trip terminates
Electricity.Some researchs afterwards consider the randomness of charging time started, but still set duration of charge as fixed value.It is existing
Region division is industrial area, shopping centre residential quarter etc. by document, and charging load is calculated according to the different charge characteristic of regional
Spatial and temporal distributions, but the space distribution of actually electric automobile is predetermined, does not simulate electric automobile because of random shifting
Dynamic and caused load spatial and temporal distributions.There is scholar to be combined power network with the network of communication lines in recent years, consider the randomness of electric automobile,
A kind of space-time model of mobile electric car load is established by random Trip chain and Markovian decision process (MDP).Should
Method relatively accurately realizes the prediction of spatial and temporal distributions, but mainly for private car, does not consider traffic, temperature etc. to electronic vapour
The environmental factor that car energy consumption has a significant impact.
It is domestic that the random device using probability statistics is concentrated mainly on to load early stage to the research of charging electric vehicle load
Annual distribution is modeled, and existing literature analyzes the charging characteristics of taxi, bus, private car etc., using Monte Carlo
Method predicts the charging total load of China's future electric automobile.Existing literature considers user's daily travel, charging starting
The stochastic variables such as time, establish the probabilistic model of charge requirement.Research afterwards starts to consider load asking in spatial distribution
Topic, existing literature study different zones charging electric vehicle characteristic by way of region division, form spatial distribution.It is existing
Document construction trip chain model, Temporal And Spatial Distribution Model is established to five big trip regions.But the studies above divides more bases to space
In function plot, it is impossible to reach required precision of the analysis space distribution to Electric Power Network Planning operation.Existing literature proposes one kind
Charging load spatial and temporal distributions Forecasting Methodology based on traffic trip matrix and cloud model, this method embody the random of user's decision-making
Property and ambiguity, preferably realize the prediction of spatial and temporal distributions, but be not involved with the different trip of all kinds of electric automobiles with filling
Electrically, it also have ignored influence of the environmental factor to load prediction.
The content of the invention
In view of this, it is an object of the invention to provide space division when a kind of meter and real-time traffic and the EV of temperature charging loads
Cloth Forecasting Methodology.All kinds of electric automobiles are classified first, analyze its space transfer and charge characteristic, it is theoretical using Trip chain
Make abstract simulation to urban traffic network with Graph-theoretical Approach.Secondly, using Monte Carlo simulation and based on Markovian decision
Random shortest path method sets up the Time-spatial diversion model of all kinds of electric automobiles, and is established according to electric automobile actual test data
Temperature and traffic energy consumption model.Finally, by taking Jincheng City as an example, the validity of this method is demonstrated.
To reach above-mentioned purpose, the present invention provides following technical scheme:
Meter and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodology, comprise the following steps:
S1:First all kinds of electric automobiles are classified, analyze its space transfer and charge characteristic;
S2:Make abstract simulation to urban traffic network with Graph-theoretical Approach using Trip chain is theoretical;
S3:All kinds of electronic vapour are set up using Monte Carlo simulation and the random shortest path method based on Markovian decision
The Time-spatial diversion model of car;
S4:Traffic and temperature energy consumption model are established according to electric automobile actual test data.
Further, the step S1 is specially:
Charging modes are fixed parking stall trickle charge, and trip mode is come out early and return late for working day, stopped over, and day off goes out
The EV types of amusement are private car A classes;
Charging modes stop over, day off steps out to be come out early and return late by way of ground fast charge, trip mode for working day
EV types be private car B classes;
Charging modes are fixed parking stall trickle charge, and trip mode is that working day traffic path is fixed, what day off did not went on a journey
EV types are officer's car A classes;
Charging modes are by way of ground fast charge, and trip mode is that working day traffic path is fixed, the EV classes that day off does not go on a journey
Type is officer's car B classes;
Charging modes are by way of bus stop quick charge, and trip mode is that traffic path is fixed with the travel time, annual
The EV types of operation are bus;
Charging modes are to be fixed by way of ground fast charge station quick charge, trip mode to set out with the off-running time, and route is random,
The whole year EV types of operation are taxi.
Further, the step S2 is specially:
Trip chain is divided into simple chain and compound catenary;Simple chain only has one for the trip purpose of every Trip chain;It is complicated
Chain is that Trip chain contains multiple trip purposes;It is divided into residence H, place of working W, purchase by resident trip survey data to trip purpose
Thing dining SE, social entertainment SR, other O, while be also charging electric vehicle place;
Wherein, the electrically-charging equipment that residence H is set is fast charge+trickle charge;
The electrically-charging equipment that place of working W is set is fast charge+trickle charge;
The electrically-charging equipment that shopping dining SE is set is fast charge;
The electrically-charging equipment that social entertainment SR is set is fast charge;
The electrically-charging equipment that other O are set is fast charge;
Urban traffic network is abstracted as by digraph using the method for graph theory, each highway in city is arranged to side,
Intersection is arranged to node, road length is arranged to the weight on side, and transportation network is abstracted as to the set of side and node.
Further, the step S3 is specially:
For each EV in region, it is known that after its starting point, using the markov with random chance factor gamma
Decision process obtains the optimal path with randomness in its space;
Wherein, the Markovian decision process is the state observed according to each moment, from available action collection
Middle to be made decisions from an action, the state of system next step is random, and its state transition probability has Ma Er can
Husband's property;Markovian decision process is:
M=(T, S, A (i), p, R)
Wherein, T is decision-making moment collection, represents to choose the time point set of action;S is state set, is represented in each decision-making
Carve, system state space collection;A (i) collects for action, represents, at any one decision-making moment, the state that policymaker observes, that is, to exist
The available action collection that state i chooses;P (j | i, it is a) state transition probability, represents under state i, a that takes action is transferred to shape
State j probability;R (i, it is a) Reward Program, represents under state i, after a that takes action, the return of policymaker's acquisition;
That is Markovian decision dynamic process is as follows:Policymaker state S (0) where T (0) moment, is selected from A (0)
One action a (0) performs, and has arrived next S (1) state by P (S (1) | S (0), a (0)) probability random transferring after execution, then
An action a (1) is performed again, and successively backward, to the last moment T (N), reaches end-state S (N);
It is the rule that a decision system selects action on each different conditions to define decision rule, is designated as f (S (0)),
f(S(1)),f(S(2))…f(N);Definition strategy is the sequence of decision rule, is designated as π=(f (S (0)), f (S (1)), f (S
(2)) ... f ((S (N))), the collection of All Policies is collectively referred to as policy class, is designated as Π;From original state S, (0) simultaneously used tactful π
∈ Π so thatMaximum or minimum value is obtained, then tactful π is optimal policy;Define value function VπFor using tactful π when
State S (i) expected returns:
Wherein E is it is expected, StFor t status, t=0 is corresponded to original state S (0), is expressed as with recursive form:
To each tactful π, its corresponding value function VπIt is a series of unique public solution of linear equations, optimal policy passes through
Recursive method obtains expected return value based on Dynamic Programming backward;
All node sets in region are considered as S, each node time instance passed through from origin-to-destination traveling is considered as decision-making
Moment t;Determine that next paths of traveling are considered as action a in each node, the path length that will be travelled is considered as return R;
If EV under steam with part probability transfer mode travel, introduce value be 0~0.5 obedience it is equally distributed with
Machine probability factor γ, represent, when policymaker is often performed with this action, there is probability γ to be transferred to another paths;By Path selection
In shortest path be combined with randomness, do not influenceed by the original state and initial action set out, for next decision-making
Remaining decision rule composition strategy be optimal policy;
(1) for private car:
Working day, private car active strokes share 3 kinds:Active strokes 1:H-W;Active strokes 2:H-W-SR/SE/O;Activity
Stroke 3:H-W+H-SR/SE/O;It is each to obey certain proportion;Working day each car Trip chain C obeys probability distribution:P (C)=pi(i
=1,2,3)
X obeys Normal probability distribution between at the beginning of every section of stroke, and its probability density function is:
Wherein σeFor Variance of Normal Distribution, μeFor normal distribution average;
All electronic private cars are considered as totality, each electric automobile is considered as a two-dimensional random variable, and it, which is combined, divides
Cloth function is:P (C ∩ x)=P (C) P (C | x), when extracting Trip chain type by the method for Monte Carlo simulation and set out
Between;
The random most short used time path computing of meter and real-time traffic:If each link length is Li, the traveling of each road of t
Speed is Vi,t, each side right value is set as by T the time required to each sidei, then havePass through band random chance factor M DP
Algorithm obtains the relevant parameter such as random most short used time path, travel speed of corresponding automobile different sections of highway, obtain private car when
Space division cloth;
(2) for officer's car
Probability distribution between at the beginning of every section of stroke is identical, equal Normal Distribution, using the illiteracy similar with private car
Special Carlow analogy method obtains the spatial and temporal distributions of officer's car;
(3) for taxi
Each taxi is 0:When 00 from company location, destination is randomly selected, using meter and real-time traffic
Obtain most short used time path in most short used time road;It is random between at work to run time return taxi company of being handed over to the next shift until two,
It is corresponding with electric taxi travel speed, obtain the spatial and temporal distributions of each electric taxi;
(4) for bus
Time of departure obedience is uniformly distributed, and for each bus, hair is extracted by the way of Monte Carlo simulation
The car time makes the time, corresponding different sections of highway travel speed, to obtain the spatial and temporal distributions of bus with the parking of each website.
Further, the step S4 is specially:
S401:Traffic energy consumption model is established according to electric automobile actual test data:
The measured data modeling of the energy consumed based on the identical mileage of different category of roads downward drivings, obtains corresponding to often
Its different periods and the specific energy consumption under the real-time jam situation of different road conditions, through street unit mileage power consumption are Eksl=0.247
+1.520/V-0.004*V+2.992*10(-5)*V;
Trunk roads unit mileage power consumption is Ezgl=-0.179+0.004*V+5.492/V;
Secondary distributor road unit mileage power consumption is Ecgl=0.21-0.001*V+1.531/V;
Branch road unit mileage power consumption is Ezl=0.208-0.002*V+1.553/V;
Wherein Eksl、Ezgl、Ecgl、EzlThrough street, trunk roads, secondary distributor road, the electric automobile unit mileage of branch road are represented respectively
Power consumption, V represent electric automobile road-section average travel speed;
For each electric automobile in region, certain position is obtained according to its Temporal And Spatial Distribution Model at moment, it is corresponding real
When traffic congestion after obtain moment EV units mileage power consumption;
S402:Temperature energy consumption model is established according to electric automobile actual test data:
The calculation formula of least square air-conditioning unlatching rate is:
Kpect=7 × 10-5·TEMP4-0.0887·TEMP3+3.673·TEMP2-56.302·TEMP+310.7681
Wherein TEMP represents temperature, KpectAir-conditioning unlatching rate is represented, draws out its functional image;
Temperature energy consumption proportionality coefficient is defined, i.e., unit distance power consumption corresponding to different temperatures consumes with benchmark after air-conditioning is opened
The ratio between electricity, the fitting formula of least square temperature proportional coefficient are:
Ktemp=10-3·(TEMP-10)2+1.2
Wherein KtempFor temperature proportional coefficient, TEMP is temperature, draws out its functional image;
According to temperatures above energy consumption model, the real time temperature at corresponding estimation range moment, electric automobile in region is obtained
To the influence situation of electric automobile unit mileage power consumption after the ratio that air-conditioning is opened, and air-conditioning unlatching.
Further, the charging Load Calculation Method of all kinds of electric automobiles:
The charging matrix of loadings PE of m rows n row is defined, each time interval is represented per a line, each row represent each node,
Matrix element represents the charging load at certain place corresponding moment, i.e. PEi,jRepresent the charge power at i-node j moment;Time interval
Quantity m, i.e. row in matrix quantity, depending on temporal resolution and calculate requirement;
(1) A classes private car and officer's car:According to the traveling spatial and temporal distributions of each electric automobile, charging modes, obtain its
The charging load of place parking site;For A classes private car, officer's car, because its place of charging is fixed, each electronic vapour
Charging load corresponds to a certain row of PE matrixes caused by car;
(2) B classes private car, officer's car:Consider that environment temperature calculates unit mileage power consumption with traffic, correspond to this
Section path length obtains this section of distance power consumption ER1;If dump energy EC is not enough to be unsatisfactory for reaching next destination, i.e., not
Meet EC>ER1, then electric automobile current location carry out a quick charge, charging duration tcr taken out using Monte Carlo mode
Take, and be limited within 30min;If the charging time started is TS, the charging end time is TE, then has:TE=TS+tcr;
It is by needing to judge that can dump energy expire during next place with A class private car calculation process difference
Sufficient distance needs;For B classes private car, officer's car, corresponded to and PE matrixes by way of the load that charged caused by charging place at it
Each row;
(3) taxi:Taxi is A class private cars for a long time, without fixed destination traveling, and taxi is by way of charging
The corresponding each row with PE matrixes of charging load caused by place;
(4) bus:Bus is the A class private cars that route is fixed, battery capacity and unit mileage power consumption are bigger, public
Car is handed in the corresponding each row with PE matrixes of load that charged caused by each website.
The beneficial effects of the present invention are:The present invention has versatility, is analyzed suitable for regional distribution, can facilitate, have
Spatial and temporal distributions situation in urban area is calculated to effect, and the load that charged when traffic changes with temperature regime can be predicted
Change, its prediction result provides foundation for researchs such as charging station planning, load schedulings.
Brief description of the drawings
In order that the purpose of the present invention, technical scheme and beneficial effect are clearer, the present invention provides drawings described below and carried out
Explanation:
Fig. 1 is Trip chain type;
Fig. 2 is Markovian decision schematic diagram;
Fig. 3 is electric automobile Time-spatial diversion schematic diagram;
Fig. 4 is taxi space transfer figure;
Fig. 5 is air-conditioning unlatching rate curve;
Fig. 6 is temperature proportional coefficient;
Fig. 7 is charging matrix of loadings;
Fig. 8 is simulation contact surface;
Fig. 9 is Jincheng City transportation network figure;
Figure 10 is each node charging load of typical case's week lower Jincheng City A, B class car;
Figure 11 is each node charging load of typical case's week lower Jincheng City bus;
Figure 12 is each node charging load of typical case's week lower Jincheng City taxi;
Figure 13 is typical case's week lower Jincheng City charging electric vehicle total load;
Figure 14 is the typical case's week of node 10, high temperature week charging load curve;
Figure 15 is the typical case's week of node 15, congestion week charging load curve;
Figure 16 is two methods load prediction comparison curves.
Embodiment
Below in conjunction with accompanying drawing, the preferred embodiments of the present invention are described in detail.
1. electric automobile is classified and traveling network
1.1 electric automobile types
Electric automobile type is more in city, and the charging modes that same type electric automobile uses also are not quite similar,
For the spatial and temporal distributions of analysis charging load, the present invention from automobile charging modes and trip feature, by electric automobile be divided into
Lower several types, as shown in table 1:
The electric automobile type of table 1
1.2 Trip chain
Trip chain can be divided into simple chain and compound catenary.Simple chain is that the trip purpose of every Trip chain only has one, more
It is fixed;Compound catenary then refers to Trip chain and contains multiple trip purposes, and flexibility is larger with randomness.It can go out to trip purpose by resident
Row survey data is divided into the H that goes home (Home), work W (Work), shopping dining SE (Shopping and Eating), social joy
Five classes of happy SR (Social and Recreational) and other affairs O (Other family/personal errands)
Type, it is charging electric vehicle place to be also simultaneously these trip purposes.The common Trip chain such as Fig. 1 institutes being related in the present invention
Show, wherein (a) (b) (c) (d) is simple chain, (e) is compound catenary.
With charging modes it is trip purpose to be mutually related, is printed and distributed according to National Development and Reform Committee《Charging electric vehicle base
Infrastructure development guide (2015-2020)》, to trip purpose charging equipment do following setting:
The electrically-charging equipment of table 2 is set
The 1.3 random walk simulations based on Markov decision theory
Urban traffic network is abstracted as by digraph using the method for graph theory.Each highway in city is arranged to side,
Intersection is arranged to node, road length is arranged to the weight on side, and with this, transportation network can be abstracted as side and node
Set.
Markovian decision process refers to that policymaker periodically or continuously observes the stochastic and dynamic with Markov property
System, sequentially make decisions.The state observed according to each moment, from an action from available action collection
Make decisions, the state of system next step is random, and its state transition probability has Markov property.One Ma Er can
Husband's decision process is made up of a five-tuple
M=(T, S, A (i), p, R):
T:Decision-making moment collection, represent to choose the time point set of action.
S:State set, represent in each decision-making moment, system state space collection.
A(i):Action collection, represents at any one decision-making moment, the state that policymaker observes, can be chosen in state i
Available action collection.
p(j|i,a):State transition probability.Represent under state i, a that takes action is transferred to state j probability.
R(i,a):Reward Program.Under expression state i, after a that takes action, the return of policymaker's acquisition.
Markovian decision dynamic process is as follows:Policymaker state S (0) where T (0) moment, the selection one from A (0)
Individual action a (0) is performed, and next S (0) state has been arrived by P (S (1) | S (0), a (0)) probability random transferring after execution.Then again
An action a (1) is performed, successively backward, to the last moment T (N), reaches end-state S (N), as shown in Figure 2.
It is the rule that a decision system selects action on each different conditions to define decision rule, is designated as f (S (0)),
f(S(1)),f(S(2))…f(N).Definition strategy is the sequence of decision rule, is designated as π=(f (S (0)), f (S (1)), f (S
(2)) ... f ((S (N))), the collection of All Policies is collectively referred to as policy class, is designated as Π.If from original state S, (0) simultaneously used
Tactful π ∈ Π so thatMaximum or minimum value is obtained, then tactful π is referred to as optimal policy.Define value function VπFor using strategy
In state S (i) expected returns during π:
Wherein St is t status, and t=0 is corresponded to original state S (0), is expressed as with recursive form:
To each tactful π, its corresponding value function VπIt is a series of unique public solution of linear equations.Above formula is referred to as
Bellman formula, optimal policy can recursive method obtains backward by the expected return value based on Dynamic Programming.
All node sets in region are considered as S, each node time instance passed through from origin-to-destination traveling is considered as decision-making
Moment t;Determine that next paths of traveling are considered as action a in each node, the path length that will be travelled is considered as return R.
By taking 4 node regions as an example:State set S={ 1,2,3,4 }, action collection A={ a, b, c, d };R (a)=65m, R (b)
=95m, R (c)=90m, R (d)=97m, π 1={ a, b }, π 2={ c, d }, if 1 is starting point, 2 be terminal, then:
R(π1)=r (a)+r (b)=160m
R(π2)=r (c)+r (d)=187m (3)
Definition action state-transition matrix P, represent to take certain take action after the probability that is shifted between each state of policymaker.If
Take action a, b, c, and d is full probability transfer, i.e., is transferred to NextState with probability 1, by taking a that takes action as an example, then P matrixes are expressed as:
When policymaker is located at node 1, execution action a will be transferred to node 2 with probability 1, when being in other nodes then with
The holding position of probability 1 is constant.Now optimal policy is that shortest path path is π 1={ a, b }.
Due to the actual generally non-shortest path of electric automobile driver trip route, it is to the current region residing for oneself
State can not possibly have a complete understanding.Therefore, the method for probability factor is introduced to handle the uncertainty of this category information.
Assuming that EV is travelled with part probability transfer mode under steam, obeyed between introducing 0~0.5 equally distributed random
Probability factor γ, when policymaker is often performed with this action, there is probability γ to be transferred to another paths.To action a in upper example
Speech, its probability transfer matrix of taking action are:
There is γ probability to be transferred to node 3 after performing action a.It is corresponding random in the case of transferable path is more than 2
Probability factor γ1~γk, and meetRandom chance factor gamma be introduced into by the shortest path in Path selection with
Randomness is combined, because the recursive algorithm backward of MDP optimizing decision processes obeys principle of optimality, i.e., no matter at the beginning of which
Which initial action is beginning state set out and employ, and the strategy of remaining decision rule composition is optimal for next decision-making
Strategy.Therefore for each EV in region, it is known that optimal using the MDP with random chance factor gamma after its starting point
Process approach can obtain the optimal path with randomness in its space.
2. electric automobile Time-spatial diversion model
2.1 private car
2.1.1 Trip chain and departure time
On weekdays, the drive active strokes of working clan share 3 kinds:Active strokes 1 are to be gone from before house on job site
Class, After Hours goes home, does not go to other places, represent to be H-W with Trip chain;Active strokes 2 are gone to work from before house, are come off duty
SR/SE/O places are gone to before in way home once, then are gone home, Trip chain is expressed as H-W-SR/SE/O;Active strokes 3 are in oneself
Before go to work, After Hours go home, then go to SR/SE/O places one to be gone home after plowing from family again, this active strokes contain two
Individual Trip chain, i.e. H-W and H-SR/SE/O.Three kinds of active strokes proportions of working day and as shown in table 3.
The Trip chain ratio of table 3
That is working day each car Trip chain C obeys following probability distribution:
P (C)=pi(i=1,2,3) (6)
X obeys Normal probability distribution between at the beginning of every section of stroke, and its probability density function is:
All electronic private cars are considered as totality, each electric automobile is considered as a two-dimensional random variable, and it, which is combined, divides
Cloth function is:
P (C ∩ x)=P (C) P (C | x) (8)
Trip chain type and departure time are extracted by the method for Monte Carlo simulation.Nonworkdays uses identical calculations side
Formula.
2.1.2 the random most short used time path computing of meter and real-time traffic
If each link length is Li, the travel speed of each road of t is Vi,t, each side right value is set as by each side
Required time Ti, then have:
Random most short used time path, the traveling of corresponding automobile different sections of highway are obtained by band random chance factor M DP algorithm
The relevant parameters such as speed, the spatial and temporal distributions of private car can be obtained.As shown in figure 3, grey circuit represents driving path, it is right respectively
Should be in each period.
2.2 officer's car
The trip rule of officer's car is relatively simple, general only to be moved between institutional settings or other places.Officer's car
Trip include two kinds of active strokes, i.e. W-W and W-SR/SE/O, at the beginning of every section of stroke between probability distribution it is identical, take
From normal distribution, the spatial and temporal distributions of officer's car are can obtain using the Monte-carlo Simulation Method similar with private car.
2.3 taxi
Daily 0:When 00, taxi starts to accept an order, driven to where client from affiliated taxi company
Ground, then the destination for going to client to go on a journey, receive new order, so circulation, as shown in Figure 4 again after arriving at.Daily 23:
Stop accepting an order after 00, return to affiliated company location and hand over to the next shift.
Each electric taxi is 0:When 00 from company location, destination is randomly selected, with private car calculating side
Method is similar, and most short used time path is obtained using the most short used time road of meter and real-time traffic.Runed at random until two between at work
The individual time of handing over to the next shift returns to taxi company.It is corresponding with electric taxi travel speed, each electric taxi can be obtained
Spatial and temporal distributions.
2.4 bus
Bus route is fixed, and time of departure obedience is uniformly distributed.For each bus, using Monte Carlo mould
The mode of plan extracts the time of departure and the parking of each website makes that the time, corresponding different sections of highway travel speed, bus can be obtained
Spatial and temporal distributions.
3 real-time traffics and temperature energy consumption model
Electric automobile energy consumption factor has a lot, and wherein temperature influences most with traffic on its unit mileage power consumption
Greatly.It is too high or too low for temperature that driver can be caused to open air-conditioning, increase electric automobile energy consumption.Automobile-shaped under different traffics
Formula speed is different, and automobile energy consumption is also different.
3.1 traffic energy consumption models
Speed, acceleration are had nothing in common with each other under different categories of roads and jam situation when vehicle travels, single caused by correspondence
Potential energy consumption is also different.Document《EV Energy Consumption modeling and continual mileage estimation research based on driving cycle》Using electronic
Automobile stand method energy consumption experimental data, has obtained the electric automobile energy consumption model of different categories of roads.
According to《Urban road engineering design specification》(CJJ37-2012), urban road should be according to road in road network
Status, communication function and to service function along the line etc., are divided into four through street, trunk roads, secondary distributor road and branch road grades, root
The data provided according to traffic department, are divided into unimpeded, substantially unimpeded, slight congestion, moderate congestion by congestion in road and seriously gather around
Stifled six kinds, the road of every kind of grade has different travel speeds under different jam situations, as shown in table 4:
Different congestion degree road speeds under 4 each category of roads of table
The measured data modeling of the energy consumed based on the identical mileage of different category of roads downward drivings, obtains corresponding to often
Its different periods and the specific energy consumption under the real-time jam situation of different road conditions, as shown in table 5, wherein Eksl、Ezgl、Ecgl、EzlRespectively
Represent through street, trunk roads, secondary distributor road, the electric automobile unit mileage power consumption (unit of branch road:KWh/km), V represents electronic
Automobile section average overall travel speed (unit:km/h).
Unit mileage power consumption under 5 each category of roads of table
For each electric automobile in region, certain position is obtained according to its Temporal And Spatial Distribution Model at moment, it is corresponding real
When traffic congestion after can obtain moment EV units mileage power consumption.
3.2 temperature energy consumption models
When outdoor temperature it is too high or too low so as to beyond the comfort zone of human body when, or environmental wet is spent under low temperature
When height need to open air-conditioning progress defrosting-defogging, electric automobile driver can typically select air-conditioning in opening vehicle.Document《Air conditioning for automobiles
Life cycle climate performance assessment models》It has studied air-conditioning unlatching rate of the domestic several big cities under different temperatures section, logarithm
According to the calculation formula that least square air-conditioning unlatching rate is fitted after processing:
Kpect=7 × 10-5·TEMP4-0.0887·TEMP3+3.673·TEMP2-56.302·TEMP+310.7681 (10)
Wherein TEMP represents temperature, KpectAir-conditioning unlatching rate is represented, draws out its functional image, as shown in Figure 5.
Document《Pure electric automobile performance estimating method and popularization feasibility study based on real example》To electric automobile air conditioner
Open the relation between unit mileage power consumption to be studied, define temperature energy consumption proportionality coefficient, be i.e. air-conditioning is different after opening
The ratio between unit distance power consumption corresponding to temperature and benchmark power consumption.Data measured by experiment in text are handled, obtained most
A young waiter in a wineshop or an inn multiplies fitting formula:
Ktemp=10-3·(TEMP-10)2+1.2 (11)
Wherein KtempFor temperature proportional coefficient, TEMP is temperature, draws out its functional image, as shown in Figure 6.
According to temperatures above energy consumption model, the real time temperature at corresponding estimation range moment, can obtain electronic in region
To the influence situation of electric automobile unit mileage power consumption after the ratio that air conditioning for automobiles is opened, and air-conditioning unlatching.
4 charging electric vehicle carry calculations
4.1 matrixes represent
Multiple nodes and moment for being related in one day is represented with the form of matrix, spatial and temporal distributions are unified, effectively receive
Collection information and it is easy to calculate and assesses, and matrix is than loop control statement more efficiently in mathematical software MatLab.Define m rows
The charging matrix of loadings PE of n row, each time interval is represented per a line, each row represent each node, and matrix element represents certain
The charging load at place corresponding moment, i.e. PEi,jRepresent the charge power at i-node j moment.The quantity m of time interval, i.e. matrix
In row quantity, depending on temporal resolution and calculate requirement.As shown in Figure 7.
4.2 all kinds of car charging Load Calculation Methods
(1) A classes private car and officer's car., can because A classes private car and A class officer's cars are fixed location trickle charge pattern
Charging carry calculation is carried out using identical method, is introduced by taking A class private cars as an example.A classes private car has the slow of fixation
Fast charging station, for its main charging modes to go home to carry out trickle charge daily, charging start time is to arrive at the house moment.Returned from car owner
Family's whole period of being set out by second day is charging interval section, until full of.Each electronic vapour can be obtained based on method above
Traveling spatial and temporal distributions, the power consumption situation of car, and then its charging load in place parking site can be obtained.To A class private savings
For car, officer's car, because its place of charging is fixed, charging load corresponds to certain of PE matrixes caused by each electric automobile
One row.
(2) B classes private car, officer's car.Because B classes private car and B class officer's cars are by way of ground fast charge pattern, can use
Identical method carries out charging carry calculation, is introduced by taking B class private cars as an example.It is gone on a journey since the departure place in path
Before, first judge whether battery dump energy EC meets the needs of this section of distance.Consider that environment temperature calculates list with traffic
Position mileage power consumption, corresponds to this section of path length and obtains this section of distance power consumption ER1.If dump energy is not enough to be unsatisfactory for reaching
Next destination, that is, be unsatisfactory for:
EC>ER1 (12)
Then electric automobile carries out a quick charge in current location, and charging duration tcr is extracted using Monte Carlo mode,
And it is limited within 60min.If the charging time started is TS, the charging end time is TE, then has:
TE=TS+tcr (13)
It is corresponding each with PE matrixes by way of the load that charged caused by charging place at it for B classes private car, officer's car
Row.
(3) taxi.Taxi is considered as A class private cars for a long time, without fixed destination traveling, and taxi is on way
The corresponding each row with PE matrixes of charging load caused by charged place.
(4) bus.The species of electric bus is many at present, and the charging modes of use also differ.It is electronic to embody
The space-time randomness of automobile distribution, it is assumed that all electric bus are the super capacity public transport by way of bus station quick charge
Car.Therefore bus is considered as the A class private cars that route is fixed, battery capacity and unit mileage power consumption are bigger, bus
In the corresponding each row with PE matrixes of load that charged caused by each website.
4.3 simulation flows are as shown in Figure 8.
5 sample calculation analysis
5.1 transportation network
As shown in figure 9, by taking Jincheng City, Shanxi Province as an example, by the abstract traffic containing 16 nodes, 20 sides of its transportation network
Network, wherein each side institute band color represents its category of roads, node 1-6 is residential block (H), and node 11-15 is trading estate
(W), node 7,8,9,10,16 is that amusement is gone shopping and rested area (SR/SE/O).According to the above-mentioned setting to city electrically-charging equipment, 6
H places and 5 W places have fast charge stake and trickle charge stake simultaneously, that is, are provided simultaneously with fast charge function and trickle charge function, 5 SR/SE/O
Place sets fast charge stake.
5.2 recoverable amounts are predicted
Domestic the year two thousand thirty car ownership prediction result is as shown in table 6.
The domestic year electric automobile recoverable amount of table 6 2030
With Jincheng City GDP compared with domestic GDP total amounts, to obtain the electric automobile recoverable amount such as institute of table 7 of the year two thousand thirty Jincheng City
Show:
7 2030 Jincheng City year of table car ownership
5.3 simulation analysis
5.3.1 parameter setting
It is 100% to set electric automobile permeability, is networked with simulating extensive electric automobile.Private car, taxi and public affairs
Business car is generally compact car, using long easy dynamic EV as standard, battery capacity 30kWh, trickle charge power 7kW, fast charge power 25kW.It is public
Friendship car is medium-and-large-sized vehicle, and using SWB612SC super capacity public transports car as standard, its capacitive energy is 5kWh, and the charging interval is
1min, charge power 300kW.Definition typical week is the degrees centigrade of temperature 25, traffic well one week.Define high temperature week
It is temperature more than 35 degree, other conditions and typical all identicals one week.Congestion week is traffic heavy congestion in peak period on and off duty, its
His condition and typical Zhou Xiang identicals one week.Typical week, high temperature week, congestion the inside of a week, nonworkdays active strokes, temperature are handed over
Logical and its corresponding electric automobile quantity provides in annex.Emulation start time assumes that all batteries of electric automobile SOC are
100%.Setting simulation step length is 1min, and simulation time is 20 weeks.
5.3.2 simulation result
(1) each node charging load of typical case's week lower Jincheng City A, B class car.
As shown in Figure 10, simulation result shows, for having the node of fast charge and trickle charge ability such as 11-16 simultaneously, charges
Load is one week to be that tendency is relatively gentle from the point of view of the cycle, peak-valley ratio is less than 30%.The main node such as 1 for undertaking fast charge function,
7th, 8,9,10 charging load fluctuations are larger, and peak-valley ratio is higher, and average peak-valley ratio is more than 75%, and highest is more than 90%.
(2) each node charging load of typical case's week lower Jincheng City bus.
As shown in figure 11, simulation result shows that super capacity public transport car is relatively equal in each website charging power load distributing daily
Even, night parking rest does not produce charging load.
(3) each node charging load of typical case's week lower Jincheng City taxi.
As shown in figure 12, simulation result shows that taxi produces charging load in each website each time, and in morning
Hand over to the next shift the period produce charging peak, it is similar to internal combustion engine taxi morning oiling characteristic.
(4) typical case's week lower Jincheng City charging electric vehicle total load.By all kinds of charging electric vehicle load superpositions of each node.
As shown in figure 13, simulation result shows, the class car of A, B two for accounting for absolute quantity determines the general trend of charging load.
(5) high temperature week lower Jincheng City charging electric vehicle load.Typical week is subjected to lateral comparison with high temperature Zhou Jiedian 10;
As shown in figure 14, simulation result shows, when temperature keeps a high position that node charging load peak can be caused to hold in one week
Continuous time increase, more typical Zhou Pingjun increase about 2.5 hours daily.
(3) under congestion serious conditions, Jincheng City charging electric vehicle load.Typical week is entered with congestion Zhou Jiedian 15
Row lateral comparison;
As shown in figure 15, simulation result shows, on the one hand node can be caused to fill after traffic congestion deterioration in city
The increase of electric load peak period, more typical Zhou Pingjun increase about 3 hours daily.On the other hand load peak can be caused to increase
Height, wherein load peak increase about 550kW.It can also be seen that traffic change also has shadow to duration of load application distribution
Ring, the peak for the load that charges shifts to an earlier date, pre-set time 15-20min.
(6) with the contrast of other method prediction result.Electric automobile quantity is arranged to 240,000 amounts, simulation step length is arranged to
Min, simulation time are arranged to 1440min (one day), and other simulation parameters are with document《Charging electric vehicle based on Trip chain is born
Lotus forecast model》(Chen Lidan, Nie Yongquan, clock celebrating [J] electrotechnics journals, 2015, (04):216-225.) it is configured.
Each node is charged, and load is superimposed to obtain region-wide total charging load condition, is contrasted with document simulation result;
As shown in figure 16, the prediction data of two methods is contrasted it can be found that there are two load peaks in one day, it is corresponding
In period at noon and the dusk period.Wherein peak value tendency is closer to the dusk period, and the inventive method period at noon peak value is bright
Aobvious to be higher than bibliography method, reason is:The spatial distribution of electric automobile is described as the phase in 5 big regions by above-mentioned literature method
Mutually transfer, majority charging moment occur in Trip chain ending phase.The inventive method simulate in the region of electric automobile when
Idle running moves, and simulation result finds to occur that fairly large electric automobile is quickly filled in job site around noon daily
Electricity.
Finally illustrate, preferred embodiment above is merely illustrative of the technical solution of the present invention and unrestricted, although logical
Cross above preferred embodiment the present invention is described in detail, it is to be understood by those skilled in the art that can be
Various changes are made to it in form and in details, without departing from claims of the present invention limited range.
Claims (6)
1. meter and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodology, it is characterised in that:This method include with
Lower step:
S1:First all kinds of electric automobiles are classified, analyze its space transfer and charge characteristic;
S2:Make abstract simulation to urban traffic network with Graph-theoretical Approach using Trip chain is theoretical;
S3:All kinds of electric automobiles are set up using Monte Carlo simulation and the random shortest path method based on Markovian decision
Time-spatial diversion model;
S4:Traffic and temperature energy consumption model are established according to electric automobile actual test data.
2. meter as claimed in claim 1 and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodologies, its feature
It is:The step S1 is specially:
Charging modes are fixed parking stall trickle charge, and trip mode is come out early and return late for working day, stopped over, and day off steps out
EV types be private car A classes;
Charging modes are to be come out early and return late by way of ground fast charge, trip mode for working day, are stopped over, the EV that day off steps out
Type is private car B classes;
Charging modes are fixed parking stall trickle charge, and trip mode is that working day traffic path is fixed, the EV classes that day off does not go on a journey
Type is officer's car A classes;
Charging modes are by way of ground fast charge, and trip mode is that working day traffic path is fixed, and the EV types that day off does not go on a journey are
Officer's car B classes;
Charging modes are by way of bus stop quick charge, and trip mode is that traffic path is fixed with the travel time, whole year operation
EV types be bus;
Charging modes are to be fixed by way of ground fast charge station quick charge, trip mode to set out with the off-running time, and route is random, annual
The EV types of operation are taxi.
3. meter as claimed in claim 1 and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodologies, its feature
It is:The step S2 is specially:
Trip chain is divided into simple chain and compound catenary;Simple chain only has one for the trip purpose of every Trip chain;Compound catenary is
Trip chain contains multiple trip purposes;It is divided into residence H, place of working W, shopping use by resident trip survey data to trip purpose
Eat SE, social entertainment SR, other O, while is also charging electric vehicle place;
Wherein, the electrically-charging equipment that residence H is set is fast charge+trickle charge;
The electrically-charging equipment that place of working W is set is fast charge+trickle charge;
The electrically-charging equipment that shopping dining SE is set is fast charge;
The electrically-charging equipment that social entertainment SR is set is fast charge;
The electrically-charging equipment that other O are set is fast charge;
Urban traffic network is abstracted as by digraph using the method for graph theory, each highway in city is arranged to side, will be handed over
Cross road mouth is arranged to node, and road length is arranged to the weight on side, and transportation network is abstracted as to the set of side and node.
4. meter as claimed in claim 1 and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodologies, its feature
It is:The step S3 is specially:
For each EV in region, it is known that after its starting point, using the Markovian decision with random chance factor gamma
Process obtains the optimal path with randomness in its space;
Wherein, the Markovian decision process is the state observed according to each moment, is selected from available action collection
Made decisions with an action, the state of system next step is random, and its state transition probability has Markov property;
Markovian decision process is:
M=(T, S, A (i), p, R)
Wherein, T is decision-making moment collection, represents to choose the time point set of action;S is state set, represent at each decision-making moment, be
System state space collection;A (i) collects for action, represents at any one decision-making moment, the state that policymaker observes, i.e., in state i
The available action collection chosen;P (j | i, it is a) state transition probability, represents under state i, a that takes action is transferred to state j's
Probability;R (i, it is a) Reward Program, represents under state i, after a that takes action, the return of policymaker's acquisition;
That is Markovian decision dynamic process is as follows:Policymaker state S (0) where T (0) moment, the selection one from A (0)
Act a (0) to perform, arrived next S (1) state by P (S (1) | S (0), a (0)) probability random transferring after execution, then held again
One action a (1) of row, successively backward, to the last moment T (N), reaches end-state S (N);
It is the rule that a decision system selects action on each different conditions to define decision rule, is designated as f (S (0)), f (S
(1)),f(S(2))…f(N);Definition strategy is the sequence of decision rule, is designated as π=(f (S (0)), f (S (1)), f (S (2)) ...
F ((S (N))), the collection of All Policies is collectively referred to as policy class, is designated as Π;From original state S (0) simultaneously used tactful π ∈ Π,
So thatMaximum or minimum value is obtained, then tactful π is optimal policy;Define value function VπIn state S during to use tactful π
(i) expected returns:
<mrow>
<msup>
<mi>V</mi>
<mi>&pi;</mi>
</msup>
<mo>=</mo>
<mi>E</mi>
<mo>&lsqb;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>R</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>S</mi>
<mi>t</mi>
</msup>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>,</mo>
<mi>&pi;</mi>
<mo>(</mo>
<mrow>
<msup>
<mi>S</mi>
<mi>t</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
Wherein E is it is expected, StFor t status, t=0 is corresponded to original state S (0), is expressed as with recursive form:
<mrow>
<msup>
<mi>V</mi>
<mi>&pi;</mi>
</msup>
<mo>=</mo>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>,</mo>
<mi>&pi;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
</munderover>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>,</mo>
<mi>S</mi>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msup>
<mi>V</mi>
<mi>&pi;</mi>
</msup>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
To each tactful π, its corresponding value function VπA series of unique public solution of linear equations, optimal policy by based on
Recursive method obtains the expected return value of Dynamic Programming backward;
All node sets in region are considered as S, each node time instance passed through from origin-to-destination traveling is considered as the decision-making moment
t;Determine that next paths of traveling are considered as action a in each node, the path length that will be travelled is considered as return R;
If EV is travelled with part probability transfer mode under steam, it is equally distributed random general to introduce the obedience that value is 0~0.5
Rate factor gamma, represent, when policymaker is often performed with this action, there is probability γ to be transferred to another paths;By in Path selection
Shortest path is combined with randomness, is not influenceed by the original state and initial action set out, surplus for next decision-making
The strategy of remaining decision rule composition is optimal policy;
(1) for private car:
Working day, private car active strokes share 3 kinds:Active strokes 1:H-W;Active strokes 2:H-W-SR/SE/O;Active strokes
3:H-W+H-SR/SE/O;It is each to obey certain proportion;Working day each car Trip chain C obeys probability distribution:P (C)=pi(i=1,
2,3)
X obeys Normal probability distribution between at the beginning of every section of stroke, and its probability density function is:
<mrow>
<msub>
<mi>f</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msqrt>
<msub>
<mi>&sigma;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo><</mo>
<mi>x</mi>
<mo>&le;</mo>
<msub>
<mi>&mu;</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mn>12</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mrow>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msqrt>
<msub>
<mi>&sigma;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>+</mo>
<mn>24</mn>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>e</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&mu;</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mn>12</mn>
<mo><</mo>
<mi>x</mi>
<mo>&le;</mo>
<mn>24</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Wherein σeFor Variance of Normal Distribution, μeFor normal distribution average;
All electronic private cars are considered as totality, each electric automobile is considered as a two-dimensional random variable, its Joint Distribution letter
Number is:P (C ∩ x)=P (C) P (C | x), Trip chain type and departure time are extracted by the method for Monte Carlo simulation;
The random most short used time path computing of meter and real-time traffic:If each link length is Li, the travel speed of each road of t
For Vi,t, each side right value is set as by T the time required to each sidei, then havePass through band random chance factor M DP algorithm
The relevant parameter such as random most short used time path, the travel speed of corresponding automobile different sections of highway is obtained, obtains the when space division of private car
Cloth;
(2) for officer's car
Probability distribution between at the beginning of every section of stroke is identical, equal Normal Distribution, using the Meng Teka similar with private car
Lip river analogy method obtains the spatial and temporal distributions of officer's car;
(3) for taxi
Each taxi is 0:When 00 from company location, randomly select destination, using meter and real-time traffic it is most short
Used time, road obtained most short used time path;It is random between at work to run time return taxi company of being handed over to the next shift until two, with electricity
Dynamic taxi travel speed is corresponding, obtains the spatial and temporal distributions of each electric taxi;
(4) for bus
Time of departure obedience is uniformly distributed, for each bus, using extracted by the way of Monte Carlo simulation dispatch a car when
Between with the parking of each website make the time, corresponding different sections of highway travel speed, obtain the spatial and temporal distributions of bus.
5. meter as claimed in claim 1 and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodologies, its feature
It is:The step S4 is specially:
S401:Traffic energy consumption model is established according to electric automobile actual test data:
The measured data modeling of the energy consumed based on the identical mileage of different category of roads downward drivings, obtains corresponding to daily not
With period and the specific energy consumption under the real-time jam situation of different road conditions, through street unit mileage power consumption is Eksl=0.247+
1.520/V-0.004*V+2.992*10(-5)*V;
Trunk roads unit mileage power consumption is Ezgl=-0.179+0.004*V+5.492/V;
Secondary distributor road unit mileage power consumption is Ecgl=0.21-0.001*V+1.531/V;
Branch road unit mileage power consumption is Ezl=0.208-0.002*V+1.553/V;
Wherein Eksl、Ezgl、Ecgl、EzlThe electric automobile unit mileage power consumption of through street, trunk roads, secondary distributor road, branch road is represented respectively
Amount, V represent electric automobile road-section average travel speed;
For each electric automobile in region, certain position is obtained according to its Temporal And Spatial Distribution Model at moment, it is corresponding to hand in real time
Moment EV units mileage power consumption is obtained after logical congestion;
S402:Temperature energy consumption model is established according to electric automobile actual test data:
The calculation formula of least square air-conditioning unlatching rate is:
Kpect=7 × 10-5·TEMP4-0.0887·TEMP3+3.673·TEMP2-56.302·TEMP+310.7681
Wherein TEMP represents temperature, KpectAir-conditioning unlatching rate is represented, draws out its functional image;
Temperature energy consumption proportionality coefficient is defined, i.e., unit distance power consumption corresponding to different temperatures and benchmark power consumption after air-conditioning unlatching
The ratio between, the fitting formula of least square temperature proportional coefficient is:
Ktemp=10-3·(TEMP-10)2+1.2
Wherein KtempFor temperature proportional coefficient, TEMP is temperature, draws out its functional image;
According to temperatures above energy consumption model, the real time temperature at corresponding estimation range moment, electric automobile air conditioner in region is obtained
To the influence situation of electric automobile unit mileage power consumption after the ratio of unlatching, and air-conditioning unlatching.
6. meter as claimed in claim 1 and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodologies, its feature
It is:The charging Load Calculation Method of all kinds of electric automobiles:
The charging matrix of loadings PE of m rows n row is defined, each time interval is represented per a line, each row represent each node, matrix
The charging load at element representation place corresponding moment, i.e. PEi,jRepresent the charge power at i-node j moment;The number of time interval
M, i.e. row in matrix quantity are measured, depending on temporal resolution and the requirement calculated;
(1) A classes private car and officer's car:According to the traveling spatial and temporal distributions of each electric automobile, charging modes, it is obtained at place
The charging load of parking site;For A classes private car, officer's car, because its place of charging is fixed, each electric automobile institute
Caused charging load corresponds to a certain row of PE matrixes;
(2) B classes private car, officer's car:Consider that environment temperature calculates unit mileage power consumption with traffic, correspond to this Duan Lu
Cheng Changdu obtains this section of distance power consumption ER1;If dump energy EC is not enough to be unsatisfactory for reaching next destination, that is, it is unsatisfactory for
EC>ER1, then electric automobile current location carry out a quick charge, charging duration tcr extracted using Monte Carlo mode,
And it is limited within 30min;If the charging time started is TS, the charging end time is TE, then has:TE=TS+tcr;
It is by needing to judge that can dump energy meet road during next place with A class private car calculation process difference
Journey needs;It is corresponding each with PE matrixes by way of the load that charged caused by charging place at it for B classes private car, officer's car
Row;
(3) taxi:Taxi is A class private cars for a long time, without fixed destination traveling, and taxi is by way of charging place
The corresponding each row with PE matrixes of caused charging load;
(4) bus:Bus is the A class private cars that route is fixed, battery capacity and unit mileage power consumption are bigger, bus
In the corresponding each row with PE matrixes of load that charged caused by each website.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710787528.1A CN107392400B (en) | 2017-09-04 | 2017-09-04 | EV charging load space-time distribution prediction method considering real-time traffic and temperature |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710787528.1A CN107392400B (en) | 2017-09-04 | 2017-09-04 | EV charging load space-time distribution prediction method considering real-time traffic and temperature |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107392400A true CN107392400A (en) | 2017-11-24 |
CN107392400B CN107392400B (en) | 2020-10-16 |
Family
ID=60347903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710787528.1A Active CN107392400B (en) | 2017-09-04 | 2017-09-04 | EV charging load space-time distribution prediction method considering real-time traffic and temperature |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107392400B (en) |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108288112A (en) * | 2018-01-30 | 2018-07-17 | 广西电网有限责任公司柳州供电局 | Region electric automobile charging station load forecasting method based on user's trip simulation |
CN108510128A (en) * | 2018-04-11 | 2018-09-07 | 华南理工大学广州学院 | A kind of region electric vehicle charging load spatial and temporal distributions prediction technique |
CN108981736A (en) * | 2018-05-28 | 2018-12-11 | 国网江苏省电力有限公司电力科学研究院 | A kind of regular electric car charge path optimization method of being gone on a journey based on user |
CN109308541A (en) * | 2018-08-02 | 2019-02-05 | 华南理工大学 | Electric vehicle charging load prediction method based on scene analysis |
CN109409939A (en) * | 2018-10-08 | 2019-03-01 | 国网天津市电力公司电力科学研究院 | For the voltage-controlled double-layer electric automobile quick charge electricity price modification method of distribution |
CN109583708A (en) * | 2018-11-08 | 2019-04-05 | 国网浙江省电力有限公司经济技术研究院 | The method for building up of multiple agent microcosmic traffic Assignment Model |
CN110466384A (en) * | 2019-07-05 | 2019-11-19 | 武汉新能源汽车工业技术研究院有限公司 | A kind of charging module group power distribution method and device |
CN110516372A (en) * | 2019-08-29 | 2019-11-29 | 重庆大学 | Meter and the electric car state-of-charge spatial and temporal distributions analogy method of Quasi dynamic traffic flow |
CN110570014A (en) * | 2019-08-07 | 2019-12-13 | 浙江大学 | Electric vehicle charging load prediction method based on Monte Carlo and deep learning |
CN110728396A (en) * | 2019-09-25 | 2020-01-24 | 国网山西省电力公司经济技术研究院 | Electric vehicle charging load comprehensive modeling method considering space-time distribution |
CN110956397A (en) * | 2019-12-03 | 2020-04-03 | 深圳供电局有限公司 | Method for quantitatively evaluating influence of network attack on voltage of power distribution network |
CN110956329A (en) * | 2019-12-02 | 2020-04-03 | 国网浙江省电力有限公司绍兴供电公司 | Load prediction method based on distributed photovoltaic and electric vehicle space-time distribution |
CN111199300A (en) * | 2018-10-31 | 2020-05-26 | 国网江苏省电力有限公司南京供电分公司 | Electric vehicle charging load space-time prediction method under vehicle-road-network mode |
CN111216730A (en) * | 2020-01-15 | 2020-06-02 | 山东理工大学 | Method, device, storage medium and equipment for estimating remaining driving range of electric automobile |
CN111291999A (en) * | 2020-02-21 | 2020-06-16 | 东南大学 | Dynamic wireless charging road section load modeling method based on hierarchical clustering |
CN111368435A (en) * | 2020-03-05 | 2020-07-03 | 浙江大学 | Uncertain demand-oriented air separation pipe network device starting and stopping and load scheduling method |
CN111428946A (en) * | 2020-04-28 | 2020-07-17 | 三峡大学 | Distributed optimal scheduling method for supply side of charging and storage station |
CN112330025A (en) * | 2020-11-06 | 2021-02-05 | 国网冀北电力有限公司张家口供电公司 | Prediction method of space-time charging load for urban electric vehicle |
CN112428834A (en) * | 2020-11-17 | 2021-03-02 | 宁波工程学院 | Monte Carlo method-based intelligent electric vehicle charging optimization method and system |
CN113094852A (en) * | 2021-03-31 | 2021-07-09 | 东北电力大学 | Electric vehicle charging load time-space distribution calculation method |
CN113191523A (en) * | 2021-01-27 | 2021-07-30 | 国电南瑞南京控制系统有限公司 | Urban electric vehicle rapid charging demand prediction method and device based on data driving mode and behavior decision theory |
CN113837498A (en) * | 2021-11-08 | 2021-12-24 | 华北电力大学 | Charging facility scale prediction method and system |
CN114037177A (en) * | 2021-11-22 | 2022-02-11 | 山东德佑电气股份有限公司 | Bus charging load optimization method in crowded traffic state based on train number chain |
CN114282377A (en) * | 2021-12-28 | 2022-04-05 | 国网四川省电力公司电力科学研究院 | Electric vehicle short-term charging load prediction method considering meteorological factors |
CN114676885A (en) * | 2022-03-02 | 2022-06-28 | 三峡大学 | Electric vehicle charging and discharging load space-time distribution prediction method |
CN114819345A (en) * | 2022-04-25 | 2022-07-29 | 重庆大学 | Electric vehicle charging load space-time distribution prediction method based on vehicle-road-station-network fusion |
CN115049167A (en) * | 2022-08-16 | 2022-09-13 | 北京市城市规划设计研究院 | Traffic situation prediction method, device, equipment and storage medium |
CN115147244A (en) * | 2022-07-21 | 2022-10-04 | 东北电力大学 | Method for achieving wind curtailment and accommodation by considering charging load-electricity price response of electric automobile |
CN115275999A (en) * | 2022-08-09 | 2022-11-01 | 东北电力大学 | Power distribution network optimal scheduling method considering electric automobile time-varying road impedance |
WO2023109099A1 (en) * | 2021-12-15 | 2023-06-22 | 重庆邮电大学 | Charging load probability prediction system and method based on non-intrusive detection |
CN117494882A (en) * | 2023-11-01 | 2024-02-02 | 吉林大学 | Urban multi-scene charging load prediction method based on vehicle operation background data |
CN117829378A (en) * | 2024-03-04 | 2024-04-05 | 华东交通大学 | Track traffic energy consumption prediction method based on space-time data |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120233062A1 (en) * | 2011-03-11 | 2012-09-13 | Kevin Terrill Cornish | Method and process for an electric vehicle charging system to automatically engage the charging apparatus of an electric vehicle |
CN106355290A (en) * | 2016-09-21 | 2017-01-25 | 四川大学 | Electric vehicle charge load prediction method and system based on Markov chain |
CN106599390A (en) * | 2016-11-23 | 2017-04-26 | 国网浙江省电力公司电动汽车服务分公司 | Method for calculating charging load considering time-space stochastic characteristics of electric taxi |
CN107067110A (en) * | 2017-04-14 | 2017-08-18 | 天津大学 | Charging electric vehicle load spatio-temporal prediction method under car Road network pattern |
-
2017
- 2017-09-04 CN CN201710787528.1A patent/CN107392400B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120233062A1 (en) * | 2011-03-11 | 2012-09-13 | Kevin Terrill Cornish | Method and process for an electric vehicle charging system to automatically engage the charging apparatus of an electric vehicle |
CN106355290A (en) * | 2016-09-21 | 2017-01-25 | 四川大学 | Electric vehicle charge load prediction method and system based on Markov chain |
CN106599390A (en) * | 2016-11-23 | 2017-04-26 | 国网浙江省电力公司电动汽车服务分公司 | Method for calculating charging load considering time-space stochastic characteristics of electric taxi |
CN107067110A (en) * | 2017-04-14 | 2017-08-18 | 天津大学 | Charging electric vehicle load spatio-temporal prediction method under car Road network pattern |
Non-Patent Citations (3)
Title |
---|
宋媛媛: "基于行驶工况的纯电动汽车能耗建模及续驶里程估算研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
陈丽丹 等: "基于出行链的电动汽车充电负荷预测模型", 《电工技术学报》 * |
陈娟 等: "基于模糊随机理论的驾驶人路径选择模型初步探讨", 《机械》 * |
Cited By (50)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108288112B (en) * | 2018-01-30 | 2020-08-21 | 广西电网有限责任公司柳州供电局 | Regional electric vehicle charging station load prediction method based on user travel simulation |
CN108288112A (en) * | 2018-01-30 | 2018-07-17 | 广西电网有限责任公司柳州供电局 | Region electric automobile charging station load forecasting method based on user's trip simulation |
CN108510128B (en) * | 2018-04-11 | 2022-04-15 | 华南理工大学广州学院 | Method for predicting space-time distribution of charging load of regional electric automobile |
CN108510128A (en) * | 2018-04-11 | 2018-09-07 | 华南理工大学广州学院 | A kind of region electric vehicle charging load spatial and temporal distributions prediction technique |
CN108981736A (en) * | 2018-05-28 | 2018-12-11 | 国网江苏省电力有限公司电力科学研究院 | A kind of regular electric car charge path optimization method of being gone on a journey based on user |
CN108981736B (en) * | 2018-05-28 | 2021-03-16 | 国网江苏省电力有限公司电力科学研究院 | Electric vehicle charging path optimization method based on user travel rule |
CN109308541A (en) * | 2018-08-02 | 2019-02-05 | 华南理工大学 | Electric vehicle charging load prediction method based on scene analysis |
CN109409939B (en) * | 2018-10-08 | 2022-02-22 | 国网天津市电力公司电力科学研究院 | Double-layer electric automobile rapid charging electricity price correction method for distribution network voltage control |
CN109409939A (en) * | 2018-10-08 | 2019-03-01 | 国网天津市电力公司电力科学研究院 | For the voltage-controlled double-layer electric automobile quick charge electricity price modification method of distribution |
CN111199300A (en) * | 2018-10-31 | 2020-05-26 | 国网江苏省电力有限公司南京供电分公司 | Electric vehicle charging load space-time prediction method under vehicle-road-network mode |
CN109583708A (en) * | 2018-11-08 | 2019-04-05 | 国网浙江省电力有限公司经济技术研究院 | The method for building up of multiple agent microcosmic traffic Assignment Model |
CN109583708B (en) * | 2018-11-08 | 2021-06-04 | 国网浙江省电力有限公司经济技术研究院 | Method for establishing multi-agent microscopic traffic distribution model |
CN110466384A (en) * | 2019-07-05 | 2019-11-19 | 武汉新能源汽车工业技术研究院有限公司 | A kind of charging module group power distribution method and device |
CN110570014A (en) * | 2019-08-07 | 2019-12-13 | 浙江大学 | Electric vehicle charging load prediction method based on Monte Carlo and deep learning |
CN110570014B (en) * | 2019-08-07 | 2022-04-01 | 浙江大学 | Electric vehicle charging load prediction method based on Monte Carlo and deep learning |
CN110516372B (en) * | 2019-08-29 | 2023-02-17 | 重庆大学 | Electric vehicle charge state space-time distribution simulation method considering quasi-dynamic traffic flow |
CN110516372A (en) * | 2019-08-29 | 2019-11-29 | 重庆大学 | Meter and the electric car state-of-charge spatial and temporal distributions analogy method of Quasi dynamic traffic flow |
CN110728396A (en) * | 2019-09-25 | 2020-01-24 | 国网山西省电力公司经济技术研究院 | Electric vehicle charging load comprehensive modeling method considering space-time distribution |
CN110956329B (en) * | 2019-12-02 | 2023-04-25 | 国网浙江省电力有限公司绍兴供电公司 | Load prediction method based on distributed photovoltaic and electric automobile space-time distribution |
CN110956329A (en) * | 2019-12-02 | 2020-04-03 | 国网浙江省电力有限公司绍兴供电公司 | Load prediction method based on distributed photovoltaic and electric vehicle space-time distribution |
CN110956397A (en) * | 2019-12-03 | 2020-04-03 | 深圳供电局有限公司 | Method for quantitatively evaluating influence of network attack on voltage of power distribution network |
CN111216730A (en) * | 2020-01-15 | 2020-06-02 | 山东理工大学 | Method, device, storage medium and equipment for estimating remaining driving range of electric automobile |
CN111291999A (en) * | 2020-02-21 | 2020-06-16 | 东南大学 | Dynamic wireless charging road section load modeling method based on hierarchical clustering |
CN111368435B (en) * | 2020-03-05 | 2022-03-22 | 浙江大学 | Uncertain demand-oriented air separation pipe network device starting and stopping and load scheduling method |
CN111368435A (en) * | 2020-03-05 | 2020-07-03 | 浙江大学 | Uncertain demand-oriented air separation pipe network device starting and stopping and load scheduling method |
CN111428946A (en) * | 2020-04-28 | 2020-07-17 | 三峡大学 | Distributed optimal scheduling method for supply side of charging and storage station |
CN111428946B (en) * | 2020-04-28 | 2022-07-01 | 三峡大学 | Distributed optimal scheduling method for supply side of charging and storage station |
CN112330025A (en) * | 2020-11-06 | 2021-02-05 | 国网冀北电力有限公司张家口供电公司 | Prediction method of space-time charging load for urban electric vehicle |
CN112428834A (en) * | 2020-11-17 | 2021-03-02 | 宁波工程学院 | Monte Carlo method-based intelligent electric vehicle charging optimization method and system |
CN113191523A (en) * | 2021-01-27 | 2021-07-30 | 国电南瑞南京控制系统有限公司 | Urban electric vehicle rapid charging demand prediction method and device based on data driving mode and behavior decision theory |
CN113094852B (en) * | 2021-03-31 | 2023-06-09 | 东北电力大学 | Electric automobile charging load time-space distribution calculation method |
CN113094852A (en) * | 2021-03-31 | 2021-07-09 | 东北电力大学 | Electric vehicle charging load time-space distribution calculation method |
CN113837498A (en) * | 2021-11-08 | 2021-12-24 | 华北电力大学 | Charging facility scale prediction method and system |
CN113837498B (en) * | 2021-11-08 | 2024-01-30 | 华北电力大学 | Charging facility scale prediction method and system |
CN114037177B (en) * | 2021-11-22 | 2024-05-14 | 山东德佑电气股份有限公司 | Bus charging load optimization method based on train number chain in crowded traffic state |
CN114037177A (en) * | 2021-11-22 | 2022-02-11 | 山东德佑电气股份有限公司 | Bus charging load optimization method in crowded traffic state based on train number chain |
WO2023109099A1 (en) * | 2021-12-15 | 2023-06-22 | 重庆邮电大学 | Charging load probability prediction system and method based on non-intrusive detection |
CN114282377A (en) * | 2021-12-28 | 2022-04-05 | 国网四川省电力公司电力科学研究院 | Electric vehicle short-term charging load prediction method considering meteorological factors |
CN114676885A (en) * | 2022-03-02 | 2022-06-28 | 三峡大学 | Electric vehicle charging and discharging load space-time distribution prediction method |
CN114676885B (en) * | 2022-03-02 | 2024-07-09 | 三峡大学 | Prediction method for charge-discharge load space-time distribution of electric automobile |
CN114819345A (en) * | 2022-04-25 | 2022-07-29 | 重庆大学 | Electric vehicle charging load space-time distribution prediction method based on vehicle-road-station-network fusion |
CN114819345B (en) * | 2022-04-25 | 2024-05-10 | 重庆大学 | Electric vehicle charging load space-time distribution prediction method based on vehicle-road-station-network integration |
CN115147244A (en) * | 2022-07-21 | 2022-10-04 | 东北电力大学 | Method for achieving wind curtailment and accommodation by considering charging load-electricity price response of electric automobile |
CN115275999A (en) * | 2022-08-09 | 2022-11-01 | 东北电力大学 | Power distribution network optimal scheduling method considering electric automobile time-varying road impedance |
CN115275999B (en) * | 2022-08-09 | 2024-07-05 | 东北电力大学 | Power distribution network optimal scheduling method considering electric automobile time-varying road impedance |
CN115049167A (en) * | 2022-08-16 | 2022-09-13 | 北京市城市规划设计研究院 | Traffic situation prediction method, device, equipment and storage medium |
CN117494882A (en) * | 2023-11-01 | 2024-02-02 | 吉林大学 | Urban multi-scene charging load prediction method based on vehicle operation background data |
CN117494882B (en) * | 2023-11-01 | 2024-05-24 | 吉林大学 | Urban multi-scene charging load prediction method based on vehicle operation background data |
CN117829378A (en) * | 2024-03-04 | 2024-04-05 | 华东交通大学 | Track traffic energy consumption prediction method based on space-time data |
CN117829378B (en) * | 2024-03-04 | 2024-05-14 | 华东交通大学 | Track traffic energy consumption prediction method based on space-time data |
Also Published As
Publication number | Publication date |
---|---|
CN107392400B (en) | 2020-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107392400A (en) | Meter and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodology | |
CN109034648B (en) | Electric vehicle cluster demand response potential evaluation method | |
CN110728396B (en) | Electric vehicle charging load comprehensive modeling method considering space-time distribution | |
Yi et al. | Energy impact evaluation for eco-routing and charging of autonomous electric vehicle fleet: Ambient temperature consideration | |
Tao et al. | Data-driven optimized layout of battery electric vehicle charging infrastructure | |
Yi et al. | Data-driven optimal charging decision making for connected and automated electric vehicles: A personal usage scenario | |
CN106599390B (en) | It is a kind of meter and electric taxi space-time stochastic behaviour charging load calculation method | |
Sayarshad et al. | Non-myopic dynamic routing of electric taxis with battery swapping stations | |
Grahn et al. | PHEV utilization model considering type-of-trip and recharging flexibility | |
CN109711630A (en) | A kind of electric car fast charge station addressing constant volume method based on trip probability matrix | |
CN108681795B (en) | Electric automobile charging load space-time prediction method under constraint of urban traffic network and user travel chain | |
Jiang et al. | Multiple agents and reinforcement learning for modelling charging loads of electric taxis | |
CN109934403A (en) | Charge load Analysis prediction technique in electric car resident region based on mathematical model | |
CN105160428A (en) | Planning method of electric vehicle fast-charging station on expressway | |
CN110674575A (en) | Electric vehicle cluster charging demand and discharging capacity model modeling method based on time sequence traveling set | |
CN105528498B (en) | Net connection intelligent electric vehicle integrated modelling and integrated control method | |
CN106875075A (en) | A kind of electric automobile charging station points distributing method based on travel behaviour | |
Acha et al. | Modelling spatial and temporal agent travel patterns for optimal charging of electric vehicles in low carbon networks | |
Su et al. | Forecast of electric vehicle charging demand based on traffic flow model and optimal path planning | |
CN108062591A (en) | Electric vehicle charging load spatial and temporal distributions Forecasting Methodology | |
CN110232219A (en) | A kind of schedulable capacity ratification method of electric car based on data mining | |
CN113112097A (en) | Electric vehicle load prediction and charging facility layout optimization method | |
Boulakhbar et al. | Electric vehicles arrival and departure time prediction based on deep learning: the case of Morocco | |
CN105160418A (en) | Charging distribution predication method based on electric vehicle application features | |
Ma et al. | A model for electric vehicle charging load forecasting based on simulated driving path |
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 |