Research Article  Open Access
Forecasting Cohesionless Soil Highway Slope Displacement Using Modular Neural Network
Abstract
The highway slope failures are triggered by the rainfall, namely, to create the disaster. However, forecasting the failure of highway slop is difficult because of nonlinear time dependency and seasonal effects, which affect the slope displacements. Starting from the artificial neural networks (ANNs) since the mid1990s, an effective means is suggested to judge the stability of slope by forecasting the slope displacement in the future based on the monitoring data. In order to solve the problem of forecasting the highway slope displacement, a displacement time series forecasting model of cohesionless soil highway slope is given firstly, and then modular neural network (MNN) is used to train it. With the randomness of rainfall information, the membership function based on distance measurement is constructed; after that, a fuzzy discrimination method to sample data is adopted to realize online subnets selection to improve the selfadapting ability of artificial neural networks (ANNs). The experiment on the sample data of Beijing city’s highway slope demonstrates that this model is superior to others in accuracy and adaptability.
1. Introduction
Forecasting the slope displacement accurately in the future becomes a very important step to judge the stability of slope. Forecasting the highway slope displacement in the future is a difficult task which involves evaluation of a large number of interrelated variables and factors and is hardly accounted for a model. Factors affecting slope stability can be divided into two main types: internal factors and external inducements. Internal factors include rock type, joints, property, and structure. External inducements include rainfalls, earthquakes, and human factors. Generally speaking, the slope failures are directly caused by the external inducements, and internal factors are the fundamentals of it. Also, most of these factors have obvious randomness, fuzziness, and variability. All factors interleaving together make the influence on the stability of slope, so evaluation criterion which is used to judge the stability of slope should be easily measured and can make clearly manifestations of the change to slope stability.
Displacement is the external expression of deformation process to mass movements of soil and rock. In the past, displacement had been used to compute the prealert, alert, and emergency phases of large rock slides; this method had been validated by collecting and analysing literature data for historical rock slope failures [1]. The forecasting model of slope displacement includes artificial neural networks (ANNs) [2–7], grey theory [8–10], wholeregion method of chaotic time series to slope deformation prediction [11]. The engineering properties of slop exhibit varied and uncertain behaviors, for the complex and imprecise physical processes associated with the formation of soil and rock [12]. So displacement considered as the external expression of slope failure shows evident nonlinear and uncertain characteristics. In order to cope with the complexity of behaviors and the spatial variability of displacement, traditional forms of forecasting models are justifiably simplified. Artificial neural networks (ANNs) refer to a class of artificial intelligence which attempt to mimic the behaviors of the human brain and nervous system, based on the data alone to determine the structure and parameters of the model. It is well suited to complex model problems where the relationship between the model variables is unknown. So, models based on ANNs have been successfully applied on virtually every problem in slope stability. However, rainfalls have important influence on the slop displacement [13–15], the forecasting displacement models of highway slop based on ANNs are focus on the time series variation of displacement, but neglected it.
Modular neural network (MNN) can solve largescale realworld problems by dividing a problem into smaller and simpler subproblems, assigning a network module to learn each of the subproblems, and recombining the individual modules into a solution to the original problem. In this paper, using rainfall information, the membership function based on distance measurement is constructed, firstly, then a fuzzy discrimination method of sample data is adopted on online subnets selection in order to improve the selfadapting ability of artificial neural networks (ANNs). The experiment with the sample data of Beijing city’s highway slope demonstrates that this model is superior to others in accuracy and adaptability.
2. Study the Relationship between Displacement and Rainfall
In order to obtain an accurate forecasting model based on MNN, a theoretical analysis of the relationship between displacement and rainfall has been made. The most common failure mode of cohesionless soil highway slope is landslide, the geological conditions is high permeability weathered layer or slope deposits which covers low hydraulic conductivity rock stratum, and the interface of them is often the potential slip surface. The result of analytic solution from Iverson [16] is used to describe the transient seepage due to rainfall. Then, an elastoplastic analysis was used to analyze the unstable slope based on MohrCoulomb yield criterion and associated flow rule.
Figure 1 shows the force diagram of cohesionless soil highway slope in plane strain condition based on MohrCoulomb yield criterion. The assumptions are as follows:(i)depth of impervious surface and initial groundwater level is denoted by and , both of them are parallel to the slope surface; (ii)the homogeneous isotropic slide body is elasticperfectly plastic material, strictly observing the MohrCoulomb yield and its flow rule; (iii)the source of groundwater recharge is only the rainfall; during this process without considering the loss such as evaporation, the soil above the groundwater level is fully saturated soil. In other words, all sliding soil masses have the same weight density (denoted by ); (iv)ignoring the initial plastic strain of sliding soil mass before raining, at first, the initial effective stress parallel to the slope surface (denoted by ) is proportional to initial effective stress perpendicular to the slope surface (denoted by ) and is calculated using where is coefficient of lateral earth pressure. To elastic material, can be calculated by poisson ratio (denoted by ) using In additional, it is assumed that intensity rainfall has a constant value during the entire process. Based on the analytical solution of Iverson, porewater pressure with depth at time (denoted by ) can be calculated by where is osmotic coefficient of fully saturated soil, is Heaviside step function, is slope angle, is saturated hydraulic conductivity, and can be calculated using Function is given by where is complementary error function.
Otherwise, should meet the requirement of Hence, the factor of slope safety (defined by ) with depth at time can be calculated using where and are effective angles of internal friction and effective cohesion of soil, and is weight density of groundwater.
As shown in Figure 1, total normal stress (denoted by ) and total shear stress (denoted by ) with depth parallel to the slope surface are calculated using Porewater pressure with depth at time (denoted by ) can be calculated using (2.9). Effective normal stresses which are parallel to the slope surface (denoted by ) and are perpendicular to the slope surfaceare (denoted by ) are given by The relationship between Mohr’s Circle of Stress and shear strength on site is shown in Figure 2.
(a) Elastic state
(b) Plastic state
Variable is defined As shown in Figure 2(a), (2.12) is satisfied: Figure 2(a) shows that the increased porewater pressure after the rainfall will shift the Mohr’s Circle of Stress to the left till it is tangent to the curve of shear strength, and it becomes yield at this site, then inequation (2.12) will become equation. However, the yield surface due to yield point does not parallel the bottom plane of slope at this moment (Figure 2(b)); direction of principal stress axis will change with increased porewater pressure until become parallel to the bottom plane of slope. During this period, stress of plastic zone can be calculated using where θ is the angle of the maximum principal stress from yaxis, and it is given as Under plane strain conditions, MohrCoulomb yield criterion can be expressed as , calculated using Using the flow rule given by it, strain increment along the critical slip surface has a translational movement, also the value of it is zero, namely, ; the relations between stress and strain increment are taken to be as in which where is shear modulus and is Lame’s constants and calculated by using For , , and , the displacement of slop can be calculated by using The slope angle, depth of impervious surface, depth initial groundwater level, weight density of groundwater, the sliding soil mass weight density, coefficient of lateral earth pressure, effective angle of internal friction, saturated hydraulic conductivity, and shear modulus were set up to , , , kN/m^{3}, kN/m^{3}, , , m^{2}/s, and . Also, the coefficient of permeability is expressed as 10^{−4} m/s and the upper limit of light rain, moderate rain, heavy rain, heavy rainstorm, severe rainstorm is , 0.025, 0.05, 0.1, and 0.3, respectively.
Figure 3 shows the effect of rainfall in slope displacement and the factor of slope safety () under different rainfall amount. There is an approximate onetoone correspondence between displacement value and with various grades of rainfall; is found to decrease markedly with increasing displacement, and the influence of rainfall on the relationship between and displacement becomes slight (Figure 3(a)). It needs about 124 hours, 53 hours, 26 hours, 14 hours, and 5 hours to light rain, moderate rain, heavy rain, heavy rainstorm and severe rainstorm to make 10 mm displacement of slop (Figure 3(b)); so rainfall has a huge impact on the displacement value.
(a)  curve
(b)  curve
During the calculation of the , some parameters of model such as physical and mechanical parameters (effective angle of internal friction and effective cohesion of soil ) are uncertain, and they are closely related to test method, soil specimens, reduction factor, and so on. So it is difficult to forecast the displacement of highway slope by this way. Because displacement has an approximate onetoone correspondence to , displacement can be used to forecast the failure of highway slop. And rainfall has important influence on the slop displacement; forecasting the displacement of cohesionless soil highway slope in the future with monitoring data based on artificial neural networks (ANNs) should consider the rainfall.
3. Methods of Study Using Modular Neural Network
The most feature of ANNs is that they are based on a selforganizing structure that resembles the biological neural system of mammalian brains, and ANNs give an efficient method to solve problems that cannot be exact solved by statistical methods for their theoretical limitations. Thus, displacement forecasting models based on ANNs ignore influence of rainfall, and the values might be unobjective and unconvincing.
MNN uses divideandconquer technique, in which a problem is divided into a set of subproblems according to the inherent class relations among training data, then gives a network module to learn each of the subproblems, recombine the individual modules into a solution to the original problem at last. Before learning displacement monitoring data, sample set of slope displacement monitoring data should be correctly classified based on the rainfall grade, after that construct the subnetwork to make the subsample data can be precisely trained, and finally recombine the subnetwork outputs into one. The total output of MNN is given by where is the network output, is the output of the ith subnetwork (denoted by ) and is the normalized output vector elements of the subnetwork, and is the number of local networks.
Figure 4 shows the network architecture of the prediction system, which consists of three layers input layer, hidden layer, and output layer, and all of them are completely connected to form a hierarchical network. Subnetworks constitute the network, state of the switches is identified by rainfall, and is the number of local networks and switches. The sample set of slope displacement monitoring data is defined as , is the sample set of rainfall, and is the number of data. and have an approximate onetoone correspondence.
The displacementtime forecasting model is established as (3.2), which allows to forecast the th displacement of slope in the future: where is delay time step.
At first, construct eigenvector for rainfall , so that LR is light rain, MR is moderate rain, HR is heavy rain, TR is terrible rainstorm, and SR is severe rainstorm. Consider fuzzy classification of rainfall grade; the displacement monitoring data of sample set are classified according to calculation of eigenvector using where is the classification number of sample set, , is clustering center of , and is the membership matrix of sample to fuzzy subset. In (3.3), and [17]. Some data may simultaneously belong to sample subsets and . After that, the back propagation algorithm is used to train the gating and local network with sample subset , and set lower bound on the error and the maximum iteration number is 2000 [5]. Final results of the MNN can be calculated by (3.1).
To the new test data , the distance measure of and is chosen as criterion to decide which local network it belongs to, calculated using (3.4) and (3.5): in which, where is the relative distance measure of and , is the average distance measure of ith sample subset, and is the number of data for ith sample subset. Displacement monitoring data is trained only by local networks in which is the smallest one based on the traditional MNN it makes the calculation accuracy of boundary samples lower to others. To solve this problem, more than one local network are chosen to training with fuzzy decision, in order to significantly enhance the accuracy of MNN model. First, the relative distance measure of among all local network normalization can be calculated using In (3.6), and .
Defining fuzzy set of relative distance measure A = {very small (denoted by ), small (denoted by ), middling (denoted by ), large ()}, the values of reflect the distance measure between and . If , then local networks should be used to train for the distance of and is short. If , then local networks cannot be used to train for the distance of and is long.
The calculation of local networks selection involves the following steps.
Step 1. Calculate between every new test data and its corresponding local networks using (3.4), (3.5), and (3.6).
Step 2. Calculate membership degree of to every fuzzy subset A with the membership degree curves as Figure 5 shows.
Step 3. Make a choice which fuzzy subset belongs to by the highest membership degree in the fuzzy command, and local networks which belong to the same set are chosen to integrate from to 1. If the local network in which actually exist, then it should be used to train the data , otherwise, select local network in which , until selecting local network in which .
Due to the symmetry of the Gaussian distribution, the center value of membership degree decreases 50% at each time by adopting Gauss function for fuzzy membership functions. The center of fuzzy subset to membership degree is 0.25/P, 0.5/P, 1/P, 2/P and the value of membership degree is 1 when or . For , all cannot meet at the same time; meanwhile, the local network corresponding to would not be the empty set. The local network corresponding to should not be chosen for selecting sequence of local network is from to 1 in Step 3.
The weight of local network that is selected is given by where is the weight of local network that is selected, is the distance measure of ith local network, and is the number of local networks that has been selected.
4. Simulation Results
Hence, in order to illustrate the effectiveness of the model validity, one of Beijing city’s highway slopes is chosen for test. Geological investigations in the collapse area in this research are based on site investigation and the evaluation of remote survey images; the location and brief geological descriptions of slope are shown as Figure 6.
(a)
(b)
Time series analysis was used to set up displacementtime forecasting model for both MNN model and ANN model [18]. And then forecast cohesionless soil highway slope displacement in the future based on the same sample data.
Setting model parameters for the classification number of sample set, the number of hidden node in multilayer feedforward neural network that is used to train local networks is 4, and delay time step which is calculated with autocorrelation method is 6 [19].
Figure 7 shows the relative error for the group of displacement monitoring data that is in the same rainfall grade, namely, ignoring influence of rainfall through selecting the displacement monitoring data in one rainfall grade. And their respective margins of error are almost the same as each other because both of the methods use back propagation algorithm to training.
Displacement data records the state of this cohesionless soil highway slope from the period May 22 to October 31, 2011. Due to the fact that rainfall series is necessary for the failure of slope, the meteorological station was built at the slope. Figure 8 shows the grade of rainfall distribution map at precipitation station of highway construction area during the whole observation period.
Figure 9 shows three curves of forecast data based on two methods and the real data during the whole observation period. The forecast data consider the influence of rainfall.
Figure 10 shows their relative error curves of displacement forecast data.
Note that in Figures 9 and 10, MNN offers a highly satisfactory method to forecast the slope displacement, which is better than the ANN model. MNN model has a moderate relative error, but ANN model shows even bigger fluctuations especially after intense rainfall.
The value of mean square error and average relative error of these models under different conditions is just as shown Table 1.

According to the abovementioned analyses, we can get the results that (1) Whichever algorithm of artificial neural networks is chosen and both models have accurate forecast data of the slope displacement. (2) When considering the influence of rainfall, the fitting situation of MNN model is better than ANN model, mean square error of MNN model is only 7.73% of ANN model’s, and average relative error of MNN model is 27.65% of ANN model’s. (3) When considering the influence of rainfall, both models’ errors are sharply decreased, and MNN model performs better than ANN model.
5. Conclusion
Displacement monitoring data of highway slop is classified by the MNN based on the same type of data with the same character, which fits the monitoring needed. Also, the method reduces the complexity of the modular neural network and its generalization capabilities could remarkably be enhanced.
When considering the influence of rainfall, the precision of MNN model is much higher than ANN model. The MNN model has a good applicability for the selflearning ability and datadriven ability than ANN model, and it can enhance the applicability of model by retraining the networks with new data at regular intervals.
Acknowledgments
This research was supported by China’s Eleventh FiveYear National Science and Technology Support Project (no. 2009BAG13A02) and Science and Technology Project of MOT (no. 2012364223300).
References
 G. B. Crosta and F. Agliardi, “Failure forecast for large rock slides by surface displacement measurements,” Canadian Geotechnical Journal, vol. 40, no. 1, pp. 176–191, 2003. View at: Publisher Site  Google Scholar
 S. H. Ni, P. C. Lu, and C. H. Juang, “A fuzzy neural network approach to evaluation of slope failure potential,” Microcomputers in Civil Engineering, vol. 11, no. 1, pp. 59–66, 1996. View at: Publisher Site  Google Scholar
 S. Lee, J. H. Ryu, K. Min, and J. S. Won, “Landslide susceptibility analysis using GIS and artificial neural network,” Earth Surface Processes and Landforms, vol. 28, no. 12, pp. 1361–1376, 2003. View at: Publisher Site  Google Scholar
 K. M. Neaupane and S. H. Achet, “Use of backpropagation neural network for landslide monitoring: A case study in the higher Himalaya,” Engineering Geology, vol. 74, no. 34, pp. 213–226, 2004. View at: Publisher Site  Google Scholar
 H. B. Wang, W. Y. Xu, and R. C. Xu, “Slope stability evaluation using Back Propagation Neural Networks,” Engineering Geology, vol. 80, no. 34, pp. 302–315, 2005. View at: Publisher Site  Google Scholar
 L. Ermini, F. Catani, and N. Casagli, “Artificial Neural Networks applied to landslide susceptibility assessment,” Geomorphology, vol. 66, no. 1–4, pp. 327–343, 2005. View at: Publisher Site  Google Scholar
 H. Guo, W. Wang, W. Guo et al., “Reliability analysis of pedestrian safety crossing in urban traffic environment,” Safety Science, vol. 50, no. 4, pp. 968–973, 2012. View at: Publisher Site  Google Scholar
 J. B. Zhao, L. Li, and Q. Gao, “Research and application of the grey theory to slope deformation prediction,” Chinese Journal of Rock Mechanics and Engineering, vol. 24, no. 2, pp. 5799–5802, 2005. View at: Google Scholar
 D. Wang, M. Zhou, and Y. Li, “Displacement analysis of slope based on wavelet transformation and grey theory,” in Proceedings of the International Conference on Electric Technology and Civil Engineering (ICETCE '11), pp. 4929–4932, Lushan, China, April 2011. View at: Publisher Site  Google Scholar
 W. Wang, Vehicle's ManMachine Interaction Safety and Driver Assistance, China Communications Press, Beijing, China, 2012.
 H. Zhiquan and F. U. Wenbin, “The wholeregion method of chaotic time series for slope deformation prediction,” Journal of North China Institute of Water Conservancy and Hydroelectric Power, vol. 30, no. 3, pp. 70–73, 2009. View at: Google Scholar
 M. B. Jaksa, The influence of spatial variability on the geotechnical design properties of a stiff, overconsolidated clay [Ph.D. thesis], The University of Adelaide, Adelaide, Australia, 1995.
 T. Glade, M. Crozier, and P. Smith, “Applying probability determination to refine landslidetriggering rainfall thresholds using an empirical “Antecedent Daily Rainfall Model”,” Pure and Applied Geophysics, vol. 157, no. 68, pp. 1059–1079, 2000. View at: Publisher Site  Google Scholar
 J. M. Gasmo, H. Rahardjo, and E. C. Leong, “Infiltration effects on stability of a residual soil slope,” Computers and Geotechnics, vol. 26, no. 2, pp. 145–165, 2000. View at: Publisher Site  Google Scholar
 C. ChienYuan, C. TienChien, Y. FanChieh, and L. ShengChi, “Analysis of timevarying rainfall infiltration induced landslide,” Environmental Geology, vol. 48, no. 45, pp. 466–479, 2005. View at: Publisher Site  Google Scholar
 R. M. Iverson, “Landslide triggering by rain infiltration,” Water Resources Research, vol. 36, no. 7, pp. 1897–1910, 2000. View at: Publisher Site  Google Scholar
 M. G. C. A. Cimino, W. Pedrycz, B. Lazzerini, and F. Marcelloni, “Using multilayer perceptrons as receptive fields in the design of neural networks,” Neurocomputing, vol. 72, no. 10–12, pp. 2536–2548, 2009. View at: Publisher Site  Google Scholar
 X. Liu, X. H. Zeng, and C. Y. Liu, “Research on artificial neural networktime series analysis of slope nonlinear displacement,” Chinese Journal of Rock Mechanics and Engineering, vol. 24, no. 19, pp. 3499–3504, 2005. View at: Google Scholar
 M. T. Rosenstein, J. J. Collins, and C. J. De Luca, “Reconstruction expansion as a geometrybased framework for choosing proper delay times,” Physica D, vol. 73, no. 12, pp. 82–98, 1994. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2012 Yanyan Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.