Data-Driven Deformation Reliability of Retaining Structures in Deep Excavations Considering Measurement Error

The deflections of retaining structures during deep excavation should be controlled below a predefined threshold in an attempt to mitigate the deformation-posed damage to the surrounding buildings and infrastructure facilities. In this paper, a new approach is presented to conducting deformation reliability assessments of retaining structures. The method is based on an extension of the classical first-passage reliability problem and is applied on the spatial scale of the retaining structures. With the proposed method, the in situ monitoring data and the associated measurement error can be easily incorporated, which accordingly improves the assessment accuracy. Moreover, the proposed method has a closed-form solution, which is beneficial for the computation efficiency and its practical application. The method was applied to assess the retaining structure reliability in a foundation excavation of a subway station in Fuzhou, China. The accuracy of the analytical results was verified through a comparison with those of Monte Carlo simulation. It was shown that the proposed method can represent well, the deformation-based safety level of retaining structures.


Introduction
The past decades have witnessed rapid growth in the number and scale of deep foundation pit projects in subway constructions around the world. Once a foundation pit accident (e.g., collapse of retaining structures) occurs, the consequence is often catastrophic, with ripple effects to the safety of adjacent buildings and underground pipelines [1][2][3][4]. In order to prevent such accidents, an important goal in the design and construction of foundation pits is to control the deformation of the soil and the retaining structures. Considering the uncertainties arising from the relevant parameters in the deformation process, the reliability approaches based on probability theory should be used to capture the randomness and variability of these parameters, yielding a probabilistic estimate of the safety level of the foundation excavations [5][6][7][8][9].
A significant number of studies are available in the literature regarding the deformation reliability assessment of retaining structures in foundation excavations. For example, Luo et al. [10] investigated the influence of statistical parameter uncertainty on the probability of serviceability failure in braced excavations based on the Bootstrap method. Wang et al. [11] conducted the probabilistic evaluation of serviceability failure in a braced excavation using the moment methods and the finite element platform Plaxis TM , and discussed the effect of spatial variability of soil parameters on the probability of serviceability failure. Goh and Kulhawy [12] calculated the reliability index of braced retaining walls by approximating the nonlinear limit state surface through a neural network model. Khosrojerdi et al. [13] evaluated the lateral deformation of geosynthetic reinforced soil walls and abutments based on some selected methods. Sayed et al. [14] studied the reliability of reinforced retaining walls based on finite element modeling. The response surface approach was utilized to represent the performance function, coupling the first order reliability method to compute the reliability index. Many of the existing studies have performed reliability assessments of retaining structures in foundation pits based on finite element software, where the computational cost could be dramatically high. Another approach is to establish a relationship between the retaining structure deformation and some key factors such as the wall deformation mode, wall stiffness, soil properties, adjacent building characteristics, and others [15][16][17][18]. However, these relationships are often empirically formulated and contain many unknown parameters to calibrate, which would become difficult in the presence of limited data and significant uncertainties. Alternatively, the direct use of in situ observed data could improve the accuracy and efficiency of structural reliability assessment [19][20][21][22][23][24][25]. For instance, Zhang et al. [22] investigated the impact of deep excavation on the adjacent metro shield tunnels using health monitoring data on the deformation. Qi and Zhou [24] developed a procedure for the Bayesian back-analysis of retaining wall reliability, making use of the deflection data monitored at multiple locations. However, the incorporation of measurement error (ME) associated with data collection [26,27] in the deformation reliability assessment of retaining structures has been limited discussed in previous studies. In fact, uncertain errors involved in the measurement of observations could deteriorate the quality of the data, and hence, the estimates of structural reliability [28][29][30].
This paper presents a method for deformation reliability analysis of retaining structures based on in situ monitoring data. The failure of a retaining structure is defined as a displacement at any depth exceeding the predefined threshold, which could differ from the engineering practice focusing on the deformations at the monitored depths only. By noting that the deformation distribution along the vertical direction of a retaining structure is a one-dimensional problem, this paper extends the classical first-passage reliability approach [31][32][33] from the temporal scale to the spatial scale. A closed-form solution is proposed for the deformation reliability assessment of retaining structures considering the potential failure risks at any depth. The impact of ME arising from data collection is explicitly incorporated in the proposed deformation reliability method. A case study is presented to demonstrate the applicability of the proposed method.

Deformation of Retaining Structures
The excavation of deep foundation pits in urban areas will inevitably result in the deformation of the soil and the retaining structures (e.g., retaining walls). Excessive deformation will trigger significant damage to adjacent buildings and infrastructure facilities. During the excavation process, the soil in the pit is removed and the earth pressure at the two sides of the retaining wall would be changed. The earth pressure on the outer side of the wall changes from static earth pressure to active earth pressure, and thus the difference in earth pressure between the inner and outer sides will lead to the horizontal deformation of the wall [34][35][36]. When the wall is supported at the top, as is often the case in engineering practice, the overall deformation of the wall along the vertical direction would behave in a "concave" way [37], as illustrated in Figure 1a (adopted from [38]).
A threshold for the lateral displacement is often used to determine the behaviour of the deformed retaining wall, which would be dependent on the excavation depth of the pit. When the deformation of the wall is monitored at some selected depths, the wall is deemed as satisfactory if all the monitored deformations do not exceed the predefined threshold, and vice versa. This method is relatively straightforward and reflects the current engineering practice for the performance evaluation of retaining structures in foundation pits, providing a simple yet useful tool for decision-makers. This approach only concerns the deformation at some certain depths (where monitored). In fact, the failure of retaining walls (i.e., the displacement exceeds the threshold value) could occur at any depth. Illustratively, as shown in Figure 1b, even if all the displacements at the monitoring points are smaller than the limit, the actual displacement curve may still exceed the threshold, leading to the failure of retaining wall. Thus, a reliability assessment approach is essentially needed moving from "point control" to "surface control". Conceptually, increasing the number of monitoring points and inclinometers will increase the accuracy of the traditional point control-based method. While this idea could improve the monitoring practice in engineering, it is often difficult to implement due to the significant constraints of costs and monitoring equipment. This fact explains the motivation of the proposed "surface control" method.
In this paper, the deformation reliability of retaining structures are considered, which is defined as the probability that the maximum displacement of the retaining structure within the monitored area does not exceed the displacement limit. Suppose that there are n pieces of monitored deformation data at different depths d 1 , d 2 , . . . d n ; namely, y 1 , y 2 , . . . y n respectively. Let y lim denote the predefined threshold for each y i . The existing approach ("point control") simply concludes that the structural behavior is satisfactory if y lim ≥ y i for ∀i = 1, 2, . . . n and unsatisfactory otherwise. Clearly, the possible occurrence of failure (displacement exceeding limit) at locations rather than {d 1 , d 2 , . . . d n } has been neglected. Also, the point control fails to take into account the ME associated with the monitored data. These disadvantages are addressed in this paper.

Mathematical Basis
The classical first-passage reliability method was originally developed for time-dependent reliability analysis based on the stochastic process theory. The basic idea of the method is to assess the probability that the structural response does not exceed the safety limit for any time point within a reference period of [0, T]. Rice [39] proposed the first-passage calculation formula, which was used in structural dynamic reliability analysis. Siegert [40] developed an approach for the first-crossing probability based on the continuous Markov process. Helstrom and Isley [41] derived the analytical solution of the first-passage time using the Markov's envelope process. Coleman [42] proposed a Poisson approximation approach to calculate the crossing frequency. The Poisson approximation establishes the relationship between the crossing rate and the structural time-dependent reliability.
However, the computational accuracy depends on the assumption that the structural response exceeds the safety limit independently with a small probability. Cramér [43] showed that the solution to time-dependent reliability based on the Poisson process assumption is asymptotically accurate when the safety limit tends to infinity. Later on, the first-crossing reliability approach was extensively applied in the estimate of structural time-dependent safety level [31,32,[44][45][46].
In this paper, the classical time-dependent reliability approach is extended, with some modifications, to estimate the deformation reliability of retaining structures, by noting that the spatial distribution of the deformation along the vertical direction is also a one-dimensional problem.
Mathematically, within a length of [0, X] of interest, the limit state function Z(x) takes the form of where x denotes the location, S(x) is the structural deformation at x, and R(x) is the deformation limit at x. If Z(x) ≥ 0 holds for all applicable x ∈ [0, X], then the retaining structure is deemed as satisfactory. Typically, R(x) can be treated as deterministic (or simply a constant in some occasions, e.g., y lim , as discussed in Section 2). Taking into account the uncertainty associated with the process S(x), the probability of Z(x) ≥ 0, ∀x ∈ [0, X] is the deformation reliability, denoted by where P( ) is the probability of the event in the bracket. Recall that the traditional "point control" approach (c.f. Section 2) simply gives where I( ) is an indicative function, which returns 1 if the statement in the bracket is true and 0 otherwise. Now we apply the first-passage reliability approach to estimate Equation (2). We assume that S(x) is a weakly stationary Gaussian process. For the purpose of further derivation, we let where E( ) denotes the mean value of the variable in the bracket. With this, Y(x) is a zero-mean stationary Gaussian process, as schematically shown in Figure 2. Thus, Equation (2) becomes The up-crossing rate of Y(x) with respect to Ω(x) at location x, ν + (x), is estimated according to whereẎ(x) andΩ(x) are the differentials of Y(x) and Ω(x) respectively. Rearranging Equation (6) yields Since Y(x) is a stationary Gaussian process, Y(x) andẎ(x) are mutually independent. Thus, it follows, where σ Y and σẎ are the standard deviations of Y(x) andẎ(x), respectively. Combining Equations (7) and (8) gives where Φ( ) is the cumulative density function (CDF) of a standard normal distribution. Assume that the upcrossings of Y(x) to Ω(x) are temporally independent and are rare, with which the Poisson process can be used to describe the occurrence of the upcrossings. Let N X be the number of upcrossings within the length [0, X], and it follows, By noting the fact that the structural deformation reliability within [0, X] equals the probability of N X = 0, one has where P(0) is the failure probability (the probability of deformation exceeding the limit) at x = 0. Specifically, as P(0) is typically small enough, it follows that ( [31][32][33]) A similar form of Equation (12) (with the terms X and x denoting time) has been extensively used for the estimate of structural time-dependent reliability [31,32,[44][45][46]. In this paper, Equation (12) is adopted with modifications to assess the deformation reliability of retaining structures.
Based on Equation (12), combining the concept of "point control", the deformation reliability, L nme (X), is formulated as follows, By comparing Equations (3) and (13), it can be seen that the proposed method in Equation (13) has moved the deformation limit state from "point control" to "surface control". Note that the subscript "nme" in Equation (13) means that L nme (X) has not (n) considered the measurement (m) error (e). The impact of ME on the deformation reliability is addressed later, in Section 4.3.

Parameter Calibration Using Observational Data
The deformation monitoring data is a set of discrete samples of a random sequence that varies with depth. Note that the relevant parameters (c.f. Ω(x) and Y(x) in Equation (13)) involved in the first-passage reliability problem are continuous along the depth x. This section presents a method for calibrating these quantities by regression analysis of the samples.

Fitting Ω (x) and Subsequently ObtainingΩ (x)
The retaining wall deformation caused by foundation pit excavation varies with the construction stage. The damage to adjacent buildings due to pit excavation at the early stages has been rarely reported but occurs mainly at the latter stages [5,47]. Therefore, the emphasis could be put on the deformation reliability assessment of retaining structures at the later stages of foundation pit excavation. In the presence of a sequence of observed deformation data, taking the depth x as the argument, the (continuous) deformation with respect to x can be fitted with an m-degree polynomial as follows, where γ i is the coefficient of the polynomial. In order to avoid Runge's phenomenon in the fitting [48] (the oscillation at the edge of the interval when using a high-order polynomial interpolation), the degree m should be carefully selected. Furthermore, one can obtain Ω (x) according to Equation (4), and furtherΩ (x) by differentiating Ω (x).

Fitting the Variance of Y(x)
With Equation (4) and a similar form of Equation (14), one can fit E(Y (x)) based on the discrete observed data y 1 , y 2 , . . . y n at depths d 1 , d 2 , . . . d n . Furthermore, employing the concept of "average variance" [49], the variance of Y(x) is computed as follows, where V ( ) represents the variance of the random variable in the bracket, and Y(d i ) is the deformation response evaluated at depth d i with the fitted function Y(x). The standard deviation of Y(x), σ Y , is then simply the square root of the variance.

Fitting the Variance ofẎ(x)
Based on the monitored deformation data y 1 , y 2 , . . . y n at depths d 1 , d 2 , . . . d n , the samples ofẎ are first approximated by With this, E(Ẏ (x)) is fitted using a polynomial expression similar to Equation (14). Subsequently, by referring to Equation (15), the variance ofẎ(x) is determined by where Ẏ (d i ) is the fitted functionẎ(x) evaluated at depth d i .

Test of Fitting Accuracy
In order to evaluate the accuracy of the fitting-based prediction of the quantities, the coefficient of determination (R-square) is a frequently-used tool [50], which is computed as follows, where y i is the fitted value, andȳ = 1 n ∑ n i=1 y i . This coefficient varies within [0, 1]. A value of R 2 approaching 1 indicates a well-fitted prediction model.

Measurement of Retaining Structure Deformation
An inclinometer system is often used to measure the deformation of retaining walls in practical engineering, which typically comprises a sensor probe, a load-bearing cable with a scaler, and a reader, as illustrated in Figure 3a,b. The sensor probe is a key component for the lateral displacement measurement in the inclinometer. A high-precision acceleration digital sensor is installed inside the probe. The tilt of the casing in four different directions, namely, A 0 , A 180 , B 0 and B 180 in Figure 3b, is measured by the probe. Two sets of guide wheels are located at the upper and lower parts of the probe respectively, with which the probe can slide up and down inside the inclinometer casing (c.f. Figure 3c).
The basic idea to measure the horizontal displacement with an inclinometer is to first determine the azimuth between the central axis of the inclinometer and the vertical line of the gravity pendulum in the probe, and then to use the geometric relationship (c.f. Figure 3d) to compute the accumulative displacement at different depths [51][52][53]. Let the inclination angle measured at a certain depth be α, and L the overstep (measurement interval), as shown in Figure 3d. With this, the horizontal displacement at this depth, δ, is computed by δ = L sin α according to the triangular relationship. Furthermore, calculating the horizontal offset from the bottom to the top, the lateral displacement at the jth measurement interval relative to the bottom, D j , is determined by where α i is the angle between the central axis of the inclinometer and the vertical line of the gravity pendulum at the ith measurement interval and δ i is the horizontal offset at the ith interval.

Sources of ME
The ME inevitably arises from the observed data in practice. In terms of the lateral displacement data monitored by an inclinometer, the ME could be caused by inclinometer accuracy and human error. The former consists of bias-shift error, sensitivity drift, rotation error, depth positioning error and others [51,55,56]. For example, the bias-shift error (also known as zero-shift error) is due to the initial checksum at the probe (c.f. A 0 and A 180 in Figure 3b) not being zero. Ideally, the sum of the readings of A 0 and A 180 should be zero with identical numbers but opposite signs. However, there is often a slight bias at the initial readings due to the imperfection of the probe. Details for other types of ME can be found elsewhere [51,55]. It was argued in [51] that the systematic error due to inclinometer accuracy, which accumulates arithmetically with the number of measurement intervals, is more important and significant than the random error (human error). As such, in this paper, only the ME due to inclinometer accuracy is considered. Taking into account its randomness by nature, the ME is represented by a random variable. The impact of ME on the deformation reliability of retaining structures is formulated in the next Section.

Deformation Reliability Incorporating ME
Recall Equation (13), which gives the deformation reliability of retaining structures without considering the potential impact of ME. Taking into account the systematic error in the measurement of lateral displacement, which dominates in ME [51], we let ∆ i denote the error associated with the ith (from the bottom to the top) measurement interval, which is a random variable representing the uncertainty of the measurement. Since the monitored data at different depths are obtained with the same inclinometer, and the ME accumulates with the number of measurement intervals, it is reasonable to assume that each ∆ i is fully correlated and follows a normal distribution. With this, it follows that where is a random variable that reflects the accuracy of the inclinometer, and k i is a scaler (i = 1, 2, . . . n) that reflects the accumulation of the ME with the number of measurement intervals. With Equation (19), where σ • denotes the standard deviation of • = ∆ i and . Let A denote the accuracy of the inclinometer. With a confidence level of 95%, it follows that Notice that the coefficient "2" in Equation (21) is explained by the 95% confidence level; that is, Combing Equations (20) and (21) gives Furthermore, in the presence of ME, the application of "point control" (c.f. Section 2 and Equation (3)) gives L pc (X) = P (y 1 + ∆ 1 ≤ y lim ∩ y 2 + ∆ 2 ≤ y lim ∩ . . . y n + ∆ n ≤ y lim ) = P (y 1 + k 1 ≤ y lim ∩ y 2 + k 2 ≤ y lim ∩ . . . y n + k n ≤ y lim ) Clearly, if y lim − y i ≥ 0 for ∀i = 1, 2, . . . n and A → 0 (i.e., infinite accuracy) in Equation (23), then L pc (X) → 1. On the other hand, if min(y lim − y i ) < 0 and A → 0, then L pc (X) → 0. These facts imply that Equation (23) is a generalized case of Equation (3). Furthermore, similar to Equation (13), the deformation reliability in the presence of ME is formulated as follows, The comparison between Equations (13) and (24) shows that the measurement accuracy of the inclinometer (A) has been explicitly incorporated in the deformation reliability of retaining structures. The applicability and accuracy of Equation (24) is demonstrated in the following.
Notice that the proposed method (Equation (13) or (24)) only needs the computational costs of calculating a one-fold integral, which enables the deformation reliability assessment to be performed efficiently compared with the traditional finite element-based approaches. Thus, the proposed method is promising to be adopted for use in practical engineering.

Background
The deep foundation pit project of a subway station considered in this Section is located in Changle District, Fuzhou, China. It is a two-story, double-column, and three-span station with a length of 293 m and a width of 21.7 m. The total construction area of the station is 17,861.2 m 2 . The station was constructed with the cut-and-cover method. The retaining structure is a 800 mm-thick continuous wall. The layout of the monitoring points is shown in Figure 4a, where the black dots labeled with "QCX" are the monitored locations. Figure 4b presents a photo from the in-situ construction, where the reinforced concrete supports are used at the top of the excavation and steel supports are used elsewhere.
The CX-06B inclinometer was used in this project to measure the deformation of retaining walls. The CX-06B inclinometer is a geotechnical instrument equipped with a high-precision accelerometer. It has been widely used in the measurements of lateral displacements for building foundations, mines, foundation pit excavation, and underground structures. The major technical features of CX-06B are as summarized in Table 1.

Deformation Reliability Assessment
We illustratively chose the monitored data at location QCX04 (c.f. Figure 4a) to perform deformation reliability assessment. The horizontal displacements measured at different depths on each day are shown in Table 2. Note that the depths in Table 2 were calculated from the bottom, which would equal 43.5 m at the top (ground surface).
The fitting of E(S(x)) is performed using a polynomial, with which Ω(x) and Y(x) in Equations (13) and (24) can be subsequently determined. The fitting accuracy is evaluated by the determination coefficient R 2 in Equation (17). Using the monitoring data from the subset of 1/4-3/4 of the depth, the dependence of R 2 on the order of polynomial (c.f. the term m in Equation (14)) is shown in Figure 5. It can be seen that when m ≥ 2, R 2 > 0.9, fairly close to 1. While a greater value of m will lead to a greater value of R 2 , the Runge's phenomenon could occur as m is sufficiently large. Thus, m was chosen as 4 for the following analyses. Furthermore, assuming that the displacement limit R (x) is 73.5 mm, Ω (x) can be determined according to Equation (4). The derivative of Ω(x) giveṡ Ω(x). The variances of both Y(x) andẎ(x) were computed with Equations (15) and (16). With this, the deformation reliability was found to be L nme (43.5) = 0.956 according to Equation (13). That is, the probability that the maximum lateral displacement of the retaining wall at location QCX04 does not exceed the predefined limit is 95.6%.
Next, the impact of ME on the deformation reliability is discussed. As introduced in Table 1, the inclinometer CX-06B used in this project has a measurement accuracy of A = 4 mm/15 m = 2.667 × 10 −4 . With this, applying Equation (23) yields L pc (43.5) = 0.518. Furthermore, with Equation (24), the deformation reliability taking into account the impact of ME is estimated as L me (43.5) = 0.956 × 0.518 = 0.495. That is, in the presence of ME, the probability that the maximum lateral displacement of the retaining wall at location QCX04 does not exceed the predefined limit decreases to 49.5%. This implies that ignoring the impact of ME could significantly overestimate the deformation reliability of retaining structures.    (14)).

Verification of the Proposed Method
In order to verify the accuracy of the proposed method in Equation (24), the Monte Carlo simulation method is used. The analytical and simulated reliabilities are expected to agree well with each other to demonstrate the applicability of Equation (24). The basic idea of the simulation-based approach is to generate a set of random lateral displacements distributed along the vertical direction of the retaining wall, and compare each of them with the predefined threshold to determine the structural state (failure or survival). For the jth simulation run (j = 1, 2, . . . q, where q is sufficiently large), the procedure is summarized as follows.
1. Let r j = 1, which becomes 0 if the jth sampled scenario of lateral deformation leads to structural failure (i.e., maximum lateral displacement exceeds the threshold). 2. Calculate the displacement distribution along the vertical direction of the retaining wall at an interval of 5 mm (excluding the monitoring points) using Equation (14). 3. Generate a set of correlated zero-mean normal random variables with a standard deviation of σ Y using the Nataf transformation method [32], and add them to the sample values from Step 2 respectively. This is to obtain a set of displacement values that are approximately continuous along the depth and have the same statistical characteristics as those used in the analytical approach. 4. Sample a realization of normally distributed with a mean value of 0 and a standard deviation of σ . Use Equation (22) to compute k i and subsequently utilize Equation (21) to calculate each ∆ i , where i = 1, 2, . . . n. 5. The sampled data at depth d i is calculated as y i + ∆ i for i = 1, 2, . . . n. 6. Combining those sampled data in Steps 3 and 5, for each of them, if the displacement exceeds the predefined limit, then set r j = 0. 7. Record r j .
With q simulation runs, the simulated deformation reliability is approximated by 1 q ∑ q j=1 r j . The analytical and simulation-based reliabilities with different displacement thresholds are presented in Figure 6 for the purposes of comparison. The number of simulation runs is q = 10 6 for each case. It can be seen that at the upper tail (e.g., ≥ 73 mm), the results from both approaches agree well with each other, which indicates the accuracy of the proposed method in Equation (24). At the lower tail, the reliability associated with the analytical approach is slightly greater than that of the simulation-based approach due to the impact of Poisson assumption in the analytical formulation (c.f. Section 3.1). Nonetheless, as the emphasis is on the behavior at the upper tail for most engineering cases, the proposed reliability assessment method in Equation (24) can be used to reasonably estimate the deformation safety level of retaining structures.

Discussions
The deformation reliability approach in Equation (13) based on the crossing rate model takes into account the potential failures at the depths rather than the monitored ones, and thus yields a reliability method that moves from "point control" to "surface control". Using the monitored data for location QCX04 in Table 2, the dependence of deformation reliability on the displacement limit is presented in Figure 7. The maximum monitored displacement is 73.35 mm according to Table 2. Correspondingly, in Figure 7, the deformation reliability equals 0 if the displacement limit is smaller than 73.35 mm. On the other hand, when R(x) ≥ 73.35 mm, the deformation reliability increases from 95.6% to 1. This means that the conventional "point control" may overestimate the deformation reliability by simply focusing on the deformations at limited points.
Furthermore, taking into account the impact of ME, the deformation reliability is also shown in Figure 7 for the purpose of comparison. Clearly, the dependence of L me in Equation (24) on the displacement threshold is smoother compared with L nme . For cases where L nme > 0 (that is, min(y lim − y i ) ≥ 0), it can be seen that L me < L nme . This can be explained by noting that in Equations (13) and (24), 2 n min i=1 y lim −y i Ad i < ∏ n i=1 I(y i ≥ y lim ) = 1, which gives L me < L nme . This comparison again suggests the important role of ME in the deformation reliability of retaining walls.

Conclusions
In this study, a new approach was developed for the deformation reliability analysis of retaining structures in deep excavations, making use of the in situ monitored data. The method expands the classical first-passage based time-dependent reliability problem into the spatial scale, and explicitly incorporates the potential impact of measurement error. The new method considers the deformation scenario of the whole structure rather than limited points, and thus moves from the conventional "point control" method to "surface control". The proposed method only involves the calculation of a one-fold integral, with which the deformation reliability assessment can be conducted efficiently compared with the traditional finite element-based approaches. An illustrative example is presented in this paper to demonstrate the applicability of the proposed method. The following conclusions can be drawn from this paper.
1. The accuracy and applicability of the proposed deformation reliability method were verified by comparing the results with those obtained from a Monte Carlo simulation. 2. If only taking into account the deformations at the depths monitored, the deformation reliability of retaining structures could be overestimated, due to the potential failures at the non-monitored points. 3. If ignoring the impact of ME, the deformation reliability could be overestimated significantly if the deformation scenario passes the "point control" examination. This suggests the importance of carefully identifying the measurement error of the equipment (inclinometer) in an attempt to achieve a reliable estimate of the safety level of retaining structures.
The proposed method can be further extended for the deformation reliability assessment of excavation surfaces containing multiple retaining structures. Further work is underway regarding the deformation reliability assessment of deep excavations on both the spatial and the temporal scales.