崔紅霞，蒲 龍

（1.浙江財經(jīng)大學東方學院 財稅學院，浙江 海寧 314408；2.浙江財經(jīng)大學 財政稅務學院，杭州 310018）

一、引言

改革開放以來，中國經(jīng)濟在較長的一段時間內(nèi)保持較高的增長速度，被譽為“中國式增長奇跡”。主流觀點認為在以政治上集權(quán)和財政上分權(quán)的中國式分權(quán)和地方官員政治晉升“競標賽”的雙重激勵下，中國地方政府為發(fā)展經(jīng)濟展開了激烈的競爭（Qian 和Weingast，1997；周黎安，2007）。盡管這種自發(fā)的政府間競爭推動了地方經(jīng)濟的快速發(fā)展，但其主要依賴資源能源的高消耗和高投入，形成以要素驅(qū)動和投資驅(qū)動為主要增長引擎的“粗放式發(fā)展模式”。在這一模式下為吸引資本流入和發(fā)展經(jīng)濟，地方政府往往選擇相對寬松的環(huán)境監(jiān)管力度進而對地區(qū)生態(tài)環(huán)境質(zhì)量造成嚴重的影響（Greenstone et al，2021）?？諝馕廴痉矫?，“霧霾”頻頻成為近年來的環(huán)保關鍵詞，2013 年持續(xù)“霧霾”天氣籠罩全國10 余個省份，同年全球空氣質(zhì)量最差的10 個城市中有7 個位于中國。在2014 年兩會期間，李克強總理明確宣布中國政府“向污染宣戰(zhàn)”，此后幾年中國空氣質(zhì)量得到明顯改觀，但截至2019 年，全國337 個地級市PM2.5 濃度均值仍為36 微克/立方米，依然顯著高于世界衛(wèi)生組織（WHO）給出的安全警戒值（10 微克/立方米）①詳見中華人民共和國生態(tài)環(huán)境部網(wǎng)站http：//www.mee.gov.cn/。。水污染現(xiàn)象也較為嚴重，我國一些城市周圍的湖泊大多處于富氧化階段，基本喪失了供水、旅游等功能。2012 年原環(huán)保部公布的重點流域水環(huán)境質(zhì)量顯示全國七大水系中，只有珠江、長江總體水質(zhì)比較良好，松花江為輕度污染，黃河、淮河為中度污染，遼河、海河為重度污染，河流污染相當嚴重。

嚴峻的生態(tài)環(huán)境顯然不利于滿足人民追求美好生活的向往，為解決這一發(fā)展不平衡問題，近年來我國陸續(xù)出臺了一系列政策文件，如《打贏藍天保衛(wèi)戰(zhàn)三年計劃》《水污染防治行動計劃》及《中華人民共和國固體廢物污染環(huán)境防治法》等，根據(jù)不同污染源明確對應的政策措施。但除從法律法規(guī)等政策角度予以矯正外，經(jīng)濟發(fā)展模式的優(yōu)化也可以有效緩解生態(tài)環(huán)境惡化。從現(xiàn)實來看，目前有相當比例的污染物排放集中來源于如石油化工、鋼鐵及印染等傳統(tǒng)制造業(yè)企業(yè)，此類企業(yè)一方面規(guī)模往往較大，且大多為地方支柱企業(yè)，能夠提供穩(wěn)定的就業(yè)和稅收保障，但另一方面，資源能源的高消耗和高投入也易形成污染物的高排放，對當?shù)厣鷳B(tài)環(huán)境帶來不利的影響。隨著我國經(jīng)濟進入以創(chuàng)新驅(qū)動為主要標志的新發(fā)展階段，優(yōu)化資源配置，轉(zhuǎn)變經(jīng)濟發(fā)展模式不僅可以有效提高經(jīng)濟增長質(zhì)量，也有利于鼓勵企業(yè)完成轉(zhuǎn)型升級，實現(xiàn)綠色發(fā)展。目前在諸多促進企業(yè)綠色發(fā)展的實踐中，生態(tài)工業(yè)園區(qū)無疑發(fā)揮著重要的作用。自2000 年前后“生態(tài)工業(yè)”概念引入我國后，原環(huán)保部開始在全國范圍內(nèi)大力推廣生態(tài)工業(yè)園區(qū)建設。生態(tài)工業(yè)園區(qū)兼具經(jīng)濟發(fā)展和生態(tài)環(huán)境保護的雙重目標，更多強調(diào)綠色、低碳、循環(huán)發(fā)展模式，契合新發(fā)展階段下生態(tài)文明建設的內(nèi)涵。但從現(xiàn)有研究來看，已有研究工業(yè)園區(qū)的相關文獻大多仍以國家級經(jīng)濟開發(fā)區(qū)或高新技術開發(fā)區(qū)為例，少數(shù)涉及生態(tài)工業(yè)園區(qū)的研究要么以定性分析為主（田金平等，2016），要么分析生態(tài)工業(yè)園區(qū)對區(qū)域經(jīng)濟增長的作用（蒲龍等，2021；Zeng et al，2021）。然而生態(tài)工業(yè)園區(qū)特色和政策發(fā)力點是以影響企業(yè)行為變化為基礎，但鮮有研究探討生態(tài)工業(yè)園區(qū)的微觀效應。生態(tài)工業(yè)園區(qū)的設立是否能夠顯著降低了企業(yè)污染物的排放水平？其影響機制又是如何？回答這些問題有利于進一步理解生態(tài)工業(yè)園區(qū)對區(qū)域生態(tài)環(huán)境影響的微觀機制，進而也為完善生態(tài)工業(yè)園區(qū)建設及提高經(jīng)濟發(fā)展質(zhì)量提供一定的思路。

鑒于此，本文以生態(tài)工業(yè)園區(qū)為例，實證檢驗其對企業(yè)污染排放的影響，構(gòu)建雙重差分模型實證研究發(fā)現(xiàn)，生態(tài)工業(yè)園區(qū)的設立會顯著降低企業(yè)工業(yè)廢水中化學需氧量的排放水平。然而對于雙重差分法而言，其需要滿足的先決條件是處理組和控制組在事前擁有共同時間變化趨勢，通過事件研究法檢驗這一條件，實證研究發(fā)現(xiàn)在生態(tài)工業(yè)園區(qū)設立之前兩組企業(yè)沒有顯著差異，表明滿足平行趨勢條件。此外，還通過變換指標度量、篩選樣本、添加額外固定效應及考慮環(huán)境規(guī)制力度等方面進行了一系列穩(wěn)健性檢驗，實證結(jié)果均表明基準回歸結(jié)論穩(wěn)健。進一步異質(zhì)性分析發(fā)現(xiàn)生態(tài)工業(yè)園區(qū)對企業(yè)化學需氧量的這一減排效應僅在東部地區(qū)企業(yè)及非國有企業(yè)樣本中顯著，機制分析則表明，生態(tài)工業(yè)園區(qū)的設立通過鼓勵企業(yè)購買更多的廢水治理設施降低了企業(yè)污染物的排放。

二、文獻綜述及研究假設

（一）文獻綜述

與本文研究相關的文獻主要包括兩個方面：工業(yè)園區(qū)的經(jīng)濟社會效應及企業(yè)減排。

首先，已有大量文獻集中討論工業(yè)園區(qū)產(chǎn)生的經(jīng)濟效應。工業(yè)園區(qū)起源于二戰(zhàn)結(jié)束后的部分西方發(fā)達國家，一般是指通過行政手段等方式，規(guī)劃出特定的區(qū)域，聚集生產(chǎn)要素和各類基礎設施，提高工業(yè)化集約程度，突出產(chǎn)業(yè)特色。工業(yè)園區(qū)能夠?qū)Ξ數(shù)亟?jīng)濟產(chǎn)生顯著影響，為進一步發(fā)揮園區(qū)的集聚效應和規(guī)模效應，地方政府也會往往提供給工業(yè)園區(qū)特定的優(yōu)惠政策（place-based policies），如稅收優(yōu)惠、融資等方面。由于能夠更快地獲得資源、基礎設施或人力資本等，工業(yè)園區(qū)也被視為影響企業(yè)選址和投資決策的重要因素之一，大多數(shù)工業(yè)園區(qū)通過特定的優(yōu)惠政策來吸引企業(yè)進駐，進而發(fā)揮規(guī)模效應帶動輻射整個地區(qū)的發(fā)展（Aggarwal，2007；Farole，2013）。雖然工業(yè)園區(qū)通過集聚效應對地區(qū)發(fā)展帶來顯著影響，但也有部分學者認為在工業(yè)園區(qū)影響企業(yè)和人口流動進而帶來額外的資源配置扭曲，此外工業(yè)活動的大量集聚也會帶來一定的擁擠成本（Busso et al，2013）。針對中國工業(yè)園區(qū)的研究，絕大多數(shù)以國家級經(jīng)濟開發(fā)區(qū)或高新技術開發(fā)區(qū)為例，如Wang（2013）以國家級經(jīng)濟開發(fā)區(qū)為例，研究發(fā)現(xiàn)開發(fā)區(qū)的設立可以有效地吸引到FDI，進而可以促進地區(qū)經(jīng)濟增長，劉瑞明和趙仁杰（2015）則以高新技術開發(fā)區(qū)為例，也發(fā)現(xiàn)了類似的結(jié)論。除檢驗對區(qū)域經(jīng)濟增長外，部分文獻還集中討論工業(yè)園區(qū)的微觀效應，即對企業(yè)生產(chǎn)率、土地利用率、TFP 及創(chuàng)新等方面的影響（Lu et al，2019；陳釗和熊瑞祥，2015；張?zhí)烊A和鄧宇銘，2020）。

與大量研究國家經(jīng)濟開發(fā)區(qū)或高新技術開發(fā)區(qū)的經(jīng)濟效應不同，現(xiàn)有研究較少關注生態(tài)工業(yè)園區(qū)對經(jīng)濟社會的影響，其中田金平等（2016）通過定性分析方式對我國生態(tài)工業(yè)園區(qū)的發(fā)展趨勢做了詳細的綜述回顧，蒲龍等（2021）通過實證分析方式，從宏觀角度分析了生態(tài)工業(yè)園區(qū)對區(qū)域經(jīng)濟增長的影響，Zeng et al（2021）則從環(huán)境保護和經(jīng)濟增長雙重目標出發(fā)，實證研究發(fā)現(xiàn)中國生態(tài)工業(yè)園區(qū)的設立可以同時實現(xiàn)上述目標，但上述文獻均未討論生態(tài)工業(yè)園區(qū)對企業(yè)行為的微觀效應。

其次，與本文研究高度相關的另一支文獻集中討論企業(yè)減排效應。環(huán)境規(guī)制是影響企業(yè)排污行為的重要因素之一，從理論上講，嚴格的環(huán)境規(guī)制可以提高企業(yè)生產(chǎn)成本，使得企業(yè)產(chǎn)出降低進而也會減少生產(chǎn)過程中所產(chǎn)生的污染物，此外，嚴格的環(huán)境規(guī)制下，企業(yè)被迫需要購買更多的污染治理設備及提高污染處理的技術水平。因此可以在一定程度上降低其污染物的排放水平（Greenstone 和Hanna，2014；姜英兵和崔廣慧，2019）。然而針對中國的實證研究往往發(fā)現(xiàn)環(huán)境規(guī)制政策對企業(yè)排污并未產(chǎn)生顯著影響，主要原因在于在中國式的分權(quán)激勵下，地方政府執(zhí)行環(huán)境規(guī)制力度較低，為吸引企業(yè)流入，會放松監(jiān)管力度（周飛舟，2006；陳碩，2010）。在低效率的執(zhí)行力度下，污染企業(yè)通過尋租、議價等方式反而會增加污染物的排放（Chichilnisky，1994；袁凱華和李后建，2015），然而即使政府提高環(huán)境規(guī)制力度，同樣會面臨發(fā)展經(jīng)濟和環(huán)境保護的權(quán)衡，因為在嚴格的環(huán)境規(guī)制下，會顯著影響地區(qū)就業(yè)、生產(chǎn)水平等方面（Liu et al，2021）。

除環(huán)境規(guī)制力度外，外資（FDI）也是影響企業(yè)污染排放的重要因素之一，已有研究發(fā)現(xiàn)FDI 對企業(yè)排污行為具有顯著的影響，但同時具有正反兩種影響方向，學術界將其總結(jié)為“污染天堂假說”和“污染光環(huán)假說”。前者主要是指在自由貿(mào)易的情況下，污染較重產(chǎn)業(yè)會從較高環(huán)境標準國家流入相對較低環(huán)境標準的國家，進而引起后者的生態(tài)環(huán)境不斷惡化（Lucas et al，1992），針對中國的部分實證研究也發(fā)現(xiàn)類似的結(jié)論（Ren et al，2014；張宇和蔣殿春，2014）。但與此相反的另一觀點則認為外資通過帶來先進的環(huán)保理念和技術，有助于改善東道國的生態(tài)環(huán)境質(zhì)量，形成“污染光環(huán)假說”（Birdsall 和Wheeler，1993；Eskeland 和Harrison，2003）。針對中國方面的研究，也有實證文獻能夠證實這一觀點（盛斌和呂越，2012；蘇丹妮和盛斌，2021）。

（二）研究假設

工業(yè)園區(qū)可以為企業(yè)集聚提供一定的便利，但由于生產(chǎn)經(jīng)營活動的大量聚集也會對當?shù)厣鷳B(tài)環(huán)境產(chǎn)生一定的影響，甚至有的工業(yè)園區(qū)一度成為污染的代名詞（田金平等，2016）。相比于其他工業(yè)園區(qū)，生態(tài)工業(yè)園區(qū)的主要目標就是要實現(xiàn)綠色、低碳、循環(huán)發(fā)展方式。因此一方面來看，生態(tài)工業(yè)園區(qū)的設立鼓勵企業(yè)利用清潔的能源，實現(xiàn)產(chǎn)業(yè)鏈上下游企業(yè)間的資源能源循環(huán)流動，降低企業(yè)的污染排放水平；另一方面來看，由于地區(qū)環(huán)境質(zhì)量直接與地方官員績效掛鉤，進一步強化了地方官員的政治激勵，這也會在一定程度上加強地方政府環(huán)境監(jiān)管力度。而根據(jù)環(huán)境規(guī)制理論，規(guī)制力度的提高會對企業(yè)產(chǎn)生“創(chuàng)新補償機制”，即進一步鼓勵企業(yè)優(yōu)化資源配置水平及提高創(chuàng)新能力，進而降低污染物排放水平?；诖耍疚奶岢鲅芯考僭O1：

生態(tài)工業(yè)園區(qū)的設立會顯著降低企業(yè)污染物排放（H1）。

但在實際分析生態(tài)工業(yè)園區(qū)對企業(yè)污染物的減排效應時，還需要進一步考慮不同地區(qū)及不同性質(zhì)企業(yè)的異質(zhì)性影響。首先，對于位于不同經(jīng)濟發(fā)展程度的企業(yè)而言，生態(tài)工業(yè)園區(qū)減排效應會存在一定差異。根據(jù)環(huán)境庫茲涅茨曲線理論，即在經(jīng)濟增長初期，環(huán)境污染隨著人均收入的增長而增長，但到一定發(fā)展階段后，會隨著人均收入的增長而降低?？傊h(huán)境污染與經(jīng)濟增長之間會呈現(xiàn)出一定的“倒U 型”關系。因此根據(jù)這一理論，認為在經(jīng)濟發(fā)達地區(qū)，對生態(tài)環(huán)境質(zhì)量的要求更高，而相對欠發(fā)達地區(qū)，為追求經(jīng)濟增長，可以選擇犧牲一定的環(huán)境質(zhì)量。其次，對于不同性質(zhì)企業(yè)而言，生態(tài)工業(yè)園區(qū)的減排效應也具有一定的差異。相對于其他企業(yè)，國有企業(yè)天然具有一定的“政治優(yōu)勢”，在面臨外生沖擊時具有一定的協(xié)調(diào)和議價能力。而對于其他類型企業(yè)而言，則更容易受到政策沖擊的影響?；诖?，本文提出研究假設2：

生態(tài)工業(yè)園區(qū)減排效應可能更多體現(xiàn)在位于東部發(fā)達地區(qū)企業(yè)樣本及非國有企業(yè)樣本中（H2）。

從理論上講，影響企業(yè)污染物排放的因素主要包括三個方面，生產(chǎn)能力、治理設施及治理效率。一般而言，生產(chǎn)能力越低、治理設施越多或治理效率越高會顯著降低企業(yè)污染物排放。但對于生態(tài)工業(yè)園區(qū)而言，因為其自身大多為國家級經(jīng)濟開發(fā)區(qū)或高新技術開發(fā)區(qū)，享有較多的優(yōu)惠政策。因此不會通過降低企業(yè)生產(chǎn)能力的方式實現(xiàn)減排效應。而對治理設施和治理效率而言，認為企業(yè)可能會更加偏重前者。一方面，根據(jù)《國家生態(tài)工業(yè)示范園區(qū)標準》，污染處理設施作為申請生態(tài)工業(yè)園區(qū)的必備條件之一，客觀需要進行大量購置；另一方面，相對于治理效率而言，購買更多治理設施更具直觀性，也更容易達到?；诖?，本文提出研究假設3：

生態(tài)工業(yè)園區(qū)的減排效應主要通過購買更多治理設施來實現(xiàn)（H3）。

三、研究設計及數(shù)據(jù)來源

（一）研究設計

本文實證檢驗生態(tài)工業(yè)園區(qū)設立是否會影響企業(yè)污染物的排放水平，由于生態(tài)工業(yè)園區(qū)的設立兼具時間-地區(qū)層面的差異。因此構(gòu)建如下的雙重差分模型：

其中：被解釋變量Yijt表示在第j市內(nèi)的第i個企業(yè)在t年的污染排放水平，對于污染物的選擇方面，主要以企業(yè)排放的化學需氧量（COD）來衡量，該指標是指在一定條件下，采用強氧化劑處理水樣時所消耗的氧化劑量，一般而言，化學需氧量越大，說明水樣中所含有機物含量越高，水體受有機物的污染越嚴重?；瘜W需氧量的度量單位為毫克/升，為消除異方差和偏態(tài)性，將企業(yè)排放的化學需氧量取對數(shù)處理，并將其作為被解釋變量。

Eparkjt為核心解釋變量，具體而言，如果城市j在第t年擁有國家級生態(tài)工業(yè)園區(qū)，那么第t年及以后Eparkjt賦值為1，否則為0。Xit和Wjt為一系列控制變量，同時包括企業(yè)層面和城市層面。企業(yè)層面控制變量中，控制了企業(yè)規(guī)模、企業(yè)資本勞動比及企業(yè)年齡，其中企業(yè)規(guī)模利用企業(yè)當年產(chǎn)值的對數(shù)值來表示，企業(yè)資本勞動比用企業(yè)總資產(chǎn)與年末就業(yè)人數(shù)之比的對數(shù)值表示，企業(yè)年齡利用企業(yè)設立年份與當年之差取對數(shù)值表示。城市層面控制變量中，控制了城市經(jīng)濟發(fā)展水平并用GDP 對數(shù)值來表示，為了檢驗是否存在“庫茲尼茨曲線”，也控制了GDP 對數(shù)值的平方項，此外，城市的人口密度及產(chǎn)業(yè)結(jié)構(gòu)均可能對企業(yè)排污產(chǎn)生影響。因此也控制了這兩個變量，其中人口密度利用年末總?cè)丝诔鞘休爡^(qū)面積表示，產(chǎn)業(yè)結(jié)構(gòu)利用第二產(chǎn)業(yè)產(chǎn)值占GDP 比重來衡量。μi和μt分別表示企業(yè)固定效應和年份固定效應；εit為隨機擾動項。考慮到核心解釋變量為城市層面變量，并且對于企業(yè)排污而言，不同行業(yè)間差異較大。因此借鑒Cameron et a（l2011）的處理，選擇雙向聚類標準誤，將模型的標準誤同時聚類到城市-行業(yè)層面上，即容許同一城市同一行業(yè)內(nèi)的企業(yè)間存在一定的相關性。

（二）數(shù)據(jù)、變量及描述性統(tǒng)計

作為繼國家級經(jīng)濟開發(fā)區(qū)和國家高新技術產(chǎn)業(yè)開發(fā)區(qū)后的新一代工業(yè)園區(qū)，生態(tài)工業(yè)園區(qū)更多強調(diào)的是綠色、低碳及循環(huán)發(fā)展模式，兼顧經(jīng)濟發(fā)展和生態(tài)建設（蒲龍等，2021）。生態(tài)工業(yè)園區(qū)建設最早可以追溯到2000 年，早期的試點主要集中在制糖、化工及冶煉等專業(yè)園區(qū)，此后原環(huán)保部著力在國家級經(jīng)濟開發(fā)區(qū)及國家高新技術產(chǎn)業(yè)開發(fā)區(qū)開展試點示范工作②《國家生態(tài)工業(yè)示范園區(qū)管理辦法》明確指出“鼓勵國家級經(jīng)濟技術開發(fā)區(qū)、國家高新技術產(chǎn)業(yè)開發(fā)區(qū)、發(fā)展水平較高的省級工業(yè)園區(qū)或其他特色園區(qū)，積極開展生態(tài)園區(qū)創(chuàng)建活動”。。截至2016 年底，由生態(tài)環(huán)境部、商務部和科技部三部委聯(lián)合命名的生態(tài)工業(yè)園區(qū)共計93 家，其中正式獲得命名的生態(tài)園區(qū)有48 家，通過規(guī)劃論證得到批準建設的生態(tài)園區(qū)有45 家，覆蓋了除西藏、青海、貴州、甘肅及海南以外的全部大陸省份。在定義處理組時沒有進一步區(qū)分正式命名生態(tài)園區(qū)和批準建設生態(tài)園區(qū)的差別，只要企業(yè)所在城市只要有一個生態(tài)園區(qū)（不管是批準建設還是正式命名），就把它定義為處理組，其余城市的企業(yè)定義為控制組③采用簡易方式處理的主要原因有兩個，首先，批準建設的生態(tài)工業(yè)園區(qū)與正式獲得命名園區(qū)基本享用一致的政策措施，包括不限于資金、招商引資、對外經(jīng)濟技術合作、服務、稅收優(yōu)惠、人才引進以及土地劃撥。其次，我們研究的企業(yè)樣本跨度為1998—2013 年，而正式獲得命名生態(tài)工業(yè)園區(qū)的最早年份在2008 年，所以如果直接以正式命名園區(qū)來定義處理組的話會造成大量的樣本損失。而最早批準建設生態(tài)園區(qū)的年份為2001 年，因此以包含正式命名和批準建設的生態(tài)工業(yè)園區(qū)來定義處理組。。

企業(yè)層面信息主要來源于兩個微觀數(shù)據(jù)庫，其中關于企業(yè)財務信息數(shù)據(jù)主要來源于《中國工業(yè)企業(yè)數(shù)據(jù)庫》，該數(shù)據(jù)庫樣本區(qū)間為1998—2013 年，涵蓋了絕大多數(shù)規(guī)模以上工業(yè)企業(yè)，也在學術界得到了廣泛應用。而企業(yè)層面污染信息數(shù)據(jù)主要來源于《中國工業(yè)企業(yè)污染排放數(shù)據(jù)庫》，該數(shù)據(jù)庫詳細披露了重點調(diào)查工業(yè)企業(yè)污染排放信息，是目前國內(nèi)所能獲得的最為詳細的企業(yè)層面污染數(shù)據(jù)庫，國內(nèi)應用較少。該數(shù)據(jù)庫是由國家統(tǒng)計局所收集的工業(yè)企業(yè)上報的原始數(shù)據(jù)，也是歷年發(fā)布的《環(huán)境統(tǒng)計年鑒》的具體數(shù)據(jù)來源。調(diào)查方法為對重點調(diào)查工業(yè)企業(yè)逐個發(fā)表填報匯總，其中重點調(diào)查單位是指篩選出的排污量占各地區(qū)排污總量85%以上的工業(yè)企業(yè)單位，該數(shù)據(jù)庫目前可獲得的樣本期間也是1998—2013 年。由于《中國工業(yè)企業(yè)數(shù)據(jù)庫》和《中國工業(yè)企業(yè)污染排放數(shù)據(jù)庫》覆蓋樣本期間相同，此外兩者報告了根據(jù)統(tǒng)一標準設定的企業(yè)身份信息。因此為合并這兩個數(shù)據(jù)庫提供了可能（陳登科，2020）。除企業(yè)層面信息外，城市層面相關數(shù)據(jù)則來自對應年份的《中國區(qū)域經(jīng)濟統(tǒng)計年鑒》和《城市統(tǒng)計年鑒》。

參照陳登科（2020）的匹配思路，在具體匹配時，首先我們按照標準的處理思路（Brandt et al，2009），將《中國工業(yè)企業(yè)數(shù)據(jù)庫》進行清洗，得到企業(yè)層面的面板數(shù)據(jù)，其次按照類似的處理思路，對《中國工業(yè)企業(yè)污染數(shù)據(jù)庫》進行清洗，然后按照企業(yè)代碼將兩個數(shù)據(jù)庫進行匹配，形成了包含污染信息和財務信息的企業(yè)數(shù)據(jù)庫。最后根據(jù)企業(yè)注冊地，將企業(yè)數(shù)據(jù)庫和城市數(shù)據(jù)庫進行合并，最終形成了本文實證研究的數(shù)據(jù)基礎，即時間跨度為1998—2013 年企業(yè)層面的面板數(shù)據(jù)。此外為防止異常值造成的干擾，對連續(xù)變量指標在1%和99%的分位數(shù)水平上進行了去尾處理，表1 匯報了基本描述性統(tǒng)計結(jié)果。

表1 描述性統(tǒng)計

四、實證結(jié)果

（一）基準回歸

根據(jù)模型（1）設定，首先評估生態(tài)工業(yè)園區(qū)對企業(yè)化學需氧量排放的影響，回歸結(jié)果見表2。

表2 基準回歸結(jié)果

表2 匯報了基準回歸結(jié)果，其中（1）列沒有添加任何控制變量，僅僅控制了企業(yè)固定效應和年份固定效應，實證結(jié)果顯示，生態(tài)工業(yè)園區(qū)的設立顯著降低了其所在城市的企業(yè)化學需氧量的排放水平。（2）列和第（3）列分別在（1）列的基礎上逐步添加了企業(yè)層面和城市層面控制變量，從表2 可以看出添加額外控制變量后，生態(tài)工業(yè)園區(qū)對企業(yè)排污的影響系數(shù)變小，但依舊在1%的顯著性水平上顯著為負，以包含全部控制變量的（3）列為例，實證結(jié)果顯示，生態(tài)工業(yè)園區(qū)的設立會顯著降低企業(yè)化學需氧量排放水平13.91%。各控制變量中，企業(yè)規(guī)模和企業(yè)年齡均對企業(yè)排污具有顯著的正向影響，而企業(yè)資本勞動比的影響不顯著。城市層面中，經(jīng)濟發(fā)展水平對污染排放的影響呈現(xiàn)出顯著倒U 型的“庫尼涅茲曲線”，即隨著城市經(jīng)濟發(fā)展會顯著提高當?shù)仄髽I(yè)的污染排放量，但當經(jīng)濟達到一定水平后，對企業(yè)排污會產(chǎn)生邊際遞減的作用。城市人口密度和產(chǎn)業(yè)結(jié)構(gòu)對企業(yè)排污具有顯著的正向影響，一般而言城市人口越密集、第二產(chǎn)業(yè)占比越高就企業(yè)排污水平也會越大，基本和現(xiàn)實吻合。

表2 的后兩列則在基準回歸的基礎上進一步展開穩(wěn)健性討論，其中（4）列在（3）列的基礎上額外控制了企業(yè)期初特征（企業(yè)可觀測到的最早年份）與年份虛擬變量的交互項，其中期初特征主要包括企業(yè)期初的規(guī)模、資本勞動比及企業(yè)年齡。通過這種方式可以剔除不同初始特征企業(yè)在以后年份可能面對的差異性沖擊，進一步降低內(nèi)生性問題帶來的干擾。按照類似的邏輯，（5）列則額外控制了城市期初特征（1998 年）與年份虛擬變量的交互項，其中城市期初特征主要包括期初的GDP 水平及對數(shù)項、人口密度和產(chǎn)業(yè)結(jié)構(gòu)等。從（4）列和（5）列的結(jié)果來看，額外控制這些期初變量與年份虛擬變量的交互項后實證結(jié)果不會發(fā)生明顯改變，即生態(tài)工業(yè)園區(qū)依然會顯著降低企業(yè)化學需氧量的排放水平。

（二）事件研究法

基準回歸結(jié)論表明生態(tài)工業(yè)園區(qū)會顯著降低企業(yè)化學需氧量的排放水平，但就雙重差分法而言，其必須滿足的一個先決條件是處理組和控制組滿足事前平行性趨勢，否則就可能存在同時影響處理組和控制組且隨時間可變的變量，進而引入內(nèi)生性問題，使得估計結(jié)論不可靠。因此，借鑒Jacobson et al，（1993）的處理方法，利用事件研究法來識別處理組和控制組在生態(tài)工業(yè)園區(qū)設立前后的動態(tài)變化，具體而言，設定實證模型如下：

其中：βk代表生態(tài)工業(yè)園區(qū)時間窗口的影響系數(shù)；Dtj0+k代表生態(tài)工業(yè)園區(qū)前后的事件窗口虛擬變量；tj0代表城市j設立生態(tài)工業(yè)園區(qū)的年份；tj0+k代表生態(tài)工業(yè)園區(qū)設立前后的各個年份，其中k取值分別是-6、-5、-4、-3、-2、-1、0、1、2、3、4+④4+表示第4 年及以后。，缺省組是t≤-7；其余控制變量和基準模型類似；標準誤仍然聚類到城市-行業(yè)層面，回歸結(jié)果如圖1 所示。

圖1 生態(tài)工業(yè)園區(qū)設立前后對企業(yè)化學需氧量的動態(tài)影響

圖1 繪制了Dtj0+k的回歸系數(shù)點估計值和置信區(qū)間，從圖1 可以清晰看出，在生態(tài)工業(yè)園區(qū)設立之前年份，企業(yè)化學需氧量排放水平?jīng)]有明顯變化，但在生態(tài)工業(yè)園區(qū)設立后，企業(yè)化學需氧量的排放呈現(xiàn)出顯著的下降趨勢。

（三）穩(wěn)健性檢驗

（1）考慮其他污染指標。在基準回歸中，主要分析企業(yè)化學需氧量排放水平的影響。但在評估水污染時，除化學需氧量外一般常用指標還包括氨氮排放量⑤氨氮排放量主要來源于污水中含氮有機物的分解，焦化、合成氨等工業(yè)廢水，大量氨氮廢水排入水體不僅引起水體富營養(yǎng)化、造成水體黑臭，給水處理的難度和成本加大，甚至對人群及生物產(chǎn)生毒害作用。。因此在表3 的（1）列中，以企業(yè)氨氮排放量對數(shù)值為被解釋變量，實證結(jié)果顯示，生態(tài)工業(yè)園區(qū)對企業(yè)氨氮排放量同樣具有顯著的負面影響。在表3 的（2）列中，直接以企業(yè)生產(chǎn)經(jīng)營過程中所排放的工業(yè)廢水為被解釋變量，實證結(jié)論顯示，生態(tài)工業(yè)園區(qū)對企業(yè)工業(yè)廢水排放也具有顯著的負向影響。因此綜合來看，變換不同被解釋變量后我們基準回歸結(jié)論依然保持穩(wěn)健。

（2）考慮排放強度?；鶞驶貧w中直接以企業(yè)排放化學需氧量取對數(shù)值作為被解釋變量，但從現(xiàn)實來看，不同規(guī)模的企業(yè)存在顯著的排放差異，即使在控制變量中已經(jīng)控制了企業(yè)規(guī)模的影響，但是在穩(wěn)健性檢驗中，進一步考慮到不同企業(yè)排放污染的系統(tǒng)性差異，利用相對值來替代基準回歸中使用的絕對值，具體而言，利用企業(yè)化學需氧量排放量與其工業(yè)總產(chǎn)值比重的對數(shù)值來衡量排污強度。實證結(jié)果見表3 的（4）列，可以生態(tài)工業(yè)園區(qū)設立對企業(yè)排污強度依然具有顯著的負向影響。

（3）刪除2010 年樣本。在基準回歸中，利用的企業(yè)樣本跨度為1998—2013 年，但部分研究也指出2010年《中國工業(yè)企業(yè)數(shù)據(jù)庫》存在指標異常等問題（譚語嫣等，2017）。因此在穩(wěn)健性檢驗部分直接刪除了企業(yè)2010 年樣本，實證結(jié)果見表3 的（3）列，可以看出剔除可能存在的異常樣本后，也不會對基準回歸產(chǎn)生顯著影響。

（4）添加額外固定效應。在基準回歸中，控制了企業(yè)固定效應和年份固定效應。因此可以吸收企業(yè)（包括城市）層面不隨時間變化的不可觀測因素和全國性沖擊影響，但就企業(yè)排污行為而言，不同行業(yè)特征也可能會導致企業(yè)排污出現(xiàn)系統(tǒng)性差異，比如部分從事冶煉、石油化工行業(yè)的企業(yè)與其他類型行業(yè)的企業(yè)在污染物排放中存在顯著差異。因此需要考慮行業(yè)層面帶來的影響。此外，不同省份針對企業(yè)排污也相應出臺了不同的政策制度。因此也需要考慮省份層面帶來的差異。因此在表3 的（4）和（5）列中，分別額外控制了企業(yè)所屬行業(yè)與年份交互項及企業(yè)所處省份與年份交互項⑥此處選擇的是企業(yè)所屬的二位數(shù)行業(yè)代碼，值得注意的是，在本文研究樣本期間內(nèi)，國民經(jīng)濟分類代碼標準發(fā)生過變化，分別是GB/T 4754-1997、GB/T 4754-2002、GB/T 4754-2011。實證處理時，按照不同代碼的具體分類標準，將全部行業(yè)代碼統(tǒng)一轉(zhuǎn)換為GB/T 4754-2002 標準代碼，使得不同年份企業(yè)所屬行業(yè)具有可比性。，進一步消除行業(yè)層面和省級層面帶來的影響，實證結(jié)果顯示，額外控制這些因素后，基準回歸結(jié)果依然保持穩(wěn)健。

（5）考慮環(huán)境規(guī)制力度帶來的影響。理論上講，地方政府的環(huán)境規(guī)制力度會顯著改變企業(yè)排污行為。因此需要剔除環(huán)境規(guī)制力度變化對企業(yè)排污帶來的影響。具體而言，利用企業(yè)所屬城市是否屬于“兩控區(qū)”范圍來考慮環(huán)境規(guī)制力度可能對企業(yè)排污造成的影響。實證結(jié)果見表3 的（6）列，可以看出，相比與基準回歸結(jié)論來看，考慮環(huán)境規(guī)制力度帶來的影響后，生態(tài)工業(yè)園區(qū)對企業(yè)化學需氧量排放的影響系數(shù)略微降低，但仍在1%的顯著性水平上顯著為負。此外是否屬于“兩控區(qū)”范圍對企業(yè)化學需氧量的排放具有顯著的負向影響，且影響系數(shù)較大，這說明環(huán)境規(guī)制力度的確會對企業(yè)排污行為產(chǎn)生顯著影響，但表3 的（6）列回歸結(jié)論同樣說明，即使考慮到環(huán)境規(guī)制力度對企業(yè)排污行為影響后，生態(tài)工業(yè)園區(qū)的設立依然對企業(yè)排污產(chǎn)生顯著影響，基準回歸結(jié)論保持穩(wěn)健。

表3 穩(wěn)健性檢驗

（四）異質(zhì)性分析

首先，我們將企業(yè)分為位于東部地區(qū)及中西部地區(qū)兩類⑦此處沒有進一步細分中部和西部地區(qū)，主要原因在于截至本文的研究期間內(nèi)（1998—2013 年）大多數(shù)生態(tài)工業(yè)園區(qū)均位于東部地區(qū)，此外少部分生態(tài)園區(qū)位于中部，鮮有生態(tài)園區(qū)位于西部地區(qū)，因此，直接僅將樣本按照區(qū)域分為東部地區(qū)和中西部地區(qū)。。實證結(jié)果見表4 前（1）列和（2）列，可以看出，東部地區(qū)生態(tài)園區(qū)的設立會對企業(yè)排污產(chǎn)生顯著的減排效應，但對于中西部地區(qū)而言，生態(tài)工業(yè)園區(qū)的設立對當?shù)仄髽I(yè)化學需氧量排放沒有產(chǎn)生顯著影響，且影響系數(shù)也遠小于中部地區(qū)。之所有產(chǎn)生這一現(xiàn)象的主要原因可能在于東部地區(qū)經(jīng)濟發(fā)展情況較高。因此也會更加注重對生態(tài)環(huán)境保護。因此生態(tài)工業(yè)園區(qū)的設立會直接加大對企業(yè)的減排效應，而對于中西部地區(qū)而言，發(fā)展經(jīng)濟仍是地方政府主要目標之一。因此生態(tài)工業(yè)園區(qū)的設立可能不會對企業(yè)排污產(chǎn)生顯著影響。

其次，按照所有權(quán)性質(zhì)將企業(yè)分為國有企業(yè)、民營企業(yè)和外資企業(yè)⑧對于企業(yè)所有權(quán)性質(zhì)的劃分主要根據(jù)《中國工業(yè)企業(yè)數(shù)據(jù)庫》中企業(yè)注冊信息來進行劃分。對于注冊類型代碼為190 的其他企業(yè)類型，則根據(jù)企業(yè)實收資本中占比最高類型進行劃分。另外外資或港澳臺企業(yè)包括外資企業(yè)和港澳臺企業(yè)。。實證結(jié)果見表4 后三列，可以看出生態(tài)工業(yè)園區(qū)的設立對民營企業(yè)和外資企業(yè)的化學需氧量排放產(chǎn)生顯著的負向影響，對國有企業(yè)的影響雖然為負，但系數(shù)遠低于非國有企業(yè)且不具有統(tǒng)計上的顯著性。出現(xiàn)這一現(xiàn)象的可能原因在于相對于其他類型企業(yè)，國有企業(yè)往往是地區(qū)支柱產(chǎn)業(yè)，能夠為地方財政提供穩(wěn)定的稅源，如煙草、電力等壟斷性行業(yè)，往往與地方政府關系較為密切，具有一定的政策議價能力。而民營企業(yè)和外資企業(yè)對政策變化比較敏感，政策執(zhí)行力度也要明顯高于國有企業(yè)。因此生態(tài)工業(yè)園區(qū)的設立雖然會減少企業(yè)污染排放水平，但其主要是對民營企業(yè)和外資企業(yè)等非國有企業(yè)的影響，對國有企業(yè)的污染排放影響不大。

表4 異質(zhì)性分析

（五）機制檢驗

最后嘗試分析生態(tài)工業(yè)園區(qū)設立如何降低了企業(yè)污染排放水平。從理論上講，影響一個企業(yè)污染物排放的因素主要有三種，企業(yè)生產(chǎn)能力和污染治理設施及污染治理效率。首先，企業(yè)生產(chǎn)能力會直接影響到污染物的排放，如果生態(tài)工業(yè)園區(qū)設立后降低企業(yè)生產(chǎn)水平，那么自然污染物排放量也會變少。其次，污染治理投入也是影響企業(yè)污染物排放水平的重要因素之一，在更加嚴格的環(huán)境規(guī)制力度下，企業(yè)不得不加大對污染治理的投入，但企業(yè)在加大對污染治理的投入時可以面臨兩種選擇，要么企業(yè)進一步加大投入力度，即購買更多的污染減排設備，要么企業(yè)在既有設備的基礎上，進一步改善治理技術水平，提高治理效率。

因此在表5 中，匯報了相應的機制檢驗，在（1）列中，以企業(yè)報告的工業(yè)總產(chǎn)值的對數(shù)值來衡量企業(yè)生產(chǎn)能力，從回歸結(jié)果來分析可知，生態(tài)工業(yè)園區(qū)的設立對企業(yè)工業(yè)總產(chǎn)值沒有發(fā)生顯著的影響，這說明生態(tài)工業(yè)園區(qū)影響企業(yè)排污并不是通過影響企業(yè)生產(chǎn)能力的機制。在表5 的后兩列中，進一步實證檢驗了污染治理投入機制，其中（2）列直接以企業(yè)當年的廢水治理設施數(shù)作為被解釋變量，（3）列以廢水治理設施的處理能力作為被解釋變量，其中處理能力是以企業(yè)每天可處理廢水噸數(shù)來衡量（實證中將該指標取對數(shù)化處理）。從實證回歸結(jié)果來看，生態(tài)工業(yè)園區(qū)的設立會顯著提高企業(yè)廢水治理設施套數(shù)，但對企業(yè)廢水治理能力沒有發(fā)生顯著的影響。綜合機制檢驗的實證結(jié)果來看，生態(tài)工業(yè)園區(qū)主要是通過鼓勵企業(yè)加大對污染治理投入力度，而非影響企業(yè)生產(chǎn)經(jīng)營能力來影響企業(yè)污染物排放水平，具體而言，生態(tài)工業(yè)園區(qū)的設立會鼓勵企業(yè)加大購買廢水治理設施，而非提高治理能力的機制對企業(yè)污染物排放產(chǎn)生減排效應。

表5 機制檢驗

五、結(jié)論及政策啟示

隨著我國經(jīng)濟發(fā)展進入新時代，迫切需要實現(xiàn)綠色經(jīng)濟增長轉(zhuǎn)型，在兼顧經(jīng)濟增長的同時盡可能實現(xiàn)資源節(jié)約和環(huán)境改善。本文以生態(tài)工業(yè)園區(qū)為例，利用企業(yè)層面污染信息，實證檢驗綠色發(fā)展對企業(yè)污染物排放的影響，實證研究發(fā)現(xiàn)：生態(tài)工業(yè)園區(qū)的設立可以有效降低企業(yè)工業(yè)廢水中化學需氧量的排放。事件研究法分析顯示，處理組和控制組企業(yè)在生態(tài)工業(yè)園區(qū)設立前化學需氧量的排放并不存在顯著差異，只有在設立之后才出現(xiàn)差異，證明了基準結(jié)論的穩(wěn)健性。此外，我們進一步考慮了多種穩(wěn)健性檢驗策略，實證結(jié)果均表明基準結(jié)論保持穩(wěn)健。通過異質(zhì)性分析發(fā)現(xiàn)，生態(tài)工業(yè)園的減排效應在東部發(fā)達地區(qū)和非國有企業(yè)樣本中更為顯著。最后機制分析結(jié)果顯示，生態(tài)工業(yè)園區(qū)的設立會通過鼓勵企業(yè)更多購買污染治理設施來實現(xiàn)減排效應。

本文研究結(jié)論衍生出的政策啟示有以下幾點。首先，本文研究結(jié)論表明，生態(tài)工業(yè)園區(qū)可以有效地減少企業(yè)污染物的排放水平，但從現(xiàn)實來看，目前我國生態(tài)工業(yè)園區(qū)數(shù)量仍然較少，且大多位于發(fā)達省份的省會城市或計劃單位市，這表明應該進一步加大對生態(tài)工業(yè)園區(qū)的建設力度，鼓勵普通地級市和中西部城市積極申報生態(tài)工業(yè)園區(qū)。其次，機制分析結(jié)果顯示目前生態(tài)工業(yè)園區(qū)主要通過鼓勵企業(yè)加大對污染治理設施數(shù)量來實現(xiàn)減排效應，而對污染治理效率方面影響并不明顯，這表明在下一步發(fā)展中應重點提升企業(yè)污染治理效率。對于污染治理設施而言，雖然短期內(nèi)加大投資的確有助于企業(yè)實現(xiàn)減排，但當污染設施數(shù)量達到一定量后，可能會產(chǎn)生一定的邊際遞減情況（如企業(yè)生產(chǎn)能力限制，導致部分污染治理設施空置等），進而帶來資源配置效率的低下，此外采購污染治理設施也會提高企業(yè)生產(chǎn)經(jīng)營成本。而污染治理能力的提升，即可以有效降低企業(yè)生產(chǎn)成本，也有助于提高企業(yè)創(chuàng)新能力，有利于企業(yè)長期發(fā)展。