# 传染病SIR模拟器：R0、疫苗接种、致死率与传播曲线分析
这个传染病SIR模拟器让你探索病原体如何在人群中传播-易感者被感染后，通过康复或死亡离开传染群体。它面向学生、科学传播者、公共卫生学习者，以及任何希望快速直观理解为什么传播或免疫方面的微小变化就能重塑疫情的人。交互式控制面板聚焦于人们最常想要测试的变量：R0、致死率、疫苗接种覆盖率、疫苗有效性、传染期和初始感染比例。图表实时更新，让你可以将易感者、感染者、康复者、免疫者和死亡曲线作为一个相互关联的疫情系统进行比较。# SIR模型的工作原理
基本SIR模型将人口分为三个主要舱室。S代表仍可能被感染的易感人群。I代表当前具有传染性、可以传播病原体的感染者。R代表因康复、获得免疫或其它原因离开传播链条、不再具有传染性的人群。该模拟器还从离开感染组的人群中追踪预估死亡数，作为重症结果的统计分支。传播速率与R0和传染期相关联。如果R0较高或人们传染时间较长，在感染组缩小之前会产生更多的新感染。如果疫苗接种从易感人群中移除足够多的人，有效再生数就会下降，疫情峰值可以变得小得多。| 控制参数 | 改变的内容 | 典型曲线效果 |
|---|---|---|
| R0 | 考虑免疫之前的传播潜能 | R0越高，感染曲线上升越快，峰值越高。 |
| 疫苗接种覆盖率 | 受保护后被移出易感人群的比例 | 覆盖率越高，Re越低，可以平缓疫情曲线。 |
| 疫苗有效性 | 在此简化模型中疫苗预防感染的程度 | 有效性越高，相同的覆盖率提供更强的保护。 |
| 传染期 | 人们保持传染性的平均时间 | 传染期越长，改变疫情时序，可能延长疫情。 |
| 致死率 | 离开感染组后被计为死亡的比例 | 致死率越高，死亡曲线上升，但不直接影响传播。 |
# R0、Re与群体免疫直觉
R0是一个感染者在完全易感人群中平均造成的二代病例数。有效再生数Re在部分人群已有免疫力、接种疫苗、隔离或其它原因无法被感染时会降低。在本模拟器中，有效疫苗接种直接降低易感性，因此展示的Re会随着受保护覆盖率的上升而下降。一个常用的群体免疫近似公式是1 - 1 / R0。对于R0为3的情况，阈值约为66.7%的有效免疫。该模拟器让这个阈值变得直观：当有效疫苗接种低于阈值时，疫情仍可能增长；当高于阈值时，传播难以维持。# 感染峰值意味着什么
感染峰值是模拟人群中同时处于传染期的最大人数。它在实际操作中通常比总感染人数更重要，因为医院、实验室、隔离项目和接触者追踪团队承受的压力来自同时存在的活跃病例。即使最终攻击率没有降为零，降低峰值也可能具有重要意义。攻击率估算的是整个模拟周期结束时的总感染比例。高攻击率意味着病原体在易感人群耗尽或被控制之前感染了大量人群。低攻击率意味着免疫、疫苗接种或弱传播阻止了广泛扩散。# 现实世界中的R0值及其对群体免疫的启示
基本再生数R0并不是某个病原体的固定生物学常数。它取决于接触模式、人口密度、文化习惯和环境因素。同一病毒在密集城市和农村地区，或在室内聚集更多的季节，可能具有不同的R0值。以下数值来自已发表研究的参考范围。| 病原体 | 典型R0范围 | 群体免疫阈值 (1 - 1/R0) | 主要传播特征 |
|---|---|---|---|
| 季节性流感 | 1.2 - 1.4 | 17% - 29% | 传染期短，具有季节性变化 |
| SARS-CoV-2 (原始株) | 2.0 - 3.0 | 50% - 67% | 症状前传播，气溶胶途径 |
| SARS-CoV-2 (Delta) | 5.0 - 8.0 | 80% - 87% | 病毒载量增加，传播更快 |
| SARS-CoV-2 (Omicron) | 8.0 - 12.0 | 87% - 92% | 免疫逃逸，上呼吸道嗜性 |
| 小儿麻痹症 | 5.0 - 7.0 | 80% - 86% | 粪口传播，长期无症状排毒 |
| 天花 | 5.0 - 7.0 | 80% - 86% | 通过全球疫苗接种运动被消灭 |
| 麻疹 | 12.0 - 18.0 | 92% - 94% | 传染性极强，空气传播，传染期长 |
| 百日咳 | 12.0 - 17.0 | 92% - 94% | 免疫力衰退导致重复感染 |
# 疫情过程中有效再生数Re如何变化
Re是疫情某一时刻的有效再生数。与假设完全易感人群的R0不同，Re考虑了免疫、疫苗接种以及任何降低易感性的因素。在此模拟器中，Re的计算公式为R0 x (1 - 受保护比例)，其中受保护比例是通过疫苗接种获得有效免疫的人口比例。模拟器头部显示的Re值会随着你调整控制参数而更新。当Re保持在1以上时，疫情持续增长。当它降至1以下时，每个感染者平均产生少于一个的新感染，疫情无法自我维持。这就是疫情控制的核心思路：通过免疫、行为改变或干预措施使Re降至并维持在1以下。# 攻击率、峰值负荷及其揭示的疫情严重程度
攻击率是在整个模拟疫情中总人口被感染的比例。这是疫情波过后最常引用的总结性指标。高攻击率意味着病原体在易感人群耗尽或控制措施阻止传播之前感染了大多数易感者。低攻击率意味着免疫、疫苗接种或传播能力本身较弱阻止了大规模感染。感染峰值-同时处于传染期的最大人数-对短期医疗压力更为重要。一个攻击率中等但峰值极高且尖锐的疫情波可能压垮医院，即使总病例数并不极端。相反，一条平缓的曲线可能在数周内达到类似的攻击率，给医疗系统留出处理病例的时间。这就是为什么公共卫生官员强调平缓曲线作为一个独立于预防所有感染的操作目标。# SIR模型中的平缓曲线与医疗容量
SIR模型中的感染曲线可以理解为同时需要医疗救助的人数。在真实疫情中，每个需要病床、氧气支持或重症监护的人都依赖于有限的资源池。当感染曲线超过可用容量时，由于系统无法提供足够的医疗服务，各种原因的死亡率都会上升。在此模型中，疫苗接种通过在被感染之前将人们从易感人群中移除来平缓曲线。通过其它措施-口罩、通风、社交距离、检测、隔离-降低R0也会在更完整的模型中降低峰值。这个模拟使机制变得可见：随着有效疫苗接种覆盖率的增加，峰值缩小、推迟或完全消失。# SIR模型背后的数学可视化
在SIR模型中，每个时间步长的新感染率取决于三个量：传播率beta、当前感染人数I以及仍处于易感状态的有效人口比例S / N。乘积beta x I x S / N被称为感染力。每一天，这个感染力度量了有多少易感者进入感染舱室。人们以康复率gamma = 1 / 传染期离开感染舱室。感染力和康复率之间的平衡决定了疫情是增长还是收缩。当beta x S / N超过gamma时，新感染超过康复，疫情扩大。当易感比例S / N下降到足够低时，gamma占主导地位，疫情收缩。参数beta在界面上并不直接可见。它通过关系式beta = R0 x gamma从R0和传染期推导得出。这就是为什么改变R0或传染期会产生相似但不完全相同的曲线形状。两个参数都影响感染力，但传染期还会拉长疫情的时间轴。# 如何在学习与教学中使用本模拟器
- 比较高R0与低R0场景：将R0设为1.5（季节性流感范围），然后设为6.0（疫苗接种前的小儿麻痹症范围）。注意即使所有其它控制参数相同，峰值高度、峰值时间和攻击率如何变化。
- 探索群体免疫阈值：从R0为3.0且无疫苗接种开始。记录攻击率。然后增加疫苗接种覆盖率直到Re降至1以下。比较覆盖率刚好低于和刚好高于阈值时的峰值和攻击率。
- 测试慢响应与快响应的效果：在不同水平设置疫苗接种覆盖率，观察峰值出现的时间。更高的覆盖率不仅降低峰值高度，通常还会推迟峰值出现，为医疗准备争取时间。
- 区分致死率与传播：将致死率设为0%，观察感染曲线。然后将致死率设为10%而不改变其它设置。感染曲线不变，但死亡人数上升。这展示了病死率和传播速度是截然不同的流行病学维度。
- 检查传染期的影响：在相同的R0下比较4天传染期与18天传染期。较长的传染期拉长曲线，推迟峰值，产生更长但更低的疫情波。
- 课堂练习-寻找阈值：让学生找出对于给定R0使Re降至1以下的最低疫苗接种覆盖率，然后将结果与公式1 - 1/R0进行比较。
# 何时以及为何使用本模拟器
- 流行病学学生：在使用微分方程或编程建模之前，将数学SIR框架与交互式曲线形状联系起来。
- 科学传播者与记者：生成图表、截图或实时解说，展示R0、疫苗接种和传染期为何对疫情轨迹至关重要。
- 公共卫生学习者：测试不同的干预组合如何改变疫情峰值和攻击率，培养对疫情应对中各种权衡的直觉。
- 任何对疫情数学感兴趣的人：无需编写代码或安装软件即可探索SIR模型。每个控制参数都会实时更新图表。