[發明專利]一種基于拉丁超立方抽樣方法進行巨災事件模擬的方法及計算機程序產品

在審

• 申請號: 202110390960.3
• 申請日: 2021-04-12
• 公開（公告）號: CN113344711A
• 公開（公告）日: 2021-09-03
• 發明（設計）人: 楊浩;袁曦;魯涵;賈凱
• 申請（專利權）人: 中國再保險（集團）股份有限公司
• 主分類號: G06Q40/08
• 分類號: G06Q40/08;G06F17/18
• 代理公司: 北京卓愛普專利代理事務所(特殊普通合伙) 11920
• 代理人: 王玉松

1.一種基于拉丁超立方抽樣方法進行巨災事件模擬的方法，其特征在于，包括模擬巨災事件序列的方法的如下步驟：

獲取參數步驟 獲取被模擬巨災事件的事件頻率和事件標識；

排序步驟 將被模擬巨災事件按照事件頻率進行排序，得到相對應的排序事件頻率和排序事件標識；

確定抽樣間隔值步驟 用排序事件頻率總和除以抽樣數N得到每次的抽樣間隔值δ；

抽樣步驟 對所述排序事件頻率按照δ的間隔依次進行抽樣，記錄每次抽樣對應的排序事件標識；

亂序步驟 將全部抽取的排序事件標識進行亂序得到模擬巨災事件序列。

2.根據權利要求1所述的方法，其特征在于，

所述獲取參數步驟還包括：獲取所述排序事件頻率的全部累積頻率WT；

所述抽樣步驟包括：

當前排序事件頻率位于與其相鄰的兩個所述累積頻率WTn之間時，記錄對應的排序事件標識，直至抽得N個所述排序事件標識。

3.根據權利要求2所述的方法，其特征在于，第一次抽樣時，在大于0并小于δ之間的范圍內任取一個事件頻率記錄為所述當前排序事件頻率初始值。

4.根據權利要求1所述的方法，其特征在于，還包括模擬巨災事件發生時間序列的方法的如下步驟：

獲取參數步驟 獲取所述被模擬巨災事件的時間發生頻率和時間標識；

排序步驟 將所述被模擬巨災事件按照所述時間發生頻率排序，得到相應的排序時間發生頻率和排序時間標識；

確定抽樣間隔值步驟 用排序時間發生頻率總和除以抽樣數N得到每次的抽樣間隔值δ_時間；

抽樣步驟 對排序后的被模擬巨災事件按照δ_時間的間隔依次進行抽樣，記錄N個每次抽樣對應的排序時間標識；

亂序步驟 將全部抽取的排序時間標識進行亂序得到模擬巨災事件發生時間序列。

5.根據權利要求4所述的方法，其特征在于，

所述獲取參數步驟還包括：獲取所述排序時間發生頻率的全部累積頻率WT_時間；

所述抽樣步驟包括：

當前排序時間發生頻率位于與其相鄰的兩個所述累積頻率WT_時間n之間時，記錄對應的排序時間標識，直至抽得N個所述排序時間標識。

6.根據權利要求1-5任一項所述的方法，其特征在于，計算所述抽樣數包括如下步驟：

取值步驟 獲取模擬年數Y；

計算概率步驟 按照式I依次計算P(k)，

其中，P(k)表示一年發生k次巨災事件的概率，k表示發生巨災事件的次數，k＝0，k＝1,k＝2,……，λ為排序事件頻率總和，e為自然常數；

計算年數步驟 按照式II依次計算Q(k)，

Q(k)＝P(k)×Y 式II

其中，Q(k)表示發生k次巨災事件的年數；k表示發生巨災事件的次數，k＝1,k＝2,……，Y表示模擬年數；

重復計算年數步驟 將每次計算得到的Q(k)與臨界值比較，當Q(k)大于等于臨界值小于1，且小于λ×Y，令Q(k)＝1，當Q(k)小于所述臨界值時，令Q(k)＝0，計算年數步驟停止；

記錄所述k為m；

求和步驟 按照式III計算N，

其中，N表示抽樣數，k表示發生巨災事件的次數，k＝1,……，m，Q(k)表示發生k次巨災事件的年數。

7.根據權利要求6所述的方法，其特征在于，所述模擬巨災事件發生時間序列的方法還包括時間換算步驟：

根據Q(k)得到k在所述Y年的次數分布；

將所述次數分布進行亂序得到年發生次數序列；

將所述年發生次數序列按照每次事件的年標識序列表示；

根據每次事件的年標識將被模擬巨災事件的時間標識按對應的年標記為換算時間標識。