有限尺度效应对流行病房室模型动力学影响的研究

《European Journal of Applied Mathematics》：An exploration of finite-size effects in the dynamics of epidemic compartmental modelling

　　

有限尺度效应（finite-size effects）是指在具有内在随机性的人或粒子交互系统中，其动力学行为偏离相应平均场方程（mean-field equations）预测的现象。尽管在实践中频繁观察到，有限尺度效应塑造和影响真实世界流行病动力学的机制仍大多未被探索。为此，研究人员在基于个体（agent-based）和平均场两种流行病建模框架下研究这些效应，尤其关注以反复暴发为特征的动力学区制。最终目标是为流行病预测与防控的公共卫生策略提供信息并改进之。作为迈向更广义目标的第一步，研究人员聚焦于一个能自然刻画反复流行波的最小流行病模型，即SIHRS房室模型（SIHRS compartmental model），该模型假设康复个体免疫力消退（waning immunity），可视为经典SIHR与SIRS房室模型的自然扩展。SIHRS框架而非用过度复杂结构掩盖底层机制，使研究人员能够隔离、识别并随时间追踪有限尺度效应。分析采用基于鞅（martingale）的早期时间步方法，推导并利用底层马尔可夫过程（Markov processes）的鞅表述。特别地，鞅表述从房室变量扩展到高阶量，这是将分析推进到短时近似之外的重要一步。分析定量且系统地理解了流行病动力学中的有限尺度效应，表明其影响由群体规模、房室尺度和时间结构共同塑造。这些发现由基于国家级与县区级COVID-19数据标定的大量数值模拟支持。

研究背景方面，在人或粒子交互系统中，基于个体的随机模型与确定性平均场常微分方程（ODEs）或偏微分方程（PDEs）之间存在固有差距：平均场方程假定充分混合的大群体，忽视个体异质性与随机性，在有限或异构结构中常失效。近几十年的观测表明，真实流行病并非总沿平均场轨迹演进，尤其在较小局部群体中个体异质性被放大。已有先驱工作仅关注无暴发或单次暴发区制下的有限尺度效应，而反复流行波（recurrent epidemic waves）中随机涨落可跨波次累积，使有限尺度效应逐渐放大且难以解析处理。为隔离并识别该效应，研究人员选取能刻画疾病波动的最小框架——SIHRS房室模型（在S、I、R基础上将R细分为住院H、康复R、死亡D，并引入免疫力消退使R可回流转为S），以此作为透明基线。该研究发表于《European Journal of Applied Mathematics》，旨在填补有限尺度效应在反复暴发区制下定量理解的空白，对公共卫生规划具有重要意义。

关键技术方法上，研究人员构建了带有独立泊松时钟（Independent Poisson Clocks, IPC）的基于个体SIHRS模型（SIHRS-IPC），利用泊松合并与分裂技术得到统计等价表述并导向Gillespie型模拟；采用鞅方法推导连续时间马尔可夫过程的鞅分解，将随机动态拆为无穷小均值积分项与无穷小方差鞅项；由此推出平均场ODE系统，并引入时间积分L²距离度量个体与确定性轨迹偏差，将其概率分解为偏置项与波动项；进一步将鞅表述推广至瞬时感染率等高阶量，结合早期时间渐近展开追踪有限尺度效应的产生与传播。数值实验涵盖国家级COVID-19数据标定及美国内华达州Carson City县区案例（使用该县真实人口与每日病例、住院数据）。

研究结果部分，第2节介绍基于个体的SIHRS-IPC模型与模拟：通过Type(I)–Type(IV)四类泊松时钟分别治理S-I接触、I的转出、H的转出及R的免疫消退，借泊松分合得到等效主权时钟率Λ_N(t)，实现高效随机模拟。第3节给出各房室变量的鞅表述并导出平均场ODE系统。第4节定义时间积分L²差异及其分解为偏置（源于非线性期望）与波动（源于过程方差）两项，偏置与波动均随N减小而增长。第5节分析偏置项：初始微小的漂移项随机偏差随时间累积，导致随机路径期望与确定性轨迹系统性分离。第6节分析波动项：证明随机过程方差上界按1/N缩放，并在非平凡时间区间内保持该逆缩放行为，该区间长度依赖房室与N。第7节Carson City案例显示，即便总体较大，小规模房室（如H）仍呈现显著随机变异性；个体随机路径较平均场更早消亡或呈不规则长期振荡。附录A.4将鞅分析用于瞬时随机感染率，推导其早期时间步表述。

讨论与结论翻译部分，研究人员指出有限尺度效应不仅取决于总人口N，还取决于各房室相对规模及系统时间演化；随机偏差可跨多个流行波累积，导致与平均场持续发散；即便总群体较大，小规模房室（如住院H）仍展现明显随机波动。理论由国家级与Carson City县区COVID-19数据模拟互补验证。研究建立了有限尺度效应在反复暴发区制下形成与演化的综合定量解析理解，明确了确定性平均场描述在此类设置中的内在局限，并表明即便总体较大，随机涨落仍可维持显著影响。所发展的高阶鞅方法不限于流行病，可推广至更广随机交互个体系统。