利用GNSS資料分析云南地區(qū)M≥5.7地震前應變場分布特征

楊建文 葉 泵 陳 佳 高 瓊 張華英

1 中國地震科學實驗場大理中心，云南省大理市濱海大道，671000 2 云南大理滇西北地殼構造活動野外科學觀測研究站，云南省大理市濱海大道，671000

地震前兆是近年來各國進行地震預報的基本途徑之一[1]。從觀測事實出發(fā)來探尋地震預報的可能性，是現(xiàn)階段經驗預測的基礎[2]。而GNSS觀測以其物理意義清晰、資料可靠、描述大地形變精度高等特點，被廣泛應用于地震預測預報。

GNSS揭示的地殼形變動態(tài)分布特征與強震的關系一直是研究熱點[3]。江在森等[4]基于GPS計算結果，對應變場空間分布與強震的關系進行研究，認為強震通常發(fā)生在最大剪應變率的高值區(qū)或邊緣地帶；洪敏等[5]基于云南省早期的28個GNSS連續(xù)站數據，對地震事件與地殼形變之間的關系進行研究認為，當部分點位運動方向背離長期趨勢運動背景時，發(fā)生地震的危險性較大；付虹等[2]對2014年發(fā)生的3次M≥6.0地震前的面應變分布特征進行研究認為，構造區(qū)局部擠壓增強可能是地震短臨階段區(qū)域構造變形最重要的特征之一。

隨著中國大陸構造環(huán)境監(jiān)測網絡、云南省政府10項重點工程、中國地震科學實驗場等項目的推進，云南地區(qū)形變觀測能力得到顯著增強，為地震機制和物理機制研究、地震異常識別等提供了可靠的數據基礎。實驗表明，強震的孕育、發(fā)生及震后調整的過程中通常會伴隨較為明顯的地殼形變，對強震前的地殼形變特征進行研究，可加深對地震孕育過程的認識。本文以2012年以來云南地區(qū)發(fā)生的7次M≥5.7地震為樣本，基于云南地區(qū)GNSS連續(xù)觀測資料，分別對7次地震前的應變場分布特征進行研究，總結其共性和差異。該工作的開展對地震前兆探測及預測預報工作具有重要意義。

1 數據和方法

GNSS資料來源于中國地震局GNSS數據產品服務平臺提供的站點位移時間序列原始數據，研究時段為2012-01～2018-12，研究范圍(21°～29°N，97°～107°E)內共包含43個GNSS連續(xù)觀測站。研究時段內共發(fā)生7次M≥5.7地震(剔除余震)，相關地震參數見表1，震中分布及GNSS站點分布見圖1。由圖1可知，7次M≥5.7地震主要位于滇西北、滇東北和滇南等地區(qū)。由表1中震源機制解可知，7次地震主要為走滑型地震。

表1 地震相關參數

圖1 GNSS連續(xù)觀測站、地震震中分布及格網劃分結果Fig.1 The distribution of GNSS continuous observation stations,earthquake epicenters and grid division

考慮到GNSS連續(xù)觀測站分布的不均勻性，首先采用克里金插值法對點位數據進行1°×1°格網化插值(圖1)，獲取均勻分布于80個格網的位移時間序列，然后通過位移時間序列求取區(qū)域應變場[6-8]。應變場具體解算方法詳見文獻[6]。

通過應變場解算可獲取3個狀態(tài)分量：反映東西與南北方向拉張與收縮的正應變εx和εy、東西與南北方向的剪應變γxy[4,6]。在此基礎上，計算9個應變參數，包括面應變、最大剪應變、主應變等[4]。其中，面應變θ和最大剪應變γmax的計算公式為[4,6]：

(1)

式(1)可確定區(qū)域內任意時段的面應變、最大剪應變參數。截取任意起止時間可獲取該時段內的累積應變量，通過滑動結束時間可進一步獲取應變場變化結果[6,8]。

2 震前應變場分布特征

應變張量是用來反映純形變的參數，能夠全面表征形變的性質與強度[4]。應變參數實際上是對不同方式形變信息的一種客觀分離，研究地殼形變與強震的關系，需要分別對多個應變參數的空間分布進行研究[4]。面應變可直接表征區(qū)域擠壓(收縮)或拉張(膨脹)的強弱，反映該區(qū)域的應變累積狀態(tài)，而最大剪應變是反映剪切應變強度的參數。震源機制解研究表明，云南地區(qū)7次M≥5.7地震多為走滑型地震，地震破裂類型多含剪切破裂，震前應變累積中包含一定量的剪切應變累積，因而強震發(fā)生前最大剪應變會表現(xiàn)出一定的異常特征。因此，分別提取面應變和最大剪應變信息可揭示云南地區(qū)M≥5.7地震前的應變場分布特征。

以7次強震前1 d為截止時間，分別計算震前6個月(180 d)的面應變、最大剪應變累積變化值，繪制應變分布圖(圖2、圖3)?？紤]到寧蒗M5.7和彝良M5.7、盈江M6.1和魯甸M6.5、魯甸M6.5和景谷M6.6地震間的時間間隔小于6個月，獲取的后者應變累積變化值中除含有后一次地震孕震變化外，可能也包含前一次地震孕育-發(fā)生的影響過程。為更好地進行應變場動態(tài)變化分析，本文分別在寧蒗M5.7和彝良M5.7、盈江M6.1和魯甸M6.5、魯甸M6.5和景谷M6.6地震間(間隔時間中點)增加面應變和最大剪應變累積的變化結果。

圖2 云南地區(qū)M≥5.7地震前面應變分布Fig.2 Surface strain distribution before M≥5.7 earthquakes in Yunnan area

圖3 云南地區(qū)M≥5.7地震前最大剪應變分布Fig.3 Maximum shear strain distribution before M≥5.7 earthquakes in Yunnan area

2.1 面應變

圖2為7次M≥5.7地震前的面應變分布圖，圖中藍色區(qū)域為擠壓變形區(qū)，紅色為拉張變形區(qū)。由圖2(a)可知，在2012-06-24寧蒗5.7級地震前，在滇西南的下關-景東-墨江一帶存在較為明顯的面應變拉張區(qū)，累積面應變達4.0×10-8，寧蒗5.7級地震發(fā)生在顯著拉張區(qū)以北的面拉張變化高梯度帶上。由圖2(c)可知，在2012-09-07彝良5.7級地震前，云南地區(qū)同時存在顯著的面擠壓和面拉張區(qū)，在滇南的臨滄-思茅-勐海和彌勒-蒙自-金平一帶分別存在2個顯著的面擠壓區(qū)，最大累積面應變達-3.5×10-8；在麗江-永勝-攀枝花-寧南一帶也存在較為明顯的面擠壓區(qū)，累積面應變達-2.0×10-8，彝良5.7級地震就發(fā)生在面擠壓變化的梯度帶上。瑞麗-騰沖-施甸一帶存在顯著的面拉張區(qū)，最大累積面應變達4.0×10-8；廣西百色一帶也存在較為明顯的面拉張區(qū)，累積面應變達3.5×10-8。由圖2(a)～2(c)可知，從寧蒗M5.7地震到彝良M5.7地震，云南地區(qū)的面應變存在較為明顯的變化過程，其中滇南地區(qū)最為顯著，在彝良5.7級地震前形成強擠壓和強拉張并存的格局。此外，下關-景東-墨江一帶的面拉張也存在明顯的減弱趨勢。由圖2(d)可知，在2013-08-31德欽5.9級地震前，滇西北及鄰區(qū)的中甸-鄉(xiāng)城一帶存在較為明顯的面應變拉張區(qū)，累積面應變達2.0×10-8，德欽5.9級地震發(fā)生在面應變拉張區(qū)的高梯度帶上；麗江-云龍-施甸-耿馬一帶存在較為明顯的面擠壓區(qū)，累積面應變達-2.0×10-8。

由圖2(e)可知，在2014-05-30盈江6.1級地震前，云南地區(qū)面應變異常較為突出，滇西南及鄰區(qū)存在強擠壓和強拉張并存的格局，在瀾滄-思茅-勐海一帶存在顯著的面拉張區(qū)，累積面應變達9.5×10-8；瑞麗-耿馬一帶存在顯著的面擠壓區(qū)，累積面應變達-4.5×10-8，盈江6.1級地震則發(fā)生在面擠壓變化的高梯度帶上。此外，越西-冕寧-鹽源和東川-鳳凰山-通海一帶也存在較為明顯的面擠壓區(qū)，滇南的金平-蒙自-文山一帶存在較為明顯的面拉張區(qū)。由圖2(g)可知，在2014-08-03魯甸6.5級地震前，滇東北及鄰區(qū)的越西-喜德-寧南-會澤-東川-彌勒一帶存在明顯的面擠壓區(qū)，累積面應變達-3.5×10-8；滇南的金平-蒙自-文山一帶存在較為明顯的面拉張區(qū)，累積面應變達3.5×10-8，魯甸6.5級地震則發(fā)生在面擠壓變化的高梯度帶上。由圖2(i)可知，在2014-10-07景谷6.6級地震前，云南地區(qū)面應變累積顯著，滇西南存在強擠壓和強拉張并存的格局，瑞麗-騰沖一帶存在明顯的面拉張區(qū)，累積面應變達5.0×10-8；臨滄-思茅-勐海一帶存在明顯的面擠壓區(qū)，累積面應變達-4.5×10-8，景谷6.6級地震發(fā)生在面擠壓的高值區(qū)。另外，會澤-東川-彌勒和文山-蒙自-金平一帶分別存在2個較為明顯的面擠壓區(qū)，最大面應變?yōu)?3.0×10-8。由圖2(e)～2(i)可知，從盈江M6.1地震到魯甸M6.5地震再到景谷M6.6地震，云南地區(qū)的面應變存在較為明顯的動態(tài)演化過程，整個過程中面擠壓變化最為明顯，冕寧-寧南-會澤-東川-彌勒和景東-瀾滄-思茅一帶的面擠壓存在由弱到強的變化過程。滇西南的施甸-耿馬-瀾滄-勐海一帶的面應變經歷了強-弱-強的變化過程。由圖2(j)可知，在2018-09-08墨江5.9級地震前，滇西南存在強擠壓和強拉張并存的格局，瑞麗-騰沖一帶存在明顯的面拉張區(qū)，累積面應變達4.5×10-8；臨滄-瀾滄-勐海一帶存在明顯的面擠壓區(qū)，累積面應變達-4.5×10-8，墨江5.9級地震發(fā)生在面擠壓變化的高梯度帶上。此外，越西-喜德-寧南-會澤和鄉(xiāng)城-中甸-麗江一帶分別存在較為明顯的面擠壓區(qū)和面拉張區(qū)。

2.2 最大剪應變

圖3為7次M≥5.7地震前最大剪應變分布圖。由圖3(a)可知，在2012-06-24寧蒗5.7級地震前，滇東北及鄰區(qū)的冕寧-寧南-東川一帶和滇南的瀾滄-思茅-墨江-金平一帶存在明顯的最大剪應變高值區(qū)，累積最大剪應變分別達3.0×10-8和5.0×10-8，寧蒗5.7級地震就發(fā)生在冕寧-寧南-東川高值區(qū)最大剪應變變化的梯度帶上。由圖3(c)可知，在2012-09-07彝良5.7級地震前，鄉(xiāng)城-木里-冕寧-喜德-馬邊、瑞麗-騰沖-施甸、彌勒-通海-新平-墨江-思茅和廣西的百色等地存在明顯的最大剪應變高值區(qū)，累積最大剪應變達5.5×10-8，彝良5.7級地震則發(fā)生在第一個高值區(qū)最大剪應變變化的高梯度帶上。由圖3(a)～3(c)可知，從寧蒗M5.7地震至彝良M5.7地震，云南地區(qū)存在顯著的最大剪應變累積過程，其中以鄉(xiāng)城-木里-冕寧-喜德-馬邊和瑞麗-騰沖-施甸等地變化最為明顯。由圖3(d)可知，在2013-08-31德欽5.9級地震前，永勝-木里-喜德-寧南-東川-通海一帶存在明顯的最大剪應變高值區(qū)，累積最大剪應變達3.5×10-8，德欽5.9級地震則發(fā)生在最大剪應變變化的高梯度帶上。由圖3(e)可知，在2014-05-30盈江6.1級地震前，在滇西南的施甸-耿馬-瀾滄-勐海一帶存在明顯的最大剪應變高值區(qū)，累積最大剪應變達10.0×10-8，盈江6.1級地震則發(fā)生在最大剪應變變化的梯度帶上。此外，會澤-東川-彌勒-文山和鄉(xiāng)城-中甸-麗江一帶也存在較為明顯的最大剪應變高值區(qū)。

由圖3(g)可知，在2014-08-03魯甸6.5級地震前，會澤-東川-鳳凰山-通海-蒙自一帶存在明顯的最大剪應變高值區(qū)，累積最大剪應變達5.0×10-8，魯甸6.5級地震則發(fā)生在最大剪應變變化的高梯度帶上。同時，鄉(xiāng)城-木里-冕寧-喜德-馬邊一帶也存在明顯的最大剪應變高值區(qū)。由圖3(i)可知，在2014-10-07景谷6.6級地震前，瑞麗-騰沖-施甸-耿馬一帶存在明顯的最大剪應變高值區(qū)，累積最大剪應變達6.0×10-8，景谷6.6級地震就發(fā)生在最大剪應變變化的高梯度帶上。另外，木里-冕寧-越西和元謀-鳳凰山-通海-蒙自-文山一帶存在較為明顯的最大剪應變高值區(qū)。由圖3(e)～3(i)可知，從盈江M6.1地震到魯甸M6.5地震再到景谷M6.6地震，云南地區(qū)的最大剪應變存在較為明顯的動態(tài)演化過程，瑞麗-騰沖-施甸-耿馬一帶的最大剪應變經歷了強-弱-強的變化過程。值得注意的是，整個變化過程中木里-冕寧-越西和元謀-鳳凰山-通海-蒙自-文山一帶的最大剪應變累積一直存在。由圖3(j)可知，在2018-09-08墨江5.9級地震前，瑞麗-騰沖-施甸-耿馬一帶存在明顯的最大剪應變高值區(qū)，累積最大剪應變達6.0×10-8，墨江5.9級地震就發(fā)生在最大剪應變變化的高梯度帶上。此外，木里-鹽源-寧南-會澤一帶也存在較為明顯的最大剪應變高值區(qū)，累積最大剪應變達4.0×10-8。

3 討 論

前文已對7次M≥5.7地震前的面應變和最大剪應變分布特征進行了分析，在此基礎上，對7次地震前的應變場信息作進一步總結(表2、表3)。

表2 云南地區(qū)M≥5.7地震前面應變相關信息

表3 云南地區(qū)M≥5.7地震前最大剪應變相關信息

由表2可知，7次M≥5.7地震前，云南地區(qū)存在顯著的面應變變化趨勢，多數情況下存在強擠壓和強拉張并存的格局，除德欽5.9級地震和魯甸6.5級地震外，其他地震前累積面應變均超過±4.0×10-8。從震前累積面應變量值與地震震級之間的對應關系看，兩者之間并無明顯關聯(lián)。從地震發(fā)生的位置看，除彝良5.7級地震和景谷6.6級地震外，其他地震均發(fā)生在面應變變化的高梯度帶上，7次M≥5.7地震中有5次發(fā)生在面擠壓區(qū)，占比71%；有2次發(fā)生在面拉張區(qū)，占比29%；3次M≥6.0地震均發(fā)生在面擠壓區(qū)。因此，面應變(特別是面擠壓)變化的高梯度帶可作為研判未來M≥5.7地震發(fā)震地點的重要區(qū)域。

由表3可知，7次M≥5.7地震前，云南地區(qū)存在顯著的最大剪應變變化趨勢，除德欽5.9級地震外，其他地震前累積最大剪應變均超過5.0×10-8。值得注意的是，震前累積最大剪應變量值與地震震級并無直接關聯(lián)。從地震發(fā)生的位置看，多數地震發(fā)生在最大剪應變變化的高梯度帶上。因此，最大剪應變變化的高梯度帶可作為研判未來M≥5.7地震發(fā)震地點的重要區(qū)域。

面應變和最大剪應變都是重要的應變場參數，基于反映不同形變信息的多種應變參數來分析應變場空間分布與強震地點的關系，可避免單一應變結果的局限性。研究表明，7次M≥5.7地震多發(fā)生在面應變(特別是面擠壓)和最大剪應變變化的高梯度帶上，這些區(qū)域通常會伴隨較大的應變場變化信息，強震發(fā)生的可能性更大。

強震的孕育與發(fā)生過程復雜多樣，不同類型和相同類型的地震在不同介質環(huán)境中的孕育過程也存在差異，目前對強震孕育過程的認知還十分有限，汶川大地震使我們意識到從大時空尺度來獲取大地震孕育過程中異常信息的重要性[9]。GNSS連續(xù)觀測站數量的增加以及觀測資料和震例的進一步積累，可為獲取強震前更加精細的應變場分布特征和總結更多震例提供理論基礎，這對于進一步探索GNSS資料在強震預測中的應用、提高異常識別能力、推進指標體系建設都具有重要意義。

4 結 語

1)7次M≥5.7地震前，云南地區(qū)存在顯著的面應變變化趨勢，多數情況下存在強擠壓和強拉張并存的格局，累積面應變一般超過±4.0×10-8，地震多發(fā)生在面應變(特別是面擠壓)變化的高梯度帶上。

2)7次M≥5.7地震前，云南地區(qū)存在顯著的最大剪應變變化趨勢，累積最大剪應變一般超過5.0×10-8，地震多發(fā)生在最大剪應變變化的高梯度帶上。

3)基于GNSS連續(xù)觀測數據獲取的面應變和最大剪應變參數物理意義明確，能定量表征區(qū)域擠壓或拉張強弱和區(qū)域剪切應變強度。其變化高梯度帶可作為研判未來M≥5.7地震發(fā)震地點的重要區(qū)域。

致謝：本文GNSS數據來源于中國地震局GNSS數據產品服務平臺(http:∥www.cgps.ac.cn)，數據分析軟件由云南省地震局洪敏高級工程師提供，在此一并表示感謝。