Some continuity estimates for ruin probability and other ruin-related quantities

基本信息

• 论文ID: 2511.12218
• 标题: Some continuity estimates for ruin probability and other ruin-related quantities
• 作者: Lazaros Kanellopoulos (University of the Aegean, Greece)
• 分类: math.PR (Probability Theory)
• 发表时间: November 18, 2025 (arXiv预印本)
• 论文链接: https://arxiv.org/abs/2511.12218

摘要

本文研究经典风险模型中破产概率的连续性性质。利用压缩积分算子的性质推导破产时赤字的连续性估计。这些结果也应用于扩散扰动的连续时间盈余过程，在该框架下，破产概率可表示为复合几何分布K与扩散项的卷积。文章推导了K的连续性不等式，并提出了该破产相关量的迭代近似方法，通过数值例子验证了结果。

研究背景与动机

核心问题

1. 稳定性问题：在保险数学中，实际的风险参数（如索赔到达强度λ和索赔额分布F）通常是未知的，需要用近似参数来估计。如何度量参数扰动对破产概率等输出量的影响？
2. 连续性估计：寻找形如 ν(ψ_α, ψ_α̃) ≤ w(δ(α, α̃)) 的不等式，其中δ度量输入参数的距离，ν度量输出函数的距离。
3. 扩散扰动模型：对于包含扩散项的风险模型，现有文献缺乏使用概率度量进行连续性分析的工作。

研究重要性

• 理论意义：完善风险理论中的稳定性分析框架，为参数估计误差的传播提供理论保证
• 实践价值：帮助精算师评估使用近似参数计算破产概率时的误差界限
• 方法创新：将Banach压缩原理和不动点理论应用于破产理论的新领域

现有方法局限性

• 强稳定性方法主要应用于无扩散的经典模型
• 显式解仅在特殊分布（指数、相位型）下存在
• 扩散扰动模型的连续性分析工具不足

核心贡献

1. 破产概率的连续性不等式（定理1）：对于具有有限(γ+1)阶矩的索赔分布，建立了基于度量ν_γ的破产概率连续性界
2. 破产时赤字的连续性界（定理2）：推导了两个破产时赤字分布之间的一致度量界，形式为： $\sup_{u≥0}|G(u,y) - \tilde{G}(u,y)| ≤ \frac{1}{c-λμ}[λQ_y(F,\tilde{F}) + |λ-\tilde{λ}|\tilde{μ}]$
3. 扩散模型中K(u)的连续性（定理3）：首次对扩散扰动模型中的复合几何分布K(u)建立连续性估计
4. 迭代近似方法（引理1）：基于Banach不动点定理提出K(u)的显式迭代逼近公式，具有明确的收敛性保证
5. 数值验证：通过多个例子验证理论界的有效性（虽然界相对宽松）

方法详解

任务定义

给定两个风险过程：

• 过程1：参数(λ, c, F)，破产概率ψ(u)，破产时赤字G(u,y)
• 过程2：参数(λ̃, c, F̃)，破产概率ψ̃(u)，破产时赤字G̃(u,y)

目标：建立形如 d_out(ψ, ψ̃) ≤ f(d_in(λ,F; λ̃,F̃)) 的不等式

其中d_in和d_out是适当选择的度量，f是递增函数满足f(0)=0。

理论框架

1. Banach空间构造

• D_γ空间：函数h满足h_γ(x) = (1+x)^γ|h(x)| ∈ D_∞（右连续左极限函数）
• 度量：ν_γ(x,y) = ∫_0^∞ (1+t)^γ|x(t)-y(t)|dt
• 完备性：(D_γ, ν_γ)是完备度量空间

2. 压缩算子

定义积分算子： $Tx(u) = \frac{λ}{c}\left[\int_u^∞ F(t)dt + \int_0^u x(u-t)F(t)dt\right]$

关键性质：

• T: D_γ → D_γ
• 压缩模数：λM_γ^X/c < 1（在净利润条件下）
• 其中 M_γ^X = ∫_0^∞(1+t)^γF(t)dt = E(X+1)^(γ+1)-1/(γ+1)

3. 不动点方法

• 破产概率ψ是算子T的唯一不动点：ψ = Tψ
• 利用三角不等式：ν_γ(ψ,ψ̃) ≤ ν_γ(Tψ,Tψ̃) + ν_γ(Tψ̃,T̃ψ̃)
• 第一项由压缩性控制，第二项度量算子差异

定理1详解（破产概率连续性）

陈述：假设EX^(γ+1) < ∞，EX̃^(γ+1) < ∞，则： $\int_0^∞(1+t)^γ|ψ(t)-\tilde{ψ}(t)|dt ≤ \frac{c}{c-λM_γ^X}\left[\frac{ν_{γ+1}(F,\tilde{F})}{γ+1} + ν_γ(F,\tilde{F})M_γ^L + \frac{|λ-\tilde{λ}|}{c}M_{γ+1}^{\tilde{X}}(1+M_γ^L)\right]$

证明要点：

1. 建立ν_γ(Tx,Ty) ≤ (λM_γ^X/c)ν_γ(x,y)（压缩性）
2. 估计ν_γ(Tψ̃,T̃ψ̃)，分解为I_1（分布差异）+ I_2（强度差异）
3. 对I_1使用换序积分：
   • ∫0^t(1+u)^γdu项给出ν{γ+1}(F,F̃)/(γ+1)
   • 卷积项利用不等式(1+z+t)^γ ≤ (1+z)^γ(1+t)^γ
4. 对I_2类似处理

特殊情况（γ=0）：退化为Kantorovich度量下的界

定理2详解（破产时赤字）

背景：破产时赤字尾部G(u,y) = ψ(u) - G(u,y)满足： $G(u,y) = \frac{λ}{c}\left[\int_0^u G(u-t,y)F(t)dt + \int_{u+y}^∞ F(t)dt\right]$

陈述： $\sup_{u≥0}|G(u,y)-\tilde{G}(u,y)| ≤ \frac{1}{c-λμ}[λQ_y(F,\tilde{F}) + |λ-\tilde{λ}|\tilde{μ}]$

其中Q_y(F,F̃) = ∫_y^∞|F(t)-F̃(t)|dt

证明策略：

• 使用一致度量ν_d(x,y) = sup_t|x(t)-y(t)|
• 固定y，将G(·,y)视为u的函数
• 验证算子T^z的压缩性（模数λμ/c）
• 估计ν_d(T^zΓ̃^z, T̃^zΓ̃^z)

定理3详解（扩散模型中的K(u)）

模型：U(t) = u + ct - S(t) + σB(t)

关键量：K(u) = Pr(L_K > u)，其中L_K是由索赔引起的记录高点的聚合损失

结构：K(u)满足瑕积分方程： $K(u) = \frac{λμ}{c}\left[A(u) + \int_0^u K(u-t)dA(t)\right]$

其中A(x)是L_o + L_c的分布函数，L_o ~ H_1(扩散项)，L_c ~ F_e(平衡分布)

连续性界： $\sup_u|K(u)-\tilde{K}(u)| ≤ \frac{1}{c-λμ}\left[λμ\left(\frac{c}{D}K(H_1,\tilde{H}_1) + \frac{|\tilde{D}-D|}{D} + \frac{K(F,\tilde{F})}{μ} + \frac{|\tilde{μ}-μ|}{μ}\right) + |λμ-\tilde{λ}\tilde{μ}|\right]$

技术难点：需要估计∫|a(t)-ã(t)|dt，其中a是H_1和H_2卷积的密度

迭代近似（引理1）

迭代格式：K_n = T_dK_，K_0 ∈ K_d

显式公式（K_0=k为常数）：

• n=1: K_1(u) = φ - (1-k)φA(u)
• n≥2: K_n(u) = φ - (1-k)φ^nA^{*n}(u) - (1-φ)∑_^{n-1}φ^iA^{*i}(u)

其中φ = 1/(1+θ)

收敛性：由BFPT保证ν_d(K_n,K) → 0

特殊情况（指数索赔，β=c/D）： $K_n(u) = φe^{-βu}S_1(βu) + \sum_{m=2}^{n-1}φ^me^{-βu}[S_{2m-1}(βu)-S_{2m-3}(βu)] + ...$

其中S_m(z) = ∑_^m z^r/r!

实验设置

数值例子设计

例1（定理1验证）：

• 分布：X ~ 混合指数(1/2·Exp(5/4) + 1/2·Exp(5/6))，X̃ ~ Exp(1)
• 参数：EX = EX̃ = 1，多组(λ,c)组合
• 度量：γ = 0, 1

例2（定理2验证）：

• 三个分布：X^(1) ~ Erlang(3,3)，X^(2) ~ Exp(1)，X^(3) ~ 混合指数
• 参数：θ = 1或4，多组(u,y)
• 比较：|G^(i)(u,y) - G^(j)(u,y)|与理论界DK2

例3（定理3验证）：

• 分布：X ~ Exp(3)，X̃ ~ 混合指数
• 扩散参数：多组(D,D̃)组合
• 比较：sup|K(u)-K̃(u)|与界DK3

例4（迭代近似）：

• 分布：X ~ Exp(β)
• 初值：k = 0.0, 0.1, ..., 1.0
• 迭代次数：n = 1,...,5
• 两组参数设置

计算方法

• 精确值：利用已知显式解（指数、Erlang情况）
• 数值积分：计算度量和界
• 迭代计算：使用引理1的显式公式

实验结果

例1结果（表1）

观察：

• γ=0, λ=5/6时：
  • 真实距离：0.0154 (c=3) → 0.0060 (c=7)
  • 理论界DK1：0.1211 → 0.0929
  • 界的宽松度：约8-15倍
• γ=1时界更宽松（约8-13倍）
• λ增大时界略微增大
• c增大时真实距离和界都减小

结论：界虽然保守但定性趋势正确

例2结果（表2）

四组比较：

1. θ=1, X^(1) vs X^(2)：真实距离0.0080-0.0832，界0.1547-0.2900
2. θ=4, X^(1) vs X^(2)：真实距离0.0034-0.0309，界0.0387-0.0735
3. θ=4, X^(1) vs X^(3)：真实距离0.0035-0.0328，界0.0413-0.0779
4. θ=4, X^(3) vs X^(2)：真实距离0.0001-0.0023，界0.0027-0.0054

发现：

• θ增大（安全负载增加）时，距离和界都显著减小
• y增大时，距离先增后减（存在最大值）
• 相似分布（X^(3)与X^(2)）的距离很小（< 0.003）
• 界的宽松度：约2-10倍，在θ=4时更紧

作者备注：界的宽松是因为需要对所有u一致成立

例3结果（表3）

扩散参数影响：

• D=1, D̃=1/10：真实距离0.0854，界0.4837（5.66倍）
• D=3, D̃=1/20：真实距离0.1334，界0.5254（3.94倍）

趋势：

• D与D̃差异越大，距离越大
• D增大时，真实距离增大
• 界的宽松度：约3-7倍

比率分析：|1-Ratio|列显示相对误差为2.94-6.75

例4结果（表4-5）

收敛速度（表4，β=2，精确值K(1)=0.3325717）：

• n=1：0.2030-0.5000（取决于k）
• n=3：0.3316-0.3336（误差<0.001）
• n=5：0.3326左右（误差<10^-5）

初值敏感性：

• k=0.0时n=5误差：0.0000008
• k=1.0时n=5误差：0.0000007
• 结论：5次迭代后对初值不敏感

表5（β=3/2，精确值K(1)=0.6573777）：

• 类似的快速收敛
• n=5时所有初值误差<3×10^-6

实用价值：3-5次迭代即可获得高精度近似

相关工作

稳定性方法

• Beirlant & Rachev (1987)：首次形式化保险数学中的稳定性问题
• Kalashnikov (2000)：强稳定性方法用于破产概率
• Enikeeva et al. (2001)：破产概率的连续性估计
• Gordienko & Vázquez-Ortega (2016)：使用压缩映射的简单连续性不等式（本文直接扩展）

多维和马尔可夫调制模型

• Benouaret & Aı̈ssani (2010)：二维经典风险模型的强稳定性
• Rusaityte (2001)：马尔可夫调制+投资模型的稳定性界
• Gordienko & Vázquez-Ortega (2018)：多维更新风险模型的连续性不等式

压缩原理应用

• Gajek & Rudź (2018)：Banach压缩原理用于regime-switch模型
• Jiang (2021)：q-尺度函数和最终破产概率
• Sánchez & Baltazar-Larios (2020)：使用BFPT的破产概率近似
• Gajek & Rudź (2025)：破产问题的不动点分析

扩散模型

• Gerber (1970)：引入扩散扰动的经典模型
• Dufresne & Gerber (1991)：系统研究三类破产概率
• Tsai (2003, 2006)：K(u)的结构和性质
• Veraverbeke (1993)：渐近行为分析

本文创新：首次将概率度量方法应用于扩散扰动模型

结论与讨论

主要结论

1. 理论贡献：
   • 建立了三个主要的连续性不等式（定理1-3）
   • 提供了可计算的误差界
   • 扩展了压缩算子方法到扩散模型
2. 方法论贡献：
   • 统一框架：选择合适的Banach空间和度量
   • 迭代近似：显式公式+收敛保证
   • 数值验证：多个例子支持理论结果
3. 实践意义：
   • 参数估计误差的定量评估
   • 近似方法的误差控制
   • 扩散模型中K(u)的快速计算

局限性

1. 界的宽松性：
   • 数值例子显示界比真实距离大3-15倍
   • 原因：需要对所有u一致成立，且使用了多次三角不等式
   • 改进空间：更精细的分析或自适应界
2. 矩条件：
   • 定理1要求EX^(γ+1) < ∞
   • 排除了某些重尾分布（如Pareto with α≤γ+1）
   • 对于重尾情况需要其他方法
3. 数值例子的局限：
   • 主要使用指数族分布（有显式解）
   • 更一般分布的验证有限
   • 高维情况未涉及
4. 扩散模型的假设：
   • 定理3假设D≥D̃和μ≥μ̃
   • 这些假设的必要性未充分讨论
   • 不满足时的处理方法未给出
5. 计算复杂度：
   • 界的计算涉及多个积分
   • 实际应用中的可行性未详细讨论

未来方向

1. 理论改进：
   • 寻找更紧的界
   • 重尾分布的处理
   • 多维模型的扩展
2. 方法扩展：
   • 其他类型的风险模型（如Lévy过程）
   • 时间依赖参数
   • 相依索赔
3. 应用研究：
   • 实际数据的参数估计
   • 稳健性分析
   • 最优近似策略

深度评价

优点

1. 数学严谨性：
   • 完整的证明细节
   • 清晰的数学框架
   • 条件假设明确
2. 方法创新：
   • 首次将概率度量系统应用于扩散模型
   • 统一处理破产概率和破产时赤字
   • 迭代近似的显式公式（引理1）
3. 结构清晰：
   • 从简单到复杂（经典→赤字→扩散）
   • 每个定理配有数值验证
   • 技术细节和直观解释平衡良好
4. 实用价值：
   • 提供可计算的误差界
   • 迭代方法收敛快（3-5次）
   • 适用于参数敏感性分析
5. 文献综述：
   • 全面覆盖相关工作
   • 清晰定位本文贡献
   • 引用恰当（32篇参考文献）

不足

1. 界的有效性：
   • 数值例子显示界较宽松（3-15倍）
   • 可能限制实际应用价值
   • 缺少界的紧性分析（是否可改进？）
2. 数值实验：
   • 主要限于指数族（易于计算精确解）
   • 缺少一般分布的蒙特卡罗验证
   • 未讨论计算成本
3. 理论局限：
   • 矩条件排除重尾情况
   • 定理3的单侧假设（D≥D̃）缺少讨论
   • 最优度量选择的指导不足
4. 应用指导：
   • 缺少如何选择γ的建议
   • 实际参数估计的集成未讨论
   • 与其他方法（如强稳定性）的比较不足
5. 写作细节：
   • 一些符号未充分解释（如M^L_γ的含义）
   • 定理2的备注提到"一致性"但未展开
   • 例4的参数选择动机不明

影响力评估

理论影响：

• 为破产理论提供了新的分析工具
• 扩展了压缩算子方法的应用范围
• 可能激发更紧界的研究

实践影响：

• 中等：界的宽松性限制直接应用
• 迭代方法有实用价值（快速收敛）
• 需要进一步工作才能用于实际风险管理

可复现性：

• 高：数值例子参数完整
• 理论证明详细
• 缺少代码但公式足够实现

领域贡献：

• 填补了扩散模型连续性分析的空白
• 提供了系统的方法论框架
• 为后续研究奠定基础

适用场景

适合的场景：

1. 参数敏感性分析：评估参数估计误差的影响
2. 近似方法验证：为数值方法提供误差界
3. 理论研究：作为更复杂模型的基准
4. 教学：展示压缩算子方法的应用

不适合的场景：

1. 精确误差控制：界可能过于保守
2. 重尾分布：矩条件不满足
3. 高维复杂模型：方法扩展非平凡
4. 实时计算：积分计算可能较慢

改进建议：

• 结合蒙特卡罗方法验证更一般情况
• 开发自适应界（根据具体参数调整）
• 提供软件包实现
• 与机器学习方法结合用于参数估计

参考文献（精选）

基础文献：

• Asmussen & Albrecher (2010): 破产概率理论的标准教材
• Dufresne & Gerber (1991): 扩散扰动模型的奠基性工作

方法论：

• Gordienko & Vázquez-Ortega (2016): 本文直接扩展的工作
• Gajek & Rudź (2018, 2025): Banach压缩原理的应用

扩散模型：

• Tsai (2003, 2006): K(u)的结构和性质
• Veraverbeke (1993): 渐近行为

总体评价：这是一篇数学严谨、方法创新的理论论文，为破产理论提供了新的分析工具。主要贡献在于将概率度量方法系统应用于扩散扰动模型，并提供了可计算的连续性界。虽然数值界相对宽松限制了直接应用价值，但作为理论基础和方法论贡献是显著的。论文适合发表在精算或应用概率的专业期刊上，对该领域的研究者有重要参考价值。未来工作应聚焦于改进界的紧性和扩展到更一般的模型。