引言

地震是地下介质受到应力作用产生破裂的自然现象，研究中强地震震源区应力状态的时空分布，对于了解地震的孕育过程具有重要意义． 陈颙(1978)通过对一些强震震例的分析研究发现，强震前震源机制趋于一致的现象时有发生，提出可用前震震源机制的一致性来描述地震活动性，同时给出了判定地震危险性的实例． 后续的相关研究也较多. 例如： 王俊国和刁桂苓(2005)的研究表明千岛岛弧M_W≥7.5地震前曾出现震源机制解与构造应力场的一致性参数降低的现象； 泽仁志玛等(2010)对2006年千岛岛弧M_W8.3地震和2001年秘鲁M_W8.4地震的研究表明，大震发生前孕震区震源一致性参数比较低，震源机制解与构造应力场的差异较小，说明前震受到了震源区应力场的统一作用，而主震发生后震源一致性参数比较高，说明后续地震的震源机制散乱，表明背景应力场的控制作用开始减弱．

2014年8月以来安徽省金寨县关庙乡发生M_L3.9震群序列，该震群序列持续时间长，震区原地1970年以来未记录到小震活动，说明该震群活动的位置不是老震区．

本文将利用P波、 SV波和SH波的极性和振幅比联合求解震源机制的方法(Snoke et al，1984； Snoke，1989)反演金寨M_L3.9震群序列中较大地震的震源机制解，并根据震源机制解的一致性参数(陈颙，1978)及震源机制解P轴方位角的变化分析该震群序列的震源一致性，同时利用震区附近3个地震台站的极性一致性和体波谱振幅相关系数(Lund，Boovarsson，2002)综合考察该震群序列的震源一致性相似程度．

1.   地震地质背景及震群序列基本情况

1.1   地震地质背景

2014年金寨M_L3.9震群序列发生在安徽省六安市金寨县关庙乡. 该地区位于大别山北麓，断裂构造复杂(姚大全等，2006)，展布有WNW向的磨子潭—晓天断裂、 梅山—龙河口断裂和NE向的落儿岭—土地岭断裂、 团风—温泉断裂(亦称商城—麻城断裂)． 金寨震群位于磨子潭—晓天断裂南侧，距其约2 km，震区附近固定地震台站相对密集. 鉴于该震群序列地震频发，安徽省地震局先后在震区附近架设了3个流动地震台，分别距震中2，5和13 km. 图 1给出了研究区内断裂和台站的分布．

图  1  金寨M_L3.9地震震中位置及震区断裂和台站分布

Figure  1.  Epicentral location of Jinzhai M_L3.9 earth-quake and distribution of faults and stations Solid triangles represent fixed seismic stations， and open triangles represent temporary ones

下载: 全尺寸图片 幻灯片

1.2   金寨震群序列基本情况

本文收集整理了2014年8月1日—11月30日安徽省地震台网记录到的金寨地震资料，包括安徽省的固定台站和所架设的3个流动台站，以及河南省、 湖北省紧邻震中区台站所记录到的地震波形数据． 在此基础上，剔除震相不清晰的资料，选取不少于3个台站记录到的M_L≥1.5地震，共计115次地震，3046条震相数据．

在对所选的115次地震进行双差重定位(Waldhauser，Ellsworth，2000)的过程中，选取震相观测记录数不少于6且震源间距小于10 km的地震事件进行组对，利用共轭梯度法反复迭代后得到了113次地震的精确震源位置，其余两次地震由于重定位组对过程中某些地震不满足条件而被剔除．

图 2给出了采用双差定位法得到的金寨震群序列地震精定位结果． 双差精定位结果显示，该序列地震的震源深度主要集中在5—6 km，占所定位地震总数的88%，表明该震群序列的震源深度较浅，分布也相对集中． 从地震的空间分布也可以看出，该序列的地震主要分布在5 km范围内，震中位置非常集中. 由于该序列的地震在空间上未出现明显的线性分布，因此很难确定其与周边断裂的关系． 从时间上可以看出，该震群序列分为两个活跃时间段： 第一个活跃时间段为8月底到9月底，序列地震活动强度相对较低； 第二个活跃时段为10月下旬到11月底，序列地震活动强度和频次均较高，震群序列活动有增强的趋势．

图  2  用双差定位法得到的金寨震群地震精定位结果

Figure  2.  Relocation result of the Jinzhai earthquake swarm sequence by using double difference location algorithm

下载: 全尺寸图片 幻灯片

2.   震源机制解反演及震源一致性参数

Snoke等提出了利用P波、 SV波和SH波的极性和振幅比联合求解震源机制的方法，简称为Snoke方法(Snoke et al，1984; Snoke，1989). 该方法利用P波、 SV波和SH波的极性，以及SV/P和SH/P振幅比等5个约束信息进行震源机制反演. 与P波极性法和垂直向SV/P振幅比方法相比，该方法约束量大为增加，反演结果更为可靠． 故本文采用Snoke方法来反演金寨震群序列中地震的震源机制解．

倪红玉等(2011)研究表明地壳速度结构模型对震源机制解的反演精度存在一定影响.经过对比，本文反演采用了表 1所示的速度结构模型，其中表层波速均为金寨、 六安和合肥地震台表层岩石实验室测试结果的平均值，其它各层速度均黄耘等(2011)给出的郯庐断裂带下扬子地块的P波速度，S波速度由v_S=v_P/1.732计算得到．

表  1  速度结构模型

Table  1.  Model of velocity structure

深度/km	层厚/km	vP/(km·s-1)	vS/(km·s-1)
0—2	2	4.98	2.94
2—5	3	5.22	3.01
5—10	5	5.59	3.23
10—15	5	6.12	3.53
15—20	5	6.46	3.73
20—25	5	6.45	3.72
25—33	8	7.03	4.06
33—40	7	8.00	4.62

下载: 导出CSV 

| 显示表格

鉴于篇幅所限，表 2仅列出了金寨震群序列中12次M_L≥3.0地震的震源机制解详细结果． 可以看出，除11月11日7时17分金寨M_L3.1地震的震源机制解为张性走滑型外，其余M_L≥3.0地震均为压性走滑型，表明该震群序列的震源机制解具有较好的一致性． 从表 2中还可以看到，该震群序列的12次M_L≥3.0地震的震源机制解P轴方位角集中在251°—278°，其中有7次地震的P轴方位角集中在256°—260°，这说明从P轴方位角也可以看出该震群序列的震源机制一致程度很好．

表  2  金寨震群M_L≥3.0地震的震源机制解

Table  2.  Focal mechanism solutions of M_L≥3.0 earthquakes of the Jinzhai earthquake swarm

发震时间 年-月-日	时：分	北纬/°	东经/°	ML	节面Ⅰ				节面Ⅱ				P轴			T轴
					方位角/°	倾角/°	滑动角/°		方位角/°	倾角/°	滑动角/°		方位角/°	倾角/°		方位角/°	倾角/°
2014-10-22	04:29	31.52	115.52	3.0	127	83	29		33	61	172		257	14		354	26
2014-10-22	13:34	31.52	115.52	3.7	125	83	19		33	71	173		258	8		350	18
2014-10-26	01:25	31.52	115.52	3.9	129	69	22		31	69	158		260	0		350	30
2014-10-26	02:16	31.52	115.51	3.7	127	83	29		33	61	172		257	14		354	26
2014-10-26	23:05	31.52	115.51	3.3	121	85	30		28	61	174		251	17		349	24
2014-10-28	19:33	31.52	115.52	3.6	138	82	24		44	67	171		269	10		3	23
2014-10-29	09:38	31.51	115.51	3.0	128	76	32		30	59	164		256	11		353	33
2014-11-06	06:22	31.51	115.53	3.6	126	82	24		32	67	171		257	10		351	23
2014-11-06	06:31	31.52	115.52	3.3	125	83	19		33	71	173		258	8		350	18
2014-11-11	03:08	31.51	115.52	3.7	141	85	30		48	61	174		271	17		9	24
2014-11-11	07:08	31.51	115.51	3.4	140	78	22		46	69	167		272	6		4	24
2014-11-11	07:17	31.52	115.51	3.0	323	81	-5		53	85	-171		278	10		188	3

下载: 导出CSV 

| 显示表格

同理，本文采用Snoke方法计算了金寨震群序列中26次M_L≥2.3地震的震源机制解，其中P轴方位角的时间变化进程如图 3所示． 可以看出，该震群序列8月底的地震P轴方位角约为278°，10月26日M_L3.9地震前后P轴方位角约为258°，11月以后的地震P轴方位角又转至278°附近． 刘泽民等(2011)利用震源机制解资料反演了东大别地区的构造应力场，认为该地区的最大主压应力轴方位角为267°，整个震群序列的震源机制解P轴方位角变化幅度约为25°，与震区构造应力场最大主压应力轴方位角的夹角最大约为10°． 本文结果与刘泽民等(2011)的结果吻合得很好． 故从P轴方位角的变化幅度及其与构造应力场的夹角可以看出，金寨序列地震震源机制解具有较好的一致性．

图  3  金寨震群序列的M-t及震源机制解P轴方位角的时序分布

Figure  3.  M-t plots and temporal distribution of the P-axis azimuth of the focal mechanism for the Jinzhai earthquake swarm sequence

下载: 全尺寸图片 幻灯片

本文使用应力场的最大和最小主压应力轴与震源机制解的P轴、 T轴夹角之和的1/2作为震源一致性参数，因此震源一致性参数最小值为0°，最大值为90°． 刘泽民等(2011)基于东大别地区的现代地震震源机制解资料，利用应力张量平均法计算得到了该地区的构造应力场特征，认为最大主压应力轴方位角为267°，倾角为5°，最小主压应力轴方位角为358°，倾角为4°． 本文将刘泽民等(2011)的东大别地区的构造应力场方向作为金寨震群附近的构造应力场方向，利用金寨震群序列26次M_L≥2.3地震的震源机制解，计算得到了该震群序列的震源一致性参数随时间的变化曲线，如图 4所示． 可以看出，该震群序列的震源一致性参数介于7°—17°，而震源一致性参数理论变化范围为0°—90°，该值越小说明震源一致性程度越高. 金寨震群序列的震源一致性参数一直处于较低值，说明其震源一致性非常好．

图  4  金寨震群序列震源一致性参数的时序分布

Figure  4.  Temporal distribution of the consistency parameters of focal mechanism for the Jinzhai earthquake swarm sequence

下载: 全尺寸图片 幻灯片

3.   体波谱振幅相关系数

3.1   原理

对于同一震源区的两次地震(或者两次地震之间的距离远小于地震与台站之间的距离)，Lund和Boovarsson(2002)提出可以由同一台站的地震记录反演两次地震直达P波、 S波的零频谱振幅值，通过计算其相关系数来描述其震源机制是否相似． 由于体波谱振幅相关系数给出了两次震中位置相近的地震事件在不同方位的零频谱值分量分布的相似性量度，因此可以用其相关系数来描述两次地震震源机制解的几何特性相似程度． 若体波谱振幅相关系数接近1，则各台记录的两次地震的波形应相似，震源机制解的几何特性也应相近； 若相关系数比较低，则两次地震的震源机制解的几何特性具有不同程度的差异． 因此，体波谱振幅相关系数可认为是其震源机制异同的归一化数值描述．

3.2   金寨震群序列谱振幅相关系数的时间变化

本文采用金寨震群震中附近波形记录较好的商城台(震中距31 km，方位角344°)、 金寨台(震中距40 km，方位角61°)和麻城台(震中距52 km，方位角218°)的波形资料，计算了该震群中8月1日—11月30日115次M_L≥1.5地震的多台谱振幅相关系数，计算过程中删除了由于发震时刻较近波形存在相互叠加的地震，具体计算步骤与程序详见刘泽民等(2010)文章． 图 5a给出了金寨震群中M_L≥1.5地震谱振幅相关系数的滑动平均时序曲线，其累积窗长为5次地震，滑动步长为1次地震. 可以看出： 震群的第一组活跃时段为8月22日—9月28日，最大地震为M_L2.6，其中8月22日—9月1日较为密集，该阶段谱振幅相关系数较高，在0.91—0.97之间波动，其平均值为0.95，9月2日—9月28日频度开始降低，谱振幅相关系数也降低，在0.71—0.92之间波动，其平均值为0.81； 第二组活跃时段从10月20日开始，最大地震为M_L3.9，震群活动水平明显上升，谱振幅相关系数在0.80—0.97之间波动，其平均值为0.90，整体较高. 图 5b给出了其中57次M_L≥2.0地震的谱振幅相关系数的滑动平均时序曲线.可以看出，谱振幅相关系数在0.86—0.97之间波动，其平均值为0.93，持续处于较高水平． 该震群序列中M_L≥1.5和M_L≥2.0地震的体波谱振幅相关系数的时序曲线，随着地震序列的谱相关系数呈起伏变化，但其整体仍处于较高水平(≥0.9)，由此可以推断该震群的震源机制相似性较好．

图  5  金寨震群的M-t及体波谱振幅相关系数时序分布 (a)M_L≥1.5；(b)M_L≥2.0

Figure  5.  M-t plots and temporal distribution of the correlation coefficient of spectral amplitude for the Jinzhai earthquake swarm sequence   (a)M_L≥1.5；(b)M_L≥2.0

下载: 全尺寸图片 幻灯片

4.   台站极性分布

震群序列台站记录到的(地震动)极性也能反映震源一致性的相似程度，若多个台站记录到的地震极性一致，则可认为该序列地震的震源机制一致，否则较离散． 本文收集了商城、 金寨、 麻城台记录的金寨震群中的115次M_L≥1.5地震的P波垂直向极性，给出了上述3个台站所记录到的地震极性时间曲线，如图 6所示. 可以看出: 商城台记录到的地震极性全部为正，金寨台除3次地震为正外，其余均为负，商城台与金寨台所记录到的地震极性表现为高度一致； 而麻城台所记录的极性则较为混乱． 麻城台的方位角为218°，由表 2所列的金寨序列较大地震的震源机制解详细结果可以看出，震源机制有一个节面的方位角集中分布在30°—50°，而麻城台正好位于该节面附近，因此导致了麻城台所记录到的地震极性较为混乱． 根据商城台与金寨台极性高度一致的现象，本文认为金寨震群整体极性分布优势明显，从时间上看，截至11月底，极性分布仍没有出现紊乱． 由于极性的一致性可以间接地反映震源机制的一致性，故本文认为该震群序列的震源一致性较好．

图  6  商城、 金寨和麻城台所记录到的金寨震群地震极性时间变化图   +1表示Pg波垂直向极性为正，0表示极性不清，-1表示Pg波垂直向极性为负

Figure  6.  Variations of seismic initial polarity recorded by the stations Shangcheng， Jinzhai and Macheng from the top down  +1 represents upcast polarity of vertical Pg wave，-1 represents downcast polarity， and 0 represents indistinct

下载: 全尺寸图片 幻灯片

5.   讨论与结论

2014年安徽省金寨县关庙乡发生显著震群序列，该震群自8月开始活动，截止到11月底，共记录到地震502次，其中M_L≥3.6地震6次，最大为M_L3.9． 双差定位研究显示，该震群序列的震源深度较浅，88%的序列地震的震源深度集中在5—6 km. 从空间分布看，该震群序列位于磨子潭—晓天断裂南侧2 km附近，序列地震集中在约5 km的范围内，空间分布方向特征不明显. 倪红玉等(2015)通过地震精定位也认为，金寨震群的空间分布无明显的展布方向，地震密集分布在2 km×2 km范围内，震源深度较浅，在3—7 km深度区间相对集中． 从时间上看，该震群序列主要划分为两个活跃时间段： 第一个活跃时间段为8月底至9月底，序列地震活动强度相对较低； 第二个活跃时段为10月下旬至11月底，序列地震活动强度、 频次均较高． 鉴于该震群序列震级有走强的趋势，本文采用多种方法综合研判了该震群序列的震源一致性．

金寨序列地震的震源机制解精确反演显示，绝大多数地震统一表现为压性走滑型，P轴方位角为251°—278°，波动幅度较小，P轴方位角与震区构造应力场最大主压应力轴方位角的最大夹角约为10°，这说明该序列地震震源机制解的一致性非常好． 本文定量计算了该序列地震的震源一致性参数，结果显示其变化范围为7°—17°，而理论变化范围为0°—90°，从震源一致性参数可以定量说明该震群序列的震源一致性非常好.倪红玉等(2015)采用Snoke方法和CAP(cut and paste)方法计算了金寨震群较大地震的震源机制解，结果表明M_L≥3.0地震震源机制解的一致性较好，节面位置基本一致，两组节面的走向分别呈NW向和NNE向，高倾角. 本文结果与倪红玉等(2015)的研究结果非常吻合．

金寨序列地震的体波谱振幅相关系数计算结果表明，虽然随着该震群序列的发展，谱振幅相关系数呈起伏波动，但整体仍处于较高水平(≥0.9)，同时震区附近台站记录到的地震极性的优势分布明显． 由此可推断该震群的震源一致性程度非常高．

金寨震区附近的东大别地区有4座大型水库，分别为梅山水库、 响洪甸水库、 佛子岭水库和磨子潭水库，均始建于20世纪五六十年代，其蓄水时间与本次震群活动时间相差甚远，加上4座水库与本次震群在空间上也存在一定距离，且震群活动前后4座水库的水位变化幅度与往年基本相同，因此本文认为本次震群活动与上述4座大型水库的蓄水无关．

此外，金寨震群附近3 km内还有两座小型电站，分别为2007年建成的黄谷滩电站(库容约100万m³)和2012年建成的关庙电站(库容约400万m³)，两座小型电站主要用于民用发电，无水位观测. 据震区政府介绍，金寨震群活动前后两电站处于正常发电水平，水位变化幅度相对较小． 虽然这两座小型电站库容非常小，蓄水时间也早于震群活动时间，但由于其距震群位置非常近，因此震群活动是否受该电站蓄水影响尚不确定．

崔子健等(2012)通过计算5个云南西北部发生的小震群序列，认为强震前小震群序列的谱振幅相关系数明显增大，接近于1.0，震源机制相似性显著； 而非前震序列谱振幅相关系数则明显较低，震源机制相似程度较弱． 黄浩和付虹(2014)使用相同的方法计算了滇西地区 2008—2011年的6个M_S≥6.0地震序列的谱振幅相关系数，认为振幅谱相关系数较高的地震序列也有可能是主余型序列，即谱振幅相关系数高的地震序列不一定都是前震系列． 因此，较低的谱振幅相关系数可以作为判别后续无更大地震的指标，但较高的相关系数可能仅仅反映了较高的应力水平，无法作为直接前震的判据． 因此，虽然金寨震群序列的体波谱振幅相关系数较高，但也只能说明其震源一致性相似，仅反映震区应力水平较高，而不能作为直接前震的判据． 同时本文认为金寨震群的震源一致性较好，也可能只反映震源区处于较高的应力状态，而在较高的应力水平下，是以发生中强地震还是以发生多次中等地震的方式释放其应变能，尚需进一步深入研究．