链上镖师 / Chain Biaoshi
CIE|灾害救援储备：实物库存与产能预留的风险稳健联合优化
用分布鲁棒机会约束模型平衡响应可靠性与经济成本
论文提出JPPR联合储备策略，将灾前实物库存与供应商生产能力预留纳入同一个两阶段风险厌恶分布鲁棒模型。
- 01：研究对象是灾害救援中的仓库、供应商、实物储备和产能储备联合配置。
- 02：模型同时处理灾后需求、实物库存存活率和供应商产能存活率的不确定性。
- 03：WMCVaR与机会约束共同控制低概率高影响短缺风险。
- 04：墨西哥湾飓风案例包含27节点网络、25个需求地点和15个候选设施位置。
- 05：样本外测试中，DRCC综合响应期可靠性达到99.3%，并保持仓库和供应商风险为零。
灾后物资为什么要把库存与产能放在一起
灾害救援中的储备决策需要同时面对响应速度、库存成本和供应中断风险。实物储备能够快速投入现场，长期采购与保管会带来较高成本；产能储备能够根据需求启动生产，响应时间依赖供应链，生产设施也可能受到灾害影响。论文将两类资源放入同一决策框架，形成联合实物与产能储备策略。
01实物储备：适合灾后初始响应，承担快速供给任务。
02产能储备：适合后续补充供给，承担弹性扩产任务。
03联合配置：同时决定仓库位置、实物数量、供应商选择和预留产能。
图1：联合储备问题的事件序列。来源：Fig. 1
从灾前准备到灾后调配：JPPR如何建模
JPPR采用两阶段决策结构。第一阶段确定仓库和供应商的设施类型、位置、实物储备数量与产能预留数量；灾害发生后，第二阶段根据实际需求、仓库库存存活率和供应商产能存活率安排运输，并计算短缺。目标函数覆盖设施建设、采购、库存持有、产能预留、运输和短缺惩罚成本。
11. 准备阶段：选择仓库与供应商位置及规模。
22. 储备阶段：配置实物库存和生产能力预留。
33. 灾害实现：观测需求与两类储备的存活率。
44. 响应阶段：分期调配物资并计算短缺。
55. 求解阶段：以Wasserstein模糊集和WMCVaR形成MIQCP。
图2：模型构建流程示意。来源：Fig. 2
模型的关键在于保留需求与储备存活率的联合不确定性，并用WMCVaR关注低概率高损失场景。原始分布鲁棒问题经过等价 reformulation 转化为可由CPLEX或Gurobi求解的混合整数二次约束规划。
墨西哥湾案例如何搭建
案例聚焦美国墨西哥湾区域的飓风救援网络，运输网络包含27个节点，其中25个节点承担需求，15个位置作为仓库或供应商候选点。研究对象为饮用水，计量单位为1000加仑。仓库和供应商均设置小型、中型、大型三种规模，历史数据选取15场飓风，并通过有放回抽样生成50个随机飓风情景。
表6：案例描述与参数设置。来源：Table 6
图3：案例研究网络。来源：Fig. 3
- Case A：仅配置实物储备。
- Case B：联合配置两类储备，忽略供应风险。
- Case C1：加入实物储备存活风险。
- Case C2：加入产能储备存活风险。
- Case D：同时考虑两类储备的供应风险。
风险情景如何改变设施与储备布局
五组案例显示，供应风险会改变仓库与供应商的空间组合。Case D选择5个仓库和8个供应商，配置19,450单位实物储备与45,000单位产能储备。考虑实物储备风险后，仓库布局向Shreveport、Atlanta、Orlando、Jackson和Jacksonville等位置移动；考虑产能风险后，供应商组合更加分散，并增加内陆节点参与。
图4：库存与供应商配置。来源：Fig. 4. Configuration of inventory and supplier.
表7：各案例结果。来源：Table 7 Case results.
Case B总成本
395.50百万美元
无供应风险设定
Case D总成本
471.99百万美元
同时考虑两类供应风险
样本外表现：DRCC如何平衡可靠性与成本
论文使用前25个样本确定第一阶段决策，再用剩余25个样本检验样本外表现，并将需求满足率乘积定义为可靠性。SAA具有较低经济成本，却在初始响应期和综合响应期呈现较低可靠性；RO提供较强保护并承担更高经济成本。DRCC处于两者之间，在风险控制、可靠性和经济成本之间取得更稳定的组合。
表8：可靠性、风险与经济成本比较。来源：Table 8
图6：不同训练样本量下SAA的样本外表现。来源：Fig. 6
DRCC可靠性
10.0%→99.3%
初始响应期至综合响应期
SAA可靠性
3.1%→67.7%
初始响应期至综合响应期
在论文给定的飓风样本和参数设置下，DRCC综合响应期可靠性达到99.3%，仓库与供应商风险均为零；其经济成本为5.66×10^7美元，低于RO模型。
风险厌恶与参数半径如何重塑储备策略
敏感性分析表明，风险厌恶权重γ从0增加到1时，仓库数量由4个增加到11个，供应商数量由5个增加到10个，经济成本由51.78百万美元上升到162.43百万美元，短缺惩罚成本由828.45百万美元降至55.28百万美元。初始响应期惩罚较高时，实物储备随γ增加得更快。
表9：不同风险厌恶权重γ下的比较结果。来源：Table 9
表13：储备成本与容量参数对最优策略和网络配置的敏感性。来源：Table 13
✓提高实物储备成本，实物储备量由24,347.5单位降至11,000单位。
✓提高仓库容量，仓库数量由9个降至2个，储备总量保持稳定。
✓提高供应商产能，供应商数量由11个降至6个。
✓提高风险厌恶程度，系统增加设施和储备投入以降低短缺损失。
计算可扩展性与模型适用范围
计算实验显示，模型在30节点和40节点合成网络上分别用707秒和1579秒获得最优解；网络规模扩大到50、60、70和80个需求节点后，计算时间达到7200秒，并出现1.20%、25.81%、26.36%和71.67%的求解间隙。结果说明，MIQCP适合中等规模网络，超大规模部署需要进一步控制模型规模或采用分解求解策略。
表14：较大规模合成实例的计算结果。来源：Table 14
- 适用条件：存在历史灾害事件样本，能够估计需求与两类储备存活率。
- 参数条件：风险厌恶权重、Wasserstein半径、短缺惩罚和设施容量需要结合业务设定校准。
- 数据边界：案例使用墨西哥湾飓风网络、饮用水需求和15场历史飓风，其他灾种与物资的外推范围由数据和参数决定。
- 计算边界：论文报告了大规模实例的时间限制和求解间隙，超大规模网络的精确求解能力有限。
论文的直接管理含义是：初始响应阶段需要保留足量实物库存，后续供给阶段可以利用分散化的产能储备；两类储备的比例应随短缺代价、存活风险和设施容量共同校准。
论文：A distributionally robust optimization approach for pre-positioning in humanitarian relief: Joint Physical and Production-Capacity Reserve strategy
作者：Zizhuo Wang，Xiang Chen，Jianghua Zhang
期刊：Computers & Industrial Engineering
链上镖师、先卜交通、先卜运筹通用论文ID
保存此推文 ID，可在粉丝群内获取相关材料，未来更多功能敬请期待
先卜科技团队：常文婕、张巧巧、王一航、夏子帆、张晓乐、祁园园、张仁义
