自年3月中旬以来,上海每日新冠确诊病例和无症状感染者的数量一路飙升。4月5日当日,新增本土新冠肺炎确诊病例例和无症状感染者例(图2)。目前疫情规模已经超过了武汉,并且还没出现拐点,上海的防疫到了最艰难的时刻。
图1:疫情笼罩下的上海上海曾是全国抗疫肺炎的模范,为何如今却成为疫情的重灾区?专家认为原因分客观和主观两方面。客观方面,奥密克戎变异株的传染性和隐蔽性较强,该变异株的传播能力比奥密克戎原始株高30%左右,而且所造成的感染大部分是轻症或无症状感染,能够逃过常规核酸检测,所以隐蔽性更强,这也给精准防控带来了困难。此外,上海是人口众多的国际大都市,国内外的贸易和人员往来太密切而复杂。主观方面,防疫人员防控意识松懈,对防疫过于自信,没能早发现早制止。当意识到疫情的严重性时,病毒早已全面蔓延开来,甚至连累了全国其它很多城市。此外,还有很多群众不配合防控也给防疫带来很大的挑战。例如有人初筛阳性,被通知要求核酸复查后却拒不配合,而是乘坐动车再转乘坐网约车回家了。
图2:上海本轮疫情走势
但值得庆幸的是,奥密克戎虽然传染性强,但致病毒性降低。很多感染者只是被影响到上呼吸道,基本没有症状。而且,上海是全国疫苗接种率最高的城市之一。因为疫苗的保护,人体免疫力增强,即便被感染,也绝大部分是轻症和无症状感染。
问题图2是2月15日-4月5日上海每日新增本土确诊和无症状人数。从图中可以看出3月15日之后,数据不断上升。将本土确诊和无症状人数合并为感染人数,表1给出了以2月15日为起点上海最近二十天(3月17日-4月5日)累计感染人数。例如表中3/17对应的数据表示2月15日至3月17日每日新增确诊和新增无症状人数的总和。
表1:上海最近二十天累计感染人数日期累计日期累计3/173///////////////////请根据表中数据,建立模型描述上海本轮疫情累计感染人数(表1中的数据)的变化规律,并对未来两日内每日新增感染人数作出预测。
模型复合增长对于一个相对封闭的城市,每日新增的本土感染病例都是被之前本土已感染病例传染的。当累计感染人数远小于城市总人口时,每日新增病例应该与累计感染人数成正比。上海目前累计感染人数约为10万,远小于上海万的总人口。因此,上海每日新增病例应该与累计感染人数成正比,或者说上海病例日增长率应该近似为一个常数。这就和固定利率的银行存款问题类似,如果我们用
表示3月17日累计病例数(相当于本金),那么3月18日的累计病例数(
)
3月17日累计病例数(
)
3月18日新增病例数(相当于利息)。当天新增病例数等于前一天累计病例数
增长率(相当于利率),所以有
其中
为日增长率。同理可以得出第
天累计病例数为
接下来我们需要验证上海的病例日增长率
是否为常数。第
+1天病例日增长率可表示为
其中
表示第
+1日新增病例数。例如,3月18日的增长率为
其中和分别为3月17、18日的累计病例数(表1)。类似地,我们可以计算得到上海近二十天累计感染人数日增长率(表2)。
表2:累计感染人数日增长率日期日增长率日期日增长率3/17--3/./././././././.3/././././././././././.从表中不难看出,日增长率基本在0.2上下波动,并且基本在0.15-0.25之间,因此可以认为日增长率为常数。我们可以用3月18日累计病例数作为模型参数
的估计,用表2中19个日增长率的算术平均值作为模型参数
的估计:
由此可得图3中的红色预测曲线,预测结果均方根误差(预测值与观测值偏差的平方和与预测次数比值的平方根)为。预测给出4月6、7日的新增病例数分别为39和。
图3:复合增长模型预测结果
当然也可以根据表1中的数据直接应用最小二乘法拟合得到模型参数:
=,
=0.拟合曲线见图3中的蓝线,相应的预测结果均方根误差为。预测给出4月6、7日的新增病例数分别为97和。
指数增长复合增长模型假设了病例数量是以天为单位增长的,但实际上病例数量在每一小时、每一分、每一秒都在增长。如果将每天的时间细分为
等份,累计病例数在每一个等份时间的增长率为
=
,则第
天的累计病例数可表示为
将上式方括号中的项定义为函数
=
。随着时间无限地细分,
将无限增大,而
将趋于0。
随着
的变化过程见表3。
表3:的根限值
0.10..0.00.
2.....
从表中不难发现,当
逐渐趋于0时,
也逐渐趋于一个常数:
这个常数被定义为自然底数e=2.
,它是重要程度不亚于
的另一个无理数。因此,第
天的累计病例数可表示为以下指数增长模型:
这就为什么像大模头这样的专家们在谈论疫情时,经常会说“呈指数级增长”。实际上,现实世界中很多数量的增长都呈现类似的规律,比如生物种群的数量增长。通过最小二乘法,我们可以拟合出指数增长的模型参数:
=,
=0.。拟合曲线见图4,相应的预测结果均方根误差为。
图4:指数增长模型预测结果
预测给出4月6、7日的新增病例数分别为97和。指数增长模型给出的结果与最小二乘法拟合得到复合增长模型的结果一致。这说明在短时期内,复合增长模型已经足够说明问题。这一点很重要,特别是当我们在中学知识范围下讨论这个问题时。
结论自年3月中旬以来,上海每日新增本土确诊和无症状人数一路飙升。根据3月17日-4月5日累计感染人数,本文先建立了复合增长模型。在此基础上,考虑到病例的增长并不是逐日进行的,本文将每天的时间进行细分。在无限细分的极限条件下,复合增长模型转变为了指数增长模型。指数增长模型给出的结果与复合增长模型的结果一致:上海的累计新冠病例数量正在以20%的日增长率增长,4月6、7日的新增病例数将达到1.8万和2.2万左右。
需要注意的是,本文的模型预测结果是建立在4月6、7日不出现拐点(即累计感染人数仍按照之前二十天的规律增长)的假设前提下得到的。
参考资料[1]国家和省市卫健委,疫情实时大数据报告,/4/5:
转载请注明:http://www.0431gb208.com/sjszlff/774.html