1. 引言
自2022年以来全国范围内陆续出现了多次较大规模疫情爆发事件。新冠疫情的爆发对全球公共健康、生命安全和经济发展都造成了巨大的危害。在大规模疫情爆发期间,我国主要采取封闭式管理方式来实现疫情的快速清零,因而对疫情期间各类人群的科学并有效的管理就显得尤为重要。虽然自新冠肺炎爆发以来,人们对于引发疫情的2019新型冠状病毒(COVID-19)的科学认知也逐步加深 [1],应对手段逐步调整,信息掌控的逐步丰富,但是大部分地区管理者尚未遇见如此大规模爆发的疫情及不同地区疫情爆发时长、人口规模与地理位置等的差异性,众多地区仍未形成较为科学的统一管理模式,使得疫情的清零周期普遍较长。应急物资配置决策收到传染病疫情时空分布、感染规模等因素的影响,而这些因素很大程度上取决于一曲采取的疫情防控措施,进而形成不同的动力学传播。传染病的防控和救治一直是我国人民普遍关心的一类问题,而对疾病的发病机制、原理以及传播特征等规律进行理论研究和定量分析 [2],对防治工作具有十分重要的意义。截至目前,国内学者对传染病防控问题的背景以及所面对的挑战都给出了自己的分析与设计,但现有研究仍有一定的不足,仍然有需要进一步进行解答的问题,例如:传染病卫生治理到底能否治理,用何种手段治理。除此之外,当前研究主要集中于定性研究,比较缺乏对问题所进行的定量研究,同时对文化及社会认知程度的研究也比较缺乏。由于传染病的特殊性,其传播过程不能通过进行实验得以重现,因此,通过数学建模进行研究就显得非常必要。数学模型以数学语言为基础,通过对实际问题进行概括,抓住主要影响因素,忽略次要影响因素,建立互相联系的数学公式,并通过求解进行定量分析,进而达到指导现实工作的目的。本文基于长春市内生活物资投放点数量的合理性并对其进行优化,建立合适的模型并带入已知数据进行求解,模型的应用可以对今后疫情防控提供有利帮助。
2. 投放点数量合理分配问题
2.1. 基本信息分析
首先根据实际数据中长春市9个区的交通网络与主要小区(1409个)的散点分布见图1,结合小区位置与路口结点信息,分析小区分布与道路分布的特点 [3]。
![](//html.hanspub.org/file/53-2622935x8_hanspub.png?20140215005915744)
Figure 1. Distribution map of traffic network and main residential areas in 9 districts of Changchun City
图1. 长春市9个区交通网络与主要小区分布图
对于各区投放点数量的合理性,我们结合各区的小区栋数、小区户数、小区人口数、街道数目、新冠疫情人数以及无症状感染者人数这些参数指标分析,可知这些指标属性为效益性指标,即属性值越大,需设置投放点的数量越多。由此可知,目前已知信息中投放量并不合理,各区生活物资投放点数量与隔离人口总数并非正相关。所以针对该问题,依据参数指标对区域进行重新评估排序是有必要的。
2.2. 指标参数的说明
我们构建指标体系来对各区进行评估得分。各区的参数指标表见表1所示。
最优解为
。
最劣解为
。
以二道区为示例,数据预处理,运用向量规范化进行以下变换
.
![](Images/Table_Tmp.jpg)
Table 1. The parameter index table of each district
表1. 各区的参数指标表
.
由此得到向量规范化属性矩阵,从而消除量纲影响。
2.3. 层次分析权重的确定
接下来通过层次分析法确定权重:
![](//html.hanspub.org/file/53-2622935x13_hanspub.png?20140215005915744)
Figure 2. Distribution map of distribution points in each district
图2. 各区投放点安排分布图
首先如图2,我们确定目标层为各投放点数量的合理安排,因素层则分别为6项参数指标。构建准则层(因素层)对目标的成对比矩阵A为:
成对比矩阵A的最大特征根为
,矩阵的阶数
。
其特征向量(权向量)为:
.
归一化后结果为:
.
一致性指标为:
.
随机一致性指标为:
一致性比率为:
.
2.4. 最终得分求解
最终得到二道区的综合得分为0.274232274,同理可以计算出其他8个区的综合得分,排序与分配的结果见表2所示。
![](Images/Table_Tmp.jpg)
Table 2. Ranking table of the scores of the nine districts
表2. 九个区的得分排序表
3. 多目标规划模型
为达到投放点数量的合理规划目标,本文以投放点数量为决策变量,构建如下的多目标规划模型,同时运用到聚类分析与集群的思想 [4]。
3.1. 目标函数
以储备数量最小化目标作为目标函数,如公式(1)所示:
, (1)
其中,
为t天的投放点i的数量,
为投放点i的储备数量,n为参与规划的目标个数,下同 [5]。
选址数量最小化目标如公式(2)所示:
, (2)
其中
是投放点i的选址系数,下同。
人均获取蔬菜包里程最小化目标如公式(3)所示:
, (3)
其中,
表示投放点i的人均获取蔬菜包里程系数,下同。
3.2. 约束条件
储备数量保障天数约束如公式(4)所示:
, (4)
其中,10代表储备数量能保障各区天数的约束值。
人均分发数量约束如公式(5)所示:
(5)
其中,31代表COVID-19疫情持续的时间,
表示各区人数约束值。
市场价格增长约束如公式(6)所示:
(6)
其中,
表示t天的生活必需品市场价格p的增长率。
下限约束如公式(7)所示:
(7)
假定在规划期内,不会有新的目标参数产生,也不会有旧的目标参数淘汰,为保证数学模型稳定运行,对各投放点的数量进行约束 [6]。
3.3. 模型转换
式5是非线性表达形式,需要将其转换成线性表达式,以进行求解。转换过程如公式(8)所示:
(8)
本题最终构建的多目标规划模型如公式(9)所示:
(9)
充分考虑未来疫情、自然灾害等特殊事件,对于政府储备物资和大规模物资分拣场所(商超)的位置与数量、规模进行合理规划 [7],得到最优的选址数量、规模及其潜在的备用场所位置结果,我们结合多目标规划、聚类分析与集群的思想,由于表格过长文中只展示一个区的结果作为实例,宽城区的50小区的最优选址数量、规模及其潜在的备用场所位置结果见表3所示。大型物资中心数量结果见表4所示。
![](Images/Table_Tmp.jpg)
Table 3. Statistical table of the number and scale of the optimal site selection of 50 communities in the broad urban area and the location of potential standby sites
表3.宽城区50个小区的最优选址数量、规模及其潜在的备用场所位置统计表
![](Images/Table_Tmp.jpg)
Table 4. Large material center number results
表4. 大型物资中心数量结果
4. 物资供应任务具体预案
4.1. 有序网络图的建立
本文依据人口密度、需求量等参数指标区别对待发放生活物资的任务,这样才能节省人力、物力的同时保证投送效率。首先我们对假定的网络格点上游(第一结点)采用枚举法,假设每个区有一个中心点,这个点到该区所有小区的距离之和为最小值,此时投送效率最高,以小区人数作为权重系数(这里将聚类准则作了调整),此时小区的人口密集度决定了中心点的偏移方向,这样可以确定9个大型站点的坐标位置;对假定的网络格点的中游(第二结点),我们仍参照第二问的区域投放量作为指标进行设置,最后依据长春市9个区交通网络与主要小区分布散点图,将一个区域所有点放在一起构成一个面,那么前面所假定的网络格点都可进行一个具体的划分类别,从而确定一个区的网络格点的全部坐标。进一步绘制有序网络图。各区投放点合理分配量不变,作出9个区的有序网络图见图3~7所示。
![](//html.hanspub.org/file/53-2622935x37_hanspub.png?20140215005915744)
Figure 3. Ordered network diagram of Chaoyang District and Erdao District
图3. 朝阳区和二道区有序网络图
根据长春市9个区的有序网络图可知:根据人口密度、需求量等参数指标区别对待发放生活物资的任务,可以节省人力、物力的同时保证投送效率。
4.2. 构建配送路径和配送时间的优化模型
本文运用优化模型,这里采用单目标规划非线性模型,并结合求解非线性规划模型的信赖域方法,提出求解模型的算法。即可得最优化卡车运送物资预案,最后对照指标评价。
![](//html.hanspub.org/file/53-2622935x38_hanspub.png?20140215005915744)
Figure 4. Ordered network diagram of the open zone and the clear moon zone
图4. 经开区和净月区有序网络图
![](//html.hanspub.org/file/53-2622935x39_hanspub.png?20140215005915744)
Figure 5. Ordered network diagram of wide urban area and green park
图5. 宽城区和绿园区有序网络图
![](//html.hanspub.org/file/53-2622935x40_hanspub.png?20140215005915744)
Figure 6. Ordered network diagram of Nanguan District and steam opening area
图6. 南关区和汽开区有序网络图
![](//html.hanspub.org/file/53-2622935x41_hanspub.png?20140215005915744)
Figure 7. Orderly network map of Changchun New Area
图7. 长春新区有序网络图
4.2.1. 决策变量分析
由于是要求得最优化的配送路径和相应的配送时间,为了能够同时照顾整体的行动力,我们区别对待不同小区位置的投送任务,把有序网络图上游九各区的中心点、中心点坐标x轴和y轴作为决策变量,利用功能模拟原理(微分模拟方法与神经网络方法)来获取配送路径、配送时间与决策变量之间的关联关系,公式为:
式中:
为第i项品种蔬菜重量;
为第t天第i项的第j项的决策变量;
为第t天第i项蔬菜单价,为路径和时间的函数关系;
为第t天第i项蔬菜投放量;
为第t天第i项蔬菜投放量的配送人量;下标:
。
4.2.2. 目标函数分析
在已经算得各分类区域得投送任务得基础上,考虑中心点到区域的所有小区距离之和为最小值,使得建立最有配送路径和配送时间的优化模型,其目标函数为:
其中:r表示小区人数;d表示中心点到小区的距离。
总投送量约束:
操作总成本约束:
购买物资总成本约束:
蔬菜种类约束:
式中:a,b分别为第t天蔬菜种类的上下界;
为采买蔬菜的总成本;
为购买物资总成本上界。
根据上述分析,可以建立起单目标规划非线性模型,求解此模型就可以得出最优化的配送路径和相应的配送时间。
4.3.3. 模型算法设计
单目标规划模型是带不等式约束的非线性优化模型。在求解过程中,采用功能模拟原理获取模型中配送路径、配送时间与其对应的决策变量关联关系的建立,采用信赖域方法对模型进行求解 [7],关键公式如公式(10)~(14)所示。
首先,令
(10)
其次,求解得出
另求得
且令:
为信赖域试探步,
为信赖域半径,
为
的近似值。最后,构造信赖域二次规划模型,可得到:
(11)
(12)
根据不等式约束信赖域方法的基本思想,模型的算法步骤如下:
1) 建立各项产量与其对应的影响因素的关联关系,运用动态规划或预测方法确定初始值
,产量的上下界
及
;
2) 确定初始点以及各个计算参数:
; (13)
3) 利用信赖域二次规划模型(SQR)求解
。试探步
的确定,采用Levemberg-Marquardt法求解:确定一个阻尼因子
,利用运算式
,求解出
,若
,则迭代停止。
4) 判断
是否正定,若正定则
,继续;否则,
,转步骤3;
5) 计算
;
6) 令
则
7) 利用BFGS方法修正产生
:
(14)
令
,转步骤3;
8) 计算得出最优解
,将其代入目标函数得到最优值。
上述步骤中,对于阻尼因子
的确定存在着一些困难,但原则上
应该选取小的正数,通常情况下可以取
。但要根据实际应用而定。在迭代过程中,
为
的缩放系数,且
(通常取4)。
5. 结论
本文以长春市疫情期间的实际数据为例,构建TOPSIS综合评价模型结合层次分析法对投放点的数量和集群思想进行合理分配。TOPSIS评价模型避免了数据的主观性,能够很好的刻画多个影响指标的综合影响力度,即适用于小样本指标,也适用于多评价单元,较为灵活方便。其次通过枚举法和聚类思想对有序网络图的上游点、中游点和下游点进行坐标确定与统计,绘制出9个区有序网络散点图;通过构建单目标非线性规划模型对卡车投送预案进行优化,单目标规划建立的模型能与实际紧密联系,本文所提出的模型可以使疫情防控人员结合实际情况对所提问题进行求解,通用性,推广性较强。
NOTES
*通讯作者。