一种凹形坡沟道中泥石流侵蚀深度计算方法技术

本发明专利技术公开了一种凹形坡沟道中泥石流侵蚀深度计算方法，包括以下步骤：通过对现有的两种泥石流侵蚀模式的分析，根据凹形坡沟道中泥石流的侵蚀特征，以及水槽实验结果，引入了沟道纵剖面曲率来描述凹形坡沟道形态变化，结合其他相关物理参数，最终得出凹形坡沟道中泥石流的侵蚀深度计算公式。该公式。综合考虑了两种侵蚀模式，能更加全面精确的计算凹形坡沟道中泥石流的侵蚀深度，对实践更具有指导意义。

【技术实现步骤摘要】
一种凹形坡沟道中泥石流侵蚀深度计算方法
：本专利技术涉及泥石流沟道基底侵蚀深度估算领域，具体涉及一种凹形坡沟道中泥石流侵蚀量估算方法。
技术介绍
：泥石流是一种常见的山区地质灾害，其沟道侵蚀作用会显著放大泥石流自身规模，扩大灾害泛滥范围。沟道松散堆积体的存在，会导致泥石流沟道纵剖面形态的改变，出现纵向上的凹点，从而影响后续泥石流灾害的运动特性和沿程侵蚀分布。已有的泥石流侵蚀计算公式中，将凹点对于泥石流侵蚀的影响归结于泥石流体与底层床体之间的附加向心力。由于附加向心力的作用，会使泥石流体与底床物质之间的正应力增加，从而底床土体不容易被携带走，侵蚀率降低。但经过专利技术人的分析，实际中沟道凹点处的侵蚀效应往往更加显著，这是因为在地形凹点处，出现了一种特殊的侵蚀类型——冲击侵蚀，这种侵蚀模式对泥石流侵蚀过程有重要影响。因此，凹形坡沟道中剪切侵蚀和冲击侵蚀并存，现有侵蚀计算方法往往局限于无限变坡模型的剪切侵蚀，难以应用于凹形坡沟道中。
技术实现思路
：为了解决上述技术问题，本专利技术公开了一种凹形坡沟道中泥石流侵蚀深度计算方法，解决现有的侵蚀计算方法适用于坡度值固定的坡道，同时对凹形坡沟道中侵蚀类型考虑不周的问题。为此，本专利技术采用如下方案：一种凹形坡沟道中泥石流侵蚀深度的计算方法，采用如下公式计算凹形坡沟道中泥石流侵蚀深度E，单位为cm，用于计算凹形坡沟道凹点处的深度：其中，泥石流初始速度为ν，m/s；泥石流流深为h，cm；泥石流的容重为ρ，N/m3；凹点处的曲率半径为Rx，m；τb是泥石流施加于凹形坡沟道的底床物质的剪应力；τres是凹形坡沟道的底床物质施加给泥石流的阻力。τb＝ρg(hcosθtanφ+(ν/Cz)2)τres＝c+ρghcosθtanφbed其中，θ为坡面角；φ为泥石流内部的摩擦角；Cz为谢才系数；c为粘聚力；φbed为底床床面摩擦角。与现有技术相比，本专利技术的效果在于：本专利技术引入了曲率半径来描述凹形坡中凹点的坡度转折陡缓程度，同时，侵蚀深度计算包括两部分，即剪切侵蚀和冲击侵蚀，能更合理的应用于凹形坡沟道中。附图说明此处所说明的附图用来提供对本专利技术实施步骤的进一步理解，构成本申请的一部分，并不构成对本专利技术实施案例的限定。图1为凹形坡沟道侵蚀初始状态。图2为泥石流侵蚀动力分析图：(a)均一坡沟道；(b)凹形坡沟道。图3为水槽实验凹形坡简图。图中A点为凹点，α为转折角，分别为15°、25°、35°，Rx为曲率半径；图4为曲率半径定义示意图；图5为凹形坡沟道中侵蚀过程示意图；图6为本专利技术函数f与曲率半径关系示意图具体实施方式以下对本专利技术的技术术语进行说明；本专利技术所指的泥石流初始速度ν，h分别是指泥石流进入沟道某一凹点处之前的龙头流速和泥位高度。本专利技术的泥石流的容重为ρ为单位体积泥石流体的重量。同时，本专利技术所述凹点，是指泥石流变坡沟道中，纵剖面坡度由陡变缓的转折点。凹点处的曲率半径Rx是刻画变坡沟道凹点的凹形程度，坡度转变越急促，说明曲率半径越小。本专利技术凹形坡沟道中泥石流侵蚀深度估算方法，具体思路如下：(1)总结现有的两种泥石流侵蚀模式现有的泥石流侵蚀率计算公式大多基于无限边坡模型，即认为沟床坡度是不变的，上覆泥石流体和下层土体之间的剪切侵蚀是主要的侵蚀方式，侵蚀动力源自层与层之间的剪切应力对τb和τres(如附图2a所示)。实际的泥石流沟道中往往不是均一坡，而是存在沿程的坡度变化。根据野外调查和实验观察，发现凹形沟道中存在另一种显著的侵蚀模式-即冲击侵蚀，侵蚀动力源自泥石流体与静止底床之间的存在着法向的正面冲击法向应力对σb和σres(如附图2b所示)。(2)开展水槽实验设计水槽实验，模拟泥石流侵蚀凹形坡松散物质，如附图2所示，水槽中A点为凹点，其转折角为α＝15°、25°、35°，泥石流的容重、速度、流深分别为ρ,v,h，凹点下游最大侵蚀深度为E。(3)引入曲率半径描述沟道形态为了得到凹形坡沟道中泥石流侵蚀深度，建立了相应的坐标系，初始的床面剖面曲线可以用b＝b(x,z)来表示，侵蚀后的床面曲线用b'＝b'(x,z)表示，如附图3所示。引入坡面角θ和曲率来刻画底床床面的形态特征。坡面角θ指沟道某一点处的坡度。(4)凹形坡沟道中泥石流侵蚀深度的确定当泥石流侵蚀凹形坡沟道中的松散物质时，侵蚀模式不仅包括剪切侵蚀，还存在另一种以冲击和碰撞为特征的冲击侵蚀模式。泥石流剪切侵蚀项E1可以由下式得到其中E1是剪切侵蚀率，τb是泥石流阻力(泥石流施加于底床物质的剪应力)，本专利技术采用Voellmy本构模型。τres是底床物质施加给泥石流体的阻力，本专利技术采用库伦阻力模型。τb＝ρg(hcosθtanφ+(ν/Cz)2)(式4)τres＝c+ρghcosθtanφbed(式5)式中φ为泥石流体内部的摩擦角，c为粘聚力，两者均为土体的基本强度参数，可通过直剪实验确定；Cz为谢才系数，其计算式为Cz＝n/R1/6，其中n为糙率系数，一般取0.02，R为水力半径，其值为h/2；φbed为底床床面摩擦角，其值一般略小于泥石流体的内摩擦角φ。通过上述公式可以计算得到本专利技术需要的τb和τres的值。在凹形坡沟道中，引入额外的一项E2来表示冲击侵蚀项。凹点上游段O1O段和下游段OO2段的坡度分别标记为θ和θ′，泥石流初始速度、流深分别为ν、h。以下游段OO2段为X轴建立了一个局部坐标系X–Z，如附图4所示，假定冲击侵蚀的强弱与泥石流泥深h、速度ν沿Z方向的分速度νz成正比，引入函数f(Rx)来描述这种正比关系。根据式曲率的定义，有上式经过简单的变形得到由此得到冲击侵蚀项E2为因此，凹形坡沟道中泥石流侵蚀深度E的计算式为通过多组实验，通过拟合确定函数f(Rx)的值，拟合图如附图5所示。因此最终的凹形坡沟道中侵蚀深度计算式为步骤(2)中，实验水槽与水平面夹角为10°。步骤(4)中凹点上游段O1O段和下游段OO2段的坡度分别标记为θ和θ′。实施例1本实施例提供一种凹形坡沟道中泥石流侵蚀深度的计算方法，用如下公式计算凹形坡沟道中泥石流侵蚀深度E，用于计算凹形坡沟道凹点处的深度，根据式10，可知:其中，泥石流初始速度为ν，m/s；泥石流流深为h，cm；泥石流的容重为ρ，N/m3。冲击侵蚀的强弱与泥石流这些相关指标有关。凹点处的曲率半径为Rx，m，为底床床面的形态固有特征。τb是泥石流施加于凹形坡沟道的底床物质的剪应力，是泥石流阻力，上覆泥石流体和下层土体之间的剪切侵蚀是主要的侵蚀方式，侵蚀动力源自层与层之间的剪切应力对τb和τres；τres是凹形坡沟道的底床物质施加给泥石本文档来自技高网...

【技术保护点】
1.一种凹形坡沟道中泥石流侵蚀深度的计算方法，其特征在于，采用如下公式计算凹形坡沟道中泥石流侵蚀深度E，单位为cm，用于计算凹形坡沟道凹点处的深度：/n

【技术特征摘要】
1.一种凹形坡沟道中泥石流侵蚀深度的计算方法，其特征在于，采用如下公式计算凹形坡沟道中泥石流侵蚀深度E，单位为cm，用于计算凹形坡沟道凹点处的深度：



其中，泥石流初始速度为ν，m/s；泥石流流深为h，cm；
泥石流容重为ρ，N/m3；
凹点处曲率半径为Rx，m；

【专利技术属性】
技术研发人员：王家鼎，李浦，胡凯衡，王飞，
申请(专利权)人：西北大学，
类型：发明
国别省市：陕西;61