Sakhone Sysavathdy 田守靜 齊龍興

摘 要：基于實際數(shù)據(jù)中不同年齡段病人的比例不一樣這一現(xiàn)象，建立一個老撾湄公血吸蟲病數(shù)學(xué)模型，模型中考慮人宿主和釘螺宿主的感染比例，并將人宿主分為不同的8個年齡段，得出了平衡點的存在性，求出基本再生數(shù)，對模型的平衡點的穩(wěn)定性進行了討論和數(shù)值模擬，最后給出最佳控制策略.

關(guān)鍵詞：老撾湄公血吸蟲病；不同年齡段；穩(wěn)定性；數(shù)值模擬；控制

中圖分類號：O175 文獻標(biāo)志碼：A

0 引言

血吸蟲病是一種慢性寄生蟲病，是由多種因素誘發(fā)的一種疾病，主要受到環(huán)境、寄生蟲、載體和宿主的影響.目前，血吸蟲病仍然是發(fā)病率和死亡率很高的一種疾病.世界衛(wèi)生組織（WHO）認為，血吸蟲病作為僅次于瘧疾的一種疾病，主要存在于第三世界的某些國家中，對世界公共衛(wèi)生與世界經(jīng)濟的發(fā)展有著極大影響[1-3].

1 老撾湄公血吸蟲病的簡介

老撾湄公血吸蟲病只發(fā)現(xiàn)于老撾最南方的占巴塞省，而其他省份沒有血吸蟲病的報道.主要集中在兩個城市：孔市和幕拉巴某市.與孔市相比，幕拉巴某市發(fā)生湄公血吸蟲病的地區(qū)較小.這兩個地區(qū)具有大量的貝類生物，其中一種就是湄公血吸蟲的感染中介.

從老撾公共衛(wèi)生占巴塞省處收集統(tǒng)計的數(shù)據(jù)來看，血吸蟲病感染率從2010年～2015年仍然很高.2010年，老撾血吸蟲病的感染率是3.52%，2011年增長至10.82%，2012年增加到17.78%，2013年下降到2.44%，但2014年又增加到9.02%，感染率最低的是2015年的1.31%.為監(jiān)測評估，需要抽取一個村子，對其進行血吸蟲病檢查，觀察村民是否有感染現(xiàn)象的出現(xiàn).具體監(jiān)測評估是：在每個城市，每隔1年或2年，在2個～4個村子里隨機抽取不同年齡層的300人～500人，對其大便進行檢查，觀察是否有寄生蟲卵出現(xiàn)，以監(jiān)測疾病感染的速度和嚴(yán)重程度.結(jié)果顯示，由于生活方式的不同，不同年齡段的人的感染率不同，如圖1所示.

2 數(shù)學(xué)模型

在中國，日本血吸蟲病及其數(shù)學(xué)模型已經(jīng)有很多相關(guān)的研究[4-8].但針對老撾湄公血吸蟲病還沒有太多的研究，尤其是利用數(shù)學(xué)模型來分析討論湄公血吸蟲病的傳播發(fā)展趨勢及其控制措施就更少；因此，本文根據(jù)湄公血吸蟲病傳播的實際情況以及來自老撾占巴賽省醫(yī)院的資料中顯示的不同年齡段的感染率（圖1）和參考文獻[9-10]，在Barbour模型[11]的基礎(chǔ)上，建立模型 （1）.首先對模型做如下假設(shè)：Pk（t）表示不同年齡段的人宿主中感染的比例，y（t）表示釘螺宿主中感染的比例.

■=aky（1-Pk）-（δH+γH）Pk=f（P1， …， P8， y）；k=1，…，8■=bfH■ekPk（1-y）-（d+θ）y=g（P1，…，P8，y）； （1）

其中（k代表不同年齡段“1：1歲～4歲，2：5歲～9歲，3：10歲～14歲，4：15歲～19歲，5：20歲～29歲，6：30歲～39歲，7：40歲～49歲，8：>49歲”），ak代表不同年齡段的病人傳染率系數(shù)，fH代表人的糞便輸出量，ek代表人的每克糞便的蟲卵數(shù)，δH代表人的死亡率，γH代表人的恢復(fù)率，b代表陽性螺的傳染率系數(shù)，d代表釘螺的自然死亡率，θ代表陽性螺的因病死亡率.根據(jù)基本再生數(shù)的生物意義可得：

R0=bfH■■.

從模型（1）可以得出平衡點的存在性.由aky（1-Pk）-（δH+γH）Pk=0可解得 ：

Pk=■

將Pk代入模型（1）的第二個方程，得：

bfH■ek■（1-y）-（d+θ）y=0，

進而得

y=0和■■=■ .

令：

F（y）=bfH■■-（d+θ），則

F（0）=bfH■■-（d+θ）=（d+θ）（R0-1）.

則有：

F（0）>0，當(dāng)R0>1時，F(xiàn)（0）<0，當(dāng)R0<1時，F(xiàn)（0）=0，當(dāng)R0=1時， 以及 F（1）=-（d+θ）<0，且

F '（y）=bfH■-■=bfH■■<0，

即F（y）關(guān)于y嚴(yán)格單調(diào)遞減.

1）若R0>1，F(xiàn)（0）>0，F(xiàn)（1）<0，F(xiàn) '（y）<0，則在（0，1）中F（y）=0只有唯一的一個y*∈（0，1）.

2）若R0<1，F(xiàn)（0）<0，F(xiàn)（1）<0，F(xiàn) '（y）<0，則在（0，1）中F（y）=0無正解.

所以，當(dāng)R0<1時，無正平衡點， 當(dāng)R0>1時，有唯一的正平衡點（ ），其中

.

由y=0可計算得無病平衡點E0（0，0，0，0，0，0，0，0，0）.下面討論E0（0，0，0，0，0，0，0，0，0）的局部穩(wěn)定性.

記系統(tǒng)（1）為：

■=a1y（1-P1）-（δH+γH）P1 = f1，■=a2y（1-P2）-（δH+γH）P2 = f2， …■=a8y（1-P8）-（δH+γH）P8 = f8，■=bfH（e1P1+…+e8P8）（1-y）-（d+θ）y = g.

可得雅可比矩陣為：

由

可得：

由式（2），解得：

λ1=λ2=…=λ7=-（δH+γH）.

由式（3）得：

λ2+[（δH+γH）+（d+θ）]λ+（δH+γH）（d+θ）-bfH■ekak=0，

其中：

a=1，b=（δH+γH）+（d+θ），c=（δH+γH）（d+θ）-bfH■ekak.

由根與系數(shù)的關(guān)系可得：

λ8+λ9=-[ （δH+γH）+（d+θ）]，λ8λ9=（δH+γH）（d+θ）-bfH■ekak=（δH+γH）（d+θ）（1-R0）.

當(dāng)R0<1時，λ8+λ9<0λ8λ9>0，可解得λ8<0；λ9<0；所以E0是局部漸近穩(wěn)定的.

當(dāng)R0>1時，λ8·λ9<0，說明λ8，λ9中必有一個正的；故E0是不穩(wěn)定的.

下面進行一些數(shù)值模擬，并對其給予分析.

3 數(shù)值模擬及控制

利用Matlab 軟件進行數(shù)值模擬.

根據(jù)參考文獻[10]，選取參數(shù)值如下：a1=1.13×10-4；a2=7.52×10-4；a3=1.18×10-3；a4=1.03×10-3；

a5=6.39×10-4；a6=3.97×10-4；a7=2.47×10-4；a8=1.39×10-4；e1=105；e2=103；e3=195；e4=170；e5=100；e6=95；

e7=75；e8=45；b=5.15×10-6；fH=1 120；δH=3.26×10-4；γH=0.003 8；d=0.02；θ=6×0.02.（時間t的單位是w）

可以計算得：

R0=bfH■■=6.402 9>1.

說明該地區(qū)湄公血吸蟲病會持續(xù)傳播下去.圖2～圖10也進一步驗證了這一理論結(jié)果.

由于老撾的國情特色，滅螺措施無法進行.為了更好的預(yù)防湄公血吸蟲病，公共衛(wèi)生部門開設(shè)了一系列的健康教育課程，以提高人們的衛(wèi)生健康意識，規(guī)范人們的日常行為，進而減少疾病發(fā)生的概率.所以在選取控制措施時，只能通過調(diào)控傳染率，即參數(shù)a1，a2，…，a8來研究并比較哪種措施最好.下面分別將8個傳播病率分別減小一半來進行比較，圖形如下（圖11～圖19）：

通過圖11～圖19可以看出，減小不同的傳染率（a1，a1，…，a8）的效果不同. 從P1～P8的圖形可以看出，除了減小相對應(yīng)的各自傳染率系數(shù)外，控制效果最顯著的是減小a3，最不明顯的是減小a8.此結(jié)論說明：在老撾，采取一定的措施使得10歲～14歲年齡段的病人的傳染率減小，對控制老撾血吸蟲病的傳播是最有效的一種策略.

4 結(jié)論

本文主要工作是建立老撾血吸蟲病數(shù)學(xué)模型，得出無病平衡點和地方病平衡點的存在性，并對平衡點的穩(wěn)定性進行了討論和數(shù)值模擬，最后給出最佳控制策略.

在老撾一些疫區(qū)，當(dāng)反復(fù)的臨床療法和環(huán)境改善相結(jié)合時，血吸蟲病感染和傳播控制有所減緩.然而，老撾湄公血吸蟲病因為中間宿主的生態(tài)特性，現(xiàn)有水平的很多干預(yù)似乎效果都不是很突出.首先，由于湄公河的水位很高，氯硝柳胺對釘螺沒有強烈的影響，滅螺劑的應(yīng)用在老撾以失敗告終.其次，由于沒有足夠的資金，當(dāng)?shù)氐纳鐓^(qū)并沒有建設(shè)完善的公共廁所和洗浴措施，通過管理人體排泄物來減少傳播的方法無法維持下去；除了外部組織的一小部分的干預(yù)，目前政府還沒有計劃用于補貼或其他形式的資金支持環(huán)境改善；因而，改善環(huán)境的舉措也是不太可行的.另外，動物宿主還有待進一步確認.目前老撾湄公血吸蟲病防治的目標(biāo)仍然是控制發(fā)病率和控制傳播途徑.

參考文獻

[1] ELBAZ T， ESMAT G. Hepatic and intestinal schistosomiasis： review [J]. Journal of Advanced Research，2013，4（5）：445-452.

[2] 中國疾病預(yù)防控制中心. 血吸蟲病防治知識問答：[EB/OL].北京：中國疾病預(yù)防控制中心，2009[2017-01-17]. http：//www.chinacdc.cn/n272442/n272530/n3226631/30541.html.

[3] 中國疾病預(yù)防控制中心. 血吸蟲病：[EB/OL].北京：中國疾病預(yù)防控制中心，2005[2017-01-17]. http：//www.chinacdc.cn/jkzt/crb/xxcb/.

[4] 薛靖波，夏尚，張瑕，等.血吸蟲病傳播控制階段湖沼型流行區(qū)傳播風(fēng)險時空分布特征研究[J]. 中國血吸蟲病防治雜志， 2016，28（6）：624-629.

[5] 齊龍興，薛夢，甘莉娟，等. 血吸蟲病在兩類易感群體中傳播的穩(wěn)定性分析[J]. 安徽大學(xué)學(xué)報（自然科學(xué)版），2015，39（6）：9-14.

[6] 陳艷艷，蔡順祥，肖瑛，等.應(yīng)用時間序列模型預(yù)測湖北省血吸蟲病流行趨勢[J].中國血吸蟲病防治雜志，2014，26（6）：613-617.

[7] 惠靜，陳蘭蓀.一個脈沖捕食-食餌系統(tǒng)的滅絕與持續(xù)生存[J]. 廣西工學(xué)院學(xué)報，2004，15（1）：1-4.

[8] 莫玉忠.一類具有時滯的控制系統(tǒng)的絕對穩(wěn)定性[J]. 廣西工學(xué)院學(xué)報，2005，16（4）：85-87.

[9] URBANI C， SINOUN M， SOCHEAT D. Epidemiology and control of Mekongii schistosomiasis[J]. Acta Tropica，2002，82（2）：157

-168.

[10] HISAKANE N， KIRINIKI M， CHIGUSA Y. The evaluation of control measures against Schistosoma mekongi in ambodia by a mathematical model[J]. Parasitogy International，2008（57）：379-385.

[11] 吳開琛. 血吸蟲病數(shù)學(xué)模型和傳播動力學(xué)及其應(yīng)用[J].中國熱帶醫(yī)學(xué)，2005，5（4）：837-844.