2013, Vol. 56
, 程鑫1 , 韩立波3 , 陈安国4 , 倪红玉1,2 2. 蒙城地球物理国家野外科学观测研究站, 合肥 230026;
3. 中国地震局地球物理研究所, 北京 100081;
4. 合肥工业大学资源与环境工程学院, 合肥 230009
2. National Geophysical Observatory at Mengcheng, Hefei 230026, China;
3. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China;
4. School of Resources and Environmental Engineering, Hefei University of Technology, Hefei 230009, China
研究表明地震前后介质应力的变化会导致岩石裂隙的张开或闭合，从而进一步影响地震波的散射特征.一些学者开始研究尾波Qc值随时间的变化特征.Chouet在1979年首次报道了尾波Qc随时间变化现象[1].随后陆陆续续的有学者报道了一些有关大震和火山喷发前尾波Qc值变化的现象，比如唐山地震前尾波Qc值变化了50%~70%[2]，Su和Aki也报道了震前尾波Qc值变化25%~29%的现象[3].其它一些学者同样发现了尾波Qc值变化的现象，但不是前兆性异常而是类似于同震响应的现象[4-8].还有一些学者发现地震前后尾波Qc值并没有出现明显的变化现象.关于尾波Qc值是否会出现时间上的变化一直存在争论，而争论主要集中在计算误差上.不支持尾波Qc值在震前发生变化的学者认为在研究尾波Qc值时由于震源机制解、震源深度和传播路径的不同使得计算结果缺乏可靠性.故而在用地震尾波Qc值研究地壳应力变化的时候如何有效减小非应力因素导致的Qc值计算结果的差异是关键问题[9].从而建议利用重复地震研究尾波Qc值进一步提高计算精度[10-11].例如，许多研究者用其研究介质变化特征[12-16].Beroza等[17]发现1989年Loma Prieta Mw6.9级地震前重复地震尾波Qc值出现5%的变化.Edoardo等[18]在研究意大利Mt. Vesuvius火山口附近地震序列时发现，在震群活动增强时，重复地震的尾波Qc值出现相应减小的现象.然而，Got等[9]和Got & Frechet[19]指出在频率域里两个重复地震虽然震中位置一致，但若破裂运动过程和辐射花样不同也会导致尾波Qc值的变化. Richard等[10]针对重复地震提出了滑动窗口谱振幅比方法，研究帕克菲尔德地区重复地震的尾波Q值变化情况，有效去除了各种干扰因素，提取尾波Qc值变化信息.本文将使用滑动窗口谱振幅比方法计算霍山窗重复地震序列尾波Qc值的变化特征，结合霍山窗小震活动分析引起尾波Qc值变化的因素.
2 滑动窗口谱振幅比方法本文使用的数据处理方法是由Richard提出的滑动窗口谱振幅比方法.由于国内还未有研究者使用该方法，本文简要介绍该方法的思路以便读者理解和参考.更详细的介绍可查阅相关参考文献[10-11].
Aki & Chouet[20]认为尾波振幅与频率有关，同时随时间衰减，可表示为：
|
(1) |
其中，S（f）表示震源谱，t是以发震时刻为起始的时间变化量，t的指数因子m表示几何衰减.计算一组重复地震尾波Qc值的变化可使用Got[9]和Beroza[17]提出的滑动窗口波谱比方法，对两个重复地震尾波部分进行如下处理：
|
(2) |
这里
|
(3) |
|
(4) |
可以看出Y（f，t）与时间和频率都有关，而W（f）只与频率有关，反映的是震源谱的变化对波形记录的影响.公式（2）表明使用这种处理方式可以将尾波Qc值的变化从包括震源谱不同在内的其它干扰源中提取出来.接下来需要处理的关键问题是如何建立方程求解ΔQc-1（f）.
滑动窗口谱振幅比方法的思路是将尾波分为M个有重叠的小窗口，每个窗口的长度为τ，窗口的中心点的时间为ti（i=1，2，…，M），具有离散频率fj（j=1，2，…，J/2-1，J是采用率）.计算第i个窗口内频率为fj时r（fj，ti）可由公式（5）获得：
|
(5) |
并通过公式（6）得到r（fj，ti）的计算误差：
|
(6) |
|
(7) |
|
(8) |
其中
|
(9) |
|
(10) |
公式（7）表示由多窗口谱分析方法[21]（Multitaper Spectral Analysis）获得的谱振幅比的自然对数，如公式（8）；公式（9）表示计算多窗口谱分析时不同本征窗口的权重；公式（10）表示在计算谱振幅比时的权重.有关多窗口谱分析的详细介绍请读者阅读相关文献[21].
有了公式（5），理论上可以建立由不同频率f和不同时间窗口t组成的关系式（2）的方程组，求解方程组获得尾波Qc值的变化.但事实上在实际观测中，表达式（6）中含有与时间无关的谱振幅比的成分，也就是需要扣除重复地震对波谱不同造成的影响.通过下面公式进行扣除，从而得到修正的波谱比：
|
(11) |
|
(12) |
公式（12）表示所有时间窗口的震源谱振幅比的加权均方差误差.另一方面，为了计算结果更稳定可靠，将频率分成L个互不重叠的频率段，中心频率为fl，每个频率段选取2V+1个频率点，建立由2V+1个方程组成的方程组，用最小二乘法求解该方程组即可求得该频率段ΔQ（fl）的值.其中V可由公式（13）获得：
|
(13) |
这里T是选取的时间窗口的长度.由2V+1个方程组成的方程组可表示如下：
|
(14) |
安徽省霍山窗位于安徽省中西部，处于地质构造复杂的北大别山地区.由于扬子克拉通与华北克拉通的陆陆碰撞，大别山地区尤其是霍山窗的地质构造比较复杂.图 1a展示霍山窗地区地形特征、断层构造特征、2008-2012年小地震分布和地震台站分布等信息.从图中可见霍山窗内有7条断裂穿过，其中磨子潭深断裂（也称青山-晓天断裂）是该区规模最大的断裂之一.该断裂走向NWW向，倾向NNE向.断裂南北分别为大别山群中深变质岩系和元古界浅变质岩系.地质学家根据大别山地质构造演化过程将大别山地区划分成多个地质构造单元[22]，而磨子潭深断裂正是两个主要构造单元（弧前复理石推覆体和反向冲断褶皱带与变质蛇绿混杂岩带）的交界线.同时也是华北地块与大别山微地块缝合带[23-25].
|
图 1 （a）霍山窗地区地质构造特征及地震、台站分布图；（b）研究区震级分布图；（c）研究区地震月频次图 Fig. 1 (a) The seismological and geological characteristics in Huoshan seismic window with small earthquake and station location; (b) Magnitude distribution of the small earthquakes; (c) Frequency distribution characteristic of the small earthquakes. |
目前已知的, 历史上安徽省共发生5级以上地震23次, 其中9次发生在霍山窗地区[26].近代小震活动也比较活跃, 且小震频次与周边地区中等强度地震有较好的对应关系[27], 故将霍山地区(115.7°E-116.5°E; 31°N-31.9°N)称为霍山震情窗, 简称霍山窗.图 1b和图 1c分别展示了霍山窗地区小地震的震级分布和月频次图.霍山窗地震的震级较小, 而时间分布极其不均匀, 有明显的活跃和平静现象.
4 数据选取及重复地震识别近年来有关重复地震的研究得到了国内许多研究者的关注[28-33].Rubin[34]定义“重复地震”为空间位置、波形和震级都具有较高的相似性的两次或多次地震事件.Schaff & Richards在研究中国境内重复地震时认为，两个地震在至少一个台站具有两个分量地震全波形相关系数大于0.8以上时，这两个地震是重复地震[35].目前国内大部分研究者都认为相关系数为0.8可以较好地识别重复地震.
考虑到霍山地区宽频带地震仪是2007年前后十五数字化改造后安装并投入使用的.我们整理了2008-2012年发生在安徽省霍山窗内的ML1.0级以上的震相数据.使用Hyp2000[36]方法进行精定位，定位结果如图 1.在地震密集的地方布设有4个台站，这4个台站的震中距较小，基本都在10 km以内.而滑动窗口谱振幅比方法需要2倍S波到时至少12 s以上，即震中距20 km以上.震中距小于这个距离的地震，尾波发育不完善影响计算可靠性.周边台站中满足条件的是六安台，其震中距为35 km左右.其它台站因震中距较大，不能较好记录霍山窗地区的小震.因而本文选取六安台作为研究对象.因为六安台震中距相对较远，为保证数据的信噪比，选取震级大于1.5（ML）的地震.去仪器响应，滤波到1~25 Hz，分别计算3个分量波形的相关系数.本文计算尾波Q值，并不需要严格意义上的重复地震，故选取三分量平均相关系数大于或等于0.75的地震对作为重复地震对.通过重复地震对的关系，搜索重复地震序列.最后只得到一组重复地震序列，其时间跨度从2008年到2012年，如图 2.
|
图 2 重复地震序列RES（Repeating Earthquake Series）的波形图.每个波形图上标注了地震的发生时间.红色虚线方框表示重复地震序列的尾波部分，起始时间为2倍S波走时. Fig. 2 The seismic wave of Repeating Earthquake Series with labeled event time. The dash frame denotes the seismic coda segment used in our research. |
地震尾波是由于地球介质的散射产生的，一般认为2倍S波走时后的波形记录即为尾波.因本文所使用方法并非计算单一地震的尾波Qc值，而是通过波形对比的方式获取重复地震对的尾波Qc值的变化.为了减小误差，更客观反映重复地震尾波Qc值的变化，在使用滑动窗口谱分析计算之前，需要将重复地震对的尾波进行对齐.对于重复地震对，一般认为将S波对齐即可认为尾波部分已对齐.本文选取S波前0.5 s和S波后3 s的时间段，并通过波形互相关方法将之对齐.以2012年12月15日的地震为参考地震，其它地震均与之进行对齐，选取该参考地震的2倍S波走时为尾波起始时间，而选取尾波振幅与背景噪音振幅相当时为尾波结束时间.根据该规则本文选取尾波长度20 s.得到各地震的尾波数据后，使用由MatLab语言编写的滑动窗口谱振幅比程序计算得到各地震对的尾波Qc-1随频率的变化曲线（图 3）.计算过程中，选取窗口长度1.5 s，滑动步长为0.25 s.图 3a展示的是每个地震与2012年地震的尾波ΔQc-1随频率变化情况，可以看出不同地震对的尾波ΔQc-1随频率的变化形态有明显区别.总体上，15 Hz以下的频段，各地震对尾波ΔQc-1有不同程度的变化，而15 Hz以上基本无变化.其中2009年的3次地震与2012年相比，在1~3 Hz和4~6 Hz频段有明显变化，10~12 Hz和13~15 Hz频段有较小变化，其它频段变化不明显.图 4展示了重复地震序列尾波Qc-1值的变化量（ΔQc-1）在各频段上随时间的变化情况.可以清楚地看出，1~3 Hz和4~6 Hz频段上，重复地震序列尾波ΔQc-1的变化经历了2008年无明显变化到2009年明显降低，2010年后恢复到基本正常的状态.其它频段ΔQc-1值基本稳定，没有明显变化.
|
图 3 每个地震相对于参考地震的尾波ΔQc-1值随频率的变化 Fig. 3 The differential variation of coda ΔQc-1 between one of the RES and the reference event |
|
图 4 重复地震序列尾波Q c-1值的变化量（ΔQc-1）在各频段上随时间的变化 Fig. 4 The time characteristics of variation of differential coda ΔQc-1 of the RES in Huoshan seismic window |
对于2009年霍山窗重复地震序列尾波ΔQc-1在低频段出现的变化现象，需要进一步分析原因.图 5绘出了霍山窗重复地震序列尾波ΔQc-1在1~3 Hz频段的变化曲线与霍山窗ML1.0级以上地震的频次的对比情况.图中黑色点线为尾波ΔQc-1变化曲线；黑色曲线为小震三月累计一周滑动频次曲线；灰色竖线为小震M-T图.可以清楚看出霍山窗小震频次在2009年之前处于较低水平，2009年出现明显升高，到2009年8月达到峰值，之后开始回落，2009年底基本恢复到正常水平.这期间霍山窗内发生了近几年最大的地震，为2009年8月30日ML 3.9级地震.但该次地震震级较小，本文作者认为，它的孕育过程不足以导致尾波Qc出现如此明显变化，而是同小震频次上升一样，可能是该地区应力场发生变化后出现的现象.
|
图 5 霍山窗重复地震序列尾波ΔQc-1在1~3 Hz频段范围的变化曲线（黑色点线）与霍山窗ML1.0级以上地震的三月累计频次（黑色曲线）的对比图，灰色竖条为ML1.0级以上地震的M-T图 Fig. 5 Comparison of the temporal variation of differential coda ΔQc-1 at frequency band 1~3Hz (black dot line) and the cumulative frequency of small earthquake with magnitude≥1.0 in three months (black curve), and the gray bar denote the M-T map |
本文用滑动窗口谱振幅比方法计算了霍山窗重复地震序列尾波ΔQc-1值的变化特征.结果发现2009年霍山窗小震活动较强时，重复地震序列的尾波Qc值发生了变化.由于霍山窗地震震级较小，在平均震中距35 km左右六安台只有震级ML≥1.5的地震有较好的信噪比.本文只找到一组重复地震序列参与计算，缺乏同时期大量数据的对比分析.然而我们使用滑动窗口谱振幅比方法计算重复地震对之间尾波Qc值的相对变化，可有效减小射线路径和震源谱的差别导致的计算误差，使计算结果相对可靠.可通过公式（15）（16）进一步推导出尾波Q值的变化量的百分比：
|
(15) |
|
(16) |
计算需要本地区尾波Qc值的绝对量.而尾波Q值在不同地区的研究结果表现出较大差异[37-40].郑先进等[41]研究安徽地区尾波Qc值发现不同台站计算出的尾波Qc值亦表现出明显差别.若以六安台的尾波Qc值研究结果（Q=120±40）作为本文计算的参考值，发现2009年霍山地区尾波Qc值增加了11%~21%.
地震尾波是岩石介质中的非均匀体的散射产生的.岩石裂隙也是非均匀体的重要部分，是导致地震波散射的重要原因.应力变化导致岩石裂隙特征如裂隙长度、宽度、流体静压力等的变化会给地震波散射造成直接影响.大量的研究结果表明震前尾波Qc值的变化不管是在量值上还是在形态上都不具有相似性.有些研究发现震前尾波Qc值减小[42-44]；而有些则发现震前尾波Qc值增加[17, 45-47]；但不管怎么样，尾波Qc值的研究表明其变化与介质应力变化是存在内在联系的.
本文研究结果发现霍山窗地区尾波Qc值在2009年出现增加现象.该现象出现后霍山窗及附近地区并没有中强地震，而是伴随着霍山窗小震活动增强现象.故我们认为霍山地区地壳应力的变化与前人研究的地震成因有关的应力变化有本质区别，更多的可能是外部应力作用的影响.
根据GPS观测到的西部测站运动速率VN和VE [48]可以认为，印度板块向N推挤的作用力在大陆内部传递的过程中，青藏亚板块的SN向缩短经过阿尔金断裂传递到塔里木和天山块体，但挤压率逐渐减小，最大压缩在拉萨和温泉之间.物质的E向滑移在喜马拉雅地区不明显，这里以挤压为主，而最大的E向滑移发生在冈底斯块体和羌塘块体.由此可见，由于西侧块体明显地向E侧向运动，对中国东部特别是大华北地区会产生影响.图 6展示的是GPS观测到的大华北及邻区水平运动速度场在2008年汶川8级特大地震前后的变化情况.可以看出汶川地震前大华北南部地区GPS水平速度场以SW方向为主，而汶川地震后变为SE方向.可见汶川8级特大地震对华北南部地区的应变有影响作用.通过统计大陆西部强震与大华北地区中强地震的对应关系也可以说明这种影响作用的存在[49].而霍山窗地质构造位置特殊，会受到应力变化的影响，进而出现小震活动增强现象[50].统计霍山窗开窗时间与西部强震发生时间的关系发现，1970年以来20次霍山窗开窗有14次在开窗前一年内，西部发生过7级以上地震.可见，大陆西部的构造运动和应力传递对霍山地区地壳应力造成一定影响，出现尾波Qc值变化，并引起霍山地区小震活动增强.
|
图 6 汶川8级地震前后大华北地区GPS水平速度场分布 （a）1997-2007；（b）2007-2009. Fig. 6 The GPS horizontal velocity field before and after the Wenchuan M8.0 earthquake |
本文通过波形互相关方法识别霍山窗重复地震对，进一步搜索重复地震序列.人工挑选一个时间跨度和采样较好的重复地震序列.用滑动窗口谱振幅比方法计算该重复地震尾波Qc值的变化，发现在2009年尾波Qc值出现明显增加，增幅可达11%~21%.分析认为霍山窗重复地震尾波Q值的变化与前人研究的震情Qc值变化有所不同，可能是受到大陆西部物质东移和应力向东传递的影响.
致谢感谢中国地震局地球物理研究所蒋长胜博士、中国地震局地震预测研究所邵志刚博士对本文研究思路提出的宝贵意见；感谢中国地震局第一监测中心杨国华老师提供的2008年汶川地震前后大华北地区GPS水平速度场图片和对本文观点提供的宝贵意见.
| [1] | Chouet B. Temporal variation in the attenuation of earthquake coda near Stone Canyon, California. Geophys. Res. Lett. , 1979, 6(3): 143-146. DOI:10.1029/GL006i003p00143 |
| [2] | Jin A, Aki K. Temporal change in coda Q before the Tangshan earthquake of 1976 and the Haicheng earthquake of 1975. J. Geophys. Res. , 1986, 91(b1): 665-673. DOI:10.1029/JB091iB01p00665 |
| [3] | Su F, Aki K. Temporal and spatial variation in coda Q-1 associated with the North Palm Springs earthquake of July 8, 1986. Pure and Applied Geophysics , 1990, 133(1): 23-52. DOI:10.1007/BF00876701 |
| [4] | Gusev A, Lemzikov V. The anomalies of small earthquake coda wave characteristics before the three large earthquakes in the Kuril-Kamchatka zone. Vulk. Seism. , 1984, 4: 76-90. |
| [5] | Sato H. Temporal change in attenuation intensity before and after the eastern Yamanashi earthquake of 1983, in central Japan. J. Geophys. Res. , 1986, 91(b2): 2049-2061. DOI:10.1029/JB091iB02p02049 |
| [6] | Peng J Y, Aki K, Chouet B, et al. Temporal change in coda associated with the Round Valley, California, earthquake of November 23, 1984. J. Geophys. Res. , 1987, 92(B5): 3507-3536. DOI:10.1029/JB092iB05p03507 |
| [7] | Tsukuda T. Coda Q before and after the 1983 Misasa earthquake of M6.2, Tottori Prefecture. Japan, Pageoph. , 1988, 128: 261-280. DOI:10.1007/BF01772600 |
| [8] | Jin A, Aki K. Spatial and temporal correlation between coda Q-1 and seismicity and its physical mechanism. J. Geophys. Res. , 1989, 94(B10): 14041-14059. DOI:10.1029/JB094iB10p14041 |
| [9] | Got J L, Poupinet G, Fréchet J. Changes in source and site effects compared to coda Q-1 temporal variations using microearthquakes doublets in California. Pure and Applied Geophysics , 1990, 134(2): 195-228. DOI:10.1007/BF00876998 |
| [10] | Aster R C, Slad G, Henton J, et al. Differential analysis of coda Q using similar microearthquakes in seismic gaps. Part 1: Techniques and application to seismograms recorded in the Anza seismic gap. Bull. Seis. Soc. Amer. , 1996, 86(3): 868-889. |
| [11] | Antolik M, Nadeau R M, Aster R C. Differential analysis of coda Q using similar microearthquakes in seismic gaps. Part 2: Application to seismograms recorded by the Parkfield high resolution seismic network. Bull. Seis. Soc. Amer. , 1996, 86(3): 890-910. |
| [12] | Niu F L, Silver P G, Daley T M. Preseismic velocity changes observed from active source monitoring at the Parkfield SAFOD drill site. Nature , 2008, 454(7201): 204-209. DOI:10.1038/nature07111 |
| [13] | Niu F L, Silver P G, Nadeau R M, et al. Migration of seismic scatterers associated with the 1993 Parkfield aseismic transient event. Nature, 2002, 1993, 426(6966): 544-548. |
| [14] | Hui L, Wen L X. Using repeated sources to quantitatively determine temporal change of medium properties: theory and an example. J. Geophys. Res. , 2012, 117(B9). DOI:10.1029/2012JB009302 |
| [15] | Sumiejski L E, Shearer P M. Temporal stability of coda Q-1 in sourthern California. Bull. Seis. Soc. Amer. , 1995, 102(2): 873-877. |
| [16] | 周龙泉, 刘桂平, 马宏生, 等. 利用重复地震观测地壳介质变化. 地震 , 2007, 27(3): 1–9. Zhou L Q, Liu G P, Ma H S, et al. Monitoring crustal media variation by using repeating earthquakes. Earthquake (in Chinese) , 2007, 27(3): 1-9. |
| [17] | Beroza G C, Cole A T, Ellsworth W L. Stability of coda wave attenuation during the Loma Prieta, California earthquake sequence. J. Geophys. Res. , 1995, 100(B3): 3977-3987. DOI:10.1029/94JB02574 |
| [18] | Edoardo D P, Francesca B, Simona P, et al. Changes in the coda decay rate and shear wave splitting parameters associated with seismic swarms at Mtt. Vesuvius, Italy. Bull. Seis. Soc. Amer. , 2004, 94(2): 439-452. DOI:10.1785/0120030141 |
| [19] | Got J L, Fréchet J. Origins of amplitude variations in seismic doublets: source or attenuation process?. Geophys. J. Int. , 1993, 114(2): 325-340. DOI:10.1111/gji.1993.114.issue-2 |
| [20] | Aki K, Chouet B. Origin of coda waves: source, attenuation, and scattering effects. J. Geophys. Res. , 1975, 80(23): 3322-3342. DOI:10.1029/JB080i023p03322 |
| [21] | Park J, Vernon F L, Lindberg C R. Multitaper spectral analysis of high-frequency seismograms. J. Geophys. Res. , 1987, 92(B12): 12657-12684. |
| [22] | 徐树桐, 江来利, 刘贻灿, 等. 大别山区(安徽部分)的构造格局和演化过程. 地质学报 , 1992, 66(1): 1–15. Xu S T, Jiang L L, Liu Y S, et al. Tectonic framework and evolution of the Dabie mountains in Anhui, Eastern China. Acta Geologica Sinica (in Chinese) , 1992, 66(1): 1-15. |
| [23] | Wang Q, Liu X, Maruyama S. Top boundary of the Dabie UHP rock, Central China. GSA Abstract with Programs , 1993, 25: 99. |
| [24] | Wang Q, Liu X, Maruyama S, et al. Top boundary of the Dabie UHP rock, Central China. J. SE. Asia Geol. , 1995, 11(4): 295-300. |
| [25] | 王清晨, 从柏林, 马力. 大别山造山带与合肥盆地的构造耦合. 科学通报 , 1997, 42(6): 575–585. |
| [26] | 陆镜元, 曹光暄, 刘庆忠, 等. 安徽省地震构造与环境分析. 合肥: 安徽科学技术出版社, 1992 : 82 -93. |
| [27] | 郑兆苾, 庆梅, 李敏莉. 霍山窗与华东中强地震关系研究. 地震学刊 , 1999, 99(2): 1–9. Zheng Z B, Qing M, Li M L. Study on correlation between Huoshan seismic window and occurrence of moderate earthquake of East China. Journal of Seismology (in Chinese) , 1999, 99(2): 1-9. |
| [28] | 李乐, 陈棋福, 钮凤林, 等. 利用"重复地震"估算丽江-宁蒗断裂带的深部滑动速率. 科学通报 , 2008, 53(23): 2925–2932. |
| [29] | 李乐, 陈棋福, 钮凤林. 唐山断裂带的板内重复地震识别及其滑动速率估算. 国际地震动态 , 2008, 11: 2. |
| [30] | 李宇彤, 吴忠良, 蒋长胜, 等. 利用辽宁区域地震台网记录分析"重复地震". 地震学报 , 2008, 30(4): 383–396. Li Y T, Wu Z L, Jiang C S, et al. Analysis on repeating earthquakes recorded by Liaoning regional seismograph network. Acta Seismologica Sinica (in Chinese) , 2008, 30(4): 383-396. |
| [31] | 蒋长胜, 吴忠良. 由"重复地震"给出的中国地震台网的定位精度估计. 中国地震 , 2005, 21(2): 147–154. Jiang C S, Wu Z L. Estimating the location accuracy of the China National Seismograph Network using repeating events. Earthquake Research in China (in Chinese) , 2005, 21(2): 147-154. |
| [32] | 蒋长胜, 吴忠良, 李宇彤. 首都圈地区"重复地震"及其在区域地震台网定位精度评价中的应用. 地球物理学报 , 2008, 51(3): 817–827. Jiang C S, Wu Z L, Li Y T. Estimating the location accuracy of the Beijing capital digital seismograph network using repeating events. Chinese J. Geophys (in Chinese) , 2008, 51(3): 817-827. |
| [33] | 郑先进, 刘泽民, 洪德全, 等. 重复地震及其在安徽霍山地区的初步应用. 国际地震动态 , 2009(4): 88. Zheng X J, Liu Z M, Hong D Q, et al. Doublet and its preliminary application in Huoshan area in Anhui. Recent Developments in World Seismology (in Chinese) , 2009(4): 88. |
| [34] | Rubin A M. Using repeating earthquakes to correct high-precision earthquake catalogs for time-dependent station delays. Bull. Seism. Soc. Amer. , 2002, 92(5): 1647-1659. DOI:10.1785/0120010180 |
| [35] | Schaff D P, Richards P G. Repeating seismic events in China. Science , 2004, 303(5661): 1176-1178. DOI:10.1126/science.1093422 |
| [36] | Fred W K. User's guide to Hypoinverse-2000, a Fortran program to solve for earthquake location and magnitudes. U. S. Geological Survey Open-File Report. 2002. |
| [37] | 陈俊华, 甘家思. 九江-瑞昌5.7级地震余震尾波Qc值分析. 大地测量与地球动力学 , 2006, 26(3): 65–68. Chen J H, Gan J S. Analysis of Qc value of coda waves about after shocks of Jiujiang-Ruichang Ms55.7 earthquake. Journal of Geodesy and Geodynamics (in Chinese) , 2006, 26(3): 65-68. |
| [38] | 杨云, 霍祝青. 苏中沿海地区尾波Qc值特征研究. 西北地震学报 , 2011, 33(2): 166–170. Yang Y, Huo Z Q. Characteristics of coda wave Qc value in the coastal area of central Jiangsu province. Northwestern Seismological Journal (in Chinese) , 2011, 33(2): 166-170. |
| [39] | 赵卫明, 杨明芝, 金延龙, 等. 宁夏中北部区域尾波Q值特性研究. 地震学报 , 2000, 22(3): 249–256. Zhao W M, Yang M Z, Jin Y L, et al. Characteristics of coda Q-value in the mid-northern part of Ningxia. Acta Seismologica Sinica (in Chinese) , 2000, 22(3): 249-256. |
| [40] | Hiramatsu Y, Hayashi N, Furumoto M, et al. Temporal changes in coda Q-1 and b value due to the static stress change associated with the 1995 Hyogo-ken Nanbu earthquake. J. Geophys. Res. , 2000, 105(B3): 6141-6151. DOI:10.1029/1999JB900432 |
| [41] | 郑先进, 刘东旺, 沈小七, 等. 安徽地区尾波Qc值特征研究. 华北地震科学 , 2011, 29(3): 8–12. Zheng X J, Liu D W, Shen X Q, et al. Research on the Qc value of Anhui area. North China Earthquake Sciences (in Chinese) , 2011, 29(3): 8-12. |
| [42] | 啜永清, 张淑亮, 苏燕, 等. 大同地震序列的尾波Q值变化特征研究. 地震 , 2004, 24(2): 66–72. Chuo Y Q, Zhang S L, Su Y, et al. Research on variation characteristics of seismic coda Q value of Datong earthquake sequence. Earthquake (in Chinese) , 2004, 24(2): 66-72. |
| [43] | 林畅彦, 古本宗充, 平松良浩, 等. 兵库县南部地震引起的尾波Q值的变化. 世界地震丛译 , 1998: 40–53. |
| [44] | 朱新运, 张帆. 浙江及邻区尾波值研究. 西北地震学报 , 2006, 28(2): 108–113. Zhu X Y, Zhang F. Research on coda Q value in Zhejiang province and neighboring area. Northwestern Seismological Journal (in Chinese) , 2006, 28(2): 108-113. |
| [45] | 黄耘, 孙业君, 王斌, 等. 2011年安徽安庆4.8级地震前尾波Q值变化研究. 国际地震动态 , 2011: 117. |
| [46] | 王伟君, 刘杰. 1999年岫岩地震序列尾波Qc的变化过程. 地震 , 2004, 24(4): 37–44. Wang W J, Liu J. Temporal process of S-Coda Qc of Xiuyan earthquake sequence in 1999. Earthquake (in Chinese) , 2004, 24(4): 37-44. |
| [47] | Sato H. Temporal changes in scattering and attenuation associated with the earthquake occurrence-a review of recent studies on coda waves. Pure App. Geophys. , 1988, 126(2-4): 465-497. DOI:10.1007/BF00879007 |
| [48] | 赖锡安, 徐菊生, 卓力格图, 等. 中国大陆主要构造块体现今运动的基本特征. 中国地震 , 2000, 16(3): 213–222. Lai X A, Xu J S, Zhuoli G T, et al. The fundamental characteristics of principal tectonic blocks present movement in Chinese mainland. Earthquake Research in China (in Chinese) , 2000, 16(3): 213-222. |
| [49] | 陈宇卫, 张军, 闫素萍, 等. 西部地区特大地震活动特征及其对大华北地震趋势的影响. 地震地质 , 2003, 25(4): 647–654. Chen Y W, Zhang J, Yan S P, et al. Characteristics of great earthquake activity in west China and its influence on seismic tendency in North China. Seismology and Geology (in Chinese) , 2003, 25(4): 647-654. |
| [50] | 陈宇卫, 张军, 庆梅, 等. "霍山地震窗"小震序列运动学参数时变特征研究. 地震 , 2007, 27(1): 26–32. Chen Y W, Zhang J, Qing M, et al. Study on the time-varying characteristics of kinematic parameters of small earthquake sequence in "Huoshan seismic window". Earthquake (in Chinese) , 2007, 27(1): 26-32. |