CN106355290A - Electric vehicle charge load prediction method and system based on Markov chain - Google Patents
Electric vehicle charge load prediction method and system based on Markov chain Download PDFInfo
- Publication number
- CN106355290A CN106355290A CN201610836365.7A CN201610836365A CN106355290A CN 106355290 A CN106355290 A CN 106355290A CN 201610836365 A CN201610836365 A CN 201610836365A CN 106355290 A CN106355290 A CN 106355290A
- Authority
- CN
- China
- Prior art keywords
- charging
- charge
- soc
- end time
- state value
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000008569 process Effects 0.000 claims abstract description 16
- 238000004088 simulation Methods 0.000 claims abstract description 13
- 230000005611 electricity Effects 0.000 claims description 16
- 238000013277 forecasting method Methods 0.000 claims description 9
- 238000000605 extraction Methods 0.000 claims description 2
- 230000006399 behavior Effects 0.000 abstract description 40
- 238000000342 Monte Carlo simulation Methods 0.000 abstract description 6
- 230000000007 visual effect Effects 0.000 abstract 1
- 101100100125 Mus musculus Traip gene Proteins 0.000 description 30
- 230000008859 change Effects 0.000 description 12
- 240000002853 Nelumbo nucifera Species 0.000 description 9
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 9
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 9
- 230000007704 transition Effects 0.000 description 7
- 238000010351 charge transfer process Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000004899 motility Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 238000004940 physical analysis method Methods 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000005309 stochastic process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Artificial Intelligence (AREA)
- Mathematical Analysis (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Primary Health Care (AREA)
- Mathematical Optimization (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Computational Mathematics (AREA)
- Quality & Reliability (AREA)
- Probability & Statistics with Applications (AREA)
Abstract
The invention discloses an electric vehicle charge load prediction method and system based on Markov chain; the method comprises the steps of setting original data information and a number of simulations; determining decision-making behaviors; accumulating charge load curves to obtain total charge load; judging whether the number of currently performed simulations is greater than a set number of simulations or not; judging whether the total charge load acting as simulation results is converged or not; if yes, outputting the simulation results; if not, returning to reset the number of simulations; based on Markov theory, stage of charge value of a power battery of an electric vehicle is considered as a stage variable in Markov chain, the method provides a visual indication for changes in the state of charge value of the power battery during one-day outgoing of a user so as to determine real-time charging need in this process, charge load is calculated through Monte Carlo simulation on such basis, and therefore load prediction results obtained by using the method are more practical.
Description
Technical field
The present invention relates to charging electric vehicle load prediction field is and in particular to a kind of be based on markovian electronic vapour
Car charging load forecasting method.
Background technology
Increasingly serious with energy crisis and environmental problem, electric automobile has become as the future automobile energy and environmental protection
Development trend." the accelerating the instruction of new-energy automobile popularization and application " put into effect according to State Council, to the year two thousand twenty, the new energy of China
Source car production and marketing is up to 5,000,000, and charging pile quantity is up to 4,500,000 it is seen that electric automobile quantity will have in the coming years
There is suitable scale.Due to the aggregation properties of charging electric vehicle, and the uncertainty of automobile user demand and behavior
With mutual difference, large-scale electric automobile accesses, and will bring huge challenge to the safe operation of electrical network and Optimized Operation, right
The operation of the planning of power system, operation and electricity market produces profound influence, therefore, sets up effective charging electric vehicle
Load forecasting model, the impact bringing for analysis charging electric vehicle load and electric automobile extensively access the regulating strategy of electrical network
Theoretical foundation is provided.
Do numerous studies for charging electric vehicle load modeling both at home and abroad, be broadly divided into three classes: based on electric automobile
Using the Monte Carlo Analogue Method of trip requirements, analyze the general of the charging electric vehicle power arriving at charging station using queuing theory
Rate analytic process and physical analysis method.Wherein, the statistics modeling method application based on Monte Carlo simulation is more universal, generally examines
Consider charging electric vehicle workload demand be subject to charging electric vehicle pattern, scale, moving law, battery behavior and electricity price regulation etc. because
The impact of element, sets up the load model based on probability statistics.Research great majority are in given starting impetus battery charge state
Set up simple model under conditions of (state of charge, soc), initiation of charge time and daily travel, ignore use
The motility of family charging behavior and randomness, load prediction subjective, not from the feature of electric automobile itself,
Real time charging demand during reflection user's trip, predicts the outcome and can not fully reflect the workload demand of reality.
Content of the invention
The technical problem to be solved is to provide one kind pre- based on markovian charging electric vehicle load
Survey method, the method is based on Markov theory, and electric automobile power battery SOC is considered as in Markov Chain
State variable, calculates charging load, therefore, the load prediction that obtains using the method by Monte Carlo simulation on this basis
Result more conforms to practical situation.
The technical scheme is that
One kind is based on markovian charging electric vehicle load forecasting method it is characterised in that including step:
(1) setting primary data information (pdi) and number realization;
(2) battery charge according to primary data information (pdi) extract unique user one day travel time and trip moment first
State value;
(3) extract single distance travelled, and calculate the battery lotus at the end of this section of stroke end time and this section of stroke
Electricity condition value;
(4) judge whether the battery charge state value at the end of described stroke is less than the threshold value setting, if less than setting
Threshold value, then execution step (5), otherwise, then execution step (7),
(5) randomly draw user's single charge duration, at the end of judging charging modes that user may take and charging
Battery charge state value and end time;Record the battery lotus at the end of this charging duration, charging modes, and charging simultaneously
Electricity condition value and end time;
(6) according to the battery charge state value at the end of this decision behavior and end time, judge that next step may be taken
Decision behavior, if next decision behavior be stop, execution step (7);If next decision behavior is to travel, execute
Step (9);
(7) randomly draw user's single parking duration, record battery charge state value and the knot at the end of this parking behavior
The bundle time;
(8) described end time and the trip end time on the same day setting are compared, if the end time is more than or equal to setting
The trip end time, then execution step (11), otherwise execution step (9);
(9) extract single distance travelled, record the battery charge state value at the end of this traveling behavior and end time;
(10) above-mentioned end time and the trip end time on the same day setting are compared, if the end time is more than or equal to setting
The trip end time, then execution step (11), otherwise execution step (4);
(11) respectively add up fill soon, the charging load curve of trickle charge, obtain total bulk charging load;
(12) whether the number of times judging at present simulation is more than the number realization setting in step (1);If so, then enter step
Suddenly (13);If it is not, then return to step (2) continues executing with;
(13) judge whether restrain as total bulk charging load of analog result;If convergence, export analog result;If no
Convergence, then return to step (1) resets number realization.
The invention has the beneficial effects as follows this method is based on Markov theory, by electric automobile power battery SOC
It is considered as the state variable in Markov Chain, automobile user trip process power battery charged shape under different decision behaviors
The change procedure of state value is described by markovian state conversion process, and the method visually reflects user one day
The state change of power battery charged state value during trip, and then determine the real time charging demand of this process, here basis
Above charging load is calculated by Monte Carlo simulation, therefore, more conform to reality using the load prediction results that the method obtains
Situation.
On the basis of technique scheme, the present invention can also do following improvement.
Further, the primary data information (pdi) in described step (1) includes:
Fill soon, the charge power of trickle charge, the battery capacity of electric automobile, the average overall travel speed of vehicle, every 100 kms consumption
Electricity;
Initial time and initial SOC that electric automobile is gone on a journey in one day first;Travel every time and terminate rear battery lotus
The threshold value of electricity condition;
The probability density function of user's single operating range;
The probability density function of the single charge duration of electric automobile;
The probability density function of single parking duration;
One day trip end time of user.
Further, the extraction single distance travelled in described step (3), and calculate this section of stroke end time and this section
Battery charge state value at the end of stroke, the calculating of described battery charge state value is drawn using equation below:
Wherein, dnFor user's single operating range, w100For the every 100 km power consumption of electric automobile, qkFor the electronic vapour of kth
The battery capacity of car, soci- for dnw1 electricity 00 electrical automobile any time battery charge state value, socj: for electric automobile at certain
Plant the battery charge state value of the subsequent time under decision behavior.
Further, judge that the charging modes detailed process that user takes is in described step (5),
Charging duration scope is soon:
Trickle charge charging duration scope is:
It is assumed that pk=5pm, then
It is a length of when user chargesThen take fast charge mode,
It is a length of when user chargesThen take trickle charge mode,
Wherein, pk、pmBe respectively fill soon, the charge power of trickle charge, tiFor charging duration, qkElectricity for kth electric automobile
Tankage.
Further, the battery charge state value at the end of the charging in described step (5) draws according to formula below:
Wherein, ai=1++ represents and fills soon, ai=1+ represents trickle charge, pk、pmBe respectively fill soon, the charge power of trickle charge, tiFor
Charging duration, qkFor the battery capacity of kth electric automobile, sociFor the battery charge state value of electric automobile any time,
socj: for the battery charge state value of electric automobile subsequent time at any time.
Further, in described step (11) respectively add up fill soon, the charging load curve of trickle charge, obtain total bulk charging bear
Lotus, concrete calculating process is as follows:
A length of during charging soon:
Trickle charge is a length of when charging:
Then the charge power that fills soon of arbitrarily t is:
The trickle charge charge power of arbitrarily t is:
Wherein, ti,1、ti,2Represent respectively fill soon, trickle charge charging duration, soci,0Proceed by the lotus of charging for electric automobile
Electricity condition, soci,1For the state-of-charge at the end of charging, t0For starting to charge up the moment,
Charging electric vehicle power is added up, can predict fill soon, trickle charge load and total load, computing formula divides
Not as follows:
Wherein, pt,1、pt,2Respectively represent t fill soon, trickle charge load, ptFor the total charging load of t electric automobile,
N1, n2 are natural number, represent the vehicle number filling soon with trickle charge respectively.
According to another aspect of the present invention, one kind is also provided to be based on markovian charging electric vehicle load prediction
System is it is characterised in that include information setup module: for inputting primary data information (pdi), and sets number realization;
First abstraction module: for extracting unique user travel time and trip first in a day from primary data information (pdi)
The battery charge state value in moment;
Second abstraction module: for extracting single distance travelled, and calculate this section of stroke end time and this section of stroke
At the end of battery charge state value;
First comparison module: for the battery charge state value at the end of relatively described stroke whether less than the threshold setting
Value;
3rd abstraction module: be used for randomly drawing user's single charge duration, and judge the charging side that user may take
Battery charge state value at the end of formula and charging and end time;Record this charging duration, charging modes, Yi Jichong simultaneously
Battery charge state value at the end of electricity and end time;
Decision-making module: for according to the battery charge state value at the end of this decision behavior and end time, judging next
Walking the decision behavior taken is to stop or travel;
4th abstraction module: be used for randomly drawing user's single parking duration, record the battery at the end of this parking behavior
SOC and time;
Second comparison module: for the trip end time on the same day of relatively described parking behavior end time and setting;
4th abstraction module: for extracting single distance travelled, record the battery charge state at the end of this traveling behavior
Value and end time;
3rd comparison module: for the trip end time on the same day of relatively described traveling behavior end time and setting;
Accumulator module: for adding up respectively to fill soon, the charging load curve of trickle charge, obtain total bulk charging load;
4th comparison module: whether the number of times for relatively simulation at present is more than the number realization setting;
Judge module: for judging whether restrain as total bulk charging load of analog result.
Brief description
Fig. 1 is electric automobile state-of-charge transfer process schematic diagram of the present invention;
Fig. 2 is based on markovian load prediction flow chart for the present invention.
Specific embodiment
Below in conjunction with accompanying drawing, the principle of the present invention and feature are described, example is served only for explaining the present invention, and
Non- for limiting the scope of the present invention.
It is illustrated in figure 1 electric automobile state-of-charge transfer process schematic diagram of the present invention, dynamic using Markov Chain description
The change procedure of power battery charge state soc:
Markov, with the naming of Russia mathematician markov, markoff process be have markov property with
Machine process, in the case of given current knowledge or information, the past (i.e. current former historic state) is for prediction in the future (i.e.
Current later to-be) it is unrelated.Can be applicable to the research neck such as the hydrology, meteorology, earthquake, economic forecasting, management decision-making
Domain.
Markov Chain (markov chain) is the discrete time stochastic process in mathematics with Markov property.It
The model of one class stochastic systems is described, it is random for referring to system in each time state in which, from current time to
The state of one time presses certain probability transfer, and to-be is only relevant with status praesenss and its transition probability, and with former
State is unrelated, i.e. markov property.
According to Markov theory, if the soc value of electric automobile any time is considered as state si, under electric automobile for the moment
The state-of-charge soc value carved is considered as state sjIt is known that sjState s only with current timeiAnd corresponding transition probability pijIt is relevant,
So the state-of-charge change during electric automobile trip in a day can be described by Markov Chain.Due to various types of
In type electric automobile, the proportion shared by private car is larger, and the charging behavior of private car user has very big motility and not true
Qualitative, during trip, the change procedure of electrokinetic cell soc is complex.Therefore, this patent is mainly used using Markov Chain
To analyze in private car one day soc value changes situation thus predicting charging load.
In figure aiDifferent values represent different decision behaviors, wherein, ai=-1 expression electric automobile during traveling;ai=1+
+、ai=1+ represents that electric automobile is filled and trickle charge soon respectively;ai=0 expression is not charged and is not also travelled.Solid arrows represent
State transfer case under certain decision behavior, heavy solid line arrows represent under current state, the possible decision behavior of next step, empty
If line arrow represents is currently not one day end-state snWhen, the possible decision behavior of next step.pijSoc value for electric automobile
Switch to the state transition probability of next state from preceding state, it has different representations as the case may be.
(1) decision behavior that user may take.State-of-charge transition probability p due to electric automobileijWith current time
State-of-charge and current time decision behavior a that electric automobile is taken within the subsequent time time periodiRelevant.Therefore, root
During going on a journey according to electric automobile, different travel behaviours (charging, not charging also does not travel and travel three kinds) are to electrokinetic cell soc
The impact that change produces, the decision behavior that may take during user is gone on a journey is classified as follows.
(2) state transition probability.User's charging behavior is to meet demand, if charged terminate rear electrokinetic cell soc
Value there is no change, then this charges just nonsensical.After a period of time it is therefore desirable to user charges, its soc value should
There is certain change, in this patent, suppose that the excursion of its soc value is 0.2≤sock< 1, by filling soon, trickle charge power different,
Can release its correspondingly fill soon, the distribution of trickle charge charging duration, derivation is as follows:
Charging duration scope soon:
Trickle charge charging duration scope:
In formula: pk、pmIt is respectively the fast, charge power of trickle charge, tiFor charging duration, qkBattery for kth electric automobile
Capacity.It is assumed that the expectation of quick charge power is about 5 times of trickle charge, it is computed obtainingTherefore, basis in turn
The current chargeable duration of user may determine that the charging modes that user may take, that is, when the chargeable duration of userWhen it is believed that private car takes fast charge mode, charging durationWhen, take is trickle charge mode.
In sum, according to taking different decision behaviors it can be deduced that the transition probability of corresponding decision, remember pijIt is from shape
State siSwitch to state sjState transition probability, its concrete representation is as follows.
In formula:Be respectively user take fill soon, the probability of trickle charge, t (ti) it is single
Secondary chargeable duration tiProbability density function, f (dn) for user's single operating range probability-distribution function.
(3) soc of subsequent time electric automobile.The soc of known current time and from current time start a period of time in
The decision behavior taken, then draw the soc value of subsequent time, and calculating process is as follows.
Work as ai=1+ or aiDuring=1++, electric automobile is charged, thenWherein charge power p
Value is:
Work as aiWhen=0, electric automobile is not charged also not travelling, then its state-of-charge does not change, i.e. socj=
soci
Work as aiWhen=- 1, electric automobile during traveling, then
In formula: pk、pmIt is respectively the fast, charge power of trickle charge, dnFor user's single operating range, w100Every for electric automobile
100km power consumption, qkBattery capacity for kth electric automobile.
To sum up can obtain the soc of subsequent time electric automobile:
It is illustrated in figure 2 the present invention and is based on markovian load prediction flow chart, the present invention provides according to existing document
Material and DOT issue American family trip survey data (national household travel survey,
Nhts) give following condition:
(1) fill soon, the charge power of trickle charge, the battery capacity of electric automobile, the average overall travel speed of vehicle, every kilometer consumption
Electricity;(2) electric automobile is gone on a journey in one day first initial time and initial state-of-charge soc value;(3) after traveling terminates every time
When soc drops to certain threshold values soc, < when 0.2, electric automobile can select to be charged in charging station;(4) electric automobile list known to
The probability density function of secondary distance travelled;(5) probability-distribution function of the single charge duration of electric automobile;(6) single stops
The probability density function of persistent period;(7) one day trip finish time t of userend.
The present invention describes the situation of change of electrokinetic cell soc during user's trip using Markov Chain, by electronic vapour
Car electrokinetic cell soc, as state variable, goes on a journey according to automobile user during being accustomed to determining unique user trip in a day
The decision behavior of day part, simulation produces user's trip sequence of a day, so that it is determined that all users charge period of a day, filling
Electrically, adopt monte carlo method that the trip sequence of user is simulated based on conditions above, process is as follows:
(1) setting primary data information (pdi) and number realization;Wherein primary data information (pdi) includes:
Fill soon, the charge power of trickle charge, the battery capacity of electric automobile, the average overall travel speed of vehicle, every 100 kms consumption
Electricity;
Initial time and initial SOC that electric automobile is gone on a journey in one day first;Travel every time and terminate rear battery lotus
The threshold value of electricity condition;
The probability density function of user's single operating range;
The probability density function of the single charge duration of electric automobile;
The probability density function of single parking duration;
One day trip end time of user.
(2) battery charge according to primary data information (pdi) extract unique user one day travel time and trip moment first
State value;
(3) extract single distance travelled, and calculate the battery lotus at the end of this section of stroke end time and this section of stroke
Electricity condition value, the calculating of wherein battery charge state value is drawn using equation below:
Wherein, dnFor user's single operating range, w100For the every 100 km power consumption of electric automobile, qkFor the electronic vapour of kth
The battery capacity of car, soci- for electric automobile any time battery charge state value, socj: for electric automobile in certain decision-making
The battery charge state value of the subsequent time under behavior.
(4) judge whether the battery charge state value at the end of described stroke is less than the threshold value setting, if less than setting
Threshold value, then execution step (5), otherwise, then execution step (7),
(5) randomly draw user's single charge duration, at the end of judging charging modes that user may take and charging
Battery charge state value and end time;Record the battery charge shape at the end of this charging duration, charging modes, and charging
State value and end time;
Wherein, judge that the charging modes that user may take are determined by formula below:
Charging duration scope is soon:
Trickle charge charging duration scope is:
It is assumed that pk=5pm, then
It is a length of when user chargesThen take fast charge mode,
It is a length of when user chargesThen take trickle charge mode,
Wherein, pk、pmBe respectively fill soon, the charge power of trickle charge, tiFor charging duration, qkElectricity for kth electric automobile
Tankage;
Battery charge state value at the end of charging draws according to formula below:
Wherein, ai=1++ represents and fills soon, ai=1+ represents trickle charge, pk、pmBe respectively fill soon, the charge power of trickle charge, tiFor
Charging duration, qkFor the battery capacity of kth electric automobile, soci- for electric automobile any time battery charge state value,
socj: for the battery charge state value of subsequent time under certain decision behavior for the electric automobile.
(6) according to the battery charge state value at the end of this decision behavior and end time, judge that next step may be taken
Decision behavior, if next decision behavior be stop, execution step (7);If next decision behavior is to travel, execute
Step (9);
(7) randomly draw user's single parking duration, record battery charge state value and the knot at the end of this parking behavior
The bundle time;
(8) described end time and the trip end time on the same day setting are compared, if the end time is more than or equal to setting
The trip end time, then execution step (11), otherwise execution step (9);
(9) extract single distance travelled, record the battery charge state value at the end of this traveling behavior and end time;
(10) above-mentioned end time and the trip end time on the same day setting are compared, if the end time is more than or equal to setting
The trip end time, then execution step (11), otherwise execution step (4);
(11) respectively add up fill soon, the charging load curve of trickle charge, obtain total bulk charging load;Concrete calculating process is such as
Under:
A length of during charging soon:
Trickle charge is a length of when charging:
Then the charge power that fills soon of arbitrarily t is:
The trickle charge charge power of arbitrarily t is:
Wherein, ti,1、ti,2Represent respectively fill soon, trickle charge charging duration, soci,0Proceed by the lotus of charging for electric automobile
Electricity condition, soci,1For the state-of-charge at the end of charging, t0For starting to charge up the moment,
Charging electric vehicle power is added up, can predict fill soon, trickle charge load and total load, computing formula divides
Not as follows:
Wherein, pt,1、pt,2Represent the filling soon of t, trickle charge load, p respectivelytBear for the total charging of t electric automobile
Lotus, n1, n2 are natural number, represent the vehicle number filling soon with trickle charge respectively.
(12) whether the number of times judging at present simulation is more than the number realization setting in step (1);If so, then enter step
Suddenly (13);If it is not, then return to step (2) continues executing with;
(13) judge whether restrain as total bulk charging load of analog result;If convergence, export analog result;If no
Convergence, then return to step (1) resets number realization.
The present invention is based on Markov theory, by automobile user trip process electrokinetic cell under different decision behaviors
The change procedure of soc is described by Markov Chain, obtains the situation of change of electrokinetic cell soc during this, and simulation is used
The one day trip sequence in family, determines charging electric vehicle period, charging modes, calculates charging load by Monte Carlo simulation, should
Method has vividly described the real time charging demand of user's trip process, is more accorded with based on the load prediction results that the method obtains
Close practical situation.
The foregoing is only presently preferred embodiments of the present invention, not in order to limit the present invention, all spirit in the present invention and
Within principle, any modification, equivalent substitution and improvement made etc., should be included within the scope of the present invention.
Claims (7)
1. one kind is based on markovian charging electric vehicle load forecasting method it is characterised in that including step:
(1) setting primary data information (pdi) and number realization;
(2) battery charge state in extract unique user one day from primary data information (pdi) travel time and trip moment first
Value;
(3) extract single distance travelled, and calculate the battery charge shape at the end of this section of stroke end time and this section of stroke
State value;
(4) whether the battery charge state value relatively at the end of described stroke is less than the threshold value setting, if less than the threshold setting
Value, then execution step (5), otherwise, then execution step (7);
(5) randomly draw user's single charge duration, the battery charge at the end of judging charging modes that user takes and charging
State value and end time;Record the battery charge state value at the end of this charging duration, charging modes, and charging simultaneously
And the end time;
(6) according to the battery charge state value at the end of this decision behavior and end time, judge the decision-making row that next step is taken
For, if next decision behavior is to stop, execution step (7);If next decision behavior is to travel, execution step (9);
(7) randomly draw user's single parking duration, record battery charge state value and the time at the end of this parking behavior;
(8) relatively described parking behavior end time and the same day setting go on a journey the end time, if the end time is more than or equal to setting
The fixed trip end time, then execution step (11), otherwise execution step (9);
(9) extract single distance travelled, record the battery charge state value at the end of this traveling behavior and end time;
(10) relatively described traveling behavior end time and the same day setting go on a journey the end time, if the end time is more than or equal to setting
The fixed trip end time, then execution step (11), otherwise execution step (4);
(11) respectively add up fill soon, the charging load curve of trickle charge, obtain total bulk charging load;
(12) whether the number of times of relatively simulation at present is more than the number realization setting in step (1);If so, then enter step
(13);If it is not, then return to step (2) continues executing with;
(13) judge whether restrain as total bulk charging load of analog result;If convergence, export analog result;If not receiving
Hold back, then return to step (1) resets number realization.
2. according to claim 1 based on markovian charging electric vehicle load forecasting method it is characterised in that
Primary data information (pdi) in described step (1) includes:
Fill soon, the charge power of trickle charge, the battery capacity of electric automobile, the average overall travel speed of vehicle, every 100 km power consumptions
Amount;
Initial time and initial SOC that electric automobile is gone on a journey in one day first;Travel every time and terminate rear battery charge shape
The threshold value of state;
The probability density function of user's single operating range;
The probability density function of the single charge duration of electric automobile;
The probability density function of single parking duration;
One day trip end time of user.
3. according to claim 1 based on markovian charging electric vehicle load forecasting method it is characterised in that
Extraction single distance travelled in described step (3), and calculate the electricity at the end of this section of stroke end time and this section of stroke
Pond SOC, the calculating of described battery charge state value is drawn using equation below:
Wherein, dnFor user's single operating range, w100For the every 100 km power consumption of electric automobile, qkFor kth electric automobile
Battery capacity, sociFor the battery charge state value of any current time of electric automobile, socjFor electric automobile in certain decision-making row
For under subsequent time battery charge state value.
4. according to claim 1 based on markovian charging electric vehicle load forecasting method it is characterised in that
Judge that the charging modes detailed process that user takes is in described step (5): charging duration scope is soon:
Trickle charge charging duration scope is:
It is assumed that pk=5pm, then
It is a length of when user chargesThen take fast charge mode,
It is a length of when user chargesThen take trickle charge mode,
Wherein, pk、pmBe respectively fill soon, the charge power of trickle charge, tiFor charging duration, qkBattery for kth electric automobile holds
Amount.
5. according to claim 1 based on markovian charging electric vehicle load forecasting method it is characterised in that
The battery charge state value at the end of charging in described step (5) draws according to formula below:
Wherein, ai=1++ represents and fills soon, ai=1+ represents trickle charge, pk、pmBe respectively fill soon, the charge power of trickle charge, tiFor charging
Duration, qkFor the battery capacity of kth electric automobile, dnFor user's single operating range, w100For the every 100 km consumptions of electric automobile
Electricity, sociFor the battery charge state value of electric automobile any time, socjFor electric automobile under certain decision behavior under
The battery charge state value in one moment.
6. according to claim 1 based on markovian charging electric vehicle load forecasting method it is characterised in that
In described step (11) respectively add up fill soon, the charging load curve of trickle charge, obtain total bulk charging load, concrete calculating process is such as
Under:
A length of during charging soon:
Trickle charge is a length of when charging:
Then the charge power that fills soon of arbitrarily t is:
The trickle charge charge power of arbitrarily t is:
Wherein, ti,1、ti,2Represent respectively fill soon, trickle charge charging duration, soci,0Proceed by the charged shape of charging for electric automobile
State, soci,1For the state-of-charge at the end of charging, t0For starting to charge up the moment,
Charging electric vehicle power is added up, predict fill soon, trickle charge load and total load, computing formula is as follows respectively:
Wherein, pt,1、pt,2Represent respectively fill soon, trickle charge load, ptFor the total charging load of electric automobile, n1, n2 are natural number,
Represent the vehicle number filling soon with trickle charge respectively.
7. one kind is based on markovian charging electric vehicle load prediction system it is characterised in that including information setting mould
Block: for inputting primary data information (pdi), and set number realization;
First abstraction module: for extracting unique user one day travel time and the trip moment first from primary data information (pdi)
Battery charge state value;
Second abstraction module: for extracting single distance travelled, and calculate this section of stroke end time and this section of stroke terminates
When battery charge state value;
First comparison module: for the battery charge state value at the end of relatively described stroke whether less than the threshold value setting;
3rd abstraction module: be used for randomly drawing user's single charge duration, and judge charging modes that user may take and
Battery charge state value at the end of charging and end time;Record this charging duration, charging modes simultaneously, and the knot that charges
Battery charge state value during bundle and end time;
Decision-making module: for according to the battery charge state value at the end of this decision behavior and end time, judging that next step is adopted
The decision behavior taking is to stop or travel;
4th abstraction module: be used for randomly drawing user's single parking duration, record the battery charge at the end of this parking behavior
State value and time;
Second comparison module: for the trip end time on the same day of relatively described parking behavior end time and setting;
4th abstraction module: for extracting single distance travelled, record battery charge state value at the end of this traveling behavior and
End time;
3rd comparison module: for the trip end time on the same day of relatively described traveling behavior end time and setting;
Accumulator module: for adding up respectively to fill soon, the charging load curve of trickle charge, obtain total bulk charging load;
4th comparison module: whether the number of times for relatively simulation at present is more than the number realization setting;
Judge module: for judging whether restrain as total bulk charging load of analog result.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610836365.7A CN106355290A (en) | 2016-09-21 | 2016-09-21 | Electric vehicle charge load prediction method and system based on Markov chain |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610836365.7A CN106355290A (en) | 2016-09-21 | 2016-09-21 | Electric vehicle charge load prediction method and system based on Markov chain |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106355290A true CN106355290A (en) | 2017-01-25 |
Family
ID=57858180
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610836365.7A Pending CN106355290A (en) | 2016-09-21 | 2016-09-21 | Electric vehicle charge load prediction method and system based on Markov chain |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106355290A (en) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107169588A (en) * | 2017-04-12 | 2017-09-15 | 中国电力科学研究院 | A kind of electric automobile charging station short-time rating Forecasting Methodology and system |
CN107392400A (en) * | 2017-09-04 | 2017-11-24 | 重庆大学 | Meter and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodology |
CN108133329A (en) * | 2017-12-29 | 2018-06-08 | 天津大学 | Consider the electric vehicle trip of charging feedback effect and charge requirement analysis method |
CN108288112A (en) * | 2018-01-30 | 2018-07-17 | 广西电网有限责任公司柳州供电局 | Region electric automobile charging station load forecasting method based on user's trip simulation |
CN109117972A (en) * | 2017-06-22 | 2019-01-01 | 南京理工大学 | A kind of charge requirement of electric car determines method |
CN109754136A (en) * | 2017-11-03 | 2019-05-14 | 蔚来汽车有限公司 | Battery equalization method and system |
CN109774492A (en) * | 2018-12-29 | 2019-05-21 | 江苏大学 | A kind of whole pure electric vehicle power distribution method based on the following driving power demand |
CN110059937A (en) * | 2019-03-29 | 2019-07-26 | 四川大学 | A kind of load modeling method of integration electric automobile full track trace space |
CN111815017A (en) * | 2020-05-29 | 2020-10-23 | 国网山东省电力公司经济技术研究院 | Electric vehicle charging load prediction method based on travel data |
CN112330025A (en) * | 2020-11-06 | 2021-02-05 | 国网冀北电力有限公司张家口供电公司 | Prediction method of space-time charging load for urban electric vehicle |
CN112671022A (en) * | 2020-12-31 | 2021-04-16 | 国网山东省电力公司泰安供电公司 | Optical storage charging station capacity optimal configuration method, system, terminal and storage medium |
CN113657456A (en) * | 2021-07-26 | 2021-11-16 | 华南理工大学 | Ordered charging method for allocating user charging based on Markov chain |
CN113887811A (en) * | 2021-10-13 | 2022-01-04 | 江苏明茂新能源科技有限公司 | Charging pile data management method and system |
CN118134228A (en) * | 2024-05-10 | 2024-06-04 | 浙江大学 | Online dispatching method for industrial park comprehensive energy system containing electric automobile load |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413180A (en) * | 2013-07-22 | 2013-11-27 | 上海电力实业有限公司 | Electric car charging load forecasting system and method based on Monte Carlo simulation method |
CN103997091A (en) * | 2014-05-23 | 2014-08-20 | 国家电网公司 | Scale electric automobile intelligent charging control method |
CN104734171A (en) * | 2015-04-16 | 2015-06-24 | 合肥工业大学 | Electric vehicle charging station modeling method for reliability assessment of power distribution network and application of electric vehicle charging station modeling method |
-
2016
- 2016-09-21 CN CN201610836365.7A patent/CN106355290A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103413180A (en) * | 2013-07-22 | 2013-11-27 | 上海电力实业有限公司 | Electric car charging load forecasting system and method based on Monte Carlo simulation method |
CN103997091A (en) * | 2014-05-23 | 2014-08-20 | 国家电网公司 | Scale electric automobile intelligent charging control method |
CN104734171A (en) * | 2015-04-16 | 2015-06-24 | 合肥工业大学 | Electric vehicle charging station modeling method for reliability assessment of power distribution network and application of electric vehicle charging station modeling method |
Non-Patent Citations (1)
Title |
---|
XU WEI等: "Optimal Allocation for Charging Piles in Multi-areas Considering Charging Load Forecasting Based on Markov Chain", 《2016 CHINA INTERNATIONAL CONFERENCE ON ELECTRICITY DISTRIBUTION》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107169588A (en) * | 2017-04-12 | 2017-09-15 | 中国电力科学研究院 | A kind of electric automobile charging station short-time rating Forecasting Methodology and system |
CN109117972A (en) * | 2017-06-22 | 2019-01-01 | 南京理工大学 | A kind of charge requirement of electric car determines method |
CN107392400A (en) * | 2017-09-04 | 2017-11-24 | 重庆大学 | Meter and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodology |
CN107392400B (en) * | 2017-09-04 | 2020-10-16 | 重庆大学 | EV charging load space-time distribution prediction method considering real-time traffic and temperature |
CN109754136A (en) * | 2017-11-03 | 2019-05-14 | 蔚来汽车有限公司 | Battery equalization method and system |
CN108133329A (en) * | 2017-12-29 | 2018-06-08 | 天津大学 | Consider the electric vehicle trip of charging feedback effect and charge requirement analysis method |
CN108133329B (en) * | 2017-12-29 | 2021-06-08 | 天津大学 | Electric automobile travel and charging demand analysis method considering charging feedback effect |
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 |
CN109774492A (en) * | 2018-12-29 | 2019-05-21 | 江苏大学 | A kind of whole pure electric vehicle power distribution method based on the following driving power demand |
CN110059937A (en) * | 2019-03-29 | 2019-07-26 | 四川大学 | A kind of load modeling method of integration electric automobile full track trace space |
CN111815017A (en) * | 2020-05-29 | 2020-10-23 | 国网山东省电力公司经济技术研究院 | Electric vehicle charging load prediction method based on travel data |
CN112330025A (en) * | 2020-11-06 | 2021-02-05 | 国网冀北电力有限公司张家口供电公司 | Prediction method of space-time charging load for urban electric vehicle |
CN112671022A (en) * | 2020-12-31 | 2021-04-16 | 国网山东省电力公司泰安供电公司 | Optical storage charging station capacity optimal configuration method, system, terminal and storage medium |
CN112671022B (en) * | 2020-12-31 | 2024-05-31 | 国网山东省电力公司泰安供电公司 | Optical storage charging station capacity optimal configuration method, system, terminal and storage medium |
CN113657456A (en) * | 2021-07-26 | 2021-11-16 | 华南理工大学 | Ordered charging method for allocating user charging based on Markov chain |
CN113657456B (en) * | 2021-07-26 | 2023-08-25 | 华南理工大学 | Ordered charging method for allocating user charging based on Markov chain |
CN113887811A (en) * | 2021-10-13 | 2022-01-04 | 江苏明茂新能源科技有限公司 | Charging pile data management method and system |
CN118134228A (en) * | 2024-05-10 | 2024-06-04 | 浙江大学 | Online dispatching method for industrial park comprehensive energy system containing electric automobile load |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106355290A (en) | Electric vehicle charge load prediction method and system based on Markov chain | |
CN105160451B (en) | A kind of micro-capacitance sensor Multiobjective Optimal Operation method containing electric vehicle | |
Zhao et al. | A deep reinforcement learning framework for optimizing fuel economy of hybrid electric vehicles | |
CN107169273A (en) | The charging electric vehicle power forecasting method of meter and delay and V2G charge modes | |
CN103246943B (en) | Vehicle operating condition multi-scale predicting method based on Markov chain | |
CN109501630A (en) | A kind of electric car charging scheme real-time recommendation method and system | |
CN111626527B (en) | Intelligent power grid deep learning scheduling method considering fast/slow charging/discharging form of schedulable electric vehicle | |
CN103903090B (en) | Electric car charging load distribution method based on user will and out-going rule | |
CN103499792B (en) | The Forecasting Methodology of available capacity of EV power battery cluster | |
CN103065199A (en) | Electric vehicle charging station load forecasting method | |
Gan et al. | A probabilistic evaluation method of household EVs dispatching potential considering users’ multiple travel needs | |
CN105160428A (en) | Planning method of electric vehicle fast-charging station on expressway | |
CN104281889B (en) | Electric automobile EV charging load multiple objective random programming method | |
CN106203720A (en) | A kind of Multiple Time Scales electric automobile cluster schedulable capacity prediction methods | |
CN106295860A (en) | A kind of electric automobile scale charge requirement Forecasting Methodology based on Monte Carlo Analogue Method | |
CN110363332A (en) | A kind of electric car charging load spatial and temporal distributions prediction technique based on individual behavior characteristic | |
CN113103905A (en) | Intelligent charging distribution adjusting method, device, equipment and medium for electric automobile | |
CN105068419B (en) | Residential area electric automobile charge/discharge control method | |
CN106786977B (en) | Charging scheduling method of electric vehicle charging station | |
CN109308541A (en) | Electric vehicle charging load prediction method based on scene analysis | |
CN105226650A (en) | Based on the micro-capacitance sensor reliability calculation method of miniature combustion engine-energy storage cooperation strategy | |
CN104102954B (en) | Optimal configuration method of distributed comprehensive energy supply system considering black start function | |
CN108876052A (en) | Electric car charging load forecasting method, device and computer equipment | |
CN106407726A (en) | Method for selecting electrical access point of electric automobile charging station by considering influence on tidal flow | |
CN107067130A (en) | A kind of quick charge station method for planning capacity based on electric automobile markov charge requirement analysis model |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170125 |