《Frontiers in Sustainable Food Systems》：Cascading failure mechanisms in the fresh agricultural product supply chain of China’s Yangtze River Delta
编辑推荐：
由于产品易腐性和企业间紧密耦合，新鲜农产品供应链极易发生级联失效。基于复杂适应系统（complex adaptive systems, CAS）和供应链韧性理论，研究人员利用企业注册记录和上市公司披露信息，构建了中国长三角的真实供应链网络，包含1,669个企业
由于产品易腐性和企业间紧密耦合，新鲜农产品供应链极易发生级联失效。基于复杂适应系统（complex adaptive systems, CAS）和供应链韧性理论，研究人员利用企业注册记录和上市公司披露信息，构建了中国长三角的真实供应链网络，包含1,669个企业节点和6,139条供需边。易腐性被直接嵌入级联失效模型：Sigmoid函数连续刻画过载和欠载失效概率，负载重新分配权重随物理距离衰减，累积负载阈值捕捉持续温和冲击下的渐进劣化。DebtRank算法结合负载权重重新识别核心风险节点；动态恢复指标刻画级联后自修复；并通过成本-收益权衡求解库存弹性系数（ε）和距离衰减系数（λ）的最优值。随机攻击和蓄意攻击下的模拟表明，网络高度异质（γ = 1.214），对随机攻击具有鲁棒性但对蓄意攻击脆弱：仅移除10%的最大度节点就会破坏约80%的连通效率，呈现非线性临界相变。累积阈值机制在持续温和冲击下带来真实但较小的额外损失，且在两种攻击类型间基本对称，而非韧性损失的主导因素；DebtRank揭示了被度中心性低估的隐藏风险节点；级联后恢复以相似速率恢复失效节点，而不论攻击类型，但连通效率主要在随机攻击后恢复，因为不区分重要性的恢复过程不太可能恢复连通性所依赖的特定枢纽；ε呈现0.35–0.45的平台型最优区间，而地理集中度（λ）的韧性收益快速饱和。四项定量策略（分级枢纽保护、适度库存缓冲、多元化应急协作和预置恢复）直接支持SDG 12和SDG 13。
近年来，重大突发事件频繁冲击全球供应链，造成巨额经济损失，而新鲜农产品因易腐性和短保质期，对冷链物流时效高度依赖，其流通损耗率远高于发达国家水平，供应链中断还会额外消耗应急冷链运输能源，分别对应联合国可持续发展目标SDG 12和SDG 13。现有级联失效研究存在四个主要缺口：一是大多依赖合成网络或小样本数据，缺乏真实大规模网络的实证基础；二是双向失效模型多采用确定性硬阈值跳变，无法刻画企业在持续温和冲击下的渐进劣化过程；三是核心风险节点识别过度依赖度中心性，忽视负载规模和冲击传播能力的联合信息；四是普遍忽略级联结束后的动态恢复能力。为此，研究人员基于复杂适应系统（CAS）和供应链韧性理论，以中国长三角真实的新鲜农产品流通企业为对象，构建了包含1,669个节点和6,139条供需边的有向加权供应链网络，提出了融合双向Sigmoid失效概率、累积负载阈值、物理距离衰减负载重分配和运输时间新鲜度损失的级联失效模型，并结合DebtRank算法和动态恢复指标，系统研究了随机与蓄意攻击下网络的韧性演化规律及关键参数最优配置。
研究得出以下主要结论：第一，网络具有极端度异质性（γ = 1.214），对随机攻击鲁棒，但对蓄意攻击极为脆弱，仅移除10%的最高度节点就导致连通效率损失约80%，呈现非线性临界相变特征。第二，累积负载阈值机制在持续温和冲击下带来真实但较小的额外损失，且在不同攻击类型间基本对称，并非韧性损失的主导因素，这源于网络负载分布极度右偏，大部分瞬时违规已足以通过基础过载/欠载通道引发失效。第三，DebtRank算法识别的核心风险节点与传统度中心性结果高度一致（Spearman ρ ≈ 0.71），但额外揭示了若干度值不高却拥有大库存或跨区域桥接位置的隐藏风险节点，单靠度指标会低估其系统性威胁。第四，网络动态恢复能力按韧性维度分化：节点数量恢复率几乎与攻击类型无关，而连通效率恢复在随机攻击后远好于蓄意攻击，原因在于盲目恢复过程难以优先恢复连通性所依赖的高中心性枢纽。第五，库存弹性系数ε在随机攻击下效果显著，且存在0.35–0.45的平台型最优区间；地理集中度系数λ的韧性收益快速饱和，在计入构建成本后最优配置趋于低集中度。第六，在多波混合攻击场景下，节点存活率Z(t)随定向攻击比例ρ增大而近似线性下降，但连通效率E(t)在ρ=0.125时即已实现总下降量的约82%，表明少量定向攻击即可促使网络结构性连通崩溃。这些结论表明，韧性提升应侧重于网络结构层面的前摄性保护，而非单纯依赖运营参数调整。
关键技术方法方面，研究人员从国家企业信用信息公示系统和上市公司年报/招股说明书获取企业注册信息及上下游供需关系，经筛查和去重得到1,669个节点和6,145条有向边，并利用CSMAR数据库交叉验证关键边。随后将网络分为供应商、产地批发商、销地批发商和零售商四个层级，基于AMap API获取真实物流距离。模型核心包括：基于Sigmoid函数的连续双向失效概率（过载/欠载）、带记忆衰减的累积负载阈值、按物理距离衰减的度-距离复合负载重分配权重（式11），以及运输时间折算的新鲜度损失系数；采用DebtRank算法对节点负载权重和损伤传播路径进行迭代求解识别核心风险节点；引入恢复动态指数（RDI）和平均恢复时间（MTTR）刻画恢复过程；通过成本-收益联合优化求解决策参数。
研究结果部分的主要发现包括：“网络拓扑性质”显示网络由229个连通分量组成，最大连通分量含814个节点，平均路径长度5.02，网络效率0.248，具有小世界性和碎片化特征；“随机与蓄意攻击的韧性差异”表明蓄意攻击高选择性破坏枢纽节点，导致连通效率急剧下降，而随机攻击下网络衰减较缓，且随攻击比例增大两者趋于一致；“关键参数敏感性”显示ε在随机攻击下显著提升韧性但在蓄意攻击下几乎无效，λ在两种场景下均提供稳定且较小的正向收益，β的影响在β≈0–1时敏感性较低，其后趋于饱和；“累积负载阈值的作用”表明该机制额外损失约0.1个百分点，且不随θ或η显著变化，说明结论对专家判断参数稳健；“DebtRank核心风险节点”揭示三个度排名14–22的隐藏枢纽，若按DebtRank目标攻击同等规模节点，其连通损失略小于纯按度排序，说明两者优化目标不同；“动态恢复能力”表明RDI(Z)与攻击类型无关，RDI(E)则显著依赖攻击类型，MTTR恒约9.4–9.5周期；“ε和λ最优阈值”通过网格搜索求得U(ε,λ)在ε≈0.43处取得峰值，ε最优区间0.35–0.45，λ最优位于低浓度端；“多波混合冲击”证实了枢纽脆弱性结论在动态攻击下依然成立。
讨论部分强调，韧性建设应前移至结构保护，包括基于DebtRank和度中心性的双指标分级枢纽保护、适度库存缓冲（避免过度冗余成本陷阱）、多元化地理分散的应急协作网络（无需刻意提高地理集中度）以及预置恢复能力。这些策略直接回应了SDG 12（负责任的消费和生产）和SDG 13（气候行动）。研究结论指出，长三角新鲜农产品供应链的脆弱性并非节点失效的线性叠加，而是微观主体适应性行为在极端异质网络中耦合涌现的临界相变现象；运营级参数调整难以应对结构性冲击，韧性提升应从网络结构层面实施差异化干预，并辅以预置的灾后恢复能力。同时，研究也承认当前网络基于公开数据构建，无法反映交易频率和品类差异，参数η、Pr、μ等基于文献而非实测校准，未来有待结合企业运营数据进一步验证。
