騰沖地區近震S波分裂研究

王雪鶴， 李永華,2*， 呼楠

1 中國地震局地球物理研究所， 北京 100081 2 中國地震局震源物理重點實驗室， 北京 100081 3 陜西省地震局， 西安 710068

0 引言

位于云南西南地區的騰沖地塊內發育第四紀活火山活動，是中國大陸最新的火山之一，其最新的一次火山噴發發生在300多年前.已有地質和地球物理研究表明，騰沖火山為緬甸微板塊向東俯沖誘導巖漿上涌所致(Lei et al., 2009；胥頤等，2012).地質調查還表明，在騰沖地區東西寬50 km，南北長90 km范圍內共有68座火山，組成了7個火山群(姜朝松，1998).針對騰沖火山群的監測顯示，該區火山區的巖漿處于一種活動狀態，并推測其作為一座活火山仍然有很大的噴發危險(姜朝松，1998).此外，該火山區地質構造復雜、地震活動頻繁，一直受到地震學家的廣泛關注.

通過地球物理方法確定火山區深部結構可為理解該火山的地質過程及演化提供重要約束.前人(謝富仁等，2001；Xu et al., 2010; Zhao et al., 2013； Tian et al., 2019)基于區域地震震源機制的應力張量反演研究了騰沖地區的構造應力場特征，結果顯示騰沖火山區東部的最大主壓應力軸方向為NE-SW向，而西部為N-S或NNW-SSE方向.但根據震源機制反演應力張量需要知道哪個節面是真正的斷層面，如果缺失這一信息，且斷層面和輔助節面有交換，應力反演就可導致不精確的結果.此外，在某種情況下，由于區域地震震源機制解的不均一性使得很難獲取可信的平均應力場特征(許忠淮,1985；Bora et al., 2018; Wan et al., 2016).已有研究表明，地震震源機制和地殼地震各向異性等地震學研究可為構造活躍區地殼變形過程的確定提供重要約束(Johnson et al., 2011; Araragi et al., 2015; Bora et al., 2018).

地震各向異性是地球上普遍存在的一種現象，而剪切波分裂現象是地震各向異性最直觀的表現形式(Crampin, 1994; Boness and Zoback, 2006).地殼中周期性薄互層、礦物晶體和裂隙的定向排列是誘發地震各向異性的主要因素(Crampin,1994).近震S波分裂主要反映了中上地殼介質的各向異性特征，一般認為中上地殼各向異性的成因主要分為兩大類，一類是構造成因(包括定向排列的礦物，巖石組構及古應力殘留的介質定向特征等等)，其產生的快波方向一般與對應的剪切組構的方向一致，慢波延遲時間代表剪切程度，這種類型中最常見的就是由大型走滑斷裂所引起的各向異性，而在火山區定向排列的巖墻也可能導致地震各向異性(Johnson et al., 2011).另外一類是應力成因，是由于介質內部受主壓應力引起的定向排列的張開的流體填充或不填充的微裂隙、閉合的可見微破裂引起，其產生的快波方向平行于區域水平最大主壓應力，慢波延遲時間代表微裂隙或破裂的密度(Boness and Zoback, 2006).

另一方面，已有研究表明，火山區往往可以觀察到顯著的地震各向異性時空變化，并將其歸結為淺層巖漿侵入引起的淺部地殼應力變化所致(Gerst and Savage, 2004; Unglert et al., 2011; Illsley-Kemp et al., 2018).由于巖漿系統中的壓力隨時間的增加/減少會導致地震各向異性快波偏振方向的改變，在一些火山區，剪切波分裂參數的監測已被建議作為預測火山爆發的一種可能方法(Gerst and Savage,2004; Unglert et al., 2011; Illsley-Kemp et al., 2018； Crampin and Gao, 2018).為此，區分騰沖火山區地震各向異性的成因，確定騰沖火山區區域應力與構造之間的相互關系，不僅可為更好地了解研究區的地殼演化提供地震學約束，還將為應用剪切波分裂這一手段監測火山活動提供堅實的基礎.

前人(毛慧玲和秦嘉政,2011; 孫長青等, 2013; Shi et al.,2012;Zhang et al.,2018)利用研究區固定臺站的波形記錄，通過接收函數和近震剪切波分裂方法對研究區地殼各向異性特征進行了分析，推測研究區中上地殼各向異性與主壓應力方向一致，但受騰沖地區固定臺站分布限制，無法反映出研究區應力橫向變化特征.

本文研究對2008-10—2019-07騰沖火山區15個固定和流動地震臺站記錄到的近震波形資料進行了剪切波分裂分析，以確定騰沖火山區的地殼應力狀態時空變化.地殼VP/VS值與地殼的物質組成、溫度狀態等密切相關，對其時空變化特征進行研究，可為巖漿活動的確定提供重要的約束(Unglert et al., 2011).為此，本文也對騰沖火山區上地殼各臺站的波速比值進行了計算.該區地殼地震各向異性和波速比分布特征的確定，不僅為理解該火山區的動力學演化過程提供了約束，也為將來應用地震各向異性手段監測火山活動奠定了基礎.

1 數據與方法

1.1 數據

本研究采用的數據包括四個部分：中國國家數字地震臺網(CNDSN)布設在騰沖火山區的火山臺站(8個)、國家臺站TNC、區域臺站MAS總共10個固定臺站從2008年10月—2019年7月的波形數據，以及2014年1月—2014年12月布設在騰沖火山區5個流動臺站記錄的波形數據.其中10個固定臺站配備了采樣率為100 Hz的BBVAS60和KS2000寬帶地震儀以及EDAS-24IP數據收集器系統.研究中收集了上述15個地震臺站記錄的發生在騰沖火山區及其周圍鄰區震級小于4級的地震，這些地震的震源深度在5～25 km之間不等，80%地震的震源深度小于15 km(圖1).

當近震S波入射到地表時，如果入射角大于臨界角(sin-1(VS/VP))時會發生全反射現象，為了避免這種全反射現象引起的波形畸變，我們需要選擇在剪切波橫波窗口內的地震波形數據進行分析.對于泊松比為0.25的泊松介質來說，其臨界角為37°，考慮到近地表低速層導致的地震波射線路徑彎曲，本文研究參照前人的做法(吳晶等，2010；Savage et al.，2016)，將入射角擴大至45°，即選取入射角小于45°的地震波形記錄，以確保所有用于分裂的波形記錄滿足剪切波窗口的要求.研究中我們主要依據中國地震臺網中心提供的地震目錄，選取了震中距<25 km，震級小于4級，震源深度大于震中距的587個近震波形數據用于剪切波分裂分析.研究中我們人工手動方法拾取了Pg，Sg波到時用于下一步的分析.

1.2 方法

本研究使用MFAST方法進行了剪切波分裂分析.MFAST是一種基于對剪切波分裂參數(快波偏振方向和延遲時間)的可能值進行網格搜索反演的全自動剪切波分裂程序(Savage et al., 2010)，已有多位學者利用該方法進行了各向異性研究(Unglert et al., 2011; Li et al., 2015; Syuhada et al., 2017; Maher and Kendall, 2018)，證明了該方法的可靠性.

MFAST方法將14個帶有特定拐角頻率的帶通濾波應用到每一個數據上以選擇最佳濾波.只有信噪比大于3的地震波形被用于下一步計算.對波形進行濾波后，MFAST在多個時間窗上執行SC91算法，在S波到時前后3 s內選擇噪聲和信號窗口，通過定義窗口的固定長度來避免個人主觀標準.傳統情況下，剪切波分裂窗口都是由我們手動選擇的，而手動選擇窗口存在一定的主觀性，因為選擇窗口的不同往往就意味著計算出的不同分裂結果.MFAST方法改善了這一問題，將選擇窗口的過程完全自動化，對多個窗口進行分裂分析，以得到在不同的窗口上都穩定的測量結果.

MFAST方法基于依賴于濾波頻段和分裂窗口選擇的SC91算法，以Teanby等(2004)的方法為基礎，采用了14個不同頻率的帶通濾波器，并基于信噪比、延遲時間和快波偏振方向誤差等一系列參數自動對分裂測量結果進行質量分級(A—D).在本研究中我們選擇網格搜索的最大延遲時間tlagmax為0.5 s.評定為A等級的結果需滿足延遲時間δt<0.8tlagmax,信噪比SNR>4,同時快波偏振方向誤差δφ<10°的條件；評定為B等級的結果需滿足延遲時間δt<0.8tlagmax, 信噪比SNR>3,同時快波偏振方向誤差δφ<25°的條件.對信噪比帶寬、濾波器、分裂窗口、聚類分析和分級原理的更具體解釋見Savage等(2010).

本文對15個臺站記錄的地震波形進行了分裂計算，最終得到了1761個分裂結果，我們依據以下標準篩選分裂結果：(1)觀察旋轉后的切向分量的能量是否達到最小(圖2b)；(2)觀察是否具有穩定的計算結果(圖2c,d)；(3)觀察校正后的粒子運動是否變為近直線形狀，校正后快慢波波形是否擬合(圖2e)；(4)如圖2f所示觀察等值線是否比較收斂(圖2f).圖2和圖 3闡述了不同質量等級的地震事件的分析過程.

圖1 研究區地質背景與地震臺站分布 紅色、黃色三角分別代表固定和流動臺站，研究區7個火山群分別用空心橢圓標示，灰點和紫點表示2008年以來所有ML≤4.0級的地震和最終用于橫波分裂分析的事件.右下方的柱狀圖所用地震事件的深度分布.兩對黑色粗箭頭代表區域主壓應力方向(Zhao et al., 2013).右上角插圖為騰沖火山區的相對位置.黑色實線代表斷層(Yang et al., 2013)，F1：怒江斷層，F2：龍嶺—瑞麗 斷層，F3：騰沖火山斷層，F4：盈江—隴川斷層.Fig.1 The geological setting and distribution of seismic stations The permanent and temporary seismic stations are showed by red triangles and yellow triangles, and the seven volcanic clusters in the study area are indicated by hollow ellipses. The gray and blue dots indicate all earthquakes with magnitude less than 4 since 2008 and events selected for shear wave splitting analysis in this study. The bar graph at the bottom right shows the number of selected events. Two pairs of thick black arrows represent the direction of regional compressive stress (Zhao et al., 2013). The inset shows the location of Tengchong Volcanic area in China. Black thin lines represent main faults (Yang et al., 2013). F1: Nujiang Fault, F2: Longling-Ruili Fault, F3: Tengchong volcano Fault, F4: Yingjiang-Longchuan Fault.

延遲時間反映了沿射線路徑的各向異性強度.我們利用歸一化延遲時間得到了各向異性百分比(A)(Savage，1999)，

(1)

其中d是震源到臺站的距離，VS是介質的橫波速度，研究區的地殼平均橫波速度為3.65 km·s-1(李永華等，2009).

2 結果

本文利用MFAST方法對15個地震臺站記錄的近震波形數據進行了橫波分裂分析，得到了1761個結果(A—D)，最終只采用了593個質量A，B級的高質量分裂結果(表1).從表1可以看出，15個臺站中總共有10個臺站的分裂結果多于10個，5個臺站的分裂結果少于10個.臺站TNC高質量的分裂結果數目最多(277個)；而流動臺站(L11,L12等)由于觀測時間較短，獲得的近震事件比較少，也沒有足夠的后方位角覆蓋(圖5)，最終的有效分裂結果比較少.但是除臺站L14外，其余每個臺站也都具有多于5個的分裂結果.

中上地殼各向異性離散現象顯著，因此對具有一定數量觀測結果的臺站進行討論才能降低隨機結果的影響.為提高討論的意義，本文只對具有10條及以上有效觀測結果的臺站加以后續分析.圖4為研究區具有10個以上有效分裂結果的臺站的快波偏振方向等面積投影玫瑰圖.從圖中可以看出多數臺站的偏振方向比較明顯，我們用統計平均的方式取得了這些臺站的優勢方向；個別臺站(MIZ，MZT,L13)得到的快波偏振方向比較分散，顯示兩個近乎正交的偏振方向，對于此類臺站需進行兩個快波偏振方向的分別統計.需要說明的是，盡管TNC、MAS等臺站具有足夠的剪切波分裂參數，但是我們并沒有發現其剪切波分裂參數隨時間存在明顯的變化.

圖2 橫波分裂分析計算示例(自動評級為A級)

圖3 橫波分裂分析計算示例(自動評級B級) 2017年1月26日QKT臺站記錄的2級地震的自動分裂結果. (a) 濾波后的E、N和Z三分量波形，垂直黑線表示手動拾取的S波到時，灰色區域代表剪切波分裂分析窗口，垂直虛線代表起始和結束時間; (b) 上面兩條波形代表向入射偏振方向和與其垂直的方向旋轉的波形，下面兩條波形是經過校正后的分量.虛線表示橫波分裂窗口的起始和結束時間范圍，從圖中可以看出切向分量上能量幾乎達到最小; (c) 從不同時間窗口確定的分裂參數，可看出該結果并不平緩穩定且存在一定的誤差; (d) 從不同的時間窗口確定的分裂參數分布，藍色十字代表最好的結果，可以看出該結果比較穩定; (e) 原始波形(左上)和校正后的波形(右上)(實線代表快波，虛線代表慢波)，和原始的的粒子運動(左下)和校正后的粒子運動(右下).如圖所示分裂校正后的波形擬合不太好，粒子運動軌跡由分裂前的形態變為近直線形態; (f) 協方差矩陣最小特征值的等值線圖，最佳結果(φ=23°，δt=0.06 s).用一個藍色的十字架標記，該等值線 不太收斂，說明結果不太穩定.Fig.3 An example of shear wave splitting analysis(rated as B automatically) Automatic splitting results for an earthquake with magnitude 2 recorded at station QKT on 26 January 2017. (a) Filtered waveform of E, N and Z component seismograms. Solid line represents the hand-picked S arrival. Dashed lines represent the minimum start and maximum end times，gray area is the time window for shear wave splitting analysis; (b) The upper two seismograms are waveforms that are rotated to the incoming polarization direction and the perpendicular direction. The bottom two seismograms are components after correction. Dashed lines represent the range of starting and ending time for shear wave splitting window. It can be seen that the energy in the tangential component is minimum after correction; (c) The splitting parameters determined from different time windows. It can be seen that the result is not stable and has certain error; (d) The distributions of splitting parameters determined from different time windows, the blue cross is the best result.The result is not stable; (e) Waveforms (top row) and particle motions (bottom row) for the original (left column) and corrected (right column) waveforms. Solid line represents fast wave, dashed line represents slow wave.The fast and slow wave is not fit after correction, the particle trajectory changed to a near-line shape after correction; (f) Contour plot of the smallest eigen value of the covariance matrix. The blue cross masks the best splitting parameter of this event (φ= 23°，δt=0.06 s). The contours are not converging, indicating that the results are not stable.

表1 各臺站基本參數與橫波分裂參數Table 1 Station parameters and results of shear-wave splitting by station

圖4 騰沖火山區具有10條以上有效記錄臺站的快波偏振方向玫瑰圖Fig.4 Rose diagrams of fast polarizations for each station with more than 10 measurements

如圖5所示，騰沖斷裂帶以西地區臺站(4個)的快波方向大多為近N-S向(如MIZ,RHT等)，位于騰沖斷裂帶西北部地區的臺站MIZ的各向異性偏振方向同與其相鄰的2個臺站MZT,L13的各向異性偏振方向都較為分散，三者其中的一個偏振方向較為一致，均為近NW-SE向.騰沖斷裂帶以東地區臺站(6個)的快波方向大多近NE-SW向(如臺站MZT,TNC,MAS等).

本研究得出的騰沖火山區的延遲時間在0.02 s到0.38 s之間變化，平均延遲時間為0.2 s.由表1可以看出，除臺站L14的延遲時間僅為0.02 s外，其他臺站的延遲時間均大于0.1 s.事實上，臺站L14只有2個高質量的測量結果，我們推測該臺站與其他臺站延遲時間結果的差異有兩種可能：(1)臺站L14附近地殼變形結構較弱，不存在明顯的各向異性；(2)有效分裂結果少，沒有全面反映臺站下方的各向異性特征.

圖5 研究區近震各向異性結果 每個臺站的快波偏振方向及大小用粗短棒表示，黑色短粗棒為本研究得出的分裂結果.紫色短粗棒為Zhang等(2018)給出的 分裂結果，藍色短粗棒為Shi等(2012)的分裂結果.Fig.5 Results of seismic anisotropy from shear wave splitting analysis of local earthquakes The polarization direction and size of the fast wave in each station are indicated by the thick short bars, and the black short thick bars are the splitting result in this study. The purple bars are the splitting result from Zhang et al. (2018). The blue bar is the splitting result from Shi et al. (2012).

前人(毛慧玲和秦嘉政，2011；太齡雪等，2015；Zhang et al.,2018；Shi et al.,2012)曾經利用云南地區固定和流動地震臺站的資料，采用近震剪切波分裂分析方法開展過各向異性研究工作.太齡雪等(2015)利用中國地震科學臺陣第一期流動地震臺站的近震波形資料，使用剪切波分裂系統分析方法SAM方法(高原等，2004)得到了云南地區的近震剪切波分裂參數，很好地揭示了研究區中上地殼的各向異性分布特征，但受觀測資料所限，其多數臺站的有效分裂結果都不足5個.Shi等(2012)同樣利用SAM方法對云南地區46個固定臺站進行了剪切波分裂分析，結果表明騰沖火山區內TNC和MAS臺站的快波偏振方向均為近N-S向，歸一化延遲時間均值為1.73 ms·km-1，與本研究得出的TNC,MAS的各向異性方向較為一致(圖5)，但延遲時間較本文得到的臺站歸一化延遲時間19.2 ms·km-1要小.毛慧玲和秦嘉政(2011)利用騰沖固定臺記錄的近震波形數據開展了剪切波分裂研究，結果顯示，騰沖臺下方的快波偏振方向主要集中在N30°E—N55°E 之間；Zhang等(2018)利用云南地區固定臺站的近震資料，采用MFAST方法開展了剪切波分裂分析，其結果顯示騰沖火山區內TNC和MAS兩個臺站的延遲時間與本文結果有較大差異，但快波偏振方向與本文較為一致(圖5).本研究得到的TNC和MAS兩個臺站的慢波延遲分別為0.23 s和0.29 s,而Zhang等(2018)采用同樣資料和方法給出的慢波延遲則不超過0.05 s.一些學者(如Gerst and Savage, 2004; Johnson et al., 2011; Keats et al., 2011; Unglert et al., 2011; Araragi et al., 2015)在新西蘭、日本等火山區開展的近震各向異性研究顯示，臺站下方的慢波延遲一般在0.1～0.3 s左右.有些研究(如Keats et al., 2011)也報道過火山噴發前慢波延遲時間降低到小于0.1 s的現象，并認為這種小的慢波延遲時間可能與巖漿活動導致的裂隙中充滿流體有關.Crampin等(2002) 認為，當孔隙流體壓力接近最大水平應力且裂縫縱橫比減小時，延遲時間將接近零.鑒于目前沒有證據表明騰沖火山區正處于火山噴發的前期，我們認為本研究的慢波延遲結果是可信的.

除了上述測量結果外，本研究在15個臺站中還獲得了93個“無效”的測量(圖6).S波在以下幾種情況下均會產生無效分裂：(1)S波穿過各向同性介質；(2)S波穿過各向異性介質時，快波偏振方向平行或者垂直于各向異性介質的快慢波方向；(3)臺站下方分布有復雜的各向異性介質，比如地下介質具有多個各向異性區域，兩個及其以上的各向異性區域被同時采樣而發生了雙折射使其各向異性相互抵消(Silver and Chan, 1991; Wustefeld and Bokelmann, 2007；強正陽和吳慶舉,2015;Saltzer et al., 2000; Menke andLevin, 2003; Levin et al., 2004).這幾類情況均會產生“無效”測量，即S波不會發生分裂且延遲時間接近為0.盡管S波的無效測量無法提供有關延遲時間的信息，但由于無效測量的事件方位與各向異性介質的快慢波的方向是平行或垂直的，因此，遠震剪切波分裂研究中，無效分裂通常也可以被用來去約束各向異性軸的方向(Silverand Chan, 1991; Wustefeld and Bokelmann, 2007).但與遠震SKS波不同的是，本文研究的近震Sg波的快波偏振方向依賴于地震震源機制，因為我們很難獲得騰沖火山區附近近震震源機制的相關信息，所以近震Sg波無效測量的意義還是不夠直觀，用無效測量去約束各向異性的方向可能存在一定誤差.

圖6 地震臺站剪切波分裂參數等面積圖 黑線代表分裂參數，藍點代表無效測量.Fig.6 Equal-area plots of splitting parameters (black lines) and null measurements (blue spots) for seismic stations Black lines represent the splitting parameters; blue spots represent the null measurements.

3 討論

3.1 各向異性強度與深度分布

由于我們使用的是近垂直入射的近震事件，因此可以利用震源深度來研究各向異性的深度分布.由圖7中可以看出快波偏振方向和延遲時間隨深度的變化并不明顯.本文又選取了有效記錄>50條的臺站RHT、TNC和MAS，分別統計了每個臺站不同深度地震事件與延遲時間的關系(圖8)，這3個臺站的結果并沒有表現出明顯的延遲時間隨深度的變化，進一步驗證了以上結論.如圖7c所示，47%的結果集中5～10 km處，各向異性強度基本都在20%以內；35%的結果集中在10～15 km處，各向異性強度在15%以內.15～20 km與20～25 km處的各向異性強度并沒有隨著地震事件數目的減少而減弱，說明各向異性強度的減弱與地震事件的數目無關，可能與震源深度的加深有關，研究區88%的結果都集中在0～15 km范圍內(圖7c)，因此推斷各向異性主要存在于淺層地殼內(<15 km)，是這一深度內各向異性疊加的結果.但是各向異性層的精確位置并非是精準的，因為各向異性強度是單位路徑上的延遲時間，反映的是整個路徑的平均各向異性(Maher and Kendall，2018).假設路徑上的各向異性是恒定的，那么騰沖火山上地殼的平均各向異性為0.016 s·km-1.

3.2 與前人Pms和P波層析成像等各向異性結果的比較

前人采用接收函數(孫長青等，2013)、遠震剪切波分裂(Shi et al.,2012;Huang et al., 2015;Zhao et al., 2013)、體波走時成像(Huang et al., 2018)等不同的手段,得到了研究區殼幔不同深度的各向異性分布特征.下面將對本文研究結果與前人各向異性研究結果進行對比分析.

Huang等(2018)利用區域地震P波走時資料開展了層析成像研究，得到了青藏高原東南緣地殼上地幔頂部不同深度的各向異性分布特征.前述研究表明，本文研究近震各向異性的主要貢獻來自15 km

圖7 各向異性快波偏振方向(a)、分裂延遲時間(b)和各向異性強度(c)隨深度的分布Fig.7 Distributions of fast wave polarization direction (a), delay time (b) and anisotropy percentage (c) of all high-quality results with depth

圖8 3個典型臺站延遲時間隨震源深度分布圖Fig.8 Distribution of delay time with focal depth for 3 seismic stations

以上的中上地殼，為此，我們將本文各向異性結果與10 km深度P波各向異性結果(Huang et al., 2018)進行了比較(圖9).10 km深度的P波各向異性結果顯示，騰沖火山區的各向異性快波偏振方向為NE-SW向，這與本研究采用近震剪切波分裂分析得到的快波偏振方向基本一致，但本文研究所揭示強烈的各向異性橫向變化在體波成像研究并沒有得到體現.產生這種差異的原因可能在于：(1)直達S波分裂研究可以給出很好的各向異性橫向變化，但其反映的各向異性是震源到接收臺站的介質各向異性的疊加效應，而不是某一深度處的各向異性特征；(2)P波成像(Huang et al., 2018)受到了到時、速度結構等因素的相互影響，而直達S波分裂受到了定位精度與S波讀取精度的影響，二者均具有一定的偏差性.

孫長青等(2013)利用固定臺站的地震資料，通過接收函數對Pms進行分裂分析，得到了川滇地區地殼的各向異性分布特征，其結果顯示騰沖地區全地殼的慢波延遲時間為0.3 s，而本文研究顯示，騰沖火山區中上地殼的平均慢波延遲時間為0.2 s，這說明騰沖火山區地殼各向異性的主要貢獻源自中上地殼.前人研究表明，騰沖地區具有高熱流、低電阻率、低波速等特征(Wang and Huangfu，2004；Xu et al., 2018; 樓海等，2002；胥頤等，2012)，并推測騰沖火山區下方10～25 km深度范圍存在巖漿囊.據此，我們可以推測，騰沖地區中上地殼介質存在顯著的變形作用與過程.

圖9 研究區不同地殼各向異性結果 每個臺站的快波偏振方向及大小用粗短棒表示.黑色短粗棒為本研究得出的分裂結果.白色細棒代表10 km深度P波各向異性成像結果(Huang et al., 2018)，紫色細棒代表Pms分裂結果 (孫長青等，2013).Fig.9 Comparisons of crustal anisotropy results from various analyses The polarization direction and size of the fast wave in each station are indicated by the thick short bars, and the black short thick bars are the splitting result in this study. The white thin bar represents the P wave anisotropic result at 10 km from body wave tomography (Huang et al., 2018), and the purple thin bar represents the result of Pms splitting (Sun et al. 2013).

3.3 地殼各向異性的成因

以往的研究表明，在遠離變形結構的區域，淺層地殼各向異性是由應力導致的裂隙的定向排列所致(Crampin,1978; Peacock et al., 1988; Rial et al., 2005)，區域應力場控制地下裂隙的方向使其平行于最大主壓應力的方向，最終使得地下各向異性介質的快波偏振方向平行于最大主壓應力方向，所以這類中上地殼的各向異性主要是由應力控制的.但在距離變形結構較近的區域，其中上地殼各向異性主要是定向排列的礦物所致，是由結構引起的各向異性.此外，在斷層附近的區域，由應力導致的中上地殼各向異性方向也有可能因受到附近斷層控制作用比較強而與斷裂帶平行(Leary et al., 1990).在前人的研究中存在類似的結果(Balfour et al., 2005; Boness and Zoback, 2006; Hurd and Bohnhoff, 2012; Li and Peng, 2017)，即斷裂帶附近區域的各向異性快波偏振方向平行于附近斷層的走向.

為了研究騰沖火山區中上地殼各向異性的主要成因，我們將本文各臺站快波偏振方向結果與斷層走向、應力分布等進行了比較(圖10)，另外，還以騰沖斷裂為中心，分別繪制了其東西兩區所有臺站記錄的地震事件的快剪切波偏振方向等面積投影玫瑰圖.在本研究中，我們使用的應力數據是前人通過反演震源機制解得到的(Zhao et al., 2013).如圖10所示，整個區域分布有兩個方向不同的主壓應力方向，大致為近N-S與近NE-SW向，表現出以騰沖斷裂帶為中心東西分區的特征.騰沖斷裂帶以西地區的應力方向主要為近N-S向，存在極少量NE-SW向，與該區大部分臺站(4個)的快波偏振方向存在很大的一致性：比如MIZ,RHT臺站的快波偏振方向均為近N-S向，與附近N-S向的應力方向幾乎平行；RST等臺站的快波偏振方向為近NE-SW向，與附近NE-SW向的應力方向也近乎平行.除此以外我們也發現，斷裂帶以西地區總體近N-S向的應力與我們在該地區記錄的所有地震事件的快波偏振方向也存在很大的一致性(圖10)，說明騰沖火山區以西地區中上地殼的各向異性主要是由應力引起的.騰沖斷裂帶以東地區的應力數據相對較少，但總體顯示出偏NE-SW向的主壓應力方向，不僅與我們在該地區記錄的所有地震事件的快波偏振方向一致，也與附近大部分臺站(7個)的快波偏振方向比較一致，比如MAS、MZT、SBT等臺站快波偏振方向均為近NE-SW向，與附近NE-SW向的主壓應力方向非常一致，說明騰沖火山區以東地區中上地殼的各向異性主要也是由應力控制的.

除此之外，位于騰沖斷裂帶北部地區的MIZ、MZT、L13臺站的各向異性偏振方向比較分散，均具有偏NW-SE向的第二個偏振方向，但該方向無論是與附近應力場方向或是與附近斷層方向都存在很大的差異，說明這類區域的各向異性可能受到了應力以及斷層的雙重影響，也有可能是由附近地表地形、地下結構構造造成的(Leary et al., 1990; Crampin and Peacock,2008)，即該區附近可能存在走向與其快波偏振方向接近的隱伏斷裂帶，有可能是隱伏斷裂帶控制了該地的各向異性.

對騰沖火山斷裂帶東西兩側中上地殼各向異性的分析表明，騰沖火山區中上地殼各向異性主要是區域主壓應力場引起的，該應力導致中上地殼微裂隙定向排列，產生平行于主壓應力方向的快波方向.

4 VP/VS

本研究中，我們根據公式(2)(Nur，1972)，利用2008-10—2019-07騰沖火山周圍10個固定臺站記錄的近震走時信息計算了VP/VS，并根據波速比與泊松比之間的關系, 將波速比轉化為泊松比(表2).

圖10 研究區各臺站偏振方向與附近主壓應力場方向的比較 每個臺站的快波偏振方向及大小用黑色粗短棒表示,粉色細棒代表最大水平主壓應力方向(Zhao et al., 2013)， 其中的玫瑰圖為騰沖斷裂帶東、西兩側所有臺站的快剪切波偏振方向等面積投影.Fig.10 Comparison of fast polarization direction and maximum main compressive stress direction The black short thick bars is the splitting results in this study, and the pink thin bars represent the maximum main compressive stress direction (Zhao et al., 2013). Rose diagrams show the fast polarizations for all seismic stations on the east and west sides of the Tengchong fault zone.

(2)

其中,T0是地震發生的起始時刻,TP是P波到時,TS是S波到時,tp=TP-T0，ts=TS-T0，在這里每條波形的波速比的值最大不能超過2.2，超過該值的記錄將會被舍棄.本文計算得該地區的平均VP/VS為1.68±0.01， 平均泊松比為0.21(表2).

已有研究表明，VP/VS值可為地殼巖石組成和流體含量的確定提供重要的約束(Christensen,1996).與全球地殼平均波速比(1.75)(Christensen,1996)和該地區地殼平均波速比(1.79～1.99)(李永華等,2009; 查小惠和雷建設, 2013)相比，本研究得到的該區平均VP/VS值相對較低.但該結果與世界其他火山區的上地殼波速比值相當(Nakajima et al., 2001; Kato et al., 2010).

通常認為，低的VP/VS意味著巖漿中的二氧化硅含量較高(55%～75%)(Christensen，1996)、流體含量小(～2%)(Nakajima et al., 2001)和/或有氣體存在.已有研究表明，騰沖火山區火山巖的SiO2

表2 固定臺站的VP/VS平均值與泊松比Table 2 Average VP/VS and Poisson′s ratio for permanent stations

含量介于48%～70%之間(向才英等，2000)、且存在豐富的地熱資源.前人研究(Wang and Huangfu, 2004；Lei et al., 2009; Xu et al., 2018; 樓海等，2002；胥頤等，2012)所揭示的騰沖地區高熱流、低電阻率、低波速等特征，被解釋為殼內巖漿囊.據此可以推測，研究區上地殼較低的VP/VS與淺表的地熱流體和/或深部充滿氣體的中酸性巖漿一致.

5 結論

本研究利用騰沖火山區15個固定和流動臺站記錄到的近震波形數據，采用MFAST方法進行了剪切波分裂分析，得到了該區近震各向異性剪切波分裂參數.由于本文研究首次使用了騰沖火山區8個火山臺站和5個流動臺站的數據，本文結果更清晰地揭示了騰沖火山區中上地殼各向異性的橫向變化特征.

研究結果顯示，騰沖火山區地震臺站下方的近震各向異性的慢波延遲時間為0.02～0.37 s，平均延遲0.2 s.各向異性強度隨深度的變化表明，0～15 km深度范圍淺層地殼介質對近震各向異性的貢獻最大.已有接收函數Pms分析(孫長青等,2013)表明，騰沖火山區地殼各向異性慢波延遲為0.3 s，據此推測研究區地殼各向異性的主要貢獻源自中上地殼.

本文結果顯示，研究區不同臺站的快波偏振方向變化很大，似乎反映了構造和區域應力場的共同作用.通過對分裂得到的快波偏振方向、應力方向和斷層走向的比較發現，以地震活動比較強烈的騰沖火山斷裂帶為界，其西側地區總體快波偏振方向呈近N-S向，而東側地區總體快波偏振方向呈NE-SW向，與Zhao等(2013)根據震源機制解給出的區域主壓應力方向基本一致，暗示研究區中上地殼各向異性主要是受主壓應力引起定向排列的裂隙所致.研究區低的地殼波速暗示騰沖火山區地殼應力場的局部變化可能是上地殼中富含氣體的中酸性巖漿膨脹所致.

此外，在騰沖火山區外圍部分臺站(如MIZ、MZT)的快波偏振方向比較分散，這意味著這類臺站附近可能存在走向與其快波方向接近的隱伏斷裂帶，即有可能是隱伏斷裂帶控制該地的各向異性，也可能因為這種臺站下上地殼各向異性受到了應力和構造的共同影響，因此表現出了較為復雜的分裂模式.盡管沒有揭示出各向異性隨時間變化這一現象，但本文圈定了最有可能受騰沖火山活動影響而導致各向異性參數變化的地震臺站，為未來騰沖火山活動監測和研究奠定了基礎.

致謝中國地震局地球物理研究所國家測震臺網數據備份中心(doi:10.11998/SeisDmc/SN)為本研究提供地震波形.文章大部分圖件使用GMT軟件繪制(Wessel and Smith,1998).兩位審稿人對提出的意見和建議對本文的提高有很大幫助，在此表示感謝.