吳家園，賴天文，潘振宇，李筱，劉華鋒，呂軍

2019年12月，武漢市出現(xiàn)了一系列非典型病毒性肺炎病例。2020年1月7日，我國確認(rèn)這類疾病的病原為一種新型冠狀病毒，1月12日世界衛(wèi)生組織將其命名為2019-nCoV，該病原感染所致的肺炎稱為新型冠狀病毒肺炎（novel coronavirus pneumonia，NCP）。被發(fā)現(xiàn)之后，該疾病在全國各地乃至周邊國家迅速擴散，各地疑似/確診病例數(shù)在短時間內(nèi)急劇上升。截至2月17日24時，我國已累計報告確診NCP病例72 436例，累計死亡病例1868例，其中，湖北累計確診59 989例（武漢42 752例），累計死亡1789例（武漢1381例）[1]。

由于疫情暴發(fā)突然，直至2020年1月23日后，全國和各地的衛(wèi)生健康委員會公布的疫情數(shù)據(jù)才日趨準(zhǔn)確[2]。由于急性傳染病的流行具有擴散快，地方性和季節(jié)性等特點，新興的大數(shù)據(jù)分析手段并不適用，需要利用傳統(tǒng)的傳染病動力學(xué)模型來分析和預(yù)測疫情的發(fā)展。在傳染病動力學(xué)模型中，基本再生數(shù)（R0）是評估傳染病傳播能力的關(guān)鍵指標(biāo)，其定義是指在一個由易感者構(gòu)成的群體中，單個已被感染的個體在喪失傳播病毒能力之前平均能感染的人數(shù)[3]。R0反映了傳染性疾病暴發(fā)的潛力和嚴(yán)重程度，數(shù)值越高表示疫情越嚴(yán)重。由于NCP疫情暴發(fā)于春節(jié)前夕，大規(guī)模的春運人口流動和人群聚集加快了疫情發(fā)展，給我國疫情的防控帶來了嚴(yán)峻的考驗。疫情暴發(fā)后，各地政府相繼啟動了重大突發(fā)公共衛(wèi)生事件Ⅰ級響應(yīng)，并采取了一系列嚴(yán)厲的防控措施，包括全面排查疑似/確診患者并隔離觀察或治療，嚴(yán)格公共場所開放管理，出行限制，禁止聚會，加強健康管理等，在一定程度上減少了人群的流動性和接觸率，有效地控制了疫情的發(fā)展。但長期的管控也給人民群眾的日常生活和各行各業(yè)的生產(chǎn)帶來了諸多不便，對我國的經(jīng)濟發(fā)展也帶來了很大的影響。目前我國的NCP疫情已經(jīng)到了防控的關(guān)鍵時期，利用合適的模型來估算現(xiàn)階段的R0以及預(yù)測其未來的趨勢，有助于及時準(zhǔn)確地評估干預(yù)措施的有效性，預(yù)測措施需要維持的時間，為部署下一階段的疫情防控提供有力的證據(jù)以及減少群眾的恐慌情緒。

本研究通過利用國家衛(wèi)生健康委員會公布的NCP疫情的官方數(shù)據(jù)，構(gòu)建NCP傳播的SIR模型，并對本次疫情的R0進(jìn)行估計和預(yù)測，旨在評估疫情的現(xiàn)狀，從而預(yù)測疫情的發(fā)展趨勢，為疫情防控提供依據(jù)與支持。

1 方法

1.1 數(shù)據(jù)來源本文數(shù)據(jù)來源于國家衛(wèi)健委截至2020年2月18日24時的官方統(tǒng)計數(shù)據(jù)[1]。除了全國的統(tǒng)計數(shù)據(jù)，我們還提取本次疫情發(fā)源地和最嚴(yán)重的湖北省，以及部分具有代表性（如累計確診人數(shù)較多，人口密度較大，或外來務(wù)工較多）的省市（廣東、河南、浙江和上海）的統(tǒng)計數(shù)據(jù)進(jìn)行分析，以全面評估和預(yù)測我國NCP疫情發(fā)展趨勢。全國及各省市自治區(qū)2018年底的人口數(shù)來源于國家統(tǒng)計局官方數(shù)據(jù)[4]，全國為14億人，湖北省為5917萬，廣東省為11 346萬，河南省為9605萬，浙江省為5737萬，上海市為2424萬。

1.2 SIR模型我們采用SIR模型來描述我國NCP疫情的發(fā)展趨勢。在SIR模型中，疫區(qū)內(nèi)的人群被分為三個部分，分別為易感人群（Susceptible，S）、感染人群（Infected，I）和移除人群（Removed，R）。其中，易感人群指尚未染病但容易被感染的人群，感染人群指已感染并可以傳播給他人的人群，移除人群指不再參與感染與被感染過程的人群（包括病愈且具有免疫能力或已死亡人群）。S(t)，I(t)和R(t)分別表示t時刻易感人群、感染人群和移出人群的數(shù)量，同時假設(shè)該傳染病滿足以下法則：①疫區(qū)內(nèi)的總?cè)丝跀?shù)（N）保持恒定；②單位時間內(nèi)，感染人群所能感染的人數(shù)與易感人群的人數(shù)成正比，感染強度為β；③單位時間內(nèi)，感染人群中移除的人數(shù)與感染人群的數(shù)量成正比，移除概率為γ。據(jù)此我們可得以下經(jīng)典SIR模型的微分方程組[5]：

其中，N=S(t)+I(t)+R(t)，β是感染人群每天平均接觸的易感人數(shù)k和傳染概率b的乘積，γ取決于可傳播的時間長度C，γ=1/C。

1.3 指數(shù)平滑法指數(shù)平滑法是由Robert G. Brown所提出，是估計時間序列趨勢性的一種常用方法，主要是通過對以往實際觀察值的差異性賦權(quán)（序列中近期的數(shù)據(jù)賦予較大的權(quán)重，遠(yuǎn)期的數(shù)據(jù)賦予較小的權(quán)重）來實現(xiàn)歷史數(shù)據(jù)的加權(quán)修勻，從而提取時間序列中的確定性信息[6]。其突出的優(yōu)點在于預(yù)測結(jié)果既反映最新的信息，也不好遺漏歷史信息，并且隨著時間的流逝，過去數(shù)據(jù)的影響逐漸減小，主要用于中短期的預(yù)測[7]。本研究將SIR模型與指數(shù)平滑法相結(jié)合，通過計算基本再生數(shù)R0預(yù)測值與真實值之間的殘差平方之和，即進(jìn)行估計與預(yù)測。本研究本研究所有計算和模型擬合過程均采用R語言（3.6.1版本）的deSolve包和sweep包來完成。

2 結(jié)果

2.1 我國當(dāng)前NCP疫情的分析與預(yù)測表1列出了從2020年1月17日至2月18日我國NCP疫情發(fā)展過程中R0的具體數(shù)據(jù)。其中，R0的最大值為2.94，見于疫情初期（1月19日），隨后，R0呈波動下降的趨勢，至2月18日，達(dá)到本觀察期內(nèi)的最小值1.31（圖1），說明我國目前所采取的疫情防控措施是積極有效的。

表1 2020年1月17日至2月18日我國新型冠狀病毒肺炎疫情的基本再生數(shù)R0

圖1 2020年1月17日至2月18日我國新型冠狀病毒肺炎基本再生數(shù)的變化趨勢

圖2顯示了指數(shù)平滑法對我國NCP基本再生數(shù)R0未來趨勢的預(yù)測結(jié)果。依據(jù)目前的數(shù)據(jù)的預(yù)測結(jié)果顯示，大約在2月底，我國NCP疫情的基本再生數(shù)R0將回落到低于1的水平，預(yù)示著疫情可能迎來“拐點”，形勢正逐漸好轉(zhuǎn)。大約在3月中下旬，基本再生數(shù)R0將回落到低于0的水平，預(yù)示著疫情將走向消亡。

2.2 各省市NCP疫情發(fā)展趨勢的分析圖3展示了各?。ㄊ?、自治區(qū)）NCP疫情基本再生數(shù)R0的變化趨勢。各?。ㄊ?、自治區(qū)）NCP疫情基本再生數(shù)R0均在疫情初期達(dá)到最高值。作為疫情起源地，湖北省疫情初期的基本再生數(shù)R0最高（6.19），急劇下降到2的水平后，保持一個緩慢下降的趨勢，說明盡管湖北省所采取的一系列疫情防控措施取得了一定的效果，但現(xiàn)階段的防疫形勢仍有壓力。而湖北省以外的各?。ㄊ?、自治區(qū)）的基本再生數(shù)R0在疫情初期達(dá)到峰值（2.43）之后，保持一個波動下降的趨勢，并且下降的速度較快。幾個累計確診病例數(shù)最多的省份，如廣東、河南、浙江，以及人口密度較大的地區(qū)，如上海市，也保持一個波動下降的趨勢，說明這些省市在疫情早期所實施的防控措施已經(jīng)起到了效果。

圖2 我國新型冠狀病毒肺炎基本再生數(shù)的變化與預(yù)測

圖3 各?。ㄊ小⒆灾螀^(qū)）新型冠狀病毒肺炎基本再生數(shù)的變化趨勢

3 討論

本研究結(jié)果顯示，全國NCP基本再生數(shù)R0在疫情初期最高，為2.94，與Imai[8]的估算結(jié)果（2.8～3.3）相吻合，略低于美國東北大學(xué)的報告結(jié)果（3.2～3.9）[9]，說明本研究的估算結(jié)果還是較為準(zhǔn)確的。與2003年同樣由冠狀病毒引起的SARS疫情相比，NCP疫情的基本再生數(shù)R0與之持平，說明NCP傳播能力與SARS相當(dāng)，同樣屬于中等傳染能力的傳染性疾病。通常將基本再生數(shù)R0區(qū)間上限作為傳染病最大傳播力的參考參數(shù)，很多已知傳染病的傳播能力都高于NCP，如寨卡病毒（1.4～6.6）[10]，中東呼吸綜合癥（2.0～6.7）[11]，天花（2.5～6.0）[12]，故可認(rèn)為NCP并不具體特別突出的傳染能力。根據(jù)我國對抗SARS的經(jīng)驗，越早采取嚴(yán)格的防控措施越有利于盡早控制疫情。

作為疫情發(fā)源地，湖北省的基本再生數(shù)R0在全國各省（市，自治區(qū)）中最高，并在疫情初期（1月中旬）達(dá)到峰值，并且傳播勢頭顯得比SARS還要兇猛，這是在無任何有效干預(yù)和防控手段的情況下疫情自然傳播的結(jié)果。自1月21日起，政府果斷采取各種防控措施，有效的遏制了NCP疫情的蔓延勢頭，基本再生數(shù)R0快速下降至2的水平，并維持波動下降的趨勢。從預(yù)測的結(jié)果來看，湖北省NCP疫情暴發(fā)在未來的一段時間內(nèi)仍繼續(xù)持續(xù)，其可能的原因在于：①湖北當(dāng)?shù)氐尼t(yī)療資源嚴(yán)重不足，不能及時發(fā)現(xiàn)、收治和隔離已感染的患者，導(dǎo)致未能及時控制感染者人數(shù)（I）；②將CT影像結(jié)果作為湖北省重災(zāi)區(qū)NCP“臨床診斷病例”的診斷依據(jù)，在一定程度上增加了新增患者的例數(shù)；③同樣具有傳播病毒能力的無癥狀感染者未及時就診，增加了傳播時間（C）。目前湖北省，尤其是武漢市，每天新增確診例數(shù)仍高達(dá)1600多例，說明此次疫情的傳染源仍未得到有效的控制。鑒于之前湖北省和武漢市醫(yī)療資源嚴(yán)重不足，許多感染者無法得到及時救治和隔離的情況，調(diào)動更多的資源支援疫情重災(zāi)區(qū)，有利于集中力量撲滅疫情，防止疫情的進(jìn)一步蔓延。

與湖北省相比，其他省份的疫情防控形勢比較樂觀，基本再生數(shù)呈持續(xù)下降的趨勢，尤其是累計確診病例數(shù)最多的幾個省份。經(jīng)歷了SARS疫情之后，各地對付急性傳染病疫情有了更豐富的經(jīng)驗。如廣東省，從發(fā)現(xiàn)首輸入性病例到啟動重大突發(fā)公共衛(wèi)生事件Ⅰ級響應(yīng)，僅僅不到一個月的時間；河南省采取的“硬核”防疫，2019年12月底已實行交通管制，21日起停止銷售活禽，新增130個定點醫(yī)院和525個發(fā)熱門診。本研究的結(jié)果也進(jìn)一步說了目前各省市所采取的防控措施是積極有效的，同時，也得益于民眾防范意識的提升以及十多年來我國醫(yī)療衛(wèi)生事業(yè)的進(jìn)步。

對于傳染病疫情的防控，關(guān)鍵在于快速將疾病再生數(shù)R0控制在1以下[13]。在傳染病動力學(xué)理論里，基本再生數(shù)R0取決于感染人群每天平均接觸的易感人數(shù)（k），傳染概率（b），以及可傳播的時間長度（C），其中k又與感染人群（I）和易感人群（S）的數(shù)量有密切關(guān)聯(lián)[14]?；诹餍胁W(xué)調(diào)查結(jié)果，我們對于新型冠狀病毒的傳播途徑有了更加清晰的認(rèn)識，包括空氣飛沫和氣溶膠傳播，接觸傳播以及可能的糞口傳播。目前政府所采取的各種防控措施，包括：集中收治和隔離感染者，可以有效減少感染人數(shù)（I）；限制出行、取消聚會，減少人群的接觸幾率，可以有效減低感染者接觸易感者的數(shù)量（k）；戴口罩、勤洗手，做好個人清潔衛(wèi)生，可以降低傳染概率（b）；出現(xiàn)疑似癥狀或者與高危人群有密切接觸后，自行隔離或接受醫(yī)學(xué)觀察，可縮短傳播時間（C）。因此，繼續(xù)保持目前的防控手段不松懈，對控制疫情的發(fā)展是非常重要的。

需要警惕的是，由于節(jié)后“返工潮”的來臨，可能對各地NCP疫情的防控帶來巨大的壓力，尤其是勞務(wù)輸入性省市，包括廣東、浙江、上海等。為防止疫情的二次爆發(fā)，積極落實早發(fā)現(xiàn)、早報告、早隔離、早治療是仍控制NCP疫情的首要手段。同時宣傳正確的生活習(xí)慣，注重個人衛(wèi)生，加強自我保護，包括正確佩戴口罩，正確洗手，合理膳食，適量運動，戒煙限酒，心理平衡；在到達(dá)目的地后，自覺在家隔離14 d，都將有助于控制疫情。

本研究的研究數(shù)據(jù)來源于公開數(shù)據(jù)，僅考慮了已確診的和有癥狀的感染者，但在病例報告過程中可能存在漏報或報告不及時等情況，同時未涵蓋無癥狀感染者，所以我們的研究存在一定的局限性和信息偏倚，可能會導(dǎo)致基本再生數(shù)R0的估算結(jié)果不能真實反映疫情的實際情況。同時，本研究的假設(shè)相對較為理想，這可能會導(dǎo)致對預(yù)測結(jié)果造成一定的影響。

綜上所述，本文通過估計NCP疫情基本再生數(shù)R0，預(yù)測疫情未來的發(fā)展趨勢。結(jié)果表明，在現(xiàn)有的高效防控措施下，我國NCP疫情正逐步得到控制，但湖北省的防控形勢仍十分嚴(yán)峻。繼續(xù)保持目前的防控手段不松懈，對控制疫情的發(fā)展非常關(guān)鍵。同時，調(diào)動更多的資源支援疫情重災(zāi)區(qū)湖北省，有利于集中力量撲滅疫情，防止疫情的進(jìn)一步蔓延。