A Comprehensive Evaluation of Resilience in Abandoned Open-Pit Mine Slopes Based on a Two-Dimensional Cloud Model with Combination Weighting

: The stability of abandoned open-pit mine slopes and their ecological environment are threatened owing to their fragile, complicated, and uncertain characteristics. This study establishes a novel evaluation indicator system for enhancing mine design and environmental protection insight. The weights in the system are assigned using a combined method, which consists of the game theory, the interval analytic hierarchy process (IAHP), and the entropy weight method (EWM). The IAHP is optimized by the improved radial movement optimal (IRMO) algorithm and the simulated annealing (SA) algorithm to ensure calculation stability and efficiency. Meanwhile, a two-dimensional cloud model (TDCM) is developed to obtain the slope resilience level and visualize the result. This comprehensive evaluation method is applied to three abandoned mine slopes in the Yellow River Basin, and the results demonstrate that the method can provide crucial insights for rational mine slope stabilization and ecological restoration.


Introduction
Mining is one of the foundational industries for a country's economic and social development, with open-pit mining being one of its primary methods.Improper mining practices have damaged mines and their ecological surroundings, forming high and steep slopes.This has increased the risk of landslides and incidents, posing significant hazards.The failure of mine slopes can cause substantial economic losses and even endanger the safety of surrounding residents [1][2][3][4].Meanwhile, as people's awareness of environmental protection gradually strengthens, mine ecological environment management has become essential to environmental protection construction.The ecological restoration of a mining area is an integral part of the ecological civilization construction process.Therefore, conducting a comprehensive evaluation of mine slopes is a complex interdisciplinary endeavor [5] requiring assessment from various perspectives, such as safety and ecological restoration suitability [6].This is of paramount significance in mine design and production processes.
In the mining industry, which focuses on safety production, reliable safety evaluation is a prerequisite to ensure the orderly and safe production of mines [7].The safety evaluation is usually viewed as a multi-criteria decision-making problem containing both qualitative and quantitative characteristics [8].A single risk evaluation method cannot satisfy the evaluation requirements of complex project targets, so a composite evaluation method combining multiple methods is often used.Researchers usually combine the AHP with other methods for safety evaluation.The authors of [8] utilized extended interval trapezoidal fuzzy soft set and the fuzzy hierarchy analysis method (FAHP) to calculate the risk factor weights of high and steep slopes in open-pit mines.The authors of [9] used the AHP, the topology object element method, and EWM to comprehensively evaluate the safety of high, steep rocky slopes, forming a multi-indicator dynamic safety evaluation method.The authors of [10] assigned causative factors to sub-factors and weightage according to the AHP.A landslide hazard map along the national highway was prepared to evaluate slope safety using the AHP model.However, the AHP involves subjective uncertainty to some extent.To mitigate the impact of subjective judgment and uncertainty, the IAHP constructs interval number judgment matrices [11], which better reflect the subjective preferences of experts and the uncertainty of the evaluation criteria.Meanwhile, the EWM [12] is an objective weighting method that determines the weights of each criterion based on the information content of the evaluation indicators, possessing advantages such as objectivity and adaptability.This paper will employ a combination of the IAHP and EWM for weighting.However, the IAHP's interval number judgment matrices may lack precision or be computationally complex [13].Therefore, to enhance the accuracy of subjective weighting and simplify computation, the IRMO algorithm will be used to establish an optimization model for objective search, obtaining optimal subjective weights for each evaluation indicator, thereby complementing both subjective and objective biases.
The current study uses a retrieval strategy to perform a keyword co-occurrence analysis through VOSviewer.The search query "mining/mine slope evaluation/assessment" is employed to query the Web of Science database.As illustrated in Figure 1, the current evaluations of slopes predominantly focus on safety or stability, with less emphasis on the ecological restoration suitability of mining slopes.However, following a series of mining activities, the surrounding ecological environment becomes fragile, and rock structures are fractured.This mining process parallels the seismic hazards experienced by urban buildings, resulting in devastating impacts.An evaluation of the suitability of ecological restoration will help to select appropriate restoration measures to ensure that the surrounding ecological environment can be effectively restored and protected [14].Therefore, the ecological restoration suitability of slopes should also be included in the evaluation object, as well as safety, to ensure the stability and health of their surrounding ecological environment.This paper combines safety and ecological evaluations of abandoned open-pit mine slopes to comprehensively evaluate slope resilience from two dimensions based on the resilience theory.Resilience was initially used to describe the mechanical properties of metals.It was later applied to systems ecology [15][16][17] to describe the ability of ecosystems to maintain normal functioning or restore balance after being damaged.The resilience evaluation encompasses both risk and recovery assessments.In 2020, China adopted the construction of "resilient cities" as a national sustainable development strategy.The continuous development of the resilience theory has made exploring cities' or systems' ability to cope with destruction a hot topic in disaster prevention and mitigation [18], spanning various applications such as deep foundation pit construction [19], buildings [20], critical infrastructure systems [21], urban systems [22,23], bridges [24], etc.Based on this, an evaluation of the slope resilience of abandoned open-pit mines can be based on two dimensions: slope safety and slope ecological restoration suitability.This has significant theoretical and practical significance for safety management and ecological environment restoration.In addition, due to the evaluation process not fully considering the impacts of evaluative fuzziness and uncertainty, the assessment results are limited.As a cognitive model based on the normal distribution function and membership function, a cloud model [25] can describe the uncertainty of quantitative information and can be used to describe qualitative concepts, achieving bidirectional transformation between the two [26][27][28].In recent years, a cloud model has been gradually introduced into the evaluation of engineering projects [29], such as subways [30][31][32], slopes [11,33], railways [34], and urban engineering construction projects [35].A TDCM is synthesized from two one-dimensional clouds, comprehensively describing the randomness and fuzziness issues under the combined influence of two factors.This paper employs the combination weighting method to derive the weights.It constructs a TDCM based on the cloud theory from both safety and ecological perspectives, comprehensively evaluating the resilience of abandoned open-pit mine slopes.
This paper addresses the inherent uncertainties in current slope evaluations by constructing a comprehensive resilience evaluation model for abandoned open-pit mine slopes.This study leverages the combination weighting method with a TDCM to enhance the accuracy of assessments, mitigating issues such as subjectivity, uncertainties, fuzziness, and randomness that may arise during the evaluation process.In terms of evaluation methods, this study proposes the subjective weight optimization method based on the IAHP-IRMO-SA.This hybrid algorithm is applied to calculate the subjective weights in the interval judgment matrix.The corresponding algorithmic calculation program is developed, leveraging the self-feedback ability of particles in the IRMO algorithm and the SA algorithm to escape from local optimal solutions based on the Metropolis criterion to accurately determine the optimal solution for subjective weights.
This paper is organized as follows: Section 2 introduces the data and methods, including the process of constructing the toughness evaluation indicator system and selecting the grading criteria.The Section 3 shows the establishment of the toughness assessment model for abandoned open-pit mine slopes based on the IAHP-IRMO-SA.The Section 4 presents a case study to verify the validity and practicality of the provided safety assessment model.The results are discussed with the main conclusions and insights in the Sections 5 and 6.In addition, due to the evaluation process not fully considering the impacts of evaluative fuzziness and uncertainty, the assessment results are limited.As a cognitive model based on the normal distribution function and membership function, a cloud model [25] can describe the uncertainty of quantitative information and can be used to describe qualitative concepts, achieving bidirectional transformation between the two [26][27][28].In recent years, a cloud model has been gradually introduced into the evaluation of engineering projects [29], such as subways [30][31][32], slopes [11,33], railways [34], and urban engineering construction projects [35].A TDCM is synthesized from two one-dimensional clouds, comprehensively describing the randomness and fuzziness issues under the combined influence of two factors.This paper employs the combination weighting method to derive the weights.It constructs a TDCM based on the cloud theory from both safety and ecological perspectives, comprehensively evaluating the resilience of abandoned open-pit mine slopes.
This paper addresses the inherent uncertainties in current slope evaluations by constructing a comprehensive resilience evaluation model for abandoned open-pit mine slopes.This study leverages the combination weighting method with a TDCM to enhance the accuracy of assessments, mitigating issues such as subjectivity, uncertainties, fuzziness, and randomness that may arise during the evaluation process.In terms of evaluation methods, this study proposes the subjective weight optimization method based on the IAHP-IRMO-SA.This hybrid algorithm is applied to calculate the subjective weights in the interval judgment matrix.The corresponding algorithmic calculation program is developed, leveraging the self-feedback ability of particles in the IRMO algorithm and the SA algorithm to escape from local optimal solutions based on the Metropolis criterion to accurately determine the optimal solution for subjective weights.
This paper is organized as follows: Section 2 introduces the data and methods, including the process of constructing the toughness evaluation indicator system and selecting the grading criteria.The Section 3 shows the establishment of the toughness assessment model for abandoned open-pit mine slopes based on the IAHP-IRMO-SA.The Section 4 presents a case study to verify the validity and practicality of the provided safety assessment model.The results are discussed with the main conclusions and insights in the Sections 5 and 6.

Basis for Indicator Selection
This paper constructs a comprehensive evaluation model from two dimensions of slope safety and ecological restoration suitability.To share a standard system of slope safety and ecological evaluation indicators to ensure objectivity, this paper evaluates the slope ecological restoration suitability side by side through the existing slope safety evaluation method.Hence, the slope safety and ecological evaluation indicators share a standard system.This study uses existing research [36][37][38][39][40][41][42][43][44][45][46][47][48][49][50][51][52] as a reference to select a comprehensive evaluation indicator for abandoned open-pit mine slopes, based on which the assessment indicator system in Table 1 is established.

D4
Degree of ecological environment vulnerability [49,52] The resilience of abandoned open-pit mine slopes is described using four aspects: geological features, slope geometric features, slope rock features, and external factors.In the proposed evaluation indicator system, geological features (primary indicator A) include the degree of weathering of the slope rock body, the damaged area of the abandoned mine, and the degree of development of the slope joints and cracks, which can reflect the influence of the geological features of the abandoned open-pit mine slopes on the safety and ecological restoration of the abandoned open-pit mine slopes.The geometric characteristics of the slope (primary indicator B) are divided into slope length and slope height, which characterize the geometric characteristics of the abandoned mine slope.Slope rock mass features (primary indicator C) includes the hardness of the slope rock, the basic quality grade of the rock, the self-stabilizing ability of the vertical slope, the classification of the integrity of the slope rock, and the cohesive force of the rock and soil body, which can reflect the rock stability of the abandoned open-pit mine slopes.The external factors (primary indicator D) include the average annual precipitation, the seismic intensity, the risk of geologic hazards, and the degree of vulnerability of ecological environments, which can reflect the impacts of the external factors of abandoned open-pit mine slopes on the safety of mine slopes and the suitability of ecological restoration.The influence of external factors on the safety and ecological restoration suitability of abandoned open-pit mine slopes can be more comprehensively reflected.

Indicator Grading Criteria
The system consists of 15 evaluation indicators.Due to the different classification standards of each indicator, it is necessary to hierarchically quantify evaluation indicators and determine a unified, comprehensive evaluation criterion for the resilience of open-pit mine slopes.Using the literature [53,54], the established evaluation indicator system was quantitatively analyzed to obtain the range of values of the secondary indicators in the index layer, as shown in Table 2.There are two types of indicators: "qualitative (QL)" and "quantitative (QN)".In this paper, the grading criteria for the evaluation indicator system are unified concerning the principles of risk grading and control [41,55].The comprehensive evaluation grade for the resilience of abandoned open-pit mine slopes is divided into four levels, Level I, Level II, Level III, and Level IV, as shown in Table 3.

Comprehensive Evaluation Model of Resilience Based on Combination Weighting and TDCM
The established model solved the indicators' optimal subjective weights using the IAHP-IRMO-SA method.The EWM was used to calculate the objective weights; the comprehensive weights were obtained based on the game theory.The TDCM was used to generate the cloud droplets and make the two-dimensional standard cloud diagram.After calculating the cloud eigenvalues of each evaluation indicator, the combined weight of these eigenvalues was calculated.Subsequently, the level of safety and ecological restoration suitability of the slopes by the degree of affiliation and proximity and the comprehensive grade were determined.The specific process can be seen in Figure 2.

Comprehensive Evaluation Model of Resilience Based on Combination Weighting and TDCM
The established model solved the indicators' optimal subjective weights using the IAHP-IRMO-SA method.The EWM was used to calculate the objective weights; the comprehensive weights were obtained based on the game theory.The TDCM was used to generate the cloud droplets and make the two-dimensional standard cloud diagram.After calculating the cloud eigenvalues of each evaluation indicator, the combined weight of these eigenvalues was calculated.Subsequently, the level of safety and ecological restoration suitability of the slopes by the degree of affiliation and proximity and the comprehensive grade were determined.The specific process can be seen in Figure 2.  The IAHP constructs the interval number judgment matrix to reduce the impacts of subjective judgment and uncertainty.This paper establishes a hierarchical model based on the IAHP, constructs the interval number judgment matrix, and introduces the concept of the degree of divergence to establish the objective optimization model.The model's objective function is the IRMO algorithm's fitness function.The IRMO algorithm is used to optimize the search and obtain the optimal subjective weights of the evaluation indicators, and the specific steps are described below.
(1) Construct a judgment matrix of intervals.The 1-9 scale method compares the importance degree of each evaluation indicator at the same level and assigns values.The interval number judgment matrix A is the relative importance of the two comparisons of indicator i and indicator j.

IAHP method
The IAHP constructs the interval number judgment matrix to reduce the impacts of subjective judgment and uncertainty.This paper establishes a hierarchical model based on the IAHP, constructs the interval number judgment matrix, and introduces the concept of the degree of divergence to establish the objective optimization model.The model's objective function is the IRMO algorithm's fitness function.The IRMO algorithm is used to optimize the search and obtain the optimal subjective weights of the evaluation indicators, and the specific steps are described below.
(1) Construct a judgment matrix of intervals.
The 1-9 scale method compares the importance degree of each evaluation indicator at the same level and assigns values.The interval number judgment matrix A = (A ij ) n×n , A ij is the relative importance of the two comparisons of indicator i and indicator j.
Set the number of intervals as a = [a l , a r ] and b = [b l , b r ]: where Equation ( 3

b. IRMO algorithm
The IRMO is a global optimization algorithm used to quickly solve the optimal value of multi-dimensional objective functions and aims at the problems of unstable search results and the insufficient search accuracy of the Radial Movement Algorithm (RMO).The data structure of the IRMO is further optimized and adjusted to enhance the self-feedback ability between particles.The search stability and accuracy are improved.Figure 3 shows the basic principle diagram of the IRMO [56].
( , ) D a b -The degree of divergence between the sum of interval numbers.When D a b is the degree of divergence from point b to interval a.The objective optimization model is established by combining the concept of the degree of divergence between interval numbers: where Equation ( 3) is the objective function in the objective optimization model; U w is the optimal subjective weight of each assessment indicator; ij W denotes the judgment range when indicator i and indicator j are compared; and w and Uj w are the optimal subjective weights of indicator i and indicator j, respectively, and is the degree of divergence of intervals ij A and ij W .When the optimal subjective weight is taken for both indicator i and indicator j, ( ,

b. IRMO algorithm
The IRMO is a global optimization algorithm used to quickly solve the optimal value of multi-dimensional objective functions and aims at the problems of unstable search results and the insufficient search accuracy of the Radial Movement Algorithm (RMO).The data structure of the IRMO is further optimized and adjusted to enhance the self-feedback ability between particles.The search stability and accuracy are improved.Figure 3 shows the basic principle diagram of the IRMO [56].In this paper, when applying the IRMO algorithm to optimize the search, the objective function denoted as Equation ( 3) in the optimization model is used as the fitness func- In this paper, when applying the IRMO algorithm to optimize the search, the objective function denoted as Equation ( 3) in the optimization model is used as the fitness function of the IRMO algorithm.The fitness function computes the fitness function value corresponding to each generation of particles, and the calculated fitness function values are compared to select the contemporary optimal solution.The solution space is reduced to one point when the algorithm calculates the last generation.The fitness function value corresponding to this point is the global optimal solution, and the parameter corresponding to it is the optimal subjective weight of the evaluation indicator w U .To accurately search and calculate the optimal subjective weight of the evaluation indicator w U , it is necessary to constrain the range of values of each variable.According to the constructed objective op-timization model, each variable must satisfy the constraints of Equation ( 4) in the objective optimization model.

c. IRMO-SA
The IRMO is an excellent heuristic algorithm with good global search ability, more stable search results, and rapid convergence, and it is not easily premature.It can quickly search out the global optimal solution, but it may fall into the local optimal situation during the search process [56].SA is based on the Metropolis criterion and can jump out of the local optimal solution while searching for a global optimal solution of the objective function [57,58].In this paper, by combining the two algorithms to improve the overall performance, a hybrid IRMO-SA algorithm is designed to solve the optimal weight solution based on the IAHP.The specific steps are described below.
(1) Construct the objective function of the IRMO-SA algorithm.
In this paper, Equation ( 3) is used as the objective function of the IRMO-SA algorithm as follows: (2) Obtain the optimal weight solution of the interval number judgment matrix based on IAHP-IRMO-SA.
In the IRMO-SA algorithm, N particles are generated in each generation, and the position of each particle in the space represents a solution vector X i,j .These solution vectors are composed of the variables w Ui (i = 1, 2, . . ., n) with the dimension of the interval number judgment matrix constructed using the IAHP method with the dimension of M. The information of all of the particles' positions is shown in Equation (7).When the IRMO-SA algorithm is applied to optimize the search, the objective function Equation ( 6) is used as the fitness function of the IRMO-SA algorithm.The fitness value of the objective function corresponding to the k-th generation of particles is calculated using the fitness function.The computed fitness values of the objective function are compared.The contemporary optimal solution is selected by using the Metropolis criterion.The solution space is reduced to one point when the algorithm calculates the last generation.The objective function fitness value corresponding to this point is the global optimal solution, and its corresponding parameter is the optimal value of the interval number judgment matrix weights w U .The variable superscripts in the matrix [X] denote the variables of the i-th indicator, and the subscripts denote the positions of the variables in the matrix.
[X] = In this paper, we set the initial parameters, the initial temperature T 0 , the annealing rate α, the termination temperature T f , and the maximum number of iterations G so that the number of iterations (k = 0, T k = T, f (Rbest)) is the contemporary optimal fitness value in the algorithm and f (Gbest) is the global optimal fitness value.The flowchart of the IRMO-SA algorithm is shown in Figure 4.
in the algorithm and f(Gbest) is the global optimal fitness value.The flowchart of the IRMO-SA algorithm is shown in Figure 4.

The Entropy Weighting Method used to Solve the Objective Weights
Due to the subjectivity of IAHP, it is necessary to introduce an objective weighting method to assign weights to each indicator objectively.Entropy is a measure of uncertainty [59].The information entropy value can determine the degree of dispersion of an indicator.The steps are described below.
To eliminate the influence of the data outline and unit of each assessment indicator, it is necessary to standardize the original data in the original sample matrix According to the meaning of the indicator, it can be divided into a positive indicator (the larger, the better) and a negative indicator (the smaller, the better).

Positive indicator:
Negative indicator: represents the j-th indicator of the i-th evaluation object; ij r is the value of the j-th indicator of the i-th evaluation object after standardization; is the minimum value of all evaluation objects in the column of the j-th indicator; and is the maximum value of the sample data of the j-th indicator evaluation object.

The Entropy Weighting Method Used to Solve the Objective Weights
Due to the subjectivity of IAHP, it is necessary to introduce an objective weighting method to assign weights to each indicator objectively.Entropy is a measure of uncertainty [59].The information entropy value can determine the degree of dispersion of an indicator.The steps are described below.
To eliminate the influence of the data outline and unit of each assessment indicator, it is necessary to standardize the original data in the original sample matrix X = (x ij ) m×n .According to the meaning of the indicator, it can be divided into a positive indicator (the larger, the better) and a negative indicator (the smaller, the better).
Positive indicator : Negative indicator : represents the j-th indicator of the i-th evaluation object; r ij is the value of the j-th indicator of the i-th evaluation object after standardization; x i min j is the minimum value of all evaluation objects in the column of the j-th indicator; and x i max j is the maximum value of the sample data of the j-th indicator evaluation object.
(2) Normalize the data to obtain the normalization matrix Y = (y ij ).
(3) Calculate the information entropy value e j for the j-th indicator.
where y ij = 0, y ij ln y ij = 0, and m is the number of evaluation indicators.
(4) Calculate the objective weight w sj for the j-th indicator.

Combination Weighting Method Based on Game Theory
This paper uses the AHP and EWM to conduct subjective and objective weighting assignments.The previously obtained subjective and objective weight values are combined through game theory to obtain the final weights of each indicator, which can largely minimize the bias between subjective and objective.The combination steps are as follows: (1) Construct the set of basis weight vectors w k .
where M represents a total of M different weighting assignment methods, and n represents the number of evaluation indicators.
(2) Construct a linear combination of weight vectors w.
β k is the linear combination weight coefficient.(3) Optimally solve the linear combination of weight coefficients β k .
Based on the game theory, the linear combination of weight coefficients β k is optimized to minimize the divergence between the linear combination of weight vectors w and the underlying set of weight vectors w k : According to the matrix differentiation property, a system of linear equations with optimized first-order derivative conditions can be obtained: By solving Equation ( 16), the optimized linear combination weight coefficients ) Normalize β k to obtain the optimal linear combination weight coefficients β * k .
(5) Calculate the combined weight w.
The optimal subjective weights obtained using the IAHP-IRMO method are combined and assigned with the objective weights determined using the EWM to obtain the comprehensive weights w of the assessment indicators: 3.2.TDCM

Basic Concepts
The cloud model not only realizes the two-way conversion of qualitative concepts and quantitative data [60] but also considers uncertainty problems, such as the ambiguity and randomness of the assessment system [61].The two-dimensional cloud model is obtained by synthesizing two one-dimensional clouds, comprehensively describing the stochasticity and ambiguity problems under the joint action of two factors.F is a twodimensional random function obeying normal distribution, and the cloud model composed of cloud drops drop(x i , y i , µ i ) is a two-dimensional normal cloud model, which is calculated as follows: x i , y i are the cloud droplet coordinates, Ex, Ey are the expectant values, En x , En y are the standard deviation values, P xi , P yi are the conditional cloud droplet coordinates, He x , He y are the hyper entropy values, and µ i is the degree of certainty.The smaller the expectant, the lower the evaluation level of the indicator; the smaller the standard deviation, the higher the acceptance of the result, and the higher the reliability of the assessment result; and the smaller the value of the hyper entropy, the smaller the dispersion of the cloud droplet, and the higher the accuracy of the evaluation result.

Two-Dimensional Cloud Modeling (1) Standard cloud
The evaluation indicator toughness level is equalized into four subintervals according to the interval [0, 10], where the j-th subinterval is denoted as [S min j , S max j ].The numerical features of the standard cloud are calculated to obtain the value intervals and numerical features of Levels I~IV: Ex ′ is the expectation of the standard cloud, En ′ is the standard deviation of the standard cloud, S max j is the upper limit value of the j-th interval, S min j is the lower limit value of the j-th interval, and He ′ is the hyper entropy of the standard cloud.The numerical characteristics of the evaluated standard cloud are shown in Table 4.The standard cloud diagram is shown in Figure 5.  (

2) Comprehensive evaluation cloud
This paper conducts a comprehensive evaluation of mine slopes regarding the two dimensions of safety and ecological restoration suitability.The two-dimensional eigenvalues of the evaluation indicators are selected as the two sets of basic variables of the twodimensional cloud model.The safety dimension and the ecological restoration suitability dimension of each evaluation indicator form a cloud droplet, which constitutes the safety evaluation cloud and the indicator's ecological restoration suitability evaluation cloud, known as the second-level cloud.Subsequently, the synthetic operation is performed, and the formula is shown below, where ', ', ' Ex En He are higher-level digital features: (3) Proximity Since the two-dimensional cloud map is displayed as a three-dimensional view, the spatial graphics will cause visual errors, so a more precise method is needed to determine the comprehensive evaluation grade of each mine.The proximity degree is introduced to make judgments as a means of obtaining accurate evaluation ratings, calculated using the following formula: (3) Proximity Since the two-dimensional cloud map is displayed as a three-dimensional view, the spatial graphics will cause visual errors, so a more precise method is needed to determine the comprehensive evaluation grade of each mine.The proximity degree is introduced to make judgments as a means of obtaining accurate evaluation ratings, calculated using the following formula: where N is the proximity of the composite level; ___ Ex is the expectant of the safety standard cloud and Ex is the expectant of the actual safety cloud; and ____ Ex ′ is the expectant of the ecological standard cloud and Ex ′ is the expectant of the actual ecological cloud.

Accuracy Analysis of IRMO-SA
To prove that the optimal subjective weights obtained using the IRMO-SA algorithm optimizing the search interval weight solution are more reliable and reasonable, this paper applies the IAHP-IRMO-SA method with the interval number eigenvalue method (IEM) to solve the third-, fourth-, and fifth-order judgment matrices of the literature [62][63][64].The authors of [62] used a hierarchical analysis and interval fuzzy hierarchical analysis to calculate the judgment matrix in the decision-making process for a flood risk assessment of a subway system; the authors of [63] calculated the interval number judgment matrix based on nonlinear fuzzy mathematics and a hierarchical analysis to evaluate the risk of sudden water on the coal mine floor; and the authors of [64] optimized the acquired research data for quantitative analyzation based on a hierarchical analysis, interval hierarchical analysis, and the interval number theory.As can be seen in Tables 5-7, the results in this paper are consistent with those in the literature.

Stability and Efficiency Verification of IRMO-SA Method
To verify the stability and efficiency of the proposed method, this paper selects the interval number judgment matrix in the literature [65] and adopts the three global optimization algorithms of RMO, IRMO, and IRMO-SA to consecutively search the interval number break matrix of the algorithm 20 times, and then it compares the results, as shown in Table 8 and Figure 6.As seen in Table 8, the maximum of the objective function fitness obtained using the IRMO-SA algorithm is 16.4865, the minimum is 16.4854, and the average is 16.4857.The standard deviation is only 0.00028, which is optimal among the three global optimization algorithms.Meanwhile, the volatility of Figure 6 indicates that the change in the results of the 20 consecutive searches of the IRMO-SA algorithm is the smallest, which means it is more stable, followed by the IRMO algorithm, and the RMO algorithm has the worst stability.Meanwhile, as seen in the volatility of Figure 6, the IRMO-SA algorithm has the smallest amplitude, which is more stable, followed by the IRMO algorithm, and the RMO algorithm is the worst.It shows that the IRMO-SA algorithm has good stability and is more advantageous when used to solve the weights of interval number judgment matrices under the IAHP.
As seen in Figure 7, the IRMO-SA algorithm initially stabilized at about 75 generations of iterations at the same number of iterations for three consecutive searches.In contrast, the RMO algorithm started to be initially stabilized at 98 generations of searches.This shows that the IRMO-SA algorithm has the fastest convergence speed, the IRMO algorithm is the second fastest, and the RMO algorithm is the slowest.The above comparison shows that the IRMO-SA algorithm has more outstanding efficiency in solving the weights of the interval number judgment matrix under the IAHP, and its convergence efficiency is more ideal.As seen in Table 8, the maximum of the objective function fitness obtained using the IRMO-SA algorithm is 16.4865, the minimum is 16.4854, and the average is 16.4857.The standard deviation is only 0.00028, which is optimal among the three global optimization algorithms.Meanwhile, the volatility of Figure 6 indicates that the change in the results of the 20 consecutive searches of the IRMO-SA algorithm is the smallest, which means it is more stable, followed by the IRMO algorithm, and the RMO algorithm has the worst stability.Meanwhile, as seen in the volatility of Figure 6, the IRMO-SA algorithm has the smallest amplitude, which is more stable, followed by the IRMO algorithm, and the RMO algorithm is the worst.It shows that the IRMO-SA algorithm has good stability and is more advantageous when used to solve the weights of interval number judgment matrices under the IAHP.
As seen in Figure 7, the IRMO-SA algorithm initially stabilized at about 75 generations of iterations at the same number of iterations for three consecutive searches.In contrast, the RMO algorithm started to be initially stabilized at 98 generations of searches.This shows that the IRMO-SA algorithm has the fastest convergence speed, the IRMO algorithm is the second fastest, and the RMO algorithm is the slowest.The above comparison shows that the IRMO-SA algorithm has more outstanding efficiency in solving the weights of the interval number judgment matrix under the IAHP, and its convergence efficiency is more ideal.
Mathematics 2024, 12, x FOR PEER REVIEW 14 of 27 number break matrix of the algorithm 20 times, and then it compares the results, as shown in Table 8 and Figure 6.As seen in Table 8, the maximum of the objective function fitness obtained using the IRMO-SA algorithm is 16.4865, the minimum is 16.4854, and the average is 16.4857.The standard deviation is only 0.00028, which is optimal among the three global optimization algorithms.Meanwhile, the volatility of Figure 6 indicates that the change in the results of the 20 consecutive searches of the IRMO-SA algorithm is the smallest, which means it is more stable, followed by the IRMO algorithm, and the RMO algorithm has the worst stability.Meanwhile, as seen in the volatility of Figure 6, the IRMO-SA algorithm has the smallest amplitude, which is more stable, followed by the IRMO algorithm, and the RMO algorithm is the worst.It shows that the IRMO-SA algorithm has good stability and is more advantageous when used to solve the weights of interval number judgment matrices under the IAHP.
As seen in Figure 7, the IRMO-SA algorithm initially stabilized at about 75 generations of iterations at the same number of iterations for three consecutive searches.In contrast, the RMO algorithm started to be initially stabilized at 98 generations of searches.This shows that the IRMO-SA algorithm has the fastest convergence speed, the IRMO algorithm is the second fastest, and the RMO algorithm is the slowest.The above comparison shows that the IRMO-SA algorithm has more outstanding efficiency in solving the weights of the interval number judgment matrix under the IAHP, and its convergence efficiency is more ideal.

Overview
This paper selects high and steep slopes in three areas of the Yellow River Basin Ecological Restoration Demonstration Project area, namely the Yanhua Quarry, Torch Quarry, and Zhubei Quarry #1, for a comprehensive resilience evaluation.The primary mining method in this area is "gourd-type" blasting at the bottom, forming high and steep traumatic surfaces with large topographic ups and downs.High and steep slopes remain from the previous mining project, with an overall slope inclination of approximately 140 • ; the land is excavated and occupied, and the vegetation is seriously destroyed.The current situation of the slopes of abandoned mines in the Yellow River Basin is shown in Figures 8 and 9.The characteristic period of the seismic reflection spectrum is 0.40 s.The peak acceleration value of ground vibration is 0.15 g, corresponding to the seismic intensity of VII levels.The magnitude is seven levels, and the comprehensive horizontal seismic coefficient is 0.035.

Overview
This paper selects high and steep slopes in three areas of the Yellow River Basin Ecological Restoration Demonstration Project area, namely the Yanhua Quarry, Torch Quarry, and Zhubei Quarry #1, for a comprehensive resilience evaluation.The primary mining method in this area is "gourd-type" blasting at the bottom, forming high and steep traumatic surfaces with large topographic ups and downs.High and steep slopes remain from the previous mining project, with an overall slope inclination of approximately 140°; the land is excavated and occupied, and the vegetation is seriously destroyed.The current situation of the slopes of abandoned mines in the Yellow River Basin is shown in Figures 8 and 9.The characteristic period of the seismic reflection spectrum is 0.40 s.The peak acceleration value of ground vibration is 0.15 g, corresponding to the seismic intensity of VII levels.The magnitude is seven levels, and the comprehensive horizontal seismic coefficient is 0.035.(1) Basic situation Table 9 lists the three quarries' basic high and steep slope conditions.Table 10 is a grading table of each quarry's ecological and environmental problems.The "Degree of Joint and Fissure Development" column describes joint and fissure development for the upper and lower sections of each quarry's rock body.The "Slope Rock Body Type" column lists the types of rock bodies found on the slope.All three quarries belong to the low mountainous and hilly area, and the compressive strength of slightly weathered and strongly weathered rocks is within the range of 50-60 MPa.

Overview
This paper selects high and steep slopes in three areas of the Yellow River Basin Ecological Restoration Demonstration Project area, namely the Yanhua Quarry, Torch Quarry, and Zhubei Quarry #1, for a comprehensive resilience evaluation.The primary mining method in this area is "gourd-type" blasting at the bottom, forming high and steep traumatic surfaces with large topographic ups and downs.High and steep slopes remain from the previous mining project, with an overall slope inclination of approximately 140°; the land is excavated and occupied, and the vegetation is seriously destroyed.The current situation of the slopes of abandoned mines in the Yellow River Basin is shown in Figures 8 and 9.The characteristic period of the seismic reflection spectrum is 0.40 s.The peak acceleration value of ground vibration is 0.15 g, corresponding to the seismic intensity of VII levels.The magnitude is seven levels, and the comprehensive horizontal seismic coefficient is 0.035.(1) Basic situation Table 9 lists the three quarries' basic high and steep slope conditions.Table 10 is a grading table of each quarry's ecological and environmental problems.The "Degree of Joint and Fissure Development" column describes joint and fissure development for the upper and lower sections of each quarry's rock body.The "Slope Rock Body Type" column lists the types of rock bodies found on the slope.All three quarries belong to the low mountainous and hilly area, and the compressive strength of slightly weathered and strongly weathered rocks is within the range of 50-60 MPa.(1) Basic situation Table 9 lists the three quarries' basic high and steep slope conditions.Table 10 is a grading table of each quarry's ecological and environmental problems.The "Degree of Joint and Fissure Development" column describes joint and fissure development for the upper and lower sections of each quarry's rock body.The "Slope Rock Body Type" column lists the types of rock bodies found on the slope.All three quarries belong to the low mountainous and hilly area, and the compressive strength of slightly weathered and strongly weathered rocks is within the range of 50-60 MPa.Hancheng City has a warm, temperate, continental, semi-arid climate, with less annual precipitation than evaporation.The winter is dry and cold with little snow, the spring is dry with little rain and wind, the summer is hot and dry with many rainstorms, and the fall cools down rapidly with cloudy and rainy days.In 2021, Hancheng City had an overall high rainfall, with more than 1100 mm of precipitation.

Determination of Indicator Weights 4.2.1. Determination of Optimal Subjective Weights using SA-RMO-IAHP
This example scoring uses the 1-9 scale method to assess safety and ecological restoration suitability.The scoring results are shown in Tables A1-A6, Table A7, Table A8, Table A9, Table A10, respectively.After constructing the judgment matrix and applying the IAHP-IRMO-SA algorithm to determine the optimal subjective weights, the optimal values of the subjective weights of safety w and ecology w' for each indicator are obtained, as shown in Figure 10.11.

Determination of Objective Weights Using EWM
The initial matrix is constructed, the entropy value of the safety evaluation indicator information is determined, and the EWM is applied to obtain the objective weights of the evaluation indicators.The quantitative values of the evaluation indicators of abandoned open-pit mine slopes and the corresponding objective weights are shown in Table 11.As shown in Figure 12, the optimal subjective weights, comprehensive objective weights have the same trend of change with the different evaluati but the magnitude of change of the objective weights is larger compared subjective weights.This is because the objective weighting method only eva dicators from the mathematical aspect, which is too absolute.In contrast, weighting method often relies on the experience of experts.The evaluatio tween the indicators can affect the indicators when the weight calculation which can show the influence degree of each indicator more objectively.Co the single subjective or objective assignment method, the comprehensive game theory combination assignment method is between the optimal subj and objective weight, and the combination assignment process considers th As shown in Figure 12, the optimal subjective weights, comprehensive weights, and objective weights have the same trend of change with the different evaluation indicators, but the magnitude of change of the objective weights is larger compared to that of the subjective weights.This is because the objective weighting method only evaluates the indicators from the mathematical aspect, which is too absolute.In contrast, the subjective weighting method often relies on the experience of experts.The evaluation process between the indicators can affect the indicators when the weight calculation is performed, which can show the influence degree of each indicator more objectively.Compared with the single subjective or objective assignment method, the comprehensive weight of the game theory combination assignment method is between the optimal subjective weight and objective weight, and the combination assignment process considers the influence of subjective human factors and combines the actual objective data, which makes the evaluation results more reasonable and reliable.

Comprehensive Evaluation and Verification
This paper comprehensively evaluates the resilience of abandoned open-pit mine slopes in terms of both safety and ecology.The stability of a slope mainly determines its safety.Ecological restoration suitability refers to the difficulty of ecologically restoring a slope.The expert scoring method was used to score the underlying indicators within each level's thresholds and calculate the evaluation cloud numerical characteristics for each underlying event concerning the natural language descriptions of each level in Table 4.This evaluation model was applied to evaluate high and steep slopes in three areas to calculate the second-level sub-safety cloud and ecological cloud, the first-level sub-safety cloud and ecological cloud digital features, and the integrated safety cloud and ecological cloud digital features.Table 12 presents the Digital characteristic parameters of cloud model for comprehensive evaluation of abandoned open-pit mine slopes.This study employed the traditional TOPSIS method [66] to validate and calculate the engineering case.Table 13 presents the results of the safety assessment and ecological suitability assessment obtained through the TOPSIS method.The ranking outcomes align with the results in this study, indicating the accuracy of the novel comprehensive evaluation method proposed for abandoned open-pit mine slope resilience in this paper.A comprehensive cloud of resilience for these slopes was generated using the MATLAB R2021b software based on the digital features of the two-dimensional evaluation cloud model from Table 12 and the standard cloud digital features from Table 4. Figure 13 illustrates the two-dimensional cloud map for the resilience assessment of the steep slopes in Zhubei Quarry #1.A preliminary analysis suggests that Zhubei Quarry #1's resilience level falls between Levels II and III.
A comprehensive cloud of resilience for these slopes was generated using the MATLAB R2021b software based on the digital features of the two-dimensional evaluation cloud model from Table 12 and the standard cloud digital features from Table 4. Figure 13 illustrates the two-dimensional cloud map for the resilience assessment of the steep slopes in Zhubei Quarry #1.A preliminary analysis suggests that Zhubei Quarry #1's resilience level falls between Levels II and III.When displayed as a three-dimensional view, using the two-dimensional cloud map to conduct a resilience assessment of abandoned open-pit mine slopes may lead to visual errors and result in overly coarse evaluation outcomes.Therefore, it is necessary to calculate proximity to more accurately characterize the resilience levels of individual indicators.Proximity is computed based on the digital features of the two-dimensional cloud model.The resilience assessment levels for steep slopes in abandoned mines can be determined accordingly, as shown in Table 14.When displayed as a three-dimensional view, using the two-dimensional cloud map to conduct a resilience assessment of abandoned open-pit mine slopes may lead to visual errors and result in overly coarse evaluation outcomes.Therefore, it is necessary to calculate proximity to more accurately characterize the resilience levels of individual indicators.Proximity is computed based on the digital features of the two-dimensional cloud model.The resilience assessment levels for steep slopes in abandoned mines can be determined accordingly, as shown in Table 14.By applying the maximum membership principle, the overall resilience assessment levels for steep slopes in the Yanhua Quarry, Torch Quarry, and Zhubei Quarry #1 are determined as follows: the Yanhua Quarry is classified as Level II, indicating that it has "relatively high resilience"; the Torch Quarry is also rated as Level II, denoting "relatively high resilience"; and Zhubei Quarry #1 is categorized as Level III, representing "moderate resilience".The resilience ranking is as follows: Zhubei Quarry #1 > Yanhua Quarry > Torch Quarry.

Discussion
This paper proposes a new method based on the TDCM with combination weighting to comprehensively evaluate slope resilience in abandoned open-pit mines.Four primary indicators (i.e., geological structural features, slope geometric characteristics, slope rock mass features, and external factors) are selected to establish the evaluation system.Fifteen secondary indicators (i.e., rock weathering degree, mine damage area, and degree of development of slope joints and fractures) are considered.The method is applied to three areas of the Yellow River Basin Ecological Restoration Demonstration Project area and compared with other methods.The results indicate that it exhibits superior stability, efficiency, and accuracy.A comparison with the actual engineering data further validates the effectiveness of the method.This section aims to further discuss the data in this paper.
(1) As indicated in Table 12, the safety ranking is as follows: Zhubei Quarry #1 > Yanhua Quarry > Torch Quarry.The ecological restoration suitability ranking is as follows: Zhubei Quarry #1 > Yanhua Quarry > Torch Quarry.By comparing the actual conditions of the mines and the data in Table 10, the safety evaluation results are generally consistent with the actual situation.The total investment in ecological restoration in Table 9 is ranked as follows: Zhubei Quarry #1 > Yanhua Quarry > Torch Quarry.This sequence indirectly reflects the ecological restoration suitability of the respective quarry slopes, aligning with the ecological assessment results obtained using the evaluation methodology employed in this study.(2) During the management process in mining areas, researchers often integrate concepts such as the comprehensive utilization of waste rock resources, the elimination of geological safety hazards, slope stability, and ecological restoration.This is achieved through practices such as slope reduction, the lowering of step slope ratios, height control, platform installation, and the creation of tiered slopes to ensure overall slope stability.The reduction in slope load and slope ratio through slope cutting can eliminate geological safety hazards and maintain slope stability.The determination of the amount of rock cutting considers the ecological restoration suitability and slope safety of the mining area, reflecting the resilience level of abandoned open-pit mine steep slopes.Based on the data in Table 9, the quantity of cutting stone ranking is as follows: Zhubei Quarry #1 > Yanhua Quarry > Torch Quarry.The evaluation results align closely with the on-site situation, validating the applicability and feasibility of the method proposed in this paper.(3) Zhubei Quarry #1 is categorized as Level III and requires intensified ecological restoration and safety management efforts.Zhubei Quarry #1 has a large area of exposed rock walls, with a significant presence of residual hills, posing challenges for ecological restoration, as illustrated in Figure 14.It is situated in a high-risk geological disaster area with severe vegetation damage.Mining activities have profoundly impacted its topography and landscape, resulting in significant soil erosion.It needs to undergo slope-cutting construction, and vegetation needs to be planted for slope surface recovery.

Conclusions
This paper proposed a new method for comprehensively evaluating slope resilience in abandoned open-pit mines.A comprehensive evaluation indicator system was established by considering the two dimensions of safety and ecological restoration suitability.Then, the IRMO-SA algorithm was employed to optimize the IAHP to enhance the calculation stability and efficiency when solving subjective weights.The TDCM was used to obtain the resilience levels.The details are as follows:

Conclusions
This paper proposed a new method for comprehensively evaluating slope resilience in abandoned open-pit mines.A comprehensive evaluation indicator system was established by considering the two dimensions of safety and ecological restoration suitability.Then, the IRMO-SA algorithm was employed to optimize the IAHP to enhance the calculation stability and efficiency when solving subjective weights.The TDCM was used to obtain the resilience levels.The details are as follows: (1) By integrating the resilience theory, this paper analyzed the safety and ecological restoration suitability of abandoned open-pit mine slopes.A comprehensive evaluation model was proposed to evaluate the safety and suitability for ecological restoration of abandoned open-pit mine slopes.(2) By integrating the IAHP, EWM, and the enhanced IRMO-SA algorithm, a combined weighting method was established.The results indicate that the proposed method can improve the objectivity and rationality of the evaluation and increase the calculation stability compared to traditional evaluation approaches.(3) The resilience of mine slopes was categorized into risk and ecological dimensions, and the TDCM was introduced.The cloud model was utilized to quantify fuzziness and uncertainty, obtain the resilience level of the mine slopes, and visualize the evaluation results.(4) The developed resilience assessment model was applied to evaluate steep slopes in three areas: the Yanhua Quarry, Torch Quarry, and Zhubei Quarry #1.The results are consistent with the actual engineering situations, demonstrating the model's applicability and feasibility.(5) This paper proposed a novel method for assessing slope resilience in abandoned open-pit mines.This method offers theoretical support for the safety management and ecological restoration of mine slopes.It assists construction personnel in taking targeted actions based on evaluation results to ensure the structural safety of abandoned mine slopes and protect their ecological environment.It also provides valuable guidance for reinforcement and ecological management projects on such slopes.

Figure 2 .
Figure 2. Flow chart of comprehensive evaluation of abandoned open-pit mine slopes.

Figure 2 .
Figure 2. Flow chart of comprehensive evaluation of abandoned open-pit mine slopes.

D
(a, b)-The degree of divergence between the sum of interval numbers.When b l = b r , D(a, b) is the degree of divergence from point b to interval a.The objective optimization model is established by combining the concept of the degree of divergence between interval numbers: ) is the objective function in the objective optimization model; w U is the optimal subjective weight of each assessment indicator; W ij denotes the judgment range when indicator i and indicator j are compared; and W ij = [ Ui and w Uj are the optimal subjective weights of indicator i and indicator j, respectively, and D(A ij , W ij ) is the degree of divergence of intervals A ij and W ij .When the optimal subjective weight is taken for both indicator i and indicator j, D(A ij , w Ui w Uj ) denotes the degree of divergence of the point w Ui w Uj to interval A ij .

Figure 4 .
Figure 4.The flowchart of the IRMO-SA algorithm.

Figure 4 .
Figure 4.The flowchart of the IRMO-SA algorithm.

Table 4 . 1 ) 27 '
Risk evaluation criteria cloud digital characterization.Mathematics 2024, 12, x FOR PEER REVIEW 12 of Ex is the expectation of the standard cloud, ' En is the standard deviation of the standard cloud, m ax j S is the upper limit value of the j-th interval, m in j S is the lower limit value of the j-th interval, and ' He is the hyper entropy of the standard cloud.The nu- merical characteristics of the evaluated standard cloud are shown in Table 4.The standard cloud diagram is shown in Figure 5.

( 2 )
Comprehensive evaluation cloudThis paper conducts a comprehensive evaluation of mine slopes regarding the two dimensions of safety and ecological restoration suitability.The two-dimensional eigenvalues of the evaluation indicators are selected as the two sets of basic variables of the two-dimensional cloud model.The safety dimension and the ecological restoration suitability dimension of each evaluation indicator form a cloud droplet, which constitutes the safety evaluation cloud and the indicator's ecological restoration suitability evaluation cloud, known as the second-level cloud.Subsequently, the synthetic operation is performed, and the formula is shown below, where Ex ′ , En ′ , He ′ are higher-level digital features:

Figure 6 .
Figure 6.Results of 20 searches of objective function fitness value by different algorithms.

Figure 6 .
Figure 6.Results of 20 searches of objective function fitness value by different algorithms.

Figure 6 .
Figure 6.Results of 20 searches of objective function fitness value by different algorithms.

Figure 8 .Figure 9 .
Figure 8. Aerial photograph of abandoned open-pit mine slopes in Yellow River Basin: (a) Yanhua Quarry; (b) Torch Quarry; (c) Zhubei Quarry #1.(The quarry area is shown within the red curve on the figures.)

Figure 8 .
Figure 8. Aerial photograph of abandoned open-pit mine slopes in Yellow River Basin: (a) Yanhua Quarry; (b) Torch Quarry; (c) Zhubei Quarry #1.(The quarry area is shown within the red curve on the figures.)

Figure 8 .
Figure 8. Aerial photograph of abandoned open-pit mine slopes in Yellow River Basin: (a) Yanhua Quarry; (b) Torch Quarry; (c) Zhubei Quarry #1.(The quarry area is shown within the red curve on the figures.)

Figure 10 .
Figure 10.Optimal values of subjective weights of comprehensive evaluation indicators.(a) Safety subjective weights; (b) ecological subjective weights.

Figure 10 .
Figure 10.Optimal values of subjective weights of comprehensive evaluation indicators.(a) Safety subjective weights; (b) ecological subjective weights.

11 .Figure 11 .
Figure 11.Weights of indicators in index layer: (a) safety weights; (b) ecological restoration suitability weights.The comparison of the security weights of the evaluation indicators at the indicator level is shown in Figure 12.

Figure 11 .
Figure 11.Weights of indicators in index layer: (a) safety weights; (b) ecological restoration suitability weights.The comparison of the security weights of the evaluation indicators at the indicator level is shown in Figure 12.

Figure 12 .
Figure 12.Comparison of weights of indicators to assess security of indicators.

Figure 13 .
Figure 13.Graph showing evaluation results of TDCM: (a) overall evaluation; (b) criterion layer evaluation; (c) index layer A evaluation; (d) index layer B evaluation; (e) index layer C evaluation; (f) index layer D evaluation.

Figure 13 .
Figure 13.Graph showing evaluation results of TDCM: (a) overall evaluation; (b) criterion layer evaluation; (c) index layer A evaluation; (d) index layer B evaluation; (e) index layer C evaluation; (f) index layer D evaluation.

Figure 14 .
Figure 14.Distribution of exposed rock walls and residual hills in Zhubei Quarry #1.

( 1 )
By integrating the resilience theory, this paper analyzed the safety and ecological restoration suitability of abandoned open-pit mine slopes.A comprehensive evaluation model was proposed to evaluate the safety and suitability for ecological restoration of abandoned open-pit mine slopes.(2)By integrating the IAHP, EWM, and the enhanced IRMO-SA algorithm, a combined weighting method was established.The results indicate that the proposed method can improve the objectivity and rationality of the evaluation and increase the calculation stability compared to traditional evaluation approaches.(3) The resilience of mine slopes was categorized into risk and ecological dimensions, and the TDCM was introduced.The cloud model was utilized to quantify fuzziness and uncertainty, obtain the resilience level of the mine slopes, and visualize the evaluation results.(4) The developed resilience assessment model was applied to evaluate steep slopes in three areas: the Yanhua Quarry, Torch Quarry, and Zhubei Quarry #1.The results are consistent with the actual engineering situations, demonstrating the model's applicability and feasibility.(5) This paper proposed a novel method for assessing slope resilience in abandoned open-pit mines.This method offers theoretical support for the safety management and ecological restoration of mine slopes.It assists construction personnel in taking targeted actions based on evaluation results to ensure the structural safety of abandoned mine slopes and protect their ecological environment.It also provides valuable guidance for reinforcement and ecological management projects on such slopes.Author Contributions: Conceptualization, L.J.; methodology, L.J.; software, P.L.; validation, P.L.; formal analysis, P.L.; investigation, P.L.; resources, P.L.; data curation, P.L.; writing-original draft preparation, P.L.; writing-review and editing, P.L. and W.Y.; visualization, P.L. and W.Y.; supervision, J.W.; project administration, P.L. and J.W. All authors have read and agreed to the published version of the manuscript.Funding: This research received no external funding.

Figure 14 .
Figure 14.Distribution of exposed rock walls and residual hills in Zhubei Quarry #1.

Table 1 .
Comprehensive evaluation indicator system of resilience in abandoned open-pit mine slopes.

Table 2 .
Grading criteria for abandoned open-pit mine slopes' comprehensive evaluation indicators.

Table 5 .
Results for third-order subjective weights.

Table 6 .
Results for fourth-order subjective weights.

Table 7 .
Results for fifth-order subjective weights.

Table 8 .
Objective function fitness values for three algorithms.

Table 8 .
Objective function fitness values for three algorithms.

Table 8 .
Objective function fitness values for three algorithms.

Table 9 .
Basic conditions of quarries.

Table 10 .
Grading table of ecological and environmental problems in quarries.

Table 11 .
Quantitative data and corresponding objective weights of indicators.

Table 11 .
Quantitative data and corresponding objective weights of indicators.

Table 12 .
Digital characteristic parameters of cloud model for comprehensive evaluation of abandoned open-pit mine slopes.

Table 14 .
Calculation of proximity of slopes.

Table 14 .
Calculation of proximity of slopes.

Table A7 .
Judgment matrix for number of ecological subjective weight intervals for criterion layer A.