□文/鄧 清 劉紋靜 高 超 成燕妮 王翼航 時 悅

（黑龍江八一農(nóng)墾大學經(jīng)濟管理學院 黑龍江·大慶）

［提要］ 基于SBM-GML 模型測算黑龍江省13 個地級市的農(nóng)業(yè)綠色全要素生產(chǎn)率，并分析其時空分異特征。通過fsQCA 方法探究農(nóng)業(yè)綠色全要素生產(chǎn)率提升路徑，結果表明：黑龍江省農(nóng)業(yè)綠色全要素生產(chǎn)率發(fā)展呈現(xiàn)分散化的特點；主要有“經(jīng)濟發(fā)展”主導“貿易開放”輔助型、“貿易開放”主導“灌溉效率”輔助型、“灌溉效率水平”主導“市場化程度”輔助型等六條提升路徑。對此，各地政府應著力提高農(nóng)業(yè)科技水平，充分發(fā)揮地區(qū)資源優(yōu)勢，優(yōu)化要素投入結構，推進農(nóng)業(yè)綠色產(chǎn)業(yè)轉型升級，加強區(qū)域合作實現(xiàn)區(qū)域均衡發(fā)展，推動提高農(nóng)業(yè)綠色全要素生產(chǎn)率。

黨的二十大報告明確提出要加快發(fā)展方式綠色轉型，積極穩(wěn)妥推進碳達峰、碳中和。農(nóng)業(yè)在實現(xiàn)經(jīng)濟社會發(fā)展綠色化、低碳化過程中具有不可替代的作用，著力提升農(nóng)業(yè)綠色發(fā)展水平將成為未來一段時期農(nóng)業(yè)發(fā)展的核心任務。黑龍江省作為中國糧食安全“壓艙石”，正面臨水土資源過度消耗、黑土地質量下降風險加大以及農(nóng)業(yè)面源污染等問題，制約著黑龍江省農(nóng)業(yè)綠色轉型。由于農(nóng)業(yè)綠色全要素生產(chǎn)率能夠全方位衡量區(qū)域農(nóng)業(yè)綠色轉型升級的綜合水平，因此研究黑龍江省農(nóng)業(yè)綠色全要素生產(chǎn)率多重影響因素的聯(lián)動作用，對于研究凝練農(nóng)業(yè)綠色全要素生產(chǎn)率的多元化提升路徑具有重要的理論和現(xiàn)實意義。

一、文獻綜述

（一）農(nóng)業(yè)綠色全要素生產(chǎn)率的測算。國外現(xiàn)有研究成果對農(nóng)業(yè)綠色全要素生產(chǎn)率的測算以數(shù)據(jù)包絡分析技術（DEA）為基礎的非參數(shù)方法為主。國內相關研究較多使用的是基于徑向、角度的DEA 模型。傳統(tǒng)的全要素生產(chǎn)率測算方法沒有考慮農(nóng)業(yè)經(jīng)濟增長的環(huán)境代價，其測算結果無法真實反映環(huán)境約束下的農(nóng)業(yè)生產(chǎn)效率，如徑向DEA 模型無法處理投入冗余和產(chǎn)出不足產(chǎn)生的松弛變量問題，角度DEA 模型需要假定投入或產(chǎn)出不變，這與現(xiàn)實并不完全相符。為克服傳統(tǒng)DEA 模型的缺陷，Chung 提出了方向性距離函數(shù)的概念，據(jù)此構建了Malmquist-Luenberger（ML）生產(chǎn)率指數(shù)，有效解決了期望產(chǎn)出和非期望產(chǎn)出聯(lián)合生產(chǎn)的環(huán)境效率評價問題。國內學者們則構建了非徑向、非角度SBM 模型，如李谷成、閆桂權使用SBM-ML 模型分析了我國農(nóng)業(yè)綠色全要素生產(chǎn)率。為解決線性規(guī)劃無可行解的問題，一些學者使用全局Malmquist-Luenberger（GML）指數(shù)測算農(nóng)業(yè)綠色全要素生產(chǎn)率，如馬國群采用SBM-GML 模型分析發(fā)現(xiàn)環(huán)境規(guī)制對農(nóng)業(yè)綠色全要素生產(chǎn)率的影響存在顯著的雙重門檻效應。這些研究成果已經(jīng)成為測算農(nóng)業(yè)綠色全要素生產(chǎn)率的重要依據(jù)。

（二）農(nóng)業(yè)綠色全要素生產(chǎn)率的影響因素。關于農(nóng)業(yè)綠色全要素生產(chǎn)率的影響因素，現(xiàn)有文獻主要從技術因素、制度因素、市場因素和農(nóng)業(yè)生產(chǎn)條件等方面展開研究。技術支撐是農(nóng)業(yè)綠色發(fā)展的基礎，其支撐作用在實證研究中得到了證實。如，Reinhard 認為技術條件既可以直接影響農(nóng)業(yè)環(huán)境效率，也可以通過投入要素間接影響環(huán)境效率。制度因素對農(nóng)業(yè)綠色發(fā)展的影響最為廣泛，幾乎在農(nóng)業(yè)發(fā)展的所有層面上都具有重要作用，并且農(nóng)業(yè)綠色發(fā)展具有一定的公共品屬性，需要一定的制度環(huán)境加以引導和規(guī)范。李谷成從宏觀層面分析了我國農(nóng)業(yè)綠色全要素生產(chǎn)率變動的制度原因。王建華從微觀層面檢驗了環(huán)境管制政策對農(nóng)業(yè)綠色發(fā)展的影響?！丁笆奈濉比珖r(nóng)業(yè)綠色發(fā)展規(guī)劃》明確提出農(nóng)業(yè)綠色發(fā)展需要堅持政府引導和市場主導的基本原則，意味著市場化將成為影響農(nóng)業(yè)綠色發(fā)展的重要因素。相關研究認為農(nóng)業(yè)市場化有助于改善農(nóng)村生態(tài)環(huán)境質量。Plassmann 強調人們必須重視市場需求在減少農(nóng)業(yè)污染方面的重要作用，否則，最復雜、最先進的技術也無法阻止農(nóng)業(yè)污染排放的增加。完善的農(nóng)業(yè)生產(chǎn)條件能夠有效緩解農(nóng)業(yè)生態(tài)環(huán)境壓力，提高農(nóng)業(yè)綠色發(fā)展水平。如，李自強對農(nóng)業(yè)灌溉系統(tǒng)的研究就在一定程度上表明了農(nóng)業(yè)生產(chǎn)條件對農(nóng)業(yè)綠色發(fā)展的重要性。

綜上所述，農(nóng)業(yè)綠色全要素生產(chǎn)率的測算方法日趨成熟，為本研究分析黑龍江省農(nóng)業(yè)綠色全要素生產(chǎn)率現(xiàn)狀提供了理論和方法。學術界從多方面研究了農(nóng)業(yè)綠色全要素生產(chǎn)率的影響因素，為本研究構建理論框架提供了重要依據(jù)。但已有研究主要探究了單個影響因素對農(nóng)業(yè)綠色全要素生產(chǎn)率的凈效應，忽視了影響因素之間的協(xié)同聯(lián)動對農(nóng)業(yè)綠色全要素生產(chǎn)率的影響。多因素聯(lián)動作用的研究有助于識別出農(nóng)業(yè)綠色全要素生產(chǎn)率的多種提升路徑，深化農(nóng)業(yè)綠色全要素生產(chǎn)率提升背后的多因素共同作用的本質認識。

二、數(shù)據(jù)來源與指標選取

（一）數(shù)據(jù)來源。本文選取2012～2021 年黑龍江省13 個地級市的面板數(shù)據(jù)。用于測算黑龍江省各地市農(nóng)業(yè)綠色全要素生產(chǎn)率的相關指標如農(nóng)業(yè)從業(yè)人員、農(nóng)業(yè)總播種面積、農(nóng)業(yè)機械總動力等來源于《黑龍江省統(tǒng)計年鑒》和各地市統(tǒng)計年鑒；前因變量如農(nóng)業(yè)機械化水平、經(jīng)濟發(fā)展水平、灌溉效率水平等指標均來源于EPS 數(shù)據(jù)庫。

（二）指標選取。農(nóng)業(yè)綠色全要素生產(chǎn)率測算，選取的投入產(chǎn)出指標如表1 所示。（表1）

農(nóng)業(yè)碳排放量總量的計算方法，本文參考劉亦文的方法，使用IPCC 碳排放系數(shù)對農(nóng)業(yè)灌溉、翻耕、化肥進行估算，考慮到數(shù)據(jù)精確性與無缺性，以農(nóng)業(yè)灌溉、翻耕面積、化肥施用折純量作為碳排放源。

（三）變量選取

1、變量定義?？紤]到農(nóng)業(yè)的復雜性與數(shù)據(jù)可得性等因素，分別參考黃秀全、熊翅新、王璐、謝賢君、陳燕君等的做法，采用農(nóng)業(yè)機械化水平、經(jīng)濟發(fā)展水平、貿易開發(fā)程度、市場化程度、灌溉效率水平等5 個變量作為農(nóng)業(yè)綠色全要素生產(chǎn)率的前因變量。將農(nóng)業(yè)綠色全要素生產(chǎn)率作為結果變量。主要變量說明見表2。（表2）

2、變量校準。采用分位數(shù)值來確定錨點的直接法對結果與前因條件進行校準。結果與前因條件均采用直接校準法，所有連續(xù)變量的完全隸屬、交叉臨界點和完全不隸屬的三個閾值，在此分別設置為95%、50%和5%三個分位數(shù)值，若對校準數(shù)據(jù)中出現(xiàn)0.5 數(shù)值的條件，將其數(shù)值替換處理為0.501，以避免對0.5 數(shù)值所屬案例進行理論分析的困難。校準結果如表3 所示。（表3）

（四）農(nóng)業(yè)綠色全要素生產(chǎn)率測算方法。本文使用SBM-GML 模型測算黑龍江省農(nóng)業(yè)綠色全要素生產(chǎn)率。

1、超效率SBM 模型。為有效解決農(nóng)業(yè)碳排放在測算過程中對效率的影響，采用Tone 提出的超效率SBM 模型，相關SBM 模型測算方法已得到有效驗證，參考焦兵、王騰、王兵等相關研究。

2、GML 指數(shù)。為在線性規(guī)劃中得到可行解，參考OH 提出的GML指數(shù)研究思路，借鑒劉鉆擴、朱寒冰等的方法，公式如下：

公式中，（xt，gx），表示投入和投入減少方向向量，（yt，gy）表示期望產(chǎn)出和期望產(chǎn)出增加方向向量，（bt，gb）表示非期望產(chǎn)出和非期望產(chǎn)出減少方向向量分別代表基于非徑向、非角度測度的當期和全域SBM 方向距離函數(shù)分別表示從t 到t+1 時期的農(nóng)業(yè)綠色全要素生產(chǎn)率、技術進步、技術效率。

三、結果與分析

（一）測算結果分析。測算并分解黑龍江省13 個地級市2012～2021年農(nóng)業(yè)綠色全要素生產(chǎn)率，結果如表4 所示。農(nóng)業(yè)綠色全要素生產(chǎn)率整體趨于平穩(wěn)，但在個別年份間存在波動，其中2015～2016 年農(nóng)業(yè)綠色全要素生產(chǎn)率為0.948，與上一年相比下降11.81%。盡管技術進步較上一年提高了7.96%，對農(nóng)業(yè)綠色全要素生產(chǎn)率產(chǎn)生正向影響，但是由于純技術效率和規(guī)模效率分別下降4.35%和14.63%，最終導致農(nóng)業(yè)綠色全要素生產(chǎn)率下降。（表4）

由表5 可知，橫向對比來看，黑龍江省內13 個地級市農(nóng)業(yè)綠色全要素生產(chǎn)率存在顯著的區(qū)域性差異。如，綏化市和雞西市等地農(nóng)業(yè)綠色全要素生產(chǎn)率的提升主要依靠強勁的技術效率和技術進步，其中鶴崗市純技術效率高達1.077；然而，齊齊哈爾市和大慶市主要是由規(guī)模效率驅動的農(nóng)業(yè)綠色全要素生產(chǎn)率提升。（表5）

表5 農(nóng)業(yè)綠色全要素生產(chǎn)率的空間差異一覽表

（二）必要性條件分析。為探求單個因素是否為結果的必要性條件，借鑒Ragin 的方法，對前因變量進行必要性檢驗，是否產(chǎn)生以結果集合作為條件集合（或者其非集）的子集，必要性檢驗結果如表6 所示。單個因素的一致性水平均不超過0.9，因此單個因素均不是必要性條件，表明其提升路徑受多因素的組態(tài)效應影響。（表6）

表6 單個因素必要性條件檢驗結果一覽表

（三）條件組態(tài)結果。為得到提升路徑條件組態(tài)，借鑒張明、杜運周等方法，將初始一致性閾值設為0.7，PRI 不低于0.7，頻數(shù)閾值設為1，得到黑龍江13 個地級市農(nóng)業(yè)綠色全要素生產(chǎn)率相應組態(tài)，存在6 條組態(tài)增長路徑。第1 條提升路徑：“經(jīng)濟發(fā)展”主導“貿易開放”輔助型；第2 條提升路徑：“經(jīng)濟發(fā)展”主導“灌溉效率”輔助型；第3 條提升路徑：“貿易開放”主導“灌溉效率”輔助型；第4 條提升路徑：“對外開放”主導“農(nóng)業(yè)機械化+市場化程度”輔助型；第5 條提升路徑：“經(jīng)濟發(fā)展水平”主導“市場化程度”輔助型；第6 條提升路徑：“灌溉效率水平”主導“市場化程度”輔助型。這6 條提升路徑的一致性都高于0.8，高于初設臨界值0.7，總體覆蓋度0.636，可以解釋超63%的案例，說明6 條路徑的解釋力度均較好。（表7）

表7 組態(tài)分析結果一覽表

（四）條件組態(tài)分析

1、經(jīng)濟發(fā)展——貿易開放型?；诖寺窂较拢崔r(nóng)業(yè)機械化程度低、市場資源配置結構未得到優(yōu)化，依靠較高的經(jīng)濟發(fā)展水平和貿易開放促進農(nóng)業(yè)綠色全要素生產(chǎn)率提高。解釋典型地區(qū)為伊春市、鶴崗市。以伊春市為例，經(jīng)濟發(fā)展水平對農(nóng)業(yè)生產(chǎn)技術進步、農(nóng)業(yè)生產(chǎn)結構等方面具有相當?shù)挠绊懥Γ?jīng)濟發(fā)展水平較高，意味著其資源配置效率較高，農(nóng)業(yè)基礎設施更加完善等，能夠有效提高農(nóng)村居民純收入；伊春市嘉蔭口岸和市政府推進貿易開放的措施政策等提高貿易開放程度。2012～2021 年伊春市農(nóng)村居民純收入由0.92 萬元增至1.94 萬元，增長110.87%，農(nóng)民有較高的收入意味著農(nóng)民所掌握的生產(chǎn)技術較高，有著較高的生活質量，在農(nóng)業(yè)生產(chǎn)、提升生產(chǎn)技術水平等方面分配更多的時間精力；貿易開放在需求端刺激農(nóng)業(yè)綠色生產(chǎn)率的提高，更為嚴格的農(nóng)產(chǎn)品出口標準在一定程度上引導農(nóng)業(yè)向著綠色化和標準化發(fā)展，從而較高的農(nóng)業(yè)經(jīng)濟發(fā)展水平和貿易開放程度的提高能夠促進伊春農(nóng)業(yè)綠色全要素生產(chǎn)率的提升。

2、經(jīng)濟發(fā)展——灌溉效率型。此路徑下，農(nóng)業(yè)綠色全要素生產(chǎn)率主要是由高經(jīng)濟發(fā)展水平主導和良好的灌溉效率水平輔助共同推動。該路徑中的典型地區(qū)為大慶市。以大慶市為例，大慶市在2019 年印發(fā)《大慶市市級農(nóng)民專業(yè)合作社示范社評定管理辦法》，著重提出加快構建現(xiàn)代農(nóng)業(yè)體系，促進農(nóng)業(yè)經(jīng)濟健康規(guī)范發(fā)展，提高農(nóng)民純收入，進一步提高其農(nóng)業(yè)經(jīng)濟發(fā)展水平，其農(nóng)民純收入由2012 年的1.15 萬元增至2021 年2.04 萬元，同比增長77.39%；灌溉效率水平10 年間同比增幅87.05%，排名升至黑龍江省第一，有效促進農(nóng)業(yè)綠色化發(fā)展。因此，農(nóng)業(yè)經(jīng)濟發(fā)展水平和灌溉效率水平的提高促進大慶市農(nóng)業(yè)綠色全要素生產(chǎn)率的提升。

3、貿易開放——灌溉效率型。此路徑下，在農(nóng)業(yè)機械化水平較為落后和市場資源配置結構未得到優(yōu)化的情況下，高貿易開放程度和灌溉效率水平的提高能夠有效促進農(nóng)業(yè)綠色全要素生產(chǎn)率的提高。該路徑中的典型地區(qū)為齊齊哈爾市。以齊齊哈爾市為例，在外貿方面，打造國際貿易供給力、營造對外貿易健康發(fā)展環(huán)境等全面擴大貿易開放程度；實施高效節(jié)水灌溉、節(jié)水增糧、灌溉基礎設施建設等一系列措施有效提高灌溉效率水平。貿易開放主導、灌溉效率輔助構成農(nóng)業(yè)綠色全要素生產(chǎn)率提升路徑。

4、貿易開放——農(nóng)業(yè)機械化+市場化程度型。此路徑下，即使經(jīng)濟發(fā)展水平相對較低，但是在高農(nóng)業(yè)機械化水平和高市場化程度因素的推動下，依靠高貿易開放程度仍能促進農(nóng)業(yè)綠色全要素生產(chǎn)率提高。該路徑下，典型地區(qū)為黑河市。以黑河市為例，2019 年國務院印發(fā)《中國（黑龍江）自由貿易試驗區(qū)總體方案》，明確提出黑龍江自貿試驗區(qū)各片區(qū)產(chǎn)業(yè)發(fā)展的功能布局和重點領域，自貿區(qū)黑河片區(qū)投入建設以來，以多元聯(lián)動發(fā)展為目標，形成數(shù)字自貿區(qū)等創(chuàng)新發(fā)展理念和“一脈串六區(qū)”空間新格局，成功全面提高了黑河市對外開放程度。黑河市統(tǒng)計年鑒數(shù)據(jù)顯示，2012～2021 年農(nóng)業(yè)機械化水平由8.38%增長至22.9%，農(nóng)業(yè)機械化水平不斷提高；市場化程度由31.46%增長至40.88%，市場化程度不斷提升。因此，在農(nóng)業(yè)機械化水平和市場化程度因素推動下，貿易開放程度的提高能促進黑河市農(nóng)業(yè)綠色全要素生產(chǎn)率的提升。

5、經(jīng)濟發(fā)展水平——市場化程度型。此路徑下，即使農(nóng)業(yè)機械化水平程度、貿易開放程度和灌溉效率水平低，但由于高經(jīng)濟發(fā)展水平和高市場化程度因素的助推，農(nóng)業(yè)綠色全要素生產(chǎn)率仍能得到提升。該形成路徑下，典型地區(qū)為佳木斯市。佳木斯市在《佳木斯市促進平臺經(jīng)濟高質量發(fā)展實施意見》中強調要抓住新經(jīng)濟發(fā)展窗口機遇，著力打造具有地方特色的平臺經(jīng)濟以帶動區(qū)塊鏈擴展和升級。佳木斯市統(tǒng)計年鑒數(shù)據(jù)顯示，2012～2021 年其農(nóng)村居民人均純收入由9，823 元增長至2.12 萬元，經(jīng)濟發(fā)展水平顯著提高；市場化程度由24.04%增長至41.24%，上升了17.2%，市場化程度不斷擴大。因此，經(jīng)濟發(fā)展水平的提高和市場化程度的擴大能促進佳木斯市農(nóng)業(yè)綠色全要素生產(chǎn)率提升。

6、灌溉效率水平——市場化程度型。此路徑下，農(nóng)業(yè)綠色全要素生產(chǎn)率主要是由高灌溉效率水平和高市場化程度共同推動。該形成路徑下，典型地區(qū)為牡丹江市。2019 年牡丹江市發(fā)布的《牡丹江市人民政府辦公室關于印發(fā)牡丹江市黑土耕地保護實施方案的通知》中強調依靠科技進步推動農(nóng)業(yè)綠色發(fā)展，改善農(nóng)業(yè)水利灌溉，提高灌溉效率水平。牡丹江市統(tǒng)計年鑒數(shù)據(jù)顯示，2012～2021 年灌溉效率水平由12.03%增長至16.33%；市場化程度由18.64%增長至30.59%。因此，灌溉效率水平的提高和市場化程度的擴大可以促進牡丹江市農(nóng)業(yè)綠色全要素生產(chǎn)率的提升。

四、主要結論及政策建議

基于SBM-GML 模型對黑龍江省13 個地級市2012～2021 年面板數(shù)據(jù)進行測算，并分析其農(nóng)業(yè)綠色全要素生產(chǎn)率的時空分異特征，通過fsQCA 方法探究農(nóng)業(yè)機械化水平、農(nóng)業(yè)經(jīng)濟發(fā)展水平、貿易開放程度、灌溉效率水平等前因變量所構建的農(nóng)業(yè)綠色全要素生產(chǎn)率提升路徑，得到如下主要結論：第一，黑龍江省農(nóng)業(yè)綠色全要素生產(chǎn)率發(fā)展呈現(xiàn)分散化的特點，主要集中于黑龍江省西部。第二，黑龍江省各地市農(nóng)業(yè)綠色全要素生產(chǎn)率的提升存在顯著的區(qū)域性差別，主要依靠經(jīng)濟發(fā)展、貿易開放、灌溉效率促進其提升。第三，黑龍江省農(nóng)業(yè)綠色全要素生產(chǎn)率共有6 條多元提升路徑，即“經(jīng)濟發(fā)展”主導“貿易開放”輔助型、“貿易開放”主導“灌溉效率”輔助型、“灌溉效率水平”主導“市場化程度”輔助型等6 條提升路徑。

基于上述研究結論，得出如下政策建議：第一，黑龍江省各地農(nóng)業(yè)綠色發(fā)展具有差異性，因地制宜，多元化施策，各地政府應著力提高農(nóng)業(yè)科技水平，進一步提高農(nóng)業(yè)生產(chǎn)技術水平與生產(chǎn)效率，解決經(jīng)濟發(fā)展與生態(tài)環(huán)境之間的矛盾關系，促進農(nóng)業(yè)綠色化發(fā)展；第二，充分發(fā)揮地區(qū)資源優(yōu)勢，優(yōu)化要素投入結構，發(fā)展地區(qū)特色農(nóng)業(yè)，推進農(nóng)業(yè)綠色產(chǎn)業(yè)轉型升級；第三，加強區(qū)域合作，優(yōu)勢互補，以農(nóng)業(yè)綠色全要素生產(chǎn)率水平較高地區(qū)帶動周邊地市，深入合作交流，實現(xiàn)區(qū)域均衡發(fā)展，整體性推動農(nóng)業(yè)綠色全要素生產(chǎn)率提高。