Collapse Risk Analysis of Deep Foundation Pits in Metro Stations Using a Fuzzy Bayesian Network and a Fuzzy AHP

Collapse risk analysis is of great significance for ensuring construction safety in foundation pits. )is study proposes a comprehensive methodology for dynamic risk analysis of foundation pit collapse during construction based on a fuzzy Bayesian network (FBN) and a fuzzy analytical hierarchy process (FAHP). Firstly, the potential risk factors contributing to foundation pit collapse are identified based on the results of statistical analysis of foundation pit collapse cases, expert inquiry, and fault tree analysis. )en, a FAHP and improved expert elicitation considering a confidence index are adopted to elicit the probability parameters of the BN. On this basis, quantitative risk reasoning and sensitivity analysis of foundation pit collapse are achieved by means of fuzzy Bayesian inference. Finally, an actual deep foundation pit in a metro station was used to illustrate a specific application of this approach, and the results were in accordance with the field observations and numerical simulation results. )e proposed approach can provide effective decision-making support for planners and engineers, which is vital to the prevention and control of the occurrence of the foundation pit collapse accidents.


Introduction
With the rapid development of urbanization around the world, metros have become a promising solution for relieving overground traffic in congested urban areas. However, metro construction is exposed to large potential risks due to various potential risk events in an uncertain environment. erefore, the number of accidents is increasing in metro engineering, especially in developing countries such as China. For example, on November 15, 2008, 21 people were killed in a collapse caused by the diaphragm wall failures in the Xianghu Station of Hangzhou Metro Line 1 [1]. Additionally, an accident that occurred in the Haizhu Square Station of Guangzhou Metro Line 2 caused five people to be trapped, and nearby buildings were cracked as a result of a diaphragm wall collapse [2]. As we can see, the frequent occurrence of foundation pit collapse accidents in metro construction has led to serious consequences and aroused the attention of all parties, including the public [3].
In recent decades, various risk assessment approaches developed based on probabilistic risk analysis (PRA) have been widely used in construction projects to avoid financial losses and personnel casualties. For example, these approaches include Monte Carlo simulations (MCS) [4,5], event tree analysis (ETA) [6,7], fault tree analysis (FTA) [8,9], and decision trees (DT) [10]. e above-mentioned methods have played a vital role in improving the risk management and control ability of large construction projects [11,12]; however, in many circumstances, these methods are incapable of coping well with uncertainties and giving satisfactory results on account of incompleteness or shortage of data [13]. Furthermore, it is often difficult or sometimes even impossible to obtain statistical data for PRA modeling, and this approach often relies on expert knowledge and experience. In recent years, the fuzzy reasoning methods [14][15][16][17] and fuzzy hybrid methods [18][19][20] have been developed based on fuzzy set theory (FST) [21] to overcome the limitations of PRA methods. Although the fuzzy reasoning and fuzzy hybrid methods can solve engineering problems under uncertainty, these methods suffer from limitations in conducting inference inversely and dynamic risk analysis. Feed-forward-like approximate reasoning approaches are strictly one way; that is, when a model is given a set of inputs, it can predict the output, but not vice versa [22][23][24]. is may have limitations on the flexibility of a safety analysis and assessment method that focuses on exploring causal relationships among risk factors. Furthermore, the aforementioned methods have a limitation in effectively achieving dynamic risk analysis and cannot achieve real-time safety control. When related parameters such as hydrological and geological parameters are changed, the above-mentioned approaches cannot accurately handle the updated features of dynamic process in project construction. erefore, effective real-time safety measures cannot be implemented with a change in environment [25].
Recently, Bayesian networks (BNs) have been developed to model the complexity of man-machine systems [26] and are widely used in many fields to perform uncertain knowledge representation and probabilistic reasoning because of the ability to represent multistate variables, dependencies among variables, and update probabilities [25,27,28]. Sousa and Einstein [29] presented a BN-based risk analysis approach for decision support in tunnel construction. Liang et al. [30] used BNs to conduct risk assessment of debris flow hazards. Khakzad et al. [31] performed a quantitative risk assessment of drilling operations based on a bow-tie and BN method. Zhou et al. [32] conducted quantitative risk analysis in foundation pit construction project based on the BN and field data.
Although BN has many advantages, it also has the limitation of using a probability measure to assess uncertainty. Traditional BNs can achieve quantitative risk analysis if exact information is sufficient and available [33], which is nearly impossible owing to lack of detailed data and incomplete knowledge, especially for construction projects [34,35]. As mentioned earlier, FST provides a powerful tool to deal with uncertainties and ambiguities associated with engineering problems [36]. erefore, the combination of fuzzy logical and BN theory may provide an alternative means to carry out risk assessment under an uncertainty environment [37]. A fuzzy BN-based safety risk analysis model was proposed for risk assessment of tunnel-induced pipeline damage by Zhang et al. [37]. Wang and Chen [38] introduced a systematic risk assessment method for metro construction using a fuzzy comprehensive BN. Similarly, Sun et al. [39] utilized a multistate fuzzy Bayesian network (FBN) to conduct risk analysis of tunnel collapse. Compared with traditional BN, FBN seems more flexible and interpretable [22]. Moreover, conventional BN can only handle precise probabilities, which in most cases are difficult to estimate, especially for large and complex projects. As FBN incorporates the advantages of both FST and BN, it can deal more effectively with uncertainty of risk data and describe the dependent and complex relationships between risk events [23,40,41]. Nevertheless, the elicitation of prior probabilities and conditional probability tables (CPTs) under data shortage or incompleteness is one of the key challenges to be solved.
In the above-mentioned studies, the probability parameters are usually determined via expert elicitation to overcome the incompleteness or lack of data, and linguistic expressions and integral value methods are usually adopted to quantify prior probabilities. However, the involvement of human judgement inevitably brings subjectivity and vagueness and thus causes the subsequent probability analysis to be unreliable and untrustworthy. In addition, the elicitation of large number of probability parameters of the BN not only puts great workload on the experts but also poses challenges to the quality or consistency of the elicited result [42]. Additionally, most of the existing researches focus mainly on risk analysis of tunnel construction, whereas few cases of foundation pit construction are reported. However, the foundation pit construction may face more potential risks than tunnel construction due to its more complex surroundings.
In construction practice, various risk factors including internal (geological, hydrological, etc.) and external ones (construction, organization, etc.) are able to lead to foundation pit collapse.
us, a systematic and integrated method is required to carry out collapse risk analysis of foundation pit in construction, and this is of great significance for ensuring construction safety. A comprehensive methodology for quantitative risk analysis is presented in this study to compensate for the shortcomings of existing models. Firstly, the potential risk factors that contribute to foundation pit collapse are identified via FTA. ereafter, a fault tree (FT) model is transformed into a corresponding BN model to conduct quantitative risk reasoning. en, the FAHP technique is adopted in this study to determine the prior probabilities and CPTs of nodes in BN, as it can overcome subjectivity and ambiguity issues and is an effective approach for handling uncertainty associated with expert evaluations. Moreover, a decomposition method is used to eliminate the number of comparisons as well as reducing possibility of inconsistency. Finally, the sensitivity performance measure is applied for sensitivity analysis to determine the critical risk factors for foundation pit collapse. e remainder of this paper is structured as follows. Section 2 presents the developed methodology for risk analysis of foundation pit collapse based on a FBN and a FAHP. A case study of an actual deep foundation pit in Qianhu Avenue Station is described in Section 3. Section 4 discusses issues to be further studied, and conclusions are given in Section 5.

The Developed Methodology
To decrease the occurrence probability of collapse risk and ensure safe construction of foundation pit, a comprehensive methodology for risk analysis of foundation pit collapse is developed based on a FBN and a FAHP. e following five steps are included in the proposed approach: (1) risk identification and FT construction based on expert knowledge and collected information; (2) FBN model establishment; (3) prior probability elicitation of the root nodes based on the FAHP and improved expert elicitation; (4) CPT elicitation of the intermediate nodes and leaf node using the FAHP, decomposition method, and improved expert elicitation; and (5) FBN-based risk reasoning and sensitivity analysis. e overall procedure for risk analysis of foundation pit collapse based on a FBN and a FAHP is shown in Figure 1.

Fuzzy Bayesian Network.
In traditional BN, the failure probability of nodes is usually treated as a precise value, but insufficient data often makes it difficult or even impossible to obtain crisp probabilities, especially for construction projects [25]. In construction engineering fields, the occurrence probabilities of basic events are usually more imprecise and vague due to deviation of measured values, construction operation errors, engineering experience, and other factors. As FST aims to cope with uncertainties caused by imprecision and fuzziness, a BN model integrated with FST is developed, and this combination is referred to as a FBN. Moreover, many studies have shown that a FBN not only can produce similar results with the same model and similar input data but also is more flexible and interpretable than a traditional BN and is the inheritance and development of a traditional BN [22,23].

Risk Identification and FT Construction.
Risk identification involves carrying out risk mechanism analysis for accidents and documenting their characteristics. In this paper, we analyze the risk mechanism of foundation pit collapse from internal factors (hydrogeological conditions) and external factors (design, construction, management factors, etc.) and reveal the causal relationships between risk factors via FTA. In addition, domain experts are also invited to query risk factors according to the construction site conditions and collected material (the literature, construction standards, and technical manuals). en, the final FT model of foundation pit collapse can be built by combining the results of FTA and expert inquiry.

Mapping the FT into a Bayesian
Network. BN construction is usually complex and is the "bottleneck" of its further application in safety risk analysis. e method of BN construction based on FT transformation is considered an effective approach to solve the above problem. Hence, BN construction based on the method of mapping the FT into a BN is adopted in the present study. Mapping from FT into the BN is based on graphical and numerical functions, and the specific mapping procedure is shown in Figure 2 [27]. In graphical mapping, the primary events, intermediate events, and the top event of the FT are directly transformed into root nodes, intermediate nodes, and the leaf node in an equivalent BN, respectively [28]. Overlapping nodes are merged into one [40]. In numerical mapping, the occurrence probabilities of the primary events in the FT are assigned to the relevant root nodes of the BN as a priori probabilities. CPTs are developed based on the Boolean gates of the FT.

Prior Probabilities and CPT Establishment.
Generally, the prior probabilities and CPTs of BN can be obtained based on a shrinkage approach [43] or hill climbing algorithms [44] if sufficient available data are collected. However, it is very difficult to present enough available data in engineering practice for BN model construction and parameterization due to the rare occurrence of foundation pit collapse and the limits of various learning algorithms. erefore, experts' domain knowledge and experience remain important sources for BN modeling, and the establishment of probability parameters of BN still relies on expert elicitation. However, the involvement of human judgement inevitably brings subjectivity and vagueness.
us, in this study, a FAHP was employed to handle the uncertainty and ambiguity of criteria and the judgement process. Moreover, an improved expert elicitation method based on a confidence index was used to guarantee survey data reliability, which takes both objectivity and subjectivity into account. An expert's judgement ability is determined according to his/her designation, educational level, and service time. e judgement ability level represented by ξ is determined based on the four indicators including expert's designation, educational level, service time, and age, and the corresponding expert judgement ability level is computed using a FAHP method. e system of expert information is illustrated in Figure 3. e expert profiles and judgement ability level calculated via the FAHP are shown in Table 1. e higher the value of ξ is, the more reliable the expert judgement will be.
In addition, the subjective reliability level represented by ψ is utilized to measure the reliability of experts' judgements on their own and is also divided into five levels, represented by "0.6, 0.7, 0.8, 0.9, and 1.0." e higher the value of ψ is, the more certain the judgement will be. When N experts are involved, the confidence index of the kth expert, indicated by θ k , can be calculated by (1). If the expert cannot make a judgement about his/her subjective reliability level or give a subjective reliability level ψ less than 0.6, the data will not be adopted.
e processes of combining FAHP, decomposition method, and improved expert elicitation to establish prior probabilities and CPTs of BN are shown in the following steps.

Generation of Root Nodes' Prior
Probabilities. In such circumstances, assume that there are n states (S 1 , S 2 , . . ., S n ) of a node N, and the probability of N at state S i is denoted by P(S i ). Generally speaking, P(S i ) is usually elicited by experts according to their experience, and it may be feasible while the number of nodes' states is not large. However, estimating probabilities to all states at a time is difficult when the number of states is large and the BN model structure is complex [45]. Moreover, it may inevitably bring uncertainty and bias during expert elicitation. As mentioned earlier, the fuzzy AHP, a decomposition method, and improved expert elicitation are adopted in this paper to elicit probability parameters of BN as follows.
Each expert individually implement a pairwise comparison according to Saaty's 1-9 scale [46], and A (k) � [a (k) ij ] n×n represents the pairwise comparison matrix obtained from the kth expert. Furthermore, each expert should also give his/her subjectivity reliability level ψ k towards his/her own judgements. en, the fuzzy positive reciprocal matrix A � [a ij ] n×n can be obtained by integrating M experts' grades using (2), where a ij is a TFN with    where M is the number of experts, θ k denotes the confidence index of the kth expert and can be obtained by (1), and θ k / M k�1 θ k represents the degree of reliability regarding the judgement data of the kth expert.
In addition, the consistency of fuzzy matrix A should be tested before the calculation of weight [47]. A new judgement matrix A must be obtained again when the above test results are inconsistent until the consistency of matrix A is verified to be acceptable.
After the consistency test, the fuzzy weight matrix can be calculated by (3)-(5). First, the geometric means of the TFNs for the ith state of node N can be obtained as By summing up the g i values, we have en, the fuzzy weight w i related to the prior probability of the ith state of node N can be computed as follows: Since the fuzzy weight w i obtained above is a fuzzy number, the "α-weighted valuation defuzzification technique" proposed by Detyniecki and Yager [48] is selected here to convert fuzzy results into crisp values. Defuzzification of a TFN w � (w l i , w m i , w u i ) is shown in where Val(w i ) denotes the transformed exact value, represents the α-weighted valuation function, and Average(F α ) is the average of the α-level set and can be obtained as follows: where u α and v α can be determined with the help of Figure 4.
Generally, we set f(α) � 1 and α � 0.5, and then, we can obtain the value of Val(w i ), as shown in (10).
Finally, P(S i ) can be obtained via normalization of w i as

CPT Establishment for a Node with One Parent.
Assuming that node N and its only parent M have n states (S N1 , S N2 , . . ., S Nn ) and m states (S M1 , S M2 , . . ., S Mm ), respectively, our aim is to derive the probabilities of each state of node N when the parent M is at state M j , i.e., P(S Ni | S Mj ) (i � 1, 2, . . ., n; j � 1, 2, . . ., m). We can judge which one (the states S Ni and S Nη (η �1, 2, . . ., n)) is more likely to occur when given the state of S Mj by inquiring experts [45]. And each expert only needs to fill in the number represented by a ij shown in Table 2. us, the corresponding comparison matrix is constructed as in the table. us, p(S Ni | S Mj ) � w ij can be obtained from the calculation of the above matrix, and the calculation procedure is the same as the generation of root nodes' prior probabilities.
Since node M has m states, m matrices should be constructed to obtain all w ij (i � 1, 2, . . ., n; j � 1, 2, . . ., m) values, which are the CPTs of the node N. erefore, we can establish the CPTs for a single-parent node based on the above procedure.

(12)
As we can see, there are a large number of state combinations of its parent nodes, and it is difficult for experts to directly estimate the above probabilities under such situations. Hence, a decomposition method is adopted here to reduce the elicitation workload by simplifying CPT elicitation [49,50]. For example, when a node A has three parents B, C, and D, its conditional probability on B, C, and D can be approximated by where α is a normalization factor to ensure that a∈A P(a | B, C, D) � 1. us, (12) can be computed as follows: where α is a normalizing constant to ensure that rough the above technique, we can see that the problem of CPT calculation for a node with multiple-parent nodes is transformed to compute the CPTs for a singleparent node. Additionally, estimating p(N � S Ni | M j � S M j P j ) individually will be much easier than directly esti-   (16). is approach can predict the potential risk of top event T and provide decision-making support for foundation pit construction.
where t � {t 1 , t 2 , . . ., t P } is states range for the leaf node T, stands for the JPD of the root nodes.

Sensitivity Analysis.
In engineering practice, decision makers usually pay more attention to which factors play an important role in foundation pit collapse. In addition, the critical checkpoints are usually determined based on expert knowledge and practical experience. Sensitivity analysis provides an alternative method that observes the relevant changes in the probability distribution of the top event by changing the probability distribution of each factor. e sensitivity performance measure (SPM) proposed by Zhang et al. [37] is applied for sensitivity analysis in this paper and can be calculated by (17). e higher the SPM (X i ) of root node x i is, the greater the contribution x i is for the risk sensitivity of top event T.
where t and x i refer to the state of the top event T and risk factor X i , respectively.  Figure 5). In addition, there are four sewage pipelines with a diameter of 800 mm on both sides of the road, only 20 m to the edge of the foundation pit. According to the investigation report, the main strata of the site comprise miscellaneous fill (1-2), silty clay (2-1), medium-coarse sand (2-4-5), round gravel (2-7), and medium weathered argillaceous siltstone (5-1-2), as shown in Figure 6. e groundwater in the station includes perched pore water and bedrock fissure water, and there is no confined water. Figure 7 shows the working site of the Qianhu Avenue Station construction. e excavation depth of foundation pit is 16 m, and an open-cut method was used in the station construction as shown in Figure 7(a). e retaining structure of the foundation pit consists of diaphragm walls and three rows of inner supports (Figure 7(b)). e three rows of supports are used to restrain the diaphragm wall deflection. e first row of support was concrete support (800 mm × 1000 mm) and the second and third rows of supports were steel supports (Φ609 mm). In addition, the excavation principle of "layering, blocking, symmetry, and balance" is adopted during the foundation pit excavation, as can be seen in Figure 6.

Risk Identification and Fault Tree Construction.
e metro construction accident data from 2004 to 2017 in China were collected through the network, accident investigation reports, and a literature review [51] and analyzed by statistical methods. We found that foundation pit collapse is the main type of accident, accounting for 44% of the total number of accidents. Moreover, foundation pit collapse is considered a large potential risk to personnel safety and surrounding facilities. erefore, foundation pit collapse risk is chosen as the top event (T) in the FT for the present study.
Due to complexities and uncertainties of the foundation pit construction process, the foundation pit collapse may be influenced by a large number of risk factors that interact with each other. erefore, a database of 48 deep foundation pit collapse cases of metro stations is analyzed to provide support for risk identification. Moreover, the results of expert inquiry and the method of FT are also used to determine the risk factors and reveal the causal relationships of risk factors. Finally, a FT diagram of foundation pit collapse for the Qianhu Avenue Station is constructed based on the above analysis. e three types of influential factors that contribute to T are revealed and analyzed in detail as follows: (1) Hydrogeological conditions (y 1 ): e geological condition (y 4 ) and hydrologic condition (y 5 ) are intrinsic risk surroundings; thus, five primary influential factors are selected, namely, thickness of soft stratum or permeable sand layer (x 1 ), accuracy/degree of site investigation (x 2 ), groundwater level (x 3 ), rainfall (x 4 ), and drainage timeliness (x 5 ) on the basis of foundation pit collapse case analyses, as shown in Table 3.      advantages of time and space saving and safety, and they can be used as waterproof curtains during construction and permanent basement walls after construction. And the inner supports are being indispensable part of the diaphragm walls as foundation pits have become deep and large, because they can effectively restrain the diaphragm wall deflection. erefore, the rationality of diaphragm wall and support design (y 6 and y 7 ) has great impact on the stability and safety of the pit. e diaphragm wall thickness (x 6 ), diaphragm wall insertion ratio (H O / H) (x 7 ), number of supports (x 8 ), and support form (x 9 ) are selected based on the cases analysis, as shown in Table 3. (3) Construction and management factors (y 3 ): e construction quality and management level play a crucial role in the pit stability, and improper construction organization and management may lead to potential risks during pit construction. ey can be further classified into three types: reliability of excavation (y 8 ), reliability of support (y 9 ), and field monitoring and management (y 10 ). en, the eight primary risk factors (x 10 , x 11 , . . ., x 17 ) related to construction and management are selected according to the cases analysis, as shown in Table 3.

Bayesian Network Modeling.
e BN model (shown in Figure 8) for collapse risk analysis of metro pit construction can be established in light of the transformation process described in Section 3.2. Each node of BN and relevant descriptions are shown in Table 3. All the nodes are considered to be multistate variables to make the results more accurate. us, the corresponding risk state of each node is divided into three levels defined by "L (low), M (medium), and H (high)," as illustrated in Table 3.

Prior Probabilities and CPT Establishment.
First, detailed information about the project's geological exploration, design, construction, and field investigation should be collected. en, domain experts are invited to express their opinions towards each node of the FBN. Clemen and Winkler [52] indicated that 3 to 5 experts are enough to achieve the most of the knowledge.
us, we invited five domain experts to elicit prior probabilities and CPTs of the FBN based on the collected information and their own experiences. e root node X 4 (rainfall) is taken as an example. e probabilities of X 4 in each state can be elicited from experts through the procedure introduced in Section 2.4.1. Each expert individually gives a pairwise comparison matrix A (k) � [a (k) ij ] 3×3 in the light of Table 2 and a subjectivity reliability level ψ k towards all states of the root node X 4 , and the fuzzy positive  3 Construction and management factors Good 0-100 Moderate 60-80 Poor 0-60 y 4 Geological condition Good 80-100 Moderate 60-80 Poor 0-60 y 5 Hydrologic condition Good 80-100 Moderate 60-80 Poor 0-60 y 6 Rationality of diaphragm wall design Reasonable 80-100 General 60-80 Unreasonable 0-60 y 7 Rationality of support design Reasonable 80-100 General 60-80 Unreasonable 0-60 y 8 Reliability of excavation High 80-100 Medium 60-80 Low 0-60 y 9 Reliability of support High 80-100 Medium 60-80 Low 0-60 y 10 Field  14 Steel support axial force <640 kN 640-1000 kN >1000 kN x 15 Overloading (multiples of the allowable value of design) <1.2 1.2-1.5 >1.5 x 16 Construction normalization 80-100 60-80 0-60 x 17 Monitoring (frequency (times/24 h))  Table 4 can be obtained by integrating the five experts' judgements. en, we can obtain the prior probabilities of node X 4 by calculating the principle eigenvector of the fuzzy matrix A � [a ij ] 3×3 (displayed in Table 4), and these probabilities are shown as follows: e prior probabilities of remaining nodes can also be determined and their prior probabilities are listed in Table 5.
After the prior probabilities of each root node are calculated, CPT elicitation for intermediate nodes and leaf node can be carried out using the FAHP and decomposition method mentioned in Section 3.3.3. e calculation of the CPT of node T is taken as an example. e leaf node T has three parent nodes, as illustrated in the FBN model, namely, y 1 , y 2 , and y 3 . First, we should obtain P(T|y 1 ), P(T|y 2 ), and P(T|y 3 ), respectively, according to the method in Section 2.4.3. An example of the calculation of the probabilities of node T when the node y 1 is at state L is presented in Table 6. Similarly, the probability of each state of T when its parent node y 1 is at states M and H can also be obtained (shown in Table 7). In addition, the probabilities of each state of node T can also be obtained under the given states of nodes y 2 and y 3 , and the results are listed in Tables 8 and 9, respectively. en, the conditional probabilities of the leaf node T can be estimated via the decomposition method as follows. For instance, when the states of nodes y 1 , y 2 , and y 3 are all H, we will have the following equations with α � 1/K: where K � P T � L y 1 � H P T � L y 2 � H P T � L y 3 � H + P T � M y 1 � H , T y 1 y 4 y 5 y 6 y 7 y 8 y 10 y 10 y 3 y 2

Root node
Intermediate node Leaf node x 5 x 6 x 7 x 8 x 9 x 10 x 11 x 12 x 13 x 14 x 15 x 16 x 17    In the same way, we can also obtain the conditional probabilities of T on the other states combination of its parent nodes (y 1 , y 2 , and y 3 ). In addition, the CPTs of each intermediate node of BN can be generated similarly. Owing to limited space, only the fuzzy CPT of T is given, as shown in Table 10.

Risk Prediction Based on Forward Reasoning.
After prior probabilities and CPTs are assigned into the established FBN model, we can predict the risk of foundation pit collapse by forward reasoning prior to foundation pit construction. is can help decision makers gain insight into the potential risks associated with the foundation pit, which is very important for decision making in the preconstruction phase. Using (16), the occurrence probabilities of foundation pit collapse with different risk levels under prior probability environments (Scenario 1 displayed in Table 11) are calculated as follows:     Table 9: e conditional probabilities of leaf node T on y 3 's various states.
Obviously, the above computation results mean that the risk state of foundation pit collapse belongs to level L (low) since P(T � L)>P(T � M)>P(T � H). However, the potential risk of foundation pit collapse shows a considerable trend moving towards level M (medium). us, the planners and engineers need to optimize the design and construction plans continuously until the collapse risk probability of level M is relatively small to ensure the safety of the foundation pit construction.
Moreover, the actual status information of root nodes will be revealed gradually as the construction progress evolves. erefore, we can achieve dynamic risk analysis of foundation pit collapse by the new evidence (current actual state of root node) entered into the FBN, which is also the unique merit of the FBN compared to other static analysis methods (FTA, ETA, etc.). When the foundation pit was excavated between YDK27 + 306.221 m and YDK27 + 322.989 m (current evaluation unit shown in Figure 6), there were several irregularities in the current evaluation unit according to the construction plan and on-site inspection. For example, considering the condition of earthwork overbreak, the excavation speed is relatively too fast, the support erection lags, and there are earthmoving vehicles passing along the periphery of the foundation pit. Specifically, the second row of steel support should be set up when the soil is excavated to the level of 11.5 m in light of the design requirements. Actually, the second row of steel support is erected until the soil is excavated to the level of 13 m, as shown in Figure 6 (overbreak part of earthwork). Furthermore, there is overloading around the foundation pit as the earthmoving vehicles pass along the periphery of the foundation pit. erefore, the new status of root nodes X 10 , X 11 , X 12 , and X 15 can be defined according to the current construction states as P(X 10 , X 11 , X 12 , X 15 � M) � 1 (Scenario 2 displayed in Table 11). en, we can calculate the occurrence probabilities of foundation pit collapse with different risk levels under current construction status through forward reasoning. Using (16), these probabilities are calculated as follows:  × P x 1 , x 2 , . . . , x 17 , y 1 , y 2 , . . . , y 10 � 0.14.

(23)
From the above calculation results, it can be seen that the risk level of foundation pit collapse corresponds to M (medium) as P(T � M) > P(T � L)>P(T � H). In addition, it should be noted that the occurrence probability of risk level H increases from 0.08 to 0.14, meaning that the potential risk of serious collapse is increasing. Hence, necessary measures must be adopted to reduce the occurrence probability of foundation pit collapse risk. Specifically, the earth excavation must be stopped immediately; in the meantime, the steel support should be set up in time. Moreover, the supervisors should strictly control the earth excavation speed and pay much more attention to the overloading condition around the foundation pit while limiting it to a reasonable range.

Comparison with Field Actual Status and Numerical
Results. To validate the reliability of the proposed approach, the field monitoring data and numerical simulation results under current construction states (Scenario 2 displayed in Table 11) are analyzed. First, the relationship between collapse risk levels and diaphragm wall deflection is established as illustrated in Table 12 in light of the related design code, the control standards, and the actual construction situation (Nanchang Metro in China). en, the monitoring data of diaphragm wall deflection under current evaluation unit are analyzed to assess the actual risk status of the foundation pit according to the relationship between collapse risk levels and diaphragm wall deflection shown in Table 12. Actually, the maximum diaphragm wall deflection in the displacement monitoring section of YDK27 + 314.605 was 28.7 mm under current evaluation unit, indicating that the risk level of the foundation pit collapse corresponds to level M. Obviously, the result obtained from the proposed approach was consistent with the field actual status of foundation pit in construction. In addition, early warning was issued on-site since the displacement rate exceeded the warning value (>3 mm/d) for consecutive two days under current construction states.
Moreover, to further illustrate the performance of foundation pit and retaining structure under current construction states, a numerical model of typical construction section under current evaluation unit was established using FLAC3D software [53] and is shown in Figure 9. As Mohr-Coulomb model can predict the retaining structure deformation and ground surface settlement of the foundation pit with reasonable accuracy [54,55], the Mohr-Coulomb yield criterion was adopted for the rock and soil material. e diaphragm wall and inner supports belong to structural member and were described by elasticity constitutive model. e horizontal constraint is applied on the left and right boundaries of the model, the front and back boundaries are fixed in the Y direction, and the full boundary constraint is applied on the bottom boundary. In addition, an interface element is arranged between the diaphragm wall and the soil layer to simulate the soil-structure interaction such as sliding and opening. e physical and mechanical parameters of material and structure applied in the numerical model are shown in Tables 13 and 14, respectively. First, the numerical model was calibrated via repeated trial based on the existing monitoring data. On this basis, the horizontal displacements of the diaphragm wall along the depth, considering the condition of earthwork overbreak and overloading around the pit, were calculated adopting the above numerical model. e overbreak depth is 1.5 m, the position of the overloading is at a distance of 5 m from the top of the pit, and the overloading value is 30 kPa. e monitoring and simulating curves of diaphragm wall deflection are shown in Figure 10. It can be seen that the maximum diaphragm wall deflection calculated by FLAC3D is 22.5 mm, and it is also indicated that the risk level of the foundation pit collapse was M according to Table 12. It is apparent that the result obtained from the proposed approach was also in accordance with the numerical result. In addition, the calculation value of maximum diaphragm wall deflection is less than the monitoring value due to only the condition of earthwork overbreak and overloading considered in numerical computation. In other words, the influence of excavation speed, timeliness of support erection, and other uncertainty factors (time-space effect, weather change, etc.) was not considered in numerical calculation. In short, the accuracy and rationality of the proposed approach were verified by comparing with field actual status and numerical results.

Critical Risk Factor Identification and Countermeasures.
In construction practice, decision makers usually pay more attention to which factors play a vital role in foundation pit collapse. However, during the process of actual construction, the critical checkpoints are usually determined based on expert knowledge and practical experience.
rough sensitivity analysis, the BN can automatically identify the critical Mathematical Problems in Engineering risk factors, and this technique is effective and efficient for controlling the risk of targeted events. e SPM values of all the root nodes (X i ; i � 1, 2, . . ., 17) in Scenario 1 are computed by (17), as shown in Figure 11. When T � H, we sort the risk factors according to the SPM value as X 11 > X 12 � X 16 > X 7 > X 6 > X 9 � X 10 � X 14 � X 15 > X 5 > X 1 � X 4 � X 8 � X 17 > X 3 > X 13 > X 2 . It is apparent that X 11 (excavation speed), X 12 (timeliness of support erection), and X 16 (construction normalization) were the top three suspected factors when serious collapse occurred. erefore, it is essential to give priority to the factors X 11 , X 12 , and X 16 during foundation pit construction. In addition, we can facilitate fault diagnosis according to the following route: X 11 ⟶ X 12 ⟶ X 16 ⟶ X 7 ⟶ X 6 ⟶ X 9 ⟶ X 10 ⟶ X 14 ⟶ X 15 . Furthermore, decision makers can adopt corresponding safety control measures in advance on the basis of the above sensitivity analysis results. For instance, importance should be attached to the restriction of excavation speed during excavation, the inner support should be erected in time, and the management of on-site construction should be strengthened in the meantime.

Discussion
Conducting accurate and quantitative risk assessment of foundation pit collapse is difficult or not even possible due to complexities and uncertainties of the foundation pit construction process. erefore, the conventional methods (FTA, ETA, BN, etc.) cannot usually present the reliable and trustworthy results for us in this case. e proposed risk analysis method is superior to the conventional methods, mainly in the following aspects: firstly, the problems of vagueness and subjectivity are well solved by using a FAHP and improved expert elicitation considering a confidence index; in the meantime, a decomposition method is adopted in this study to reduce the elicitation workload and the possibility of inconsistency by simplifying CPT elicitation; secondly, since the occurrence probability of foundation pit collapse risk can be calculated quantitatively using both the prior probability and the actual fault state revealed during construction of the root nodes, the proposed method can achieve real-time deductive reasoning, dynamic risk control, and management during construction; and finally, the critical risk factors can be identified automatically through    sensitivity analysis to help decision makers judge the key checkpoints accurately. e proposed model was applied before and during construction of the foundation pit in Qianhu Avenue Station, and the design and construction schemes were continuously optimized to reduce the collapse risk during foundation pit construction to a maximum extent. e results were in accordance with the field observations and numerical simulation results, indicating that the model proposed in this paper is of certain feasibility and reliability.
e Qianhu Avenue Station construction went smoothly and no collapses occurred throughout foundation pit excavation. In addition, the results of sensitivity analysis indicated that X 11 (excavation speed), X 12 (timeliness of support erection), and X 16 (construction normalization) were the top three risk factors influencing the foundation pit collapse in Qianhu Avenue Station.
Nevertheless, the approach proposed in the present study still has some limitations. For example, this approach relies greatly on domain experts during the establishment of prior probabilities and CPTs of the FBN. Although the subjective bias and uncertainty are reduced to some extent by using the FAHP and improved expert elicitation, some subjective bias still remains. Moreover, the risk factors considered in this study are determined based on the results of statistical analysis of foundation pit collapse cases, expert inquiry, and fault tree analysis. us, some other risk factors may not be taken into consideration due to the shortage of foundation pit cases. Future work will focus on developing a data-driven BN model by collecting a large number of   foundation pit cases along with field monitoring data to further reduce subjectivity and enhance the robustness of the proposed model.

Conclusions
Foundation pit construction is associated with large potential risks owing to various risk events in a complicated uncertain environment. As foundation pit collapse in metro construction can cause heavy casualties and financial losses, it has aroused the attention of all parties. However, uncertain construction environments and data shortages make it difficult to accurately predict the potential collapse risk in foundation pit construction. is study proposes a systematic model for dynamic risk analysis of foundation pit collapse during construction based on a FBN and a FAHP. In terms of priori probabilities and CPTs establishment, a FAHP and decomposition method were adopted to reduce subjectivity and vagueness. In addition, improved expert elicitation considering a confidence index was used to guarantee survey data reliability. A case study concerning an actual deep foundation pit in a metro station was utilized to illustrate a specific application of this model, and the result was in accordance with the field observations and numerical simulation results.
rough the FBN-based risk analysis, the potential risks of foundation pit collapse were predicted quantitatively to help decision makers adopt necessary safety measures in advance; additionally, the critical risk factors were determined automatically to contribute to judging the key checkpoints accurately and effectively. e proposed approach can provide effective decision-making support for planners and engineers and guarantee on-site construction safety of foundation pit. Meanwhile, the results of this study can be used in the field case with similar retaining structure of foundation pit.

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

Conflicts of Interest
e authors wish to confirm that there are no known conflicts of interest associated with this publication.