動(dòng)力學(xué)模型在預(yù)測深圳市寶安區(qū)結(jié)核病的應(yīng)用

張磊 劉艷紅

（深圳市寶安區(qū)慢性病防治院 廣東 深圳 518133）

動(dòng)力學(xué)模型在預(yù)測深圳市寶安區(qū)結(jié)核病的應(yīng)用

張磊 劉艷紅

（深圳市寶安區(qū)慢性病防治院 廣東 深圳 518133）

目的：預(yù)測深圳市寶安區(qū)結(jié)核病的發(fā)病趨勢(shì)。方法：采用動(dòng)力學(xué)模型對(duì)寶安區(qū)2010－2013年的肺結(jié)核數(shù)據(jù)建立動(dòng)力微分學(xué)方程組，并對(duì)擬合結(jié)果進(jìn)行評(píng)價(jià)，對(duì)2014－2016年新發(fā)涂陽肺結(jié)核數(shù)據(jù)分別進(jìn)行了預(yù)測。結(jié)果：預(yù)測2014－2016年結(jié)核病的發(fā)病人數(shù)分比為1355，1259和1170。結(jié)論：寶安區(qū)的結(jié)核病防治工作卓有成效，結(jié)核病已進(jìn)入下降通道

肺結(jié)核病；預(yù)測；動(dòng)力學(xué)模型

結(jié)核病是嚴(yán)重危害人類健康的慢性傳染病。自上世紀(jì)八十年代后期以來，由于人口流動(dòng)加劇、HIV/AIDS流行和耐藥性結(jié)核病增多等因素的影響，結(jié)核病在全球重新蔓延，己成為全球重大的公共衛(wèi)生問題和社會(huì)問題。本研究利用SEIR倉室的理論和方法，充分考慮流動(dòng)人口、耐藥結(jié)核病和合并HIV感染等因素，建立能反映結(jié)核病動(dòng)力學(xué)特征的數(shù)學(xué)模型。通過對(duì)模型動(dòng)力學(xué)特征的定性、定量分析和數(shù)值模擬，預(yù)測未來結(jié)核病的發(fā)病情況，顯示疾病的發(fā)展過程，揭示流行規(guī)律，為人們防制決策提供理論基礎(chǔ)和數(shù)據(jù)支持。

1.材料和方法

1.1 資料來源

由于結(jié)核病屬于國家法定乙類傳染病，寶安區(qū)實(shí)行結(jié)核病歸口管治，因此從寶安區(qū)慢性病防治院收集了2010－2014年肺結(jié)核病發(fā)病人數(shù)，其中發(fā)病人數(shù)包括了涂陽肺結(jié)核和涂陰肺結(jié)核患者。由于2011年12月成立龍華新區(qū)，因此本文相關(guān)數(shù)據(jù)均剔除上述龍華新區(qū)相關(guān)數(shù)據(jù)。

1.2 分析方法

1.2.1 模型假設(shè)

（1）人群對(duì)結(jié)核病普遍易感，患結(jié)核病的幾率一致；

（2）結(jié)核病治療后的傳染性降低，和普通人群一致；

（3）若合并HIV感染則轉(zhuǎn)為無法治愈者；

（4）不考慮人口出生對(duì)結(jié)核病的影響；

（5）本地區(qū)的人口總量基本保持不變；

（6）兒童接種卡介苗對(duì)結(jié)核病的發(fā)病無影響。

1.2.2 參數(shù)說明

S：易感人群，即普通人群；

I1：去向不明或治療失敗患者；

I2：發(fā)病人群中治愈者或成功完成療程患者；

R：耐藥者；

N：研究總?cè)巳海ㄒ赘腥巳?，發(fā)病人群；

ε：未治療或去向不明者占患病者比重；

C：每年每名患者平均C次能引起感染的接觸次數(shù)；

r：易感人群受感染的傳染率

1.2.3 建立結(jié)核病在人群中的模型

1.2.4 研究內(nèi)容

對(duì)上述初步模型建立結(jié)核病動(dòng)力微分學(xué)方程組；確定關(guān)鍵參數(shù)，對(duì)研究地區(qū)的動(dòng)力學(xué)方程進(jìn)行驗(yàn)證、優(yōu)化；根據(jù)模型，預(yù)測研究地區(qū)的中短期（2年左右）的趨勢(shì)。

1.3 統(tǒng)計(jì)軟件

采用matlab對(duì)微分學(xué)方程組編程求解方程組，模擬結(jié)核病在各個(gè)人群中的轉(zhuǎn)歸。

2.結(jié)果

2.1 建立的結(jié)核病動(dòng)力微分學(xué)方程組

2.2 模型初值的確定

2010 年深圳市寶安區(qū)人口265.24萬，發(fā)現(xiàn)結(jié)核病人1940例，成功治療959例，成功治療率50.31，治療失敗的86例，耐藥10例。

2.3 模型擬合結(jié)果

經(jīng)過計(jì)算擬合，發(fā)現(xiàn)病人發(fā)現(xiàn)數(shù)有下降趨勢(shì)，成功治療數(shù)上升，耐藥數(shù)也有下降，擬合結(jié)果良好，具體詳情見表1。

表1 2010年－2013年實(shí)際值與模擬值統(tǒng)計(jì)表

2.4 模型預(yù)測

利用建立的模型，預(yù)測2014－2016寶安區(qū)結(jié)核病的相關(guān)數(shù)據(jù)，具體見表2。

表2 2014年－2016年預(yù)測結(jié)核病模型相關(guān)結(jié)果統(tǒng)計(jì)表

3.討論

目前對(duì)結(jié)核病的流行病學(xué)研究方法主要是描述性研究、分析性研究、實(shí)驗(yàn)性研究，而動(dòng)力學(xué)研究卻相對(duì)較少。宋妮［1］在結(jié)核病傳播的數(shù)學(xué)模型的動(dòng)力學(xué)基礎(chǔ)上，利用Hurwitz判據(jù)和復(fù)合矩陣?yán)碚?，討論了具有常?shù)移入的結(jié)核病傳播的數(shù)學(xué)模型的動(dòng)力學(xué)性質(zhì)，研究了具有常數(shù)移入的結(jié)核病模型的地方病平衡點(diǎn)的局部和全局漸進(jìn)穩(wěn)定性，得到了具有常數(shù)移入的結(jié)核病模型的地方病平衡點(diǎn)是局部和全局漸進(jìn)穩(wěn)定的充分條件。賀曉新［2］根據(jù)中國結(jié)核病流行病學(xué)規(guī)律，針對(duì)Aruma模型的局限性，建立聯(lián)立遞推方程模型。利用所建模型對(duì)北京市1979－－－2000年的疫情發(fā)展過程進(jìn)行模擬，并與實(shí)際進(jìn)行比較以評(píng)價(jià)模型的擬合優(yōu)度。朱帥［3］利用Jacobian矩陣性態(tài)和索引文獻(xiàn)的結(jié)論.證明了無病平衡點(diǎn)的全局漸近穩(wěn)定性和地方病平衡點(diǎn)的全局漸近穩(wěn)定性.Sungim Whanga［4］運(yùn)用流行病學(xué)動(dòng)力學(xué)理論，開發(fā)了結(jié)核病傳播時(shí)態(tài)序列的SEIR模型，利用該模型，他們提出了優(yōu)化的結(jié)核病控制策略，同時(shí)考慮了結(jié)核病發(fā)病率的降低與成本控制。針對(duì)感染結(jié)核桿菌人數(shù)的增加、耐多藥結(jié)核病和在吸毒人員中HIV多種嚴(yán)重問題，Lebcir，R M［5］采用了系統(tǒng)動(dòng)力學(xué)模型來模擬、評(píng)估結(jié)核病的流行，提出了多種建設(shè)性的意見。

本次利用2010年的初始數(shù)據(jù)作為模型初值，利用matlab軟件，建立了微分學(xué)方程組，模擬結(jié)核病在人群中的轉(zhuǎn)歸，經(jīng)過2011年－2013年的模擬，絕大多數(shù)指標(biāo)擬合良好，顯示模型在中短期的預(yù)測方面有較好的擬合效果。2014－2016年的預(yù)測結(jié)果發(fā)現(xiàn)結(jié)核病已進(jìn)入下降通道，發(fā)病人數(shù)和成功治療人數(shù)均向有利的方向發(fā)展，提示寶安區(qū)的結(jié)防工作卓有成效。

但是和任何一種時(shí)間序列的預(yù)測模型一樣，模型建立后都需要及時(shí)補(bǔ)充和更新數(shù)據(jù)，使模型得以不斷改進(jìn)，從而使預(yù)測更加符合實(shí)際。

［1］宋妮，薛亞套，何志本，等 .具有常數(shù)移人的結(jié)核病模型穩(wěn)定性分析［J］.中北大學(xué)學(xué)報(bào)（自然科學(xué)版），2009，（3）：197－201.

［2］賀曉新、金水高、張立興，等 .結(jié)核病疫情發(fā)展預(yù)測模型的建立及評(píng)價(jià)［J］.中國防癆雜志，2004，26（1）：10－15.

［3］朱帥，馬純，等.一類具有常數(shù)輸人的SLIT的TB模型的分析［J］.太原師范學(xué)院學(xué)報(bào)（自然科學(xué)版），2007，6（3）：114－115，119.

［4］Sungim Whanga，Sunhwa Choib，Eunok Jungb.A dynamic model for tuberculosis transmission and optimal treatment strategies in South Korea［J］.Journal of Theoretical Biology，2011，（1）：120－131.

［5］Lebcir，R M1；Atun，R A2；Coker，R J.System Dynamic simulation of treatment policies to address colliding epidemics of tuberculosis，drug resistant tuberculosis and injecting drug users driven HIV in Russia［J］.Journal of the Operational Research Society，2010，（8）：1238－1248.

Dynamical Model in Predicting Tuberculosis in Bao＇an District of Shenzhen City

Zhang Lei，Liu Yanhong Bao＇an Chronic Diseases Prevent and Cure Hospital，Shenzhen，Guangdong 518133

Objective：To predict the trends of the incidence of tuberculosis in Bao＇an District of Shenzhen City.Methods：The dynamical model was used to build momentum in differential equations on the basis of the tuberculosis data from 2010to 2013in Bao＇an district，and the results of the fitting was evaluated，the data of new smear positive pulmonary tuberculosis from 2014to 2016was predicted.Results：The predicted number of tuberculosis from 2014to 2016would be 1355，1259and 1170respectively.Conclusion：TB control in Bao＇an District was effective，and tuberculosis has been declined.

TB；predict；dynamical model

R181

B

1009－6019（2015）08－0025－02

張磊，男，江蘇人，1981年，碩士研究生，主治醫(yī)師，從事傳染性疾病控制工作。

深圳市科技計(jì)劃項(xiàng)目（醫(yī)療衛(wèi)生類）：201203266