Multi-factors fusion method of debris flow prediction based on genetic programming
-
摘要: 泥石流是一种多发的地质灾害,常对人民生命财产安全带来极大的威胁,其暴发不仅与降雨有关,还与众多地质环境因子相关。本文以流域面积、松散物质比率、沟床平均坡度为地质因子,以最大小时雨强(T)和总降雨量(R)的乘积作为降雨指数,在获取的泥石流地质因子和降雨指数因子综合样本库的基础上,采用遗传规划法建立了泥石流临界降雨指数智能预测模型,克服了以往以雨量为单一指标的预警模型的弊端,模型验证结果显示,泥石流预测精度高、适应性强。Abstract: Debris flow is a frequent geological disaster, which often poses a great threat to the safety of people’s lives and property. Outbreak of debris flow is not only related to rainfall, but also related to many geological and environmental factors. In this paper, the watershed area, ratio of loose materials and the average slope of the gully bed are taken as the geological factors,the maximum hourly rainfall intensity (T) and the total rainfall (R) are taken as the rainfall index, the sample database is established by means of geological factors and rainfall index, genetic programming is used to establish a prediction model for the critical rainfall index of debris flow, which overcomes the shortcomings of the previous warning model that used rainfall as a single indicator, model verification results show that the model has high warning accuracy and strong adaptability, which can realize timely warning.
-
Keywords:
- debris flow /
- geological factors /
- rainfall index /
- genetic programming /
- intelligent prediction
-
2010年8月8日凌晨,舟曲县城北侧三眼峪沟和罗家峪沟同时暴发特大山洪泥石流,城区三分之一被淹,共造成1435人死亡,330人失踪,直接经济损失超过10亿元[1-2],给舟曲县城居民生命财产造成了巨大损失,也给当地生产生活带来严重困难。灾害引起了党中央、国务院、中央军委及全国人民的高度关注,启动国家二级救灾应急响应,同时批复专项资金进行治理。其中,重力式拦挡坝是舟曲泥石流治理中的最主要工程之一。
周龙茂等[3]认为拦挡坝在泥石流治理中发挥着重要作用,但同时拦挡坝也是最易遭受破坏、失去防灾功能的泥石流防治构筑物,因此,在泥石流设计中,为了拦挡坝不被破坏,往往设计得很保守,王念秦等[4]提出这种设计容易造成两个极端现象:保守,造成资金浪费;冒进,防治工程失败。要做到既能保证拦挡坝安全,又能将投资最小化,就要求对泥石流拦挡坝抗冲击力验算方法提出新要求。传统的泥石流冲击力计算经验公式[5]只能通过大量试算表述结果,不能表述过程。而将三维有限元数值分析方法应用到泥石流拦挡坝稳定性验算中[6-7],能将过程和结果同时呈现,即通过分步加载的方法,逐步呈现拦挡坝的位移情况和抗冲击力过程中的破损情况。
关于拦挡坝的数值模拟研究还比较少,本文将考虑损伤的混凝土本构模型与有限元计算方法相结合,对舟曲泥石流混凝土拦挡结构进行力学分析,最终确定了拦挡坝的抗冲击力合理区间,以期为泥石流治理工程的设计提供借鉴。
1. 模型建立及相关参数的设定
杨东旭等[8]、许海亮等[9]、张睿骁等[10]认为冲击力是破坏防治工程构筑物的主要作用力之一,其大小与泥石流流量、流速、容重等有关。泥石流冲击力是泥石流防治工程设计的重要参数,分为流体整体冲击力和个别石块的冲击力两种,在设计中取两种计算结果较高者为设计依据。文章采用《泥石流灾害防治工程设计规范》(DZT 0239—2004)中的经验公式[5]作为数值计算结果的参考和验证。
流体整体冲击力计算公式:
$$f = K\frac{{{\gamma_C}}}{g}v_c^2$$ (1) 式(1)中:f−冲击力/Pa;
K−系数,取2.5;
${\gamma_C} $ −泥石流重度/(t·m−3);g−重力加速度,取9.8 m·s−2;
vc−断面处泥石流流速(m·s−1)。
个别石块的冲击力计算公式:
$${F_b} = \sqrt {\frac{{48EJ{V^2}W}}{{g{L^3}}}} \cdot \sin \alpha $$ (2) 式(2)中:Fb−泥石流大石块冲击力/(t·m−2);
E−工程构件弹性模量/(t·m−2);
J−工程构件界面中心轴的惯性矩/m4;
V−石块运动速度(m·s−1);
W−石块重量/t;
L−构件长度/m;
α−石块运动方向与构件受力面的夹角/(°)。
泥石流具体参数和计算结果见表1,编号和《甘肃省舟曲县三眼峪沟泥石流灾害设计报告》[1]中保持一致。
表 1 泥石流冲击力计算参数及结果Table 1. Debris flow impact calculation parameters and results编号 ${\gamma _c}$/(t·m−3) ${v_c}$/(m·s−1) $L$/m $W$/t $V$/(m·s−1) $\sin \alpha $ $E$/(t·m−2) $J$/m4 $f$/(t·m−2) ${F_b}$/(t·m−2) 大1号坝 2.09 6.56 5.5 162 10.55 0.946 2.8 274.63 22.49 7.14 大2号坝 2.09 6.60 2.8 81 7.53 0.946 2.8 274.63 22.76 7.34 大3号坝 2.03 6.61 8.6 259.2 13.20 0.946 2.8 216.00 22.17 19.95 大4号坝 2.03 8.67 4.3 94.5 9.33 0.946 2.8 421.88 38.15 19.62 小2号坝 2.03 6.56 3.1 83.7 7.92 0.946 2.8 421.88 21.84 4.28 小3号坝 2.13 9.86 5.8 129.6 10.84 0.946 2.8 274.63 51.77 8.98 小4号坝 2.13 5.49 2.8 75.6 7.53 0.946 2.8 421.88 16.05 14.17 小6号坝 2.13 8.12 4.7 121.5 9.76 1.000 2.8 512.00 35.11 37.01 主1号坝 2.13 6.53 11.2 361.8 15.06 0.946 2.8 421.88 22.71 44.35 主2号坝 2.13 6.16 7.5 234.9 12.32 0.946 2.8 343.00 20.21 12.44 计算区域按地质资料分高程、分区域模拟。模型向上游及下游分别延伸至坝体厚度的2倍,模型高度方向自坝基向下延伸坝体垂直部分的2倍。以大2号坝为例进行数值模拟,砼坝体和地基土数值分析具体参数见表2。
分析中将泥石流流体的冲击力P简化为静力加载到坝体侧面,计算坝体的极限抗冲击能力,简化计算力学模型如图1所示。
表 2 混凝土坝和地基碎石土参数Table 2. Concrete dam and gravel soil parameter名称 坝体 沟床碎石土 弹性模量/GPa 泊松比 损伤阈值 拉压强度比 弹性模量/MPa 泊松比 黏聚力/kPa 内摩擦角/(°) 剪胀角/(°) 取值 24.0 0.2 2×10−4 0.15 240 0.2 5.0 40 40 2. 重力式拦挡坝抗冲击力数值分析
数值模拟使用有限元软件ABAQUS进行计算分析。冯帅等[11]认为数值计算出的泥石流的极限抗冲压力大于经验公式计算出的泥石流整体冲击压力。因此,本文分别按泥石流流体高度h=H/2坝高(工况1)及h=H(工况2)两种工况进行分析,将数值计算出的抗冲击力限定在一合理区间。计算中加载每一荷载增量后均计算至收敛,并记录坝体的最大位移。加载至破坏时(计算不收敛,坝体位移不断增大)的压力P与坝体最大位移曲线由倾斜直线变为水平线,坝体所能承受的最大冲击力Pu可由压力P与坝体最大位移曲线水平段的和坐标求出。应力以拉为正,以压为负,应力的单位为Pa,长度单位为m,位移单位为mm,其他单位均采用国际单位制。具体模型边界和网格划分图2所示。
2.1 工况1条件下重力式拦挡坝抗冲击力数值分析
大坝按泥石流冲击高度h=H/2坝高计算,在拦挡坝的一侧施加泥石流冲击力,荷载增量取值100 kPa,每加一次荷载,计算至收敛,并且记录一次坝体的最大位移;一直持续加载至破坏时,即计算不收敛且坝体位移不断增大时停止计算。
周勇等[12]采用结构动力学的方法,建立了泥石流冲击荷载与拦挡坝的动力方程,提出拦挡坝的坝顶处有最大的位移,为本次工程力学计算提供了一种思路。将泥石流冲击力与相应荷载下坝体位移进行统计,形成图3所示冲击力与坝体最大位移关系曲线,会发现坝体水平位移随着拦挡坝冲击力增大呈对数曲线递增,泥石流流体高度h=H/2工况条件下施加的最大冲击力Pu=3500 kPa(357.14 t/m2)。
坝体损毁过程如图4所示,冲击力达到800 kPa(水平位移1 mm)时坝体开始出现初始损伤;冲击力达到1100 kPa(水平位移1.4 mm)时两侧坝肩和基础局部都出现较明显损伤;冲击力达到1400 kPa时(水平位移2 mm),泄水孔和泄水涵洞边缘出现局部损伤;冲击力达到2700 kPa(水平位移4.8 mm)时,基础出现大面积损伤,沿正面泄水涵洞和泄水孔形成纵向损伤;冲击力达到3200 kPa(水平位移7.2 mm)时,坝肩、基础及坝体正中损伤贯通,损伤区呈“W”型,坝体已基本失去功能,在工程实际应用中已达到破坏极限;冲击力达到3500 kPa时,坝体大面积损伤,超过坝体总面积的2/3,水平位移高达14 mm,整体性降低或消失,这只是一种模拟现象,在工程实际应用中泥石流物质沿坝体破坏处流通,坝体已不存在整体位移现象。
2.2 工况2条件下重力式拦挡坝极限抗冲压数值分析
大坝按泥石流冲击高度h=H坝高计算,荷载增量取值125 kPa,每加一次荷载,计算至收敛,同时记录一次坝体的最大位移,将泥石流冲击力与相应荷载下坝体位移进行统计,形成图5所示冲击力与坝体最大位移关系曲线,会发现坝体水平位移随着拦挡坝冲击力增大也呈对数曲线递增,泥石流流体高度h=H工况条件下施加的最大冲击力Pu=2490 kPa(254.08 t/m2)。
坝体损毁过程如图6所示,冲击力达到375 kPa(水平位移1.7 mm)时坝体开始出现初始损伤;冲击力达到500 kPa(水平位移2.3 mm)时两侧坝肩和基础局部都出现较明显损伤;冲击力达到1250 kPa(水平位移6.6 mm)时,泄水孔和泄水涵洞边缘出现局部损伤;冲击力达到1500 kPa(水平位移8.4 mm)时,基础出现局部损伤,沿正面泄水涵洞和泄水孔形成纵向损伤,两侧坝肩损伤较严重;冲击力达到2000 kPa(水平位移14.8 mm)时,严重损伤区呈“W”型,坝肩、基础和坝体中心部位损伤基本贯通,坝体已基本失去功能,在工程实际应用中已达到破坏极限;冲击力达到2490 kPa时,坝体大面积损伤,超过坝体总量的2/3,水平位移高达40.4 mm,整体性降低或消失,这也只是一种模拟现象,在工程实际应用中泥石流物质沿坝体破坏处流通,坝体已不存在整体位移。
3. 坝体冲击力安全验算
同样方法计算三眼峪及各支沟泥石流重力式拦挡工程冲击力,得到表3,结合表1可以看出,泥石流经验公式计算的单宽冲击力均小于工况1数值模拟验算结果,均大于工况2数值模拟验算结果,与工况1和2的平均值接近。2012年建成至今,大部分坝体库容淤积过半,部分甚至已淤满,说明拦挡坝经受住了各种泥石流冲击破坏的考验。
表 3 重力坝冲击力数值模拟验算结果对比表(单位:t/m2)Table 3. Comparison of results of numerical simulation of impact of gravity dam (unit: t/m2)编号 经验公式
计算冲击力h=H/2
最大冲击力h=H
最大冲击力计算平均值 大1号坝 22.49 27.00 17.55 22.28 大2号坝 22.76 27.36 17.78 22.57 大3号坝 22.17 26.64 17.32 21.98 大4号坝 38.15 45.84 29.80 37.82 小2号坝 21.84 26.28 17.08 21.68 小3号坝 51.77 62.16 40.40 51.28 小4号坝 16.05 19.32 12.56 15.94 小6号坝 37.01 44.52 28.94 36.73 主1号坝 44.35 53.28 34.63 43.96 主2号坝 20.21 24.36 15.83 20.10 4. 讨论
4.1 泥石流冲击高度对坝体影响
将2种工况进行比较,工况1的冲击力达到800 kPa时坝体开始出现初始损伤,最大冲击力为3500 kPa,而工况2的冲击力达到375 kPa时桩体开始出现初始损伤,最大冲击力为2490 kPa。说明冲击高度对坝体的影响较大,随高度增加,达到初损的冲击力荷载几乎成倍数减少,而最大冲击力也减少1000 kPa。这说明在拦挡坝设计中泄水涵洞和泄水孔的预留很重要,为减少拦挡坝的冲击破坏,应尽量选择低坝,同时在不影响坝体安全和停淤功能的基础上,应多布设泄水涵洞和泄水孔,降低坝前壅水位,最大可能避免或减少高水位过流。
4.2 拦挡坝损伤过程对拦挡坝设计的指导意义
拦挡坝的损伤从两侧坝肩开始,再到拦挡坝中间部位的泄水涵洞及泄水孔边缘,然后到基础,最后坝肩、中间部位和基础形成“W”型的贯通破坏,其破损部位按先后顺序依次为“坝肩—泄水涵洞及泄水孔—基础—“W”型贯通”4个过程。设计时应重点考虑这几处薄弱环节,要相应的进行专门的加固处理。
4.3 拦挡坝冲击力设计的合理范围
从安全和经济方面考虑,在拦挡坝设计中冲击力的考虑应该取工况1和工况2之间值较合理,工况1存在风险,工况2偏保守,而工况1和2的中间值接近经验公式计算冲击力。从表3中可知,三眼峪设计中的冲击力选择也是工况1和2的平均值,从而保障了拦挡坝的安全运行。
4.4 可作为现有泥石流设计理论的有效补充
由于真实模拟泥石流重力式拦挡坝的室内大型实验难度比较大,野外测定随机性太大,故本文采用数值模拟的方式来分析三眼峪沟拦挡坝的受力情况。结果表明数值模拟对分析问题有一定的指导意义,可以和现有的泥石流设计理论结合,为以后工程设计提供安全对比,但是不能代替物理实验和理论计算。
5. 结论
(1)本文通过有限元软件ABAQUS进行拦挡坝数值模拟计算和分析,比较2种工况条件发现:随着泥石流冲击高度增加,达到初损的冲击力荷载成倍数减少,而最大冲击力也减少1000 kPa;拦挡坝的损伤从两侧坝肩开始,再到拦挡坝中间部位的泄水涵洞及泄水孔边缘,形成“W”型的贯通破坏。
(2)通过与经验公式计算冲击力比较,发现拦挡坝设计中冲击力选择工况1和2的平均值较合理,可以为工程设计提供安全对比。
-
[1] 崔鹏, 杨坤, 陈杰. 前期降雨对泥石流形成的贡献:以蒋家沟泥石流形成为例[J]. 中国水土保持科学, 2003, 1(1):11-15. [CUI P, YANG K, CHEN J. Relationship between occurrence of debris flow and antecedent precipitation:taking the Jiangjia gully as an example[J]. Science of Soil and Water Conservation, 2003, 1(1):11-15.(in Chinese)]
[2] 韦方强, 胡凯衡, 陈杰. 泥石流预报中前期有效降水量的确定[J]. 山地学报, 2005, 23(4):453-457. [WEI F Q, HU K H, CHEN J. Determination of effective antecedent rainfall for debris flow forecast[J]. Journal of Mountain Science, 2005, 23(4):453-457.(in Chinese)]
[3] 丛威青, 潘懋, 李铁锋, 等. 降雨型泥石流临界雨量定量分析[J]. 岩石力学与工程学报, 2006, 25(增刊1):2808-2812.[CONG W Q, PAN M, LI T F, et al. Quantitative analysis of critical rainfall-triggered debris flows[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25 (Sup1):2808-2812.(in Chinese)]
[4] 田冰, 王裕宜, 洪勇. 泥石流预报中前期降水量与始发日降水量的权重关系:以云南省蒋家沟为例[J]. 水土保持通报, 2008, 28(2):71-75. [TIAN B, WANG Y Y, HONG Y. Weighted relation between antecedent rainfall and process precipitation in debris flow prediction-A case study of Jiangjia gully in Yunnan Province[J]. Bulletin of Soil and Water Conservation, 2008, 28(2):71-75.(in Chinese)]
[5] 刘艳辉, 唐灿, 李铁锋, 等. 地质灾害与降雨雨型的关系研究[J]. 工程地质学报, 2009, 17(5):656-661. [LIU Y H, TANG C, LI T F, et al. Statistical relations between geo-hazards and rain-type[J]. Journal of Engineering Geology, 2009, 17(5):656-661.(in Chinese)]
[6] 倪化勇, 王德伟. 基于雨量(强)条件的泥石流预测预报研究现状、问题与建议[J]. 灾害学, 2010, 25(1):124-128. [NI H Y, WANG D W. Present status, problem and advice on the research of prediction and forecasting of debris flow based on rainfall condition[J]. Journal of Catastrophology, 2010, 25(1):124-128.(in Chinese)]
[7] 梁光模, 姚令侃. 确定暴雨泥石流临界雨量的研究[J]. 路基工程, 2008(6):3-5.[LIANG G M, YAO L K. Study on determining the critical rainfall of rainstorm mudflow[J]. Subgrade Engineering, 2008 (6):3-5.(in Chinese)]
[8] 白利平, 孙佳丽, 南赟. 北京地区泥石流灾害临界雨量阈值分析[J]. 地质通报, 2008, 27(5):674-680. [BAI L P, SUN J L, NAN Y. Analysis of the critical rainfall thresholds for mudflow in Beijing, China[J]. Geological Bulletin of China, 2008, 27(5):674-680.(in Chinese)]
[9] 王海芝. 北京山区基于历史资料的泥石流临界雨量研究[J]. 城市地质, 2008, 3(1):18-21. [WANG H Z. A study of the critical rainfall volume for mudflows based on his-torical data in the mountainous area of Beijing[J]. Urban Geology, 2008, 3(1):18-21.(in Chinese)]
[10] 涂剑, 马超, 杨海龙. 北京山区暴雨泥石流激发雨量条件[J]. 中国水土保持科学, 2017, 15(5):103-110. [TU J, MA C, YANG H L. Rainfall condition of triggering debris flows in Beijing mountain regions[J]. Science of Soil and Water Conservation, 2017, 15(5):103-110.(in Chinese)]
[11] 贾三满, 路璐, 翟淑花, 等. 北京山区泥石流预警阈值初步研究[J]. 城市地质, 2016, 11(3):1-7. [JIA S M, LU L, ZHAI S H, et al. Preliminary study on early warning threshold of debris flow in Beijing mountainous area[J]. Urban Geology, 2016, 11(3):1-7.(in Chinese)]
[12] 王军, 李学峰. 基于Logistic回归的泥石流灾害预警模型[J]. 厦门理工学院学报, 2018, 26(3):73-77. [WANG J, LI X F. Logistic regression for early warning of debris flow hazards[J]. Journal of Lujiang University, 2018, 26(3):73-77.(in Chinese)]
[13] 何朝阳, 许强, 巨能攀, 等. 基于降雨过程自动识别的泥石流实时预警技术[J]. 工程地质学报, 2018, 26(3):703-710. [HE C Y, XU Q, JU N P, et al. Real-time early warning technology of debris flow based on automatic identification of rainfall process[J]. Journal of Engineering Geology, 2018, 26(3):703-710.(in Chinese)]
[14] 王之君, 李仁年, 拓万全. 基于历史观测数据的泥石流突变起动模式研究——以陇南山区白龙江流域为例[J]. 灾害学, 2019, 34(2):140-144. [WANG Z J,LI R N, TUO W Q.Historical field observation data-based study on the catastrophic initiation mode of debris flow:a case of the Bailong river basin in Longnan mountain areas[J]. Journal of Catastrophology, 2019, 34(2):140-144.(in Chinese)]
[15] 王进孝.土石流临界降雨线极限状态之研究-基因演算最佳化类神经网路[D].台湾:国立台湾大学,2004.[WANG J X. Study on the critical line of debris flow based on genetic evolutional neural networks[D].Taiwan:National Taiwan University, 2004.(in Chinese)] -
期刊类型引用(8)
1. 李晓玮,刘翠娜. 基于AHP法山区公路边坡稳定性评价及危害性分析. 河北地质大学学报. 2024(02): 61-65+134 . 百度学术
2. 林毅,马强,张凯,李岚,李嘉宁,陆井龙,肇毓. 有效降雨量在辽宁省内公路水毁风险评价中的应用. 中国地质灾害与防治学报. 2024(04): 85-92 . 本站查看
3. 谷梦瑶,徐新胜,何雨辰. 面向变电站机器人巡检的加权嵌套决策树数据质量协同评价方法. 智慧电力. 2024(09): 65-71 . 百度学术
4. 陈红旗,方志伟,祁小博. 地质灾害危害程度分级标准优化研究. 灾害学. 2023(04): 18-22 . 百度学术
5. 胡志瑞,范宏涛,马欣,黄玮,赵大江,李鹏,李小琼. 高分辨率格网模型在地质灾害风险普查数据坐标转换中的应用研究. 宁夏大学学报(自然科学版). 2023(03): 275-281 . 百度学术
6. 任静. 流域超标准洪涝灾害数据库的构建研究. 农业灾害研究. 2023(11): 302-304 . 百度学术
7. 惠春,宁立新,程昌秀. 中国东南沿海城市海啸灾害风险评估. 自然灾害学报. 2022(01): 49-59 . 百度学术
8. 杨志全,杨溢,朱颖彦. “地质灾害防治”课程在人才培养中的目标定位及面临的主要问题与解决对策探讨. 大学. 2020(33): 53-54 . 百度学术
其他类型引用(6)
计量
- 文章访问数: 370
- HTML全文浏览量: 171
- PDF下载量: 442
- 被引次数: 14