文澤宙，熊 磊

（1.湖南科技大學(xué) 商學(xué)院，湖南 湘潭 411201；2.重慶理工大學(xué) 財(cái)會(huì)研究與開(kāi)發(fā)中心，重慶 400054）

近來(lái)年，隨著我國(guó)工業(yè)化和城市化進(jìn)程的快速推進(jìn)，城市大氣污染特別是霧霾污染問(wèn)題日趨惡化，而且覆蓋范圍仍有不斷擴(kuò)張之勢(shì)。2018年全國(guó)338個(gè)地級(jí)及以上城市中，城市環(huán)境空氣質(zhì)量超標(biāo)比重為64.2%，平均超標(biāo)天數(shù)達(dá)20.7%，其中首要污染物為PM2.5和PM10等①數(shù)據(jù)來(lái)源于2018年《中國(guó)生態(tài)環(huán)境狀況公報(bào)》。。一部名為《穹頂之下》的紀(jì)錄片更是將“霧霾”一詞帶進(jìn)了尋常百姓家，也讓公眾逐漸了解到霧霾污染給人們身心健康帶來(lái)的損害，打破“霧霾鎖城”的困境成為新時(shí)代人們最為迫切的期許。與所有環(huán)境污染類(lèi)似，霧霾污染同樣具有典型的負(fù)外部性，市場(chǎng)機(jī)制在治霾問(wèn)題上存在失靈，因此需要政府強(qiáng)有力的規(guī)制政策予以干預(yù)［1］。隨著我國(guó)一系列霧霾治理相關(guān)手段和政策的運(yùn)用和實(shí)施，霧霾污染情況已有所好轉(zhuǎn)，但由于霧霾污染的空間溢出效應(yīng)、各地區(qū)各自為戰(zhàn)等問(wèn)題的存在，使得多地出現(xiàn)霧霾反撲或PM2.5反彈的現(xiàn)象，要想從根本上治理霧霾污染仍任重道遠(yuǎn)。2018年習(xí)近平在全國(guó)生態(tài)環(huán)境保護(hù)大會(huì)上明確要求：“要以空氣質(zhì)量明顯改善為剛性需求，強(qiáng)化聯(lián)防聯(lián)控，基本消除重污染天氣，還老百姓藍(lán)天白云、繁星閃爍”。可見(jiàn)，理清環(huán)境規(guī)制對(duì)城市霧霾污染的異質(zhì)性空間影響，對(duì)于科學(xué)制定和實(shí)施行之有效的環(huán)境規(guī)制政策，徹底打贏藍(lán)天保衛(wèi)戰(zhàn)，進(jìn)而統(tǒng)籌推進(jìn)我國(guó)生態(tài)文明建設(shè)和經(jīng)濟(jì)高質(zhì)量發(fā)展具有重大理論和現(xiàn)實(shí)意義。

霧霾最早爆發(fā)于英美等工業(yè)化國(guó)家，是由于PM等細(xì)顆粒物在空中集聚導(dǎo)致大氣消光強(qiáng)度增加和能見(jiàn)度降低的一種污染現(xiàn)象［2］。隨著全球各國(guó)對(duì)霧霾污染危害性認(rèn)識(shí)的加深，關(guān)于霧霾污染的相關(guān)研究如雨后春筍般涌現(xiàn)。對(duì)于霧霾污染成分的研究，學(xué)術(shù)界普遍認(rèn)同PM2.5是其最主要組成部分，但受制于PM2.5數(shù)據(jù)的缺乏，使用該數(shù)據(jù)在全國(guó)中小尺度區(qū)域展開(kāi)探討的研究成果還十分有限［3］。學(xué)者們關(guān)注的焦點(diǎn)是PM10、總懸浮顆粒物（TSP）、空氣質(zhì)量指數(shù)（AQI）等相對(duì)容易獲得的大氣污染物［4-6］。文獻(xiàn)研究表明：霧霾污染會(huì)影響居民健康［7］、降低政府公信力［8］、削弱城市吸引力［9］等。陳詩(shī)一等［10］也指出霧霾污染嚴(yán)重阻礙了中國(guó)經(jīng)濟(jì)高質(zhì)量發(fā)展，可見(jiàn)霧霾污染已成為新時(shí)代一個(gè)重大的民生問(wèn)題。與此同時(shí)，學(xué)者們針對(duì)霧霾污染的成因進(jìn)行了積極探索，如李賓等［11］以湘潭市城區(qū)為例分析了交通擁堵對(duì)霧霾的影響，得出交通擁堵造成汽車(chē)尾氣排放增加進(jìn)而加大了霧霾污染的結(jié)論。邵帥等［12］在研究城市化與霧霾污染的關(guān)系時(shí)發(fā)現(xiàn)兩者呈現(xiàn)顯著的正相關(guān)性，即城市化加劇了霧霾污染。東童童等［13］在運(yùn)用中國(guó)省級(jí)面板數(shù)據(jù)進(jìn)行研究后指出，工業(yè)勞動(dòng)與資本集聚是造成我國(guó)霧霾污染程度加重的重要推手。

為了應(yīng)對(duì)日益嚴(yán)重的霧霾污染問(wèn)題，學(xué)者們從不同視角探討了環(huán)境規(guī)制對(duì)霧霾污染的影響情況。劉晨躍等［14］從中介效應(yīng)入手，研究發(fā)現(xiàn)環(huán)境規(guī)制通過(guò)直接和間接雙重路徑共同作用于霧霾治理。王書(shū)斌等［15］基于企業(yè)投資偏好視角，運(yùn)用門(mén)檻回歸模型驗(yàn)證了不同環(huán)境規(guī)制對(duì)霧霾脫鉤的影響效應(yīng)，發(fā)現(xiàn)環(huán)境行政管制和環(huán)境污染監(jiān)管比環(huán)境經(jīng)濟(jì)規(guī)制對(duì)霧霾脫鉤效應(yīng)更好。黃壽峰［16］從影子經(jīng)濟(jì)視角，運(yùn)用動(dòng)態(tài)半?yún)?shù)分析方法實(shí)證研究環(huán)境規(guī)制對(duì)霧霾污染的影響，發(fā)現(xiàn)環(huán)境規(guī)制并未有效直接抑制霧霾污染，而且其間接作用反而會(huì)通過(guò)非正式部門(mén)惡化環(huán)境。此外，呂守軍等［17］通過(guò)構(gòu)建地方政府環(huán)境規(guī)制決策模型進(jìn)行研究，結(jié)果表明即使有政府制度安排，地方利益集團(tuán)的博弈均衡仍會(huì)偏離社會(huì)福利最大化目標(biāo)，從而導(dǎo)致PM2.5治理困局。

學(xué)者們關(guān)于環(huán)境規(guī)制與霧霾污染的研究已取得豐富成果，但現(xiàn)有研究多是集中于全國(guó)整體或省級(jí)數(shù)據(jù)的分析，來(lái)自地級(jí)及以上城市的經(jīng)驗(yàn)證據(jù)較少，且基本未考慮環(huán)境規(guī)制的跨區(qū)域影響。此外，雖有Cole等［1，18］個(gè)別學(xué)者考慮了異質(zhì)性環(huán)境規(guī)制作用于霧霾污染的效應(yīng)，但將環(huán)境規(guī)制區(qū)分為正式和非正式2個(gè)方面以探討其對(duì)霧霾污染影響的文獻(xiàn)仍十分缺乏。有鑒于此，本研究在充分考慮環(huán)境規(guī)制的空間溢出效應(yīng)和霧霾污染的空間流動(dòng)性基礎(chǔ)上，選擇中國(guó)城市這一霧霾污染“重災(zāi)區(qū)”作為研究對(duì)象，運(yùn)用空間計(jì)量方法試圖分別從正式、非正式環(huán)境規(guī)制2個(gè)角度揭示霧霾污染減排效果，以期為多種環(huán)境規(guī)制工具的優(yōu)化組合、各地區(qū)霧霾污染聯(lián)防聯(lián)控機(jī)制的合理制定提供一定依據(jù)。

1 研究設(shè)計(jì)與數(shù)據(jù)說(shuō)明

1.1 雙變量LISA模型

為考察空間中不同地理屬性觀測(cè)值之間的相關(guān)關(guān)系，需要運(yùn)用探索性空間數(shù)據(jù)分析方法（ESDA），其旨在通過(guò)測(cè)度莫蘭指數(shù)（Moran’s I）從數(shù)學(xué)上解析地理學(xué)第一定律，即空間上任何要素均存在相關(guān)性，距離近的要素比距離遠(yuǎn)的要素相關(guān)性更高。而傳統(tǒng)單變量莫蘭指數(shù)主要反映的是空間中同一地理屬性觀測(cè)值之間的相關(guān)關(guān)系［19］，聚焦本研究所考察主題，通過(guò)對(duì)傳統(tǒng)單變量莫蘭指數(shù)進(jìn)行改造，構(gòu)造雙變量莫蘭指數(shù)，這樣既能彌補(bǔ)不同地理屬性觀測(cè)值之間關(guān)系的割裂，同時(shí)更能達(dá)到研究目標(biāo)。

雙變量全局空間自相關(guān)（bivariate global moran’s I）的表達(dá)式為：

雙變量局部空間自相關(guān)（bivariate local moran’s I）的表達(dá)式為：

式中：xi、yj分別代表城市i、j的地理屬性觀測(cè)值，x為（正式或非正式）環(huán)境規(guī)制，y為霧霾污染；這里wij代表鄰接權(quán)重矩陣，即城市i和j不相鄰或者為同一地區(qū)時(shí)wij＝0，城市i和j相鄰時(shí)wij＝1；n代表樣本城市總數(shù)；S2代表樣本方差分別代表城市i屬性A與城市j屬性B觀測(cè)值的方差標(biāo)準(zhǔn)化值。

1.2 空間杜賓模型

對(duì)于空間問(wèn)題的研究需要構(gòu)建反映空間因素的空間計(jì)量模型。本研究在空間杜賓模型的基礎(chǔ)上，通過(guò)模型的識(shí)別選擇最優(yōu)形式，進(jìn)而檢驗(yàn)正式、非正式環(huán)境規(guī)制對(duì)霧霾污染的空間影響。空間杜賓模型的具體表達(dá)形式為：

式中：Xit和Yit分別代表t時(shí)期第i個(gè)城市的解釋變量與被解釋變量觀測(cè)值；Kit、Kjt分別代表t時(shí)期第i和j個(gè)城市的控制變量觀測(cè)值；ρ表示被解釋變量的空間自回歸系數(shù)；wij同樣選用鄰接權(quán)重矩陣；α、β、γ、η為待估參數(shù)向量；μi為空間效應(yīng)；νt為時(shí)間效應(yīng)；εit為干擾項(xiàng)。

1.3 變量選取及數(shù)據(jù)來(lái)源

除了前文所闡述的以霧霾污染作為被解釋變量、正式或非正式環(huán)境規(guī)制作為解釋變量，還引入其他可能對(duì)霧霾污染產(chǎn)生一定干擾作用的因素作為控制變量，并與解釋變量和被解釋變量一起納入到回歸分析當(dāng)中。各變量的具體說(shuō)明如下：

F1=-0.018X1+0.298X2-0.259X3-0.271X4-0.087X5+0.249X6-0.049X7+0.139X8-0.009X9+0.054X10+0.113X11

被解釋變量：由于PM2.5被認(rèn)為是霧霾元兇，因此選用PM2.5來(lái)近似衡量城市霧霾污染程度。囿于我國(guó)從2012年開(kāi)始空氣質(zhì)量監(jiān)測(cè)才正式進(jìn)入“PM2.5時(shí)代”，且僅有部分城市有PM2.5統(tǒng)計(jì)數(shù)據(jù)［20］，導(dǎo)致直接運(yùn)用國(guó)內(nèi)統(tǒng)計(jì)數(shù)據(jù)難以達(dá)到研究目標(biāo)。因此，借鑒已有相關(guān)文獻(xiàn)對(duì)該數(shù)據(jù)的處理方法，使用ArcGIS軟件對(duì)哥倫比亞大學(xué)社會(huì)經(jīng)濟(jì)數(shù)據(jù)和應(yīng)用中心（SEDAC）公布的全球PM2.5衛(wèi)星監(jiān)測(cè)柵格數(shù)據(jù)進(jìn)行解析①當(dāng)前哥倫比亞大學(xué)社會(huì)經(jīng)濟(jì)數(shù)據(jù)和應(yīng)用中心（SEDAC）公布的最新PM2.5數(shù)據(jù)僅到2016年。，從而可得到樣本期內(nèi)中國(guó)城市年均PM2.5濃度具體值。

解釋變量：根據(jù)環(huán)境規(guī)制行為主體的不同，將環(huán)境規(guī)制細(xì)分為正式環(huán)境規(guī)制（FER）和非正式環(huán)境規(guī)制（IFER），力求揭示差異化環(huán)境規(guī)制手段對(duì)城市霧霾污染的空間影響。其中，正式環(huán)境規(guī)制的主體是政府和環(huán)保部門(mén)等官方組織，其主要借助強(qiáng)制性的法律法規(guī)等手段勒令或倒逼企業(yè)進(jìn)行霧霾污染等環(huán)境治理，由于環(huán)境治理程度與綠化覆蓋率存在高度正相關(guān)性，故借鑒鄺嫦娥等［21］的做法，使用市轄區(qū)建成區(qū)綠化覆蓋率作為正式環(huán)境規(guī)制的替代變量；非正式環(huán)境規(guī)制的主體則是公眾和環(huán)保NGO等民間組織或團(tuán)體，其主要通過(guò)舉報(bào)、信訪或輿論監(jiān)督等途徑達(dá)到保護(hù)環(huán)境之目的。已有文獻(xiàn)表明：環(huán)境規(guī)制強(qiáng)度會(huì)隨著居民收入水平的變化而改變［22］，因?yàn)榫用袷杖氲奶岣邥?huì)改善生活條件，使得其逐漸由物質(zhì)追求轉(zhuǎn)向精神追求和對(duì)美好生活環(huán)境的需求，居民保護(hù)環(huán)境的積極性更加高漲，故選用城鎮(zhèn)居民人均可支配收入來(lái)衡量非正式環(huán)境規(guī)制。

控制變量：參照已有學(xué)者的研究成果，從霧霾污染源和霧霾治理積極性2個(gè)方面選取對(duì)霧霾污染具有重要影響的控制變量［23-24］。這里所討論的霧霾污染源主要是人為污染源，包括產(chǎn)生工業(yè)廢氣的工業(yè)污染源、產(chǎn)生汽車(chē)尾氣的交通運(yùn)輸污染源以及產(chǎn)生工地?fù)P塵等的其他污染源。首先，工業(yè)污染源的形成與外商投資（FI）和產(chǎn)業(yè)結(jié)構(gòu)（IS）密切相關(guān)，因?yàn)橥馍掏顿Y的質(zhì)量、結(jié)構(gòu)以及環(huán)境準(zhǔn)入門(mén)檻不盡相同，會(huì)對(duì)當(dāng)?shù)仂F霾污染產(chǎn)生不同程度的影響［25］，即可能會(huì)因?yàn)橐M(jìn)外資而增加（或減少）地方工業(yè)污染源，進(jìn)而加?。ɑ蚓徑猓┏鞘徐F霾污染。而產(chǎn)業(yè)結(jié)構(gòu)的高級(jí)化和合理化調(diào)整也會(huì)直接改變工業(yè)污染源的結(jié)構(gòu)和布局，導(dǎo)致霧霾污染總體排放量發(fā)生變化［26］，分別以GDP中實(shí)際使用外商投資總額占比、第三產(chǎn)業(yè)增加值占比來(lái)衡量外商投資和產(chǎn)業(yè)結(jié)構(gòu)；其次，交通運(yùn)輸污染源是霧霾的主要來(lái)源之一，交通條件（TC）的不斷改善導(dǎo)致公共和私人交通工具從而交通運(yùn)輸污染源顯著增加，促使霧霾污染物排放激增，選取公路里程數(shù)表征交通條件；最后，房產(chǎn)建設(shè)（HC）過(guò)程中產(chǎn)生的工地?fù)P塵也是霧霾污染的一大來(lái)源，以全市房地產(chǎn)開(kāi)發(fā)投資完成額作為房產(chǎn)建設(shè)的替代變量。另一方面，霧霾治理的積極性與受教育水平（EL）有較大聯(lián)系，因?yàn)榻逃降母叩驮谝欢ǔ潭壬嫌绊懼姾推髽I(yè)對(duì)待環(huán)境保護(hù)的重視程度，同樣對(duì)霧霾污染起作用，這里選用每萬(wàn)人在校大學(xué)生人數(shù)作為教育水平的衡量指標(biāo)。

鑒于數(shù)據(jù)的連續(xù)性和可獲得性，本研究樣本時(shí)序?yàn)?005—2016年，涵蓋中國(guó)262個(gè)城市12年的面板數(shù)據(jù)。除前文已作說(shuō)明的中國(guó)城市年均PM2.5濃度具體值，其他數(shù)據(jù)主要源自歷年的《中國(guó)城市統(tǒng)計(jì)年鑒》和《中國(guó)區(qū)域統(tǒng)計(jì)年鑒》，部分缺失數(shù)據(jù)采用線性插值法估算而得，同時(shí)所有變量均作了對(duì)數(shù)化處理，以減輕極端觀測(cè)值對(duì)模型估計(jì)所帶來(lái)的影響。表1顯示了相關(guān)變量的描述性統(tǒng)計(jì)結(jié)果。

表1 相關(guān)變量的描述性統(tǒng)計(jì)

2 實(shí)證結(jié)果分析

2.1 空間自相關(guān)分析

在研究環(huán)境規(guī)制對(duì)城市霧霾污染的空間影響之前，需要先對(duì)其空間相關(guān)性是否顯著存在進(jìn)行檢驗(yàn)，因此下面借助Geoda軟件中的單變量和雙變量LISA模塊進(jìn)行空間自相關(guān)分析。

表2分別給出了2005—2016年城市環(huán)境規(guī)制與霧霾污染的全局空間自相關(guān)檢驗(yàn)下單變量和雙變量全局Moran’s I統(tǒng)計(jì)值。單變量全局Moran’s I值顯示，樣本期內(nèi)城市霧霾污染、正式環(huán)境規(guī)制和非正式環(huán)境規(guī)制的Moran’s I值均至少在5%水平上顯著為正，且霧霾污染的Moran’s I值（均值為0.759）最高，非正式環(huán)境規(guī)制的Moran’s I值（均值為0.458）次之，正式環(huán)境規(guī)制的Moran’s I值（均值為0.151）最低，這表明三者在空間分布上并非處于隨機(jī)狀態(tài)，而是表現(xiàn)出顯著的空間正相關(guān)性，但相關(guān)程度存在一定差異。由于霧霾污染的跨區(qū)域流動(dòng)性比環(huán)境規(guī)制更強(qiáng)，故更容易發(fā)生空間聚集效應(yīng)，使得霧霾污染物呈塊狀分布，空間自相關(guān)性自然更為明顯。進(jìn)一步分析雙變量全局Moran’s I值，發(fā)現(xiàn)樣本期內(nèi)“正式環(huán)境規(guī)制—霧霾污染”的Moran’s I值通過(guò)了1%的顯著性水平檢驗(yàn)，“非正式環(huán)境規(guī)制—霧霾污染”的Moran’s I值除了2016年未通過(guò)顯著性檢驗(yàn)以外其他年份至少在10%水平下顯著，且絕大多數(shù)時(shí)期“正式環(huán)境規(guī)制—霧霾污染”的Moran’s I值大于“非正式環(huán)境規(guī)制—霧霾污染”的Moran’s I值，這說(shuō)明在空間上262個(gè)城市間正式、非正式環(huán)境規(guī)制與霧霾污染的地理分布存在正向依賴(lài)關(guān)系。總的來(lái)看，無(wú)論是正式、非正式環(huán)境規(guī)制和霧霾污染的單變量全局自相關(guān)檢驗(yàn)，還是環(huán)境規(guī)制與霧霾污染之間的雙變量全局自相關(guān)檢驗(yàn)，都印證了地理學(xué)第一定律的正確性，反映出不同區(qū)域間各要素均存在空間關(guān)聯(lián)，因此環(huán)境規(guī)制與霧霾污染的空間關(guān)聯(lián)屬性必須納入到后續(xù)的回歸分析之中。

表2 環(huán)境規(guī)制與霧霾污染的單變量和雙變量全局Moran’s I統(tǒng)計(jì)值

雙變量全局自相關(guān)檢驗(yàn)可能隱藏了局部個(gè)體中也許存在的“非典型”特征，即與全局自相關(guān)檢驗(yàn)不一致的情況，故為了進(jìn)一步顯示正式（非正式）環(huán)境規(guī)制與霧霾污染局部個(gè)體的空間集聚和分異特征，對(duì)其雙變量局部自相關(guān)進(jìn)行檢驗(yàn)。表3列舉了2005、2010和2016年3個(gè)觀測(cè)時(shí)點(diǎn)上，5%顯著性水平下環(huán)境規(guī)制和霧霾污染的雙變量LISA聚類(lèi)結(jié)果。

從表3可知：環(huán)境規(guī)制與城市霧霾污染的聚類(lèi)類(lèi)型以H-H（集聚區(qū)）和L-L（蕭條區(qū)）為主①，雖然部分城市無(wú)顯著的空間相關(guān)性，但在整體上呈現(xiàn)出空間正相關(guān)關(guān)系，驗(yàn)證了前文得到的全局自相關(guān)檢驗(yàn)的結(jié)論。從3個(gè)觀測(cè)時(shí)點(diǎn)來(lái)看，顯著區(qū)域正式環(huán)境規(guī)制與霧霾污染的H-H、H-L（孤島區(qū)）型城市數(shù)量有所增加，L-L、L-H（空心區(qū)）型城市數(shù)量有所減少，但由于H-H型城市數(shù)量遠(yuǎn)高于其他聚類(lèi)類(lèi)型城市數(shù)量，使得H-H和L-L型的正相關(guān)城市數(shù)量之和始終占據(jù)主導(dǎo)，即表現(xiàn)出顯著的空間正相關(guān)性；而顯著區(qū)域非正式環(huán)境規(guī)制與霧霾污染的H-H型城市數(shù)量雖然也有些許上漲，但是L-L型城市數(shù)量下降較多，且H-L和L-H型的負(fù)相關(guān)城市數(shù)量不斷上升，從2005年的正相關(guān)城市數(shù)量占主導(dǎo)到2010、2016年負(fù)相關(guān)城市數(shù)量實(shí)現(xiàn)反超，這也進(jìn)一步說(shuō)明了為什么全局自相關(guān)檢驗(yàn)中2010、2016年的Moran’s I值雖然為正但并不是非常顯著。

2.2 空間杜賓模型估計(jì)分析

通過(guò)對(duì)正式、非正式環(huán)境規(guī)制與霧霾污染的全局以及局部空間自相關(guān)分析，可以基本認(rèn)定正式、非正式環(huán)境規(guī)制對(duì)霧霾污染確實(shí)存在空間影響。關(guān)于具體影響效應(yīng)如何、正式和非正式環(huán)境規(guī)制對(duì)霧霾污染的空間影響是否存在差異等問(wèn)題，下面通過(guò)STATA軟件和空間計(jì)量方法對(duì)式（3）所構(gòu)建的空間杜賓模型進(jìn)行回歸分析。環(huán)境規(guī)制對(duì)霧霾污染空間影響的相關(guān)模型估計(jì)結(jié)果如表4所示。

表3 環(huán)境規(guī)制和霧霾污染的雙變量LISA聚類(lèi)結(jié)果

表4 模型估計(jì)結(jié)果

續(xù)表（表4）

在回歸之前需要先對(duì)模型進(jìn)行初步識(shí)別。從表4中非空間靜態(tài)面板模型的拉格朗日乘數(shù)（LM）和穩(wěn)健的拉格朗日乘數(shù)（robust LM）結(jié)果可以看出：正式和非正式環(huán)境規(guī)制空間滯后的LM、robust LM和空間誤差的LM統(tǒng)計(jì)量均通過(guò)了1%顯著性水平，而兩者空間誤差的robust LM檢驗(yàn)均未通過(guò)顯著性檢驗(yàn)，表明正式和非正式環(huán)境規(guī)制對(duì)霧霾污染都存在顯著的空間影響效應(yīng)，但是空間滯后模型要優(yōu)于空間誤差模型。進(jìn)一步通過(guò)LR和Wald檢驗(yàn)結(jié)果可知，空間杜賓模型同時(shí)優(yōu)于空間滯后和空間誤差模型。因此，后續(xù)基于空間杜賓模型分析環(huán)境規(guī)制對(duì)霧霾污染的空間影響。與此同時(shí)，還選用非空間靜態(tài)面板模型與空間杜賓模型進(jìn)行對(duì)比分析。根據(jù)Hausman統(tǒng)計(jì)量可以確定非空間靜態(tài)面板模型和空間杜賓模型選用固定效應(yīng)較為適宜，再結(jié)合空間杜賓模型的時(shí)間效應(yīng)、空間效應(yīng)的LR統(tǒng)計(jì)量（限于篇幅，未詳細(xì)列出），最終正式環(huán)境規(guī)制選擇時(shí)空雙固定效應(yīng)空間杜賓模型、非正式環(huán)境規(guī)制選擇時(shí)間固定效應(yīng)空間杜賓模型進(jìn)行實(shí)證分析。

從不考慮空間效應(yīng)的非空間靜態(tài)面板模型估計(jì)結(jié)果來(lái)看，正式、非正式環(huán)境規(guī)制均對(duì)霧霾污染產(chǎn)生負(fù)向影響，環(huán)境規(guī)制每提升1%，霧霾污染至少降低0.01%，但正式環(huán)境規(guī)制作用不顯著，非正式環(huán)境規(guī)制則通過(guò)了5%的顯著性水平檢驗(yàn)。總的來(lái)看，環(huán)境規(guī)制政策的制定和實(shí)施抑制了霧霾污染的增加，但正式、非正式環(huán)境規(guī)制效用大小的不協(xié)調(diào)也反映出當(dāng)前異質(zhì)性環(huán)境規(guī)制工具組合還有一定的優(yōu)化空間，合理搭配使用不同環(huán)境規(guī)制工具，對(duì)于有效治霾可以起到事半功倍的效果。

從考慮空間效應(yīng)的空間杜賓模型估計(jì)結(jié)果來(lái)看，正式、非正式環(huán)境規(guī)制的估計(jì)系數(shù)分別為-0.028 6和-0.140 9，且均在1%水平上顯著，意味著正式、非正式環(huán)境規(guī)制確實(shí)會(huì)顯著改善環(huán)境霧霾污染程度。與非空間靜態(tài)面板模型估計(jì)結(jié)果比較可以發(fā)現(xiàn)：通過(guò)把空間因素納入到回歸分析當(dāng)中，將提高環(huán)境規(guī)制的顯著性，這里既包括環(huán)境規(guī)制估計(jì)系數(shù)顯著性的增加（系數(shù)從不顯著或一般顯著變?yōu)榉浅ｏ@著），也包括環(huán)境規(guī)制估計(jì)系數(shù)絕對(duì)值的增加（系數(shù)絕對(duì)值分別增加0.019 6和0.113 1）。同時(shí)在考慮空間因素之后，回歸模型的擬合優(yōu)度明顯增強(qiáng)，進(jìn)一步說(shuō)明運(yùn)用空間計(jì)量模型進(jìn)行回歸分析是合理且必要的。這是因?yàn)檎健⒎钦江h(huán)境規(guī)制存在空間溢出效應(yīng)，一個(gè)地區(qū)實(shí)施環(huán)境規(guī)制在改善本地區(qū)霧霾污染情況的同時(shí)也會(huì)對(duì)鄰近地區(qū)產(chǎn)生反饋效應(yīng)，即在一定程度上減少了其對(duì)鄰近地區(qū)的霧霾溢出，同樣鄰近地區(qū)環(huán)境規(guī)制的實(shí)施也有利于減輕本地區(qū)霧霾污染程度，正是這種環(huán)境規(guī)制空間溢出效應(yīng)的存在導(dǎo)致了“1+1＞2”的治霾效果，即環(huán)境規(guī)制的實(shí)際治霾效用要遠(yuǎn)大于僅考慮本地區(qū)環(huán)境規(guī)制作用的情況。此外，無(wú)論是非空間靜態(tài)面板模型還是空間杜賓模型中，非正式環(huán)境規(guī)制系數(shù)的絕對(duì)值要大于正式環(huán)境規(guī)制系數(shù)，一個(gè)可能的原因是因?yàn)檎江h(huán)境規(guī)制政策存在時(shí)滯性，使得其從制定到實(shí)施需要花費(fèi)相當(dāng)長(zhǎng)的時(shí)間，等到政策真正實(shí)施之際霧霾污染情況可能已經(jīng)發(fā)生變化，事先設(shè)計(jì)的正式環(huán)境規(guī)制也就不再適宜，降低了環(huán)境規(guī)制的時(shí)效性。表4還給出了空間杜賓模型估計(jì)的滯后項(xiàng)空間系數(shù)，正式與非正式環(huán)境規(guī)制的rho值均顯著為正，表明在大氣循環(huán)作用下，城市霧霾污染在空間上呈現(xiàn)正相關(guān)態(tài)勢(shì)。從最一般層面來(lái)說(shuō)，鄰近地區(qū)霧霾污染濃度每提高1%，將引起本地區(qū)霧霾污染濃度至少顯著提高0.787 4%，顯現(xiàn)出“局部俱樂(lè)部”現(xiàn)象，這與前文的分析基本一致。

從控制變量來(lái)看，當(dāng)前的外商直接投資和交通條件會(huì)顯著提高霧霾污染程度，這是因?yàn)楦鞯卣疄榱颂岣弑镜亟?jīng)濟(jì)績(jī)效通常會(huì)加大招商引資力度，而對(duì)于外商投資的質(zhì)量沒(méi)有很好的甄別，這就使得一些高排放、高污染企業(yè)有機(jī)會(huì)轉(zhuǎn)移到環(huán)境規(guī)制強(qiáng)度較低的地區(qū)，自然也就造成該地霧霾污染加重，使經(jīng)濟(jì)績(jī)效在短期內(nèi)得到改善而環(huán)境績(jī)效卻遭受了長(zhǎng)期的毀滅性打擊。交通條件的改善誘發(fā)了汽車(chē)購(gòu)買(mǎi)量的激增，從而間接加劇了大氣污染中霧霾污染物的含量。產(chǎn)業(yè)結(jié)構(gòu)系數(shù)僅在正式環(huán)境規(guī)制占主導(dǎo)時(shí)顯著促進(jìn)了霧霾污染程度，正如祝麗云等［26］的研究指出，我國(guó)效率低但占比卻相當(dāng)高的低端服務(wù)業(yè)嚴(yán)重加劇了霧霾污染程度。教育水平、房產(chǎn)建設(shè)的估計(jì)系數(shù)并未通過(guò)顯著性檢驗(yàn)，說(shuō)明當(dāng)前教育對(duì)于霧霾等環(huán)境保護(hù)方面的宣傳尚不充分，跟其他因素相比，房產(chǎn)建設(shè)也并未達(dá)到顯著擴(kuò)大霧霾污染的程度。此外，比較非空間靜態(tài)面板模型與空間杜賓模型的結(jié)果可知：控制變量在納入空間因素后影響大小、方向、顯著程度等發(fā)生了一系列變化，且由于后者擬合效果更好，這既說(shuō)明了使用空間計(jì)量方法的合理性和常規(guī)面板模型的不再適用，同時(shí)也說(shuō)明在各城市聯(lián)系日益緊密、城市間空間效應(yīng)日趨明顯的今天，任何因素僅考慮本地的影響是遠(yuǎn)遠(yuǎn)不夠的，還要考慮其他城市因素對(duì)本地的影響，否則政策的實(shí)施將難以達(dá)到預(yù)計(jì)效果。

2.3 穩(wěn)健性檢驗(yàn)

在262個(gè)城市中，可能存在個(gè)別地區(qū)的相鄰城市未被納入到樣本中的情況，導(dǎo)致運(yùn)用鄰接權(quán)重矩陣進(jìn)行的實(shí)證分析可能存在一定誤差，下面分別構(gòu)建地理權(quán)重矩陣和經(jīng)濟(jì)權(quán)重矩陣對(duì)上述模型重新進(jìn)行回歸，以檢驗(yàn)原模型估計(jì)結(jié)果是否具有穩(wěn)健性。其中，地理權(quán)重矩陣借鑒了王喬等［27］的思路，利用城市經(jīng)緯度計(jì)算出城市之間的距離，再以2城市間的距離除以所有城市的距離之和即可得到所需地理權(quán)重矩陣的相應(yīng)位置元素。而經(jīng)濟(jì)權(quán)重矩陣是按照林光平等［28］的做法，使其主對(duì)角元素為零，非對(duì)角元素為城市間人均實(shí)際GDP差值的絕對(duì)值的倒數(shù)。以2個(gè)新構(gòu)建的權(quán)重矩陣為基礎(chǔ)進(jìn)行回歸所得到的穩(wěn)健性檢驗(yàn)結(jié)果如表5所示。不難發(fā)現(xiàn)，無(wú)論是基于地理權(quán)重矩陣還是基于經(jīng)濟(jì)權(quán)重矩陣的空間杜賓模型回歸結(jié)果中，正式環(huán)境規(guī)制和非正式環(huán)境規(guī)制對(duì)霧霾污染均表現(xiàn)出顯著的抑制作用，且空間系數(shù)顯著為正，這與原模型估計(jì)結(jié)果相一致，說(shuō)明回歸結(jié)果是穩(wěn)健可靠的。同時(shí)，從模型擬合情況也可以看出：基于鄰接權(quán)重矩陣的空間杜賓模型擬合優(yōu)度更高，說(shuō)明原回歸結(jié)果能夠更好地揭示環(huán)境規(guī)制對(duì)城市霧霾污染的空間影響。

表5 穩(wěn)健性檢驗(yàn)結(jié)果

3 結(jié)論與建議

大面積霧霾污染已逐漸變得常態(tài)化，如何打破霧霾鎖城現(xiàn)狀是新時(shí)代背景下生態(tài)文明建設(shè)的重要內(nèi)容。本文中以城市間普遍存在的空間關(guān)聯(lián)特征為切入點(diǎn)，運(yùn)用2005—2016年中國(guó)262個(gè)城市的面板數(shù)據(jù)，借助空間計(jì)量方法探討環(huán)境規(guī)制對(duì)霧霾污染的空間影響情況。通過(guò)研究可以得出以下結(jié)論：

1）首先對(duì)環(huán)境規(guī)制和霧霾污染進(jìn)行了空間自相關(guān)分析，全局空間自相關(guān)檢驗(yàn)的單變量與雙變量Moran’s I值顯示，霧霾污染與正式、非正式環(huán)境規(guī)制無(wú)論是單個(gè)變量還是兩兩之間均存在顯著的空間正相關(guān)性，局部空間自相關(guān)檢驗(yàn)也得出了類(lèi)似結(jié)論，整體而言正式、非正式環(huán)境規(guī)制對(duì)霧霾污染均存在空間影響，正式環(huán)境規(guī)制與霧霾污染的空間關(guān)聯(lián)度略高于非正式環(huán)境規(guī)制。

2）其次構(gòu)建空間杜賓模型檢驗(yàn)正式、非正式環(huán)境規(guī)制對(duì)霧霾污染的空間影響，發(fā)現(xiàn)城市間霧霾污染本身具有顯著的空間依賴(lài)性，正式、非正式環(huán)境規(guī)制對(duì)霧霾污染均存在顯著的空間溢出效應(yīng)，但是非正式環(huán)境規(guī)制對(duì)霧霾污染治理的效果要高于正式環(huán)境規(guī)制。此外，外商直接投資和交通條件均會(huì)顯著提高霧霾污染程度，產(chǎn)業(yè)結(jié)構(gòu)僅在正式環(huán)境規(guī)制占主導(dǎo)時(shí)會(huì)顯著提高霧霾污染，而教育水平和房產(chǎn)建設(shè)則對(duì)霧霾污染的影響不明顯。

3）最后分別構(gòu)建地理權(quán)重矩陣和經(jīng)濟(jì)權(quán)重矩陣進(jìn)行穩(wěn)健性檢驗(yàn)，結(jié)果表明原模型估計(jì)結(jié)果是穩(wěn)健的。

基于以上結(jié)論，課題組認(rèn)為要完善環(huán)境規(guī)制，促進(jìn)不同環(huán)境規(guī)制主體共同治霾，同時(shí)由于空間因素作用，霧霾污染的區(qū)域性特征日益突出，進(jìn)行區(qū)域協(xié)作成為徹底鏟除霧霾污染的必然選擇。具體建議如下：

1）應(yīng)進(jìn)一步優(yōu)化和豐富環(huán)境規(guī)制工具，對(duì)于正式環(huán)境規(guī)制可以采取更多的綠色激勵(lì)政策而相應(yīng)減少不合理的費(fèi)用型規(guī)制，增加正式規(guī)制過(guò)程的透明性，避免“尋租”行為導(dǎo)致的規(guī)制效用降低。同時(shí)政府應(yīng)引導(dǎo)公眾和環(huán)保NGO參與治霾，充分發(fā)揮非正式環(huán)境規(guī)制的輔助作用。

2）環(huán)境規(guī)制政策的制定應(yīng)因地制宜、精準(zhǔn)治霾，對(duì)于不合時(shí)宜的規(guī)制手段要及時(shí)予以革新，既要通過(guò)設(shè)置大氣污染物特別排放限值、優(yōu)化能源消費(fèi)結(jié)構(gòu)等手段加大對(duì)高霧霾污染地區(qū)的治理，也要提高低霧霾污染地區(qū)的外商投資質(zhì)量和產(chǎn)業(yè)結(jié)構(gòu)高級(jí)化與合理化，避免形成新的高霧霾污染中心。

3）大力推進(jìn)區(qū)域聯(lián)防聯(lián)控，建立統(tǒng)一的跨區(qū)域聯(lián)合預(yù)警和環(huán)境信息共享平臺(tái)，以城市群為單位展開(kāi)試點(diǎn)，逐漸推廣到全國(guó)所有城市，由點(diǎn)及面循序漸進(jìn)，形成全國(guó)各城市政企民共同治霾的良好局面，打贏藍(lán)天保衛(wèi)戰(zhàn)指日可待。