內蒙古呼和浩特—包頭地區地震目錄的完整性分析①

韓曉明，曹孟娜，張 帆

（內蒙古自治區地震局，內蒙古 呼和浩特 010010）

0 引言

在地震活動性研究中地震目錄的完整性至關重要。由于受到臺站布局、觀測儀器性能、觀測人員分析能力等不確定因素的影響，地震目錄往往是不完整和非均勻的，即地震目錄中除了天然地震之外，還摻雜著些許人為地震。用這種非純天然的地震目錄去研究某一區域的地震活動性是不合適的，有必要對其中的人為地震事件加以剔除，保證地震活動性研究結果是建立在地震目錄的完整性和均勻性的基礎上。關于人為地震產生的原因，Habermann和傅征祥都做了系統歸納，主要有3種因素：①臺網監測能力的變化；②地震監測報告的缺失；③由臺站分布和震級測定算法不同而引起的震級漂移。前兩種因素主要影響地震目錄中的小震級事件的完整程度。因此，需要確定一個震級閥值來保證在不同時期之間的地震目錄都是完整的，這個震級閥值就是最小完整性震級Mc（Minimum magnitude of completeness），它是指特定尺度的時—空體內100%的地震事件可被檢測到的最低震級［1］，是評價地震目錄的一個重要參數。

國內外關于地震目錄完整性的研究由來已久，且方法眾多。概括來講，主要研究最小完整性震級在時間上的變化規律，空間上的變化規律，或者同時研究其時間和空間變化。而在各自的時間范圍和空間范圍選取上又有著明顯的不同。作者分析認為，應主要依據各地區斷裂構造展布、地震活動空間分布、以及臺網改造時間等因素而選取合適的時間或空間范圍。焦遠碧等根據我國地震臺網監測能力繪制了全國最小完整性震級空間分布等值線圖［2］；黃瑋瓊等分析了中國大陸不同地震區（帶）歷史地震的平均發生率和震級—頻度分布關系，給出了中國大陸分區地震目錄（MS≥4.7和5以上）完整性的時間分布［3－4］。Stefan Wiemer等用震級—頻度分布FMD （frequency－ magnitude distribution）方法研究了阿拉斯加、美國西部、日本地區的地震目錄最小完整性震級的空間分布［5］。蘇有錦分析了川滇地區1970－2001年的區域地震目錄，給出了該地區的最小完整性震級及主要地震活動帶（區）最小完整性震級的時、空分布［6］。

1 研究區域及資料選取

內蒙古地區地域遼闊，臺站分布稀疏、非均勻性顯著。部分地區監測能力很弱，少有地震記錄，對其進行地震目錄的完整性分析意義不大。相對這些臺網監測能力偏低的地區，呼和浩特—包頭地區（以下簡稱呼包地區）是內蒙古地震監測能力最強的區域?！熬盼濉逼陂g，呼包地區建成了由1個臺網中心、7個臺站組成的呼和浩特數字遙測地震臺網，基本覆蓋了呼包地區；“十五”期間，呼和浩特數字遙測臺網進一步擴建調整，對區域地震臺站進行數字化改造，主要包括觀測設備、傳輸信道等的更新和升級，使該區域的臺網技術水平、監測能力、觀測精度和觀測數據的利用率都得到大幅提升。既保證該區域有足夠多的地震樣本數目參與計算，也能客觀的對其地震目錄的完整性進行精細分析。

表1 呼和浩特—包頭地區數字測震臺網臺站基本參數表

根據全國地震目錄，并結合內蒙古地震臺網目錄，經過一致性處理，選取時間范圍為2001年1月1日－2010年4月30日，空間范圍為40°～42°N，108°～114°E，震級下限 ML1.0，以2個月時間窗長進行了時間分布統計；以震級間隔為ML0.1進行了震級分布統計（圖3）。共計1 072個地震事件參與計算，震級主要集中于ML1.5～2.5之間，最大震級ML5.0（圖1、圖2、圖3（b））。從圖3（a）可以看出，2003－2005年呼包地區地震活動水平相對較高，2006－2008年地震活動水平最低，2009年以來略有回升?？紤]到選取的空間范圍較小，因此本文只做最小完整性震級在時間上的變化分析，未做空間變化分析。

圖1 呼包地區地震事件的時間分布（ML≥1.0，2001－01－01－2010－04－30）Fig.1 Time distribution of seismic events in Hohhot－Baotou area（ML≥1.0，2001－01－01－2010－04－30）.

圖2 呼包地區斷裂構造及地震事件的空間分布（ML≥1.0，2001－01－01－2010－04－30）Fig.2 Spatial distribution of faults and seismic events in Hohhot－Baotou area（ML≥1.0，2001－01－01－2010－04－30）.

圖3 呼包地區地震事件的統計（ML≥1.0，2001－01－01－2010－04－30）Fig.3 Statistics of seismic events number for every two years and for every 0.1magnitude in Hohhot－Baotou area（ML≥1.0，2001－01－01－2010－04－30）.

2 呼包地區地震目錄的完整性分析

2.1 最小完整性震級Mc的確定及時間變化分析

基于Rydelek等［1］對最小完整性震級Mc的具體定義，在上述選取的時—空—強范圍內，根據地震震級—頻度公式logN＝a－bM，對2001年以來呼包地區ML≥1.0地震進行掃描計算，掃描窗長為100個地震事件，步長為50個地震事件。計算結果顯示，2001年1月－2010年4月，呼包地區平均最小完整性震級Mc＝1.8，基本浮動范圍為ML1.6～2.2，即最小完整性震級的最大值為2001年度的ML2.2，最小值為2010年度的 ML1.6（圖4、圖5）。實際上，由于地震臺站數目隨時間在不斷變化（總體上不斷增加，但由于個別時段的整改停測，所以有時候數目會表現為減少）以及地震分析方法的改進，地震目錄最小完整性震級Mc隨著時間的流逝應該是非均勻性減小的。如圖5中所示，呼包地區的Mc值隨著時間流逝的總體變化趨勢雖是逐漸減小的，但也有個別時段的變化是相反的，2005－2006年的Mc值就有著轉折升高現象，這主要取決于該時段內呼包地區地震活動特點、臺站布局、地震信噪比和人為因素等。

圖4 呼包地區地震事件的最小完整性震級Mc（ML≥1.0，2001－01－01－2010－04－30）Fig.4 The minimum integrality magnitude of seismic events in Hohhot－Baotou area（ML≥1.0，2001－01－01－2010－04－30）.

圖5 呼包地區Mc隨時間的變化形態（ML≥1.0，2001－01－01－2010－04－30）Fig.5 Morphological changes of Mcwith time in Hohhot－Baotou area（ML≥1.0，2001－01－01－2010－04－30）.

除了定量的分析Mc值隨著時間的變化趨勢，還試圖利用震級序號圖對Mc的時序變化做定性分析。震級序號圖方法是按照地震發生的先后順序對地震事件進行排序，由地震序號和對應的震級來討論地震目錄完整性的變化。這種方法主要排除地震叢集事件等因素的影響，可以反映特定時空范圍內區域臺網監測能力的短期或者不連續的變化。本文選取2001年1月1日至2010年4月30日呼和浩特—包頭地區的小震事件作為研究對象，震級下限ML＝1.0，按照一天內地震發生的時間先后進行編號排序，構成震級－序號分布圖。可能由于本文選定的時空范圍過小，地震樣本數目沒有很好的滿足計算要求，圖6中所示的震級序號圖效果較差，用于表征各個震級—時間范圍內地震數目的色澤分布不均勻、不密實，未能直觀的表達呼包地區2001年以來臺網監測能力的實際變化，這里僅作參考。

圖6 呼包地區地震事件的震級序號圖（ML≥1.0，2001－01－01－2010－04－30）Fig.6 The magnitudes and event sequence diagram in Hohhot－Baotou area（ML≥1.0，2001－01－01－2010－04－30）.

2.2 RS檢驗結果

最小完整性震級雖然能夠在給定的臺網監測能力下計算出來，但由于地震臺網所在地理、社會、人文環境下的自然噪聲和文化噪聲的影響在個別時段內是非常強烈的，這就可能導致理論上應該能被記錄到（大于最小完整性震級）的小震級事件被這些噪聲背景所淹沒，從而反過來影響最小完整性震級的確定。有鑒于此，Rydelek等［1］提出了一種考慮噪聲背景進行地震目錄完整性檢驗的方法，即Rydelek－Sacks檢驗（簡稱RS檢驗），這種方法的實現必須滿足以下2個前提條件：（1）在任何震級檔的地震事件都必須服從泊松分布；（2）白天的噪聲水平要強于夜間的噪聲水平。

在給定的震級M0水平內，一個地震的發震時間可以在一個24小時的“時鐘”上用單位振幅的矢量來表示，它具有一個相位角，即一個地震事件可以表示成一個相位矢量，那么選取的時空范圍內的地震事件按照發生的先后順序就組成一個相位矢量軌跡圖，也即所有地震事件的相位相加得到一個總相位矢量，其模為R，它反應了地震總體上的發生規律。按照第一個假設，即如果地震的發生是隨機的，其相位矢量軌跡應該做布朗運動，那么得到一個模大于R的總相位矢量的概率［1］：

其中N是地震總數。因此，如果R超過了某個臨界值RC（在“時鐘”上表現為臨界半徑），那么所考慮的目錄便具有周期約為24小時的周期性。對于95%的置信水平，理論上有［1］

如果R≥RC，則可以認為地震目錄受周期調制，因而小于等于給定震級M0的地震目錄都是不完整的。即可以根據所記錄到的地震事件是否可被優勢周期約為24小時的周期過程調制來判斷地震目錄是否完整，假如地震目錄出現周期性，則認為是不完整的。

圖7是呼包地區2001年以來ML≥1.0地震事件的RS檢驗結果，這個結論與圖5中得出的結果是相一致的。圖7顯示，2001－2010年，呼包地區1.0≤ML＜2.2地震事件的相位矢量R≥RC，即這個震級范圍內的地震目錄是不完整的；ML≥2.2的地震事件的相位矢量R＜RC，即ML≥2.2的地震目錄是完整的。即如果選取2001年以來完整的地震目錄，震級下限只能取Mc的最大值ML2.2，而不能取最小Mc值1.6，這樣才能保證2001－2010年期間各個時段內的地震目錄都是完整的。

圖7 呼包地區地震目錄完整性的RS檢驗（ML≥1.0，2001－01－01－2010－04－30）Fig.7 The RS test for integrity of the earthquake catalog in Hohhot－Baotou area（ML≥1.0，2001－01－01－2010－04－30）.

3 結論和討論

首先需要說明的一個問題是，臺網監測能力為某一震級，是指當這個地區發生這一震級以上的地震時，至少有4個臺站能記錄到震相完整且清晰可辨的地震波形。因此，最小完整性震級與臺網監測能力所指的震級有一定的相關性，即最小完整性震級能模糊的反映某一地區的臺網監測能力，但與精確表達臺網監測能力的特定震級在本質上還是有區別的。二者基本上屬于理論值和實際值之間的關系，可以通過臺網監測能力震級來佐證所計算的最小完整性震級Mc的準確性。而根據內蒙古地區測震臺網監測能力空間分布（圖8），呼包地區的監測能力應該能達到ML1.5～2.0，這個震級范圍基本與本文計算的最小完整性震級范圍（ML1.6～2.2）相符，也與通過RS檢驗得出的結果相吻合。

圖8 內蒙古自治區地震監測能力示意圖（內蒙古自治區地震局臺網中心提供）Fig.8 Schematic diagram of seismic monitoring capability in Inner Mongolia Autonomous Region（from Earthquake Administration of Inner Mongolia Autonomous Regional Network Center）.

本文著重分析Mc值在時間上的變化規律，通過Mc時序曲線（圖5）可以發現各時段對應的臺網監測能力變化和地震活動特點，但其中的變化也包含著些許自然和人為因素的干擾，所以不能作為精確指標對臺網監測能力和區域地震活動性做出評價和判斷。

目前關于最小完整性震級Mc的計算，都是在基于文中所述的2個假設前提下進行的，特別是根據第一個假設而衍生出很多的計算方法，如完整震級范圍的EMR（Entire－magnitude－range method）方法、最大曲率方法 MAXC（Maximum curvature－method）、擬合度檢測法 GFT（Goodness－of－fit test）和b值穩定性判定Mc的 MBS（Mcby bvalue stability）方法等。為了使Mc的測算和評估更加精確，還需綜合利用多種方法加以驗證，這需要在今后的工作中逐步完成。

致謝：中國地震局地球物理研究所蔣長勝博士提供了RS檢驗和震級序號圖計算程序，在此深表謝意！

［1］Rydelek P A，I S Sacks.Testing the completeness of earthquakecatalogs and the hypothesis of self－similarity［J］.Nature，1989，337：251－253.

［2］焦遠碧，吳開統，楊滿棟.我國地震臺網監測能力及臺網觀測條件質量評定［J］.中國地震，1990，6（4）：1－71

［3］黃瑋瓊，李文香，曹學鋒.中國大陸地震資料完整性研究之一：以華北地區為例［J］.地震學報，1994，16（3）：274－280.

［4］黃瑋瓊，李文香，曹學鋒.中國大陸地震資料完整性研究之二：分區地震資料基本完整的起始年分布圖象［J］.地震學報，1994，16（4）：423－432.

［5］Stefan Wiemer，Max Wyss.Minimum Magnitude of Completeness in Earthquake Catalogs：Examples from Alaska，the Western United States，and Japan［J］.BSSA，2000，90（4），859－869.

［6］蘇有錦，李永莉，李忠華，等.川滇地區區域地震目錄完整性最小震級分析［J］.地震研究，2003，26（增刊）：10－16.

［7］劉杰，陳棋福，陳顳.華北地區地震目錄完全性分析［J］.地震，1996，16（1）：59－67.

［8］何宗海.中國南北地震帶中北部最小完整性震級的探討［J］.西北大學學報（自然科學版），1994，24（5）：411－416.

［9］Cao A M，S S Gao.Temporal variation of seismic b－valuesbeneath northeastern Japan Island arc［J］.Geophys.Res.Lett.，2002，29：no.9，doi 10.1029／2001GL013：775.

［10］Ogata Y，K Katsura.Analysis of temporal and spatial heterogeneityof magnitude frequency distribution inferred from earthquakecatalogs［J］.Geophys.J.Int.，1993，113：727－738.

［11］李建彪，甘衛軍，冉永康，等.青藏高原東部構造塊體的運動學及變形特征分析［J］.西北地震學報，2006，28（2）：97－103.

［12］韓曉明，薛丁，劉文廣，等.喀喇沁地區震群活動與周圍區域中強地震的對應關系研究［J］.西北地震學報，2009，31（4）：380－384.