2. 中国地震局地震预测研究所, 北京市复兴路63号, 100036
大量研究证实[1-5],地壳内部介质普遍存在各向异性。EDA理论[6]认为,各向异性的来源主要是地壳内存在大量充满流体、呈近垂直平行排列的微裂隙,在偏应力作用下,这些裂隙不均匀地成行排列,从而表现出地震各向异性。当剪切波在定向排列的微裂隙中传播时,会分裂成分别以不同速度前进、波振动方向近乎垂直的两列波[7-18],依传播的快慢分为快波和慢波,其中,快波的偏振方向跟应力方向、活动断层走向有关,慢波的延迟时间则跟地壳介质的裂隙密度和各向异性程度有关[19-21]。改变应力场会改变微裂隙的几何参数,造成剪切波分裂参数的变化[22-28]。因此,追踪研究剪切波分裂参数的时空信息变化,对探索地震孕育的过程和发震机理具有应用前景[1, 29-33]。
2013-10-01乳山发生ML3.8地震之后余震不断,形成震群活动,目前仍在持续。截止到2016-03-01,共发生ML<0.0地震3 640次,0.0≤ML≤0.9地震7 108次,1.0≤ML≤1.9地震1 444次,2.0≤ML≤2.9地震228次,3.0≤ML≤3.9地震26次,4.0≤ML≤4.9地震3次,ML≥5.0地震1次(图 1给出了乳山地震序列的M-t图),并于2014-01-07、04-04、09-16和2015-05-22在乳山又先后发生ML4.7、ML4.5、ML4.1和ML5.0地震,造成较大的社会影响。基于当前地震活动的特点,为了更好地对这个震群进行详细研究,从2014-05起,山东省地震局在震中附近十几km范围内架设18个流动地震台,组成临时台阵。这些地震台站记录了大量近场地震数字记录,为研究乳山地震近场区的地震学特征提供了丰富的数据资料。
本文利用乳山地区的乳山台(RSH)和18个流动地震台组成的台阵记录到的波形记录,分析乳山地震近场区的各向异性参数, 研究该观测区域的应力环境特征,为研究乳山地震近场区的地震学特征及乳山震群的发震机理提供参考依据。
1 资料和方法截至2016-03-01,山东数字地震台网共记录到12 450个乳山地震事件。由于地震数量太多,本文随机抽取150次地震用于RSH台剪切波分裂分析。所抽取的样本满足以下条件:1)P波与S波的到时差ΔtS-P≤1.7 s; 2)-0.9≤ML≤1.0。根据山东数字地震台网和乳山台阵所定位的地震,其震源深度为6~14 km。流动台剪切波分裂使用的资料为2014-06~2015-12乳山台阵记录到的地震事件,选取剪切波窗口内共90条地震记录用于计算分析。
对18个流动台记录的波形数据开展分析,发现一些流动台仪器参数配置可能有错误,只有DSDC、XHJC、PUSC、BSTZ、HYSZ、CXZX和GJDT 7个流动台记录的数据可以用于剪波分裂分析,故本文采用SAM方法[34-35]对这7个流动台及乳山台(RSH)记录到的240条波形记录进行处理,台站分布如图 2所示。该方法包括相关函数计算、时间延迟校正和偏振分析检验3个部分。
下面结合具体震例来阐述数据处理流程。图 3为乳山台记录到的一个地震事件波形,震级为ML0.0。对地震波形数据进行滤波、相关函数计算、时间延迟校正,如果经过时间延迟校正后接近线性偏振,说明计算结果可靠。图 4和图 5分别给出了时间延迟校正前后的质点偏振图及剪切波波形。
从图 4、图 5看出,经过偏振方向及时间延迟校正后,剪切波质点呈线性偏振,说明所求的剪切波分裂参数是正确可靠的。
2 结果与讨论乳山地区位于胶东地块东南部,南邻黄海。从全球构造来看,其位于环太平洋构造带上,受西太平洋板块俯冲作用和郯庐断裂左行走滑运动的影响。区域内断裂构造十分发育,按发育方向主要有NNE、NE、NW和NS方向,NNE和NE方向断裂最为发育[36-37]。
一些学者曾对胶东地区及山东地区背景应力场作过研究。张玲等[38]对胶东地区地震震源机制解进行计算分析,得出胶东地区背景应力场方向为ENE;山长仑等[39]利用P、SH、SV波的初动及振幅比获得山东及附近区域132次地震的震源机制解,对该区域地壳应场进行分析,得出山东及附近区域背景应力场的方向为ENE;郑建常等[40]使用Gephart和Forsyth的FMSI方法,利用山东地区387次中小地震震源机制解数据,得出胶东地区最大主压应力场为ENE。
本文通过对乳山地震近场区分布的7个流动地震台和乳山台(RSH)记录到的地震序列波形数据进行分析计算,得到8个台站的剪切波分裂参数,如表 1所示。从这8个台站快剪切波偏振方向的等面积玫瑰图(图 6、图 7)可以看出,所有台站的快剪切波的偏振优势方向比较明显;台站RSH、PUSC、HYSZ、XHJC的快剪切波优势偏振方向为近NE方向,该优势偏振方向与2次显著性乳山地震震源机制解的P轴方向基本一致,这2次显著性地震分别是乳山ML4.7和ML4.5(图 8), 与该地区背景主压应力场的方向总体上一致;台站GJDT位于海岛上,其快剪切波的偏振方向比较离散,可能与不规则的地形和复杂的局部结构有关[29, 41-42]。
台站DSDC的快剪切波偏振优势方向为近NW方向,此方向与乳山地震序列震中的时间演化方向一致[43],如图 9所示。对比乳山地震序列中较大地震的震源机制(图 8)发现,震源机制的NW向节面走向和地震序列展布的方向吻合。据岩石资料及野外地质调查,乳山地区成矿后的断裂十分发育,主要有NW、NNE和NE向断裂,其中NW向断裂形成较晚,切割其他方向的断层,以左行压扭性为主,具多期活动特点[44]。另外,该地震序列附近的乳山断裂,其产状与序列展布的方向和CAP计算结果无法吻合。综合上述分析推测,DSDC台附近存在近NW向的断裂构造,乳山震群可能是该断裂活动的结果,该断裂可能为乳山地震序列的发震断层。
乳山地区基底构造为栖霞复背斜东段,中生代以来, 随着强烈的岩浆活动, 东西向构造基底被破坏, 在其上叠加了一系列大规模的NNE、NE向断裂构造[44-45],因此,该地区存在EW向断裂构造的基底。台站BSTZ和CXCZ的快剪切波偏振优势方向为近EW方向,与地区背景应力场的方向不同,认为附近可能存在近EW向的断裂构造。
乳山地区分布着许多规模不等、方向不一的断裂构造,这样复杂的地质构造背景造成该地区构造应力非常复杂,其构造应力场受不同的地质构造影响,因此乳山地区地壳介质地震各向异性受到区域构造应力和局部构造活动的双重影响。
3 结语本文利用乳山台阵和乳山台(RSH)记录到的240条波形,通过剪切波分裂分析,计算得出乳山地震近场区的各向异性参数,并结合震源机制和地震序列展布方向,推测出乳山地震序列的发震断层,获得该区域的应力环境构造特征。
本文研究结果显示,乳山地区地壳介质地震各向异性受到区域构造应力和局部构造活动的双重影响,台站RSH、PUSC、HYSZ、XHJC的快剪切波优势偏振方向为近NE方向,与该地区背景主压应力场方向总体一致,与该研究区内的2次显著性乳山地震的震源机制P轴的方向基本一致;GJDT台的快剪切波的优势偏振方向比较离散,这与该台站位于海岛上有关,也可能与不规则的地形和复杂的局部结构有关;台站DSDC的快剪切波偏振优势方向与该地区背景主压应力场的方向不一致,为近NW方向,该方向与乳山地震序列震中的时间演化方向一致,与CAP计算结果NW节面走向一致,据此推测DSDC台附近隐伏近NW方向的活动断裂,地震序列可能是该断裂的活动结果,该断裂可能为乳山地震序列的发震断层;BSTZ台和CXZX台的快剪切波偏振优势方向为近EW方向,认为BSTZ台和CXZX台附近可能存在近EW向的活动断裂。
[1] |
Crampin S, Chastin S, Gao Y. Shear-Wave Splitting in a Critical Crust:Ⅲ-Preliminary Report of Multi-Variable Measurements in Active Tectonics[J]. J Appl Geophys, 2003, 54: 265-277 DOI:10.1016/j.jappgeo.2003.01.001
(0) |
[2] |
Silver P G. Seismic Anisotropy Beneath the Continents:Probing the Depths of Geology[J]. Annual Rev Earth Planet Sci, 1996, 24: 385-432 DOI:10.1146/annurev.earth.24.1.385
(0) |
[3] |
Woodhouse J H, Giardini D, Li X D. Evidence for Inner-Core Anisotropy from Free Oscillations[J]. Geophys Res Lett, 1986, 13: 1549-1552 DOI:10.1029/GL013i013p01549
(0) |
[4] |
Morelli A, Dziewonski A M, Woodhouse J H. Anisotropy of the Inner-Core Inferred from PKIKP Travel Times[J]. Geophys Res Lett, 1986, 13: 1545-1548 DOI:10.1029/GL013i013p01545
(0) |
[5] |
Creager K C. Anisotropy of the Inner Core from Differential Travel Times of the Phases PKP and PKIKP[J]. Nature, 1992, 356: 309-314 DOI:10.1038/356309a0
(0) |
[6] |
Crampin S. Seismic Wave Propagation Through a Cracked Solid:Polarization as A Possible Dilatancy Diagnostic[J]. Geophys J R Astron Soc, 1978, 53: 467-496 DOI:10.1111/gji.1978.53.issue-3
(0) |
[7] |
Crampin S, Evans R. Analysis of Records of Local Earthquakes:The Turkish Dilatancy Projects (TDP1 and TDP2)[J]. Geophys J R Astron Soc, 1985, 83: 1-16 DOI:10.1111/j.1365-246X.1985.tb05152.x
(0) |
[8] |
Crampin S, Booth D C. Shear-Wave Polarizations Near the North Anatolian Fault-Ⅱ.Interpretation in Terms of Crack-Induced Anisotropy[J]. Geophys J R Astron Soc, 1985, 83: 75-92 DOI:10.1111/j.1365-246X.1985.tb05157.x
(0) |
[9] |
Crampin S, Booth D C, Krasnova M A, et al. Shear-Wave Polarizations in the Peter the First Range Indicating Crack-Induced Anisotropy in a Thrust-Fault Regime[J]. Geophys J R Astron Soc, 1986, 84: 401-412 DOI:10.1111/j.1365-246X.1986.tb04362.x
(0) |
[10] |
Crampin S, Lovell J H. A Decade of Shear-Wave Splitting in the Earth's Crust:What Use Can We Make of It?and What Should Use Do Next[J]. Geophys J Int, 1991, 107: 387-408 DOI:10.1111/gji.1991.107.issue-3
(0) |
[11] |
Evans R, Beamish D, Crampin S, et al. The Turkish Dilatancy Project (TDP3):Multidisciplinary Studies of A Potential Earthquake Source Region[J]. Geophys J R Astron Soc, 1987, 91: 265-286 DOI:10.1111/j.1365-246X.1987.tb05227.x
(0) |
[12] |
Chen T C, Booth D C, Crampin S. Shear-Wave Polarizations Near the North Anatolian Fault-Ⅲ.Observations of Temporal Changes[J]. Geophys J R Astron Soc, 1987, 91: 287-311 DOI:10.1111/j.1365-246X.1987.tb05228.x
(0) |
[13] |
Daley T M, McEvilly T V, Majer E L. Multiply-Polarized Shear-Wave VSPs from the Cajon Pass Drillhole[J]. Geophys Res Lett, 1988, 15: 1001-1004 DOI:10.1029/GL015i009p01001
(0) |
[14] |
Kaneshima S, Ito H, Sugihara M. Shear Wave Polarization Anisotropy Observed in a Rift Zone in Japan[J]. Tectonophysics, 1989, 157: 281-300 DOI:10.1016/0040-1951(89)90145-5
(0) |
[15] |
Buchbinder G G R. Shear-Wave Splitting and Anisotropy in the Charlevoix Seismic Zone, Quebec[J]. Geophys Res Lett, 1985, 12(7): 425-428 DOI:10.1029/GL012i007p00425
(0) |
[16] |
高原, 冯德益, 郑治真, 等. 大同-阳高Ms5.8级地震前S波分裂异常变化[J]. 华北地震科学, 1993, 11(2): 1-13 (Gao Yuan, Feng Deyi, Zheng Zhizhen, et al. Abnormal Variations of S-Wave Splitting before Datong-Yanggao Ms5.8 Earthquake[J]. North China Earthquake Sciences, 1993, 11(2): 1-13)
(0) |
[17] |
高原, 郑斯华, 孙勇. 唐山地区地壳裂隙各向异性[J]. 地震学报, 1995, 17(3): 283-293 (Gao Yuan, Zheng Sihua, Sun Yong. Seismic Anisotropy of the Crust in Tangshan Area[J]. Acta Seismologica Sinica, 1995, 17(3): 283-293)
(0) |
[18] |
姚陈, 王培德, 陈运泰. 卢龙地区S波偏振与上地壳裂隙各向异性[J]. 地球物理学报, 1992, 35(3): 305-315 (Yao Chen, Wang Peide, Chen Yuntai. Shear-Wave Polarization and Crack Induced Anisotropy of Upper Crust in Lulong, North China[J]. Acta Geophysica Sinica, 1992, 35(3): 305-315 DOI:10.3321/j.issn:0001-5733.1992.03.005)
(0) |
[19] |
Crampin S. Calculable Fluid-Rock Inter Actions[J]. J Geol Soc, 1999, 156: 501-514 DOI:10.1144/gsjgs.156.3.0501
(0) |
[20] |
Gao Y, Wu J, Fu K Y, et al. Shear-Wave Splitting in the Crust in North China:Stress, Faults and Tectonic Implications[J]. Geophys J Int, 2011, 187(2): 642-654 DOI:10.1111/gji.2011.187.issue-2
(0) |
[21] |
Shi Y, Gao Y, Tai L, et al. The Shear-Wave Splitting in the Crust and the Upper Mantle Around the Bohai Sea, North China[J]. Journal of Asian Earth Sciences, 2015, 111: 505-516 DOI:10.1016/j.jseaes.2015.06.015
(0) |
[22] |
Crampin S, Zatsepin S V. Modelling the Compliance of Crustal Rock-Ⅱ.Response to Temporal Changes before Earthquakes[J]. Geophys J Int, 1997, 129: 495-506 DOI:10.1111/gji.1997.129.issue-3
(0) |
[23] |
Gao Y, Crampin S. Observations of Stress Relaxation before Earthquakes[J]. Geophys J Int, 2004, 157: 578-582 DOI:10.1111/gji.2004.157.issue-2
(0) |
[24] |
Zhang Z J, Teng J W, Badal J, et al. Construction of Regional and Local Seismic Anisotropic Structures from Wide-Angle Seismic Data:Crustal Deformation in the Southeast of China[J]. J Seismol, 2009, 13: 241-252 DOI:10.1007/s10950-008-9124-0
(0) |
[25] |
高原, 郑斯华, 王培德. 海南省东方地区1992年小震群剪切波分裂研究[J]. 地球物理学报, 1996, 39(2): 221-232 (Gao Yuan, Zheng Sihua, Wang Peide. Shear Wave Splitting Study on Small Earthquake Swarm of 1992 in Dongfang of Hainan, South China[J]. Acta Geophysica Sinica, 1996, 39(2): 221-232 DOI:10.3321/j.issn:0001-5733.1996.02.009)
(0) |
[26] |
石玉涛, 高原, 吴晶, 等. 云南地区地壳介质各向异性:快剪切波偏振特性[J]. 地震学报, 2006, 28(6): 574-585 (Shi Yutao, Gao Yuan, Wu Jing, et al. Seismic Anisotropy of the Crust in Yunnan, China:Polarizations of Fast Shear-Waves[J]. Acta Seismologica Sinica, 2006, 28(6): 574-585 DOI:10.3321/j.issn:0253-3782.2006.06.002)
(0) |
[27] |
吴晶, 高原, 陈运泰, 等. 首都圈西北部地区地壳介质地震各向异性特征初步研究[J]. 地球物理学报, 2007, 50(1): 209-220 (Wu Jing, Gao Yuan, Chen Yuntai, et al. Seismic Anisotropy in the Crust in Northwestern Capital Area of China[J]. Chinese Journal of Geophysics, 2007, 50(1): 209-220 DOI:10.3321/j.issn:0001-5733.2007.01.027)
(0) |
[28] |
吴晶, 高原, 蔡晋安, 等. 华夏地块东南部地壳地震各向异性特征初步研究[J]. 地球物理学报, 2007, 50(6): 1748-1756 (Wu Jing, Gao Yuan, Cai Jinan, et al. Preliminary Study on Seismic Anisotropy in the Crust in Southeast of Cathaysia Block[J]. Chinese Journal of Geophysics, 2007, 50(6): 1748-1756 DOI:10.3321/j.issn:0001-5733.2007.06.015)
(0) |
[29] |
赖院根, 刘启元, 陈九辉, 等. 首都圈地区横波分裂与地壳应力场特征[J]. 地球物理学报, 2006, 49(1): 189-196 (Lai Yuangen, Liu Qiyuan, Chen Jiuhui, et al. Shear Wave Splitting and the Features of the Crustal Stress Field in the Capital Circle[J]. Chinese Journal of Geophysics, 2006, 49(1): 189-196 DOI:10.3321/j.issn:0001-5733.2006.01.025)
(0) |
[30] |
Crampin S, Volti T, Stefansson R. A Successfully Stress-Forecast Earthquake[J]. Geophys J Int, 1999, 138: 1-5 DOI:10.1046/j.1365-246x.1999.00807.x
(0) |
[31] |
Crampin S, Peacock S. A Review of the Current Under Standing of Shear-Wave Splitting and Common Fallaciesin Interpretation[J]. Wave Motion, 2008, 45: 675-722 DOI:10.1016/j.wavemoti.2008.01.003
(0) |
[32] |
苗庆杰, 刘希强. 山东乳山震群剪切波分裂参数的时间演化特征分析[J]. 地震学报, 2016, 38(2): 220-231 (Miao Qingjie, Liu Xiqiang. Temporal Variation Characteristics of Shear-Wave Splitting for the Rushan Earthquake Swarm of Shandong Province[J]. Acta Seismologica Sinica, 2016, 38(2): 220-231)
(0) |
[33] |
Crampin S. Developing Stress-Monitoring Sites Using Cross-Hole Seismology to Stress-Forecast the Times and Magnitudes of Future Earthquake[J]. Tectonophysics, 2001, 338: 233-245 DOI:10.1016/S0040-1951(01)00079-8
(0) |
[34] |
高原, 刘希强, 梁维, 等. 剪切波分裂系统分析方法(SAM)软件系统[J]. 中国地震, 2004, 20(1): 101-107 (Gao Yuan, Liu Xiqiang, Liang Wei, et al. Systematic Analysis Method of Shear-Wave Splitting:SAM Software System[J]. Earthquake Research in China, 2004, 20(1): 101-107 DOI:10.3969/j.issn.1001-4683.2004.01.012)
(0) |
[35] |
高原, 石玉涛, 梁维, 等. 剪切波分裂系统分析系统SAM (2007):软件系统[J]. 中国地震, 2008, 24(4): 345-353 (Gao Yuan, Shi Yutao, Liang Wei, et al. Systematic Analysis Method of Shear-Wave Splitting SAM (2007):Software System[J]. Earthquake Research in China, 2008, 24(4): 345-353 DOI:10.3969/j.issn.1001-4683.2008.04.004)
(0) |
[36] |
杨峰华, 张均. 试论胶东地区金成矿体系的统一[J]. 地质找矿论丛, 2002, 17(4): 225-233 (Yang Fenghua, Zhang Jun. The Unification of Gold Deposit-Forming System in Jiaodong Region[J]. Contributions to Geology and Mineral resources Research, 2002, 17(4): 225-233 DOI:10.3969/j.issn.1001-1412.2002.04.002)
(0) |
[37] |
高建伟, 赵国春, 毛小红, 等. 山东乳山金青顶金矿成矿构造和应力场研究[J]. 现代地质, 2011, 25(6): 1099-1107 (Gao Jianwei, Zhao Guochun, Mao Xiaohong, et al. Research of Metallotectonics and Tectonic Stress Field of Jinqingding Gold Deposit, Shandong Province[J]. Geoscience, 2011, 25(6): 1099-1107 DOI:10.3969/j.issn.1000-8527.2011.06.008)
(0) |
[38] |
张玲, 周翠英, 王峰吉, 等. 山东地区各分区地震应力场特征分析[J]. 华北地震科学, 2004, 22(4): 12-15 (Zhang Ling, Zhou Cuiying, Wang Fengji, et al. Characteristics of Stress Field in Each Subregion of Shandong Area[J]. North China Earthquake Sciences, 2004, 22(4): 12-15 DOI:10.3969/j.issn.1003-1375.2004.04.004)
(0) |
[39] |
山长仑, 李霞, 范培乐, 等. 山东及附近区域地震断层错动性质与地壳应力场特征[J]. 地震地磁观测与研究, 2013, 34(5/6): 32-39 (Shan Changlun, Li Xia, Fan Peile, et al. The Properties of Earthquake Fault Slip and Features of Crustal Stress Field in Shandong and Nearby Regions[J]. Seismological and Geomagnetic Observation and Research, 2013, 34(5/6): 32-39)
(0) |
[40] |
郑建常, 王鹏, 李冬梅, 等. 使用小震震源机制解研究山东地区背景应力场[J]. 地震学报, 2013, 35(6): 773-784 (Zheng Jianchang, Wang Peng, Li Dongmei, et al. Tectonic Stress Field in Shandong Region Inferred from Small Earthquake Focal Mechanism Solutions[J]. Acta Seismologica Sinica, 2013, 35(6): 773-784 DOI:10.3969/j.issn.0253-3782.2013.06.001)
(0) |
[41] |
孙进, 石玉涛, 赵博, 等. 张家口-渤海地震活动带中东段地壳剪切波分裂研究[J]. 地震, 2013, 33(3): 24-32 (Sun Jin, Shi Yutao, Zhao Bo, et al. On Shear Wave Splitting in the Zhangjiakou-Bohai Seismic Belt[J]. Earthquake, 2013, 33(3): 24-32)
(0) |
[42] |
苗庆杰, 刘希强, 石玉燕, 等. 基于地壳介质各向异性分析山东地区构造应力特征[J]. 中国地震, 2016, 32(1): 63-71 (Miao Qingjie, Liu Xiqiang, Shi Yuyan, et al. Tectonic Stress Analysis based on the Crustal Seismic Anisotropy in the Shandong Area[J]. Earthquake Research in China, 2016, 32(1): 63-71 DOI:10.3969/j.issn.1001-4683.2016.01.006)
(0) |
[43] |
高太忠, 赵伦山, 杨敏之. 山东牟平乳金矿带成矿演化机理探讨[J]. 大地构造与成矿学, 2001, 25(2): 155-160 (Gao Taizhong, Zhao Lunshan, Yang Minzhi. Gold Mineralization and Its Evolution in The Mouping-Rushan Gold Ore Belt, Shandong Province, China[J]. Geotectonica et Metallogenia, 2001, 25(2): 155-160 DOI:10.3969/j.issn.1001-1552.2001.02.007)
(0) |
2. Institute of Earthquake Science, CEA, 63 Fuxing Road, Beijing 100036, China