中國市域生豬生產(chǎn)時(shí)空演變特征及影響因素

王如如，房艷剛,2,3*，劉建志

（1. 東北師范大學(xué)地理科學(xué)學(xué)院，吉林 長春 130024；2. 東北師范大學(xué)長白山地理過程與生態(tài)安全教育部重點(diǎn)實(shí)驗(yàn)室，

吉林 長春 130024；3. 東北師范大學(xué)中國東北研究院，吉林 長春 130024）

二戰(zhàn)以來，在快速人口城鎮(zhèn)化背景下，眾多新興發(fā)展中國家城鄉(xiāng)居民對(duì)于動(dòng)物源性食品（ADF）的需求不斷增加，從而推動(dòng)了牲畜產(chǎn)量的增長，這一過程被稱為“畜牧革命”[1]。“畜牧革命”在中國表現(xiàn)最為顯著[2]。中國人均肉類消費(fèi)量從1978年的8.90 kg增長到2017年的62.26 kg，總增長率高達(dá)599.55%。肉類消費(fèi)的上升倒逼農(nóng)業(yè)結(jié)構(gòu)轉(zhuǎn)型，構(gòu)成了中國隱形農(nóng)業(yè)革命的一部分[3]。從肉類消費(fèi)結(jié)構(gòu)上看，豬肉在中國肉類消費(fèi)中居主導(dǎo)地位（豬肉占比60%以上）[4-5]，生豬市場潛力巨大。然而，現(xiàn)階段中國生豬生產(chǎn)正面臨著規(guī)?；c環(huán)保和疫病防治的矛盾、貿(mào)易政策變動(dòng)與供給波動(dòng)的復(fù)雜關(guān)聯(lián)。一方面，隨著生豬規(guī)模化生產(chǎn)的快速發(fā)展，環(huán)境污染問題日益突出[6]，疫病風(fēng)險(xiǎn)引發(fā)的產(chǎn)業(yè)脆弱性持續(xù)加大，土地利用規(guī)劃背景下的土地資源脅迫問題愈加嚴(yán)重[7]；另一方面，我國生豬綜合競爭力明顯低于發(fā)達(dá)國家[8]，國外生豬產(chǎn)品的大量涌入給國內(nèi)生豬生產(chǎn)帶來了巨大沖擊。新形勢下，中美貿(mào)易摩擦對(duì)大豆等蛋白飼料進(jìn)口量的限制進(jìn)一步提高了生豬生產(chǎn)成本。2019年以來，受“豬周期”下行、非洲豬瘟疫情沖擊、地方環(huán)保限養(yǎng)等多因素疊加影響，我國生豬產(chǎn)能明顯下降，生豬價(jià)格大幅度上漲[9]。因此，在消費(fèi)需求變化、可持續(xù)發(fā)展及中美貿(mào)易摩擦等背景下，基于生豬生產(chǎn)的時(shí)空格局演變特征探索我國生豬產(chǎn)業(yè)的空間布局優(yōu)化策略，對(duì)穩(wěn)定生豬市場供給、促進(jìn)生豬產(chǎn)業(yè)可持續(xù)發(fā)展、提升農(nóng)業(yè)綜合競爭力具有重要意義。

國內(nèi)外關(guān)于生豬生產(chǎn)的相關(guān)研究主要涉及生豬生產(chǎn)模式[10]、環(huán)境規(guī)制[11]、生豬質(zhì)量安全[12-13]、市場波動(dòng)[14-15]、空間布局[16-17]及其影響因素[18]等。對(duì)于生豬生產(chǎn)空間布局的研究，學(xué)者多采用基尼系數(shù)、地域集中度指數(shù)、區(qū)域重心分析法、空間自相關(guān)等來揭示生豬生產(chǎn)的區(qū)域差異[19-20]，研究尺度涉及國家、區(qū)域、省域單元。在空間布局方面，趙俊偉等[21]認(rèn)為生豬生產(chǎn)的空間集聚具有一定的階段性特征，時(shí)空分布呈現(xiàn)“集-散-集”的變化特征；劉勍等[19]認(rèn)為生豬生產(chǎn)區(qū)域布局呈現(xiàn)向北向西移動(dòng)的趨勢；但王歡和喬娟[17]認(rèn)為東南地區(qū)仍是生豬生產(chǎn)優(yōu)勢區(qū)，生豬生產(chǎn)區(qū)域布局呈現(xiàn)向東向北移動(dòng)的趨勢。同時(shí)，還有學(xué)者指出隨著規(guī)?；陌l(fā)展，生豬生產(chǎn)區(qū)域?qū)⑾螯S淮海平原和東北地區(qū)等糧食主產(chǎn)區(qū)以及其他邊遠(yuǎn)地區(qū)轉(zhuǎn)移[22]。在影響因素方面，大部分學(xué)者提出生豬生產(chǎn)受相鄰地區(qū)和本地區(qū)因素的共同影響[17,21]。此外，還有學(xué)者通過對(duì)生豬產(chǎn)地移動(dòng)的分析發(fā)現(xiàn)生豬生產(chǎn)的自然性布局條件減弱，經(jīng)濟(jì)性布局條件不斷增強(qiáng)[16,23]。總的來看，目前我國生豬生產(chǎn)時(shí)空演變特征及影響因素的研究成果豐碩。但是，此類研究大多基于省域尺度，對(duì)市域尺度生豬生產(chǎn)的研究鮮見。由于各省域內(nèi)部自然和社會(huì)經(jīng)濟(jì)條件差異明顯，省域尺度研究對(duì)于生豬生產(chǎn)時(shí)空演變規(guī)律和驅(qū)動(dòng)機(jī)制的刻畫相對(duì)不夠細(xì)致精準(zhǔn)。在研究內(nèi)容上，多數(shù)研究集中在對(duì)生豬平均生產(chǎn)水平影響因素的探究，忽略了不同生豬生產(chǎn)水平影響因素的差異。此外，此類研究對(duì)于生豬生產(chǎn)空間布局優(yōu)化的探討相對(duì)較少。

鑒于此，本文以市域?yàn)榛狙芯繂卧?，首先運(yùn)用有序樣品聚類法探尋中國生豬生產(chǎn)發(fā)展演化的關(guān)鍵節(jié)點(diǎn)，確定研究截面；其次，基于標(biāo)準(zhǔn)差橢圓和探索性空間數(shù)據(jù)分析等空間分析方法探索其時(shí)空演變特征；最后利用分位數(shù)回歸模型對(duì)時(shí)空格局演化機(jī)制進(jìn)行定量剖析，以期為我國生豬養(yǎng)殖業(yè)的空間布局優(yōu)化及制定生豬養(yǎng)殖業(yè)可持續(xù)發(fā)展政策提供科學(xué)參考。

1 研究方法與數(shù)據(jù)來源

1.1 研究方法

1.1.1 有序樣品聚類法 有序樣品聚類法又稱最優(yōu)分割法，是針對(duì)有序樣本的聚類分析方法。該方法在分類時(shí)不打破樣本順序[24]，基于類內(nèi)離差平方和(W)最小且類間離差平方和(B)最大的原則，選取最優(yōu)分割點(diǎn)，對(duì)連續(xù)樣本數(shù)據(jù)劃分階段，從而達(dá)到分類的目的。

式中：n表示樣本數(shù)量；L表示包含n個(gè)樣本的序列被劃分的類別數(shù)目；nk表示第k類含有的樣品數(shù)量；Xkj表示第k類第j個(gè)樣本對(duì)應(yīng)的指標(biāo)；表示第k類指標(biāo)的均值；為所有樣本的均值。

1.1.2 標(biāo)準(zhǔn)差橢圓 標(biāo)準(zhǔn)差橢圓是通過以中心、長軸、短軸、方位角為基本參數(shù)定量描述地理要素的中心性、展布性、方向性、空間形態(tài)等，從而揭示地理要素空間分布特征的方法[25]。本文基于ArcGIS 10.2軟件構(gòu)建生豬出欄密度標(biāo)準(zhǔn)差橢圓，分析其空間分布特征。

1.1.3 分位數(shù)回歸 分位數(shù)回歸用于研究自變量與因變量條件分位數(shù)之間的關(guān)系，使用殘差絕對(duì)值的加權(quán)平均作為最小化目標(biāo)函數(shù)，在消除極端值影響的同時(shí)還能提供因變量條件分布的完整特性[26]。相比OLS，分位數(shù)回歸具有兩大優(yōu)勢：一是對(duì)離群值與偏態(tài)分布的不敏感性；二是對(duì)于因變量條件分布完整特征的反映[27]。本研究中，采用分位數(shù)回歸方法對(duì)生豬出欄密度進(jìn)行回歸估計(jì)以反映不同生豬出欄水平影響因素的差異。具體的計(jì)量模型如下：

式中：Y表示生豬出欄密度；X表示自變量；i表示市域；m表示自變量的個(gè)數(shù)；μ表示隨機(jī)誤差項(xiàng)；β0表示常數(shù)項(xiàng)，βi為待估計(jì)的解釋變量的系數(shù)。

分位數(shù)回歸系數(shù)的估計(jì)通過使殘差絕對(duì)值的加權(quán)平均最小化得出，即：

式中：yi代表因變量的第i個(gè)觀測值；xi代表自變量的第i個(gè)觀測值；θ是估計(jì)時(shí)人為設(shè)定的分位數(shù)值，本研究設(shè)定為10，25，50，75，90；β是對(duì)應(yīng)x的偏回歸系數(shù)。將(3)式設(shè)定的模型代入(4)式的模型可估算出y相應(yīng)分位數(shù)上的回歸系數(shù)。根據(jù)模型需要保證數(shù)據(jù)的穩(wěn)定性，本文對(duì)數(shù)值大于1的因子進(jìn)行對(duì)數(shù)變換[28]。

1.2 數(shù)據(jù)來源與指標(biāo)選擇

首先運(yùn)用有序樣品聚類法選取研究截面。本文借助Matlab 9.1軟件平臺(tái)，采用有序樣品聚類法對(duì)1978—2017年中國生豬出欄密度進(jìn)行分析，首先得到1988、1994、2000、2008年四個(gè)最優(yōu)分割點(diǎn)。囿于數(shù)據(jù)的可獲取性，本文不考慮1988、1994年。此外，相關(guān)研究表明，2000年是我國生豬出欄量增長率由下降到穩(wěn)步提升的重要節(jié)點(diǎn)，2008年是增長率由穩(wěn)步提升到波動(dòng)下降的重要節(jié)點(diǎn)[21]。因此，本文選取2000、2008和2017年作為研究截面。

本文所用統(tǒng)計(jì)數(shù)據(jù)主要來源于相應(yīng)年份的《中國統(tǒng)計(jì)年鑒》與各省統(tǒng)計(jì)年鑒。以上兩個(gè)年鑒缺失的數(shù)據(jù)獲取于各地級(jí)市統(tǒng)計(jì)年鑒與統(tǒng)計(jì)公報(bào)。生豬生產(chǎn)空間格局的研究范圍為中國大陸315個(gè)市域行政單元，按2016年最新市域單元?jiǎng)澐智闆r對(duì)部分市域數(shù)據(jù)進(jìn)行歸并處理。數(shù)據(jù)分析過程中，對(duì)出入較大的數(shù)據(jù)做了剔除處理，對(duì)異常值或缺失數(shù)據(jù)的市域單元，用插值或趨勢外推法進(jìn)行估計(jì)。對(duì)于影響因素研究，由于2000年自變量數(shù)據(jù)的可獲取性較差，本文僅選取2008、2017年兩個(gè)時(shí)間截面，并將前述范圍缺失數(shù)據(jù)的市域剔除，最終選取了297個(gè)市域行政單元。氣溫、高程與耕地?cái)?shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http://www.resdc.cn/）；限制開發(fā)區(qū)名錄來源于各省主體功能區(qū)規(guī)劃。

本文使用的生豬出欄密度是地區(qū)生豬出欄總量除以其行政區(qū)面積。生豬生產(chǎn)時(shí)空格局演變是區(qū)域本底因素和內(nèi)外部驅(qū)動(dòng)因素交互作用、相互影響的結(jié)果。依據(jù)指標(biāo)選擇的科學(xué)性和可得性原則，本文構(gòu)建如表1的影響因素測度指標(biāo)體系。

表1 生豬生產(chǎn)時(shí)空格局影響因素測度指標(biāo)體系Table 1 Measurement index system of factors influencing the spatio-temporal pattern of hog production

具體指標(biāo)的解釋如下：1）區(qū)域本底因素。①自然資源稟賦是農(nóng)業(yè)區(qū)域形成發(fā)展的自然基礎(chǔ)，對(duì)于生豬養(yǎng)殖而言主要包括飼料資源、本底資源、自然條件[29]。飼料資源包括地均玉米產(chǎn)量（Core）、地均大豆產(chǎn)量（Soy），對(duì)生豬這種耗糧型產(chǎn)業(yè)的影響較大[17]；本底資源包括生豬養(yǎng)殖業(yè)所需的水土資源，以人均耕地面積（Land）、人均水資源量（Water）進(jìn)行表征，人均耕地面積越大，既能夠滿足生豬生產(chǎn)規(guī)?；瘜?duì)土地資源的需求[21]，也有利于負(fù)荷更多的畜禽污染物[20]，人均水資源量越多對(duì)生豬這一耗水型產(chǎn)業(yè)的支撐作用越強(qiáng)；自然條件選取年平均氣溫（Temp）、高程（Elev），一般來講，過冷和過熱都會(huì)提高生豬養(yǎng)殖成本，高程表征的地形地勢條件將會(huì)影響生豬的規(guī)模化生產(chǎn)。②區(qū)位條件能夠影響生豬運(yùn)輸、交易成本以及資金、技術(shù)、信息等現(xiàn)代生產(chǎn)要素的傳遞[29]，采用單位面積公路里程（Road）進(jìn)行表征。③環(huán)境條件主要體現(xiàn)在區(qū)域主體功能對(duì)生豬生產(chǎn)的約束。本文選取區(qū)域內(nèi)限制開發(fā)區(qū)面積占比（Rest）來衡量區(qū)域環(huán)境容納能力。

2）內(nèi)外部驅(qū)動(dòng)因素。①經(jīng)濟(jì)發(fā)展水平對(duì)生豬生產(chǎn)的影響復(fù)雜。一方面，經(jīng)濟(jì)發(fā)展水平的提升能夠?yàn)樯i養(yǎng)殖業(yè)提供更多的資金與技術(shù)支持，另一方面，當(dāng)經(jīng)濟(jì)發(fā)展水平達(dá)到一定高度后，產(chǎn)業(yè)結(jié)構(gòu)升級(jí)對(duì)生豬生產(chǎn)具有排斥和置換作用。因此，兩者可能呈倒“U”關(guān)系。本文選取人均GDP（Agdp）表征經(jīng)濟(jì)發(fā)展水平。②區(qū)域勞動(dòng)力、土地等生產(chǎn)要素快速非農(nóng)化會(huì)壓縮生豬生產(chǎn)的要素投入。本文選擇鄉(xiāng)村非農(nóng)就業(yè)率（Nona）表征鄉(xiāng)村非農(nóng)化水平。③市場需求一定程度上取決于消費(fèi)者數(shù)量與城鄉(xiāng)居民收入水平，對(duì)生豬生產(chǎn)發(fā)展具有推動(dòng)作用。本文選取區(qū)域年末總?cè)丝冢≒opu）及城鄉(xiāng)居民人均可支配收入（Inco）進(jìn)行表征。

2 中國市域生豬出欄密度的時(shí)空格局及演化

2.1 中國市域生豬出欄密度的時(shí)序演變特征

2000—2017 年中國生豬出欄密度總體呈現(xiàn)快速增長的趨勢，總增長率達(dá)到32.46%。2000—2008年，生豬出欄密度由0.65頭/hm2增加到0.81頭/hm2，增幅達(dá)24.48%；2008—2017年，生豬出欄密度由0.81頭/hm2增加到0.86頭/hm2，增幅為6.41%。這表明中國生豬出欄密度在2000—2008年增長速度較快，2008—2017年增長速度趨緩。從核密度曲線來看（圖1），2000—2017年生豬出欄低密度區(qū)所占比重明顯下降，整體縮小了中國生豬生產(chǎn)的地區(qū)差距。這一現(xiàn)象通過泰爾指數(shù)的對(duì)比更加明顯，2000—2017年中國市域生豬出欄密度的泰爾指數(shù)由0.417下降到0.409，再到0.407，表明中國生豬生產(chǎn)的地區(qū)差異不斷縮小，“追趕效應(yīng)”顯現(xiàn)。

圖1 2000—2017年全國生豬出欄密度的核密度曲線Fig. 1 Kernel density curve of the density of hog slaughterhouses of China from 2000 to 2017

2.2 中國市域生豬出欄密度的時(shí)空演變特征

2000—2017 年中國市域生豬出欄密度區(qū)域差異顯著，以“胡煥庸線”為界呈現(xiàn)出西低東高的空間格局。四川盆地、湘江流域、黃淮海平原等糧食主產(chǎn)區(qū)和哈大沿線、珠江流域下游等經(jīng)濟(jì)較發(fā)達(dá)區(qū)為生豬生產(chǎn)高值區(qū)（>2.21頭/hm2，圖2）。

2000—2008 年中國各地區(qū)生豬出欄密度普遍增長（圖3），增長區(qū)域占比達(dá)75.56%。增長幅度較大的區(qū)域主要分布在東北地區(qū)、秦巴山區(qū)、川西高原、云南高原和橫斷山區(qū)等，此類地區(qū)基期生豬出欄密度普遍較小，對(duì)整體的生豬出欄分布格局未造成較大影響；負(fù)增長區(qū)域散布在河北、寧夏、廣西中東部、豫皖蘇交界和蘇南等地區(qū)。該時(shí)段生豬生產(chǎn)分布重心由湖北省孝感市向西北方向移動(dòng)36.21 km到達(dá)湖北省隨州市南部，生豬生產(chǎn)“北進(jìn)西移”現(xiàn)象顯現(xiàn)（圖4）。2008—2017年，全國生豬出欄密度變化率區(qū)域差異顯著，增長區(qū)域占比達(dá)56.83%，主要分布在黃淮海平原，黃土高原，珠江流域中上游以及西藏、青海與甘肅等地，而減少區(qū)域主要分布在東北地區(qū)，長江流域的四川、重慶、湖南、長江三角洲，東南沿海以及部分大城市（如北京、天津、鄭州等地，圖3）。該時(shí)段生豬生產(chǎn)分布重心由隨州市南部向東北方向偏移了13.54 km。2000—2017年，生豬生產(chǎn)分布重心總體上向北偏西方向移動(dòng)39.10 km，且“北進(jìn)”速度快于“西移”，呈現(xiàn)出“南豬北移”的演變格局（圖4）。

圖2 2000—2017年中國市域生豬出欄密度時(shí)空格局圖Fig. 2 Spatio-temporal pattern of the density of hog slaughterhouses in China’s cities during 2000—2017

圖3 2000—2017年中國市域生豬出欄密度變化率圖Fig. 3 The change rate of the density of hog slaughterhouses in China’s cities during 2000—2017

分地區(qū)來看，2000—2017年東北地區(qū)生豬出欄密度偏低（2017年僅達(dá)到0.78頭/hm2），呈現(xiàn)總體增長但先升后降的態(tài)勢，2000—2017年總體增幅為56.78%，其中2000—2008年上升89.61%，2008—2017年下降17.32%；黃淮海平原生豬出欄密度較高（2017年達(dá)到3.43頭/hm2），2000—2017年增幅達(dá)37.55%，除一些特大城市（如北京、天津）外，其余地區(qū)生豬出欄密度呈持續(xù)增長態(tài)勢；長江流域生豬出欄密度較高，2000—2017年整體呈現(xiàn)先升后降的演化趨勢，其中上游地區(qū)表現(xiàn)出較大幅度的先升后降趨勢，中游地區(qū)平穩(wěn)增加（[15%，35%）），下游地區(qū)的皖南山區(qū)增長率持續(xù)增大，而長江三角洲地區(qū)基本呈持續(xù)下降趨勢；湘江流域平均生豬出欄密度極高（2017年達(dá)到了3.46頭/hm2），但在2008—2017年明顯下降，降幅達(dá)20.64%；珠江流域上游地區(qū)生豬出欄密度處于中等水平（2017年為1.58頭/hm2），但始終保持較強(qiáng)的增長態(tài)勢，2000—2017年增長率高達(dá)81.10%；東南沿海地區(qū)生豬出欄密度處于中等水平（2017年達(dá)到1.59頭/hm2），但2008—2017年出現(xiàn)下降趨勢，降幅達(dá)9.12%。

圖4 2000—2017年中國市域生豬出欄密度標(biāo)準(zhǔn)差橢圓曲線圖Fig. 4 The standard deviation ellipse of the density of hog slaughterhouses in China’s cities during 2000—2017

2.3 中國市域生豬出欄密度空間關(guān)聯(lián)特征

2000—2017 年，中國生豬出欄密度全局空間自相關(guān)Moran's I值都通過了1%水平下的顯著性檢驗(yàn)，且全為正值，表明中國市域生豬出欄密度在空間分布上呈現(xiàn)顯著的空間集聚現(xiàn)象。2000—2017年，Moran’s I值下降態(tài)勢明顯，由2000年的0.56下降到2008年的0.52和2017年的0.37，一定程度上反映出自然本底、距離等地理因素對(duì)生豬生產(chǎn)分布的作用弱化。

進(jìn)一步對(duì)中國市域生豬出欄密度局部空間自相關(guān)進(jìn)行測度，結(jié)果顯示（圖5）：2000—2017年高高集聚區(qū)的市域主要分布在四川盆地、湘江流域和黃淮海平原。2000—2008年三大高高集聚區(qū)的數(shù)量沒有發(fā)生明顯變化。2008—2017年，四川盆地、湘江流域的生豬出欄密度高高集聚區(qū)收縮，兩者的數(shù)量分別由11個(gè)減少到8個(gè)，7個(gè)減少到5個(gè)，而黃淮海平原高高集聚區(qū)的數(shù)量由9個(gè)大幅度擴(kuò)大到24個(gè)。2000—2017年，低低集聚區(qū)的市域主要分布在西藏、青海、寧夏、甘肅、黑龍江、內(nèi)蒙古大部以及吉林東部山區(qū)；高—低和低—高類型區(qū)數(shù)量較少，“高地現(xiàn)象”和“洼地現(xiàn)象”都不明顯。

3 中國市域生豬出欄密度時(shí)空格局演化的影響因素

本文對(duì)影響因子進(jìn)行方差膨脹因子檢驗(yàn)（VIF）以消除多重共線性的影響，選取VIF值小于5的最佳影響因素組合進(jìn)行基于bootstrap方法的分位數(shù)回歸（QR）?？v向來看，2008與2017年兩個(gè)時(shí)間截面通過顯著性檢驗(yàn)的影響因子及其影響強(qiáng)度存在差異；橫向來看，不同出欄水平影響因子的影響強(qiáng)度也存在差異，具體如下（表2）。

1）飼料資源。①地均玉米產(chǎn)量對(duì)生豬出欄密度起顯著正向影響。隨著分位數(shù)的提高，回歸系數(shù)不斷減小，地均玉米產(chǎn)量對(duì)生豬出欄密度的正向促進(jìn)作用減弱。這主要是由于生豬出欄密度高值區(qū)大部分位于糧食主產(chǎn)區(qū)，地均玉米產(chǎn)量較大（如黃淮海平原）；部分生豬出欄密度高值區(qū)也可以通過域外玉米的調(diào)入彌補(bǔ)本地玉米等飼料資源的不足。②2008、2017年地均大豆產(chǎn)量對(duì)生豬出欄密度起顯著正向影響，回歸系數(shù)隨條件分位數(shù)的增大而減?。ㄈ?017年由第10分位數(shù)的0.223減小到第90分位數(shù)的0.141），因此，地均大豆產(chǎn)量同樣對(duì)生豬出欄密度較低地區(qū)的促進(jìn)作用較強(qiáng)。2017年地均大豆產(chǎn)量的正向影響程度有所減弱，主要由于近年來國外進(jìn)口大豆對(duì)國內(nèi)的補(bǔ)充（2017年大豆進(jìn)口量是2008年的2.55倍），生豬生產(chǎn)對(duì)國內(nèi)大豆的依賴程度減弱。當(dāng)前，中美貿(mào)易摩擦使大豆等蛋白飼料的進(jìn)口量短期內(nèi)大量減少，對(duì)生豬養(yǎng)殖業(yè)飼料供給產(chǎn)生較大影響。

圖5 2000—2017年中國市域生豬出欄密度LISA圖Fig. 5 LISA cluster diagrams of the density of hog slaughterhouses in China’s cities during 2000—2017

表2 QR模型結(jié)果匯總表Table 2 Result summary table of QR model

2）水土資源。①人均耕地面積與生豬出欄密度呈負(fù)相關(guān)關(guān)系，與預(yù)期不符。這主要因?yàn)樵谌司孛娣e較小的地區(qū)，農(nóng)戶在種植業(yè)中獲得的收入較低，傾向于選擇種植和養(yǎng)殖相結(jié)合的方式提高收入，從而促進(jìn)生豬養(yǎng)殖業(yè)的發(fā)展。橫向上看，隨著分位數(shù)的提高，回歸系數(shù)的絕對(duì)值不斷變大?？v向上看，2017年各分位數(shù)上回歸系數(shù)的絕對(duì)值明顯大于2008年，負(fù)相關(guān)程度增強(qiáng)，主要是因?yàn)樵?017年，隨著種植業(yè)比較收益的進(jìn)一步下降，部分人均耕地面積較小的地區(qū)發(fā)展生豬養(yǎng)殖業(yè)的意愿更大，生豬養(yǎng)殖大幅度上升（如黃淮海平原）。②2008年人均水資源量的回歸系數(shù)顯著為正，符合預(yù)期方向，但2017年回歸系數(shù)不顯著。2008—2017年由于環(huán)境規(guī)制、農(nóng)業(yè)政策等因素影響，長江流域、東北地區(qū)等人均水資源量豐富區(qū)生豬養(yǎng)殖大幅減少，同時(shí)受南水北調(diào)、集約化生產(chǎn)等因素影響，黃淮海平原等缺水區(qū)生豬產(chǎn)業(yè)供水壓力有所緩解，養(yǎng)殖量也不斷增加。

3）自然條件。①年平均氣溫在各分位數(shù)上對(duì)生豬出欄密度呈顯著正向影響。隨著分位數(shù)的提高，年平均氣溫的回歸系數(shù)基本下降（如2017年由第25分位數(shù)的0.110下降到第90分位數(shù)的0.053），表明其正向促進(jìn)作用隨著生豬出欄密度的增加而變小。這主要因?yàn)樯i出欄密度較高的區(qū)域大部分位于氣溫較適宜的地區(qū)（基于生豬出欄密度的空間格局可知）。②高程對(duì)生豬出欄密度呈負(fù)向影響。縱向上看，2008年各分位數(shù)上的回歸系數(shù)不顯著，而2017年回歸系數(shù)基本顯著，負(fù)向影響開始顯現(xiàn)；橫向上看，2017年回歸系數(shù)的絕對(duì)值隨分位數(shù)的提高呈增長趨勢（如由第25分位數(shù)的-0.067減少到第90分位數(shù)的-0.161）。2017年生豬生產(chǎn)規(guī)模化水平高于2008年（2017年出欄量為500頭以上的生豬養(yǎng)殖場占出欄量為50頭以上的比重是2008年的1.56倍），中、高分位數(shù)規(guī)模化水平高于低分位數(shù)（2017年生豬生產(chǎn)高值區(qū)年出欄500頭以上的規(guī)模養(yǎng)殖場占全國的60%以上），因此，隨著規(guī)?；潭鹊募由?，高程的負(fù)向影響凸顯。

4）環(huán)境條件。限制開發(fā)區(qū)面積占比與生豬出欄密度呈負(fù)相關(guān)關(guān)系，與預(yù)期相符。橫向上看，回歸系數(shù)隨條件分位數(shù)的增大而減小?？v向上看，2017年各分位數(shù)上回歸系數(shù)的絕對(duì)值小于2008年，負(fù)向影響程度有所減弱，這主要是因?yàn)?008—2017年長江流域、東北平原等限制開發(fā)區(qū)面積占比較小地區(qū)的生豬出欄密度不斷減小。2008年之后相繼頒布的《水污染防治行動(dòng)計(jì)劃》、《農(nóng)業(yè)部關(guān)于促進(jìn)南方水網(wǎng)地區(qū)生豬養(yǎng)殖布局調(diào)整優(yōu)化的指導(dǎo)意見》等使長江流域環(huán)境規(guī)制壓力加大，成為其生豬出欄密度下降的重要原因。

5）區(qū)位條件。單位面積公路里程與生豬出欄密度呈顯著正相關(guān)關(guān)系?？v向上看，2017年各分位數(shù)上的回歸系數(shù)小于2008年，正向影響程度有所減弱，主要是由于2017年我國交通基礎(chǔ)設(shè)施明顯提升，對(duì)生豬生產(chǎn)的促進(jìn)作用弱化。橫向上看，2017年回歸系數(shù)隨條件分位數(shù)的提高而增加，主要是因?yàn)樯i出欄密度高值區(qū)生豬商品率較高，對(duì)交通運(yùn)輸條件的要求較高。

6）經(jīng)濟(jì)規(guī)模。整體來看，人均GDP與生豬出欄密度之間呈倒“U”關(guān)系，符合預(yù)期。具體而言，2008年隨著分位數(shù)水平的提升，人均GDP的正向促進(jìn)作用逐漸減弱（回歸系數(shù)由第10分位數(shù)的0.373下降到第90分位數(shù)的0.040），而2017年人均GDP的回歸系數(shù)不顯著且基本為負(fù)。2017年人均GDP普遍高于2008年，高分位數(shù)的人均GDP普遍高于低分位數(shù)。因此，在經(jīng)濟(jì)發(fā)展水平相對(duì)較低的時(shí)期和地區(qū)，人均GDP的提高可以促進(jìn)當(dāng)?shù)厣i生產(chǎn)，當(dāng)人均GDP達(dá)到一定水平后，其正向作用不斷減小，而負(fù)向作用開始顯現(xiàn)，如在東南沿海地區(qū)，經(jīng)濟(jì)發(fā)展引發(fā)的產(chǎn)業(yè)轉(zhuǎn)型升級(jí)使生豬養(yǎng)殖業(yè)不斷被置換出去。

7）鄉(xiāng)村非農(nóng)化水平。如表2所示，隨著分位數(shù)的提高，2008年鄉(xiāng)村非農(nóng)化水平的回歸系數(shù)顯著為正（第10～75分位數(shù)），且不斷降低，而2017年回歸系數(shù)顯著為負(fù)（第25～90分位數(shù)），且絕對(duì)值不斷變大，因此鄉(xiāng)村非農(nóng)化水平對(duì)生豬生產(chǎn)的影響也具有倒“U”趨勢。在非農(nóng)化水平較低時(shí)，當(dāng)?shù)胤寝r(nóng)產(chǎn)業(yè)的發(fā)展具有反哺功能，可為生豬生產(chǎn)提供資金、技術(shù)等多方面支持，促進(jìn)其生產(chǎn)經(jīng)營；隨著城鎮(zhèn)化和工業(yè)化進(jìn)程加快，非農(nóng)產(chǎn)業(yè)的發(fā)展會(huì)壓縮生豬養(yǎng)殖業(yè)勞動(dòng)力、土地等要素投入，抑制生豬養(yǎng)殖業(yè)的發(fā)展，如四川、湖南等地農(nóng)業(yè)勞動(dòng)力大量轉(zhuǎn)移[31]，生豬養(yǎng)殖量也呈明顯下降趨勢。

8）市場需求。區(qū)域年末總?cè)丝趯?duì)生豬出欄密度呈正向顯著影響，且這種正向影響在中、低生豬出欄水平體現(xiàn)的較為明顯（如2008年的第10～50分位數(shù)），這主要是因?yàn)樯i出欄密度高的地區(qū)生豬商品率較高，更多面向域外市場。結(jié)果顯示，城鄉(xiāng)居民人均可支配收入在各分位數(shù)上基本未通過顯著性檢驗(yàn)，可能是因?yàn)槌青l(xiāng)居民人均可支配收入增加對(duì)豬肉的需求作用具有階段性、本地生豬市場的擴(kuò)大也可以通過其他區(qū)域得到滿足、收入增加對(duì)當(dāng)?shù)剞r(nóng)戶生豬生產(chǎn)意愿的作用復(fù)雜等。

4 結(jié)論與討論

4.1 結(jié)論

1）2000—2017 年，中國市域生豬出欄密度顯著增長，總增長率達(dá)32.46%。四川盆地、湘江流域、黃淮海平原等糧食主產(chǎn)區(qū)和哈大沿線、珠江流域下游等部分經(jīng)濟(jì)較發(fā)達(dá)區(qū)為生豬生產(chǎn)高值區(qū)。2008—2017年全國生豬出欄密度變化率區(qū)域差異顯著，增長幅度較大的區(qū)域主要分布在黃土高原以及西藏、青海與甘肅等，而減少幅度較大的區(qū)域主要分布在東北地區(qū)，長江流域中上游、長江三角洲以及部分大城市。2000—2017年，“南豬北移”的生豬生產(chǎn)區(qū)域演變格局顯現(xiàn)。

2）中國市域生豬出欄密度時(shí)空格局演變受區(qū)域本底因素與內(nèi)外部驅(qū)動(dòng)因素綜合作用。不同生豬出欄水平的影響因素、各因素的影響強(qiáng)度及作用方向存在差異?？v向來看，隨著全球化發(fā)展、交通等基礎(chǔ)設(shè)施完善、規(guī)?；M(jìn)程加快，飼料資源、水資源、區(qū)位條件等因素的作用減弱，地形條件等的影響加?。粰M向來看，隨著生豬出欄水平的增加，年平均氣溫、飼料資源等對(duì)生豬出欄密度的直接作用減弱，區(qū)位條件的作用增強(qiáng)。內(nèi)外部驅(qū)動(dòng)因素對(duì)生豬出欄密度的影響較為復(fù)雜，其中經(jīng)濟(jì)發(fā)展水平、鄉(xiāng)村非農(nóng)化水平與生豬出欄密度呈倒“U”關(guān)系。

4.2 討論

綜合考慮各地生豬生產(chǎn)水平、發(fā)展趨勢以及資源稟賦、社會(huì)經(jīng)濟(jì)發(fā)展水平、環(huán)境承載力等因素，本研究對(duì)我國重要的生豬生產(chǎn)主產(chǎn)區(qū)和潛力發(fā)展區(qū)提出以下差別化的分區(qū)優(yōu)化策略。

1）東北地區(qū)水土資源富集，飼料資源豐富且成本較低，是我國生豬生產(chǎn)的潛力增長區(qū)。未來應(yīng)加大生豬產(chǎn)業(yè)布局，提升玉米等飼料就地轉(zhuǎn)化率，同時(shí)加大大豆蛋白飼料作物的種植，建設(shè)一批高標(biāo)準(zhǔn)種養(yǎng)結(jié)合生產(chǎn)基地，提升農(nóng)業(yè)綜合競爭力；積極推行經(jīng)濟(jì)高效的糞污收集處理利用的無害化設(shè)施裝備，解決冷涼氣候?qū)S便資源化利用不足等問題；此外，努力補(bǔ)齊畜產(chǎn)品加工流通短板，提升生豬市場營銷能力。

2）隨著“南豬北移”的趨勢不斷強(qiáng)化，黃淮海平原逐漸成為我國生豬供給的核心區(qū)域，但同時(shí)也面臨著水土資源短缺加劇的挑戰(zhàn)，傳統(tǒng)要素與經(jīng)營模式支撐下的生豬生產(chǎn)增長潛力有限。未來應(yīng)充分發(fā)揮傳統(tǒng)糧食種植和牲畜養(yǎng)殖比較優(yōu)勢，鼓勵(lì)種養(yǎng)結(jié)合的生豬養(yǎng)殖模式[32]；在穩(wěn)定現(xiàn)有生產(chǎn)規(guī)模的基礎(chǔ)上，依靠技術(shù)進(jìn)步實(shí)現(xiàn)資源的高效利用，提高集約化水平，促進(jìn)生豬產(chǎn)業(yè)的可持續(xù)發(fā)展。

3）長江流域作為我國第一大流域，水熱條件優(yōu)越，農(nóng)業(yè)資源稟賦高，人口密集，市場廣闊，是我國重要的生豬主產(chǎn)區(qū)，但由于勞動(dòng)力非農(nóng)化進(jìn)程加快和環(huán)境治理壓力加大，生豬生產(chǎn)呈明顯下降趨勢。未來應(yīng)精準(zhǔn)評(píng)估區(qū)內(nèi)生豬生產(chǎn)的環(huán)境承載力差異，在長江中上游地區(qū)應(yīng)遵循“共抓大保護(hù)、不搞大開發(fā)”的策略，調(diào)減生豬養(yǎng)殖規(guī)模；在長江入?？诘拳h(huán)境承載力較高的地區(qū)應(yīng)保持適度規(guī)模的生豬生產(chǎn)，加快推動(dòng)生豬生產(chǎn)向綠色化發(fā)展。

4）云貴地區(qū)生豬出欄密度整體大幅度增加，但受地形條件影響，目前散養(yǎng)戶居多，大規(guī)模養(yǎng)殖戶發(fā)展較慢，同時(shí)由于區(qū)位條件的限制，遠(yuǎn)離消費(fèi)市場。未來應(yīng)適度提高生豬規(guī)模化水平，完善交通等基礎(chǔ)設(shè)施，加快物流產(chǎn)業(yè)發(fā)展，擴(kuò)大生豬市場；同時(shí)大力推廣滇南小耳豬、撒壩豬等地方優(yōu)良豬種。

5）東南沿海地區(qū)在水熱條件、環(huán)境承載力、交通區(qū)位等方面優(yōu)勢明顯，在當(dāng)前生豬供需矛盾加大的背景下，應(yīng)穩(wěn)定生豬生產(chǎn)布局和現(xiàn)有生產(chǎn)規(guī)模，重點(diǎn)通過發(fā)展適度規(guī)模化養(yǎng)殖提高勞動(dòng)生產(chǎn)率，避免生豬養(yǎng)殖業(yè)的過快下滑；同時(shí)要瞄準(zhǔn)高端市場，強(qiáng)化品牌建設(shè)，提升生豬產(chǎn)業(yè)競爭力。

本文僅選取典型截面數(shù)據(jù)開展生豬生產(chǎn)時(shí)空格局演化特征及其影響因素的探究，一方面對(duì)各地區(qū)生豬生產(chǎn)動(dòng)態(tài)規(guī)律的揭示不夠充分，另一方面對(duì)經(jīng)濟(jì)要素對(duì)于生豬生產(chǎn)分布的復(fù)雜影響刻畫不足。未來應(yīng)進(jìn)行面板數(shù)據(jù)的深入分析，更加細(xì)致地揭示各地區(qū)生豬生產(chǎn)發(fā)展演化特征及驅(qū)動(dòng)機(jī)制。此外，本文對(duì)于養(yǎng)殖技術(shù)、規(guī)?；?、全球化等對(duì)生豬生產(chǎn)的影響探討不足，未來應(yīng)充分利用農(nóng)業(yè)大數(shù)據(jù)等多源數(shù)據(jù)進(jìn)一步開展深入研究。