• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    傳染病類突發(fā)公共衛(wèi)生事件風(fēng)險(xiǎn)評(píng)估與應(yīng)對(duì)

    2020-05-26 12:03:49趙成珍梁循王軍禮
    中國流通經(jīng)濟(jì) 2020年5期
    關(guān)鍵詞:突發(fā)公共衛(wèi)生事件持續(xù)時(shí)間

    趙成珍 梁循 王軍禮

    摘要:突發(fā)公共衛(wèi)生事件會(huì)沖擊社會(huì)經(jīng)濟(jì)運(yùn)行,特別是具備傳染性時(shí)這種沖擊的程度和范圍會(huì)更深更廣,而通過建立模型來模擬傳染病類突發(fā)公共衛(wèi)生事件,將為傳染病評(píng)估和應(yīng)對(duì)提供積極有效的依據(jù)。故以SEIR模型為基準(zhǔn),針對(duì)SEIR模型存在的不足,將重大傳染病事件分為三個(gè)階段發(fā)展三階段模型,將相關(guān)參數(shù)由靜態(tài)參數(shù)改進(jìn)為動(dòng)態(tài)參數(shù),分類模擬封城以及傳染病存在復(fù)發(fā)可能性的情形,并利用2020年新冠肺炎疫情數(shù)據(jù)進(jìn)行動(dòng)態(tài)模擬。對(duì)經(jīng)典SEIR模型以及改進(jìn)的三階段模型分別進(jìn)行參數(shù)估計(jì)和系統(tǒng)仿真發(fā)現(xiàn),在無人為干預(yù)的自然傳播情況下,傳染病感染數(shù)量峰值會(huì)出現(xiàn)在1~2個(gè)月后;加強(qiáng)隔離、加大醫(yī)療資源投入能顯著減小傳播規(guī)模;改進(jìn)的三階段模型能對(duì)傳染病類事件起到很好的預(yù)測(cè)作用,可為未來傳染病評(píng)估和應(yīng)對(duì)提供有效依據(jù)。由此,為有效應(yīng)對(duì)傳染病類突發(fā)公共衛(wèi)生事件,建議從加強(qiáng)隔離、加大醫(yī)療資源投入、疫情后期持續(xù)關(guān)注等角度著手完善我國傳染病應(yīng)急防控體系,在傳染病發(fā)生時(shí)能夠?qū)ο嚓P(guān)人員做到盡早隔離,提高治愈率,減少傳染數(shù)量,從源頭上降低突發(fā)公共衛(wèi)生事件對(duì)經(jīng)濟(jì)的負(fù)面影響。

    關(guān)鍵詞:傳染病模型;突發(fā)公共衛(wèi)生事件;沖擊程度;沖擊范圍;持續(xù)時(shí)間

    中圖分類號(hào):F064.1文獻(xiàn)標(biāo)識(shí)碼:A文章編號(hào):1007-8266(2020)05-0084-11

    收稿日期:2020-04-16

    基金項(xiàng)目:教育部人文社會(huì)科學(xué)項(xiàng)目“供給側(cè)改革、債務(wù)率調(diào)整與中國宏觀經(jīng)濟(jì)波動(dòng)”(19YJC790189);北京市教育委員會(huì)社會(huì)科學(xué)項(xiàng)目“創(chuàng)新和完善宏觀調(diào)控體系研究——基于供給側(cè)改革視角”(SM202010037006);北京物資學(xué)院校級(jí)項(xiàng)目“新型冠狀病毒肺炎傳染病模型構(gòu)建及在城市疫情防控中的應(yīng)用研究”(2020YQYJ12)

    作者簡介:趙成珍(1982—),男,山東省濱州市人,北京物資學(xué)院經(jīng)濟(jì)學(xué)院教師,博士,主要研究方向?yàn)榻?jīng)濟(jì)系統(tǒng)的數(shù)量建模等;梁循(1965—),男,北京市人,中國人民大學(xué)信息學(xué)院教授,博士,主要研究方向?yàn)樯鐣?huì)計(jì)算、人工智能、網(wǎng)絡(luò)金融;王軍禮(1978—),男,河南省平輿縣人,國務(wù)院發(fā)展研究中心公共管理與人力資源研究所研究員,北京大學(xué)能源經(jīng)濟(jì)與可持續(xù)發(fā)展研究中心研究員,理學(xué)博士,主要研究方向?yàn)榻?jīng)濟(jì)系統(tǒng)的量化分析。

    一、引言

    重大傳染病疫情、重大食物和職業(yè)中毒等突然發(fā)生、可能或已經(jīng)嚴(yán)重影響公眾健康的事件統(tǒng)稱突發(fā)公共衛(wèi)生事件。突發(fā)公共衛(wèi)生事件具有嚴(yán)重的社會(huì)經(jīng)濟(jì)危害性。在目前人口流動(dòng)速度不斷加快、生態(tài)環(huán)境不斷惡化的情況下,突發(fā)公共衛(wèi)生事件時(shí)有發(fā)生,給相關(guān)國家造成了一定的損失。及時(shí)有效處理各類突發(fā)公共衛(wèi)生事件,是維護(hù)國家和社會(huì)穩(wěn)定、保障人民生命安全的重要手段。

    在傳染病類突發(fā)公共衛(wèi)生事件中,由于傳染病具有很強(qiáng)的傳染性,且隨著經(jīng)濟(jì)和政治的全球化、交通的日益便利,人員流動(dòng)速度和商品貿(mào)易數(shù)量不斷增加,不僅加大了突發(fā)事件在區(qū)域之間的關(guān)聯(lián)度,而且加劇了傳染病的快速傳播。2003年的非典型肺炎(SARS)、2009年的H1N1流感、2014年的埃博拉以及2020年的新型冠狀病毒感染肺炎(以下簡稱“新冠肺炎”)等重大傳染病疫情暴發(fā)后,經(jīng)濟(jì)最先受到?jīng)_擊的就是旅游業(yè)、餐飲業(yè)、零售業(yè)及交通運(yùn)輸業(yè)。比如,根據(jù)國家統(tǒng)計(jì)局?jǐn)?shù)據(jù),在2003年5月非典型肺炎疫情的高峰期,我國內(nèi)地入境旅游人數(shù)同比下降31%,餐飲營業(yè)額同比下降15.5%,鐵路客運(yùn)量同比下降62.5%,社會(huì)消費(fèi)品零售總額增速降至4.3%,此外確診者及其密切接觸人群的治療和隔離也造成了勞動(dòng)力規(guī)模的下降。重大傳染病一旦暴發(fā),會(huì)威脅經(jīng)濟(jì)貿(mào)易的發(fā)展和人類的生命安全,其影響的范圍和嚴(yán)重性一般遠(yuǎn)高于其他類型的突發(fā)公共衛(wèi)生事件,因此世界各國一直都在關(guān)注傳染病的防治問題,我國在應(yīng)對(duì)傳染病方面也付出了巨大的努力。中華人民共和國成立之前,鼠疫、天花、霍亂、血吸蟲病等反復(fù)暴發(fā)。中華人民共和國成立之后,不斷制定和完善法律法規(guī),堅(jiān)持?jǐn)U大對(duì)公共衛(wèi)生事業(yè)的投入,取得了顯著成就,不僅消滅了霍亂、天花,戰(zhàn)勝了非典型肺炎等突發(fā)傳染病,流行傳染病的結(jié)構(gòu)也發(fā)生了變化,過去流行的傳染病發(fā)病率不斷下降。

    重大傳染病疫情發(fā)生時(shí),一個(gè)國家擁有一套完善的能夠有效管理突發(fā)公共衛(wèi)生事件的應(yīng)急防控體系是非常重要的。而相對(duì)于非典型肺炎疫情等其他突發(fā)公共衛(wèi)生事件,新冠肺炎是近年來我國乃至全球感染人數(shù)最多、波及范圍最廣、對(duì)經(jīng)濟(jì)影響更大的突發(fā)重大傳染病,在眾多突發(fā)公共衛(wèi)生事件中比較具有代表性。而且,對(duì)新冠肺炎疫情的應(yīng)對(duì)盡管體現(xiàn)了我國自2003年非典型肺炎疫情之后應(yīng)急防控體系極大的完善和進(jìn)步,但同時(shí)也暴露了我國應(yīng)急防控體系中依然存在的問題和不足。因此,本研究利用此次新冠肺炎疫情數(shù)據(jù),在傳統(tǒng)SEIR模型基礎(chǔ)上改進(jìn)傳染病動(dòng)力學(xué)模型,通過模型模擬疫情擴(kuò)散發(fā)展過程和趨勢(shì),既可為以后傳染病類突發(fā)公共衛(wèi)生事件防控措施的制定和效果評(píng)估提供科學(xué)依據(jù),也可為目前我國應(yīng)急防控體系存在的問題提供行之有效的建議,以便快速及時(shí)有效地應(yīng)對(duì)突發(fā)公共衛(wèi)生事件,避免問題持續(xù)發(fā)酵導(dǎo)致更嚴(yán)重的經(jīng)濟(jì)損失。

    二、文獻(xiàn)概述

    通過建模求解研究傳染病的歷史最早可以追溯到1760年,貝努利(Bernoulli D)[ 1 ]通過建立模型研究接種疫苗即牛痘對(duì)天花病毒傳播的影響。哈默(Hamer W H)[ 2 ]通過建立離散時(shí)間模型了解麻疹流行的復(fù)發(fā)情況,并創(chuàng)造性地假設(shè)了發(fā)病率,即每單位時(shí)間的新增病例數(shù)取決于易感者與確診者數(shù)量的乘積。羅斯(Ross R)[ 3 ]通過建立傳染病的微分方程模型,研究瘧疾在蚊蟲和人類間的傳播行為,結(jié)果發(fā)現(xiàn),若將蚊蟲數(shù)量控制在臨界值以下,瘧疾的流行將得到控制,羅斯還憑借對(duì)瘧疾在蚊子與人類間傳播動(dòng)力學(xué)的研究獲得1902年第二屆諾貝爾醫(yī)學(xué)獎(jiǎng)。隨后,其他的確定性傳染病模型也不斷被提出。20世紀(jì)初,數(shù)理流行病學(xué)逐漸興起,對(duì)傳染病模型的研究愈發(fā)深入。克馬克(Kermack W O)等[ 4 ]通過對(duì)黑死病等的研究,將被研究地區(qū)的人群分為三類,分別是易感者(Suscep? tible)、確診者(Infective)、康復(fù)者(Recoverer),提出了著名的SIR倉室模型,適用于描述通過病毒進(jìn)行傳播的傳染病。隨后,克馬克等[ 5 ]又在SIR倉室模型基礎(chǔ)上創(chuàng)建SIS倉室模型,提出了著名的閾值定理。但是,由于很多傳染病在發(fā)病前通常會(huì)有一定時(shí)間的潛伏期,比如艾滋病的潛伏期一般為2年到12年,狂犬病的潛伏期一般為1~3個(gè)月,多則可達(dá)幾十年,因此考慮到潛伏期的因素,后來的研究在SIR模型的基礎(chǔ)上不斷完善,增加了一個(gè)新的傳染病群體,即潛伏者(Exposed),提出了SEIR模型。

    近年來,國內(nèi)學(xué)者結(jié)合我國各時(shí)期傳染病數(shù)據(jù),多次應(yīng)用SEIR模型,并對(duì)之進(jìn)行改進(jìn)。徐恭賢等[ 6 ]通過建立SEIR模型對(duì)有潛伏期且治愈后終身免疫的非典型肺炎疫情進(jìn)行研究,實(shí)驗(yàn)數(shù)據(jù)表明,模型模擬結(jié)果與疫情實(shí)際情況比較吻合。王蓮花等[ 7 ]利用總?cè)丝诰哂蟹N族動(dòng)力學(xué)且具有一般傳染率的SEIR模型證明,地方病具有全局漸進(jìn)穩(wěn)定性。王樹忠等[ 8 ]研究了潛伏期和感染期均具有傳染性的SEIR模型,并最終得出結(jié)論,即對(duì)于潛伏期和感染期均具有傳染性的傳染病,只有控制潛伏者和確診者人數(shù),才能有效控制疫情傳播。吳昊澄等[ 9 ]擬合SEIR模型研究浙江省的學(xué)校流行性腮腺炎疫情,對(duì)比不同條件下應(yīng)急接種所達(dá)到的效果發(fā)現(xiàn),學(xué)校越早實(shí)行應(yīng)急措施,越早應(yīng)急接種,控制效果越好。崔景安等[ 10 ]根據(jù)手足口病具有潛伏期的特征,應(yīng)用SEIR模型預(yù)測(cè)北京市手足口病流行趨勢(shì),為手足口病預(yù)防控制提供參考。李冬梅等[ 11 ]建立的SEIR模型假設(shè)潛伏期具有常數(shù)輸入率,證明疾病模型僅存在地方病平衡點(diǎn),其基于我國2009年甲型流感病毒實(shí)際數(shù)據(jù)的模擬結(jié)果表明,嚴(yán)格控制處于潛伏期的病毒攜帶者進(jìn)入可有效控制和消除疫情。黃璜等[ 12 ]分別從自然傳播和采取不同隔離防控措施的角度構(gòu)建諾如病毒SEIR模型,為防止醫(yī)院暴發(fā)性感染諾如病毒提供了理論支持。王昕煒等[ 13 ]基于SEIR倉室模型,在潛伏者向確診者轉(zhuǎn)化過程中引入時(shí)滯參數(shù)刻畫潛伏期特性,發(fā)現(xiàn)在最優(yōu)疫苗接種策略情況下可有效控制傳染病傳播。

    總體來看,目前盡管SEIR模型已被應(yīng)用于很多流行傳染病的分析,但很少有實(shí)證研究探討影響面如此之廣的傳染病事件,且現(xiàn)有研究在估計(jì)參數(shù)時(shí)大多用整個(gè)過程估計(jì)出一個(gè)值,不符合現(xiàn)實(shí)不同階段的情況。本研究將重大傳染病事件分為三個(gè)階段,并對(duì)參數(shù)的估計(jì)加以改進(jìn),把死亡率和康復(fù)率作為時(shí)變量進(jìn)行疫情趨勢(shì)模擬,更符合現(xiàn)實(shí)情況,可為實(shí)施相關(guān)措施以解決影響社會(huì)經(jīng)濟(jì)發(fā)展的源頭問題提供思路。

    三、研究設(shè)定與模型構(gòu)建及模擬

    本研究以SEIR模型為基準(zhǔn),在此基礎(chǔ)上進(jìn)行改進(jìn)。本研究將傳染病流行范圍內(nèi)的人劃分為四類:一是易感者(記為St),代表t時(shí)刻未感染但有可能被感染的人群;二是潛伏者(記為Et),代表t時(shí)刻已經(jīng)被感染但尚未表現(xiàn)出發(fā)病癥狀且具有傳染能力的人群;三是確診者(記為It),代表t時(shí)刻已經(jīng)感染病毒并表現(xiàn)出相應(yīng)的發(fā)病癥狀,且具有傳染能力的人群;四是康復(fù)者(記為Rt),代表t時(shí)刻已經(jīng)康復(fù)且不會(huì)再被感染的人群,不再影響系統(tǒng)動(dòng)力學(xué)行為。

    (一)無人為干預(yù)時(shí)傳染病模型構(gòu)建與模擬

    1.模型假設(shè)

    借鑒傳統(tǒng)經(jīng)典模型的假定,加之傳染性疾病暴發(fā)時(shí)間較短,且從流行到結(jié)束一般不會(huì)持續(xù)很長時(shí)間,因此本研究進(jìn)行以下假設(shè):

    (1)假設(shè)總?cè)丝跒镹,保持不變,不考慮人口的遷入和遷出,忽略人口自然出生和死亡對(duì)N的影響,即:

    2.模型傳播機(jī)制

    無人為干預(yù)情況下的模型傳播機(jī)制參見圖1。

    為單獨(dú)觀測(cè)康復(fù)者人數(shù)變化,傳播機(jī)制中的R僅代表康復(fù)者。σ代表單位時(shí)間內(nèi)潛伏者轉(zhuǎn)化為確診者的比例,即潛伏期的倒數(shù)。γ代表確診者的單日康復(fù)率,簡稱康復(fù)率或治愈率。v代表確診者的單日死亡率,簡稱死亡率。僅進(jìn)行普通治療,不防控即自然傳播情況下的經(jīng)典SEIR模型微分動(dòng)力學(xué)方程為:

    3.參數(shù)設(shè)定

    本研究分析的傳染病流行范圍是以湖北省為例,因此經(jīng)典模型的參數(shù)估計(jì)基于湖北省疫情數(shù)據(jù)計(jì)算得到死亡率v和康復(fù)率γ。圖2展示了湖北省2020年1月24日至2020年3月30日的疫情情況。由圖2可以看出,2020年2月2日之前,由于樣本量較小,數(shù)據(jù)特征不穩(wěn)定,如單日死亡率和單日康復(fù)率波動(dòng)較大,但2月2日之后數(shù)據(jù)特征趨于穩(wěn)定。因此,選取2020年2月2日作為實(shí)證分析的起始點(diǎn),以在一定程度上減少評(píng)估誤差。

    2020年2月1日到2月2日疫情詳細(xì)數(shù)據(jù)參見表1。其中,單日死亡率指單日新增死亡人數(shù)與留存確診人數(shù)之比。比如,2020年2月2日單日死亡率的計(jì)算過程為:(294-249)/6 738=0.007。單日康復(fù)率的計(jì)算方法與之類似。

    在經(jīng)典模型中,死亡率v為由實(shí)際數(shù)據(jù)計(jì)算得到的平均死亡率0.007,康復(fù)率γ為平均康復(fù)率0.007。通過在萬德(Wind)金融終端檢索,得到2018年湖北省戶籍人口為6 172.9萬人。設(shè)N的初始值N0為6 172.9萬人,I的初始值I0為2020年2月2日的留存確診人數(shù)8 565人,R的初始值R0為2020年2月2日的康復(fù)人數(shù)215人。將I視為潛伏期為0天的人群,設(shè)潛伏期在0~14天平均分布,則潛伏者E的人數(shù)為確診者人數(shù)的14倍,即E的初始值E0估計(jì)為11.991萬人。易感者S的初始值S0= N0-E0-I0-R0,即6 160.031萬人。從理論上講,一名患者可以傳染2~3人,記為Cn,本研究取Cn為2,潛伏期平均為7天,潛伏者的非線性傳染系數(shù)β2估計(jì)為Cn/7,即0.286。由于確診者會(huì)伴有咳嗽及發(fā)熱等明顯癥狀,且一般會(huì)在醫(yī)院接受治療,密切接觸其他人群的機(jī)會(huì)遠(yuǎn)小于潛伏者,假設(shè)其接觸其他人群的機(jī)會(huì)為潛伏者的0.02倍,則確診者非線性傳染系數(shù)β1的估計(jì)值為0.006。自然傳播情況下模型主要參數(shù)設(shè)置參見表2。

    4.模擬分析

    現(xiàn)實(shí)中,在疫情發(fā)生初期,由于目前的傳染病預(yù)防機(jī)制存在一定時(shí)滯,比如需要對(duì)某種病毒是否存在人傳人的情況等進(jìn)行醫(yī)學(xué)上的論證,相關(guān)防控措施一般不會(huì)在疫情初期就達(dá)到極其嚴(yán)格的程度。因此,在發(fā)展初期疫情一般遵循自然傳播過程,此時(shí)人們尚未意識(shí)到該病的傳染性,不會(huì)采取嚴(yán)格的防控措施,醫(yī)療資源也不會(huì)因此增加。根據(jù)假設(shè)的自然傳播情況下的初始條件以及估計(jì)的參數(shù)值,本研究利用Matlab軟件對(duì)經(jīng)典SEIR模型的微分動(dòng)力學(xué)方程進(jìn)行數(shù)值求解,仿真得到新型冠狀病毒自然傳播的情況。發(fā)展趨勢(shì)具體參見圖3。

    從圖3的仿真結(jié)果可以看出,在自然傳播階段,如果不采取防控措施,且醫(yī)療處理水平維持不變,即康復(fù)率不變,新冠肺炎的蔓延十分迅速,且感染人群數(shù)量龐大,達(dá)到千萬級(jí)別。潛伏者數(shù)量在40天左右達(dá)到最高峰,之后開始下降;確診者數(shù)量在60天左右達(dá)到最高峰,之后開始下降,但要持續(xù)很長時(shí)間才能下降到較低的水平,一直到300天的時(shí)候還有200萬左右的確診者和2萬多名潛伏者。可以看出,如果放任傳染病發(fā)展,將導(dǎo)致勞動(dòng)力供給大幅下降且持續(xù)時(shí)間很長,給經(jīng)濟(jì)帶來長期嚴(yán)重負(fù)面影響,因此有必要分析防控措施能否有效阻止疫情蔓延。

    (二)積極防控條件下傳染病模型構(gòu)建與模擬

    1.模型假設(shè)與傳播機(jī)制

    假設(shè)對(duì)傳染病疫情進(jìn)行積極防控和治療等,如相關(guān)部門在采取一系列措施對(duì)易感者、潛伏者、確診者進(jìn)行不同程度隔離以降低傳染率的同時(shí),積極組織研究新的更有效的治療方案,加大醫(yī)療資源投入,如增加醫(yī)院數(shù)量、增加醫(yī)護(hù)人員數(shù)量、增加收治率、提高康復(fù)率,則傳染病經(jīng)典SEIR模型的傳播機(jī)制將變?yōu)槿鐖D4所示的傳播機(jī)制。

    其中,α1代表易感者和潛伏者的暴露程度,由對(duì)易感者和潛伏者流動(dòng)和接觸的限制程度決定。α2代表確診者的暴露程度,由醫(yī)護(hù)治療期間醫(yī)院內(nèi)部對(duì)確診者隔離治療的程度決定。α1和α2的值越低,代表隔離措施越嚴(yán)格。g代表醫(yī)療資源提升程度,收治率越高,醫(yī)療方案越有效,g值越大,即治療康復(fù)率越高,死亡率越低。根據(jù)積極防控條件下的疫情傳播機(jī)制,SEIR模型的微分動(dòng)力學(xué)方程可以表示為:

    2.模型結(jié)果分析

    本研究通過調(diào)整參數(shù)α1、α2、g的大小來模擬對(duì)易感者、潛伏者、確診者不同隔離程度以及不同醫(yī)療資源強(qiáng)度下疫情的發(fā)展趨勢(shì),采用科學(xué)方法評(píng)估防控措施的效果。當(dāng)參數(shù)α1分別取1、0.5、0.1時(shí),即針對(duì)易感者和潛伏者的隔離措施越來越嚴(yán)格時(shí),潛伏者、確診者數(shù)量的發(fā)展趨勢(shì)可分別由圖5、圖6看出;當(dāng)參數(shù)α2分別取1、0.5、0.1時(shí),即針對(duì)確診者的隔離措施越來越嚴(yán)格時(shí),潛伏者、確診者數(shù)量的發(fā)展趨勢(shì)可分別由圖7、圖8看出;當(dāng)參數(shù)g分別取1、5、10時(shí),即醫(yī)療資源投入越來越大時(shí),如增加醫(yī)護(hù)人員、增加床位、研發(fā)出特效藥等,潛伏者、確診者數(shù)量的發(fā)展趨勢(shì)可分別由圖9、圖10看出。

    在圖5和圖6中,α1代表對(duì)易感者和潛伏者的隔離程度,1代表不隔離,0.5代表隔離程度中等,0.1代表隔離程度高。由圖5可以看出,當(dāng)α1=1時(shí),代表自然傳播情況下潛伏者數(shù)量的發(fā)展趨勢(shì),潛伏者數(shù)量先經(jīng)歷40天的上升期,最高峰時(shí)超過1 000萬人,150天后才下降至16萬人;當(dāng)α1=0.5時(shí),即采取某些隔離措施使有效接觸率降低一半時(shí),潛伏者數(shù)量自隔離措施實(shí)施起開始下降,到51天時(shí)便降至1萬人;當(dāng)α1=0.1時(shí),即采取更加嚴(yán)厲的隔離措施時(shí),數(shù)量也是自隔離措施實(shí)施起開始下降,只需要18天,潛伏者數(shù)量便可降至1萬人。從三種情況下潛伏者數(shù)量發(fā)展趨勢(shì)的仿真結(jié)果可以看出,針對(duì)易感者和潛伏者的隔離措施越嚴(yán)格,潛伏者數(shù)量越能得到有效控制。也就是說,對(duì)易感者和潛伏者采取隔離措施可對(duì)整個(gè)疫情控制起到正向作用。

    由圖6可以看出,對(duì)易感者和潛伏者采取隔離措施,對(duì)控制確診者數(shù)量也起到了很好的作用。在自然傳播情況下,確診者數(shù)量在66天時(shí)達(dá)到最高峰,超過3 000萬人,且300天后才降至200萬人左右;如果采取某些隔離措施使易感者和潛伏者的有效接觸率降低一半,即α1=0.5時(shí),確診者數(shù)量在32天時(shí)便達(dá)到最高峰,僅17萬多人,并于300天后降至2.7萬人;如果對(duì)易感者和潛伏者采取比較嚴(yán)格的隔離措施,即α1=0.1時(shí),則在18天時(shí)就達(dá)到最高峰,僅10萬多人,并于300天后降至2 000多人。從確診者數(shù)量的發(fā)展趨勢(shì)也可以看出,對(duì)易感者和潛伏者采取隔離措施可對(duì)疫情控制起到比較明顯的作用。一方面,可使確診者數(shù)量更早地達(dá)到峰值,有效縮短疫情持續(xù)時(shí)間;另一方面,可有效降低確診者數(shù)量。但同時(shí)我們也發(fā)現(xiàn),在僅僅采取隔離措施的情況下,疫情后半階段仍然會(huì)持續(xù)較長時(shí)間。因此,對(duì)易感者和潛伏者采取隔離措施有助于縮小傳染病導(dǎo)致的勞動(dòng)力供給缺口,但經(jīng)濟(jì)受疫情影響的持續(xù)時(shí)間仍然較長。

    在圖7和圖8中,α2代表對(duì)確診者的隔離程度,1代表不隔離,0.5代表隔離程度中等,0.1代表隔離程度強(qiáng)。針對(duì)確診者采取不同程度的隔離措施,可對(duì)潛伏者和確診者數(shù)量控制起到一定作用,但效果不是特別明顯。這主要是因?yàn)椋环矫?,確診者數(shù)量相對(duì)于潛伏者和易感者數(shù)量來說還是很少的;另一方面,確診者一般都在醫(yī)院接受治療,這本身就已經(jīng)自然限制了其接觸大量人群的機(jī)會(huì)。因此,對(duì)確診者繼續(xù)加強(qiáng)隔離的防控效果不如對(duì)易感者和潛伏者繼續(xù)加強(qiáng)防控的效果明顯。但是,考慮到確診者傳染的人群以醫(yī)護(hù)人員為主,而醫(yī)護(hù)人員被傳染對(duì)醫(yī)療資源來說是一種損失,因此加強(qiáng)對(duì)確診者的隔離,從另外一個(gè)層面來看,是一種對(duì)醫(yī)療資源的加強(qiáng)。本研究的實(shí)證分析也表明,醫(yī)療資源的加強(qiáng)對(duì)疫情控制效果明顯,因此對(duì)確診者采取隔離措施也比較重要。

    在圖9和圖10中,g代表醫(yī)療補(bǔ)給程度,1代表醫(yī)療資源不變,5代表醫(yī)療資源補(bǔ)給程度中等,10代表醫(yī)療資源補(bǔ)給程度高。當(dāng)加大醫(yī)療資源投入時(shí),潛伏者數(shù)量的變化趨勢(shì)主要表現(xiàn)為峰值的降低,盡管起到了一定的作用,但效果不是很明顯。確診者數(shù)量的變化趨勢(shì)一方面表現(xiàn)為峰值的降低,且效果特別明顯,說明盡管加大醫(yī)療資源投入對(duì)潛伏者數(shù)量控制效果一般,卻可以十分有效地控制確診者數(shù)量暴發(fā)的程度;另一方面表現(xiàn)為疫情后半階段持續(xù)時(shí)間的大幅縮短,有助于更快地控制疫情。可見,采取積極的醫(yī)療措施不僅有助于縮小疫情導(dǎo)致的勞動(dòng)力供給缺口,而且可以縮短經(jīng)濟(jì)受疫情影響的持續(xù)時(shí)間。

    如果綜合運(yùn)用以上疫情防控措施,即在對(duì)易感者、潛伏者、確診者采取嚴(yán)格隔離措施的同時(shí),加大醫(yī)療資源投入,令模型參數(shù)α1=0.1,α2=0.1,g= 10,則疫情防控效果如圖11所示。可以看出,潛伏者數(shù)量自開始就一直下降,在18天時(shí)便降至1萬人以下,確診者數(shù)量在9天時(shí)達(dá)到最高峰65 000多人,45天可降至1萬人。對(duì)比圖3所示的自然傳播情況,理想假設(shè)下的積極治療加隔離防控措施,既可以通過降低感染數(shù)量縮小勞動(dòng)力供給缺口,又可以通過縮短疫情持續(xù)時(shí)間更早解決沖擊經(jīng)濟(jì)的源頭問題。

    (三)三階段傳染病模型構(gòu)建與模擬

    1.三階段模型的數(shù)據(jù)基礎(chǔ)

    模型把人群分為易感者、潛伏者、確診者、康復(fù)者四類,這樣在設(shè)定初始值時(shí)可能會(huì)導(dǎo)致易感者數(shù)量初始值偏大。從此次湖北省新冠肺炎疫情看,最初暴發(fā)區(qū)域僅限于武漢,如果直接以湖北省總?cè)丝跍p去其他三類人群數(shù)量的差值作為易感者數(shù)量,顯然不符合現(xiàn)實(shí)。SEIR模型的相關(guān)應(yīng)用研究在進(jìn)行相關(guān)參數(shù)估計(jì)時(shí),一般將疫情的整個(gè)過程當(dāng)成一個(gè)階段來估計(jì),但考慮到現(xiàn)實(shí)中疫情會(huì)受到媒體報(bào)道以及政府緊急防控政策等的沖擊和影響,參數(shù)估計(jì)會(huì)隨疫情階段的變化而改變,因此本研究在分析疫情趨勢(shì)時(shí),嘗試對(duì)疫情進(jìn)行分階段考慮,把人員流動(dòng)趨于穩(wěn)定的階段作為分析的起始點(diǎn)。這樣既可以獲得更加符合現(xiàn)實(shí)情況的參數(shù)估計(jì),又由于在疫情擴(kuò)張期,隨著時(shí)間的推移,易感者覆蓋范圍越來越廣直至覆蓋全省,可以在一定程度上糾正易感者數(shù)量初始值偏大的問題。

    根據(jù)本次疫情發(fā)展的時(shí)間軸,本研究將之大致劃分為三個(gè)階段。第一階段為自然傳播期,從2019年12月8日第一例患者出現(xiàn)開始至2020年1月19日武漢通報(bào)一夜新增136名患者為止,在這個(gè)階段公眾對(duì)本次疫情幾乎沒有防備,符合自然傳播的特征。2020年1月20日,鐘南山赴武漢考察,經(jīng)媒體報(bào)道,公眾得知新冠肺炎會(huì)人傳人,加之當(dāng)時(shí)正值春節(jié)返鄉(xiāng)時(shí)段,短時(shí)間內(nèi)武漢有大量人口出城,致使易感區(qū)域逐漸擴(kuò)大到湖北全省。從國家衛(wèi)生健康委員會(huì)統(tǒng)計(jì)的密切接觸者數(shù)據(jù)(圖12)可以發(fā)現(xiàn),新增密切接觸者數(shù)量從2020年1月20日開始至1月29日一直上升,之后趨于平緩。據(jù)此,本研究將2020年1月20日至1月29日這個(gè)時(shí)間段定義為疫情發(fā)展的第二階段。在此階段,人員大量流動(dòng),醫(yī)療資源緊缺,加之新建醫(yī)院、新增醫(yī)護(hù)人員都需要一定時(shí)間,相應(yīng)的管理應(yīng)對(duì)措施也不能一步到位,死亡率和康復(fù)率并未開始改善。2020年1月23日,武漢封城,全國各地陸續(xù)采取嚴(yán)格隔離措施,1月30日以后政府防控措施的作用開始顯現(xiàn),本研究將此時(shí)點(diǎn)開始往后的時(shí)間段定義為第三階段。

    2020年2月7日晚湖北省新型冠狀病毒感染肺炎疫情防控指揮部在例行新聞發(fā)布會(huì)上匯報(bào)的情況顯示,武漢市28家定點(diǎn)醫(yī)療機(jī)構(gòu)床位增加到8 895張,火神山醫(yī)院1 000張病床已經(jīng)交付使用,雷神山醫(yī)院1 600張病床計(jì)劃于2月8日交付使用。截至2020年2月6日,國家及各省市共派遣99支醫(yī)療隊(duì)9 277人支援武漢??梢钥闯觯诖穗A段,隨著醫(yī)療資源的不斷補(bǔ)充,死亡率和康復(fù)率逐步改善(具體參見表3),因此對(duì)模型中的死亡率和康復(fù)率也要進(jìn)行相應(yīng)調(diào)整,以使模擬效果更符合現(xiàn)實(shí)情況。

    2.三階段模型構(gòu)建與參數(shù)擬合

    為驗(yàn)證改進(jìn)模型在疫情預(yù)測(cè)模擬方面的有效性,實(shí)證分析起始點(diǎn)的選擇應(yīng)盡量靠前,因?yàn)樵皆邕M(jìn)行比較準(zhǔn)確的預(yù)測(cè),越可及時(shí)為疫情防控提供合理參考依據(jù)。此外,也要考慮疫情數(shù)據(jù)特征的穩(wěn)定性。根據(jù)前面的分析,在疫情發(fā)展的第三階段選擇預(yù)測(cè)的起始點(diǎn)相對(duì)比較合理。還有,我們對(duì)疫情階段的劃分是以重大事件的發(fā)生為依據(jù)的,而從具體事件的發(fā)生到疫情數(shù)據(jù)的體現(xiàn)存在一定的時(shí)滯,這個(gè)特點(diǎn)從表3的疫情數(shù)據(jù)也可看出,死亡率和康復(fù)率在2020年2月7日以前變化幅度相對(duì)較大。綜合以上分析,最終選取2020年2月7日作為模型預(yù)測(cè)模擬的起始點(diǎn)。

    第三階段的傳播機(jī)制仍然如圖4所示,但在此階段,隨著醫(yī)療資源的不斷補(bǔ)充,死亡率和康復(fù)率均呈現(xiàn)出穩(wěn)定的改善趨勢(shì),因此本研究將模型中的康復(fù)率γ和死亡率v設(shè)為時(shí)變量,初始值分別為表3中的0.010和0.004。改進(jìn)的模型方程為:

    假設(shè)N的初始值N0為湖北省人口6 172.9萬人,I的初始值I0為2020年2月7日留存確診者數(shù)量20 662人,R的初始值R0為2020年2月7日康復(fù)者數(shù)量832人。由于第三階段隔離措施已經(jīng)加強(qiáng),因此潛伏者E的數(shù)量不再等于2020年2月7日確診者數(shù)量的14倍,本研究將此階段E的初始值E0估計(jì)為2020年1月30日留存確診者數(shù)量的14倍,即60 676人。易感者數(shù)量S的初始值S0=N0-E0-I0-R0,即61 646 830人。

    改進(jìn)的模型參數(shù)仍然取α1=0.1,α2=0.1,除死亡率和康復(fù)率外,其他參數(shù)仍然參照表2,醫(yī)療資源強(qiáng)度系數(shù)g=1,表示暫時(shí)沒有特效藥,但醫(yī)護(hù)資源穩(wěn)步改善,表現(xiàn)為模型中的康復(fù)率γ和死亡率v變成了時(shí)變量。

    3.實(shí)證結(jié)果與模型有效性評(píng)價(jià)

    按照改進(jìn)的模型和優(yōu)化的參數(shù)設(shè)計(jì),本研究得出的仿真結(jié)果如圖13所示。現(xiàn)實(shí)中此階段的防控措施比較嚴(yán)格,醫(yī)療資源也源源不斷地支援湖北,在這樣的假設(shè)下模擬結(jié)果顯示,潛伏者數(shù)量立即得到控制,一直處于下降狀態(tài)。得益于嚴(yán)格的隔離措施,潛伏者數(shù)量在初期可以保持較快的下降速度,后期隨著基數(shù)的減少,下降幅度趨于平緩,最終完全得到控制;確診者數(shù)量由潛伏者轉(zhuǎn)化而來,繼續(xù)經(jīng)歷了一段上升期,但10天后便達(dá)到最高峰54 000人左右,之后一直下降,這主要得益于隔離措施的嚴(yán)格(使?jié)摲邤?shù)量不斷減少)和醫(yī)療資源的加強(qiáng)(使康復(fù)率不斷提高);康復(fù)者數(shù)量增長較快,但隨著確診者數(shù)量的不斷減少,康復(fù)者數(shù)量的增長最終也趨于平緩??傮w來看,疫情自防控措施實(shí)施起就可以快速往好的方向發(fā)展,并在50天后慢慢進(jìn)入尾聲。留存確診者數(shù)量預(yù)測(cè)數(shù)據(jù)與實(shí)際數(shù)據(jù)對(duì)比結(jié)果參見圖14,說明改進(jìn)的模型具備預(yù)測(cè)疫情的可靠性。

    4.對(duì)封城措施的模擬分析

    一個(gè)很有效的控制疫情的辦法就是盡早隔離。現(xiàn)實(shí)結(jié)果顯示,疫情暴發(fā)初期,疫情發(fā)源地并未得到較好的隔離處置,待到采取暫停公共交通等封城措施時(shí),已有大量潛伏者和確診者流出武漢,這在理論上直接增加了易感者的數(shù)量。如果在疫情暴發(fā)初期就對(duì)武漢采取封城措施,疫情控制效果能否更好?本研究利用改進(jìn)的SEIR模型對(duì)此進(jìn)行模擬仿真。從國家衛(wèi)生健康委員會(huì)官網(wǎng)發(fā)布的通知可以看到,截至2020年1月19日22時(shí),武漢市累計(jì)報(bào)告新冠肺炎病例198例,已康復(fù)出院25例,死亡3例?!段錆h統(tǒng)計(jì)年鑒》公布的2018年末武漢市常住人口為1 108.1萬人。假設(shè)疫情初期即2020年1月19日之前疫情沒有蔓延到武漢市以外,如果當(dāng)時(shí)對(duì)武漢采取封城措施,將疫情控制在武漢市內(nèi)進(jìn)行處理,則模型中總?cè)丝贜的初始值N0為武漢市常住人口1 108.1萬人,I的初始值I0為170人,R的初始值R0為25人。潛伏者E的初始值E0估計(jì)為確診者的14倍,即2 380人。易感者S的初始值S0= N0-E0-I0-R0,即11 078 425人。假設(shè)采取嚴(yán)格的隔離措施,即模型參數(shù)仍取α1=0.1,α2=0.1,由于早期的死亡率和康復(fù)率尚不穩(wěn)定,不具有一般規(guī)律性,因此本次模擬對(duì)死亡率和康復(fù)率的估計(jì)仍然采用前文對(duì)第三階段的時(shí)變?cè)O(shè)計(jì)。其他參數(shù)仍然參照表2,醫(yī)療資源強(qiáng)度系數(shù)g=1,表示暫時(shí)沒有特效藥,但醫(yī)護(hù)資源情況穩(wěn)步改善。模擬得出的疫情發(fā)展情況參見圖15。

    可以發(fā)現(xiàn),隔離以及醫(yī)療資源等防控措施影響著疫情的走勢(shì),但前期對(duì)疫情暴發(fā)地采取的封城措施決定了易感者的基數(shù)。在模擬封城10~11天后,確診者數(shù)量達(dá)到峰值1 600多人,50天后,確診者數(shù)量僅剩48人,可以說疫情完全得到了控制。但是,如果沒有及時(shí)封城,導(dǎo)致整個(gè)湖北省都成為易感區(qū),疫情完全得到控制要經(jīng)歷更長的時(shí)間,且影響范圍更廣。因此,及時(shí)采取封城措施一方面可以縮小疫情對(duì)經(jīng)濟(jì)的影響范圍,另一方面可以更快地解決源頭問題,推動(dòng)經(jīng)濟(jì)更早地進(jìn)入恢復(fù)期。

    5.對(duì)傳染病存在復(fù)發(fā)情況的模擬分析

    盡管我國在2003年應(yīng)對(duì)非典型肺炎疫情時(shí)已經(jīng)積累了一些經(jīng)驗(yàn),如要對(duì)患者進(jìn)行隔離,但在應(yīng)對(duì)本次傳染病疫情時(shí),卻因多了潛伏者具備傳染性這個(gè)特征而沒有及時(shí)重視對(duì)潛伏者的隔離。因此,在面對(duì)類似重大傳染病疫情時(shí),需要在總結(jié)以往防控經(jīng)驗(yàn)的基礎(chǔ)上,多考慮一些有可能出現(xiàn)的傳染病特征。比如,近期有證據(jù)表明,部分患者在康復(fù)后仍有可能核酸檢測(cè)呈陽性,因此我們有必要考慮傳染病復(fù)發(fā)的可能,并通過實(shí)證分析得出傳染病存在復(fù)發(fā)可能性時(shí)的疫情趨勢(shì)。

    我們?cè)谛鹿诜窝讉魅静∧P偷幕A(chǔ)上添加這一假設(shè),即康復(fù)者存在一定比例的復(fù)發(fā)可能性,并再次具有傳染性,再次成為潛伏者。假設(shè)這一比例為k,可以得出如圖16所示的存在復(fù)發(fā)可能性時(shí)的傳播機(jī)制。

    考慮到疫情逐漸得到控制,醫(yī)護(hù)資源不斷改善,我們假設(shè)k隨時(shí)間遞減,設(shè)k的初始值為0.01,即初始時(shí)每天有1%的康復(fù)者會(huì)復(fù)發(fā)變?yōu)闈摲撸@個(gè)比例隨著時(shí)間的推移到100天時(shí)線性遞減為0。在圖11所示的積極治療防控情況下,對(duì)存在復(fù)發(fā)和不存在復(fù)發(fā)兩種情況進(jìn)行仿真對(duì)比,結(jié)果參見圖17。

    可以看出,如果存在一定比例的康復(fù)者復(fù)發(fā),即使比例不大,也會(huì)延長疫情持續(xù)的時(shí)間,進(jìn)而阻礙社會(huì)經(jīng)濟(jì)的盡快恢復(fù)。因此,傳染病疫情防控工作不能忽視對(duì)康復(fù)者的后續(xù)關(guān)注。一方面,康復(fù)者短期內(nèi)要繼續(xù)做好自我隔離,盡量減少可能的傳染;另一方面,醫(yī)療團(tuán)隊(duì)要對(duì)傳染病是否存在復(fù)發(fā)可能性進(jìn)行必要的研究。

    四、結(jié)論與建議

    (一)結(jié)論

    本研究以SEIR模型為基準(zhǔn),對(duì)處于傳染病流行范圍的人進(jìn)行劃分,分為易感者、潛伏者、確診者、康復(fù)者四類;對(duì)傳統(tǒng)經(jīng)典SEIR模型存在的不足(即估計(jì)參數(shù)時(shí)多以整個(gè)過程估計(jì)出一個(gè)值,不符合現(xiàn)實(shí)中不同階段的情況)進(jìn)行改進(jìn),將重大傳染病事件分為三個(gè)階段,發(fā)展三階段模型;對(duì)相關(guān)參數(shù)加以改進(jìn),把死亡率和康復(fù)率作為時(shí)變量,更符合現(xiàn)實(shí)情況。通過對(duì)經(jīng)典SEIR模型及改進(jìn)的三階段模型分別進(jìn)行參數(shù)估計(jì)和系統(tǒng)仿真,可得出以下三個(gè)主要結(jié)論:

    1.基于對(duì)經(jīng)典SEIR模型的應(yīng)用,在無人為干預(yù)的自然傳播情況下,傳染病感染數(shù)量峰值會(huì)出現(xiàn)在1~2個(gè)月后

    基于湖北疫情數(shù)據(jù)進(jìn)行參數(shù)擬合與動(dòng)態(tài)仿真,結(jié)果發(fā)現(xiàn),在無人為干預(yù)的自然傳播情況下,潛伏者數(shù)量在40天左右時(shí)達(dá)到峰值,確診者數(shù)量在60天左右時(shí)達(dá)到峰值,且要持續(xù)很長時(shí)間才能降到較低水平。

    2.加強(qiáng)隔離、加大醫(yī)療資源投入能顯著減小傳播規(guī)模

    基于經(jīng)典SEIR模型,通過調(diào)整隔離程度和醫(yī)療資源強(qiáng)度參數(shù),模擬不同隔離程度與醫(yī)療資源投入下疫情的發(fā)展趨勢(shì)。結(jié)果發(fā)現(xiàn),隨著對(duì)易感者和潛伏者隔離程度的提高,傳染病整體的傳播規(guī)模顯著減小;隨著醫(yī)療資源投入的加大,如增加醫(yī)護(hù)人員、增加床位、研發(fā)出特效藥等,傳染病的傳播規(guī)模顯著減小。通過對(duì)湖北疫情進(jìn)行系統(tǒng)擬合發(fā)現(xiàn),隔離措施嚴(yán)格時(shí),潛伏者數(shù)量自隔離措施實(shí)施起開始下降,到51天時(shí)便降至1萬人;隔離措施更嚴(yán)格時(shí),只要18天,潛伏者數(shù)量便可降至1萬人。同樣,加大醫(yī)療資源投入后,確診者數(shù)量在9天時(shí)達(dá)到峰值65 000多人,45天時(shí)可降至1萬人。

    3.改進(jìn)的三階段模型能對(duì)傳染病類事件起到很好的預(yù)測(cè)作用,可為未來傳染病評(píng)估和應(yīng)對(duì)提供有效依據(jù)

    改進(jìn)的三階段模型將死亡率和康復(fù)率由靜態(tài)參數(shù)變?yōu)閯?dòng)態(tài)參數(shù)進(jìn)行動(dòng)態(tài)仿真,使得運(yùn)用改進(jìn)后的模型預(yù)測(cè)結(jié)果與實(shí)際情況擬合度較高,結(jié)果誤差很小,說明模型預(yù)測(cè)精度較高,具有可靠性。模型動(dòng)態(tài)仿真結(jié)果表明,加強(qiáng)隔離、加大醫(yī)療資源投入可使疫情快速往好的方向發(fā)展,并于50天后慢慢進(jìn)入尾聲。模擬封城情形發(fā)現(xiàn),模擬封城10~ 11天后確診者數(shù)量達(dá)到峰值,50天后確診者基本消失,說明盡早封城可有效控制疫情規(guī)模,進(jìn)而縮小經(jīng)濟(jì)受沖擊范圍。模擬傳染病存在復(fù)發(fā)可能性的情形發(fā)現(xiàn),若存在一定比例的康復(fù)者復(fù)發(fā),即使比例不大,也會(huì)延長疫情持續(xù)時(shí)間。

    (二)建議

    結(jié)合本研究模型分析結(jié)果,應(yīng)從嚴(yán)格隔離、加強(qiáng)醫(yī)療資源投入、疫情后期持續(xù)關(guān)注三個(gè)角度,對(duì)新冠肺炎疫情中依然存在的問題提出合理建議,完善我國傳染病應(yīng)急防控體系,這有助于今后及時(shí)有效解決類似突發(fā)公共衛(wèi)生事件,從源頭上減少疫情對(duì)經(jīng)濟(jì)的負(fù)面影響。

    1.對(duì)傳染病類突發(fā)公共衛(wèi)生事件盡早隔離

    由模型分析可以看出,阻斷傳播途徑是降低感染數(shù)量的重要手段。要對(duì)突發(fā)公共衛(wèi)生事件做到早發(fā)現(xiàn)早處理,就要對(duì)傳染病信息保持足夠警惕和實(shí)時(shí)關(guān)注,如果不能第一時(shí)間采取有效措施消除隱患,傳染病一旦暴發(fā),波及范圍和損失就會(huì)不斷擴(kuò)大。建議完善傳染病直報(bào)系統(tǒng),擴(kuò)展傳染病直報(bào)系統(tǒng)職責(zé)范圍,方便新發(fā)傳染病信息及時(shí)上報(bào),幫助政府決策。

    重大突發(fā)傳染病疫情中全國各地醫(yī)護(hù)人員普遍缺乏醫(yī)療物資的現(xiàn)象反映了我國醫(yī)用應(yīng)急物資儲(chǔ)備不足的現(xiàn)狀,國家應(yīng)設(shè)立專門機(jī)構(gòu)在國內(nèi)建立應(yīng)急醫(yī)療物資儲(chǔ)備體系,動(dòng)態(tài)管理,定期更換補(bǔ)充應(yīng)急物資,靈活管理應(yīng)急物品儲(chǔ)備,在儲(chǔ)備物品臨近保質(zhì)期或保質(zhì)期限超過一半時(shí)派發(fā)到一線醫(yī)院使用,避免浪費(fèi)。疫情一旦發(fā)生,由主管的地方政府和部門及時(shí)將應(yīng)急物資送往醫(yī)院等抗疫一線機(jī)構(gòu),協(xié)助做好醫(yī)護(hù)人員和患者的隔離,避免醫(yī)護(hù)人員交叉感染。

    2.增加投入提高康復(fù)率是降低傳染數(shù)量的重要手段

    由模型分析結(jié)果可以看出,提高康復(fù)率是降低感染數(shù)量的又一重要手段,因此需要加強(qiáng)對(duì)醫(yī)療資源的集中整合。目前,我國一線二線城市與大部分中小城市、農(nóng)村及偏遠(yuǎn)地區(qū)的公共衛(wèi)生服務(wù)水平差距較大,政府對(duì)農(nóng)村及偏遠(yuǎn)地區(qū)公共衛(wèi)生管理的財(cái)政投入仍然較少,公共衛(wèi)生人才主要集中于大中城市,面對(duì)類似此次新冠肺炎的傳染病等突發(fā)公共衛(wèi)生事件,農(nóng)村及偏遠(yuǎn)地區(qū)的醫(yī)療服務(wù)水平明顯不足,基層疾病防控機(jī)構(gòu)衛(wèi)生應(yīng)急專業(yè)能力相對(duì)不足,高素質(zhì)復(fù)合型人才短缺,處理重大傳染病疫情的專業(yè)隊(duì)伍缺乏,難以有效應(yīng)對(duì)。因此,需要加大對(duì)農(nóng)村及偏遠(yuǎn)地區(qū)醫(yī)療方面的財(cái)政投入,提高農(nóng)村及偏遠(yuǎn)地區(qū)應(yīng)對(duì)傳染病等突發(fā)公共衛(wèi)生事件的能力,滿足其開展公共衛(wèi)生管理的基本需求。

    *感謝北京物資學(xué)院碩士研究生趙輝、趙淑夢(mèng)對(duì)本研究的貢獻(xiàn)。

    參考文獻(xiàn):

    [1]BERNOULLI D.Essai dune nouvelle analyse de la mor? talitécausée par la petite vérole et des avantages de linocu? lation pour la prévenir[D].Paris:AcadéMie Royale Des Sci? ences,1976.

    [2]HAMER W H.The milroy lectures on epidemic diseases in England——the evidence of variability and of persistency of type[J].The lancet,1906,167(4 305):569-574.

    [3]ROSS R.The mathematics of malaria[EB/OL].[2020-04-05]. https://www.bmj.com/content/1/2626/1023.1.

    [4]KERMACK W O,MCKENDRICK A G.A contribution to the mathematical theory of epidemics[J].Proceedings of the royal society of london,1927,115(772):700-721.

    [5]KERMACK W O,MCKENDRICK A G.Contributions to the mathematical theory of epidemics-II. the problem of ende? micity[J].Bulletin of mathematical biology,1991,53(1-2):57-87.

    [6]徐恭賢,馮恩民,王宗濤,等.SARS流行病的SEIR動(dòng)力學(xué)模型及其參數(shù)辨識(shí)[J].黑龍江大學(xué)自然科學(xué)學(xué)報(bào),2005(4):459-462.

    [7]王蓮花,剛毅,張鳳琴.具有常數(shù)輸入的SEIR模型的穩(wěn)定性分析[J].數(shù)學(xué)的實(shí)踐與認(rèn)識(shí),2009,39(18):208-214.

    [8]王樹忠,劉曉宇,李冬梅.一類潛伏期和染病期均有傳染力的SEIR模型的穩(wěn)定性分析[J].哈爾濱理工大學(xué)學(xué)報(bào),2010,15(2):71-75.

    [9]吳昊澄,王臻,何凡,等.基于SEIR模型的流行性腮腺炎暴發(fā)疫情控制效果評(píng)價(jià)[J].中國預(yù)防醫(yī)學(xué)雜志,2013,14(1):72-75.

    [10]崔景安,葉萌,宋國華,等.北京市手足口病的流行趨勢(shì)預(yù)測(cè)[J].生物數(shù)學(xué)學(xué)報(bào),2014,29(1):131-135.

    [11]李冬梅,桂春羽,溫盼盼.一類潛伏期具有常數(shù)輸入率的SEIR模型在流感防控中的應(yīng)用[J].數(shù)學(xué)的實(shí)踐與認(rèn)識(shí),2015,45(12):160-166.

    [12]黃璜,蘭玲,單旭征,等.SEIR模型在諾如病毒醫(yī)院感染的應(yīng)用研究[J].中華醫(yī)院感染學(xué)雜志,2018,28(14):2 205-2 208.

    [13]王昕煒,彭海軍,鐘萬勰.具有潛伏期時(shí)滯的時(shí)變SEIR模型的最優(yōu)疫苗接種策略[J].應(yīng)用數(shù)學(xué)和力學(xué),2019,40(7):701-712.

    Risk Assessment and Response of Public Health Emergencies in Infectious Diseases

    ZHAO Cheng-zhen1,LIANG Xun2and WANG Jun-li1

    (1.School of Economics,Beijing Wuzi University,Beijing 101149,China;2.School of Information,Renmin University of China,Beijing 100872,China;3.Institute of Public Management and Human Resources,Development Research Center of the State Council,Beijing 100010,China)

    Abstract:Major public health events will impact the socio-economic operation,especially when it is contagious,the extent and scope of such impact will be more serious. Establishing the model can help us to simulate public health emergencies in infectious diseases;and it will provide us with the effective basis for the evaluation and countermeasures of that. Based on the SEIR model and considering its deficiency,the authors develop a three-step development model,change related parameters from static ones to dynamic ones,simulate the situation of lockdown and reappearance possibility,and carry out the dynamic simulation based on data of new coronavirus outbreak in 2020. It is found that:under the situation of no human intervention,the peak number of infections will appear in one or two months;strengthening isolation and increasing medical investment will significantly reduce the scale and scope of outbreak;and the improved three-step model can play its role in forecasting infection related issues,and it can provide us with effective basis for the evaluation and countermeasures of the infection in the future. Therefore,to effectively cope with the public health emergencies,we should,first,perfect out prevention system from aspects such as restricting isolation,increasing investment in medical resources and continuously focusing on this issues;second,we should isolating related person as early as possible with the happening of diseases and improve the recovery rate to reduce the number of infections,and reduce to fundamentally weaken the negative impact on economy.

    Key words:infectiousdiseasemodel;publichealthemergency;levelofimpact;extendofimpact;thedurationoftheoutbreak

    猜你喜歡
    突發(fā)公共衛(wèi)生事件持續(xù)時(shí)間
    佳木斯市突發(fā)公共衛(wèi)生事件醫(yī)療救治體系及運(yùn)行機(jī)制研究
    基于突發(fā)公共衛(wèi)生事件中的衛(wèi)生應(yīng)急人力資源管理分析
    突發(fā)公共衛(wèi)生事件中微博、微信的議題呈現(xiàn)異同
    科技傳播(2016年16期)2017-01-03 21:54:06
    學(xué)校突發(fā)公共衛(wèi)生事件防控體系的建立及意義探索
    突發(fā)公共衛(wèi)生事件風(fēng)險(xiǎn)管理的新思考
    The 15—minute reading challenge
    科研院所突發(fā)公共衛(wèi)生事件防控工作的探析
    科技視界(2015年27期)2015-10-08 15:25:42
    基于SVD的電壓跌落持續(xù)時(shí)間檢測(cè)新方法
    極寒與北極氣壓變動(dòng)有關(guān),持續(xù)時(shí)間不確定
    俄語體與持續(xù)時(shí)間結(jié)構(gòu)組合規(guī)律的認(rèn)知語義闡釋
    99久久中文字幕三级久久日本| 亚洲欧美一区二区三区国产| 中国三级夫妇交换| 国产女主播在线喷水免费视频网站| 亚洲欧美一区二区三区黑人 | 一级爰片在线观看| 乱码一卡2卡4卡精品| 97超碰精品成人国产| 免费大片黄手机在线观看| 九九在线视频观看精品| av在线观看视频网站免费| 又大又黄又爽视频免费| 天堂8中文在线网| 国产片特级美女逼逼视频| 亚洲精品乱码久久久v下载方式| 另类亚洲欧美激情| 成人免费观看视频高清| 精品久久久久久久久亚洲| 多毛熟女@视频| 亚洲成人手机| 嘟嘟电影网在线观看| av专区在线播放| 国产白丝娇喘喷水9色精品| 只有这里有精品99| 制服人妻中文乱码| 国产 精品1| av在线app专区| 精品人妻熟女毛片av久久网站| 伊人亚洲综合成人网| 亚洲怡红院男人天堂| 免费大片18禁| 丰满饥渴人妻一区二区三| av专区在线播放| 国产有黄有色有爽视频| 一本色道久久久久久精品综合| 亚洲欧美日韩卡通动漫| 国产无遮挡羞羞视频在线观看| 秋霞伦理黄片| 久久热精品热| 国产精品三级大全| 日韩 亚洲 欧美在线| av在线老鸭窝| 欧美三级亚洲精品| 各种免费的搞黄视频| 男男h啪啪无遮挡| 三级国产精品片| 青春草视频在线免费观看| 久久久久久久久久人人人人人人| 国产一区二区在线观看av| 天美传媒精品一区二区| 免费看av在线观看网站| 久久亚洲国产成人精品v| 欧美精品一区二区大全| 亚洲精品乱久久久久久| 18禁在线播放成人免费| 国国产精品蜜臀av免费| av.在线天堂| 亚洲少妇的诱惑av| 美女视频免费永久观看网站| 国产老妇伦熟女老妇高清| 热re99久久国产66热| 丰满迷人的少妇在线观看| 日韩大片免费观看网站| 国产国语露脸激情在线看| 婷婷色麻豆天堂久久| 永久网站在线| 亚洲丝袜综合中文字幕| 亚洲精品一区蜜桃| 亚洲av不卡在线观看| 在线免费观看不下载黄p国产| 亚洲欧美清纯卡通| a 毛片基地| 18禁观看日本| 国产又色又爽无遮挡免| 在线天堂最新版资源| 中文欧美无线码| 黑人高潮一二区| 国产黄片视频在线免费观看| a 毛片基地| 久久精品国产鲁丝片午夜精品| 91精品一卡2卡3卡4卡| 人体艺术视频欧美日本| 街头女战士在线观看网站| 国产精品人妻久久久影院| 午夜免费鲁丝| 欧美日韩av久久| 国产乱人偷精品视频| 国产伦理片在线播放av一区| 日本黄大片高清| 看非洲黑人一级黄片| 国产国拍精品亚洲av在线观看| 日日摸夜夜添夜夜爱| 水蜜桃什么品种好| 在线观看人妻少妇| 国产精品一区二区三区四区免费观看| 日韩av不卡免费在线播放| 美女国产视频在线观看| 日本爱情动作片www.在线观看| 午夜激情av网站| 十八禁高潮呻吟视频| 欧美激情国产日韩精品一区| 国产免费一级a男人的天堂| 日韩伦理黄色片| 狂野欧美激情性xxxx在线观看| av电影中文网址| 亚洲av不卡在线观看| 亚洲精华国产精华液的使用体验| 日本黄大片高清| 亚洲伊人久久精品综合| 人人妻人人添人人爽欧美一区卜| 精品一区二区三区视频在线| 国产乱人偷精品视频| 91精品三级在线观看| 男的添女的下面高潮视频| 精品一区二区免费观看| 精品少妇久久久久久888优播| 免费日韩欧美在线观看| 中文字幕制服av| 国产精品不卡视频一区二区| 三级国产精品欧美在线观看| 国产成人91sexporn| 啦啦啦中文免费视频观看日本| 91在线精品国自产拍蜜月| 在线观看美女被高潮喷水网站| 亚洲av电影在线观看一区二区三区| 久久久精品94久久精品| 欧美少妇被猛烈插入视频| 免费久久久久久久精品成人欧美视频 | 午夜免费鲁丝| 91精品伊人久久大香线蕉| 丝袜美足系列| 高清欧美精品videossex| 亚洲四区av| 五月天丁香电影| 水蜜桃什么品种好| 丝袜在线中文字幕| 一个人看视频在线观看www免费| av播播在线观看一区| 日韩成人av中文字幕在线观看| 久热久热在线精品观看| 九九在线视频观看精品| 男女高潮啪啪啪动态图| 亚洲人成网站在线播| 婷婷色av中文字幕| 午夜福利视频精品| 国产精品99久久久久久久久| 天堂俺去俺来也www色官网| 男女高潮啪啪啪动态图| 日本色播在线视频| 在线观看免费日韩欧美大片 | 亚洲成色77777| 狠狠精品人妻久久久久久综合| 国产亚洲最大av| 国产精品99久久久久久久久| 免费少妇av软件| 国产男女内射视频| 午夜福利在线观看免费完整高清在| 国产精品一区二区在线观看99| 久久鲁丝午夜福利片| 国产不卡av网站在线观看| 少妇的逼水好多| 国产午夜精品一二区理论片| h视频一区二区三区| 久热这里只有精品99| 少妇的逼水好多| 久久毛片免费看一区二区三区| 亚洲图色成人| 简卡轻食公司| 韩国高清视频一区二区三区| 久久久久久久久久人人人人人人| 97在线视频观看| 国产69精品久久久久777片| 黄色一级大片看看| 少妇的逼好多水| 韩国高清视频一区二区三区| 美女脱内裤让男人舔精品视频| xxx大片免费视频| 熟女人妻精品中文字幕| 国产成人精品福利久久| 一级a做视频免费观看| 国产亚洲一区二区精品| 国产av码专区亚洲av| 国产国拍精品亚洲av在线观看| 97精品久久久久久久久久精品| 伦精品一区二区三区| 国产国语露脸激情在线看| 一本大道久久a久久精品| 午夜视频国产福利| 久久久久精品久久久久真实原创| 狂野欧美白嫩少妇大欣赏| 亚洲,欧美,日韩| 亚洲国产av影院在线观看| 大香蕉97超碰在线| 99国产精品免费福利视频| 久久久久精品久久久久真实原创| 最近的中文字幕免费完整| 欧美激情国产日韩精品一区| 大码成人一级视频| 嘟嘟电影网在线观看| 99国产综合亚洲精品| 免费久久久久久久精品成人欧美视频 | 校园人妻丝袜中文字幕| 一级爰片在线观看| 欧美bdsm另类| 少妇 在线观看| 久久鲁丝午夜福利片| 日韩一本色道免费dvd| 午夜福利在线观看免费完整高清在| 精品一品国产午夜福利视频| av女优亚洲男人天堂| 精品久久国产蜜桃| 国产熟女欧美一区二区| 在线观看www视频免费| 91久久精品电影网| 狠狠婷婷综合久久久久久88av| 欧美最新免费一区二区三区| 黑丝袜美女国产一区| 欧美三级亚洲精品| 亚洲av中文av极速乱| 国产精品无大码| 草草在线视频免费看| 国产黄片视频在线免费观看| 国产日韩欧美视频二区| 精品一品国产午夜福利视频| 狠狠婷婷综合久久久久久88av| 国产成人91sexporn| 亚洲综合色惰| 免费av中文字幕在线| 亚洲欧美日韩另类电影网站| 中文字幕人妻丝袜制服| 91精品伊人久久大香线蕉| 人妻 亚洲 视频| 日本wwww免费看| 国产视频内射| 国产精品 国内视频| 丝袜美足系列| 免费看光身美女| 一边亲一边摸免费视频| 伊人久久国产一区二区| 婷婷色综合大香蕉| 秋霞在线观看毛片| 99视频精品全部免费 在线| 中国国产av一级| 在线观看人妻少妇| 草草在线视频免费看| 亚洲欧美日韩另类电影网站| 成人毛片60女人毛片免费| 亚洲,一卡二卡三卡| 午夜免费观看性视频| 久久久久精品久久久久真实原创| 免费看av在线观看网站| 亚洲精品乱码久久久久久按摩| 午夜91福利影院| 国产成人精品一,二区| 久久久久久久久久人人人人人人| 看非洲黑人一级黄片| 中文字幕最新亚洲高清| 日韩欧美精品免费久久| 久热久热在线精品观看| 性高湖久久久久久久久免费观看| 免费人妻精品一区二区三区视频| 五月玫瑰六月丁香| 精品国产一区二区久久| 久久久久人妻精品一区果冻| 国产黄色视频一区二区在线观看| 亚洲精华国产精华液的使用体验| 精品人妻偷拍中文字幕| 国产国拍精品亚洲av在线观看| 国模一区二区三区四区视频| 国产精品麻豆人妻色哟哟久久| 欧美激情极品国产一区二区三区 | 婷婷色综合大香蕉| av在线播放精品| 欧美 日韩 精品 国产| 欧美激情 高清一区二区三区| 欧美精品一区二区免费开放| 国产高清有码在线观看视频| 久热久热在线精品观看| 夫妻性生交免费视频一级片| 亚洲怡红院男人天堂| 亚洲欧洲精品一区二区精品久久久 | 在线天堂最新版资源| 91午夜精品亚洲一区二区三区| 午夜激情久久久久久久| 国内精品宾馆在线| 少妇精品久久久久久久| 99视频精品全部免费 在线| 菩萨蛮人人尽说江南好唐韦庄| 男男h啪啪无遮挡| 乱码一卡2卡4卡精品| 肉色欧美久久久久久久蜜桃| 久久精品国产自在天天线| 日日爽夜夜爽网站| 亚洲图色成人| 欧美xxⅹ黑人| 国产亚洲欧美精品永久| 最近手机中文字幕大全| 成人综合一区亚洲| 22中文网久久字幕| 亚洲av电影在线观看一区二区三区| 另类精品久久| 精品人妻熟女毛片av久久网站| 少妇被粗大猛烈的视频| 久久人人爽av亚洲精品天堂| 国产亚洲精品第一综合不卡 | 久久99热这里只频精品6学生| 秋霞伦理黄片| 免费久久久久久久精品成人欧美视频 | a 毛片基地| 在线天堂最新版资源| 毛片一级片免费看久久久久| 九草在线视频观看| 国产av国产精品国产| 久热这里只有精品99| 亚洲精品色激情综合| 老熟女久久久| 少妇人妻久久综合中文| 久久精品国产a三级三级三级| 国产精品国产三级专区第一集| 亚洲国产精品成人久久小说| 国产免费福利视频在线观看| 在现免费观看毛片| 精品熟女少妇av免费看| 国产精品久久久久成人av| 国产精品一区二区在线观看99| 精品少妇久久久久久888优播| 日韩成人伦理影院| 9色porny在线观看| 人人妻人人爽人人添夜夜欢视频| 在线观看美女被高潮喷水网站| 午夜免费男女啪啪视频观看| av一本久久久久| 男女无遮挡免费网站观看| 七月丁香在线播放| 国产日韩欧美亚洲二区| 亚洲激情五月婷婷啪啪| 国产免费视频播放在线视频| 久热久热在线精品观看| 国产在线免费精品| 精品一区二区免费观看| 亚洲精品国产av成人精品| 老女人水多毛片| 精品人妻熟女毛片av久久网站| 久久久午夜欧美精品| 一级片'在线观看视频| 中文乱码字字幕精品一区二区三区| 国产不卡av网站在线观看| 日本午夜av视频| 午夜福利视频精品| 成人二区视频| 精品国产露脸久久av麻豆| 日韩av免费高清视频| 妹子高潮喷水视频| 成年人免费黄色播放视频| 在线观看三级黄色| 啦啦啦啦在线视频资源| 国产男女超爽视频在线观看| 亚洲综合色网址| 999精品在线视频| 免费av不卡在线播放| xxx大片免费视频| 国产乱来视频区| 国产黄色免费在线视频| 人妻系列 视频| 久久久午夜欧美精品| 亚洲精品成人av观看孕妇| 91午夜精品亚洲一区二区三区| 亚洲欧美一区二区三区国产| 亚洲一区二区三区欧美精品| 麻豆乱淫一区二区| 欧美日韩视频精品一区| 午夜精品国产一区二区电影| 国产一区二区在线观看av| 高清毛片免费看| 国产一区二区在线观看日韩| 亚洲精品乱码久久久v下载方式| 又大又黄又爽视频免费| 涩涩av久久男人的天堂| 精品久久蜜臀av无| 久久这里有精品视频免费| 国国产精品蜜臀av免费| 熟女av电影| 天堂8中文在线网| 在线精品无人区一区二区三| 国产综合精华液| 国产探花极品一区二区| av卡一久久| 亚洲成色77777| 丰满少妇做爰视频| 国产一区二区三区综合在线观看 | 久久精品国产亚洲网站| 另类亚洲欧美激情| 日韩伦理黄色片| 日本欧美国产在线视频| 国产探花极品一区二区| 成人国语在线视频| 在线亚洲精品国产二区图片欧美 | 热99国产精品久久久久久7| 九九爱精品视频在线观看| 五月开心婷婷网| 久久久久视频综合| 亚洲综合精品二区| 午夜精品国产一区二区电影| 久久久久久久久久久免费av| 日本vs欧美在线观看视频| 亚洲人与动物交配视频| 日本黄色日本黄色录像| 国产伦精品一区二区三区视频9| 夜夜骑夜夜射夜夜干| 人妻一区二区av| 亚洲精品色激情综合| 国产精品蜜桃在线观看| 国产熟女欧美一区二区| 男女高潮啪啪啪动态图| 精品亚洲成a人片在线观看| 蜜臀久久99精品久久宅男| 美女cb高潮喷水在线观看| 亚洲人成网站在线观看播放| 亚洲av成人精品一二三区| 午夜激情福利司机影院| 人体艺术视频欧美日本| 精品一区在线观看国产| 在线观看人妻少妇| 国产极品天堂在线| 观看美女的网站| 欧美成人午夜免费资源| 人成视频在线观看免费观看| 亚洲av电影在线观看一区二区三区| av在线app专区| 国产免费福利视频在线观看| 91久久精品国产一区二区成人| 一级毛片aaaaaa免费看小| 91午夜精品亚洲一区二区三区| 久久人人爽av亚洲精品天堂| 男人添女人高潮全过程视频| 亚洲欧洲国产日韩| 狂野欧美激情性bbbbbb| 久久久久久久久久人人人人人人| 免费大片18禁| h视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 久久久精品区二区三区| 九九在线视频观看精品| 一二三四中文在线观看免费高清| 精品一品国产午夜福利视频| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产欧美日韩在线播放| 日韩不卡一区二区三区视频在线| 精品一区二区三区视频在线| 亚洲精品视频女| 国产精品久久久久久精品古装| 99热国产这里只有精品6| av黄色大香蕉| 高清av免费在线| 亚洲精华国产精华液的使用体验| 日日摸夜夜添夜夜爱| 人体艺术视频欧美日本| 亚洲欧美一区二区三区国产| 国产熟女午夜一区二区三区 | 午夜激情久久久久久久| 亚洲国产精品999| 国产片内射在线| 免费人成在线观看视频色| 色吧在线观看| 极品少妇高潮喷水抽搐| 亚洲国产av影院在线观看| 十八禁高潮呻吟视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品国产av在线观看| 啦啦啦中文免费视频观看日本| 成人亚洲欧美一区二区av| 亚洲国产色片| h视频一区二区三区| 免费大片18禁| 欧美日韩综合久久久久久| 下体分泌物呈黄色| 91aial.com中文字幕在线观看| 免费黄网站久久成人精品| 久久久久视频综合| 亚洲国产精品999| 丝袜美足系列| 交换朋友夫妻互换小说| 美女视频免费永久观看网站| 中国三级夫妇交换| 欧美老熟妇乱子伦牲交| 国产精品久久久久久av不卡| 91久久精品电影网| 亚洲国产欧美日韩在线播放| 爱豆传媒免费全集在线观看| 91精品一卡2卡3卡4卡| 在线观看免费高清a一片| 国产精品久久久久久久电影| 女性被躁到高潮视频| 制服丝袜香蕉在线| 午夜91福利影院| 97超碰精品成人国产| 2018国产大陆天天弄谢| 国产在线一区二区三区精| 2022亚洲国产成人精品| 热99国产精品久久久久久7| 国产成人精品在线电影| 中国国产av一级| 午夜影院在线不卡| 精品人妻熟女毛片av久久网站| 久久久精品94久久精品| 青春草视频在线免费观看| 在线观看三级黄色| 国产免费又黄又爽又色| 精品久久蜜臀av无| 黄色欧美视频在线观看| 高清av免费在线| 美女中出高潮动态图| 欧美精品国产亚洲| 精品一区二区三卡| 欧美精品国产亚洲| 中文字幕久久专区| 国产国拍精品亚洲av在线观看| 日韩视频在线欧美| 国产欧美另类精品又又久久亚洲欧美| 91成人精品电影| 国产欧美另类精品又又久久亚洲欧美| 少妇精品久久久久久久| 黑人高潮一二区| 一级爰片在线观看| 一边亲一边摸免费视频| 在线观看免费视频网站a站| 精品久久国产蜜桃| 午夜免费男女啪啪视频观看| 18禁在线播放成人免费| 十八禁网站网址无遮挡| 91精品国产国语对白视频| 美女脱内裤让男人舔精品视频| 国产老妇伦熟女老妇高清| 18在线观看网站| 国产精品久久久久久精品古装| 99热国产这里只有精品6| 成人综合一区亚洲| 国产av国产精品国产| 我的老师免费观看完整版| 欧美 日韩 精品 国产| 又粗又硬又长又爽又黄的视频| 麻豆成人av视频| 午夜福利影视在线免费观看| 国产黄频视频在线观看| 久久久国产欧美日韩av| 人人澡人人妻人| 肉色欧美久久久久久久蜜桃| 另类精品久久| 国产成人精品婷婷| 国产在线免费精品| 狂野欧美白嫩少妇大欣赏| 日韩制服骚丝袜av| a级片在线免费高清观看视频| 亚洲,一卡二卡三卡| 成人毛片a级毛片在线播放| 蜜桃久久精品国产亚洲av| 国国产精品蜜臀av免费| 成人亚洲欧美一区二区av| a级毛色黄片| 各种免费的搞黄视频| 国产探花极品一区二区| 美女国产高潮福利片在线看| 日韩不卡一区二区三区视频在线| 少妇被粗大的猛进出69影院 | 精品一区二区三区视频在线| 人人妻人人澡人人看| 国产又色又爽无遮挡免| 能在线免费看毛片的网站| 国产精品成人在线| 午夜影院在线不卡| 久久 成人 亚洲| 伊人亚洲综合成人网| 99久久人妻综合| 在线观看国产h片| 夜夜看夜夜爽夜夜摸| 999精品在线视频| 国产精品免费大片| 成人亚洲精品一区在线观看| 日韩一区二区三区影片| 日韩中文字幕视频在线看片| 日本黄色日本黄色录像| 精品人妻熟女av久视频| 精品一区二区免费观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲欧洲国产日韩| 在线播放无遮挡| 国产亚洲精品久久久com| 五月伊人婷婷丁香| 少妇的逼水好多| 最近手机中文字幕大全| 欧美日韩精品成人综合77777| 在线观看免费视频网站a站| 18禁在线无遮挡免费观看视频| av播播在线观看一区| 亚洲三级黄色毛片| 好男人视频免费观看在线| 91国产中文字幕| 久久精品熟女亚洲av麻豆精品| 嫩草影院入口| 国产一区二区三区av在线| 一级毛片aaaaaa免费看小| av国产久精品久网站免费入址| 高清黄色对白视频在线免费看| 亚洲精品成人av观看孕妇| 成人手机av| 狠狠精品人妻久久久久久综合| 亚洲欧洲国产日韩| 丝袜喷水一区| 中文精品一卡2卡3卡4更新| 亚洲精品久久成人aⅴ小说 | 亚洲国产精品成人久久小说| 国产极品天堂在线| 中文字幕免费在线视频6| 各种免费的搞黄视频| 亚洲精品一二三|