劉 帥，張航宇，蔡文靜

(1.中國(guó)社會(huì)科學(xué)院財(cái)經(jīng)戰(zhàn)略研究院，北京 100006；2.北京大學(xué)中國(guó)農(nóng)業(yè)政策研究中心，北京 100871；3.滄州市文化廣電和旅游局，河北 滄州 061011)

黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展已上升為國(guó)家戰(zhàn)略。黃河流域流經(jīng)中國(guó)9個(gè)省份，是中國(guó)第二長(zhǎng)河，連接青藏高原、黃土高原和華北平原，流域內(nèi)煤炭、石油等資源豐富，既是重要的經(jīng)濟(jì)地帶，也是重要的生態(tài)屏障。黃河流域各省份2019年總?cè)丝?.22億，占全國(guó)總?cè)丝诘谋壤^30%；國(guó)內(nèi)生產(chǎn)總值24.74萬億，占全國(guó)的25.1%。黃河流域是中國(guó)重要的農(nóng)產(chǎn)品產(chǎn)區(qū)，糧食和肉類產(chǎn)量占全國(guó)的1/3 左右。黃河流域各省份2019年第一產(chǎn)業(yè)產(chǎn)值占全國(guó)的29.6%，超過第二產(chǎn)業(yè)(26.2%)和第三產(chǎn)業(yè)占比(23.7%)。黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展是習(xí)近平總書記親自謀劃、親自部署、親自推動(dòng)的重大國(guó)家戰(zhàn)略，已納入“十四五”規(guī)劃和二○三五年遠(yuǎn)景目標(biāo)。黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展也是貫徹綠色發(fā)展理念和“綠水青山就是金山銀山”理念的重要實(shí)踐載體。

農(nóng)業(yè)綠色全要素生產(chǎn)率(total factor productivity，TFP)指在傳統(tǒng)TFP中考慮環(huán)境污染等非期望產(chǎn)出[1]，統(tǒng)籌兼顧農(nóng)業(yè)發(fā)展、資源節(jié)約和環(huán)境保護(hù)，將其納入一個(gè)統(tǒng)一分析框架[2]。國(guó)內(nèi)外對(duì)傳統(tǒng)TFP的研究比較豐富，包括理論內(nèi)涵[3]、測(cè)算方法[4-5]、實(shí)證測(cè)算[6]、影響因素[7]、時(shí)空特征[8]等。在農(nóng)業(yè)領(lǐng)域，國(guó)內(nèi)外學(xué)者對(duì)農(nóng)業(yè)TFP也進(jìn)行了深入研究。例如，ANG等[9]測(cè)算了美國(guó)農(nóng)業(yè)TFP，并對(duì)其進(jìn)行了分解；王璐等[10]利用微觀數(shù)據(jù)計(jì)算中國(guó)農(nóng)戶農(nóng)業(yè)TFP，并分析了其結(jié)構(gòu)變遷問題；張樂等[11]、李欠男等[12]則從宏觀角度測(cè)算了中國(guó)或省際農(nóng)業(yè)TFP。國(guó)內(nèi)外學(xué)者對(duì)農(nóng)業(yè)綠色TFP也進(jìn)行了有益探索，測(cè)算農(nóng)業(yè)綠色TFP的關(guān)鍵在于指標(biāo)體系和方法，筆者對(duì)國(guó)內(nèi)外相關(guān)文獻(xiàn)進(jìn)行了梳理，部分代表作如表1所示。

從評(píng)價(jià)指標(biāo)體系上來看，投入方面，多數(shù)研究涵蓋了各類農(nóng)業(yè)生產(chǎn)的投入要素，但也存在較大差異，如能源、用水、資本等投入指標(biāo)并沒有得到廣泛的使用，一些研究可能受限于研究條件等，甚至簡(jiǎn)化了投入指標(biāo)。產(chǎn)出方面，期望產(chǎn)出主要有廣義農(nóng)業(yè)(即農(nóng)、林、牧、漁業(yè))和狹義農(nóng)業(yè)2類，也有針對(duì)特定農(nóng)業(yè)子行業(yè)的研究。非期望產(chǎn)出主要有農(nóng)業(yè)碳排放和農(nóng)業(yè)面源污染2種。整體上看，農(nóng)業(yè)綠色TFP測(cè)算指標(biāo)體系的構(gòu)建主觀性較強(qiáng)，學(xué)界對(duì)此尚未達(dá)成共識(shí)。從研究方法上看，多數(shù)學(xué)者是利用數(shù)據(jù)包絡(luò)分析(data envelopment analysis，DEA)的基本原理，納入非期望產(chǎn)出，運(yùn)用各種擴(kuò)展的模型進(jìn)行測(cè)算，部分學(xué)者也對(duì)傳統(tǒng)模型進(jìn)行了改進(jìn)。從研究對(duì)象上看，既有跨國(guó)研究，也有針對(duì)某一國(guó)的研究；國(guó)內(nèi)學(xué)者對(duì)中國(guó)及省級(jí)層面的研究較多，也有對(duì)特定區(qū)域進(jìn)行研究。遺憾的是，雖然農(nóng)業(yè)在黃河流域占據(jù)重要地位，但關(guān)注黃河流域農(nóng)業(yè)綠色TFP問題的文獻(xiàn)較少。另外，可能由于測(cè)算指標(biāo)、研究方法等不同，不同學(xué)者的研究結(jié)論也相差甚遠(yuǎn)。如展進(jìn)濤等[20]研究表明，中國(guó)農(nóng)業(yè)綠色TFP在2000—2015年間年均下降0.14%；而侯孟陽(yáng)等[21]的測(cè)算結(jié)果顯示，2000年之后中國(guó)農(nóng)業(yè)綠色TFP穩(wěn)定上升。在農(nóng)業(yè)碳排放和面源污染等指標(biāo)的核算上，不同學(xué)者也采用了不同方法[22-23]，事實(shí)上兩者的核算本身也是研究的難點(diǎn)。綜合來看，雖然國(guó)內(nèi)外對(duì)農(nóng)業(yè)綠色TFP有了一定研究，但仍有一些不足有待完善。

基于以上分析，筆者以黃河流域省份為研究對(duì)象，采用考慮非期望產(chǎn)出的SBM模型及多種空間方法，對(duì)農(nóng)業(yè)綠色TFP進(jìn)行測(cè)算和分析。該研究的創(chuàng)新點(diǎn)有：測(cè)算指標(biāo)方面，在總結(jié)借鑒現(xiàn)有研究的基礎(chǔ)上，構(gòu)建了更為完善科學(xué)的農(nóng)業(yè)綠色TFP測(cè)算指標(biāo)體系，力圖測(cè)算結(jié)果準(zhǔn)確、科學(xué)；在指標(biāo)量化上，選擇燈光數(shù)據(jù)衡量農(nóng)業(yè)產(chǎn)值，借鑒聯(lián)合國(guó)政府間氣候變化專門委員會(huì)(IPCC)提供的研究方法核算農(nóng)業(yè)部門碳排放數(shù)據(jù)，較傳統(tǒng)數(shù)據(jù)更客觀；在研究對(duì)象上，黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展是國(guó)家戰(zhàn)略，研究黃河流域農(nóng)業(yè)綠色TFP對(duì)促進(jìn)黃河流域農(nóng)業(yè)綠色發(fā)展和生態(tài)保護(hù)都具有重要參考價(jià)值。

1 研究方法、指標(biāo)體系與數(shù)據(jù)來源

1.1 考慮非期望產(chǎn)出的SBM模型

DEA是一種基于數(shù)學(xué)規(guī)劃的運(yùn)籌學(xué)方法，在效率測(cè)算中得到了國(guó)內(nèi)外學(xué)者的廣泛應(yīng)用。經(jīng)典的DEA模型如CCR[24]、BCC[25]，都是徑向和分段的，容易對(duì)效率值出現(xiàn)高估現(xiàn)象。TONE[26]提出了一種改進(jìn)模型，即非徑向、非角度(non-radial and non-oriented)的基于松弛(slacks-based measure，SBM)的效率測(cè)算方法，有效解決了上述問題。為了解決非期望產(chǎn)出的問題，TONE[27]進(jìn)一步提出了包括非期望產(chǎn)出(undesirable outputs)的SBM模型。根據(jù)研究實(shí)際，筆者采用考慮非期望產(chǎn)出的SBM模型對(duì)農(nóng)業(yè)綠色TFP進(jìn)行測(cè)算，該模型為

(1)

(2)

在實(shí)際測(cè)算中，還應(yīng)考慮規(guī)模報(bào)酬問題。設(shè)U和L分別為規(guī)模報(bào)酬eλ的上下限，則L≤eλ≤U，e=(1,…,1)。當(dāng)兩者都為1時(shí)，表明模型是可變規(guī)模報(bào)酬；當(dāng)L=0、U=1時(shí)，模型是規(guī)模報(bào)酬遞減的；當(dāng)L=1、U=∞時(shí)，模型是規(guī)模報(bào)酬遞增的。

農(nóng)業(yè)生產(chǎn)的規(guī)模報(bào)酬回報(bào)問題在學(xué)術(shù)界尚未達(dá)成一致。傳統(tǒng)經(jīng)濟(jì)理論認(rèn)為，農(nóng)業(yè)生產(chǎn)的規(guī)模報(bào)酬是不變的[28]，但仇童偉等[29]則認(rèn)為傳統(tǒng)農(nóng)業(yè)卷入分工經(jīng)濟(jì)后，顯著提高了農(nóng)業(yè)規(guī)模報(bào)酬。王嫚嫚等[30]卻發(fā)現(xiàn)，部分研究顯示糧食生產(chǎn)中存在規(guī)模報(bào)酬遞減的情況。因此，選擇更為一般的規(guī)模報(bào)酬(general returns to scale，GRS)開展分析。在GRS情形下，可以控制規(guī)模報(bào)酬回報(bào)率在某個(gè)區(qū)間范圍內(nèi)，更符合農(nóng)業(yè)生產(chǎn)的規(guī)律。筆者參照COOPER等[31]的建議，將范圍設(shè)定為L(zhǎng)=0.8、U=1.2。

對(duì)于期望產(chǎn)出和非期望產(chǎn)出，兩者的比例也是模型應(yīng)該考慮的問題。一般來說，非期望產(chǎn)出是伴隨著期望產(chǎn)出產(chǎn)生的，通常假定為兩者同比例增減。COOPER等[31]的實(shí)證研究表明，比例不同只會(huì)造成在前沿面上的投影變化，而最終的結(jié)果排序不受影響。王寶義等[13]的研究表明，不同比例下計(jì)算結(jié)果的大小和趨勢(shì)基本相同。因此該研究確定兩者的比例為1∶1。

綜上，該研究選擇基于一般規(guī)模報(bào)酬的考慮非期望產(chǎn)出的SBM模型，其中期望產(chǎn)出與非期望產(chǎn)出比例為1∶1。

1.2 Kernel 密度估計(jì)

Kernel密度估計(jì)是一種非參數(shù)的密度制圖方法，通過連續(xù)的密度曲線刻畫隨機(jī)變量的分布形態(tài)。隨機(jī)變量x的Kernel 密度估計(jì)方程為

(3)

(4)

帶寬決定了核密度函數(shù)的平滑程度，帶寬越大，方差越小，估計(jì)的偏差越大。因此需要選擇合適的帶寬，該研究采用下式?jīng)Q定帶寬。

(5)

(6)

式(5)～(6)中，σ為方差；span為四分位的跨度。

1.3 空間自相關(guān)分析

探索性空間分析是一種針對(duì)空間數(shù)據(jù)的數(shù)量分析方法，主要用來描述空間分布、識(shí)別空間離群點(diǎn)、判斷空間管理模式等，空間自相關(guān)是其中最常用的分析方法?？臻g自相關(guān)通常分為全域空間自相關(guān)和局域空間自相關(guān)，前者描述某個(gè)變量在整個(gè)研究區(qū)域的空間分布模式，反映的是空間依賴性，后者描述某個(gè)變量在整個(gè)研究區(qū)域內(nèi)與鄰近單元的相關(guān)程度，反映的是空間異質(zhì)性。筆者主要研究黃河流域農(nóng)業(yè)綠色TFP的空間分布特征，因此使用全域自相關(guān)進(jìn)行分析。Moran′sI是應(yīng)用最廣泛的空間自相關(guān)統(tǒng)計(jì)量，公式為

(7)

(8)

1.4 空間馬爾科夫鏈分析

馬爾科夫鏈(Markov chain)是基于隨機(jī)過程理論，通過構(gòu)造狀態(tài)轉(zhuǎn)移概率矩陣，分析事件發(fā)展變化特征。農(nóng)業(yè)綠色TFP從t期到t+1期狀態(tài)的轉(zhuǎn)移概率為

Pij=nij/ni，

(9)

(10)

式(9)～(10)中，P為轉(zhuǎn)移概率；X為農(nóng)業(yè)綠色TFP;nij為在該時(shí)間區(qū)間內(nèi)從狀態(tài)i變?yōu)闋顟B(tài)j的次數(shù)；ni為初始狀態(tài)j出現(xiàn)的總次數(shù)。

將農(nóng)業(yè)綠色TFP分為高、中、低3類，即3種狀態(tài)，由此可以構(gòu)造3×3狀態(tài)轉(zhuǎn)移概率矩陣。將“空間滯后”的概念引入傳統(tǒng)馬爾科夫鏈，構(gòu)造空間馬爾科夫鏈，以此考察空間因素對(duì)狀態(tài)轉(zhuǎn)移概率的影響。某地區(qū)的空間滯后類型取決于農(nóng)業(yè)綠色TFP初始空間滯后算子?？臻g滯后算子是指該地區(qū)周圍鄰居觀察值的空間加權(quán)。對(duì)于地區(qū)i，空間滯后算子(L)的計(jì)算公式為

(11)

式(11)中，xi為觀測(cè)值；wij為空間權(quán)重矩陣中的對(duì)應(yīng)元素。對(duì)于空間權(quán)重矩陣，此處采用地理鄰接0-1矩陣。

將空間滯后類型分為高、中、低3類，將傳統(tǒng)馬爾科夫鏈3×3狀態(tài)轉(zhuǎn)移概率矩陣分解為3個(gè)3×3狀態(tài)轉(zhuǎn)移概率矩陣。通過比較傳統(tǒng)和空間馬爾科夫鏈對(duì)應(yīng)概率，可以分析鄰接地區(qū)對(duì)當(dāng)?shù)剞r(nóng)業(yè)綠色TFP狀態(tài)轉(zhuǎn)移的影響。

1.5 指標(biāo)體系與數(shù)據(jù)來源

以2004—2017年黃河流域各省份為研究范圍，黃河上游包括青海、四川、甘肅、寧夏和內(nèi)蒙古，中游包括陜西和山西，下游包括河南和山東。筆者的研究對(duì)象是狹義農(nóng)業(yè)，即種植業(yè)。該研究在借鑒并總結(jié)已有研究的基礎(chǔ)上，構(gòu)建了一個(gè)更為完善的指標(biāo)體系，并在部分指標(biāo)數(shù)據(jù)的衡量上有所創(chuàng)新(表2)。

表2 農(nóng)業(yè)綠色TFP測(cè)算指標(biāo)體系

投入指標(biāo)包括勞動(dòng)、土地等，涵蓋了現(xiàn)有研究中不同指標(biāo)體系的投入指標(biāo)，同時(shí)還包含了多數(shù)研究忽略的能源和資本2類投入要素。在指標(biāo)數(shù)據(jù)上，個(gè)別指標(biāo)因沒有細(xì)化到種植業(yè)層次，采用相應(yīng)的總量乘以種植業(yè)在生產(chǎn)總值中的占比來衡量。期望產(chǎn)出是種植業(yè)總產(chǎn)值，數(shù)據(jù)來自國(guó)家統(tǒng)計(jì)局網(wǎng)站。非期望產(chǎn)出有農(nóng)業(yè)碳排放和農(nóng)業(yè)面源污染2類。對(duì)于農(nóng)業(yè)碳排放，現(xiàn)有研究多根據(jù)化肥、農(nóng)藥等使用量進(jìn)行推算，容易造成對(duì)碳排放的低估。該研究借鑒SHAN等[32]的方法，采用IPCC分部門排放核算方法計(jì)算農(nóng)業(yè)部門的碳排放數(shù)據(jù)，公式為

(12)

式(12)中，C為碳排放量，百萬t；s為部門；v為能源種類；E為能源消費(fèi)量，需折算成萬t標(biāo)準(zhǔn)煤；VNC為能源低位發(fā)熱量，億J·t-1；CCE為碳排放系數(shù)，t·萬億J-1；FCO為碳氧化因子，一般取100%。能源消費(fèi)量數(shù)據(jù)來源于《中國(guó)能源統(tǒng)計(jì)年鑒》，以各省的統(tǒng)計(jì)年鑒作為補(bǔ)充。其他各變量參考IPCC提供的系數(shù)和GB/T 2589—2020《綜合能耗計(jì)算通則》。

由于種植業(yè)面源污染主要來源于化肥，因此面源污染主要包括總氮和總磷[33]?？偟暮怂惴绞綖檗r(nóng)業(yè)氮肥和復(fù)合肥中含氮量之和乘以氮肥流失率；總磷的核算方式為農(nóng)業(yè)磷肥和復(fù)合肥中含磷量之和乘以磷肥流失率；復(fù)合肥中的含氮(磷)量參考陳同斌等[34]的研究，均取值15%。氮肥流失率和磷肥流失率則因地區(qū)差異而有所不同，具體數(shù)值參考賴斯蕓[35]的研究。

2 結(jié)果分析

2.1 時(shí)間序列分析

從圖1可以看出，2004—2017年黃河流域整體農(nóng)業(yè)綠色TFP先下降后回升，但仍未恢復(fù)至最初水平，因此整體農(nóng)業(yè)綠色TFP略有下降。上、中、下游地區(qū)農(nóng)業(yè)綠色TFP具有分化現(xiàn)象，下游地區(qū)先下降后回升，波動(dòng)幅度較大且同樣未回升到初始水平。上游地區(qū)整體上處于緩慢下降的態(tài)勢(shì)。中游地區(qū)在2005年下降較明顯，之后一直小幅波動(dòng)至2013年，在2014年出現(xiàn)較大的提升，在2017年略微超過了2004年。中游地區(qū)在2005和2014年分別出現(xiàn)了跳躍式下降和回升，這區(qū)別于其他游段，進(jìn)一步分析發(fā)現(xiàn)主要是由于山西省出現(xiàn)了明顯變化，且中游樣本少，導(dǎo)致中游整體波動(dòng)較大。對(duì)比3個(gè)游段的農(nóng)業(yè)綠色TFP可知，下游的效率值一直最高，2010年前上游的農(nóng)業(yè)綠色TFP水平高于中游，但在2014年后中游實(shí)現(xiàn)了反超。不同游段之間的差距也在變大。

圖1 黃河流域及其上、中、下游農(nóng)業(yè)綠色TFP變化趨勢(shì)Fig.1 The changing trend of agricultural green TFP in Yellow River Basin

在研究時(shí)段內(nèi)，農(nóng)業(yè)綠色TFP的變化大體可分為3個(gè)階段。第1階段是2004—2008年，農(nóng)業(yè)綠色TFP整體處于下滑態(tài)勢(shì)，中游和上游地區(qū)下降比較明顯，下游地區(qū)則先大幅下降繼而開始回升，在2008年又出現(xiàn)較大降幅；第2階段是2009—2013年，這段時(shí)間整體農(nóng)業(yè)綠色TFP比較平穩(wěn)，呈小幅上下波動(dòng)，下游地區(qū)則在波動(dòng)中下降且降幅比較明顯。第3階段是2014—2017年，這段時(shí)間整體農(nóng)業(yè)綠色TFP出現(xiàn)了穩(wěn)步回升，中游上升最明顯，其他游段也出現(xiàn)了明顯提高，但仍低于2004年水平。

具體到每個(gè)省份來看，多數(shù)測(cè)算結(jié)果都保持比較穩(wěn)定的態(tài)勢(shì)，但也出現(xiàn)了個(gè)別極值，少數(shù)省份效率值差別過大，與現(xiàn)有研究存在差異。出現(xiàn)這種現(xiàn)象的原因可能是：第一，該研究的測(cè)算指標(biāo)體系中加入了非期望產(chǎn)出農(nóng)業(yè)面源污染，與現(xiàn)有的一些研究只包括一種非期望產(chǎn)出不同，筆者通過測(cè)算不包含農(nóng)業(yè)面源污染的農(nóng)業(yè)綠色TFP發(fā)現(xiàn)，并未出現(xiàn)上述問題；第二，傳統(tǒng)的單階段DEA方法受到環(huán)境因素和隨機(jī)噪聲的影響；第三，DEA對(duì)異常值比較敏感，測(cè)算結(jié)果還與研究對(duì)象、時(shí)間范圍、研究尺度等多方面因素有關(guān)。

從圖2可以明顯看出，上游省份中的四川、甘肅和內(nèi)蒙古農(nóng)業(yè)綠色TFP比較低，而同處上游的寧夏和青海農(nóng)業(yè)綠色TFP較高，且寧夏在黃河流域省份中農(nóng)業(yè)綠色TFP最高，可見上游地區(qū)農(nóng)業(yè)綠色TFP內(nèi)部差異較大。從中游往下，農(nóng)業(yè)綠色TFP逐漸升高。從圖中的分布可知，黃河流域農(nóng)業(yè)綠色TFP分布極不均勻，存在較大的省際差異。

圖2 黃河流域各省份農(nóng)業(yè)綠色TFP對(duì)比Fig.2 Comparison of agricultural green TFP of different provinces in Yellow River Basin

黃河流域各省份農(nóng)業(yè)結(jié)構(gòu)差異較大，西部省份農(nóng)業(yè)中畜牧業(yè)占比高，而東部則以種植業(yè)為主，因此進(jìn)一步將黃河流域分為牧區(qū)省份和農(nóng)區(qū)省份，分析兩者在綠色TFP方面的差異。當(dāng)該省份畜牧業(yè)占農(nóng)業(yè)比值超過30%時(shí)，則認(rèn)為其為牧區(qū)省份。牧區(qū)省份包括青海、四川、寧夏、內(nèi)蒙古和山西，其余為農(nóng)區(qū)省份?？梢钥闯觯羺^(qū)省份均位于中西部。由圖3可見，2類地區(qū)農(nóng)業(yè)綠色TFP在研究時(shí)段內(nèi)表現(xiàn)出先下降后回升的變化趨勢(shì)；農(nóng)區(qū)的農(nóng)業(yè)綠色TFP一直高于牧區(qū)，但從2013年開始兩者的差距明顯縮小。進(jìn)一步分析發(fā)現(xiàn)，牧區(qū)省份內(nèi)部的農(nóng)業(yè)綠色TFP差距明顯大于農(nóng)區(qū)省份，圖4展示了農(nóng)區(qū)和牧區(qū)農(nóng)業(yè)綠色TFP的方差和極差(最大值與最小值的差值)變化，可以看出牧區(qū)的內(nèi)部分化和差異明顯大于農(nóng)區(qū)。

圖3 黃河流域農(nóng)區(qū)與牧區(qū)農(nóng)業(yè)綠色TFP變化趨勢(shì)

圖4 黃河流域農(nóng)區(qū)與牧區(qū)農(nóng)業(yè)綠色TFP的方差與極差變化趨勢(shì)Fig.4 The variance and range trend of agricultural green TFP of farming and animal husbandry areas in Yellow River Basin

2.2 空間自相關(guān)分析

為了探索黃河流域省份農(nóng)業(yè)綠色TFP的空間相關(guān)性，基于地理鄰接0-1矩陣和地理距離矩陣計(jì)算全局Moran′sI值，結(jié)果見表3。從結(jié)果來看，Moran′sI值都沒有通過顯著性檢驗(yàn)，在多數(shù)年份下為負(fù)，若干年份在鄰接矩陣下為正。從Moran′sI值的變化趨勢(shì)來看，2種矩陣下指數(shù)的變化一致。以上結(jié)果表明，黃河流域農(nóng)業(yè)綠色TFP以分散態(tài)勢(shì)為主，空間效應(yīng)較弱，整體上沒有形成一個(gè)有機(jī)整體，存在“各自為政”的現(xiàn)象；整體的變化趨勢(shì)也不穩(wěn)定，說明黃河流域農(nóng)業(yè)綠色TFP在空間協(xié)作上沒有連貫性，缺乏整體性和一體化思維。與長(zhǎng)江經(jīng)濟(jì)帶相比，黃河流域的整體通航能力較弱，同時(shí)由于黃河河道呈“幾”字型，客觀上導(dǎo)致黃河流域各省份之間的空間關(guān)聯(lián)較弱。黃河流域環(huán)境治理是一個(gè)系統(tǒng)工程，需要流域內(nèi)各省份協(xié)同推進(jìn)。黃河流域農(nóng)業(yè)綠色TFP在空間上的分散對(duì)于整個(gè)流域的綠色發(fā)展起到了阻礙作用。

表3 Moran′s I值測(cè)算結(jié)果

2.3 時(shí)間演變分析

采用Kernel 密度估計(jì)黃河流域農(nóng)業(yè)綠色TFP的動(dòng)態(tài)演進(jìn)情況，部分代表年份的變化見圖5。從分布位置上看，2004、2016和2017年呈明顯的雙峰分布，這種分布形態(tài)下農(nóng)業(yè)綠色TFP呈現(xiàn)兩極分化態(tài)勢(shì)。其中2004年主峰位于右側(cè)，即農(nóng)業(yè)綠色TFP較高的省份居多；2016年主峰位于左側(cè)，主峰高度增高，同時(shí)次鋒向左偏移，峰的寬帶變?。?017年則沒有明顯的主峰，2個(gè)波峰高度和寬帶比較接近。2005—2015年分布形態(tài)穩(wěn)定，為單峰分布，主峰位于左側(cè)。這段時(shí)間出現(xiàn)了比較明顯的右拖尾現(xiàn)象，表明這期間農(nóng)業(yè)綠色TFP較低，但有個(gè)別省份農(nóng)業(yè)綠色TFP水平較高，從而拖長(zhǎng)了曲線分布，也表明農(nóng)業(yè)綠色TFP存在極化現(xiàn)象。進(jìn)一步分析發(fā)現(xiàn)，在2005—2015年間，主峰高度先上升后下降，而寬帶表現(xiàn)為逐漸擴(kuò)大。雙峰分布意味著高水平TFP和低水平TFP集聚現(xiàn)象并存，而單峰則表現(xiàn)為低水平集聚。從2004—2017年的整體變化可以看出，黃河流域農(nóng)業(yè)綠色TFP經(jīng)歷了雙峰—單峰—雙峰的變化，極化趨勢(shì)呈現(xiàn)出兩極化—單極化—兩極化的變化，而整個(gè)階段的地區(qū)差異一直存在。

為了方便展示，僅列出部分代表年份。

2.4 馬爾科夫鏈分析

通過求解馬爾科夫鏈轉(zhuǎn)移矩陣考察黃河流域農(nóng)業(yè)綠色TFP的時(shí)空變化情況。根據(jù)農(nóng)業(yè)綠色TFP的整體分布情況，將黃河流域農(nóng)業(yè)綠色TFP劃分為低、中、高3個(gè)類型，2個(gè)分界點(diǎn)分別是0.25和0.75，3個(gè)類型分別用1、2、3進(jìn)行表示。將研究區(qū)間劃分為2004—2010年和2011—2017年2個(gè)時(shí)間段進(jìn)行考察。傳統(tǒng)馬爾科夫轉(zhuǎn)移矩陣如表4所示，矩陣中的數(shù)字為發(fā)生轉(zhuǎn)換的概率。如在2004—2010年間的矩陣中，第1列第1行即表示在t期處于低水平的省份有87.5%的轉(zhuǎn)換沒有發(fā)生類型變化，而有12.5%的概率是從低水平向上轉(zhuǎn)移為中等水平，沒有低水平狀態(tài)直接上升為高水平狀態(tài)。據(jù)此分析可知，2004—2010年間，農(nóng)業(yè)綠色TFP處在中等水平的省份，有12.5%的概率向下轉(zhuǎn)移至低水平狀態(tài)，18.75%的概率向上轉(zhuǎn)移至高水平狀態(tài)，而更多的是同狀態(tài)轉(zhuǎn)移(68.75%)，即仍保留在中等水平。農(nóng)業(yè)綠色TFP處在高水平的省份，有71.43%的概率仍保持在高水平，有28.57%的概率會(huì)向下轉(zhuǎn)移為中等水平。2011—2017年間，概率轉(zhuǎn)移矩陣發(fā)生了變化。低水平的省份向上轉(zhuǎn)移為中等水平的概率下降(4.76%)，更多的會(huì)保留在低水平(95.23%)，仍沒有省份直接從低水平跨越至高水平。對(duì)于中等水平省份來說，在該區(qū)間沒有發(fā)生向下轉(zhuǎn)移，向上轉(zhuǎn)移至高水平的概率下降至14.29%，但保留在中等水平的概率提高至85.71%。高水平省份不再發(fā)生向下轉(zhuǎn)移，而是全部保留在了高水平狀態(tài)。對(duì)角線上的概率大于非對(duì)角線上的概率，這表明存在“俱樂部收斂”，即同一狀態(tài)水平的省份具有趨同態(tài)勢(shì)。綜合來看，低水平狀態(tài)向上轉(zhuǎn)移的概率在下降，且較難直接跨越至高水平；而高水平狀態(tài)向下轉(zhuǎn)移的概率變小，且一般不會(huì)降為低水平；中等水平省份會(huì)發(fā)生向上或向下轉(zhuǎn)移，但向下的概率不斷降低。

表4 傳統(tǒng)馬爾科夫鏈概率轉(zhuǎn)移矩陣

為了克服傳統(tǒng)馬爾科夫鏈轉(zhuǎn)移矩陣缺乏空間因素的缺陷，在傳統(tǒng)馬爾科夫鏈概率轉(zhuǎn)移矩陣的基礎(chǔ)上引入空間權(quán)重矩陣。以各地初始年份的空間滯后算子將黃河流域省份劃分為鄰近低、中、高水平3類，分別用Ⅰ、Ⅱ、Ⅲ表示。由表5可知，第一，空間格局對(duì)于黃河流域農(nóng)業(yè)綠色TFP具有重要影響。通過對(duì)比空間與傳統(tǒng)2類馬爾科夫鏈概率轉(zhuǎn)移矩陣，發(fā)現(xiàn)概率發(fā)生了變化，且不同空間類型的變化也不盡相同。第二，在鄰接類型Ⅰ下，即鄰接地區(qū)為低水平，2004—2010年間低水平和中等水平的類型均保持不變，而高水平向下轉(zhuǎn)移為中等水平的概率為40%，保持高水平的概率為60%，與傳統(tǒng)馬爾科夫鏈概率轉(zhuǎn)移矩陣相比，向上轉(zhuǎn)移和保持高水平的概率都明顯下降。與此相反的是，當(dāng)鄰接類型為Ⅲ時(shí)，低水平能夠向上躍升為高水平(2.5%)，而中等水平向上遷移和高水平向下遷移的概率都出現(xiàn)下降，2011—2017年的結(jié)果也表現(xiàn)出以上特征。第三，空間馬爾科夫鏈概率轉(zhuǎn)移表現(xiàn)出的“效率鎖定”現(xiàn)象為“俱樂部收斂”提供了空間解釋?！靶舒i定”現(xiàn)象指鄰接地區(qū)為低水平時(shí)，向上轉(zhuǎn)移的概率下降；而鄰接地區(qū)為高水平時(shí)，向下轉(zhuǎn)移的概率也下降。黃河流域農(nóng)業(yè)綠色TFP表現(xiàn)出高低各自集聚的“俱樂部收斂”現(xiàn)象。

表5 空間馬爾科夫鏈概率轉(zhuǎn)移矩陣

3 結(jié)論

以2004—2017年黃河流域省份為研究對(duì)象測(cè)算黃河流域農(nóng)業(yè)綠色TFP，并對(duì)其空間格局和動(dòng)態(tài)演進(jìn)進(jìn)行研究。在總結(jié)借鑒已有研究的基礎(chǔ)上，提出了一個(gè)更為完善的農(nóng)業(yè)綠色TFP測(cè)算指標(biāo)體系，運(yùn)用考慮非期望產(chǎn)出的SBM模型，使用空間自相關(guān)分析、Kernel 密度估計(jì)和馬爾科夫鏈概率分析等方法，分析了黃河流域農(nóng)業(yè)綠色TFP的空間格局和動(dòng)態(tài)演進(jìn)規(guī)律，主要得出以下結(jié)論：第一，黃河流域農(nóng)業(yè)綠色TFP整體上并不高，有很大的提升空間，且上、中、下游具有明顯差異。2004—2017年變化趨勢(shì)整體上為略微下降，在此過程中上下波動(dòng)起伏，各省份間差異較大。第二，Kernel 密度估計(jì)結(jié)果表明，黃河流域農(nóng)業(yè)綠色TFP經(jīng)歷了雙峰—單峰—雙峰的變化趨勢(shì)，極化趨勢(shì)呈現(xiàn)兩極化—單極化—兩極化的變化。從這種變化來看，黃河流域的農(nóng)業(yè)綠色TFP并不穩(wěn)定，沒有形成比較明顯和穩(wěn)定的發(fā)展態(tài)勢(shì)。第三，空間自相關(guān)分析發(fā)現(xiàn)，黃河流域農(nóng)業(yè)綠色TFP處于空間分化狀態(tài)，整體的空間關(guān)聯(lián)性較弱，沒有形成統(tǒng)一的有機(jī)整體。結(jié)合馬爾科夫鏈分析可知，高水平和低水平農(nóng)業(yè)綠色TFP地區(qū)各自集聚，從而使空間分化更加明顯，形成“俱樂部收斂”。第四，由傳統(tǒng)和空間馬爾科夫鏈分析可知，低水平地區(qū)較難躍升為高水平地區(qū)，而保持不變的概率較大。在空間滯后的作用下，“效率鎖定”現(xiàn)象明顯，即鄰接地區(qū)為低水平時(shí)，向上轉(zhuǎn)移的概率下降；而鄰接地區(qū)為高水平時(shí)更容易保持在高水平。

從以上分析可以看出，黃河流域農(nóng)業(yè)綠色TFP具有以下3個(gè)特點(diǎn)：第一，整體效率偏低。黃河流域在整體農(nóng)業(yè)綠色TFP約為0.433，屬于較低水平。從投入-產(chǎn)出體系分析，主要原因在于農(nóng)業(yè)投入多，而農(nóng)業(yè)產(chǎn)出多的同時(shí)非期望產(chǎn)出也在增大。第二，地區(qū)差異明顯。黃河流域跨東、中、西三大板塊，流經(jīng)我國(guó)三大階梯。各省份具有不同的地域特色，經(jīng)濟(jì)發(fā)展水平、要素稟賦、環(huán)境氣候等不盡相同，各地農(nóng)業(yè)種植結(jié)構(gòu)、技術(shù)水平也存在差異，從而導(dǎo)致了農(nóng)業(yè)綠色TFP的地區(qū)差異。地區(qū)差異明顯造成“一條腿長(zhǎng)、一條腿短”的失衡狀態(tài)，從而抑制了整體效率的提升。第三，空間關(guān)聯(lián)偏弱。黃河流域整體的通航能力較弱，雖然各省份都依黃河而興農(nóng)業(yè)，但彼此的關(guān)聯(lián)卻不緊密，不利于形成合理的農(nóng)業(yè)空間布局。

4 政策啟示

第一，著力提升農(nóng)業(yè)綠色TFP。從投入-產(chǎn)出體系角度看，加大使用農(nóng)藥、化肥、農(nóng)膜等重要的農(nóng)業(yè)生產(chǎn)資料，會(huì)同時(shí)帶來農(nóng)業(yè)污染。因此，應(yīng)加強(qiáng)落實(shí)農(nóng)藥、化肥零增長(zhǎng)行動(dòng)方案。改進(jìn)施肥方式，提高肥料利用率，減少不合理投入。大力推廣新型農(nóng)藥，提升裝備水平，加快轉(zhuǎn)變病蟲害防控方式，實(shí)現(xiàn)農(nóng)藥減量控害，保障農(nóng)業(yè)生產(chǎn)安全、農(nóng)產(chǎn)品質(zhì)量安全和生態(tài)環(huán)境安全。同時(shí)應(yīng)加大農(nóng)業(yè)科技投入，特別是結(jié)合種業(yè)振興方案，提高農(nóng)作物育種質(zhì)量，提高農(nóng)業(yè)全要素生產(chǎn)率。

第二，應(yīng)樹立黃河流域整體性思維，協(xié)同推進(jìn)農(nóng)業(yè)綠色發(fā)展。一方面，黃河流域農(nóng)業(yè)綠色TFP本身具有較大提升空間，省際差異較大，空間分化現(xiàn)象嚴(yán)重；另一方面，空間分布格局對(duì)農(nóng)業(yè)綠色TFP的動(dòng)態(tài)演化具有重要影響。黃河流域各省份間的關(guān)聯(lián)性和協(xié)同性較差，應(yīng)運(yùn)用整體性思維，打破行政壁壘，推動(dòng)黃河流域區(qū)域協(xié)調(diào)發(fā)展。

第三，落實(shí)“綠水青山就是金山銀山”的綠色發(fā)展理念。農(nóng)業(yè)在黃河流域省份占比相對(duì)較高，農(nóng)業(yè)是基礎(chǔ)產(chǎn)業(yè)，但同時(shí)也面臨嚴(yán)峻的生態(tài)考驗(yàn)。因此在黃河流域農(nóng)業(yè)發(fā)展上，要應(yīng)用“綠水青山就是金山銀山”的理念系統(tǒng)規(guī)劃農(nóng)業(yè)發(fā)展。特別是要加強(qiáng)農(nóng)業(yè)面源污染治理，加強(qiáng)污染監(jiān)測(cè)監(jiān)管，種養(yǎng)結(jié)合，發(fā)展循環(huán)綠色農(nóng)業(yè)。針對(duì)碳達(dá)峰碳中和的遠(yuǎn)期目標(biāo)，農(nóng)業(yè)領(lǐng)域應(yīng)提前做好謀劃。

第四，因地制宜發(fā)展特色農(nóng)業(yè)，化“空間差異”為合理梯度，避免同質(zhì)化惡性競(jìng)爭(zhēng)。根據(jù)各省份的地理、氣候、要素稟賦等因素，制定適合當(dāng)?shù)貙?shí)際的農(nóng)業(yè)綠色發(fā)展規(guī)劃，避免“一刀切”思維。各地區(qū)應(yīng)根據(jù)當(dāng)?shù)剞r(nóng)業(yè)發(fā)展階段、農(nóng)業(yè)特點(diǎn)制定農(nóng)業(yè)綠色發(fā)展方案。對(duì)于農(nóng)業(yè)綠色TFP較高的省份，可在保持穩(wěn)定增長(zhǎng)的基礎(chǔ)上，發(fā)展高值、特色農(nóng)業(yè)，探索農(nóng)業(yè)多元經(jīng)營(yíng)。對(duì)于農(nóng)業(yè)TFP較低的省份，應(yīng)將重點(diǎn)放在提高農(nóng)業(yè)產(chǎn)出上，加強(qiáng)農(nóng)業(yè)污染治理，促進(jìn)農(nóng)業(yè)可持續(xù)綠色發(fā)展?？梢岳命S河流域的空間差異，因地制宜地進(jìn)行農(nóng)業(yè)空間分工，形成各地區(qū)聯(lián)合發(fā)展的合力。