Two-dimensional Dynamics Simulation of Two-phase Debris Flow

To investigate the movement mechanism of debris flow, a two-dimensional, two-phase, depthintegrated model is introduced. The model uses Mohr-Coulomb plasticity for the solid rheology, and the fluid stress is modeled as a Newtonian fluid. The interaction between solid and liquid phases, which plays a major role in debris flow movement, is assumed to consist of drag and buoyancy forces. The applicability of drag force formulas is discussed. Considering the complex interaction between debris flow and the bed surface, a combined friction boundary condition is imposed on the bottom, and this is also discussed. To solve the complex model equations, a numerical method with second-order accuracy based on the finite volume method is proposed. Several numerical experiments are performed to verify the feasibilities of model and numerical schemes. Numerical results demonstrate that different solid volume fractions substantially affect debris flow movement.

Coupled Model of Two-phase Debris Flow,Sediment Transport and Morphological Evolution

The volume fraction of the solid and liquid phase of debris flows, which evolves simultaneously across terrains, largely determines the dynamic property of debris flows. The entrainment process significantly influences the amplitude of the volume fraction. In this paper, we present a depth-averaged two-phase debris-flow model describing the simultaneous evolution of the phase velocity and depth, the solid and fluid volume fractions and the bed morphological evolution. The model employs the Mohr–Coulomb plasticity for the solid stress, and the fluid stress is modeled as a Newtonian viscous stress. The interfacial momentum transfer includes viscous drag and buoyancy. A new extended entrainment rate formula that satisfies the boundary momentum jump condition （Iverson and Ouyang, 2015） is presented. In this formula, the basal traction stress is a function of the solid volume fraction and can take advantage of both the Coulomb and velocity-dependent friction models. A finite volume method using Roe’s Riemann approximation is suggested to solve the equations. Three computational cases are conducted and compared with experiments or previous results. The results show that the current computational model and framework are robust and suitable for capturing the characteristics of debris flows.

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.

Protective effects of baffles with different positions,row spacings,heights on debris flow impact

The baffle effectively slowed down debris flow velocity,reduced its kinetic energy,and significantly shortened the distance of debris flow movement.Consequently,they are widely used for protection against natural hazards such as landslides and mudslides.This study,based on the threedimensional DEM(Discrete Element Method),investigated the impact of different baffle positions on debris flow protection.Debris flow velocity and kinetic energy variations were studied through single-factor experiments.Suitable baffle positions were preliminarily selected by analyzing the influence of the first-row baffle position on the impact force and accumulation mass of debris flow.Subsequently,based on the selected baffle positions and four factors influencing the effectiveness of baffle protection(baffle position(P),baffle height(h),row spacing(S_(r)),and angle of transit area(α)),an orthogonal design was employed to further explore the optimal arrangement of baffles.The research results indicate that the use of a baffle structure could effectively slow down the motion velocity of debris flows and dissipate their energy.When the baffle is placed in the transit area,the impact force on the first-row baffle is greater than that when the baffle is placed in the deposition area.Similarly,when the baffle is placed in the transit area,the obstruction effect on debris flow mass is also greater than that when the baffle is placed in the deposition area.Through orthogonal experimental range analysis,when the impact on the first row of baffles is used as the evaluation criterion,the importance of each influencing factor is ranked asα>P>S_(r)>h.When the mass of debris flow behind the baffle is regarded as the evaluation criterion,the rank is changed to P>α>S_(r)>h.The experimental simulation results show that the optimal baffle arrangement is:P_(5),S_(r)=16,α=35°,h=9.

Debris flow runout behaviors considering the influences of densely populated buildings

Debris flows pose serious risks to communities in mountainous areas,often resulting in large losses of human life and property.The impeding presence of urban buildings often affects the runout behavior and deposition of debris flows.But the impact of different building densities and sizes on debris flow dynamics has yet to be quantified to guide urban planning in debris flow risk zones.This study focused on a debris flow that occurred in Zhouqu County,Gansu Province,China on August 7th,2010,which was catastrophic and destroyed many buildings.The FLO-2D software was used to simulate this debris flow in two scenarios,i.e.the presence and the absence of buildings,to obtain debris-flow intensity parameters.The developed model was then used to further analyze the influence of large buildings and narrow channels within the urban environment.The simulation results show that considering the presence of buildings in the simulation is essential for accurate assessment of debris flow intensity and deposition distribution.The layout of buildings in the upstream urban area,such as large buildings or parallel buildings which form narrow channels,can affect the flow velocity and depth of debris flow heading towards downstream buildings.To mitigate damage to downstream buildings,the relative spacing(d/a)between upstream and downstream buildings should not exceed a value of two and should ideally be even lower.These findings provide valuable insights for improving the resistance of mountainous cities to urban debris flows.

Evaluation of frictional pressure drop correlations for air-water and air-oil two-phase flow in pipeline-riser system

Accurate prediction of the frictional pressure drop is important for the design and operation of subsea oil and gas transporting system considering the length of the pipeline. The applicability of the correlations to pipeline-riser flow needs evaluation since the flow condition in pipeline-riser is quite different from the original data where they were derived from. In the present study, a comprehensive evaluation of 24prevailing correlation in predicting frictional pressure drop is carried out based on experimentally measured data of air-water and air-oil two-phase flows in pipeline-riser. Experiments are performed in a system having different configuration of pipeline-riser with the inclination of the downcomer varied from-2°to-5°to investigated the effect of the elbow on the frictional pressure drop in the riser. The inlet gas velocity ranges from 0.03 to 6.2 m/s, and liquid velocity varies from 0.02 to 1.3 m/s. A total of885 experimental data points including 782 on air-water flows and 103 on air-oil flows are obtained and used to access the prediction ability of the correlations. Comparison of the predicted results with the measured data indicate that a majority of the investigated correlations under-predict the pressure drop on severe slugging. The result of this study highlights the requirement of new method considering the effect of pipe layout on the frictional pressure drop.

Debris Fan Produced by Failure of Canyon-Blocking Pyroclastic Flows

Ash-rich pyroclastic flows from the cataclysmic eruption of Mount Mazama (~7700 yr. B. P.), Cascade volcanic arc, Oregon, entered and blocked the narrow, bedrock-lined canyon of the Williamson River approximately 35 to 44 km from the source volcano. The blockage impounded a body of water which then released producing four stratigraphic units in the downstream debris fan. The four stratigraphic units are a boulder core comprised of locally sourced bedrock boulders and three sand-rich units including a fine-grained sand unit, a sandy pumice gravel (±basalt/hydrovolcanic tuff) unit, and a pumice pebble-bearing, crystal-rich sand unit. Hand-drilled auger holes up to ~1.6 m deep were used to obtain samples of the sand-rich units. Units were delimited using surface and down-hole observations, composition and texture, estimated density, statistical parameters of grain size, and vertical and lateral distribution of properties. Overtopping followed by rapid incision into the ash-rich pyroclastic flows progressively cleared the canyon, but a bedrock knickpoint near the head of the canyon limited the volume of debris available for transport to about 0.04 km3 to 0.08 km3. Co-deposition of bedrock boulders and lithic-rich sand was followed by rapid deposition with minimal reworking of remobilized pyroclastics. Continued draining of the impounded lake sent hyperconcentrated flows onto the debris fan depositing pumice-rich gravels that graded upward to crystal-rich sands.

Dynamic-based model for calculating the boulder impact force in debris flow

The boulder impact force in debris flow is generally calculated by static methods such as the cantilever beam models.However,these methods cannot describe the dynamic scenario of boulder collision on structures,so the inertia and damping effects of the structures are not involved causing an overestimation on the boulder impact force.In order to address this issue,a dynamic-based model for calculating the boulder impact force of a debris flow was proposed in this study,and the dynamic characteristics of a cantilever beam with multiple degrees of freedom under boulder collision were investigated.By using the drop-weight method to simulate boulders within debris flow,seven experiments of drop-weight impacting the cantilever beam were used to calibrate the error of the dynamicbased model.Results indicate that the dynamic-based model is able to reconstruct the impact force history on the cantilever beam during impact time and the error of dynamic-based model is 15.3%in calculating boulder impact force,significantly outperforming the cantilever beam model’s error of 285%.Therefore,the dynamic-based model can overcome the drawbacks of the static-based models and provide a more reliable theoretical foundation for the engineering design of debris flow control structures.

Discharge evolution law of debris flow based on a sharp bend physical modeling test

For the basins with debris flow development,its channel terrain exhibits a tortuous shape,which significantly restricts the movement of debris flows and leads to severe erosion effects on the concave bank.Therefore,this study aims to analyze the shear force of debris flows within the bend channel.We established the relationship between the shear force and bend curvature through laboratory experiments.Under the long-term erosion by debris flows,the curvature radius of bends gradually increases,however,when this increasing trend reaches an equilibrium state with the intensity of debris flow discharge,there will be no significant change in curvature radius.In general,the activity pattern and discharges of debris flows would remain relatively stable.Hence,we can infer the magnitude of debris flow discharges from the terrain parameters of the bend channel.

Socio-scientific quantification of the comprehensive benefits of debris flow mitigation measures for villages in western Sichuan, China

Debris flow hazards seriously threaten thesafety and sustainable development of mountainousareas. Numerous debris flow mitigation measures havebeen implemented worldwide;however, acomprehensive assessment of the specific disasterreduction effects of these measures and their economic,social and ecological benefits is yet to be performed.The western region of Sichuan Province frequentlysuffers from geohazards such as debris flow, and thegovernment has adopted many mitigation measures.This study assessed the benefits of debris flowmitigation measures and identified the key influencingfactors via a field-based study conducted in 81 villagesin western Sichuan province, China. A framework forthe evaluation of the benefits of rural debris flowmitigation measures was constructed andquantitatively evaluated using a survey. Snowballsampling was performed to recruit 81 village leadersand 468 farmers. The results showed that managementand engineering measures were the main methodsused to mitigate debris flow;ecological measures wereauxiliary. The average satisfaction scores of farmers forthese three types of measures were 4.07, 3.90, and 3.56,respectively (as measured on a five-point Likert scale).In contrast, in terms of the benefits of these mitigationmeasures, only a small proportion of villages (11.11%)obtained a high level of comprehensive benefits fromthe debris flow mitigation measures, while the majority(88.89%) received medium to low-level benefits. Toimprove this situation, we further studied and foundthat the main factors that restricted villages fromachieving high-level comprehensive benefits were theunpredictable nature of debris flows, labour forceoutflow and remoteness. Effective control measures, agood economic environment and strong governmentassistance were reported as crucial factors forimproving these comprehensive benefits. This studyprovides socio-scientific references for decisionmakingon rural debris flow mitigation measures while keeping villages at the centre of economic development.

Interaction Mechanisms between Natural Debris Flow and Rigid Barrier Deflectors:A New Perspective for Rational Design and Optimal Arrangement

Rigid barrier deflectors can effectively prevent overspilling landslides,and can satisfy disaster prevention requirements.However,the mechanisms of interaction between natural granular flow and rigid barrier deflectors require further investigation.To date,few studies have investigated the impact of deflectors on controlling viscous debris flows for geological disaster prevention.To investigate the effect of rigid barrier deflectors on impact mechanisms,a numerical model using the smoothed particle hydrodynamics(SPH)method with the Herschel–Bulkley model is proposed to simulate the interaction between natural viscous flow and single/dual barriers with and without deflectors.This model was validated using laboratory flume test data from the literature.Then,the model was used to investigate the influence of the deflector angle and multi-barrier arrangements.The optimal configuration of multi-barriers was analyzed with consideration to the barrier height and distance between the barriers,because these metrics have a significant impact on the viscous flow pile-up,run-up,and overflow mechanisms.The investigation considered the energy dissipation process,retention efficiency,and dead-zone formation.Compared with bare barriers with similar geometric characteristics and spatial distribution,rigid barriers with deflectors exhibit superior effectiveness in preventing the overflow and overspilling of viscous debris flow.Recommendations for the rational design of deflectors and the optimal arrangement of multi-barriers are provided to mitigate geological disasters.

Dynamic response of buildings under debris flow impact

This study employs the smoothed particle hydrodynamics–finite element method(SPH–FEM) coupling numerical method to investigate the impact of debris flow on reinforced concrete(RC)-frame buildings. The methodology considers the variables of debris flow depth and velocity and introduces the intensity index IDV(IDV = DV) to evaluate three different levels of debris flow impact intensity. The primary focus of this study is to investigate the dynamic response and failure mechanism of RC-frame buildings under debris flow impact, including structural failure patterns, impact force and column displacement. The results show that under a highintensity impact, a gradual collapse process of the RCframe building can be observed, and the damage mode of the frame column reflects shear failure or plastic hinge failure mechanism. First, the longitudinal infill walls are damaged owing to their low out-of-plane flexural capacity;the critical failure intensity index IDV value is approximately 7.5 m2/s. The structure cannot withstand debris flows with an intensity index IDV greater than 16 m2/s, and it is recommended that the peak impact force should not exceed 2100 k N. The impact damage ability of debris flow on buildings mostly originates from the impact force of the frontal debris flow, with the impact force of the debris flow body being approximately 42% lower than that of the debris flow head. Finally, a five-level classification system for evaluating the damage status of buildings is proposed based on the numerical simulation and investigation results of the disaster site.

Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Modelling of debris-flow susceptibility and propagation: a case study from Northwest Himalaya

The geological and geographical position of the Northwest Himalayas makes it a vulnerable area for mass movements particularly landslides and debris flows. Mass movements have had a substantial impact on the study area which is extending along Karakorum Highway(KKH) from Besham to Chilas. Intense seismicity, deep gorges, steep terrain and extreme climatic events trigger multiple mountain hazards along the KKH, among which debris flow is recognized as the most destructive geohazard. This study aims to prepare a field-based debris flow inventory map at a regional scale along a 200 km stretch from Besham to Chilas. A total of 117 debris flows were identified in the field, and subsequently, a point-based debris-flow inventory and catchment delineation were performed through Arc GIS analysis. Regional scale debris flow susceptibility and propagation maps were prepared using Weighted Overlay Method(WOM) and Flow-R technique sequentially. Predisposing factors include slope, slope aspect, elevation, Topographic Roughness Index(TRI), Topographic Wetness Index(TWI), stream buffer, distance to faults, lithology rainfall, curvature, and collapsed material layer. The dataset was randomly divided into training data(75%) and validation data(25%). Results were validated through the Receiver Operator Characteristics(ROC) curve. Results show that Area Under the Curve(AUC) using WOM model is 79.2%. Flow-R propagation of debris flow shows that the 13.15%, 22.94%, and 63.91% areas are very high, high, and low susceptible to debris flow respectively. The propagation predicated by Flow-R validates the naturally occurring debris flow propagation as observed in the field surveys. The output of this research will provide valuable input to the decision makers for the site selection, designing of the prevention system, and for the protection of current infrastructure.

Regulation effect of the grille spacing of a funnel-type grating water–sediment separation structure on the debris flow performance

The size of pores or the grille spacing of water–sediment separation structures directly affects their regulation effect on the debris flow performance.A suitable pore size or grille spacing can effectively improve the water–sediment separation ability of the structure.The new funnel-type grating water–sediment separation structure(FGWSS)combines vertical and horizontal structures and provides a satisfactory water–sediment separation effect.However,the regulation effect of the grille spacing of the structure on the debris flow performance has not been studied.The regulation effect of the structure grille spacing on the debris flow performance is studied through a flume test,and the optimal structure grille spacing is obtained.An empirical equation of the relationship between the relative grille spacing of the structure and the sediment separation rate is established.Finally,the influence of the water–sediment separation structure on the regulation effect of debris flows is examined from two aspects:external factors(properties of debris flows)and internal factors(structural factors).The experimental results show that the gradation characteristics of solid particles in debris flows constitute a key factor affecting the regulation effect of the structure on the debris flow performance.The optimum grille spacing of the FGWSS matches the particle size corresponding to the material distribution curves d85~d90 of the debris flow.The total separation rate of debris flow particles is related to the grille spacing of the structure and the content of coarse and fine particles in the debris flow.

Comparison of debris flow susceptibility assessment methods:support vector machine,particle swarm optimization,and feature selection techniques

The selection of important factors in machine learning-based susceptibility assessments is crucial to obtain reliable susceptibility results.In this study,metaheuristic optimization and feature selection techniques were applied to identify the most important input parameters for mapping debris flow susceptibility in the southern mountain area of Chengde City in Hebei Province,China,by using machine learning algorithms.In total,133 historical debris flow records and 16 related factors were selected.The support vector machine(SVM)was first used as the base classifier,and then a hybrid model was introduced by a two-step process.First,the particle swarm optimization(PSO)algorithm was employed to select the SVM model hyperparameters.Second,two feature selection algorithms,namely principal component analysis(PCA)and PSO,were integrated into the PSO-based SVM model,which generated the PCA-PSO-SVM and FS-PSO-SVM models,respectively.Three statistical metrics(accuracy,recall,and specificity)and the area under the receiver operating characteristic curve(AUC)were employed to evaluate and validate the performance of the models.The results indicated that the feature selection-based models exhibited the best performance,followed by the PSO-based SVM and SVM models.Moreover,the performance of the FS-PSO-SVM model was better than that of the PCA-PSO-SVM model,showing the highest AUC,accuracy,recall,and specificity values in both the training and testing processes.It was found that the selection of optimal features is crucial to improving the reliability of debris flow susceptibility assessment results.Moreover,the PSO algorithm was found to be not only an effective tool for hyperparameter optimization,but also a useful feature selection algorithm to improve prediction accuracies of debris flow susceptibility by using machine learning algorithms.The high and very high debris flow susceptibility zone appropriately covers 38.01%of the study area,where debris flow may occur under intensive human activities and heavy rainfall events.

Glacial debris flow susceptibility mapping based on combined models in the Parlung Tsangpo Basin,China

Machine learning(ML)-based prediction models for mapping hazard(e.g.,landslide and debris flow)susceptibility have been widely developed in recent research.However,in some specific areas,ML models have limited application because of the uncertainties in identifying negative samples.The Parlung Tsangpo Basin exemplifies a region prone to recurrent glacial debris flows(GDFs)and is characterized by a prominent landform featuring deep gullies.Considering the limitations of the ML model,we developed and compared two combined statistical models(FA-WE and FA-IC)based on factor analysis(FA),weight of evidence(WE),and the information content(IC)method.The final GDF susceptibility maps were generated by selecting 8 most important static factors and considering the influence of precipitation.The results show that the FA-IC model has the best performance.The areas with a very high susceptibility to GDFs are primarily located in the narrow valley section upstream,on both sides of the valley in the middle and downstream of the Parlung Tsangpo River,and in the narrow valley section of each tributary.These areas encompass 86 gullies and are characterized as"narrow and steep".

Integrated numerical simulation of hydraulic fracturing and production in shale gas well considering gas-water two-phase flow

Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.

Arbitrary High-Order Fully-Decoupled Numerical Schemes for Phase-Field Models of Two-Phase Incompressible Flows

Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.

Numerical Simulation of Oil-Water Two-Phase Flow in Low Permeability Tight Reservoirs Based on Weighted Least Squares Meshless Method

In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.