Study on extracting method of single slope surface shape based on DEM:taking Wanzhou district of three gorges reservoir area as example

This paper focused on the extracting method of single slope shape. Applying the software of ARCGIS 9.0, DEM (digital elevation model) was established. From the DEM, topographic characters, including valleys and ridges were extracted. Some valley lines were extended in order to intersect with the ridges nearby. All slope were divided into different slope surface, which enclosed by valleys and ridges. And the slope surface shapes were defined three types, Line Slope, Upper-concave and lower-convex slope, Upper-convex and lower-concave slope, according to their functions. And the judge formula of single slope surface shape was brought forward. Taking Wanzhou District as test area, it indicates that auto-extracting method of single slope surface shape has high precision relatively. This study can provide references to the studies of region geological disasters prevention and cure.

Mathematical Model of the Identical Slope Surface

The formation of the identical slope surface and the method of construction are discussed. On the basement of building the parameter equation of variable-radius circle family envelope, the frequently used parameter equation of the identical slope surface of the top of taper moving along column helix, horizental arc and line is built. The equation can be used to construct the identical slope surface’s contours, gradient lines and three dimensional figures correctly.

Simulation on slope uncertainty derived from DEMs at different resolution levels:a case study in the Loess Plateau

Slope is one of the crucial terrain variables in spatial analysis and land use planning, especially in the Loess Plateau area of China which is suffering from serious soil erosion. DEM based slope extracting method has been widely accepted and applied in practice. However slope accuracy derived from this method usually does not match with its popularity. A quantitative simulation to slope data uncertainty is important not only theoretically but also necessarily to applications. This paper focuses on how resolution and terrain complexity impact on the accuracy of mean slope extracted from DEMs of different resolutions in the Loess Plateau of China. Six typical geomorphologic areas are selected as test areas, representing different terrain types from smooth to rough. Their DEMs are produced from digitizing contours of 1:10,000 scale topographic maps. Field survey results show that 5 m should be the most suitable grid size for representing slope in the Loess Plateau area. Comparative and math-simulation methodology was employed for data processing and analysis. A linear correlativity between mean slope and DEM resolution was found at all test areas, but their regression coefficients related closely with the terrain complexity of the test areas. If taking stream channel density to represent terrain complexity, mean slope error could be regressed against DEM resolution (X) and stream channel density (S) at 8 resolution levels and expressed as (0.00158+0.031S-0.0325)X-0.0045S2-0.155S+0.1625, with a R2 value of over 0.98. Practical tests also show an effective result of this model in applications. The new development methodology applied in this study should be helpful to similar researches in spatial data uncertainty investigation.

Methods for the construction of DEMs of artificial slopes considering morphological features and semantic information

As human activities increase,artificially modified terrain is increasingly widely distributed in road,hydrological,and urban construction.Artificially modified terrain plays an important role in protecting from geological disasters and in the planning and design of urban landscapes.Compared with natural slopes,artificial slopes have obvious morphological characteristics.Traditional modeling methods are no longer suitable for digital elevation model(DEM)modeling of artificial slopes because they often seriously distort the DEM results.In this paper,from the perspective of morphological characteristics,artificial slopes are divided into two types,namely,regular slopes and irregular slopes,based on whether the top and bottom lines of the artificial slope are parallel.Then,according to the morphological characteristics of the two types of slopes,the following DEM construction methods are designed:the first method(perpendicular+inverse distance weighted)is suitable for regular slopes,and the second method(perpendicular+high-accuracy surface modeling)is suitable for irregular slopes.Finally,a DEM construction test is carried out using the artificial slopes in the study area.The results show that for the regular and irregular slopes in the study area,the construction method proposed in this paper has significant advantages in morphological accuracy over the traditional method(triangulated irregular network),and the elevation accuracy method is also superior to the traditional method(using this method,the mean error and standard deviation error of the regular slope DEM are 0.08 m and 0.13 m,respectively,and those of the irregular slope DEM are 0.08 m and 0.06 m).In addition,the top lines and bottom lines can be included in the DEM construction of the background area after processing the elevation information of the boundary line to realize a smooth transition in the boundary between the artificial slope and the background area.

Comparison of Slope Length Factor Extraction in Hillslope Soil Erosion Model with Different DEM Resolutions

In this study, non-cumulative slope length(NCSL) calculation method and spatial analytical calculation(SAC) method were respectively applied to extract slope length and slope length factor from 10 sample areas, which are located in Ansai County, north Shaanxi Province. The comparison of computation precision between variable DEM resolutions showed that NCSL was superior to SAC entirely. And the results were best when the DEM resolutions were 5 and 10 m. Besides, the results of slope length factor were nearly the same under the two conditions. So DEM of 10 m resolution can be used to extract slope length.

Slope deformation partitioning and monitoring points optimization based on cluster analysis

The scientific and fair positioning of monitoring locations for surface displacement on slopes is a prerequisite for early warning and forecasting.However,there is no specific provision on how to effectively determine the number and location of monitoring points according to the actual deformation characteristics of the slope.There are still some defects in the layout of monitoring points.To this end,based on displacement data series and spatial location information of surface displacement monitoring points,by combining displacement series correlation and spatial distance influence factors,a spatial deformation correlation calculation model of slope based on clustering analysis was proposed to calculate the correlation between different monitoring points,based on which the deformation area of the slope was divided.The redundant monitoring points in each partition were eliminated based on the partition's outcome,and the overall optimal arrangement of slope monitoring points was then achieved.This method scientifically addresses the issues of slope deformation zoning and data gathering overlap.It not only eliminates human subjectivity from slope deformation zoning but also increases the efficiency and accuracy of slope monitoring.In order to verify the effectiveness of the method,a sand-mudstone interbedded CounterTilt excavation slope in the Chongqing city of China was used as the research object.Twenty-four monitoring points deployed on this slope were monitored for surface displacement for 13 months.The spatial location of the monitoring points was discussed.The results show that the proposed method of slope deformation zoning and the optimized placement of monitoring points are feasible.

Probabilistic assessment of rock slopes stability using the response surface approach——A case study

This paper examines the stability condition of a jointed rock slope in the south western region of Saudi Arabia using deterministic and probabilistic approaches,under both dry and wet conditions.The study area is characterized by complex geology in rugged terrains.The stability analysis is carried out using the code FLAC3D to generate a 3-dimensional,ubiquitous joint model,to determine the influence of the dominant,unfavourable discontinuity orientation with respect to the slope face.The deterministic analysis is first implemented using the mean values of the selected random variables,namely the dip,dip direction and friction angle of the dominant discontinuity set,and the stability condition is assessed with a factor of safety based on the classical frictional joint constitutive model.A Box-Behnken design(BBD)approach is then adopted to create the surface response function as a second order polynomial for the factor of safety.To do so,fifteen FLAC3D models are generated in accordance with the BBD.Based on this,10,000 simulations of different slope realizations are carried out using Monte-Carlo simulation technique,and the probability of unsatisfactory of performance of the rock slope is assessed.It is shown that the probabilistic approach provides more insight and confidence in the stability condition of the rock slope,both under dry and steady state heavy rainfall conditions.A discussion is presented on the significance of accepting lower safety factors when heavy rainfall conditions are encountered.

A Method Combining Numerical Analysis and Limit Equilibrium Theory to Determine Potential Slip Surfaces in Soil Slopes

This paper describes a precise method combining numerical analysis and limit equilibrium theory to determine potential slip surfaces in soil slopes. In this method, the direction of the critical slip surface at any point in a slope is determined using the Coulomb’s strength principle and the extremum principle based on the ratio of the shear strength to the shear stress at that point. The ratio, which is considered as an analysis index, can be computed once the stress field of the soil slope is obtained. The critical slip direction at any point in the slope must be the tangential direction of a potential slip surface passing through the point. Therefore, starting from a point on the top of the slope surface or on the horizontal segment outside the slope toe, the increment with a small distance into the slope is used to choose another point and the corresponding slip direction at the point is computed. Connecting all the points used in the computation forms a potential slip surface exiting at the starting point. Then the factor of safety for any potential slip surface can be computed using limit equilibrium method like Spencer method. After factors of safety for all the potential slip surfaces are obtained, the minimum one is the factor of safety for the slope and the corresponding potential slip surface is the critical slip surface of the slope. The proposed method does not need to pre-assume the shape of potential slip surfaces. Thus it is suitable for any shape of slip surfaces. Moreover the method is very simple to be applied. Examples are presented in this paper to illustrate the feasibility of the proposed method programmed in ANSYS software by macro commands.

A Comparative Visualization Strategy of Probabilistic Susceptibility Maps for Analyzing Different Types of Slope Failures

This paper presents a comparative visualization strategy of slope failure susceptibility maps for analyzing different types of simultaneous occurrences of slope failures. Through the SEM （structural equation modeling）, slope failure susceptibility maps are produced by using causal factors （i.e., geographical information, satellite remotely sensed data）. As for a conventional pair-wise comparative procedure, the differences between susceptibility maps are delineated on difference maps, that can be, however, applied for evaluating differences only between pairs of susceptibility maps. One of the strong requirements from specialists working on slope stability evaluation is a comparative and visualization strategy of susceptibility maps with respect to ＂different types of simultaneous slope failures＂, for which the discussion is insufficient in the previous research activities for constructing the quantitative models for slope failure hazard mapping. As a measure, a color composite map based on susceptibility maps has been produced. The combination of assigning susceptibility maps to RGB-color planes is determined based on an index of ＂NCCT （normalized correlated color temperature）＂ which represents the relationship between chromaticity and human visual perception. Through the cases examined, the result indicates that the proposed color composite map, as a heuristic visualization strategy, is useful for simultaneously evaluating the hazardous areas affected by ＂different types of slope failures＂.

Comprehensive Analysis Method of Slope Stability Based on the Limit Equilibrium and Finite Element Methods and Its Application

To study the safety and stability of large slopes, taking the right side slope of the Yuxi’an tunnel of the Yuchu Expressway Bridge in Yunnan Province as an example, limit equilibrium and finite element analysis were applied to engineering examples to calculate the stability coefficient of the slope before and after excavation in the natural state. After comparative analysis, it was concluded that the former had a clear mechanical model and concept, which could quickly provide stability results;the latter could accurately determine the sliding surface of the slope and simulate the stress state changes of the rock and soil mass. The stability coefficients calculated by the two methods were within the stable range, but their values were different. On this basis, combined with the calculation principles, advantages and disadvantages of the two methods, a comprehensive analysis method of slope stability based on the limit equilibrium and finite element methods was proposed, and the rationality of the stability coefficient calculated by this method was judged for a slope case.

Characteristics of runoff processes and nitrogen loss via surface flow and interflow from weathered granite slopes of Southeast China

Rainfall intensity and slope gradient are two of the most important factors affecting the variations of runoff nitrogen(N).However,the effects of slope gradient and rainfall intensity on N loss via surface flow and interflow on weathered granite slopes are poorly understood.In this study,12 artificial rainfalls(three rainfall intensities and four slope gradients)were simulated to investigate the coupling loss characteristics of surface flow–interflow–total nitrogen(TN),nitrate nitrogen(NO_3^--N)and ammonia nitrogen(NH_4^+-N)on weathered granite slopes.The results show that slope gradient has a greater impact on the surface flow when the rainfall intensity is relatively large.The effect gradually weakens with the decrement of rainfall intensity.The interflow yield increases firstly with the prolongation of rainfall duration,then tends to be stable and finally decreases.The total surface flow percentage increases with rainfall intensity while it decreases with increasing slope gradient with a range of 10.88%-71.47%.The TN loss concentration of the surface flow continually decreases with rainfall duration while that of the interflow shows different fluctuations.However,the TN loss loads of both surface flow and interflow increase with increasing rainfall intensity and slope gradient.The NO_3^--N concentration of interflow is much higher than that of the surface flow.The NH_4^+-N concentration is always less than that of NO_3^--N with no significant difference between surface flow and interflow.The percentages of the TN,NO_3^--N,and NH_4^+-N total loss load and concentration of surface flow and interflow were analyzed.The results show that N loss via both surface flow and interflow occurs mainly in the form of NO_3^--N.Most of the N loss is caused by interflow which is the preferential path of runoff nutrient loss.These findings provide data support and underlying insights for the control of runoff and N loss on the weathered granite slopes.

Reliability analysis of slopes considering spatial variability of soil properties based on efficiently identified representative slip surfaces

Slope reliability analysis considering inherent spatial variability(ISV)of soil properties is timeconsuming when response surface method(RSM)is used,because of the"curse of dimensionality".This paper proposes an effective method for identification of representative slip surfaces(RSSs)of slopes with spatially varied soils within the framework of limit equilibrium method(LEM),which utilizes an adaptive K-means clustering approach.Then,an improved slope reliability analysis based on the RSSs and RSM considering soil spatial variability,in perspective of computation efficiency,is established.The detailed implementation procedure of the proposed method is well documented,and the ability of the method in identifying RSSs and estimating reliability is investigated via three slope examples.Results show that the proposed method can automatically identify the RSSs of slope with only one evaluation of the conventional deterministic slope stability model.The RSSs are invariant with the statistics of soil properties,which allows parametric studies that are often required in slope reliability analysis to be efficiently achieved with ease.It is also found that the proposed method provides comparable values of factor of safety(FS)and probability of failure(Pf)of slopes with those obtained from direct analysis and lite rature.

Analysis of Earthquake-Triggered Failure Mechanisms of Slopes and Sliding Surfaces

Earthquake-induced landslides along the Dujiangyan-Yingxiu highway after the Ms 8.0 Wenchuan earthquake in 2008 were investigated. It was found that： （1） slopes were shattered and damaged during the earthquake and open tension cracks formed on the tops of the slopes; （2） the upper parts of slopes collapsed and slid, while the lower parts remained basically intact, indicating that the upper parts of slopes would be damaged more heavily than the lower parts during an earthquake. Large-scale shaking table model tests were conducted to study failure behavior of slopes under the Wenchuan seismic wave, which reproduced the process of deformation and failure of slopes. Tension cracks emerged at the top and upper part of model, while the bottom of the model remained intact, consistent with field investigations. Depth of the tension crack at the top of model is 32 cm, i.e., 3.2 m compared to the prototype natural slope with a height of 14 m when the length scale ratio （proto/model） is lo. Acceleration at the top of the slope was almost twice as large as that at the toe when the measured accelerations on shaking table are 4.85 m/s2 and 6.49 m/s2, which means that seismic force at the top of the slope is twice the magnitude of that at the toe. By use of the dynamic-strength-reduction method, numerical simulation was conducted to explore the process and mechanism of formation of the sliding surface, with other quantified information. The earthquake-induced failure surfaces commonly consist of tension cracks and shear zones. Within 5 mfrom the top of the slope, the dynamic sliding surface will be about 1 m shallower than the pseudo-static sliding surface in a horizontal direction when the peak ground acceleration （PGA） is 1 m/s2; the dynamic sliding surface will be about 2 m deeper than the pseudo-static sliding surface in a horizontal direction when the PGA is lo m/sL and the depths of the dynamic sliding surface and the pseudo-static sliding surface will be almost the same when the PGA is 2 m/s2. Based on these findings, it is suggested that the key point of anti-seismic design, as well as for mitigation of post-earthquake, secondary mountain hazards, is to prevent tension cracks from forming in the upper part of the slope. Therefore, the depth of tension cracks in slope surfaces is the key to reinforcement of slopes. The depth of the sliding surface from the pseudo-static method can be a reference for slope reinforcement mitigation.

A Simple Monte Carlo Method for Locating the Three-dimensional Critical Slip Surface of a Slope

Based on the assumption of the plain-strain problem, various optimization or random search methods have been developed for locating the critical slip surfaces in slope-stability analysis, but none of such methods is applicable to the 3D case. In this paper, a simple Monte Carlo random simulation method is proposed to identify the 3D critical slip surface. Assuming the initial slip to be the lower part of a slip ellipsoid, the 3D critical slip surface is located by means of a minimized 3D safety factor. A column-based 3D slope stability analysis model is used to calculate this factor. In this study, some practical cases of known minimum safety factors and critical slip surfaces in 2D analysis are extended to 3D slope problems to locate the critical slip surfaces. Compared with the 2D result, the resulting 3D critical slip surface has no apparent difference in terms of only cross section, but the associated 3D safety factor is definitely higher.

Statistical Distribution of Surface Slope in A 3-D Ocean Wave Field

A joint probability density function (PDF) for surface slopes in two arbitrary directions is derived on the basis of Longuet Higgins's linear model for three-dimensionol (3-D) random wave field. and the correlation moments of surface slopes. as parameters in the PDF, are expressed in terms of directional spectrum of ocean waves. So long as the directional spectrum model is given, these parameters are determined. Since the directional spectrum models proposed so far are mostly parameterized by the wind speed and fetch, this allows for substituting these parameters with thc wind speed and fetch. As an example, the wind speed and fetch are taken to be 14 m ' s and 200 km, and the Hasselmann and Donclan directional spectra are, respectively, use to compute these parameters. Some novel results a reobtained. One of the increasing interesting results is that the variances of surface slope in downwind and cross-wind directions determined by the Donclan directional spectra are close to those measured by Cox and Munk (1954). Some discussions are made on these results.

3D Identification and Stability Analysis of Key Surface Blocks of Rock Slope

Complicated geological structures make it difficult to analyze the stability of rock slopes, such as faults, weak intercalated layers or joint fissures. Based on 3D geological modeling and surface block identifying methods, an integrated methodology framework was proposed and realized to analyze the stability of surface blocks in rock slopes. The surface blocks cut by geological structures, fissures or free faces could be identified subjected to the four principles of closure, completeness, uniqueness and validity. The factor of safety(FOS)of single key block was calculated by the limit equilibrium method. If there were two or more connected blocks, they were defined as a block-group. The FOS of a block-group was computed by the Sarma method. The proposed approach was applied to an actual rock slope of a hydropower project, and some possible instable blocks were demonstrated and analyzed visually. The obtained results on the key blocks or block-groups provide essential information for determining potential instable region of rock slopes and designing effective support scheme in advance.

Comparative Study on Response Surfaces for Reliability Analysis of Spatially Variable Soil Slope

This paper focuses on the performance of the second-order polynomial-based response surfaces on the reliability of spatially variable soil slope. A single response surface constructed to approximate the slope system failure performance function G（X） （called single RS） and multiple response surfaces constructed on finite number of slip surfaces （called multiple RS） are developed, respectively. Single RS and multiple RS are applied to evaluate the system failure probability pf for a cohesive soil slope together with Monte Carlo simulation （MCS）. It is found thatpe calculated by single RS deviates significantly from that obtained by searching a large number of potential slip surfaces, and this deviation becomes insignificant with the decrease of the number of random variables or the increase of the scale of fluctuation. In other words, single RS cannot approximate G（X） with reasonable accuracy. The value of pc from multiple response surfaces fits well with that obtained by searching a large number of potential slip surfaces. That is, multiple RS can estimate G（X） with reasonable accuracy.

Comparative study of mining methods for reserves beneath end slope in flat surface mines with ultra-thick coal seams

The paper aims to identify a reasonable method for mining ultra-thick coal seams in an end-slope in surface mine, With a case study of Heidaigou surface coal mine(HSCM), the paper conducted a comparative research on three mining methods, namely Underground Mining Method(UMM), Highwall Mining System(HMS) and Local Steep Slope Mining Method(LSSMM). A model was firstly established to simulate the impact that UMM and HMS exert on monitoring points and surface deformation. The way that stripping and excavation amount varies with different slope angle, and the corresponding end slope stability were analyzed in the mode of LSSMM. Then a TOPSIS model was established by taking into account six indicators such as recovery ratio, technical complexity and adaptability, the impact on surface mining production, production safety and economic benefits. Finally, LSSMM was determined as the best mining method for mining ultra-thick coal seams in end slope in HSCM.

Static and dynamic analysis on slope stability using a DFN-DEM approach on the right abutment of the Karun 4 dam

Uncertainty in input fracture geometric parameters during analysis of the stability of jointed rock slopes is inevitable and therefore the stochastic discrete fracture network(DFN)-distinct element method(DEM) is an efficient modeling tool. In this research, potentially unstable conditions are detected in the right abutment of the Karun 4 dam and downstream of the dam body as a case study. Two extreme states with small and relatively large block sizes are selected and a series of numerical DEM models are generated using a number of validated DFN models. Stability of the rock slope is assessed in both static and dynamic loading states. Based on the design basis earthquake(DBE) and maximum credible earthquake(MCE) expected in the dam site, histories of seismic waves are applied to analyze the stability of the slope in dynamic earthquake conditions. The results indicate that a MCE is likely to trigger sliding of rock blocks on the rock slope major joint. Furthermore, the dynamic analysis also shows a local block failure by the DBE, which can consequently lead to slope instability over the long term. According to the seismic behavior of the two models, larger blocks are prone to greater instability and are less safe against earthquakes.

Seismic stability analysis of slopes with pre-existing slip surfaces

In analyzing seismic stability of a slope with upper bound limit analysis method, the slip surface is often assumed as a log-spiral or plane slip surface. However, due to the presence of a weak layer and unfavorable geological structural surface or a bedrock interface with overlying soft strata, the preexisting slip surface of the slope may be irregular and composed of a series of planes rather than strictly logspiral or plane shape. A computational model is developed for analyzing the seismic stability of slopes with pre-existing slip surfaces. This model is based on the upper bound limit analysis method and can consider the effect of anchor bolts. The soil or rock is deemed to follow the Mohr-Coulomb yield criterion. The slope is divided into multiple block elements along the slip surface. According to the displacement compatibility and the associated flow rule, a kinematic velocity field of the slope can be obtained computationally. The proposed model allows not only calculation of the rate of external work owing to the combined effect of self-weight and seismic loading, but also that of the energy dissipation rate caused by the slip surface, interfaces of block elements and anchorage effect of the anchors. Considering a direct relationship between the rate of external work and the energy dissipation rate, the expressions of yield acceleration and permanent displacement of anchored slopes can be derived. Finally, the validity of this proposed model is illustrated by analysis on three typical slopes. The results showed that the proposed model is more easily formulated and does not need to solve complex equations or time consuming iterations compared with previous methods based on the conditions of force equilibrium.