1. 引言
1.1. 研究背景
传染病在人群以及其它种群中的传播是一个很复杂的过程,利用动力学的方法建立传染病的数学模型,并通过数学模型对传染病进行定量与定性的分析和研究以判断某种传染病在某一地区是否会蔓延持续下去。1926年Kermack和Mckendrick [1] 研究了1665~1666年伦敦的黑死病以及1906年孟买的瘟疫,构建了具有重要意义的仓室模型。基于SIR的简单模型对各类流行病的特点建立了各种升级版的SIR模型,模型预言的结果在流行病预防和控制决策过程中也会给出重要的参考价值。1932年他们又在文献 [2] 提出了SIS模型,并提出了阈值理论,Hethcote [3] [4] 对这些工作进行了系统的总结。文献 [5] [6] 建立相应的传染病动力学模型,从模型对应的方程角度分类,建立的模型分别为常微分方程模型 [7] [8] 、随机微分方程模型 [9] 、偏微分方程模型 [10] 等。建立更符合实际的传染病模型能更好地预测传染病的发病趋势,通过对一个更符合现实情况的模型进行研究,文献 [11] 对不同传染病进行预测同时找到更合适的防控措施来控制疾病的蔓延,这将大大地减少人们的患病率或致死率等,对人类健康的发展具有重要的作用。
1.2. 研究现状及问题
Gumel等 [12] 认为传染病绝大多数传播是由无症状感染者造成的,因此传染病模型至少应包含无症状传播这一重要特征。Mohsen等 [13] 提出了一个考虑移民和隔离影响的传染病模型,研究表明正确实施隔离措施有助于控制疾病传播。Kamara等 [14] 建立了潜伏期人群传播的SEIR数学模型,证明了无病平衡点和地方病平衡点的局部稳定性和全局稳定性。Zamir等 [15] 考虑了传染病的感染和消除,并对所建立的模型的再生数进行了灵敏度分析,找到影响传播的最敏感的参数。到目前为止,对疫苗接种的研究也在陆续展开,占主导地位的研究方法也是两种:调查和建立数学模型。例如kate等人利用ODE模型 [16] 对比了5种年龄结构的优先次序策略,发现优先给20至49岁的成年人接种疫苗可最大限度地降低发病率,但在大多数情况下优先给60岁以上的老人接种疫苗可最大限度地降低死亡率。类似地,Brody等人 [17] 也通过仓室模型得出了优先给大于60岁的老年人接种疫苗能最大程度降低死亡率的结果。Shen等人 [18] 研究了接种疫苗在缓解疫情方面的有效性,还有一些研究则侧重于得到使传播得到控制的最佳疫苗覆盖率以及接种疫苗对再生数的影响。现如今,在流行病爆发时,主要的医疗控制手段有注射疫苗、隔离防治、药物治疗等。疾病爆发后,群众也会自发地采取少去人群密集的公共区域、购买增强对疾病抵抗力的药物等措施。但对于一些新的传染病或是以及存在的但又发生新的变异的传染病,传统的治疗方案可能作用有限,需要医疗专家找到新的方案,这期间耗费的时间我们也需要考虑。当今,注射疫苗是控制传染病最直接最有效的方法,但一般而言,一种新的传染病出现需要几个月的时间研发、量产新的疫苗。这些在我们研究疾病内在机理、采取控制措施、构建数学模型时都应该考虑进去(有些因素影响力较小,在构建模型时可以省去)。建立合适的模型,结合相关理论知识,传染病专家可以发现疾病的传播规律和发展趋势并依此找到更好的预防和控制方案,从而最大程度地减少人们的生命财产损失。
1.3. 研究意义
有部分传染病是属于单链RNA病毒,可以感染动物或人类。RNA病毒具有遗传不稳定性的特点,由单链RNA病毒引起的传染病的传播可以通过接触和飞沫、空气以及污染物等方式传播,其在人与人之间的主要传播方式为呼吸道飞沫和接触。此外,物体因接触病毒而被污染,然后间接传播给易感者。虽然被污染的物体表面有助于病毒的传播,但这并不被认为是病毒传播的主要方式。这种单链RNA病毒不仅给全球公共卫生治理带来威胁,而且对社会的正常生活产生重大影响。目前单链RNA病毒的变体对现有疫苗也构成了严重威胁。特别地,隔离在防止传染病的进一步传播方面发挥着特殊作用。采取包括隔离在内的一系列措施,能控制传播、降低该病的病死率。单链RNA病毒及其变种不仅给全球卫生体系带来了前所未有的压力,也从根本上挑战了对病毒和疾病的认识。需要通过建模相关科学研究和各方面进行控制来抑制传播,所以对此类传染病的研究仍需继续。
2. 相关基础知识
为方便阅读理解本章节对文中所涉及的专业名词及引文定理等做简要阐释。
2.1. 基本再生数的计算方法——下一代矩阵法
基本再生数 [19] (Basic Reproduction Number,通常缩写为R0)是流行病学中用来衡量传染病传播能力的重要指标之一。它表示在人群中,一个患者平均会传染多少个新的个体。具体来说,基本再生数描述了在一个完全易感人群中,单个感染个体在其传染期间(通常是感染后可传染的时间段)会传播给多少个新的易感个体。如果一个疾病的基本再生数大于1,那么这种疾病就有可能形成疫情并在人群中传播。相反,如果基本再生数小于1,那么疾病传播会逐渐减弱并最终消失。基本再生数的计算通常涉及疾病的传播速率、感染概率、潜伏期等因素。它对于评估疾病的传播情况、制定预防控制策略以及预测疫情发展趋势都具有重要意义。下一代矩阵法可以用来计算基本再生数。
下一代矩阵法 [20] (Next Generation Matrix Method)是一种流行病学中用于估计基本再生数(R0)的方法。它是基于传统的传染病传播模型的改进,主要由David Champredon和James Watmough在2015年提出。下一代矩阵法的基本思想是将传染病的传播过程分解成若干个组成部分,然后通过矩阵的运算来描述这些组成部分之间的转化关系,从而计算出基本再生数。这种方法相对于传统的微分方程模型更为直观和灵活,尤其适用于复杂的传染病模型。
具体来说,下一代矩阵法的步骤如下:
设系统为
代表第i个染病仓室的新感染率,
,
代表个体通过所有方式转移到i仓室的速率,
代表个体离开i仓室的速率。
设
是系统的无病平衡点,则有
为无病平衡点
处的导数
其中F和V是
矩阵,且
则
的谱半径是系统的基本再生数,即
。
通过这种方法,可以更清晰地理解传染病在人群中的传播过程,并且可以用于评估不同干预措施对基本再生数的影响,从而指导公共卫生策略的制定和调整。
2.2. 平衡点局部渐近稳定的判别方法——Hurwitz判据
Hurwitz判据(Hurwitz Criterion)是一种用于判断实系数多项式的根的位置的方法。该判据是基于多项式的系数来确定多项式根的位置,从而判断多项式是否具有稳定的实根。
对于一个实系数多项式
,Hurwitz判据可以用以下步骤进行:
1) 构造Hurwitz矩阵:构造一个
的矩阵,其中第一行为
,第二行为
,依此类推,直到第n行为
;
2) 计算主子式:计算Hurwitz矩阵的所有
主子式(subdeterminants),其中
;
3) 判断稳定性:如果所有主子式的值都大于零,则多项式的所有实根都是负的;如果有任何一个主子式的值小于等于零,则多项式至少有一个正实根。
Hurwitz判据的主要思想是利用多项式的系数来间接地判断多项式根的位置,而不需要直接求解根。这种方法在控制论、信号处理以及工程数学等领域中有着广泛的应用,特别是在判断系统稳定性方面。
3. 模型的建立
本文考虑RNA单链病毒引起的传染病,它们在爆发时,速度之快、传播范围之广并且有较长的潜伏期,而有关部门采取的可控措施为隔离控制和注射疫苗。患病者被发现症状后立刻接受治疗,对部分病情较为严重的病例采取隔离治疗的方案。对于新出现的传染病和一些新变异的传染病,疫苗从研发到上市一般需要6个月的时间。但对于一般传染病,其疫苗之前就已经上市,疾病爆发初期就可以投入使用。对于日常生活中一般的传染病,患者初期染病后不会立刻进入隔离控制阶段,以普通流感为例,患病初期症状较轻,服用一般抗病毒药物即可。若是病情加重引起发烧、肺炎等,可能会住院隔离。对这类传染病模型就有其局限性。而本文参考了Sanjoy Basu等人 [21] 的SLETIR模型,建立SLEIRVW模型。
其中S是易感人群,L是隔离人群(S中的一部分),E是病毒在潜伏期是的感染人群,I是感染人群,R是康复者,V是接种疫苗的人群,W是环境中所含病毒量,
是输入者,
是易感人群转移到隔离人群的比率,
是隔离人群转移到易感人群的比率,d是自然死亡率,
是有效接触率,e1是易感人群接种比例,e2是接种后感染但在潜伏期的比率,e3是接种后感染的概率,
是潜伏期到感染的变化率,
是感染但不传播的概率,
是感染治疗后康复的比率,A是天生免疫者,m是每个仓室病毒释放量,
是消灭率。该系统中不考虑因病死亡人口且有一定的初始人口,模型的传播流程图如下图1。
Figure 1. A warehouse diagram of an infectious disease model with environmental virus load and vaccine impact
图1. 一类带有环境病毒量和疫苗影响的传染病模型仓室图
因为
所以有
当
时,
为了不使一般性,可行域为
4. 平衡点和再生数的计算
通过计算得到,无病平衡点
,其中
。
下面用下一代矩阵算法得出基本再生数,首先有
,
将
,
分别关于E,I,L,V,W分别求导,并将无病平衡点
代入上式可以得到
,
继续计算可以得到
其特征值矩阵为
所以得到基本再生数
令系统右端为0得到
通过计算可以得到地方病平衡点,当
时,地方病平衡点
,
其中
,
,
设参数值
,
,
,
,
,
,
,
,
,
4.1. 无病平衡点的稳定性
定理1:当
时,系统的无病平衡点
在
内是局部渐近稳定的,当
时不稳定。
证:在无病平衡点处,系统的Jacobian矩阵为
矩阵
的特征方程为
显然得到
,
.其中计算得到
由上式可知
,
,
,
,
。当
时,满足以下条件
(i)
;(ii)
;(iii)
;
(iv)
。
则无病平衡点处的Jacobian矩阵的七个特征值均有负实部。通过Hurwitz判据得到系统的无病平衡点
在可行域
内是局部渐近稳定的。
4.2. 地方病平衡点稳定性
定理2:当
时,系统的地方病病平衡点
在
内是局部渐近稳定的,当
时不稳定。
证:在无病平衡点处,系统的Jacobian矩阵为
矩阵
的特征方程为
显然得到
,
。其中计算得到
,
,
由上式可知
,
,
,
,
。当
时,满足以下条件
(i)
;(ii)
;(iii)
;
(iv)
。
则地方病平衡点处的Jacobian矩阵的七个特征值均有负实部。通过Hurwitz判据得到系统的地方病平衡点
在可行域
内是局部渐近稳定的。
5. 总结
研究传染病模型具有多方面的现实意义,尤其是在流行病学、公共卫生和医学领域:
1) 预测疾病传播趋势:传染病模型可以帮助预测疾病在人群中的传播趋势,包括感染人数、疫情爆发时间和地点等,为政府和卫生部门制定相应的防控策略提供科学依据。2) 评估干预措施效果:通过模拟不同的干预措施(如隔离、口罩佩戴、疫苗接种等),可以评估这些措施对疫情传播的影响,从而指导决策者采取最有效的防控措施。3) 优化资源分配:传染病模型可以帮助优化医疗资源的分配,包括病床、医护人员和医疗设备等,以应对疫情期间可能出现的患者激增情况。4) 研究疾病传播机制:通过建立传染病模型,可以深入研究疾病的传播机制、潜伏期、传播途径等重要参数,有助于更好地理解疾病的发生和传播规律。5) 制定应对策略:传染病模型的研究结果可以为制定更有效的疫情应对策略提供科学依据,包括早期诊断、密切接触者追踪、社交隔离等措施。6) 预警和监测:建立传染病模型可以帮助建立疾病预警系统,及时监测疫情变化,提高对突发传染病的应对能力。本文分析了一类带有环境病毒量和疫苗影响的SLEIRVW传染病动力学模型,不仅考虑了疫苗对系统的影响,还将考虑了隔离仓室以及潜伏期人群,计算得到了模型的基本再生数
,并证明了当
时,系统的无病平衡点
在
内是局部渐近稳定的,当
时,系统的地方病病平衡点
在
内是局部渐近稳定的。从基本再生数的式子中可以看出降低接触率以及提高疫苗接种效果能使得基本再生数
减小,即能够有效防治该类传染病,这也为合理防控提供了建议。当然后续还要继续研究,通过最优控制策略和数值模拟,通过现实数据进行拟合与验证,以保证模型的真实性以及能更好的运用。总体来说,研究传染病模型对于预测和控制疾病传播、保障公众健康、优化医疗资源分配等方面具有重要的现实意义,可以有效应对传染病疫情带来的挑战。