中國大陸特征地震震級分布規(guī)律*

陳中天 郭 星 潘 華 李金臣

1) 中國北京100081中國地震局地球物理研究所 2) 中國北京100082環(huán)境保護部核與輻射安全中心

?

中國大陸特征地震震級分布規(guī)律*

陳中天1)郭 星2),*潘 華1)李金臣1)

1) 中國北京100081中國地震局地球物理研究所 2) 中國北京100082環(huán)境保護部核與輻射安全中心

根據(jù)特征地震的震級分布符合正態(tài)分布的基本假設， 對特征地震的震級分布規(guī)律進行定量化研究． 考慮到原地復發(fā)型大地震的震級數(shù)據(jù)稀少， 提出利用大地震的同震位移轉換得到震級數(shù)據(jù)， 再對不同觀測點所得到的震級數(shù)據(jù)進行歸一化處理， 使所有歸一化震級數(shù)據(jù)均服從均值為0的同分布， 進而可以統(tǒng)計得到反映特征地震震級變異性的標準差σ． 最后基于廣泛搜集的中國大陸54個特征地震的同震位移， 利用本文給出的統(tǒng)計方法得到一個通用的標準差σ=0.08， 為進一步研究特征地震的發(fā)生概率模型提供了重要依據(jù)．

特征地震 震級 同震位移 小樣本

引言

Aki(1984)和Schwartz和Coppersmith(1984)以彈性回跳理論(Reid， 1910)為基本物理假設， 在分析圣安德烈斯斷層上反復發(fā)生的大地震資料基礎上， 提出了“特征地震”的概念， 即特征地震是指在某些斷層長期活動的過程中， 反復發(fā)生的并往往表現(xiàn)出相似破裂長度、 位錯分布和震級大小的多次大地震．

考慮到特征地震的震級存在不確定性， 一般認為特征地震震級的變化區(qū)間大約在0.5個震級單位(聞學澤， 1995； 楊明， 劉百篪， 2000； 陳立春， 2002)， 但相關的定量化研究鮮有報道． 本文基于特征地震模型和特征地震震級符合正態(tài)分布的基本假定(Schwartz， Coppersmith， 1984； Youngs， Coppersmith， 1985)， 對特征地震震級分布規(guī)律進行定量化研究． 本研究不討論其它有別于特征地震的大地震破裂模型， 那些不符合特征地震模型的斷層源也不予考慮．

欲研究特征地震的震級分布規(guī)律， 需對不同斷層源上特征地震的震級進行歸一化處理， 再對歸一化的數(shù)據(jù)進行統(tǒng)計分析． 然而， 對于某個斷層源， 很難獲得兩個以上特征地震的精確震級數(shù)據(jù)， 但考慮到同震位移的大小直接與震級大小相關， 我國很多斷層上可以觀測到大量同震位移數(shù)據(jù)， 因此， 本文提出利用特征地震的同震位移數(shù)據(jù)轉換得到歸一化震級數(shù)據(jù)， 進而通過對歸一化震級數(shù)據(jù)的統(tǒng)計分析來對特征地震的震級分布作定量化研究．

此外， 考慮到每個觀測點上的地震事件個數(shù)均非常有限， 基于多組小樣本數(shù)據(jù)的統(tǒng)計分布與實際分布也可能存在一定的差異， 本研究對這種差異性進行定量化分析， 并對每個歸一化震級數(shù)據(jù)進行修正， 使所有震級數(shù)據(jù)符合同分布．

本文廣泛搜集了中國大陸54個特征地震的同震位移數(shù)據(jù)， 擬將同震位移數(shù)據(jù)轉換為歸一化無量綱震級數(shù)據(jù)， 并對每個無量綱震級數(shù)據(jù)進行修正； 然后通過統(tǒng)計分析， 計算得到歸一化震級分布的標準差σ．

1 利用同震位移數(shù)據(jù)研究特征地震的震級分布規(guī)律

一般定量化研究特征地震的震級分布規(guī)律， 需對不同斷層源上的特征地震震級進行歸一化處理， 再對歸一化的數(shù)據(jù)進行統(tǒng)計分析得到． 然而， 大地震的復發(fā)周期一般都很長， 人類有較精確的大地震震級記載的時間又很短， 任何斷層源上都很難得到兩個以上特征地震的精確震級數(shù)據(jù)． 但是， 很多斷層上可以通過地質手段觀測得到大量同震位移數(shù)據(jù)， 而同震位移的大小直接與大地震的震級大小相關， 因此可以考慮利用斷層上的同震位移數(shù)據(jù)來研究特征地震的震級分布規(guī)律．

1.1 利用特征地震的同震位移數(shù)據(jù)確定歸一化震級數(shù)據(jù)

地震矩與斷層運動所釋放的能量直接相關(Kanamori， 1977)， 其中地震矩可以表示為

(1)

式中：M0為地震矩， 單位為N·m；μ為斷層附近的巖石剪切模量, 單位為N/m2；A為斷層面積， 單位為m2;D為平均同震位移， 單位為m．

地震矩M0與震級M存在以下經驗關系：

lgM0=1.5M+C，

(2)

式中C為常數(shù)項， 根據(jù)郭星(2014)對震級-地震矩關系的統(tǒng)計分析可知， 當震級為我國面波震級MS時C=8.61， 當震級為矩震級MW時C=9.05．

聯(lián)合式(1)和式(2)， 可得

(3)

圖1 特征地震模型(引自Schwartz， Coppersmith， 1984)Fig.1 Characteristic earthquake model (after Schwartz， Coppersmith， 1984)

若基于Schwartz和Coppersmith(1984)提出的特征地震模型(圖1)， 對于同一斷層或同一斷層段， 式(3)中的lg(μA)為一常數(shù)， 則震級M與平均同震位移D之間的關系可表示為

(4)

(5)

對于任一平均同震位移數(shù)據(jù)Dk， 根據(jù)式(4)和式(5)可以得到一個歸一化無量綱震級數(shù)據(jù)Mk， 即

(6)

大量古地震研究(錢洪等， 1990； 王華林， 1995； 萬永魁等， 2014)表明： 大地震在斷層上各處的同震位移符合特征地震模型， 在斷層或斷層中間段的同震位移較大， 至兩端則趨于零； 對于斷層上反復發(fā)生的多次特征地震事件， 不同位置處的同震位移數(shù)據(jù)則保持固定的比例關系不變． 因此， 對于式(6)， 并不要求Dk必須為平均同震位移， 僅保證在歸一化無量綱震級數(shù)據(jù)計算過程中所使用的同一斷層上的同震位移數(shù)據(jù)是出于同一觀測點即可．

1.2 歸一化震級數(shù)據(jù)的不確定性處理

由于每個觀測點上的地震事件個數(shù)非常有限， 而基于多組小樣本數(shù)據(jù)的統(tǒng)計分布與實際分布可能存在一定的差異(郭星， 潘華， 2015)． 為此， 本文設定初始正態(tài)分布的均值μ=0， 標準差σ=0.01， 利用蒙特卡羅方法給出了分別對應于每組抽取2次、 3次、 4次、 5次地震同震位移數(shù)據(jù)的Mk值分布(圖2)， 隨機抽取的次數(shù)為2萬次. 結果表明： 每組抽取樣本數(shù)不同， 所對應的Mk值分布的標準差也是不同的． 這樣對于同震位移數(shù)據(jù)量不同的觀測點， 所統(tǒng)計得到的Mk值分布就不是獨立同分布的， 而不是同分布的數(shù)據(jù)很難放在一起進行統(tǒng)計分析． 因此， 在對實際歸一化震級數(shù)據(jù)Mk進行統(tǒng)計分析的過程中， 需要對數(shù)據(jù)進行修正， 使多組小樣本震級數(shù)據(jù)盡量滿足同分布．

為了對每一個歸一化無量綱震級數(shù)據(jù)進行修正， 本文定義修正系數(shù)ρ為初始分布的標準差與多組小樣本數(shù)據(jù)統(tǒng)計分布的標準差之比． 在最終的統(tǒng)計分析之前， 需要根據(jù)該組同震位移數(shù)據(jù)的個數(shù)， 對計算得到的每個歸一化無量綱震級數(shù)據(jù)進行修正， 即乘上一個修正系數(shù)ρ． 經過反復檢驗， 結果顯示： 初始分布標準差不同， 計算所得到的修正系數(shù)ρ也會有一定差異， 雖然這種差異無規(guī)律， 但差異很小， 對最終的計算結果影響也很?。?因此本文直接選取初始正態(tài)分布(μ=0，σ=0.01)計算得到一組修正系數(shù)：ρ(2)=1.41，ρ(3)=1.22，

圖2 抽取2次(a)， 3次(b)， 4次(c)， 5次(d)地震事件的Mk值初始

ρ(4)=1.15，ρ(5)=1.12．

2 中國大陸地區(qū)特征地震震級分布的標準差計算

2.1 數(shù)據(jù)選取

通過對斷層的實際觀測可以獲得大量同震位移數(shù)據(jù)， 但并不是所有同震位移數(shù)據(jù)均屬同一級別的特征地震． 例如， 海原斷裂帶分為全段破裂、 中段破裂、 西段破裂等3種破裂方式， 其中全段破裂的震級高， 同震位移大， 而中段和西段破裂則屬于次一級破裂(冉永康， 鄧起東， 1998)． 此外， 并非所有大地震的發(fā)生均滿足準周期模式， 對于一些大地震的復發(fā)規(guī)律異?；虿幻鞔_的斷層， 大地震的發(fā)生是否符合特征地震模型尚需進一步考證． 基于以上考慮， 本研究確定了選取大地震同震位移數(shù)據(jù)的基本原則如下：

1) 斷層源或斷層分段上大地震的復發(fā)必須符合準周期模型；

2) 同震位移數(shù)據(jù)須屬于同一個斷層源或斷層分段， 不可以是不同級別的特征地震， 若無法明確是否為同一級別的大地震， 則不可選用；

3) 每個斷層上的一組同震位移數(shù)據(jù)必須為同一觀測點所獲取．

根據(jù)以上原則， 本文在中國大陸16個不同斷層上選取了共計54個同震位移數(shù)據(jù)， 具體列于表1， 其中夏墊斷裂上有來源于不同觀測點的兩組數(shù)據(jù)， 分別作為單獨的數(shù)據(jù)集參與統(tǒng)計分析．

表1 本文所用的中國大陸16個斷層上的同震位移數(shù)據(jù)

續(xù)表1

圖3 本文所統(tǒng)計的54個特征地震的Mk值分布直方圖Fig.3 Histogram of Mk values for 54 characteristic earthquakes

2.2 計算結果

基于本文搜集的54個同震位移數(shù)據(jù)， 根據(jù)式(6)可以計算得到54個歸一化無量綱震級數(shù)據(jù)， 并對其進行修正， 得到歸一化無量綱震級數(shù)據(jù)Mk的分布圖， 如圖3所示．

假設特征地震的震級分布滿足正態(tài)分布， 對54個修正后的歸一化無量綱震級數(shù)據(jù)進行統(tǒng)計分析， 結果顯示其均值為0， 標準差σ為0.08． 圖3的直方圖直觀地顯示出，Mk的分布接近于正態(tài)分布， 而且此數(shù)據(jù)集可通過顯著性水平為0.05的正態(tài)分布假設檢驗． 考慮到同震位移數(shù)據(jù)量較少， 要獲取更為精確的特征地震震級分布， 尚需積累更多更精確的同震位移觀測資料．

根據(jù)郭星(2014)對特征地震復發(fā)周期的變異性研究， 若取強震復發(fā)周期變異系數(shù)大于0.35， 假設地震矩累積率是恒定的， 強震復發(fā)周期的變異性均由特征地震的震級不確定性所造成， 則可計算得到特征地震震級分布的標準差應大于0.1． 根據(jù)本文的統(tǒng)計結果標準差為0.08， 則可推測出斷層上的地震矩累積率并不是恒定不變的．

3 討論與結論

研究特征地震的震級分布規(guī)律， 可以為研究大地震復發(fā)的地震矩累積和釋放過程提供重要依據(jù)． 考慮到很難直接獲取足夠多的特征地震震級數(shù)據(jù)， 本文提出利用同震位移數(shù)據(jù)研究特征地震震級的分布規(guī)律． 由于未涉及其它有別于特征地震的大地震破裂模式， 本研究具有一定的局限性， 需要同特征地震模型一起不斷完善和更新．

由于每個斷層觀測點上的同震位移數(shù)據(jù)量非常有限， 基于多組小樣本數(shù)據(jù)的統(tǒng)計分布與實際分布可能存在一定的差異， 因此本文定義初始分布的標準差與多組小樣本數(shù)據(jù)分布的標準差之比為修正系數(shù)ρ． 在最終統(tǒng)計分析之前， 需要根據(jù)該組同震位移數(shù)據(jù)的個數(shù)， 對每個計算所得到的歸一化無量綱震級數(shù)據(jù)進行修正， 即乘以修正系數(shù)ρ， 使得修正后的數(shù)據(jù)盡可能符合同分布.

本文在中國大陸16個不同斷層上選取了共計54個同震位移數(shù)據(jù)， 將同震位移數(shù)據(jù)轉換為歸一化無量綱震級數(shù)據(jù)， 并對每個無量綱震級數(shù)據(jù)進行了修正． 基于特征地震的震級分布滿足正態(tài)分布的基本假設， 統(tǒng)計得到其標準差為0.08， 據(jù)此推測斷層上地震矩累積率并不是恒定不變的， 而對于地震矩累積率變異性的定量化研究將在下一步工作中開展．

柴熾章， 廖玉華， 張文孝， 許文俊， 申旭輝， 田勤儉． 2001． 靈武斷裂晚第四紀古地震及其破裂特征[J]． 地震地質， 23(1)： 15--23．

Chai C Z， Liao Y H， Zhang W X， Xu W J， Shen X H， Tian Q J． 2001． Late Quaternary paleoearthquakes and their rupture features along the Lingwu fault[J]．SeismologyandGeology， 23(1)： 15--23 (in Chinese)．

陳立春． 2002． 河套斷陷帶的古地震、 強震復發(fā)規(guī)律和未來可能強震地點[D]． 北京： 中國地震局地質研究所： 49--67．

Chen L C． 2002．Paleoearthquakes，theLawofStrongEarthquakeRecurrenceandPotentialSitesfortheOccurrenceofFutureStrongEarthquakesintheHetaoFault-DepressionZone[D]． Beijing： Institute of Geology， China Earthquake Administration： 49--67 (in Chinese)．

陳立春， 冉勇康， 常增沛． 2003． 色爾騰山山前斷裂得令山以東段晚第四紀活動特征與古地震事件[J]． 地震地質， 25(4)： 555--565．

Chen L C， Ran Y K， Chang Z P． 2003． Characteristics of Late Quaternary faulting and paleoseismic events on the east of Delingshan segment of the Sertengshan piedmont fault[J]．SeismologyandGeology， 25(4)： 555--565 (in Chinese)．

鄧起東， 汪一鵬， 廖玉華， 張維岐， 李孟鑾． 1984． 斷層崖崩積楔及賀蘭山山前斷裂全新世活動歷史[J]． 科學通報, 29(9)： 557--560．

Deng Q D， Wang Y P， Liao Y H， Zhang W Q， Li M L． 1984． Colluvial wedges and Holocene activity along the rang-front fault of Helan Shan[J]．ChineseScienceBulletin, 29(9)： 557--560．

杜鵬， 柴熾章， 廖玉華， 宋方敏， 孟廣魁， 雷啟云， 王銀． 2009． 賀蘭山東麓斷裂南段套門溝—榆樹溝段全新世活動與古地震[J]． 地震地質, 31(2)： 256--264．

Du P， Chai C Z ， Liao Y H， Song F M， Meng G K， Lei Q Y， Wang Y． 2009． Study on Holocene activity of the south segment of the eastern piedmont fault of Helan mountains between Taomengou and Yushugou[J]．SeismologyandGeology, 31(2)： 256--264 (in Chinese)．

董治平， 何文貴， 戴華光． 2002． 包代河斷裂古地震初步研究[J]． 西北地震學報, 24(3)： 225--229．

Dong Z P， He W G， Dai H G． 2002． Preliminary research on paleoearthquakes of the Baodaihe fault[J]．NorthwesternSeismologicalJournal, 24(3)： 225--229 (in Chinese)．

郭星． 2014． 強震復發(fā)的隨機特征滑動模型及其應用方法研究[D]． 北京： 中國地震局地球物理研究所： 44--70．

Guo X． 2014．StochasticCharacteristic-SlipModelandItsApplicationMethod[D]． Beijing： Institute of Geophysics， China Earthquake Administration： 44--70 (in Chinese)．

郭星， 潘華． 2015． 強震復發(fā)間隔變異系數(shù)的一種計算方法[J]． 地震學報， 37(3)： 411--419．

Guo X， Pan H． 2015． A method for computing the aperiodicity parameter of the strong earthquake recurrence interval [J]．ActaSeismologicaSinica， 37(3)： 411--419 (in Chinese)．

江娃利， 侯治華， 肖振敏， 謝新生． 2000． 北京平原夏墊斷裂齊心莊探槽古地震事件分析[J]． 地震地質, 22(4)： 413--422．

Jiang W L， Hou Z H， Xiao Z M， Xie X S． 2000． Study on paleoearthquakes of Qixinzhuang trench at the Xiadian fault， Beijing plain [J]．SeismologyandGeology, 22(4)： 413--422 (in Chinese)．

江娃利， 謝新生， 王煥貞， 馮西英， 張康富． 2003． 山西大同盆地恒山北緣斷裂全新世古地震活動[J]． 中國地震, 19(1)： 8--19．

Jiang W L， Xie X S， Wang H Z， Feng X Y， Zhang K F． 2003． Holocene paleoseismic activities along the northern piedmont fault of Hengshan mountain， Datong Basin， Shanxi Province[J]．EarthquakeResearchinChina, 19(1)： 8--19(in Chinese)．

錢洪， 羅灼禮， 聞學澤． 1990． 鮮水河斷裂帶上特征地震的初步研究[J]． 地震學報, 12(1)： 22--29．

Qian H， Luo Z L， Wen X Z． 1990． Preliminary study of characteristic earthquakes in the Xianshuihe fault zone[J]．ActaSeismologicaSinica, 12(1)： 22--29 (in Chinese)．

冉勇康， 鄧起東， 楊曉平， 張晚霞， 李如成， 向宏發(fā)． 1997． 1679年三河—平谷8級地震發(fā)震斷層的古地震及其重復間隔[J]． 地震地質, 19(3): 193--201．

Ran Y K， Deng Q D， Yang X P， Zhang W X， Li R C， Xiang H F． 1997． Paleoearthquakes and recurrence interval on the seismogenic fault of 1679 Sanhe-PingguM8 earthquake， Hebei and Beijing[J]．SeismologyandGeology, 19(3): 193--201 (in Chinese)．

冉勇康， 鄧起東． 1998． 海原斷裂的古地震及特征地震破裂的分級性討論[J]． 第四紀研究， (3): 271--277．

Ran Y K， Deng Q D． 1998． Paleoearthquakes along Haiyuan fault and discussion of grading on rupture of large earthquakes[J]．QuaternarySciences， (3)： 271--277 (in Chinese)．

湯勇， 胡朝忠， 田勤儉， 王林， 楊攀新， 熊仁偉． 2014． 云南龍蟠—喬后斷裂劍川段古地震初步研究[J]． 地震, 34(3)： 117--124．

Tang Y， Hu C Z， Tian Q J， Wang L， Yang P X， Xiong R W． 2014． A preliminary study of paleo-earthquake in the Jianchuan section of Longpan-Qiaohou fault zone， Yunnan Province[J]．Earthquake， 34(3)： 117--124 (in Chinese)．

萬永魁， 沈軍， 于曉輝， 戴訓也， 王雷， 邵博， 李祥． 2014． 北京平原區(qū)夏墊活動斷裂滑動速率及古地震復發(fā)間隔[J]． 防災科技學院學報, 16(3)： 38--45．

Wan Y K， Shen J， Yu X H， Dai X Y， Wang L， Shao B， Li X． 2014． The slipping rates and ancient earthquake recurrence interval at Xiadian fault， Beijing plain[J].JournalofInstituteofDisasterPrevention, 16(3)： 38--45 (in Chinese)．

謝新生， 趙晉泉， 江娃利， 安衛(wèi)平， 李自紅， 張大衛(wèi)， 程新源， 高樹義， 王春和， 孫昌斌， 閆成國． 2007． 山西交城斷裂帶西張?zhí)讲廴率拦诺卣鹧芯縖J]． 地震地質, 29(4)： 744--755．

Xie X S， Zhao J Q， Jiang W L， An W P， Li Z H， Zhang D W， Cheng X Y， Gao S Y， Wang C H， Sun C B， Yan C G． 2007． Study on Holocene paleoearthquake in Xizhang trench on the Jiaocheng fault zone， Shanxi Province[J]．SeismologyandGeology, 29(4)： 744--755 (in Chinese)．

王華林． 1995． 1668年郯城8.5級地震斷裂的全新世滑動速率、 古地震和強震復發(fā)周期[J]． 西北地震學報， 17(4)： 1--12．

聞學澤． 1995． 活動斷裂地震潛勢的定量評估[M]． 北京： 地震出版社： 1--150．

Wen X Z． 1995．QuantitativeEstimatesofSeismicPotentialonActiveFaults[M]． Beijing： Seismological Press： 1--150 (in Chinese)．

楊明， 劉百篪． 2000． 時間相依的地震危險性概率分析[J]． 西北地震學報， 22(1)： 10--15．

Yang M， Liu B C． 2000． Probabilistic seismic hazard analysis depending on time[J]．NorthwesternSeismologicalJournal， 22(1)： 10--15 (in Chinese)．

張世民， 王丹丹， 劉旭東， 張國宏， 趙俊香， 羅明輝， 任俊杰， 王瑞， 張英禮． 2008． 北京南口—孫河斷裂晚第四紀古地震事件的鉆孔剖面對比與分析[J]． 中國科學: D輯， 38(7)： 881--895．

Zhang S M， Wang D D， Liu X D， Zhang G H， Zhao J X， Luo M H， Ren J J， Wang R， Zhang Y L． 2008． Using borehole core analysis to reveal Late Quaternary paleoearthquakes along the Nankou-Sunhe fault， Beijing[J]．ScienceinChina:SeriesD， 51(8)： 1154--1168．

Aki K． 1984． Asperities， barriers， characteristic earthquakes and strong motion prediction[J]．JGeophysRes， 89(B7)： 5867--5892．

Kanamori H． 1977． The energy release in great earthquakes[J]．JGeophysRes， 82(20): 2981--2987．

Papadimitriou E， Wen X Z， Karakostas V， Jin X S． 2004． Earthquake triggering along the Xianshuihe fault zone of western Sichuan， China[J]．PureApplGeophys， 161(8)： 1683--1707．

Reid H F． 1910．TheMechanicsoftheEarthquake，theCaliforniaEarthquakeofApril18， 1906[R]． Washington： State Earthquake Investigation Commission， Carnegie Institution of Washington： 43--47．

Schwartz D P， Coppersmith K J． 1984． Fault behavior and characteristic earthquakes： Examples from the Wasatch and San Andreas fault zones[J]．JGeophysRes， 89(B7)： 5681--5698．

Youngs R R， Coppersmith K J． 1985． Implications of fault slip rates and earthquake recurrence models to probabilistic seismic hazard estimates[J]．BullSeismolSocAm， 75(4)： 939--964．

The magnitude distribution law of characteristic earthquakes in Chinese mainland

Chen Zhongtian1)Guo Xing2),*Pan Hua1)Li Jinchen1)

1)InstituteofGeophysics，ChinaEarthquakeAdministration，Beijing100081，China2)NuclearandRadiationSafetyCenter，MinistryofEnvironmentProtection，Beijing100082，China

Based on the fundamental assumption that the magnitude distribution of characteristic earthquakes obeys the normal distribution， the quantification research are carried out on the magnitude distribution of characteristic earthquakes． Given sparse magnitude data of large earthquakes recurrence in situ， this paper uses magnitude data converted by coseismic displacement data， and each magnitude obeys identical distribution after normalization， then the standard deviationσwhich reflects the variability of the normalized magnitude can be got by statistic analysis． Using the approach proposed in this paper， analysis of coseismic displacements of 54 characteristic earthquakes in Chinese mainland suggests a generic value ofσ=0.08， which provides an important basis for further studying occurrence probabilistic model of characteristic earthquakes.

characteristic earthquake； magnitude； coseismic displacement； small sample

國家科技支撐項目(2012BAK15B01-08)資助.

2016-02-17收到初稿， 2016-06-15決定采用修改稿.

10.11939/jass.2016.06.009

P315.08

A

陳中天， 郭星， 潘華， 李金臣. 2016. 中國大陸特征地震震級分布規(guī)律. 地震學報, 38(6): 898--905. doi:10.11939/jass.2016.06.009.

Chen Z T, Guo X, Pan H, Li J C. 2016. The magnitude distribution law of characteristic earthquakes in Chinese mainland.ActaSeismologicaSinica, 38(6): 898--905. doi:10.11939/jass.2016.06.009.

*通訊作者 e-mail: guoxing@chinansc.cn