Background
With the continuous improvement of the urbanization level of China and the continuous increase of the automobile holding capacity, the congestion pressure of urban traffic is higher and higher. At the same time, the demand for urban traffic tends to be hierarchical and diverse. The modern tramcar as a public transportation system between the urban rapid rail transit and the traditional public transportation mode has wide development prospect in China.
However, trams have unique driving characteristics compared to social vehicles and traditional public transport. Conventional signal control methods will adversely affect their operating efficiency and safety. Therefore, on the basis of analyzing the running characteristics of the tramcar, a corresponding intersection signal control method is needed to improve the running efficiency of the tramcar.
In the field of traffic signal control, maximizing the green wave bandwidth of a road is a common coordination control strategy, and a MAXBAND algorithm is taken as a typical strategy. The passive priority strategy is widely applied to the coordination control of the main line signals of the tramcar, and the tramcar can continuously pass through a plurality of intersections in a green wave band by coordinating the phase difference of each intersection, so that the parking times and delay are effectively reduced.
On the basis of the trunk green wave, the network green wave model introduces intersection trunk crossing constraint and network outer ring closed-loop constraint, but the range of a feasible domain is further narrowed by additionally adding constraint conditions, and sometimes even no feasible solution exists.
In summary, scholars at home and abroad have obtained certain research results in the aspect of green wave algorithm and tramcar signal priority, but still have some defects: (1) at present, the coordination control method of the tramcar traffic signals is mostly researched on the basis of MAXBAND, the latest result is only the main line green wave control of the tramcar stopping on a single-line straight line, and the consideration of the tramcar line network forming is lacked. (2) The development of a network green wave control model is not mature, and particularly, an effective method is lacked to solve the problem of over-tight constraint conditions, so that a global optimal solution is difficult to obtain.
Disclosure of Invention
The purpose of the invention is as follows: the existing tramcar green wave model can only carry out green wave optimization on a single-line tramcar, and the consideration of line network formation is lacked; the existing network green wave model often has the situation that the feasible region is narrow or even is empty, and the optimal solution is difficult to obtain. Aiming at the defects of the prior art, the invention aims to provide a tramcar network green wave coordination control method and device, wherein tramcar constraints are introduced into a network green wave model, and 0-1 integer variables are used for representing whether a road green wave is interrupted. When necessary, the model breaks the green wave of partial road sections, and relaxes related constraint conditions to obtain the maximum value of the weighted sum of the network green wave bandwidths, so that the urban streetcars and social vehicles can efficiently and smoothly run, and the parking times and delay are reduced.
The technical scheme is as follows: in order to achieve the purpose, the invention adopts the following technical scheme:
a tram network green wave coordination control method comprises the following steps:
(1) acquiring road flow information related to a network and single-point timing information of each intersection, and determining road section geometric information, vehicle information and variable upper and lower limits in the network;
(2) constructing a mixed integer linear programming model for tramcar network green wave coordination control; the mixed integer linear programming model aims at maximizing the weighted sum of bidirectional green wave bandwidths of the social vehicles, and constraint conditions of the model comprise variable upper and lower limit constraints, intersection phase sequence combination constraints, social vehicle green wave constraints, tramcar green wave constraints, social vehicle and tramcar interaction constraints, intersection trunk line crossing constraints and network outer ring closed loop constraints; the method comprises the steps that 0-1 variable representing whether a road section green wave is interrupted or not is introduced into social vehicle green wave constraint, tramcar green wave constraint, social vehicle and tramcar interaction constraint and network outer ring closed loop constraint, green wave bandwidth is forced to be 0 when the green wave is interrupted, green wave time difference constraint in the social vehicle green wave constraint and tramcar green wave constraint, constraint for avoiding green wave bands from touching red light time, bidirectional minimum bandwidth constraint, social vehicle and tramcar interaction constraint and product of integer added in the network outer ring closed loop constraint and the introduced 0-1 variable representing whether the road section green wave is interrupted or not serve as penalty items;
(3) and calculating the global optimal solution of the mixed integer linear programming model by using a branch-and-bound method to obtain an optimal control scheme, wherein the optimal control scheme comprises the optimal cycle length, the phase sequence mode of each intersection and the phase difference of each intersection.
Preferably, the road flow rate ratio is taken as a bandwidth weight in the model objective function; the mark j indicates the jth road, the mark i indicates the ith intersection in the uplink direction, and
at the ith intersection, which represents the upstream direction of the trunk line j, the road flow rate ratio calculation formula is:
wherein:
and
respectively indicate intersections
And
the flow rate ratio between the up and down leg; v
i jAnd
respectively indicate intersections
And
the traffic volume between the uplink and downlink road sections; SF
i jAnd
respectively indicate intersections
And
the up and down leg saturation flow rates.
Preferably, the objective function of the mixed integer linear programming model is expressed as:
wherein: d is a trunk set; n is
jIs the number of intersections of trunk j, n
j-1 is the number of segments of trunk j;
and
are respectively an intersection
And
and the green wave bandwidth of the upstream and downstream social vehicles is reduced.
Preferably, the social vehicle green wave constraint comprises:
green wave time difference constraint:
wherein:
and
respectively in the green wave system of social vehicles
And
the time difference between the time center moments of the uplink red light and the downlink red light;
is an intersection
The time difference from the time center moment of the downlink red light to the time center moment of the nearest uplink red light;
intersection in trunk line j in green wave system of representative social vehicle
And
integer variables of road segment phase constraints;
indicating intersection
And
a binary variable of 0-1 indicating whether the green wave of the road section is interrupted or not, wherein 1 indicates the interruption and 0 indicates the non-interruption; m is a positive number not less than 1000;
is an intersection
The time difference between the right side of the uplink red light time and the center line of the green wave band of the uplink social vehicle closest to the right side of the uplink red light time;
is an intersection
The time difference between the left side of the downlink red light time and the center line of the green wave band of the downlink social vehicle closest to the left side of the downlink red light time;
and
are respectively an intersection
And
travel time of the social vehicles between ascending and descending;
and
respectively indicate intersections
Red light time in up and down directions;
to avoid the constraint set by the time when the green wave band of the variable social vehicle touches the red light:
social vehicle minimum bandwidth constraints:
wherein: bgminRepresenting a social vehicle bandwidth minimum;
when the green wave is broken, the green bandwidth is forced to 0:
preferably, the social vehicle green wave constraint further comprises:
and (3) balancing and constraining the bandwidth of the uplink green wave and the downlink green wave:
social vehicle travel time constraints:
wherein: z is the reciprocal of period C;
and
respectively indicate intersections
And
length of the section between the upper and lower directions.
Preferably, the tramcar green wave constraint includes:
green wave time difference constraint:
wherein:
and
respectively in the green wave system of the tramcar
And
the time difference between the time center moments of the uplink red light and the downlink red light;
is an intersection
The time difference from the time center moment of the downlink red light to the time center moment of the nearest uplink red light;
for a tramcar green wave system, a junction in a trunk line j
And
integer variables of road segment phase constraints;
is an intersection
The time difference between the right side of the uplink red light time and the central line of the green wave band of the uplink tramcar closest to the right side of the uplink red light time;
is an intersection
The time difference between the left side of the downlink red light time and the center line of the green wave band of the downlink tramcar closest to the left side of the downlink red light time;
and
are respectively an intersection
And
travel time of the tramcar between ascending and descending;
and
respectively indicate intersections
Red light time in up and down directions;
the tramcar green wave band can not touch the restriction of the red light:
wherein:
and
are respectively an intersection
And
the green wave bandwidth of the tramcar is increased between the ascending and descending;
is an intersection
The tramcar emptying time of (a); v. of
tminAnd v
tmaxRespectively representing the lower limit and the upper limit of the running speed of the tramcar; d represents the length of the tramcar body;
indicating intersection
The width of (d);
minimum bandwidth constraint of the tramcar:
wherein: btminRepresenting the minimum value of the bandwidth of the tramcar;
when the green wave is broken, the green bandwidth is forced to 0:
preferably, the tramcar green wave constraint further comprises:
tramcar travel time constraint:
wherein:
and
respectively indicate intersections
And
the number of tramcar stops in the ascending direction and the descending direction of the middle road section;
and
respectively indicate intersections
And
length of the section between the upper and lower directions; a represents the average acceleration of the tramcar; b represents the average deceleration of the tramcar; tau is
minRepresenting the lower limit of the stopping time of the tramcar;
and
are respectively an intersection
And
the running time of the tramcar between the ascending and the descending;
and
are respectively an intersection
And
the stopping time of the h stopping station of the tramcar between the ascending and the descending;
bidirectional bandwidth balance constraint of the tramcar:
preferably, the social vehicle and tram interaction constraint is expressed as:
wherein:
and
are respectively societyIn the green wave system of the vehicle, the intersection
And
the time difference between the time center moments of the uplink red light and the downlink red light;
and
respectively in the green wave system of the tramcar
And
the time difference between the time center moments of the uplink red light and the downlink red light; m is
iAnd
respectively represent intersections in trunk j
And
and integer variables of uplink and downlink interactive constraints between the two.
Preferably, the network outer ring closed-loop constraint is a constraint generated by the mutual correlation of phase differences of all intersections when a road provided with green waves forms a closed loop in a network structure; the method of relaxing the constraint is to add corresponding penalty terms to the left of the constraint inequality for all roads in the closed loop
Based on the same inventive concept, the tram network green wave coordination control device provided by the invention comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the computer program realizes the tram network green wave coordination control method when being loaded into the processor.
Has the advantages that: the tramcar network green wave coordination control method provided by the invention can maximize the weighted sum of the green wave bandwidth of social vehicles under the condition of meeting the tramcar passing requirement. Under necessary conditions, the model allows partial road sections not to be provided with green waves, corresponding constraint conditions are relaxed, the feasible region range is expanded, the green wave bandwidth weighted sum of the social vehicles is realized as far as possible, as many social vehicles as possible directly pass through the intersection without stopping, and the global optimal effect of urban road traffic operation is achieved.
Detailed Description
The invention is further described with reference to the following figures and specific embodiments.
The embodiment of the invention discloses a tramcar network green wave coordination control method, which mainly comprises the following steps:
and S1, acquiring the road flow information related to the network and the single-point timing information of each intersection, and determining the geometric information of the road section, the vehicle information and the upper and lower limits of the variable in the network. The method specifically comprises the following steps:
collecting road flow information: acquiring the traffic flow demand of each road by using flow acquisition equipment (including but not limited to electromagnetic induction, ultrasonic induction, video monitoring and the like), and then obtaining the flow rate ratio of each road;
determining a single-point timing scheme of each intersection in the network: calculating the key flow rate ratio of each turn according to the flow and the canalization condition of each direction of each intersection, and then obtaining the green signal ratio;
determining the values of the following input parameters of the model: the system comprises a road section length, an intersection size, upper and lower speed limits of the tramcar and the social vehicle, upper and lower cycle limits, the length of the tramcar, the acceleration and deceleration of the tramcar, the lower bandwidth limit of the tramcar and the social vehicle, the number of stop stations of the tramcar on the road section and the lower stop limit of the tramcar.
S2, constructing a mixed integer linear programming model for tramcar network green wave coordination control: wherein the objective function of the model is the weighted sum of the two-way green wave bandwidth of the maximized social vehicles; the constraint conditions include: constraint conditions of upper and lower variable limits; combining constraint conditions of phase sequences of intersections; social vehicle green wave constraints; the tramcar green wave constraint condition; social vehicles and trams are subjected to interactive constraint conditions; network intersection trunk line crossing constraint; network outer ring closed loop constraint; the method comprises the steps that 0-1 variable representing whether a road section green wave is interrupted or not is introduced into social vehicle green wave constraint, tramcar green wave constraint, social vehicle and tramcar interaction constraint and network outer ring closed loop constraint, green wave bandwidth is forced to be 0 when the green wave is interrupted, green wave time difference constraint in the social vehicle green wave constraint and tramcar green wave constraint, constraint for avoiding green wave bands from touching red light time, bidirectional minimum bandwidth constraint, social vehicle and tramcar interaction constraint and network outer ring closed loop constraint are increased by the product of an integer and the introduced 0-1 variable representing whether the road section green wave is interrupted or not are used as penalty terms;
and S3, calculating a global optimal solution of the mixed integer linear programming model for the tramcar network green wave coordination control by using a branch and bound method to obtain an optimal control scheme, wherein the optimal control scheme comprises an optimal cycle length, a phase sequence mode of each intersection and a phase difference of each intersection. And the optimized tramcar network green wave control can be realized by inputting the optimal control scheme into a signal control machine of a local traffic control bureau.
On the basis of a classic network green wave model, relevant restraint of the tramcar is introduced, and the running requirements of the tramcar and social vehicles are met; and introducing a variable 0-1 for indicating whether the green wave of the road section is interrupted or not, allowing part of the road sections not to be provided with the green wave band, and expanding the feasible region range of the model to pursue the maximum weighted sum of green wave bandwidth of the social vehicles.
Step S1 specifically includes: and after the flow of each road is obtained, calculating a road flow rate ratio according to the road traffic capacity, wherein the road flow rate ratio is used as a bandwidth weight in the model objective function. The upper mark j represents the jth road, the lower mark i represents the ith intersection in the uplink direction, and the order
At the ith intersection, representing the upstream direction of trunk j, the road flow rate ratio is calculated:
wherein:
indicating intersection
And
the flow rate ratio of the upstream (downstream) section therebetween;
and
respectively indicate intersections
And
the traffic volume between the uplink and downlink road sections;
indicating intersection
And
inter-up (down) leg saturation flow rate;
calculating the split ratio of each turn according to the phase setting of each intersection
Indicating intersection
Red light time in up (down) direction, unit: cycles; by using
Indicating intersection
Left turn green time in up (down) direction, unit: cycles;
the values of the following variables are set as input parameters for the model:
vtmin(vtmax) Represents the lower limit (upper limit) of the running speed of the tramcar, and the unit is m/s;
vgmin(vgmax) The lower limit (upper limit) of the running speed of the social vehicles at the intersection is represented by m/s;
Cmin(Cmax) Indicating green wave system signal periodLower limit (upper line), unit: s;
d represents the length of the tramcar body, and the unit is as follows: m;
a represents the average acceleration of the tramcar, and the unit is: m/s2;
b represents the average deceleration of the tramcar, and the unit is: m/s2;
btminRepresents the minimum value of the bandwidth of the tramcar, and the unit is as follows: cycles;
bgminrepresents the social vehicle bandwidth minimum, in units: cycles;
and
respectively indicate intersections
And
the number of tramcar stops in the ascending direction and the descending direction of the middle road section;
indicating intersection
The width of (d);
and
respectively indicate intersections
And
length of the section between the upper and lower directions;
τminrepresents the lower limit of the stop time of the tramcar, and the unit is as follows: s;
in step S2, the tramcar network green wave coordination control method is expressed as a mixed integer linear programming problem, and the maximum sum of the bidirectional green wave bandwidth weights of the social vehicles is taken as an objective function:
wherein:
d is a trunk set;
njis the number of intersections of trunk j, nj-1 is the number of segments of trunk j;
setting the value range of the model variable, namely an upper limit and a lower limit:
wherein:
z is the inverse of period C, and Cmin≤C≤Cmax;
Is an intersection
The right side (left side) of the up (down) red light time and the closest up (down) social vehicle greenTime difference of the center line of the wave band, unit: cycles;
is an intersection
Time difference between the right side (left side) of the up (down) red light time and the center line of the up (down) tram green wave band closest to the right side, unit: cycles;
is an intersection
And
green bandwidth of the inter-ascending (descending) social vehicle, unit: cycles;
is an intersection
And
green wave bandwidth of the ascending (descending) tramcar, unit: cycles;
is an intersection
A binary variable 0-1 representing the phase sequence setting, the value of which is 0 or 1;
is an intersection
The tramcar emptying time, unit: cycles;
is an intersection
And
travel time of the social vehicles between ascending (descending), unit: cycles;
is an intersection
And
travel time of the ascending (descending) tramcar, unit: cycles;
is an intersection
And
the unit of the running time of the ascending (descending) tramcar is as follows: cycles;
is an intersection
And
the stopping time of the h stopping station of the ascending (descending) tramcar between the first and the second stop stations, the unit: cycles;
indicating intersection
And
a binary variable of 0-1 (1: interrupted, 0: not interrupted) whether the green wave of the middle road section is interrupted or not;
inf represents positive infinity.
Adding intersection phase sequence combination constraints:
wherein:
is an intersection
The time difference from the time center moment of the downlink red light to the time center moment of the nearest uplink red light is positive when the time center moment of the uplink red light is positioned on the right side of the time center moment of the downlink red light, and negative otherwise;
the following social vehicle green wave constraints were added:
green wave time difference constraint:
wherein:
and
respectively in the green wave system of social vehicles
And
the red light time is selected according to the time difference between the time center moments of the uplink red light and the time center moments of the downlink red light, and the red light time is met by the left side (right side) of the green wave band of the same uplink (downlink) social vehicle
First up (down) red light time of and
first up (down) red light time, unit: cycles;
intersection in trunk line j in green wave system of representative social vehicle
And
inter-road phase constraintAn integer variable of (a);
in order to avoid the time when the green wave band of the variable social vehicle touches the red light, the following constraints are set:
the minimum bandwidth constraint of the social vehicles, the bidirectional bandwidth of the social vehicles should not be less than the minimum bandwidth requirement:
when the green wave is interrupted, the width of the green wave band is forced to be 0
In order to ensure that the green wave bandwidth of the uplink and downlink social vehicles can reflect the traffic demands of the social vehicles in two directions to a certain extent, the uplink and downlink green wave bandwidth balance constraint is added:
social vehicle travel time constraints:
the following tram green wave constraints were added:
green wave time difference constraint:
wherein:
and
respectively in the green wave system of the tramcar
And
the red light time is selected according to the time difference between the time center moments of the uplink red light and the time center moment of the downlink red light, and the red light time is met by the left side (right side) of the green wave band of the same uplink (downlink) tramcar
First up (down) red light time of and
first up (down) red light time, unit: cycles;
for a tramcar green wave system, a junction in a trunk line j
And
integer variables of road segment phase constraints;
the tramcar can conflict with the left-turn traffic flow in the direction, so that the green wave band of the tramcar can not touch the red light on the upper line nor the red light on the lower line. In addition, a certain emptying time is reserved for the tramcar, so the following constraints are set:
the tramcar travel time is the sum of the travel time and the station stop time:
time of flight
The following constraints can be satisfied:
the tram stop time lower limit needs satisfy the passenger demand of going up and down:
the bidirectional bandwidth balance constraint of the tramcar is characterized in that in order to ensure that the service levels of the tramcar in the uplink direction and the tramcar in the downlink direction are the same, the uplink bandwidth and the downlink bandwidth are also the same:
the minimum bandwidth constraint of the tramcar, the bidirectional bandwidth of the tramcar should not be less than the lowest bandwidth requirement:
when the green wave is broken, the green bandwidth is forced to 0:
social vehicle and tram interaction constraint is added,
and
the time value of the phase difference is an integral multiple relation of the period, namely:
wherein:
m
iand
respectively represent intersections in trunk j
And
and integer variables of uplink and downlink interactive constraints between the two.
Adding network intersection constraints:
wherein:
if there are two crossed main lines, the numbers are alpha and beta, and the crossing point is the ith direction of the main line alpha
αEach signal control intersection
Is the i-th direction of the main line beta
βEach signal control intersection
Then
And
respectively indicate intersections
Left turn green time in up and down direction;
and
respectively indicate intersections
Left turn green time in up and down direction;
is an intersection
A binary variable 0-1 representing the phase sequence setting, the value of which is 0 or 1;
is an intersection
A binary variable 0-1 representing the phase sequence setting, the value of which is 0 or 1;
is a crossing
The time center of the ascending red light arrives at the intersection
Time difference of the time center moment of the up red light. When crossing
The time center of the ascending red light is positioned at the intersection
When the time center of the ascending red light is right, the value is positive, otherwise, the value is negative;
adding network outer ring closed-loop constraint, namely when a road provided with green waves forms a closed loop in a network structure, the constraint generated by the mutual correlation of the phase differences of all intersections needs to be constructed respectively aiming at a specific network structure, and the method for relaxing the constraint is that for all roads in the closed loop, corresponding punishment items are added on the left side of a constraint inequality. For the network structure shown in fig. 2, the closed loop constraint at the outer ring of the network is as follows:
wherein:
m1,2,3,4the integer variables corresponding to the outer-loop constraints composed of the main lines 1,2,3 and 4 in fig. 2. For an outer ring closed loop consisting of the main lines α,.,. beta, a similar relaxed constraint can be constructed as described above. Integer variable m in network outer ring closed loop constraints1,2,3,4Will be composed ofα,...,βAnd (4) substitution.
Based on the same inventive concept, the tram network green wave coordination control device provided by the embodiment of the invention comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, and the computer program realizes the tram network green wave coordination control method when being loaded to the processor.
The following describes the scheme of the embodiment of the present invention in detail with reference to the specific scenario shown in fig. 3. In this scenario, the physical conditions of the road network, and the traffic information of the intersection approach and road segment are all known. The first road network has 8 intersections, and each intersection has 4 entrance roads. Fig. 3 shows only the key roads in the road network, for a total of 9, and the other roads are not shown. The flow rate of each road is 1000pcu/h, and the traffic capacity is 2000 pcu/h. The proportion of left turn and straight going of each inlet passage is 1: 4.
according to the traffic flow demand of each road, the flow rate ratio of each road is obtained, and the flow rate ratio is shown in a table 1;
TABLE 1 Main road traffic, traffic capacity and flow rate ratio of road network
Calculating the key flow rate ratio of each turn according to the flow and the canalization condition of each direction of each intersection, and then obtaining the green signal ratio; the steering flow of the intersection entrance road is shown in a table 2; the proportion of the left turn green light in the trunk direction and the proportion of the straight red light in the trunk direction at the intersection obtained by the method are shown in a table 3.
TABLE 2 intersection node entry lane steering flow
TABLE 3 proportion of the cycle of the left turn green light and the straight red light at the intersection
Determining the values of the following input parameters of the model: the system comprises a road section length, an intersection size, upper and lower speed limits of the tramcar and the social vehicle, upper and lower cycle limits, the length of the tramcar, the acceleration and deceleration of the tramcar, the lower bandwidth limit of the tramcar and the social vehicle, the number of stop stations of the tramcar on the road section and the lower stop limit of the tramcar.
The length of each road segment from fig. 2 is 500 meters except for the length of the h-b road segment of 1000 meters. The width of each intersection is 30 m; the upper limit and the lower limit of the speed of the social vehicle are both 16.7m/s, and the upper limit and the lower limit of the speed of the tramcar are both 10.75 m/s; the upper limit of the period is 150s, and the lower limit is 90 s; the length of the tramcar is 15m, and the acceleration and deceleration of the tramcar are both 3m/s2(ii) a The lower limit of the bandwidth of the tramcar is 0.05cycles, and the lower limit of the bandwidth of the social vehicle is 0.15 cycles; the number of tramcar stop stations on each road is 1; the lower limit of the stopping time of the tramcar is 5 s.
And constructing the mixed integer linear programming model for the tramcar network green wave coordination control according to the information, and calculating the global optimal solution of the mixed integer linear programming model for the tramcar network green wave coordination control by using a branch-and-bound method. And the finally obtained target function global maximum value is 3 cycles.
Table 4 shows the signal control scheme for each intersection in the model output result. The phase difference between adjacent intersections and the green bandwidth optimization results are shown in table 5.
TABLE 4 signalization scheme at each intersection
TABLE 5 phase difference and Green Bandwidth between Adjacent intersections
As can be seen from table 5, when the sum of the green bandwidth weights of the social vehicles is maximized, only 5 roads are provided with green bands, and 4 roads are not provided with green bands. The space-time diagram is plotted according to the solution of the model, as shown in fig. 4.
Intersections and green wave bands of the main lines {1,2,3,4,5} in the upward (downward) direction are drawn in this order from bottom to top (top to bottom). The red light time of the up (down) straight phase at each intersection is indicated by the horizontal solid line (dashed line). The intersection non-trunk direction phase duration is indicated by a bold solid line. Different intersections are marked with different colors. The social vehicle green wave is indicated by a light oblique line, and the tramcar green wave is indicated by a dark oblique line.