华北地区小震精定位及构造意义

张广伟 雷建设， 谢富仁 郭永霞 兰从欣

1）中国北京100085中国地震局地壳应力研究所

2）中国北京100081中国地震局地球物理研究所

3）中国北京100080北京市地震局

华北地区小震精定位及构造意义

张广伟1）雷建设1），谢富仁1）郭永霞2）兰从欣3）

1）中国北京100085中国地震局地壳应力研究所

2）中国北京100081中国地震局地球物理研究所

3）中国北京100080北京市地震局

收集了华北地区2001年10月—2009年9月发生在该区内的7 519次地震事件被166个数字化地震台站记录到的74 181条P波和74 465条S波到时资料，采用双差定位法对这些小震进行了重新定位，最终得到了5 511次小震的定位结果.结果表明，重新定位后的结果较好地改进了原定位的精度.从平面上来看，重新定位后的小震更加集中于断裂带附近，条带状更为明显.从震源深度分布来看，研究区内重新定位后小震的震源深度多集中在5—15km的范围内，表明研究区孕震层基本位于地壳的中上部，但不同地震带上小震分布的优势深度存在明显差异.山西地震带上的小震多集中在上地壳，而张家口—渤海地震带陆区部分包括唐山断裂带及邢台地震段上的小震多集中于中上地壳.本研究结果还获得了前人研究结果未揭示的小震分布规律：① 在唐山断裂带北段东侧存在一个近北北东走向的小震活动；② 渤海海域地震带上的小震多处于10km深度处，且有些断层上的小震展示出从地表至10km深度处的垂直分布；③ 太行山前的小震随深度向太行山下方有逐渐变深的趋势，表征太行山山前断裂带东西两侧孕震层有显著差异.这些结果对于认识华北地区地震断层及深部动力学过程均有重要意义.

小地震 地震带 双差定位 华北地区

引言

本文的研究区域为35.5°N—42.5°N，110°E—121°E（图1）.由于受印度板块的北向俯冲和太平洋板块的西向俯冲作用（Yin，2000；Zhao，2004；Liuetal，2004；Lei，Zhao，2005，2006），该区构造环境较为复杂，不仅分布有太行山隆起、燕山隆起和华北平原，而且还存在有很多活动断裂（如邢台断裂带、三河平谷断裂带和唐山断裂带等）和历史强震（如1679年三河平谷8级地震、1966年邢台MS7.2地震、1976年唐山MS7.8地震等）.特别是，近年来该区小震不仅活动频繁，且分布广泛.小震在平面上的分布不仅与断层的走向密切相关，而且在深度上可能与断层的深部几何展布有关，因此小震通常作为认识断层存在和活动的重要标志.所以小震位置的精确程度直接关系到我们对研究区构造活动的认识.

双差定位法（Waldhauser，Ellsworth，2000）被认为是提高地震特别是小震震源位置精度的一种有效手段.由于该方法主要采用两次相近地震到同一台站的走时残差之差，从而能够有效地降低由于地壳速度结构的不确定性而引起的误差.因此，国内外不同学者运用该方法在不同研究区开展了大量研究工作，在地壳孕震层的深度范围和断层走向等方面均获得了富有成效的研究成果（杨智娴，陈运泰，2004；Fukuyamaetal，2003；朱艾斓等，2005；Baietal，2006；李乐等，2007）.就本研究的华北地区而言，也已有很多学者采用该方法开展了类似的研究工作.例如，杨智娴和陈运泰（2004）对1998年1月—1999年3月发生在张北—尚义交界地区MS6.2地震震中附近的87次余震进行重新精定位，揭示出张北—尚义地震的发震构造为北北东向断层；朱艾斓等（2005）对在1980—2000年发生在我国首都圈地区（39°N—41°N，115°E—118°E）的2 098次小震进行重新定位，通过小震线性分布刻画出顺义震群及延怀盆地一带几条倾角较陡的隐伏活动断裂；李乐等（2007）对1980—2004年首都圈及邻区 （35°N—42°N，110°E—120°E）内的地震开展重新定位，发现地震表现为与区域构造更为密切的“井”字形活动分布，并揭示出首都圈地区潜在的地震空区和陡倾角的隐伏断层；于湘伟等（2010b）利用1993—2004年发生在华北地区（35°N—42°N，111°E—120°E）的6 771次地震进行重定位，揭示出新河断裂下方存在两条NE走向、SW倾斜和不同深度的地震活动带.然而，这些研究工作所用的资料尽管有相同之处，但无论是所用资料的时间段还是空间分布均存在较大差别，尤其是他们所用资料包含大量的模拟资料且资料时间仅至2004年.在随后5年时间研究区内密集的数字化地震台站积累了大量高质量到时资料，这些资料对于提高小震定位精度具有重要意义.

研究区内小震定位的精确程度，可能直接影响到我们对断层面相关参数的认识（万永革等，2008；王福昌等，2008），也会直接关系到研究区内地震成像结果的质量.尽管前人对该区做过大量的成像研究工作，但他们均是将走时残差折衷到波速不均匀性和震源参数不确定性等参数反演中，因此无论是获得的速度横向不均匀性还是震源参数的精确程度均具有一定局限（于湘伟等，2003；Huang，Zhao，2004；齐诚等，2006；Leietal，2008，2011）.在本研究中，我们认为有必要收集2001—2009年华北地区数字化记录以来积累的大量高质量到时资料，对该区内的小震进行重新精确定位，以更好地认识地震断层深部结构形态，同时也为进一步成像提供更为可靠的基础资料.

图1 研究区域构造图及历史强震分布Fig.1 Tectonic map and historical strong earthquakes in study region

1 资料与方法

本研究收集了2001年10月—2009年9月9年时间段内发生在华北地区内的7 519次地震事件被166个数字化地震台站所记录到的74 181条P波和74 465条S波到时资料.这些地震台站主要集中于北京和天津地区，在河北省分布虽较为稀疏但却相对均匀，在山西北部、辽宁和山东等地也有少量分布（图2）.为减少本研究中小震震相可能的误判和到时拾取误差，在重新定位前使用一维速度模型（Leietal，2008）和采用HYPODD中计算走时方法，获得理论走时（图3a，c）.该一维速度模型是综合了前人的研究成果（国家地震局《深部物探成果》编写组，1986；嘉世旭等，2005）获得.尽管如此，我们还是比较了前人该区所用的一维速度模型（于湘伟等，2003），结果发现不同速度模型并不会对定位结果产生有意义的影响.最后，将观测走时减去理论走时，获得了我们所需的绝对走时残差.多数P波走时残差集中在±3.0s以内，S波残差分布在－5.0—4.0s，只有少数资料的到时残差超出了±5s；且震中距超过200km，走时残差发生明显畸变.为了有更多可靠资料用于定位，本研究选择残差为±5s、震中距在200km范围内的到时资料（图3b，d），且每个事件至少被5个台站记录到.经过条件筛选，符合条件的事件数为7 357个，P波和S波的到时数分别为67 815条和67 489条（图3）.在重定位前，将资料集中无初始震源深度的小震深度设为0km，这些地震总共263次.由于数量较少，占地震总数不到4%，因而这种假定不会影响我们定位结果的主要结构特点.由小震的初始位置分布（图4a）可知，这些小震的空间分布不均匀，主要集中在张渤（张家口—渤海）地震带、山西地震带、太行山山前断裂带及邢台地震段等主要地震带上，而在其它地区则较为稀疏；在深度方向上，有些小震深度达到30km，且震源呈现层状密集分布.

图4 重新定位前后地震参数对比（a）重新定位前的小震震中平面图、侧面图和小震深度分布频度图.图中圆圈表示获得重定位的地震，“十”字表示未获得重定位的地震，F表示地震频度Fig.4 Comparison of earthquake parameters before and after relocation（a）Small earthquakes before relocation in map and side views，and histogram of the earthquakes with depth.Circles stand for the relocated earthquakes，while crosses stand for the earthquakes that cannot be relocated，Fstands for earthquake frequency

双差法（Waldhauser，Ellsworth，2000）是目前较为通用的提高地震定位精度的一种有效技术方法.该方法主要采用绝对到时残差之差，即相对走时残差或称为双差来修定地震位置.这种相对走时残差的使用主要是基于两个震源之间的距离远小于事件到台站的距离，认为两个事件传播到台站的射线路径几乎相同，从而可以有效地降低由于速度模型的不确定性对定位结果的影响，也会有意义减少到时资料拾取时的系统误差.

两个相邻事件i和j到相同台站k的走时残差之差（双差drijk）可表示为

图4 重新定位前后地震参数对比（b）重新定位后的小震震中平面图、侧面图和小震深度分布频度图；（c）重定位前后深度分布对比.图中圆圈表示获得重定位的地震，F表示地震频度Fig.4 Comparison of earthquake parameters before and after relocation（b）Small earthquakes after relocation in map and side views，and histogram of the earthquakes with depth；（c）Comparison of focal depths between before and after relocation.Circles stand for the relocated earthquakes，Fstands for earthquake frequency

由于我们假定地震i与地震j到台站的射线路径几乎相同，即∫则式（1）可简化为

另外，双差drikj也可以表示为

最终，我们可以将式（4）反演问题归结为一般线性方程组

式中，A为雅可比矩阵，y为待求解未知量（发震时刻、震中经度、纬度和震源深度），b为观测量（到时差之差）.实际计算中所用的地震事件和观测资料均较多，因而式（5）中A为大型稀疏矩阵，因此我们采用共轭梯度法（Paige，Saunders，1982）求解式（5），最终得到阻尼最小二乘解.有关该方法的细节，可参阅Waldhauser和Ellsworth（2000）文章.

2 相关参数的选取

采用双差法对华北地区的小震进行重新定位之前，我们必须对研究区内的小震进行配对，即要求配对的地震事件在一定的空间范围内能够搜索到至少一次地震事件.这个空间范围主要根据研究区内的地震密集程度和所需定位精度来确定，即通常所说的搜索半径.该搜索半径对定位结果有一定影响：若该值过小，事件簇中能够建立联系的事件对就减少，能够定位的事件数也会随之减少；反之，若该值过大，超出了震源区的速度结构异常，便会直接影响定位结果的精度.根据本研究区内的地震分布特点，为得到更高精度的定位结果，且有更多地震能获得重定位，我们将初始确定事件对的搜索半径设为10km，发现有7 154个事件能够进行配对并参与重新定位.

由于S波到时资料拾取的精度相对于P波要低些，因此为减少S波资料的拾取误差和保证S波的走时残差发挥足够作用，这里我们将P波到时残差的权重取为1.0，而S波到时残差权重取为0.5.以观测报告中的定位结果作为初始值，对地震进行重定位.在重定位过程中进行两轮运算：在第一轮的8次迭代中，采用标准偏差的5倍作为截断值对地震数据进行剔除；在第二轮的10次迭代中，采用标准偏差的3倍作为截断值对地震数据进行剔除.在迭代过程中，对不同的距离阀值进行反复比较，结果发现选取9km时定位结果最佳.当均方根残差达到由数据的噪音水平决定的阀值或迭代次数达到最大值时（Waldhauser，Ellsworth，2000），迭代过程停止.最终我们获得了5 511个高质量的地震事件震源位置（图4b），其均方根残差由重定位前的0.44s下降到重定位后的0.18s.在东西、南北及深度上的平均相对误差分别为0.53，0.49km和0.57km，且三方向相对误差均在1.0km以内的事件达4 485个，占重定位事件的80%.尽管定位后的一些事件的相对误差超过1.0km，但这些事件空间分布较为零散，说明双差法对地震簇定位具有明显优势.

3 结果与分析

图4b展示出我们获得的重新定位后的5 511个地震事件震中分布的平面图和经度方向及纬度方向的侧面图.尽管从平面图上重定位前后地震位置在多数地区难以观察到明显变化，但在较密集的地震条带上的地震位置变化还是很明显的，如唐山地震带南段和邢台地震段重定位后地震表现更为聚集；从侧面图也可以清楚地看出重定位前后明显变化（图4a，b）.重定位前地震多集中在5—10km深度范围内，特别在10km深度上地震更为集中；但在重定位后地震多分布在5—15km深度范围内，且震源的深度分布形态更接近于正态分布（图4c）.由重定位后地震分布平面图（图4b）可以看出，地震主要集中在山西地震带、邢台—文安地震带和张渤地震带等地震条带上，而在其它地区则相对稀疏.在山西地震带，地震多集中在盆地边缘，表征山西地震带的发震构造背景是由一系列的断陷盆地组成（刘光勋，阎凤忠，1995）；邢台—文安地震带中的邢台地震段小震主要集中分布在新河断裂与束鹿盆地的交界带；张渤地震带地震多集中分布于NE向和NW向断裂带上及一些断裂交汇处，其中在夏垫及唐山断裂带上小震的条带状分布非常明显.为更清晰地认识和理解该区内发震断层性质，本文主要集中于山西地震带、邢台—文安地震带、张渤地震带陆区段、太行山山前地震带和渤海海域地震带等几个纵剖面来深入探讨研究区内小震的空间分布特征以及它们所蕴涵的深层次含义.需要说明的是，在重定位过程中有较零散分布的2 008次地震被去除或未能获得重定位，与双差法优势于成簇地震定位有关.另外，这些丢失的地震多分布在断裂带以外（图4a），因而不会影响本文着重讨论断裂带的结构形态.

3.1 山西地震带

山西地震带位于稳定的鄂尔多斯地块与华北盆地之间，受太平洋板块西向俯冲和印度板块北向逆冲的共同作用，使得该地震带或构造带为一张性地震活动带（Yin，2000；Liuet al，2004；Zhao，2004；Lei，Zhao，2005，2006）.由于该地震带所处的特殊地理位置和其所具有的特定构造特点，因而无论是地质学还是地球物理学的研究均倍受关注（Heetal，2004；王秀文等，2010；李传金等，2010）.地质学研究结果表明，该地震带为一条北北东向的地堑型活动构造带，主要受新生代断陷带控制（刘光勋，阎凤忠，1995）；而地球物理学研究表明，山西地震带上的震源机制解多为走滑正断层，地堑两侧的应力作用多以张性为主，可能表明深部有热物质上涌（Heetal，2004；王秀文等，2010）.该断裂带为我国强震活动集中带，以强度大、频度高和破坏严重为特点.以上世纪末为例，该地震带曾发生过1989年大同MS6.1地震，1991年忻州MS5.1地震和大同MS5.8地震，以及1999年大同—阳高MS5.6地震.这些地震的发生可能标志着山西地震带的地震活动进入一个新的活动期（王秀文等，2010）.因此，近年来震群的频繁活动，为了解山西地震带的活动特征提供了有利条件.

从重定位后的震源深度分布（图5a）可以看出，小震主要发生在上地壳，震源深度集中在5—12km.表明山西地震带的发震层深度较浅，而重定位前的小震零散的分布在0—15km的范围内.小震的分布在六棱山山前、系舟山山前等地区表现为较垂直分布，可能暗示了断裂带的倾角较陡.从小震平面分布可以看出，在太原盆地北缘、忻定盆地系舟山山前断裂带以及大同盆地东南缘等地区小震分布较为密集，而在太原盆地的交城断裂和太谷断裂带小震分布较为稀疏.小震主要集中在太原凸起区，该段小震密集主要与田庄断裂和新城—亲贤断裂的活动密切相关.这样的震群分布特点与前人研究结果较为一致（侯廷爱等，1997）.系舟山山前断裂是忻定盆地东南边界断裂，沿该断裂带的小震有着与断裂走向NE向相一致的走向.大同盆地附近的地震分布不均匀，西北边界的口泉断裂带上小震分布较少，但在大同盆地东南缘的六棱山北麓断裂带附近小震较多，其走向主要为NE、NNE向，可能暗示出了六棱山北麓断裂带的走向.

图5 重定位后小震震源深度纵剖面图.（a）山西地震带；（b）邢台—文安地震带；（c）张渤地震带.除剖面上方的地形外，本图与图6—10中的纵横轴均成比例Fig.5 Vertical profiles of the relocated small earthquakes.（a）Shanxi seismic zone；（b）Xingtai－－Wenan seismic zone；（c）Zhang－－Bo seismic zone.No vertical exaggerationexcept for the topography above vertical cross sections in Fig.6—10

3.2 邢台—文安地震带

邢台—文安地震带位于华北平原，根据地震密集程度可分为邢台段和文安段.其中，邢台段横穿束鹿盆地、晋县盆地、新河凸起和宁晋凸起，走向大致为NE向的地震带（徐锡伟等，2000）.该断裂带曾发生过1966年邢台MS6.8和MS7.2地震，造成8 000余人死亡和巨大经济损失.研究认为，邢台段中的新河断裂是该次地震的发震断裂（许桂林，范国胜，1986；徐杰等，1988）.另外，2006年文安MS5.1地震也发生在邢台—文安地震带中的文安段上，尽管在震中附近没有断裂对应，但在震中以东方向小震较多，且与地震带走向基本一致（刁桂苓等，2006）.这些中强震的发生，可能说明邢台—文安地震带近年来较为活跃.通过重新定位，发现重定位后邢台段小震较为密集，震源主要集中分布在5—20km的范围内；到北部的文安段，震源深度有所增加，表明该段的发震层有所加深（图5b）.Lei等（2008，2011）认为文安地震前震源区确实存在明显的低速异常，可能与西太洋板块西向深俯冲至太行山前并在地幔转换带中滞留、脱水和软流圈热物质上涌等（Zhao，2004；Lei，Zhao，2005，2006）动力学过程相关.

为更清晰地比较重定位前后邢台段与文安段震源区附近小震分布的差异性，以及分析定位后小震所蕰含的物理意义，我们将定位前后的小震分布分别展示于图6和图7.由图6中AA′剖面可知，邢台段震源区重定位前小震零散地分布在0—30km范围内（图6a），而重定位后小震多集中在5—18km之间（图6b），表明该深度层为邢台段的优势发震层；由BB′剖面可知，重定位后的小震分布在新河断裂呈现出上陡下缓的“铲状”形态（图6d），而重定位前的小震分布并未显示出这样的分布特征（图6c）.这些结果均表明，小震的重定位可以清晰地反映出断层深部结构形态和活动态势.

图6 邢台地震带小震重定位前（a，c）、后（b，d）震源深度纵剖面对比图.六角星表示1966年邢台MS7.2地震Fig.6 Comparison between small earthquakes before（a，c）and after（b，d）the relocation around Xingtai seismic zone along the vertical cross sections.Hexagram represents the 1966Xingtai MS7.2earthquake

通过对比文安段以东的南北向纵剖面和穿过文安主震震源区东西向纵剖面的重定位前后结果（图7），可以看出重定位前后小震分布的深度范围变化不大，但重定位后能清晰分辨出文安震源区附近3个明显的小震群：第一个位于文安地震的东北方向，表现为明显的双层发震层（图7b）；第二个位于文安地震正东方向，表现为从地表近垂直连续延伸至大约25km的深度（图7d）；第三个恰好位于文安主震震源区上方，但该震群小震最少（图7d）.这些结果说明，文安地震震源区存在强烈的横向不均匀性.

3.3 张渤地震带

张渤地震带位于华北平原与燕山山脉之间的过渡带.该地震带由20多条断续分布的断裂组成，单条断裂长20—40km，倾向不一，倾角较陡，具有正断兼左旋走滑的运动方式，被认为是一条第四纪以来的活动断裂带（方颖，张晶，2009），华北北部地区的强震活动明显地受该构造断裂带所控制.重力资料结果表明，该断裂带的中部和东南部为不同走向线性异常的交接带，西北部是大兴安岭—太行山重力梯度扭曲带（徐杰等，1998）.航磁资料的分析结果表明，沿断裂带主要表现为NNW至近EW向的局部正异常及少量NW向局部负异常和NE向局部正负异常组成的NW向航磁异常条带，也是不同走向的磁场交接带（徐杰等，1998）.从NW与NE向两组断裂的相互交切组合关系来看，张渤带自西而东分为张北—怀来、南口—三河、天津—塘沽、渤海中部和蓬莱—烟台等5个构造交汇段（徐杰等，1998）.直达P波和莫氏面反射波PmP联合高分辨率成像结果显示，张渤带自东南向西北有由宽变窄的趋势（Leietal，2008）.这一结果与前人分析的地质、地貌和小震活动结果（北京市地震局①北京市地震局.2006.北京市2007年度地震趋势研究报告.）所显示的变化模式相类似.

图7 文安地震带小震重定位前（a，c）、后（b，d）纵剖面对比图.六角星表示2006年文安MS5.1地震Fig.7 Comparison of small earthquakes before（a，c）and after（b，d）the relocation around Wenan seismic zone along the vertical cross sections.Hexagram represents the 2006Wenan MS5.1earthquake

张渤地震带小震活动最为频繁，这为我们认识该地震带提供了有利条件.通过对比重定位前后该地震带上小震的分布，可以看出重定位前后小震的震源深度有着十分明显的变化，重定位前震源分布在0—30km，重定位后主要集中在5—20km（图5c），表明震源深度优势分布主要集中在中上地壳.从重定位后的结果可以清楚地看出，张渤地震带自NW向SE小震活动不仅在水平方向上有明显分段不连续特点，即在115.3°E和116.8°E处均出现了不同程度的明显不连续现象，而且在深度方向也表现出显著的分段特点.在西北部的涿鹿、延怀盆地一带较浅，而在中部的昌平顺义一带则明显加深，到夏垫、宝坻和唐山一带又逐渐变浅.这种由西向东的变化趋势，与前人的研究结果（朱艾斓等，2005）相类似，但在昌平顺义一带的震源深度分布有所差异.前人的结果显示出该段小震震源深度可达25km左右，而本研究震源深度却仅为20km.这一结果与层析成像结果显示的该段低速异常体顶部界面深度（Leietal，2008，2011）相类似.说明本研究的结果可能更符合该段的孕震构造，因而认为本研究相对前人的定位结果有了一定意义的改进.

为进一步认识张渤地震带，我们还详细地分析了唐山断裂带的小震分布（图8）.在重新定位前，小震沿断裂带方向无明显变化趋势（图8a），在唐山断裂带北段垂直于断裂方向也无法判定有两条断裂（图8e），但重定位后小震分布呈现出明显的结构特点，且主要分布在20km深度以上（图8b，d，f）.沿断层走向，多数小震集中在唐山断裂带北段，在宁河段震源深度较大，而在唐山主震段小震震源深度较浅（均在15km深度之上），呈现出明显上隆的分布格局（图8b）.前人对唐山断裂带的研究成果（孙若昧，刘福田，1995）也揭示出类似的发震构造.人工地震测深结果显示，在唐山主震附近中地壳15—20km深度范围内存在一低速层，唐山主震震中的位置便位于这个低速层的顶面附近（赖晓玲等，2006）；天然地震层析成像结果也显示出唐山主震就发生在高低速异常边界，且接近于高速异常区（Leietal，2008，2011）.在垂直于断层走向的纵剖面上，重定位后小震的分布相对之前更加集中于断裂带（图8d）.唐山断裂带北段呈现出两条明显的断层，且东边断层可能略深于西边断层（图8f）.最近研究结果显示出的唐山地震带存在3个明显震群分布特点，即唐山、卢龙（迁安）和滦县震群（于湘伟等，2010a）.本文定位结果不仅证实了该结构特点，还揭示出了在唐山断裂带北段东侧存在一个近北北东走向的断层（图8f）.

图8 唐山断裂带重定位前（a，c，e）、后（b，d，f）震源剖面对比图.六角星表示1976年唐山MS7.8地震Fig.8 Comparison of small earthquakes before（a，c，e）and after（b，d，f）the relocationalong three vertical cross sections around Tangshan seismic zone along the verticalcross sections.Hexagram represents the 1976Tangshan MS7.8earthquake

万永革等（2008）采用首都圈数字地震台网资料编制的唐山地震观测报告，运用双差法得到1 896个高精度的定位结果，然后应用这些小震资料反演获得了唐山主震北段的断层参数（表1）.为证实我们结果的可靠性，本研究采用万永革等（2008）反演断层参数的算法，使用与其类似的地震事件，获得了该区一个新的断层参数（表1）.结果显示，本研究结果与前人的结果有较好的一致性，但本研究结果的标准差明显降低，说明本研究获得的华北地区小震的定位精度相对于前人有明显提高.这显然对进一步深入理解研究区内断层结构形态有着重要意义.

表1 唐山断裂北段断层参数反演结果与万永革等（2008）结果对比Table 1 Comparison between previous and present result on northern segment of the Tangshan fault

图9 太行山山前断裂带重定位前（a，c，e）、后（b，d，f）小震分布纵剖面对比图Fig.9 Comparison between the small earthquakes before（a，c，e）and after（b，d，f）the relocation around Taihang mountain front in vertical cross sections

3.4 太行山山前及渤海海域地震带

太行山山前断裂带位于华北断块区的中部，为太行山断块隆起区与华北平原裂陷盆地区的构造分界线，由黄庄—高丽营、徐水、保定—石家庄、邯郸、汤东和汤西等10多条NE－－NNE向断裂组成（徐杰等，2000）.从小震的分布可以看出，在保定段、宁晋段及磁县段小震较多，在此我们由北而南展示出这3个地区的小震分布纵剖面（图9）.通过重定位前后的对比，可以看出重定位后小震更加密集（图9b，d）.在保定段小震近乎垂直的密集分布暗示出该段存在两个倾角较为陡立的断层；在宁晋段及磁县段，断裂带西侧的小震展现出向太行山下方震源深度逐渐加深的趋势.图9中的虚线是采用CRUST2.0模型（Bassinetal，2000）所展示的莫霍面深度.该模型是由CRUST5.1模型（Mooneyetal，1998）改进后获得.尽管该模型的莫霍面分辨率不高，但在太行山东西两侧仍表现出了由东向西逐渐变深的趋势，这与曾融生等（1995）认为莫霍面深度由断裂带东侧的32—34km增加到西侧的40—42km相一致.太行山山前断裂带两侧地壳厚度的显著变化，可能导致了断裂带东西两侧孕震层深度的变化，这可能是造成震源深度逐渐加深的主要原因之一.尽管野外地质调查均显示，太行山山前断裂多数表现为正断层类型，但可能由于晚第三纪以来的区域构造作用，主压应力轴逐渐转变为ENE－－WSW（徐杰等，2000）.这与震源机制解联合滑动资料及水压致裂结果反演获得的近EW向水平主压应力方向（Wan，2010）较为接近.因此，认为研究区内当前的近EW水平压应力，也有可能为导致太行山山前断裂带西侧震源深度加大的原因之一.

渤海位于张渤地震带的最东段，既有NNE向右旋走滑断裂带，又有WNW向左旋走滑断裂带.渤海海域的地震活动可能与由太平洋板块西向俯冲而形成的ENE向挤压应力和渤海上地幔隆起产生的构造应力的共同作用密切相关（苗庆杰等，2010）.前人研究结果表明，渤海地区上地幔上隆，其上隆中心最薄处地壳只有28km厚（赖晓玲等，2007）.其下方可能存在地幔热柱（滕吉文等，1997）.渤海海域自有记录以来，曾发生过9次6级以上强震，其中7级以上地震有4次，最近一次为1969年7月18日MS7.4地震（滕吉文等，1997）.初始小震深度基本集中在10km这个深度上（图10a，c，e），重定位以后有少量小震从地表延伸至10km深度（图10b，d，f）.说明该处断层在这个深度范围均表现为活动状态.本研究结果显示发震层位于上地壳脆性层内.环文林等（1989）对渤海地区小震的定位结果显示，渤海地震均发生在地壳以内，但震源深度多在地壳上部的10—20km范围内.这些结果表明，本项研究有意义地改进了前人定位结果，这种改进可能与我们采用双差定位法，从而克服了速度模型不确定性和减少了到时资料的读取误差密切相关.

图10 渤海海域重定位前（a，c，e）、后（b，d，f）小震纵剖面对比图Fig.10 Comparison between small earthquakes before（a，c，e）and after（b，d，f）the relocation along vertical cross sections in Bohai Sea region

4 讨论与结论

本研究运用双差定位法对我国华北地区2001—2009年发生的地震进行了重新定位，获得5 511个定位后的地震事件.这些事件在东西、南北及深度方向上的平均相对误差分别为0.53，0.49km和0.57km，均方根残差由重定位前的0.44s下降到重定位后的0.18s.研究区内震源深度集中分布在5—15km，进一步证实了研究区的发震层在中上地壳.重定位后小震分布更加集中于断裂带，呈现出更加明显的条带状分布.

张渤地震带陆区段地震多集中在20km深度之上，但无论是在深度方向上还是横向方向上均呈现出明显的不均匀性，这与前人的地球物理结果相类似.然而，本项研究结果显示，渤海海域的地震多集中在10km深度以上，这无疑扭转了前人认为多在10—20km深度范围的认识.山西地震带的小震主要分布在上地壳，且小震集中在盆地边缘，表征山西地震带的发震构造是由一系列的断陷盆地组成.太行山山前断裂西侧的小震随深度向太行山下方有逐渐加深的趋势，反映了太行山山前断裂东西两侧地壳厚度的显著变化.

重定位后邢台—文安地震带的邢台段小震在新河断裂处呈现出上陡下缓的“铲状”形态；文安震源区附近存在3个明显的小震群：第一个位于文安地震的东北方向，表现为明显的双层发震层；第二个位于文安地震正东方向，表现为从地表近垂直连续延伸至大约25km深度；第三个恰好位于文安主震震源区上方，但该震群小震最少.这些结果可能说明邢台—文安地震带存在着强烈横向不均匀性.

通过使用与前人相同的反演方法和基本相同的资料集，发现唐山断裂带北段的断层参数与前人基本相同，但标准差明显降低.说明本项研究结果相对前人的结果有了一定意义上的改进.这是因为本研究还揭示出了在唐山断裂带北段东侧存在一个近北北东走向的小震群.因此，本项研究进一步深化了我们对该区地震断层深部结构形态的认识，同时也为获得高质量的成像结果提供了更为可靠的基础资料.

刁桂苓，张四昌，孙佩卿，孙佩卿，张跃刚，高景春，毛国良.2006.2006年7月4日文安5.1级地震［J］.地震地质，28（3）：497－－502.

方颖，张晶.2009.张家口—渤海断裂带分段活动性研究［J］.地震，29（3）：136－－140.

国家地震局《深部物探成果》编写组.1986.中国地壳上地幔深部探测成果［M］.北京：地震出版社：1－－407.

侯廷爱，程新原，孟雁英，阎维彰.1997.太原盆地震群的几个特征［J］.山西地震，（3）：12－－17.

环文林，汪素云，常向东，武宦英，李群.1989.渤海地震活动特征［J］.地震研究，12（1）：1－－9.

嘉世旭，齐诚，王夫运，陈棋福，张先康，陈颙.2005.首都圈地壳网格化三维结构［J］.地球物理学报，48（6）：1316－－1324.

赖晓玲，孙译，刘志.2006.华北强震区地震测深研究［J］.大地测量与地球动力学，26（1）：55－－62.

赖晓玲，李松林，孙译.2007.渤海及邻区3次7级以上地震的深部构造背景［J］.大地测量与地球动力学，27（1）：31－－33.

李传金，徐佩芬，孙勇军，贾建喜.2010.用远震接收函数研究山西地区地壳厚度变化：“晋中坡折带”及其地质意义探讨［J］.地球物理学报，53（5）：1143－－1148.

李乐，陈棋福，陈颙.2007.首都圈地震活动构造成因的小震精定位分析［J］.地球物理学进展，22（1）：24－－34.

刘光勋，阎凤忠.1995.从山西地震带看大同—阳高地震［J］.山西地震，（1）：3－－6.

苗庆杰，周翠英，马玉香，苗春兰，李霞，刘霞.2010.1969年渤海7.4级地震研究述评［J］.华北地震科学，28（1）：9－－15.

齐诚，赵大鹏，陈颙，陈棋福，王宝善.2006.首都圈地区地壳P波和S波三维速度结构及其与大地震的关系［J］.地球物理学报，49（3）：805－－815.

孙若昧，刘福田.1995.京津唐地区地壳结构与强震的发生：Ⅰ.P波速度结构［J］.地球物理学报，38（5）：599－－607.

滕吉文，张中杰，张秉铭，杨顶辉，万志超，张慧.1997.渤海地球物理场与深部潜在地幔热柱的异常构造背景［J］.地球物理学报，40（4）：468－－480.

万永革，沈正康，刁桂苓，王福昌，胡新亮，盛书中.2008.利用小震分布和区域应力场确定大震断层面参数方法及其在唐山地震序列中的应用［J］.地球物理学报，51（3）：569－－583.

王福昌，万永革，胡顺田.2008.粒子群算法在主震断层面参数估计中的应用［J］.地震研究，31（2）：149－－154.

王秀文，宋美琴，杨国华，李斌，杨晓泓.2010.山西地区应力场变化与地震的关系［J］.地球物理学报，53（5）：1127－－1133.

徐杰，方仲景，杨理华.1988.1966年邢台7.2级地震的构造背景和发震构造［J］.地震地质，10（4）：51－－59.

徐杰，宋长青，楚全芝.1998.张家口—蓬莱断裂带地震构造特征的初步探讨［J］.地震地质，20（2）：146－－154.

徐杰，高战武，宋长青，孙建宝.2000.太行山山前断裂带的构造特征［J］.地震地质，22（2）：111－－122.

徐锡伟，于贵华，王峰，顾梦林，孙振国，刘保金，尤惠川.2000.1966年邢台地震群的发震构造模型：新生断层形成？先存活断层摩擦粘滑？［J］.中国地震，16（4）：364－－378.

许桂林，范国胜.1986.华北盆地北部断裂活动特征及其与地震的关系［J］.华北地震科学，4（1）：37－－47.

杨智娴，陈运泰.2004.用双差地震定位法再次精确测定1998年张北—尚义地震序列的震源参数［J］.地震学报，26（2）：115－－120.

于湘伟，陈运泰，王培德.2003.京津唐地区中上地壳三维P波速度结构［J］.地震学报，25（1）：1－－14.

于湘伟，陈运泰，张怀.2010a.京津唐地区中小地震重新定位［J］.地震学报，32（3）：257－－269.

于湘伟，张怀，陈运泰.2010b.华北地区地震重新定位结果分析［J］.大地测量与地球动力学，30（2）：29－－33.

曾融生，孙为国，毛桐恩，林中洋，胡鸿翔，陈光英.1995.中国大陆莫霍界面深度图［J］.地震学报，17（3）：322－－326.

朱艾斓，徐锡伟，胡平，周永胜，林元武，陈桂华，甘卫军.2005.首都圈地区小震重新定位及其在地震构造研究中的应用［J］.地质评论，51（3）：268－－274.

Bai L，Kawasaki1I，Zhang T Z，Ishikawa Y.2006.An improved double－difference earthquake location algorithm using sP phases：application to the foreshock and aftershock sequences of the 2004earthquake offshore of the Kii peninsula，Japan（MW＝7.5）［J］.EarthPlanetsSpace，58：823－－830.

Bassin C，Laske G，Masters G.2000.The current limits of resolution for surface wave tomography in North America［J］.EOSTransAGU，81：F897.

Fukuyama E，Ellsworth W L，Waldhauser F，Kubo A.2003.Detailed fault structure of the 2000western Tottori，Japan，earthquake sequence［J］.BullSeismSocAmer，93（4）：1468－－1478.

He J K，Cai D S，Li Y X，Gong Z.2004.Active extension of the Shanxi rift，north China：does it result from anticlockwise block rotations？［J］.TerraNova，16（1）：38－－42.

Huang J L，Zhao D P.2004.Crustal heterogeneity and seismotectonics of the region around Beijing，China［J］.Tectonophysic，385（1－4）：159－－180.

Lei J S，Xie F R，Lan C X，Xing C Q，Ma S Z.2008.Seismic images under the Beijing region inferred from P and PmP data［J］.PhysEarthPlanetInter，168（3－4）：134－－146

Lei J S，Zhao D P，Xie F R，Liu J.2011.An attempt to detect temporal variations of crustal structure in the source area of the 2006Wen－An earthquake in North China［J］.JAsianEarthSci，40（4）：958－－976.

Lei J S，Zhao D P.2005.P－wave tomography and origin of the Changbai intraplate volcano in Northeast Asia［J］.Tectonophysics，397（3－4）：281－－295.

Lei J S，Zhao D P.2006.Global P－wave tomography：On the effect of various mantle and core phases［J］.PhysEarth PlanetInteri，154（1）：44－－69.

Liu M，Cui X J，Liu F T.2004.Cenozoic rifting and volcanism in eastern China：A mantle dynamic link to the Indo－Asian collision？［J］.Tectonophysics，393（1－4）：29－－42.

Mooney W，Laske G，Masters G.1998.CRUST 5.1：A global crustal model at 5°×5°［J］.JGeophysRes，103（B1）：727－－747.

Paige C C，Saunders M A.1982.LSQR：Sparse linear equations and least squares problems［J］.ACMTransactionson MathematicalSoftware，8（2）：195－－209.

Waldhauser F，Ellsworth W L.2000.A double difference earthquake location algorithm：method and application to the North Hayward fault，California［J］.BullSeismSocAmer，90（6）：1353－－1368.

Wan Y G.2010.Contemporary tectonic stress field in China［J］.EarthquakeSci，23（4）：377－－386.

Yin A.2000.Mode of Cenozoic east－west extension in Tibet suggesting a common origin of rifts in Asia during the Indo－Asian collision［J］.JGeophysRes，105（B9）：21745－－21759.

Zhao D P.2004.Global tomographic images of mantle plumes and subducting slabs：insight into deep Earth dynamics［J］.PhysEarthPlanetInteri，146（1－2）：3－－34.

张广伟 中国地震局地壳应力研究所.2008年吉林大学应用地球物理专业毕业，获学士学位；2011年中国地震局地壳应力研究所毕业，获理学硕士学位.现主要从事小震精定位等方面的研究工作.

Precise relocation of small earthquakes occurred in North China and its tectonic implication

Zhang Guangwei1）Lei Jianshe1），Xie Furen1）Guo Yongxia2）Lan Congxin3）
1）InstituteofCrustalDynamics，ChinaEarthquakeAdministration，Beijing100085，China
2）InstituteofGeophysics，ChinaEarthquakeAdministration，Beijing100081，China
3）EarthquakeAdministrationofBeijingMunicipality，Beijing100080，China

In the present study we collected 74 181P－wave and 74 465S－wave arrival times from 7 519local earthquakes in North China recorded by 166seismic stations during October 2001to September 2009，and relocated these earthquakes by using double－difference location algorithm，finally obtained 5 511re－located earthquakes.Our result shows that the locations of the relocated earthquakes have been improved significantly over those of original ones.In map view the relocated earthquakes are more concentrated on obvious strips around the faults，while most of them occurred at 5—15km depths，suggesting the seismogenic layer in the region being in the mid－upper crust.However，different seismic zones show obvious contrasts in predominant earthquake focal depths.Most small earthquakes around the Shanxi seismic zone occurred in upper crust，while those in continental part of the Zhang－－Bo seismic zone and in Xingtai seismic zone mainly occurred in mid－upper crust.The present study obtained some results not revealed by previous studies：①Small earthquakes occurred in NNE－－SSW trend on east side of the northern segment of the Tangshan seismic zone.② Most relocated earthquakes in sea region of the Bohai sea zone are at about 10km depth，with a few earthquakes having a depth from a few to 10km.③Small earthquakes around the fault zone in front Taihang mountain gets deeper westward under the mountain，suggesting a deeper seismogenic layer beneath the mountain.These results are helpful for understanding fault seismic activity and deep mantle dynamics of North China region.

small earthquake；seismic zone；double－difference location algorithm；North China region

10.3969/j.issn.0253－3782.2011.06.001

P315.2

A

张广伟，雷建设，谢富仁，郭永霞，兰从欣.2011.华北地区小震精定位及构造意义.地震学报，33（6）：699－－714.

Zhang Guangwei，Lei Jianshe，Xie Furen，Guo Yongxia，Lan Congxin.2011.Precise relocation of small earthquakes occurred in North China and its tectonic implication.ActaSeismologicaSinica，33（6）：699－－714.

中央级公益性科研院所基本科研业务专项（ZDJ2007－1）和北京市自然科学基金项目（8092028）资助.

2010－11－24收到初稿，2011－04－26决定采用修改稿.

e－mail：jshlei＿cj＠hotmail.com