2019年2月2日新疆塔城<i>M</i><sub>S</sub>5.2地震房屋震害特点与发震构造

黄帅堂; 马建; 常想德; 阿里木江·亚力昆; 任静

doi:10.11899/zzfy20200211

2019年2月2日新疆塔城M_S5.2地震房屋震害特点与发震构造

doi: 10.11899/zzfy20200211

1.
新疆维吾尔自治区地震局, 乌鲁木齐 830011
2.
中国地震台网中心, 北京 100045

基金项目:

中国地震局地震预测研究所基本科研业务费专项重点项目(援疆项目) 2016IES0103

地震应急青年课题 CEAEDEM201916

新疆地震局局基金课题 201910

详细信息

作者简介:
黄帅堂, 男, 生于1990年。硕士研究生。主要从事活动构造、构造地貌方面研究。E-mail:shuaith@126.com

计量
- 文章访问数: 75
- HTML全文浏览量: 68
- PDF下载量: 11
- 被引次数: 0
出版历程
- 收稿日期: 2019-02-19
- 刊出日期: 2020-06-20

Seismic Disaster Characteristics of Buildings and Seismogenic after the Xinjiang M_S5.2 Earthquake on February 2, 2019

1.
Earthquake Agency of Xinjiang Uygur Autonomous Region, Urumqi 830011, China
2.
China Earquake Networks Center, Beijing 100045, China

摘要

摘要: 2019年2月2日新疆塔城市发生M_S5.2地震，根据现场震害调查确定本次地震极震区烈度为Ⅵ度。地震造成震中附近20世纪八、九十年代建造的个别土木结构房屋出现中等程度以上破坏，少数轻微破坏(含旧裂缝加宽)，多数基本完好；其他结构类型房屋多数基本完好。本次地震震级小、有感范围广，震区土木结构房屋缺少抗震措施、年久失修、地基处理不当是造成房屋出现破坏的主要原因。综合分析等震线分布及深部构造特征，初步判定本次地震发震构造为塔克台断裂向塔城盆地迁移、扩展的隐伏断裂。
- 塔城M_S5.2地震 /
- 房屋破坏 /
- 震害特点 /
- 发震构造
Abstract: An M_S5.2 earthquake occurred on February 2, 2019 in Tacheng, Xinjiang. The earthquake intensity is VI degree. Based on a field survey and a related data analysis, we reached two conclusions. First, the earthquake didn't cause large-scale collapse and destruction. The earthquake caused a moderate degree of damage to individual soil-wood in the 1980s and 1990s near the epicenter, and a few minor damages (including the increase of old cracks); most of the other types of buildings are basically intact. Although the earthquake has a small magnitude it has a wide sensible area. Lack of seismic construction measures, disrepair for years and improper foundation treatment is the main cause of damage for earthquake-stricken houses. Second, the seismogenic structure causing this earthquake was located at Tacheng basin, a hidden fault that migrated and expanded from the Taktai fault to the Tacheng Basin.
- Tacheng M_S5.2 earthquake /
- Housing destruction /
- Characteristics of seismic damage /
- Seismogenic structure

HTML全文

引言

中缅油气管道工程是我国重大能源战略项目，全长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断裂变形观测子系统能够有效观测破坏性地震引起的断层变形量。

1. GNSS断裂变形观测方法及观测系统设计

1.1 精密单点定位原理

精密单点定位算法采用双频相位和伪距观测数据，利用国际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.2 基于PPP技术的断裂位错观测方法

在断裂两盘各布设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} $ ）分别为观测点A、B的初始空间坐标，（ $ {x}_{}^{A},{y}_{}^{A},{z}_{}^{A} $ ）、（ $ {x}_{}^{B},{y}_{}^{B},{z}_{}^{B} $ ）分别为观测点A、B同一时刻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 GNSS断裂位错观测系统设计及建设

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级大地震。经实地考察，选择A、B点作为断裂两盘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

| 显示表格

2. 观测数据及其处理

2.1 观测数据

率先建成的镇安、定西岭观测站于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

下载: 全尺寸图片幻灯片

2.2 观测站测点间基线长度

确定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

下载: 全尺寸图片幻灯片

3. 结果分析

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

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

3.1 误差分析

本文利用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所示。单点定位误差值越小，精度越高。系统观测值N、E向的误差均＜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

| 显示表格

3.2 观测结果的有效性

观测的最终目的是对天然气管道进行保护，观测结果是站点的坐标变动，如果观测精度低于管道容许压缩应变范围，观测失去意义（张慧敏，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），以上因子带来的影响体现为观测值与真值的差值∆_x、∆_y、∆_z进行坐标变换后得到的对应值，即站心坐标系下对应于观测站本身水平和高程方向的定位误差（代桃高，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可知，N、E向误差收敛值均处于毫米级水平，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

| 显示表格

3.3 GNSS管道观测阈值

在系统实际应用中，由于偶然因素多，加上管道自身存在微小变形，会导致起始坐标、起始间距发生变化，综合各种因素，对于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

| 显示表格

4. 结语

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

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

图 1 塔城地震烈度、调查点及震区场地条件分布

Figure 1. The seismic intensity, survey point distribution and site conditions map of the Tacheng earthquake

下载: 全尺寸图片幻灯片

图 2 土木结构房屋受损情况

(a)禾斯阔甫尔村“干打垒”房屋中等程度破坏，后墙外闪，出现纵横墙间裂缝；(b)阿热散村老旧土木结构房屋后纵墙外闪；(c)、(d)阿热散村土木结构房屋已有裂缝加宽；(e)农二队土木结构房屋门洞角裂缝加宽、墙皮脱落；(f)科克塔尔村“干打垒”房屋出现纵横墙间裂缝

Figure 2. Damage of soil-wood by the earthquake

下载: 全尺寸图片幻灯片

图 3 其他结构房屋受损情况

(a)巴斯博孜达克东村老旧砖木结构房屋裂缝加宽；(b)农二队砖混结构房屋墙体裂缝加宽；(c)农二队砖混结构房屋预制板顶连接处出现裂缝；(d)滑雪场砖混结构房屋门洞角出现斜裂缝；(e)恰夏镇中心学校小学教学楼框架结构轻微破坏，填充墙墙体门洞角出现斜裂缝；(f)恰夏镇中心学校小学教学楼框架结构轻微破坏，廊洞部出现横向裂缝

Figure 3. Damage of other houses by the earquake

下载: 全尺寸图片幻灯片

图 4 震区安居富民房完好

Figure 4. Inact comfortable house by the earthquake

下载: 全尺寸图片幻灯片

图 5 区域地震构造

Figure 5. Regional seismic structure diagram

下载: 全尺寸图片幻灯片

图 6 塔克台断裂断错地貌

Figure 6. Offset geomorphology of Taktai fault

下载: 全尺寸图片幻灯片

图 7 塔城M_S5.2地震发震构造

Figure 7. Seismogenic model for the Tacheng M_S5.2 event

下载: 全尺寸图片幻灯片