Free access
TECHNICAL PAPERS
Jul 15, 2009

Highway Project Level Life-Cycle Benefit/Cost Analysis under Certainty, Risk, and Uncertainty: Methodology with Case Study

Publication: Journal of Transportation Engineering
Volume 135, Issue 8

Abstract

One of the key steps in the highway investment decision-making process is to conduct project evaluation. The existing project level life-cycle cost analysis approaches for estimating project benefits maintain limited capacity of probabilistic risk assessments of input factors such as highway agency costs, traffic growth rates, and discount rates. However, they do not explicitly address cases where those factors are under uncertainty with no definable probability distributions. This paper introduces an uncertainty-based methodology for highway project level life-cycle benefit/cost analysis that handles certainty, risk, and uncertainty inherited with input factors for the computation. A case study is conducted to assess impacts of risk and uncertainty considerations on estimating project benefits and on network-level project selection. First, data on system preservation and expansion, usage, and candidate projects for state highway programming are used to compute project benefits using deterministic, risk-based, and uncertainty-based analysis approaches, respectively. Then, the three sets of estimated project benefits are implemented in a stochastic optimization model for project selection. Significant differences are revealed with and without uncertainty considerations.

Introduction

Risk and Uncertainty Considerations in Highway Project Level Life-Cycle Benefit/Cost Analysis

One of the key steps in the highway investment decision-making process is to realistically estimate project level life-cycle costs and benefits of different types of highway projects. Different highway facilities such as pavements and bridges have different useful service lives. In order to compare the merit of different projects on an equal basis, the life-cycle cost analysis approach needs to be adopted to evaluate the total economic worth of the initial construction cost and discounted future maintenance and rehabilitation costs in the facility life cycle.
Over the last decade, the life-cycle cost analysis approach has been widely used for evaluating highway pavement and bridge projects. As related to pavement project evaluation, the Federal Highway Administration (FHwA) made a concerted effort for the use of life-cycle cost analysis in highway pavement design (FHwA 1998). Hicks and Epps (1999) explored alternative pavement life-cycle design strategies with a logical comparison between conventional mixtures and the mixture containing asphalt rubber pavement materials. Wilde et al. (1999) introduced a life-cycle cost analysis framework for rigid pavement design. Abaza (2002) developed an optimal life-cycle cost analysis model for flexible pavements. Falls and Tighe (2003) enhanced life-cycle cost analysis through the development of cost models using the Alberta roadway maintenance and rehabilitation analysis application. Labi and Sinha (2005) and Peshkin et al. (2005) studied systematic preventive maintenance and the optimum timing strategies to achieve minimum pavement life-cycle costs. Chan et al. (2008) evaluated life-cycle cost analysis practices in Michigan. For bridge project evaluation, Purvis et al. (1994) performed life-cycle cost analysis of bridge deck protection and rehabilitation. Mohammadi et al. (1995) introduced the concept of incorporating life-cycle costs into highway bridge planning and design. Hawk (2003) developed a bridge life-cycle cost analysis software tool for bridge project evaluation.
In recent years, researchers began to utilize the risk-based life-cycle cost analysis approach to establish mathematical expectations of highway project benefits. For instance, Tighe (2001) performed a probabilistic life-cycle cost analysis of pavement projects by incorporating mean, variance, and probability distribution for typical construction variables, such as pavement structural thickness and costs. Reigle et al. (2005) incorporated risk considerations into the pavement life-cycle cost analysis model. Setunge et al. (2005) developed a methodology for risk-based life-cycle cost analysis of alternative rehabilitation treatments for highway bridges using Monte Carlo simulation.

Motivation and Study Objectives

Because of lacking of pertinent information, in many instances it might not be possible to establish a meaningful probability distribution to possible outcomes of a specific input factor such as construction, rehabilitation, and maintenance costs and traffic growth in risk-based life-cycle benefit/cost analysis. That is, the input factors are under uncertainty with no definable probability distributions. Consequently, the mathematical expectation of the input factor cannot be established. Further, risk and uncertainty inherited with input factors for project level life-cycle benefit/cost analysis may vary from project to project. Some projects may only involve risk cases for some input factors, whereas other projects may only experience uncertainty cases for some input factors. In more general situations, a project may face mixed cases of certainty, risk, and uncertainty concerning all input factors for the computation. This necessitates developing a new uncertainty-based methodology for highway project level life-cycle benefit/cost analysis that could rigorously handle such general situations.
In this paper, we will first introduce a methodology for highway project level life-cycle benefit/cost analysis that considers certainty, risk, and uncertainty associated with input factors for the computation. A case study is then conducted to assess impacts of risk and uncertainty considerations in estimating project level life-cycle benefits and on network-level project selection. Discussions of usefulness of the proposed methodology and directions to its refinement are provided in the last section.

Framework of the Proposed Methodology

The section starts with the discussion of common agency cost and user cost categories for pavement and bridge facilities, respectively. It then introduces a project level life-cycle cost analysis approach for computing agency costs and user costs, as well as estimating overall project level life-cycle benefits for pavements and bridges. Next, risk and uncertainty issues associated with input factors for the computation are addressed. The last part of this section provides a generalized framework for uncertainty-based highway project level life-cycle benefit/cost analysis where the input factors are under certainty, risk, and uncertainty.

Pavement and Bridge Life-Cycle Agency Costs and User Costs

In this study, the pavement or bridge life-cycle is defined as the time interval between two consecutive construction events. Maintenance and rehabilitation treatments are performed within the pavement or bridge life-cycle. The pavement and bridge life-cycle agency cost and user cost components are briefly discussed in the following:

Pavement Life-Cycle Agency Costs

Cost analysis is a cardinal element of any highway project life-cycle benefit/cost analysis. All costs incurred over pavement life-cycle including those of construction, rehabilitation, and maintenance treatments need to be included into the analysis.

Bridge Life-Cycle Agency Costs

Bridge agency costs are primarily involved with costs of bridge design and construction/ replacement, deck and superstructure rehabilitation and replacement, and maintenance treatments.

Pavement/Bridge Life-Cycle User Costs

User costs are incurred by highway users in the pavement or bridge life cycle. User cost components mainly include costs of vehicle operation, travel time, vehicle crashes, and vehicle air emissions (FHwA 2000; AASHTO 2003). Each user cost component consists of two cost categories: user cost under normal operation conditions and excessive user cost due to work zones (FHwA 1998).

Pavement/Bridge Life-Cycle Activity Profiles and User Cost Profiles

Pavement/Bridge Life-Cycle Activity Profiles

The pavement or bridge life-cycle activity profile refers to the frequency, timing, and magnitude of construction, rehabilitation, and maintenance treatments within its life cycle. A typical life-cycle activity profile represents the most cost-effective way of implementing strategically coordinated treatments to achieve the intended service life. In practice, pavement life-cycle activity profiles are determined using preset time intervals for treatments and condition triggers for treatments, respectively. Many state transportation agencies currently use preset time intervals because of lacking consensus in condition trigger values and consistency in pavement condition data. With respect to bridge life-cycle activity profiles, the preset time interval approach is also commonly used. Table 1 lists the typical frequency and timing of major treatments in pavement and bridge service lives used by the FHwA, American Association of State Highway and Transportation Officials (AASHTO), and state transportation agencies (FHwA 1987, 1991; Gion et al. 1993; INDOT 2002; AASHTO 2003).
Table 1. Typical Frequency and Timing of Major Treatments in Pavement and Bridge Life Cycles
FacilityMaterial typeService life (year)Treatment frequencyTiming
PavementFlexible 40Thin overlay +15th year
    Thick HMA overlay30th or 33rd year
 Rigid 40PCC joint sealing +7th year
    PCC joint sealing +15th year
    PCC repair techniques +23rd year
    Thick HMA overlay +30th year
    HMA crack sealing37th year
    PCC overlay +30th year
    PCC joint sealing35th year
BridgeConcreteChannel beam35Deck rehabilitation20th year
  T-beam/girder70Deck rehabilitation +20th, 55th year
    Superstructure replacement35th year
  Slab60Deck rehabilitation30th, 45th year
 PrestressedBox-beam65Deck rehabilitation +20th, 50th year
 Concrete  Deck replacement35th year
  Box girder50Deck rehabilitation20th, 35th year
 SteelBox-beam/girder70Deck rehabilitation +20th, 55th year
    Deck replacement35th year
  Truss80Deck rehabilitation +25th, 65th year
    deck replacement40th year
The life-cycle agency costs for each type of pavements or bridges can be quantified on the basis of the proposed life-cycle activity profile as Table 1. For a specific pavement or bridge project, the construction, rehabilitation, and maintenance costs in the pavement or bridge life-cycle can be estimated using historical data on the unit rates of construction, rehabilitation, and maintenance treatments multiplied by the project size. A geometric growth rate represented by a constant percent of annual growth can be used to establish annual routine maintenance costs for future years based on the first year routine maintenance cost within an interval between two major treatments.

Pavement/Bridge Life-Cycle Annual User Cost Profiles

For each user cost component, the first-year user costs under normal operation conditions within an interval between two major treatments can be calculated. A geometric growth rate can be used for estimating annual user costs in future years within the same interval based on the first year user costs. The excessive user costs caused by project work zones such as delay costs need to be considered for the year involving major treatments.

Estimation of Project Level Life-Cycle Benefits

The typical life-cycle activity profile for pavements or bridges represents the most cost-effective investment strategy to manage pavement or bridge facilities. If any needed treatment fails to be timely implemented as per the typical life-cycle activity profile, an early termination of the service life is expected. As such, the typical life-cycle activity profile can be used as the base case activity profile and the case with early service-life termination can be considered as an alternative case activity profile. For each type of pavements or bridges, the reduction in life-cycle agency costs of the base case activity profile compared with the alternative case activity profile can be computed as project level life-cycle agency benefits of timing implementing the needed project. Similarly, the decrease in life-cycle user costs according to the base case activity profile against the alternative case activity profile can be estimated as the project level life-cycle user benefits.
Fig. 1 illustrates an example of base case and alternative case activity profiles for the steel-box beam bridge and the method for estimating project level life-cycle agency benefits and user benefits by keeping the typical life-cycle activity profile for the bridge. For the base case life-cycle activity profile, agency costs in the T -year bridge service life consist of initial bridge construction cost CCON in year 0, first deck rehabilitation cost CDECKREH1 in year t1 , deck replacement cost CDECKREP in year t2 , second deck rehabilitation cost CDECKREH2 in year t3 , and annual routine maintenance costs. The annual routine maintenance costs between two major treatments in the bridge life-cycle will gradually increase over time due to the combined effect of higher traffic demand, aging materials, climate conditions, and other nonload related factors. Different geometric gradient growth rates are used for intervals between year 0 and t1 , t1 and t2 , t2 and t3 , and t3 and T , respectively.
Fig. 1. Illustration of base case and alternative case life cycles for the steel-box beam bridge
For the alternative life-cycle activity profile, it is assumed that the deck replacement project (with the cost of CPROJECT ) is actually implemented y1 years after year t2 as the base case profile, namely, CDECKREP in year t2 is replaced by CPROJECT in year t2+y1 . This will defer the second deck rehabilitation by y1 years. Because of postponing deck replacement and the second deck rehabilitation, the bridge service life may experience an early termination of y2 years. As for the annual routine maintenance costs, different geometric gradient growth rates are used for intervals between year 0 and t1 , t1 and t2+y1 , t2+y1 and t3+y1 , and t3+y1 and Ty2 , correspondingly. In particular, the annual routine maintenance cost profiles for the base case and alternative case profiles are identical from year 0 to year t2 . The project level life-cycle agency benefits are estimated as the reduction in bridge life-cycle agency costs quantified according to the base case activity profile compared with the alternative case activity profile.
The primary user cost items include vehicle operating costs, travel time, vehicle crashes, and vehicle air emissions. For each user cost item, the base case and alternative case annual user cost profiles in bridge life cycle follow a pattern similar to the profile of annual routine maintenance costs in bridge life cycle. In either the base case profile or alternative case profile, the first-year user cost amounts immediately after the major treatments including bridge construction, first deck rehabilitation, deck replacement, and second deck rehabilitation are directly computed on the basis of the unit user cost in constant dollars per vehicle mile of travel (VMT) and the annual VMT. The unit user cost per VMT is estimated according to average travel speed and roadway condition as in Li and Sinha (2004). Geometric growth rate is then applied to the first-year user cost amount for each interval between two major treatments to establish the annual user cost amounts for subsequent years within the interval. Additional work zone related costs are estimated using the procedures in FHwA (1998, 2000) and AASHTO (2003), and added to the annual user 235 cost amounts for the years in which major treatments are implemented. This ultimately establishes the base case and alternative case annual user cost profiles for vehicle operating costs, travel time, vehicle crashes, and vehicle air emissions, respectively.
For each user cost item, the annual user cost profiles for the base case and alternative case are identical from year 0 to t2 and are different for the remaining years in the bridge life cycle. The travel demand in terms of annual VMT for a specific year after year t2 could be different between the base case and alternative case due to the fact that the traffic volume, i.e., annual average daily traffic (AADT) and/or travel distance associated with the bridge might change for the two cases. The consumer surplus concept is employed to separately compute the user benefits by comparing the base case and alternative case annual user cost profiles for intervals from year t2 to t2+y1 , t2+y1 to t3 , t3 to t3+y2 , t3+y2 , Ty2 , and Ty2 to T . The total project level life-cycle user benefits are the aggregation of individual user benefit items associated with reductions in vehicle operating costs, travel time, vehicle crashes, and vehicle air emissions in the bridge life cycle. With equal weights assigned for agency benefits and user benefits, the total project level life-cycle benefits by keeping the typical life-cycle activity profile for the bridge are established by combining the two sets of benefits.

Estimation of Project Level Life-Cycle Benefits in Perpetuity

The project level life-cycle benefits in perpetuity can be quantified on the basis of the base case and alternative life-cycle activity profiles. As the base case life-cycle activity profile represents the most cost-effective investment strategy, investment decisions are always made with the intention to keep abreast of the base case life-cycle activity profile. For the base case life-cycle activity profile in perpetuity, the base case typical facility life-cycle is assumed to be repeated an infinite number of times. For the alternative case life-cycle activity profile in perpetuity, early termination of service life may occur in the first life-cycle, in the first and second life cycles or in the first several life cycles. After experiencing early service life terminations, the base case typical facility life cycle is expected to be resumed back for the subsequent life cycles in perpetuity horizon. This is because that the base case life-cycle profile represents the most cost-effective investment strategy that the decision maker always aims to achieve. Without loss of generality, the alternative case life-cycle profile in perpetuity in this study adopts early terminations for the first two life cycles and the base case life-cycle profile is used for subsequent life cycles in perpetuity horizon. The reduction in project level life-cycle agency costs between the base case and the alternative case life-cycle activity profiles in perpetuity is computed to establish project level life-cycle agency benefits in perpetuity.
Similarly, the reduction in project level life-cycle user costs between the base case and the alternative case life-cycle annual user cost profiles in perpetuity for vehicle operating costs, travel time, vehicle crashes, and vehicle air emissions can be separately computed and summed up to establish project level life-cycle user benefits in perpetuity. With equal weights considered for agency benefits and user benefits, they can be directly added to establish overall project level life-cycle benefits in perpetuity.

Risk Considerations in Estimating Project Level Life-Cycle Benefits

Primary Input Factors under Risk Considerations

Project construction, rehabilitation, and maintenance costs may not remain as predicted. Traffic demand may not follow the projected path. Discount rate may fluctuate over time during the pavement or bridge life-cycle. Such variations will in turn result in changes in the overall project level life-cycle benefits. In this study, the unit rates of project construction, rehabilitation, and maintenance treatments, traffic growth rates, and discount rates are primary input factors considered for probabilistic risk assessments.

Selection of Probability Distributions for the Input Factors under Risk Considerations

The minimum and maximum values of above input factors under risk considerations are bounded by non-negative values. For each of the risk factors, the distribution of its possible outcomes could be either symmetric or skewed. Such distribution characteristics can be readily modeled by the Beta distribution that is continuous over a finite range and also allows for virtually any degree of skewness and kurtosis. The Beta distribution has four parameters: lower bound (L), upper bound (H), and two shape parameters α and β , with density function given by
f(xα,β,L,H)=Γ(α+β)(xL)α1(Hx)β1Γ(α)Γ(β)(HL)α+β1(LxH)
(1)
where the Г -functions serve to normalize the distribution so that the area under the density function from L to H is exactly 1.
The mean and variance of the Beta distribution are given as
μ=αα+βandσ2=αβ(α+β)2(α+β+1)
(2)

Using Simulation for Probabilistic Risk Assessments

Simulation is essentially a rigorous extension of sensitivity analysis that uses randomly sampled values from the input probability distribution to calculate discrete outputs. Two types of sampling techniques are commonly used to perform simulations. The first type is the Monte Carlo sampling technique that uses random numbers to select values from the probability distribution. The second type is the Latin Hypercube sampling technique where the probability scale of the cumulative distribution curve is divided into an equal number of probability ranges. The number of ranges used is equal to the number of iterations performed in the simulation. The Latin Hypercube sampling technique is likely to achieve convergence in fewer iterations as compared to those of the Monte Carlo sampling technique (FHwA 1998).

Uncertainty Considerations in Estimating Project Level Life-Cycle Benefits

As a practical matter, the input factors under risk considerations may not be readily characterized using reliable probability distributions. Consequently, a meaningful mathematical expectation for each factor cannot be established and this invalidates risk-based analysis. Shackle’s model introduced herein is well suited to handle each input factor under uncertainty where no probability distribution can be readily established for a number of possible outcomes (Shackle 1949).
In general, Shackle’s model overcomes the limitation of inability to establish the mathematical expectation of possible outcomes of each input factor for project level life-cycle benefit/cost analysis according to the following procedure. First, it uses degree of surprise as a measure of uncertainty associated with the possible outcomes in place of probability distribution. Then, it introduces a priority index by jointly evaluating each known outcome and the associated degree of surprise pair. Next, it identifies two outcomes of the input factor maintaining the maximum priority indices, one on the gain side and the other on the loss side from the expected outcome X(E) . The expected outcome could be the average value or the mode of all known possible outcomes, but it is not the mathematical expectation as outcome probabilities are unknown. The two outcomes need to be standardized to remove the associated degrees of surprise. The absolute deviations of two outcomes relative to the expected outcome are terms as standardized focus gain xSFG and standardized focus loss xSFL from the expected outcome X(E) . This model yields a triple xSFL,X(E),xSFG for each input factor under uncertainty. More details of Shackle’s model are in Ford and Ghose (1998), Young (2001), and Li and Sinha (2004, 2006).
To simplify the application of Shackle’s model for uncertainty-based analysis, the grand average of simulation outputs from multiple iterations of replicated simulation runs can be used as the expected outcome X(E) for an input factor under uncertainty
X(E)=m=1Mn=1NXiM×N
(3)
where Xi=a simulation output representing a possible outcome; N=number of iterations in each simulation run; and M=number of replicated simulation runs. If higher valued outcomes are preferred for an input factor, the absolute deviation of the average value of simulation outputs that are lower than the expected outcome can used as standardized focus loss value xSFL and the absolute deviation of the average value of simulation outputs that are equal or higher than the expected outcome can used as standardized focus gain value xSFG for the input factor under uncertainty
xSFL=m=1Mn=1NrXiM×NrX(E)
(4)
xSFG=m=1M(n=1NXin=1NrXi)M×(NNr)X(E)
(5)
where Nr=number of simulation outputs in the rth simulation run such that Xi<X(E) if a higher outcome value is preferred for the input factor.
In some cases, lower outcome values are preferred for an input factor such as the discount rate. The Nr for computing the standardized focus loss value xSFL and the standardized focus gain value xSFG thus refers to number of simulation outputs in the rth simulation run such that Xi>X(E) .
As an extension of Shackle’s model dealing with the input factor under uncertainty, a decision rule is introduced to help compute a single value X for the input factor based on the triple xSFL,X(E),xSFG that can be used for estimating project benefits. Assuming that the decision-maker only tolerates loss from the expected outcome for the input factor under uncertainty by ΔX and if higher outcome values are preferred, the decision rule is set as
X={X(E),ifxSFLΔXX(E)xSFL[1ΔXX(E)],otherwise}
(6)
When lower outcome values are preferred for an input factor, the decision rule is revised to
X={X(E),ifxSFLΔXX(E)+xSFL[1+ΔXX(E)],otherwise}
(7)
If the standardized focus loss xSFL from the expected outcome X(E) does not exceed ΔX , the expected outcome value will be utilized for the input factor for the computation. This will produce an identical input factor value for both uncertainty-based and risk-based analyses. If the standardized focus loss xSFL from the expected outcome X(E) exceeds ΔX , a penalty is applied to derive a unique value for the input factor. Different tolerance levels ΔX ’s may be used for different input factors under uncertainty.

Generalized Framework for Uncertainty-Based Project Level Life-Cycle Benefit/Cost Analysis

Fig. 2 shows a generalized framework for uncertainty-based highway project level life-cycle benefit/cost analysis with input factors under mixed cases of certainty (the input factor is purely deterministic with single value), risk (the input factor has a number of possible outcomes with a known probability distribution), and uncertainty (the input factor has a number of possible outcomes with unknown probabilities). If an input factor is under certainty, the single value of the factor can be used for the computation. If an input factor is under risk, the mathematical expectation of the factor can be utilized for the computation. If an input factor is under uncertainty, the single value of the factor determined according to the decision rule extended from Shackle’s model can be adopted for the computation.
Fig. 2. Generalized framework for estimating overall project benefits under certainty, risk, and uncertainty
By using values of input factors determined under certainty, risk or uncertainty, the proposed framework helps establish project level life-cycle agency benefits and user benefits concerning decrease in agency costs, reduction in vehicle operating costs, shortening of travel time, decrease in vehicle crashes, and cutback of vehicle air emissions in perpetuity horizon, respectively. The combination of certainty, risk, and uncertainty cases for input factors may vary by project benefit item for the same project and may also vary for different types of highway projects.

Impacts of the Proposed Methodology on Estimating Project Benefits

Comparison of Estimated Project Benefits for Project Level Impact Assessments

Project level impact assessments compare project level life-cycle benefits separately estimated using the deterministic, risk-based, and the uncertainty-based project level life-cycle cost analysis approaches. For the application of deterministic project level life-cycle benefit/cost analysis, project benefits are calculated by assuming that all input factors are under certainty and each input factor has a single value. These values are directly used for the computation.
For the application of risk-based project level life-cycle benefit/cost analysis, project benefits are calculated by assuming that input factors regarding unit rates of construction, rehabilitation, and maintenance treatments, traffic growth rates, and discount rates are all under risk. The remaining input factors such as pavement or bridge service life and timing of treatments are still treated as being under certainty with single values. Monte Carlo simulations are executed to establish the grand average values of simulation outputs as mathematical expectations of input factors under risk. The single values of input factors under certainty and the grand average values of input factors under risk are used for the computation.
For the application of the uncertainty-based methodology, project benefits are calculated by assuming that the input factors regarding unit rates of construction, rehabilitation, and maintenance treatments, traffic growth rates, and discount rates are all under uncertainty or under mixed cases of risk and uncertainty. The remaining input factors are still considered under certainty with single values. For the input factor under risk, the grand average value as the mathematical expectation is established using Monte Carlo simulation outputs. For the input factor under uncertainty, the grand average value of simulation outputs is adjusted according to the preset decision rule. The single values of input factors under certainty, the grand average values of input factors under risk, and the adjusted grand average values of input factors under uncertainty are used for the computation.

Comparison of Project Selection for Network Level Impact Assessments

In order to assess the network level impacts of adopting different approaches for project benefit estimation, the three sets of project benefits computed using the deterministic, risk-based, and uncertainty-based approaches are separately applied to a stochastic optimization model for network-level project selection. The network level impacts are assessed by cross comparison of the overall benefits of selected projects and consistency matching rates of project selection using the three different sets of project benefits with the actual project selection and programming practice. This section briefly discusses the stochastic optimization formulation for finding the optimal subset of highway projects from all candidate projects to achieve maximized total project level life-cycle benefits where there is stochasticity in the available budget.
Consider a state transportation agency that carries out highway network-level project selection over a future project implementation period of tΩ years. The agency makes first round of investment decisions many years prior to project implementation using an estimated budget for all years. With time elapsing, updated budget information on the first few years of the multiyear project selection and programming period becomes available that motivates the agency to refine the investment decisions. In each refined decision-making process, the annual budget for the first few years that can be accurately determined is treated as a deterministic value, while the budget for the remaining years without accurate information is still handled as a stochastic budget.

Stochastic Optimization Model for Network Level Project Selection

Assuming that the multiyear budget is refined Ω times and each time an increasing number of years with accurate budget information from the first year is obtained, Ω -decision stages are therefore involved. Without loss of generality, a discrete probability distribution of budget possibilities can be assumed for each year where no accurate budget information is available. The stochastic model with - Ω -stage budget recourses can be formulated as a deterministic equivalent program that combines first stage decisions using the initial budget estimate with expected values of recourse functions for the remaining (Ω-1) stages (Birge and Louveaux 1997).
The stochastic model with Ω -stage budget recourses formulated as a Multichoice Multidimensional Knapsack Problem (MCMDKP) is shown below:
MaximizeAT.X1+ω=2ΩEξω[Qω(Xω(p),ξω)]
(8)
Stage 1
SubjecttoCktT.X1E(Bkt1)
(9)
X1=decision vector with 01 integer elements.
Stage 2
Eξ2[Q2(X2(p),ξ2)]=max{AT.X2(p)Bkt2(p)=E(Bkt2)}
(10)
SubjecttoCktT.X2(p)Bkt2(p)
(11)
X1+X2(p)1
(12)
X1 and X2(p)=decision vectors with 01 integer elements.
Stage L
EξL[QL(XL(p),ξL)]=max{AT.XL(p)BktL(p)=E(BktL)}
(13)
SubjecttoCktT.XL(p)BktL(p)
(14)
X1+X2(p)++XL(p)1
(15)
X1,X2(p),,XL(p)=decision vectors with 01 integer elements.
Stage Ω
EξΩ[QΩ(XΩ(p),ξΩ)]=max{AT.xΩ(p)BktΩ(p)=E(BktΩ)}
(16)
SubjecttoCktT.XΩ(p)BktΩ(p)
(17)
X1+X2(p)++XL(p)++XΩ(p)1
(18)
X1,X2(p),,XL(p),,XΩ(p)=decision vectors with 01 integer elements.
In the objective function as Eq. (8), the first term is for the overall project benefits in the first stage decisions using the initial budget estimate and the second term is for the expected value of overall project benefits for the remaining (Ω1) stages with budget recourse decisions. Eqs. (9) - (17) are budget constraints for the optimization model. The constraints are imposed by highway asset management program k and by project implementation year t for the multiyear project implementation period. Eqs. (10) - (16) compute the expected values of maximized total project benefits under budget constraints at stages 2, L , and Ω , respectively. At each stage, the budget used for the computation is determined by selecting a possible budget that maintains the least amount of deviations from the expected budget for that stage. For instance, the stage L budget for the computation is chosen from all possible budgets BktL(1) , BktL(2),,BktL(pL) such that Σk=1KΣt=1M[BktL(p)E(BktL)]2 is minimized. Eqs. (12) - (18) ensure that each highway project can be selected at most once in the multistage recourse decision process.
In the above optimization model, budget constraints are imposed by highway asset management program k and by year t for the multiyear project implementation period. In current practices, available budgets for different highway asset management programs such as pavement preservation program and bridge preservation program are not transferable across programs. For each asset management program, the yearly constrained multiyear budgets in the optimization model can however be treated as a cumulative budget for all years combined. For the option of cumulative budget constraints, the notations BktL(p) , p(BktL(p)) , and E(BktL) in Eqs. (9) - (17) are replaced by Σt=1MBktL(p) , p(Σt=1MBktL(p)) and E(Σt=1MBktL) , where E(Σt=1MBktL)=Σp=1p1[p(Σt=1MBktL(p))Σt=1MBktL(p)] (L=1,2,,Ω) . In this study, the optimization model is applied using the yearly constrained budget scenario and the cumulative budget scenario, respectively.

Model Solution

For the purpose of the present paper, the solution algorithm developed based on the LaGrangian relaxation technique was implemented using a customized computer code.

Case Study

A case study was conducted to examine the impacts of using deterministic, risk-based, and uncertainty-based project level life-cycle cost analysis approaches on computing the benefits of individual highway projects. The computed project benefits were used to assess the network-level impacts of adopting different project level life-cycle cost analysis approaches on project selection results.

Data Sources

Data Collection and Processing for Highway Project Benefit Estimation

For assessing the project level impacts of using deterministic, risk-based, and uncertainty based project level life-cycle cost analysis approaches for project benefit estimation, historical data on the Indiana state highways for period 1990–2006 were collected to establish the base case life-cycle activity profiles and annual user cost profiles for different types of pavements and bridges. The data items collected mainly included project type and size; unit rates of construction, rehabilitation, and maintenance treatments; unite rates of vehicle operating costs, travel time, crashes, and air emissions; traffic volume and growth rates; discount rates, etc. Table 2 presents Beta distribution parameters established for those factors on the basis of historical data.
Table 2. Input Values of Factors for Risk- and Uncertainty-Based Project Benefit Analysis
Input factorsMeanStandarddeviationBeta distribution parameters
L H α β
Flexible pavement costConstruction1,353,537694,614588,3853,165,8402.494.50
(1990, $/lane-mile)Rehabilitation155,287509,87929,1471,119,8632.564.50
 Resurfacing52,93819,68926,364101,6022.564.50
 Routine maintenance13849942,1862.274.50
Rigid pavement costConstruction1,334,841763,709674,2992,947,1732.254.50
(1990, $/lane-mile)Rehabilitation383,704242,26057,9522,052,8962.414.50
 Routine maintenance32320441,9813.104.50
All pavement costPreventive4,1206,54418621,9992.564.50
(1990, $/lane-mile)maintenance      
Concrete bridge costDeck62420.13872.394.50
(1990, $ft2 )Superstructure110820.23722.394.50
 Substructure115920.13722.394.50
Steel bridge costDeck86590.47342.174.50
(1990, $ft2 )Superstructure171750.47342.174.50
 Substructure206990.47342.174.50
Annual routine maintenance growth3%1%1%5%4.504.50
Annual traffic growth2%1%1%3%4.504.50
Discount rate4%1%3%5%4.504.50
Furthermore, 11-year data on 7,380 candidate projects (grouped into 5,068 contracts) proposed for Indiana state highway programming during 1996–2006 were collected for applying the deterministic, risk-based, and uncertainty based project level life-cycle cost analysis approaches for project benefit estimation. For each pavement or bridge project, base case and alternative case life-cycle activity profiles and annual user cost profiles were established. As described in the proposed methodology, the agency benefits and user benefits associated with reduction in vehicle operating costs, shortening of travel time, decrease in vehicle crashes, and cutback of vehicle air emissions for each project were separately estimated by comparing the respective base case and alternative case life-cycle profiles. For the application of the deterministic life-cycle cost analysis approach, the single values of all input factors were utilized for estimating the project level life-cycle benefits.
For the application of the risk-based life-cycle cost analysis approach, Beta distribution parameter values for the input factors regarding unit rates of construction, rehabilitation, and maintenance treatments; traffic growth rates; and discount rates were applied in 10 simulation runs, each with 1,000 iterations using the RISK software, Version 4.5 (Palisade Corp., Ithaca, N.Y., 2007). The Latin Hypercube stratified sampling technique was used in the simulations to reach faster convergence. The grand average of simulation runs for each risk factor was adopted for computing the mathematical expectations of agency benefits and user benefits.
When conducting risk-based analysis, it was found that project benefits related to decrease in agency costs, reduction in vehicle operating costs, and cutback of vehicle air emissions were not very sensitive to the variations of simulation outputs of the input factors under risk. However, travel time and vehicle crash reductions varied considerably with the simulation outputs of the factors. For this reason, the project user benefits concerning travel time and vehicle crash reductions were further estimated using the uncertainty-based analysis approach. Specifically, the grand average values of simulation runs for unit rates of construction, rehabilitation, and maintenance treatments, traffic growth rates, and discount rates were adjusted according to the preset decision rules as the proposed methodology for uncertainty-based analysis. The adjusted values were used to compute the benefits of travel time and vehicle crash reductions under uncertainty.

Data Collection and Processing for Network Level Highway Project Selection

The three sets of project level life-cycle benefits estimated for the 7,380 candidate projects were used to assess the network-level impacts of using deterministic, risk-based, and uncertainty-based project level life-cycle cost analysis approaches for estimating project benefits on project selection results.
Additional data on available budgets by highway asset management program and by project implementation year for period 1996–2006 were collected. The annual average budget was approximately 700 million dollars with 4% increment per year. The initially estimated budget for the project implementation period was found to have being updated three times by the Indiana Department of Transportation (DOT). This provided 4-stage budget recourses in the application of the stochastic optimization model for project section. The budget adjustments were mainly made on pavement preservation, bridge preservation, system expansion, and maintenance programs, with changes varying from 32 to +60% .

Summary of Estimated Project Level Life-Cycle Benefits

Table 3 lists project level life-cycle benefits of some pavement and bridge projects. On average, the present worth amounts of project level life-cycle benefits estimated using deterministic, risk-based, and uncertainty-based analysis approaches for the 7,380 projects are 4.18, 7.14, and 6.64 million dollars per project (in 1990 constant dollars), respectively. The average benefit-to-cost ratios are 3.24, 5.54, and 5.16, correspondingly. The significant difference between the project benefits estimated using the deterministic analysis approach and risk-based analysis approach is mainly attributable to large standard deviations of input factors considered for probabilistic risk assessments. The comparable results of project benefits computed using the risk-based analysis approach and uncertainty-based analysis approach are intuitive. This is because the grand average of simulation outputs for each input factor under uncertainty is adjusted only if the deviation between the grand average as the expected outcome and standardized focus loss value exceeds the preset threshold level. The input factor values for risk-based and uncertainty-based analyses will be identical if no adjustment is made.
Table 3. Project Level Life-Cycle Benefits of Some Pavement and Bridge Projects Computed Using Deterministic, Risk-Based, and Uncertainty-Based Analysis Approaches (1990 Constant Dollars)
ContractNo.LetyearLanesLength(miles)AADTWork typeProjectcostProject benefits estimated under
CertaintyRiskUncertainty
12021200040.1169,200Bridge widening2,291,0006,959,43411,703,26411,703,264
12040200040.5032,630Pavement resurfacing4,620,0004,776,3196,927,6696,365,844
12077200022.063,170Pavement resurfacing3,000,0009,436,80415,545,50115,545,501
12158199923.7016,770Added travel lanes750,0003,036,2535,405,6214,806,134
20694199621.343,420Flexible pave. replace51,00043,704131,989131,989
21743199640.4025,310Pavement rehabilitation696,0001,271,5741,878,3751,878,375
217491998213.634,190Pavement resurfacing11,573,00038,024,31963,943,22563,943,225
21825199642.5311,150Pavement rehabilitation151,000504,5741,033,2741,505,738
21931199640.782,664Rigid pavement replace196,000705,235736,046736,046
21944199629.461,100Pavement rehabilitation131,000239,334353,545353,545
22026199620.158,291Bridge widening108,000267,380299,746254,516
22032199646.3012,274Pavement resurfacing754,0001,743,1882,753,2592,559,337
22044199621.1013,994Pavement resurfacing2,757,0006,169,0676,773,2425,702,627
22119199840.1027,700Pavement rehabilitation264,000445,933658,734658,734
22264199621.137,843Pavement resurfacing1,226,0003,566,5667,164,6116,450,209

Comparisons of Project Selection Results

Comparison of Total Benefits of Selected Projects

Fig. 3 illustrates the total benefits of projects selected using the optimization model based on three sets of estimated project benefits (deterministic, risk-based, and uncertainty-based), two types of budgets (deterministic and stochastic), and two budget constraint scenarios (yearly constrained and cumulative). Regardless of budget types and budget constraint scenarios, the total benefits of selected projects are the lowest for project benefits estimated using the deterministic analysis approach and are the highest for project benefits computed using the risk-based analysis approach.
Fig. 3. Comparison of total benefits of selected projects using deterministic and stochastic budgets under yearly constrained and cumulative budget scenarios (1996–2006)
Despite approaches used for computing project benefits and types of budgets used in the optimization model, the project selection using the cumulative budget scenario generally yielded higher total benefits. The results are not unexpected. The cumulative budget scenario does not have year-by-year budget restrictions as those added to the yearly constrained budget scenario. This entails more flexibility to the optimization model in conducting project selection, leading to increases in the total project benefits.

Comparison of Number of Selected Contracts

Table 4 presents the comparison of contracts selected using the three sets of project benefits, two types of budgets, and two budget constraint scenarios. The matching rates were established in reference to the contracts being authorized by the Indiana DOT. One match is counted if a contract is both selected in the optimization model application and also authorized by the Indiana DOT.
Table 4. Summary of Consistency in Contract Selection Results under Different Extents of Risk and Uncertainty Considerations
YearNo. ofcontractsIndianaDOTauthorizedYearly constrained budgetCumulative budgetAll methodsmatched withIndiana DOT
DeterministicStochasticDeterministicStochastic
MD MR MU MD MR MU MD MR MU MD MR MU No.%
199646444343339038843739039443941141443941241531972
199741235838733834438633634339037037239036937425070
199842927540835136340935336141337537741437737718768
199941132337632233338132233238835235238835135220363
200061057857650651657950451458254454458654654641672
200141841239534835839634335639336336339336036628970
200242242139934334339833934340237337340637337729169
200346946143737338144037137544441341444641341831568
200464964860851953161552152861257858061357858146371
200540840638033733938433734038735535938935736428269
200637637535530230735330330335733333635933433825969
Total5,0684,7004,7543,8714,2034,7783,8964,1894,8074,6254,4844,8234,6604,508  
Total Match with Indiana DOT4,4003,8283,8894,4213,8173,8774,4514,1294,1454,4664,1314,1683,274 
% Match with Indiana DOT94%81%83%94%81%82%95%88%88%95%88%89% 70%
Note: MD , MR , and MU=project benefits estimated using deterministic based, risk-based, and uncertainty-based analysis approaches, respectively.
For the deterministic budget, the average matching rates for the three sets of estimated project benefits and two budget constraint scenarios are 81–95%. Irrespective of using project benefits estimated by the deterministic, risk-based or uncertainty-based life-cycle cost analysis approach, the use of cumulative budget constraint scenario in the optimization model for project selection resulted in the selection of a higher number of contracts and with a higher matching rate. The net increases in the matching rates for the cumulative budget scenario as opposed to the yearly constrained budget scenario are 1% for deterministic project benefits, 7% for risk-based project benefits, and 5% for uncertainty-based project benefits, respectively. The relative increases in the matching rates resulted from the use of the cumulative budget scenario versus the yearly constrained budget scenario are 1%94%=1.1% for deterministic based project benefits, 7%81%=9% for risk-based project benefits, and 5%83%=6% for uncertainty-based project benefits, correspondingly.
For the stochastic budget, the average matching rates for the three sets of estimated project benefits and two budget constraint scenarios also range from 81–95%. The use of cumulative budget constraint scenario in the optimization model for project selection resulted in the selection of a higher number of contracts and with a higher matching rate. The increases in the matching rates for the cumulative budget scenario as opposed to the yearly constrained budget scenario are 1% for deterministic based project benefits, 7% for risk-based project benefits, and 7% for uncertainty-based project benefits, respectively. The relative increases in the matching rates are 1%94%=1.1% , 7%81%=9% , and 7%82%=8.5% , correspondingly.
Irrespective of budget types and budget constraint scenarios, the use of project benefits estimated by the deterministic life-cycle cost analysis approach for project selection produced a higher percentage of matching rate as compared to matching rates established for project benefits estimated by risk-based and uncertainty-based analysis approaches. The matching rates for project benefits estimated using the uncertainty-based analysis approach are slightly higher than those of the project benefits computed by the risk-based analysis approach. In particular, increases in the matching rates are 2% for yearly constrained deterministic budget, 2% for yearly constrained stochastic budget, 0 percent for cumulative deterministic budget, and 1 percent for cumulative stochastic budget, respectively. The relative increases in the matching rates are 2%81%=2.5% , 2%81%=2.5% , 0%82%=0% , and 1%88%=1.1% , accordingly.
Without regard to using different approaches for project benefit estimation and employing different types of budgets and budget constraint scenarios in the optimization model for project selection, the average matching rate between projects selected using the optimization model and actually authorized by the Indiana DOT for the eleven-year analysis period is 70%. After removing this portion of matching rate invariant to approaches used for project benefit analysis and types of budgets and budget constraint scenarios used in the optimization model for project selection, the relative increases in the matching rates of project selection resulted from the use of uncertainty-based analysis approach versus the risk-based analysis approach for project benefit estimation are 2%(81%70%)=18% for yearly constrained deterministic budget, 2%(81%70%)=18% for yearly constrained stochastic budget, 0%(82%70%)=0% for cumulative deterministic budget, and 1%(88%70%)=9% for cumulative stochastic budget, accordingly.

Summary and Conclusions

This paper proposed a new methodology for highway project level life-cycle benefit/cost analysis that handles certainty, risk, and uncertainty inherited with input factors for the computation. A case study was conducted to assess the impacts of risk and uncertainty considerations in estimating project level life-cycle benefits and on the results of network-level project selection.
The case study results revealed that using project level life-cycle benefits estimated by the proposed uncertainty-base analysis approach yielded a higher percentage of matching rate with the actual programming practice of the Indiana DOT as compared to the matching rate of using the project benefits computed by the risk-based analysis approach. The relative increase in matching rate with uncertainty considerations is up to 2.5%. After removing the portion of matching rate invariant to approaches used for project benefit estimation and types of budgets and budget constraint scenarios considered in the optimization model for project selection, the relative increase in the matching rate is as high as 18 percent. The difference is quite significant. The proposed methodology offers a means for transportation agencies to explicitly address uncertainty issues in project level life-cycle benefit/cost analysis that would enhance the existing risk-based life-cycle cost analysis approach.

Recommendations

Application of the proposed methodology requires collecting a large amount of data. This may limit the methodology application primarily to state and large-scale local transportation agencies that maintain sufficient historical data on highway system preservation, expansion, operations, and expenditures. In addition, the customized Beta distribution parameters need to be updated over time to reflect changes in the values of input factors for the analysis. Moreover, the equally assigned weights for project level life-cycle agency benefits and user benefits may be adjusted to assess the impact of such changes on the estimation of project benefits and on the results of network-level project selection.

Notation

The following symbols are used in this paper:
A
=
vector of benefits of N contracts, A=[a1,a2,,aN] ;
ai
=
benefits of contract i , i=1,2,,N ;
BktL(p)
=
pth possibility of budget for program category k in year t in stage L ;
Ckt
=
vector of costs of N contracts using budget from program category k in year t , Ckt=[c1kt,c2kt,cNkt] ;
cikt
=
cost of contract i using budget from highway asset management program k in year t ;
E(BktL)
=
expected budget in stage L , where E(BktL)=Σp=1PL[P(BktL(p))BktL(p)] ;
EξL[Q(XL(p),ξL)]
=
mathematical expectation of the recourse function in stage L ;
i
=
1,2,,N ;
k
=
1,2,,K for different highway asset management programs, which typically include bridge preservation, pavement preservation, safety improvements, roadside improvements, system expansion, state park highway facilities, ITS installations, and maintenance programs;
L
=
1,2,,Ω ;
p
=
1,2,,pL , where pL=sL.sL+1.sΩ ;
p(BktL(p))
=
probability of having budget scenario BktL(p) occur in stage L ;
Q(XL(p),ξL)
=
recourse function in stage L ;
t
=
1,2,,M .
XL(p)
=
decision vector using budget BktL(p) in stage L , XL(p)=[X1,X2,,XN] ;
xi
=
decision variable of contract i , i=1,2,,N , (xi=01) ; and
ξL
=
randomness associated with budget and decision space comprised of all possible combinations in the values of decision variables in stage L .

Acknowledgments

The writers acknowledge the Midwest Regional Universities Transportation Center at the University of Wisconsin—Madison for financial support of this research. This paper does not necessarily reflect the views of project sponsor.

References

AASHTO. (2003). User benefit analysis for highways, American Association of State Highway and Transportation Officials, Washington, D.C.
Abaza, K. A. (2002). “Optimum flexible pavement life-cycle analysis model.” J. Transp. Eng., 128(6), 542–549.
Birge, J. R., and Louveaux, F. (1997). Introduction to stochastic programming, Springer, New York.
Chan, A., Keoleian, G., and Gabler, E. (2008). “Evaluation of life-cycle cost analysis practices used by the Michigan Department of Transportation.” J. Transp. Eng., 134(6), 236–245.
Falls, L. C., and Tighe, S. (2003). “Improving LCCA through the development of cost models using the Alberta roadway rehabilitation and maintenance analysis.” Transportation Research Board 82nd Annual Meeting, Washington, D.C.
Federal Highway Administration (FHwA). (1987). Bridge management systems, Demonstration Project No. 71. U.S. Department of Transportation, Washington, D.C.
Federal Highway Administration (FHwA). (1991). Pavement management systems, U.S. Department of Transportation, Washington, D.C.
Federal Highway Administration (FHwA). (1998). Life-cycle cost analysis in pavement design—In search of better investment decisions, U.S. Department of Transportation, Washington, D.C.
Federal Highway Administration (FHwA). (2000). Highway economic requirements system, U.S. Department of Transportation, Washington, D.C.
Ford, J. L., and Ghose, S. (1998). “Lottery designs to discriminate between Shackle’s theory: Expected utility theory and non-expected utility theories.” Annals of operations research, Kluwer Academic, Norwell, Mass.
Gion, L. C., Gough, J., Sinha, K. C., and Woods, R. E. (1993). User’s manual for the implementation of the Indiana bridge management system, Purdue Univ., West Lafayette, Ind.
Hawk, H. (2003). “Bridge life cycle cost analysis.” NCHRP Rep. 483, National Academies, Washington, D.C.
Hicks, R. G., and Epps, J. A. (1999). “Life cycle cost analysis of asphalt-rubber paving materials.” Final Rep. Volumes 1 and II to Rubber Pavements Association, Tempe, Ariz.
INDOT. (2002). Indiana Department of Transportation design manual, Indiana Department of Transportation, Indianapolis.
Labi, S., and Sinha, K. C. (2005). “Life-cycle evaluation of flexible pavement preventive maintenance.” J. Transp. Eng., 131(10), 744–751.
Li, Z., and Sinha, K. C. (2004). “A methodology for multicriteria decision making in highway asset management.” Transportation Research Record. 1885, 79–87.
Li, Z., and Sinha, K. C. (2006). “Application of Shackle’s model for highway project evaluation under uncertainty.” Proc., 9th Int. Conf. on Applications of Advanced Technology in Transportation, Chicago, ASCE, Reston, Va., 67–73.
Mohammadi, J., Guralnick, S. A., and Yan, L. (1995). “Incorporating life-cycle costs in highway-bridge planning and design.” J. Transp. Eng., 121(5), 417–424.
Peshkin, D. G., Hoerner, T. E., and Zimmerman, K. A. (2005). “Optimal timing of pavement preventive maintenance treatment applications.” NCHRP Rep. No. 523, National Academies, Washington, D.C.
Purvis, R. L., Babaei, K., Clear, K. C., and Markow, M. J. (1994). “Life-cycle cost analysis for protection and rehabilitation of concrete bridges relative to reinforcement corrosion.” SHRP-S-377, National Academies, Washington, D.C.
Reigle, J. A., and Zaniewski, J. P. (2005). “Risk-based life-cycle cost analysis for project-level pavement management.” Transportation Research Record. 1816, 34–42.
Setunge, S., et al. (2005). “Whole of life cycle cost analysis in bridge rehabilitation.” Rep. No. 2002-005-C-03, RMIT Univ., Melbourne, Australia.
Shackle, G. L. S. (1949). Expectation in economics, 2nd Ed., Cambridge University Press, Cambridge, U.K.
Tighe, S. (2001). “Guidelines for probabilistic pavement life cycle cost analysis.” Transportation Research Record. 1769, 28–38.
Wilde, J. W., Waalkes, S., and Harrison, R. (1999). “Life cycle cost analysis of Portland cement concrete pavements.” FHWA/TX-00/0-1739-1, Univ. of Texas at Austin, Tex.
Young, R. A. (2001). Uncertainty and the environment: Implications for decision making and environmental policy, Edward Elgar, Cheltenham, U.K.

Information & Authors

Information

Published In

Go to Journal of Transportation Engineering
Journal of Transportation Engineering
Volume 135Issue 8August 2009
Pages: 516 - 526

History

Received: May 20, 2008
Accepted: Nov 10, 2008
Published online: Jul 15, 2009
Published in print: Aug 2009

Permissions

Request permissions for this article.

Authors

Affiliations

Assistant Professor, Dept. of Civil, Architectural and Environmental Engineering, Illinois Institute of Technology, Chicago, IL 60616 (corresponding author). E-mail address: [email protected]
Sunil Madanu
Graduate Research Assistant, Dept. of Civil, Architectural and Environmental Engineering, Illinois Institute of Technology, Chicago, IL 60616.

Metrics & Citations

Metrics

Citations

Download citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited by

View Options

Media

Figures

Other

Tables

Share

Share

Copy the content Link

Share with email

Email a colleague

Share