錢(qián)小仕 王福昌 盛書(shū)中
(防災(zāi)科技學(xué)院,河北 065201)
大量的地震實(shí)際資料分析表明,地震震級(jí)-頻度關(guān)系并非總是滿足G-R公式,在許多情況下震級(jí)高端或低端出現(xiàn)明顯的 “掉頭”和“擺尾”現(xiàn)象,說(shuō)明震級(jí)分布為負(fù)指數(shù)分布僅在有限震級(jí)區(qū)時(shí)成立。實(shí)際資料特別是強(qiáng)震數(shù)據(jù)對(duì)G-R關(guān)系的偏離,一定程度上影響了利用該關(guān)系式推測(cè)未來(lái)強(qiáng)震發(fā)生危險(xiǎn)性的可信度。極值理論表明對(duì)充分大的閾值,隨機(jī)變量超過(guò)閾值的超出量的極限分布為廣義帕累托(GPD)分布,為研究強(qiáng)震震級(jí)分布提供了統(tǒng)計(jì)理論依據(jù)?;趶V義帕累托分布給出若干地震活動(dòng)性參數(shù)的估計(jì)方法,包括強(qiáng)震震級(jí)分布、地震平平復(fù)發(fā)周期、重現(xiàn)水平和期望重現(xiàn)震級(jí)、地震危險(xiǎn)性概率和潛在震級(jí)上限估計(jì)等,并對(duì)云南地區(qū)歷史地震資料進(jìn)行分析,研究該地區(qū)強(qiáng)震震級(jí)分布特征。
設(shè)X1,X,2,…,Xn為地震震級(jí)隨機(jī)變量列,假設(shè)它們相互獨(dú)立且服從同一分布F(x),Balkema,de Haan(1974)和 Pickands(1975)指出,若最大震級(jí)近似服從廣義極值分布則對(duì)充分大的震級(jí)閾值u,震級(jí)超出量X-u近似服從廣義帕累托分布且兩者具有相同的形狀參數(shù)ξ。
若假設(shè)最大地震震級(jí)分布為廣義極值分布,根據(jù)上述理論,可以推得若干個(gè)基于廣義帕累托分布的地震活動(dòng)性參數(shù)估計(jì)。
(1)地震震級(jí)分布的估計(jì):
其中Nu為n個(gè)地震中超過(guò)給定發(fā)震震級(jí)閾值u的次數(shù),?()F y為3參數(shù),,hξ σ′的廣義帕累托分布,其中
(2)平平復(fù)發(fā)周期估計(jì): 假設(shè)地震數(shù)據(jù)為日觀測(cè)數(shù)據(jù),一年按365天計(jì)算,若震級(jí)分布為3參數(shù)廣義帕累托分布G(x;u,σ~,ξ),則發(fā)生震級(jí)為x的地震平平復(fù)發(fā)周期為
(3)重現(xiàn)水平估計(jì):對(duì)給定概率p,不超過(guò)概率p的發(fā)震震級(jí)(即p-分平數(shù))為
若假設(shè)地震數(shù)據(jù)為日觀測(cè)數(shù)據(jù),則T年重現(xiàn)期的重現(xiàn)水平相當(dāng)于的分平數(shù)。
(4)地震危險(xiǎn)性概率: 未來(lái)T年發(fā)震震級(jí)超過(guò)x的概率
(5)潛在震級(jí)上限估計(jì):當(dāng)形狀參數(shù)ξ<0時(shí),震級(jí)分布有上限
以云南地區(qū)(21o~29oN,97o~106oE)1930-2010年M≥2.0地震目錄為基礎(chǔ)資料,進(jìn)行實(shí)證分析。通過(guò)觀察樣本平平超出量函數(shù)是否呈直擬和不同閾值選取下參數(shù)估計(jì)的穩(wěn)定狀況,選取閾值u=5.0。采用極大似然法估計(jì)廣義帕累托分布的參數(shù),得和的標(biāo)準(zhǔn)差分別為0.0329和0.0603,于是可得ξ的95%置信區(qū)時(shí)為[-0.351,-0.222]。從ξ的95%置信區(qū)時(shí)可以看出,ξ為負(fù)值,所以云南地區(qū)震級(jí)有上限,上限估計(jì)為Mmax=8.3。對(duì)比云南地區(qū)所給地震目錄資料范圍內(nèi)的歷史最大震級(jí)(8.0級(jí)),說(shuō)明廣義帕累托分布不僅可以作為確定潛在震源區(qū)最大震級(jí)上限的一種有效方法。此外,在工程抗震中,潛在震級(jí)上限的確定總是與特定建筑結(jié)構(gòu)的使用年限相聯(lián)系,即給定年限內(nèi)可能出現(xiàn)的最大地震,這相當(dāng)于要估計(jì)高分平數(shù)。模擬計(jì)算表明,高分平數(shù)估計(jì)在一定范圍內(nèi)受閾值選取的影響較小,表現(xiàn)穩(wěn)定。另一方面在實(shí)際應(yīng)用中還可以根據(jù)建筑結(jié)構(gòu)的抗震要求選擇相應(yīng)的高分平數(shù)作為上限估計(jì),因此利用高分平數(shù)估計(jì)潛在震級(jí)上限不僅具備使用的靈活性而且還具有一定的估計(jì)穩(wěn)定性。