An Improved Monte Carlo Method Based on Neural Network and Fuzziness Analysis: A Case Study of the Nanpo Dump of the Chengmenshan Copper Mine

+e landslide of dump is a man-made geological disaster which will bring great harm to the surrounding people and environment, and probabilistic reliability analysis is commonly used to analyze the probability of slope landslide or whether protective measures should be taken. Monte Carlo simulation is the most commonly used method, but there are some problems, such as low efficiency, statistical ambiguity of small samples, and the fuzzy transition interval of the stability criterion. +is paper proposes an improved Monte Carlo method that uses an improved bootstrap method to process small samples of geotechnical data, employs ELM (extreme learning machine) based on PSO (particle swarm optimization) to fit the limit equilibrium method function, and constructs the safety factor membership function of the dump site considering the fuzzy transition interval. +is method was applied to an example slope of the dump site in Chengmenshan, Jiangxi. Comparing the analysis result with the result of the traditional MCS (Monte Carlo Search) method, it was found that after adding the safety factor membership function, the result was closer to the actual situation of the dump site, and the probability of failure and reliability index values were closer to those of the dangerous state; after the original function was replaced by the PSO-ELMmodel, the efficiency of theMCSmethod was greatly improved while the results maintained high consistency with the original results; the MCS method combined with the bootstrap method not only simulated the fuzzy uncertainty of the original sample statistics and distribution type but also expressed the reliability index and probability of failure as a two-sided confidence interval with a certain confidence level. +e above conclusion proves the effectiveness and superiority of this method compared with the original MCS method.


Introduction
As a man-made geological disaster with the rise of mining industry, dump landslide not only threatens the life and property safety of mine personnel but also has a huge impact on the surrounding environment. Dump sites in mining areas have been studied by domestic and foreign scholars for stability analysis. At present, the most commonly used methods [1] in engineering are the limit equilibrium method and numerical analysis based on the finite element method, obtaining a certain safety factor value of the analyzed slope and clarifying the sliding surface or part that is most likely to be damaged. However, it often ignores widespread uncertainty of slopes' stability analysis (such as the uncertainty of the calculation parameters of the rock and soil, the geometric model, and the functional relationship between the slopeinfluencing factors and the safety factor), so the effect is not satisfactory. Reliability analysis introduces probability theory into the slope stability evaluation system and obtains the slope failure probability and reliability index based on consideration of the above uncertain factors, so as to achieve a more reasonable and scientific evaluation of the stability of the slope. In recent years, it has been widely valued and studied by scholars. e reliability analysis methods of soil slopes can be roughly divided into the first-order second moment method [2][3][4][5][6][7][8], Monte Carlo method [9][10][11][12][13][14], and response surface method [11,[15][16][17][18][19][20][21][22][23]. Among them, the Monte Carlo method is favored in soil slope reliability analysis due to its accuracy and ease of operation (no matter how complicated the limit state equations and calculations are, according to statistical principles, only enough simulation sampling times and random number sequences are needed to obtain an accurate failure probability value), but some problems also exist.
Firstly, the original limit equilibrium iterative calculation method that links the soil parameters and the slope safety factor requires a certain amount of time, and the Monte Carlo method has to repeat the process hundreds of thousands or even millions of times, so a huge amount of time is needed. Many scholars use models from other fields to simulate and approximate the functional relationship between slope safety factors and soil parameters to solve the problem of low efficiency. He [24] used sample data obtained by deterministic calculations to train a support vector machine (SVM) model to approximate the functional function in the reliability analysis and then combined it with the Monte Carlo method to calculate the probability of slope instability. After the application, it was proved that the accuracy of this method was better than the (first-order reliability method) FORM method, and the efficiency was higher than the ordinary MCS method. Tan et al. [25] proposed a method to replace the original functional function in reliability analysis with an SVM and (radial basis function) RBF model fitted based on two new sample selection methods and combined different examples to illustrate the accuracy and efficiency of the two. Su [26] constructed a Monte Carlo method based on Gaussian regression response surface and verified the effectiveness and efficiency of the method through three slope calculation examples. In addition, artificial bee colony evolution algorithm [27], artificial neural network [28], Kriging [29][30][31], and vector projection [32] have been used to construct the functional relationship between input parameters and output safety factors, alleviating the problem of low efficiency to a certain extent.
Secondly, the traditional MCS method also has the problem of fuzziness. Fuzziness is mainly reflected in two aspects: the traditional reliability analysis using Z < 0 as the criterion to judge slope failure is too arbitrary, ignoring the fuzzy interval of the intermediate transition between stability and instability, and the variability of the sample mean and standard deviation caused by insufficient sample test data will further cause fuzzy uncertainty of the soil parameter distribution types and statistics and ultimately affect the reliability calculation results. In response, some scholars try to solve these problems from the perspective of combining fuzzy mathematics theory and mathematical statistics methods. Habibagahi and Meidani [33], Xu et al. [34], Jia and He [35], Lou [36], Xu [37], Anvar et al. [38], and others successively established the membership function between safety factor and slope stability based on fuzzy mathematics theory. e membership function form has undergone an evolution from linear to nonlinear, triangular to trapezoidal and then to ridge distribution, and it has become closer to the reality of slope stability. However, the ways to determine the undetermined coefficients of the membership function are relatively scarce and subjective. It is necessary to collect as many corresponding slope examples as possible to help judge the rationality of the membership function. Most and Knabe [39], Luo et al. [40], and Tang [41] used the well-known bootstrap theory in statistics to expand the sample in order to reduce the turbulence of the reliability calculation caused by the variability of the parameter sample mean, standard deviation, and distribution, and the reliability index can be more reasonably characterized as a confidence interval with a certain confidence level. is will not only fully simulate the variability of sample statistics, distribution types, and reliability calculation results but also make the slope reliability calculation results more reasonable and true. However, the current combination of bootstrap and Monte Carlo methods is rarely applied, because each bootstrap subsample generated by sampling requires a Monte Carlo calculation, which brings a huge computational burden, so it is necessary to combine some measures to improve the efficiency of the Monte Carlo method.
e Chengmenshan Copper Mine in Jiujiang, Jiangxi Province, due to site selection constraints, had to construct a dump in the lake area near the mining area where the base bearing capacity is weak. In addition, the height of the heap load is relatively high, so it is very easy for progressive landslide instability to be produced over time. Wang [1] conducted a comprehensive reliability analysis on the traditional MCS used in the Chengmenshan dump site (natural unsupported state); the average value of the calculation results under different limit equilibrium functions is P f � 6.65%, β � 1.71, and the results of different methods are highly variable. According to the 1997 US Army Corps of Engineers Index [42], the situation is between unsatisfactory and poor, which is slightly different from the actual situation of the Nanpo dump site (support measures are needed to prevent landslide). e reason might be that the statistical uncertainty caused by insufficient rock and soil samples and the fuzziness of the traditional MCS method were ignored.
According to the question above, an improved Monte Carlo method was proposed that uses improved bootstrap methods to process small samples of geotechnical data, applies particle swarm optimization (PSO)-optimized extreme learning machine (ELM) to fit limit equilibrium method function, and constructs the safety factor membership function of the dump site considering the fuzzy transition interval. It was applied to the reliability analysis of the second-stage dump site of Chengmenshan Copper Mine, and the results were compared with those of the traditional MCS method for the purpose of verifying effectiveness and rationality of the improved Monte Carlo method.

Overview of the Dump.
Chengmenshan Copper Mine is located in Chengmen Town, Jiujiang, Jiangxi. e secondstage dump site is located at the Nanpo ravine and Dachengmen valley on the southeast side of the slope. It is divided into two parts: Nanpo (south slope) and Chengmengou dump sites. e climate in the mining area has four distinct seasons, humid, hot spring/summer, cold, and dry autumn/winter, the annual average temperature is 17°C, and the annual average precipitation is 1420 mm. e basement of the dump site on the south slope is generally a slippery stratum. When the overlying waste material accumulates to a certain extent, under the pressure of its own weight and external factors, it is easy for a slip surface to be produced along the slippery bottom layer to the slope surface area and cause a landslide accident. So far, many bottom heaves of dump slope have been found. e second-stage dump site of the mine has six characteristic sections within the delineated area. Section E-E in the eastern part of the Nanpo dump site was selected for analysis (natural unsupported state). A cross-sectional view is shown in Figure 1, and the statistical characteristics of the 20 experimental datasets and 14 expanded datasets are shown in Table 1.
In the process of parameter sampling, the cohesion and internal friction angle of rock and soil mass meet the linear negative correlation, and the correlation coefficient is taken as −0.5.

Results of Traditional Monte Carlo Slope Reliability
Analysis (Limit Equilibrium) Method. For geotechnical slope engineering, the state function of the slope's stability can be expressed as where the value of the state function Z is the safety factor and X 1 , X 2 , X 3 , . . . , X m are m-many random variables with a certain distribution, which are generally the key factors affecting the stability of a slope such as cohesion and severity of rock and soil. e number of times when the safety factor value Z ≤ 1 is counted as M and the total number of simulations is N, and then the probability of failure can be obtained according to the law of large numbers: e reliability index can be expressed as Using different limit equilibrium methods to calculate the safety factor, it was found that the values obtained by the Janbu simplified method were the lowest, so this method was used for subsequent calculation. Slide is a two-dimensional slope stability analysis software based on limit equilibrium method produced by Canadian Rocscience company and is highly praised in engineering applications due to its advantages of practicality, high efficiency, and accuracy. e reliability analysis results obtained by "slide" software using the MCS method are shown in Figure 2, the failure probability of the Nanpo slope dump of the copper mine is P f � 15.25% and the reliability index is β � 1.026 (also expressed by RI), and the average value of safe factors is 1.045. e results showed that the slope is basically in a stable state at the end of dumping, but there is a lack of sufficient safety factor reserve, so landslide accidents are very easy to occur under static conditions. Compared with the actual situation, the calculation results tend to be conservative.

Problems
Problem 1. It was found that the corresponding probability of failure varies greatly when using different limit equilibrium methods for MCS calculation. By the Janbu simplified method, it was 15.25%, and by the Bishop simplification method and Spencer's method, it was about 1.1%, and the relative error reached 92.4%. e reason should be that the basis for judging the state of the slope based on whether the safety factor is greater than 1 does not conform to the actual situation. e safety factor calculated by the simplified method for the Nanpo slope dump of the Chengmenshan Copper Mine is basically around the critical value 1; although calculated by other methods (such as the Bishop simplified method), it is only 0.05 higher than the Janbu simplified method. e distribution is mostly on the right side of the critical value 1, so there is a problem of large differences in failure probability (as shown in Figure 3). And the accuracy of the calculation results under this criterion will also be affected. So, it is necessary to introduce a stable membership function μ (z) containing the safety factor of the intermediate transition state to improve the judgment method of the slope state in traditional reliability analysis.

Problem 2.
e limit equilibrium program realized in MATLAB not only needs to input the sliding surface coordinates, number of blocks, trial radius, soil parameters, and boundary coordinates and other parameters in advance [43] but also needs to simplify the geometric model boundary to a neatly planned boundary and combine some optimized algorithms (such as genetic algorithm [44]) to help search for minimum safety factors. In this paper, the results of running the MCS method in MATLAB are P f � 15.8% and β � 0.998, and the time consuming is 63005.44 s. e calculation efficiency of the original MCS method has been criticized, and the bootstrap sampling method needs to repeat the whole process of Monte Carlo simulation tens of thousands or even hundreds of thousands of times, which will take a huge amount of time, so other fast and accurate methods are needed to replace the iterative calculation of the limit equilibrium method in MATLAB.
Problem 3. Due to various reasons in actual slope engineering, the experimental data of rock and soil parameters are often limited, so is this dump (there are only 20 sets of data). e statistical uncertainty caused by small sample may lead to deviations in the subsequent Monte Carlo calculation results. e following part of this article will focus on methods and applications to solve the above problems.

Improved Bootstrap Method.
Bootstrapping is an effective method to solve the statistical uncertainty of small sample data. e idea of the bootstrap method is to randomly sample the initial samples with replacement to obtain a large number of bootstrap subsamples (the sample size is the same as the original sample) that contain the original sample information, then calculate the estimated value of statistics and Akaike information criterion (AIC) value for each subsample, get their optimal probability density distribution according to the AIC [45], and finally perform Monte Carlo calculation to  obtain their reliability index, probability of failure, and corresponding confidence interval. Assume the sample of soil parameters X � (c i , φ i , c i ), i � 1, 2, . . . , N , then randomly sample with replacement N times as shown in Figure 4 to obtain a bootstrap subsample B j � B 1,j , B 2,j , . . . , B N,j , and repeat this step M times to obtain M subsamples. e theoretical basis and good convergence of the bootstrap method have long been proved by scholars, but it also has the problem that the sampling range is small and the probability distribution is concentrated on a small number of points for small samples (sample size of 10-30), which causes the calculation results to deviate from the true distribution. Some scholars [46][47][48] proposed an improved    Mathematical Problems in Engineering 5 bootstrap method for this problem: first arrange the original sample in ascending order to get the order statistic . , x n ), then, for each observation x i , a larger and a smaller statistic are randomly selected in its neighborhood U i according to the rule of uniform distribution to expand the sampling range and obtain information outside the observation point: e bootstrap method can continue to be used for parameter estimation based on the improved samples X′ � (x 1 ′ , x 2 ′ , . . . , x n ′ ) obtained by formula (4). Since the original sample size of 20 is less than 30, thus is a small sample, the size can be expanded by formula (4) before bootstrap sampling, but not every observation value needs to be expanded in this way, especially when it is very close to (or even the same as) an adjacent value. Use MATLAB for random permutation and combination after the six groups of data were, respectively, expanded, and the expanded sample statistics are shown in Table 1.
One of four types of functions commonly used in the distribution of rock and soil parameters, normal, lognormal, extreme value I-type, and Wilber distribution, is selected by calculate the AIC value of the original sample (the probability distribution function with the smallest AIC value can be considered as the one that best fits the probability distribution of the test data). It is found that the AIC value corresponding to the normal distribution is the smallest and the optimal probability distributions of the six types of data of the initial sample are all normal distributions, which conforms to the description in the literature [1]. e subsample mean, standard deviation, and optimal probability density function distribution type are recorded after each bootstrap sampling and used as the statistics and distribution type of random variables in the subsequent Monte Carlo calculation.

Original sample
Bootstrap sub-sample

PSO-ELM Model. Extreme learning machine (ELM) is a new type of feedforward single hidden layer neural network.
In addition to inheriting the good self-organization and selfadapting ability of general neural network algorithms, it also has fewer adjustable parameters, faster speed, better prediction accuracy, strong versatility, and other advantages. e input weight ω l×n and hidden layer threshold b l×n are randomly selected, and then the input quantity X n×Q , output weight β l×m , and excitation function g(x) are combined to obtain the expression of the output t j : where β � H + T.
In ELM, the input weight hidden layer threshold is randomly selected, and the excitation function and input are known.
erefore, the training of the extreme learning machine is to calculate the output weight through equation (7), and then use it to perform predictions.
PSO is an intelligent optimization algorithm that simulates the foraging behavior of a bird colony to achieve swarm optimization. is method treats the solution of each problem as a particle. e distance between the particle's spatial position and the target is used as the fitness value, and each particle has a specific moving direction and speed. e PSO algorithm updates the speed and fitness in every search by tracking the individual's extreme values (Pbest, the optimal position of fitness values calculated from the positions of individuals) and the group's extreme values (Gbest, the optimal position of fitness of all particles in the population) until the termination requirements are met. e termination condition can be the number of iterations, the extreme value error of two consecutive searches below a specific tolerance value, or a mixture of both. e extreme learning machine algorithm optimized by particle swarm optimization is shown in Table 2.

Membership Function of Safety Factor in Dump Sites.
Monte Carlo reliability analysis considers that there are only three states in a slope: limit equilibrium state (Z � 1), unstable state (Z < 1), and stable state (Z > 1). However, the corresponding relationship between the safety factor and the stable situation in the actual soil slope is fuzzy and uncertain, which means the slope has a fuzzy interval of intermediate transition between completely unstable and completely stable. Even if the safety factor is greater than 1 or even 1.2, it may become unstable; less than 1 also has the possibility of stability. erefore, the way the Monte Carlo method estimates slope failure probability based on the proportion of statistical failure times to total simulation times may be biased, so it is necessary to introduce the stability membership function μ(Z) of the safety factor to judge the stability of soil slope equation (2) can be improved as As the name implies, the membership function is the degree to which the safety factor belongs to the stable state of the slope, μ(Z) ⟶ 0. Slope becomes more unstable, μ(Z) � 0.5. e probability of slope instability and stability is 0.5. At this time, the fuzziness is the strongest and the stability state is the most difficult to judge. At μ(Z) ⟶ 1, the slope becomes more stable. e membership functions commonly used to characterize the stability of geotechnical engineering structures include ridge distribution, quadratic parabolic distribution, Mathematical Problems in Engineering and reduced semitrapezoidal distribution. Each of these membership functions can also be divided into three types: descending, rising, and intermediate. e membership function of the safety factor for stability should express the degree to which the safety factor belongs to the concept of stability, and additionally, the slope of the graph curve should reflect the strength of the fuzziness that affects the judgment of the steady state. e smaller the safety factor, the closer μ (Z) is to 0, the easier it is to judge that the slope is in an unstable state. Similarly, the larger the safety factor, the closer μ (Z) is to 1, and the slope is likely to be in a stable state. e fuzziness is the strongest when the safety factor and μ (Z) are in the middle, so the curve slope of the membership function first increases and then decreases, with the characteristics of slow at both ends and steep in the middle. Based on the above characteristics and engineering application principles, the raised-ridge-shaped distribution membership function was selected as the stable membership function of the slope safety factor. e function expression is as follows: e steps to determine the specific form in the membership function are as follows: Step 1: collect the corresponding relationship between the actual stability of the dump site and the value of the safety factor. We collected 48 sets of statistical data on the stability factors of different dumps in China, as shown in Table 3. In the data, the highest value of safety factor of unstable slope is 1.21, and the lowest value is 0.92; therefore, it can be assumed that a dump slope with safety factor greater than 1.21 must be stable, and a dump slope with safety factor less than 0.92 can no longer remain stable. e median value of 1.065 is used to divide the safety factor into two intervals. e stable state of dump sites changes most obviously near 1.065, so the stability state of a dump slope with a safety factor of 1.065 is the fuzziest, and the corresponding membership degree should be 0.5.
Step 2: preliminarily determine the undetermined coefficient of the membership function. According to the analysis given in Step 1, a � 0.92, c � 1.21, b � (a + c)/2 � 1.065, and substituting into equation (9), the membership function of the safety factor of raised-ridge-shaped distribution can be expressed as Figure 5 shows a ridge-shaped distribution function diagram of the relationship between the safety factor and the steady state membership based on the collected data. e point with a membership value of 1 corresponds to the stable state of the dump, and the point with a membership value of 0 indicates an unstable state. e function diagram can fit the data points of the dump sites well and follow the slope membership function law mentioned above.
Step 3: use the degree of conformity between the statistical law of the dump sites and the safety factor interval, with the membership interval corresponding to the fuzzy language value, to judge whether the membership function determined in Step 1 and 2 is practical.
ere are nine commonly used fuzzy linguistic values to describe a certain state; however, due to the small amount of collected samples, the number of samples in each fuzzy subset would be too small, which could cause fluctuations in statistical results, so it is necessary to reduce the fuzzy language subset describing slope stability. In this paper, five fuzzy linguistic values of stable π 1 , basically stable π 2 , critical π 3 , basically unstable π 4 , and unstable π 5 were selected to describe the stability state of the dump slope, and the membership interval calculation formula of fuzzy linguistic value in [37] was reduced to 5, as shown in equation (11), where a is a constant used to separate the interval, a ∈ [0.5, 1], and here we take a � 0.51: Table 2: Particle swarm optimization-extreme learning machine (PSO-ELM) algorithm program.
Algorithm: algorithmic flow of PSO-ELM (1) Obtain the training and testing datasets (2) Begin ELM training (3) Set ELM parameters randomly (4) Use the mean square error (MSE) as the fitness function (5) Initialize PSO population (Inipop) (6) Calculate the fitness value of each candidate solution (7) S � global best solution (8) For i � 1 to maximum iteration number do (9) For i � 1 to P do (10) Update the velocity and position of the ith particle (11) Evaluate the fitness of the ith particle 1 (2) Update the personal best solution of the ith particle (13) S � current global best solution (14) End for (15) End for (16) End (17) Obtain the optimal input weights and hidden biases of hidden layer neurons using S (18) Use the optimal input weights and hidden layer neurons for ELM test 8 Mathematical Problems in Engineering e safety factor interval corresponding to the fuzzy language value of the dump slope stability can be obtained by substituting each end value of the membership interval into equation (12). Table 4 lists the statistical results of the dump sites' samples, the membership degree interval, and the safety factor interval, where n i represents the number of slopes in the dump in a stable state in the π i safety factor interval and N i is the number of slope instances in the corresponding safety factor interval. If the membership function is selected reasonably, the value of n i /N i increases to 1 as the membership interval approaches 1: It can be seen from Table 4 that the n i /N i value is in good agreement with the membership interval corresponding to the ith fuzzy language value. e closer it is to the steady state, the closer its value is to 1. erefore, it can be considered that the slope stability membership function proposed in this paper according to formula (10) can reflect the overall law of the slope stability of the dump site with the change of the safety factor, which can be used to improve the criterion of the subsequent Monte Carlo calculation state function. e flowchart of specific application steps of methods is shown in Figure 6.

Bootstrap Statistical Uncertainty Simulation.
is paper used MATLAB to perform bootstrap sampling on the original sample and selected N B � 10 4 as the number of subsamples (also the number of samples) based on ensuring good convergence and computational efficiency. Box-andwhisker plots of 20 randomly selected bootstrap subsamples and the original sample are shown in Figure 7. It can be seen that the generated subsamples perfectly retain the data information of the original sample.
For this procedure, 10,000 simulations generated 10,000 subsamples, corresponding to the mean, standard deviation, and AIC value of 10,000 simulated samples. Figures 8 and 9 and Table 5 show the distribution of the subsamples' mean and standard deviation, and the number of times that the four common probability density functions are identified as the optimal distribution. It is not difficult to see that except for bedrock cohesion, for the other indicators, the distribution of the mean and standard deviation of the bootstrap subsamples is close to that of the original sample, and the optimal probability distribution function of the subsamples obtained by the bootstrap method is mostly the same as the original sample, which is normal distribution, only a few of the optimal probability distributions of the subsample are the other indicators.
is shows that the bootstrap subsamples can effectively reflect the basic characteristics of the original sample data, which means that the original limited data sample has statistical uncertainty and AIC value variability. is problem may affect the selection of subsequent  Mathematical Problems in Engineering distribution parameters and probability distribution function types, and the bootstrap method can simulate such variability.

Realizing the PSO-ELM-MCS Model in MATLAB.
is paper uses the Janbu simplified method to calculate the safety factor of the sample to form the training and test sets and train the PSO-ELM model to achieve the purpose of replacing the original limit equilibrium function. e specific steps are as follows: Step 1: determine training and test sets. In addition to the 34 sets of original samples, 44 sets of samples randomly selected from the bootstrap subsamples were substituted in the Slide software to calculate the safety factor to form 70 sets of training samples and 8 sets of test samples.
Step 2: set initial parameters. Six rock and soil parameters corresponded to one safety factor value in each group of samples, so there were six input layer units of the PSO-ELM model and one output layer unit. If the number of hidden layer units is too high, it will lead to excess performance, and if the number is too low, it will affect the prediction accuracy. After many calculations, we decided to set the number to 25. e learning factor C 1 � C 2 � 1.49, the sigma function was selected as the excitation function, and the linearly decreasing weight that can account for both local and global search was selected as the inertial weight.
Step 3: train and test models. e 70 sets of training samples obtained in Step 1 were combined with the PSO model to start training, and the weights and thresholds obtained after training were assigned to the ELM model, then the test samples were substituted into    Figure 10 shows a comparison between the predicted and real values of the PSO-ELM model, as well as the predicted value of ordinary ELM. e predicted values of the safety factor obtained by the two models were not much different from the true value, but the PSO-ELM model was better than the ELM model in terms of the root mean square error (RMSE) of prediction accuracy or the coefficient of determination (R 2 ), which means particle swarm algorithm optimization of extreme learning machine does improve accuracy. e reliability analysis result obtained by running the PSO-ELM-MCS method in MATLAB was as follows: probability of failure P f � 15.02%, reliability index RI � 1.0356, and average safety factor 1.0454, which are highly consistent with the results obtained by the Slide software in the previous paper, indicating that the PSO-ELM model can perfectly fit the original Monte Carlo method to express the functional of the relationship between soil parameters and safety factors and reliability indicators. Meanwhile, the time to obtain the safety factor by the limit equilibrium method and the PSO-ELM model in MATLAB was statistically compared, and it was found that the efficiency of the latter was significantly better. e PSO-ELM model only required 0.0741 s, while the limit equilibrium method took 12.03 s. Figure 11 shows the calculation results of the improved PSO-ELM-MCS model after adding the slope stability membership function as the Monte Carlo calculation stable state criterion and expresses the relationship between the safety factor value and the membership degree. Comparing the probability of failure P f � 36.85% and reliability index β � 0.34 with the calculation results of the PSO-ELM-MCS model (P f � 15.02%, β � 1.0356), it can be found that the failure probability of the dump site is large improved, and reliability is significantly decreased. Comparing the reliability calculation results obtained by different limit equilibrium methods (shown in Table 6), it is found that probability of failure and reliability indicators maintain good consistency under different limit equilibrium methods, which is obviously different from the previous results obtained by using ordinary Monte Carlo reliability analysis (the relative error dropped from more than 90% to about 3%). at is to say, the results of Monte Carlo calculation with the stability membership function of the safety factor added are closer to reality than the traditional Monte Carlo reliability calculation, and it will not have the illusion of stability for simulators and engineers.

Monte Carlo Calculation Combined with PSO-ELM, Bootstrap Method, and Fuzzy Membership Function.
is paper calculated the probability of failure P f,i and reliability index β i (i � 1, 2, 3, . . . , N B ) of each bootstrap subsample and obtained the distribution probability density function and confidence interval (using the reliability index and the 5% and 95% quantile values of failure probability as the upper and lower limits of the 90% two-sided confidence interval), and the results as shown in Figure 12 and Table 7.
It can be seen that the reliability index of the dump site basically varies between −0.5 and 1, and the probability is the highest near 0.25, with an average value of 0.243, lower than the average value of 0.34 without bootstrap sampling. e probability of failure varies between 0.1 and 0.8, with the highest probability near 0.4, and the average value is 0.4075, slightly higher than the 0.3664 without bootstrap sampling.
is shows that the reliability index and failure probability obtained by the ordinary MCS method may have variability, and the bootstrap method can simulate this variability. In addition, compared to the MCS method, which generally only obtains a single reliability index and failure probability value, the MCS method combined with bootstrap can express the reliability index and probability of failure as a twosided confidence interval with a certain confidence level, which can reflect the actual reliability level of the dump slope.

Discussion and Conclusion
is paper proposed an improved Monte Carlo method that uses improved bootstrap methods to process small samples of geotechnical data, applies particle swarm optimization (PSO)-optimized extreme learning machine (ELM) to fit limit equilibrium method function, and constructs the safety factor membership function of the dump site considering the fuzzy transition interval, the results of each method combination are compared and analyzed, and the conclusions are as follows: e MCS method using the PSO-ELM model to replace the original functional function can greatly improve efficiency on the basis of ensuring accuracy. e time consumed changed from the original 63005.44 s to 45.312 s, and the probability of failure (15.25% and 15.02%) and the reliability index (1.027 and 1.036) are almost the same, which proves that the artificial intelligence and neural network algorithms have good self-learning abilities and are effective at fitting implicit functional functions. After the safety factor membership function of the dump site was added to the Monte Carlo method, the fluctuation of reliability index and probability of failure of the Nanpo dump site under different limit equilibrium methods was effectively improved, and the results showed an inclination toward slope instability (the failure probability of 36.64% and reliability index of 0.34 are significantly lower than the 15.02% and 1.036 without the membership function), which is closer to the situation at the actual dump site. is shows that the addition of a fuzzy transition zone of safety factor can take into account the uncertainty and gradualness of the slope in the process of failure, reduce the fluctuation of calculation results to a certain extent, and make it more real and effective. Although the MCS method combined with the bootstrap method has a lot of repeated calculation work, the time changed from 44.519 s to 20275.23 s, but time was exchanged for accuracy, and the bootstrap subsample was generated to simulate the uncertainty brought by the small sample size of the experimental data (uncertainty of sample mean, standard deviation, optimal probability density distribution, reliability index, and failure probability), making the reliability calculation result change from a single value to a confidence interval with fuzzy characteristics, which is closer to reality.
e average values of reliability analysis obtained by the improved Monte Carlo method based on bootstrap in the second-phase dump site of the mine are P f � 40.75% and β � 0.243, which indicates a hazardous situation that requires immediate protective measures. e dump site immediately implemented preventive measures of antislip steel rail piles plus spoil grading and classified discharge to avoid the aggravation of displacement and the occurrence of landslides.
On summary, the method proposed in this paper can be used in the reliability analysis of the dump, and the effect is better than the traditional MCS method.

Limitation
However, the methods and research in this paper still have some shortcomings that could be improved. Such as the calculation process of MCS is not simplified and still consumes a lot of time, and the fuzzy membership function of the dump safety coefficient determined in this paper may still be subjective because of insufficient data collection and requires a professional authority to collect thousands of waste dump examples to reduce uncertainty. Follow-up studies will be carried out to improve the above deficiencies.

Conflicts of Interest
e authors declare that they have no conflicts of interest.