張文婷 季靈運(yùn)* 朱良玉 蔣鋒云 徐曉雪

1)中國地震局第二監(jiān)測(cè)中心，西安 710054 2)國家遙感中心地質(zhì)災(zāi)害研究部，北京 100036

圖1 伽師地震區(qū)域構(gòu)造Fig.1 Tectonic setting map of the Jiashi earthquake.a 黑色細(xì)線表示斷層(F1柯坪斷裂；F2奧茲格爾他烏斷裂；F3托克散阿塔能拜勒褶皺-逆斷裂帶；F4科克布克三山褶皺-逆斷裂帶；F5奧依布拉克褶皺-逆斷裂帶；F6喀拉鐵克斷裂；F7邁丹斷裂；F8皮羌斷裂；F9巴什托普斷裂)，灰色粗線表示二級(jí)塊體劃分，白色方框?yàn)樯?、降軌SAR數(shù)據(jù)覆蓋范圍，震源機(jī)制解數(shù)據(jù)來自GCMT(1976—2020年)(1)https: ∥www.globalcmt.org/。和USGS(2020年)(2)https: ∥earthquake.usgs.gov/earthquakes/。，藍(lán)色箭頭表示GPS水平速度場(chǎng)；b 黑色實(shí)線表示二級(jí)塊體劃分，紅色圓圈表示1900年以來發(fā)生的7級(jí)以上地震，紅色方框表示圖a的位置

表 1 不同機(jī)構(gòu)提供的震源機(jī)制解Table 1 Focal mechanism solutions of Jiashi MS6.4 earthquake

0 引言

據(jù)中國地震臺(tái)網(wǎng)測(cè)定，2020年1月19日21時(shí)27分新疆喀什地區(qū)伽師縣(39.83°N，77.21°E)發(fā)生了6.4級(jí)地震，震中位于南天山西段的柯坪推覆體附近(圖1)，距離伽師縣56km，震源深度約為16km(孟令媛等，2020；聶曉紅等，2020)，此次地震是自1997—2003年伽師6～7級(jí)震群事件以來發(fā)生的又一次強(qiáng)烈地震。多家機(jī)構(gòu)對(duì)此次地震的震源機(jī)制解進(jìn)行了解算，得到的結(jié)果均顯示此次地震為一次逆沖型地震(表1)。目前已有學(xué)者基于InSAR技術(shù)獲取了此次地震的同震形變場(chǎng)(溫少妍等，2020)，并反演了斷層面滑動(dòng)分布(李成龍等，2020)。本文在此基礎(chǔ)上優(yōu)化了斷層面幾何構(gòu)形及滑動(dòng)分布反演方法，獲取了此次地震更精細(xì)的斷層面滑動(dòng)分布，然后對(duì)區(qū)域GPS面應(yīng)變率、構(gòu)造動(dòng)力背景和區(qū)域深淺部構(gòu)造進(jìn)行了綜合分析，探討了2020年新疆伽師6.4級(jí)地震與1997—2003年伽師6～7級(jí)震群的相互關(guān)系。

1 區(qū)域構(gòu)造背景

2020年新疆伽師6.4級(jí)地震的震中位于南天山前陸褶皺沖斷帶，其北部為天山褶皺帶，西南為西昆侖-帕米爾高原，東部和南部為剛性的塔里木盆地(喬學(xué)軍等，2007)(圖1)。受印度板塊與歐亞板塊持續(xù)碰撞的影響，該區(qū)域的構(gòu)造活動(dòng)主要表現(xiàn)為天山地塊逆沖于塔里木地塊之上(田勤儉等，2006)，GPS觀測(cè)到該區(qū)域的地殼縮短率為6～8mm/a(王琪等，2001)。在這樣的逆沖環(huán)境下，伽師地區(qū)中強(qiáng)地震頻發(fā)，僅過去的50a內(nèi)，距此次震中50km的范圍內(nèi)就發(fā)生了13次具有多種震源機(jī)制類型的6級(jí)以上地震(周仕勇等，2006；張竹琪等，2008)。

從區(qū)域相對(duì)較小空間尺度的構(gòu)造背景分析，2020年新疆伽師6.4級(jí)地震的震源區(qū)位于南天山山前的柯坪推覆體。新生代期間，特別是第四紀(jì)以來，天山地區(qū)強(qiáng)烈而持久的再生造山作用不僅形成了分布式逆斷裂作用的格局，使地殼近SN向縮短和垂向隆升，還在天山南北麓形成逆斷裂-褶皺帶或背馱式逆斷裂推覆體系，引起天山山脈不斷地向外側(cè)擴(kuò)展，逆斷裂-褶皺帶的主體活動(dòng)期也從山體一側(cè)的中新世早期向前陸盆地逐漸變新(鄧起東等，2000)?？缕和聘搀w的形成及擴(kuò)展也與滑脫面的伸展有關(guān)，天山主體沿高角度逆斷裂向南、北兩側(cè)前陸盆地推擠的過程中，在盆地下順寒武系底界面形成了低角度滑脫斷裂，該滑脫斷裂的形成和不斷向盆地方向的伸展，導(dǎo)致之上多排褶皺-逆斷裂帶的逐步形成與推覆體的擴(kuò)展形成了柯坪推覆體(宋方敏等，2006)。推覆體由一系列近NEE-SWW走向、向S逆沖的薄皮式?jīng)_斷帶組成(Allenetal.，1994)，EW長(zhǎng)約300km，SN寬60～140km，北以邁丹斷裂與南天山晚古生代造山帶為界，南以柯坪斷裂與塔里木盆地為界，在新生代天山造山帶變形向盆地?cái)U(kuò)展的過程中扮演著極其重要的角色(張耀，2019)?？缈缕簲鄩K的人工地震剖線揭示位于基底及其蓋層間存在一個(gè)總體向天山傾斜的滑脫層，具有上陡下緩的 “犁式”特點(diǎn)(楊曉平等，2008)，塔里木盆地的基底沿滑脫面自北向南逆沖推覆，推覆體南淺北深，滑脫面上部發(fā)育高角度斷坡斷層，這些斷層或出露地表或隱伏其下，與滑脫層形成多排近EW向的疊瓦狀鏟式斷層(李杰等，2015)。受地殼橫向縮短不均勻性的影響，柯坪推覆體被NNW向的皮羌 “類轉(zhuǎn)換斷層”(F8)分割為東、西2段，此次地震就發(fā)生在柯坪推覆體的西段，該段由多排走向NEE—近EW向的倒轉(zhuǎn)背斜褶皺、褶皺前緣斷裂及其間背馱式盆地或谷地等組成，且具有明顯的構(gòu)造分段性特征，從南到北依次由柯坪塔格逆斷裂-背斜(F1)、奧茲格爾他烏逆斷裂-背斜(F2)、托克散阿塔能拜勒逆斷裂-背斜(F3)、科克布克三山逆斷裂-背斜(F4)和奧依布拉克逆斷裂-背斜(F5)等組成(徐錫偉等，2006)，每個(gè)逆斷裂-背斜以單面山的形式出現(xiàn)，自新生代晚期向塔里木盆地逆沖推覆，地殼縮短量約為45km，地殼縮短率為37%(楊曉平等，2006)，其中早期形成的褶皺-逆斷裂帶在北側(cè)，晚期形成的褶皺-逆斷裂帶在南側(cè)，反映推覆體由北向南的擴(kuò)展(宋方敏等，2006)。地貌測(cè)量結(jié)果顯示，柯坪推覆體西段的逆斷裂-褶皺帶活動(dòng)性從北向南逐漸增強(qiáng)，表現(xiàn)出南天山多排逆斷層在前展模式中的活動(dòng)強(qiáng)度由后排構(gòu)造向前排構(gòu)造轉(zhuǎn)移的特點(diǎn)。從區(qū)域地層厚度分析，認(rèn)為柯坪逆斷裂-褶皺帶西段滑脫層的埋深約為9km，東段滑脫層的埋深約為5km(楊曉平等，2008)，向N緩傾的滑脫面在天山根部會(huì)聚到邁丹斷裂上(吳傳勇等，2014)。2020年伽師6.4級(jí)地震的發(fā)生為深入理解柯坪推覆體的動(dòng)力學(xué)過程帶來了新的啟示。

2 2020年新疆伽師6.4級(jí)地震的InSAR同震形變場(chǎng)

此次新疆伽師6.4級(jí)地震的震中區(qū)域人口稀少，傳統(tǒng)大地測(cè)量觀測(cè)站點(diǎn)分布稀疏，難以獲得高分辨率的同震形變場(chǎng)(圖1)。InSAR技術(shù)作為常用的影像大地測(cè)量方法，已被廣泛應(yīng)用于同震形變場(chǎng)的獲取及發(fā)震斷層幾何參數(shù)及其滑動(dòng)分布的反演(Massonnetetal.，1998；喬學(xué)軍等，2007；單新建等，2009；張國宏等，2010；宋小剛等，2015；季靈運(yùn)等，2017；屈春燕等，2017)?；诖?，本文利用InSAR技術(shù)獲取了2020年伽師6.4級(jí)地震的同震形變場(chǎng)，并以此為約束反演其發(fā)震斷層的幾何參數(shù)及斷層面的滑動(dòng)分布特征。

本文利用升、降軌Sentinel-1 SAR衛(wèi)星數(shù)據(jù)獲取了伽師6.4級(jí)地震的同震形變場(chǎng)(表2)。InSAR數(shù)據(jù)處理在GAMMA商業(yè)軟件平臺(tái)上進(jìn)行，首先基于美國NASA發(fā)布的30m空間分辨率的SRTM DEM地形數(shù)據(jù)模擬并消除地形相位，得到研究區(qū)的InSAR同震干涉圖(圖2)，然后基于Delaunay三角網(wǎng)的最小費(fèi)用網(wǎng)絡(luò)流算法對(duì)差分干涉圖進(jìn)行相位解纏，解纏起始點(diǎn)選取在遠(yuǎn)離震中的穩(wěn)定點(diǎn)。使用遠(yuǎn)場(chǎng)數(shù)據(jù)進(jìn)行二次多項(xiàng)式擬合去除形變場(chǎng)中的殘余軌道誤差(Fialkoetal.，2006)，基于地形相關(guān)法估計(jì)并去除大氣相位延遲誤差，最后根據(jù)SAR衛(wèi)星參數(shù)將相位轉(zhuǎn)化為位移，得到此次地震的升、降軌LOS向同震形變場(chǎng)(圖3a，d)。

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

圖2 中的干涉條紋變化范圍顯示此次地震形變場(chǎng)EW向跨度約40km，SN向跨度約20km，圖中1個(gè)色周對(duì)應(yīng)2.8cm的位移量，條紋從正變化為負(fù)代表地表沿LOS向(雷達(dá)視線方向)靠近衛(wèi)星運(yùn)動(dòng)，干涉圖的整體相干性較高且震源區(qū)不存在由于形變梯度過大導(dǎo)致的失相干，表明同震破裂未達(dá)到地表。形變區(qū)域內(nèi)有柯坪和奧茲格爾他烏2條逆沖活動(dòng)斷裂，干涉條紋的變化顯示升、降軌形變場(chǎng)在柯坪斷裂以北(斷層上盤)為顯著上升，在奧茲格爾他烏斷裂以有北少量下沉，2條軌道的形變趨勢(shì)基本一致，表明此次地震以垂向運(yùn)動(dòng)為主導(dǎo)，符合逆沖型破裂引起的形變特征。最終獲取的同震形變(圖3a，d)結(jié)果顯示，升軌結(jié)果LOS向最大抬升量為5.9cm，最大沉降量為3.7cm；降軌結(jié)果LOS向最大抬升量約為6.4cm，最大沉降量約為2cm，與前人的結(jié)果較為一致(李成龍等，2020；溫少妍等，2020)。結(jié)合InSAR形變場(chǎng)分布及已有斷層的空間分布，認(rèn)為此次地震為典型的逆沖型地震破裂事件，其發(fā)震斷層為柯坪斷裂。

圖2 InSAR同震差分干涉圖升(a)、降軌(b)Fig.2 Coseismic ascending(a)and descending(b) interferograms from InSAR.

3 2020年新疆伽師6.4級(jí)地震的同震滑動(dòng)分布

3.1 均勻滑動(dòng)模型

由于缺少此次地震發(fā)震斷層的先驗(yàn)信息，在反演斷層面精細(xì)化滑動(dòng)分布前，需采用均勻滑動(dòng)模型對(duì)發(fā)震斷層的位置、震源深度、走向、傾角等幾何參數(shù)進(jìn)行反演。本文采用一種基于貝葉斯方法的均勻滑動(dòng)反演方法，該反演方法假設(shè)地震由單一斷層面滑動(dòng)引起，發(fā)震斷層面是一個(gè)嵌入在介質(zhì)均勻且各向同性的彈性半空間中的矩形面(Okada，1985)，其滑移量和滑動(dòng)角都均勻分布。反演時(shí)給定較大范圍的上、下界并隨機(jī)給出相應(yīng)的初始值，采用馬爾科夫鏈蒙特卡洛方法對(duì)參數(shù)空間進(jìn)行采樣，基于貝葉斯方法構(gòu)建后驗(yàn)概率函數(shù)，通過迭代不斷對(duì)參數(shù)空間進(jìn)行優(yōu)化直至收斂，進(jìn)而計(jì)算出最佳擬合參數(shù)(Bagnardietal.，2018)。圖3 所示的藍(lán)色矩形為參數(shù)出現(xiàn)的頻數(shù)，紅色豎線為最優(yōu)解。

為提高計(jì)算效率，在反演前利用四叉樹分割方法對(duì)原始的高分辨率InSAR觀測(cè)值進(jìn)行降采樣，將采樣的形變梯度閾值設(shè)為2.8cm，最終得到1i327個(gè)升軌觀測(cè)值和1i680個(gè)降軌觀測(cè)值，既保留了形變細(xì)節(jié)特征又能減少數(shù)據(jù)冗余。反演得到的發(fā)震斷層參數(shù)如表3 所示，說明此次地震是以傾滑運(yùn)動(dòng)為主的逆沖型地震。圖3c、3f為觀測(cè)值與模擬值間的差值，可以看出均勻滑動(dòng)模型反演結(jié)果的殘差較小，表明了反演結(jié)果的可靠性(圖4)。

圖3 均勻滑動(dòng)模型參數(shù)直方圖Fig.3 Frequency histograms of uniform slip model parameters.

表 3 均勻滑動(dòng)反演的斷層參數(shù)Table 3 Fault parameters derived by uniform slip inversion

圖4 均勻滑動(dòng)模型的反演結(jié)果Fig.4 Simulation results of the uniform-slip model.a、d 升軌和降軌的InSAR同震形變場(chǎng)；b、e 相應(yīng)的模擬同震形變場(chǎng)；c、f 殘差圖

3.2 分布式滑動(dòng)模型

為了獲取更為精細(xì)的斷層面滑動(dòng)分布，本文以均勻滑動(dòng)模型反演得到的發(fā)震斷層幾何參數(shù)為基礎(chǔ)，引入基于約束最小二乘的分布式滑動(dòng)反演方法，計(jì)算伽師地震發(fā)震斷層面的精細(xì)滑動(dòng)分布。

考慮到震源區(qū)的構(gòu)造背景，本文在前人研究的基礎(chǔ)上(李成龍等，2020)，對(duì)分布式滑動(dòng)反演的斷層面幾何構(gòu)型進(jìn)行了改進(jìn)。首先將均勻滑動(dòng)模型反演得到的斷層長(zhǎng)度擴(kuò)展至40km?？缈缕簲鄩K的人工地震剖面揭示區(qū)域基底及其蓋層間存在一個(gè)總體向天山傾斜的滑脫層，具有上陡下緩的 “犁式”特點(diǎn)(楊曉平等，2008)，滑脫面上部發(fā)育高角度斷坡斷層，這些斷層或出露地表或隱伏其下，與滑脫層形成多排近EW向的疊瓦狀鏟式斷層(李杰等，2015)。因此，根據(jù)均勻滑動(dòng)反演的發(fā)震斷層參數(shù)及柯坪斷裂在地表出露位置間的幾何關(guān)系，將斷層面設(shè)置為鏟狀，其傾角從地表到地表之下4.25km深處設(shè)為57.8°，4.25～8km深處設(shè)為10.7°，這樣的斷層幾何參數(shù)設(shè)置與深部地球物理結(jié)果一致。

最后，將斷層面剖分為930個(gè)邊長(zhǎng)為2km的等邊三角形單元，基于約束最小二乘方法計(jì)算每個(gè)三角形單元的滑動(dòng)量，通過最小化模型和觀測(cè)值之間的殘差確定最佳參數(shù)。通過權(quán)衡數(shù)據(jù)擬合程度與斷層滑動(dòng)粗糙度，對(duì)相鄰斷層片的滑動(dòng)量施加平滑約束。

圖5 分布式滑動(dòng)模型的擬合結(jié)果Fig.5 Simulation results of the distributed-slip model.a、d 升軌和降軌的InSAR同震形變場(chǎng)；b、e 相應(yīng)的模擬同震形變場(chǎng)；c、f 殘差圖

與均勻滑動(dòng)模型結(jié)果相比，分布式滑動(dòng)模型結(jié)果的殘差明顯減小(圖5)。圖6 為獲取的斷層滑動(dòng)分布圖，結(jié)果顯示斷層面的最大傾滑量約為0.29m，位于地下約5km深處，最大走滑量約為0.03m。計(jì)算得到同震的累計(jì)地震矩為1.73×1018N·m，合矩震級(jí)MW6.1。

圖6 伽師地震的同震滑動(dòng)分布Fig.6 Slip distributions for the Jiashi earthquake.a 走滑分量；b 傾滑分量

4 討論

4.1 基于現(xiàn)今GPS的區(qū)域應(yīng)變率場(chǎng)分析

南天山前陸褶皺沖斷帶是大陸內(nèi)部典型的逆沖斷裂帶，地質(zhì)結(jié)果和GPS觀測(cè)數(shù)據(jù)表明天山南、北兩側(cè)以逆斷層-褶皺為代表的薄皮推覆構(gòu)造自晚第四紀(jì)以來活動(dòng)非常強(qiáng)烈，這些構(gòu)造帶以水平向的地殼縮短和垂向的構(gòu)造抬升吸收調(diào)節(jié)了天山地區(qū)至少一半的縮短變形(鄭文俊等，2019)。此次伽師地震所處的柯坪推覆體作為塔里木盆地周緣惟一一個(gè)以寒武系為滑脫層的薄皮沖斷帶，在新生代天山造山帶變形向盆地內(nèi)部擴(kuò)展的過程中扮演著極其重要的角色(張耀，2019)。因此，對(duì)柯坪推覆體及周邊區(qū)域的形變特征進(jìn)行研究有助于分析此次地震的發(fā)震背景。

本文基于GPS水平速度場(chǎng)(Wangetal.，2020)和Shen等(2015)提出的應(yīng)變計(jì)算方法，處理并獲取了柯坪地塊及其周邊地區(qū)的面應(yīng)變率場(chǎng)(圖7)。研究區(qū)域的面應(yīng)變率整體為負(fù)值，說明該區(qū)域以擠壓逆沖變形為主，與區(qū)域所處的逆沖構(gòu)造背景一致，區(qū)域內(nèi)發(fā)生的中強(qiáng)地震也以逆沖型居多。值得注意的是，在塔里木盆地與天山交界帶西側(cè)存在一處面應(yīng)變率顯著降低的區(qū)域(圖7 中的藍(lán)色區(qū)域)，說明此處受到擠壓更為強(qiáng)烈，此次伽師地震的震中位于該區(qū)域中應(yīng)變率梯度變化較大處。此外，在該低值區(qū)東側(cè)，為了調(diào)節(jié)由于柯坪推覆體前展變形速率差異較大引起的橫向擠壓縮短不均勻性，派生出NNW向的皮羌 “類轉(zhuǎn)換斷層”，該斷層由柯坪塔格推覆構(gòu)造系西段各推覆體前緣逆斷裂東端的NNW向弧形拐彎段串聯(lián)而成，是一條不連續(xù)的、代表兩側(cè)不等量地殼縮短的邊界(徐錫偉等，2006)，該斷裂歷史上也發(fā)生過6.8級(jí)地震，表明其現(xiàn)今可能處于活動(dòng)狀態(tài)。

圖7 伽師地震及其鄰近地區(qū)的面應(yīng)變率分布圖Fig.7 The plane strain rate field around the Jiashi earthquake and its surrounding areas.

圖8 柯坪推覆體的深、淺部構(gòu)造剖面圖Fig.8 The shallow and deep tectonic profile of the Kalpin nappe.a 地形與垂直于柯坪斷裂走向的GPS速率剖面；b 小地震剖面與深部構(gòu)造解譯。AA′剖面的位置見圖1，黑色曲線代表F1、F2、F6、F7 4條斷層和深部滑脫面，小地震數(shù)據(jù)時(shí)間范圍為1980—2015年

4.2 2020年新疆伽師6.4級(jí)地震的發(fā)震斷層與構(gòu)造動(dòng)力背景分析

柯坪推覆體由多條近EW向的逆斷裂-背斜帶組成，沿著統(tǒng)一的滑脫層由北向南逆沖推覆(張耀，2019)，其滑脫面發(fā)育在地下5～9km深處，向N緩傾的滑脫面在天山根部會(huì)聚到邁丹斷裂上(吳傳勇等，2014)，推覆體南淺北深，總體向N傾(曲國勝等，2003)。根據(jù)前文分析，此次地震的發(fā)震斷層為鏟型結(jié)構(gòu)，從地表至地下深部由高傾角逐漸過渡為低傾角，震源深度約為5km，印證了柯坪地塊為薄皮推覆運(yùn)動(dòng)的地學(xué)認(rèn)識(shí)(陳杰等，2001)。

圖8 指示了此次地震的發(fā)震構(gòu)造動(dòng)力背景，沿過震中的剖線AA′(位置見圖1)繪制出柯坪推覆體的構(gòu)造剖面圖(田勤儉等，2006；李杰等，2015)，將剖面線兩側(cè)50km范圍內(nèi)的GPS水平向速率(Wangetal.，2020)投影到垂直于柯坪斷裂的方向，同時(shí)提取了剖面線兩側(cè)10km范圍內(nèi)的小地震精定位數(shù)據(jù)及地形數(shù)據(jù)(圖8)。可以看出，垂直于柯坪斷裂的GPS運(yùn)動(dòng)速率自北向南呈遞減趨勢(shì)，速率的最大差值約為5mm/a，說明柯坪推覆體受SN向擠壓變形較大。因此，此次地震可能是由柯坪地塊向塔里木盆地推覆的過程中前方?jīng)_斷席發(fā)生錯(cuò)斷所引起的，且發(fā)震斷層位于柯坪推覆體前緣小地震較稀疏的區(qū)域。值得注意的是，西南天山的構(gòu)造變形不僅分配在推覆體前緣的新生代逆斷裂-褶皺帶上，通過山根處的邁丹斷裂(F7)全新世活動(dòng)證據(jù)發(fā)現(xiàn)推覆體根部斷裂現(xiàn)在仍強(qiáng)烈活動(dòng)，吸收了一部分縮短變形(吳傳勇等，2014)，從圖7 可以看出其小地震活動(dòng)也較為稀疏，該區(qū)的地震危險(xiǎn)性不容忽視。

4.3 2020年新疆伽師6.4級(jí)地震與1997—2003年伽師6～7級(jí)震群相互關(guān)系的討論

1997—2003年，南天山前陸環(huán)境不同斷裂帶相繼破裂引發(fā)了3組6～7級(jí)強(qiáng)烈地震(圖1)，震群主要位于塔里木盆地邊緣的麥蓋提斜坡帶與喀什凹陷的過渡邊界。由于該區(qū)域位于塔里木盆地西北端、帕米爾弧型構(gòu)造東北部以及天山褶皺帶柯坪塊體的交會(huì)部位，地質(zhì)構(gòu)造環(huán)境和受力因素復(fù)雜(高國英等，2004)，因此地震的震源機(jī)制復(fù)雜多變，既有逆沖型，也有走滑和正斷型(趙翠萍等，2008)(圖1)。其中，1997年1月—2003年1月伽師震源區(qū)地震的震源機(jī)制解以走滑型為主，且有一定數(shù)量的地震為正斷型，與區(qū)域所處的逆沖環(huán)境不一致。而2003年巴楚-伽師地震的震源機(jī)制解則以逆斷型為主，與區(qū)域震源機(jī)制解類型較為一致(徐錫偉等，2006)。

地震是新構(gòu)造運(yùn)動(dòng)的表現(xiàn)形式之一，地震發(fā)生的地點(diǎn)、強(qiáng)度與地殼最新的構(gòu)造變動(dòng)樣式有著密切的聯(lián)系。已有的地震精確定位、震源機(jī)制解、地殼—上地幔三維S波速度結(jié)構(gòu)和各種探測(cè)成果顯示，1997年伽師強(qiáng)震群的主體發(fā)震構(gòu)造為NNW向下蘇洪隱伏斷層，該斷層為右旋走滑的NNW向 “類轉(zhuǎn)換斷層”，起到分割、調(diào)節(jié)西大東小的地殼縮短邊界轉(zhuǎn)換作用(劉啟元等，2000；徐錫偉等，2006)。而阿圖什1996年6.9級(jí)和1998年6.0級(jí)地震、2003年巴楚-伽師地震、柯坪塊體發(fā)生的多次6級(jí)地震和此次伽師6.4級(jí)地震，無論是震源斷錯(cuò)性質(zhì)還是主壓應(yīng)力方向或震源破裂特征都很一致，反映出擠壓環(huán)境條件下的地殼縮短作用，表現(xiàn)出柯坪塊體長(zhǎng)期以來的斷層運(yùn)動(dòng)特性。

圖9 柯坪推覆體及鄰近區(qū)域的構(gòu)造背景示意圖Fig.9 Schematic map showing regional tectonic background of the Kalpin nappe and surrounding regions.

從1997—2003年伽師6～7級(jí)地震震群與此次伽師6.4級(jí)地震所處的構(gòu)造環(huán)境來看，均發(fā)生在南天山逆沖推覆構(gòu)造體系的前陸盆地環(huán)境中。對(duì)柯坪推覆體的構(gòu)造動(dòng)力環(huán)境進(jìn)行簡(jiǎn)化，如圖8 所示。可以看出，在整體逆沖推覆的構(gòu)造背景中，存在主逆斷層(柯坪斷裂)以及在前陸盆地下方發(fā)育的隱伏逆沖斷層(F9)，隱伏斷裂帶的形成可能是由塔里木盆地北緣的非均勻變形及現(xiàn)今構(gòu)造運(yùn)動(dòng)所引起(劉啟元等，2000)，隨著逆沖推覆帶向盆地發(fā)展，一些斷層向前伸展到前陸盆地下方，被較厚的新生代沉積層覆蓋，形成隱伏斷層，這些隱伏斷層也具有孕育中強(qiáng)地震的能力(張竹琪等，2008)。在逆沖推覆體上，由于巖石力學(xué)屬性存在差異且區(qū)域構(gòu)造應(yīng)力場(chǎng)在不同部位也存在一定的差異性，故在大的逆沖推覆背景下可能會(huì)伴生發(fā)育規(guī)模相對(duì)較小的正斷層，甚至還伴生發(fā)育與主逆斷層走向近垂直的橫向走滑斷層(皮羌斷裂、下蘇洪隱伏斷層)，這些伴生正斷層和走滑斷層也都具備發(fā)生中強(qiáng)地震的能力。因此，根據(jù)區(qū)域構(gòu)造動(dòng)力背景及區(qū)域構(gòu)造分布進(jìn)行綜合分析，認(rèn)為1997—2003年伽師6～7級(jí)地震震群與此次伽師6.4級(jí)地震的發(fā)生機(jī)理應(yīng)為南天山前陸環(huán)境中不同級(jí)別、不同性質(zhì)的斷裂相繼破裂形成不同類型的破裂事件，應(yīng)為同源事件，共同受控于南天山的逆沖推覆作用。

5 結(jié)論

(1)本文利用升、降軌Sentinel-1 SAR衛(wèi)星數(shù)據(jù)獲取了2020年新疆伽師6.4級(jí)地震的InSAR同震形變場(chǎng)，并以InSAR形變場(chǎng)為約束，反演獲取了斷層面的精細(xì)滑動(dòng)分布。結(jié)果顯示，此次地震的發(fā)震斷層為柯坪斷裂。斷層面最大傾滑量約為0.29m，位于地下約5km深處，通過同震釋放能量計(jì)算得到的震級(jí)為MW6.1，是一次以逆沖為主兼具少量走滑運(yùn)動(dòng)的地震破裂事件。

(2)2020年新疆伽師6.4級(jí)地震是發(fā)生在南天山前陸環(huán)境中的一次典型逆沖破裂事件，與1997—2003年在該區(qū)域發(fā)生的多次6～7級(jí)地震一起構(gòu)成系列強(qiáng)震群。從區(qū)域構(gòu)造動(dòng)力背景和區(qū)域深淺部構(gòu)造綜合分析，認(rèn)為1997—2003年伽師6～7級(jí)地震震群與此次伽師6.4級(jí)地震的發(fā)生機(jī)理應(yīng)為南天山前陸環(huán)境中不同級(jí)別、不同性質(zhì)的斷裂相繼破裂形成不同類型的破裂事件，共同受控于南天山的逆沖推覆作用。

致謝中國地震局地球物理研究所房立華研究員為本研究提供了小地震數(shù)據(jù)；歐洲航天局提供了Sentinel-1衛(wèi)星影像。在此一并表示感謝!