CN113902182B - 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
- CN113902182B CN113902182B CN202111150148.XA CN202111150148A CN113902182B CN 113902182 B CN113902182 B CN 113902182B CN 202111150148 A CN202111150148 A CN 202111150148A CN 113902182 B CN113902182 B CN 113902182B
- Authority
- CN
- China
- Prior art keywords
- charging
- energy consumption
- constraint
- robust
- vehicle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005265 energy consumption Methods 0.000 title claims abstract description 96
- 238000005457 optimization Methods 0.000 title claims abstract description 67
- 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 34
- 230000014509 gene expression Effects 0.000 claims description 29
- 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 the robustness of an obtained charging optimization result and avoiding the limitation that the traditional uncertain optimization method presupposes probability distribution of uncertain parameter clothes. 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 energy consumption uncertainty 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 electronics & Transportation Review,2017,100: 115-. 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, and Pan Jing Hui, an electric bus ordered charging strategy [ J ] based on double-layer optimization, a power grid and clean energy source, 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 basic electric bus system [ J ]. Transportation Research Part E: logics and Transportation Review,2020, 142 ] based on given capacity of a power distribution network, the number of chargers and bus operation information, so as to optimize the charging scheduling and management of the electric bus rapid charging system and effectively reduce 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 compiled fleet charging strategy 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 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 electric bus is out of the station, so that the delayed departure of the electric bus and the like may occur, 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 using a robust optimization technology to depict uncertainty of operation energy consumption, on one hand, the robustness of the obtained charging optimization result is guaranteed, namely, the obtained solutions are feasible as long as the 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 clothes 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 per vehicle m required to execute shift kIn a single interval, i.e.WhereinIs a nominal value for 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 budget uncertainty set definition is shown in equation (1). Wherein Z is shown as formula (2).
Constructed uncertainty set U of operating energy consumption buget Can be defined as:
wherein r mk The 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 should be noted that 0 ≦ r mk K ≦ k, i.e. r mk The 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 r mk K, which means that the worst case occurs, the charging scheme is most conservative.
Step (2) establishing 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 obtains a charging scheme with higher robustness and optimality.
(2.1) parameter definition
Collection
M is the number 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 }
K m Class set of electric public transport vehicles m, K m ={1,2,…,A m In which A is m Is 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 is mk Is the starting time, s, of the shift k executed by the vehicle leaving the charging station mk Is the moment when the shift k is finished and returns to the charging station
parameter(s)
f t The 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
q mk The remaining battery power of the vehicle m returning to the charging station after the shift k is executed (unit: kW. h)
q max Battery capacity (unit: kW. h)
q min Minimum threshold allowed by battery capacity (unit: kW. h)
Δe t Power consumption of electric bus per minute (unit: kW. h)
P max Total 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
p a Maximum allowable charging power for battery of electric bus (unit: kW)
p b Maximum charging power allowed by charging gun (unit: kW)
p c Upper limit of charging power (unit: kW)
d mt If the bus m supplements the electric quantity at the moment t, d mt 1 is ═ 1; otherwise d mt =0
(2.2) decision variables
p mt Required charging function at time t of m for vehicleRate, integer variable, (unit: kW)
U mt Auxiliary binary variable for charge continuity control, for determining whether a change in the state of charge occurs from time t to t +1, U mt ∈{0,1}
V mt Auxiliary binary variable for charge continuity control, for determining whether a change in the state of charge has occurred from time t to t-1, V mt ∈{0,1}
Uncertainty set | z of robust model constraint (19) mi Dual variable corresponding to | less than or equal to 1
Uncertainty set | z of robust model constraint (20) mi Dual variable corresponding to | less than or equal to 1
Uncertainty set | z of robust model constraint (21) mk Dual 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:
U mt ∈{0,1} (13)
V mt ∈{0,1} (14)
0≤p mt ≤p c (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 the electricity rates of the consumed electric quantities of all vehicles at all times. The constraint condition formula (4) indicates that the electric bus can not be supplemented with electric quantity during the operationTherefore, 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) b And the maximum power p acceptable for the battery 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 simultaneously served is not more than the total power available for the charging stations; 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. Constraint expressions (8.1) to (8.3) and constraint expressions (9.1) to (9.5) both indicate charging continuity control, in which once charging of the vehicle is started, the vehicle reaches a predetermined amount of electricity or stops for a predetermined time, and charging is not performed intermittently. Note that the form of the constraints (9.1) - (9.5) is relatively complex, since the charging period after the last shift has been performed includes two parts, namelyThe 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. The constraint equations (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 can be converted into min/max forms, such as 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 w i Are respectivelyAnd 0. ltoreq. z i A 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 is also 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 m m To 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:
x mr ∈{0,1} (35)
in the main problem, c mr Refers to the cost of the charging schedule r for the electric bus m. x is the number of mr Is a binary variable, x if vehicle m selects charging schedule r mr 1, otherwise x mr =0。p mtr And d mtr And (2.1) p in the parameter definition mt And d mt The expressed meanings are consistent, and the lower subscript R is only used as a scheme set R for distinguishing the vehicles m m Different schemes 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 R m One charging scheme of (1) is matched; restraint (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 p max μ, the total power available for 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 T e 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 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), and 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 δ mr A dual variable being a constraint equation (32); for charging power at time t, w t A dual variable of constraint formula (33); for the charging gun at time t,. epsilon t Is a dual variable of the constraint equation (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.
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 R m Updating 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 electric bus fleet robust charging optimization model is established, the influence of multiple constraints such as time-of-use electricity price, charging station resource constraint, uncertain energy consumption requirements of vehicles, operation time and the like on a charging scheme of the electric bus fleet is fully considered, actual operation conditions are 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 dealt 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: (1) 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 considering the charging under uncertain operation energy consumption, and further carrying out linearization process on the robust optimization model; (3) the design column generation algorithm decomposes an original robust model, and the model is quickly solved by utilizing the independence of main problems and sub problems, 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 a single 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 consumption buget Can be defined as:
from the above equation, the key parameters for determining the set of uncertainty in energy consumption includeP and f mk 。
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 integrally higher per minute probably due to traffic jam or work rush hour and the like, 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 the vehicles in different periods of time in table 1 are collated according to the analysis results of the operation energy consumption of the electric buses. And accumulating the power consumption per minute according to the time period of the operation shift of the vehicle 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 c mk Arrival time s mk Nominal value of cumulative energy consumptionByAnd calculating to obtain.
The residual parameter ρ directly determines the fluctuation range of the uncertainty interval. R mk Is an important parameter that affects the degree to which a solution is conserved, r mk The larger the resulting solution is, the more conservative. Note that the parameter r mk In the embodiment, the difference of the vehicle m is not distinguished, only the difference of the shift is distinguished, and f mk Is satisfied as an integer between 0-k. Since the number of vehicle operation shifts to which this example relates is at most 8, f mk The total number of values of (2) is 8. For vehicles with less than 8 operating shifts, the vehicles only need to adopt Gamma consistent with the number of shifts mk The number of (2). Such as 3 vehicles in a shift, a r actually used in a trial calculation mk The parameters comprising only f m1 ,Г m2 And r m3 . The present embodiment adopts p of 20%, r m1 =1,Г m2 =2,Г m3 =3,Г m4 =4,Г m5 =5,Г m6 =6,Г m7 =7,Г m8 These 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 VALUE METER
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.2kW & h, wherein the energy consumption with low electricity price accounts for 75% and is about 2880kW & h. According to the fact that the available total power resource in the low-electricity-price time period from 0:00 to 8:00 is 3040 kW.h, the resource utilization rate of the low-electricity price can be calculated to be as high as 95%, and the result of the robust optimization model can fully utilize the resource of the low-electricity price, reduce charging cost and improve 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 battery power distribution corresponding to the vehicle can be obtained by considering the condition of determining the energy consumption and the condition of deviation of the energy consumption (rho is 20%). 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 charge range with a minimum charge of 50kW · h above the charge safety threshold, in a solution to 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 relatively easy to operate, the battery capacities of the robust model and the deterministic model do not fall below the safety threshold of 20%, but at the time of the departure, the vehicle B2 cannot be fully charged according to the charging schedule of the deterministic model, the battery capacity is charged only to 84kW · h, the departure may be delayed, and the charging schedule of the next day may be affected even if the departure is not delayed. 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 range. 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 consumption buget Can be defined as:
where M is the total number of vehicles, r mk The 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 less than or equal to r mk K is less than or equal to k, i.e. r mk The maximum value of (a) is the maximum number of shifts of which the energy consumption deviates from the nominal value simultaneously; when the deviation value of the energy consumption of the vehicle m in all shifts k is maximum, namely, each shift r mk K, the worst case occurs, with the most conservative charging scheme;
step (2), establishing a robust optimization model for orderly charging of electric bus fleets 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 the operation energy consumption defined in the step (1);
(2.1) parameter definition
Collection of
M is the number set of electric buses served by the charging station, wherein M is {1,2, …, M }
T is a set of times divided a day, {1,2, …,1440/Δ T }
K m Shift set of vehicle m, K m =P1,2,…,A m In which A is m Is the total number of shifts of the vehicle m
The set of operating hours for m shift k of the vehicle,where cmk is the starting time of the vehicle leaving the charging station to execute shift k, s mk Is the moment when the shift k is finished and returns to the charging station
parameter(s)
f t : electricity price corresponding to time t, unit: Yuan/kW.h
Δ t: unit time, unit: minute (min)
q mk : the remaining battery power, unit, of the vehicle m returning to the charging station after performing shift k: kW. h
G mk : electric quantity value supplemented after vehicle m performed shift k, unit: kW.h
q max : battery capacity, unit: kW.h
q min : minimum threshold allowed for battery charge, unit: kW.h
Δe t : the power consumption of the electric bus per minute, unit: kW.h
P max : total power supplied by the charging station, in units: kW (power of kilowatt)
μ: loss rate from charging station to charging gun power transfer process
N: total number of charging guns in charging station
p a : maximum charging power allowed by the battery of the electric bus, unit: kW (kilo power)
p b : maximum charging power allowed by the charging gun, unit: kW (power of kilowatt)
p c : upper limit value of charging power, unit: kW (power of kilowatt)
d mt : if the vehicle m is replenished with electricity at time t, d mt 1; otherwise d mt =0
(2.2) decision variables
p mt : required charging power at time t of vehicle m, integer variable, unit: kW (power of kilowatt)
U mt : 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 +1 mt ∈{0,1}
V mt : auxiliary binary variable for charge continuity control, for determining whether a change in state of charge, V, has occurred from time t to t-1 mt ∈{0,1}
Uncertainty concentration of robust model constraint conditional expression (19)Corresponding dual variable
Uncertainty set | z of robust model constraint conditional expression (19) mi Dual variable corresponding to | less than or equal to 1
Uncertainty concentration of robust model constraint conditional expression (20)Corresponding dual variable
Uncertainty set | z of robust model constraint conditional expression (20) mi Dual variable corresponding to | less than or equal to 1
Uncertainty concentration of robust model constraint conditional expression (21)Corresponding dual variable
Uncertainty set | z of robust model constraint conditional expression (21) mk Dual 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:
U mt ∈{0,1}(13)
V mt ∈{0,1} (14)
0≤p mt ≤p c (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 operation, 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) b And the maximum power p acceptable for the battery a (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, it 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 supplied and the amount of electricity consumed are consistent, that is, the electric bus must be fully charged before the next day of departure; constraint equations (13) - (15) give the types of variables;
(2.4) converting the determined charging optimization model into a robust optimization model for orderly charging of electric bus fleets under 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 w i Are respectivelyAnd 0. ltoreq. z i A 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 a robust optimization model for orderly charging of the electric bus fleet under uncertain operation energy consumption 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 the total charging power and the charging gun in the original robust model are covered, and the task of the resource limit constraints is to satisfy the charging scheme set R of the uncertain energy consumption per vehicle m m Selecting 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:
x mr ∈{0,1}(35)
in the main problem, c mr Refers to the cost of the charging schedule r for the vehicle m; x is the number of mr Is a binary variable, x if vehicle m selects charging schedule r mr 1, otherwise x mr =0;p mtr The required charging power of the vehicle at the moment t of m in the charging scheme r is expressed as unit: kW; d is a radical of mtr Means that in the charging scheme r, the vehicle m supplements the electric quantity at the moment t, then d mtr 1, otherwise mtr 0; the lower corner mark R is only used as a scheme set R for distinguishing the vehicles m m Different charging schemes r; equation (31) aims to minimize the charging of the final selected scheme of the electric busThe total cost is consistent with the target 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 R m One 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 p max Mu, 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 the charging scheme r, the sub-problem of the 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 δ mr A dual variable being a constraint equation (32); for charging power at time t, w t A dual variable of constraint formula (33); for the charging gun at time t,. epsilon t A 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 sub-problems of different vehicles m, and if the objective function of the sub-problem is less than 0, adding a charging scheme set R m Updating the scheme set of the main problem and returning to the step (3.2); if the objective functions of the subproblems are all greater than 0, stoppingAnd solving to obtain the lowest charging cost of the robust optimization model for orderly charging of the electric bus fleet under the condition of considering the uncertain operation energy consumption of the energy consumption.
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 CN113902182A (en) | 2022-01-07 |
CN113902182B true 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) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114997478A (en) * | 2022-05-19 | 2022-09-02 | 郑州天迈科技股份有限公司 | New energy bus charging scheduling method based on hybrid multi-target particle swarm algorithm |
CN115330267B (en) * | 2022-10-11 | 2023-04-07 | 北京建筑大学 | Battery charging and replacing facility layout method, device, equipment and medium based on demand behaviors |
CN115688394B (en) * | 2022-10-18 | 2023-12-26 | 上海科技大学 | V2G distribution robust optimization method considering multiple uncertainties of power grid |
CN115689310B (en) * | 2022-11-09 | 2024-06-04 | 东南大学 | Robust evaluation method for resource allocation economy of urban pure electric bus system |
CN115952985B (en) * | 2022-12-21 | 2023-08-18 | 大连理工大学 | Mixed scheduling method of module vehicle and bus in multi-line multi-shift traffic system |
CN116485157B (en) * | 2023-06-16 | 2023-08-22 | 四川国蓝中天环境科技集团有限公司 | Electric bus charging plan optimization method considering charging station vehicle queuing |
Citations (5)
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 |
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 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111008723B (en) * | 2019-09-24 | 2023-12-22 | 华北电力大学 | Optimization method for design of distributed energy PEV charging station |
-
2021
- 2021-09-29 CN CN202111150148.XA patent/CN113902182B/en active Active
Patent Citations (5)
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 |
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 (2)
Title |
---|
基于多层混合效应的电动汽车能耗估计模型;刘锴等;《武汉理工大学学报(交通科学与工程版)》;20151215(第06期);全文 * |
考虑时空间限制的电动汽车充电站布局模型;孙小慧等;《地理科学进展》;20120615(第06期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113902182A (en) | 2022-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
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. | Integrating plug-in electric vehicles into power grids: A comprehensive review on power interaction mode, scheduling methodology and mathematical foundation | |
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 | |
CN112488444B (en) | Electric vehicle fast and slow synchronous ordered charging scheduling method and electric quantity settlement method | |
Yu et al. | A real time energy management for EV charging station integrated with local generations and energy storage system | |
Mohammed et al. | Interruptible charge scheduling of plug-in electric vehicle to minimize charging cost using heuristic algorithm | |
CN116151486B (en) | Multi-time-scale random optimization method and device for photovoltaic charging station with 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 | |
Jaiswal et al. | Optimal load management of plug-in electric vehicles with demand side management in vehicle to grid application | |
Ni et al. | Hierarchical optimization of electric vehicle system charging plan based on the scheduling priority | |
Saldanha et al. | Control strategies for smart charging and discharging of plug-in electric vehicles | |
CN113799640A (en) | Energy management method suitable for microgrid comprising electric vehicle charging pile | |
Deng et al. | Hierarchical operation management of electric vehicles for depots with PV on-site generation | |
Kim | Smart charging architecture for between a plug-in electrical vehicle (PEV) and a smart home | |
Amin et al. | Renewable energy maximization for pelagic islands network of microgrids through battery swapping using deep reinforcement learning | |
Huang et al. | Recharging and transportation scheduling for electric vehicle battery under the swapping mode | |
Wang et al. | Admission and scheduling mechanism for electric vehicle charging with renewable energy | |
Diaz-Londono et al. | Coordination strategies for electric vehicle chargers integration in electrical grids | |
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 | |
Liu et al. | Operation-area-constrained adaptive primary frequency support strategy for electric vehicle clusters |
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 |