CN101826258B - Method for predicting simple accidents on freeways - Google Patents

Method for predicting simple accidents on freeways Download PDF

Info

Publication number
CN101826258B
CN101826258B CN2010101458112A CN201010145811A CN101826258B CN 101826258 B CN101826258 B CN 101826258B CN 2010101458112 A CN2010101458112 A CN 2010101458112A CN 201010145811 A CN201010145811 A CN 201010145811A CN 101826258 B CN101826258 B CN 101826258B
Authority
CN
China
Prior art keywords
section
curve
traffic
vertical
highway
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
CN2010101458112A
Other languages
Chinese (zh)
Other versions
CN101826258A (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 CN2010101458112A priority Critical patent/CN101826258B/en
Publication of CN101826258A publication Critical patent/CN101826258A/en
Application granted granted Critical
Publication of CN101826258B publication Critical patent/CN101826258B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention relates to a method for predicting simple accidents on freeways, which can be used for predicting the occurrence number of traffic accidents on some road sections of freeways and belongs to the field of traffic safety. The traffic accidents are not only concerned with speeds and speed differences, but also are closely associated with traffic volume and horizontal and longitudinal line shapes. The invention establishes relational models comprising intersection angles of the accidents and horizontal curves, gradients of longitudinal slopes and the like by utilizing the method of statistical regression. Practice proves that the accident number predicted with the method can be preferably coincident with the actually occurred accident number, and therefore, the method provides reliable theories and basis for reforming accident areas, preventing the occurrence of the traffic accidents, and reducing the severities of accidents.

Description

Method for predicting simple accidents on freeways
Technical field
The present invention is a kind of method for predicting simple accidents on freeways, can be used to predict the generation number of certain road section traffic volume accident of highway, belongs to the traffic safety field.
Background technology
Statistics shows that highway all is the highest highway of safe coefficient, and on the contrary in China in worldwide.According to Japan's statistics, the accident rate of its highway per 100 km only is 1/2 or 1/3 of a common road, and the magnitude of traffic flow is about 10 times of common road; The statistics of the U.S. is that the accident rate of its highway per 100 km only is 1/10 of a common road, also is far superior to Japan.And in China, the magnitude of traffic flow of highway is basic to similar abroad, even is lower than a lot of countries, but according to Public Security Department of China Ministry of Public Security statistics, the per 100 km accident rate but is more than 4 times of common road.If factors such as the volume of traffic are taken into account, the traffic safety status of China's highway will be more troubling again.Therefore, be necessary to set up the perfect traffic accidents Forecasting Methodology of a cover, for the prediction and the prevention of China's traffic accidents provides theoretical foundation.
According to the data of American National major accident research institute (NCSS), traffic mortality and travelling speed gradient (Δ V)) biquadratic be directly proportional.Approximate function is as follows:
Death = ( ΔV 114.24 ) 4 - - - ( 1 )
In the formula: Death-traffic mortality, %;
Δ V-running velocity gradient, km/h;
People such as Pei Yulong of Harbin Institute of Technology carry out regretional analysis to vehicle speed standard difference and hundred million truck kilometer accident rates on seven highways of China, have obtained the relational model of hundred million truck kilometer accident rates and vehicle speed standard deviation.Model shows: the speed of a motor vehicle distributes discrete more, and accident rate is high more, and model is as follows:
AR=9.583e 0.0553σ(2)
In the formula: AR---hundred million truck kilometer accident rates;
The standard deviation of σ---the speed of a motor vehicle (km/h).
The Du Boying of Tongji University carries out on the basis of analysis-by-synthesis at accident rate and the operating speed to Foreign Expressway, set forth on the highway traffic accident prediction method based on operating speed.
The computing formula of highway mortality ratio is as follows:
I Death = [ V ( 2 E - 09 ) × V 3.667 - ( 1 - V / V ‾ ) ] · [ ΔV 114.24 ] 4 - - - ( 3 )
I in the formula Death---traffic mortality (inferior/10 6Veh*km);
V---travelling speed (km/h);
Δ V---velocity gradient, the i.e. difference of section operating speed and average operating speed (km/h);
Figure GSA00000084388500022
---average operating speed (km/h).
As far back as 1964, Solomon studied the speed of a motor vehicle and the relation of safety with regard to beginning, and the relation between speed and the accident is discussed mostly, does not consider the influence of the volume of traffic.Make a general survey of both at home and abroad, in the influence research to speed and safety, major part all is in the relation of discussing between travelling speed, velocity contrast and the accident, not have to quantize the variation such as the volume of traffic, so conclusion and inequality have stronger cogency.Because the special national conditions and the highway of China can not directly be indiscriminately imitated external speed of a motor vehicle theoretical model, also need to carry out concrete theoretical research and engineering practice, constantly sum up the speed of a motor vehicle safety theory of China's highway, propose reasonably to solve countermeasure.
Summary of the invention
The present invention has taken all factors into consideration the influence to accident of horizontal curve factor, vertical section factor, the volume of traffic on the basis that the domestic communication operation characteristic is analyzed, and has obtained effect preferably through actual verification.
Use the brief traffic hazard frequency of formula prediction highway step:
Road is a three-dimensional entity.General said route is meant the locus of center line of road.The projection of route on surface level is called the planimetric map of route.Vertically dissecing again along center line, the row expansion then is the skiagraph of route.The turning point of two slope sections is relaxed with one section curve for the ease of driving on the vertical section, is called vertical curve.
(1) determine prediction highway condition, collect road information, comprising:
The average corner of horizontal curve:
Ave _ angle = Σ | α i | n - - - ( 4 )
The parametric representation meaning:
a iThe corner of i bar horizontal curve in the expression highway section;
N represents to comprise in the highway section horizontal curve number;
Horizontal curve figure as shown in Figure 1
2. vertical curve factor:
A basic variable supposing every vertical curve is V (i), and its unit is the variation of every 100m vertical curve gradient.
Figure GSA00000084388500031
V ( i ) = | g i - g i + 1 | L xi - - - ( 5 )
g iThe expression gradient, g i=tan θ;
The weight of vertical curve i:
Figure GSA00000084388500033
VC: the slope change value after the weighting
VC=∑ iWV(i)×V(i) (6)
The weighting gradient of vertical section:
Ave_slope=∑ iWG(k)×|g k|(7)
In the highway section on the k/weight in descending highway section
Figure GSA00000084388500034
g k: the gradient in k highway section, highway section
Vertical curve figure as shown in Figure 2
3. cart number percent is formed in traffic;
(2) calculate exposure variable EXPO
EXPO=AADT*365*L*10 -6*Y
The parametric representation meaning:
AADT: annual average daily traffic
L: road section length
Y: prediction continues the time
(3) gather the cart ratio, utilize Stata9.0 software to carry out statistical study, adopt the regression analysis of rejecting backward: set up full model earlier, judgement according to index of correlation (| Z| value minimum) among the output result, one of each rejecting least meets the variable that enters model, no longer contains till the independent variable that does not meet criterion in regression equation.Therefore, remove incongruent independent variable successively and return again, obtain brief hazard model at last:
λ i=EXPO·EXP(-2.676614+0.0071095·Ave_angle+0.737331·VC+0.2539619·Ave?slope+6.14963·Truck)(8)
The parametric representation meaning:
λ i: i section prediction accident number;
EXPO: expose variable;
Y: prediction continues the time;
The L:i road section length;
The average corner of horizontal curve in the Ave_angle:i highway section;
VC:: vertical curve index, the slope change value after the weighting;
Average_slope: vertical curve index, the weighting gradient of vertical section;
Truck: cart ratio.
The present invention has taken all factors into consideration the influence to accident of horizontal curve factor, vertical section factor, the volume of traffic on the basis that the domestic communication operation characteristic is analyzed, and has obtained effect preferably through actual verification.
Description of drawings:
Fig. 1 horizontal curve synoptic diagram; Fig. 2 vertical curve synoptic diagram; Fig. 3 horizontal curve example schematic; Fig. 4 vertical curve example schematic.
Embodiment:
As Fig. 3, Figure 4 shows that the flat vertical instance graph of a certain highway, this highway annual average daily traffic is 2500, to the annual investigation statistics of this highway, the cart proportion is 13%.
The first step: calculate the average corner of horizontal curve
Ave _ angle = Σ | α i | n = 5 18 π + 1 3 π 2 = 11 36 π
Second step:
Vertical curve L 1Length=10000*[2%-(3%)]=500m
Vertical curve L 2Length=8000*4.5%=360m
This road section length=150+500+300+360+280=1590m
V ( 1 ) = V ( j ) = | g j - g j + 1 | L j = 2 % - ( - 3 % ) 5 = 1 %
Figure GSA00000084388500043
V ( 2 ) = V ( i ) = | g i - g i + 1 | L i = 4.5 % 3.6 = 1.25 %
Figure GSA00000084388500051
VC: the slope change value after the weighting
VC=∑ iWV(i)×V(i)=V(1)*WV(1)+V(2)*WV(2)=0.54%
Longitudinal gradient i in the highway section 1=-3% shared weight:
Figure GSA00000084388500052
Longitudinal gradient i in the highway section 2=2% shared weight:
Figure GSA00000084388500053
Longitudinal gradient i in the highway section 3=-2.5% shared weight:
Figure GSA00000084388500054
The weighting gradient of vertical section:
Ave_slope=∑ jWG(k)×|g k|=g1*WG(1)+g2*WG(2)+g3*WG(3)=3%*0.094+2%*0.189+2.5%*0.176=0.011
The 3rd step: calculate cart number percent
Know by condition: Truck=13%
The 4th step: expose variable and calculate
EXPO=AADT*365*L*10 -6*Y=2500*365*3*10 -6*1=2.74
The 5th step: accident number prediction
λ i=EXPO·EXP(-2.676614+0.0071095·Ave_angle+0.737331·VC+0.2539619·Ave_slope+6.14963·Truck%)
=2.74*EXP(-2.676614+0.0071095*0.96+0.737331*0.0054+0.2539619*0.011+6.14963*13%
=2.74*0.155
=1
Use respectively the Chongqing Chengdu-Chongqing at a high speed, the upper bound at a high speed, long ten thousand at a high speed, Yunnan sieve (gateway of the village) rich (rather) at a high speed, the road traffic of many highways such as the Beijing-Tianjin pool is according to top, the simple accidents on freeways forecast model verified the result is as showing
Highway Road section length (km) The volume of traffic (/ day) Cart ratio (%) The prediction duration (year) The prediction traffic hazard The actual traffic accident
Luo Fu (K40+791~K41+390) 1 4300 46.2% 1 5 5
Luo Fu (k38+318~K40+39) 1.720 4300 46.2% 1 8 7
Chengdu-Chongqing (k8+500~K9+764) 1.264 12473 48% 2 7 6
Chengdu-Chongqing (k2+100~K3+456) 1.56 17506 48% 2 4 4
Long by ten thousand (K152+327~K154+772) 2.444 2310 32% 2 8 9
Long by ten thousand (K182+108~K183+507) 1.398 1500 32% 2 10 11
The upper bound (K26+870~K27+900) 1.03 10748 15.7% 2 17 15
The upper bound (K24+789~K26+930) 2.414 30380 25.7% 2 29 25
The Beijing-Tianjin pool (K6~K8) 2 18169 27.44% 4 66 64
The Beijing-Tianjin pool (K49~K56) 7 19400 43.11% 4 25 20
As can be seen from the table, prediction accident number and actual accidents number can reasonablely coincide, thereby the trouble-saving theoretical foundation that provides is provided.

Claims (1)

1. method for predicting simple accidents on freeways is characterized in that carrying out according to the following steps:
Route is meant the locus of center line of road; The projection of route on surface level is called the planimetric map of route; Vertically dissecing again along center line, the row expansion then is the skiagraph of route; The turning point of two slope sections is relaxed with one section curve for the ease of driving on the vertical section, is called vertical curve;
(1) determine prediction highway condition, collect road information, comprising:
1.1 the average corner of horizontal curve:
Ave _ angle = Σ | α i | n
The parametric representation meaning:
α iThe corner of i bar horizontal curve in the expression highway section;
N represents to comprise in the highway section horizontal curve number;
1.2. vertical curve factor:
A basic variable supposing every vertical curve is V (i), and its unit is the variation of every 100m vertical curve gradient;
Figure FSB00000591028000012
V ( i ) = | g i - g i + 1 | L xi
g iThe expression gradient, g i=tan θ;
The weight of vertical curve i:
VC: the slope change value after the weighting
VC=∑ iWV(i)×V(i)
The weighting gradient of vertical section:
Ave_slope=∑ iWG(k)×|g k|
In the highway section on the k/weight in descending highway section
Figure FSB00000591028000015
g k: the gradient in k highway section, highway section
1.3. traffic is formed, cart number percent;
(2) calculate exposure variable EXPO
EXPO=AADT*365*L*10 -6*Y
The parametric representation meaning:
AADT: annual average daily traffic
L: road section length
Y: prediction continues the time
(3) utilize Stata9.0 software to carry out statistical study, adopt the regression analysis of rejecting backward: obtain brief hazard model at last:
λ i=EXPO·EXP(-2.676614+0.0071095·Ave_angle+0.737331·VC
+0.2539619·Ave_slope+6.14963·Truck)
The parametric representation meaning:
λ i: i section prediction accident number;
EXPO: expose variable;
Y: prediction continues the time;
The L:i road section length;
The average corner of horizontal curve in the Ave_angle:i highway section;
VC: vertical curve index, the slope change value after the weighting;
Average_slope: vertical curve index, the weighting gradient of vertical section;
Truck: the cart ratio is obtained by historical data.
CN2010101458112A 2010-04-09 2010-04-09 Method for predicting simple accidents on freeways Expired - Fee Related CN101826258B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010101458112A CN101826258B (en) 2010-04-09 2010-04-09 Method for predicting simple accidents on freeways

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010101458112A CN101826258B (en) 2010-04-09 2010-04-09 Method for predicting simple accidents on freeways

Publications (2)

Publication Number Publication Date
CN101826258A CN101826258A (en) 2010-09-08
CN101826258B true CN101826258B (en) 2011-12-07

Family

ID=42690162

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010101458112A Expired - Fee Related CN101826258B (en) 2010-04-09 2010-04-09 Method for predicting simple accidents on freeways

Country Status (1)

Country Link
CN (1) CN101826258B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102103715A (en) * 2010-11-18 2011-06-22 上海海事大学 Negative binomial regression-based maritime traffic accident investigation analysis and prediction method
CN105096586B (en) * 2014-05-15 2018-06-05 交通运输部公路科学研究所 Expressway in Plain accident forecast method based on traffic flow character parameter
CN105701579A (en) * 2016-03-08 2016-06-22 北京工业大学 Prediction method for predicting traffic accidents on basic section of dual-lane secondary road in plateau area
CN105608902B (en) * 2016-03-28 2017-10-03 辽宁省交通科学研究院 A kind of expressway traffic accident Multiple trauma identification system and method
CN106530720B (en) * 2016-12-28 2018-12-14 吉林大学 A kind of identification of Expressway Road traffic safety stain section and method for early warning
CN107403053A (en) * 2017-08-07 2017-11-28 中建隧道建设有限公司 Measurement and Data Processing method and device for shield-tunneling construction
CN108320514A (en) * 2018-04-04 2018-07-24 哈尔滨工业大学 The analysis method that accident rate is influenced based on the Tobit highway route indexs returned
CN110264724B (en) * 2019-07-09 2020-08-25 哈尔滨工业大学 Interactive highway traffic accident prediction method
CN112419719B (en) * 2020-11-18 2022-06-07 济南北方交通工程咨询监理有限公司 Method and system for evaluating traffic operation safety of highway
CN114627643B (en) * 2022-02-28 2023-07-14 青岛海信网络科技股份有限公司 Highway accident risk prediction method, device, equipment and medium
CN116959248B (en) * 2023-07-22 2024-04-05 同济大学 Expressway accident risk optimization method based on interactive self-interpretation model

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101016052A (en) * 2007-01-25 2007-08-15 吉林大学 Warning method and system for preventing deviation for vehicle on high standard highway

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004094444A (en) * 2002-08-30 2004-03-25 Tokio Marine Research Institute Information processing method for preventing traffic accident

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101016052A (en) * 2007-01-25 2007-08-15 吉林大学 Warning method and system for preventing deviation for vehicle on high standard highway

Also Published As

Publication number Publication date
CN101826258A (en) 2010-09-08

Similar Documents

Publication Publication Date Title
CN101826258B (en) Method for predicting simple accidents on freeways
Ehsani et al. Modeling of vehicle fuel consumption and carbon dioxide emission in road transport
CN105608902B (en) A kind of expressway traffic accident Multiple trauma identification system and method
CN103413263B (en) Travel time index entropy traffic circulation evaluation method
CN101246514B (en) City fast road intercommunicated overpass simulation design system and method for establishing design model
CN101976501B (en) Principal component analysis and neural network based port road safety prediction method
CN103606269B (en) A kind of control method improving traffic efficiency of freeway construction area
CN104732075A (en) Real-time prediction method for urban road traffic accident risk
CN101246513A (en) City fast road intercommunicated overpass simulation design system and selection method
CN101930670A (en) Method for predicting social vehicle running time on bus travel road section
CN103234544A (en) Methods for building power consumption factor model and estimating following-up driving range of electric car
CN109543882B (en) Bus network density calculation method based on optimal bus average station spacing
Lai et al. Optimization of rail transit alignments considering vehicle dynamics
CN102201021A (en) Expressway aided design system
CN103632535B (en) Judgment method for section pedestrian crossing signal lamp arrangement
CN104575050A (en) Express way ramp intelligent inducing method and device based on floating vehicles
CN103235861A (en) Anticipation and evaluation method for implementation of passenger and cargo separating road
CN103646374B (en) A kind of routine bus system running comfort index calculation method
CN104091440B (en) Based on the urban traffic blocking Cost Estimation of category of roads and type of vehicle
Garber et al. Incorporating crash risk in selecting congestion-mitigation strategies: Hampton Roads area (Virginia) case study
CN103942950A (en) Method for predicating degree of reliability of traffic circulation of arterial highway under snow and ice environments
CN106067078A (en) Bus berth distribution optimization method for double platforms
Jung et al. Evaluation of lateral driving safety for elderly drivers in high-speed curve driving conditions
CN109190935B (en) Highway charging station planning method considering service area and vehicle accident
Barnes et al. Impact of peak and off-peak tolls on traffic in san francisco–oakland bay bridge corridor in california

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent for invention or patent application
CB03 Change of inventor or designer information

Inventor after: He Yulong

Inventor after: Zhong Liande

Inventor after: Sun Xiaoduan

Inventor after: Chen Yongsheng

Inventor after: Liu Xiaoming

Inventor after: Zhang Jie

Inventor after: Sun Mingling

Inventor after: Xu Ting

Inventor after: Lian Jia

Inventor before: He Yulong

Inventor before: Wang Yidai

Inventor before: Lian Jia

Inventor before: Wang Chao

Inventor before: Zhong Liande

Inventor before: Sun Xiaoduan

Inventor before: Chen Yongsheng

Inventor before: Bao Xingjian

Inventor before: Xu Ting

Inventor before: Zhang Jie

Inventor before: Hou Shuzhan

Inventor before: Wang Huarong

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: HE YULONG ZHONG LIANDE SUN XIAODUAN CHEN YONGSHENG BAO XINGJIAN XU TING ZHANG JIE HOU SHUZHAN WANG HUARONG WANG YIYI LIAN JIA WANG CHAO TO: HE YULONG ZHONG LIANDE SUN XIAODUAN CHEN YONGSHENG LIU XIAOMING ZHANG JIE SUN MINGLING XU TING LIAN JIA

C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20111207

Termination date: 20130409