潘 睿，盛書中，萬永革

（防災(zāi)科技學(xué)院，河北 三河 065201）

0 引言

據(jù)中國地震臺網(wǎng)中心測定，北京時間2009年7月9日19時19分，云南省姚安縣發(fā)生了6.0級地震（震中25.6°N，101.1°E），震源深度約10km，本次地震造成重大人員傷亡和財產(chǎn)損失。本次地震前在其附近地區(qū)自2000年以來，已經(jīng)發(fā)生了4次中等強(qiáng)度地震（圖1）。然而，該區(qū)域的1∶20萬區(qū)域地質(zhì)資料和先前的地震地質(zhì)資料都未能顯示或發(fā)現(xiàn)任何活動斷裂構(gòu)造的形跡［1］，反映本次地震斷層可能為地下隱伏斷層。

假定可以用一個平面來模擬地震發(fā)震斷層面（即數(shù)學(xué)模型），即大多數(shù)余震會叢集發(fā)生在此主震斷層面上或附近，因此，根據(jù)余震空間分布就可以較精確地勾畫出主震斷層面的形狀和位置（即發(fā)震斷層的走向和傾角）。王鳴等［2］采用這種原則和方法，利用1989年大同—陽高地震余震精確定位參數(shù)，反演出該地震主震斷層面參數(shù)。萬永革等［3］采用模擬退火全局搜索和高斯—牛頓局部搜索相結(jié)合的反演方法，利用2002年4月1日至2006年5月31日發(fā)生在唐山地震破裂區(qū)的“余”震精確定位資料，求得了唐山地震、灤縣地震、寧河及盧龍斷裂帶的斷層面走向、傾角、位置及滑動角等參數(shù)。該方法在斷層面參數(shù)研究上得到了較為廣泛的應(yīng)用［4－6］，本研究將該方法應(yīng)用于2009年姚安地震發(fā)震斷層面參數(shù)的確定，為后續(xù)相關(guān)研究提供參考。

圖1 研究區(qū)域地震構(gòu)造簡圖［8］

1 數(shù)據(jù)的篩選及反演結(jié)果

1.1 數(shù)據(jù)情況

姚安地震后，云南省地震局快速地在余震區(qū)布設(shè)了流動數(shù)字臺網(wǎng)，該臺網(wǎng)和固定臺站記錄到了大量余震，王長在等［7］對這些余震資料進(jìn)行了雙差層析成像研究，給出了大量精準(zhǔn)定位的余震資料，為本文研究提供了基礎(chǔ)。本研究使用的747個余震的分布情況見圖2，震級分布和震源深度分布見圖3。

1.2 斷層面參數(shù)反演結(jié)果

圖2 余震震中分布圖

圖3 姚安地震基本要素分布圖

從震中分布圖2可見，余震震中分布具有較強(qiáng)的叢集性，雖有少數(shù)余震分布較為偏遠(yuǎn)，但是總體形態(tài)上并沒有削弱余震叢集性，所以本文研究依據(jù)余震叢集性原則選取圖2中四邊形內(nèi)的余震反演斷層面參數(shù)。反演結(jié)果見圖4，發(fā)震斷層面的走向為296.6°，傾角為86.6°；由圖4c和4d可見，余震大體上對稱地分布在發(fā)震斷層面兩側(cè)；根據(jù)萬永革等［3］所提的斷層頂點(diǎn)確定原則，確定斷層面的4個頂點(diǎn) 的 位 置 分 別 為：（25.53°N，101.10°E，2.42km），（25.54° N， 101.11°E， 18.47km），（25.59°N，101.00°E，18.47km）和（25.58°N，101.00°E，2.42 km）。

1.3 斷層反演結(jié)果分析

為分析我們反演所得結(jié)果的合理性，我們搜集了其他專家學(xué)者對本次姚安地震斷層面參數(shù)的研究結(jié)果（表1），秦雙龍等［8］利用P波初動資料解出地震斷層面參數(shù)，張建國等［1］結(jié)合流動數(shù)字地震觀測的資料和地震序列觀測記錄的雙差定位處理，并借鑒相關(guān)地震勘探的研究成果獲得地震斷層面參數(shù)。

圖4 本次地震斷層面參數(shù)結(jié)果

表1 2009年云南姚安MS6.0級地震斷層面參數(shù)表

由表1可見，本文所獲取的走向和USGS以及秦雙龍等［8］基本一致，與張建國等［1］存在7.85°的差異；本研究所得的傾角與其它研究結(jié)果存在10°左右的差異。為了判斷該差異是否在合理的范圍內(nèi)，我們搜集了其他14組地震斷層面參數(shù)研究結(jié)果（表2），并對這14組研究結(jié)果進(jìn)行統(tǒng)計分析，分析14組地震斷層面參數(shù)研究結(jié)果的差異范圍，從而為本文研究結(jié)果是否合理提供判斷依據(jù)。

表2 地震斷層面參數(shù)研究結(jié)果差異統(tǒng)計表

續(xù)表2

圖5 同一地震斷層走向（a）和傾角（b）研究結(jié)果最大差異統(tǒng)計圖

根據(jù)圖5可見，盡管這14次地震斷層面的走向和傾角研究結(jié)果的最大差異較為離散，但是總體上不影響其作為參照來考查本研究結(jié)果的合理性，可見我們給出的姚安地震主斷層面參數(shù)與其他研究者給出結(jié)果的差異是可以接受的，即相對于這14次地震斷層面參數(shù)研究結(jié)果的差異，本研究結(jié)果與其他研究結(jié)果的差異在可接受的差異范圍內(nèi)，故本研究所得的姚安地震斷層面參數(shù)是可信的。

2 驗證斷層面參數(shù)結(jié)果穩(wěn)定性

圖6 斷層走向（a）和斷層傾角（b）隨時間變化曲線圖

當(dāng)中等強(qiáng)度及以上地震發(fā)生后，需要及時地對地震發(fā)震斷層面參數(shù)進(jìn)行研究，以便及時判斷災(zāi)情，為震后應(yīng)急救援等提供依據(jù)?？梢姡鸷罂焖佾@取斷層面參數(shù)具有重要意義，那么本方法是否能在震后快速給出可靠的發(fā)震斷層面參數(shù)呢？刁桂苓等［24］對蘆山地震斷層的研究表明，地震序列早期的地震目錄可以確定斷層面參數(shù)，但僅應(yīng)用于限定極震區(qū)的范圍。針對該問題并結(jié)合本次地震，我們將探討本次地震后利用精定位余震反演斷層面方法何時能反演出可靠的斷層面參數(shù)，即對反演結(jié)果趨于穩(wěn)定所需的時間。

依據(jù)精定位余震的發(fā)震時刻，本研究按余震發(fā)生的時間順序來選擇余震，選擇距離首個余震發(fā)生后24h內(nèi)的余震作為震后第1天的數(shù)據(jù)，選擇震后48h內(nèi)的數(shù)據(jù)為第2天的數(shù)據(jù)，依次類推，選取了震后12 d（288h）內(nèi)的余震數(shù)據(jù)，用這12d數(shù)據(jù)對斷層面參數(shù)分別進(jìn)行反演，反演結(jié)果見圖6。

由圖6可見，從第2天起，斷層走向和傾角反演結(jié)果趨于穩(wěn)定狀態(tài)；第1天所發(fā)生的余震數(shù)量占全部余震數(shù)量的14.1%，數(shù)據(jù)量略少，所以導(dǎo)致其反演所得地震斷層面參數(shù)不穩(wěn)定；第2天累積的余震數(shù)量占全部余震數(shù)量的40%左右，所得發(fā)震斷層面參數(shù)開始趨于穩(wěn)定。在本震例中我們使用震后2d內(nèi)的精定位余震就可以獲取穩(wěn)定地震主斷層面參數(shù)，這說明使用本文所用反演方法，利用主震后短時間內(nèi)發(fā)生的余震精定位結(jié)果可以快速地反演出主斷層面參數(shù)，反映該反演方法可能為震后應(yīng)急提供快速反演斷層面參數(shù)的技術(shù)支持，為今后其它地震發(fā)震斷層面參數(shù)的快速獲取提供新的思路。

3 討論與結(jié)論

根據(jù)秦雙龍等［8］以及張建國等［1］的對該地震震源機(jī)制解研究，他們均認(rèn)為本次姚安地震主斷裂走向是NWW 向，傾角陡立。精定位余震震中分布明顯呈NWW 向密集條帶狀（圖2），反演出發(fā)震斷層面與余震關(guān)系情況（圖4），我們可以推斷本次地震發(fā)震構(gòu)造為馬尾箐斷裂，這也與秦雙龍等［8］的研究結(jié)果相符合。本研究反演出的發(fā)震斷層傾角約86°，與其得出斷層傾角陡立相符合。

本文將利用余震分布確定發(fā)震斷層面的方法應(yīng)用于2009年姚安地震，得到以下結(jié)果：1）2009年姚安地震發(fā)震斷層面的走向和傾角分別為296.6°和86.6°，并推斷本次地震發(fā)震構(gòu)造為所在區(qū)域馬尾箐斷裂；2）本震例中使用震后2d內(nèi)余震資料便可以獲得可靠的發(fā)震斷層面參數(shù)，這說明利用大地震后短時間內(nèi)所發(fā)生的余震資料，可以較為快速地確定地震斷層面參數(shù)，為今后其它地震發(fā)震斷層面參數(shù)的快速獲取提供新的思路。本研究是在余震精準(zhǔn)定位的基礎(chǔ)上進(jìn)行的，余震若為速報結(jié)果，或是其它震例，何時能獲得穩(wěn)定的斷層面參數(shù)？還有待于進(jìn)一步研究。

［1］ 張建國，劉麗芳，李西，等.姚安、大姚中強(qiáng)震區(qū)地震構(gòu)造初析［J］.地震地質(zhì)，2009，31（3）：536－542.

［2］ 王鳴，王培德.1989年10月18日大同—陽高地震的震源機(jī)制和發(fā)震構(gòu)造［J］.地震學(xué)報，1992，14（4）：407－415.

［3］ 萬永革，沈正康，刁桂岺，等.利用小震分布和區(qū)域應(yīng)力場確定大震斷層面參數(shù)方法及其在唐山地震序列的應(yīng)用［J］.地球物理學(xué)報，2008，51（3）：793－804.

［4］ 楊超群，孟凡順，萬永革.采用精確定位小震資料反演伽師地震斷層面［J］.地球物理學(xué)進(jìn)展，2013，28（6）：2865－2871.

［5］ 盛書中，萬永革，王未來，等.2010年玉樹MS7.1級地震發(fā)震斷層面參數(shù)的確定［J］.地球物理學(xué)進(jìn)展，2014，29（4）：1555－1562.

［6］ Wang X S，F(xiàn)eng X D，Xu X W，et al.Fault plane parameters of Shanhe－Pinggu M 8earthquake in 1679determined using present－day small earthquakes［J］.Earthq Sci，2014，27：607－614.

［7］ 王長在，吳建平，房立華，等.2009年姚安地震序列及震源區(qū)三維P波速度結(jié)構(gòu)研究［J］.地震學(xué)報，2011，33（2）：123－133.

［8］ 秦雙龍，張建國，廖麗霞，等.2009年姚安6.0級地震震源機(jī)制與發(fā)震構(gòu)造的分析研究［J］.內(nèi)陸地震，2012，26（1）：52－61.

［9］ 龍思勝.2001年四川雅江6.0級地震序列的破裂特征及發(fā)震構(gòu)造［J］.中國地震，2004，20（1）：1－11.

［10］ 高國英，聶曉紅，夏愛國，等.2003年伽師6.8級地震序列特征和震源機(jī)制的初步研究［J］.中國地震，2004，20（2）：179－186.

［11］ 王衛(wèi)民，李麗，趙連鋒，等.2003年2月24日新疆伽師MS6.5級地震震源破裂過程研究［J］.地球物理學(xué)報，2005，48（2）：343－351.

［12］ 趙翠萍，陳章立，鄭斯華，等.2003年9月27日中、俄、蒙邊界MS7.9地震震源機(jī)制及破裂過程［J］.2005，27（3）：237－249.

［13］ 郭志，高星，王衛(wèi)民，等.2006年12月6日臺灣南部濱海MS7.2級地震破裂過程研究［J］.地球物理學(xué)報，2008，51（4）：1103－1113.

［14］ 郭祥云，陳學(xué)忠，李艷娥，等.2008年5月12日四川汶川8.0級地震與部分余震的震源機(jī)制解［J］.地震，2010，30（1）：50－60.

［15］ 鄭勇，馬宏生，呂堅，等.汶川地震強(qiáng)余震MS5.6的震源機(jī)制解及其與發(fā)震構(gòu)造的關(guān)系［J］.中國科學(xué)：地球科學(xué)，2009，34（4）：413－426.

［16］ 蘇金蓉，郭志.西南天山—帕米爾2008年10月5日MS6.7級地震震源機(jī)制研究［J］.地球物理學(xué)報，2013，56（2）：504－512.

［17］ 張輝，王熠熙.2012年5月3日金塔—阿拉善盟5.4級地震震源機(jī)制解［J］.西北地震學(xué)報，2012，34（2）：205－206.

［18］ 謝祖軍，金筆凱，鄭勇，等.近遠(yuǎn)震波形反演2013年蘆山地震震源參數(shù)［J］.中國科學(xué)：地球科學(xué)，2013，43（6）：1010－1019.

［19］ 呂堅，王曉山，蘇金蓉，等.蘆山7.0級地震序列的震源位置與震源機(jī)制解特征［J］.地球物理學(xué)報，2013，56（5）：1753－1763.

［20］ 羅鈞，趙翠萍，周連慶.2013年8月香格里拉—德榮MS5.9地震序列震源機(jī)制與應(yīng)力場特征［J］.地球物理學(xué)報，2015，58（2）：424－435.

［21］ 趙凌云，張輝，陳俊華，等.基于CAP方法得2013湖北巴東5.1級地震震源機(jī)制與發(fā)震構(gòu)造背景研究［J］.地震工程學(xué)報，2014，36（4）：1014－1018.

［22］ 趙旭，黃志斌，房立華，等.2014年云南盈江6.1地震震源機(jī)制研究［J］.中國地震，2014，30（3）：462－473.

［23］ 劉麗芳，徐甫坤.利用精確定位余震資料確定2014年云南魯?shù)?.5級地震的斷層面參數(shù)［J］.地震研究，2014，37（4）：489－494.

［24］ 刁桂苓，馮向東，趙翠萍，等.蘆山地震的震源斷層及在應(yīng)急救援中應(yīng)用的可能性［J］.國際地震動態(tài)，2013，11：32－33.