具有接種與潛伏期的SVEIR模型穩定性研究

摘" 要：該文建立一類具有接種與潛伏期的SVEIR倉室模型，利用再生矩陣的方法求出模型的基本再生數，通過基本再生數確定模型的傳播動力學。當R0≤1時，系統的無病平衡點是全局漸近穩定的；當R0＞1時，系統的地方病平衡點是全局漸近穩定的。

關鍵詞：傳染病模型；基本再生數；Lyapunov泛函；全局穩定性；平衡點

中圖分類號：R183" " " " 文獻標志碼：A" " " " "文章編號：2095-2945（2023）16-0016-04

Abstract： In this paper， a kind of SVEIR compartmental model with inoculation and incubation period is established， the basic reproduction number of the model is obtained by the method of reproducing matrix， and the propagation dynamics of the model is determined by the basic reproduction number. If R0≤1， then the disease free equilibrium is globally asymptotically stable， and if R0gt;1， then the endemic equilibrium is globally asymptotically stable.

Keywords： models of epidemic diseases; basic reproduction number; Lyapunov functional; global stability; equilibrium point

隨著社會經濟的發展和人類生活水平的提高，各種社會問題也相繼出現?？股乇粸E用的問題隨著醫療水平的不斷提高而愈演愈烈，工業廢棄物（廢氣、粉塵、廢水和廢渣等）的排放導致人類生存環境惡化，城市化的迅速發展導致城市人口密度增加，諸多因素導致了微生物的生態平衡發生破壞。一些傳染病依然危害著人類的健康，例如：艾滋病、埃博拉出血熱、新型肝炎等。影響傳染病傳播的因素有很多，包括媒介報道、氣象因素、空氣污染物等[1-2]。

傳染病動力學是對傳染病進行理論性定量研究的一種重要方法。其根據疾病的發生、傳播、發展規律，以及與之有關的社會等因素，建立能反映傳染病動力學特性的數學模型。Kermack 和Mckendrick 描述了一類由病毒或細菌直接感染的傳染病模型，該模型中包括3類人群：易感者（Ｓ）、感染者（Ｉ）、治愈者（Ｒ）。其他關于傳染病動力學的文章可以參考文獻[3-5]

接種對預防傳染病是很重要的。通常情況下，不同疾病的接種過程也是不同的。例如：乙肝疫苗要接種3次，并且每2次之間的時間間隔是固定的[6]。文獻[7]中認為，對易感者和新生兒接種都是可行的，并且接種以后的免疫作用是暫時的。有些傳染病是具有潛伏期的，不同疾病的潛伏時長也有所不同。

1" 模型介紹

本文建立了一類具有接種和潛伏期的傳染病模型，模型中共有5個倉室：易感者（Ｓ）、接種者（Ｖ）、潛伏期患者（Ｅ）、感染者（Ｉ）和治愈者（Ｒ）。其傳播機制圖如圖1所示，根據傳播機制圖得如下模型

式中：？撰表示出生率，？滋表示自然死亡率，p表示易感者的接種率，？茁和？滓？茁（０lt;？滓lt;1）分別表示易感者和接種者的感染率，？酌表示潛伏期患者的轉移率，？孜表示治愈率，d表示因病死亡率。

下面來分析系統的正性和有界性。將式（1）改寫成如下形式

X′=Ｇ（X），（2）

式中：X＝（Ｓ，Ｖ，Ｅ，Ｉ，Ｒ）∈R5，

。（3）

容易得到

。（4）

由文獻[8]中的引例2可知，所有的t≥0時，X（t）∈R■■均成立。

將式（1）中的所有方程相加，得

（Ｓ＋Ｖ＋Ｅ＋Ｉ＋Ｒ）′=？撰-？滋（Ｓ＋Ｖ＋Ｅ＋Ｉ＋Ｒ）-dI≤？撰-？滋（Ｓ＋Ｖ＋Ｅ＋Ｉ＋Ｒ）。（5）

？贅＝（Ｓ＋Ｖ＋Ｅ＋Ｉ＋Ｒ）∈R■■：Ｓ＋Ｖ＋Ｅ＋Ｉ＋Ｒ≤■是式（1）的一個正的不變集。

2" 基本再生數和平衡點的存在性

2.1" 無病平衡點的存在性和基本再生數

令式（1）的右端等于零，則存在一個無病平衡點P0（Ｓ0，Ｖ0，０，０，０），其中

S0＝■，V0＝■

利用基本再生矩陣的方法[9]，求出了模型的基本再生數，即

2.2 地方病平衡點的存在性

設Ｐ*（Ｓ*，Ｖ*，Ｅ*，Ｉ*，Ｒ*）為式（1）的地方病平衡點，滿足下列方程

解方程可得：當Ｒ０gt;1時，式（8）在區間（０，■）上存在唯一的正根。此時，式（1）存在唯一的地方病平衡點Ｐ*（Ｓ*，Ｖ*，Ｅ*，Ｉ*，Ｒ*）。其中

3 平衡點的全局穩定性

3.1 無病平衡點的全局穩定性

定理1：當Ｒ０≤1時，式（1）的無病平衡點Ｐ０是全局漸近穩定的。

證明：無病平衡點Ｐ0（Ｓ0，Ｖ0，０，０，０）中的Ｓ0，Ｖ0滿足下列方程

。（10）

那么式（1）可被改寫為

。（11）

定義Lyapunov函數為

Ｖ１＝（S-S0-S0ln■）+（V-V0-V0ln■）+E+■I 。（12）

V1沿著式（11）的全導數為

（13）

顯然，當Ｒ０≤1時，V1′≤0，式（1）得出的無病平衡點P0是全局漸近穩定的。

3.2" 地方病平衡點的全局穩定性

定理2：當Ｒ０gt;1時，式（1）得出的地方病平衡點P *是全局漸近穩定的。

證明：地方病平衡點P*（Ｓ*，Ｖ*，Ｅ*，Ｉ*，Ｒ*）滿足下列方程

。（14）

記■＝x，■＝y，■＝z，■＝u，那么式（1）可被改寫為

定義Lyapunov函數為

Ｖ２＝Ｓ*（x-1-lnx）+V*（y-1-lny）+E*（z-1-lnz）+■I*（u-1-lnu）。（16）

Ｖ２沿著式（15）的全導數為

由參考文獻[7]可知，V2′≤0。當R0gt;1時，式（1）的地方病平衡點P*是全局漸近穩定的。

4" 結論

本文提出了具有接種與潛伏期的傳染病模型，通過再生矩陣法得到了模型的基本再生數。通過構造Lyapunov函數，證明了當R0≤1時，無病平衡點是全局漸近穩定的，傳染病最終消失；當R0gt;1時，地方病平衡點是全局漸近穩定的，傳染病將持續存在。通過基本再生數可以看出接種對疾病暴發的影響，建議人們盡早接種疫苗，有效預防流行病。

參考文獻：

[1] SOORYANARAIN H， ELANKUMARAN S. Environmental role in inflfluenza virus outbreaks[J]. Annu. Rev.Anim. Biosci， 2015，3（1）：347-373.

[2] ECCLES R. An explanation for the seasonality of acute upper respiratory tract viral infections[J]. Acta oto-laryngologica， 2002，122（2）：183-191.

[3] 王素霞，董玲珍，王曉燕，等.具接種的隨機SIQR流行病模型中平衡點的穩定性分析[J].數學的實踐與認識，2019，49（14）：140-149.

[4] 郭曉君，李大治．一類具有預防接種且帶隔離項的SIQR傳染病模型的穩定性分析[J]．南通大學學報（自然科學版），2008，7（4）：87-93.

[5] 王賀橋，周艷麗，王美娟，等．具有連續預防接種的雙線性傳染率SIQR流行病模型[J].上海理工大學學報，2007，29（2）：113-116.

[6] GABBUTI A， ROMANO L， BLANC P， et al. Long-term immunogenicity of hepatitis B vaccination in a cohort of Italian healthy adolescents[J]. Vaccine， 2007，25：3129-3132.

[7] LI J， YANG Y， ZHOU Y. Global stability of an epidemic model with latent stage andvaccination[J]. Nonlinear Anal. RWA.， 2011，12：2163-2173.

[8] YANG X， CHEN L. Permanence and positive periodic solution for the single-speciesnonautonomous delay diffusive model[J]. Comput. Math. Appl.， 1996，32：109-116.

[9] VAN D D P， WATMOUGH J. Reproduction numbers and sub-threshold endemicequilibria for compartmental models of disease transmission[J]. Math. Biosci.，2002，180： 29-48.