• ISSN 1673-5722
  • CN 11-5429/P

软弱地层综合管廊三维非线性地震反应特性研究

潘强 赵丁凤 曹泽伟 钱健 周正龙

郑舒元,王建雄,崔建文,2022. 中缅油气管道跨断裂地震变形的GNSS观测. 震灾防御技术,17(4):751−763. doi:10.11899/zzfy20220417. doi: 10.11899/zzfy20220417
引用本文: 潘强,赵丁凤,曹泽伟,钱健,周正龙,2023. 软弱地层综合管廊三维非线性地震反应特性研究. 震灾防御技术,18(1):44−52. doi:10.11899/zzfy20230106. doi: 10.11899/zzfy20230106
Zheng Shuyuan, Wang Jianxiong, Cui Jianwen. GNSS Monitoring of Cross Fault Seismic Deformation of China Myanmar Oil and Gas Pipeline Project[J]. Technology for Earthquake Disaster Prevention, 2022, 17(4): 751-763. doi: 10.11899/zzfy20220417
Citation: Pan Qiang, Zhao Dingfeng, Cao Zewei, Qian Jian, Zhou Zhenglong. Three-dimensional Nonlinear Seismic Response Characteristics of Utility Tunnel under Soft-weak Strata[J]. Technology for Earthquake Disaster Prevention, 2023, 18(1): 44-52. doi: 10.11899/zzfy20230106

软弱地层综合管廊三维非线性地震反应特性研究

doi: 10.11899/zzfy20230106
基金项目: 国家自然科学基金(52208349)
详细信息
    作者简介:

    潘强,男,生于 1990 年。工程师。主要从事地基处理与地震工程研究。E-mail:354992238@qq.com

    通讯作者:

    周正龙,男,生于 1988 年。博士,高级工程师。主要从事地震工程研究。E-mail:zhouzhen.glong@163.com

Three-dimensional Nonlinear Seismic Response Characteristics of Utility Tunnel under Soft-weak Strata

  • 摘要: 软弱地层综合管廊在强震作用下易发生破坏。考虑土体非线性及水的影响,以苏通GIL综合管廊工程为依托,采用三维线性梁单元模拟螺栓与钢筋,对地震动进行幅值标定,建立三维精细化有限元模型,并根据不同地震动输入方式及强度,从衬砌应力分布、张开量及结构损伤角度分析综合管廊非线性地震反应特性。研究结果表明,45°共轭方向与拱腰位置处管廊受力较大,横向、纵向地震动作用下,随着地震动强度的增加,应力增加明显;横向地震动输入对管廊环间张开量的影响较小,当地震动达到峰值附近时,张开量增长明显,并在一定范围内波动;横向地震动对管廊结构损伤的影响较大,峰值加速度达0.3 g时,管廊内部结构在两端连接处及中间支撑处连接点出现拉压变形。
  • 中缅油气管道工程是我国重大能源战略项目,全长1 726.8 km,其包括原油管道与天然气管道。两管道起自缅甸西海岸皎漂,经云南省瑞丽市58号界碑入境,在贵州安顺市分离,天然气管道向南经贵阳市到达贵港市。管道沿线地质背景复杂,地震活动强烈,穿越喜马拉雅地震带、滇西南地震带、鲜水河滇东地震带、右江地震带和长江中游地震带,历史上曾发生过8级地震1次、7~7.9级地震16次、6~6.9级地震65次,对管道站址影响烈度最大达8度以上。在云南境内,有3条曾发生过7级和7级以上地震的断裂与管道相交,分别为龙陵-瑞丽断裂、凤仪-定西岭断裂、小江断裂,在管线近场有18条活动断裂通过,这些断裂具备发生强烈地表错位的地震构造条件(高涵等,2020),存在未来发生强震的可能,对管道安全营运造成潜在威胁。

    强烈地震对地下管道有破坏作用(Wang等,1991),这种作用主要有2种表现形式,即强烈振动引起的管道连接脱落和地面变形引起的管体变形破坏。1976年河北唐山大地震的地震动和地表破裂对秦皇岛至北京的输油管道造成了多处变形、管体破损(王亮亮等,2003);1985年墨西哥地震的地震动造成了管线管体破碎与接头脱落;2002年美国阿拉斯加发生7.9级地震,管道穿越的断层破裂长达1 m以上,但由于采取了特殊的抗震措施,避免了管道破裂(Hall等,2003);2016年日本福冈7.3级地震中,地下水管由于地震破裂而变形损坏,泄露的水损毁了地面道路(孟晨,2018)。为减轻地震对管道破坏造成的影响,多破坏性地震地区油气管道工程重视管道的抗震措施,如在阿拉斯加油气管道穿越断层时,在管道下方采用了滑轮设计,为管道横向和纵向移动保留了空间,避免了在2002年7.9级地震中受到损坏(郭守德等,2019)。此外,布设管道地震观测系统,在破坏性地震发生后能够快速判断管道状况,及时采取措施,这也是有效减轻油气管道地震灾害的方法。在中缅油气管道沿线建立地震观测系统,当管道沿线一定范围内发生破坏性地震时,对于观测管道遭受的地震动强度及管道穿越的断裂变形情况具有重要意义。

    中缅油气管道地震观测系统由3套子系统构成,即观测管道地震动强度的强震动观测子系统、在管道与3条断层相交处观测断裂变形的GNSS(Global Navigation Satellite System)断裂变形观测子系统和观测管道变形的光纤管道变形观测子系统。本文主要介绍GNSS断裂变形观测子系统。

    GNSS是高精度全球卫星定位导航系统,已大量应用于建筑物沉降、区域性地貌变化、地壳变形观测中(帅向华等,2011李希亮等,2014刘昌伟等,2019李兆隆等,2019)。采用相对定位技术,GNSS可获得毫米级的定位精度,但需在观测点一定距离外布设参考站,而精密单点定位(Precise Point Positioning,PPP)无须布设参考站,数据处理更方便(张小红等,2015邹冠华,2018)。中缅油气管道断裂变形观测子系统在管道穿越的断裂两侧各布设1台GNSS接收机,采用PPP技术精确定位同一时刻断裂两侧观测点的位置,即可确定断层两侧的相对位错,再结合油气管道容许变形,可估计穿越断层管道可能的变形情况。利用已建成的2条断层GNSS断裂变形观测子系统获取的观测数据,本文研究了跨断层GNSS观测的数据质量及观测误差,结果表明,本项目建立的GNSS断裂变形观测子系统能够有效观测破坏性地震引起的断层变形量。

    精密单点定位算法采用双频相位和伪距观测数据,利用国际GPS服务组织(IGS)或其他分析中心发布的精密星历和钟差同时解算观测站坐标、钟差和对流层延迟及模糊度参数,PPP可利用单台双频GPS接收机的伪距和载波相位观测值进行精密绝对定位。

    对于双频GNSS接收机,PPP的伪距和载波相位消电离层组合观测方程可表示为:

    $$ P=\rho +c\left({\delta t}_{{\rm{r}}}-{\delta t}^{{\rm{s}}}\right)+{M}_{{\rm{W}}}{\xi }_{{\rm{T}}}+{\varepsilon }_{{\rm{P}}} $$ (1)
    $$ \varphi =\rho +c\left({\delta t}_{{\rm{r}}}-{\delta t}^{{\rm{s}}}\right)+\lambda N+{M}_{{\rm{W}}}{\xi }_{{\rm{T}}}+{\varepsilon }_{{\text{φ}} } $$ (2)

    式中, $ P、\varphi $ 分别为消电离层的伪距和载波相位; $ \rho $ 为卫星至接收机的几何距离; $ {\delta t}_{{\rm{r}}}\mathrm{和}{\delta t}^{{\rm{s}}} $ 分别为接收机和卫星的钟误差,后者可利用精密卫星钟差产品消除; $ {M}_{{\rm{W}}} $ 为对流层湿分量投影函数; $ {\xi }_{{\rm{T}}} $ 为信号传播方向的天顶对流层延迟; $ {\varepsilon }_{{\rm{P}}} $ $ {\varepsilon }_{\text{φ}} $ 为测量噪声;λ为消电离层组合观测值的波长;为消电离层组合观测值的浮点模糊度。

    基于上述2个观测方程,综合考虑卫星端、接收机端及信号传输过程中有关误差的影响,使用各种改正模型对定位过程改正发展的数据解算方法,可在4颗以上卫星接收数据的情况下,获取观测的位置信息,其静态定位精度可达厘米级,可满足工程观测的需要。

    在断裂两盘各布设1个GNSS观测点,利用PPP技术确定同一时刻2个观测点的相对位置,在断裂发生位错的情况下,由2个观测点相对位置变化可确定断裂的位错量:

    $$ D=L-{L}_{0} $$ (3)
    $$ {L}_{0}=\sqrt{{{(x}_{0}^{A}-{x}_{0}^{B})}^{2}+{{(y}_{0}^{A}-{y}_{0}^{B})}^{2}+{{(z}_{0}^{A}-{z}_{0}^{B})}^{2}} $$
    $$ L=\sqrt{{{(x}_{}^{A}-{x}_{}^{B})}^{2}+{{(y}_{}^{A}-{y}_{}^{B})}^{2}+{{(z}_{}^{A}-{z}_{}^{B})}^{2}} $$

    式中,D为断裂发生的位错量,( $ {x}_{0}^{A},{y}_{0}^{A},{z}_{0}^{A} $ )、( $ {x}_{0}^{B},{y}_{0}^{B},{z}_{0}^{B} $ )分别为观测点AB的初始空间坐标,( $ {x}_{}^{A},{y}_{}^{A},{z}_{}^{A} $ )、( $ {x}_{}^{B},{y}_{}^{B},{z}_{}^{B} $ )分别为观测点AB同一时刻GNSS的定位坐标, $ {L}_{0}、L $ 分别为断裂发生位错前、后2个观测点间的基线长度。

    当不考虑断裂的竖向位错时,有:

    $$ {L}_{0}=\sqrt{{{(x}_{0}^{A}-{x}_{0}^{B})}^{2}+{{(y}_{0}^{A}-{y}_{0}^{B})}^{2}} $$ (4)
    $$ L=\sqrt{{{(x}_{}^{A}-{x}_{}^{B})}^{2}+{{(y}_{}^{A}-{y}_{}^{B})}^{2}} $$ (5)

    将位错量分解到3个空间坐标方向上,可确定3个方向的位错量,有:

    $$ \left\{\begin{aligned}&X={(x}^{A}-{x}_{0}^{A})-{(x}^{B}-{x}_{0}^{B})\\ &Y={(y}^{A}-{y}_{0}^{A})-{(y}^{B}-{y}_{0}^{B})\\ &Z={(z}^{A}-{z}_{0}^{A})-{(z}^{B}-{z}_{0}^{B})\end{aligned}\right. $$ (6)
    1.3.1   观测点选址

    在云南境内管道将穿越多条断裂,根据地震地质资料,选择龙陵-瑞丽断裂、凤仪-定西岭断裂及小江断裂西支作为断裂位错观测对象。历史上这3条断裂均发生过7级以上的强震,仍存在再次发生7级以上破坏性地震的可能。

    为保证观测结果能够真实反映断裂两盘的位错情况,观测站址应选择在断层两盘坚固完整的基岩上,不受地表水土流失、滑坡等因素影响。站址场地应开阔、无明显遮挡物。经实地考察,根据管道穿越断裂处地形地貌、断裂两盘情况,3个断裂位错GNSS观测子系统观测点布设如图1图3所示。

    图 1  龙陵-瑞丽跨断层观测系统布设示意图
    Figure 1.  Layout of Longling-Ruili cross fault monitoring system
    图 2  凤仪-定西岭跨断层观测系统布设示意图
    Figure 2.  Layout of Fengyi-Dingxiling cross fault monitoring system
    图 3  小江断裂西支跨断层观测系统布设示意图
    Figure 3.  Layout of fault monitoring system across the west branch of Xiaojiang fault

    (1)镇安观测站

    镇安GNSS观测站距龙陵县城17 km,距龙陵县镇安镇5 km,位于龙陵-瑞丽断裂上。龙陵-瑞丽断裂经龙陵、芒市、遮放、瑞丽盆地,延入缅甸中央盆地,全长约135 km,其北端与高黎贡断裂呈弧形相接。1976年5月29日,该断裂曾发生7.3级和7.4级大地震。经实地考察,选择AB点作为断裂两盘GNSS观测点场址,两点间距159 m,如图1所示。

    (2)定西岭观测站

    弥渡定西岭GNSS观测站位于弥渡县红岩镇定西岭,距弥渡红岩镇14 km,处在凤仪-定西岭断裂上。凤仪-定西岭断裂属红河大断裂带,红河断裂带是青藏高原东缘1条重要的大地构造和活动块体边界,该断裂及其周边地区地震活动较强烈,1925年曾发生凤仪7级地震。实地考察显示,在管道与断裂交汇处,断裂分为2支(图2),为观测交汇处断裂的位错,在被断裂分割的3个地块上各选择1个GNSS观测点,构成如图2所示的三角形测量体系,其两两间距分别为300、321、258 m。

    (3)小碑当观测站

    小碑当GNSS观测站位于寻甸县羊街镇小碑当村,距寻甸县城17 km,距羊街镇6 km,处在小江断裂带西支上。小江断裂带地处青藏高原东南侧,近SN走向,全长约400 km,系川滇菱形地块与华南地块的边界,其北起云南巧家,向南经东川分为东、西支,近平行延伸至抚仙湖,然后呈扫帚状向南撇开为多条断层。其构造运动十分活跃,地震活动频繁并诱发多种类型的地质灾害,1833年小江断裂带西支曾发生嵩明8级大地震。小江断裂与前述2条断裂相比,规模更大,在管道与断裂交汇处,断裂破碎带宽近300 m,两盘稳定地块间距应≥300 m,按相关规范的规定,选择的观测点如图3所示,观测点间距400 m。

    1.3.2   GNSS断裂位错观测子系统建设

    GNSS断裂位错观测子系统由2~3个跨断裂GNSS固定观测点观测墩、GNSS测量系统、供电系统、通信系统组成。

    (1)观测墩

    按照《全球导航卫星系统连续运行基准站网技术规范》(GB/T 28588—2012)(中华人民共和国国家质量监督检验检疫总局等,2012),GNSS混凝土观测墩应建在基岩上,如场地基岩埋深>20 m,则混凝土观测墩建在地下20 m深度处。观测墩高出地面3~5 m。建成的观测墩和观测房如图4图5所示,观测房内包含1个观测墩。

    图 4  跨断层变形观测系统GNSS固定观测墩
    Figure 4.  GNSS fixed observation pier of cross fault deformation monitoring system
    图 5  跨断层变形观测系统GNSS观测房(含观测墩)
    Figure 5.  Cross fault deformation monitoring system GNSS observation room (Including observation pier)

    (2)供电和通信

    供电系统采用太阳能直流供电,太阳能供电系统以最大连续阴雨10 d计算,可在连续10 d无日照天气时为设备供电。

    获取的观测数据通过4G移动网络远程传输到昆明数据中心,并可实现对设备的远程操控。

    (3)GNSS测量系统

    GNSS测量系统由GNSS天线和信号接收机组成,其技术性能指标如表1所示,采用S10接收机。

    表 1  GNSS测量系统技术性能指标
    Table 1.  Technical performance indexes of GNSS measurement system
    指标 基本性能
    通道 ≥24(并行)
    测量方式 独立完整的码与载波相位测量、高精度的多重相关 L1/L2 伪距测量
    采样 最大20 Hz
    RTK网络 VRS、FKP、MAC
    数据存储 ≥32 MB
    通信 RJ45、RS232,无线802.11 b
    功耗 2 W
    电源 10.5~28 V 直流,带压电保护功能。
    工作温度 −30 ℃~70 ℃
    静态精度(MonNET 后处理) 水平:3 mm+0.5 ppm RMS
    垂直:5 mm+1 ppm RMS
    快速静态基线精度(后处理)
    水平:5 mm+0.5 ppm
    垂直:10 mm+1 ppm
    动态定位:RTK
    水平:10 mm+1 ppm
    垂直:20 mm+1 ppm
    天线类型 L1/L2 零相位微对中天线,抑径板可以减少多路径干扰
    天线增益 50 dB,内置低噪音放大器
    使用温度 −40 ℃~70 ℃
    天线包装 防水、密封,高技术材料外壳机械强度高
    下载: 导出CSV 
    | 显示表格

    率先建成的镇安、定西岭观测站于2019年9月5日投入运行。数据连续采集,采样间隔为15 s,数据每2 h生成1个观测文件,至2019年10月7日,共获取5 400条观测数据,期间由于故障等原因,2个观测站数据均存在部分缺失。数据文件格式为RINEX2.11,系统常驻可观测GPS系统卫星可达8颗,BDS卫星可观测到15颗,受数据处理方法的限制,本文处理的数据均源自GPS系统。

    以仪器安装完成后测定的观测点初始位置( $ {X}_{0} $ $ {Y}_{0} $ $ {Z}_{0} $ )为参考点,龙陵观测站和定西岭观测站观测点PPP定位结果相对于测点初始位置偏移随时间的变化分别如图6图7所示。由图6图7可知,在3个方向上,各时刻的观测值围绕着初始位置波动,反映出各时刻确定的坐标位置有差异,在基础未发生位移的情况下,可将其视为GNSS观测的干扰背景,这种干扰包括观测数据解算算法本身包含的误差、大气环境等产生的影响。

    图 6  龙陵观测站观测点PPP定位结果相对于测点初始位置偏移随时间的变化
    Figure 6.  Variation of PPP positioning results of Longling observation station relative to the initial position offset of measurement points with time
    图 7  定西岭观测站观测点PPP定位结果相对于测点初始位置偏移随时间的变化
    Figure 7.  Variation of PPP positioning result of Dingxiling observation station relative to the initial position offset of the measurement point with time

    确定GNSS观测站2个观测点同一时刻空间坐标后,即可按式(3)确定两测点间基线长度相对于初始长度的变化,当变化量超过测量误差时,则可认为两测点间的断裂两盘发生了错位。

    不考虑竖向的变化,水平面内两观测点间基线长度变化可表示为:

    $$ D=L-{L}_{0}=\sqrt{{({x}^{A}-{x}^{B})}^{2}+{({y}^{A}-{y}^{B})}^{2}}-\sqrt{{{(x}_{0}^{A}-{x}_{0}^{B})}^{2}+{{(y}_{0}^{A}-{y}_{0}^{B})}^{2}} $$ (7)

    对于竖向点间距的变化,有:

    $$ Z={(z}^{A}-{z}_{0}^{A})-{(z}^{B}-{z}_{0}^{B}) $$ (8)

    2个观测站观测点间水平面内基线长度在不同时刻的变化如图8图9所示,2个观测站观测点间竖向间距在不同时刻的变化如图10图11所示。

    图 8  龙陵观测站观测点间水平面内基线长度的变化
    Figure 8.  Change of baseline length in horizontal between observation points of Longling observation station
    图 9  定西岭观测站观测点间水平面内基线长度的变化
    Figure 9.  Change of baseline length in horizontal between observation points of Dingxiling observation station
    图 10  龙陵观测站观测点间竖向间距的变化
    Figure 10.  Change of vertical distance between two measuring points at Longling observation station
    图 11  定西岭观测站观测点间竖向间距的变化
    Figure 11.  Change of vertical point spacing between observation points of Dingxiling observation station

    图8图11可知,2个观测站观测点间基线长度随时间变化,在所处理的时间段内,这种变化的最大值在水平向与竖向均≤10 cm,且无突出干扰值。由于在该时段内,未发生大的地震事件,可认为2个观测站间距在该波动范围内为正常状态。假定不发生地震时,2个观测站基线长度变化结果稳定在该数值范围内,则需解决的问题是观测时产生的测量误差满足管道变形安全要求。

    GNSS观测数据解算误差有多种来源,如卫星轨道摄动、电离层、对流层、大气含水量、多路径效应影响等(李征航等,2010),目前虽发展了消除这些影响的方法,但在观测条件不利的情况下,如降雨、降水等,误差会相应增大。

    本文利用G-Nut/Anubis(以下简称Anubis)对布设的5个观测站数据从多路径效应、单点定位精度、信噪比、数据利用率及周跳方面进行数据质量分析(张涛等,2017陈秀德等,2018康朝虎等,2018)。

    (1)多路径效应

    多路径效应是指接收机收到卫星发射的信号以外,同时接收到观测站附近物体反射的信号现象,接收到的不同路径信号会与原本接收到的卫星信号叠加,从而产生时延效应,带来误差,且该效应受环境的影响较大。不同观测站观测数据的多路径误差如图12所示。多路径均方根值受观测站周围环境影响,值越小表示数据质量越高。GPS系统观测的多路径误差多在40 cm上下浮动,GLONASS系统多路径误差较大。

    图 12  定西岭、龙陵观测站观测数据的多路径误差
    Figure 12.  Multipath error of observation data at Dingxiling and Longling observation stations

    (2)单点定位精度

    不同观测站观测数据的单点定位误差如图13所示。单点定位误差值越小,精度越高。系统观测值NE向的误差均<1 m,U向误差集中于2 m左右。

    图 13  定西岭、龙陵观测站单点定位误差
    Figure 13.  Single point positioning error of Dingxiling and Longling observation stations

    (3)信噪比

    信噪比为载波信号与噪声的比值,影响信噪比的因素较多,接收机本身、观测环境及多路径效应均会对其造成一定影响,该值不仅可用于评估观测数据质量,同时也可作为GNSS接收机性能优劣的参考(布金伟等,2017姚文敏等,2019)。不同观测站观测数据的信噪比如图14所示。由图14可知,观测数据信噪比基本在40 dB左右,各信道信噪比相差较小,较稳定。

    图 14  定西岭、龙陵观测站观测数据信噪比
    Figure 14.  Signal to noise ratio of observation data of Dingxiling and Longling observation stations

    (4)数据完整性

    观测数据完整性是观测质量的直接决定因素,由于各种偶然因素及传输过程中的损耗,导致接收机实际接收的观测数据少于理论数据量。周跳也称为整周跳变,发生的次数越多,表示丢失的观测数据越多,可反映卫星信号失锁导致的整周计数偏差情况。对图14观测数据进行统计,可知GPS系统的数据可利用率远高于GLONASS系统,而其周跳次数经长时统计,所有台站周跳频次维持在较低水平,能够保证观测精度及数据可靠性(表2)。

    表 2  数据利用率及周跳频次统计
    Table 2.  Statistics of data utilization and weekly hop frequency
    测站 数据利用率/% 周跳频次
    MDGCDS 100 0.006 2
    MDGCF 100 0.012 0
    MDGCDX 83 0.002 8
    LLGCD 100 0.004 4
    LLGCF 100 0.005 6
    下载: 导出CSV 
    | 显示表格

    观测的最终目的是对天然气管道进行保护,观测结果是站点的坐标变动,如果观测精度低于管道容许压缩应变范围,观测失去意义(张慧敏,2009孟立朋等,2016)。

    中缅油气管道跨越断层时一般在管沟中平直铺设,荷载垂直于管道轴向且沿轴向均匀分布,管道尺寸参数如表3所示。

    表 3  管道尺寸参数
    Table 3.  Design parameters of different pipeline projects
    管道项目 外径/m 壁厚/m 设计压力/MPa
    天然气管道 0.010 6 0.022 9 15
    原油管道 0.813 0 0.028 6 15
    成品油管道 0.219 1 0.009 5 15
    下载: 导出CSV 
    | 显示表格

    根据理论,管道受到地震影响时,由于管道自身的抗震设计,不会直接发生破裂,但如果管道受地震影响过大,其变形量超出容许范围时,会发生不可逆的损害,管道容许变形长度对于GNSS观测系统警报阈值计算至关重要。确定管道容许变形长度时分为2种情况,一种为管道受拉状态下产生的最大位移,另一种为管道受压状态下的最大位移。当管道受到拉伸、压缩时管道最大容许长度变化及土壤与管道外表面之间单位长度上的摩擦力根据《输油(气)钢质管道抗震设计规范》(SY/T 0450—2004)(国家发展和改革委员会,2004)中的相关规定计算。

    考虑管道变形时,还需考虑弯曲变形对管道的危害,通过计算轴向运行弯曲变形能力,可得到管道能够承受的沉降偏差最大值,根据《输气管道工程设计规范》(GB 50251—2015)(中华人民共和国住房和城乡建设部,2015)中的相关规定计算。

    不同管道项目容许拉伸量及容许压缩量如表4所示。

    表 4  管道容许变形计算
    Table 4.  Calculation of pipeline ultimate deformation
    管道项目 容许拉伸量/cm 容许压缩量/cm 容许弯曲量/cm
    天然气管道 41.13 27.57 8.4
    原油管道 132.59 70.56 152.9
    成品油管道 66.08 20.04 63.77
    下载: 导出CSV 
    | 显示表格

    由于PPP精度随观测时长而收敛,进行5、6 h以上时长收敛,可达厘米级精度(李浩军等,2010赵爽,2018淦邦,2020)。

    具体定位精度的评价因素较多,主要包括平面位置精度因子(HDOP)、高程精度因子(VDOP)、空间位置精度因子(PDOP)、接收机钟差精度因子(TDOP)及几何精度因子(GDOP),以上因子带来的影响体现为观测值与真值的差值xyz进行坐标变换后得到的对应值,即站心坐标系下对应于观测站本身水平和高程方向的定位误差(代桃高,2017王何鹏等,2019),转换公式为:

    $$ \left[\begin{array}{c}E\\ N\\ U\end{array}\right]=\left[\begin{array}{*{20}{c}}-\sin\lambda & \cos\lambda &0\\ -\cos\lambda \sin\varphi & -\cos\lambda \sin\varphi & \cos\varphi \\ \cos\lambda \cos\varphi & \sin\lambda \cos\varphi & \sin\varphi \end{array}\right]\left[\begin{array}{c}\Delta_x\\ \Delta_ y\\ \Delta_z\end{array}\right] $$ (9)

    式中, $ \lambda $ $\varphi$ 分别为该点所在的经、纬度。

    大理弥渡和保山龙陵观测站GNSS观测精度如图15图16所示。由图15可知,采样时间前段数据值波动较大,这是因为台站安置初期会出现基座沉降导致的安置点位变化及PPP初期带来的较大观测误差,PPP观测须经过长时间收敛才会趋于稳定值,如图15图16所示,布设观测站观测结果的均方根误差基本在第400个历元后开始收敛并趋于稳定,取第400个历元后的误差均值为误差收敛值。

    图 15  定西岭观测站基线误差解算结果
    Figure 15.  Baseline error calculation results of Dingxiling observation station
    图 16  龙陵观测站基线误差解算结果
    Figure 16.  Baseline error calculation results of Longling observation station

    处理后PPP定位误差的收敛值如表5所示。由表5可知,NE向误差收敛值均处于毫米级水平,U向误差收敛值由于GNSS本身的局限性而处于厘米级水平,将基线误差收敛值与管道容许拉伸量和容许压缩量进行对比,可知基线误差收敛值均小于3种管道容许拉伸量和容许压缩量,说明GNSS观测精度能够满足油气管道检测需求。

    表 5  PPP基线误差收敛值
    Table 5.  Convergence value of PPP baseline error
    方向 定西岭观测站基线
    误差收敛值/cm
    龙陵观测站基线
    误差收敛值/cm
    N 0.514 0.412
    E 0.381 0.379
    U 3.341 4.314
    下载: 导出CSV 
    | 显示表格

    在系统实际应用中,由于偶然因素多,加上管道自身存在微小变形,会导致起始坐标、起始间距发生变化,综合各种因素,对于GNSS管道观测的阈值,不能单纯用极限变形量与极限误差收敛值评估,应结合数据的摇摆及观测误差进行考虑,可按下式计算:

    $$ {\sigma }_{{\rm{t}}}=\left[{\Delta L}_{{\rm{t}}}\right]+{\sigma }_{g}+{\sigma }_{i} $$ (10)
    $$ {\sigma }_{{\rm{c}}}=\left[{\Delta L}_{{\rm{c}}}\right]+{\sigma }_{g}+{\sigma }_{i} $$ (11)
    $$ {\sigma }_{{\rm{h}}}={h}_{\max}+{\sigma }_{g}+{\sigma }_{i} $$ (12)

    式中, $ {\sigma }_{{\rm{t}}} $ 为GNSS管道拉伸观测阈值, ${\sigma }_{{\rm{c}}}$ 为GNSS管道压缩观测阈值, $ {{\sigma }}_{{\rm{h}}} $ 为GNSS管道偏差阙值, $ {\sigma }_{g} $ 为第g条基线的误差收敛值, $ {\sigma }_{i} $ 为第i个观测点一时间序列上实测水平距或高程差与真实水平距和高程差间差值的标准差。

    不同观测站观测值与真实值差值的标准差如表6所示。观测时长内每条基线间的警报阈值计算结果如表7表9所示。

    表 6  不同观测站观测值与真实值差值的标准差
    Table 6.  Difference standard deviation between observed values and true values at different stations
    基线段 标准差/cm
    定西岭观测站A-C 0.058
    定西岭观测站B-C 0.060
    定西岭观测站A-B 0.058
    龙陵观测站A-B 0.029
    下载: 导出CSV 
    | 显示表格
    表 7  N向GNSS管道观测阈值
    Table 7.  Observation threshold of GNSS pipeline in horizontal N direction
    基线段 拉伸阈值/cm 压缩阈值/cm 弯曲阈值/cm
    天然气管道 原油管道 成品油管道 天然气管道 原油管道 成品油管道 天然气管道 原油管道 成品油管道
    定西岭观测站A-C 41.702 133.162 66.652 28.142 71.132 20.612 8.972 153.472 64.342
    定西岭观测站B-C 41.704 133.164 66.654 28.144 71.134 20.614 8.974 153.474 64.344
    定西岭观测站A-B 41.702 133.162 66.652 28.142 71.132 20.612 8.972 153.472 64.342
    龙陵观测站A-B 41.571 133.031 66.521 28.011 71.001 20.481 8.943 153.443 64.313
    下载: 导出CSV 
    | 显示表格
    表 8  E向GNSS管道观测阈值
    Table 8.  Observation threshold of GNSS pipeline in horizontal Edirection
    基线段 拉伸阈值/cm 压缩阈值/cm 弯曲阈值/cm
    天然气管道 原油管道 成品油管道 天然气管道 原油管道 成品油管道 天然气管道 原油管道 成品油管道
    定西岭观测站A-C 41.569 133.029 66.519 28.007 70.999 20.479 8.839 153.339 64.209
    定西岭观测站B-C 41.571 133.031 66.521 28.011 71.001 20.481 8.841 153.341 64.211
    定西岭观测站A-B 41.569 133.029 66.519 28.009 70.999 20.479 8.839 153.339 64.209
    龙陵观测站A-B 41.538 132.998 66.488 27.978 70.968 20.448 8.808 153.308 64.178
    下载: 导出CSV 
    | 显示表格
    表 9  U向GNSS管道观测阈值
    Table 9.  Observation threshold of GNSS pipeline in vertical U direction
    基线段 拉伸阈值/cm 压缩阈值/cm 弯曲阈值/cm
    天然气管道 原油管道 成品油管道 天然气管道 原油管道 成品油管道 天然气管道 原油管道 成品油管道
    定西岭观测站A-C 44.529 135.989 69.479 30.969 73.959 23.439 11.799 156.299 67.169
    定西岭观测站B-C 44.531 135.991 69.481 30.971 73.961 23.441 11.801 156.301 67.171
    定西岭观测站A-B 44.529 135.989 69.479 30.969 73.959 23.439 11.799 156.299 67.169
    龙陵观测站A-B 45.473 136.933 70.423 31.913 73.93 24.383 12.743 157.243 68.113
    下载: 导出CSV 
    | 显示表格

    本文针对中缅油气管道工程范围内的2条断裂带,利用2019年沿断裂布设的GNSS固定台站观测数据,确保观测数据在误差允许范围内对断裂带变形趋势进行长期观测与分析,以降低中缅油气管道受地震的影响。当精度达到容许拉伸应变误差允许范围内时,布设GNSS观测站对断裂带进行有针对性的观测切实可行,相较于从IGS特征站观测数据内插得到断裂带数据更直接快捷,能够直观反映断裂带变形趋势,且规避了参考站受地震波影响导致观测结果产生的误差,从而保证准确获得地震位移波形。

    本研究还存在部分问题,PPP观测精度虽随观测时长而收敛,且观测时长满足其收敛要求,但由于观测精度与拉伸应变误差较接近,不能完全保证突发的偶然误差能够得到处理。有研究表明采用BDS与GPS双系统乃至引入3个以上系统能够显着提高PPP观测精度,而本项目使用了GPS单系统,仍有改进空间;进行精度论证时,仅评价了NEU向上的精度差,由于NEU向均为一维方向,并未进行三维精度论证,结果不够严谨,需改进。

  • 图  1  有限元模型

    Figure  1.  3D finite element model

    图  2  修正的Davidenkov模型剪应力-剪应变关系曲线

    Figure  2.  Shear stress-shear strain curve for the modified Davidenkov model

    图  3  输入地震动时程和傅里叶谱

    Figure  3.  Input time history of ground motion and Fourier spectrum

    图  4  管廊监测点布置

    Figure  4.  Utility tunnel monitoring points

    图  5  管廊横断面应力分布(Mises应力)

    Figure  5.  Utility tunnel cross-sectional stress distribution (Mises stress)

    图  6  管廊横断面内力分布

    Figure  6.  Internal force distribution in the cross section of the utility tunnel

    图  7  管片纵向环间张开量

    Figure  7.  The amount of opening between the longitudinal rings of the utility tunnel

    图  8  管廊环间张开量时程曲线

    Figure  8.  Time history curve of opening amount between utility tunnel rings

    图  9  管廊结构损伤因子云图

    Figure  9.  Damage factors of utility tunnel structure

    表  1  土体参数

    Table  1.   Soil parameters

    土层名称最大动剪切
    模量Gmax/kPa
    泊松比参数A参数B参考剪应变γr参数β参考围压/kPa参数C1参数C2门槛剪应变γth参数m参数n
    1粉质黏土混粉土层58 6440.410.45950.000 3770.0031 0001.0510.1430.000 20.345668.9
    2粉土层123 0320.410.45950.000 5080.0031 0000.9210.1630.000 20.345668.9
    1粉细砂层176 4000.410.45950.000 4780.0031 0000.8130.1850.000 20.345668.9
    2细砂层202 7520.410.45950.000 5410.0031 0000.8130.1850.000 20.345668.9
    1中粗砂层219 9780.410.45950.000 6170.0031 0000.8130.1850.000 20.345668.9
    下载: 导出CSV

    表  2  混凝土参数

    Table  2.   Concrete parameters

    混凝土强度
    等级
    密度/
    (kg·m−3
    泊松比弹性模量/GPa膨胀角/(°)偏心率双轴与单轴抗压
    强度之比fb0/fc0
    屈服形态影响
    参数K
    黏性系数
    C302 3000.230300.11.160.6670.000 5
    C602 5000.236380.11.160.6670.0005
    下载: 导出CSV

    表  3  钢筋与螺栓参数

    Table  3.   Rebar and bolt parameters

    材料密度/(kg·m−3弹性模量/GPa泊松比屈服强度/MPa
    HRB400钢筋7 8002000.3400
    10.9级螺栓7 8002000.3900
    下载: 导出CSV
  • 韩俊艳, 王小强, 郭之科等, 2022. 纵向非一致激励下埋地管道的振动台试验研究. 岩石力学与工程学报, 41(6): 1256—1266 doi: 10.13722/j.cnki.jrme.2021.0515

    Han J. Y. , Wang X. Q. , Guo Z. K. , et al. , 2022. Shaking table test of buried pipelines under longitudinal non-uniform excitation. Chinese Journal of Rock Mechanics and Engineering, 41(6): 1256—1266. (in Chinese) doi: 10.13722/j.cnki.jrme.2021.0515
    黄雨, 于淼, Bhattacharya S. , 2013.2011年日本东北地区太平洋近海地震地基液化灾害综述. 岩土工程学报, 35(5): 834—840

    Huang Y. , Yu M. , Bhattacharya S. , 2013. Review on liquefaction-induced damages of soils and foundations during 2011 of the Pacific Coast of Tohoku Earthquake (Japan). Chinese Journal of Geotechnical Engineering, 35(5): 834—840. (in Chinese)
    蒋录珍, 李杰, 陈隽, 2015. 非一致地震激励下综合管廊接头响应数值模拟. 世界地震工程, 31(2): 101—107

    Jiang L. Z. , Li J. , Chen J. , 2015. Numerical simulation of utility tunnel joint effects under non-uniform earthquake excitation. World Earthquake Engineering, 31(2): 101—107. (in Chinese)
    钱七虎, 陈晓强, 2007. 国内外地下综合管线廊道发展的现状、问题及对策. 地下空间与工程学报, 3(2): 191—194 doi: 10.3969/j.issn.1673-0836.2007.02.001

    Qian Q. H. , Chen X. Q. , 2007. Situation, problems and countermeasures of utility tunnel’ development in China and abroad. Chinese Journal of Underground Space and Engineering, 3(2): 191—194. (in Chinese) doi: 10.3969/j.issn.1673-0836.2007.02.001
    史晓军, 陈隽, 李杰, 2008, 地下综合管廊大型振动台模型试验研究. 地震工程与工程振动, 28(6): 116—123

    Shi X. J., Chen J., Li J., 2008. Shaking table test on underground utility tunnel. Earthquake Engineering and Engineering Dynamics, 28(6): 116—123. (in Chinese)
    施有志, 华建兵, 阮建凑等, 2018. 地下综合管廊地震动力响应三维数值分析. 工程地质学报, 26(3): 785—793 doi: 10.13544/j.cnki.jeg.2017-163

    Shi Y. Z. , Hua J. B. , Ruan J. C. , et al. , 2018. Three dimensional numerical analysis of dynamic response of underground utility tunnels during earthquake. Journal of Engineering Geology, 26(3): 785—793. (in Chinese) doi: 10.13544/j.cnki.jeg.2017-163
    汤鹏, 杨明, 庄海洋等, 2022. 越江电力地下综合管廊结构横向抗震性能研究. 防灾减灾工程学报, 42(3): 507—515 doi: 10.13409/j.cnki.jdpme.202006025

    Tang P. , Yang M. , Zhuang H. Y. , et al. , 2022. Lateral seismic performance of the utility tunnel crossing the Yangzi River. Journal of Disaster Prevention and Mitigation Engineering, 42(3): 507—515. (in Chinese) doi: 10.13409/j.cnki.jdpme.202006025
    袁勇, 朱力, 禹海涛等, 2022. 综合管廊抗震分析研究进展综述. 隧道建设(中英文), 42(11): 1821—1831

    Yuan Y. , Zhu L. , Yu H. T. , et al. , 2022. Review on seismic analysis of utility tunnels. Tunnel Construction, 42(11): 1821—1831. (in Chinese)
    仉文岗, 韩亮, 陈志雄等, 2020. 双仓综合管廊抗震性能模型试验研究. 岩土工程学报, 42(1): 100—108

    Zhang W. G. , Han L. , Chen Z. X. , et al. , 2020. Model tests on seismic performance of double-box underground utility tunnel. Chinese Journal of Geotechnical Engineering, 42(1): 100—108. (in Chinese)
    中华人民共和国住房和城乡建设部, 2011. GB 50010—2010 混凝土结构设计规范. 北京: 中国建筑工业出版社.

    Ministry of Housing and Urban-Rural Development of the People’s Republic of China, 2011. GB 50010—2010 Code for design of concrete structures. Beijing: China Architecture & Building Press. (in Chinese)
    中华人民共和国住房和城乡建设部, 国家市场监督管理局, 2021. GB/T 51438—2021 盾构隧道工程设计标准. 北京: 中国建筑工业出版社.

    Ministry of Housing and Urban-Rural Development of the People’s Republic of China, State Administration for Market Regulation, 2021. GB/T 51438—2021 Standard for design of shield tunnel engineering. Beijing: China Architecture & Building Press. (in Chinese)
    Chen G. X. , Zhao D. F. , Chen W. Y. , et al. , 2019. Excess pore-water pressure generation in cyclic undrained testing. Journal of Geotechnical and Geoenvironmental Engineering, 145(7): 04019022. doi: 10.1061/(ASCE)GT.1943-5606.0002057
    Chen G. X. , Wang Y. Z. , Zhao D. F. , et al. , 2021. A new effective stress method for nonlinear site response analyses. Earthquake Engineering & Structural Dynamics, 50(6): 1595—1611.
    Duan X. , Dong Q. , Ye W. J. , 2019. Experimental study on seismic performance of prefabricated utility tunnel. Advances in Civil Engineering, 2019: 8968260.
    Ruan B. , Zhao K. , Wang S. Y. , et al. , 2019. Numerical modeling of seismic site effects in a shallow estuarine bay (Suai Bay, Shantou, China). Engineering Geology, 260: 105233. doi: 10.1016/j.enggeo.2019.105233
    Wang H. F. , Lou M. L. , Zhang R. L. , 2017. Influence of presence of adjacent surface structure on seismic response of underground structure. Soil Dynamics and Earthquake Engineering, 100: 131—143. doi: 10.1016/j.soildyn.2017.05.031
    Zhang J. H. , Yuan Y. , Yu H. T. , 2019. Shaking table tests on discrepant responses of shaft-tunnel junction in soft soil under transverse excitations. Soil Dynamics and Earthquake Engineering, 120: 345—359. doi: 10.1016/j.soildyn.2019.02.013
  • 期刊类型引用(3)

    1. 熊政辉,曲哲,宁超列,卜春尧,尹建明,薄涛. 我国地震巨灾保险试点实践对比研究. 灾害学. 2025(01): 48-52+159 . 百度学术
    2. 李云仙,子芝兰. 基于AHP-FCE的地震巨灾保险运行效果评价——以大理州政策性农房地震巨灾保险为例. 保险理论与实践. 2024(10): 33-55 . 百度学术
    3. 景冰冰,施唯. 2023年土耳其M_S7.8地震灾区损失与保险研究. 震灾防御技术. 2023(03): 495-504 . 本站查看

    其他类型引用(6)

  • 加载中
图(9) / 表(3)
计量
  • 文章访问数:  123
  • HTML全文浏览量:  58
  • PDF下载量:  20
  • 被引次数: 9
出版历程
  • 收稿日期:  2022-10-19
  • 刊出日期:  2023-03-31

目录

/

返回文章
返回