一种河流漂木总量的计算方法与流程

1.本发明涉及漂木灾害技术领域，具体涉及一种河流漂木总量的计算方法。


背景技术：

2.漂木是指存在于山地系统中由滑坡崩塌、河岸侵蚀、火灾、风力作用、自然死亡或人为砍伐等各种因素产生的大量枯倒木。在河流运动过程中，漂木多以单个或批量的形式出现，会造成巨大危害，特别是在洪水期间，漂木的堆积或堵塞对桥梁、涵洞等河流基础设施造成严重威胁。漂木在结构物前的堆积将堵塞过水通道，减小过流断面，导致回水上升，从而导致洪水泛滥面积增大、持续时间延长；并伴随河床淤高、河道撕裂和局部冲刷侵蚀，对结构稳定性十分不利。因此，漂木量的估计对于评价漂木危害及河道防护结构的设计说至关重要的作用。如何估计流域内漂木储量，是研究漂木灾害效应的基础，对于山区小流域暴雨洪峰灾害风险评估及山区工程设施有效防灾减灾具有至关重要作用，对判断桥涵上游漂木数量以及桥涵孔径设计有重要指导性意义。
3.现阶段流域内漂木储量大多采用经验公式来获取漂木源的估计，存在准确性不够的问题。


技术实现要素：

4.本发明所要解决的技术问题是提供一种从汇集机制出发，可以提高计算准确性的河流漂木总量的计算方法。
5.本发明解决其技术问题所采用的技术方案是：一种河流漂木总量的计算方法，所述计算方法包括：
6.来自河岸坡面的体积输入量lw
in,hs
的计算和和来自河岸侵蚀的体积输入量lw
in,fc
的计算；
7.所述lw
in,hs
的计算包括：坡面林木的死亡砍伐量im的计算、外界因素作用导致的林木倾倒量if的计算以及滑坡崩塌泥石流引起的物源补充量is的计算；
8.所述lw
in,fc
的计算包括：洪水诱发的河岸侵蚀量i
be
的计算以及洪水运动过程中对堆积木材的挟带量ie的计算。
9.进一步地，所述im的计算包括：确定河岸单位面积坡面上成活和死亡林木的体积b
l
、坡面林木的死亡率m、坡面林木的平均高度h以及坡面林木因自然死亡或零星砍伐而倾倒于河道中的体积或长度分数；
10.所述im＝(b
l
×m×h×
pm)n1；其中，n1取1或2。
11.进一步地，所述if的计算包括：确定外界因素作用前单位面积河岸坡面上成活的林木的体积bf、一定时期内的坡面上林木倾倒率tf、成活林木的平均高度hf以及坡面上林木因外界因素作用而倾倒于河道中的体积或长度分数pf；
12.所述if＝(bf×
tf×
hf×
pf)n2；其中，n2取1或2。
13.进一步地，所述is的计算包括：确定河岸坡面崩塌区域内成活和死亡的林木体积ss
、与河道相交的滑坡崩塌或泥石流区的个数ns、滑坡崩塌或泥石流区域面积as以及滑坡崩塌或泥石流的平均重现期ts；
14.所述is＝(ss×ns
×as
×
t
s-1
)rc，其中，rc为夹带比例。
15.进一步地，所述i
be
的计算包括：确定河岸侵蚀前单位面积坡面上成活和死亡的林木体积b
l
、河岸(横向距离)侵蚀率e以及河道侧坡面上林木因河岸侵蚀而倾倒于河道中的体积或长度分数p
be
；
16.所述i
be
＝(b
l
×e×
p
be
)n3，其中，n3取1或2。
17.本发明的有益效果是：本发明将河流流域内漂木汇集来源划分为坡面源和河岸源两大类，分别基于不同汇集机制建立评估模型，实现了对漂木源的定量计算，可以更为准确地计算出河流流域内可能产生的漂木总量。该计算方法对于准确估计山区小流域内漂木储量，构建漂木在桥涵处的堆积堵塞预测分析模型，量化漂木堵塞桥涵诱发的回水上升高度具有重要作用。
附图说明
18.图1是本发明的计算流程图；
具体实施方式
19.下面结合附图和实施例对本发明进一步说明。
20.如图1所示，本发明的一种河流漂木总量的计算方法，所述计算方法包括：来自河岸坡面的体积输入量lw
in,hs
的计算和来自河岸侵蚀的体积输入量lw
in,fc
的计算。
21.所述lw
in,hs
的计算包括：河岸坡面上因自然死亡和零星砍伐带来的死亡砍伐量im的计算，外界因素作用导致的林木倾倒量if的计算以及滑坡崩塌泥石流引起的物源补充量is的计算；
22.所述lw
in,fc
的计算包括：确定洪水诱发的河岸侵蚀量i
be
的计算以及洪水运动过程中对堆积木材的挟带量ie的计算。
23.计算出上述参数后，可以计算出河流可以产生的漂木总量lw
in
，lw
in
＝im+if+is+i
be
+ie。
24.本发明将河流流域内漂木汇集来源划分为坡面源和河岸源两大类，分别基于不同汇集机制建立评估模型，实现了对漂木源的定量计算，可以更为准确地计算出河流流域内可能产生的漂木总量。该计算方法对于准确估计山区小流域内漂木储量，构建漂木在桥涵处的堆积堵塞预测分析模型，量化漂木堵塞桥涵诱发的回水上升高度具有重要作用。
25.本发明中的外界因素包括但不限于火灾因素、暴风雨因素和风力因素。
26.im可通过如下方式计算：确定河岸单位面积坡面上成活和死亡林木的体积b
l
、坡面林木的死亡率m、坡面林木的平均高度h以及坡面林木因自然死亡或零星砍伐而倾倒于河道中的体积或长度分数pm；所述im＝(b
l
×m×h×
pm)n1；其中，n1取1或2。当河道两侧坡面均有森林覆盖时n1取2，单侧森林覆盖时n1取1。
27.一般情况下，0《pm《0.5。pm的计算需要考虑河道宽度、林木高度、林木的倾倒角度及其与河道的距离，可采用sickle和gregory提出的衰减模型来进行估计。例如：采用该方法得到对于15m宽的河道河道，平均高度为50m的林木其平均pm约为0.1，而对于5m宽的河道，
其平均pm约为0.05。此外，pm值随距离河道的距离而急剧下降，且与树高有关。在不适用随机衰减模型假设的地形或在相对较短的长度范围内进行研究的地形当中，可使用实地测量来定义上述参数值。
28.河岸坡面上林木的平均高度h、河岸单位面积坡面上成活和死亡林木的体积b
l
可基于lidar dsm和实地调查进行联合估计，取得每个径阶的株数及树高曲线(可通过林分生长曲线表或richards生长曲线模型获取)，根据每个径阶的中值和树高曲线上查出的该径阶的平均高值，即可从二元材积表中查出各径阶的单株平均材积，乘以株数得到该径阶的林木材积。各径阶林木材积相加即为该区域内林分体积。坡面林木的死亡率m可通过当地林业部门获取数据。
29.if可通过如下方式计算：确定外界因素作用前单位面积坡面上成活的林木体积bf、一定时期(如：一年、两年、一月)内坡面林木倾倒率tf、成活林木的平均高度hf以及坡面林木因外界因素作用而倾倒于河道中的体积或长度分数pf；所述if＝(bf×
tf×
hf×
pf)n2；其中，n2取1或2。当河道两侧坡面均有森林覆盖时n2取2，单侧森林覆盖时n2取1。
30.hf、bf也可基于lidar dsm和实地调查进行联合估计，方法同h及b
l
。tf也可通过基于多期遥感数据和长期定位观察获取数据。一般情况下，0《pf《0.5，pf的确定方法与也可采用sickle和gregory提出的衰减模型来进行估计，同pm。
31.is可通过如下方式计算：确定河岸坡面崩塌区域内成活和死亡的林木体积ss、与河道相交的滑坡崩塌或泥石流区的个数ns、滑坡崩塌或泥石流区域面积as以及滑坡崩塌或泥石流的平均重现期ts；所述is＝(ss×ns
×as
×
ts)rc，其中，rc为夹带比例。一般认为坡面产生崩塌后，整个崩塌区域内林木全部进入河道河道，即rc＝1。
32.ss也可基于lidar dsm和实地调查进行联合估计，方法同b
l
和bf。ns和as也可通过基于多期遥感数据和长期定位观察获取数据。ts可依据该河流流域历史数据。
33.i
be
可通过如下方式计算：确定河岸侵蚀前单位面积坡面上成活和死亡的林木体积b
l
、河岸侵蚀率e以及河道侧坡面上林木因河岸侵蚀而倾倒于河道中的体积或长度分数p
be
；所述i
be
＝(b
l
×e×
p
be
)n3，其中，n3取1或2。当河道两侧坡面均有森林覆盖时n3取2，单侧森林覆盖时n3取1。河岸侵蚀率e主要指河道宽度方向(横向距离)
34.一般来讲，p
be
要比pf和pm的值要大，因为所有被河岸侵蚀的林木都紧挨着河道，而被侵蚀的林木则倾向于向河道内倾倒，所以有，0《p
be
《1.0。p
be
也可采用sickle和gregory提出的衰减模型来进行估计。b
l
计算方法如上所述。e也可基于多期遥感数据和长期定位观察获取数据。
35.对于ie，由于该部分漂木源输入缺乏调查数据且很难量化，因此一般情况下认为ie＝0。