王弘燁，李翠霞

（東北農(nóng)業(yè)大學(xué)經(jīng)濟(jì)管理學(xué)院，黑龍江 哈爾濱 150030）

奶業(yè)的可持續(xù)發(fā)展是中國(guó)農(nóng)業(yè)供給側(cè)結(jié)構(gòu)性改革的重要任務(wù)和突破口，是實(shí)現(xiàn)鄉(xiāng)村振興戰(zhàn)略“產(chǎn)業(yè)興旺”目標(biāo)的重要抓手，對(duì)中國(guó)農(nóng)業(yè)的市場(chǎng)環(huán)境穩(wěn)定具有不可忽略的地位[1]。2018年6月頒布的 《國(guó)務(wù)院辦公廳關(guān)于推進(jìn)奶業(yè)振興保障乳品質(zhì)量安全的意見》確立了高質(zhì)量、綠色、可持續(xù)發(fā)展的指導(dǎo)思想，提出要全面振興奶業(yè)，提高奶業(yè)競(jìng)爭(zhēng)力，實(shí)現(xiàn)奶業(yè)現(xiàn)代化。當(dāng)下，原料乳的生產(chǎn)不但面臨著日益嚴(yán)重的資源緊缺問題，而且養(yǎng)殖產(chǎn)生的糞污造成的環(huán)境污染日益嚴(yán)峻。根據(jù)《第二次全國(guó)污染源普查公報(bào)》顯示，2017年畜禽養(yǎng)殖業(yè)排放化學(xué)需氧量1 000.53萬t，氨氮11.09萬t，總氮59.63萬t，總磷11.97萬t，是農(nóng)業(yè)污染源的首要排放源。其中，奶牛糞便日排放量顯著高于其他畜禽品種，一頭成年母牛日糞便產(chǎn)生量分別是生豬、肉羊、肉雞和蛋雞的3.57倍、13.46倍、218.57倍和116.67倍。。以上種種問題制約著奶業(yè)的高質(zhì)量、健康、可持續(xù)發(fā)展，嚴(yán)重影響奶業(yè)振興目標(biāo)的實(shí)現(xiàn)[2-4]。如何化解奶業(yè)可持續(xù)發(fā)展面臨的諸多矛盾，實(shí)現(xiàn)資源、環(huán)境、經(jīng)濟(jì)之間的協(xié)調(diào)發(fā)展，已成為政府、學(xué)術(shù)界共同關(guān)注的焦點(diǎn)問題，亟待解決。

對(duì)于原料乳生產(chǎn)空間布局的研究，學(xué)者多采用區(qū)位基尼系數(shù)、奶牛養(yǎng)殖重心模型、區(qū)域重心分析法、空間自相關(guān)等方法來揭示原料乳生產(chǎn)的區(qū)域差異[5-7]，研究范圍涉及國(guó)家、區(qū)域、省域單元。在原料乳生產(chǎn)空間布局方面：袁艷云等[5]認(rèn)為京津冀地區(qū)原料乳生產(chǎn)與乳制品加工業(yè)布局匹配度較高,但與飼料原料種植布局匹配度較低。王鈿等[6]認(rèn)為中國(guó)奶業(yè)的地理集聚程度具有一定的階段性特征，呈現(xiàn)“減弱—增強(qiáng)—減弱”的發(fā)展態(tài)勢(shì)。薛聰?shù)萚7]認(rèn)為中國(guó)奶牛養(yǎng)殖重心呈現(xiàn)出先向東南方向移動(dòng)再向西南方向移動(dòng)的演變趨勢(shì)。為保障中國(guó)奶業(yè)可持續(xù)發(fā)展，關(guān)鍵是實(shí)現(xiàn)原料乳生產(chǎn)過程中資源、產(chǎn)業(yè)、環(huán)境、市場(chǎng)與產(chǎn)能之間由低級(jí)到高級(jí)、由無序到有序的整體協(xié)作演化[8]。促進(jìn)原料乳生產(chǎn)協(xié)調(diào)發(fā)展的途徑主要有兩個(gè)：一是合理布局中國(guó)奶牛養(yǎng)殖區(qū)域。當(dāng)下，中國(guó)奶業(yè)生產(chǎn)區(qū)域間存在非均衡增長(zhǎng)特征[9]，優(yōu)勢(shì)區(qū)域主要分布在內(nèi)蒙古、河北、山東、陜西、新疆等北方省份[10]。各地區(qū)奶牛養(yǎng)殖的路徑依賴、飼料資源、財(cái)政政策、乳企數(shù)量、技術(shù)進(jìn)步等因素均會(huì)導(dǎo)致奶牛養(yǎng)殖布局變遷[11-12]。二是促進(jìn)區(qū)域內(nèi)原料乳生產(chǎn)協(xié)調(diào)發(fā)展。周鵬和樊斌[13]通過構(gòu)建乳業(yè)產(chǎn)業(yè)系統(tǒng)指標(biāo)體系，實(shí)證檢驗(yàn)得出中國(guó)乳業(yè)產(chǎn)業(yè)系統(tǒng)整體協(xié)同度水平較低。各奶業(yè)產(chǎn)區(qū)區(qū)域特征明顯，不同省份奶業(yè)發(fā)展優(yōu)勢(shì)差異較大,且原料乳生產(chǎn)受相鄰地區(qū)和本地區(qū)因素的共同影響[14-15]。

綜上所述，上述研究成果為本文研究奠定了堅(jiān)實(shí)基礎(chǔ)，提供了有益的經(jīng)驗(yàn)借鑒，但同時(shí)也存在著進(jìn)一步深化的空間：一是尚未構(gòu)建原料乳生產(chǎn)系統(tǒng)，明晰其耦合協(xié)調(diào)機(jī)理。二是尚未測(cè)度系統(tǒng)耦合協(xié)調(diào)水平，探究其空間關(guān)聯(lián)格局并分析其驅(qū)動(dòng)因素。鑒于此，本文在明晰原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)機(jī)理的基礎(chǔ)上，運(yùn)用耦合協(xié)調(diào)度模型測(cè)算中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平，分析其時(shí)空演變規(guī)律，運(yùn)用空間相關(guān)性檢驗(yàn)方法分析耦合協(xié)調(diào)發(fā)展的空間關(guān)聯(lián)格局，并進(jìn)一步采用空間杜賓模型對(duì)耦合協(xié)調(diào)發(fā)展的驅(qū)動(dòng)因素進(jìn)行分析，以期為促進(jìn)中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展提供理論支撐與經(jīng)驗(yàn)借鑒，為政府相關(guān)部門運(yùn)用差異化措施促進(jìn)各地區(qū)奶業(yè)可持續(xù)發(fā)展提供量化決策支持。

1 原料乳生產(chǎn)系統(tǒng)構(gòu)成要素及耦合協(xié)調(diào)發(fā)展機(jī)理

由系統(tǒng)論可知[16]，系統(tǒng)是若干要素以一定結(jié)構(gòu)形式聯(lián)結(jié)構(gòu)成的具有某種功能的有機(jī)整體，由經(jīng)濟(jì)地理學(xué)中的區(qū)域經(jīng)濟(jì)差異理論可知[17]，受區(qū)域經(jīng)濟(jì)發(fā)展水平、資源稟賦、生態(tài)環(huán)境和奶牛養(yǎng)殖布局等因素的影響，原料乳生產(chǎn)協(xié)調(diào)發(fā)展水平存在區(qū)域差異。本文借鑒中國(guó)科學(xué)院可持續(xù)發(fā)展研究組設(shè)計(jì)的中國(guó)可持續(xù)發(fā)展指標(biāo)體系的系統(tǒng)層劃分方式，將原料乳生產(chǎn)系統(tǒng)分為資源稟賦、產(chǎn)業(yè)基礎(chǔ)、環(huán)境抗逆、市場(chǎng)需求和原料乳產(chǎn)能等5個(gè)子系統(tǒng)，將原料乳生產(chǎn)過程看作一個(gè)以資源稟賦為前提條件，以產(chǎn)業(yè)基礎(chǔ)為輔助支撐，以環(huán)境抗逆為重要保障，以市場(chǎng)需求為動(dòng)力源泉，以原料乳產(chǎn)能為目標(biāo)導(dǎo)向構(gòu)成共生界面，耦合而成的具有供給原料乳功能的復(fù)雜有機(jī)整體，其在功能方面表現(xiàn)了資源稟賦子系統(tǒng)的承載能力、產(chǎn)業(yè)基礎(chǔ)子系統(tǒng)的推動(dòng)能力、環(huán)境抗逆子系統(tǒng)的容納能力、市場(chǎng)需求子系統(tǒng)的帶動(dòng)能力和原料乳產(chǎn)能子系統(tǒng)的回應(yīng)能力（圖1）。原料乳生產(chǎn)系統(tǒng)是一個(gè)開放的、非平衡的、具有非線性相互作用的動(dòng)態(tài)漲落系統(tǒng)，資源稟賦等5個(gè)子系統(tǒng)之間形成了相互作用、相互補(bǔ)充、共同發(fā)展的耦合共生關(guān)系。本文將5個(gè)子系統(tǒng)間相互作用、協(xié)調(diào)發(fā)展的程度界定為原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度。

圖1 原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展機(jī)理Fig. 1 Coupling and coordinated development mechanism of dairy production system

2 研究方法與指標(biāo)體系

2.1 研究方法

2.1.1 耦合協(xié)調(diào)度模型 耦合度和協(xié)調(diào)度分別反映系統(tǒng)之間相互作用的強(qiáng)弱程度和相互作用中良性耦合程度的大小，從而判斷協(xié)調(diào)狀況的好壞，可以表征各功能之間是在高水平上相互促進(jìn)還是低水平上相互制約[18]。本文用耦合協(xié)調(diào)度描述中國(guó)原料乳生產(chǎn)系統(tǒng)的各子系統(tǒng)間和諧一致發(fā)展的程度。為了耦合協(xié)調(diào)度的計(jì)算結(jié)果更科學(xué)、更準(zhǔn)確，本文采用極差法[19]對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，采用熵值法[20]測(cè)度各子系統(tǒng)各指標(biāo)權(quán)重，計(jì)算原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度過程如下。

1）計(jì)算各子系統(tǒng)的綜合發(fā)展指數(shù)：

式中：q表示原料乳生產(chǎn)系統(tǒng)子系統(tǒng)個(gè)數(shù)；m表示指標(biāo)個(gè)數(shù)，ij表示第i個(gè)地區(qū)的第j個(gè)指標(biāo)，xij表示標(biāo)準(zhǔn)化處理后的數(shù)值為各子系統(tǒng)中各指標(biāo)的權(quán)重，為各子系統(tǒng)綜合發(fā)展指數(shù)，Uq越大，表示該子系統(tǒng)綜合發(fā)展水平越高[21]。

2）本文采用容量耦合模型[22]計(jì)算耦合度：

3）構(gòu)建耦合協(xié)調(diào)度模型[23]，計(jì)算耦合協(xié)調(diào)度：

式中：D為原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度，D∈[0,1]；C為式（2）計(jì)算得到的耦合度，C∈[0,1]；T代表原料乳生產(chǎn)系統(tǒng)的綜合發(fā)展指數(shù)，T∈[0,1]；wq代表各子系統(tǒng)權(quán)重，反映了各子系統(tǒng)的重要程度。借鑒已有研究成果[24]，將耦合協(xié)調(diào)度劃分為3個(gè)階段10個(gè)類型（表1）。

表1 耦合協(xié)調(diào)度評(píng)價(jià)標(biāo)準(zhǔn)和類型Table 1 Evaluation criteria and types of the coupling coordination

2.1.2 空間相關(guān)性檢驗(yàn) 本文采用Moran’sI[25]刻畫原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度空間分布與集聚特征，其計(jì)算公式如下：

式中：Moran’sI∈[-1, 1]，當(dāng)Moran’sI為正時(shí)觀測(cè)值呈現(xiàn)空間正相關(guān)性，反之為負(fù)相關(guān)性。Moran’sI絕對(duì)值越大，表明空間相關(guān)性越強(qiáng)。若Moran’sI趨近于0，則表示數(shù)據(jù)之間相互獨(dú)立，在空間上隨機(jī)分布。xi和xj分別表示區(qū)域i、區(qū)域j地理單元上的觀測(cè)值，Wij為地理單元相互之間0-1鄰接矩陣的權(quán)重矩陣，其中，i和j分別表示空間單元編號(hào)，即中國(guó)31個(gè)省份，0-1空間權(quán)重矩陣的對(duì)角線為0，假定同單元之間不相鄰，當(dāng)Wij=1表示相鄰，當(dāng)Wij=0表示不相鄰，和S2分別為均值和平方差。

通過構(gòu)造正態(tài)統(tǒng)計(jì)量Z檢驗(yàn)Moran’sI顯著性水平：

式中：E(I)為Moran’sI的期望值，Var(I)為Moran’sI的方差。

2.1.3 空間杜賓模型 區(qū)域農(nóng)村人力資本水平、市場(chǎng)化水平等因素或?qū)υ先樯a(chǎn)系統(tǒng)耦合協(xié)調(diào)關(guān)系產(chǎn)生影響，但同一因素在不同地區(qū)可能會(huì)表現(xiàn)出差異化的影響效果，為了考慮這種空間異質(zhì)性問題，本文采用空間杜賓模型[26]（Spatial Durbin Model，SDM）進(jìn)一步深入分析，其一般模型形式如下：

式中：Y表示因變量；ρ表示空間回歸系數(shù)；W為空間鄰接權(quán)重矩陣；WY表示因變量的空間滯后項(xiàng)；X為解釋變量；β表示各解釋變量的回歸系數(shù)；WX表示解釋變量的空間滯后項(xiàng)；δ為解釋變量空間滯后系數(shù)；ε為隨機(jī)誤差向量?？臻g杜賓模型同時(shí)考察驅(qū)動(dòng)因素對(duì)本地區(qū)和鄰近地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的影響。

2.2 指標(biāo)選取

2.2.1 原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度評(píng)價(jià)指標(biāo)體系 原料乳生產(chǎn)耦合協(xié)調(diào)度評(píng)價(jià)指標(biāo)體系既是一個(gè)多層次多因素的復(fù)雜體系，又是綜合反映各子系統(tǒng)發(fā)展?fàn)顩r的指標(biāo)集合體。本文基于《全國(guó)奶業(yè)發(fā)展規(guī)劃（2016—2020）》的總體思路、基本原則和發(fā)展目標(biāo)，考慮奶牛養(yǎng)殖分布格局的影響因素，并借鑒吳強(qiáng)等[27]構(gòu)建的生豬生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度的評(píng)價(jià)指標(biāo)體系，選擇資源稟賦、產(chǎn)業(yè)基礎(chǔ)、環(huán)境抗逆、市場(chǎng)需求和原料乳產(chǎn)能5個(gè)方面，共16個(gè)指標(biāo)組成原料乳生產(chǎn)耦合協(xié)調(diào)度評(píng)價(jià)指標(biāo)體系。其中，環(huán)境規(guī)制強(qiáng)度借鑒周建軍等[28]研究成果，用規(guī)制系數(shù)反映，其計(jì)算公式為：其中A表示各省域面積，π代表圓周率。指標(biāo)的具體含義、正逆向性等具體見表2。

表2 原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度評(píng)價(jià)指標(biāo)Table 2 Evaluation indexes of the coupling degree of the coordination of raw milk production system

2.2.2 驅(qū)動(dòng)因素選取 原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展受多種因素的共同影響和制約，本文基于資源環(huán)境經(jīng)濟(jì)學(xué)、生態(tài)經(jīng)濟(jì)學(xué)和發(fā)展經(jīng)濟(jì)學(xué)的理論分析，具體驅(qū)動(dòng)因素指標(biāo)選擇見表3。

表3 原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的驅(qū)動(dòng)因素Table 3 Driving factors of the coupled and coordinated development of raw milk production system

1）財(cái)政支農(nóng)。政府的財(cái)政支持是奶牛良種繁育、奶牛養(yǎng)殖技術(shù)升級(jí)的重要資金來源。政府出臺(tái)的奶牛養(yǎng)殖補(bǔ)貼政策有利于促進(jìn)良種的普及和養(yǎng)殖規(guī)模的擴(kuò)張，進(jìn)而影響原料乳生產(chǎn)效率。

2）對(duì)外開放水平。隨著對(duì)外開放程度不斷提高，由于中國(guó)資源稟賦等存在劣勢(shì)，國(guó)內(nèi)外奶價(jià)存在差異，驅(qū)動(dòng)乳制品大量進(jìn)口，對(duì)各地區(qū)原料乳生產(chǎn)造成極大沖擊[29]。

3）農(nóng)村人力資本水平。農(nóng)村人力資本水平的提升能促進(jìn)知識(shí)等高端要素流動(dòng)，實(shí)現(xiàn)奶牛養(yǎng)殖的技術(shù)進(jìn)步[30]，是奶牛養(yǎng)殖業(yè)向更高層級(jí)演進(jìn)的關(guān)鍵因素，本文用農(nóng)村居民人均受教育年限衡量，各地區(qū)農(nóng)村勞動(dòng)力不同受教育程度的人口比重乘以對(duì)應(yīng)的受教育年限,文盲或者半文盲為0年、小學(xué)為6年、初中為9年、高中及中專為12年、大專及以上為16年。

4）市場(chǎng)化水平。市場(chǎng)環(huán)境是技術(shù)創(chuàng)新有效發(fā)生的基礎(chǔ)性條件，市場(chǎng)化進(jìn)程的加快有利于發(fā)揮市場(chǎng)機(jī)制對(duì)要素資源配置的決定性作用[31]，促進(jìn)要素向更高生產(chǎn)率、更高創(chuàng)新水平的企業(yè)及產(chǎn)業(yè)流動(dòng)，升級(jí)企業(yè)及產(chǎn)業(yè)的創(chuàng)新激勵(lì)機(jī)制，推動(dòng)技術(shù)創(chuàng)新，進(jìn)而實(shí)現(xiàn)資源有效配置，降低環(huán)境污染程度，減少產(chǎn)能過剩，推動(dòng)經(jīng)濟(jì)綠色增長(zhǎng)。

5）科技創(chuàng)新水平。科技創(chuàng)新可以推動(dòng)奶牛養(yǎng)殖業(yè)的技術(shù)進(jìn)步，從而實(shí)現(xiàn)原料乳生產(chǎn)效率的提升。不同地區(qū)的科技創(chuàng)新水平存在差異性會(huì)導(dǎo)致原料乳生產(chǎn)效率的不同，且會(huì)引起原料乳生產(chǎn)由效率高的地區(qū)向效率低的地區(qū)緩慢轉(zhuǎn)移，最終達(dá)到均衡狀態(tài)。

6）產(chǎn)業(yè)集聚度。產(chǎn)業(yè)集聚通過實(shí)現(xiàn)內(nèi)、外部規(guī)模經(jīng)濟(jì)影響原料乳生產(chǎn)效率。最佳規(guī)模理論認(rèn)為產(chǎn)業(yè)集聚和原料乳生產(chǎn)效率之間是倒“U”型關(guān)系。

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

基于數(shù)據(jù)的價(jià)值性、可得性以及研究時(shí)段的完整性，本文選取中國(guó)（除港、澳、臺(tái)地區(qū)外）31省（市、區(qū)）作為研究樣本，選取2010—2019年作為研究時(shí)段。其中飼料生產(chǎn)能力、防疫檢驗(yàn)?zāi)芰?、道路通達(dá)性、環(huán)境規(guī)制強(qiáng)度、消費(fèi)市場(chǎng)潛力、城鎮(zhèn)化水平、城鎮(zhèn)消費(fèi)市場(chǎng)能力、鄉(xiāng)村消費(fèi)市場(chǎng)能力、科技創(chuàng)新水平、對(duì)外開放水平指標(biāo)的測(cè)算數(shù)據(jù)來源于《中國(guó)統(tǒng)計(jì)年鑒》，奶牛布局指數(shù)、規(guī)模飼養(yǎng)能力、技術(shù)水平、乳制品加工能力、原料乳生產(chǎn)能力指標(biāo)的測(cè)算數(shù)據(jù)來源于《中國(guó)奶業(yè)統(tǒng)計(jì)資料》，污染治理投資能力、資源化利用能力指標(biāo)的測(cè)算數(shù)據(jù)來源于《中國(guó)環(huán)境統(tǒng)計(jì)年鑒》，農(nóng)業(yè)機(jī)械化水平、財(cái)政支農(nóng)、農(nóng)村人力資本水平、奶業(yè)區(qū)位熵指標(biāo)的測(cè)算數(shù)據(jù)來源于《中國(guó)農(nóng)村統(tǒng)計(jì)年鑒》，市場(chǎng)化水平的測(cè)算數(shù)據(jù)來源于《中國(guó)市場(chǎng)化指數(shù)數(shù)據(jù)庫》。

3 結(jié)果與分析

3.1 中國(guó)原料乳生產(chǎn)系統(tǒng)滯后子系統(tǒng)動(dòng)態(tài)變化

選取樣本期的初始時(shí)間（2010年）、中間時(shí)間（2013年、2016年）和結(jié)束時(shí)間（2019年）為節(jié)點(diǎn)，從原料乳生產(chǎn)系統(tǒng)滯后子系統(tǒng)的動(dòng)態(tài)變化來看（表4），資源稟賦子系統(tǒng)滯后省份的數(shù)量變化為4→4→4→3，產(chǎn)業(yè)基礎(chǔ)子系統(tǒng)滯后省份的數(shù)量變化為2→3→2→1，環(huán)境抗逆子系統(tǒng)滯后省份的數(shù)量變化為8→6→7→6，市場(chǎng)需求子系統(tǒng)滯后省份的數(shù)量變化為0→0→0→0，原料乳產(chǎn)能子系統(tǒng)滯后省份的數(shù)量變化為17→18→18→21。根據(jù)滯后子系統(tǒng)動(dòng)態(tài)變化可知，中國(guó)乳制品市場(chǎng)需求整體向好，資源稟賦滯后呈縮小趨勢(shì)，產(chǎn)業(yè)基礎(chǔ)較為鞏固，環(huán)境抗逆能力緩慢增強(qiáng)，原料乳產(chǎn)能亟待提升。由此可見，加大污染防治力度是促進(jìn)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的有力抓手，不斷提升原料乳產(chǎn)能是實(shí)現(xiàn)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的重中之重。

表4 2010、2013、2016和2019年原料乳生產(chǎn)系統(tǒng)滯后子系統(tǒng)Table 4 Raw milk production system lagging subsystem in 2010, 2013, 2016 and 2019

3.2 中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的時(shí)空演變特征

3.2.1 時(shí)序演變分析 為了判斷中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的時(shí)空演變特征，根據(jù)中國(guó)31個(gè)?。ㄊ?、區(qū)）的耦合協(xié)調(diào)度，運(yùn)用Arcgis 10.8軟件，將2010年、2013年、2016年、2019年中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度進(jìn)行空間可視化（圖2）。

圖2 中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度的時(shí)空演變過程圖Fig. 2 Spatial and temporal evolution of the coupling coordination of China’s raw milk production system

2010—2019年原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度呈上漲趨勢(shì)。其中，西藏、寧夏上升最為明顯，但導(dǎo)致兩個(gè)地區(qū)上升的原因不同，西藏是由于產(chǎn)業(yè)基礎(chǔ)不斷發(fā)展的帶動(dòng)，寧夏則由于原料乳產(chǎn)能快速提升的促進(jìn)，使得兩地區(qū)逐步優(yōu)化子系統(tǒng)匹配程度。僅有北京、海南、江西3個(gè)地區(qū)系統(tǒng)耦合協(xié)調(diào)度明顯下降，分析發(fā)現(xiàn)這些地區(qū)資源稟賦供給趨弱，原料乳產(chǎn)能下降與市場(chǎng)需求的快速提升不匹配，導(dǎo)致耦合協(xié)調(diào)度降低。

從各地區(qū)平均耦合協(xié)調(diào)度來看，分別有9.68%、22.58%、67.74%的地區(qū)處于協(xié)調(diào)發(fā)展階段、協(xié)調(diào)過渡階段、失調(diào)衰退階段。耦合協(xié)調(diào)度較高的地區(qū)主要為山東、河北、內(nèi)蒙古等原料乳主產(chǎn)區(qū)，這些地區(qū)自然資源稟賦占優(yōu)，奶業(yè)產(chǎn)業(yè)基礎(chǔ)堅(jiān)實(shí)。耦合協(xié)調(diào)度較低的地區(qū)主要為北京、上海等經(jīng)濟(jì)發(fā)達(dá)地區(qū)及西藏、青海、貴州等西部地區(qū)。經(jīng)濟(jì)發(fā)達(dá)地區(qū)的原料乳生產(chǎn)系統(tǒng)失調(diào)主要受當(dāng)?shù)刳吘o的環(huán)保政策和第三產(chǎn)業(yè)為主的經(jīng)濟(jì)結(jié)構(gòu)影響，奶業(yè)發(fā)展受到政策因素等多方面的限制；西部地區(qū)系統(tǒng)失調(diào)的原因主要來自產(chǎn)業(yè)基礎(chǔ)薄弱，無力滿足奶業(yè)發(fā)展現(xiàn)實(shí)需求。3.2.2 空間特征分析 為更加直觀的看出2010—2019年原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的時(shí)空演變特征。本文運(yùn)用Arcgis 10.8軟件，分別繪制出2010年、2013年、2016年、2019年原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度的空間趨勢(shì)分析圖（圖3）。其中，X軸指向東方，Y軸指向北方，分別代表了系統(tǒng)耦合協(xié)調(diào)度在中國(guó)東西（綠色曲線）、南北（藍(lán)色曲線）方向的投影。

圖3 中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度的空間趨勢(shì)圖Fig. 3 Spatial trend of the coupling coordination of China’s raw milk production system

從東西方向上投影看，趨勢(shì)線由“從西向東平滑上升”的直線型，逐步演變?yōu)椤皷|西高，中部低”的“凹”型，趨勢(shì)線愈發(fā)平緩，說明東西方向各地區(qū)之間系統(tǒng)耦合協(xié)調(diào)度差距逐漸縮小，從南北方向投影看，中部偏北持續(xù)抬升，趨勢(shì)線由“從南向北平滑上升”趨勢(shì)逐步演變?yōu)椤巴埂毙?。趨?shì)線愈發(fā)陡峭，說明南北方向上各地區(qū)之間系統(tǒng)耦合協(xié)調(diào)度差距逐漸擴(kuò)大。南北方向趨勢(shì)線較東西方向始終更為陡峭，說明系統(tǒng)耦合協(xié)調(diào)度的南北方向差異始終大于東西方向差異。

3.3 中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的空間關(guān)聯(lián)格局

3.3.1 空間自相關(guān)分析 盡管各地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度有明顯差異，但仍顯露出區(qū)域性耦合協(xié)調(diào)度接近的現(xiàn)象，如冀豫魯、黑吉遼、粵桂湘等區(qū)域內(nèi)耦合協(xié)調(diào)度差異較小。為此，本文采用Moran’sI測(cè)算耦合協(xié)調(diào)度的空間自相關(guān)性。由表5可知，2010—2019年Moran’sI位于0.290 8-0.349 7之間,呈現(xiàn)“M”型趨勢(shì)，且在1%的水平上通過了顯著性檢驗(yàn)，表明中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度存在正向的空間自相關(guān)性，即相鄰地區(qū)的系統(tǒng)耦合協(xié)調(diào)度表現(xiàn)出空間集聚特征。各地區(qū)原料乳生產(chǎn)系統(tǒng)的耦合協(xié)調(diào)度不僅與各子系統(tǒng)匹配程度有關(guān)，還與周邊鄰近地區(qū)的原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展?fàn)顩r存在正相關(guān)性。

表5 中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度的Moran’s I及其顯著性檢驗(yàn)Table 5 Moran’s I and its significance test for the coupling coordination of China’s raw milk production system

3.3.2 空間集聚特征 為進(jìn)一步詳細(xì)考察原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度的空間集聚特征，本文運(yùn)用MATLAB R2022a軟件繪制了原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度的Moran’sI散點(diǎn)圖（圖4）。結(jié)果顯示：位于Ⅰ、Ⅲ象限的地區(qū)占大部分，在2010、2013、2016和2019年分別占總地區(qū)的74.19%、70.97%、70.97%和70.97%，呈現(xiàn)高—高集聚和低—低集聚特征，表明原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度存在顯著的正向空間自相關(guān)性，同質(zhì)性溢出效應(yīng)明顯。

圖4 中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度的Moran’s I散點(diǎn)圖Fig. 4 Moran’s I scatter plot of the coupling coordination of China’s raw milk production system

3.4 中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的驅(qū)動(dòng)因素分析

3.4.1 空間效應(yīng)診斷檢驗(yàn) 中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度具有較強(qiáng)的空間相關(guān)性特征，因此研究原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的驅(qū)動(dòng)因素，必須要考慮空間因素?？臻g面板計(jì)量模型嵌套了空間效應(yīng)和時(shí)間效應(yīng)，可以盡量避免因時(shí)空因素的存在而使回歸結(jié)果出現(xiàn)偏差。本文運(yùn)用Stata16.0軟件對(duì)空間計(jì)量模型進(jìn)行識(shí)別檢驗(yàn)（表6）。首先，Robust LM統(tǒng)計(jì)量均通過了1%水平的顯著性檢驗(yàn),表明應(yīng)該選擇空間杜賓模型（SDM）來研究。其次，Wald和LR統(tǒng)計(jì)量均通過了1%水平的顯著性檢驗(yàn)，表明SDM不能退化為SAR和SEM。再次，Hausman統(tǒng)計(jì)量通過1%的顯著性檢驗(yàn)，表明固定效應(yīng)較為合適。最后，SDM的固定效應(yīng)包括空間固定效應(yīng)、時(shí)間固定效應(yīng)和時(shí)空雙固定效應(yīng)。時(shí)間和空間固定效應(yīng)聯(lián)合顯著性的LR統(tǒng)計(jì)量通過了1%水平的顯著性檢驗(yàn)，表明時(shí)空雙固定效應(yīng)優(yōu)于時(shí)間固定效應(yīng)和空間固定效應(yīng)。因此，本文選擇時(shí)空雙固定效應(yīng)的SDM模型來測(cè)度原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展驅(qū)動(dòng)因素的影響。

表6 空間計(jì)量模型識(shí)別檢驗(yàn)結(jié)果Table 6 Spatial econometric model identification test results

3.4.2 空間溢出效應(yīng)測(cè)度與分析 對(duì)耦合協(xié)調(diào)發(fā)展驅(qū)動(dòng)因素的面板數(shù)據(jù)進(jìn)行時(shí)空雙固定效應(yīng)的SDM估計(jì)結(jié)果顯示（表7）：從回歸系數(shù)來看，對(duì)外開放水平（O）、農(nóng)村人力資本水平（E）、市場(chǎng)化水平（M）、科技創(chuàng)新水平（S）、產(chǎn)業(yè)集聚度（L）五個(gè)變量均通過顯著性檢驗(yàn)，且符號(hào)均為正，表明上述5個(gè)變量對(duì)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展均有顯著正向影響。從解釋變量的空間滯后項(xiàng)看，W×農(nóng)村人力資本水平、W×市場(chǎng)化水平分別在1%、5%的水平上通過顯著性檢驗(yàn)，且符號(hào)為正；W×產(chǎn)業(yè)集聚在5%的水平上通過顯著性檢驗(yàn)，且符號(hào)為負(fù)，表明農(nóng)村居民人均受教育年限越長(zhǎng)、市場(chǎng)化水平越高的地區(qū)，對(duì)周邊鄰近地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的帶動(dòng)作用越強(qiáng)。產(chǎn)業(yè)集聚度越高對(duì)周邊鄰近地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的阻礙作用越明顯。

表7 時(shí)空固定效應(yīng)的面板SDM估計(jì)結(jié)果Table 7 Panel SDM spatial effects for the spatiotemporal fixed effects

3.4.3 空間溢出效應(yīng)分解 由于空間自回歸系數(shù)ρ在1%的水平上顯著為正，模型中各解釋變量的估計(jì)系數(shù)值并不能直接解釋該變量對(duì)因變量的影響效應(yīng)。因此，需要采用偏微分方法對(duì)空間效應(yīng)進(jìn)行分解[32]，以便更好地解釋各個(gè)驅(qū)動(dòng)因素的直接效應(yīng)、間接效應(yīng)，以及對(duì)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展影響的總效應(yīng)（表8）。

表8 時(shí)空固定效應(yīng)的面板SDM空間效應(yīng)分解Table 8 Decomposition of the panel SDM spatial effects for the spatiotemporal fixed effects

1）直接效應(yīng)。直接效應(yīng)是驅(qū)動(dòng)因素對(duì)本地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平的影響。對(duì)外開放水平（O）、農(nóng)村人力資本水平（E）、市場(chǎng)化水平（M）、科技創(chuàng)新水平（S）、產(chǎn)業(yè)集聚度（L）的彈性系數(shù)為正，均通過顯著性檢驗(yàn)。對(duì)外開放水平的提升可以引入國(guó)外先進(jìn)技術(shù)，刺激本土企業(yè)開展研發(fā)，同時(shí)進(jìn)出口貿(mào)易的發(fā)展在短期內(nèi)補(bǔ)足了本地區(qū)原料乳產(chǎn)能不足的困境，優(yōu)化原料乳產(chǎn)能子系統(tǒng)，進(jìn)而提升本地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平。農(nóng)村人力資本累積對(duì)技術(shù)的發(fā)明創(chuàng)造及引進(jìn)、消化、吸收起促進(jìn)作用，高素質(zhì)、專業(yè)型人才能夠有效促進(jìn)前沿技術(shù)轉(zhuǎn)化為現(xiàn)實(shí)生產(chǎn)力，促進(jìn)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展。科技創(chuàng)新水平是推動(dòng)奶業(yè)結(jié)構(gòu)調(diào)整升級(jí)的內(nèi)生動(dòng)力，在供給端，通過更新原料乳生產(chǎn)技術(shù)，促進(jìn)勞動(dòng)、資本生產(chǎn)率的提高；在需求端，通過延伸乳制品品類，創(chuàng)造新的乳制品消費(fèi)需求，推動(dòng)乳制品產(chǎn)業(yè)發(fā)展，是促進(jìn)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展的重要因素。市場(chǎng)化水平的彈性系數(shù)最大，說明市場(chǎng)化水平相較于其他驅(qū)動(dòng)因素更具有直接影響力，一方面市場(chǎng)化水平的提高是一種制度紅利，市場(chǎng)化水平提升有助于促進(jìn)要素市場(chǎng)成熟，提升市場(chǎng)資源配置效率，進(jìn)而推動(dòng)本地區(qū)的奶業(yè)結(jié)構(gòu)調(diào)整。另一方面市場(chǎng)化水平高的地區(qū)，企業(yè)和政府的利益共同體聯(lián)系強(qiáng)度減弱，企業(yè)不得不在高度競(jìng)爭(zhēng)的市場(chǎng)中不斷強(qiáng)化自身創(chuàng)新能力，進(jìn)而提升原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平。產(chǎn)業(yè)集聚可以提升本地區(qū)奶業(yè)專業(yè)化水平，使生產(chǎn)資料交易費(fèi)用降低，提高原料乳生產(chǎn)效率，進(jìn)而形成省際差異。通過產(chǎn)業(yè)鏈各環(huán)節(jié)組織程度的集中，引領(lǐng)生產(chǎn)要素在本地區(qū)內(nèi)集聚，正向影響本地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平。財(cái)政支農(nóng)（F）對(duì)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度影響為正，但并不顯著。盡管財(cái)政支農(nóng)力度越來越大，但在低違規(guī)違約成本情況下，部分養(yǎng)殖主體拿到資金后繼續(xù)從事高投入、高污染、低產(chǎn)出的原料乳生產(chǎn)行為，造成了財(cái)政資金低效率，所以財(cái)政支農(nóng)對(duì)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平的影響不顯著。

2）間接效應(yīng)。間接效應(yīng)反映的是驅(qū)動(dòng)因素對(duì)鄰近地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平的影響，也稱之為溢出效應(yīng)。農(nóng)村人力資本水平（E）、市場(chǎng)化水平（M）的彈性系數(shù)為正，在1%水平上顯著。隨著農(nóng)業(yè)人力資本的提高，流動(dòng)規(guī)模逐漸增加，通過知識(shí)溢出，農(nóng)業(yè)人力資本在一定區(qū)域內(nèi)形成“1+1＞2”的集聚效應(yīng)，降低了省際間奶業(yè)信息傳播與知識(shí)交流成本，強(qiáng)化了局部區(qū)域技術(shù)創(chuàng)新意識(shí)，促進(jìn)了鄰近地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展。市場(chǎng)化在一定區(qū)域內(nèi)能充分地調(diào)動(dòng)生產(chǎn)要素，實(shí)現(xiàn)資源合理配置并提供資金支持，市場(chǎng)化水平越高，開放性、競(jìng)爭(zhēng)性特征越顯著，越能帶動(dòng)鄰近地區(qū)耦合協(xié)調(diào)水平的提升。產(chǎn)業(yè)集聚度（L）的彈性系數(shù)為-0.022，在5%水平上顯著。說明產(chǎn)業(yè)集聚度能阻礙鄰近地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展。由于產(chǎn)業(yè)集聚引領(lǐng)生產(chǎn)要素集聚，表現(xiàn)出空間虹吸效應(yīng)，形成局部區(qū)域內(nèi)奶業(yè)增長(zhǎng)的中心點(diǎn)，導(dǎo)致鄰近地區(qū)子系統(tǒng)間失衡，抑制鄰近地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平的提升。而本地的財(cái)政支農(nóng)、對(duì)外開放水平、科技創(chuàng)新水平對(duì)鄰近地區(qū)系統(tǒng)耦合協(xié)調(diào)水平無顯著影響。

3）總效應(yīng)?？傂?yīng)反映的是直接效應(yīng)和間接效應(yīng)對(duì)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平的綜合影響程度。農(nóng)村人力資本水平（E）、市場(chǎng)化水平（M）的彈性系數(shù)為正，在1%水平上顯著。未來應(yīng)該繼續(xù)推進(jìn)市場(chǎng)化進(jìn)程，激發(fā)市場(chǎng)主體活力，充分釋放市場(chǎng)潛力。同時(shí)發(fā)揮人力資源優(yōu)勢(shì)，加強(qiáng)農(nóng)村實(shí)用人才隊(duì)伍建設(shè)，推動(dòng)奶業(yè)可持續(xù)發(fā)展。此外，財(cái)政支農(nóng)、對(duì)外開放水平、科技創(chuàng)新水平、產(chǎn)業(yè)集聚度對(duì)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平的總體影響力不高，總效應(yīng)未通過顯著性檢驗(yàn)。

4 結(jié)論與對(duì)策建議

4.1 結(jié)論

本文基于2010—2019年中國(guó)31?。ㄊ小^(qū)）面板數(shù)據(jù)，以原料乳生產(chǎn)系統(tǒng)為研究對(duì)象，采用耦合協(xié)調(diào)度模型、空間相關(guān)性檢驗(yàn)方法以及空間杜賓模型對(duì)中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平進(jìn)行測(cè)度，并探究其時(shí)空演變特征及驅(qū)動(dòng)因素，以期為中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平的提高提出有效建議。研究結(jié)果表明：

1）在時(shí)間維度上，中國(guó)各地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度逐年穩(wěn)定上升，大部分地區(qū)處于失調(diào)衰退階段和協(xié)調(diào)過渡階段，市場(chǎng)需求子系統(tǒng)與原料乳產(chǎn)能子系統(tǒng)為主要滯后子系統(tǒng)。

2）在空間維度上，原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)度呈現(xiàn)“東高西低，北高南低”態(tài)勢(shì)，南北方向差異逐漸擴(kuò)大，東西方向差異逐漸減少。系統(tǒng)耦合協(xié)調(diào)度存在顯著的正向空間自相關(guān)性，呈現(xiàn)高—高集聚和低—低集聚特征，同質(zhì)溢出性效應(yīng)明顯。

3）空間杜賓模型檢驗(yàn)與估計(jì)結(jié)果表明：直接效應(yīng)層面，對(duì)外開放水平、農(nóng)村人力資本水平、市場(chǎng)化水平、科技創(chuàng)新水平、產(chǎn)業(yè)集聚度對(duì)地區(qū)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展具有顯著促進(jìn)作用；間接效應(yīng)層面，農(nóng)村人力資本水平和市場(chǎng)化水平具有顯著空間正溢出效應(yīng)，產(chǎn)業(yè)集聚度具有顯著空間負(fù)溢出效應(yīng)；總效應(yīng)層面，農(nóng)村人力資本水平和市場(chǎng)化水平對(duì)中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)發(fā)展具有顯著促進(jìn)作用。

4.2 對(duì)策建議

根據(jù)研究結(jié)果，中國(guó)原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平存在顯著的空間分異形態(tài)和空間溢出現(xiàn)象，因此在制定和實(shí)施奶業(yè)發(fā)展相關(guān)規(guī)制政策時(shí)需注意以下方面：

1）原料乳主產(chǎn)省份應(yīng)以降本、增效、綠色發(fā)展為原則，促進(jìn)種養(yǎng)結(jié)合、草畜配套。加強(qiáng)奶業(yè)綠色發(fā)展資金保障和監(jiān)管，落實(shí)生態(tài)保護(hù)補(bǔ)償政策，推動(dòng)奶業(yè)可持續(xù)發(fā)展。建立健全現(xiàn)代奶業(yè)生產(chǎn)經(jīng)營(yíng)體系，支持有條件的奶農(nóng)發(fā)展乳制品加工，提高奶業(yè)增值收益和抗市場(chǎng)風(fēng)險(xiǎn)能力。破除省份間市場(chǎng)壁壘，促進(jìn)生產(chǎn)要素自由流動(dòng)，輻射周邊省份，帶動(dòng)?xùn)|部、西部其他省份協(xié)同發(fā)展。

2）東部發(fā)達(dá)省份應(yīng)堅(jiān)持市場(chǎng)導(dǎo)向，持續(xù)推進(jìn)市場(chǎng)化進(jìn)程。鼓勵(lì)企業(yè)進(jìn)行獨(dú)立研究，推動(dòng)企業(yè)“走出去”，通過引進(jìn)、消化、再吸收，推進(jìn)奶牛養(yǎng)殖向信息化、智能化、數(shù)字化邁進(jìn)，提高奶業(yè)綜合競(jìng)爭(zhēng)力。同時(shí)積極鼓勵(lì)優(yōu)秀畢業(yè)生、具有工作技術(shù)經(jīng)驗(yàn)的農(nóng)民工和科技人員返鄉(xiāng)進(jìn)行創(chuàng)業(yè)創(chuàng)新活動(dòng)，逐步提升農(nóng)村人力資本對(duì)鄰近省份原料乳生產(chǎn)系統(tǒng)耦合協(xié)調(diào)水平的空間溢出效應(yīng)，整體提高奶業(yè)發(fā)展水平。

3）西部生態(tài)資源豐富的省份，應(yīng)充分利用自然資源優(yōu)勢(shì)，筑牢產(chǎn)業(yè)基礎(chǔ)，著眼于利益分配合理化和提升奶農(nóng)地位，推動(dòng)奶業(yè)產(chǎn)業(yè)鏈調(diào)整完善。強(qiáng)化技術(shù)服務(wù)，建立專家指導(dǎo)服務(wù)制度，加強(qiáng)對(duì)項(xiàng)目的技術(shù)指導(dǎo)服務(wù)，推動(dòng)建立以企業(yè)為主體、科研院所為支撐、產(chǎn)學(xué)研結(jié)合的奶業(yè)科技創(chuàng)新體系。并進(jìn)一步深化西部對(duì)外開放，利用一帶一路建設(shè)工程，提升與其他國(guó)家貿(mào)易合作水平，推動(dòng)西南和西北原料乳生產(chǎn)系統(tǒng)協(xié)調(diào)發(fā)展，形成奶業(yè)可持續(xù)發(fā)展的長(zhǎng)期動(dòng)力。