Simulation of two-phase debris flow scouring bridge pier

Debris flows are typical two-phase flows, which commonly accompany erosion in mountainous areas, and may destroy bridge engineering by scouring. In this study, a physically-based two-phase model is applied for the simulation of debris flow scouring of bridge pier. In this model, the shear stress of debris flow on an erodible bed is considered to be a function of the solid shear stress, fluid shear stress, and volume fraction; accordingly, the erosion is incorporated into the two-phase model. Using a highaccuracy computational scheme based on the finite volume method, the model is employed for simulating a dynamic debris flow over an erodible bed. The numerical results are consistent with the experimental data, and verify the feasibility of the two-phase model. Moreover, a simple numerical test is performed to exhibit the fundamental behaviour of debris flow scouring of bridge pier, which shows that the degree of erosion on each side of the pier is higher compared to other areas. The scouring depth is influenced by the variations of solid volume fraction and velocity of debris flow and pier width.     

A quasi single-phase model for debris flows and its comparison with a two-phase model

A depth-averaged quasi single-phase mixture model is proposed for debris flows over inclined bed slopes based on the shallow water hydrosediment-morphodynamic theory with multi grain sizes. The stresses due to fluctuations are incorporated based on analogy to turbulent flows, as estimated using the depth-averaged k-? turbulence model and a modification component. A fully conservative numerical algorithm, using wellbalanced slope limited centred scheme, is deployed to solve the governing equations. The present quasi single-phase model using four closure relationships for the bed shear stresses is evaluated against USGS experimental debris flow and compared with traditional quasi single-phase models and a recent physically enhanced two-phase model. It is found that the present quasi single-phase model performs much better than the traditional models, and is attractive in terms of computational cost while the two-phase model performs even better appreciably.     

Velocity profile of debris flow based on quadratic rheology model

The quadratic rheology model considers the yield stress, viscous stress, turbulent stress and disperse stress, so it is used in this study to derive the velocity profile of debris flows. The quadratic model with the parabolic eddy viscosity was numerically solved, and an analytical solution was derived for the quadratic model with a constant eddy viscosity. These two solutions were compared with the Arai-Takahashi model that excluded the viscous stress and the yield stress. The three models were tested by using 17 experiment cases of debris flows over rigid beds. The results prove that the quadratic model with parabolic and constant eddy viscosities is applicable to muddy and granular flows, whereas the Arai-Takahashi model tends to overestimate the flow velocity near the water surface if a plug-like layer exists. In addition, the von Karman constant and the zero-velocity elevation in the three models are related to sediment concentration. The von Karman constant decreases first and then increases as the sediment concentration increases. The zero-velocity elevation is below the bed surface, likely due to the invalidity of the non-slip boundary condition for the debris flows over fixed beds.     

降雨模式对震后泥石流起动模式影响的试验研究: 以九寨天堂沟为例

降雨过程中降雨强度的变化会影响土体渗透率及饱和过程, 从而改变土体的力学性质, 影响泥石流起动模式及破坏规模。为探究不同降雨模式对震后泥石流起动机制的影响, 自制了小比例模型槽, 结合可控雨型的降雨模拟系统, 进行了人工降雨诱发泥石流的室内模型试验; 基于不同降雨模式下泥石流的起动过程分析, 对坡体内部含水率和孔隙水压力的变化规律进行了研究。研究结果表明: 递增型降雨模式下泥石流发生突然, 呈整体滑坡转化为泥石流起动模式, 坡体破坏规模最大; 递减型降雨模式下表现为后退式溃散失稳起动模式; 均匀型降雨模式下则表现为溯源侵蚀起动模式; 中峰型降雨模式下以局部滑坡转化为泥石流起动模式; Ⅴ型降雨模式下则由坡面侵蚀加剧转化为泥石流启动模式, 破坏规模最小。研究结果可以为九寨沟地区泥石流的预报预警提供参考。      

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.     

泥石流危险范围预测模型及在昆明东川城区的应用

结合泥石流危险范围模型实验数据,运用多元回归分析方法探讨了泥石流危险范围预测,并进行了误差分析。以昆明市东川城区后山3条泥石流沟为例,运用该模型对其危险范围进行了预测分析,为东川城区泥石流防灾提供了科学依据。     

RECENT DEVELOPMENTS IN DEBRIS FLOW RESEARCH IN ITALY

ＲＥＣＥＮＴＤＥＶＥＬＯＰＭＥＮＴＳＩＮＤＥＢＲＩＳＦＬＯＷＲＥＳＥＡＲＣＨＩＮＩＴＡＬＹＭａｒｃｈｉＬｏｒｅｎｚｏ；ＴｅｃｃａＰｉａＲ．（ＩｎｓｔｉｔｕｔｅｆｏｒＰｒｅｖｅｎｔｉｏｎｏｆＨｙｄｒｏｌｏｇｉｃａｌａｎｄＧｅｏｌｏｇｉｃａｌＨａｚａｒｄ...     

Magnitude-frequency relationship of debris flows in the Jiangjia Gully,China

The magnitude-frequency(MF) relationship of debris flows is the basis for engineering designs and risk quantification. However, because of the lack of debris flow monitoring data, research progress in this area has been relatively slow. The MF relationship of debris flows in Jiangjia Gully, Yunnan Province was evaluated based on a regression analysis of 178 debris flow events that occurred from 1987-2004. The magnitude-cumulative frequency(MCF) relationship of the debris flows in the Jiangjia Gully is consistent with the linear logarithmic transformation function. Moreover, observed data for debris flows in Hunshui Gully of Yunnan Province and Huoshao Gully, Liuwan Gully, and Niwan Gully of Gansu Province were used to verify the function. The results showed that the MCF relationship of highfrequency debris flows is consistent with the power law equation, although the regression coefficients in the equation are considerably different. Further analysis showed a strong correlation between the differences in the constants and the drainage area and daily maximum precipitation.     

Hybrid simulation of the initiation and runout characteristics of a catastrophic debris flow

On 13 August 2010, a catastrophic debris flow with a volume of 1.17 million m³ occurred in Xiaojiagou Ravine near Yingxiu town of Wenchuan county in Sichuan Province, China. The main source material was the landslide deposits retained in the ravine during the 2008 Wenchuan earthquake. This paper describes a two-dimensional hybrid numerical method that simulates the entire process of the debris flow from initiation to transportation and finally to deposition. The study area is discretized into a grid of square zones. A two dimensional finite difference method is then applied to simulate the rainfall-runoff and debris flow runout processes. The analysis is divided into three steps; namely, rainfall-runoff simulation, mixing water and solid materials, and debris flow runout simulation. The rainfall-runoff simulation is firstly conducted to obtain the cumulative runoff near the location of main source material and at the outlet of the first branch. The water and solid materials are then mixed to create an inflow hydrograph for the debris flow runout simulation. The occurrence time and volume of the debris flow can be estimated in this step. Finally the runout process of the debris flow is simulated. When the yield stress is high, it controls the deposition zone. When the yield stress is medium or low, both yield stress and viscosity influence the deposition zone. The flow velocity is largely influenced by the viscosity. The estimated yield stress by the equation, τ _y = ρghsin θ, and the estimated viscosity by the equation established by Bisantino et al. (2010) provide good estimates of the area of the debris flow fan and the distribution of deposition depth.     

Geomorphic characteristics of hillslope and channelized debris flows: A case study in the Shitou area of central Taiwan

The data on the hillslope and channelized debris flows in the Shitou area of central Taiwan occurred during Typhoons Toraji and Nali in 2001 were applied in this paper. The geomorphic parameters, including the flow length, gully gradient, drainage area and form factor of the debris flows were determined by spatial analysis using a Geographic Information System （GIS） based on the data derived from field investigation, aerial photographs, and topographical maps. According to such determined geomorphic parameters, the threshold conditions and empirical equations, such as the relationship between the gully gradient and drainage area and that between gully length and drainage area and topographic parameter, are presented and used to distinguish the geomorphic characteristics between the channelized and hillslope debris flows.     

Effects of Roughness Elements Distribution on Overland Flow Resistance

Roughness elements are various in a mountain area; they include gravel and ground surface vegetation that often result in surface friction drag to resist overland flows. The variation and characteristics of flow resistance strongly impact the overland flow process and watershed floods. In view of the universal existence of natural vegetation, such as Chlorophytum malayense(CM) or Ophiopogon bodinieri(OB), and the sand-gravel bed of the river channel, it is important to understand the role of different types of roughness elements in flow resistance. This study was performed to investigate and compare through flume experiments the behaviors of overland flow resistance by the reaction of multi-scale configuration of different roughness elements. The result showed that the resistance coefficient gradually reduced versus the increase of flow rate in unit width and tended to be a constant when q = 3.0 l/s.m, Fr = 1.0, and Re = 4000 for slopes of 6 to 10 degrees. The gap of the vegetated rough bed and the gravel rough bed is limited to the same as the gap of the two types of vegetation, CM and OB. It was noted that the vegetation contributed to the increase in form resistance negatively and may lead to the mean resistance on decrease. To classify the flow pattern, the laminar flows were described by DarcyWeisbach's equation. In the study the f-Re equation of vegetated bed was developed with f ?5000 Re.The friction coefficient for laminar flows can be regarded as the critical value for identifying the transformation point of the flow pattern.     

Changes in runout distances of debris flows over time in the Wenchuan earthquake zone

A large number of debris flows occurred in the Wenchuan earthquake zone after the 12 May 2008 earthquake.The risks posed by these debris flows were rather high.An appropriate model is required to predict the possible runout distance and impacted area.This paper describes a study on the runout characteristics of the debris flows that occurred in the Wenchuan earthquake zone over the past four years.A total of 120 debris flows are analyzed.Separate multivariate regression models are established for the runout distances of hill-slope debris flows and channelized debris flows.The control variables include type of debris flow,debris flow volume,and elevation difference.Comparison of the debris flows occurring before and after the earthquake shows that the runout distance increased after the earthquake due to sufficient material supply and increased mobility of the source materials.In addition,the runout distances of annual debris flow events in 2008,2010 and 2011 are analyzed and compared.There is a tendency that the runout distance decreases over time due to the decreasing source material volume and possible changes of debris flow type.Comparison between the debris flows in the earthquake zone and the debris flows in Swiss Alps,Canada,Austria,and Japan shows that the former have a smaller mobility.     

Effects of loose deposits on debris flow processes in the Aizi Valley,southwest China

Loose deposits, rainfall and topography are three key factors that triggering debris flows.However, few studies have investigated the effects of loose deposits on the whole debris flow process.On June 28, 2012, a catastrophic debris flow occurred in the Aizi Valley, resulting in 40 deaths.The Aizi Valley is located in the Lower Jinsha River,southwestern Sichuan Province, China. The Aizi Valley debris flow has been selected as a case for addressing loose deposits effects on the whole debris flow process through remote sensing, field investigation and field experiments. Remote sensing interpretation and laboratory experiments were used to obtain the distribution and characteristics of the loose deposits, respectively. A field experiment was conducted to explore the mechanics of slope debris flows, and another field investigation was conducted to obtain the processes of debris flow formation, movement and amplification. The results showed that loose deposits preparation, slope debris flow initiation,gully debris flow confluence and valley debris flow amplification were dominated by the loose deposits.Antecedent droughts and earthquake activities may have increased the potential for loose soil sources in the Aizi Valley, which laid the foundation for debris flow formation. Slope debris flow initiated under rainfall, and the increase in the water content as well as the pore water pressure of the loose deposits were the key factors affecting slope failure. The nine gully debris flows converged in the valley, and the peak discharge was amplified 3.3 times due to a blockage and outburst caused by a large boulder. The results may help in predicting and assessing regional debris flows in dry-hot and seismic-prone areas based on loose deposits, especially considering large boulders.     

Influence of housing and urban development on debris flow flooding and deposition

Debris flows form deposits when they reach an alluvial fan until they eventually stop.However,houses located in the alluvial fan might affect the debris flow flooding and deposition processes.Few previous studies have considered the effects of houses on debris flow flooding and deposition.This study conducted model experiments and numerical simulations using the Kanako2D debris flow simulator to determine the influence of houses on debris flow flooding and deposition.The model experiments showed that when houses are present,the debris flow spreads widely in the cross direction immediately upstream of the houses,especially when the flow discharge is large or the grain size is small.Houses located in the alluvial fan also influence the deposition area.The presence of houses led to flooding and deposition damage in some places and reduced the damage in others.The simulation also demonstrated the influence of houses.Both the model experiment and the simulation showed that houses change the flooding and deposition areas.     

Sediment movement in periodic alternating current

The present paper summarizes the results of previous studies, including the structure and principle of the rotary ring flume for researching the fine sediment movement, the mechanism of the flume, method of eliminating the influence of the centrifugal force and sediment movement experiments with periodic alternating current. Also included are the experiment-based relationship among sediment concentration, bed shear stress and silt carrying capacity, a proposed erosion-deposition function and bed erosion-deposition calculation together with the results of verification.     

Geomorphologic analysis and physico-dynamic characteristics of Zhatai-Gully debris flows in SW China

Zhatai gully is a typical debris flow channel in Butuo county of Sichuan province, southwestern China. The geomorphologic features are analyzed and the physical-dynamic characteristics are discussed on the basis of field investigation and laboratory tests. Geomorphologic analysis indicates that Zhatai-gully drainage in relation to debris flow can be divided into source area, transport area, and deposition area. The source area has a steep slope and has very limited vegetation cover, which favors runoff, allowing loose solid materials to be mobilized easily and rapidly. In the transport area, there are many small landslides, lateral lobes, and loose materials distributed on both banks. These landslides are active and constantly providing abundant source of soils for the debris flows. In the deposition area, three old debris-flow deposits of different ages can be observed. The dynamic calculation shows that within the recurrence intervals of 50 and 100 years, debris flow discharges are 155.77m³/s and 1y8.19m³/s and deposition volumes are 16.39 x 10⁴ m³ and 18.14 x 10⁴ m³, respectively. The depositional fan of an old debris flow in the outlet of the gully can be subdivided into six layers. There are three debris flow deposits on left and two on the right side of the gully. Grain-size tests of sediments from the soil, gulley bed deposits, and the fresh and old debris flow deposits showed that high amounts of clay and fine gravel were derived from the soil in the source area whereas much of the gravel fraction were sourced from the gully bed deposits. Comprehensive analysis indicates that Zhatai gully is viscous debris-flow gully with moderate to high frequency and moderate to large magnitude debris flows. The risk of a debris flow disaster in Zhatai-gully is moderate and poses a potential threat to the planned hydroelectric dam. Appropriate engineering measures are suggested in the construction and protection of the planned hydroelectric station.     

Influences of the Wenchuan Earthquake on sediment supply of debris flows

The 5.12 Wenchuan Earthquake and the subsequent rainstorms induced a large number of landslides, which later were transformed into debris flows. To evaluate the effect of the earthquake on the sediment supply of debris flows, eight debris flow basins near Beichuan City, Sichuan Province, China were chosen as the study area. The area variations of the debris flow source after the Wenchuan Earthquake and the subsequent rainstorm are analyzed and discussed in this paper. Interpretations of aerial photographs (after the 5.12 Wenchuan Earthquake) and SPOT5 images (after the rainstorm event of September 24, 2008) as well as field investigations were compared to identify the transformation of landslide surface in the study area, indicating that the landslide area in the eight debris flow basins significantly increased. The loose sediment area on the channel bed increased after the rainstorm event. In order to estimate the relationship of the landslide area with the rainfall intensity in different return periods, a model proposed by Uchihugi was adopted. Results show that new landslide area induced by heavy rainfall with 50-year and 100-year return period will be 0.87 km2 and 1.67 km2, respectively. The study results show the Wenchuan earthquake had particular influences on subsequent rainfall-induced debris flow occurrence.     

Risk assessment of highways affected by debris flows in Wenchuan earthquake area

Since the Wenchuan earthquake in China on May 12^th, 2008, highways in earthquake-affected areas have been frequently interrupted by debris flows. We analyzed the hazard effect modes and damage processes along highways and developed three key indexes, scale of debris flows, deposits on highways and river blockage, to describe quantitatively the highway disasters. By combining the empirical methods and the actual terrain conditions, we proposed new methods to determine the value of hazard indexes. In addition, we used the economic value and resistance of highway as vulnerability assessment indexes, then determined the specific subindexes for the subgrade, bridges and culverts, and developed a way for the quantified vulnerability zoning. Moreover, we proposed the assessment and mapping methods for highway risk. The risk is described into 5 grades: extremely low risk, low risk, middle risk, high risk and extremely high risk. We applied these methods in a case study carried out on provincial highway S303 from Yingxiu Town to Wolong Town, in Wenchuan County. Analysis of debris flow risk for the whole highway, showed that the total length of highway in extremely low risk area was 28.26 km, 4.83 km in low risk area, 8.0 km in middle risk area, 3.65 km in high risk area, and 3.06 km in extremely high risk area. The assessment results are consistent with the field survey data which reflected the disaster situation. This risk method can be used objectively to evaluate the debris-flow risk along highways, and is useful for highway reconstruction in mountainous areas suffering from active debris flows.     

Debris flow susceptibility assessment using a probabilistic approach: A case study in the Longchi area,Sichuan province,China

The Longchi area with the city of Dujiangyan, in the Sichuan province of China, is composed of Permian stone and diorites and Triassic sandstones and mudstones intercalated with slates. An abundance of loose co-seismic materials were present on the slopes after the May 12, 2008 Wenchuan earthquake, which in later years served as source material for rainfall-induced debris flows or shallow landslides. A total of 48 debris flows, all triggered by heavy rainfall on 13th August 20l0, are described in this paper. Field investigation, supported by remote sensing image interpretation, was conducted to interpret the co-seismic landslides in the debris flow gullies. Specific characteristics of the study area such as slope, aspect, elevation, channel gradient, lithology, and gully density were selected for the evaluation of debris flow susceptibility. A score was given to all the debris flow gullies based on the probability of debris flow occurrence for the selected factors. In order to get the contribution of the different factors, principal component analyses were applied. A comprehensive score was obtained for the 48 debris flow gullies which enabled us to make a susceptibility map for debris flows with three classes. Twenty-two gullies have a high susceptibility, twenty gullies show a moderate susceptibility and six gullies have a low susceptibility for debris flows.     

Debris flows introduced in landslide deposits under rainfall conditions: The case of Wenjiagou gully

Debris flows often occur in landslide deposits during heavy rainstorms. Debris flows are initiated by surface water runoff and unsaturated seepage under rainfall conditions. A physical model based on an infinitely long, uniform and void-rich sediment layer was applied to analyze the triggering of debris-flow introduced in landslide deposits. To determine the initiation condition for rainfall-induced debris flows, we conducted a surface water runoff and saturated-unsaturated seepage numerical program to model rainfall infiltration and runoff on a slope. This program was combined with physical modeling and stability analysis to make certain the initiation condition for rainfall-introduced debris flows. Taking the landslide deposits at Wenjiagou gully as an example, the initiation conditions for debris flow were computed. The results show that increase height of surface-water runoff and the decrease of saturated sediment shear strength of are the main reasons for triggering debris-flows under heavy rainfall conditions. The debris-flow triggering is affected by the depth of surface-water runoff, the slope saturation and shear strength of the sediment.