王　輝　曹建玲　洪順英　徐岳仁　荊　鳳

1)中國地震局地震預(yù)測重點實驗室，北京　100036 2)中國地震局地震預(yù)測研究所，北京　100036

?

2008年和2014年2次新疆于田M7地震之間的黏彈性應(yīng)力轉(zhuǎn)移

王輝1，2)曹建玲1，2)洪順英2)徐岳仁2)荊鳳2)

1)中國地震局地震預(yù)測重點實驗室，北京100036 2)中國地震局地震預(yù)測研究所，北京100036

青藏高原西北邊緣地區(qū)在2008年和2014年先后發(fā)生了2次7級以上強震，這2次地震發(fā)生在阿爾金斷裂帶西南端，對于認識區(qū)域動力背景和潛在地震危險性有重要意義。文中采用分層黏彈性模型和最大應(yīng)力策略研究了這2次地震之間的相互關(guān)系。計算結(jié)果表明，2008年于田地震的同震應(yīng)力擾動觸發(fā)了后續(xù)的余震活動，該地震產(chǎn)生的同震應(yīng)力擾動對2014年于田地震的發(fā)生影響不顯著；但是，此次地震6a后震后變形所導致的應(yīng)力擾動有利于2014年于田地震的發(fā)生。2014年于田地震發(fā)生后，阿爾金斷裂帶西南段上的同震和震后庫倫應(yīng)力顯著增加。考慮到阿爾金斷裂帶西南段上次強震的離逝時間和應(yīng)力積累-釋放過程，該斷裂帶上的地震危險性值得進一步關(guān)注。

于田地震分層黏彈性模型庫倫破裂應(yīng)力阿爾金斷裂帶西南段潛在地震危險性

0　引言

新疆于田地區(qū)位于青藏高原西北邊緣，該地區(qū)是青藏塊體和塔里木塊體交會的地區(qū)。柯崗斷裂帶、阿爾金斷裂帶、康西瓦斷裂帶等大型斷裂帶在此交會，構(gòu)造活動強烈。雖然該地區(qū)缺乏地震歷史記載，但是在過去數(shù)年中，該地區(qū)先后發(fā)生了2008年3月21日M7.3地震和2014年2月12日M7.3地震(圖1)。這2次地震分別發(fā)生在阿爾金斷裂帶西端的不同分支上，宏觀震中相距約200km，震源機制類型分別屬于正斷層型和走滑兼正傾滑型。如此短的時間間隔中、如此近的空間范圍內(nèi)發(fā)生的2次M7地震為研究強震之間的相互作用及其對區(qū)域潛在地震危險性的影響提供了很好的科學實例(Zhuetal.，2015)。

圖1　于田地區(qū)活動構(gòu)造簡圖Fig. 1　Simplified tectonic map in the Yutian region，Xinjiang.紅色圓點代表1970 年以來的7級以上強震

于田地區(qū)平均海拔4,000m，該地區(qū)人煙稀少，因此對該區(qū)域2次M7地震的研究相對較少。徐錫偉等(2011)對2008年于田地震的地表破裂帶進行了野外調(diào)查，發(fā)現(xiàn)此次地震的發(fā)震斷層為位于阿爾金斷裂西南尾端的張性構(gòu)造，該地震表明青藏高原的昆侖塊體、柴達木塊體、祁連塊體存在著E向滑移(Xuetal.，2013)?；贗nSAR資料和地震波資料，許多人對2008年于田地震破裂帶上的滑動分布進行了研究(Elliottetal.，2010; Furuyaetal.，2011；張國宏等，2011；劉智榮等，2012)。利用這些成果，不少人對這2次地震的動力學背景、應(yīng)力擾動等做出了初步研究，并對區(qū)域潛在地震危險性進行了判定(萬永革等，2010; 程惠紅等，2014; 任俊杰等，2014; 宋金等，2014; 周云等，2015)。這些結(jié)果大多基于彈性假設(shè)，沒有考慮中、下地殼黏彈性效應(yīng)的影響，由此得到的結(jié)果存在一定的局限性。本文采用分層黏彈性模型研究2008年和2014年2次于田地震所造成的區(qū)域應(yīng)力調(diào)整以及這2次地震的相互關(guān)系，得到的結(jié)果能為認識區(qū)域潛在地震危險性提供科學依據(jù)。

1　方法與模型

1.1計算方法

庫倫應(yīng)力變化ΔCFS可以根據(jù)如下公式計算：

其中，Δτ為剪應(yīng)力變化；μ′為視摩擦系數(shù)，通常取為0.4；Δσ為正應(yīng)力變化。Δτ定義當其方向與特定方向一致時為正；Δσ定義壓應(yīng)力為正(Stein，1999)。根據(jù)庫倫破裂準則，上述公式表明庫倫應(yīng)力變化增加的地方更加趨近于破裂，而庫倫應(yīng)力變化減少的地方相對更加安全(Freed，2005)。

前人提出了利用位錯模型計算半無限彈性空間中位錯所造成的應(yīng)力變化的具體表達式(Okada，1992)，并在研究單個大地震的同震應(yīng)力變化(McCloskeyetal.，2005; Parsonsetal.，2008; Todaetal.，2008；萬永革等，2010)以及一系列地震所造成的區(qū)域應(yīng)力擾動(Steinetal.，1997；Papadimitriouetal.，2004; 沈正康等，2004; 王輝等，2008)中得到了廣泛應(yīng)用。在此基礎(chǔ)上，考慮到中、下地殼-上地幔黏彈性松弛效應(yīng)的震后變形模型也在震后應(yīng)力調(diào)整研究中得到了應(yīng)用(沈正康等，2003; Freed，2005; Luoetal.，2010; Heetal.，2011; Shanetal.，2013；徐晶等，2013)。

本研究中采用半無限空間的彈性/非彈性成層模型(Wangetal.，2006)計算2008年和2014年2次于田地震的同震與震后應(yīng)力調(diào)整。與經(jīng)典彈性位錯模型(Okada，1992)相比，彈性/非彈性成層模型可以考慮中、下地殼—上地幔介質(zhì)流變所導致的上地殼應(yīng)力積累，進而研究與時間相關(guān)的震后變形過程。為了研究應(yīng)力調(diào)整與后續(xù)地震活動的關(guān)系，主要關(guān)注庫倫應(yīng)力變化量≥0.01MPa的區(qū)域范圍(Stein，1999)。

研究地震所造成的應(yīng)力擾動通常選擇下次地震的震源深度作為觀測深度?？紤]到地震定位深度的誤差，同時小震往往集中發(fā)生在1個深度范圍內(nèi)，我們采用最大庫倫應(yīng)力策略(the maximum-Coulomb-stress strategy)，分別計算每個節(jié)點上孕震層不同深度上的ΔCFS，并取其中的最大值作為該節(jié)點上的ΔCFS結(jié)果，由此得到的ΔCFS能夠更好地避免計算誤差，以解釋區(qū)域的小震活動。綜合前人結(jié)果，取于田地區(qū)上地殼孕震層的深度為0～15km(胥頤等，2006; Huangetal.，2011)。

1.2模型設(shè)置

對2008年和2014年2次于田地震的研究主要基于地震學和InSAR等的結(jié)果。震源機制結(jié)果表明，2008年于田地震主要為正斷兼走滑型地震，而2014年于田地震主要表現(xiàn)為走滑兼正斷型地震(表1)。野外調(diào)查結(jié)果表明，2008年于田地震的地表破裂帶呈SN—SN—NNE向展布，全長約31km(徐錫偉等，2011)。根據(jù)地表破裂類型及其組合特征，該次地震的地表破裂帶可以分為3個基本段，地表破裂帶測量到的最大左旋走滑位移為1.8m，最大垂直位移約2.0m。InSAR觀測的同震形變場分析結(jié)果表明，該地震的震源破裂帶長45km，破裂面上的最大滑動量超過4m，反演誤差<0.1m(Elliottetal.，2010)(圖2)。

表1　2008年和2014年2次于田地震的震源參數(shù)

Table1　Source parameters of the 2008 and 2014 Yutian earthquakes

事件時間經(jīng)度/(°)緯度/(°)走向/(°)傾向/(°)滑動角/(°)來源12008-03-2081.3735.4320352286GCMT82.4735.4921969292USGS22014-02-1282.5936.2724282-4GCMT

2014年2月12日的于田M7.3地震位于2008年地震宏觀震中東北方向約200km處。此次地震同樣發(fā)生在阿爾金斷裂帶南端。利用衛(wèi)星遙感資料，發(fā)現(xiàn)該地震的地表破裂帶可識別長度約為9km，破裂面上未見明顯的水平位錯(徐岳仁等，2015)。根據(jù)遠場P波波形和SH波波形資料反演得到的震源過程結(jié)果顯示，該次于田地震震級為MW6.9，震源深度10km，震源破裂面長度達到85km，地震滑動以左旋走滑為主并具有正傾滑分量，最大值為2.10m(周云等，2015)。

由于巖石圈變形在短期內(nèi)表現(xiàn)為彈性，長期變形更接近于流體，可以采用Maxwell體作為巖石圈變形的一階近似。地震學資料顯示，于田地區(qū)的莫霍面深度大約為62km(Huangetal.，2011)。與此同時，根據(jù)強震震后回彈估計的青藏高原中、下地殼-上地幔的黏滯系數(shù)≤1×1018Pa,s(Ryderetal.，2007；張晁軍等，2008; 孫玉軍等，2013)。綜合區(qū)域地震波速結(jié)構(gòu)(Mechieetal.，2012)，設(shè)置模型材料參數(shù)如表2。

表2　分層黏彈模型的模型參數(shù)

Table2　Parameters of layered viscoelastic model

深度/kmVp/km·s-1Vs/km·s-1ρ/kg·m-3η/Pas上地殼0～276.03.452760中、下地殼27～406.13.5528001×1018上地幔40～S607.14.1532001×1020

2　2008年于田地震對2014年于田地震的影響

2.1同震應(yīng)力擾動

為了方便與前人結(jié)果進行比較，首先根據(jù)2008年于田地震的震源機制解計算了均勻位錯造成的同震庫倫應(yīng)力變化。震源破裂帶的參數(shù)參考文獻(萬永革等，2010)，地震破裂面長度根據(jù)經(jīng)驗公式(Wellsetal.，1994)計算為57.5km，滑動量為1.88m，斷層寬度為25.6km。圖3a給出了深度為8km處最優(yōu)面方向上的ΔCFS結(jié)果。計算結(jié)果顯示，2008年于田地震導致破裂帶兩側(cè)斷層上的同震ΔCFS減少，地震破裂帶端點附近斷層上的同震ΔCFS增加。具體而言，于田地區(qū)的康西瓦斷裂帶、大紅柳灘斷裂帶、泉水溝斷裂帶和郭扎錯斷裂帶東段附近的同震ΔCFS減少超過0.01MPa，而阿爾金斷裂帶、郭扎錯斷裂帶西段和龍木措斷裂帶附近的同震ΔCFS增加超過了0.01MPa。雖然上述結(jié)果與前人結(jié)果較為一致(萬永革等，2010)，但是該結(jié)果與后續(xù)余震的空間分布差別較大。根據(jù)中國地震臺網(wǎng)統(tǒng)一地震目錄(http: ∥data.earthquake.cn/datashare/datashare_tyml_query.jsp)，選取了2009年1月1日至2013年12月31日的區(qū)域小震，發(fā)現(xiàn)2008年于田地震后的余震主要集中在破裂面兩側(cè)附近。統(tǒng)計結(jié)果顯示，在2008年于田地震同震ΔCFS量值>0.01MPa的空間范圍內(nèi)，只有大約19.2%的小震發(fā)生在同震ΔCFS增加的地方。由此可見，雖然根據(jù)經(jīng)驗公式和均勻位錯模型能夠比較快速地計算出較大空間范圍內(nèi)的區(qū)域應(yīng)力擾動情況，但是其與后續(xù)余震活動的空間分布差別較大，均勻位錯模型的結(jié)果具有一定的局限性。

結(jié)合地震波和InSAR地表形變觀測，Elliott 等(2010)給出了2008年于田地震的斷層滑動模型(圖2)。采用Elliot的斷層模型作為斷層位錯輸入，計算了2008年于田地震在0～15km深度內(nèi)各個深度的同震ΔCFS，然后選取不同深度ΔCFS的最大值作為區(qū)域整個孕震層中最大ΔCFS的空間分布(圖3b)。比較圖3a和圖3b可以看出，非均勻位錯模型計算的同震ΔCFS量值>0.01MPa的空間范圍略小于均勻模型結(jié)果。非均勻位錯模型計算的同震ΔCFS顯著影響區(qū)的空間分布形態(tài)與均勻模型的結(jié)果大致相同，都是斷層面兩側(cè)同震ΔCFS減少，斷層面2個端點處同震ΔCFS增加。然而，非均勻位錯模型計算的同震ΔCFS增加的區(qū)域面積明顯大于均勻位錯模型，特別是非均勻位錯模型計算的同震ΔCFS在接近地震破裂面的地區(qū)為正值，正好解釋了地震破裂面附近較頻繁的余震活動。同樣，將余震空間分布與同震ΔCFS空間分布進行比較，可以發(fā)現(xiàn)于田地區(qū)2009年以來的地震有53%發(fā)生在同震ΔCFS增加超過0.01MPa的區(qū)域，2008年于田地震的發(fā)生顯著影響了破裂帶附近區(qū)域的地震活動。

圖3　均勻位錯模型和非均勻位錯模型計算的2008年于田地震同震ΔCFS空間分布Fig. 3　Coseismic ΔCFS produced by the 2008 Yutian earthquake.黃色圓點為2009年1月1日至2013年12月31日的小震活動，黑色線條為區(qū)域主要活動斷裂帶，綠色線條為地震破裂帶；a 均勻位錯模型計算的最優(yōu)面上的同震ΔCFS，b 非均勻位錯模型計算的最優(yōu)面上的同震ΔCFS

圖4　非均勻位錯模型計算的2008年于田地震同震ΔCFS 誤差空間分布Fig. 4　Error of coseismic ΔCFS produced by uneven rupture model of the 2008 Yutian earthquake.Green lines show the seismic rupture.

為了考察同震ΔCFS計算的可靠性，進一步根據(jù)Elliot模型的誤差計算了2008年于田地震同震ΔCFS的誤差分布(圖4)。圖4 的結(jié)果顯示，同震ΔCFS的誤差分布與同震ΔCFS的空間分布模式大致一致，同震ΔCFS誤差最大的區(qū)域集中在地震破裂面的兩側(cè)附近，地震破裂帶南段上較大的滑動誤差導致其附近相對較大的同震ΔCFS的誤差。另外，由于地震破裂面滑動誤差主要為0.05～0.10m，遠小于反演的斷層滑動量，同震ΔCFS誤差擾動的量值也不到同震ΔCFS量值的10%。

比較均勻位錯模型和非均勻位錯模型的結(jié)果可以發(fā)現(xiàn)，非均勻斷層滑動模型計算得到的同震ΔCFS分布比簡單均勻模型的結(jié)果更好地解釋實際觀測到的余震空間分布。另外，無論是均勻模型還是非均勻模型的同震ΔCFS計算結(jié)果均表明，2008年于田地震在2014年于田地震震源區(qū)附近所造成的同震ΔCFS變化遠<0.01MPa，2008年于田地震所產(chǎn)生的同震應(yīng)力擾動對2014年于田地震的發(fā)生影響不顯著。

2.2黏彈性應(yīng)力擾動

雖然同震ΔCFS的結(jié)果顯示2008年于田地震對2014年于田地震的發(fā)生沒有顯著影響，但是由于青藏高原地區(qū)的中、下地殼相對較軟(Ryderetal.，2007；張晁軍等，2008; 孫玉軍等，2013)，柔軟中、下地殼的黏性變形會導致上地殼的應(yīng)力集中。因此，2008年于田地震發(fā)生6a后的震后變形影響范圍有可能擴展到2014年于田地震的震源區(qū)附近。

圖5　2008年于田地震造成5a后的ΔCFS空間分布Fig. 5　Postseismic ΔCFS produced by the 2008 Yutian earthquake in 2014.黃色圓點為2009年1月1日至2013年12月31日的小震活動，綠色線條為地震破裂帶；a分層黏彈模型計算的最優(yōu)面上的震后ΔCFS，b投影到區(qū)域斷層分段上的震后ΔCFS

采用1.2節(jié)中的模型和參數(shù)設(shè)置，計算了2008年于田地震6a后，中下地殼黏彈性松弛效應(yīng)所造成的區(qū)域應(yīng)力擾動(圖5)。圖5a給出了2008年于田地震6a后附近區(qū)域整個孕震層中最大ΔCFS的空間分布。與同震ΔCFS的空間分布相比，2008年于田地震6a后變形所造成的應(yīng)力擾動空間范圍進一步擴大，應(yīng)力擾動>0.01MPa的區(qū)域擴展到距離2008年于田地震地震破裂面近200km的地方。對比2008年于田地震后的ΔCFS空間分布和2009—2013年間的小震分布，可以發(fā)現(xiàn)，在震后應(yīng)力擾動>0.01MPa的區(qū)域內(nèi)發(fā)生的小震大約74%集中在ΔCFS增加的地區(qū)，只有約26%的小震集中在ΔCFS減少的地區(qū)。

根據(jù)計算的震后應(yīng)力場和前人提供的斷層幾何與運動學參數(shù)(萬永革等，2010)，進一步計算了震后ΔCFS影響顯著區(qū)域內(nèi)，主要活動斷裂帶分段上的震后ΔCFS(圖5b)。計算結(jié)果顯示，2008年于田地震導致康西瓦斷裂帶東段上的震后ΔCFS增加最大，約為0.07MPa。普魯斷裂帶上從東到西的震后ΔCFS增加0.01～0.05MPa不等，貢嘎斷裂帶上的震后ΔCFS減少0.13～0.59MPa。郭扎錯斷裂帶東段上的震后ΔCFS減少0.11MPa，而其西段上的震后ΔCFS增加0.04～0.10MPa。龍木措斷裂帶上震后ΔCFS增加0.03MPa。泉水溝和大紅柳灘斷裂帶上的震后ΔCFS減少0.02～0.04MPa。值得注意的是，2008年于田地震的震后應(yīng)力調(diào)整使得2014年于田地震破裂面上的震后ΔCFS增加0.011MPa。

對比2008年于田地震同震ΔCFS和震后ΔCFS的空間分布圖像可以看出，雖然2008年的于田地震所造成的同震ΔCFS對2014年于田地震的影響不顯著，但是震后應(yīng)力松弛效應(yīng)導致的應(yīng)力轉(zhuǎn)移有利于6a后2014年于田地震的發(fā)生。

圖6　2014年于田地震造成的震后ΔCFS空間分布Fig. 6　Coseismic and postseismic ΔCFS produced by the 2014 Yutian earthquake.黃色圓點為2014年2月12日至2014年12月31日的小震活動，綠色線條為地震破裂帶；a分層黏彈模型計算的最優(yōu)面上的同震ΔCFS，b投影到區(qū)域斷層分段上的同震ΔCFS，c分層黏彈模型計算的5a后的最優(yōu)面上的震后ΔCFS，d投影到區(qū)域斷層分段上的震后ΔCFS

3　2014年于田地震所造成的應(yīng)力擾動

2014年于田地震的發(fā)生可能再次擾動了區(qū)域的應(yīng)力應(yīng)變分布，從而影響后續(xù)地震活動。根據(jù)地震破裂過程給出的破裂面滑動分布結(jié)果(周云等，2015)，同樣計算了2014年于田地震所造成的同震ΔCFS和震后ΔCFS應(yīng)力擾動(圖6)。圖6a給出了2014年于田地震所造成的同震ΔCFS，圖6b給出了投影到區(qū)域主要斷裂分段上的同震ΔCFS，圖6c是5a以后黏彈性松弛導致的震后ΔCFS，圖6d則給出了5a以后主要斷裂分段上的震后ΔCFS。

計算結(jié)果顯示，2014年于田地震的發(fā)生主要使得阿爾金斷裂帶西南段的同震ΔCFS增加，與地震破裂帶平行的附近區(qū)域同震ΔCFS減少區(qū)與增加區(qū)相互間隔。從余震分布的情況來看，在同震ΔCFS變化較顯著的區(qū)域內(nèi)，大約76.1%的余震發(fā)生在同震ΔCFS增加的地區(qū)。這些余震主要集中在地震破裂帶附近，并且呈現(xiàn)出南、北2個余震集中區(qū)。與反演的地震破裂面滑動分布結(jié)果相對比，可以看出余震主要集中在地震破裂面滑動分布相對較小的區(qū)域。也就是說，2014年于田地震后的余震活動可能補償了主震破裂面上的同震滑動分布。這與2008年汶川M8.0地震的主震滑動與余震活動分布模式類似(陳九輝等，2009)。

同樣地，計算了2014年于田地震在區(qū)域主要活動斷裂帶分段上造成的同震ΔCFS。結(jié)果表明，2014年于田地震使得康西瓦斷裂帶東段上的同震ΔCFS增加0.16～0.85MPa，普魯斷裂帶上的同震ΔCFS減少0.42MPa，貢嘎斷裂帶上的同震ΔCFS減少0.42MPa?？缕簲嗔褞贤穰FS減少0.03MPa，郭扎錯斷裂帶上的同震ΔCFS主要表現(xiàn)為減少，大約為0.04～0.06MPa。與前人計算類似的結(jié)果證明了本文結(jié)果的可靠性(周云等，2015)。

震后黏彈性松弛將使2014年于田地震的震后變形影響范圍進一步擴大。預(yù)計2014年于田地震后5a，黏彈性松弛效應(yīng)將使震后ΔCFS量值>0.01MPa的范圍大約達到顯著同震形變范圍的1.5倍，但是其形態(tài)大致保持不變?？滴魍邤嗔褞|段、中段上的震后ΔCFS變化增強，達到0.35～0.98MPa，貢嘎斷裂帶上的震后ΔCFS減少0.03～0.04MPa，普魯斷裂帶上的震后ΔCFS減少0.21～0.23MPa?？缕簲嗔褞险鸷螃FS減少0.03～0.05MPa，郭扎錯斷裂帶上的同震ΔCFS主要表現(xiàn)為減少，大約為0.03～0.08MPa。2014年于田地震震后變形的影響還將擴展到其南邊的一些斷裂帶上。

4　討論

由于具有明確的物理含義，利用位錯模型研究震后ΔCFS并判斷區(qū)域余震趨勢對于潛在地震危險性的判定有重要意義。前人利用相對較簡單的均勻彈性位錯模型計算2008年和2014年2次于田地震的同震ΔCFS，得到的結(jié)果能夠?qū)τ诤罄m(xù)余震的快速判定提供指導( 萬永革等，2010；宋金等，2014)。然而，均勻彈性位錯模型結(jié)果在近場與實際的余震活動符合程度較差。即使采用相對完善的震源破裂模型作為輸入能夠得到更加符合實際情況的結(jié)果，但是彈性模型并不能考察震后松弛所導致的應(yīng)力轉(zhuǎn)移(周云等，2015)。

除了同震ΔCFS與余震分布的關(guān)系研究之外，震后ΔCFS對后續(xù)地震的影響也很顯著。同一個斷裂帶上或者鄰近地區(qū)的強震活動往往表現(xiàn)為叢集活動，2次大地震的時間間隔也在數(shù)年至數(shù)百年的尺度，如此長時間內(nèi)的地殼流變效應(yīng)不可忽略。本文的研究也表明，2008年于田地震的同震ΔCFS對2014年于田地震的發(fā)生影響不大，如果僅憑同震ΔCFS對區(qū)域潛在地震危險性進行判定可能會帶來誤導。事實上，震后黏彈性松弛所導致的應(yīng)力轉(zhuǎn)移使得2008年于田地震在2014年于田地震震源區(qū)產(chǎn)生超過了0.01MPa的ΔCFS，從而2008年于田地震的震后變形有利于2014年于田地震的發(fā)生。雖然黏滯系數(shù)的選取會影響震后變形的計算結(jié)果，但是下地殼黏滯系數(shù)計算結(jié)果的誤差往往在2～3個數(shù)量級(張晁軍等，2008)，數(shù)值試驗表明較低的下地殼黏滯系數(shù)更有利于震后變形的擴散。無論如何，震后黏性松弛所導致的應(yīng)力轉(zhuǎn)移需要特別引起注意。

在前人結(jié)果的基礎(chǔ)上，本文采用多種方法考慮ΔCFS的計算誤差。首先，選取整個孕震層不同深度上ΔCFS的最大值作為最終的ΔCFS，避免了后續(xù)地震震源深度誤差對ΔCFS計算結(jié)果的影響，最大ΔCFS與余震分布的對比也更加合理。其次，根據(jù)位錯模型的誤差給出了同震ΔCFS的誤差分布，結(jié)果表明同震ΔCFS的誤差主要集中在斷裂帶兩側(cè)附近，誤差的量值小于ΔCFS變化絕對值的10%。由此，本文得到的區(qū)域最優(yōu)面方向的ΔCFS空間分布圖像相對可靠。最后，前人的研究表明，對ΔCFS計算影響最大的因素是接收斷層的傾角誤差，其次是斷層模型、接收斷層的走向和滑動角方向的誤差(Wangetal.，2014)。受到于田地區(qū)活動構(gòu)造研究局限性的限制，該區(qū)域各個接收斷裂帶上的同震ΔCFS和震后ΔCFS只能提供相對的定性參考，其結(jié)果還需要根據(jù)后續(xù)斷層幾何學與運動學的研究進一步完善。

圖7　2008年于田地震的震后變形與和2014年于田地震同震影響所造成的區(qū)域應(yīng)力擾動Fig. 7　Regional comprehensive stress disturbance produced by 2008 and 2014 Yutian earthquakes.

5　結(jié)論

本文采用分層黏彈性模型詳細分析了2008年和2014年2次M7于田地震所造成的同震和震后應(yīng)力擾動，并研究了2008年于田地震對2014年于田地震的影響，得到了如下結(jié)論：

(1)2008年于田地震的發(fā)生主要使得破裂面附近區(qū)域的同震ΔCFS增加，絕大部分余震發(fā)生在同震ΔCFS增加的地區(qū)。此次地震使得阿爾金斷裂帶西端的ΔCFS增加，但是其對2014年于田地震的影響并不顯著；

(2)中、下地殼的黏性松弛效應(yīng)導致2008年于田地震的震后應(yīng)力調(diào)整，有利于2014年于田地震的發(fā)生；

(3)2014年于田地震的發(fā)生進一步增加了阿爾金斷裂帶西段的ΔCFS，考慮到該斷裂帶上較高的活動速率和過去較長時間的地震平靜，該斷裂帶上的潛在地震危險性值得進一步關(guān)注。

致謝感謝審稿專家提出的寶貴意見。Elliot博士提供了2008年于田地震的震源破裂參數(shù)，王衛(wèi)民博士提供了2014年于田地震的震源破裂參數(shù)，汪榮江博士提供了計算軟件，在此一并致以謝忱!

陳九輝，劉啟元，李順成，等. 2009. 汶川MS8.0地震余震序列重新定位及其地震構(gòu)造研究 [J]. 地球物理學報，52(2): 390—397.

CHEN Jiu-hui，LIU Qi-yuan，LI Shun-cheng，etal. 2009. Seismotectonic study by relocation of the WenchuanMS8.0 earthquake sequence [J]. Chinese Journal of Geophysics，52(2): 390—397(in Chinese).

程惠紅，龐亞瑾，董培育，等. 2014. 于田2008年和2014年2次MS7.3地震孕育的應(yīng)力環(huán)境 [J]. 地球物理學報，57(10): 3238—3246. doi: 10.6038/cjg20141012.

CHEN Hui-hong，PANG Ya-jin，DONG Pei-yu，etal. 2014. Analysis of the stress environment of the 2008 and 2014 YutianMS7.3 earthquake [J]. Chinese Journal of Geophysics，57(10): 3238—3246(in Chinese).

劉智榮，黃靜宜，白相東，等. 2012. 基于梯度下降法反演于田MW7.1地震同震滑動分布 [J]. 大地測量與地球動力學，32(6): 8—12.

LIU Zhi-rong，HUANG Jing-yi，BAI Xiang-dong，etal. 2012. Coseismic slip distribution of YutianMW7.1 earthquake based on steepest descend method [J]. Journal of Geodesy and Geodynamics，32(6): 8—12(in Chinese).

馬宏生，張國民，劉杰，等. 2006. 中國大陸活動地塊邊界帶強震活動特征的研究 [J]. 地震地質(zhì)，28(1): 48— 60.

MA Hong-sheng，ZHANG Guo-min，LIU Jie，etal. 2006. Research on the characteristics of large earthquake activity on the active tectonic boundaries in Chinese mainland [J]. Seismology and Geology，28(1): 48— 60(in Chinese).

任俊杰，王信國，程佳. 2014. 2014年于田7.3級地震的發(fā)震構(gòu)造及動力學背景的初步分析 [J]. 震災(zāi)防御技術(shù)，9(1): 74—82.

REN Jun-jie，WANG Xin-guo，CHENG Jia. 2014. Preliminary analysis of seismogenic fault and geodynamic background of 2014 YutianMS7.3 earthquake [J]. Technology for Earthquake Disaster Prevention，9(1): 74—82(in Chinese).

沈正康，萬永革，甘衛(wèi)軍，等. 2003. 東昆侖活動斷裂帶大地震之間的黏彈性應(yīng)力觸發(fā)研究 [J]. 地球物理學報，46(6): 786—795.

SHEN Zheng-kang，WAN Yong-ge，GAN Wei-jun，etal. 2003. Viscoelastic triggering among large earthquakes along the east Kunlun fault system [J]. Chinese Journal of Geophysics，46(6): 786—795(in Chinese).

沈正康，萬永革，甘衛(wèi)軍，等. 2004. 華北地區(qū)700年來地殼應(yīng)力場演化與地震的關(guān)系研究 [J]. 中國地震，20(3): 211—228.

SHEN Zheng-kang，WAN Yong-ge，GAN Wei-jun，etal. 2004. Crustal stress evolution of the last 700 years in North China and earthquake occurrence [J]. Earthquake Research in China，20(3): 211—228(in Chinese).

宋金，周龍泉. 2014. 2014 年于田MS7.3 地震產(chǎn)生的靜態(tài)庫侖應(yīng)力變化及對周邊斷層的影響 [J]. 中國地震，30(2): 168—177.

SONG Jin，ZHOU Long-quan. 2014. The static stress triggering effects related with the YutianMS7.3 earthquake [J]. Earthquake Research in China，30(2): 168—177(in Chinese).

孫玉軍，董樹文，范桃園，等. 2013. 中國大陸及鄰區(qū)巖石圈三維流變結(jié)構(gòu) [J]. 地球物理學報，56(9): 2936—2946. doi: 10.6038/cjg20130908.

SUN Yu-jun，DONG Shu-wen，F(xiàn)AN Tao-yuan，etal. 2013. 3D rheological structure of the continental lithosphere beneath China and adjacent regions [J]. Chinese Journal of Geophysics，56(9): 2936—2946(in Chinese).

萬永革，沈正康，盛書中，等. 2010. 2008年新疆于田7.3級地震對周圍斷層的影響及其正斷層機制的區(qū)域構(gòu)造解釋 [J]. 地球物理學報，53(2): 280—289. doi: 10.3969/j.issn.0001-5733.2010.02.006.

WAN Yong-ge，SHEN Zheng-kang，SHENG Shu-zhong，etal. 2010. The mechanical effects of the 2008MS7.3 Yutian，Xinjiang earthquake on the neighboring faults and its tectonic origin of normal faulting mechanism [J]. Chinese Journal of Geophysics，53(2): 280—289(in Chinese).

王輝，劉杰，石耀霖，等. 2008. 鮮水河斷裂帶強震相互作用的動力學模擬研究 [J]. 中國科學(D輯) 38(7): 808—818.

WANG Hui，LIU Jie，SHI Yao-lin，etal. 2008. Dynamic simulation of interactions between major earthquakes on the Xianshuihe fault zone [J]. Science in China(Ser D): Earth Sciences，51(10): 1388—1400.

徐晶，邵志剛，馬宏生，等. 2013. 鮮水河斷裂帶庫侖應(yīng)力演化與強震間關(guān)系 [J]. 地球物理學報，56(4): 1146—1158. doi: 10.6038/cjg20130410.

XU Jing，SHAO Zhi-gang，MA Hong-sheng，etal. 2013. Evolution of Coulomb stress and stress interaction among strong earthquakes along the Xianshuihe fault zone [J]. Chinese Journal of Geophysics，56(4): 1146—1158(in Chinese).

徐錫偉，譚錫斌，吳國棟，等. 2011. 2008年于田MS7.3地震地表破裂帶特征及其構(gòu)造屬性討論 [J]. 地震地質(zhì)，33(2): 462— 471. doi: 10.3969/j.issn.0253-4967.2011.02.019.

XU Xi-wei，TAN Xi-bin，WU Guo-dong，etal. 2011. Surface rupture features of the 2008 YutianMS7.3 earthquake and its tectonic nature [J]. Seismology and Geology，33(2): 462— 471(in Chinese).

胥頤，劉建華，劉福田，等. 2006. 天山-帕米爾結(jié)合帶的地殼速度結(jié)構(gòu)及地震活動研究 [J]. 地球物理學報，49(2): 469— 476.

XU Yi，LIU Jian-hua，LIU Fu-tian，etal. 2006. Crustal velocity structure and seismic activity in the Tianshan-Pamir conjunctive zone [J]. Chinese Journal of Geophysics，49(2): 469— 476(in Chinese).

徐岳仁，陳立澤，申旭輝，等. 2015. 基于GF-1衛(wèi)星影像解譯2014年新疆于田MS7.3地震同震地表破裂帶 [J]. 地震，35(2): 61—71.

XU Yue-ren，CHEN Li-ze，SHEN Xu-hui，etal. 2015. Interpretating coseismic surface rupture zone of the 2014 YutianMS7.3 earthquake using GF-1 satellite images [J]. Earthquake，35(2): 61—71(in Chinese).

張晁軍，曹建玲，石耀霖. 2008. 從震后形變探討青藏高原下地殼黏滯系數(shù) [J]. 中國科學(D輯)，38(10): 1250—1257.

ZHANG Chao-jun，CAO Jian-ling，SHI Yao-lin. 2009. Studying the viscosity of lower crust of Qinghai-Tibet plateau according to post-seismic deformation [J]. Science in China(Ser D)，52(3): 411— 419.

張國宏，屈春燕，單新建，等. 2011. 2008年MS7.1于田地震InSAR同震形變場及其震源滑動反演 [J]. 地球物理學報，54(11): 2753—2760. doi: 10.3969/j.issn.0001-5733.2011.11.005.

ZHANG Guo-hong，QU Chun-yan，SHAN Xin-jianetal. 2011. The coseismic InSAR measurements of 2008 Yutian earthquake and its inversion for source parameters [J]. Chinese Journal of Geophysics，54(11): 2753—2760(in Chinese).

周云，王衛(wèi)民，熊林，等. 2015. 2014年2月12日MW6.9于田地震震源破裂過程及對周圍斷層的應(yīng)力影響 [J]. 地球物理學報，58(1): 184—193. doi: 10.6038/cjg20150116.

ZHOU Yun，WANG Wei-min，XIONG Lin，etal. 2015. Rupture process of 12 February 2014，YutianMW6.9 earthquake and stress change on nearby faults [J]. Chinese Journal of Geophysics，58(1): 184—193(in Chinese).

Elliott J R，Walters R J，England P C，etal. 2010. Extension on the Tibetan plateau: Recent normal faulting measured by InSAR and body wave seismology [J]. Geophysical Journal International，183: 503—535. doi: 10.1111/j.1365-246X.2010.047,54.x.

Freed A M. 2005. Earthquake triggering by static，dynamic，and postseismic stress transfer [J]. Annual Review of Earth and Planetary Sciences，33: 335—367. doi: 10.1146/annurev.earth.33.092203.122,505.

Furuya M，Yasuda T. 2011. The 2008 Yutian normal faulting earthquake(MW7.1)，NW Tibet: Non-planar fault modeling and implications for the Karakax Fault [J]. Tectonophysics，511: 125—133. doi: 10.1016/j.tecto.2011.09.003.

He J，Vernant P，Chery J，etal. 2013. Nailing down the slip rate of the Altyn Tagh Fault [J]. Geophysical Research Letters，40: 5382—5386. doi: 10.1002/2013GL057497.

He J，Xia W，Lu S，etal. 2011. Three-dimensional finite element modeling of stress evolution around the Xiaojiang fault system in the southeastern Tibetan plateau during the past～500 years [J]. Tectonophysics，507: 70—85. doi: 10.1016/j.tecto.2011.05.009.

Huang G C D，Roecker S W，Levin V. 2011. Lower-crustal earthquakes in the West Kunlun range [J]. Geophysical Research Letters，38(L01314). doi: 10.1029/2010GL045893.

Luo G，Liu M. 2010. Stress evolution and fault interactions before and after the 2008 great Wenchuan earthquake [J]. Tectonophysics，491: 127—140. doi: 10.1016/j.tecto.2009.12.019.

McCloskey J，Nalbant S S，Steacy S. 2005. Earthquake risk from co-seismic stress [J]. Nature，434(3): 291.

Mechie J，Yuan X，Schurr B，etal. 2012. Crustal and uppermost mantle velocity structure along a profile across the Pamir and southern Tien Shan as derived from project TIPAGE wide-angle seismic data [J]. Geophysical Journal International，188: 385— 407. doi: 10.1111/j.1365-246X.2011.052,78.x.

Okada Y. 1992. Internal deformation due to shear and tensile faults in a half-space [J]. Bulletin of the Seismological Society of America，82(2): 1018—1040.

Papadimitriou E，Wen X，Karakostas V，etal. 2004. Earthquake triggering along the Xianshuihe fault zone of western Sichuan，China [J]. Pure and Applied Geophysics，161: 1683—1707.

Parsons T，Ji C，Kirby E. 2008. Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan Basin [J]. Nature，454(7203): 509—510. doi: 10.1038/nature07177.

Ryder I，Parsons B，Wright T J，etal. 2007. Post-seismic motion following the 1997 Manyi(Tibet)earthquake: InSAR observations and modelling [J]. Geophysical Journal International，169(3): 1009—1027.

Shan B，Xiong X，Wang R，etal. 2013. Coulomb stress evolution along Xianshuihe-Xiaojiang fault system since 1713 and its interaction with Wenchuan earthquake，May 12，2008 [J]. Earth and Planetary Science Letters，377-378: 199—210. doi: 10.1016/j.tifl.2013.06.044.

Stein R S. 1999. The role of stress transfer in earthquake occurrence [J]. Nature，402(6762): 605— 609.

Stein R S，Barka A A，Dieterich J H. 1997. Progressive failure on the North Anatolian Fault since 1939 by earthquake stress triggering [J]. Geophysical Journal International，128(3): 594— 604.

Toda S，Lin J，Meghraoui M，etal. 2008. 12 May 2008M=7.9 Wenchuan，China，earthquake calculated to increase failure stress and seismicity rate on three major fault systems [J]. Geophysical Research Letters，35(L17305): doi: 10.1029/2008GL034903.

Wang H，Liu M，Cao J，etal. 2011. Slip rates and seismic moment deficits on major active faults in mainland China [J]. Journal of Geophysical Research-Solid Earth，116(B02405). doi: 10.1029/2010JB007821.

Wang J，Xu C，F(xiàn)reymueller J T，etal. 2014. Sensitivity of Coulomb stress change to the parameters of the Coulomb failure model: A case study using the 2008MW7.9 Wenchuan earthquake [J]. Journal of Geophysical Research-Solid Earth，119(4): 3371—3392. doi: 10.1002/2012JB009860.

Wang R，Lorenzo-Martin F，Roth F. 2006. PSGRN/PSCMP-A new code for calculating co- and post-seismic deformation，geoid and gravity changes based on the viscoelastic-gravitational dislocation theory [J]. Computers & Geosciences，32(4): 527—541. doi: 10.1016/j.cageo.2005.08.006.

Wells D L，Coppersmith K J. 1994. New empirical relationships among magnitude，rupture length，rupture width，rupture area，and surface displacement [J]. Bulletin of the Seismological Society of America，84(4): 974—1002.

Xu X，Tan X，Yu G，etal. 2013. Normal- and oblique-slip of the 2008 Yutian earthquake-Evidence for eastward block motion，northern Tibetan plateau [J]. Tectonophysics，584: 152—165. doi: 10.1016/j.tecto.2012.08.007.

Zhu S，Miao M. 2015. How did the 2013 Lushan earthquake(MS=7.0)trigger its aftershocks？Insights from static Coulom stress change calculations [J]. Pure and Applied Geophysics，172(10): 2481—2495. doi: 10.100,7/s00024-015-1064-3.

VISCOELASTIC STRESS TRANSFER BETWEEN 2008 AND 2014 YUTIANM7 EARTHQUAKES，XINJIANG

WANG Hui1，2)CAO Jian-ling1，2)HONG Shun-ying2)XU Yue-ren2)JING Feng2)

1)KeyLaboratoryofEarthquakePrediction，ChinaEarthquakeAdministration，Beijing100036，China2)InstituteofEarthquakeScience，ChinaEarthquakeAdministration，Beijing100036，China

Two earthquakes with magnitude larger than 7.0 occurred in 2008 and 2014 on the southwestern end of the Altyn Tagh Fault，which is located in the northwestern borderland of Tibetan plateau. Occurrences of these two earthquakes provide important insights into regional geodynamics and potential seismic risk. Layered viscoelastic model is employed in the paper to study the interaction between these two events. We find that most of aftershocks were triggered by coseismic stress produced by the 2008 Yutian earthquake，and the effect of this earthquake is insignificant on the occurrence of the 2014 Yutian earthquake. However，stress transfer by viscoelastic relaxation of postseismic deformation is in favor of occurrence of the 2014 Yutian earthquake. The coseismic and postseismic stress transfer produced by the 2014 Yutian earthquake leads to stress increasing on the western segment of the Altyn Tagh Fault. Since the occurrence time of the last major earthquake on the western segment of the Altyn Tagh Fault is tens of years ago，it should have accumulated large moment deficit on the fault segment. The Altyn Tagh Fault should be considered as a fault with high potential seismic risk.

Yutian earthquake，layered viscoelastic model，ΔCFS，western segment of the Altyn Tagh Fault，potential seismic risk

10.3969/j.issn.0253- 4967.2016.03.011

2015-03-30收稿，2016-04-15改回。

中國地震局地震預(yù)測研究所基本科研業(yè)務(wù)專項(2014IES0102)與國家自然科學基金(41104058，41104057)共同資助。

P315.72+7

A

0253-4967(2016)03-0646-14

王輝，男，1976年生，研究員，主要研究方向為地殼形變與地球動力學，電話： 010-88015551，E-mail: wanghui500@gmail.com。