Hopf bifurcation, bistability and transient grazing in apiecewise-smooth hospital-overload model

技术摘要：分段光滑医院过载模型中的霍普夫分岔、双稳态与瞬态掠过现象

问题陈述
本文针对流行病学建模中的一个关键操作性差距提出了探讨：即感染是否能侵入人群（由基本再生数 $R_0$ 决定）与医疗系统是否能吸收爆发期间产生的重症病例，这两者之间是否存在区别。经典模型通常假设，如果地方病平衡点与常规入院量相兼容，则系统是安全的。然而，在快速爆发的流行波期间，重症病例可能会在入院前在等待状态中累积，从而导致即使长期平衡是可持续的，也会出现暂时的过载。作者旨在分析由容量限制的医院入院所产生的非线性动力学结构，特别关注从需求驱动型向容量驱动型入院机制的转变。

方法论
研究采用了一个连续分段光滑的流行病学模型（SEIQRH），其中入院通量被定义为需求驱动通量（$aq$）与容量驱动通量（$c(b-h)_+$）的最小值。这种表述创建了一个切换曲面$\Sigma$，在该曲面上向量场是连续的，但雅可比矩阵是不连续的。

关键方法论组成部分包括：

• 无量纲化表述： 模型通过缩放将参数减少到无量纲组，从而将医院容量 $b$ 孤立为一个主要的分岔参数。
• 谱约简： 作者识别出两个单侧雅可比矩阵（用于正常护理的 $J_N$ 和用于容量受限护理的 $J_C$）共有的一个精确且与机制无关的特征对。这使得两个机制的六阶特征多项式均可简化为五阶多项式，从而便于进行显式的 Routh–Hurwitz 稳定性分析。
• 局部与全局分岔分析：
  • 局部： 计算第一李雅普诺夫系数以确定霍普夫分岔的临界性。
  • 全局： 使用单射击法（single-shooting formulation）结合弗洛凯乘子（Floquet multipliers）对周期轨道进行数值延拓，以追踪局部霍普夫分支与有限振幅吸引子之间的连接。
  • 瞬态分析： 对“掠过”（grazing）事件进行渐近分析，即正常护理下的疫情轨迹切向接触切换曲面的现象。
• 严格界限： 利用 Metzler 比较系统推导出关于首次切换时间（即首次达到容量限制时刻）的全局有效下界。
• 验证： 通过使用符号计算（SymPy）和高精度算术（mpmath）独立复现结果，以验证霍普夫系数和吸引子的性质。

主要贡献与结果

1. 谱约简与稳定性：
   研究表明，两个单侧雅可比矩阵共享一个对应于人口更替率（$-m$）的精确特征对。这使得稳定性分析可以简化为五阶多项式。对于基准参数，容量受限的平衡点在容量阈值 $b_H \approx 0.2263 b_{crit}$ 处发生亚临界霍普夫分岔。正的第一李雅普诺夫系数（$\ell_1 \approx 42.44$）证实了新出现的局部周期轨道是不稳定的。

2. 双稳态与全局动力学：
   数值延拓揭示了不稳定的局部霍普夫分支通过一个较低的周期折叠点（$b^-_{LPC}$）连接到一个稳定的、有限振幅的分段光滑周期分支。在流行病平衡容量阈值（$b_{crit}$）之外存在第二个折叠点（$b^+_{LPC}$）。这种结构创造了两个截然不同的双稳态区间：

   • $b^-_{LPC} < b < b_H$：稳定平衡点与稳定周期吸引子并存。
   • $b_{crit} < b < b^+_{LPC}$：稳定平衡点与稳定周期吸引子并存。
     这暗示医院在平衡状态下的充足性并不保证动力学上的安全性；系统的状态高度依赖于初始条件和扰动。

3. 瞬态掠过与定标律：
   作者定义了一个瞬态阈值 $b_{trans}$，代表正常护理爆发轨迹产生的最大入院压力。当 $b < b_{trans}$ 时，轨迹会掠过切换曲面。在唯一且非退化极大值的假设下，论文推导出了关于过载持续时间（$\Delta \tau$）和累积入院赤字的特定定标律：

   • 过载持续时间随 $(b_{trans} - b)^{1/2}$ 定标。
   • 累积入院赤字随 $(b_{trans} - b)^{3/2}$ 定标。
     这些结果强调，系统在接近瞬态阈值时对容量短缺的响应是奇异且非线性的。

4. 严格的切换时间界限：
   Metzler 比较论证提供了一个关于首次容量限制时间的严格且全局有效的下界。虽然有用，但论文指出，随着易感人群的耗尽，该界限与真实切换时间之间的相对差距会增大。

5. 经验应用（英格兰国民医疗服务体系 NHS England）：
   论文将一个简化的存量-流量版本的模型应用于 2020/21 年冬季的 NHS England 数据。虽然该模型成功追踪了入院量与床位占用量之间的滞后关系（验证了存量-流量机制），但作者明确指出，仅凭现有数据无法识别潜在的等待类以及具体的容量阈值（$b_{crit}, b_{trans}$）。

意义与主张
该论文声称其主要意义在于将三个动力学上截然不同的问题区分开来：入侵、地方病充足性以及瞬态过载。它认为，平衡状态下的医院充足性并不保证瞬态或动力学上的安全性。

• 非线性复杂性： 该模型揭示了容量限制的入院可以产生复杂的动力学，包括亚临界霍普夫分岔和双稳态，而这些现象是平滑非线性近似或简单阈值模型无法捕捉到的。
• 操作含义： 存在一个瞬态阈值 $b_{trans}$，它可能显著高于地方病阈值 $b_{crit}$（在基准情况下，$b_{trans} \approx 53.5 b_{crit}$），这表明仅基于地方病平衡计算进行的规划是不充分的。增量容量必须考虑到瞬态峰值压力。
• 技术严谨性： 该工作对边界-平衡转换（持续类型）进行了严格分类，并为连续分段光滑系统中的掠过事件建立了精确的渐近定标律，将其与冲击型（Filippov）系统的平方根定标区分开来。

作者对政策建议保持谨慎，指出基准参数仅具说明性，完整的校准需要有关入院、占用、出院、人员配备和等待时间的关联数据，而这些数据在经验检查中并未获得。其结果旨在作为一种结构性演示，展示容量约束如何改变流行病动力学，而非针对特定地区的校准预测。