劉艷慧 朱良玉 季靈運(yùn) 徐 晶 張文婷 徐曉雪

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

鮮水河斷裂位于青藏高原東緣，是一條分割川滇塊體和巴顏喀拉塊體的弧形左旋走滑大型活動斷裂帶[1](圖1)，在青藏高原東部晚新生代地殼變形中起著重要的調(diào)節(jié)作用。由于獨(dú)特的地理位置與構(gòu)造背景，鮮水河斷裂帶成為中國大陸西部最主要的強(qiáng)震活動帶之一，1725年至今共發(fā)生6級以上地震22次，其中7級以上地震8次[2]，僅中南段(道孚-石棉)于近300 a中就發(fā)生過10次6.5級以上地震[3]，顯示出在近代具有強(qiáng)烈的活動性。2014-11-22康定發(fā)生的6.3級地震是距今最近的一次強(qiáng)震，打破了鮮水河斷裂帶30 a未曾有大震發(fā)生的記錄[4]，也為認(rèn)識地震的震后形變機(jī)制、探索區(qū)域流變結(jié)構(gòu)提供機(jī)會，進(jìn)而為判定后續(xù)地震危險(xiǎn)性提供重要依據(jù)。

圖中黑色實(shí)線表示斷層，紅色圓圈表示5級以上地震，紅色五角星表示康定地震震中，藍(lán)色箭頭表示GPS水平速度場(ITRF2008參考框架下)，藍(lán)色框表示InSAR數(shù)據(jù)范圍；左下角小圖中黑色實(shí)線表示二級塊體劃分，紅色圓點(diǎn)表示1990年后7級以上地震，紅色方框表示大圖位置圖1 鮮水河斷裂帶構(gòu)造背景Fig.1 Tectonic setting of the Xianshuihe fault zone

巖石圈介質(zhì)的力學(xué)結(jié)構(gòu)(尤其是下地殼和上地幔的流變學(xué)結(jié)構(gòu))對地殼形變特征起著決定性作用?；诘卣鹫鸷笳硰椥运沙诂F(xiàn)象，學(xué)者們找到一種認(rèn)識巖石圈流變結(jié)構(gòu)的可行方法：一系列周期性的地震造成區(qū)域應(yīng)力場改變，觸發(fā)地殼中層粘彈性應(yīng)力松弛響應(yīng)，從而產(chǎn)生與時(shí)間相關(guān)的地表形變；運(yùn)用地震周期模型模擬這些與時(shí)間相關(guān)的震后地表形變，就可以反演該地區(qū)的流變結(jié)構(gòu)。例如，在青藏高原內(nèi)部，Ryder等[5-6]利用1997年瑪尼7.6級地震震后形變，運(yùn)用地震周期模型，通過瞬時(shí)形變速率與深部應(yīng)力的粘彈性松弛的聯(lián)系推斷出該地區(qū)存在一個(gè)低粘度(粘度系數(shù)為4×1018Pa·s)的中地殼層。

前人[7-8]通過地質(zhì)研究和地球物理測深及成像方法均揭示出青藏高原東緣地區(qū)存在與中下地殼物質(zhì)通道流密切相關(guān)的低速低阻異常，并推測在青藏高原東部地殼內(nèi)的高導(dǎo)層(HCL)具有較低的粘度。然而，前人的研究結(jié)果只是定性描述青藏高原東緣區(qū)域下地殼存在低阻異常，缺乏對其物理參數(shù)的具體研究。康定地震發(fā)生后，不少學(xué)者[9-10]應(yīng)用不同時(shí)間跨度、不同空間展布的大地測量資料對康定地震震后形變開展初步研究，但是該區(qū)域流變結(jié)構(gòu)仍然沒有得到很好的約束。因此，本文以2014年康定6.3級地震震后InSAR形變資料作為約束，采用三層粘彈性地震周期模型及遺傳算法反演鮮水河斷裂南東段康定地震震區(qū)深度介質(zhì)粘彈性層的厚度及粘滯系數(shù)，探索該區(qū)域深部流變結(jié)構(gòu)，為后續(xù)該地區(qū)的地震活動預(yù)測研究提供幫助。

1 數(shù)據(jù)與方法

1.1 數(shù)據(jù)介紹

本文采用2015～2019年兩個(gè)相鄰軌道的降軌Sentinel-1影像數(shù)據(jù)來獲取康定地震的震后形變，軌道號分別為33和135(表1)，東西向覆蓋約 460 km，南北向覆蓋約 480 km(圖1)。InSAR 數(shù)據(jù)處理采用 GAMMA 商業(yè)軟件平臺，引入空間分辨率為 30 m 的 SRTM 數(shù)字高程模型作為外部地形信息以提高配準(zhǔn)精度；利用兩軌法生成差分干涉圖，使用 DEM 模擬并消除地形相位；考慮到研究區(qū)域植被茂盛、氣候多變，InSAR 相干性差，需要對干涉圖進(jìn)行多視及濾波處理，提高相位相干性；然后進(jìn)行相位解纏、軌道精化及去大氣等處理；最后采用 π-RATE 軟件包[11]計(jì)算鮮水河斷裂帶兩個(gè)軌道的平均地殼形變速率場，并將雷達(dá)視線方向形變速率轉(zhuǎn)換到斷層平行方向，得到康定地震沿?cái)鄬悠叫蟹较虻恼鸷笮巫兯俣葓?圖2)。

表1 Sentinel-1數(shù)據(jù)參數(shù)

圖2 研究區(qū)InSAR震后形變速度場Fig.2 Post-seismic deformation velocity field of InSAR in the studied area

1.2 方法

與時(shí)間相關(guān)的震后形變可以用粘彈性地震周期模型來解釋，在該類模型中，地震周期同震階段產(chǎn)生的應(yīng)力在深度上通過粘彈性過程得到松弛，從而在整個(gè)地震周期的震間期產(chǎn)生與時(shí)間有關(guān)的地表速度。不同學(xué)者研究瑪尼震后形變所用的地震周期模型大多是在均質(zhì)半無限Maxwell流變粘彈性層上方有彈性上地殼的二層模型的基礎(chǔ)上構(gòu)建的。由于大地電磁和地震波數(shù)據(jù)表明鮮水河斷裂下方存在軟流層[12-13]，本文采用三層模型(彈性上地殼和彈性上地幔之間存在一個(gè)Maxwell粘彈性中地殼層)[6](圖3)構(gòu)建粘彈性地震周期模型，以約束康定地區(qū)下方中地殼層的厚度和粘度。為計(jì)算三層模型在整個(gè)地震周期內(nèi)的時(shí)變速度，DeVries等[14]發(fā)展了一個(gè)地震間變形隨時(shí)間演化的二維半解析解，以響應(yīng)限于中地殼層的粘彈性應(yīng)力松弛的無限序列周期地震。使用傳播矩陣方法計(jì)算位錯(cuò)引起的應(yīng)力和位移，同時(shí)要求層界面處的應(yīng)力連續(xù)，并在粘彈性三層模型的傅立葉域可以找到對單一地震響應(yīng)的解。

圖3 三層模型幾何示意圖Fig.3 Geometric sketch of three-layer model

由于時(shí)間序列的位移量是對能引起斷層偏移的一系列地震的響應(yīng)，假設(shè)具有復(fù)發(fā)間隔T的周期地震序列，基于麥克斯韋流變?yōu)榫€性的事實(shí)，可以求出單個(gè)地震的速度v1，從而得到N次前震的響應(yīng)：

(1)

式中，t是距離上次地震的時(shí)間。隨著地震次數(shù)N的增加，地震間速度分布接近一個(gè)周期穩(wěn)定狀態(tài)，在N趨于無窮大時(shí)，速度v表示為：

v=Δu(C1e-β1t+C2e-β2t)

(2)

式中， Δu為同震偏移，系數(shù)C1、C2和弛豫時(shí)間β1、β2表示粘彈性中地殼層厚度、彈性上地殼層厚度、同震破裂延伸的深度、Maxwell時(shí)間、地震復(fù)發(fā)周期和波數(shù)構(gòu)成的復(fù)雜函數(shù)，模型實(shí)現(xiàn)代碼已由DeVries等人公開，可在大地測量(GPS、InSAR)數(shù)據(jù)的約束下，采用遺傳算法反演出該地區(qū)巖石圈的流變結(jié)構(gòu)。

遺傳算法是由生物進(jìn)化論原理發(fā)展而來的搜索優(yōu)化算法，該算法有利于避免最終解陷入適應(yīng)度局部最大處，能搜索到全局最大處，因而被廣泛應(yīng)用于地學(xué)參數(shù)研究領(lǐng)域。本文使用遺傳算法計(jì)算時(shí)初始種群數(shù)為50，交叉概率為0.85，變異概率為0.2，目標(biāo)函數(shù)為模型值與觀測值的殘差加權(quán)平方和：

(3)

2 反演結(jié)果

為探索與康定地震震后速率相適應(yīng)的介質(zhì)參數(shù)，我們利用上述三層粘彈性周期模型進(jìn)行反演，其中粘滯系數(shù)的范圍設(shè)置為1×(1016～1019) Pa·s，粘彈性層的范圍設(shè)置為5～15 km。除了大地測量觀測值，此模型還需要地震復(fù)發(fā)周期、長期斷層滑動率和有效斷層閉鎖深度等參數(shù)。歷史地震、同震滑動和古地震的地表破裂調(diào)查結(jié)果表明，康定斷裂的地震平均復(fù)發(fā)周期約為200 a[15]?？刀ǖ貐^(qū)分為3個(gè)分支斷裂,分別是雅拉河斷裂、色拉哈斷裂和折多塘斷裂，單條斷裂的滑動速率小于10 mm/a，色拉哈斷裂滑動速率約為5.5 mm/a[1]，3條斷層滑動速率之和約為10mm/a[16]。因此，本文將地震復(fù)發(fā)周期設(shè)定為200 a，斷層滑動速率設(shè)定為10 mm/a。GPS和InSAR估算結(jié)果表明，康定段的閉鎖深度為10～15 km[17-19]，在此范圍內(nèi)，我們通過測試最終確定閉鎖深度為12 km。另外，通過分析不同彈性層厚度的擬合結(jié)果，設(shè)定彈性層厚度為12 km。

模型最佳模擬結(jié)果如圖4所示。圖4(a)顯示觀測到的地震后速度(藍(lán)色點(diǎn))和每個(gè)模型的最佳擬合速度曲線(紅色曲線為三層粘彈性周期模型，黑色曲線為彈性模型)，可以看出，與彈性模型相比，粘彈性模型能更好地解釋該區(qū)域的震后形變特征；灰色曲線表示彈性層不同厚度(5 km、10 km、15 km、20 km、25 km)的擬合曲線，結(jié)果表明，彈性層越薄，形變速率越大，但當(dāng)彈性層過薄時(shí)，該模型不能約束康定地震震后形變；紅色曲線顯示了當(dāng)彈性層厚度為12 km時(shí)的最優(yōu)結(jié)果，得到粘滯系數(shù)為9.9×1017Pa·s，粘彈性層厚度為5 km。結(jié)果表明，斷層附近的擬合效果較好，距離斷層越遠(yuǎn)擬合效果越差，而且斷層南西盤總體上比北東盤擬合效果要好。當(dāng)距離斷層大于30 km時(shí)，斷層南西盤速度曲線誤差明顯，擬合效果較差。斷層北東盤整體擬合效果不佳，這說明同一粘滯系數(shù)并不能同時(shí)解釋斷層兩側(cè)的形變特征。對其模型殘差(圖4(b))分析發(fā)現(xiàn)，南西盤距斷層10 km以內(nèi)和北東盤距斷層5 km以內(nèi)的數(shù)據(jù)擬合殘差較小，最小殘差約為0.01 mm/a；距斷層越遠(yuǎn)擬合殘差越大，最大殘差為5.64 mm/a，且斷層左側(cè)的擬合殘差遠(yuǎn)小于右側(cè)的擬合殘差。圖4(c)為遺傳算法收斂曲線，黑點(diǎn)為最優(yōu)值，藍(lán)色為平均值?？梢钥闯?，模型能夠快速收斂，在20代收斂于1.6 mm/a，說明反演模型的參數(shù)設(shè)置合理。

圖4 不同模型的模擬結(jié)果Fig.4 Simulation results of different models

鑒于同一粘滯系數(shù)不能同時(shí)解釋斷層兩側(cè)的形變特征，我們將兩側(cè)分別反演，結(jié)果如圖5所示，上圖(圖(a1)、(b1)、(c1))為斷層左側(cè)(南西盤)數(shù)據(jù)的擬合結(jié)果，下圖(圖(a2)、(b2)、(c2))為斷層右側(cè)(北東盤)數(shù)據(jù)的擬合結(jié)果。圖5(a)(藍(lán)色點(diǎn)表示參與模型計(jì)算的約束數(shù)據(jù)，紅色曲線為3層粘彈性周期模型最優(yōu)擬合曲線，灰色曲線表示彈性層不同厚度(5 km、10 km、15 km、20 km、25 km)的擬合曲線，黑色曲線為彈性模型擬合曲線)為觀測到的地震后速度和模型最佳擬合速度曲線。相比圖4(a)，兩側(cè)分開擬合效果更好，適應(yīng)康定震后形變約束的南西盤最佳粘滯系數(shù)為8.7×1017Pa·s，其相對應(yīng)的粘彈性層厚度為5 km，距斷層3～30 km范圍內(nèi)擬合效果較好，小于3 km和大于30 km時(shí)擬合效果較差(圖(a1))；北東盤最佳粘滯系數(shù)為1.2×1018Pa·s，其相對應(yīng)的粘彈性層厚度為6 km，距斷層小于15 km時(shí)擬合效果不佳，可能是因?yàn)樾巫兯俾首兓^快不能很好地約束模型，而大于15 km時(shí)擬合效果明顯優(yōu)于整體擬合(圖(a2))。圖5(b)為模型擬合殘差，南西盤距斷層過近或太遠(yuǎn)時(shí)擬合殘差較大，但擬合殘差平均值小于1 mm/a(圖(b1))；北東盤距斷層稍遠(yuǎn)(大于30 km)時(shí)，擬合殘差較小，擬合殘差平均值同樣小于1 mm/a(圖(b2))。圖5(c)為遺傳算法收斂曲線，黑點(diǎn)為最優(yōu)值，藍(lán)色為平均值?？梢钥闯觯Ｐ途軌蚩焖偈諗?，分別收斂于1.2 mm/a和1.9 mm/a，說明反演模型的參數(shù)設(shè)置合理。

圖5 斷層兩側(cè)分別反演的最佳模擬結(jié)果Fig.5 Optimal simulation results on both sides of fault

3 討 論

3.1 震后形變特征分析

為了弄清該三層模型中不同參數(shù)對震后形變速度的影響程度及方式，我們嘗試每個(gè)參數(shù)不同值的正演，結(jié)果如圖6所示。假設(shè)斷層以10 mm/a的速度滑動，地震復(fù)發(fā)周期為100 a。當(dāng)彈性上地殼厚度H=20 km、粘彈性下地殼厚度h=10 km、粘滯系數(shù)eta=1×1018Pa·s時(shí)，隨著有效閉鎖深度D從1 km變化到20 km，模型產(chǎn)生震后變形的速度越來越快，速度變化量越來越大(圖6(a))；將有效閉鎖深度D固定為20 km，將彈性上地殼厚度H從20 km變化到60 km，模型產(chǎn)生震后變形的速度越來越慢，速度變化量越來越小(圖6(b))；同理，將粘彈性下地殼厚度h從1 km變化到40 km，模型產(chǎn)生震后變形的速度越來越快，速度變化量并不明顯(圖6(c))；將粘滯系數(shù)從 1×1016Pa·s變化到1×1020Pa·s，模型產(chǎn)生震后變形的速度越來越快，但速度變化量越來越小，很快達(dá)到最大速度值(圖6(d))。實(shí)驗(yàn)表明，隨著有效閉鎖深度增大，形變速率越來越大；上地殼彈性層厚度越大，形變速率越小；下地殼粘彈性層厚度越大，形變速率越大；下地殼粘彈性層粘滯系數(shù)越小越容易引起變形。

圖6 不同參數(shù)調(diào)試的斷層平行速度剖面Fig.6 Parallel velocity profiles of faults adjusted by different parameters

3.2 斷層兩側(cè)介質(zhì)差異

鮮水河斷裂帶位于川滇、松潘-甘孜和揚(yáng)子3個(gè)地塊的交界區(qū)，受新生代以來印度板塊與歐亞板塊之間持續(xù)的消減和碰撞作用，早期構(gòu)造活動可追溯到47～42 Ma，并于18～12 Ma分階段性地侵入折多山巖體，在巖漿侵位后，于5.5～3 Ma階段發(fā)生韌性左行剪切走滑[20-21]。因此，研究區(qū)斷層南西盤以三疊系巖石類型為主的區(qū)域出現(xiàn)了新近-古近系侵入巖，而斷層北東盤地層巖性較為復(fù)雜，從元古到新生代均有出露，以元古界、泥盆系、石巖-二疊系、二疊紀(jì)侵入巖以及三疊紀(jì)等巖石類型為主。本文以康定地震震后形變速度為約束，反演鮮水河斷裂南東段下地殼粘滯系數(shù)，研究表明，南西盤最佳粘滯系數(shù)(8.7×1017Pa·s)小于北東盤最佳粘滯系數(shù)(1.2×1018Pa·s)，證明鮮水河斷裂帶南東段兩側(cè)上地殼物質(zhì)存在顯著的橫向介質(zhì)差異。結(jié)合研究區(qū)的深部速度結(jié)構(gòu)特征[4, 13]發(fā)現(xiàn)，鮮水河斷裂帶南東段東側(cè)呈現(xiàn)出高速異常，西側(cè)呈現(xiàn)出低速異常，表明鮮水河斷裂帶南東段兩側(cè)上地殼物質(zhì)存在顯著的橫向介質(zhì)差異，而康定MS6.3地震就發(fā)生在該高、低速異常區(qū)的分界線上。

4 結(jié) 語

1)數(shù)值實(shí)驗(yàn)結(jié)果表明，隨著有效閉鎖深度增大，形變速率越來越大；上地殼彈性層厚度越大，形變速率越?。幌碌貧ふ硰椥詫雍穸仍酱?，形變速率越大；下地殼粘彈性層粘滯系數(shù)越小越容易引起變形。

2)以康定震后InSAR形變資料作為約束，采用三層粘彈性地震周期模型及遺傳算法反演鮮水河斷裂南東段的康定地震震源區(qū)深度介質(zhì)粘彈性層的厚度及粘滯系數(shù)。結(jié)合鮮水河斷裂區(qū)域地質(zhì)構(gòu)造背景和地球物理觀測表明，康定區(qū)域下地殼粘彈性層的厚度為5 km，粘滯系數(shù)為9.9×1017Pa·s，且斷層兩側(cè)下地殼流變存在顯著的橫向不均一性；將斷層兩側(cè)分別反演得出，斷層南西盤下地殼粘滯系數(shù)(8.7×1017Pa·s)略小于斷層北東盤下地殼粘滯系數(shù)(1.2×1018Pa·s)。