新能源示范城市政策對企業(yè)綠色創(chuàng)新影響研究
——基于我國A股上市公司樣本的PSM-DID分析

彭 順

(重慶工商大學 長江上游經(jīng)濟研究中心，重慶 400067)

一、引言

為應對日益凸顯的環(huán)境污染和國際碳減排壓力，我國開始放棄以“高能耗、高污染、高排放”為特征的傳統(tǒng)發(fā)展方式，積極探索可再生能源的創(chuàng)新發(fā)展模式。自2007年始，我國陸續(xù)出臺了《可再生能源中長期發(fā)展規(guī)劃》《可再生能源發(fā)展“十二五”規(guī)劃》，對新能源發(fā)展進行了部署。2012年國家能源局發(fā)布了《關于申報新能源示范城市和產(chǎn)業(yè)園區(qū)的通知》，正式啟動新能源示范城市和產(chǎn)業(yè)園區(qū)的申報工作。繼而2014年國家能源局發(fā)布《關于公布創(chuàng)建新能源示范城市(產(chǎn)業(yè)園區(qū))名單(第一批)的通知》，確定了81個城市和8個產(chǎn)業(yè)園區(qū)為第一批創(chuàng)建新能源示范城市和產(chǎn)業(yè)園區(qū)。

企業(yè)作為實體單位，在新能源示范城市建設過程中發(fā)揮著重要作用。企業(yè)在新能源發(fā)電、建造5G基站、組建大數(shù)據(jù)中心等過程中的創(chuàng)新研發(fā)活動無疑為新能源高質(zhì)量發(fā)展和城市能源轉型注入了強大的動力。而能源轉型不是一蹴而就的，為穩(wěn)步發(fā)展、穩(wěn)中向前，需要對新能源示范城市可再生能源發(fā)展政策(以下簡稱新能源示范城市政策)影響下的企業(yè)綠色創(chuàng)新潛力進行深入評估，以加強企業(yè)綠色創(chuàng)新與新能源示范城市政策的有效銜接，促進可再生能源健康持續(xù)發(fā)展。

在此背景下，本文擬在新能源示范城市這一“準自然實驗”條件下，采用傾向得分匹配模型(PSM)與雙重差分法(DID)相結合的方法評估企業(yè)的綠色創(chuàng)新潛力，即研究新能源示范城市政策對企業(yè)綠色創(chuàng)新影響的大小、機制和時效，以及異質(zhì)性企業(yè)對政策的處置效應。本研究對于新能源示范城市進一步建設和發(fā)展具有一定的指導意義，同時為后續(xù)新能源示范城市的政策制定提供借鑒與參考。

二、文獻綜述

可再生能源發(fā)展政策屬于一種環(huán)境規(guī)制政策。環(huán)境規(guī)制最早來源于經(jīng)濟學原理的市場失靈外部性理論。新古典學派代表Pigou[1]將“外部經(jīng)濟”進一步擴展為“外部經(jīng)濟”與“外部不經(jīng)濟”，提出了“庇古稅”，即通過補貼或征稅來應對環(huán)境外部性，此為政府干預型環(huán)境規(guī)制；Eskeland和Jimenez[2]將政府的環(huán)境規(guī)制進一步劃分為市場激勵型、命令控制型以及資源協(xié)議或環(huán)境認證等形式的“軟”手段監(jiān)管型。此后對環(huán)境規(guī)制的分類研究則繞不開這三種類型。相對于命令控制型，市場激勵型環(huán)境規(guī)制往往被視為更加靈活的監(jiān)管方式[3-4]，具有以下幾種優(yōu)勢：一是具有市場化信息成本節(jié)約優(yōu)勢與長期激勵減污技術發(fā)展優(yōu)勢[5]；二是使企業(yè)的創(chuàng)新靈活性增加、生產(chǎn)力提高[6]；三是使企業(yè)和產(chǎn)業(yè)的邊際減排成本相等[7]。新能源示范城市政策兼有市場激勵(自愿申請)與命令控制(排放約束)兩種類型環(huán)境規(guī)制的特點，而企業(yè)對環(huán)境監(jiān)管的主動自愿行為則屬于“軟”手段監(jiān)管型。

近些年國內(nèi)外學者熱衷于研究環(huán)境規(guī)制對企業(yè)綠色創(chuàng)新的影響。Johnstone等[8]研究發(fā)現(xiàn)，可交易能源許可證制度有助于誘發(fā)與化石能源直接競爭的行業(yè)的技術創(chuàng)新；Calel和Dechezleprêtre[9]檢驗了歐盟排放交易體系對企業(yè)技術變化的影響，結果發(fā)現(xiàn)對企業(yè)低碳技術的促進作用水平大約為10%；Weiss和Anisimova[10]對瑞典高污染行業(yè)企業(yè)開展實證檢驗，結果顯示，靈活的命令控制型環(huán)境規(guī)制和經(jīng)濟激勵政策工具能夠通過提高能源使用效率顯著產(chǎn)生創(chuàng)新補償效應。錢麗等[11]的研究表明，資源環(huán)境約束能夠正向影響中國企業(yè)綠色創(chuàng)新能力；于連超等[12]發(fā)現(xiàn)，環(huán)境稅能夠顯著促進企業(yè)綠色創(chuàng)新，并呈現(xiàn)滯后性特征；李青原和肖澤華[13]發(fā)現(xiàn)，排污收費“倒逼”了企業(yè)綠色創(chuàng)新能力，“倒逼”效應體現(xiàn)在外部壓力和內(nèi)部激勵上。

基于政策試點的示范城市類環(huán)境規(guī)制對企業(yè)綠色創(chuàng)新影響研究的文獻也較為豐富。鄧玉萍等[14]研究發(fā)現(xiàn)，《萬家企業(yè)節(jié)能低碳行動實施方案》對企業(yè)綠色創(chuàng)新有著積極的影響；徐佳和崔靜波[15]研究發(fā)現(xiàn)，低碳城市政策能夠有效促進高碳行業(yè)中國有企業(yè)和非國有企業(yè)的綠色技術創(chuàng)新；鐘昌標等[16]評估了低碳試點政策對企業(yè)綠色創(chuàng)新的積極促進作用；曹廷求等[17]研究發(fā)現(xiàn)，《綠色信貸指引》的實施對企業(yè)綠色創(chuàng)新活動具有正向影響，能夠促進重污染企業(yè)積極進行綠色轉型；宋德勇等[18]評估了智慧城市建設促進綠色技術創(chuàng)新的顯著程度；金環(huán)等[19]研究發(fā)現(xiàn)，電商示范城市建設試點能夠顯著促進企業(yè)綠色技術創(chuàng)新，而且對民營企業(yè)綠色技術創(chuàng)新的作用更顯著。

通過眾多文獻發(fā)現(xiàn)，驗證“波特假說”(1)波特假說：環(huán)境規(guī)制可能會促使企業(yè)進行更多的創(chuàng)新活動以減少污染，而這些創(chuàng)新活動會間接地促進企業(yè)生產(chǎn)力的提高，從而抵消由環(huán)境規(guī)制而增加的成本，并且提升企業(yè)在市場上的盈利能力和相對競爭力。中環(huán)境規(guī)制對企業(yè)綠色創(chuàng)新的影響近些年已成為研究焦點，對政策試點的深入研究也成為學者新的探索領域，尤其是政策試點與企業(yè)綠色創(chuàng)新的關系更是成為研究重點。但是這些研究在使用模型時多以控制變量和增減變量來探討政策試點的處理效應，未從縱向角度研究其動態(tài)效應，也未能深度挖掘效應的持續(xù)時間。而且，綜觀近年文獻，尚未發(fā)現(xiàn)針對新能源示范城市這一政策試點與企業(yè)綠色創(chuàng)新關系的研究。基于此，本文將彌補以上不足，主要貢獻如下：第一，采用面板傾向得分匹配模型(PSM)結合雙重差分法(DID)，驗證新能源示范城市政策試點能否誘發(fā)企業(yè)綠色創(chuàng)新；第二，以政策時點的縱向視角觀察新能源示范城市政策的處理效應，構建政策實施后三期的動態(tài)效應模型，以識別新能源示范城市政策對每一期企業(yè)綠色創(chuàng)新的效應以及持續(xù)期，此處理方法有助于更深入發(fā)掘新能源示范城市政策對企業(yè)綠色創(chuàng)新的時點縱向影響及相互作用規(guī)律；第三，創(chuàng)造性地運用面板分位數(shù)結合雙重差分法評估企業(yè)綠色創(chuàng)新，以數(shù)值橫向水平估計企業(yè)綠色創(chuàng)新能力，用以考察新能源示范城市政策與企業(yè)綠色創(chuàng)新之間的定量關系；第四，從所有制、所屬行業(yè)、四大審計的企業(yè)異質(zhì)性角度，探索新能源示范城市政策對企業(yè)綠色創(chuàng)新的影響效應及機制。

三、研究模型與樣本

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

本文以2003-2019年A股上市公司為研究對象，其中，剔除了金融保險業(yè)上市公司、非正常交易上市公司(包括ST，ST*以及PT)以及相關數(shù)據(jù)缺失的上市公司。本文數(shù)據(jù)主要包括兩個部分：一是企業(yè)創(chuàng)新數(shù)據(jù)，其中公司發(fā)明專利和實用新型專利的專利分類號信息從中國研究數(shù)據(jù)服務平臺(CNRDS)獲??；二是公司其他特征數(shù)據(jù)，其中財務數(shù)據(jù)主要來源于國泰安數(shù)據(jù)庫。將上述數(shù)據(jù)匹配后，最終獲得11 135個年度觀測值。

關于新能源示范城市內(nèi)企業(yè)的識別。借鑒向?qū)捇⒑完戙慬20]識別開發(fā)區(qū)內(nèi)企業(yè)的方法，即如果企業(yè)地址信息顯示企業(yè)位于新能源示范城市內(nèi)，則將該企業(yè)識別為新能源示范城市企業(yè)；反之，則為非新能源示范城市企業(yè)。

關于政策時點的識別。由于新能源示范城市的創(chuàng)建是2012年申報、2014年公布名單，無法獲取每個城市準確的準備時間，只能選定一個統(tǒng)一的政策時點。因此，借鑒評估政策效應的雙重差分法的相關文獻以及研究新能源示范城市的文獻[21-23]，選取政策時點為公布名單的2014年。

(二)PSM樣本選取

Rosenbaum和Rubin[24-25]經(jīng)研究提出了一種通過得分匹配獲得“反事實”的估計方法，這種方法直接促成了傾向匹配得分模型(PSM)的產(chǎn)生。在該模型的應用中，首先是盡可能多地找到實驗組與對照組之間存在差異的變量，組成匹配變量集，記為X≡(x1，x2，…，xn)，然后將分組變量Treated(0-1)對變量集X進行參數(shù)回歸，常用回歸方式為logit(probit)，具體如式(1)所示：

logit(Treated)=Xβ+ε。

(1)

采用式(1)可得回歸系數(shù)矩陣β。

通過式(2)可計算出得分值ps。

(2)

式中，Xi代表個體i的匹配變量集。

式(2)的本質(zhì)是將多維的匹配變量擬合成一個一維且介于[0，1]之間的變量psi，即個體i進入實驗組的概率。得到psi后，可以根據(jù)特定的匹配方式，將實驗組與對照組匹配在一起，最后將實驗組與對照組進行對比，以識別出二者之間的差異。

但是該得分方法針對的是截面數(shù)據(jù)的個體，對于面板數(shù)據(jù)而言，忽略時間效應顯然是不可取的。由此，石大千等[26]在研究中采用了逐期匹配的方法，先將面板數(shù)據(jù)拆分為n>1期截面數(shù)據(jù)，依據(jù)式(1)和式(2)對每期截面數(shù)據(jù)中所有匹配變量進行逐期logit(probit)回歸，計算出每個個體每期的psi，t；然后分別對每期的psi，t進行匹配，得到每期的匹配樣本；最后將每期的樣本合并成面板數(shù)據(jù)。本文將此方法記為PSM-DID。謝申祥等[27]在石大千等人的研究基礎上將PSM進一步優(yōu)化，即將每期截面數(shù)據(jù)的得分值用等時間權重加權獲得一個綜合得分值，根據(jù)這個綜合得分值匹配出合適的個體樣本，把這些樣本按時間序列組合形成面板數(shù)據(jù)。本文將此方法記為PSMs-DID。另外，該方法還對特殊變量進行了處理，因為在匹配前后特殊變量必須滿足一致性，即同一行業(yè)相匹配或者同一企業(yè)性質(zhì)相匹配。

參考Chen等[28]的研究，確定公司規(guī)模(Size)、資產(chǎn)收益率(ROA)、資產(chǎn)負債率(Lev)、總資產(chǎn)周轉率(ATO)、現(xiàn)金流比率(Cashflow)、托賓Q值(TobinQ)、公司成立年限(FirmAge)為PSM方法的匹配變量。

(三)模型設計

為評估新能源示范城市政策對企業(yè)綠色創(chuàng)新的影響，對國家能源局2014年發(fā)布的新能源示范城市做一次準自然實驗。通過構建試點城市與非試點城市、政策前與政策后的交互項，剔除一些不隨個體變化或時間變化的因素，以合理評估新能源示范城市的政策效應。借鑒陳林和伍海軍[29]關于雙重差分法的研究成果與建議，本文建立如下計量模型：

Yi，t=α+βTreatedi×Timet+∑ρControli，t+λi+σt+ηc+εi，t。

(3)

式中，Yi，t表示試點城市某企業(yè)在某年的綠色創(chuàng)新變量，以授權專利等數(shù)據(jù)表征；Treatedi為試點城市虛擬變量，若某企業(yè)為新能源示范城市企業(yè)，為“1”，反之為“0”；Timet為新能源示范城市政策實施時間虛擬變量，2014年之后為“1”，之前為“0”；Treatedi×Timet是Treatedi和Timet的交互項，也稱did虛擬變量，其系數(shù)β正是本文最關心的政策效應；Controli，t是選取的一組控制變量；λi，σt和ηc代表個體、時間以及行業(yè)截距項的固定效應；εi，t表示計量模型隨機擾動項。

(四)變量描述性分析

表1為各變量描述性統(tǒng)計結果。本文被解釋變量為企業(yè)綠色專利授權情況，包括企業(yè)綠色專利授權量(GrePatNum)與企業(yè)綠色專利授權比率(GrePatRat)，以及企業(yè)綠色專利授權量(GrePatNum)的替代變量——綠色發(fā)明專利授權量(InPatNum)和綠色實用新型專利授權量(UmPatNum)，而且對企業(yè)綠色專利授權量作1+N的對數(shù)處理?？刂谱兞繛楣疽?guī)模(Size)、資產(chǎn)收益率(ROA)、資產(chǎn)負債率(Lev)、總資產(chǎn)周轉率(ATO)、現(xiàn)金流比率(Cashflow)、托賓Q值(TobinQ)、公司成立年限(FirmAge)。其中，公司規(guī)模(Size)參照規(guī)模經(jīng)濟理論[30-32]選?。还境闪⒛晗?FirmAge)參照組織慣性理論[33]選?。黄渌?個變量是財務報表中的重要指標，屬于被廣泛接受的企業(yè)績效衡量標準，代表企業(yè)市場價值和資產(chǎn)重置價值之間的關系。依據(jù)企業(yè)特征并結合樣本，選取特殊變量為企業(yè)的所有制(SOE)、所屬行業(yè)(Industry)、四大審計(Big4)，即以是否為國有企業(yè)、是否屬高能耗行業(yè)、是否請四大審計公司進行審計來代表企業(yè)的異質(zhì)性。

表1 各變量描述性統(tǒng)計結果

四、實證結果分析

(一)政策的唯一性檢驗

要證明新能源示范城市政策對企業(yè)綠色創(chuàng)新具有影響力，根據(jù)“反事實”的研究框架，一定要排除掉“非政策時點”與“非政策試點個體”的干擾，由此，必須對政策時點的唯一性與政策試點個體的唯一性逐個進行檢驗。

1.政策時點的唯一性——平行趨勢檢驗

從政策時點上考量，有效應用雙重差分法的重要假設是，如果沒有新能源示范城市政策的影響，隨著時間推移，實驗組和對照組的企業(yè)綠色創(chuàng)新變化趨勢不具有系統(tǒng)性差異。因此，需要檢驗平行趨勢或者共同趨勢的假設是否成立，即檢驗實驗組和對照組的企業(yè)綠色創(chuàng)新在政策實施之前必須具有相同的發(fā)展趨勢。如果不滿足這一條件，那么雙重差分法得出的政策效應β則不完全是真實的政策效應。

企業(yè)綠色專利授權情況的平行趨勢如圖1所示，橫軸表示年度，圖1(a)縱軸表示加總到城市層面的企業(yè)綠色專利授權量的平均值，圖1(b)縱軸則表示企業(yè)綠色專利授權比率的平均值，以新能源示范城市名單發(fā)布的2014年為政策前與政策后的分界線。在政策前(2003-2014年)，隨著年度的變化，對于企業(yè)綠色專利授權量均值，其實驗組總體略高于對照組，兩者的波動大致呈平行趨勢；實驗組的企業(yè)綠色專利授權比率均值與對照組呈現(xiàn)“你追我趕”的波動趨勢。在政策后(2014-2018年)，實驗組的企業(yè)綠色專利授權量均值攀升速度明顯加快，與對照組之間的差距逐步拉大；同樣，對于綠色專利授權比率均值，其實驗組與對照組之間的差距也增大，并且再未出現(xiàn)被“超越”現(xiàn)象。因此，在政策時點2014年之后，確實使原本具有平行趨勢的實驗組與對照組之間的差距逐漸拉大，雙重差分法的平行趨勢假說從圖形上得到支持。

(a)

(b)圖1 平行趨勢檢驗-時間趨勢圖

2.政策試點個體的唯一性——安慰劑檢驗

從政策試點個體上考量，需要證明隨著時間的推移，沒有依據(jù)新能源示范城市名單選擇的個體不會出現(xiàn)顯著性變化。為避免出現(xiàn)人為的設定，本文采取隨機抽樣的方式進行實驗組的選取。對抽取的樣本按照公式(3)進行回歸，得到系數(shù)β，t值與置信區(qū)間；依此重復多次后，得到隨機抽樣的估計值區(qū)間及其分布；最后與原模型的估計值比較差異。這種檢驗政策試點個體唯一性的方法稱為安慰劑檢驗法。

原模型的估計系數(shù)為0.049，t值為2.45。采用1 000次隨機抽樣，安慰劑檢驗結果如圖2所示。隨機數(shù)采用的是正態(tài)分布函數(shù)，大樣本下系數(shù)均值為0，90%水平下的置信區(qū)間為[-1.65，1.65]。從圖2中可以看出，實際試點城市樣本下的估計值明顯異于隨機抽樣試點個體的結果，原模型穩(wěn)健性的結果說明政策試點個體具有唯一性，即新能源示范城市名單中試點城市的唯一性得到支持。

(a)

(b)

(c)圖2 安慰劑檢驗結果

(二)企業(yè)綠色創(chuàng)新的回歸結果分析

表2為企業(yè)綠色創(chuàng)新的回歸結果。第1-2欄為沒有進行PSM匹配的回歸結果，全樣本量為11 135；第3-4欄為參照了石大千等[26]逐期PSM匹配方法的回歸結果，通過近鄰1∶3匹配將樣本量由11 135篩選至8 573；第5-6欄為借鑒了謝申祥等[27]使用的方法先獲得逐期得分再以時間權重進行匹配的回歸結果，通過近鄰1∶3匹配將樣本量由11 135篩選至8 409。表2中的模型均采用了Gaure[34]，Guimaraes和Portugal[35]提出的固定效應的回歸方法，將個體、時間、行業(yè)作為固定截距項。

模型回歸結果顯示，綠色專利授權量與綠色專利授權比率在DID，PSM-DID，PSMs-DID三個模型中的did系數(shù)均大于0，且至少在10%顯著性水平上顯著。這說明無論是以綠色專利授權量(GrePatNum)衡量的綠色創(chuàng)新絕對水平，還是以綠色專利授權比率(GrePatRat)衡量的綠色創(chuàng)新相對水平，新能源示范城市政策均對試點城市企業(yè)綠色創(chuàng)新產(chǎn)生了積極的推動作用。

表2中綠色創(chuàng)新絕對水平模型回歸的did系數(shù)分別為0.049，0.039，0.040，顯著性水平依次為1%，10%，5%，可見PSMs-DID回歸結果的大小處于另外兩個結果之間。而通過傾向得分匹配方法可知，PSMs-DID將平衡面板數(shù)據(jù)匹配后依然是平衡面板數(shù)據(jù)，而PSM-DID匹配后的數(shù)據(jù)不一定是平衡面板數(shù)據(jù)，可見PSMs-DID的理論機制準確且結果穩(wěn)健。因此在后文的回歸分析中，主要采用PSMs-DID模型。

表2 企業(yè)綠色創(chuàng)新的回歸結果

由表2可知，綠色創(chuàng)新相對水平的回歸結果均比綠色創(chuàng)新絕對水平更加穩(wěn)健，平均顯著性水平更高。但是由于知識產(chǎn)權在最近幾年得到了進一步保障，專利授權量呈現(xiàn)高速增長趨勢，如果用相對水平去衡量綠色創(chuàng)新程度不免會遺漏綠色專利授權的基數(shù)及其增加值。而代表綠色創(chuàng)新絕對水平的綠色專利授權量卻切實地用到了環(huán)境污染的治理中，相較而言，解釋力度更高。因此下文以綠色專利授權量(GrePatNum)為主要被解釋變量。

為進一步保證回歸結果的穩(wěn)健性，需要考察新能源示范城市政策誘發(fā)上市公司企業(yè)綠色創(chuàng)新活動的類別。本文以綠色發(fā)明專利授權量(InPatNum)和綠色實用新型專利授權量(UmPatNum)作為企業(yè)綠色專利授權量的替代變量，來驗證不同專利類別的綠色創(chuàng)新水平。為消除異方差，對這兩種綠色專利類別變量進行對數(shù)處理。

表3為兩種綠色專利類別變量的回歸結果，其中模型的回歸參照了表2中的PSM處理方式與固定個體、時間、行業(yè)的回歸方法。觀察表3可知，在三個模型中綠色發(fā)明專利授權量(InPatNum)的did系數(shù)都不顯著，而綠色實用新型專利授權量(UmPatNum)的did系數(shù)都顯著。綠色發(fā)明專利授權量的did系數(shù)在0.01左右，綠色實用新型專利授權量的did系數(shù)在0.05左右，比較兩者系數(shù)可知，綠色實用新型專利授權量比綠色發(fā)明專利授權量更受政策的影響。綠色實用新型專利具有研發(fā)周期短、消耗資源少和實用性強等特點[36]，而綠色發(fā)明專利比綠色實用新型專利的研發(fā)周期要長。本文觀測的年份為2003-2019年，從新能源示范城市開始創(chuàng)建的2014年至2019年，僅有5年時間，綠色發(fā)明專利授權量不顯著也能被解釋。由此得出結論，新能源示范城市創(chuàng)建后的5年中，新能源示范城市政策主要影響了綠色實用新型專利授權量。

表3 兩種綠色專利類別變量的回歸結果

(三)政策時點對企業(yè)綠色創(chuàng)新的動態(tài)效應分析

從縱向維度政策時點上繼續(xù)考察其對企業(yè)綠色創(chuàng)新的動態(tài)效應。參照蔣靈多和陸毅[37]的做法，納入政策時點2014年前的-3，-2，-1期以及2014年后的+1，+2，+3，+4期，生成7個時點虛擬變量，再與Treated變量相乘，生成did虛擬變量，然后將去掉的2014年的時點虛擬變量即0期作為政策效應0值的對照，最后在加入控制變量和固定個體、時間、行業(yè)的條件下，分別將企業(yè)綠色專利授權量和企業(yè)綠色專利授權比率對所有生成的did虛擬變量進行回歸。根據(jù)回歸結果中的did虛擬變量系數(shù)畫出政策時點的動態(tài)效應圖，即以直觀圖的方式呈現(xiàn)新能源示范城市政策在不同年度對企業(yè)綠色創(chuàng)新的動態(tài)效應，如圖3所示。

(a)

(b)注：實線為did虛擬變量系數(shù)，虛線為系數(shù)的上下置信區(qū)間圖3 政策時點的動態(tài)效應圖

從圖3(a)可以看出，對于企業(yè)綠色專利授權量，其交互項的系數(shù)在政策時點2014年前-3，-2期時幾乎為0，-1期時雖不接近于0值，但為負數(shù)，說明在政策時點前實驗組和對照組之間不存在顯著差異，這進一步滿足了平行趨勢的假設；而在政策時點2014年后+1期時交互項的系數(shù)顯著為正，隨后在+2，+3，+4期時一直保持為正，說明新能源示范城市政策在實施后產(chǎn)生了一個長期的正效應，企業(yè)的綠色創(chuàng)新受其長期影響。由圖3(b)可以看出，以企業(yè)綠色專利授權比率為被解釋變量的政策時點動態(tài)效果與圖3(a)大致相同，并且都表明政策實施后交互項的系數(shù)在+2期時開始降低，由此便產(chǎn)生一個疑問：新能源示范城市政策效應的持續(xù)期有多長？

為了解答上述疑問，本文借鑒Bernile等[38]研究的沖擊的動態(tài)效應方法。該研究在分析多期災難對基金波動性的影響時，對災難發(fā)生前1年與災難發(fā)生后第2年、災難發(fā)生前1年與災難發(fā)生后第3年、災難發(fā)生前1年與災難發(fā)生后第4年三個子樣本分別作了多期DID回歸，以觀察災難發(fā)生對基金波動性的動態(tài)效應。依據(jù)這種樣本處理方法，本文選取2013年與2015年的樣本、2013年與2016年的樣本、2013年與2017年的樣本，分別使用企業(yè)綠色專利授權量與企業(yè)綠色專利授權比率進行回歸，結果如表4所示，其中2013年的樣本為政策時點的對照。

表4 政策對企業(yè)綠色創(chuàng)新的動態(tài)效應回歸結果

由表4可知，政策時點2014年后的企業(yè)綠色專利授權量和企業(yè)綠色專利授權比率的did系數(shù)逐年降低，由政策第2年(2015年)的0.152，0.247下降到政策第3年(2016年)的0.117，0.210，再下降到政策第4年(2017年)的幾乎為0且不顯著，與對照年份2013年的企業(yè)綠色創(chuàng)新水平差異不明顯，可見政策對企業(yè)綠色創(chuàng)新的影響效應隨著時間推移不斷降低，并且在第4年幾乎消失，說明政策效應的持續(xù)期大致為3年。

(四)企業(yè)綠色專利授權量的分位數(shù)回歸

在研究政策效應縱向時點維度的基礎上，本文對政策效應的橫向維度即綠色專利授權量的大小進行研究。不同水平的綠色專利授權量對新能源示范城市政策有著怎樣的反應？本文通過構建面板分位數(shù)的計量模型并采用Machado和Silva[39]提出的面板分位數(shù)估計方法來解釋這個問題，模型如下：

Qv(τ|xi.t)=αi+β(τ)Treatedi×Timet+∑ρ(τ)Controli.t+εi.t。

(4)

式中，Qv(τ|xi.t)為綠色專利授權量的τ分位數(shù)，αi為截距項，β(τ)為did的τ分位的系數(shù)，ρ(τ)為控制變量的τ分位的系數(shù)，εi.t為隨機擾動項。

在處理樣本時，由于綠色專利授權量的0值較多，故采用嚴格的1∶1近鄰逐期PSM方法對綠色專利授權量為0值的樣本進行匹配，最終獲得0值的樣本量為3 000，大于0值的樣本量為1 700。顯然這些數(shù)據(jù)不能獲得全分位的分位數(shù)，只能從0.6分位開始，每0.1作為一個分位點。

不同綠色專利授權量水平下的企業(yè)綠色創(chuàng)新面板分位數(shù)回歸結果如表5所示。在控制變量、個體固定的條件下，企業(yè)綠色專利授權量的did系數(shù)從0.6分位數(shù)至0.9分位數(shù)逐漸增加，并且顯著性水平也有逐漸顯著的趨勢。這進一步說明了回歸結果的穩(wěn)健性，即在不同分位數(shù)下，did系數(shù)的正負、政策效應的影響方向均未發(fā)生變化。

表5 企業(yè)綠色創(chuàng)新面板分位數(shù)回歸結果

為了使面板分位數(shù)在采用分位數(shù)回歸方法時，自變量分組數(shù)值能夠跟隨因變量分組而發(fā)生變化，故將政策虛擬變量did與企業(yè)變量交叉相乘組成交互項，來檢驗政策對企業(yè)變量的調(diào)節(jié)效應。選取企業(yè)規(guī)模(Size)、現(xiàn)金流比率(Cashflow)、公司成立年限(FirmAge)分別代表企業(yè)的發(fā)展屬性、企業(yè)的盈利屬性、企業(yè)成長屬性并與政策虛擬變量did組成交互項，對企業(yè)綠色專利授權量的分位數(shù)進行回歸。

did調(diào)節(jié)效應的面板分位數(shù)回歸結果如表6所示。政策虛擬變量did與企業(yè)規(guī)模(Size)的交互項代表了受到政策沖擊的企業(yè)其規(guī)模對自身綠色創(chuàng)新的影響。觀察系數(shù)可知，在不同綠色專利授權量分位數(shù)下，did×Size的系數(shù)保持在0.113左右，可見企業(yè)規(guī)模對企業(yè)綠色創(chuàng)新影響非常穩(wěn)定，說明對于不同的綠色專利授權量分位數(shù)，政策對不同規(guī)模的企業(yè)的綠色創(chuàng)新影響程度基本相似。在三組面板分位數(shù)回歸結果中did×Cashflow的系數(shù)數(shù)值最大，大于1，且從0.6分位數(shù)至0.9分位數(shù)逐步上升，由1.177上升至1.406，這說明在不同綠色專利授權量分位數(shù)下，現(xiàn)金流比率的綠色創(chuàng)新邊際效應非常明顯，且顯著提升。而did×FirmAge的系數(shù)在此模型中是負向顯著，并且政策試點內(nèi)企業(yè)綠色專利授權量分位數(shù)增加，公司成立年限對綠色創(chuàng)新抑制作用卻隨之下降?？赡艿脑蚴牵簩τ?014年后成立的公司，國家對其綠色發(fā)展有更高的審核要求，樣本本身的偏差導致了負向作用，而當變量達到一定水平，企業(yè)年限的正向作用對樣本偏差的負向作用會有一定程度的緩解。

表6 did調(diào)節(jié)效應的面板分位數(shù)回歸結果

五、異質(zhì)性分析

因為表3涉及的三個回歸模型中PSMs-DID的回歸結果更穩(wěn)健，所以下文對企業(yè)的所有制異質(zhì)性、所屬行業(yè)異質(zhì)性和四大審計異質(zhì)性的分析仍采用此模型，結果如表7所示。

(一)所有制異質(zhì)性

國有企業(yè)通常會積極響應國家政策，發(fā)揮帶頭作用，同時國有企業(yè)獲得政策信息的時點一般會早于非國有制企業(yè)。因此，下文驗證新能源示范城市試點中的國有企業(yè)是否在綠色創(chuàng)新中起到了引領作用。具體做法為：首先按照企業(yè)所有制類型篩選數(shù)據(jù)，然后對各個類型的數(shù)據(jù)進行回歸。

表7中第1-2欄分別為非國有企業(yè)、國有企業(yè)在PSMs-DID中的回歸結果。對比這兩個結果，發(fā)現(xiàn)國有企業(yè)did系數(shù)為正向顯著，而非國有企業(yè)為負向不顯著。相較于表2中的did系數(shù)0.040，表7中國有企業(yè)的did系數(shù)0.058有明顯提升，表明在新能源示范城市政策影響下國有企業(yè)在綠色創(chuàng)新中確實起到了帶頭作用。

(二)所屬行業(yè)異質(zhì)性

若某地級市準備申請新能源示范城市或已被確定為新能源示范城市，則這座城市中屬于高能耗行業(yè)的企業(yè)受影響最大，由此可以驗證高能耗行業(yè)的企業(yè)綠色創(chuàng)新的政策效應，并進一步驗證政策的實施是否起到了有的放矢的作用。依據(jù)高能耗行業(yè)的分類，將采礦業(yè)，制造業(yè)，電力、熱力、燃氣及水產(chǎn)和供應業(yè)，建筑業(yè)，交通運輸、倉儲和郵政業(yè)五類行業(yè)歸為高能耗行業(yè)，其他行業(yè)為低能耗行業(yè)。表7中第3-4欄為低能耗與高能耗行業(yè)的企業(yè)綠色創(chuàng)新回歸結果。對比這兩個結果發(fā)現(xiàn)，高能耗行業(yè)相比低能耗行業(yè)，企業(yè)綠色創(chuàng)新的貢獻度更大。低能耗行業(yè)的did系數(shù)為10%顯著水平負向顯著，而高能耗行業(yè)為1%顯著水平正向顯著，可見新能源示范城市政策對高能耗行業(yè)的企業(yè)確實起到了約束作用，而且新能源示范城市政策的實施具有高度針對性。

表7 所有制、所屬行業(yè)、四大審計異質(zhì)性的回歸結果

(三)四大審計異質(zhì)性

財務狀況好的公司會請四大審計公司進行審計，而四大審計公司一般也會審計財務好的公司，這是互為因果的關系。財務狀況好的公司是否會主動承擔企業(yè)責任、積極響應政府號召呢？自2006年中國成為全球第一大碳排放國后，企業(yè)的競爭必定是綠色、環(huán)保、低碳、生態(tài)的競爭已成為共識，誰在這些方面的創(chuàng)新走在前面，誰就有可能在未來立于主導地位。財務狀況好的公司有充足的資金擔負綠色專利研究的費用，以備公司的長遠發(fā)展；而財務狀況不佳的公司不會有閑置資金去搞綠色研發(fā)，更不會邀請四大審計公司進行審計。因此，本文將四大審計作為企業(yè)是否有充足資金的指標，而只有在確定企業(yè)有閑置資金時，政策的實施是否能有效促進企業(yè)綠色創(chuàng)新才能得到進一步驗證。

非四大審計與四大審計的回歸結果如表7中第5-6欄所示。二者的樣本量之比為7 579∶810，然而，樣本量少的四大審計企業(yè)對綠色創(chuàng)新具有顯著性，非四大審計企業(yè)不具有顯著性。四大審計企業(yè)的did系數(shù)為0.161，相較于表2中的回歸結果0.040，足足增加了3倍。由此進一步得出結論：企業(yè)的資金是否充足是影響政策對企業(yè)綠色創(chuàng)新產(chǎn)生效應大小的關鍵因素。

六、結論與對策建議

本文基于2003-2019年我國A股上市公司的面板數(shù)據(jù)，采用傾向得分匹配模型(PSM)和雙重差分法(DID)相結合的方法，檢驗評估了新能源示范城市政策對企業(yè)綠色創(chuàng)新的影響效應，研究結論如下。第一，新能源示范城市政策誘發(fā)了試點城市內(nèi)企業(yè)綠色創(chuàng)新行為，企業(yè)綠色創(chuàng)新絕對水平和相對水平顯著提升。第二，從綠色專利類別來看，相對于綠色發(fā)明專利，新能源示范城市政策更顯著影響企業(yè)對綠色實用新型專利的研發(fā)。第三，新能源示范城市政策在其實施的第4年(2017年)對企業(yè)綠色創(chuàng)新的政策效應幾乎減弱為零，說明政策效應的持續(xù)期大致為3年。第四，在企業(yè)綠色專利授權量的0.6分位數(shù)至0.9分位數(shù)的水平下，政策效應did系數(shù)隨著分位數(shù)的增大而顯著增大。在企業(yè)調(diào)節(jié)變量的政策效應中，企業(yè)發(fā)展屬性(企業(yè)規(guī)模)的did系數(shù)較為穩(wěn)定，而企業(yè)盈利屬性(現(xiàn)金流比率)的did系數(shù)值最大，且隨分位數(shù)的增大而升高，是企業(yè)綠色創(chuàng)新的重要影響因素。第五，國有企業(yè)比非國有企業(yè)研發(fā)創(chuàng)新更加積極，新能源示范城市政策對前者影響更顯著；高能耗行業(yè)的企業(yè)受新能源示范城市政策的影響比非高能耗行業(yè)的企業(yè)大，說明新能源示范城市政策非常具有針對性；四大審計企業(yè)的政策效應did系數(shù)比非四大審計企業(yè)大，說明企業(yè)的資金充足是影響政策對企業(yè)綠色創(chuàng)新產(chǎn)生效應的關鍵因素。

基于上述結論，本文提出如下對策建議：首先，新能源示范城市政策對企業(yè)綠色創(chuàng)新活動具有正向促進作用，表明推行新能源示范城市政策符合我國可再生能源的發(fā)展策略，因此，政府可以考慮進一步擴大新能源示范城市政策的實施范圍，貫徹“綠水青山就是金山銀山”的新時代綠色發(fā)展理念；其次，新能源示范城市應鼓勵企業(yè)加大對綠色發(fā)明專利的研發(fā)投入，長期的環(huán)境治理必須依靠綠色發(fā)明專利與綠色實用新型專利的支持，而研發(fā)投入是綠色創(chuàng)新的基本動力，對此政府也應給予適當補貼，以便更加有效推動企業(yè)的自主創(chuàng)新；再次，新能源示范城市的企業(yè)綠色創(chuàng)新熱度維持在3年左右，故應以3年為一個周期做好新能源示范城市政策的宣傳工作，使政策效應更持久；最后，對于不同類型的企業(yè)，地方政府應因“類”制宜地制定環(huán)境規(guī)制政策。