黑龍江省農(nóng)業(yè)綠色全要素生產(chǎn)率提升研究

［提要］基于SBM-GML模型測算黑龍江省13個地級市的農(nóng)業(yè)綠色全要素生產(chǎn)率，并分析其時空分異特征。通過fsQCA方法探究農(nóng)業(yè)綠色全要素生產(chǎn)率提升路徑，結(jié)果表明：黑龍江省農(nóng)業(yè)綠色全要素生產(chǎn)率發(fā)展呈現(xiàn)分散化的特點；主要有“經(jīng)濟發(fā)展”主導“貿(mào)易開放”輔助型、“貿(mào)易開放”主導“灌溉效率”輔助型、“灌溉效率水平”主導“市場化程度”輔助型等六條提升路徑。對此，各地政府應(yīng)著力提高農(nóng)業(yè)科技水平，充分發(fā)揮地區(qū)資源優(yōu)勢，優(yōu)化要素投入結(jié)構(gòu)，推進農(nóng)業(yè)綠色產(chǎn)業(yè)轉(zhuǎn)型升級，加強區(qū)域合作實現(xiàn)區(qū)域均衡發(fā)展，推動提高農(nóng)業(yè)綠色全要素生產(chǎn)率。黨的二十大報告明確提出要加快發(fā)展方式綠色轉(zhuǎn)型，積極穩(wěn)妥推進碳達峰、碳中和。農(nóng)業(yè)在實現(xiàn)經(jīng)濟社會發(fā)展綠色化、低碳化過程中具有不可替代的作用，著力提升農(nóng)業(yè)綠色發(fā)展水平將成為未來一段時期農(nóng)業(yè)發(fā)展的核心任務(wù)。黑龍江省作為中國糧食安全“壓艙石”，正面臨水土資源過度消耗、黑土地質(zhì)量下降風險加大以及農(nóng)業(yè)面源污染等問題，制約著黑龍江省農(nóng)業(yè)綠色轉(zhuǎn)型。由于農(nóng)業(yè)綠色全要素生產(chǎn)率能夠全方位衡量區(qū)域農(nóng)業(yè)綠色轉(zhuǎn)型升級的綜合水平，因此研究黑龍江省農(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ù)包絡(luò)分析技術(shù)（DEA）為基礎(chǔ)的非參數(shù)方法為主。國內(nèi)相關(guān)研究較多使用的是基于徑向、角度的DEA模型。傳統(tǒng)的全要素生產(chǎn)率測算方法沒有考慮農(nóng)業(yè)經(jīng)濟增長的環(huán)境代價，其測算結(jié)果無法真實反映環(huán)境約束下的農(nóng)業(yè)生產(chǎn)效率，如徑向DEA模型無法處理投入冗余和產(chǎn)出不足產(chǎn)生的松弛變量問題，角度DEA模型需要假定投入或產(chǎn)出不變，這與現(xiàn)實并不完全相符。為克服傳統(tǒng)DEA模型的缺陷，Chung提出了方向性距離函數(shù)的概念，據(jù)此構(gòu)建了Malmquist-Luenberger（ML）生產(chǎn)率指數(shù)，有效解決了期望產(chǎn)出和非期望產(chǎn)出聯(lián)合生產(chǎn)的環(huán)境效率評價問題。國內(nèi)學者們則構(gòu)建了非徑向、非角度SBM模型，如李谷成、閆桂權(quán)使用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)率的影響存在顯著的雙重門檻效應(yīng)。這些研究成果已經(jīng)成為測算農(nóng)業(yè)綠色全要素生產(chǎn)率的重要依據(jù)。（二）農(nóng)業(yè)綠色全要素生產(chǎn)率的影響因素。關(guān)于農(nóng)業(yè)綠色全要素生產(chǎn)率的影響因素，現(xiàn)有文獻主要從技術(shù)因素、制度因素、市場因素和農(nóng)業(yè)生產(chǎn)條件等方面展開研究。技術(shù)支撐是農(nóng)業(yè)綠色發(fā)展的基礎(chǔ)，其支撐作用在實證研究中得到了證實。如，Reinhard認為技術(shù)條件既可以直接影響農(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ā)展的重要因素。相關(guān)研究認為農(nóng)業(yè)市場化有助于改善農(nóng)村生態(tài)環(huán)境質(zhì)量。Plassmann強調(diào)人們必須重視市場需求在減少農(nóng)業(yè)污染方面的重要作用，否則，最復雜、最先進的技術(shù)也無法阻止農(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)狀提供了理論和方法。學術(shù)界從多方面研究了農(nóng)業(yè)綠色全要素生產(chǎn)率的影響因素，為本研究構(gòu)建理論框架提供了重要依據(jù)。但已有研究主要探究了單個影響因素對農(nóng)業(yè)綠色全要素生產(chǎn)率的凈效應(yīng)，忽視了影響因素之間的協(xié)同聯(lián)動對農(nóng)業(yè)綠色全要素生產(chǎn)率的影響。多因素聯(lián)動作用的研究有助于識別出農(nóng)業(yè)綠色全要素生產(chǎn)率的多種提升路徑，深化農(nóng)業(yè)綠色全要素生產(chǎn)率提升背后的多因素共同作用的本質(zhì)認識。二、數(shù)據(jù)來源與指標選?。ㄒ唬?shù)據(jù)來源。本文選取2012～2021年黑龍江省13個地級市的面板數(shù)據(jù)。用于測算黑龍江省各地市農(nóng)業(yè)綠色全要素生產(chǎn)率的相關(guā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ā)展水平、貿(mào)易開發(fā)程度、市場化程度、灌溉效率水平等5個變量作為農(nóng)業(yè)綠色全要素生產(chǎn)率的前因變量。將農(nóng)業(yè)綠色全要素生產(chǎn)率作為結(jié)果變量。主要變量說明見表2。（表2）2、變量校準。采用分位數(shù)值來確定錨點的直接法對結(jié)果與前因條件進行校準。結(jié)果與前因條件均采用直接校準法，所有連續(xù)變量的完全隸屬、交叉臨界點和完全不隸屬的三個閾值，在此分別設(shè)置為95%、50%和5%三個分位數(shù)值，若對校準數(shù)據(jù)中出現(xiàn)0.5數(shù)值的條件，將其數(shù)值替換處理為0.501，以避免對0.5數(shù)值所屬案例進行理論分析的困難。校準結(jié)果如表3所示。（表3）（四）農(nóng)業(yè)綠色全要素生產(chǎn)率測算方法。本文使用SBM-GML模型測算黑龍江省農(nóng)業(yè)綠色全要素生產(chǎn)率。1、超效率SBM模型。為有效解決農(nóng)業(yè)碳排放在測算過程中對效率的影響，采用Tone提出的超效率SBM模型，相關(guān)SBM模型測算方法已得到有效驗證，參考焦兵、王騰、王兵等相關(guān)研究。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)率、技術(shù)進步、技術(shù)效率。三、結(jié)果與分析（一）測算結(jié)果分析。測算并分解黑龍江省13個地級市2012～2021年農(nóng)業(yè)綠色全要素生產(chǎn)率，結(jié)果如表4所示。農(nóng)業(yè)綠色全要素生產(chǎn)率整體趨于平穩(wěn)，但在個別年份間存在波動，其中2015～2016年農(nóng)業(yè)綠色全要素生產(chǎn)率為0.948，與上一年相比下降11.81%。盡管技術(shù)進步較上一年提高了7.96%，對農(nóng)業(yè)綠色全要素生產(chǎn)率產(chǎn)生正向影響，但是由于純技術(shù)效率和規(guī)模效率分別下降4.35%和14.63%，最終導致農(nóng)業(yè)綠色全要素生產(chǎn)率下降。（表4）由表5可知，橫向?qū)Ρ葋砜?，黑龍江省?nèi)13個地級市農(nóng)業(yè)綠色全要素生產(chǎn)率存在顯著的區(qū)域性差異。如，綏化市和雞西市等地農(nóng)業(yè)綠色全要素生產(chǎn)率的提升主要依靠強勁的技術(shù)效率和技術(shù)進步，其中鶴崗市純技術(shù)效率高達1.077；然而，齊齊哈爾市和大慶市主要是由規(guī)模效率驅(qū)動的農(nóng)業(yè)綠色全要素生產(chǎn)率提升。（表5）表5農(nóng)業(yè)綠色全要素生產(chǎn)率的空間差異一覽表（二）必要性條件分析。為探求單個因素是否為結(jié)果的必要性條件，借鑒Ragin的方法，對前因變量進行必要性檢驗，是否產(chǎn)生以結(jié)果集合作為條件集合（或者其非集）的子集，必要性檢驗結(jié)果如表6所示。單個因素的一致性水平均不超過0.9，因此單個因素均不是必要性條件，表明其提升路徑受多因素的組態(tài)效應(yīng)影響。（表6）表6單個因素必要性條件檢驗結(jié)果一覽表（三）條件組態(tài)結(jié)果。為得到提升路徑條件組態(tài)，借鑒張明、杜運周等方法，將初始一致性閾值設(shè)為0.7，PRI不低于0.7，頻數(shù)閾值設(shè)為1，得到黑龍江13個地級市農(nóng)業(yè)綠色全要素生產(chǎn)率相應(yīng)組態(tài)，存在6條組態(tài)增長路徑。第1條提升路徑：“經(jīng)濟發(fā)展”主導“貿(mào)易開放”輔助型；第2條提升路徑：“經(jīng)濟發(fā)展”主導“灌溉效率”輔助型；第3條提升路徑：“貿(mào)易開放”主導“灌溉效率”輔助型；第4條提升路徑：“對外開放”主導“農(nóng)業(yè)機械化+市場化程度”輔助型；第5條提升路徑：“經(jīng)濟發(fā)展水平”主導“市場化程度”輔助型；第6條提升路徑：“灌溉效率水平”主導“市場化程度”輔助型。這6條提升路徑的一致性都高于0.8，高于初設(shè)臨界值0.7，總體覆蓋度0.636，可以解釋超63%的案例，說明6條路徑的解釋力度均較好。（表7）表7組態(tài)分析結(jié)果一覽表（四）條件組態(tài)分析1、經(jīng)濟發(fā)展——貿(mào)易開放型?；诖寺窂较?，即農(nóng)業(yè)機械化程度低、市場資源配置結(jié)構(gòu)未得到優(yōu)化，依靠較高的經(jīng)濟發(fā)展水平和貿(mào)易開放促進農(nóng)業(yè)綠色全要素生產(chǎn)率提高。解釋典型地區(qū)為伊春市、鶴崗市。以伊春市為例，經(jīng)濟發(fā)展水平對農(nóng)業(yè)生產(chǎn)技術(shù)進步、農(nóng)業(yè)生產(chǎn)結(jié)構(gòu)等方面具有相當?shù)挠绊懥?，?jīng)濟發(fā)展水平較高，意味著其資源配置效率較高，農(nóng)業(yè)基礎(chǔ)設(shè)施更加完善等，能夠有效提高農(nóng)村居民純收入；伊春市嘉蔭口岸和市政府推進貿(mào)易開放的措施政策等提高貿(mào)易開放程度。2012～2021年伊春市農(nóng)村居民純收入由0.92萬元增至1.94萬元，增長110.87%，農(nóng)民有較高的收入意味著農(nóng)民所掌握的生產(chǎn)技術(shù)較高，有著較高的生活質(zhì)量，在農(nóng)業(yè)生產(chǎn)、提升生產(chǎn)技術(shù)水平等方面分配更多的時間精力；貿(mào)易開放在需求端刺激農(nóng)業(yè)綠色生產(chǎn)率的提高，更為嚴格的農(nóng)產(chǎn)品出口標準在一定程度上引導農(nóng)業(yè)向著綠色化和標準化發(fā)展，從而較高的農(nóng)業(yè)經(jīng)濟發(fā)展水平和貿(mào)易開放程度的提高能夠促進伊春農(nóng)業(yè)綠色全要素生產(chǎn)率的提升。2、經(jīng)濟發(fā)展——灌溉效率型。此路徑下，農(nóng)業(yè)綠色全要素生產(chǎn)率主要是由高經(jīng)濟發(fā)展水平主導和良好的灌溉效率水平輔助共同推動。該路徑中的典型地區(qū)為大慶市。以大慶市為例，大慶市在2019年印發(fā)《大慶市市級農(nóng)民專業(yè)合作社示范社評定管理辦法》，著重提出加快構(gòu)建現(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、貿(mào)易開放——灌溉效率型。此路徑下，在農(nóng)業(yè)機械化水平較為落后和市場資源配置結(jié)構(gòu)未得到優(yōu)化的情況下，高貿(mào)易開放程度和灌溉效率水平的提高能夠有效促進農(nóng)業(yè)綠色全要素生產(chǎn)率的提高。該路徑中的典型地區(qū)為齊齊哈爾市。以齊齊哈爾市為例，在外貿(mào)方面，打造國際貿(mào)易供給力、營造對外貿(mào)易健康發(fā)展環(huán)境等全面擴大貿(mào)易開放程度；實施高效節(jié)水灌溉、節(jié)水增糧、灌溉基礎(chǔ)設(shè)施建設(shè)等一系列措施有效提高灌溉效率水平。貿(mào)易開放主導、灌溉效率輔助構(gòu)成農(nóng)業(yè)綠色全要素生產(chǎn)率提升路徑。4、貿(mào)易開放——農(nóng)業(yè)機械化+市場化程度型。此路徑下，即使經(jīng)濟發(fā)展水平相對較低，但是在高農(nóng)業(yè)機械化水平和高市場化程度因素的推動下，依靠高貿(mào)易開放程度仍能促進農(nóng)業(yè)綠色全要素生產(chǎn)率提高。該路徑下，典型地區(qū)為黑河市。以黑河市為例，2019年國務(wù)院印發(fā)《中國（黑龍江）自由貿(mào)易試驗區(qū)總體方案》，明確提出黑龍江自貿(mào)試驗區(qū)各片區(qū)產(chǎn)業(yè)發(fā)展的功能布局和重點領(lǐng)域，自貿(mào)區(qū)黑河片區(qū)投入建設(shè)以來，以多元聯(lián)動發(fā)展為目標，形成數(shù)字自貿(mào)區(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è)機械化水平和市場化程度因素推動下，貿(mào)易開放程度的提高能促進黑河市農(nóng)業(yè)綠色全要素生產(chǎn)率的提升。5、經(jīng)濟發(fā)展水平——市場化程度型。此路徑下，即使農(nóng)業(yè)機械化水平程度、貿(mào)易開放程度和灌溉效率水平低，但由于高經(jīng)濟發(fā)展水平和高市場化程度因素的助推，農(nóng)業(yè)綠色全要素生產(chǎn)率仍能得到提升。該形成路徑下，典型地區(qū)為佳木斯市。佳木斯市在《佳木斯市促進平臺經(jīng)濟高質(zhì)量發(fā)展實施意見》中強調(diào)要抓住新經(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ā)布的《牡丹江市人民政府辦公室關(guān)于印發(fā)牡丹江市黑土耕地保護實施方案的通知》中強調(diào)依靠科技進步推動農(nóng)業(yè)綠色發(fā)展，改善農(nóng)業(yè)水利灌溉，提高灌溉效率水平。牡丹江市統(tǒng)計年鑒數(shù)據(jù)顯示，2012～2021年灌溉效率水平由12.03%增長至16.33%；市場化程度由18.64%增長至30.59%。因此