Discussion of Relationship between the Wenchuan Earthquake and Lushan Earthquake from the Viewpoint of Coulomb Failure Stress Change and Spatial Distribution of Aftershocks
-
摘要: 本文以龙门山及周边地区为研究对象，考虑区域地质构造差异、主要活动断裂带、地表附加重力影响，建立能反映地表起伏和岩石圈分层结构的龙门山地区三维粘弹性有限元模型。以地壳水平运动速率观测值为约束条件重建研究区现今构造背景应力场，在此基础上分别模拟了汶川地震和芦山地震的发生机理。通过分析同震库仑破裂应力变化与余震空间分布的关系，探讨了2次地震主震对余震的触发作用以及汶川地震对芦山地震的影响。研究表明，汶川地震和芦山地震的余震大部分由其主震触发，汶川地震对芦山地震的余震有约6.78%的触发作用。汶川地震的同震库仑破裂应力在芦山地震主震位置的增加值约为0.016MPa，如果龙门山断裂带南段库仑破裂应力年累积速率按照0.4×10-3-0.6×10-3MPa·a-1计算，汶川地震使芦山地震提前了约27-40年。计算还表明汶川地震和芦山地震的发生使鲜水河断裂带南段和虎牙断裂的库仑破裂应力增加，这些断裂带在未来发生地震的可能性增加。Abstract: Taking the Longmenshan and its surrounding area as the research area, with consideration of the difference of regional geological structure, the main active fault zone, additional surface gravity, we constructed the 3-D viscoelastic finite element model of the Longmenshan area with irregular topography and the layered lithosphere structure. Using the observed values of crustal horizontal movement as constraint condition, the present tectonic background stress field in the research area is reconstructed. Based on it, the occurrence of Wenchuan earthquake and Lushan earthquake are simulated respectively. Through the analysis the relationship between coseismic Coulomb stress change and the spatial distribution of aftershocks, the trigger action of main shock to aftershock and the influence of Wenchuan earthquake on Lushan earthquake is investigated. The results show that the most aftershocks of Wenchuan earthquake and Lushan earthquake are triggered by their main shock, and only 6.78% aftershock of Lushan earthquake is triggered by Wenchuan earthquake. The coseismic Coulomb stress of the location of Lushan earthquake caused by Wenchuan earthquake is about 0.016MPa. Assuming tha the annual cumulative rate of Coulomb stress in the southern segment of the Longmenshan fault zone is 0.4-0.6×10-3MPa·a-1, the occurrence of Wenchuan earthquake causes the occurrence of Lushan earthquake nearly 27-40 years earlier than it should. Our results also show that the Wenchuan earthquake and Lushan earthquake cause the Coulomb stress increase in the southern segment of the Xianshuihe fault zone and Huya fault, which increases the possibility of earthquake in this area in the future.
-
引言
中缅油气管道工程是我国重大能源战略项目，全长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 ℃ 天线包装 防水、密封，高技术材料外壳机械强度高 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 time2.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 station3. 结果分析
由图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 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 根据理论，管道受到地震影响时，由于管道自身的抗震设计，不会直接发生破裂，但如果管道受地震影响过大，其变形量超出容许范围时，会发生不可逆的损害，管道容许变形长度对于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 由于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龙陵观测站基线
误差收敛值/cmN 0.514 0.412 E 0.381 0.379 U 3.341 4.314 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 表 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 表 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 表 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 4. 结语
本文针对中缅油气管道工程范围内的2条断裂带，利用2019年沿断裂布设的GNSS固定台站观测数据，确保观测数据在误差允许范围内对断裂带变形趋势进行长期观测与分析，以降低中缅油气管道受地震的影响。当精度达到容许拉伸应变误差允许范围内时，布设GNSS观测站对断裂带进行有针对性的观测切实可行，相较于从IGS特征站观测数据内插得到断裂带数据更直接快捷，能够直观反映断裂带变形趋势，且规避了参考站受地震波影响导致观测结果产生的误差，从而保证准确获得地震位移波形。
本研究还存在部分问题，PPP观测精度虽随观测时长而收敛，且观测时长满足其收敛要求，但由于观测精度与拉伸应变误差较接近，不能完全保证突发的偶然误差能够得到处理。有研究表明采用BDS与GPS双系统乃至引入3个以上系统能够显着提高PPP观测精度，而本项目使用了GPS单系统，仍有改进空间；进行精度论证时，仅评价了N、E、U向上的精度差，由于N、E、U向均为一维方向，并未进行三维精度论证，结果不够严谨，需改进。
-
图 1 研究区域 (根据邓起东等，2011)
Figure 1. Research area of this study (after Deng et al., 2011)
图 4 地壳水平运动速率模拟值与观测值的对比
Figure 4. Comparison of the simulated and the observed values of the crustal horizontal movement velocity
图 5 龙门山及周边地区主要断裂带上的库仑破裂应力年累积速率
Figure 5. Cumulative Coulomb failure stress of year in the main fault zone in Longmenshan and its surrounding area
图 6 汶川地震的同震库仑破裂应力变化与余震的空间分布
红色实心圆圈表示汶川地震主震，白色圆圈表示余震
Figure 6. Coseismic Coulomb stress change of the Wenchuan earthquake and the spatial distribution of aftershock
图 7 芦山地震的同震库仑破裂应力变化与余震的空间分布
红色实心五角星表示芦山地震主震，白色圆圈表示余震
Figure 7. Coseismic Coulomb stress change of the Lushan earthquake and the spatial distribution of aftershock (Solid red star represents the main shock of Lushan earthquake, the white circle represents the aftershock)
图 8 汶川地震的同震库仑破裂应力在芦山地震破裂面上的投影
红色实心圆圈表示汶川地震，红色实心五角星表示芦山地震
Figure 8. The projection of the coseismic Coulomb stress of the Wenchuan earthquake on the fracture plane of the Lushan earthquake
图 9 汶川地震和芦山地震共同产生的同震库仑破裂应力变化与余震的空间分布
红色实心五角星表示芦山地震，白色圆圈表示余震
Figure 9. Coseismic Coulomb stress change of the Wenchuan earthquake and Lushan earthquake as well as the spatial distribution of aftershock
图 10 汶川地震和芦山地震共同引起的龙门山地区主要活动断层的库仑破裂应力变化
Figure 10. Coulomb stress change of the main active faults in Longmenshan and its surrounding area caused by the Wenchuan earthquake and Lushan earthquake
表 1 研究区介质分层材料参数
Table 1. Material parameters of the layered medium in the research area
分层 深度/km 四川盆地 青藏高原东缘 E/1010·Pa υ η/1022·Pa·s E/1010·Pa υ η/1022·Pa·s 地表 0-0.5 3.75 0.21 0.8 3.75 0.21 0.8 上地壳 0.5-20 7.2 0.246 9.96 7.39 0.25 10.3 下地壳 20-Moho面 12.5 0.253 0.925 11.9 0.254 2.02 岩石圈 Moho面-100 17.5 0.265 0.05 17.50 0.265 0.05 软流圈 100-200 17.5 0.265 0.5 17.50 0.265 0.5 
下载: 导出CSV
表 2 汶川地震和芦山地震参数
Table 2. Earthquake parameters of the Wenchuan earthquake and Lushan earthquake
事件 经度 纬度 走向 倾角 滑动角 最大水平滑动量/m 破裂长度/km 数据来源 汶川 103.4°E 31.0°N 231° 35° 138° 4.9 > 300 中国地震台网中心；Global CMT；徐锡伟等，2008 芦山 102.89°E 30.31°N 212° 44° 92° 1.3 46.7 USGS；曾祥方等，2013；徐锡伟，2013；张勇等，2013
下载: 导出CSV
表 3 研究区内主要活动断层参数及汶川地震和芦山地震共同引起的库仑破裂应力变化
Table 3. Parameters of the main active faults in research area and the Coulomb stress change caused by the Wenchuan earthquake and Lushan earthquake
断裂名称 走向/° 倾角/° 滑动角/° 库仑破裂应力变化/MPa 龙门山断裂带东北段 225 60 180 -4.2-3.2 龙门山断裂带南段 220 60 90 -2-1.8 鲜水河断裂带南段 142-159.5 90 0-45 0.005-0.15 龙日坝断裂 205-229 60 135 -0.067-0 岷江断裂 180 45-60 45-70 -2.96-0 虎牙断裂 150 75 45 0~0.08
下载: 导出CSV
-
曹建玲, 石耀霖, 张怀, 王辉, 2009.青藏高原GPS位移绕喜马拉雅东构造结顺时针旋转成因的数值模拟.科学通报, 54(2):224-234. http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB200902017.htm 陈立春, 冉勇康, 王虎等, 2013.芦山地震与龙门山断裂带南段活动性.科学通报, 58(20):1925-1932. http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201320007.htm 陈连旺, 陆远忠, 刘杰等, 2001.1966年邢台地震引起的华北地区应力场动态演化过程的三维粘弹性模拟.地震学报, 23(5):480-491. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXB200105003.htm 陈连旺, 张培震, 陆远忠等, 2008.川滇地区强震序列库仑破裂应力加卸载效应的数值模拟.地球物理学报, 51(5):1411-1421. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200805015.htm 陈运泰, 杨智娴, 张勇, 刘超, 2013.从汶川地震到芦山地震.中国科学:地球科学, 43(6):1064-1072. http://www.cnki.com.cn/Article/CJFDTOTAL-ZSCD201701022.htm 董培育, 程惠红, 曾祥方, 石耀霖, 2013.四川芦山MS 7.0级地震导致周边断层的应力变化.科技导报 (北京), 31(12):19-22. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=kjdb201312015&dbname=CJFD&dbcode=CJFQ 邓起东, 陈社发, 赵小麟, 1994.龙门山及其邻区的构造和地震活动及动力学.地震地质, 16(4):389-403. http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ404.013.htm 邓起东, 陈桂华, 朱艾斓, 2011.关于2008年汶川MS 8.0地震震源断裂破裂机制几个问题的讨论.中国科学:地球科学, 41(11):1559-1576. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=jdxk201111002&dbname=CJFD&dbcode=CJFQ 房立华, 吴建平, 王未来等, 2013.四川芦山MS 7.0级地震及其余震序列重定位.科学通报, 58(20):901-1909. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=kxtb201320004&dbname=CJFD&dbcode=CJFQ 李传友, 宋方敏, 冉勇康, 2004.龙门山断裂带北段晚第四纪活动性讨论.地震地质, 26(2):248-258. http://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ200402006.htm 李玉江, 陈连旺, 陆远忠, 詹自敏, 2013a.汶川地震的发生对周围断层稳定性影响的数值模拟.地球科学——中国地质大学学报, 38(2):398-410. http://www.cnki.com.cn/Article/CJFDTOTAL-DQKX201302023.htm 李玉江, 陈连旺, 杨树新, 2013b.基于应变能变化的芦山强震同震效应的数值模拟.震灾防御技术, 8(4):361-369. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20130403&journal_id=zzfyjs 李玉江, 陈连旺, 刘少峰, 杨树新, 荆燕, 2014.芦山地震的发生对周围断层影响的数值模拟.地球学报, 35(5):627-634. doi: 10.3975/cagsb.2014.05.13 刘杰, 易桂喜, 张致伟等, 2013.2013年4月20日四川芦山M 7.0级地震介绍.地球物理学报, 56(4):1404-1407. doi: 10.6038/cjg20130434 缪淼, 朱守彪, 2013.2013年芦山MS7.0地震产生的静态库仑应力变化及其对余震空间分布的影响.地震学报, 35(5):619-631. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDW201310008003.htm 单斌, 熊熊, 郑勇, 刁法启, 2009.2008年5月12日MW7.9汶川地震导致的周边断层应力变化.中国科学 (D辑), 39(5):537-545. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=jdxk200905001&dbname=CJFD&dbcode=CJFQ 单斌, 熊熊, 郑勇等, 2013.2013年芦山地震导致的周边断层应力变化及其与2008年汶川地震的关系.中国科学:地球科学, 43(6):1002-1009. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK201306008.htm 邵志刚, 周龙泉, 蒋长胜, 马宏生, 张浪平, 2010.2008年汶川MS 8.0地震对周边断层地震活动的影响.地球物理学报, 53(8):1784-1795. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dqwx201008005&dbname=CJFD&dbcode=CJFQ 石耀霖, 曹建玲, 2008.中国大陆岩石圈等效粘滞系数的计算和讨论.地学前缘, 15(3):82-95. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY200803006.htm 史翔, 冉勇康, 陈立春, 王虎, 刘瑞春, 2009.龙门山中央断裂北川-邓家一带古地震初步研究.第四纪研究, 29(3):494-501. http://www.cnki.com.cn/Article/CJFDTOTAL-DSJJ200903010.htm 万永革, 吴忠良, 周公威等, 2000.几次复杂地震中不同破裂事件之间的"应力触发"问题.地震学报, 22(6):568-576. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXB200006001.htm 万永革, 沈正康, 盛书中, 徐晓枫, 2009.2008年汶川大地震对周围断层的影响.地震学报, 31(2):128-139. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXB200902002.htm 王椿镛, 楼海, 吕智勇等, 2008.青藏高原东部地壳上地幔S波速度结构——下地壳流的深部环.中国科学 (D辑), 38(1):22-32. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200801003.htm 王辉, 刘杰, 石耀霖, 张怀, 张国民, 2008.鲜水河断裂带强震相互作用的动力学模拟研究.中国科学 (D辑), 38(7):808-818. http://www.cnki.com.cn/Article/CJFDTOTAL-JDXK200807003.htm 王敏中, 王炜, 武际可, 2011.弹性力学教程.北京:北京大学出版社. 吴建平, 明跃红, 王椿镛, 2006.川滇地区速度结构的区域地震波形反演研究.地球物理学报, 49(5):1369-1376. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200605015.htm 徐晶, 邵志刚, 马宏生, 张浪平, 2013.鲜水河断裂带库仑应力演化与强震间关系.地球物理学报, 56(4):1146-1158. doi: 10.6038/cjg20130410 徐锡伟, 闻学泽, 叶建青等, 2008.汶川MS 8.0地震地表破裂带及其发震构造.地震地质, 30(3):597-629. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=dzdz200803003&dbname=CJFD&dbcode=CJFQ 徐锡伟, 陈桂华, 于贵华等, 2013a.芦山地震发震构造及其与汶川地震关系讨论.地学前缘, 20(3):11-20. http://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201303002.htm 徐锡伟, 闻学泽, 韩竹军等, 2013b.四川芦山7.0级强震:一次典型的盲逆断层型地震.科学通报, 58(20):1887-1893. http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201320002.htm 许冲, 徐锡伟, 2014.2013年芦山地震滑坡空间分布样式对盲逆断层构造的反映.科学通报, 59(11):979-986. http://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201411005.htm 杨强, 党亚民, 2010.利用GPS速度场估算青藏高原地壳韧性层等效粘滞系数分布的研究.测绘学报, 39(5):497-502. http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201005012.htm 杨兴悦, 陈连旺, 杨立明, 李玉江, 谭佩, 2013.巴颜喀拉块体强震动力学过程数值模拟.地震学报, 35(3):304-314. http://cdmd.cnki.com.cn/Article/CDMD-85403-1013107357.htm 张培震, 徐锡伟, 闻学泽等, 2008.2008年汶川8.0级地震发震断裂的滑动速率、复发周期和构造成因.地球物理学报, 51(4):1066-1073. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX200804017.htm 周光泉, 刘孝敏, 1996.粘弹性理论.合肥:中国科学技术大学出版社. Freed A. M., 2005. Earthquake triggering by static, dynamic, and postseismic stress transfer. Annual Review of Earth and Planetary Sciences, 33:335-367. doi: 10.1146/annurev.earth.33.092203.122505 Harris R. A., 1998. Introduction to special section:Stress triggers, stress shadows, and implications for seismic hazard. Journal of Geophysical Research, 103(B10):24347-24358. doi: 10.1029/98JB01576 Jia K., Zhou S. Y., Zhuang J. C., et al., 2014. Possibility of the independence between the 2013 Lushan earthquake and the 2008 Wenchuan earthquake on Longmen Shan fault, Sichuan, China. Seismological Research Letters, 85(1):60-67. doi: 10.1785/0220130115 Liu M., Luo G., Wang H., 2014. The 2013 Lushan earthquake in China tests hazard assessments. Seismological Research Letters, 85(1):40-43. doi: 10.1785/0220130117 Luo G., Liu M., 2010. Stress evolution and fault interactions before and after the 2008 Great Wenchuan earthquake. Tectonophysics, 491:127-140. doi: 10.1016/j.tecto.2009.12.019 Parsons T., Ji C., Kirby E., 2008. Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan basin. Nature, 454(7203):509-510. doi: 10.1038/nature07177 Ran K. Y., Chen W. S., Xu X. W., et al., 2013. Paleoseismic events and recurrence interval along the Beichuan-Yingxiu fault of Longmenshan fault zone, Yingxiu, Sichuan, China. Tectonophysics, 584:81-90. doi: 10.1016/j.tecto.2012.07.013 Toda S., Lin J., Meghraousi M., et al. 2008. 12 May 2008 M=7.9 Wenchuan, China, earthquake calculated to increase failure stress and seismicity rate on three major fault systems. Geophysical Research Letters, 35(17). https://darchive.mblwhoilibrary.org/bitstream/handle/1912/3372/2008GL034903.pdf;sequence=1 Wang Y. Z., Wang F., Wang M., et al., 2014. Coulomb stress change and evolution induced by the 2008 Wenchuan earthquake and its delayed triggering of the 2013 MW6.6 Lushan earthquake. Seismological Research Letters, 85(1):52-59. doi: 10.1785/0220130111 Wells D. L., Coppersmith K. J., 1994. New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement. Bulletin of the Seismological Society of America, 84(4):974-1002. http://www.academia.edu/1151208/New_Empirical_Relationships_among_Magnitude_Rupture_Length_Rupture_Width_Rupture_Area_and_Surface_Displacement 期刊类型引用(3)
1. 熊政辉，曲哲，宁超列，卜春尧，尹建明，薄涛. 我国地震巨灾保险试点实践对比研究. 灾害学. 2025(01): 48-52+159 . 
百度学术2. 李云仙，子芝兰. 基于AHP-FCE的地震巨灾保险运行效果评价——以大理州政策性农房地震巨灾保险为例. 保险理论与实践. 2024(10): 33-55 . 百度学术3. 景冰冰，施唯. 2023年土耳其M_S7.8地震灾区损失与保险研究. 震灾防御技术. 2023(03): 495-504 . 本站查看其他类型引用(6)
-
下载:
点击查看大图
计量
- 文章访问数: 137
- HTML全文浏览量: 46
- PDF下载量: 18
- 被引次数: 9
