Probabilistic inverse-analysis and reliability prediction of rainfall-induced landslides for slope with multi-source information
-
摘要:
概率反分析是推断不确定土体参数统计特征的重要手段,可以使边坡可靠度评估更接近工程实际。然而目前的概率反分析很少使用多源信息(包括监测数据、观测信息和边坡服役记录),因为这通常涉及数千个随机变量和高维似然函数的评估。因此融合多源信息对空间变异土体参数进行概率反分析进而预测降雨条件下的边坡可靠度是一项具有挑战性的难题。文章将改进的基于子集模拟的贝叶斯更新(mBUS)方法与自适应条件抽样(aCS)算法相结合,构建了空间变异土体参数概率反分析和边坡可靠度预测的框架,并以某一公路边坡为例验证了该框架的有效性。研究结果表明:通过融合多源信息所获得的土体参数后验统计特征与现场观测结果基本吻合;用更新后的土体参数预测得到2004年9月12日该边坡在暴雨工况下的失效概率为23.1%,符合实际边坡失稳情况,说明在此框架下可以充分利用多源信息解决高维概率反分析问题。
Abstract:Probabilistic inverse-analysis is an essential approach to infer statistical characteristics of uncertain soil parameters, making the slope reliability assessment closer to engineering reality. However, current probabilistic inverse analysis rarely integrates multi-source information, including monitored data, field observation information, and slope survival records. Conducting the probabilistic inverse-analysis of spatially varying soil parameters and slope reliability prediction under rainfalls by integrating the multi-source information is a challenging issue due to the involvement of thousands of random variables and the evaluation of high-dimensional likelihood functions. In this paper, a modified Bayesian updating with subset simulation (mBUS) method is combined with adaptive conditional sampling (aCS) algorithm to establish a framework for probabilistic inverse analysis of spatially variable soil parameters and reliability prediction of slopes. The effectiveness of this framework is validated using a highway slope as a case study. The research results show that the posterior statistical characteristics of soil parameters obtained by integrating multi-source information are in good agreement with field observation results. Additionally, the probability of slope failure under heavy rainfall on September 12, 2004 with the updated soil parameters is 23.1 %, which is in line with the actual slope instability. Within this framework, multi-source information can be fully utilized to address high-dimensional probabilistic inverse analysis problems.
-
0. 引言
近年来,随着国家对天然气管道工程建设的重视力度不断提高,我国的天然气管道工程建设步伐不断加快,目前国家干线管道系统已初具规模,并加速形成“主干互联、区域成网”的全国天然气基础网络[1]。
由于管道线路较长,通常跨越不同地区,沿线地质地貌条件复杂,因此在建设过程中难以避免会穿越不良地质地段,在不稳定因素的影响下极有可能诱发地质灾害,其中滑坡是山区油气管道工程中常见的地质灾害类型[2]。在滑坡作用下,管道容易产生变形甚至破损,导致燃气泄漏,严重威胁沿线民众的正常生活。因此,如何对管道滑坡进行有效的支护并保护管道安全值得关注。
花管微型桩是一种新型压力灌浆式抗滑桩,集中了钢花管桩与微型桩的优点,在20世纪80年代开始在国内外应用并很快在基坑支护、滑坡治理等领域被广泛采用。目前,国内外学者对花管微型桩加固滑坡的承载性能已经开展了部分研究工作,但针对管道—滑坡体系支护性能研究方面还未涉及。潘锋[3]提出了一种注浆钢花管桩加固滑坡的理论计算公式;陈强等[4]通过离心模型试验研究发现注浆钢管桩在控制滑坡土体位移方面的有益作用;Wang等[5]针对钢花管桩微型桩的双注浆技术展开了研究,并通过试验证明了此工艺可以有效提高微型桩承载性能。
螺纹微型桩来源于地基基础工程中使用的螺纹桩,因其利用“螺丝钉比钉子更牢固”的原理[6],使螺纹桩拥有较好的承载性能,但螺纹桩目前只在地基工程中应用较多,还未被作为支挡结构应用在滑坡治理工程中。国内外对其进行的研究多集中在竖向承载性能方面,如Krasiński[7]采用数值模拟方法对螺纹桩承载机理进行了研究,并通过现场载荷试验验证了方法的可靠性;方崇等[8]通过静载试验分析了螺杆桩的竖向承载力传递特性与受力特征;叶阳升等[9]对高速铁路中的螺杆桩复合地基进行了原位测试试验,研究获得了直杆端与螺纹段在不同受力状态下的桩侧摩阻力关系;Malik等[10]进行了螺纹桩与直杆桩承载性能的对比研究,结果证明螺纹桩轴径比为2~4.1的螺旋桩端承载力比类似桩轴径的直杆桩高2~12倍;孟振等[11]通过模型试验手段针对砂土地基中的螺杆桩承载特性展开了研究。
为探究花管微型桩与螺纹微型桩加固管道—滑坡体系下桩体变形特征、桩体两侧土压力的空间分布规律以及在两种不同微型桩加固保护下管道的变形特点及破坏模式,进一步分析对比两种不同微型桩的支护性能,本文以中贵天然气管道K558+700滑坡工程为依托,以千斤顶模拟滑坡水平推力[12],开展7种不同水平加载压力工况下的室内大型推桩模型试验,并对桩顶位移、桩体应变、桩体两侧土压力、天然气管道应变以及管道两侧的土压力进行测试。进一步分析桩体弯矩以及桩体两侧峰值土压力空间分布规律,确定不同微型桩桩体变形特征,并结合管道变形破坏情况,探究两种微型桩在管道滑坡下的承载特性,总结试验成果,为我国管道地段滑坡灾害治理提供理论参考。
1. 工程概况
中贵天然气管道K558+700滑坡治理工程位于甘肃省陇南市成县黄陈镇中湾村。2020年8月12—18日中湾村附近连降暴雨,导致中贵天然气管道K558+700处斜坡发生明显滑动变形,管道滑坡区域现场见图1。共产生2处滑坡及1处滑塌体。其中,中贵天然气管道位于H1滑坡中部,横坡敷设通过,管径1016 mm,管道走向186°,滑坡区管道埋深2.2~4.1 m。H1滑坡纵长约190 m,宽约186 m,滑坡平面面积约2.25×104 m2,滑体平均厚度约8.5 m,滑体总体积约19.13×104 m2,属于中型土质滑坡。该滑坡地表变形强烈,尤其是与管道斜交的乡村水泥道路已完全损毁,受滑坡推挤作用,威胁管道安全运营。
2. 试验设计
2.1 相似设计
本次试验设计以中贵天然气管道K558+700处H1滑坡典型断面为工程原型,参考试验模型箱尺寸,确定试验模型的几何相似比尺为:
(1) 式中:CL——几何相似常数;
LP——原型尺寸;
LM——模型尺寸。
根据模型试验中各因素对现象影响的大小,抓住其起主要作用的因素,略去其次要因素的原则。以微型桩几何尺寸(L)、密度(ρ)、重力加速度(g)为主要控制参数,其相似比分别为30∶1、1∶1、1∶1,其余参数相似比可根据Buckinghamπ定理导出,如表1所示。
表 1 相似比设计Table 1. Design of similarity constant物理量 相似比 物理量 相似比 几何尺寸 CL=30 变形模量 质量密度 摩擦角 重度 Cy=1 黏聚力 应变 时间 位移 重力加速度 2.2 相似材料及配比设计
制作与原型形状完全相同的螺纹桩、花管桩较为困难,因此对其桩体模型进行简化,采用阻氧双色PP-R管材通过注水泥浆模拟桩体,其中花管桩模型为在PP-R管材长度方向上每隔8 cm环绕管道螺旋打孔并采用压力灌注水泥浆模拟;螺纹桩模型为在PP-R管材外径螺旋缠绕塑料软管并灌注水泥浆进行模拟;承台采用硬质木板模拟;天然气管道采用直径20 cm的PVC管模拟。制作完成的模型桩见图2。
滑坡岩土体相似材料较为复杂,在配比过程中不易满足推导得到的相似关系,因此试验选取滑体、滑带及基岩材料时,以最易影响原型滑坡岩土体稳定性的若干参数 (容重、黏聚力、内摩擦角、弹性模量等)为基础,参考其他学者已得到的部分研究成果[13],通过正交配比设计,并经直剪试验和三轴试验检验参数取值的准确合理性,确定使用红粉土、石英砂、水泥、石膏、水的混合物模拟Ⅳ级基岩,具体质量配比为70∶30∶5∶3∶10;滑体采用红粉土、石英砂、水的混合物模拟,具体质量配比为70∶20∶10;滑动带采用石英砂、土、滑石粉、水的混合物模拟,具体质量配比为27∶52∶35∶15。得到的模型材料相关力学参数及其与原型材料的对比见表2。
表 2 模型材料与原型材料相关物理性质参数Table 2. Physical property parameters related to model material and prototype material物理力学参数 重度/(kN·m−3) 内摩擦角/(°) 黏聚力/kPa 弹性模量/MPa 滑体 原型 19.0 30.0 40.0 / 模型 19.0 30.5 1.2 / 滑带 原型 18.5 25.0 30.0 / 模型 18.0 25.0 1.0 / 基岩 原型 27.0 / / 5000 模型 26.7 / / 160 2.3 模型设计
根据相似设计,试验使用的模型箱尺寸为 4 m×2 m×2 m(长×宽×高),箱体两侧由钢板及有机玻璃组成,箱体内填筑土体,从上往下分三层:滑体、滑带及基岩层。箱体内共布置4组微型桩群模型,其中左侧2组为花管桩,右侧2组为螺纹桩;桩后滑体内埋置有天然气管道模型。箱体后侧设置有反力墙,反力墙上安装有液压千斤顶,千斤顶通过作用在承压板上从而给滑坡模型施加水平推力。模型整体示意图如图3。
微型桩群桩模型单桩布置形式为5行×3列,共布置4组微型群桩,以模型箱中线为分界线分为左右两侧,左侧2组为花管桩,右侧2组为螺纹桩;单桩桩体模型长160 cm,直径2.5 cm。微型桩群桩埋深160 cm,从桩顶到桩底依次穿过滑体、滑带与基岩层,其中滑体层厚75 cm,滑体坡面角40°,滑带厚度5 cm,基岩层厚100 cm;天然气管道布置在桩前滑体中,距离桩体20 cm,埋深40 cm。桩顶采用承台方式将群桩连接为整体,单桩中心间距为7 cm,承台长38 cm,宽24 cm,厚2 cm,如图4所示。
2.4 传感器布设位置及监测设备
试验数据采用DH3816N静态应变测试系统采集,该系统具有60个采集通道,可同时采集应变、土压力等测试数据。
2.4.1 百分表与应变片布设位置
选取模型箱体中间两组微型桩群桩作为研究对象,在桩顶承台设置百分表测试桩顶位移,并在群桩中选取不同位置的典型单桩作为应变测试桩粘贴应变片,百分表设置位置及测试单桩的位置如图5所示。
测试单桩为2#花管桩及3#螺纹桩群桩中的第二列单桩,并对其进行编号,从桩后到桩前依次为1、2、3号测试单桩,每一测试单桩在桩体两侧沿桩体深度间隔一定距离粘贴应变片,每根测试桩体上共粘贴2×6=12个应变片,所有测试桩体共布设12×6=72个应变片,具体布设位置如图6所示(以1号测试单桩为例,其余测试桩相同)。
天然气管道山侧与河侧对称粘贴3组应变片,一组位于花管桩侧中部,一组位于螺纹桩侧中部,以监测在两种不同微型桩群桩支护作用下天然气管道的变形并作为对比,进一步优化天然气管道滑坡的支护方案设计;另一组粘贴在天然气管道接口处,以监测管道接口处的变形,进一步研究在滑坡作用下天然气管道接口处的处理方法。天然气管道应变片布设位置见图7。
2.4.2 土压力布设位置
为了研究花管微型桩与螺纹微型桩支护下桩周土压力空间分布规律,在2#群桩与3#群桩前后沿桩体深度间隔一定距离对称布置土压力盒,布设位置如图8,共布设6×2×2=24个土压力盒,以监测在滑坡水平推力作用下两种微型桩群桩前后不同深度位置土压力的变化,以进一步分析两种不同桩体变形特征。
在天然气管道前后两侧对称布设土压力盒,布设位置分别位于花管桩支护侧、螺纹桩支护侧及管道接口附近,共布设3×2=6个土压力盒,以监测天然气管道附近土体的应力变化。具体布设位置如图9。
2.5 加载工况设计
试验采用油压千斤顶施加水平荷载模拟滑坡推力,加载方法为慢速维持加载法分级加载[14],初始加载压力为1 MPa,每级加载压力增加0.5 MPa,试验一共设置7组加载工况。每级加载后观察桩顶百分表读数稳定后进行下级加载,具体加载工况见表3。
表 3 各级加载工况Q值Table 3. Q values of loading conditions at all levels工况 加载压力/MPa 工况 加载压力/MPa 1 1.0 5 3.0 2 1.5 6 3.5 3 2.0 7 4.0 4 2.5 − − 3. 试验结果分析
3.1 桩顶水平位移
由图10荷载下桩顶水平位移数据分析可知,在各级水平荷载作用下,螺纹桩侧3#与4#群桩顶水平位移明显大于花管桩侧1#与2#群桩,在4 MPa加载压力作用下螺纹桩侧位移量为花管桩侧位移量的两倍之多,说明在相同水平推力荷载下花管桩抗滑承载力更强;在2 MPa以下水平荷载作用下桩顶位移变化微弱,2~3 MPa时桩顶位移开始有小幅度增长,当荷载达到3 MPa时,桩顶位移变化量出现转折,增长幅度明显加大。根据桩顶位移容许值为30~50 mm,以1.67 mm(相当于原型的50 mm)为作为试验模型桩顶位移容许值,螺纹桩侧桩顶位移量在水平荷载为2.5 MPa达到容许值,花管桩侧桩顶位移量在水平荷载为3 MPa达到容许值,表明在本次试验条件下花管微型桩的加固效果优于螺纹微型桩。
3.2 微型桩群桩两侧土压力分析
3.2.1 桩体两侧土压力变化时程特性
根据图11可以看出,除个别测点外,在加载时段微型桩山侧土压力变化呈阶梯状上升,反映了随着加压荷载的增大,微型桩靠山侧土体应力在不断增加,并且有明显分级现象。
对比两侧不同桩体类型的土压力变化可知,花管桩侧不同埋深的土压力值普遍高于螺纹桩侧,花管桩侧在TH2测点达到最大值为46.09 kPa,相比螺纹侧在TL2测点达到的最大值20.54 kPa,花管桩侧土体最大应力高出螺纹桩侧两倍以上,结合螺纹桩桩顶最大位移与花管桩桩顶位移比值为2.4,说明花管桩侧土体受挤压程度大,分析为土拱效应[15]导致,在花管微型桩支护作用下桩后土体应力得不到有效释放,进一步表明花管桩的抗弯刚度高于螺纹桩。
由图11不同测点深度的土压力变化曲线可以看出,螺纹微型桩与花管微型桩的中上部位置土压力值最大,桩顶次之,桩底一般较小。
由图12,河侧螺纹桩测点4因土压力盒出现故障导致采集数据失真,故不作为分析对象。由图12可知,桩体河侧土压力变化趋势与山侧几乎保持一致,在加载阶段都有明显的上升趋势且各级加载土压力有明显的上升幅度加大的现象,卸载阶段开始下降并最终稳定在零值附近。
两侧桩体都在测点3(滑带附近位置)的最大级加载阶段达到峰值土压力,且螺纹桩侧峰值土压力(38.81 kPa)大于花管桩侧峰值土压力(32.57 kPa),高出比值约120%,对比山侧峰值土压力,花管桩河侧峰值土压力比山侧减少41%,螺纹桩河侧峰值土压力增加89%,由于抗滑桩承担并抵消了一部分滑坡推力,桩体河侧的土压力明显小于山侧,但花管桩两侧土压力差值明显大于螺纹桩两侧土压力差值,表明花管桩在承担滑坡推力作用方面发挥了良好的效果。
3.2.2 桩体两侧峰值土压力分析
提取桩体两侧各测点在各级加压荷载下土压力的峰值土压力进行单独分析,如图13绘制沿桩身深度不同位置随荷载改变的峰值土压力变化曲线分布图,峰值土压力基本随加压荷载的增大而增大,花管桩侧在4 MPa荷载下0.4 m桩深位置达到最大值46.09 kPa,螺纹桩侧在4 MPa荷载下0.4 m桩深位置达到最大值20.54 kPa。
花管桩侧峰值土压力沿桩深分布图大体呈“S”曲线形,且在各级荷载下分布规律保持一致,桩体中上部(位于滑带与桩顶中部)土压力值最大,在滑带附近土压力值有所减小,又在桩体中下部(位于滑带与桩底中部)有小幅度增大,往桩底位置开始减小。螺纹桩侧峰值土压力沿桩深分布图大体呈双“S”曲线形,同样在桩体中上部位置土压力值最大,滑带附近减小,桩体中下部土压力值最小,桩底土压力值开始增大。总的来说,花管桩侧土压力分布规律为桩中>桩顶>桩底,螺纹桩侧土压力分布规律为桩中>桩底>桩顶。
滑面以上的土压力由推力荷载的传递引起,滑面以下的土压力主要由微型桩群的变形引起[16]。根据上述土压力分布规律,由此说明,花管桩在滑面以上桩中上部位置所受的推力荷载最大,滑面以下桩体中下部位置变形最大;螺纹桩在滑面以上桩中上部位置所受的推力荷载最大,滑面以下桩底位置变形最大。
根据图14可知,花管桩侧在4 MPa荷载下0.65 m桩深位置达到最大值32.57 kPa,螺纹桩侧在4 MPa荷载下0.65 m桩深位置达到最大值38.81 kPa。
花管桩河侧峰值土压力分布形式与山侧分布形式大体呈现相似,土压力大小分布规律仍为桩中>桩顶>桩底,但土压力最大值出现位置下移到滑面附近,说明在花管微型桩支护作用下推力荷载在土体的分布形式发生了改变,由于花管桩刚度大于桩周土体,土体应力传递到桩体上产生应力集中,桩后土体应力重分布,最大水平应力作用点下移。螺纹桩侧河侧峰值土压力相比山侧分布形式产生了较大差异,土压力大小分布规律变为桩中>桩顶>桩底,桩底的土压力相比山侧有很大幅度减小。与花管桩侧类似,由于桩前土体应力重分布,土压力最大值出现位置下移到滑面附近且比山侧土压力最大值高出89%,说明螺纹微型桩并没有有效承担推力荷载,变形较大,在推力荷载作用桩前土体应力仍较大。
3.3 桩身弯矩分析
选取花管微型桩1号测试单桩与螺纹微型桩1号测试单桩作为研究对象,根据桩身两侧应变数据可由以下公式计算测点弯矩值[17]:
(2) 式中:M——测试截面的弯矩/(kN·m);
E——桩身材料的弹性模量/MPa;
I——测试界面的惯性矩/m4;
ε1——微型桩后的桩身应变;
ε2——微型桩前的桩身应变;
d——桩身直径/m。
图15—16为两种桩体各截面测点在分级加压荷载下的弯矩分布图,图中弯矩正值表示桩后侧受拉。
从图15可以看出,在1~3.5 MPa加压荷载下,弯矩分布特征几乎保持一致,花管桩桩身弯矩沿深度方向呈“M”形分布,桩顶与桩底弯矩值较小,桩身中上部与中下部出现较大的负弯矩,桩身离模拟滑面以上5 cm位置处产生最大正弯矩;在4 MPa加载压力作用下,弯矩分布曲线有所偏移,桩身中上部位置负弯矩极具增大,最大正弯矩位置下移到桩中位置(模拟滑面处)且桩身正弯矩长度范围增大到0.5 m,滑面以下负弯矩最大值减小,且在桩底出现正弯矩。说明在1~3.5 MPa推力荷载下,花管桩桩体变形特征一致,桩后桩身中上部与中下部位置受压变形,桩后桩身离模拟滑面0~10 cm段受拉变形;在4 MPa推力荷载下,桩身弯曲变形加剧,在桩深0.4 m 处出现最大负弯矩−71.96 kN·m,在桩深0.8 m处(模拟滑面位置)出现最大正弯矩16.47 kN·m。花管桩设计时对桩身变形较大处可开展进一步的优化工作。
由图16可知,螺纹桩桩身弯矩分布沿深度方向呈“S”形,桩体正负弯矩位置在模拟滑面附近大致呈旋转对称分布,滑面以上大部分区段为负弯矩,滑面以下为正弯矩;螺纹桩在桩深0.4处出现最大负弯矩−64.44 kN·m,在桩深0.9m处(模拟滑面以下10 cm)出现最大正弯矩65.85 kN·m,相比花管桩,最大负弯矩相差不大,最大正弯矩高出近4倍,表明在相同推力荷载工况下,螺纹微型桩变形程度大于花管微型桩。
3.4 管道受力变形特征分析
3.4.1 管道两侧土压力变化规律及分布特性
为了研究天然气管道在两种不同微型桩支护作用下的变形特征,通过埋设土压力盒监测管道两侧土体应力变化,图17为管道山侧与河侧各测点的土压力变化时程曲线。其中TG1/TG1'测点位于花管桩支护侧,TG3/TG3'测点位于螺纹桩支护侧,TG2/TG2'测点位于管道接口位置两侧。
根据图17可以看出,管道山侧与河侧各测点土压力在加载阶段都有明显的随时间增大的趋势。山侧各测点变化趋势大致相同,在初期加载土压力值增长表现平稳,而在后期加载阶段土压力出现跳跃式增长,并在最后加载阶段达到峰值,位于螺纹桩支护侧的TG3测点峰值土压力达22.99 kPa,为花管桩支护侧的1.73倍;结合图18各级荷载下管道山侧峰值土压力分布曲线可知,管道山侧峰值土压力随荷载的增大而增大,并且螺纹桩支护侧峰值土压力整体大于花管桩支护侧及管中接口位置。说明花管桩在推力荷载下承载力表现优于螺纹桩,减小了传递到管道的滑坡推力,而螺纹桩在较大推力荷载下由于抗弯性能不足,变形严重,破坏后不能有效承担并抵消滑坡推力,传递到桩前管道的应力较大,从而导致管道易产生变形破坏。
由图17,管道河侧土压力变化趋势与山侧存在较大差异,在加载阶段整体呈波动性增长,测点最大峰值土压力在螺纹桩支护侧为12.29 kPa,最小峰值土压力在花管桩支护侧为1.54 kPa,比值7.98∶1。结合图19管道河侧峰值土压力分布曲线,在分级推力荷载作用下,峰值土压力规律表现为:花管桩支护侧<管中接口位置<螺纹桩支护侧。
3.4.2 管道弯矩分析
由图20分级荷载管道弯矩分布图可以看出,在较小推力荷载作用下(1~3 MPa),管道接口附近弯矩值整体较小,花管桩支护侧呈现负弯矩,螺纹桩支护侧呈现正弯矩,绝对值大小相差不大,说明在较小推力荷载下,管道变形程度较小;在较大推力荷载作用下(3.5~4 MPa),管道接口处弯矩急剧增大,最大值达−317.75 kN·m,说明在较大推力荷载下管道接口位置变形严重甚至破坏,在设计滑坡区管道工程时应注意管道接口位置的加固;另一方面,对比花管桩支护侧与螺纹桩支护侧管道弯矩可知,弯矩分布曲线明显向螺纹桩侧偏移,并整体大于花管桩支护侧,说明在本试验条件下花管微型桩支护作用下天然气管道的变形程度低于螺纹微型桩支护,在滑坡推力作用下 ,花管微型桩对管道的保护效益更突出,花管桩更适用于作为管道滑坡区域的支挡结构。
4. 结论
(1)花管微型桩山侧及河侧峰值土压力沿桩深分布形式基本相似,大体呈“S”曲线形,桩后土体土拱效应明显,且在各级荷载下分布形式大致保持一致,总体来说花管桩侧土压力分布规律为桩中>桩顶>桩底,滑带附近的桩体周围土压力较大,在抗滑桩设计工作中应重点考虑优化。
(2)螺纹桩山侧峰值土压力沿桩深分布图大体呈双“S”曲线形,河侧峰值土压力相比山侧分布形式产生了较大差异,桩底的土压力相比山侧减小很大幅度;随外部荷载的增加桩周土压力增加幅度较大,表明螺纹微型桩在横向承载性能方面有所欠缺。
(3)花管桩桩身弯矩沿深度方向呈“M”形分布,桩身离模拟滑面以上5 cm位置处产生最大正弯矩;螺纹桩桩身弯矩分布沿深度方向呈“S”形,桩体正负弯矩位置在模拟滑面附近大致呈旋转对称分布,滑面以上大部分区段为负弯矩,滑面以下为正弯矩;在相同推力荷载工况下,螺纹微型桩变形程度大于花管微型桩。
(4)在滑坡作用下花管微型桩可以有效减小传递到管道的坡体应力,在一定程度上预防管道受力破坏;而螺纹桩在较大推力荷载下抗弯性能不足,变形严重,破坏后不能有效承担滑坡推力,传递到桩前管道的应力较大,从而导致管道变形程度更为强烈;在本试验条件下,花管微型桩对管道的保护效益突出,更适用于作为管道-滑坡区域的支挡结构。
-
表 1 回填层岩土体参数先验信息
Table 1 Prior statistical characteristics of soil parameters for the backfill layer
土体类型 土体参数 μ σ COV 概率分布 回填层 饱和渗透系数/(m·h−1) 0.108 0.0648 0.6 对数正态分布 黏聚力/kPa 10 3 0.3 对数正态分布 内摩擦角/(°) 15 3 0.2 对数正态分布 饱和体积含水量 0.315 — — — 饱和重度/(kN·m−3) 20 — — — -
[1] 文海家, 张岩岩, 付红梅, 等. 降雨型滑坡失稳机理及稳定性评价方法研究进展[J]. 中国公路学报,2018,31(2):15 − 29. [WEN Haijia, ZHANG Yanyan, FU Hongmei, et al. Research status of instability mechanism of rainfall-induced landslide and stability evaluation methods[J]. China Journal of Highway and Transport,2018,31(2):15 − 29. (in Chinese with English abstract)] WEN Haijia, ZHANG Yanyan, FU Hongmei, et al . Research status of instability mechanism of rainfall-induced landslide and stability evaluation methods[J]. China Journal of Highway and Transport,2018 ,31 (2 ):15 −29 . (in Chinese with English abstract)[2] 胡康,任光明,常文娟,等. 基于节理不确定性的可靠度分析——以西藏某岩质边坡为例[J]. 中国地质灾害与防治学报,2022,33(2):53 − 60. [HU Kang, REN Guangming, CHANG Wenjuan, et al. Reliability analysis based on joint uncertainty:A case study of a rock slope in Tibet[J]. The Chinese Journal of Geological Hazard and Control,2022,33(2):53 − 60. (in Chinese with English abstract)] HU Kang, REN Guangming, CHANG Wenjuan, et al . Reliability analysis based on joint uncertainty: A case study of a rock slope in Tibet[J]. The Chinese Journal of Geological Hazard and Control,2022 ,33 (2 ):53 −60 . (in Chinese with English abstract)[3] 陈忠源, 戴自航, 简文彬. 基于因子权重反分析的新近失稳土质边坡稳定性评价云模型[J]. 中国地质灾害与防治学报,2023,34(4):125 − 133. [CHEN Zhongyuan, DAI Zihang, JIAN Wenbin. Cloud model for stability evaluation of recently failed soil slopes based on weight inversion of influencing factors[J]. The Chinese Journal of Geological Hazard and Control,2023,34(4):125 − 133. (in Chinese with English abstract)] CHEN Zhongyuan, DAI Zihang, JIAN Wenbin . Cloud model for stability evaluation of recently failed soil slopes based on weight inversion of influencing factors[J]. The Chinese Journal of Geological Hazard and Control,2023 ,34 (4 ):125 −133 . (in Chinese with English abstract)[4] 王继玲, 周维博, 孙梨梨, 等. 石川河富平地下水库渗透系数空间变异性研究[J]. 水文地质工程地质,2023,50(3):34 − 43. [WANG Jiling, ZHOU Weibo, SUN Lili, et al. Study on the spatial vriability of hydraulic conductivity of underground reservoir in Fuping section of Shichuan River[J]. Hydrogeology & Engineering Geology,2023,50(3):34 − 43. (in Chinese with English abstract)] WANG Jiling, ZHOU Weibo, SUN Lili, et al . Study on the spatial vriability of hydraulic conductivity of underground reservoir in Fuping section of Shichuan River[J]. Hydrogeology & Engineering Geology,2023 ,50 (3 ):34 −43 . (in Chinese with English abstract)[5] YANG Haoqing,ZHANG Lulu,PAN Qiujing,et al. Bayesian estimation of spatially varying soil parameters with spatiotemporal monitoring data[J]. Acta Geotechnica,2021,16(1):263 − 278. DOI: 10.1007/s11440-020-00991-z
[6] CAO Zijun,WANG Yu,LI Dianqing. Site-specific characterization of soil properties using multiple measurements from different test procedures at different locations–A Bayesian sequential updating approach[J]. Engineering Geology,2016,211:150 − 161. DOI: 10.1016/j.enggeo.2016.06.021
[7] DIAZDELAO F A,GARBUNO-INIGO A,AU S K,et al. Bayesian updating and model class selection with Subset Simulation[J]. Computer Methods in Applied Mechanics and Engineering,2017,317:1102 − 1121. DOI: 10.1016/j.cma.2017.01.006
[8] STRAUB D,PAPAIOANNOU I. Bayesian updating with structural reliability methods[J]. Journal of Engineering Mechanics,2015,141(3):04014134. DOI: 10.1061/(ASCE)EM.1943-7889.0000839
[9] 刘贤, 揭鸿鹄, 蒋水华, 等. 融合历史降雨下斜坡稳定性观测信息的可靠度分析[J]. 地球科学,2023,48(5):1865 − 1874. [LIU Xian, JIE Honghu, JIANG Shuihua, et al. Slope reliability analysis incorporating observation of stability performance under a past rainfall event[J]. Earth Science,2023,48(5):1865 − 1874. (in Chinese with English abstract)] LIU Xian, JIE Honghu, JIANG Shuihua, et al . Slope reliability analysis incorporating observation of stability performance under a past rainfall event[J]. Earth Science,2023 ,48 (5 ):1865 −1874 . (in Chinese with English abstract)[10] 曹子君, 胡超, 苗聪, 等. 基于分层贝叶斯学习的滨海软土地层高效识别方法[J]. 地球科学,2023,48(5):1730 − 1741. [CAO Zijun, HU Chao, MIAO Cong, et al. Efficient identification method of coastal soft soil stratum based on hierarchical Bayesian learning[J]. Earth Science,2023,48(5):1730 − 1741. (in Chinese with English abstract)] CAO Zijun, HU Chao, MIAO Cong, et al . Efficient identification method of coastal soft soil stratum based on hierarchical Bayesian learning[J]. Earth Science,2023 ,48 (5 ):1730 −1741 . (in Chinese with English abstract)[11] DEPINA I,OGUZ E A,THAKUR V. Novel Bayesian framework for calibration of spatially distributed physical-based landslide prediction models[J]. Computers and Geotechnics,2020,125:103660. DOI: 10.1016/j.compgeo.2020.103660
[12] PAPAIOANNOU I,BETZ W,ZWIRGLMAIER K,et al. MCMC algorithms for subset simulation[J]. Probabilistic Engineering Mechanics,2015,41:89 − 103. DOI: 10.1016/j.probengmech.2015.06.006
[13] ZHANG J,ZHANG L M,TANG W H. Slope reliability analysis considering site-specific performance information[J]. Journal of Geotechnical and Geoenvironmental Engineering,2011,137(3):227 − 238. DOI: 10.1061/(ASCE)GT.1943-5606.0000422
[14] AU S K,BECK J L. Estimation of small failure probabilities in high dimensions by subset simulation[J]. Probabilistic Engineering Mechanics,2001,16(4):263 − 277. DOI: 10.1016/S0266-8920(01)00019-4
[15] WANG C H,FANG Li,CHANG D T T,et al. Back-analysis of a rainfall-induced landslide case history using deterministic and random limit equilibrium methods[J]. Engineering Geology,2023,317:107055. DOI: 10.1016/j.enggeo.2023.107055
[16] 蒋水华, 李典庆, 周创兵, 等. 考虑参数空间变异性的非饱和土坡可靠度分析[J]. 岩土力学,2014,35(9):2569 − 2578. [JIANG Shuihua, LI Dianqing, ZHOU Chuangbing, et al. Reliability analysis of unsaturated slope considering spatial variability[J]. Rock and Soil Mechanics,2014,35(9):2569 − 2578. (in Chinese with English abstract)] JIANG Shuihua, LI Dianqing, ZHOU Chuangbing, et al . Reliability analysis of unsaturated slope considering spatial variability[J]. Rock and Soil Mechanics,2014 ,35 (9 ):2569 −2578 . (in Chinese with English abstract)[17] FREDLUND D G,XING Anqing,HUANG Shangyan. Predicting the permeability function for unsaturated soils using the soil-water characteristic curve[J]. Canadian Geotechnical Journal,1994,31(4):533 − 546. DOI: 10.1139/t94-062
[18] YANG Haoqing,ZHANG Lulu,XUE Jianfeng,et al. Unsaturated soil slope characterization with Karhunen–Loève and polynomial chaos via Bayesian approach[J]. Engineering with Computers,2019,35(1):337 − 350. DOI: 10.1007/s00366-018-0610-x
[19] HUANG J,GRIFFITHS D V. Determining an appropriate finite element size for modelling the strength of undrained random soils[J]. Computers and Geotechnics,2015,69:506 − 513. DOI: 10.1016/j.compgeo.2015.06.020
[20] LIU Xin,WANG Yu,KOO R C H,et al. Development of a slope digital twin for predicting temporal variation of rainfall-induced slope instability using past slope performance records and monitoring data[J]. Engineering Geology,2022,308:106825. DOI: 10.1016/j.enggeo.2022.106825
[21] ZHANG J,ZHANG L M,TANG W H. Bayesian framework for characterizing geotechnical model uncertainty[J]. Journal of Geotechnical and Geoenvironmental Engineering,2009,135(7):932 − 940. DOI: 10.1061/(ASCE)GT.1943-5606.0000018
-
期刊类型引用(2)
1. 杨晓华,王东清,袁帅,张毅. 基于均匀化理论与上限分析的膨胀土滑坡稳定性分析. 水文地质工程地质. 2024(02): 172-182 . 百度学术
2. 刘苗. 影响某天然气管道滑坡发育特征及变形监测分析. 中国石油和化工标准与质量. 2024(09): 152-154 . 百度学术
其他类型引用(1)