垂直形變趨勢累計率在川滇地區(qū)的應(yīng)用1

張　超 王世進(jìn) 王　寧 王生文 徐佳垚

?

垂直形變趨勢累計率在川滇地區(qū)的應(yīng)用1

張超 王世進(jìn) 王寧 王生文 徐佳垚

（中國地震局第一監(jiān)測中心，天津 300180）

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

垂直形變趨勢累計率跨斷層測量精密水準(zhǔn)形變異常強(qiáng)度

引言

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

基于跨斷層測量的研究基本上都是通過計算斷層的三維運動情況，然后對其運動規(guī)律進(jìn)行發(fā)掘，發(fā)現(xiàn)異常，這些方法的最終目的都是為了尋找或確定某條斷層或某個場地附近斷層的運動變化機(jī)理，進(jìn)而開展地震預(yù)測預(yù)報工作。本文在該領(lǐng)域?qū)＜已芯康幕A(chǔ)上，深入研究垂直形變趨勢累計率（陳兵等，2000），對其反映的物理機(jī)制以及區(qū)域應(yīng)變場在中強(qiáng)震發(fā)生前后的變化規(guī)律進(jìn)行了進(jìn)一步研究，并得到了一些有用的結(jié)論。

1　基本原理

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

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

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

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

其中為參與計算測線的總觀測年數(shù)。

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

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

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

依據(jù)（4）、（5）式，計算垂直形變趨勢性的量化標(biāo)準(zhǔn)，即（6）式，垂直形變趨勢累計率。

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

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

2　跨斷層觀測場地情況

基于以上原理，本文以川滇地區(qū)為實驗區(qū)，利用布設(shè)在重要斷裂帶上的跨斷層場地的水準(zhǔn)觀測數(shù)據(jù)開展實驗，并進(jìn)行了分析。

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

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

續(xù)表

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

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

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

3　實驗分析

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

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

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

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

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

4　討論

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

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

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

陳兵，江在森，胡斌，2000．?dāng)鄬有巫冓厔葑兓芯俊怪毙巫冓厔堇鄯e率c值．地震地磁觀測與研究，21（1）：50—55．

焦青，周俊萍，1998．跨斷層垂直位移速率的動態(tài)演化特征與強(qiáng)震的關(guān)系．地震，18（3）：265—273．

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

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

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

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

張晶，黎凱武，武艷強(qiáng)等，2011．?dāng)鄬踊顒訁f(xié)調(diào)比在地震預(yù)測中的應(yīng)用．地震，31（3）：19—26．

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

張興飛，1992．跨斷層形變測量動態(tài)數(shù)據(jù)模型的研究．武漢測繪科技大學(xué)學(xué)報，17（2）：79—85．

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

鄭勇，馬宏生，呂堅等，2009．汶川地震強(qiáng)余震（S≥5.6）的震源機(jī)制解及其與發(fā)震構(gòu)造的關(guān)系．中國科學(xué)D輯：地球科學(xué)，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基金項目 科技部科技基礎(chǔ)性工作專項（2015FY210400），中國地震局第一監(jiān)測中心主任基金（FMC2015006），中國地震局監(jiān)測、預(yù)測、科研三結(jié)合課題（CEA-JC/3JH163703）共同資助

2016-07-11

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

張超，王世進(jìn)，王寧，王生文，徐佳垚，2017．垂直形變趨勢累計率在川滇地區(qū)的應(yīng)用．震災(zāi)防御技術(shù)，12（2）：329—337． doi：10.11899/zzfy20170209