传染病模拟练习题与答案

传染病模拟练习题与答案某沿海城市A市人口约200万，2023年11月出现一种新型呼吸道病毒HRV-23疫情。经疾控中心监测，该病毒潜伏期为2-5天（平均3.5天），传染期为出现症状前1天至症状消失后2天（平均传染期5天）。传播途径以飞沫传播为主，气溶胶传播为辅，基本再生数R0初期为3.2。人群普遍易感，无天然免疫力。当地医疗机构日均可处理重症患者200例，重症率约为感染总数的5%（重症定义为需住院治疗的肺炎或呼吸衰竭）。截至12月1日，A市报告首发病例（输入性病例），12月3日报告第2例本地感染病例，12月5日累计报告10例确诊病例（均为轻症，无重症）。假设：①所有病例均被及时发现并隔离（隔离后无传染性）；②未采取其他防控措施时，每个感染者日均接触25人，接触者感染概率为4%；③忽略人口流动与死亡；④初始易感者数量等于总人口（200万）。问题1：根据12月1日至12月5日的疫情数据，计算该阶段病毒的实际再生数Rt，并分析其与R0存在差异的可能原因（需列出计算过程）。问题2：若A市未采取任何防控措施（即不隔离病例），假设R0保持3.2不变，使用SIR模型（易感者S、感染者I、康复者R）模拟疫情发展。已知初始时S0=2,000,000，I0=1（12月1日首发病例），R0=0，传染率β=R0/(平均传染期×S0)，恢复率γ=1/平均传染期（单位：天⁻¹）。请：（1）写出SIR模型的微分方程组；（2）计算β和γ的具体数值（保留4位小数）；（3）估算疫情达到峰值时的感染者数量（提示：峰值时dI/dt=0，且S=γ/β）。问题3：A市于12月10日启动“3天居家健康监测”政策（即发现病例后，其密切接触者需居家监测3天，期间接触人数降至日均5人，感染概率降至1%）。假设12月10日当天有100例感染者（均未被隔离），密切接触者数量为每个感染者的前3天接触人群（即每个感染者在传染期内接触的25人/天×3天=75人）。计算该政策实施后，这100例感染者在剩余传染期内（假设剩余传染期为2天）可能引发的续发病例数（需考虑政策前后的传播差异）。问题4：A市计划通过疫苗接种降低疫情风险，疫苗保护率为80%（即接种者感染概率降至20%）。若目标是将有效再生数Rt控制在1以下（Rt=R0×(1-疫苗接种率×保护率)），计算需要达到的最低疫苗接种率（保留2位小数）。问题5：12月15日，A市累计报告500例确诊病例，其中重症25例（符合5%重症率）。假设重症患者需住院10天，且12月15日后新增感染者中重症率保持5%，日均新增感染者数量为200例（未采取新措施时）。若医院重症床位仅150张，请评估12月20日（5天后）是否会出现重症床位短缺，并说明理由（需计算具体数值）。答案1：实际再生数Rt的计算需基于传播代际数据。首发病例（12月1日）为第1代，第2例本地病例（12月3日）为第2代（间隔2天），12月5日累计10例中，除首发病例和第2例，剩余8例为第3代（假设第2代病例在12月3日出现后，传染期为5天，可传播至12月8日，但12月5日已报告10例，说明第3代病例集中在12月4-5日）。传播代际时间（serialinterval）取平均潜伏期3.5天（因二代病例多在一代病例出现症状后潜伏期内感染）。截至12月5日，第1代1例，第2代1例，第3代8例（10-1-1=8）。Rt=（第n代病例数）/（第n-1代病例数）。第2代Rt=1/1=1；第3代Rt=8/1=8（但需考虑时间间隔）。更准确的计算应使用指数增长模型：新增病例数=I0×e^(rt)，其中t为时间（天）。12月1日I=1，12月5日I=10，t=4天（12月1日到12月5日共4天），则10=1×e^(4r)，解得r=ln(10)/4≈0.5756天⁻¹。Rt=r×平均传染期+1（近似公式），平均传染期5天，故Rt≈0.5756×5+1≈3.88（或用Rt=1+r×T，T为传染期，更准确的公式为Rt=e^(rT)，则Rt=e^(0.5756×5)=e^2.878≈17.7，显然不合理，说明早期数据波动大，应采用代际法）。更合理的方法是：首发病例（12月1日）传染期为12月1日-12月6日（症状前1天至症状消失后2天，假设症状出现在12月3日，则传染期为12月2日-12月5日，共4天），实际在12月3日前接触的人可能感染第2例（12月3日），12月3-5日接触的人感染后续8例（12月5日报告）。因此，首发病例在传染期内引发1（第2例）+8（第3代）=9例，Rt=9/1=9。但初期因病例少、防控措施（如病例被隔离）导致实际传播受限，而R0=3.2是无干预时的理论值，差异原因可能包括：①首发病例被隔离后传染期缩短（实际传染期可能小于5天）；②密切接触者被追踪并采取防护措施（如戴口罩）降低感染概率；③早期病例症状轻，传播力弱于平均水平。答案2：（1）SIR模型微分方程组：dS/dt=-βSIdI/dt=βSIγIdR/dt=γI（2）R0=βS0/γ（因初始时S≈S0），故β=R0×γ/S0。已知平均传染期=1/γ=5天，故γ=1/5=0.2天⁻¹。代入R0=3.2，S0=2,000,000，得β=3.2×0.2/2,000,000=0.64/2,000,000=3.2×10⁻⁷人⁻¹·天⁻¹（保留4位小数为3.2000×10⁻⁷）。（3）疫情峰值时dI/dt=0，即βSI=γI→S=γ/β。代入γ=0.2，β=3.2×10⁻⁷，得S=0.2/(3.2×10⁻⁷)=625,000人。此时易感者剩余625,000人，感染过的人数为S0-S=2,000,000-625,000=1,375,000人。根据SIR模型，峰值时感染者数量I=(S0SR)，但R=γ∫Idt，近似可用I峰值≈(R01)/R0×S0（当R0>1时），但更准确的公式为I峰值=S0(γ/β)(R01)/R0×(γ/β)（推导复杂）。另一种方法是，当S=γ/β=625,000时，总感染人数约为S0S=1,375,000，而峰值时I≈(βS0/γ1)×γ/β（近似），代入得I≈(3.21)×625,000/3.2≈2.2×195,312.5≈429,687人（实际需数值模拟，此处为估算）。答案3：政策实施前（12月10日前），每个感染者日均接触25人，感染概率4%，剩余传染期假设为2天（原传染期5天，已过去3天）。但政策于12月10日启动，密切接触者为前3天接触的75人（25人/天×3天），其中：前1天（12月9日）接触的25人：未被政策影响，感染概率4%，传染期剩余2天（12月10-11日），但12月10日已启动政策，可能这部分人已被纳入监测，假设政策对所有密切接触者生效。政策实施后（12月10-11日），每个感染者剩余2天传染期，接触人数降至5人/天，感染概率1%。计算续发病例数：（1）前3天密切接触者（75人）：假设这些人在被监测期间（3天）不会再传播，但需计算他们被当前100例感染者感染的概率。每个感染者接触75人，感染概率为：政策前接触的25人/天×3天=75人，其中前2天（12月8-9日）接触的50人可能在政策前已感染（但题目假设政策12月10日启动，故前3天接触的75人包括12月7-9日接触的人群，此时尚未隔离，感染概率为4%）。但更准确的是，每个感染者在传染期内的传播分为两部分：政策实施前已传播的部分（12月10日前）和政策实施后剩余的2天。假设100例感染者在12月10日时，已处于传染期的第3天（传染期共5天，已过去3天，剩余2天），则：前3天（12月7-9日）：每天接触25人，感染概率4%，已传播的病例数=100×3天×25人/天×4%=100×3×25×0.04=300例（但这些病例可能已被报告或处于潜伏期）。剩余2天（12月10-11日）：接触5人/天，感染概率1%，续发病例数=100×2天×5人/天×1%=100×2×5×0.01=10例。因此，总续发病例数≈前3天已传播的300例（但可能已被隔离）+剩余2天的10例=10例（因前3天接触者可能已被纳入监测，实际新增为政策后的10例）。答案4：有效再生数Rt=R0×(1-疫苗接种率×保护率)≤1。已知R0=3.2，保护率=80%=0.8，设接种率为p，则：3.2×(10.8p)≤110.8p≤1/3.2≈0.31250.8p≥10.3125=0.6875p≥0.6875/0.8≈0.8594，即85.94%。因此，最低疫苗接种率需达到85.94%。答案5：12月15日重症患者25例，占用床位25张（住院10天，故12月15日入院的25例需占用床位至12月25日）。12月15日后，日均新增感染者200例，重症率5%，即日均新增重症=200×5%=10例。12月15-20日共5天，新增重症=10×5=50例，累计重症=25+50=75例。但需考虑已住院的重症患者是否出院：12月15日入院的25例，住院10天，出院时间为12月25日，故12月20日时仍未出院；12月16日入院的10例（12月16日新增重症），出院时间为12月26日，12月20日仍住院；同理，12月17日入院的10例（出院12月27日）、12月18日（出院12月28日）、12月19日（出院12月29日）、12月20日（出