地球物理学报  2016, Vol. 59 Issue (6): 2126-2137   PDF    
2014年新疆于田MS7.3地震序列重定位及其发震构造初步研究
唐明帅 , 王海涛 , 李艳永 , 热依木江 , 孔祥艳 , 魏芸芸 , 上官文明 , 魏斌     
新疆维吾尔自治区地震局, 乌鲁木齐 830011
摘要: 2014年2月12日在新疆于田县发生了MS7.3地震,主震前一天在震区发生了MS5.4前震,震后余震活动频繁,由于震区台站十分稀疏和不均匀、地壳速度结构复杂,台网常规定位结果精度有限,很难从中获得序列的空间分布特征和活动趋势的正确认识.本文首先利用位于震区附近的于田地震台5年记录的远震波形数据,采用接收函数方法研究了震区附近的地壳结构,建立了震源区的地壳速度模型.在此基础上,联合震相到时和方位角对2014年于田MS7.3地震序列(从2014年02月11日-2014年04月30日,共计577次地震)进行了重新绝对定位.结果显示,(1) 重定位后的前震和主震震中位置明显向地表破裂带及其附近的阿尔金分支断裂(南肖尔库勒断裂和阿什库勒-肖尔库勒断裂)靠近,两者相距5.4 km,主震位置为36.076°N、82.576°E,震源深度为22 km, 前震位置为36.055°N、82.522°E,震源深度为19 km;(2) 本文重定位结果显示,余震序列沿NEE-SWW展布,优势分布长度约73 km、宽度约16 km,平均震源深度为14.8 km,其中77%的余震分布在地表破裂带的西南端,这部分余震中少数沿阿什库勒-肖尔库勒断裂分布,绝大多数沿北东东向的南肖尔库勒断裂分布,位于地表破裂带东北端的余震沿阿什库勒-肖尔库勒断裂分布,但发生在地表破裂带的余震极少;重定位后,位于地表破裂带西南侧的震中分布由台网目录的近南北向变为北东向,与地表破裂带、南肖尔库勒断裂和阿什库勒-肖尔库勒断裂走向一致;(3) 沿重定位剖面的地震分布,可推断位于地表破裂带西南段的南肖尔库勒断裂与位于北东段的阿什库勒-肖尔库勒断裂倾向反向,南肖尔库勒断裂的倾向为SE,阿什库勒-肖尔库勒断裂的倾向为NW,这与本次地震野外考察得到的断裂性质一致.综合重定位结果、地表破裂带分布、震源机制解、南肖尔库勒断裂和阿什库勒-肖尔库勒断裂的性质认为,2014年于田MS7.3地震的发震构造为阿尔金断裂西南尾段的两条分支断裂——南肖尔库勒断裂和阿什库勒-肖尔库勒断裂.
关键词: 2014年于田MS7.3地震序列      联合震相到时和方位角绝对定位法      重新定位      发震构造     
Relocation of the 2014 Yutian, Xinjiang, MS7.3 earthquake sequence and a preliminary study of its seismogenic structure
TANG Ming-Shuai, WANG Hai-Tao, LI Yan-Yong, REYI MuJiang, KONG Xiang-Yan, WEI Yun-Yun, SHANGGUAN Wen-Ming, WEI Bin     
Earthquake Administration of Xinjiang Uygur Autonomous Region, Urumqi 830011, China
Abstract: On February 12, 2014, an MS7.3 earthquake occurred in Yutian County, Xinjiang, China. One day before this earthquake there was an MS5.4 foreshock occurred in the same place. After the mainshock, the seismicity of aftershocks was highly active. Study on the characteristics of this earthquake sequence may help us to strengthen our understanding of the seismogenic environment in this area. However, due to the sparse distribution of the seismic network and the complex crustal structure in the source region, it is hard to obtain a precise distribution of the earthquakes to further analyze the seismic sequence spatial distribution and the earthquake occurrence tendency in this source region. Using the teleseismic waveforms data recorded by station Yutian near the 2014 Yutian epicenter region in Xinjiang, we investigated the crustal thickness, Poisson's ratio and the depth of crustal discontinuity beneath station Yutian by grid-search-stacking method of receiver function, further obtained the crustal S-wave velocity structure. Finally, we established the crustal velocity model in the source region for relocating the earthquakes. Based on the crustal velocity model, we relocated the 557 earthquakes which belong to the 2014 Yutian, Xinjiang, MS7.3 earthquake sequence (February 11-April 30, 2014) by the combined seismic phase arrival time and azimuth absolute location algorithm. The following characteristics can be found in the relocated result. (1) The epicenters of the relocated mainshock and the foreshock are much closer to the surface rupture and the Altyn Tagh branch faults (the south Xor Kol fault and the Ashikule-Xor Kol fault) near surface fracture, and the two epicenters are only 5.4 km apart from each other. The location of the mainshock is 82.576°E, 36.076°N and the focal depth is 22 km. The location of the foreshock is 82.522°E, 36.055°N and the focal depth is 19 km. (2) The relocated distribution of the earthquakes are consistent with the surface rupture of the mainshock. The aftershock zone is mainly along NEE-SWW direction with a spatial scale of approximately 73 km long and 16 km wide. Strong variations can be found in the spatial distribution of the aftershocks southwest to the rupture fault, where the aftershock zone changes it strike direction from N-S to northeastward, which is consistent with the strike of the south Xor Kol fault and the Ashikule-Xor Kol fault, which are branches of the Altyn Tagh Fault that is close to the rupture fault. Around 77% aftershocks occurred at the southwest end of the surface rupture zone, mainly along the south Xor Kol fault, and few aftershocks occurred along the surface rupture zone, and the aftershocks at the northeast end of the surface rupture zone mainly occurred on the Ashikule-Xor Kol fault. (3) According to earthquakes distribution along profiles, the dip direction of the south Xor Kol fault (along SE direction) is opposite to that of the Ashikule-Xor Kol fault (along NW direction), which is confirmed by the field surveys. Taking the characters of the south Xor Kol fault and the Ashikule-Xor Kol fault, and the distributions of relocation results, and focal mechanism solution and surface rupture zone into consideration, we speculate that the seismogenic faults of 2014 Yutian MS7.3 earthquake are the south Xor Kol fault and the Ashikule-Xor Kol fault, which are two branches of the Altyn Tagh fault zone..
Key words: 2014 Yutian MS7.3 earthquake sequence      Combined seismic phase and azimuth absolute location algorithm      Relocation      Seismogenic structure     
1 引言

北京时间2014年2月12日下午17时19分,在新疆和田地区于田县发生了MS7.3强烈地震(简称2014年于田7.3地震),主震前一天在震区发生了MS5.4前震,震后余震活动频繁,震中区的烈度为Ⅸ度,等震线长轴呈北东东向分布(新疆维吾尔自治区地震局,2014).此次地震位于2008年3月21日于田MS7.3地震震中东侧约110 km,位于我国巴颜喀拉块体西边界的阿尔金断裂带西南段.

2014年于田MS7.3地震发生后,不同研究单位或研究小组开展了相关的地震学研究,给出了此次地震的震源位置(图 1表 1)、震源机制解(周云等,2015)、破裂过程(张勇等,2014)和地震分布特征(张广伟等,2014王俊等,2014王晓欣等,2014吴传勇等,2014房立华等,2015),这些结果对于理解本次于田地震的发生机理具有重要意义.多家机构给出的震源机制结果认为,本次于田地震显示出左旋走滑性质.震源破裂过程结果表明,此次地震破裂主要集中深部(张勇等,2014).此次于田MS7.3地震震区周边地震台站十分稀疏且基本分布在震中北边,除主震外,只有极少数余震才能被位于震中南边的个别地震台站记录到,图 2为参加地震定位的近台(包括临时台)分布.加之震区地壳结构、活动构造十分复杂,这些因素导致国内外各机构或不同研究小组给出主震和前震的震中平面位置相差较大,分布在阿尔金断裂带西南段的不同分支附近(图 1表 1);图 3为新疆地震台网基于塔吉克3400走时表,利用单纯型定位方法获得此次地震序列的震中和深度分布情况,可以看出此次地震的余震主要集中在主震南部,该部分的余震呈现近南北向分布、跨越了不同的分支断裂,近南北向的分布特征与主震的震源机制解明显不一致.

图 1 不同研究机构或小组给出的主震与前震的震中分布 图中的断层和地表破裂数据李海兵等(2014),虚线为推测断裂,其中f1为南硝尔库勒断裂,f2为阿什库勒—硝尔库勒断裂,f3为北硝尔库勒断裂. Fig. 1 The epicenter distribution of mainshock and foreshock by different research institutions and groups The fault data and rupture data are from Li(2014).Dotted line are inferred faults;f1 is the south Xor Kol fault,f2 is the Ashikule-Xor Kol fault,f3 is the north Xor Kol fault.
表 1 不同研究机构或小组给出的前震和主震定位结果(表 1中带*深度表示矩心深度) Table 1 The results of location of mainshock and foreshock by different research institutions and groups(the H value with asterisk are centroid depth)

地震定位的精度不仅受台网布局、可用定位震相以及射线路径上速度结构等因素影响,还与定位方法密切相关.因此,要获得对此次地震及其序列的准确位置和深入研究,必须选用合理的地壳速度结构和适合台站稀疏地区地震定位的方法.由于位于新疆南部的阿尔金断裂带西南地区自然条件十分恶劣,在2008年以前该地区没有地震台站,目前对该区域的地壳结构研究只有大区域的地震层析成像研究,其分辨率明显不足.远震体波接收函数方法是现阶段国内外运用远震波形数据获得地震台站下方间断面深度的最有效手段(Langston,1979).我们在对2008年于田MS7.3地震及其序列重新定位的基础上认识到,在台网稀疏的情况下,可靠的震相方位角参与地震定位能明显提高地震的定位精度(Bratt and Bache,1988唐明帅等,2010).本文基于震区附近的于田地震台5年记录的远震波形数据,应用接收函数方法研究了台站下方的地壳结构,建立了用于地震定位的地壳速度模型,在此基础上,通过精细分析震相到时和部分台站记录的P波方位角,采用联合震相到时和方位角的定位方法(Schweitzer,2001)对2014年于田MS7.3地震序列进行了重新绝对定位,获得了更加精确的地震定位结果,对余震的空间分布特征和发震构造进行了分析.本文的定位结果对深入研究和正确认识2014年于田MS7.3地震的发震断层及其性质、地表破裂组合特征和地震破裂过程等科学问题提供了更准确的基础数据和约束条件.

2 地震定位资料分析

2014年2月12日于田MS7.3地震震中周围的地震台站分布十分稀疏.当时距震中100 km内的地震台站只有1个于田(YUT)台(约56 km),500 km内的地震台站有8个且都位于震中的北面.主震发生之后,于13日23点和14日18点,新疆地震局在震区北面架设了流动1(LD1)、流动2(LD2)两个流动地震台,两个流动台到主震的距离分别为92 km、82 km(图 2).

图 2 参加地震定位的近台分布 Fig. 2 The distribution of near-seismic stations in study region

根据新疆地震局目录结果,本文选取了震源区ML≥1.6(时间段为2014年2月11日—4月30日、有3组完整的震相或者记录地震台站≥4个)的577次地震进行重新定位,其中447次地震的定位台数介于3~5个.在重定位过程中,根据每次重定位地震的台站分布情况,我们分析了部分台站记录地震波形的方位角,并让这些方位角参与地震定位,增强了对震源信息的约束;同时我们结合当时新疆数字地震台网记录实际波形的情况,对ML≥3.6且定位地震台站大于11个的48次地震做了重新选台.2014年于田主震定位的台站分布,除了应用新疆台网19个台站(YUT、MIF、YJB、QMO、TAZ、HTTZ、SSL、YCH、BCH、AKS、RUQ、KUC、BAC、LTA、 KOL、KSH、KMS、SSH、YMS)记录的数据外,还使用了西藏狮泉河(SQHE)、改则(GZE)、双湖(SUH)地震台,青海沱沱河(TTH)、小黑山(XHS)、乌图美仁(WTM)地震台的数据,这些数据对于改善台网布局具有重要作用,共有25个地震台参与了主震重定位,台站布局最大空隙角为67°.

3 震区附近速度结构研究

在地震定位中,震源区的速度结构是影响地震定位精度的重要因素(Aki and Richards,1980).由于2014年于田MS7.3地震震区附近的地壳结构研究程度很低.因此,本文利用于田地震台从2009年1月—2013年12月记录的远震波形数据,应用接收函数的H-κ叠加方法(Zhu and Kanamori,2000)和线性反演方法(Randall,1989Ammon et al.,1990Herrmann,2004)研究了该台下方的地壳结构.于田台的接收函数波形中,在滞后初至P波约3.0 s处存在一组能量较强且连续的震相,我们推断其为地壳内某界面的一次转换波Ps震相,由该间断面产生的多次转换波PpPs、PsPs+PpSs震相也较明显,图 4A中的红线分别标出了该间断面产生不同震相的理论到时.参考震区附近的人工源地震勘探结果(李秋生等,2000贺日政等,2001),结合区域地质情况,选取6.2 km·s-1和5.5 km·s-1分别作为于田台的地壳和壳内间断面之上的平均P波速度,用于接收函数的H-κ叠加.本文选取位于于田台东南方位的接收函数获得于田台下方的地壳厚度为60.8 km、地壳波速比为1.65(图 4b),于田台下方的壳内间断面深度为24.8 km、地壳波速比为1.64(图 4a).鉴于于田台下方地壳结构的复杂性,以及到2014年于田地震震区存在一定距离(约56 km),为了更好地获得对震源区地壳结构的认识,我们应用和达法(李善邦,1981)计算了于田地震序列中定位台数介于3~5个的447次地震的波速比,获得从震区到台站(YUT、LD1、LD2、YJB、MIF)的平均波速比为1.70.基于上述对于田台地壳结构和震区波速比的研究、结合于田台接收函数的线性反演结果,我们建立了用于2014年于田MS7.3地震序列定位的地壳速度结构如图 4C.

图 3 初始震中分布图和震源深度沿经度方向的剖面 数据来自新疆地震台网的分析结果.圆点表示余震震中(大小与震级成正比),不同圆点的颜色代表不同的震源深度.
红色五角星表示主震和前震的震中位置,色标表示震源深度的范围.断裂数据同图 1.
Fig. 3 The distributions of initial epicenter and focal depth profiles along the longitude The data are from observation reports of Xinjiang Earthquake Administration.Red stars represent the locations of mainshock and foreshock.Circles represent aftershocks with magnitude proportional to its size.The color of the circles changes with depth from blue to red.The fault is same as Fig.1.
图 4 于田台接收函数结果及定位采用的速度模型 (a)于田台下方壳内间断面的H-κ叠加振幅;(b)于田台下方Moho面的H-κ叠加振幅,图中右上角的数字分别表示间断面厚度和平均波速比、字母表示台站名,dH代表间断面厚度误差、dκ代表波速比误差.(A)—(B)为(a)—(b)对应的接收函数剖面图,红线为根据H-κ叠加结果对应间断面产生的转换波(Ps)和多次波(PpPs,PsPs+PpSs)的理论到时;(C)为定位采用的速度模型. Fig. 4 The results of receiver functions at station Yutian and velocity model used in the study (a)The results of inter-crustal discontinuity H-κ grid-stacking-search beneath station YUT;(b)The results of Moho H-κ grid-stacking-search,the numbers in the right above are the discontinuity thickness and wave velocity ratio,the letter is station name,dH is the error of discontinuity thickness,dκ is the error of wave velocity ratio.(A)—(B)are corresponding receiver functions profiles of(a)—(b),the red lines are arrival times of Ps,PpPs,PsPs+PpSs phases for the discontinuity converted waves by the results of H-κ grid-stacking-search.(C)is velocity model used in the study.
4 联合震相到时和方位角的线性定位方法

为了最大限度地让一切可用的数据参与地震定位,Schweitzer开发了HYPOSAT绝对定位方法(Schweitzer,2001),该方法是基于经典的Geiger法(Geiger,19102012),除采用传统的震相(Pn、Sn、Pg、Sg等)到时外,还可应用震相到时差、台站方位角和地震射线参数等来参与定位.其可采用的地球结构速度模型包括使用者给定的具体区域模型和通用的全球模型.Hyposat定位方法基于方程(1)求解的.

(1)

其中ti是走时数据;Δti是走时残差;dtj是同一台站两种震相的走时差;Δdtj是两震相走时差的残差;pk是观测射线参数;Δpk是射线参数的残差;Zb是观测后方位角;ΔZb是观测后方位角的残差;δt0是一次迭代震源时间的改变量;δφ是一次迭代震源纬度改变量;δλ是一次迭代震源经度改变量;δz0是一次迭代震源深度改变量.多次迭代求解方程组,求得震源参数.剔除计算过程中残差过大的数据.

5 定位结果与分析

基于本文建立的速度模型,我们采用HYPOSAT地震定位方法,引入P波震相观测后方位角,设置置信水平为95%,对选取的577次地震进行了重定位.重定位后的误差椭圆平均半长轴为3.27 km、平均半短轴为2.16 km、平均取向为131.5°,走时残差为0.397 s.

5.1 前震、主震定位结果和分析

我们获得2014年于田MS7.3地震及其前震的震源位置信息,并与国内外不同权威机构定位结果进行了比较(表 1图 1).其中前震位置为36.055°N、82.522°E,震源深度为19 km,对应的误差椭圆半长轴为3.77 km,半短轴为2.8 km、取向为172.2°,走时残差为0.890 s.主震位置为36.076°N、82.576°E,震源深度为22 km,对应的误差椭圆半长轴为3.46 km、半短轴为3.32 km、取向为167.9°,走时残差为1.000 s.本文主震的震源深度为22 km,符合地震一般发生在中上地壳的规律.

表 1可知,本文的主震震中纬度与其他研究机构或小组给出的主震震中纬度的最大差为0.171°、最大经度差为0.109°,小于其他不同研究机构或小组给出主震震中之间的最大纬度差0.315°、最大经度差0.133°.从图 1可知,各机构给出的前震和主震位置有较为明显的差异,其震中分别分布在阿尔金断裂带的不同分支附近.本文重定位后的前震和主震震中更靠近地表破裂带及其附近的阿尔金分支断裂(南肖尔库勒断裂和阿什库勒—肖尔库勒断裂),与新疆测震台网测定的结果相比,前震纬度向北移动了约0.113°、主震纬度向南移动了约0.047°,均位于地表破裂带西南端西北,两者的震中相距5.4 km,远小于新疆地震目录给出的前震与主震的震中距离(20.2 km).震后实地考察发现本次地震的破裂带沿阿尔金断裂带西南段的两条近平行的分支断裂南肖尔库勒断裂和阿什库勒—肖尔库勒断裂分布,以左旋走滑性质的破裂为主,这种破裂性质与南肖尔库勒断裂和阿什库勒—肖尔库勒断裂的运动特征一致(李海兵等,20142015),也与前震和主震的震源机制结果一致(周云等,2015李金和王琼,2015).2014年于田MS7.3地震位于阿尔金断裂带尾端多条分支断裂的构造区,本文重定位后前震和主震震中介于南肖尔库勒断裂和阿什库勒—肖尔库勒断裂之间,靠近阿什库勒—肖尔库勒断裂,到南肖尔库勒断裂的垂直距离为5.5 km.综合上述分析,我们推断此次于田MS7.3地震的发震断层为阿尔金断裂西南段的分支断裂——南肖尔库勒断裂和阿什库勒—肖尔库勒断裂.

5.2 重定位的余震震中空间分布特征

精确的震源位置不仅可以为发震构造的判定提供依据,还是研究震源区孕震环境、地震活动性、深部构造以及地球内部物理学等的重要基础数据(张国民等,2002).此外,还可为地震序列的相对定位提供可靠的参考位置(Waldhauser and Ellsworth,2000).震后野外考察表明此次于田地震的地表破裂带沿南肖尔库勒断裂分布长约15 km,沿阿什库勒—肖尔库勒断裂分布长约10 km,两条断裂间不连续分布右阶雁行状张剪裂隙长5 km(李海兵等,20142015).重定位后的余震分布图像(图 5)与重定位前的结果(图 3)相比,重定位后的余震震中分布明显向南肖尔库勒断裂和阿什库勒—肖尔库勒断裂靠近,条带性更加明显,整体沿NEE-SWW走向展布,与地表破裂带、南肖尔库勒断裂和阿什库勒—肖尔库勒断裂的走向一致.重定位后的余震震中分布长度约73 km,但其分布并不连续,沿地表破裂带余震很少,而在地表破裂带之外的西南端和北东端呈两个条带分布(图 5).其中位于地表破裂带北东的余震约占余震总数的23%,余震分布相对均匀,有少数震级大于4.0的余震,震中分布长度约为22 km,起始点紧邻地表破裂带东端,余震总体走向和位置与地表破裂带东段延长线一致.位于地表破裂带西南的余震数量和强度都显著大于地表破裂带西北的余震规模,数量约占余震总数的77%,震中分布长度约为31 km,起始点紧邻破裂带西端,总体走向和位置与地表破裂带西段延长线一致.余震数量在主震北东和西南两侧存在差异,表明两侧断裂特征存在差异.

图 5 重新定位后的震中分布及震源深度剖面位置 AB剖面上的白色三角为剖面的中心位置,LK为沿断层走向的剖面,图中主震震源机制周云等(2015) Fig. 5 The distributions of relation epicenter and the location of depth profiles White triangle on the profile AB indicates the center of the section,LK represent the fault strike.The focal mechanisms for main shock is from Zhou(2015)

重定位前后变化最显著的是位于地表破裂带西南端的余震分布.台网目录给出这部分地震的震中为近南北向展布,跨越了阿尔金断裂带西南段尾端的多条分支断裂(图 3),重新定位后的这部分余震分布更加密集,其震中分布的优势方向呈北东向展布,少量分布在阿什库勒—肖尔库勒断裂两侧,大多数分布在南肖尔库勒断裂两侧,南盘偏多(图 5),符合野外考察(李海兵等,2015)获得的断层性质(走滑兼逆冲),与走滑型地震在断层两侧的余震分布相对均匀、逆冲型地震的余震主要分布在断层上盘(Chang et al.,2000房立华等,2013黄媛等,2008)的结论吻合.位于地表破裂带北东端的余震主要沿阿什库勒—肖尔库勒断裂展布.从主震、余震空间分布可以看出,本次地震破裂具有双向破裂的特点,即从主震分别向北东、西南两侧扩展.根据重定位序列分布求得的地震断层面倾角(万永革等,2008)为88.77°,与阿尔金走滑断裂西南端倾角较陡的情况吻合,对此次地震序列的震源机制特征分析表明,震源机制解给出的节面倾角90°,位于地表破裂带西南端的地震有走滑、正断和逆断型地震(李金和王琼,2015).综合以上分析,我们推断2014年于田MS7.3地震的发震断层为左旋走滑兼有逆断分量的南肖尔库勒断裂和左旋走滑兼有伸展分量的阿什库勒—肖尔库勒断裂.

我们把距断层走向剖面(图 5中的LK剖面)6km范围内的地震定位结果投影到断层面静态滑动量分布图(图 6)上,可以看出在主震破裂区静态滑动量大的区域,余震分布很少且震级较小,余震主要分布在滑动量较小甚至为零的区域,表明主破裂区的大部分应力在主震过程中得以释放,而在主震没有破裂的区域,可能是介质强度相对较大,主震破裂导致的应力增强不足以引起同震破裂,但引起了应力积累,最终以余震的方式释放(王卫民等,2005赵翠萍等,2008冀战波等,2014).

图 6 主震破裂滑动量(来自张勇等,2014)与地震在破裂面上的分布 圆点表示地震,大小与震级成正比,五角星同图 2 Fig. 6 The distribution of static fault slip and earthquake sequence(Zhang et al,2014) Red stars represent the locations of mainshock and foreshock.Circles represent aftershocks with magnitude proportional to its size.
5.3 重定位的震源深度分布及其意义分析

2014年于田地震发生后,国内外不同机构给出了此次于田MS7.3地震及其前震的震源深度(表 1),大部分结果在10~12 km之间.对于前震和主震,GCMT给出的矩心深度分别为21.7 km和18.3 km,本文重定位给出的震源深度分别为19 km和22 km.图 3图 8a显示,地震目录给出于田地震序列的震源深度主要集中在5~10 km(约占73%),分布在0~5 km的约占18%.本文重定位后的地震序列平均震源深度为14.8 km(图 7图 8b),约80%的地震事件分布在6~25 km的深度范围内,约11%分布在0~5 km,大部分余震的震源深度浅于主震的震源深度,其中地震序列中的46次ML≥3.5余震中有30次地震的深度位于15~26 km,在主震上方的地震极少,表明主震破裂使得破裂区能量释放较为充分,而在地表破裂带两侧地震分布在深度上较为均匀(图 6图 7).应用PTD方法测得此次地震序列106次地震的平均震源深度为10.6 km,其中有70次地震的震源深度介于10~20 km(宋秀青等,2014).使用双差定位方法对这次地震序列进行重定位得到的震源深度主要集中在4~12 km(房立华等,2015).使用CAP方法反演获得此次地震序列前震、主震、及早期MS≥3.5共18次地震的平均震源深度为15.6 km,大部分地震的震源深度为15~20 km(李金和王琼,2015).上述不同研究得到震源深度的差异主要是由于台站分布稀疏,缺乏足够近的震相数据,导致地震定位时对深度的约束较差(唐明帅和王海涛,2011),以及使用不同的速度模型和定位方法时,可能导致较大的深度差异(房立华等,2013).

图 7 震源深度分布剖面 剖线位置如图 5,沿AB剖线每个地震到剖线的距离小于8 km,到其余剖线的距离小于5 km. Fig. 7 The distributions of focal depth along different profiles Profile′s locations and legends are the same as in Fig.5.The distance between each earthquake and the axis of the cross-section is less than 5 km,along profile AB less than 8 km.
图 8 重定位前(a)、后(b)的震源深度分布柱状图 Fig. 8 Histogram of focal depth distribution before(a)and after(b)earthquake relocation

图 7(bce)显示,本文重定位后地震在深度方向具有近垂直分布的特征,这与阿尔金走滑断裂倾角较陡的特征一致.图 5图 7(bc))显示分布于南肖尔库勒断裂南边的地震多于北边的,同时可推断南肖尔库勒断裂的倾向为SE,这与本次地震野外考察(李海兵等,2015)得到的结论(沿南硝尔库勒断裂分布的地表破裂带不仅具有左旋走滑位移,而且反映出挤压性质,具有逆冲分量,断裂南盘略有上升,因此,推测地震断裂面应向南陡倾)一致.图 5图 7e显示分布于阿什库勒—硝尔库勒断裂南边的地震多于北边的,从图 7e可推断阿什库勒—肖尔库勒断裂的倾向为NW,与本次地震野外考察(李海兵等,2015)给出的结论(沿阿什库勒—肖尔库勒断裂分布的同震地表破裂带不仅具有左旋走滑性质,而且还具有伸展性质,沿断裂发育较多的张性裂缝,并且断裂北盘略带有下降,反映出地震断裂面向NNW方向陡倾,反映出走滑-伸展性质)一致.上述结论说明本文重定位给出了地震序列可靠的空间分布特征.

本文重定位后的震源深度分布下界约为26 km(图 7),而震源深度的下界面有可能取决于脆-韧性转换带的深度.震区附近的壳内间断面深度约为24.8 km,与重定位后的震源深度的下界一致,我们推断在震区地壳内约25 km深度位置可能存在脆-韧转换带.这一结果对于确定阿尔金断裂西南段地壳内的发震层厚度,阐明地震的成因和机制以及地震危险性分析提供了重要的约束条件(张国民和李丽,2003).而位于帕米尔东北缘的西克尔地震台附近的壳内间断面和震源深度的下界面的深度约为21 km(赵翠萍,2006唐明帅等,2014),说明位于新疆不同区域,可能存在不同深度的脆-韧转换带.研究表明(李海兵等,2015),2014年于田MS7.3地震震区的断裂与郭扎错和龙木错断裂构成“阿尔金断裂”向西南方向的延伸部分,表明这些断裂的发育成熟度低于阿尔金主断裂,这可能在一定程度上造成了地壳强度上的差异(宋美琴等,2012),从而造成了位于地表破裂带北东方向余震的深度下限略深于西南方向余震深度这一现象.

6 结论

本文首先研究了位于震源区附近于田地震台下方的地壳结构,建立了震源区的地壳速度模型.在此基础上,联合震相到时和方位角对2014年于田MS7.3地震序列进行了重定位,获得以下结论.

(1) 本文选取位于于田台东南方位的接收函数获得台站下方的地壳厚度为60.8 km、地壳波速比为1.65,壳内间断面深度为24.8 km、波速比为1.64,以及相应的速度结构.

(2) 重定位后的前震和主震震中位置明显向地表破裂带、阿什库勒—肖尔库勒断裂和南肖尔库勒断裂靠近,两者相距5.4 km,前震位置为36.055°N、82.522°E,震源深度为19 km;主震位置为36.076°N、82.576°E,震源深度为22 km.

(3) 重定位后的震中优势沿NEE-SWW方向展布,与此次地震地表破裂带和附近的断裂走向一致.余震序列长度约73 km、宽度约16 km,平均震源深度为14.8 km.77%的余震分布在地表破裂带的西南部、沿南肖尔库勒断裂分布;23%的余震分布在地表破裂带东北部、沿阿什库勒—肖尔库勒断裂分布.位于地表破裂带西南部的震中分布格局变化较大,由台网目录的近南北向变为北东向,与地表破裂带、南肖尔库勒断裂和阿什库勒—肖尔库勒断裂走向一致.

(4) 根据重定位后沿剖面的地震分布(图 7),可推断位于西南段的南肖尔库勒断裂与北东段的阿什库勒—肖尔库勒断裂倾向反向,南肖尔库勒断裂的倾向为SE,阿什库勒—肖尔库勒断裂的倾向为NW,这与本次地震野外考察得到的断层性质一致.

(5) 综合重新定位结果、地表破裂带分布、震源机制解、南肖尔库勒断裂和阿什库勒—肖尔库勒断裂的性质分析认为,2014年于田MS7.3地震的发震构造为阿尔金断裂西南段分支断裂——南肖尔库勒断裂和阿什库勒—肖尔库勒断裂.

致谢

本文大部分图件使用GMT软件绘制(Wessel and Smith,1991),特此申明.在本文成文过程中,万永革教授提供了利用小震分布确定大震断层面参数的计算程序,张勇博士提供了此次地震的破裂过程数据,葛粲博士和冀战波博士给予了无私帮助,在利用定位结果分析断层性质时与吴传勇博士进行了有益的讨论,在此对他们深表感谢.感谢审稿专家对本文细致的审阅和提出富有建设性的修改意见,提高了文章的质量.

参考文献
Aki K, Richards P G. Quantitative Seismology: Theory and Methods. San Francisco: W.H.Freeman and Co, 1980 : 28 -119.
Ammon C J, Randall G E, Zandt G. 1990. On the nonuniqueness of receiver function inversions. J.Geophys.Res. , 95(B10): 15303–15318.
Bratt S R, Bache T C. 1988. Locating events with a sparse network of regional arrays. Bull.Seism.Soc.Am. , 78(2): 780–798.
Chang C H, Wu Y M, Shin T C, et al. 2000. Relocation of the 1999 Chi-Chi Earthquake in Taiwan. Terr.Atmos.Oceanic Sci. , 11(3): 581–590.
Fang L H, Wu J P, Wang W L, et al. 2013. Relocation of the mainshock and aftershock sequences of MS7.0 Sichuan Lushan earthquake. Chin.Sci.Bull. , 58(28): 3451–3459. doi: 10.1007/s11434-013-6000-2.
Fang L H, Wu J P, Wang W L, et al. 2015. Relocation of the 2014 MS7.3 earthquake sequences in Yutian,Xinjiang. Chinese J.Geophys. , 58(3): 802–808. doi: 10.6038/cjg20150310.
Geiger L C.1910.Herdbestimmung bei Erdbeben aus den Ankunftszeiten.Nachrichten von der Gesellschaft der Wissenschaften zu G ttingen,Mathematisch-Physikalische Klasse,331-349.
Geiger L C.2012.Probability Method for the Determination of Earthquake Epicenters from the Arrival Time Only(in Chinese).Zhao Z H(translated).China Earthquake Networks Center.
He R Z, Gao R, Li Q S, et al. 2001. Corridor gravity fields and custal density structures in Tianshan(Dushanzi)-west Kunlun(Quanshuigou)GGT. Acta Geoscientia Sinica (in Chinese) , 22(6): 553–558.
Herrmann R B.2004.Computer programs in seismology.Saint Louis University,71-90, http://www.eas.slu.edu/eqc/eqccps.html.
Huang Y, Wu J P, Zhang T Z, et al. 2008. Relation of the MS8.0 Wenchuan earthquake and its aftershock sequence. Science in China(Series D) (in Chinese) , 38(10): 1242–1249.
Ji Z B, Zhao C P, Wang Q, et al. 2014. Source rupture process of Yutian,Xinjiang,MS7.3 earthquake on 21 March 2008. Acta Seismologica Sinica (in Chinese) , 36(3): 339–349.
Langston C A. 1979. Structure under Mount Rainier,Washington,inferred from teleseismic body waves. J.Geophys.Res. , 84(B9): 4749–4762.
Li H B, Sun Z M, Pan J W, et al. 2014. Field study of the 12 February 2014 Yutian MS7.3 earthquake: A special surface rupture zone. Acta Geoscientica Sinica (in Chinese) , 35(3): 391–394.
Li H B, Pan J W, Sun Z M, et al. 2015. Seismogenic structure and surface rupture characteristics of the 2014 MS7.3 Yutian earthquake. Acta Geologica Sinica (in Chinese) , 89(1): 180–194.
Li J, Wang Q. 2015. The analysis on the characteristics of the focal mechanisms for the MS7.3 earthquake sequence on February 12,2014 in Yutian. Earthquake Research in China (in Chinese) , 31(1): 110–120.
Li Q S, Lu D Y, Gao R, et al. 2000. Explosion seismic probing across the West Kunlun-Tarim contact zone. Science in China(Series D) (in Chinese) , S: 16–21.
Li S P. Earthquake China. (in Chinese) Beijing: Earthquake Press, 1981 : 120 -121.
Paul Wessel, Walter H, Smith F. 1991. Free Software Helps Map and Display Data. EOS , 72(41): 441–448.
Randall G E. 1989. Efficient calculation of differential seismograms for lithospheric receiver functions. Geophys.J.Int. , 99(3): 469–481.
Schweitzer J. 2001. HYPOSAT-An enhanced routine to locate seismic events. Pure Appl.Geophys. , 158(1): 277–289.
Song M Q, Zheng Y, Ge C, et al. 2012. Relocation of small to moderate earthquakes in Shanxi Province and its relation to the seismogenic structures. Chinese J.Geophys. , 55(2): 513–525. doi: 10.6038/j.issn.0001-5733.2012.02.014.
Song X Q, Miao F J, Liu S Q, et al. 2014. The depth determination of the 2014 Yutian MS7.3 earthquake sequence,Xinjiang. Earthquake Research in China (in Chinese) , 30(2): 198–207.
Tang M S, Wang H T, Duan T S. 2010. Relocating the seismic sequence of Yutian MS7.3 earthquake in 2008 using the data of Xinjiang Hetian seismic array. Inland Earthquake (in Chinese) , 24(3): 227–235.
Tang M S, Wang H T. 2011. F-detection algorithm and its preliminary application to seismic surface reflected wave phase identification. Acta Seismologica Sinica (in Chinese) , 33(6): 776–787.
Tang M S, Zheng Y, Ge C, et al. 2014. Study on crustal structure in the northeastern Pamir region by P receiver functions. Chinese J.Geophys , 57(10): 3176–3188. doi: 10.6038/cjg20141007.
Waldhauser F, Ellsworth W L. 2000. A double-difference earthquake location algorithm: Method and application to the Northern Hayward Fault,California. Bull.Seism.Soc.Am. , 90(6): 1353–1368.
Wan Y G, Shen Z K, Diao G L, et al. 2008. An algorithm of fault parameter determination using distribution of small earthquakes and parameters of regional stress field and its application to Tangshan earthquake sequence. Chinese J.Geophys. (in Chinese) , 51(3): 793–804.
Wang J, Song X Q, Chen X J, et al. 2014. Study on precision positioning the Yutian,Xinjiang MS7.3 earthquake. Earthquke Research in China , 30(2): 188–197.
Wang W M, Li L, Zhao L F, et al. 2005. Rupture process of Jiashi,Xinjiang earthquake(MS6.5)of Feb.24,2003. Chinese Journal of Geophysics (in Chinese) , 18(2): 343–351.
Wang X X, Liu R F, Zou L Y, et al. 2014. Determination of the parameters for the 12 February 2014 Yutian MS7.3 earthquake in Xinjiang Uygur Autonomous Region. Acta Seismologica Sinica (in Chinese) , 36(3): 522–528.
Wu C Y, Zhang Z Q, Zhao C P, et al. 2014. 2014 Yutian MS7.3 earthquake: Structural response of the Bayankala tectonic-block to eastward extrusion. Chinese J.Geophys , 57(10): 3226–3237. doi: 10.6038/cjg20141011.
Earthquake Administration of Xinjiang Uygur Autonomous Region. 2014. Disaster management of the 12 February 2014 Yutian MS7.3 earthquake. Earthquake Administration of Xinjiang Uygur Autonomous Region : 1–134.
Zhang G M, Wang S Y, Li L, et al. 2002. Focal depth research of earthquakes in mainland China: Implication for tectonics. Chin.Sci.Bull. , 47(12): 969–974.
Zhang G M, Li L. 2003. Rheology of crustal media and a related seismogenic model. Seismology and Geology (in Chinese) , 25(1): 1–10.
Zhang G W, Lei J S, Sun C Q. 2014. Relocation of the 12 February 2014 Yutian,Xinjiang,mainshock(MS7.3)and its aftershock sequence. Chinese J.Geophys , 57(3): 1012–1020. doi: 10.6038/cjg20140330.
Zhang Y, Xu L S, Chen Y T, et al. 2014. Fast inversion for the rupture process of the 12 February 2014 Yutian MW6.9 earthquake: Discussion on the impacts of focal mechanisms on rupture process inversions. Acta Seismologica Sinica (in Chinese) , 36(2): 159–164.
Zhao C P.2006.Seismological studies on the characteristics of Jashi source region from 1997 to 2003[Ph.D.thesis](in Chinese).Beijing: Institute of Geophysics,China Earthquake Administration,1-58.
Zhao C P, Chen Z L, Zheng S H. 2008. Source rupture process of 3 Jiashi MS6 events(1998-2003)and its correlation with the aftershock activity. Chinese Journal of Geophysics , 51(4): 1093–1102. doi: 10.3321/j.issn:0001-5733.2008.04.018.
Zhou Y, Wang W M, Xiong L, et al. 2015. Rupture process of 12 February 2014,Yutian MW6.9 earthquake and stress change on nearby faults. Chinese J.Geophys. , 58(1): 184–193. doi: 10.6038/cjg20150116.
Zhu L P, Kanamori H. 2000. Moho depth variation in Southern California from teleseismic receiver functions. J.Geophys Res , 105: 2969–2980.
房立华, 吴建平, 王未来, 等. 2013. 四川芦山MS7.0级地震及其余震序列重定位. 科学通报 , 58(20): 1901–1909.
房立华, 吴建平, 王未来, 等. 2015. 2014年新疆于田MS7.3级地震序列重定位. 地球物理学报 , 58(3): 802–808.
Geiger L C.2012.只由到时确定地震震中的概率方法.赵仲和译.中国地震局台网中心.
贺日政, 高锐, 李秋生, 等. 2001. 新疆天山(独山子)-西昆仑(泉水沟)地学断面地震与重力联合反演地壳构造特征. 地球学报 , 22(6): 553–558.
黄媛, 吴建平, 张天中, 等. 2008. 汶川8.0级大地震及其余震序列重定位研究. 中国科学D辑:地球科学 , 38(10): 1242–1249.
冀战波, 赵翠萍, 王琼, 等. 2014. 2008年3月21日新疆于田MS7.3地震破裂过程研究. 地震学报 , 36(3): 339–349.
李海兵, 孙知明, 潘家伟, 等. 2014. 2014年于田MS7.3地震野外调查-特殊的地表破裂带. 地球学报 , 35(3): 391–394.
李海兵, 潘家伟, 孙知明, 等. 2015. 2014年于田MS7.3地震地表破裂特征及其发震构造. 地质学报 , 89(1): 180–194.
李金, 王琼. 2015. 2014年2月12日于田7.3级地震序列震源机制特征分析. 中国地震 , 31(1): 110–120.
李秋生, 卢德源, 高锐, 等. 2000. 横跨西昆仑-塔里木接触带的爆炸地震探测. 中国科学: D辑 , 30(增刊): 16–21.
李善邦. 中国地震. 北京: 地震出版社, 1981 : 120 -121.
宋美琴, 郑勇, 葛粲, 等. 2012. 山西地震带中小震精确位置及其显示的山西地震构造特征. 地球物理学报 , 55(2): 513–525.
宋秀青, 缪发军, 刘双庆, 等. 2014. 2014年新疆于田MS7.3地震序列的震源深度测定. 中国地震 , 30(2): 198–207.
唐明帅, 王海涛, 段天山. 2010. 利用和田地震台阵数据对2008年于田7.3级地震序列重新定位. 内陆地震 , 24(3): 227–235.
唐明帅, 王海涛. 2011. F检测算法及其在识别地震地表反射波震相中的初步应用. 地震学报 , 33(6): 776–787.
唐明帅, 郑勇, 葛粲, 等. 2014. 帕米尔东北缘地壳结构的P波接收函数研究. 地球物理学报 , 57(10): 3176–3188.
万永革, 沈正康, 刁桂苓, 等. 2008. 利用小震分布和区域应力场确定大震断层面参数方法及其在唐山地震序列中的应用. 地球物理学报 , 51(3): 793–804.
王俊, 宋秀青, 陈向军, 等. 2014. 新疆于田MS7.3地震主震精定位研究. 中国地震 , 30(2): 188–197.
王卫民, 李丽, 赵连锋, 等. 2005. 2003年2月24日新疆伽师MS6.5级地震震源破裂过程研究. 地球物理学报 , 48(2): 343–351.
王晓欣, 刘瑞丰, 邹立晔, 等. 2014. 2014年2月12日新疆于田MS7.3地震参数测定. 地震学报 , 36(3): 522–528.
吴传勇, 张竹琪, 赵翠萍, 等. 2014. 2014年新疆于田MS7.3级地震:巴颜喀喇地块侧向挤出的构造响应. 地球物理学报 , 57(10): 3226–3237.
新疆维吾尔自治区地震局.2014.2014年2月12日新疆于田7.3级地震灾害损失评估.新疆维吾尔自治区地震局,1-134.
张国民, 汪素云, 李丽, 等. 2002. 中国大陆地震震源深度及其构造含义. 科学通报 , 47(9): 663–668.
张国民, 李丽. 2003. 地壳介质的流变性与孕震模型. 地震地质 , 25(1): 1–10.
张广伟, 雷建设, 孙长青. 2014. 2014年2月12日新疆于田MS7.3级地震主震及余震序列重定位研究. 地球物理学报 , 57(3): 1012–1020.
张勇, 许力生, 陈运泰, 等. 2014. 2014年2月12日于田MW6.9地震破裂过程初步反演:兼论震源机制对地震破裂过程反演的影响. 地震学报 , 36(2): 159–164.
赵翠萍.2006.1997-2003年新疆伽师震源区特征的地震学方法研究[博士论文].北京:中国地震局地球物理研究所,1-58.
赵翠萍, 陈章立, 郑斯华. 2008. 1998-2003年伽师三次不同类型MS6级地震震源破裂过程及短期内余震活动特征. 地球物理学报 , 51(4): 1093–1102.
周云, 王卫民, 熊林, 等. 2015. 2014年2月12日MW6.9于田地震震源破裂过程及对周围断层的应力影响. 地球物理学报 , 58(1): 184–193.