引言

断层面参数是描述断层构造和地震机制的重要参数,该参数既体现了断层构造的性质,又为地震发震构造判定提供依据,在地球物理学、地质学等学科中具有重要地位。确定断层参数的方法主要有地质学方法、地球物理学方法等,其中地质学方法主要是通过浅层的地质信息推测断层产状,虽然比较直观,但断层浅部出露与深部在构造形态上可能存在较大差异,存在一定的局限性;地球物理学方法包括天然地震法、地震勘探、电法勘探等,其中地震勘探、电法勘探等物探方法多用于城市活断层探测,天然地震法确定地震断层面参数的方法主要有体波和面波联合法、面波波形及小震确定层面参数等。地球物理学方法所采用的资料包含深部信息,能更好地刻画深部断层形态。

现今,随着地震台网的密集以及地震定位精度的提高,越来越多的学者开始使用小震分布确定断层面参数(王鸣等,1992Ouillon等,2008万永革等,2008王福昌等, 2012, 2013盛书中等,2014),但研究所使用的地震数据均为双差定位后的精定位地震目录数据。自2009年1月1日开始,由国家地震台网和31个区域地震台网组成的覆盖中国的地震监测台网初步建成,通过统一编目系统(黄文辉等,2016),实现了国家地震台网和区域地震台网的统一编目。通过改变国家地震台的数据上传方式,形成新的地震目录编辑方法,统一将地震事件进行震相删选和重新定位,大大提高地震参数的测定精度。本文尝试基于实例地震数据,直接使用中国地震台网统一地震目录来确定断层面参数,并将研究结果与前人利用精定位数据得到的结果进行比较,以验证该方法的可行性。

1976年河北唐山发生了MS 7.8地震,许多学者对唐山地震的发震构造及破裂做了大量研究(虢顺民等,1977Butler等,1979陈运泰等,1979李钦祖等,1980张之立等,1980王景明等,1981Nábělek等,1987尤惠川等,2002),且震后震源区小震频发,故该地区积累了丰富的小震资料。随着地震定位方法在地震学中的应用及小震定位精度的提高,对唐山地区小震分布的研究也越来越多,张宏志等(2008)采用双差定位法对唐山地震震区中小地震重新定位,重新定位后的唐山断裂南段走向为NNE向,断裂北段转为NE向,滦县断裂附近区域地震分布呈“丁”字形,宁河断裂地震分布无明显优势方向;万永革等(2008)将小震分布确定断层面参数的方法应用于唐山地震序列,把唐山地震序列分为宁河断裂段、唐山断裂南段、唐山断裂北段、卢龙断裂段和滦县断裂段,对唐山地震序列定量研究,利用模拟退火算法和高斯迭代算法相结合的算法给出各段地震断层面的走向、倾角、位置及其误差。因此,本文选择地震资料丰富且研究程度较高的唐山地震进行研究及实例分析。

1 资料情况

本文使用的地震数据来源于国家科技基础条件平台——国家地震科学数据共享中心1提供的中国地震台网统一地震目录,研究区域为117.2°—119.2°E、39°—40°N(图 1)。选取2009年1月1日—2018年10月10日发生的4250次地震事件,震源深度集中分布在3—20km,震级主要为ML 2.5以下(图 2)。


图 1 唐山地区小震分布 Fig. 1 Distribution of small earthquakes in the Tangshan area

图 2 小震深度(a)和震级分布图(b) Fig. 2 Histogram of depth (a) and magnitude (b) of small earthquakes

1 http://data.carthquake.cn

2 研究方法

本研究采用万永革等(2008)提出的利用小震拟合断层面参数的方法。该方法基于2个基本假设:小震均发生在断层面及其附近区域;发震断层面可近似为1个或多个平面。基于小震震源位置到该平面的距离最小,建立求解断层面参数的数学模型,采用模拟退火全局搜索和高斯-牛顿局部搜索相结合的方法,给出全局最优断层面参数及其误差,同时还给出了断层面的顶点坐标,更直观地展示断层形态。该方法被广泛地应用于断层面参数的确定(Zhou等,2010李迎秋等,2011刘白云等,2012杨超群等,2013盛书中等,2014Wang等,2014潘睿等,2015)。为了方便结果的对比分析,断层面参数拟合所选取的各段数据范围(图 1中矩形方框)同万永革等(2008)的研究,即将唐山地震序列分为5段拟合。

3 研究结果及分析

对各段断层面参数进行了分段拟合,拟合结果与万永革等(2008)的结果对比见表 1。其中,第1段为宁河断裂段,该段走向为246.6°,倾角为81.8°,走向和倾角的标准差相对较大,倾角与万永革等(2008)的研究结果相差15.8°,主要原因为该段小震的数量少,丛集性较差,因此反演断层参数结果最差。第2段为唐山断裂南段,该段走向为213.4°,倾角为81.9°。第3段为唐山断裂北段,该段走向为231.4°,倾角为89.1°,断层近乎直立(图 3(c)),小震集中分布在断层面4km范围内(图 3(d)),由于该段小震数目最多,反演得到走向和倾角的标准差最小,与万永革等(2008)的研究结果差值也最小。第4、5段分别为卢龙断裂段和滦县断裂段,其中卢龙断裂段的走向和倾角分别为46.1°和89.3°,断层面近乎直立,滦县断裂段的走向和倾角分别为125.1°和76.2°,结果与万永革等(2008)的研究结果相近。由表 1可见,各段断层顶点坐标的反演结果与万永革等(2008)的结果较为一致,其原因可能为两者反演断层面参数时选择的地震资料范围一致;断层面的深度分布与万永革等(2008)的结果相比均较浅,本文给出的断层面上边界均为4km左右,万永革等(2008)给出的断层面上边界分布更深,为6—10km。

表 1 运用小震资料求得的唐山地震序列各段断层面走向、倾角、标准差和位置 Table 1 Fault plane parameters determined by using small earthquake for segments in Tangshan earthquake

图 3 唐山断裂北段小震拟合结果 Fig. 3 The fitting result by using small earthquakes along the north segment of Tangshan fault

前人对唐山地震断层面作了大量研究,根据P波初动(张之立等,1980)和野外地质调查(虢顺民等,1977尤惠川等,2002)得到断层走向为N30°E;根据面波资料(Butler等,1979)得到断层面走向为N40°E;根据大地测量资料(陈运泰等,1979)和卫星资料(王景明等,1981)得到断层走向为N50°E。上述结果表明唐山地震破裂非常复杂,实际断层面并非1个简单的几何面,而是错综复杂的破裂体。陈运泰等(1979)的研究结果表明唐山地震总体走向为N49°E,初始破裂的走向为N30°E,即破裂起始于N30°E走向的南段断裂,而后转向N50°E走向的北段断裂(万永革等,2008),与本文给出的断层南段走向213.4°、北段走向231.4°基本一致。对于断层倾角,李钦祖等(1980)张之立等(1980)给出总体倾向为120°,倾角为90°;Nebelek等(1987)给出断层为西北倾向,南段倾角为78°,北段倾角为80°;尤惠川等(2002)通过野外地质调查给出总体断层倾向向西,倾角70°—80°,与本文给出的断层南段倾角81.9°、北段倾角89.1°较为一致。杜晨晓等(2010)根据前人研究结果得到滦县地震的震源参数,其断层走向为N120°E,倾角为80°,与本文给出的滦县段反演结果较为一致。

本文结果与万永革等(2008)的研究结果存在一定差异,为判定每段差异是否在可接受的范围内,搜集了18组由不同学者和机构给出的同一地震断层面参数(表 2),对其差异进行统计可在一定程度上反映对断层面参数的约束能力,并将其作为参考,探讨本文研究结果的合理性。各组参数中走向和倾角的最大、最小差异值如图 4所示。

表 2 地震断层面参数研究结果差异统计 Table 2 Statistical results of the differences in seismic fault plane parameters

图 4 同一地震断层走向(a)和倾角(b)差异值统计 Fig. 4 The differences of fault strike (a) and dip (b) of the same earthquake by different studies

图 4可见,不同学者和机构给出的断层走向和倾角的差异值虽有一定的离散性,但其差异范围可作为参考。本文得到的走向和倾角的差异值均在18组结果的最大及最小差异范围内,故利用中国地震台网统一地震目录计算得到的结果是可靠的,证明中国地震台网统一地震目录可以用于断层面参数的确定。

4 结论与讨论

本文各段断层的顶点坐标(表 1)与万永革等(2008)的结果一致,但断层深度有一定的差距。本文结果显示各段断层上边界均为4km左右,万永革等(2008)给出的断层上边界分布更深,造成差异的主要原因可能是万永革等(2008)认为5km以上的小震是深部破裂引发的沉积层破裂,故在确定发震断层面参数时未考虑5km以上的小震,而本文在进行断层面参数反演时,将所有的地震事件都考虑在内,因此反演的断层面上边界分布较浅。

胡新亮等(2001)运用小孔径数字地震台网对唐山地区的地震进行重新定位,通过对比以往地震目录中给出的震源深度,表明唐山地区的地震发生在地壳浅层。于湘伟等(2010)采用双差定位法对华北地区的地震重新定位,精定位后的结果显示83%的地震震源深度位于0—15km,与其他学者的研究结果一致(胡幸平等,2013赵博等,2013李红光等,2015),上述研究表明华北地区的地震主要发生在中上地壳。王椿镛等(2017)给出华北地壳厚度为32—35km,其中上地壳厚度为10—12km,中地壳为8—10km。本文采用的地震数据震源深度主要分布于3—20km,位于华北地壳的中上部,与华北地区地震震源深度分布一致(图 2(a)),由此表明统一地震目录数据给出的地震深度范围是合理的,因此,反演出的断层面上、下边界位置是可靠的。

本文基于中国地震台网统一地震目录提供的地震资料,应用小震确定断层面参数法确定了唐山地震序列的断层面参数,研究中数据分段及其选取范围均参考了万永革等(2008)的研究。对所得的各断层面参数进行对比分析,表明本文结果与万永革等(2008)研究结果的差异在可接受范围内,证明了中国地震台网统一地震目录可以用于断层面参数的确定。

随着地震台网的加密布设,其地震定位能力和定位精度均显著提高。地震活跃地区大量的地震定位数据为拟合断层的几何形态奠定了基础,今后,国家台网统一地震目录可广泛地应用于活断层的发现及其形态的确定。

致谢: 感谢审稿专家提出的宝贵修改意见以及国家地震科学数据共享中心提供的数据。
参考文献
陈运泰, 林邦慧, 王新华, 等, 1979. 用大地测量资料反演的1976年唐山地震的位错模式[J]. 地球物理学报, 22(3): 201-217. DOI:10.3321/j.issn:0001-5733.1979.03.001
杜晨晓, 谢富仁, 张扬, 等, 2010. 1976年MS 7.8唐山地震断层动态破裂及近断层强地面运动特征[J]. 地球物理学报, 53(2): 290-304. DOI:10.3969/j.issn.0001-5733.2010.02.007
郭祥云, 陈学忠, 李艳娥, 2010. 2008年5月12日四川汶川8.0级地震与部分余震的震源机制解[J]. 地震, 30(1): 50-60.
郭志, 高星, 王卫民, 2008. 2006年12月26日台湾南部滨海MS 7.2级地震破裂过程研究[J]. 地球物理学报, 51(4): 1103-1113. DOI:10.3321/j.issn:0001-5733.2008.04.019
虢顺民, 李志义, 程绍平, 等, 1977. 唐山地震区域构造背景和发震模式的讨论[J]. 地质科学, 12(4): 305-320.
胡新亮, 刁桂苓, 高景春, 等, 2001. 华北、西南一些地区地震发生在地壳浅部的证据[J]. 地震学报, 23(4): 427-435. DOI:10.3321/j.issn:0253-3782.2001.04.011
胡幸平, 崔效锋, 2013. 华北地区中部地震精定位与构造应力场研究[J]. 震灾防御技术, 8(4): 351-360. DOI:10.3969/j.issn.1673-5722.2013.04.002
黄文辉, 康英, 苏柱金, 等, 2016. 全国统一编目系统设计与实现[J]. 地震地磁观测与研究, 37(4): 170-175.
李红光, 王利亚, 孙刚, 等, 2015. 华北地区中小地震重新定位和地震活动特征研究[J]. 地震, 35(1): 28-37.
李启雷, 李玉丽, 马丽, 等, 2016. 2016年青海门源MS 6.4地震震源机制与震源深度计算[J]. 地震研究, 39(S1): 55-61.
李钦祖, 张之立, 靳雅敏, 等, 1980. 唐山地震的震源机制[J]. 地震地质, 2(4): 59-67.
李晓峰, 2017. 2016年1月21日青海门源MS 6.4地震震源机制解及发震构造初步探讨[J]. 地震工程学报, 39(4): 657-661. DOI:10.3969/j.issn.1000-0844.2017.04.0657
李迎秋, 刘鑫, 万永革, 2011. 利用智利地震震源分布确定智利地区主震断层面[J]. 防灾科技学院学报, 13(2): 51-56. DOI:10.3969/j.issn.1673-8047.2011.02.010
刘白云, 袁道阳, 张波, 等, 2012. 1879年武都南8级大地震断层面参数和滑动性质的厘定[J]. 地震地质, 34(3): 415-424. DOI:10.3969/j.issn.0253-4967.2012.03.003
刘丽芳, 徐甫坤, 2014. 利用精确定位余震资料确定2014年云南鲁甸6.5级地震的断层面参数[J]. 地震研究, 37(4): 489-494. DOI:10.3969/j.issn.1000-0666.2014.04.001
龙思胜, 2004. 2001年四川雅江6.0级地震序列的破裂特征及发震构造[J]. 中国地震, 20(1): 1-11. DOI:10.3969/j.issn.1001-4683.2004.01.001
罗钧, 赵翠萍, 周连庆, 2015. 2013年8月香格里拉德钦-得荣MS 5.9地震序列震源机制与应力场特征[J]. 地球物理学报, 58(2): 424-435.
吕坚, 王晓山, 苏金蓉, 等, 2013. 芦山7.0级地震序列的震源位置与震源机制解特征[J]. 地球物理学报, 56(5): 1753-1763.
潘睿, 盛书中, 万永革, 2015. 2009年姚安地震断层面参数的确定[J]. 华北地震科学, 33(2): 63-68. DOI:10.3969/j.issn.1003-1375.2015.02.011
盛书中, 万永革, 王未来, 等, 2014. 2010年玉树MS 7.1地震发震断层面参数的确定[J]. 地球物理学进展, 29(4): 1555-1562.
苏金蓉, 郭志, 2013. 西南天山-帕米尔2008年10月5日MS 6.7级地震震源机制研究[J]. 地球物理学报, 56(2): 504-512.
万永革, 沈正康, 刁桂苓, 等, 2008. 利用小震分布和区域应力场确定大震断层面参数方法及其在唐山地震序列中的应用[J]. 地球物理学报, 51(3): 793-804. DOI:10.3321/j.issn:0001-5733.2008.03.020
王椿镛, 吴庆举, 段永红, 等, 2017. 华北地壳上地幔结构及其大地震深部构造成因[J]. 中国科学:地球科学, 47(6): 684-719.
王福昌, 靳志同, 钱小仕, 等, 2012. 由余震分布确定大地震子断层及其参数的模糊聚类方法[J]. 地震学报, 34(6): 793-803. DOI:10.3969/j.issn.0253-3782.2012.06.006
王福昌, 万永革, 钱小仕, 等, 2013. 由地震分布丛集性给出断层参数的一种新方法[J]. 地球物理学报, 56(2): 522-530.
王景明, 郑文俊, 陈国顺, 等, 1981. 唐山地震地面主破裂带及地震成因探讨[J]. 地震研究, 4(4): 437-450.
王鸣, 王培德, 1992. 1989年10月18日大同-阳高地震的震源机制和发震构造[J]. 地震学报, 14(4): 407-415.
吴宝峰, 2017. 2017年11月18日西藏米林6.9级地震震源机制[J]. 地震地磁观测与研究, 38(6): 26-29. DOI:10.3969/j.issn.1003-3246.2017.06.005
谢祖军, 金笔凯, 郑勇, 等, 2013. 近远震波形反演2013年芦山地震震源参数[J]. 中国科学:地球科学, 43(6): 1010-1019.
谢祖军, 郑勇, 姚华建, 等, 2018. 2017年九寨沟MS 7.0地震震源性质及发震构造初步分析[J]. 中国科学:地球科学, 48(1): 79-92.
杨超群, 孟凡顺, 万永革, 2013. 采用精确定位小震资料反演伽师地震断层面[J]. 地球物理学进展, 28(6): 2865-2871.
杨宜海, 范军, 花茜, 等, 2017. 近震全波形反演2017年九寨沟M 7.0地震序列震源机制解[J]. 地球物理学报, 60(10): 4098-4104. DOI:10.6038/cjg20171034
易桂喜, 龙锋, 梁明剑, 等, 2017. 2017年8月8日九寨沟M 7.0地震及余震震源机制解与发震构造分析[J]. 地球物理学报, 60(10): 4083-4097. DOI:10.6038/cjg20171033
尤惠川, 徐锡伟, 吴建平, 等, 2002. 唐山地震深浅构造关系研究[J]. 地震地质, 24(4): 571-582. DOI:10.3969/j.issn.0253-4967.2002.04.012
于湘伟, 张怀, 陈运泰, 2010. 华北地区地震重新定位结果分析[J]. 大地测量与地球动力学, 30(2): 29-33.
张宏志, 刁桂苓, 陈祺福, 等, 2008. 1976年唐山7.8级地震震区现今地震震源机制分析[J]. 地震研究, 31(1): 1-6.
张辉, 王熠熙, 2012. 2012年5月3日金塔-阿拉善盟5.4级地震震源机制解[J]. 西北地震学报, 34(2): 205-206.
张之立, 李钦祖, 谷继成, 等, 1980. 唐山地震的破裂过程及其力学分析[J]. 地震学报, 2(2): 111-129.
赵博, 高原, 石玉涛, 2013. 用双差定位结果分析华北地区的地震活动[J]. 地震, 33(1): 12-21. DOI:10.3969/j.issn.1000-3274.2013.01.002
赵翠萍, 陈章立, 郑斯华, 等, 2005. 2003年9月27日中、俄、蒙边界MS 7.9地震震源机制及破裂过程研究[J]. 地震学报, 27(3): 237-249. DOI:10.3321/j.issn:0253-3782.2005.03.001
赵旭, 黄志斌, 房立华, 2014. 2014年云南盈江MS 6.1地震震源机制研究[J]. 中国地震, 30(3): 462-473. DOI:10.3969/j.issn.1001-4683.2014.03.019
郑勇, 马宏生, 吕坚, 等, 2009. 汶川地震强余震MS(≥ 5.6)的震源机制解及其与发震构造的关系[J]. 中国科学D辑:地球科学, 39(4): 413-426.
Nábělek J., Chen W. P., Ye H., 1987. The Tangshan earthquake sequence and its implications for the evolution of the North China basin[J]. Journal of Geophysical Research, 92(B12): 12615-12628. DOI:10.1029/JB092iB12p12615
Ouillon G., Ducorbier C., Sornette D., 2008. Automatic reconstruction of fault networks from seismicity catalogs:three-dimensional optimal anisotropic dynamic clustering[J]. Journal of Geophysical Research, 113(B1): B01306.
Wang X. S., Feng X. D., Xu X. W., et al, 2014. Fault plane parameters of Sanhe-Pinggu M 8 earthquake in 1679 determined using present-day small earthquakes[J]. EarthqButler R., Stewart G. S., Kanamori H.. 1979. The July 27. 1976 Tangshan, China earthquake-a complex sequence of intraplate events. Bulletin of the Seismological Society of America, 69(1): 207-220.
Nábělek J., Chen W. P., Ye H., 1987. The Tangshan earthquake sequence and its implications for the evolution of the North China basin[J]. Journal of Geophysical Research, 92(B12): 12615-12628. DOI:10.1029/JB092iB12p12615
Ouillon G., Ducorbier C., Sornette D., 2008. Automatic reconstruction of fault networks from seismicity catalogs:three-dimensional optimal anisotropic dynamic clustering[J]. Journal of Geophysical Research, 113(B1): B01306.
Wang X. S., Feng X. D., Xu X. W., et al, 2014. Fault plane parameters of Sanhe-Pinggu M 8 earthquake in 1679 determined using present-day small earthquakes[J]. Earthquake Science, 27(6): 607-614. DOI:10.1007/s11589-014-0099-3
Zhou C. Y., Diao G. L., Geng J., et al, 2010. Fault plane parameters of Tancheng M 81/2 earthquake on the basis of present-day seismological data[J]. Earthquake Science, 23(6): 567-576. DOI:10.1007/s11589-010-0756-0