布魯氏菌病人—畜動力學建模及穩定性分析

郝博　張瑞明

摘要：根據布魯氏菌病傳播的特點，建立了具有病菌傳染項的人畜動力學模型。在該模型中，動物種群分為易感項、潛伏期項和感染項，人類則被分為高危人群和低危人群；感染個體主要是由染病的動物或者環境中的病菌傳染。從理論上求出了系統的基本再生數，分析了各平衡點局部穩定性并且通過極限系統和構造適當的Lyapunov函數，證明了無病平衡點和惟一的正平衡點是全局漸近穩定的。

關鍵詞：布魯氏菌病；傳染病模型；基本再生數；局部穩定性；全局穩定性

中圖分類號：O157 文獻標識碼：A 文章編號：0439-8114（2015）24-6324-04

DOI：10.14088/j.cnki.issn0439-8114.2015.24.059

Abstract： According to the characteristics of Brucellosis，a dynamic model with Brucella in environment was established. In this model，the animal population was divided into susceptible，latent，infected. Human was divided into high-risk group and low-risk group；Susceptible individuals can contract the disease in two ways： infected animals and Brucella in environment. The basic reproductive number was calculated theoretically， the local stability of each equilibrium point was analyzed， then through the limit system and constructing appropriate Lyapunov functions to prove that the disease-free equilibrium and the unique positive equilibrium was globally asymptotically stable.

Key words：Brucellosis；epidemiological model；basic reproduction number；local stability；global stability

布魯氏菌病，是由布魯氏菌（Brucella）以及相同菌屬所引起的人畜共患傳染病，其主要侵害生殖系統而引起流產、不孕、睪丸炎等[1]。布魯氏菌病通常是由牛布魯氏流產桿菌、小型反芻動物布魯氏菌病波狀熱、綿羊布魯氏菌、豬布氏桿菌以及犬種布魯氏菌所引起的[2]。布魯氏菌病不僅會造成嚴重的經濟損失，而且還會感染人體對人造成的傷害。2012年中國31個省市報道有布魯氏菌病患者，全國共有患者39 515人，無論人數還是發病范圍都達到近些年的高峰，其中內蒙古、黑龍江、山西、河北等我國最主要的牧區患病人數最多，但目前尚未見人與人水平傳染的報道[3]。

傳染病數學模型的理論研究涉及常微分方程組、時滯微分方程組、偏微分方程組、一隨機微分方程等方面[4]。A？觙nseba等[5]建立了具有易感者、染病者、環境中布魯氏菌（SIC）的羊群布魯氏菌病動力學模型，模型考慮了直接接觸傳染和環境中布魯氏菌的間接傳染。聶靜等[6，7]根據奶牛布魯氏菌病的傳播特征及其規律建立了具有潛伏期和間接傳染的SEIV模型，從理論上求出了系統的基本再生數，分析了各平衡點的穩定性。李明濤等[8]根據布魯氏菌病具有年齡特征，建立了具有階段結構的羊群布魯氏菌病動力學模型，證明了平衡點的全局漸近穩定性。

1 模型的建立

羊、牛、豬是中國最主要的家畜，需求量龐大，這3種家畜的主要養殖地又是中國布魯氏菌病的主要疫區，存在大量的相互調入。布魯氏菌病的傳播速率很快，而人和家畜對其都易感，因此對于患病動物必須及時撲殺。布魯氏菌對外界的理化因素具有一定的抵抗力，可以生存較長的時間。此外，布魯氏菌可侵入呼吸、消化、生殖系統黏膜以及損傷甚至未損傷完整皮膚等[1，9]。人患病與職業有十分密切的關系，并且農村發病率遠遠高于城市，所以將人群分為高危人群和低危人群。不同動物以及不同規模農場存在差異，動物之間和動物與人之間的接觸存在差異，因此它們之間的接觸與總的群體數量形成一定的比例關系。

3 小結與討論

本研究主要根據布魯氏菌病本身的特點及其傳播特性，建立了具有外界輸入，潛伏期和環境中病菌傳染項的人畜動力學模型，通過基本再生數R0的表達形式可以看出，動物和病菌是布魯氏菌病的傳染源貢獻度，定理1證明了R0<1時無病平衡點是全局穩定的，隨著時間的推移，所以種群的患病者都趨于0，也就是疾病最終消亡。定理3證明了R0>1時正平衡點是全局穩定的，患病者數量最終趨于穩定的正值，也就是形成地方病。系統的無病平衡點和正平衡點都是全局漸近穩定的，環境中的布魯氏菌也是布魯氏菌病傳播的一個重要因素，在實際中應該采取控制引入量，及時有效地殺菌和患病動物捕殺相結合的方式才能使疾病得到有效控制。

參考文獻：

[1] 金寧一，胡仲明，馮書章，等.新編人獸共患病學[M].北京：科學出版社，2007.

[2] 文學忠，于瑞華，姜秋杰.布魯氏菌病近況[J].吉林畜牧獸醫，2007（5）：20-23.

[3] 毛景東，王景龍，楊艷玲.布魯氏菌病的研究進展[J].中國畜牧獸醫，2011，38（1）：222-226.

[4] 楊 偉.傳染病動力學的一些數學模型及其分析[D].上海：復旦大學，2010.

[5] A？魵NSEBA B， BENOSMAN C，MAGAL P. A model for ovine brucellosis incorporating direct and indirect transmission[J]. J Biol Dyn，2010，4（1）：2-11.

[6] 聶 靜.中國奶牛布魯氏菌病動力學建模及防控措施研究[D].廣州：中北大學，2013.

[7] 聶 靜，萬 輝，張 娟.奶牛布魯氏菌病的動力學分析[J].中北大學學報（自然科學版），2013，34（2）：93-97.

[8] 李明濤，孫桂全，靳 禎.具有階段結構的羊群布魯氏菌病動力學分析[J].中國科技論文在線精品論文，2014，7（1）：52-57.

[9] MATOPE G，BHEBHE E.Risk factors for Brucella spp. infection in small holder household herds[J]. Epidemiol Infect，2011，139：157-164.

[10] VAN DEN DRIESSCHE P，WATMOUGH J.Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission[J]. Math Biosci，2002，180：29-48.

[11] 周義倉，靳 禎，秦林軍.常微分方程及其應用[M].北京：科學出版社，2003.