●張銀銀 劉云飛

一、引言

2014 年浙江首次提及特色小鎮(zhèn)，之后各地特色小鎮(zhèn)建設(shè)如火如荼，2016 年和2017 年住房城鄉(xiāng)建設(shè)部分別公布了兩批中國特色小鎮(zhèn)建設(shè)名單，共計403 個，省級特色小鎮(zhèn)創(chuàng)建達到2000多個。特色小鎮(zhèn)建設(shè)是推進新型城鎮(zhèn)化的重要組成部分，黨的十八大和十八屆三中全會中提出，“把生態(tài)文明理念和原則全面融入城鎮(zhèn)化全過程，”國家“十三五”規(guī)劃和《國務(wù)院關(guān)于深入推進新型城鎮(zhèn)化建設(shè)的若干意見》明確支持加快發(fā)展中小城市和特色鎮(zhèn)的政策導(dǎo)向。2016 年國家發(fā)改委發(fā)布《關(guān)于加快美麗特色?。ǔ牵╂?zhèn)建設(shè)的指導(dǎo)意見》，要求城鎮(zhèn)建設(shè)優(yōu)化生態(tài)環(huán)境，提升城鎮(zhèn)發(fā)展品質(zhì)，加快建設(shè)美麗特色新型小（城）鎮(zhèn)。2020 年9 月，國家發(fā)改委再次發(fā)文《關(guān)于促進特色小鎮(zhèn)規(guī)范健康發(fā)展的意見》，要求“準(zhǔn)確把握特色小鎮(zhèn)發(fā)展內(nèi)涵，發(fā)揮要素成本低、生態(tài)環(huán)境好、體制機制活等優(yōu)勢，打造經(jīng)濟高質(zhì)量發(fā)展的新平臺、新型城鎮(zhèn)化建設(shè)的新空間、城鄉(xiāng)融合發(fā)展的新支點、傳統(tǒng)文化傳承保護的新載體?！薄巴七M特色小鎮(zhèn)多元功能聚合，打造宜業(yè)宜居宜游的新型空間?！碧厣℃?zhèn)建設(shè)對生態(tài)環(huán)境建設(shè)的影響到底怎樣？學(xué)者們主要從特色小鎮(zhèn)/城鎮(zhèn)化/城市化與生態(tài)環(huán)境協(xié)調(diào)或耦合的關(guān)系研究（方創(chuàng)琳①，劉海猛②；2006，2019）、協(xié)調(diào)關(guān)系的識別（梁振民，2019）③、耦合模式及判別（劉耀彬，2005）④、生態(tài)環(huán)境承載力（唐志強，2019）⑤等方面進行理論研究。特色小鎮(zhèn)的研究主要聚焦于概念（王小章，2016）⑥、特征（衛(wèi)龍寶等，2016）⑦、空間組織形式（盛世豪，張偉明，2016）⑧指標(biāo)衡量體系（吳一洲等，2016）、建設(shè)經(jīng)驗（郁建興等，2016）⑨、定位與功能聚合（王博雅等，2020；易開剛等，2019）⑩?、治理（姚尚建，2017）、產(chǎn)業(yè)鏈空間分置（莊晉財?shù)龋?018）?、空間分布特征（時浩楠等，2019）?、特色指標(biāo)體系（王長松等，2019）等?，現(xiàn)有的文獻對兩者關(guān)系實證研究還鮮有人研究。本文將基于288 個城市中國特色小鎮(zhèn)建設(shè)為研究對象，運用PSM-DID 方法，2009 至2018 年10 年的數(shù)據(jù)進行探索與分析，具有重要的研究意義。

二、特色小鎮(zhèn)建設(shè)對生態(tài)環(huán)境建設(shè)的影響分析

特色小鎮(zhèn)和生態(tài)環(huán)境建設(shè)是相互耦合的開放復(fù)雜系統(tǒng)，蘊含著對立統(tǒng)一的規(guī)律。兩者之間的關(guān)系本質(zhì)上是城鎮(zhèn)化推進或經(jīng)濟發(fā)展與生態(tài)環(huán)境之間的關(guān)系，針對該問題生態(tài)學(xué)、環(huán)境科學(xué)、經(jīng)濟學(xué)、管理學(xué)等各學(xué)科領(lǐng)域都有深入研究。西方學(xué)者主要是圍繞人類活動與自然環(huán)境之間的關(guān)系展開分析與討論，包括社會生態(tài)系統(tǒng)、可持續(xù)生計研究、人口發(fā)展環(huán)境研究等。在我國，1979 年吳傳鈞提出人地關(guān)系地域系統(tǒng)，1984 年馬世駿等學(xué)者提出社會—經(jīng)濟—自然復(fù)合生態(tài)系統(tǒng)，1999 年方創(chuàng)琳提出人地耦合圈理論。國內(nèi)外相關(guān)研究理論為本論文研究奠定了堅實的理論基礎(chǔ)。

特色小鎮(zhèn)是集“產(chǎn)城人”的空間載體，突破了行政區(qū)劃與園區(qū)的空間限制，注重宜居與宜創(chuàng)融合，工業(yè)化與信息化融合，區(qū)域特色與產(chǎn)業(yè)前沿融合，一二三產(chǎn)業(yè)融合（“四個融合”）。本質(zhì)上看，特色小鎮(zhèn)建設(shè)的初衷是要實現(xiàn)小鎮(zhèn)建設(shè)與生態(tài)環(huán)境的空間耦合。耦合過程是人類與自然之間相互作用、相互影響，并趨于協(xié)調(diào)統(tǒng)一的過程，一般經(jīng)歷無序狀態(tài)、中間狀態(tài)和有序狀態(tài)三個階段。無序狀態(tài)一般是是城鎮(zhèn)化系統(tǒng)包括人口、經(jīng)濟、社會、基礎(chǔ)設(shè)施等子系統(tǒng)與生態(tài)環(huán)境系統(tǒng)包括大氣、土壤、能源、生物等子系統(tǒng)之間雖然在空間與時間上相互作用，但是對立多于統(tǒng)一的狀態(tài)。兩個系統(tǒng)之間通過人流、物質(zhì)流、能流、信息流等中介發(fā)生作用，但是脅迫大于促進作用，生態(tài)環(huán)境惡化就是無序狀態(tài)的表象。中間狀態(tài)是城鎮(zhèn)化系統(tǒng)與生態(tài)環(huán)境系統(tǒng)在時間、空間、組織等多維度多要素的相互作用，不同組織的利益主體之間權(quán)衡經(jīng)濟社會發(fā)展與生態(tài)環(huán)境保護之間的利益。有序狀態(tài)是在長期的多方利益博弈之后，達到一種相互協(xié)調(diào)統(tǒng)一的暫時狀態(tài)，當(dāng)生態(tài)系統(tǒng)或經(jīng)濟社會系統(tǒng)發(fā)展內(nèi)外部干擾，如突發(fā)性自然災(zāi)害、突破性技術(shù)進步、金融危機等，破壞了能量守恒或超出閾值之后，又將進入新一輪的演化過程。實證研究有利于更加準(zhǔn)確地判斷特色小鎮(zhèn)建設(shè)與生態(tài)環(huán)境系統(tǒng)之間關(guān)系處于哪個階段，便于在特色小鎮(zhèn)建設(shè)中查找和解決問題。

三、研究設(shè)計

（一）實證模型的設(shè)定

1.雙重差分法。本文擬采用PSM-DID 方法檢驗特色小鎮(zhèn)建設(shè)對生態(tài)環(huán)境的影響。雙重差分法（Difference-in-Difference，DID）常用于評估政策影響效果，并且能夠較好地克服模型評估中的內(nèi)生性問題。本文的政策影響包括特色小鎮(zhèn)建立前后趨勢變化，同時包括時間效應(yīng)和政策處理效應(yīng)，適用于雙重差分法。特色小鎮(zhèn)建設(shè)與設(shè)立開發(fā)區(qū)等政策相似，同一地區(qū)在不同年份會新增1 個或多個建設(shè)單位，本文所研究的是2016 年、2017 年住房城鄉(xiāng)建設(shè)部分別公布了兩批特色小鎮(zhèn)，第一批127 個特色小鎮(zhèn)，第二批為276 個特色小鎮(zhèn)，這為我們采用雙重差分進行“準(zhǔn)自然實驗”提供良好的基礎(chǔ)。具體來說，在樣本中將2016 年第一批獲得127 個特色小鎮(zhèn)的地級市作為處理組，其余沒有獲批特色小鎮(zhèn)的地級市作為對照組；因2017 年公布特色小鎮(zhèn)與2016 年僅差1 年，時間間隔較短，2017 年將兩批合計中獲得特色小鎮(zhèn)的地級市作為處理組，其余沒有獲批特色小鎮(zhèn)的地級市作為對照組。同時，本文為了更好解釋是否特色小鎮(zhèn)越多的地級市對環(huán)境的影響越顯著，將2017 年累加獲批2 個以上的地級市作為處理組，其余作為對照組進行分析。本文擬先采用單時期的雙重差分法進行測算，公式如下：

其中，下標(biāo)i 和t 分別表示地區(qū)和年份。Yit是生態(tài)環(huán)境變化的指示性變量，在模型中由lnwastew 表示取對數(shù)后的工業(yè)廢水排放量；lnsul 表示取對數(shù)后的工業(yè)二氧化硫排放量；lnsmem表示取對數(shù)后的工業(yè)煙（粉）塵排放量；utea 表示一般工業(yè)固體廢物綜合利用率；wadi 表示污水處理廠集中處理率。模型實證模擬過程中發(fā)現(xiàn)，部分被解釋變量在不同時期檢驗結(jié)果存在較大差異，因此，在分析結(jié)果中將只顯示影響顯著的被解釋變量。Tsxz 表示特色小鎮(zhèn)建設(shè)前后的虛擬指示性變量，tsxza 為2016 年第一批建設(shè)的特色小鎮(zhèn)，tsxzb 為2017 年兩批累計建設(shè)的特色小鎮(zhèn)，tsxzc 為2017 年后兩批共計某城市已有兩個（含）以上特色小鎮(zhèn)建設(shè)。若Tsxz 等于1 為實驗組，代表設(shè)立為特色小鎮(zhèn)的城市；若Tsxz 等于0 為控制組，為沒有設(shè)立特色小鎮(zhèn)的城市。Time 為時間指示性變量，time 等于1 為2016 年之后，等于 0 為 2016 年之前；timee 等于 1 為 2017 年之后，等于 0 為2017 年之前。Xit表示其他控制變量。模型檢驗過程中主要關(guān)注β1的系數(shù)，若β1為正，代表特色小鎮(zhèn)建設(shè)可能會導(dǎo)致污染排放量的增長，或者廢物廢水處理率提升；若β1為負(fù)，代表特色小鎮(zhèn)建設(shè)可能會導(dǎo)致污染排放量的下降，或者廢物廢水處理率降低。

2.傾向得分匹配法（PSM）。DID 方法能夠較好解決政策評估中的內(nèi)生性問題，但是，由于城市眾多差異較大，盡可能選取一批與處理組特征相類似的城市作為對照組，以消除樣本的選擇偏差。具體操作如下：首先計算傾向得分值（weight），建立一個回歸模型，選取評價兩組相似度的指標(biāo)，包括人均GDP（取對數(shù)）、人口數(shù)（取對數(shù)）、地方財政預(yù)算支出（取對數(shù)），因變量設(shè)為虛擬變量，處理組設(shè)為1、對照組設(shè)為0，計算出某城市作為處理組的概率。其次，根據(jù)weight，對每個確定為處理組的城市，從對照組中尋找最相近的城市作為對照組。

3.PSM-DID 方法。PSM 能避免樣本的選擇性偏差，DID 能夠解決內(nèi)生性并評估政策處理效應(yīng)，兩者結(jié)合能夠更加準(zhǔn)確反映政策效果?？刂谱兞縓it 包括：（1）政府支出（lngov）。采用取對數(shù)后的地方財政預(yù)算支出表示，一般政府支出越多，地方政府投資越大，可能影響污染排放量；同時，若政府將部分支出用于污染治理和防治又會促進廢水廢物處理率提升。（2）地區(qū)人口規(guī)模（lnpeople）。地區(qū)人口規(guī)模越大，污染排放量越顯著，本文采用取對數(shù)后的年末總?cè)丝诒硎尽＃?）地區(qū)經(jīng)濟發(fā)展水平（lnpergdp）。一般經(jīng)濟發(fā)展水平越高，產(chǎn)業(yè)的高級化程度越顯著，污染的治理能力較強，本文采用取對數(shù)后的人均國內(nèi)生產(chǎn)總值表示。（4）就業(yè)水平（lnemp）。一般情況而言，就業(yè)水平越高，污染的排放量越顯著，本文采用取對數(shù)的單位從業(yè)人員數(shù)表示。（5）基礎(chǔ)設(shè)施建設(shè)（lnintt）?；A(chǔ)設(shè)施建設(shè)越完善的地區(qū)，環(huán)境污染防治的能力越強，本文采用取對數(shù)后的國際互聯(lián)網(wǎng)用戶數(shù)來表示。

（二）樣本選擇與數(shù)據(jù)來源

本文所使用的2009-2018 年數(shù)據(jù)來源于《中國城市統(tǒng)計年鑒》（2010—2019）。統(tǒng)計年鑒中主要顯示的是291 個城市的相關(guān)數(shù)據(jù)，但是，由于巢湖市、三沙市、海東市的數(shù)據(jù)嚴(yán)重缺失，在模型應(yīng)用時將這3 地剔除，剩余288 個城市10 年約2880 組相關(guān)數(shù)據(jù)。數(shù)據(jù)的統(tǒng)計性描述見表1。

四、實證結(jié)果及分析

（一）平行趨勢檢驗

在使用雙重差分法進行檢驗實驗效果時，必須滿足平行趨勢假設(shè)，也就意味著實驗組和控制組在接受實驗之前具有一致的發(fā)展趨勢。如果實驗組和控制組之間存在系統(tǒng)性差異，直接影響估計結(jié)果的可靠性，因此，為了保證模型的適用性，本文對實驗組和控制組的環(huán)境污染與防治的相關(guān)指標(biāo)進行了同趨勢檢驗。如圖1—圖2 所示，筆者經(jīng)過多次模擬嘗試之后發(fā)現(xiàn)，并非所有的被解釋變量都顯著，只有wadi（污水處理廠集中處理率）、lnwastew（取對數(shù)后的工業(yè)廢水排放量）、lnsul（取對數(shù)后的工業(yè)二氧化硫排放量）顯著（此變量在修正標(biāo)準(zhǔn)誤之后也不顯著），其他解釋變量不顯著，因此，只列出前兩個被解釋變量的平行趨勢。如圖所示，兩個被解釋變量在特色小鎮(zhèn)建設(shè)之前的3~4 年都保持相同趨勢。wadi 圖示在2014 至2016 年趨勢幾乎重合，在2017 年分異后又交叉；lnwastew 圖示2015 至2016 年間趨勢相似，2017 年后從折線圖上看相似，需要通過實證分析進一步驗證差異；lnsul 圖示 2015 至 2017 年間幾乎重合，2017 年之后又出現(xiàn)差異。

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

圖1 wadi 的平行趨勢檢驗

圖2 lnwastew 的平行趨勢檢驗

（二）實證檢驗結(jié)果

在PSM-DID 基準(zhǔn)回歸過程中，首先對288 個城市根據(jù)傾向得分尋找最相近的城市做對照組，然后再用DID 方法測度特色小鎮(zhèn)建設(shè)前后對生態(tài)環(huán)境的影響，在匹配中可能存在無法匹配和單個城市重復(fù)匹配的可能性，所以每組的數(shù)據(jù)值存在一定的差異。

1.特色小鎮(zhèn)建設(shè)對生態(tài)環(huán)境的影響。3 個模型都控制了時間固定效應(yīng)和城市固定效應(yīng)，模型（1）中tsxza×time 的系數(shù)為-2.006，在10%的置信水平顯著，表明在排除其他因素影響后，2016年第一批特色小鎮(zhèn)建設(shè)使得污水處理廠集中處理率下降2.006個單位。為了進一步驗證特色小鎮(zhèn)建設(shè)對生態(tài)影響是正還是負(fù)，筆者嘗試發(fā)現(xiàn)，2016 年第一批特色小鎮(zhèn)建設(shè)對工業(yè)廢水排放量（lnwastew）、工業(yè)二氧化硫排放量（lnsul）、工業(yè)煙（粉）塵排放量（lnsmem）、工業(yè)固體廢物綜合利用率（utea）影響并不顯著，所以在表中并未列出。為了深入分析，因兩批特色小鎮(zhèn)只相隔1 年，用兩批特色小鎮(zhèn)匯總數(shù)據(jù)再進行模型檢驗，模型（2）顯示出特色小鎮(zhèn)建設(shè)對取對數(shù)的工業(yè)廢水排放量（lnwastew）的影響，在10%的置信水平顯著，影響系數(shù)為0.0821，說明兩批合計的特色小鎮(zhèn)建設(shè)增加了工業(yè)廢水排放量，在排除其他因素的影響后，因變量變動1 個單位，工業(yè)廢水排放量增加8.21%。筆者還嘗試兩批特色小鎮(zhèn)中超過2 個特色小鎮(zhèn)的城市篩選統(tǒng)計入模型（3），發(fā)現(xiàn)特色小鎮(zhèn)建設(shè)數(shù)量的增加對生態(tài)環(huán)境相關(guān)指標(biāo)的影響并不顯著。

2.控制變量的影響。人均GDP（lnpergdp）對各項指標(biāo)的影響是正向的，一方面經(jīng)濟發(fā)展較快的地區(qū)的排污量可能增加，另一方面治理生態(tài)環(huán)境的投入和成本也會增加。例如，在模型（2）和模型（3）中l(wèi)npergdp 顯著為正，說明經(jīng)濟發(fā)展水平越高，污染排放量越高；在模型（1）中系數(shù)顯著為正，說明經(jīng)濟發(fā)展水平越高的地區(qū)污染治理能力越強。地方財政預(yù)算支出（lngov）的增加，會加大對環(huán)境治理的投入力度，使得工業(yè)廢水排放量和工業(yè)二氧化硫的排放量等指標(biāo)下降。人口數(shù)量增加或變動會影響工業(yè)廢水排放量增加或變動?；ヂ?lián)網(wǎng)等基礎(chǔ)設(shè)施建設(shè)一方面會提高污水集中處理率，另一方面基礎(chǔ)設(shè)施投入也會增加污染排放量。從控制變量的結(jié)果可以看出，經(jīng)濟的快速發(fā)展會導(dǎo)致污染排放量增加，同時，也增強了污染治理能力，呈現(xiàn)邊發(fā)展邊治理的趨勢。

表2 特色小鎮(zhèn)建設(shè)對生態(tài)環(huán)境影響的PSM-DID 基準(zhǔn)回歸結(jié)果

3.穩(wěn)健性檢驗。我國東中西部地區(qū)經(jīng)濟發(fā)展水平、資源稟賦、經(jīng)濟地理條件存在巨大差異，同樣的政策是否在不同的地區(qū)效果存在差異？在表3 中呈現(xiàn)的是特色小鎮(zhèn)建設(shè)對東中西部地區(qū)工業(yè)廢水排放量影響的穩(wěn)健性檢驗，結(jié)果表明：特色小鎮(zhèn)建設(shè)尤其對東部地區(qū)的工業(yè)廢水排放量影響顯著，結(jié)果與表2 所示是一致的，東部地區(qū)特色小鎮(zhèn)建設(shè)對工業(yè)廢水排放量的影響最大。

表3 特色小鎮(zhèn)建設(shè)對東中西部地區(qū)工業(yè)廢水排放量的影響

五、結(jié)論與建議

由實證檢驗結(jié)果可知，特色小鎮(zhèn)建設(shè)與生態(tài)環(huán)境之間的關(guān)系還處于無序狀態(tài)，可能會導(dǎo)致污染排放量增加，地方政府面臨更大的污染治理壓力。同時，經(jīng)濟發(fā)展水平越高的地區(qū)，污染排放量越多。主要原因有以下兩點：一是特色小鎮(zhèn)都有明確的主導(dǎo)產(chǎn)業(yè)如裝備制造、科創(chuàng)等，產(chǎn)業(yè)的培育與快速發(fā)展期會占用大量的資源，可能導(dǎo)致特色小鎮(zhèn)建設(shè)帶來的污染排放量增加。二是特色小鎮(zhèn)建設(shè)是城鎮(zhèn)多元功能的聚合，在完善配套設(shè)施、拓展發(fā)展空間過程中，可能會出現(xiàn)小鎮(zhèn)建設(shè)與環(huán)境保護之間的沖突。特色小鎮(zhèn)在全國推進是從2016 年開始，本文中考察時間范圍到2018年，正處于特色小鎮(zhèn)快速發(fā)展期，難免會出現(xiàn)環(huán)境沖突。為促進特色小鎮(zhèn)高質(zhì)量發(fā)展，應(yīng)充分把握特色小鎮(zhèn)建設(shè)內(nèi)涵，科學(xué)定位特色小鎮(zhèn)主導(dǎo)產(chǎn)業(yè)，推動產(chǎn)業(yè)綠色化和高級化發(fā)展；科學(xué)規(guī)劃利用城鎮(zhèn)空間，因地制宜優(yōu)化配套設(shè)施布局，保護修復(fù)生態(tài)環(huán)境，打造宜居宜業(yè)的特色小鎮(zhèn)。

注釋：

①方創(chuàng)琳，楊玉梅.城市化與生態(tài)環(huán)境交互耦合系統(tǒng)的基本定律[J].干旱區(qū)地理，2006（02）:1

②劉海猛，方創(chuàng)琳，李詠紅.城鎮(zhèn)化與生態(tài)環(huán)境“耦合魔方”的基本概念及框架[J].地理學(xué)報，2019（08）:1489

③梁振民，陳才.城市化與生態(tài)環(huán)境非協(xié)調(diào)性耦合關(guān)系識別研究[J].安徽大學(xué)學(xué)報（自然科學(xué)版），2019（01）：72

④劉耀彬，宋學(xué)鋒.城市化與生態(tài)環(huán)境耦合模式及判別[J].地理科學(xué)，2005（08）：408

⑤唐志強.城鄉(xiāng)融合視野下農(nóng)村轉(zhuǎn)移人口質(zhì)量與轉(zhuǎn)移地區(qū)生態(tài)環(huán)境承載力耦合研究[J].2019（04）:70

⑥王小章.特色小鎮(zhèn)的“特色”與“一般”[J].浙江社會科學(xué)，2016（03）：46

⑦衛(wèi)龍寶，史新杰.浙江特色小鎮(zhèn)建設(shè)的若干思考與建議[J].浙江社會科學(xué)，2016（03）：28

⑧盛世豪，張偉明.特色小鎮(zhèn)：一種產(chǎn)業(yè)空間組織形式[J].浙江社會科學(xué)，2016（03）：36

⑨郁建興，張蔚文，高翔.浙江省特色小鎮(zhèn)建設(shè)的基本經(jīng)驗與未來[J].浙江社會科學(xué)，2017（06）：143

⑩王博雅，張車偉，蔡翼飛.特色小鎮(zhèn)的定位與功能再認(rèn)識——城鄉(xiāng)融合發(fā)展的重要載體[J].2020（01）：143

?易開剛，厲飛芹.特色小鎮(zhèn)的功能聚合與戰(zhàn)略選擇——基于浙江的考察[J].經(jīng)濟理論與經(jīng)濟管理，2019（08）：104

?莊晉財，盧文秀，華賢宇.產(chǎn)業(yè)鏈空間分置與特色小鎮(zhèn)產(chǎn)業(yè)培育[J].學(xué)習(xí)與實踐，2018（08）：36

?時浩楠，楊雪云.國家級特色小鎮(zhèn)空間分布特征[J].干旱區(qū)資源與環(huán)境，2019（03）：39

?王長松，賈世奇.中國特色小鎮(zhèn)的特色指標(biāo)體系與評價[J].南京社會科學(xué)，2019（02）：79