中國工業(yè)水污染排放的空間格局及結(jié)構(gòu)演變研究

石敏俊 鄭丹++雷平 袁靜沛

摘要 我國工業(yè)企業(yè)污染物排放的空間格局正在發(fā)生變化。現(xiàn)有的統(tǒng)計數(shù)據(jù)不能提供區(qū)分區(qū)域和行業(yè)層面的工業(yè)企業(yè)污染物排放數(shù)據(jù)，無法刻畫地市尺度上各個工業(yè)部門污染排放的空間格局及其演變。工業(yè)企業(yè)污染物排放的空間認知是制定環(huán)境污染控制和區(qū)域產(chǎn)業(yè)發(fā)展規(guī)劃的科學基礎(chǔ)。本文聚焦于工業(yè)水污染，構(gòu)建了新的工業(yè)企業(yè)水污染排放的空間估算方法，估算分行業(yè)分區(qū)域的工業(yè)水污染排放量，來分析地市尺度上中國工業(yè)水污染排放的空間格局及其演變。本文構(gòu)建的空間估算方法首先將全國分行業(yè)污染排放量依據(jù)產(chǎn)值分解到省級和地市層面，再依據(jù)企業(yè)規(guī)模校正系數(shù)、生產(chǎn)工藝校正系數(shù)和環(huán)境管制能力校正系數(shù)對未考慮區(qū)域差異的排放強度進行校正，最后采用適時修正法對污染物排放量的估算結(jié)果進行部門和區(qū)域的平衡調(diào)整。以基于空間估算方法估算的工業(yè)COD為例，剔除個別異常值后，2010年省級尺度重點行業(yè)工業(yè)COD排放量估算結(jié)果的平均誤差為12.7%，數(shù)據(jù)質(zhì)量優(yōu)良，顯示本估算方法較強的適用性。根據(jù)本文的結(jié)果，發(fā)現(xiàn)2005年至2010年間我國工業(yè)COD的減排主要由造紙及紙制品業(yè)的技術(shù)效應(yīng)貢獻，其中又以西部地區(qū)造紙及紙制品業(yè)貢獻突出，其次是東部地區(qū)，中部地區(qū)略有貢獻，東北部地區(qū)貢獻為負值。其他對工業(yè)COD減排作出較大貢獻的行業(yè)包括食品加工業(yè)和化學原料及化學品制造業(yè)，但不同產(chǎn)業(yè)的區(qū)域差異顯著?；诳沙掷m(xù)性發(fā)展的視角，技術(shù)效應(yīng)與規(guī)模效應(yīng)同樣顯著的產(chǎn)業(yè)應(yīng)當作為區(qū)域的主導產(chǎn)業(yè)予以重點發(fā)展。

關(guān)鍵詞工業(yè)水污染；空間格局；工業(yè)COD；技術(shù)效應(yīng)；規(guī)模效應(yīng)

中圖分類號X24

文獻標識碼A文章編號1002-2104（2017）05-0001-07DOI：10.12062/cpre.20170309

工業(yè)企業(yè)污染物排放超標是水環(huán)境惡化的主要原因。水環(huán)境的主要污染指標為化學需氧量（COD）、五日生化需氧量和總磷。近年來，工業(yè)企業(yè)污染物排放格局正在發(fā)生變化。一方面，“十一五”期間我國大力推進環(huán)境治理，二氧化硫、化學需氧量等污染物排放總量趨于減少；另一方面，隨著部分制造業(yè)企業(yè)開始從沿海地區(qū)向中西部地區(qū)轉(zhuǎn)移[1-4]，制造業(yè)分布格局正在發(fā)生變化，這也直接影響到工業(yè)企業(yè)污染物排放的空間格局。據(jù)國家環(huán)保部統(tǒng)計，2002年至2010年，西部地區(qū)工業(yè)廢氣排放量占全國的比例從22%上升到28%，工業(yè)固體廢棄物產(chǎn)生量的比例從26%上升到31%，工業(yè)COD排放量的比例從21%上升到27%。工業(yè)企業(yè)污染物排放的空間認知是環(huán)境污染控制和區(qū)域產(chǎn)業(yè)發(fā)展規(guī)劃制定的科學基礎(chǔ)，然而，目前的環(huán)境統(tǒng)計尚不能提供分區(qū)域和分行業(yè)的工業(yè)企業(yè)污染物排放數(shù)據(jù)，無法刻畫地市尺度上各個工業(yè)部門污染物排放的空間格局及其演變。據(jù)此，本文聚焦于工業(yè)水污染，在對分行業(yè)分區(qū)域的工業(yè)水污染物排放量進行空間估算的基礎(chǔ)上，構(gòu)建地市尺度分行業(yè)工業(yè)水污染排放強度數(shù)據(jù)庫，并應(yīng)用于分區(qū)域工業(yè)水污染排放的空間格局及結(jié)構(gòu)演變的分析。本研究拓展的數(shù)據(jù)挖掘技術(shù)，可為微觀尺度的環(huán)境政策制定和區(qū)域環(huán)境規(guī)劃提供科學參考。

1工業(yè)水污染排放空間估算方法

已有的空間估算方法大多是通過挖掘可獲得數(shù)據(jù)變量和估算對象之間的內(nèi)在聯(lián)系，構(gòu)建數(shù)學模型，據(jù)此在有限的數(shù)據(jù)限制下對未知變量的空間分布進行估算[5-8]。由于工業(yè)企業(yè)污染物排放的空間格局既與工業(yè)企業(yè)生產(chǎn)能力的空間分布密切相關(guān)，也受到工業(yè)企業(yè)的工藝設(shè)備和技術(shù)水平、污染物處理能力、區(qū)域環(huán)境管治政策等因素的影響，后者主要是影響工業(yè)企業(yè)的污染物排放強度，因此，基于特定變量的空間分布來估算工業(yè)企業(yè)污染物排放的空間格局存在較大的困難。

因此，本文根據(jù)可獲取的有限信息，構(gòu)建了一種新的空間估算方法，盡可能準確地估算分行業(yè)分區(qū)域的工業(yè)水污染物排放量。以中國環(huán)境統(tǒng)計年鑒提供的全國范圍的分行業(yè)工業(yè)水污染物排放量和省區(qū)或地市尺度的工業(yè)水污染物排放總量為控制點，通過估算工業(yè)企業(yè)污染物排放強度的區(qū)域差異，分解得到分行業(yè)分地市的工業(yè)企業(yè)水污染物排放量。

1.1工業(yè)水污染排放估算步驟

工業(yè)企業(yè)污染物排放強度區(qū)域差異的影響因素主要有企業(yè)規(guī)模差異、生產(chǎn)工藝差異、環(huán)境管治能力差異三個方面[9-12]。因此，本文首先將全國分行業(yè)的工業(yè)污染物排放總量根據(jù)各省或各地市的分行業(yè)產(chǎn)值，分解到各省或各地市的各行業(yè)。其次，從企業(yè)規(guī)模差異、生產(chǎn)工藝差異和環(huán)境管治能力差異三個方面，分別構(gòu)建企業(yè)規(guī)模校正系數(shù)、投產(chǎn)時間校正系數(shù)和環(huán)境管治能力校正系數(shù)，對未考慮排放強度的區(qū)域差異、按照產(chǎn)值分解的污染物排放量初始值進行校正。最后，對污染物排放量的估算結(jié)果進行平衡調(diào)整，以保證估算得到的分行業(yè)污染物排放量加和與分地區(qū)污染物排放量加和相等。平衡調(diào)整采用R.A.S.法（適時修正法）進行。R.A.S.法是通過構(gòu)造行比例和列比例，使得行和與列和等于控制量。

在控制產(chǎn)品種類、原料和工藝等污染排放影響因素的條件下，不同規(guī)模的企業(yè)具有不同的排污系數(shù)。每個地區(qū)不同規(guī)模企業(yè)的占比存在差異，基于不同規(guī)模企業(yè)的排污系數(shù)比值和各類規(guī)模企業(yè)的比例，可構(gòu)建企業(yè)規(guī)模校正系數(shù)。投產(chǎn)時間早晚可以反映生產(chǎn)工藝水平差異，污染物排放國家標準中對不同投產(chǎn)時間的企業(yè)規(guī)定不同排污限值，因此基于不同投產(chǎn)時間的排污限值比值和不同投產(chǎn)時間企業(yè)的比例，可構(gòu)建投產(chǎn)時間校正系數(shù)。環(huán)境管治能力與區(qū)域經(jīng)濟發(fā)展水平、污染物集中處理能力、產(chǎn)業(yè)集聚程度密切相關(guān)[13-14]。因此，環(huán)境管治能力校正系數(shù)由污水處理廠集中處理能力、經(jīng)濟發(fā)展水平和產(chǎn)業(yè)集聚程度三部分組成。污染物集中處理能力采用污水處理廠集中處理所去除的工業(yè)污染量占工業(yè)污染物總?cè)コ康谋嚷剩?jīng)濟發(fā)展水平采用人均GDP，產(chǎn)業(yè)集聚程度采用從業(yè)人員占比指標。每個指標采用與全國（全省）平均水平的比值來體現(xiàn)區(qū)域差異。根據(jù)產(chǎn)業(yè)集聚程度指標比重略輕、污染物集中處理能力指標和經(jīng)濟發(fā)展水平指標比重較高的原則賦予不同的指標權(quán)重，將校正后的估算結(jié)果與測試集反復比較，確定相關(guān)指標權(quán)重和參數(shù)。

在具體估算過程中，先是省區(qū)尺度的空間估算，然后是地市尺度的空間估算。地市尺度的空間估算以省級尺度的分行業(yè)估算結(jié)果作為控制點，這樣可以將地市尺度的估算誤差控制在本省范圍之內(nèi)。為保證估算結(jié)果與實測數(shù)據(jù)一致，在進行估算有效性檢驗后，對于存在實測數(shù)據(jù)的地市以實測數(shù)據(jù)替換估算數(shù)據(jù)，并再次進行平衡調(diào)整。

我國工業(yè)污染物排放主要集中在少數(shù)污染密集型行業(yè)。2010年，農(nóng)副食品加工業(yè)、食品制造業(yè)、飲料制造業(yè)、紡織業(yè)、造紙及紙制品業(yè)和化學原料及化學制品制造業(yè)6個行業(yè)合計排放了69.45%的工業(yè)COD。為簡化估算工作量，工業(yè)COD的空間估算以農(nóng)副食品加工業(yè)、食品制造業(yè)、飲料制造業(yè)、紡織業(yè)、造紙及紙制品業(yè)和化學原料及化學制品制造業(yè)6個行業(yè)作為重點行業(yè)，其余行業(yè)歸并為其他行業(yè)。因此，省級尺度上估算了大陸地區(qū)31個省級行政單元6個重點行業(yè)和其他行業(yè)的污染排放，地市尺度上估算了332個地級行政單元6個重點行業(yè)和其他行業(yè)的污染排放，4個直轄市、海南省和西藏自治區(qū)作為一個空間單元進行估算。

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

工業(yè)污染物排放量和污水處理的相關(guān)數(shù)據(jù)來自《中國環(huán)境年鑒》、《中國城市統(tǒng)計年鑒》、《中國城市發(fā)展報告》、各省區(qū)統(tǒng)計年鑒和環(huán)保統(tǒng)計資料；區(qū)域經(jīng)濟的相關(guān)數(shù)據(jù)來自《中國區(qū)域經(jīng)濟統(tǒng)計年鑒》；工業(yè)企業(yè)的相關(guān)數(shù)據(jù)來自中國工業(yè)企業(yè)數(shù)據(jù)庫，其中2010年工業(yè)企業(yè)產(chǎn)值根據(jù)2009年工業(yè)企業(yè)產(chǎn)值結(jié)構(gòu)和價格指數(shù)推算而來；工業(yè)企業(yè)分行業(yè)污染物排放國家標準來自環(huán)保部[15]；工業(yè)企業(yè)分行業(yè)排污系數(shù)來自《第一次全國污染源普查工業(yè)污染源產(chǎn)排污系數(shù)手冊》[16]?；谝延械臄?shù)據(jù)，可以得到分產(chǎn)業(yè)的COD排放量（見表1）。

1.3估算結(jié)果的有效性檢驗

借鑒數(shù)據(jù)挖掘方法的思路，以公布的部分省份分行業(yè)工業(yè)水污染物排放量和部分地市工業(yè)水污染物排放量的實測數(shù)據(jù)作為測試集，對估算結(jié)果進行了有效性檢驗。

在2010年工業(yè)COD排放量估算結(jié)果有效性檢驗中，省級單元的測試集樣本量為30個，地市單元的測試集樣本量為144個。剔除浙江的紡織業(yè)、新疆的紡織業(yè)和山東的化學原料及化學制品制造業(yè)3個異常值，2010年省級尺度重點行業(yè)工業(yè)COD排放量估算結(jié)果的平均誤差為12.70%。2010年地市尺度重點行業(yè)工業(yè)COD排放量估算結(jié)果，剔除異常值后，排放強度低、經(jīng)濟發(fā)達區(qū)域的加權(quán)平均誤差為26.7%，中等排放強度區(qū)域的加權(quán)平均誤差為22.1%，高排放強度區(qū)域的加權(quán)平均誤差為23.5%。

在2005年工業(yè)COD排放量估算結(jié)果有效性檢驗中，省級單元的測試集樣本量為18個，地市單元的測試集樣本量為134個。剔除山東的農(nóng)副食品制造業(yè)和新疆的飲料制造業(yè)2個異常值，2005年省級尺度重點行業(yè)工業(yè)

COD排放量估算結(jié)果的平均誤差為7.85%。2005年地市尺度重點行業(yè)工業(yè)COD排放量估算結(jié)果，剔除異常值后，排放強度低、經(jīng)濟發(fā)達區(qū)域的加權(quán)平均誤差為28.32%，中等排放強度區(qū)域的加權(quán)平均誤差為27.29%，高排放強度區(qū)域的加權(quán)平均誤差為19.76%。

2工業(yè)水污染排放的空間格局及結(jié)構(gòu)演變：以工業(yè)COD為例

2.1基礎(chǔ)研究模型

區(qū)域環(huán)境質(zhì)量的影響因素主要有規(guī)模效應(yīng)、技術(shù)效應(yīng)和結(jié)構(gòu)效應(yīng)三個方面[13]。其中，規(guī)模效應(yīng)是指在其他條件不變的情況下，經(jīng)濟規(guī)模擴大或縮小將伴隨著污染物排放的增加或減少。技術(shù)效應(yīng)是指經(jīng)濟發(fā)展和研發(fā)支出增加推動技術(shù)進步，從而產(chǎn)生兩方面的影響：一是在其他條件不變時，技術(shù)進步提高生產(chǎn)率，改善資源利用效率，降低單位產(chǎn)出的要素投入；二是清潔技術(shù)不斷開發(fā)并取代舊技術(shù)，可以有效地循環(huán)利用資源，降低單位產(chǎn)出的污染排放。規(guī)模效應(yīng)導致污染負荷加劇，技術(shù)效應(yīng)則會抵消污染負荷。

基于污染排放恒等式，可構(gòu)建方程如方程1。

其中COD為工業(yè)COD排放總量，下標指不同區(qū)域或產(chǎn)業(yè)的COD排放量。i、j分別代表區(qū)域和產(chǎn)業(yè)；GDP為區(qū)域地區(qū)生產(chǎn)總值，SGDP為單位地區(qū)生產(chǎn)總值排放強度；GP為工業(yè)總產(chǎn)值，SGP為單位工業(yè)總產(chǎn)值排放強度。

根據(jù)現(xiàn)有的統(tǒng)計數(shù)據(jù)發(fā)布，現(xiàn)在僅有全國層面分行業(yè)的COD排放量，即僅有方程1a有數(shù)據(jù)基礎(chǔ)。更加深入的分區(qū)域經(jīng)濟增長環(huán)境效益研究1c和分區(qū)域經(jīng)濟結(jié)構(gòu)環(huán)境效益研究只有在本文數(shù)據(jù)的基礎(chǔ)上才有可能。由于沒有公布可比較的分區(qū)域COD排放數(shù)據(jù)，甚至1b的區(qū)域COD排放比較都較困難。

全微分方程1d，可得COD排放的時間微分方程2，并進而推導方程3，即COD排放演變貢獻模型。

方程3中d代表差分，反映了不同時期變量的變化量。右側(cè)第一式反映的區(qū)域i產(chǎn)業(yè)j的COD排放規(guī)模效應(yīng)，第二式反映了區(qū)域i產(chǎn)業(yè)j的COD排放技術(shù)效應(yīng)。在第二部分本文推導了332個城市的7個行業(yè)（6個細分行業(yè)和“其他”），因此我們可以在地級城市層面研究332×7=2 324組4 648個規(guī)模效應(yīng)和技術(shù)效應(yīng)對總COD貢獻指標，但這顯然太繁瑣且偏離了本文的主題?；谝陨?，下文首先介紹基于本文數(shù)據(jù)開展的COD排放空間格局與產(chǎn)業(yè)結(jié)構(gòu)演變，其次討論影響COD排放空間格局與產(chǎn)業(yè)結(jié)構(gòu)演變的規(guī)模效應(yīng)和技術(shù)效應(yīng)。

2.2工業(yè)水污染的空間格局與產(chǎn)業(yè)結(jié)構(gòu)演變

根據(jù)公布的全國行業(yè)數(shù)據(jù)，2005年至2010年，我國COD排放量從556萬t下降到435萬t，降幅21.8%。在減排總量中，本文研究的6個行業(yè)減排貢獻率超過85%（見表1），證明了本文選擇的重點行業(yè)的合理性。值得說明的是，模型3計算出來的減排貢獻率指標符號與COD排放量變動方向一致，即本文提到的減排貢獻概念是一個正向概念，某領(lǐng)域COD減排量在減排總量中比重越大，貢獻率越高。

表1顯示，6個產(chǎn)業(yè)中，造紙及紙制品業(yè)貢獻了最多的COD減排量，貢獻率達到60.70%，其次是食品加工業(yè)（16.56%）以及化學原料及化學品制造業(yè)（10.58%）。表1的數(shù)據(jù)還顯示飲料制造業(yè)與紡織業(yè)的COD排放與大勢相反，2005年至2010年間減排貢獻率分別為-3.72%和-1.07%。分產(chǎn)業(yè)COD減排總量研究離開行業(yè)規(guī)模的變化是沒有意義的，因為無法判斷究竟是排放規(guī)模還是排放強度的變化造成了這一差異。當引入行業(yè)規(guī)模后，我們發(fā)現(xiàn)，從2005年至2010年，所有研究行業(yè)的規(guī)模都在擴張，減排的規(guī)模效應(yīng)貢獻全部為負值，真正推動COD減排的是排放強度所反映的技術(shù)效應(yīng)，尤其是造紙及紙制品業(yè)，在規(guī)模擴張超過1倍的背景下，其技術(shù)效應(yīng)造成的COD減排貢獻率達到191.99%，不僅抵消了規(guī)模效應(yīng)造成的增排，還卓有余裕。2005年至2010年間，食品加工業(yè)和化學原料與化學制品制造業(yè)的技術(shù)效應(yīng)也貢獻了較高的COD減排，貢獻率分別達到73.13%和57.04%。表1的數(shù)據(jù)顯示，飲料制造業(yè)和紡織業(yè)的COD排放變動趨勢差異不僅是行業(yè)規(guī)模擴張帶來的規(guī)模效應(yīng)的原因，也是研究年份期間排放強度下降緩慢、技術(shù)效應(yīng)貢獻有限的原因。類似的行業(yè)還有食品制造業(yè)，這些行業(yè)的共同特征為低技術(shù)效應(yīng)和低規(guī)模效應(yīng)共存，二者的內(nèi)在聯(lián)系，值得進一步的探究。

雖然基于全國行業(yè)數(shù)據(jù)能得出一些有意思的結(jié)論，但這些結(jié)論的實踐指導意義仍然值得商榷，原因不僅在中國巨大的空間差異，還在中國產(chǎn)業(yè)空間布局處于不斷重構(gòu)之中。尤其是近十年來，中國制造業(yè)在沿海和內(nèi)陸間的空間轉(zhuǎn)移是否與環(huán)境規(guī)制有關(guān)，產(chǎn)業(yè)空間重構(gòu)是否同步帶來了污染物排放的空間重構(gòu)？什么樣的污染企業(yè)從沿海遷往了中西部？這些企業(yè)是采用了更加消耗資源環(huán)境的技術(shù)？還是借助搬遷機會，在一個更高的平臺生產(chǎn)？這些問題才是更加值得關(guān)心的問題。應(yīng)用本文拓展的分行業(yè)分區(qū)域數(shù)據(jù)，我們可以開展更深入的研究?！糐P〗

表2描述了COD排放的空間結(jié)構(gòu)演變，表3分行業(yè)描述了COD減排的空間差別化貢獻。表2的數(shù)據(jù)證實中國2005年至2010年的COD減排中，東部地區(qū)發(fā)揮了巨大的作用，東部地區(qū)的企業(yè)實現(xiàn)了超過40%的減排。從產(chǎn)業(yè)占全國的比重看，東部地區(qū)的減排效率似乎還有待于進一步的提高。表2的數(shù)據(jù)也緩解了社會各界對西部經(jīng)濟發(fā)展中的環(huán)境污染擔心。至少從COD排放視角，2005年至2010年全國的COD減排西部的貢獻率達到31.93%，相比西部地區(qū)的工業(yè)總產(chǎn)值比重，其環(huán)境效率表現(xiàn)在四大板塊中堪稱首位。中部和東北地區(qū)的COD空間排放結(jié)構(gòu)基本與產(chǎn)業(yè)結(jié)構(gòu)變動一致。

表3提供了進一步的事實。如前文所述，COD減排最大的動力來自造紙及紙制品業(yè)。分區(qū)域、分產(chǎn)業(yè)的數(shù)據(jù)顯示，造紙和紙制品在COD減排中的良好表現(xiàn)主要是西部地區(qū)和東部地區(qū)，在東北，造紙及紙制品業(yè)的COD減排貢獻甚至是負值。在減排貢獻排名第二的食品加工業(yè)，西部地區(qū)也表現(xiàn)優(yōu)良。東部地區(qū)在化學原料及化學制品制造業(yè)COD減排表現(xiàn)優(yōu)良，而同期西部地區(qū)減排貢獻率為負值。此外，西部地區(qū)飲料制造業(yè)COD排放貢獻率為-4.93%，顯著影響了本階段COD減排的成績。

總體看，基于數(shù)據(jù)挖掘技術(shù)，本部分重點討論了COD污染的空間格局與產(chǎn)業(yè)結(jié)構(gòu)演變，發(fā)現(xiàn)西部地區(qū)表現(xiàn)優(yōu)良，造紙及紙制品業(yè)表現(xiàn)優(yōu)良，但也存在區(qū)域和產(chǎn)業(yè)例外。造成這些現(xiàn)象的原因究竟是區(qū)域?qū)用娴漠a(chǎn)業(yè)技術(shù)效應(yīng)還是規(guī)模效應(yīng)，將在下文繼續(xù)討論。

2.3工業(yè)水污染空間格局與產(chǎn)業(yè)結(jié)構(gòu)演變的規(guī)模與技術(shù)效應(yīng)

從可持續(xù)性發(fā)展的視角看，一個地區(qū)的主導產(chǎn)業(yè)應(yīng)當至少具備兩種特征，首先該產(chǎn)業(yè)應(yīng)當具有良好的資源環(huán)境效率，其次還應(yīng)當有良好的成長性。2005年至2010年，在本文所研究的6個行業(yè)中，東部地區(qū)符合以上兩個特征的是造紙及紙制品業(yè)以及化學原料及化學制品制造業(yè)（見表4）。雖然東部地區(qū)造紙及紙制品業(yè)期間規(guī)模急劇擴張，規(guī)模效應(yīng)的COD排放貢獻率為-40.07%，但排放強度的迅速下降使得其技術(shù)效應(yīng)的COD排放貢獻率高達59.67%，從而該行業(yè)的東部COD排放不增反減，表現(xiàn)突出。經(jīng)濟與環(huán)境效益表現(xiàn)更加突出的是西部地區(qū)的造紙及紙制品業(yè)和食品加工業(yè)。與這兩個產(chǎn)業(yè)相比，飲料制造業(yè)與化學原料及化學品制造業(yè)雖然也得到高速增長，但從技術(shù)效應(yīng)的同比變化相比，至少從COD排放角度，飲料制造業(yè)和化

注：2005年山東、新疆的各行業(yè)COD排放強度以及黑龍江的農(nóng)副食品制造業(yè)、造紙及紙制品業(yè)和化學原料及化學制品制造業(yè)COD排放強度為實測數(shù)據(jù)，其余為估算結(jié)果；2010年黑龍江、浙江、山東、廣西和新疆的各行業(yè)工業(yè)COD排放強度為實測數(shù)據(jù)，其余為估算結(jié)果。“-”表示該省份該行業(yè)沒有規(guī)模以上企業(yè)。

學原料及化學品制造業(yè)不是西部地區(qū)的環(huán)境優(yōu)勢產(chǎn)業(yè)。與東部和西部相比，中部地區(qū)和東北地區(qū)相關(guān)產(chǎn)業(yè)都無突出表現(xiàn)。

整體看，不同產(chǎn)業(yè)在不同區(qū)域的技術(shù)效應(yīng)和規(guī)模效應(yīng)差異顯著，尤其是技術(shù)效應(yīng)。同一產(chǎn)業(yè)在不同區(qū)域顯著差異的資源環(huán)境利用效率利用本研究的數(shù)據(jù)拓展技術(shù)得到充分展示?；谶@一技術(shù)，在表4中，我們可以比較直觀地判斷各地的環(huán)境優(yōu)勢產(chǎn)業(yè)，從而得出直觀的政策建議。這一技術(shù)的價值不僅體現(xiàn)在實踐中，更可以為理論界所用，為探討諸如不同行業(yè)的資源環(huán)境利用效率空間差異等熱點問題提供數(shù)據(jù)基礎(chǔ)。

表3的數(shù)據(jù)來自地市級城市的匯總數(shù)據(jù)，因此，相關(guān)的研究還可以深入到省際層面和地市級層面。表5提供了分省分行業(yè)排放強度。根據(jù)表5，結(jié)合相關(guān)行業(yè)發(fā)展情況，我們可以得出比表2—表4更加豐富、微觀的結(jié)論。

3結(jié)論與啟示

由于分行業(yè)分區(qū)域的環(huán)境污染物排放未見統(tǒng)計數(shù)據(jù)發(fā)布，因此我國分行業(yè)污染物排放的空間格局及其演變一直是研究的難點。本文基于工業(yè)企業(yè)水污染物排放量的空間估算方法，在地市級層面估算了分行業(yè)分區(qū)域的COD排放數(shù)據(jù)。這一數(shù)據(jù)挖掘技術(shù)及其形成的數(shù)據(jù)可以為相關(guān)實證分析和政策研究提供新的數(shù)據(jù)基礎(chǔ)，并推動環(huán)境經(jīng)濟研究的深化。

基于本文估算方法估算得到的數(shù)據(jù)，分析了分區(qū)域工業(yè)水污染排放的空間格局及結(jié)構(gòu)演變，估算了2005年和2010年分區(qū)域分行業(yè)的工業(yè)企業(yè)COD排放的空間格局及其結(jié)構(gòu)演變，重點研究了我國工業(yè)COD排放的6個主要行業(yè)，得到如下結(jié)論。

第一，2005—2010年間我國工業(yè)COD減排主要是依賴西部和東部地區(qū)造紙及紙制品業(yè)的技術(shù)效應(yīng)，其中西部地區(qū)尤為顯著。同期為我國工業(yè)COD減排做出較大貢獻的還有東部地區(qū)的化學原料及化學品制造業(yè)、食品加工業(yè)；西部地區(qū)的食品加工業(yè)；中部地區(qū)的造紙及紙制品制造業(yè)。從可持續(xù)發(fā)展的視角，這些產(chǎn)業(yè)部門無論從增長速度還是減排技術(shù)方面，都可以確立為所在區(qū)域的主導產(chǎn)業(yè)。同期東北沒有技術(shù)效應(yīng)和規(guī)模效應(yīng)都表現(xiàn)突出的產(chǎn)業(yè)。

第二，從整體上看，同一產(chǎn)業(yè)在不同區(qū)域的資源環(huán)境效率差異顯著，其原因值得深入研究。

第三，本文基于數(shù)據(jù)挖掘技術(shù)的污染排放空間估算方法可以在更加微觀的空間尺度和更廣泛的行業(yè)與污染物種類層面得到進一步運用。（編輯：李琪）

參考文獻（References）

[1]石敏俊，楊晶，龍文，等.中國制造業(yè)分布的地理變遷與驅(qū)動因素[J].地理研究，2013，32（9）：1708-1720. [SHI Minjun， YANG Jing， LONG Wen， et al. Changes in geographical distribution of Chinese manufacturing sectors and its driving forces [J]. Geographical research， 2013， 32（9）： 1708-1720.]

[2]HE C，WANG J. Regional and sectoral differences in the spatial restructuring of Chinese manufacturing industries during the postWTO period[J]. Geo journal， 2012， 77（3）：361-381.

[3]吳愛芝，孫鐵山，李國平.中國紡織服裝產(chǎn)業(yè)的空間集聚與區(qū)域轉(zhuǎn)移[J].地理學報，2013，68（6）：775-790. [WU Aizhi， SUN Tieshan， LI Guoping. Spatial agglomeration and regional shift of textile and garment industry in China [J]. Acta geographica sinica， 2013， 68（6）： 775-790.]

[4]毛琦梁，董鎖成，王菲，等.中國省區(qū)間制造業(yè)空間格局演變[J].地理學報，2013，68（4）：435-448. [MAO Qiliang， DONG Suocheng， WANG Fei， et al. Evolving spatial distribution of manufacturing industries in China [J]. Acta geographica sinica， 2013， 68（4）： 435-448.]

[5]WANG Jinfeng，REIS B Y，HU Maogui. Area disease estimation based on sentinel hospital records[J].Plos one，2011，6（8）：1-8.

[6]LI Lianfa，WANG Jinfeng，WU Jun. A spatial model to predict the incidence of neural tube defects[J].BMC public health，2012（12）： 951-961.

[7]WANG Jinfeng，HU Maogui，XU Chengdong，et al. Estimation of citywide air pollution in Beijing[J].Plos one， 2013， 8（1）：1-6.

[8]戚偉，李穎，劉盛和，等. 城市晝夜人口空間分布的估算及其特征——以北京市海淀區(qū)為例[J]. 地理學報， 2013， 68（10）：1344-1356. [QI Wei， LI Ying， LIU Shenghe， et al. Estimation of urban population at daytime and nighttime and analyses of their spatial pattern： a case study of Haidian District， Beijing [J]. Acta geographica sinica， 2013， 68（10）： 1344-1356.]

[9]趙海霞，曲福田，郭忠興.環(huán)境污染影響因素的經(jīng)濟計量分析：以江蘇省為例[J].環(huán)境保護，2006（4）：43-49. [ZHAO Haixia， QU Futian， GUO Zhongxing. Environmental pollution and its affecting factors： a theoretical modle and its application to Jiangsu， China [J]. Environmental protection， 2006（4）： 43-49.]

[10]MATTHEW A C，ROBERT J R. ELLIOTT K S. Industrial characteristics， environmental regulations and air pollution： an analysis of the UK manufacturing sector[J].Journal of environmental economics and management，2005， 50（1）：121-143.

[11]JIANG L，LIN C，LIN P. The determinants of pollution levels： firmlevel evidence from Chinese manufacturing [J]. Journal of comparative economics，2014， 42（1）：118-142.

[12]傅京燕.產(chǎn)業(yè)特征、環(huán)境規(guī)制與大氣污染排放的實證研究——以廣東省制造業(yè)為例[J].中國人口·資源與環(huán)境，2009，19（2）：73-77. [FU Jingyan. An empirical analysis on industrial characteristics， environmental regulation and air pollution： Guangdongs manufacturing industry [J]. China population，resources and environment， 2009，19（2）： 73-77.]

[13]GROSSMAN G M，KRUEGER A B. Environmental impacts of a North American free trade agreement [J]. Social science electronic publishing， 1991，8（2）：223-250.

[14]COPELAND B A，TAYLOR M S. NorthSouth trade and the environment [J].Quarterly journal of economics. 1994，109（3）：755-787.

[15]環(huán)境保護部.水污染物排放標準[S].[2013-10-10].http：//kjs.mep.gov.cn/hjbhbz/bzwb/ shjbh /swrwpfbz/.[Ministry of Environmental Protection. Standards for the discharge of water pollutants[S]. [2013-10-10].http：//kjs.mep.gov.cn/hjbhbz/bzwb/ shjbh /swrwpfbz/.]

[16]國務(wù)院第一次全國污染源普查領(lǐng)導小組辦公室.第一次全國污染源普查工業(yè)污染源產(chǎn)排污系數(shù)手冊[S].[2013-10-15].http：//cpsc.mep.gov.cn/.[Office of first national survey of pollution source. Handbook for pollutant discharge coefficients of industrial pollutant sources[S]. [2013-10-15].http：//cpsc.mep.gov.cn/.]

作者簡介：石敏俊，教授，博導，主要研究方向為區(qū)域經(jīng)濟、綠色發(fā)展。Email： mjshi@ruc.edu.cn。

基金項目：國家自然科學基金項目“資源環(huán)境約束下中國經(jīng)濟中長期發(fā)展前景與綠色發(fā)展模式”（批準號：71173212），“環(huán)境規(guī)制政策對產(chǎn)業(yè)空間分布的影響研究：基于環(huán)境便利性的視角”（批準號：71503241），“產(chǎn)業(yè)空間布局的資源環(huán)境效應(yīng)”（批準號：41471112）；國家社會科學基金重大項目“中國經(jīng)濟綠色發(fā)展的理論內(nèi)涵、實現(xiàn)路徑與政策創(chuàng)新”（批準號：15ZDC006）。