CN115099702A - Electric bus daytime running charging optimization method based on Lagrange relaxation algorithm - Google Patents
Electric bus daytime running charging optimization method based on Lagrange relaxation algorithm Download PDFInfo
- Publication number
- CN115099702A CN115099702A CN202210874670.0A CN202210874670A CN115099702A CN 115099702 A CN115099702 A CN 115099702A CN 202210874670 A CN202210874670 A CN 202210874670A CN 115099702 A CN115099702 A CN 115099702A
- Authority
- CN
- China
- Prior art keywords
- station
- vehicle
- charging
- line
- bus
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000009977 dual effect Effects 0.000 claims description 12
- 230000008878 coupling Effects 0.000 claims description 10
- 238000010168 coupling process Methods 0.000 claims description 10
- 238000005859 coupling reaction Methods 0.000 claims description 10
- 230000000295 complement effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 238000011478 gradient descent method Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 description 5
- 238000005265 energy consumption Methods 0.000 description 4
- 239000000446 fuel Substances 0.000 description 3
- 150000001875 compounds Chemical class 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 238000010792 warming Methods 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
- B60L53/00—Methods of charging batteries, specially adapted for electric vehicles; Charging stations or on-board charging equipment therefor; Exchange of energy storage elements in electric vehicles
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- Biodiversity & Conservation Biology (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Mechanical Engineering (AREA)
- Transportation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Primary Health Care (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Power Engineering (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
Abstract
On the premise of ensuring a vehicle scheduled schedule, an electric bus daytime running charging optimization method based on a Lagrange relaxation algorithm aims at minimizing the total charging time of a system, obtains an integer optimization original problem model of a full-line network level charging scheme based on the electric quantity state of a discrete battery, converts the integer optimization original problem model into a Lagrange relaxation problem, decomposes the integer optimization original problem model into a plurality of subproblems, converts the integer optimization original problem model into a plurality of compact subproblems, obtains the optimal solution of a reduced original problem model P by solving the optimal solution of the plurality of compact subproblems, accurately describes the charging problem of a large-scale bus network, reduces the number of decision variables and constraint conditions in the model step by step, reduces the problem scale, improves the solving efficiency and carries out fast and efficient solution.
Description
Technical Field
The invention relates to an electric bus daytime running and charging optimization method based on a Lagrange relaxation algorithm.
Background
Urbanization and motorization lead to the rapid rise of vehicle ownership, which causes a great deal of environmental (such as noise and pollutant emission) and social (such as energy consumption and traffic jam) problems, the energy consumption generated in the traffic field reaches more than 25% of the total global energy consumption, the emission of greenhouse gases accounts for 27%, and the use of fuel vehicles is considered to be a serious threat to global warming.
Electric vehicles are used as substitutes for fuel vehicles, and are of great help to reduce the negative impact of urban traffic on climate and improve air quality conditions. Over the past decades, many cities around the world have provided significant ongoing investment in urban traffic electrification, with significant success particularly in public traffic electrification.
Compared with a fuel bus, the electric bus has the greatest defect that the limited battery capacity cannot meet the requirement of operation energy consumption all day long, and electricity needs to be supplemented in the daytime operation process; the non-linear relationship between the battery charging time and the electric quantity state brings a challenge to the control of the electric quantity of the bus journey.
However, in the existing research, the non-linear charging process of the electric bus and the electric quantity state tracking in the process of the journey are not considered in a fusion manner, and a quick and effective solving method is not formed by coordination of the bus operation schedule and the charging time sequence in a large-scale network.
Disclosure of Invention
The technical problem to be solved by the invention is to provide an electric bus daytime running and charging optimization strategy based on a Lagrange relaxation algorithm, according to a vehicle operation schedule and road network information, the electric quantity state of a vehicle from each station is tracked and recorded by depicting a nonlinear process of battery charging, an electric bus daytime running and charging optimization model based on the Lagrange relaxation algorithm is established, and an efficient solving algorithm is designed.
In order to solve the technical problem, the inventionThe technical scheme adopted by the invention is as follows: an electric bus daytime running charging optimization method based on a Lagrange relaxation algorithm comprises the following specific steps: s1, obtaining electric bus data and line operation data containing the number of lines, constructing a full-line network-level charging scheme integer optimization model based on the electric quantity state of a discrete battery, and obtaining an original problem model P containing coupling constraint; s2, converting the coupling constraint into an objective function of an original problem model P by introducing a Lagrange multiplier, so as to obtain a Lagrange relaxation problem P' (mu, lambda) with a new constraint condition; s3, decomposing the Lagrangian relaxation problem P' (mu, lambda) into a plurality of sub-problems P according to the line number l ' new decision variables are introduced, i.e. whether the vehicles of the respective line are charged at the respective station, so that a plurality of said sub-questions P are presented l ' conversion to multiple compact sub-problems P l "; s4, solving a plurality of compact sub-problems P l ", solving the target function of P' (mu, lambda), i.e. solving the target function of the original problem model P, and converting the optimal solution for solving the lower bound of the original problem model P into the solution of the Lagrangian dual problem P LD Solving the Lagrangian dual problem P by using a subgradient descent method LD And after each iterative solution, updating a Lagrange multiplier according to the current solution, and calculating the upper-bound optimal solution of the original problem model P by adopting a k-greedy algorithm until the upper-bound optimal solution of the original problem model P is equal to the lower-bound optimal solution, wherein the algorithm reaches a termination condition, the optimal solution of the original problem model P is obtained, namely, a daytime running charging scheme of each electric bus is obtained, and the optimization of daytime running charging of the electric buses is completed.
In some embodiments, the electric bus data includes the maximum charging times of each bus at each station, the minimum electric quantity protection value of each bus, and the electric bus electric quantity state, and the route operation data includes the route station, the arrival and departure schedule of each station, and the number of charging piles at each station.
In some embodiments, the objective function of the original problem model P in step S1 is specifically:
s.t.
wherein,
s.t. represents the constraint condition is;
n and L respectively represent the position and the number of bus stations, L represents the L-th bus line, L belongs to L, N represents the N-th bus station, N belongs to N, and N l Represents the total number of stations of route l;
a l,j scheduled arrival time, d, at jth stop for vehicles representing route l l,j The vehicles on the line l are scheduled to leave the station at the jth station for a predetermined time, j e { 1, …, N l -1 }, discretizing the planning time interval T into intervals in units of Δ tmin, T, T', T "each being a certain time point within the planning time interval T, the charging scheme being determined by a decision variable u l,j tt’ Or u l,j t’t” Is expressed when the vehicle is in time interval t, t']Or [ t', t "]When charging u is positive at jth station of line l l,j tt’ =1 or u l,j t’t” =1;
S is a discrete set, S represents an electric bus electric quantity state SOC, S = {0,1,2, …, Q }, each value S e S in the set S represents an electric quantity state reached after the vehicle starts to be charged from the SOC =0 state for S time intervals Δ t, i.e., S =0 represents an electric quantity reached by the battery after the vehicle is charged for Δ t duration, S =2 represents an electric quantity reached by the battery after the vehicle is charged for 2 Δ t duration, time required for the battery to be fully charged from the state of SOC =0 is H hours, Q represents an electric quantity state Q, i.e., a full state, and Q =60H Δ t;
station j =1, …, N for all routes of route l l Of variable v l,j,s Whether the electric quantity state of the vehicle is s or not when the vehicle leaves the jth station is shown, when the variable is 1, the electric quantity state is shown to be s, and when the variable is 0, the electric quantity state is not shown to be sQuantity v l,1,Q =1 denotes the state of charge of the vehicle on route l leaving the 1 st station is Q, variable v l,j-1,s Indicating whether the electric quantity state of the vehicle when leaving the j-1 th station is s or not;
g is an electric public transport network, G = (N, L), and the minimum electric quantity protection value of the vehicle of each line L is G l ;N l Represents the total number of stations of the route l, N (1, j) represents the jth station of the route l, N (1, j) is equal to N l ;
D l,j Representing the amount of power consumed by a vehicle on route l for its journey from jth to j +1 th stop, D l,j-1 Represents the amount of power consumed by the vehicle on route l for its trip from the jth-1 st station to the jth station;
the maximum number of times of charging permitted by the vehicles on the line l at each station is U l And the number of charging piles of each station N belonging to N is M n 。
In some embodiments, the lagrangian relaxation problem P' (μ, λ) is specifically:
s.t.
wherein, mu l,j And λ n,t Is the lagrange multiplier described in step S2.
In some embodiments, the sub-problem has | L | sub-problems P l Is a compound ofA variable 0 to 1 anda bar-constrained 0-1 integer programming model, where | S | represents the number of elements in a discrete set S, a l,n Indicating the scheduled arrival time of the vehicle at the n stations on the route l, d l,n Indicating the scheduled departure time of the vehicle on the line l at the n stations, the l sub-question P l ' described as:
s.t.
wherein the charging scheme is determined by a decision variable u l,n tt’ Or u l,n t’t” Denotes when the vehicle is in time interval t, t']Or [ t', t "]Then, u is charged at the nth station of the line l l,n tt’ =1 or u l,n t’t” =1。
In some embodiments, the sub-problem P is presented to l ' New decision variables are introduced into,Represents [ t, t +1 ]]When the vehicle on line l is charging at its jth station, and the vehicle on line l is not charging at jth station,wherein j =1, …, N l ,t∈{a l,j ,…,d l,j H, solve the subproblem P l ' equivalent transformation compact sub-problem P l ", compact sub-problem P l "have | L | number, L th compact sub-problem P l "may be described as:
s.t.
wherein, variableRepresenting the battery state of charge of the vehicle on line l as it leaves the jth station,representing the battery state of charge of the vehicle on line i when it leaves the 1 st station,representing the battery state of charge of the vehicle on line i as it leaves the j-1 st station,is a new decision variable, which when it is 1, indicates [ t, t +1 ]]The vehicle on line l is charging at its j-1 st station.
In some embodiments, the step S4 is specifically: step S41, let the number of charging times U of the vehicle on the route l at the jth station l,j Number of vehicles charged simultaneously m in each station t interval =0 n,t = 0; step S42, traversing interval t epsilon [ a ∈ ] l,j +1,d l,j ]If at allAnd isIn time, let U l,j + 1; traversal interval t e [ a ] l,j ,d l,j ],(ii) a Step S43, and obtaining,N is N, T is T, the target function of the Lagrangian relaxation problem P' (mu, lambda) is
(ii) a The original problem model P has a target function formula of(ii) a Step S44, respectively using mu and lambda to represent vectors formed by Lagrange multipliers in the Lagrange relaxation problem P '(mu, lambda), and enabling L (mu, lambda) to represent the optimal solution corresponding to the Lagrange relaxation problem P' (mu, lambda), so that the solution of the lower bound of the original problem model P can be converted into the solution of the Lagrange dual problem P LD :Step S45, utilizing a sub-gradient algorithm to solve the Lagrangian dual problem P LD And (6) solving.
In some embodiments, step S45 specifically includes: step S01, the Lagrangian dual problem P is processed LD Initializing, setting the current iteration number i =0, initializing Lagrangian multipliers mu and lambda, and setting an initial lower bound lb = - ∞ and an initial upper bound ub = + ∞; step S02, update multipliers mu, lambda and solve compact subproblems P using a commercial solver l ", obtain the optimal solution,(ii) a Step S03, calculating U by using the steps S41-S44 l,j (i) And m n,t (i) And obtaining P l "objective function value obj of LR (i) Solving the original problem by using a k-greedy algorithm and obtaining an objective function value obj corresponding to a feasible solution (i) Recording the current iteration result, if the constraint is relaxed, namely U l,j (i) ≤U l And m n,t (i) ≤M n If both are satisfied, the upper and lower bounds are checked and updated: if obj LR (i) If lb is greater than lb, update lower bound lb = obj LR (i) (ii) a If obj (i) If the sub is less than the sub, the upper bound ub is updated to be = obj (i) (ii) a Step S04, updating the lagrangian multiplier by using the secondary gradient, making the update step of the multiplier δ =1/2(n +1), and updating the multiplier by using the following formula: mu.s l,j (i+1) =max{μ l,j (i) +δ·△μ l,j (i) ,0},λ n,t (i+1) =max{λ n,t (i) +δ·△λ n,t (i) 0} in which
Step S05, judging whether the algorithm termination condition is satisfied, if all the complementary relaxation conditions are satisfied, the μ l,j (i+1) (U l,j (i) -U l ) =0 and λ n,t (i+1) (m n,t (i) -M n ) =0, the algorithm is terminated, and at this time, the optimal solution corresponding to the lagrangian relaxation problem is the optimal solution of the original problem, and the superscripts (i) and (i +1) respectively represent the ith iteration and the (i +1) th iteration; when lb = ub, the algorithm is terminated, and the optimal solution corresponding to the Lagrangian relaxation problem is the optimal solution of the original problem; when the iteration number I exceeds the maximum iteration number I max When the algorithm is ended, outputting a target function and an optimal solution corresponding to the upper and lower bounds; if the three conditions are not met, the step S01-the step S05 are circulated, the iteration times are +1 in sequence, the optimal solution of the original problem is obtained until one of the three conditions is met, the daytime running charging scheme of each electric bus is obtained, and the daytime running charging of the electric buses is optimized.
The scope of the present invention is not limited to the specific combinations of the above-described features, and other embodiments in which the above-described features or their equivalents are arbitrarily combined are also intended to be encompassed. For example, the above features and the technical features (but not limited to) having similar functions disclosed in the present application are mutually replaced to form the technical solution.
Due to the application of the technical scheme, compared with the prior art, the invention has the following advantages: the invention provides an electric bus daytime running charge optimization method based on a Lagrange relaxation algorithm, which aims at minimizing the total charging time of a system on the premise of ensuring a preset schedule of a vehicle, obtains a full-line network level charging scheme integer optimization original problem model based on the electric quantity state of a discrete battery, converts the integer optimization original problem model into a Lagrange relaxation problem, decomposes the Lagrange relaxation problem into a plurality of subproblems, and finally converts the integral optimization original problem model into a plurality of compact subproblems.
Drawings
FIG. 1 is a schematic flow diagram of an electric bus daytime running charging optimization strategy based on a Lagrange relaxation algorithm;
FIG. 2 is a flow chart of the k-greedy algorithm of the present invention.
Detailed Description
As shown in fig. 1, the electric bus daytime running charging optimization method based on the lagrangian relaxation algorithm includes the following steps:
step S1: acquiring electric bus data and line operation data containing line quantity, wherein the electric bus data comprise the maximum charging times of each bus at each station and the minimum electric quantity protection value of the bus, the line operation data comprise route stations, arrival and departure schedules of each station and the number of charging piles of each station, a full-line network-level charging scheme integer optimization model based on the electric quantity state of a discrete battery is constructed, and an original problem model P containing coupling constraint is obtained, and specifically:
s.t.
wherein s.t. represents the constraint condition of,
andfor coupling constraint, G is an electric public transportation network, G = (N, L), N and L respectively represent the position of a bus station and the number of bus lines, L represents the L-th bus line, and L belongs to L and N l Represents the total number of stations of route l, N (1, j) represents the jth station of the route, N (1, j) is equal to N l ,a l,j Scheduled arrival time, d, at jth stop for vehicles representing route i l,j Scheduled departure time, D, for vehicles representing route l at jth stop l,j Vehicles representing route l from their jth station to jth +1Amount of electricity consumed by journey of station, D l,j-1 Representing the amount of power consumed by a vehicle on line l for its trip from the jth-1 st station to the jth station, j e { 1, …, N l -1 }, the planning period T is discrete as intervals in units of Δ tmin, T, T', T "are all time points within the planning period T, the electric bus state of charge SOC is represented as a discrete set S = {0,1,2, …, Q }, the time required for the battery to charge from the state of SOC =0 to full charge is H hours, Q = 60H/. DELTA.t, each value S ∈ S in the set represents the state of charge of the vehicle after charging from the state of SOC =0 for S time intervals Δ T, i.e. S =0 represents the amount of charge of the battery after the vehicle has charged for a time interval Δ T, S =2 represents the amount of charge of the battery after the vehicle has charged for 2 time intervals Δ T; since the charging process of the battery is nonlinear, the increment of the battery charge between each state S e S is inconsistent, i.e. the state-of-charge difference between S and S +1 is not necessarily equal to the difference between S-1 and S; in order to prolong the service life of the battery, the minimum electric quantity protection value of the vehicle of each line l is G l I.e. the vehicle SOC at the end of each journey cannot be lower than G l The maximum number of charging times allowed by the vehicles on the route l at each station is U l Considering the consumption of the battery electric quantity in the journey, the vehicle needs to be charged when passing through the station to meet the electric quantity requirement of the subsequent journey, but because the time from the vehicle to the station is different and the number of charging piles at the same station is limited, the vehicle coming first can be interrupted from charging and let the vehicle coming later and having an urgent schedule charge, the vehicle continues to be charged after the vehicle leaves, but the charging times of the vehicle are calculated every time the vehicle is plugged and pulled out, and the number of the charging piles of each station N belonging to N is M n The charging scheme is defined by a decision variable u l,j tt’ Or u l,j t’t” Is expressed when the vehicle is in time interval t, t']Or [ t', t "]When charging u is positive at jth station of line l l,j tt’ =1 or u l,j t’t” =1, station j =1, …, N for all routes of route l l Of variable v l,j,s Indicating whether the state of charge of the vehicle when leaving the jth station is s, and when the variable is 1,that is, when the state of charge is s and the variable is 0, the state of charge is not s and the variable v is l,1,Q =1 denotes the state of charge of the vehicle on route l leaving the 1 st station is Q, variable v l,j-1,s Indicating whether the state of charge of the vehicle when leaving the j-1 th station is s or not.
The original problem model P has two types of decision variables, which are respectively a charging schedule of the electric bus and an optimal battery state when the electric bus leaves a station, and the original problem model P has three types of constraints, which are respectively constraints related to a route schedule, battery electric quantity state constraints and coupling constraints.
Step S2, by introducing Lagrange multiplier mu l,j And λ n,t And mu l,j ≥0,λ n,t And ≧ 0, converting the coupling constraint condition into an objective function of the original problem P, thereby obtaining a Lagrangian relaxation problem P '(μ, λ) with a new constraint condition, wherein the Lagrangian relaxation problem P' (μ, λ) is specifically:
s.t.
step S03, apply the constant term in the above-mentioned lagrangian relaxation problem P' (μ, λ) objective function, i.e. the constant termRemoving and then decomposing the Lagrangian relaxation problem P' (mu, lambda) into | L | subproblems according to the number of lines, wherein each subproblem is respectively the decision optimization of each line charging scheme, and the first subproblem P l ' may be described as: first sub-problem P l ' described as:
s.t.
the charging scheme is determined by a decision variable u l,n tt’ Or u l,n t’t” Is expressed when the vehicle is in time interval t, t']Or [ t', t "]Then, u is charged at the nth station of the line l l,n tt’ =1 or u l,n t’t” =1, sub-problem P l Is a compound ofA variable 0 to 1 andthe strip constrained 0-1 integer programming model is extremely difficult to solve in a large-scale problem when the running schedule of each line is long in duration, so that new decision variables are continuously introduced,Represents [ t, t +1 ]]When the vehicle on line i is charging at its jth station, otherwiseWherein j =1, …, N l ,t∈{a l,j ,…,d l,j }; variables ofRepresenting the battery state of charge of the vehicle on line l as it leaves the jth station,representing the battery state of charge of the vehicle on route i when it leaves the 1 st station,representing the battery state of charge of the vehicle on line i as it leaves the j-1 st station,is a new decision variable, which when it is 1, indicates [ t, t +1 ]]Then the vehicle on line l is charging at the j-1 st station; thereby solving a plurality of the sub-questions P l ' equivalent transformation compact sub-problem P l ", compact sub-problem P l "is described as:
s.t.
step S4, solving a plurality of compact sub-problems P by using a sub-gradient descent method l ", each iteration is followed byThe Lagrange multiplier is updated by the previous solution until the algorithm reaches a termination condition, and the optimal solution of the original problem model P is obtained, and the method specifically comprises the following steps:
step S41, let the number of charging times U of the vehicle on the route l at the jth station l,j =0, number m of vehicles charged simultaneously in each station t interval n,t =0;
Step S42, traversing interval t epsilon [ a ∈ l,j +1,d l,j ]If, ifAnd is provided withWhen making U l,j + 1; traversal interval t e [ a ] l,j ,d l,j ],;
Step S43, and obtaining
The objective function of the Lagrangian relaxation problem P' (μ, λ) is then
Is equal to
Step S44, respectively using μ and λ to represent vectors formed by lagrangian multipliers in the lagrangian relaxation problem P '(μ, λ), and making L (μ, λ) represent an optimal solution corresponding to the lagrangian relaxation problem P' (μ, λ), so that the solution of the lower bound of the original problem model P can be converted into a solution of a lagrangian dual problem:
step S45, the Lagrangian dual problem P is processed LD Initializing, setting the current iteration number i =0, initializing Lagrangian multipliers mu and lambda, and setting an initial lower bound lb = - ∞ and an initial upper bound ub = + ∞; updating multipliers mu, lambda and solving compact subproblems P using a commercial solver l ", obtain the optimal solution,Calculating U by using the steps S41-S44 l,j (i) And m n,t (i) And obtaining P l "objective function value obj of LR (i) Solving the original problem by using a k-greedy algorithm and obtaining an objective function value obj corresponding to a feasible solution (i) The k-greedy algorithm solution is shown in fig. 2, the current iteration result is recorded, and if the constraint is relaxed, U is obtained l,j (i) ≤U l And m n,t (i) ≤M n If both are satisfied, the upper and lower bounds are checked and updated: if obj LR (i) If lb is greater than lb, update lower bound lb = obj LR (i) (ii) a If obj LR (i) If lb is greater than lb, update lower bound lb = obj LR (i) (ii) a Updating the Lagrange multiplier by using the secondary gradient, enabling the updating step of the multiplier to be delta =1/2(n +1), and updating the multiplier by adopting the following formula:
wherein
Judging whether the algorithm termination condition is satisfied, if all the complementary relaxation conditions are satisfied, the mu value is l,j (i+1) (U l,j (i) -U l ) =0 and λ n,t (i+1) (m n,t (i) -M n ) =0, the algorithm is terminated, and the optimal solution corresponding to the Lagrange relaxation problem is the optimal solution of the original problem at this time; when lb = ub, the algorithm is terminated, and the optimal solution corresponding to the Lagrangian relaxation problem is the optimal solution of the original problem; when the iteration number I exceeds the maximum iteration number I max When the algorithm is ended, outputting a target function and an optimal solution corresponding to an upper boundary and a lower boundary; if the three conditions are not met, the step S01-the step S05 are circulated, the iteration times are +1 in sequence, the optimal solution of the original problem is obtained until one of the three conditions is met, the daytime running charging scheme of each electric bus is obtained, and the daytime running charging of the electric buses is optimized.
The embodiment is as follows: in the embodiment, public transport network and daytime operation data (including station positions, vehicle travel schedules and the like) disclosed by a website are used, a planning time range T is divided into 600 time intervals delta T, each time interval represents 1min, the vehicle battery capacity Cap =100, and the battery minimum electric quantity protection value G l And (5). Because of the original problem P and the transformed sub-problem P l ' and compact sub-problem P l "all are described as integer programming models, and all can be solved directly by using commercial solvers (e.g., Cplex, Gurobi, etc.).
To highlight the effectiveness of the protocol of the present application, the examples compare: 1) directly solving the original problem P using Gurobi, and 2) solving the compact sub-problem P therein using steps S41-S45 l "optimal solution and solution efficiency using Gurobi for both methods, the results are shown in Table 1:
in Table 1, the value of N, L, M is calculated n 、U l Respectively represent: the total number of stations, the total number of lines and the number M of charging piles of each station N belonging to N in the network n The maximum number of charging times allowed by the vehicles on each route l at each station is U l (ii) a The objective function column represents an objective function value corresponding to the optimal solution, and the time consumed for solving is the CPU calculation time in seconds. The result shows that the scheme provided by the application can effectively solve the problem of optimizing the daytime running and charging of the electric bus, and the solving efficiency is obviously superior to that of a solution scheme directly using a Gurobi solver; and when the number of stations and the number of lines respectively reach 300/150, Gurobi cannot solve due to large problem scale, and the scheme still shows excellent solving capability, so that the scheme is more stable and efficient than the conventional method.
The above embodiments are merely illustrative of the technical ideas and features of the present invention, and the purpose thereof is to enable those skilled in the art to understand the contents of the present invention and implement the present invention, and not to limit the protection scope of the present invention. All equivalent changes and modifications made according to the spirit of the present invention should be covered in the protection scope of the present invention.
Claims (8)
1. An electric bus daytime running charging optimization method based on a Lagrange relaxation algorithm is characterized by comprising the following steps: the method comprises the following specific steps:
s1, obtaining electric bus data and line operation data containing the number of lines, constructing a full-line network-level charging scheme integer optimization model based on the electric quantity state of a discrete battery, and obtaining an original problem model P containing coupling constraint;
s2, converting the coupling constraint into an objective function of an original problem model P by introducing a Lagrange multiplier, so as to obtain a Lagrange relaxation problem P' (mu, lambda) with a new constraint condition;
s3, decomposing the Lagrangian relaxation problem P' (mu, lambda) into a plurality of sub-problems P according to the line number l ' new decision variables are introduced, i.e. whether the vehicles of the respective line are charged at the respective station, so that a plurality of said sub-questions P are presented l ' converting into multiple compact sub-problems P l ”;
S4, solving a plurality of compact sub-problems P l ", the objective function of P' (mu, lambda) is solved, that is, the objective function of the original problem model P is solved,
converting the optimal solution for solving the lower bound of the original problem model P into the solution of the Lagrangian dual problem P LD Solving the Lagrangian dual problem P by using a sub-gradient descent method LD And updating the Lagrange multiplier according to the current solution after each iteration solution,
calculating the upper-bound optimal solution of the original problem model P by adopting a k-greedy algorithm,
and obtaining the optimal solution of the original problem model P until the optimal solution of the upper bound of the original problem model P is equal to the optimal solution of the lower bound, and obtaining the optimal solution of the original problem model P, namely obtaining the daytime running charging scheme of each electric bus to complete the optimization of the daytime running charging of the electric buses.
2. The electric bus daytime running charging optimization method based on the Lagrangian relaxation algorithm as claimed in claim 1, wherein: the electric bus data comprise the maximum charging times of each bus at each station, the minimum electric quantity protection value of the bus and the electric bus electric quantity state, and the line operation data comprise the approach stations, the arrival and departure timetables of the stations and the charging pile number of the stations.
3. The electric bus daytime running charging optimization method based on the Lagrangian relaxation algorithm as claimed in claim 2, wherein the method comprises the following steps: target of original problem model P in step S1The function is specifically:
s.t.
wherein,
s.t. represents the constraint condition is;
n and L respectively represent the position and the number of bus stations, L represents the L-th bus line, L belongs to L, N represents the N-th bus station, N belongs to N, and N l Represents the total number of stations of route l;
a l,j scheduled arrival time, d, at jth stop for vehicles representing route l l,j The scheduled departure time of the vehicle representing the line l at the jth station, j epsilon { 1, …, N l -1 }, discretizing the planning time interval T into intervals in units of Δ tmin, T, T', T "each being a certain time point within the planning time interval T, the charging scheme being determined by a decision variable u l,j tt’ Or u l,j t’t” Denotes when the vehicle is in time interval t, t']Or [ t', t "]When charging u is positive at jth station of line l l,j tt’ =1 or u l,j t’t” =1;
S is a discrete set, S represents an electric bus electric quantity state SOC, S = {0,1,2, …, Q }, each value S e S in the set S represents an electric quantity state reached after the vehicle starts to be charged from the SOC =0 state for S time intervals Δ t, i.e., S =0 represents an electric quantity reached by the battery after the vehicle is charged for Δ t duration, S =2 represents an electric quantity reached by the battery after the vehicle is charged for 2 Δ t duration, time required for the battery to be fully charged from the state of SOC =0 is H hours, Q represents an electric quantity state Q, i.e., a full state, and Q =60H Δ t;
station j =1, …, N for all routes of route l l Of variable v l,j,s Whether the electric quantity state of the vehicle is s or not when the vehicle leaves the jth station is shown, when the variable is 1, the electric quantity state is shown to be s, and when the variable is 0, the electric quantity state is not shown to be s l,1,Q =1 denotes the state of charge of the vehicle on route l leaving the 1 st station is Q, variable v l,j-1,s Indicating whether the electric quantity state of the vehicle when leaving the j-1 th station is s or not;
g is an electric public transport network, G = (N, L), and the minimum electric quantity protection value of the vehicle of each line L is G l ;N l Represents the total number of stations of the route l, N (1, j) represents the jth station of the route l, N (1, j) is epsilon N l ;
D l,j Representing the amount of power consumed by a vehicle on route l for its journey from jth to j +1 th stop, D l,j-1 Represents the amount of power consumed by the vehicle on route l for its trip from the jth-1 st station to the jth station;
the maximum allowable charging times of the vehicles on the line l at each station are U l And the number of charging piles of each station N belonging to N is M n 。
5. The electric bus daytime running charging optimization method based on the Lagrangian relaxation algorithm as claimed in claim 4, wherein the method comprises the following steps: the sub-problem has | L | pieces, sub-problem P l ' is a composition containingA variable 0 to 1 and0-1 integer programming model of strip constraintType, where | S | represents the number of elements in a discrete set S, a l,n Indicating a scheduled arrival time of a vehicle at n stations on route l, d l,n Indicating the scheduled departure time of the vehicle on the line l at the n stations, the l sub-question P l ' described as:
s.t.
wherein the charging scheme is determined by a decision variable u l,n tt’ Or u l,n t’t” Denotes when the vehicle is in time interval t, t']Or [ t', t "]While being on the lineCharging u at nth station of l l,n tt’ =1 or u l,n t’t” =1。
6. The electric bus daytime running charging optimization method based on the Lagrangian relaxation algorithm as claimed in claim 5, wherein: to the sub-question P l ' Internally introducing New decision variables,Represents [ t, t +1 ]]When the vehicle on line l is charging at its jth station, and the vehicle on line l is not charging at jth station,wherein j =1, …, N l ,t∈{a l,j ,…,d l,j Will sub-problem P l ' equivalent transformation compact sub-problem P l ", compact sub-problem P l "has | L | number, the L-th compact sub-problem P l "may be described as:
s.t.
wherein, variableRepresenting the battery state of charge of the vehicle on line l as it leaves the jth station,representing the battery state of charge of the vehicle on route i when it leaves the 1 st station,representing the battery state of charge of the vehicle on line i as it leaves the j-1 st station,is a new decision variable, which when it is 1, indicates [ t, t +1 ]]The vehicle on line l is charging at its j-1 st station.
7. The electric bus daytime running charging optimization method based on the Lagrangian relaxation algorithm as claimed in claim 6, wherein: the step S4 specifically includes: step S41, let the number of times U of charging the vehicle on the route l at the jth station l,j Number of vehicles charged simultaneously m in each station t interval =0 n,t = 0; step S42, traversing interval t epsilon [ a ∈ l,j +1,d l,j ]If at allAnd isWhen making U l,j + 1; traversal interval t e [ a ] l,j ,d l,j ],(ii) a Step S43, and obtaining,N is N, T is T, the target function formula of the Lagrangian relaxation problem P' (mu, lambda) is
(ii) a The target function formula of the original problem model P is(ii) a Step S44, respectively using mu and lambda to represent vectors formed by Lagrange multipliers in the Lagrange relaxation problem P '(mu, lambda), and enabling L (mu, lambda) to represent the optimal solution corresponding to the Lagrange relaxation problem P' (mu, lambda), so that the solution of the lower bound of the original problem model P can be converted into the solution of the Lagrange dual problem P LD :Step S45, utilizing a sub-gradient algorithm to solve the Lagrangian dual problem P LD And (6) solving.
8. The electric bus daytime running charging optimization method based on the Lagrangian relaxation algorithm as claimed in claim 7, wherein the method comprises the following steps: step S45 specifically includes: step S01, the Lagrangian dual problem P is processed LD Initializing, setting the current iteration number i =0, initializing Lagrangian multipliers mu and lambda, and setting an initial lower bound lb = - ∞ and an initial upper bound ub = + ∞; step S02, updating multipliers mu and lambda and solving compact subproblems P by using a commercial solver l ", obtain the optimal solution,(ii) a Step S03, calculating U by using the steps S41-S44 l,j (i) And m n,t (i) And obtaining P l "objective function value obj of LR (i) Solving the original problem by using a k-greedy algorithm and obtaining an objective function value obj corresponding to a feasible solution (i) Recording the current iteration result, if the constraint is relaxed, i.e. U l,j (i) ≤U l And m n,t (i) ≤M n If both are satisfied, the upper and lower bounds are checked and updated: if obj LR (i) If lb is greater than lb, update lower bound lb = obj LR (i) (ii) a If obj (i) If ub is less than ub, update the upper bound ub = obj (i) (ii) a Step S04, updating the lagrangian multiplier by using the secondary gradient, making the update step of the multiplier δ =1/2(n +1), and updating the multiplier by using the following formula: mu.s l,j (i+1) =max{μ l,j (i) +δ·△μ l,j (i) ,0},λ n,t (i+1) =max{λ n,t (i) +δ·△λ n,t (i) 0} in which
Step S05, judging whether the algorithm termination condition is satisfied, if all the complementary relaxation conditions are satisfied, the μ l,j (i+1) (U l,j (i) -U l ) =0 and λ n,t (i+1) (m n,t (i) -M n ) =0, the algorithm is terminated, and at this time, the optimal solution corresponding to the lagrangian relaxation problem is the optimal solution of the original problem, and the superscripts (i) and (i +1) respectively represent the ith iteration and the (i +1) th iteration; when lb = ub, the algorithm is terminated, and the optimal solution corresponding to the Lagrangian relaxation problem is the optimal solution of the original problem; when the iteration number I exceeds the maximum iteration number I max Time, algorithmTerminating, and outputting a target function and an optimal solution corresponding to the upper and lower bounds; if the three conditions are not met, the step S01-the step S05 are repeated, the iteration times are sequentially +1, and the optimal solution of the original problem is obtained until one of the three conditions is met, namely the daytime running charging scheme of each electric bus is obtained, and the optimization of the daytime running charging of the electric buses is completed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210874670.0A CN115099702B (en) | 2022-07-25 | 2022-07-25 | Electric bus daytime running charging optimization method based on Lagrange relaxation algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210874670.0A CN115099702B (en) | 2022-07-25 | 2022-07-25 | Electric bus daytime running charging optimization method based on Lagrange relaxation algorithm |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115099702A true CN115099702A (en) | 2022-09-23 |
CN115099702B CN115099702B (en) | 2022-11-15 |
Family
ID=83298144
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210874670.0A Active CN115099702B (en) | 2022-07-25 | 2022-07-25 | Electric bus daytime running charging optimization method based on Lagrange relaxation algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115099702B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115729193A (en) * | 2022-11-24 | 2023-03-03 | 沈阳建筑大学 | Steelmaking-continuous casting tundish batch plan optimization method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113361073A (en) * | 2021-05-10 | 2021-09-07 | 上海工程技术大学 | Inventory path joint optimization method based on improved Lagrange relaxation algorithm |
CN113887032A (en) * | 2021-09-26 | 2022-01-04 | 云南电网有限责任公司 | Electric automobile ordered charging and discharging control method based on Lagrange distributed algorithm |
-
2022
- 2022-07-25 CN CN202210874670.0A patent/CN115099702B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113361073A (en) * | 2021-05-10 | 2021-09-07 | 上海工程技术大学 | Inventory path joint optimization method based on improved Lagrange relaxation algorithm |
CN113887032A (en) * | 2021-09-26 | 2022-01-04 | 云南电网有限责任公司 | Electric automobile ordered charging and discharging control method based on Lagrange distributed algorithm |
Non-Patent Citations (1)
Title |
---|
田小鹏 等: "基于拉格朗日分解的电动汽车充电路径优化", 《兰州交通大学学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115729193A (en) * | 2022-11-24 | 2023-03-03 | 沈阳建筑大学 | Steelmaking-continuous casting tundish batch plan optimization method |
Also Published As
Publication number | Publication date |
---|---|
CN115099702B (en) | 2022-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wu et al. | A novel fast-charging stations locational planning model for electric bus transit system | |
US11087422B2 (en) | Method and device for determining a configuration for deployment of a public transportation system | |
CN108955711B (en) | Navigation method applied to intelligent charging and discharging of electric automobile | |
CN103915869B (en) | A kind of Intelligent charging system of electric automobile based on mobile device and method | |
CN111311122B (en) | Charging station construction method for smart city | |
El-Taweel et al. | Integrated utility-transit model for optimal configuration of battery electric bus systems | |
CN109489676B (en) | Electric vehicle charging navigation method considering power grid information and charging station information | |
CN107274665A (en) | Bus transport capacity resource method and system for planning | |
CN114282821B (en) | Scheduling method, system and equipment for shared electric automobile | |
CN111047120B (en) | Electric automobile charging load prediction method under road-electricity coupling network | |
CN114462864B (en) | Electric bus route vehicle dispatching method under influence of charging facility sharing strategy | |
CN115099702B (en) | Electric bus daytime running charging optimization method based on Lagrange relaxation algorithm | |
Sivagnanam et al. | Minimizing energy use of mixed-fleet public transit for fixed-route service | |
CN112097783B (en) | Electric taxi charging navigation path planning method based on deep reinforcement learning | |
CN112116125A (en) | Electric vehicle charging navigation method based on deep reinforcement learning | |
He et al. | Expansion planning of electric vehicle charging stations considering the benefits of peak‐regulation frequency modulation | |
Zhou et al. | Electric bus charging scheduling on a bus network | |
Ruiz et al. | An optimal battery charging and schedule control strategy for electric bus rapid transit | |
Su et al. | Cost minimization control for electric vehicle car parks with vehicle to grid technology | |
CN113222248B (en) | Automatic taxi-driving charging pile selection method | |
Liu et al. | An MHO approach for electric bus charging scheme optimization based on energy consumption estimation | |
CN113011686A (en) | Low-carbon operation method for electric taxi fleet and battery replacement operator | |
Indira | Optimised business model for effective charging of batteries in electric vehicles | |
Binal Vansola et al. | Set cover model-based optimum location of electric vehicle charging stations | |
Solntsev et al. | Improving the Efficiency of Electric Buses in a City Transport Company |
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 | ||
CP03 | Change of name, title or address |
Address after: 35 / F, block a, Suzhou City Life Plaza, 251 pinglong Road, Gusu District, Suzhou City, Jiangsu Province 215000 Patentee after: Jiangsu Boyuxin Information Technology Co.,Ltd. Country or region after: China Address before: 35 / F, block a, Suzhou City Life Plaza, 251 pinglong Road, Gusu District, Suzhou City, Jiangsu Province 215000 Patentee before: Suzhou BOYUXIN Transportation Technology Co.,Ltd. Country or region before: China |
|
CP03 | Change of name, title or address |