[发明专利]沟谷型泥石流沟道流通区三维地形演化分析方法

摘要

本发明公开了一种沟谷型泥石流沟道流通区三维地形演化分析方法，属于地质灾害防灾减灾技术领域，本发明方法在泥石流流通区沿沟道中泓线每隔一定距离采集沟道横断面地形高程，基于离散化思想将各个横断面进行等距离分段，通过测得每个分段上任意时间对应的流速和泥深，计算每个横断面分段的沟道侵蚀速率及泥石流过程中的侵蚀深度，获取每个横断面泥石流侵蚀后的地形高程，最后将相邻横断面地形等高点通过直线连接，可得到沟道流通区沟床在泥石流冲刷作用下的地形演化，本发明方法不仅克服了传统方法在泥石流侵蚀计算中经验参数取值的主观性和不确定性，同时可得到沟道流通区地形经泥石流冲刷后的三维地形，可为泥石流防灾减灾工作提供参考。

主权项

1.一种沟谷型泥石流沟道流通区三维地形演化分析方法，其特征在于，包括以下步骤：步骤一、沟道横断面地形高程Z(i)的采集；通过现场调查确定沟谷型泥石流流通区范围，并沿沟道中泓线每隔一定距离L获取沟道横断面；在横断面方向上将各横断面水平等距离分为n段，并获取每段沟床表面的高程为Z(i)，单位m，其中，i表示沟床横断面分段标号，i＝1，2，…，n；步骤二、计算任意时间每个横断面对应分段的泥石流瞬态侵蚀速率E(i,t)；计算横断面对应分段i在任意时间t的泥石流流速v(i,t)和泥深h(i，t)，并按照如下计算公式计算横断面对应分段i的瞬态侵蚀速率E(i,t)：其中，为g重力加速度，取值为9.8N/kg，θ为横断面位置的沟道沿泥石流流动方向的坡度，B_D为沟床土体的孔隙水压力系数，为沟床土体的内摩擦角，为泥石流颗粒物质的内摩擦角，β为泥石流流速沿深度方向的分布系数；步骤三、计算任意时间每个横断面对应分段的沟床侵蚀后的地表高程Z′(i)；按照下式计算泥石流对横断面分段i的侵蚀深度dsc及沟床横断面侵蚀后的地表高程Z′(i)；Z'(i)＝Z(i)‑dsc(i)其中，T为泥石流过程的总历时；步骤四、泥石流侵蚀后的地形等高线构建，得到沟道演化后的三维地形；计算每个横断面在泥石流侵蚀作用下的地形变化后，相邻两个横断面地形等高点通过直线连接，可得到流通区沟道在泥石流侵蚀作用下的地形等高线。