基于SIQR的新冠肺炎疫情动态传播模型研究.pdf
《基于SIQR的新冠肺炎疫情动态传播模型研究.pdf》由会员分享,可在线阅读,更多相关《基于SIQR的新冠肺炎疫情动态传播模型研究.pdf(5页珍藏版)》请在咨信网上搜索。
1、第 卷第期湖北工业大学学报 年 月V o l N o J o u r n a l o fH u b e iU n i v e r s i t yo fT e c h n o l o g yA u g 收稿日期 第一作者钞寅康(),男,湖北襄阳人,湖北工业大学硕士研究生,研究方向为系统仿真与设计.通信作者龚立雄(),男,湖北仙桃人,工学博士,湖北工业大学副教授,研究方向为生产建模.文章编号 ()基于S I Q R的新冠肺炎疫情动态传播模型研究钞寅康,龚立雄,雷彬文(湖北工业大学机械工程学院,湖北 武汉 )摘要以武汉新冠疫情为案例,运用系统动力学理论,基于传统S I R模型进行改进,增加了潜伏期、
2、感染率、规避系数、隔离率、治愈率、死亡率等参数,建立有限区域内的新冠肺炎疫情传播演化S I Q R系统动力学模型,研究和分析疫情传播路径以及传播影响因素.参考相关政府部门对武汉的疫情采取了强力控制措施,有效阻断疫情传播等实际情况,对新冠疫情的传播演化机理进行分析和验证,证明了本文S I Q R模型的优越性.通过多因素影响分析及讨论,结果表明:规避系数大于时能明显遏制疫情扩散;病毒潜伏期的增加会提高疫情不可控性;隔离率的增加对疫情有抑制作用但同时会增加疫情防控时长;疫苗覆盖率达到 以上能形成人群免疫屏障.最后在实验分析基础上,提出新冠疫情的预防及控制方案,对于重大公共疫情的控制与防治具有重大意义
3、.关键词S I Q R模型;传播演化;规避系数;疫苗覆盖率;系统动力学 中图分类号T P 文献标识码A新型冠状病毒肺炎是一种传染性很强通过呼吸道传染的疾病,基于新冠疫情的发展趋势,目前主要从两方面进行专业领域的研究,一种从疫情的控制角度,主要针对病毒的平均潜伏期以及基本再生数 ;另一种则是基于S I R模型以疫情传播的建模和预测为主 .其中,阿曼多等利用S I R模型预测新冠肺炎疫情的演变,预测的数据与具体病例数据有很好的拟合.任中贵提出了基于S I R模型的C OV I D 疫情数据分析,其中感染率和移除率这两个重要参数的拟合结果说明及时采取封城等防疫措施对医疗救助能力的提升十分重要.或者采
4、用S E I R模型进行数据拟合 ,然后进行适当的假设来推测疫情后期发展趋势,为疫情防控提供参考.S E I R模型相当于在S I R模型的基础上增加了传染病的潜伏期.赵成珍在文章中对S E I R模型改进为三阶段模型然后对现实数据进行拟合;丁振华提出了一种新的新冠肺炎传播非线性动力学模型;刘建国 在基于空间交互模型的基础上同时构建S E I R模型;董章功 提出了一种基于传统的传染病动力模型S E I R和差分整合移动平均自回归模型A R I MA构建的S E I R A R I MA混合模型;以上模型都在S E I R模型基础上有所改进.但是S I R模型和S E I R模型二者均不能完全
5、适用于湖北武汉的疫情发展的实际情况,考虑到 年月 日采取封城措施,对感染者进行隔离,在模型中增加隔离者提高模型的仿真水平,提高仿真结果与实际疫情的拟合度,给疫情防控提供更有价值的参考.在上述背景下,本文基于系统动力学分析,构建了新冠肺炎疫情扩散的S I R模型,并在传统的S I R模型上进行改进和优化,改进为S I Q R模型.最后,以武汉疫情为数据样本,使用改进后的S I Q R模型进行仿真建模实验,并通过实验寻求预防和控制新冠疫情传播的方法.结果表明:该模型能有效地复现新冠病毒传播演化过程,为重大公共疫情的预测与防控提供指导依据.S I Q R模型与方法S I Q R模型构建本次疫情传播迅
6、速与传染病研究中S I R模型的特点相类似,但S I R模型在此次疫情中并不完全适用.以往的疫情没有如此果断且迅速的采取大规模隔离措施,因此本文在S I R模型的基础上增加隔离人群进行改进,使用S I Q R模型进行有限区域内的疫情传播演化仿真.改进后的S I Q R模型中,治愈的传染个体拥有免疫抗体,设传染率(感染系数)为,恢复率(恢复系数)为,S I Q R模型状态转换图如图所示.S I Q R模型将有限区域内的人口分为以下四类:)染病者i(t),t时刻已被感染且具有传染力的人数;)易感者s(t),t时刻还未染病但可能被传染SIRQ1Q2?2?1?21-?1?2?图S I Q R模型状态转
7、换图的人数;)隔离者q(t),t时刻已被进行隔离措施的人数;)恢复者r(t),t时刻已从染病者中移出的人数.设总人口为N(t),则有N(t)s(t)i(t)q(t)r(t).S I Q R模型假设S I Q R模型基于以下四个假设:)假设感染者的传染能力较强,但凡与易感染者接 触 就 有 被 感 染 的 风 险.假 设 易 感 染 总 数 为s(t),在t时刻单位时间内,被病人感染的比例系数为,那么s(t)i(t)就表示在t时刻内被病人传染的人数.)假设人口基数不变,忽略人口的流动、死亡和出生等种群动力因素.即N(t)K.)假设易感者和感染者分别按一定概率进行隔 离,在t时 刻 进 行 隔 离
8、 的 人 数 为s(t)i(t).)假设病人数量与从染病者中移出的人数在t单位时间内成正相关,且相关系数为,则单位时间内移出者的数量为 i(t)q(t).基于以上四个假设条件,感染机制如S(i)Q(i)I(i)S(j)Q(j)I(j)I(i)Q(i)R(i)()由以上所有假设可知:若对感染者和易感者采取一定的隔离措施,感染者的恢复率增长为 i(t),易感个体的下降率为i(t)s(t),感染者的增长率达到i(t)s(t)i(t)q(t)i(t).易感者从患病到移出的过程可用微分方程表示,如dS(t)dt I(t)S(t)NS(t),dI(t)dt I(t)S(t)NI(t)Q(t)N()I(t)
9、,dQ(t)dtS(t)I(t)Q(t)N,dQ(t)dtS(t)Q(t)Q(t),dR(t)dt I(t)Q(t)()S I Q R模型参数设定及应用模型参数设定本案例基于系统动力学进行S I Q R新冠疫情传播演化仿真,时间以离散数据的形式表示,分别使用易感者、感染者、隔离者和治愈者来区分个体的不同状态.隔离者又分为在家隔离者和医院收治隔离者两类群体.移出者分为感染死亡/隔离死亡/治愈者.感染者区分为入院治疗和重症隔离两种情况,两者将不同程度地影响死亡人数和治愈人数.模型公式及参数设定如表.表S I Q R新冠疫情传播演化模型的相关计算公式感染易感染者(感染者潜伏者接触感染率接触人数隔离者
10、隔离者患病率患者接触人次)/总人口易感染者隔离易感染者隔离率感染者隔离感染者感染者隔离率隔离者患病隔离者隔离者患病率/潜伏期隔离者未患病隔离者(隔离者患病率)隔离者治愈隔离者(隔离者死亡率)隔离者死亡隔离者隔离者死亡率感染者死亡感染者感染者死亡率感染者治愈感染者(感染者死亡率)S I Q R模型的验证本文数据来源于中国湖北省卫生健康委员会官网(h t t p:w i wh u b e i g o v c n)公布的湖北省 年 月 日 年月的疫情数据进行实验.然后基于S I Q R模型对新冠疫情传播演化机理进行仿真和推演,可以得到封城前后的数据参数如表所示.表封城前后数据/(封城前)/(封城后)
11、总人口 接触人数 患者接触感染率 潜伏者接触感染率 患病周期潜伏期规避系数隔离率院内治愈率 院外治愈率相比于传统S I R模型,增加了潜伏期、规避系数、隔离率等参数.经仿真求解,得出封城前和封城后治愈人数、感染人数和死亡人数三者的曲线,如图、图所示.通过图、图的分析,可以将我国新冠疫情累计确诊数量的拟合情况大致分为三个阶段:)第一阶段为 月 日月 日,疫情初期防控力度小,传染途径广,感染人数呈指数增长,并产生了大量的二、三代感染者,但此时医院资源供给相对充足,感染者治愈率较高.)第二阶段月 日月日,随着防护意湖北工业大学学报 年第期05101520253.02.52.01.51.00.5?10
12、4?图封城前曲线010203040506070809087654321?104?图封城后曲线识及政府管控力度增加,聚集感染事件减少,治愈系数已经大于感染系数,感染人数达到峰值后迅速下降,但由于前期大量感染者的涌入造成医疗资源陷入匮乏,导致治愈者人率并未达到峰值,且死亡人数还在增长.)第三阶段月日后,封城物理阻断了传播途径,治愈人数呈指数增长,感染者人数逐渐趋近于零,死亡人数曲线已趋于平缓.本文模型中考虑到武汉及时的“封城”行为,以及政府在多方面实施强有力措施,使新冠肺炎传播在时间、空间上具有明显的异质性,而本次仿真中得到的指数增长,可以推论是因为病毒平均潜伏期较长,在有限区域内对特定人群(中老
13、年人等免疫力较差的群体)的感染性强,导致疫情爆发速度快,医疗体系在短时间内不堪重负、瘫痪所造成的.多因素的影响研究与讨论为了更加贴近武汉疫情传播真实现状,进行有效模型的构建,本文在模型构建过程中基于政府以及群众不同的角度在疫情发展的不同阶段,采取的不同措施进行分析.通过添加规避系数以及隔离率两大因素,主要体现政府采取措施的力度以及群众防疫意识的强弱;通过添加病毒潜伏期和疫苗覆盖率两大因素,主要从病毒分析,体现病毒对医疗以及社会经济的影响.规避系数对疫情传播演化的影响规避系数的定义是将该地区各方面防疫措施对新冠肺炎疫情扩散程度的影响引入模型,是对该地区的政府是否采取措施,措施力度是否严厉,以及医
14、疗资源,人口素养等多方面因素的综合.0102030405060708090100500000040000003000000200000010000000.10.50.9?/图封城后规避系数不同对感染人数的影响武汉封城举措,影响规避系数,图中规避系数越大,感染者人数越少.其中分为三个小阶段:):视为该地区没有采取防疫风险等级制度的相关措施,区域内人群活动自由,频繁发生交叉感染,隔离点缺乏,患者无法及时得到治疗和隔离.疫情管控呈现出失控的态势.):该地区采取了对应防疫风险等级制度的相关措施,内防扩散、外防输入,在疫情爆发过程中及时隔离了二、三代感染者以及其密切接触者,并及时保障了足够的医疗资源.)
15、:该地区在发现疫情时及时进行了严格的隔离措施,人口都有较高的防控素养,并配有足够的医疗资源保障,疫情在该地区几乎很难爆发.隔离率对疫情传播演化的影响本模型中,隔离率不同时对感染者数量的影响如图所示.0102030405060708090100500000040000003000000200000010000000.20.40.8?/%?/图隔离率不同时对感染者数量的影响由图可知,隔离率对疫情扩散程度的影响至关重要,隔离率越高,感染者人数越少,隔离率的提高能有效遏制疫情传播,但会给医疗体系带来巨大负担,从而造成较晚的峰点到来和较长的疫情防控时间,给社会经济发展和人们健康造成沉重打击.平均潜伏期对
16、疫情传播演化的影响本模型在计算感染者人数的参数时,将本次疫情备受关注的新冠肺炎潜伏期引入模型,发现当病毒在区域内人群当中平均潜伏期不同时,会造成不同的疫情扩散情况,如图所示.由图分析可知,当其他参数相同时,病毒潜伏期越短,感染者数量峰值越低,但峰点到来较慢,疫情持续时间更长,对地区影响更持久;当病毒潜伏期越长,感染者数量峰值越高,但峰点到来较快,疫情持续时间更短,对地区造成短时间内较大的冲击性第 卷第期钞寅康,等基于S I Q R的新冠肺炎疫情动态传播模型研究影响.0102030405060708090100500000040000003000000200000010000003114?/?/
17、d图潜伏期不同时对感染人数的影响新冠疫苗接种对疫情传播演化的影响研究疫苗覆盖率因素过程中,将已经接种疫苗设定为免疫人群,S I Q R模型中易感染者人数会除去已经接种疫苗人数.例如,当接种疫苗率为 时,易感染者人数在原先的值总人口感染者人数基础上除去接种疫苗的人数即 的总人口.易感染者人数(总人口感染者人数)总人口),规避系数变为.当疫苗覆盖率为、时,感染者,易感染者、治愈者的人数变化曲线如图.171319253137434955616773798591971031091151211271331396000500040003000200010000?0?20/%?40/%?6/%0?8/%0?
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 SIQR 肺炎 疫情 动态 传播 模型 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。