(華僑大學 經(jīng)濟與金融學院,福建 泉州 362021)

工業(yè)環(huán)境效率是工業(yè)經(jīng)濟可持續(xù)發(fā)展的關(guān)鍵,也是新時期工業(yè)產(chǎn)業(yè)轉(zhuǎn)型和環(huán)境保護面臨的重大挑戰(zhàn)。一方面，對外貿(mào)易、FDI作為技術(shù)來源的重要途徑,是提高工業(yè)環(huán)境效率的關(guān)鍵動力之一;另一方面,根據(jù)“向底線賽跑”論斷和“污染天堂”假說,對外貿(mào)易、FDI會使發(fā)達國家或地區(qū)的污染性企業(yè)遷移到落后國家或地區(qū),進而造成污染、降低環(huán)境效率。就目前的發(fā)展看,整體上我國仍是發(fā)展中國家；局部上,我國地域遼闊,各地區(qū)經(jīng)濟發(fā)展差異較大。從東至西,內(nèi)部區(qū)域可分為相對發(fā)達地區(qū)、中等發(fā)達地區(qū)和相對落后地區(qū),空間分布等級明顯,且不同地區(qū)的區(qū)位優(yōu)勢導致貿(mào)易活動在各地區(qū)存在明顯差異。因此,本文從對外貿(mào)易和FDI角度出發(fā)探討對我國工業(yè)環(huán)境效率影響的空間差異具有重要的理論價值和實際意義。

對外貿(mào)易對環(huán)境效率影響研究最先是由Grossman、Krueger建立相應的分析框架,并將國際貿(mào)易環(huán)境影響分解為規(guī)模效應、結(jié)構(gòu)效應和技術(shù)效應3個方面[1]。此后,國內(nèi)外許多學者開展了相關(guān)研究,但結(jié)果不盡相同。其中,國外Dua、Esty的觀點認為,對外貿(mào)易的全球化競爭會使各國以環(huán)境為代價,降低環(huán)境效率[2]。另一類研究表明對外貿(mào)易有利于環(huán)境效率的提高[3,4]。在國內(nèi),對外貿(mào)易與環(huán)境污染、對外貿(mào)易與能源效率以及對外貿(mào)易與工業(yè)環(huán)境效率都是學者討論的熱點,相關(guān)研究也存在較大的爭議。一類研究認為，對外貿(mào)易有利于能源環(huán)境效率或工業(yè)環(huán)境效率的提高,有利于環(huán)境保護[5-7];另一類研究則截然相反,認為發(fā)達國家的貿(mào)易轉(zhuǎn)移使我國成為環(huán)境污染的“避難所”,對外貿(mào)易惡化了我國的環(huán)境[8]。對FDI與環(huán)境效率的研究以“污染天堂”假說為代表,該觀點在國內(nèi)外研究均存在較大爭議。在我國,相關(guān)研究可以分為三類：第一類認為FDI對環(huán)境效率具有顯著的促進作用,否定了“污染天堂”假說[9-11];第二類認為FDI對我國環(huán)境效率的作用不顯著,原因在于外商投資一方面給本地工業(yè)企業(yè)帶來了先進技術(shù),促進了科技合作,另一方面又引進了發(fā)達國家高耗能、高污染的產(chǎn)業(yè),兩方面的共同作用給工業(yè)環(huán)境效率帶來正、負兩方面影響,兩種影響相互抵消[12];第三類認為我國存在著“污染天堂”假說的證據(jù)[13]。

從已有的研究來看,對外貿(mào)易和FDI均可作為市場開放程度的相關(guān)指標,可認為是技術(shù)進步的重要來源,但是研究結(jié)果大相徑庭。筆者認為,其中一個重要的原因是由于技術(shù)手段的限制,這些研究未能很好地考慮空間因素,從空間異質(zhì)性的角度出發(fā)對該問題進行考量。目前,隨著相關(guān)研究和技術(shù)手段的不斷發(fā)展,空間探索性分析、時空地理加權(quán)回歸模型、面板分位回歸等方法為本研究提供了良好的研究基礎(chǔ),因此可從空間異質(zhì)性出發(fā),繼續(xù)考察這一問題,其研究結(jié)果將更具有說服力。在“一帶一路”倡議的發(fā)展背景下,我國對外貿(mào)易和外商直接投資格局正在發(fā)生變化,關(guān)于對外貿(mào)易和FDI對工業(yè)環(huán)境效率影響的探究對相關(guān)環(huán)境政策的制定提供參考與借鑒。

1 研究方法及數(shù)據(jù)來源

1.1 GTWR模型

傳統(tǒng)的回歸模型在估計區(qū)域數(shù)據(jù)時往往忽略了數(shù)據(jù)的空間結(jié)構(gòu),僅能得到平均影響,Fotheringham提出的地理加權(quán)回歸(GWR)模型考慮了空間的非平穩(wěn)性,解釋變量的影響隨著區(qū)域的不同而異,是變量研究空間異質(zhì)性的重要工具。在實際運用中,GWR模型也存在一定的缺陷。如該模型只能研究截面數(shù)據(jù),缺乏時間維度考慮,即使在研究中可研究不同時間節(jié)點上的情況來說明問題,但相較于面板數(shù)據(jù)節(jié)點數(shù)據(jù)具有一定的偶然性,因此仍存在一定的偏差。Huang等[14]提出的時空地理加權(quán)(GTWR)模型考慮了時間維度的變化,克服了地理加權(quán)回歸模型的不足。本文將采用GTWR模型來解決空間異質(zhì)性問題,模型表達式為:

(1)

式中,(yi;xi1,xi2,…,xid)為被解釋變量和解釋變量在時空坐標(ui,vi,ti)處的觀測值。通常,對i、j的空間距離可采用各個省域的經(jīng)度和緯度數(shù)據(jù)來計算。t為時間距離;εi～N(0,σ2)、βk(ui,vi,ti)為系數(shù)函數(shù),是函數(shù)βk(u,v,t)在i點的值。

利用時空加權(quán)回歸模型的局部線性估計方法可求得各回歸系數(shù)在觀測點i處的估計值:

=[XTW(ui,vi,ti)]-1XTW(ui,vi,ti)Y

(2)

權(quán)重W(ui,vi,ti)設(shè)為觀測點i到其他觀測點的距離的函數(shù),通常采用高斯距離函數(shù)。此外,帶寬選擇對GTWR模型的精度影響很大,常采用交叉確認法確認模型的空間帶寬和時間帶寬。在實際運用中,可按照AIC準則,即AIC最小時對應的為最優(yōu)帶寬。分位回歸模型的估計結(jié)果可反映不同層次中參數(shù)的變化狀況,因此可用來研究異質(zhì)性問題[15],且在估計中是根據(jù)樣本信息按不同分位點分割數(shù)據(jù)進行估計,有利于歸類分析。本文將在GTWR模型的基礎(chǔ)上繼續(xù)使用面板分位回歸方法對空間異質(zhì)性問題進行研究。

1.2 面板分位回歸

分位回歸模型最早由RogerKoenker、Gilbert Bassett于1987年提出,它克服了最小二乘法的一些局限,能很好地反應各分位點的影響。但一直以來分位回歸模型停留在對截面數(shù)據(jù)的估計,未能發(fā)揮面板數(shù)據(jù)的優(yōu)點,直到2004年Koenker提出固定效應面板分位回歸模型。目前面板分位回歸模型仍處于發(fā)展階段,該模型可充分發(fā)揮面板數(shù)據(jù)模型和分位數(shù)回歸優(yōu)點,在控制個體差異的基礎(chǔ)上，可有效地分析解釋變量對不同分位點上被解釋變量條件分布的影響程度。相較于時空地理加權(quán)回歸模型,面板分位回歸模型依據(jù)某一變量的統(tǒng)計數(shù)字特征對研究對象進行了歸類,能說明處于不同數(shù)字特征類別對象的空間差異,具有鮮明的層次性,可作為空間異質(zhì)性研究工具?，F(xiàn)有的面板分位回歸方法主要有三大類:一階差分分位回歸法、固定效應變換分位回歸法和引進虛擬變量的懲罰法[16],本文在研究中只考察研究對象的層次差異,以探討我國工業(yè)環(huán)境效率的空間異質(zhì)性問題。以固定效應分位回歸模型為分析方法,具體表達式為:

(3)

1.3 數(shù)據(jù)來源

本文研究數(shù)據(jù)來源于相關(guān)年份的《中國統(tǒng)計年鑒》、《中國對外貿(mào)易統(tǒng)計年鑒》、《中國環(huán)境統(tǒng)計年鑒》以及各省的統(tǒng)計年鑒。

2 對外貿(mào)易、FDI與工業(yè)環(huán)境效率的空間格局分析

2.1 對外貿(mào)易空間分布格局分析

利用各省貿(mào)易進出口總額數(shù)據(jù)對我國對外貿(mào)易空間格局進行分析?？紤]到一方面進出口貿(mào)易總額中包含價格因素的影響，另一方面各省進出口貿(mào)易總額在逐年增加,絕對數(shù)很難反應出各省的貿(mào)易實力,因此在分析中采用相對數(shù)當年各省進出口貿(mào)易總額占當年全國貿(mào)易總額的比值,以解決上述問題。運用ArcGIS軟件繪制出2000年、2008年和2015年對外貿(mào)易情況的空間分布圖(見圖1—4)。

圖1 2000年對外貿(mào)易空間分布情況

分析圖1—4可見：首先,我國對外貿(mào)易省際差距較大,呈現(xiàn)出顯著的空間非均衡發(fā)展。我國對外貿(mào)易東中西層級分布明顯,貿(mào)易活動主要集中在東部地區(qū),尤其是廣東、江蘇、上海、浙江和山東5個省市,西部地區(qū)的貿(mào)易強度較弱,尤其是甘肅、青海、寧夏和西藏4個省區(qū)。其次,從時間演變過程看,整個貿(mào)易活動呈現(xiàn)出了“強者愈強,弱者愈弱”的局面,整個貿(mào)易活動在逐步向東部沿海地區(qū)集聚。圖中呈現(xiàn)了各省進出口總額占全國進出口總額比重的空間分布情況,在時間維度上發(fā)現(xiàn)西部各省進出口總額所占比重明顯減少,中部各省也出現(xiàn)了不同程度的縮減。第三,東西區(qū)域內(nèi)存在差異顯著的省份。西部地區(qū)整體上貿(mào)易活動強度較弱,但四川省是在整個西部地區(qū)中經(jīng)濟發(fā)展較好的省份,新疆為我國向西開放的“橋頭堡”和“樞紐站”,貿(mào)易活動顯著強于西部其他省區(qū)。整體上,東部地區(qū)貿(mào)易活動較集中,但福建省的貿(mào)易強度和貿(mào)易演變過程較緩慢,與鄰近沿海省市差距較大。

圖2 2008年對外貿(mào)易空間分布情況

圖3 2015年對外貿(mào)易空間分布情況

圖4 2000年FDI空間分布情況

2.2 FDI空間格局分析

本文運用外商直接投資實際利用額數(shù)據(jù)對FDI空間分布格局進行分析研究，同時采用相對數(shù)對各省當年FDI實際利用額占全國當年FDI實際利用額的比值進行了分析,空間分布結(jié)果見圖5、圖6。從圖5、圖6可見:①2000年、2008年和2015年我國FDI空間非均衡性顯著,整體上東部地區(qū)各省市FDI水平顯著高于中部地區(qū)和西部地區(qū),區(qū)域?qū)蛹壊町愝^大,但也不乏個別“異?！钡氖∈?如2000年、2008年和2015年上海市的實際利用外資額分別為31.6億美元、100.84億美元和184.59億美元,占全國實際利用外商投資額的比重甚微,西部地區(qū)的四川省利用外商直接投資額較大,明顯強于西部其他省份。②整體上來說，各省、市、自治區(qū)實際利用外商投資額在逐漸增加,且呈現(xiàn)出自東向西逐漸擴展的發(fā)展趨勢,中部地區(qū)的FDI水平變動較顯著,成為外商直接投資的重點區(qū)域,尤其是河南、安徽、江西、湖南4個省。

圖5 2008年FDI空間分布情況

圖6 2015年FDI空間分布情況

2.3 工業(yè)環(huán)境效率及其空間格局分析

環(huán)境效率是提高環(huán)境質(zhì)量的關(guān)鍵所在,工業(yè)環(huán)境效率又是提高整體環(huán)境效率的重中之重。對我國工業(yè)環(huán)境效率的測度和空間分布格局的分析一方面是對我國工業(yè)環(huán)境現(xiàn)狀的具體認識,另一方面可與對外貿(mào)易、FDI空間分布格局進行對比,初步判斷三者間存在的相關(guān)關(guān)系。①建立工業(yè)環(huán)境效率測算指標體系。在工業(yè)環(huán)境效率的測算中,相關(guān)指標體系已較健全,在投入和產(chǎn)出變量的選擇中,成本即為投入變量,收益為產(chǎn)出變量,環(huán)境污染物可考慮為成本,因此可作為投入處理。借鑒相關(guān)研究[11,17]可建立的指標體系見表1。其中,根據(jù)現(xiàn)有相關(guān)文獻研究,工業(yè)能源消費量選取煤炭消費量為具體指標;工業(yè)資本投入采用永續(xù)盤存法測算,借鑒焦兵等[10]的研究,把規(guī)模以上工業(yè)企業(yè)固定資產(chǎn)凈值年均余額之差當作新增固定資產(chǎn)投資,參照張軍等[18]的成果,折舊率取9.6%,并且用固定資產(chǎn)投資價格指數(shù)進行平減;工業(yè)總產(chǎn)值以2000年為基期,運用各年工業(yè)出廠價格指數(shù)進行平減,以剔除價格因素的影響。②運用上述超效率DEA模型測算我國工業(yè)環(huán)境效率。超效率模型與傳統(tǒng)DEA模型具有相同的表達式,不同的是在計算第n個決策單元效率時,將第n個決策單元的投入和產(chǎn)出用所有單元決策單元投入和產(chǎn)出的線性組合代替。相較于傳統(tǒng)的DEA模型,超效率DEA模型的優(yōu)點在于不但能判斷決策單元是否有效,而且還可進行決策單元間的效率比較,根據(jù)效率計算結(jié)果對決策單元進行排名。結(jié)果中若效率值小于1,其意義同傳統(tǒng)DEA模型計算結(jié)果,表示相對無效;若測算結(jié)果為有效,則其效率值大于等于1。當值大于1時,表明在1的基礎(chǔ)上再擴大若干,決策單元仍能保持相對最優(yōu)效率。如若決策單元效率為1.6,則表示決策單元規(guī)模再擴大60%仍可保持相對最優(yōu)效率[19]。

表1 工業(yè)環(huán)境效率測算指標體系

根據(jù)以上測算結(jié)果,可以對我國工業(yè)環(huán)境效率的空間分布及演化過程做進一步的分析,見圖7—9。其中,由于香港特別行政區(qū)、澳門特別行政區(qū)、臺灣地區(qū)、海南省、西藏自治區(qū)部分關(guān)鍵數(shù)據(jù)缺乏,因此在測算中予以了剔除,但是為了完整性仍然保留了底圖。

圖7 2000年我國工業(yè)環(huán)境效率空間分布狀況

圖8 2008年我國工業(yè)環(huán)境效率空間分布狀況

圖9 2015年我國工業(yè)環(huán)境效率空間分布狀況

從圖7—9可見:①我國整體工業(yè)環(huán)境效率東部地區(qū)明顯高于中西部地區(qū),尤其是北京、上海、浙江、福建等省市。由于經(jīng)濟相對發(fā)達、技術(shù)條件優(yōu)越、環(huán)境治理投資充裕,工業(yè)效率處于明顯的領(lǐng)先地位。②從整體演變過程看,東部地區(qū)有進一步提高環(huán)境效率的趨勢,中西部地區(qū)與東部地區(qū)的環(huán)境效率差距并沒有隨時間的推移而逐漸縮小,而是呈現(xiàn)出繼續(xù)擴大的發(fā)展態(tài)勢。③黑龍江省屬于我國的老工業(yè)基地,經(jīng)濟發(fā)展相對落后,但工業(yè)環(huán)境效率一直處于較高水平,該結(jié)果與宋馬林等[6]的研究結(jié)果一致。西部地區(qū)中,廣西、重慶和新疆工業(yè)環(huán)境效率相對較高。④存在著環(huán)境效率變化顯著的省份。西部地區(qū)的四川省環(huán)境效率有了顯著提高,中部地區(qū)的河南省工業(yè)環(huán)境效率提高相對顯著,東部地區(qū)的天津市、福建省、廣東省、山東省的環(huán)境效率則出現(xiàn)下降,尤其是天津和福建兩地。

綜合上述研究發(fā)現(xiàn)：整體分布上，我國對外貿(mào)易、FDI與工業(yè)環(huán)境效率有基本一致的狀態(tài)。在演變路徑上,我國工業(yè)環(huán)境效率空間分布與對外貿(mào)易空間分布呈現(xiàn)出相同的變化趨勢,而與FDI的變化路徑相反。在部分省、市、自治區(qū)中,如四川的工業(yè)環(huán)境效率變化與FDI呈現(xiàn)出了一致的變化趨勢?？梢?，我國對外貿(mào)易和FDI與工業(yè)環(huán)境效率間存在著一定的相關(guān)關(guān)系。根據(jù)已有研究,對對外貿(mào)易、FDI和工業(yè)環(huán)境效率間的關(guān)系進行猜想:①對外貿(mào)易作為技術(shù)來源的關(guān)鍵因素,對我國工業(yè)環(huán)境效率存在著顯著的促進作用,貿(mào)易活動愈活躍的地區(qū),對外貿(mào)易對工業(yè)環(huán)境的促進作用愈明顯;在貿(mào)易活動貧乏的地區(qū),對外貿(mào)易對工業(yè)環(huán)境的促進作用不顯著。②在經(jīng)濟發(fā)達地區(qū),由于自身發(fā)展良好,會對進入的外資企業(yè)篩選,阻止具有污染性或污染嚴重投資,因此FDI對工業(yè)環(huán)境效率具有一定的促進作用。在投資環(huán)境優(yōu)越的地區(qū),FDI是主要技術(shù)來源,可促進地區(qū)工業(yè)環(huán)境效率的提高;在經(jīng)濟相對落后的地區(qū),FDI技術(shù)效應和環(huán)境成本效應相互抵消,結(jié)果正負均可能存在。通過GTWR模型和分位回歸模型對該判斷進行驗證,并對其他相關(guān)因素對我國工業(yè)環(huán)境效率影響的空間異質(zhì)性做一分析。

3 模型分析

3.1 控制變量

本文利用超效率DEA模型對工業(yè)環(huán)境效率進行了測度，在此將探討對外貿(mào)易與FDI對工業(yè)環(huán)境效率的影響和表現(xiàn)出的區(qū)域差異。將工業(yè)環(huán)境效率作為被解釋變量,對外貿(mào)易和FDI作為主要解釋變量。其中,對外貿(mào)易用各地區(qū)進出口貿(mào)易總額與地區(qū)生產(chǎn)總值的比值為具體指標,FDI用實際利用外商直接投資額為具體指標。根據(jù)相關(guān)研究,模型中加入以下控制變量:①市場化程度(MS)。市場化程度既會降低交易成本,也會刺激企業(yè)創(chuàng)新,有利于工業(yè)環(huán)境效率的提高[20]。具體采用工業(yè)私營企業(yè)銷售產(chǎn)值占規(guī)模以上工業(yè)銷售總產(chǎn)值的比重度量。②技術(shù)投入(R&D)。已有研究表明,技術(shù)進步會降低能源消耗,提高能源利用率,有利于工業(yè)環(huán)境效率的提高[21],而R&D是技術(shù)進步的主要推動力。技術(shù)投入的具體指標選擇規(guī)模以上工業(yè)企業(yè)R&D內(nèi)部經(jīng)費支出,部分年份只統(tǒng)計了大中型工業(yè)企業(yè)R&D數(shù)據(jù),因此部分指標用大中型工業(yè)企業(yè)數(shù)據(jù)進行替代。③工業(yè)化水平(IZ)。在環(huán)境效率較低的地區(qū),工業(yè)化水平拉動經(jīng)濟增長的正面效應大于負面效應,因此可能會出現(xiàn)工業(yè)化水平對環(huán)境效率的正向影響,環(huán)境效率較高的地區(qū),工業(yè)化水平可能會對環(huán)境效率產(chǎn)生負面影響[6],具體的工業(yè)化水平指標用地區(qū)工業(yè)增加值占GDP的比重表示,各變量的描述性統(tǒng)計見表2。

表2 模型變量描述性統(tǒng)計

注:數(shù)據(jù)中不包括香港與澳門特別行政區(qū)、臺灣地區(qū)、海南省和西藏自治區(qū)。

3.2 GTWR模型分析

根據(jù)上述模型理論和變量,估計具體模型,在估計中核函數(shù)選取高斯核函數(shù),帶寬依據(jù)AIC最小原則進行了選取,估計參數(shù)的描述性統(tǒng)計見表3。

表3 GTWR模型參數(shù)描述性統(tǒng)計

從表3可見,各變量對工業(yè)環(huán)境影響的區(qū)域差異不可忽視。首先主要解釋變量對外貿(mào)易、FDI。整體對外貿(mào)易和FDI對我國工業(yè)環(huán)境效率的影響都呈現(xiàn)出一定的梯度分布,其中FDI的梯度分布更明顯。對外貿(mào)易對工業(yè)環(huán)境的負向影響主要分布在北京、吉林、遼寧、黑龍江、陜西、甘肅、青海7個省市,這些地區(qū)除北京市外均為經(jīng)濟欠發(fā)達地區(qū)。北京為發(fā)達地區(qū),本身技術(shù)水平較高,對外貿(mào)易在很大程度上為技術(shù)輸出,存在一定的負向影響。其他6個地區(qū)除經(jīng)濟欠發(fā)達外,屬于我國傳統(tǒng)的老工業(yè)基地,主要發(fā)展重工業(yè),技術(shù)水平較低,進出口貿(mào)易主要以重工業(yè)產(chǎn)品為主,這種進出口的技術(shù)效應低于其增加產(chǎn)出而引起的環(huán)境成本,因此對工業(yè)環(huán)境效率的影響呈現(xiàn)出了負向作用。相反,對外貿(mào)易對環(huán)境效率影響為正的地區(qū)主要分布在東部沿海和中部地區(qū),其中福建、湖北、山西3省尤為明顯。FDI對工業(yè)環(huán)境效率影響為負的地區(qū)主要分布在西部地區(qū)和中部地區(qū)各省,其顯著性有待采用面板分位回歸模型進一步檢驗。

分析其他控制變量可見,市場結(jié)構(gòu)對工業(yè)環(huán)境效率影響為負的地區(qū)主要有北京、天津、吉林、山東和寧夏,對其他地區(qū)影響為正向;R&D產(chǎn)生負向影響的主要有山西和新疆等省區(qū);工業(yè)化水平對環(huán)境效率產(chǎn)生負向影響的主要有北京、天津、河北、山西、甘肅、寧夏、青海和新疆8個省、市、自治區(qū),在工業(yè)化水平較低的地區(qū),工業(yè)化對環(huán)境效率的影響為負。而在工業(yè)化水平較高的地區(qū),工業(yè)化對環(huán)境效率的影響為正,較符合工業(yè)發(fā)展的實際情況,可認為工業(yè)化水平有利于工業(yè)技術(shù)的提高。

3.3 面板分位回歸分析

對變量繼續(xù)進行面板分位回歸分析,估計結(jié)果見表4。

表4 面板分位回歸估計結(jié)果

注:***、**、*分別表示在1%、5%、10%顯著水平下顯著；括號內(nèi)為估計系數(shù)的標準誤。

分析分為回歸估計結(jié)果可見,對外貿(mào)易對工業(yè)環(huán)境效率存在著顯著的正向影響,且隨著分位點的提高,這種影響在逐漸增強,表明在工業(yè)環(huán)境效率加高的地區(qū),對外貿(mào)易對工業(yè)環(huán)境效率的影響越大,對外貿(mào)易是工業(yè)技術(shù)進步的重要來源;FDI對工業(yè)環(huán)境效率的影響梯度顯著,在中位數(shù)及以下。即環(huán)境效率較低的省、市、自治區(qū),FDI對工業(yè)環(huán)境效率的影響為負;在環(huán)境效率較高的地區(qū),FDI對工業(yè)環(huán)境效率的影響為正。結(jié)合本文的工業(yè)環(huán)境效率的空間分布規(guī)律發(fā)現(xiàn),影響為負的地區(qū)主要為中西部地區(qū),影響為正的地區(qū)主要為東部沿海地區(qū),與上述GTWR模型的估計結(jié)果一致,但估計結(jié)果的統(tǒng)計顯著水平上不顯著。佟連軍等[11]的觀點認為,外商投資一方面給本地工業(yè)企業(yè)帶來了先進技術(shù),促進了科技合作;另一方面又引進了發(fā)達國家高耗能、高污染的產(chǎn)業(yè)。兩方面的共同作用給工業(yè)環(huán)境效率帶來正、負兩方面的影響,兩種影響相互抵消,總體上外商直接投資對工業(yè)環(huán)境效率的影響并不顯著在一定意義上成立。值得注意的是,盡管外商直接投資對工業(yè)環(huán)境效率的影響不顯著,但在工業(yè)環(huán)境效率較低的中西部地區(qū),其影響為負,在工業(yè)環(huán)境效率較高的東部地區(qū),其影響為正。即在西部地區(qū)外商直接投資的環(huán)境成本效應加大,高污染、高耗能產(chǎn)業(yè)引進較多;在東部地區(qū)外商直接投資的技術(shù)效應較大,外商投資引進以技術(shù)為主,因此不能否認“污染天堂”假說在我國存在。其他變量如市場結(jié)構(gòu)、技術(shù)投入、工業(yè)化水平對工業(yè)環(huán)境效率的影響均為正向,存在著一定的分布梯度,與GTWR模型的估計結(jié)果基本一致。因此,各地在改善環(huán)境、提高工業(yè)環(huán)境效率時需注重上述各因素的影響。在我國工業(yè)環(huán)境效率的分析中,空間異質(zhì)性問題不可忽視。

綜上所述,對外貿(mào)易、FDI與我國工業(yè)環(huán)境效率間的關(guān)系為:從我國整體和工業(yè)效率分布類別看,對外貿(mào)易對我國工業(yè)環(huán)境效率的影響為正，經(jīng)濟越發(fā)達的地區(qū),對外貿(mào)易對工業(yè)環(huán)境效率的影響越大；從各省看,在技術(shù)落后的重工業(yè)發(fā)展地區(qū)和技術(shù)輸出型地區(qū),對外貿(mào)易對工業(yè)環(huán)境效率的影響為負。FDI對工業(yè)環(huán)境效率的影響存在梯度分布，整體看,技術(shù)效應和環(huán)境成本效應的相互作用,在各省的影響不顯著。

4 結(jié)論及啟示

本文運用空間描述統(tǒng)計、超效率DEA模型對我國對外貿(mào)易、FDI及工業(yè)環(huán)境效率的空間分布進行了分析。結(jié)果表明:①從整體分布看,我國對外貿(mào)易活動與工業(yè)外商直接投資活動主要集中在東部地區(qū),目前仍表現(xiàn)為“東強西弱”;從時間演變看,對外貿(mào)易活動表現(xiàn)為繼續(xù)向東集聚,呈現(xiàn)出“強者愈強,弱者愈弱”的發(fā)展態(tài)勢,外商直接投資活動逐漸在向中西部轉(zhuǎn)移擴散。此外,在對外貿(mào)易活動和外商直接投資活動中,個別省份異于周邊地區(qū)。②我國工業(yè)環(huán)境效率表現(xiàn)為明顯的梯度分布,東部地區(qū)明顯高于中西部地區(qū),但整體環(huán)境效率由東向西在逐漸改善。在此演變過程中，東、中、西部存在工業(yè)環(huán)境效率變化較顯著的省份。③對外貿(mào)易、FDI對我國工業(yè)環(huán)境效率的影響存在一定分布梯度,FDI的影響更顯著。從我國整體和工業(yè)效率分布類別看,對外貿(mào)易對我國工業(yè)環(huán)境效率的影響為正,且經(jīng)濟越發(fā)達的地區(qū),對外貿(mào)易對工業(yè)環(huán)境效率的影響越大。但從各省份看,在技術(shù)落后的重工業(yè)發(fā)展地區(qū)和技術(shù)輸出型地區(qū),對外貿(mào)易對工業(yè)環(huán)境效率的影響為負。④在經(jīng)濟發(fā)達、工業(yè)環(huán)境效率較高的地區(qū),FDI對工業(yè)環(huán)境效率具有一定的促進作用,可提高地區(qū)工業(yè)環(huán)境效率；在經(jīng)濟相對落后、工業(yè)環(huán)境效率較低的地區(qū),FDI會導致環(huán)境進一步惡化,抑制工業(yè)環(huán)境效率的提高。盡管FDI的技術(shù)效應和環(huán)境成本效應在相互作用中的正負作用不顯著,但不能否認“污染天堂”假說在我國的存在。此外,市場結(jié)構(gòu)、技術(shù)投入、工業(yè)化水平對我國工業(yè)環(huán)境效率均存在影響,且影響的分布梯度明顯。優(yōu)化市場結(jié)構(gòu)、合理化技術(shù)投資和提高工業(yè)化水平是工業(yè)環(huán)境效率須考慮的問題。在我國整體工業(yè)環(huán)境效率的研究中,工業(yè)環(huán)境效率的空間異質(zhì)性不容忽視。

基于上述結(jié)論,首先在“一帶一路”倡議的發(fā)展背景下,我國堅持“走出去”和“引進來”并重的發(fā)展策略,我國整體貿(mào)易格局、外商直接投資格局正在發(fā)生著變化,各地區(qū)在加大對外開放程度時需注意對環(huán)境的潛在影響,合理引進發(fā)展成熟的高清潔技術(shù),提高區(qū)域的整體工業(yè)技術(shù)水平,保護環(huán)境。注重發(fā)揮市場的作用,建立完善的市場機制,確保低效、污染的國有企業(yè)退出市場,高效、清潔的民營企業(yè)得到平等的發(fā)展機會。其次,在“中國制造2025”的發(fā)展背景下,我國制造業(yè)產(chǎn)業(yè)格局正發(fā)生顯著的變化,東部制造業(yè)產(chǎn)業(yè)正在逐步向中西部轉(zhuǎn)移,轉(zhuǎn)移的產(chǎn)業(yè)往往以勞動密集型和資源密集型為主。一方面,相對于中西部地區(qū),這些轉(zhuǎn)移產(chǎn)業(yè)具有較高的技術(shù)效率,在一定程度上能提高工業(yè)環(huán)境效率;另一方面,長期來看資源密集型產(chǎn)業(yè)的發(fā)展會大量開發(fā)當?shù)刈匀毁Y源,造成環(huán)境破壞,降低工業(yè)環(huán)境效率,因此承接的產(chǎn)業(yè)需根據(jù)地區(qū)發(fā)展狀況及時轉(zhuǎn)型或進行技術(shù)引進,同時當?shù)卣畱⑼晟频沫h(huán)境產(chǎn)權(quán)市場與法律體系。