摘要
2014年2月12日,在新疆于田县发生了里氏7.3级地震.在该地震震中附近,前人研究证明发育了大量规模不同的活动断层(如康西瓦断裂与贡嘎错断裂等).根据地震触发理论,地震发生后因地壳同震变形会导致其周边不同性质断裂破裂应力发生变化,进而影响其地震的潜在危险性.本文利用地震远场波形记录,反演了该地震滑动模型.之后,根据弹性无限半空间位错理论,计算了该地震在近场范围内活动断裂上的同震应力变化.其目的在于讨论于田地震引起的附近断裂上的库仑应力变化以及这些活动断裂可能潜在的地震危险性.在地震发生后,从国际地震学联合会(IRIS)地震数据中心,下载了震中距离介于30°~90°的地震远场波形记录,为保证台站方位角分布均匀,从中挑选了27个不同方位角的高信噪比地震记录参与理论地震图的生成和波形反演过程.我们采用广义射线理论计算生成远场理论地震波形数据.每个子断层参数的反演则利用基于全局化反演的快速模拟退火反演方法.在有限断层反演过程中,我们采用了强调波形拟合的相关误差函数作为待反演的目标函数,拟合的断层参数使目标函数为最小.之后,根据弹性无限半空间位错理论,以库仑破裂准则为基础,结合反演得到的地震震源机制解和地震位错模型,计算该地震引起的近场断层面上库仑应力的变化.由远场波形计算结果可以看到,于田地震的震源深度为10km,地震断层的倾角约71.9°,破裂面上最大的同震位移达到210cm,以左旋走滑为主并具有正倾滑分量,地震能量主要在前15s内释放.由此得到该地震的地震矩为2.91×1019 N·m,地震震级为Mw6.9.于田地震引发的余震,大致分布在三个区域内:普鲁断裂北部、康西瓦断裂东部和贡嘎错断裂中部.弹性应力计算结果表明,于田地震导致阿尔金断裂西段、普鲁断裂中段、康西瓦断裂东段和贡嘎错断裂中段的静态库仑应力明显增加,其中以康西瓦断裂东段和贡嘎错断裂中段应力增量为最大,分别达到了0.05 MPa和0.04 MPa.大量研究证明,当地震所导致的库仑应力变化大于0.01 MPa时将具有明显的地震触发作用.根据本文结果,2014年于田Mw6.9地震使普鲁断裂、贡嘎错断裂和康西瓦断裂上的库仑应力增量均超过了触发阈值,具有被触发出地震的潜在危险.因此,在以后的地震学研究中,应加强对该三条断裂地震危险性的研究和监测.此外,近6年以来,研究区域发生了3次6级以上的地震.这些地震均沿着贡嘎错断裂,由南西向北东迁移,逐步靠近阿尔金断裂,并且逐渐由正倾滑型地震转变为走滑型地震.阿尔金断裂的走滑速率达到了9mm·a-1,所以,尽管本次地震导致的阿尔金断裂库仑应力增量小于0.01 MPa,阿尔金的地震危险性也应该加强监测.
On February 12 2014,an Ms7.3earthquake hit Yutian,Xinjiang,the epicenter is at the intersection of Karakax fault and Gonggacuo fault.Earthquake triggering theory indicates that Coulomb stress on nearby faults will change because of crustal coseismic slip after earthquake,and it will affect earthquake potential risk.In this paper,we estimate focal mechanism and rupture process with far field seismic wave data,and calculate coseismic stress change on nearby faults around epicenter.The purpose is to discuss Coulomb stress change and seismic potential hazard on these faults caused by Yutian earthquake.After the earthquake,we download seismic far field wave data of which epicenter is between 30°to 90°from IRIS.We select27 high SNR(signal to noise ratio)seismic record to make theoretical seismogram.We use generalized ray theory to get synthetic seismic wave map,and sub fault parameter inverse process is based on simulated annealing algorithms.In the way of changing fault parameter to fit actual and synthetic wave form,optimum solution of every sub fault is found.Based on fault rupture inversion model and infinite elastic half space theory,we calculate Coulomb stress change on nearby faults caused by the earthquake.The inversion results show that source depth of earthquake is 10 km,dip angle is 71.9°,largest coseismic displacement is 210 cm.Seismic moment is 2.91×1019 N·m,the main seismic energy is released in former 10 second.Aftershocks are mainly distributed in three regions:north Pulu fault,east Karakax fault and centre Gonggacuo fault.Stress increased significantly on western segment of Altyn fault,central part of Pulu fault,eastern segment of Karakax fault and central segment of Gonggacuo fault.Among them the largest stress change is 0.05 MPa on Karakax fault and 0.04 MPa on Gonggacuo fault.Former research shows that,Coulomb stress change caused by earthquake larger than 0.01 MPa will dramatically increase seismic risk on faults.In our research,stress change on Pulu fault,Gonggacuo fault and Karakax fault exceed triggering threshold,so all these three faults have seismic risk.In the near future seismic potential hazard on these faults should be closely monitored.In the past 6years,3 moderately strong earthquakes happen in the study area,epicenters migrate from Gonggacuo fault to Altyn fault,from southwest to northeast.Though stress change on Altyn fault is much lower,in consideration of 9mm·a-1 slip rate,seismic risk on Altyn fault should cause enough attention.
引文
Barka A.1999.The 17august 1999Izmit earthquake.Science,285(5435):1858-1859.
Chan C H,Stein R S.2009.Stress evolution following the 1999Chi-Chi,Taiwan,earthquake:consequences for afterslip,relaxation,aftershocks and departures from Omori decay.Geophysical Journal International,177(1):179-192.
Deng J S,Sykes L R.1997.Stress evolution in southern Californiaand triggering of moderate-,small-,and micro-size earthquakes.Journal of Geophysical Research:Solid Earth(1978—2012),102(B11):24411-24435.
Deng Q D,Zhang P Z,Ran Y K,et al.2002.Basic feature of Chinaactive structures.Science in China(D Series)(in Chinese),32(12):1020-1030.
Gonzalez-Huizar H,Velasco A A.2011.Dynamic triggering:Stressmodeling and a case study.Journal of Geophysical Research:Solid Earth(1978—2012),116(B2),doi:10.1029/2009JB007000.
Hao J L,Wang W M,Yao Z X,2011.Source process of the 2011Mw9.0Tohuko Japan earthquake.Science China Earth Science(in Chinese),41(6):745-749,doi:10.1007/s11430-011-4241-y.
Harris R A,Simpson R W.1998.Suppression of large earthquakesby stress shadows:A comparison of Coulomb and rate-and-statefailure.Journal of Geophysical Research:Solid Earth(1978—2012),103(B10):24439-24451.
He J,ChéryJ.2008.Slip rates of the AltynTagh,the Kunlun and theKarakorum faults(Tibet)from 3Dmechanical modeling.Earthand Planetary Science Letters,274(1-2):50-58.
He J,PeltzerG.2010.Poroelastic triggering in the 9-22January 2008Nima-Gaize(Tibet)earthquake sequence.Geology,38(10):907-910.
He J K,VernantP,ChéryJ,et al.2013.Nailing down the slip rateof the AltynTagh fault.Geophysical Research Letters,40(20):5382-5386.
Hong S Y,Shen X H,Shan X J,et al.The calculation and analysisof the co-seismic deformation field of Yutian Ms7.3earthquakebasing on the ascending and descending orbit ASAR data.Remote Sensing for Land&Resources(in Chinese),2010,84(1):98-102.
Jia K,Zhou S,Wang R.2012.Stress interactions within the strongearthquake sequence from 2001to 2010in the Bayankala blockof eastern Tibet.Bulletin of the Seismological Society ofAmerica,102(5):2157-2164.
Jiang Tingchen,et al.Deformation analysis of Yutianearthquakebased on ScanSARinterferometry.Journal of Geodesy andGeodynamics,2012,(6):13-16
Johnson K M,Segall P,Yu S B.2005.A viscoelastic earthquakecycle model for Taiwan.Journal of Geophysical Research:Solid Earth(1978—2012),110(B10),doi:10.1029/2004JB003516.
King G C P,Stein R S,Lin J.1994.Static stress changes and thetriggering of earthquakes.Bulletin of the Seismological Societyof America,84(3):935-953.
Lasserre C,PeltzerG,CrampéF,et al.2005.Coseismic deformationof the 2001 Mw=7.8Kokoxiliearthquake in Tibet,measured bysynthetic aperture radar interferometry.J.Geophys.Res.,110:B12408,doi:10.1029/2004JB003500.
Li H B,Van der Woerd J,Sun Z M,et al.2012.Co-seismic andcumulative offsets of the recent earthquakes along the Karakaxleft-lateral strike-slip fault in western Tibet.Gondwana Research,21(1):64-87.
Lin J,Stein R S.2004.Stress triggering in thrust and subductionearthquakes and stress interaction between the southern SanAndreas and nearby thrust and strike-slip faults.Journal ofGeophysical Research:Solid Earth(1978—2012),109(B2),doi:10.1029/2003JB002607.
Liu C L,Zheng Y,Ge C,et al.2013.Rupture process of the M7.0Lushan earthquake.Science in China(Ser D),43(6):1020-1026.
MériauxA S,Ryerson F J,Tapponnier P,et al.2004.Rapid slipalong the central AltynTagh Fault:Morphochronologic evidencefrom Cherchen He and SulamuTagh.J.Geophys.Res,109(B6):B06401,doi:10.1029/2003JB002558.
McCloskey J,Nalbant S S,Steacy S.2005.Indonesian earthquake:Earthquake risk from co-seismic stress.Nature,434(7031):291-291.
Parsons T,Ji C,Kirby E.2008.Stress changes from the 2008Wenchuan earthquake and increased hazard in the Sichuanbasin.Nature,454(7203):509-510.
Robinson D P,Brough C,Das S.2006.The Mw7.8,2001Kunlunshanearthquake:Extreme rupture speed variabilityand effect of faultgeometry.J.Geophys.Res.,111(B8):B08303,doi:10.1029/2005JB004137.
Savage J C.2000.Viscoelastic-coupling model for the earthquakecycle driven from below.Journal of Geophysical Research:Solid Earth(1978—2012),105(B11):25525-25532.
Segall P,Harris R.1987.Earthquake deformation cycle on the SanAndreas fault near Parkfield,California.Journal of GeophysicalResearch:Solid Earth(1978—2012),92(B10):10511-10525.
Shan B,Xiong X,Zheng Y,et al.Stress changes on major faultscaused by Mw7.9 Wenchuan earthquake,May 12,2008.Sci.China Earth Sci.,2009,52(5):593-601.
Shan B,Xiong X,Zheng Y,et al.Stress changes on major faultscaused by 2013 Lushan earthquake,and its relationship with2008 Wenchuan earthquake.Science China(Earth Sciences),2013,56(7):1169-1176,doi:10.1007/s11430-013-4642-1.
Shen Z K,Wang M,Li Y X,et al.2001.Crustal deformation alongthe AltynTagh fault system,western China,from GPS.Journal of Geophysical Research:Solid Earth(1978—2012),106(B12):30607-30621.
Shimazaki K,Nakata T.1980.Time-predictable recurrence modelfor large earthquakes.Geophysical Research Letters,7(4):279-282.
Stein R S,King G C P,Lin J.1992.Change in failure stress on thesouthern San Andreas fault system caused by the 1992magnitude=7.4 Landers earthquake.Science,258(5086):1328-1332.
Stein R S,King G C P,Lin J.1994.Stress triggering of the 1994M=6.7 Northridge,California,earthquake by its predecessors.Science,265(5177):1432-1435.
Stein R S,Barka A A,Dieterich J H.1997.Progressive failure onthe North Anatolian fault since 1939 by earthquake stresstriggering.Geophysical Journal International,128(3):594-604.
Wan Y G,Shen Z K,Sheng S Z,et al.2010.The mechanicaleffects of the 2008 M7.3 Yutian,Xinjiang earthquake on theneighboring faults and itstectonic origin of normal faultingmechanism.ChineseJournal of Geophysics,53(2):280-289.doi:10.3969/j.issn.000i-5733.2010.02.006.
Wang W M,Zhao L F,Li J,et al.2008.Rupture process of theMs8.0 Wenchuan earthquake of Sichuan,China.Chinese J.Geophys.(in Chinese),51(5):1403-1410.
Yao Z X,Ji C.The inverse problem of finite fault study in timedomain.Chinese J.Geophysics.(in Chinese),1997,40(5):691-701.
Zhang G H,Shan X J,Li W D.The Coulomb failure stress changeassociated with the Ms8.0 Wenchuan earthquake and the riskprediction of its surrounding faults.Seismicity and Geology(inChinese),2008,30(4):935-944.
Zhang G W,Lei J S,Sun C Q.2013.Relocation of the 12February2014Yutian,Xinjiang,mainshock(Ms7.3)and its aftershocksequence.Chinese J.Geophys.(in Chinese),57(3):1012-1020,doi:10.6038/cjg20140330.
Ziv A,Rubin A M.2000.Static stress transfer and earthquaketriggering:No lower threshold in sight?.Journal of GeophysicalResearch:Solid Earth(1978—2012),105(B6):13631-13642.
邓起东,张培震,冉勇康等.2002.中国活动构造基本特征.中国科学(D辑),32(12):1020-1030.
郝金来,王卫民,姚振兴.2011.2011年3月11日日本东北沿海Mw9.0级地震震源破裂过程.中国科学:地球科学,41(6):745-749.
洪顺英,申旭辉,单新建等.2010.基于升降轨ASAR的于田Ms7.3级地震同震形变场信息提取与分析.国土资源遥感,(4):98-102.
蒋廷臣,李陶,刘经南.2012.基于宽幅SAR干涉测量的新疆于田7.3级地震形变分析.大地测量与地球动力学,32(6):13-16.
刘成利,郑勇,葛粲等.2013.2013年芦山7.0级地震的动态破裂过程.中国科学:地球科学,43(6):1020-1026
单斌,熊熊,郑勇等.2009.2008年5月12日Mw7.9汶川地震导致的周边断层应力变化.中国科学D辑:地球科学,39(5):537-545.
单斌,熊熊,郑勇等.2013.2013年芦山地震导致的周边断层应力变化及其与2008年汶川地震的关系.中国科学:地球科学,43(6):1002-1009.
万永革,沈正康,盛书中等.2010.2008年新疆于田7.3级地震对周围断层的影响及其正断层机制的区域构造解释.地球物理学报,53(2):280-289,doi:10.3969/j.issn.0001-5733.2010.02.006.
王卫民,赵连锋,李娟等.2008.四川汶川8.0级地震震源过程.地球物理学报,51(5):1403-1410,doi:10.3321/j.issn:0001-5733.2008.05.013
姚振兴,纪晨.1997.时间域内有限地震断层的反演问题.地球物理学报,40(5):691-701.
张国宏,单新建,李卫东.2008.汶川Ms8.0地震库仑破裂应力变化及断层危险性初步研究.地震地质,30(4):935-944.
张光伟,雷建设,孙长青.2014.2014年2月12日新疆于田Ms7.3级地震主震及余震序列重定位研究.地球物理学报,57(3):1012-1020,doi:10.6038/cjg20140330.