2017 年伊朗MW7.3 地震的震后余滑分布及其對2018 年MW6.0 地震的觸發(fā)影響＊

蔣子琴 楊瑩輝 陳 強徐 倩 徐 浪 黃小梅

1） 中國成都 611756 西南交通大學地球科學與環(huán)境工程學院

2） 中國成都 610059 成都理工大學地質(zhì)災害防治與地質(zhì)環(huán)境保護國家重點實驗室

引言

2017 年11 月12 日（世界時）伊朗邊境的薩波爾扎哈布鎮(zhèn)發(fā)生MW7.3 地震，該地震為二十一世紀以來在扎格羅斯山脈發(fā)生的震級最大的一次地震（Taymazet al，2018），地震震中位于（45.96°E，34.91°N）（USGS，2017）。2018 年8 月25 日該地區(qū)的賈萬魯?shù)驴h發(fā)生了一次MW6.0 地震（USGS，2018），震中位于（46.24°E，34.61°N）。兩次地震均位于扎格羅斯造山帶的高扎格羅斯斷層附近，震中僅相距約42 km，發(fā)震時間相隔不到一年 （圖1）。扎格羅斯山脈由亞歐板塊與阿拉伯板塊長期相互碰撞擠壓而形成，該區(qū)域孕育了豐富多樣的斷層系統(tǒng)，強烈的板塊活動造就了該區(qū)域頻繁的地震活動（Jahaniet al，2009；Hatzfeld，Molnar，2010）。

薩波爾扎哈布地震發(fā)生后，一些研究人員（楊百存等，2018；Yanget al，2018b； Heet al，2019； Yanget al，2019）基于合成孔徑雷達（synthetic aperture radar，縮寫為SAR）數(shù)據(jù)對該地震的地表形變及發(fā)震斷層展開了一系列研究。上述研究的結(jié)果表明：薩波爾扎哈布地震由近南北走向的低傾角斷層破裂導致，同震運動以逆沖兼走滑為主，且強震發(fā)生后該震區(qū)震后形變顯著。但是，關(guān)于發(fā)震斷層的震后運動特征未見詳細分析。而震后斷層活動是揭示主震引起的應力變化和余震觸發(fā)機制的關(guān)鍵信息，因此本文擬收集覆蓋研究區(qū)域的Sentinel-1 衛(wèi)星升、降軌SAR 數(shù)據(jù)，利用短基線集（small baseline subset，縮寫為SBAS）干涉合成孔徑雷達(interferometric synthetic aperture rada，縮寫為InSAR）技術(shù)（Moraet al，2003）提取薩波爾扎哈布地震后283 天的地表時序形變，以期揭示該地震同震斷層的運動性質(zhì)，解析其震后形變機制，為進一步討論震后效應對后續(xù)地震的觸發(fā)作用提供依據(jù)。

地震發(fā)生后會引起周圍區(qū)域應力的增加或卸載，而這可能會觸發(fā)或延緩后續(xù)地震的發(fā)生，例如1990 年菲律賓MW7.7 地震被認為觸發(fā)了其后的MW6.0 和MW6.5 強余震事件（李健等，2017），而2010 年MS7.1 玉樹地震并未觸發(fā)2011 年MS5.2 囊謙地震，囊謙地震更多應考慮背景地震活動性的影響（單斌等，2012）。萬永革等（2009）計算了汶川大地震在周圍斷層上的庫侖破裂應力變化，其結(jié)果對該地區(qū)的地震危險性分析具有重要的參考價值。有效地判定前震對后續(xù)地震的觸發(fā)效應對于認識某區(qū)域地震活動規(guī)律、評價未來地震的危險性具有重要作用，因此本文擬根據(jù)彈性位錯理論（Okada，1985）分別反演獲得薩波爾扎哈布地震震后余滑模型和賈萬魯?shù)碌卣鸬耐饠鄬踊瑒?，并通過計算薩波爾扎哈布地震的同震及震后余滑對賈萬魯?shù)碌卣鸢l(fā)震斷層的庫侖應力作用，探討兩次地震間的觸發(fā)關(guān)系，以期為深入探究扎格羅斯造山帶的地震活動特性提供參考。

1 薩波爾扎哈布地震的震后形變和余滑反演

1.1 數(shù)據(jù)獲取與處理

本文觀測數(shù)據(jù)為Sentinel-1 衛(wèi)星的C 波段合成孔徑雷達影像（European Space Agency，2014），其數(shù)據(jù)獲取周期短、分辨率高等優(yōu)點有利于短期時序監(jiān)測。針對薩波爾扎哈布地震的震后形變，本文選取覆蓋研究區(qū)域的12 景升軌影像和13 景降軌影像，利用短基線集干涉測量技術(shù)獲取震后283 天的時序形變。Sentinel-1A 影像的覆蓋范圍如圖1 所示，表1 列出了Sentinel-1A 影像的詳細參數(shù)。

表1 Sentinel-1A 衛(wèi)星的SAR 影像參數(shù)Table 1 SAR image parameters of Sentinel-1A satellite

圖1 2017 年和 2018 年兩次伊朗地震震源區(qū)的構(gòu)造背景及余震分布Fig. 1 Tectonic settings of the source regions of the two Iran earthquakes in 2017 and 2018 and aftershock distribution

為減少時空失相干的影響，升、降軌SAR 影像干涉圖集的時間基線均設置為150 天，垂直基線均設置為150 m，最終獲得27 組升軌干涉對和29 組降軌干涉對。SAR 影像處理過程中，引入精密定軌星歷（precise orbit ephemerides，縮寫為POD）數(shù)據(jù)以減少軌道誤差影響，使用空間分辨率約90 m 的數(shù)字高程模型（digital elevation model，縮寫為DEM）扣除地形相位影響，并采用最小費用流（minimum cost flow，縮寫為MCF）方法進行干涉相位解纏，繼而通過地理編碼和形變解算得到時序結(jié)果。

地殼孔隙介質(zhì)的彈性回彈、同震斷層面上的余滑以及下地殼與上地幔間的黏彈性松弛為三種主要的震后形變效應（Maroneet al，1991；Peltzeret al，1998；Jónssonet al，2003；Hsuet al，2006）。Feng 等（2018）關(guān)于主震后一個月內(nèi)累計余滑的研究表明，孔隙彈性回彈作用對薩波爾扎哈布地震的影響較小，而黏彈性松弛效應主要對遠場區(qū)域產(chǎn)生影響，且其效果在震后長時間跨度下更加顯著（Zhaoet al，2017；賀克鋒等，2019）。由于本研究中地表監(jiān)測時段較短，且主要關(guān)注近場的震后形變，因此可忽略黏彈性松弛效應，所監(jiān)測到的地表形變僅由余滑效應所致，將被用于反演薩波爾扎哈布地震的震后余滑分布。

為了提高反演效率，采用四叉樹方法（Lohman，Simons，2005）分別對升、降軌形變場進行降采樣。盡管震后余滑效應可能持續(xù)數(shù)月至數(shù)年時間，但相對于漫長的地震周期而言，這一過程仍是短暫的，因此可以采用與同震斷層相同的方法進行反演獲取余滑模型。本文使用“兩步法”，以Yang 等（2018b）研究中所估計的同震斷層幾何模型為初始斷層參數(shù)，首先假設在斷層均勻滑動的情形下，利用模擬退火算法搜索最優(yōu)的斷層幾何參數(shù)，然后將斷層離散為尺寸較小的矩形源，通過線性反演獲得斷層滑動的精細分布。由于震后不同時期的形變具有相似的地下活動過程，依據(jù)上述方法，本文首先聯(lián)合震后升軌83 天、降軌79 天的形變進行反演，獲得一組最優(yōu)斷層幾何參數(shù)，再以觀測周期內(nèi)不同階段的震后形變?yōu)榧s束，最終反演得到各階段對應的余滑分布。

1.2 結(jié)果與分析

以震后第一期的觀測數(shù)據(jù)為參考，將計算所得的其它影像與參考影像的成像時間差視為震后天數(shù)，得到震后的各期形變場，如圖2 所示。可見：形變場時序圖顯示薩波爾扎哈布地震發(fā)生后可提取到顯著的地表形變，其量級和分布均清晰可辨。震后35 天（2017 年12 月17 日），升軌觀測場上已經(jīng)識別到明顯的形變信號，其衛(wèi)星視線（line of sight，縮寫為LOS）向的抬升形變量達到50 mm 左右；升、降軌觀測形變的累積形變量和形變分布范圍均隨著時間推移而不斷增大，升軌形變以LOS 向抬升為主，降軌形變則兼具明顯的LOS 向抬升和沉降信號，形變分布呈“雙橢圓”狀；InSAR 觀測的震后地表形變分布均具有不對稱性，且抬升形變量明顯高于沉降形變量，呈現(xiàn)出類似逆沖斷層典型的地下活動特征。

圖2 2017 年薩波爾扎哈布MW7.3 地震的震后形變場（a） 升軌視線向形變量；（b） 降軌視線向形變量Fig. 2 Post-seismic deformation field of the Sarpol Zahab MW7.3 earthquake in 2017（a） Ascending LOS deformation；（b） Descending LOS deformation

為進一步分析震后形變特征，在主要形變區(qū)選取四個點繪制其時序形變累積圖，如圖3a所示，可見：累積形變量持續(xù)增大；震后初期形變增長較快，隨著時間的推移，形變增長速率呈下降趨勢；震后約220 天時，形變增長明顯減緩；形變序列整體符合對數(shù)函數(shù)變化規(guī)律。

圖3 點形變序列（a）及震后余滑的地震矩張量M0 和矩震級MW 隨時間的變化（b）Fig. 3 Point deformation sequences （a） and variation of the seismic moment tensor M0 and the moment magnitude MW of the afterslip with post-seismic time （b）

根據(jù)震后每一期累積形變進行余滑反演，結(jié)果如圖4 所示，圖中余滑模型的滑移量對應于震后該時期的累積滑移量。反演結(jié)果顯示：滑動以逆沖運動為主，深度上集中分布于沿斷層傾向10—45 km 處，且呈現(xiàn)出主滑區(qū)余滑量隨時間延長逐漸增加的特點；震后180 天左右，在深度50 km 的區(qū)域有一小塊聚簇的滑動，而后此處的累積滑移量逐漸增加，斷層活動向深部延伸也對應其地表形變速率逐漸降低，但累積余滑和形變量仍持續(xù)增長；震后283 天的余滑最顯著，滑移量達到0.7 m，位于斷層平面右下角區(qū)域。

圖4 2017 年薩波爾扎哈布MW7.3 地震的震后余滑分布Fig. 4 The afterslip distribution of Sarpol Zahab MW7.3 earthquake in 2017

利用反演的各期余滑模型計算對應的地震矩和矩震級，結(jié)果如圖3b 所示?？梢姡赫鸷?83 天的余滑釋放的地震矩為1.58×1019N·m，大約相當于一次MW6.74 地震；地震矩和矩震級均呈現(xiàn)出先快速增加再趨于平緩的指數(shù)函數(shù)變化規(guī)律。在時間演化上，該指數(shù)變化規(guī)律與圖3a 的時序累計形變變化規(guī)律形成對應，印證了余滑反演結(jié)果的合理性。

基于反演的最優(yōu)余滑模型，分別模擬計算了升、降軌InSAR 的LOS 向地表形變，并將觀測值與模擬值作差得到殘差分布，如圖5 所示。結(jié)果顯示：模擬形變與觀測形變整體上具有較高的相似度，殘差基本小于1.5 cm；主形變區(qū)的殘差量級較小，而形變場的西南和東北方向部分區(qū)域存在較明顯的殘差，主要來源于觀測誤差。由于研究區(qū)域位于山區(qū)，植被覆蓋茂密，易造成體散射失相關(guān)，且山區(qū)地形起伏較大，因此這些區(qū)域殘存部分大氣延遲誤差，使得時序解算的結(jié)果受到影響。

圖5 2017 年薩波爾扎哈布MW7.3 地震震后余滑的模擬形變場和殘差。（a） 升軌形變場Fig. 5 The simulated post-seismic deformation fields and the residuals based on the afterslip of the Sarpol Zahab MW7.3 earthquake in 2017。 （a） The ascending deformation field

圖5 2017 年薩波爾扎哈布MW7.3 地震震后余滑的模擬形變場和殘差。（b） 降軌形變場Fig. 5 The simulated post-seismic deformation fields and the residuals based on the afterslip of the Sarpol Zahab MW7.3 earthquake in 2017。（b） The descending deformation field

2 薩波爾扎哈布地震對2018 年MW6.0 地震的觸發(fā)影響

2.1 接收斷層的幾何參數(shù)

探究地震間的觸發(fā)關(guān)系需獲得后續(xù)地震事件準確的斷層幾何模型。本文利用地震發(fā)生前后的四景Sentinel-1B 衛(wèi)星SAR 影像，提取后續(xù)賈萬魯?shù)翸W6.0 地震的升、降軌同震地表形變。SAR 影像的詳細參數(shù)列于表2，圖1 顯示所使用的SAR 影像完整地覆蓋了地震區(qū)域。采用“二軌法”差分干涉技術(shù)處理影像數(shù)據(jù)，并引入SRTM4 （Shuttle Radar Topography Mission 4）DEM 數(shù)據(jù)去除地形相位，經(jīng)差分干涉、濾波、相位解纏、地理編碼等步驟，獲得了升、降軌觀測的LOS 向地表形變場，如圖6a 所示。

表2 Sentinel-1B 衛(wèi)星SAR 影像參數(shù)Table 2 SAR image parameters of the Sentinel-1B satellite

由圖6a 可見：InSAR 觀測的最大抬升形變量為5.5 cm，最大沉降形變量為4.5 cm，整體形變條紋比較清晰，說明同震斷層破裂可能未延伸至地表；升軌和降軌觀測值均圍繞震中呈現(xiàn)出“花瓣”狀的形變條紋，升軌觀測值主要位于西側(cè)的兩個“花瓣”上，而降軌主要形變則位于東側(cè)的兩個“花瓣”，觀測結(jié)果表明升軌西側(cè)上、下兩瓣形變分別為抬升和沉降，降軌東側(cè)上、下兩瓣形變分別為沉降和抬升。根據(jù)上述地表形變場的分布特征推測該地震是由近南北走向的走滑斷層運動所致。

同樣采用“兩步法”反演獲得賈萬魯?shù)碌卣鸬臄鄬踊瑒幽Ｐ?，如圖7 所示。結(jié)果顯示：本次地震斷層走向為355.6°，傾角為89.4°，滑動以右旋走滑為主兼具部分正傾滑運動，斷層的最大滑移量達1.2 m，位于地下約6.5 km 處；滑動集中分布于沿斷層面4—10 km 的區(qū)域，最淺層的子斷層滑動不顯著，表明此次地震基本未破裂至地表，符合同震形變觀測結(jié)果。根據(jù)所獲斷層滑動模型正演所得到的模擬地表形變場（圖6b），并與觀測值作差，得到殘差分布（圖6c），結(jié)果顯示升軌東側(cè)與降軌西側(cè)的部分區(qū)域存在部分殘差，這可能是由于主震區(qū)位于山區(qū)，地形起伏使得InSAR 干涉過程中存在部分大氣延遲誤差，對反演精度造成了影響?？傮w而言，本文估計的斷層滑動模型能較好地擬合InSAR 觀測形變場，反演殘差較小，說明斷層參數(shù)的可靠性高。

圖6 2018 年賈萬魯?shù)翸W6.0 地震升軌（上）、降軌（下）的同震形變場及其模擬殘差（a） 觀測形變場；（b） 模擬形變場；（c） 殘差Fig. 6 Coseismic deformation fields and simulation residuals in ascending （upper） and descending （lower） orbits of the Javanrud MW6.0 earthquake in 2018（a） Observed InSAR deformation field；（b） Simulated deformation field；（c） Residual

圖7 InSAR 反演所得的賈萬魯?shù)翸W6.0 地震的斷層滑動模型（a） 斷層滑動空間分布；（b） 斷層平面上的滑動分布，黃色星形代表2018 年賈萬魯?shù)翸W6.0 地震震源在斷層面上的投影位置Fig. 7 Fault slip model of the Javanrud MW6.0 earthquake constrained by InSAR observations（a） 3D view of the fault slip；（b） Slip distribution on the fault plane，where the yellow star represents the projected location of the hypocenter of the 2018 MW6.0 Javanrud earthquake

2.2 應力觸發(fā)關(guān)系

為了研究薩波爾扎哈布地震及其震后效應對賈萬魯?shù)碌卣鸬挠|發(fā)影響，本文結(jié)合Yang 等（2018b）估計的薩波爾扎哈布地震的同震斷層模型和本文反演的震后余滑模型，計算同震及震后作用在賈萬魯?shù)碌卣饠鄬悠矫嫔系膸靵銎屏褢ψ兓?，其中接收斷層參?shù)根據(jù)2.1 節(jié)InSAR 反演所獲的同震斷層幾何參數(shù)設定。計算庫侖應力時，摩擦系數(shù)大多取值于 ［ 0.1，0.8 ］ 區(qū)間內(nèi)（Deng，Sykes，1997；Yanget al，2018a）。本研究測試了該區(qū)間內(nèi)的不同摩擦系數(shù)，分別計算得到相應的庫侖破裂應力變化，結(jié)果表明摩擦系數(shù)的改變對于應力變化的正負值空間分布影響不大，但對應力數(shù)值有一定影響，這與萬永革等（2000）的數(shù)值試驗結(jié)果相同。由于最佳摩擦系數(shù)值難以確定，故本文選取中間值0.4 進行計算。

圖8a 和8b 展示了2017 年薩波爾扎哈布MW7.3 地震在2018 年賈萬魯?shù)翸W6.0 地震斷層平面上的庫侖應力變化場，可見：前者對后者的震源區(qū)域造成了明顯的正值庫侖應力傳輸，正庫侖應力區(qū)域主要位于震源東北方向，且應力增加的幅值相對較大，約為0.12 MPa，震中位置的庫侖應力增加量為0.04 MPa，達到了以往研究中所提出的大約0.01 MPa 的靜態(tài)庫侖應力觸發(fā)閾值（Kinget al，1994；Lin，Stein，2004），表明薩波爾扎哈布地震對于賈萬魯?shù)碌卣鸬陌l(fā)生具有一定的觸發(fā)效應。圖8c 和8d 展示了薩波爾扎哈布地震的震后余滑在賈萬魯?shù)碌卣饠鄬悠矫嫔系膽ψ兓瘓觯梢姡赫鸷笥嗷瑢Z萬魯?shù)碌卣鹫鹪磪^(qū)亦存在正值庫侖應力傳輸，這是由于余滑是同震斷層的繼續(xù)滑動作用，但是相較于同震影響，震后余滑導致的庫侖應力增加量級相對較小，約為0.008 MPa，正庫侖應力區(qū)域主要位于震源東南方向，震中位置的庫侖應力增加值為0.002 MPa。庫侖應力的計算結(jié)果說明賈萬魯?shù)碌卣鹬饕芩_波爾扎哈布主震的觸發(fā)，而震后余滑對其影響較小，總體而言，薩波爾扎哈布地震及其震后運動均對賈萬魯?shù)碌卣鸬陌l(fā)生具有觸發(fā)作用。

圖8 2017 年薩波爾扎哈布MW7.3 地震及其震后余滑在2018 年賈萬魯?shù)翸W6.0 地震發(fā)震斷層上觸發(fā)的庫侖應力變化ΔCFS（a，b） MW7.3 地震對 MW6.0 地震觸發(fā)的庫侖應力變化ΔCFS；（c，d） MW7.3 地震的震后余滑對MW6.0 地震觸發(fā)的庫侖應力變化ΔCFS；（e，f） 本文計算使用的2017 年MW7.3 地震同震斷層模型（Yang et al，2018b）。圖中黃色和紅色星形分別代表2018 年MW6.0 和2017 年MW7.3 地震震源在各自斷層面上的投影位置Fig. 8 The Coulomb failure stress change ΔCFS on the fault plane of 2018 MW6.0 earthquake triggered by the 2017 MW7.3 earthquake and its afterslip（a，b） ΔCFS on the seismogenic fault plane of the MW6.0 earthquake triggered by the MW7.3 earthquake；（c，d） ΔCFS on the fault plane of the MW6.0 earthquake triggered by the MW7.3 earthquake afterslip；（e，f） The coseismic fault of the 2017 MW7.3 earthquake used in the calculation of ΔCFS （Yang et al，2018b）。 The yellow and red stars represent the projected location of the hypocenter of the 2018 MW6.0 and the 2017 MW7.3 events on their fault planes

3 討論與結(jié)論

本文利用SBAS-InSAR 技術(shù)提取的震后形變場（圖2）與已有研究（楊百存等，2018；Yanget al，2019）提取的同震形變結(jié)果具有較高的相似性，可見在本研究監(jiān)測時段內(nèi)，斷層在震后呈現(xiàn)出與同震效應類似的活動特征，即薩波爾扎哈布地震的主震斷層被錯動后持續(xù)滑動所產(chǎn)生的余滑效應為震后地表形變的主要成分，相對于同震形變場，震后地表形變的位置明顯向西南方向發(fā)生偏移，對應于同震斷層錯動方向。

反演估計的余滑模型以逆沖運動為主，這與薩波爾扎哈布地震的同震斷層運動特征相似（Yanget al，2018b；張慶云等，2020）。將震后283 天的余滑分布與已有的同震滑動結(jié)果（Yanget al，2018b）繪制在一起進行對比，如圖9 所示，可見：余滑位于同震滑動的西南方區(qū)域，兩者在空間分布上具有明顯的相關(guān)性，震后283天的余滑累積滑移量僅為0.7 m 左右，而同震滑動最大值可達4.99 m，說明余滑相對于同震滑動并不明顯，且余滑模型（圖6）的傾角較小，主余滑區(qū)位置較淺。上述分析表明，震后余滑作用使同震斷層進一步加劇滑動并擴大了滑動區(qū)域，但整體上并未改變同震斷層滑動分布的特征。

圖9 2017 年MW7.3 地震的同震及震后斷層滑動分布Fig. 9 The co- and post-seismic fault slip distribution of the 2017 MW7.3 earthquake

本文基于InSAR 形變反演了賈萬魯?shù)碌卣鸬耐饠鄬幽Ｐ停⑵渥鳛閼τ嬎愕慕邮諗鄬?。計算的庫侖應力變化表明，薩波爾扎哈布地震及其震后余滑對賈萬魯?shù)碌卣鸬陌l(fā)生具有不同程度的觸發(fā)作用，其中同震在斷層面上造成的庫侖應力增值達0.12 MPa，震后余滑則繼續(xù)對斷層面應力狀態(tài)產(chǎn)生影響，根據(jù)圖8b 和8d 中正、負應力變化分布可知，余滑會增大或減小斷層面上部分區(qū)域的應力狀態(tài)。主震位錯產(chǎn)生的庫侖應力變化是強余震活動的重要原因，但并非全部原因（郝平等，2004），故本文從地質(zhì)構(gòu)造背景方面進一步探究導致地震發(fā)生的可能因素。

賈萬魯?shù)碌卣鹞挥谠窳_斯山脈褶皺逆沖帶的西北部，該區(qū)域的斷裂分布（圖1）顯示研究區(qū)主要受到北東—南西向的擠壓應力作用（文磊等，2015），反演的斷層面滑動分布（圖7）顯示其同震運動與區(qū)域背景構(gòu)造運動趨勢相關(guān)。因此，需要完整地解讀震間、同震和震后這一地震周期中的應力傳輸過程，并結(jié)合研究區(qū)域的地質(zhì)活動背景，才能更好地評估地震發(fā)生的危險性。本文僅計算了靜態(tài)應力觸發(fā)作用，而冀戰(zhàn)波等（2014）針對2008 年于田MS7.3 地震對后續(xù)地震的庫侖應力觸發(fā)研究認為，對于近場區(qū)域的余震，動態(tài)庫侖應力能夠更好地解釋其分布。盡管本文結(jié)果已顯示薩波爾扎哈布地震對賈萬魯?shù)碌卣鸬撵o態(tài)應力觸發(fā)產(chǎn)生了影響，但同時考慮動態(tài)庫侖應力的影響，或許更能解釋后續(xù)地震的成因，本文將會在下一步研究中將其效應考慮在內(nèi)。

強震后伴隨的大量余震可能會對震后形變場產(chǎn)生影響。萬永革等（2005）計算了地震斷層面及附近余震對震后位移場的影響，結(jié)果顯示大量余震發(fā)生后可產(chǎn)生能被GPS 觀測到的地表位移。由于目前尚未獲得大量包含詳細震源機制的余震信息，在InSAR 震后形變場中也未觀測到明顯的余震相關(guān)的形變信號，因此本文未考慮余震的影響。根據(jù)已有研究的余震地表位移場的分析，多數(shù)余震的位移場方向與主震基本一致，且余震釋放能量遠小于主震（萬永革等，2005，2009），本文觀測到的震后形變模式基本與主震相同，可以推測震后形變場中的少量余震位移并不會改變本研究中主震和余滑對賈萬魯?shù)碌卣鸬挠|發(fā)影響。

除余震位移場之外，黏彈性松弛效應也已被證明在短期的震后形變中存在一定影響（Guoet al，2019），但本研究也顯示，余滑在斷層近場起主導作用，通過中遠場的觀測數(shù)據(jù)才能更好地判讀黏彈性松弛模型的擬合結(jié)果。本文的研究范圍主要在近場，且目前的觀測結(jié)果中未發(fā)現(xiàn)明顯的遠場形變信號，故反演中未考慮震后的黏彈性效應。此外，Yang 等（2019）將薩波爾扎哈布地震后220 天的二維形變與其三維同震形變分量進行對比的結(jié)果顯示，二者在垂直和東西方向上的變形趨勢相同，表明在這一時期震后形變的主要成分是余滑效應，本文結(jié)果與其相符，同時也說明黏彈性松弛形變的量級較小。相應地，我們也認為短期內(nèi)該地區(qū)的黏彈性效應所產(chǎn)生的庫侖應力變化影響不大，如果針對長時期的地震觸發(fā)研究，黏彈性效應將是一個必須考慮的因素。上述黏彈性效應造成的地表形變雖然并不顯著，但在實際余滑反演中，可能會對余滑量級造成一定程度的高估，但顯而易見的是，本研究中余滑所造成的地表形變?nèi)匀徽贾鲗У匚?，因此，本文反演的余滑模型仍具有較高的可靠性。