□于宏娟

（濟(jì)南大學(xué)，山東 濟(jì)南 250000）

近年來，黃河流域的發(fā)展質(zhì)量備受關(guān)注。產(chǎn)業(yè)轉(zhuǎn)型升級(jí)與綠色可持續(xù)發(fā)展滯后已經(jīng)成為黃河流域發(fā)展不平衡、自身發(fā)展不充分的突出表現(xiàn)[1]。綠色可持續(xù)是逐步實(shí)現(xiàn)社會(huì)經(jīng)濟(jì)與生態(tài)保護(hù)協(xié)調(diào)發(fā)展的有效路徑，是破解資源渠道、環(huán)境約束的主流趨勢(shì)，更是黃河流域轉(zhuǎn)型升級(jí)的重要方向，也必將在城市轉(zhuǎn)型升級(jí)中發(fā)揮不可替代的角色。

1 文獻(xiàn)綜述

在研究尺度上，主要從區(qū)域和行業(yè)兩個(gè)方面進(jìn)行研究，從區(qū)域視角來看，涉及全球和全國(guó)范圍。從全球范圍來看，李小勝和安慶賢（2012）[2]研究發(fā)現(xiàn)，考慮環(huán)境管制后APEC 的TFP 增長(zhǎng)水平高于不考慮環(huán)境情況下的數(shù)值。從全國(guó)范圍來看，Rusiawan W.等（2015）[3]運(yùn)用增長(zhǎng)核算法測(cè)算印尼的全要素生產(chǎn)率，發(fā)現(xiàn)全要素生產(chǎn)率對(duì)生產(chǎn)效率和碳排放有影響。張毅（2021）[4]研究非期望產(chǎn)出的超效率SBM 模型探究長(zhǎng)江經(jīng)濟(jì)帶、長(zhǎng)三角地區(qū)的綠色全要素生產(chǎn)率。此外，國(guó)內(nèi)的研究區(qū)域從全國(guó)各個(gè)省份到各個(gè)經(jīng)濟(jì)帶，比如長(zhǎng)三角、長(zhǎng)江經(jīng)濟(jì)帶、京津冀等[5-7]。從產(chǎn)業(yè)視角來看，國(guó)外學(xué)者Ball V.E.等（2004）[8]研究發(fā)現(xiàn)，美國(guó)農(nóng)業(yè)GTFP 的年均增長(zhǎng)率低于傳統(tǒng)TFP，國(guó)內(nèi)學(xué)者的研究主要涉及工業(yè)、制造業(yè)、農(nóng)業(yè)、物流業(yè)、旅游產(chǎn)業(yè)等各個(gè)產(chǎn)業(yè)。

在GTFP 測(cè)算方法上，主要有參數(shù)法和非參數(shù)法兩大類。參數(shù)法主要是指隨機(jī)前沿生產(chǎn)函數(shù)（SFA）。數(shù)據(jù)包絡(luò)分析法（DEA）是一種非參數(shù)研究方法，相對(duì)而言，該方法更加適用于效率計(jì)算，也是國(guó)內(nèi)外學(xué)者偏好的一種運(yùn)算方法。胡立和等（2020）[9]基于數(shù)據(jù)包絡(luò)分析方法和非期望產(chǎn)出的超效率SBM模型探究長(zhǎng)江經(jīng)濟(jì)帶、長(zhǎng)三角地區(qū)的綠色全要素生產(chǎn)率。

梳理上述研究結(jié)果可以發(fā)現(xiàn)，學(xué)術(shù)界鮮有對(duì)黃河流域各地級(jí)市綠色全要素生產(chǎn)率的研究。因此文章對(duì)黃河流域77 個(gè)地級(jí)市測(cè)算2008—2018 年綠色全要素生產(chǎn)率，研究黃河流域綠色全要素生產(chǎn)率的時(shí)空演化特征，具有重要意義。

2 研究方法

2.1 指標(biāo)選取與數(shù)據(jù)來源

2.1.1 投入變量

投入變量分別為黃河流域77 個(gè)地級(jí)市的資本投入、勞動(dòng)力投入和能源投入，其中，計(jì)算資本投入的根據(jù)張軍等（2004）[10]研究的永續(xù)盤存法，如式（1）。

式中：K表示資本存量，d表示折舊率，I表示新增社會(huì)固定資產(chǎn)投資，i代表城市，t代表時(shí)間。投資指標(biāo)以2008 年投資價(jià)格指數(shù)為基期平減，數(shù)據(jù)主要來源于《中國(guó)城市統(tǒng)計(jì)年鑒》。

借鑒單豪杰（2008）[11]的算法，折舊率定為10.96%。勞動(dòng)力投入是根據(jù)每個(gè)地方級(jí)城市年底的就業(yè)人數(shù)來衡量的，數(shù)據(jù)來自《中國(guó)勞動(dòng)統(tǒng)計(jì)年鑒》。在能源投入方面，文章選取了社會(huì)各界關(guān)注的二氧化碳排放量，這也是對(duì)現(xiàn)有研究的完善。

參考陳詩一（2009）[12]的做法，并根據(jù)《中國(guó)能源統(tǒng)計(jì)年鑒》可知，碳排放量=相關(guān)能源消耗量×相關(guān)碳排放系數(shù)×氧化系數(shù)×發(fā)熱量。

2.1.2 產(chǎn)出變量

產(chǎn)出變量包括期望產(chǎn)出和非期望產(chǎn)出，這與綠色經(jīng)濟(jì)發(fā)展中應(yīng)考慮的其他因素更為一致。期望產(chǎn)出采用以2008 年為基期平減的各時(shí)期實(shí)際GDP。參考張毅（2021）的研究方法，非期望產(chǎn)出以工業(yè)廢氣排放、工業(yè)廢水排放和工業(yè)煙（粉）塵排3 個(gè)指標(biāo)來衡量。選取黃河流域地區(qū)2008—2018 年的數(shù)據(jù)作為研究樣本。數(shù)據(jù)來源于《中國(guó)環(huán)境統(tǒng)計(jì)年鑒》《中國(guó)城市統(tǒng)計(jì)年鑒》，缺失數(shù)據(jù)采取插值法補(bǔ)齊。

2.2 研究區(qū)域與研究方法

2.2.1 研究區(qū)域

根據(jù)黃河流域水利委員會(huì)對(duì)黃河流域的自然劃分，黃河一共流經(jīng)9 個(gè)省（自治區(qū)），文章考慮到各個(gè)地區(qū)的經(jīng)濟(jì)關(guān)聯(lián)度，未將黃河流域流經(jīng)四川省的區(qū)域考慮進(jìn)來，去除了屬于長(zhǎng)江流域的城市（漢中市、安康市、商洛市、隴南市），借鑒師博等（2021）[13]的劃分方法，將研究區(qū)域限定在黃河流域的77 個(gè)城市，具體見表1。

表1 黃河流域77 個(gè)地級(jí)市

2.2.2 研究方法

文章參考Tone K（2002）[14]代表研究，考慮非期望產(chǎn)出和松弛投入，運(yùn)用非徑向非角度Super-SBM模型計(jì)算綠色全要素生產(chǎn)率，該模型很好地考慮了松弛變量的問題，具體假設(shè)如下。

在規(guī)模報(bào)酬不變的條件下，綠色全要素生產(chǎn)率的決策單元（DMU）有n個(gè)，記DMUj（j=1，2，…，n），每個(gè)DMU 有m種投入，記xi（i=1，2，…，m）；共得到w種產(chǎn)出，其中期望產(chǎn)出記為yr（r=1，2，…，w1），非期望產(chǎn)出記為zt（t=1，2，…，w2），則Super-SBM 模型如式（2）、式（3）。

3 黃河流域城市綠色全要素生產(chǎn)率的時(shí)空演變分析

3.1 黃河流域綠色全要素生產(chǎn)率時(shí)序演變分析

通過非期望產(chǎn)出的Super-SBM模型，測(cè)算黃河流域77 個(gè)地級(jí)市的綠色全要素生產(chǎn)率，通過算數(shù)平均計(jì)算出整體和黃河上游、中游和下游的綠色全要素生產(chǎn)率值，將其變化繪制成圖1。

從沿黃流域整體來看，整體的綠色發(fā)展不太穩(wěn)定，呈波動(dòng)狀態(tài)，但近年來有上升趨勢(shì)。

2008—2009 年，國(guó)際金融危機(jī)暴發(fā)，各國(guó)經(jīng)濟(jì)受到影響，政府采取措施導(dǎo)致資源錯(cuò)配，綠色全要素生產(chǎn)率受到影響呈下降趨勢(shì)。

“十二五”期間，2011 年我國(guó)首次提出綠色發(fā)展理念，各級(jí)政府也積極開展和落實(shí)綠色發(fā)展，因此在2011—2012 年的綠色全要素生產(chǎn)率上漲幅度較大。

2012—2015 年總體呈衰退趨勢(shì)，這與黃河流域城市的發(fā)展模式息息相關(guān)，第二產(chǎn)業(yè)占比較大，而且多為能源消耗型加工企業(yè)，這也是黃河流域綠色發(fā)展落后的原因之一。

截至2015 年，黨的十八屆五中全會(huì)明確提出五大發(fā)展理念，從2015 年開始，整個(gè)黃河流域的綠色全要素生產(chǎn)率呈上升趨勢(shì)[15]。

根據(jù)圖1 可知，對(duì)比沿黃流域上游、中游和下游3 個(gè)流域的綠色全要素生產(chǎn)率，發(fā)現(xiàn)沿黃河流域總體的發(fā)展不均衡，呈黃河上游＞黃河下游＞黃河中游的狀態(tài)。黃河上游地域?qū)拸V，土地資源較為豐富，主要發(fā)展的產(chǎn)業(yè)是農(nóng)牧業(yè)，黃河上游是重要的生態(tài)保護(hù)屏障，受國(guó)家政策支持，因此綠色全要素生產(chǎn)率較高。下游地區(qū)有著兩個(gè)省會(huì)城市，分別為濟(jì)南市、鄭州市，另外，沿海城市較多，經(jīng)濟(jì)發(fā)展水平較高，更加注重保護(hù)環(huán)境，因此綠色全要素生產(chǎn)率較為平衡[16]。中游地區(qū)能源資源較為豐富，資源利用效率較低，生態(tài)環(huán)境較為脆弱，因此綠色全要素生產(chǎn)率較低。

圖1 2008—2018 年沿黃流域綠色全要素生產(chǎn)率變動(dòng)趨勢(shì)圖

3.2 黃河流域城市綠色全要素生產(chǎn)率空間演變分析

為了進(jìn)一步探究沿黃流域77 個(gè)城市綠色發(fā)展的程度，文章采用ArcGIS 10.2 可視化黃河流域77 個(gè)地級(jí)市的綠色全要素生產(chǎn)率，并根據(jù)自然斷裂點(diǎn)分類法，將黃河流域2008 年、2018 年的綠色全要素生產(chǎn)率分為4 檔，分別為低水平區(qū)、中低水平區(qū)、中高水平區(qū)和高水平區(qū)。

由圖2 可以看到，黃河流域綠色全要素生產(chǎn)率分布不均勻，高水平區(qū)主要分布于黃河流域的邊緣地帶，說明沿黃河流域邊緣地帶的城市綠色全要素生產(chǎn)率較高，綠色發(fā)展水平優(yōu)于中間地帶。黃河流域綠色全要素生產(chǎn)率的空間分布特征在不同年份發(fā)生了較大的變化，可以看出各個(gè)年份的低水平區(qū)、中低水平區(qū)、中高水平區(qū)和高水平區(qū)的地級(jí)市數(shù)量變化，由2008 年的31 個(gè)、28 個(gè)、14 個(gè)和4 個(gè)，變?yōu)?018 年的39 個(gè)、25 個(gè)、8 個(gè)和5 個(gè)。2008 年高水平區(qū)主要集中在黃河上游，中下游發(fā)展情況不樂觀；2018 年黃河中下游綠色群要素生產(chǎn)率逐漸往中高水平發(fā)展，集中在黃河流域周邊地區(qū)。

圖2 2008 年和2018 年黃河流域綠色全要素生產(chǎn)率空間分布

4 結(jié)論與建議

4.1 結(jié)論

基于黃河流域77 個(gè)地級(jí)城市2008—2018 年的面板數(shù)據(jù)，采用Super-SBM模型測(cè)算綠色全要素生產(chǎn)率，得出結(jié)論如下。

一是從時(shí)序演變圖來看，沿黃河流域整體的綠色全要素生產(chǎn)率呈波動(dòng)趨勢(shì)，綠色發(fā)展水平不穩(wěn)定。黃河流域綠色發(fā)展水平呈黃河上游＞黃河下游＞黃河中游的態(tài)勢(shì)。

二是從空間演變圖來看，沿黃河流域整體效率不高，效率值分布不均勻，高水平區(qū)主要分布在黃河流域邊緣地帶[17]。

4.2 建議

一是堅(jiān)持落實(shí)綠色發(fā)展理念，走高質(zhì)量發(fā)展道路。黃河流域所擁有的自然資源、文化精華是其他流域不可代替的，黃河中游擁有大量資源，多數(shù)城市為能源消耗型城市，應(yīng)加快轉(zhuǎn)型，減少經(jīng)濟(jì)增長(zhǎng)對(duì)能源消耗的依賴，實(shí)現(xiàn)低碳、綠色發(fā)展，促進(jìn)當(dāng)?shù)亟?jīng)濟(jì)快速增長(zhǎng)[18]。

二是要進(jìn)一步實(shí)現(xiàn)經(jīng)濟(jì)增長(zhǎng)，充分發(fā)揮中心城市輻射作用，帶動(dòng)黃河中下游城市發(fā)展，推動(dòng)區(qū)域經(jīng)濟(jì)一體化進(jìn)程。針對(duì)不同類型的城市采取不同的措施，濟(jì)南、西安、鄭州等地經(jīng)濟(jì)發(fā)展較好的城市要加強(qiáng)研發(fā)綠色科技，優(yōu)化升級(jí)產(chǎn)業(yè)結(jié)構(gòu)，助力黃河流域發(fā)展。

三是要大力發(fā)展數(shù)字經(jīng)濟(jì)，以數(shù)字經(jīng)濟(jì)推動(dòng)生產(chǎn)和消費(fèi)等領(lǐng)域的發(fā)展，利用數(shù)字化和智能化的手段改造資源消耗型城市，優(yōu)化資源配置效率，進(jìn)一步提升黃河流域城市綠色發(fā)展水平[19]。