汶川地震震前與同震斷層滑移的比較分析1

梁向洲 郝圣旺

?

汶川地震震前與同震斷層滑移的比較分析1

梁向洲 郝圣旺

（燕山大學(xué)建筑工程與力學(xué)學(xué)院，河北秦皇島066004）

本文利用GPS觀測的1999—2007年汶川震前3期地表變形數(shù)據(jù)和2008年汶川同震地表變形數(shù)據(jù)，結(jié)合地震位錯理論，通過高斯變換和坐標(biāo)旋轉(zhuǎn)建立斷層模型，運用遺傳算法，反演了龍門山斷裂帶斷層震前3期和同震滑動參數(shù)。結(jié)果表明龍門山斷層震前3期平均走滑位移為-5.39mm，傾向位移為2.66mm，與同震斷層滑移相比較，發(fā)現(xiàn)震前斷層的滑移趨勢與同震斷層滑移一致，均為逆沖兼右旋的擠壓運動。比較震前3期逆沖方向的滑移量，發(fā)現(xiàn)逆沖滑移有加速的現(xiàn)象。并根據(jù)震前和同震的斷層滑動量估算了汶川地震復(fù)發(fā)周期。

GPS 斷層滑移 震前 同震 汶川地震

引言

地震通常都伴隨斷層的滑移，長期構(gòu)造運動產(chǎn)生地殼變形，當(dāng)變形超過斷層的承受臨界值時，斷層就會發(fā)生滑動破壞。斷層滑移是理解震源機制和發(fā)震機理的重要參數(shù)（柳暢等，2014）。目前對斷層滑移的研究大多集中于同震滑移，由于地表斷層破裂數(shù)據(jù)不能反映斷層的真實滑移，李志才等（2009）和王閱兵等（2012）根據(jù)GPS觀測的地表同震變形數(shù)據(jù)，通過對斷層同震滑移的反演計算，分析地震的發(fā)震類型和破壞特征，并據(jù)此理解震源機制。震前斷層滑動是與地震的孕育和發(fā)生直接關(guān)聯(lián)的物理現(xiàn)象。所以，震前斷層滑移與同震參數(shù)的比較分析是理解和進(jìn)行地震預(yù)測的一個具有明顯物理背景的分析方法。

基于此，本文以汶川地震為背景，利用中國地殼運動觀測中心提供的1999年至2001年、2001年至2004年、2004年至2007年共3期的153個固定站點GPS測量數(shù)據(jù)（水平精度2mm，垂直精度6mm）分析震前地表變形的變化趨勢，采用遺傳算法反演震前3期斷層滑移參數(shù)，并將震前斷層滑移參數(shù)與同震測量的38個GPS站點數(shù)據(jù)反演得到的同震斷層參數(shù)進(jìn)行比較，分析震前斷層滑移與同震斷層滑移的關(guān)聯(lián)。

1 斷層位置與GPS測站布置

汶川地震震后地質(zhì)考查表明（焦青等，2008；劉靜等，2008；李勇等，2008），此次地震由中央斷裂帶、山前江油-灌縣斷裂帶和山后茂汶-汶川斷裂帶等多條斷層同時參與地震斷層破裂過程，但中央斷裂帶為主要斷裂帶，其走滑和傾滑位移較大，山前和山后斷裂帶傾滑運動很小，可以忽略，僅考慮其走滑運動，因此本文將中央斷裂帶作為模型研究對象。斷層位置與GPS測點分布如圖1所示，根據(jù)斷層位置，將其擬合為一條直線，將斷層視為平面，且斷層面上的位移是均勻分布的。資料顯示（陳運泰，2008），斷層破裂長度在300km左右，通過地質(zhì)考察等方式（李海兵等，2008；張培震等，2008）得到汶川地震主斷裂帶斷層近地面傾角為60°—70°，斷層在20km處斷層隨深度向下逐漸變緩。根據(jù)這些背景資料，本文選取斷層長度=300km，傾角=60°，斷層深度=20km作為模型形狀參數(shù)。

2 震前3期斷層滑移參數(shù)

2.1 震前斷層參數(shù)計算方法

斷層滑移參數(shù)的計算通常采用反演算法（徐果明，2003）。采用圖2所示的斷層位錯模型，圖中豎直方向坐標(biāo)軸的原點=0定義在地表平面，為震源深度。斷層長、寬分別為和，傾角為。軸平行于斷層走向。1、2、3分別代表斷層上盤相對于下盤的走滑、傾滑和張引位錯分量。

為了提高斷層位移反演的精度，本文采用所有GPS測量點的3個方向測量值和模型值差的絕對值作為目標(biāo)函數(shù)值，即：

在整個算法中的關(guān)鍵步驟是基于位錯模型（Steketee，1958；Okada，1985；1992）計算斷層滑動參數(shù)1、2、3對應(yīng)的地表變形u0、u0、u0。根據(jù)點源位錯模型，斷層內(nèi)一點位移對地表產(chǎn)生的影響可以表示為：

由點源位錯計算公式（2）沿斷層面積分，可得有限矩形位錯公式。

2.2 模型驗證

下面，通過一個模擬斷層來校驗本文算法。已知斷層滑移等參數(shù)列于表1，其中為斷層長度，為斷層寬度，為斷層深度，為斷層傾角，1、2、3為斷層3個方向的滑移。據(jù)此計算地表位移分量，選取其中5組數(shù)據(jù)（表2）進(jìn)行分析。然后，基于這5組數(shù)據(jù)，利用遺傳算法反演斷層滑動參數(shù)，結(jié)果列于表3。計算結(jié)果與給出的斷層參數(shù)一致。為更進(jìn)一步說明計算效果，圖3示出了反演計算中目標(biāo)函數(shù)演化過程。可以看出，目標(biāo)函數(shù)收斂迅速，較好地趨近0且并未有波動趨勢，表明本文算法的穩(wěn)定性和較好的收斂效果。

表1 模擬斷層參數(shù)

表2 斷層正演結(jié)果

表3 斷層反演結(jié)果

2.3 震前三期斷層滑移特征

將GPS測量的震前3期地表位移數(shù)據(jù)進(jìn)行高斯換算和坐標(biāo)旋轉(zhuǎn)，把大地坐標(biāo)系轉(zhuǎn)換到以斷層下緣線為水平軸的平面坐標(biāo)系中。利用變換后的地表位移以及上述計算模型反演龍門山斷層1999年到2007年震前3期滑動位移的累積量和斷層同震滑移量。震前3期的時間（以年為單位）為：第一期1999.771—2001.678，時間間隔1.907年；第二期為2001.678—2004.583，時間間隔2.905年，第三期為2004.583—2007.579，時間間隔2.996年。

從震前3期的目標(biāo)函數(shù)演化過程（圖4）可以看出，3期結(jié)果均收斂較好，其中3期的目標(biāo)函數(shù)值分別為0.0231mm、0.0261mm、0.0283mm，平均每個測點每個方向的誤差小于1×10-4mm，同震目標(biāo)函數(shù)值為3.93×10-2m，每個測點每個方向的誤差小于10-3m，震前3期和同震反演誤差均小于GPS測量數(shù)據(jù)的誤差，故反演結(jié)果可靠。

同震斷層附近地表位移如圖5所示，同震反演斷層走滑位移為-3.029m，傾向位移為4.725m。震后地質(zhì)考察結(jié)果顯示（吳珍漢等，2008；劉健等，2012），中央斷裂帶中，映秀至北川斷裂帶平均逆沖位移為6.4m，北川至青川斷裂帶逆沖為5m，映秀至北川破裂帶垂直位移量從0.2m到11m不等，伴隨的右旋平均水平位移量為2—3m。所以，同震斷層反演結(jié)果與地質(zhì)觀測結(jié)果相符合。陳運泰等（2008）通過地震波反演得出，斷層平均逆沖位錯為5m，斷層平均走滑位錯為2m，與本文的結(jié)果基本一致。

另一方面，根據(jù)斷層滑移參數(shù)可以計算地震矩：

計算得矩震級W=7.985。所以，本文反演的結(jié)果與根據(jù)地震波計算的矩震級為W=7.9（張勇等，2008）基本相當(dāng)。

為了與同震滑移數(shù)據(jù)進(jìn)行對比，了解震前滑移與地震的關(guān)聯(lián)，表4給出了震前3期的斷層滑移的反演結(jié)果?？梢钥闯?，震前3期總的平均走滑速率為-5.39mm/a，傾向速率為2.66mm/a，斷層有逆沖和右旋走滑的趨勢。由于走滑位移會受到斷層兩端的約束，而傾向位移不受地表約束更易發(fā)生滑動，所以從震前斷層滑移參數(shù)分析可知，斷層的傾向滑移以逆沖為主。這與震后觀測結(jié)果及上述同震反演結(jié)果（逆沖和右旋運動）一致。這個結(jié)果可能預(yù)示著，可以通過震前斷層的運動趨勢來推測地震發(fā)生時斷層的破壞形式。

進(jìn)一步地，斷層震前各期的反演結(jié)果（表4）表明，震前第一、二、三期的年均走向滑移分別為5.76mm/a、5.44mm/a、4.97mm/a，逐漸減小。而傾向的年均滑移分別為2.29mm/a、2.66mm/a和3.02mm/a，顯示上升趨勢。這表明斷層的矢量運動方向從以走滑為主的右旋運動逐漸偏轉(zhuǎn)到以傾向滑移為主的逆沖運動。該趨勢預(yù)示著斷層在未來地震時可能發(fā)生以逆沖為主的破壞，這與汶川地震時斷層發(fā)生的巨大逆沖滑移相一致。

表4 龍門山斷層震前3期反演結(jié)果

龍門山斷層震前和同震的滑移參數(shù)可以反映地震的能量累計和釋放量，進(jìn)而能夠估算汶川地震的復(fù)發(fā)周期。龍門山斷裂帶晚第四紀(jì)以來地貌學(xué)的研究結(jié)果顯示（馬保起等，2005），茂汶-汶川斷裂、北川-映秀斷裂和江油-灌縣斷裂晚第四紀(jì)逆沖滑動速率的總和約為1mm/a左右。龍門山斷裂帶歷史上沒有超過7級的地震，6級以上地震只有3次，由于地震震級相差一級釋放的能量相差很大，所以忽略復(fù)發(fā)周期內(nèi)地震釋放的能量，由于臨近地震，后2期的逆沖滑動有加速趨勢，故將第一期逆沖滑動速率2.29mm/a作為長期積累量，每年平均釋放量為1mm/a，用積累量與釋放量之差計算凈積累量為1.29mm/a，利用斷層平均滑動位移4.725m與斷層的滑動速率的凈積累量之比計算得到地震復(fù)發(fā)周期為3700年，地質(zhì)學(xué)估計的地震復(fù)發(fā)周期為2000—10000年（謝富仁等，2008；李玉江等，2012）。

3 結(jié)論

（1）本文計算了震前3期斷層滑移累積量，斷層平均走滑速率為-5.39mm/a，傾滑向速率為2.66mm/a，與同震相比，斷層運動形式一致，都為逆沖兼右旋的滑動。

（2）比較震前3期的滑移速率，從第一、二期到第二、三期逆沖滑移速度增加，說明逆沖方向有加速的趨勢，斷層滑移加速方向和地震發(fā)生時斷層的滑動形式一致。

（3）根據(jù)震前、同震斷層滑移和地質(zhì)資料，利用斷層同震逆沖滑移與震前斷層滑動速率的凈積累量推算了汶川地震的復(fù)發(fā)周期，結(jié)果表明汶川地震復(fù)發(fā)周期為3700年。

陳運泰，2008．汶川特大地震的震級和斷層長度．科技導(dǎo)報，26（10）：26—27．

陳運泰，許力生，張勇等，2008．2008年5月12日汶川特大地震震源特性分析報告．北京：中國地震局地球物理所．

焦青，楊選輝，許麗卿，2008．汶川8.0級地震前后龍門山斷裂活動特征淺析．大地測量與地球動力學(xué)，28（4）：7—11，31．

雷英杰，2014．MATLAB遺傳算法工具箱及應(yīng)用．2版．西安：西安電子科技大學(xué)出版社．

李海兵，付小方，VAN DER WOERD J.等，2008．汶川地震（S8.0）地表破裂及其同震右旋斜向逆沖作用．地質(zhì)學(xué)報，82（12）：1623—1643．

李勇，周榮軍，Densmore A. L.，2008．映秀-北川斷裂的地表破裂與變形特征．地質(zhì)學(xué)報，82（12）：1688—1706．

李玉江，陳連旺，2012．汶川S8.0地震大震復(fù)發(fā)周期的研究進(jìn)展．地球物理學(xué)進(jìn)展，27（2）：455—463．

李志才，張鵬，金雙根，蔣志浩，溫?fù)P茂，2009．基于GPS觀測數(shù)據(jù)的汶川地震斷層形變反演分析．測繪學(xué)報，38（2）：108—113．

劉健，熊探宇，趙越，張永雙，陳群策，2012．龍門山活動斷裂帶運動學(xué)特征及其構(gòu)造意義．吉林大學(xué)學(xué)報（地球科學(xué)版），42（S2）：320—330．

劉靜，張智慧，文力等，2008．汶川8級大地震同震破裂的特殊性及構(gòu)造意義——多條平行斷裂時活動的反序型逆沖地震事件．地質(zhì)學(xué)報，82（12）：1707—1722．

柳暢，石耀霖，朱伯靖，程惠紅，楊小林，2014．地殼流變結(jié)構(gòu)控制作用下的龍門山斷裂帶地震發(fā)生機理．地球物理學(xué)報，57（2）：404—418．

馬保起，蘇剛，侯治華等，2005．利用岷江階地的變形估算龍門山斷裂帶中段晚第四紀(jì)滑動速率．地震地質(zhì)，27（2）：234—242．

王閱兵，金紅林，付廣裕，2012．Yabuki&Matsu'ura方法在汶川W7.9地震反演中的應(yīng)用．地震，32（2）：121—128．

吳珍漢，張作辰，2008．四川汶川S8.0級地震的地表變形與同震位移．地質(zhì)通報，27（12）：2067—2075．

謝富仁，張永慶，張效亮，2008．汶川S8.0級地震發(fā)震構(gòu)造大震復(fù)發(fā)間隔估算．震災(zāi)防御技術(shù)，3（4）：337—344．

徐果明，2003．反演理論及其應(yīng)用．北京：地震出版社．

張培震，徐錫偉，聞學(xué)澤等，2008．2008年汶川8.0級地震發(fā)震斷裂的滑動速率、復(fù)發(fā)周期和構(gòu)造成因．地球物理學(xué)報，51（4）：1066—1073．

張勇，馮萬鵬，許力生，周成虎，陳運泰，2008．2008年汶川大地震的時空破裂過程．中國科學(xué) D輯：地球科學(xué)，38（10）：1186—1194．

周明，1999．遺傳算法原理及應(yīng)用．北京：國防工業(yè)出版社．

Okada Y., 1985. Surface deformation due to shear and tensile faults in a halfspace. Bulletin of the Seismological Society of America, 75(4):1135—1154.

Okada Y., 1992. Internal deformation due to shear and tensile faults in a half-space. Bulletin of the Seismological Society of America, 92(2):1018—1040.

Steketee J. A., 1958. On volterra's dislocations in a semi-infinite elastic medium. Canadian Journal of Physics, 36(2):192—205.

Comparsion of Pre-earthquake and Coseismic Fault Slip of the 2008 Wenchuan Earthquake

Liang Xiangzhou and Hao Shengwang

(School of Civil Engineering and Mechanics, Yanshan University, Qinhuangdao 066004, China)

In order to get the acknowledge about the pre-earthquake fault slip and its relationship with coseismic slip of the 2008 Wenchuan Earthquake, this paper inverts the GPS measurements from 1999 to 2007 and the coseismic measurementsto infer the pre-earthquake and coseismic slip characteristics of Longmenshan fault. The earthquake dislocation theory and the genetic algorithm are used in the inversion. It is shown that the mean strike-slip from 1999 to 2007 is about 5.39 mm, and the dip slip is almost 2.66 mm. This pre-earthquake slip tendency is consistent with the coseismic reverse and right-lateral slip. The results of the dip slips occurred in 1999 to 2001, 2001 to 2004 and 2004 to 2007 before the 2008 Wenchuan Earthquake indicatethat the dip slip is accelerating with time. This acceleration slip may be a precursor of earthquake. Finally, based on the inversion results of the pre-earthquanke and coseismic fault slip, we estimate the recurrence period of Wenchuan earthquake.

GPS；Fault slip；Pre-earthquake；Coseisimic；Wenchuan Earthquake

國家自然科學(xué)基金（11672258）、河北省自然科學(xué)基金（D2015203398）和973計劃（2013CB834100）

2016-04-11

梁向洲，男，生于1991年。碩士研究生。主要從事地震斷層滑移與地震前兆的研究。 E-mai：liangxiangzhou@foxmail.com

梁向洲，郝圣旺，2017．汶川地震震前與同震斷層滑移的比較分析．震災(zāi)防御技術(shù)，12（1）：107—114．doi：10.11899/zzfy20170111