孫 赫 季靈運 朱良玉 趙 強

1 中國地震局第二監(jiān)測中心，西安市西影路316號，710054

?

新疆于田MS7.3地震同震與震后形變機制研究

孫 赫1季靈運1朱良玉1趙 強1

1 中國地震局第二監(jiān)測中心，西安市西影路316號，710054

利用InSAR技術(shù)獲取2008-03-21新疆于田MS7.3地震的同震和震后形變場。同震分布式滑動反演結(jié)果表明，同震斷層最大滑動量達5.4 m，主要分布在南部斷層的0～5 km深度附近，地震以正斷錯動為主，兼有左旋走滑分量。震后形變結(jié)果表明，發(fā)震斷層北段兩側(cè)存在差異性運動，最大累積差異形變在震后782 d達15 cm。進一步分析表明，震后斷層余滑可能是震后形變的主要機制。余滑反演結(jié)果表明，震后2 a斷層余滑量相對較小，滑移區(qū)范圍明顯減小且均位于淺部區(qū)域，北部斷層能量釋放較徹底，南部仍存有少量能量，整體能量基本釋放完。

于田地震；InSAR；同震形變場；震后形變場；余滑

2008-03-21新疆于田MS7.3地震[1]震中位于西昆侖地震帶與阿爾金地震帶交匯區(qū)，地表破裂長度達30多km，在正斷層的基礎上存在左旋走滑位移[2-7]。區(qū)域構(gòu)造如圖1所示。

紅色矩形框為本文研究區(qū)域；黑色線為區(qū)域內(nèi)斷層分布；紅色實線為野外地質(zhì)考察確定的地表破裂[2]；地形數(shù)據(jù)為90 m分辨率的SRTM數(shù)據(jù)；藍色圓圈為NEIC震后5個月的余震分布圖1 于田地震區(qū)域構(gòu)造示意圖Fig.1 Local faults of the 2008 Yutian earthquake

地震發(fā)生后，許多專家對其發(fā)震機制與發(fā)震斷層構(gòu)造進行多方面研究，其中包括野外地質(zhì)考察[2]、地質(zhì)構(gòu)造與地殼形變研究[5]以及基于GPS、InSAR技術(shù)的探測與模擬[6-9]等。本文利用SAR數(shù)據(jù)獲取的InSAR形變場反演于田地震的同震與震后斷層滑動分布。對比分析后發(fā)現(xiàn)，震后余滑量相對較小，能量基本釋放完全。

1 D-InSAR同震形變結(jié)果與分析

1.1 反演算法簡介

本文利用Wang等[10]提出的SBIF(sensitivity based iterative fitting)反演程序。在反演處理之前，采用四叉樹方法將InSAR技術(shù)獲取的形變場數(shù)據(jù)進行降采樣[11]，以提高運算速度。反演計算時，采用InSAR技術(shù)獲取的形變場作為約束。發(fā)震斷層走向通過震后野外實地考察確定，傾角根據(jù)ISC提供的震源機制解和文獻[7]綜合確定。斷層主要包括南北兩段，具體位置與方向如圖2中黑色實線所示。北段走向為196.4°，傾角為52°，南段走向為漸變，由北向南從155°到230°，傾角為52°。為了得到較合理的反演結(jié)果，利用數(shù)據(jù)擬合程度與斷層滑動粗糙度曲線之間的關(guān)系確定該斷層的平滑因子。圖3為于田地震同震斷層應力面的粗糙度與數(shù)據(jù)擬合曲線，分析斷層應力面平滑條件，得到同震的平滑因子為0.05。2008年于田地震震源深度為19 km，而且大部分影響深度為地表以下10 km處，故設置反演模型深度為20 km。

1.2 同震形變結(jié)果與分析

利用D-InSAR技術(shù)分別對ENVISAT SAR數(shù)據(jù)和ALOSPALSAR數(shù)據(jù)(表1)進行處理，并采用90 m分辨率的SRTM數(shù)據(jù)作為外部DEM數(shù)據(jù)，得到于田地震同震形變場(圖2(a)和2(b)，形變方向為LOS方向)。研究時間域內(nèi)，在冰雪覆蓋區(qū)及部分地勢陡峭的山區(qū)出現(xiàn)失相干現(xiàn)象。加之衛(wèi)星軌道入射角差異，升降軌數(shù)據(jù)獲取的干涉圖信息不盡相同。升軌PALSAR數(shù)據(jù)獲取的最大形變區(qū)位于上盤，沉降量達110 cm，下盤沉降約為50 cm。降軌ENVISAT數(shù)據(jù)在最大形變區(qū)失相干，上、下盤最大沉降與抬升均為40 cm左右。

圖2 2008年于田地震形變結(jié)果(包括ENVISAT與PALSAR兩類數(shù)據(jù))Fig.2 The deformation results of 2008 Yutian earthquake including ENVISAT and PALSAR

表1 SAR影像數(shù)據(jù)參數(shù)

反演過程中，首先對獲得的InSAR形變場進行降采樣，以提高運算速度；然后對ENVISAT和PALSAR同震形變場進行共同約束反演，得到模擬形變場(圖2(e)～圖(f))。殘差結(jié)果顯示，差值在±10 cm之間，精度比較可靠?；瑒臃植冀Y(jié)果(圖4)顯示，滑移主要集中在地下0～5 km深度，南段斷層的滑移量比較大，最大達5.4 m；北段斷層滑移量約為4 m，集中分布在地下10～18 km深度。斷層滑動分布證實了于田地震的正斷層破裂機制，同時兼有一定的左旋走滑分量。

圖4 于田地震同震反演的斷層面滑動分布Fig.4 The inversion slip distribution of coseismic for Yutian earthquake

2 震后形變場獲取與機制探討

為獲取震后形變信息，選用9景L波段的ALOS PALSAR數(shù)據(jù)，時間跨度為2008-05-26～2010-07-17。其中,垂直基線的閾值為±2 500 m，時間基線的閾值為35～800 d。通過影像數(shù)據(jù)的自由組合，選取11個相干性較高的干涉對進行計算，詳細參數(shù)見表2，時-空基線關(guān)系如圖5所示。

采用SBAS技術(shù)對表2中的干涉對進行處理，得到于田地震震后累積形變場，如圖6所示。2008-05-26～2010-07-17發(fā)震斷層北段兩側(cè)存在差異形變，累積差異形變約為15 cm，表明震后斷層存在活動。為了定量分析斷層兩側(cè)地表形變變化，分析圖6中剖線在各時間段內(nèi)的形變變化(圖7)。圖7顯示，震后46 d內(nèi)基本沒有形變，184 d時發(fā)生形變，直到506 d后形變基本穩(wěn)定，說明斷層滑移基本停止。圖8給出了基于指數(shù)函數(shù)和對數(shù)函數(shù)的斷層北段附近a點的形變時間序列擬合曲線，結(jié)果表明，震后2 a內(nèi)形變逐年增加，累積超過5 cm。

表2 選取干涉對列表信息

圖5 干涉對基線組合Fig.5 Baseline combination of all the interferomeric pairs

圖6 于田地震震后2008～2010年累積形變(虛線為形變時間序列剖面位置)Fig.6 The postseismic cumulative deformation of Yutian earthquake from 2008 to 2010(Dotted lines represent profile position for the time series of deformation)

圖6～8的震后累積形變及震后形變時間序列結(jié)果表明，震后斷層北段存在余滑。為了得到震后斷裂帶的滑動情況，同樣采用上文的反演策略與算法，以獲取的震后InSAR地殼形變?yōu)榧s束，對震后InSAR降采樣形變場(圖9(a))進行反演，得到模擬形變場(圖9(b))與斷裂帶滑動分布(圖10)。研究區(qū)域內(nèi)殘差結(jié)果差異在±3 cm內(nèi)。由圖10知，震后2 a內(nèi)，斷層滑動量相對很小，最大滑動量為0.7 m，南北兩段斷層的滑移區(qū)均集中在地下0～5 km深度。結(jié)合同震與震后斷層滑動分布結(jié)果，北段斷層能量由深至淺進行釋放，震后2a基本釋放完；南段斷層震后滑動區(qū)范圍明顯減小，同時余滑量也明顯減弱，由同震的5.4 m降到震后的0.7 m。

圖7 2008～2010年剖線上各時段累積形變(剖線位置見圖6)Fig.7 Accumulative deformation of profile line in different stages from 2008 to 2010 (profile position in Fig.6)

圖8 2008-2010年于田地震震后a點形變時間序列Fig.8 Time series of postseismic deformation in LOS for point a of Yutian earthquack in different stages from 2008 to 2010

圖9 InSAR降采樣形變場、形變場模擬結(jié)果及殘差Fig.9 Resampled deformation, inversed deformation and residual for inversed data

圖10 于田地震震后斷層面的滑動分布Fig.10 The inversion slip distribution of postseismic for Yutian earthquake

震例研究表明，震后地殼形變的機制主要有3種，即震后余滑、孔隙彈性回彈、震后粘彈性松弛[12]。震后余滑產(chǎn)生的形變主要集中在斷層附近，余滑一般主要發(fā)生于斷裂中部并且局限在上部4 km[13]; 孔隙彈性回彈產(chǎn)生孔隙彈性回彈形變，主要影響近場；而震后粘彈性松弛影響的范圍比較寬廣。在震后1～2 a，3種因素都有可能產(chǎn)生震后形變。但是，分析于田地震的震后形變場發(fā)現(xiàn)，形變主要發(fā)生在發(fā)震斷層北段，并且沿斷層兩側(cè)形變差異明顯，比較符合震后余滑機制。后續(xù)反演分析亦以余滑機制為主，但是也可能有其他兩種機制的貢獻。

3 結(jié) 語

采用InSAR技術(shù)獲取的同震及震后形變場顯示，2008年于田地震符合正斷層的形變機制兼有左旋走滑分量。利用SBIF反演方法分別對獲取的形變場進行模擬反演計算，得到相應的模擬形變場及斷層滑動分布。同震南段斷層活動強烈，分布在地下深度0～15 km，最大滑移量為5.4 m。震后2 a時間內(nèi)發(fā)震斷層北段兩側(cè)形變差異顯著，震后形變的主要機制可能是由震后斷層余滑所致。震后斷層活動較同震弱，能量基本釋放完全，最大余滑量約為0.7 m。南北兩段斷層的滑移區(qū)均位于地下0～5 km的淺部區(qū)域。

[1] 尹光華，蔣靖祥，吳國棟. 2008年3月21日于田7.4級地震的構(gòu)造背景[J]. 干旱區(qū)地理，2008，31(4)：543-549(Yin Guanghua，Jiang Jingxiang，Wu Guodong. Tectonic Background of the MS7.4 Earthquake at Yutian on March 21, 2008[J]. Arid Land Ceography，2008，31(4)：543-549)

[2] 徐錫偉，譚錫斌，吳國棟，等. 2008年于田MS7.3地震地表破裂帶特征及其構(gòu)造屬性討論[J]. 地震地質(zhì)，2011，33(2)：462-471(Xu Xiwei，Tan Xibin，Wu Guodong，et al. Surface Ruture Features of the 2008 Yutian MS7.3 Earthquack and Its Tectonic Nature[J]. Seismology and Geology，2011，33(2)：462-471)

[3] 洪順英，申旭輝，單新建，等. 基于升降軌ASAR的于田MS7.3級地震同震形變場信息提取與分析[J]. 國土資源遙感，2010，87(4)：98-102(Hong Shunying，Shen Xuhui，Shan Xinjian，et al. The Calculation and Analysis of the Coseismic Deformation Field of Yutian MS7.3 Earthquake Basing on the Ascending and Descending Orbit ASAR Date[J]．Remote Sensing for Land and Resources，2010，87(4)：98 -102)

[4] 張廣偉，雷建設，孫長青．2014年2月12日新疆于田MS7.3級地震主震及余震序列重定位研究[J]．地球物理學報，2014，57(3)：1 012-1 020(Zhang Guangwei，Lei Jianshe，Sun Changqing. Relocation of the 12 February 2014 Yutian, Xinjiang, Mainshock (MS7.3) and Its Aftershock Sequence[J]. Chinese Journal of Geophysics，2014，57(3)：1 012-1 020)

[5] 萬永革，沈正康，盛書中，等. 2008年新疆于田7.3級地震對周圍斷層的影響及其正斷層機制的區(qū)域構(gòu)造[J]. 地球物理學報，2010，53(2)：280-289(Wan Yongge，Shen Zhengkang， Sheng Shuzhong，et al. The Mechanical Effects of the 2008 MS7.3 Yutian，Xinjiang Earthquake on the Neighboring Faults and Its Tectonic Origin of Normal Faulting Mechanism[J]. Chinese Journal of Geophysics，2010，53(2): 280-289)

[6] 王凡，陳為濤，王敏，等. 利用GPS觀測資料分析2008年于田MS7.3地震的同震位移及震后形變[J]. 地球物理學報，2011，54(9)：2 250-2 255(Wang Fan，Chen Weitao，Wang Min，et al. Coseismic and Postseismic Deformation of MS7.3 Yutian Earthquake Derived from GPS Data[J]. Chinese Journal of Geophysics，2011，54(9)：2 250-2 255)

[7] 張國宏，屈春燕，單新建，等. 2008年MS7.1于田地震InSAR同震形變場及其震源滑動反演[J].地球物理學報，2011，54(11)：2 753-2 760(Zhang Guohong，Qu Chunyan，Shan Xinjian，et al. The Coseimic InSAR Measurements of 2008 Yutian Earthquake and Its Inversion for Source Parameters[J]. Chinese Journal of Geophysics，2011，54(11): 2 753-2 760)

[8] Lasserre G，Peltzer C，Crampe F，et al. Coseismic Deformation of the 2001 MW7.8 Kokoxili Earthquake in Tibet，Measured by Synthetic Aperture Radar Interferometry[J]. Journal of Geophysical Research Atmospheres，2005，110(B12)：1-18

[9] Shen Z K，Sun J B，Zhang P Z，et al. Slip Maxima at Fault Junctions and Rupturing of Barriers During the 2008 Wenchuan Earthquake[J]. Nature Geoscience，2009，2(10)：718-724

[10]Wang R，Xia Y，Grosser H，et al. The 2003 Bam(SE Iran)Earthquake：Precise Source Parameters from Satellite Radar Interferometry[J]. Geophysical Journal International，2004，159(3)：917-922

[11]Lohman R B，Simons M. Some Thoughts on the Use of InSAR Data to Constrain Models of Surface Deformation：Noise Structure and Data Downsampling[J]. Geochemistry Geophysics Geosystems，2005，6(1)：359-361)

[12]譚凱，王琪，王曉強，等. 震后形變的解析模型和時空分布特征[J]. 大地測量與地球動力學，2005，25(4)： 23-26(Tan Kai，Wang Qi，Wang Xiaoqiang，et al. Analyticmodels and Space-Time Distribution of Postseismic Deformation[J]. Journal of Geodesy and Geodynamics，2005，25(4)：23-26)

[13]Scholz C H，Wyss M，Smith S W. Seismic and Aseismic Slip on the San Andreas Fault[J]. Journal of Geophysical Research Atmospheres，1969，74(8)：2 049-2 069

About the first author:SUN He，assistant engineer，majors in crustal deformation and earthquake prediction，E-mail：sunhele@163.com.

Coseismic Deformation and Postseismic Deformation Mechanism Yutian MS7.3 Earthquake

SUNHe1JILingyun1ZHULiangyu1ZHAOQiang1

1 Second Crust Monitoring and Application Center, CEA,316 Xiying Road，Xi’an 710054，China

On March 21th, 2008, a MS7.3 earthquake occurred in Yutian of Xinjiang. In this paper, InSAR technology is used to obtain the coseismic and postseismic deformation. The inversion of coseismic slip distribution shows that the largest slip achieves 5.4 m and that the slip is dominant in the south of the fault, near 0-5 km depth beneath the ground surface. The earthquake is a normal faulting event with a slight sinistral strike-slip component. Postseismic deformation indicates that there are differences in movement in the north section of seismogenic fault, and the largest cumulative deformation difference is about 15 cm in 782 d after the earthquake. According to the further analysis, the afterslip of postseismic fault may be the main mechanism for deformation. The degree of postseismic afterslip is relatively small in two years and the slip region decreases greatly and for the most part is located in shallow areas. The energy of north fault is relatively exhaustive released, while there is little energy remaining in south part of the fault.

Yutian earthquake；InSAR；coseismic deformation；postseismic deformation；afterslip

Special Fund for Earthquake Research of CEA, No.201508009.

2015-12-15

項目來源：中國地震局地震行業(yè)科研專項(201508009)。

孫赫，助理工程師，主要從事地殼形變與地震預測研究，E-mail：sunhele@163.com。

10.14075/j.jgg.2016.12.004

1671-5942(2016)012-1052-04

P315

A