Integration Optimization and Effect Analysis of Earthquake Flow Observation Shield
-
摘要: 为提升地震流动观测质量，适应恶劣条件下野外长期工作，研制集传感器、供电、通讯网络三大模块于一体的流动观测仪器整体防护罩。利用2套GL-PS2一体化短周期地震仪进行对比观测，一套安装于南京台室内基岩观测墩，另一套配备防护罩浅层掩埋于户外。结果表明新型研制的观测仪整体集成防护罩具有以下优点：①保温性能好，内部环境日温差小于3℃，月温度变化小于5℃。②防水性好，经历多次暴雨大雪天气后，防护罩内部干燥，无水滴或结露现象。③台基噪声略优于室内观测墩。④2台设备相关系数垂直向优于水平向，夜晚优于白天。⑤地震记录波形相关系数均在0.6以上，震级越大，相关性越高，且信噪比差距较小。Abstract: In order to improve the earthquake flow observation quality, a new seismograph shield is developed. Integrated optimization of sensor, power supply and communication network is able to achieve overall protection, which is adapted to long-term field work under severe conditions. Comparative observation with two GL-PS2 integrated short period seismograph has been done, one is installed on the bedrock observation pier of Nanjing seismic station, the other is shallowly buried out of doors with integrated optimization of shield. The results show that the newly developed observer Shield and integrated protection have the following advantages: ①good thermal insulation, internal environment daily temperature difference less than 3℃, monthly temperature change less than 5℃. ②good water resistance, and the interior of the protective cover is dry without water drop or condensation after heavy rain and heavy snow. ③the noise with integrated protection is slightly better than that of indoor observation pier. ④the correlation coefficient of the two is better in the vertical direction than in the horizontal direction, and better in the night than in the day. ⑤the correlation coefficients of the waveforms of seismic records are all above 0.6, the larger the magnitude, the higher the correlation, and the difference between the signal to noise ratio is small.
-
引言
跨断层形变测量是通过观测布设在断层两侧控制点的相对关系，以监测断层（上、下盘）的相对运动。由于控制点布设的位置，跨断层测量能更容易、更直接地观测到断层的相对运动情况。很多学者对跨断层测量进行了研究。张超（1981）将断层上、下盘作为刚性块体，从几何学的角度进行了研究，利用基线观测数据、水准观测数据得到断层三维运动情况。薄万举等（1998）在此基础上，对八宝山断裂带上的跨断层场地数据进行了处理，对该地区的断层运动情况进行分析，发现在一定区域内跨断层测量对异常形变监测和震情预测研究有重要意义。张兴飞（1992）针对跨断层形变测量的特点，构建了基于混合回归理论的动态数据处理模型，该模型对地震前兆有较好的反映。焦青等（1998）利用华北和南北地震带跨断层流动水准和定点水准观测资料，研究了断层垂直位移速率异常的空间分布以及震源区的逼近情况，发现通过研究断层垂直位移速率预测地震是一种可行方法。张希等（2013）将自适应控制理论引入断层流动形变分析中，通过构建慢时变前兆系统结构，对5个典型场地资料进行分析，发现其异常识别、异常分类的判定较原始曲线更为直观，与以往的形变指标相比也有所改进。张晶等（2011）从断层受力运动方向入手，对断层在一定外力作用下的运动情况进行了不同的组合，形成断层活动协调比理论，分析了断层活动性质和应变积累状态，对地震中长期预测有一定的参考意义。
基于跨断层测量的研究基本上都是通过计算断层的三维运动情况，然后对其运动规律进行发掘，发现异常，这些方法的最终目的都是为了寻找或确定某条断层或某个场地附近断层的运动变化机理，进而开展地震预测预报工作。本文在该领域专家研究的基础上，深入研究垂直形变趋势累计率（陈兵等，2000），对其反映的物理机制以及区域应变场在中强震发生前后的变化规律进行了进一步研究，并得到了一些有用的结论。
1. 基本原理
垂直形变趋势累计率是通过垂直形变时变曲线建立的一种可较明确地量化趋势变化的判定方法。此方法以断层总体速率变化绝对量为基础，再根据断层总的垂直活动量判断断层运动趋势，而且具有表征孕震应力场时空动态演化的物理意义。
首先通过公式（1）获得断层的相对运动速率（上盘相对下盘，上升为正）：
$${v_i} = \frac{{{y_{i + 1}} - {y_i}}}{{{t_{i + 1}} - {t_i}}}$$ (1) i表示观测值序列号，ti为观测序列所对应的时间，yi表示ti时刻上盘相对下盘观测的高差值，vi表示ti时刻相邻2期断层的相对垂直形变速率。
为了排除干扰因素，更真实地计算断层上、下盘的相对运动速度，本文在计算断层相对运动速率时，首先将不同年份（序号为j）相同月份（序号为k）的速率值vj, k取平均得到vk，然后再在此基础上，用每一年每一期的速率值（vj, k）减去每一期对应的平均值vk获取vj, k'，称之为总体性去周期，其中j为年。具体见式（2）、（3）：
$${\bar v_k} = \frac{1}{m}\sum\limits_{j = 1}^m {{v_{j,k}}} $$ (2) 其中m为参与计算测线的总观测年数。
$${v'_{j,k}} = {v_{j,k}} - {\bar v_k}$$ (3) 该想法启发于大地水准面的确定。与相邻年同期观测数据做差或较短观测期间不同时间段（年）中相同期（年中月）取平均值获取相对运动速率的传统方法相比，这样做能够更好地消除因环境等因素引起的长短周期影响。
总体去周期后，那么以起算时刻（序号为s，场地起始观测年份 < 起算年份 < 最后观测年份）至计算时刻（序号为m、k0所对应年份和相应的月份）这一时间段断层总体垂直速率变化量va, k0如下：
$${v_{{\rm{a}},{k_0}}} = \sum\limits_{j = s}^{m - 1} {\sum\limits_{k = 1}^{12} {{{v'}_{j,k}} + \sum\limits_{k = 1}^{{k_0}} {{{v'}_{m,k}}} } } $$ (4) 其中，一般情况下s=m-10。
将每期垂直形变速率取绝对值求和，得到（5）式，vd, k0为此时间段垂直形变变化幅度。
$${v_{{\rm{d}},{k_0}}} = \sum\limits_{j = s}^{m - 1} {\sum\limits_{k = 1}^{12} {\left| {{{v'}_{j,k}}} \right|} } + \sum\limits_{k = 1}^{{k_0}} {\left| {{{v'}_{m,k}}} \right|} $$ (5) 依据（4）、（5）式，计算垂直形变趋势性的量化标准，即（6）式，垂直形变趋势累计率。
$${D_{\rm{c}}} = \frac{{{v_{{\rm{a}},{k_0}}}}}{{{v_{{\rm{d}},{k_0}}}}} = \frac{{\sum\limits_{j = s}^{m - 1} {\sum\limits_{k = 1}^{12} {{{v'}_{j,k}} + \sum\limits_{k = 1}^{{k_0}} {{{v'}_{m,k}}} } } }}{{\sum\limits_{j = s}^{m - 1} {\sum\limits_{k = 1}^{12} {\left| {{{v'}_{j,k}}} \right|} } + \sum\limits_{k = 1}^{{k_0}} {\left| {{{v'}_{m,k}}} \right|} }}$$ (6) 如果断层运动速率vi为恒值，那么va, k0 = vd, k0 =0，，即Dc为0/0型。为使Dc与实际物理情况相符，此情况下规定Dc=1。
由（6）式可以发现Dc∈[-1, 1]，那么确定垂直形变趋势起算时间后，如果Dc > 0，则断层运动趋势与起算时刻的运动趋势一致，且断层活动性为压性；如果Dc < 0，则断层运动趋势与起算时刻的运动趋势相反，且断层活动性为张性；如果Dc=0，说明断层恢复到起算时的状态，即垂直形变无变化。|Dc|越接近1，则可以认为断层形变异常强度越高；如果|Dc|越接近0则反之。在分析过程中需要注意的是，获取的Dc值是基于一定的起止时间段的，因此其仅反映此时间段内断层的活动性，而且可能会与一些块体总体（长趋势）活动性不尽相同。
另外，速率绝对幅度的大小及积累时间的长短与相应的地震强度和复发周期有关。强震孕育到一定程度后，会出现闭锁、微破裂、预滑移等现象，必然对断层形变的大小和方向产生影响，从而反映到Dc值的异常变化。
2. 跨断层观测场地情况
基于以上原理，本文以川滇地区为实验区，利用布设在重要断裂带上的跨断层场地的水准观测数据开展实验，并进行了分析。
以下为川滇地区跨断层流动场地水准资料观测情况。各个场地、各条测线的观测情况如表 1，川滇地区跨断层流动观测场地分布情况如图 1。
表 1 川滇地区在测流动跨断层场地情况表Table 1. List of the across fault monitoring in the Sichuan-Yunnan region场地名称 测段 所跨断层 开始观测时间 复测周期 侏倭 AB、AC 鲜水河断裂 1980.4 6期/年 格篓坝子 AB、AC 鲜水河断裂 1974.8 6期/年 虚墟 BA、CA 鲜水河断裂 1973.2 6期/年 沟普 BA、CA 鲜水河断裂 1983.5 6期/年 龙灯坝 BA、CA、DA 鲜水河断裂 1985.5 6期/年 老乾宁 13、53、54 鲜水河断裂 1979.9 6期/年 折多塘 BE 鲜水河断裂 2013.11 6期/年 团结 12、13、43 鲜水河断裂 1986.5 6期/年 棉蟹支线 13、23、57（停）、78（停） 安宁河则木河断裂 1990.3 6期/年 安顺场 BA、CA 鲜水河断裂 1977.8 6期/年 冕宁支线 13、23、54 安宁河则木河断裂 1981.9 6期/年 宁南支线 54 安宁河则木河断裂 1972.6 6期/年 汤家坪 AB、AC 安宁河则木河断裂 1981.9 6期/年 汤家坪新 BC 安宁河-则木河断裂 2012.9 6期/年 尔乌 31 安宁河则木河断裂 1981.9 6期/年 蒲江 BC 蒲江-新津-成都-德阳断裂 1977.7 6期/年 双河 12 龙门山断裂 1977.7 6期/年 灌县 34 龙门山断裂-灌县断裂 1978.2 6期/年 七盘沟 CA 龙门山断裂 1986.5 6期/年 荞窝 BC 安宁河-则木河断裂 2010 6期/年 西昌 BA、BC 安宁河-则木河断裂 1974.6 6期/年 宝兴 BX 盐井-五龙断裂 2008.8 6期/年 映秀 BC 龙门山断裂-映秀断裂 1985.5 6期/年 剑川 21、31 剑川断裂 1982.2 12期/年 丽江 1x-2x、43 丽江断裂 1982.2 12期/年 永胜 14、1-2x、3-2x 程海断裂 1982.2 12期/年 楚雄 1x-2、32、35 楚雄断裂 1982.2 12期/年 羊街 1x-x5 小江断裂-西支断裂 1982.2 12期/年 宜良 2x-1x、2x-3x、4-3x 小江断裂东支断裂带南羊街断裂 1997.7 12期/年 峨山 12、34 曲江断裂 1982.3 12期/年 通海 12、32、34 曲江断裂 1982.3 12期/年 建水 21、31 越屏-建水断裂 1982.3 12期/年 石屏 21、31 石屏-建水断裂 1982.3 12期/年 下关① 12、32 洱海断裂 1982.2 12期/年 ① 说明：例如下关测段12，即表示从上盘控制点1到下盘控制点2；2x-1x表示上盘2号新控制点到下盘新1号控制点，x代表新点。 
图 1 川滇地区跨断层流动观测场地及断层分布图Figure 1. The distribution of faults and across fault monitoring sites in the Sichuan-Yunnan region上述场地中部分测线因场地环境改变或观测墩（点）遭破坏等情况而停测或更改观测路线的测段，以及未跨过断层的测段并未列出。另外由于一些新建场地的观测时间距现在较近，故这些场地的相关信息也未列出。
川滇地区跨断层观测场地多数建于20世纪80年代，主要分布在鲜水河断裂带、龙门山断裂带、安宁河-则木河断裂带以及小江断裂带等几大断裂带上，其场地分布见图 1。场地位置用五角星标注。
3. 实验分析
川滇地区的观测场地主要位于南北地震带上，而南北地震带是我国地震高发地区，拥有丰富的震例（刘辛中等，2015），因此作者选取该地区地震为样本进行分析。另外为弥补云南、四川跨断层场地观测周期不同以及因外界环境等因素造成的数据不全对计算的影响，在数据预处理过程中对数据进行了3次样条插值处理。其中，为获取地震前后川滇地区垂直形变趋势较为真实的变化情况，尽量选择发震时间前后各个场地都有观测数据的时间作为实验时间，对受外界条件影响或在核实异常过程中认为存在质量问题的观测数据不予采用，也尽量避免因插值造成的不必要误差。
由于篇幅原因，基于第2章所列观测场地的观测资料仅对汶川地震、芦山地震前后Dc值变化所反映的规律分别进行讨论，其中起算时间分别选为1998、2003年。
图 2中，（a）与（d）分别为汶川地震发生前4个月、震后2个月的垂直形变趋势累计率，（b）与（c）为地震发生前2个月、震后1个月的垂直形变趋势累计率。可以看出，汶川地震发生在垂直形变趋势变化较大的区域附近，而地震前后，发震地附近的垂直形变趋势发生较大改变（矩形框示意范围），反映其附近断层垂直形变异常强度改变较大。地震发生前2个月形变趋势表现为运动平缓，即形变异常强度变小，而震后则逐渐恢复到与震前相同的情况。作者认为，这种情况说明发震地附近震前断层的“正常”运动趋势受阻，导致断层附近出现更为明显的闭锁现象，为地震的发生提供了必要条件。相关文献（赵静等，2012）也给出了类似的观点，进一步说明跨断层水准能够较有效地反映地震前兆。实验过程中，作者发现在汶川地震发生后2个月（图 2（d））川滇地区的垂直形变异常强度与汶川地震发震前（图 2（b））的垂直形变异常强度极为相似，而通过查询中国地震台网中心发布的地震目录可以发现，在汶川地震发生2个多月后，即2008年7月24日，汶川主震震中东北方约250km处的广元发生了6.0级强余震（郑勇等，2009；温瑞智等，2009）。需要说明的是，川滇地区跨断层观测数据的获取时间主要集中在每个观测月的中上旬，故图 2（d）显示的为广元地震发生前该地区的垂直形变异常强度，这从另一个方面（形变学）说明跨断层手段所反映的垂直形变异常强度具有较好的映震能力。
图 2 汶川地震（2008.5.20）前后垂直形变速率趋势变化Figure 2. The vertical deformation trend variation before and after the Wenchuan earthquake图 3中，（a）、（b）分别为该地区发震前6个月和1个月的垂直形变趋势图，（c）、（d）为震后1个月和2个月的垂直形变趋势图。从图中可以发现芦山地震的发震地同样处于垂直形变趋势变化较大（即垂直形变异常强度值变化较大，主要集中在震中东南方向，矩形框示意范围）的地区附近。川滇地区2012年9月份与2013年5月、6月的垂直形变趋势累计率几乎保持一致，而该地区2013年3月较2012年9月的垂直形变累计率发生了较大改变，即断层垂直形变异常强度值改变较大，随即2014年4月20日在该区域便发生MS 6.7芦山地震。
图 3 芦山地震（2013.4.20）前后垂直形变速率趋势变化Figure 3. The vertical deformation trend variation before and after the Lushan earthquake总的来说，垂直形变趋势累计率在以下2方面有较好的效果：① 从地震发生的位置上来看，发震地一般都处于断层形变异常强度值改变较大的地区，地震发生前断层附近的垂直形变异常强度有突然变缓或加剧的趋势；② 从地震前后应变能变化的角度研究，地震发生前应变能大量累积，地震发生后得到释放，但还处于平衡自适应阶段，震后一段时间应变能达到平衡，使其与震前相当，较好地反映了地震孕震、发震以及震后平静期的应变能变化过程（许昭永等，2015）。
4. 讨论
通过分析川滇地区的震例，得到了几方面认识：从时间角度看，地震发生前后，川滇地区跨断层场地附近垂直形变异常强度几乎维持不变；从空间位置角度看，地震（中强震）主要发生在垂直形变异常强度较大地区或者其边缘及附近；从震级角度来看，较大震级的地震在震前的异常形变强度值变化比较小震级地震更为明显，但目前缺少更为确切的依据，需要开展进一步研究。总体来说，垂直形变趋势累计率在时间和空间上对中短期地震预测预报能够发挥不错的效果，但对于短临预报，由于数据采集频率以及跨断层场地分布情况的影响，该手段的有效性有待进一步验证。另外，该方法需要足够的监测场地密度，否则仅从单个跨断层场地的垂直形变趋势变化速率以及趋势变化幅度分析，会因缺少断层运动的整体性判断而更加难以获得真实的空间异常变化特征，因此对地震预测预报的帮助可能较小。地震发生过程是应变能积累-释放-再积累的过程，因此在选用垂直形变趋势累计率方法时，起算时间的选择也是不可忽略的因素之一。
致谢: 感谢2位审稿人对本文的悉心指导，提出的建议对本文质量提升有很大帮助。本文所使用的跨断层观测数据为一线观测人辛苦所得以及兄弟单位慷慨提供，在此一并感谢。 -
图 1 防护罩集成优化示意图
注： 1-保温层；2-防护罩；3-4G路由器；4-O型密封圈；5-太阳能板葛兰头；6-水平仪；7-覆土；8-GPS、GPRS、WIFI、葛兰头；9-太阳能控制器；10-锂电池；11-铝制隔板；12-地震仪；13-太阳能板
Figure 1. Integrated optimization schematic of shield
表 1 2套设备标定结果
Table 1. Calibration results of two seismographs
通道 周期/s（变化率/%） 阻尼（变化率/%） 电压灵敏度/V/（m·s-1）（变化率/%） 基岩_ud 2.1084（5.4185) 0.716816（1.3885） 998.34（0.17） 基岩_ew 2.0783（3.9128） 0.703904（0.4379） 1000.93（0.09） 基岩_ns 2.0415（2.0763） 0.702160（0.6846） 1004.08（0.41） 集成_ud 1.9972（0.1402） 0.691597（2.1786） 1001.66（0.17） 集成_ew 2.0229（1.1467） 0.704364（0.3728） 997.49（0.25） 集成_ns 2.0457（2.2867） 0.701876（0.7248） 999.39（0.06） 
下载: 导出CSV
表 2 2台设备平均台基噪声
Table 2. Average noises of two seismographs
地震计 噪声平均值×10-8/m·s-1 平均噪声counts值（p-p） 噪声等级 基岩_ud 5.80101 1101 Ⅱ 基岩_ew 7.80126 1481 Ⅱ 基岩_ns 9.60356 1823 Ⅱ 集成_ud 5.08683 965 Ⅱ 集成_ew 7.43162 1411 Ⅱ 集成_ns 7.85839 1492 Ⅱ
下载: 导出CSV
表 3 不同地震相关系数统计
Table 3. Statistics of correlation coefficient with different earthquakes
震级M 地震个数 相关系数R 0.5—0.6 0.6—0.7 0.7—0.8 0.8—0.9 0.9以上 1.0—2.0 UD — 2 1 — — EW 1 2 — — — NS — 3 — — — 2.0—3.0 UD — 1 6 — — EW — 3 4 — — NS — 2 5 — — 3.0—4.0 UD — 2 5 2 — EW — 4 4 1 — NS — 3 5 1 — 4.0—5.0 UD — — 3 7 — EW — — 5 5 — NS — — 6 4 — 5.0—6.0 UD — — — 5 — EW — — 2 3 — NS — — 2 3 — 6.0—7.0 UD — — — — 3 EW — — — 2 1 NS — — — 2 1 7.0以上 UD — — — — 2 EW — — — — 2 NS — — — 1 1
下载: 导出CSV
表 4 地震波形信噪比偏差统计
Table 4. Deviation statistics of signal to noise ratio of seismic waveform
震级M 地震个数 SNR1-SNR2 -0.5—0 0—0.5 0.5—1 1—1.5 1.5以上 1.0—2.0 UD 3 — — — — EW 2 1 — — — NS 3 — — — — 2.0—3.0 UD 4 2 1 — — EW 3 4 — — — NS 4 3 — — — 3.0—4.0 UD — 2 5 2 — EW — 5 4 — — NS — 6 3 — — 4.0—5.0 UD — 4 5 1 EW — 1 5 4 — NS — 3 3 4 — 5.0—6.0 UD — 2 3 — EW — 1 2 3 — NS — 2 1 2 — 6.0—7.0 UD — — — 3 — EW — — 1 2 — NS — — 1 2 — 7.0以上 UD — — — 2 1 EW — — — 3 — NS — — — 2 —
下载: 导出CSV
-
陈颙, 周华伟, 葛洪魁, 2005. 华北地震台阵探测计划. 大地测量与地球动力学, 25(4): 1-5. https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB200504000.htm 葛洪魁, 陈海潮, 欧阳飚等, 2013. 流动地震观测背景噪声的台基响应. 地球物理学报, 56(3): 857-868. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201303016.htm 花鑫升, 石金虎, 谭雅丽等, 2018. 浅层地震勘探资料揭示汤东断裂特征. 震灾防御技术, 13(2): 276-283. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20180203&journal_id=zzfyjs 姜旭东, 杨建思, 徐志强等, 2009. 应对巨大地震的应急流动观测系统. 地震地磁观测与研究, 30(5): 59-65. https://www.cnki.com.cn/Article/CJFDTOTAL-DZGJ200905010.htm 李世成, 崔建文, 张潜等, 2011.2009年姚安Ms6.0强震动应急观测. 地震研究, 34(3): 322-330. https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB201102004.htm 李世成, 林国良, 崔建文, 2017. 云南强震动流动观测发展及相关问题讨论. 地震研究, 40(4): 629-637. doi: 10.3969/j.issn.1000-0666.2017.04.015 任叶飞, 温瑞智, 周宝峰, 2011. 强震动记录编目——以汶川余震流动观测记录为例. 震灾防御技术, 6(4): 406-415. doi: 10.3969/j.issn.1673-5722.2011.04.006 孙宏志, 2016. 地震计智能恒温系统的设计与实现. 地震工程学报, 38(6): 1004-1009. https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201606024.htm 孙路强, 柳艳丽, 刘文兵等, 2014. 沉积层地震观测噪声的对比试验. 震灾防御技术, 9(4): 932-936. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20140422&journal_id=zzfyjs 王恒知, 俞言祥, 何斌等, 2018. 高邮-宝应4.9级地震地面运动特征分析. 地震地质, 40(6): 1349-1361. https://www.cnki.com.cn/Article/CJFDTOTAL-DZDZ201806011.htm 温瑞智, 周正华, 李小军等, 2009. 汶川Ms8.0地震的强余震流动观测. 地震学报, 31(2): 219-225. https://www.cnki.com.cn/Article/CJFDTOTAL-DLXB201910010.htm 杨建思, 2008. 汶川地震所见的零星感受与启示. 国际地震动态, (7): 40-45. doi: 10.3969/j.issn.0253-4975.2008.07.007 张云鹏, 李孝宾, 王伟涛等, 2017. 云南宾川地震信号发射台的流动观测数据服务系统及数据质量评估. 地震研究, 40(4): 525-533. doi: 10.3969/j.issn.1000-0666.2017.04.003 郅红魁, 程楠, 郭慧等, 2019. 地震计恒温低湿智能监控仪研制. 地震地磁观测与研究, 40(1): 128-133. https://www.cnki.com.cn/Article/CJFDTOTAL-DZGJ201901019.htm Burdick S., Li C., Martynov V., et al., 2008. Upper mantle heterogeneity beneath North America from travel time tomography with global and USArray transportable array data. Seismological Research Letters, 79(3): 384-392. doi: 10.1785/gssrl.79.3.384 Yang Y. J., Ritzwoller M. H., Lin F. C., et al., 2008. Structure of the crust and uppermost mantle beneath the western United States revealed by ambient noise and earthquake tomography. Journal of Geophysical Research: Solid Earth, 130(B12): B12310. doi: 10.1029/2008JB005833/full -
下载:
点击查看大图
计量
- 文章访问数: 130
- HTML全文浏览量: 41
- PDF下载量: 4
- 被引次数: 0
