青海地區地震目錄最小完整性震級研究*

王培玲，胡 玉

(青海省地震局，青海西寧810001)

0 引言

青海省地處青藏高原北部地區，是一個地震頻發區。受自然條件限制，測震和前兆監測臺站分布不均勻。多年來，地震學方法在地震預測預報研究中發揮著重要的作用，利用地震目錄進行統計分析和參數計算占主導地位。但是，基于地震目錄進行的地震學預報方法在地震活動時、空、強圖像變化特征的研究中，存在許多問題，如地震目錄最小震級 (蘇有錦等，2003;焦遠碧等，1990)和地震學參數間的相關性問題 (韓渭賓等，2001)等。

地震目錄最小完整性震級MC是指地震目錄中該震級以上的地震是完整的、沒有遺漏的。地震目錄的完整性由實際地震監測能力決定。在使用地震目錄時忽略不同地區、不同時間段內的監測能力，籠統地選取一個起始震級進行地震學方法的計算研究，會影響計算結果的實用性、普適性和科學性 (沈建文，蔡長青，1998)。

本文以青海地區 (32°～40°N，89°～104°E)為研究區，開展不同統計時段最小完整性震級的研究，給出青海地區及各區不同時間段內，地震目錄的最小完整性震級。

圖1 2007年前青海省模擬地震觀測臺站分布圖Fig.1 Distribution of analogue seismic observation stations in Qinghai Province before 2007

1 青海省地震監測歷史

青海省的地震監測工作最早可追溯到1957年，由中國科學院地球物理研究所在西寧市建立了配有513型地震儀的地震臺。1964年為監視我國第一次核試驗，分別在烏蘭、格爾木、門源等地建立了地震臺，在完成核試驗監視任務之后，停止觀測。1972年青海省地震隊成立后，先后建成了門源、西寧、烏蘭、香日德、瑪多和諾木洪地震臺。1976～1981年陸續建成了大武、格爾木、玉樹、德令哈、湟源、樂都等測震臺，初步形成了地震觀測臺網。由于青海地域遼闊，地質構造復雜，受自然環境的制約，測震臺 (點)總體分布稀疏且極度不平衡，即東北部較多，西部較少，西南部更少。以現有的地震臺站可監控72×104km2范圍地域內的地震，其監控能力遠遠不夠，降低了觀測資料的科學性和可用性，不能滿足地震監測、預報和科研工作的需要。“九五”期間對格爾木、西寧、花土溝、湟源和都蘭地震臺進行了數字化觀測技術的改造。由于改擴建的臺項少，無法組網，不能發揮其應有的功能，青海省仍然以模擬地震臺網為主，用于日常地震觀測的臺站僅有8個 (圖1)，監測能力仍然很低 (青海省地震局，2005)。2007年9月青海省建成了數字地震觀測臺網。截至2013年底青海數字地震觀測臺網中心，能夠接收青海及鄰區省份共47個臺站 (圖2)的數字地震觀測資料，大大地提高了青海地震監測能力，也提高了西部和西南部的地震監測能力。

圖2 2007年后青海省數字地震觀測臺站分布圖Fig.2 Distribution of digital seismic observation stations in Qinghai Province after 2007

圖3 1970年以來青海地區地震事件震級—序號圖Fig.3 Magnitudes and related order number diagram of earthquakes in Qinghai region since 1970

2 地震目錄最小完整性震級分析

青海地區使用的中小地震目錄始于1970年，由于社會經濟文化發展不平衡、臺站布局和建設不完善，造成了地震目錄中不同時間段和不同地區的中小地震的缺失，即地震目錄的不完整。因此在開展地震學方法研究工作前，首先對所使用的地震目錄進行完整性分析 (何宗海，1994;劉麗芳等，2012;李志海，2011;秦娟，2012)，以保證后續使用地震目錄開展統計和參數計算研究的正確性和科學性。

由于G-R關系計算最小完整性震級 (MC)的方法基于簡單的分布假定，計算中便于操作，因此本文主要利用該方法對青海地區內不同時間段、不同區域的地震目錄進行最小完整性震級的計算研究。還基于G-R關系對計算MC的其他方法，如定性評估的“震級—序號”法 (Ogata et al.，1991)、定量評估的“最大曲率”MAXC方法及擬合度分別為90%和95%的GFT方法 (Wiemer S，Wyss，2000)進行了計算研究和佐證。

本文以1970年以來青海地區 (32°～40°N，89°～104°E)小震目錄為例，進行MC計算，研究其時序變化特征 (馮建剛等，2012;王鵬等，2011)

2.1 “震級—序號”法

利用“震級—序號”法按地震發生時間的先后順序排序，考察不同震級的地震數密度分布來定性分析1970～2013年青海地區小震目錄的MC。從圖3可以看出，在1970～1980年間，青海地區小震主要集中在2.2～3.5級之間，1980～1990年間主要集中在2.0～3.0級之間，1990～2008年主要集中在1.5～2.5級之間，2008年至今主要集中在1.0～2.0級之間。由此看出，隨著監測能力的提高，特別是2007年9月青海省數字地震臺網正式運行后，M≤2地震記錄明顯增多，目錄的最小完整性震級MC在1.0～2.0之間。

2.2 多種方法定量分析

利用MAXC方法和GFT方法進行多方法的定量分析。其中MAXC方法是將震級—頻度曲線一階導數最大值對應的震級作為MC。GFT方法是通過搜索給定的實際與理論震級—頻度分布下的擬合程度百分比來確定MC。這里分別稱為 GFT-90%和GFT-95%。由于不同方法同時使用，結果不唯一，產生優選問題，設定優先級為 GFT-95% ＞GFT-90% ＞MAXC，選取可計算且優先級高的MC-Best為最終結果。

圖4給出了最大曲率法MAXC(紅色)的計算結果，擬合度分別為90%(藍色)和95%(藍綠色)的檢測結果，利用多方法的定量分析，得出最佳最小完整性震級MC-Best(黑色)，即根據現有的地震目錄，能計算得到的最接近理論的線性震級—頻度關系的擬合，從圖中可以看出2007年以后，擬合度95%的曲線 (藍綠色)和擬合度90%的曲線都存在，也就是說2007年后的地震目錄在進行多種方法的定量分析時，能夠滿足設定的優先級 (GFT-95% ＞GFT-90% ＞MAXC)，最優起算震級為M≤2。由此進一步說明青海省數字地震臺網正式運行前 (1970～2007年)，由于臺站稀少且分布不均勻，監測能力比較低，很難達到和滿足設定的優先級別。而2008年以后，由于數字地震臺網的運行，臺站數量的增多和分布的改善，使得青海地區的監測能力得到了提高。

圖4 青海地區最小完整性震級M C時序變化曲線 (多方法)Fig.4 Temporal variation curve of the minimum completed magnitudes(M C)in Qinghai region(multi-methods)

圖5 1970～2013年青海地區小震目錄G-R關系擬合曲線Fig.5 The fitted curve of G-R relation of small earthquakes catalogue in Qinghai region from 1970 to 2013

2.3 G-R關系分析

基于地震目錄和統計地震學的MC計算方法主要是基于震級不小于MC的地震在震級—頻度分布上滿足 G-R關系的假定 (Gutenberg，Richter，1944)，即

式中，M是地震震級，N是研究區內震級大于等于M的總次數，a和b為統計常數。當震級—頻度分布中，能夠滿足G-R關系時，對應的最小起始震級即為MC，并認為這些地震的記錄是完整的。

選取1970～2013年間青海地區小震目錄進行G-R關系計算，從圖5中雖然可以看出M≤1.3地震明顯偏離直線關系，但是卻不能將1.3作為地震目錄的MC。因為在1970～2013年間，青海省臺網進行了多次調整和架設，直到2007年9月青海數字地震臺網投入運行后，地震臺站數量成倍增長、布局得到改善，提高了地震臺網分布密度，但是實際地震監測能力仍達不到該水平。這個震級下限反映了40多年來青海地區的地震活動優勢地區或時段的監測能力，而可能將監測能力較差地區或時段的實際情況掩蓋掉。所以本文在進行G-R關系分析時，將根據不同時段，對不同區域進行分析，結合相關系數R，選取擬合效果最好的起算震級作為該時段內的最小完整性震級MC。

圖6 1970～2013年青海地區小震目錄不同最小完整性震級相關系數變化曲線Fig.6 Correlation coefficent variation curve of different minimum completed magnitudes of small earthquake catalog in Qinghai region from 1970 to 2013

根據1970～2013年青海省小震目錄，通過G-R關系擬合得出不同最小完整性震級的相關系數可以看出 (圖6)，當起算震級為2.0時，得到的相關系數最高，G-R關系擬合效果最好，因此這一時間段內，青海地區小震目錄的最小完整性震級可以選為2.0級。

從2000～2013年，青海省小震目錄的G-R關系擬合曲線 (圖7)可以看出，當M≤1.5時，樣本點明顯偏于直線下方，因此取最小完整性震級為1.5。由此得出了不同時間段、不同區域內，小震目錄的最小完整性震級 (表1)。

從表1中可以看出，不同時間段、不同區域地震目錄的最小完整性震級存在明顯差異，如同一時間段內唐古拉地震帶與祁連地震帶相比，唐古拉地震帶最小完整性震級明顯高于祁連地震帶，這也反映出唐古拉地震帶的實際監測能力明顯低于祁連地震帶。自2008年后唐古拉地震帶上由原先的沒有地震監測臺增加到5個臺，監測能力得到提高。

圖7 2000～2013年青海地區小震目錄G-R關系擬合曲線Fig.7 The fitted curve of G-R relation of small earthquakes in Qinghai region from 2000 to 2013

表1 不同時間段、不同地區地震目錄最小完整性震級Tab.1 Minimum completed magnitude of earthquake catalogue in different time at different area

3 討論與結論

本文介紹了青海省的地震觀測發展歷史和現狀，對青海省地震觀測臺網臺站布設、地震記錄進行了梳理。但由于青海省地震觀測臺網布設歷史記載不詳實，無法準確掌握臺站布設增減的具體時間，只能籠統的將青海省地震觀測發展時期按模擬記錄和數字記錄加以區分，即2008年前為模擬地震臺網觀測，2008年后為數字地震臺網觀測。從地震震級序號法和多種方法定量分析，都可以明顯的看出2008年前后，青海省地震目錄的最小完整性震級在1.0～2.0之間。而利用G-R關系計算出2008年以后地震目錄的最小完整性震級為1.5。由此可以看出數字地震觀測臺網大大地提高了青海省的地震監測能力，有些區域不同時間段內的最小完整性震級存在明顯的差異，如庫瑪地震帶內1990～2013年和2000～2013年兩個時間段內，雖然不存在地震臺網架設增減問題，但是，這兩個時間段相比，可能2000～2013年是M≥1.3地震的活動優勢活動時間段。而1990～2013年M≤1.5地震相對比較少。由此認為本文最終給出的最小完整性震級是合理的。

筆者根據日常地震預報工作中選取時間段的實際使用情況，人為地對地震目錄進行時間段的劃分可能不夠嚴謹，但是這樣劃分僅為了方便日常預報工作中利用地震目錄進行測震學參數時間掃描計算和地震危險性分析等工作要求。

本文計算程序由蔣長勝博士提供，在此表示衷心感謝!

馮建剛，蔣長勝，韓立波，等.2012.甘肅測震臺網監測能力及地震目錄完整性分析［J］.地震學報，34(5):646-658.

韓渭賓，易桂喜，辛華.2001.與頻度或能量有關的各種地震學指標的相關性研究［J］.地震學報，23(6):582-587.

何宗海.1994.中國南北地震帶中北部最小完整性震級的探討［J］.西北大學學報，24(5):411-416.

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

李志海，蔣長勝，黃瑜，等.2011.新疆地區地震目錄最小完整性震級和臺網科學布局研究［J］.地震學報，33(6):1-13.

劉麗芳，李志海，蔣長勝.2012.云南地區地震目錄最小完整性震級研究［J］.地震研究，35(4):491-498.

秦娟，蔡輝騰，李光，等.2012.重慶地區小地震完整性分析［J］.華南地震，32(1):64-70.

青海省地震局.2005.青海省地震監測志［M］.北京:地震出版社，9-10.

沈建文，蔡長青.1998.起算震級與衰減不確定性校正［J］.中國地震，14(4):36-43.

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

王鵬，鄭建常，趙金花，等.2011.山東地區最小完整性震級的確定—R-S檢驗［J］.華北地震科學，29(3):19-22.

Gutenberg R.，Richter C.F..1944.Frequency of Earthquakes in California［J］.BSSA，34:185-188.

Ogata Y.，Imoto M.，Katsura K..1991.13-D Spatial Variation of b-values of Magnitude-frequency Distribution beneath the Kantodistrict，Japan［J］.Geophys J.Int.，104(1):135-146.

Wiemer S.，Wyss M..2000.Minimum Magnitude of Complete ness Repor-ting in Earthquake Catalogs:Examples from Alaska，the Western United States，and Japan［J］.Bull Seism.Soc.Am.，90(4):859-869.