Deformation Prediction of a Deep Foundation Pit Based on the Combination Model of Wavelet Transform and Gray BP Neural Network

)e purpose of this study was to predict the deformation of a deep foundation pit based on a combination model of wavelet transform and gray BP neural network. Using a case of a deep foundation pit, a combination model of wavelet transform and gray BP neural network was used to predict the deformation of the deep foundation pit. )e results show that compared with the traditional gray BP neural network model, the relative error of the combination model of wavelet transform and gray BP neural network was reduced by 2.38%. )is verified that the combined model has high accuracy and reliability in the prediction of foundation pit deformation and also conforms to the actual situation of the project. )e research results can provide a valuable reference for foundation pit deformation monitoring.


Introduction
At present, with the vigorous development of large-scale infrastructure in China, the excavation process of foundation pits will not only cause the deformation of foundation pits themselves, but also cause the deformation and displacement of adjacent buildings, rail transit, and underground comprehensive pipe corridors. In serious cases, it will not only cause huge economic losses, but also bring serious social consequences [1,2]. erefore, how to establish a suitable deformation prediction model to monitor the excavation deformation of the foundation pit and ensure the construction safety of the excavation process plays an important role. However, because there are many influencing factors of deformation in the process of foundation pit excavation construction, and since these factors have randomness and uncertainty, it is difficult for this to be expressed accurately with a simple mathematical model [3][4][5]. For this reason, many experts and scholars have put forward many prediction methods for the research of foundation pit excavation deformation prediction models.
Recently, the prediction methods for foundation pit deformation mainly include artificial neural networks, gray correlation degree methods, fuzzy comprehensive evaluations, Kalman filtering models, and hybrid combination models [6][7][8].
However, the above methods have achieved certain results in the deformation prediction of foundation pit, but there are also shortcomings. e convergence speed of artificial neural network is relatively slow, which makes local optimization easy, making the prediction result error larger; the gray correlation degree theory can easily reduce the prediction accuracy of the model when dealing with the situation of large data fluctuation; the fuzzy comprehensive evaluation model still has great difficulty in choosing which membership function. erefore, many scholars also try to use the hybrid combination prediction model to predict the deformation of foundation pits. Xu et al. [9] made use of the fact that the stability of the signal data decomposed by wavelet transform is better than that of the original signal and established the combination model of wavelet gray series and parallel series, respectively, to effectively predict the deformation of foundation pits. Liu et al. [10,11] used the combination of gray model and BP neural network to predict the deformation of deep foundation pits under the condition of few samples and insufficient information, but the final fitting results were too different from the measured results. Zheng et al. [12] established the prediction model of a wavelet BP neural network. rough the comparison with the prediction results of the BP neural network model, it is concluded that the prediction results of the wavelet BP neural network model are more accurate than those of the BP neural network model. e above models use the advantages of various deformation prediction models to predict the deformation of deep foundation pits using the hybrid combination model method [13,14], which improves the accuracy of prediction results. In fact, foundation pit deformation is based on timeseries prediction, including the trend term and the random term [15]. Among them, the trend term reflects the main law of foundation pit deformation; it is the main basis of deformation prediction and belongs to the nonstationary series. e random term belongs to noise series and has certain stability, which is one of the main reasons affecting the accuracy of the prediction method. erefore, these mixed combination model methods do not fully consider or distinguish the characteristics and influence of the trend term and the random term in the prediction process, and they have many shortcomings, such as being relatively simple and having their own characteristics and application occasions, so they cannot fully mine the original data information, and the prediction accuracy needs to be improved [16]. erefore, we tried to synthesize the advantages of three kinds of deformation prediction models using a wavelet denoising to remove and impute the outliers of the original signal data to obtain a more real and reliable data signal after denoising [17][18][19]. Using the GM (1, 1) model [20], which has a better effect in processing incomplete signals and few samples, and using the BP neural network with strong computing ability and the advantage of a good error correction ability [21,22], the prediction model of foundation pit deformation based on wavelet denoising gray BP neural network was established.
In this study, a new gray BP neural network method based on wavelet denoising is proposed. e prediction model method is widely used in the deformation prediction of deep foundation pit excavation. Firstly, the effective and reliable data of foundation pit deformation were obtained by wavelet denoising; secondly, the GM (1, 1) model was used to predict foundation pit deformation, and the predicted value was taken as the input sample value of BP neural network; and finally, the expected output prediction result of foundation pit deformation was obtained by learning and training. e contribution of this study is twofold: (1) e innovation of this study is to propose a new gray BP neural network method based on wavelet denoising. Compared with the traditional gray BP neural network model, this method can improve the accuracy and stability of deformation prediction of deep foundation pit, more accurately predict the future deformation of deep foundation pit, and be more in line with the engineering practice. (2) It can provide a useful reference for deformation monitoring and prediction of deep foundation pit excavation construction in the future and ensure the safety of deep foundation pit excavation construction.
e rest of this study is organized as follows. e related work is described in Section 2. e details of the proposed combined forecasting model method are shown in Section 3. e application case evaluation is presented in Section 4. e conclusion is drawn in Section 5.

Denoising Method of Deformation Signal Based on
Wavelet Transform. Foundation pit deformation monitoring data are a signal including noise, and its mathematical expression is as follows: where s (x) represents deformation monitoring data, which include real deformation and random noise (or system error); f (x) represents the real deformation, and e (i) is the random noise signal.
In the analysis of deformation monitoring data, to make the monitoring data more real, it is necessary to deal with the noise. Generally speaking, the noise signal is usually a highfrequency signal, so the real signal is a stable low-frequency signal.
e denoising principle of wavelet analysis is to reduce part of e (i), the noise signal, which is essentially the denoising of the signal or monitoring data s (x) so that the real deformation signal f (x) is separated from the signal s (x). In general, this method consists of three steps: wavelet decomposition, threshold quantization, and high-frequency coefficient reconstruction. In this study, the principle of wavelet multiresolution analysis was used for signal processing [23].
Suppose that space V 0 is approximated by space L 2 (R), the deformation signal f (x) is regarded as belonging to the approximate space a, and L 2 (R) is decomposed into orthogonal components V 1 and W 1 , then space V 1 is further decomposed into orthogonal components V 2 and W 2 , and so on; the approximate space V 0 of space L 2 (R) is decomposed into infinite orthogonal spaces: (2) where ⊕ is the sum of two approximate spaces of decomposition, and then, any signal f(t) can be decomposed into the sum of infinite orthogonal spaces, which can be expressed as follows: 2 Mathematical Problems in Engineering where scale function ϕ N, is the standard orthogonal basis of W j and L 2 (R). e whole energy of all finite signal spaces (L 2 (R) � ⊕ ∞ j�− ∞ W j ) is covered by the closure requirements of space W j , which is complete and satisfies the requirements of the Taking the two-level decomposition as an example, the structure of the signal wavelet multiresolution analysis tree is shown in Figure 1. In Figure 1, S represents the original signal, and A j and D j represent the decomposed low-frequency signal and high-frequency signal on scale j, respectively. en, the low-frequency signal is further decomposed into the approximation part of low frequency and high frequency. e relationship can be expressed as follows: It is worth noting that the different choices of wavelet function and threshold processing will lead to different signal denoising effects. At present, the Daubechies (dbN) wavelet and the Symlet wavelet are commonly used for signal denoising. e Symlet wavelet is based on the improvement of the dbN wavelet function. Its expression is as follows: e Symlet wavelet and the Daubechies wavelet have symmetry, which is suitable for image processing. Both wavelet functions can reduce phase distortion in signal analysis and reconstruction.
In this study, there are four commonly used threshold selection methods.
where N is the sampling length of the signal.

Min/Max reshold (Minimaxi).
Its principle is to minimize the maximum risk of the estimate, producing an extreme value of the minimum mean square error. Its threshold algorithm is as follows:

Unbiased Risk Estimation reshold (Rigrsure).
According to the risk curve Rish, the minimum risk point is denoted as k min , and then, the unbiased risk estimation threshold is defined as follows:

Heursure reshold (Heursure).
Heuristic threshold combines unbiased risk valuation threshold and fixed threshold, and it is defined as follows: If we can choose an appropriate threshold to denoise the signal containing noise by the threshold method, then we can achieve the purpose of denoising without losing truth.
To measure the effect of signal denoising, we adopted two quality evaluation methods. e first method is the rootmean-square error method (RMSE), which refers to the mean square error of the decomposed and reconstructed signal and the original signal. Its expression is as follows: where N is the signal length, F (i) is the original signal, and f(i) is the signal after wavelet denoising. RMSE represents the difference between the original signal and the denoised signal. In practical application, when the RMSE value is larger, the denoised signal has a low approximation to the original signal, and the denoising effect is poor; otherwise, the denoising effect is good. Another method is the signal-to-noise ratio (SNR) method, which is the ratio of signal value to noise value. It is generally believed that the higher the SNR, the better the denoising effect. It is defined as follows: Figure 1: A tree corresponding to multiresolution analysis.
where P S is the power of the original signal and P n is noise power.

e Gray Modeling
Method. e gray GM (1, 1) model is a common prediction model, which is composed of firstorder differential equations of first-order variables. Assume the original discrete sequence is as follows: A new sequence is generated by accumulating a sequence at a time: Further, the first-order differential equation is established for equation (13), and the gray GM (1, 1) model is as follows: where a and u are gray parameters, which can be solved by the least-squares method: where After solving a, the result is substituted into equation (14) to obtain: For cumulative generation, the prediction model is restored to: where x (0) (2), . . . , x (0) (n) is the signal data part for back testing; x (0) (n + 1), x (0) (n + 2), . . . is the data part of the predicted signal. Finally, the residual and relative residual are calculated as follows: e accuracy of the gray GM (1, 1) model can be evaluated by the posterior variance ratio and small error probability.
e posterior variance ratio and small error probability are, respectively: where S 1 and S 2 are the variances of model sequence x (0) (k) and residual sequence e (0) (k), respectively. e is the mean value of e (0) (k). e smaller the C, the more discrete the original sequence, and the larger the variance; the smaller the residual, the higher the accuracy of the simulation model. e model accuracy level can be expressed as follows: the level of Max {P, the level of C}, as shown in Table 1.

BP Neural Network Model.
A BP neural network is a multilayer feedforward network with error propagation (the sum of squares of network error reaches the minimum value). A BP neural network model is shown in Figure 2.
Assume the following: input node x i , hidden node y j , output node z k , network connection weight w ij of input node and hidden layer node, and network connection weight v jk of hidden layer node and output node. e expected value of the output node is t k . e output of hidden layer node is as follows: e calculated output of the output node is as follows: where net j � i w ij x i − Q j and net k � j v jk y j − Q k e error of the output node is as follows: If E meets the requirements, the sample meets the expected output; if E does not meet the requirements, the network connection weight needs to be corrected. where η is the training times, and α is the momentum factor. To improve the generalization ability of BP network, we adopted the early termination method. At the same time, the foundation pit monitoring data were divided into three subsets: the first subset was the training sample set, which was used to calculate the gradient and correct the weight and threshold of the network. e second subset was the validation sample set, which was used to minimize process fitting and monitor the error of the validation sample set during the training process. e third subset was the test sample set, which was used to test the final solution to confirm the actual predictive ability of the network.
is study studied the use of a three-layer BP neural network (input layer, hidden layer, and output layer composition structure) for foundation pit deformation prediction.
e specific selection of the parameters of the foundation pit deformation prediction model was to use the previous 12 monitoring data to construct training samples. If there were too many nodes in the input layer, the calculation cycle and workload would increase, and it would also have a certain impact on the prediction accuracy of the model. erefore, the first n-4 displacement values of the original data were used as the training sequence of the BP neural network, and the last four actual measurements were used as the predicted values (verification samples). Certain monitoring data of the designated measuring point were the output, and the eight monitoring data before the measuring point were the input. e eight input data and one output data of the measuring point constituted a training sample. In this way, 11 sets of training samples could be constructed using the previous 12 monitoring data.
e BP neural network model of foundation pit deformation prediction had 8 input variables and 1 output variable, and the hidden layer had a three-layer network model with 11 nodes. On the premise that the training accuracy could be guaranteed, it was enough to ensure that the hidden layer node was larger than the input layer node; otherwise, it would affect the training speed. In the calculation of the network topology, the tansig transfer function was used to build the network, the gradient descent training was used as the training function, the training accuracy was set to 0.0001, and the maximum number of training steps was 10 000. e training simulation calculation was performed to obtain the corresponding fitting value of the foundation pit deformation.

Proposed Combination Model of Wavelet Transform and Gray BP Neural Network
Based on wavelet analysis and gray BP neural network modeling method, we constructed a series combination model of wavelet transform and gray BP neural network. e series combination model uses the gray model to predict the data sequence after wavelet denoising and then takes the predicted value as the input sample value of BP neural network. After learning and training, the expected output prediction result is finally obtained. e specific modeling steps are as follows: (1) After preprocessing the original data sequence, the signal data sequence after wavelet denoising is obtained. (2) e processed signal data are judged to see whether they meet the conditions of gray modeling, and the nonconformed signal data series is transformed. (3) e signal data after denoising are predicted by the gray GM (1, 1) model. (4) According to the characteristics of the BP neural network, the prediction value of the gray GM (1, 1) model is taken as the input signal data sequence of the network. (5) By learning and training BP neural network, the predicted value of the expected output is finally obtained.
In the second point of the above steps, it is worth noting that the purpose of modeling was deformation prediction. To improve the simulation accuracy of post-prediction, it was necessary to delete old information and add new information consecutively; in particular, when the quantitative change accumulated to a certain degree and qualitative change occurred, we removed the old data that could no longer reflect the current characteristics of the system. erefore, from the perspective of prediction, the metabolic model was the most ideal prediction model.
After appropriate changes, the expression of the GM (1, 1) model can be written as follows:

(25)
Considering the accuracy of fitting and prediction, the applicable conditions of the model are as follows [24]: When − a < − 0.5, the GM (1, 1) model is used for medium-and long-term prediction. When 0.5 < − a − < 0.8, the GM (1, 1) model can generally be used for short-term forecasting. When 0.8 < − a < − 1, the GM (1, 1) model should use residual correction. When − a > 1, the GM (1, 1) model should not be used. e gray modeling feasibility test was performed on the original data series, generally including grade ratio test, smoothness, and smoothness ratio test.
Secondly, the smoothness and smoothness ratio of data series x (0) are tested; the condition for x (0) to be a smooth discrete function is as follows:( and ε is as small as possible. If the original data sequence x (0) can meet the conditions of the grade ratio test, smooth ratio, and smoothness test, it can be modeled in gray. If the above conditions are not met, the gray modeling conditions are not met, and the data need to be transformed. en, the rank ratio of the transformed sequence of the original data sequence x (0) must be within the coverage range.
In this study, the horizontal displacement data of foundation pit pressure were used as the input data, and the initial value of the original data sequence was used as the initial value of the cumulative sequence. Because the gray model needed to be trained first, and several measured data needed to be reserved to verify the accuracy of the model, the first n− 4 data of the original data were used as training samples, and the last four data were used to test the network. ey were then substituted into the already built gray GM (1, 1) model for training. e whole research process is summarized, and the technical route is given as follows (Figure 3).

Application Studies
In this section, a real-time dataset from the deformation measurement of a foundation pit is used to evaluate the effectiveness of the proposed method.

General Situation of Deep Foundation Pit Engineering and Description of Monitoring Data.
e basement of the project was two floors, the excavation depth of the foundation pit was about 10 m, and the area of the foundation pit was about 10408 m 2 . e foundation pit is located at main traffic arteries, and the excavation of the foundation pit has a safety impact on the surrounding buildings and roads. e site survey results of the foundation pit showed that the main soil layers within the influence range of the excavation depth of the foundation pit on the site were as follows: ① miscellaneous fill, with a thickness of 3.50-6.10 m; ② silt, with a thickness of 3.50-14.30 m; ③ silty soil, with a thickness of 1.70-5.30 m; ④ silt (including mud), with a thickness of 2.90-15.10 m; ⑤ silty clay, with a thickness of 1.20-3.60 m; and ⑥ coarse medium sand (including mud), with a thickness of 8.00-12.30 m. e hydrological conditions within the influence range of the excavation depth of the foundation pit were the following: the aquifer in the site was mainly the pore confined water in layers ③ and ⑤. e confined water level was about 3.73 m, and the permeability coefficient was K m � 10.23 m/d. e buried depth of the stable water level of the borehole measured during the survey was about 2.50 m. e supporting structure of foundation pit was a single row of SMW piles with two concrete supports. e first support bottom elevation was − 3.20 m, and the second support bottom elevation was − 7.7 m. e protection level of foundation pit was level I. In this project, 18 horizontal displacement observation points (S1-S18) were arranged along the edge of the support structure to measure the horizontal displacement (pile top). e layout of the measurement points is shown in Figure 4. e observation time was 152 days. e whole construction process of the foundation pit excavation was divided into the following four working conditions: (1) e first working condition: initial stress solution of undisturbed soil and initial displacement return to zero (30 days) (2) e second working condition: construction of pile foundation enclosure, excavation of the first layer of soil (the thickness of the soil layer was 1 m), construction of purlin, and the first concrete support (52 days) (3) e third working condition: construction of the second concrete support and excavation of the second soil layer (the thickness of the soil layer was 5 m) (40 days) (4) Fourth working condition: excavation of the third layer of soil (the thickness of the soil layer was 4 m) to the bottom of the foundation pit (30 days) Because the space was limited, and the middle part of the long side of the foundation pit was in a weak position, its monitoring data played a control role in the deformation of the support structure. erefore, to verify the effectiveness of the method in this study, the S18 point was selected as the research object.
e original records of horizontal displacement observation of pile foundation are shown in Table 2.
It is worth noting that to ensure the reliability of the actual monitoring data of the deformation of the foundation pit, in the model prediction and processing stage, the network model adopted a multistep rolling real-time prediction method, taking the current horizontal displacement of the pile top as the new input and replacing the horizontal displacement of the previous period; this method was used to continuously replace the input data, and it carried out a rolling dynamic forecast for the next time. In addition, the GM (1, 1) model uses the original certain information to predict and analyze unknown information. erefore, it is more sensitive to the original data. According to the new information principle, we know that the more newly acquired signals have the smallest noise error. At the same time, a single-sequence GM (1, 1) model can be constructed according to the time-accumulated deformation, and to meet the prediction requirements of the GM (1, 1) model, the data series must be positive. erefore, this study presents another expression for the horizontal deformation data measurement, which is expressed in terms of cumulative horizontal displacement and is beneficial in reducing the error; this is suitable for the original data sequence with a small number of calculations and a small sample, and the calculation accuracy can also meet the requirements. Mathematical Problems in Engineering

Horizontal Displacement Data (Pile Top) Based on
Wavelet Denoising. A different selection of wavelet function and threshold processing will lead to different denoising effects. In this study, the threshold value and the different parameters of wavelet basis function were selected to make a comparative study of the data effect after denoising, and the best model of wavelet denoising for the horizontal displacement of S18 point pile top was determined, and then, the next stage of analysis and prediction was carried out. For the evaluation of a different threshold selection, unified scale � SLN, soft threshold, Daubechies' 3 wavelets, and one-level decomposition were used. A fixed threshold (sqtwolog), unbiased risk estimation threshold (rigrsure), and heuristic threshold (heursure) were used to denoise the original data, and different filtering effects were obtained. e evaluation results of three different threshold selections are shown in Table 3.
According to the comparison of Table 3, in relation to the definition of root-mean-square error (RMSE) and the signalto-noise ratio (SNR), the smaller the RMSE value, the higher the approximation between the denoised signal and the original signal, and the better the denoising effect; the larger the signal-to-noise ratio, the better the denoising effect. From these two evaluation indexes, the denoising effect based on the unbiased risk estimation threshold (rigrsure) method is relatively better. erefore, to predict more accurately, we chose the principle of rigrsure for wavelet denoising.
We performed a comparison and evaluation of different wavelet basis functions. ere are many choices of wavelet basis functions, and the denoising effect of different wavelet basis functions is different. In this study, only Db6 and Sym6 wavelets with orthogonality and compactness were selected. We used the principle of the rigrsure threshold, scale � SLN, and one-level decomposition to denoise the original data; the evaluation results are shown in Table 4.
It can be seen from Table 4 that the root-mean-square error of the Db6 wavelet function was relatively smaller, and the signal-to-noise ratio was larger.
erefore, Db6 was chosen as the wavelet basis function of wavelet denoising. Finally, we decided to use the Db6 wavelet function, one-level decomposition, rigrsure function, and scale � SLN to denoise the original observation data and obtained the denoised data as shown in Table 2.

Predicting the Deformation Based on the Gray BP Network Model Based on the Measured Data.
According to the method introduced in the second part, the 8-step gray GM (1, 1) model was selected to predict the 9-12 horizontal displacement data using the measured horizontal displacement data of the pile top in the period of 1-8. e fitting value test and prediction value test of the 8-step gray GM (1, 1) model are shown in Table 5.
It can be seen from Table 5 that the posteriori difference ratio of the fitting value of the GM gray (1, 1) model was 0.31, less than or equal to 0.35, and its model accuracy grade was grade I; its prediction result was very good, and its reliability was high. In addition, it can be seen from the test results that, with the increase in observation periods, the relative error was also gradually increasing, mainly because the data change trend in this period suddenly slowed down, and the prediction error of the gray model itself also gradually increased with the increase in observation periods, which makes the prediction accuracy not high. erefore, the BP neural network prediction was carried out according to the GM (1, 1) fitting results. Considering the small amount of data, we took the GM (1, 1) fitting value of the pile top horizontal displacement data from the first period to the eighth period as the input sample, took the    Table 6. It can be seen from Table 6 that the accuracy of the gray BP combined model was improved to a certain extent compared with the single gray GM (1, 1), and the average relative error of the final model prediction was reduced by 1.47%, indicating that the prediction of the gray BP model is relatively reliable and effective.

Deformation Prediction of Gray BP Neural Network
Based on Wavelet Denoising. In the same way, according to the method in the previous section, considering the small amount of data, we used the GM (1, 1) fitting prediction result of the horizontal displacement of the pile top (wavelet denoising data) from the first period to the eighth period as the output value is used to predict the horizontal displacement of the BP neural network in periods 9-12, as shown in Table 7.
As can be seen from Table 7, with the increase in observation periods, the relative error had not simply increased. e relative error in the 10th and 11th periods was smaller and better, while the relative error in the 9th and 12th periods was larger. is shows that the deformation data after wavelet denoising can better reflect the real horizontal displacement of the pile top.

Comparison and Analysis of Gray BP Network Prediction Model before and after Wavelet
Denoising. Based on the above analysis results, the comparison between the measured data and the gray BP prediction data before and after wavelet denoising is shown in Table 8.
It can be seen from Table 8 that the prediction results of the gray BP model before and after wavelet denoising were larger than the actual observation data, which was due to a sudden change in the original observation data. From the comparison of the two groups of data, the average relative error of the prediction results of the gray BP model based on wavelet denoising was 4.32%, while the average relative error of the prediction results of the gray BP neural network model based on the original monitoring data was 6.70%. erefore, for the prediction data of the S18 pile top horizontal displacement, the gray BP model based on wavelet denoising was more accurate and reliable than the data without denoising.
In addition, it can be seen from Table 8 that with the increase in excavation depth of the foundation pit, the change in horizontal displacement increased. However,   during the predicted observation period from the 9th to the 12th (during the fourth stage of construction), although the horizontal displacement value of the pile top increased gradually, the increase was relatively gentle. is is because the horizontal displacement of the pile top tended to be stable with the completion of the excavation Earth pressure release and the formation of the second support. e supporting structure system formed by retaining the pile, crown beam, and support played a good role in controlling the horizontal deformation of the pile top. All the predicted results are also in line with engineering practice.

Conclusion
is research combines the advantages of wavelet transform, a gray model, and a BP neural network prediction model. We established a gray BP neural network prediction model based on wavelet transform and conducted deformation prediction research on deep foundation pit engineering. Studies have shown that the gray BP neural network model based on wavelet transform has an average relative error of 4.32% for the prediction results of deep foundation pit deformation, which is 2.38% less than the traditional gray BP neural network model. e prediction accuracy was improved, and the reliability of the method was verified. is method can provide a useful reference for deformation monitoring and the prediction of the deep foundation pit excavation construction in the future and ensure the safety of the deep foundation pit excavation construction process.
Due to the modeling granularity, our modeling only considered some key information, but some details were not considered, such as the influence of the settlement of the foundation around the foundation pit, the change in groundwater, and the change in the internal force of the concrete internal support. However, in fact, if these random factors are considered, the relative error of deformation prediction of foundation pits will be caused. erefore, in future work, we will study how to use other real-time monitoring data to predict the deformation of deep foundation pits. We will also apply other algorithms, including advanced learning algorithms, to search for the optimal solution more quickly. In addition, how to provide a clearer visual expression model structure to guide the safety of foundation pit excavation and construction will be another direction of future foundation pit deformation prediction research work.

Data Availability
e simulation data used to support the findings of this study are included within the article.

Conflicts of Interest
e author(s) declare(s) that there are no conflicts of interest regarding the publication of this study.