臧 陽,韓顏顏,孟令媛
(1.中國地震局地球物理研究所,北京 100081;2.中國地震臺網(wǎng)中心,北京 100045)
緬甸弧俯沖帶是印度板塊與緬甸塊體相互作用的地區(qū),緬甸弧的形成以及周邊兩大板塊間的運(yùn)動和動力學(xué)模式歷來是地學(xué)工作者的研究熱點(diǎn)。Meade(2007)通過塊體模型模擬了印度與亞洲碰撞區(qū)的運(yùn)動模式,認(rèn)為印度板塊在NNE向運(yùn)動的同時(shí),呈現(xiàn)自西向東的俯沖;Liu和Bird(2008)利用大地測量結(jié)果、斷層滑動速率和區(qū)域構(gòu)造應(yīng)力場建立了青藏高原東南部的運(yùn)動與動力模式;張浪平等(2013)以地震活動參數(shù)為基礎(chǔ)數(shù)據(jù),準(zhǔn)定量地給出了緬甸弧區(qū)域板塊間接觸的三維幾何形態(tài)。前人的研究結(jié)果表明:印度板塊向緬甸塊體東側(cè)強(qiáng)烈擠壓俯沖作用導(dǎo)致了中源地震的頻發(fā),而東側(cè)縫合帶即實(shí)皆斷裂的水平走滑運(yùn)動則是緬甸弧大部分淺源強(qiáng)震發(fā)生的主要原因(Maurinetal,2010;Hurukawa,Maung,2011)。
印度板塊與緬甸塊體的相互作用不僅在緬甸弧俯沖帶附近表現(xiàn)出劇烈的強(qiáng)震活動,其對我國川滇地區(qū)的構(gòu)造環(huán)境和強(qiáng)震活動同樣有著深遠(yuǎn)的影響。川滇地區(qū)位于印度板塊向緬甸塊體俯沖的前緣,是印度板塊東向擠壓作用直接影響的地區(qū),同時(shí)受到青藏高原物質(zhì)東向擠出作用的共同影響,構(gòu)造上表現(xiàn)為一系列近似平行于緬甸弧的斷裂,如鮮水河斷裂、紅河斷裂等(蘇有錦,秦嘉政,2001)。胡家富等(2008)聯(lián)合地震活動與速度結(jié)構(gòu)分析認(rèn)為:印度板塊從94°E附近開始俯沖,且縫合于實(shí)皆斷裂,并且給出了印度板塊向東俯沖的動力學(xué)模型,指出我國滇西附近地區(qū)是板塊俯沖作用直接影響的區(qū)域。地震活動性方面,由于川滇地區(qū)的地形、地貌和地質(zhì)構(gòu)造格局均受到印度板塊俯沖作用的影響和控制,因此2個(gè)區(qū)域的地震活動存在密切的相關(guān)性,劉建中(1992)認(rèn)為川滇地區(qū)地震活躍期的開始與結(jié)束較緬甸弧附近地區(qū)相應(yīng)的地震活躍期滯后,滯后時(shí)間為2~3年。
研究緬甸弧地區(qū)強(qiáng)震對川滇地區(qū)主要斷裂的應(yīng)力加載作用對于深入理解緬甸弧強(qiáng)震觸發(fā)機(jī)制及其對川滇地區(qū)的應(yīng)力影響具有重要意義。該研究現(xiàn)階段仍以同震靜態(tài)庫侖應(yīng)力觸發(fā)為主要手段(孟令媛等,2015),缺少對震后累積應(yīng)力作用的系統(tǒng)研究。本文在前人研究的基礎(chǔ)上,選取Burgers體構(gòu)建粘彈性介質(zhì)模型,分別研究緬甸弧俯沖帶附近中、淺源強(qiáng)震在川滇地區(qū)主要斷裂產(chǎn)生的同震及震后累積庫侖應(yīng)力,嘗試通過緬甸弧地區(qū)震源機(jī)制解特征給出該地區(qū)中源地震可能的發(fā)震機(jī)制,并通過分析累積庫侖應(yīng)力計(jì)算結(jié)果,系統(tǒng)地描述緬甸弧地區(qū)強(qiáng)震對川滇地區(qū)的應(yīng)力影響。
本文考慮緬甸弧地區(qū)強(qiáng)震應(yīng)力觸發(fā)作用對川滇地區(qū)的影響,建立粘彈性應(yīng)力觸發(fā)模型,該模型考慮到了上地幔物質(zhì)對地震能量的運(yùn)移作用。由于存在粘彈性松弛效應(yīng),這種應(yīng)力傳遞作用存在時(shí)間上的延遲,可以適用于研究幾百年尺度的長期應(yīng)力累積影響(徐晶等,2013;孟令媛等,2016)。
地震的發(fā)生主要是由脆性巖石破裂后的摩擦運(yùn)動造成的,假設(shè)地球介質(zhì)為半無限空間,且空間內(nèi)各向同性,若已知斷層位錯(cuò)的幾何參數(shù)和錯(cuò)動大小,可求得斷層面上庫侖破裂應(yīng)力變化為:
ΔCFS=Δτ+f(Δσn+Δp)
(1)
式中:Δτ為剪切應(yīng)力變化;Δσn為正應(yīng)力變化;f為斷層面介質(zhì)摩擦系數(shù);Δp為孔隙壓力變化,調(diào)節(jié)著斷層面上的有效正應(yīng)力。Δp造成的摩擦系數(shù)的減少可用f′=f(1-B)表示,其中B為Skemptons系數(shù),取值范圍為0~1。當(dāng)巖石應(yīng)力改變遠(yuǎn)遠(yuǎn)快于巖石中Δp的擴(kuò)散時(shí),Δp可以通過B值在式(1)中得到體現(xiàn),且轉(zhuǎn)換為:
ΔCFS=Δτ+f′Δσn
(2)
式中:f′包含了孔隙流體和斷層面上介質(zhì)的特性,一般取值為0.2~0.8(Kingetal,1994),本文取值為0.7。
本文收集了1969—2018年緬甸弧附近、緬甸弧俯沖帶及東側(cè)實(shí)皆斷裂附近地區(qū),在全球矩心矩張量工程(Global Centroid Moment Tensor Project,簡稱GCMT,)有震源機(jī)制解記錄的所有淺源及中源地震。包括MW4.0~4.9地震24次(淺源12次,中源12次),MW5.0~5.9地震127次(淺源72次,中源55次),MW6.0~6.9地震17次(淺源10次,中源7次),MW7.0~7.9地震2次(淺源0次,中源2次),共170次(淺源94次,中源76次)地震??傮w地震震級范圍為MW4.7~7.2,震源深度范圍為6~157 km(圖1),中源地震數(shù)量從70~100 km隨震源深度的增加逐漸增多,在100~110 km深度達(dá)到峰值,隨后數(shù)量逐漸減少。其中,MW≥6.0地震信息列于表1。
表1 1969—2018年緬甸弧地區(qū)MW≥6.0地震信息
圖1 1969—2018年緬甸弧地區(qū)MW≥4.7地震震源深度分布圖Fig.1 Distribution of MW≥4.7 events along depth near Burma arc from 1969 to 2018
以Crust1.0模型為基礎(chǔ)計(jì)算得到川滇地區(qū)上、中、下地殼的平均厚度分別為20,15和10 km??紤]地殼及上地幔介質(zhì)的滯彈性,根據(jù)楊強(qiáng)和黨亞民(2010)利用GPS速度場估算的青藏高原地殼韌性層等效粘滯系數(shù)分布的研究,青藏高原及周邊區(qū)域脆性-韌性轉(zhuǎn)換面一般位于中地殼,等效粘滯系數(shù)范圍在中地殼為1019~1022Pa·s,在下地殼為1017~1020Pa·s,在Moho面則降至1016~1018Pa·s,該結(jié)果與劉代志(1992)、Flesch等(2005)、石耀霖和曹建玲(2008)以及張晁軍等(2008)的研究基本一致。考慮到粘滯松弛存在短期和長期效應(yīng),采用Burgers體模型,取楊強(qiáng)和黨亞民(2010)研究結(jié)果中各層粘滯系數(shù)范圍的低值賦值為η1,高值賦值為η2,并結(jié)合Crust1.0的速度和密度模型,最終建立研究區(qū)簡化分層模型,如表2所示。
表2 研究區(qū)介質(zhì)模型參數(shù)
本文以鄧起東(2007)匯總的中國活動構(gòu)造數(shù)據(jù)作為基礎(chǔ),將川滇地區(qū)共51條主要斷裂、129個(gè)斷層段、2 112個(gè)斷層離散點(diǎn)作為接收斷層模型,使用Wang等(2006)的PSGRN/PSCMP程序,計(jì)算格林函數(shù)以及水平分層粘彈性介質(zhì)中同震和震后庫侖應(yīng)力。由于川滇地區(qū)強(qiáng)震多為淺源地震,本文主要考慮緬甸弧強(qiáng)震對川滇地區(qū)主要活動斷裂淺部的應(yīng)力影響,在計(jì)算格林函數(shù)的過程中,將觀測深度統(tǒng)一設(shè)置為地下10 km。
研究表明:中、深源地震的成因明顯區(qū)別于淺源地震,即隨著震源深度的增加,地下溫度和壓力不斷增加使得脆性和摩擦行為被強(qiáng)烈抑制,巖石表現(xiàn)為可抑制地震的韌性行為(David,John,1960)。中源地震的誘發(fā)機(jī)制主要是脫水或與流體相關(guān)的致脆以及塑性剪切失穩(wěn)(邵同賓,嵇少丞,2015),而其破裂機(jī)制是利用傳統(tǒng)巖石破裂模型計(jì)算中源地震庫侖應(yīng)力的理論基礎(chǔ)。
緬甸弧地區(qū)震源機(jī)制解資料顯示(圖2),沿其東側(cè)實(shí)皆斷裂主要分布著一系列淺源地震,其中實(shí)皆斷裂中段以純走滑型地震為主,與實(shí)皆斷裂右旋走滑的運(yùn)動性質(zhì)相符合;實(shí)皆斷裂北側(cè)和南側(cè)則分布著一系列逆沖型地震,震源機(jī)制解P軸方向整體為NE向,反映了實(shí)皆斷裂南北兩側(cè)的應(yīng)力環(huán)境以NE向的擠壓為主。緬甸弧俯沖帶西側(cè)分布著一系列淺源地震,震源機(jī)制解類型較為復(fù)雜,總體上以21.5°N附近為分割線,其北側(cè)震源機(jī)制以走滑為主,南側(cè)以正斷為主,反映了該區(qū)域淺部較為復(fù)雜的應(yīng)力環(huán)境(圖3)。緬甸弧地區(qū)中源地震震源機(jī)制以逆沖型為主,隨著震源深度的增加,這一特征更為明顯。
圖2 緬甸弧地區(qū)MW≥5.0(a)、MW≥6.0(b)地震震源機(jī)制解
圖3 緬甸弧地區(qū)MW≥5.0地震震源機(jī)制解隨震源深度和緯度分布
值得關(guān)注的現(xiàn)象是,緬甸弧地區(qū)中源地震震源機(jī)制解B軸方向沿緯度具有一定的分布規(guī)律,北側(cè)B軸在水平面上投影方向以NWW-SEE方向?yàn)橹?,南?cè)轉(zhuǎn)為SWW-NEE方向。假設(shè)緬甸弧板塊俯沖方向垂直于俯沖帶邊界,且整體自西向東俯沖,單元俯沖板片俯沖方向在水平面上的投影見圖4a。由圖可見,緬甸弧俯沖邊界整體呈半圓弧狀,向東側(cè)彎曲,在24°N以北區(qū)域(圖4a中粉色區(qū)域),板塊整體俯沖方向?yàn)镹W-SE向,在23°N以南區(qū)域(圖4a中藍(lán)色區(qū)域)為SW-NE向,23~24°N區(qū)域(圖4a中綠色區(qū)域)為近W-E向。由于俯沖邊界的弧形彎曲,致使相鄰單元板片俯沖方向出現(xiàn)多處相交,以23°~24°N區(qū)域最為顯著。假設(shè)俯沖到70~160 km以下的板片仍為脆性巖石結(jié)構(gòu),隨著俯沖深度的增加,俯沖方向所在延長線相交,意味著相鄰單元板片之間在垂直于俯沖方向上的相互擠壓作用力逐漸增強(qiáng),最終可能導(dǎo)致脆性板片內(nèi)部形成以俯沖方向?yàn)樽呦虻哪鏇_型破裂,理論上逆沖地震震源機(jī)制解B軸方向或斷層面走向應(yīng)與俯沖方向一致。
出現(xiàn)上述破裂機(jī)制的可能性隨著俯沖深度的增加不斷增大,以位于緬甸弧俯沖帶上震源深度≥100 km的所有MW≥5.0地震為樣本,對斷裂面走向進(jìn)行分區(qū)域統(tǒng)計(jì)。如圖4b所示,在24°N以北區(qū)域,地震破裂面走向以SE向?yàn)橹?,?3°N以南區(qū)域則以SW向?yàn)橹鲗?dǎo),在23~24°N區(qū)域則多為近EW方向。該統(tǒng)計(jì)結(jié)果表明:斷層面走向與俯沖方向存在高度一致性,支持上述關(guān)于俯沖帶中源地震破裂機(jī)制的假設(shè),這些中源地震仍為脆性板片內(nèi)部巖石的擠壓破裂,在計(jì)算庫侖應(yīng)力時(shí)可沿用常規(guī)淺源地震的計(jì)算方法。
圖4 緬甸弧俯沖帶邊界及單元板片俯沖方向在水平面的投影(a)和不同緯度范圍中源地震斷裂面走向比例分布圖(b)
另外,在震源稍淺的區(qū)域還發(fā)育有P軸方向與俯沖方向一致的逆沖型地震,傾角一般較小,這些地震與發(fā)生在俯沖帶的大多數(shù)淺源地震類似,可能是由上馱板塊與下伏板片之間摩擦運(yùn)動所導(dǎo)致,如2011年3月11日日本MW9.0地震(李圣強(qiáng)等,2013)以及2015年4月25日尼泊爾MW7.8地震(劉靜等,2015),主震震源機(jī)制解均表現(xiàn)出相似的特征,盡管2次地震分別為洋-陸碰撞和陸-陸碰撞模式形成的大型逆沖型地震,但二者的震源機(jī)制解均揭示了俯沖板片克服上馱板塊一側(cè)的擠壓力導(dǎo)致地震發(fā)生的物理機(jī)制。
通過理論模型實(shí)驗(yàn)考察不同類型單次強(qiáng)震對川滇地區(qū)主要斷裂的影響。實(shí)驗(yàn)需要建立理論地震模型,計(jì)算川滇地區(qū)主要斷裂的同震庫侖應(yīng)力大小。計(jì)算時(shí),地震震級統(tǒng)一取為MW7.0,破裂尺寸及位移距離根據(jù)Wells和Coppersmith(1994)研究中的經(jīng)驗(yàn)公式,在已知矩震級和震源機(jī)制解的前提下進(jìn)行換算,破裂方式均假設(shè)為雙側(cè)均勻破裂。
根據(jù)前文對緬甸弧地區(qū)強(qiáng)震震源機(jī)制解特征以及推測的中源地震發(fā)震機(jī)制的分析,首先考察典型中源強(qiáng)震對川滇地區(qū)的影響,震源模型經(jīng)緯度、深度以及震源機(jī)制解均參考真實(shí)發(fā)生在緬甸弧的MW≥6.0地震。測試的中源地震模型分別為:B軸沿俯沖方向的逆沖型地震,P軸沿俯沖方向的逆沖型地震其參數(shù)見表3。同震庫侖應(yīng)力計(jì)算結(jié)果如圖5a,b所示,2次中源地震對川滇地區(qū)的應(yīng)力加載作用量級均較小,最大僅為0.001~0.01 bar,其中1號地震模型應(yīng)力加載區(qū)主要集中在小滇西地區(qū)(圖5a),2號地震模型集中在滇西南地區(qū)(圖5b),由于中源地震模型具有較深的震源深度,因此對川滇地區(qū)的應(yīng)力加載作用十分有限。
實(shí)皆斷裂運(yùn)動性質(zhì)較明確,區(qū)域構(gòu)造環(huán)境相對簡單,其上發(fā)生的淺源強(qiáng)震震源機(jī)制較為穩(wěn)定,且與川滇地區(qū)距離最近。選取發(fā)生在實(shí)皆斷裂中段的右旋走滑型淺源地震以及實(shí)皆斷裂北段逆沖型淺源地震(表3)作為地震模型,計(jì)算得到川滇地區(qū)同震庫侖應(yīng)力如圖5c,d所示。由于震源深度的差異,且與川滇地區(qū)距離更近,3,4號地震模型對川滇地區(qū)的應(yīng)力加載作用明顯較中源地震強(qiáng),最高值達(dá)0.01~0.1 bar,其中3號地震模型對云南西部地區(qū)存在明顯加載作用,庫侖應(yīng)力最高值位于小滇西附近(圖5c),4號地震模型對小滇西大部分區(qū)域則為應(yīng)力卸載,但對滇東北地區(qū)至鮮水河斷裂帶均有較顯著的應(yīng)力加載作用,且明顯強(qiáng)于3號地震模型(圖5d)。
表3 地震模型震源參數(shù)
圖5 1號(a),2號(b),3號(c)和4號(d)地震模型對川滇地區(qū)主要斷裂同震庫侖應(yīng)力計(jì)算結(jié)果
考慮1969—2018年緬甸弧地區(qū)MW≥6.0的9次中源地震對川滇地區(qū)的影響,從1969年10月17日第1次中源地震開始(表1),以每一次中源地震的發(fā)震時(shí)刻作為時(shí)間切片計(jì)算該時(shí)刻在川滇地區(qū)累積的庫侖應(yīng)力,終止時(shí)間為2018年12月31日。如圖6所示,1988年8月6日MW7.2地震由于震級大且震源深度相對較淺(92 km),震后造成川滇地區(qū)應(yīng)力出現(xiàn)大幅變化,最大累積庫侖應(yīng)力達(dá)0.01~0.1 bar,位于小滇西—滇西北附近地區(qū)(圖6c)。該地震發(fā)生后,緬甸弧地區(qū)陸續(xù)發(fā)生多次MW≥6.0中源強(qiáng)震,但由于震級均小于MW7.0,單次地震影響有限(參考1,2號模型地震計(jì)算結(jié)果),隨著時(shí)間的推移,川滇地區(qū)累積庫侖應(yīng)力整體呈現(xiàn)減小趨勢。截至2018年底,最大累積庫侖應(yīng)力僅在小滇西和滇西北局部地區(qū)達(dá)到0.001~0.01 bar,其余地區(qū)均小于0.001 bar(圖6d~j)。因此,僅當(dāng)在緬甸弧地區(qū)發(fā)生的中源地震的震級較強(qiáng)(MW≥7.0),且震源深度較淺(<100 km時(shí)),才可能在川滇地區(qū)產(chǎn)生一定程度的應(yīng)力加載,一般情況下,緬甸弧地區(qū)中源強(qiáng)震的發(fā)生在川滇地區(qū)累積的庫侖應(yīng)力量級不會超過0.01 bar。
圖6 緬甸弧地區(qū)中源強(qiáng)震對川滇地區(qū)主要斷裂累積庫侖應(yīng)力時(shí)間進(jìn)程(以9次中源地震發(fā)生后第1天及2018年12月31日作為時(shí)間切片)
圖7為1969—2018年緬甸弧地區(qū)淺源地震在川滇地區(qū)累積的庫侖應(yīng)力計(jì)算結(jié)果,計(jì)算方法同中源地震。盡管緬甸弧地區(qū)淺源地震震級整體低于中源地震(表1),但由于距離優(yōu)勢,對川滇地區(qū)的影響明顯強(qiáng)于中源地震。1992年6月15日一次MW6.3地震就能夠造成小滇西地區(qū)累計(jì)庫侖應(yīng)力的明顯增強(qiáng)(圖7d),隨后1994年5月29日MW6.5地震的發(fā)生,使小滇西地區(qū)最大累積庫侖應(yīng)力達(dá)到0.01~0.1 bar(圖7e)。此后緬甸弧地區(qū)發(fā)生的淺源地震震級均≤MW6.8,但與中源地震結(jié)果相反,淺源地震所造成的累積庫侖應(yīng)力在不斷增強(qiáng),截至2018年12月31日,盡管最大累積庫侖應(yīng)力始終未能超過0.1 bar,但小滇西大范圍地區(qū)庫侖應(yīng)力大小達(dá)到0.01~0.1 bar。
圖7 緬甸弧地區(qū)淺源強(qiáng)震對川滇地區(qū)主要斷裂累積庫侖應(yīng)力時(shí)間進(jìn)程(以9次淺源地震發(fā)生后第1天及2018年12月31日作為時(shí)間切片)
圖8為1969—2018年緬甸弧地區(qū)19次MW≥6.0地震在川滇地區(qū)的累積庫侖應(yīng)力,由于中源地震的影響較小,截至2018年12月31日與單純考慮淺源強(qiáng)震的結(jié)果基本一致,說明從應(yīng)力觸發(fā)的角度來看,緬甸弧地區(qū)強(qiáng)震對川滇地區(qū)主要斷裂的應(yīng)力加載作用以淺源地震占主導(dǎo),中源地震除震級≥MW7.0,且震源深度較淺等特殊情況外,在川滇地區(qū)所造成的庫侖應(yīng)力變化一般可忽略不計(jì)。
圖8 1969—2018年緬甸弧地區(qū)所有MW≥6.0強(qiáng)震對川滇地區(qū)主要斷裂累積庫侖應(yīng)力計(jì)算結(jié)果(截至2018年12月31日)
根據(jù)中國地震臺網(wǎng)中心全球7級強(qiáng)震目錄,1900—2000年緬甸弧俯沖帶共發(fā)生5次MS≥7.0地震,統(tǒng)計(jì)川滇地區(qū)在后續(xù)1年以及1~2年內(nèi)發(fā)生的所有MS≥6.0地震,發(fā)現(xiàn)這些地震的震中位置在空間上形成明顯分界:震后1年內(nèi)地震分布在滇西南—小滇西—川滇藏交界附近,距緬甸弧相對較近;而震后1~2年地震分布在龍門山斷裂—川滇交界中部—滇南附近地區(qū),距緬甸弧相對較遠(yuǎn)(圖9)。緬甸弧地區(qū)中源強(qiáng)震的發(fā)生,對川滇地區(qū)的地震活動性似乎存在一定影響,且隨時(shí)間的推移影響的距離逐漸增大,因此,從地震活動性的統(tǒng)計(jì)學(xué)角度來看,緬甸弧中源地震與川滇地區(qū)強(qiáng)震在時(shí)空上具有明顯相關(guān)性。另外,川滇地區(qū)在緬甸弧地區(qū)中源強(qiáng)震后1年內(nèi)發(fā)生的MS≥6.0地震的空間位置,與圖5a,b中累積庫侖應(yīng)力超過0.001 bar的區(qū)域具有一定的一致性,可能反映了中源地震的應(yīng)力影響在其中發(fā)揮著一定的作用。
圖9 1900—2018年緬甸弧附近MS≥7.0中源地震與后續(xù)2年內(nèi)川滇地區(qū)MS≥6.0地震統(tǒng)計(jì)關(guān)系Fig.9 Statistical relationship between intermediateMS≥7.0 earthquakes near Burma arc and theMS≥6.0 earthquakes in the Sichuan-Yunnan region occurred within the subsequent 2 years from 1900 to 2018
本文收集了1969—2018年緬甸弧附近地區(qū)共170次(淺源94次、中源76次)地震的震源機(jī)制解,通過分析這些地震的空間位置及震源機(jī)制解特征,得到以下結(jié)論:
(1)緬甸弧東側(cè)實(shí)皆斷裂中段,地震震源機(jī)制以右旋走滑型為主,與斷裂運(yùn)動性質(zhì)相符,實(shí)皆斷裂南北兩側(cè)多發(fā)逆沖型地震,主壓應(yīng)力方向?yàn)榻麼E向。
(2)緬甸弧中部多為發(fā)生在俯沖板片上的逆沖型中源地震,震源深度較深時(shí),震源機(jī)制解B軸方向與單元板片俯沖方向一致,可能是相鄰單元板片之間相互擠壓導(dǎo)致的脆性破裂;震源深度較淺時(shí),部分地震P軸方向與俯沖方向一致,且傾角較小,可能反映了下伏板片與上馱板塊之間的摩擦運(yùn)動。
(3)緬甸弧西側(cè)淺源地震震源機(jī)制解以走滑及正斷型為主,總體上以21.5°N附近為分割線,其北側(cè)走滑型占主導(dǎo),南側(cè)正斷型占主導(dǎo),反映了該區(qū)淺部較為復(fù)雜的應(yīng)力環(huán)境。
在對緬甸弧地區(qū)中源地震發(fā)震機(jī)制進(jìn)行合理假設(shè)的基礎(chǔ)上,通過建立緬甸弧至川滇地區(qū)粘彈性介質(zhì)模型,以川滇地區(qū)51條主要斷裂、129個(gè)斷層段、2 112個(gè)斷層離散點(diǎn)作為接收斷層模型,分別計(jì)算了單次地震模型,緬甸弧地區(qū)中源、淺源及19次MW≥6.0地震對川滇地區(qū)的累積庫侖應(yīng)力,根據(jù)計(jì)算結(jié)果得到以下結(jié)論:
(1)緬甸弧地區(qū)單次中源地震對川滇地區(qū)應(yīng)力加載作用主要集中在小滇西及滇西南地區(qū),MW7.0地震引起的最大同震庫侖應(yīng)力范圍為0.001~0.01 bar。
(2)實(shí)皆斷裂單次淺源地震對川滇地區(qū)的應(yīng)力加載作用主要集中在小滇西及滇西北地區(qū),MW7.0地震引起的最大同震庫侖應(yīng)力范圍為0.01~0.1 bar。
(3)緬甸弧地區(qū)中源地震對川滇地區(qū)累積庫侖應(yīng)力影響較小,僅當(dāng)?shù)卣鹫鸺壿^強(qiáng)(MW≥7.0),且震源深度較淺(<100 km)時(shí),才可能產(chǎn)生一定應(yīng)力加載作用,一般情況下,中源地震在川滇地區(qū)的累積庫侖應(yīng)力不超過0.01 bar。
(4)緬甸弧地區(qū)淺源地震在川滇地區(qū)累積的庫侖應(yīng)力明顯強(qiáng)于中源地震,從應(yīng)力觸發(fā)角度來看,緬甸弧地區(qū)強(qiáng)震對川滇地區(qū)的應(yīng)力加載作用以淺源地震占主導(dǎo),中源地震一般可忽略不計(jì)。
(5)截至2018年12月31日,1969—2018年緬甸弧地區(qū)所有MW≥6.0地震對川滇地區(qū)的最大累積庫侖應(yīng)力小于0.1 bar,應(yīng)力加載作用最強(qiáng)的區(qū)域集中在小滇西—滇西北附近地區(qū)。
盡管緬甸弧地區(qū)中源強(qiáng)震在統(tǒng)計(jì)上對后續(xù)川滇地區(qū)地震活動具有某種程度的影響,但正如前面所提到的,其對川滇地區(qū)的應(yīng)力加載作用非常有限,一般不會超過0.01 bar。根據(jù)以往研究,當(dāng)庫侖應(yīng)力超過0.1 bar時(shí),才會對斷裂具有明顯的加載作用,引發(fā)強(qiáng)震(Stein,Lisowski,1983),而即使考慮1969—2018年所有淺源和中源地震的影響,對川滇地區(qū)的累積庫侖應(yīng)力最大值依然小于0.1 bar。因此,緬甸弧地區(qū)中源強(qiáng)震的應(yīng)力加載作用,可能是后續(xù)在川滇地區(qū)地震活動的影響因素之一,當(dāng)同時(shí)滿足較高的背景應(yīng)力等其他因素時(shí)才會導(dǎo)致地震的最終發(fā)生。另外,中源地震的發(fā)生可能反應(yīng)了更大尺度上印度板塊向歐亞板塊推擠作用的增強(qiáng),這可能是促進(jìn)川滇地區(qū)后續(xù)強(qiáng)震活動的一個(gè)重要因素。對比緬甸弧地區(qū)中源、淺源強(qiáng)震對川滇地區(qū)的應(yīng)力影響,表明淺源地震要顯著強(qiáng)于中源地震,因此在關(guān)注緬甸弧地區(qū)中源地震對川滇地區(qū)影響的同時(shí),還應(yīng)對淺源強(qiáng)震給予足夠的重視。