動(dòng)物疫病防控策略EBCR分析法及其實(shí)證

廖 祺

（中南財(cái)經(jīng)政法大學(xué) 統(tǒng)計(jì)與數(shù)學(xué)學(xué)院，武漢 430073）

0 引言

動(dòng)物疫病給社會(huì)各個(gè)層面所造成的影響是巨大的。如何防控動(dòng)物疫病既是一個(gè)技術(shù)問(wèn)題，也是一個(gè)經(jīng)濟(jì)問(wèn)題。動(dòng)物疫病防控策略是指為了實(shí)現(xiàn)動(dòng)物疫病防控計(jì)劃的目標(biāo)，減輕動(dòng)物傳染病對(duì)人類社會(huì)的影響而采取的技術(shù)投入、制度安排等針對(duì)性的預(yù)防和控制方案或措施。對(duì)動(dòng)物疫病防控策略進(jìn)行經(jīng)濟(jì)學(xué)評(píng)估已經(jīng)逐漸成為各國(guó)制定動(dòng)物衛(wèi)生政策的一項(xiàng)基礎(chǔ)工作，它是指根據(jù)疫病所處的疫情階段，對(duì)當(dāng)期正在實(shí)施的策略與下一階段防控的方案策略進(jìn)行評(píng)估，從而為當(dāng)前防控策略是否延續(xù)、調(diào)整或終結(jié)提供良好的實(shí)證基礎(chǔ)和理論支撐，是逐步推進(jìn)防控工作的關(guān)鍵，也只有更加科學(xué)地安排和使用各種防控資金，才能提高有限衛(wèi)生資源的使用效率，更好地實(shí)現(xiàn)疫病的防控，因此，動(dòng)物疫病防控策略的經(jīng)濟(jì)學(xué)評(píng)估具有重要的實(shí)踐指導(dǎo)意義。

本文從動(dòng)物疫病防控資金的效用和效果分析的角度，將實(shí)施策略可能帶來(lái)的預(yù)期風(fēng)險(xiǎn)損失納入評(píng)估體系，從而構(gòu)建了EBCR分析法以評(píng)估動(dòng)物疫病的防控策略，是對(duì)這類風(fēng)險(xiǎn)決策問(wèn)題的初步探索。

1 方法的提出

動(dòng)物疾病防控策略經(jīng)濟(jì)學(xué)評(píng)估必須借助經(jīng)濟(jì)學(xué)的有關(guān)理論和分析方法，以病理學(xué)、免疫學(xué)、公共衛(wèi)生學(xué)、傳染病學(xué)及各種有關(guān)自然學(xué)科理論為依據(jù)，對(duì)動(dòng)物疫病防控過(guò)程中的投入產(chǎn)出、生產(chǎn)消費(fèi)以及風(fēng)險(xiǎn)損失等相關(guān)問(wèn)題進(jìn)行定性和定量分析。國(guó)內(nèi)外有關(guān)研究主要圍繞防控策略效果展開(kāi)，學(xué)者們通常對(duì)防控措施實(shí)施產(chǎn)生的直接結(jié)果、效用、收益和福利進(jìn)行分析。綜合分析近30年來(lái)動(dòng)物衛(wèi)生經(jīng)濟(jì)的研究軌跡，主要特點(diǎn)是：微觀研究較多而宏觀研究較少，單方面的研究多而綜合性的研究少，決策支持研究多而行為研究少，應(yīng)用研究較多而理論研究較少，沒(méi)有形成學(xué)科體系[1]。

在傳統(tǒng)的動(dòng)物疫病防控經(jīng)濟(jì)分析中，一般采用的是成本收益率法，即BCR（Benefit and Cost Ratio）分析法，主要通過(guò)對(duì)動(dòng)物疫病防控基礎(chǔ)成本投入和防控的短期效果進(jìn)行評(píng)估，但是動(dòng)物疫病防控決策一旦被實(shí)施就意味著至少產(chǎn)生兩方面的經(jīng)濟(jì)影響：一方面是實(shí)施具體策略的防控成本，另一方面是采用不同防控策略又會(huì)產(chǎn)生不同的風(fēng)險(xiǎn)損失，而直接影響著防控的產(chǎn)出實(shí)際效果。將風(fēng)險(xiǎn)分析納入經(jīng)濟(jì)學(xué)評(píng)估體系正是研究的發(fā)展趨勢(shì)。許多學(xué)者利用動(dòng)力學(xué)模型對(duì)包括動(dòng)物疫病在內(nèi)的傳染病進(jìn)行大量仿真研究，為動(dòng)物疫病防控策略的風(fēng)險(xiǎn)評(píng)估奠定了較好的基礎(chǔ)[2]。本文將風(fēng)險(xiǎn)概率導(dǎo)致的風(fēng)險(xiǎn)損失納入效益評(píng)估中，稱為擴(kuò)展的成本收益率法——EBCR（Extended Benefit and Cost Ratio）分析法。

EBCR分析法是以疫病防控基礎(chǔ)成本評(píng)估、暴發(fā)損失評(píng)估、暴發(fā)風(fēng)險(xiǎn)損失評(píng)估為動(dòng)物疫病防控策略經(jīng)濟(jì)學(xué)評(píng)估的基礎(chǔ)，以動(dòng)力學(xué)模型以及動(dòng)物衛(wèi)生經(jīng)濟(jì)學(xué)主要經(jīng)濟(jì)學(xué)評(píng)價(jià)方法中最常用的成本收益分析法為主要方法。評(píng)估的具體內(nèi)容由三項(xiàng)構(gòu)成：（1）基礎(chǔ)防控成本評(píng)估；（2）暴發(fā)損失評(píng)估；（3）風(fēng)險(xiǎn)損失評(píng)估。

我國(guó)對(duì)動(dòng)物疫病的防控管理實(shí)行動(dòng)物疫病分類管理制度。針對(duì)不同級(jí)別的動(dòng)物疫病，我國(guó)現(xiàn)有的主要防控策略一般可分為三種，即以撲殺為主的綜合防控策略、撲殺結(jié)合免疫的綜合防控策略和以免疫為主的綜合防控策略。針對(duì)這三種防控策略如何從經(jīng)濟(jì)性、實(shí)際效果等方面進(jìn)行評(píng)估和策略的轉(zhuǎn)換呢？本文以亞洲Ⅰ型口蹄疫為例進(jìn)行實(shí)證研究。

2 模型構(gòu)建、變量選取和數(shù)據(jù)說(shuō)明

2.1 EBCR模型的構(gòu)建

在對(duì)甲、乙兩種策略進(jìn)行經(jīng)濟(jì)學(xué)評(píng)估時(shí)，首先對(duì)兩種策略的不同防控措施成本（稱為基礎(chǔ)防控成本）進(jìn)行分析評(píng)估；其次，為了評(píng)估另一部分成本，即風(fēng)險(xiǎn)損失，可通過(guò)動(dòng)力學(xué)模型對(duì)不同策略行為導(dǎo)致的疫情影響進(jìn)行預(yù)測(cè)，然后在對(duì)該種疫病歷史暴發(fā)損失進(jìn)行經(jīng)濟(jì)學(xué)評(píng)估的基礎(chǔ)上，分析暴發(fā)損失與疫點(diǎn)數(shù)的關(guān)系，對(duì)風(fēng)險(xiǎn)損失進(jìn)行預(yù)測(cè)評(píng)估；最后，基于以上數(shù)據(jù)結(jié)論，使用成本收益率法進(jìn)行是否應(yīng)該退出免疫的決策評(píng)估。

假定當(dāng)前防控方案為甲策略，假設(shè)性解決方案為策略乙，采取乙策略的效益主要來(lái)自于減少的甲策略基礎(chǔ)防控成本，乙策略下的風(fēng)險(xiǎn)損失與甲策略下的風(fēng)險(xiǎn)損失之差表現(xiàn)為乙策略的成本，即，如果乙策略的風(fēng)險(xiǎn)損失大于甲策略的風(fēng)險(xiǎn)損失，則乙策略的成本構(gòu)成除了包括相應(yīng)的基礎(chǔ)防控措施成本，還包括由于實(shí)施乙策略而增加的風(fēng)險(xiǎn)損失；反之，則乙策略的收益除了包括減少的甲策略基礎(chǔ)防控成本，還包括由于實(shí)施乙策略而減少的風(fēng)險(xiǎn)損失。

Δ風(fēng)險(xiǎn)損失=乙策略風(fēng)險(xiǎn)損失-甲策略風(fēng)險(xiǎn)損失

在一定考察時(shí)期范圍m年內(nèi)：

若Δ風(fēng)險(xiǎn)損失＞0，則乙策略的EBCR為：

若Δ風(fēng)險(xiǎn)損失＜0，即相比較于甲策略，采取乙策略可以減少風(fēng)險(xiǎn)損失，減少的風(fēng)險(xiǎn)損失體現(xiàn)為采取乙策略的收益，乙策略的EBCR為：

EBCR＞1的策略為可行策略，即如果成本收益率EBCR=B/C＞1，乙策略的防控效益高于成本，該防控措施可以接受。但是數(shù)年考察期內(nèi)，在不同年份采取乙策略，其EBCR值有所差別，EBCR值為最大的年份被認(rèn)為對(duì)應(yīng)策略轉(zhuǎn)換最優(yōu)時(shí)機(jī)。

2.2 變量的選取

甲策略的基礎(chǔ)防控成本A1包括：疫病監(jiān)測(cè)調(diào)查成本a1和年均應(yīng)急儲(chǔ)備成本a2，前者由疫病監(jiān)測(cè)成本a11和疫病調(diào)查成本a12構(gòu)成，后者受到應(yīng)急儲(chǔ)備成本a21與理論儲(chǔ)備量a22這兩方面主要因素的影響。

乙策略的基礎(chǔ)防控成本B1包括：免疫計(jì)劃實(shí)施成本b1和疫病監(jiān)測(cè)調(diào)查成本b2，前者由疫苗成本b11和落實(shí)成本b12構(gòu)成，后者由疫病監(jiān)測(cè)成本b21和疫病監(jiān)測(cè)成本b22構(gòu)成。

甲、乙策略的風(fēng)險(xiǎn)損失分別為a3和b3。

2.3 數(shù)據(jù)來(lái)源及說(shuō)明

本文的數(shù)據(jù)主要來(lái)源于以下統(tǒng)計(jì)資料：

（1）動(dòng)物疫病疫情數(shù)據(jù)。國(guó)外動(dòng)物疫病疫情數(shù)據(jù)來(lái)自于世界動(dòng)物衛(wèi)生組織網(wǎng)站；國(guó)內(nèi)口蹄疫疫點(diǎn)數(shù)據(jù)來(lái)自《獸醫(yī)公報(bào)》和農(nóng)業(yè)部獸醫(yī)局。

（2）監(jiān)測(cè)數(shù)據(jù)。監(jiān)測(cè)方案參照《2015年國(guó)家動(dòng)物疫病監(jiān)測(cè)與流行病學(xué)調(diào)查計(jì)劃》；監(jiān)測(cè)數(shù)據(jù)來(lái)自農(nóng)業(yè)部獸醫(yī)局、《獸醫(yī)公報(bào)》和國(guó)家口蹄疫參考實(shí)驗(yàn)室。

（3）畜牧業(yè)養(yǎng)殖數(shù)據(jù)。牛、羊、駱駝的養(yǎng)殖數(shù)據(jù)來(lái)自《農(nóng)村統(tǒng)計(jì)年鑒》、《中國(guó)農(nóng)業(yè)統(tǒng)計(jì)年鑒》、《中國(guó)畜牧業(yè)統(tǒng)計(jì)年鑒》和國(guó)家統(tǒng)計(jì)局網(wǎng)站。

（4）強(qiáng)制免疫數(shù)據(jù)。口蹄疫強(qiáng)制免疫數(shù)據(jù)參照《2016年國(guó)家動(dòng)物疫病強(qiáng)制免疫計(jì)劃》、《2015—2022年中國(guó)畜禽疫苗市場(chǎng)評(píng)估及未來(lái)發(fā)展趨勢(shì)研究報(bào)告》。

（5）其他數(shù)據(jù)。疫苗價(jià)格、免疫副反應(yīng)率、疫苗注射費(fèi)用和儲(chǔ)存費(fèi)用數(shù)據(jù)來(lái)自文獻(xiàn)資料、現(xiàn)場(chǎng)調(diào)查和專家訪談；口蹄疫病原及感染相關(guān)參數(shù)來(lái)自專家意見(jiàn)、參考文獻(xiàn)，或由數(shù)據(jù)分析、參數(shù)估計(jì)和數(shù)據(jù)擬合得到。

本文分邊境與內(nèi)陸省份兩個(gè)區(qū)域進(jìn)行研究。根據(jù)疫情特點(diǎn)，內(nèi)陸地區(qū)包括：甘肅、青海、寧夏、四川、重慶、陜西、貴州、湖南、湖北、河南、山西、北京、安徽、江西、海南、河北、天津、山東、江蘇、上海、浙江、福建、廣東；邊境地區(qū)包括：內(nèi)蒙古、黑龍江、吉林、遼寧、廣西、云南、新疆、西藏。

3 實(shí)證分析——以亞洲I型口蹄疫的免疫策略為例

3.1 當(dāng)前防控策略問(wèn)題分析

口蹄疫（Foot and Mouth Disease，F(xiàn)MD）以感染偶蹄動(dòng)物為主，是因口蹄疫病毒引起的熱性、急性、高度傳染性動(dòng)物疫病，屬于一類動(dòng)物疫病，亞洲Ⅰ型口蹄疫是已知七個(gè)主要的口蹄疫類型之一。2005年，亞洲Ⅰ型口蹄疫曾擴(kuò)散到中國(guó)、俄羅斯和蒙古，但2010年至今，我國(guó)再?zèng)]有報(bào)告發(fā)生亞洲Ⅰ型臨床病例，也未監(jiān)測(cè)到病原學(xué)陽(yáng)性，我國(guó)亞洲Ⅰ型持續(xù)多年無(wú)疫的狀態(tài)[3]，易感動(dòng)物個(gè)體和群體抗體合格率長(zhǎng)期維持在80%以上，已經(jīng)符合《OIE陸生動(dòng)物衛(wèi)生法典》關(guān)于免疫無(wú)口蹄疫的規(guī)定。

我國(guó)《國(guó)家中長(zhǎng)期動(dòng)物疫病防治規(guī)劃》明確提出亞洲Ⅰ型口蹄疫下一階段防治目標(biāo)是達(dá)到凈化狀態(tài)，2020年實(shí)現(xiàn)全國(guó)非免疫無(wú)疫。因此，一方面，當(dāng)前采取的強(qiáng)制免疫為主的綜合防控策略已經(jīng)與凈化計(jì)劃的目標(biāo)不相匹配，停止強(qiáng)制免疫已成為防控計(jì)劃階段轉(zhuǎn)換的必然要求；另一方面，一種疫病實(shí)現(xiàn)從控制到消滅計(jì)劃通常歷時(shí)十多年，期間需要大量的經(jīng)費(fèi)作為支撐[4]。因此，本文所要討論的是，針對(duì)亞洲Ⅰ型口蹄疫，如何有計(jì)劃的對(duì)存欄動(dòng)物停止實(shí)施強(qiáng)制免疫，采取以退出免疫為主，加強(qiáng)監(jiān)測(cè)這種假設(shè)性策略方案從經(jīng)濟(jì)學(xué)角度分析是否可行。

3.2 退出強(qiáng)制免疫策略基礎(chǔ)防控成本

退出強(qiáng)制免疫策略基礎(chǔ)防控成本主要由監(jiān)測(cè)調(diào)查成本和年均應(yīng)急儲(chǔ)備成本構(gòu)成。

退出免疫情況下，為防控亞洲Ⅰ型口蹄疫，政府對(duì)邊境和內(nèi)陸地區(qū)必須投入的年均基礎(chǔ)成本經(jīng)核算為192354.12萬(wàn)元和292425.15萬(wàn)元。結(jié)果如下頁(yè)表1所示。

3.3 繼續(xù)強(qiáng)制免疫策略基礎(chǔ)防控成本

繼續(xù)強(qiáng)制免疫策略基礎(chǔ)防控成本主要由監(jiān)測(cè)調(diào)查成本和免疫計(jì)劃實(shí)施成本構(gòu)成。繼續(xù)強(qiáng)制免疫策略基礎(chǔ)成本核算統(tǒng)計(jì)結(jié)果見(jiàn)表2：

表1 退出免疫策略基礎(chǔ)成本

表2 繼續(xù)強(qiáng)制免疫策略基礎(chǔ)成本

4 風(fēng)險(xiǎn)評(píng)估

4.1 養(yǎng)殖場(chǎng)分類

本文以養(yǎng)牛場(chǎng)為例，將養(yǎng)牛場(chǎng)按照染病狀態(tài)分為7類，即易感場(chǎng)類（S）、感染場(chǎng)類、康復(fù)場(chǎng)類（R）、免疫場(chǎng)類和病毒攜帶場(chǎng)類（C）。其中，依據(jù)是否表現(xiàn)臨床癥狀，將感染場(chǎng)分為暫無(wú)癥狀類養(yǎng)殖場(chǎng)（E）、無(wú)癥狀養(yǎng)殖場(chǎng)（Y）。依據(jù)免疫后抗體水平，免疫場(chǎng)類分為具有部分抗體養(yǎng)殖場(chǎng)（V1）、完全具有抗體養(yǎng)殖場(chǎng)（V2）

4.2 模型流程圖

模型描述了口蹄疫病毒在不同類別養(yǎng)牛場(chǎng)之間的傳播過(guò)程以及不同類別養(yǎng)牛場(chǎng)數(shù)量的動(dòng)態(tài)變化情況，從而揭示口蹄疫傳播隨時(shí)間變化的內(nèi)在演化規(guī)律，流程圖見(jiàn)圖1。

圖1 口蹄疫在各類養(yǎng)殖場(chǎng)中傳播的流程圖

假設(shè)免疫有效率為k，單位時(shí)間內(nèi)kS個(gè)牛場(chǎng)進(jìn)入免疫倉(cāng)室V1。免疫保護(hù)期為6個(gè)月，單位時(shí)間內(nèi)將有λV2個(gè)場(chǎng)失去免疫力，從V2類進(jìn)入S類。易感場(chǎng)被暫無(wú)癥狀場(chǎng)感染成為暫無(wú)癥狀場(chǎng)（）。潛伏期過(guò)后，暫無(wú)癥狀場(chǎng)E，表現(xiàn)出臨床癥狀、進(jìn)入無(wú)癥狀類Y或者關(guān)閉。

攜帶類養(yǎng)殖場(chǎng)C可以成為無(wú)癥狀類養(yǎng)殖場(chǎng)Y，或者攜帶類養(yǎng)殖場(chǎng)C。無(wú)癥狀場(chǎng)Y和攜帶類C能夠恢復(fù)成為康復(fù)類養(yǎng)殖場(chǎng)R，康復(fù)類養(yǎng)殖場(chǎng)R也存在關(guān)閉的可能性。

具有部分抗體養(yǎng)殖場(chǎng)V1能夠轉(zhuǎn)化為完全具有抗體養(yǎng)殖場(chǎng)V2；或者受到病毒感染（）而進(jìn)入Y。具有完全抗體養(yǎng)殖場(chǎng)V2過(guò)了疫苗保護(hù)期后轉(zhuǎn)化為易感養(yǎng)殖場(chǎng)S。

4.3 微分方程模型

由流程圖1可以建立微分方程（1）：

4.4 系統(tǒng)演化

模型（1）描述了亞洲Ⅰ型口蹄疫從2005年4月份暴發(fā)以來(lái)的演化過(guò)程。以2007年各類養(yǎng)殖場(chǎng)數(shù)量為變量S、E、Y、C、R、V1和V2的初始條件，各類養(yǎng)殖場(chǎng)將按照模型（1）進(jìn)行演化。如果實(shí)施免疫退出計(jì)劃，即模型（1）中k（免疫率）為0，那么模型（1）將變?yōu)槟Ｐ停?）：

如果在時(shí)間點(diǎn)t0實(shí)施免疫退出計(jì)劃，原本按照模型（1）傳播的亞洲Ⅰ型口蹄疫，將轉(zhuǎn)化為按照模型（2）傳播。

4.5 參數(shù)估計(jì)

模型（1）中的參數(shù)主要來(lái)源于四個(gè)方面：一是文獻(xiàn)查詢；二是調(diào)查數(shù)據(jù)仿真；三是專家經(jīng)驗(yàn)；四是假設(shè)。在這個(gè)模型中，假設(shè)的參數(shù)只有1個(gè)即p，用于描述易感動(dòng)物感染口蹄疫病毒后不表現(xiàn)臨床癥狀的比例。參數(shù)值、實(shí)際意義和來(lái)源見(jiàn)表3：

表3 模型（1）中參數(shù)值和說(shuō)明

4.6 數(shù)據(jù)擬合

數(shù)值模擬包括參數(shù)估計(jì)和模型驗(yàn)證、基本再生數(shù)計(jì)算、亞洲Ⅰ型口蹄疫發(fā)生概率和規(guī)模計(jì)算以及參數(shù)敏感性分析。

采用MCMC參數(shù)估計(jì)方法，運(yùn)用函數(shù)D˙=(1-p)δE.擬合模型數(shù)值解和疫情數(shù)據(jù)。

通過(guò)擬合累計(jì)疫點(diǎn)數(shù)據(jù)，可以估計(jì)數(shù)參數(shù)向量M,q,a.。擬合全國(guó)疫點(diǎn)數(shù)據(jù)，獲取M1,q1,a1參數(shù)值；擬合內(nèi)地省份和邊境省份疫點(diǎn)數(shù)據(jù)，獲取M2,q2,a2，M3,q3,a3參數(shù)值，見(jiàn)表4。

表4 M1,q1,a1參數(shù)估計(jì)結(jié)果及置信區(qū)間

4.7 基本再生數(shù)

基本再生數(shù)分為免疫措施下基本再生數(shù)R0和退出免疫情況下基本再生數(shù)R1。

R0通過(guò)以下矩陣計(jì)算：

則基本再生數(shù)為：R0=ρ(FV-1)，其中ρ為譜半徑。

基本再生數(shù)詳見(jiàn)表5和下頁(yè)表6：

表5 免疫與免疫退出后基本再生數(shù)

表6 內(nèi)陸和邊境地區(qū)傳播閾值計(jì)算結(jié)果

5 風(fēng)險(xiǎn)損失值評(píng)估

5.1 動(dòng)物疫病暴發(fā)損失評(píng)估方法

動(dòng)物疫病的暴發(fā)對(duì)地區(qū)或者社會(huì)所造成的影響是多方面的。動(dòng)物疫病主要影響主體包括政府、社會(huì)公眾、生產(chǎn)企業(yè)、養(yǎng)殖場(chǎng)（戶）、銷售企業(yè)等；從影響層面看，疫病暴發(fā)主要造成了不同程度的社會(huì)、經(jīng)濟(jì)和環(huán)境層面的損失[5]。

損失評(píng)估常用的方法是收益現(xiàn)值法和現(xiàn)行市價(jià)法。在動(dòng)物產(chǎn)品加工和養(yǎng)殖銷售市場(chǎng)上，根據(jù)現(xiàn)行市價(jià)法分析計(jì)算的條件，可以用現(xiàn)行市價(jià)法對(duì)動(dòng)物產(chǎn)品加工和養(yǎng)殖銷售的各種交易行為進(jìn)行價(jià)值估算。對(duì)因動(dòng)物疫病所造成的政府損失、企業(yè)損失、公共衛(wèi)生損失、消費(fèi)者損失中的直接醫(yī)療費(fèi)用的計(jì)算，還有養(yǎng)殖場(chǎng)（戶）損失中的動(dòng)物免疫費(fèi)用、動(dòng)物產(chǎn)品價(jià)值損失、動(dòng)物死亡損失等的計(jì)算，都可以用現(xiàn)行市價(jià)法中的直接法進(jìn)行相應(yīng)的損失估計(jì)。本文實(shí)證分析中，對(duì)動(dòng)物疫病損失評(píng)估時(shí)，對(duì)人感染疫病死亡所造成的環(huán)境污染的環(huán)境損失評(píng)估、人力資本損失評(píng)估，采用收益現(xiàn)值法進(jìn)行計(jì)算。

5.2 動(dòng)物疫病損失評(píng)估指標(biāo)體系

在構(gòu)建具體評(píng)估指標(biāo)時(shí)，對(duì)指標(biāo)的選擇不僅要依據(jù)重要性原則，還應(yīng)考慮數(shù)據(jù)的可獲取性，盡量保障每項(xiàng)指標(biāo)都能獲取數(shù)據(jù)資料。根據(jù)口蹄疫暴發(fā)傳播機(jī)制和其流行病學(xué)特點(diǎn)，以及《口蹄疫防治技術(shù)規(guī)范》中的危機(jī)處理要求，構(gòu)建口蹄疫暴發(fā)損失評(píng)估指標(biāo)體系，如表7所示：

表7 口蹄疫暴發(fā)損失評(píng)估指標(biāo)體系

5.3 歷史暴發(fā)損失值評(píng)估

我國(guó)亞洲Ⅰ型口蹄疫暴發(fā)始于2005年4月，止于2009年6月。經(jīng)統(tǒng)計(jì)，邊境地區(qū)撲殺死亡數(shù)量為2432頭，而內(nèi)陸地區(qū)為5724頭。由亞洲Ⅰ型口蹄疫導(dǎo)致的死亡數(shù)，兩個(gè)地區(qū)由于亞洲Ⅰ型口蹄疫造成的平均損失值分別為2.72萬(wàn)元/頭和2.08萬(wàn)元/頭。根據(jù)前文指標(biāo)體系，經(jīng)統(tǒng)計(jì)核算，兩區(qū)域亞洲Ⅰ型口蹄疫歷史暴發(fā)總損失值如表8所示：

表8 2005—2009年亞洲Ⅰ型口蹄疫各區(qū)域歷史暴發(fā)損失值Li（萬(wàn)元）

5.4 風(fēng)險(xiǎn)損失值評(píng)估結(jié)果

對(duì)邊境和內(nèi)陸兩個(gè)地區(qū)暴發(fā)疫點(diǎn)數(shù)與相應(yīng)暴發(fā)損失值進(jìn)行回歸分析得出，邊境地區(qū)與內(nèi)陸地區(qū)亞洲Ⅰ型口蹄疫暴發(fā)疫點(diǎn)數(shù)x與損失值y之間相關(guān)系數(shù)分別為0.9749和0.8126，反映出損失值y與疫點(diǎn)數(shù)x線性相關(guān)程度較高。

邊境和內(nèi)陸地區(qū)在不同時(shí)間采取退出免疫，風(fēng)險(xiǎn)損失值評(píng)估結(jié)果如表9和表10所示：

表9 邊境地區(qū)退出免疫時(shí)間與暴發(fā)損失規(guī)模

表10 內(nèi)陸地區(qū)退出免疫時(shí)間與暴發(fā)損失規(guī)模

6 EBCR分析結(jié)果

退出免疫決策的成本主要是由退出免疫后的防控措施而付出的直接基礎(chǔ)成本構(gòu)成。收益則主要體現(xiàn)在退出免疫措施減少的大量的強(qiáng)制免疫基礎(chǔ)成本投入；減少的損失值視為動(dòng)物疫病防控決策的收益，由于退出免疫造成風(fēng)險(xiǎn)損失值的增加，這里視為收益的一種減少。因此，在一定考察時(shí)期范圍m年內(nèi)，第n（n=1,2,3,...，m）年退出免疫策略的EBCR為：

EBCR＞1的決策為可行決策，設(shè)定數(shù)年考察期內(nèi)，不同年份采取該策略EBCR大小有所差別，EBCR最大的決策也為時(shí)間上退出免疫的最優(yōu)選擇。

從年均基礎(chǔ)動(dòng)物疫病防控成本的比較上看，采取強(qiáng)制免疫措施每年投入的直接成本較高。

由于強(qiáng)制免疫下b3n=0，a3n＞0,即Δ風(fēng)險(xiǎn)損失=b3n-a3n＜0,則退出免疫作為乙策略，其增加的風(fēng)險(xiǎn)屬于成本的增量。第n年(n=1,2，3，...，6)退出免疫決策：

即在短期評(píng)估內(nèi)，將減少的強(qiáng)制免疫防控成本投入視為退出免疫決策的收益。以2017—2022年風(fēng)險(xiǎn)分析結(jié)果為參照，若在不同年份退出免疫，其EBCR不同，例如，邊境地區(qū)若2017年初采取退出免疫措施，則：

若2018年退出免疫，即2017年仍強(qiáng)制免疫，則：

以此類推，在2017—2022年考察期范圍內(nèi)，邊境地區(qū)在每年初采取退出免疫的EBCR值依次為：3.08、1.91、1.20、0.75、0.43、0.19。內(nèi)陸地區(qū)在每年初采取退出免疫策略的EBCR值依次為：3.21、2.25、1.35、0.81、0.45、0.19。

7 結(jié)束語(yǔ)

針對(duì)亞洲Ⅰ型口蹄疫防控策略的經(jīng)濟(jì)學(xué)評(píng)估結(jié)果顯示，從2017—2022年，不論是邊境地區(qū)或是內(nèi)陸地區(qū)，在2017—2019年三年內(nèi)采取退出免疫策略的EBCR值均大于1，退出免疫策略具有可行性，且越早采取退出免疫措施，該策略的EBCR值越高。建議我國(guó)盡早采取退出免疫策略，開(kāi)展凈化措施，保障2020年實(shí)現(xiàn)全國(guó)非免疫無(wú)疫目標(biāo)。

[1]魏穎.衛(wèi)生經(jīng)濟(jì)學(xué)與衛(wèi)生經(jīng)濟(jì)管理[M].北京:人民衛(wèi)生出版社,1998.

[2]謝仲倫,楊楠,黃保續(xù)等.2004.SARS控制策略仿真研究:防控措施效果比較[J].系統(tǒng)仿真學(xué)報(bào),2016,(12).

[3]鄧光明.口蹄疫流行態(tài)勢(shì)及國(guó)家防控策略[J].獸醫(yī)導(dǎo)刊,2016,(17).

[4]韋欣捷,陳雯雯,林萬(wàn)龍等.發(fā)達(dá)國(guó)家動(dòng)物疫病防控財(cái)政支持政策及啟示[J].農(nóng)業(yè)經(jīng)濟(jì)問(wèn)題,2011,(7).

[5]鄭雪光,滕翔雁,朱迪國(guó)等.全球重大動(dòng)物疫病狀況及影響[J].中國(guó)動(dòng)物檢疫,2014,(1).