1. 引言
近年来,为了更好地研究传染病在空间上的演化情况以及更好地预测传染病的流行趋势,反应扩散传染病模型是众多学者广泛关注的研究内容,并已取得了很多成果。作为反应扩散传染病模型的一个重要研究内容之一便是行波解,由于空间平移的不变性,行波解不但可以很好地揭示传染病的传播过程,震荡现象,还能预测未来的传染规律和发展趋势,具有十分重要的实践意义和理论意义。其中,行波解的存在性和唯一性的证明已取得了丰硕的成果,可参见文献 [1] - [8] 。除此之外,行波解稳定与否可以直观反应传染病的传播形态能否保持不变且不受干扰,对传染病的预防和控制也具有很强的指导意义。因此对行波解稳定性的研究也是当前的研究热点和难点之一。
在研究反应扩散方程行波解稳定性时,代表性的研究方法就是Mei等人对于加权能量法的研究,他们在文献 [9] 中利用加权能量估计的方法和比较原则证明非局部Nicholson非线性方程行波解的全局指数稳定性。Yang等人分别在文献 [10] 和 [11] 中采用加权能量法分别研究单调系统和非拟单调系统的行波解的稳定性。Hsu等人的文献 [12] 针对一类带时滞的反应扩散模型,研究其行波解的存在性和稳定性,他们的研究结果表明围绕行波解的初始扰动属于适当的加权Sobolev空间时,模型的行波解是指数稳定的。Li和Weng [13] 也采取同样的加权能量法和比较原理研究SIS传染病模型的行波解在大波速情况下的指数稳定性。Su和Zhang [14] 研究一类带时滞和非拟单调的反应扩散方程的行波解的稳定性,采用傅里叶变化再结合加权能量法建立恰当的能量方程,证明了无论是在单调还是非单调的情况下,方程的行波解都
以收敛率
在
空间上是全局指数稳定的,其中
是某个正常数。更多相关文献还可以参
见 [15] [16] [17] [18] 。
最近,将空间扩散系统进行空间离散化得到扩散性格微分系统引起了研究者越来越多的关注。格微分系统是偏微分系统的离散化形式,与空间连续系统相比,其动力学行为更加复杂,可以更好的刻画实际问题,并且在数值计算,参数估计等方面具有明显的优势。目前研究者大部分都是针对三维以下的离散系统进行行波解稳定性方面的研究。Zhang和Tian [19] 建立一个一维格上的两物种Lotka-Volterra竞争系统,并得出模型大波速行波解的全局稳定性。Yang和Zhang [20] 研究非单调离散扩散方程的行波解的稳定性和唯一性,利用加权能量方法和非线性Halanay不等式,证明当初始振动很小的时候,所有的行波解是时间指数稳定的。Hsu和Lin [21] 研究离散的SIS传染病模型的行波解稳定性,借助加权能量法和比较原理得出当初始扰动在加权Sobolev空间任意大时,所有非临界波速行波解的全局指数稳定性。Su和Zhang [22] 还研究了一类格上的三种群Lotka-Volterra竞争系统波解的稳定性,利用类似的比较原理和加权能量的方法,得到大波速波前解的指数渐近稳定性。Zhao等人在文献 [23] 中研究一类格上的考虑年龄结构的二维人口模型,研究行波解的存在性和唯一性,其后利用挤压技术,建立非最小波前解的稳定性。其他类似的参考文献可参看 [24] [25] [26] 。
本文针对一类具有多种传播途径的霍乱传染病进行研究,这是一种由于摄入受霍乱弧菌污染的食物或水引起的急性肠道感染疾病,如果未经合理诊治,感染者可在短期内死亡。19世纪时期,源自印度恒河三角洲的最初宿主使霍乱在世界各地的蔓延。霍乱曾经在全球范围引起过七次世界大流行,之前的六次大流行使得各大洲百万人口因此丧生。最近的第七次霍乱大流行于1961年始于南亚,1971年蔓延至非洲,1991年至到美洲。据世界卫生组织估计,全球每年因霍乱导致300万至500万例感染者,有28,000到142,000人因此死亡,不仅是国际卫生条例规定的国际检疫传染病之一,也是我国规定的甲类传染病之一。近期比较严重的两次霍乱疫情,一次是于2008年~2009年在津巴布韦爆发,引起超过9万八千例感染者,最终统计有4200人死亡,尽管在国际机构的帮助下,疫情于2009年结束,但于2018年又一次卷土重来,造成2000多例感染者;另外一次2010年海地在地震之后也爆发严重的霍乱疫情,影响了82万多人,造成9792人因此丧生。在最近几十几年间,霍乱还陆陆续续在各国各地爆发,尤其是一些不发达或者发展中的国家,例如刚果(2008年),伊拉克(2008年),越南(2009年),尼日利亚(2010年),墨西哥(2013年),南苏丹(2014年)和索马里(2017年)。
霍乱因其传播方式与众不同而受到研究者的广泛关注,其具有多种传播途径,不但通过人与人之间直接传播,还通过人与被污染的水源进行非直接传播。Codeco [27] 改进了标准的传染病SIR模型,首次计算在水源环境中霍乱弧菌的浓度,但该模型只假设不洁水源为唯一感染途径。Hartley等人 [28] 随后在Codeco工作的基础上,充分考量到了霍乱多途径传播的特点,将霍乱病菌分为高感染阶段和低感染阶段,结合两个新的环境元素构建一个高维霍乱传染病模型。Mukandavire等人 [29] 再将Hartley的模型进行一些简化,采用非线性发生率来描述人与被污染的水源之间的传播,重点刻画传播过程中的饱和效应。Tien和Earn [30] 构造并分析一个具有多种传播途径的高维霍乱模型,对该模型模型进一步拓展分析和数值模拟。Shuai和Driessche [31] 依照病毒毒性的强弱把水源中的病毒分为不同类型,同时也把受感染的人群按照对应的病毒传染性的强弱划分为有限个类型进行分析。Liao等人 [32] 研究带有非局部扩散项的霍乱传染病模型行波解的存在性问题,利用Schauder不动点定理,构造上下解得到行波解的存在性,其次构造Lyapunov函数结合Lebesgue控制收敛原理,得到行波解在
处的渐近行为。更多的相关文献可以参见 [33] [34] [35] [36] [37] 。
本文拟针对一类格上的高维霍乱传染病模型的行波解进行稳定性分析。据我们所知,目前针对三维以上的格微分系统的行波解稳定性问题,结论还比较少。本文的主要结构如下。首先建立格上的高维模型系统,给出相应的初值问题的解的存在唯一性和比较原理,接下来结合加权能量方法和比较原理讨论解的能量指数衰减估计,再利用嵌入定理证明大波速行波解的全局渐近指数稳定性,最后是全文总结。
2. 霍乱传染病模型
霍乱传染病的传播比较复杂,不但通过人与人之间进行直接传播,还通过人与不洁水源之间进行非直接传播。我们在文献 [32] 的基础上,建立带扩散项的霍乱传染病模型如下
(1)
(2)
(3)
(4)
(5)
和
分别代表易感者,接种者,感染者和移出者在t时刻x处的密度,
表示水源在t时刻x处的霍乱病菌浓度。模型中参数Λ为输入率,
和
分别表示环境与人之间传播和人与人之间传播的传染率系数,
为霍乱病菌的衰减率,
为死亡率,
为染病者的复原率,
为病原体的去除率。所有的参数都为正数
为正扩散率系数。
注意到系统方程中R的独立性,为了简化计算,在后面的分析中只考虑方程组(1)~(4)即可。显然模型(1)~(4)具有无病平衡点
,以及地方病平衡点
,且表达式为
求雅可比矩阵并代入初始值,得基本再生数为:
。
3. 行波解稳定性分析
本节对系统(1)~(4)进行空间离散化得到如下四维格微分模型
(6)
(7)
(8)
(9)
其中
。则模型系统(6)~(9)的解
为行波解。
再令
,
为波速,将
代入系统(6)~(9),可得如下的行波方程
(10)
(11)
(12)
(13)
且系统的行波解要满足下面的边界条件
以及
接下来本节利用加权能量法讨论系统模型(10)~(13)的行波解是指数稳定的。在进行稳定性证明之前,先介绍一些必要的记号。
设
是一个一般常数,
是特殊的一些常数。设I是一个区间,特殊情况下取
。
为定义在I上的平方可积函数区间。
为Sobolev空间,是定义在区间I上的L2-函数
的集合,其导数形式
且属于
。
为L2-空间上赋予加权函数
的函数空间,其范数定义为:
。
为赋予加权函数的Sobolev空间,其范数定义为:
令
是一个正常数,Β是一个Banach空间。再定义
为
上Β-值连续函数空间。
为
上的Β-值L2-函数空间。再定义加权函数如下:
,
其中
,
为一个充分大的数使得
。并令
以及
接下来证明本文主要的稳定性定理如下
定理1 对任意给定的行波解
,
则存在一个
以及
使得
。
以及
,
其中C和
为正常数。
为方便起见,再定义
其中
,
,
。
本节分以下三步来证明定理3.1。
第一步.
收敛到
。
引理1 (比较原理)假设
,
分别是系统(10)~(13)带有相应初值
和
的解,如果对任意的
满足
那么对任意的
,都有
通过计算,
满足
(14)
且满足初始条件
,
,
。
易知
,则有
引理2 存在一个正常数
使得
.
证. 首先将式(14)中四个等式分别乘以
,
,再由
可得
(15)
其中
由Cauchy-Schwarz不等式
知
和
,
.
将不等式(15)关于
和t在
上积分,由
可知,在无穷远处
,
从而得到
(16)
将(16)中不等式相加可得以下结论:
.
其中
引理3 对任意的
存在一个
,
,在满足下列假设条件(17)的情况下
(17)
使得对所有的
都有
,
。
证. 分以下四种情况进行证明
Case 1当
时,
此时可计算
Case 2.当
时,
,以及
此时可计算
Case 3.当
此时
,
Case 4.当
此时
,
,可计算出
综上即可得出当
充分小时,总存在足够小的正常数
,使得
。证毕。
进一步分别对(15)关于
进行微分得
(18)
再将(15)两边同时乘以
,
得
(19)
引理4 对任意的正实数
,成立估计式
证. 再次利用Cauchy-Schwarz不等式从而得到
.
以及
.
由
可知
,所以在
上关于
和t积分可得
(20)
.
将(20)中不等式相加得
.
以及
利用和引理4相同的方法可以得出如下引理(此处省略证明):
引理5 当
充分小时,存在一个正常数
使得对所有的
,都有
。
引理6 对任意的
,存在一个
和一个正常数
使得
再对不等式(15)右边的最后一项进行估计。利用行波解
的性质,可知对任意的
,
都是有界的。由Young-不等式
,当
时,存在一个正常数
,有:
则
结合上述结论,容易得到以下一致先验估计。
引理7 对任意
,存在一个正常数
和
,
,使得
由标准的Sobolev嵌入不等式,对充分大的
,
。但是对于任意区间
,可得
。于是可得下面的
估计。
引理8当
时,对充分大的
,
存在正常数
,使得
为了将引理8所得到的结论扩展到整个区间
上的估计,故需要证明在
处
的收敛性。引理9 对所有的
,总存在一个正常数
,使得
。
证. 首先易知
,以及
,再基于
的有界性,当
时,由上述不等式变为
(21)
将不等式(21)相加可得
其中
,在
上积分,存在一个常数
使得
。
因此,必存在一个正常数
使得
结合引理8和引理9,可证明出当
时,
的
收敛性。
因此,可得结论
其中
。
第二步.证明
收敛到
。
用类似于第一步的方法即可马上证明出第二步的结论,故此省略。
第三步.证明
收敛到
。
由挤压原理,易得以下结论:
引理 10存在一些正常数C,有
4. 总结
本文研究一类离散扩散霍乱传染病模型行波解的稳定性问题,是liao和Yang在文章 [32] 中对霍乱模型行波解的存在性进行研究之后,作进一步延续。本文主要证明模型的大波速行波解是关于时间t的全局指数渐近稳定的,即行波解能保持稳定不受外界环境的干扰。这些有效的理论结果不仅对霍乱传染病模型适用,还可以推广到所有的水源性传染病模型行波解的研究中。而小波速行波解的稳定性研究以及在非拟单调情况下的稳定性,都将是以后的工作重点之一。
基金项目
重庆市基础研究与前沿探索项目(cstc2020jcyj-msxmX0394)。
NOTES
*通讯作者。