CN106548632A - A kind of adjacent upstream and downstream signalized intersections queue length method of estimation - Google Patents

A kind of adjacent upstream and downstream signalized intersections queue length method of estimation Download PDF

Info

Publication number
CN106548632A
CN106548632A CN201610901122.7A CN201610901122A CN106548632A CN 106548632 A CN106548632 A CN 106548632A CN 201610901122 A CN201610901122 A CN 201610901122A CN 106548632 A CN106548632 A CN 106548632A
Authority
CN
China
Prior art keywords
downstream
phase state
queue length
signalized intersections
upstream
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
Application number
CN201610901122.7A
Other languages
Chinese (zh)
Other versions
CN106548632B (en
Inventor
唐克双
姚佳蓉
李克平
孙剑
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tongji University
Original Assignee
Tongji University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tongji University filed Critical Tongji University
Priority to CN201610901122.7A priority Critical patent/CN106548632B/en
Publication of CN106548632A publication Critical patent/CN106548632A/en
Application granted granted Critical
Publication of CN106548632B publication Critical patent/CN106548632B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/065Traffic control systems for road vehicles by counting the vehicles in a section of the road or in a parking area, i.e. comparing incoming count with outgoing count

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Traffic Control Systems (AREA)

Abstract

The present invention relates to a kind of adjacent upstream and downstream signalized intersections queue length method of estimation, it is characterised in that comprise the following steps:1) vehicle Flow Detection data and the position of fixed point detector are obtained, and is calculated within each signal period, the initial queue on Through Lane between adjacent upstream and downstream signalized intersections;2) four kinds of phase states according to adjacent upstream and downstream signalized intersections in signalized intersections signal timing dial each cycle of data acquisition;3) the initial queue length and phase state according to adjacent upstream and downstream signalized intersections, the queue length of the adjacent upstream and downstream signalized intersections in each signal period is calculated respectively, and maximum queue length and the correspondence moment of upstream and downstream signalized intersections are obtained respectively, compared with prior art, the present invention has the advantages that association measurement analysis, effectiveness height, the suitability are wide.

Description

A kind of adjacent upstream and downstream signalized intersections queue length method of estimation
Technical field
The present invention relates to traffic information field, estimates more particularly, to a kind of adjacent upstream and downstream signalized intersections queue length Method.
Background technology
Skeleton of the major urban arterial highway as transportation network, carries the sensible work(of long range rapid transit and commuting traffic Can, jam is but most frequent and serious.City letter control major trunk roads traffic congestion is not only reflected as the independence in section itself Attribute, and show influence each other between section, it is interrelated on, in the case where congestion is particularly queued up and is overflowed, section Between have obvious space-time characteristic.On the other hand, the detection technique and equipment of diversification, it is city such as to pinpoint coil, electricity police etc. The monitoring of major trunk roads traffic behavior and congestion warning provide data qualification.The coverage rate and transmission frequency of detector is to traffic shape The estimation of state has a significant impact, but the main traffic detection data source of current most of small and medium-sized cities is still single low frequency (1-min) detector is pinpointed, therefore research is intended to the testing conditions based on low precision and is relatively accurately estimated by traffic shock wave theory The queue length of signal-control crossing.
At present, the research estimated with regard to queue length both at home and abroad is mainly based upon accumulation traffic flow input/output model and friendship The method of logical ripple theoretical model.Accumulation input/output model is only applicable to situation of the queue length not less than detector location, and The model cannot obtain queue length and change over the Evolution in Spatial Dimension.Passed through based on the theoretical model of traffic shock wave The formation and dissipation of the calculating reduction vehicle queue of traffic shock wave velocity of wave, but the precision and coverage rate of the examined equipment of input of model Affect larger.In terms of data qualification, estimate mostly to employ high frequency based on the queue length of fixed point detection data in recent years Detection data (in units of s or ms), or using the detection technique of two kinds of coil, video, Floating Car etc. or more, utilize The method of multisource data fusion is to improve the accuracy of estimation, but such detection environment is not universal in reality.Therefore, build Be based on low frequency testing conditions queue length method of estimation for the queue up monitoring overflowed of signal-control crossing shows with important Sincere justice.
The content of the invention
The purpose of the present invention be exactly in order to overcome defect that above-mentioned prior art is present and provide a kind of association measurement analysis, Effectiveness is high, the adjacent upstream and downstream signalized intersections queue length method of estimation that the suitability is wide.
The purpose of the present invention can be achieved through the following technical solutions:
A kind of adjacent upstream and downstream signalized intersections queue length method of estimation, comprises the following steps:
1) vehicle Flow Detection data and the position of fixed point detector are obtained, and are calculated within each signal period, it is adjacent on The initial queue on Through Lane between downstream signal crossing;
2) four kinds according to adjacent upstream and downstream signalized intersections in signalized intersections signal timing dial each cycle of data acquisition Phase state, including first phase state:Upstream red light downstream red light, second phase state:Upstream green light downstream red light, the 3rd Phase state:Upstream green light downstream green light and the 4th phase state:Upstream red light downstream green light;
3) the initial queue length and phase state according to adjacent upstream and downstream signalized intersections, calculates each signal week respectively The queue length of the adjacent upstream and downstream signalized intersections in the phase, and the maximum queue length of stream signal crossing is obtained respectively Lj+1,m, downstream signal crossing maximum queue length Lj,m, stream signal crossing there is moment of maximum queue length tj+1,mOccurs the moment t of maximum queue length with downstream signal crossingj,m
Described step 1) in, adjacent upstream and downstream signalized intersections initial queue length QinitialFor:
QC=Σ qab*Ra+Σqab*ga
Wherein, QCFor adjacent upstream and downstream signalized intersections, qabIt is the second level stream of the turn-on interval within a signal period Amount, RaFor the duration at red light interval in a signal period, gaFor the duration at green light interval in a signal period, c is to be estimated The traffic capacity of one signal period of signalized intersections.
Described step 3) specifically include following steps:
31) the assembly ripple velocity of wave queued up in red light interval is obtained according to the second level flow of downstream signal crossing;
32) judge whether downstream signal crossing overflows according to occupation rate threshold value, if it is not, then adjacent upstream and downstream signal cross The queue length of mouth independently can be calculated, and carry out step 33), if so, then judge that the queue length of stream signal crossing is received To the impact of downstream signal crossing, and carry out step 34);
33) obtain queuing assembly ripple and evanescent wave meets the moment, the moment corresponding queue length is maximum queue length;
34) the corresponding queuing of every kind of phase state is calculated according to the phase state of different adjacent upstream and downstream signalized intersections Length, and the maximum queue length obtained in a signal period that added up by sequential.
Described step 32) in, if current downstream road section exit ramp detector occupation rate is more than or equal to occupation rate threshold value O, Then judge that downstream signal crossing is overflowed, if current downstream road section exit ramp detector occupation rate is less than occupation rate threshold valueThen Judge that downstream signal crossing is not overflowed, occupation rate threshold valueCalculating formula be:
Wherein,To pinpoint the detection meansigma methodss of detector vehicle flowrate, LvFor average effective vehicle commander, vfFor free stream velocity.
Described step 33) in, the calculating formula of the maximum queue length of upstream and downstream signalized intersections is:
Qinitial4/S01*tmax
Wherein, ω1To assemble ripple velocity of wave, tmaxAssemble ripple to queue up and evanescent wave meets the moment, seconds of the q for Through Lane Level flow, kjFor jam density, vfFor free stream velocity, ω4For the assembly ripple velocity of wave for secondary queuing occur, S0For saturation volume rate.
Described step 34) in,
When phase state is in first phase state, then downstream signal crossing is in first phase state finish time Queue length is maximum, is spacing L of adjacent upstream and downstream signalized intersectionstotal, during the spilling that downstream signal intersection vehicles are queued up Carve tsFor:
Lj,m=Ltotal
tj,m=t1+D1
Wherein, t1For the start time of first phase state, Lj,t1It is downstream intersection in first phase state start time Queue length;
Stream signal crossing is maximum in the queue length of first phase state finish time:
tj+1,m=t1+D1
Lj+1,m=Qinitial4/S0+Lj+1,t11(t2-t1)
Wherein, D1For the persistent period of first phase state, ω4For the assembly ripple velocity of wave for secondary queuing occur, S0For saturation Flow rate, Lj+1,t1For stream signal crossing first phase state start time queue length, t2Open for second phase state Begin the moment.
Described step 34) in,
When phase state is in second phase state, then row of the downstream signal crossing in the phase state finish time Team leader's degree is maximum, is spacing L of adjacent upstream and downstream signalized intersectionstotal, the spilling moment that downstream signal intersection vehicles are queued up tsFor:
Lj,m=Ltotal
tj,m=t2+D2
Wherein, t2For the start time of second phase state, Lj,t2It is downstream intersection in second phase state start time Queue length, ω1To assemble ripple velocity of wave;
When upstream crossing produces the maximum moment queued up when overflowing after the moment, then stream signal crossing is second The queue length of phase state finish time is maximum, then have:
Lj+1,m=Qinitial4/S0+Lj+1,t21*D2
tj+1,m=t2+D2
Wherein, ω4For the assembly ripple velocity of wave for secondary queuing occur, S0For saturation volume rate, D2For continuing for second phase state Time, Lj+1,t2For stream signal crossing second phase state start time queue length;
When upstream crossing produced the maximum moment queued up before the moment is overflowed, the maximum row of stream signal crossing Team leader's degree takes the maximum queuing L before overflowingj+1,m" and the queue length L of second phase state finish timej+1,m' higher value, Then have:
Lj+1,m"=Qinitial4/S01*tj+1,m
Lj+1,m'=ω4(t34-ts)+ω1(t2+D2-t34)
Lj+1,m=max { Lj+1,m″,Lj+1,m′}
Wherein, ω2For evanescent wave velocity of wave, GsBright moment, t are opened for green light34To assemble ripple and sailing out of the moment that ripple meets, ω3 To sail out of ripple velocity of wave, tsFor the spilling moment that downstream signal intersection vehicles are queued up.
Described step 34) in,
When phase state is in third phase state, then downstream signal crossing is in third phase state finish time Queue length is maximum, is spacing L of adjacent upstream and downstream signalized intersectionstotal, during the spilling that downstream signal intersection vehicles are queued up Carve tsFor:
Lj,m=Ltotal
tj,m=t3+D3
Wherein, t3For the start time of third phase state, Lj,t3It is downstream intersection in third phase state start time Queue length, D3For the persistent period of third phase state, ω1To assemble ripple velocity of wave;
The maximum queuing time t of upstream crossingj+1,mWith length Lj+1,mRespectively:
Lj+1,m2(tj+1,m-t3)-Ltotal
Wherein, Lj+1,t3For upstream crossing third phase state start time queue length, t34Hand over for stream signal Prong assemble ripple and stream signal crossing sail out of the moment that ripple meets, S0For saturation volume rate, ω2For evanescent wave velocity of wave.
If queuing spilling situation is not solved at the end of third phase state, there is secondary queuing in stream signal crossing, Then the maximum queue length of upstream crossing is:
Lj+1,m2D4+Qinitial4/S0+Lj+1,se
Lj+1,se4*Qj+1,se/S0
Qmax=Qj,t+Qj+1,t+q1*(tj+1,m-t3)
Wherein, Lj+1,seFor the secondary queue length of stream signal crossing, ω2For evanescent wave velocity of wave, ω3To sail out of ripple ripple Speed, ω4For the assembly ripple velocity of wave for secondary queuing occur, S0For saturation volume rate, D4For the persistent period of the 4th phase state, Qj+1,se The secondary queuing vehicle number for being stream signal crossing at the end of green light, QmaxIt is the maximum started counting up from downstream stop line Queuing vehicle number, Qj+1,t、Qj,tRespectively queuing vehicle number of the upstream and downstream signalized intersections in third phase state start time, t4For the start time of the 4th phase state, q1For the second level data on flows that stream signal intersection exit road detector is obtained.
Described step 34) in,
When phase state is in four phase states, if queuing up at the end of a upper phase state still overflow, under Trip signalized intersections vehicle queue in four phase states starts to reduce from crossing spacing, the maximum of downstream signal crossing Spacing of the queue length for upstream and downstream signalized intersections, if at the end of the green light of downstream signal crossing failing to dissipate completely, Remaining vehicle becomes secondary queuing, and stream signal crossing forms secondary queuing by the vehicle that a upper phase state does not dissipate, and Queue up in the 4th phase state and upstream accumulate, maximum is reached in the 4th phase state finish time queue length, then had:
Lj+1,m2D4+Qinitial4/S0+Lj+1,se
Lj,m=Ltotal
Lj,se4*Qj,se/S0
Wherein, LtotalFor the spacing of adjacent upstream and downstream signalized intersections, Qj,seTerminate in green light for downstream signal crossing When secondary queuing vehicle number, Qj,tFor S0For saturation volume rate, ω2For evanescent wave velocity of wave, ω3To sail out of ripple velocity of wave, t3For the 3rd Phase state start time, D4For the persistent period of the 4th phase state, Lj,seSecondary queuing for stream signal crossing is long Degree.
Compared with prior art, the present invention has advantages below:
First, association measurement analysis:The present invention is different from the situation that existing queuing method of estimation only considers single crossing, will Two neighboring crossing is used as holistic approach object, it is contemplated that the impact queued up to upstream crossing is overflowed in section, to crossing Between traffic relatedness carried out quantitative analysiss.
2nd, effectiveness is high:Recognition methodss and data correcting method when the present invention proposes queuing overflow detector, it is ensured that Detector is queued up in the case of occupied and to estimate the effectiveness of model.
3rd, the suitability is wide:The present invention uses the low frequency fixed point testing equipment that data transmission frequencies are 1min, to current Vehicle Detection the well adapting to property of environment of the most of small and medium-sized cities of China, universality are stronger.And under low frequency condition still Higher estimated accuracy can be so obtained, the reliability of the adjustment model is preferable.
Description of the drawings
Fig. 1 is double sections detection environment schematic in the present invention.
Fig. 2 is the emulation road network figure that adopts in the embodiment of the present invention.
Fig. 3 is the phase state schematic diagram that Adjacent Intersections are present.
Fig. 4 is traffic oscillogram of the first phase state in the case of spilling.
Fig. 5 is traffic oscillogram of the second phase state in the case of spilling.
Fig. 6 is traffic oscillogram of the third phase state in the case of spilling.
Fig. 7 is traffic oscillogram of the 4th phase state in the case of spilling.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.
Embodiment:
By adopting double sections low frequency detection method, (data transmission frequencies pinpoint testing equipment for the low frequency of 1min to the present invention Detected) queue length of the adjacent upstream and downstream crossing within each signal period is obtained, comprise the following steps:
1) vehicle Flow Detection data and the position of fixed point detector are obtained, and are calculated within each signal period, it is adjacent on The initial queue on Through Lane between downstream signal crossing, adjacent upstream and downstream signalized intersections initial queue length QinitialFor:
QC=Σ qab*Ra+Σqab*ga
Wherein, QCFor total vehicle number that signalized intersections were reached within a signal period, qabIt is a signal period The second level flow of interior turn-on interval, RaFor the duration at red light interval in a signal period, gaFor green light in a signal period The duration at interval, c are the traffic capacity of one signal period of signalized intersections to be estimated;
2) four kinds according to adjacent upstream and downstream signalized intersections in signalized intersections signal timing dial each cycle of data acquisition Phase state, including first phase state:Upstream red light downstream red light, second phase state:Upstream green light downstream red light, the 3rd Phase state:Upstream green light downstream green light and the 4th phase state:Upstream red light downstream green light, as shown in Figure 3;
3) the initial queue length and phase state according to adjacent upstream and downstream signalized intersections, calculates each signal week respectively The queue length of the adjacent upstream and downstream signalized intersections in the phase, and obtain its maximum queue length and maximum queue length occur Moment;
Comprise the following steps:
31) the assembly ripple velocity of wave queued up in red light interval is obtained according to the second level flow of upstream and downstream signalized intersections, works as row When team exceedes entrance driveway detector, then it is used to calculate using the flow of exit ramp detector and queues up, examine beyond exit ramp when queuing up Then think when surveying device that section queuing occurs overflows, the necessity of double sections detector is:Each track stream of entrance driveway detector Amount data can be used to calculate turning rate, and each track flow of exit ramp is correspondingly modified, and exclude vehicle lane change behavior Affect;
32) judge whether downstream signal crossing overflows according to occupation rate threshold value, occupation rate threshold valueCalculating formula be:
Wherein,To pinpoint the detection meansigma methodss of detector vehicle flowrate, LvFor average effective vehicle commander, vfFor free stream velocity; If it is not, then the queue length of adjacent upstream and downstream signalized intersections independently can be calculated, and carry out step 33), if so, then in judgement The queue length of trip signalized intersections is overflowed by downstream signal crossing and is affected, and carries out step 34);
33) obtain queuing assembly ripple and evanescent wave meet the moment, the moment corresponding queue length is maximum queue length, Upstream and downstream maximum queue length LmaxCalculating formula be:
Lmax=Qinitial4/S01*tmax
Wherein, ω1To assemble ripple velocity of wave, tmaxAssemble ripple to queue up and evanescent wave meets the moment, seconds of the q for Through Lane Level flow, kjFor jam density.
34) the corresponding queuing of every kind of phase state is calculated according to the phase state of different adjacent upstream and downstream signalized intersections Length, and added up by sequential and obtain the maximum queue length of Adjacent Intersections in a signal period;
As shown in figure 4, when phase state is in first phase state, then downstream signal crossing is in first phase state The queue length of finish time is maximum, is spacing L of adjacent upstream and downstream signalized intersectionstotal, downstream signal intersection vehicles row The spilling moment t of teamsFor:
Lj,m=Ltotal
Wherein, t1For the start time of first phase state, Lj,tWhen starting for first phase, the queuing of downstream intersection is long Degree;Lj,mFor the queue length that downstream intersection is maximum.
In red light phase place, not queued up to overflow by downstream intersection is affected, therefore upstream intersects for upstream intersection Mouth is queued up and accumulation can arrive first phase state finish time always, queuing of the upstream crossing in first phase state finish time Length is maximum:
tj+1,m=t1+D1
Lj+1,m=Qinitial4/S0+Lj+1,t1(t2-t1)
Wherein Di(i=1,2,3,4) represent the persistent period of i-th phase state, ω4To assemble ripple velocity of wave, with ω1Number Value is equal.
As shown in figure 5, when phase state is in second phase state, then downstream signal crossing is tied in the phase state The queue length at beam moment is maximum, is spacing L of adjacent upstream and downstream signalized intersectionstotal, the queuing of downstream signal intersection vehicles Spilling moment tsFor:
Lj,m=Ltotal
Upstream crossing overflows into " passive red light " phase place due to the queuing in downstream, and first upstream crossing is not being received The maximum moment t for queuing up is produced in the case of downstream influencesj+1,mCalculated, if tj+1,mOverflowing moment tsAfterwards, upstream Crossing is queued in tsCan start afterwards to increase, it is maximum in the queue length of second phase state finish time:
Lj+1,m=Qinitial4/S0+Lj+1,t1*D2
tj+1,m=t2+D2
If tj+1,mOverflowing moment tsBefore, i.e., upstream crossing just generated maximum queuing before the spilling of downstream, Then ripple is assembled in queue clearance to downstream and the moment t that ripple meets is sailed out of in upstream34Can start afterwards to increase, the maximum of upstream crossing Queue length takes the maximum higher value queued up with the queue length of second phase state finish time before overflowing:
Lj+1,m=Qinitial4/S01*tj+1,m
Lj+1,m'=ω4(t34-ts)+ω1(t2+D2-t34)
Lj+1,m=max { Lj+1,m,Lj+1,m'}
Wherein, GsGreen light for each signal period opens bright moment, ω3To sail out of ripple velocity of wave.
As shown in fig. 6, when phase state is in third phase state, then downstream signal crossing is tied in the phase state The queue length at beam moment is maximum, is spacing L of adjacent upstream and downstream signalized intersectionstotal, the queuing of downstream signal intersection vehicles Spilling moment tsFor:
Lj,m=Ltotal
Overflow upstream and downstream crossing after occurring and may be regarded as a section, in downstream intersection evanescent wave and upstream crossing collection The moment t that knot ripple meetsj+1,mProduce the maximum queue length of upstream crossing:
Lj+1,m2(tj+1,m-t3)-Ltotal
If queuing spilling situation is not also solved at the end of the phase state, upstream crossing occurs remaining queuing, The then secondary queue length of upstream crossing is calculated as follows:
Qmax=Qj,t+Qj+1,t+q1*(tj+1,m-t3)
Lj+1,se4*Qj+1,se/S0
Wherein, QmaxIt is the maximum queuing vehicle number started counting up from downstream stop line, Qj+1,t、Qj,tRespectively upstream and downstream is handed over Queuing vehicle number of the prong in the phase state start time, Qj+1,seThe secondary queuing for being upstream crossing at the end of green light Vehicle number, Lj+1,seThe secondary queue length for being upstream crossing at the end of green light.
As shown in fig. 7, when phase state is in four phase states, the vehicle queue of downstream intersection is between crossing Away from starting to reduce, if failing to dissipate at the end of the green light of downstream completely, remaining queue length is calculated:
Lj,m=Ltotal
Lj,se4*Qj,se/S0
The vehicle that upstream crossing is not dissipated can form secondary queuing, and the initial queue for adding the cycle is ranked length It is cumulative, at the end of the phase state, queue length is maximum:
Lj+1,se4*Qj+1,t/S0
Lj+1,m2D4+Qinitial4/S0+Lj+1,se
4) phantom that letter is drained is set up with VISSIM, precision test is carried out to queue length method of estimation.
The present invention sets up the phantom that letter is drained with VISSIM, carries out precision to queue length method of estimation and tests Card, Fig. 2 is the emulation section at Qingdao City Foochow South Road, and research section includes trunk roads, secondary distributor road, cross junction and T-shaped The road conditions such as crossing.
Modeling builds section model according to road geometry data and timing data first in VISSIM, and employing is actually obtained The detector data on flows for taking demarcates the traffic composition of section model, non-intersection speed distribution and Intersections timing scheme etc. Parameter, to ensure that the traffic behavior for emulating road network is coincide with virtual condition as much as possible.Wherein, pinpoint detector data to adopt The data of 7 thirty of morning peak to 8 thirty on November 23rd, 2015, the position of detector are as shown in Figure 1.Finally, run VISSIM Emulation road network, by the section flow for contrasting the flow for verifying each entrance driveway in crossing and fixed point detector data collection meter, really The concordance of emulation road network flow and actual road network flow is protected.
According to the time varying characteristic of actual traffic stream, basic, normal, high and segmentation four kinds of vehicles of change are provided with phantom Input type, respectively 600veh/lane/h, 1000veh/lane/h, 1200veh/lane/h and 600-1200-1200- For each vehicle input type, choose 35,40,42 3 random seeds carries out emulation testing to 600veh/lane/h., and totally 12 Group experimental group.Every time simulation time is 4200s, the simulation calculation time of preheating time and 3600s including 600s.Evaluate As a result the numerical value that 3600s calculates the time and obtains only is taken, the mean absolute difference (MAE) and average percent difference of maximum queue length is calculated (MAPE), as a result as shown in table 1.
N is number of detection cycles,For measured value, Y is value of calculation.
1 queue length estimation error statisticses table of table
From simulation results, in embodiment, the mean absolute difference of queue length is 13.58m, about two cars Length, mean percent difference are 10.97%, therefore the mean accuracy of the queue length method of estimation is more than 89%.

Claims (10)

1. a kind of adjacent upstream and downstream signalized intersections queue length method of estimation, it is characterised in that comprise the following steps:
1) vehicle Flow Detection data and the position of fixed point detector are obtained, and is calculated within each signal period, adjacent upstream and downstream The initial queue on Through Lane between signalized intersections;
2) four kinds of phase places according to adjacent upstream and downstream signalized intersections in signalized intersections signal timing dial each cycle of data acquisition State, including first phase state:Upstream red light downstream red light, second phase state:Upstream green light downstream red light, third phase State:Upstream green light downstream green light and the 4th phase state:Upstream red light downstream green light;
3) the initial queue length and phase state according to adjacent upstream and downstream signalized intersections, is calculated in each signal period respectively Adjacent upstream and downstream signalized intersections queue length, and obtain the maximum queue length L of stream signal crossing respectivelyj+1,m、 The maximum queue length L of downstream signal crossingj,m, stream signal crossing there is the moment t of maximum queue lengthj+1,mWith under There is the moment t of maximum queue length in trip signalized intersectionsj,m
2. a kind of adjacent upstream and downstream signalized intersections queue length method of estimation according to claim 1, it is characterised in that Described step 1) in, adjacent upstream and downstream signalized intersections initial queue length QinitialFor:
Q i n i t i a l = Q c - c , Q c > c 0 , Q c ≤ c
QC=∑ qab*Ra+∑qab*ga
Wherein, QCFor adjacent upstream and downstream signalized intersections, qabIt is the second level flow of the turn-on interval within a signal period, Ra For the duration at red light interval in a signal period, gaFor the duration at green light interval in a signal period, c is signal to be estimated The traffic capacity of one signal period of crossing.
3. a kind of adjacent upstream and downstream signalized intersections queue length method of estimation according to claim 2, it is characterised in that Described step 3) specifically include following steps:
31) the assembly ripple velocity of wave queued up in red light interval is obtained according to the second level flow of downstream signal crossing;
32) judge whether downstream signal crossing overflows according to occupation rate threshold value, if it is not, then adjacent upstream and downstream signalized intersections Queue length independently can be calculated, and carry out step 33), if so, then judge that the queue length of stream signal crossing is subject to down The impact of trip signalized intersections, and carry out step 34);
33) obtain queuing assembly ripple and evanescent wave meets the moment, the moment corresponding queue length is maximum queue length;
34) the corresponding length of queuing up of every kind of phase state is calculated according to the phase state of different adjacent upstream and downstream signalized intersections Degree, and the maximum queue length obtained in a signal period that added up by sequential.
4. a kind of adjacent upstream and downstream signalized intersections queue length method of estimation according to claim 3, it is characterised in that Described step 32) in, if current downstream road section exit ramp detector occupation rate is more than or equal to occupation rate threshold valueUnder then judging Trip signalized intersections overflow, if current downstream road section exit ramp detector occupation rate is less than occupation rate threshold valueThen judge downstream Signalized intersections do not overflow, occupation rate threshold valueCalculating formula be:
O ‾ = L v * q ‾ v f
Wherein,To pinpoint the detection meansigma methodss of detector vehicle flowrate, LvFor average effective vehicle commander, vfFor free stream velocity.
5. a kind of adjacent upstream and downstream signalized intersections queue length method of estimation according to claim 3, it is characterised in that Described step 33) in, the calculating formula of the maximum queue length of upstream and downstream signalized intersections is:
Qinitial4/S01*tmax
ω 1 = q - 0 k j - q / v f
Wherein, ω1To assemble ripple velocity of wave, tmaxAssemble ripple to queue up and evanescent wave meets the moment, second level streams of the q for Through Lane Amount, kjFor jam density, vfFor free stream velocity, ω4For the assembly ripple velocity of wave for secondary queuing occur, S0For saturation volume rate.
6. a kind of adjacent upstream and downstream signalized intersections queue length method of estimation according to claim 3, it is characterised in that Described step 34) in,
When phase state is in first phase state, then queuing of the downstream signal crossing in first phase state finish time Length is maximum, is spacing L of adjacent upstream and downstream signalized intersectionstotal, the spilling moment t that downstream signal intersection vehicles are queued ups For:
t s = t 1 + L t o t a l - L j , t 1 ω 1
Lj,m=Ltotal
tj,m=t1+D1
Wherein, t1For the start time of first phase state, Lj,t1For downstream intersection first phase state start time row Team leader's degree;
Stream signal crossing is maximum in the queue length of first phase state finish time:
tj+1,m=t1+D1
Lj+1,m=Qinitial4/S0+Lj+1,t11(t2-t1)
Wherein, D1For the persistent period of first phase state, ω4For the assembly ripple velocity of wave for secondary queuing occur, S0For saturated flow Rate, Lj+1,t1For stream signal crossing first phase state start time queue length, t2Start for second phase state Moment.
7. a kind of adjacent upstream and downstream signalized intersections queue length method of estimation according to claim 3, it is characterised in that Described step 34) in,
When phase state is in second phase state, then downstream signal crossing is long in the queuing of the phase state finish time Degree is maximum, is spacing L of adjacent upstream and downstream signalized intersectionstotal, the spilling moment t that downstream signal intersection vehicles are queued upsFor:
t s = t 2 + L t o t a l - L j , t 2 ω 1
Lj,m=Ltotal
tj,m=t2+D2
Wherein, t2For the start time of second phase state, Lj,t2For downstream intersection second phase state start time row Team leader's degree, ω1To assemble ripple velocity of wave;
When upstream crossing produces the maximum moment queued up when overflowing after the moment, then stream signal crossing is in second phase The queue length of state finish time is maximum, then have:
Lj+1,m=Qinitial4/S0+Lj+1,t21*D2
tj+1,m=t2+D2
Wherein, ω4For the assembly ripple velocity of wave for secondary queuing occur, S0For saturation volume rate, D2For second phase state it is lasting when Between, Lj+1,t2For stream signal crossing second phase state start time queue length;
When upstream crossing produced the maximum moment queued up before the moment is overflowed, the maximum length of queue of stream signal crossing Degree takes the maximum queuing L before overflowingj+1,m" and the queue length L of second phase state finish timej+1,m' higher value, then have:
t j + 1 , m = L j + 1 , t + ω 2 G s - ω 1 t 2 ω 2 - ω 1
Lj+1,m"=Qinitial4/S01*tj+1,m
t 34 = L j + 1 , m + ω 4 t s + ω 3 t j + 1 , m ω 3 + ω 4
Lj+1,m'=ω4(t34-ts)+ω1(t2+D2-t34)
Lj+1,m=max { Lj+1,m″,Lj+1,m′}
Wherein, ω2For evanescent wave velocity of wave, GsBright moment, t are opened for green light34To assemble ripple and sailing out of the moment that ripple meets, ω3To sail From ripple velocity of wave, tsFor the spilling moment that downstream signal intersection vehicles are queued up.
8. a kind of adjacent upstream and downstream signalized intersections queue length method of estimation according to claim 3, it is characterised in that Described step 34) in,
When phase state is in third phase state, then queuing of the downstream signal crossing in third phase state finish time Length is maximum, is spacing L of adjacent upstream and downstream signalized intersectionstotal, the spilling moment t that downstream signal intersection vehicles are queued ups For:
t s = t 3 + L t o t a l - L j , t 3 ω 1
Lj,m=Ltotal
tj,m=t3+D3
Wherein, t3For the start time of third phase state, Lj,t3For downstream intersection third phase state start time row Team leader's degree, D3For the persistent period of third phase state, ω1To assemble ripple velocity of wave;
The maximum queuing time t of upstream crossingj+1,mWith length Lj+1,mRespectively:
t j + 1 , m = L t o t a l + L j + 1 , t 3 + ω 2 * t 3 - ω 1 * t 34 ω 2 - ω 1
Lj+1,m2(tj+1,m-t3)-Ltotal
Wherein, Lj+1,t3For upstream crossing third phase state start time queue length, t34For stream signal crossing Assemble ripple and stream signal crossing sail out of the moment that ripple meets, S0For saturation volume rate, ω2For evanescent wave velocity of wave.
9. a kind of adjacent upstream and downstream signalized intersections queue length method of estimation according to claim 8, it is characterised in that If queuing spilling situation is not solved at the end of third phase state, there is secondary queuing in stream signal crossing, then upstream is handed over The maximum queue length of prong is:
Lj+1,m2D4+Qinitial4/S0+Lj+1,se
Lj+1,se4*Qj+1,se/S0
Q j + 1 , s e = Q m a x - S 0 * L t o t a l + ω 2 * t 3 + ω 3 * t 4 ω 2 + ω 3
Qmax=Qj,t+Qj+1,t+q1*(tj+1,m-t3)
Wherein, Lj+1,seFor the secondary queue length of stream signal crossing, ω2For evanescent wave velocity of wave, ω3To sail out of ripple velocity of wave, ω4For the assembly ripple velocity of wave for secondary queuing occur, S0For saturation volume rate, D4For the persistent period of the 4th phase state, Qj+1,seFor Secondary queuing vehicle number of the stream signal crossing at the end of green light, QmaxIt is the maximum row started counting up from downstream stop line Team's vehicle number, Qj+1,t、Qj,tRespectively queuing vehicle number of the upstream and downstream signalized intersections in third phase state start time, t4 For the start time of the 4th phase state, q1For the second level data on flows that stream signal intersection exit road detector is obtained.
10. a kind of adjacent upstream and downstream signalized intersections queue length method of estimation according to claim 3, its feature exist In described step 34) in,
When phase state is in four phase states, if queuing up at the end of a upper phase state still overflow, downstream letter Number crossing vehicle queue in four phase states starts to reduce from crossing spacing, the maximum queuing of downstream signal crossing Spacing of the length for upstream and downstream signalized intersections, it is if at the end of the green light of downstream signal crossing failing to dissipate completely, remaining Vehicle becomes secondary queuing, and stream signal crossing forms secondary queuing by the vehicle that a upper phase state does not dissipate, and Four phase states are queued up and are upstream accumulated, and reach maximum in the 4th phase state finish time queue length, then have:
Lj+1,m2D4+Qinitial4/S0+Lj+1,se
Lj,m=Ltotal
Q j , s e = Q j , t - S 0 * ω 2 * t 3 + ω 3 * D 4 ω 2
Lj,se4*Qj,se/S0
Wherein, LtotalFor the spacing of adjacent upstream and downstream signalized intersections, Qj,seIt is downstream signal crossing at the end of green light Secondary queuing vehicle number, Qj,tFor S0For saturation volume rate, ω2For evanescent wave velocity of wave, ω3To sail out of ripple velocity of wave, t3For third phase State start time, D4For the persistent period of the 4th phase state, Lj,seFor the secondary queue length of stream signal crossing.
CN201610901122.7A 2016-10-17 2016-10-17 A kind of adjacent upstream and downstream signalized intersections queue length estimation method Active CN106548632B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610901122.7A CN106548632B (en) 2016-10-17 2016-10-17 A kind of adjacent upstream and downstream signalized intersections queue length estimation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610901122.7A CN106548632B (en) 2016-10-17 2016-10-17 A kind of adjacent upstream and downstream signalized intersections queue length estimation method

Publications (2)

Publication Number Publication Date
CN106548632A true CN106548632A (en) 2017-03-29
CN106548632B CN106548632B (en) 2019-01-25

Family

ID=58369024

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610901122.7A Active CN106548632B (en) 2016-10-17 2016-10-17 A kind of adjacent upstream and downstream signalized intersections queue length estimation method

Country Status (1)

Country Link
CN (1) CN106548632B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108053649A (en) * 2017-12-14 2018-05-18 迈锐数据(北京)有限公司 A kind of vehicle queue length detecting system, method and apparatus
CN108765981A (en) * 2018-06-20 2018-11-06 昆明理工大学 A kind of real-time queue length prediction technique of divided lane
CN108765941A (en) * 2018-05-29 2018-11-06 重庆大学 A kind of signalized intersections vehicle arriving rate method of estimation
CN108806282A (en) * 2018-06-01 2018-11-13 浙江大学 Track group maximum queue length method of estimation based on sample travel time information
CN109166313A (en) * 2018-09-10 2019-01-08 南京市公安局交通管理局 A kind of spilling method for early warning according to car data excessively
CN110047292A (en) * 2019-05-29 2019-07-23 招商局重庆交通科研设计院有限公司 Road section congestion warning method
CN110164128A (en) * 2019-04-23 2019-08-23 银江股份有限公司 A kind of City-level intelligent transportation analogue system
EP3625786A4 (en) * 2018-08-06 2020-04-15 Beijing Didi Infinity Technology and Development Co., Ltd. Systems and methods for determining traffic conditions
CN111724600A (en) * 2020-06-30 2020-09-29 公安部交通管理科学研究所 Real-time queuing length analysis method and system based on edge calculation
CN113284336A (en) * 2021-04-04 2021-08-20 北方工业大学 Traffic state characterization method for calculating upstream and downstream intersection co-scheduling

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020077742A1 (en) * 1999-03-08 2002-06-20 Josef Mintz Method and system for mapping traffic congestion
CN102819957A (en) * 2012-06-05 2012-12-12 浙江大学 Road section queue overflow identification method based on upstream coil detectors
CN105023433A (en) * 2015-07-01 2015-11-04 重庆大学 Method for predicting range influenced by abnormal traffic event of highway

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020077742A1 (en) * 1999-03-08 2002-06-20 Josef Mintz Method and system for mapping traffic congestion
CN102819957A (en) * 2012-06-05 2012-12-12 浙江大学 Road section queue overflow identification method based on upstream coil detectors
CN105023433A (en) * 2015-07-01 2015-11-04 重庆大学 Method for predicting range influenced by abnormal traffic event of highway

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张强: "《非等周期协调信号交叉口间车辆排队模式及测算方法》", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 *
钱喆: "《过饱和交通状态下的信号控制关键技术研究》", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108053649A (en) * 2017-12-14 2018-05-18 迈锐数据(北京)有限公司 A kind of vehicle queue length detecting system, method and apparatus
CN108765941A (en) * 2018-05-29 2018-11-06 重庆大学 A kind of signalized intersections vehicle arriving rate method of estimation
CN108806282A (en) * 2018-06-01 2018-11-13 浙江大学 Track group maximum queue length method of estimation based on sample travel time information
CN108765981A (en) * 2018-06-20 2018-11-06 昆明理工大学 A kind of real-time queue length prediction technique of divided lane
EP3625786A4 (en) * 2018-08-06 2020-04-15 Beijing Didi Infinity Technology and Development Co., Ltd. Systems and methods for determining traffic conditions
CN109166313A (en) * 2018-09-10 2019-01-08 南京市公安局交通管理局 A kind of spilling method for early warning according to car data excessively
CN109166313B (en) * 2018-09-10 2021-06-25 南京市公安局交通管理局 Overflow early warning method according to vehicle passing data
CN110164128A (en) * 2019-04-23 2019-08-23 银江股份有限公司 A kind of City-level intelligent transportation analogue system
CN110164128B (en) * 2019-04-23 2020-10-27 银江股份有限公司 City-level intelligent traffic simulation system
CN110047292A (en) * 2019-05-29 2019-07-23 招商局重庆交通科研设计院有限公司 Road section congestion warning method
CN111724600A (en) * 2020-06-30 2020-09-29 公安部交通管理科学研究所 Real-time queuing length analysis method and system based on edge calculation
CN113284336A (en) * 2021-04-04 2021-08-20 北方工业大学 Traffic state characterization method for calculating upstream and downstream intersection co-scheduling

Also Published As

Publication number Publication date
CN106548632B (en) 2019-01-25

Similar Documents

Publication Publication Date Title
CN106548632B (en) A kind of adjacent upstream and downstream signalized intersections queue length estimation method
CN106652458B (en) Online urban road path forms Time Estimation Method based on virtual vehicle trajectory reconstruction
CN103680157B (en) A kind of vehicle queue's overflow pre-judging method towards city bottleneck road
Mousa Analysis and modeling of measured delays at isolated signalized intersections
CN104157152B (en) A kind of traffic signal optimization control method for intersection vehicle queue overflow status
CN106097730A (en) The method of estimation of a kind of section vehicle queue length, Apparatus and system
CN101894477B (en) Self-locking control method for controlling road network traffic through urban signal lamps
CN101739824B (en) Data fusion technology-based traffic condition estimation method
CN104778834A (en) Urban road traffic jam judging method based on vehicle GPS data
CN106710215B (en) Bottleneck upstream lane grade traffic status prediction system and implementation method
Zheng et al. A methodological framework of travel time distribution estimation for urban signalized arterial roads
CN101783073B (en) Signalized intersection delayed measuring method based on bisection detector
CN106571030A (en) Queuing length prediction method in multi-source traffic information environment
CN106530749A (en) Signal control intersection queuing length estimation method based on single section low frequency detection data
CN112885089B (en) Main line green wave intelligent diagnosis method based on multi-dimensional indexes
CN103942957A (en) Method for calculating signalized intersection vehicle queuing length under saturation condition
CN104750963A (en) Intersection delay time estimation method and device
CN107240264A (en) A kind of non-effective driving trace recognition methods of vehicle and urban road facility planing method
CN106530695A (en) Urban trunk road vehicle travel time real-time prediction method based on Internet of vehicles
CN106355882A (en) Traffic state estimation method based on in-road detector
CN106023602B (en) Mountainous City signalized intersections delay estimation method
CN106846809A (en) A kind of farthest point calculating method of queuing based on vehicle parking number of times
CN105070060B (en) A kind of urban road traffic state method of discrimination based on public transport vehicle-mounted gps data
Chepuri et al. Examining performance of an urban corridor using microscopic traffic simulation model under mixed traffic environment in India
CN107067723A (en) A kind of method of estimation of Urban road hourage

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