四川木里小震群活動特征研究1

蔡一川 戴仕貴 魏婭玲

（四川省地震局，成都 610041）

四川木里小震群活動特征研究1

蔡一川 戴仕貴 魏婭玲

（四川省地震局，成都 610041）

本文通過對四川木里地區地震震源位置的重新確定，反演了地震震源機制解，結果顯示：①重定位后的小震群震中分布成帶狀，地震震源深度分布在0—12km范圍內；②經過對震群空間分布進行仔細分析，認為其發震構造是小金河斷裂西側的一條NWW向分支斷裂；③3次4.0級以上地震震源參數存在明顯差異，淺源逆沖事件表現出受垂直方向應力（重力）作用的特征，走滑事件表現出與區域應力構造活動有關。

震群活動 地震條帶 雙差定位方法 震源機制

蔡一川，戴仕貴，魏婭玲，2015.四川木里小震群活動特征研究.震災防御技術，10（4）：933—940.doi：10.11899/zzfy20150411

引言

地震震群活動是地殼變動的結果，是地殼運動增強的客觀反映，當大范圍構造應力出現增強時，首先會在介質相對較薄弱的部位發生顯著地震或以震群活動的方式表現出來。本文中的“小震群”意指在較短時間內局部區域內密集發生的一組中小地震（MS≤5.0）。正在活動中的震群，究竟是“大震”的前兆震群，還是并不伴有后續大震發生一般性的震群活動，對其作出快速、準確判別，對減輕地震災害和維護社會穩定具有重要的現實意義。在我國，基于小震群活動開展的大震預測研究一直是眾多研究人員關注的重點（趙根模，1990；柯龍生等，1977；朱麗霞，1999；啜永清等，2000；Lund等，2002；劉正榮，2004；王俊國等，2005；崔子健等，2011）。

四川木里縣地處青藏高原南緣，是橫斷山脈在四川境內表現最典型的地帶，地質結構復雜；其大地構造位置屬于松潘-甘孜地槽區南緣的兩個I級構造單元揚子準地臺與松潘-甘孜地槽的槽臺結合部；小金河斷裂帶為槽臺結合的接觸性斷裂，位于木里縣城南14km處，此斷裂由一系列東西向展布、北傾、中部向南弧形突出的斷裂組成。2013年7月至12月期間，四川區域臺網地震觀測資料揭示，在四川木里地區發生了大量微震和有感小震群活動，11月22日發生的3次ML≥4.0地震造成木里縣部分房屋開裂，給當地居民的生活帶來了不小的影響，因此很有必要對該震群地震活動特征進行分析研究。定位結果表明，這些地震以震群形式集中發生在四川木里縣項腳鄉，得到的可定位地震竟達64 538次，最大震級為ML4.5，震中集中分布在幾公里范圍內。目前，由于研究者們對此震群成因尚存有不同認識，因此本文采用雙差定位方法對此震群進行重新定位，用CAP方法對ML≥4.0地震求解震源機制解，對震群地震頻次與錦屏一級水庫水位關系進行分析，試圖對此小震群活動特征進行一些初步的探索和分析。

1 資料與方法

1.1 臺站和資料的選取

四川木里震群發生前，震中集中區域附近250km范圍內正常運行的固定臺站22個（圖1，紅色三角形）、水庫臺站37個（圖1，綠色五角星）。在臺站觀測儀器的配置上，22個固定臺站全部配備寬頻帶記錄儀（四川臺17個、云南臺5個），37個水庫臺站全部配備短周期記錄儀。四川、云南區域臺網固定臺站與四川水庫企業臺站并網監測，構成的地震監測網對木里震群區形成了很好的方位覆蓋，保證了地震定位結果的可靠性。其中，記錄Sg與Pg到時差小于4s的臺站有5個，近臺方位分布較好，在地震事件定位中發揮了重要作用。

圖1 重定位臺站分布Fig.1 The distribution of stations for the earthquake relocation

2013年7月至12月期間，并網觀測到的可定位地震64538次，其中ML4.0—4.9地震4次，ML3.0—3.9地震38次，ML2.0—2.9地震489次，ML1.0—1.9地震5952次，為研究木里震群活動特征提供了資料保證。在進行初始編目定位時，采用了LOC3D、Hypo2000和單純形定位方法，由于不同定位方法的速度模型和參數存在差別，因此定位結果可能出現系統偏差。在研究中為了獲得較為可靠的初始震源位置，選擇ML≥1.0、記錄臺站數≥5的地震事件，統一采用Hypo2000定位方法（Wang等，2007）和林向東等（2013）的速度模型進行重定位。林向東等（2013）的模型是根據川西地區人工地震測深剖面和波形擬合結果，并做適當調整后建立的速度模型（見圖2）。

圖2 研究地區地殼速度模型Fig. 2 The crustal velocity model used in the study region

1.2 方法基本原理及參數設置

在小震定位方面采用了雙差定位方法，該方法的基本原理是：如果兩個地震震源之間的距離小于事件到臺站的距離，那么震源區和這個臺站之間的整個射線路徑幾乎相同，則可以通過地震事件兩兩組對的方法在一定程度上消除地殼速度結構橫向不均勻性帶來的定位誤差。近年來，雙差定位法在我國得到了廣泛應用（李永紅等，2011；馮紅武等，2013；張廣偉等，2013）。本研究中采用適合木里地區的速度模型開展雙差定位。由于S波到時資料拾取精度相對P波要低些，因此在定位過程中將P波權重設為0.5—1.0，而S波權重設為0.1—0.5。為了獲得較高精度的定位結果，設定事件對的距離閥值為800m，事件對與臺站之間的最大距離設為200km，最小連接數和最小觀測數為8。這些設置在本例的實際應用中很好地排除了干擾，既保留了P波震相相對較高的精度，又保證了足夠多的數據量。

在震源機制的研究中，國內外學者經過系統研究積累了大量寶貴的經驗，已經取得了大量重要成果（呂堅等，2011；易桂喜等，2012；趙凌云等，2010）。其中，Zhao等（1994）和Zhu等（1996）提出的CAP（Cut And Paste）方法獲取震源機制解，在國內外諸多研究成果中得到了充分證明。其優點體現在反演震源機制解時使用的是近震寬頻帶波形資料，彌補了利用P波初動求解震源機制解受臺站數量限制的缺陷。

2 震群重定位結果分析

為了了解木里小震群的活動特征，本文盡可能多地收集研究區域內所有地震的各類震相數據，首先采用Hypo2000定位方法進行初始定位，然后采用雙差定位方法重新定位。共獲得6094個重定位結果，定位成功率為94.5%，E-W、N-S和U-D三個方向大致的定位誤差分別為0.21、0.2和0.38km。

圖3 木里小震群重定位前、后震中分布Fig.3 The epicenter distribution before and after relocation of the earthquake swarm in the Muli area

圖3（a）是初始定位震中分布圖，結果顯示地震震中集中在小金斷裂帶西側的木里項腳一帶，分布成團狀，沒有明顯的幾何特征，僅有少數微震彌散在其四周。重定位后，ML≥1.0地震震中分布顯示出在平面上高度集中，在西北-東南向呈條帶分布，對密集區進行放大，獲得重定位震中分布如圖3（b）所示，結果發現小震群分布呈“劍狀”幾何特征，“劍柄”位于東南端，“劍尖”指向西北。通過對震群空間分布判定，木里震群并不發生在小金主斷裂帶上，而是發生在小金斷裂西側的NWW向分支斷裂上。依照圖3（b）對小震群序列按NW和NE方向進行了剖面劃分，并將地震震源深度投影到AA'和BB'兩個剖面上（見圖4）。圖4（a）是地震條帶沿直線AA'的震源擴展剖面圖，顯示木里小震群總體震源深度較淺，在地震條帶AA'的東南段，地震震源深度分布范圍為0—12km，西北段地震震源深度集中分布在7—10km范圍內；地震條帶總長約8km，空間幾何圖形就像一個橫著的大寫字母“L”。圖4（b）是地震條帶沿直線BB'的震源擴展剖面圖，震源深度沿BB'剖面總體展布傾角較徒，由最小二乘擬合得到的平均展布傾角約86°，略傾向北，震源深度呈南淺北深的特征。

圖4 木里小震群重定位后深度剖面圖Fig.4 The depth profile map after relocation of the earthquake swarm in the Muli area

3 較大地震震源機制反演結果分析

首先，選取四川、云南區域臺網震中距小于250km的寬頻帶儀器數據資料，共有21個臺（圖1中的紅色三角形），在0.05—20Hz的頻率范圍內這些臺站儀器的幅頻特性曲線都是平直的，可獲得較好的臺站儀器記錄。然后，再從這些臺站中篩選出信噪比較高、P波初動明顯、臺站方位分布良好的臺站記錄參與反演計算。在研究時間段內，木里地區共發生ML≥4.0地震4次（見表1），表1中的事件D屬于疊加地震，故僅對其它3次地震進行了反演，且獲得了較可靠結果。反演結果顯示，3次ML≥4.0地震震源深度（表2）與精定位結果給出的震源深度值很接近（表1）。

表1 四川木里ML≥4.0地震重定位結果Table 1 The relocation results of ML≥4.0 earthquakes in the Muli area

表2 四川木里ML≥4.0地震震源機制解反演結果Table 2 The inversion results of ML≥4.0 earthquake in the Muli area

續表

采用CAP方法反演獲得了四川木里震群中3次ML≥4.0地震震源機制解，圖5給出了3次4級以上地震震源機制反演的深度擬合誤差分布圖，其震源機制解參數列于表2。結果顯示，2013年10月22日18時事件A的震源機制反演結果在2—8km震源深度范圍內比較穩定，震源深度在4km處獲得最小深度擬合誤差，這表明該地震屬淺源地震，錯動類型為逆沖，震源機制P軸方位近NS向，其傾角近水平（1°）；2013年11月22日03時事件B的震源機制反演結果在20km深度范圍內比較一致，震源深度在7km處獲得最小深度擬合誤差，錯動類型表現為走滑，其深度值較事件A的深度值增加了3km，震源機制P軸方位為N50°W向，傾角近水平（1°）；2013年11月22日07時事件C的震源機制反演結果在20km深度范圍內比較穩定，震源深度在8km處獲得最小深度擬合誤差，錯動類型也為走滑，其深度值較事件A的深度值增加了4km，震源機制P軸方位則表現為近N25°W向，傾角為6°。綜上所述，從地震震源機制解節面、力軸參數和錯動類型來看，這3次4級以上地震震源參數存在明顯差異，淺源逆沖事件A受垂直方向應力作用明顯，走滑事件B和C與區域構造活動有關，說明該區域發生的小震群活動具有一定的復雜性。

圖5 木里地區ML≥4.0地震震源機制隨不同震源深度變化圖Fig.5 The variation of the focal mechanism solution with different depth of ML≥4.0 earthquakes in the Muli area

4 震群活動與錦屏一級水庫蓄水的關系

由北向南穿行的雅礱江在四川木里縣與鹽源縣交界處的白碉一帶，突然來了一個近于180°的大轉折，向北繞行約150km，形成一個“幾”字形的錦屏大河灣后，然后才繼續向南流去。木里小震群位于錦屏一級水庫的上游，與庫大壩之間的距離約為50km，屬于錦屏一級水庫庫區范圍（見圖1）。2012年10月8日，錦屏一級水電站庫區開始下閘蓄水，到2012年底水位從海拔（下同）1652m升至1713m，完成第一階段蓄水，水位上升61m；2013年6月中旬開始第二階段蓄水，7月18日水位升至1800m的水庫死水位，完成第二階段蓄水，在第一階段蓄水的基礎上，水位上升87m，蓄水總量約24.2億立方米；2013年9月15日開始第三階段蓄水，至10月15日水位上升至1840m，水位再次上升40m（圖6a）。

圖6（b）是木里震群地震頻次5日滑動均值，從圖中可以看出，研究區域水庫地震類型屬于直接響應類型，在第一階段蓄水期間，研究區域地震頻次沒有發生變化；在第二階段蓄水期間，庫水位上升至一個新的高度，研究區地震頻次快速增加，由平靜期5次/日，迅速升高至400次/日；在第三階段蓄水期間，庫水位上升至最高1840m，地震頻次反而出現了階段性下降，較大地震的發生略滯后于最高庫水位，在庫水最高位后的6—40天內，發生了4次4.0級以上地震，地震頻次在11月22日達到了最高值1528次/日，隨后又逐漸減少。圖5顯示出地震頻次與庫水位上升有明顯關聯特征，從而可以判定庫水位加載過程所導致的裂隙或斷層面上的應力變化，成為庫區震群活動的直接影響因素。

圖6 錦屏一級水庫水位與地震頻次的關系Fig.6 The relationship between water level of Jinping reservoir and frequency of the earthquake swarms

5 小結與討論

通過對四川木里地區小震群震源位置的重新確定、反演較大地震震源機制解、小震群活動頻次與錦屏一級水庫庫水位上升關系的研究分析，得出以下結論：

（1）木里小震群地震頻次與錦屏一級水庫庫水位的上升有明顯關聯特征，庫水位加載過程所導致的斷層面上的應力變化是小震群活動的直接影響因素。

（2）該小震群有一個明顯的地震條帶，長約8km，地震震源深度分布在0—12km范圍內。經過對震群空間分布進行仔細分析，認為發震構造與小金河斷裂西側的一條NWW向分支斷裂有關，并不發生在NW向小金河主斷裂帶上。

（3）在震群較大的地震震源機制的研究中，分析出了兩種錯動類型，其中，發震時間靠前的淺源逆沖斷層錯動表現出受垂直方向應力（重力）作用的特征；發震時間相對較晚的走滑斷層錯動明顯與區域構造應力場（或地應力場）變化引發的構造活動增強有關。

根據上述結論可知，木里小震群活動頻次與錦屏一級水庫庫水位上升關系密切，看似僅為一般性的水庫誘發地震；但是，該地區小震群活動分布明顯成帶，且絕大多數地震震源深度較一般性水庫誘發地震震源深度偏深，表現出區域應力有增強的現象。因此，正在活動中的四川木里小震群究竟是前兆震群，還是一般性的普通震群，目前尚存爭議。

啜永清，王趙麗，孟燕英等，2000.用初動符號的變化識別前震和震群的研究.中國地震，16（3）：242—255.

崔子健，李志雄，陳章立等，2011.小震群序列類型判定方法研究的現狀及相關問題研究的趨勢.國際地震動態，（6）：32—40.

馮紅武，張元生，劉旭宙等，2013.2013年甘肅岷縣漳縣MS6.6地震及其余震序列重定位.地震工程學報，35（3）：443—447.

柯龍生，李用普，章愛娣等，1977.前震序列與一般震群活動規律的探討.地球物理學報，20（2）：1—16.

劉正榮，2004.地震預報方法.北京：地震出版社.

林向東，葛洪魁，徐平，Douglas Dreger，蘇金蓉，王寶善，武敏捷，2013.近場全波形反演：蘆山7.0級地震及余震矩張量解.地球物理學報，56（12）：4037—4047.

李永紅，胡新亮，許萍等，2011.以雙差定位方法對郯城8?級地震震中附近現代小震重新定位.西北地震學報，33（1）：71—75.

呂堅，鄭勇，馬玉虎等，2011.2010年4月14日青海玉樹MS4.7、MS7.1、MS6.3地震震源機制解發震構造研究.地球物理學進展，26（5）：1600—1606.

王俊國，刁桂苓，2005.千島島弧大震前哈佛大學矩心矩張量（CMT）解一致性的預測意義.地震學報，27（2）：178—183.

易桂喜，龍鋒，張致偉，2012.汶川MS8.0地震余震震源機制時空分布特征.地球物理學報，55（4）：1213—122.

趙根模，1990.前震和普通震群震源深度分布的比較.中國地震，6（4）：57—66.

朱麗霞，1999.1997年青海省震群活動的某些特征.西北地震學報，21（3）：67—71.

張廣偉，雷建設，2013.四川蘆山7.0級強震及其余震序列重定位.地球物理學報，56（5）：1764—1771.

趙凌云，鄧津，陳俊華等，2010.基于CAP方法的震源機制研究.長江科學院院報，27（5）：81—84.

Lund B.， B?evarsson R.， 2002. Correlation of microearthquake body-wave spectral amplitudes. Bul1. Seism. Soc. Am.， 92 （6）：2419—2433.

Wang C.Y.， Han W.B.， Wu J.P. et al.， 2007. Crustal structure beneath the eastern margin of the Tibetan Plateau and its tectonic implications. J. Geophys. Res.， 112， b07307， doi: 10.1029/2005JB003873.

Zhao L.S.， Helmberger D.V.， 1994. Source estimation from broadband regional seismograms.Bull. Seis. Soc. Amer.，84: 91—104.

Zhu L.P.， Helmberger D.V.， 1996. Advancement in source estimation techniques using broadband regional seismograms. Bull. Seis. Soc. Amer.， 86: 1634—164l.

The Active Characteristics of Small Earthquake Swarms in the Muli Area of Sichuan Province

Cai Yichuan， Dai Shigui and Wei Yaling

（Earthquake Administration of Sichuan Province， Chengdu 610041， China）

In this article， we relocated the seismic source location of the Sichuan Muli areas earthquakes， and inverted the focal mechanism of those larger earthquakes. The relationship between the water level of the Jinping reservoir and the frequency of the earthquake swarm was analyzed. The results showed that: 1）The epicenter of the relocated small earthquake swarm distributes into seismic zone， and the earthquake focal depth is in the range of 0—12 km. （2）By analyzing the earthquake swarm spatial distribution， we found that the earthquake swarms are generated from one branch fault on the west Xiaojin fault；（3）The focal mechanism of three earthquakes with magnitude greater than 4.0 is significantly different， in which the shallow source thrust events are effected by the vertical stress， and the strike-slip events are related to regional tectonic activity.

Earthquake swarm activities；Seismic zone；Double difference location method；Focal mechanism

中國地震局測震臺網青年骨干培養專項（專題編號：20140320、20140321）

2014-11-22

蔡一川，男，生于1979年。工程師。主要從事系統運維工作。E-mall：caiyichuan197910@sina.com