Next Article in Journal
Impact of Land Use Change on the Water Environment of a Key Marsh Area in Vientiane Capital, Laos
Next Article in Special Issue
Land Subsidence Due to Groundwater Exploitation in Unconfined Aquifers: Experimental and Numerical Assessment with Computational Fluid Dynamics
Previous Article in Journal
Soil-Matric-Potential-Based Irrigation Scheduling to Increase Yield and Water Productivity of Okra
Previous Article in Special Issue
Comparative Analysis of Tree-Based Ensemble Learning Algorithms for Landslide Susceptibility Mapping: A Case Study in Rize, Turkey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification Method of River Blocking by Debris Flow in the Middle Reaches of the Dadu River, Southwest of China

1
College of Earth Science and Engineering, Shandong University of Science and Technology, Qingdao 266590, China
2
School of Tourism and Urban-Rural Planning, Chengdu University of Technology, Chengdu 610081, China
3
Chengdu Center of China Geological Survey, Chengdu 610081, China
4
College of Water Resource and Hydropower, Sichuan University, Chengdu 610065, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(24), 4301; https://doi.org/10.3390/w15244301
Submission received: 26 October 2023 / Revised: 11 December 2023 / Accepted: 13 December 2023 / Published: 18 December 2023
(This article belongs to the Special Issue Risk Analysis in Landslides and Groundwater-Related Hazards)

Abstract

:
Debris flow is a typical natural disaster in the middle reaches of the Dadu River, which seriously threatens the safety of life and property of local residents. However, there is currently a lack of a comprehensive analysis methods applicable to the blockage of river channels by debris flow in the Dadu River basin, limiting disaster prevention and mitigation in this area. Based on previous large-scale model tests carried out in the middle reaches of the Dadu River, the debris flows are divided into dam-type debris flows and submerged debris flows. The calculation formulas for the maximum travel distance of the two kinds of debris flows entering the river are obtained via theoretical derivation. The formulas for calculating the length and volume of debris flow accumulation are derived, and the relationship between the debris flow loss coefficient and river blocking degree in the middle part of the Dadu River is analyzed. An identification method of river blocking by debris flow is put forward in this study. By calculating the maximum blocking degree, S (the ratio of the maximum driving distance of the debris flow to the width of the river), and the volume of the source materials needed to form a debris flow dam under the conditions that the debris flow does not reach the opposite bank (V1), reaches the opposite bank but does not block the river (V2), and reaches the opposite bank (V3), the form of debris flow blocking the river is distinguished. When S = 1, V > V3, complete blockage occurs; when S = 1, V > V2, the river is mostly blocked; when S < 1, V > V1, the river is half-blocked. This study established an identification method of river blocking by debris flow, providing a basis for early warning for river blocking and disaster prevention in the middle reaches of the Dadu River.

1. Introduction

Debris flow is a special type of torrent that carries a large amount of sediment, rocks, and other materials, and that is caused by short-term heavy precipitation in mountainous areas and snow-melting [1,2,3,4]. Furthermore, after debris flows enter rivers, they can trigger secondary disasters, such as dammed lakes, which can threaten hydropower stations, transportation infrastructure, and the lives and properties of local residents. Debris flow disasters have caused enormous losses worldwide. In May 1970, a debris flow disaster occurred in the town of Yange, Peru. The flow speed was nearly 100 m/s, and 5000 × 104 m3 of solid material was washed away, resulting in nearly 20,000 deaths [5]. On 15 December 1999, heavy rain hit all parts of Venezuela, causing dozens of gullies to develop large-scale debris flows at the same time, and resulting in 30,000 casualties and enormous economic losses of up to 10 billion US dollars [6]. On 19 August 2019, Wenchuan County in the Sichuan Province of China suffered multiple debris flow disasters caused by heavy precipitation; approximately 446,000 people were affected, and this event caused nearly 15.89 billion yuan of direct economic losses [7,8].
The Dadu River is located in northwestern Sichuan Province. The middle reaches of the Dadu River have deep river valleys, strong tectonic activity, and complex stratigraphic lithology, resulting in the frequent occurrence of debris flow disasters in this region. On 23 July 2009, a debris flow occurred in Xiangshui gully, Kangding city, and completely blocked the Dadu River, causing 16 deaths, 38 missing, and considerable property losses [9]. On 12 July 2012, a debris flow occurred in Tangjia gully, Shimian County, Yaan city, and blocked the Dadu River, resulting in two deaths and five missing persons [10]. On 4 July 2013, the Xiongjia gully debris flow completely blocked the Zhuma River and caused 19 deaths [11].
Hence, identifying the risk of river blocking by debris flows is important for disaster prevention and mitigation in the middle reaches of the Dadu River. At present, there are four main methods for studying the mechanism of river blocking by debris flows: (1) experimental studies using flume tests [12,13,14,15,16]; (2) statistical models based on field investigations [17,18,19]; (3) numerical simulations; and (4) theoretical analysis [20,21,22]. However, there are few practical studies on the identification of river blocking by debris flows in the Dadu River basin, which serves as one of the major hydropower bases in China. Therefore, a more applicable and comprehensive identification method is needed for the Dadu River to reduce the loss of life and property caused by debris flow disasters in the basin.
In this study, the interaction mechanism between the debris flow and the water flow after the debris flow enters the river is analyzed based on the field investigation in the middle reaches of the Dadu River, and the debris flows occurred in this area are divided into a barrier type and a submerged type according to whether the debris flow into the river is below the water surface and the characteristics of the debris flow. This study takes the three-dimensional spatial relationship between tributaries and main streams into account, formulas for calculating the travel distance of different types of debris flows are established, and finally, a comprehensive identification method for debris flow interception under different models is proposed.

2. Study Area

The study area is the middle reaches of the Dadu River (Kangding–Shimian section), covering an area of 6183.33 km2. The study area is crisscrossed by water systems and numerous tributaries, such as the Moxi River, Wasigou, Nanya River, Anshun River, and Tianwan Gou. The valleys are deep, the maximum elevation difference is 5200 m, as illustrated in Figure 1. The middle reaches of the Dadu River are located in the transition zone from the Sichuan Basin to the Qinghai–Tibet Plateau. The middle reaches of the Dadu River are divided into three sections, i.e., the Kangding section, the Luding section, and the Shimian section. The middle reaches of the Dadu River are located in the transition zone from the Sichuan Basin to the Qinghai–Tibet Plateau. Affected by the southeast and southwest monsoon and the cold air of Qinghai–Tibet Plateau, the vertical climate difference is obvious. The areas below 1800 m altitude in the study area have tropical monsoon climate and experience local heavy rainfall. There are six tributaries in this area, and each tributary may have different confluence characteristics. The average annual total flow in this section is 330.56 × 108 m3, the average annual flow is 1218 m3/s, and the maximum peak flow is 6050 m3/s. The location of the study area is shown in Figure 1.
The middle section of the main stream of the Dadu River flows along the Dadu River fault zone, which is a zone of abrupt changes in topography and landforms that constitutes the eastern boundary of the Qinghai–Tibet Plateau. The middle reaches of the Dadu River are located on the eastern margin of the Qinghai–Tibet Plateau, and strata other than Cretaceous and Cambrian strata are exposed. Macroscopically, the middle reaches of the Dadu River are located at a tectonic intersection with complex structures. The north side is the Sichuan–Qinghai block, the west side is the Sichuan–Yunnan block, and the east side is the South China block. The NW-trending Xianshuihe fault, the NE-trending Longmenshan fault and the N–S-trending Anninghe fault constitute a Y-shaped structure in the Shimian section.
Field investigation shows that there are a total of 441 debris flows in the study area, and these debris flows can be divided into four categories based on the area of the debris flow deformation region, i.e., extra-large, large, medium and small (Figure 2).

3. Type of Interaction between the Debris Flow in a Branch Channel and the Main River

Debris flows tend to accumulate on flatter terrain, and this is caused by two factors. First, when the debris flows rush out from a gully, the gravitational potential energy is reduced because the topographic slope decreases; second, the debris flow usually spreads in a fan shape after it stops, leading to an increase in friction. As the middle reaches of the Dadu River have flat topography, debris flows are usually deposited in the river channel. After a debris flow enters the main river, some source materials may not flow downstream together with the river since the debris flow deposits too quickly or the river does not have sufficient capacity to scour and transport the material [23,24]. These source materials can resist river erosion and be deposited in the river channel, thus gradually forming a natural dam and eventually leading to river blockage. This is defined as a dam-type debris flow in this study, as illustrated in Figure 3a. In addition, if the water level of the main river is high, making it impossible for a debris flow entering the main river to form a dam-type debris flow, the debris flow will rush into the main river and be submerged, forming a submerged-type debris flow, as illustrated in Figure 3b.
The flow field of the debris flow is disordered, and the intermixing effect is difficult to describe quantitatively. To facilitate its theoretical description, the mixing effect at the interface between debris flow and water flow is often ignored. Instead, it is considered a Newtonian interaction with different densities. Considering the different viscosities of the debris flow, it is difficult for turbulent flow at the interface between water and the debris flow to disrupt the fluid structure of the high-viscosity debris flows because debris flows with higher viscosity have higher structural force. However, debris flows with higher structural forces can be easily impacted by water flow due to their poor structure, resulting in more intense mixing. In practice, compared with viscous debris flows, dilute debris flows are less likely to block rivers due to their susceptibility to scouring [25]. Therefore, only cohesive debris flows are considered in this study.
For viscous debris flows, there is an interface between the debris flow and the water flow in the river, as the material exchange between them is ignored (Figure 3). When the debris flow rushes into the river, it is subjected to a bypass resistance at the debris flow head and a frictional resistance from the riverbed at the bottom surface of the debris flow. When the debris flow enters the river in a submerged manner, its top surface is subject to gravitational and shear forces from the river flow, the specific direction of which varies depending on the direction of their mutual movement.

4. Calculation Model for Dam-Type Debris Flow

4.1. Equations of the Motion of a Debris Flow in a Variable-Slope Channel

Takahashi derived the equation for the motion of debris flows in variable-slope channels based on the conservation of momentum theorem [26].
d d t [ 1 2 ( h + h f ) y γ d g v d ] = 1 2 ( h + h f ) y γ d sin θ + γ d g q T v u cos ( θ u θ ) + 1 2 g h u 2 cos θ u cos ( θ u θ ) [ ( σ p ) c d u K a + ρ ] F
where t refers to the time; g represents the gravitational acceleration; h and h f stand for the simplified shape parameters for debris flow head in y–z two-dimensional plane (Figure 4); γ d is the gravity of debris flow; ρ σ refers to the density of sand; ρ means the density of water; c d u is the volume specific concentration of the debris flow; F is the ground friction force on the debris flow; and K a = tan 2 ( 45 ° ϕ / 2 ) ( ϕ is the angle of dynamic friction of the debris flow materials) is the coefficient of the active earth pressure.
Without considering the change in direction due to the x-directional flow pressure (Figure 5), the equation can be used to calculate the y-directional movement distance of the debris flow in the channel with a variable slope.
After the debris flow enters the river, its head and body are resisted by both the riverbed and the river flow. Without considering the material exchange between the river flow and the debris flow, the water pressure at the confluence is simplified to hydrostatic pressure, a hydrostatic pressure of the river is added to Equation (1), and the kinetic equation of the motion of debris flow can be obtained, as follows:
d d t [ 1 2 ( h + h f ) y γ d g v d ] = 1 2 ( h + h f ) y γ d sin θ + γ d g q T v u cos ( θ u θ ) + 1 2 g h u 2 cos θ u cos ( θ u θ ) [ ( σ p ) c d u k a + ρ ] F f
where f is the resistance of the river flow to the debris flow, which is determined as follows:
f = C D h u ρ v u 2 8 cos θ + 1 2 g h u 2 ρ cos θ
where the first term on the right side is the drag force of the debris flow when moving in the water, and the second term on the right side is the hydrostatic pressure. Among them, the value of the drag coefficient, CD, is related to the shape of the debris flow and the Reynolds number of the main river flow, which is generally in the range of 0.5–2.5 in engineering practices. Liu considered the shape of the debris flow head and recommended the values of the drag force coefficient, as listed in Table 1 [27].
From the law of conservation of mass, it follows that
d d t [ 1 2 ( h + h f ) y ] = q T
Integrating this and combining the initial value condition x ( t = 0 ) = 0 , Equation (2) gives
d v d t = v t + Q t R
where
R = ( γ d γ w ) g c d u n cos θ ( γ d γ w ) c d u + γ w Q = v u cos ( θ u θ ) + cos θ u cos ( θ u θ ) [ ( γ d γ w ) c d u k a + γ w ] g h u 2 [ ( γ d γ w ) c d u + γ w ] v u γ w cos θ g h u 2 [ ( γ d γ w ) c d u + γ w ] v u C D γ w v u cos θ 8 [ ( γ d γ w ) c d u + γ w ]
Integrating Equations (1)–(4) from 0 to t and substituting the boundary conditions yields v(t = 0) = 0, giving:
y = 1 4 R t 2 + Q t
Substitute yields the formula for the maximum transport distance, as follows:
L = Q 2 R

4.2. Analysis of Debris Flow Movements Considering the Spatial Relationship of the Channel

4.2.1. Maximum Travel Distance at the Cross-Section of the Channel

Equation (6) is a method for calculating the travel distance of a debris flow in a channel with variable slope in a two-dimensional plane, which is applicable when the debris flow intersects the river at a right angle. However, debris flows often rush into rivers at a certain angle. Therefore, the movement calculation model is simplified to give a reasonable travel distance calculation for different spatial relationships.
As shown in Figure 6, after the debris flow enters the river, the movement track of the debris flow head is actually curved due to the fluid and bed forces. The complex changes in the head shape during curved movement will have a significant impact on both the resistance and dynamics equations. To simplify the process, the x- and y-directions are considered to be separate and without interaction, and the calculation of the y-direction travel distance still uses Equation (6), where the angle of the channel is calculated based on the spatial relationship in the plane, resulting in
cos θ = tan 2 β + 1 tan 2 β + sin 2 α + 1
where α is the longitudinal slope angle of the river, and β is the angle between the debris flow and the river channel.
The equation of the motion for the x-directional debris flow is
1 2 ρ 0 x h f a = 1 2 g h f 2 ρ 0 cos θ x + 1 2 c d u h f x g ( ρ σ ρ ) sin θ x f 2
where ρ 0 is the average density of the debris flow, θ x is the slope angle in the x-direction of the river channel, and f 2 is the internal resistance of the debris flow.
The angle of slope in the x-direction of the river channel can be calculated with the following equations:
sin θ x = tan β sin α 1 + tan 2 β ( 1 + sin 2 α )
cos θ x = 1 + tan 2 β 1 + tan 2 β ( 1 + sin 2 α )
The composition of debris flows is complex. Therefore, to establish a simple and practical resistance model, this study adopts structural two-phase resistance models to consider the resistance term in debris flows and mainly considers the frictional stress of solid-phase particles, which are independent of shear deformation, to obtain the velocity component in the river width direction as follows:
v l = ( 1 2 R t + Q ) sin β + a t cos ( π β )
The farthest distance is reached when the speed is 0, and the time is
t l = 2 Q sin β R sin β a cos ( π β )
Integrating over v l gives a maximum distance of
l max = ( 1 4 R t l 2 + Q t l ) sin β + 1 2 a t l 2 cos ( π β )

4.2.2. Deposit Morphology of Debris Flows

In Figure 7, the length of the debris flow in the x-direction is still unknown, but it can be seen that the shape and size in the x-direction can be calculated by determining the slope values of the front slope and the back slope, and finally the three-dimensional accumulation dimensions of the debris flow can be obtained.
As the debris flow enters the river, if the slope of the deposit surface is large, the surface flow can erode and further reshape the deposit slope, until forming a stable slope with a critical slope. Kuang derived the calculation equation of the critical slope values γ as follows [28]:
tan γ = C ( ρ σ ρ ) tan ϕ C ( ρ σ ρ ) + ρ [ 1 + 0.52 ( q J S W 2 g d m 3 ) 1 3 ]
where C is the volumetric concentration of debris flow deposits after the debris flows stop moving, d m is the mean particle size of the debris flow materials and q J S W is the single width flow rate of the surface flow over the deposit.

5. Calculation Model for Submerged-Type Debris Flow

In the analysis of the dam-type debris flow, the forces in the plane of the movement direction are mainly taken into account when calculating the transport distance and determining the degree of river blocking, but the upper water pressure on the debris flow is not considered, so it is not suitable for describing the submerged-type debris flow. Therefore, the dynamic movement characteristics of submerged-type debris flows will be analyzed in this section.

5.1. Underwater Movement Characteristics of Debris Flows

Underwater cohesive debris flows have an interface at their top during movement and deposition. Walker divided the deposit profile of underwater debris flows into four gradient processes [29]. (1) The uppermost part is a turbulent layer, in which deposit particles are disorganized and where multiple transport mechanisms (rolling, suspension, dragging, etc.) coexist. (2) When the debris flow travels a certain distance, the deposit begins to show a reverse progression. At this point, the mixing effect of the water reduces the viscosity of the upper layer of the debris flow body and concentrates the coarse particles toward the middle. (3) In the third state, water flow mixing continues to strengthen, and the turbulent flow ripples through the entire debris flow body. The larger particles settle, and the smaller particles cease to move as the flow velocity decreases, thus creating a positive progressive structure. (4) As the debris flow continues to travel and deposit, the deposit profile regressive layers become stable, and the traction of the upper part of the debris flow body increases and gradually becomes the dominant movement.
After entering the main river, the debris flow moves to the opposite bank, and the debris flow body deposits to form a submerged dam, which affects the flow movement, so whether the debris flow body can travel to the opposite bank is an important factor affecting this kind of river blocking mode. However, the underwater motion process of debris flows is complicated, and the external force forms and internal deposition modes are different at different stages. To facilitate theoretical analysis, it is necessary to simplify the process.

5.2. Movement Model of the Debris Flow Head

5.2.1. Equations of the Motion in the Two-Dimensional Plane

After the debris flow rushes into the main river, the front of the debris blow has a certain morphology. The end of the debris flow can be considered to form an interface with the water flow since the mixing of the debris flow by the water flow is ignored. Hence, the debris flow and the water flow can be analyzed as two separate fluids in motion and coupled according to their interaction characteristics.
The motion pattern of the debris flow head is illustrated in Figure 8, based on the model tests and field studies. Disregarding the material exchange between the debris flow and the riverbed, the motion equation of the mixed flow can be obtained as follows:
t ( ρ d u ) + y ( ρ d u 2 ) + z ( ρ d u w ) = ρ d g sin θ p y z τ z x
t ( ρ d w ) + y ( ρ d u w ) + z ( ρ d w 2 ) = ρ d g cos θ p z y τ z x
where p and τ z x are the static pressure and shear stress at the interface between the debris flow and the water flow, respectively.

5.2.2. Basic Equation of the Debris Flow Head

Substituting the momentum equations and mass equations into the model shown in Figure 8, and integrating the mass conservation equations for solids and mixed flows from the riverbed to the interface between the debris flow and water flow, while the volume concentration of solid particles can be considered to be close to zero as it is small at the interface, gives:
t 0 h C V d z + x 0 h u C V d z = w s C V b
t 0 h ρ d d z + x 0 h ρ d u d z = ρ d w e ρ s w s C V b
where u p , w p are the y- and z-directional flow rates at the debris flow and water flow interfaces, respectively; w b is the z-directional flow rate at the riverbed; and CVb and CVp are the volume concentrations of debris flow materials at the water flow interface and at the riverbed bottom, respectively.
For Equations (17) and (18), only the debris flow head is integrated in the y-direction, and assuming that the debris flow body does not change over time, y f t = u f , y e t = u e = v e ( v e is the flow rate of the debris flow body), Equations (17) and (18) can be simplified as
t y e y f 0 h C V d z d y = y e y f w s C V b d y
t y e y f 0 h ρ d d z d y = y e y f ρ d w e d y y e y f ρ σ w s C V b d y
t y e y f 0 h u ρ d d z d y = ρ g R sin θ x e x f 0 h C V d z d y 1 2 ρ g R cos θ C V h ¯ 2 0 h f p d z + 0 h e p d z x e x f z p d y + x e x f z b d y

5.3. Travel Distances of the Debris Flow Head

Since the measured data of submerged-type debris flows are limited, accurate characterization analysis on the debris flow head is difficult. Liu [27] studied the characteristics of submerged flows using flume tests and proposed a simplified method for the basic equations of debris flow heads based on the test results.
By ignoring the head speed, Equations (19) and (20) can be simplified as
d d t ( p 0 A ¯ ) = 0
d d t ( c A ¯ ) = 0
Simplifying the pressure and friction terms in Equation (21), a simplified equation can be obtained:
d d t ( ρ d U A ¯ ) = ρ g R sin θ C A ¯ + ρ g R cos θ C v h ¯ 2 / 2 ρ C D U f 2 A z f p ρ U 2 P p / 2 f b ρ U 2 P b / 2
where f b , f p are the coefficients of friction at the intersection and riverbed, respectively; and P P , P b are the wetted lengths at the intersection and the riverbed, respectively.
Assuming that the average head height is constant and considering that there is a functional relationship between the head geometry and head height yields
L f = ξ f h max , L e = ξ e h max A ¯ = ξ A ¯ h max 2 , A z = ξ A z h max 2 P p = ξ p h max , P b = ξ b h max
After substitution of these relations into Equations (20)–(24), these equations can be simplified as
d d t ( ρ d h max 2 ) = 0 , d d t ( C h max 2 ) = 0
ξ A ¯ d d t ( ρ d U h max 2 ) = ρ g R sin θ C ξ A ¯ h max 2 + 1 2 ρ g R cos θ C v h max 2 ρ C D ( U + ξ f d h max d t ) 2 ξ A z h max 2 1 2 f p ρ U 2 ξ p h max 1 2 f b ρ U 2 ξ b h max
The changes in density and volume concentration of the mixed flow are not considered when the velocity of the river water is low, and the analytical solution to Equations (25) and (26) is
U = a ( a + c U 0 ) e 2 a c t ( a c U 0 ) c ( a c U 0 ) + ( a + c U 0 ) e 2 a c t
where
a = [ g R sin θ C V 0 ξ A ¯ h 0 2 + 1 2 g R cos θ C V 0 h e 2 ] / [ ξ A ¯ ( 1 + R C V 0 h 0 2 ) ] c = [ C D ξ A ¯ h 0 2 + 1 2 f p ξ p h 0 + 1 2 f b ξ b h 0 ] / [ ξ A ¯ ( 1 + R C V 0 h 0 2 ) ]
The velocity sequence { U i } of the debris flow is calculated at interval Δ T . By using Equation (27), the maximum travel distance of the debris flow along the river after the debris flow entered the river can be calculated as
l max = i = 1 n U i Δ T i sin β

6. Identification Method of River Blocking

The motion equations and maximum travel distance for dam-type and submerged-type debris flows are given in this study. A comprehensive method to determine the degree of river blocking by a debris flow is proposed herein, taking into account two factors: the length of the debris flow deposit and the amount of material.
(1) Maximum blocking degree
According to Equation (28), the maximum travel distance lmax can be calculated, and then, the maximum blocking degree S can be defined as
S = l max L width   of   main   river
where S [ 0 , 1 ] .
Although this index can reflect the degree of river blocking by debris flow, it is still insufficient to represent the actual river blocking capacity because only the maximum travel distance is discussed, and the water interception capacity of the section is not considered. Therefore, the material conditions needed for river blocking should be considered in order to establish a more reasonable criterion.
(2) Source volume
Zhou et al. [30] proposed a formula for calculating the total amount of debris flow deposition in the main river:
V b = ( 1 2 tan 14 + 1 2 tan ϕ ) H 2 L
This equation is applicable when the river is completely blocked, but it is no longer applicable when a local blockage actually occurs in practice. This study improves the formula based on the model illustrated in Figure 9 by replacing the critical slope value calculated in Equation (14) with the 14° assumed in Equation (30), and assuming the minimum height of the debris flow head is hf, the lateral accumulation slope of the debris flow along the river is determined as
tan ψ = H h f l f
By establishing a local coordinate axis, the calculation formula of the height of different parts of the debris flow head can be obtained as
h ( x ) = H x tan ψ
where x ( 0 , l f ) .
Equation (28) calculates the transport distance after the debris flow enters the river. In this study, the transport distance after the debris flow entered the river (lmax) is classified according to the correlation between the transport distance after the debris flow entered the river and the opposite bank, which is divided into three categories in total, as described below.
Equation (28) is used to calculate the travel distance of the debris flow after merging with the river. In this study, the travel distance (lmax) of debris flows after merging with rivers is classified according to the relative relationship between the debris flow and the opposite bank, which is divided into three categories as follows:
First, l max ( l f , L ) , if the debris flow has not reached the opposite bank, the volume of its locally blocked dam can be calculated by the following equation:
V b 1 = ( 1 2 tan γ + 1 2 tan ϕ ) H 2 ( l max l f ) + ( 1 2 tan γ + 1 2 tan ϕ ) 0 l f h 2 ( x ) d x = ( 1 2 tan γ + 1 2 tan ϕ ) ( H 2 ( 2 l max l f ) l f 2 H tan ψ + l f 3 3 tan 2 ψ )
Second, l max [ L , L + l f ) , if the front of the debris flow head reaches the opposite bank but has not yet formed a complete river-blocking dam, the volume of the dam can be calculated by the following formula:
V b 2 = ( 1 2 tan γ + 1 2 tan ϕ ) H 2 ( l max l f ) + ( 1 2 tan γ + 1 2 tan ϕ ) 0 L l max + l f h 2 ( x ) d x = ( 1 2 tan γ + 1 2 tan ϕ ) ( H 2 L ( L l max + l f ) 2 H tan ψ + ( L l max + l f ) 3 3 tan 2 ψ )
Third, if l max L + l f , the debris flow can be considered to be blocking the river, and the volume of the dam is the same as in Equation (31), replacing only the critical slope:
V b 3 = ( 1 2 tan γ + 1 2 tan ϕ ) H 2 L
As the debris flow is scoured and mixed by the water flow as it enters the main river, some of the solid materials in the debris flow are carried away by the water flow. Therefore, the total amount of actual materials required to form a debris flow dam needs to be greater than the amount of debris flow that is deposited in the main river. The definition of the loss coefficient is therefore defined as follows:
r = V w V w + V b
where V w is the volume of solids washed away by the river and V b is the volume of debris flow solids forming the dam.
From the above equation, it can be seen that the determination of the loss coefficient is of great significance to the analysis of the material amount required for debris flow blockage. Based on the field model tests in the Dadu River basin, a correlation between measured loss coefficients and river blocking coefficients is obtained, as shown in Figure 10.
As presented in Figure 10, there is a significant negative correlation between the loss coefficient and the blocking degree. The loss coefficient is difficult to obtain in practice, and the maximum blocking degree coefficient obtained from Equation (29) can be used to estimate the loss coefficient by the correlation presented in Figure 10.
The amount of materials required to form the dam can be calculated using Equations (31)–(33) and combined with the loss coefficient to give the total amount of material:
V n = V b n 1 r
where n is the classification of debris flow travel distance intervals; V b 1 , V b 2 , V b 3 are the amount of debris flow materials that are deposited in the main river under conditions where the head of the debris flow does not reach the opposite bank, reaches the opposite bank but does not blocked, and reaches the opposite bank and is blocked, respectively; and V1, V2, and V3 are the total amount of source materials required to form a debris flow dam under conditions where the debris flow head does not reach the opposite bank, reaches the opposite bank but does not block the river, and reaches the opposite bank and blocks the river, respectively.
On the basis of the maximum blocking degree coefficient and the amount of material sources required, the identification criterion of the river blocking by debris flow is proposed as follows:
When S = 1, V > V3, completely blocked;
When S = 1, V > V2, mostly blocked;
When S < 1, V > V1, half blocked.

7. Discussion

Currently, the assessment of debris flow hazard intensity and extent primarily rely on empirical formulas, statistical analogies, and numerical analyses based on dynamic processes. While empirical formulas and statistical analogies can estimate the velocity and travel distance of debris flows, they often exhibit significant errors. For the movement distance of a single debris flow event, the error in empirical formulas may exceed 50% (Guo et al., 2014). Thus, we established the method for the identification of river blocking by debris flow through the modeling test carried out in the Dadu River. This method has an advantage over other methods in that it is easier to compute, provides a quicker assessment of the debris flow, and expedites the provision of emergency relief. Nevertheless, the loss coefficient of this method is calculated through model experiments, and it is necessary to supplement the debris flow data in the Dadu River area to improve its accuracy.
There are many softwares, such as Mass flow2.5, flow 3D 2022R1 and ANSYS CFX2021, which can be used to simulate debris flow blocking rivers, among which the Massflow software demonstrates high computational efficiency, supports secondary development, and has been applied in the numerical analysis of real geological hazards [31]. Subsequently, we will employ the comprehensive analysis method via the software and examine the case of debris flow in the Dadu River in order to verify its precision and suitability.

8. Conclusions

In this study, an identification method for river blocking by debris flows in the middle reaches of the Dadu River is proposed, and the following conclusions can be drawn:
  • The calculation formulas for the maximum travel distance of the two kinds of debris flows entering the river are obtained through theoretical derivation.
  • The formulas for calculating the length and volume of debris flow accumulation are derived, and the relationship between the debris flow loss coefficient and river blocking degree in the middle part of the Dadu River is analyzed.
  • Based on the relationship between the maximum blocking degree coefficient and the amount of material sources needed, the identification criterion of river blocking by debris flow is proposed.
  • This identification method can swiftly identify a complete dam blockage in a river, but more examples are needed to adjust its loss coefficient.

Author Contributions

Conceptualization, G.F. and Z.S.; methodology, Z.S. and Y.C.; software, Y.C.; validation, G.F.; formal analysis and investigation, G.F. and Z.S.; resources, G.F.; writing—original draft preparation, Y.C. and D.L.; writing—review and editing, G.F. and D.L.; visualization, Y.C.; supervision, Z.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (U20A20111).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, H.; Bai, X.; Zhai, X.; Zhao, J.; Zhu, X.; Li, C.; Liu, K.; Wang, Q. Research on extension evaluation method of mudslide hazard based on analytic hierarchy Process–Criteria importance through intercriteria correlation combination assignment of game theory ideas. Water 2023, 15, 2961. [Google Scholar] [CrossRef]
  2. Pei, Y.; Qiu, H.; Yang, D.; Liu, Z.; Ma, S.; Li, J.; Cao, M.; Wufuer, W. Increasing landslide activity in the Taxkorgan River Basin (eastern Pamirs Plateau, China) driven by climate change. Catena 2023, 223, 106911. [Google Scholar] [CrossRef]
  3. Pei, Y.; Qiu, H.; Zhu, Y.; Wang, J.; Yang, D.; Tang, B.; Wang, F.; Cao, M. Elevation dependence of landslide activity induced by climate change in the eastern Pamirs. Landslides 2023, 20, 1115–1133. [Google Scholar] [CrossRef]
  4. Yang, D.; Qiu, H.; Ye, B.; Liu, Y.; Zhang, J.; Zhu, Y. Distribution and Recurrence of Warming-Induced Retrogressive Thaw Slumps on the Central Qinghai-Tibet Plateau. JGR Earth Surf. 2023, 128, e2022JF007047. [Google Scholar] [CrossRef]
  5. Cluff, L.S. Peru earthquake of May 31, 1970; Engineering geology observations. Bull. Seismol. Soc. Am. 1971, 61, 511–533. [Google Scholar] [CrossRef]
  6. Lyon, B. Enhanced Seasonal Rainfall in Northern Venezuela and the Extreme Events of December 1999. J. Clim. 2003, 16, 2302–2306. [Google Scholar] [CrossRef]
  7. Yang, F.; Fan, X.; Siva Subramanian, S.; Dou, X.; Xiong, J.; Xia, B.; Yu, Z.; Xu, Q. Catastrophic debris flows triggered by the 20 August 2019 rainfall, a decade since the Wenchuan earthquake, China. Landslides 2021, 18, 3197–3212. [Google Scholar] [CrossRef]
  8. Li, Y.; Liu, X.N.; Gan, B.R.; Wang, X.K.; Yang, X.G.; Li, H.B.; Long, X.Y.; Zhou, J.W. Formation-Evolutionary Mechanism Analysis and Impacts of Human Activities on the 20 August 2019 Clustered Debris Flows Event in Wenchuan County, Southwestern China. Front. Earth Sci. 2021, 9, 616113. [Google Scholar] [CrossRef]
  9. Tang, C.; Yang, Y.; Su, Y.; Ding, J.; Huang, W. The disastrous 23 July 2009 debris flow in Xiangshui Gully, Kangding County, Southwestern China. J. Mt. Sci. 2011, 8, 131–139. [Google Scholar] [CrossRef]
  10. Xie, H.; Liu, W.M.; Zhao, J.H.; Hu, K.H. Characteristics of Tangjiagou Debris Flow in Shimian of Sichuan in July 14, 2012. J. Earth Sci. Environ. 2013, 35, 90–97. [Google Scholar]
  11. Song, Z.; Deng, R.G.; Chen, Z.S. Analysis of motion process and accumulation zone characteristics of a typical debris flow blocking river case: A case study of the debris flow in Xiongjia Gully. J. Nat. Disasters 2016, 25, 74–80. [Google Scholar]
  12. Best, L.J. Sediment transport and bed morphology at river channel confluences. Sedimentology 1988, 35, 481–498. [Google Scholar] [CrossRef]
  13. Leite Ribeiro, M.; Blanckaert, K.; Roy, A.G.; Schleiss, A.J. Flow and sediment dynamics in channel confluences. J. Geophys. Res. Solid Earth 2012, 117, F01035. [Google Scholar] [CrossRef]
  14. Zheng, H.; Shi, Z.; Hanley, K.J.; Peng, M.; Guan, S.; Feng, S.; Chen, K. Deposition characteristics of debris flows in a lateral flume considering upstream entrainment. Geomorphology 2021, 394, 107960. [Google Scholar] [CrossRef]
  15. Wang, X.; Yan, X.; Duan, H.; Liu, X.; Huang, E. Experimental study on the influence of river flow confluences on the open channel stage–discharge relationship. Hydrol. Sci. J. 2019, 64, 2025–2039. [Google Scholar] [CrossRef]
  16. Yu, B.; Yang, C.Y.; Yu, M. Experimental study on the critical condition of river blockage by a viscous debris flow. Catena 2022, 213, 106198. [Google Scholar] [CrossRef]
  17. Zhu, P.Y.; Cheng, Z.L.; You, Y. Research on causes of river blocking by sediment delivery of Peilonggou Gully debris now in the Sichuan-Xizang Highway. J. Nat. Disasters 2000, 9, 80–83. [Google Scholar]
  18. Li, H.; Qi, S.; Chen, H.; Liao, H.; Cui, Y.; Zhou, J. Mass movement and formation process analysis of the two sequential landslide dam events in Jinsha River, Southwest China. Landslides 2019, 16, 2247–2258. [Google Scholar] [CrossRef]
  19. Huang, H.; Xie, Z.S.; Shi, S.W.; Yu, T. Characteristics and Countermeasures of River-blocking Debris Flow of Haermu Gully in Wenchuan Area After Earthquake. Bull. Soil Water Conserv. 2015, 35, 327–337. [Google Scholar]
  20. Zhang, S.; Wang, F.W.; Li, R. First insight into the catastrophic Atami debris fow induced by a rain gush on 3 July 2021 in Shizuoka, Japan. Landslides 2022, 19, 527–532. [Google Scholar] [CrossRef]
  21. Gong, X.L.; Chen, X.Q.; Chen, J.G.; Song, D.R. Effects of material composition on deposition characteristics of runoff-generated debris flows. Landslides 2023, 20, 2603–2618. [Google Scholar] [CrossRef]
  22. Ikeya, H. Debris flow and its countermeasures in Japan. Bull. Eng. Geol. Environ. 1989, 40, 15–33. [Google Scholar] [CrossRef]
  23. Prochaska, A.B.; Santi, P.M.; Higgins, J.D.; Cannon, S.H. Debris-flow runout predictions based on the average channel slope (ACS). Eng. Geol. 2008, 98, 29–40. [Google Scholar] [CrossRef]
  24. Fang, M.; Qi, X. Simulation and Prediction Algorithm for the Whole Process of Debris Flow Based on Multiple Data Integration. Water 2023, 15, 2778. [Google Scholar] [CrossRef]
  25. Wang, F.; Wang, J.D.; Chen, X.Q.; Chen, J.G. The influence of temporal and spatial variations on phase separation in debris flow deposition. Landslides 2019, 16, 497–514. [Google Scholar] [CrossRef]
  26. Takahashi, T.; Yoshida, K. Study on the deposition of debris flow, part l-deposition due to abrupt change of bed slope. Annu. Disaster Prev. Res. Inst. Kyoto Univ. 1997, 22, 315–328. [Google Scholar]
  27. Liu, C.R. Critera of Blocking Large River by Debris Flow and Disaster Reduction Countermeasures in Highways along Rivers; Southwest Jiaotong University: Chengdu, China, 2014; pp. 56–66. [Google Scholar]
  28. Kuang, S.F. Study on behaviors and deposit processes of debris flow at the confluence. J. Sediment Res. 1995, 1, 1–15. [Google Scholar]
  29. Walker, R.G. Turbidities and associated coarse clastic deposits. Facies Models 1976, 3, 25–36. [Google Scholar]
  30. Zhou, B.F.; Li, D.J.; Lv, R.R.; Luo, D.F.; Yang, Q.X. Direction for the Mitigations of Debris Flow; Science Press: Beijing, China, 1993. [Google Scholar]
  31. Ouyang, C.; He, S.; Tang, C. Numerical analysis of dynamics of debris flow over erodible beds in Wenchuan earthquake-induced area. Eng. Geol. 2015, 194, 62–72. [Google Scholar] [CrossRef]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Water 15 04301 g001
Figure 2. Illustration of the study area: (a) distribution of the debris flows; and (b) debris flow density in the study area.
Figure 2. Illustration of the study area: (a) distribution of the debris flows; and (b) debris flow density in the study area.
Water 15 04301 g002
Figure 3. Two typical motion types of debris flows in the middle reaches of the Dadu River after entering into the main river channel.
Figure 3. Two typical motion types of debris flows in the middle reaches of the Dadu River after entering into the main river channel.
Water 15 04301 g003aWater 15 04301 g003b
Figure 4. Diagram of the motion of the debris flow.
Figure 4. Diagram of the motion of the debris flow.
Water 15 04301 g004
Figure 5. Diagram of the debris flow in plane in a channel with a variable slope.
Figure 5. Diagram of the debris flow in plane in a channel with a variable slope.
Water 15 04301 g005
Figure 6. Illustration of the spatial relationship between the debris flow and the river.
Figure 6. Illustration of the spatial relationship between the debris flow and the river.
Water 15 04301 g006
Figure 7. Diagram of the x-directional forces.
Figure 7. Diagram of the x-directional forces.
Water 15 04301 g007
Figure 8. Illustration of the motion of the debris flow head.
Figure 8. Illustration of the motion of the debris flow head.
Water 15 04301 g008
Figure 9. Calculation of debris flow deposits.
Figure 9. Calculation of debris flow deposits.
Water 15 04301 g009
Figure 10. Correlation between the loss coefficient and the blocking degree based on the field model tests in the Dadu River basin.
Figure 10. Correlation between the loss coefficient and the blocking degree based on the field model tests in the Dadu River basin.
Water 15 04301 g010
Table 1. Values for the drag coefficients CD.
Table 1. Values for the drag coefficients CD.
Conditions of the Debris Flow Entering the RiverCD
Gentle riverbed, low velocity and low turbulence0.5–1.0
Sloping riverbed, medium velocity and medium turbulence1.0–1.5
Steeply sloping areas, high velocity and turbulence1.5–2.5
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Song, Z.; Fan, G.; Chen, Y.; Liu, D. Identification Method of River Blocking by Debris Flow in the Middle Reaches of the Dadu River, Southwest of China. Water 2023, 15, 4301. https://doi.org/10.3390/w15244301

AMA Style

Song Z, Fan G, Chen Y, Liu D. Identification Method of River Blocking by Debris Flow in the Middle Reaches of the Dadu River, Southwest of China. Water. 2023; 15(24):4301. https://doi.org/10.3390/w15244301

Chicago/Turabian Style

Song, Zhi, Gang Fan, Yanni Chen, and Darui Liu. 2023. "Identification Method of River Blocking by Debris Flow in the Middle Reaches of the Dadu River, Southwest of China" Water 15, no. 24: 4301. https://doi.org/10.3390/w15244301

APA Style

Song, Z., Fan, G., Chen, Y., & Liu, D. (2023). Identification Method of River Blocking by Debris Flow in the Middle Reaches of the Dadu River, Southwest of China. Water, 15(24), 4301. https://doi.org/10.3390/w15244301

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop