高速远程滑坡裹气流态化模型试验研究

为分析高速远程滑坡运动过程中碎屑流的裹气流态化特性,以汶川地震触发的谢家店子滑坡为原型,运用自行设计的流化床试验装置,开展了一系列高速远程滑坡碎屑流裹气流态化特性的试验研究。结果表明：(1) 对于颗粒粒径变化范围广(0.1～7 mm者均有)、颗粒形状极不规则、且底部细颗粒含量高的类滑坡碎屑,其床层压降随着供风量的增加并非呈线性增长趋势,而是表现出明显的气泡控制性特征,随堆积体内部气泡和气垫层位置的不同,堆积体各部位床层压降相差较大;(2) 当供风量较低时,碎屑堆积体整体稳定性较好,在堆积体表面仅可见有限气体逸散孔的分布;当供风量达到一定值后,堆积体内部流态化特性增强,可见砂质密云、气孔、纵向脊等一系列与高速远程滑坡实例中相类似现象的出现;(3) 对于滑坡类不规则颗粒,其间容易形成架桥现象,增强碎屑流整体稳定性,使气流不易扩散,在堆积体下方形成高压气垫层,产生空气擎托现象。     

颗粒材料底面动摩擦系数特征研究

离散态颗粒物质具有明显不同于普通固体的界面摩擦特性,而摩擦系数是界面摩擦特性的主要表征参数之一。通过倾斜仪开展不同级配条件下颗粒材料的滑动摩擦试验,基于视频图像解析以及函数拟合方法,建立滑动位移与滑动时间的最佳函数拟合关系,分析滑动过程的加速度并推算底面动摩擦系数,研究颗粒粒径、质量配比等级配因素对颗粒材料底面动摩擦系数的影响。研究结果表明:(1)各级配颗粒材料的平均底面动摩擦系数随着运动时间的增加均呈线性减小趋势;(2)对于单粒径材料,与粗颗粒相比,细颗粒具有较大的底面动摩擦系数;(3)对于双粒径材料,随着细颗粒含量的增加,颗粒材料的平均底面动摩擦系数先急剧降低至最小值(细颗粒含量≤40%),后急剧增加(细颗粒含量40% ~60%),最终增加趋势明显变缓(细颗粒含量≥60%)。     

高速远程滑坡颗粒流研究进展

颗粒流是高速远程滑坡物质演化过程的一个重要阶段,也是从细观角度揭示其超常运动特性的重要手段.颗粒流所采用的主要研究方法以及取得的重要理论成果,可以为高速远程滑坡动力学机理的研究提供重要的技术手段和理论支持.本文聚焦颗粒流研究进展,从颗粒流基本概念、流态特征及流变本构模型、颗粒流粒径分选机制等方面进行了系统梳理;进而,从理论、实验及数值计算模型3个方面对高速远程滑坡研究中颗粒流理论及方法的应用进行了系统性述评;在此基础上,提出了从颗粒流角度研究高速远程滑坡动力学机理涉及的关键科学问题:高速远程滑坡高流动性的起源涉及哪些物理过程？如何量化和模拟其多分散性和破碎过程？如何量化描述颗粒尺寸分布的时空演变及其与流动的耦合？如何从其沉积特征中探究流动的传播机制？针对这些问题,从基于沉积学特征的颗粒流物理力学过程、考虑尺度效应的颗粒流动力学特性研究、基于颗粒流力学过程的滑坡运动机理及其本构模型、新技术新方法的应用4个方面提出下一步应重点开展的研究工作.      

重力流的运动特点和类型

作者以鄂尔多斯地区奥陶系丰富的重力流沉积为例,借助非牛顿流体力学的基本原理,对重力流运动的特点和类型作了进一步的探讨,重力流运动的特点是:具层移(流)运动、有旋性、挤压性、涡旋性、惯性、颗粒上浮性、颗粒下沉、颗粒碰撞、剪切作用和下切作用。重力流可分为:滑动流、混合流、碎石(屑)流、颗粒流、液化流和浊流。     

滑坡-碎屑流物理模型试验及运动机制探讨

滑坡-碎屑流由于高速、远程的特点常常引发灾难性事故,其复杂的运动机制导致预测致灾范围非常困难。通过开展室内模型试验,研究了碎屑粒径、滑床糙率和挑坎对运动特性的影响。试验结果发现,滑坡碎屑运动距离受控于前端碎屑,且随着碎屑的粒径增大而增加,增加滑床糙率、挑坎均可使碎屑的运动距离减小。在前人研究成果的基础上结合碎屑材料的力学特性探讨了滑坡-碎屑流出现流态化的原因和高速远程机制,即高速运动中颗粒间的作用力远小于完整岩体,因此颗粒间的“黏聚力”不能维持滑坡体的整体性,同时致使滑坡体与滑床接触的过程中传递至滑坡体内部的摩阻力减少,从而导致碎屑滑坡的远程结果。     

颗粒粒径对滑坡碎屑流动力特征及能量转化的影响——以四川省三溪村滑坡为例

滑坡碎屑流的能量演化机制涉及复杂的碰撞、摩擦和能量转化,对滑坡灾害的防治具有重要意义。以四川省三溪村滑坡为例,采用离散元法建模,研究了不同土颗粒粒径下滑坡的运动堆积特征、能量演化过程及其对建筑物的影响。结果表明,在不同粒径条件下,滑坡的堆积特征变化不大,但它们对能量转换和建筑物的冲击力有显著的影响。粒径越大,滑坡启动速度越快,峰值动能速度越高,碰撞耗散能量越大,摩擦耗散能量越少。粗粒土中颗粒间距越大,颗粒间的碰撞效应越明显,有利于能量传递,因此,对房屋的冲击力越大。因此,在模拟过程中不能忽视颗粒粒径对滑坡碎屑流动力学特征的影响。这些研究结果揭示了不同粒径土粒在滑坡运动过程中的能量演化机制,为能量演化对建筑冲击的影响提供了初步的认识,可为滑坡碎屑流的防治提供指导依据。     

高速远程滑坡运动堆积过程中的能量传递机制

高速远程滑坡往往引发灾难性事故,开展运动堆积过程定量研究,对于探究滑坡发生机理及预测致灾范围具有重要意义.基于室内物理模型试验,通过PIV技术分析高速摄像机在试验过程中拍摄的照片,获取了运动过程中滑体颗粒的水平速度、竖直速度与位移等运动参数.从滑体颗粒群和不同位置处单体颗粒角度,分析高速远程滑坡的运动演化规律.结果显示:（1）滑体颗粒前端出现高速区,该高速区随着滑坡停止具有一定的保持性.颗粒间存在明显的碰撞现象;（2）从不同位置单体颗粒来看,前部颗粒位移量最大,速度波动频繁,碰撞频次最高,能量多次补充;中部颗粒位移量其次,速度有波动过程,但不及前部频繁;后部颗粒位移量最小,速度基本呈递降趋势,能量逐渐减小.结合重庆鸡尾山滑坡以及Black Rapids Glacier滑坡实例分析,揭示了高速远程滑坡运动堆积过程中滑体颗粒间存在碰撞及能量传递现象,从而进一步探究高速远程滑坡形成机制,在监测预防、灾害治理等方面具有现实指导意义.      

论崩塌滑坡—碎屑流高速远程问题

高速远程崩塌滑坡—碎屑流具有规模大、速度快、滑程远、多态化、常转向、冲程多、冲击性和摧毁性等特征。崩塌或滑坡高速远程实质上是其解体后的碎屑流运动形式。碎屑流高速远程与崩塌滑坡规模、物质成分结构、地形高差、沟道形态和引发因素及运动路经的环境等因素密切相关。崩塌滑坡变形破坏形式一般显示为蠕动—拉裂—剪断—滑移—冲出—解体—碎屑流化的过程。崩塌滑坡形成机理主要基于残余强度、蠕变作用和孔隙水压力等理论认识进行解释。碎屑流运动机理主要立足于势能动能转化传递、气体浮托和颗粒流运动理论予以解释。动势能转化、气体浮托作用和颗粒流运动三种解释是层次不同、相互补充的关系,不是彼此独立的。基于成年人在复杂地形下能够奔跑逃生,崩塌滑坡—碎屑流前锋的运动速度5m/s作为高速运动的下限值是比较合理的。崩塌滑坡—碎屑流区域的前后缘高差(H)与前后缘水平距离(L)的比值小于0.4或L/H值大于2.5可作为其远程运动的判据。崩塌滑坡—碎屑流成灾模式包括直接压覆、解体推挤、碰撞冲击、气浪吹袭、激流涌浪、堰塞湖淹没与滑坡坝溃决—洪水泥石流等多种形式。     

基于PFC3D的鱼鳅坡滑坡运动过程分析

以贵州省开阳县鱼鳅坡滑坡为研究对象,采用颗粒流离散元（PFC3D）对其破坏运动过程进行数值模拟。采用Ball-Wall建模方法建立滑坡模型,对滑坡不同关键部位颗粒进行位移、速度监测,阐明其破坏运动特征。结果表明,降雨为鱼鳅坡滑坡的直接诱发因素。该滑坡在破坏初始阶段以蠕滑变形为主,随着变形量的增加,滑坡体不断挤压坡脚,滑坡岩土体到达应力平衡极限,坡脚产生剪切破坏,并向上牵引发展,滑坡发生整体滑动,斜坡变形破坏模式为蠕滑-拉裂,按照力学条件为牵引式破坏。滑坡滑动最高时速12.4 m/s,最大滑移80 m,滑动阶段持续50 s。研究成果可为对该类滑坡影响范围预测,以及工程措施的制定具有一定的参考意义。     

基于PFC2D的刘涧滑坡破坏运动过程模拟

以陕西省洛南县刘涧滑坡为研究对象,采用颗粒流离散元法对其破坏运动过程进行数值模拟。首先通过双轴数值试验对滑坡饱和土体进行细观参数标定,并与室内试验中饱和土体宏观力学参数进行对比,经验证该细观参数能应用到滑坡的破坏运动分析中,进而引入颗粒流（PFC2D）程序中平行黏结模型,采用ball-wall建模方法建立滑坡模型,对滑坡不同关键部位颗粒进行位移、速度监测,阐明其破坏运动特征。模拟结果表明,降雨为刘涧滑坡的直接诱发因素,斜坡变形破坏模式为由坡脚开挖引起的自前缘向后部牵引-孔隙水压力诱发的后部向前缘推移式滑塌。总体特征为上部推移,中部剪切,下部牵引;滑坡滑动最高时速13.4 m/s,最大滑移170 m,滑动阶段持续25 s。利用颗粒流法对滑坡的破坏运动过程模拟具有较好的适用性,可为工程决策提供依据。      

Experimental study of the dynamic behavior and segregation of density-bidisperse granular sliding masses at the laboratory scale

From the understanding of dynamics and processes of rapid granular flows and the granular-segregation mechanism in gravity-driven flow, we can clarify the particle-composition structure in the downstream areas of avalanches in geophysical contexts, such as landslides, rock falls, and snow-slab avalanches. Such dynamics also provide a basis for geophysical studies. This study experimentally investigates the dynamic behavior and segregation phenomena of a density-bidisperse, rapid, granular flow down a quasi-2D, rough, inclined rectangular chute. Particles with two density ratios are used to investigate the mechanism of density-induced segregation, and four chute-inclination angles are tested to examine the influence of driving forces. The dynamics of the mixture flow—which includes the flow-depth evolution, stream-wise and depth-wise velocity profiles, shear rate, and granular temperature in the upper high-shear band of the flow—are obtained from particle image velocimetry (PIV) measurements. The two-dimensional concentration distributions of the particles in the stream-wise direction are also obtained using 2D image processing to determine the segregation state. In the upstream region, the variation in the concentration of heavier particles is defined as the strength of the density-induced segregation state, S_d. Our results indicate that the mixture-flow parameter—particularly the shear rate and the granular temperature in the upper high-shear band—crucially influence the strength of particle segregation in granular avalanches. In the upstream region, a higher shear rate and a higher granular temperature in the upper high-velocity band result in a smaller drag force in the mixture flow, causing stronger density-induced particle segregation. These results well describe the entire processes of dense granular flows, from upstream initiation to the downstream steady state. Therefore, they reveal the structure of the mixed flow in the depth direction and are expected to explain various gravity-driven mixture granular flows.

     

From solid to granular gases: the steady state for granular materials

This paper aims at extending the well‐known critical state concept, associated with quasi‐static conditions, by accounting for the role played by the strain rate when focusing on the steady, simple shear flow of a dry assembly of identical, inelastic, soft spheres. An additional state variable for the system, the granular temperature, is accounted for. The granular temperature is related to the particle velocity fluctuations and measures the agitation of the system. This state variable, as is in the context of kinetic theories of granular gases, is assumed to govern the response of the material at large strain rates and low concentrations. The stresses of the system are associated with enduring, frictional contacts among particles involved in force chains and nearly instantaneous collisions. When the first mechanism prevails, the material behaves like a solid, and constitutive models of soil mechanics hold, whereas when inelastic collisions dominate, the material flows like a granular gas, and kinetic theories apply. Considering a pressure‐imposed flow, at large values of the normal stress and small values of the shear rate, the theory predicts a nonmonotonic shear rate dependence of the stress ratio at the steady state, which is likely to govern the evolution of landslides. Copyright © 2013 John Wiley & Sons, Ltd.     

岩石碎屑流运移堆积过程数值模拟

岩石高速远程崩滑是一种特殊的、危害性很大的地质灾害,崩滑后期形成的高速岩石碎屑流具有很高的运移能力,破坏力极强。以岩石碎屑流为研究对象,重点模拟岩石碎屑流形成后,碎屑流的运移堆积过程和最大运移堆积距离。利用二维颗粒离散元数值模拟方法,建立碎屑流二维运动堆积模型,分析了岩石碎屑流初始高度、体积、加速斜面坡度、堆积底面摩擦系数以及堆积底面地表起伏程度对于岩石碎屑流最大运移堆积距离的影响。实验结果显示:岩石碎屑流运移堆积过程体现出了碎屑流物质的离散性和流动性; 碎屑流最大堆积距离随高度、体积、斜面坡度增大而增大; 随堆积底面摩擦系数、堆积底面地表起伏程度增大而减小。     

青藏高原古高速远程滑坡沉积学特征研究

高速远程滑坡运动学机理是国际工程地质领域亟待解决的重大前沿性关键科学问题.为探索高速远程滑坡的运动学机理,以青藏高原不同构造背景下的三大高速远程滑坡为研究对象,通过现场工程地质调查与分析,详细探讨了滑坡运动路径上所揭露出的各种表面与剖面沉积学特征,结果显示:(1)沿滑坡运动路径上依次可见大型堆积平台、纵向脊、横向脊、堆...     

岩质滑坡−碎屑流的运动距离计算公式研究

滑坡?碎屑流的远程运动距离是碎屑流体所能够达到的最大堆积距离,是灾害预警和评估的重要指标。通过总结已有碎屑流运动距离研究成果,从岩体势能入手开展研究,并结合量纲分析,首先建立了运动距离与势能之间的基本方程。其后,采用4种颗粒材料开展岩质碎屑流滑槽试验,研究碎屑体体积V、滑移区坡度?、碎屑粒径d以及最大垂直运动距离H等对碎屑流运动距离L的影响,通过逐步拟合回归,建立了基于势能的岩质滑坡?碎屑流最大水平运动距离的计算公式。最后,采用汶川地震触发的38个岩质碎屑流以及17个其他典型岩质滑坡碎屑流数据对计算公式进行了验证,结果表明,考虑该运动距离计算公式具有较好的可靠性,能够为山区滑坡?碎屑流灾害预警工作提供理论指导。     

A numerical study of particle motion and two‐phase interaction in aeolian sand transport using a coupled large eddy simulation – discrete element method

Wind‐blown sand movement, considered as a particle‐laden two‐phase flow, was simulated by a new numerical code developed in the present study. The discrete element method was employed to model the contact force between sand particles. Large eddy simulation was used to solve the turbulent atmospheric boundary layer. Motions of sand particles were traced in the Lagrangian frame. Within the near‐surface region of the atmospheric boundary layer, interparticle collisions will significantly alter the velocity of sand. The sand phase is quite dense in this region, and its feedback force on fluid motion cannot be ignored. By considering the interparticle collision and two‐phase interaction, four‐way coupling was achieved in the numerical code. Profiles of sand velocity from the simulations were in good agreement with experimental measurements. The mass flux shows an exponential decay and is comparable to reported experimental and field measurements. The turbulence intensities and shear stress of sand particles were estimated from particle root‐mean‐square velocities. Distributions of slip velocity and feedback force were analysed to reveal the interactions between sand particles and the continuous fluid phase.     

不同粒度颗粒在非恒定管流中跟随性试验

管道中固液两相流水击对管道和输送系统可能产生严重的破坏,而固液两相在这种非恒定流中的运动特性是计算最大水击压力变化的重要依据。采用粒子图像测速技术（Particle Image Velocimetry,PIV）,通过试验研究水击发生时,水平圆管中不同平均流速、颗粒粒径条件下,流体介质和粗颗粒在管道断面的速度分布以及粗颗粒跟随性的变化规律。研究结果表明:①在水击发生的不同时刻,圆管流中粗颗粒的流速在管道断面分布呈不规则的抛物线型分布,主要表现为靠近管道壁面底部的颗粒流速略小于靠近管道顶部流速,当颗粒粒径大于1.5 mm,平均流速小于2.5 m/s时,粗颗粒表现出明显的沉降特性;②粗颗粒的跟随性与颗粒受力有密切关系,其中颗粒速度与流体速度的变化量是影响颗粒受力的重要参数;③基于试验数据拟合得到了水击条件下粗颗粒跟随性系数k的经验公式,并分析了颗粒粒径、管道直径、两相流平均流速以及水击发生时间等不同参数对粗颗粒跟随性系数的影响,公式计算值与实测值之间的误差在5%以内。     

Process dependence of grain size distributions in rock avalanche deposits

Rock avalanches are a form of hazardous long-runout landslide and leave fragmented deposits of complex sedimentology that, if studied in detail, can provide insight into their emplacement processes. Complexity arises due to the myriad overlapping factors known to contribute to the final deposit fabric, such as source structures, lithology (i.e. material properties), topographic feedback, substrate interaction and emplacement processes (i.e. internal factors), as well as our reliance on (un)suitable exposures. Herein, we present sedimentological data from two carbonate rock avalanche deposits (Tschirgant in Austria and Flims in Switzerland), where changes in lithology can be eliminated from the causal equation due to their largely mono-mineralic composition. We further eliminated the effects of external influences such as topography or substrate interactions by detailed facies mapping of the deposit interior. Since sedimentary properties locally vary within less than 1-m² outcrop area, emplacement processes are the only causes that remain to explain the different fabrics. Characteristic (fractal) grain size distributions of three distinctive sub-facies in the interior of these, and other, rock avalanche deposits—jigsaw-fractured, fragmented, and shear zone facies—can be linked to specific processes acting during emplacement. We suggest that a heterogeneously distributed and progressively increasing particle breakage in the moving granular mass best explains the ranges of fractal dimensions and associated features for the respective sub-facies, from simple breakage along pre-existing planes, through dynamic fragmentation which locally minimises coordination number, to zones of shear concentration. No exotic emplacement mechanisms (such as air-layer lubrication or fluidised substrates) are required to produce these features; continued, heterogeneous degrees of fragmentation of an initially intact source rock best explains the sedimentary record of rock avalanches.     

微粒对多晶冰流变行为的影响(Ⅰ):应力下的流变行为

研究了应力下微粒对多晶冰流变行为的影响.结果表明:微粒分布在晶界或同时分布在晶界和晶内都将增加多晶冰的流变速率.当微粒只分布在晶界时,含微粒浓度为1 wt.%的多晶冰有最高的最小流变速率;当微粒同时分布在晶内与晶界时,多晶冰的最小流变速率基本与微粒浓度无关.微粒提高流变速率的原因来源于提高了流变过程中的位错密度,分布在晶界和晶内的微粒可以分别通过阻碍晶界滑移,发展内应力及Frank-Read位错增值机制提高流变过程中的位错密度.