垂直形變趨勢累計率在川滇地區的應用1

張　超 王世進 王　寧 王生文 徐佳垚

?

垂直形變趨勢累計率在川滇地區的應用1

張超 王世進 王寧 王生文 徐佳垚

（中國地震局第一監測中心，天津 300180）

本文在研究形變趨勢累計率的基礎之上，對數據處理方法進行了局部改進，將跨斷層場地布設較多的川滇地區作為實驗區，利用精密水準測量分析該地區的垂直形變異常強度，并結合川滇地區震例開展研究。分析發現強震發生前后一段時期，震中附近形變異常強度幾乎保持一致，但在震前（1、2個月左右）會發生較大改變，震后逐步恢復到與震前相同的情況；強震主要發生在形變異常強度值變化較大區域附近。垂直形變趨勢累計率在中短期的映震效果較好。

垂直形變趨勢累計率跨斷層測量精密水準形變異常強度

引言

跨斷層形變測量是通過觀測布設在斷層兩側控制點的相對關系，以監測斷層（上、下盤）的相對運動。由于控制點布設的位置，跨斷層測量能更容易、更直接地觀測到斷層的相對運動情況。很多學者對跨斷層測量進行了研究。張超（1981）將斷層上、下盤作為剛性塊體，從幾何學的角度進行了研究，利用基線觀測數據、水準觀測數據得到斷層三維運動情況。薄萬舉等（1998）在此基礎上，對八寶山斷裂帶上的跨斷層場地數據進行了處理，對該地區的斷層運動情況進行分析，發現在一定區域內跨斷層測量對異常形變監測和震情預測研究有重要意義。張興飛（1992）針對跨斷層形變測量的特點，構建了基于混合回歸理論的動態數據處理模型，該模型對地震前兆有較好的反映。焦青等（1998）利用華北和南北地震帶跨斷層流動水準和定點水準觀測資料，研究了斷層垂直位移速率異常的空間分布以及震源區的逼近情況，發現通過研究斷層垂直位移速率預測地震是一種可行方法。張希等（2013）將自適應控制理論引入斷層流動形變分析中，通過構建慢時變前兆系統結構，對5個典型場地資料進行分析，發現其異常識別、異常分類的判定較原始曲線更為直觀，與以往的形變指標相比也有所改進。張晶等（2011）從斷層受力運動方向入手，對斷層在一定外力作用下的運動情況進行了不同的組合，形成斷層活動協調比理論，分析了斷層活動性質和應變積累狀態，對地震中長期預測有一定的參考意義。

基于跨斷層測量的研究基本上都是通過計算斷層的三維運動情況，然后對其運動規律進行發掘，發現異常，這些方法的最終目的都是為了尋找或確定某條斷層或某個場地附近斷層的運動變化機理，進而開展地震預測預報工作。本文在該領域專家研究的基礎上，深入研究垂直形變趨勢累計率（陳兵等，2000），對其反映的物理機制以及區域應變場在中強震發生前后的變化規律進行了進一步研究，并得到了一些有用的結論。

1　基本原理

垂直形變趨勢累計率是通過垂直形變時變曲線建立的一種可較明確地量化趨勢變化的判定方法。此方法以斷層總體速率變化絕對量為基礎，再根據斷層總的垂直活動量判斷斷層運動趨勢，而且具有表征孕震應力場時空動態演化的物理意義。

首先通過公式（1）獲得斷層的相對運動速率（上盤相對下盤，上升為正）：

表示觀測值序列號，t為觀測序列所對應的時間，y表示t時刻上盤相對下盤觀測的高差值，v表示t時刻相鄰2期斷層的相對垂直形變速率。

為了排除干擾因素，更真實地計算斷層上、下盤的相對運動速度，本文在計算斷層相對 運動速率時，首先將不同年份（序號為）相同月份（序號為）的速率值取平均得到，然后再在此基礎上，用每一年每一期的速率值（）減去每一期對應的平均值獲取，稱之為總體性去周期，其中為年。具體見式（2）、（3）：

其中為參與計算測線的總觀測年數。

該想法啟發于大地水準面的確定。與相鄰年同期觀測數據做差或較短觀測期間不同時間段（年）中相同期（年中月）取平均值獲取相對運動速率的傳統方法相比，這樣做能夠更好地消除因環境等因素引起的長短周期影響。

總體去周期后，那么以起算時刻（序號為，場地起始觀測年份<起算年份<最后觀測年份）至計算時刻（序號為、0所對應年份和相應的月份）這一時間段斷層總體垂直速率變 化量如下：

其中，一般情況下=-10。

依據（4）、（5）式，計算垂直形變趨勢性的量化標準，即（6）式，垂直形變趨勢累計率。

如果斷層運動速率v為恒值，那么==0，即c為0/0型。為使c與實際物理情況相符，此情況下規定c=1。

另外，速率絕對幅度的大小及積累時間的長短與相應的地震強度和復發周期有關。強震孕育到一定程度后，會出現閉鎖、微破裂、預滑移等現象，必然對斷層形變的大小和方向產生影響，從而反映到c值的異常變化。

2　跨斷層觀測場地情況

基于以上原理，本文以川滇地區為實驗區，利用布設在重要斷裂帶上的跨斷層場地的水準觀測數據開展實驗，并進行了分析。

以下為川滇地區跨斷層流動場地水準資料觀測情況。各個場地、各條測線的觀測情況如表1，川滇地區跨斷層流動觀測場地分布情況如圖1。

表1　川滇地區在測流動跨斷層場地情況表

續表

①說明：例如下關測段12，即表示從上盤控制點1到下盤控制點2；2x-1x表示上盤2號新控制點到下盤新1號控制點，x代表新點。

上述場地中部分測線因場地環境改變或觀測墩（點）遭破壞等情況而停測或更改觀測路線的測段，以及未跨過斷層的測段并未列出。另外由于一些新建場地的觀測時間距現在較近，故這些場地的相關信息也未列出。

川滇地區跨斷層觀測場地多數建于20世紀80年代，主要分布在鮮水河斷裂帶、龍門山斷裂帶、安寧河-則木河斷裂帶以及小江斷裂帶等幾大斷裂帶上，其場地分布見圖1。場地位置用五角星標注。

3　實驗分析

川滇地區的觀測場地主要位于南北地震帶上，而南北地震帶是我國地震高發地區，擁有豐富的震例（劉辛中等，2015），因此作者選取該地區地震為樣本進行分析。另外為彌補云南、四川跨斷層場地觀測周期不同以及因外界環境等因素造成的數據不全對計算的影響，在數據預處理過程中對數據進行了3次樣條插值處理。其中，為獲取地震前后川滇地區垂直形變趨勢較為真實的變化情況，盡量選擇發震時間前后各個場地都有觀測數據的時間作為實驗時間，對受外界條件影響或在核實異常過程中認為存在質量問題的觀測數據不予采用，也盡量避免因插值造成的不必要誤差。

由于篇幅原因，基于第2章所列觀測場地的觀測資料僅對汶川地震、蘆山地震前后c值變化所反映的規律分別進行討論，其中起算時間分別選為1998、2003年。

圖2中，（a）與（d）分別為汶川地震發生前4個月、震后2個月的垂直形變趨勢累計率，（b）與（c）為地震發生前2個月、震后1個月的垂直形變趨勢累計率。可以看出，汶川地震發生在垂直形變趨勢變化較大的區域附近，而地震前后，發震地附近的垂直形變趨勢發生較大改變（矩形框示意范圍），反映其附近斷層垂直形變異常強度改變較大。地震發生前2個月形變趨勢表現為運動平緩，即形變異常強度變小，而震后則逐漸恢復到與震前相同的情況。作者認為，這種情況說明發震地附近震前斷層的“正常”運動趨勢受阻，導致斷層附近出現更為明顯的閉鎖現象，為地震的發生提供了必要條件。相關文獻（趙靜等，2012）也給出了類似的觀點，進一步說明跨斷層水準能夠較有效地反映地震前兆。實驗過程中，作者發現在汶川地震發生后2個月（圖2（d））川滇地區的垂直形變異常強度與汶川地震發震前（圖2（b））的垂直形變異常強度極為相似，而通過查詢中國地震臺網中心發布的地震目錄可以發現，在汶川地震發生2個多月后，即2008年7月24日，汶川主震震中東北方約250km處的廣元發生了6.0級強余震（鄭勇等，2009；溫瑞智等，2009）。需要說明的是，川滇地區跨斷層觀測數據的獲取時間主要集中在每個觀測月的中上旬，故圖2（d）顯示的為廣元地震發生前該地區的垂直形變異常強度，這從另一個方面（形變學）說明跨斷層手段所反映的垂直形變異常強度具有較好的映震能力。

圖3中，（a）、（b）分別為該地區發震前6個月和1個月的垂直形變趨勢圖，（c）、（d）為震后1個月和2個月的垂直形變趨勢圖。從圖中可以發現蘆山地震的發震地同樣處于垂直形變趨勢變化較大（即垂直形變異常強度值變化較大，主要集中在震中東南方向，矩形框示意范圍）的地區附近。川滇地區2012年9月份與2013年5月、6月的垂直形變趨勢累計率幾乎保持一致，而該地區2013年3月較2012年9月的垂直形變累計率發生了較大改變，即斷層垂直形變異常強度值改變較大，隨即2014年4月20日在該區域便發生S6.7蘆山地震。

總的來說，垂直形變趨勢累計率在以下2方面有較好的效果：①從地震發生的位置上來看，發震地一般都處于斷層形變異常強度值改變較大的地區，地震發生前斷層附近的垂直形變異常強度有突然變緩或加劇的趨勢；②從地震前后應變能變化的角度研究，地震發生前應變能大量累積，地震發生后得到釋放，但還處于平衡自適應階段，震后一段時間應變能達到平衡，使其與震前相當，較好地反映了地震孕震、發震以及震后平靜期的應變能變化過程（許昭永等，2015）。

4　討論

通過分析川滇地區的震例，得到了幾方面認識：從時間角度看，地震發生前后，川滇地區跨斷層場地附近垂直形變異常強度幾乎維持不變；從空間位置角度看，地震（中強震）主要發生在垂直形變異常強度較大地區或者其邊緣及附近；從震級角度來看，較大震級的地震在震前的異常形變強度值變化比較小震級地震更為明顯，但目前缺少更為確切的依據，需要開展進一步研究。總體來說，垂直形變趨勢累計率在時間和空間上對中短期地震預測預報能夠發揮不錯的效果，但對于短臨預報，由于數據采集頻率以及跨斷層場地分布情況的影響，該手段的有效性有待進一步驗證。另外，該方法需要足夠的監測場地密度，否則僅從單個跨斷層場地的垂直形變趨勢變化速率以及趨勢變化幅度分析，會因缺少斷層運動的整體性判斷而更加難以獲得真實的空間異常變化特征，因此對地震預測預報的幫助可能較小。地震發生過程是應變能積累-釋放-再積累的過程，因此在選用垂直形變趨勢累計率方法時，起算時間的選擇也是不可忽略的因素之一。

致謝：感謝2位審稿人對本文的悉心指導，提出的建議對本文質量提升有很大幫助。本文所使用的跨斷層觀測數據為一線觀測人辛苦所得以及兄弟單位慷慨提供，在此一并感謝。

薄萬舉，謝覺民，郭良遷，1998．八寶山斷裂帶形變分析與探討．地震，18（1）：63—68．

陳兵，江在森，胡斌，2000．斷層形變趨勢變化研究——垂直形變趨勢累積率c值．地震地磁觀測與研究，21（1）：50—55．

焦青，周俊萍，1998．跨斷層垂直位移速率的動態演化特征與強震的關系．地震，18（3）：265—273．

劉辛中，劉志廣，朱爽等，2015．安寧河-則木河斷裂帶斷層活動特征及其與地震活動的關系．震災防御技術，10（2）：253—261．

溫瑞智，周正華，李小軍等，2009．汶川S8.0地震的強余震流動觀測．地震學報，31（2）：219—225．

許昭永，文麗敏，石寶文等，2015．川滇菱形地塊的應變能積累釋放周期和強震預測．地震學報，37（5）：774—786．

張超，1981．跨斷層位移測量與局部斷層活動．地震，1（5）：24—25．

張晶，黎凱武，武艷強等，2011．斷層活動協調比在地震預測中的應用．地震，31（3）：19—26．

張希，唐紅濤，賈鵬等，2013．跨斷層形變的慢時變結構分析與異常識別初探．武漢大學學報·信息科學版，38（6）：710—715．

張興飛，1992．跨斷層形變測量動態數據模型的研究．武漢測繪科技大學學報，17（2）：79—85．

趙靜，江在森，武艷強等，2012．汶川地震前龍門山斷裂帶閉鎖程度和滑動虧損分布研究．地球物理學報，55（9）：2963—2972．

鄭勇，馬宏生，呂堅等，2009．汶川地震強余震（S≥5.6）的震源機制解及其與發震構造的關系．中國科學D輯：地球科學，39（4）：413—426．

The Application of the Accumulation Rate of Vertical Deformation Trend in the Sichuan-Yunnan Region

Zhang Chao, Wang Shijin, Wang Ning, Wang Shengwen and Xu Jiayao

(The First Crust Monitoring and Application Center, CEA, Tianjin 300180, China)

Based on the study of the cumulative rate of vertical deformation trend, the data processing method is modified partially, and the Sichuan-Yunnan region with more cross-fault site is taken as experimental area. The intensity of vertical deformation anomalies from precision leveling is analyzed. The results show that the intensity of deformation anomalies near the epicenter remains almost the same in a period before and after the occurrence of strong earthquakes, but the great change may take place one and two months before and after the earthquake. Then, it would gradually recover to the original situation. We also found that strong earthquakes occur mainly near the regions where the values of anomalous intensity change greatly. The reflect effect of accumulative rate of the vertical deformation tendency is good in the short to medium term.

The accumulation rate of vertical deformation trend; Cross-over fault measurement; Precise leveling; Intensity of deformation anomalies

1基金項目 科技部科技基礎性工作專項（2015FY210400），中國地震局第一監測中心主任基金（FMC2015006），中國地震局監測、預測、科研三結合課題（CEA-JC/3JH163703）共同資助

2016-07-11

張超，男，生于1987年。碩士研究生，助理工程師。主要從事區域地殼形變數據處理與方法等研究。E-mail：875011229@qq.com

張超，王世進，王寧，王生文，徐佳垚，2017．垂直形變趨勢累計率在川滇地區的應用．震災防御技術，12（2）：329—337． doi：10.11899/zzfy20170209