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

Road traffic accident multi-happening section identifying method Download PDF

Info

Publication number
CN100454354C
CN100454354C CNB200610113504XA CN200610113504A CN100454354C CN 100454354 C CN100454354 C CN 100454354C CN B200610113504X A CNB200610113504X A CN B200610113504XA CN 200610113504 A CN200610113504 A CN 200610113504A CN 100454354 C CN100454354 C CN 100454354C
Authority
CN
China
Prior art keywords
accident
road
traffic
accident rate
critical value
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.)
Expired - Fee Related
Application number
CNB200610113504XA
Other languages
Chinese (zh)
Other versions
CN1932923A (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

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
Figure C20061011350400051
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 } - { 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
Figure C20061011350400062
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

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 γ 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,
γ-high accident number of times critical value,
λ-Parameter for Poisson Distribution, K is a nonnegative integer,
Calculate high accident rate critical value R, be shown below:
P { Y &GreaterEqual; R } = 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
Wherein: the Y-accident rate;
The high accident rate critical value of R-;
N is a nonnegative integer, R=4;
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 all greater than level of signifiance α corresponding down critical value γ and R, 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 CN1932923A (en) 2007-03-21
CN100454354C true 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)

Families Citing this family (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
CN102402715A (en) * 2010-09-13 2012-04-04 方正国际软件有限公司 Method and device for presenting incident hotspot region
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
CN103366551A (en) * 2012-04-05 2013-10-23 郭海锋 Road traffic safety evaluation method
CN102855395B (en) * 2012-08-16 2015-05-20 长安大学 Method for distinguishing road black spot
CN103955596B (en) * 2014-03-14 2017-09-22 安徽科力信息产业有限责任公司 A kind of accident focus synthetic determination method based on traffic accident acquisition technique
CN104392076B (en) * 2014-12-16 2017-10-31 东南大学 City road network pedestrian traffic accident black spot recognition method
CN106935030A (en) * 2017-03-31 2017-07-07 青岛海信网络科技股份有限公司 A kind of expressway safety hidden danger section recognition methods and device
CN106875687B (en) * 2017-04-24 2020-01-14 哈尔滨工业大学 Sliding window method-based automatic identification method for accident multi-occurrence road sections
CN115240407B (en) * 2022-06-10 2024-01-12 深圳市综合交通与市政工程设计研究总院有限公司 Method and device for identifying black spots of traffic accidents, electronic equipment and storage medium

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
109国道宁夏石中段交通事故多发点研究. 胡江碧,刘璟,彭利人,任福田,刘小明.北京工业大学学报,第29卷第3期. 2003
109国道宁夏石中段交通事故多发点研究. 胡江碧,刘璟,彭利人,任福田,刘小明.北京工业大学学报,第29卷第3期. 2003 *
基于数理统计原理的交通事故多发点识别. 姜华平,许洪国.济南交通高等专科学校学报,第9卷第3期. 2001
基于数理统计原理的交通事故多发点识别. 姜华平,许洪国.济南交通高等专科学校学报,第9卷第3期. 2001 *
道路交通事故多发点质量控制鉴别法的改进. 裴玉龙.哈尔滨工业大学学报,第38卷第1期. 2006
道路交通事故多发点质量控制鉴别法的改进. 裴玉龙.哈尔滨工业大学学报,第38卷第1期. 2006 *

Also Published As

Publication number Publication date
CN1932923A (en) 2007-03-21

Similar Documents

Publication Publication Date Title
CN100454354C (en) Road traffic accident multi-happening section identifying method
WO2020244288A1 (en) Method and apparatus for evaluating truck driving behaviour based on gps trajectory data
CN102360525B (en) Discriminant analysis-based high road real-time traffic accident risk forecasting method
CN103077610A (en) Road trip time estimating method and system
Cafiso et al. Driving performance, alignment consistency, and road safety: real-world experiment
CN103903433A (en) Real-time dynamic judgment method and device for road traffic state
CN108922168B (en) A kind of mid-scale view Frequent Accidents road sentences method for distinguishing
CN104408925A (en) Array radar based intersection running state evaluation method
CN101465059A (en) Early-warning system for inspecting and judging urban road traffic safety situation
CN113870564B (en) Traffic jam classification method and system for closed road section, electronic device and storage medium
Afghari et al. Bayesian latent class safety performance function for identifying motor vehicle crash black spots
CN103593973A (en) Urban road traffic situation assessment system
CN106526149B (en) A kind of Pavement Condition prediction technique based on be open to traffic duration and the volume of traffic
CN100481153C (en) Method for automatically inspecting highway traffic event based on offset minimum binary theory
CN108922170B (en) Intersection safety evaluation method based on electronic police snapshot data
CN113436432A (en) Method for predicting short-term traffic risk of road section by using road side observation data
CN103971519A (en) System and method of using traffic conflicts for judging accident-prone sections
CN103778782B (en) A kind of traffic behavior division methods based on semi-supervised learning
Tarko et al. Estimating the frequency of crashes as extreme traffic events
Lim et al. Identifying High-Crash-Risk Intersections: Comparison of Traditional Methods with the Empirical Bayes–Safety Performance Function Method
CN104361349A (en) Car inspection device and toll data fusion based abnormal traffic state identification method and system
CN105894325A (en) Method for quantitatively analyzing power market index system
CN104809869A (en) Viaduct entrance ramp traffic state identifying method based on lifting height limiting rack
CN104050730A (en) Method and apparatus for processing vehicle reported information
Das et al. Investigating the safety impact of raised pavement markers on freeways in Louisiana

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