Final report for ITS Center project: Incident capacity estimation
UVA Center for Transportation Studies
A Research Project Report
For the National ITS Implementation Research Center
A U.S. DOT University Transportation Center
CHARACTERIZATION OF ACCIDENT CAPACITY REDUCTION Ling
Qin
Department of Civil Engineering
Email: lq9s@virginia.edu
Dr. Brian L. Smith
Department
of Civil Engineering
Email: briansmith@Virginia.EDU
Smart Travel Lab Report No. STL-2001-02University of Virginia
Center for Transportation Studies at the University of Virginia produces outstanding transportation professionals, innovative research results and provides important public service. The Center for Transportation Studies is committed to academic excellence, multi-disciplinary research and to developing state-of-the-art facilities. Through a partnership with the Virginia Department of Transportation’s (VDOT) Research Council (VTRC), CTS faculty hold joint appointments, VTRC research scientists teach specialized courses, and graduate student work is supported through a Graduate Research Assistantship Program. CTS receives substantial financial support from two federal University Transportation Center Grants: the Mid-Atlantic Universities Transportation Center (MAUTC), and through the National ITS Implementation Research Center (ITS Center). Other related research activities of the faculty include funding through FHWA, NSF, US Department of Transportation, VDOT, other governmental agencies and private companies.
Disclaimer: The contents of this report reflect the views of the authors, who are responsible for the facts and the accuracy of the information presented herein. This document is disseminated under the sponsorship of the Department of Transportation, University Transportation Centers Program, in the interest of information exchange. The U.S. Government assumes no liability for the contents or use thereof.
ABSTRACT
Incidents are a major cause of urban highway congestion. Incidents include any event that temporarily reduces roadway capacity, such as accidents, debris, disabled vehicles, and hazardous material spills. Incident capacity reduction will be used in the incident management systems, advanced traveler information systems, queuing analysis, and computer simulation models. A study conducted in 1970 estimated that an accident or disabled vehicle blocking one lane out of three lanes will reduce traffic flow by an average of 50 percent, an accident blocking two lanes out of three lanes will reduce traffic flow by an average of 79 percent, and an accident or disabled vehicle blocking shoulder lane(s) out of three lanes will reduce traffic flow by an average of 33 percent. However, very little other research has comprehensively addressed the impact of incidents on capacity. The premise of this project is that the incident capacity reduction is best modeled as a random variable, not a deterministic value, as is the current practice.
Extensive traffic flow and incident data for the Hampton Roads region of Virginia in the Smart Travel Laboratory provides us the opportunity to model incident capacity reduction as a random variable. Capacity under prevailing conditions can be estimated by calibrating a speed-flow and /or density-flow curve for a given highway. The peak of this curve defines capacity. When an incident occurs and a bottleneck is formed, the reduced capacity of the roadway is reached and can be measured directly as incident capacity. Incident capacity reduction can be computed as the difference between these two values over the capacity under prevailing conditions, and then modeled as a random variable.
This research focuses on estimating accident capacity reductions with one lane and two lanes out of three lanes blocked, and modeling them as random variables based on the traffic flow and accident data for the Hampton Roads region. The results indicate that accident capacity reduction with one lane out of three lanes blocked can be modeled as Beta distribution with an average of 63 percent, which is fairly higher than the result of previous research (50 percent), and accident capacity reduction with two lanes out of three lanes blocked can be modeled as Beta distribution with an average of 77 percent which is slightly lower than the result of previous research (79 percent).
Table of Contents
ABSTRACT ⅰ
TABLE OF CONTENTS ⅲ
LIST OF FIGURES ⅵ
LIST OF TABLES ⅷ
GLOSSARY ⅹ
ACKNOWLEDGEMENTS ⅺ
CHAPTER 1: INTRODUCTION 1
1.1.1 Overview 1
1.1.2 Importance for Queuing Analysis 2
1.1.3 Importance for Computer Simulation Models 5
CHAPTER 2: LITERATURE REVIEW 11
CHAPTER 3: METHODOLOGY 16
3.1 Definition of Capacity 16
3.2 Capacity Estimation Methods 17
3.2.1 Estimation of Capacity under Prevailing Conditions 17
3.2.2 Accident Capacity Estimation 21
3.3 Modeling Accident Capacity Reduction as a Random Variable 24
3.3.1 Overview 24
3.3.2 Activity Ⅰ: Hypothesizing Families of Distributions 24
3.3.2.1 Summary Statistics 25
3.3.2.2 Histogram 25
3.3.3 Activity Ⅱ: Estimating Parameters 27
3.3.4 Activity Ⅲ: Determining How Representative the Fitted Distributions Are 28
3.3.4.1 Heuristic Procedures 29
3.3.4.2 Goodness-of-Fit Test 30
3.4 Summary 32
CHAPTER 4: DATA ANALYSIS AND RESULTS 33
4.1 Study Site 33
4.2 Data Collection 36
4.2.1 Incident Data 36
4.2.2 Traffic Flow Data 39
4.2.3 Relationship between the Incident and the Traffic Flow Database 41
4.3 Data Analysis and Results 43
4.3.1. Capacity under Prevailing Conditions 44
4.3.2. Accident Capacity Reduction 47
4.3.3. Modeling Accident Capacity Reduction as a Random Variable 51
4.3.3.1 One Lane out of Three Lanes Blocked by Accidents 52
4.3.3.2 Two Lanes out of Three Lanes Blocked by Accidents 61
4.3.3.3 Shoulder Lanes out of Three Lanes Blocked by Accidents 67
4.4 Summary 68
CHAPTER 5: CONCLUSIONS 70
5.1 Conclusions 70
5.2 Contributions of This Work 72
5.3 Suggested Future Research 73
5.3.1 Improving the Quality of the Data from the Smart Travel Laboratory 73
5.3.2 Future Research Efforts 73
5.3.2.1 Improving the Methodology Used in This Research 73
5.3.2.2 Improving Traffic Condition Forecasting 74
5.3.2.3 Further Research on the Impact of Incidents on Traffic Flow 74
5.4 Summary 75
REFERENCES 76
LIST OF FIGURES
Figure 1.1: The Basic Queuing Process 3
Figure 3.1: Speed-Density Diagram 18
Figure 3.2: Speed-Flow Diagram 19
Figure 3.3: Flow-Density Diagram 19
Figure 3.4: An Example of Speed-flow Plot 21
Figure 3.5: An Example of Accident Flow 22
Figure 3.6: Moving Average of Five Successive Accident Flows 23
Figure 4.1: Hampton Roads Location Designations 34
Figure 4.2: Hampton Roads Station Designations 34
Figure 4.3: The Relationships between Incident and Traffic Flow Databases 42
Figure 4.4: Speed-Flow Diagram of Station 69 45
Figure 4.5: Speed-Flow Diagram of Station 4 45
Figure 4.6: Accident Traffic Flow of 1999-04712 48
Figure 4.7: Accident Traffic Flow of 2000-26305 49
Figure 4.8: Accident Traffic Flow of 1999-17390 at Station 131, 138 and 139 50
Figure 4.9: Histogram of Accident Capacity Reduction for One Lane out of Three Lanes Blocked 54
Figure 4.10: Density/Histogram Overplot for Accident Capacity Reduction with One Lane out of Three Lanes Blocked 58
Figure 4.11: Distribution-Function-Differences Plot for Accident Capacity Reduction with One Lane out of Three Lanes Blocked 59
Figure 4.12: Histogram of Accident Capacity Reduction for Two Lanes out of Three Lanes Blocked 63
Figure 4.13: Density/Histogram Overplot for Accident Capacity Reduction with Two Lanes out of Three Lanes Blocked 65
Figure 4.14: Distribution Function Differences Plot for Accident Capacity Reduction with Two Lanes out of Three Lanes Blocked 65
LIST OF TABLES
Table 1.1: Estimated Traffic Characteristics under the Incident Situation 4
Table 1.2: Sensitivity Analysis of Estimated Traffic Characteristics under the Incident Situation 5
Table 2.1: Percentage of Freeway Capacity Available under Incident Conditions 15
Table 4.1: Relationship between Location and Station ID 35
Table 4.2: Summary of Hampton Roads’ Incident Data 37
Table 4.3: Number of Lanes Blocked by Disabled Vehicles 38
Table 4.4: Number of Lanes Blocked by Accidents 38
Table 4.5: Duration of Disabled Vehicles and Accidents 39
Table 4.6: Estimation of Capacity under Prevailing Conditions for Selected Stations 46
Table 4.7: Estimation of capacity under Prevailing Conditions for All Stations 47
Table 4.8: Data Summary of Accident Capacity Reduction for One Lane out of Three Lanes Blocked 53
Table 4.9: The Selected Distributions and Estimated Parameters of Accident Capacity
Reduction for One Lane out of Three Lanes Blocked 56
Table 4.10: Chi-Square Test Calculations 60
Table 4.11: The Critical Values of Chi-Square Test (DF=7) 60
Table 4.12: Data Summary of Accident Capacity Reduction for Two Lanes out of Three Lanes Blocked 62
Table 4.13: The Selected Distribution and Estimated Parameters of Accident Capacity
Reduction for Two Lanes out of Three Lanes Blocked 64
Table 4.14: Chi-Square Test Calculations 66
Table 4.15: The Critical Values of Chi-Square Test (DF=5) 67
Table 4.16: Data Summary of Accident Capacity Reduction for Shoulder Lanes out of Three Lanes Blocked 68
Table 5.1: Comparison of Mean Values of This Research with Goolsby’s Results 70
Table 5.2: Comparison of the Methodology Used in This Research with the Methodology Used in Goolsby’s Study 71
GLOSSARY
Incident: any non-recurrent event that causes reduction of roadway capacity, such as accidents, debris, disabled vehicles, and hazardous material spills.
Capacity: the maximum hourly rate at which persons or vehicles can reasonably be expected to traverse a point or uniform section of a lane or roadway during a given time period under prevailing roadway, traffic, and control conditions.
Volume: the total number of vehicles that cross a point on the highway during a period of time, normally one hour
Flow: the number of vehicles passing a specific point or short section in a given period of time in a single lane. It is expressed as vehicle per hour per lane.
Speed: is defined as the average rate of motion and is expressed in miles per hour (mi/hr).
Time mean speed: the arithmetic mean of the speeds of vehicles passing a point on a highway during an interval of time.
Density: the number of vehicles occupying a section of roadway in a single lane. It is expressed as vehicle per mile per lane.
Occupancy: measure of the percentage of time for which a vehicle is detected over a detector on a highway during an interval of time.
Optimum Speed: the speed at which the flow reaches the maximum value (i.e. capacity)
Optimum Density: the density at which the flow reaches the maximum value (i.e. capacity)
Chapter 1: Introduction
Motivation
Overview
Incidents include any event that temporarily reduces roadway capacity, such as accidents, debris, disabled vehicles, and hazardous material spills (Faradyne, 2000). Incidents are a major cause of urban highway congestion. It is estimated that incidents account for 60 percent of the vehicle-hours lost to congestion (Cambridge Systematic, Inc, 1990). Incident management systems are in place in many cities to reduce the time lost due to incidents. Advanced traveler information systems are also used to provide travelers with important traffic information concerning incident congestion.
Queuing analysis and computer simulation models are important analysis techniques used in incident management systems and advanced traveler information systems. Queuing analysis can be used to estimate traffic characteristics under incident situations, including the estimation of the maximum queue length, average queue length, maximum individual delay, average individual delay, and total delay. Computer simulation models can be used to study transportation systems, transportation plans, and management strategies and evaluate the performances in the laboratory rather than in the field. Queuing analysis can be regarded as a simplified simulation model.
Remaining incident capacity is a key input of queuing analysis and computer simulation models. It is estimated that an accident or disabled vehicle blocking one out of three lanes will reduce traffic flow by an average of 50 percent, an accident blocking two out of three lanes will reduce traffic flow by an average of 79 percent, and an accident or disabled vehicle blocking shoulder lane out of three lanes will reduce traffic flow by 33 percent (Goolsby, 1970). This result is widely used by practitioners. However, very little other research has comprehensively addressed the impact of incidents on capacity. Furthermore, the stochastic characteristics of incident capacity reduction have not been investigated thoroughly. The premise of this project is that the capacity reduction caused by incidents is best modeled as a random variable, not a deterministic value, as is the current practice.
Capacity can logically be considered as a random variable rather than a deterministic value due to the variations in traffic control, weather, and other conditions. Even under ideal conditions, capacity is not a constant due to the variations in driver and vehicle characteristics. Also, incident capacity reduction is a random variable rather than a deterministic value due to the variations in incident characteristics (e.g., duration, extent, time of day, and background volume). Modeling incident capacity reduction as a random variable should provide for a more realistic estimation of incident characteristics.
Importance for Queuing Analysis
"Queuing analysis involves the mathematical study of queue that is a common
phenomenon that occurs whenever the current demand for a service exceeds the current capacity to provide that service" (Hillier and Lieberman, 1986). The basic queuing process is shown in Figure 1.1.
Queuing System
Customer Served
Served Customer
Figure 1.1: The Basic Queuing Process (Source: Frederick S. Hillier and Gerald J. Lieberman, Introduction to Operations Research, 1986)
Queuing analysis can be used to estimate traffic characteristics under the incident situation, including the estimation of the maximum queue length, average queue length, maximum individual delay, average individual delay, and total delay. The inputs include normal capacity (m ), the traffic demand (l ) when the incident occurred, the incident capacity (m R) and the duration (tR). It can be seen that incident capacity is one of the key inputs. Table 1.1 is an example of estimated traffic characteristics under the incident situation based on the deterministic queuing analysis.
Table 1.1: Estimated Traffic Characteristics under the Incident Situation
| Estimated Traffic characteristics | Equations |
| Time Duration in Queue (hours) tQ |
tR(m - m R)/ (m - l ) |
| Number of Vehicles Queued (vehicles) NQ |
l tQ |
| Maximum Queue Length (vehicles) QM |
tR(l - m R) |
| Average Queue Length (vehicles) QQ |
tR(l - m R)/2 |
| Maximum Individual Delay (minutes) dM |
60tR(l - m R)/ l |
| Average Individual Delay (minutes) dQ |
30tR(l - m R)/ l |
| Total Delay (vehicles) TD | tR tQ(l - m R)/2 |
Source: Adolf D. May, Traffic Flow Fundamentals, 1990
A sensitivity analysis study can be undertaken to assess the effect of the estimation of incident capacity reduction on the estimated traffic characteristics under the incident situation. "Sensitivity analysis involves investigating the effect on the model estimation caused by making changes in the values of the model parameters" (Hillier and Lieberman, 1986). For example, consider a three-lane directional freeway
with a total capacity of 6000 vehicles per hour. Then, assume that during the middle of the day, the traffic demand is at 80 percent of capacity. Assume an incident blocks one out of three lanes and reduces traffic flow by 46, 48, 50, 52 and 54 percent, and the reduction in capacity lasts for about 45 minutes. Table 1.2 shows the results of sensitivity analysis.
The results of sensitivity analysis indicate that the estimation of incident capacity reduction affects the estimated traffic characteristics under the incident situation significantly. Small changes in incident capacity reduction (two percent) modify the estimation of total delay significantly (11 or more percent), and thus, affect the estimation of total system performance. It is very important, therefore, to estimate incident capacity reduction accurately.
Table 1.2: Sensitivity Analysis of Estimated Traffic Characteristics Under the Incident Situation
| Estimated Traffic Characteristics |
Incident Capacity Reduction (%) | ||||
| 46% |
48% |
50% |
52% |
54% | |
|
Time Duration in Queue (hours) |
1.73 |
1.80 |
1.88 |
1.95 |
2.03 |
|
Number of Vehicles Queued (vehicles) |
8280 |
8640 |
9000 |
9360 |
9720 |
|
Maximum Queue Length (vehicles) |
1170 |
1260 |
1350 |
1440 |
1530 |
|
Average Queue Length (vehicles) |
585 |
630 |
675 |
720 |
765 |
|
Maximum Individual Delay (minutes) |
14.63 |
15.75 |
16.88 |
18.00 |
19.13 |
|
Average Individual Delay (minutes) |
7.31 |
7.88 |
8.44 |
9.00 |
9.56 |
|
Total Delay (vehicles) |
1009.13 |
1134.00 |
1265.63 |
1404.00 |
1549.13 |
"Computer simulation models incorporate queuing analysis, car-following theory, shock wave analysis, and other analytical techniques into a framework for simulating complex components or systems of interactive components" (May, 1990). CORSIM, a widely used microscopic simulation model, models the "rubbernecking" phenomenon in addition to the physical blockage on traffic flow by incidents. Drivers of vehicles in adjacent lanes tend to slow down to see what is happening as they pass an incident. This phenomenon, namely "Rubbernecking" phenomenon, results in lower speeds and then lower capacity, and explains the additional capacity reduction beyond the corresponding physical loss.
The rubberneck factor indicates the percentage of capacity reduction of each remaining lane that is not blocked during the incident. For example, if one lane out of three lanes was completely blocked and each of the remaining two lanes was influenced by a rubberneck factor of 10 percent, the capacity of the link would be reduced as shown in the following equation:
RC = (100%) (1/3) + (10%) (1/3) +(10%) (1/3)
RC = 40%
1/3 indicates the portion of volume normally carried in each lane and 10% indicates the capacity reduction in each remaining lane that is not blocked during the incident. In this example, the roadway’s capacity would be reduced by 40% during an incident.
This is a somewhat simplified model of incident capacity reduction. The calculated result, 40 percent, is less than Goolsby ‘s commonly accepted value of 50 percent capacity reduced by an accident blocking one lane out of three lanes. If analysts want to decide on an appropriate rubberneck factor, they need to collect data and perform simulations to test several alternative rubberneck factors, just as Cragg (1994) did in simulation analysis of route diversion strategies for freeway incident management. Furthermore, incident capacity reduction is modeled as a deterministic value, not a random variable in CORSIM.
As for the necessity of modeling incident capacity reduction as a random variable instead of a deterministic value in simulation models, here is an example to illustrate the danger of replacing the probability distribution of a random variable by its mean in simulation models. This example is derived from the book of simulation modeling and analysis (Law and Kelton, 2000). This example considers an incident that creates a bottleneck. It is assumed that the vehicle arrives at the bottleneck every one minute, and the vehicle departs the bottleneck every 0.99 minute. Furthermore, the interarrival times and departing times are assumed to be random variables with exponential distributions. Then based on the M/M/1 queuing analysis model, the average number of vehicles waiting in the queue can be calculated as (l /m )2/[1-(l /m )], where l is the mean interarrival rate (vehicles per time interval), and m is the mean departing rate (vehicles per time interval). The result indicates that the long-run average number of vehicles in the queue is approximately 98.
If we replace the probability distributions of the interarrival times and departing times by its mean, that is, we assume that each interarrival time is exactly one minute and each departing time is exactly 0.99 minutes, then no vehicle ever waits in the queue. The comparison of the results indicates that the variability of the probability distributions, rather than just their means, has a significant impact on the estimation of the congestion level under the incident situation.
1.2 Research Objective, Tasks and Scope
The objective of this research is to conduct a comprehensive study on the impact of incidents on capacity, based on the traffic flow and incident data for the Hampton Roads region of Virginia in the Smart Travel Laboratory of the University of Virginia. This objective is fulfilled through the following tasks:
1. Estimate capacities under prevailing conditions by calibrating speed-flow curves for several segments of Hampton Roads to provide a frame-of-reference.
2. Measure Incident capacity as the 10-minute minimum oversaturated flow in the bottleneck created by an incident.
3. Calculate the absolute value of incident capacity reduction as the difference of incident capacity and the capacity under prevailing conditions and then, calculate the percentage value of incident capacity reduction as the absolute value of incident capacity reduction over the capacity under prevailing conditions and model it as a random variable following a probability distribution.
A preliminary study indicated that few cases with measurable capacity reduction for disabled vehicles could be found based on the incident and traffic flow data for the Hampton Roads region in the Smart Travel Laboratory. It is probably because that most disabled vehicle event occurred on shoulder lanes with short durations. The statistics of durations and number of lanes blocked by different types of incidents are given in Chapter Four. The preliminary study also indicates that most cases with measurable capacity reduction are one lane, two lanes and shoulder lanes out of three lanes blocked by accidents. Based on the preliminary study, this research focuses on estimating accident capacity reduction for one lane, two lanes or shoulder lanes out of three lanes blocked and modeling it as a random variable.
One thing needs to be mentioned here is that, the percentage value of accident capacity reduction is modeled as a random variable and presented as the result in this research instead of the absolute value of accident capacity reduction. The advantage of modeling and presenting the percentage value instead of the absolute value of the accident capacity reduction is that, the accident capacity can be estimated for the other study sites given the capacity under prevailing conditions.
1.3 Organization of the Thesis
The remainder of this thesis includes a literature review, the description of methodology, data analysis and results, and the conclusions of this work.
Chapter Two describes a literature review about incident management systems, previous work on this topic, and related contents in the Highway Capacity Manual (2000).
Chapter Three first describes the existing definition of capacity and capacity estimation methods. After that, our efforts to clarify the definition of accident capacity and to develop the methodology of accident capacity estimation are discussed. Then, the methodology of modeling the accident capacity reduction as a random variable is presented.
Chapter Four describes the study site, data collection, process of calculation and modeling, and the results. This chapter discusses the methodology in more detail.
Chapter Five summarizes the conclusions and contributions of the thesis and discusses suggested further research based on this work.
1.4 Summary
This chapter describes the motivation, objective, tasks, scope and organization of this thesis. First, the importance of estimating incident capacity reduction accurately and modeling it as a random variable is discussed. Then, the objectives, main tasks, and scope of this research are described. Finally, the organization of the thesis is presented.
Chapter 2: Literature Review
2.1 Incident Management Systems
As discussed in Chapter One, incidents are a major cause of urban highway congestion. Additional fuel consumption and air pollution are commonly associated with congestion caused by incidents. The other serious problems caused by incidents include the risk of secondary crashes, and the danger of incident responders working at the scene. According to a study conducted by Minnesota Department of Transportation (1982), secondary crashes accounts for 13 percent of all crashes occurred during peak hours. Furthermore, incident responders are vulnerable and exposed to injury. According to the statistics of the IACP (1998), in 1997, nearly 40 percent of all law enforcement officers who died on duty died in traffic.
Incident management systems are now in place in many cities to minimize the congestion caused by incidents, and improve the safety of motorists, crash victims, and incident responders. For these purposes, agencies responding to incidents, such as the Police, the Fire, 911 Dispatch, the Towing and Recovery, the Emergency Medical Service (EMS), and transportation agencies, need to coordinate effectively and efficiently. In addition, the main activities of incident management systems, including detection, verification, motorist information, response, site management, traffic management and clearance, need to be performed in a systematic way. The good performance of an individual activity cannot guarantee the good performance of the whole system, which is actually the objective of an incident management system.
Some benefits of an effective incident management system can be quantified, such as improved response time, reduced delay, improved air quality, reduced occurrence of secondary incidents, and improved safety of responders, crash victims and other motorists. The other benefits are qualitative. For example, the cooperation of response agencies is improved. The public gets to know, understand, and participate in the systems. Drivers feel more comfortable knowing the operation of such systems. These qualitative benefits are also very important.
One important thing related to this research needs to be mentioned here. Most incident management systems emphasize the response and clearance activities, which are proven to be effective. Except these activities, some consideration must be given to traffic, including not only how traffic affects the incident, but also how the incident affects traffic. The knowledge of how incident affects traffic is very important for good traffic management within incident management systems, which can help fulfill the purpose of releasing congestions. For example, this research computes the accident capacity reduction and models it as a random variable. The results can be used in queuing analysis and computer simulation models to estimate the number of vehicles queued. Then management strategies, such as route diversion, can be adopted to reduce the number of vehicles queued, and therefore, release congestions to some extents.
2.2 Previous Work on This Topic
Goolsby conducted a study on the influence of incidents on freeway quality of
service in 1970. The Gulf Freeway in Houston was selected for the study because of the extensive surveillance system existing there. The 6.5-mile study section has three lanes in each direction. According to the author, an accurate log of freeway incidents including accidents and disabled vehicles was maintained on weekdays from 6:00am to 6:00pm for two years (1968 -1969) on the study section. This study collected the volume in 1-minute interval in the bottleneck created by an incident.
A total of 517 1-minute volume counts in the bottleneck were available for 27 incidents. 312 1-minute volumes under normal conditions were collected downstream of the study site to provide a frame-of-reference traffic volume. Based on these data, the author concluded that an accident or disabled vehicle blocking one lane out of three lanes reduces flow by 50 percent. An accident blocking two lanes out of three lanes reduces flow by 79 percent. An accident or disabled vehicle blocking shoulder lane(s) out of three lanes reduces flow by 33 percent.
This study initiated the research on incident capacity reduction and revealed the necessity of comprehensive research on this topic. A critical review of the paper reveals the following limitations of this study:
1. This study was conducted in 1970 and traffic conditions have changed significantly in the following 31 years. Current estimation of incident capacity reduction might be different from Goolsby’s results.
2. This study only counted 312 1-minute volumes under normal conditions downstream of the study site, not the capacity under prevailing conditions, to provide a frame-of-reference. The incident capacity reduction might be overestimated.
3. This study only used the 517 1-minute volume counts through the bottleneck created by an incident, and the measurement of incident capacity has not been well described. The incident capacity might be overestimated or underestimated because of the unstable characteristics of traffic flow rates using short measurement intervals (1-minute).
4. This study modeled the incident capacity reduction as a deterministic value, not a random variable, and the small sample size (27 incidents) made it impossible to model the incident capacity reduction as a random variable.
2.3 Related Contents in the Highway Capacity Manual (2000)
According to the Highway Capacity Manual (2000), the behavior of traffic streams during and immediately after the occurrence of an incident is not well understood. The relationships among speed, density, and flow may be discontinuous at the point of capacity and the maximum queue departing flow rate may be less than capacity under stable flow. Various observations of freeway queue departing flow rate range from 1,800 pcphpl to 2,400 pcphpl.
Estimation of percentage of freeway capacity available under incident conditions is addressed in the Highway Capacity Manual (2000). Table 2.1 illustrates the results. Unfortunately, the research work supporting this information has not been found through the literature review, including the references provided in the Highway Capacity Manual (2000). The estimated values of incident capacity reduction with one lane and two lanes out of three lanes blocked from this table are very consistent with Goolsby’s results. It is possible that Goolsby’s work in 1970 is one of the researches supporting this information.
Table 2.1: Percentage of Freeway Capacity Available under Incident Conditions
| Number of Freeway Lanes in Each Direction | Shoulder Disablement | Shoulder Accident | Lanes Blocked | ||
|
1 |
2 |
3 | |||
|
2 |
0.95 |
0.81 |
0.35 |
0.00 |
N/A |
|
3 |
0.99 |
0.83 |
0.49 |
0.17 |
0.00 |
|
4 |
0.99 |
0.85 |
0.58 |
0.25 |
0.13 |
|
5 |
0.99 |
0.87 |
0.65 |
0.40 |
0.20 |
|
6 |
0.99 |
0.89 |
0.71 |
0.50 |
0.25 |
|
7 |
0.99 |
0.91 |
0.75 |
0.67 |
0.36 |
|
8 |
0.99 |
0.93 |
0.78 |
0.63 |
0.41 |
Source: Highway Capacity Manual (2000)
2.4 Summary
This chapter gives a literature review about incident management systems, previous work on this topic, and related contents in the Highway Capacity Manual (2000). This literature review helps understand the background, motivation, and objective of this research furthermore.
Chapter 3: Methodology
3.1 Definition of Capacity
According to the Highway Capacity Manual (2000), the capacity of a facility is defined as the maximum hourly rate at which persons or vehicles can reasonably be expected to traverse a point or uniform section of a lane or roadway during a given time period under prevailing roadway, traffic, and control conditions. Also, the Highway Capacity Manual (2000) categorizes traffic flow within a basic freeway segment into three flow types: undersaturated flow, queue discharge flow, and oversaturated flow. Each flow represents different conditions on the freeway.
Undersaturated Flow represents traffic flow that is unaffected by upstream or downstream conditions. This flow is generally defined within a speed range of 55 to 75 mph at low to moderate flow rates and a range of 45 to 65 mph at high flow rates.
Queue Discharge flow represents traffic flow that has just passed through a bottleneck and is accelerating back to the free-flow speed of the freeway. This flow type is generally defined within a narrow range of flows, 2000 to 2300 pcphpl, with speeds ranging from 35 mph up to the free-flow speed of the freeway section.
Oversaturated Flow represents traffic flow that is influenced by the effects of a downstream bottleneck. Traffic flow in the congested regime can vary over a broad range of flows and speeds depending on the severity of the bottleneck.
According to these definitions and categorizations, the accident capacity can be defined as the minimum 15-minute oversaturated flow at the upstream of a bottleneck created by an accident. Because oversaturated flow upstream of a bottleneck means that the capacity level has been reached in the bottleneck, it is possible to make a reliable capacity estimate. Also, given the variation in observed capacities, analysts may wish to use appropriate measurement interval to get the minimum 15-minute oversaturated flow as accident capacity to reduce the risk of capacity overestimates or underestimates.
According to the Highway Capacity Manual (2000), the time period used in most capacity analyses is 15-minute, which is considered to be the interval during which stable flow exists. In order to provide a quantitative understanding of the impacts of short (less than 15-minute) measurement intervals on traffic flow rate, Smith (2001) conducted a study and found that stable flow rates may be calculated using measurement intervals as short as 10-minute, and that statistically significant improvements in stability can be achieved by adding 2-minute to any measurement interval. According to the Highway Capacity Manual (2000) and Smith’s results (2001), the 10-minute interval is used in this work to estimate accident capacity reduction. 10-minute intervals are also used in this work to estimate capacity under prevailing conditions to provide a frame-of-reference.
3.2 Capacity Estimation Methods
3.2.1 Estimation of Capacity under Prevailing Conditions
Capacity under prevailing conditions can be estimated by calibrating a speed-
flow and /or density-flow curve for a given segment of highway. The peak of this curve defines capacity. This kind of methods is based on the fundamental models describing speed-flow, flow-density, and speed-density relationships. Theoretical speed-density, speed-flow, and flow –density diagrams are shown in Figure 3.1, Figure 3.2 and Figure 3.3. A linear speed-density relationship is assumed to simplify these models. Flow is defined as the number of vehicles passing a specific point or short section in a given period of time in a single lane. It is expressed in vehicles per hour per lane. Speed is defined as the average rate of motion and is expressed in miles per hour. Density is defined as the number of vehicles occupying a linear section of roadway in a single lane. It is expressed in vehicles per mile per lane. qm is capacity under prevailing conditions that could be estimated empirically or calculated theoretically. U0 is optimum speed corresponding to qm. k0 is optimum density corresponding to qm.

Figure 3.1: Speed-Density Diagram

Figure 3.2: Speed-Flow Diagram

Figure 3.3: Flow-Density Diagram
Speed-density, speed-flow, and flow–density diagrams shown in Figure 3.1, Figure 3.2 and Figure 3.3 can be expressed as theoretical models in mathematical forms. Efforts have been devoted to calculate the capacity value from the theoretical models. But the stochastic nature of the observations near capacity makes it very difficult. Instead of this kind of analytical approach, the empirical approach can be used to estimate the capacity value.
Venkatanarayayana (2001) developed the speed-flow plots and estimated the capacities under prevailing conditions for several segments of Hampton Roads through an empirical approach in the Smart Travel Laboratory during the summer of 2000. According to Venkatanarayayana (2001), on plotting speed- flow diagram, a near parabola bounded by a top line can be observed for most of the segments of Hampton Roads considered. Figure 3.4 shows the speed-flow plot with 142,000 10-minute observations from June 15th, 1998 to July 29th, 1999 on the eastbound I-64. The top line exists there because of the maximum speed limits on the roads (65 mph). Also, a number of points can be seen inside the parabola because different roadway, traffic, and control conditions exist even on the uniform segments of the Hampton Roads. The process that Venkatanarayayana (2001) used to estimate the capacity for several segments of Hampton Roads during a given time period is discussed in Chapter Four in detail.
The traffic flow rate can be taken to the nearest integer, in multiples of 50, because it is impossible and unnecessary to measure traffic flow rate exactly and errors exist inherently. In Figure 3.4, a lot of points can be seen on the "peak" of the curve and form a "peak area". The traffic flow rates within this "peak area" vary from 2450 vphpl to 2550 vphpl. To reduce the risk of underestimates or overestimates, the integer, in multiple of 50, within the "peak area", 2500 vphpl is regarded as the capacity under prevailing conditions from June 15th, 1998 to July 29th, 1999 on the eastbound of I-64.

Figure 3.4: An Example of Speed-flow Plot
3.2.2 Accident Capacity Estimation
When an accident occurs and a bottleneck is formed, the capacity of the roadway
is reached. There is no space for more traffic flow and traffic demand exceeds the capacity at that time. Thus accident capacity can be measured directly as bottleneck traffic flow according to the definition of capacity as the maximum hourly rate at which persons or vehicles can reasonably be expected to traverse a point or uniform section of a lane or roadway during a given time period under prevailing roadway, traffic, and control conditions (Highway Capacity Manual, 2000).
Figure 3.5 provides an example of accident flow. This accident occurred on I-564 and blocked the first lane. The duration period was recorded in the incident database from 6:29am to 6:44am. Traffic volume was collected every two minutes and transformed to an hourly rate. Figure 3.5 plots traffic flow from 30-minute before the accident to 30-minute after the accident. It shows that traffic flow fell at 6:18am, and was restored at 6:48am. When this kind of pattern appears, the judgment can be made that a bottleneck is formed and the capacity of the roadway is reached. For this accident situation, traffic flow began to fall before the beginning time recorded. It is reasonable because some accidents might not be detected and recorded immediately after it occurred.

Figure 3.5: An Example of Accident Flow
The variation of traffic flow can be seen clearly in this figure because the time interval in this practice is 2-minute. There exists the risk of overestimates or underestimates of accident capacity because of the variation of traffic flow. To reduce the effect of the variation of traffic flow on accident capacity estimation, the moving average of five successive flows is calculated and then the minimum value is selected as an accident capacity. As discussed before, 10-minute interval is used in this work to estimate accident capacity reduction. This is why the moving average of five successive flows is used to estimate the accident capacity. For example, the moving average of five successive flows can be calculated as:
Mt = (Xt-2 + Xt-1 +Xt + Xt+1 + Xt+2) / 5
Where Mt is the moving average of five successive flows
Xt is the traffic flow at time t
Xt+1 is traffic flow at about two minutes after time t
Xt-1 is traffic flow at about two minutes before time t
The moving average of five successive accident flows of the previous example is shown in Figure 3.6. The minimum value (1953 veh/hour) can be regarded as the accident capacity for this example. To simplify this task, analysts could base on his/her own judgment to select the five successive flows and calculate the average as the accident capacity. The criteria is that the average of these five successive flows is the minimum of all the average of these five successive flows during the period from 30-minute before the accident to 30-minute after the accident. For the illustrated example, accident capacity could be calculated as the average of traffic flow from 6:26am to 6:36am, which is 1953 veh/hour.
![]()
Figure 3.6: Moving Average of Five Successive Accident Flows
3.3 Modeling Accident Capacity Reduction as a Random Variable
3.3.1. Overview
The absolute value of accident capacity reduction can be calculated as the difference of the estimation of the capacity under prevailing conditions and accident capacity. Then the percentage value of accident capacity reduction can be calculated as the accident capacity reduction over the capacity under prevailing conditions. The following work is to model accident capacity reduction (refer to the percentage value) as a random variable, not a deterministic value. Chapter One has explained the importance of modeling accident capacity reduction as a random variable. This section will describe how to choose the distribution that best represent the observed data set. Because accident capacity reduction is continuous, this section will only address selecting the continuous distribution that best represent the observed data set.
The activities involve in choosing the "correct" distribution include 1. Families of distributions that might be representative of the observed data set are hypothesized on the basis of the summary statistics, histogram of the observed data set, the shapes of the hypothesized distributions, and the other information, 2. The parameter values are specified for these candidate distributions, and 3. The hypothesized distributions are evaluated on how representative they are for the observed data set and the best- fitted distribution is selected.
3.3.2. Activity Ⅰ: Hypothesizing Families of Distributions
3.3.2.1 Summary Statistics
Summary statistics, including minimum observation, maximum observation, mean, median, variance, coefficient of variation, and skewness, can be used to suggest families of distributions that might fit the observed data set. The coefficient of variation, calculated as the standard deviation divided by the mean, can be used to measure the variability or dispersion of a distribution as an alternative to variance. The skewness is a measure of the symmetry of the distribution. The skewness of a symmetric distribution is zero. The skewness of a distribution that has a longer "right tail" than "left tail" is positive. The skewness of a distribution that has a longer "left tail" than "right tail" is negative.
By comparing the summary statistics of the observed data set and characteristics of the theoretical distributions, families of distributions may be hypothesized. For example, the coefficient of variation of the exponential distribution is equal to one. The skewness of the normal distribution is equal to zero. The comparison between the summary statistics of the observed data set and characteristics of the theoretical distributions can at least exclude some theoretical distributions from the hypothesized family of distributions. For example, the normal distribution cannot represent the observed data set with the skewness different from zero apparently. Also, the exponential distribution cannot represent the observed data set with the coefficient of variation different from one apparently.
3.3.2.2 Histogram
A histogram provides a graphical estimation of the plot of the density function corresponding to the fitted distribution of the observed continuous data set. To make a histogram, the range of the values in the data set is broken into k adjacent intervals [b0, b1) , [b1, b2) , … , [bk-1, bk]. These intervals have the same interval width D b = bj – bj-1. The definition of the histogram function h(x) is as follows:
0 if x < b0
h(x) = hj if bj-1 £ x < bj for j = 1,2, …k
0 if x£ b0
where hj is the proportion of the Xi’s that are in the jth interval [bj-1, bj].
The sample size and the interval width affect whether the histogram is representative of the observed data set. If the sample size is too small, the histogram will often be "ragged" regardless of how the interval width is chosen. If the interval width is too small, the histogram will also be "ragged". If the interval width is too big, it is possible that too much information in the observed data set is left out and the histogram cannot be representative. There is no definitive guide for choosing the interval width. Neiswanger (1943) proposed that the number of intervals should be between 10 and 25. Sturges (1926) proposed that the number of intervals should be estimated as follows:
Range
I =
1 + (3.322) logN
where I = number of intervals
Range = largest observed value minus smallest observed value
N = number of observations
Law and Kelton (2000) regarded such rules as not very useful and recommended trying several different values of interval width and choosing the smallest one that gives a "smooth" histogram. These rules and recommendations will be considered comprehensively in this research.
3.3.3. Activity Ⅱ: Estimating Parameters
In activity Ⅰ,
families of distributions have been hypothesized to represent the observed data
set. The following work is to estimate the parameter values to obtain the completely
specified distributions from the observed data set. Maximum-likelihood estimator
(MLEs) is a widely used technique for estimating parameter values for a given
distribution. In the discrete case, given the observed data set
X1, X2, …, Xn, suppose a discrete distribution
has been hypothesized for this data set with one unknown parameter q . Let
denote
the probability mass function for this distribution. The likelihood function is
defined as
L(q )
=
... ![]()
Which is the joint probability mass function when the data are independent. The
parameter q is then estimated as the value q that maximizes L(q ) over all other permissible values of q .
In the
continuous case, density function for the continuous distribution,
instead of probability mass function for the discrete distribution, is used here
because the probability for a continuous random variable equal to any fixed number
is zero. Given the observed data set X1, X2, …, Xn,
suppose a continuous distribution has been hypothesized for this data set with
one unknown parameter q . Let
denote
the density function for this distribution. The likelihood function is defined
as
L(q ) =
…![]()
The parameter q is then estimated as the value q that maximizes L(q ) over all other permissible values of q . If the hypothesized distribution has several unknown parameters, these parameters can be estimated as the values that jointly maximizes L(q ) over all other permissible values of these parameters.
3.3.4 Activity Ⅲ: Determining How Representative the Fitted Distributions Are
In activity Ⅱ, one or more completely specified distributions have been
determined with the estimated parameter values. The following work is to examine how well these distributions represent the true underlying distribution of the observed data set. Furthermore, the distribution that provides the best-fit needs to be determined if several distributions are evaluated as "representative". Heuristic procedures and goodness-of-fit hypothesis tests can be used to determine how representative the fitted distributions are.
3.3.4.1 Heuristic Procedures
The Density/Histogram overplot and the Distribution-Function-Differences Plot are two widely used heuristic procedures for comparing fitted distributions with the true underlying distribution. For continuous observed data set, the Density/Histogram overplot can be made by plotting the density function of the fitted distribution over the histogram. If the fitted distribution is a good representation for the observed data set and the sample size is sufficiently large, the density function of the fitted distribution and the histogram should look very similar.
The Density/Histogram overplot compares the individual probabilities of the fitted distribution and the true underlying distribution, while the Distribution-Function-Differences plot compares the cumulative probability of the fitted distribution and the true underlying distribution. The empirical cumulative distribution for the observed data set is defined as:
Number of Xi’s £ X
Fn(x) =
N
The Distribution-Function-Differences plot is a plot of the differences between the
fitted
cumulative distribution
and
the empirical cumulative distribution Fn(x) for the observed data set.
If the fitted distribution is a good representation for the observed data set
and the sample size is sufficiently large, the Distribution-Function-Differences
plot should be close to a horizontal line at height zero.
3.3.4.2 Goodness-of-fit Test
While
heuristic procedures can provide graphic comparison between the fitted distribution
and the true underlying distribution, goodness-of-fit test can assess formally
whether the observed data set is an independent sample from a particular fitted
distribution (Law and Kelton, 2000). In other words, given the observed data set
X1, X2, …, Xn, and the density function of the
fitted distribution
,this statistical hypothesis test can be used to test the following null hypothesis:
H0:
The observed data set X1, X2, …, Xn is an independent
sample from the distribution with density function ![]()
According to Law and Kelton (2000), the t-statistics is computed as the measure of the "distance" between the fitted distribution and the observed data set. Different goodness-of-fit tests, such as the Chi-Square test, the Kolmogorov-Smirnov test, and the Anderson-Darling test, use different methods to calculate the T value. Decision is made whether to reject the null hypothesis based on the comparison of T and a critical value t (a ). The probability that T is greater than t (a ) is equal to a when the null hypothesis is true. Typically, the level a is set to be 0.05 to 0.1. The critical value t (a ) depends on a , the sample size n, the test itself, and sometimes, the fitted distribution. For a test of level a , if T is greater than t(a ), reject the null hypothesis. Otherwise, do not reject the null hypothesis.
Failure to reject the null hypothesis should not be interpreted as "accepting the null hypothesis as being true" (Law and Kelton, 2000). That is to say, goodness-of-fit tests can only assist determining whether to reject the null hypothesis. Whether to accept the null hypothesis cannot be determined through these tests. If the sample size is small, it is difficult for goodness-of-fit tests to examine the disagreements between the true underlying distribution of the observed data set and the fitted distribution based on the limited information provided by the short data set.
This research
uses the chi-square test to assess whether the observed data set is an independent
sample from a particular fitted distribution. To compute the Chi-Square test statistic
in the continuous case, the entire range of the fitted distribution is first broken
into k adjacent intervals [b0, b1) , [b1, b2)
, … , [bk-1, bk]. Then given the density function
of the fitted distribution,
the proportion pj of the Xi’s that would fall in the jth
interval is computed as:

Finally, the chi-square test statistic is computed as:
![]()
The
c 2 test statistic for large n is approximately a c 2 distribution.
H0 is rejected if
.
The number of degrees of freedom of the c 2 distribution is calculated
as:
N=(I - 1) – P
where N = number of degrees of freedom
I = number of intervals being compared
P = number of parameters estimated
3.4 Summary
This chapter describes the methodology used in this research. First, the existing definition of capacity and capacity estimation methods is presented. After that, the definition of accident capacity and methodology of estimating capacity under prevailing conditions and estimating accident capacity used in this research are discussed. Finally, the methodology of modeling the accident capacity reduction as a random variable is described.
Chapter 4: Data Analysis and Results
4.1 Study Site
The Hampton Roads region of Virginia is selected for this research because the Smart Travel Laboratory of the University of Virginia has incident data dating from March 28th, 1992 and traffic flow data dating from June 28th, 1998 for this region. The Smart Travel Laboratory is connected to the Hampton Roads Smart Traffic Center that is the main freeway data resource for this laboratory. For the purpose of traffic management, the Hampton Roads are divided into several segments, which are named as "locations" in the incident database. The locations with incident data include: 564-01, E64-01, E64-02, E-64-03, W64-01, W64-02, W64-03, W64-04, W64-05, W64-06, W64-07, and w64-08. The length of each segment is about one to three miles. Accident data was recorded by locations in the incident database. The location designations can be seen in Figure 4.1.
Traffic flow data was recorded by inductive loop detectors installed at 203 places on the interstate facilities of the region, referred to as "stations". In general, there are several stations within each location in each direction, and each station consists of a single loop detector installed in each of the travel lanes. Traffic flow data was recorded by station ID (the identification number of the station) in the traffic flow database. The relationship between location and station ID is the key to relate traffic flow data to incident data. The station designations can be seen in Figure 4.2, and the relationship between location and station ID can be seen in Table 4.1.

Figure 4.1: Hampton Roads Location Designations

Figure 4.2: Hampton Roads Station Designations
Table 4.1: Relationship between Location and Station ID
| Location |
Station ID |
| 564-01 |
E 132,135,136 W 131,138,139 |
| E64-01 |
E 3,8 W 4,6 |
| E64-02 |
E 22,24,26,68 W 21,23,30,44 |
| E64-03 |
E 31 W 15 |
| W64-01 |
E 19,39,43,47,51 W 17,36,40,46,54 |
| W64-02 |
E 58,60 W 56,62 |
| W64-03 |
E 67, 71 W 65, 69 |
| W64-04 |
E 76, 83 W 80,81,85 |
| W64-05 |
E 87,92 W 91,96 |
| W64-06 |
E 98,105,111 W 104,108,117 |
| W64-07 |
E 120,123 W 122 |
| W64-08 |
E 126 W 125 |
This table only lists the stations on mainlines and does not include the stations on ramps. This research only uses the traffic flow data corresponding to the stations on mainlines. Traffic flow on ramps is often very low. When an accident occurred on ramps, the bottleneck often cannot be established and the methodology of accident capacity reduction discussed in Chapter Three cannot be used.
Station ID is corresponding to location and direction. For example, station 132 is on the eastbound of the location 564-01, while station 131 is on the westbound of the location 564-01. There is no information of the exact place where the accident occurred within the location. There are several stations within one location. It is possible that the traffic flow recorded by one station was affected by the accident more significantly than the traffic flow recorded by the other stations. It is difficult to decide which station best provides the information of traffic flow affected by accident to estimate accident capacity. This difficult situation needs to be considered during the process of accident capacity estimation and will be discussed later on.
4.2 Data Collection
4.2.1 Incident Data
The Hampton Roads’ incident data in the Smart Travel Laboratory describes the type, location, lanes blocked, beginning time, duration and the other information of each reported incident that occurred at the locations described in section 4.1 dating from March 28th, 1992. The type of incidents, number of each type of incidents, and percentages of each type of incident dating from March 28th, 1992 until February 18th, 2001 (the end date for data analysis in this research) are summarized in Table 4.2. It can be seen that disabled vehicles account for 72.9 percent of reported incidents and accidents account for 8.2 percent of reported incidents. The other types, such as Abandoned, Bridge, Debris, TEOC, Tunnel, and VMS Change, are special cases reported for the Hampton Roads Smart Traffic Center. They are not within the scope of this research.
Table 4.2: Summary of Hampton Roads’ Incident Data in the Smart Travel Laboratory
| Type of incident | Total number | Percentage |
| Disabled |
68410 |
72.9% |
|
Accident |
7698 |
8.2% |
|
Abandoned |
6958 |
7.4% |
|
Bridge |
2945 |
3.1% |
|
Debris |
2895 |
3.1% |
|
TEOC |
1843 |
2.0% |
|
Tunnel |
1734 |
1.8% |
|
Condition Chg |
752 |
0.8% |
|
VMS Change |
254 |
0.3% |
|
Other |
404 |
0.4% |
|
Total |
93893 |
100% |
The Number of Lanes blocked by disabled vehicles and accidents, number of incidents of each category, and corresponding percentage are summarized in Table 4.3 and Table 4.4. The total number of lanes varies from two to four. There are more than one hundred descriptions of which lane or lanes were blocked by the incidents in the incident database. Those that cannot be categorized into one, two, or shoulder lanes blocked are categorized into "other".
Table 4.3: Number of Lanes Blocked by Disabled Vehicles
| Number of lanes blocked | Total Number | Percentage |
| One |
1513 |
0.22% |
|
Two |
54 |
0.10% |
|
Shoulder |
64518 |
94% |
|
Other |
2323 |
3.4% |
|
Total |
68408 |
100% |
Table 4.4: Number of Lanes Blocked by Accidents
| Number of lanes blocked | Total Number | Percentage |
| One |
2063 |
26.8% |
|
Two |
1593 |
20.7% |
|
Shoulder |
2827 |
36.7% |
|
Other |
1275 |
16.6% |
|
Total |
7698 |
100% |
The durations of incidents, number of incidents within each category of duration, and corresponding percentages of incidents dating from March 28th, 1992 until February 18th, 2001 are summarized in Table 4.5. There are three categories for the durations of incidents in this table: short, 0-15 minutes; medium, 15-30 minutes; and long, 30+ minutes. This categorization is widely used in practice. Also, only disabled vehicles and accidents are considered because the other types of incidents are not within the scope of this research, as mentioned before.
Table4.5: Duration of Disabled Vehicles and Accidents
| Duration |
0-15 minutes |
15-30 minutes |
30+ minutes |
Total Number | |||
| Number |
Percentage |
Number |
Percentage |
Number |
Percentage | ||
| Disabled |
48210 |
70.5% |
12556 |
18.4% |
7644 |
11.2% |
68410 |
|
Accident |
1811 |
23.5% |
1549 |
20.1% |
4338 |
56.4% |
7698 |
It can be seen that most disabled vehicles occurred on shoulder lanes with short durations, while most of the accidents blocked one lane, two lanes and shoulder lanes with long durations. This may be the reason that few cases with significant capacity reduction for disabled vehicles can be found, and most cases with significant capacity reduction are one lane, two lanes and shoulder lanes out of three lanes blocked by accidents, as discussed in Chapter One. The scope of this research is thus limited to estimating accident capacity reduction for one lane, two lanes or shoulder lanes out of three lanes blocked and modeling it as a random variable.
4.2.2 Traffic Flow Data
The Hampton Roads’ traffic flow data in the Smart Travel Laboratory records date, station ID, traffic volume, traffic occupancy, time mean speed, collect length and the other information corresponding to the station ID from June 28th, 1998 and is updated every two minutes for the Hampton Roads region. The traffic volume is the number of vehicles on all the lanes of the freeway counted during the collection length. The collection lengths are between 100 and 255 seconds with the average of slightly more than 120 seconds. Traffic flow rate in vphpl is calculated by the formula below.
Traffic flow rate (vphpl) =
(volume /collect length in sec /number of lanes) * 3600 sec/hour
Thus, through this calculation the traffic volume in vehicle per second is transformed to traffic flow rate in vehicle per hour per lane. The stable traffic flow is calculated as 10-minute average traffic flow, as discussed in Chapter Three.
One important thing needs to be mentioned here. There is no information about the types of vehicles in the Hampton Roads’ traffic flow data. As we know, the heavy vehicles other than passenger cars, such as trucks, buses, and recreational vehicles, affects traffic flow, and thus capacity of the roadway. The Highway Capacity Manual (2000) presents tables to determine the numerical value of the heavy vehicle factor that is used to consider the reduction in capacity due to the presence of heavy vehicles.
Subramanyan (2000) presented the results of the estimation of the truck percentage at four study sites on I-64 and I-44 based on the field data in his thesis. The truck percentage on the westbound of I-64 was estimated as 10 percent, while the truck percentage on the eastbound of I-64 was estimated as 8 percent. The other two study sites are not related to this research. This research does not consider the heavy vehicle factor because the estimation of capacity under prevailing conditions and accident capacity reduction will all decrease due to the heavy vehicle factor. The percentage value of accident capacity reduction will not change so much when considering the heavy vehicle factor.
Time mean speed is defined as the arithmetic mean of the speeds of vehicles passing a point on a highway during an interval of time. Individual speeds are recorded passing a point, and are arithmetically averaged. Another speed measure, space mean speed, is defined as the length of the segment divided by the average running time of vehicles to traverse the segment. This research uses time mean speed as the speed measure. Occupancy is defined as the measure of the percentage of time for which a vehicle is detected over a detector on a highway during an interval of time. This research is not concerned about the occupancy.
4.2.3 Relationship between the Incident and the Traffic Flow Database
Four databases stored in the Smart Travel Laboratory are used in this research, including (1) HR.Inc_Roadway, which describes the location where an incident occurred; (2) HR.Incident, which stores information about each incident, uniquely identified by the identification number of incidents TMS_CALL_NUMBER; (3) HR.Station_Data, which records that traffic flow data corresponding to the station ID; This table stores information about each station; and (4) HR.Station_Info, which describes the relationship between location and station ID. One thing needs to be mentioned here to eliminate confusions. In some cases one keyword has different field names in the different databases. For example, the field "Location" in (1) HR.Inc_Roadway and the field "Interstate" in (4) HR.Station_Info all refer to the segments of Hampton Roads, while the field "Direction" in (1) HR.Inc_Roadway and the field "Lane_Type" in (4) HR.Station_Info all refer to roadway direction. This needs to be taken care of when connecting databases by keywords. Please refer to database documentation in the Smart Travel Laboratory to get more detailed information about these databases. The relationships among these four databases are shown in Figure 4.3.
Figure 4.3: The Relationships between Incident and Traffic Flow Databases
The traffic flow from 30 minutes before the accident occurred to 30 minutes after the accident cleared at the special station corresponding to the location where the accident occurred can be selected using the following steps:
Step1: Select incidents occurred on the specified location by the fields "Location" and "Direction" from (1) HR.Inc_Roadway. For example, all incidents occurred on the eastbound ("Direction") of 564-01("Location") can be selected from (1) HR.Inc_Roadway.
Step2: Select accidents from the incidents selected at Step One by the fields "Type" in (2) HR.Incident. This can be realized because (1) HR.Inc_Roadway and (2) HR.Incident are related through TMS_CALL_NUMBER. Then decide the beginning time and the ending time for each selected accidents by the field "Begin" and "Duration" in (2) HR.Incident. The ending time is equal to the beginning time plus the duration in minute.
Step3: Find the stations corresponding to the specified location from (4) HR.Station_Info. For example, Station 132, 135, and 136 are on the eastbound of the location 564-01 (Please refer to table 4.1). Then the traffic flow from 30 minutes before the accident occurred to 30 minutes after the accident cleared at the special station corresponding to the location where the accident occurred can be selected from (3) HR.Station_Data by the joint information of the station ID, beginning time, and ending time decided in the previous steps.
4.3 Data Analysis and Results
4.3.1. Capacity under Prevailing Conditions
Using the method of estimating of capacity under prevailing conditions discussed in Chapter Three, capacities for several segments of Hampton Roads under prevailing conditions were estimated using the following steps:
Step1: Several segments of Hampton Roads were selected to estimate capacities under prevailing conditions. These segments were selected because each of them is believed to have uniform characteristics. These segments include I564 West Bound, I64 (south of I564 jct) East Bound, I64 (south of I564 jct) West Bound, and I64 (south of I264 jct) West Bound (Please refer to figure 4.1). Then station 139, 39, 69, and 4 were selected to estimate capacities for the corresponding segments under prevailing conditions.
Step2: Five successive traffic flow values were summarized to obtain volumes over 10-minute. The volumes over 10-minute were then converted to flow rates in vehicle per hour per lane. The time mean speeds were averaged within every 10-minute interval. Speed-flow diagram was plotted using about one year’s traffic flow data for each station. On plotting speed-flow diagram, a near parabola bounded by a top line was observed for most of the basic sections considered. The speed-flow plots of station 69 and 123 are shown as expected in Figure 4.4 and Figure 4.5.
Step3: The capacities of station 139, 39, 69, and 4 were first estimated. The results are shown in Table 4.6.

Figure 4.4: Speed-Flow Plot of Station 69

Figure 4.5: Speed-Flow Plot of Station 123
Table 4.6: Estimation of Capacity under Prevailing Conditions for Selected Stations
| Station ID | Location Description |
No. of lanes |
Max flow (vphpl) |
Data from (date) |
Data to (date) | Total data points (‘000) |
| 139 |
I564 West Bound |
3 |
1950 |
8/8/99 |
8/16/00 |
101 |
|
39 |
I64 (south of I564 jct) East Bound |
3 |
2500 |
6/15/98 |
7/29/99 |
142 |
|
69 |
I64 (south of I564 jct) West Bound |
3 |
2450 |
8/2/99 |
7/31/00 |
110 |
|
4 |
I64 (south of I264 jct) West Bound |
3 |
2250 |
8/2/99 |
7/31/00 |
140 |
Source: Venkatanarayayana’s Working Paper at the Smart Travel Laboratory
Step4: The estimation of capacity under prevailing conditions of station 139, 39, 69, and 4 can be used as the estimation of capacity for the corresponding locations under prevailing conditions, and then the estimation of capacity related to each station can be decided by the relationship between locations and stations. The results are shown in Table 4.7. The capacity under prevailing conditions is related to each station in Table 4.1 to provide a frame-of-comparison because estimation of accident capacity reduction is related to one of the stations, not limited to station 139, 39, 69, and 4.
Table 4.7: Estimation of Capacity under Prevailing Conditions for All Stations
| Location |
Station id |
No. of lanes |
Maximum flow (veh/hour/lane) |
| 564-01 |
W 131,138,139 |
W 3 |
W 1950 |
|
E64-01 |
W 4,6 |
W 3 |
W 2250 |
|
E64-02 |
W 21,23,30,44 |
W 3 |
W 2250 |
|
E64-03 |
W 15 |
W 3 |
W 2250 |
|
W64-01 |
E 19,39,43,47,51 W 17,36,40,46,54 | E 3 W 3 |
E 2500 W 2400 |
| W64-02 |
E 58,60 W 56,62 | E 3 W 3 |
E 2500 W 2400 |
| W64-03 |
E 67, 71 W 65, 69 | E 3 W 3 |
E 2500 W 2450 |
| W64-04 |
E 76, 83 W 80,81,85 | E 3 W 3 |
E 2500 W 2400 |
| W64-05 |
E 87,92 W 91,96 | E 3 W 3 |