Research on Slope Deformation Prediction Based on Fractional-Order Calculus Gray Model

Civil Engineering College, Chongqing ree Gorges University, Chongqing 404100, China Jiangxi Engineering Research Center of Water Engineering Safety and Resources Efficient Utilization, Nanchang 330099, China Changjiang River Scientific Research Institute of Changjiang Water Resources Commission, Wuhan 430010, China College of Environment and Civil Engineering, Chengdu University of Technology, Chengdu 610059, China School of River and Ocean Engineering, Chongqing Jiao-tong University, Chongqing 400074, China


Introduction
Slope disaster prediction and forecasting is the frontier problem and research focus of current environmental geology research and is also one of the critical aspects for people to understand the natural environment [1,2], with regard to geological balance of the surface, which has wide practical significance and important theoretical value.Deformation is the macroscopic and visible physical information in the process of slope evolution, and predicting slope deformation can better understand the evolution behavior of slope.erefore, scholars have done a lot of research on slope deformation prediction and model.Generally speaking, based on the existing deformation monitoring data, nonlinear prediction is often discussed from the perspective of time series statistics or formation mechanism.Nowadays, it is more popular to apply the mathematical model for nonlinear comprehensive prediction combined with the existing data, which can further interpret the internal mechanism of slope instability, further predict the trend of deformation and failure, and provide relevant measures and treatment methods for engineering.In addition, there have been many successful cases for guiding practical projects.
At present, the slope deformation prediction is based on slope displacement monitoring data, reservoir water level change, rainfall monitoring, etc., and therefore, using certain mathematical methods or means to obtain the slope deformation prediction model is one of the main technical difficulties in the slope engineering.
e mechanical knowledge still cannot solve the problem of slope deformation prediction properly.erefore, mathematical methods or approaches are adopted to predict the deformation of the slope, in order to obtain the deformation trend.Foreign scholars started related research earlier.For example, Jibson [3] adopted a regression model to estimate the displacement of slope under the action of strong vibration and achieved the expected results.Kaund et al. [4] employed the backpropagation neural network model to predict slope displacement and groundwater height.At the same time, Jin et al. [5], Tan et al. [6], Zhao et al. [7], Li et al. [8], Peng et al. [9], and other domestic scholars proposed the models of slope deformation prediction in terms of slope instability and deformation mechanism, thus obtaining the optimal combination weight coefficients according to forecast effectiveness, residual error, and other parameters and carried out fitting to learn deformation prediction in order to enhance the prediction accuracy; Zhang et al. [10] and Ma [11] used the vector machine in mathematics to predict the slope deformation of Baima River and Ancient Tree House, which is almost the same as the actual measurement.Moreover, in respect of the prediction of landslide displacement, Huang, Cao, Zhu, Zhou, Xie, and Cai [12][13][14][15][16][17] adopted some mathematical models and intelligent methods to predict the landslides in Sichuan and ree Gorges Reservoir Area in China, which provides relevant nonengineering measures for landslide prevention and control and achieved good results.ese methods used in slope or landslide deformation generally have their own advantages and disadvantages, and often different methods get different prediction results.e main differences lie in terms of the selection singleness of methods, the effectiveness of monitoring data, and the acceptability of error judgment.
Using mathematical models and intelligent evaluation methods to accurately predict slope deformation and at the same time integrating the mechanical theory interpretation with the deformation prediction has important theoretical significance and practical engineering value.e fractionalorder calculus method is an extension of the classical calculus method with memory and genetic functions, especially in the complex system description; it boasts the advantages of simple and fast modeling, easy to implementing, and accurate positioning of description and has become one of the most important tools for mathematical modeling of complex mechanical and physical processes.ere are relatively few studies in this area, and among the engineering application foundations, there are a few reports on fractional creep and fractional plasticity.Wang et al. [18], Chen et al. [19], He et al. [20], and Wu et al. [21] came up with the fractional calculus theory to propose a nonlinear creep model for rock objects and conducted related experimental research; Tang et al. [22] put forward the creep damage model of loses with variable parameters based on fractional calculus.
ese studies suggest that the fractional-order calculus method can be used in geotechnical engineering research.In addition, the gray theory has the advantages of moderate sample size, quantitative qualitative combination, small computational workload, high prediction accuracy, and high accuracy.erefore, it is a mathematical model commonly used in deformation prediction.Whether the advantages of both gray theory and fractional-order calculus can be combined and used in practical engineering prompt the necessity of research.In this paper, the fractional-order calculus method is used to derive the fractional calculus gray prediction model based on the time series of deformation data.In order to verify its rationality, taking the monitoring data of a deep foundation pit in Chongqing, Southwest China, as the engineering case example, the slope deformation prediction calculation and analysis are carried out to provide technical support for the slope deformation treatment and maintenance.

Derivation of Fractional-Order Calculus Gray Model of Slope Deformation Prediction.
Assume that the slope deformation monitoring data are X * � (X * (1), X * (2), • • • , X * (n)), and use the appendix Lemma 1 and Lemma 2 to establish a fractional gray difference formula: where α (1) Derivate parameters a and b according to the least squares method (see Appendix Lemma 4 for details): e time response of the descendant into the fractional gray model (see Appendix Lemma 5 for details) is where t is the response time and the unit is the day.e above formula is a fractional gray model for slope deformation prediction.
2 Advances in Civil Engineering

Implementation of Fractional Calculus Gray Model for
Slope Deformation Prediction.According to the above mathematical lemma and derivation, the detailed steps of the fractional-order calculus gray model of slope deformation are as follows: (1) First of all, process the original displacement data, and subtract the initial value from the existing time series data so that the initial value is 0. is is to directly apply the solution formula of Lemma 5.In this step, it should be noted that the postprocessing time series data cannot be negative because usually the cumulative displacement for slope deformation prediction is adopted, and thus, it can be satis ed.(2) Determine the order of the fractional order, and apply the parameters of the lemma gray model.(3) Substitute the parameters into the solution of the albinism di erential equation of the gray model, adding the initial values to all the simulated values, thus obtaining the predicted value of the slope deformation.
In the above derivation, the derivative order p value needs to be arti cially given.In order to obtain a reasonable p value, an optimization algorithm (such as the gradient method, genetic algorithm, and particle swarm algorithm) is used to solve the problem.erefore, the solution of the p value can be expressed as min where X (exp) (t) is the monitoring data and 1 represents rst-norm.
In order to facilitate visual reading and make it easy to use, the Matlab toolbox is called to compile the corresponding program, and the implementation process is shown in Figure 1.

Engineering Example of Slope Deformation
Prediction of a Deep Foundation Pit

Survey of Engineering Geology of the Deep Foundation Pit
Slope. e deep foundation pit slope project of an international garden in Chongqing is located in the Northern New District of Chongqing.It covers an area of about 1,200 mu (1 mu ≈ 666.66 m 2 ), of which the pure water surface area of the lake is about 160 mu. e terrain falls on the lakeside like a gentle slope, wherein Region B covers 300 mu, which is located in the central lakeside hinterland of the project.
is region belongs to the subtropical monsoon climate area, its annual average rainfall is 1100 mm, and the rainfall is mostly concentrated from May to September, taking up 70% of the annual rainfall; the long-time average annual temperature is 18.6 °C, the extreme minimum temperature is 4.5 °C, and the extreme maximum temperature is 42.5 °C; the long-time average annual relative humidity is 80%, and the annual distribution is the largest in December, namely, 87%, and the smallest in August, namely, 74%; and the frost period generally lasts 10∼20 days, and fog days are up to 20∼35 days, and the number of sunshine hours is as many as 1384.2∼1542.8hours.
e area is a type of denuded hilly landform, with no faults passing through it.
e occurrence of rock formation is SW310 °∠9 °in the exposed area of the bedrock near the area, and there are two sets of ssures in the rock mass.e structural ssures are not developed, with no faults passing through it, and the development degree of ssures is simple.
e slope angle of the terrain is 2∼28 °; the average slope angle is about 22 °; the topography and geomorphology are relatively complex; the dip angle of the rock layer is gentle, with no faults passing through it; the ssures are not developed; the rock mass takes on the shape of the medium thick to thick layer; the lithology and soil layer combination are binary combination, and the exposed bedrocks are the mud rocks and sand stones of Middle Jurassic Upper Shaximiao Formation (J2s); the quaternary soil layer includes el-dlQ (Q4el + dl), Qml (Q4ml), and Qapl (Q4al + pl); the soil layer is 0.50 to 4.7 m thick (2.7 m on average), and the dip angle of the geosynthetic interface is 3∼14 degrees; and the local geotechnical interface is forward and nonemptive, and dip size is ≥20 °, accounting for 28.8% of the planning area.e relationship between the rock mass penetrating the structural plane and the oblique (edge) is relatively complicated, and the in uence of surface water and groundwater on the rock mass is small, the basic intensity of earthquakes of the VI degree.Human activities that a ect the geological environment are not strong.
ere are no bad geological phenomena such as dangerous rock collapse, landslide, debris ow, and ground collapse, and thus, the geological environment is complex to some degree.
According to the engineering geological conditions, three kinds of protection methods are used to support and protect the slope of the deep foundation pit, such as protection treatment of the lattice anchor rod or ordinary slideresistant pile or slide-resistant pile with prestressed anchor cable by combining with the C25 shotcrete panel.In the process of construction, monitor the horizontal and vertical displacements of the slope top, monitor the anchor rod and prestressed anchor cable press, and make a macroscopic inspection of the Earth's surface.

Acquisition and Analysis of Deformation Monitoring
Data. e safety grade of the foundation pit slope in Region B of an international garden in Chongqing is Grade I, the whole excavation panorama of the deep foundation pit slope is shown in Figure 2(a), and the monitoring plane graph is shown in Figures 2(b) and 2(c).During the construction process, before the foundation pit is excavated to 5 m, monitor at the frequency of once every 7 days; before the whole foundation pit is excavated from 5 m to 10 m, monitor at the frequency of once every 3 days; and when the excavation of the whole foundation pit slope exceeds 10 m, monitor at the frequency of once every day.e horizontal displacement and the displacement of the pile top are tested with a total station instrument by means of the photoelectrode coordinate method, and the datum point needed to monitor the deformation observation point accords with second-order accuracy.At the time of settlement observation, the error in the observation station's height difference of the observation point is ≤0.5 mm and in the observation of plane displacement, the error in the coordinates of observation points is ≤3 mm.Two observation points are set on Building No. Six, the ird Phase of the southern slope of the foundation pit in this region, and they are named as JC1 and JC2. e measured deformation values of the JC1 and JC2 monitoring points are all 5 m in foundation pit depth, and the monitoring time is concentrated from June to September, during which the seasonal rainfall is relatively obvious.erefore, the change in the range of the reservoir water level in the central lakeside hinterland is relatively large, and the actual monitoring data are shown in Figure 3.

Analysis of Deformation Prediction
Results.Select the first 50 time series data of the JC-1 monitoring point to build the model, and predict the displacement changes of the following 26 time series; and select the first 30 time series data of the JC-2 monitoring point to build the model, and predict the displacement changes of the following 15 time series.Among them, the JC-1 monitoring point is located at the top of the deep foundation pit, and the deep foundation pit of the JC-2 monitoring point is located at the bottom position.e sudden change of the monitoring point JC-1 at 10∼20 days is mainly due to the local temporary loading effect at the top of the slope.
According to the measured data of JC-1 and JC-2 monitoring points, using the fractional calculus gray prediction model and gray prediction model, the p value was obtained by the genetic algorithm in the calculation, the p value of JC-1 was taken as 0.5, and the p value of JC-2 was taken as 0.6, and the results are shown in Figure 4. From the deformation prediction of JC-1 and JC-2 monitoring points, the deformation prediction value calculated by the fractional-order calculus gray model is smaller than the measured value on the whole, the deformation prediction value calculated by the gray model is larger than the measured value on the whole, the overall variation trends of the deformation prediction values calculated by the fractionalorder calculus gray model and gray model are the same, and there is no apparent singularity.With the increase of the time series of the deformation monitoring point, the absolute error value of the deformation prediction of the fractional calculus gray model and the measured value are gradually increased, and the absolute error value margin of the deformation prediction of the gray model and the measured value are getting larger and larger; the gradually accumulated value is getting bigger and bigger.When calculating the fractional-order calculus gray model, the minimum absolute error of the JC-1 monitoring point is 2.02 mm, appearing in the first half of monitoring time series prediction, the maximum absolute error of it is 3.98 mm, appearing in the last phase of monitoring time series prediction, and the front-to-rear difference reaches 1.96 mm; the minimum absolute error of the JC-2 monitoring point is 0.125 mm, appearing in the middle section of monitoring time series prediction, the maximum absolute error of it is 0.709 mm, appearing in the last phase of monitoring time series prediction, and the front-to-rear difference reaches 0.584 mm; the variation trends and laws of absolute errors of JC-1 and JC-2 monitoring points are basically the same.With the increase of monitoring time series, the absolute error increases gradually, and the difference between the maximum and minimum values of the absolute error is gradually on the rise as well.When calculating the gray model, the minimum absolute error of JC-1 is 3.79 mm, appearing at the first time series point of monitoring time series prediction, and the maximum absolute error of it is 8.58 mm, appearing at the last time series point of monitoring time series prediction, and the front-to-rear difference reaches 4.79 mm; the minimum absolute error of JC-2 is 0.467 mm, appearing at the first time series point of monitoring time series prediction, the maximum absolute error of it is 1.843 mm, appearing at the last time series point of monitoring time series prediction, and the front-to-rear difference reaches 1.376 mm; the variation trends and laws of absolute errors of JC-1 and JC-2 monitoring points are basically the same.Comparing the deformation prediction of the fractional-order calculus gray model and gray model, the absolute error calculated by the gray model is monotonically increasing, while the absolute error calculated by the fractional-order calculus gray model first increases, then decreases, and finally increases.e absolute error of the gray model is one time larger than that of the fractionalorder calculus gray model.In order to make a clearer comparison, the relative errors of the measured data and the predicted calculation data of JC-1 and JC-2 monitoring points are calculated.As shown in Figure 5, the variation trends and laws of relative errors of JC-1 and JC-2 monitoring points are basically the same as that of the absolute errors which will not be elaborated here in detail.

Discussion
(1) Analysis of the influence of the p value on prediction results.e method for determining the 4 Advances in Civil Engineering p value is given above.e impact of different p values on the model prediction results is further analyzed.Figure 6 shows the sensitivity analysis of different p values for the settlement results of the model using JC-2 monitoring points as an example.It can be seen from Figure 6 that different p values have a greater impact on the model prediction results.In general, as the p value increases, the overall value of the prediction result decreases.

6
Advances in Civil Engineering (2) Di erent from the traditional gray model, the fractional-order gray model introduces the derivative order p, which makes the traditional gray model be a special case of the fractional-order gray model.erefore, the fractional-order gray model can better t the original sequence with strong nonlinearity.From the calculation example in this paper, the prediction accuracy of the fractional-order gray model is obviously improved compared with the existed gray model, but it still has a certain gap from the real monitoring value.In addition, the determination of the p value requires the aid of an optimization algorithm, requires more calculation amount of the model than that of the traditional gray model.

Conclusions
(1) Slope deformation is closely related to soil properties and engineering geological environment.
Rainwater erosion and slope load may directly accelerate slope deformation.erefore, the analysis of the development trend of slope deformation needs to consider the internal and external factors.
It is often combined with the existing monitoring data to use mathematical models to carry out nonlinear calculations and predict the trend of change, which can provide technical support for slope management and promote the disaster prevention and reduction.(2) With the fractional-order calculus method, based on the deformation monitoring data time series, from the theoretical derivation point of view, the fractional-order calculus gray model is proposed after analysis and derivation, and the applicable conditions and ranges are nonnegative sequences, and the initial value is 0. e steps of using the fractional-order calculus gray prediction model are given.Using the Matlab language to compile the corresponding program, these corresponding functions can be realized one by one.
(3) Based on the established fractional-order calculus gray model method, the monitoring data of a deep foundation pit slope in Chongqing, Southwest China, was selected as an example to verify the model.e calculation results show that the fractional-order calculus method has better prediction e ect than that of the gray model, and the latter has a maximum relative error of 37.96%, while the former has a maximum relative error of 17.63%.
ese results show that the fractional-order gray model has a certain application prospect in the prediction of mountain slope deformation.Compared with the measured data, the relative errors are within the acceptable range, which can provide reference for similar projects.
e inverse of Equation (10b) is in the form ( Proof completed.

Figure 2 :Figure 5 :Figure 3 :Figure 4 :
Figure 2: Deep foundation pit slope, monitoring layout, and engineering geological section: (a) whole panorama of deep foundation pit slope; (b) monitoring floor plan of foundation pit slope B district treatment project; (c) Building No. Six engineering geological section.

1 Figure 6 :
Figure 6: Parameter p sensitivity analyses of the JC-2 monitoring point.