中國動物疫病經濟損失的空間格局演化研究*

李燕凌，王繼紅，蘇青松

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

一、引 言

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

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

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

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

二、中國動物疫病經濟損失的區(qū)域差異

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

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

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

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

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

續(xù)表1

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

三、中國動物疫病經濟損失的空間分布分析

（一）動物疫病經濟損失的空間分布分析方法

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

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

2.SDE主要參數(shù)

計算公式如下。

3.標準差橢圓分析結果

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

（二）2012～2016年主要動物疫病流行數(shù)據(jù)及SDE模型建構

1.主要動物疫病流行數(shù)據(jù)

根據(jù)農業(yè)部《獸醫(yī)公報》數(shù)據(jù)，統(tǒng)計2012～2016年全國31省區(qū)禽流感、新城疫、豬藍耳、豬瘟及口蹄疫感染病例數(shù)，制作五種主要動物疫病流行數(shù)據(jù)。

2.主要動物疫病SDE計算結果

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

表2 2012～2016年五種疫病標準差橢圓分析計算結果

（三）2012～2016年主要動物疫病SDE模型結果分析

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

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

通過確定2012～2016年禽流感疫情發(fā)生的中心城市與每年方位角轉變得出分析結果：五年來禽流感疫情多覆蓋在東南地區(qū)，傳播方向為“自西北向東南再向西北”回旋趨勢，2014年有向東北和偏南部傳播趨勢。疫情覆蓋區(qū)域呈“縮小—擴大—再縮小—再擴大”反反復復模式，疫情控制稍有疏忽便會呈現(xiàn)回溫之勢?？梢娂词挂咔榈玫娇刂迫圆荒苄傅　?/p>

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

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

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

3.豬藍耳疫情空間分布格局分析

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

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

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

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

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

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

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

四、中國動物養(yǎng)殖業(yè)經濟損失差異的影響因素

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

（一）變量定義

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

表3 動物疫病流行風險區(qū)域差異的影響因素

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

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

2.動物飼養(yǎng)生產水平

本文主要通過三個指標評價地區(qū)動物飼養(yǎng)生產水平，研究各指標對動物疫病流行風險的影響。一是規(guī)?；B(yǎng)殖程度，用動物規(guī)?；B(yǎng)殖戶數(shù)比來衡量，比例越高規(guī)?；B(yǎng)殖程度越高；二是動物生產力水平，由出欄率衡量，出欄率越高動物生產水平越高；三是養(yǎng)殖密度，由人均動物飼養(yǎng)量衡量，人均動物飼養(yǎng)量越高養(yǎng)殖密度越高。

3.動物疫病防控管理

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

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

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

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

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

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

表4 多元回歸分析

（四）回歸結果分析

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

五、結論及政策啟示

（一）結論

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

（二）政策啟示

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

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

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

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

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

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