郝博　張瑞明

摘要：根據(jù)布魯氏菌病傳播的特點(diǎn)，建立了具有病菌傳染項(xiàng)的人畜動(dòng)力學(xué)模型。在該模型中，動(dòng)物種群分為易感項(xiàng)、潛伏期項(xiàng)和感染項(xiàng)，人類則被分為高危人群和低危人群；感染個(gè)體主要是由染病的動(dòng)物或者環(huán)境中的病菌傳染。從理論上求出了系統(tǒng)的基本再生數(shù)，分析了各平衡點(diǎn)局部穩(wěn)定性并且通過極限系統(tǒng)和構(gòu)造適當(dāng)?shù)腖yapunov函數(shù)，證明了無病平衡點(diǎn)和惟一的正平衡點(diǎn)是全局漸近穩(wěn)定的。

關(guān)鍵詞：布魯氏菌病；傳染病模型；基本再生數(shù)；局部穩(wěn)定性；全局穩(wěn)定性

中圖分類號(hào)：O157 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：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）以及相同菌屬所引起的人畜共患傳染病，其主要侵害生殖系統(tǒng)而引起流產(chǎn)、不孕、睪丸炎等[1]。布魯氏菌病通常是由牛布魯氏流產(chǎn)桿菌、小型反芻動(dòng)物布魯氏菌病波狀熱、綿羊布魯氏菌、豬布氏桿菌以及犬種布魯氏菌所引起的[2]。布魯氏菌病不僅會(huì)造成嚴(yán)重的經(jīng)濟(jì)損失，而且還會(huì)感染人體對(duì)人造成的傷害。2012年中國(guó)31個(gè)省市報(bào)道有布魯氏菌病患者，全國(guó)共有患者39 515人，無論人數(shù)還是發(fā)病范圍都達(dá)到近些年的高峰，其中內(nèi)蒙古、黑龍江、山西、河北等我國(guó)最主要的牧區(qū)患病人數(shù)最多，但目前尚未見人與人水平傳染的報(bào)道[3]。

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

1 模型的建立

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

3 小結(jié)與討論

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

參考文獻(xiàn)：

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

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

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

[4] 楊 偉.傳染病動(dòng)力學(xué)的一些數(shù)學(xué)模型及其分析[D].上海：復(fù)旦大學(xué)，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] 聶 靜.中國(guó)奶牛布魯氏菌病動(dòng)力學(xué)建模及防控措施研究[D].廣州：中北大學(xué)，2013.

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

[8] 李明濤，孫桂全，靳 禎.具有階段結(jié)構(gòu)的羊群布魯氏菌病動(dòng)力學(xué)分析[J].中國(guó)科技論文在線精品論文，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] 周義倉(cāng)，靳 禎，秦林軍.常微分方程及其應(yīng)用[M].北京：科學(xué)出版社，2003.