基于熵权-正态云模型的泥石流灾害易发性评价

近年来,由于极端天气事件频发和工程建设活动加剧,泥石流灾害对人类的威胁愈加严峻。因此,在综合分析泥石流灾害发育特征及形成机理的基础上,进行泥石流灾害易发性评价对防灾减灾工程活动具有指导意义。选取沟谷岸坡坡度、沟床纵坡比降、植被覆盖率、单位面积固体物源储量、汇水面积、雨季降雨量建立评价指标体系,采用熵权法确定评价指标权重,将泥石流灾害易发性分为高易发、中易发、低易发和不易发四级,建立基于正态云模型的泥石流灾害易发性评价方法,并以西秦岭地区5条泥石流沟为例验证了该评价模型的合理性。进一步将该模型用于陕西吴堡井沟泥石流灾害易发性评价,结果表明井沟泥石流具有中易发性,评价结果与实际相符。     

基于Fisher判别法的台风暴雨泥石流预测模型

对台风暴雨泥石流发生的可能性进行定量预测,有助于减少危险区内的人员伤亡、降低经济损失。以台湾地区南投县陈有兰溪流域的47条泥石流沟为研究对象,从泥石流形成所需的地形地貌、物源和降雨条件中,初步选取台风暴雨泥石流发生的影响因子,包括沟床平均坡度、有效流域面积、形状系数、主沟长度、岩性、崩滑比和平均雨强。根据因子重要性排序结果,选择崩滑比和平均雨强作为模型的预测因子,基于Fisher判别法建立了台风暴雨泥石流预测模型。采用随机取样技术,选取70%的数据用于构建模型,剩余30%的数据用于验证模型。以精确度、准确率、漏报率和误报率指标,定量评价模型的预测效果,并确定最优的预测模型。结果表明:基于Fisher判别法构建的台风暴雨泥石流预测模型,综合考虑了泥石流形成所需的物源条件和降雨条件,弥补了降雨阈值模型仅依靠降雨资料分析的不足,预测效果更好。     

Size of a debris flow deposition: model experiment approach

 The prediction of the dangerous extent of a debris flow deposition is of vital importance, but difficult to achieve. Precise prediction of the depositional boundary of a debris flow event is impossible, but the size of a debris flow deposition could provide some estimates of the area, length, width, and thickness of a debris flow deposition. Based on in situ depositional experiments performed on a debris flow creek just after debris flows, a rule of thumb expressed by a group of equations containing the multiple-variate nonlinear functions is proposed in this paper. The interrelationships between the size and the causation also are discussed, and some empirical formulae to calculate the causative parameters for different regions are presented. Received: 24 April 1995 · Accepted: 21 June 1995     

北京市泥石流易发区降雨预警阈值研究

泥石流灾害合理的雨量预警阈值不仅与历史泥石流灾害发生时的降雨量有关,且与研究区域的气候、地形地貌、地质、植被等密切相关。论文采用雨场分割法和GIS技术研究了影响泥石流启动的降雨和地质背景两大因素,在对北京市泥石流灾害易发分区的基础上,结合北京地区已发生的82起泥石流的易发性分区和雨量值,提出了不同泥石流易发等级条件下的雨量预警阈值。研究成果已经在2015年7月17日北京房山区西区沟泥石流预警中成功应用,为泥石流区域预警预报提供了一种新的思路。     

Regional debris flow susceptibility analysis based on principal component analysis and self-organizing map: a case study in Southwest China

A method was developed to analyze the susceptibilities of 541 regional basins affected by debris flows at the Wudongde Dam site in southwest China. Determining susceptibility requires information on source material quantity and occurrence frequency. However, the large number of debris flows can hinder the individual field investigation in a each small basin. Factors that may trigger debris flows can be identified using remotely sensed interpretation information. Susceptibility analysis can then be conducted based on these factors. In this study, SPOT5 satellite imagery, digital elevation models (DEM), a lithology distribution map, and rainfall monitoring data were used to identify 12 debris flow trigger factors: basin relief ratio, slope gradient in the initiation zone, drainage density, downslope curvature of the main channel, vegetation coverage, main channel aspect, topographic wetness index, Melton’s ruggedness number, lithology, annual rainfall, form factor, and cross-slope curvature of the transportation zone. Principal component analysis was used to obtain the eight principal components of these factors that contribute to susceptibility results. Then, a self-organizing map method was adopted to analyze the principal components, which resulted in a debris flow susceptibility classification. Field validation of 26 debris flow basins was used to evaluate the errors of the susceptibility classification, as well as assess the causes of such errors. The study found that principle component analysis and self-organizing map methodologies are good predictors of basin susceptibility to debris flows.     

基于黏弹塑性本构模型的泥石流数值模拟

泥石流是一个世界范围的地质灾害问题,多年来受到广泛而深入的研究。大多数工作致力于分离地研究岩土材料破坏前的失稳机制和失稳后的快速流动扩展。尝试以一种连续方法模拟完整的泥石流过程,仅用一个本构模型描述泥石流的起始、扩展和停滞。首先讨论了重力增大和基质吸力降低诱发泥石流的力学机制,描述了泥石流过程中岩土材料力学性质的弹塑性–黏性转变现象,对泥石流进行了黏弹塑性的全程模拟,并对挡土墙的冲击效应进行了启发性研究,证明了研究中所使用的本构模型可以准确地描述岩土材料的固流转化现象,数值方法可以很好地对均质连续介质进行大位移模拟。     

兰州市区泥石流趋势分析的系统熵模型

本文依据泥石流灾害体形态确定的地貌发育程度构建了泥石流灾害的评价模型,并将熵理论综合应用于泥石流评价模型,应用该模型对兰州地区山谷泥石流趋势进行了分析,并指出了评价模型的不足.     

运用Voellmy模型方法模拟泥石流运动堆积动力过程

泥石流流体具有速度快、破坏力强的特征,严重威胁着山区城镇的建设、民众生命以及财产安全。如何针对泥石流的运动堆积过程进行模拟一直是泥石流研究的热点和难点。基于近似Voellmy解的连续介质理论,建立泥石流动力模型,采用数值模拟仿真的方法,以一典型的泥石流沟为研究实例进行分析,通过模拟泥石流运动、堆积的动力过程,分析其最大速度、压力、过流量、堆积高度以及堆积面积,为后期泥石流危险性评价区划及防治提供有力支撑。     

基于HEC-HMS模型的不同雨型泥石流流量变化特征

为了探究不同雨型条件下泥石流流量变化特征,通过搜集汶川震区典型泥石流降雨数据,将其概化为三次峰值早到型、三次峰值型、三次峰值晚到型三种雨型。以HEC-HMS水文模型构建高家沟流域模型,在获得流域清水流量结果基础上采用雨洪修正法计算不同雨型下泥石流流量。结果表明三次峰值早到型、三次峰值型、三次峰值晚到型峰值清水流量分别为33.5,41.5,45.8 m3/s,泥石流峰值流量分别为166.83,206.67,228.08 m3/s,误差为-25.6%、-7.8%、1.7%;三种雨型下流量从上游至沟口以线性方式演变,且随降雨峰值推迟其演变速率不断增大,而增长幅度逐渐减小;泥石流暴发出现在峰值降雨前后,属于降雨激发型泥石流;三种雨型随降雨峰值推迟,流量增长阶段所需时间越长,分别为5.5,6,9 h,而衰退阶段所需时间越短,分别为14.5,8,2.5 h。研究表明该模型能为缺少降雨监测数据地区泥石流研究提供技术支持。     

基于模型试验的泥石流坡面物源启动预警模型

为进一步探明银洞子沟物源区坡面物源启动机理与降雨的相关性,在进行详细的现场考察与工程勘察工作后选取了典型模型体,并采用20余组人工降雨物理试验方法将降雨强度与坡度设为控制变量,研究了4种坡形、5种雨强条件下,坡面松散物源的失稳机制及破坏模式。试验结果揭示了地表变形与地下物理力学参数变化的定量响应关系,并基于试验发生破坏的临界雨量,建立了银洞子沟传统I-D预警模型。之后提炼可靠预警参数(坡度、深部体积含水率),通过数学回归分析方法并采用Exponential模型,得出IGD、IGM多参数新型预警判别式,实现了传统I-D模型的有效修正,并具一定的可靠性和实用价值。     

Exploring the impact of introducing a physical model into statistical methods on the evaluation of regional scale debris flow susceptibility

Natural Hazards - Regional scale debris flow susceptibility is widely evaluated by statistical methods. However, the initiation mechanism of debris flow is not considered, which leads to the...     

泥石流的二维数学模型

泥石流是在重力作用下,由砂粒石块和水等组成的固液混合物,是一种发生于山区的复杂的地质灾害现象。泥石流主要是由暴雨诱发引起的,它沿着复杂的三维地形高速流动,具有流体流动的特性。为了模拟泥石流的运动规律,预测降雨诱发的泥石流的到达距离和泛滥范围,减少和避免泥石流引起的灾害,把泥石和雨水组成的固液混合物假定为遵循均匀、连续、不可压缩的、非定常的牛顿流体运动规律。基于质量守恒方程和Naiver-stokes方程,采用深度积分方法,推导出了一个模拟泥石流运动的二维数学模型。所有方程式可用有限差分法来求解。结合GIS,该模型可用于预测泥石流的流动距离和泛滥范围,以及泛滥范围内的危险房屋和路段,也可以用于泥石流灾害的风险性分析。     

泥石流汇流槽可靠度分析

汇流槽是泥石流治理工程中常用的工程措施之一。文章基于汇流槽的倾倒破坏、滑移破坏及地基破坏模式，通过对汇流槽设计影响强烈的岩土参数、几何特征、荷载因素等敏感因子的求解方法及分布特征分析，从汇流槽的抗倾稳定、抗滑稳定和基底应力要求三方面建立极限状态方程。对于每个极限状态方程，在泰勒级数对非线性极限状态方程线性化后，运用一次二阶矩的理论，先假定失效点P^＊，求解出相应的均值和方差。然后根据可靠的定义，得到相应的可靠度指标，通过迭代法求解真正失效点的可靠度指标βi。最后根据3个状态方程相互独立的假定，得出整个结构的可靠度指标β。该方法在平川泥石流防治工程中得到了成功的运用。这对于可靠度设计方法在泥石流防治工程中的运用进行了有益的探讨。     

The RUSLE erosion index as a proxy indicator for debris flow susceptibility

Debris flows represent dangerous occurrences in many parts of the world. Several disasters are documented due to this type of fast-moving landslides; therefore, natural-hazard assessment of debris flows is crucial for safety of life and property. To this aim, much current work is being directed toward developing geotechnical-hydraulic models for the evaluation of debris flow susceptibility. A common base for such current models is parameterization of background predisposing and triggering factors such as inherent characteristics of geo-materials, topography, landscape and vegetation cover, rainfall regime, human activities, etc. which influence the occurrence of these processes on slopes. The same factors are also taken into account in soil erosion prediction models. Consequently, it seems worth investigating the effectiveness of the soil erosion index as debris flows susceptibility indicator. To this aim, a logistic regression analysis was carried out between the erosion index assessed by means of the Revised Universal Soil Loss Equation (RUSLE) model and the inventory of debris flows that have occurred in an area in Sicily (Southern Italy). Model assumptions were verified and validated by means of a series of statistical tools. Different possible scenarios were also evaluated by considering hypothetical changes in soil erosion rate under different rain erosivity conditions. Notwithstanding the rough approximations in model data collection, the outcomes appear encouraging.     

基于不同侵蚀模型的高速崩滑碎屑流动力过程模拟分析

为实现高速崩滑碎屑流沿程侵蚀动力过程的模拟分析,采用连续理论方法对NomashRiver碎屑流动力过程进行了数值模拟。其中,在连续理论模型中采用3种不同侵蚀速率模型,并采用HLLC近似Riemann解对有限体积数值离散控制体单元的界面通量进行了计算求解。致灾范围及运动时间的计算结果均与实际灾害情况吻合良好,验证了计算模拟的正确有效性,并对最终堆积深度、运动速度和侵蚀区域侵蚀深度进行了分析讨论。结果显示:采用McDougall侵蚀模型得到的最终堆积平均深度和最大深度与实测情况较为接近;每个时刻采用Medina侵蚀模型得到的最大速度值最大,其次是采用McDougall侵蚀模型的结果,最后是采用Pitman侵蚀模型的结果;采用McDougall侵蚀模型得到的侵蚀深度分布较为连续,其最大值8.1m与估测值8m比较接近,采用Medina侵蚀模型和Pitman侵蚀模型得到的侵蚀深度结果则较为分散,其最大值分别为10.9m和8.6m。     

基于HBP本构模型的泥石流动力过程SPH数值模拟

本构模型是描述泥石流流变特性的关键,也是决定其动力过程数值模拟准确性的核心问题之一。泥石流流体属多相混合物,现有的研究已证实其存在剪切增稠或剪切变稀的现象,传统基于Bingham及Cross线性本构关系的数值模型难以准确描述泥石流流变特性。文中探讨了Bingham模型在低剪应变率下的数值发散问题,在光滑粒子流体动力学（SPH）方法框架上建立了整合Herschel-Bulkley-Papanastasiou(HBP）本构关系的稀性泥石流动力过程三维数值模型。相比传统基于浅水波假设的二维数值模型,所述方法从三维尺度建立SPH形式下的泥石流浆体纳维?斯托克斯方程并进行数值求解,可获取泥石流速度场时空分布及堆积形态,同时采用HBP本构关系描述泥石流流变特性,能在确保数值收敛的前提下反映泥石流流体在塑性屈服过渡段及大变形状态下应力?应变的非线性变化。为验证提出方法的合理性,结合小型模型槽实验观测进行了对比,结果表明数值模拟与实测结果基本吻合。     

泥石流成因机理的非饱和土力学理论研究

泥石流是一种具有较强破坏力的自然山地灾害。对于它的预报研究历来为人们所重视，并建立了很多雨量预报模型。然而，这些雨量预报模型的预报时间很短，往往只能在灾害发生前几十分钟作出预报。论文应用非饱和土强度理论对降雨型泥石漉的成园机理进行了研究，提出降雨型泥石流的形成过程可以划分为2个阶段：第一个阶段与前期实效降雨量有关；第二个阶段与短历时强降雨有关。并对各个阶段降雨作用机理以及固体松散物质的力学性质变化特征进行探讨。为预先判断在降雨条件下，会不会发生泥石流以及所需要的降雨量和雨型提供依据。     

Residual luminescence in modern debris flow deposits from western Taiwan: A single grain approach

In an active fold-and-thrust belt, sediment is commonly transported in the form of debris and colluviums, and the ages of such deposits can be important for seismic hazard assessment. These sediments are mainly transported by high concentration flows that travel very short distances before burial, and therefore have very few opportunities for daylight exposure and bleaching, which limits the application of luminescence dating in such environments. Luminescence studies on modern debris flow deposits and channel sediments in central Taiwan indicate a mixed population of bleached and unbleached grains. The well bleached grains constitute a small proportion of the sediments and could only be identified using the single grain dating technique; doing single aliquot analysis may result in overestimation of the age of such sediments by 3–25 ka. The data reveal that outwash sediments are better bleached than debris flow sediments. Our results show that quartz that has experienced partial resetting can be used for luminescence dating only through the use of the single grain technique.     

Substantiation of debris flow velocity from super-elevation: a numerical approach

The use of super-elevations that a forced vortex flow leaves on the valley walls of a curved flume is a plausible approach toward estimating debris flow velocities in earthquake-induced geo-hazard studies. The centrifugal force of a speeding flow is responsible for a higher flow depth on the outer bend. However, in reality, a flow is not steady, and only the highest flow-marks are left at the outer and inner bends of the flow, which can lead to an inaccurate estimation of the actual velocity. Seeing the real scenario of the field, a series of numerical flume tests using smoothed particle hydrodynamics (SPH) is conducted to validate the estimation of debris flow velocities from flow-marks. Velocities estimated from flow-marks are lower than real velocities near the source region, but they converge to real velocities as the distance to the source increases. Based on several simulations, a best-fit line is proposed for adjusting debris flow velocity from mud-marks, and it is used to estimate flow velocities of the well-documented debris event called “Shiraito river debris flow,” which happened near the rim of the Hakone Crater, Kanagawa Prefecture, Japan, ensuing from the 1923 Great Kanto earthquake.     

Assessment of dynamic impact force of debris flow in mountain torrent based on characteristics of debris flow

Landslides and debris flows that occur around residential areas are considered, globally, as significant disasters that cause damage to human life and property. With terrain slope defining the flow characteristics of debris flows, flow depth, flow velocity, and impact force vary by time and distance. In particular, when a structure is located in the flow path of debris flows, the flow characteristics of debris flows vary by terrain slope and direction angle. To simulate the flow characteristics of these debris flows, the simulation results obtained by FLO-2D were analyzed with six-stage conditions for the research area. In the analysis, the flow depth, flow velocity, and impact force were estimated on the basis of the outlet of the research area in the presence and absence of structure(s) at certain distances. With this, the variation of the impact force in accordance with the variation of the flow depth of the debris flows was highly similar to the simulation results obtained by FLO-2D, when the correction index (α) of the suggested dynamic impact force equation was 0.3–0.4. There were sections where the estimated value of the impact force was overestimated near the outlet, and it was judged that the fixed values of the terrain factors (width, roughness coefficient, slope, etc.) caused the impact force to be overestimated. However, the correlation analysis showed that the correlation index was above the normal ranges in the suggested dynamic impact force equation for debris flows with the application of the terrain factors.