1. 引言
众所周知,世界正在经历一场阿片类药物危机,美国作为超级大国影响尤为严重,此次危机不仅仅对美国社会健康问题带来巨大的负面影响,同时还会对美国现有的法律体系问题和美国经济的重要部门带来不可磨灭的影响和复杂深刻的挑战。根据美国毒品管制署所提供的数据报告,美国联邦、州地方毒品犯罪的主要是受麻醉止痛药(合成阿片)和吗啡的影响,而美国各州之间受到这两类阿片药物的影响程度是因地而异的。
更详细地说,如果阿片类药物的危机蔓延到美国的所有横截面,将会对美国政治、经济局势带来十分恶劣的影响。基于上诉阿片药物危机的问题,本文立足美国的肯塔基州、弗吉尼亚州、西弗吉尼亚州、宾夕法尼亚州、俄亥俄州五个州的各个县阿片类药物影响,根据美国毒品管制署与美国人口普查所知的社会经济数据,综合得出五个州的阿片药物在未来的影响情况,并对美国阿片药物的控制提供控制性的积极建议。
因此,我们用五大州的Drug Report数据占五个州总Drug Reports数据的比例建立时序散点图,对筛选后剩下州的时序散点图进行四次拟合得到拟合曲线,再通过截距判断最早发生的州位置;然后用该州的各县Drug Report数据占总县Drug Report的比例,用制图软件建立三维曲面图,通过峰值的高低判断县位置。研究阿片药物的传播和扩散,确定阿片药物的最开始的传染源以及阈值水平,利用传染并传播模型SIR对问题进行分析。本文将对整个阿片药物的传播系统抽象成为复杂网络,在该网络上,节点表示美国五个州各个县的个体,边表示美国五个州各个县之间的相互联系,表示节点之间的相互关系,由于阿片药物使用的传播方式与传染病传播及其相似,因此建立SIR模型针对复杂网络上阿片药物的进行科学的演化,从而预测未来阿片药物发展的具体位置与时间。考虑不同阶段网络上阿片药物的传播及传播阈值的影响,引进相对感染力和传统时间尺度的概念,分析无标度网络的相对传染里的阈值小于小世界的网络阈值,从而设置合适的阈值来控制传播向着有利于降低阿片药物传播的方向发展。
根据美国人口普查所提供的社会经济数据,分析数据对阿片药物使用和上瘾的趋势是否存在影响,综合考虑社会与个体的关系以及个体处在不同阶段的差异。受到文献 [1] [2] [3] [4] [5] 启发,我们基于SIRS提出的阿片药物多感染传播阶段的SInR模型。感染阿片药物的个体在不同的阶段具有非均匀的感染力,同时引进感染尺度与恢复尺度的概率,从复杂网络结构、网络规格和相对的感染力方面进行仿真研究。得出个体有着不同传播的影响因素,并研究多因素在阿片药物传播网络行为规律,确定有效性和参数,最后得出控制合适的阈值可以有效的控制阿片药物传播的最终影响力,进而得到相应控制策略。
2. 假设和名词符号说明
2.1. 基本假设
均匀性假设,美国各大州各个县的个体自身不存在差异,个体对于阿片类药品的接触、使用和上瘾的概率是等效的,阿片药物的传播是等概率接触、传播。假设本文中阿片药物的传播方式对一些实际的传播特性和网络结构的影响是可以忽略不计的。假设美国总人数是不变的,并且感染病毒的时间尺度远小于个体的生命周期,从而不考虑个体的出生与自然死亡。假设美国国家毒品管制署提供的犯罪实验室向缉毒局提供的报告和美国人口普查局提供的数据是正确无误的,与实践情况没有任何的偏差,并且数据和报告具有权威性和说服力。国各个州及各个县之间合成阿片中药物报告数较小的合成药物对推测五个州及各个县最早出现阿片类药物的影响是可以忽略不计的,最早出现阿片类药物的州及其具体位置可根据阿片药物报告数中数目在百以上来推测。
2.2. 基本假设
阿片类药物:缓解疼痛的一种药物,容易上瘾。
吗啡:一系列吗啡类毒品的总称,同时也是一种非法生产、容易高度上瘾的一种镇痛药物吗啡。
合成类阿片:由人工合成的阿片。
2.3. 符号说明
为阿片药物在基本复杂网络中的感染率;
为阿片药物在基本复杂网络中的恢复率;
为阿片药物在基本复杂网络的传播效率(传播临界值即阈值);
为阿片药物的易感者与免疫者的转化率;
为阿片药物在第n阶段时传播感染的感染率;
为阿片药物在第n阶段时传播感染的恢复率;
为阿片药物在S状态的易染者个体和处于I状态的邻居的j权重。
3. 模型建立及求解
3.1. 时序散点模型与三维曲面模型
用2010~2017年美国Drug Reports数据的比例建立时序散点图,然后对筛选后剩下州的时序散点图进行四次拟合得到拟合曲线,再通过截距判断最早发生的州位置;然后用2010~2017年该州的各县Drug Report数据占总县Drug Report的比例,用制图软件建立位置、时间、比例三维曲面图,通过峰值的高低判断县位置。并通过图形分析随着时间的推移,报告的五个州及其县之间的合成阿片和吗啡事件(病例)的特点。
首先,分析吗啡时序散点图1与图2,可以看出KENTUCKY、VIRGINIA、WEST VIRGINIA三个州Drug Report数据占五个州总Drug Reports数据的比例每年都相对比较低,大多数仅占10%以内,而PENNSYLVANIA和OHIO两州Drug Report数据占五个州总Drug Reports数据的比例高达35%以上。因此,从概率上来讲KENTUCKY、VIRGINIA、WEST VIRGINIA三个州不太可能为最早发生吗啡使用的具体位置。由拟合曲线图可以清晰的看到PENNSYLVANIA的截距在0.5左右,而OHIO的截距在0.35左右,两者相差0.15,从而可以确定PENNSYLVANIA州为最早发生吗啡使用的位置。
![](//html.hanspub.org/file/6-2620873x17_hanspub.png)
Figure 1. PENNSYLVANIA morphine time series discrete fitting map
图1. PENNSYLVANIA吗啡时序离散拟合图
![](//html.hanspub.org/file/6-2620873x18_hanspub.png)
Figure 2. OHIO morphine time series discrete fitting map
图2. OHIO吗啡时序离散拟合图
然后运用相同的方法,由于VIRGINIA和WEST VIRGINIA两个州Drug Report数据占五个州总Drug Reports数据的比例与PENNSYLVANIA、OHIO、KENTUCKY三个州相比要低的多,所以仅对PENNSYLVANIA、OHIO、KENTUCKY三个州所对应的散点图进行四次拟合,如图3,图4,图5所示。
因此可以确定五大州中最早发生的为PENNSYLVANIA。再将PENNSYLVANIA州的各个县建立位置、时间、比例三维立体曲面模型,如图6所示,判断PENNSYLVANIA的具体位置。
![](//html.hanspub.org/file/6-2620873x19_hanspub.png)
Figure 3. KY discrete fitting diagram of synthetic opioid time series
图3. KY合成阿片时序离散拟合图
![](//html.hanspub.org/file/6-2620873x20_hanspub.png)
Figure 4. PENNSYLVANIA discrete fitting diagram of synthetic opioid time series
图4. PENNSYLVANIA合成阿片时序离散拟合图
![](//html.hanspub.org/file/6-2620873x21_hanspub.png)
Figure 5. OHIO discrete fitting diagram of synthetic opioid time series
图5. OHIO合成阿片时间序列离散拟合图
![](//html.hanspub.org/file/6-2620873x22_hanspub.png)
Figure 6. PENNSYLVANIA three-dimensional surface model
图6. PENNSYLVANIA三维立体曲面模型
3.2. 小世界网络模型与无标度网络模型
基于上诉分析,将阿片药物传播系统抽象成基本复杂网络。对于基本的复杂网络,小世界性和无标度性,其中小世界网络具有较短的平均距离和较大的聚集系数;同时在基恩复杂网络中大部分的节点只有少数点链路,同时某些节点却拥有着与其他节点大量的链接,表现在无标度性上的就是幂律分布,综合考虑阿片药物传播系统类似网络的增长特性、同时系统中新节点更倾向与那些具有高度的“大”节点相连接的优先连接特性。因此,阿片药物传播系统类似的基本复杂网络主要包括小世界网络与无标度网络两种网络模型。
复杂网络中的多种形式的传播现象都可以看作是服从某一规律的传播行为。实际的复杂网络系统的构成元素众多、组织因素庞大、这导致信息在网络上的传播规律过程是十分复杂的,因此还需构建小世界网络模型和无标度网络模型来解释传播行为规律并且解决文中的相关问题,并将本部分中阿片药物的扩散传播抽象为感染传播模型。
3.3. SIR模型
SIR模型是用来描述疾病传播过程中被感染个体在治愈后具有免疫力,因而成为健康个体并且不会被传染,或者难以治愈而死亡,从而移除传播系统不再感染其他个体。
SIR模型中,文中所要研究的美国五大州各个县的个体主要被分为三类:易染者S (易染个体)、感染者I (染病个体)、恢复者R (免疫个体)。在初始条件下,通常认为系统中一个或者少数几个处于感染状态,其余个体处于易感染的状态。一般情况下,SIR模型传播规则如图7所示。
记
分别表示易感者、感染者、免疫者在美国五大州各个县总个体所占的比例,
表示各初始状态下人数的初值,且均为负数,因此得到了SIR模型的微分方程如下:
因此,在基本复杂网络中的k节点条件下,设
分别表示在k节点时的易感者、感染者、免疫者在感染点与恢复点在t时刻所占的各个比例:
此时我们取用k节点的传统时间尺度关系式:
因此SIR模型的微分方程可变化可得:
初始情况下,当
时,可取用
,可得
,
并构建
辅助函数:
进而可以推导:
即:
因此本文令
,则可以求出:
因此在本题中只需满足
,即可知道:
而根据题目中的已知条件,我们必须满足:
因此可知:
从而可知,在问题一条件下的阈值是:
综上所述,在基本复杂网络中,阿片药物的传播存在一个有限的正的传播临界值
,如果在阿片药物传播过程中,阿片药物的有效传播率大于这个临界值,传播感染将会持续稳定的进行下去,最终情况下,阿片药物的传播将处于一个稳定的平衡状态,即称该基本复杂网路处于激活相态;如果阿片药物的有效传播率小于这个临界值,那么感染源的个数将会受到限制,阿片药物的传播无法大范围的进行传播,此时的基本复杂网络处于吸收相态。因此在阿片药物传播的小世界网络与无标度网络中,存在阈值
,将阿片药物传播网络的激活相态和吸收相态分开。
基于上诉对阈值
的确定,选择合适阈值,针对阿片药物传播系统利用SIR模型演化5年,本文发现,在将来的阿片药物传播感染在美国PENNSYLVANIA 州持续传播,最后并趋近平稳状态。如图8,图9所示。
![](//html.hanspub.org/file/6-2620873x49_hanspub.png)
Figure 9. SIR evolutionary synthetic drugs
图9. SIR演化合成药物
3.4. SIR模型
基于第二部分中美国人口普查所提供的社会经济数据,本文在第一部分SIR模型的基础上,考虑在阿片药物传播过程中,免疫者可能会丧失免疫能力转化成为易感者,因此本文对模型进行优化,引入SIRS模型,考虑转化概率
的影响。
此时,阿片药物传播系统的传播机理如图10所示:
![](//html.hanspub.org/file/6-2620873x51_hanspub.png)
Figure 10. Dissemination mechanism of opiate drug delivery system
图10. 阿片药物传播系统的传播机理
其微分方程表达式:
由于在阿片药物传播系统中,个体与个体之间的联系并不是单一的相互联系,同时对于个体所在的状态阶段的不同也会有着不同的感染力。
3.5. SInR模型
由于在阿片药物传播系统中未考虑个体之间的差异,而根据美国人口普查所提供的社会经济数据,因此美国五个州各个县中个体及其周围的邻居接触的概率不同从而导致阿片药物传播的概率也不同,同时由于社会经济因素的影响导致美国五个州各个县的个体自身的免疫能力的不同从而导致阿片药物传播概率的不同。
在阿片药物传播系统中,个体在每一个阶段都具有不同的感染传播概率,基于此,本文建立具有多感染传播阶段的SInR模型进行求解,SInR模型的工作机理如图11:
其微分方程表达式为:
在上诉公式中,
表示阿片药物在第n阶段时产波感染个体的比例,且满足下列条件:
若此时在S状况下,阿片药物传播系统复杂网络中,易染者个体i将以概率为
被阿片药物感染成I状况(即由S变成
),而满足权重关系式:
根据上诉权重关系,可以确定阿片药物的状态阶段,利用SInR模型不断进行演化,最终发现直到
,阿片药物的传播系统保持稳定状态,即假设具有免疫系统的美国五大州各个县的个体也不会被感染。
综上所述,根据SInR模型的演化,本文发现阿片合成物易于获得和使用广泛是滥用现象出现的主要因素,过去的数据中可以明显的发现18~40岁的年龄段的人群占有阿片上瘾的数量不断增加。现有大量的含有阿片的药物流通在市场上,现行的制度没有明确上瘾的发生条件,处方阿片类药物使用的广泛增加和随之而来的误用造成的最紧迫问题。现在市场上有大量且容易获得的吗啡加剧了对于阿片的上瘾。由于合成阿片的易获得使阿片上瘾有着低龄化的趋势。
因此,综合上诉的分析,我们制定以下策略:
1) 首先加强对含有阿片药物的管制,控制每个患者定期的用量避免上瘾。
2) 严厉打击市场上所有的毒品,防止人们受到毒品的影响。
而根据上诉所建立的模型,计算确定在该策略下阿片药物的有效传播效率,将该效率与第一部分与第二部分时阿片药物感染传播的临界阈值进行对比,可发现,我们制定策略后的传播效率低于临界阈值,即阿片药物传播抽象的基本复杂网络呈现吸收相态,则说明我们策略有效。若要使制定策略成功,必须使制定策略后阿片药物的有效传播效率低于临界阈值
即可,即
即可。
4. 结论
本文中所建立的模型充分利用了阿片药物在传播过程中与传染病传染相类似的特点以及阿片药物的传播规律,抽象基本复杂网络,充分结合SIR模型与SIRS模型的传染模型特性,并在多感染源阶段的条件下,建立与之对应的SInR模型。文中所建立的模型将阿片药物传播时所经受的多种具有不同阶段性和不同传播效率的产波感染力的影响,同时利用阿片药物传播阈值的影响,引入阿片药物的相对感染传播能力与传播时间尺度的概念,对整个阿片药物传播系统进行对比仿真分析,实验数据与实际中阿片药物传播与传染的影响相一致。
同时,文中所建立的模型并未确定阿片药物传播系统抽象处的复杂网络的最大适用范围,以及确定实际运用的网络模型,可能与实际阿片药物传播系统网络存在一定差距。因此,基于本文模型,构建具有多边权重的复杂网络,充分考虑在阿片药物传播过程中个体与周围环境、个体与个体之间的影响关系。
基金项目
同济大学浙江学院第七届教改项目(项目编号:0118037)。
NOTES
*通讯作者。