CN113902182A - Robust charging optimization method for electric bus fleet considering energy consumption uncertainty - Google Patents
Robust charging optimization method for electric bus fleet considering energy consumption uncertainty Download PDFInfo
- Publication number
- CN113902182A CN113902182A CN202111150148.XA CN202111150148A CN113902182A CN 113902182 A CN113902182 A CN 113902182A CN 202111150148 A CN202111150148 A CN 202111150148A CN 113902182 A CN113902182 A CN 113902182A
- Authority
- CN
- China
- Prior art keywords
- charging
- energy consumption
- constraint
- robust
- power
- 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
- 238000005265 energy consumption Methods 0.000 title claims abstract description 94
- 238000005457 optimization Methods 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000005611 electricity Effects 0.000 claims abstract description 26
- 230000009977 dual effect Effects 0.000 claims description 28
- 230000014509 gene expression Effects 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 5
- 239000013589 supplement Substances 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000001502 supplementing effect Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 238000000844 transformation Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000007726 management method Methods 0.000 description 9
- 230000003111 delayed effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
Abstract
A robust charging optimization method for an electric bus fleet based on uncertain operation energy consumption is characterized by firstly utilizing a robust optimization technology to construct a budget uncertain set to depict uncertainty of operation energy consumption, ensuring robustness of an obtained charging optimization result and avoiding limitation that a traditional uncertain optimization method assumes probability distribution of uncertain parameter clothes in advance. Secondly, a robust charging optimization model of the electric bus fleet is built, the time-of-use electricity price, charging station resource constraints, uncertain energy consumption requirements of vehicles, the influence of multiple constraints such as operation time on the charging scheme of the electric bus fleet is fully considered, and the actual operation condition is fitted. And finally, solving according to a frame design column generation algorithm of the robust model, decomposing the robust optimization model into a main problem and a sub problem for iterative solution, and accelerating the solving efficiency. The invention provides an economic and reliable charging scheme for the charging management of the electric bus fleet, obviously reduces the charging cost, effectively copes with energy consumption fluctuation, and improves the economy and stability of the electric bus charging system.
Description
Technical Field
The invention belongs to the technical field of urban electric public system management and uncertain problem optimization, and relates to a robust optimization method for orderly charging of an electric bus fleet facing to energy consumption fluctuation influence, in particular to a robust charging optimization method for the electric bus fleet considering energy consumption uncertainty.
Background
The electric public transport system is used as the backbone of green traffic, and provides a good solution for solving the problems of environmental pollution, energy shortage and the like. Compared with the traditional fuel vehicle, the electric vehicle has the difficulties of mileage limitation, battery technology and the like, so that the charging management of the electric vehicle becomes a key link for ensuring the safe and stable operation of a bus fleet. However, the operation power consumption of the electric bus is greatly uncertain due to the influence of weather, road environment, traffic conditions, passenger capacity and other factors, and the power consumption of the same shift may have a deviation of about 5% -20% or more, thereby making the actual charging operation very challenging.
Through literature and patent search, it can be found that in current research results, a charging optimization strategy of an electric bus generally performs charging management based on a determined power consumption parameter, and ignores the uncertainty of energy consumption in a driving process. Wang et al, in the literature [ Wang Y, Huang Y, Xu J, et al, optimal retrieval scheduling for urban electric buses [ J ]. Transportation Research Part E logics & Transportation view,2017,100: 115. 132 ], propose a modeling framework for optimizing electric bus charging plans, which determines charging plans and operation decisions, while minimizing the total annual cost. The method is characterized in that a discrete binary particle swarm algorithm is adopted in a document (Wu Xiao Mei, Von Qijin, rigor, permissive garden, Pan Jing, electric bus ordered charging strategy [ J ] based on double-layer optimization, power grid and clean energy, 2021,37(01): 119-. He et al propose a network modeling framework in a document [ He Y, Liu Z, Song Z. optimal charging scheduling and management for a fast-charging bus system [ J ]. Transportation Research Part E: logics and Transportation Review,2020, 142 ] based on given power distribution network capacity, number of chargers, and bus operation information, to optimize the charging scheduling and management of the electric bus rapid charging system, effectively reducing the total charging cost including both power demand cost and energy cost. The charging optimization based on the determined power consumption parameters can simplify the problem, but the reliability of the thus-formulated charging strategy for the fleet is poor. This is because the power consumption of the electric bus directly determines the required charging duration. Ignoring the uncertainty of the electric bus energy consumption can therefore result in suboptimal or even infeasible solutions for the public transportation charging system. For example, when the power consumption of the electric bus increases during a shift, the electric bus is executed according to the original charging strategy, and various situations may occur, such as delayed departure caused by the fact that the battery power is lower than the safety threshold of the battery in the operation process or the battery is insufficient to supplement the power before the bus leaves, and the reliability and the safety of the public transportation system are damaged.
Disclosure of Invention
Based on the problem of neglecting the uncertainty of energy consumption in the charging optimization of the electric bus fleet, the invention provides the robust charging optimization method of the electric bus fleet in consideration of the uncertainty of the energy consumption. Firstly, a budget uncertain set is constructed by utilizing a robust optimization technology to depict uncertainty of operation energy consumption, on one hand, the robustness of the obtained charging optimization result is ensured, namely, the obtained solutions are feasible as long as parameters are disturbed in the given uncertain set, and on the other hand, the limitation that the traditional uncertain optimization method assumes the probability distribution of uncertain parameter service in advance is avoided. And secondly, establishing a robust optimization model for charging of the electric bus fleet under uncertain operation energy consumption, wherein time-of-use electricity price, charging station resource limitation and fluctuation of operation energy consumption are fully considered, so that the final optimization result meets the uncertain energy consumption for charging of the electric bus fleet and the charging cost is minimum. And finally, solving according to a frame design column generation algorithm of the robust model, decomposing the robust optimization model into a main problem and a sub problem for iterative solution, and accelerating the solving efficiency.
The technical scheme of the invention is as follows:
the robust charging optimization method of the electric bus fleet considering the uncertainty of energy consumption comprises the following steps:
step (1) of constructing an uncertainty set of operation energy consumption
And in consideration of the fluctuation of the energy consumption of the vehicle m in the operation process, a budget uncertain set is introduced to deal with the uncertainty of the operation energy consumption. Uncertain energy consumption of each vehicle m required to execute shift kIn a single interval, i.e.WhereinIs a nominal value of the uncertain energy consumption,the maximum deviation value representing the uncertain energy consumption is generally assumed to be p times the nominal value, i.e.The set of pre-computed uncertainties is defined as shown in equation (1). Wherein Z is shown as formula (2).
Constructed uncertainty set U of operating energy consumptionbugetCan be defined as:
wherein rmkThe energy consumption uncertainty budget for executing shift k is a parameter reflecting the uncertainty level, and the optimality and the robustness of the solution can be weighed. It is noted that 0 ≦ rmkK ≦ k, i.e. rmkThe maximum value of (a) is the maximum number of shifts for which the energy consumption deviates from its nominal value at the same time. When the deviation of the energy consumption of the vehicle m from the value of k is maximum for all shifts, i.e. per rmkK, which means that the worst case occurs, the charging scheme is most conservative.
Step (2), establishing an electric bus fleet charging robust optimization model under uncertain operation energy consumption
And based on the defined uncertain energy consumption set, converting the electric bus charging determined optimization model into a robust optimization model for orderly charging of the electric bus fleet under uncertain operation energy consumption. The method can ensure the feasibility of the charging optimization result without knowing the specific distribution of uncertain energy consumption, and obtain a charging scheme with higher robustness and optimality.
(2.1) parameter definition
Collection
M is the set of electric buses served by the charging station, wherein M is {1,2, …, M }
T is a set of times divided in one day, T ═ 1,2, …,1440/Δ T }
KmClass set of electric public transport vehicles m, Km={1,2,…,AmIn which A ismIs the total number of shifts of the vehicle m
M shift K of vehicle (K belongs to K)m) The set of operating time instants of (c),wherein c ismkIs the starting time, s, of the shift k executed by the vehicle leaving the charging stationmkIs the moment when the shift k is finished and returns to the charging station
parameter(s)
ftThe price of electricity at time t (unit: yuan/kW. h)
Delta t is unit time, (unit: minute), the invention selects 5 minutes as one unit time
qmkThe remaining battery power of the vehicle m returning to the charging station after the shift k is executed (unit: kW. h)
qmaxBattery capacity (unit: kW. h)
qminMinimum threshold allowed by battery power (unit: kW. h)
ΔetPower consumption of electric bus per minute (unit: kW. h)
PmaxTotal power supplied by charging station (unit: kW)
Mu loss rate of power transmission process from charging station to charging gun
N total number of charging guns in charging station
paMaximum allowable charging power for battery of electric bus (unit: kW)
pbMaximum charging power allowed by charging gun (unit: kW)
pcUpper limit of charging power (unit: kW)
dmtIf the bus m supplements the electric quantity at the moment t, d mt1 is ═ 1; otherwise dmt=0
(2.2) decision variables
pmtRequired charging power at time t of m for a vehicle, an integer variable, (unit: kW)
UmtAuxiliary binary variable for charge continuity control, for determining whether a change in state of charge occurs from time t to t +1, Umt∈{0,1}
VmtAuxiliary binary variable for charge continuity control for determining whether or not a change occurs in the state of charge from time t to t-1, Vmt∈{0,1}
Uncertainty set | z of robust model constraint (19)miDual variable corresponding to | less than or equal to 1
Uncertainty set | z of robust model constraint (20)miDual variable corresponding to | less than or equal to 1
Uncertainty set | z of robust model constraint (21)mkDual variable corresponding to | less than or equal to 1
(2.3) establishing a determined charging optimization model taking the minimum total charging cost of the electric bus fleet as an objective function, wherein the model is as follows:
an objective function:
constraint conditions are as follows:
Umt∈{0,1} (13)
Vmt∈{0,1} (14)
0≤pmt≤pc (15)
in the formula, the objective function formula (3) is intended to minimize the total charging cost of the electric bus fleet, and is the sum of electricity prices of all vehicles at all times of consumed electricity. The constraint condition formula (4) indicates that the electric bus cannot be replenished with electricity during operation, so the charging power is 0. Ensuring that the charging power can not exceed the maximum power p of the charging gun through the constraint condition formula (5)bAnd the maximum power p acceptable for the batterya. The constraint condition formula (6) and the constraint condition formula (7) are constraint conditions of charging station resource limitation, and the constraint condition formula (6) requires that the total charging power sum of the electric buses which are served simultaneously is not more than the total power available for the charging station; the constraint condition formula (7) requires that the number of the electric buses charged simultaneously is at most the total number of the charging guns. The constraint expressions (8.1) to (8.3) and the constraint expressions (9.1) to (9.5) both represent charging continuity control in which the vehicle, once charged, reaches a predetermined amount of electricity or stops for a predetermined time, and is not intermittently charged. Note that the form of constraint conditions (9.1) - (9.5) is complicated because the charging period after the last shift is performed includes two parts, i.e., the charging period after the last shift is performedThe charging continuity of the two-part join requires special handling. The constraint condition formula (10) indicates that the battery electric quantity value after each electricity supplementing should not exceed the battery capacity; the constraint equation (11) ensures that the amount of power replenished is not below the minimum threshold after completion of the normal operation of the next shift. The constraint (12) means that the amount of electricity supplied and the amount of electricity consumed are identical, i.e. the electric bus must be fully charged before the next day of departure. ConstrainingConditional expressions (13) to (15) give the types of variables.
(2.4) converting the determined charging optimization model into a charging robust optimization model under the condition of uncertain operation energy consumption
In order to convert the deterministic charging optimization model into a robust model that takes into account the uncertain energy consumption, it is necessary to replace the constraints (10) - (12) with the following constraints (16) - (18), when the model contains uncertain energy consumption parameters.
The three constraints of the constraint equations (16) - (18) are difficult to handle, and the min/max can be converted into the equations (19) - (21).
(2.5) conversion of robust optimization model
To further linearize the constraints (19) - (21), the transformations were derived from Bertsimas' robust optimization theory and dual theory. The specific process is as follows:
equation (22) is first obtained by unifying equations (19) - (21). The equivalent formula of equation (22) is equation (23).
The right end of equation (23) is the maximization problem, as in equation (24), and its dual problem is as shown in equation (25).
Wherein v and wiAre respectivelyAnd 0. ltoreq. ziA dual variable of ≦ 1. From the strong dual theorem of linear programming, it can be seen that the problem of equation (24) is bounded, and then equation (25) is bounded as its dual problem, and the two objective function values are consistent. It can therefore be deduced that equation (22) is equivalent to equation (26) and, in turn, equivalent to equation (27).
Finally, a peer constraint of constraints (19) - (21) is obtained, namely the following three sets of equations (28.1) - (30.4), each set of equations containing 4 subformulae.
Thus, the final robust model is transformed to include: an objective function (3), constraint expressions (4) - (9), constraint expressions (28.1) - (30.4), and constraint expressions (13) - (15).
Step (3) solving a design column generation algorithm
And (3.1) decomposing the final charging robust optimization model in the step (2) into a main problem and a sub problem, and accelerating the solving process.
The main problems are as follows: resource limit constraints (6) - (7) of total charging power and charging power in the original robust model are covered, and the task is to satisfy the charging scheme set R of uncertain energy consumption per bus mmTo select a generated charging schedule to minimize the overall charging cost of the fleet. At the same time, the selected strategies all need to meet the resource constraints of the charging station. The mathematical expression is as follows:
xmr∈{0,1} (35)
in the main problem, cmrRefers to the cost of the charging schedule r for the electric bus m. x is the number ofmrIs a binary variable, x if vehicle m selects charging schedule r mr1, otherwise xmr=0。pmtrAnd dmtrAnd p in the definition of (2.1) parametersmtAnd dmtThe expressed meanings are consistent, and the lower subscript R is only used as a scheme set R for distinguishing the vehicles mmDifferent scheme r. Equation (31) aims to minimize the total charge cost of the final selected scheme of the electric bus, consistent with the original robust model objective function. Constraints (32) are constraints on buses that represent each bus only with set RmOne charging scheme of (1) is matched; in the constraint (33), for each time T epsilon T, the accumulated power of all the charging strategies at the time T is required to be less than or equal to pmaxμ, the total power available at the electric bus charging station. The constraint (34) ensures that the number of electric buses simultaneously served by all selected charging strategies at any one time te T does not exceed the total number of charging guns. Equation (35) indicates the variable type.
The sub-problems are: each vehicle is treated as a separate subproblem, the purpose of which is to find a new charging scheme that improves the current optimal solution of the main problem on the basis of the dual variables provided by the main problem. Under policy r, the sub-problem for vehicle M, M ∈ M can be expressed as:
constraint expressions (4) - (5), constraint expressions (8) - (9), constraint expressions (28.1) - (30.4), and constraint expressions (13) - (15).
The objective function of the subproblem is equation (36), consisting of dual variables of charge cost and charge resource, where δmrA dual variable being a constraint equation (32); for charging power at time t, wtA dual variable of constraint formula (33); for the charging gun at time t,. epsilontIs a dual variable of the constraint equation (34). Other constraints are consistent with the original robust model, and the generated model is used for generatingThe electrical scheme still needs to guarantee operation time limitation, charging continuity, battery capacity limitation, uncertain energy consumption requirements and the like.
And (3.2) solving the main problem in the step (3.1) to obtain dual variables for providing the sub-problem.
(3.3) respectively solving the subproblems of different electric buses m, and if the objective function of the subproblems is less than 0, adding a charging scheme set RmUpdating the scheme set of the main problem and returning to the step (3.2); and if the objective functions of the sub-problems are all larger than 0, stopping solving to obtain the lowest charging cost of the robust optimization model of the electric bus fleet under the consideration of the energy consumption uncertainty.
The invention has the beneficial effects that:
the invention provides an electric bus fleet robust charging optimization method based on uncertain operation energy consumption. Firstly, the fluctuation of the operation energy consumption is described by adopting a budget uncertain set, so that the defects of traditional uncertain optimization methods such as random planning and fuzzy planning are overcome, the optimization result is insensitive to parameter disturbance, and the method is very suitable for public transport systems with higher requirements on stability. Secondly, an optimized robust charging model of the electric bus fleet is established, the time-of-use electricity price, charging station resource constraints, uncertain energy consumption requirements of vehicles, the influence of multiple constraints such as operation time on the charging scheme of the electric bus fleet is fully considered, the actual operation condition is fitted, and the practicability is high. And finally, designing a column generation algorithm to solve the model. The column generation algorithm has the advantages of accurate solving result, high calculating speed and the like in the aspect of solving the large-scale mixed integer programming model. In conclusion, the invention provides an efficient and reliable charging scheme for the charging management of the electric bus fleet, and when the actual energy consumption is deviated from the nominal value, the robust charging scheme with uncertain energy consumption can still normally operate in a given battery range, so that the fluctuation of energy consumption is effectively coped with, the safe operation of the bus system is ensured, the charging cost is reduced, and the economical efficiency of the management level is improved.
Drawings
FIG. 1 is a schematic flow diagram of the process of the present invention.
Fig. 2 is an analysis graph of the operation energy consumption per minute of the electric bus in the embodiment of the present invention.
Fig. 3 is a distribution diagram of energy consumption of the electric bus fleet in different electricity price periods in the embodiment of the invention.
Fig. 4 is a graph comparing battery power of the electric bus B1 determination model and robust model pair in the embodiment of the present invention.
Fig. 5 is a graph comparing battery power of the electric bus B2 determination model and robust model pair in the embodiment of the present invention.
Detailed Description
The following describes in detail embodiments of the present invention with reference to examples, and simulates the effects of the present invention.
Fig. 1 is a schematic flow chart of the method of the present invention, as shown in fig. 1, which mainly comprises three steps: constructing an uncertainty set of operation energy consumption; (2) establishing a robust optimization model for charging of the electric bus fleet under uncertain operation energy consumption, wherein the robust optimization technology is mainly used for converting the determined charging optimization model of the electric bus fleet into the robust optimization model for charging under the uncertain operation energy consumption, and further carrying out linearization process on the robust optimization model; (3) and a design column generation algorithm decomposes an original robust model, and a model is quickly solved by utilizing the independence of the main problem and the sub-problem, so that a reliable and economic management scheme is obtained for charging the electric bus fleet.
The specific flow of this embodiment is as follows:
1. method for constructing uncertainty set of energy consumption by adopting certain city operation energy consumption data
As described in step (1), the vehicles m each need to perform an indeterminate energy consumption for shift kIn one interval, i.e.WhereinIs a nominal value of the uncertain energy consumption, representing the maximum deviation value, generally assumed to be p times the nominal value, i.e.Constructed uncertainty set U of energy consumptionbugetCan be defined as:
from the above equation, the key parameters for determining the set of uncertainty in energy consumption includeP and fmk。
In this embodiment, the operation energy consumption historical data of the pure electric bus in a certain city is selected to perform numerical fitting, as shown in fig. 2. The power consumption per minute has fluctuation along with the change of time, the minimum value is 0.142 kW.h, and the maximum value can reach 0.310 kW.h. The power consumption is higher per minute due to traffic jam or work peak, and the average value is close to 0.28 kW.h; the energy consumption value of 17:00-19:30 is also higher and is close to 0.25 kW.h; the operation energy consumption is relatively stable in other periods, and is basically about 0.18 kW.h.
For convenience, the average energy consumption per minute values of different periods of the vehicle in table 1 are collated according to the analysis result of the operation energy consumption of the electric bus. And accumulating the power consumption per minute according to the time period of the vehicle operation shift to obtain the power consumption of the shift.
TABLE 1 average energy consumption per minute values for different periods
I.e. when vehicle m executes shift k, i.e. departure time cmkArrival time smkNominal value of cumulative energy consumptionByAnd calculating to obtain.
The residual parameter ρ directly determines the fluctuation range of the uncertainty interval. RmkIs an important parameter that affects the degree to which a solution is conserved, rmkThe larger the resulting solution is, the more conservative. Note that the parameter rmkIn the embodiment, the difference of the vehicle m is not distinguished, only the difference of the shift is distinguished, and fmkIs satisfied as an integer between 0-k. Since the number of vehicle operation shifts to which the present example relates is at most 8, rmkThe total number of values of (2) is 8. For vehicles with less than 8 operating shifts, only the Gamma consistent with the number of the shift is neededmkThe number of (2). Such as 3 vehicles in a shift, a r actually used in a trial calculationmkThe parameters comprising only fm1,Гm2And rm3. The present embodiment adopts p of 20%, rm1=1,Гm2=2,Гm3=3,Гm4=4,Гm5=5,Гm6=6,Гm7= 7,Гm8These parameters verify the utility of the effectiveness of the electric bus charging robusta optimization model that accounts for energy consumption uncertainty.
2. Setting experimental parameters according to the electric bus fleet charging robust optimization model under uncertain operation energy consumption established in the step (2) and the step (3)
In the experimental example, data of 64 vehicles on 3 lines actually operated in a certain city are selected, including a plurality of parameters such as an operation schedule, charging station resources, battery capacity, time-of-use electricity price and the like, and are specifically shown in tables 2 and 3. The time-of-use electricity price means that the charging price is changed in different time periods. According to the three distribution states of the electricity price, the electricity price is divided into a low-valley electricity price, a flat-peak electricity price and a high-peak electricity price. The price is the lowest between 00:00 and 8:00 and is 0.3 yuan/kW.h.
TABLE 2 parameters of the model for determination of electric bus charging optimization
TABLE 3 time-share electricity price table
Wherein the operation schedule needs to be converted into time, that is, one time Δ t every 5 minutes, 1440 minutes a day is discretized into 288 times t. The first time in the morning, 6:00, is set to be t ═ 1. A sample trip table for an electric bus is shown in table 4.
TABLE 4 example journey table of electric bus
3. Optimization results of robust charging model
The method is adopted to carry out robust charging optimization of the electric bus fleet in consideration of energy consumption uncertainty, and the optimization result shows that the lowest charging cost of 64 fleets is 1539.4 yuan. The total energy consumption required by the electric bus fleet is 3840.2 kW.h, wherein the energy consumption at low power price accounts for 75% and is about 2880 kW.h. According to the low-electricity-price time period of 0: 00-8: 00, the available total power resource is 3040 kW.h, the resource utilization rate of the low-electricity price can be calculated to be as high as 95%, which shows that the result of the robust optimization model can fully utilize the low-electricity-price resource, reduce the charging cost and improve the economy of charging management. Meanwhile, since the charging resource of low power rates is almost occupied, a part of the power consumption needs to be distributed in the flat rate period, about 24% (fig. 3).
In order to further verify the anti-interference performance and stability of the charging scheme obtained by the robust charging optimization method, the condition of determining the energy consumption and the condition of deviation of the energy consumption are considered (rho is 20%), and the battery electric quantity distribution corresponding to the vehicle can be obtained respectively. Fig. 4 and fig. 5 respectively show the comparison of the charging scheme of the optimization model and the charging scheme of the robust optimization model on the electric bus battery power distribution when the energy consumption has a deviation of 20%.
Fig. 4 shows that in the worst case, the electric bus B1 can operate normally in a given battery power range with a minimum power of 50kW · h above the power safety threshold, under the solution of the robust model. Under the solution of the deterministic model, the battery use range of the electric bus exceeds a given range, the lowest electric quantity is only 7kW & h, and the safe operation of the electric bus is seriously influenced.
Fig. 5 illustrates that even if the vehicle is operating with ease, and the battery capacities of the robust model and the deterministic model do not fall below the safety threshold of 20%, at the time of the outbound, vehicle B2 cannot be fully charged according to the deterministic model charging schedule, and the battery capacity is only charged to 84kW · h, which may require a delayed departure, i.e., a departure that is not delayed, may affect the next day's charging schedule. Therefore, when the worst case occurs, the electric bus charging system has two problems in determining the solution of the energy consumption model. First, in extreme cases where energy consumption deviates significantly from the nominal value, the electric bus may have its battery level below a safe threshold before returning to the charging station, which may reduce battery life beyond its given safe level. In severe cases, the battery power may be completely exhausted, thereby causing a safety hazard in the operation process. Second, the electric bus will need to spend longer charging time at the charging station due to the increase of energy consumption, delaying normal operation.
In conclusion, the robust charging optimization method for the electric bus fleet, which considers the uncertainty of energy consumption, can remarkably reduce the charging cost, effectively cope with the fluctuation of energy consumption, and improve the economy and stability of the electric bus charging system.
Claims (1)
1. The robust charging optimization method of the electric bus fleet considering the uncertainty of energy consumption is characterized by comprising the following steps of:
step (1) of constructing an uncertainty set of operation energy consumption
Considering the fluctuation of the energy consumption of the vehicle m in the operation process, a budget uncertain set is introduced to deal with the uncertainty of the operation energy consumption; uncertain energy consumption of each vehicle m required to execute shift kIn a single interval, i.e.WhereinIs a nominal value of the uncertain energy consumption,maximum deviation value representing uncertain energy consumption, p times nominal value, i.e.The budget uncertainty set is defined as shown in formula (1); wherein Z is shown as formula (2);
constructed uncertainty set U of operating energy consumptionbugetCan be defined as:
where M is the total number of vehicles, ΓmkThe energy consumption uncertainty budget for executing shift k is a parameter reflecting the uncertainty level, and the optimality and the robustness of the solution can be balanced; and 0 is more than or equal to gammamkK or less, i.e. gammamkThe maximum value of (1) is energy consumption at the same timeThe maximum number of shifts from its nominal value; when the energy consumption deviation of the vehicle m is maximum for all shifts k, i.e. each ΓmkK, the worst case occurs, with the most conservative charging scheme;
step (2), establishing an electric bus fleet charging robust optimization model under uncertain operation energy consumption
Converting the electric bus charging determination optimization model into a robust optimization model for orderly charging of the electric bus fleet under uncertain operation energy consumption based on the uncertainty set of operation energy consumption defined in the step (1);
(2.1) parameter definition
Collection
M: the charging station serves a set of electric public transport vehicles, M ═ 1,2
T: a set of times divided a day, T ═ 1,2
Km: class set of electric public transport vehicles m, Km={1,2,...,AmIn which A ismIs the total number of shifts of the vehicle m
The set of operating moments for m shifts k of the vehicle,wherein c ismkIs the starting time, s, of the shift k executed by the vehicle leaving the charging stationmkIs the moment when the shift k is finished and returns to the charging station
parameter(s)
ft: electricity price corresponding to time t, unit: Yuan/kW.h
Δ t: unit time, unit: minute (min)
qmk: the remaining battery power, unit, of the vehicle m returning to the charging station after performing shift k: kW.h
Gmk: electric quantity value supplemented after vehicle m performed shift k, unit: kW.h
qmax: battery capacity, unit: kW.h
qmin: minimum threshold allowed for battery charge, unit: kW.h
Δet: electric power consumption of the electric bus per minute, unit: kW.h
Pmax: supplied by charging stationsTotal power, unit: kW (power of kilowatt)
μ: loss rate from charging station to charging gun power transfer process
N: total number of charging guns in charging station
pa: maximum charging power allowed by the battery of the electric bus, unit: kW (power of kilowatt)
pb: maximum charging power allowed by the charging gun, unit: kW (power of kilowatt)
pc: upper limit value of charging power, unit: kW (power of kilowatt)
dmt: if the bus m supplements the electric quantity at the moment t, dmt1 is ═ 1; otherwise dmt=0
(2.2) decision variables
pmt: required charging power at time t of vehicle m, integer variable, unit: kW (power of kilowatt)
Umt: an auxiliary binary variable for charge continuity control, for determining whether a change, U, in the state of charge has occurred from time t to t +1mt∈{0,1}
Vmt: auxiliary binary variable for charge continuity control, for determining whether a change in state of charge, V, has occurred from time t to t-1mt∈{0,1}
Uncertainty set | z of robust model constraint (19)miDual variable corresponding to | less than or equal to 1
Uncertainty set | z of robust model constraint (20)miDual variable corresponding to | less than or equal to 1
Uncertainty set | z of robust model constraint (21)mkDual variable corresponding to | less than or equal to 1
(2.3) establishing a determined charging optimization model taking the minimum total charging cost of the electric bus fleet as an objective function, wherein the model is as follows:
an objective function:
constraint conditions are as follows:
Umt∈{0,1} (13)
Vmt∈{0,1} (14)
0≤pmt≤pc (15)
in the formula, the objective function formula (3) aims to minimize the total charging cost of the electric bus fleet, and is the sum of the electricity prices of all vehicles at all times; the constraint condition formula (4) indicates that the electric bus cannot be supplemented with electric quantity during the operation period, so that the charging power is 0; ensuring that the charging power can not exceed the maximum power p of the charging gun through the constraint condition formula (5)bAnd the maximum power p acceptable for the batterya(ii) a The constraint condition formula (6) and the constraint condition formula (7) are constraint conditions of charging station resource limitation, and the constraint condition formula (6) requires that the total charging power sum of the electric buses which are served simultaneously is not more than the total power available for the charging station; the constraint condition formula (7) requires that the number of the electric buses charged simultaneously is at most the total number of the charging guns; the constraint expressions (8.1) - (8.3) and the constraint expressions (9.1) - (9.5) both represent charging continuity control, that is, once the vehicle starts charging, reaches a predetermined amount of electricity or stops for a predetermined time, and charging is not performed intermittently; the constraint condition formula (10) indicates that the battery electric quantity value after each electricity supplementing should not exceed the battery capacity; the constraint condition formula (11) ensures that the battery capacity of the supplemented electric quantity is not lower than the lowest threshold value after the normal operation of the next shift is finished; the constraint equation (12) means that the amount of electricity to be replenished and the amount of electricity to be consumed areThe electric buses are consistent, namely, the electric buses are fully charged before the electric buses are out of the station on the next day; constraint equations (13) - (15) give the types of variables;
(2.4) converting the determined charging optimization model into a charging robust optimization model under the condition of uncertain operation energy consumption
In order to convert the determined charging optimization model into a robust model considering uncertain energy consumption, constraint conditional expressions (10) - (12) are replaced by constraint conditional expressions (16) - (18), and the model contains uncertain energy consumption parameters;
constraints (16) - (18) may transform the min/max form, as shown in equations (19) - (21):
(2.5) conversion of robust optimization model
To further linearize the constraint equations (19) - (21), the transformations are derived from Bertsimas' robust optimization theory and dual theory; the specific process is as follows:
firstly, unifying the formulas (19) - (21) to obtain a formula (22); the equivalent formula of the formula (22) is formula (23);
the right hand end of equation (23) is the maximization problem, as in equation (24), and its dual problem is as shown in equation (25):
wherein v and wiAre respectivelyAnd 0. ltoreq. ziA dual variable of less than or equal to 1; according to the strong dual theorem of linear programming, the problem of the formula (24) is feasible and bounded, and then the problem of the formula (25) as the dual problem is also feasible and bounded, and the objective function values of the two are consistent; it can therefore be deduced that equation (22) is equivalent to equation (26), and in turn equivalent to equation (27):
finally, equivalent constraints of constraint conditions (19) - (21), namely equations (28.1) - (30.4), are obtained, each set of equations comprising 4 sub-equations:
thus, the final robust model is transformed to include: an objective function (3), constraint expressions (4) - (9), constraint expressions (28.1) - (30.4), and constraint expressions (13) - (15);
step (3) solving a design column generation algorithm
(3.1) decomposing the final charging robust optimization model in the step (2) into a main problem and a sub-problem, and accelerating the solving process;
the main problems are as follows: resource limit constraints (6) - (7) of total charging power and charging power in the original robust model are covered, and the task is to satisfy the charging scheme set R of uncertain energy consumption per bus mmSelecting a generated charging scheme to minimize the overall charging cost of the fleet; meanwhile, the selected strategies all need to meet the resource constraint of the charging station; the mathematical expression is as follows:
xmr∈{0,1} (35)
in the main problem, cmrThe cost of the charging scheme r of the electric bus m is indicated; x is the number ofmrIs a binary variable, x if vehicle m selects charging schedule rmr1, otherwise xmr=0;pmtrAnd dmtrAnd (2.1) p in the parameter definitionmtAnd dmtThe expressed meanings are consistent, and the lower subscript R is only used as a scheme set R for distinguishing the vehicles mmDifferent schemes r in (1); the formula (31) aims to minimize the total charging cost of the scheme finally selected by the electric bus, and is consistent with the objective function of the original robust model; the constraint expression (32) is a constraint condition for buses, and represents that each bus is only connected with the set RmOne charging scheme of (1) is matched; in the constraint condition formula (33), for each time T epsilon T, the accumulated power of all charging strategies at the time T is required to be less than or equal to pmaxMu, i.e. the total power available for electric bus-charging stations; the constraint condition formula (34) ensures that the number of the electric buses simultaneously served by all the selected charging strategies at any time te ∈ T does not exceed the total number of the charging guns; formula (35) indicates the variable type;
the sub-problems are: each vehicle is treated as a separate subproblem, with the aim of finding a new charging scheme that improves the current optimal solution of the main problem on the basis of the dual variables provided by the main problem; under strategy r, the sub-problem of vehicle M, M ∈ M can be expressed as:
constraint expressions (4) - (5), constraint expressions (8) - (9), constraint expressions (28.1) - (30.4), constraint expressions (13) - (15);
the objective function of the sub-problem is equation (36), consisting of dual variables of charging cost and charging resource, where δmrA dual variable being a constraint equation (32); for charging power at time t, wtA dual variable of constraint formula (33); for the charging gun at time t,. epsilontA dual variable being constraint (34); other constraints are consistent with the original robust model, and the generated charging scheme still needs to guarantee operation time limitation, charging continuity, battery capacity limitation, uncertain energy consumption requirements and the like;
(3.2) solving the main problem in the step (3.1) to obtain dual variables for providing to the sub-problems;
(3.3) respectively solving the subproblems of different electric buses m, and if the objective function of the subproblems is less than 0, adding a charging scheme set RmUpdating the scheme set of the main problem and returning to the step (3.2); and if the objective functions of the sub-problems are all larger than 0, stopping solving to obtain the lowest charging cost of the robust optimization model of the electric bus fleet under the consideration of the energy consumption uncertainty.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111150148.XA CN113902182B (en) | 2021-09-29 | 2021-09-29 | Robust charging optimization method for electric bus fleet considering energy consumption uncertainty |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111150148.XA CN113902182B (en) | 2021-09-29 | 2021-09-29 | Robust charging optimization method for electric bus fleet considering energy consumption uncertainty |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113902182A true CN113902182A (en) | 2022-01-07 |
CN113902182B CN113902182B (en) | 2022-09-20 |
Family
ID=79189274
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111150148.XA Active CN113902182B (en) | 2021-09-29 | 2021-09-29 | Robust charging optimization method for electric bus fleet considering energy consumption uncertainty |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113902182B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115330267A (en) * | 2022-10-11 | 2022-11-11 | 北京建筑大学 | Battery charging and replacing facility layout method, device, equipment and medium based on demand behaviors |
CN115688394A (en) * | 2022-10-18 | 2023-02-03 | 上海科技大学 | V2G distribution robust optimization method considering multiple uncertainties of power grid |
CN115689310A (en) * | 2022-11-09 | 2023-02-03 | 东南大学 | Robust evaluation method for resource allocation economy of urban pure electric bus system |
CN115952985A (en) * | 2022-12-21 | 2023-04-11 | 大连理工大学 | Mixed scheduling method of module vehicle and bus in multi-line multi-shift traffic system |
CN116485157A (en) * | 2023-06-16 | 2023-07-25 | 四川国蓝中天环境科技集团有限公司 | Electric bus charging plan optimization method considering charging station vehicle queuing |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110009152A (en) * | 2019-04-03 | 2019-07-12 | 东南大学 | A kind of consideration electricity turns gas and probabilistic regional complex energy system operation robust Optimal methods |
CN111008723A (en) * | 2019-09-24 | 2020-04-14 | 华北电力大学 | Optimization method for design of distributed energy PEV charging station |
CN111654036A (en) * | 2020-05-18 | 2020-09-11 | 天津大学 | Two-stage robust optimization scheduling method for power distribution network considering energy storage quick charging station |
CN111845426A (en) * | 2020-07-01 | 2020-10-30 | 大连理工大学 | Pure electric bus charging power distribution and optimization method based on column generation framework |
CN112347615A (en) * | 2020-10-20 | 2021-02-09 | 天津大学 | Power distribution network hybrid optimization scheduling method considering light storage and fast charging integrated station |
CN112465323A (en) * | 2020-11-19 | 2021-03-09 | 大连理工大学 | Short-term robust scheduling method for cascade hydropower station coupling daily electric quantity decomposition and day-ahead market bidding |
-
2021
- 2021-09-29 CN CN202111150148.XA patent/CN113902182B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110009152A (en) * | 2019-04-03 | 2019-07-12 | 东南大学 | A kind of consideration electricity turns gas and probabilistic regional complex energy system operation robust Optimal methods |
CN111008723A (en) * | 2019-09-24 | 2020-04-14 | 华北电力大学 | Optimization method for design of distributed energy PEV charging station |
CN111654036A (en) * | 2020-05-18 | 2020-09-11 | 天津大学 | Two-stage robust optimization scheduling method for power distribution network considering energy storage quick charging station |
CN111845426A (en) * | 2020-07-01 | 2020-10-30 | 大连理工大学 | Pure electric bus charging power distribution and optimization method based on column generation framework |
CN112347615A (en) * | 2020-10-20 | 2021-02-09 | 天津大学 | Power distribution network hybrid optimization scheduling method considering light storage and fast charging integrated station |
CN112465323A (en) * | 2020-11-19 | 2021-03-09 | 大连理工大学 | Short-term robust scheduling method for cascade hydropower station coupling daily electric quantity decomposition and day-ahead market bidding |
Non-Patent Citations (3)
Title |
---|
刘锴等: "基于多层混合效应的电动汽车能耗估计模型", 《武汉理工大学学报(交通科学与工程版)》 * |
刘锴等: "电动汽车充电站布局优化研究", 《城市交通》 * |
孙小慧等: "考虑时空间限制的电动汽车充电站布局模型", 《地理科学进展》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115330267A (en) * | 2022-10-11 | 2022-11-11 | 北京建筑大学 | Battery charging and replacing facility layout method, device, equipment and medium based on demand behaviors |
CN115688394A (en) * | 2022-10-18 | 2023-02-03 | 上海科技大学 | V2G distribution robust optimization method considering multiple uncertainties of power grid |
CN115688394B (en) * | 2022-10-18 | 2023-12-26 | 上海科技大学 | V2G distribution robust optimization method considering multiple uncertainties of power grid |
CN115689310A (en) * | 2022-11-09 | 2023-02-03 | 东南大学 | Robust evaluation method for resource allocation economy of urban pure electric bus system |
CN115689310B (en) * | 2022-11-09 | 2024-06-04 | 东南大学 | Robust evaluation method for resource allocation economy of urban pure electric bus system |
CN115952985A (en) * | 2022-12-21 | 2023-04-11 | 大连理工大学 | Mixed scheduling method of module vehicle and bus in multi-line multi-shift traffic system |
CN115952985B (en) * | 2022-12-21 | 2023-08-18 | 大连理工大学 | Mixed scheduling method of module vehicle and bus in multi-line multi-shift traffic system |
CN116485157A (en) * | 2023-06-16 | 2023-07-25 | 四川国蓝中天环境科技集团有限公司 | Electric bus charging plan optimization method considering charging station vehicle queuing |
CN116485157B (en) * | 2023-06-16 | 2023-08-22 | 四川国蓝中天环境科技集团有限公司 | Electric bus charging plan optimization method considering charging station vehicle queuing |
Also Published As
Publication number | Publication date |
---|---|
CN113902182B (en) | 2022-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zheng et al. | Integrating plug-in electric vehicles into power grids: A comprehensive review on power interaction mode, scheduling methodology and mathematical foundation | |
CN113902182B (en) | Robust charging optimization method for electric bus fleet considering energy consumption uncertainty | |
Wu et al. | Hierarchical operation of electric vehicle charging station in smart grid integration applications—An overview | |
Zheng et al. | A novel real-time scheduling strategy with near-linear complexity for integrating large-scale electric vehicles into smart grid | |
Yang et al. | Optimal dispatching strategy for shared battery station of electric vehicle by divisional battery control | |
Kumar et al. | V2G capacity estimation using dynamic EV scheduling | |
Yu et al. | A real time energy management for EV charging station integrated with local generations and energy storage system | |
CN110598904B (en) | Vehicle network energy interaction optimization method considering renewable energy consumption in market environment | |
Zheng et al. | Smart charging algorithm of electric vehicles considering dynamic charging priority | |
Ni et al. | Hierarchical optimization of electric vehicle system charging plan based on the scheduling priority | |
Jaiswal et al. | Optimal load management of plug-in electric vehicles with demand side management in vehicle to grid application | |
Li et al. | Multi-objective optimal operation of centralized battery swap charging system with photovoltaic | |
Saldanha et al. | Control strategies for smart charging and discharging of plug-in electric vehicles | |
Athulya et al. | Electric vehicle recharge scheduling in a shopping mall charging station | |
Marasciuolo et al. | Chance-constrained calculation of the reserve service provided by EV charging station clusters in energy communities | |
Yin et al. | Optimal scheduling of electric vehicle integrated energy station using a novel many-objective stochastic competitive optimization algorithm | |
Das et al. | Game theoretical energy management of EV fast charging station with V2G capability | |
Diaz-Londono et al. | Coordination strategies for electric vehicle chargers integration in electrical grids | |
Marasciuolo et al. | Two-stage scheduling of electrical vehicle charging station clusters in a community of DC microgrids | |
CN110570098A (en) | Electric automobile charging and battery replacing station control method considering battery replacing demand and photovoltaic uncertainty | |
Yu et al. | Optimal battery swapping and charging strategy considering on-site solar generation | |
Arcos-Aviles et al. | Fuzzy control-based energy management system for interconnected residential microgrids using the forecasts of power generation and load demand | |
Chang et al. | A coordinated charging strategy for PV-Assisted charging station of electric vehicles based on charging service price | |
Uko | Optimization of Vehicle to Grid System in a Power System with Unit Commitment | |
Aziz et al. | Vehicle to grid (v2g) for peak shaving: New trend, benefits, and issues |
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 |