基于形变仪器同震形变波衰减规律的研究1

龚丽文 陈丽娟 刘 琦 郑许东



基于形变仪器同震形变波衰减规律的研究1

龚丽文1）陈丽娟1）刘 琦2）郑许东1）

1）重庆市地震局，重庆 401147 2）中国地震局地震预测研究所，北京 100036

水管仪 垂直摆 同震形变波 衰减因子 振幅-震级关系

引言

近年来，地震前兆观测已经积累了大量的观测资料，对于观测资料的应用和预报成果的产出迫在眉睫。形变台网主要记录形变演化，在日常观测中发现，同一区域不同仪器记录的形变信息可能有较大差异，特别是位于同一观测硐室不同倾斜仪所记录到的倾斜趋势完全不一致，这对数据的物理解释和应用造成一定困难，但其固体潮记录和同震形变波记录都明显优于同台的其它应变仪器。对于倾斜仪的同震记录有不同角度的研究，杨玲英等（2012）和方宏芳等（2010）针对同台不同仪器的同震响应作了对比分析，认为摆式仪器的同震响应幅度大于同台水管仪的原因是前者的阻尼远小于后者；张利兵等（2013）对比不同台同一类仪器的同震响应，认为其同震响应具有一定差异性，可能受台站构造环境和仪器参数影响；邱永平（2011）对比分析了同套仪器对不同巨大远震的同震响应；吴利军等（2016）利用体应变数据，统计并分析地震波最大振幅、震级及震中距的关系；吕品姬等（2010）和曹喜等（2014）分别研究了台站倾斜仪的映震能力，地倾斜类仪器一般倾向于对周边的中强地震具有短临异常指示意义。

同震形变波是受地震波激发、形变观测仪在多测点观测到的地震时倾斜、应变或应力的波动现象（牛安福等，2005），研究震时形变波动与地震震源的关系能为识别地震短临前兆及短临预测提供重要的方法途径（冯德益等，1993，1994）。倾斜仪记录到的同震形变波集中在面波，其主要参数包括面波延迟时间、最大变形幅度及同震持续时间等，它的震荡时间和响应幅度都与震级大小和震中距有一定关系（牛安福等，2006）。面波是由体波产生的派生波，能量强、传播远，具有很强的频散现象，其演化产生的转动分量比平动分量更富有高频信息（李宏男等，2001）。同震形变波能量的传递与演化对地震危险区的触发能力之间的关系仍需进一步研究，不同区域的特大远震所产生的同震形变波对区域应力场的影响幅度是判断震后区域危险性的重要依据，而观测仪器所记录到的最大振幅是量化同震形变波能量的最直接的依据（Uang等，1990）。因此，本文以黔江台倾斜仪记录到的442次同震形变波数据为基础，对比分析该台2套倾斜仪同震形变波的记录能力，统计同震形变波的振幅-震级（-）关系，为后续分析同震形变波对区域应力场影响能力提供一定依据。

1 台站简介及数据来源

黔江台地处黔江区仰头山风景区内（29.55°N，108.80°E），海拔1020m，该区域地处渝东南武陵山区，地貌受地质拼迭的控制，以低、中山为主，山脉走向多与构造线方向一致。台站位于NE向的彭水基地断裂和黔江断裂之间，距台站4km处有第四纪仰头山断层通过。1856年在上述断层附近发生6¼级的黔江小南海地震（图1（a））。黔江台2007年改建为形变综合台，主要形变仪器包括SS-Y伸缩仪、DSQ水管仪和VS垂直摆，3套仪器均安装在20世纪70年代人工开凿的山洞内，该洞距洞口50m处分岔为两支洞，夹角80°，NS向支洞长约29m，EW向支洞长约42m，仪器布设见图1（b）。山洞顶部覆盖厚度超过70m，旁覆厚度超过500m，洞内干燥，洞温15.1℃，保温效果好。山洞基岩为二叠系栖霞灰岩，岩石完整、致密，各方面均符合设计要求。虽然该台2套倾斜仪的工作原理不同（图1（c）），但是均能记录到明显的固体潮，同震记录丰富，数据连续性较好，数据质量高，便于分析结果的对比，提高分析结果的可信度。

图1 黔江台站及仪器概况

本文选取该台2套倾斜仪在2010—2015年内记录到的442次分布全球的地震（据月报和日志文件统计）数据，统计不同测项对不同地震响应的分量振幅。由于仪器采样率的限制（分采样）和测项方向的差异性，使得仪器记录到的同震最大振幅具有一定误差，但本文选择的地震样本数量足够多，在一定程度上减小了这种误差；而最大振幅是倾斜仪在同震形变波的激发下仪器的震动响应程度，是同震形变波能量的一个量度，并非地面倾斜的真实角度，不具有矢量性。因此，本文采用倾斜仪最大振幅模量来量化同震形变波的振幅大小，并用振幅模量之和来消减测项因方向记录的差异，即：

式中，MAX、MIN为同震最大、最小倾斜幅度。

本文通过对初始分量振幅作简单的预处理，最大程度降低误差，从而提高统计数据的可靠性。

图2 同震形变波波形图

2 同震形变波波形对比分析

刘琦等（2013）利用S变换对高采样率钻孔应变同震形变波形作了详细的研究，认为仪器的采样率是能否记录到完整波形的先决条件，10Hz采样率的形变仪器能够记录到比较完整的地震形变波波形。由于该台3套仪器均为分采样，不能记录到完整的波形，导致记录到的近震主要以短期的脉冲信号为主（图2（a）、（b）），部分同震形变波还伴随有阶跃信号，能量较集中，振幅较大。随着震中距的增加，地震波发生频散，能量不断衰减，同震振幅不断减小，面波与体波逐渐分离，同震形变波持续时间变长（图2（c）、（d））。当震中距达到一定距离后，形变仪器记录到巨大远震引起的震荡波（图2（c）、（d）），主要包括球型震荡和环型震荡（张致伟等，2008；杨跃文等，2010）。

鉴于垂直摆与水管仪的观测物理量一致，但由于观测原理和基线长度不同，导致仪器对同震形变波信号的记录能力存在差异，加上采样率为分采样，记录到的波形与真实波形有一定出入。因此，本文通过对比分析2套仪器同震形变波振幅，来研究同震形变波的衰减规律和仪器对高频信号的记录能力。

3 同震形变波振幅统计分析

一般情况下，定点形变仪器对100km以内的4级地震、200km以内的5级地震、300km以内的6级地震及全球的7、8级地震都有响应（吕品姬等，2010）。通过该台的日志记录和月报资料，统计出2010—2015年442次地震，其中倾斜类仪器（垂直摆434次，水管仪408次）明显多于应变类仪器（伸缩仪168次），可能是地震面波产生的转动分量比平动分量更富有高频成分，使得倾斜仪器更容易记录到同震波动（李宏男等，2001）。

同震幅度的大小受很多因素影响，包括震级、震中距、发震构造、震源机制、台址环境及仪器的灵敏度等，其中震级为主要因素，震级越大，释放的能量越多，同震幅度越大（牛安福等，2005）。因此，本文通过量取2套倾斜仪所有同震记录的最大振幅，按仪器分类统计，并分别作振幅-震级图，结果如图3所示。由图可知：①2套仪器均能记录到周边4级以上同震，其中垂直摆对中小地震的记震能力强于水管仪；②对中强地震的记录，垂直摆振幅分布比较离散，水管仪的同震记录更稳定；③随着震级的增大，2套仪器的振幅均呈指数增长，这种关系与震级-能量的关系趋于一致。

图3 同震响应幅度对比图

对于同一同震形变波信号，2套仪器记录具有一定的差异性，除采样率可能在一定程度上影响观测结果外，还有可能是因为同震响应的最大幅度与台站仪器的固有频率有较强的关系，垂直摆自振周期较小，能捕捉到更多的高频信息，而水管仪内部液体的性质在一定程度上决定了该仪器相对较稳定的特性。

4 同震形变波的振幅-震级关系

同震形变波的振幅与地震传播的能量成正相关，而随着传播距离的增加，地震波不断发生散射、衰减，使其能量不断衰减，其衰减规律不仅与震源参数（震级和震中距）有关，还与传播介质参数（传播途径及传播介质的力学性质）和观测场地环境（观测环境和构造背景）等有关（公茂盛等，2003）。由于2套仪器处于同一观测硐室，其观测场地环境相同。不同地震的震源参数和传播介质参数差异性较大，且很难量化处理。但仪器记录地震的分布相对较集中，主要分布在环太平洋地震带上，按地震的分布特征将其分成6个研究区域，包括中国川滇地区及邻区、中国东海-中国台湾-菲律宾地区、日本岛弧地区、苏门答腊-印度尼西亚地区、所罗门群岛地区以及阿根廷-秘鲁-智利地区；位于同一个区域的地震，其构造环境相似，震中距、传播途径及传播介质参数相近，即仪器的同震幅度理论上只与震级大小有关，从而简化了研究对象，便于进一步研究振幅-震级关系。

图4 不同地震对应倾斜仪的幅度分布图

表1 不同研究区域振幅-震级关系表

续表

注：为震中距；i为最小震级；为地震记录数；c表示垂直摆；s表示水管仪。

图5 不同研究区域振幅-震级关系图（一）

图5 不同研究区域振幅-震级关系图（二）

表2 不同研究区域的拟合参数

注：、为拟合参数；R为拟合方差；为衰减因子；c表示垂直摆；s表示水管仪。

从图5可以看出，除中国川滇地区及邻区的-关系相对较离散外，其它地区的曲线拟合程度均较高，特别是在震中距大于2000km的区域地震，水管仪记录到-值的拟合程度均在0.8以上。随着震中距增加（智利地区除外），其值呈指数下降，而值呈缓慢的上升（表2），从图6（a）中可以看出，2套倾斜仪均显示ln与呈线性关系，且参数相近，即-关系式可以简化成单一参变量的指数关系式，该参变量能指示不同区域地震能量的衰减程度。因此，本文将其定义为衰减因子，原-关系式即可写成：

式中，为振幅，为震级，c和s分别是垂直摆和水管仪的同震形变波衰减因子。

将最大振幅和震级分别按仪器代入公式（2）、（3），分别计算出不同区域的衰减系数，如图6（b）显示：衰减系数随着震中距的增加呈现缓慢上升趋势，特别是印尼（＝4000km）至所罗门（＝6000km）处，其衰减系数增长速度较快，而印尼地区的地震震级较大，衰减较小，该区域发生巨大远震对我国地震危险区的影响较大；当震中距达到一定距离后（＞6000km），其衰减系数趋于3，这可能是地震波的传播途径发生变化导致的。从衰减系数随震中距的变化特征来看，垂直摆与水管仪的观测结果具有很强的一致性，特别是智利地震，由于其拟合程度较高，衰减系数都趋于3，更加说明拟合结果的可靠性。

图6 指数系数b与lna关系图（a）及衰减系数k与震中距关系图（b）

当震中距达到一定范围后，地震波的传播途径开始发生明显变化，特别是震中距远远大于地球半径的智利地区，地震产生的同震形变波在该地区不再单纯地通过地表传递（沿地表传播时间大于形变波延迟时间），而是通过地球内部不断折射到达台站，其能量衰减也明显小于沿地表传波，这是由于沿地球内部传播以体波传播为主，速度较快、频散弱、能量消耗小，从而使得观测区接收到的能量较大，可能会触发区域应力场的调整。

5 讨论

本文针对黔江台2套倾斜仪记录到的442次地震进行了详细的统计分析，按仪器分别对每次地震的振幅和形态做了详细的对比分析，得出以下结论：

（1）2套倾斜仪灵敏度高，能记录到丰富的同震形变波信息。但由于仪器观测原理不同，记录到的高频信息具有一定差异性，其中垂直摆由于自振周期较小，波谱信息记录较丰富，在提高采样率的前提下，有利于研究波形的时频特征；而水管仪由于基线较长，记录到的震动信号较稳定，振幅的拟合程度较高，有利于研究同震形变波的能量衰减规律。

（2）大量震例显示2套仪器记录到的同震形变波振幅具有一定的规律性，即振幅随着震级增加呈指数上升，而指数系数与ln呈显著的反比线性关系；随着震中距的增加，同震持续时间增长，-曲线拟合程度随之增高。2套仪器的统计结果趋于一致，可信度较高。

蔡骞，2013．伸缩仪和水管仪测试与实验研究．武汉：中国地震局地震研究所，1—53．

曹喜，董海龙，张彩艳等，2014．嘉峪关水管倾斜仪观测资料映震能力分析．地震地磁观测与研究，35（3—4）：174—177．

崔庆谷，杨星，杨跃文，2014．水管仪、垂直摆及水平摆高频段记录数据差异的对比．大地测量与地球动力学，34（4）：175—178．

方宏芳，刘序俨，张凯，2010．漳州地震台伸缩仪和体应变仪同震响应分析．地震研究，33（3）：287—291．

冯德益，聂永安，郭瑞芝等，1993．构造波及其在地震预报中的应用．地壳形变与地震，13（1）：1—9．

冯德益，吴国有，陈化然等，1994．地震波动力学特征变化指标在短期地震预报中的应用．地震，（1）：12—22．

公茂盛，谢礼立，章文波，2003．地震动输入能量衰减规律的研究．地震工程与工程振动，23（3）：15—24．

李宏男，孙立晔，2001．地震面波产生的地震动转动分量研究．地震工程与工程振动，21（1）：15—23．

刘琦，张晶，晏锐等，2013．高采样率四分量钻孔应变同震响应分析．中国地震，29（1）：57—67．

吕品姬，陈志遥，赵斌等，2010．定点倾斜观测映震能力综述．大地测量与地球动力学，30（S2）：50—56．

牛安福，张晶，吉平，2005．强地震引起的同震形变响应．内陆地震，19（1）：1—7．

牛安福，吉平，高福旺等，2006．印尼强地震引起的同震形变波．地震，26（1）：131—137．

邱永平，2011．宁波台动水位与垂直摆对印尼8.9和日本9.0级地震的同震响应．国际地震动态，（4）：28—33．

吴利军，冯琼松，张波，2016．体应变地震波最大振幅与地震震级和震中距的统计关系探讨．震灾防御技术，11（3）：592—599．

吴艳霞，2005．VS新型垂直摆倾斜仪的研究．武汉：武汉理工大学，1—58．

熊先保，林立峰，杨婕，2013．水管倾斜仪与摆式倾斜仪抗干扰能力分析．华南地震，33（2）：34—40．

杨玲英，崔庆谷，毛先进等，2012．云龙地震台水管倾斜仪与水平摆倾斜仪同震响应对比分析．地震研究，35（1）：48—52．

杨跃文，刘冰冰，秦嘉政等，2010．丽江台水管仪频谱分析及同震震后效应研究．防灾减灾学报，26（1）：61—65．

张利兵，刘仕锦，阳光等，2013．姑咱台和乡城台VS型垂直摆倾斜仪对比分析．四川地震，（4）：25—31．

张致伟，程万正，阮祥，2008．用SSY伸缩仪资料研究巨大远震引起的振荡波及谱．大地测量与地球动力学，28（5）：27—33．

Uang C. M., Bertero V. V., 1990. Evaluation of seismic energy in structures. Earthquake Engineering & Structural Dynamics, 19(1): 77—90.

Study on Attenuation Law of Co-seismic Deformation Waves Based on Instruments

Gong Liwen1), Chen Lijuan1), Liu Qi2)and Zheng Xudong1)

1) Chongqing Earthquake Agency, Chongqing 401147, China 2) Institute of Earthquake Forecasting, China Earthquake Administration, Beijing 100036, China

For 442 co-seismic deformation waves recorded by two sets of inclinometers in Qianjiang station, we conducted comparison and analysis on the spectrum features, and the statistics relationship between amplitudeand magnitudeThe results are as follows: ① Vertical pendulous, with smaller vibration period and damping coefficient, which can record abundant spectral information, is suitable for the time-frequency analysis for spectrum. Water-tube tilt, with longer baseline, bigger damping coefficient, and higher anti-interference capacity, which amplitude record is relatively stable, is suitable for analysis of energy attenuation. ② The amplitude- magnitude relationship is exponential, and the fitting degree of-curve increases with the epicenter distance. Whereas the index coefficientis inversely proportional with the ln. ③ The coefficient () is the main factor used to quantify the energy attenuation of co-seismic deformation waves, and the earthquakes in different areas are with different attenuation coefficients that increase with the epicenter distance (≤6000km). The above results are of significance for constructing model of the earthquake induced by huge earthquake, as well as for determining its parameter settings.

Water-tube tilt; Vertical pendulous; Co-seismic deformation waves; Attenuation coefficient; Amplitude-magnitude relationship

龚丽文，陈丽娟，刘琦，郑许东，2018．基于形变仪器同震形变波衰减规律的研究．震灾防御技术，13（2）：399—409．

10.11899/ zzfy20180215

2017年度震情跟踪定向工作任务（2017010231、162202）

2017-08-10

龚丽文，男，生于1988年。助理工程师。主要研究方向：定点形变监测及前兆异常。E-mail：glw777@126.com