揚(yáng)州市禽流感風(fēng)險(xiǎn)評(píng)估模型的研究進(jìn)展探析

郭志榮 徐步

摘 要 揚(yáng)州水網(wǎng)分布密集而廣泛，養(yǎng)殖以水禽鴨、鵝為主，同時(shí)揚(yáng)州處在候鳥東部遷徙線上，境內(nèi)的高郵湖、邵伯湖是候鳥的重要棲息地，也是候鳥和家禽的交匯點(diǎn)。因此，禽流感入侵揚(yáng)州市的風(fēng)險(xiǎn)很大?；诖?，針對(duì)揚(yáng)州市特定的生態(tài)系統(tǒng)和養(yǎng)殖環(huán)境，通過文獻(xiàn)分析和專家建議，確定了揚(yáng)州市禽流感發(fā)生的風(fēng)險(xiǎn)因子，采用層次分析法建立揚(yáng)州市高致病性禽流感傳播風(fēng)險(xiǎn)評(píng)估模型，以期實(shí)現(xiàn)揚(yáng)州市禽流感疫情發(fā)生風(fēng)險(xiǎn)的定量評(píng)估。

關(guān)鍵詞 禽流感;風(fēng)險(xiǎn)評(píng)估模型;層次分析法;揚(yáng)州

中圖分類號(hào)：S858.3 文獻(xiàn)標(biāo)識(shí)碼：B DOI：10.19415/j.cnki.1673-890x.2019.20.076

自2013年3月在我國(guó)首次發(fā)現(xiàn)人感染H7N9禽流感病例以來，禽流感成為國(guó)人聞之色變的疫情，給國(guó)人造成的心理恐慌不亞于2003年的SARS疫情。進(jìn)入2017年，禽流感在我國(guó)高發(fā)勢(shì)頭更加明顯，截止到2017年4月24日，已有22個(gè)省份報(bào)告H7N9病例511例，死亡188人，其中江蘇省尤為突出，感染發(fā)病人數(shù)和死亡人數(shù)均居全國(guó)第一，感染發(fā)病人數(shù)達(dá)124人，其中35人死亡[1]。2018年，我國(guó)內(nèi)地H5+H7流感疫情共有12起，其中上半年多是北方蛋雞場(chǎng)發(fā)生H7N9疫情，而下半年多是南方家禽場(chǎng)發(fā)生H5N6疫情，全國(guó)共撲殺家禽20多萬羽[2]。禽流感疫情導(dǎo)致人們心理恐慌，嚴(yán)重影響了社會(huì)生活的穩(wěn)定。

近幾年，揚(yáng)州市發(fā)生了兩次較大的禽流感疫情，發(fā)生地都在江都區(qū)的規(guī)?；仪蒺B(yǎng)殖場(chǎng)。2015年8月3日，江都區(qū)一養(yǎng)殖場(chǎng)出現(xiàn)H5N2亞型高致病性禽流感疫情，該養(yǎng)殖場(chǎng)飼養(yǎng)的鵝有23395只發(fā)病，其中死亡3106只。2018年11月22日，江都區(qū)一家禽養(yǎng)殖場(chǎng)發(fā)生H5N6亞型高致病性禽流感疫情，該養(yǎng)殖場(chǎng)存欄家禽11810羽，發(fā)病1200羽，死亡320羽。疫情發(fā)生后，江都區(qū)疾控防治中心嚴(yán)格按照疫情防治規(guī)范要求，對(duì)養(yǎng)殖場(chǎng)周邊家禽進(jìn)行了撲殺和無害化處理。2015年8月的疫情，撲殺或無害化處理養(yǎng)殖場(chǎng)及其周邊鵝50 252只;2018年11月的疫情，共撲殺或無害化處理周邊家禽31 329羽[3]。這兩次疫情嚴(yán)重影響了揚(yáng)州市家禽養(yǎng)殖業(yè)發(fā)展。

揚(yáng)州市地處江蘇省中間地帶，是一個(gè)畜禽養(yǎng)殖、流通大市。2016年，揚(yáng)州市家禽出欄4315萬只，存欄1443萬只，產(chǎn)值達(dá)100多億元，而農(nóng)業(yè)總產(chǎn)值478.8億元，家禽占農(nóng)業(yè)比重超過20%，同時(shí)，每年又有超過3000萬只家禽通過流通領(lǐng)域進(jìn)入或運(yùn)離揚(yáng)州市[4]。而且揚(yáng)州市也是我國(guó)主要的羽絨加工基地，一旦感染禽流感疫情，后果不堪設(shè)想。

客觀上，揚(yáng)州市的水網(wǎng)分布密集而廣泛，養(yǎng)殖又以水禽鴨、鵝為主;經(jīng)過揚(yáng)州市的高郵湖，連接江蘇、安徽兩省，水域總面積達(dá)700 km2，野生水鳥眾多;揚(yáng)州處在候鳥東部遷徙線上，境內(nèi)的高郵湖、邵伯湖是候鳥的重要棲息地，也是候鳥和家禽的交匯點(diǎn)。綜上所述，揚(yáng)州市引起禽流感入侵的風(fēng)險(xiǎn)很大。

禽流感不但影響我國(guó)家禽業(yè)的發(fā)展，更事關(guān)公眾的身體健康和生命安全，因此對(duì)禽流感的風(fēng)險(xiǎn)評(píng)估與應(yīng)急控制的研究對(duì)完善現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)體系、保障公眾健康具有重要意義?？紤]到揚(yáng)州市獨(dú)有的養(yǎng)殖環(huán)境和生態(tài)體系，應(yīng)建立高致病性禽流感發(fā)生風(fēng)險(xiǎn)的定量評(píng)估模型，有效評(píng)估疫情發(fā)生的風(fēng)險(xiǎn)程度。

1 國(guó)內(nèi)外研究現(xiàn)狀

數(shù)學(xué)建模是傳染病研究中的重要方法。數(shù)學(xué)建模就是針對(duì)傳染病的傳播特點(diǎn)，收集現(xiàn)有傳染病發(fā)生的數(shù)據(jù)，在此基礎(chǔ)之上提出數(shù)學(xué)模型，然后運(yùn)用數(shù)學(xué)理論分析出影響傳染病傳播的關(guān)鍵因素，并借助計(jì)算機(jī)模擬技術(shù)揭示傳染病傳播的流行趨勢(shì)，從而有效預(yù)防和控制傳染病傳播。

1760年，瑞士數(shù)學(xué)家Daniel Bernoulli[5]最早使用數(shù)學(xué)方法研究傳染病天花的傳播途徑。1926年，Kemaek與Mekendriek[6]在分析了1665—1666年黑死病在倫敦的傳播規(guī)律和1906年瘟疫在孟買的傳播特點(diǎn)之后，構(gòu)造了著名的SIR倉(cāng)室模型，從此揭開了傳染病數(shù)學(xué)建模研究的新篇章。近年來，對(duì)高致病性禽流感傳播的數(shù)學(xué)建模研究成果屢見不鮮。2007年，Iwami等[7]建立SI-SIR模型，分析一類變異后可以在人與人、人與禽之間傳播的禽流感病毒，利用靈敏度分析和數(shù)值模擬提出了有效的防控措施，即隔離染病家禽與人的接觸，同時(shí)增加對(duì)染病禽類的檢疫。2009年，S.C.Inyama[8]建立了一類人-禽流感傳播的數(shù)學(xué)模型，計(jì)算了該模型無病平衡點(diǎn)的穩(wěn)定性，提出了非常有助于控制和減少疾病傳播的措施，即在禽流感爆發(fā)區(qū)域，地方政府應(yīng)該杜絕家禽的轉(zhuǎn)移，只要發(fā)現(xiàn)家禽感染禽流感就要盡快撲殺，同時(shí)合理的補(bǔ)償養(yǎng)殖戶的損失。2017年，Chen[9]建立了SEQS模型，通過對(duì)模型的研究提出建議，即禽流感疫情一旦發(fā)生，要立即對(duì)染病者進(jìn)行隔離，同時(shí)加大檢疫外來禽類的力度，必要時(shí)要停止外來禽類的輸入和本地禽類的輸出。

2013年以來，我國(guó)的禽流感疫情雖然整體可控，但由于禽流感病毒變異非?？?，傳播迅速，嚴(yán)重危害禽類養(yǎng)殖業(yè)，同時(shí)造成社會(huì)恐慌。因此，結(jié)合地域特點(diǎn)以及禽流感的染病機(jī)理，應(yīng)建立相應(yīng)的數(shù)學(xué)模型來研究其傳播途徑，以便于獲得有效措施防控禽流感。

2 揚(yáng)州市禽流感風(fēng)控模型項(xiàng)目進(jìn)展

2.1 調(diào)研情況

課題組從禽流感發(fā)生與傳播研究現(xiàn)狀、禽流感風(fēng)險(xiǎn)分析、禽流感影響研究、揚(yáng)州市禽流感研究進(jìn)展四個(gè)方面展開文獻(xiàn)搜集，了解了研究的前沿情況。在查閱文獻(xiàn)、分析文獻(xiàn)之后，聯(lián)系業(yè)內(nèi)專家，結(jié)合揚(yáng)州市的生態(tài)系統(tǒng)和養(yǎng)殖環(huán)境，設(shè)計(jì)了相應(yīng)的調(diào)查問卷，先后在邵伯湖區(qū)、高郵馬棚進(jìn)行了為期1個(gè)月的實(shí)地調(diào)查與農(nóng)戶訪談，在野外隨機(jī)采訪了31戶鴨農(nóng)。

根據(jù)整理分析調(diào)查資料，初步得出揚(yáng)州市禽流感的發(fā)生風(fēng)險(xiǎn)與候鳥、家禽、交通貿(mào)易及水稻田種植有一定關(guān)系。在此基礎(chǔ)之上，課題組結(jié)合揚(yáng)州市自然環(huán)境與經(jīng)濟(jì)社會(huì)特征，確定市禽流感發(fā)生風(fēng)險(xiǎn)因素。

2.2 實(shí)驗(yàn)情況

課題組聯(lián)合江蘇省家禽研究所的力量，通過利用簡(jiǎn)單隨機(jī)的抽樣方法，對(duì)揚(yáng)州市2017年9月到2018年1月采集的3個(gè)規(guī)模養(yǎng)雞場(chǎng)、30戶左右的散養(yǎng)戶，3個(gè)活禽市場(chǎng)、6個(gè)家禽屠宰點(diǎn)作為數(shù)據(jù)采集點(diǎn)的289份血清樣品進(jìn)行血清樣本檢測(cè)和病原學(xué)特征檢測(cè)，同時(shí)利用統(tǒng)計(jì)學(xué)方法對(duì)規(guī)模養(yǎng)殖場(chǎng)禽流感的致病因子進(jìn)行分析。

病原學(xué)檢測(cè)結(jié)果表明，各監(jiān)測(cè)場(chǎng)點(diǎn)主要以H9亞型禽流感病原為主，且每年10月至次年5月是病毒的活躍期。從不同區(qū)縣檢測(cè)結(jié)果來看，各區(qū)縣均檢測(cè)到了禽流感病毒，其中，江都區(qū)和高郵市為禽流感可能流行的主要區(qū)域。

3 構(gòu)建揚(yáng)州市禽流感風(fēng)險(xiǎn)評(píng)估模型的計(jì)劃

針對(duì)揚(yáng)州市特定的家禽養(yǎng)殖環(huán)境以及家禽交易、物流情況，本課題擬建立高致病性禽流感發(fā)生風(fēng)險(xiǎn)評(píng)估模型框架，以期實(shí)現(xiàn)評(píng)估禽流感疫情的發(fā)生風(fēng)險(xiǎn)。

3.1 禽流感風(fēng)控模型的研究思路

依據(jù)科學(xué)原則，審慎確定風(fēng)險(xiǎn)因素，利用層次分析法（AHP），建立高致病性禽流感發(fā)生的風(fēng)險(xiǎn)評(píng)估模型[1]。

3.2 禽流感風(fēng)控模型的構(gòu)建

建立高致病性禽流感免疫預(yù)防風(fēng)險(xiǎn)評(píng)估模型主要運(yùn)用層次分析法，其步驟有以下4點(diǎn)。1）建立評(píng)估體系。依據(jù)層次分析法的建模要求，把高致病性禽流感免疫預(yù)防控制作為目標(biāo)層，風(fēng)險(xiǎn)因素作為準(zhǔn)則層，子風(fēng)險(xiǎn)因素作為方案層，建立相應(yīng)的層次結(jié)構(gòu)評(píng)估指標(biāo)體系。2）確定指標(biāo)的單層權(quán)重。層次結(jié)構(gòu)體系確定了上、下層指標(biāo)元素間的隸屬關(guān)系。對(duì)于同層各元素，以相鄰上層有聯(lián)系的元素為準(zhǔn)，使用1～9比率標(biāo)度法衡量相對(duì)重要程度，并分別進(jìn)行兩兩比較，聘請(qǐng)8～10位長(zhǎng)期從事禽畜病學(xué)的專家，運(yùn)用評(píng)分辦法判斷其相對(duì)重要性或優(yōu)劣程度，從而構(gòu)造判斷矩陣。對(duì)各層的因素進(jìn)行對(duì)比分析，求解各因素的相對(duì)權(quán)重，其最大特征值對(duì)應(yīng)的特征向量歸一化后，即為該層次相應(yīng)元素（對(duì)于上一層次某一元素）的權(quán)重，依此進(jìn)行單層次評(píng)估研究。3）確定各指標(biāo)的整體權(quán)重。每一指標(biāo)的綜合權(quán)重值為所屬各層權(quán)重值之和，也表示該指標(biāo)的整體性綜合評(píng)估結(jié)果。4）建立評(píng)估模型。確定免疫預(yù)防風(fēng)險(xiǎn)等級(jí)，利用各指標(biāo)的整體權(quán)重計(jì)算總體風(fēng)險(xiǎn)程度，建立高致病性禽流感免疫預(yù)防風(fēng)險(xiǎn)的定量評(píng)估模型。

4 結(jié)語

對(duì)揚(yáng)州市禽流感風(fēng)險(xiǎn)評(píng)估與應(yīng)急控制模型進(jìn)行研究，能夠有效根據(jù)當(dāng)前養(yǎng)殖環(huán)境和免疫措施評(píng)判本區(qū)域高致病性禽流感發(fā)生風(fēng)險(xiǎn)，使揚(yáng)州市的決策者能夠評(píng)估價(jià)格快速制定出預(yù)防措施，減少財(cái)產(chǎn)損失，保護(hù)人身安全，使風(fēng)險(xiǎn)降至最低。

參考文獻(xiàn)：

[1] 王琦梅，劉社蘭，陳恩富.人感染H7N9禽流感流行病學(xué)研究進(jìn)展[J].中華預(yù)防醫(yī)學(xué)雜志，2017，51（2）：183-187.

[2] 搜狐新聞.農(nóng)業(yè)農(nóng)村部發(fā)布兩起高致病性禽流感疫情 如何應(yīng)對(duì)？[DB/OL].（2018-11-15）[2019-06-10].http：//m.sohu.com/a/277695280_649918.

[3] 中國(guó)新聞網(wǎng).江蘇揚(yáng)州江都區(qū)發(fā)生家禽H5N6亞型高致病性禽流感疫情[DB/OL].（2018-11-23）[2019-06-10].http：//sports.chinanews.com/sh/2018/11-23/8683795.shtml.

[4] 揚(yáng)州市統(tǒng)計(jì)局.2016年揚(yáng)州市國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)[N/OL].揚(yáng)州日?qǐng)?bào)，（2017-03-15）_[2019-06-10].http：//tjj.yangzhou.gov.cn/yztjj/tjjtjgb/201703/050b1fc83dbf4737b0e83114ecb87371.shtml.

[5] D.Bernoulli. Essaidune nouvelle analyse de la mortalitecausee par la petite verole et desavantages de linoculation pour al prevenir， in Memoires de Mathematiques et de physique[J].Paris：Academic Royale des Sciences. 1760：1-45.

[6] Kermack WO， MeKendriek AG. A Contribution to the mathematical theory of epidemics[J].Proceedings of the Royal Society of London. Series A， Containing Papers of a Mathematical and Physical Chearacter（1905-1934）. 1926， 115（772）：700-721.

[7] Shigo Iwami， Yasuhiro Takeuchi， Xianning Liu. Avian-human in?uenza epidemic model[J].Mathematical Biosciences.2007（207）：1-25.

[8] S.C.Inyama. Mathematical for bird ?u disease transimission with no bird migration[J]. globaljournal of mathematical sciences. 2009， 8（2）：75-81.

[9] Yong-xue Chen. Global analysis of an SEQS avian in?uenza model[J].journal of Fujian Agriculture and Forestry University（Natural Science Edition）.2010， 39（2）：173-176.

（責(zé)任編輯：劉昀）