Abstract

The objective of the present study is to develop models to predict the deterioration of pavement distress of the urban road network. Genetic programming (GP) has been used to develop five models for the prediction of pavement distress: Model 1 for the cracking progression, Model 2 for the ravelling progression, Model 3 for the pothole progression, Model 4 for the rutting progression, and Model 5 for the roughness progression. The data have been collected from the roads of Patiala City, Punjab, India; during the years 2012–2015, the network of 16 roads have been selected for the data collection purposes. The data have been divided into two sets, that is, training dataset (data collected during the years 2012 and 2013) and validation dataset (data collected during the years 2014 and 2015). The two fitness functions have been used for the evaluation of the models, that is, coefficient of determination (R2) and root mean square error (RMSE), and it is inferred that GP models predict with high accuracy for pavement distress and help the decision makers for adequate and timely fund allocations for preservation of the urban road network.

1. Introduction

India’s road network is the second largest in the world, with the total road length increased from 0.399 million kilometers as on 31 March 1951 to 4.698 million kilometers as on 31 March 2011 [1]. The urban roads clutch huge traffic and bear the economic loss due to poor condition of the roads that accumulates to a gigantic sum. The total transportation cost comprised two basic components, that is, road cost and user cost. Researcher’s studies [2, 3] on various aspects of road pavement and transportation have indicated that the road user cost constitutes a major part (around 80–90% of the total transportation cost) which could be saved substantially with the help of a good decision support system. The utilization of urban roads and problems faced while maintenance are totally different than highway road network. In urban roads, the trip length is small, the number of trips and diversion is large, and the intensity of traffic is heavy and nonuniform. In the last twenty years, the maintenance management systems have improved significantly due to the advances in technologies. Many researchers are working on the development of pavement maintenance management system (PMMS) and the effect of road surface deformations [417]. PMMS is one of such tools that assist engineers to maintain the roads in their serviceable conditions [17]. Mathew et al. [18] made efforts to develop the deterioration models for performance of rural roads using artificial neural network (ANN) and compared with regression techniques. They concluded that an ANN model proved to be more suitable than the conventional regression models because of its ability to adjust to a changing environment. Fwa et al. [19] presented a genetic algorithm- (GA-) based procedure for solving multiobjective network level pavement maintenance programming problems. They adopted the concepts of Pareto optimal solution set and rank-based fitness evaluation and two methods for selecting an optimal solution. Chan et al. [20] applied GA for optimizing the resource allocation for pavement maintenance programming. They also compared the performance of the various constraint handling methods, namely, prioritized resource allocation method, decode-repair method, and penalty method. Terzi [21] used ANN in modeling the Pavement Service Rating (PSR) for the flexible pavements. They developed an ANN model using input variables which were longitudinal cracking, total cracking, patching, rut depth, slope variance, and output as panel data of PSR. The model was trained and tested using 74 datasets obtained from AASHO test results. The comparison of PSR had been done with the ANN model and PCI had been obtained from AASHO equation. The values of coefficient of determination (R2) for the ANN model were 0.83 and 0.82 for the training set and testing set, respectively, and for PSI, the values of R2 were 0.99 and 0.87 for the training and testing set, respectively. Reddy and Veeraragavan [22] developed deterioration models and growth curves for rebound deflection, rut depth, and crack area for two types of pavement surfacing. They developed predictive models for deflection, rut depth, and crack area progression with traffic. Bianchini and Bandini [23] proposed a neurofuzzy model to predict the performance of flexible pavements using the parameters routinely collected by agencies to characterize the condition of an existing pavement. Their proposed hybrid model for predicting pavement performance was characterized by multilayer feed-forward neural networks that led the reasoning process of the IF-THEN fuzzy rules. Khraibani et al. [24] presented a nonlinear mixed-effects model for the evaluation and prediction of pavement deterioration. Luo [25] proposed the application of an auto regression method to pavement performance modeling to improve the predictive accuracy of predictions when there are only limited or incomplete data available.

The literature reveals that there is successful implementation of PMMS for maintaining the road network and making optimal use of resources in developed countries. However, still attempts are made to transfer such systems to developing countries from developed ones based on techniques, methodologies, and resources that are improper to local condition [26]. Successful implementation of PMMS for different road categories in developing countries like India will lead to optimal use of resources and better and safer roads. The objective of the present study is to develop models for urban road network of Patiala (Punjab, India) using genetic programming (GP) to predict the pavement distress deterioration.

2. Data Collection

An efficient decision support system depends on accuracy, integrity, reliability, and the most important completeness of the database. The data have been collected on the selected 16 roads of Patiala, and the database has been organized in such a manner that the data will be retrieved easily. A detailed description of the collected data and the procedure of data collection and equipments used for the PMMS have been presented in the subsections.

2.1. Identification and Selection of Patiala City Road Sections

The very first step is to identify the road network and select the road sections. In the present study, Patiala City (Punjab, India) road network has been identified and 16 road sections have been selected for developing pavement distress models. Figure 1 shows a map of urban road networks of Patiala city, and Table 1 consists of inventory data of selected road sections. The selected road sections are homogeneous within them but vary considerably from each other in terms of traffic and pavement width.

2.2. Traffic Volume Data

In the present study, traffic volume data have been collected from Municipal Corporation, Patiala. Table 2 consists of the traffic volume of each road section, and Table 3 shows volume-to-capacity ratio of the selected road network. The Equivalent Standard Axle (ESA) repetition in the analysis year (in millions) has been calculated as per IRC 37 [31].

2.3. Structural Evaluation of Pavements

The structural evaluation has been undertaken to assess the pavement’s structural ability to receive wheel loads plying over it, that is, measurement of rebound deflection. The normal practice is to use the Benkelman beam deflection (BBD) method for evaluating the structural condition of the flexible pavement as per the procedure laid down in IRC [27]. The dual rear wheels of the truck are centered above the selected point. The probe of the Benkelman beam is inserted between the dual wheels and placed on the selected point. The initial reading (D0) is recorded when the rate of deformation of the pavement is equal or less than 0.025 mm/minute. The truck is slowly driven a distance of 2.70 m. An intermediate reading (Di) is recorded when the rate of recovery of the pavement is equal or less than 0.025 mm/minute. The truck is driven forward a further 9 m. The final reading (Df) is recorded when the rate of recovery of pavement is equal to or less than 0.025 mm/minute. The rebound deflections have been calculated for all the selected road sections by taking the initial, intermediate, and final deflection readings of the BBD test. The rebound deflection value at any point is given by the following two conditions:(i)If (DiDf)  0.025 mm,then rebound deflection = 2 (D0Df).(ii)If (DiDf) > 0.025 mm,then rebound deflection = 2 (D0Df) + 2.91 [2 (DiDf)].

Pavement temperature has been recorded at least once every hour by inserting thermometer in the standard hole and filling up the hole with glycerol. The rear axle weight of the truck is 8170 kg with a tyre pressure of 5.6 kg/cm2, and spacing between the tyre walls is 35 mm. Test points are taken at a distance of 60 cm from the pavement edge if the lane width is less than 3.5 m and at 90 cm distance from the pavement edge for wider lanes. For divided four lane stretch, the measurement points should be at a distance of 1.5 m from the pavement edge. Since the deflections measured by the Benkelman beam (BB) are influenced by the pavement temperature and seasonal variations in different climates, pavement temperature and soil subgrade details have also been collected at all observation points for making subsequent corrections to the deflection values. The characteristic deflection value is obtained as a sum of the mean and standard deviation for all the corrected rebound deflection values for a particular section. From the BBD values determined as above, the modified structural number (MSN) for all the pavement sections has been calculated [30]:

Figure 2 shows the BBD survey on UR02, and Figure 3 presents the progressive variations of BBD observed during the years 2012–2015 on the selected road sections of Patiala City, that is, UR01 to UR16.

2.4. Functional Evaluation of Pavements

Functional evaluation of pavements consists of collection of road condition data related to surface distress including pavement distress in terms of cracks, potholes, rutting, and ravelling. In the present study, cracking, potholes, ravelling, rutting, and roughness have been measured for all the selected road sections for the development of prediction models for pavement distress. The two other parameters, skid resistance and mean texture depth, have been also measured. Skid resistance has been measured in terms of skid resistance value and sideways force coefficient (SFC). But, there has not been any noticeable decrement in the value of SFC from the year 2012 to 2015. Similarly, there has not been any significant change in the mean texture depth measured using the “Sand Patch Method.”

2.4.1. Cracking Measurement

A number of representative test sections of length 100 m have been selected for cracking measurements for each pavement section. Cracking (alligator, longitudinal, and transverse) has been visually inspected. In case of alligator cracking, the area covered under the distress has been marked in the form of a rectangular box with a chalk on the ground and is measured with a tape. In case of longitudinal and transverse cracks, effective width has been taken as 50 cm and actual length has been measured. Cracked area is expressed as the percentage of the total pavement area. Figure 4 shows the cracks present on the road section UR09, and Figure 5 shows the progressive variation of the cracking area for the road sections UR01 to UR16 during the years 2012–2015.

2.4.2. Pothole Measurement

Pothole means open cavity on the pavement. Measurement of potholes for each road section has been done visually. Minimum diameter of 150 mm and minimum depth of 25 mm of pothole have been taken [28]. Figure 6 shows the pothole on the section UR04, and Figure 7 shows the progressive variations of potholes for the road sections UR01 to UR16 during the years 2012–2015.

2.4.3. Ravelling Measurement

Ravelling is the loss of material from wearing surface. This distress type is associated with thin surfacing, such as, surface dressing, seal coat, and premix carpet. The affected area has been measured by taking into account the area enclosed in regular geometric shapes, and then it has been expressed as the percentage of the total pavement area. Figure 8 shows the ravelling on the section UR13, and Figure 9 shows the progressive variations of ravelling for the road sections UR01 to UR16 during the years 2012–2015.

2.4.4. Rutting Measurement

Rutting refers to permanent deformation on the pavement surface in the wheel path. Permanent deformation is commonly defined as distortions within the individual layers of the pavement structure. Rutting is caused by structural and/or material failure that can lead to hazardous driving conditions. It is very common in flexible pavements. The rut depth has been measured with 2 m straight edge placed in the transverse direction. Twelve observations have been carried out for each selected road section. The maximum rut depth (in mm) has been taken in each observation slot. The rut depths have been averaged to get the mean rut depth of the selected road section. Figure 10 shows the rutting on the section UR01, and Figure 11 shows the progressive variation of rutting for the road sections UR01 to UR16 during the years 2012–2015.

2.4.5. Roughness Measurement

Road roughness refers to surface irregularities in the longitudinal direction. Roughness is an important pavement evaluation parameter because it affects not only ride quality but also vehicle delay costs, fuel consumption, and maintenance costs. Roughness has been measured with the fifth wheel bump integrator or simply known as roughometer. Unevenness Index (in cm/km) has been calculated for each section with the following equation:

Roughness has been calculated in the form of International Roughness Index (IRI, m/km). Above Unevenness Index (UI) value has been converted into International Roughness Index (IRI in m/km) by using the following equation [28]:where is the Unevenness Index in mm/km. is the International Roughness Index in m/km. Figure 12 shows bump integrator test (roughness test) in progress on the section UR11, and Figure 13 shows the progressive variations of IRI for the road sections UR01 to UR16 during the years 2012–2015.

3. Methods

GP has commenced as an endeavor to discover how computers could be trained to solve problems without being explicitly programmed to do so. GP is an extension of genetic algorithms proposed by Koza [29] who defines GP as a domain independent problem solving impends in which computer programs are evolved to solve, or approximately solve, problems based on Darwinian principles of reproduction and analogs of naturally occurring genetic operations such as reproduction, crossover, and mutation. The five major prelude steps for the basic version of GP require the user to specify (a) the set of terminals for each stem of the to-be-evolved program; (b) the set of primitive functions for each stem of the to-be-evolved program; (c) the fitness measure; (d) certain parameters for controlling the run; and (e) the termination criterion and method for designating the result of the run.

3.1. Steps to Execute GP

The execution steps of GP are as follows:(1)Select the dataset that has to be used as initial population (generation 0) with the available functions and terminals.(2)Apply the following steps iteratively on the population until the termination criterion has been satisfied:(a)Execute each program in the population and ascertain its fitness (explicitly or implicitly) using the problem’s fitness measure.(b)Select one or two individual(s) from the population based on the probability of fitness measure to participate in the genetic operations in (c).(c)Generate new individual(s) for the population with the help of following genetic operations:(i)Reproduction: copy the selected individual to the new population.(ii)Crossover: create new children for the new population by recombining randomly chosen parts from two selected individuals.(iii)Mutation: create one new child for the new population by randomly mutating a randomly chosen part of one selected individual.(iv)Architecture-altering operations: choose operation from the available repertoire of such operations and create one new child for the new population by applying the selected architecture-altering operation to one selected individual.(3)When termination criterion has been satisfied, the best program in the population produced during the run has been designated as the result of the experiment. If the experiment is successful, the result may be a solution (or approximate solution) to the problem.

3.2. Development of the GP Model

The control parameters as given in Koza [29] have been investigated in this study. The population size () and the number of children produced (λ) have been taken as 200 and 250, respectively. The larger is the number of generations, the greater are the chances of evolving a solution; hence, the number of generations has been taken as 100000. The values for the parameters, namely, crossover rate and mutation rate, have been selected as 0.80 and 0.04, respectively, after experimentations. The values of other parameters, that is, training percentage, selection method, and tournament size of substitution, are taken as 70, “tournament” and 3, respectively, for the development of the GP model. Function set (+, −, ∗, /, sqrt(), tanh(), exp(), log(), fabs(), pow(x, y), pow(x, 2)) has been taken for generation of model equations. The architectural details of the selected GP model are given in Table 4.

3.2.1. Models Generated Using GP

In the GP model development, the addition is chosen as the linking utility. The prediction models generated using the GP are listed in Table 5. It consists of five models: Model 1, Model 2, Model 3, Model 4, and Model 5 for cracking, ravelling, potholes, rutting, and roughness, respectively.

4. Results and Discussions

The results obtained from this approach have been used for the quantitative assessment of the model’s predictive abilities and are presented in Table 6. The dataset gathered during the years 2012 and 2013 has been used for training models; that is, for the prediction of the year 2013 pavement distress, the dataset of the year 2012 has been used. The postregression fits for the training dataset of all the models (Model 1, Model 2, Model 3, Model 4, and Model 5) are presented in Figures 1418, respectively. It can be observed from Figures 1418 that the some deviation can be noticed between the observed and predicted values because of the modeling approach adopted and because of the assumed normal construction quality for the bituminous surfacing considering the use of optimum binder content in bituminous base and wearing courses as per the routine practice for maintaining the urban roads in India. But, this deviation is quite reasonable for the urban road sections of different age and traffic loading conditions.

4.1. Validation of GP Models

To further test the efficacy and reliability of the models, the datasets collected during the years 2014 and 2015 have been used for the validation purposes. The results of proposed models for validation dataset of year 2014 and year 2015 are presented in Table 7. For the prediction of pavement distress of the year 2014, the data gathered during the year 2013 has been used, and for the prediction of pavement distress of year 2015, the data gathered during the year 2014 has been used. The postregression fits for the validation datasets (2014–2015) of all the models (Model 1, Model 2, Model 3, Model 4, and Model 5) are presented in Figures 1928, respectively. All the GP models have been predicting with good competence for both the datasets (2014 and 2015). The results are revealing that these GP models can be applied successfully to the roads of Patiala City, Punjab, India.

5. Conclusions

(i)The objective of this study was to explore the applicability of GP models for pavement distress. Five models have been developed for the cracking progression (Model 1), ravelling progression (Model 2), pothole progression (Model 3), rutting progression (Model 4), and roughness progression (Model 5). The data collected during the year 2012 to year 2015 on the 16 selected roads of Patiala, Punjab, India, have been used for the development of these models.(ii)We can say that GP has built a good correlation between observed and predicted values and the parameters (like cracking, ravelling, pothole, rutting, and roughness) that have been taken in the present study are sufficient for the prediction of pavement distress.(iii)The deviations between observed and predicted values are only because of the modeling approach adopted and because of the assumed normal construction quality for the bituminous surfacing considering the use of optimum binder content in bituminous base and wearing courses, but this deviation is quite reasonable for the urban road sections of different age and traffic loading conditions.(iv)These models can be applied to other city roads with a range of traffic varying from 0.1 to 0.85 equivalent standard axle repetitions in the initial analysis year and MSN in the range of 2 to 4.5.(v)As an outcome, we can say that GP models provide good results for the prediction of pavement distress and may serve as a predictive model for the prediction of pavement distress for all the cities having homogeneous traffic and climate conditions like Patiala, Punjab, India.

Conflicts of Interest

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