張 碩，武媛媛，武子健

（1.河北工業(yè)大學 經濟管理學院，天津300130；2.吉林大學 生物與農業(yè)工程學院，長春130012）

一、引言

黨的十九大提出，中國經濟已經由高速增長階段轉向高質量發(fā)展階段，高質量發(fā)展是貫徹“創(chuàng)新、協調、綠色、開放、共享”新發(fā)展理念的根本體現，而綠色發(fā)展則是建設現代化經濟體系的必然要求。“綠色經濟”一詞早在1989 年被英國經濟學家David Pearce 提出，此后學者們展開對綠色發(fā)展的研究。[1][2]近年來，諸多學者圍繞區(qū)域綠色發(fā)展進行了較為深入的研究。在綠色發(fā)展的內涵方面，胡鞍鋼和周紹杰認為綠色發(fā)展觀是第二代可持續(xù)發(fā)展觀，強調經濟系統(tǒng)、自然系統(tǒng)和社會系統(tǒng)間的系統(tǒng)性、整體性和協調性。[3]車磊等認為綠色發(fā)展的核心理念在于低碳、循環(huán)和可持續(xù)發(fā)展，其關鍵在于提升綠色發(fā)展效率。[4]在綠色發(fā)展效率測度方面，隨機前沿方法、數據包絡分析法等被學者們廣泛運用。[5]在指標體系構建方面，資本、勞動力、能源消耗多被納入指標體系，產出指標則多選擇國內生產總值來衡量。[6][7]

河北省環(huán)抱北京、天津兩大直轄市，是中國北方重要省份，承擔了疏解北京非首都核心功能、擴大環(huán)境容量生態(tài)空間等重要任務。然而，當前對河北省綠色發(fā)展的研究較少，缺乏理論和實踐支撐。有鑒于此，現先運用熵值法計算出工業(yè)環(huán)境成本指數，構建河北省綠色經濟效率評價指標體系，最后運用DEA 模型和Malmquist 指數測度2010—2018年河北省11 個地級市綠色經濟效率，并對其時空演化特征進行分析，以期為河北省可持續(xù)發(fā)展提供參考。

二、綠色經濟效率指標體系構建與處理

（一）投入指標選取與數據處理

現以投入—產出理論為基礎，并考慮數據的科學性與可得性，構建河北省綠色經濟效率評價指標體系，見表1。

表1 河北省綠色經濟效率評價指標體系

在勞動力投入方面，選取城鎮(zhèn)單位從業(yè)人員數作為勞動力就業(yè)人數。在資本投入方面，采用永續(xù)盤存法估算河北省各地歷年資本存量，以2010 年為基期，用2010 年固定資產投資額除以2010—2018 年固定資產投資額的年平均增長率與折舊率之和作為2010 年的資本存量。其中，折舊率δ 取9.5%，每年的固定資產投資額均通過固定資產投資指數、價格指數先折算到基期，具體計算步驟與公式如下。

式（1）中，K0為基期固定資本存量，I0為基期固定資產投資額，g 為固定資產價格年平均增長率，δ為固定資產折舊率。

式（2）中，Kt為某地第t 年的固定資本存量，Kt-1為該地第t-1 年的固定資本存量，δ 為折舊率9.5%，It為該地第t 年的固定資產投資額，已折算到基期。

在資源投入方面，基于數據的可得性和科學性，選取了工業(yè)用電瓦數作為衡量資源投入的標準。

（二）產出指標選取與數據處理

在產出方面，引用張軍和佃杰采取的綠色GDP計算方法（綠色GDP=GDP/環(huán)境成本綜合指數），[8]其反映的是單位環(huán)境成本帶來的經濟產出值，若綠色GDP 值越大，則經濟增長質量越高，若綠色GDP值越小，則經濟增長質量越低，GDP 均折算到基期。

基于數據的可獲得性，現選取各市工業(yè)廢水排放量、工業(yè)煙塵排放量、工業(yè)二氧化硫排放量作為基本指標，運用熵值法構建工業(yè)環(huán)境成本指數，具體計算步驟與公式如下。

首先，對各市工業(yè)廢水排放量、工業(yè)煙塵排放量、工業(yè)二氧化硫排放量指標進行無量綱化處理，設某城市某指標為xij，標準化后值為yij，則：

其次，運用熵值法確定各指標權重。因為熵值法要涉及到對數計算，而在無量綱處理時存在標準化后值為0 的情況，為避免對數計算無意義，先對數據進行處理，為保證平移后對原始數據影響盡可能小，故選擇向右平移0.0001 個單位，公式為：

計算第j 個指標下i 地指標比重：

計算第j 個指標的熵值：

計算第j 個指標的權重：

最后，計算綜合指數：

（三）數據來源

現所需的統(tǒng)計數據均來源于2011—2019 年的《中國統(tǒng)計年鑒》《中國城市統(tǒng)計年鑒》《河北經濟年鑒》，對于缺失的2018 年固定資產投資額數據，根據《中國城市統(tǒng)計年鑒2019》發(fā)布的各市環(huán)比增長率計算后補充。2016 年石家莊工業(yè)用電數據缺失，采用前后年份平均值進行補齊。

三、發(fā)展評價與演化分析

（一）河北省工業(yè)環(huán)境成本指數變化情況

現選取工業(yè)廢水排放量、工業(yè)煙塵排放量和工業(yè)二氧化硫排放量作為指標，運用熵值法計算得出2010—2018 年河北省城市工業(yè)環(huán)境成本指數，如表2 所示。

經濟發(fā)展最好的唐山市環(huán)境成本指數最高，說明該市發(fā)展對環(huán)境的外部性較大。排名比較靠前的邯鄲、石家莊等城市也大致與經濟發(fā)展狀況相符，說明與省內城市比較，這些城市的發(fā)展在一定程度上仍是以犧牲環(huán)境為代價的。廊坊的工業(yè)環(huán)境成本排在第十位，但2018 年國內生產總值排在全省第六位，這與其地理位置、產業(yè)結構、京津冀協同發(fā)展政策等因素相關。

表2 2010—2018 年河北省城市工業(yè)環(huán)境成本指數及排名

從時間維度上來看，除個別年份外，河北省各城市工業(yè)環(huán)境成本指數的變化呈現出很強的一致性，經濟較發(fā)達的城市工業(yè)環(huán)境成本指數一直在全省處于靠前地位，具體如圖1 所示。

圖1 2010—2018 年河北省各城市工業(yè)環(huán)境成本指數變化趨勢圖

（二）基于Malmquist 指數的效率分解

在考察經濟效率方面，DEA 模型被學者廣泛采用?，F借助DEAP2.1 軟件，計算2010—2018 年河北省11 個城市的Malmquist 生產率及其分解指數。假設規(guī)模報酬可變，投入導向型。得到2017—2018 年河北省綠色經濟效率的Malmquist 生產率及其分解指數，如表3 所示。

表3 2017—2018 年河北省綠色經濟效率情況

從城市層面來看，共有四個城市的技術效率變化高于全省平均水平，廊坊市最高，七個城市低于全省平均水平，說明各城市之間綠色經濟效率差異明顯，在現有技術條件下，資源配置能力普遍較低。各城市純技術效率差異不大，集中在0.96—1.2 之間，其中，有七個地級市達到或超過1。在規(guī)模效率方面，河北省有三個地級市超過了1，廊坊最高為2.285，有八個地級市接近1，說明河北省各城市規(guī)模效應發(fā)揮較好。

為了進一步考察更長時期內各市綠色經濟效率的區(qū)域差異，得到河北省11 城市2010—2018 年技術效率變化、技術進步、純技術效率、規(guī)模效率、全要素生產率，如表4 所示。

表4 河北省各城市綠色經濟效率情況

根據表4 中的數據可知，2010—2018 年河北省全要素生產率變化為-3.9%，說明生產力水平和生產效率有下降情況。其中，技術效率變化指數為0.542，技術進步變化指數為1.786，說明技術進步貢獻了大部分的全要素生產率。從地級市層面來看，大部分城市全要素生產率均小于1，表明各城市綠色發(fā)展效率尚有較大提升空間。其中衡水最高，表明衡水市發(fā)展付出的環(huán)境代價最小，效率最高。部分城市的全要素生產率變化與地理位置、政策和產業(yè)結構緊密相關。從技術進步變化指數來看，各地級市均超過1。從純技術效率來看，各城市差異較小，石家莊、承德、張家口、秦皇島、唐山、邢臺、衡水均達到或大于1，廊坊、保定、滄州、邯鄲小于1，生產規(guī)模尚有調整的空間。綜合以上分析可知，技術進步貢獻了河北省大部分的全要素生產率，純技術效率和規(guī)模效率不高是全要素生產率下降的最主要原因。

（三）時空演化

從時間維度來看，近年來河北省全要素生產率變化不大，在1 上下浮動，如表5 所示，圖2 中的折線也說明了這一趨勢。由于全要素生產率在1 上下浮動，導致技術效率變化指數和技術進步變化指數呈現關于1 對稱的分布形態(tài)，二者在2015—2016年呈現出較大差異。

表5 2010—2018 年河北省綠色經濟效率變動情況

圖2 2011—2018 年河北省綠色經濟效率Malmquist指數分解結果

進一步分解技術效率變化指數來考察，如圖3 所示。純技術效率在1 上下浮動，變化幅度不大。技術效率變化和規(guī)模效率變化趨勢具有很明顯的一致性，二者呈現出先上升再下降，又上升的變化趨勢。

圖3 2011—2018 年河北省綠色經濟發(fā)展技術效率變化指數分解

為了從空間維度考察河北省綠色經濟效率變化規(guī)律與區(qū)域差異，現借助ArcGIS10.6 軟件，選取2011—2012 年、2014—2015 年、2017—2018 年作為考察時間節(jié)點，繪制河北省地級市綠色經濟技術效率、純技術效率和規(guī)模效率空間分布圖。其中，運用自然斷點法將效率值劃分為高水平、中等水平、較低水平和低水平四類。

圖4 河北省綠色經濟技術效率分布格局及演化

圖4為河北省綠色經濟技術效率分布，2011—2012 年技術效率高水平城市為張家口和衡水，技術效率低水平城市為石家莊和邢臺；2014—2015 年技術效率高水平城市為石家莊和邯鄲，技術效率低水平城市為衡水和秦皇島；2017—2018 年技術效率高水平城市為廊坊，技術效率低水平城市為滄州。總體來看，河北省綠色經濟技術效率變化規(guī)律性不強，冀北地區(qū)大致處于高水平或中等水平，中部城市處于較低水平，總體呈現向好發(fā)展趨勢。

圖5 為河北省純技術效率分布，2011—2012 年純技術效率高水平的城市為石家莊，純技術效率低水平城市為邢臺和滄州；2014—2015 年純技術效率高水平城市為滄州，純技術效率低水平城市為承德、廊坊、石家莊和衡水；2017—2018 年純技術效率高水平城市為廊坊，純技術效率低水平城市為邢臺、邯鄲和滄州?？傮w來看，近年來河北省綠色經濟純技術效率整體呈現出下降趨勢。

圖5 河北省純技術效率分布格局及演化

圖6 河北省規(guī)模效率分布格局及演化

圖6為河北省綠色經濟規(guī)模效率分布，2011—2012 年規(guī)模效率高水平的城市為張家口和衡水，規(guī)模效率低水平城市為石家莊和邢臺；2014—2015 年規(guī)模效率高水平的城市有邯鄲，規(guī)模效率低水平城市有秦皇島；2017—2018 年規(guī)模效率高水平的城市有廊坊，規(guī)模效率低水平城市有滄州。整體來看，河北省北部區(qū)域規(guī)模效率水平較高，中部和南部地區(qū)處水平較低。

四、結論與對策建議

現運用DEA—Malmquist 指數測度了2010—2018 年河北省地級市綠色經濟效率和全要素指數的變化，并在此基礎上分析了其現狀和時空變化趨勢。研究發(fā)現，一是經濟發(fā)展水平高的地區(qū)存在環(huán)境成本高的情況，說明河北省當前仍存在以環(huán)境為代價謀發(fā)展的現象，尚未實現高質量發(fā)展。二是河北省綠色經濟效率發(fā)展水平不高，仍有較大提升空間。三是技術進步貢獻了大部分的全要素生產率，而規(guī)模效率和技術效率貢獻不足，仍需提高規(guī)模效應，提高資源配置效率。

鑒于此研究結論，針對綠色發(fā)展提出以下建議。一是積極融入京津冀協同發(fā)展，環(huán)北京區(qū)域發(fā)揮好相應作用。繼續(xù)建設張家口“京張體育文化旅游休閑帶”，推動“旅游+”融合發(fā)展。立足“首都水源涵養(yǎng)功能區(qū)和生態(tài)環(huán)境支撐區(qū)”功能定位，著力推動交通、生態(tài)、產業(yè)一體化建設。二是貫徹落實“創(chuàng)新、協調、綠色、開放、共享”新發(fā)展理念。推動唐山、石家莊等經濟較發(fā)達地區(qū)產業(yè)結構升級，優(yōu)化和完善環(huán)境規(guī)制政策，將生態(tài)指標納入政府績效考核，多主體參與環(huán)境治理，實現高質量發(fā)展。推動滄州、邯鄲等地積極承接北京產業(yè)轉移，提高資源配置效率和技術效率，進而提高全要素生產率。三是深入實施創(chuàng)新驅動發(fā)展戰(zhàn)略。大力支持技術創(chuàng)新，發(fā)展優(yōu)勢產業(yè)，打造創(chuàng)新集群，營造良好的創(chuàng)新生態(tài)。推動大數據、人工智能等領域創(chuàng)新發(fā)展，促進行業(yè)數字化加速轉型，提高效率和優(yōu)化經濟結構。充分發(fā)揮地理位置優(yōu)勢，打造跨區(qū)域優(yōu)勢產業(yè)鏈。圍繞雄安新區(qū)集聚創(chuàng)新資源、提升創(chuàng)新能力、打造創(chuàng)新平臺、建設創(chuàng)新高地。四是因地制宜，發(fā)揮各市特色和資源優(yōu)勢促進高質量發(fā)展，各地應發(fā)揮本地發(fā)展?jié)摿?，實現優(yōu)勢互補。