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

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

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

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

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

1.1 GTWR模型

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

(1)

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

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

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

(2)

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

1.2 面板分位回歸

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

(3)

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

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

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

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

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

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

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

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

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

圖4 2000年FDI空間分布情況

2.2 FDI空間格局分析

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

圖5 2008年FDI空間分布情況

圖6 2015年FDI空間分布情況

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

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

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

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

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

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

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

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

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

3 模型分析

3.1 控制變量

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

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

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

3.2 GTWR模型分析

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

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

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

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

3.3 面板分位回歸分析

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

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

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

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

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

4 結(jié)論及啟示

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

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