鄭寶華　劉東皇

摘 要：基于2001-2015年的省級面板數(shù)據(jù)，在利用曼奎斯特指數(shù)計算工業(yè)TFP基礎(chǔ)上，構(gòu)建空間計量模型，實(shí)證分析了我國工業(yè)TFP對區(qū)域環(huán)境污染的空間溢出效應(yīng)。研究發(fā)現(xiàn)，中國省域環(huán)境污染存在顯著的正向空間相關(guān)性，區(qū)域環(huán)境污染程度受其他與之相近空間地區(qū)的影響；工業(yè)TFP的提升對區(qū)域內(nèi)的環(huán)境污染改善有顯著促進(jìn)作用，對區(qū)域外的環(huán)境污染改善溢出效應(yīng)亦顯著，大力促進(jìn)創(chuàng)新驅(qū)動，改變經(jīng)濟(jì)增長方式對區(qū)域環(huán)境生態(tài)有明顯的改善作用；另外經(jīng)濟(jì)發(fā)展水平向好對區(qū)域內(nèi)環(huán)境污染改善有顯著促進(jìn)作用；第二產(chǎn)業(yè)占比對本地區(qū)和相鄰地區(qū)環(huán)境污染程度增加具有顯著的推動作用；外商投資水平有利于本地區(qū)的環(huán)境改善。

關(guān)鍵詞： 工業(yè)全要素生產(chǎn)率； 空間溢出； 環(huán)境污染

中圖分類號：F424.6 文獻(xiàn)標(biāo)識碼：A 文章編號：2095-7394（2018）01-0079-09

我國工業(yè)經(jīng)歷了長達(dá)三十多年的高速發(fā)展，增長總量舉世矚目。但是以往工業(yè)經(jīng)濟(jì)增長主要是依靠要素的大量投入獲得的，而不是主要依靠全要素生產(chǎn)率（TFP）的提升，這種主要依靠大量要素投入的工業(yè)增長帶來了嚴(yán)重的環(huán)境污染問題，例如多地出現(xiàn)的霧霾現(xiàn)象和大區(qū)域的水污染等，因此，這種犧牲環(huán)境的增長模式不可持續(xù)。

當(dāng)前我國經(jīng)濟(jì)增長方式正從要素驅(qū)動向創(chuàng)新驅(qū)動轉(zhuǎn)變，更加重視生態(tài)文明建設(shè)和全要素增長率的提升，而其中全要素增長率之所以得到重視，一方面，是因?yàn)槲覈?jīng)濟(jì)已經(jīng)面臨中等收入階段，社會也進(jìn)入老齡化，人口紅利將逐漸消失，經(jīng)濟(jì)增長不能只靠資本等要素投入，需要向依托技術(shù)進(jìn)步的內(nèi)生增長模式轉(zhuǎn)型[1]；另一方面，TFP在經(jīng)濟(jì)增長中比例的提升，也是技術(shù)效率提升與技術(shù)進(jìn)步的結(jié)果，通過技術(shù)效率提升和技術(shù)進(jìn)步可以減少能源、物質(zhì)等要素的投入，減輕環(huán)境污染，工業(yè)產(chǎn)品制造過程中釋放的污染物是環(huán)境污染的主要來源，在保持工業(yè)經(jīng)濟(jì)增長的同時，工業(yè)TFP的提升有利于環(huán)境改善。在此情形下，研究工業(yè)TFP和環(huán)境污染之間的關(guān)系顯得較有必要。

一、文獻(xiàn)綜述

所謂全要素生產(chǎn)率（TFP）是指在經(jīng)濟(jì)增長的貢獻(xiàn)來源中，一部分來自資本、勞動力等投入要素的增長，余下的則歸因?yàn)榧夹g(shù)進(jìn)步和效率的改善，這一部分即為反映經(jīng)濟(jì)增長質(zhì)量的全要素生產(chǎn)率（TFP），在很大程度上決定了一個國家或者地區(qū)的經(jīng)濟(jì)增長質(zhì)量。目前關(guān)于TFP的研究主要集中于以下兩個方面：一是關(guān)于TFP的測算，張軍、施少華（2003）[2]、郭慶旺、賈俊雪（2005）[3]、傅曉霞、吳利學(xué)（2006）[4]、李靜、孟令杰等（2006）[5]、李征（2016）[6]等都用各自的方法測算了各個階段的區(qū)域全要素生產(chǎn)率及特征；Huang et al（1997）[7]、科埃利（Coelli）、普拉薩達(dá)·拉奧（Prasada Rao）（2003）[8]、朱鐘棣、李小平（2005）[9]、涂正革（2007）[10]、楊勇（2008）[11]、劉興凱、張誠（2010）[12]、楊汝岱（2015）[1]則對部門或制造業(yè)等行業(yè)的TFP進(jìn)行了研究。二是關(guān)于TFP的影響因素分析，米勒（Miller）（2000）[13]考察了貿(mào)易開放度和人力資本對全要素生產(chǎn)率的影響；顏鵬飛、王兵（2004）[14]認(rèn)為人力資本和制度因素對全要素生產(chǎn)率的提高均有重要影響；趙偉、張萃（2008）[15]研究了我國制造業(yè)集聚對全要素生產(chǎn)率的影響；劉秉鐮等（2010）[16]采用空間面板計量方法研究了我國交通基礎(chǔ)設(shè)施與全要素生產(chǎn)率之間的關(guān)系，蔣殿春、王曉嬈（2015）[17]研究了我國 R&D 結(jié)構(gòu)對全要素生產(chǎn)率的影響。

關(guān)于環(huán)境的研究一直離不開經(jīng)濟(jì)增長因素。李國璋、孔令寬（2008）[18]對中國環(huán)境庫茲涅茨曲線是否存在及其形狀進(jìn)行了檢驗(yàn)；彭水軍、包群（2006）[19]，曹光輝等（2006）[20]等探究了中國環(huán)境質(zhì)量改善的拐點(diǎn)；還有些學(xué)者檢驗(yàn)了具體的環(huán)境效應(yīng)，如經(jīng)濟(jì)規(guī)模的增長、結(jié)構(gòu)的變化、對外貿(mào)易、FDI以及技術(shù)的進(jìn)步等對環(huán)境質(zhì)量的影響情況，黃菁（2009） [21]運(yùn)用改進(jìn)的Divisia指數(shù)分解方法分析我國工業(yè)污染物變化的各種影響效應(yīng)；鄧柏盛、宋德勇（2008）[22]研究了對外貿(mào)易對環(huán)境質(zhì)量的影響；白俊紅、呂曉紅（2015）[23]研究了FDI質(zhì)量和環(huán)境污染改善之間的關(guān)系。

目前關(guān)于工業(yè)TFP與環(huán)境污染之間的研究并不多，有一些學(xué)者運(yùn)用新的研究方法，將環(huán)境作為要素投入或者產(chǎn)出，測算環(huán)境TFP或者低碳TFP，如楊俊、邵漢華（2009） [24]；吳軍、笪鳳媛、張建華（2010） [25]；鄭寶華、謝忠秋（2011）[26]等。聶國卿、尹向飛等（2010）[27]用湖南省數(shù)據(jù)研究了工業(yè)TFP對環(huán)境污染的影響。總體來看，關(guān)于工業(yè)TFP與環(huán)境污染之間的研究目前主要存在兩個局限。一是研究工業(yè)TFP對環(huán)境污染之間影響關(guān)系的文獻(xiàn)較少，已有的研究只是利用了個別省份內(nèi)部的數(shù)據(jù)，缺乏對國家區(qū)域內(nèi)的研究；二是在具體的實(shí)證研究中，忽視了工業(yè)TFP對環(huán)境污染空間溢出的影響。

本文關(guān)注中國工業(yè)TFP對環(huán)境污染的影響。一般來說，區(qū)域工業(yè)TFP的提高，可以通過區(qū)域工業(yè)的技術(shù)進(jìn)步、資源配置效率的改善等途徑實(shí)現(xiàn)，而無論工業(yè)生產(chǎn)中的技術(shù)進(jìn)步、資源配置效率的改善都會有利于生產(chǎn)中能源、原材料等物質(zhì)的投入減少，從而有利于降低環(huán)境污染物的排放。從研究方法方面，由于本文采用的是工業(yè)企業(yè)的分地區(qū)數(shù)據(jù)，這就不能忽略變量在地理空間上的溢出效應(yīng)?？臻g溢出效應(yīng)，是指一個變量發(fā)生變化時，不僅會對本區(qū)域目標(biāo)產(chǎn)生影響，而且會對區(qū)域之外目標(biāo)產(chǎn)生影響。而對于本文中的企業(yè) TFP 這一解釋變量來說，地區(qū)間企業(yè) TFP 水平差異會帶來競爭效應(yīng)、學(xué)習(xí)效應(yīng)和知識溢出效應(yīng)，使得其在空間上存在明顯的溢出效應(yīng)性。同樣，一些環(huán)境污染物能夠很容易的跨區(qū)域流動也呈現(xiàn)出非常明顯的空間相關(guān)性。如果我們在計量建模時忽略了這些特性，而僅將各個地區(qū)視作一個獨(dú)立的樣本，采用經(jīng)典的計量模型進(jìn)行分析，勢必會使得估計結(jié)果產(chǎn)生偏誤，也不能客觀地反映經(jīng)濟(jì)事物之間的空間聯(lián)系。

按照已有學(xué)者對空間計量經(jīng)濟(jì)學(xué)的研究，可以通過設(shè)置空間權(quán)重的模型來反映空間數(shù)據(jù)之間的相關(guān)關(guān)系。[28]本文亦將利用空間計量模型來考察中國工業(yè)TFP與環(huán)境污染的相關(guān)關(guān)系。

二、研究方法與數(shù)據(jù)來源

（一）空間相關(guān)性檢驗(yàn)?zāi)Ｐ?/p>

空間自相關(guān)分析方法由全局空間自相關(guān)和局部空間自相關(guān)兩部分構(gòu)成，本文用空間自相關(guān)指數(shù)Moran's I來檢驗(yàn)中國各地區(qū)的環(huán)境污染的空間相關(guān)性是否存在。Moran's I指數(shù)公式：

[I=ni=1nj=1nWij（Xi-X）（Xj-X）i=1nj=1nWiji=1n（Xi-X）2=i=1nj=1nWij（Xi-X）（Xj-X）S2i=1nj=1nWij]

[S2=1ni=1n（Xi-X）2]，[X=1ni=1nXi]

式中：[n]為空間單元數(shù)，[Wij]為空間權(quán)重矩陣，[Xi]和[Xj]分別是地區(qū)[i] 和[j] 的實(shí)際產(chǎn)出觀測值， [X]是觀測值的平均值，[S2]為離差平方和的均值。Moran's I的取值范圍為[-1 1]。若該指數(shù)小于0，則目標(biāo)區(qū)域存在空間負(fù)相關(guān)性時；若該指數(shù)大于則0，則目標(biāo)區(qū)域存在空間正相關(guān)性時，空間相關(guān)性越明顯則指數(shù)值越大；若此指數(shù)等于0，則目標(biāo)區(qū)域空間分布相互獨(dú)立。空間權(quán)重矩陣[Wij]采用地理距離權(quán)重矩陣進(jìn)行設(shè)定，用兩地區(qū)省會城市間的公路里程的算術(shù)平均值平方的倒數(shù)來構(gòu)造其元素。另外，區(qū)域間環(huán)境污染的空間集聚特征還可以通過繪制區(qū)域環(huán)境污染 Moran 散點(diǎn)圖來直觀反映。

（二）空間計量模型與溢出效應(yīng)分解

為了實(shí)證考察工業(yè)TFP對環(huán)境污染的空間影響，本文擬建立如下空間計量模型：

[Y=αLn+ρWY+βX+θWX+ε] （2）式（2）中，[Y]為環(huán)境污染排放，[X]為各工業(yè)TFP，同時，在實(shí)證研究過程中，加入了經(jīng)濟(jì)發(fā)展水平、人力資本存量、產(chǎn)業(yè)結(jié)構(gòu)及外商投資四項控制變量。 [α]為常數(shù)項，[Ln]為[n]×1階單位矩陣，[n] 為省份個數(shù)，[ε]為誤差項；[W] 為空間權(quán)重矩陣， [WY]和 [WX]分別考慮了被解釋變量和解釋變量的空間相關(guān)性和空間依賴性。在空間計量模型的估計結(jié)果中，如果[ρ]≠0，則[WY]的回歸系數(shù)[ρ]、 [WX]的回歸系數(shù)[θ]以及 [X]的 回 歸 系 數(shù) [β]就不能直接解釋變量的空間溢出效應(yīng)。為了對空間計量模型的回歸系數(shù)作出合理解釋，可以借鑒Lesage and Pace著作提出的的空間模型偏微分方法，可將（2）式改為：

[（In-ρW）Y=aLn+βX+θWX+ε]

[Y=r=1kSr（W）Xr+V（W）Lna+V（W）ε] （3）

[Sr（W）=V（W）（Inβr+Wθr）] （4）

[V（W）=（In-ρW）-1] （5）

其中，[X]為[n×k]矩陣，即樣本容量為[n]的[k]列解釋變量，[Xr]為第[r]個解釋變量；[In]表示[n]階單位矩陣，[θr]表示[WX]的第[r]個變量的估計系數(shù)。將（3）式展開可得：

[Y1Y2 ?Y3=Sr（W）11 Sr（W）12 … Sr（W）1nSr（W）21 Sr（W）22 … Sr（W）2n ? ? ? ?Sr（W）n1 Sr（W）n2 … Sr（W）nnX1rX2r ?Xnr] [+（In-ρW）-1+（ε+Lna）] （6）

[?Yi?Xjr=Sr（W）ij] [?Yi?Xir=Sr（W）ii]

其中，[Sr（W）ij]衡量的是區(qū)域[j]的第[r]個解釋變量對區(qū)域[i]被解釋變量的影響；[Sr（W）ii]衡量的是區(qū)域[i]的第[r]個解釋變量對區(qū)域[i]被解釋變量的影響。某個地區(qū)解釋變量的變化將不僅影響本地區(qū)的解釋變量，而且影響其他地區(qū)的被解釋變量。將前者稱為直接效應(yīng)，后者稱為間接效應(yīng)。

（三）變量設(shè)計

1.中國工業(yè)TFP

由于數(shù)據(jù)包絡(luò)分析法是一種非參數(shù)方法，無需設(shè)立特定的函數(shù)具體形式，只需通過測算產(chǎn)出與投入相對效率的前沿面，本文亦選用由Fare等提出的DEA-Malmquist指數(shù)方法來計算工業(yè)TFP，我國30個省級區(qū)域（西藏因數(shù)據(jù)少排除在基本決策單元之外）為基本決策單元（decision making unit， DMU），通過構(gòu)造每一時期的最佳實(shí)踐前沿來度量各個省區(qū)全要素生產(chǎn)率變化，要素投入（In-put）包含勞動力、資本，而總產(chǎn)出（Output）則把工業(yè)總產(chǎn)值作為指標(biāo)。

為了剔除價格變化的影響，用工業(yè)品出廠價格指數(shù)將各年的名義總產(chǎn)值轉(zhuǎn)化為以 2001年價格表示的實(shí)際總產(chǎn)值。勞動力指標(biāo)來自于工業(yè)企業(yè)全部從業(yè)人員年平均人數(shù)。資本存量用工業(yè)固定資產(chǎn)凈值年均余額表示，數(shù)據(jù)通過固定資產(chǎn)投資價格指數(shù)平減為以 2001年價格表示的數(shù)據(jù)。

2.環(huán)境污染（EP）

很多研究主要使用工業(yè)“三廢”指標(biāo)中的某個指標(biāo)或全部指標(biāo)來衡量環(huán)境污染程度，本文和白俊紅（2015）[23]一致選擇二氧化硫作為環(huán)境污染的代表指標(biāo)。其原因在于：一是在中國面臨的環(huán)境污染問題中，空氣污染最引人注目，其是大氣污染物的主要成分，而且對人的身體有直接危害性，當(dāng)作為一種重要化合物構(gòu)成酸雨后，對植物危害性也嚴(yán)重，因而是造成環(huán)境污染的一個重要因素；二是二氧化碳的分省數(shù)據(jù)相比其它數(shù)據(jù)來說較為完備且容易得到，方便進(jìn)行研究。具體來說，本文通過考察二氧化硫排放強(qiáng)度，即單位 GDP的工業(yè)二氧化硫排放量作為指標(biāo)進(jìn)行分析。

3.其它控制變量

本文還設(shè)計了其它控制變量以降低變量遺漏所造成的估計結(jié)果偏差。

（1）經(jīng)濟(jì)發(fā)展水平（EC）。環(huán)境問題始終是離不開經(jīng)濟(jì)發(fā)展水平。自英國工業(yè)革命以來，世界主要國家的經(jīng)濟(jì)發(fā)展歷程表明，經(jīng)濟(jì)發(fā)展初期總是伴隨著環(huán)境污染，當(dāng)前許多發(fā)展中國家依然以犧牲環(huán)境來換取較快的經(jīng)濟(jì)增長，而許多發(fā)達(dá)國家已經(jīng)越過環(huán)境庫茲尼茨曲線的拐點(diǎn)，其經(jīng)濟(jì)增長模式已對環(huán)境污染影響不大。本文采用地區(qū)GDP來反映地區(qū)的經(jīng)濟(jì)發(fā)展水平。

（2）人力資本（HC）。人力資本反映了凝聚在人體之中的知識、技能和熟練程度等。地區(qū)人力資本水平越高，意味著該地區(qū)人們的知識素質(zhì)和技能儲備水平也越高，從而也越有利于采用集約化的生產(chǎn)方式，進(jìn)而也有益于環(huán)境改善。人力資本水平數(shù)據(jù)采用6 歲及以上年齡人口的平均受教育年限，不同等級的教育年限設(shè)定為小學(xué)6 年、初中9年、高中12年、大專及以上16年，具體計算為總教育年限與總?cè)丝诘谋戎怠?/p>

（3）產(chǎn)業(yè)結(jié)構(gòu)（IS）。與第一和第三產(chǎn)業(yè)相比，環(huán)境污染物主要來源于第二產(chǎn)業(yè)的污染排放。自1992年來，我國第二產(chǎn)業(yè)占GDP 的比重一直保持在40%以上，因而第二產(chǎn)業(yè)占比成為影響環(huán)境污染的一個不容忽視的因素。本文用各地區(qū)第二產(chǎn)業(yè)產(chǎn)值/地區(qū)國內(nèi)生產(chǎn)總值來衡量。

（4）外商投資（FDI）： 以工業(yè)企業(yè)總產(chǎn)值中外商投資和港澳臺商投資占比。一種說法認(rèn)為，外商直接投資可以給東道國帶來先進(jìn)的生產(chǎn)技術(shù)和管理方式，有利于環(huán)境污染的減少，外資的產(chǎn)業(yè)發(fā)展和技術(shù)溢出效應(yīng)等也會促進(jìn)環(huán)境改善；還有些學(xué)者認(rèn)為了，外資會把落后的、污染嚴(yán)重產(chǎn)業(yè)帶到東道國，成為所謂的“污染天堂”。本文選擇實(shí)際外商直接投資利用額占GDP比重來衡量。

鑒于數(shù)據(jù)的可獲得性，本文選取 2001-2015年中國大陸30 個省份（西藏除外）的面板數(shù)據(jù)。為了確保數(shù)據(jù)的真實(shí)可靠，本文利用各個數(shù)據(jù)庫和年鑒對所選數(shù)據(jù)進(jìn)行了充分的核對、篩選；為了提高估計的準(zhǔn)確性，除了前文提到以外，本文用固定資產(chǎn)投資價格指數(shù)對相應(yīng)的貨幣量進(jìn)行平減指數(shù)計算調(diào)整成基期為 2001年的數(shù)值；為了減少異方差和統(tǒng)計偏誤，本文對相關(guān)變量進(jìn)行了自然對數(shù)處理。

四、實(shí)證分析

（一）省域TFP測度

本文采取2001-2015 年面板數(shù)據(jù)和DEA-Malmquist指數(shù)方法測出各地各期的工業(yè)TFP，表1是2001-2015年全國的TFP測度結(jié)果。

由表1分析得知，總體上，我國工業(yè)FP在2001-2004年期間呈下降趨勢，但2004年以后我國工業(yè)TFP逐漸上升，說明2004年以后，我國工業(yè)經(jīng)濟(jì)發(fā)展中的創(chuàng)新驅(qū)動力在增強(qiáng)。

（二）環(huán)境污染空間相關(guān)性檢驗(yàn)

根據(jù)所測算的各期省域環(huán)境污染，在地理距離權(quán)重矩陣下，采用Moran's I指數(shù)檢驗(yàn)了考察期內(nèi)我國各省環(huán)境污染的空間相關(guān)性，表2報告了檢驗(yàn)結(jié)果，圖1描述了其演變趨勢。

從表2和圖1結(jié)果可以看出，Moran's I指數(shù)均顯著為正，整體上一直趨于增長態(tài)勢，并且這種地理空間分布上的空間依賴趨勢自2008年以后表現(xiàn)的更為顯著。表明各區(qū)域的低環(huán)境污染分布并非是完全隨機(jī)的，表現(xiàn)出較強(qiáng)的空間集群特征，其空間聯(lián)系的特征是： 環(huán)境污染相似的地區(qū)趨于集聚，也就是說環(huán)境污染較高的省區(qū)相對趨于和環(huán)境污染較高的省區(qū)相鄰近；較低的省區(qū)相對趨于和較低的省區(qū)相鄰近。值得提出的是，在分析我國環(huán)境污染演化過程中，不能排除全局空間自相關(guān)性的存在。

區(qū)域環(huán)境污染的Moran's I指數(shù)已經(jīng)表明，區(qū)域環(huán)境污染的全局空間自相關(guān)性，用 Moran 散點(diǎn)圖可以進(jìn)一步說明區(qū)域環(huán)境污染分布的局部特征，如圖2 所示。受篇幅所限，本文僅以 2015年結(jié)果為例說明。

（三）空間計量模型設(shè)定

根據(jù)上文對各省空間環(huán)境污染相關(guān)性檢驗(yàn)，各空間單元存在集群現(xiàn)象，隨著時間推移，這種現(xiàn)象更為顯著。為了檢驗(yàn)中國工業(yè)TFP對環(huán)境污染的空間溢出效應(yīng)，需進(jìn)一步構(gòu)建空間計量模型進(jìn)行，采取合適方法進(jìn)行估計。首先建立非空間經(jīng)典計量模型來檢驗(yàn)工業(yè)TFP對區(qū)域環(huán)境污染的影響，然后再納入空間因素建立空間面板計量模型。經(jīng)豪斯曼檢驗(yàn) （Hausman test），經(jīng)典計量模型應(yīng)采用固定效應(yīng) （Fixed effects） 模型。

式（7）為區(qū)域環(huán)境污染非空間面板數(shù)據(jù)模型的估計結(jié)果：

其中，調(diào)整后的[R2]=0.9112，[F]=312.6920，[Log-L]=419.512 3， [D-W]=0.521 7。（7）式括號中的數(shù)值為顯著性概率。通過回歸結(jié)果，TFP、經(jīng)濟(jì)發(fā)展水平與產(chǎn)業(yè)結(jié)構(gòu)的回歸系數(shù)通過了 1% 的顯著性水平檢驗(yàn)，且前面兩者的系數(shù)均為負(fù)，這表明工業(yè)TFP當(dāng)前對環(huán)境污染水平有顯著的抑制作用，同時，表明近十幾年來我國的經(jīng)濟(jì)發(fā)展模式向環(huán)境友好型轉(zhuǎn)型；產(chǎn)業(yè)結(jié)構(gòu)的系數(shù)為正，說明我國第二產(chǎn)業(yè)的污染治理技術(shù)還是處于低端水平，而且第二產(chǎn)業(yè)中傳統(tǒng)重化工產(chǎn)業(yè)還是占有很大比重；而人力資本存量和外商投資系數(shù)雖然為負(fù)，但是并不顯著。

從非空間經(jīng)典計量模型的估計結(jié)果還可以看出，變量之間存在顯著的相關(guān)性 （[D-W]統(tǒng)計量值為0.5217）。由此，我們應(yīng)該建立區(qū)域環(huán)境污染的空間計量模型，以消除變量之間的空間相關(guān)性。我們經(jīng)常使用的空間計量模型主要包括空間滯后模型 （SLM） （ 也稱空間自回歸模型，SAR） 、空間誤差模型（SEM）和空間杜賓模型（ SDM） 三種。本文對這三種模型進(jìn)行固定效應(yīng)和隨機(jī)效應(yīng)檢驗(yàn)。

首先根據(jù) AIC 準(zhǔn)則，在固定效應(yīng)與隨機(jī)效應(yīng)之間進(jìn)行選擇，AIC 值越小則模型的解釋能力越強(qiáng)。其次，對比每一模型的Log likeli-hood與[R2]值，它們的值越大說明模型的擬合優(yōu)度越高。第三步，對空間杜賓模型是否可以轉(zhuǎn)化為空間滯后模型和空間誤差模型進(jìn)行檢驗(yàn)，如果檢驗(yàn)結(jié)果拒絕了可以轉(zhuǎn)化的原假設(shè) （兩個原假設(shè)為[H0]：[θ=0] 和[H0]：[θ+δρ=0] ），則選擇杜賓模型。按照上述步驟，分別對 SLM、SDM、SEM 進(jìn)行了固定效應(yīng)和隨機(jī)效應(yīng)回歸，表3報告了回歸結(jié)果。根據(jù)遴選原則，文中選擇了杜賓模型的固定效應(yīng)模型，并將其作為最終的解釋模型，進(jìn)行空間溢出效應(yīng)分解和分析。

從表3的估計結(jié)果，系數(shù)[ρ]值顯著不為零，且回歸系數(shù)均保持在0.5左右，表明不能直接用5個解釋變量的回歸系數(shù)來解釋各自的經(jīng)濟(jì)意義。因此，本文采用空間回歸模型偏微分方法，將TFP及其他四個控制變量對環(huán)境污染的空間溢出效應(yīng)進(jìn)行分解，表4報告了分解結(jié)果。

1.工業(yè)TFP的空間溢出效應(yīng)。實(shí)際結(jié)果顯示，工業(yè)TFP對環(huán)境污染的降低存在溢出效應(yīng)，且直接效應(yīng)、間接效應(yīng)和總效應(yīng)的系數(shù)分別通過1%、10%和5%的顯著性檢驗(yàn)。工業(yè)TFP每提高1%，直接促進(jìn)區(qū)域環(huán)境污染減少0.951%，間接減少環(huán)境污染0.127 1%，總體上減少環(huán)境污染1.078 1%，表明工業(yè)TFP的提高能夠有效降低環(huán)境污染。

2.其他控制變量的溢出效應(yīng)。根據(jù)表4的分解結(jié)果，經(jīng)濟(jì)發(fā)展水平每提高1%，直接減少環(huán)境污染0.167 1%，且通過了5%的顯著性檢驗(yàn)，間接增加0.073 2%，且并不顯著，總體上減少環(huán)境污染0.093 9%，通過了10%的顯著性檢驗(yàn)。人力資本存量每提高1%，直接減少環(huán)境污染0.132 0%，間接增加環(huán)境污染0.073 2%，總體上減少環(huán)境污染0.067%，但是都不顯著。產(chǎn)業(yè)結(jié)構(gòu)變量每提高1%，直接增加環(huán)境污染0.291 1%，間接增加環(huán)境污染0.012 3%，總體增加環(huán)境污染0.303 4%，且分別通過1%、10%和5%的顯著性檢驗(yàn)。外商投資水平每增加1%，直接減少環(huán)境污染0.113 5%，且通過10%的顯著性檢驗(yàn)，間接減少0.042 1%，總體上減少環(huán)境污染0.155 6%，但是間接效應(yīng)和總體效應(yīng)都不顯著。

五、結(jié)論

在當(dāng)前我國經(jīng)濟(jì)增長方式從要素驅(qū)動向創(chuàng)新驅(qū)動轉(zhuǎn)變、加強(qiáng)生態(tài)文明建設(shè)的背景下，研究工業(yè)全要素增長率（TFP）與環(huán)境污染的關(guān)系具有重要的理論價值和現(xiàn)實(shí)意義。本文利用2001-2015年的省級面板數(shù)據(jù)，通過構(gòu)建空間動態(tài)面板模型，實(shí)證測算了工業(yè)TFP以及經(jīng)濟(jì)發(fā)展水平、人力資本、產(chǎn)業(yè)結(jié)構(gòu)、外商投資等影響因素對環(huán)境污染的直接效應(yīng)、間接效應(yīng)以及總效應(yīng)。研究結(jié)論如下：

（ 1） 各省級區(qū)域環(huán)境污染具有顯著的空間相關(guān)性，即各省級區(qū)域環(huán)境污染程度不僅受到本區(qū)域工業(yè)TFP等因素的影響，同時受到來自鄰近省級區(qū)域的工業(yè)TFP等因素的空間影響，區(qū)域間環(huán)境污染之間存在一定的自相關(guān)性；（2） 在考慮各省級區(qū)域空間相關(guān)性前提下，并控制了其他解釋變量的影響作用以后，本區(qū)域工業(yè)TFP每提高1% ，能夠顯著的促進(jìn)本區(qū)域環(huán)境污染程度減少 0.951 0%，并帶動區(qū)域間環(huán)境污染程度減少0.127 1% ，表明工業(yè)TFP對環(huán)境污染的減少具有顯著的直接溢出效應(yīng)和間接溢出效應(yīng)；（3）本區(qū)域的經(jīng)濟(jì)發(fā)展水平提高能夠促進(jìn)本區(qū)域的環(huán)境污染程度降低，對環(huán)境污染程度降低的總效應(yīng)顯著，說明當(dāng)前經(jīng)濟(jì)發(fā)展向環(huán)境友好型發(fā)展，間接效應(yīng)雖不顯著但系數(shù)為正，其原因是區(qū)域經(jīng)濟(jì)發(fā)展，使當(dāng)?shù)氐漠a(chǎn)業(yè)結(jié)構(gòu)進(jìn)一步優(yōu)化，而將一些污染嚴(yán)重的產(chǎn)業(yè)轉(zhuǎn)移到鄰近地區(qū)，而造成鄰近地區(qū)的環(huán)境污染加重了；（4）產(chǎn)業(yè)結(jié)構(gòu)水平對本地區(qū)和相鄰地區(qū)環(huán)境污染程度增加具有明顯的促進(jìn)作用，說明在傳統(tǒng)重化工產(chǎn)業(yè)還是占有很大比重的第二產(chǎn)業(yè)里，其污染治理技術(shù)還需進(jìn)一步提高；（5）外商投資水平對環(huán)境污染程度的直接效應(yīng)系數(shù)為負(fù)且顯著，說明外資企業(yè)的增加帶來了一定的先進(jìn)環(huán)境治理技術(shù)和環(huán)境理念，很多外資企業(yè)是出口導(dǎo)向型的，為了達(dá)到出口較嚴(yán)格的環(huán)境標(biāo)準(zhǔn)，也迫使產(chǎn)業(yè)鏈相關(guān)企業(yè)提高了環(huán)境質(zhì)量。

本文主要探討了工業(yè)TFP對區(qū)域環(huán)境污染的空間溢出影響，實(shí)證結(jié)果說明一個地區(qū)的工業(yè)TFP提高能夠有效通過技術(shù)進(jìn)步或配置效率的提高降低本地的環(huán)境污染，而且還能帶動周邊地區(qū)的技術(shù)進(jìn)步或配置效率提高，從而一定程度上降低周邊地區(qū)的環(huán)境污染。同時，也說明國家大力促進(jìn)創(chuàng)新驅(qū)動，改變經(jīng)濟(jì)增長方式，未來將會對區(qū)域環(huán)境生態(tài)有明顯的改善作用。本研究通過實(shí)證得到了一些有益的結(jié)論，但也存在一些不足，例如我國區(qū)域經(jīng)濟(jì)發(fā)展水平存在差異，工業(yè)TFP對環(huán)境污染的影響也可能存在地區(qū)差異，文章未將考察對象分地區(qū)研究，稍顯不足。作者將在后續(xù)研究中進(jìn)一步探討這些問題。

參考文獻(xiàn)：

[1] 楊汝岱.中國制造業(yè)企業(yè)全要素生產(chǎn)率研究[J].經(jīng)濟(jì)研究，2015（2）：61-74.

[2] 張軍， 施少華.中國經(jīng)濟(jì)全要素生產(chǎn)率變動：1952—1998[J]. 世界經(jīng)濟(jì)文匯， 2003（2）： 17-24.

[3] 郭慶旺， 賈俊雪.中國全要素生產(chǎn)率的估算： 1979—2004[J]. 經(jīng)濟(jì)研究， 2005（5）： 51-60.

[4] 傅曉霞， 吳利學(xué). 技術(shù)效率， 資本深化與地區(qū)差異——基于隨機(jī)前沿模型的中國地區(qū)收斂分析[J]. 經(jīng)濟(jì)研究， 2006（10）： 52-61.

[5] 李靜， 孟令杰， 吳福象. 中國地區(qū)發(fā)展差異的再檢驗(yàn)： 要素積累抑或 TFP[J]. 世界經(jīng)濟(jì)， 2006， 29（1）： 12-22.

[6] 李征. 中國全要素生產(chǎn)率演變研究[D]. 吉林大學(xué)，2016.

[7] HUANG Y， MENG X. China's Industrial Growth and Efficiency： A comparison between the state and the TVE sectors [J]. Journal of the Asia Pacific Economy， 1997， 2（1）： 101-121.

[8] COELLI T J， PRASADA RAO D S. Total Factor Productivity Growth in Agriculture： A Malmquist Index Analysis of 93 Countries， School of Economics， University of Queensland[R]. Working Paper Series， 2003.

[9] 朱鐘棣， 李小平. 中國工業(yè)行業(yè)資本形成、全要素生產(chǎn)率變動及其趨異化： 基于分行業(yè)面板數(shù)據(jù)的研究[J]. 世界經(jīng)濟(jì)， 2005， 28（9）： 51-62.

[10] 涂正革. 全要素生產(chǎn)率與區(qū)域經(jīng)濟(jì)增長的動力——基于對 1995—2004 年 28 個省市大中型工業(yè)的非參數(shù)生產(chǎn)前沿分析[J]. 南開經(jīng)濟(jì)研究， 2007（4）： 14-36.

[11] 楊勇. 中國服務(wù)業(yè)全要素生產(chǎn)率再測算[J]. 世界經(jīng)濟(jì)， 2008（10）： 46-55.

[12] 劉興凱， 張誠. 中國服務(wù)業(yè)全要素生產(chǎn)率增長及其收斂分析[J]. 數(shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究， 2010（3）： 55-67.

[13] MILLER S M，UPADHTAY M P. The effects of openness， trade orientation， and human capital on total factor productivity[J].Journal of development economics， 2000， 63（2）： 399-423

[14] 顏鵬飛， 王兵. 技術(shù)效率， 技術(shù)進(jìn)步與生產(chǎn)率增長： 基于 DEA 的實(shí)證分析[J]. 經(jīng)濟(jì)研究， 2004（9）： 88-96.

[15] 趙偉， 張萃. 中國制造業(yè)區(qū)域集聚與全要素生產(chǎn)率增長[J]. 上海交通大學(xué)學(xué)報（哲學(xué)社會科學(xué)版）， 2008， 16（5）： 52-56.

[16] 劉秉鐮， 武鵬， 劉玉海. 交通基礎(chǔ)設(shè)施與中國全要素生產(chǎn)率增長——基于省域數(shù)據(jù)的空間面板計量分析[J]. 中國工業(yè)經(jīng)濟(jì)， 2010（3）： 54-64.

[17] 蔣殿春， 王曉嬈. 中國 R&D 結(jié)構(gòu)對生產(chǎn)率影響的比較分析[J]. 南開經(jīng)濟(jì)研究， 2015 （2）： 59-73.

[18] 李國璋，孔令寬.環(huán)境庫茨涅茨曲線在中國的適用性[ J].廣東社會科學(xué)，2008（2）：37- 43.

[19] 彭水軍，包群.中國經(jīng)濟(jì)增長與環(huán)境污染[ J].世界經(jīng)濟(jì)，2006（5）：48- 58.

[20] 曹光輝，汪鋒，張宗益，鄒暢.我國經(jīng)濟(jì)增長與環(huán)境污染關(guān)系研究[J].中國人口、資源與環(huán)境，2006（1）：25- 29.

[21] 黃菁.環(huán)境污染與工業(yè)結(jié)構(gòu)[ J].統(tǒng)計研究，2009（12）：68-73.

[22] 鄧柏盛，宋德勇.我國對外貿(mào)易、FDI與環(huán)境污染之間關(guān)系的研究[J].國際貿(mào)易問題，2008（4）：101- 108.

[23] 白俊紅，呂曉紅. FDI質(zhì)量與中國環(huán)境污染的改善[J].國際貿(mào)易問題，2015（8）：72-83.

[24] 楊俊，邵漢華.環(huán)境約束下的中國工業(yè)增長狀況研究——基于Malmquist-Luenberger指數(shù)的實(shí)證分析[J].數(shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究，2009（9）：64-78.

[25] 吳軍，笪鳳媛，張建華.環(huán)境管制與中國區(qū)域生產(chǎn)率增長[J].統(tǒng)計研究，2010（1）：83-89.

[26] 鄭寶華、謝忠秋.基于低碳經(jīng)濟(jì)的中國區(qū)域全要素生產(chǎn)率研究[J].經(jīng)濟(jì)學(xué)動態(tài)，2011（10）：38-41.

[27] 聶國卿，尹向飛. 湖南省工業(yè)全要素生產(chǎn)率對環(huán)境污染的影響研究 [J].湖南大學(xué)學(xué)報，2010（11）：61-66.

[28] 李婧，譚清美，白俊紅.中國區(qū)域創(chuàng)新生產(chǎn)的空間計量分析——基于靜態(tài)與動態(tài)空間面板模型的實(shí)證研究[J]. 管理世界，2010（7）：43-55.

Industry TFP， Spatial Spillover and Environmental Pollution in China

——Based on Spatial Measurement on Provincial Unit

ZHENG Bao-hua，LIU Dong-huang

（School of Business，Jiangsu University of Technology，Changzhou 213001， China）

Abstract： On the basis of provincial panel data during 2001-2015， the paper calculates industrial TFP based on Malmquist method， constructs a spatial econometric model， analyses spatial spillover effect of China's industrial TFP on regional environment pollution. The study finds that the provincial environmental pollution has significant positive spatial correlation between the regions； The industrial TFP has a significant role in the regional environmental pollution as well as outside region environmental pollution， the spillover effect is significant， Innovation driven and the change of growth mode has obvious effect on improving the regional ecological environment. The economic development level have a significant role in promoting environmental improvement in the region； Second industry accounts for a significant role in promoting environmental pollution in the region and adjacent areas； Foreign investment plays a significant role in promoting the regional environmental improvement.

Key words： industry TFP；spatial spillover； environmental pollution

責(zé)任編輯 趙文清