泥石流固体物质起动临界深度与最大冲出总量测算方法、泥石流规模预报方法

技术领域

本发明涉及泥石流活动规模的测算方法，特别是涉及一种测算泥石流起动时可裹挟的固体物质临界深度的方法，一种泥石流最大冲出总量测算方法，以及基于此的泥石流活动规模预报方法，属于地质灾害防治技术、地质灾害监测预警技术领域。

背景技术

泥石流形成是泥石流物源区固体物质在降水与地形条件耦合作用下随水起动汇集的结果。在一次降水条件下，有多少固体物质能够在水力条件下起动，是衡量泥石流活动规模的有效方法。该指标通常也可称为“最大冲出总量”，具体指一次降水条件下单场泥石流活动所输移的固体物质体积量。

在泥石流防治技术领域，针对上述指标量的测算，现有技术有两种方案：其一、综合分析各类地形参数、固体物质性质参数在一定水力条件下可能参与泥石流活动的固体物质数量，即“动储量”。例如，采用滑塌面积率、不良地质体累计分布长度与主沟总长的比值对松散物质的活动性进行判定，采用灰色系统模型对流域内动储量进行定量预测，利用泥石流流量、输沙量及其与沟床纵坡、岩性之间的关系反演建立了松散固体物质临界集中量模型，综合面积滑塌率、不良地质体长度比、流域内爆发泥石流碎屑物的最低标准三个指标，对泥石流易发程度和规模进行判别，以及依据泥石流历时冲刷深度推导动储量的计算公式等。其二、综合泥石流固体物质冲出总量的多项影响因素建立经验模型。例如，基于流域面积、沟道坡度、沟道高差、物源储量、岩性特征等参数对泥石流固体物质冲出总量进行研究，建立相关预测模型。也有研究结果是在经验模型基础上进一步给出最大冲出量与泥石流物源区总物质储量间的简洁经验关系，例如，将流域内总物源的10％作为水源条件为百年一遇降水条件下单场泥石流的最大冲出量，以及在20年一遇条件下，灾区泥石流沟的最大冲出量与沟道内松散物质储量之间存在很强的幂函数关系。

现有技术在测算最大冲出总量时，存在缺陷在于多是从定性或是半经验半定量的角度进行，尚无统一的度量标准。由此造成因研究区差异而使各类测算的度量精度差异较大，以及各类测算方案普适性不高、应用明显受限。并且对于根据已有样本建立的经验模型，若发生泥石流规模超过已有样本规模，则估算量将严重不足。这进一步导致，现有测算方法在针对偏离常态的黑天鹅式的泥石灾害的规模估计不足。而实际灾害防治中，这类灾害事件又是最需要借助科学测算数据制订灾后应急与救援方案的。

发明内容

本发明的目的就是针对现有技术的不足，提供一种有普遍适用性的泥石流最大冲出总量测算方法。

为实现上述目的，本发明首先提供一种泥石流固体物质起动临界深度测算方法，其技术方案如下：

一种泥石流固体物质起动临界深度测算方法，测算泥石流起动时可裹挟的固体物质临界深度，其特征在于：首先，现场调查获取泥石流物源区基本参数；其次，根据物源区水源条件与固体物质条件判断物源区泥石流起动属于非饱和渗流－径流起动模式，或者是非饱和渗流起动模式，最后，计算确定泥石流固体物质起动临界深度，具体：对于非饱和渗流－径流模式，依式1计算确定泥石流物源起动临界深度H

式中，H

单位m，

H

c—物源区固体物质粘聚力，单位kN/m

γ

h—物源区径流深度，单位m，物源区基本参数确定，

γ

n—物源区固体物质孔隙度，物源区基本参数确定，

θ—物源区固体物质所处斜坡坡度，单位°，物源区基本参数确定，

上述泥石流固体物质起动临界深度测算方法是测算在一次降水条件下泥石流起动时可裹挟的固体物质临界深度。方案原理在于：泥石流固体物质起动临界深度H

上述方法中，现场调查包括了针对工程所在物源区现场的各种测绘、测量、模拟实验测试，以及历史灾害记录获取，以及有参照借鉴作用的经验数据获取等。

上述方法中，物源区径流深度h、物源区固体物质粘聚力c、物源区固体物质内摩擦角

根据本发明泥石流固体物质起动临界深度测算方法可应用于泥石流防治工程设计。例如，根据测算得到的泥石流物源起动临界深度H

上述泥石流固体物质起动临界深度测算方法在泥石流防治工程设计中的应用。

以上述泥石流固体物质起动临界深度测算方法为基础，本发明同时提供一种泥石流最大冲出总量测算方法，其技术方案如下：

一种利用上述泥石流固体物质起动临界深度测算方法实现的泥石流最大冲出总量测算方法，用于测算一次降雨引发泥石流的冲出总量，其特征在于：依如下步骤实施：

步骤S1、现场调查获取泥石流物源区基本参数；

步骤S2、根据泥石流物源区基本参数划分物源区斜坡坡体，确定每一斜坡坡体规格参数，并将泥石流物源区基本参数配置至每一斜坡坡体；

步骤S3、对于每一斜坡坡体，根据水源条件与固体物质条件判断物源区泥石流起动属于非饱和渗流－径流起动模式，或者是非饱和渗流起动模式，计算每一斜坡坡体的泥石流物源起动临界深度H

步骤S4、对于每一斜坡坡体，依式3计算确定泥石流可起动的固体物质体积量V

V

式中，V

D

L

Z

H

步骤S5、依式4计算物源区内一次降雨引发泥石流的冲出总量V

式中，V

上述泥石流最大冲出总量测算方法，首先将根据地形数据将物源区划分为不同的斜坡坡体作为计算单元，其次，针对每一斜坡坡体计算泥石流物源起动临界深度H

本发明同时提供上述泥石流最大冲出总量测算方法在泥石流灾害监测预警中的应用，以及上述泥石流最大冲出总量测算方法在泥石流防治工程设计中的应用。

泥石流最大冲出总量是衡量泥石流活动规模大小简洁且有效的指标。在泥石流起动条件中，地形条件与固体物质特征条件相对稳定，降水条件则存在较高的变化率。因而，对于本发明泥石流最大冲出总量测算方法，若将测算原理中利用降水条件具体为实时动态降水条件，则测算结果是随实时降水条件变化的泥石流的冲出总量V

一种利用上述泥石流最大冲出总量测算方法实现的泥石流活动规模预报方法，其特征在于：根据降雨实时监测数据确定物源区径流深度h、物源区固体物质粘聚力c、物源区固体物质内摩擦角

与现有技术相比，本发明的有益效果是：本发明泥石流固体物质起动临界深度测算方法与最大冲出总量测算方法，是全定量测算方法。测算方法从泥石流起动的机理分析着手，将影响固体物质移动形成泥石流的三大基本条件分别分解为可定量化的量值，再通过具有明晰物理意义的动力学模型建立测算方案。该测算方案的度量依据可靠性高，能够针对性解决现有技术受限于某一特定物源区对象，以及对偏离常值的意外极端泥石流灾害事件测算不足的缺陷，具有更广的适用性。本发明还提供基于泥石流最大冲出总量测算方法建立的泥石流活动规模预报方法，能够将测算方法科学、普适性的有益效果带入泥石流险情预警。

附图说明

图1是泥石流固体物质起动临界深度测算方法技术路线示意图。

图2a、图2b是2013年“7.26”泥石流前5天降雨量及暴发过程中的分钟雨量。

具体实施方式

下面结合附图，对本发明的优选实施例作进一步的描述。

实施例一

如图1～图2所示，采用本发明方法测算四川省都江堰市白沙河流域锅圈岩泥石流沟2013年7月26日降雨事件引发的泥石流中的相关量值。

图1是泥石流固体物质起动临界深度测算方法技术路线示意图。

1、开展现场调查

现场调查包括了针对工程所在物源区现场的各种测绘、测量、模拟实验测试，以及历史灾害记录获取，以及有参照借鉴作用的经验数据获取等。

锅圈岩泥石流沟位于都江堰市区以北约10km处，是深溪沟左岸的一条支沟，隶属于白沙河水系。该泥石流沟分布于龙门山断裂带的中部偏南区段，是汶川地震的核心震区，地震烈度较高，具体等级为XI度。流域(包括沟道和两侧斜坡)特征为：整个流域面积约0.15km

物源区地形数据：现场调查或GIS数据提取的DEM数据。

物源区固体物质性质参数：现场取样，完成必要的现场或室内实验测量确定泥石流源区固体物质饱和容重γ

物源区水动力条件参数：从锅圈岩沟雨量站的雨量监测数据获取降雨过程数据。数据类型为实时监测数据见图2a、图2b所示“7.26”泥石流前5天降雨量及暴发过程中的分钟雨量。

2、划分计算单元

根据物源区斜坡坡体特征划分斜坡坡体作为计算单元。在物源区内划分出4个斜坡坡体(表1)，再将泥石流物源区基本参数配置到每一斜坡坡体(表2)。

表1斜坡坡体规格参数

表2斜坡坡体基本参数

由于本实施例采用实时降雨数据，因而用于测算泥石流规模的物源区径流深度h、物源区固体物质粘聚力c、物源区固体物质内摩擦角

以实时监测降雨数据为基础，采用Richards水分运动方程与VG导水率模型将水源条件与源区固体物质进行耦合，借助hydrus2D软件模拟物源区固体物质的降雨入渗特征过程。当各计算单元坡脚开始出现大量积水(软件模拟水头高度＝0.6m)时，判断为固体物质达到最大临界深度时刻，标记该时刻为时刻t

3、测算每一斜坡坡体泥石流物源起动临界深度

根据模拟结果中的边界通量是否有表层径流数据存在，判别每一斜坡坡体的起动模式。若存在表层径流数据，则属于非饱和渗流－径流起动模式；反之，则属于非饱和渗流起动模式。根据数据判别显示，物源区内4个斜坡坡体的泥石流起动属于非饱和渗流起动模式。将各斜坡坡体的各参数代入式2，计算得到每一斜坡坡体的泥石流固体物质起动临界深度H

4、测算每一斜坡坡体可起动的固体物质体积量V

将各参数代入式3计算得到可起动的固体物质体积量V

表3斜坡坡体H

5、测算物源区内本次降雨条件引发泥石流的冲出总量V

根据表3数据，依式4计算有，物源区内“7.26”降事件引发泥石流的冲出总量V

6、结果验证

将测算结果V