引言

在中国大陆东西部中间有一条狭长的过渡地带(99°E—106°E)，它不仅是地势、 重力和地壳厚度的梯度带，而且其地震活动水平之高也为世人瞩目，这就是人们常说的南北地震带(毛可等，1997)． 在某种程度上，南北地震带可作为青藏地块向东运动并挤压华北地块的代表地带． 当印度板块运动加剧引起我国大陆西部地块地震活动高潮后，东部华北地块随之也相应出现地震活动高潮． 20世纪，一系列大地震沿青藏地块北部边界构造带(南北地震带中北段)连续发生，包括1920年海原M_S8.5、 1927年古浪M_S8.0、 1933年叠溪M_S7.5、 1954年山丹M_S7.3和民乐M_S7.0、 1976年松潘M_S7.2和平武M_S7.2等地震，说明青藏地块这一条边界构造带在50年左右的时间内活动加强，从而集中释放出巨大的地震能量． 与此同时华北地块也相继发生了1937年菏泽M_S7.0、 1966年邢台M_S7.2、 1969年渤海M_S7.4、 1975年海城M_S7.3和1976年唐山M_S7.8以及一系列M_S≥6.0地震． 朱兆才(1994)研究发现，1966—1976年青藏地块与华北地块发生的M_S≥7.0地震有相互跳迁的呼应关系． 因此，继2008年以来南北地震带中北段相继发生了2008年5月12日汶川M_S8.0和2013年4月20日芦山M_S7.0地震后，接下来华北地块的地震又该如何发展呢？

1998年1月10日张北M_S6.2地震至今，华北地块M6.0地震平静已达17年多． 2011年1月12日南黄海M_S5.0地震打破了华北地块近5年半时间的M5.0地震平静格局，之后陆续发生了2013年1月23日辽宁灯塔M_S5.1、 2013年4月21日黄海M_S5.0、 2013年4月22日内蒙科尔沁M_S5.3和2013年5月18日黄海M_S5.1等一系列M_S≥5.0地震． 这种M5.0地震密集连发的现象为历史所罕见，且上述M_S≥5.0地震在空间上明显呈带状分布，使得华北地块M6.0地震的形势越发复杂而严峻． 因此，本文从动力边界和活动地块的角度，侧重于研究南北地震带中北段M_S≥7.0地震对华北地块M_S≥6.0地震的影响，以期对未来华北地块M_S≥6.0地震的发震地点和趋势判定工作提供借鉴和参考．

1.   动力背景分析

大量地震的震源机制解和地应力实测结果表明(邓起东等，1979； 丁旭初，张文涛，1988； 杨树新等，2012)，我国大陆存在近似于辐射状的区域应力场，其辐射中心位于青藏断块东部，自西向东，从近N-S方向逐步顺时针旋转至NNE- -SSW，NE- -SW，ENE- -WSW及NW- -SE方向． 在青藏高原地区地壳厚度为46—56 km，急剧变化梯度带向东凸出，凸出地带前方两侧分布有鄂尔多斯盆地和四川盆地. 前缘的龙门山一翼出现一系列ENE向挤压为主的右旋逆断层，而六盘山一翼则存在NW走向挤压为主的左旋逆断层，表明青藏地块在此对两翼的挤压和向东的滑动趋势． 杨树新等(2012)研究表明，南北地震带位于中国大陆东西部两个一级地应力场分区的结合部位，受到印度板块、 太平洋板块和菲律宾板块的联合作用．

新构造时期以来，华北地块的应力场是太平洋板块向欧亚大陆俯冲以及与印度板块碰撞作用下青藏地块向我国东部大陆挤压的联合作用下形成的(邓起东等，2002)，华北地震的发生也是两大板块运动共同作用于华北地块的结果(朱兆才，1994)． 无论是单次大地震还是利用单台记录的小地震资料所得的华北地块地震震源机制解均较为一致，即发震应力场的主压应力方向为ENE-WSW，发震断层一般均为具有大倾角的走滑断层(吴佳翼等，1979)． 由此也说明，华北地块不仅受太平洋板块西向俯冲作用的影响，也是印度板块自西南向东北作用于我国大陆地壳的结果． 地质学和地球物理学的研究结果(嵇少丞等，2008)也认为，青藏地块在近SN向压缩、 抬升和地壳增厚的同时，也存在着EW向的延伸，并对华北和华南地块产生推挤作用．

华北地台是中国最大的古老克拉通，位于中亚造山带与秦岭-大别-苏鲁造山带之间，如图 1所示． 可以看出，华北地块的东边界为郯庐断裂带，北边界为西拉木伦断裂． 秦岭-大别-苏鲁、 松潘-甘孜、 祁连和昆仑等造山带的交汇部位及其附近区域，即南北地震带的中北段、 龙门山-六盘山一带是印度板块作用下青藏地块挤压华北地块的典型作用地带． 因此该地区强烈地震活动所产生的应力扰动和调整，必然会对中国东部特别是大华北地块的地震活动产生影响．

图  1  华北地块和秦岭-大别造山带区域地质简图(改绘自洪大卫等， 2003)

Figure  1.  The regional geological simplified map of North China block and Qinling-Dabie orogenic belt (revised after Hong et al， 2003)

下载: 全尺寸图片 幻灯片

2.   地震活动分析

考虑到地震资料的完整性，本文所用资料引自1900年以来中国地震台网中心的全国M5.0地震目录库. 据统计，1900年以来我国大陆东部地区(30°N—44°N、 104°E—125°E)共发生M_S≥6.0地震(不计余震)34次，其中南北地震带中北段M_S≥7.0地震9次(表 1)． 大量震源机制解压应力轴所标示的我国现代构造应力场压应力分布方向与利用GPS观测得到的我国现今水平地壳运动方向一致性表明(汪素云等，1996； 周硕愚等，2001)： 我国大陆和海域现代构造应力场与现地壳水平运动具有统一性和确定的方向性； 印度板块向北推挤青藏高原，以及高原东部被挤出的构造地块往NE- -SE方向的滑动挤压，是我国东部海域及其邻区现代构造应力场形成的主要动力(徐杰等，2012)． 因此，为了更直观地反映印度板块作用下青藏地块对华北地块的挤压作用，本文以我国东部辐射状应力场的近似中心(35°N、 100°E)为动力源点，如图 2中三角形所示，研究华北地块M_S≥6.0地震震中与该简化动力源点的距离及其与地震发生时间的关系. 该距离由下式确定(张少泉，1977):

式中： λ和φ为地震震中经度和纬度(λ₀=100°，φ₀=35°)； Δ为地震震中与简化动力源点的距离，单位为km.

图  2  中国大陆及其邻区震源机制主压应力轴和GPS运动模型速度矢量分布图(引自周硕愚等，2001)

Figure  2.  Distribution of main compression stress of focal mechanism solutions and GPS velocity vectors in Chinese mainland and its adjacent areas (after Zhou et al， 2001)

下载: 全尺寸图片 幻灯片

各地震震中到动力源点距离及相关信息见表 1．可以看出，南北地震带中北段发生M_S≥7.0地震后，华北地块后续随之发生一系列M_S≥6.0地震，且每组首发的M_S≥6.0地震距南北地震带中北段最近一次M_S≥7.0地震的时间间隔为1个月至11.8年，平均为3.5年． 每组M_S≥6.0地震的 持续活跃时间基本接近，依次为16.1，14.5，13.5和18.6年，期间华北地块发生M_S≥6.0地震的次数为4—7次．

表  1  1900年以来我国大陆东部地区M_S≥6.0地震及其到动力源点距离

Table  1.  M_S≥6.0 earthquakes in eastern part of Chinese mainland since the year 1900 and their distances to dynamic source point

组号	序号	发震时间	北纬/°	东经/°	MS	发震时间	Δ/km	Δt	MS≥6.0地震持续活跃时间
		年-月-日
		1910-01-08	35.0	122.0	6.7	山东青岛	2214
		1917-01-24	31.3	116.2	6.2	安徽霍山	1699
		1917-05-28	39.7	124.0	6.1	辽宁丹东	2436
	*1	1920-12-16	36.7	104.9	8.5	宁夏海原	526
	2	1922-09-29	39.2	120.5	6.5	渤海	2089
	3	1929-01-14	40.7	111.3	6.0	内蒙古呼和浩特	1284
①	4	1932-04-06	31.4	115.0	6.0	湖北麻城	1577	1.8年	16.1年
	5	1940-01-19	42.7	121.3	6.0	内蒙古奈曼旗	2259
	6	1942-07-09	43.5	121.9	6.0	内蒙古通辽	2342
	7	1944-12-19	39.7	124.3	6.8	鸭绿江口	2465
	*8	1927-05-23	37.7	102.2	8.0	甘肃古浪	372
	*9	1933-08-25	32.0	103.7	7.5	四川叠溪	502
	10	1934-01-21	41.1	108.3	6.2	内蒙古五原	1063
②	11	1937-08-01	35.3	115.2	7.0	山东菏泽	1529	5个月	14.5年
	12	1945-09-23	39.5	119.0	6.2	河北滦县	1950
	13	1948-05-23	37.7	121.8	6.0	山东威海	2196
	*14	1954-02-11	39.0	101.5	7.3	甘肃山丹	469
	*15	1954-07-31	38.8	104.2	7.0	甘肃民乐	594
	16	1966-03-22	37.5	115.1	7.2	河北邢台	1533
	17	1967-03-27	38.5	116.5	6.3	河北河间	1688
③	18	1969-07-18	38.2	119.4	7.4	渤海	1965	11.8年	13.5年
	19	1975-02-04	40.7	122.8	7.3	辽宁海城	2340
	20	1976-04-06	40.2	112.1	6.2	内蒙和林格尔	1330
	21	1976-07-28	39.5	118.2	7.8	河北唐山	1874
	22	1979-07-09	31.5	119.3	6.0	江苏溧阳	2001
	*23	1976-08-16	32.6	104.1	7.2	四川松潘	494
	*24	1976-08-23	32.5	104.3	7.2	四川平武	517
	25	1976-09-23	39.8	106.5	6.2	内蒙古磴口	837
④	26	1979-08-25	41.2	108.0	6.0	内蒙古五原	1048	1个月	18.6年
	27	1989-10-19	39.9	113.8	6.1	山西大同	1472
	28	1996-05-03	40.7	109.7	6.4	内蒙古包头	1149
	29	1998-01-10	41.2	114.8	6.2	河北张北	1615
⑤	*30	2008-05-12	31.0	103.5	8.0	四川汶川	570	2.2年
	*31	2013-04-20	30.3	103.0	7.0	四川芦山	606
注：*所示地震为南北地震带中北段MS≥7.0地震，Δ为各地震震中与简化动力源点的距离，Δt为南北地震带中北段MS≥7.0地震与华北地块首发MS≥6.0地震的时间间隔.

下载: 导出CSV 

| 显示表格

根据表 1中我国大陆东部地区历次M_S≥6.0地震的发震时间及震中与动力源点的距离Δ，可以得到1920年以来南北地震带中北段M_S≥7.0地震发生后，华北地块M_S≥6.0地震的时空演化图像，如图 3所示．可以看出，20世纪南北地震带中北段M_S≥7.0地震发生后，华北地块M_S≥6.0地震有随时间从南北地震带附近大体向东迁移的规律．从时空进程来看，大体分为4组，各组地震均以南北地震带中北段M_S≥7.0地震开始，后续华北地区陆续发生一系列M_S≥6.0地震，且前3组中每组地震“序列”的斜率基本相同，每组地震“序列”的迁移视速度约为80 km/a; 第四组地震“序列”的迁移视速度略有不同，约为53.4 km/a.从图 3中还可以看出，1920年以来，华北地区有3次M_S≥6.0地震的发生不符合该迁移规律，依次为1922年9月29日渤海M_S6.5、1976年4月6日内蒙古和林格尔M_S6.2和1996年5月3日内蒙古包头M_S6.4地震(图 3中编号为2，20和28的地震)．

图  3  南北地震带中北段M_S≥7.0地震与华北地块M_S≥6.0地震时空演化图

①—⑤分别对应表1中的组号； 1—31分别对应表1中的地震序号A区： 环渤海地区； B区： 晋冀蒙地区； C区： 南北地震带中北段

Figure  3.  The spatio-temporal evolution of M_S≥7.0 earthquakes in central-northern segment of North-South Seismic Belt and M_S≥6.0 earthquakes in North China block

①—⑤ correspond to the group number in Table 1； 1—31 correspond to the earthquake number in Table 1.A zone: Circum-Bohai area; B zone: Shanxi-Hebei-Inner Mongolia region; C zone: Central-northern segment of North-South Seismic Belt

下载: 全尺寸图片 幻灯片

图 4给出了表 1中南北地震带中北段每组M_S≥7.0地震和后续华北地块M_S≥6.0地震的空间分布图． 可以看出，南北地震带中北段发生M_S≥7.0地震后，后续华北地块M_S≥6.0地震随着时间的推移，由西向东迁移． 首先在南北地震带北段的内蒙地区附近活动，接着在晋、 冀、 蒙及其交界地区，最后在郯庐断裂带的渤海海峡及其附近区域活动． 由此也可以说明，华北地震更多的与印度板块作用下青藏地块向我国大陆东部的挤压作用相关． 图 4d中环渤海地区未发生M_S≥6.0地震，该组相对于动力源点的迁移距离最短，而这组M_S≥6.0地震的持续活动时间为4组中最长的(18.6年)，因此造成该组的迁移速率相对较慢(图 3)．

图  4  南北地震带中北段M_S≥7.0地震和华北地块后续M_S≥6.0地震空间分布图

1—29分别对应表1中的地震序号; 图(a)- -(d)分别对应表1中组号①—④

Figure  4.  The spatial distribution of M_S≥7.0 earthquakes in central-northern segment of North-South Seismic Belt and the subsequent M_S≥6.0 earthquakes in North China block

1—29 correspond to the earthquake number in Table 1； Figs.(a)- -(d) correspond to the group number ①—④ in Table 1

下载: 全尺寸图片 幻灯片

结合图 1，3和4不难发现，南北地震带中北段发生M_S≥7.0地震后，华北地区M_S≥6.0地震首先在地块内部发生，接着在地块边界及北部的西拉木伦断裂和东部的郯庐断裂带附近发生，多数发生在两条断裂的交汇区域(环渤海地区)． 图 3中不符合迁移规律的3次例外地震，分别对应图 4a中的2号地震、 图 4c中的20号地震和图 4d中的28号地震．直观地看，均是发震次序出现了倒置，即相对南北地震带中北段M_S≥7.0地震而言，华北地块M_S≥6.0地震本应是由西至东、由近及远的发生，但上述3次地震却是在距离动力源点更 远处的地震发生后，再次回到距动力源点相对较近的位置发生的．

进一步研究发现，将图 4中华北地块M_S≥6.0地震按每组地震发生的顺序排列，计算相邻两次地震的经纬度差，即顺次用后面地震的经纬度减去前面地震的经纬度．若一次地震在前一次地震的北侧发生，即为北向移步，此时纬向迁移步长为正；反之为南向移步，则纬向步长为负．同样，若一次地震在前一次地震的东侧发生，即为东向移步，此时经向迁移步长为正；反之为西向移步，则经向迁移步长为负． 图 5给出了1920年以来华北地块每次发生M_S≥6.0地震相对于前一次M_S≥6.0地震的经纬向迁移步长． 由图 5a的纬向迁移步长可以看出，华北地块M_S≥6.0地震在纬向上南北跳迁明显(即步长值正负变化均有)，但60%的M_S≥6.0地震发生在(39°N±1.5°)区域内． 而从图 5b的经向迁移步长可以看出，华北地块85%的M_S≥6.0地震发生在前一次地震的东侧(即步长值为正)，距前一次地震平均为2.7°的位置． 因此，从华北地块M_S≥6.0地震经纬向的迁移来看，同样也可反映出来自西南方向的印度板块的推挤作用．

图  5  华北地块M_S≥6.0地震的纬向(a)和经向(b)迁移步长

Figure  5.  The migration step of the M_S≥6.0 earthquakes in North China block in latitude direction (a) and in longitude direction (b)

下载: 全尺寸图片 幻灯片

3.   讨论与结论

板内地震的动态迁移现象反映了板内应力和应变的传递过程． 同样，南北地震带中北段M_S≥7.0地震后华北地块M_S≥6.0地震的东向迁移图像，也反映了印度板块运动作用力向中国大陆的传递过程． 曾融生等(1994)根据地震波层析成像研究认为，印度地壳自南向北向西藏地块的下地壳挤入，在西藏北部塔里木地区受阻不能继续向北挤入，而滇青川地区上地幔的塑性大，从而使西藏地壳由西向东俯冲． 沙海军等(2009)研究认为，中国大陆特大地震的孕育与Ⅰ级和Ⅱ级活动地块(张培震等，2003)的相互作用密切相关，即当中国大陆发生特大地震后，将改变其所在活动地块边界带上的应力分布，进而改变区域构造应力场的分布，对相关活动地块的地震活动起调整作用． 因此可以认为，南北地震带M_S≥7.0 地震的发生主要是由于印度板块推挤运动的影响． 在其作用下，相邻的华北地块的应力场也随之调整，进而产生了华北地块M_S≥6.0地震的系列东向迁移活动．

在分析和讨论地震定向迁移现象的成因时，首先把迁移现象看作是一种“波”的传播过程(Mogi，1968; Savage，1971; Wood，Allen，1973; Scholz，1977)．既然短时间内地壳的岩石圈结构不可能发生很大的变化，因此只要动力源是相对稳定的，那么其作用传播过程也应该是稳定的(马宗晋等，1992)．为了检验华北地震的东向迁移主要来自印度板块运动这一结论，本文尝试改变青藏地块东向作用力代表点的经纬度坐标，再分别计算华北地块M_S≥6.0地震震中到动力源点的距离，并分别绘制各地震震中到动力源点的距离与地震发生时间的关系，如图 6所示． 可以看出，华北地区各组地震“序列” 也大体呈现东向迁移规律，其迁移速率在76km／a附近波动，具有一定的稳定性．依据马宗晋等(1992)的观点，图 3和图 6中相对稳定的迁移图像说明华北地块所受的力源是稳定的．吴云等(1999)基于国家攀登计划建立的全国 GPS网21个测点的复测结果，分析认为印度板块、 太平洋板块和菲律宾板块共同控制着中国大陆的现今地壳运动和形变，但印度板块似乎起主导作用． GPS资料的数值模拟结果显示： 三者对中国大陆驱动的边界作用强度之比约为4 ∶ 1.25 ∶ 1(肖兰喜等，2003； 肖兰喜，2009)； 南北地震带在中国大陆的现今地壳运动和形变过程中起重要调节作用，沿35°N线似乎存在一个调节带； 三大板块对我国大陆的动力作用是动态稳定的，直观表现为推挤速度的稳定． 对比图 3与图 6，发现图 3中地震“序列”的东向近乎匀速的迁移规律似乎更清晰，说明将中国东部应力场的辐射中心点(35°N、 100°E)作为简化的动力源点具有一定的区域代表性． 这一结果在一定程度上与上述观点互为佐证．

图  6  基于不同动力源点的南北地震带中北段M_S≥7.0地震与华北地块M_S≥6.0地震时空演化图

(a) 动力源点(30°N、 100°E); (b) 动力源点(35°N、 105°E)

Figure  6.  The spatio-temporal evolution of M_S≥7.0 earthquakes in central-northern segment of North-South Seismic Belt and M_S≥6.0 earthquakes in North China block based on different power points

(a) Dynamic source point (30°N， 100°E); (b) Dynamic source point (35°N， 105°E)

下载: 全尺寸图片 幻灯片

至于图 3中出现的不符合“序列”迁移规律的3次例外地震，作者也作了大量的分析研究． 由于华北地块受动力环境控制的同时又受局部和阶段性孕震环境的影响，其应力和应变场既体现出整体的一致性(东向平移)，又显示出动态的阶段性变化特征(陈蓥，2011)． 其中值得一提的是，张跃刚和胡新康(2005)通过GPS资料研究认为，华北地块大体有向东平移的趋势，但1995—1996年华北地块边界的相对运动较为突出，鄂尔多斯地块与其它地块间的运动呈现明显的分离，由最初的向东运动逐渐转为向西运动，由此认为向东平移是华北地块水平运动的背景场，某些地块在某些时段出现的背离运动是区域构造运动的暂时调整，可能与区域孕震过程有关． 这或许在某种程度上可以解释1996年5月3日包头M_S6.4地震的“西向迁移”，当然这只是一种可能，其具体原因有待进一步深入研究． 至于1922年渤海M_S6.5地震和1976年和林格尔M_S6.2地震，由于暂时未收集到GPS等相关资料，在此不便妄加评述．

综上，众多的研究结果表明，南北地震带强震对华北地块的应力场和形变场均具有一定的影响(嵇少丞等，2008)，但这种影响不一定以触发地震的方式呈现，也可能是以其它方式释放能量． 换言之，南北地震带强震对华北地块地震的发生有影响，但不是充分必要条件． 马宗晋等(1992)研究认为，地震的定向迁移活动存在“快”、 “ 慢”两种过程，即这种迁移是“慢”过程在定向动力附加作用的影响下，产生的群体优势定向的现象，很难用一次或多次动力作用过程的复杂变化加以解释． 因此上述出现的个别不符合常规的“迁移”图像也正说明了地震迁移的复杂性． 其实，华北地块M_S≥6.0地震也并不一定全是南北地震带强震的影响，太平洋板块的俯冲作用也是其中影响因素之一(吴佳翼等，1979； 孙文斌等，1985)．

上述的地震迁移图像只能说明在一定时间段内(1900年以来)，印度板块对华北地区的作用力更强，因为强震的孕育过程需要几十年—几百年的时间尺度． 陈宇卫等(2003)研究结果表明，印度板块的北向推挤作用力使我国大陆西侧地块有明显向东的运动趋势． 在其影响下，大华北M_S≥6.0地震与西部M_S≥7.5特大地震具有同步活动的特点，但滞后于西部地震活动时间约15年左右． 本文研究结果也表明，华北地块M_S≥6.0地震主要受控于印度板块作用下青藏地块向我国大陆东部挤压的影响，在其作用下，南北地震带中北段发生M_S≥7.0地震后，华北地块M_S≥6.0地震将发生一系列的东向迁移，空间上60%的M_S≥6.0地震的发震地点在(39°N±1.5°)区域内． 根据以往各组地震所得的迁移视速度平均约为80 km/a，华北地块首发M_S≥6.0地震距南北地震带中北段最近一次M_S≥7.0地震的时间间隔约为1个月至11.8年，平均为3.5年. 但实际表 1中第一、 二、 四组“序列”中，华北地块首发M_S≥6.0地震距南北地震带中北段最近一次M_S≥7.0地震的时间间隔均不超过2年. 众所周知，南北地震带中北段近年开始了新一组的M_S≥7.0地震活动，先后发生了2008年汶川M_S8.0和2013年芦山M_S7.0地震，且芦山M_S7.0地震的发生距今已有2.2年． 因此根据本文的研究结果，未来1—3年内华北地块存在发生M_S≥6.0地震的可能，晋冀蒙交界和环渤海及其附近地区值得重点关注．