Abstract

Pavement maintenance and rehabilitation (M&R) plan for maintaining the pavement quality in an acceptable level has direct influence on the required budget. Deterministic budgeting is an unrealistic assumption, so, in this study, a two-stage stochastic model using integer programming is developed to address uncertainty in budgeting. Another aim of this study is to develop an executive model that considers a broad range of parameters at network level maintenance and rehabilitation planning. While having too many details in planning problems makes them more complicated, some restrictions called “technical constraints” were considered to reduce solution time of solving procedure as well as improve M&R activities assignment efficiency. Comparing results of the stochastic model with a deterministic model for a case study revealed that the two-stage stochastic model led to increased total cost compared to the deterministic one due to considering probability in budgeting. However, the developed model provides several M&R plans that are compatible with budget variation.

1. Introduction

Network level pavement management is complex due to the influential factors such as network size, available budget, pavement performance, distress and maintenance and rehabilitation (M&R) history, traffic and climatic conditions, nonuniform structural properties of pavement and construction operation [13]. Each of the aforementioned factors changes during the years, so it is crucial to estimate them within the life-span of pavement. On the other hand, the deviations between the predicted and actual values are inevitable. Thus, uncertainty in the predictions of the influential factors, in addition to their values change with elapse of time, makes the pavement management problems further complex.

Among the influential factors in Pavement Management System (PMS), the one dealing with the limited budget is of great importance. Taking limited budget into account, Wang, Zhang, and Machemehl proposed an optimization process based on the linear integer programming to select the best M&R strategies at the network level for a 5-year period of planning [4]. Wu, Flintsch, and Chowdhury focused on development of a decision-making model to optimize budgeting of the short-term preventive M&R actions [5]. They utilized goal planning methods for multiobjective problems and Analytic Hierarchy Process (AHP) for prioritization under multicriteria. The model was to achieve two objectives: minimizing the M&R costs and maximizing the pavement quality during its service life. Priya, Srinivasan, and Veeraragavan formulated the problems regarding determination of the appropriate types of M&R actions and the time for implementation of them using an effective optimization method [6]. Jesus et al. used a simple but practical model to optimize M&R management with limited budget using linear integer programing [7]. The model determined the budget assignments required to achieve the predetermined goals. Chai et al. established a correlation between the budget assignments and pavement quality indices for a road network in Queensland and, based on the correlation, a model developed [8]. In their study, an index, namely, pavement sustainability index, indicating budget deficit was introduced.

Single-objective optimization for M&R work planning deals with a function in which one independent variable is optimized. However, its main drawback is its suboptimal response compared to that obtained from multiobjective planning [9]. Nevertheless, various pavement management studies have developed and utilized multiobjective optimization models. A research developed a biobjective deterministic optimization model which simultaneously satisfied the objectives of both minimization of total maintenance costs and maximization of performance of the road network [10]. Also, an optimization methodology for county paved roads has been devised that identifies the best mix of preservation projects within budget, maximizing traffic (passengers and trucks traffic) on treated roads, maximizing the weighted average PSI, and minimizing risk [11].

There are also various optimization approaches for developing decision-making models in PMS. The Analytic Hierarchy Process (AHP) is the simplest and most popular method for multiobjective M&R action prioritization [12, 13]. Others are the life-cycle cost analysis for probabilistic model using fuzzy logic [14] and the application of Artificial Neural Network (ANN) for solving multiobjective optimization problems in PMS [1520]. Moreover, an optimization tool based on a hybrid Greedy Randomized Adaptive Search Procedure (GRASP) has been developed. This tool facilitates the design of optimal maintenance programs subject to budgetary and technical restrictions and explores the effects of different budgetary scenarios on overall network conditions [21].

Budget restrictions have always been the main challenge of PMS to deal with. The PMS models can be developed based on available budget, with the assumption that the available budget is equal to the predicted budget. However, due to the political and economic risks involved in any financial system, there is uncertainty in budget allocations to organizations. Lack of budget can lead to postponing M&R actions implementation in some years during the planning period, which negatively affects the whole predictions. Therefore, deterministic budget cannot be a realistic assumption for planning purposes.

Most studies on PMS have considered pavement deterioration as nondeterministic variable [2326]. Durango-Cohen took pavement deterioration models as nondeterministic models and evaluated the M&R policies efficiently when imprecise pavement deterioration models are available [27]. Some researchers found the uncertainty in pavement deterioration models as function of uncertainty in pavement structural design, traffic and climatic conditions, and pavement age [2832]. However, it is important not to restrict the uncertainty to the pavement deterioration [23]. Kuhn [33] considered the uncertainty in the return of investment and Chootinan et al. [34] took the predicted traffic volumes as nondeterministic variable. Yet, there are few studies that considered future budget as uncertain [22, 35].

Gao and Zhang used robust optimization for budgeting M&R at project level considering uncertainty in the annual budget [36]. The approach efficiency was examined using a 20-year pavement performance data and it was found reasonable. Li developed a probabilistic model to allocate nondeterministic budget to the highway network [37]. The model was formulated as probabilistic multiobjective and multidimensional knapsack problem and the results indicated that the uncertainty has significant positive effect on the efficiency of investment.

It seems that the application of uncertainty in development of M&R planning is still immature [38]. Two-stage stochastic programming is useful tool to develop planning models [3944]. Uncertainty in budgeting justifies application of two-stage stochastic planning in PMS. In the present study, a model for M&R work planning is developed using two-stage stochastic approach with nondeterministic budget. It renders a multi-optional M&R solution under the nondeterministic budget circumstance. Moreover, at network level PMS, consideration of more decomposed condition indicators and consequently extension of details are important study areas that need further investigation in models development. In previous studies, due to the complexity of the problem, especially in large-scale issues or long-term planning periods, less attention has been paid to increasing the number of details. Although the dimensions of the networks expanded as the studies developed, the number of M&R treatments used is still constant and limited [23]. Unlike most models that use four M&R strategies at the most, in this research, eight M&R categories were defined using four pavement quality indicators reflecting thermal distress, structural, skid, and roughness conditions and another that reflects a specific combination of them. Clearly, this level of detail and its application in pavement M&R planning have led to increased model accuracy as well as a higher degree of problem difficulty and complexity. As such, some restrictions called “technical constraints” were defined to reduce solution time of solving procedure. This advantage is applicable for huge networks or long-term planning durations and helps avoid expending too much time in solving M&R planning issues. Therefore, the development of a linear integer two-stage stochastic programming model, which can be used for large-scale networks or for long-term planning periods while also addressing the issue of a greater number of details and considering nondeterministic budget parameter, is one of the innovations of this study.

2. Literature

In this section, a comparison of the applied methodologies mentioned above is summarized in Table 1. As can be seen, there are few studies that utilize more than one condition indicator and at the same time consider uncertainty of budgeting in pavement M&R planning. Accordingly, this study attempts to develop a two-stage stochastic model involving a larger number of condition indicators which allow the network level M&R work planning to become closer to project-level one.

3. Research Approach

In this research, formulation of integer programing at the network level is developed based on a two-stage stochastic model on the condition of nondeterministic budget. GAMS software was employed to solve the developed model using the CPLEX solver. A small network with 10 sections and a planning period of 3 years is considered as a case study. The results of deterministic and probabilistic models were compared for the case study to assess how budget uncertainty affects the results of M&R planning.

3.1. Model Development

Indices, parameters, and variables of the developed method were described in Tables 2, 3, and 4, respectively. t, n, m, and s are the main indices of the model corresponding to the year, section ID, M&R action ID, and budgeting scenario. Index b is an auxiliary index that was used to round the pavement quality indicator values with precision of 0.1. Notably, t addresses the end of a year in the planning period.

In accordance with Table 4, the goal of the model is to find the binary variables of x and y when, with application of the restrictions of the problem, the objective function is minimum. 1 represents implementation of an M&R action and 0 means no M&R application. Moreover, it can be found in Table 4 that there are four normalized indicators of (qf, qt, qs, and qr) as well as a compound index of qo to quantify pavement condition with respect to structural condition, thermal distresses, skid resistance, roughness, and a linear combination of them, respectively. Notably, the values of these indicators were normalized to be between 0 and 1.

As given in Table 5, there are 8 categories of M&R actions with an ID number attributed to them.

Generally, there is administration cost (AC) and user cost (UC) included in the work planning optimization models. The AC is M&R actions expenses that the administration pays while UC is the sum of vehicle operation cost (VOC), delay cost (DC), and crash cost (CC) that users of the road pay indirectly [45]. In this research, VOC, which depends on the pavement quality, and DC, which is due to treatment type and the time duration it blocks the road, were used in the developed model. It should be noted that a few efforts have been made in literature to determine how pavement condition affects accidents rate and crash costs and what these effects are. Therefore, in this study, due to absence of required data and models related to accidents rate and crash costs affected by pavement condition, the crash cost effects are not considered in the modeling.

Since the deviation in conditions in each scenario leads to the decision change, in the two-stage stochastic model of this article, a cost-based parameter (Ca) was defined to take the cost due to the difference between the first and second stage decision variables into account. Table 6 lists all types of the influential costs on objective functions that were utilized for cost analysis.

Equation (1) expresses the objective function (ob) where, in addition to minimization of all costs, maximization of pavement network condition at the end of analysis period is going to be achieved. Since the objective function is to be minimized, consumed pavement financial value was defined and considered in it, which should be minimized in order to have maximized pavement condition at the end of the analysis period. Simply, by multiplying a transformer term , the utility related parameter (va: pavement financial value) was converted to the cost related variable (Va: consumed pavement financial value) which caused the objective function to be formulated as (1) in the form of a single minimizing function. Actually, in this way, a multiobjective function has been changed to a single-objective function.Equations (2) to (45) present all of the constraints of the model. Equations (2) to (7) indicate the constraints of auxiliary variables of pavement condition. These equations convert the continuous values of the input parameters corresponding to condition indicators (between 0 and 1) into discretized values with increment of 0.1. Equations (8) to (11) define the constraints of pavement deterioration rate and yield the condition deterioration rate for each condition indictor.Subsequently, (12) shows the constraints of pavement condition at each year and was defined to compute the values of pavement condition indictors at the end of each year.Equation (13) is the constraint of the overall condition computation and gives the linear relation between the overall condition (qo) and condition indicators (q: qf, qt, qs and qr).With the aid of (14) to (16), the calculation of M&R costs based on distress intensity level is possible. As such, these equations define the constraints related to the intensity of distress for localized maintenance application.The constraint for budget allocation (17) relates to the restriction on the annual budget.Equation (18) is the constraint that determines the applications of localized maintenance or corrective rehabilitation which is done based on cost comparison. The equation imposes a restriction that recommend a corrective M&R action instead of several localized preventive M&R actions if the costs of the former is less than the latter.Equations (19) to (27) are used for user cost computations. The constraints of auxiliary variables for vehicle operation cost are defined by (19) to (23), and the constraints of vehicle operation cost calculation are defined by (24) to (27).Equations (28) and (29) are the constraints to avoid the application of more than one M&R action from a M&R policy during each year. These restrictions do not allow applying more than one strategy from Table 5 within a year.Equation (30) is the constraint to correlate between first and second stage decision variables.As well, extra constraints are defined according to the existence criteria in Table 7 which are named technical constraints ((31) to (44)). These constraints determine the appropriate limits with using some threshold values which create the zone for ignoring M&R action.In order to have realistic M&R planning, since each scenario is based on the variation in annual budget, it is crucial to have the same M&R actions for those scenarios with similar budgets in past years. In other words, budget of each year is effective in the M&R assignment of the same and coming years and has no influence on the past years’ M&R assignment. Accordingly, (45) presents the restriction titled “nonanticipative constraint in scenarios” (this formulation is developed based on sequence in budgeting order given in Table 9).

3.2. Problem Solving Procedure

This study develops an integer programming-based model involving a multiobjective function. The multiobjective function was changed to a single-objective function by consideration of the transformer term which converts the utility related parameter to the cost related variable. The General Algebraic Modeling System (GAMS) software was used to solve the linear integer programming problem. GAMS is a high-level modeling system for mathematical programming and optimization. GAMS is specifically designed for modeling linear, nonlinear, and mixed integer optimization problems. The settings of the software were modified to meet the specifications of the solving algorithm of the work planning problem. The CPLEX solver and the MIP (Mix Integer Programming) option were set to solve the model since this study dealt with a linear integer programming problem. CPLEX is a GAMS solver that allows users to combine the high-level modeling capabilities of GAMS with the power of CPLEX optimizers. Notably, this software was used because it is widely used for accurate problem solving of linear models; however, any software with the ability to solve mathematical models can be employed for the purpose of this research.

3.3. Case Study

Since the runtime for problem solving is directly related to the length of the analysis period and the size of the pavement network, a small pavement network and short analysis period have been selected for evaluating the model and comparing the approaches. A 3-year data collection on pavement conditions and M&R implementation was carried out covering 10 pavement network sections in Mashhad, Iran. In order to be concise and since the aim of this article is to compare results of probabilistic model in condition of budget uncertainty with those of deterministic model, the budget specifications were brought here and values for the other parameters of the models were not described (S1).

In this case study, two budget modes were considered for each year. Table 8 gives the predicted budget for each year as well as the probability of it. If the predicted budget is not provided, 30% of relative deficit was estimated for each year’s budget. The budget values were predicted by inquiry from undertaken administration, but the related probabilities were estimated based on an obvious assumption; if the forecast is longer, the probability of occurrence of the event will be relatively reduced. As there are two budget allocations for each year of the planning period, subsequently, 8 scenarios exist for the planning period of 3 years (23 = 8). The scenarios are defined according to the whole situations can be possible for budgeting at each year. Accordingly, the budgets and their probabilities for each scenario could be determined as given in Table 9. Notably, the sum of probabilities is equal to 1.

4. Results Analysis

In the developed two-stage stochastic model, parameter ‘Ca’ is defined as the cost of inequality of the first stage variable and the second stage variable. Due to the decision change, the value of Ca is added to the costs. Ca is the cost that is paid by decision maker to mitigate the effect of budget variation on planning and total cost due to the scenario change. Thus, by minimizing the total cost, the decision made in the second stage has less variation compared to the first stage. On the other hand, the greater the value of Ca, the less the variation at the second stage in comparison with first stage and vice versa. Even so, if the value of Ca is greater than a boundary threshold, the increase in Ca has no influence on the decision and, subsequently, on the costs. The boundary threshold is the point where the decisions of the first and second stages are the same and that is why (based on the line plotted in Figure 1) the result of the model shows no sensitivity to the budget variation after that point (for the case study, Ca = 60000).

The value of Ca is dependent on the conditions of each project. Ca values of the case study are presented in Table 10 as a percentage of the total budget and the developed model was solved for each Ca value. Also, the total cost for each Ca was computed by finding the response value of the objective function (ob according to (1)) using the developed model.

Based on the results of Table 10, the total cost versus Ca is plotted in Figure 1. It can be observed that, with Ca increase, the total cost increases up until a point that the increase in Ca has no impact on total cost. As mentioned earlier, the reason of it is the same values for the first and second stage decision variables () from that point on. As shown in Figure 1, Ca of 60000 is approximately corresponding to the boundary threshold for the case study of this research. For the Ca less than the boundary threshold, first, there is smooth total cost decrease with the Ca reduction. Then, there is turning or break-even point on the line where the slope of the line increases dramatically with Ca decrease. This point is considered the best value for Ca because the values less than it lead to the high sensitivity of the results of the model to the variation of Ca. Also, as the Ca approaches zero, it can be inferred that the decision and plan alteration do not impose any additional cost which is contrary to the reality. Therefore, in the case study, 891 is the appropriate value for the Ca based on Figure 1.

Once Ca is determined, the effects of budget uncertainty on the results of M&R planning can be assessed with comparing the results obtained from the two-stage stochastic and the deterministic models. Notably, the deterministic model was developed similar to the probabilistic model except that the “s” index was eliminated from formulations. Consequently, some terms in the objective function and some restrictions were omitted or changed in the deterministic model. More precisely, in addition to using “x” as the decision variable instead of “y” in all of the equations, (30) and (45) were eliminated and objective function was changed as (46) in the deterministic model.Finally, parameter values corresponding to the first scenario were used to solve the deterministic model. According to (30), since the first stage decision variable is equal to the second stage decision variable in the first scenario (), the values obtained for first scenario could correspond to the first stage decision variable. Therefore, in order to make the comparison between the deterministic and stochastic models possible, the results of the first stage decision variable of the two-stage stochastic model were compared to the decision variable of deterministic model.

Moreover, in order to assess the impact of technical constraints on the results, the deterministic model was also solved without any technical constraints addressed by (31) to (44). As such, adding technical constraints to the deterministic model resulted in a 91% decrease in runtime for solving the problem. Therefore, considering the technical constraints in the model reduces the solution time of solving procedure. This advantage is applicable for huge networks or long-term planning durations and helps avoid expending too much time in solving M&R planning issues.

Table 11 presents the results of M&R actions assignment for the deterministic and stochastic models. Each number in Table 11 refers to M&R action category ID presented in Table 5. “Noting” shows that any M&R action was not assigned to the related section at the investigated year. It can be seen that 23 percent of the assignments in the stochastic model, i.e., 7 out of 30, are different from deterministic ones.

Based on the definitions given in Table 2, Figure 2 compares the component costs of the deterministic and stochastic models, while Figure 3 makes the comparison between the total costs of them. The total cost of the stochastic model is higher than the deterministic one due to the uncertainty. It can be inferred that considering the probability in budgeting variation within the modeling process and, subsequently, alteration in M&R plan, results in an increase in the costs at the network level during the planning period. Instead of this, the stochastic model provides several M&R plans that are compatible with budget variation.

Figure 4 illustrates the ratios of the effective costs to the total cost for both models. It can be seen that a significant portion of the total cost is due to the vehicle operation cost (Cu) indicating the higher sensitivity of it in both models. Thus, minimizing Cu has more impact on minimization of the total cost. As shown in Figures 2 and 3, although the M&R operation cost (Co) of deterministic model is higher than those of stochastic model, the lower cost of Cu of deterministic model leads to lessened total cost compared to the stochastic model.

Figure 5 shows each year’s consumed budget as well as the average consumed budget within the planning period for deterministic and stochastic models. In order to better distinguish the cost values in each year, Figure 6 gives each year’s effective costs separately for deterministic and stochastic models. It can be observed that each model consumed the most out of the available budget for first year. In other words, more budgets are used in the initial years of the planning period. This proves that the objective functions of the models have more sensitivity to the initial years of the planning period. The reason of it is the high influence of the vehicle operation cost (Cu) and, subsequently, the pavement condition on the total cost. Therefore, the model suggests that the maximum M&R operation costs (Co) in the initials years reach the desired pavement condition in order to reduce the Cu costs substantially. In the last year of the planning period, the pavement condition is approximately in a desired level and, consequently, the Cu does not increase rapidly as any dramatic deterioration exists. Therefore, in the last year of the planning period, M&R operation cost has greater effect on the results of the model compared to the Cu and the model intends to reduce M&R operation cost (Co). Subsequently, the total cost reduces gradually in the last year of the planning period.

Table 12 yields the total costs of each scenario for deterministic and stochastic model. Also, the percentage increases in the total costs in accordance with the deterministic model are given. It can be observed that the total costs of the stochastic model are higher than those of the deterministic one. For scenario 1, it is 3%, and in the worst scenario that there is deficit in all of years of the planning period, it is obtained as 42%. According to Table 9, the maximum budget deficit in the stochastic model can be 2674 (difference between 8914 and 6240) thousand USD by comparing the best and worst case scenarios (the and scenarios). In accordance with Table 12, this amount resulted in an increase of 7374 (difference between 26807 and 19433) thousand USD in the total cost which is approximately three times greater than the maximum budget deficit. This indicates the importance of funding in pavement management.

Although using the two-stage stochastic model increases the total cost, it has some advantages. Table 13 presents the results of M&R assignments in each planning year for each scenario. Since the different conditions of budgeting are considered in the planning process, a flexible M&R planning is resulted for the budget variations.

Due to the Nonanticipative constraint considered in scenarios, it can be seen that in all scenarios with equal budgets at prior years, similar M&R actions were assigned. Such M&R assignment gives a plan that can be pursued at any year of the planning period according to budgets which had been implemented in prior years.

5. Summary and Conclusion

An M&R planning model which considers nondeterministic budget was developed in this study based on a linear integer two-stage stochastic programming approach. The model can be used for large-scale networks or for long-term planning periods, although it addresses a great number of details. As such, the effects of budget uncertainty on results of M&R planning were assessed by comparing the results obtained from the two-stage stochastic and the deterministic models. From the comparison of the results of the deterministic and stochastic models, following conclusions can be drawn:(i)Taking the probabilistic effect of the budgeting within modeling process and, subsequently, alteration in M&R planning into account results in an increase in the costs of the planning period. It was computed as 3% for the case study. However, considering the various conditions for budgeting in the two-stage stochastic model leads to flexible M&R actions planning for budget variation so that it can be followed at any year of the planning period regardless of having the exact expected budget in that year.(ii)The annual budget deficit leads to the increased costs during the planning period. The total financial loss could reach 3 times as much as the total deficit for the case study.(iii)Both deterministic and stochastic models have the most sensitivity to the vehicle operation cost (Cu). Thus, minimization of the Cu has the greatest impact on the minimization of the total cost.(iv)Both deterministic and stochastic models used the most of the available budget in the first year of the planning period. In other words, the total cost has greater sensitivity to the initial years of the planning period, so providing sufficient budget for the initial years is of much importance.

This study can be helpful to other researchers in area relevant to pavement M&R planning since it employs a linear integer programming model. There is not any limitation for application of the model if the parameters values are fully available. However, completely gathering data for a huge network or for a long time planning period could impose considerable costs to undertaking administrations. It would be worthwhile for other researchers and scholars to develop this model by using heuristic and metaheuristic methods like genetic algorithm (GA) and others which also reduce runtime in solving the problem. In addition to budgeting, this study can be considered by other researchers as it takes into account uncertainty in some other parameters such as pavement performance deterioration, traffic estimations, and other similar items. Moreover, determining pavement condition effects on accidents rate and related crash costs and, consequently, considering the crash cost effects in the modeling can be taken into account by other researches. Also, a comparison of the model, not only with the deterministic approach but also with other approaches such as robust optimization or fuzzy mathematical programming that allow for inclusion of the uncertainty of budgeting, could be greatly noteworthy for researchers.

Data Availability

The GAMS software input data used in this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Supplementary Materials

An example of coding in the GAMS software that contains all the input data about the model parameters. (Supplementary Materials)