CN106920395A - A kind of traffic impedance computation method based on parameter calibration - Google Patents

A kind of traffic impedance computation method based on parameter calibration Download PDF

Info

Publication number
CN106920395A
CN106920395A CN201710263992.0A CN201710263992A CN106920395A CN 106920395 A CN106920395 A CN 106920395A CN 201710263992 A CN201710263992 A CN 201710263992A CN 106920395 A CN106920395 A CN 106920395A
Authority
CN
China
Prior art keywords
section
parameter
flow
parameter calibration
speed
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
CN201710263992.0A
Other languages
Chinese (zh)
Other versions
CN106920395B (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.)
HANGZHOU COMPREHENSIVE TRANSPORTATION RESEARCH CENTER
Original Assignee
HANGZHOU COMPREHENSIVE TRANSPORTATION RESEARCH CENTER
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 HANGZHOU COMPREHENSIVE TRANSPORTATION RESEARCH CENTER filed Critical HANGZHOU COMPREHENSIVE TRANSPORTATION RESEARCH CENTER
Priority to CN201710263992.0A priority Critical patent/CN106920395B/en
Publication of CN106920395A publication Critical patent/CN106920395A/en
Application granted granted Critical
Publication of CN106920395B publication Critical patent/CN106920395B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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/0108—Measuring and analyzing of parameters relative to traffic conditions based on the source of data
    • G08G1/012—Measuring and analyzing of parameters relative to traffic conditions based on the source of data from other sources than vehicle or roadside beacons, e.g. mobile networks
    • 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
    • G08G1/0129—Traffic data processing for creating historical data or processing based on historical data

Landscapes

  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Traffic Control Systems (AREA)

Abstract

The invention discloses a kind of traffic impedance computation method based on parameter calibration.The present invention is to be processed link proportion function, is converted into simple regression problem, then obtains parameter to be calibrated using least square method.Objectivity, scientific shortcoming instant invention overcomes Traditional Man empirical method, and data are acquired higher, possess stronger actual operation.

Description

A kind of traffic impedance computation method based on parameter calibration
Technical field
The present invention relates to a kind of traffic impedance computation method based on parameter calibration, belong to traffic programme technical field.
Background technology
Impedance Function can assess the congested in traffic influence to journey time, and this is the basis of traffic assignation and equilibrium analysis And premise.Impedance Function includes intersection impedance function and link proportion function.It is domestic in terms of the research of link proportion function Many scholars and research institution are proposed different Mathematical Modelings, and wherein the value of model parameter is generally according to artificial experience pair Numerical value is judged, and is adjusted according to final calculation result, and objectivity, science have been short of.
The content of the invention
The defect of technical elements is being solved in view of above-mentioned existing traffic Impedance Function, the present invention is based on magnanimity traffic flow letter Breath, proposes a kind of traffic impedance computation method based on parameter calibration.
The acquisition of basic data and treatment, including obtain the section flow and speed for parameter calibration；Described section The whole section section that flow is detected by radio telephone mobile system each bicycle road flow conversion in 5 minutes is formed, described The velocity amplitude that is detected using radio telephone mobile system of speed.
It is determined that intending the continuous stream or interruption flowpath segment demarcated, the shifting of the section working day radio telephone of continuous month is chosen The data that dynamic system detectio is arrived, to parameter calibration.
Build section Impedance Function model：
In formula：T represents the link travel time between two intersections；tfWhen the expression volume of traffic is 0, between two intersections Link travel time；V represents section automobile traffic amount；C represents section actual capacity；α, β represent to be calibrated respectively Impedance Function parameter.
Parameter calibration is carried out to building section Impedance Function model, specifically：
S1, logarithmetics treatment is carried out to section Impedance Function：
For certain section LiHave：
Section Impedance Function after taking the logarithm is expressed as：
In formula：V represents the sections of road speed between two intersections, vfRepresent section free stream velocity.
Ifβ=k,Ln α=b, then have y=kx+b；It is translated into simple regression point Analysis problem.
S2, calibrating parameters C, vf
According to speed, the corresponding relation of flow detection value, speed-flow scatter diagram is drawn, scatter diagram is fitted, marked Make parameter C and vf。
S3, obtain parameter alpha to be calibrated, β using least square method.
Beneficial effects of the present invention：The present invention is handed on the basis of BRP typical case's roadlock models according to magnanimity such as flow, speeds Through-flow information, proposition is a kind of based on the traffic impedance computation method that parameter calibration is carried out using least square method, overcomes tradition The objectivity of artificial experience method, scientific shortcoming, and data are acquired higher, possess stronger actual operation.
Brief description of the drawings
Fig. 1 is the literary west section matched curve eastwards of (Gu Cuilu-Xueyuan Road) 4 tracks all the way schematic diagram.
Fig. 2 is the literary section of (Xueyuan Road-Gu Cuilu) 4 track east orientation West Road all the way matched curve schematic diagram.
Fig. 3 is Wang Jianglu (Jiang Chenglu-Qiu Taolu) 3 tracks west section matched curve eastwards schematic diagram.
Fig. 4 is the track east orientation West Roads of Wang Jianglu (Qiu Taolu-Jiang Chenglu) 3 section matched curve schematic diagram.
Specific embodiment
Basic thought of the invention is to be processed link proportion function, is converted into simple regression problem, is then utilized Least square method obtains parameter to be calibrated.
Basic step of the invention is as follows：
C1, the acquisition of basic data and treatment.
C2, section Impedance Function model construction and parameter calibration.
Step c1, the acquisition of basic data and treatment.Specifically include:
Obtain section flow, the speed data for parameter calibration：
Road traffic delay can be divided into Uninterrupted and interrupted flows.Typically continuous stream is divided three classes：It is highway, quick Road/viaduct, ring road；Typically stopped by between and be divided into four classes according to category of roads：Trunk roads, secondary distributor road, main branch road, secondary branch Road.Different according to form of fracture is six classes by interruption flow point：The non-separation of central strip machine, the non-line of central strip machine, center point Every Computer method, center do not separate the non-separation of machine, center do not separate the non-line of machine, center do not separate Computer method.Continuous stream and Between stop data source it is as follows：
(1) flow value of radio telephone mobile system RTMS (radio telephone mobile system) detections.Will Each bicycle road flow conversions of whole section section 5min that RTMS is detected are whole section vehicle flowrate, and expand calculation for 1h flows. Unit：pcu/h.
(2) velocity amplitude of RTMS detections.The velocity amplitude of detection is one-to-one relationship with flow value.Unit：km/h.
(3) data volume.It is determined that intending the continuous stream or interruption flowpath segment demarcated, the section working day of continuous one month is chosen RTMS detection datas, to parameter calibration.
Step c2, section Impedance Function model construction and parameter calibration.Specifically include：
C21, selection section Impedance Function model
Current transcad does not receive User Defined Impedance Function form, uses the Bureau of Public Roads after modification parameter BPR functions：
In formula：T --- the link travel time (min) between two intersections；
tf--- when the volume of traffic is 0, the link travel time (min) between two intersections；
V --- section automobile traffic amount (/h)；
C --- section actual capacity (/h)；
α, β --- Impedance Function parameter.
C22, parameter calibration
C221, logarithmetics treatment is carried out to BRP function formulas
Before carrying out parameter calibration to link proportion function, logarithmetics treatment is carried out to BRP function formulas, obtained：
For certain section LiHave：
Therefore, the BRP functions after taking the logarithm can also be expressed as：
In formula：T, tf, V, C are constants, and implication is ibid；
V --- the sections of road speed between two intersections, km/h；
vf--- section free stream velocity, km/h；
Ifβ=k,Ln α=b, then have y=kx+b.Simple regression point can be converted into Analysis problem.Section travel speed v and link counting V can be detected by RTMS and obtained in formula, therefore want calibrating parameters α, β, Firstly the need of calibrating C and vfValue.
C222, calibrating parameters C, vf。
According to RTMS speed, the corresponding relation of flow detection value, speed-flow scatter diagram is drawn, scatter diagram is intended Close, parameter can be calibrated：Section bicycle road traffic capacity C, section free stream velocity vf。
C223, obtain parameter alpha to be calibrated, β using least square method
The general principle of least square method is that the residual sum of squares (RSS) for making model reaches minimum value, can be expressed as with formula：
Q (b, k) is residual sum of squares (RSS) in formulaB the and k values of minimum are reached, is designated as
Least square method is really to find appropriate k and b so that function Q values reach minimum in formula, are asked according in calculus The principle of extreme value is rightWithDerivation respectively, obtains：
Solve equation availableWithValue, such as following formula：
In formula：
α, the value of β are can obtain after trying to achieve parameter k and b：
α=eb, β=k
Embodiment：
1st, simulation scale and data acquisition
(1) simulation scale
Road network scope comprising Hangzhou around field at a high speed within and lower sand form a team formed a team with Yuhang and the west of a city future technology City.The scope covers Hangzhou city Travel Range, it is sufficient to emulate the flow distribution of Hangzhou through street net.Road in simulation scale Net covering high speed, through street, secondary distributor road and branch road.
(2) data to parameter calibration are obtained
Uninterrupted and interrupted flows section is used to the data source and data volume of parameter calibration：
Flow.Each bicycle road flow conversions of whole section section 5min that RTMS is detected are whole section vehicle flowrate, And expand calculation for 1h flows.Unit：pcu/h.
Speed.The velocity amplitude of detection is one-to-one relationship with flow value.Unit：km/h.
Data volume.It is determined that intending the section demarcated, the section working day RTMS detection data of continuous month is chosen, be used to Parameter calibration.This chooses the RTMS data in 11 days~July 11 June in 2015.
(3) the OD data for model testing are obtained
Rely on the real-time mobile phone signaling data of Zhejiang movement, base station about 52000, mobile subscriber in research range of the present invention Quantity more than 5,000,000, average daily signaling bar number more than 500,000,000, every average daily signaling bar number 102 of mobile phone, the signaling of maximum ratio Type is location updating.Data acquisition time is the August of 4 days to 2015 July in 2015 4.
(4) the resident trip OD based on mobile phone signaling data estimates
Using cellphone subscriber as sample data, the OD distributions of permanent resident population are extrapolated according to the following formula.
In formula：OD is distributed for permanent resident population OD；Od is that the od drawn using mobile phone user data is distributed；a (average ownership) is the per capita ownership of cellphone subscriber, unit：Portion/people；p(penetration rate of Mobile phone) it is mobile phone permeability；M (market share) is the occupation rate of market of China Mobile；d(detection Probability) for mobile subscriber's mobile phone is detected probability.
2nd, parameter calibration
(1) speed-flow scatter diagram, calibrating parameters C, v are drawnf
According to RTMS speed, the corresponding relation of flow detection value, speed-flow scatter diagram is drawn, scatter diagram is intended Close, parameter can be calibrated：Section bicycle road traffic capacity C, section free stream velocity vf.Matched curve figure such as accompanying drawing 1 is to Fig. 4 It is shown；As shown in Table 1, a cutout parameter calibration result is as shown in Table 2 for continuous stream parameter calibration result.
(2) with least square method calibrating parameters α, β
Parameter alpha, β, continuous stream calibration result in demarcating continuous stream section and interruption flowpath segment BRP functions with least square method As shown in Table 1, interruption flows calibration result as shown in Table 2.
The continuous stream parameter calibration result of table one
Title Category of roads α β Traffic capacity C
Highway Highway 85 0.75 6 3500/2,4500/3,6000/4
Through street Through street 80 0.8 5.5 3500/2,4500/3,6000/4
The lower ring road of central strip Ring road 55 0.75 9 1600/1,3000/2
The lower ring road that center does not separate Ring road 55 0.75 9 1600/1,3000/2
The upper ring road of central strip Ring road 55 0.75 9 1600/1,3000/2
The upper ring road that center does not separate Ring road 55 0.75 9 1600/1,3000/2
Grade separation reversing loop Ring road 55 0.75 9 1600/1,3000/2
Stop parameter calibration result between table two
3rd, the inspection of roadlock model
Evening peak is obtained according to methods such as microwave bayonet socket data, manual research altogether in the range of the macroscopic artificial of Hangzhou at present The actual traffic data on flows of 626 sections of period.These actual measurement road data distribution of grades are in expressway/through street, trunk Road.Secondary distributor road and branch road.Each grade flow is in actual measurement total amount, and proportion is as shown in the table.
The simulating traffic of table three and measured discharge comparing result
From table three it can be found that under current simulation result, evening peak period, the actual measurement of 616 one-way road sections is total Amount and model are more or less the same, error 7.06%, and the flow accounting in each grade road is close with actual, a 6% maximum left side of error The right side, precision substantially conforms to require.

Claims (1)

1. a kind of traffic impedance computation method based on parameter calibration, it is characterised in that the method is comprised the following steps：
The acquisition of basic data and treatment, including obtain the section flow and speed for parameter calibration；Described section flow The whole section section detected by radio telephone mobile system each bicycle road flow conversion in 5 minutes is formed, described speed The velocity amplitude that degree is detected using radio telephone mobile system；
It is determined that intending the continuous stream or interruption flowpath segment demarcated, the section working day radio telephone of continuous month movement system is chosen The data that system is detected, to parameter calibration；
Build section Impedance Function model：
t = t f [ 1 + α ( V i C i ) β ]
In formula：T represents the link travel time between two intersections；tfWhen the expression volume of traffic is 0, the section between two intersections Running time；V represents section automobile traffic amount；C represents section actual capacity；α, β represent roadlock to be calibrated respectively Function parameter；
Parameter calibration is carried out to building section Impedance Function model, specifically：
S1, logarithmetics treatment is carried out to section Impedance Function：
l n ( t t f - 1 ) = β l n ( V C ) + l n α
For certain section LiHave：
t t f = L i v i L i v f = v f v i
Section Impedance Function after taking the logarithm is expressed as：
l n ( v f v - 1 ) = β l n ( V C ) + l n α
In formula：V represents the sections of road speed between two intersections, vfRepresent section free stream velocity；
Ifβ=k,Ln α=b, then have y=kx+b；Simple regression analysis is translated into ask Topic；
S2, calibrating parameters C, vf
According to speed, the corresponding relation of flow detection value, speed-flow scatter diagram is drawn, scatter diagram is fitted, calibrated Parameter C and vf；
S3, obtain parameter alpha to be calibrated, β using least square method.
CN201710263992.0A 2017-04-21 2017-04-21 A kind of traffic impedance computation method based on parameter calibration Expired - Fee Related CN106920395B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710263992.0A CN106920395B (en) 2017-04-21 2017-04-21 A kind of traffic impedance computation method based on parameter calibration

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710263992.0A CN106920395B (en) 2017-04-21 2017-04-21 A kind of traffic impedance computation method based on parameter calibration

Publications (2)

Publication Number Publication Date
CN106920395A true CN106920395A (en) 2017-07-04
CN106920395B CN106920395B (en) 2019-04-12

Family

ID=59567587

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710263992.0A Expired - Fee Related CN106920395B (en) 2017-04-21 2017-04-21 A kind of traffic impedance computation method based on parameter calibration

Country Status (1)

Country Link
CN (1) CN106920395B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107730920A (en) * 2017-10-23 2018-02-23 淮阴工学院 A kind of dynamically changeable lane control method based on spike nail light
CN110807926A (en) * 2019-11-05 2020-02-18 武汉理工大学 Road impedance prediction method and device based on hybrid traffic
CN111275965A (en) * 2020-01-20 2020-06-12 交通运输部科学研究院 Real-time traffic simulation analysis system and method based on internet big data
CN113111271A (en) * 2021-04-20 2021-07-13 智慧足迹数据科技有限公司 Travel OD data sample expansion method and device, computer equipment and storage medium
CN113393690A (en) * 2021-08-17 2021-09-14 长沙理工大学 Hybrid traffic balance distribution method considering automatic driving special lane

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070183479A1 (en) * 2006-01-17 2007-08-09 Airnet Communications Corporation Method to calibrate RF paths of an FHOP adaptive base station
CN101930670A (en) * 2010-08-12 2010-12-29 东南大学 Method for predicting social vehicle running time on bus travel road section
JP5100845B2 (en) * 2008-10-28 2012-12-19 株式会社パスコ Road measuring device and road measuring method
CN104318773A (en) * 2014-11-04 2015-01-28 杭州市综合交通研究中心 Traffic jam determining method based on traffic jam space-time total amount
CN106327871A (en) * 2016-09-06 2017-01-11 华南理工大学 Highway congestion forecasting method based on historical data and reservation data

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070183479A1 (en) * 2006-01-17 2007-08-09 Airnet Communications Corporation Method to calibrate RF paths of an FHOP adaptive base station
JP5100845B2 (en) * 2008-10-28 2012-12-19 株式会社パスコ Road measuring device and road measuring method
CN101930670A (en) * 2010-08-12 2010-12-29 东南大学 Method for predicting social vehicle running time on bus travel road section
CN104318773A (en) * 2014-11-04 2015-01-28 杭州市综合交通研究中心 Traffic jam determining method based on traffic jam space-time total amount
CN106327871A (en) * 2016-09-06 2017-01-11 华南理工大学 Highway congestion forecasting method based on historical data and reservation data

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107730920A (en) * 2017-10-23 2018-02-23 淮阴工学院 A kind of dynamically changeable lane control method based on spike nail light
CN110807926A (en) * 2019-11-05 2020-02-18 武汉理工大学 Road impedance prediction method and device based on hybrid traffic
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
CN113111271A (en) * 2021-04-20 2021-07-13 智慧足迹数据科技有限公司 Travel OD data sample expansion method and device, computer equipment and storage medium
CN113111271B (en) * 2021-04-20 2023-08-29 智慧足迹数据科技有限公司 Trip OD data sample expansion method and device, computer equipment and storage medium
CN113393690A (en) * 2021-08-17 2021-09-14 长沙理工大学 Hybrid traffic balance distribution method considering automatic driving special lane

Also Published As

Publication number Publication date
CN106920395B (en) 2019-04-12

Similar Documents

Publication Publication Date Title
CN106920395B (en) A kind of traffic impedance computation method based on parameter calibration
CN104217250B (en) A kind of urban rail transit new line based on historical data opens passenger flow forecasting
CN103745089A (en) Multi-dimensional public transport operation index evaluation method
CN103440411B (en) Set up based on exposed population group area the method for traffic noise pollution model of acoustic environment functional areas
CN106327871B (en) A kind of crowded prediction technique of highway of fusion historical data and reservation data
CN103699775B (en) Urban road traffic guidance strategy automatic generation method and system
CN105070057B (en) Method for monitoring real-time road conditions of road
CN104298540B (en) A kind of underlying model parameter correcting method of traffic simulation software
CN106683404A (en) Method for obtaining bus passenger flow OD (origin-destination) matrix through mobile phone location technology
CN109686091B (en) Traffic flow filling algorithm based on multi-source data fusion
CN101281684A (en) Area traffic control system for synergism operation of inducement and zone control of display panel
CN102968901A (en) Method for acquiring regional congestion information and regional congestion analyzing device
CN107767011B (en) Pedestrian characteristic acquisition system and service level dynamic evaluation method for railway station
CN106898142B (en) A kind of path forms time reliability degree calculation method considering section correlation
Rosik et al. Regional dispersion of potential accessibility quotient at the intra-European and intranational level. Core-periphery pattern, discontinuity belts and distance decay tornado effect
CN105678034A (en) Pedestrian crossing time model and signalized intersection pedestrian crosswalk width optimizing method
Sun et al. Impacts of rainfall and catchment characteristics on bioretention cell performance
Liu et al. Recent practices in regional modeling of nonmotorized travel
McAndrews et al. How do the definitions of urban and rural matter for transportation safety? Re-interpreting transportation fatalities as an outcome of regional development processes
CN106571032A (en) OD calibration method using mobile phone signaling big data and dynamic traffic assignment
CN105551241B (en) A kind of real-time jamming analysis method based on FCD and EP multi-source datas
CN106803227A (en) A kind of public transport OD analysis methods of the main volume of the flow of passengers based on public transport vehicle-mounted gps data
Xiong et al. Spatial heterogeneity in the effects of built environments on walking distance for the elderly living in a mountainous city
CN110889086B (en) Method for evaluating influence of shared electric leasing automobile on urban automobile exhaust emission
CN106530700A (en) Method for judging traffic congestion based on fixed detector

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190412

Termination date: 20200421