《Annals of Actuarial Science》：Arbitrage-free catastrophe reinsurance valuation for?compound dynamic contagion claims
摘要：研究人员研究了在时变（time-inhomogeneous）参数设定下的复合动态传染过程（Compound Dynamic Contagion Process, CDCP）——即同时包含自激（self-excited）跳跃与外部激发（externally excited）跳跃的Cox过程且其强度服从Ornstein-Uhlenbeck型均值回复带跳扩散——应用于巨灾（catastrophe）索赔到达建模，并在不完全保险市场中通过Esscher变换确定等价鞅测度（Equivalent Martingale Measure, EMM）以获得再保险合同的无套利（no-arbitrage）保费。研究人员推导了时变CDCP联合过程（λt,Nt,Ct,Mt,Λt,t）在物理测度P及风险中性测度P~下的无穷小生成元（infinitesimal generator）A与A~，给出强度λt、索赔计数Nt与累计索赔Ct一阶矩的显式表达式（含时齐与时间依存情形），证明使保险人盈余过程{Rt}成为P~-鞅的测度变换之存在性，进而得到巨灾再保险仲裁免费（arbitrage-free）保费等于E~[CT∣F0]。该研究将Dassios & Zhao (2011) 的标准DCP推广至时变参数情形，并将Sondermann (1991)–Dassios & Jang (2003) 的再保险无套利框架扩展至具双跳跃 contag效应与随机强度衰减的动态传染模型。
论文解读——《Annals of Actuarial Science》：时变参数复合动态传染过程(CDCP)下巨灾再保险的无套利定价
一、研究背景与意义
传统再保险定价常采用净保费原理或期望值原理，未充分考虑索赔过程的自激传染效应（如一次巨灾引发后续余震或连带损失）、外部随机冲击导致的强度跳变，以及强度本身的均值回复（mean-reversion）与随时间变化的暴露水平。经典Cox过程（Cox with shot-noise intensity）仅含外部激发跳跃，Hawkes过程仅含自激效应，均无法同时刻画巨灾险中两类传染机制。Dassios & Zhao (2011) 提出的复合动态传染过程（Compound Dynamic Contagion Process, CDCP）弥补了这一缺陷，但其原模型假设时齐（time-homogeneous）参数，与实际保险组合随季节、组合构成、风险环境变化而时变的情形不符。此外，再保险市场是不完全市场，不存在唯一等价鞅测度，需借助Esscher变换选取合理的风险调整测度以实现无套利（no-arbitrage / actuarial no-arbitrage）定价。本研究将CDCP推广至时变参数（time-inhomogeneous CDCP, 允许a(t),ρ(t),μG(t),μH(t;t)时间依存），严格推导其无穷小生成元（infinitesimal generator）A及Esscher变换后测度P~下的A~，获得各关键量的矩显式解，并基于使保险人净盈余过程为P~-鞅给出巨灾再保险的无套利保费公式，为不完备市场中 catastrophe reinsurance 的精算公允定价提供理论依据。
二、主要关键技术方法
研究人员采用以下方法开展研究：(1) 建立时变参数CDCP模型——强度λt满足带自激跳（分布G(y;t)，均值μG(t)）与外部激发跳（分布H(x;t)，均值μH(t)）的均值回复SDE，κ(t)=δ?μG(t)；(2) 应用随机过程无穷小生成元（infinitesimal generator）技术，对试函数f(λ,n,c,m,Λ,t)作用A得生成元表达式，通过Dynkin公式导出E[λt∣λ0]满足的一阶线性ODE并求解获时变期望；(3) 构造指数F-鞅ζt=eK(t)θNteB(t)λte?νCte?ΛtψMt/E0[?]，其中B(t),K(t)满足由Af=0导出的非线性ODE，以此作为Radon–Nikodym导数dP~/dP∣FT=ζT实现Esscher变换（Esscher transform）；(4) 利用测度变换下生成元关系A~f~=A(f~e?Bλ)/e?Bλ推导P~下参数变化（自激率乘θj^(ν)，外部激发率乘ψ，强度漂移修正），验证使盈余{Rt}成P~-鞅的条件等价于再保险无套利保费E~[CT∣F0]=P~下累计索赔期望。
三、研究结果（按论文章节小结）
2.1 时齐标准CDCP（Standard Time-Homogeneous CDCP）
研究人员给出时齐CDCP联合过程(λt,Nt,Ct,Mt,Λt,t)在物理测度下无穷小生成元Af（式2.2），其中漂移项δ(a?λ)?f/?λ表均值回复，自激跳项λ[∫∫f(λ+y,n+1,c+ξ,m,Λ,t)dG(y)dJ(ξ)?f(?)]描述claim触发强度跳增，外部激发跳项ρ[∫f(λ+x,n,c,m+1,Λ,t)dH(x)?f(?)]描述外生冲击。取试函数f=λ,f=n,f=c结合Dynkin公式，获命题1：E[λt∣λ0]=κρμH+aδ+(λ0?κρμH+aδ)e?κt（κ=δ?μG>0为平稳条件），E[Nt∣λ0]=μ1t+(λ0?μ1)(1?e?κt)/κ，E[Ct∣λ0]=μJE[Nt∣λ0]，其中μ1=(ρμH+aδ)/κ。
2.2 时变（Time-Inhomogeneous）CDCP
研究人员将a→a(t),ρ→ρ(t),G(y)→G(y;t)(μG(t)),H(x)→H(x;t)(μH(t))，解SDE得λt=λ0e?δt+δ∫0ta(s)e?δ(t?s)ds+∑衰减跳项（式2.5）。给出时变情形下生成元A（式2.9），其中a(t),ρ(t),G(y;t),H(x;t)替代常数。取f=λ得Aλ=?κ(t)λ+ρ(t)μH(t)+a(t)δ，经Dynkin公式导出ODE dμλ(t)/dt=?κ(t)μλ(t)+ρ(t)μH(t)+a(t)δ，解得命题2：μλ(t)=E[λt∣λ0]=λ0e?∫0tκ(s)ds+e?∫0tκ(s)ds∫0te∫0sκ(u)du{ρ(s)μH(s)+a(s)δ}ds；进而E[Nt∣λ0]=∫0tμλ(s)ds，E[Ct∣λ0]=μJ∫0tμλ(s)ds。
3 保险市场与无套利（Insurance Market and No-Arbitrage）
研究人员定义再保险策略?u∈[0,1]（比例分保），保险人净盈余Rt=Pt?Ct（Pt为已收总保费）。定义套利再保险策略为终端增益Gt,T(?)≥0a.s.且E[Gt,T(?)∣Ft]>0。引述Fundamental Theorem：无套利?存在等价鞅测度P~使{Rt}为P~-鞅。若R0=0，则E~[Ct∣F0]=E~[Pt∣F0]=无套利再保险保费。
3.1 等价鞅测度（Equivalent Martingale Measures）
定义P~～P使{Rt}为P~-鞅（Definition 5）。因市场不完全，用Esscher变换选P~。
定理1（Martingale for CDCP）
取f(λt,Nt,Ct,Mt,Λt,t)=eK(t)θNteB(t)λte?νCte?ΛtψMt，令Af=0得B(t),K(t)满足非线性ODE：B′(t)?δB(t)+θj^(ν)g^(?B(t))+??1=0，K′(t)+aδB(t)+ρ[ψh^(?B(t))?1]=0，其中g^,h^,j^为G,J,H的Laplace变换。此时{f(λt,…,t)}t∈T是(F,P)-局部鞅。
定理2（解的存在唯一性）
在θ,ψ≥1,ν<0(j^(ν)>1),?=?(θj^(ν)?1)<0,δ>θj^(ν)μG下，B(t)=G?1(t)唯一存在（G(B)=∫bBdu/f1(u),f1(B)=δB?θj^(ν)(g^(?B)?1)），K(t)=?aδ∫0tB(s)ds+ρ∫0t[1?ψh^(?B(s))]ds唯一。
取ζt=f(λt,…,t)/E0[f(λt,…,t)]为F-鞅>0，dP~/dP∣FT=ζT定义EMM。
3.2 等价鞅测度下无穷小生成元（Infinitesimal Generator under EMM）
引理1给出A~f~(λ,0)=A(f~(λ,0)e?B(t)λt)/e?B(t)λt。
推论1（生成元形式）
在P~下CDCP生成元为：
A~f=?f/?t+δ(a(t)?λ)?f/?λ+λ?f/?Λ
+λθj^(0)[∫∫f(λ+y,n+1,c+ξ,m,Λ,t)dG(y)dJ(ξ)?f(?)]
+ρψ[∫f(λ+x,n,c,m+1,Λ,t)dH(x)?f(?)]
即自激跳率从λ变为λ?θ（因j^(0)=1当取ν=0对应只调整claim计数Nt的Esscher参数），外部激发率从ρ(t)变为ρ(t)ψ，其余漂移项不变。取f=λ得A~λ=?κ(t)λ+θρ(t)μH(t)+a(t)δ，据此可计算E~[λt],E~[Nt],E~[Ct]，无套利再保险保费即为E~[CT∣λ0]。
四、讨论与结论
研究人员将Dassios & Zhao的CDCP从时齐推广至时变参数，严格建立了该过程在物理测度与Esscher变换后风险中性测度下的无穷小生成元，获得强度、索赔次数及累计索赔一阶矩的显式（含积分形式）表达式。通过指数鞅进行Esscher变换，证明在不完全再保险市场中可选取等价鞅测度使保险人盈余为鞅，从而巨灾再保险无套利保费等于风险中性测度下累计索赔的条件期望。测度变换后自激跳跃率被θ倍缩放、外部激发率被ψ倍缩放，体现了Esscher参数对跳跃强度的风险调整。该框架兼顾了巨灾险的自激与外源传染特征、强度均值回复衰减及时变承保暴露，为catastrophe reinsurance的精算公允定价与solvent capital评估提供较通用的随机模型与无套利定价理论基础。论文结论对再保险定价实务中引入动态传染效应与随时间变化的业务参数具指导意义，并为后续研究拓展至多类型再保险（excess-of-loss, stop-loss）及随机利率情形奠定基础。
