CN107038862A - A kind of traffic assignation modification method under the conditions of tail number restricted driving - Google Patents
A kind of traffic assignation modification method under the conditions of tail number restricted driving Download PDFInfo
- Publication number
- CN107038862A CN107038862A CN201710306621.6A CN201710306621A CN107038862A CN 107038862 A CN107038862 A CN 107038862A CN 201710306621 A CN201710306621 A CN 201710306621A CN 107038862 A CN107038862 A CN 107038862A
- Authority
- CN
- China
- Prior art keywords
- traffic
- tail number
- restricted driving
- distribution
- assignation
- 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
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0137—Measuring and analyzing of parameters relative to traffic conditions for specific applications
Abstract
A kind of traffic assignation modification method under the conditions of being restricted driving the invention discloses tail number.This method include obtaining the basic road network information in newest city, obtain initial right of way table for traffic assignation, obtain traffic zone information for traffic assignation, obtain the newest each trip mode OD distribution matrixs in city, obtain motor vehicle tail number restricted driving scheme for test, determine twice the ratio of OD distribution matrixs used in traffic assignation, carry out eight steps of automobile traffic amount that each section is distributed under multipath capacity limit iteration traffic assignation, output tail number restricted driving scheme.The present invention has taken into full account the existing motor vehicle tail number restricted driving scheme of China, foundation is used as using restricted motor vehicle tail number number, the trip mode OD distribution matrixs limited by tail number are split as unrestricted tail number part and limited tail number part, and carry out traffic assignation respectively.The inventive method can be used in the virtual traffic test platform of city testing the quantization of motor vehicle tail number restricted driving strategy implement effect.
Description
Technical field：
A kind of traffic assignation modification method under the conditions of being restricted driving the present invention relates to tail number, belongs to urban road network traffic point
With technical field.
Background technology：
Motor vehicle tail number restricted driving measure is a kind of distinctive traffic system of China, originally controls plan as a kind of traffic administration
Purpose when being slightly born is to alleviate urban traffic blocking, and Beijing implements during being earliest Beijing Olympic Games test event in 2007
Motor vehicle tail number odd-and-even license plate rule.Recently as the swift and violent growth of the major Urban vehicles poputations of China, transport need by
Gradually close to even greater than transportation supplies, most intuitively performance is exactly the traffic congestion of urban road, and many cities are private except putting into effect
Outside family's car limit purchase policy, it is also desirable to more directly alleviate traffic congestion by tail number restricted driving measure, and tail number restricted driving measure has
The estimation of effect property is essential link.
Urban road network method of traffic assignment international at present is very ripe, but distinctive for China
Traffic assignation modification method under the conditions of the restricted driving of motor vehicle tail number, or research blank.Driver under the conditions of being restricted driving due to tail number
Path selection result and be not carried out during tail number restricted driving measure be it is completely different, while tail number restrict driving after point mode OD distribution
Respective change can occur for matrix, and still carrying out tail number restricted driving measure validity test using the method for traffic assignment not being corrected is
It is inappropriate, therefore more than ten cities of tail number restricted driving measure were had been carried out at present, it is only capable of implementing tail number restricted driving by actual
Link counting situation after scheme carries out scheme evaluation, and then tail number restricted driving scheme is adjusted accordingly, and can not shift to an earlier date
Validity test is carried out to the tail number restricted driving scheme to be implemented.
Under smart city, wisdom traffic construction background, many cities will all set up city virtual traffic system testing and put down
Platform, and the test of tail number restricted driving scheme will turn into urban traffic control control program and quantify one of function of test platform.According to
The basic road network information in newest city, traffic zone information, point mode OD distribution matrixs data and the tail number for test are restricted driving
Objective for implementation, action time and the sphere of action of scheme, are modified to the idiographic flow of OD distribution matrixs and traffic assignation, obtain
Traffic distribution result under the conditions of being restricted driving to tail number, you can quantization is carried out to the tail number restricted driving scheme implementation result tested and commented
Estimate, it is to avoid directly implement inappropriate tail number restricted driving scheme and brought inconvenience to the traffic trip of city dweller.Therefore, of the invention
Propose the traffic assignation modification method under the conditions of a kind of tail number is restricted driving.
The content of the invention
The purpose of the present invention is the Path selection result of driver under the conditions of being restricted driving for existing tail number with tail number is not carried out
It is completely different during restricted driving measure, but asking there is presently no the traffic assignation modification method under the conditions of being restricted driving for tail number
Topic, proposes the traffic assignation modification method under the conditions of a kind of tail number restricted driving, has taken into full account that each big city of China is existing motor-driven
Tailstock restricted driving scheme, is modified to method of traffic assignment, compensate for traffic assignation modification method under the conditions of tail number is restricted driving
Blank, is surveyed to accurate in the virtual traffic test platform of city, the quick quantization for realizing motor vehicle tail number restricted driving strategy implement effect
Examination has greater significance.
Above-mentioned purpose is achieved through the following technical solutions：
A kind of traffic assignation modification method under the conditions of tail number restricted driving, this method comprises the following steps：
A kind of traffic assignation modification method under the conditions of tail number restricted driving, this method comprises the following steps：
(H) the basic road network information in newest city, including intersection numbering Node_ID (dimensionless), section numbering are obtained
Link_ID (dimensionless), road section length Ll(unit is m) and section practical capacity Cl(unit is pcu/h)；
(I) the initial right of way table W for traffic assignation is obtained, its method is each section numbering in first obtaining step (A)
Desin speed S corresponding to Link_IDD(unit is km/h) and Maximum speed limit SL(unit is km/h), then using formula (1)
Calculating obtains the free flow journey time T corresponding to each section numbering Link_IDl0(unit is s), these sections numbering Link_
Journey time corresponding to ID is initial right of way table W：
Tl0=3600 × { 0.001 × Ll/min{SD,SL}} (1)；
(J) the traffic zone information distributed for urban transportation, the numbering Cell_ID comprising each traffic zone are obtained
Node serial number Node_ID and section numbering Link_ID that (dimensionless), each traffic zone are included；
(K) the newest each trip mode OD distribution matrixs in city are obtained；
(L) obtain the region motor vehicle tail number restricted driving scheme for test, comprising in the daily restricted driving time in working day by
The motor vehicle tail number number N of limitationL(unit is individual, value is 1,2 ..., step (D) 9), corresponding to tail number restricted driving scheme
In trip mode (trip mode restricted driving by tail number), traffic zone in the completely covered step (C) in tail number restricted driving region
The numbering Cell_ID of (by tail number restricted driving cell)；
Or：
The urban parts major trunk roads motor vehicle tail number restricted driving scheme for test is obtained, during comprising being restricted driving daily in working day
Interior restricted motor vehicle tail number number NL(unit is individual, value is 1,2 ..., 9), corresponding to tail number restricted driving scheme
The numbering Link_ by tail number restricted driving section that the trip mode restricted driving by tail number, tail number restricted driving scheme in step (D) are covered
ID”；
(M) it is restricted motor vehicle tail number number N according to what is determined in step (E)LIt is determined that dividing for the first time with second of traffic
Ratio P with OD distribution matrixs used respectively1And P2, value rule such as formula (2)：
(N) multipath-capacity restraint traffic assignment method under the conditions of being restricted driving using tail number is iterated traffic assignation；
(I) the automobile traffic amount V that each section is distributed is exportedl。
Traffic assignation modification method under the conditions of a kind of described tail number restricted driving, the city described in step (D) is respectively gone on a journey
Mode OD distribution matrixs, including：Private minibus trip mode OD distribution matrixs CarOD (unit is pcu/h), regular public traffic go out
Line mode OD distribution matrixs BusOD (unit is pcu/h), track traffic trip mode OD distribution matrixs RailOD (unit behave/
H), motorcycle trip mode OD distribution matrixs MotorOD (unit is pcu/h), trip of taxi mode OD distribution matrixs
TaxiOD (unit is pcu/h), non-motor vehicle trip mode OD distribution matrixs BikeOD (unit behaviour/h) and walking trip side
It is several in formula OD distribution matrixs WalkOD (unit behaviour/h).
It will be restricted driving in traffic assignation modification method under the conditions of a kind of described tail number restricted driving, the step (F) by tail number
Trip mode OD distribution matrixs be split as two parts of unrestricted tail number part and limited tail number part, and carry out traffic twice point
Match somebody with somebody, the unrestricted tail number part OD distribution matrixs of the first sub-distribution, the second sub-distribution is limited tail number OD distribution matrixs.
The tail number limit used in traffic assignation modification method under the conditions of a kind of described tail number restricted driving, the step (G)
Multipath-capacity restraint traffic assignment method under the conditions of row, its iteration traffic assignation process is comprised the steps of：
(G12) traffic assignation iterations is designated as K, and the OD amounts ratio of the secondary iteration loading of kth (k=1,2 ... ... K) is Qk,
AndMake k=1；
(G13) the OD distribution matrixs of each trip mode in step (D) are multiplied by Q respectivelyk, obtain being used for kth time iteration friendship
The OD distribution matrixs that the reduction of fractions to a common denominator is matched somebody with somebody, start kth time iteration；
(G14) for the trip mode restricted driving by tail number obtained in step (E), obtained in the step corresponded to (G2)
OD distribution matrixs be multiplied by P1The OD distribution matrixs (unrestricted tail number part) for first time traffic assignation are obtained, are corresponded to
The step of (G2) in obtained OD distribution matrixs be multiplied by P2Carry out the transfer amendment of OD distribution matrixs again afterwards, that is, obtain being used for second
The OD distribution matrixs (limited tail number part) of secondary traffic assignation；
Or：
For the trip mode restricted driving by tail number obtained in step (E), the OD obtained in the step corresponded to (G2)
Distribution matrix is multiplied by P1The OD distribution matrixs (unrestricted tail number part) for first time traffic assignation are obtained, by its corresponding step
Suddenly the OD distribution matrixs obtained in (G2) are multiplied by P2Obtain OD distribution matrixs (the limited tail number portion for second of traffic assignation
Point)；
(G15) for the trip mode do not restricted driving in step (D) by tail number, shifted using OD distribution matrixs in step (G3)
The diverted traffic produced during amendment updates the trip mode OD distribution matrixs do not restricted driving by tail number, with reference to right of way table W and renewal
The trip mode OD distribution matrixs do not restricted driving afterwards by tail number carry out kth time iteration traffic assignation, and the volume of traffic after distribution is tired out
Meter is loaded into corresponding road section automobile traffic amount VlOn (unit is pcu/h)；
(G16) according to each section automobile traffic amount V obtained in step (G4)l, using Bureau of Public Road's roadlock
Function model updates the journey time T in each section in right of way table Wl(unit is s), i.e. formula (3)：
Tl=Tl0[1+α(Vl/Cl)β] (3)
Wherein α=0.15, β=4；
(G17) after using the unrestricted tail number part OD distribution matrixs obtained in step (G3) and step (G5) updated
Right of way table W carries out first time traffic assignation, and is loaded into corresponding road section automobile traffic amount V by the volume of traffic after distribution is accumulativel
On；
(G18) according to each section automobile traffic amount V obtained in step (G6)l, right of way table W is updated using formula (3)
In each section journey time Tl；
(G19) by by the traffic zone information matches obtained in the numbering Cell_ID of tail number restricted driving cell and step (C),
Obtain the section numbering Link_ID in all sections that these traffic zones are included, and by the right of way table W obtained in step (G7)
In these sections number corresponding journey time TlValue be changed to any number not less than 99999, obtain tail number restricted driving bar
Right of way table W under partL；
Or：
By in the right of way table W obtained in step (G7) by the corresponding journey time T of tail number restricted driving section numbering Link_IDl
Value be changed to any number not less than 99999, obtain the right of way table W under the conditions of tail number is restricted drivingL；
(G20) using the tail number obtained in the limited tail number part OD distribution matrixs and step (G8) obtained in step (G3)
Right of way table W under the conditions of restricted drivingLCarry out second of traffic assignation, and be loaded into corresponding road section machine by the volume of traffic after distribution is accumulative
Motor-car volume of traffic VlOn；
(G21) according to each section automobile traffic amount V obtained in step (G9)l, right of way table W is updated using formula (3)
In each section journey time Tl；
(G22) make k=k+1, if k≤K, return to step (G2), otherwise iteration traffic assignation terminate.
OD moments of distribution are carried out in traffic assignation modification method under the conditions of a kind of described tail number restricted driving, the step (G3)
The process of battle array transfer amendment is comprised the steps of：
(G3d) the distribution volume of traffic that note traffic occurs between cell i and road traffic simulation cell j is Qij(i ∈ { Cell_ID },
J ∈ { Cell_ID }), remember QijTransfer amount be QTij, for all friendships corresponding to the cell number Cell_ID that is restricted driving by tail number
Logical generating capacity and traffic attraction, Schilling QTij=Qij, then make Qij=0；
(G3e) for remaining Q in OD distribution matrixsij, calculate the rate γ that detours under the conditions of tail number is restricted drivingij；
(G3f) OD amount transfer functions Q is combinedTij=F (γij,Qij), calculate remaining QijTransfer amount QTij, then make Qij=
Qij-QTij, that is, complete the transfer amendment of OD distribution matrixs.
Tail number limit is calculated in traffic assignation modification method under the conditions of a kind of described tail number restricted driving, the step (G3b)
The rate γ that detours under the conditions of rowijThe step of it is as follows：
(G3b1) by by the traffic obtained in the numbering Cell_ID of tail number restricted driving cell and step (C) described in claim 1
Cell information is matched, and obtains the section numbering Link_ID in all sections that these traffic zones are included, and by claim 4
Corresponding journey time T is numbered in these sections in newest right of way table W under middle current procedureslValue be changed to be not less than 99999
Any number, obtain tail number restrict driving under the conditions of right of way table WL；
(G3b2) the right of way table W under the conditions of being restricted driving respectively according to right of way table W and tail numberLOD is calculated to the shortest path between ij
Weigh DijAnd DLij；
(G3b3) the rate γ that detours is calculated using formula (4)ij：
γij=DLij/Dij (4)。
Traffic assignation modification method under the conditions of a kind of described tail number restricted driving, OD amounts described in the step (G3c) turn
Move function QTij=F (γij,Qij) expression formula be formula (5)：
Wherein β1、β0For coefficient to be calibrated, γ0To specify all volume of traffic under confidence level not change trip mode
Maximum transfer rate, γ1The minimum rate of transform of trip mode, γ are changed for all volume of traffic under specified confidence level0And γ1Can
Directly obtained by survey, β can be demarcated using formula (6)1And β0：
Updated in traffic assignation modification method under the conditions of a kind of described tail number restricted driving, the step (G4) not by tail number
The method of the trip mode OD distribution matrixs of restricted driving is：By obtained all diverted traffic QTijAccording to what is do not restricted driving by tail number
Same OD is split to an ij person-time volume of traffic ratio in trip mode OD distribution matrixs and each superposition is entered, and is reconverted into
Respective original unit's volume of traffic is the trip mode OD distribution matrixs do not restricted driving after being updated by tail number.
Traffic assignation modification method under the conditions of a kind of described tail number restricted driving, private minibus, regular public traffic, motorcycle
The conversion coefficient respectively θ mutually changed with original unit's volume of traffic with person-time volume of traffic of taxi these four trip modescar=
1.3、θbus=40.0, θmotor=1.0, θtaxi=1.3, original unit's volume of traffic is multiplied by conversion coefficient and obtains a person-time volume of traffic, people
The secondary volume of traffic divided by conversion coefficient are to be changed into original unit's volume of traffic；These three trip modes of track traffic, non-motor vehicle and walking
A person-time volume of traffic be original unit's volume of traffic, without conversion.
Beneficial effect：
A kind of tail number proposed by the present invention restrict driving under the conditions of traffic assignation modification method, by by going out for being limited by tail number
Line mode OD distribution matrixs are split as unrestricted tail number part and limited tail number part, and carry out traffic assignation respectively, by region
Or the tail number restricted driving scheme of urban major trunk roads is integrated well in terms of OD distribution matrixs, right of way table and assigning process three
Into the related procedure of traffic assignation, both ensure that method of traffic assignment is easy to computational, but conscientiously, consider tail exactly
Effective object, action time and the action space of number restricted driving scheme.The present invention has taken into full account China's existing machine in each big city
Motor-car tail number restricted driving scheme, is modified to method of traffic assignment, compensate for traffic assignation modification method under the conditions of tail number is restricted driving
Blank, to accurate in the virtual traffic test platform of city, the quick quantization for realizing motor vehicle tail number restricted driving strategy implement effect
Test has greater significance.
Brief description of the drawings
Fig. 1 is the overview flow chart of the inventive method.
Embodiment
With reference to embodiment, the present invention is furture elucidated, it should be understood that following embodiments are only used for
The bright present invention rather than limitation the scope of the present invention.
Embodiment 1：
Method described in the present embodiment can realize that motor vehicle tail number restricted driving scheme in region is implemented, and comprise the concrete steps that：
It is as shown in Figure 1 the overview flow chart of the present invention, is further illustrated with reference to embodiment.
(A) the basic road network information in newest city, including intersection numbering Node_ID (dimensionless), section numbering are obtained
Link_ID (dimensionless), road section length Ll(unit is m) and section practical capacity Cl(unit is pcu/h)；
Embodiment chooses a simplified city, and its basic road network information is as follows：
Intersection information:Node_ID=1,2 ... ..., 24；
Road section information：It is shown in Table 1.
Acquired road section information in the embodiment of table 1
(B) the initial right of way table W for traffic assignation is obtained, its method is each section numbering in first obtaining step (A)
Desin speed S corresponding to Link_IDD(unit is km/h) and Maximum speed limit SL(unit is km/h), then using formula (1)
Calculating obtains the free flow journey time T corresponding to each section numbering Link_IDl0(unit is s), these sections numbering Link_
Journey time corresponding to ID is initial right of way table W.
The free flow journey time T obtained by data and use formula (1) calculating in the present embodiment needed for this stepl0It is shown in Table
1, then the initial right of way table W of embodiment be shown in Table 2.
The initial right of way table W of the present embodiment of table 2
| Section numbering Link_ID | Right of way (s) |
| 1 | 28 |
| 2 | 18 |
| 3 | 14 |
| 4 | 24 |
| …… | …… |
| 8 | 33 |
| …… | …… |
| 11 | 46 |
| …… | …… |
| 76 | 34 |
(C) the traffic zone information distributed for urban transportation, the numbering Cell_ID comprising each traffic zone are obtained
Node serial number Node_ID and section numbering Link_ID that (dimensionless), each traffic zone are included；
Traffic zone information in the present embodiment is shown in Table 3.
The traffic zone information of the present embodiment of table 3
(D) the newest each trip mode OD distribution matrixs in city are obtained；
Three kinds of Passenger Traveling Choice individual minibus, regular public traffic and non-motor vehicle in the present embodiment, newest city is each
Trip mode OD distribution matrixs are shown in Table 4-6.
CarOD (pcu/h) in the present embodiment of table 4
| Cell_ID | 1 | 2 | 3 | 4 | 5 | Traffic generating capacity adds up to |
| 1 | 3000 | 1000 | 5000 | 1000 | 1000 | 11000 |
| 2 | 1000 | 3000 | 5000 | 1000 | 1000 | 11000 |
| 3 | 5000 | 5000 | 1000 | 5000 | 5000 | 21000 |
| 4 | 1000 | 1000 | 5000 | 3000 | 1000 | 11000 |
| 5 | 1000 | 1000 | 5000 | 1000 | 3000 | 11000 |
| Road traffic simulation amount adds up to | 11000 | 11000 | 21000 | 11000 | 10500 | 65000 |
BusOD (pcu/h) in the present embodiment of table 5
| Cell_ID | 1 | 2 | 3 | 4 | 5 | Traffic generating capacity adds up to |
| 1 | 70 | 30 | 70 | 30 | 20 | 220 |
| 2 | 30 | 70 | 70 | 20 | 30 | 220 |
| 3 | 70 | 70 | 100 | 70 | 70 | 380 |
| 4 | 30 | 20 | 70 | 70 | 30 | 220 |
| 5 | 20 | 30 | 70 | 30 | 70 | 220 |
| Road traffic simulation amount adds up to | 220 | 220 | 380 | 220 | 220 | 1260 |
BikeOD (people/h) in the present embodiment of table 6
(E) obtain the region motor vehicle tail number restricted driving scheme for test, comprising in the daily restricted driving time in working day by
The motor vehicle tail number number N of limitationL(unit is individual, value is 1,2 ..., step (D) 9), corresponding to tail number restricted driving scheme
In trip mode (trip mode restricted driving by tail number), traffic zone in the completely covered step (C) in tail number restricted driving region
The numbering Cell_ID of (by tail number restricted driving cell)；
Region motor vehicle tail number restricted driving scheme in the present embodiment is traffic zone 3 on weekdays 07:00~09:00、17:
00~19:00 for 2 tail numbers of private minibus restricted driving, i.e. NL=2, the trip mode restricted driving by tail number be private minibus,
It is 3 by tail number restricted driving cell number Cell_ID.
(F) it is restricted motor vehicle tail number number N according to what is determined in step (E)LIt is determined that dividing for the first time with second of traffic
Ratio P with OD distribution matrixs used respectively1And P2, value rule such as formula (2)；
P in the present embodiment1=80%, P2=20%.
(G) traffic assignation is carried out using multipath-capacity restraint traffic assignment method, comprised the following steps that：
(G1) traffic assignation iterations is designated as K, and the OD amounts ratio of the secondary iteration loading of kth (k=1,2 ... ... K) is Qk,
AndMake k=1；
K=5 in the present embodiment, Q1=30%, Q2=25%, Q3=20%, Q4=15%, Q5=10%.
(G2) the OD distribution matrixs of each trip mode in step (D) are multiplied by Q respectivelyk, obtain being used for kth time iteration traffic
The OD distribution matrixs of distribution, start kth time iteration；
The OD distribution matrixs that the present embodiment is used for the 1st iteration are shown in Table 7-9.
The CarOD (pcu/h) of the 1st iteration of table 7
| Cell_ID | 1 | 2 | 3 | 4 | 5 | Traffic generating capacity adds up to |
| 1 | 900 | 300 | 1500 | 300 | 300 | 3300 |
| 2 | 300 | 900 | 1500 | 300 | 300 | 3300 |
| 3 | 1500 | 1500 | 300 | 1500 | 1500 | 6300 |
| 4 | 300 | 300 | 1500 | 900 | 300 | 3300 |
| 5 | 300 | 300 | 1500 | 300 | 900 | 3300 |
| Road traffic simulation amount adds up to | 3300 | 3300 | 6300 | 3300 | 3300 | 19500 |
The BusOD (pcu/h) of the 1st iteration of table 8
| Cell_ID | 1 | 2 | 3 | 4 | 5 | Traffic generating capacity adds up to |
| 1 | 21 | 9 | 21 | 9 | 6 | 66 |
| 2 | 9 | 21 | 21 | 6 | 9 | 66 |
| 3 | 21 | 21 | 30 | 21 | 21 | 114 |
| 4 | 9 | 6 | 21 | 21 | 9 | 66 |
| 5 | 6 | 9 | 21 | 9 | 21 | 66 |
| Road traffic simulation amount adds up to | 66 | 66 | 114 | 66 | 66 | 378 |
The BikeOD (people/h) of the 1st iteration of table 9
| Cell_ID | 1 | 2 | 3 | 4 | 5 | Traffic generating capacity adds up to |
| 1 | 600 | 150 | 150 | 150 | 18 | 1068 |
| 2 | 150 | 600 | 150 | 18 | 150 | 1068 |
| 3 | 150 | 150 | 150 | 150 | 150 | 750 |
| 4 | 150 | 18 | 150 | 600 | 150 | 1068 |
| 5 | 18 | 150 | 150 | 150 | 600 | 1068 |
| Road traffic simulation amount adds up to | 1068 | 1068 | 750 | 1068 | 1068 | 5022 |
(G3) for the trip mode restricted driving by tail number obtained in step (E), obtained in the step corresponded to (G2)
OD distribution matrixs be multiplied by P1The OD distribution matrixs (unrestricted tail number part) for first time traffic assignation are obtained, are corresponded to
The step of (G2) in obtained OD distribution matrixs be multiplied by P2Carry out the transfer amendment of OD distribution matrixs again afterwards, that is, obtain being used for second
The OD distribution matrixs (limited tail number part) of secondary traffic assignation；
For CarOD in the present embodiment, unrestricted tail number part OD distribution matrixs are shown in Table 10, are multiplied by P2The use obtained afterwards
11 are shown in Table respectively in the OD distribution matrixs of transfer amendment；
The unrestricted tail number part CarOD (pcu/h) of table 10
| Cell_ID | 1 | 2 | 3 | 4 | 5 | Traffic generating capacity adds up to |
| 1 | 720 | 240 | 1200 | 240 | 240 | 2640 |
| 2 | 240 | 720 | 1200 | 240 | 240 | 2640 |
| 3 | 1200 | 1200 | 240 | 1200 | 1200 | 5040 |
| 4 | 240 | 240 | 1200 | 720 | 240 | 2640 |
| 5 | 240 | 240 | 1200 | 240 | 720 | 2640 |
| Road traffic simulation amount adds up to | 2640 | 2640 | 5040 | 2640 | 2640 | 15600 |
Table 11 is used for the CarOD (pcu/h) for shifting amendment
| Cell_ID | 1 | 2 | 3 | 4 | 5 | Traffic generating capacity adds up to |
| 1 | 180 | 60 | 300 | 60 | 60 | 660 |
| 2 | 60 | 180 | 300 | 60 | 60 | 660 |
| 3 | 300 | 300 | 60 | 300 | 300 | 1260 |
| 4 | 60 | 60 | 300 | 180 | 60 | 660 |
| 5 | 60 | 60 | 300 | 60 | 180 | 660 |
| Road traffic simulation amount adds up to | 660 | 660 | 1260 | 660 | 660 | 3900 |
(G3a) the distribution volume of traffic that note traffic occurs between cell i and road traffic simulation cell j is Qij(i ∈ { Cell_ID },
J ∈ { Cell_ID }), remember QijTransfer amount be QTij, for all friendships corresponding to the cell number Cell_ID that is restricted driving by tail number
Logical generating capacity and traffic attraction, Schilling QTij=Qij, then make Qij=0；
In the present embodiment, for the traffic zone 3 in table 11, Schilling QT1*=Q1*、QT2*=Q2*、QT*1=Q*1、QT*2=
Q*2, then make Q1*=0, Q2*=0, Q*1=0, Q*1=0, herein in subscript " * " take 1 successively, 2 ..., 5；
(G3b) for remaining Q in OD distribution matrixsij, calculate the rate γ that detours under the conditions of tail number is restricted drivingij；
(G3b1) by by the traffic obtained in the numbering Cell_ID of tail number restricted driving cell and step (C) described in claim 1
Cell information is matched, and obtains the section numbering Link_ID in all sections that these traffic zones are included, and by claim 3
Corresponding journey time T is numbered in these sections in newest right of way table W under middle current procedureslValue be changed to be not less than 99999
Any number, obtain tail number restrict driving under the conditions of right of way table WL；In the 1st iterative process of the present embodiment, this step it is current
Newest right of way table W is table 2, and the right of way of the Link_ID that cell 3 is included in table 2 is changed into 99999, W is obtainedL。
(G3b2) the right of way table W under the conditions of being restricted driving respectively according to right of way table W and tail numberLOD is calculated to the shortest path between ij
Weigh DijAnd DLij；
In the present embodiment under current procedures, the shortest path power calculated according to table 2 between OD pairs is shown in Table 12, according to right of way table WL
Calculate the power of the shortest path between OD pairs and be shown in Table 13：
The shortest path of table 12 weighs Dij(s)
| Cell_ID | 1 | 2 | 3 | 4 | 5 |
| 1 | 236 | 418 | —— | 418 | 568 |
| 2 | 418 | 236 | —— | 568 | 418 |
| 3 | —— | —— | —— | —— | —— |
| 4 | 418 | 568 | —— | 236 | 418 |
| 5 | 568 | 418 | —— | 418 | 236 |
The shortest path of table 13 weighs DLij(s)
| Cell_ID | 1 | 2 | 3 | 4 | 5 |
| 1 | 236 | 418 | —— | 418 | 1412 |
| 2 | 418 | 236 | —— | 1412 | 418 |
| 3 | —— | —— | —— | —— | —— |
| 4 | 418 | 1412 | —— | 236 | 418 |
| 5 | 1412 | 418 | —— | 418 | 236 |
(G3b3) the rate γ that detours is calculated using formula (4)ij；
The rate of detouring obtained in the present embodiment under current procedures is shown in Table 14：
Table 14 detours rate γij
| Cell_ID | 1 | 2 | 3 | 4 | 5 |
| 1 | 1.0 | 1.0 | —— | 1.0 | 2.5 |
| 2 | 1.0 | 1.0 | —— | 2.5 | 1.0 |
| 3 | —— | —— | —— | —— | —— |
| 4 | 1.0 | 2.5 | —— | 1.0 | 1.0 |
| 5 | 2.5 | 1.0 | —— | 1.0 | 1.0 |
(G3c) OD amounts transfer function QTij=F (γij,Qij), calculate remaining QijTransfer amount QTij, then make Qij=Qij-
QTij, that is, complete the transfer amendment of OD distribution matrixs；
OD amount transfer functions Q described in step (G3c)Tij=F (γij,Qij) expression formula be formula (5)：
Wherein β1、β0For coefficient to be calibrated, γ0To specify all volume of traffic under confidence level not change trip mode
Maximum transfer rate, γ1The minimum rate of transform of trip mode, γ are changed for all volume of traffic under specified confidence level0And γ1Can
Directly obtained by survey, β can be demarcated using formula (6)1And β0：
In this step of the present embodiment, γ is got under 0.95 confidence level0=1.2, γ1=3.0, according to formula (6)
(5) calculated, the results are shown in Table 15,16：
Table 15 is limited tail number part CarOD (pcu/h)
| Cell_ID | 1 | 2 | 3 | 4 | 5 |
| 1 | 180 | 60 | 0 | 60 | 17 |
| 2 | 60 | 180 | 0 | 17 | 60 |
| 3 | 0 | 0 | 0 | 0 | 0 |
| 4 | 60 | 17 | 0 | 180 | 60 |
| 5 | 17 | 60 | 0 | 60 | 180 |
Table 16CarOD diverted traffic (pcu/h)
| Cell_ID | 1 | 2 | 3 | 4 | 5 |
| 1 | 0 | 0 | 60 | 0 | 43 |
| 2 | 0 | 0 | 30 | 43 | 0 |
| 3 | 300 | 300 | 60 | 300 | 300 |
| 4 | 0 | 43 | 180 | 0 | 0 |
| 5 | 43 | 0 | 60 | 0 | 0 |
(G4) for the trip mode do not restricted driving in step (D) by tail number, shifted using OD distribution matrixs in step (G3)
The diverted traffic produced during amendment updates the trip mode OD distribution matrixs do not restricted driving by tail number, with reference to right of way table W and renewal
The trip mode OD distribution matrixs do not restricted driving afterwards by tail number carry out kth time iteration traffic assignation, and the volume of traffic after distribution is tired out
Meter is loaded into corresponding road section automobile traffic amount VlOn (unit is pcu/h)；
By the diverted traffic between in table 16 each OD pairs in the present embodiment, according to the people between table 8 and same OD pairs of table 9
Secondary volume of traffic ratio is split and each superposition is entered, and is reconverted into respective original unit's volume of traffic, such as QT15=
Person-time volume of traffic ratio in 43pcu/h, table 8 and table 9 between same OD pairs is (6 × 40=240)：18=40：3, then by 40
Person-time diverted traffic be converted to the 1pcu/h regular public traffic trip volume of traffic and be superimposed to the Q of table 815In, by 3 person-times of transfer
The volume of traffic is that 3 person-times of the non-motor vehicle trip volume of traffic is superimposed to the Q of table 915In.BusOD and BikeOD after renewal are shown in respectively
Table 17 and table 18.
BusOD after the renewal of table 17
| Cell_ID | 1 | 2 | 3 | 4 | 5 | Traffic generating capacity adds up to |
| 1 | 21 | 9 | 21 | 9 | 7 | 67 |
| 2 | 9 | 21 | 21 | 7 | 9 | 67 |
| 3 | 21 | 21 | 30 | 21 | 21 | 114 |
| 4 | 9 | 7 | 21 | 21 | 9 | 67 |
| 5 | 7 | 9 | 21 | 9 | 21 | 67 |
| Road traffic simulation amount adds up to | 67 | 67 | 114 | 67 | 67 | 382 |
BikeOD after the renewal of table 18
| Cell_ID | 1 | 2 | 3 | 4 | 5 | Traffic generating capacity adds up to |
| 1 | 600 | 150 | 150 | 150 | 21 | 1071 |
| 2 | 150 | 600 | 150 | 21 | 150 | 1071 |
| 3 | 150 | 150 | 150 | 150 | 150 | 750 |
| 4 | 150 | 21 | 150 | 600 | 150 | 1071 |
| 5 | 21 | 150 | 150 | 150 | 600 | 1071 |
| Road traffic simulation amount adds up to | 1071 | 1071 | 750 | 1071 | 1071 | 5034 |
In current procedures, the 1st iteration for carrying out automobile traffic amount using table 2 and table 17 is distributed, and by after distribution
The volume of traffic is accumulative to be loaded into corresponding road section automobile traffic amount VlOn.
(G5) according to each section automobile traffic amount V obtained in (G4)l, right of way table W Zhong Ge roads are updated using formula (3)
The journey time T of sectionl；
Right of way table W is shown in Table 19 after the renewal that this step of this example is obtained.
The right of way table W of table 19
(G6) after using the unrestricted tail number part OD distribution matrixs obtained in step (G3) and step (G5) updated
Right of way table W carries out first time traffic assignation, and is loaded into corresponding road section automobile traffic amount V by the volume of traffic after distribution is accumulativel
On；
In this step of the present embodiment, first time traffic assignation is carried out using table 10 and table 19, and by the volume of traffic after distribution
It is accumulative to be loaded into corresponding road section automobile traffic amount VlOn.
(G7) according to each section automobile traffic amount V obtained in step (G6)l, updated using formula (3) in right of way table W
The journey time T in each sectionl；
(G8) by by the traffic zone information matches obtained in the numbering Cell_ID of tail number restricted driving cell and step (C), obtain
The section numbering Link_ID in all sections included to these traffic zones, and by the right of way table W obtained in step (G7)
Number corresponding journey time T in these sectionslValue be changed to any number not less than 99999, obtain tail number restricted driving condition
Under right of way table WL；
Right of way table W under the conditions of the tail number restricted driving that current procedures are obtainedLIt is shown in Table 20.
The right of way table W of table 20L
| Section numbering Link_ID | Right of way (s) |
| 1 | 30 |
| 2 | 19 |
| 3 | 15 |
| 4 | 26 |
| …… | …… |
| 8 | 99999 |
| …… | …… |
| 11 | 99999 |
| …… | …… |
| 76 | 37 |
(G9) using the tail number obtained in the limited tail number part OD distribution matrixs and step (G8) obtained in step (G3)
Right of way table W under the conditions of restricted drivingLCarry out second of traffic assignation, and be loaded into corresponding road section machine by the volume of traffic after distribution is accumulative
Motor-car volume of traffic VlOn；
Traffic assignation is carried out using table 15 and table 20 in the present embodiment current procedures, and by the volume of traffic after distribution it is accumulative plus
It is downloaded to corresponding road section automobile traffic amount VlOn.
(G10) according to each section automobile traffic amount V obtained in step (G9)l, right of way table W is updated using formula (3)
In each section journey time Tl；
Right of way table W is shown in Table 21 after the renewal obtained in the present embodiment current procedures.
The right of way table W of table 21
| Section numbering Link_ID | Right of way (s) |
| 1 | 35 |
| 2 | 22 |
| 3 | 17 |
| 4 | 30 |
| …… | …… |
| 8 | 42 |
| …… | …… |
| 11 | 58 |
| …… | …… |
| 76 | 41 |
(G11) make k=k+1, if k≤K, return to step (G2), otherwise iteration traffic assignation terminate；
K is made to take 2,3,4,5 successively, return to step (G2) continuation the 2nd, 3,4,5 iteration traffic assignations complete the 5th and changed
After traffic assignation, iteration terminates.
(H) the automobile traffic amount V that each section is distributed is exportedl。
By 5 iteration traffic assignations in the present embodiment, the automobile traffic amount V that final each section is distributed is obtainedl。
In this example, it can be seen that by eight steps of the inventive method, can quickly and easily complete region motor-driven
Point mode OD distribution matrixs amendment and traffic assignation process under the conditions of the tailstock number restricted driving, the logicality of the modification method are strong, easy
Realized in sequencing.
Embodiment 2：
The present embodiment and the difference of embodiment 1 be, the method described in the present embodiment, it is possible to achieve only in city
The motor vehicle tail number restricted driving scheme that part major trunk roads are implemented, be specifically：
The step (E) includes work to obtain the urban parts major trunk roads motor vehicle tail number restricted driving scheme for testing
Restricted motor vehicle tail number number N in the daily restricted driving time in dayL(unit be it is individual, value be 1,2 ..., 9), tail number limit
Trip mode (trip mode restricted driving by tail number), tail number restricted driving scheme in step (D) corresponding to row scheme are covered
The numbering Link_ID in section (by tail number restricted driving section)；
Urban parts major trunk roads motor vehicle tail number restricted driving scheme in the present embodiment for section 4,8,42,46 on weekdays
07:00~09:00、17:00~19:00 for 2 tail numbers of private minibus restricted driving, i.e. NL=2, the trip side restricted driving by tail number
Formula is private minibus, is 4,8,42,46 by the numbering Link_ID values in tail number restricted driving section.
The step (G3) for " for the trip mode restricted driving by tail number obtained in step (E), the step corresponded to
(G2) the OD distribution matrixs obtained in are multiplied by P1Obtain OD distribution matrixs (the unrestricted tail number portion for first time traffic assignation
Point), the OD distribution matrixs obtained in the step (G2) corresponded to are multiplied by P2Obtain the OD distributions for second of traffic assignation
Matrix (limited tail number part) "；
For CarOD in the present embodiment, unrestricted tail number part OD distribution matrixs are limited tail number part OD distributions with table 10
Matrix is with table 11.
The step (G4) is " for the trip mode do not restricted driving in step (D) by tail number, with reference to right of way table W and step
(G2) the trip mode OD distribution matrixs do not restricted driving by tail number obtained in carry out kth time iteration traffic assignation, and by after distribution
The volume of traffic accumulative be loaded into corresponding road section automobile traffic amount VlOn "；
In the present embodiment current procedures, the 1st iteration for carrying out automobile traffic amount using table 2 and table 8 is distributed, and will be divided
The volume of traffic after matching somebody with somebody adds up to be loaded into corresponding road section automobile traffic amount VlOn.
The step (G8) for " by the right of way table W obtained in step (G7) by Link_ID pairs of tail number restricted driving section numbering
The journey time T answeredlValue be changed to any number not less than 99999, obtain the right of way table W under the conditions of tail number is restricted drivingL”。
Right of way table W under the conditions of the tail number restricted driving that current procedures are obtainedLIt is shown in Table 21.
The right of way table W of table 21L
| Section numbering Link_ID | Right of way (s) |
| 1 | 31 |
| 2 | 20 |
| 3 | 15 |
| 4 | 99999 |
| …… | …… |
| 8 | 99999 |
| …… | …… |
| 11 | 48 |
| …… | …… |
| 76 | 36 |
Remaining step is same as Example 1.
In this example, it can be seen that equally can quickly and easily be completed only in urban parts master by the inventive method
The traffic assignation process under conditions of motor vehicle tail number restricted driving scheme is implemented in arterial highway, now without entering to a point mode OD distribution matrixs
Row amendment, the logicality of the modification method is strong, be easy to sequencing realization.
It should be pointed out that above-mentioned embodiment is only intended to clearly illustrate example, and not to embodiment
Limit, there is no necessity and possibility to exhaust all the enbodiments.Each part being not known in the present embodiment
Realized with prior art.For those skilled in the art, the premise of the principle of the invention is not being departed from
Under, some improvements and modifications can also be made, these improvements and modifications also should be regarded as protection scope of the present invention.
Claims (9)
1. the traffic assignation modification method under the conditions of a kind of tail number restricted driving, it is characterised in that this method comprises the following steps：
(A) the basic road network information in newest city, including intersection numbering Node_ID (dimensionless), section numbering Link_ are obtained
ID (dimensionless), road section length Ll(unit is m) and section practical capacity Cl(unit is pcu/h)；
(B) the initial right of way table W for traffic assignation is obtained, its method is each section numbering Link_ID in first obtaining step (A)
Corresponding desin speed SD(unit is km/h) and Maximum speed limit SL(unit is km/h), is then calculated using formula (1) and obtained
Free flow journey time T corresponding to each section numbering Link_IDl0(unit is s), corresponding to the numbering Link_ID of these sections
Journey time be initial right of way table W：
Tl0=3600 × { 0.001 × Ll/min{SD,SL}} (1)；
(C) the traffic zone information distributed for urban transportation, the numbering Cell_ID (nothings comprising each traffic zone are obtained
Dimension), each traffic zone the node serial number Node_ID and section numbering Link_ID that are included；
(D) the newest each trip mode OD distribution matrixs in city are obtained；
(E) the region motor vehicle tail number restricted driving scheme for test is obtained, comprising being restricted in working day in the daily restricted driving time
Motor vehicle tail number number NL(unit is individual, value is 1,2 ..., in step (D) 9), corresponding to tail number restricted driving scheme
Traffic zone is (by tail in the completely covered step (C) in trip mode (trip mode restricted driving by tail number), tail number restricted driving region
Number restricted driving cell) numbering Cell_ID；
Or：
The urban parts major trunk roads motor vehicle tail number restricted driving scheme for test is obtained, comprising in the daily restricted driving time in working day
Restricted motor vehicle tail number number NL(unit is individual, value is 1,2 ..., step 9), corresponding to tail number restricted driving scheme
(D) the numbering Link_ID " by tail number restricted driving section that the trip mode restricted driving by tail number, tail number restricted driving scheme in are covered；
(F) the OD distribution matrixs by tail number restricted driving trip mode determined in step (E) are split as into two parts to be respectively used to twice
In traffic assignation, motor vehicle tail number number N is restricted according to what is determined in step (E)LIt is determined that dividing for the first time with second of traffic
Ratio P with OD distribution matrixs used respectively1And P2, value rule such as formula (2)：
(G) multipath-capacity restraint traffic assignment method under the conditions of being restricted driving using tail number is iterated traffic assignation；
(H) the automobile traffic amount V that each section is distributed is exportedl。
2. the traffic assignation modification method under the conditions of a kind of tail number restricted driving according to claim 1, it is characterised in that step
(D) each trip mode OD distribution matrixs in city described in, including：Private minibus trip mode OD distribution matrixs CarOD is (single
Position be pcu/h), regular public traffic trip mode OD distribution matrixs BusOD (unit is pcu/h), track traffic trip mode OD divide
Cloth matrix RailOD (unit behaviour/h), motorcycle trip mode OD distribution matrixs MotorOD (unit is pcu/h), taxi
(unit is by trip mode OD distribution matrixs TaxiOD (unit is pcu/h), non-motor vehicle trip mode OD distribution matrixs BikeOD
People/h) and walking trip mode OD distribution matrixs WalkOD (unit behaviour/h) in it is several.
3. the traffic assignation modification method under the conditions of a kind of tail number restricted driving according to claim 1, it is characterised in that described
The trip mode OD distribution matrixs restricted driving by tail number are split as unrestricted tail number part and limited tail number part two in step (F)
Part, and traffic assignation twice is carried out, the unrestricted tail number part OD distribution matrixs of the first sub-distribution, the second sub-distribution is limited tail number OD
Distribution matrix.
4. the traffic assignation modification method under the conditions of a kind of tail number restricted driving according to claim 1, it is characterised in that described
Multipath-capacity restraint traffic assignment method under the conditions of the tail number restricted driving used in step (G), its iteration traffic assignation process
Comprise the steps of：
(G1) traffic assignation iterations is designated as K, and the OD amounts ratio of the secondary iteration loading of kth (k=1,2 ... ... K) is Qk, andMake k=1；
(G2) the OD distribution matrixs of each trip mode in step (D) are multiplied by Q respectivelyk, obtain being used for kth time iteration traffic assignation
OD distribution matrixs, start kth time iteration；
(G3) for the trip mode restricted driving by tail number obtained in step (E), the OD obtained in the step corresponded to (G2)
Distribution matrix is multiplied by P1The OD distribution matrixs (unrestricted tail number part) for first time traffic assignation are obtained, by its corresponding step
Suddenly the OD distribution matrixs obtained in (G2) are multiplied by P2Carry out the transfer amendment of OD distribution matrixs again afterwards, that is, obtain being used for second of friendship
The OD distribution matrixs (limited tail number part) that the reduction of fractions to a common denominator is matched somebody with somebody；
Or：
For the trip mode restricted driving by tail number obtained in step (E), the OD distributions obtained in the step corresponded to (G2)
Matrix Multiplication is with P1Obtain the OD distribution matrixs (unrestricted tail number part) for first time traffic assignation, the step corresponded to
(G2) the OD distribution matrixs obtained in are multiplied by P2Obtain the OD distribution matrixs (limited tail number part) for second of traffic assignation；
(G4) for the trip mode do not restricted driving in step (D) by tail number, corrected using OD distribution matrixs transfer in step (G3)
When the diverted traffic that produces update the trip mode OD distribution matrixs do not restricted driving by tail number, with reference to right of way table W and after updating not
The trip mode OD distribution matrixs restricted driving by tail number carry out kth time iteration traffic assignation, and by the volume of traffic after distribution it is accumulative plus
It is downloaded to corresponding road section automobile traffic amount VlOn (unit is pcu/h)；
Or：
For the trip mode do not restricted driving in step (D) by tail number, with reference to being obtained in right of way table W and step (G2) not by tail number
The trip mode OD distribution matrixs of restricted driving carry out kth time iteration traffic assignation, and are loaded into the volume of traffic after distribution is accumulative pair
Answer section automobile traffic amount VlOn；
(G5) according to each section automobile traffic amount V obtained in step (G4)l, using Bureau of Public Road's Impedance Function mould
Type updates the journey time T in each section in right of way table Wl(unit is s), i.e. formula (3)：
Tl=Tl0[1+α(Vl/Cl)β] (3)
Wherein α=0.15, β=4；
(G6) right of way using the unrestricted tail number part OD distribution matrixs and step (G5) that are obtained in step (G3) after updated
Table W carries out first time traffic assignation, and is loaded into corresponding road section automobile traffic amount V by the volume of traffic after distribution is accumulativelOn；
(G7) according to each section automobile traffic amount V obtained in step (G6)l, right of way table W Zhong Ge roads are updated using formula (3)
The journey time T of sectionl；
(G8) by by the traffic zone information matches obtained in the numbering Cell_ID of tail number restricted driving cell and step (C), this is obtained
The section numbering Link_ID in all sections that a little traffic zones are included, and by these in the right of way table W obtained in step (G7)
Number corresponding journey time T in sectionlValue be changed to any number not less than 99999, obtain under the conditions of tail number restricts driving
Right of way table WL；
Or：
By in the right of way table W obtained in step (G7) by the corresponding journey time T of tail number restricted driving section numbering Link_IDlValue more
Any number not less than 99999 is changed to, the right of way table W under the conditions of tail number is restricted driving is obtainedL；
(G9) restricted driving using the tail number obtained in the limited tail number part OD distribution matrixs and step (G8) obtained in step (G3)
Under the conditions of right of way table WLCarry out second of traffic assignation, and be loaded into corresponding road section motor vehicle by the volume of traffic after distribution is accumulative
Volume of traffic VlOn；
(G10) according to each section automobile traffic amount V obtained in step (G9)l, right of way table W Zhong Ge roads are updated using formula (3)
The journey time T of sectionl；
(G11) make k=k+1, if k≤K, return to step (G2), otherwise iteration traffic assignation terminate.
5. the traffic assignation modification method under the conditions of a kind of tail number restricted driving according to claim 4, it is characterised in that described
The process that the transfer amendment of OD distribution matrixs is carried out in step (G3) is comprised the steps of：
(G3a) the distribution volume of traffic that note traffic occurs between cell i and road traffic simulation cell j is Qij(i ∈ { Cell_ID }, j ∈
{ Cell_ID }), remember QijTransfer amount be QTij, for all traffic corresponding to the cell number Cell_ID that is restricted driving by tail number
Generating capacity and traffic attraction, Schilling QTij=Qij, then make Qij=0；
(G3b) for remaining Q in OD distribution matrixsij, calculate the rate γ that detours under the conditions of tail number is restricted drivingij；
(G3c) OD amount transfer functions Q is combinedTij=F (γij,Qij), calculate remaining QijTransfer amount QTij, then make Qij=Qij-
QTij, that is, complete the transfer amendment of OD distribution matrixs.
6. the traffic assignation modification method under the conditions of a kind of tail number restricted driving according to claim 5, it is characterised in that described
The rate γ that detours under the conditions of tail number is restricted driving is calculated in step (G3b)ijThe step of it is as follows：
(G3b1) by by the traffic zone obtained in the numbering Cell_ID of tail number restricted driving cell and step (C) described in claim 1
Information matches, obtain the section numbering Link_ID in all sections that these traffic zones are included, and by claim 4 when
Corresponding journey time T is numbered in these sections in newest right of way table W under preceding steplValue be changed to not less than 99999 appoint
Meaning numerical value, obtains the right of way table W under the conditions of tail number is restricted drivingL；
(G3b2) the right of way table W under the conditions of being restricted driving respectively according to right of way table W and tail numberLCalculate OD and D is weighed to the shortest path between ijij
And DLij；
(G3b3) the rate γ that detours is calculated using formula (4)ij：
γij=DLij/Dij (4)。
7. the traffic assignation modification method under the conditions of a kind of tail number restricted driving according to claim 5, it is characterised in that described
OD amount transfer functions Q described in step (G3c)Tij=F (γij,Qij) expression formula be formula (5)：
Wherein β1、β0For coefficient to be calibrated, γ0The maximum of trip mode is not changed for all volume of traffic under specified confidence level
The rate of transform, γ1The minimum rate of transform of trip mode, γ are changed for all volume of traffic under specified confidence level0And γ1It can pass through
Survey is directly obtained, and β can be demarcated using formula (6)1And β0：
8. the traffic assignation modification method under the conditions of a kind of tail number restricted driving according to claim 5, it is characterised in that described
The method for the trip mode OD distribution matrixs do not restricted driving by tail number is updated in step (G4) is：By obtained all diverted traffics
QTijAccording to not split simultaneously to an ij person-time volume of traffic ratio by same OD in the trip mode OD distribution matrixs that tail number is restricted driving
Each superposition is entered, and is reconverted into OD points of the trip mode do not restricted driving after respective original unit's volume of traffic is updated by tail number
Cloth matrix.
9. the traffic assignation modification method under the conditions of a kind of tail number restricted driving according to claim 8, it is characterised in that：It is private
Minibus, regular public traffic, person-time volume of traffic of motorcycle and taxi these four trip modes and original unit's volume of traffic mutually turn
The conversion coefficient changed respectively θcar=1.3, θbus=40.0, θmotor=1.0, θtaxi=1.3, original unit's volume of traffic is multiplied by conversion
Coefficient is to obtain a person-time volume of traffic, and person-time volume of traffic divided by conversion coefficient are to be changed into original unit's volume of traffic；Track traffic, non-maneuver
A person-time volume of traffic for these three trip modes of car and walking is original unit's volume of traffic, without conversion.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710306621.6A CN107038862B (en) | 2017-05-04 | 2017-05-04 | Traffic distribution correction method under tail number restriction condition |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710306621.6A CN107038862B (en) | 2017-05-04 | 2017-05-04 | Traffic distribution correction method under tail number restriction condition |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107038862A true CN107038862A (en) | 2017-08-11 |
| CN107038862B CN107038862B (en) | 2020-02-14 |
Family
ID=59536847
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710306621.6A Active CN107038862B (en) | 2017-05-04 | 2017-05-04 | Traffic distribution correction method under tail number restriction condition |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107038862B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108269399A (en) * | 2018-01-24 | 2018-07-10 | 哈尔滨工业大学 | A kind of high ferro passenger forecast method based on the anti-push technologies of network of highways passenger flow OD |
| CN108399465A (en) * | 2017-10-12 | 2018-08-14 | 东南大学 | The OD distribution methods of area traffic control strategy are implemented in a kind of reply |
| CN108399466A (en) * | 2017-10-12 | 2018-08-14 | 东南大学 | A kind of OD distribution methods being accustomed in view of resident trip |
| CN108765941A (en) * | 2018-05-29 | 2018-11-06 | 重庆大学 | A kind of signalized intersections vehicle arriving rate method of estimation |
| CN109325617A (en) * | 2018-09-04 | 2019-02-12 | 青岛海信网络科技股份有限公司 | A kind of urban traffic status prediction technique and device |
| CN111275965A (en) * | 2020-01-20 | 2020-06-12 | 交通运输部科学研究院 | Real-time traffic simulation analysis system and method based on internet big data |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103440764A (en) * | 2013-08-19 | 2013-12-11 | 同济大学 | Urban road network vehicle travel path reconstruction method based on vehicle automatic identification data |
| CN103593535A (en) * | 2013-11-22 | 2014-02-19 | 南京洛普股份有限公司 | Urban traffic complex self-adaptive network parallel simulation system and method based on multi-scale integration |
| CN104183119A (en) * | 2014-08-19 | 2014-12-03 | 中山大学 | Real-time traffic flow distribution prediction system based on road section OD backstepping |
-
2017
- 2017-05-04 CN CN201710306621.6A patent/CN107038862B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103440764A (en) * | 2013-08-19 | 2013-12-11 | 同济大学 | Urban road network vehicle travel path reconstruction method based on vehicle automatic identification data |
| CN103593535A (en) * | 2013-11-22 | 2014-02-19 | 南京洛普股份有限公司 | Urban traffic complex self-adaptive network parallel simulation system and method based on multi-scale integration |
| CN104183119A (en) * | 2014-08-19 | 2014-12-03 | 中山大学 | Real-time traffic flow distribution prediction system based on road section OD backstepping |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108399465A (en) * | 2017-10-12 | 2018-08-14 | 东南大学 | The OD distribution methods of area traffic control strategy are implemented in a kind of reply |
| CN108399466A (en) * | 2017-10-12 | 2018-08-14 | 东南大学 | A kind of OD distribution methods being accustomed in view of resident trip |
| CN108399465B (en) * | 2017-10-12 | 2021-09-07 | 东南大学 | OD distribution method for implementing regional traffic management strategy |
| CN108399466B (en) * | 2017-10-12 | 2021-09-07 | 东南大学 | OD distribution method considering resident travel habits |
| CN108269399A (en) * | 2018-01-24 | 2018-07-10 | 哈尔滨工业大学 | A kind of high ferro passenger forecast method based on the anti-push technologies of network of highways passenger flow OD |
| CN108269399B (en) * | 2018-01-24 | 2020-12-25 | 哈尔滨工业大学 | High-speed rail passenger flow demand prediction method based on road network passenger flow OD reverse thrust technology |
| CN108765941A (en) * | 2018-05-29 | 2018-11-06 | 重庆大学 | A kind of signalized intersections vehicle arriving rate method of estimation |
| CN108765941B (en) * | 2018-05-29 | 2022-03-11 | 重庆大学 | Signalized intersection vehicle arrival rate estimation method |
| CN109325617A (en) * | 2018-09-04 | 2019-02-12 | 青岛海信网络科技股份有限公司 | A kind of urban traffic status prediction technique and device |
| CN111275965A (en) * | 2020-01-20 | 2020-06-12 | 交通运输部科学研究院 | Real-time traffic simulation analysis system and method based on internet big data |
| CN111275965B (en) * | 2020-01-20 | 2021-02-05 | 交通运输部科学研究院 | Real-time traffic simulation analysis system and method based on internet big data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107038862B (en) | 2020-02-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107038862A (en) | A kind of traffic assignation modification method under the conditions of tail number restricted driving | |
| Wu et al. | Energy-optimal speed control for electric vehicles on signalized arterials | |
| Davidson et al. | Autonomous vehicles: what could this mean for the future of transport | |
| Werner et al. | Passenger car equivalencies of trucks, buses, and recreational vehicles for two-lane rural highways | |
| CN105216782B (en) | Plug-in hybrid-power automobile energy management method based on energy predicting | |
| CN104866654B (en) | A kind of construction method of integrated urban dynamic traffic emulation platform | |
| CN107330559A (en) | A kind of hybrid customization public bus network planing method of many terminus multi-vehicle-types | |
| CN107392400A (en) | Meter and real-time traffic and the EV of temperature charging load spatial and temporal distributions Forecasting Methodology | |
| CN105205557A (en) | Design method for conventional urban public transit network | |
| CN105946857A (en) | Parallel plug-in hybrid electric vehicle (PHEV) energy management method based on intelligent transportation system | |
| CN105857312B (en) | A kind of highway heavy truck speed travels optimization method | |
| CN112309119B (en) | Urban traffic system capacity analysis optimization method | |
| CN107921886A (en) | Method for calculating the set point for being used for the fuel and power consumption for managing hybrid moto vehicle | |
| DE102014224758A1 (en) | METHOD AND DEVICE FOR PROGNOSTICATING ELECTRIC VEHICLE ENERGY RECEIPT | |
| CN102855386A (en) | Footway and non-motorized vehicle lane evaluating method and device | |
| Arasan et al. | Characteristics of trips by foot and bicycle modes in Indian city | |
| CN105912748A (en) | Method of evaluating contribution degree of each local mode of vehicle body structure | |
| CN108171975B (en) | Urban automobile running speed prediction method based on road section and intersection distribution | |
| CN113936475A (en) | Traffic control method and system for optimizing multi-road resources | |
| Schücking et al. | Influencing factors on specific energy consumption of EV in extensive operations | |
| DE102017213984A1 (en) | Method for operating a navigation device for a motor vehicle | |
| Unni et al. | Energy Consumption Analysis for the Prediction of Battery Residual Energy in Electric Vehicles | |
| CN103426294B (en) | A kind of division methods based on the other multi-layer Fuzzy of emergency traffic flow priority | |
| CN117690301B (en) | Expressway diversion induction method considering induction compliance rate | |
| CN107124294A (en) | A kind of In-vehicle networking evaluation method based on data fusion |
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 |