Murree, located in the sub-Himalayan Mountains, with geographical position extending from 33°52′ to 33°59′ North and 73°24′ to 73°31′ East is the worst slide-affected area of Pakistan. This area is recurrently affected by landslides and causes severe damages to land, life lines, houses, livestock and human life. This article has tried to determine the socio-economic determinants of landslide risk perception in Murree hills of Pakistan. Information was collected from 200, randomly selected sample households through questionnaire survey. The questionnaire consisted of questions regarding socio-economic characteristics of the respondents and risk perceptions about landslides. Socio-economic variables affecting landslide risks were also determined through relevant literature. A binary logistic regression was used to find out determinants of landslide risk perception in the study area. Results revealed that out of five variables, three variables including educational level, location and past experience had significant impact on landslide risk perception. The study recommends important policy implications for preparedness and mitigation plans against landslides in the study area.
Citation: Said Qasim, Muhammad Qasim, Rajendra P. Shrestha, Amir Nawaz Khan. Socio-economic determinants of landslide risk perception in Murree hills of Pakistan[J]. AIMS Environmental Science, 2018, 5(5): 305-314. doi: 10.3934/environsci.2018.5.305
Related Papers:
[1]
Ma Jun, Han Xinyu, Xu Qian, Chen Shiyou, Zhao Wenbo, Li Xiaoke .
Reliability-based EDM process parameter optimization using kriging model and sequential sampling. Mathematical Biosciences and Engineering, 2019, 16(6): 7421-7432.
doi: 10.3934/mbe.2019371
[2]
Charles Wiseman, M.D. .
Questions from the fourth son:
A clinician reflects on immunomonitoring, surrogate markers and systems biology. Mathematical Biosciences and Engineering, 2011, 8(2): 279-287.
doi: 10.3934/mbe.2011.8.279
[3]
Dongning Chen, Jianchang Liu, Chengyu Yao, Ziwei Zhang, Xinwei Du .
Multi-strategy improved salp swarm algorithm and its application in reliability optimization. Mathematical Biosciences and Engineering, 2022, 19(5): 5269-5292.
doi: 10.3934/mbe.2022247
[4]
Zulqurnain Sabir, Hafiz Abdul Wahab, Juan L.G. Guirao .
A novel design of Gudermannian function as a neural network for the singular nonlinear delayed, prediction and pantograph differential models. Mathematical Biosciences and Engineering, 2022, 19(1): 663-687.
doi: 10.3934/mbe.2022030
[5]
Zulqurnain Sabir, Muhammad Asif Zahoor Raja, Aldawoud Kamal, Juan L.G. Guirao, Dac-Nhuong Le, Tareq Saeed, Mohamad Salama .
Neuro-Swarm heuristic using interior-point algorithm to solve a third kind of multi-singular nonlinear system. Mathematical Biosciences and Engineering, 2021, 18(5): 5285-5308.
doi: 10.3934/mbe.2021268
[6]
Konstantin Weise, Erik Müller, Lucas Poßner, Thomas R. Knösche .
Comparison of the performance and reliability between improved sampling strategies for polynomial chaos expansion. Mathematical Biosciences and Engineering, 2022, 19(8): 7425-7480.
doi: 10.3934/mbe.2022351
[7]
Xuewei Wang, Xiaohu Zhao, Fei Li, Qiang Lin, Zhenghui Hu .
Sample entropy and surrogate data analysis for Alzheimer’s disease. Mathematical Biosciences and Engineering, 2019, 16(6): 6892-6906.
doi: 10.3934/mbe.2019345
[8]
Zulqurnain Sabir, Muhammad Asif Zahoor Raja, Abeer S. Alnahdi, Mdi Begum Jeelani, M. A. Abdelkawy .
Numerical investigations of the nonlinear smoke model using the Gudermannian neural networks. Mathematical Biosciences and Engineering, 2022, 19(1): 351-370.
doi: 10.3934/mbe.2022018
[9]
Nawaz Ali Zardari, Razali Ngah, Omar Hayat, Ali Hassan Sodhro .
Adaptive mobility-aware and reliable routing protocols for healthcare vehicular network. Mathematical Biosciences and Engineering, 2022, 19(7): 7156-7177.
doi: 10.3934/mbe.2022338
[10]
Xiaoye Zhao, Yinlan Gong, Lihua Xu, Ling Xia, Jucheng Zhang, Dingchang Zheng, Zongbi Yao, Xinjie Zhang, Haicheng Wei, Jun Jiang, Haipeng Liu, Jiandong Mao .
Entropy-based reliable non-invasive detection of coronary microvascular dysfunction using machine learning algorithm. Mathematical Biosciences and Engineering, 2023, 20(7): 13061-13085.
doi: 10.3934/mbe.2023582
Abstract
Murree, located in the sub-Himalayan Mountains, with geographical position extending from 33°52′ to 33°59′ North and 73°24′ to 73°31′ East is the worst slide-affected area of Pakistan. This area is recurrently affected by landslides and causes severe damages to land, life lines, houses, livestock and human life. This article has tried to determine the socio-economic determinants of landslide risk perception in Murree hills of Pakistan. Information was collected from 200, randomly selected sample households through questionnaire survey. The questionnaire consisted of questions regarding socio-economic characteristics of the respondents and risk perceptions about landslides. Socio-economic variables affecting landslide risks were also determined through relevant literature. A binary logistic regression was used to find out determinants of landslide risk perception in the study area. Results revealed that out of five variables, three variables including educational level, location and past experience had significant impact on landslide risk perception. The study recommends important policy implications for preparedness and mitigation plans against landslides in the study area.
β: Target reliability; X: Random variable; μX: The mean ofX; d: Deterministic design variable; GMSE: Generalized mean square cross-validation error; SORA: Sequential optimization and reliability assessment; SVM: Support vector machine; DDV: Deterministic design variable; SVR: Support vector regression; MCC: Multiple correlation coefficient; RBDO: Reliability-based design optimization; MCS: Monte Carlo simulation; DRM: Dimensionality reduction methods; RBF: Radial basis functions; RAAE: Relative average absolute error; LHD: Latin hypercube design; RMSE: Root mean squared error; FFD: Full factorial design; RRMSE: Relative root mean squared error; KKT: Karush Kuhn Tucker; RMAE: Relative maximum absolute error; VF: Variable fidelityFORM: First-order reliability methods; LF: Low fidelity; CBS: Constraint boundary sampling; HF: High fidelity; EGRA: Efficient global reliability analysis; LAS: Local adaptive sampling; EGO: Efficient global optimization; MPP: Most possible point; PRS: Polynomial response surfaces; RIA: Reliability index approach; EFF: Expected feasibility function; PMA: Performance measure approach
1.
Reliability-based design optimization
To improve the production efficiency and reduce manufacturing cost, design optimization of structure or system has become an inevitable trend in modern product. Due to the various uncertainty in engineering, the traditional deterministic design optimization method cannot ensure the practical application effect. On the other hand, in methods based on safety factor, the selection of safety factor depends on the manufacturing capacity and quality control means of the enterprise, which does not provide a quantitative measure of the safety margin and is not quantitatively linked to the influence of different design variables [1,2]. Furthermore, with emergence of new processes, new methods and new production conditions, the safety factor is difficult to determine, which also makes it difficult to cope with the complex environment of design, manufacture and use. RBDO is a natural extension of deterministic optimization method, where various uncertainty such as material properties, load conditions, manufactures, assembly and surrounding environment are considered. By considering the random distribution of various factors that affect the output response (stress, strain, vibration, fatigue, etc.) in the design optimization stage, RBDO can adapt to the complex environment in practical engineering, so that the production efficiency and manufacturing cost can be balanced.
1.1. RBDO model
In RBDO, the optimal solution that meets the minimum reliability level is searched. The RBDO model is defined as follows [3,4]:
In Eq. (1), d is the deterministic design variable with the upper bound dU and lower bound dL; X is the random variable, whose mean is μX with the upper bound μUX and lower bound μLX; f(d,μX) is the objective function; gi(d,X) is the performance function, whose value less than 0 is regarded as a failure event; βti is the target reliability; Φ(−βti) is the maximum allowable failure probability.
Because of the uncertainty of random variable X, the output response of the performance function at any point is also uncertain, so the failure probability P(gi(d,X)≤0) can be defined as follows:
FG(g)=P(gi(d,X)≤0)=∫gi(X)≤0⋯∫fX(X)dX
(2)
In Eq. (2), fX(X) is the joint probability density function of random variable X. As can be seen, Eq. (2) is a multi-dimensional integral problem, which is difficult to solve in engineering. Therefore, a variety of reliability analysis methods are developed to avoid solving the multi-dimensional integral problem directly [5,6,7].
1.2. Reliability analysis
Reliability analysis is an important part of RBDO, which is used to evaluate the reliability level of current design. Commonly used reliability analysis methods can be divided into traditional analytical methods [8,9] and numerical simulation methods [10,11]. The traditional analytic method calculates the failure probability by expanding the performance function at a specific point in a Taylor series or approximating the probability density function in Eq. (2).
Traditional analytical methods include methods based on probability density function approximation [12] and methods based on most probable (failure) point, MPP [13]. The former evaluates the probability density of performance function by assuming the output distribution type. The latter calculates the failure probability by expanding the performance function at a specific point in a Taylor series. MPP based methods include first-order reliability methods (FORM) [14,15], second-order reliability methods (SORM) [16,17,18] and MPP point-based dimensionality reduction methods (DRM) [19]. The traditional analytical methods are efficient, but for complex performance function, the calculation accuracy is poor.
Different from traditional analytical methods, the numerical simulation methods evaluate the responses of a large number of samples, then the failure probability is calculated by the ratio of failure sample number to total sample number [20]. With the characteristics of easy realization and high accuracy [21], Monte Carlo Simulation (MCS) is the most commonly used numerical simulation method, which is not affected by distribution of random variables and types of performance functions. But accurate reliability analysis often requires a lot of samples, and the response values at these samples usually need computer simulation or real experiments to obtain. Hence, the computational cost of MCS is usually very high [22]. To improve the computational efficiency of MCS, the important sampling (IS) strategy can be adopted. By relocating the random sampling center to increase the probability of the sample point falling into the failure region, the IS strategy reduces the number of sample points needed to calculate the failure probability, thus the computational costs can be reduced heavily [23,24]. Besides the IS strategy, other methods to improve the efficiency of MCS include Line Sampling [25,26], Directional Simulation [27,28] and Subset Simulation [29,30] et al.
In surrogate-assisted RBDO, MCS is most commonly used to calculate the failure probability and corresponding gradient. In MCS, a group of random samples are generated according to the probability distribution of random variables, and then their responses are evaluated to judge whether the failure event occurs. The ratio of the failure number to the total number is regarded as the failure probability. Using MCS to calculate the failure probability is as follows:
The test pointsxj,i=1,⋅⋅⋅,N in Eq. (4) are the same as that used for estimating the failure probability P(g(X)≤0). In other words, the calculation of ∂Pj(μkX)∂μX does not require additional performance function evaluation. Details about the Eq. (4) are in references [31]. It's worth noting that, if the random variable obeys other distribution type, Rosenblatt transformation [32] or Nataf transformation [33] can be used to transform them to the normal distribution.
1.3. Coupling treatment of reliability analysis and design optimization
From the mathematical model of RBDO in Eq. (1), it can be seen that the solution process of RBDO is divided into two stages: reliability analysis stage to calculate whether the failure probability of the current design point meets the reliability requirement and design optimization stage to find the optimal solution that meets the probability constraint. To handle the coupling relationship between the two stages, three types of methods are formed: double loop methods, single loop methods and decoupling methods. The characteristics of these three methods are compared in Table 1. Generally, the double loop methods are easy to implement, but due to the multiplication of the number of iterations in both design optimization and reliability analysis loops, the computational cost of double loop methods is very high. In single loop methods, the reliability analysis loop is approximated by its Karush–Kuhn–Tucker (KKT) optimality conditions. Therefore, the efficiency of single-loop methods is very high for linear and moderate nonlinear performance functions. But due to the large error caused by approximation, the single loop methods have bad performance for highly nonlinear problem. In decoupling methods, reliability analysis and design optimization are performed sequentially, therefore these methods can achieve the balance of accuracy and efficiency.
Table 1.
Comparison of double loop methods, single loop methods and decoupling methods.
Reliability index approach (RIA) [34,35] and performance measure approach (PMA) [36,37,38] are the most commonly used double loop methods. Because reliability analysis is needed at each iteration point of design optimization, the calculation efficiency of the double loop methods is very low [39]. Using the Karush Kuhn Tucker (KKT) condition instead of the probabilistic constraint to directly solve the RBDO problem is the most commonly used single loop methods [40,41,42]. Although the single loop method has high efficiency, Aoues shows that the stability of the single loop strategy is insufficient through standard examples [43].
Different from the double loop methods and single loop methods, the decoupling methods use the deterministic constraints to approximate the probability constraint sequentially. Li proposed one of the earliest decoupling strategies to solve RBDO problem by linear approximation of reliability index [47]. Then some scholars introduced MPP [48] and developed approximate decoupling strategies such as adaptive sequential linear programming algorithms [49,50] to balance accuracy, efficiency and convergence. In addition, Wu put forward the equivalent probabilistic constraint based on safety factor, which simplifies the solving process of performance function value and improves the solving efficiency [51]. Du proposed the well-known sequential optimization and reliability assessment (SORA) method, which decouples the RBDO process into a sequential single loop process to improve the efficiency of probabilistic design [52]. Other decoupling methods in RBDO include the convex linear method [55], saddle point approximation method [56], optimal shifting vector [4] etc.
In addition, to alleviate the computational time caused by solution in both physical space and normalized space, hybrid method (HM) [57,58], improved hybrid method (IHM) [59,60], optimum safety factor (OSF) method [61,62], robust hybrid method (RHM) [61,62,63,64,65] and hybrid modified method [66] are developed. These methods solved the optimization problem and the reliability problem simultaneously in the hybrid design space (HDS) [63]. Therefore, the huge calculation time can be reduced.
1.4. Optimization algorithms in RBDO
In reliability analysis or design optimization, various optimization algorithms are used to calculate the MPP and design point. Commonly used optimization algorithms include sequential approximate programming (SAP), active set method, trust region algorithm, penalty function method and so on.
In addition, evolutionary methods and swarm intelligence algorithms are also used in RBDO. Garakani [67] used adaptive Metropolis algorithm to obtain data and combined multi-gene genetic programming (MGGP) model with MCS to calculate the failure probability. Yadav established the robust design optimization (RDO) and RBDO model of trusses, and then solved them by particle swarm optimization (PSO) algorithm [68]. Tong [69] proposed a new hybrid reliability analysis algorithm based on PSO-optimized Kriging model. By comparing the calculation results of PSO-optimized Kriging model with the original model, it is verified that this method can greatly reduce the number of sample points and improve the solution efficiency. Aiming at minimizing maintenance cost under reliability constraints, Chen proposed a collaborative PSO-Dynamic Programming (co PSO-DP) method for multi-stage and multi scheme maintenance [70]. In order to reduce the weight and cost of crane metal structure (CMS), Fan established a discrete double loop RBDO model [71]. In the inner loop, Latin hypercube sampling (LHS) and MCS are used to analyze the reliability, and ant colony optimization (ACO) algorithm is used in the outer loop to solve the optimization results satisfying the deterministic and probabilistic constraints. Antos [72] combined ant colony algorithm with first-order reliability method to form nested optimization cycle, and applied it to the RBDO of geosynthetics reinforced earth retaining wall. Okasha used the improved firefly algorithm to solve the RBDO problem of truss structures with discrete design variables, which can obtain a reasonable solution with a reasonable amount of computations [73].
In evolutionary methods and swarm intelligence algorithms, there are no special requirements on the mathematical properties (such as convexity, continuity or explicit definition) of objective function and constraints in the optimization process. Therefore, evolutionary methods and swarm intelligence algorithms may be the research hotspot of RBDO in the future.
2.
Research on surrogate models
The evaluation of implicit performance function is an important part of reliability analysis and RBDO. If the real response value of implicit performance function is called directly during reliability analysis and design optimization, the computational cost is unacceptable. Therefore, the introduction of surrogate models in reliability analysis and RBDO has become a research hotspot.
2.1. Commonly used Surrogate Models
Constructed by interpolation or fitting of representative sample points, surrogate models have the characteristics of simple structure and high computational efficiency. Therefore, surrogate models have been widely used to approximate the implicit performance function in engineering. Most commonly used surrogate models include polynomial response surfaces (PRS) [36], Kriging approximation models [74,75,76], radial basis functions (RBF) [77], and support vector machine (SVM) or support vector regression (SVR).
Polynomial response surface model is also called polynomial regression model or polynomial approximation model, which is simple and easy to implement. Thus, any implicit function can be expressed by polynomial approximation. Commonly used second order polynomial response surface model is formulated as follows [36,78]:
ˆg(x)=a0+∑ni=1aixi+∑ni=1∑nj≥1aijxixj
(5)
In Eq. (5), xi is the ith component of the n dimension design variable; a0 is the constant term; ai is the coefficient of first term; aij is the coefficient of quadratic term (also known as cross term).
Kriging model consists of a polynomial response function model and a Gaussian random process, which is formulated as follows:
ˆy(x)=F(β,x)+Z(x)
(6)
In Eq. (6), F(β,x) is the global polynomial response surface model with weight coefficientβ; Z(x) is a nonparametric random process model with 0 as mean value and σ2 as variance [79,80]. In Kriging model, an indicator of a priori uncertainty prediction can be provided, thus it is widely used in engineering problems [81].
Taking the Euclidean distance between the predicted points and the sample point as the design variable, RBF model is constructed by linear combination of radial symmetric kernel functions. The RBF model is expressed as follows [77,82]:
∧y(x)=∑nk=1ωkψk(‖x−xk‖)
(7)
In Eq. (7), n is the number of sample points; ωk is the weight coefficient; ψ(‖x−xk‖) is the basis function centered on sample xk; ‖x−xk‖ is the Euclidean distance between the predicted points and the sample points. For highly nonlinear problems, good fitting performance can usually be observed using RBF model [77].
Support vector machine (SVM) is a supervised learning model based on structural risk minimization criteria. SVR model is a specific form of support vector machine in the prediction field. SVR realizes the approximate of the low-dimensional spatial problem by mapping the parameter combination vector to the high-dimensional space and constructing the regression function in the high-dimensional space [83]. SVR model is expressed as follows:
∧p(x)=〈w⋅Φ(x,xk)〉+b
(8)
In Eq. (8), 〈w⋅Φ(x,xk)〉 is the inner product of hyperparameterwand support vector regression kernel functionΦ(x,xk); bis the undetermined coefficient.
In addition to the above models, some achievements are also made in variable fidelity (VF) approximation models and ensemble approximation models. Variable fidelity model (also termed as multi-fidelity model) is an approximate model constructed by fusion of high and low fidelity samples. In VF model, the global approximation trend of implicit function is reflected by low fidelity model, and the accuracy is guaranteed by introducing a small number of high-fidelity samples [84,85,86,87]. As shown in Figure 1, the variable fidelity model fv(x) shows good fitting performance to the high-fidelity model fh(x) by fusing the low fidelity model fl(x) and four high fidelity samples.
The VF model is expected to fuse the advantages of LF and HF models and thus become a new choice in engineering. Nevertheless, how to determine the variable fidelity modeling framework and how to quantitatively compare the computational cost of HF model, LF model and VF model still need to be researched [88].
Ensemble model is another choice which expects to combines the advantage of multiple single surrogate models [89,90,91,92]. The modeling strategies of ensemble model can be divided into two categories: weighting and voting [90]. The former is based on the prediction accuracy of single models, and a larger weight coefficient is assigned to the model with higher accuracy. The weight coefficients of each single model in the ensemble model is added up to 1 to ensure the accuracy of the prediction value of the ensemble model [89,92]. Voting strategy can be regarded as a special weighting strategy, which is usually used in the case of multiple iterations in modeling or optimization. The core is to set the weight of the model with the highest precision to 1, and the weights of other models to 0. In single surrogate models, the modeling performance usually relies on the selection of samples. But the ensemble model can avoid these shortcomings, thus it is able to achieve a robust modeling effect [90].
2.2. Selection of samples
Construction of surrogate model depends on the selection of samples, which can be divided into one-shot sampling methods and sequential sampling methods [94]. Commonly used one-shot sampling methods include full factorial design (FFD) [94,95], orthogonal design [96,97], uniform design [98,99], Latin hypercube design (LHD) [100,101], etc. Sequential sampling method is to establish the initial surrogate model by one-shot sampling first, and then update the surrogate model by adding samples until the required accuracy is achieved.
In full factorial design, all levels of each factor are combined to carry out the design of experiments. However, with the increase of design variables and levels, the number of samples required will increase exponentially, which cannot be acceptable in engineering. Therefore, this method is only suitable for the case of fewer design factors and levels. In order to reduce the computational cost of full factorial design, some representative samples instead of all FFD samples can be selected in orthogonal design or uniform design. Furthermore, Latin hypercube design can also be adopted to ensure the representativeness of samples and reduce the repeatability of samples through the layered treatment in the whole design space. Figure 2 shows the sample distribution of 2-factor-3-level full factor experiment design (9 samples), 3-factor-3-level orthogonal experiment design (9 samples) and 9-point Latin hypercube design.
Figure 2.
Sample distribution of one-shot sampling methods.
One-shot sampling methods is easy to implement, but the accuracy of the surrogate model cannot be guaranteed. Increasing the number of samples will increase the computational cost and even may cause overfitting. On the contrary, fewer samples will affect the accuracy of surrogate model. Therefore, it is necessary to add sequential samples to update the surrogate model in reliability analysis and RBDO. Sequential sampling methods is a strategy which improves the model accuracy by dividing the modeling process into several stages. The flowchart of sequential sampling methods is shown in Figure 3. First, the initial surrogate model is built by one-shot sampling, and then the surrogate model is updated by adding samples according to some criteria until the required accuracy is achieved.
Figure 4 further describes the process of sequential sampling. f(x) is the true performance function which needs to be approximated. f0ˆ(x) is the initial surrogate model constructed by initial samples (black circles) generated from one-shot sampling method. Using the sequential samples (black boxes) to update the surrogate model, the final surrogate mode fsˆ(x) with high approximating accuracy is acquired.
The selection of sequential samples can be divided into following methods: error-based methods, methods based on distance to existing samples and methods of sampling near constrained boundaries (Seen in Table 2). The error-based methods can guarantee the accuracy of the surrogate model, but it is difficult to accurately evaluate the modeling error with limit samples. The method based on the distance to the existing samples can evenly distribute the sample points in the design domain. But for highly nonlinear or discontinuous problems, these methods may have poor performance. Sampling near the constraint boundary can improve the fitting accuracy of the constraint boundary, but if the constraint boundary is approximated inaccurately, many samples will be wasted. In engineering, these three sampling methods are often combined to obtain the comprehensive advantages [102,103,104,105].
Table 2.
Sequential sampling methods.
Methods
Error-based methods
Methods based on distance to existing samples
Methods of sampling near constrained boundaries
Advantages
Guarantee the modeling precision
Ensure the uniformity of sample distribution
Ensure the fitting accuracy of the constrained boundary
Disadvantages
Accurate error is difficult to assess
Poor effect for highly nonlinear or discontinuous problem
Sampling at inaccurate constraint boundary may cause waste of samples
The accuracy of the surrogate model reflects the fitting degree between the predicted response function and the real response function. The following error indexes can be used in accuracy evaluation of surrogate model [116]:
(1) Multiple correlation coefficient R2
R2 is used to evaluate the relative error of the surrogate model in the whole design space [117]. The larger the value of R2, the higher the accuracy of the surrogate model is.
R2=1−∑mi=1(yi−ˆyi)2∑mi=1(yi−ˉyi)2=1−MSEVar
(9)
(2) Relative Average Absolute Error (RAAE)
RAAE is used to evaluate the global relative error of the surrogate model. The smaller the value of RAAE, the smaller the global relative error is.
RAAE=∑mi=1|yi−ˆyi|m×STD
(10)
(3) Root Mean Squared Error (RMSE)
RMSE is used to evaluate the closeness between the predicted value and the true value [118]. The smaller the value of RMSE, the higher the accuracy of the surrogate model is.
RMSE=√{∑mi=1(yi−ˆyi)2}/m
(11)
(4) Relative Root Mean Squared Error (RRMSE)
RRMSE is used to indicate the degree of accuracy improvement of the surrogate model. The smaller the value of RRMSE, the higher the degree of accuracy improvement of the surrogate model is.
RRMSE=1STD√∑mi=1(yi−ˆyi)2m
(12)
(5) Relative Maximum Absolute Error (RMAE)
RMAE is used to evaluate the local prediction error of the surrogate model. The smaller the RMAE value, the smaller the local prediction error of the surrogate model is.
RMAE=maxi=1,⋯,m|yi−ˆyi|STD
(13)
(6) Generalized Mean Square Cross-validation Error (GMSE)
GMSE is used to evaluate the prediction error of the surrogate model with limited samples. The smaller the value of GMSE, the smaller the prediction error of the surrogate model is.
GMSE=1m∑mj=1(fj−ˆfj(−j))2
(14)
In equations (9)-(14), m is the number of samples used for model accuracy evaluation; yi is the actual response value of performance function, and ˆy is the predicted value obtained by the surrogate model; ˉy=1mm∑yiyi; ˆfj(−j) is the predicted value at point x(j) by the surrogate model constructed using the remainingm−1sample points except point (x(j),fj). MSE is the mean square error of the experimental data, Var is the variance of the experimental data, and STD is the standard deviation of the experimental data. MSE, Var and STD are calculated as follows:
MSE=∑mi=1(yi−ˆyi)2m,Var=∑mi=1(yi−ˉy)2m,STD=√Var
(15)
3.
Surrogate-assisted RBDO methods
In recent decades, many researches have been performed in surrogate-assisted RBDO methods, especially in sequential sampling based RBDO methods. The sequential sampling methods for surrogate model are reviewed in previous section. Different from the sequential sampling methods in surrogate model, some characteristics of RBDO solution process are considered. For example, sequential sampling methods on improving the accuracy of failure probability calculation [103,119], using MPP for sequential sampling [120], using the samples in the key area of reliability analysis or RBDO [122] and using MCS samples to build surrogate model [22] are also considered in surrogate-assisted RBDO.
The sequential sampling methods in surrogate-assisted RBDO can be divided into global modeling methods and local modeling methods. It's worth noting that, the global modeling methods can be used in RBDO or other optimization problems, but the local modeling methods can only be used in RBDO. In global modeling methods, accurate performance function boundaries are first approximated in the whole design space (Figure 5(a)), and then RBDO is conducted to find the optimal design. Instead, accurate fitting is only needed in the key area of RBDO (usually near the iterative design point) in local modeling methods (Figure 5(b)). The core of local modeling methods is to update the surrogate model and carry out the design optimization sequentially, as shown in Figure 5 (b)[102].
Figure 5.
Sample distribution of global modeling method and local modeling method [102].
The sample distribution of global modeling methods and local modeling methods are compared in Figure 5. The dotted line is the real performance function curve f(x1,x2); the solid line is the prediction curve of the surrogate model ˆf(x1,x2); the red "*" is the sample used to construct the surrogate model, and the blue "+" is the optimization iterative design point of RBDO; the green ellipse in Figure 5 (b) is the local sampling region. It can be seen that the global modeling methods fit the boundary of the whole performance function accurately, while the local modeling methods only fit the key area of the optimization iteration.
3.1. RBDO methods based on global modeling
The global modeling methods in RBDO can be divided into the global modeling methods based on one-shot sampling and the global modeling methods based on sequential sampling. Using one-shot sampling, it is difficult to balance the solution precision and computational cost, so current global modeling methods are mostly based on sequential sampling. The basic flowchart of global modeling method in RBDO is shown in Figure 6. After the initial sample set and design point are determined, the initial surrogate model is constructed. And then the surrogate model is updated by global sequential sampling until the accuracy satisfies the requirement. After that, the accurate global surrogate model is used to replace the implicit objective function and performance function in RBDO. Finally, the RBDO problem is solved by double loop methods, single loop methods or decoupling methods.
Figure 6.
The flowchart of global modeling methods in RBDO.
Following describes the commonly used global modeling methods: Constraint boundary sampling (CBS). Constraint boundary sampling is developed to realize the accurate approximation of the constraint boundary in the whole design space [113]. New samples are mainly selected on the approximate constraint boundary. When the accuracy of the global surrogate model satisfies the requirement, RBDO is carried out. The sequential sampling criteria in CBS are as follows:
In CBS, samples locating near the predicted constrained boundary will be preferred.
3.2. RBDO method based on local modeling
The global modeling method can fit the constraint boundary in the whole design space well, but there are always some inactive samples in RBDO (Such as the samples located on the lower left quarter in Figure 5(a)), which will inevitably cause waste of computational resources. Therefore, the local modeling methods to fit the key region of design space is proposed. The flowchart of local modeling methods is shown in Figure 7. Different from global sampling methods, the sequential samples are mainly selected in a small local region around the current design which is calculated by double loop methods, single loop methods or decoupling methods. Therefore, the RBDO process and sequential sampling are conducted sequentially.
Figure 7.
The flowchart of local modeling methods in RBDO.
Following is a brief introduction of the representative local modeling method LAS (local adaptive sampling method). In LAS, the samples are preferentially selected near the constraint boundaries around the current design point [102]. The size of local sampling region is determined according to the target reliability and nonlinearity of the performance functions, which is shown as:
Where nc is used to quantify the nonlinearity of the performance functions; ∇ˆgi(x) is the predicted gradient value at x for the ith probabilistic constraint. (x1,⋯,xM) are the M testing points evenly located within the βt-sphere region. Nis the number of probabilistic constraints. Seen from Eq. (17), a larger local sampling region is assigned to performance function with larger target reliability and higher nonlinearity.
After the local sampling region is determined, MSE sampling criterion and CBS criterion are introduced to select the sequential samples. The detailed sequential sampling criterion in LAS is as follows:
In Eq. (18), if no performance function boundary is contained in the local sampling region, the prediction error and distance to the existing samples are used to select new samples. Otherwise, CBS sampling criterion in Eq. (16) is adopted.
3.3. Numerical example
To compare the performance of global modeling method and local modeling method, a two-dimensional RBDO problem with three probabilistic constraints is introduced, which is formulated as follows:
Figure 8 shows the locations of sample points when using CBS and LAS to solve this problem. Black "o" is the initial samples using LHS method and red "x" is the sequential samples; blue "+" is the design point. d0 is the initial design point and dopt is the optimal design point. As can be seen, both the CBS method and LAS method can obtain the optimal design dopt. The surrogate model using CBS method has good fit performance in the whole design space and most sequential sample points are located near the constraint boundary (Seen in Figure 8(a)). Differently, most of the sequential sample points in LAS are concentrated in the small local region around current iterative design points (Blue ellipse in Figure 8(b)). Therefore, the fitting capability of surrogate model is weak in the upper right quarter of the design space. Because sequential sampling is avoided in the meaningless region of RBDO process, LAS has higher efficiency. It's worth noting that, if the target reliability changes, the sequential sampling process of LAS method will be conducted again. Therefore, the surrogate model constructed by LAS cannot be reused (Seen in Figure 8(c)).
Figure 8.
The distribution of sample points for Numerical example.
Table 3 shows the comparison of global modeling methods and local modeling methods. In global modeling methods, an accurate global surrogate model is built in the whole design space, therefore this model can also be used when the target reliability changes. But when the target reliability is fixed, the efficiency is usually lower than local modeling methods (Seen in Figure 8). On the other hand, the local modeling method is more efficient for problem with fixed target reliability. But when the target reliability changes, the local sampling region and corresponding sequential sample points will change, which means different surrogate models must be built for problems with different target reliability.
Table 3.
Comparison between global modeling methods and local modeling methods.
Method
Global modeling methods
Local modeling methods
Advantages
Ensure the modeling accuracy in the whole design domain; No need to remodel when the target reliability changes
High utilization rate of samples and high fitting accuracy of key areas
Disadvantages
Modeling efficiency is usually lower than local modeling methods
Remodeling is needed when the target reliability changes
In this paper, the surrogate-assisted RBDO methods are reviewed. Commonly used reliability analysis methods, RBDO methods, surrogate models and corresponding sampling and accuracy evaluation methods are summarized and discussed. Then the surrogate-assisted RBDO methods are divided into the global modeling methods and local modeling methods. The advantages and dis advantages of these two methods are compared.
(1) This paper introduces the mathematical model of RBDO, summarizes the existing reliability analysis methods and RBDO methods. Different RBDO methods are classified and compared.
(2) This paper summarizes the commonly used surrogate model, sampling method and accuracy evaluation methods of surrogate models. Different sampling methods are classified and compared.
(3) This paper divided the surrogate-assisted RBDO methods into the global modeling methods and local modeling methods. Using a representative numerical example, the advantages and disadvantages of these two modeling methods are compared.
Although many surrogate-assisted RBDO methods have been developed, there are still some research space to further strengthen the application of these methods in engineering. Future researches can be focused on the following aspects:
(1) Quantification of modeling error of surrogate model
Generally, the number of samples used to build the surrogate model is finite, thus the prediction accuracy of the surrogate model cannot be well guaranteed. Meanwhile, in existing researches, the modeling error of surrogate model is mainly used in sequential sampling. This strategy can achieve the optimal RBDO solution for general functions, but it may have large error for multimodal functions, or even cannot meet the reliability requirements. Therefore, it is of great significance to quantify the modeling error and construct its conservative model as an alternative to the implicit performance function, which finally achieve the performance improvement on the premise of ensuring the reliability requirements.
(2) Determination of high-fidelity samples and low fidelity samples in variable fidelity model
Commonly used variable fidelity methods are based on the situation that the function form (often from empirical formula, etc.) of the low fidelity model is known. In view of the unknown function form of low fidelity model, it is very important to allocate the high and low fidelity samples reasonably in RBDO process. In addition, the study of the variable fidelity method based on the surrogate model will also be a useful supplement.
(3) Construction of surrogate model for high dimensional RBDO problem with limited samples
The existence of "dimension disaster" makes the computational cost of existing RBDO methods unacceptable, which not only exists in the construction of surrogate model, but also in the process of reliability analysis and design optimization. In order to acquire an accurate and efficient RBDO solution, the problem dimension must be limited. Effective variable selection methods, high dimensional model representation methods, high dimensional reliability analysis methods and high dimensional optimization methods will be future research hotspots.
(4) Construction of classification model to balance the precision and cost in RBDO
Besides using surrogate models to accurately approximate the implicit objective and performance functions, introduction of the classification method such as support vector machine is another viable option. However, the existing support vector machine method is commonly confronted with precision problem because of lack of real function response information. Therefore, it is necessary to develop a more effective classification model to balance the precision and cost in RBDO.
Acknowledgments
This research was funded by National Natural Science Foundation of China (Grant No. 51905492), Key Scientific and Technological Research Projects in Henan Province (Grant No. 202102210089, 202102210274), Innovative Research Team (in Science and Technology) in University of Henan Province (Grant No. 20IRTSTHN015).
References
[1]
Solana MC, Kilburn CR (2003) Public awareness of landslide hazards: the Barranco de Tirajana, Gran Canaria, Spain. Geomorphology 54: 39–48. doi: 10.1016/S0169-555X(03)00054-0
[2]
Nadim F, Kjekstad O, Peduzzi P, et al. (2006) Global landslide and avalanche hotspots. Landslides 3: 159–173. doi: 10.1007/s10346-006-0036-1
[3]
Ho MC, Shaw D, Lin S, et al. (2008) How do disaster characteristics influence risk perception? Risk Anal 28: 635–643. doi: 10.1111/j.1539-6924.2008.01040.x
[4]
Papathoma-Köhle M, Neuhäuser B, Ratzinger K, et al. (2007) Elements at risk as a framework for assessing the vulnerability of communities to landslides. Nat Hazard Earth Eys 7: 765–779. doi: 10.5194/nhess-7-765-2007
[5]
Kamp U, Growley BJ, Khattak GA, et al. (2008) GIS-based landslide susceptibility mapping for the 2005 Kashmir earthquake region. Geomorphology 101: 631–642. doi: 10.1016/j.geomorph.2008.03.003
[6]
Khattak GA, Owen LA, Kamp U, et al. (2010) Evolution of earthquake-triggered landslides in the Kashmir Himalaya, northern Pakistan. Geomorphology 115: 102–108. doi: 10.1016/j.geomorph.2009.09.035
[7]
Khan AN (2001) Impact of landslide hazards on housing and related socio-economic characteristics in Murree (Pakistan). Pak Econ Soc Rev: 57–74.
[8]
Archer DR, Fowler HJ (2008) Using meteorological data to forecast seasonal runoff on the River Jhelum, Pakistan. J Hydrol 361: 10–23. doi: 10.1016/j.jhydrol.2008.07.017
[9]
Farooq S, Malik M (1996) Landslide Hazard Management and Control in Pakistan-A Review: International Centre for Integrated Mountain Development (ICIMOD).
[10]
Pearce A (1987) Plan for demonstration in Tehsil Murree for improving landslide-stability by reforestation and drainage improvement. Consultant's Report to FAO/UNDP Project PAK/78/036.
[11]
Abbasi A, Khan M, Ishfaq M, et al. (2002) Slope failure and landslide mechanism in Murree area, North Pakistan. Geol Bull Univ Peshawar 35: 125–137.
[12]
Calvello M, Papa MN, Pratschke J, et al. (2016) Landslide risk perception: a case study in Southern Italy. Landslides 13: 349–360. doi: 10.1007/s10346-015-0572-7
[13]
Khan AN (1995) Landslide Hazards and Policy-Response in Pakistan: A Case Study of Murree. Commission on Science and Technology for Sustainable Development in the South: 35.
[14]
Owen LA, Kamp U, Khattak GA, et al. (2008) Landslides triggered by the 8 October 2005 Kashmir earthquake. Geomorphology 94: 1–9. doi: 10.1016/j.geomorph.2007.04.007
[15]
Khan AN, Collins AE, Qazi F (2011) Causes and extent of environmental impacts of landslide hazard in the Himalayan region: a case study of Murree, Pakistan. Nat Hazards 57: 413–434. doi: 10.1007/s11069-010-9621-7
[16]
Rahman A, Khan AN, Collins AE (2014) Analysis of landslide causes and associated damages in the Kashmir Himalayas of Pakistan. Nat Hazards 71: 803–821. doi: 10.1007/s11069-013-0918-1
[17]
EM-DAT. The OFDA/CRED International Disaster Database.
[18]
Myers CA, Slack T, Singelmann J (2008) Social vulnerability and migration in the wake of disaster: the case of Hurricanes Katrina and Rita. Popul Env 29: 271–291. doi: 10.1007/s11111-008-0072-y
[19]
Eidsvig UM, McLean A, Vangelsten BV, et al. (2014) Assessment of socioeconomic vulnerability to landslides using an indicator-based approach: methodology and case studies. Bull Eng Geol Environ 73: 307–324. doi: 10.1007/s10064-014-0571-2
[20]
Kuhlicke C, Scolobig A, Tapsell S, et al. (2011) Contextualizing social vulnerability: findings from case studies across Europe. Nat Hazards 58: 789–810. doi: 10.1007/s11069-011-9751-6
[21]
Shirley WL, Boruff BJ, Cutter SL (2012) Social vulnerability to environmental hazards. Hazards Vulnerability and Environmental Justice: Routledge. pp. 143–160.
[22]
Lin S, Shaw D, Ho MC (2008) Why are flood and landslide victims less willing to take mitigation measures than the public? Nat Hazards 44: 305–314. doi: 10.1007/s11069-007-9136-z
[23]
Siagian TH, Purhadi P, Suhartono S, et al. (2014) Social vulnerability to natural hazards in Indonesia: driving factors and policy implications. Nat Hazards 70: 1603–1617. doi: 10.1007/s11069-013-0888-3
[24]
Chen W, Cutter SL, Emrich CT, et al. (2013) Measuring social vulnerability to natural hazards in the Yangtze River Delta region, China. Int J Disaster Risk Sci 4: 169–181. doi: 10.1007/s13753-013-0018-6
[25]
Yoon DK (2012) Assessment of social vulnerability to natural disasters: a comparative study. Nat Hazards 63: 823–843. doi: 10.1007/s11069-012-0189-2
[26]
Zhou Y, Li N, Wu W, et al. (2014) Assessment of provincial social vulnerability to natural disasters in China. Nat Hazards 71: 2165–2186. doi: 10.1007/s11069-013-1003-5
[27]
Zou LL, Wei YM (2010) Driving factors for social vulnerability to coastal hazards in Southeast Asia: results from the meta-analysis. Nat Hazards 54: 901–929. doi: 10.1007/s11069-010-9513-x
[28]
Schmidtlein MC, Shafer JM, Berry M, et al. (2011) Modeled earthquake losses and social vulnerability in Charleston, South Carolina. Appl Geogr 31: 269–281. doi: 10.1016/j.apgeog.2010.06.001
[29]
Rufat S, Tate E, Burton CG, et al. (2015) Social vulnerability to floods: Review of case studies and implications for measurement. IJDRR 14: 470–486.
[30]
Kellens W, Zaalberg R, Neutens T, et al. (2011) An analysis of the public perception of flood risk on the Belgian coast. Risk Analysis31: 1055–1068.
[31]
Sjöberg L, Moen BE, Rundmo T (2004) Explaining risk perception. An evaluation of the psychometric paradigm in risk perception research 10: 665–612.
[32]
Raaijmakers R, Krywkow J, van der Veen A (2008) Flood risk perceptions and spatial multi-criteria analysis: an exploratory research for hazard mitigation. Nat Hazards 46: 307–322. doi: 10.1007/s11069-007-9189-z
[33]
Burn DH (1999) Perceptions of flood risk: a case study of the Red River flood of 1997. Water Resour Res 35: 3451–3458. doi: 10.1029/1999WR900215
[34]
Ludy J, Kondolf GM (2012) Flood risk perception in lands "protected" by 100-year levees. Nat Hazards 61: 829–842. doi: 10.1007/s11069-011-0072-6
[35]
Sjöberg L (2000) Factors in risk perception. Risk Anal 20: 1–12. doi: 10.1111/0272-4332.00001
[36]
Armaş I, Avram E (2009) Perception of flood risk in Danube Delta, Romania. Nat Hazards 50: 269–287. doi: 10.1007/s11069-008-9337-0
[37]
Government of Pakistan (GOP) (1999). 1998 district census report of Rawalpindi. Population census organization of Pakistan, Islamabad.
[38]
Yamane T (1973) Statistics: An introductory analysis.
[39]
Salvati P, Bianchi C, Fiorucci F, et al. (2014) Perception of flood and landslide risk in Italy: a preliminary analysis. Nat Hazard Earth Eys 14: 2589–2603. doi: 10.5194/nhess-14-2589-2014
[40]
Landeros-Mugica K, Urbina-Soria J, Alcántara-Ayala I (2016) The good, the bad and the ugly: on the interactions among experience, exposure and commitment with reference to landslide risk perception in México. Nat Hazards 80: 1515–1537. doi: 10.1007/s11069-015-2037-7
[41]
Sudmeier-Rieux K, Jaquet S, Derron MH, et al. (2012) A case study of coping strategies and landslides in two villages of Central-Eastern Nepal. Appl Geogr 32: 680–690. doi: 10.1016/j.apgeog.2011.07.005
[42]
Pilgrim NK (1999) Landslides, Risk and Decision making in Kinnaur District: Bridging the Gap between Science and Public Opinion. Disasters 23: 45–65. doi: 10.1111/1467-7717.00104
[43]
Larsen M (2008) Rainfall-triggered landslides, anthropogenic hazards, and mitigation strategies. Adv Geosci 14: 147–153. doi: 10.5194/adgeo-14-147-2008
[44]
Kreibich H, Thieken AH, Petrow T, et al. (2005) Flood loss reduction of private households due to building precautionary measures--lessons learned from the Elbe flood in August 2002. Nat hazard earth sys 5: 117–126. doi: 10.5194/nhess-5-117-2005
[45]
Grothmann T, Reusswig F (2006) People at risk of flooding: why some residents take precautionary action while others do not. Nat Hazards 38: 101–120. doi: 10.1007/s11069-005-8604-6
[46]
Miceli R, Sotgiu I, Settanni M (2008) Disaster preparedness and perception of flood risk: A study in an alpine valley in Italy. J Environ Psychol 28: 164–173. doi: 10.1016/j.jenvp.2007.10.006
[47]
Botzen W, Aerts J, Van Den Bergh J (2009) Dependence of flood risk perceptions on socioeconomic and objective risk factors. Water Resour Res 45.
[48]
Siegrist M, Gutscher H (2006) Flooding risks: A comparison of lay people's perceptions and expert's assessments in Switzerland. Risk Anal 26: 971–979. doi: 10.1111/j.1539-6924.2006.00792.x
This article has been cited by:
1.
Zia Ud Din Taj, Ahmad Bilal, Muhammad Awais, Shuaib Salamat, Messam Abbas, Adnan Maqsood,
Design exploration and optimization of aerodynamics and radar cross section for a fighter aircraft,
2023,
133,
12709638,
108114,
10.1016/j.ast.2023.108114
2.
Tiansheng Miao, Jiayuan Guo, Guanghua Li, He Huang,
Inversion-based identification of DNAPLs-contaminated groundwater based on surrogate model of deep convolutional neural network,
2023,
23,
1606-9749,
129,
10.2166/ws.2022.437
3.
Abdelwahid Boutemedjet, Smail Khalfallah, Mahfoudh Cerdoun, Abdelkader Benaouali,
A multi-objective optimization strategy based on combined meta-models: Application to a wind turbine,
2024,
0954-4089,
10.1177/09544089231223859
4.
Roger Ballester Claret, Ludovic Coelho, Christian Fagiano, Cédric Julien, Didier Lucor, Nicolò Fabbiane,
Reliability based optimisation of composite plates under aeroelastic constraints via adapted surrogate modelling and genetic algorithms,
2024,
347,
02638223,
118461,
10.1016/j.compstruct.2024.118461
5.
Chongjian Yang, Junle Yang, Yixiao Qin,
Research on Comparative of Multi-Surrogate Models to Optimize Complex Truss Structures,
2024,
28,
1226-7988,
2268,
10.1007/s12205-024-0196-3
6.
Zhenglei He, Mengna Hong, Hongze Zheng, Jinfeng Wang, Qingang Xiong, Yi Man,
Towards low-carbon papermaking wastewater treatment process based on Kriging surrogate predictive model,
2023,
425,
09596526,
139039,
10.1016/j.jclepro.2023.139039
7.
R. Allahvirdizadeh, A. Andersson, R. Karoumi,
Improved dynamic design method of ballasted high-speed railway bridges using surrogate-assisted reliability-based design optimization of dependent variables,
2023,
238,
09518320,
109406,
10.1016/j.ress.2023.109406
8.
Yupeng Cui, Yang Yu, Siyuan Cheng, Mingxiu Wei, Yu Pan, Zewei Dong,
Multi-performance reliability-based concept-detailed co-design of offshore structures using modified conjugate FR algorithm,
2024,
300,
00298018,
117275,
10.1016/j.oceaneng.2024.117275
9.
Debiao Meng, Hengfei Yang, Shiyuan Yang, Yuting Zhang, Abílio M.P. De Jesus, José Correia, Tiago Fazeres-Ferradosa, Wojciech Macek, Ricardo Branco, Shun-Peng Zhu,
Kriging-assisted hybrid reliability design and optimization of offshore wind turbine support structure based on a portfolio allocation strategy,
2024,
295,
00298018,
116842,
10.1016/j.oceaneng.2024.116842
Said Qasim, Muhammad Qasim, Rajendra P. Shrestha, Amir Nawaz Khan. Socio-economic determinants of landslide risk perception in Murree hills of Pakistan[J]. AIMS Environmental Science, 2018, 5(5): 305-314. doi: 10.3934/environsci.2018.5.305
Said Qasim, Muhammad Qasim, Rajendra P. Shrestha, Amir Nawaz Khan. Socio-economic determinants of landslide risk perception in Murree hills of Pakistan[J]. AIMS Environmental Science, 2018, 5(5): 305-314. doi: 10.3934/environsci.2018.5.305