发明名称 泥石流起动临界清水流量测算方法及起动监测方法
摘要 本发明公开了一种沟道侵蚀型泥石流起动的临界清水流量测算方法及起动监测方法。针对现有技术只适用于坡面侵蚀型泥石流起动临界流量测算的缺陷,本发明提供了一种沟道侵蚀型泥石流起动的临界清水流量测算方法。本方法基于沟道侵蚀型泥石流起动时刻沟道粗化层中值粒径的起动概率与沟道地形特征数据的对应关系的原理,通过统计分析泥石流起动时固体物质中值粒径的起动概率,结合起动概率计算方法测算泥石流起动清水流量。本发明进一步提供了沟道侵蚀型泥石流起动降雨量阈值测算方法,以及利用雨量阈值测算监测实现的泥石流起动监测系统。本发明方法原理可靠,计算过程科学简便,能够应用于缺资料地区沟道侵蚀性泥石流起动临界降雨量的确定。
申请公布号 CN103453944B 申请公布日期 2016.01.06
申请号 CN201310403573.4 申请日期 2013.09.06
申请人 中国科学院、水利部成都山地灾害与环境研究所 发明人 王涛;陈宁生;邓明枫;杨成林
分类号 G01F1/00(2006.01)I;G01W1/14(2006.01)I 主分类号 G01F1/00(2006.01)I
代理机构 成都赛恩斯知识产权代理事务所(普通合伙) 51212 代理人 张帆
主权项 一种沟道侵蚀型泥石流起动临界清水流量测算方法,其特征在于:依如下步骤实施:步骤S1、确定沟道侵蚀型泥石流形成区特征过流断面根据调查测绘确定泥石流沟道内泥石流形成区范围,在所述沟道泥石流形成区内的上段区域选择一段沟道顺直、沟床比降均一的沟段作为特征过流沟段,以所述特征过流沟段的中间断面作为特征过流断面;步骤S2、确定地形地质基本参数数据调查测绘测量确定地形地质基本参数数据,所述地形地质基本参数数据包括:特征过流断面松散物质粗化层中值粒径d<sub>50</sub>,特征过流断面宽度B、固体物质容重γ<sub>s</sub>、特征过流沟段比降i、沟道糙率n;步骤S3、确定沟道侵蚀型泥石流起动临界清水流量Q联立式1~式5测算确定沟道内泥石流起动时的临界清水流量Q:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>P</mi><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mn>1</mn><msqrt><mn>2</mn><mi>&pi;</mi></msqrt></mfrac><msubsup><mo>&Integral;</mo><mrow><mo>-</mo><mo>&infin;</mo></mrow><mi>&omega;</mi></msubsup><mi>exp</mi><mo>[</mo><mo>-</mo><mfrac><msup><mi>x</mi><mn>2</mn></msup><mn>2</mn></mfrac><mo>]</mo><mi>dx</mi></mrow>]]></math><img file="FDA0000843640250000013.GIF" wi="745" he="187" /></maths>        式1式中,P—沟道松散固定物质粗化层中值粒径起动概率,取值73%~85%,ω—计算颗粒起动概率的积分上限,依式2计算确定,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&omega;</mi><mo>=</mo><mn>5.03</mn><msqrt><mfrac><mrow><msub><mi>&gamma;</mi><mi>s</mi></msub><mo>-</mo><mi>&gamma;</mi></mrow><mi>&gamma;</mi></mfrac><mi>g</mi><msub><mi>d</mi><mn>50</mn></msub></msqrt><mo>&times;</mo><msup><mrow><mo>[</mo><mfrac><mi>h</mi><msub><mi>d</mi><mn>50</mn></msub></mfrac><mo>]</mo></mrow><mn>0.167</mn></msup><msup><mi>V</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><mn>2.7</mn></mrow>]]></math><img file="FDA0000843640250000012.GIF" wi="967" he="187" /></maths>       式2式中,γ<sub>s</sub>—特征过流断面固体物质容重,g/cm<sup>3</sup>,步骤S2确定,γ—水的容重,g/cm<sup>3</sup>,取常值,g—重力加速度,m/s<sup>2</sup>,取常值,d<sub>50</sub>—特征过流断面粗化层松散物质中值粒径,m,步骤S2确定,V—特征过流断面平均流速,m/s,依式1、2联立计算确定,h—沟道侵蚀型泥石流起动临界水深,m,依式4计算确定;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>V</mi><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><msup><mi>R</mi><mfrac><mn>2</mn><mn>3</mn></mfrac></msup><msup><mi>i</mi><mfrac><mn>1</mn><mn>2</mn></mfrac></msup></mrow>]]></math><img file="FDA0000843640250000023.GIF" wi="324" he="141" /></maths>             式3式中,R—特征过流断面水力半径,m,i—特征过流沟段比降,%,步骤S2确定,n—沟道糙率,%,步骤S2确定;<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>h</mi><mo>=</mo><mfrac><mi>RB</mi><mrow><mi>B</mi><mo>-</mo><mn>2</mn><mi>R</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000843640250000022.GIF" wi="367" he="132" /></maths>         式4式中,B—特征过流断面宽度,m,步骤S2确定,Q=VBh   式5式中,Q—沟道侵蚀型泥石流起动临界清水流量,m<sup>3</sup>/s,B—特征过流断面宽度,m,步骤S2确定。
地址 610041 四川省成都市武侯区人民南路四段九号