2019年12月31日，武汉市卫生健康委员会首次通报了27例“不明原因感染肺炎”^[1]。2020年1月7日确定导致该肺炎的病原体为新型冠状病毒^[2]，WHO在2月11日将疾病正式命名为新型冠状病毒肺炎（COVID-19）^[3]。COVID-19在人群中具有较强的传染性，可能高于严重急性呼吸综合征（severe acute respiratory syndrome，SARS）和中东呼吸综合征（middle east respiratory syndrome，MERS）^[4-6]，我国将其纳入法定乙类传染病管理，并采取甲类传染病的预防控制措施^[7]。截至2020年2月25日24：00，全国34个省（自治区、直辖市、特别行政区）累计报告COVID-19确诊病例78 190例，其中广东省1 347例，是湖北省外累计病例数最多的省份^[8-9]。已有研究显示，各省份疫情的严重程度受到湖北省人口流入的影响，广东省疫情输入风险较高^[10]。尤其是广东省当前处于全面复工复产阶段，人口输入规模较大，输入疫情的风险可能增高。本研究结合21个城市的人口迁入和输入疫情的省份报告病例数情况，对广东省各地级市疫情的输入风险进行评估和预警，为各市疫情研判和防控提供支撑。

对象与方法

1.数据来源：从国家卫生健康委员会及各省级卫生健康委员会官方网站收集各省份和广东省每个地级市每日报告新增病例数（不含临床确诊病例），时间截至2020年2月25日。我国30个省份（不含中国香港地区、中国澳门地区和中国台湾地区数据）迁入至广东省21个地级市的迁入指数来源于百度迁徙（http://qianxi.baidu.com），表征人口迁入规模，时间为2020年1月10日至2月25日。各省份的人口数据来源于2019年中国统计年鉴。根据2019年广东省统计年鉴^[11]，可将广东省21个地级市划分为4个区域：珠三角地区包括广州、深圳、珠海、佛山、江门、东莞、中山、惠州和肇庆；粤东地区包括汕头、汕尾、潮州和揭阳；粤西地区包括湛江、茂名和阳江；粤北地区包括韶关、河源、梅州、清远和云浮。

2.疫情输入风险指数计算：基于前期的研究^[10]，假设迁入的人群符合来源地的人群分布，且发病风险与来源地发病率相一致，则人口迁入导致的疫情输入风险可由迁入规模、来源地发病人数和人口进行表示。由于2020年后全国首次报告COVID-19病例的时间为1月10日，故计算1月10日至2月25日的疫情输入风险指数，计算公式：

式中，Risk_in，t表示第t天的疫情输入风险指数，case_i，t表示第t天i省的新增病例数，pop_i表示i省的人口数，iMI_i，t表示第t天由i省迁入的人口规模指数。进行风险来源分析时，对疫情输入风险指数进行自然对数标准化，公式如下：

式中，SRisk_in表示标化疫情输入风险指数，Risk_in表示疫情输入风险指数。

根据滞后时间计算每个城市疫情输入风险指数的滑动累加值，综合21个地级市所有滑动累加风险指数，取四分位数作为阈值，划分风险预警等级，由小到大分别为“低风险”“中风险”“中高风险”和“高风险”。

3.统计学分析：采用Excel 2019软件对数据进行初步整理，分析广东省累计报告病例和每日新增病例的时间变化趋势；采用R 3.6.1软件描述病例空间分布特征，计算和分析每个地级市的输入风险指数、分析输入风险不同滞后天数与新增病例数的关系，并计算每个城市未来的疫情输入风险。采用双侧检验，以P＜0.05为差异有统计学意义。

结果

1.疫情基本情况：截至2020年2月25日24：00，广东省累计报告确诊病例1 347例，主要聚集在珠三角地区，占90.0%（1 212/1 347）；除云浮市外，其他市均报告病例，其中深圳市报告病例数最多（417例），其次为广州市（346例）。

2.疫情输入风险：2020年1月10日至2月25日，广东省平均每天的疫情输入风险指数为44.03，其中深圳、东莞、广州的风险排前3位，分别为11.72、10.25和7.25（表 1）。将每个市来源于不同省份的风险指数进行对数标准化，见图 1。除了湛江的疫情输入最高风险来自海南省以外，其他市均来自湖北省。此外，位于珠三角地区和粤北地区的城市来自湖南省、江西省和广西壮族自治区输入疫情风险也较高；粤东地区受到福建省疫情输入的影响较大；粤西地区的风险主要来自海南省。

3.疫情输入风险指数与报告病例数相关性：分别对广东省滞后0~9 d的疫情输入风险指数和每日新增病例数进行相关性分析，所有相关系数均具有统计学意义，且滞后4 d时，二者的相关系数最高（r＝0.73）。因此，下面的风险预警用滞后4 d的风险指数。见图 2。

4.疫情输入风险预警：各城市未来4 d的累计风险预警等级时空分布见图 3，东莞、深圳、广州和中山近期持续具有高输入风险，其他城市风险出现降低趋势。以2月25日后的4 d为例，东莞、深圳、中山、广州、佛山和惠州风险等级为高，4 d的累计疫情输入风险指数分别为38.85、21.59、11.67、11.25、6.19和5.92，且单日风险整体呈升高趋势。低风险城市有6个，分别为云浮、汕尾、阳江、梅州、茂名和潮州。每个城市输入风险来源排名前5的省份中，湖北仍是最高，其次为湖南、广西和重庆。见图 4。

讨论

COVID-19作为一种新发传染病，具有较强的人际传播能力^[12]。控制传染源和阻断传播途径是最主要的防控策略。已有研究显示，人口流动对疫情传播风险具有重要的影响^[10]。广东省作为外来人口大省，面临输入疫情的风险更高。因此，广东省制定了“外防输入，内防扩散”的防控策略，并采取了一系列严格的防控措施，目前每天报告病例数在个位数上波动，已取得了阶段性的成效。本研究基于人口迁徙数据，对广东省21个地级市“外防输入”进行风险评估，筛选重要输入风险地，为各市风险预警和精准防控提供了技术支撑。

截至2020年2月25日，广东省累计报告1 347例COVID-19病例，占湖北省以外33个省（直辖市、自治区、特别行政区）病例总数的10.36%（1 347/13 003）^{[8, 13]}。90%的病例主要集中在珠三角地区，这可能与珠三角地区外来人口较多有关。人口数据显示，2018年年末珠三角地区的常住人口为6 300.99万，户籍人口为3 628.04万^[11]。有文献报道，截至2月7日广东省近80%的病例均为输入性病例^[6]。疫情输入高风险城市主要分布在珠三角地区，对每个城市的疫情输入风险来源进行分析，发现湛江外其他城市的疫情输入最高风险均来自湖北省。这可能与湖北省的高发病率和人口流入多均有关。据百度迁徙大数据显示，在2020年1月10日至2月25日广东省的总人口迁入指数中，湖北省占2.86%，排第8位；其中1月24日前累计迁入占比5.08%，排第5位^[14]，提示早期湖北省输入广东省各市的疫情比较严重。此外，广东省的邻近省份如湖南、江西等省也是各城市输入疫情的来源。

本研究分析了广东省不同滞后时间的疫情输入风险指数与每日新增病例数的相关性，发现滞后时间为4 d的相关系数最高。这与Guan等^[15]基于全国1 099例确诊病例报告的潜伏期中位数相一致，提示风险指数可能对未来4 d的疫情发展具有较好的预警作用。基于此，进一步计算了广东省21个城市的滑动累计风险指数，并综合所有城市的4 d滑动累计风险指数划分风险等级。在未来4 d的风险预警中，东莞、深圳、广州和中山市具有高输入风险；进一步对风险来源进行分析，发现尽管各市有所不同，但主要来源仍是湖北省和广东省的邻近省份。提示在未来短期时间内，外防输入的策略仍需有力落实，尤其是外来人口众多的城市（如东莞、深圳、广州等）更需继续保持疫情控制力度，对工厂、企业、工地、出租屋等复工复产的重点场所要采取精细化防控措施，以及安排就业人员错峰出行避免在公共交通工具发生大规模人群聚集，对来自高风险省份的迁入人员加强筛查和健康管理，最大程度降低疫情输入风险。

本研究存在局限性。首先，疫情输入风险指数基于人口迁入来源地的每日报告病例数计算，未考虑无症状感染者，可能低估了输入疫情的风险；第二，输入风险指数的只考虑输入性疫情，对由输入病例导致的本地暴发预警能力有限；第三，流入广东的人群可能并不是来源地人群的一个无偏样本，一定程度上影响风险评估和预警的准确性。

利益冲突  所有作者均声明不存在利益冲突

志谢 感谢各级卫生健康委员会公开COVID-19疫情数据；感谢百度地图团队公开发布百度迁徙指数；向奋战在新型冠状病毒肺炎疫情防治前线的广大工作人员致敬！