李燕凌，王繼紅，蘇青松

(1.湖南農(nóng)業(yè)大學(xué)公共管理與法學(xué)學(xué)院，長(zhǎng)沙 410128；2.中國(guó)農(nóng)業(yè)大學(xué)人文與發(fā)展學(xué)院，北京 100093)

一、引 言

2013 年以來(lái)中國(guó)較大范圍流行H7N9 禽流感，到2018 年流行非洲豬瘟，再次引發(fā)公眾與學(xué)界對(duì)動(dòng)物養(yǎng)殖業(yè)經(jīng)濟(jì)損失空間分布及其影響因素的高度關(guān)注。董雯（2018）、趙璐（2014）、段宜嘉（2018）等學(xué)者研究表明，動(dòng)物疫病傳播空間（地理）分布呈現(xiàn)區(qū)域性特點(diǎn)。部分疫病遍及全球，而部分疫病僅局部地區(qū)流行。李燕凌（2015）研究指出全部動(dòng)物疫病分布均呈區(qū)域性特點(diǎn)。從疫病本身而言，動(dòng)物傳染病和非傳染病各有分布規(guī)律，受病因特點(diǎn)、宿主特征和環(huán)境條件、自然條件及社會(huì)條件差異等影響，表現(xiàn)出反地域性特點(diǎn)。正是動(dòng)物疫病傳播的區(qū)域性差異，客觀上反映了養(yǎng)殖業(yè)經(jīng)濟(jì)損失的空間分布不同，為利用其空間分布規(guī)律規(guī)避養(yǎng)殖業(yè)經(jīng)濟(jì)損失風(fēng)險(xiǎn)提供了機(jī)遇。因此，研究動(dòng)物疫病的地理分布特征，闡明其引起差異原因，有助于制定科學(xué)的經(jīng)濟(jì)避損及疫病傳播防控對(duì)策。

國(guó)內(nèi)外學(xué)者對(duì)動(dòng)物疫病流行風(fēng)險(xiǎn)造成的養(yǎng)殖業(yè)經(jīng)濟(jì)損失和區(qū)域差異有較深入研究。1991年世界動(dòng)物衛(wèi)生組織（OIE）將風(fēng)險(xiǎn)分析列入《陸生動(dòng)物衛(wèi)生法典》中，推動(dòng)動(dòng)物疫病風(fēng)險(xiǎn)分析在動(dòng)物衛(wèi)生管理決策中應(yīng)用。Cheney Shreve等（2016）從人類(lèi)、社會(huì)和環(huán)境多角度探討2001年英國(guó)口蹄疫疫情迅速擴(kuò)散的原因，并強(qiáng)調(diào)疫病風(fēng)險(xiǎn)管理對(duì)預(yù)防和控制疫病流行的重要性。MustafaAl-Zoughool 等（2015）根據(jù)瘋牛病流行病學(xué)特征和疫病傳播特點(diǎn)，構(gòu)建瘋牛病流行風(fēng)險(xiǎn)數(shù)學(xué)模型，估算瘋牛病流行風(fēng)險(xiǎn)和未來(lái)走向。Yang（2010）綜述了數(shù)理統(tǒng)計(jì)模型等風(fēng)險(xiǎn)評(píng)估方法在媒價(jià)傳播性疫病中的應(yīng)用前景。我國(guó)動(dòng)物衛(wèi)生分析體系的構(gòu)建工作起步較晚，但發(fā)展迅速。周曉農(nóng)等（2014）構(gòu)建區(qū)域疫病流行現(xiàn)狀、潛在傳播風(fēng)險(xiǎn)和機(jī)構(gòu)工作能力等指標(biāo)，計(jì)算疫病流行風(fēng)險(xiǎn)評(píng)估指數(shù)對(duì)風(fēng)險(xiǎn)區(qū)域分類(lèi)，進(jìn)而采用ArcGIS 軟件繪制疫病流行風(fēng)險(xiǎn)地圖。李滋睿等（2010）提出動(dòng)物疫病流行指數(shù)計(jì)算公式，根據(jù)2004～2008五年間重大動(dòng)物疫病流行數(shù)據(jù)計(jì)算各區(qū)域疫病流行指數(shù)，將動(dòng)物疫病的流行區(qū)域分成嚴(yán)重流行區(qū)、較重流行區(qū)、中度流行區(qū)、散發(fā)區(qū)和潔凈區(qū)五個(gè)流行風(fēng)險(xiǎn)等級(jí)。

國(guó)外針對(duì)日益頻繁、損害威脅嚴(yán)重的動(dòng)物疫情公共危機(jī)演化機(jī)理的復(fù)雜性研究不足，特別是對(duì)動(dòng)物疫病流行風(fēng)險(xiǎn)造成的養(yǎng)殖業(yè)經(jīng)濟(jì)損失的區(qū)域差異研究文獻(xiàn)較少。國(guó)內(nèi)學(xué)者主要集中在區(qū)域差異影響因素和研究方法兩方面。李鵬等（2014）認(rèn)為動(dòng)物疫病發(fā)生和流行具有不確定性和可控性。不確定性主要體現(xiàn)在疫病發(fā)生和流行遵循一定的自然規(guī)律，不以人的意志為轉(zhuǎn)移，可控性體現(xiàn)在疫病發(fā)生和流行受復(fù)雜的自然因素和社會(huì)因素影響，傳染源、傳播途徑和易感動(dòng)物三個(gè)基本環(huán)節(jié)在特定情況下通過(guò)一定手段可切斷。對(duì)動(dòng)物疫病的流行具體受哪些自然因素和社會(huì)因素影響，導(dǎo)致形成區(qū)域差異，梁琛等（2009）總結(jié)地理生態(tài)環(huán)境對(duì)人及動(dòng)物健康的影響長(zhǎng)遠(yuǎn)，特別是疫源地地理生態(tài)環(huán)境對(duì)疫病空間分布和流行趨勢(shì)影響巨大。楊林生等（2010）認(rèn)為地理生態(tài)環(huán)境對(duì)動(dòng)物疫病流行造成影響。廢水大量排放會(huì)污染和破壞環(huán)境，畜禽廢棄物吸納能力也將大大衰減，使疫病病原體的傳播幾率增大。梁小珍等（2013）認(rèn)為，地區(qū)動(dòng)物飼養(yǎng)生產(chǎn)狀況不僅體現(xiàn)該地區(qū)養(yǎng)殖傳統(tǒng)和養(yǎng)殖業(yè)發(fā)展水平，一定程度上反映由養(yǎng)殖業(yè)發(fā)展趨勢(shì)引起的疫病流行風(fēng)險(xiǎn)變動(dòng)。王華等（2013）認(rèn)為地區(qū)動(dòng)物養(yǎng)殖水平越高，動(dòng)物疫病發(fā)生和流行風(fēng)險(xiǎn)也越大。與之截然相反觀點(diǎn)認(rèn)為，地區(qū)動(dòng)物飼養(yǎng)生產(chǎn)水平越高、規(guī)?；曫B(yǎng)比例越高，動(dòng)物疫病發(fā)生和流行風(fēng)險(xiǎn)越低。劉芳等（2011）認(rèn)為防疫人員素質(zhì)是疫病流行風(fēng)險(xiǎn)區(qū)域差異的重要影響因素，應(yīng)不斷加強(qiáng)防疫隊(duì)伍建設(shè)。

綜上，選取指數(shù)研究法、ArcGIS標(biāo)準(zhǔn)差橢圓分析法對(duì)2012～2016年中國(guó)動(dòng)物養(yǎng)殖業(yè)經(jīng)濟(jì)損失的空間格局演化展開(kāi)定量分析，利用STATA回歸分析法探究其影響因素。

二、中國(guó)動(dòng)物疫病經(jīng)濟(jì)損失的區(qū)域差異

由于中國(guó)養(yǎng)殖業(yè)經(jīng)濟(jì)損失主要由動(dòng)物疫病引起，動(dòng)物疫情公共危機(jī)若得不到及時(shí)有效管控，會(huì)引發(fā)食品安全、市場(chǎng)供求失衡、社會(huì)恐慌等系列社會(huì)和經(jīng)濟(jì)問(wèn)題，為針對(duì)動(dòng)物疫病造成的經(jīng)濟(jì)損失加以區(qū)劃，須構(gòu)建可全面反映該地區(qū)動(dòng)物疫病流行的指標(biāo)。分析動(dòng)物疫病流行情況，首先了解動(dòng)物疫病發(fā)生次數(shù)、發(fā)病動(dòng)物數(shù)、病死動(dòng)物數(shù)和該地區(qū)動(dòng)物總頭數(shù)①數(shù)據(jù)來(lái)源于中國(guó)農(nóng)業(yè)部《獸醫(yī)公報(bào)》。，這些是動(dòng)物疫病流行病學(xué)基本指標(biāo)，可用其構(gòu)建出反映各地區(qū)動(dòng)物疫病流行情況總指數(shù)。本文通過(guò)計(jì)算動(dòng)物疫病流行指數(shù)表示某個(gè)地區(qū)動(dòng)物疫病流行情況。借鑒李滋睿等（2010）動(dòng)物疫病流行指數(shù)公式如下：

式（1）中，E是指某種疫病在某一地區(qū)流行指數(shù)，Ni、Oi、Di、Ki分別指某種病在某地區(qū)發(fā)生次數(shù)、發(fā)病動(dòng)物數(shù)、死亡動(dòng)物數(shù)和撲殺動(dòng)物數(shù)，是指某種動(dòng)物疫病在全國(guó)總發(fā)生次數(shù)、總發(fā)病動(dòng)物數(shù)、總死亡動(dòng)物數(shù)和總撲殺動(dòng)物數(shù)， Ai和分別指當(dāng)年某地區(qū)和全國(guó)動(dòng)物總飼養(yǎng)量②數(shù)據(jù)來(lái)源于《中國(guó)畜牧獸醫(yī)年鑒》。。用平均流行指數(shù)表示一個(gè)地區(qū)所有動(dòng)物疫病的流行情況，計(jì)算公式如下：

式（2）中，E指某地區(qū)動(dòng)物疫病平均流行指數(shù)，E1，E2，…，En分別指該地區(qū)發(fā)生各種動(dòng)物疫病流行指數(shù)。在我國(guó)動(dòng)物疫病防治實(shí)踐中，重大動(dòng)物疫病對(duì)畜牧業(yè)生產(chǎn)和社會(huì)公共衛(wèi)生等領(lǐng)域危害極大，因此，防治重大動(dòng)物疫病的流行是獸醫(yī)工作重點(diǎn)。目前，對(duì)我國(guó)畜牧業(yè)造成重大危害的動(dòng)物疫病有口蹄疫、禽流感、新城疫、豬瘟、豬蘭耳病等數(shù)十種。為從總體上把握我國(guó)動(dòng)物疫病經(jīng)濟(jì)損失的空間分布，選取口蹄疫、禽流感、新城疫、豬瘟、豬蘭耳病五種最主要?jiǎng)游镆卟?，作為重大疫病區(qū)劃的代表病種。選取理由如下：第一，這五種病是目前我國(guó)畜牧業(yè)威脅最大的動(dòng)物疫病，且近幾年在我國(guó)不同區(qū)劃均有大流行。雖然2014年爆發(fā)的小反芻獸疫影響較大，全國(guó)22省市均有發(fā)生，但農(nóng)業(yè)部公開(kāi)數(shù)據(jù)多為數(shù)省合并病例數(shù)，且2012、2013年幾乎無(wú)病例數(shù)報(bào)告，故不便于采用分省數(shù)據(jù)比較分析。第二，豬（口蹄疫、豬瘟、豬藍(lán)耳病）、牛（口蹄疫）及禽類(lèi)（禽流感、新城疫）為五種患疫病動(dòng)物主體，占畜牧業(yè)80%以上。第三，雖然不同動(dòng)物疫病發(fā)病規(guī)律和影響因素存在一定差別，但本文不是從單種疫病研究其流行病學(xué)規(guī)律，而是考慮動(dòng)物疫病流行在地理分布上的經(jīng)濟(jì)、社會(huì)及地理因素影響。因豬、牛及禽類(lèi)在我國(guó)各省區(qū)均有養(yǎng)殖，可較全面反映動(dòng)物疫病全國(guó)流行特點(diǎn)和區(qū)域性特征。因此，開(kāi)展動(dòng)物養(yǎng)殖業(yè)經(jīng)濟(jì)損失的空間格局變化研究，忽略個(gè)別動(dòng)物疫病病種之間差異。第四，口蹄疫、禽流感、新城疫、豬瘟、豬蘭耳病分別是豬、牛及禽類(lèi)的主要危害疫病，故成為重點(diǎn)防疫疫病。

通過(guò)計(jì)算2012～2016年各省區(qū)三種疫病流行指數(shù)，取平均值得到各省區(qū)動(dòng)物疫病流行指數(shù)，反映省際動(dòng)物疫病連續(xù)5年流行狀況（見(jiàn)表1）。因此，可用動(dòng)物疫病流行指數(shù)這一指標(biāo)對(duì)全國(guó)動(dòng)物疫病進(jìn)行區(qū)劃。

表1 31?。ㄖ陛犑?、自治區(qū)）2012～2016年重大動(dòng)物疫病流行指數(shù)

續(xù)表1

世界糧農(nóng)組織（FAO）和世界動(dòng)物衛(wèi)生組織（OIE）在制定動(dòng)物疫病風(fēng)險(xiǎn)分級(jí)制度中，將動(dòng)物疫病流行風(fēng)險(xiǎn)分為極高風(fēng)險(xiǎn)、高風(fēng)險(xiǎn)、中等風(fēng)險(xiǎn)和低風(fēng)險(xiǎn)四類(lèi)。結(jié)合動(dòng)物疫病風(fēng)險(xiǎn)分級(jí)制度并參考李滋睿等（2010）分級(jí)方法計(jì)算出2012～2016年中國(guó)大陸31省區(qū)動(dòng)物疫病平均流行指數(shù)，再將31省區(qū)分為極高風(fēng)險(xiǎn)區(qū)域（E≥8）、高風(fēng)險(xiǎn)區(qū)域（8≥E≥4）、中等風(fēng)險(xiǎn)區(qū)域（4≥E≥1）和低風(fēng)險(xiǎn)區(qū)域（1≥E≥0）四個(gè)不同風(fēng)險(xiǎn)等級(jí)的動(dòng)物疫病流行風(fēng)險(xiǎn)區(qū)域。由表1 可知，天津、江西、重慶、貴州、云南、西藏、陜西、甘肅、寧夏、新疆10省區(qū)為極高風(fēng)險(xiǎn)區(qū)域；浙江、湖北、廣西、青海4省區(qū)為高風(fēng)險(xiǎn)區(qū)域；北京、河北、吉林、黑龍江、上海、江蘇、安徽、河南、湖南、廣東、四川、遼寧12省區(qū)為中等風(fēng)險(xiǎn)區(qū)域；山西、內(nèi)蒙古、福建、山東、海南的5省區(qū)為低風(fēng)險(xiǎn)區(qū)域。

三、中國(guó)動(dòng)物疫病經(jīng)濟(jì)損失的空間分布分析

（一）動(dòng)物疫病經(jīng)濟(jì)損失的空間分布分析方法

1.疫情空間分布格局描述

為探索2012～2016年疫情波及范圍和空間分布格局，利用ArcGIS10.1軟件中標(biāo)準(zhǔn)差橢圓工具對(duì)每一年病例分布開(kāi)展標(biāo)準(zhǔn)差橢圓分析。SDE方法通過(guò)以中心、長(zhǎng)軸、短軸、方位角為基本參數(shù)的空間分布橢圓，定量描述研究對(duì)象的空間分布整體特征。具體而言，空間分布橢圓以地理要素空間分布的平均中心為中心，分別計(jì)算其在X方向和Y方向上的標(biāo)準(zhǔn)差，以此定義包含要素分布的橢圓軸。使用該橢圓可查看要素分布是否被拉長(zhǎng)，由此確定特定方向。SDE方法可從全局與空間角度定量解釋地理要素空間分布的中心性、展布性、方向性及空間形態(tài)等特征。橢圓空間分布范圍表示地理要素空間分布的主體區(qū)域，其中中心表示地理要素在二維空間上分布的相對(duì)位置，方位角反映其分布主趨勢(shì)方向（即正北方向順時(shí)針旋轉(zhuǎn)到橢圓長(zhǎng)軸的角度），長(zhǎng)軸表征地理要素在主趨勢(shì)方向上的離散程度。

2.SDE主要參數(shù)

計(jì)算公式如下。

3.標(biāo)準(zhǔn)差橢圓分析結(jié)果

將標(biāo)準(zhǔn)差數(shù)設(shè)為一個(gè)標(biāo)準(zhǔn)差，通過(guò)以平均中心、長(zhǎng)軸、短軸、方位角為基本參數(shù)的空間分布橢圓，定量描述2012～2016年病例點(diǎn)的空間分布整體特征，分析結(jié)果可生成一個(gè)包含每年病例點(diǎn)空間分布特性的橢圓面。通過(guò)對(duì)標(biāo)準(zhǔn)差橢圓分析結(jié)果中平均中心、長(zhǎng)軸、短軸、方位角等基本參數(shù)變化的解讀，可發(fā)現(xiàn)動(dòng)物疫病病例點(diǎn)的空間擴(kuò)散趨勢(shì)。

（二）2012～2016年主要?jiǎng)游镆卟×餍袛?shù)據(jù)及SDE模型建構(gòu)

1.主要?jiǎng)游镆卟×餍袛?shù)據(jù)

根據(jù)農(nóng)業(yè)部《獸醫(yī)公報(bào)》數(shù)據(jù)，統(tǒng)計(jì)2012～2016年全國(guó)31省區(qū)禽流感、新城疫、豬藍(lán)耳、豬瘟及口蹄疫感染病例數(shù)，制作五種主要?jiǎng)游镆卟×餍袛?shù)據(jù)。

2.主要?jiǎng)游镆卟DE計(jì)算結(jié)果

將2012～2016年各省五種主要疫病感染病例數(shù)利用ArcGIS10.1軟件中標(biāo)準(zhǔn)差橢圓分析工具，計(jì)算平均中心、長(zhǎng)軸、短軸、方位角等基本參數(shù)定量描述每年病例點(diǎn)空間分布整體特征。橢圓空間分布范圍表示地理要素空間分布的主體區(qū)域，其中中心表示地理要素二維空間分布的相對(duì)位置，平均中心以X和Y坐標(biāo)確定2012～2016年各年份五種疫情爆發(fā)中心點(diǎn)所在地區(qū)，長(zhǎng)、短軸表征地理要素主趨勢(shì)方向的離散程度，方位角反映其分布主趨勢(shì)方向（即正北方向順時(shí)針旋轉(zhuǎn)到橢圓長(zhǎng)軸的角度），分析計(jì)算結(jié)果見(jiàn)表2。

表2 2012～2016年五種疫病標(biāo)準(zhǔn)差橢圓分析計(jì)算結(jié)果

（三）2012～2016年主要?jiǎng)游镆卟DE模型結(jié)果分析

1.禽流感疫情空間分布格局分析

2012～2016年橢圓平均中心及覆蓋范圍呈現(xiàn)由中心向東北、東南方向移動(dòng)的軌跡。疫情分布總體表現(xiàn)為“西北-東北，東南”的空間分布格局。2012年禽流感疫情標(biāo)準(zhǔn)差橢圓平均中心在甘肅省天水市（X 坐標(biāo)=105.88 m，Y 坐標(biāo)=34.40 m），標(biāo)準(zhǔn)差橢圓方位角為135.36°。2013 年在湖北省恩施土家族苗族自治州（X坐標(biāo)=109.00 m，Y坐標(biāo)=30.15 m），標(biāo)準(zhǔn)差橢圓方位角為35.26°。2014年在河南省南陽(yáng)市（X坐標(biāo)=112.30 m，Y坐標(biāo)=32.80 m），標(biāo)準(zhǔn)差橢圓方位角為48.70°。2015年在安徽省安慶市（X坐標(biāo)=116.10 m，Y坐標(biāo)=30.30 m），標(biāo)準(zhǔn)差橢圓方位角為60.45°。2016年在青海省黃南藏族自治州（X坐標(biāo)=101.40 m，Y坐標(biāo)=35.10 m），標(biāo)準(zhǔn)差橢圓方位角為124.80°。

通過(guò)確定2012～2016年禽流感疫情發(fā)生的中心城市與每年方位角轉(zhuǎn)變得出分析結(jié)果：五年來(lái)禽流感疫情多覆蓋在東南地區(qū)，傳播方向?yàn)椤白晕鞅毕驏|南再向西北”回旋趨勢(shì)，2014年有向東北和偏南部傳播趨勢(shì)。疫情覆蓋區(qū)域呈“縮小—擴(kuò)大—再縮小—再擴(kuò)大”反反復(fù)復(fù)模式，疫情控制稍有疏忽便會(huì)呈現(xiàn)回溫之勢(shì)。可見(jiàn)即使疫情得到控制仍不能懈怠。

2.新城疫疫情空間分布格局分析

2012～2016年橢圓平均中心及覆蓋范圍呈現(xiàn)由中心向東南方向移動(dòng)的軌跡。疫情分布總體表現(xiàn)為“中部-東南”的空間分布格局。2012年新城疫疫情標(biāo)準(zhǔn)差橢圓平均中心在陜西省漢中市（X坐標(biāo)=107.57 m，Y坐標(biāo)=32.64 m），標(biāo)準(zhǔn)差橢圓方位角為131.14°。2013年在湖南省湘西土家族苗族自治州（X坐標(biāo)=109.60 m，Y坐標(biāo)=27.98 m），標(biāo)準(zhǔn)差橢圓方位角為113.13°。2014年在湖南省邵陽(yáng)市（X坐標(biāo)=110.48 m，Y坐標(biāo)=26.70 m），標(biāo)準(zhǔn)差橢圓方位角為44.90°。2015年在湖南省邵陽(yáng)市城步苗族自治縣（X坐標(biāo)=110.30 m，Y坐標(biāo)=26.30 m），標(biāo)準(zhǔn)差橢圓方位角為47.84°。2016年在湖南省常德市（X坐標(biāo)=111.40 m，Y坐標(biāo)=29.30 m），標(biāo)準(zhǔn)差橢圓方位角為81.84°。

通過(guò)確定2012～2016年新城疫疫情發(fā)生的中心城市與每年方位角轉(zhuǎn)變得出分析結(jié)果：五年來(lái)新城疫疫情有四年平均中心發(fā)生在湖南，僅2012年發(fā)生在陜西，但距湖南較近。可見(jiàn)2012～2016年新城疫疫情大多由湖南向周邊省份傳播，傳播范圍不廣，多集中在湖南省周邊省份，覆蓋范圍呈現(xiàn)“由湖南省向周邊城市”發(fā)散的圓形區(qū)域，且2012～2016年覆蓋面積幾乎相同，因此，控制新城疫疫情應(yīng)嚴(yán)抓疫情發(fā)源地，特別是湖南等周邊地區(qū)，目前嚴(yán)控新城疫疫病的傳染源是遏制疫病的最好方法。

3.豬藍(lán)耳疫情空間分布格局分析

2012～2016年橢圓平均中心及覆蓋范圍集中在東南部，呈現(xiàn)由湖南、湖北向四周擴(kuò)散的軌跡。疫情分布總體表現(xiàn)為“湖南、湖北-周邊各省份”的空間分布格局。2012年豬藍(lán)耳疫情標(biāo)準(zhǔn)差橢圓平均中心在湖北省武漢市（X坐標(biāo)=114.20 m，Y坐標(biāo)=30.36 m），標(biāo)準(zhǔn)差橢圓方位角為82.70°。2013年在湖北省恩施土家族苗族自治州（X 坐標(biāo)=110.16 m，Y 坐標(biāo)=30.65 m），標(biāo)準(zhǔn)差橢圓方位角為61.80°。2014 年在河南省信陽(yáng)市（X 坐標(biāo)= 115.39 m，Y 坐標(biāo)= 31.90 m），標(biāo)準(zhǔn)差橢圓方位角為171.36°。2015 年在江西省九江市（X 坐標(biāo)= 116.03 m，Y 坐標(biāo)= 28.90 m），標(biāo)準(zhǔn)差橢圓方位角為81.10°。2016 年在貴州省銅仁市（X 坐標(biāo)= 109.30 m，Y 坐標(biāo)= 27.68 m），標(biāo)準(zhǔn)差橢圓方位角為69.00°。

通過(guò)確定2012～2016年豬藍(lán)耳疫情發(fā)生的中心城市與每年方位角轉(zhuǎn)變得出分析結(jié)果：五年來(lái)豬藍(lán)耳疫情有兩年平均中心發(fā)生在湖北，其余三年均在湖北周邊城市回旋。疫病覆蓋范圍呈現(xiàn)“擴(kuò)大—縮小—再擴(kuò)大”趨勢(shì)，呈圓圈狀向周邊發(fā)散。因此，要控制豬藍(lán)耳疫情應(yīng)加大湖北等周邊省份疫情防治，從源頭控制疫情，防止繼續(xù)擴(kuò)散。

4.豬瘟疫情空間分布格局分析

2012～2015年橢圓平均中心及覆蓋范圍集中在南部，呈現(xiàn)由南向東北方向逐漸轉(zhuǎn)移，到2016年又轉(zhuǎn)移到南部的軌跡。疫情分布總體呈現(xiàn)為“南部-東北-南部”的空間分布格局。2012年豬瘟疫情標(biāo)準(zhǔn)差橢圓平均中心在貴州省黔西南州（X 坐標(biāo)=105.14 m，Y 坐標(biāo)=26.10 m），標(biāo)準(zhǔn)差橢圓方位角為62.10°。2013 年在貴州省遵義市（X 坐標(biāo)= 107.67 m，Y 坐標(biāo)= 28.57 m），標(biāo)準(zhǔn)差橢圓方位角為49.70°。2014年在廣西壯族自治區(qū)百色市（X坐標(biāo)=105.32 m，Y坐標(biāo)=24.88 m），標(biāo)準(zhǔn)差橢圓方位角為137.70°。2015 年在河北省唐山市（X 坐標(biāo)=118.30 m，Y 坐標(biāo)=39.57 m），標(biāo)準(zhǔn)差橢圓方位角為53.62°。2016 年在四川省宜賓市（X 坐標(biāo)=105.30 m，Y 坐標(biāo)=28.20 m），標(biāo)準(zhǔn)差橢圓方位角為43.10°。

通過(guò)確定2012～2016年豬瘟疫情發(fā)生的中心城市與每年方位角轉(zhuǎn)變得出分析結(jié)果：五年來(lái)疫病發(fā)生主要區(qū)域在中南方向，2015年有“中南—東北”擴(kuò)散趨勢(shì)，在北方疫情可控，2016年又移回到中南部地區(qū)。可見(jiàn)，豬瘟在北方易于控制在中南部則易于傳播。五年內(nèi)豬瘟疫情覆蓋面積呈現(xiàn)“擴(kuò)大—縮小—再擴(kuò)大”趨勢(shì)。

5.口蹄疫疫情空間分布格局分析

因2012～2016年口蹄疫疫情發(fā)生地區(qū)較為離散，無(wú)法分析其分布規(guī)律。

四、中國(guó)動(dòng)物養(yǎng)殖業(yè)經(jīng)濟(jì)損失差異的影響因素

動(dòng)物疫病發(fā)生和流行受復(fù)雜的自然因素和社會(huì)因素影響，呈現(xiàn)區(qū)域性差異。綜合相關(guān)文獻(xiàn)（Cheney Shreveetal，2016；MustafaAl-Zoughooletal，2015），采用面板數(shù)據(jù)模型，從地理生態(tài)環(huán)境、動(dòng)物飼養(yǎng)生產(chǎn)水平和動(dòng)物疫病防控管理等方面展開(kāi)因素分析。

（一）變量定義

借鑒王東陽(yáng)等（2009）研究，從地理生態(tài)環(huán)境、動(dòng)物飼養(yǎng)生產(chǎn)水平和動(dòng)物疫病防控管理等方面設(shè)置相關(guān)變量，利用前文計(jì)算獲取的動(dòng)物疫病分布主趨勢(shì)方向結(jié)果（標(biāo)準(zhǔn)差橢圓方位角）作為被解釋變量，設(shè)為Y，影響動(dòng)物養(yǎng)殖業(yè)經(jīng)濟(jì)損失空間分布結(jié)構(gòu)變化的各個(gè)影響因素為解釋變量，分別設(shè)為森林覆蓋率（x1）、單位土地面積廢水排放量（x2）、動(dòng)物規(guī)?；B(yǎng)殖戶(hù)數(shù)比例（x3）、出欄率（x4）、人均動(dòng)物飼養(yǎng)量（x5）、畜牧獸醫(yī)站建設(shè)比例（x6）、技術(shù)人員比例（x7）。相關(guān)變量及其定義見(jiàn)表3。

表3 動(dòng)物疫病流行風(fēng)險(xiǎn)區(qū)域差異的影響因素

1.地理生態(tài)環(huán)境

動(dòng)物疫病流行與地理生態(tài)環(huán)境因素緊密相關(guān)。影響流行過(guò)程的地理生態(tài)環(huán)境因素主要包括地理環(huán)境、土壤、水、空氣等，對(duì)傳染源、傳播途徑和易感動(dòng)物三個(gè)環(huán)節(jié)的作用錯(cuò)綜復(fù)雜，如一定的地理?xiàng)l件（海、河，高山等）對(duì)傳染源的轉(zhuǎn)移產(chǎn)生一定限制，成為天然隔離屏障。我國(guó)在自然屏障較好的海南島、四川盆地、遼東半島等六省五區(qū)建立無(wú)規(guī)定動(dòng)物疫病示范區(qū)，充分考慮地域自然屏障狀況。因此，本研究選取森林覆蓋率、單位土地面積廢水排放量作為地理生態(tài)環(huán)境因素的基礎(chǔ)指標(biāo)。

2.動(dòng)物飼養(yǎng)生產(chǎn)水平

本文主要通過(guò)三個(gè)指標(biāo)評(píng)價(jià)地區(qū)動(dòng)物飼養(yǎng)生產(chǎn)水平，研究各指標(biāo)對(duì)動(dòng)物疫病流行風(fēng)險(xiǎn)的影響。一是規(guī)模化養(yǎng)殖程度，用動(dòng)物規(guī)?；B(yǎng)殖戶(hù)數(shù)比來(lái)衡量，比例越高規(guī)?；B(yǎng)殖程度越高；二是動(dòng)物生產(chǎn)力水平，由出欄率衡量，出欄率越高動(dòng)物生產(chǎn)水平越高；三是養(yǎng)殖密度，由人均動(dòng)物飼養(yǎng)量衡量，人均動(dòng)物飼養(yǎng)量越高養(yǎng)殖密度越高。

3.動(dòng)物疫病防控管理

2014年農(nóng)業(yè)部發(fā)布的《加強(qiáng)重大動(dòng)物疫病防控延伸績(jī)效管理指標(biāo)體系》中將動(dòng)物疫病防控管理的績(jī)效分為疫病監(jiān)測(cè)、應(yīng)急處置、衛(wèi)生監(jiān)督管理、防治能力建設(shè)、協(xié)調(diào)經(jīng)費(fèi)落實(shí)情況等九方面。這些疫病防控延伸績(jī)效管理指標(biāo)直接指明鄉(xiāng)鎮(zhèn)畜牧獸醫(yī)站的主要職責(zé)和工作內(nèi)容。鄉(xiāng)鎮(zhèn)畜牧獸醫(yī)站的建設(shè)情況、技術(shù)水平狀況直接影響地區(qū)動(dòng)物防疫工作開(kāi)展情況。因動(dòng)物疫病防控管理影響因素眾多且相互關(guān)聯(lián)，相關(guān)數(shù)據(jù)收集限制極大。因此，通過(guò)鄉(xiāng)鎮(zhèn)畜牧獸醫(yī)建設(shè)情況和疫病防控技術(shù)水平評(píng)價(jià)地區(qū)動(dòng)物疫病防控管理水平。其中，用畜牧獸醫(yī)站建設(shè)比例衡量鄉(xiāng)鎮(zhèn)畜牧獸醫(yī)站建設(shè)情況，用技術(shù)人員比例衡量疫病防控技術(shù)水平。

（二）數(shù)據(jù)來(lái)源

被解釋變量由下列數(shù)據(jù)來(lái)源計(jì)算所得，解釋變量的數(shù)據(jù)主要來(lái)自2012～2016 年《獸醫(yī)公報(bào)》、2012～2016年《中國(guó)畜牧獸醫(yī)年鑒》和2012～2016年《中國(guó)統(tǒng)計(jì)年鑒》。其中，動(dòng)物疫病發(fā)生次數(shù)、發(fā)病動(dòng)物數(shù)、死亡動(dòng)物數(shù)、撲殺動(dòng)物數(shù)來(lái)自2012～2016年《獸醫(yī)公報(bào)》；存欄量、出欄量、動(dòng)物飼養(yǎng)規(guī)模（戶(hù)）數(shù)、鄉(xiāng)鎮(zhèn)畜牧獸醫(yī)站數(shù)量、技術(shù)人員數(shù)量等來(lái)自2012～2016年《中國(guó)畜牧獸醫(yī)年鑒》；森林覆蓋率、有效灌溉面積、廢水排放總量、地區(qū)年末人口數(shù)和鄉(xiāng)鎮(zhèn)級(jí)區(qū)劃數(shù)等來(lái)自2012～2016年《中國(guó)統(tǒng)計(jì)年鑒》。

（三）面板數(shù)據(jù)模型建立

根據(jù)動(dòng)物疫病分布的主趨勢(shì)方向（標(biāo)準(zhǔn)差橢圓方位角）Y與各個(gè)影響因素的單因素分析，建立被解釋變量Y與解釋變量森林覆蓋率（x1）、單位土地面積廢水排放量（x2）、動(dòng)物規(guī)?；B(yǎng)殖戶(hù)數(shù)比例（x3）、出欄率（x4）、人均動(dòng)物飼養(yǎng)量（x5）、畜牧獸醫(yī)站建設(shè)比例（x6）、技術(shù)人員比例(x7)的關(guān)系，構(gòu)建面板數(shù)據(jù)模型如下：

其中，e為誤差項(xiàng)。進(jìn)一步采用STATA軟件分析，各項(xiàng)回歸系數(shù)見(jiàn)表4。

表4 多元回歸分析

（四）回歸結(jié)果分析

從回歸模型分析結(jié)果看，按照各因素對(duì)動(dòng)物疫病分布的主趨勢(shì)方向影響程度大小依次為單位土地面積廢水排放量、畜牧獸醫(yī)站建設(shè)比例、森林覆蓋率、出欄率、技術(shù)人員比例、動(dòng)物規(guī)?；B(yǎng)殖戶(hù)數(shù)比例、人均動(dòng)物飼養(yǎng)量，且顯著性P 均小于0.1，說(shuō)明回歸方程變量之間線性關(guān)系顯著。其中單位土地面積廢水排放量、動(dòng)物規(guī)?；B(yǎng)殖戶(hù)數(shù)比例、出欄率、畜牧獸醫(yī)站建設(shè)比例與動(dòng)物疫病分布的標(biāo)準(zhǔn)差橢圓方位角成正方向變動(dòng)，即單位土地面積廢水排放量、動(dòng)物規(guī)?；B(yǎng)殖戶(hù)數(shù)比例、出欄率、畜牧獸醫(yī)站建設(shè)比例越大標(biāo)準(zhǔn)差橢圓方位角越大，動(dòng)物疫病傳播方向越向東南至西北方向傳播，而其他因素提高均會(huì)不同程度導(dǎo)致動(dòng)物疫病標(biāo)準(zhǔn)差橢圓方位角減小，即森林覆蓋率、人均動(dòng)物飼養(yǎng)量、技術(shù)人員比例越大動(dòng)物疫病越從東北至西南方向傳播。

五、結(jié)論及政策啟示

（一）結(jié)論

本文運(yùn)用數(shù)理統(tǒng)計(jì)和空間分析方法，選擇五種主要?jiǎng)游镆卟〔》N，從一般意義上探討中國(guó)動(dòng)物養(yǎng)殖業(yè)經(jīng)濟(jì)損失的空間格局演化，結(jié)論得出2012～2016年中國(guó)養(yǎng)殖業(yè)經(jīng)濟(jì)損失不相同，根據(jù)動(dòng)物疫病流行指數(shù)將全國(guó)31省區(qū)劃分為四個(gè)等級(jí)的風(fēng)險(xiǎn)區(qū)域，其中天津、江西、重慶等10省區(qū)為較嚴(yán)重地區(qū)，且損失嚴(yán)重地區(qū)通常連片聚集，還多發(fā)于邊界地區(qū)，呈現(xiàn)“邊界-核心”的結(jié)構(gòu)并日益明顯。中國(guó)養(yǎng)殖業(yè)經(jīng)濟(jì)損失的空間格局演化表現(xiàn)為其中心由中部向東南方向轉(zhuǎn)移，擴(kuò)散區(qū)域具有聯(lián)動(dòng)性。2012～2016 五年間經(jīng)濟(jì)損失嚴(yán)重區(qū)域的標(biāo)準(zhǔn)差橢圓覆蓋范圍整體由北向南遷移，呈現(xiàn)“擴(kuò)大-縮小-再擴(kuò)大”離散特征。中國(guó)養(yǎng)殖業(yè)經(jīng)濟(jì)損失的空間格局演化方向受地理生態(tài)環(huán)境、動(dòng)物飼養(yǎng)生產(chǎn)水平、動(dòng)物疫病防控管理影響。

（二）政策啟示

根據(jù)中國(guó)養(yǎng)殖業(yè)經(jīng)濟(jì)損失嚴(yán)重區(qū)域及其空間格局演化規(guī)律，對(duì)動(dòng)物疫病發(fā)生的極高風(fēng)險(xiǎn)區(qū)域、高風(fēng)險(xiǎn)區(qū)域、中等風(fēng)險(xiǎn)區(qū)域和低風(fēng)險(xiǎn)區(qū)域分別采用相對(duì)應(yīng)的防控措施。

（1）盡快制定和實(shí)施重大動(dòng)物疫病區(qū)域化撲滅計(jì)劃。視我國(guó)重大動(dòng)物疫病在各個(gè)區(qū)域的發(fā)生情況，進(jìn)行分區(qū)可撲滅性選擇，秉承“先局部、后全國(guó)，先重災(zāi)區(qū)，再輕微疫區(qū)”的原則，制定出重大動(dòng)物疫病分區(qū)撲滅計(jì)劃。該計(jì)劃首先要針對(duì)我國(guó)的極高風(fēng)險(xiǎn)和高風(fēng)險(xiǎn)區(qū)域，確認(rèn)撲滅計(jì)劃的整體目標(biāo)及分階段目標(biāo)，明確具體實(shí)施步驟和過(guò)程，確保撲滅經(jīng)費(fèi)到達(dá)重大動(dòng)物疫病疫區(qū)，在基礎(chǔ)設(shè)施和人員分配上一定要合理、充足，且工作人員要做好保護(hù)措施。

（2）進(jìn)一步加強(qiáng)重大動(dòng)物疫病區(qū)域性應(yīng)急反應(yīng)能力。要相對(duì)疫情進(jìn)行有效處理，制定一套重大動(dòng)物疫情應(yīng)急反應(yīng)機(jī)制必不可少，也是疫情處理的前期工作，因此要求政府要進(jìn)一步優(yōu)化突發(fā)重大動(dòng)物疫情應(yīng)急反應(yīng)的工作程序，制定出一套可行性應(yīng)急預(yù)案，提高應(yīng)急預(yù)案的區(qū)域針對(duì)性。另外，面對(duì)突發(fā)的重大動(dòng)物疫情應(yīng)建立重大動(dòng)物疫情的快速反應(yīng)體系，制定出可在短時(shí)間內(nèi)撲滅疫情的應(yīng)急反應(yīng)系統(tǒng)。為提高政府對(duì)于動(dòng)物疫情的應(yīng)急反應(yīng)能力，建議在國(guó)家、省、市、縣分層次建立動(dòng)物疫情應(yīng)急反應(yīng)工作小組，要求工作小組均為訓(xùn)練有素的動(dòng)物疫病專(zhuān)業(yè)人員，在小組中明確責(zé)任，做到權(quán)責(zé)分明。

（3）進(jìn)一步完善無(wú)規(guī)定動(dòng)物疫病區(qū)建設(shè)工作。根據(jù)前文對(duì)于重大動(dòng)物疫病的區(qū)劃，政府適當(dāng)調(diào)整無(wú)規(guī)定動(dòng)物疫病區(qū)域，在重大動(dòng)物疫區(qū)低風(fēng)險(xiǎn)區(qū)域建設(shè)無(wú)疫區(qū)，也可在隔離條件較好的中等風(fēng)險(xiǎn)區(qū)域適當(dāng)建立無(wú)疫區(qū)。在不適宜建設(shè)無(wú)疫區(qū)的高風(fēng)險(xiǎn)和極高風(fēng)險(xiǎn)區(qū)域或畜牧業(yè)大省可由一些大規(guī)模養(yǎng)殖戶(hù)建設(shè)相對(duì)安全的禽畜隔離區(qū)域，既可擴(kuò)大出口,也能為我國(guó)畜牧業(yè)養(yǎng)殖模式轉(zhuǎn)變提供示范作用。

（4）打破區(qū)域界限，改變主體單一治理狀況，利用省份間相互合作，充分調(diào)動(dòng)各區(qū)域的社會(huì)資源、最大限度地發(fā)揮市場(chǎng)及社會(huì)組織力量，提高動(dòng)物疫病防控效果。對(duì)于無(wú)疫區(qū)的地方政府，要求向省中央提供真實(shí)、準(zhǔn)確的動(dòng)物疫病流行病學(xué)數(shù)據(jù)，以便中央整體統(tǒng)一監(jiān)測(cè)，大力推進(jìn)我國(guó)無(wú)疫區(qū)的具體工作，為其他疫病區(qū)域治理提供借鑒。

（5）進(jìn)一步完善重大動(dòng)物疫病區(qū)域防控法規(guī)等，根據(jù)《中華人民共和國(guó)動(dòng)物防疫法》，我國(guó)政府需要制定和完善疫情防控工作的具體過(guò)程，其中包括：疫情監(jiān)測(cè)、疫情公布、檢疫監(jiān)督、獸醫(yī)衛(wèi)生管理、疫情追溯、市場(chǎng)準(zhǔn)入、區(qū)域化管理、應(yīng)急反應(yīng)、撲殺賠償?shù)确矫娴囊?guī)章制度。同時(shí)，加大對(duì)違犯動(dòng)物疫病檢疫、不遵守強(qiáng)制免疫、疫情隱瞞報(bào)告等行為的處罰力度，提高動(dòng)物防疫執(zhí)法權(quán)威。從根本上控制動(dòng)物疫病發(fā)生和傳播，從而最大限度降低中國(guó)養(yǎng)殖業(yè)經(jīng)濟(jì)損失。