王 茉,劉俊利

(西安工程大學 理學院,陜西 西安 710048)

0 引 言

傳染病一直是人類致力于攻克的難題,為了預防和控制疾病的傳播,學者們常常建立微分方程數(shù)學模型研究分析傳染病特性[1-2]。麻疹是一種由副黏病毒感染引起的急性呼吸道疾病,其傳染性極強,在人口密集區(qū)域且未普及接種疫苗地區(qū)易流行傳播,約2～3年形成一次大流行。其中兒童期為高發(fā)年齡段,且大多數(shù)是5歲以下兒童,發(fā)病率可高達95%以上,是全球幼兒死亡的主要原因之一。臨床表現(xiàn)為發(fā)熱、上呼吸道炎癥、眼結膜炎以及皮膚上現(xiàn)斑丘疹等,常見的并發(fā)癥有中耳炎、肺炎、麻疹腦炎,嚴重者可危及生命。潛伏期最長可達21天,并且麻疹發(fā)病具有較為明顯的季節(jié)性規(guī)律,其中春季和冬季為高發(fā)季節(jié)[3-8]。目前接種常規(guī)麻疹疫苗是預防麻疹疾病的關鍵公共衛(wèi)生戰(zhàn)略,但在實際生活中,往往存在漏種、不及時接種以及自身免疫等因素使接種效果受到影響,因此對麻疹病情的傳播及治療還需要做進一步的深入研究。

目前,國內外諸多學者對于麻疹疾病的流行病學特點都已經進行了較為全面的分析。其中,考慮到接種一劑疫苗的時效性問題,2015年佘連兵[9]建立了一類具有二次接種的麻疹傳染病模型,并對模型的全局穩(wěn)定性進行了分析。由于麻疹患病者的治療情況較為復雜,2017年Garba等[10]研究了關于疫苗接種和麻疹治療相結合的確定性模型,定性評估了接種和治療相結合的方式對于麻疹染病人群的綜合影響。隨著二次接種的普及率越來越高,Li等[11]針對我國麻疹疫苗接種的二劑常規(guī)劑量情況進行了建模和分析,并通過數(shù)值模擬驗證了二劑疫苗接種率對麻疹疾病控制的影響。Stephane等[12]建立并研究了一類通過補充常規(guī)免疫的偏微分方程,得出了通過提高麻疹疫苗接種覆蓋率可控制麻疹傳播的結論。在文獻[12]的研究基礎上,本文考慮了麻疹分人群接種的傳播影響,建立并分析一類帶有分層接種的常微分方程麻疹模型,證明了模型平衡點的穩(wěn)定性以及疾病的一致持久性。

1 模型建立

將總人口N(t)分為6類:S(t)、I(t)、R(t)、V1(t)、V2(t)、V3(t)分別表示t時刻人群中的易感者、染病者、恢復者、接種易感者、接種染病者和接種恢復者,模型如下:

(1)

式中:A表示人口輸入率;β表示感染率;r表示恢復率;μ表示自然死亡率;k1、k2分別表示第1次、第2次接種疫苗覆蓋率;τ1、τ2分別表示第1次、第2次接種疫苗的有效性;θ1表示由于接種原因,相對于染病者,接種染病者對易感者的相對傳染性;θ2表示相對于易感者的傳染,染病者和接種染病者對接種易感者的相對傳染性;并且0≤τ1≤1、0≤τ2≤1、0≤θ1≤1、0≤θ2≤1。

定義集合Ω如下:

假定模型(1)中所有參數(shù)都是非負的,且A>0,μ>0,根據(jù)文獻[13]中定理2.1的證明,可得系統(tǒng)的解是非負的。將模型(1)中的6個方程相加得總人口N(t)=S(t)+I(t)+R(t)+V1(t)+V2(t)+V3(t)滿足方程

因此Ω為模型(1)的正向不變集。

2 平衡點和基本再生數(shù)

模型(1)總存在唯一的無病平衡點

根據(jù)文獻[14-15],有

(2)

計算得到模型(1)的基本再生數(shù)為

R0=[Aβ(μ+τ2k2)(μ+r+k1θ1)+Aβk1θ1θ2(1-τ1)(k1+μ+r)]/

[(k1+μ+r)(μ+r)(k1+μ)(μ+τ2k2)]

3 平衡點的局部穩(wěn)定性

定理1如果R0<1,則模型(1)的無病平衡點E0是局部漸近穩(wěn)定的;如果R0>1,則E0不穩(wěn)定。

證明模型(1)在平衡點E0處的雅可比矩陣為

J(E0)=

得

則E0處的特征方程為

(λ+μ)(λ+k1+μ)2(λ+μ+τ2k2)·

[(λ-x1+k1+μ+r)(λ-θ1x2+μ+r)-

θ1x1(x2+k1)]=0

因此,該特征方程存在4個負特征根:λ1=-μ、λ2=-μ-τ2k2、λ3=λ4=-k1-μ。其余特征根滿足方程

λ2+aλ+b=0

其中

由 Hurwitz 判據(jù)[16]可知, 當R0<1時, 有a>0,b>0, 從而方程λ2+aλ+b=0 的根均具有負實部, 此時無病平衡點E0是局部漸近穩(wěn)定的。當R0>1 時, 得b<0, 即方程λ2+aλ+b=0 有一個正實根, 此時無病平衡點E0是不穩(wěn)定的, 從而定理1結論成立。

4 疾病的持久性

用s(A)表示n×n矩陣A的穩(wěn)定性模,定義為

s(A)max{Rez:z定義為A的一個特征值}

由式(2)可知,顯然F-V是不可約的,并且具有非負的非對角元,則s(F-V)是(F-V)的一個簡單特征值,且對應正的特征向量[17]。

引理1[13]下列等價條件成立:

R0>1?s(F-V)>0,R0<1?s(F-V)<0

定理2如果R0>1,則模型(1)至少存在一個正平衡點,且?ε>0,使得模型(1)具有非負初值條件的任意解(S(t)，I(t)，R(t)，V1(t)，V2(t)，V3(t))滿足

證明由模型(1)得

則

則?η>0,?t0>0,當t>t0時,有

S(t)>m1-η

(3)

當t>t0時,由式(3)有

則

當t>t0時,由式(3)有

rV2-μV3≥τ1k1S-μV3

則

定義

X0={(S,I,R,V1,V2,V3)∈X,

I+V2>0,?t≥0}

?X0=XX0

易證模型(1)關于X和X0是正不變的。

令Φ(t)=(S(t),I(t),R(t),V1(t),V2(t),V3(t)),定義

M?={Φ(0);Φ(t)∈?X0,t≥0}

下證

M?={(S,0,R,V1,0,V3):S≥0,

R≥0,V1≥0,V3≥0}

(4)

為了證明式(4)成立,假設Φ(0)∈M?,只需證明?t≥0,有I=0,V2=0即可。若不成立,則?t1≥0,使得I(t1)+V2(t1)>0,因此Φ(t1)∈X0這與Φ(0)∈M?相矛盾。 用ω(Φ(0))表示模型(1)從Φ(0)∈X出發(fā)的解的ω-極限集。令

P=U(ω(Φ(0)):Φ(0)∈M?}

在M?={(S,0,R,V1,0,V3):S≥0,R≥0,V1≥0,V3≥0}上有

(5)

其中Φ(t)是初值屬于X0的模型(1)的任意解,只需證明

Ws(E0)∩X0=?

(6)

其中Ws(E0)是E0的穩(wěn)定流形。假設式(6)不成立,則在X0中存在一個解Φ(t),滿足

(7)

則存在充分小的δ>0,當T1>0時,有

(8)

當t≥T1時,根據(jù)模型(1)的第2個和第5個方程,由式(8)得

即

其中

5 數(shù)值模擬

通過數(shù)值模擬來驗證模型(1)理論結果的正確性。從文獻[10]中得到參數(shù)τ1=0.85,τ2=0.98,文獻[19]中得到參數(shù)A=1 340 000,k1=0.629 1,k2=0.8,r=1/0.633 3。參數(shù)以月為單位,假設人的平均壽命是840個月,因此死亡率u=1/840≈0.001 2?？紤]參數(shù)θ1=0.01,θ2=0.05，可得以下結論。

當β=0.000 008 5時,求得R0=8.212 5>1,圖1分別表示染病者和接種染病者與時間的關系,顯示了麻疹疾病是持續(xù)存在的,與定理2結論一致。

當β=0.000 000 85時,無病平衡點為E0=(2.1×106,0,0,2.6×105,0,1.1×109),R0=0.821 3<1。圖2分別表示染病者與接種染病者與時間的關系,表明此時無病平衡點E0是局部漸近穩(wěn)定的,麻疹疾病絕滅,與定理1結論一致。

圖3為部分相關系數(shù)(PRCC)的結果R0對每個參數(shù)的依賴性。從圖3可知,β和A對R0有較大的正的影響,r和k1對R0有較大的負的影響。相比較而言,μ,θ1和τ1對R0的影響較小,即不敏感。另外,麻疹疾病的感染率β對R0的影響最大,其次是恢復率r,然后是第一次接種覆蓋率k1和人口輸入率A,這些因素在很大程度上影響著疾病的存在和持續(xù)傳播,因此,為了減少麻疹染病者人數(shù)的增加并達到控制疾病流行的目的,可以采取的有效措施有:①盡量減少麻疹染病者和普通人群以及從未接種過麻疹疫苗的易感者的接觸,從而達到降低感染率的目的;②在麻疹患者染病的整個潛伏期過程中,需盡量縮短治愈周期,控制治療時間在潛伏期范圍內,以期望達到高效的治愈率,從而避免了麻疹疾病再次傳播;③全面擴大新生兒首次接種覆蓋率,第一次接種的覆蓋率對基本再生數(shù)的影響要大于第二次接種覆蓋率,因此新生兒出生后的第一劑麻疹疫苗覆蓋率需盡可能達到100%。

6 結 語

本文在文獻[12]的基礎上考慮了分人群接種對麻疹疾病傳播動力學的影響,通過計算得到了模型的無病平衡點和基本再生數(shù),證明了無病平衡點的局部漸近穩(wěn)定性和疾病的一致持久性,并進行數(shù)值模擬驗證。結果表明：感染率、接種覆蓋率以及恢復率對基本再生數(shù)的影響較大,對閾值的敏感性較強。由此總結得到控制麻疹流行的3個有效措施,從而達到控制麻疹疾病流行的目的。