畢雪昊 鄒 偉

1 引 言

鄉(xiāng)村振興戰(zhàn)略是中國現(xiàn)代化進(jìn)程中促進(jìn)農(nóng)業(yè)農(nóng)村發(fā)展的綱領(lǐng)性戰(zhàn)略，產(chǎn)業(yè)繁榮發(fā)展、環(huán)境治理優(yōu)美是鄉(xiāng)村振興實現(xiàn)的重要表現(xiàn)。產(chǎn)業(yè)興旺是振興鄉(xiāng)村的首要要求，也是鄉(xiāng)村發(fā)展的經(jīng)濟(jì)基礎(chǔ)。在廣大農(nóng)村地區(qū)，種植業(yè)與養(yǎng)殖業(yè)是兩大農(nóng)業(yè)支柱產(chǎn)業(yè)，既肩負(fù)著國家糧食安全的重任，也是保障城鄉(xiāng)居民食品消費需求的基礎(chǔ)。實現(xiàn)農(nóng)村產(chǎn)業(yè)興旺需要促進(jìn)和落實種植業(yè)和養(yǎng)殖業(yè)的良好、安全發(fā)展（陳政等，2020［1］）。生態(tài)宜居是鄉(xiāng)村振興的關(guān)鍵也是鄉(xiāng)村環(huán)境需求和綠色生態(tài)產(chǎn)業(yè)發(fā)展導(dǎo)向。近年來養(yǎng)殖業(yè)發(fā)展迅猛，但伴隨著大量養(yǎng)殖廢棄物的產(chǎn)生及不當(dāng)排放，引發(fā)了水體污染、農(nóng)業(yè)面源污染，對鄉(xiāng)村生態(tài)環(huán)境造成嚴(yán)重破壞（譚永風(fēng)等，2021［2］；Jones 等,2019［3］）。養(yǎng)殖廢棄污染問題根源于種養(yǎng)業(yè)發(fā)展不協(xié)調(diào)。不合理的產(chǎn)業(yè)空間布局導(dǎo)致的“種養(yǎng)分離”不僅會增加養(yǎng)殖業(yè)廢棄物還田成本強化排污風(fēng)險，而且不利于種植業(yè)化肥減量和有機(jī)肥資源化利用，影響農(nóng)業(yè)面源污染防治。因此有必要治理“種養(yǎng)分離”問題，促進(jìn)種植業(yè)與養(yǎng)殖業(yè)協(xié)調(diào)發(fā)展，而系統(tǒng)認(rèn)知種養(yǎng)殖業(yè)空間協(xié)調(diào)性，正是促進(jìn)協(xié)調(diào)發(fā)展的基礎(chǔ)。

為了對養(yǎng)殖業(yè)污染進(jìn)行有效管控，促進(jìn)產(chǎn)業(yè)綠色可持續(xù)發(fā)展，中央及地方頒布了環(huán)境規(guī)制。《畜禽規(guī)模養(yǎng)殖污染防治防治條例》（2014）、《畜牧法》（2015）、《水污染防止行動計劃》（2015）、《畜禽養(yǎng)殖禁養(yǎng)區(qū)劃定技術(shù)指南》（2016）、《畜禽規(guī)模養(yǎng)殖場糞污資源化利用設(shè)施建設(shè)規(guī)范》（2018）、《關(guān)于促進(jìn)畜牧業(yè)高質(zhì)量發(fā)展的意見》（2020）等法律法規(guī)和政策文件都對養(yǎng)殖的污染治理做出了安排。通過與種植業(yè)有效結(jié)合，將廢棄物資源化利用就地還田，是降低養(yǎng)殖污染物排放最有效且經(jīng)濟(jì)可行的措施之一（付強等，2020［4］）。環(huán)境規(guī)制會倒逼養(yǎng)殖產(chǎn)業(yè)轉(zhuǎn)移，形成新的產(chǎn)業(yè)布局，而種養(yǎng)業(yè)發(fā)展協(xié)調(diào)程度是對這一影響的直接響應(yīng)，并會對養(yǎng)殖污染治理效能形成面上反饋。探究環(huán)境規(guī)制對種養(yǎng)業(yè)協(xié)調(diào)發(fā)展的影響關(guān)系，既是對中國農(nóng)業(yè)發(fā)展布局研究的進(jìn)一步深入，也是對區(qū)域產(chǎn)業(yè)融合研究的拓展，對廣泛理解環(huán)境規(guī)制效用，深刻認(rèn)識產(chǎn)業(yè)發(fā)展的政策響應(yīng)機(jī)制，完善環(huán)境政策有著重要的現(xiàn)實意義。

已有研究在種養(yǎng)產(chǎn)業(yè)發(fā)展及其影響因素方面已形成了一定的成果，但仍存在可以深入分析的空間。在種植業(yè)與養(yǎng)殖業(yè)產(chǎn)業(yè)發(fā)展方面，現(xiàn)有文獻(xiàn)基于市級、省級層面的數(shù)據(jù)進(jìn)行測算，對產(chǎn)業(yè)布局和變遷規(guī)律進(jìn)行了總結(jié)（黃愛軍，1995［5］；熊學(xué)振等，2022［6］），在種植業(yè)發(fā)展方面，學(xué)者們總結(jié)了不同區(qū)域種植業(yè)的空間演進(jìn)特征。欒軍強等（2021）［7］指出山東省糧食供需重心西向趨勢明顯，劉冬等（2021）［8］針對陜西省各個地級市的不同農(nóng)作物進(jìn)行了空間格局分析，liu 等（2018）［9］在對中國1980-2011 年間的種植結(jié)構(gòu)的時空變化進(jìn)行了刻畫和總結(jié)。有關(guān)養(yǎng)殖業(yè)空間格局的研究同樣較為成熟。針對肉羊（樊宏霞等，2012［10］）、生豬（張園園等，2019［11］）、蛋雞（曹光喬等，2010［12］）等具體的畜禽養(yǎng)殖產(chǎn)業(yè)，學(xué)者們對其空間集聚規(guī)律都進(jìn)行了細(xì)致的梳理。此外，李洋等（2021）［13］著眼于種養(yǎng)業(yè)結(jié)構(gòu)優(yōu)化，基于比較優(yōu)勢理論，從運用334 個地級市數(shù)據(jù)，對中國種植業(yè)和養(yǎng)殖業(yè)生產(chǎn)格局進(jìn)行分析，指出當(dāng)前種養(yǎng)業(yè)結(jié)構(gòu)的合理化有待加強。有關(guān)種養(yǎng)產(chǎn)業(yè)結(jié)合的影響因素，學(xué)者們也進(jìn)行了深入探索。黃煒虹等（2017）［14］基于農(nóng)戶調(diào)查數(shù)據(jù)指出經(jīng)濟(jì)收益會影響種植和養(yǎng)殖聯(lián)合生產(chǎn)的意愿與行為。從范圍經(jīng)濟(jì)的角度考慮，與單一從事種植業(yè)或養(yǎng)殖業(yè)相比，種養(yǎng)結(jié)合式的生產(chǎn)有著更高的效率，農(nóng)戶對生產(chǎn)效率的追求也是促進(jìn)種養(yǎng)業(yè)協(xié)調(diào)發(fā)展的主要原因之一（錢忠好和李友藝，2020［15］；侯國慶等，2022［16］）。進(jìn)一步地，因為效率由生產(chǎn)要素的投入和組合決定，生產(chǎn)要素獲取的成本也成為了種養(yǎng)業(yè)能否結(jié)合的重要依據(jù)（彭艷玲等，2019［17］；姜天龍等，2022［18］）。此外，農(nóng)戶的個體特征也會導(dǎo)致其對是否從事種養(yǎng)結(jié)合生產(chǎn)有著不同的偏好（Huang 等, 2022［19］）。綜上可以發(fā)現(xiàn)，在產(chǎn)業(yè)發(fā)展空間格局的研究尺度方面，當(dāng)前已有研究偏向于從省域、市域等相對宏觀的研究尺度進(jìn)行了空間分析（李洋等，2021［13］）。側(cè)重單一產(chǎn)業(yè)發(fā)展空間形態(tài)的文獻(xiàn)較為充分，將多產(chǎn)業(yè)發(fā)展進(jìn)行聯(lián)合討論的研究還需要深入。在影響因素探析方面，從政策層面探究種養(yǎng)業(yè)協(xié)調(diào)發(fā)展響應(yīng)這一議題有著更廣泛的討論空間（周建軍等，2018［20］）。已有部分文獻(xiàn)（冀名峰等，2019［21］；潘丹，2015［22］；孟祥海等，2022［23］；孟祥海等，2019［24］）從農(nóng)業(yè)補貼等經(jīng)濟(jì)政策進(jìn)行分析，這部分文獻(xiàn)主要側(cè)重于在經(jīng)濟(jì)收益的影響框架下進(jìn)行探討，從環(huán)境政策角度出發(fā)的影響分析能夠為該方面的研究提供新的視角。對于種植業(yè)和養(yǎng)殖業(yè)而言，其有效連接的關(guān)鍵在于養(yǎng)殖業(yè)產(chǎn)生的廢棄物能夠及時資源化利用，為種植業(yè)提供養(yǎng)分。不同于其他產(chǎn)業(yè)，農(nóng)業(yè)利潤率極低（盛來運，2007［25］），出于成本控制的原則，在小范圍區(qū)域內(nèi)的廢棄物就地還田更能保證種養(yǎng)結(jié)合的經(jīng)濟(jì)可持續(xù)性，因而較小的研究尺度更能精確反應(yīng)產(chǎn)業(yè)協(xié)調(diào)發(fā)展程度。因此本文將基于以往研究提供的經(jīng)驗，進(jìn)一步縮小研究尺度，以縣域單位作為研究對象，對我國種養(yǎng)業(yè)發(fā)展情況進(jìn)行空間分析。在方法上，本文將采用一致性系數(shù)衡量各縣域單位的種養(yǎng)業(yè)協(xié)調(diào)發(fā)展程度，以期為種養(yǎng)結(jié)合相關(guān)研究提供更多的文獻(xiàn)作用與研究參考。同時，環(huán)境規(guī)制政策是重要的農(nóng)業(yè)政策之一，污染防治、資源利用、綠色發(fā)展是種養(yǎng)業(yè)協(xié)調(diào)發(fā)展的主要目的和重要意義，也是環(huán)境政策實施的出發(fā)點。分析環(huán)境規(guī)制對種養(yǎng)業(yè)協(xié)調(diào)發(fā)展的影響能夠拓寬該領(lǐng)域的研究內(nèi)容，提供新的研究視角，因此本文還將重點討論環(huán)境規(guī)制對種養(yǎng)協(xié)調(diào)發(fā)展的影響，并通過環(huán)境成本理論對影響效果進(jìn)行解釋。基于以上，本文以縣域空間為研究尺度，在分析種植業(yè)與養(yǎng)殖業(yè)產(chǎn)業(yè)發(fā)展空間格局的基礎(chǔ)上，總結(jié)二者協(xié)調(diào)發(fā)展規(guī)律，并探究環(huán)境規(guī)制下養(yǎng)殖業(yè)的產(chǎn)業(yè)轉(zhuǎn)移效應(yīng)及環(huán)境規(guī)制對種養(yǎng)業(yè)協(xié)調(diào)發(fā)展的影響，可以為進(jìn)一步促進(jìn)種養(yǎng)業(yè)協(xié)調(diào)發(fā)展，改善農(nóng)業(yè)生產(chǎn)生活環(huán)境提供實證依據(jù)。

2 數(shù)據(jù)與研究方法

2.1 數(shù)據(jù)來源

本文以糧食產(chǎn)量表征區(qū)域種植業(yè)發(fā)展水平，以肉類總產(chǎn)量表征區(qū)域養(yǎng)殖業(yè)發(fā)展水平，相關(guān)數(shù)據(jù)來源于2002-2020 年的《中國縣域統(tǒng)計年鑒》、各省（市、自治區(qū)）、市（自治州）、縣統(tǒng)計年鑒以及各地的經(jīng)濟(jì)社會發(fā)展統(tǒng)計公報。參考以往研究（虞祎等，2011［26］）以環(huán)境規(guī)制政策個數(shù)來表征環(huán)境規(guī)制變量，相關(guān)數(shù)據(jù)來源于中國法律數(shù)字圖書館①https://law.cnki.net/的法律法規(guī)庫。分別通過設(shè)置檢索關(guān)鍵詞“畜禽 污染”、“畜禽 糞污”、“畜禽 廢棄物”，檢索后獲取各地的養(yǎng)殖業(yè)環(huán)境規(guī)制個數(shù)，并根據(jù)行政管轄關(guān)系，將國家、各省、各地級市頒布的環(huán)境規(guī)制政策數(shù)量加至各縣級單位單獨頒布的政策數(shù)量上，形成對該變量的賦值。本文研究單元為縣級行政區(qū)，考慮到行政區(qū)劃的變化，以2019 年的縣級行政區(qū)劃為基準(zhǔn)，并對之前的縣級行政單位進(jìn)行修正處理。

2.2 研究方法

2.2.1 區(qū)域重心分析法

重心是源于物理力學(xué)中的概念，指空間存在一點使其各個方向上的力量對比保持相對均衡，區(qū)域重心被廣泛運用到各種空間研究領(lǐng)域。本文通過區(qū)域重心模型分析縣域種植業(yè)與養(yǎng)殖業(yè)的產(chǎn)業(yè)空間分布格局和轉(zhuǎn)移過程，模型基本設(shè)定如下：

其中，Xj、Yj分別表示第j 年產(chǎn)業(yè)重心所在地位置的經(jīng)度和緯度；xi、yi分別表示第i 個縣級行政單位空間質(zhì)心的經(jīng)度和緯度；Tij表示第i 個縣級單位第j 年的產(chǎn)業(yè)產(chǎn)量。

基于不同時期的產(chǎn)業(yè)重心坐標(biāo)可以計算出重心移動距離，公式如下：

其中，D（t+k）-t表示t 至t+k 時期的產(chǎn)業(yè)重心移動距離；l 為地球表面坐標(biāo)單位（度）轉(zhuǎn)化為平面距離單位（千米）的系數(shù)；Xt、Yt分別表示t 時期重心所在地的經(jīng)緯坐標(biāo)。

2.2.2 空間重心重疊性與變動一致性

采用空間重心重疊性系數(shù)，從靜態(tài)視角衡量種植業(yè)和養(yǎng)殖業(yè)產(chǎn)業(yè)重心的重合程度，具體計算公式如下（謝坤等，2021［27］）：

其中，（Xp，Yp）和（Xb，Yb）分別表示種植業(yè)與養(yǎng)殖業(yè)的空間重心坐標(biāo)；S 表示種植業(yè)產(chǎn)業(yè)重心與養(yǎng)殖業(yè)產(chǎn)業(yè)中心的空間距離，單位為千米。S 值越小，表明二者空間重心重疊性越高，則種植業(yè)與養(yǎng)殖業(yè)空間重合度越高，反之，空間重合度越低，S 為0 意味著種養(yǎng)業(yè)重心完全重疊，空間分布一致。

采用變動一致性系數(shù)，從動態(tài)視角衡量種植業(yè)與養(yǎng)殖業(yè)的重心轉(zhuǎn)移規(guī)律的一致性程度。一致性程度可以通過種養(yǎng)業(yè)二者重心相較前一個時間點產(chǎn)生位移的矢量角度θ 表現(xiàn)，角度越小說明二者重心變動越一致。由于θ ∈［0,180］，因此將其余弦值設(shè)置為變動一致性系數(shù)（C），取值范圍為［-1,1］，余弦值越大說明種養(yǎng)業(yè)重心變動越一致，反之說明變動越不一致，當(dāng)C 為-1 時，說明二者重心變化方向完全相反，當(dāng)C 為1 是說明遷移方向相同。根據(jù)余弦定理，變動一致性系數(shù)如下：

其中，Δ 表示與前一期相比的坐標(biāo)變化量；其余變量含義同前式。

2.2.3 協(xié)調(diào)度指數(shù)

為分析各地區(qū)養(yǎng)殖業(yè)與種植業(yè)發(fā)展在空間上的協(xié)調(diào)性，借鑒人口與經(jīng)濟(jì)空間一致性系數(shù)（封志明和劉曉娜，2013［28］）來表征養(yǎng)殖業(yè)與種植業(yè)的產(chǎn)業(yè)集中度的關(guān)系，即種養(yǎng)協(xié)調(diào)度指數(shù)，具體計算公式如下：

其中，Ki表示某年i 縣的種養(yǎng)業(yè)協(xié)調(diào)度，K 小于1 表明養(yǎng)殖業(yè)產(chǎn)業(yè)集聚低于種植業(yè)，K 大于1 表明養(yǎng)殖業(yè)產(chǎn)業(yè)集聚高于種植業(yè)，K 等于1 表明養(yǎng)殖業(yè)與種植業(yè)一致，發(fā)展較為協(xié)調(diào)；Bi、Pi分別表示該年i 縣養(yǎng)殖業(yè)產(chǎn)量與種植業(yè)產(chǎn)量。

參考相關(guān)研究（謝坤等，2021［27］），將K 取值劃分為5 段，每段代表不同的種養(yǎng)業(yè)發(fā)展類別。當(dāng)協(xié)調(diào)度指數(shù)小于0.5，意味著該區(qū)域種植業(yè)發(fā)展程度遠(yuǎn)高于養(yǎng)殖業(yè)，二者協(xié)調(diào)度很差，種植業(yè)有機(jī)肥施用情況不樂觀；當(dāng)K ∈（0.5,0.8］時，說明種植業(yè)發(fā)展程度略高于養(yǎng)殖業(yè)；協(xié)調(diào)度指數(shù)位于0.8 至1.2 之間表明種養(yǎng)業(yè)協(xié)調(diào)度較好，兩產(chǎn)業(yè)發(fā)展較均衡；K ∈（1.2,2］意味著養(yǎng)殖業(yè)發(fā)展程度略高于種植業(yè)，養(yǎng)殖產(chǎn)生的廢棄物數(shù)量略多，有一定的環(huán)境污染風(fēng)險；當(dāng)協(xié)調(diào)指數(shù)超過2 時，說明養(yǎng)殖業(yè)產(chǎn)業(yè)集聚遠(yuǎn)高于種植業(yè)，二者發(fā)展不協(xié)調(diào)。

表1 種養(yǎng)業(yè)協(xié)調(diào)度分級

2.2.4 面板回歸模型

為考察環(huán)境規(guī)制與種養(yǎng)協(xié)調(diào)度的關(guān)系，設(shè)定如下的面板回歸模型進(jìn)行分析：

其中，i 和t 分別表示第i 個縣級行政單位和第t 年；E 為環(huán)境規(guī)制變量，用環(huán)境規(guī)制政策個數(shù)表征；X為5 個控制變量，根據(jù)農(nóng)業(yè)生產(chǎn)特征，并參考相關(guān)文獻(xiàn)，從經(jīng)濟(jì)、社會和技術(shù)等維度選取控制變量：經(jīng)濟(jì)因素包括地區(qū)生產(chǎn)總值（單位：億元）、第一產(chǎn)業(yè)增加值（單位：億元），社會因素包括年末總?cè)丝冢▎挝唬喝f人）、人口密度（單位：人/平方公里），技術(shù)因素包括農(nóng)業(yè)機(jī)械總動力（單位：萬千瓦特）；α、λ、β均為變量的系數(shù)；L 為常數(shù)項；μi、νt、εit分別為地區(qū)固定效應(yīng)、時間固定效應(yīng)和殘差項。

3 種養(yǎng)業(yè)發(fā)展協(xié)調(diào)度的空間演變

3.1 種植業(yè)與養(yǎng)殖業(yè)重心變化特征

從圖1 和表2 可以發(fā)現(xiàn)，養(yǎng)殖業(yè)生產(chǎn)重心表現(xiàn)出分布范圍小、遷移速度慢、遷移方向多樣的特點。整體上看，2001-2019 年期間，養(yǎng)殖業(yè)生產(chǎn)重心相對集中，主要分布在河南省的西南部，集中在方城縣、社旗縣、唐河縣、驛城區(qū)和泌陽縣境內(nèi)。在研究期內(nèi)，養(yǎng)殖業(yè)重心坐標(biāo)由東經(jīng)113°7′22″，北緯32°41′3″（河南省南陽市唐河縣）變動至東經(jīng)113°48′2″，北緯 32°58′53″（河南省駐馬店市驛城區(qū)），整體移動了70.86km，向東移動59.36km，向北移動38.70km。在重心遷移的每個階段內(nèi)，遷移距離基本在40 公里內(nèi)，僅2016-2019年間移動距離較長。每階段遷移方向雖略有差別，但是整體上我國養(yǎng)殖業(yè)重心是在向東向北轉(zhuǎn)移，與具體養(yǎng)殖產(chǎn)業(yè)的轉(zhuǎn)移規(guī)律，如“南豬北養(yǎng)”、“西豬東進(jìn)”、“南牛北移”較為相符（趙俊偉等，2019［29］；孫爽軒等，2020［30］），說明在這一時期內(nèi)，東部地區(qū)和北部地區(qū)的養(yǎng)殖業(yè)生產(chǎn)總量在全國范圍內(nèi)的比重有一定的提升。

圖1 2001-2019 年中國種植業(yè)生產(chǎn)重心與養(yǎng)殖業(yè)生產(chǎn)重心動態(tài)變化軌跡

圖2 2001-2019 年中國縣域種植業(yè)與養(yǎng)殖業(yè)重心重疊性與一致性

表2 2001-2019 年種植業(yè)與養(yǎng)殖業(yè)重心遷移結(jié)果

種植業(yè)生產(chǎn)重心呈現(xiàn)出分布范圍廣、遷移速度快，遷移方向在研究期內(nèi)較為一致的特點。從整體上來看，種植業(yè)重心在研究期內(nèi)相對分散，分布在河南的通許縣、滑縣、鄢陵縣、舞陽縣以及山東省莘縣。從整體上看，在研究期內(nèi)，種植業(yè)重心坐標(biāo)由東經(jīng)113°50′48″，北緯33°32′33″（河南省漯河市舞陽縣）變動至東經(jīng)115°34′17″，北緯35°56′54″（山東省聊城市莘縣），向東移動了127.88km，向北移動了285.49km，整體移動了312.82km。在重心遷移的各階段中除了2013-2016 年間，其余階段種植業(yè)生產(chǎn)重心遷移距離均超過40 公里，遷移速度較養(yǎng)殖業(yè)更快。然而除2010 至2016 年階段，重心均沿著東北方向遷移，說明近20 年我國的種植業(yè)在快速地向東向北轉(zhuǎn)移，相比之下西部和南部的種植業(yè)生產(chǎn)占全國的比重有所下降。

綜合二者的生產(chǎn)重心遷移路徑來看，2001 年至2019 年間兩產(chǎn)業(yè)重心基本集中在河南省境內(nèi)，但是存在一定的偏離，二者產(chǎn)業(yè)轉(zhuǎn)移方向在整體上較為同向，即向東向北。造成這一結(jié)果可能的原因是，東部和北部地區(qū)的土地資源更為豐富，為種植業(yè)發(fā)展提供了較好的要素條件（徐春春等，2013［31］），同時，種植業(yè)的產(chǎn)品即為養(yǎng)殖業(yè)的生產(chǎn)資料，飼料供應(yīng)是影響?zhàn)B殖業(yè)的重要原因之一，因而養(yǎng)殖業(yè)的變遷方向會呈現(xiàn)出與種植業(yè)變遷的相對一致性（衛(wèi)龍寶和張菲，2012［32］；王歡和喬娟，2017［33］）。

3.2 種養(yǎng)業(yè)重心重疊性與一致性

2001-2019 年間，中國種植業(yè)與養(yǎng)殖業(yè)空間分布存在一定的空間錯位，但是二者在空間遷移的方向上較為相似。從種植業(yè)與養(yǎng)殖業(yè)重心之間的距離來看，二者重心錯位距離在117.41km 至371.12 公里之間波動，表現(xiàn)出“增-減-增”的變化態(tài)勢，2010 前，種植業(yè)與養(yǎng)殖業(yè)的重心距離表現(xiàn)出逐漸遞增的趨勢，全國層面上二者的空間一致性在逐步降低。2010-2013 年，二者重心重疊性增大，表現(xiàn)為空間上的距離縮短。而2013年后，種植業(yè)與養(yǎng)殖業(yè)生產(chǎn)重心距離進(jìn)一步拉大，研究期內(nèi)二者的生產(chǎn)重心重疊性小，整體上呈現(xiàn)出一定的空間錯位。從種植業(yè)與養(yǎng)殖業(yè)重心遷移的變動方向來看，二者變動方向較為相似，整體上都是向東向北進(jìn)行遷移，但遷移速度差距大導(dǎo)致了重心重疊性小。在2001-2004 年間以及2007 至2010 年階段，中國種植業(yè)與養(yǎng)殖業(yè)生產(chǎn)重心基本保持了同向遷移，一致性程度高。在2010-2013 年間，二者生產(chǎn)重心遷移方向相反，但是相向而行，因而這一階段的重心重疊性也較高。其余階段的重心一致性變化量較小，均在0 值附近徘徊?？偟膩碚f，種植業(yè)與養(yǎng)殖業(yè)重心變動一致性指數(shù)較高，但由于重心遷移速度不同，仍存在一定的空間錯位。

3.3 種養(yǎng)業(yè)協(xié)調(diào)度的時空特征

圖3 為不同等級協(xié)調(diào)度的縣級行政單位數(shù)量分別在全國、東部地區(qū)、中部地區(qū)和西部地區(qū)的占比①根據(jù)地理位置及經(jīng)濟(jì)發(fā)展發(fā)展情況，將我國大陸地區(qū)劃分為東、中、西部分，其中東部地區(qū)包括北京、天津、河北、遼寧、上海、江蘇、浙江、福建、山東、廣東、海南等11 個?。ㄊ校?中部地區(qū)包括山西、吉林、黑龍江、安徽、江西、河南、湖北、湖南，共8 個省級行政區(qū)；西部地區(qū)包括四川、重慶、貴州、云南、西藏、陜西、甘肅、青海、寧夏、新疆、廣西、內(nèi)蒙古，共12 個省級行政區(qū)。?？偟膩砜?，在2016 年及之前中國種養(yǎng)協(xié)調(diào)度指數(shù)空間分布較為均衡且養(yǎng)殖業(yè)與種植業(yè)聚集較一致的縣級行政單位占比較多。2019 年的數(shù)據(jù)顯示養(yǎng)殖業(yè)聚集程度超過種植業(yè)的區(qū)縣數(shù)量遠(yuǎn)超種植業(yè)相對聚集的區(qū)縣，這可能是由于2017-2018 年間非洲豬瘟、禽流感等畜禽疫病頻發(fā)，對國內(nèi)養(yǎng)殖產(chǎn)業(yè)形成了較大的沖擊，導(dǎo)致養(yǎng)殖場大規(guī)模搬遷，養(yǎng)殖業(yè)更加集中。分地區(qū)來看，東部和西部地區(qū)的協(xié)調(diào)度分布變化趨勢和全國趨勢類似，同時長期以來，東西部地區(qū)的養(yǎng)殖業(yè)產(chǎn)量集中度基本高于種植業(yè)產(chǎn)量集中度，說明在這兩個地區(qū)養(yǎng)殖業(yè)發(fā)展更優(yōu)。中部地區(qū)則呈現(xiàn)出相反的分布。在研究期內(nèi)，中部地區(qū)種植業(yè)產(chǎn)量聚集程度更高，這可能是由于中部省份，如黑龍江、吉林、河南、安徽等地，地勢平坦，土壤肥沃、自然資源豐富，更適宜種植類產(chǎn)業(yè)發(fā)展，因而在中部省份的種植業(yè)集中度相對較高。在東部地區(qū)，發(fā)展種植業(yè)所需的土地、耕地資源在數(shù)量上更為緊張，相比之下種植業(yè)集中度略遜于養(yǎng)殖業(yè)。由于氣候、地形、水肥等自然因素，部分西部省份的耕地資源相對匱乏，但是牧草地等養(yǎng)殖用地資源豐富，導(dǎo)致此區(qū)域的養(yǎng)殖業(yè)發(fā)展較種植業(yè)更為發(fā)達(dá)。

圖3 全國及各地區(qū)不同等級種養(yǎng)業(yè)協(xié)調(diào)度占比

4 環(huán)境規(guī)制對種養(yǎng)業(yè)協(xié)調(diào)度的影響

環(huán)境規(guī)制對種養(yǎng)業(yè)協(xié)調(diào)度的影響主要是通過其對養(yǎng)殖業(yè)的作用實現(xiàn)的，因此采用面板回歸模型，運用STATA15.0 軟件估計分別估計環(huán)境規(guī)制對養(yǎng)殖業(yè)及種養(yǎng)業(yè)協(xié)調(diào)度的影響效應(yīng)，回歸結(jié)果見表3。

4.1 環(huán)境規(guī)制對養(yǎng)殖產(chǎn)業(yè)的影響

模型1 為環(huán)境規(guī)制影響?zhàn)B殖業(yè)產(chǎn)量的回歸結(jié)果，由于種植業(yè)的產(chǎn)品及副產(chǎn)物為養(yǎng)殖提供了生產(chǎn)的投入要素，因此控制變量的選擇在式（7）的基礎(chǔ)上增加種植業(yè)發(fā)展水平這一變量，與環(huán)境規(guī)制變量一起共同組成該模型的自變量。由結(jié)果可知環(huán)境規(guī)制變量的回歸系數(shù)在1%的顯著性水平下為負(fù)，說明地區(qū)環(huán)境政策實施個數(shù)越多、環(huán)境管制力度越強對養(yǎng)殖業(yè)發(fā)展越不利。事實上，長期以來養(yǎng)殖業(yè)生產(chǎn)帶來的固液氣廢棄物排放一直對周圍的生態(tài)環(huán)境和人居環(huán)境產(chǎn)生著惡劣的影響。在實施環(huán)境規(guī)制政策之前，養(yǎng)殖主體以幾乎零環(huán)境成本的運作模式進(jìn)行生產(chǎn)。伴隨著規(guī)制的執(zhí)行，環(huán)境成本開始成為養(yǎng)殖成本中的主要部分，其比重也隨著規(guī)制的嚴(yán)格程度的加深而不斷加重，給生產(chǎn)增加了負(fù)擔(dān)。由于地區(qū)間環(huán)境規(guī)制實施的差異，養(yǎng)殖產(chǎn)品的社會生產(chǎn)價格因此發(fā)生大幅增長，規(guī)制嚴(yán)格地區(qū)養(yǎng)殖廠商的生產(chǎn)利潤因此受到壓縮，這會迫使其從當(dāng)?shù)厥袌鲛D(zhuǎn)出，造成地區(qū)產(chǎn)量的下降。“赫-俄理論”認(rèn)為地區(qū)會根據(jù)自身生產(chǎn)要素稟賦的豐裕程度來決定生產(chǎn)方式。可以將環(huán)境要素納入養(yǎng)殖生產(chǎn)過程，那么地區(qū)間環(huán)境要素稟賦的差距就會引起養(yǎng)殖業(yè)的遷入和遷出。發(fā)展養(yǎng)殖業(yè)需要的豐裕的環(huán)境要素稟賦有兩層內(nèi)涵（虞祎等，2011［26］），一是地區(qū)具有豐富的自然資源和社會條件以消納養(yǎng)殖產(chǎn)生的廢棄物，也即“真實的環(huán)境豐?！保欢堑貐^(qū)放松環(huán)境管制，默許廢棄物和污染排放，此為“虛假的環(huán)境豐裕”。對于環(huán)境政策緊收的地區(qū)來說，適宜養(yǎng)殖業(yè)的環(huán)境稟賦是相對不足的，不僅不具有發(fā)展養(yǎng)殖業(yè)的環(huán)境比較優(yōu)勢，且嚴(yán)格實施的環(huán)境規(guī)制也增加了在當(dāng)?shù)剡M(jìn)行養(yǎng)殖生產(chǎn)的比較成本。此外，“污染天堂理論”也同樣指出，污染型產(chǎn)業(yè)傾向于從環(huán)境規(guī)制嚴(yán)格的地區(qū)向環(huán)境規(guī)制寬松的地區(qū)轉(zhuǎn)移（Walter 和Ugelow，1979［34］），養(yǎng)殖業(yè)同樣服從這一規(guī)律?？偟膩碚f，環(huán)境規(guī)制是一類管制措施，環(huán)境政策的收緊會使得養(yǎng)殖戶或養(yǎng)殖企業(yè)不可避免地面對合規(guī)成本的增加?；诂F(xiàn)實和理論的分析，執(zhí)行環(huán)境規(guī)制會對養(yǎng)殖產(chǎn)業(yè)形成擠出效應(yīng)，然而這種擠出效應(yīng)是否會進(jìn)一步影響種養(yǎng)業(yè)協(xié)調(diào)則需要進(jìn)行更深入地探究。

從控制變量的回歸結(jié)果來看，各變量均對因變量有著顯著的影響。（1）種植業(yè)發(fā)展的影響在1%水平下顯著為正，說明種植產(chǎn)業(yè)發(fā)達(dá)能夠帶動當(dāng)?shù)仞B(yǎng)殖業(yè)發(fā)展?？赡艿脑蚴?，種植業(yè)為養(yǎng)殖業(yè)發(fā)展提供了飼料這一重要的生產(chǎn)投入元素，因此在種植業(yè)產(chǎn)量較高能夠促進(jìn)本地區(qū)養(yǎng)殖業(yè)發(fā)展。（2）經(jīng)濟(jì)因素。經(jīng)濟(jì)因素包括地區(qū)生產(chǎn)總值與第一產(chǎn)業(yè)地區(qū)生產(chǎn)總值，二者均對養(yǎng)殖業(yè)發(fā)展有著顯著地正向促進(jìn)。這一結(jié)果反應(yīng)了養(yǎng)殖業(yè)會隨著經(jīng)濟(jì)水平的提高而發(fā)展，與其他學(xué)者的研究結(jié)果相似（Han 等，2020［35］）?？赡艿脑蚴牵B(yǎng)殖場或者養(yǎng)殖企業(yè)的建設(shè)需要一定的投資，相比于其他農(nóng)業(yè)產(chǎn)業(yè)，養(yǎng)殖業(yè)更需要資本投入，因此經(jīng)濟(jì)發(fā)達(dá)地區(qū)更具備養(yǎng)殖業(yè)發(fā)展所需的資本條件（楊少華和王凱，2022［36］）。（3）社會因素。地區(qū)人口變量在0.01的顯著性水平上表現(xiàn)出正向促進(jìn)作用，而人口密度變量呈現(xiàn)出顯著的負(fù)向影響?？赡艿脑蚴?，地區(qū)人口數(shù)量表明地區(qū)對肉類產(chǎn)品的需求量，人口數(shù)量越多說明當(dāng)?shù)仞B(yǎng)殖產(chǎn)品的市場需求量越大，有利于區(qū)域養(yǎng)殖產(chǎn)業(yè)的發(fā)展（Zhang 等，2022［37］）。相關(guān)政策的出臺可以對人口密度負(fù)向影響?zhàn)B殖產(chǎn)業(yè)發(fā)展做出了解釋。為了避免養(yǎng)殖污染對周邊居民的影響，《畜禽養(yǎng)殖污染防治管理辦法》、《畜禽養(yǎng)殖禁養(yǎng)區(qū)劃定技術(shù)指南》等政策文件均指出要在人口密度大的地區(qū)實施禁養(yǎng)、限養(yǎng)，鼓勵養(yǎng)殖小區(qū)、養(yǎng)殖場搬遷至低人口密度地區(qū)。因此在人口密度較高的地區(qū)養(yǎng)殖業(yè)發(fā)展往往受到限制。（4）技術(shù)因素。技術(shù)因素以農(nóng)業(yè)機(jī)械總動力進(jìn)行表征，結(jié)果表明農(nóng)業(yè)機(jī)械變量在10%的顯著性水平下負(fù)向影響?zhàn)B殖業(yè)發(fā)展。該變量的影響系數(shù)較小且顯著性水平比其他變量低，說明其對因變量的解釋十分有限。這可能是由于變量選取的代表性不足，導(dǎo)致其解釋力較差。

4.2 環(huán)境規(guī)制對種養(yǎng)業(yè)協(xié)調(diào)度的影響

根據(jù)前文對種養(yǎng)協(xié)調(diào)度的設(shè)定，協(xié)調(diào)度為1 時表明種植業(yè)與養(yǎng)殖業(yè)協(xié)調(diào)發(fā)展程度達(dá)到頂峰，因此有關(guān)環(huán)境規(guī)制對協(xié)調(diào)度影響的回歸估計，應(yīng)當(dāng)通過非線性的形式進(jìn)行表達(dá)。將因變量種養(yǎng)業(yè)協(xié)調(diào)度分為兩個部分，一部分為K≤1，另一部分則是K>1，分別進(jìn)行估計，模型2 和模型3 分別展示了兩階段的結(jié)果。當(dāng)協(xié)調(diào)度指數(shù)小于等于1 時候，環(huán)境規(guī)制變量對其發(fā)揮顯著的正向影響，而協(xié)調(diào)度指數(shù)大于1 時，環(huán)境規(guī)制的回歸系數(shù)顯著為負(fù)。這一結(jié)果說明當(dāng)種植業(yè)產(chǎn)量集聚程度高于養(yǎng)殖業(yè)時，執(zhí)行環(huán)境政策會促進(jìn)地區(qū)養(yǎng)殖業(yè)比重的增加，而當(dāng)養(yǎng)殖業(yè)產(chǎn)量集聚程度高于種植業(yè)時，環(huán)境規(guī)制會促進(jìn)種植業(yè)比重增加，即環(huán)境規(guī)制的實施會修正過高的種植業(yè)或養(yǎng)殖業(yè)產(chǎn)業(yè)集聚，促進(jìn)種養(yǎng)業(yè)協(xié)調(diào)發(fā)展。環(huán)境成本理論可以這一結(jié)果進(jìn)行解釋。前文已述，環(huán)境規(guī)制實施會產(chǎn)生成本效應(yīng)，提高養(yǎng)殖業(yè)生產(chǎn)的成本，主要包括合規(guī)成本和治污成本兩方面（沈坤榮等，2017［38］）。對于養(yǎng)殖主體而言，實施相對寬松的環(huán)境政策的地區(qū)是規(guī)避合規(guī)成本的“污染天堂”，在此地從事養(yǎng)殖業(yè)的生產(chǎn)主體會面臨相對較低的環(huán)境規(guī)制合規(guī)成本（Walter 和 Ugelow，1979［34］）。隨著環(huán)境規(guī)制政策的實施，治污成本成為養(yǎng)殖生產(chǎn)主體的重點關(guān)注對象。養(yǎng)殖業(yè)的污染來源于畜禽糞污和廢棄物，根據(jù)《畜禽養(yǎng)殖業(yè)污染防治技術(shù)政策》，養(yǎng)殖業(yè)治污主要有能源化和肥料化兩種方式。其中糞污肥料化利用更為主流（彭艷玲等，2019［17］）?！秶鴦?wù)院辦公廳關(guān)于加快推進(jìn)畜禽養(yǎng)殖廢棄物資源化利用的意見》中鼓勵養(yǎng)殖糞污近地、就地還田消納，以實現(xiàn)肥料化利用。種植業(yè)發(fā)達(dá)的地區(qū)具備更強的糞污消納條件，有著更強的環(huán)境稟賦優(yōu)勢。因此在環(huán)境規(guī)制的影響下，養(yǎng)殖主體會傾向于將產(chǎn)業(yè)遷移至種植業(yè)發(fā)達(dá)地區(qū)，使得種養(yǎng)同步集聚，實現(xiàn)種養(yǎng)業(yè)協(xié)調(diào)發(fā)展。

從控制變量的回歸結(jié)果來看，僅地區(qū)生產(chǎn)總值與地區(qū)人口對二者協(xié)調(diào)發(fā)展有著顯著的影響。當(dāng)協(xié)調(diào)度指數(shù)小于等于1 時，地區(qū)生產(chǎn)總值的回歸系數(shù)顯著為負(fù)，當(dāng)協(xié)調(diào)度指數(shù)大于1 時，其回歸系數(shù)顯著為正，說明經(jīng)濟(jì)發(fā)達(dá)程度對種養(yǎng)業(yè)協(xié)調(diào)發(fā)展有負(fù)向影響?？赡艿脑蚴?，在經(jīng)濟(jì)發(fā)達(dá)地區(qū)產(chǎn)業(yè)專業(yè)化程度相對較高，因此不利于二者共同發(fā)展。對于地區(qū)人口這一變量而言，協(xié)調(diào)度指數(shù)小于等于1 時其表現(xiàn)出顯著的正向影響，而指數(shù)超過1 時其表現(xiàn)出顯著的負(fù)向影響，說明人口越多的地區(qū)種養(yǎng)業(yè)協(xié)調(diào)發(fā)展程度越高?？赡艿脑蚴?，地區(qū)人口越多，對種植業(yè)與養(yǎng)殖業(yè)產(chǎn)品的需求量越大，因而會促進(jìn)當(dāng)?shù)貎蓚€產(chǎn)業(yè)的共同發(fā)展。

5 結(jié)論與啟示

本文以2001-2019 年縣域?qū)用娴臄?shù)據(jù)為基礎(chǔ)，對中國種養(yǎng)業(yè)協(xié)調(diào)發(fā)展的時空特征進(jìn)行分析，探討了環(huán)境規(guī)制對其的影響。具體結(jié)論如下：

第一，2001-2019 年，種植業(yè)和養(yǎng)殖業(yè)生產(chǎn)重心主要集中在河南省境內(nèi)。養(yǎng)殖業(yè)生產(chǎn)重心遷移范圍相對較小，遷移速度較慢，研究期內(nèi)整體的重心遷移方向為東北方向，而種植業(yè)生產(chǎn)重心分布較為分散，重心遷移速度較快，整體上同樣朝東北方向遷移。

第二，由于種植業(yè)和養(yǎng)殖業(yè)在生產(chǎn)重心遷移速度上具有一定的差距，盡管遷移方向大致相同，但是整體上在研究期內(nèi)兩個產(chǎn)業(yè)的空間重疊性較差，存在一定的錯位分布情況。

第三，在全國層面上，2016 年前種植業(yè)與養(yǎng)殖業(yè)發(fā)展在空間上整體較為平衡，發(fā)展協(xié)調(diào)性有所改善。在地區(qū)層面上，中部地區(qū)表現(xiàn)出更高的種植業(yè)產(chǎn)業(yè)聚集程度，東西部地區(qū)的養(yǎng)殖業(yè)聚集程度普遍高于種植業(yè)。

第四，環(huán)境規(guī)制會抑制地區(qū)養(yǎng)殖業(yè)的發(fā)展，同時會促進(jìn)地區(qū)種養(yǎng)業(yè)協(xié)調(diào)。此外，經(jīng)濟(jì)發(fā)達(dá)程度、地區(qū)人口數(shù)量會對當(dāng)?shù)仞B(yǎng)殖業(yè)發(fā)展產(chǎn)生促進(jìn)作用，而人口密度則會抑制養(yǎng)殖業(yè)發(fā)展。對于種養(yǎng)協(xié)調(diào)發(fā)展而言，地區(qū)生產(chǎn)總值和人口數(shù)量會分別對其產(chǎn)生負(fù)向和正向影響。

上述研究結(jié)論可以帶來如下啟示：

首先，盡管環(huán)境規(guī)制的實施抑制了養(yǎng)殖業(yè)發(fā)展但是在促進(jìn)種養(yǎng)業(yè)協(xié)調(diào)方面卻發(fā)揮著重要的作用，這意味著環(huán)境政策在調(diào)整我國農(nóng)業(yè)空間布局和產(chǎn)業(yè)結(jié)構(gòu)方面具有一定的效力。因當(dāng)繼續(xù)執(zhí)行環(huán)境政策，倒逼種植業(yè)與養(yǎng)殖業(yè)結(jié)合并協(xié)調(diào)發(fā)展。通過制定養(yǎng)殖廢棄物還田規(guī)范、廢棄物消納農(nóng)田配套標(biāo)準(zhǔn)等，形成對養(yǎng)殖主體的指導(dǎo)與約束，在減少養(yǎng)殖污染的同時實現(xiàn)農(nóng)田有機(jī)利用。

其次，當(dāng)前種養(yǎng)分離的專業(yè)化生產(chǎn)模式不僅引發(fā)養(yǎng)殖業(yè)排污和環(huán)境污染，更是造成了種植業(yè)養(yǎng)分和肥料獲取的單一性，助長化肥過量使用，抑制土壤有機(jī)質(zhì)堆積，違背了可持續(xù)農(nóng)業(yè)的愿景。可以通過優(yōu)化種植結(jié)構(gòu)，增加農(nóng)地對本地養(yǎng)殖廢棄物的消納能力，以及飼料供應(yīng)能力。一方面實現(xiàn)種植業(yè)對養(yǎng)殖業(yè)原料供應(yīng)與環(huán)境保持的助力，另一方面促進(jìn)養(yǎng)殖業(yè)對種植業(yè)的有機(jī)質(zhì)反哺。在優(yōu)化種植結(jié)構(gòu)的過程中，增強糧食安全保障以及飼料糧供應(yīng)能力，同時減輕養(yǎng)殖業(yè)生產(chǎn)成本，形成種養(yǎng)業(yè)之間良好的互動循環(huán)，實現(xiàn)農(nóng)業(yè)結(jié)構(gòu)的優(yōu)化布局，促進(jìn)種養(yǎng)協(xié)調(diào)發(fā)展。

最后，耕地資源是確保養(yǎng)殖業(yè)環(huán)境問題得到有效解決的寶貴資源，也是保障糧食持續(xù)供應(yīng)的基礎(chǔ)條件。從多方面多角度進(jìn)行耕地保護(hù)，禁止耕地“非農(nóng)化”利用。對于違規(guī)占用耕地進(jìn)行種樹造林的地區(qū)實行退林還耕，將被毀林地等恢復(fù)成耕地并復(fù)墾種糧，增強糧食生產(chǎn)和供應(yīng)能力。通過生態(tài)補償、土地整理、高標(biāo)準(zhǔn)農(nóng)田建設(shè)等方式，提高土地利用主體耕地保護(hù)意識及能力，為種養(yǎng)協(xié)調(diào)發(fā)展鞏固土地基礎(chǔ)。