一類具有不同時滯的布魯氏菌病模型的穩(wěn)定性分析

吳嫚　夏米西努爾·阿布都熱合曼

摘要：建立一類同時具有離散時滯和分布時滯的布魯氏菌病模型，研究該模型無病平衡點和地方病平衡點的存在性和局部漸近穩(wěn)定性，并通過構(gòu)造合適的Lyapunov泛函得出各類平衡點的全局漸近穩(wěn)定性.最后，通過數(shù)值模擬對結(jié)果進行解釋，并探究時滯效應(yīng)對疾病傳播過程的影響.

關(guān)鍵詞：布魯氏菌病; 時滯; 穩(wěn)定性; Lyapunov泛函

中圖分類號：O175 文獻標志碼：A 文章編號：1001-8395（2023）05-0616-07

布魯氏菌病是布魯氏菌屬的革蘭氏陰性菌引起的一種人畜共患疾病，牛羊等家畜是布魯氏菌的主要宿主.易感動物可能在與傳染性動物之間的直接接觸中感染疾病，動物患布魯氏菌病通常的癥狀是流產(chǎn)、睪丸炎等，這會對生殖和生育、新生兒存活和產(chǎn)奶量產(chǎn)生不利影響，造成經(jīng)濟損失.布魯氏菌能在流產(chǎn)的胎兒、肉類、乳制品、塵埃、糞便、水、泥漿和土壤中長期存活，易感動物因接觸到環(huán)境中的布魯氏菌而感染也成為主要疾病傳播途徑之一[1].人類常常因為接觸到被感染動物污染的環(huán)境、與感染動物直接接觸、食用未經(jīng)殺菌的牛奶或乳制品、未經(jīng)高溫殺菌的肉類和動物肝臟等感染布魯氏菌病.人類布魯氏菌病通常表現(xiàn)為急性發(fā)熱性疾病，包含虛弱、嗜睡以及發(fā)燒，可以持續(xù)并發(fā)展為慢性致殘性疾病，伴有嚴重并發(fā)癥，例如骨關(guān)節(jié)并發(fā)癥、胃腸道并發(fā)癥和呼氣道并發(fā)癥等.布魯氏菌病在人與人之間的傳播極其罕見.

目前，越來越多的學(xué)者關(guān)注并深入研究布魯氏菌病，更多關(guān)于布魯氏菌病的數(shù)學(xué)模型被建立并用于分析疾病的傳播規(guī)律[2-6].考慮到布魯氏菌能在環(huán)境中長期存活這一因素，眾多學(xué)者建立了具有間接傳播的布魯氏菌病模型[7-9].另外，由于布魯氏菌在宿主體內(nèi)的潛伏期是不可忽略的，許多帶時滯的流行病模型被建立并研究[10].關(guān)于布魯氏菌病的時滯模型主要研究了潛伏期、免疫期以及檢測和剔除感染性動物所需的時間這一類的時滯.Hou 等[11]在布病方面做了很多工作，建立了一個動物布魯氏菌病傳播的通用 SEIB 動力學(xué)模型，考慮了動物布魯氏菌病的一般發(fā)生率和潛伏期離散時滯τ并分析了模型平衡態(tài)的動態(tài)行為; 文獻[12]建立了一個一般的動態(tài)模型及其相應(yīng)的離散時滯模型，發(fā)現(xiàn)時滯對離散時滯模型的平衡點穩(wěn)定性是沒有影響的.此外，文獻[13]建立了一個具有分布時滯的動物布魯氏菌病的一般動力學(xué)模型，通過計算和數(shù)值模擬得出在滿足假設(shè) （H2） 的前提條件下分布時滯不會改變系統(tǒng)的動態(tài)行為這一結(jié)論; 最后數(shù)值模擬結(jié)果表明在不滿足假設(shè) （H2） 的前提下系統(tǒng)將經(jīng)歷周期振蕩.Lolika等[14]建立了一個包含2個離散延遲的布魯氏菌病模型，分析了模型的動力學(xué)性質(zhì).分析和數(shù)值結(jié)果表明在一定的條件下這2種時滯可以使系統(tǒng)不穩(wěn)定，并導(dǎo)致 Hopf 分支.前面提到布魯氏菌能在一定條件的環(huán)境中存活較長的時間，這一性質(zhì)會導(dǎo)致環(huán)境中的布魯氏菌是一段時間內(nèi)的累積量，因此，本文將引入一個分布時滯來表示環(huán)境中布魯氏菌的量.基于上訴考慮，本文將建立一個同時具有離散時滯和分布時滯的布魯氏菌病模型.

1模型建立

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

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

4數(shù)值模擬

5總結(jié)

本文建立了一個具有不同時滯的動物布魯氏菌病SIB動力學(xué)模型.在一般生物學(xué)意義的假設(shè)下，首先證明了系統(tǒng)（1）解的正性和有界性.然后通過構(gòu)造適當(dāng)?shù)?Lyapunov 泛函，分析了系統(tǒng)（1）平衡點的穩(wěn)定性，發(fā)現(xiàn)平衡點的動力學(xué)性質(zhì)依賴于基本再生數(shù)R0：如果R0≤1，無論初始值如何，動物布魯氏菌病最終都會滅絕; 如果R0>1，動物布魯氏菌病的傳播是持續(xù)的，并最終達到地方性穩(wěn)定狀態(tài).這些結(jié)果表明，時滯對系統(tǒng)（1）的動態(tài)行為沒有影響.

參考文獻

[1] MCDERMOTT J J， GRACE D， ZINSSTAG J. Economics of brucellosis impact and control in low-income countries[J]. Revue Scientifique et Technique，2013，32（1）：249-261.

[2] YANG C Y， LOLIKA P O， MUSHAYABASA S， et al. Modeling the spatiotemporal variations in brucellosis transmission[J]. Nonlinear Analysis：Real World Applications，2017，38：49-67.

[3] LI M T， SUN G Q， WU Y F， et al. Transmission dynamics of a multi-group brucellosis model with mixed cross infection in public farm[J]. Applied Mathematics and Computation，2014，237：582-594.

[4] SUN G Q， ZHANG Z K. Global stability for a sheep brucellosis model with immigration[J]. Applied Mathematics and Computation，2014，246：336-345.

[5] KANG Y Y， ABDURAHMAN X. Dynamic behavior of a brucellosis infection model with the effect of infected dogs[J]. IAENG International Journal of Applied Mathematics，2020，50（3）：636-642.

[6] MENG K， ABDURAHMAN X. Study of the brucellosis transmission with multi-stage[J]. Communications in Mathematical Science，2018，20：1-14.

[7] HOU Q， QIN H Y. Global dynamics of a multi-stage brucellosis model with distributed delays and indirect transmission[J]. Mathematical Biosciences and Engineering，2019，16（4）：3111-3129.

[8] AINSEBA B， BENOSMAN C， MAGAL P. A model for ovine brucellosis incorporating direct and indirect transmission[J]. Journal of Biological Dynamics，2010，4（1）：2-11.

[9] 康瑤瑤，夏米西努爾·阿布都熱合曼. 垂直傳染與環(huán)境因素對野牛布魯氏菌傳播的影響研究[J]. 數(shù)學(xué)的實踐與認識，2020，50（18）：139-146.

[10] COOKE K L， VAN-DEN-DRIESSCHE P. Analysis of an SEIRS epidemic model with two delays[J]. Journal of Mathematical Biology，1996，35（2）：240-260.

[11] HOU Q， ZHANG F. Global dynamics of a general brucellosis model discrete delay[J]. Journal of Applied Analysis and Computation， 2016，6（1）：227-241.

[12] HOU Q， WANG T. Global stability and a comparison of SVEIP and delayed SVIP epidemic models with indirect transmission[J]. Communications in Nonlinear Science and Numerical Simulation，2017，43：271-281.

[13] HOU Q， QIN H Y. Global dynamics of a multi-stage brucellosis model with distributed delays and indirect transmission[J]. Mathematical Biosciences and Engineering，2019，16（4）：3111-3129.

[14] LOLIKA P O， STEADY M. Dynamics and stability analysis of a brucellosis model with two discrete delays[J]. Discrete Dynamics in Nature and Society，2018，2018：1-20.

Stability Analysis of a Brucellosis Model with Different Time Delays

WU Man，ABDURAHMAN Xamxinur（College of Mathematics and Systems Science， Xinjiang University， Urumqi 830046， Xinjiang）

Abstract：In this paper， a Brucellosis model with both discrete and distributed time delays is established. The existence and local asymptotic stability of the disease-free equilibrium and endemic equilibrium are studied. The global asymptotic stability of all kinds of equilibrium is obtained by constructing suitable Lyapunov functional. Finally， numerical simulation was used to explain the results， and the influence of the delay effect on the process of disease transmission was explored.

Keywords：Brucellosis; delay; stability; Lyapunov functional

2020 MSC：34D20; 34D23

（編輯 陶志寧）