CN1932923A - Road traffic accident multi-happening section identifying method - Google Patents

Road traffic accident multi-happening section identifying method Download PDF

Info

Publication number
CN1932923A
CN1932923A CNA200610113504XA CN200610113504A CN1932923A CN 1932923 A CN1932923 A CN 1932923A CN A200610113504X A CNA200610113504X A CN A200610113504XA CN 200610113504 A CN200610113504 A CN 200610113504A CN 1932923 A CN1932923 A CN 1932923A
Authority
CN
China
Prior art keywords
accident
road
traffic
frequency
signifiance
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
CNA200610113504XA
Other languages
Chinese (zh)
Other versions
CN100454354C (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.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CNB200610113504XA priority Critical patent/CN100454354C/en
Publication of CN1932923A publication Critical patent/CN1932923A/en
Application granted granted Critical
Publication of CN100454354C publication Critical patent/CN100454354C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Traffic Control Systems (AREA)

Abstract

The invention is concerned with the distinguishing method for the road with more traffic accident. The invention uses double target of alpha-gamma to distinguish the road with more traffic accident, assigns the remarkable level alpha, according to the Poisson probability formula or the frequency accumulation chart, after computes the critical value gamma of the high accident number and the critical value R of the high accident rate, compares the gamma and the R of the un-distinguish road to get the distinguishing result of the road with more traffic accident.

Description

Road traffic accident multi-happening section identifying method
Technical field
The present invention relates to the traffic safety field, more relate to a kind of complete reliable road traffic accident multi-happening section identifying method in this field.
Background technology
The first step in road safety improvement project also is a step of most critical, is exactly the position of determining the accident-prone road section of urgent need improvement.Present most widely used method is the casualty data statistical analysis method, and different according to casualty data statistical treatment discriminant criterion and method of discrimination be divided into various method again, but these methods all have separately applicable elements and deficiency.
The absolute number method is divided into accident frequency method and accident rate method again.The accident frequency method is not considered factors such as the volume of traffic, severity of injuries, uses separately the reflection accident situation is had one-sidedness, needs a large amount of and long investigation statistics; Simultaneously the accident number of times how much lack theoretic tightness, dividing the accident-prone road section, influenced by artificial subjective factor bigger.The accident rate method has been considered the influence of the volume of traffic and two key factors of road section length on the basis of accident frequency method, but does not find the foundation of tight division accident-prone road section yet in theory.
The probability statistics distribution is theoretical tight, but practical condition harshness.The one, the accident number that counts by average road section length is on some roads and do not meet Poisson distribution, therefore can not be applied widely; In addition, can not determine a numerical value stereotypedly, consider the difference that the whole nation is regional, regional.
The cumulative frequency method is not considered the order of severity of volume of traffic factor and accident, though the catastrophe point with the cumulative frequency matched curve has certain mathematical theory foundation as the standard that defines the accident-prone road section simultaneously, but this method only is suitable for the similar matched curve of road conditions has the situation of flex point (sudden change), can't explain its practical significance in practice.
Summary of the invention
The purpose of this invention is to provide a kind of road traffic accident multi-happening section identifying method, it is more complete, reliable to make road traffic accident multi-happening section differentiate.
The present invention is a kind of discrimination method of road traffic accident multi-happening section, it is characterized in that: it may further comprise the steps:
1) utilize the checkout equipment that is arranged on the road to detect road traffic accident number and the traffic data of collecting highway section to be identified, the traffic hazard number is meant that highway section to be identified is in a measurement period and the corresponding traffic hazard number of times of highway mileage number, with the volume of traffic of traffic hazard number of times, be the road accident rate in this highway section divided by this highway section;
2) utilize computing machine to read road traffic accident number and traffic data, utilize statistical analysis software differentiation traffic hazard number and accident rate whether to meet Poisson distribution, can utilize the Poisson distribution Chi-square Test module (Chi-square) among the existing SAS of statistical software to realize;
3) when traffic hazard number and accident rate meet Poisson distribution, given level of signifiance α calculates high accident number of times critical value γ and high accident rate critical value R under the level of signifiance α according to Poisson new probability formula (1);
P { X ≥ γ } = 1 - Σ K = 0 γ - 1 λ K K ! e - λ = Σ k = r ∞ λ K K ! e - λ ≤ α - - - ( 1 )
Wherein: X-accident frequency,
γ-accident number of times critical value,
λ-Parameter for Poisson Distribution, K is a nonnegative integer,
When traffic hazard number and accident rate do not meet Poisson distribution, statistical method according to routine, calculate the frequency of identical accident number and identical accident rate generation respectively, and according to accident number and accident rate from small to large order respectively with its frequency accumulative total addition, obtain accumulative total accident frequency and accumulative total accident rate frequency, behind the given level of signifiance α, cumulative frequency pairing accident number and the accident rate equivalent with (1-α) are critical value γ, R;
4) traffic hazard number of times and the road accident rate with highway section to be identified compares with level of signifiance α corresponding down critical value γ and R, if traffic hazard number of times and road accident rate are all greater than level of signifiance α corresponding down critical value γ and R, perhaps the two one of greater than the respective value under the level of signifiance α, judge that promptly this highway section is the accident-prone road section.
The described checkout equipment that is arranged on the road is unit, MetroCount5600 roadside or the vehicle magnetic reflection magnitude of traffic flow and analyser NC97.
The present invention uses the two indexs of α-γ to differentiate road traffic accident multi-happening section, and its beneficial effect is mainly reflected in:
1. when accident number and accident rate meet Poisson distribution, consider the influence of the volume of traffic simultaneously, avoided the subjectivity of only differentiating with the accident number of times, thereby avoided the highway section erroneous judgement of the low accident rate of high accident number of times is decided to be the accident-prone road section, strengthened discriminating simultaneously the low high accident rate of accident number of times highway section.
2 when accident number and accident rate do not meet Poisson distribution, and this method has been widened traffic hazard cumulative frequency method.
Description of drawings
Fig. 1 is an operational flowchart of the present invention;
Traffic hazard cumulative frequency figure among Fig. 2 embodiment 2.
Embodiment
Embodiment 1 accident number and accident rate meet the example of Poisson distribution
The traffic hazard number is meant highway section to be identified in a measurement period and the corresponding traffic hazard number of times of highway mileage number, is statistical unit with the accident number that takes place on the long 1Km highway section in this example, and road type is an expressway.The division of other road types is as shown in table 1.The corresponding second class road type urban road, when differentiating road traffic accident multi-happening section with " length less than L and the highway section of the number 〉=2 that had an accident " replacements " length 1Km highway section ", the constant multiple decision criteria of putting of urban road that is of other steps.L wherein is meant the city intersection spacing: subsidiary road 30m, major trunk roads 60m.
Table 1 road type
Type Roadway characteristic
1 2 3 4 Expressway comprises that highway, Class I highway, secondary expressway urban road do not comprise suburbs and counties under the jurisdiction of a large city's highway, it is characterized in that crossing density is big by two, Class III highway level Four, substandard highway
1) given level of signifiance α=0.05 time, the Poisson new probability formula of calculating high accident number of times critical value γ is as shown in Equation (2):
P { X &GreaterEqual; 3 } = 1 - P { X = 0 } - P { X = 1 } - P { X = 2 } = 1 - &Sigma; K = 0 2 ( 0.8 ) k K ! e 0.8 = 0.0474 < 0.05 - - - ( 2 )
Wherein: X-accident frequency;
γ-high accident number of times critical value;
λ 1-Parameter for Poisson Distribution, the total number of accident of the interior generation of institute's road type one-period of study obtains λ divided by the type road total kilometrage in the Free Region 1Estimated value.λ in this example 1Estimated value is 0.8.
So γ=3, promptly accident number of times critical value is 3.Then the highway section of all accident generation number 〉=3 just needs further to investigate the accident rate situation, thereby determines the accident-prone road section fully.
2) high accident rate critical value R is calculated in given α=0.01 time, as shown in Equation (3):
P { Y &GreaterEqual; 4 } = 1 - P { Y = 0 } - P { Y = 1 } - P { Y = 2 } - P { Y = 3 } = 1 - &Sigma; N = 0 3 ( 0.8 ) N N ! e 0.8 = 0.0091 < 0.01 - - - ( 3 )
Wherein: the Y-accident rate;
The high accident rate critical value of R-;
λ 2-Parameter for Poisson Distribution, the interior accident rate of institute's road type one-period of studying obtains λ divided by the type road total kilometrage in the Free Region 2Estimated value.
λ in this example 2Estimated value be 0.8.So R=4.
This example adopts for avoiding accident rate numerical value less than normal Form.
Accident number of times 〉=3 when highway section to be identified, accident rate 〉=4 o'clock, we just judge this accident-prone road section, highway section.Avoided the subjectivity of only differentiating with the accident number of times; Consider the influence of the volume of traffic simultaneously, avoided the highway section erroneous judgement of the low accident rate of high accident number of times is decided to be the accident-prone road section, strengthened discriminating simultaneously the low high accident rate of accident number of times highway section.
Embodiment 2 accident numbers and accident rate do not meet the example of Poisson distribution:
Do not meet the situation of Poisson distribution for accident number and accident rate, be divided into some little highway sections waiting to study road by every 1Km length, find out the path section number that the number that has an accident equates, be that a certain accident number is waiting to study the frequency that takes place on the road, try to achieve the frequency of this accident number with this accident frequency divided by total path section number, frequency is added up addition from small to large according to the accident number get traffic hazard and count cumulative frequency.Behind the given level of signifiance α, the cumulative frequency pairing accident number equivalent with (1-α) is critical value γ.The accident number is got accident rate divided by the volume of traffic, find out the path section number that accident rate equates, try to achieve the accident rate observing frequency, from small to large the observing frequency cumulative addition is got the road accident rate cumulative frequency according to accident rate divided by total path section number.Behind the given level of signifiance α, the cumulative frequency pairing accident rate equivalent with (1-α) is critical value R.
The accident number of domestic 109 national highways in Ningxia is analyzed the Poisson distribution that does not show on the probability statistics by statistics, its accident statistics and calculate as shown in table 2ly, and Fig. 2 is traffic hazard number of times frequency accumulative total figure.
Given level of signifiance α=0.05,1-α=0.95 is critical value γ with transverse axis cumulative frequency 95% pairing accident number of times promptly, from table 2 and Fig. 2, gets γ=15.
Given level of signifiance α=0.10,1-α=0.90 is critical value γ with transverse axis cumulative frequency 90% pairing accident number promptly, from table 2 and Fig. 2, gets γ=13.
The computing method of high accident rate critical value R are similar.
After the given level of signifiance 0.05, the highway section of the number of times that has an accident≤15 accounts for 95% of sum, thus these highway sections, accident number of times>15 as abnormity point, compare the position of finally definite accident-prone road section again with high accident rate critical value R.This method has been widened the cumulative frequency method, no longer is subject to the situation that the cumulative frequency curve has flex point.
Domestic 109 national highway accident statistics and the reckoners in table 2 Ningxia
Accident is counted X 0 1 2 3 4 5 6 7 8 9
Accident frequency f 8 19 20 41 31 20 18 13 10 8
Frequency (%) 3.3 7.9 8.3 17 12.9 8.3 7.5 5.4 4.1 3.3
Cumulative frequency (%) 3.3 11.2 19.5 36.5 49.4 57.7 65.2 71 75 78
Accident is counted X 10 11 12 13 14 15 16 17 18 19
Accident frequency f 10 7 4 4 6 6 4 4 4 4 ∑f=241
Frequency (%) 4.1 2.9 1.7 1.7 2.5 2.5 1.7 1.7 1.6 1.6
Cumulative frequency (%) 82.1 85 86.7 88.4 90.9 93.4 95.1 96.8 98.4 100

Claims (2)

1. the discrimination method of a road traffic accident multi-happening section is characterized in that, it may further comprise the steps:
1) utilizes the checkout equipment that is arranged on the road to detect road traffic accident number and the traffic data of collecting highway section to be identified,, obtain the road accident rate in this highway section the volume of traffic of traffic hazard number of times divided by this highway section;
2) utilize computing machine to read road traffic accident number and traffic data, and differentiate the traffic hazard number and whether accident rate meets Poisson distribution;
3) when traffic hazard number and accident rate meet Poisson distribution, given level of signifiance α calculates high accident number of times critical value γ and high accident rate critical value R under the level of signifiance α according to Poisson new probability formula (1);
P { X &GreaterEqual; &gamma; } = 1 - &Sigma; K = 0 &gamma; - 1 &lambda; K K ! e - &lambda; = &Sigma; k = r &infin; &lambda; K K ! e - &lambda; &le; &alpha; - - - ( 1 )
Wherein: X-accident frequency,
γ-accident number of times critical value,
λ-Parameter for Poisson Distribution, K is a nonnegative integer,
When traffic hazard number and accident rate do not meet Poisson distribution, calculate the frequency of identical accident number and identical accident rate generation respectively, and according to accident number and accident rate from small to large order respectively with its frequency accumulative total addition, obtain accumulative total accident frequency and accumulative total accident rate frequency, behind the given level of signifiance α, cumulative frequency pairing accident number and the accident rate equivalent with (1-α) are critical value γ, R;
4) during given level of signifiance α, the traffic hazard number of times and the road accident rate in highway section to be identified are compared with level of signifiance α corresponding down critical value γ and R, if traffic hazard number of times and road accident rate are all greater than level of signifiance α corresponding down critical value γ and R, perhaps the two one of greater than the respective value under the level of signifiance α, judge that promptly this highway section is the accident-prone road section.
2. the discrimination method of a kind of road traffic accident multi-happening section according to claim 1 is characterized in that the described checkout equipment that is arranged on the road is unit, MetroCount5600 roadside or the vehicle magnetic reflection magnitude of traffic flow and analyser NC97.
CNB200610113504XA 2006-09-29 2006-09-29 Road traffic accident multi-happening section identifying method Expired - Fee Related CN100454354C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB200610113504XA CN100454354C (en) 2006-09-29 2006-09-29 Road traffic accident multi-happening section identifying method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB200610113504XA CN100454354C (en) 2006-09-29 2006-09-29 Road traffic accident multi-happening section identifying method

Publications (2)

Publication Number Publication Date
CN1932923A true CN1932923A (en) 2007-03-21
CN100454354C CN100454354C (en) 2009-01-21

Family

ID=37878732

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB200610113504XA Expired - Fee Related CN100454354C (en) 2006-09-29 2006-09-29 Road traffic accident multi-happening section identifying method

Country Status (1)

Country Link
CN (1) CN100454354C (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101833610A (en) * 2010-04-09 2010-09-15 北京工业大学 Accident black-spot identification optimizing method
CN102034009A (en) * 2010-12-20 2011-04-27 东南大学 Equivalent road accident number method-based identifying equipment for accident-prone sections
CN102254430A (en) * 2011-06-03 2011-11-23 东南大学 Method for distinguishing accident-prone road section by using traffic conflicts
CN102360525A (en) * 2011-09-28 2012-02-22 东南大学 Discriminant analysis-based high road real-time traffic accident risk forecasting method
CN102402715A (en) * 2010-09-13 2012-04-04 方正国际软件有限公司 Method and device for presenting incident hotspot region
CN102855395A (en) * 2012-08-16 2013-01-02 长安大学 Method for distinguishing road black spot
CN103366551A (en) * 2012-04-05 2013-10-23 郭海锋 Road traffic safety evaluation method
CN103955596A (en) * 2014-03-14 2014-07-30 安徽科力信息产业有限责任公司 Accident hotspot comprehensive judging method based on traffic accident collection technology
CN104392076A (en) * 2014-12-16 2015-03-04 东南大学 Urban road network pedestrian traffic accident black spot identification method
CN106875687A (en) * 2017-04-24 2017-06-20 哈尔滨工业大学 A kind of accident-prone road section automatic identification method based on sliding window method
CN106935030A (en) * 2017-03-31 2017-07-07 青岛海信网络科技股份有限公司 A kind of expressway safety hidden danger section recognition methods and device
CN115240407A (en) * 2022-06-10 2022-10-25 深圳市综合交通与市政工程设计研究总院有限公司 Traffic accident black point identification method and device, electronic equipment and storage medium

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101833610A (en) * 2010-04-09 2010-09-15 北京工业大学 Accident black-spot identification optimizing method
CN102402715A (en) * 2010-09-13 2012-04-04 方正国际软件有限公司 Method and device for presenting incident hotspot region
CN102034009A (en) * 2010-12-20 2011-04-27 东南大学 Equivalent road accident number method-based identifying equipment for accident-prone sections
CN102034009B (en) * 2010-12-20 2012-07-18 东南大学 Equivalent road accident number method-based identifying equipment for accident-prone sections
CN102254430A (en) * 2011-06-03 2011-11-23 东南大学 Method for distinguishing accident-prone road section by using traffic conflicts
CN102360525B (en) * 2011-09-28 2013-10-16 东南大学 Discriminant analysis-based high road real-time traffic accident risk forecasting method
CN102360525A (en) * 2011-09-28 2012-02-22 东南大学 Discriminant analysis-based high road real-time traffic accident risk forecasting method
CN103366551A (en) * 2012-04-05 2013-10-23 郭海锋 Road traffic safety evaluation method
CN102855395A (en) * 2012-08-16 2013-01-02 长安大学 Method for distinguishing road black spot
CN103955596A (en) * 2014-03-14 2014-07-30 安徽科力信息产业有限责任公司 Accident hotspot comprehensive judging method based on traffic accident collection technology
CN103955596B (en) * 2014-03-14 2017-09-22 安徽科力信息产业有限责任公司 A kind of accident focus synthetic determination method based on traffic accident acquisition technique
CN104392076A (en) * 2014-12-16 2015-03-04 东南大学 Urban road network pedestrian traffic accident black spot identification method
CN106935030A (en) * 2017-03-31 2017-07-07 青岛海信网络科技股份有限公司 A kind of expressway safety hidden danger section recognition methods and device
CN106875687A (en) * 2017-04-24 2017-06-20 哈尔滨工业大学 A kind of accident-prone road section automatic identification method based on sliding window method
CN106875687B (en) * 2017-04-24 2020-01-14 哈尔滨工业大学 Sliding window method-based automatic identification method for accident multi-occurrence road sections
CN115240407A (en) * 2022-06-10 2022-10-25 深圳市综合交通与市政工程设计研究总院有限公司 Traffic accident black point identification method and device, electronic equipment and storage medium
CN115240407B (en) * 2022-06-10 2024-01-12 深圳市综合交通与市政工程设计研究总院有限公司 Method and device for identifying black spots of traffic accidents, electronic equipment and storage medium

Also Published As

Publication number Publication date
CN100454354C (en) 2009-01-21

Similar Documents

Publication Publication Date Title
CN1932923A (en) Road traffic accident multi-happening section identifying method
WO2020244288A1 (en) Method and apparatus for evaluating truck driving behaviour based on gps trajectory data
CN108550263B (en) Expressway traffic accident cause analysis method based on fault tree model
CN111460551B (en) Simulation, classification and compilation method of two-dimensional vehicle overrun load spectrum
CN102663252B (en) Combined type pavement usability performance evaluation method for underground road
CN114783183A (en) Monitoring method and system based on traffic situation algorithm
CN103413263A (en) Travel time index entropy traffic circulation evaluation method
CN110807919A (en) Urban road network traffic operation situation evaluation method based on vehicle passing data
CN103593973A (en) Urban road traffic situation assessment system
CN113935090B (en) Random traffic flow fine simulation method for bridge vehicle-induced fatigue analysis
Biswas et al. Assessment of level-of-service on urban arterials: a case study in Kolkata metropolis
CN108922170B (en) Intersection safety evaluation method based on electronic police snapshot data
CN113313357A (en) Traffic road safety evaluation method based on Gaussian process regression analysis
CN101075377A (en) Method for automatically inspecting highway traffic event based on offset minimum binary theory
CN102789690A (en) Illegal vehicle identifying method and system
CN109101734B (en) Prediction method for continuous rigid frame bridge downwarping risk
Lim et al. Identifying High-Crash-Risk Intersections: Comparison of Traditional Methods with the Empirical Bayes–Safety Performance Function Method
CN113919731A (en) Express way intersection area safety evaluation method based on traffic conflicts
CN115762140B (en) Expressway traffic accident risk prediction method considering variable heterogeneity
Das et al. Investigating the safety impact of raised pavement markers on freeways in Louisiana
CN116246491B (en) Track feature mining-based rear-end collision risk road section dynamic identification and control method
CN113920725B (en) Urban suburban main road accident prediction method based on mesoscopic level and storage medium
Coburn et al. Injury outcomes and costs for work zone crashes
Wang et al. Identification of risk sources of abnormal driving vehicles of expressway in port city
Vijay et al. Examining Crash Variable Based on Collision Type for PredictingCrash Severity on Urban Highways

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20090121

Termination date: 20140929

EXPY Termination of patent right or utility model