Cramér–Lundberg 模型
让我从最简单的保险风险模型开始,为您系统性地讲解 Cramér–Lundberg 模型。这个模型是经典风险理论的核心,用来描述保险公司盈余随时间变化的过程。
第一步:模型的基本组成部分
考虑一个保险公司,其盈余过程 \(U(t)\) 由三个基本要素构成:
- 初始资本:用 \(u\) 表示,\(u \geq 0\),是公司在 \(t=0\) 时刻的初始准备金。
- 保费收入:假设保费以常数速率 \(c > 0\) 连续流入。这意味着在时间区间 \((0, t]\) 内,保费总收入为 \(ct\)。
- 索赔过程:索赔的发生是随机的。假设索赔按照泊松过程到达。具体来说,在时间区间 \((0, t]\) 内发生的索赔次数 \(N(t)\) 服从参数为 \(\lambda t\) 的泊松分布,即 \(P(N(t)=k) = e^{-\lambda t} (\lambda t)^k / k!\),其中 \(\lambda > 0\) 是单位时间的平均索赔次数。
每一个索赔的金额也是一个随机变量。设第 \(i\) 次索赔的金额为 \(Y_i\),我们假设 \(\{Y_i\}\) 是独立同分布的随机变量序列,其共同的分布函数为 \(F(y)\),且 \(F(0)=0\),期望值 \(\mathbb{E}[Y_1] = \mu < \infty\)。此外,索赔额序列 \(\{Y_i\}\) 与索赔发生过程 \(N(t)\) 相互独立。
第二步:盈余过程的数学定义
综合以上三个要素,保险公司在时刻 \(t\) 的盈余 \(U(t)\) 可以定义为:
\[U(t) = u + ct - \sum_{i=1}^{N(t)} Y_i \]
其中:
- \(u\) 是初始资本。
- \(ct\) 是到时刻 \(t\) 为止累计收取的保费。
- \(\sum_{i=1}^{N(t)} Y_i\) 是到时刻 \(t\) 为止累计支付的索赔总额。
这个过程 \(\{U(t), t \geq 0\}\) 就称为 Cramér–Lundberg 模型(或称经典复合泊松风险模型)。它是一个随机过程,其样本路径(即一次具体的实现)是初始值为 \(u\) 的直线(斜率为 \(c\)),然后在随机的时间点(索赔发生时刻)向下跳跃,跳跃的幅度就是随机的索赔额 \(Y_i\)。
第三步:安全负载条件与破产概率
保险公司运营的核心是保持 solvent(有偿付能力),即盈余始终为正。破产(Ruin)被定义为盈余首次降至零或零以下。
- 破产时间:\(\tau = \inf\{ t > 0: U(t) < 0 \}\)。如果对所有 \(t \geq 0\) 都有 \(U(t) \geq 0\),则定义 \(\tau = \infty\)。
- 最终破产概率:这是风险理论研究的核心,定义为 \(\psi(u) = P(\tau < \infty | U(0)=u)\),即在无限时间范围内,从初始资本 \(u\) 出发最终发生破产的概率。
为了保证公司有长期持续经营的可能(即 \(\psi(u) < 1\)),保费收入必须足够覆盖平均的索赔支出。这引出了安全负载条件:
\[c > \lambda \mu \]
这个不等式的意义是:单位时间的平均保费收入 \(c\) 必须大于单位时间的平均索赔支出 \(\lambda \mu\)。\(\theta = \frac{c}{\lambda \mu} - 1\) 称为相对安全负载。当 \(\theta > 0\) 时,模型是“有利可图的”,且可以证明此时 \(\psi(u) < 1\)。如果 \(c \leq \lambda \mu\),则 \(\psi(u) = 1\),破产是必然事件。
第四步:破产概率的积分微分方程与Cramér–Lundberg近似
破产概率 \(\psi(u)\) 满足一个重要的积分微分方程。通过考虑在第一个微小时间间隔 \(dt\) 内可能发生的情况(无索赔或有索赔),可以推导出:
\[c \psi'(u) = \lambda \psi(u) - \lambda \int_0^u \psi(u-y) dF(y) - \lambda (1 - F(u)) \]
这是一个更新型的方程。
在满足安全负载条件 \(c > \lambda \mu\) 的前提下,一个里程碑式的结果是 Cramér–Lundberg 不等式:
存在一个正常数 \(R\),称为调节系数或Lundberg指数,使得:
\[\psi(u) \leq e^{-R u}, \quad \forall u \geq 0 \]
这个 \(R\) 是方程
\[\lambda + cR = \lambda M_Y(R) \]
的唯一正根,其中 \(M_Y(r) = \mathbb{E}[e^{r Y_1}]\) 是索赔额 \(Y_1\) 的矩母函数。这个不等式给出了破产概率的一个指数型上界,直观地显示了初始资本 \(u\) 越大,破产概率的上界越小。
更进一步,在一定的正则条件下(例如索赔额分布是“非格点的”,且矩母函数在某个邻域内有定义),破产概率存在一个精确的渐近表达式,即 Cramér–Lundberg 近似:
\[\psi(u) \sim C e^{-R u}, \quad \text{当} \ u \to \infty \]
其中 \(C\) 是一个正常数。这意味着对于大额初始资本,破产概率主要表现出指数衰减行为。
第五步:模型的扩展与应用
基本的 Cramér–Lundberg 模型有许多推广,以贴近更复杂的现实:
- 带投资收益的模型:将盈余的一部分用于投资,比如投资于一个布朗运动或几何布朗运动描述的金融市场,使模型变为带有随机投资的扩散-跳跃过程。
- 带利息力的模型:考虑资金的时间价值,在模型中引入一个常数的利息力 \(\delta\),此时盈余会按指数增长。
- 带再保险的模型:保险公司通过再保险(如比例再保险、超额损失再保险)将部分风险转移给再保险公司,这会影响索赔额的分布和保费收入。
- 扰动模型:在模型中增加一个扩散项(布朗运动)来反映除大额索赔外的、连续的小额波动(如日常运营费用、投资收益的连续波动等),即 \(dU(t) = c dt - dS(t) + \sigma dW(t)\),其中 \(S(t)\) 是累计索赔过程,\(W(t)\) 是标准布朗运动。
Cramér–Lundberg 模型及其理论构成了非寿险精算数学、风险管理以及排队论等领域的重要基石,它提供了一个分析机构(如保险公司、银行)在随机冲击下的生存能力和资本需求的严谨数学框架。