# A Two-Stage Model for Project Optimization in Transportation Infrastructure Management System

Authors: Zhang Chen; Liyuan Liu; Li Li; Hui Li
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

Mathematical optimization is very important for project decision in the Transportation Infrastructure Management System (TIMS). However, it has not been widely employed in TIMS due to poor performance of conventional optimization models in calculation speed and practical application. Therefore, it is necessary to improve the performance of optimization models. According to the process of decision-making in transportation management, a novel two-stage project optimization model, including budget allocation and project distribution, was proposed in this paper. Moreover, the methods of dynamic programming (DP) and genetic algorithm (GA) were applied to obtain an effective solution. The findings indicate that the new optimization method can provide a satisfactory and reasonable maintenance schedule for transportation infrastructure maintenance agencies whose routine management will benefit from the newly proposed model.


---

## Body

## 1. Introduction

Transportation infrastructure includes roads, bridges, tunnels, airports, railways, and seaports. As a result, Transportation Infrastructure Management System (TIMS) correspondingly covers many subsystems, among which Pavement Management System (PMS) and Bridge Management System (BMS) are the most important ones. Whatever kind of subsystem it is, project optimization is a key element in the process of decision-making for the infrastructure management. Specifically, project optimization [1] refers to finding an optimal maintenance strategy with maximized benefit through arranging pavement maintenance reasonably in terms of time and space in the planning period. Most of existing research efforts on project optimization of TIMS focused on PMS.Since decision supporting system was introduced into the second generation of PMS, project optimization has been paid much attention by pavement researchers and management agencies. Currently, there are mainly two categories of project optimization methods for network-level pavement management system, namely, prioritization method and mathematical optimization [2, 3]. Prioritization method is to carry out project selection based on some principles prescreened and then to determine the maintenance strategy for each year in the planning period. In 1994, Hass et al. [4] summarized the characteristics of different prioritization methods, among which the two most popular methods are based on infrastructure performance parameters and economic analysis parameters. The prioritization method based on performance parameters is adopted to arrange road maintenance projects by California, US, where road roughness, damaged condition, and average daily traffic volume are considered as influencing factors to develop the prioritization principles. Hudson uses rainfall, number of freeze-thaw cycles, and damage severity as parameters to establish an expression for prioritization through regression or variance analysis so as to determine the priority index of every section requiring reconstruction in the road network. While prioritization based on performance parameters is convenient for calculation and of high pertinence, the results may be far from the economic optimum. Thus the prioritization method based on economic analysis parameters is relatively advanced regarding economic optimum. In 1980s, the prioritization method based on economic analysis parameters was very popular [5–8]. The State of Washington, US, used the total cost, including initial project construction cost, maintenance cost, user operating expense, delay cost, and pavement salvage, as the prioritization indicator in their PMS while the performance index was adopted in UK. At the same time, Hass et al. used equivalent annual cost in Canada, which refers to the ratio of the project construction cost and the corresponding expected lifetime to sequence. Among the PMSs adopting the prioritization method, the most typical one is PAVER in US. PAVER uses the method of benefit-cost-increment to determine the priorities of maintenance projects according to budget optimization. Another representative one is the Pavement Maintenance Decision Support System of Shanghai (PMDSSS), China. According to different managerial preferences of road administrators and agencies, this PMDSSS has developed eight different prioritization principles taking into account either performance parameters or economic analysis parameters or both, with different combinations of traffic volume, pavement damage condition, riding quality, structural capability, and economic indicators. The principles greatly help policy makers to formulate large or medium maintenance plans and long-term rehabilitation strategies in the planning years. However, the tradeoff between the maintenance strategy and the time has not been included in the prioritization method, which leads to huge disparity between the calculated result and the actual optimal solution. Therefore, many researchers have been, instead, focusing on the mathematical optimization, which refers to considering each project in the planning period, the possible maintenance plan, and the implementation schedule through mathematical calculations.The idea of optimization using mathematical method was proposed in 1970. The mathematical optimization for PMS falls into two categories: static and dynamic. Integer programming is mostly used in the former. For instance, PMS in Denmark [9], HDM-III in World Bank [10], RAMS of Texas, US, and PARS of Ontario, Canada, all adopt 0-1 integer programming [11, 12]. PMS in the State of Indiana, US [13, 14], uses the number of workday needed for maintenance as the decision variable, and constraints from budget, manpower, machine, and material are also included to build up the integral programming model. In the dynamic programming, PMS of Arizona, US [15], was the first to successfully introduce Markov decision process to the network-level pavement management system. When improving the PAVER system, Feighan et al. [16] managed to develop the strategy for maintenance and reconstruction that minimizes the cost of road network through dynamic programming. In 1994, Liu and Yao [17] used Markov decision process to minimize maintenance cost. In 1995, Zou [18] formulated a dynamic model to predict the pavement performance and the adaptive iteration algorithm. In addition, to solve the problems of pavement maintenance strategy, he adopted the analytical hierarchy process model and the heuristic optimization technique. In 2001, Nunoo [19] optimized the maintenance plan of integrated pavement using shuffled complex evolution algorithm. In 2003, Chan et al. [20, 21] attempted to use genetic algorithm to develop a multiobjective programming model to optimize the strategy for highway pavement maintenance. Meanwhile, they adopted heuristic algorithm to yield the result. Ferreira et al. [22] also developed a probabilistic segment-linked optimization model together with a genetic algorithm heuristic with the objective of minimizing total discounted cost of M&R actions. While mathematical optimization can produce the optimal calculation results, the large number of factors to be considered, huge data processing, and barely satisfactory calculation speed greatly limit its practical application. In order to apply mathematic optimization to project optimization, the key point is to find a model or algorithm which can save computing resources as well as meeting the practical needs.Compared to the existing research studies, the main contributions of this paper are listed as follows. Firstly, the project optimization is separated into two independent and interrelated processes, namely, budgets allocation and project distribution, and a new two-stage model is developed. Secondly, the technique of dynamic programming (DP) and genetic algorithm (GA) are applied to solve the model and yield an effective solution. Finally, the new programming method is verified to be effective through the case study in Shanghai and the poor calculation speed and the practical application limitations of conventional methods are improved.


## 2. Integer Programming in Conventional Project Optimization

Integer programming (IP) in the conventional project optimization model can be stated as follows:(1)Objectivemax⁡⁡Z=∑t=1T∑i=1N∑j=1mXijt*Bijt(2)Subjecttozz∑t=1T∑i=1N∑j=1mXijt*Cijt≤A(3)∑j=1mXijt=1(i=1,2,…,N;t=1,2,…,T)(4)Xijt={1iftreatmentjisappliedinsegmentiinyeart0otherwise(5)Bijt,Cijt,A,St,Vt>0,
where Z is total maintenance benefit; Bijt, Cijt are, respectively, the benefit and cost caused by implementing treatment j in project i in year t; A is the total budget in planning period; T is the length of planning period, normally 5 or 10 years; N is the amount of total road units; m is the total number of treatments for each project.The model above is a large integer programming problem in which the set of feasible solutions is very huge. Therefore, it is necessary to simplify the model in order to calculate conveniently as well as meeting practical needs. Equation (1) is transformed to the following form:
(6)∑t=1T∑i=1N∑j=1mXijt*Bijt=∑i=1N∑j=1mXij1*Bij1+⋯+∑i=1N∑j=1mXijk*Bijk+⋯+∑i=1N∑j=1mXijT*BijT.
If the expression ft(Yt)=∑i=1N∑j=1mXijt*Bijt (Yt is the budget of year t) is considered as the maintenance benefit of year t, then the total maintenance benefit in planning years should be equal to the summation of maintenance benefit of each year. Therefore,
(7)∑t=1T∑i=1N∑j=1mXijt*Bijt=f1(Y1)+⋯+ft(Yt)+⋯+fT(YT)=∑t=1Tft(Yt),
where Yt is budget of year t and ft(Yt) is maintenance benefit of year t.It is shown in (1) and (7) that solution of project optimization can be approached with the following steps. The first step is budget optimizing to get the budget allocation in each year in the planning period. The second step is project portfolio (including optimization of project scheme and schedule) to get the optimal benefit maintenance strategy in each year based on the budget allocation of first step. After several iterations of optimization, selection, and comparison through the repetition of above two steps, the optimal maintenance strategy under the optimal budget allocation can be achieved finally. Actually, this method reflects well the actual iterative decision-making process of the government agencies, which is “budget allocation, project arrangement, budget adjustment, and project adjustment.” Through this iterative method, project optimization is divided into two relatively simple processes, namely, budget optimization and project distribution.


## 3. Two-Stage Optimization Approach

According to the above analysis, project optimization can be divided into two stages. The first one is about how to allocate budgets. The second one determines maintenance projects in each year based on the budget allocation in the first stage. The first stage is defined as budget allocation model and the second as project distribution model, which build up the two-stage optimization approach.


### 3.1. Budget Allocation Model

Budget allocation model is, under certain amount of fund, to find a reasonable budget allocation for each year in planning period in order to maximize the benefit of budget. The model focuses on the reasonable way of budget allocation.In the process of allocating budgets, pavement management agencies need to consider not only the performance of road network but also the budget limits (maximal and minimal budgets for each year are expressed withSt and Vt, resp.). Therefore, the budget allocation model can be designed as follows:
(8)Objectivemax⁡∑t=1Tft(Yt)=f1(Y1)+f2(Y2)ssssssssddssssssssssssssssssss+⋯+ft(Yt)+⋯+fT(YT)(9)Subjectto∑t=1TYt≤A(10)Yt≤St(t=1,2,…,T)(11)Yt≥Vt(t=1,2,…,T).
In fact, there are a large number of maintenance projects in the road network for the planning period, and each project has several treatments. Different treatments bring different maintenance benefits. Therefore, it is necessary to pose some limitations to budget allocation considering that pavement management agencies may utilize fuzzy decisions unconsciously in the process of decision-making which makes their decisions rather reasonable. Firstly, budget change is discontinuous when management agencies are allocating or adjusting budgets; therefore the budget change is set as integer times of the smallest unit. Secondly, maximal maintenance benefit in the planning period under certain budget allocation schedule equals the summation of maximal benefits in each year:
(12)∑t=1T‍ft(Yt)=max⁡(f1(Y1))+⋯+max⁡(ft(Yt))+⋯+max⁡(fT(YT))=∑t=1T‍max⁡(ft(Yt)).


### 3.2. Project Distribution

Project distribution model is about how to arrange maintenance projects under given budget allocation schedule(Y1,Y2,…,YT) in each year in order to obtain the maximal maintenance benefit. It can be designed as follows:
(13)Objective∑t=1T(max⁡∑i=1n∑j=1m(XijtBijt))(14)Subjectto∑i=1n∑j=1m(XijtCijt)≦Yt,(t=1,2,…,T)(15)∑j=1mXijt≤1(i=1,2,…,n;t=1,2,…,T),Xijt=0or1,
where n is the total amount of projects in road network; m is the total amount of treatments for each project; Boolean variable Xijt=1 (or 0) means treatment j is applied (or not) to project i in year t.


### 3.3. Relation of Budget Allocation Model and Project Distribution Model

The results of the budget allocation model and the project distribution model are interactional. At first, annual budget will be provided through the budget allocation model. Then, the maximal maintenance benefit and the project schedule will be constructed through the project distribution model according to the result of the first model and then a feedback will be offered to the first model to judge whether the solution is the best or not. If not, the budget allocation model will optimize again and produce a new array of budgets in each year, then the corresponding total maintenance benefit and maintenance schedule will be obtained through the project distribution model. Through iterations of the two models, the optimal strategy will be obtained in the planning years.


## 3.1. Budget Allocation Model

Budget allocation model is, under certain amount of fund, to find a reasonable budget allocation for each year in planning period in order to maximize the benefit of budget. The model focuses on the reasonable way of budget allocation.In the process of allocating budgets, pavement management agencies need to consider not only the performance of road network but also the budget limits (maximal and minimal budgets for each year are expressed withSt and Vt, resp.). Therefore, the budget allocation model can be designed as follows:
(8)Objectivemax⁡∑t=1Tft(Yt)=f1(Y1)+f2(Y2)ssssssssddssssssssssssssssssss+⋯+ft(Yt)+⋯+fT(YT)(9)Subjectto∑t=1TYt≤A(10)Yt≤St(t=1,2,…,T)(11)Yt≥Vt(t=1,2,…,T).
In fact, there are a large number of maintenance projects in the road network for the planning period, and each project has several treatments. Different treatments bring different maintenance benefits. Therefore, it is necessary to pose some limitations to budget allocation considering that pavement management agencies may utilize fuzzy decisions unconsciously in the process of decision-making which makes their decisions rather reasonable. Firstly, budget change is discontinuous when management agencies are allocating or adjusting budgets; therefore the budget change is set as integer times of the smallest unit. Secondly, maximal maintenance benefit in the planning period under certain budget allocation schedule equals the summation of maximal benefits in each year:
(12)∑t=1T‍ft(Yt)=max⁡(f1(Y1))+⋯+max⁡(ft(Yt))+⋯+max⁡(fT(YT))=∑t=1T‍max⁡(ft(Yt)).


## 3.2. Project Distribution

Project distribution model is about how to arrange maintenance projects under given budget allocation schedule(Y1,Y2,…,YT) in each year in order to obtain the maximal maintenance benefit. It can be designed as follows:
(13)Objective∑t=1T(max⁡∑i=1n∑j=1m(XijtBijt))(14)Subjectto∑i=1n∑j=1m(XijtCijt)≦Yt,(t=1,2,…,T)(15)∑j=1mXijt≤1(i=1,2,…,n;t=1,2,…,T),Xijt=0or1,
where n is the total amount of projects in road network; m is the total amount of treatments for each project; Boolean variable Xijt=1 (or 0) means treatment j is applied (or not) to project i in year t.


## 3.3. Relation of Budget Allocation Model and Project Distribution Model

The results of the budget allocation model and the project distribution model are interactional. At first, annual budget will be provided through the budget allocation model. Then, the maximal maintenance benefit and the project schedule will be constructed through the project distribution model according to the result of the first model and then a feedback will be offered to the first model to judge whether the solution is the best or not. If not, the budget allocation model will optimize again and produce a new array of budgets in each year, then the corresponding total maintenance benefit and maintenance schedule will be obtained through the project distribution model. Through iterations of the two models, the optimal strategy will be obtained in the planning years.


## 4. Solving the Two-Stage Model

The key to applying the mathematic optimization method into practice is to find a reasonable solving method which ensures certain accuracy and meets practical requirements. In the objective function of (8), it is difficult to be expressed for fi(Yi) by explicit functions, which causes difficulty in effectively solving it with normal algorithm. Given that genetic algorithm (GA) is able to search the global optimal solution in a complicated space while its objective functions are not necessary to be explicit functions, GA is chosen to solve the budget allocation model (8).In the budget allocation model, decision (maintenance strategy) of each stage (planning year) is the function of pavement network condition given by the decision of last stage. That is, each stage will influence the next stage through the network condition resulting from the different treatment strategy. Therefore, under a given budget allocation(Y1,Y2,…,YT), the formulation of pavement network maintenance strategy in planning years is a multistage problem in which all stages have mutual connections, namely, a classic dynamic programming problem. Therefore, the method of dynamic programming is employed to solve the projects distribution model (13).


### 4.1. GA for Budget Allocation Model

Before the calculation with GA, some parameters including length of chromosome, species group scale, intersection rate, and mutation rate have to be set in advance according to the scale of pavement network, length of planning period, and budget. Then, the budget allocation model can be solved in the following process, as shown in Figure1.Figure 1
Calculation process of budget allocation with GA.(1)
Initial species: each chromosome in the species represents a budget allocation mode.
(2)
Fitness function: the fitness function of chromosome is derived from the maximal maintenance benefit under a budget allocation mode. The function is computed using the project allocation model.
(3)
Genetic manipulation of chromosome: through the genetic manipulation of chromosome, the budget allocation with superior maintenance benefit can be inherited and those inferior ones will be eliminated.
(4)
Termination of iteration: certain iteration times can be set as the termination condition of GA.In the above calculation process, the fitness function of chromosomes, which are not subject to the budget constraints, can be limited by penalty function or be given a very small fitness value. The chromosome with the biggest fitness value gets to be obtained by GA and the corresponding budget allocation manner is the final solution of the budget allocation model.


### 4.2. DP for Project Distribution Model

Project distribution model is mainly used to arrange project implementation in each year under a given optimization result of budget allocation. The output results should include a list of maintenance projects in each planning year. Meanwhile, through the project distribution model, the fitness value can be achieved for the above-mentioned chromosome which can serve as the basis of calculation for further optimization.As discussed above, the determination of maintenance strategy of pavement network in each stage is a multistage decision-making problem in which each stage is interactional. The process has the following characters: at the beginning of each stage, optimal decision of each stage is only related to road network condition, but not related to the decision of the previous stage in which the road network condition is known. The combination of optimal decision in each stage is the best strategy for the planning period. On the other hand, the process of decision-making is a multistage chain in which each stage does not influence the following stages. Therefore, the dynamic programming is a powerful tool to solve this kind of problem.


#### 4.2.1. Calculation Process of DP

According to the characteristics of DP and project distribution model, the process of DP in project distribution can be described as shown in Figure2.Figure 2
DP process of project distribution model.In Figure2, the whole planning period of T years can be divided into T stages. Vector Yt (pavement network condition in year t, budget of year t) represents initial condition of each stage. The decision can be made according to Ut(Yt)=Ut (pavement network condition in year t, budget of year t). Finally the maximal maintenance benefit, benefit in this stage, can be calculated based on the Ut(Yt) and the condition transition function Yt+1=F(Yt,Ut(Yt)).


#### 4.2.2. The Calculation of Maximal Maintenance Benefit (Determination ofUt(Yt))

Calculation model of maximal maintenance benefitUt(Yt) in each year can be devised as follows:
(16)max⁡∑i=1n∑j=1m(XijBij)(17)s.t∑i=1n∑j=1m(XijCij)≤y(18)∑j=1mXij≤1,Xij=0or1,
where n is the number of projects; m is the number of treatments for each project; y is the total maintenance budget; Xij=1 (or 0) means treatment j is selected (or not) in project i.The following is a recursion equation built up by the DP method, in which there are 3 treatments for each project. The cost of each treatment is defined as integer and described from small to large aswi1, wi2, and wi3. At the same time, the corresponding benefit is vi1, vi2, and vi3, respectively:
(19)f(i,m)={f(i+1,m)0<m<wi1max⁡⁡(f(i+1,m),iiif(i+1,m-wi1)+vi1)wi1<m<wi2max⁡⁡(f(i+1,m),f(i+1,m-wi1)+vi1,f(i+1,m-wi2)+vi2)wi2<m<wi3max⁡⁡(f(i+1,m),f(i+1,m-wi1)+vi1,f(i+1,m-wi2)+vi2,iiif(i+1,m-wi3)+vi3)Xi>m>wi3.
Then constraints can be formulated as follows:
(20)f(n,m)={0m<wn1vn1wn2<m<wn1max⁡(vn1,vn2)wn2<m<wn3max⁡(vn1,vn2,vn3)m>wn3.
The function of f(i,m) can be solved through the backward deduction method and the optimal solution for each stage can be achieved through the backtracking method.


## 4.1. GA for Budget Allocation Model

Before the calculation with GA, some parameters including length of chromosome, species group scale, intersection rate, and mutation rate have to be set in advance according to the scale of pavement network, length of planning period, and budget. Then, the budget allocation model can be solved in the following process, as shown in Figure1.Figure 1
Calculation process of budget allocation with GA.(1)
Initial species: each chromosome in the species represents a budget allocation mode.
(2)
Fitness function: the fitness function of chromosome is derived from the maximal maintenance benefit under a budget allocation mode. The function is computed using the project allocation model.
(3)
Genetic manipulation of chromosome: through the genetic manipulation of chromosome, the budget allocation with superior maintenance benefit can be inherited and those inferior ones will be eliminated.
(4)
Termination of iteration: certain iteration times can be set as the termination condition of GA.In the above calculation process, the fitness function of chromosomes, which are not subject to the budget constraints, can be limited by penalty function or be given a very small fitness value. The chromosome with the biggest fitness value gets to be obtained by GA and the corresponding budget allocation manner is the final solution of the budget allocation model.


## 4.2. DP for Project Distribution Model

Project distribution model is mainly used to arrange project implementation in each year under a given optimization result of budget allocation. The output results should include a list of maintenance projects in each planning year. Meanwhile, through the project distribution model, the fitness value can be achieved for the above-mentioned chromosome which can serve as the basis of calculation for further optimization.As discussed above, the determination of maintenance strategy of pavement network in each stage is a multistage decision-making problem in which each stage is interactional. The process has the following characters: at the beginning of each stage, optimal decision of each stage is only related to road network condition, but not related to the decision of the previous stage in which the road network condition is known. The combination of optimal decision in each stage is the best strategy for the planning period. On the other hand, the process of decision-making is a multistage chain in which each stage does not influence the following stages. Therefore, the dynamic programming is a powerful tool to solve this kind of problem.


### 4.2.1. Calculation Process of DP

According to the characteristics of DP and project distribution model, the process of DP in project distribution can be described as shown in Figure2.Figure 2
DP process of project distribution model.In Figure2, the whole planning period of T years can be divided into T stages. Vector Yt (pavement network condition in year t, budget of year t) represents initial condition of each stage. The decision can be made according to Ut(Yt)=Ut (pavement network condition in year t, budget of year t). Finally the maximal maintenance benefit, benefit in this stage, can be calculated based on the Ut(Yt) and the condition transition function Yt+1=F(Yt,Ut(Yt)).


### 4.2.2. The Calculation of Maximal Maintenance Benefit (Determination ofUt(Yt))

Calculation model of maximal maintenance benefitUt(Yt) in each year can be devised as follows:
(16)max⁡∑i=1n∑j=1m(XijBij)(17)s.t∑i=1n∑j=1m(XijCij)≤y(18)∑j=1mXij≤1,Xij=0or1,
where n is the number of projects; m is the number of treatments for each project; y is the total maintenance budget; Xij=1 (or 0) means treatment j is selected (or not) in project i.The following is a recursion equation built up by the DP method, in which there are 3 treatments for each project. The cost of each treatment is defined as integer and described from small to large aswi1, wi2, and wi3. At the same time, the corresponding benefit is vi1, vi2, and vi3, respectively:
(19)f(i,m)={f(i+1,m)0<m<wi1max⁡⁡(f(i+1,m),iiif(i+1,m-wi1)+vi1)wi1<m<wi2max⁡⁡(f(i+1,m),f(i+1,m-wi1)+vi1,f(i+1,m-wi2)+vi2)wi2<m<wi3max⁡⁡(f(i+1,m),f(i+1,m-wi1)+vi1,f(i+1,m-wi2)+vi2,iiif(i+1,m-wi3)+vi3)Xi>m>wi3.
Then constraints can be formulated as follows:
(20)f(n,m)={0m<wn1vn1wn2<m<wn1max⁡(vn1,vn2)wn2<m<wn3max⁡(vn1,vn2,vn3)m>wn3.
The function of f(i,m) can be solved through the backward deduction method and the optimal solution for each stage can be achieved through the backtracking method.


## 4.2.1. Calculation Process of DP

According to the characteristics of DP and project distribution model, the process of DP in project distribution can be described as shown in Figure2.Figure 2
DP process of project distribution model.In Figure2, the whole planning period of T years can be divided into T stages. Vector Yt (pavement network condition in year t, budget of year t) represents initial condition of each stage. The decision can be made according to Ut(Yt)=Ut (pavement network condition in year t, budget of year t). Finally the maximal maintenance benefit, benefit in this stage, can be calculated based on the Ut(Yt) and the condition transition function Yt+1=F(Yt,Ut(Yt)).


## 4.2.2. The Calculation of Maximal Maintenance Benefit (Determination ofUt(Yt))

Calculation model of maximal maintenance benefitUt(Yt) in each year can be devised as follows:
(16)max⁡∑i=1n∑j=1m(XijBij)(17)s.t∑i=1n∑j=1m(XijCij)≤y(18)∑j=1mXij≤1,Xij=0or1,
where n is the number of projects; m is the number of treatments for each project; y is the total maintenance budget; Xij=1 (or 0) means treatment j is selected (or not) in project i.The following is a recursion equation built up by the DP method, in which there are 3 treatments for each project. The cost of each treatment is defined as integer and described from small to large aswi1, wi2, and wi3. At the same time, the corresponding benefit is vi1, vi2, and vi3, respectively:
(19)f(i,m)={f(i+1,m)0<m<wi1max⁡⁡(f(i+1,m),iiif(i+1,m-wi1)+vi1)wi1<m<wi2max⁡⁡(f(i+1,m),f(i+1,m-wi1)+vi1,f(i+1,m-wi2)+vi2)wi2<m<wi3max⁡⁡(f(i+1,m),f(i+1,m-wi1)+vi1,f(i+1,m-wi2)+vi2,iiif(i+1,m-wi3)+vi3)Xi>m>wi3.
Then constraints can be formulated as follows:
(20)f(n,m)={0m<wn1vn1wn2<m<wn1max⁡(vn1,vn2)wn2<m<wn3max⁡(vn1,vn2,vn3)m>wn3.
The function of f(i,m) can be solved through the backward deduction method and the optimal solution for each stage can be achieved through the backtracking method.


## 5. Model Verification through Case Study

The model and algorithm are tested and verified through a case study with the data collected from asphalt pavements of 12 districts in Shanghai in 2004. The planning period lasts 10 years, the total number of sections surveyed is 867, the total road length is 238,332 meters, and the total area is 3,093,293 square meters.(1)
Parameters of funds are provided as follows:
(a)
total budget: 160 million RMB,(b)
minimum investment for each year: 13 million RMB,(c)
maximum investment for each year: 17 million RMB,(d)
the step change of investment for each year: 100 thousand RMB.(2)
Parameters of GA are provided as follows:
(a)
reproduction rate: 0.7,(b)
mutation rate: 0.05,(c)
crossover rate: 0.3.The calculation software is programmed with Microsoft Visual Basic 6.0. The model proposed and developed in this paper (hereafter referred to as “new method”) is compared with the current pavement maintenance method (hereafter referred to as “conventional method”), in which maintenance would be conducted once the PCI is lower than 75, currently adopted by pavement management agencies in Shanghai. The calculation results are listed in Tables 1 and 2.Table 1
Results of new method.
Year
No. of segments to be treated
Length of road to be treated (m)
Area of pavement to be treated (m2)
Total maintenance cost (¥10,000)
Average PCI of pavement network
Average PII of pavement network
2005
60
17859
187210
1320
89.6
483649.2
2006
56
17234
223318
1470
88.3
361548
2007
63
17298
221853
1540
87.1
571428.7
2008
81
23026
303050
1450
86.8 (minimal)
888703.5
2009
89
23103
307059
1700
87
782662.4
2010
112
28812
334866
1670
88.3
780185.2
2011
85
22393
310335
1670
88.8
755266.7
2012
80
22567
328590
1660
89.3
675452.1
2013
98
22582
272024
1700
90.6
520031.3
2014
75
22688
319542
1580
91.1 (minimal)
504012.8
Total
799
217562
2807847
15760
88.69
6322940Table 2
Results of conventional method.
Year
No. of segments to be treated
Length of road to be treated (m)
Area of pavement to be treated (m2)
Total maintenance cost (¥10,000)
Average PCI of pavement network
Average PII of pavement network
2005
78
24152
269404
2440
90.27
572719
2006
43
13433
173309
1188
88.58 (minimal)
302004
2007
73
20923
257218
2145
87.88
621756
2008
133
37972
461683
2820 (maximal)
89.26
1131875
2009
121
33224
446845
2810
90.77
962347
2010
118
29992
366621
2356
92.51
749037
2011
115
31383
435860
2559
94.38
864322
2012
79
19932
279639
1370
95.09
551289
2013
63
15851
224198
1039
95.13 (maximal)
340663
2014
34
8620
130474
559 (minimal)
94
150450
Total
857
235482
3045251
19286
91.78 (average)
6246469In the tables, the maintenance benefit of the pavement network is indicated with PII (Pavement Improvement Index). PII mainly considers restoration of pavement indicators, interval time between two maintenance services, and social impact. The comparison of the data in Tables1 and 2 is shown in Figures 3 and 4.Figure 3
Maintenance cost comparison of road network.Figure 4
Average PCI comparison of road network.Based on the case study results, the following observations can be drawn.(1)
The total maintenance benefit increases (the PII improves from 6246469 to 6322940 after optimization with the new method), while the average PCI of road network in the planning period decreases after adopting the new method. At the same time, the total cost in the planning period is significantly reduced (the total cost decreases from¥19286 to ¥15760, a drop of 18%). That is, pavement management agencies will achieve better maintenance benefit through spending fewer resources with the new method.(2)
As shown in Figure3, through the conventional method, budget demand in each planning year is significantly different (i.e., the maximal budget demand is ¥28,200,000 in the fourth year, and the minimal demand is ¥5,590,000 in the tenth year. The variance of the investment in each year is up to 2447). That seriously violates the rule of a stable financial plan of local finance department. Through the new method, the budget of each planning year is well-controlled. The difference of each year’s budget effectively decreases (i.e. the maximal budget demand is ¥17,000,000 in the fifth and ninth year, and the minimal demand is ¥13,200,000 in the first year. The variance of the investment in each year is lowered down to 387). It well reflects the rule of a stable financial plan of finance department.(3)
Figures4 and 5 indicate that the PCI and the PII of road network in each year, similar to the maintenance budget, change remarkably when adopting the conventional method. On the contrary, the new method can significantly mitigate the variation of PCI and PII. The variance of PCI decreases from 8.3 to 4.5 and the variance of PII decreases from 928403 to 505346.In general, pavement management agencies will get more maintenance benefit with fewer funds through the new model. Meanwhile, the amount of investment and the variability of pavement network PCI in each planning year can be effectively controlled to meet the rule of a stable project plan for pavement management agencies. Consequently, the new model and the solving method proposed in this paper are effective and practical for the infrastructure management.Figure 5
PII comparison.


## 6. Conclusion

Choosing a proper model and algorithm is very critical for methods of mathematic optimization. To solve the problems of poor calculation speed and practical application limitation of the conventional mathematic method, a two-stage model consisting of the budget allocation model and the project distribution model is developed in this paper. In addition, the model is solved through dynamic programming and genetic algorithm. Furthermore, it is verified to be effective by practical data in a case study. The findings indicate that the new optimization method can provide a satisfactory and reasonable maintenance schedule for transportation infrastructure maintenance agencies whose routine management will benefit from the newly proposed model.



---

*Source: 914515-2014-07-13.xml*