SIR模型疫情计算器

使用确定性传染病模型,模拟易感、感染和康复人群的变化。

SIR疫情模拟
输入各类人群的初始人数、每日传播率和康复率,以及模拟时长。

了解SIR疫情模型

SIR模型是研究传染病传播的基础数学框架。它将封闭人群划分为三类:可能感染的易感者、能够传播疾病的感染者,以及假定已具备免疫力且不再具有传染性的康复者。模型名称源于这三类人群的英文首字母S、I和R。每个时间步中,人群由易感状态转为感染状态,再由感染状态转为康复状态;总人口保持不变。 传播率用beta表示。每个模拟日的新增感染人数等于beta乘以易感人数和感染人数,再除以总人口。康复率用gamma表示,新增康复人数等于gamma乘以感染人数。相应的微分方程为dS/dt = -beta SI/N、dI/dt = beta SI/N - gamma I和dR/dt = gamma I。本计算器使用步长为一天的欧拉法近似这些连续方程,并限制每日转移人数不超过来源人群的人数。 在最简单的SIR假设下,beta与gamma之比就是基本再生数。大于一表示当几乎所有人都易感时,感染人数最初可能增长;小于一则表示下降。随着易感人数减少,有效再生数会降低,因此即使基本再生数大于一,疫情也可能达到峰值。结果中的峰值是模拟期间最大的感染人数及其出现日期。最终各类人数描述的是所选时间终点的模型状态,并不一定代表疫情彻底结束时的情况。 现实疫情很少满足所有SIR假设。人群存在年龄差异、接触网络、出生、死亡、迁移、疫苗接种、行为变化、不同的传染期、漏报及再次感染。beta和gamma也可能随时间改变。一天的数值步长适合学习和快速情景分析,但精度不及更小步长的数值求解器或拟合后的统计模型。请勿将输出用作预测或公共卫生决策依据。可靠的预测需要经过校准的数据、不确定性区间、敏感性分析,以及熟悉病原体和监测系统的流行病学专家审核。

SIR模型示例

以下示例说明传播率和康复率假设如何影响疫情曲线。

初始情景预期变化解读
N 10,000;S 9,990;I 10;beta 0.3;gamma 0.1;100天R0为3,初期感染增长社区疫情示例
N 1,000;S 990;I 10;beta 0.08;gamma 0.1;60天感染人数总体下降R0低于1
N 50,000;S 39,990;I 10;R 10,000;beta 0.2;gamma 0.1增长慢于全体易感的人群既有免疫力降低实际传播

如何运行SIR模拟

  1. 输入总人口以及初始易感、感染和康复人数,确保三者之和等于总人口。
  2. 以小数输入假设的每日传播率和康复率。
  3. 选择要模拟的整数天数。
  4. 运行模型,比较最终各类人数、感染峰值和再生数。

SIR模型常见问题

beta和gamma分别代表什么?

在这个简化模型中,beta控制可能引发感染的接触速率。gamma是每日康复率,其倒数近似于平均传染期。

R0大于一意味着什么?

在模型假设下,每名感染者最初平均引发超过一例新感染。这表示存在增长潜力,并非确定的预测。

为什么三类人数之和必须等于总人口?

封闭的SIR人群中,每个人必须且只能属于一类。总数不符会凭空增加或遗漏人数,违反方程中的守恒关系。

模型是否包含疫苗接种或再次感染?

这个基础版本既不包含疫苗接种导致的人群转移,也不包含免疫力丧失。更复杂的分室模型可以加入这些机制。

按天模拟的精度如何?

每日欧拉步长提供了便于理解的近似结果,适合学习和情景比较。科研或政策工作应采用经过校准的模型、更精细的数值方法及不确定性分析。