[發(fā)明專利]一種基于地震儀實(shí)測(cè)及數(shù)值模擬近實(shí)時(shí)烈度圖生成方法

在審

• 申請(qǐng)?zhí)枺?/td>: 201911408975.7
• 申請(qǐng)日: 2019-12-31
• 公開（公告）號(hào): CN111060962A
• 公開（公告）日: 2020-04-24
• 發(fā)明（設(shè)計(jì)）人: 周立;陳志高;張偉;范濤
• 申請(qǐng)（專利權(quán)）人: 武漢地震科學(xué)儀器研究院有限公司;中國地震局地震研究所;南方科技大學(xué)
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28;G01V1/30;G01V1/34
• 代理公司: 北京棘龍知識(shí)產(chǎn)權(quán)代理有限公司 11740
• 代理人: 謝靜

1.一種基于地震儀實(shí)測(cè)及數(shù)值模擬近實(shí)時(shí)烈度圖生成方法，其特征在于包括以下步驟：

1)預(yù)先設(shè)置好地球介質(zhì)模型；

2)通過地震儀獲取地震數(shù)據(jù)；

3)通過將獲取的地震數(shù)據(jù)經(jīng)計(jì)算得到儀器地震烈度值；

4)通過地震預(yù)警或速報(bào)快速在震后數(shù)秒或數(shù)十秒獲取震源基本參數(shù)；

5)通過1)和4)運(yùn)用數(shù)值模擬方法快速計(jì)算震區(qū)附近強(qiáng)地面運(yùn)動(dòng)信息得到該區(qū)域的理論烈度值；

6)通過處理3)地震儀器實(shí)測(cè)的儀器烈度值和5)數(shù)值模擬計(jì)算得到的理論烈度值，當(dāng)理論烈度和儀器烈度存在差異時(shí)，用儀器烈度校正理論烈度；

7)在得到震中區(qū)各網(wǎng)格點(diǎn)的烈度值后，可繪制出烈度等值線圖，其主要根據(jù)地震行業(yè)規(guī)范，對(duì)插值后得到的柵格圖進(jìn)行配色處理，得到最終的地震烈度圖。

2.根據(jù)權(quán)利要求1所述一種基于地震儀實(shí)測(cè)及數(shù)值模擬近實(shí)時(shí)烈度圖生成方法，其特征在于，在步驟3中地震儀烈度值的計(jì)算方法如下：

3-1)首先通過對(duì)地震的研究，統(tǒng)計(jì)得到儀器烈度I_i與峰值地面加速度PGA以及烈度與峰值地面速度PGV的經(jīng)驗(yàn)關(guān)系：

當(dāng)Ⅴ≤I_i≤Ⅷ時(shí)，I_i與峰值地面加速度PGA的經(jīng)驗(yàn)關(guān)系公式：

I_i＝a·log(PGA-b)

其中a，b為擬合得到的常數(shù)，我們擬合的全球數(shù)據(jù)得到a＝3.659；b＝1.662；在特定區(qū)域時(shí)，a，b值會(huì)隨本地的介質(zhì)不同有所變化；

當(dāng)Ⅷ＜I_i≤Ⅻ時(shí)，I_i與峰值地面速度(PGV)的經(jīng)驗(yàn)關(guān)系公式：

I_i＝a·log(PGV)+b

其中a，b為擬合得到了常數(shù)，我們擬合的全球數(shù)據(jù)得到a＝3.466；b＝2.351；在特定區(qū)域時(shí)，a，b值會(huì)隨本地的介質(zhì)不同有所變化；

3-2)根據(jù)統(tǒng)計(jì)儀器烈度I_i與峰值地面加速度PGA以及烈度與峰值地面速度PGV的經(jīng)驗(yàn)關(guān)系公式，最終獲得地震的儀烈度值。

3.根據(jù)權(quán)利要求1或2所述一種基于地震儀實(shí)測(cè)及數(shù)值模擬近實(shí)時(shí)烈度圖生成方法，其特征在于，在步驟4中所述的震源基本參數(shù)包含震源位置、發(fā)震時(shí)刻、震級(jí)、斷層面解，其中所述的震源位置包含經(jīng)緯度以及深度值。