A Prediction Model of Structural Settlement Based on EMD-SVR-WNN

,


Introduction
With the rapid development of information technology, automatic monitoring of structural deformation has become an important way to ensure structural safety [1][2][3]. Based on the monitoring data, the prediction of structural deformation to prevent potential safety hazards during construction and operation of large buildings (such as subway, tunnel, and bridge) has become a research hotspot of structural deformation monitoring in recent years [4]. As the scale of the project increases, higher requirements are put forward on the construction and operation safety of the structure. e processing of the real-time monitoring in the field of structural deformation monitoring also becomes especially important. According to the accurate and effective analysis of the real-time structural monitoring data, the structural deformation prediction model can be effectively established, which is of great significance to ensure the structural safety.
At present, the prediction models in the field of structural deformation prediction can be roughly divided into single prediction models and combined prediction models. e single prediction methods such as regression analysis method, time series analysis [5], grey system theory [6], and artificial intelligence method [7][8][9][10] are commonly used in the structural deformation prediction [11]. Li et al. [12] applied the GM (1, 1) grey system theory model to the prediction of subgrade deformation and explained the application conditions of the theory with engineering examples. However, the grey system requires the data sequence to have nonnegative characteristics, and the data accumulation has the rule of grey index, which limits the wide application of Grey eory in structural deformation prediction. Zhenget al. [13] proposed the applicability principle of predicting the deformation of highway high slope based on stepwise regression analysis and verified the proposed method is an effective prediction method. Time series analysis requires the data to be linear and stable, and the deformation monitoring data in practical engineering are complex and nonlinear, which will affect the prediction accuracy. Chen and Wei [14] introduced fuzzy time series into deformation analysis and prediction. Experimental results show that fuzzy time series can be effectively applied to deformation prediction, and the prediction accuracy is improved. Liu and Hao [15] used the wavelet neural network to establish the prediction model. e result shows that this method has the great nonlinear function approaching capability. However, there are some shortcomings in the applicability and generalization of the single model. In view of the problems of a single model, some improved models [16,17] and combined models have achieved better results.
In order to further improve the prediction accuracy of structural deformation monitoring data, researchers have proposed a hybrid model that combines two or more prediction models, which can correspondingly optimize the processing of nonstationary time series signals. Liu et al. [18] proposed a Least-Square Support Vector Machine (LS-SVM) model based on the phase space reconstruction for forecasting nonlinear time series of dam deformation, and the experimental results indicate that the forecasting performance of the proposed method is significantly superior to that of the traditional multiple regression method. Jingzhou et al. [19] established the Kalman-ARIMA-GARCH (Generalized Autoregressive Conditional Heteroskedasticity) model to predict bridge structural deformation, and the results showed that the modified model can effectively improve the prediction accuracy. Su et al. [20] established a deformation prediction model by combining support vector machine, phase space reconstruction, wavelet analysis, and Particle Swarm Optimization (PSO). e results show that the method has high modeling efficiency and prediction accuracy. Jiang et al. [21] proposed the Least-Squares-Support Vector Machine-Markov Chain (LS-SVM-MC) model, which was based on the Markov Chain predicted by Least-Square Support Vector Machine (LS-SVM) to correct the error, and applied the model to the deformation analysis of an arch dam. e results show that the model has high prediction accuracy. Zhao et al. [22] proposed an enhanced empirical wavelet transform (MSCEWT) based on the maximum-minimum length curve method. is method is used to decompose the signal into a series of intrinsic mode functions (IMFs). Cao et al. [23] proposed a combined prediction model for dam deformation prediction. e model uses the EMD method to decompose the original signal and then uses the correlation vector machine and ARIMA model to predict the decomposed components, respectively. Ren et al. [24] used the EMD method to decompose the original deformation monitoring data into subcomponents with different characteristics and then used GA algorithm to optimize the parameters of the WNN model to predict the values of different components, and finally fuse the predicted values of each subsequence. In summary, when predicting deformation monitoring data, decomposing the signal into internal components with different characteristics is of great significance to the improvement of prediction accuracy.
Aiming at the poor applicability of a single model and the characteristics of deformation data, this paper proposes a combined prediction model based on Empirical Mode Decomposition, Support Vector Regression, and Wavelet Neural Network (EMD-SVR-WNN). Firstly, the EMD model is used to decompose data into different components with different physical characteristics. EMD decomposition can effectively decompose monitoring data into trend and random fluctuation term according to energy distribution. Secondly, the Support Vector Regression (SVR) prediction model [25] can overcome the shortcomings caused by highdimensional model construction in dealing with nonlinear problems.
irdly, the Wavelet Neural Network (WNN) prediction model [26,27] is used to predict the trend term. e WNN prediction model combines the advantages of wavelet transform and neural network model, and it combines good time-frequency localization characteristic and self-learning functions of neural network theory, which makes it have the best ability of function approximation and fault tolerance and overcome the shortcomings of slow learning speed and low network training success rate in traditional neural networks. Finally, the prediction results of the SVR model and WNN model are fused to obtain the final model prediction results. e results show that the EMD-SVR-WNN model proposed in this paper can achieve good results and has high prediction accuracy, and it is an effective structural deformation prediction model, which can provide important decision support for structural safety prevention. e rest of this paper is organized as follows. Section 2 reviews the related work in the prediction of structural settlement and gives a brief introduction to the related methodology.
e EMD-SVR-WNN prediction model is introduced in Section 3. Section 4 discusses the prediction results and performance analysis. Section 5 makes the conclusion of this paper and gives some recommendations.

Methodology
2.1. Empirical Mode Decomposition. Empirical Mode Decomposition (EMD) is an adaptive signal decomposition method, which can adaptively decompose nonlinear and nonstationary time series into trend and random terms, including several Intrinsic Mode functions (IMF) with the same characteristic scale and a residual component. IMF must satisfy two conditions at the same time when decomposing [28][29][30]. Firstly, the total number of extreme points is equal to or at most one different point from the total number of zero crossing points in all time periods of decomposing signals. Secondly, in any time period, the mean value of the upper envelope composed of local maximum points and the lower envelope composed of local minimum points is 0. For any structural deformation sequence, the specific decomposition procedure is as follows: (1) Find the extreme point of the original signal x(n), and all the extracted local maximum and minimum points were interpolated by cubic spline 2 Advances in Civil Engineering interpolation to obtain the upper (U max (n)) and lower (V min (n)) envelopes. (2) Calculate the mean value of the upper and lower envelope, and equation (1) is obtained: (1) (3) Calculate the difference between the signal and the mean value of the envelope, and equation (2) is obtained: (4) Determine whether T(n) satisfies the above two IMF conditions. If the condition is met, take T(n) as the first component of x(n), i.e., c 1 . If the above condition is not met, T(n) will be taken as the new original signal to continue the iteration until the above conditions are met. (5) Calculate the value of the residual component, it can be noted as follows: (6) Take the residual component r 1 (n) as the original signal and repeat steps (1) to (5). en, continuously decompose the 1, 2..., m IMF components that satisfy the condition, and the iteration stops until the residual component is a monotonic function. It can be noted as follows: 2.2. Support Vector Regression. Support Vector Machine (SVM) is a pattern recognition method based on statistical theory [31][32][33], which was first used to solve the classification problem in pattern recognition. Support Vector Regression (SVR) is the application of support vector in the field of function regression. SVR is a typical nonlinear forward feed network with hidden units, which can realize regression prediction of time series. e structure settlement data is time series data with random, nonstationary, and complex characteristics. e SVR model can overcome the defects caused by the construction of high-dimensional model. And it also has a solid statistical theory, and fewer data samples can be used to train data [34]. erefore, the SVR model is widely used to deal with nonlinear and nonstationary time series regression prediction problems. e main idea of support vector regression is to map data to a high-dimensional feature space through nonlinear mapping and then conduct linear regression on the data in this high-dimensional space [35]. Given a set of data points (x 1 , y 1 ), (x 2 , y 2 ), ..., (x m , y m ) for prediction model, where x i ∈ R n , y i ∈ R n , and m is the total number of the sampling points. A linear regression function can be noted as follows: where ω is weight vector and b is bias value; the ω and b can be obtained by minimizing the following function: where ‖ω‖ 2 is the expression that controls the function capacity, C is the penalty coefficient, ε is the insensitive loss function [33], and ξ i , ξ * i are, respectively, the upper and lower training error subject to the ε. e optimization problem of the above can be expressed as the following nonlinear regression function: where α i and β i are Lagrange operators and e kernel functions commonly used at present include polynomial kernel function, radial basis kernel function, and sigmoid kernel function.
Different kernel functions have different advantages; choosing the appropriate kernel function is very important for constructing nonlinear regression prediction models. In the field of machine learning, the radial basis kernel function can effectively realize the nonlinear mapping from input space to high-dimensional feature space. e convergence domain of the radial basis kernel function is wide and the performance is stable, and it is suitable for arbitrary distributed samples. erefore, combined with the characteristics of structural settlement deformation data, the radial basis kernel function is selected as the kernel function of the support vector regression model in this paper.

Wavelet Neural Network.
Wavelet neural network (WNN) is based on the topology of BP neural network, the Wavelet basis function is selected for the transmission signals of hidden layer nodes in the Neural Network. WNN combines wavelet transform with neural network, and the advantages of them are fully utilized. e structure diagram of the wavelet neural network model is shown in Figure 1, where x n is the input parameter of the wavelet neural network, y m is the output parameter of the wavelet neural network, and g i (x) is the wavelet basis function. When the signal h(x) is input in the wavelet neural network, the wavelet basis function g is used for fitting, and equation (8) is obtained:

Advances in Civil Engineering
where h(x) is the fitted function, N is the number of wavelet bases, w n is the weight of input layer and hidden layer, a n is the scale factor of the wavelet basis function, and b n is the translation factor of wavelet basis function. e prediction process of the wavelet neural network model can be divided into two parts: training process and prediction process. e training process of the network is to find the parameters that make the network converge. e prediction process is to use the parameters obtained by network training to obtain the predicted value of the whole network. e training process of wavelet neural network is similar to the neural network method, which combines the forward propagation of signals with the backpropagation of errors. Wavelet neural network has strong learning ability and function approximation ability, and it overcomes the shortcomings of traditional neural network such as slow learning speed and low success rate of network training [36,37]. In view of these characteristics of wavelet neural network, this paper selects the wavelet neural network model to predict the trend term of monitoring data.

EMD-SVR-WNN Prediction
Model. e EMD-SVR-WNN prediction model proposed in this paper is based on empirical mode decomposition, support vector regression, and wavelet neural network model. e EMD decomposition model decomposes structural deformation data into different components with different physical characteristics and reduces the nonstationary characteristic of deformation monitoring data to a large extent. en, different prediction models are constructed according to the internal characteristics of different components. SVR prediction model overcomes the shortcomings of neural network model, such as overfitting and easy to fall into local optimum. It is considered as an alternative method of artificial neural network [38] and has great advantages in dealing with nonstationary and nonlinear time series. erefore, the SVR model is selected to predict the random term of the decomposed monitoring data. WNN prediction model combines the advantages of wavelet transform and BP neural network and has strong self-learning ability and localization of wavelet transform [39]. erefore, the WNN model is adopted to predict the trend term of decomposed monitoring data.
In actual projects, monitoring data are generally obtained by automatic monitoring equipment. It is assumed that the original deformation monitoring data can be expressed as a two-dimensional matrix X, where m is the number of sensors and n is the number of sampling points for each sensor; then, the original data matrix can be expressed as X ∈ R n×m . Suppose the kth column of the original data is selected for analysis, which is recorded as en, it is decomposed by EMD: where p represents the number of Intrinsic Mode Functions (IMF) of different feature scales and r n represents the residual component. In the process of signal decomposition, the frequency and energy distribution of the original signal will change. In order to further determine the trend term and random term of the signal, analyze the energy changes of different IMF components and extract the features of each different IMF component. e specific energy characteristic matrix solution steps are as follows: (1) Firstly, the energy of p IMF components is calculated according to (2) en, the energy of each IMF component is normalized according to (3) e energy characteristic matrix can be expressed as follows: Assuming that the energy feature matrix values of the IMF components of the front f layer are less than the energy threshold ε, the IMF components of the front f layer are used as the random term X r , and the remaining IMF components are used as the trend term X t : x 2 x n Input layer Hidden layer Output layer … … … Figure 1: Wavelet neural network model structure. 4 Advances in Civil Engineering SVR and WNN models are used to predict the random term X r and trend term X t , respectively. e final model prediction result can be obtained by fusing the prediction result X r of the random term and the prediction result X t of the trend term: If d m is the target output value, the mapping relationship between target input X ⌢ and target output Y ⌢ is R z ⟶ R, and the data matrix X ⌢ and Y ⌢ is as follows: In formula (15), nis the length of the data sequence and z is the embedding dimension. e specific steps and flow charts of the EMD-SVR-WNN prediction model proposed in this paper are shown in Figure 2. e implementation process of the method is as follows: Step1: by performing EMD decomposition on the structural settlement monitoring data, low-frequency components and high-frequency components are obtained. e value of the energy characteristic matrix is used to determine whether the IMF component belongs to trend X t or random data X r , according to formulas (10)- (13).
Step2: the random term X r uses the SVR prediction model to predict; select the first z data samples in the sample as the training set to predict n − z data, and the random term prediction result is recorded as X r .
Step 3: the trend term X t uses the WNN prediction model to predict; select the first z data samples in the sample as the training set to predict n − z data, and the trend term prediction result is recorded as X t .
Step 4: the prediction results of random and trend term are fused to obtain the final prediction results according to formula (14).

Performance Evaluation Criteria.
e evaluation criteria to measure the error of the prediction results included root mean square error (RMSE), mean absolute error (MAE), mean square percentage error (MSPE), and mean absolute percentage error (MAPE). e RMSE and MAPE are selected as evaluation criteria of EMD-SVR-WNN model in this paper. e formulas are as follows: where y i is the real value of the structural deformation, y i is the predicted value of the structural deformation, and N is the number of predicted samples.

Prediction Results and Performance Analysis
To evaluate the effectiveness of the proposed model, the experiment was carried out using the ground settlement data of a subway station in TianJin, and the settlement data contained 1000 sampling points. It can be seen from Figure 3 that the monitoring data reflects the actual characteristics of cumulative vertical displacement of the surface. In this paper, the cumulative variation of surface settlement of 800 sampling points in a section of a construction subway is selected for experimental testing. e first 500 data are selected as training samples to predict the last 300 data.

Analysis of EMD-SVR-WNN Model
Results. In order to illustrate the necessity of EMD decomposition, the amplitude-frequency characteristics of the settlement monitoring data to be tested are analyzed firstly, and the Fast Fourier Transform method is used to transform the time domain into the frequency domain. It can be seen from Figure 4 that the settlement monitoring data is mainly concentrated in the low-frequency part. is phenomenon shows that the data N Y Structural settlement monitoring data EMD decomposition  In order to further analyze the characteristics of the settlement monitoring data, the power spectral density of the data is analyzed, as shown in Figure 5. rough the analysis of Figure 5, this shows that the set of settlement data can be decomposed into multiple components with different frequencies. At the same time, because the set of settlement data changes slowly and is nonlinear and unstable, selecting a suitable decomposition method can reduce the instability of the signal and extract the characteristic information of the signals at different frequencies. Data prediction based on the abovementioned signal decomposition can greatly improve the prediction accuracy.
EMD decomposition method is suitable for processing nonlinear and nonstationary signals. e settlement signals are decomposed into high-frequency components and lowfrequency components by EMD decomposition, and then different components are predicted by appropriate prediction methods, so as to reduce the impact of the instability of settlement data on the prediction results. e trend and random terms of data are obtained by EMD decomposition. e decomposition results are shown in Figure 6. It can be seen from Figure 6 that the EMD decomposition divides the data into 5 layers, of which the first few layers show obvious high-frequency characteristics, and the higher the level of the data, the lower the frequency. At the same time, the data after EMD decomposition becomes more stable and smooth, and the data is adaptively divided into trend items and random items. is paper extracts the energy features of the abovementioned 5-layer IMF components and normalizes the energy values of different IMF components to construct an energy feature vector matrix. e constructed energy feature vector is shown in Table 1. It can be seen from the energy characteristic values in Table 1 that the energy value of the 5th layer IMF component accounts for the largest proportion and can be regarded as the main energy of the signal. From the perspective of signal decomposition, the trend term, that is, the low-frequency component contains the main energy of the signal, and the high-frequency component accounts for the secondary energy of the signal. erefore, in this article, the IMF5 component is used as the low-frequency component, and the remaining IMF components are superimposed as the highfrequency component of the data sequence. e result is shown in Figures 7 and 8.
After the EMD decomposition of the deformation monitoring data, the data is divided into four IMF components and one residual IMF5. e components after EMD decomposition show more obvious physical characteristics. e analysis and prediction of settlement deformation data can achieve better results on the basis of the EMD decomposition. In this paper, the first four IMF components are firstly predicted by the SVR model. After superimposing the first four layers components, the first 500 data are selected as training samples, while the last 300 data are selected as test samples, and the predicted results are compared with the real values.
e predicted results of SVR model are shown in Figure 9.
It can be seen from Figure 9 that the predicted value of the SVR model is very close to the real value of the  monitoring data, and only the prediction results at the extreme points have a certain deviation. e maximum absolute error between the predicted value and the real value is 0.3090 mm, and the minimum absolute error between the predicted value and real value is 0.00007 mm. e result shows that the trend of the predicted deformation value is basically consistent with the real value. After the prediction of the random term, the trend term is predicted by the WNN model. e first 500 data of the fifth layer component are selected as training samples to predict the last 300 data. e predicted results of the WNN model are shown in Figure 10.
It can be seen from Figure 10 that the predicted value of the WNN model is almost identical with the real monitoring value, and the prediction accuracy is high. e maximum absolute error between the predicted value and the real value is 0.0898 mm, and the minimum absolute error between the predicted value and the real value is 0.00003 mm.
Finally, the prediction results of SVR and WNN are fused to obtain the final prediction results of the EMD-SVR-WNN model. e final prediction results are shown in Figure 11.
It can be seen from Figure 11 that the prediction results of the proposed model are very close to real values, and the prediction accuracy is very high. It can be concluded from the experimental results that the RMSE of the proposed   Advances in Civil Engineering model is 0.0531 mm and the precision reaches 99.19%, which indicates that the proposed model is an effective settlement prediction model. e result also indicates that the EMD-SVR-WNN model can be well applied to the field of structural deformation prediction and achieve higher prediction accuracy.

Performance Analysis.
In order to evaluate the effectiveness of our proposed approach, the SVR, WNN, and DBN-SVR models are selected to compare with the EMD-SVR-WNN model proposed in this paper. In the SVR model, the nuclear function is configured as "rbf," the number of iterations is 10,000, and the penalty factor is taken as 0.01. In the WNN model, the number of iterations is 1000, the number of the hidden layer nodes is 3. In the DBN-SVR model, the number of network layers in the DBN model is set as 3, the number of iterations is 200, the kernel function of the SVR classifier is "rbf," the number of iterations is 10,000 and the penalty factor is 0.01. e computer configuration for the experiment is shown Table 2.
Under the experimental environment and the parameter conditions in this paper, the performance comparison of different prediction models is shown in Table 3 and the comparison of different prediction models with real values is shown in Figure 12.
It can be seen from Figure 12 that the SVR model has the worst prediction performance; the DBN-SVR and WNN model are better than the SVR model, while they show weakness compared with the EMD-SVR-WNN model. e prediction value of the proposed model is almost coincided with the monitoring data.
e experimental results show that the proposed EMD-SVR-WNN prediction model has better prediction performance, and it is an effective model for structural deformation prediction.
It can be seen from Table 3 that the prediction accuracy of the EMD-SVR-WNN model, respectively, improves 87.10%, 81.46%, and 63.35% compared with the SVR, WNN,    Advances in Civil Engineering and DBN-SVR models. e average prediction accuracy is 77.30% higher than other prediction models. Computation complexity is also an important criterion for model evaluation. It usually refers to the resources required at runtime after being written into an executable program. e higher the complexity of the model, the longer the running time; running time is used to discuss computation complexity in this paper. DBN-SVR and WNN models have the longest and shortest running time, respectively, but the accuracy is low. e proposed model has a relatively short running time while ensuring prediction accuracy. e experimental results show that the proposed EMD-SVR-WNN prediction model has higher prediction accuracy than other prediction models, and it is an effective structural deformation prediction model.

Conclusions and Recommendations
Structural deformation is inevitable during the construction process, so it is of great significance to avoid accidents through structural deformation prediction. Most of the existing structural deformation prediction models lack the mining of internal features within the data, and there is a lack of discussion on the applicability of the prediction model. Aiming at the above problems, this paper proposes a combined prediction model based on Empirical Mode Decomposition, Support Vector Regression, and Wavelet Neural Network (EMD-SVR-WNN). Based on the analysis of the measured data at Tianjin Station, the following conclusions can be drawn: (1) e data is decomposed into different components with different physical characteristics through the EMD decomposition method. According to the different characteristics of different signal components, the appropriate prediction model is selected, which can effectively improve the accuracy of deformation prediction. (2) According to the characteristics of random items and trend items, the SVR model has great advantages in dealing with nonlinear problem and can overcome the defects caused by high-dimensional model construction. e WNN prediction model combines the advantages of wavelet transform and neural network model. SVR and WNN are finally selected for prediction, respectively. (3) Finally, the prediction results of the SVR model and WNN model are fused to obtain the final model prediction results. Combined with practical engineering, the results show that the EMD-SVR-WNN model proposed in this paper can achieve good results and has high prediction accuracy; this is of great significance to the study of structural deformation in engineering.
Data Availability e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.   Advances in Civil Engineering 9