1、第 卷 第 期 年 月人 民 长 江 ，收稿日期：基金项目：国家管网集团西南管道有限责任公司科技研发项目（）作者简介：葛 华，男，高级工程师，主要从事管道工程设计工作。：通信作者：李 虎，男，助理工程师，硕士，主要从事管道完整性研究工作。：文章编号：（）引用本文：葛华，方迎潮，朱建平，等 穿越管道滑坡地表位移预警阈值及预警模型研究 人民长江，（）：穿越管道滑坡地表位移预警阈值及预警模型研究葛华，方 迎 潮，朱 建 平，李虎，王庆，时 建 辰（国家管网集团西南管道有限责任公司，四川 成都；四川省地质工程勘察院集团有限公司，四川 成都）摘要：为了解决管道穿越滑坡区地表位移预警阈值难以确定的问题，运
2、用弹性地基梁模型和滑坡位移控制方法推导了滑坡区管道位移和应力力学模型，获取了管道能承受的最大位移；通过管道最大位移和地表位移的耦合关系，定量化确定了滑坡地表位移预警阈值计算方法；结合地下水预警阈值与管道应力阈值，构建了管道滑坡多因素耦合预警模型；并将其运用到贵州 管道滑坡灾害。结果表明：滑坡灾害地表位移蓝、黄、红色预警阈值分别为 ，地下水各级阈值为 ，。现场监测数据显示，地表位移达到 ，触发蓝色预警，地下水位高度为 ，达到黄色预警，综合预警级别达到关注级。根据地表位移动态预警、地下水预警与管道应力预警构建的管道滑坡灾害动态耦合预警模型，能在一定程度上为管道滑坡灾害提供合理的预警，保障管道安全运
3、行，提高预警成功率。关 键 词：地表位移；预警阈值；预警模型；力学模型；滑坡；管道中图法分类号：文献标志码：引 言管道作为油气的主要运输手段，承担着中国 的原油和 的天然气运输，助力着中国的经济发展。由于中国地形地貌特点，长输油气管道不可避免地穿越滑坡区，其中尤以西南管道公司所辖管道为甚。资料显示，西南管道公司所辖管道中，山区管道占比超过，管道沿线地质环境条件复杂，其中穿越大量的滑坡区，滑坡灾害发育引发管道周围土体发生变形，从而会导致埋地管道产生弯曲、压缩、扭曲、拉裂、局部屈曲等破坏行为。国内外计算滑坡区管道位移力学模型常将滑坡作为力荷载形式并作为主动力，在此基础上，梁政运用弹性地基梁理论对管
4、道横穿滑坡时的受力行为进行了推导；邓道明等首次得出拉、压当量轴力下管道的内力及变形的计算公式；郝建斌等假设管后土体为刚性楔形体，并对横穿滑坡时管道所受的推力进行了推导。许多学者根据滑坡位移塑性变形、双参数模型、土体位移四次分布假设等对滑坡管道的位移进行了理论推导。但是，以上滑坡管道位移力学模型的研究多以力荷载作为主动力，不便于体现管道极限位移与土体位移之间的关系，对滑坡监测工程有着一定的局限性。最早对滑坡灾害预警预报的方法是在 世纪 年代，包括通过观察地表变形时间、地下水变化异常、地表裂缝的变化和动物的异常表现等状况判断滑坡是否即将失稳的现象预报法，以及斋藤法、位移时间曲线变化判断法、统计数学
5、模型法和参数预报法等位移预报模型。不少学者针对管道与滑坡的相互关系提出了不同的监测预警模型，日本学者 第 期 葛 华，等：穿越管道滑坡地表位移预警阈值及预警模型研究利用分布式滑坡概念模型估计强降雨条件对浅层滑坡的影响；刘人杰以川气东送普光 宜昌段为例，开展了滑坡与管道变形的监测预警研究；俞政运用灰色模型对忠武输气管道滑坡段进行了预警研究；李岩岩以西南某天然气管道为例建立了滑坡管道预警模型；冷建成等根据管道应力监测数据，采用 神经网络作为管道应力变化预警模型；陈星明通过监测管道的变形确定预警等级；许强等提出了一种改进的切线角及对应的滑坡预警判据，解决了位移与时间量纲不统一的问题。但以上研究均未对
6、管道滑坡灾害的地表位移和地下水预警阈值指出定量化阈值计算方法，也未形成多因素下的预警模型。本文将根据弹性地基梁模型，以弹性地基梁模型推导管道位移力学模型，结合管道最大位移和地表位移的耦合关系，以及地下水与滑坡稳定性的关系，确定地表位移和地下水预警阈值计算方法，并与管道应力监测预警相结合建立管道滑坡灾害多因素耦合预警模型。横穿滑坡管道极限位移计算众多学者针对管道极限位移计算多采用均布力荷载、二次型力荷载或楔形体力荷载形式计算管道的变形特征，在实际监测工程中，以力的形式计算得到的管道变形特征并不能完全表示管道屈服时能达到的极限位移值，不便于监测预警。以弹性地基梁模型为基础推导管道力学方程，省略了滑
7、坡推力的计算，运用管道位移与应力之间的相互关系进行推导，能为监测预警提供管道极限位移值，为预警阈值提供基础。因此，本次研究运用半无限连续弹性地基梁模型推导管道在横穿滑坡作用时能承受的极限位移。非滑坡区管道力学行为由于非滑坡区管道为小变形，且管道受到力始终处于平衡状态，可将此段看作是半无限连续弹性地基梁 模型，如图 所示，则此段的控制微分方程如式（）所列。图 管道穿越滑坡示意 （）式中：为管道的弯曲刚度；为管道一点距 点的距离，；为管道挠度，；为弹性地基系数。令 ，且有（），（），可解得 （）。设 点的弯矩为，则可得到方程的解为 （）（）则根据挠度、转角、弯矩之间的相互关系可以求得在 点处的挠度
8、（）、转角（）、弯矩（）和剪力（）如下：（）（）（）（）（）横穿滑坡区管道位移当管道横穿滑坡时，可将管道变形看作是对称变形，邓道明等将管道简化为梁，假定滑坡区管道只受均匀滑坡推力，因此，本次研究在前人的理论基础上，运用弹性地基梁模型推导管道力学模型，将管道受到的滑坡作用简化为受上部土体的均布荷载，可将管道看作是受均布荷载的弹性梁模型。因挠度的四阶导数等于横向荷载，则将 段设为 （）由于管道的连续性，则在 点处的力学行为与非滑坡区管道力学状态相同，又因为滑坡主滑面 点处其转角为，且滑坡不断蠕变过程中管道上 点的位移逐渐变大，点的最大位移即为管道跨中极限位移值，则边界条件为（）（）（）式中：为跨中
9、 点位移，；为 点到 点的水平距离，。联立式（）（）可解得 段挠度方程和在 点的弯矩为 ()（）（）（）由上一节可知，段与 段对称，则整段的挠度方程可由分段函数表示为 人 民 长 江 年()，()（）（）（）（），（）根据应变关于挠度的关系 ，假设管道在弹性范围内，则管道满足胡克定律，由于管道一点应力可看作由该点处的位移控制，则管线由弯曲引起的轴向应力为()，()（）（），（）由于管道在正常运行内压下工作，忽略温差的影响，因此管道总轴向应力 为弯曲产生的轴向应力 与内压产生的轴向应力 之和。（）式中：为管道泊松比；为内压，；为管道外径，；为管道壁厚，。因此，通过式（）可确定管道弯曲产生的最大应
10、力，结合式（）和（）通过迭代算法，即可求出当总轴向应力达到屈服时的最大位移量，并可通过不同 值计算得到滑坡位移分布。滑坡地表位移预警阈值根据以上滑坡动态变化过程，可知滑坡地表位移是造成管道破坏的最主要因素，针对缓慢型变形滑坡，将管道应力与滑坡形变相结合，运用式（），（），（）可计算得出管道能承受的最大位移值，并进行简化，以最不利方式考虑，即管道与滑坡看作是协同变形，管道的变形位移值就等于土体变形位移值。将管道挠度达到 的最大挠度作为滑坡形变指标红色级预警阈值，即 红，此红色级预警表示管道在此时刻已经受力非常大，极有可能在极短的时间内发生屈服，达到管道破坏。因此，此红色级预警时间点应是管道最危险
11、的时刻，维护工作需在此阶段之前结束，若未结束则维修人员需立刻退场，对管道进行关阀等相关措施，保障人员与管道安全。为了安全考虑，以管 土协同变形为最不利状态，则滑坡形变位移预警的蓝色级预警以滑坡开始运动为指标，以 地表位移监测设备为例，设备的精度为，且考虑到监测受干扰等异常情况，监测数据误差符合 原则（为设备精度），则取滑坡地表位移监测数据达到 作为蓝色预警阈值。蓝色预警时应时刻关注灾害体变形情况，做好应急处理准备措施，以及在此刻考虑灾害的危险性等。利用地表位移监测数据可以得出一段时间内的地表变形速度，以 月以内的平均形变速度确定，若监测工程较为紧急，则可适当降低变形速度获取周期。以管道挠度达到
12、红色级预警时的挠度（）与地表变形速度 相除，得到滑坡以此速度进行发育时管道能达到破坏的时间 红 ，按照抢修时间和人员调动总时间，反算出黄色预警时间点 黄，如下式：黄 （）根据式（）计算得出的黄色预警时间，结合地表位移、时间 和变化速度 的关系，即可得到黄色预警时间下地表位移预警阈值，如式（）所示。黄色预警时则应是灾害点处理开始时间点，需要预留时间对灾害点进行处理和维护，防止灾害发生，减小灾害发生所产生的损失。黄 黄（）（）因此，以监测数据达到 作为地表位移蓝色预警阈值 蓝，以式（）计算得到的 黄作为地表位移黄色预警阈值，以管道挠度达到 的最大挠度 红作为滑坡形变指标红色级预警阈值，如表 所列。
13、表 滑坡地表位移预警级别划分标准 形变指标地表位移蓝色级蓝 黄黄色级黄 红红色级红 滑坡地下水位预警阈值由滑坡动态变化过程可知，地下水与滑坡稳定性有着密切的联系。在降雨条件下，外界环境发生改变，雨水入渗到斜坡体内部，改变了其原有的地下水深度及岩土体物理力学状态，滑坡体地下水位变高，降低了滑面的强度指标，从而也使抗滑力减小，当稳定系数 突破临界值 时，斜坡体极有可能产生破坏，出现变形，进而可能演变为滑坡灾害。有许多研究都对地下水对滑坡稳定性的关系进行了分析。第 期 葛 华，等：穿越管道滑坡地表位移预警阈值及预警模型研究研究的主要问题是需要解决抗滑力与地下水深度之间的关系。在地震作用下，降雨型滑坡
14、稳定性计算模型如图 所示。图 降雨型斜坡地下水计算模型 考虑地震作用和降雨条件，将潜在滑坡体简化为二维平面问题。由于雨水下渗作用，潜水面上升，潜水面高度变大使得滑坡稳定性降低，进而可能导致滑坡产生。假设斜坡堆积体的厚度为，地下水位高度为，斜坡的滑面倾角为，近似取坡度值。根据莫尔 库伦定律以及极限平衡理论计算斜坡体的稳定系数，计算公式如下：（）（）式中：为单位斜坡体的抗滑力，；为单位斜坡下滑力，；为斜坡体的有效黏聚力，；，分别为坡体物质和水的密度大小，；为重力加速度值，；为地震加速度，；表示斜坡体的铅直厚度值，；表示水位；为坡体有效内摩擦角；为斜坡的倾角。根据稳定性系数 的计算公式可知，斜坡的稳
15、定性 与水位 呈负相关关系，当 不断增大时，逐渐减小。运用式（）反算滑坡体临界水位，得到稳定性系数与滑坡水位关系式（）：（）（）依据 滑坡防治工程勘查规范知稳定性系数 为不稳定状态；为欠稳定状态；为基本稳定状态；为稳定状态。根据不同稳定性系数代入式（）计算得出滑坡地下水位，以 ，和 计算得到的地下水位高度，和 分别作为蓝色、黄色和红色级预警阈值，得到以地下水位为监测预警的各级预警阈值，如表 所列。表 地下水位预警级别划分标准 形变指标稳定性系数水位深度蓝色级 黄色级 红色级 多因素耦合预警模型 滑坡耦合预警模型将滑坡地表位移和地下水高度相结合构建滑坡本体二维预警模型，以滑坡地表位移预警为主，以
16、地下水位预警为辅，得到滑坡耦合预警模型，如表 所列。表 滑坡地表位移与地下水位的二维预警模型 地下水位不同地表位移对应的预警级别蓝 黄黄 红红 蓝色级黄色级红色级 蓝色级黄色级红色级 黄色级红色级红色级 滑坡与管道耦合预警模型管道穿越滑坡区时，其主要关注对象为管道是否达到破坏，需要对管道本体应力进行预警，以便了解管道应力的动态变化。参照相关规范，将管道应力 预警分级规定为管道附加应力的，作为蓝色预警、黄色预警和红色预警阈值，即管道应力预警级别划分如表 所列。表 管道应力预警级别划分 应力指标管道允许附加应力蓝色级 黄色级 红色级 将滑坡本体预警与管道应力预警相互耦合，以直接反映管道动态变化的应力预警为主，以诱发管道破坏的滑坡本体预警为辅，建立管道滑坡灾害动态耦合预警模型，如表 所列。若现场监测方式仅有一种或两种，则根据相应的设备进行预警阈值设置。表 管道滑坡灾害二维预警模型 管道应力预警级不同滑坡本体预警级别对应的动态耦合预警级别蓝色级黄色级红色级蓝色级关注级警示级警报级黄色级关注级警示级警报级红色级警示级警报级警报级 此模型可广泛应用于管道穿越滑坡灾害的预警当 人 民 长 江 年中，