◎ 王家庭 趙 麗 孫 哲 王 璇

我國(guó)區(qū)域城市化與環(huán)境污染關(guān)系的空間計(jì)量研究

◎ 王家庭 趙 麗 孫 哲 王 璇

結(jié)合我國(guó)29個(gè)省區(qū)2000-2010年的面板數(shù)據(jù)，運(yùn)用空間計(jì)量方法實(shí)證研究了區(qū)域城市化與環(huán)境污染的關(guān)系。實(shí)證結(jié)果表明：（1）工業(yè)污染、生活污染與區(qū)域城市化的關(guān)系依舊呈現(xiàn)倒N型；（2）其他影響因素對(duì)環(huán)境污染的影響程度和方向各不相同；（3）目前中國(guó)的大部分省區(qū)均處于隨城市化率增加環(huán)境污染加劇的階段。

區(qū)域城市化 環(huán)境污染 空間面板數(shù)據(jù)模型

一、引 言

隨著經(jīng)濟(jì)的發(fā)展，人民生活水平逐漸提高，城市化進(jìn)程快速發(fā)展。我國(guó)城市化水平自1992年起，以每年1.2%的平均速度增長(zhǎng)。到2010年底，全國(guó)城市化率達(dá)到49.95%的水平。然而隨之而來(lái)的現(xiàn)實(shí)問(wèn)題就是日益嚴(yán)重的環(huán)境污染。城市化過(guò)程中產(chǎn)生了大量的工業(yè)污染物和生活污染物，當(dāng)這些污染物的排放總量超過(guò)了環(huán)境的自凈能力時(shí)，就會(huì)對(duì)生態(tài)環(huán)境造成大范圍的破壞。

從區(qū)域的角度來(lái)看，各地區(qū)經(jīng)濟(jì)發(fā)展過(guò)程不同，城市化進(jìn)程也有所不同，在此過(guò)程中產(chǎn)生的環(huán)境污染問(wèn)題也就不盡相同。東部地區(qū)經(jīng)濟(jì)發(fā)展較早，城市化進(jìn)程較快，也就伴隨著繁雜的環(huán)境污染問(wèn)題，中西部地區(qū)相對(duì)來(lái)說(shuō)，經(jīng)濟(jì)水平較低，城市化還在不斷建設(shè)的過(guò)程中，所產(chǎn)生的環(huán)境問(wèn)題相對(duì)少而簡(jiǎn)單。然而，一個(gè)地區(qū)的污染并不會(huì)只局限在該地區(qū)內(nèi)，而是會(huì)向周邊省區(qū)滲透，造成周邊省區(qū)的環(huán)境質(zhì)量變差。因此，從空間視角研究區(qū)域城市化與環(huán)境污染之間的關(guān)系問(wèn)題，具有十分重要的現(xiàn)實(shí)意義。

目前，國(guó)內(nèi)外大量的學(xué)者已對(duì)城市化過(guò)程中的環(huán)境污染問(wèn)題進(jìn)行了研究。

國(guó)外方面的相關(guān)研究主要有：Howard（1898）[1]嘗試使用合理的規(guī)劃方法，研究城市化和城市生態(tài)環(huán)境之間的協(xié)調(diào)發(fā)展。Meadows（1972）[2]提出了世界城市化前景的“有極限增長(zhǎng)”的預(yù)測(cè)，使得世界各國(guó)開(kāi)始恐慌城市化進(jìn)程中的資源消耗和環(huán)境污染問(wèn)題。Shafik和Bandyopadhyay（1992）[3]研究了不同的污染物和人均國(guó)民生產(chǎn)總值兩個(gè)變量的對(duì)數(shù)形式之間的關(guān)系。Grossman和Krueger（1995）[4]研究了城市空氣污染、水污染和人均國(guó)內(nèi)生產(chǎn)總值之間的關(guān)系。Vernon（2003）[5]則從城市人口的角度，研究了城市人口規(guī)模對(duì)生態(tài)環(huán)境的影響。

國(guó)內(nèi)方面的相關(guān)研究主要有：姜乃力（1999）[6]和史愛(ài)玲等（1999）[7]分析了城市化造成的“熱島效應(yīng)”、“渾濁島效應(yīng)”、“雨島效應(yīng)”和局部環(huán)流等氣候影響。盛學(xué)良等（2001）[8]分析了城市化引起的城市生態(tài)環(huán)境資源問(wèn)題。楊文舉（2002）[9]提出城市化進(jìn)程和生態(tài)環(huán)境建設(shè)具有明顯的階段性、相互作用性和相互聯(lián)系性。劉耀彬等（2008）[10][11]分析了江西省城市化發(fā)展過(guò)程中的環(huán)境污染排放效應(yīng)。盛光耀（2009）[12]認(rèn)為城市化的發(fā)展模式應(yīng)該與資源環(huán)境相適應(yīng)，并考慮環(huán)境政策的影響。

由此可見(jiàn)，國(guó)外相關(guān)研究起步較早，視角廣泛，方法比較系統(tǒng)化。相比而言，由于中國(guó)城市化特殊性強(qiáng)，國(guó)內(nèi)相關(guān)研究發(fā)展起步較晚，但發(fā)展迅速。整體來(lái)看，研究城市化過(guò)程中的環(huán)境污染分析多為定性分析，且未從空間角度來(lái)分析區(qū)域之間的相互影響。

基于此，本文以空間面板數(shù)據(jù)模型為切入點(diǎn)，研究各省區(qū)作為單位個(gè)體的區(qū)域城市化與環(huán)境污染之間的關(guān)系，考慮相鄰區(qū)域間的相互影響。首先利用Moran’s I指數(shù)評(píng)判環(huán)境污染的區(qū)域間的相關(guān)性，之后建立計(jì)量模型，嘗試確定我國(guó)區(qū)域城市化與環(huán)境污染之間的關(guān)系的曲線形狀并量化區(qū)域間的空間相關(guān)性，試圖解答怎樣的環(huán)境污染是合適的、環(huán)境污染的主要影響因素及其影響方向等問(wèn)題。

二、理論分析

城市化水平與環(huán)境污染物的排放，伴隨著時(shí)間的變化呈現(xiàn)出不同的特征，本文將其分為四個(gè)階段，所得曲線如圖1所示。其中細(xì)線表示城市化水平變化曲線，粗線表示環(huán)境污染水平變化曲線。

（一）城市化與環(huán)境污染的初始發(fā)展階段

圖1 城市化水平與環(huán)境污染物排放量的時(shí)期變化圖

在0至T1時(shí)期，城市化水平較低，處于0至10%之間，屬于城市化進(jìn)程的初期階段。在此階段，產(chǎn)業(yè)結(jié)構(gòu)以農(nóng)業(yè)為主，大部分人口集中于農(nóng)村地區(qū)，城市人口低速增加，甚至出現(xiàn)回流現(xiàn)象。聚集的城市的人口數(shù)量相對(duì)較少，產(chǎn)生的生活廢物較少。而有限的收入水平，使得居民的消費(fèi)結(jié)構(gòu)主要以生活必需品為主，對(duì)環(huán)境造成的污染是相當(dāng)有限的。同時(shí)，工業(yè)生產(chǎn)水平低，資源使用較少，投入的有限決定了非期望產(chǎn)出（也就是環(huán)境污染物）的有限。總體來(lái)看，環(huán)境污染物排放處于較低水平。而由于此時(shí)的經(jīng)濟(jì)發(fā)展也屬于初期階段，不穩(wěn)定因素較多，也就造成了環(huán)境污染的波動(dòng)性圖形。在某一時(shí)刻，達(dá)到波動(dòng)的最低點(diǎn)，對(duì)應(yīng)城市化水平為X1*，環(huán)境污染物排放量水平為Y1*。

（二）城市化與環(huán)境污染快速增長(zhǎng)階段

在T1至T2時(shí)期，城市化發(fā)展速度加快，城市化總體水平迅速提高，達(dá)到30%至70%之間。在此過(guò)程中，伴隨著經(jīng)濟(jì)水平的快速增長(zhǎng)。經(jīng)濟(jì)的發(fā)展帶來(lái)技術(shù)進(jìn)步，農(nóng)村勞動(dòng)力生產(chǎn)效率提高，農(nóng)村出現(xiàn)剩余勞動(dòng)力；同時(shí)，城市中的工業(yè)迅速壯大擴(kuò)張，需要更多的勞動(dòng)力投入。于是農(nóng)村人口更多更快地往城市地區(qū)聚集，生活污染逐漸增多。居民生活水平也快速提高，開(kāi)始不僅滿足于必需品消費(fèi)，傾向于改善衣、住、行的消費(fèi)，而且更多的是一種數(shù)量的提高。這種消費(fèi)結(jié)構(gòu)促進(jìn)了相關(guān)產(chǎn)品的生產(chǎn)，從而產(chǎn)生較多的工業(yè)污染物。而社會(huì)總生產(chǎn)也進(jìn)入高速化階段，投入生產(chǎn)的資源急速增加，在大部分地區(qū)采用的是耗能型技術(shù)的前提下，環(huán)境污染，無(wú)論是工業(yè)污染，還是生活污染，都進(jìn)入迅速攀升階段。

（三）城市化緩慢增長(zhǎng)，環(huán)境污染達(dá)到極值階段

在T2至T3時(shí)期，城市化進(jìn)程逐漸進(jìn)入穩(wěn)定階段，城市化水平達(dá)到75%以上，經(jīng)濟(jì)也相應(yīng)進(jìn)入穩(wěn)速發(fā)展階段。城市人口增長(zhǎng)速度減緩，城市居民素質(zhì)提高，開(kāi)始注重生活質(zhì)量的提高，傾向于節(jié)能環(huán)保型產(chǎn)品的消費(fèi)，生活廢物排放傾向于減少。同時(shí)，這種消費(fèi)結(jié)構(gòu)的改變相應(yīng)地改變了生產(chǎn)結(jié)構(gòu)，工業(yè)環(huán)境污染物排放量減少。此時(shí)的生產(chǎn)技術(shù)也趨于成熟，更多地區(qū)的技術(shù)傾向于節(jié)能環(huán)保型，但由于環(huán)境污染具有一定的粘性，排放量會(huì)繼續(xù)增長(zhǎng)，直至達(dá)到一點(diǎn)，經(jīng)濟(jì)發(fā)展帶來(lái)的技術(shù)效應(yīng)大于人口、資源在空間上集聚產(chǎn)生的規(guī)模效應(yīng)，污染物排放開(kāi)始減少，這一點(diǎn)也就是污染物排放最多的臨界點(diǎn)，最高值點(diǎn)。此時(shí)對(duì)應(yīng)的城市化水平為X2*，環(huán)境污染物排放量水平為Y2*。

（四）城市化均衡，環(huán)境污染緩慢減少階段

在T3時(shí)期之后，進(jìn)入城市化進(jìn)程的均衡階段，城市化水平在此階段達(dá)到80%到90%之間，經(jīng)濟(jì)發(fā)展也相應(yīng)地進(jìn)入均衡狀態(tài)。人口和經(jīng)濟(jì)不再向城市聚集，農(nóng)村和城市的發(fā)展不再具有足夠大的差異性，使得資源出現(xiàn)大量的單向轉(zhuǎn)移。城市人口比率、生產(chǎn)水平、資源投入以及技術(shù)水平都相應(yīng)進(jìn)入穩(wěn)態(tài)。此時(shí)的污染物排放跨過(guò)臨界值，進(jìn)入緩慢減少階段。

將城市化水平與環(huán)境污染分別作為X軸和Y軸，建立直角坐標(biāo)系，以更直接地表示二者關(guān)系，所得結(jié)果如圖2所示。

三、指標(biāo)選取及相關(guān)說(shuō)明

（一） 被解釋變量的選取

圖2 城市化與環(huán)境污染的關(guān)系圖

從現(xiàn)實(shí)情況看，工業(yè)污染和生活污染是環(huán)境污染的主要組成部分，其影響因素會(huì)有所不同，于是本文將環(huán)境污染分為工業(yè)污染與生活污染兩部分進(jìn)行細(xì)化研究。

1. 工業(yè)污染的變量選取。選取工業(yè)廢水排放量（萬(wàn)噸）、工業(yè)廢氣排放量（萬(wàn)噸）、工業(yè)二氧化硫排放量（萬(wàn)噸）、工業(yè)煙塵排放量（萬(wàn)噸）、工業(yè)固體廢棄物產(chǎn)生量（萬(wàn)噸）五個(gè)工業(yè)污染物相關(guān)指標(biāo)。

2. 生活污染的變量選取。選取生活二氧化硫排放量（萬(wàn)噸）、生活煙塵排放量（萬(wàn)噸）、生活污水排放量（萬(wàn)噸）、生活垃圾清運(yùn)量（萬(wàn)噸）四個(gè)相關(guān)指標(biāo)。

同時(shí)，對(duì)工業(yè)污染物指標(biāo)體系和生活污染物指標(biāo)體系通過(guò)熵值法計(jì)算降維，分別得到工業(yè)污染和生活污染綜合指數(shù)，定義為Y1和Y2，以保證較為全面地反映環(huán)境污染水平。

（二）解釋變量的選取

1.工業(yè)污染回歸方程解釋變量的選取

蔡孝箴（1998）[13]認(rèn)為，城市化過(guò)程除了人口城市化為主要代表外，還應(yīng)包括產(chǎn)業(yè)城市化、生活方式城市化和地域城市化三個(gè)方面的考慮。我們以控制變量形式將這三種城市化加入模型之中。此外，從區(qū)域自身發(fā)展條件和外在發(fā)展條件考慮，加入?yún)^(qū)域技術(shù)水平指標(biāo)和貿(mào)易開(kāi)放程度指標(biāo)。

（1）選取各省區(qū)城鎮(zhèn)人口比重（%）作為本文中的主要解釋變量——人口城市化指標(biāo)，也就是普遍意義的城市化，在回歸模型中定義為X1。

（2）選取各省區(qū)工業(yè)產(chǎn)值占地區(qū)生產(chǎn)總值的比率（%）①作為產(chǎn)業(yè)城市化指標(biāo)，定義為X2。工業(yè)產(chǎn)值在生產(chǎn)總值中的比重，也就是工業(yè)化率，對(duì)環(huán)境污染起正向影響作用。

（3）選取各省區(qū)非生活必需品消費(fèi)性支出占城鎮(zhèn)居民家庭平均每人全年消費(fèi)性支出的比值（%）作為生活方式城市化指標(biāo)，定義為X3。居民用于非生活必需品的消費(fèi)，初期對(duì)環(huán)境污染起加速作用。而伴隨著經(jīng)濟(jì)的發(fā)展，居民收入水平進(jìn)一步提高，開(kāi)始注重生活環(huán)境質(zhì)量的提高，傾向于購(gòu)買(mǎi)節(jié)能環(huán)保型產(chǎn)品，減少了環(huán)境污染物的排放。

（4）選取人口密度（人/平方公里）作為地域城市化的指標(biāo)，定義為X4。人口密度大的地區(qū)也就意味著更多的擴(kuò)大再生產(chǎn)，更多的消費(fèi)，更多的日常生活活動(dòng)，也就伴隨著更多的污染物排放。

（5）選取各地區(qū)用于科學(xué)技術(shù)的財(cái)政支出在財(cái)政總支出中所占比重（%）作為區(qū)域內(nèi)部技術(shù)水平指標(biāo)，定義為X5。技術(shù)進(jìn)步會(huì)從各方面降低污染物的排放。

（6）選取各地區(qū)進(jìn)出口總額與地方生產(chǎn)總值之比作為貿(mào)易開(kāi)放程度指標(biāo)，定義為X6。依據(jù)比較優(yōu)勢(shì)理論，國(guó)際貿(mào)易使得環(huán)境污染嚴(yán)重的產(chǎn)業(yè)由發(fā)達(dá)國(guó)家移入發(fā)展中國(guó)家，或是由環(huán)境管制嚴(yán)格的國(guó)家轉(zhuǎn)移到管制寬松的國(guó)家，以此規(guī)避環(huán)境污染約束。國(guó)際貿(mào)易使得環(huán)境污染產(chǎn)生“替代效應(yīng)”，也就是“污染天堂”假說(shuō)。我們可以認(rèn)為，貿(mào)易開(kāi)放度越高，本省區(qū)的污染越嚴(yán)重。

2. 生活污染方程解釋變量的選取

造成生活污染的影響因素，我們可以從和居民生活息息相關(guān)的衣食住行等方面進(jìn)行考慮。在工業(yè)污染方程中，涉及到的四個(gè)城市化的指標(biāo)中，人口城市化、生活方式城市化和地域城市化指標(biāo)也會(huì)對(duì)生活造成污染，而產(chǎn)業(yè)城市化指標(biāo)——工業(yè)化率通過(guò)影響生活垃圾的排放影響生活污染指數(shù)的變化。由此我們繼續(xù)沿用這四個(gè)指標(biāo)。

除此之外，根據(jù)數(shù)據(jù)的可行性和可得性，考慮另外兩個(gè)因素的影響。

（1）選取私人汽車(chē)擁有量（萬(wàn)輛），定義為X7。私人汽車(chē)擁有量越多，對(duì)大氣環(huán)境造成的污染越嚴(yán)重。

（2）選取餐飲業(yè)企業(yè)營(yíng)業(yè)額（億元），定義為X8。餐飲業(yè)企業(yè)營(yíng)業(yè)額，反映了居民外出就餐的消費(fèi)程度，而餐飲業(yè)是造成生活垃圾和污水排放的主要因素。營(yíng)業(yè)額越多，表明食物垃圾和廢水排放越嚴(yán)重。

（三）相關(guān)說(shuō)明

本文選取時(shí)間范圍為2000年至2010年，選取截面為全國(guó)29個(gè)省區(qū)（西藏和海南地區(qū)除外）②。2000年至2003年的城市化率指標(biāo)來(lái)源于對(duì)應(yīng)年份的《中國(guó)全面建設(shè)小康社會(huì)進(jìn)程統(tǒng)計(jì)監(jiān)測(cè)報(bào)告》，2003年與2004年單位地區(qū)生產(chǎn)總值能耗來(lái)源于《中國(guó)能源統(tǒng)計(jì)年鑒》相關(guān)數(shù)據(jù)計(jì)算得到，其他指標(biāo)均來(lái)源于對(duì)應(yīng)年份的《中國(guó)統(tǒng)計(jì)年鑒》。

四、實(shí)證分析

（一） 基于Moran’s I指數(shù)的空間相關(guān)性的分析

1. Moran’s I 指數(shù)理論簡(jiǎn)要介紹

Moran’s I，又稱(chēng)為空間自相關(guān)指數(shù)，取值范圍為－1≤I≤1。I＜0表示負(fù)向的空間相關(guān)性，說(shuō)明不同的屬性值傾向于集聚；I＞0表示正向的空間相關(guān)性，說(shuō)明相似的屬性值傾向于集聚；I=0表示屬性值的分布與區(qū)位數(shù)據(jù)的分布相互獨(dú)立，不存在空間相關(guān)性。I的絕對(duì)值越大，表示空間的聚集度越高。

如果需要探討局部空間集聚的存在與否，確定每個(gè)區(qū)域?qū)θ值目臻g相關(guān)性的貢獻(xiàn)大小，空間自相關(guān)的全局評(píng)估掩蓋了多大程度上的局部不穩(wěn)定性，就要運(yùn)用局部空間自相關(guān)分析。本文選用Moran’s I散點(diǎn)圖進(jìn)行分析。Moran’s I散點(diǎn)圖分為四個(gè)象限，對(duì)應(yīng)相鄰單元之間的四種局部空間相關(guān)性形式：第一象限對(duì)應(yīng)高觀測(cè)值的區(qū)域單元被同樣的高觀測(cè)值的區(qū)域單元包圍的形式（簡(jiǎn)稱(chēng)高-高，H-H）；第二象限表示低觀測(cè)值的區(qū)域單元被高觀測(cè)值的區(qū)域包圍（簡(jiǎn)稱(chēng)低-高，L-H）；第三象限表示低觀測(cè)值的區(qū)域單元被低觀測(cè)值的區(qū)域單元所包圍（簡(jiǎn)稱(chēng)低-低，L-L）；第四象限表示高觀測(cè)值的區(qū)域單元被低觀測(cè)值的區(qū)域所包圍（簡(jiǎn)稱(chēng)高-低，H-L）。

2. 基于Moran’s I指數(shù)值的環(huán)境污染全局空間相關(guān)性分析

由于Moran’s I指數(shù)值反映截面數(shù)據(jù)的性質(zhì)，而本文涉及兩個(gè)環(huán)境污染綜合指數(shù)11年的數(shù)據(jù)，比較龐雜。嘗試對(duì)每個(gè)環(huán)境污染綜合指數(shù)的11年數(shù)據(jù)進(jìn)行平均化處理，形成截面數(shù)據(jù)，然后計(jì)算Moran’s I指數(shù)，并進(jìn)行分析，結(jié)果如表1所示。

表1 我國(guó)各省區(qū)的區(qū)域環(huán)境污染水平的全局Moran’s I指數(shù)

表1結(jié)果顯示，工業(yè)和生活環(huán)境污染綜合指數(shù)的Moran’s I指數(shù)的數(shù)值都大于0，說(shuō)明我國(guó)的區(qū)域環(huán)境污染從整體上來(lái)看，并不是每個(gè)省區(qū)隨機(jī)分布并獨(dú)立存在的，而是存在著明確的正向空間依賴關(guān)系。環(huán)境污染嚴(yán)重的區(qū)域傾向于集聚在一起，一個(gè)區(qū)域的環(huán)境污染水平受到相鄰區(qū)域的同向影響。工業(yè)污染綜合指數(shù)的Moran’s I的數(shù)值較大，而生活綜合污染指數(shù)的Moran’s I的數(shù)值較小，說(shuō)明工業(yè)污染的空間相關(guān)性強(qiáng)于生活污染的空間相關(guān)性。

3. 基于Moran’s I散點(diǎn)圖的環(huán)境污染局部空間相關(guān)性分析

為了更進(jìn)一步判定我國(guó)各省區(qū)環(huán)境污染的高低屬性，進(jìn)行局部空間相關(guān)性分析。

（1）工業(yè)污染局部空間相關(guān)性分析

由圖3的散點(diǎn)圖（左圖）可以看到，大多數(shù)省區(qū)集中于第一象限和第三象限，也就是工業(yè)環(huán)境污染集聚的高-高和低-低類(lèi)型占主導(dǎo)。具體省區(qū)的分布情況如表2所示。位于第一象限（高-高類(lèi)型）的省區(qū)主要有黑龍江、遼寧（東北老工業(yè)基地省區(qū)）；山西、河南、湖南、湖北、安徽、江西（中部工業(yè)省區(qū)）；河北、山東、上海、江蘇、浙江、福建、廣東（東部沿海省區(qū)）。東北老工業(yè)基地省區(qū)長(zhǎng)久以來(lái)的產(chǎn)業(yè)基礎(chǔ)都是工業(yè)；而中部工業(yè)省區(qū)有的憑借資源優(yōu)勢(shì)、有的依靠東部地區(qū)的產(chǎn)業(yè)轉(zhuǎn)移，近些年來(lái)工業(yè)經(jīng)濟(jì)崛起；而涉及到的東部沿海省區(qū)都是輕紡工業(yè)發(fā)達(dá)的地區(qū)。以工業(yè)為主的產(chǎn)業(yè)結(jié)構(gòu)決定了這些省區(qū)高污染的地位。處于第三象限（低-低類(lèi)型）的省區(qū)為新疆、甘肅、寧夏、青海、云南、廣西，這些都是經(jīng)濟(jì)不發(fā)達(dá)的西部省區(qū)。位于第二象限（低-高類(lèi)型）的省區(qū)多為東部省區(qū)的腹地，本身經(jīng)濟(jì)發(fā)展水平不高，污染程度較低，但是有受東部地區(qū)重污染影響的可能性。

圖3 工業(yè)污染綜合指數(shù)的Moran’s I散點(diǎn)圖及顯著性水平圖

從圖3的顯著性水平圖（右圖）來(lái)看，在工業(yè)污染集聚的低-低類(lèi)型中，新疆和青海通過(guò)了顯著性水平檢驗(yàn)，也就是這種特征最為明顯的省區(qū)。這和其經(jīng)濟(jì)發(fā)展水平密切相關(guān)，與其他省區(qū)相比，這兩個(gè)城市的工業(yè)化水平更低，工業(yè)污染排放較少，而周邊省區(qū)也大多為工業(yè)發(fā)展水平不高的省區(qū)，也就形成了以這兩個(gè)省區(qū)為核心的低工業(yè)污染圈。而在高-高類(lèi)型中，通過(guò)了顯著性水平的省區(qū)為河北、山西、山東、河南、江蘇、安徽、上海。這七個(gè)省區(qū)中，有重工業(yè)省區(qū)，也有輕工業(yè)省區(qū)，這就決定了其高污染區(qū)域核心的地位，與周邊省區(qū)形成高污染圈，也就不足為奇了。而四川成為工業(yè)污染高-低特性的核心，也就是說(shuō)四川相對(duì)于周邊省區(qū)來(lái)說(shuō)，工業(yè)化程度較高，吸引著周邊省區(qū)資源向四川的集聚，同時(shí)也產(chǎn)生了污染排放的集聚，形成了“高-低”特性。

（2）生活污染局部空間相關(guān)性分析

由圖4的散點(diǎn)圖（左圖）可以看出，多數(shù)省區(qū)依然位于第一象限和第三象限，也就是生活環(huán)境污染集聚的高-高和低-低類(lèi)型占主導(dǎo)。具體來(lái)看，如表3的區(qū)域象限分布，位于第一象限（高-高類(lèi)型）的省區(qū)中東部省區(qū)有6個(gè)，中部省區(qū)5個(gè)，西部省區(qū)1個(gè)，東北部省區(qū)2個(gè)。地區(qū)的生活污染程度，一方面決定于該區(qū)域的經(jīng)濟(jì)發(fā)展水平。經(jīng)濟(jì)水平越高，人民的可支配收入越高，用于消費(fèi)的部分也就越多，相應(yīng)的就會(huì)造成更多的環(huán)境污染，比如上海、江蘇、廣東等。另一方面決定于當(dāng)?shù)氐奈幕諊m然經(jīng)濟(jì)水平不一定發(fā)達(dá)，但是當(dāng)?shù)厝嗣駜A向于更多的生活消費(fèi)，也會(huì)造成較多的生活污染，比如湖南、河南等。處于第三象限（低-低類(lèi)型）的省區(qū)為新疆、甘肅、青海，都是經(jīng)濟(jì)不發(fā)達(dá)的西部省區(qū)，當(dāng)?shù)鼐用窬窒抻跐M足溫飽的程度，不會(huì)產(chǎn)生很多的生活污染。位于第二象限（低-高類(lèi)型）的省區(qū)可分為兩類(lèi)。一種是北京、天津的經(jīng)濟(jì)發(fā)達(dá)地區(qū)，這些地域由于經(jīng)濟(jì)發(fā)展水平較高，當(dāng)?shù)卣约熬用褚呀?jīng)開(kāi)始重視污染問(wèn)題，產(chǎn)生了環(huán)保意識(shí)，人為地減少生活污染的排放。另一種是經(jīng)濟(jì)發(fā)展水平一般，且當(dāng)?shù)鼐用癫蛔⒅厣钯|(zhì)量的地區(qū)，但這些地區(qū)都有受到周邊高生活污染省區(qū)影響的可能性。

表2 我國(guó)各省區(qū)工業(yè)污染綜合指數(shù)的空間相關(guān)模式

圖4 生活污染綜合指數(shù)的Moran’s I散點(diǎn)圖及顯著性水平圖

由圖4的顯著性水平圖（右圖），可以看出，在生活污染集聚的低-低類(lèi)型中，新疆和青海通過(guò)了顯著性水平檢驗(yàn)，成為特征核心省區(qū)。一方面與這些省區(qū)的經(jīng)濟(jì)發(fā)展水平相關(guān)，另一方面也與當(dāng)?shù)氐娘嬍沉?xí)慣相關(guān)，傾向于產(chǎn)生較少的生活污染。而周邊省區(qū)也大多為生活污染排放較少的省區(qū)，也就形成了以這兩個(gè)省區(qū)為核心的低生活污染圈。而在高-高類(lèi)型中，通過(guò)顯著性水平檢驗(yàn)的省區(qū)為山東、江蘇、安徽、上海、江西、福建。也就是說(shuō)，以這六個(gè)省區(qū)為核心，形成了高污染區(qū)域。而四川成為生活污染高-低特性的核心，也就是說(shuō)四川相對(duì)于周邊省區(qū)來(lái)說(shuō)，有較多的生活污染的排放，這與當(dāng)?shù)氐娘嬍沉?xí)慣是密切相關(guān)的，形成了“高-低”特性。

（二）傳統(tǒng)固定效應(yīng)模型的估計(jì)與分析

利用傳統(tǒng)固定效應(yīng)模型——LSDV（Least Square Dummy Variables），也就是普通面板數(shù)據(jù)模型，進(jìn)行回歸分析并檢驗(yàn)空間相關(guān)性的存在與否。

設(shè)定本文基本模型為：

回歸結(jié)果如表4所示。

分析表4的數(shù)據(jù)，可以看出：

工業(yè)污染的方程中，人口城市化指標(biāo)的一次方、二次方和三次方的系數(shù)通過(guò)了1%的顯著性水平檢驗(yàn)，且與工業(yè)環(huán)境污染物的關(guān)系，呈現(xiàn)倒N型。而其他五個(gè)控制變量，四個(gè)通過(guò)了不同程度的顯著性檢驗(yàn)。

生活污染的方程中，人口城市化一次方、二次方和三次方的系數(shù)均未通過(guò)顯著性檢驗(yàn)，控制變量中三個(gè)通過(guò)了不同程度的顯著性檢驗(yàn)，我們依然能夠觀察到生活污染伴隨區(qū)域城市化發(fā)展的變化曲線呈現(xiàn)倒N型趨勢(shì)。

表3 我國(guó)各省區(qū)生活污染綜合指數(shù)的空間相關(guān)模式

表4 基于LSDV的估計(jì)結(jié)果和空間相關(guān)性檢驗(yàn)

LMerr統(tǒng)計(jì)量用于檢驗(yàn)空間誤差模型中的空間自相關(guān)系數(shù)的顯著性，LMlag統(tǒng)計(jì)量用于檢驗(yàn)空間滯后模型中的空間自回歸系數(shù)的顯著性。LMerr和LMlag統(tǒng)計(jì)量都是通過(guò)截面數(shù)據(jù)的最小二乘估計(jì)（OLS估計(jì)）得到的。Anselin et al（2008）[14]指出，將空間權(quán)重矩陣做相應(yīng)的變化后，這兩個(gè)統(tǒng)計(jì)量能夠繼續(xù)應(yīng)用于面板數(shù)據(jù)模型中兩個(gè)模型使用的判斷?？臻g誤差模型或空間滯后模型的選擇方法為：如果LMerr統(tǒng)計(jì)量的顯著性優(yōu)于LMlag統(tǒng)計(jì)量，則選取模型SEM；如果LMlag優(yōu)于LMerr，則選取SLM。

由表4可以看到，在兩個(gè)方程中，LMlag、Robust LMlag、LMerr和Robust LMerr統(tǒng)計(jì)量均顯著不為0，且LMlag統(tǒng)計(jì)值通過(guò)了顯著性檢驗(yàn)，也就說(shuō)明工業(yè)環(huán)境污染和生活環(huán)境污染均存在著明顯的空間自相關(guān)性，也就是工業(yè)環(huán)境污染和生活環(huán)境污染都會(huì)有臨近區(qū)域間的相互影響。再根據(jù)空間滯后模型與空間誤差模型的選擇標(biāo)準(zhǔn)進(jìn)行判斷，兩個(gè)方程中，LMerr和Robust LMerr統(tǒng)計(jì)量均未通過(guò)顯著性檢驗(yàn)，而LMlag和Robust LMlag統(tǒng)計(jì)量全通過(guò)了顯著性檢驗(yàn)，而且，LMlag和Robust LMlag統(tǒng)計(jì)量大于對(duì)應(yīng)的LMerr和Robust LMerr統(tǒng)計(jì)量，因此工業(yè)污染方程和生活污染方程均選擇空間滯后模型進(jìn)行回歸，測(cè)度相鄰區(qū)域間的相互影響程度。

（三）空間面板數(shù)據(jù)模型的估計(jì)與分析

1. 模型介紹

Anselin（1998）[15]和Elhorst（2003）[16]在傳統(tǒng)的面板數(shù)據(jù)模型的基礎(chǔ)上，明確引入考慮空間相關(guān)性的空間滯后誤差項(xiàng)和空間滯后因變量。Anselin（1998）提出了隨機(jī)效應(yīng)的空間誤差模型。Elhorst（2003）把空間面板數(shù)據(jù)模型分成了四個(gè)小類(lèi)：空間固定效應(yīng)模型、空間固定系數(shù)模型、空間隨機(jī)效應(yīng)模型、空間隨機(jī)系數(shù)模型，給出了四個(gè)模型的對(duì)數(shù)似然函數(shù)，還分析了極大似然估計(jì)量的漸進(jìn)性質(zhì)。

空間計(jì)量模型包括空間誤差模型（SEM，又稱(chēng)為空間自相關(guān)模型）和空間滯后模型（SLM，又稱(chēng)空間自回歸模型）。何江和張馨之（2006）[17]提出，回歸目標(biāo)集中于研究某些特定個(gè)體時(shí)，選擇固定效應(yīng)模型更好；何江和張馨之（2006），王火根和沈利生（2007）[18]，吳擁政和顏日初（2009）[19]，這三組學(xué)者的研究都表明，空間面板研究文獻(xiàn)一般均選擇固定效應(yīng)模型。基于前人的研究成果，本文選擇空間固定效應(yīng)模型進(jìn)行回歸分析。

空間誤差固定效應(yīng)模型（Fixed Effects Spatial Error Model）的設(shè)定如下：

空間滯后固定效應(yīng)模型（Fixed Effects Spatial Lag Model）的設(shè)定如下：

上述兩個(gè)模型的矩陣匯總表示形式如下：

2. 回歸結(jié)果分析

對(duì)兩個(gè)方程的空間滯后固定效應(yīng)模型進(jìn)行回歸，結(jié)果如表5所示。

由表5可以看出，加入空間相關(guān)性的考慮之后，兩個(gè)模型的可決系數(shù)都有一定程度的提高，這說(shuō)明考慮空間相關(guān)性在一定程度上提高了模型的擬合優(yōu)度。Y1方程中，除了城市化率指標(biāo)的系數(shù)依舊通過(guò)顯著性檢驗(yàn)之外，其他控制變量的系數(shù)也通過(guò)了一定水平的顯著性檢驗(yàn)。Y2方程中，城市化指標(biāo)的一次方項(xiàng)和三次方項(xiàng)的系數(shù)通過(guò)10%的顯著性水平檢驗(yàn)，控制變量中通過(guò)顯著性檢驗(yàn)的指標(biāo)也增多。變量系數(shù)顯著性提高，也從一定程度上說(shuō)明采用空間滯后固定模型的正確性?？傮w來(lái)看，兩個(gè)方程中污染物與城市化指標(biāo)的關(guān)系，依舊全部呈現(xiàn)倒N型，這與LSDV估計(jì)結(jié)果相同。從理論來(lái)上講，倒N型曲線是符合環(huán)境發(fā)展特征的，只是與倒U型曲線比較來(lái)說(shuō)，在城市化初期有較大的波動(dòng)。

表5 工業(yè)污染和生活污染的空間固定效應(yīng)模型估計(jì)結(jié)果

兩種污染物的模型中，自回歸系數(shù) 的估計(jì)值均為正，并且都通過(guò)1%或5%水平下的顯著性檢驗(yàn)，說(shuō)明各省區(qū)間的工業(yè)污染物排放和生活污染物排放確實(shí)具有較強(qiáng)的正向相互影響作用。為了更好的解釋環(huán)境污染的影響因素，我們更詳細(xì)地分析每一個(gè)模型方程：

（1）工業(yè)污染方程

第一個(gè)方程中，工業(yè)污染綜合指數(shù)作為被解釋變量，與區(qū)域城市化指標(biāo)的呈現(xiàn)的倒N型曲線中，兩個(gè)轉(zhuǎn)折點(diǎn)處對(duì)應(yīng)的區(qū)域城市化率為X1*=20.63%和X2*=80.70%。

具體來(lái)說(shuō)，區(qū)域城市化率低于20.63%的省市，還處于城市發(fā)展的初期階段；而區(qū)域城市化率介于20.63%和80.70%之間，表明各省區(qū)的工業(yè)環(huán)境污染隨著區(qū)域城市化加劇而惡化；區(qū)域城市化率超過(guò)80.70%的省區(qū)，則進(jìn)入工業(yè)環(huán)境污染隨區(qū)域城市化深入而緩和的階段。

根據(jù)中國(guó)2010年各省區(qū)的城市化率指標(biāo)，各省市在工業(yè)污染與城市化率呈現(xiàn)的倒N型關(guān)系中的所處階段的分布如圖5所示。粗實(shí)線表示區(qū)域城市化率為20.63%的基準(zhǔn)線，對(duì)應(yīng)的為極小值處的X1*的區(qū)域城市化水平；細(xì)虛線表示區(qū)域城市化率為80.70%的基準(zhǔn)線，對(duì)應(yīng)的是極大值處的X2*的區(qū)域城市化水平。除了北京在2010年的人口城市化率達(dá)到85.93%，上海達(dá)到88.86%，超過(guò)了第二個(gè)轉(zhuǎn)折點(diǎn)，進(jìn)入工業(yè)污染隨著人口城市化率的增加而減少的階段，其他省區(qū)依然處于工業(yè)廢水排放量隨人口城市化率增加而增加的階段。

由此可以看出，北京和上海兩個(gè)直轄市已經(jīng)跨過(guò)以污染換增長(zhǎng)的階段，經(jīng)濟(jì)發(fā)展進(jìn)入平穩(wěn)發(fā)展階段，更具有長(zhǎng)期持續(xù)發(fā)展的特性，成為“經(jīng)濟(jì)發(fā)展與環(huán)保兩手抓”的代表。

對(duì)于其他控制變量：伴隨著工業(yè)產(chǎn)值比率增加，各省區(qū)的工業(yè)污染加劇，這是與事實(shí)相符合的；伴隨著非生活必需品消費(fèi)性支出比率增加，工業(yè)污染加劇，這說(shuō)明目前居民生活水平的提高很大程度上表現(xiàn)于家用電器等生活用品的量的增加，從而促進(jìn)生產(chǎn)的增多，資源的消耗；伴隨著人口密度的增加，工業(yè)污染在減少；伴隨著區(qū)域內(nèi)部技術(shù)的進(jìn)步，工業(yè)污染也在減少；伴隨著貿(mào)易開(kāi)放度的增加，工業(yè)污染在減少。

表示工業(yè)污染的空間相關(guān)性的自回歸系數(shù) 的估計(jì)值為正，并且通過(guò)1%水平下的顯著性檢驗(yàn)，說(shuō)明各省區(qū)間的工業(yè)污染物排放具有較強(qiáng)的正向空間相關(guān)性。也就是說(shuō)，本省區(qū)的工業(yè)污染物排放受相鄰地區(qū)工業(yè)污染排放的影響，且影響方向?yàn)橥?。自回歸系數(shù)的估計(jì)值為0.3941，表明相鄰省區(qū)的工業(yè)污染綜合指數(shù)變化1%，則會(huì)導(dǎo)致本省區(qū)工業(yè)污染綜合指數(shù)同向變動(dòng)0.3941%。

圖5 基于工業(yè)污染的2010年區(qū)域城市化率分布情況

（2）生活污染方程

在第二個(gè)方程中，區(qū)域生活污染綜合指數(shù)作為被解釋變量，與區(qū)域城市化率呈現(xiàn)的倒N型曲線中，轉(zhuǎn)折點(diǎn)依次為19.10%和74.32%。

參照中國(guó)2010年各省區(qū)的人口城市化率，各省區(qū)在生活污染和城市化率呈現(xiàn)的倒N型關(guān)系中的所處階段的具體分布如圖6所示。粗實(shí)線表示區(qū)域城市化率為19.10%的基準(zhǔn)線，細(xì)虛線表示區(qū)域城市化率為74.32%的基準(zhǔn)線。

具體來(lái)看，北京（城市化率為85.90%）、天津（城市化率為79.20%）和上海（城市化率為88.86%）三個(gè)直轄市，已越過(guò)第二個(gè)轉(zhuǎn)折點(diǎn)，進(jìn)入?yún)^(qū)域生活污染排放隨著區(qū)域城市化率增加而減少的階段，大多數(shù)省區(qū)仍處在兩者正相關(guān)的部分。直轄市在生活污染治理中的標(biāo)桿作用在這里可以看到，在注重經(jīng)濟(jì)建設(shè)的同時(shí)，關(guān)注精神文明的建設(shè)，居民環(huán)保觀念較強(qiáng)。

在其他控制變量中，產(chǎn)業(yè)結(jié)構(gòu)的系數(shù)未通過(guò)顯著性檢驗(yàn)，說(shuō)明從現(xiàn)實(shí)來(lái)講，生活污染與工業(yè)化的關(guān)系較?。浑S著非生活必需品消費(fèi)性支出比率的增加，生活污染增加；隨著人口密度的增加，生活污染減少，說(shuō)明除了少數(shù)省區(qū)人口密度過(guò)度集中之外，人口密度還未對(duì)居民生活造成負(fù)面影響，而是規(guī)模效應(yīng)依然占著主導(dǎo)地位；隨著私人汽車(chē)擁有量的增多，生活污染在增加，這是與現(xiàn)實(shí)相符合的；隨著餐飲業(yè)營(yíng)業(yè)額的增加，生活污染也在加劇。

表示生活污染的空間相關(guān)性的自回歸系數(shù) 的估計(jì)值也為正，并且通過(guò)5%水平下的顯著性檢驗(yàn)，說(shuō)明各省區(qū)間的生活污染物排放具有較強(qiáng)的正向空間相關(guān)性。也就是說(shuō)，本省區(qū)的生活污染物排放受相鄰地區(qū)生活污染排放的影響，且影響方向?yàn)橥?。自回歸系數(shù)的估計(jì)值為0.2124，表明相鄰省區(qū)的生活污染綜合指數(shù)變化1%，則會(huì)導(dǎo)致本省區(qū)生活污染綜合指數(shù)同向變動(dòng)0.2124%。

圖6 基于生活污染的2010年區(qū)域城市化率分布情況

五、結(jié)論及政策建議

通過(guò)實(shí)證分析，我們得到以下結(jié)論：（1）全國(guó)29個(gè)省區(qū)的環(huán)境污染大多為同屬性聚集。工業(yè)污染和生活污染的空間相關(guān)性多為“高-高”、“低-低”聚集類(lèi)型；無(wú)論是工業(yè)污染還是生活污染，其“高-高”空間集聚類(lèi)型的核心區(qū)域多為東部沿海省區(qū)，“低-低”空間集聚類(lèi)型的核心區(qū)域?yàn)槲鞑渴^(qū)；四川則是“高-低”污染類(lèi)型的核心區(qū)域。（2）環(huán)境污染與區(qū)域城市化的關(guān)系呈現(xiàn)倒N型曲線。工業(yè)污染和生活污染都隨著區(qū)域城市化率增加，呈現(xiàn)出倒N型發(fā)展趨勢(shì)，且工業(yè)污染的倒N型曲線更明顯；并且大多數(shù)省份都處在隨著城市化率的增加，污染物增多的階段。這與我國(guó)目前現(xiàn)狀基本相符。（3）各省區(qū)的環(huán)境污染與其他控制變量間存在明顯的相互關(guān)系。工業(yè)污染與工業(yè)化率、消費(fèi)結(jié)構(gòu)、貿(mào)易開(kāi)放度呈現(xiàn)正相關(guān)性，與人口密度、技術(shù)進(jìn)步呈負(fù)相關(guān)；生活污染與消費(fèi)結(jié)構(gòu)、私人汽車(chē)擁有量、餐飲業(yè)企業(yè)營(yíng)業(yè)收入呈現(xiàn)正相關(guān)關(guān)系，與人口密度呈現(xiàn)負(fù)相關(guān)關(guān)系。（4）各省區(qū)的環(huán)境污染在空間維度相互影響。各省區(qū)間的工業(yè)污染物排放和生活污染物排放都具有正向空間相關(guān)性，且工業(yè)污染的空間維度相關(guān)性更強(qiáng)一些。

為了實(shí)現(xiàn)城市化與環(huán)境污染的同步、協(xié)調(diào)、均衡發(fā)展，提出以下建議：

1. 聚焦城鎮(zhèn)人口發(fā)展，合理調(diào)整人口政策。中國(guó)的大部分省區(qū)目前還處于隨著區(qū)域城市化率提高，環(huán)境污染加劇的階段，并且在很長(zhǎng)的一段時(shí)間都將處于這種狀態(tài)。地方政府不能為了盡快實(shí)現(xiàn)環(huán)境污染隨城市人口比重增加而減少的情況，盲目地采取促進(jìn)城市人口增加的政策。人口的增長(zhǎng)應(yīng)該和經(jīng)濟(jì)社會(huì)的發(fā)展相協(xié)調(diào)。過(guò)度的城市人口，會(huì)給環(huán)境造成更大的壓力，乃至最后造成整個(gè)環(huán)境系統(tǒng)承載能力的崩潰。因此，必須聚焦城鎮(zhèn)人口變化，在人口發(fā)展偏離自然增長(zhǎng)時(shí)，要運(yùn)用一定的政策手段給予合理控制，使人口恢復(fù)自然增長(zhǎng)，保證穩(wěn)定地越過(guò)轉(zhuǎn)折點(diǎn)。

2. 發(fā)展“綠色工業(yè)”，優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)。目前，環(huán)境污染隨著工業(yè)產(chǎn)值比率的提升而惡化。這是因?yàn)樵谖覈?guó)目前的工業(yè)結(jié)構(gòu)中，依然是高污染的重工業(yè)占主要地位。為了控制環(huán)境污染，就要逐步淘汰工業(yè)結(jié)構(gòu)中的高污染的產(chǎn)業(yè)，積極引進(jìn)高新技術(shù)，改造傳統(tǒng)產(chǎn)業(yè)，建立起低消耗、低污染、高效益的工業(yè)結(jié)構(gòu)。大力倡導(dǎo)生態(tài)經(jīng)濟(jì)、循環(huán)經(jīng)濟(jì)和低碳經(jīng)濟(jì)。另一方面，要大力發(fā)展具有高產(chǎn)出、高就業(yè)、低消耗和低污染等特點(diǎn)的第三產(chǎn)業(yè)，如金融、商貿(mào)、科教等，使其逐步替代工業(yè)在生產(chǎn)總值中的頭等地位，促進(jìn)整個(gè)經(jīng)濟(jì)向更合理的結(jié)構(gòu)轉(zhuǎn)變。

3. 加強(qiáng)產(chǎn)業(yè)轉(zhuǎn)移質(zhì)量評(píng)估，促進(jìn)貿(mào)易開(kāi)放的技術(shù)效應(yīng)?！拔廴咎焯谩奔僬f(shuō)在本文的實(shí)證中不成立，但這并不意味著我們可以肆無(wú)忌憚地接受?chē)?guó)際貿(mào)易帶來(lái)的產(chǎn)業(yè)轉(zhuǎn)移。在整體趨勢(shì)上，中國(guó)以其發(fā)展中國(guó)家的地位，無(wú)法趨避資源消耗型、污染密集型產(chǎn)業(yè)的轉(zhuǎn)移進(jìn)入，但可以從中選優(yōu)。地方政府應(yīng)盡量發(fā)揮產(chǎn)業(yè)轉(zhuǎn)移中外來(lái)的先進(jìn)環(huán)境管理理念和污染防治技術(shù)；更多地接受高端制造業(yè)、高新技術(shù)產(chǎn)業(yè)、新能源和節(jié)能環(huán)保等產(chǎn)業(yè)的轉(zhuǎn)移，實(shí)現(xiàn)區(qū)域經(jīng)濟(jì)發(fā)展與減少污染排放的雙贏目標(biāo)；在合理擴(kuò)大貿(mào)易開(kāi)放度的同時(shí)，加強(qiáng)對(duì)產(chǎn)業(yè)轉(zhuǎn)移的質(zhì)量評(píng)估，盡量減少以污染換增長(zhǎng)的發(fā)展模式。

4. 關(guān)注地區(qū)間相互影響，促進(jìn)地區(qū)間和諧發(fā)展。我國(guó)地區(qū)間的環(huán)境污染具有較強(qiáng)的正向空間相關(guān)性。也就是說(shuō)，在實(shí)現(xiàn)本地的經(jīng)濟(jì)發(fā)展的同時(shí)有可能犧牲周邊地區(qū)的環(huán)境質(zhì)量。為了全國(guó)整體環(huán)境水平的提高，地區(qū)政府要避免“短視”行為，不能只關(guān)注當(dāng)?shù)氐慕?jīng)濟(jì)發(fā)展，深刻牢記“一榮俱榮、一損俱損”的思想，高度重視相鄰地區(qū)間的相互影響性，做到相鄰省區(qū)間相互關(guān)照、和諧發(fā)展。在制定經(jīng)濟(jì)規(guī)劃、環(huán)境規(guī)劃和社會(huì)規(guī)劃時(shí)，可采取跨地區(qū)合作的方式，減緩地區(qū)間經(jīng)濟(jì)發(fā)展與環(huán)境污染的矛盾，促進(jìn)地區(qū)間的共同和諧發(fā)展。

注釋?zhuān)?/p>

①?gòu)闹袊?guó)目前的發(fā)展來(lái)看，工業(yè)化多集中于城市發(fā)展，農(nóng)村地區(qū)工業(yè)化發(fā)展較低，相對(duì)所占份額較少；同時(shí)由于統(tǒng)計(jì)數(shù)據(jù)中缺少單獨(dú)的城市工業(yè)化指標(biāo)，因此用各省區(qū)全省的工業(yè)化指標(biāo)代替城市地區(qū)的工業(yè)化指標(biāo)。

②西藏部分?jǐn)?shù)據(jù)缺失，海南由于其地理位置與其他地區(qū)不相鄰，無(wú)法考慮空間相關(guān)性，因此將二者剔除。

[1]Howard ． E 著,金經(jīng)元譯．明日的田園生活[M]．北京:北京城市出版社, 1998．

[2]Meadows D H, Meadows D L, Randers J, et al． The Limits to Growth[M]．New York: Universe Books, 1972．

[3]Shafik N, Bandyopadhyay S． Economic growth and environmental quality: time-series and cross-country evidence[R]． World Bank Policy Research Working Paper, 1992, No．904．

[4]Grossman G, Krueger A． Environmental Impacts of A North American Free Trade Agreement[R]． NBER Working Paper, 1991, No．3914．

[5]Vernon Henderson． The Urbanization Process and Economic Growth: the So-what Question[J]． Journal of Economic Growth, 2003(8): 47-71．

[6]姜乃力．城市化對(duì)大氣環(huán)境負(fù)面影響及其對(duì)策[J] 遼寧城鄉(xiāng)環(huán)境科技,1999(02): 63-66．

[7]史愛(ài)玲,閆慶松．城市化對(duì)環(huán)境的影響和對(duì)策[J]．山東環(huán)境,1999(01): 35-36．

[8]盛學(xué)良,董雅文,John Toon．城市化對(duì)生態(tài)環(huán)境的影響和對(duì)策[J]．環(huán)境導(dǎo)報(bào),2001 (06): 8-9．

[9]楊文舉,孫海寧．淺析城市化進(jìn)程中的生態(tài)環(huán)境問(wèn)題[J]．生態(tài)經(jīng)濟(jì),2002 (03): 31-34．

[10]劉耀彬．城市化與生態(tài)環(huán)境耦合機(jī)制及調(diào)控研究[M]．北京:經(jīng)濟(jì)科學(xué)出版社,2007: 134-147．

[11]劉耀彬．城市化與資源環(huán)境相互關(guān)系的理論與實(shí)證研究[M]．北京:中國(guó)財(cái)政經(jīng)濟(jì)出版社,2007: 30-36．

[12]盛光耀．城市化模式與資源環(huán)境的關(guān)系[J] 城市問(wèn)題,2009 (01): 11-17．

[13]蔡孝箴．城市經(jīng)濟(jì)學(xué)[M]．天津:南開(kāi)大學(xué)出版社,1998． 50-53．

[14]Anselin L, Gallo J． L, Jayet H． Spatial Panel Econometrics[A]． Mátyás L．,Sevestre P． The Econometrics of Panel Data[C]． Berlin Heidelberg:Springer, 2008．

[15] Anselin L． Spatial econometrics: methods and models[M]． Dordrecht: Kluwer Academic Publishers,1988．

[16]Elhorst J P． Specification and Estimation of Spatial Panel Data Models[J] International Regional Science Review． 2003,26(3): 244-268．

[17]何江,張馨之．中國(guó)區(qū)域經(jīng)濟(jì)增長(zhǎng)及其收斂性:空間面板數(shù)據(jù)分析[J]．南方經(jīng)濟(jì),2006(5): 44-52．

[18]王火根,沈利生．中國(guó)經(jīng)濟(jì)增長(zhǎng)與能源消費(fèi)空間面板分析[J]．?dāng)?shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究,2007(12): 98-149．

[19]吳擁政,顏日初．中部六省地級(jí)市區(qū)金融發(fā)展與經(jīng)濟(jì)增長(zhǎng)的空間面板分析[J]．統(tǒng)計(jì)與信息論壇, 2009(3): 63-71．

A Spatial Econometric Study of the Relationship between China’s Regional Urbanization and Environmental Pollution

Wang Jiating, Zhao Li, Sun Zhe, Wang Xuan

This paper has analyzed the relationship between regional urbanization and environmental pollution by using spatial econometric method, based on panel data of China’s 29 provinces and autonomous regions from 2000 to 2010. The result shows that: (1) The relationship between regional urbanization and industrial pollution, living pollution is still an inverted N type; (2) Influence degree and direction on environmental pollution varies from other different influencing factors; (3) At present, most of China’s provinces and autonomous regions are in the stage that environmental pollution grows with the increasing urbanization rate.

regional urbanization; environmental pollution; spatial panel data model

TU984.115

10.3969/j.issn.1674-7178.2013.03.001

王家庭，南開(kāi)大學(xué)中國(guó)城市與區(qū)域經(jīng)濟(jì)研究中心副教授，經(jīng)濟(jì)學(xué)博士，碩士生導(dǎo)師，研究方向?yàn)槌鞘信c區(qū)域經(jīng)濟(jì)、土地與房地產(chǎn)經(jīng)濟(jì)；趙麗，南開(kāi)大學(xué)中國(guó)城市與區(qū)域經(jīng)濟(jì)研究中心碩士生，研究方向?yàn)槌鞘信c區(qū)域經(jīng)濟(jì)；孫哲，南開(kāi)大學(xué)中國(guó)城市與區(qū)域經(jīng)濟(jì)研究中心碩士生，研究方向?yàn)槌鞘信c區(qū)域經(jīng)濟(jì)；王璇，南開(kāi)大學(xué)經(jīng)濟(jì)學(xué)碩士，天津國(guó)際郵輪母港有限公司經(jīng)濟(jì)師，研究方向?yàn)槌鞘信c區(qū)域經(jīng)濟(jì)。

（責(zé)任編輯：盧小文）

國(guó)家社會(huì)科學(xué)基金項(xiàng)目《快速城市化進(jìn)程中我國(guó)城市蔓延的成本-收益評(píng)估與治理模式構(gòu)建研究》（12BJY048）成果。