Situation(背景)
某码头设有 c 个泊位供货轮停靠装卸。由于泊位有限,到港货轮需排队等待。保险公司计划推出一款保险产品:货轮到港时缴纳保费 P,若实际等待时间超过承诺时限 T,保险公司支付赔偿金 C(C ≫ P)。保费须与承诺时间关联,需为保险公司构建定价模型,探讨盈亏平衡与利润最大化条件。
Task(任务)
- 将码头作业抽象为数学模型,导出等待时间概率分布
- 建立保险公司期望利润函数,分析盈亏平衡条件
- 针对保费与承诺时间的不同耦合关系,构造多种定价模型并求解利润最大化
- 通过数值实验验证理论推导,分析参数敏感性
Action(行动)
排队论模型:M/M/c 系统
将码头作业抽象为 M/M/c 排队系统——货轮到达服从 Poisson 过程(到达率 λ),每个泊位服务时间服从负指数分布(服务率 μ),FCFS 排队。在稳定性条件 ρ = λ/(cμ) < 1 下,导出稳态概率分布与 Erlang C 公式(排队等待概率):
$$P_Q = \frac{\dfrac{a^c}{c!} \cdot \dfrac{1}{1-\rho}}{\displaystyle\sum_{n=0}^{c-1} \dfrac{a^n}{n!} + \dfrac{a^c}{c!} \cdot \dfrac{1}{1-\rho}}$$
等待时间超出承诺时限 T 的概率为:
P(W > T) = PQ ⋅ e−(cμ − λ)T
保险利润模型
单船期望利润:
π(P, T) = P − C ⋅ PQ ⋅ e−(cμ − λ)T
盈亏平衡条件:P ≥ C ⋅ PQ ⋅ e−(cμ − λ)T。
三种定价模型
模型 A:精算公平定价 + 安全附加
P(T) = (1 + θ) ⋅ C ⋅ PQ ⋅ e−(cμ − λ)T
保险公司始终盈利,无约束最优为 T → 0。该模型只考虑保险人视角,未纳入需求侧约束。
模型 B:指数衰减定价(市场驱动型)
P(T) = p0 e−γT
利润最大化的一阶条件给出最优承诺时间:
$$T^* = \frac{1}{(c\mu - \lambda) - \gamma} \ln\left(\frac{(c\mu - \lambda) C P_Q}{\gamma p_0}\right)$$
极大值条件:保费衰减率 γ 必须低于系统过剩服务率 cμ − λ。
模型 C:线性定价
P(T) = p0 − kT
最优承诺时间:
$$T^* = \frac{1}{c\mu - \lambda} \ln\left(\frac{(c\mu - \lambda) C P_Q}{k}\right)$$
二阶条件始终满足。
市场约束
引入两项市场合理性约束: - 最大附加系数 Lmax = 3:保费不得超过期望赔付的 3 倍 - 最低赔付概率 pmin = 0.5%:保单实际赔付概率不得低于 0.5%
数值实验
以基准参数 (λ = 8, μ = 3, c = 3, C = 100000, p0 = 8000, γ = 0.5) 开展系统的数值实验,使用 Python(numpy/scipy/matplotlib)实现排队系统仿真、三种定价模型求解、参数扫描与可视化。
Result(结果)
核心发现
1. 约束收窄策略空间
无约束 FOC 给出的内点解(指数模型 TFOC* = 5.99 天)不满足市场约束。引入 Lmax 与 pmin 后,可行区间收窄至 [4.61, 5.06] 天(宽度仅 0.45 天),最优解移至 pmin 约束边界。
2. 保费衰减率 γ 的可行区间远比理论预期狭窄
理论区间为 γ ∈ [0.33, 1.0],但实际可行区间仅为 γ ∈ [0.25, 0.519],宽度为理论区间的 27%。γ > 0.519 时盈亏平衡下界超越 pmin 约束上界,可行区间闭合。
3. 利润随泊位数 c 单调递增
c 从 3 增到 6 时,单船利润从 131 元跃升至 4780 元(增幅约 36 倍)。机制:c 增大同时压低 PQ(赔付基准缩小)和抬高 cμ − λ(风险衰减加快),而保费仅以 γ 速率衰减——两者衰减速率之差持续扩大。
4. 到达率 λ 对利润的系统性压制
λ 增大只产生单方向不利影响(PQ 上升),缺乏对冲机制。λ ≥ 8.5 时指数模型丧失可行解,该阈值先于系统稳定性极限 λ = 9 出现——市场约束比系统稳定性设定了更严苛的准入门槛。
5. 利润的数学本质
$$\pi(T) = \underbrace{p_0 e^{-\gamma T}}_{\text{保费端}} - \underbrace{C P_Q e^{-(c\mu-\lambda)T}}_{\text{赔付端}}$$
利润来源于两项指数衰减的速率差。保险公司的盈利空间不取决于保费的绝对水平,而取决于保费与赔付随 T 下降的”速度差”。
三种定价模型对比
| 模型 | 最优 T* | 单船利润 | 特点 |
|---|---|---|---|
| 精算公平 | 0 天(无约束)/ 5.07 天(约束) | 23925 元 / 约束下更低 | 始终盈利,未考虑需求侧 |
| 指数衰减 | 5.06 天(约束边界) | 131 元 | 可行区间极窄,利润微薄 |
| 线性 | 3.65 天 | 4097 元 | T 较小时保费高,利润显著占优 |
可视化
数值实验生成 7 张分析图:三种定价模型利润曲线对比、γ 敏感性、λ 系统性影响、(T, λ) 利润热力图、泊位数影响、γ 连续扫描、保费构成分解。
项目文件
完整数学建模报告(含公式推导、数值实验、7 张可视化图表、Python 代码)位于项目目录。
代码由三个 Python
文件协同完成:insurance_model.py(模型库)、plot_analysis.py(可视化)、compute_tables.py(数值表格)。