Abstract

Due to the uncertainty and variability of various factors affecting the pavement performance, the change in pavement performance cannot be completely determined. In addition, this uncertainty is not accurately reflected by the pavement performance prediction model. In particular, the gray GM (1, 1) model is very suitable due to it is ability to better predict the existing situation of a domestic asphalt pavement along with the actual performance of a road surface of the “small sample, poor information” gray system. In this regard, the gray GM (1, 1) model is being increasingly used to forecast the performance of an asphalt pavement. When a gray GM (1, 1) model is used to predict the performance of an asphalt pavement, the condition number of the GM (1, 1) model matrix is too large, which, in turn, leads to the deviation of calculation and even wrong results in some cases. This study analyzed the reason for a large condition number of the GM (1, 1) model matrix. Combined with the numerical characteristics of the pavement condition index (PCI) and pavement quality index (PQI), this study focused on the annual, monthly, and daily attenuations of PCI and PQI to the condition number of the GM (1, 1) model matrix. Accordingly, we propose a method to forecast the performance of an asphalt pavement using the monthly attenuation of PCI and PQI. The PCI and PQI in Hunan Province in recent years have been predicted, and the findings reveal that the prediction GM (1, 1) model for the monthly attenuation of PCI and PQI not only effectively lowered the condition number of the matrix but also ensured that the relative error was small.

1. Introduction

The performance of a road surface encompasses five aspects: driving quality, driving safety, structural strength, damage status, and appearance. Factors that affect pavement performance include traffic volume, structural strength, structural composition, material characteristics, and environmental factors. To that end, Fernandes conducted a study on a new asphalt mixture using waste materials and found that it can be used for road-paving works in order to improve the performance without compromising on human and environmental safety [1]. Correspondingly, Camargo and Bernucci carried out a study in order to characterize the field performance of Brazil’s first pavement rehabilitation job using a field-blended rubber asphalt mixture and attempted to predict its future performance [2]. Similarly, Kim et al. conducted a research to determine the effects of temperature segregation on the density and mechanical properties of Louisiana asphalt mixtures [3], whereas Fostinelli performed a cross-sectional study with the purpose of evaluating occupational exposure to polycyclic aromatic hydrocarbons (PAHs) in workers involved in the pavement construction of a new highway in Northern Italy, where modified bitumen was utilized as a binder for hot mix asphalt [4]. Analogously, Gardete adopted an approach that can be potentially very cost effective in enhancing the performance of the bituminous mixtures and achieving better/more durable pavements [5]. On the contrary, Zadshir et al. investigated the impact of three modifiers (one petroleum-based and two bio-based) on the thermomechanical properties of an oxidized asphalt binder [6]. Meanwhile, Yousefdoost et al. presented the outcomes of the project’s first phase wherein a database of mechanical properties of 28 typical dense-graded Australian asphalt mixes was established, in addition to the rheological properties of their binder constituents [7]. Correspondingly, Tanakizadeh and Shafabakhsh examined existing micromechanical models such as the inverse rule of mixtures, Hashin–Shtrikman model, generalized self-consistent scheme (GSCS), and 4-phase model, with a view to determine the dynamic shear modulus of aged asphalt mastics [8]. Braswell et al. examined permeable pavement performance when built over a low-conductivity clay soil [9], while Santos et al. posited a full process-based comparative life cycle assessment (LCA) in order to obtain a better understanding of the environmental impact of reducing mixing temperature by incorporating warm mix technologies, namely, chemical additive-based and foam-based, along with different rates of recycling (0% and 50% RAP) [10].

The prediction of pavement performance worldwide can be broadly classified as deterministic, probabilistic, intelligent, or other types of models. Due to different conditions, many countries and regions have proposed different forms of predictive models. For example, the empirical statistical regression method, time-series analysis method, autoregressive moving average model method with difference, adaptive control Markov model method, normal Markov chain method, artificial neural network method, BP neural network method, gray methods, gray-Markov combination method, matter element method, and gray matter element method are used to predict pavement performance [1117]. Owing to the underlying uncertainty and variability of various factors affecting the pavement performance, it is not possible to completely determine the change in pavement performance; this uncertainty, however, is not accurately reflected by the prediction model. Since Professor Deng proposed the gray system theory [18], it is being increasingly used in order to predict the performance of an asphalt pavement. Gray theory is a method to study the problem of “less data and poor information” uncertainty. It mainly realizes the correct description and effective monitoring of the system’s operational behavior and evolution law by extracting valuable information from the generation and development of “partial” known information. During the past decade, major universities and scientific research institutes within the highway industry have awarded doctoral or master thesis to studies related to gray system theory in order to predict the performance of an asphalt pavement every year [1932]. The gray GM (1, 1) model is particularly suitable because it can effectively predict the existing scenario of a domestic asphalt pavement and the actual performance of a road surface of the “small sample, poor information” gray system. However, the ill-conditioned predicament of the GM (1, 1) model matrix limits its computational accuracy. More specifically, the condition number of a gray system model matrix is larger, which poses a challenge. Many scholars do not consider the influence of ill-conditioned matrix on prediction results whilst using gray models to make predictions for asphalt pavement performance. After Professor Deng proposed the gray system theory, some scholars have attempted to analyze why the matrix was ill-conditioned or what caused the matrix to be ill-conditioned [3338]. Against this backdrop, the main objective of this paper is to analyze the condition number of the GM (1, 1) model matrix which causes matrix ill-conditioning. Taking the GM (1, 1) model as an illustration, the rationale behind a large condition number of the GM (1, 1) matrix was analyzed by combining the numerical characteristics of pavement performance indicators, such as the pavement condition index (PCI) and pavement quality index (PQI). Notably, the size and impact of the PCI and PQI’s annual, quarterly, monthly, half-month, ten-day, and daily attenuations on the GM (1, 1) model matrix condition number were evaluated in this study. In addition, a method is postulated to predict the performance of an asphalt pavement based on the monthly attenuation of PCI and PQI. This method not only is more accurate than the existing prediction model of asphalt pavement performance but also effectively reduces the condition number of the gray model matrix.

2. Research Significance

Premised on the GM (1, 1) prediction model of asphalt pavement performance, this paper analyzes the size and influence of the condition number of the annual attenuation value, quarterly attenuation value, monthly attenuation value, half-month attenuation value, ten-day attenuation value, and daily attenuation value of PCI and PQI on the model matrix. Importantly, the GM (1, 1) model matrix condition number formed by the monthly attenuation value is the smallest after comparing the GM (1, 1) model matrix condition numbers that are formed by the PCI and PQI’s annual attenuation value, quarterly attenuation value, monthly attenuation value, half-month attenuation value, ten-day attenuation value, and daily attenuation value. To that end, a GM (1, 1) model for predicting the performance of an asphalt pavement with the PCI and PQI’s monthly attenuation value as a nonnegative smooth sequence is established, which then reduces the ill-conditioned number of the GM (1, 1) model matrix. Additionally, the relative error of the monthly attenuation value GM (1, 1) prediction models of the PCI and PQI is also smaller than that of other asphalt pavement performance prediction models.

2.1. Reason for a Large Condition Number of GM (1, 1) Model Matrix

Definition 1. Set as a nonnegative smooth sequence ; is ’s 1-AGO sequence, denotes the neighboring mean generation sequence, whereas is considered a gray differential equation [18]. It is also known as the definition of the GM (1, 1) model.

Theorem 1. Let be a GM (1, 1) model parameter, and

Then, the least-squares estimation parameter sequence of the gray differential equation isthat is, . The coefficient matrix is as follows [16]:

Reference [34] makes the assumption that, concerning the situation of , when , the maximum condition number value of is signified by . When , the condition number of is predicated on the extent to which is close to .

When , the matrix witnesses a severe drift. Notably, the gray system studies small samples; however, in actual production and life, the value of varies significantly with actual production and life issues. If or above, the value of is also larger. In other words, the condition number of the matrix is correspondingly larger, which, in turn, leads to the morbid matrix. When , the condition number of is predicated on the extent to which is close to n − 1. According to example 6.1 in reference [18],

Currently, the condition number of the GM (1, 1) model matrix is 249098.48. The criteria of an ill-conditioned matrix’s conditional number size are as follows: good condition (<10), mild ill condition (10, 100), strong ill condition (100, 1000), and severe ill condition (>1000) [38]; meanwhile, is severely ill. Different ill-conditioned degrees of the matrix will result in an unstable or even wrong solution of the model.

2.2. Number of GM (1, 1) Matrix Conditions with Different Attenuations of PCI and PQI

Section 2.1 points out that if or above, when , then the value of is also larger. In terms of the performance indicators of an asphalt pavement, the road surface condition index (PCI) and the driving quality index (PQI) are usually 100 points. Table 1 depicts the pavement performance indexes of three representative lines of main highways in Hunan over the past few years [31].

In the event that the detected values of PCI and PQI are directly used as the nonnegative smoothing sequences within the GM (1, 1) model, the resulting number of matrix conditions is very large. Table 2 illustrates the condition number of the GM (1, 1) model matrix obtained by directly using the values of PCI and PQI.

Considering the PCI of line 1 as the sequence, the GM (1, 1) model matrix is , whereas the condition number is 1002500. According to Table 2, the condition number of the matrix obtained by directly using the values of PCI and PQI is far beyond the range of severe ill condition (>1000).

2.2.1. Condition Number of Matrix Based on the Annual Attenuations of PCI and PQI

Asphalt pavement performance indicators are inclusive of PCI and PQI. In the initial period without any maintenance, their values reduce each year from the beginning of operation; in other words, the road surface condition and driving quality deteriorate by the year. Meanwhile, the annual attenuations of PCI and PQI can be defined as the PCI and PQI values for the previous year minus those pertaining to the current year. Taking the example of the period from 2013 to 2014, these attenuations can be obtained by subtracting the PCI and PQI values for 2014 from the PCI and PQI values for 2013. Owing to the year-by-year decay nature of PCI and PQI values, this approach ensures the annual attenuation satisfies the requirement that is a nonnegative smooth sequence in the GM (1, 1) model. Meanwhile, Table 3 depicts the PCI and PQI’s annual attenuations of three representative lines of the main highways in Hunan Province listed in Table 2.

After obtaining the annual attenuations of PCI and PQI, the annual attenuation was used as the nonnegative smooth sequence of in the GM (1, 1) model. The model parameter was solved using equations (1) and (2). Table 4 illustrates the condition number of the GM (1, 1) model matrix that was obtained using the PCI and PQI’s annual attenuations.

The PCI’s annual attenuation of line 3 denotes the sequence. The GM (1, 1) model matrix is , whereas the condition number is 1318.2. Importantly, the condition number of the GM (1, 1) model matrix derived using the annual attenuation is partially in a severely ill-conditioned state (>1000) and partly in a relatively strong ill-conditioned state (100, 1000).

2.2.2. Condition Number of Matrix Based on the Monthly Attenuations of PCI and PQI

The monthly attenuations of both PCI and PQI can be defined as the annual attenuations of PCI and PQI evenly distributed to 12 months per year; that is, monthly attenuation of PCI and PQI = annual attenuation of PCI and PQI/12. Table 5 illustrates the monthly attenuations of PCI and PQI for three typical representative line pavements in Hunan Province.

After obtaining the monthly attenuations of PCI and PQI, the monthly attenuation was considered the nonnegative smooth sequence within the GM (1, 1) model. The model parameter was solved using equations (1) and (2). Table 6 illustrates the condition number of the GM (1, 1) model matrix that was obtained using the monthly attenuations of PCI and PQI.

The monthly attenuation of PQI of line 2 represents the sequence. The GM (1, 1) model matrix is , whereas the condition number is 21.02. It may be important to note that the condition number of the matrix obtained using the monthly attenuation lies in the slightly ill-conditioned range (10, 100).

2.2.3. Condition Number of Matrix Based on the Daily Attenuations of PCI and PQI

The daily attenuations of PCI and PQI can be defined as the annual attenuations of PCI and PQI evenly distributed to 365 days per year; that is, monthly attenuation of PCI and PQI = annual attenuation of PCI and PQI/365. Table 7 depicts the daily attenuations of both PCI and PQI for three typical representative line pavements in Hunan Province.

After obtaining the daily attenuations of PCI and PQI, the daily attenuation was utilized as the nonnegative smooth sequence within the GM (1, 1) model, whereas the model parameter was solved using equations (1) and (2). Table 8 illustrates the condition number of the GM (1, 1) model matrix that was obtained using the daily attenuations of PCI and PQI.

The monthly attenuation of PQI of line 3 signifies the sequence. The GM (1, 1) model matrix is , whereas the condition number is 4406.3. The condition numbers of the matrix obtained using the daily attenuations are in the severely ill-conditioned range (>1000).

A comparison between the data shown in Tables 38 revealed that when taking the annual attenuation of PCI or PQI, the condition number of the GM (1, 1) model matrix is in the strong and severely ill-conditioned range. Meanwhile, when considering the daily attenuation, the matrix condition number lies in the severely ill-conditioned range. Finally, the matrix condition number is in the slightly ill-conditioned range when taking the monthly attenuation.

According to the authors’ calculations, if the quarterly, half-month, and ten-day attenuations of PCI or PQI are used, the resulting condition number of the matrix is higher than that of the monthly attenuation.

For this reason, the monthly attenuation of PCI or PQI was used in this paper as the initial nonnegative smoothing sequence in the GM (1, 1) model in order to forecast the performance of an asphalt pavement.

3. Prediction Model of Asphalt Pavement Performance Using Reduced Condition Number of Matrix

Considering the matrix calculated in Section 2 using different values of PCI or PQI, it can be easily observed that, since the GM (1, 1) model matrix is a symmetric matrix, the value corresponding to the second row and second column of the matrix is n − 1. The calculation theory of condition numbers demonstrates that it is only when the other three numbers of the matrix and n − 1 are an order of magnitude that the condition number will be reduced below the level of mild ill condition. To illustrate this, in the matrix that was obtained from the annual attenuation, n − 1 has a value of 5, which is only a single digit, whereas the others are at least a hundred, sometimes even reaching a thousand 1329.4. In the matrix obtained from the monthly attenuation, n − 1 has a value of 4, which again is a single digit. Additionally, the other three values are all single digits as well; therefore, the condition number is not large. During the usage of the daily attenuation, the magnitude of difference between the other three numbers and n − 1 became larger, thereby increasing the condition number accordingly.

Using the monthly attenuation of PCI or PQI, the GM (1, 1) model used for predicting the performance of an asphalt pavement is as follows:

Let be the annual value of a nonnegative smoothing sequence of order . denotes the 1-AGO sequence of ; refers to the close neighboring mean generating sequence of ; and signifies the gray differential equation. Let be the GM (1, 1) model parameter, and

The least-squares estimation parameter sequence of the gray differential equation is

Both equations (5) and (6) are GM (1, 1) prediction models for asphalt pavement performance; they are capable of reducing the condition number of the matrix. Each value in the sequence denotes the monthly attenuation of PCI or PQI.

4. Prediction of Asphalt Pavement Performance of the Main Highway in Hunan Province

In the trunk highways of Hunan Province, the evaluation of the performance of asphalt pavements was carried out in accordance with the standard “Highway Performance Assessment Standards” (JTG H20-2007) (China). For the roads of Grade II and below Grade II, the road pavement technical condition evaluation index adopts the pavement quality index (PQI), which consists of two individual indicators, namely, the pavement condition index (PCI) and the pavement driving quality index (RQI).

For the trunk highways, the pavement damage and the flatness technical condition are detected by automatic detection equipment, that is, the multifunctional rapid comprehensive detection vehicle. The equipment can not only identify various types of pavement diseases specified in the “Highway Performance Assessment Standards” (JTG H20-2007) but also meet the requirements of the specifications. After the detection vehicle is detected, the pavement damage image is recognized by the artificial assistant software, and the PCI is obtained. The pavement flatness data are analyzed and processed by special software to obtain the driving quality index (RQI) The detection vehicle and software to obtain PCI and RQI are shown in Figure 1.

The pavement quality index (PQI) is calculated according to the following equation:where PCI is the pavement condition index and RQI is the pavement driving quality index.

The pavement condition index (PCI) is calculated according to the following equations:where DR is the pavement breakage rate (%), the sum of the various damaged areas and the percentage of the surveyed pavement area; is the i-type area of a pavement damage; is the survey of the pavement area; is the weight of the i-type damage to the pavement; and is the type of pavement damage in the item considering the degree of damage.

The pavement smoothness is calculated by the pavement driving quality index (RQI), according to the following equation, and the test data are stored for a long time in units of 100 m:where IRI is International Flatness Index, is 0.026 for primary roads and 0.0185 for other grades of highways, and is 0.65 for the first-class roads and 0.58 for other grades of roads.

The provincial highway of Hunan Province has witnessed rapid-paced development over the past decade. With the drastic increase in traffic pressure, the pavement condition index of main highways has decreased year by year. The main highways have entered extensive maintenance and overhaul periods in recent years. The determination of maintaining the asphalt pavement and timing of major repairs, as well as the allocation of maintenance funds, necessitates a clearer understanding about the future changes in pavement conditions.

The representative lines of the main trunk highway are secondary (technical level) and heavy traffic (traffic level), and the overhaul plan is divided into “white to black” and “black to black.” The construction of each layer is carried out after the previous structural layer meets the requirements, but the influence of different original pavements on the performance of the asphalt pavement will be different. Therefore, when the representative lines of the main trunk highway are selected, they are classified according to the overhaul plan. Different structural materials are one of the main factors affecting the performance of an asphalt pavement and should also be classified when selecting representative lines of the main trunk highway. Due to numerous trunk highways in Hunan Province and the complicated maintenance history, the most general situation was selected for comprehensive statistics.

The PCI and PQI of three typical representative routes were predicted in this study based on the monthly attenuations of PCI or PQI. Tables 9 and 10 illustrate pavements of the Hunan main highway based on the monthly attenuations of PCI and PQI.

To illustrate this, the predicted monthly attenuation of PCI of line 1 was calculated as X(k + 1) = −8.7404 × exp(−0.075274 × k) + 8.8954. The predicted value is X(k + 1) − X(k). Suppose X(6) = −8.7404 × exp(−0.07527 × k) + 8.8954 = 2.8964 and X(5) = −8.7404 × exp(−0.075274 × k) + 8.8954 = 2.4274. Then, X(6) − X(5) = 0.4690. It denotes the predicted monthly attenuation of PCI of line 1 in 2018. Similarly, according to the aforementioned method, the predicted monthly attenuations of PCI of lines 2 and 3 were calculated after the prediction calculation. Among them, the predicted monthly attenuation of PCI of line 2 was calculated as X(k + 1) = −4.0328 × exp(−0.13542 × k) + 4.2336. Likewise, the predicted monthly attenuation of PCI of line 3 was calculated as X(k + 1) = −36.1908 × exp(−0.012024 × k) + 36.3374.

The predicted monthly attenuation of PQI of line 1 was calculated as X(k + 1) = −8.0464 × exp(−0.058269 × k) + 8.2189. Similarly, the predicted monthly attenuation of PQI of line 2 was calculated as X(k + 1) = −9.2576 × exp(−0.041277 × k) + 9.3943. Likewise, the predicted monthly attenuation of PQI of line 3 was calculated as X(k + 1) = 38.6564 × exp(0.0079912 × k) − 38.5298.

The predicted monthly attenuation value of PCI or PQI was then multiplied by 12 and converted into the annual attenuation. Subsequently, as illustrated in Tables 11 and 12, the predicted PCI or PQI of the current year was obtained by subtracting the measured PCI or PQI values of the previous year.

Table 13 depicts that the predicted relative error of line 2 in this study is comparable with that of the standard decay equation [11] postulated by the Sun Lijun team of Tongji University as well as the method proposed in the document [31]. The findings demonstrate that the prediction capability in this paper is better than that of the other two methods. As compared to the prediction errors of lines 1 and 3, the model proposed in this paper is also much better than Sun Lijun model as well as the model postulated in the literature [31]. However, only the comparative results of line 2 are shown in Table 13 due to space constraints.

During the last part of early 2018, the relative error of PCI prediction of line 2 was shown to be 0.33 in the Sun Lijun model, 0.92 in reference [31], and 0 in this paper. For the same period, the relative error of PQI prediction of line 2 was 0.16 in the Sun Lijun model, 0.63 in reference [31], and −0.028 in this paper. Notably, if the absolute value of relative error is compared with the PCI and PQI of all the other years, this model forecasts the absolute value of relative error, most of which are much smaller than those of the models proposed by Sun Lijun and reference [31].

5. Conclusion

The predicament of predicting the performance of an asphalt pavement is confronted during the precise design and maintenance of durable asphalt pavements. The foundation for solving this problem is determining the type of model to be used for prediction in addition to the capability of the model. This paper focuses on the GM (1, 1) model which is extensively used in asphalt pavement performance prediction and draws the following conclusions:(1)The direct phenomenon that renders the matrix conditional number excessively large in the gray GM (1, 1) model is too large for the difference between the other three numbers in the matrix and n − 1.(2)Using the original values of PCI and PQI as the nonnegative smoothing sequence in the GM (1, 1) model, the generated matrix entails a large number of conditions and forms an ill-conditioned matrix.(3)After comparing the condition number of the GM (1, 1) matrix formed by the annual, quarterly, monthly, half-month, ten-day, and daily attenuations of PCI and PQI, the condition number of the GM (1, 1) matrix formed by the monthly attenuation was found to be the minimum.(4)The GM (1, 1) model for forecasting the performance of an asphalt pavement using PCI and PQI is premised on the monthly attenuation as the nonnegative smoothing sequence.(5)The relative error of the monthly attenuation of the GM (1, 1) prediction model is also smaller than that of the other models that have been proposed by references.

Since PCI and PQI adopt the 100-point system in China, when using the annual PCI and PQI, the condition number of the matrix is large in the GM (1, 1) model prediction. When using monthly attenuations of PCI and PQI, the condition number of the GM (1, 1) model matrix is small. So this method is applicable to the analysis of PCI and PQI raw data of the pavement for the Chinese Hunan Province conditions.

The GM (1, 1) prediction model based on the monthly attenuations of PCI and PQI not only reduces the number of matrix conditions but also improves the accuracy of prediction and lowers the relative prediction error. To put it succinctly, this model provides a new and effective way of predicting the performance of an asphalt pavement.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The project was supported by the National Key Research and Development Program of China (Grant No. 2017YFC0805300), Research Foundation of Education Bureau of Hunan Province, China (Grant No. 15B010), and Open Fund of State Engineering Laboratory of Highway Maintenance Technology, Changsha University of Science & Technology (Grant No. kfj140102).