京津冀環(huán)境污染時(shí)空演進(jìn)研究

李從欣 李國(guó)柱 崔文靜

基金項(xiàng)目：基金項(xiàng)目河北省社會(huì)科學(xué)基金項(xiàng)目《京津冀環(huán)境污染的時(shí)空演進(jìn)及防治對(duì)策研究》（HB18YJ013）。

作者簡(jiǎn)介：作者簡(jiǎn)介李從欣（1974—），女，河北石家莊人，博士，教授、碩士生導(dǎo)師，主要研究方向?yàn)榻y(tǒng)計(jì)方法及應(yīng)用、生態(tài)[JP2]統(tǒng)計(jì);李國(guó)柱（1971—），男，河北滄州人，博士，教授、碩士生導(dǎo)師，主要研究方向?yàn)閿?shù)量經(jīng)濟(jì)、能源環(huán)境統(tǒng)計(jì);崔文靜[JP]（1981—），女，河北石家莊人，碩士，副教授，主要研究方向?yàn)楹暧^經(jīng)濟(jì)。

[摘要]使用2004—2017年京津冀13個(gè)城市的面板數(shù)據(jù)，采用時(shí)空地理加權(quán)回歸模型，對(duì)京津冀環(huán)境污染時(shí)空演進(jìn)影響因素異質(zhì)性進(jìn)行了分析。研究結(jié)果表明，產(chǎn)業(yè)結(jié)構(gòu)和產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染時(shí)空演進(jìn)影響的回歸系數(shù)各不相同，存在顯著的時(shí)空異質(zhì)性。產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染時(shí)空演進(jìn)影響差異巨大，不同城市產(chǎn)業(yè)結(jié)構(gòu)調(diào)整存在巨大的潛力。產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染時(shí)空演進(jìn)的影響存在顯著的異質(zhì)性，在利益和社會(huì)責(zé)任選擇方面，外部環(huán)境壓力政策更有利企業(yè)承擔(dān)社會(huì)責(zé)任。

[關(guān)鍵詞]關(guān)鍵詞環(huán)境污染;時(shí)空演進(jìn);異質(zhì)性;時(shí)空地理加權(quán)回歸模型

[中圖分類(lèi)號(hào)]中圖分類(lèi)號(hào)F127;X321[文獻(xiàn)標(biāo)識(shí)碼]A文獻(xiàn)標(biāo)志碼[文章編號(hào)]1673-0461（2020）08-0061-05

京津冀由于政治地位、經(jīng)濟(jì)發(fā)展水平不同，產(chǎn)業(yè)結(jié)構(gòu)和工業(yè)經(jīng)濟(jì)存在著顯著的差異。在社會(huì)發(fā)展進(jìn)程中，工業(yè)經(jīng)濟(jì)發(fā)展是生態(tài)環(huán)境的主要污染源，工業(yè)“三廢”（廢水、廢氣和固體廢棄物）對(duì)環(huán)境的影響成為政府各部門(mén)和學(xué)術(shù)界關(guān)注的熱點(diǎn)。其中，工業(yè)廢水由于含有有毒、有害物質(zhì)，對(duì)環(huán)境危害更大，直接威脅生態(tài)安全和人類(lèi)健康?！吨攸c(diǎn)流域水污染防治規(guī)劃（20162020年）》指出海河流域的北京、 天津、廊坊、石家莊、邢臺(tái)、保定、滄州等地環(huán)境污染問(wèn)題相對(duì)突出，要求加大對(duì)這些地區(qū)的造紙、 焦化、 印染、皮革等產(chǎn)業(yè)結(jié)構(gòu)和布局的調(diào)整力度。因此探索京津冀環(huán)境污染時(shí)空演進(jìn)及影響因素的異質(zhì)性，對(duì)提升京津冀生態(tài)環(huán)境具有積極的現(xiàn)實(shí)意義。

一、文獻(xiàn)綜述

對(duì)環(huán)境污染時(shí)空演進(jìn)及影響因素異質(zhì)性的研究主要集中在以下兩個(gè)方面：

第一，環(huán)境污染影響因素的研究。工業(yè)排放、汽車(chē)尾氣、燃煤排放等因素是造成環(huán)境污染的主要因素。第二產(chǎn)業(yè)主要由工業(yè)組成，工業(yè)污染是環(huán)境污染的主要來(lái)源[1]。莊汝龍等 [2]和沈曉梅等 [3]發(fā)現(xiàn)產(chǎn)業(yè)結(jié)構(gòu)與環(huán)境污染存在正相關(guān)關(guān)系，第二產(chǎn)業(yè)比重促進(jìn)工業(yè)廢水排放。部分國(guó)內(nèi)學(xué)者研究了產(chǎn)業(yè)集聚對(duì)環(huán)境污染的影響，閆逢柱等[4]、劉習(xí)平等[5]通過(guò)實(shí)證分析等發(fā)現(xiàn)產(chǎn)業(yè)集聚有利于環(huán)境質(zhì)量的改善，減少環(huán)境污染;周明生等[6]發(fā)現(xiàn)產(chǎn)業(yè)集聚與環(huán)境污染存在顯著的倒“U”型關(guān)系;徐輝等[7]發(fā)現(xiàn)產(chǎn)業(yè)專(zhuān)業(yè)化集聚對(duì)環(huán)境污染有先抑制后促進(jìn)的作用;劉軍等[8]發(fā)現(xiàn)產(chǎn)業(yè)集聚增加了環(huán)境污染。本文進(jìn)行分析時(shí)，除了選取產(chǎn)業(yè)結(jié)構(gòu)作為核心解釋變量外，還選取了產(chǎn)業(yè)盈利能力作為核心指標(biāo)，盈利是企業(yè)持續(xù)經(jīng)營(yíng)的目標(biāo)之一，隨著獲利能力的提高，企業(yè)是否更愿意承擔(dān)社會(huì)責(zé)任，進(jìn)行綠色技術(shù)創(chuàng)新，為減輕環(huán)境污染作出貢獻(xiàn)

第二，影響因素異質(zhì)性分析使用方法的選擇。Du等[11]使用空間計(jì)量分析了PM25影響因素的異質(zhì)性。這類(lèi)方法需要對(duì)樣本進(jìn)行空間自相關(guān)檢驗(yàn)，計(jì)算空間權(quán)重矩陣，其得出來(lái)的回歸系數(shù)是固定的;李丹丹等[12] 采用地理加權(quán)回歸模型（GWR）對(duì)中國(guó)省域碳排放影響因素時(shí)空異質(zhì)性進(jìn)行了研究，結(jié)果證明影響因素存在明顯的時(shí)空異質(zhì)性;李建豹等[13]使用地理加權(quán)回歸模型（GWR）對(duì)中國(guó)省域碳排放影響因素的空間異質(zhì)性進(jìn)行了分析，得出的結(jié)論與李丹丹等人的結(jié)論相同，這種分析方法得出的系數(shù)不再是固定的，不同區(qū)域回歸模型的系數(shù)隨空間變化而變化，但是時(shí)間維度系數(shù)的變化卻沒(méi)有體現(xiàn)出來(lái)。

通過(guò)對(duì)上述文獻(xiàn)的總結(jié)，這些學(xué)者之所以有不同的結(jié)論，與他們選取的指標(biāo)、樣本，采用的方法有很大關(guān)系。京津冀環(huán)境污染的變化既有時(shí)間上的，也有空間上的，本文進(jìn)行分析時(shí)采用的是時(shí)空地理加權(quán)回歸模型，回歸系數(shù)不僅隨空間變化，而且隨著時(shí)間維度的變化，回歸系數(shù)也會(huì)不同。通過(guò)對(duì)京津冀環(huán)境污染時(shí)空演進(jìn)影響因素異質(zhì)性研究，探索產(chǎn)業(yè)結(jié)構(gòu)和產(chǎn)業(yè)盈利能力在不同時(shí)間和空間層面上對(duì)環(huán)境污染的影響機(jī)理，為京津冀因地制宜地制定減輕環(huán)境污染的政策提供較為科學(xué)的數(shù)據(jù)支持。

二、變量說(shuō)明和理論模型的構(gòu)建

（一）變量說(shuō)明

解釋變量：環(huán)境污染（Enp）。對(duì)環(huán)境污染的測(cè)度主要有兩種方法：一是環(huán)境污染綜合指數(shù)。如許和連等[14]和郭四代等[15]使用熵值法測(cè)算環(huán)境污染，測(cè)算結(jié)果為指數(shù)值越大，表明環(huán)境污染愈嚴(yán)重。這種方法雖然把廢水、廢氣和固體廢棄物等都考慮進(jìn)去了，但進(jìn)行影響因素分析時(shí)，卻無(wú)法分析不同因素對(duì)水污染、空氣污染的影響狀況。二是單一指標(biāo)。這類(lèi)指標(biāo)是年鑒公布的指標(biāo)，較為客觀，如金春雨等[16]和李從欣等[17]使用二氧化硫排放量測(cè)量環(huán)境污染。本文主要分析的是工業(yè)廢水排放對(duì)環(huán)境造成的污染，因此使用工業(yè)廢水排放量（萬(wàn)噸）作為環(huán)境污染的代理變量。

參考許正松等[18]和戴宏偉等[19]的文獻(xiàn)，本文以產(chǎn)業(yè)結(jié)構(gòu)（Ins）和產(chǎn)業(yè)盈利能力（Inp）作為核心解釋變量。工業(yè)是京津冀地區(qū)經(jīng)濟(jì)發(fā)展的主要支柱產(chǎn)業(yè)，工業(yè)類(lèi)型主要集中在煉鋼、造紙、紡織以及化學(xué)品制造等相關(guān)方面，這些行業(yè)最主要的特征是高排污、高排放。產(chǎn)業(yè)結(jié)構(gòu)調(diào)整一直是京津冀協(xié)同發(fā)展的主題。產(chǎn)業(yè)結(jié)構(gòu)通過(guò)第二產(chǎn)業(yè)產(chǎn)值占地區(qū)生產(chǎn)總值的比重（%）測(cè)量，產(chǎn)業(yè)盈利能力通過(guò)規(guī)模以上工業(yè)企業(yè)利潤(rùn)總額除以工業(yè)企業(yè)個(gè)數(shù)（萬(wàn)元/個(gè)）測(cè)量。

控制變量：人口密度（Pod）、經(jīng)濟(jì)發(fā)展水平（Pgdp）、科技支出強(qiáng)度（Tsi）和利用外資程度（Dfu）。人口密度2004—2015年數(shù)據(jù)來(lái)自年鑒，2016年和2017年的數(shù)據(jù)通過(guò)年末人口數(shù)除以行政區(qū)域土地面積（人/平方公里）測(cè)量，經(jīng)濟(jì)發(fā)展水平以人均地區(qū)生產(chǎn)總值（元/人）測(cè)量，科技支出強(qiáng)度用科學(xué)技術(shù)支出除以地區(qū)生產(chǎn)總值測(cè)量，利用外資程度用當(dāng)年實(shí)際使用外資金額（以當(dāng)年匯率折算成人民幣）除以地區(qū)生產(chǎn)總值測(cè)量。

（二）理論模型的構(gòu)建

本文研究的內(nèi)容是京津冀空氣污染時(shí)空演進(jìn)影響因素的異質(zhì)性，使用的數(shù)據(jù)是面板數(shù)據(jù)，采用的研究方法是Huang等（2010）[20]的時(shí)空地理加權(quán)回歸模型（GTWR，Geographically and Temporally Weighted Regression）。該方法既可以分析京津冀空氣污染時(shí)空演進(jìn)影響因素隨地理位置變動(dòng)的異質(zhì)性，又可以分析影響因素隨時(shí)間變動(dòng)的異質(zhì)性。模型的具體形式如公式（1）：

（三）數(shù)據(jù)來(lái)源

工業(yè)廢水排放量、第二產(chǎn)業(yè)占地區(qū)生產(chǎn)總值的比重、規(guī)模以上工業(yè)企業(yè)個(gè)數(shù)、規(guī)模以上工業(yè)企業(yè)利潤(rùn)總額、人口密度、人均地區(qū)生產(chǎn)總值、科學(xué)技術(shù)支出和當(dāng)年實(shí)際使用外資金額等數(shù)據(jù)來(lái)源于2005—2018年的《中國(guó)城市統(tǒng)計(jì)年鑒》。數(shù)據(jù)分析使用的軟件是ArcGIS102和SPSS220。

三、京津冀環(huán)境污染時(shí)空演進(jìn)影響因素時(shí)空異質(zhì)性實(shí)證研究

（一）京津冀環(huán)境污染時(shí)空演進(jìn)影響因素異質(zhì)性研究

工業(yè)廢水排放量為因變量，產(chǎn)業(yè)結(jié)構(gòu)和產(chǎn)業(yè)盈利能力為核心變量，代入公式（4）進(jìn)行估計(jì)，由于參數(shù)估計(jì)結(jié)果較多，無(wú)法一一列出，本文列出部分結(jié)果如表1所示。

從表1可知，模型的判定系數(shù)為0806 6，調(diào)整的判定系數(shù)是0800 0，模型擬合程度較好。產(chǎn)業(yè)結(jié)構(gòu)、產(chǎn)業(yè)盈利能力、人口密度、經(jīng)濟(jì)發(fā)展水平、科技支出強(qiáng)度、利用外資程度對(duì)環(huán)境污染時(shí)空影響的五欄參數(shù)估計(jì)結(jié)果數(shù)值明顯不同。根據(jù)對(duì)應(yīng)的P值全部為0，小于005，說(shuō)明六個(gè)影響因素對(duì)空氣污染時(shí)空演進(jìn)存在明顯的異質(zhì)性。產(chǎn)業(yè)結(jié)構(gòu)和產(chǎn)業(yè)盈利能力的回歸系數(shù)取最小值時(shí)，可以降低環(huán)境污染，而其他四欄回歸系數(shù)都會(huì)使環(huán)境污染上升。

（二）產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染時(shí)空演進(jìn)異質(zhì)性分析

圖1是產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染時(shí)空演進(jìn)影響系數(shù)估計(jì)結(jié)果圖。從時(shí)間維度上可知，產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染影響呈“UM”型的是邯鄲和邢臺(tái)，產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染影響呈“M”型的是衡水、石家莊和保定，產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染影響呈倒“草帽”型的是北京、廊坊、滄州和張家口，產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染影響呈倒“N”型的是天津、唐山、秦皇島和承德。從空間維度上可知，2005年，京津冀13個(gè)城市產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染影響程度差距非常大。京津冀西南方位城市的產(chǎn)業(yè)結(jié)構(gòu)回歸系數(shù)是負(fù)值，東北方位城市的產(chǎn)業(yè)結(jié)構(gòu)回歸系數(shù)是正值，且數(shù)值差異非常明顯，西南方位的石家莊、保定、衡水、邢臺(tái)和邯鄲的產(chǎn)業(yè)結(jié)構(gòu)有利于環(huán)境污染降低，而東北方位的城市如秦皇島、承德、唐山和天津卻不利于環(huán)境污染的降低。到2010年，京津冀13個(gè)城市產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染影響程度差距雖然縮小，但此時(shí)13個(gè)城市的產(chǎn)業(yè)結(jié)構(gòu)均造成環(huán)境污染上升。2015年，京津冀13個(gè)城市的產(chǎn)業(yè)結(jié)構(gòu)的回歸系數(shù)的最大值和最小值均小于2010年產(chǎn)業(yè)結(jié)構(gòu)回歸系數(shù)的最大值和最小值，說(shuō)明2015年產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染的影響整體上比2010年有所下降。到2017年，13個(gè)城市的產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染的回歸系數(shù)與2015年相比均出現(xiàn)下降趨勢(shì)，下降最明顯的是邢臺(tái)和邯鄲，這兩個(gè)城市的產(chǎn)業(yè)結(jié)構(gòu)使環(huán)境污染降低。

圖題圖12004—2017年產(chǎn)業(yè)結(jié)構(gòu)對(duì)京津冀環(huán)境污染時(shí)空演進(jìn)影響異質(zhì)性特征

（三）產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染時(shí)空演進(jìn)異質(zhì)性分析

圖2是產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染時(shí)空演進(jìn)影響系數(shù)估計(jì)結(jié)果圖。從時(shí)間維度上可知，產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染影響呈明顯“M”型的是邯鄲、邢臺(tái)和石家莊，產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染影響呈平緩“M”型的是秦皇島和衡水，產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染影響呈倒“草帽”型的是北京、天津、廊坊和承德，產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染影響呈“W”型的是張家口、唐山、保定和滄州。從空間維度上可知，2005年產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染時(shí)空演進(jìn)影響系數(shù)差距較大，邯鄲、邢臺(tái)和石家莊三城市的產(chǎn)業(yè)盈利能力有利于環(huán)境污染的降低;到2010年時(shí)，13個(gè)城市的產(chǎn)業(yè)盈利能力回歸系數(shù)雖然存在顯著差異，但它們之間的差距與2005年相比顯著縮小;2010年到2015年，京津冀13個(gè)城市產(chǎn)業(yè)盈利能力的回歸系數(shù)差距幾乎在相同的范圍內(nèi);到2015年，產(chǎn)業(yè)盈利能力回歸系數(shù)的差距又開(kāi)始擴(kuò)大，除秦皇島的產(chǎn)業(yè)盈利能力能降低環(huán)境污染外，其它城市均使環(huán)境污染上升，使環(huán)境污染上升幅度最大的三個(gè)城市分別是邯鄲、邢臺(tái)和石家莊;到2017年，產(chǎn)業(yè)盈利能力回歸系數(shù)最大的三個(gè)城市仍然是邯鄲、邢臺(tái)和石家莊。

圖題圖22004—2017年產(chǎn)業(yè)盈利能力對(duì)京津冀環(huán)境污染[JZ]時(shí)空演進(jìn)影響異質(zhì)性特征

（四）穩(wěn)健性檢驗(yàn)

為了保證分析結(jié)果的穩(wěn)健性，本文以工業(yè)二氧化硫排放量為因變量，產(chǎn)業(yè)結(jié)構(gòu)、產(chǎn)業(yè)盈利能力為核心變量，人口密度、經(jīng)濟(jì)發(fā)展水平、科技支出強(qiáng)度和利用外資程度為控制變量進(jìn)了檢驗(yàn)，檢驗(yàn)結(jié)果與前面分析結(jié)果大體一致。由于篇幅有限，沒(méi)有一一列出，通過(guò)穩(wěn)健性檢驗(yàn)，本文的研究結(jié)論是可靠的。

四、研究結(jié)論與討論

本文使用2004—2017年京津冀13個(gè)城市的面板數(shù)據(jù)，采用GTWR模型，分析了產(chǎn)業(yè)結(jié)構(gòu)和產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染時(shí)空演進(jìn)影響的異質(zhì)性，目標(biāo)是探索各影響因素在不同時(shí)間和空間上對(duì)環(huán)境污染影響的差異性。通過(guò)實(shí)證分析得出了以下結(jié)論。

第一，通過(guò)時(shí)空地理加權(quán)回歸模型可以發(fā)現(xiàn)，產(chǎn)業(yè)結(jié)構(gòu)和產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染時(shí)空演進(jìn)影響的回歸系數(shù)各不相同，存在顯著的時(shí)空異質(zhì)性。產(chǎn)業(yè)結(jié)構(gòu)回歸系數(shù)的最小值和最大值差距接近2 000，說(shuō)明不同時(shí)間和不同空間上產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染的影響差異非常大，產(chǎn)業(yè)結(jié)構(gòu)調(diào)整存在著很大的潛力，京津冀不同城市應(yīng)當(dāng)實(shí)行差異化的產(chǎn)業(yè)結(jié)構(gòu)調(diào)整政策。產(chǎn)業(yè)盈利能力回歸系數(shù)的最小值和最大值差距在20左右，說(shuō)明不同時(shí)間和不同空間上產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染的影響差異不是太大，不同地區(qū)不同年份產(chǎn)業(yè)盈利能力回歸系數(shù)是負(fù)值的占總樣本的126%，說(shuō)明京津冀13個(gè)城市在承擔(dān)社會(huì)責(zé)任方面還有所欠缺，不同城市的相關(guān)領(lǐng)導(dǎo)部門(mén)應(yīng)加大引導(dǎo)力度，使企業(yè)的盈利目標(biāo)和社會(huì)責(zé)任同等重要。

第二，產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染時(shí)空演進(jìn)的影響存在顯著的異質(zhì)性。從2005年、2010年和2015年產(chǎn)業(yè)結(jié)構(gòu)回歸系數(shù)來(lái)看，京津冀西南方位的城市產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染的變動(dòng)幅度最大，由降低環(huán)境污染轉(zhuǎn)變?yōu)槭弓h(huán)境污染快速上升。京津冀東北方位的城市產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染的影響整體趨于下降的趨勢(shì)，這些城市主要位于北京的東邊。說(shuō)明當(dāng)環(huán)境污染成為全社會(huì)防控主題時(shí)，北京及其東邊的城市能快速的做出反應(yīng)，率先降低環(huán)境污染。2015年后京津冀環(huán)境污染的壓力由東北方位的城市轉(zhuǎn)向西南方位的城市，隨著京津冀環(huán)境污染治理一體化的提出，在京津冀環(huán)境污染共同治理壓力下，西南方位的城市如石家莊、邢臺(tái)和邯鄲，產(chǎn)業(yè)結(jié)構(gòu)對(duì)環(huán)境污染的影響在2015年后出現(xiàn)下降趨勢(shì)。

第三，產(chǎn)業(yè)盈利能力對(duì)環(huán)境污染時(shí)空演進(jìn)的影響存在顯著的異質(zhì)性。石家莊、邢臺(tái)和邯鄲在2005年時(shí)產(chǎn)業(yè)盈利能力回歸系數(shù)正向效應(yīng)排在前三位，在降低環(huán)境污染方面屬于佼佼者，到了2017年這三個(gè)城市的產(chǎn)業(yè)盈利能力回歸系數(shù)負(fù)向效應(yīng)排在前三位，在增加環(huán)境污染方面又成為領(lǐng)頭羊，石家莊、邢臺(tái)和邯鄲的企業(yè)在承擔(dān)社會(huì)責(zé)任方面有所欠缺。2010年到2015年京津冀13個(gè)城市產(chǎn)業(yè)盈利能力回歸系數(shù)雖然有顯著異質(zhì)性，但變化區(qū)間不大，主要跟“十二五”節(jié)能減排有很大關(guān)系，在限電、限產(chǎn)的強(qiáng)大壓力下，所有城市的企業(yè)社會(huì)責(zé)任增強(qiáng)，愿為降低能耗，減輕環(huán)境污染做貢獻(xiàn)。2016年后，節(jié)能減排政策有所松動(dòng)，大部分城市的產(chǎn)業(yè)盈利能力回歸系數(shù)開(kāi)始上升，利益驅(qū)動(dòng)超過(guò)社會(huì)責(zé)任。因此在利益和社會(huì)責(zé)任選擇方面，外部環(huán)境壓力政策更有利于企業(yè)承擔(dān)社會(huì)責(zé)任。

[參考文獻(xiàn)]

[1][ZK（#]張可，汪東芳.經(jīng)濟(jì)集聚與環(huán)境污染的交互影響及空間溢出[J].中國(guó)工業(yè)經(jīng)濟(jì)，2014（6）：70-82.

[2]莊汝龍，宓科娜，梁龍武.中國(guó)工業(yè)廢水排放格局及其驅(qū)動(dòng)因素[J].長(zhǎng)江流域資源與環(huán)境，2018，27（8）：1765-1775.

[3]沈曉梅，胡凱莉，盛前，等.技術(shù)進(jìn)步視角下產(chǎn)業(yè)結(jié)構(gòu)對(duì)工業(yè)廢水排放的門(mén)檻效應(yīng)研究——基于省際面板數(shù)據(jù)的實(shí)證檢驗(yàn)[J].水利經(jīng)濟(jì)，2020，38（1）：8-13，85.

[4]閆逢柱，蘇李，喬娟.產(chǎn)業(yè)集聚發(fā)展與環(huán)境污染關(guān)系的考察——來(lái)自中國(guó)制造業(yè)的證據(jù)[J].科學(xué)學(xué)研究，2011（1）：79-83.

[5]劉習(xí)平，宋德勇.城市產(chǎn)業(yè)集聚對(duì)城市環(huán)境的影響[J].城市問(wèn)題，2013（3）：9-14.

[6]周明生，王帥.產(chǎn)業(yè)集取是導(dǎo)致區(qū)域環(huán)境污染的“兇手”嗎：來(lái)自京津冀地區(qū)的證據(jù)[J].經(jīng)濟(jì)體制改革，2018（5）：185-190.

[7]徐輝，楊燁.人口和產(chǎn)業(yè)集聚對(duì)環(huán)境污染的影響——以中國(guó)的100 個(gè)城市為例[J].城市問(wèn)題，2017（1）：53-60.

[8]劉軍，程中華，李廉水.產(chǎn)業(yè)聚集與環(huán)境污染[J].科研管理，2016（6）：134-140.

[9]李從欣，李國(guó)柱.我國(guó)區(qū)域科技創(chuàng)新經(jīng)濟(jì)效益影響因素研究[J].廣西社會(huì)科學(xué)，2017（7）：76-80.

[10][ZK（#]LI [JP3]C X，LI G Z，LI Z N.The threshold effect of the efficiency of science and technological services on regional environmental governance in China[J].Growth and change.2019，50（6）：1-17.[JP]

[11]DU Y Y，SUN T S，PENG J，et al.Direct and spillover effects of urbanization on PM2.5 concentrations in China s top three urban agglomerations[J].Journal of cleaner production，2018（190）：72-83.

[12]李丹丹，劉銳，陳動(dòng).中國(guó)省域碳排放及其驅(qū)動(dòng)因子的時(shí)空異質(zhì)性研究[J].中國(guó)人口·資源與環(huán)境，2013，23（7）：84-92.

[13]李建豹，黃賢金，吳常艷，等.中國(guó)省域碳排放影響因素的空間異質(zhì)性分析[J].經(jīng)濟(jì)地理，2015，35（11）：21-28.

[14]許和連，鄧玉萍.外商直接投資導(dǎo)致了中國(guó)的環(huán)境污染了嗎？——基于中國(guó)省際面板數(shù)據(jù)的空間計(jì)量研究[J].管理世界，2012（2）：30-43.

[15]郭四代，張華，郭杰，等.基于空間計(jì)量模型的中國(guó)環(huán)境污染評(píng)價(jià)及影響因素分析[J].生態(tài)學(xué)雜志，2018，37（2）：471-481.

[16]金春雨，吳安兵.工業(yè)經(jīng)濟(jì)結(jié)構(gòu)、經(jīng)濟(jì)增長(zhǎng)對(duì)環(huán)境污染的非線性影響[J].中國(guó)人口·資源與環(huán)境，2017，27（10）：64-73.

[17]李從欣，李子寧，李國(guó)柱.區(qū)域環(huán)境污染的時(shí)空演進(jìn)研究[J].統(tǒng)計(jì)與管理，2019（9）：39-42.

[18]許正松，孔凡斌.經(jīng)濟(jì)發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)與環(huán)境污染：基于江西省的實(shí)證分析[J].當(dāng)代財(cái)經(jīng)，2014（8）：15-20.

[19]戴宏偉，回瑩.京津冀霧霾污染與產(chǎn)業(yè)結(jié)構(gòu)、城鎮(zhèn)化水平的空間效應(yīng)研究[J].經(jīng)濟(jì)理論與經(jīng)濟(jì)管理，2019（5）：4-19.

[20] HUANG B，WU B，BARRYC M.Geographically and temporally weighted regression for modeling spatiotemporal variation in house prices[J].International journal of geographical information science，2010，24（3）：383-401.

Study on the Spatiotemporal Evolutionof Environmental Pollution in BeijingTianjinHebei：

Spatiotemporal Heterogeneity Based on Influence Factors

Li Congxin，Li Guozhu，Cui Wenjing

（School of Economics and Trade，Hebei GEO University，Shijiazhuang 050031，China）

Abstract：

Serious environmental pollution in BeijingTianjinHebei region directly threatens the ecological security and human health.Based on the panel data from 2004 to 2017，this paper analyzes the heterogeneity of influence factors on the spatiotemporal evolution of environmental pollution in BeijingTianjinHebei region using geographically and temporally weighted regression （GTWR）.The results show that the regression coefficients of industrial structure and industrial profitability on the spatiotemporal evolution of environmental pollution are different，and there is significant spatiotemporal heterogeneity.Industrial structure has great influence on the spatial and temporal evolution of environmental pollution，and there is great potential for industrial structure adjustment in different cities.The influence of industrial profitability on the temporal and spatial evolution of environmental pollution has significant heterogeneity.In terms of the choice of interests and social responsibility，the external environmental pressure policy is more conducive to corporate social responsibility.

Key words： environmental pollution;spatiotemporal evolution;heterogeneity;geographically and temporally weighted regression model

責(zé)任編輯（責(zé)任編輯：張麗陽(yáng)）