張涵博
巨災保險是巨災風險管理中的重要工具。通過巨災保險,可以將自然災害的損失風險進行有效的分散和轉(zhuǎn)移,促進防災減災和災后重建工作的開展,對于一國的經(jīng)濟社會繁榮穩(wěn)定有著非常重要的作用。我國是一個自然災害頻發(fā)的國家,面臨著除了火山爆發(fā)之外所有的巨災風險,巨災保險對我國有著尤為重要的意義。
巨災風險有著與常規(guī)風險截然不同的特點:從發(fā)生頻率來看,常規(guī)風險通常發(fā)生頻率較高,但巨災風險的發(fā)生頻率極低;從損失規(guī)模來看,常規(guī)風險的損失一般只與單個風險標的的價值有關,而在巨災風險中通常某一地區(qū)的所有風險標的都受影響;從空間和時間分布來看,常規(guī)風險的發(fā)生和分布比較分散,但巨災風險通常會在某一時間、某一地區(qū)集中發(fā)生,造成高額損失。上述特點使得巨災保險相對于普通保險在產(chǎn)品開發(fā)和定價方面有著更大的困難:常規(guī)風險的頻發(fā)性使得采用大數(shù)定律等傳統(tǒng)保險精算方法估計風險損失的波動性較小,而巨災風險的低頻高危性使得傳統(tǒng)方法對風險損失的估計與實際情況可能會有較大出入。要衡量巨災風險的大小,必須通過巨災建模(Catastrophe Modeling),在經(jīng)驗數(shù)據(jù)的基礎上構(gòu)建巨災損失模型。
災損曲線是反映巨災的損失程度與發(fā)生頻率之間關系的曲線,是巨災損失模型的核心組成部分。它可以確定某一損失規(guī)模巨災事件的發(fā)生頻率,從而求出巨災總損失的數(shù)學期望值并實現(xiàn)巨災保險的定價。超概率曲線(Exceedance Probability Curve,EP)是災損曲線的主要形式之一。本文將首先介紹超概率曲線構(gòu)建的一般原理,進而提出一種在經(jīng)驗數(shù)據(jù)的基礎上通過計算機模擬構(gòu)建巨災超概率曲線的方法,最后將探究超概率曲線在巨災保險中的應用及意義。
超概率曲線是概率論中進行累積頻率分析(Cumulative Frequency Analysis)的重要工具,它是用圖形表示在給定時期內(nèi)或給定事件中超過一定損失水平事件的發(fā)生概率,通常用縱軸表示概率,橫軸表示損失水平。對于某種災害損失,超概率是指在某段時間內(nèi)或某些災害事件中,發(fā)生大于特定損失規(guī)模災害事件的概率。當損失為0時,超過這一損失的概率為1;隨著損失規(guī)模的提高,超概率將隨之下降。
一般而言,在構(gòu)建超概率曲線時,服從以下一些基本假設:
1.有一系列自然災害損失事件Ei;
2.對于每起事件,都有發(fā)生概率pi和對應的損失Li;
3.某一年中災害發(fā)生的次數(shù)并不局限于一次,可以是多次;
4.所有災害事件都是相互獨立的,服從 Bernoulli分布,即對于任意的m和n,都有P(Em/En)=P(Em),每起災害損失事件的發(fā)生時間、發(fā)生頻率、損失程度等參數(shù)都不受其他災害損失事件發(fā)生與否的影響。
在以上假設下,每一起災害損失事件Ei都是一個獨立隨機變量,其概率分布函數(shù)定義如下。
P(Ei發(fā)生)=pi
P(Ei不發(fā)生)=(1-pi)
若災害不發(fā)生,則損失為零。所以某年中某一災害事件Ei的期望損失為:
如果我們用AAL(Average Annual Loss)來表示所有災害事件在一年中的總期望損失,則其表達式為:
假設某一年中只發(fā)生了一次災害事件,則某一損失水平的超概率EP(Li),即超過該損失水平的災害損失事件的發(fā)生概率,就相當于1減去這一年中所有低于該損失水平的災害事件都沒有發(fā)生的概率。它可以由以下式子求得:
將相應的損失額度Li和相應的EP(Li)繪制在圖上,即可得到這一年中巨災損失的超概率曲線,如下圖:
圖:超概率曲線示例圖
在圖中,橫軸Li表示損失水平(橫軸上的數(shù)字是假設的損失值),縱軸EP(Li)表示其所對應的超概率。在損失水平為0時,所有損失事件的損失額度都會超過0,此時對應的超概率是1。一般而言,隨著損失水平的提高,相應的超概率也越來越低,但其下降幅度將逐步減緩。即超概率曲線滿足以下性質(zhì):
在以上基本原理中,若要構(gòu)建年度災害損失超概率曲線并求得AAL,必須能夠確定在某一年中,某一特定的災害損失事件Ei發(fā)生的確切概率 pi。理論上,pi可以通過統(tǒng)計多年的經(jīng)驗數(shù)據(jù)獲得。但在實際中,一方面,由于我國巨災監(jiān)控體系是在建國之后才逐步建立和完善的,歷史數(shù)據(jù)并不全面;另一方面,改革開放后經(jīng)濟社會快速發(fā)展使得經(jīng)濟體制轉(zhuǎn)型前后的災害損失數(shù)據(jù)缺乏可比性,因此中國巨災損失歷史數(shù)據(jù)非常有限,通過簡單統(tǒng)計這些數(shù)據(jù)得到的pi并不準確。在構(gòu)建超概率曲線時應當借助計算機生成虛擬的災害事件來模擬災害的損失情況,以彌補歷史數(shù)據(jù)的不足。本部分將提出一種在經(jīng)驗數(shù)據(jù)的基礎上通過計算機模擬構(gòu)建年度巨災損失超概率曲線的方法。首先根據(jù)歷史數(shù)據(jù)得到單起災害事件中損失的分布情況,然后建立每年某種巨災事件發(fā)生次數(shù)的分布,最后由計算機模擬年度災害損失情況以構(gòu)建年度巨災損失超概率曲線。
第一步:建立單起災害事件損失分布
一般而言,對于地震、臺風這樣的自然災害,我國每年都會發(fā)生數(shù)十次,在二十年的時間里能夠積累到數(shù)百個有效的災害損失數(shù)據(jù)。如果這些災害事件的損失都服從同樣的分布,數(shù)百起事件已經(jīng)足以構(gòu)成大樣本,通過統(tǒng)計這些數(shù)據(jù)就能找到單起災害事件的損失分布。在第一步中,我們把某種自然災害的總體看成一個隨機變量,把所收集的該種災害損失事件看作是巨災損失總體的樣本。在構(gòu)建單起災害事件損失分布時接受以下假設:
1.收集的每起該種災害損失事件Ei都是該種災害損失隨機變量ξ的一個樣本,且樣本容量足夠大;
2.所有災害損失事件Ei都服從同樣的分布;
3.所有災害損失事件Ei都是相互獨立的;
4.在每一起災害損失事件中,對于某一可能的損失水平Li,有對應的損失概率pi。
在以上假設下,某一時間段內(nèi)我國所有的該種災害損失事件相當于在同樣條件下重復進行的數(shù)百次實驗,構(gòu)成了一個獨立試驗序列概型。由于每次試驗的結(jié)果Ei與其他各次試驗結(jié)果無關,因此這一系列重復實驗構(gòu)成了n重Bernoulli實驗。
由于假設所有災害損失是一個隨機變量總體ξ,可以設(x1,x2,…xn,)是總體 ξ的樣本 Ei的觀察值,即每起災害損失事件中的經(jīng)濟損失數(shù)額。將它們按大小排列為,令
則Fn(x)的圖形就構(gòu)成了累積頻率曲線。它是跳躍式上升的一條階梯曲線。若觀測值不重復,則每一躍度為1/n;若觀測值有重復,則按照1/n的倍數(shù)跳躍上升。
根據(jù)Bernoulli大數(shù)定律,在獨立試驗序列中,當試驗次數(shù)n無限增加時,事件A的頻率k/n(k是n次試驗中事件A發(fā)生的次數(shù)),收斂于它的概率P(A)=p。即,對于任意給定的ε>0,有
也就是說,當試驗在不變的條件下重復進行很多次時,隨機事件的頻率在它的概率附近擺動,并將隨著實驗次數(shù)的增加而不斷趨近其概率。在樣本的累積頻率曲線Fn(x)中,對于任何實數(shù)x,F(xiàn)n(x)等于樣本的n個觀察值中不超過x的個數(shù)除以樣本容量n,因此Fn(x)可以作為未知的整體ξ的分布函數(shù)Fξ(x)的一個近似。樣本容量n越大,近似得越好。在數(shù)百起災害損失事件的大樣本下,可以近似地認為樣本累積頻率曲線Fn(x)等同于災害損失整體分布曲線Fξ(x)。由于Fn(x)實質(zhì)上是所有低于某一損失水平x的災害事件發(fā)生的概率,因此用1-Fn(x)可得出所有高于x的災害損失事件的發(fā)生概率,即單起災害損失的超概率曲線。
需要說明的是,并不是所有的災害損失數(shù)據(jù)都能夠符合前提假設:
第一,對于最重要的“同分布”假設,災害損失的原始數(shù)據(jù)很可能無法滿足:隨著經(jīng)濟社會的發(fā)展,一個地區(qū)的財富積累會越來越多。在其他條件不變的情況下,同樣破壞程度的災害事件造成的損失將逐步擴大。再加上通貨膨脹因素,嚴格來說,在長時間跨度中,災害損失并不符合同分布假設。在這種情況下應對損失數(shù)據(jù)進行一定的處理,以消除通貨膨脹或經(jīng)濟增長等外部因素的影響。
第二,自然災害的特點可能使得“獨立”這一假設在現(xiàn)實中不一定成立。同一類型的不同災害事件之間,以及不同類型的不同災害事件之間,的確可能存在一定的因果關系,并不一定是完全獨立的。由于科學界對自然災害的形成和相互作用機理仍未充分研究,因此只能暫時假定所有災害事件都是相互獨立的。
第三,某些自然災害的發(fā)生頻率較低,難以滿足“樣本容量足夠大”的假設。在我國,海嘯、火山爆發(fā)發(fā)生的可能性微乎其微,由于樣本容量太小,無法通過上述方法來找出這兩種自然災害的損失分布情況。
此外,對于一些極端損失事件,在必要時應當作為極端值舍棄。如汶川地震,其損失占2008年當年地震直接經(jīng)濟損失的99.16%,占 1990~2008年中國地震總損失的90%以上,遠超出其他地震的損失規(guī)模(鄭通彥,2009)。極端大的損失對應著非常小的概率,為了維持數(shù)據(jù)的可處理性,根據(jù)小概率事件的實際不可能性原理,在構(gòu)建地震的單起災害損失分布時應當去除汶川地震這樣的極端事件。
經(jīng)過以上數(shù)據(jù)處理后可建立單起災害事件損失分布模型。通過它可以得到在每起災害事件中,損失水平超過Li的概率EP(Li),由此構(gòu)建單起災害損失事件超概率曲線。
第二步:構(gòu)建年度災害發(fā)生次數(shù)分布模型
與構(gòu)建單起災害損失分布模型的假設類似,我們假定每年某種災害的發(fā)生次數(shù)也是相互獨立并服從同一分布的。一般而言,每年某種自然災害(如地震、臺風、洪水等)的發(fā)生次數(shù)往往呈現(xiàn)出圍繞著平均值上下波動的現(xiàn)象。如果有相當長的年度災害發(fā)生次數(shù)數(shù)據(jù)(如20年左右),可以驗證其是否服從正態(tài)分布。
首先,繪出年度災害發(fā)生次數(shù)數(shù)據(jù)的分布直方圖,并求出其峰度與偏度。如果直方圖呈現(xiàn)出正態(tài)分布的形態(tài),并且這些數(shù)據(jù)的峰度和偏度趨近正態(tài)分布的特征(峰度為3,偏度為0),則可以近似地認為年度災害發(fā)生次數(shù)數(shù)據(jù)基本符合正態(tài)分布。
為了進一步檢驗樣本數(shù)據(jù)是否服從正態(tài)分布,可以通過統(tǒng)計軟件來進行相關檢驗:通過Eviews軟件可以求出樣本數(shù)據(jù)的Jarque-Bera統(tǒng)計量,通過 SAS或SPSS軟件也可以求出樣本數(shù)據(jù)的Shapiro-Wilk或Kolmogorov-Smirnov統(tǒng)計量。如果統(tǒng)計量的值大于0.05,則接受檢驗原假設,即樣本數(shù)據(jù)服從正態(tài)分布。
若樣本數(shù)據(jù)的期望為μ,標準差為σ,樣本數(shù)據(jù)具有較大的樣本容量,符合正態(tài)分布的特征,且能夠通過正態(tài)分布的相關檢驗,則可近似地認為年度災害發(fā)生次數(shù)分布服從期望為μ、標準差為σ的正態(tài)分布。
第三步:通過計算機模擬構(gòu)建年度災害損失超概率曲線
在前面兩步中,我們分別建立了單起災害損失分布和年度災害發(fā)生次數(shù)分布。接下來,通過計算機隨機抽取每年災害發(fā)生次數(shù)和每次災害損失,經(jīng)過大量的抽取模擬后(如1000年),可以得到長時間段的年度災害損失狀況,以此可構(gòu)建年度災害損失分布曲線。
在構(gòu)建某種年度災害損失超概率曲線時,我們假設:1.任意一起災害的損失都服從第一步中得到的單起災害事件損失分布;2.任意一年中某種災害發(fā)生次數(shù)都服從第二步中得到的年度災害發(fā)生次數(shù)分布,即期望為μ、標準差為σ的正態(tài)分布。
在以上假設下,采用以下算法模擬某種災害年度損失總額:
步驟一,根據(jù)年度災害發(fā)生次數(shù)分布生成第i年災害發(fā)生次數(shù)ni;
步驟二,根據(jù)單起災害事件損失分布生成ni次災害損失事件,每次分別有kij的損失,構(gòu)成第i年中的災害損失序列;
步驟三,將這ni次災害事件的損失加總,構(gòu)成第i年的災害損失總額
步驟四,將步驟一至步驟三重復m次,得到m年中每年的年度災害損失總額序列:K1,K2,…,Ki…,Km-1,Km。
最后,采用第一步中的方法,使用這m個模擬的年度災害損失總額構(gòu)建它們的分布模型,就得到了年度災害損失超概率曲線。在計算機的輔助下可以模擬出數(shù)千年的年度災害損失情況,由此得到的超概率曲線比簡單通過幾十年的災害損失記錄得到的超概率曲線更加準確。
以超概率曲線為代表的災損曲線確立了災害損失的大小與發(fā)生概率之間的關系,是巨災保險費率厘定、風險評估、承保決策及風險轉(zhuǎn)移過程中必不可少的重要工具。在前文構(gòu)建的超概率曲線的基礎上,本部分將探究災損曲線在巨災保險中的應用。
第一,災損曲線為巨災保險的費率厘定提供了依據(jù)。在厘定保險的純保費時,需要滿足的基本條件是收取保費的現(xiàn)值應等同于損失賠付期望的現(xiàn)值。由于災害事件發(fā)生頻率低,通過歷史經(jīng)驗數(shù)據(jù)很難準確估計出未來損失的期望。以美國1994年北嶺地震為例,在地震前美國地震保險的低費率(每千美元保額的保費約為2美元)導致北嶺地震造成的125億美元巨額保險索賠,超過了此前20年間美國保險公司所收取的地震保險保費總和(朱文杰,2006),因此合理的費率對巨災保險至關重要。在建立合適的災損曲線后,保險公司能夠掌握災害損失大小與發(fā)生概率之間的關系,通過計算出年度總期望損失(AAL)而得到的巨災保險費率將更加符合實際情況。
第二,災損曲線使保險公司能夠較準確地評估其所承保的巨災風險的規(guī)模。保險公司對于災害賠付的承受能力是有一定限度的,為了確保公司經(jīng)營的穩(wěn)定性和連續(xù)性,對承保風險的評估非常必要。根據(jù)超概率曲線可以得到一定損失規(guī)模災害事件的重現(xiàn)周期,從而評估巨災風險的大小。在超概率曲線示例圖中,損失為100的災害事件對應的超概率約為0.052,這意味著損失規(guī)模超過100的災害事件的重現(xiàn)周期是1÷0.052≈19.2 年。如果一家保險公司設定其所能接受的最大風險是重現(xiàn)期為50年的災害事件,即 EP=0.02,那么根據(jù)超概率曲線示例圖,在該巨災保險中最大可能損失約為132。
第三,災損曲線可以為保險公司巨災保險的承保決策提供依據(jù)。為了確保安全穩(wěn)健經(jīng)營,保險公司會盡量避免出現(xiàn)償付能力不足的情況。假如一家保險公司可以接受償付能力不足情況的最大出現(xiàn)概率是p0,用L代表賠付損失,n代表承保保單的數(shù)量,a代表每份保單收取的保費,S代表現(xiàn)有盈余,那么保險公司經(jīng)營的約束條件是:
在超概率曲線中,根據(jù)EP=p0可以求出對應的L。若保險公司的保費收入與盈余之和無法達到L,那么此時就無法滿足以上約束條件,繼續(xù)承保該巨災保險會使公司面臨償付能力不足的風險。
第四,災損曲線可用于巨災風險轉(zhuǎn)移決策。巨災風險具有發(fā)生頻率低、損失規(guī)模大、較易集中發(fā)生的特點,單一一家保險公司往往無法單獨承受巨災保險的巨額賠付。承保巨災保險的保險人可選擇再保險或巨災債券等風險轉(zhuǎn)移措施,以進一步分散巨災風險并減少需要持有的自有資本,但需支付再保費或債券利息而降低保險公司的預期收益。為了在保證償付能力的前提下實現(xiàn)收益最大化,保險公司必須維持自留風險和分出風險的平衡。在評估巨災風險轉(zhuǎn)移措施時,災損曲線具有以下作用:一方面,通過災損曲線,原保險人和再保險人可以計算自留風險和分出風險的期望賠付,進而確定再保險的保費或巨災債券的利息;另一方面,保險人可以通過災損曲線計算出在滿足償付能力約束條件時自身所能承受的最大風險,進而決定分出額的大小。因此,災損曲線可以幫助保險公司做出合適的巨災風險轉(zhuǎn)移決策。
總體而言,災損曲線作為巨災模型的核心組成部分,是評估巨災風險規(guī)模的重要工具,對巨災保險的產(chǎn)品開發(fā)、運營管理有著決定性的作用。應用計算機模擬使得災損曲線能夠擺脫歷史數(shù)據(jù)不足帶來的困擾,增強了災損曲線的準確性,并使得巨災保險的定價和保險公司的決策更加合理。