曹穎， 黃江培* ， 錢佳威， 付虹

1 云南省地震局， 昆明 650224 2 中國科學(xué)技術(shù)大學(xué)地球和空間科學(xué)學(xué)院， 合肥 230026

0 引言

據(jù)中國地震臺(tái)網(wǎng)測(cè)定，2014年8月3日16∶ 30在云南省昭通市魯?shù)榭h發(fā)生了MS6.5地震(震中103.3°E，27.1°N，震源深度12 km)(CENC，http:∥www.cenc.ac.cn).此次地震震級(jí)雖不大，但地震極震區(qū)最高烈度可達(dá)Ⅸ度，且由于震源淺、震區(qū)人口密度高，是云南省近年來典型的“小震”致大災(zāi)事件(李西等，2018)，引起國內(nèi)外的高度關(guān)注.魯?shù)榈卣鸬难芯拷Y(jié)果，讓我們對(duì)魯?shù)榈卣鸬恼鹪刺匦?、運(yùn)動(dòng)學(xué)特征及發(fā)震構(gòu)造等有了深刻的認(rèn)識(shí).重定位結(jié)果顯示，魯?shù)榈卣鹩嗾鹦蛄械恼鹬蟹植贾饕獮榻麰W向和SES向的共軛型分布，震源深度主要分布于13～15 km，主震位于破裂底部(張廣偉等，2014；房立華等，2014；王光明等，2018)；結(jié)合震源機(jī)制解和震源破裂過程研究(張廣偉等，2014；徐錫偉等，2014；張勇等，2014，2015；許力生等，2014)，可知此次地震為一次高傾角左旋走滑型地震，是在NW向主壓應(yīng)力與NE向主張應(yīng)力相互作用下發(fā)生的以NNW向破裂為主的共軛破裂事件.此次地震周邊地質(zhì)構(gòu)造復(fù)雜，大多數(shù)研究(徐錫偉等，2014；李西等，2014，2018；房立華等，2014；張勇等，2014，2015；程佳等，2016，付芮等，2015)在結(jié)合了野外考察、余震序列分布、震源機(jī)制解、地震烈度分布、庫倫應(yīng)力變化和震源破裂過程等研究結(jié)果后認(rèn)為發(fā)震斷裂為NNW向的包谷垴—小河斷裂.以往也有對(duì)魯?shù)榈卣鸢l(fā)生前后地下速度結(jié)構(gòu)的研究，Riaz等(2017)基于雙差層析成像方法給出了魯?shù)榈卣鹬苓厖^(qū)域的地殼層析成像結(jié)果，發(fā)現(xiàn)魯?shù)榈卣鹫鹪磪^(qū)低P波速度異常被高P波速度區(qū)域所包圍，大多數(shù)地震發(fā)生在低P波速度和低泊松比區(qū)域，受主震破裂的限制.趙小艷和孫楠(2014)利用震源位置和速度結(jié)構(gòu)的聯(lián)合反演得出了魯?shù)榈卣鹩嗾鹬饕植荚赑波速度較高的區(qū)域.戴宗輝(2016)利用地震背景噪聲互相關(guān)方法得到了魯?shù)榈卣鹎昂蟮南鄬?duì)速度變化，發(fā)現(xiàn)較大幅度的同震波速變化主要集中于則木河斷裂—小江斷裂及蓮峰斷裂區(qū)域.龐衛(wèi)東等(2017)利用背景噪聲自相關(guān)方法得到了魯?shù)榈卣鹎昂笳鹪磪^(qū)周邊5個(gè)臺(tái)站的相對(duì)速度變化，發(fā)現(xiàn)魯?shù)榈卣鹁哂忻黠@的同震速度變化，推測(cè)可能是青藏高原東南緣與華南地塊相互作用而造成了研究區(qū)地殼構(gòu)造應(yīng)力變化.雖然前人已經(jīng)做了很多與地殼速度結(jié)構(gòu)變化有關(guān)的研究工作，但是均缺乏高精度空間分辨率的結(jié)果，所以提供魯?shù)榈卣鸢l(fā)生前后高精度的地殼波速變化的空間分布變化，有助于我們更加深入的了解地震孕育、發(fā)生和恢復(fù)的過程.

技術(shù)的進(jìn)步，尤其是地震臺(tái)網(wǎng)空間密度的增加，使監(jiān)測(cè)波速變化的靈敏度增加了一個(gè)數(shù)量級(jí)以上.監(jiān)測(cè)波速變化有使用固定的人工震源和被動(dòng)震源的方法.但是使用人工震源時(shí)，對(duì)所有危險(xiǎn)區(qū)域進(jìn)行監(jiān)測(cè)的費(fèi)用太高且在某些地區(qū)禁止這樣做，因此大多研究采用被動(dòng)震源的方法及數(shù)據(jù).而使用被動(dòng)震源數(shù)據(jù)的一種方法是使用剪切波分裂方法，這是通過各向異性介質(zhì)傳播的結(jié)果(Crampin et al.，1990；高原和吳晶，2008；常利軍等，2010).一種方法是比較具有相似位置和機(jī)制的重復(fù)地震的尾波中的散射波(Poupinet et al.，1984；Schaff and Beroza，2004；Taira et al.，2008；Zhao and Peng，2009).還有近10年發(fā)展迅速的基于背景噪聲的方法，該方法依賴于臺(tái)站的分布，已經(jīng)成功應(yīng)用于研究不同區(qū)域的波速變化(Brenguier et al.，2008；Obermann et al.，2014；劉志坤和黃金莉，2010).上述方法均不能得到速度變化的高精度空間分布，而能得到空間分布的一種方法是將不同時(shí)間段的速度模型相減(例如：Chiarabba et al.，2009；Patanè et al.，2006；Westman et al.，2012)，然而Julian和Foulger(2010)指出這種把不同時(shí)間段層析成像反演的速度模型相減的結(jié)果代表真實(shí)的速度變化的假設(shè)是值得懷疑的，因?yàn)橛捎诘卣鹞恢玫淖匀蛔兓鸬牡卣鹕渚€分布變化，即使速度結(jié)構(gòu)沒有變化，層析成像結(jié)果也會(huì)不同，即使地震位置沒有變化和地震臺(tái)站分布是固定的，但由于隨機(jī)觀測(cè)誤差，在反演出的模型中也會(huì)出現(xiàn)誤差.更好的方法是同時(shí)反演多個(gè)數(shù)據(jù)集，所以Julian和Foulger(2010)提出了一種依賴時(shí)間的成像方法，該方法通過最小化模型的差異和不同時(shí)間段的到時(shí)殘差來同時(shí)反演多個(gè)數(shù)據(jù)集，該方法對(duì)正演數(shù)據(jù)進(jìn)行了測(cè)試，與傳統(tǒng)的直接相減的方法相比，Julian和Foulger(2010)的方法能夠得到更可靠的速度隨時(shí)間的變化，但在一定程度上可能還會(huì)受到不同時(shí)間段數(shù)據(jù)分布不均勻的影響.Qian等(2018)提出了一種新的時(shí)移層析成像方法，該方法基于雙差層析成像方法(Zhang and Thurber，2003)，通過使用不同時(shí)間段的事件對(duì)構(gòu)建的到時(shí)差，可以克服部分的觀測(cè)誤差的影響和不同時(shí)間段的地震射線分布差異的影響.在反演過程中通過使用復(fù)雜的迭代策略從而解決了射線分布差異所產(chǎn)生的假性變化，為了消除不需要的模型差異，一個(gè)時(shí)間段的速度模型被用作另一個(gè)時(shí)間段的初始模型，因此這種新的時(shí)移層析成像可以確定可靠的波速隨時(shí)間的變化.

在本研究中，我們利用基于雙差層析成像的時(shí)移層析成像方法，使用2008年1月至2017年12月的地震數(shù)據(jù)進(jìn)行時(shí)移層析成像反演來得到2014年云南魯?shù)镸S6.5地震發(fā)生前后震源區(qū)及其周邊地區(qū)的地下介質(zhì)P波速度變化，并詳細(xì)分析了P波速度隨時(shí)間和空間的演化規(guī)律，進(jìn)而探討引起這種變化的原因.

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

1.1 基于雙差層析成像的時(shí)移層析成像方法

雙差層析成像方法通過利用絕對(duì)到時(shí)和相對(duì)到時(shí)同時(shí)得到地震速度模型和地震事件定位，與常規(guī)的地震層析成像方法相比，雙差層析成像方法能得到更加精確的源區(qū)域的速度結(jié)構(gòu)(Zhang and Thurber，2003，2006)，雙差層析成像方法能表示為

(1)

(2)

在雙差層析成像方法中，用來構(gòu)建事件對(duì)的事件均來自同一時(shí)間段，即(2)式中的事件i和事件j來自同一時(shí)間段，為了求解兩個(gè)時(shí)間段之間的速度變化，則要求事件i和事件j屬于兩個(gè)不同的時(shí)間段，即事件i發(fā)生在時(shí)間段1內(nèi)，事件j發(fā)生在時(shí)間段2內(nèi)，我們將(2)式改變?yōu)?/p>

(3)

(2)式和(3)式的不同點(diǎn)在于慢度參數(shù)δu1和δu2分別是時(shí)間段1和時(shí)間段2的慢度參數(shù)，并且事件i和事件j分屬于不同的時(shí)間段，這里為屬于不同時(shí)間段的事件構(gòu)建的到時(shí)差.

對(duì)于基于雙差層析成像方法的時(shí)移層析成像方法，首先利用雙差層析成像方法使用時(shí)間段1的數(shù)據(jù)來得到3D的速度模型，然后將所得到的3D速度模型作為初始模型使用(3)式反演時(shí)間段2相對(duì)于時(shí)間段1的速度變化δu2.在這個(gè)過程中，僅使用不同時(shí)間段的事件對(duì)構(gòu)建到時(shí)差，時(shí)間段1的慢度u1保持固定，到時(shí)差用來反演時(shí)間段2相對(duì)于時(shí)間段1的慢度變化δu2，同時(shí)也反演兩個(gè)時(shí)間段的事件位置.在計(jì)算程序中采用偽彎曲射線追蹤方法(Um and Thurber，1987)發(fā)現(xiàn)射線并計(jì)算走時(shí)，選取兩個(gè)時(shí)間段的地震匹配相對(duì)到時(shí)，對(duì)來自時(shí)間段1中的事件i進(jìn)行射線追蹤時(shí)用的速度模型是初始速度模型，來自時(shí)間段2中的事件j進(jìn)行射線追蹤時(shí)用的是每次迭代反演更新后的速度模型.模型表示為固定的三維節(jié)點(diǎn)，速度值用線性插值方法來進(jìn)行插值(Thurber，1983).

新的時(shí)移層析成像方法不僅如傳統(tǒng)層析成像方法一樣能提供高空間分辨率，同時(shí)能夠改進(jìn)傳統(tǒng)層析成像方法的一些缺陷，即能克服部分的觀察誤差的影響和不同時(shí)期地震射線分布的差異的影響，得出較為可靠的不同時(shí)間段之間的速度變化.因此這種新方法可以確定可靠的速度隨時(shí)間變化的空間分布.

1.2 數(shù)據(jù)

為了更好的監(jiān)測(cè)余震的發(fā)生，保證余震序列數(shù)據(jù)的完整性，云南省地震局在魯?shù)榈卣鹫鹪磪^(qū)附近架設(shè)了3個(gè)流動(dòng)臺(tái)，并加入了昭通市小孔徑測(cè)震臺(tái)網(wǎng)的2個(gè)臺(tái)站以及中國地震局地球物理研究所架設(shè)的巧家地震臺(tái)陣的5個(gè)臺(tái)站，加上原有的云南省區(qū)域地震臺(tái)網(wǎng)以及四川和貴州區(qū)域臺(tái)網(wǎng)的共享臺(tái)站，使得魯?shù)榈卣鹫鹪磪^(qū)及其附近的地震監(jiān)控能力大大提升.在魯?shù)榈卣鸢l(fā)生前，該區(qū)域的地震主要由云南區(qū)域臺(tái)網(wǎng)和四川、貴州區(qū)域臺(tái)網(wǎng)的共享臺(tái)站所記錄，臺(tái)站的分布會(huì)隨著臺(tái)網(wǎng)規(guī)劃和大地震的發(fā)生而產(chǎn)生變化.

本研究所使用的數(shù)據(jù)是由云南區(qū)域數(shù)字地震臺(tái)網(wǎng)所記錄到的地震事件的觀測(cè)報(bào)告數(shù)據(jù)，由于魯?shù)榈卣鸬挠嗾鸱植技性谝粋€(gè)小區(qū)域內(nèi)，射線分布少且不均勻，為此我們?cè)黾恿擞嗾饏^(qū)以外區(qū)域的地震來提高地震的射線覆蓋度.選取的研究區(qū)域?yàn)轸數(shù)榈卣鹩嗾饏^(qū)及其附近(26.6—27.6°N，103—104°E)，時(shí)間段為2008年1月1日至2017年12月31日，我們選取了由42個(gè)臺(tái)站(圖1)所記錄到的5623個(gè)地震事件(圖2)，每個(gè)事件由震中距小于200 km的5個(gè)以上臺(tái)站所記錄.一般來說，5個(gè)以上臺(tái)站所記錄到的地震事件的震相都較為清晰，拾取精度小于0.1 s.

為了確定魯?shù)榈卣鹫鹪磪^(qū)及其附近的震后速度變化，我們根據(jù)數(shù)據(jù)集劃分了五個(gè)地震數(shù)目大致相同的時(shí)間段(表1)，每個(gè)時(shí)段的P波射線數(shù)目大致相同.其中P1是魯?shù)榈卣鸢l(fā)生前，P2是同震期，P3是震后1—4個(gè)月，P4和P5分別是震后約1年和約3年的兩個(gè)時(shí)間段，這兩個(gè)時(shí)間段的地震數(shù)目大致相同.為了檢查觀測(cè)報(bào)告中震相數(shù)據(jù)的可靠性，繪制了五個(gè)時(shí)間段的P波震相走時(shí)曲線，并且使用所選初始一維速度模型(表2)繪制了理論走時(shí)曲線，如圖3所示.從圖中可看出P波走時(shí)在所選取的距離范圍內(nèi)呈現(xiàn)出明確的線性趨勢(shì)，并且震相的離散度較小，表明觀測(cè)震相數(shù)據(jù)的可靠性較高.觀測(cè)震相走時(shí)與理論走時(shí)曲線吻合較好，說明初始一維速度模型精度較高.圖4為五個(gè)時(shí)段內(nèi)的P波二維平面射線分布圖，可看出五個(gè)時(shí)段內(nèi)的魯?shù)榈卣鹫鹪磪^(qū)的射線密度分布均較密集.

表1 五個(gè)時(shí)段的數(shù)據(jù)Table 1 Data of the five periods

圖1 2014魯?shù)镸S6.5地震震源區(qū)附近的臺(tái)站分布黑框表示研究區(qū).Fig.1 The distribution of stations near Ludian 2014 MS6.5 earthquake source regionBlack frame denote research area.

圖2 研究區(qū)域活動(dòng)斷裂、地震分布及網(wǎng)格劃分紅點(diǎn)表示地震，“叉”字表示網(wǎng)格節(jié)點(diǎn)，黑線表示斷裂.Fig.2 Geological tectonic settings, seismic stations used in this study and grid nodes division in the studied areaRed dots denote earthquakes, crisscrosses denote grid nodes, black lines denote faults.

2 反演步驟

2.1 雙差地震層析成像方法

由時(shí)移層析成像方法的技術(shù)路線，在劃分完時(shí)間段后，首先使用雙差層析成像方法對(duì)每個(gè)時(shí)段的數(shù)據(jù)進(jìn)行反演，以得到每個(gè)時(shí)段的三維P波速度結(jié)構(gòu).首先進(jìn)行地震對(duì)匹配，選擇地震對(duì)之間的最大距離為20 km，每個(gè)地震最多可與10個(gè)地震組成地震對(duì)，最終構(gòu)建了相對(duì)到時(shí)數(shù)據(jù)，如表1所示.

在離散化的模型空間中進(jìn)行反演，選取(27.1°N，103.4°E)為坐標(biāo)原點(diǎn)，每個(gè)時(shí)段均劃分了相同的網(wǎng)格.綜合考慮每個(gè)時(shí)段內(nèi)的地震分布劃分網(wǎng)格，在反演之前進(jìn)行了大量的分辨率測(cè)試，以尋求最佳的網(wǎng)格分布.分別測(cè)試了3 km×3 km、5 km×5 km、10 km×10 km的網(wǎng)格間隔.由于地震分布不均勻，在地震分布密集的魯?shù)榈卣鹫鹪磪^(qū)的橫向分辨率達(dá)3 km×3 km，震源區(qū)外圍區(qū)域地震較少，橫向分辨率可達(dá)5 km×5 km，研究區(qū)邊緣地震分布很少，分辨率為10 km×10 km.如圖2所示，垂直向網(wǎng)格位于0、2、5、7、10、16、24、30、46 km深度.準(zhǔn)確的參考地殼模型是獲得穩(wěn)定結(jié)果的關(guān)鍵，所以在本研究中結(jié)合多個(gè)地殼模型(陳飛，2017；王光明等，2018；房立華等，2014；Riaz et al.，2017)建立了初始一維速度模型，其中陳飛(2017)是通過使用重力數(shù)據(jù)與面波數(shù)據(jù)聯(lián)合反演得到的全國三維速度結(jié)構(gòu)，王光明等(2018)和房立華等(2014)均是對(duì)魯?shù)榈卣鸺捌溆嗾疬M(jìn)行了重定位，Riaz等(2017)則是對(duì)魯?shù)榈卣鹫鹪磪^(qū)進(jìn)行了雙差層析成像.在最后得到的模型中，Moho面的深度范圍在30～46 km，這與鄭勇等(2013)通過背景噪聲成像得到的結(jié)果相似，每個(gè)時(shí)間段均采用相同的初始一維速度模型，如表2所示.

圖3 五個(gè)時(shí)段的P波震相走時(shí)Fig.3 P phase travel-time curve of the five periods

圖4 五個(gè)時(shí)段P波二維射線分布紅色點(diǎn)表示地震，黑色三角形表示臺(tái)站.Fig.4 Distribution of 2-D P wave ray paths for the five periodsRed dots denote earthquakes, black triangles denote stations.

表2 一維P波速度模型表Table 2 1-D P wave velocity model

在反演中，為了反演出較小尺度的異常并能找到地質(zhì)上可行特征的反演模型而系統(tǒng)地測(cè)試了不同組合的平滑參數(shù)(η)和阻尼參數(shù)(ε)，最后設(shè)置了一個(gè)較小的平滑參數(shù)(η=6)，然后測(cè)試了不同的ε，并選擇一個(gè)可以確保條件數(shù)在150至200之間的ε，保持條件數(shù)的數(shù)值范圍較低但不太窄是很重要的，因?yàn)檎鎸?shí)的條件數(shù)與使用LSQR所計(jì)算的條件數(shù)不同，可能更小(Zhang and Thurber，2007).所以最后5個(gè)時(shí)間段的ε選取會(huì)有一定差別.對(duì)5個(gè)時(shí)段的數(shù)據(jù)均進(jìn)行了16次迭代以確保最終的結(jié)果能夠達(dá)到令人滿意的收斂，前6次迭代，重點(diǎn)放在擬合絕對(duì)數(shù)據(jù)上，包括3個(gè)速度與地震位置的聯(lián)合反演及3個(gè)重定位的反演.而在后10次迭代中，重點(diǎn)放在到時(shí)差數(shù)據(jù)上，包括7個(gè)速度與地震位置的聯(lián)合反演及3個(gè)重定位的反演.反演后5個(gè)時(shí)段到時(shí)殘差的均方根都有顯著的降低，如表3所示，表明了所獲得的最終模型的穩(wěn)定性.

圖5 五個(gè)時(shí)段不同深度處的P波速度棋盤測(cè)試結(jié)果(a) 時(shí)段1； (b) 時(shí)段2； (c) 時(shí)段3； (d) 時(shí)段4； (e) 時(shí)段5.Fig.5 Checkerboard resolution test of P wave velocity for the five periods(a) P1; (b) P2; (c) P3; (d) P4; (e) P5.

表3 五個(gè)時(shí)段反演前后觀測(cè)到時(shí)差與理論到時(shí)差之間的均方根殘差變化Table 3 The RMS residuals between observed and predicted differential travel times based on 1D model and 3D model for the five periods

2.2 時(shí)移層析成像方法

同時(shí)在本文中對(duì)5個(gè)時(shí)段的反演結(jié)果的分辨率進(jìn)行棋盤測(cè)試，在初始一維速度模型上分別增加了±5%的正負(fù)相間速度擾動(dòng).如圖5所示，在此選取了5、10、16 km深度剖面的分辨率測(cè)試結(jié)果，可以看出P2～P5這4個(gè)時(shí)段內(nèi)的三個(gè)深度的魯?shù)榈卣鹫鹪磪^(qū)的P波速度均恢復(fù)較好，而P1時(shí)間段由于射線數(shù)目較少，在魯?shù)榈卣鹫鹪磪^(qū)恢復(fù)的相對(duì)較差.為了進(jìn)一步量化模型分辨率，我們使用了衍生加權(quán)總和(DWS，Thurber and Eberhart-Philips，1999)對(duì)模型分辨率進(jìn)行補(bǔ)充.DWS值代表給定網(wǎng)格節(jié)點(diǎn)附近所有射線長(zhǎng)度(乘以數(shù)據(jù)權(quán)重)的總和，一般來說，DWS值是解釋可接受分辨率區(qū)域的有效工具.總的來說，DWS=100(圖6中的白色等值線)包圍的基本是棋盤測(cè)試恢復(fù)良好的區(qū)域(圖5).圖6為最后反演出的5個(gè)時(shí)間段的5、10、16 km深度的P波速度分布圖，可看出在5個(gè)時(shí)段的不同深度處均能較好的表現(xiàn)出P波速度.

在使用雙差層析成像方法計(jì)算出5個(gè)時(shí)段的P波速度模型后，接著用時(shí)移層析成像方法得到相鄰時(shí)間段的P波速度變化.在這里首先構(gòu)建兩個(gè)時(shí)段的事件對(duì)的到時(shí)差，為了盡量多的構(gòu)建分屬于兩個(gè)時(shí)段的兩個(gè)事件的事件對(duì)，要求地震對(duì)之間的最大距離為50 km，一個(gè)時(shí)間段內(nèi)的每個(gè)事件最多可以與30個(gè)事件相連接.從P1到P2，構(gòu)建了平均空間距離為7 km的616007個(gè)到時(shí)差，挑選出了70470個(gè)分屬于兩個(gè)時(shí)段的事件對(duì)的到時(shí)差；從P2到P3，構(gòu)建了平均空間距離為4.2 km的667502個(gè)到時(shí)差，最后挑選出250680個(gè)符合要求的到時(shí)差；從P3到P4，我們構(gòu)建了平均空間距離為4.8 km的462766個(gè)到時(shí)差，最后挑選出符合要求的到時(shí)差164765個(gè)；從P4到P5，我們構(gòu)建了平均空間距離為7.3 km的665468個(gè)到時(shí)差，最后挑選出了206155個(gè)符合要求的到時(shí)差.

對(duì)于時(shí)移層析成像方法，使用由雙差層析成像方法所得到的第一個(gè)時(shí)間段的速度模型作為初始模型，重定位結(jié)果作為初始位置，然后使用兩個(gè)相鄰時(shí)間段所構(gòu)建的到時(shí)差來得到速度隨時(shí)間的變化.為了確保最后的結(jié)果能很好的收斂，反演過程經(jīng)過了32次迭代，反演后到時(shí)殘差的均方根變化如表4所示，4個(gè)到時(shí)殘差的均方根均有顯著的下降，說明最后得到的速度變化與兩個(gè)相鄰時(shí)間段內(nèi)的事件非常合適.

表4 不同時(shí)段利用時(shí)移層析成像方法反演前后觀測(cè)到時(shí)差與理論到時(shí)差之間的均方根殘差變化Table 4 The RMS residuals between observed and predicted differential travel times based on initial 3D model and final 3D model for the different periods

同時(shí)還進(jìn)行了棋盤測(cè)試，與3D層析成像的棋盤測(cè)試相類似，在第一個(gè)時(shí)間段的初始3D模型中加入了±5%的正負(fù)相間速度擾動(dòng)，圖7顯示了4個(gè)相鄰時(shí)段棋盤測(cè)試的恢復(fù)結(jié)果.P1與P2，地震分布和所使用臺(tái)站差異較大，雖然盡量選取多的地震對(duì)，但是仍然過于分散，所以分辨率恢復(fù)的并不好，這在一定程度上會(huì)影響反演結(jié)果的精度.其余時(shí)段由于地震分布和所使用臺(tái)站差異不大，所以分辨率相對(duì)較好.但是棋盤測(cè)試恢復(fù)出來的模型分辨率與反演實(shí)際數(shù)據(jù)獲得的速度模型之間沒有完全的相互映射，這是因?yàn)檫@兩個(gè)模型具有不同的速度(Share et al.，2019)，所以同時(shí)也使用了DWS值(圖7中的白色虛線)來解釋可接受分辨率的區(qū)域，如圖7所示.由于在魯?shù)榈卣鹫鹪磪^(qū)有相對(duì)多數(shù)量的射線覆蓋，該區(qū)域棋盤恢復(fù)的相對(duì)較好，圖中也顯示了DWS值大于300的區(qū)域基本覆蓋了震源區(qū)，所以接下來主要討論魯?shù)榈卣鹫鹪磪^(qū)的P波速度變化.

圖6 五個(gè)時(shí)段不同深度處的P波速度和地震分布黑點(diǎn)表示地震，白色五角星表示主震，白色實(shí)線表示DWS值大于100的區(qū)域. (a) 時(shí)段1； (b) 時(shí)段2； (c) 時(shí)段3； (d) 時(shí)段4； (e) 時(shí)段5.Fig.6 Distribution of P-wave velocity and earthquakes at different depths for the five periodsBlack dots denote earthquakes, white five-pointed stars denote the main shock, white lines denote the areas that their DWS value are greater than 100. (a) P1; (b) P2; (c) P3; (d) P4; (e) P5.

圖7 不同時(shí)段之間利用時(shí)移層析成像方法得到的不同深度處的棋盤測(cè)試結(jié)果白色虛線表示DWS值大于300的區(qū)域. (a) 時(shí)段2相對(duì)于時(shí)段1； (b) 時(shí)段3相對(duì)于時(shí)段2； (c) 時(shí)段4相對(duì)于時(shí)段3； (d) 時(shí)段5相對(duì)于時(shí)段4.Fig.7 Recovered checkerboard models for different periods by using time-lapse tomographyWhite dotted lines denote the area that their DWS values are greater than 300. (a) P2-P1; (b) P3-P2; (c) P4-P3; (d) P5-P4.

3 2014年云南魯?shù)榈卣鸬乃俣茸兓?/h2>

根據(jù)分辨率和DWS值的分布情況主要討論魯?shù)榈卣鹫鹪磪^(qū)P波速度變化，所以在圖8中主要顯示了分辨率較好的魯?shù)榈卣鹫鹪磪^(qū)在5、10、16 km深度處的不同時(shí)間段之間的P波速度變化.由圖8可以清楚的看出魯?shù)榈卣鸢l(fā)生前后的P波速度的降低及恢復(fù).其中圖8a是由P1和P2兩個(gè)時(shí)段的數(shù)據(jù)得到的魯?shù)榈卣鹜鹚俣茸兓?，可以看出每個(gè)深度處的震源區(qū)的速度均下降，下降的幅度不大，約在0.2%以下，大多為百分之零點(diǎn)零幾，其中16 km深度處的速度下降幅度小于5 km和10 km深度處.圖8b是P2和P3兩個(gè)時(shí)段的數(shù)據(jù)得到的震后1—4個(gè)月相對(duì)同震期的速度變化，可看出震源區(qū)的速度大幅下降，下降幅度大多在0.2%～3%內(nèi)，這個(gè)范圍與前人對(duì)相同量級(jí)地震發(fā)生前后速度變化研究結(jié)果相一致.例如Li等(1998)計(jì)算的2004年P(guān)arkfieldM6地震發(fā)生后速度下降了2.5%；Nishimura等(1998)計(jì)算的1998年Intra-PlateM6.1地震發(fā)生后速度下降了0.3%～10%.相比P2-P1，P3-P2的速度下降的幅度有很大的變化，達(dá)到了最大的速度降幅，說明魯?shù)榈卣鸢l(fā)生后速度并沒有立即下降到最大，而是在震后1—4個(gè)月期間速度下降到了最大.這與龐衛(wèi)東等(2017)利用背景噪聲自相關(guān)方法得到的魯?shù)榈卣鸢l(fā)生后最大波速降出現(xiàn)在震后1個(gè)月的時(shí)間相一致，也與劉志坤和黃金莉(2010)計(jì)算的汶川地震、王俊等(2020)計(jì)算的蘆山地震發(fā)生后的最大波速降出現(xiàn)在震后1—4個(gè)月內(nèi)相符合.

圖8 不同深度處的不同時(shí)段之間的P波速度變化和地震分布(a) 時(shí)段2相對(duì)于時(shí)段1； (b) 時(shí)段3相對(duì)于時(shí)段2； (c) 時(shí)段4相對(duì)于時(shí)段3； (d) 時(shí)段5相對(duì)于時(shí)段4.Fig.8 Distribution of temporal P-wave velocity changes and earthquakes at different depths for different periods(a) P2-P1; (b) P3-P2; (c) P4-P3; (d) P5-P4.

圖9 不同時(shí)段之間的P波速度變化及地震事件沿垂直剖面AA′和BB′的分布圖(剖面位置見圖2)紅色虛線表示DWS大于300的區(qū)域. (a) 時(shí)段2相對(duì)于時(shí)段1； (b) 時(shí)段3相對(duì)于時(shí)段2； (c) 時(shí)段4相對(duì)于時(shí)段3； (d) 時(shí)段5相對(duì)于時(shí)段4.Fig.9 Distribution of temporal P-wave velocity changes along the cross sections AA′, BB′ for different periods(section location is in figure 2)Red dotted lines denote the areas that their DWS values are greater than 300. (a) P2-P1; (b) P3-P2; (c) P4-P3; (d) P5-P4.

在震后速度下降達(dá)到最大值后，震源區(qū)開始愈合，速度開始上升，由P4-P3和P5-P4的速度變化可以清楚的看出震源區(qū)在時(shí)空上的愈合過程.愈合從淺層開始，首先是5 km深度處的震源區(qū)的速度開始上升，上升幅度大約為0.1%以下，而10 km和16 km深度處的速度仍然是在下降，說明淺層區(qū)域開始愈合.隨著時(shí)間發(fā)展，愈合延伸至更深的深度，5 km和10 km深度處的震源區(qū)的速度也上升，16 km深度處的震源區(qū)的速度雖然在下降，但是有一部分區(qū)域的速度在上升，說明愈合在16 km深度處也在進(jìn)行，上升的幅度約在0.2%～2%，說明基本在向震前的水平進(jìn)行愈合.

總的來看，2014年魯?shù)榈卣鸢l(fā)生后，震源區(qū)的速度變化特征為，同震期速度下降，但沒有下降至最大，而是在震后1—4個(gè)月內(nèi)速度下降至最大，然后愈合過程開始，愈合從淺層逐漸發(fā)展至深層.值得注意的一點(diǎn)是，由于我們的觀測(cè)同時(shí)對(duì)應(yīng)于空間和時(shí)間的平均值，因此不應(yīng)過度解釋所有變化幅度的數(shù)值大小.

由于魯?shù)榈卣鹩嗾鸱植汲使曹椥头植?，所以我們還做了如圖2所示的兩個(gè)剖面的速度變化示意圖(圖9)，其中AA′剖面主要沿著地震序列的SES方向穿過震源區(qū)，BB′剖面主要沿著地震序列的近EW向穿過震源區(qū)，這樣我們就能看到兩個(gè)余震優(yōu)勢(shì)分布方向上的速度變化的時(shí)空分布.在同震期，AA′和BB′兩個(gè)剖面下方均為速度下降，且靠近AA′剖面的A′端即SES端的區(qū)域和靠近BB′剖面的B′端即近E向端的區(qū)域的速度下降幅度略大，這可能是因?yàn)檫@兩個(gè)區(qū)域是魯?shù)榈卣鸢l(fā)生后余震最初集中發(fā)生的區(qū)域.在震后1—4月，兩個(gè)剖面下方速度下降劇烈，其中靠近A′端的區(qū)域和靠近B端的區(qū)域的速度下降最為劇烈，這是因?yàn)檎鸷笥嗾鸱植及l(fā)生了變化，不再局限于主震附近，近EW向分支的余震繼續(xù)往W向發(fā)展，SES向分支的余震則仍然往SES向發(fā)展(王光明等，2018)，對(duì)于余震發(fā)展的方向的解釋，Riaz等(2017)和趙小艷和孫楠(2014)均認(rèn)為是由于震源區(qū)的北部存在一個(gè)高速異常，對(duì)余震有阻礙作用.兩個(gè)剖面的震后愈合恢復(fù)過程均是從淺層開始，AA′剖面余震分布較為集中的SES端的區(qū)域及BB′剖面余震最開始發(fā)生的靠近近E端的區(qū)域的下方至5 km深度的速度上升，隨著時(shí)間發(fā)展，AA′剖面的SES端區(qū)域下方速度上升幅度增大，延伸至約8 km深度，同時(shí)BB′剖面下方速度上升的區(qū)域延伸至近W端，即后來余震發(fā)展的方向，且該區(qū)域下方速度上升幅度較大，同時(shí)速度上升的區(qū)域擴(kuò)大至約15 km深度.

兩個(gè)剖面的速度變化的發(fā)展過程同時(shí)也反映了余震分布的發(fā)展過程.結(jié)合房立華等(2014)和王光明等(2018)對(duì)魯?shù)榈卣鹩嗾鹦蛄械闹囟ㄎ谎芯拷Y(jié)果，可以得出：魯?shù)榈卣鸢l(fā)生后余震主要分布在主震附近及其東南側(cè)，這些區(qū)域下方的速度開始下降，下降幅度較小，此后主震西側(cè)開始出現(xiàn)余震，而東南側(cè)的余震分布范圍基本穩(wěn)定，余震呈現(xiàn)出不對(duì)稱的共軛分布特征，主震西側(cè)及東南側(cè)這兩個(gè)區(qū)域下方的速度下降劇烈，速度下降至最大后愈合過程開始，余震最初出現(xiàn)的主震附近及其東南側(cè)區(qū)域下方的淺層的速度上升，隨著時(shí)間的發(fā)展，愈合延伸至余震后來發(fā)展的區(qū)域，即主震西側(cè)，且該區(qū)域下方速度上升幅度增大，整個(gè)近EW向分支區(qū)域下方速度上升的深度擴(kuò)展至約15 km深度處，而余震分布穩(wěn)定的主震東南側(cè)區(qū)域下方速度上升幅度增大，范圍擴(kuò)大至約8 km.

4 討論

一般認(rèn)為，引起地震發(fā)生前后波速變化的物理機(jī)制有以下幾種：地震斷層區(qū)的破壞與愈合(Li et al.，1998；Vidale and Li，2003)，強(qiáng)地面運(yùn)動(dòng)引起的地表破壞(Peng and Ben-Zion，2006)，斷層區(qū)及周邊地殼介質(zhì)的應(yīng)力變化(Wegler et al.，2009；Brenguier et al.，2008；Niu et al.，2008；Pei et al.，2019)，地下水位變化(Sens-Sch?nfelder and Wegler，2006)等(劉志坤和黃金莉，2010).而要準(zhǔn)確地判定引起地震波速變化的物理機(jī)制是非常困難的.魯?shù)榈卣鸢l(fā)生后，對(duì)魯?shù)榈卣鹫鹪催\(yùn)動(dòng)學(xué)、余震分布特征、應(yīng)力變化以及地質(zhì)構(gòu)造等方面的研究有很多，可以將本研究所得到的結(jié)果與已有的地質(zhì)與地球物理研究成果相結(jié)合，從而進(jìn)一步分析造成速度變化的可能機(jī)制.

一般來說，速度變化能反映強(qiáng)地面運(yùn)動(dòng)所造成的巖石損傷，預(yù)計(jì)主要發(fā)生在頂部幾百米(Pei et al.，2019)，而本研究在中上地殼均發(fā)現(xiàn)了P波速度的明顯變化，因此大地震的強(qiáng)地面運(yùn)動(dòng)引起的淺部地殼介質(zhì)破壞對(duì)本研究得到的波速變化結(jié)果影響不大，不是主要因素.

魯?shù)榈卣鸢l(fā)生后，在由中國地震局組織的野外應(yīng)急科學(xué)考察中并沒有發(fā)現(xiàn)明顯的地震地表破裂現(xiàn)象(徐錫偉等，2014)，但是李西等(2018)基于詳細(xì)的野外地質(zhì)地貌調(diào)查發(fā)現(xiàn)了自龍頭山鎮(zhèn)以南謝家營(yíng)盤—光明村—王家坡一段長(zhǎng)約8 km斷續(xù)展布的地表破裂帶，呈左行右階排列，伴有30～35 cm走滑分量，局部段落伴有15～25 cm傾滑分量.這段地表破裂帶主要位于魯?shù)榈卣鹩嗾饚ES向分支上，而主震附近及近EW向分支上并沒有觀察到地表破裂，而在本研究中整個(gè)震源區(qū)均有明顯的波速變化，表明本研究得到的波速變化對(duì)整個(gè)震源區(qū)的介質(zhì)性質(zhì)變化都有所反映.并且地表破裂帶應(yīng)該僅對(duì)地表淺層區(qū)域有較大影響，對(duì)深層區(qū)域的介質(zhì)影響不大，而本研究得到的波速變化從淺層至深層均有明顯變化，所以判斷地震斷層區(qū)的破壞并不是影響本研究所得到的波速變化的主要原因，但是對(duì)淺層區(qū)域及主震東南方向區(qū)域介質(zhì)變化具有一定的影響.

破裂區(qū)深度內(nèi)的結(jié)構(gòu)變化和開放裂縫導(dǎo)致應(yīng)力場(chǎng)的變化可能發(fā)生在不同的深度(Pei et al.，2019)，通過觀察圖8發(fā)現(xiàn)同震期的5 km和10 km深度處及震后1—4個(gè)月的10 km和16 km深度處均呈現(xiàn)出按象限分布的現(xiàn)象，且在以前的很多研究中(Wegler et al.，2009；Cheng et al.，2010；劉志坤和黃金莉，2010；王俊等，2020)均發(fā)現(xiàn)了這一現(xiàn)象.為了探究波速變化呈象限分布是否是由地震的應(yīng)力變化引起的，查找了相關(guān)的研究結(jié)果.繆淼和朱守彪(2016)計(jì)算了魯?shù)榈卣痨o態(tài)庫倫應(yīng)力變化，發(fā)現(xiàn)兩個(gè)方向的余震分布都主要集中在主震產(chǎn)生的庫倫應(yīng)力增大區(qū)，付芮等(2015)也發(fā)現(xiàn)余震主要位于沿破裂面及主破裂以西的庫侖應(yīng)力增強(qiáng)區(qū).對(duì)比圖8發(fā)現(xiàn)5 km和10 km深度處及震后1—4個(gè)月的10 km及16 km深度處P波速度下降的區(qū)域基本為余震主要分布區(qū)，說明P波速度下降區(qū)主要對(duì)應(yīng)于應(yīng)力增大區(qū)，表明魯?shù)榈卣鹬髡甬a(chǎn)生的應(yīng)力變化與P波速度變化之間有緊密的聯(lián)系，這也從另一角度反映出，地震引起的應(yīng)力變化是造成魯?shù)榈卣鹚俣茸兓囊粋€(gè)重要原因.并且在第3節(jié)內(nèi)容中提到魯?shù)榈卣鸢l(fā)生后最大速度下降發(fā)生在震后1—4個(gè)月，并沒有發(fā)生在同震期，而速度的下降及上升在空間上的分布變化均與余震發(fā)生的空間變化有關(guān)，這表明地震波速變化與余震的時(shí)空分布有關(guān)，可能是余震的動(dòng)態(tài)和靜態(tài)應(yīng)力造成震源區(qū)介質(zhì)物理屬性的改變從而影響了地震波速的變化.

魯?shù)榈卣鸢l(fā)生在南北地震帶南段，青藏高原東南緣鮮水河—小江斷裂系以東走向NE向昭通—蓮峰斷裂帶西段，北臨NNW至近NS向大涼山斷裂和馬邊—鹽津斷裂，屬青藏高原東南部川滇菱形塊體與華南地塊之間I級(jí)塊體邊界構(gòu)造帶(徐錫偉等，2014).龐衛(wèi)東等(2017)認(rèn)為魯?shù)榈卣鹗乔嗖馗咴瓥|南緣與華南地塊相互作用而造成地殼內(nèi)部構(gòu)造應(yīng)力變化的結(jié)果；程佳等(2016)認(rèn)為魯?shù)榈卣鹗艿搅酥饕獊碜杂趧t木河斷裂帶上的應(yīng)力作用并積累地震能量；洪敏利用應(yīng)變場(chǎng)時(shí)序分析法(洪敏等，2014)計(jì)算出的魯?shù)榈卣鸢l(fā)生前的云南地區(qū)應(yīng)變場(chǎng)分布發(fā)現(xiàn)滇東北地區(qū)出現(xiàn)明顯的擠壓活動(dòng)增強(qiáng)(1)云南省地震局形變測(cè)量中心 2018年年度會(huì)商報(bào)告..所以區(qū)域應(yīng)力場(chǎng)的變化對(duì)魯?shù)榈卣鸬陌l(fā)生有著很大的影響，而魯?shù)榈卣饘?duì)后續(xù)余震具有觸發(fā)作用，從而對(duì)震源區(qū)介質(zhì)的物理性質(zhì)變化有著一定的影響.

5 結(jié)論

運(yùn)用基于雙差層析成像的時(shí)移層析成像方法能夠得到高精度的地震波速變化的時(shí)空分布.在本研究中我們利用云南區(qū)域數(shù)字地震臺(tái)網(wǎng)所記錄到的地震資料，使用該方法來得到了2014年魯?shù)镸S6.5地震發(fā)生前后的P波速度變化的時(shí)空分布，結(jié)果表明：2014年魯?shù)榈卣鸢l(fā)生后，同震期的震源區(qū)速度下降，但沒有下降至最大，而是在震后1—4個(gè)月內(nèi)速度下降至最大，接著P波速度開始上升，開始愈合，愈合過程從淺層逐漸發(fā)展至深層，并逐漸恢復(fù)至震前水平.由于魯?shù)榈卣鹩嗾饏^(qū)呈共軛形態(tài)，所以我們做了分別穿過兩個(gè)余震分支的垂直剖面，發(fā)現(xiàn)P波速度變化與余震分布變化相一致，余震最初分布的主震附近及其東南側(cè)區(qū)域下方的速度最先下降，然后主震西側(cè)開始出現(xiàn)余震，則主震西側(cè)及東南側(cè)這兩個(gè)區(qū)域下方的速度下降劇烈，愈合過程則是余震最初出現(xiàn)的區(qū)域下方開始，之后延伸至余震后來發(fā)展的區(qū)域.

進(jìn)一步將所得到的結(jié)果和已有的地質(zhì)與地球物理觀測(cè)相結(jié)合來分析引起這種變化的可能機(jī)制，認(rèn)為強(qiáng)地面運(yùn)動(dòng)引起的地表破壞及地震斷層區(qū)的破壞并不是影響本研究所得到的波速變化的主要原因，但是對(duì)淺層區(qū)域介質(zhì)變化具有一定的影響.由于在空間上魯?shù)榈卣饝?yīng)力變化與波速變化之間有較緊密的聯(lián)系，且波速變化與余震分布變化相一致，所以我們認(rèn)為地震引起的應(yīng)力變化是造成魯?shù)榈卣鸩ㄋ僮兓囊粋€(gè)重要原因，余震的動(dòng)態(tài)和靜態(tài)應(yīng)力會(huì)造成震源區(qū)介質(zhì)物理屬性的改變從而影響了地震波速的變化.并且由以往的研究也可得出區(qū)域構(gòu)造應(yīng)力場(chǎng)變化對(duì)魯?shù)榈卣鸬陌l(fā)生有很大作用，從而對(duì)震源區(qū)地下物理性質(zhì)的改變?cè)斐捎绊?

致謝感謝中國科技大學(xué)地球和空間科學(xué)學(xué)院胡景博士在本研究過程中提供的幫助，感謝云南省地震局的王光明助理研究員和洪敏高工對(duì)文章的熱心幫助，感謝三位匿名審稿專家提出的寶貴意見和建議.