Three-dimensional stability calculation method for high and large composite slopes formed by mining stope and inner dump in adjacent open pits

The 2D limit equilibrium method is widely used for slope stability analysis.However,with the advancement of dump engineering,composite slopes often exhibit significant 3D mechanical effects.Consequently,it is of significant importance to develop an effective 3D stability calculation method for composite slopes to enhance the design and stability control of open-pit slope engineering.Using the composite slope formed by the mining stope and inner dump in Baiyinhua No.1 and No.2 open-pit coal mine as a case study,this research investigates the failure mode of composite slopes and establishes spatial shape equations for the sliding mass.By integrating the shear resistance and sliding force of each row of microstrip columns onto the bottom surface of the strip corresponding to the main sliding surface,a novel 2D equivalent physical and mechanical parameters analysis method for the strips on the main sliding surface of 3D sliding masses is proposed.Subsequently,a comprehensive 3D stability calculation method for composite slopes is developed,and the quantitative relationship between the coordinated development distance and its 3D stability coefficients is examined.The analysis reveals that the failure mode of the composite slope is characterized by cutting-bedding sliding,with the arc serving as the side interface and the weak layer as the bottom interface,while the destabilization mechanism primarily involves shear failure.The spatial form equation of the sliding mass comprises an ellipsoid and weak plane equation.The analysis revealed that when the coordinated development distance is 1500 m,the error rate between the 3D stability calculation result and the 2D stability calculation result of the composite slope is less than 8%,thereby verifying the proposed analytical method of equivalent physical and mechanical parameters and the 3D stability calculation method for composite slopes.Furthermore,the3D stability coefficient of the composite slope exhibits an exponential correlation with the coordinated development distance,with the coefficient gradually decreasing as the coordinated development distance increases.These findings provide a theoretical guideline for designing similar slope shape parameters and conducting stability analysis.

Numerical three-dimensional modeling of earthen dam piping failure

A physically-based numerical three-dimensional earthen dam piping failure model is developed for homogeneous and zoned soil dams.This model is an erosion model,coupled with force/moment equilibrium analyses.Orifice flow and two-dimensional(2D)shallow water equations(SWE)are solved to simulate dam break flows at different breaching stages.Erosion rates of different soils with different construction compaction efforts are calculated using corresponding erosion formulae.The dam's real shape,soil properties,and surrounding area are programmed.Large outer 2D-SWE grids are used to control upstream and downstream hydraulic conditions and control the boundary conditions of orifice flow,and inner 2D-SWE flow is used to scour soil and perform force/moment equilibrium analyses.This model is validated using the European Commission IMPACT(Investigation of Extreme Flood Processes and Uncertainty)Test#5 in Norway,Teton Dam failure in Idaho,USA,and Quail Creek Dike failure in Utah,USA.All calculated peak outflows are within 10%errors of observed values.Simulation results show that,for a V-shaped dam like Teton Dam,a piping breach location at the abutment tends to result in a smaller peak breach outflow than the piping breach location at the dam's center;and if Teton Dam had broken from its center for internal erosion,a peak outflow of 117851 m'/s,which is 81%larger than the peak outflow of 65120 m3/s released from its right abutment,would have been released from Teton Dam.A lower piping inlet elevation tends to cause a faster/earlier piping breach than a higher piping inlet elevation.

Three-dimensional critical slip surface locating and slope stability assessment for lava lobe of Unzen volcano

Even Unzen volcano has been declared to be in a state of relative dormancy,the latest formed lava lobe No.11 now represents a potential slope failure mass based on the latest research.This paper concentrates on the stability of the lava lobe No.11 and its possible critical sliding mass.It proposes geographic information systems (GIS) based three-dimensional (3D) slope stability analysis models.It uses a 3D locating approach to identify the 3D critical slip surface and to analyze the 3D stability of the lava lobe No.11.At the same time,the new 3D approach shows the effectiveness in selecting the range of the Monte Carlo random variables and locating the critical slip surface in different parts of the lava lobe No.11.The results are very valuable for judging the stability of the lava lobe and assigning the monitoring equipments.

Rapid and robust slope failure appraisal using aerial photogrammetry and 3D slope stability models

Slope failures are an inevitable aspect of economic pit slope designs in the mining industry.Large open pit guidelines and industry standards accept up to 30%of benches in open pits to collapse provided that they are controlled and that no personnel are at risk.Rigorous ground control measures including real time monitoring systems at TARP(trigger-action-response-plan)protocols are widely utilized to prevent personnel from being exposed to slope failure risks.Technology and computing capability are rapidly evolving.Aerial photogrammetry techniques using UAV(unmanned aerial vehicle)enable geotechnical engineers and engineering geologists to work faster and more safely by removing themselves from potential line-of-fire near unstable slopes.Slope stability modelling software using limit equilibrium(LE)and finite element(FE)methods in three dimensions(3D)is also becoming more accessible,user-friendly and faster to operate.These key components enable geotechnical engineers to undertake site investigations,develop geotechnical models and assess slope stability faster and in more detail with less exposure to fall of ground hazards in the field.This paper describes the rapid and robust process utilized at BHP Limited for appraising a slope failure at an iron ore mine site in the Pilbara region of Western Australia using a combination of UAV photogrammetry and 3D slope stability models in less than a shift(i.e.less than 12 h).

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.

Three-dimensional mixed convection stagnation-point fow past a vertical surface with second-order slip velocity

This study is concerned with the three-dimensional(3D)stagnation-point for the mixed convection flow past a vertical surface considering the first-order and secondorder velocity slips.To the authors’knowledge,this is the first study presenting this very interesting analysis.Nonlinear partial differential equations for the flow problem are transformed into nonlinear ordinary differential equations(ODEs)by using appropriate similarity transformation.These ODEs with the corresponding boundary conditions are numerically solved by utilizing the bvp4c solver in MATLAB programming language.The effects of the governing parameters on the non-dimensional velocity profiles,temperature profiles,skin friction coefficients,and the local Nusselt number are presented in detail through a series of graphs and tables.Interestingly,it is reported that the reduced skin friction coefficient decreases for the assisting flow situation and increases for the opposing flow situation.The numerical computations of the present work are compared with those from other research available in specific situations,and an excellent consensus is observed.Another exciting feature for this work is the existence of dual solutions.An important remark is that the dual solutions exist for both assisting and opposing flows.A linear stability analysis is performed showing that one solution is stable and the other solution is not stable.We notice that the mixed convection and velocity slip parameters have strong effects on the flow characteristics.These effects are depicted in graphs and discussed in this paper.The obtained results show that the first-order and second-order slip parameters have a considerable effect on the flow,as well as on the heat transfer characteristics.

Methods applied in Australian industry to evaluate coal mine slope stability

Strata failure is a principal hazard in open cut coal mining as it has the potential to cause multiple fatalities.Prior to the excavation of any slope,a geotechnical assessment should review the likely slope performance,including the risk of slope failure.Controls to manage this risk to an acceptable level should accompany the geotechnical analysis.A survey of 43 practising geotechnical engineers indicated that kinematic and 2D limit equilibrium analyses were the methods most commonly applied to analyse excavated slope stability.While these methods are well established and widely applied in the broad rock engineering disciplines(e.g.civil,hard rock),a recent review of over 60 slope failures suggests these methods have limited suitability for modelling the complex failure mechanisms observed in excavated coal mine slopes.Kinematic techniques do not adequately capture the rock mass component of excavated slope failure and do not provide a geospatial location of potential failure,while 2D limit equilibrium techniques do not adequately capture the 3D mechanisms of excavated slope failures.Methods which do consider the 3D mechanisms of slope failure are summarized for industry consideration and application.

Stochastic seismic slope stability assessment using polynomial chaos expansions combined with relevance vector machine

This paper presents probabilistic assessment of seismically-induced slope displacements considering uncertainties of seismic ground motions and soil properties.A stochastic ground motion model representing both the temporal and spectral non-stationarity of earthquake shakings and a three-dimensional rotational failure mechanism are integrated to assess Newmark-type slope displacements.A new probabilistic approach that incorporates machine learning in metamodeling technique is proposed,by combining relevance vector machine with polynomial chaos expansions(RVM-PCE).Compared with other PCE methods,the proposed RVM-PCE is shown to be more effective in estimating failure probabilities.The sensitivity and relative influence of each random input parameter to the slope displacements are discussed.Finally,the fragility curves for slope displacements are established for sitespecific soil conditions and earthquake hazard levels.The results indicate that the slope displacement is more sensitive to the intensities and strong shaking durations of seismic ground motions than the frequency contents,and a critical Arias intensity that leads to the maximum annual failure probabilities can be identified by the proposed approach.

Linear spatial instability analysis in 3D boundary layers using plane-marching 3D-LPSE

It is widely accepted that a robust and efficient method to compute the linear spatial amplified rate ought to be developed in three-dimensional （3D） boundary layers to predict the transition with the e^N method, especially when the boundary layer varies significantly in the spanwise direction. The 3D-linear parabolized stability equation （3D- LPSE） approach, a 3D extension of the two-dimensional LPSE （2D-LPSE）, is developed with a plane-marching procedure for investigating the instability of a 3D boundary layer with a significant spanwise variation. The method is suitable for a full Mach number region, and is validated by computing the unstable modes in 2D and 3D boundary layers, in both global and local instability problems. The predictions are in better agreement with the ones of the direct numerical simulation （DNS） rather than a 2D-eigenvalue problem （EVP） procedure. These results suggest that the plane-marching 3D-LPSE approach is a robust, efficient, and accurate choice for the local and global instability analysis in 2D and 3D boundary layers for all free-stream Mach numbers.

Experimental investigation on the invert stability of operating railway tunnels with different drainage systems using 3D printing technology

In recent years, the invert anomalies of operating railway tunnels in water-rich areas occur frequently,which greatly affect the transportation capacity of the railway lines. Tunnel drainage system is a crucial factor to ensure the invert stability by regulating the external water pressure(EWP). By means of a threedimensional(3D) printing model, this paper experimentally investigates the deformation behavior of the invert for the tunnels with the traditional drainage system(TDS) widely used in China and its optimized drainage system(ODS) with bottom drainage function. Six test groups with a total of 110 test conditions were designed to consider the design factors and environmental factors in engineering practice,including layout of the drainage system, blockage of the drainage system and groundwater level fluctuation. It was found that there are significant differences in the water discharge, EWP and invert stability for the tunnels with the two drainage systems. Even with a dense arrangement of the external blind tubes, TDS was still difficult to eliminate the excessive EWP below the invert, which is the main cause for the invert instability. Blockage of drainage system further increased the invert uplift and aggravated the track irregularity, especially when the blockage degree is more than 50%. However, ODS can prevent these invert anomalies by reasonably controlling the EWP at tunnel bottom. Even when the groundwater level reached 60 m and the blind tubes were fully blocked, the invert stability can still be maintained and the railway track experienced a settlement of only 1.8 mm. Meanwhile, the on-site monitoring under several rainstorms further showed that the average EWP of the invert was controlled within 84 k Pa, while the maximum settlement of the track slab was only 0.92 mm, which also was in good agreement with the results of model test.

Linear stability theory with the equivalent spanwise wavenumber correction in 3D boundary layers

The prediction on small disturbance propagation in complex three-dimensional(3D) boundary layers is of great significance in transition prediction methodology, especially in the aircraft design. In this paper, the linear stability theory(LST) with the equivalent spanwise wavenumber correction(ESWC) is proposed in order to accurately predict the linear evolution of a disturbance in a kind of boundary layer flow with a vital variation in the spanwise direction. The LST with the ESWC takes not only the scale of the mean flow with the significant variation but also the wavenumber evolution of the disturbance itself. Compared with the conventional LST, the results obtained by the new method are in excellent agreement with those of the numerical simulations. The LST with the ESWC is an effective method on the prediction of the disturbance evolution in 3D boundary layers, which improves the prediction of the LST in the applications to complex 3D boundary layers greatly.

Slope Stability Considering the Top Building Load

Slope stability is one of the most important subjects of geotechnics. The slope top-loading plays a key role in the stability of slopes in hill slope areas. When the building load is too large or the point of action from the shoulder is too close, the shear stress of the slope will be significantly greater than its shear strength, resulting in reduced slope stability. Therefore, it is of great importance to study the relationship between the building load and the stability of the slope. This study aims to analyze the influence of different building loads applied at different distances on the top of the slope and deduces their effects on the slope stability. For this purpose, a three-dimensional slope model under different building loads with different distances to the slope shoulder was established using the finite-difference analysis software Flac3D. The results show that the loads applied at different distances on the top of the slope have different effects on the slope stability. The slope factor of safety (FOS) increases with the increase of the distance between the top-loading and the slope shoulder;it varies from 1.37 to 1.53 for the load P = 120 KPa, 1.27 to 1.53 for the load P = 200 KPa, and from 1.18 to 1.44 for P = 300 KPa, resulting in the decrease of the coincidence area between the load-deformation and the potential sliding surface. The slope is no longer affected by the potential risk of sliding at approximately 20 m away from the slope shoulder.

基于3Dslope的露天矿边坡稳定性计算

针对F3断层切割平庄西露天矿西北帮边坡形成三维滑体,采用二维计算方法已经不能满足计算要求,因此采用3Dslope边坡计算软件对西北帮边坡进行稳定性计算.结果表明:通过分析西北帮边坡的影响因素,确定潜在滑坡模式,计算由F3断层和弱层以及假设的侧滑面控制的三维滑体边坡稳定系数.从而确定滑体的大小,控制F3断层影响范围,采取相应的治理措施,该成果为安全生产提供技术支持.

降雨作用下抗滑桩与三维土工网垫植草护坡坡体的稳定分析

针对抗滑桩与三维土工网垫植草护坡在降雨作用下的稳定性问题,采用有限元强度折减法,利用ABAQUS建立抗滑桩与三维土工网垫植草护坡坡体的三维有限元模型,以坡底水平位移突变情况作为边坡失稳的判断标准,分别比较了不同工况下护坡安全系数随降雨时长的变化规律,以及降雨作用下抗滑桩桩位、桩长对抗滑桩与三维土工网垫植草护坡坡体的稳定性影响,并与未考虑降雨时的影响规律进行对比。结果表明:降雨时长对三维土工网垫植草工况的边坡稳定性影响较为明显;抗滑桩桩长对边坡稳定性影响较大,桩长宜取8 m;无降雨作用下将桩位置于坡中时边坡整体获得最大安全系数,降雨作用下将抗滑桩置于相对位置(L 0.4=4.8 m)时能够起到更好的抗滑效果;抗滑桩与三维土工网垫植草在降雨作用时加固护坡的效果比无降雨时更好,且优于原边坡。

西气东输管道工程FA014+2～FA014+3号桩跨越场地的边坡稳定性评价

通过野外地质调查确定研究区的区域地震活动、工程地质条件和受力边界条件 ,采集样品进行测试获取岩体的基本物理力学指标 ,在此基础上进行不同条件下 (天然、饱水、饱水 +震动力 )的边坡三维模拟定量分析 ,评价边坡的稳定性 ,提出工程建议。

基于DEM的三维微条柱法计算滑移型塌岸安全系数的研究

为了对边坡进行三维稳定性分析,在传统的二维条分法的基础上进行拓展,基于高分辨率DEM构建边坡三维模型,将二维条分法中的土条拓展为三维条柱法中的条柱,通过分析每个条柱上的受力状况,根据极限平衡理论,以GIS二次开发的方式,将三维毕肖普法引入GIS软件平台,实现了三维圆弧滑动面的构造和三维安全系数的计算,该方法实现了三维模型构建与三维安全系数计算过程的自动化和一体化。应用本方法,以金沙江溪洛渡库区老木沟滑移型塌岸体为例,进行三维稳定系数的计算,与二维计算结果相比,结果偏于安全,主要的原因在于将滑动面视为球形滑动面,因此需要在程序中加入控制性结构面或控制性节点,使其与实际情况更加接近。

边坡稳定三维分析的极限平衡方法

介绍对边坡稳定分析三维极限平衡方法的一个新解法 ,它是二维Spencer法在三维条件下的扩展。该法保证了滑坡体三个方向的静力平衡 ,同时 ,还增加了一个整体力矩平衡条件 ,行界面的条间力不再假定为水平 ,条底的剪力方向也不假定为平行于主滑平面。考题和工程实例说明了其可行性和实用性。

空间三维滑坡敏感性分区工具及其应用

对于滑坡敏感性分区目前有三种方法:定性法、统计法和基于岩土定量模型的确定性方法。定性法基于对滑坡敏感性或灾害评估的人为判断;统计法用一个来源于结合了权重因子的预测函数或指标;而确定性法,或者说是物理定量模型法以质量、能量和动量守恒定律为基础。二维确定性模型广泛用于土木工程设计,而无限边坡模型(一维)也用于滑坡灾害分区的确定性模型。文中提出了一个新的基于GIS(地理信息系统)的滑坡敏感性分区系统,这个系统可用于从复杂地形中确认可能的危险三维(3-D)滑坡体。所有与滑坡相关的空间数据(矢量或栅格数据)都被集成到这个系统中。通过把研究区域划分为边坡单元并假定初始滑动面是椭球的下半部分,并使用Monte Carlo随机搜索法,三维滑坡稳定性分析中的三维最危险滑面是三维安全系数最小的地方。使用近似方法假定有效凝聚力、有效摩擦角和三维安全系数服从正态分布,可以计算出滑坡失稳概率。3DSlopeGIS是一个计算机程序,它内嵌了GIS Developer kit(ArcObjects of ESRI)来实现GIS空间分析功能和有效的数据管理。应用此工具可以解决所有的三维边坡空间数据解问题。通过使用空间分析、数据管理和GIS的可视化功能来处理复杂的边坡数据,三维边坡稳定性问题很容易用一个友好的可视化图形界面来解决。将3DSlopeGIS系统应用到3个滑坡敏感性分区的实例中:第一个是一个城市规划项目,第二个是预测以往滑坡灾害对临近区域可能的影响,第三个则是沿着国家主干道的滑坡分区。基于足够次数的Monte Carlo模拟法,可以确认可能的最危险滑坡体。这在以往的传统边坡稳定性分析中是不可能的。

西气东输管道工程EJ065桩附近边坡稳定性分析

西气东输管道工程EJ0 6 5桩附近边坡陡峻 ,并处于强震区 ,必须进行墩位地基所在岩体和岩壁边坡稳定性评价 (包括震动条件下 ) ,以评价跨越墩位的适宜性。通过野外地质调查确定研究区的区域地震活动、工程地质条件和受力边界条件 ,采集样品进行测试获取岩体的基本物理力学指标 ,作者进行了不同条件下 (天然、饱和、震动力 )的边坡三维模拟定量分析 ,评价了边坡的稳定性 。

露天转地下开采岩体稳定性三维数值模拟

露天转地下开采围岩应力和工程尺度同时影响着地下工程和露天边坡的稳定性.通过三维数值模拟,揭示了露天边坡内地下开采采场周围和边坡的力学环境,探讨了围岩移动变形、应力分布和破坏机理,分析了边坡稳定性状况.研究表明,在扰动边坡下进行地下开采,坡脚处的局部弧形破坏区将进一步恶化,但不会影响边坡的整体稳定性.由于边坡的卸荷作用,导致采场上覆岩层成拱机制减弱,采空区覆岩存在整体垮冒的可能性.