课件_SARS的传播学习.pptx

SARS 的传播(2003 年 A 题)第1页/共91页一、一、SARS 的传播的传播 SARS(严重急性呼吸道综合症，俗称：非典型肺炎)是 21 世纪第一个在世界范围内传播的传染病。SARS 的爆发和蔓延给我国的经济发展和人民生活带来了很大影响，我们从中得到了许多重要的经验和教训，认识到定量地研究传染病的传播规律、为预测和控制传染病蔓延创造条件的重要性。请你们对 SARS 的传播建立数学模型，具体要求如下：第2页/共91页 (1)对附件 1 所提供的一个早期的模型，评价其合理性和实用性。(2)建立你们自己的模型，说明为什么优于附件 1 中的模型；特别要说明怎样才能建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型，这样做的困难在哪里？对于卫生部门所采取的措施做出评论，如：提前或延后 5 天采取严格的隔离措施，对疫情传播所造成的影响做出估计。附件 2 提供的数据供参考。第3页/共91页 (3)收集 SARS 对经济某个方面影响的数据，建立相应的数学模型并进行预测。附件 3 提供的数据供参考。(4)给当地报刊写一篇通俗短文，说明建立传染病数学模型的重要性。以下只考虑建立 SARS 的传播模型。第4页/共91页二、微分方程二、微分方程 在研究一些涉及到变化规律的问题，特别是所研究的问题中涉及变量的变化率时，我们就可考虑微分方程模型。许多自然现象以及社会、经济、工程等领域中的问题，如传染病的蔓延，种群的相互竞争，经济增长的预测等，均可以通过微分方程模型来描述。第5页/共91页 所谓微分方程微分方程，就是表示未知函数、未知函数的导数与自变量之间的关系的方程。未知函数是一元函数的微分方程，称之为常微分方程常微分方程；未知函数是多元函数的微分方程，称之为偏微分方程偏微分方程。第6页/共91页 下面列出的等式都是微分方程。Malthus人口方程：虎克定律第7页/共91页 牛顿万有引力方程：波动方程第8页/共91页 热传导方程：势方程或 Laplace 方程 第9页/共91页 许多自然现象以及社会、经济、工程等领域中的问题，如地震破坏程度的估计，传染病的蔓延，种群的相互竞争，经济增长的预测等，其内在规律和发展趋势的描述，不能通过实验的方式来实现，必须通过机理分析的方法用微分方程模型来表示。第10页/共91页 这些问题又不像高等数学课程中所谓的微分方程应用题那样：假设条件给定，求解的结果就是问题的答案，并且答案是唯一确定的；而是要分析具体情况或进行类比才能给出假设条件，做出不同的假设，就会得到不同的方程(模型)；结论也不是确定的、唯一的，求解的结果还要用来解释实际现象并接受检验。第11页/共91页三、三、传染病传播模型传染病传播模型 人们不可能去做传染病传播的试验以获取数据，从医疗卫生部门得到的资料也是不完全和不充分的。不同类型的传染病的传播过程有其各自不同的特点，弄清这些特点需要相当多的病理知识，这里更不可能从医学的角度来分析各种传染病的传播，所以，我们只能按照一般的传播机理建立模型。第12页/共91页 传染病传播问题和自然科学中一些已经有确定规率的问题不同，不可能立即对它做出恰当的假设，建立完善的模型，只能先做出最简单的假设，建立模型，得出结果，分析是否符合实际，然后针对其不合理或不完善处，进行修改或补充假设，逐步得到较为合理的模型。第13页/共91页 模型模型 1(SI 模型)假设条件 (1)人群分为易感染者(Susceptible)和已感染者(Infective)两类，以下简称健康者和病人。时刻 t 这两类人在总人数中所占的比例分别记为s(t)和 i(t)。(2)在疾病传播期内所考察地区的总人数 N不变，既不考虑生死，也不考虑迁移，并且时间以天为计量单位。第14页/共91页 (3)每个病人每天有效接触的平均人数是常数，称为日接触率日接触率。当病人与健康者有效接触时，使健康者受感染变为病人。第15页/共91页 根据假设，每个病人每天可使 s(t)个健康者变为病人。因为病人数为 Ni(t)，所以每天共有 Ns(t)i(t)个健康者被感染，即病人数 Ni(t)的增加率为 Ns(t)i(t)。于是得到人员流程图如下第16页/共91页进而有 再设初始时刻(t=0)病人的比例为 i0，则由 s(t)+i(t)=1，得到初值问题 第17页/共91页 初值问题的解为 第18页/共91页 可画出 i(t)t 和 di/dt i 的图形为 i(t)t 的图形第19页/共91页di/dt i 的图形第20页/共91页 于是可知：当 t 时，i1，即所有人终将被传染，全变为病人(见下图及公式)。这显然不符合实际情况。其原因是模型中没有考虑到病人可以治愈，人群中的健康者只能变成病人，病人不会再变成健康者。第21页/共91页 然而，这个模型在传染病流行的前期还是可用的，可用它来预报传染病高潮的到来：当 i=1/2 时，di/dt 达到最大值(di/dt)m，这个时刻为这时病人增加得最快，可以认为是医院得门诊量最大的一天，预示着传染病高潮的到来，是医疗卫生部门关注的时刻(见 di/dt i 图)。第22页/共91页 还可以看出，tm 与 成反比。因为日接触率 表示给定地区的卫生水平，越小卫生水平越高，所以改善保健设施、提高卫生水平可以推迟传染病高潮的到来。第23页/共91页 模型模型 2(不考虑出生和死亡的 SIS 模型)有些传染病如伤风、痢疾等治愈后免疫力很低，可以假定无免疫性，于是病人被治愈后变成健康者，健康者还可以被感染再变成病人，所以在 SI 模型的基础上，增加一个假设条件就会得到 SIS 模型。第24页/共91页 假设条件 (1)人群分为易感染者(Susceptible)和已感染者(Infective)两类，以下简称健康者和病人。时刻 t 这两类人在总人数中所占的比例分别记为 s(t)和 i(t)。(2)在疾病传播期内所考察地区的总人数 N不变，既不考虑生死，也不考虑迁移，并且时间以天为计量单位。第25页/共91页 (3)每个病人每天有效接触的平均人数是常数，称为日接触率日接触率。当病人与健康者有效接触时，使健康者受感染变为病人。(4)每天被治愈的病人数占病人总数的比例为常数，称为日治愈率日治愈率。病人被治愈后称为仍可被感染的健康者，1/称为这种传染病的平均传染期平均传染期。第26页/共91页如果考虑到假设条件(4)，则人员流程图如下 于是有第27页/共91页 记初始时刻的病人的比例 i0(i0 0)，从而 SI 模型可以修正为我们称之为 Bernolli(贝努里)方程的初值问题，其解析解为第28页/共91页其中 =/。由 和 1/的含义可知，是整个传染期内每个病人有效接触的平均人数，称为接触数接触数。于是有第29页/共91页我们画出 di/dt i 和 i t 的图形为 di/dt i 的图形（1）第30页/共91页i(t)t 的图形（1）第31页/共91页di/dt i 的图形（1）第32页/共91页i(t)t 的图形（1）第33页/共91页 模型模型 3(考虑出生和死亡的 SIS 模型)当传染病的传播周期比较长时，若不考虑出生和死亡因素显然不妥，接下来考虑带有出生和死亡情况的 SIS 模型。假设条件 (1)人群分为易感染者(Susceptible)和已感染者(Infective)两类，以下简称健康者和病人。时刻 t 这两类人在总人数中所占的比例分别记为 s(t)和 i(t)。第34页/共91页 (2)在疾病传播期内所考察地区的总人数为N，总认为人口的出生率与死亡率相同，并且新生婴儿全为易感染者。记平均出生率为，则人口的平均寿命为 1/。(3)每个病人每天有效接触的平均人数是常数，称为日接触率日接触率。当病人与健康者有效接触时，使健康者受感染变为病人。第35页/共91页 (4)每天被治愈的病人数占病人总数的比例为常数，称为日治愈率日治愈率。病人被治愈后称为仍可被感染的健康者，1/称为这种传染病的平均传染期平均传染期。在上述的假设条件下，人员流程图如下第36页/共91页于是有第37页/共91页 记初始时刻的健康者和病人的比例分别是 s0(s0 0)和 i0(i0 0)，从而考虑出生和死亡的 SIS 模型为第38页/共91页而由 s+i=1 有 ds/dt=di/dt，于是，上式的第二个方程变为恒等式，从而模型简化为 如果令 =/(+)，则 仍表示整个传染期内每个病人有效接触的平均人数，即接触数接触数。于是，以下的求解与讨论与不考虑出生和死亡的 SIS 模型相同。第39页/共91页 模型模型 4(不考虑出生和死亡的 SIR 模型)许多传染病如天花、流感、肝炎、麻疹等治愈后均有很强的免疫力，所以病愈的人既非健康者(易感染者)，也非病人(已感染者)，它们已经退出传染系统。第40页/共91页 模型的假设条件为 (1)人群分为健康者、病人和病愈免疫的移出者移出者(Removed)三类，三类人在总人数 N 中占的比例分别为 s(t)，i(t)和 r(t)。(2)病人的日接触率为，日治愈率为，传染期接触数为 =/。(3)在疾病传播期内所考察地区的总人数 N不变，既不考虑生死，也不考虑迁移，并且时间以天为计量单位。第41页/共91页 在上述的假设条件下，人员流程图如下 于是有 s(t)+i(t)+r(t)=1第42页/共91页 记初始时刻的健康者和病人的比例分别是s0(s0 0)和 i0(i0 0)(不妨设移出者的初始值 r0=0)，于是得到 SIR 模型为如下的初值问题第43页/共91页而由 s+i+r=1 有 dr/dt=di/dt ds/dt，于是，上式的第三个方程变为恒等式，从而模型简化为 上述的初值问题无法求出解析解，只能通过数值解法求出数值解。第44页/共91页 例如，取 =1，=0.3，i(0)=0.02，s(0)=0.98，则求得数值解如下表。相应的 i(t)、s(t)曲线和 i s 曲线如下图。t012345678i(t)0.0200 0.0390 0.0732 0.1285 0.2033 0.27950.3312 0.3444 0.3247s(t)0.9800 0.9525 0.9019 0.8169 0.6927 0.54380.3995 0.2839 0.2027t91015202530354045i(t)0.2863 0.2418 0.0787 0.0223 0.0061 0.00170.0005 0.00010s(t)0.1493 0.11450.0543 0.0434 0.0408 0.04010.0399 0.0399 0.0398第45页/共91页SIR 模型的i(t)、s(t)曲线 第46页/共91页SIR 模型的 i s 曲线第47页/共91页 在实际应用 SIR 模型时，模型中的参数经常通过一些统计资料来估计。事实上，能够求出解析解的微分方程模型是非常有限的，所以人们经常利用定定性理论性理论从方程本身推出解的相关性质。对于上述的 SIR 模型，就可以采用相轨线分析相轨线分析的方法，来获得i(t)、s(t)的一般变化规律。(参教案，略)第48页/共91页 模型模型 5(考虑出生和死亡的 SIR 模型)模型的假设 (1)人群分为健康者、病人和病愈免疫的移出者(Removed)三类，三类人在总人数 N 中占的比例分别为 s(t)，i(t)和 r(t)。(2)病人的日接触率为，日治愈率为，传染期接触数为 =/。第49页/共91页 (3)在疾病传播期内所考察地区的总人数为N，总认为人口的出生率与死亡率相同，并且新生婴儿全为易感染者。记平均出生率为，则人口的平均寿命为 1/。在上述的假设条件下，人员流程图如下第50页/共91页此时由假设条件有s(t)+i(t)+r(t)=1 第51页/共91页 记初始时刻的健康者和病人的比例分别是s0(s0 0)和 i0(i0 0)(不妨设移出者的初始值 r0=0)，于是得到考虑出生和死亡的 SIR模型如下 第52页/共91页而由 s+i+r=1 有 dr/dt=di/dt ds/dt，于是，上式的第三个方程变为恒等式，从而模型简化为 第53页/共91页 采用相轨线分析，可以证明：若 1，则 i=0，s=1；若 1，则 i=ie，s=se，于是(ie,se)=(1/,(1)/)第54页/共91页四、四、SARS 的的传播模型传播模型 SARS 作为一种传染病，它的发生、传播、控制与消失有着传染病的一般规律。而它作为一种新发生的传染病，又具有特殊性传染性强、一旦感染就有生命危险。由于人们初期对 SARS 认识不足，导致传染范围过广，速度过快，这和以往的传染病不同。但我们仍可根据传染病建模原理，建立 SIR 及其推广模型。它们的基本思路相同，差异在于人口分类的多少，关键在于参数的确定。第55页/共91页 1.问题的分析问题的分析 建立 SARS 模型，需要我们考虑以下几个方面。了解 SARS 传播机理和传播状况，并且给出建模原理、方法、思路和框图。模型中的人群至少有三类：易感者 S、确诊病人 I 和退出者 R(包括治愈与死亡)，也可以再加入潜伏者 Z、隔离者 Q、疑似病人 L 等，要弄清楚他们之间的关系。第56页/共91页 模型应包含对于传传染染率率、治治愈愈率率和死死亡亡率率等重要概念的清晰表述。模型分析和计算中要给出上述参参数数的的估估计计方方法法和估估计计值值。还可包括平均治愈天数、隔离率和潜伏期等参数。模型的结果应该提供预测值(用数量或曲线来说明高峰期和持续时间)和隔离措施的效果(包括提前和推迟控制时间的影响，隔离人数多少的影响及遗漏病人的影响)。对于结果的分析应包括误差分析及模型和方法的通用性分析。第57页/共91页 2.基本假设基本假设 在了解 SARS 传播机理的基础上，不难做出如下的合理假设。单位时间内感染人数与现有的感染者成比例，称之为传染率传染率；单位时间内治愈人数与现有的感染者成比例，称之为治愈率治愈率；单位时间内死亡人数与现有的感染者成比例，称之为死亡率死亡率；第58页/共91页 SARS 患者治愈恢复后不再被感染；可以忽略各类人口的自然死亡率，忽略迁移的影响。上述的三个比例系数传染率、治愈率和死亡率是三个最关键的参数，它们可以是常数，也可以是时间或人口的函数。第59页/共91页 传染率的含义是每天每个 SARS 感染者传染的人数，其确定的原则是：当天新增 SARS病人人数除以当天 SARS 感染者人数，再进行曲线拟合即可。治愈率和死亡率是 SARS 患者每天治愈和死亡所占的比例，可以通过当天 SARS 感染治愈和死亡人数除以当天 SARS 感染人数，再进行曲线拟合即可。第60页/共91页 3.SARS 模型的模型的假设、参数说明及建立假设、参数说明及建立 3.1 假设假设 地区人口总数 N 可视为常数，流入人口等于流出人口。据人口所处的健康状态可将人群分为健康者、SARS 病人、疑似 SARS 病人、隐性感染者、退出者(被治愈者、免疫者和死亡者)等情况。在政府的强制措施下，人口基本不流动，故无病源的流入和流出，避免了交叉感染，降低了感染基数。第61页/共91页 隔离的人断绝了与外界的联系，不具有传染性。SARS 康复者二度感染的概率为 0。国家完善了监控手段，加强了对 SARS 病毒监控的力度，故可假设所有感染 SARS 病毒的人群都进入了 SARS 病人类和疑似类。由于对 SARS 病原体的研究不够深入，无有效药物可以使人体免疫，同时 SARS 病毒感染后，大量繁殖，破坏免疫系统，故不可免疫。第62页/共91页3.2 参数的设定和符号说明参数的设定和符号说明s(t):t 时刻健康者（不包括病愈和死亡）在总 体人群中的比例；i(t)：t 时刻SARS病人在总体人群中的比例；l(t)：t 时刻疑似病人在总体人群中的比例；r(t)：t 时刻被治愈者、死亡者和免疫者在总体 人群中的比例之和；z(t)：t 时刻潜伏病人在总体人群中的比例；第63页/共91页：SARS病人日接触率：每个病人每天有效接触（足以使健康者受感染变为病人）的平均人数.：日治愈率：每天被治愈的病人占病人总数的比例.：日转化率：每天危险群体中的疑似病人被确诊为SARS患者的比例.：每天危险群体中的疑似病人被确诊为SARS免疫者的比例.第64页/共91页：日死亡率：每天SARS病人死亡的数量和当天病人总数量的比值.：疑似感染率：每天健康者被感染为疑似病人的比例.：隐性感染者被确诊为病人占总隐性患者的比例.：隐性感染者被确诊为免疫者占总隐性患者的比例.：健康者转换为隐性感染者占总健康者的比例.第65页/共91页 3.3 SARS 传播传播模型模型 SARS 的的 SIR 模型模型 SARS 作为一种传染病，符合传染病的传播机理。当传染病处于初始阶段，人们对其认识不足，对病人没有进行分类(疑似、确诊)，也无控制措施，因此处于一种自然传播状态，可以将总人口分成三类：易感者、病人和退出者时，可直接套用前面的 SIR 模型处理问题。第66页/共91页 由假设，每个病人每天可使 1s(t)个健康者变为病人，因为病人人数为 Ni(t)，所以每天共有 1Ns(t)i(t)个健康者被感染，于是 1Nsi 就是病人数 Ni 的增加率。又因为每天治愈率为，死亡率为，所以每天有 Ni 个病人被治愈，有 Ni 个病人死亡。于是，病人的变化率为第67页/共91页 由于 s(t)+i(t)+r(t)=1，且对于移出者有故 SARS 传染病的 SIR 模型如下：第68页/共91页 SARS 的的 SILR 模型模型 随着人们对 SARS 病情的了解，政府采取强有力的预防措施，医院也对病人进行了分类(疑似、确诊)。故此时我们可将人群分为四类：易感者、病人、疑似病人及移出者，其人员框图和建立 SILR 模型如下：第69页/共91页第70页/共91页第71页/共91页 而由 s(t)+i(t)+l(t)+r(t)=1，有dr/dt=di/dt ds/dt dl/dt上式中的第四个方程变为恒等式，于是得到简化三维模型：第72页/共91页 SARS 的的 SILZR 模型模型 我们也可考虑易感者中包含处于潜伏期的病人的情况。此时我们将人群分为五类：易感者、潜伏期病人、确诊病人、疑似病人及退出者，其人员框图和建立的 SILZR 模型如下：第73页/共91页第74页/共91页第75页/共91页 而由 s(t)+i(t)+l(t)+r(t)+z(t)=1，有dr/dt=di/dt ds/dt dl/dt dz/dt，于是上述系统简化为如下的四维系统第76页/共91页 3.4 模型求解模型求解 参数的确定和分析参数的确定和分析 (1)参数、的确定第77页/共91页 利用 EXCEL 电子表格处理附件 2 中所给的数据得：=0.055076，=0.038183 =0.002443 注注：也可根据附件 2 中所给数据，画出上述各个参数的统计曲线，再根据统计曲线画出其对应的概率分布，取其概率平均值作为参数值。第78页/共91页 (2)参数 1 和 2 的确定 确定 1。显然，从我们建立的模型无法得到 s、i 的解析解，因为所得到的模型均为非线性的。为了解决这个问题，我们采用 MATLAB软件中“龙格库塔”方法求它们的数值解。具体确定方法如下：第79页/共91页 第一步，通过附件中给的实际统计数据算出每一天的 s 和 i，做出它们与时间的函数关系，如图 所示。图 根据实际数据拟合的 i(t)图象第80页/共91页 图 根据实际数据拟合的 i(t)图象第81页/共91页 第二步，对 1 取一组数，分别画出由通过龙格库塔法解出的模型数值解随时间变化的图象，如图 所示。第三步，将这组图象与由实际数据所得图象相比较、调试。我们发现当 1 1.5 时，理论图形与实际图形有最佳的吻合，如图 所示。第82页/共91页图 根据数值解作出 i(t)图像(1 1.5)第83页/共91页 从两个图形可看出，它们的高峰期、缓解期和平稳期曲线相当符合，具有相同的发展趋势。在 t0,10 的 SARS 初期范围内，曲线变化各不相同。这主要是因为在 4 月 20 日之前，没有相关数据的统计和报道。由于数据的不全，根据边界值画出来的曲线与通过数值解得到的曲线相比较，不能准确反映 SARS 发生初期时的趋势，所以边界值应该去掉，而通过数值解模拟的曲线可以得到之前的发展趋势。第84页/共91页 通过对 SARS 蔓延期特点的分析，如图，在符合所给数据反映的规律基础上，还能够模拟缺乏数据的 SARS 初始状态，所以曲线是合理的。第85页/共91页 确定 2。与确定 1 的情形类似，先根据实际数据画出图形，如图 所示。第86页/共91页 然后再对 2 取一组数，分别画出通过模型解出的数值解随时间变化的图象，将这组图象与由实际数据所得图象相比较、调试。发现当2 1.0 时，理论图形与实际图形有最佳的吻合。如图。第87页/共91页整个曲线反映了疑似患者在 SARS 的过程中的变化规律第88页/共91页 注注：(1)至于模型()和()中的其他参数,a,b,e 的确定类似于上述参数。(2)可以利用相平面、相轨线对结果进行分析与检验。(3)类似前面的传染病模型中内容，可对SARS 的传播过程进行分析，即各参数的各种变化，如何导致疾病的蔓延、高峰、缓解、控制等。(4)模型的改进，如随机偏微分方程组随机偏微分方程组模型等。第89页/共91页谢谢第90页/共91页感谢您的观看！第91页/共91页