Intermittent slipping of landslide regulated by dilatancy evolution and velocity-weakening friction law: an efficient numerical scheme

When a block of dense sandy soil moves downhill, the shear-induced soil dilatancy along the basal shear boundary produces a negative value of excess pore pressure that increases the basal frictional resistance. Dilatancy angle,Ψ, the degree to which the basal soil dilates due to the shear, normally evolves during slope failure. A study by other researchers shows that if Ψ is constant, the block of dense soil will remain stable(or unstable) sliding when the velocity-weakening rate of the basal friction coefficient of the block is small(or large) enough. Moreover, during unstable sliding processes, the block of dense soil exhibits "periodic" patterns of intermittent slipping. Here, we used a more efficient and accurate numerical scheme to revisit that study. We expanded their model by assuming Ψ evolves during slope failure. Consequently, we acquired completely different results. For instance, even though the velocity-weakening rate of the friction coefficient is fixed at the same smaller(or larger) value that those researchers use, the stable(or unstable) steady states of landslide they predict will inversely change to unstable(or stable) when Ψ decreases(or increases) with the increase of slide displacement to a value small(or large) enough. Particularly, in unstable processes, the soil block exhibits "aperiodic" styles of intermittent slipping, instead of "periodic". We found out that the stick states appearing later last longer(or shorter) in the case of decreasing(or increasing) Ψ. Moreover, because the basic states of landslides with impacts of dilatancy evolution are not steady nor periodic, traditional stability-analysis methods cannot be "directly" used to analyze the stability of such landslides. Here, we broke through this technical problem to a degree. We showed that combining a concept called "quasi-steady-state approximation" with a traditional stability-analysis technique can qualitatively predict the instability onset of the landslides. Through this study, we demonstrated that the combination of Chebyshev collocation(CC) and 4~(th)-order Runge-Kutta methods is more accurate and efficient than the numerical scheme those researchers use.     

地下采动山体变形活动的模型实验

本文在分析横山山体地质条件及因地下采空而诱发出的山体变形破坏特征的基础上,对采空区地表山体的变形破坏进行了工程地质力学模拟研究。结果表明,横山山体之下3号煤层的采动坍陷是引起山体变形破坏的直接原因;毗邻山体的电厂区地表隆起变形破坏,是因采空区上覆岩层在坍陷时所产生的向电厂方向的蠕动,致使厂区地基土遭受挤压和变形的结果;随着3号煤层采动的结束及其上覆岩层坍陷过程的完成,电厂区将会恢复到稳定状态。 通过分析,山体岩层倾角增大时有无地下采空影响的山体变形破坏的不同特征,得出了同时具有地表侵蚀临空面和地下采空临空面的山体与仅有地表侵蚀临空面之斜玻的山体变形破坏是不同的。从不同模型材料所得实验结果的一致性,验证了此模拟实验在理论上和实际中都有很好的相似性。     

Quantifying spatial scale of positive and negative terrains pattern at watershed-scale: Case in soil and water conservation region on Loess Plateau

The positive and negative terrains (P-N terrains) widely distributed across China’s Loess Plateau constitute the dual structure characteristic of loess landforms. Analysis of loess P-N terrains at the watershed scale can serve to elucidate the structural characteristics and spatial patterns of P-N terrains, which benefits a better understanding of watershed evolution and suitable scales for loess landform research. The Two-Term Local Quadrat Variance Analysis (TTLQV) is calculated as the average of the square of the difference between the block totals of all possible adjacent pairs of block size, which can be used to detect both the scale and the intensity of landscape patches (e.g., plant/animal communities and gully networks). In this study, we determined the latitudinal and longitudinal spatial scale of P-N terrain patterns within 104 uniformly distributed watersheds in our target soil and water conservation region. The results showed that TTLQV is very effective for examining the scale of P-N terrain patterns. There were apparently three types of P-N terrain pattern in latitudinal direction (i.e., Loess Tableland type, Loess Hill type, and Transitional Form between Sand and Loess type), whereas there were both lower and higher values for P-N terrain pattern scales in all loess landforms in the longitudinal direction. The P-N terrain pattern also clearly presented anisotropy, suggesting that gully networks in the main direction were well-developed while others were relatively undeveloped. In addition, the relationships between the first scales and controlling factors (i.e., gully density, nibble degree, watershed area, mean watershed slope, NDVI, precipitation, loess thickness, and loess landforms) revealed that the first scales are primarily controlled by watershed area and loess landforms. This may indicate that the current spatial pattern of P-N terrains is characterized by internal force. In selecting suitable study areas in China’ Loess Plateau, it is crucial to understand four control variables: the spatial scale of the P-N terrain pattern, the watershed area, the main direction of the watershed, and the loess landforms.     

黄土高原边坡特征与破坏规律的分区研究

为了研究黄土高原自然边坡的特征及破坏规律,按山系与水系或水系的分水岭、地貌单元、地层岩性特征等条件,将黄土高原划分为8个区:临洮—永靖区、天水—通渭区、兰州—会宁区、陇东区、靖边—安塞区、隰县—离石区、甘泉—吉县区和汾渭区.根据极限状态边坡的4个野外判别标准,测量了8个区510个自然极限状态黄土边坡断面,分区采用指数模型回归边坡坡高与坡宽的相关关系,计算各区20、50、100m坡高的边坡稳定系数和失效概率.结果表明:黄土高原的边坡特征与破坏形式具有分区特征,且南北差异性明显.临洮—永靖区边坡坡高与坡宽呈线性关系,表明该区边坡坡度不随坡高变化,边坡稳定性受内摩擦角控制;兰州—会宁区和靖边—安塞区高坡陡,低坡缓,高坡不稳定,易发生错落式滑坡;天水—通渭区、甘泉—吉县区和汾渭区高坡缓,低坡陡,稳定性计算结果显示高坡和低坡都较为稳定,但由于地层结构和地貌的特点,高边坡易发生低速蠕变型滑坡或高速远程滑坡;陇东区边坡整体上较为稳定;隰县—离石区受黏粒含量较高的Q1地层控制,高边坡稳定性较差;50m左右坡高的黄土边坡稳定性对强度指标内聚力、内摩擦角的敏感度都高,易于失稳.     

Effects of gully erosion and gully filling on soil degradation in the black soil region of Northeast China

Gully erosion has caused soil degradation and even reduced soil productivity.However,only few studies on the effects of gully erosion and artificial controlling measures on soil degradation in the Black Soil Region of Northeast China are available.Thus,this study explores the relationships between gully erosion,gully filling and soil parameters.Two sets of soil samples were collected in the field at:(1) 72 sample points in the gully erosion study area,60 sample points in the ephemeral and classical gully erosion area(3,518 m2),12 sample points in the deposition zone(443 m2),(2)10 reference points along a slope unaffected by gully erosion representing the original situation before the gully was formed.All soil samples were analyzed for gravel content(GC),soil organic matter(SOM),total nitrogen(TN),available nitrogen(AN),available phosphorus(AP),and available potassium(AK).The soil property values on unaffected slope were fitted by the polynomial curves as the reference values in no gully erosion area.The interpolated soil property values in gully eroded study area were compared with these polynomial curves,respectively,and then,changes of soil property values were analyzed.Gully erosion caused an increase in GC and a decrease in SOM,TN,AN,AP and AK.The change of GC,SOM,TN,AN,AP,AK was 8.8%,-9.04 g kg-1,-0.92 g kg-1,-62.28 mg kg-1,-29.61 mg kg-1,-79.68 mg kg-1.The soil property values in the study area were below optimal values.Thus,we concluded that gully erosion and gully filling caused both on-site and off-site soil degradation.Soil degradation area was 0.65 % of the cultivated land.In addition,it was proved that gully filling were an improper soil and water conservation measure,which seems to exacerbate the problem.Thus,it is suggested that soil where soil is deep is moved to fill the gully,and then the area around the filled gullies should be covered by grass for preventing the formation and development of the gully.     

南梁煤矿多煤层开采地下水系统演化规律

煤层开采对所在矿区地下水系统有着重要影响。以往研究单煤层开采对地下水系统的影响较多,而对多煤层的影响研究甚少,特别对于我国西部缺水矿区。以南梁煤矿为例,运用地下水系统演化理论和岩石力学模拟等,对该矿井水文地质结构、矿井涌水变化规律、矿井地下水流场演变、矿井地下水化学成分变化等方面进行了综合分析研究,重构了多煤层开采条件下南梁矿井地下水系统流动模型,初步揭示了矿井水化学成分的演化机理。研究结果表明,2^-2煤单煤层开采时,顶板导裂带的最大发育高度为42.1 m,而2^-2和3^-1煤层重复开采时则增大为83.1 m,相应地应力、位移、塑性区范围后者也比前者增大许多。这揭示出多煤层重复采动明显增大了顶板导裂带的发育高度,加剧了矿井水文地质结构变异,进一步地,导水裂隙带改变了天然地下水渗流路径,沟通了不同含水层之间的水力联系,增强了地下水流动速度和水文地球化学作用,整体扩大了地下水流动系统的规模,从降雨入渗→导裂带渗流→各煤层涌水→井底水仓排水构成了一个自然-人工复合地下水流动模式。研究成果可以为南梁煤矿的矿井水害防治及水资源高效利...     

Gully erosion regionalization of black soil area in northeastern China

Gully erosion is the frequent and main form of soil erosion in the black soil area of the northeastern China, which is one of the most important commodity grain production bases in China. It is encroaching upon the fertile farmland there. Regionalization of gully erosion can reveal the spatial distribution and regularity of the development of gully erosion. Based on the eco-geographical regional background features of the black soil area, this study combined the regionalization with influencing factors of the development of gully erosion. GIS spatial analysis, geostatistical analysis, spatial statistics, reclassification, debris polygon processing and map algebra methods were employed. As a result, the black soil area was divided into 12 subregions. The field survey data on type, length, volume and other characteristics indicators of gully erosion were used to calibrate the results. Then the features of every subregion, such as where the gully erosion is, how serious it is, and why it happens and develops, were expounded. The result is not only an essential prerequisite for gully erosion surveys and monitoring, but also an important basis for gully erosion prevention.     

Local scour and the laws of scour pit’s shape downstream of debris flow sabo dam

The erosion shape and the law of development of debris flow sabo dam downstream is a weak part in the study on debris flow erosion. The shape and development of scour pit have an important effect on the stability and safety of debris flow sabo dam, which determines the foundational depth of the dam and the design of protective measures downstream. Study on the scouring law of sabo dam downstream can evaluate the erosion range and reasonably arrange auxiliary protective engineering. Therefore, a series of flume experiments are carried out including different debris flow characteristics （density is varying from 1.5 t/m3 to 2.1 t/m~） and different gully longitudinal slopes. The result shows that the scour pit appears as an oval shape in a plane and deep in the middle while superficial at the ends in the longitudinal section, the position of the maximum depth point moves towards downstream with an increase of flume slope angle. The maximum depth of scour pit is mainly affected by the longitudinal slope of gully, density of debris flow, and the characteristics of gully composition （particle size and the viscosity of soil）. The result also indicates that the viscosity of soil will weaken the erosion extent. The interior slopes of scour pit are different between the upstream and the downstream, and the downstream slope is smaller than the upper one. For the viscous and non-viscous sands with the same distribution of gradation, the interior slope of non- viscous sand is smaller than the viscous sand.According to tbe regression analysis on the experimental data, the quantitative relationship between the interior slope of scour pit, slope of repose under water and the longitudinal slope of gully is established and it can be used to calculate the interior slope of scour pit. The results can provide the basis for the parameter design of the debris flow control engineering foundation.     

地下采动顺向缓倾构造山坡的滑移机理

采空区地表山体侧向变动,不同于一般天然山坡,也与采空区一般上覆岩层的变形破坏有异;它是二者复合机理的效应。本文在分析考察了毗邻电厂的横山顺倾构造山体、剖析了地下采空情况认为:山体侧向变动中,软弱夹层有决定性作用,变动范围、速率和规模与地下采空有关;并利用地质力学模型试验和数值模拟,探索了采动引起山体应力场及变动规律。结果表明,山体岩层的变形、位移、破坏、由直接顶板向地表发展,采空坍陷诱发了软弱夹层的蠕滑,则产生山体侧向滑移;电厂区地表隆起变形是山体侧向滑移挤压地基土的反映。通过现场实际调研、变形观测资料分析与数值模拟和模型试验的对比研究,提出了“坍落拱梁”的成生效应、挤压蠕滑效应、失稳效应;揭露了顺倾构造山体在采空影响下,具有地表、地下的“复合临空面”的“复合应力场”中“复合变动”的“复合机理”;并提出这种山体侧向变动机理的典型地质模式,借以论证山体稳定性。     

An estimation model for the fragmentation properties of brittle rock block due to the impacts against an obstruction

Mountain hazards with large masses of rock blocks in motion – such as rock falls, avalanches and landslides – threaten human lives and structures. Dynamic fragmentation is a common phenomenon during the movement process of rock blocks in rock avalanche, due to the high velocity and impacts against obstructions. In view of the energy consumption theory for brittle rock fragmentation proposed by Bond, which relates energy to size reduction, a theoretical model is proposed to estimate the average fragment size for a moving rock block when it impacts against an obstruction. Then, different forms of motion are studied, with various drop heights and slope angles for the moving rock block. The calculated results reveal that the average fragment size decreases as the drop height increases, whether for free-fall or for a sliding or rolling rock block, and the decline in size is rapid for low heights and slow for increasing heights in the corresponding curves. Moreover, the average fragment size also decreases as the slope angle increases for a slidingrock block. In addition, a rolling rock block has a higher degree of fragmentation than a sliding rock block, even for the same slope angle and block volume. Finally, to compare with others' results, the approximate number of fragments is estimated for each calculated example, and the results show that the proposed model is applicable to a relatively isotropic moving rock block.     

ArcGIS环境下DEM的坡长计算与误差分析

坡长是水土保持、土壤侵蚀和环境评价等研究中的基本因子。在GIS环境下,坡长的提取是基于格网DEM进行的,其结果受到DEM精度、DEM结构、流向提取算法和栅格距离计算方式的影响。在确定分辨率和精度的DEM下,坡长计算误差主要来自流向提取算法和栅格距离计算方式。本文对这两种DEM坡长计算误差源进行了详细的分析,指出集成到ArcGIS软件中的坡长计算模型D8存在较大缺陷,尽管栅格距离计算方式对坡长影响最大误差为8.9%,并可通过统计方式实现整体坡长的修正,但D8算法的单流向特性使得水流方向变得确定和不连续,进而导致坡面单元坡长产生较大误差,因此,利用D8算法进行数值模拟计算需要考虑其精度问题。同时应提高坡长计算的准确性,需要发展更为完善的流向计算方法和坡长计算模型。     

岩质边坡平面滑动锚固方向角的三维优化方法

以岩质边坡中常见的平面滑动为研究对象,研究其最优锚固方向角的计算方法。将锚索自由段单位长度能提供的最大抗滑增量作为目标控制变量,将坡面和滑动面特征参数与锚索设计参数作为优化控制自变量,通过坐标系转换得到的线性方程组对锚索的预拉力进行分解,并根据锚索支护时的三维模型建立了锚索自由段长度的优化公式,进而推导了用于锚固方向角三维优化的新计算方程。在该方程的基础上借助MATLAB软件中的fmincon函数对锚索加固方向不受限制时的锚固方向角进行了优化。最后通过锚固方向角敏感性分析与工程实例分析相结合,证明了推荐的最优锚固方向角计算方法的有效性与先进性。新方法有效地解决了边坡坡面与滑动面走向存在夹角时最优锚固方向角的求解问题,可进一步提高锚索的锚固效益,降低边坡的支护费用。      

Deformation features and failure mechanism of steep rock slope under the mining activities and rainfall

Underground mining activities and rainfall have potential important influence on the initiation and reactivation of the slope deformations, especially on the steep rock slope. In this paper, using the discrete element method (UDEC), numerical simulation was carried out to investigate deformation features and the failure mechanism of the steep rock slope under mining activities and rainfall. A steep rock slope numerical model was created based on a case study at the Wulong area in Chongqing city, China. Mechanical parameters of the rock mass have been determined by situ measurements and laboratory measurements. A preliminary site monitoring system has been realized, aiming at getting structure movements and stresses of unstable rock masses at the most significant discontinuities. According to the numerical model calibrated based on the monitoring data, four types of operation conditions are designed to reveal the effect of mining excavation and extreme rainfall on the deformation of the steep rock slope.     

Wind tunnel study of multiple factors affecting wind erosion from cropland in agro-pastoral area of Inner Mongolia, China

In this paper,the process of wind erosion on two kinds of soil from the agro-pastoral area of Inner Mongolia are studied using wind tunnel experiments,considering the wind speed,blown angle of wind and soil moisture content.The results showed that the modulus of soil wind erosion increases with an increase of wind speed.When the wind speed exceeds a critical value,the soil wind erosion suddenly increases.The critical speed for both kinds of soil is within the range of 7-8m·s-1.For a constant wind speed,the rate of soil wind erosion changes from increasing to falling at a critical soil slope.The critical slope of loam soil and sandy loam soil is 20° and 10°,respectively.Soil moisture content has a significant effect on wind erosion.Soil wind erosion of both soils decreases with an increase of the soil water content in two treatments,however,for treatment two,the increasing trends of wind erosion for two soils with the falling of soil water content are no significant,especially for the loam soil,and in the same soil water content,the wind erosion of two soils in treatment one is significantly higher than treatment two,this indicates reducing the disturbance of soil surface can evidently control the soil wind erosion.     

泥质砂岩残积土边坡降雨冲刷特性

泥质砂岩残积土作为一种结构性很强的特殊土,具有崩解性强、抗冲蚀性差以及扰动性极大的特点,对工程建设有较大影响。为了探究泥质砂岩残积土边坡降雨冲刷机理,设计了边坡降雨冲刷试验,通过现场三维激光扫描技术测试分析了其表面冲刷效应;利用高密度电法进一步明确了泥质砂岩残积土边坡的入渗特性、表面冲刷演化机制及冲刷破坏机理。结果表明：冲刷试验的最初阶段,降水入渗强且主要向坡脚处运移,坡表未形成明显的细沟;冲刷试验中期,坡脚处土体最先达到饱和而形成坡面径流,细沟贯通扩大形成小规模冲槽以及片蚀区;冲刷试验后期,坡面中部和坡脚处土体冲蚀严重,坡脚处的冲槽向上部延伸,片蚀区扩大,导致表层土体结构发生变化,渗透性差异明显;泥质砂岩残积土坡体降雨冲刷主要划分为表层溅蚀、下层潜蚀和细沟贯通3个阶段,坡面土体流失主要发生在最后一个阶段,细沟率达到最高值16.9%,细沟贯通率也高达0.74。研究结果可以为深入探讨泥质砂岩残积土边坡冲蚀防护和研究冲蚀防护机理提供基础资料。     

公路膨胀土路堑边坡防护方案探讨

针对膨胀土路堑边坡的公路工程特点，总结以往膨胀土路堑边坡的防护措施，根据以往膨胀土路堑边坡的防护成败资料和经验。以保持边坡土体天然含水量状态的相对稳定，保持边坡土体结构的相对完整性。保持边坡土体足够的抗剪强度和防护工程应能适应边坡膨胀土体可能产生的膨胀变形与膨胀力不遭破坏。首先用DAH混合溶液改良膨胀土边坡的表层土，使DAH混合溶液与膨胀土产生离子交换发生化学反应，从而使表层膨胀土体改变性能，成为正常土，然后结合以往路堑边坡处理方法进行综合处理。实践证明，此方法不但充分利用了膨胀土．节约了工程投资，而且又绿化、美观了边坡，使其与周围环境相协调，从而为膨胀土地区公路路堑边坡处理提供了一种切实可行的防护措施。     

Mechanism of downcutting erosion of debris flow over a movable bed

The phenomenon of debris flow is intermediate between mass movement and solid transport. Flows can be sudden, severe and destructive. Understanding debris flow erosion processes is the key to providing geomorphic explanations, but progress has been limited because the physical-mechanical properties, movement laws and erosion characteristics are different from those of sediment-laden flow. Using infinite slope theory, this research examines the process and mechanism of downcutting erosion over a moveable bed in a viscous debris flow gully. It focuses specifically on the scour depth and the critical slope for viscous debris flow,and formulas for both calculations are presented.Both scour depth and the critical conditions of downcutting erosion are related to debris flow properties(sand volume concentration and flow depth) and gully properties(longitudinal slope,viscous and internal friction angle of gully materials,and coefficient of kinetic friction). In addition, a series of flume experiments was carried out to characterize the scouring process of debris flows with different properties. The calculated values agreed well with the experimental data. These theoretical formulas are reasonable, and using infinite slope theory to analyze down cutting erosion from viscous debris flow is feasible.     

黄土高原丘陵沟壑区切沟侵蚀与地形关系分析——以纸坊沟流域为例

沟蚀是土壤侵蚀研究的主要内容之一,地形地貌是沟蚀的一个重要影响因子。本文以安塞纸坊沟流域作为研究区域,选择土地利用方式、土壤类型、坡度坡长因子、平面曲率、坡向和地形湿度指数6个因子,通过因子内切沟所占比重/整个研究区切沟所占比重计算各个因子的权重值,通过空间叠加分析土壤侵蚀敏感性,并通过重分类的方法把土壤侵蚀敏感性分为基本无侵蚀、轻度侵蚀、中度侵蚀、强度侵蚀、剧烈侵蚀5个等级,来研究切沟侵蚀与地形的关系。结果表明：切沟多发生在坡度坡长较大、地表湿度较高的林草地区域和更容易发生侵蚀的黄绵土区域,并且多分布在阴坡的凹面;对比分析切沟侵蚀和土壤侵蚀敏感性,切沟大多分布在中等侵蚀敏感性以上的区域,约占总切沟的90%;实验权重值对验证区冲沟的响应精度为82.43%（中度侵蚀及其以后阶段）,与实际值90.53%相差不大,说明此种方法对黄土丘陵沟壑区具有一定适用性,对黄土丘陵沟壑区水土保持工作有重要意义。     

济南保泉供水数据库的建立及其功能简介

济南保泉供水数据库的建成,为及时预报地下水动态及优化调整地下水开采方案提供了一个信息管理系统。文中介绍了该数据库的建库情况、开发功能模块及在保泉供水工作中的应用。     

Deep profile soil moisture distribution characteristics for different microtopographies on the semi-arid Loess Plateau,China

Soil moisture is a limiting factor for vegetation restoration on the Loess Plateau, China. Micro-topography may cause heterogeneities in the distribution of soil moisture, but little is known about its effect on deep soil moisture. Our study aims to explore the distribution and impact of soil moisture within the upper 10 m of soil for different microtopographies. Taking undisturbed slope as the control, five micro-topographies were selected. Soil moisture over a depth of 0-10 m from 2017 to 2019 was investigated, and soil particle size and soil organic matter were measured. Variance analysis and multiple comparisons were used to analyze the difference in soil moisture for different microtopographies and multiple-linear regression was used to analyze the influence of micro-topography on soil moisture. There are significant differences in soil moisture within the different layers underlying the examined micro-topographies, while the inter-annual variation in soil water storage for the selected microtopographies increase with increased rainfall. The depth of influence of micro-topographic vegetation on soil moisture exceeded 1000 cm for a gully(GU), 740 cm for a sink hole(SH), 480 cm for a scarp(SC), 360 cm for an ephemeral gully(EG) and 220 cm for a platform(PL). Micro-topography will cause the heterogeneous distribution of soil moisture in the shallower layers, which changes the vegetation distribution differences between micro-topographies. This may be the survival strategy of herbaceous vegetation in response to climate change in the Loess Plateau. For future vegetation restoration efforts, we need to pay attention to the influence of microtopography on soil moisture.