地震危險性分析中地震時空統計分布模型研究

徐偉進

（中國地震局地球物理研究所，北京 100081）

（作者電子信箱，徐偉進：wjxuwin＠163．com）

近幾十年來，隨著地震臺網建設的完善，世界各國均積累了大量的地震目錄資料。如何將地震目錄資料直接應用于地震危險性分析中是目前地震學研究的熱點問題之一。已有學者提出多種空間光滑方法將地震目錄資料直接應用于地震危險性分析中，然而，對于與地震危險性分析密切相關的地震時空統計分布的研究卻相對滯后，缺乏堅實的理論基礎。為了探討如何科學合理地將地震目錄資料直接應用于地震危險性分析中，我們根據中國大陸的地震目錄資料，在對地震目錄進行完整性分析和采用不同除叢方法進行余震刪除的基礎上，對地震發生的時間、空間統計分布特征做了系統、深入的研究，并定量研究了余震時空叢集對地震危險性分析的影響。

首先，我們采用定量統計方法，對我國儀器記錄地震和歷史破壞性地震的完整起始年限做了分析，得到了完整的地震目錄，以此作為研究地震時空統計分布特征的基礎資料。

對于時間分布特征研究，我們以地震區為研究單元，采用時空窗法和基于ETAS模型的隨機除叢法 兩種余震刪除方法對地震序列中的余震進行了刪除。使用χ2檢驗、K－S檢驗、游程檢驗對刪除余震后的主震目錄進行了一般時間獨立性檢驗，然后又進一步使用擴散熵分析及時間K－function法對主震目錄的時間相關性和時間叢集性進行了嚴格的統計檢驗，結果表明，我國大陸地震事件的發生具有時間獨立性，即地震發生在時間上遵從泊松過程，且是典型的非平穩泊松過程。

在地震空間分布特征的研究上，我們采用 “Cellular Seismology”法并輔以Wilcoxon－test和Spearman Rank Correlation test等統計檢驗方法，對 “未來地震發生在已發生地震附近”這一假設進行了驗證分析。通過對我國8個地震區和22個地震帶的檢驗研究，結果表明“未來地震發生在已發生地震附近”是一個真實、可測量的物理現象，即空間分布的地震活動性資料可以作為未來地震潛在震源區的界定，這也意味著我國地震空間分布呈現叢集模式。進一步地，我們使用K－function空間點模式分析工具，以Neyman－Scott空間叢集模型為基本模型，分析了中國大陸地區地震點空間分布模型，結果表明，我國地震空間分布符合廣義Thomas模型。

最后，我們定量研究了余震時空叢集對地震危險性分析的影響。我們采用蒙特卡洛方法，以ETAS模型為地震活動性模型模擬地震目錄，以基于空間光滑地震活動性模型的地震危險性分析方法計算地震峰值加速度，分析了余震對地震危險性分析結果的影響。結果表明，余震對50年超越概率10%水平下峰值加速度的影響均值為6%左右，最大為18%左右，并且超越概率水平越高，余震影響越大。此外，文中我們還提出了一種在地震危險性分析中考慮“強余震產生的地震動值超過主震產生的地震動值”對地震危險性影響的計算方法。計算結果表明，強余震對某些場點的地震危險性分析結果的影響可能達到60%以上。

本文的研究可以為我國直接采用地震目錄資料進行地震危險性分析提供理論基礎；進一步加深對我國地震時空分布特征的認識和了解；明晰余震時空叢集對地震危險性分析結果的影響。