地方政府竞争下环境规制对区域碳生产率的非线性影响
——基于门槛特征与空间溢出视角
李珊珊，罗良文
(中南财经政法大学 经济学院，武汉 430074)
内容提要 ：立足于中介效应视角，本文从技术创新、产业结构、环境分权三个方面探讨环境规制对区域碳生产率影响的理论机制，并以此为基础，进一步运用2008-2015年中国30个省市区面板数据考察环境规制对区域碳生产率的门槛效应及空间溢出效应。实证结果显示：地方政府竞争背景下，环境规制对区域碳生产率的影响会随技术创新、产业结构以及环境分权变化而发生根本性转变。换言之，在不同技术创新、产业结构、环境分权水平条件下，环境规制对区域碳生产率的净效应分别呈现V型、倒V型、倒V型门槛特征；邻近地区间环境规制与区域碳生产率均存在局部聚集的空间特征，表现为明显的正向溢出效应。为实现2030年减排承诺的目标，除选择适宜的环境规制方式、力度以外，还需要从技术创新、产业结构、环境分权及区域协调等多个方面，优化环境规制的区域碳生产率效应的实施效果。
关键词 ：环境规制；碳生产率；门槛特征；空间面板杜宾模型
2015-2016年期间，中国碳生产率呈现新的变化趋势，即自1998年以来首次能源碳排放量与经济增长的脱钩。但根据《全球碳预算报告》预测，2017年能源消费碳排放量将逆势反弹3.5%，这一反弹与煤炭、石油、天然气的消费增长密切相关，反映了高碳“锁定效应”背景下经济增长对化石能源的依赖。为实现中国在《巴黎协定》中承诺2030年单位GDP碳强度下降60%-65%的目标，考虑到环境的双重负外部性，单纯依靠市场力量推动低碳经济转型远远不够，需要来自环境规制的激励与约束。本文通过考察环境规制对区域碳生产率影响的中介效应，试图改善现有环境规制的减排效应，为寻求摆脱高碳排放的路径依赖，更好地服务于低碳经济转型提供参考。
一 、文献综述
国内外学者关于环境规制影响的研究成果颇丰，主要集中于环境规制对经济增长的影响，或者主要考虑环境规制对生态环境质量的影响，然而，在经济增长和生态环境改善的两难困境中，单纯从某一方面探讨环境规制的影响无助于困境的根本解决。
环境规制与低碳经济关系的国内外文献主要包括“绿色悖论效应”、“倒逼效应”以及动态关系等三种观点。其中，“绿色悖论效应”的观点最早由Sinn(2008)提出，指明环境规制促进碳排放水平增长的可能性路径，包括环境规制设计失效、高碳能源供给的短期反弹以及时间滞后等路径。随后一些文献立足于高碳能源供给反弹机制的路径支持“绿色悖论效应”观点(Gerlagh, 2010; Michielsen, 2014)，并按碳排放增长的时效区分为“强绿色悖论效应”与“弱绿色悖论效应”两种(Van der Werf and Di Maria，2012)；支持“倒逼效应”的观点认为高额碳税不仅能降低长期碳排放总量，预期的改变还能减少短期碳排放(Hoel，2010)，表现为当清洁能源成本相对于化石能源更低时，环境规制会倒逼高能耗、高排放企业运用成本更低廉的清洁能源(Van der Ploeg and Withagen，2012)。早期国内的文献认为环境规制与低碳经济之间无明显关系(何小钢等，2011；张先锋等，2014)，近年来文献大多支持两者关系动态变化的观点，对环境规制的碳减排路径及效应关注较多，认为环境规制与低碳经济之间呈现倒U型变化趋势，而这一拐点的跨越与环境规制水平(张华，2014)、产业结构(徐盈之等，2015)以及技术创新(李巍等，2017)密切相关。
进一步梳理国内外学者关于碳生产率的文献发现，研究主要集中在碳生产率测度方法的构建及其影响因素方面。从研究方法的角度对碳生产率影响因素的文献进行梳理，将相关研究分为两类。一是以指数分解技术为基础的碳生产率因素分解研究。相关文献研究因指数分解技术的差异而存在分解因素的差别，如何建坤、苏明山(2009)将碳生产率增长分解为经济结构、能源技术效率和能源结构变动因素，Meng and Niu(2012)则运用LMDI因素分解技术将碳生产率分解为技术创新和产业结构调整因素。二是探讨碳生产率与特定因素相关关系的研究。包含碳生产率与某些特定因素相关关系的研究，如分别针对碳生产率与碳排放总量(谌伟等，2010)、人均GDP(张成等，2013)、能源效率(吴晓华、李磊，2014)之间动态关系的研究，而结合环境规制因素进行考察的文献结论存在较大差异，分别得出环境规制对区域碳生产率的影响呈现先抑制后促进的效应(刘传江等，2015)，环境规制会促进碳生产率增长(李小平等，2016)。
然而，对于环境规制促进碳生产率变动的诱导机制以及研究结论差异的内在原因，鲜有文献涉及。此外，若在环境规制对区域碳生产率的影响研究中忽略空间特征，则会忽略各区域碳生产率的空间集聚、溢出特征等重要影响因素，解决这一问题还需构建空间计量模型。本文有如下研究贡献：第一，从理论上探讨技术创新、产业结构、环境分权等中介变量的诱导机制，并在两者相关关系的基础上，运用非线性门槛面板模型考察环境规制对碳生产率影响的门槛特征；第二，结合地方政府竞争背景，在空间面板杜宾模型中纳入环境规制、区域碳生产率的空间交互项，进一步考察门槛特征下环境规制对区域碳生产率的空间溢出效应。
二 、环境规制影响区域碳生产率的中介效应机制
(一)环境规制对区域碳生产率的影响机制
式(1)、(2)中，t 为时间，i 为生产者数量，σ 为两种中间产品的替代弹性，σ >1，α 为资本产出份额，0<α <1，Y 为最终产品的产出，将劳动力总投入量标准化为1，则L ct +L pt ≤1，A 为技术进步，K 为中间产品物质资本投入。
1.假设
(1)假设国家有三个生产部门，分别为中间产品物资资本供应部门、中间产品生产部门、最终产品生产部门，劳动、资本要素可以自由流动，其中，最终产品市场、中间产品市场是完全竞争市场，而中间产品物质资本供应商提供与技术创新融合的机器设备，具有市场垄断势力。
古龙说：“爱笑的女孩，运气都不会太差。”想想的确如此。好脾气像十月，秋高气爽，让人心旷神怡；坏脾气则像六月，说变就变，让人防不胜防，避之不及。
(2)假设整个生产过程分为两阶段，第一阶段生产中间产品，其生产非低碳技术设备，也可引进低碳技术设备，其中，引进低碳技术设备进行生产的中间产品，用X c 表示，其生产过程不产生碳排放；另一种为运用非低碳技术设备生产的中间产品，用X p 表示，其生产过程产生与产量同比例的碳排放，并将中间产品生产函数形式设定为Cobb-Douglas函数以及资本增加型技术进步。第二阶段运用中间产品生产最终产品，生产函数形式设定类似。
3.后备猪的初次配种，一般在初情期后的第3个情期开始配种,此阶段后备母猪已达到体成熟，能满足怀孕的需要，采取一个情期多次配种，提高产仔数。
最终产品生产函数：
(1)
中间产品生产函数：
(2)
本部分在Acemoglu (2012)分析框架的基础上拓展，将环境规制内生化，构建污染品生产者通过低碳、非低碳技术产品相对生产规模改变以及区域碳生产率变化的理论模型。主要思路为，假定中间产品区分为低碳技术中间产品和非低碳技术中间产品两种，其中低碳技术中间产品生产促进产量的增加而不产生碳排放，而非低碳技术中间产品生产在提高产量的同时同比例增加碳排放，则环境规制水平的提高会改变非低碳技术中间产品与低碳技术中间产品的相对利润，进而通过中间产品相对生产规模变化影响区域碳生产率变动。
(3)根据Acemoglu(2012)的模型设立中间产品的技术创新函数，生产者最初选择非低碳技术创新模式，而未考虑环境规制的压力实施低碳技术自主研发创新。
中间产品技术创新函数：
而低碳技术中间产品生产者利润π ct 的最大化条件为：
(3)
初始状态下两种技术创新模式关系为：
(4)
式(3)、(4)中，γ 为技术存量增加率，η 为技术创新成功的概率，S 为技术创新活动参与的技术人员数，r 为中间产品的物质资本的投入成本，类似地，将技术人员总数标准为1，则S ct +S pt ≤1。
(4)假设中间产品X p 的生产会产生环境负外部性，在产权界定明晰的前提下，中间产品X p 的生产需要支付碳排放的机会成本τ ，对应现实中的环境规制成本。
2.模型构建
研究组患者接受少腹逐瘀汤治疗，组方包括干姜3 g、五灵脂6 g、延胡索10 g、蒲黄10 g、小茴香6 g、当归15 g、赤芍10 g、没药10 g、肉桂6 g以及川芎10 g，上述药物用水煎服，指导患者少腹逐瘀汤月经来潮服用。一天2次，一次200 ml，连续口服5天，3个月经周期为1个疗程[4-5]。
lnCP it =α 0+α 1lnER i,t -1·I (M ≤β 1)+α 2lnER i,t -1·I (β 1<M ≤β 2)+…+α n lnER i,t -1·I (β n-1 <M ≤β n )+α n+1 lnER i,t -1·I (M >β n )+X it ·γ +ε it
该公司为株式会社捷太格特在中国的独资公司，机床产品业务包含外圆磨床、无心磨床、卧式加工中心、模具加工中心、专机、超精密加工机、程序控制器（PLC）以及数控产品等。公司的机床制造立足于发动机行业高运转率、高生产效率、低故障和低运行成本的原则，追求设备的长期稳定性、高精度和高刚性。
(5)
1.低碳技术创新
(6)
(7)
求解非低碳技术中间产品物质资本供应商的利润最大化条件为：
maxπ mt =m t K pt -α 2r p K pt
(8)
将式(6)代入式(8)，并对π mt 求关于非低碳技术中间产品价格m t 的最优化条件：
m t =αr p
李白、鲍照都创作了为封建女性的悲剧命运代言的诗歌，如《白头吟》和《王昭君》。鲍照借女子失宠写正直之士不容于世，如《选诗补注》所说“此殆明远为人所间，见弃于君，故借是题以喻所怀。”失宠女子的命运中依稀有作者自己的身影，表面上是代言，实为自己抒愤。
(9)
将上述式(7)、式(8)、式(10)代入非低碳技术设备供应商利润π mt ：
2017年7月7 —11日，内蒙古自治区兴安盟出现强降水天气过程，期间伴有短时强降水、冰雹、雷雨大风等强对流天气。7日08：00至11日08：00，3个站累积降雨量超过50.0 mm，分别为中旗吐列毛都镇罕查干55.9 mm、阿尔山五岔沟镇南沟55.8 mm、中旗吐列毛都54.0 mm。累积降雨量为25.0～50.0 mm的站点有55个，主要分布在中西部地区。主要降水时段在7日白天至8日白天。7日11：00左右，降水从南部地区开始，随着时间推移雨带逐渐向北推进。7日白天降水集中在中南部地区，夜间降水集中在中北部地区，8日白天之后则以分散性的雷阵雨天气为主，9日以后降水强度和范围明显减弱。
(10)
A jt =(1+γη j S jt )A j,t -1,j =c ,p
准确搜集气象信息是实现气象预报准确率提高的基础。因此，必须加强对地面综合气象观测的建设。在建设过程中，应该将气象预测作为建设目的和中心任务，利用最新技术协调各方力量，实现综合气象观测能力的提高。在对相关体系进行建设的过程中，气象工作人员必须熟悉气象知识和当地气象条件。气象信息搜集应严格依据相关规范，对于搜集到的气象信息，要综合汇总分析，运用专业软件进行计算，以提高结果的准确性。
(11)
与π mt 的推导过程类似，得出低碳技术设备供应商利润π nt ：
(12)
接下来，将低碳技术和非低碳技术供应商利润π nt /π mt 进行对比，得出环境规制水平、低碳与非低碳技术物质资本相对成本、中间产品相对价格、相对技术创新、劳动力相对投入等因素的影响，如下：
π nt /π mt =[(1-τ t )(1-q t )]-1/(1-α) (r c /r p )-α/(1-α )(p ct /p pt )1/(1-α) (A ct /A pt )α/(1-α )(L ct /L pt )
(13)
由式(13)可知，相对技术创新、低碳与非低碳技术中间产品投入比、环境规制水平是影响低碳技术和非低碳技术供应商利润、相对生产规模及区域碳生产率的主要因素，因此，本文选择低碳技术创新、产业结构、环境分权作为中介变量，并分别考察其中介效应机制。
在化肥施用上，应以氮、磷为主，其中氮肥占总施肥量的60%、磷肥占80%、钾肥占20%。具体来说，亩施尿素31.4 kg、过磷酸钙96 kg、硫酸钾8.7 kg，或二铵25 kg、尿素21.6 kg，配施生物菌肥 50～60 kg。
(二)中介效应
分别求对π pt 求关于K pt 、L pt 的最优化得到：
进行磨料质量百分浓度对线材表面去除效率的影响实验时，将实验中高压泵的压力设定为30MPa,流量为80L/min；所用的高压管道直径为16mm,宝石喷嘴直径为0.45mm,磨料聚焦管直径为1.5mm、长度为76mm,磨料流量为20mL/s；喷嘴对线材的入射角度为水平向下倾斜45°,喷嘴距线材 30mm；采用直径为6mm的45号钢线材为加工材料,加工过程中线材匀速移动,速度为20mm/s；采用80目的石榴石砂[10],其中磨料的质量百分浓度分别为40%、50%、60%、70%。
当低碳技术存量较低时，非低碳技术物质资本供应商利润大于低碳技术物质资本供应商利润，中间产品为非低碳技术中间产品；当低碳技术存量累积到某一阶段时，低碳技术物质资本供应商利润与非低碳技术物质资本供应商利润相等时，低碳技术中间产品和非低碳技术中间产品并存；当低碳技术存量累积到某一门槛值时，低碳技术中间产品生产规模相对扩大。由此可预期，较低的低碳技术存量难以促进低碳技术中间产品生产，此时环境规制约束下低碳技术创新投入对区域碳生产率的贡献不明显，主要表现为环境规制的“成本遵循效应”(Jaffe and Stavins, 1995)，只有当低碳技术突破某一门槛值，环境规制约束才能发挥低碳技术创新的正向促进作用，可能部分削弱或完全抵消“成本遵循效应”，表现为环境规制的“波特效应”(Porter and Van der Linde，1995)。
工程岩体在地质作用以及人为扰动作用下不仅表现为短时的弹塑性力学特性，而且在一定应力应变条件下呈现与时间效应相关的变形和破坏性质，这种性质即为岩体的流变特性[1]。一般情况下，巷道围岩处于复杂的地质环境中，尤其深部围岩在应力作用下，呈现出明显的流变性质[2-3]，但是经过国内外学者的大量研究，验证了地质环境中水因素对岩石变形影响较大，主要包括：可逆物理作用、不可逆化学作用和基本的力学作用等[4]。而各种影响因素的组合作用，使得岩石内部结构体系发生劣化，甚至是物质组成改变，导致岩石内部裂隙大量生成发育，最终贯通为宏观裂纹；同时，也会造成岩石力学参数降低，抵抗破坏能量削弱[5]。
2.产业结构
S it =S (KS it ,IS it )且T it =T (ER it ,PAT it ,HR it )
3.环境分权
环境分权本质上通过中央与地方之间环境相关财权、事权分配，促进对地方政府环境治理努力程度的激励与约束。环境集权下由中央统一供给环境公共物品，能有效约束地方政府环境治理的“搭便车”行为，并实现环境公共物品生产的规模经济(Gray and Shadbegian，2004)，当环境集权水平持续下降，较高的环境分权水平会适度弱化中央对地方政府环境治理努力程度的约束，但所给予的地方政府环境管理方式的自由，有利于地方政府通盘考虑本地区经济基础、环境质量，灵活选择相对高效的碳减排治理方案(Magnani，2000；Millimet，2003)，此时环境分权水平适度提高可能会强化地方政府环境治理努力程度的激励；过高的环境分权则意味着地方政府环境管理自主权与管控能力，容易形成地方环境管理部门与污染企业之间的合谋，导致为隐瞒真实的污染排放而篡改数据的行为(聂辉华等，2007；梁平汉等，2014)，对于短期危害小且空间外溢性强的碳排放而言，更容易导致对地方政府环境治理努力程度的软化。因此，环境规制的区域碳生产率效应与环境分权水平密切相关。
三 、面板模型 、变量与数据来源
(一)门槛面板模型
门槛效应是指环境规制对区域碳生产率的促进并不是自动进行的，其正面促进作用需要突破一个或几个关键因素，否则环境规制效应难以得到发挥。根据第三部分推导，本文选取的关键因素为低碳技术存量、产业结构、环境分权，作为内生的门槛变量，运用Hansen(1999)提出的门槛面板模型考察环境规制驱动区域碳生产率提升的过程中是否存在门槛效应。考虑到环境治理投资从投入到效果显现存在一定的时滞，本文选用环境规制变量的滞后一期进行考察，构建模型如下：
在上述假设基础上，求解非低碳技术与低碳技术中间产品生产商、物质资本供应商利润最大化条件，非低碳技术产品生产者减排通过引进低碳技术设备进行生产。若对非低碳技术产品以从价税τ 方式征收环境税，同时对低碳技术产品的低碳技术设备成本以补贴率q进行补贴。其中，非低碳技术中间产品生产者利润π pt 的最大化条如下：
(14)
式(14)中，β 1、β 2…β n 为变量M 的门槛值，变量M 包括技术创新、产业结构以及环境分权，I (·)为示性函数，ER 为环境规制，CP 为碳生产率，X 为控制变量的集合，ε 为随机干扰项，其余为相应变量的参数。其中，对于式(14)中控制变量X 的选择，根据kaya恒等式(15)及结构、技术创新等式(16)：
C it =Y it ×S it ×T it
(15)
适度的环境规制成本能倒逼污染密集型企业的低碳化转型(Porter and Van der Linde，1995)，或促使企业迁至环境规制较为宽松的地区。对污染密集型企业依赖程度较低的地区，内化的环境规制成本较低，无法扭转低碳技术与非低碳技术中间产品的相对利润及生产规模，进而对区域碳生产率变动的影响不显著；对污染密集型企业依赖程度较高的地区，内化的环境规制成本较高，会明显改变低碳技术与非低碳技术中间产品相对利润及生产规模，有利于强化环境规制对区域碳生产率的正面效应；对于少数过度依赖污染密集型产业的地区，产业低碳化转型也可能导致该地区经济增速大幅放缓，为满足地方经济、就业的发展需求，地方政府可能运用财政补贴、税收减免等方式部分抵消“环境遵循成本”压力，导致环境规制的区域碳生产率效应弱化。因此，环境规制的区域碳生产率效应会随产业结构调整而动态变化。
(16)
⟹CP it =Y it /C it =1/[S (KS it ,IS it )×T (ER it ,PAT it ,HR it )]
5) 可靠性技术研究。研究多余度结构的磁通切换电动机，提出切实可行的余度结构方案；研发模块化、集成化的控制器，优化软、硬件结构，提高系统抗干扰能力和可靠性。
(17)
式(15)-(17)中，Y 为产出，S 为结构变量，包括要素禀赋结构KS 、产业结构IS ，T 为技术水平，由环境治理投资ER 、技术创新PAT 、人力资本HR ，此外还与对外开放程度FDIS 与TRS 有关，上述变量为区域碳生产率影响模型控制变量X 的集合。
(二)空间面板杜宾模型
结合中国式分权下地方政府竞争的背景，各地区环境规制可能存在相互模仿或策略性行为，或各地区碳生产率可能表现出空间聚集特征，实证研究中若忽略核心变量的空间相关性，则会导致估计结果的严重偏误。为验证环境规制、碳生产率的空间相关性，本文考察2008-2015年期间环境规制、碳生产率的moran’ I指数及其显著性水平，结果显示，环境规制、碳生产率的moran’ I指数分别为0.315、0.422，均达到1%的显著性水平，这表明实证模型中需纳入两者的空间交互特征。空间面板杜宾模型能同时纳入因变量与自变量的空间交互项，因此本文选用空间面板杜宾模型进行参数估计，模型如下所示：
(18)
式(18)中，空间权重矩阵w ij 分别为地理邻接权重矩阵与地理距离权重矩阵，其中地理邻接权重矩阵相邻为1，不相邻为0，而地理距离权重矩阵为区域之间质心距离平方的倒数表示。
(三)变量与数据来源
区域碳生产率CP ，运用各地区GDP 与能源消费碳排放的比值表示，其中碳排放包括煤、汽油、柴油、天然气消费的直接碳排放，还包括电力消费的间接碳排放；环境规制，运用各地区环境治理投资额占GDP 的比重表示；此外，结构控制变量，要素禀赋结构KS 用各地区资本存量与就业比值表示，其中资本存量借鉴张军等(2004)的方法按永续盘存法进行估算，产业结构IS 用各地区第三产业产值占GDP 的比重表示；技术控制变量，技术创新PAT 用各地区每万人拥有的专利授权项数表示，人力资本HR 用各地区劳动力人均受教育年限表示；开放程度相关变量，对外开放程度FDIS 用各地区FDI 占GDP 比值表示，TRS 用各地区进出口贸易总额占GDP 比值表示；环境分权ED ，借鉴祁毓等(2014)的思路，运用环境保护系统人员的地区分布刻画环境分权水平，具体公式为
等式右边第二项为经济规模影响因子，引入该因子能剔除经济规模变化的偏差。
考虑到地方官员政绩考核体系变化对环境规制实施的影响，如2007年底我国首次将单位GDP 能耗降低的约束性目标纳入地方官员考核体系，占考核总分的40%左右(齐晔，2013)，为推动碳强度目标下降的承诺，于2014年进一步将单位GDP 碳排放指标纳入地方官员考核体系。因此，本文以2008-2015年中国30个省市区面板数据为样本，能更为准确地反映门槛变量对环境规制效应的影响方向与影响程度。各变量数据的描述性统计如表1所示。
表1 变量数据描述性统计
注：由Stata 14.0软件整理得到，保留三位小数。
表2 门槛效应估计与检验结果
注：各统计值为bootstrap法反复抽样1000次所得；***、**、*分别代表1%、5%、10%的显著性水平。
四 、模型估计与结论分析
(一)门槛效应检验
本文运用Stata 14.0 软件进行分析，首先根据bootstrap方法的F统计值和P值确定技术创新、产业结构、环境分权的门槛数量，然后运用门槛面板模型估计其门槛值，以及该门槛值下环境规制变量的显著性水平。表2中技术创新的单一门槛均达到1%的显著性水平，门槛值为1.749，而双重门槛、三重门槛效应在统计上不显著；产业结构的单一门槛达到1%的显著性水平，门槛值为3.723，而双重门槛、三重门槛效应在统计上不显著；环境分权的单一门槛、双重门槛均达到1%的显著性水平，门槛值分别为-0.176、0.482，而三重门槛效应在统计上不显著。这一结论与技术创新、产业结构、环境分权中介效应的理论分析一致。
以技术创新、产业结构以及环境分权的单一门槛值为界，根据门槛变量均值对样本进行分组：低技术创新水平区域[-0.194, 2.414]，高技术创新水平区域(2.414, 3.210]；低产业结构偏向区域[3.498, 3.592]，高产业结构偏向区域(3.592, 4.339]；低度环境分权水平区域[-0.630, -0.196]，中度环境分权水平区域(-0.196, 0.482]，高环境分权水平区域(0.482, 0.744]。不同的样本分组具体对应的省市区如表3所示。
表3 按门槛变量门槛值的样本分组结果
(二)空间面板杜宾模型估计结果
根据技术创新、产业结构以及环境分权水平，对样本进行分组，分别考察不同区域环境规制对碳生产率影响的异质性，其中按环境分权门槛分类的三组样本，高度环境分权水平组仅山西省，对此，本文将高度与中度环境分权水平组合并为中高度环境分权水平组，如表4、表5所示。对比基于地理邻接权重与地理距离权重的估计结果，不同技术创新、产业结构以及环境分权水平的区域，核心变量的估计结果基本一致，仅在估计系数大小及显著性方面存在差异，保证了估计结果的稳健性。
对于表4、表5的模型AI、AII估计结果而言，技术创新水平的差异会明显影响环境规制的碳生产率效应，呈现“V型”的单一门槛特征，具体来看，在技术创新水平较低的区域，即每万人拥有的专利授权数小于e 2.414(约11项)时，环境规制变量的估计系数分别为-0.019、-0.024，对应5%、15%的显著性水平，表明环境规制强化会抑制区域碳生产率增长，此时环境规制成本无法刺激企业创新或引进节能减排技术，反而因挤占企业生产成本而导致技术创新投入以及能源利用效率的降低；在技术创新水平较高的区域，即每万人拥有的专利授权数大于e 2.414(约11项)时，环境规制变量的估计系数均为正，但不显著，说明环境规制对区域碳生产率的净效应不明显，这是由于随着技术创新能力的积累，环境规制的“波特效应”开始凸显，逐步抵消环境规制“成本遵循效应”的负面影响，前文所述技术创新的中介效应机制得以印证。从分组对应的区域来看，高技术创新水平区域仅占样本总量的1/5，包括北京、天津、上海、江苏、浙江、广东等6个地区，技术创新的整体水平较低，限制了技术创新对环境规制“波特效应”的正面激励，存在较大的提升潜力。
总之，基于可靠性思维推出的状态监控借助于各种传感器、信号采集、信号处理和信息分析等技术和手段，对来自于汽车、发动机系统、设备、部件、组件与相关的外部数据或信号进行检测和观察，并通过人机对话判断汽车、发动机技术状态，制定决策，确定维修内容，这是一门综合性与工程技术实践与管理。
对于表4、表5的模型BI、BII的估计结果而言，产业结构对环境规制碳生产率的影响效应是显著的，呈现为“倒V型”的单一门槛特征，具体来看，在产业结构偏向较低的区域，即第三产业产值比重小于36.31%时，环境规制变量的估计系数分别为0.028、0.024，均达到1%的显著性水平，说明环境规制强化有利于区域碳生产率增长，原因可能在于第三产业发展滞后，反映出该区域对高能耗产业的依赖程度较高，环境规制内化的高成本会促使高能耗产业的低碳转型，或向环境规制较为宽松地区的迁移，进而改善该区域碳生产率；在产业结构偏向较高的区域，即第三产业产值比重大于36.31%时，环境规制变量的估计系数均为负值，但不显著，此时该地区对高能耗产业的依赖程度较低，环境规制成本内化对产业结构低碳化转型的动力很有限。从分组对应的区域来看，低产业结构偏向区域占样本总量的1/5，包括河北、安徽、江西、河南、陕西、青海等6个地区，产业结构调整的中介效应存在一定程度的贡献。
表4 空间面板杜宾模型估计结果(地理邻接权重矩阵)
注：***、**、*、o分别代表1%、5%、10%及15%的显著性水平，括号内为对应的p值。
对于表4、表5的模型CI、CII的估计结果而言，随着环境分权水平的提高，环境规制对碳生产率的净效应发生了根本性变化，呈现出“倒V型”的门槛特征，具体来看，在低度环境分权水平的区域，即环境分权水平低于e -0.196(0.822)时，环境规制变量的估计系数分别达到5%、10%的显著性水平，表明环境规制强化有利于提升碳生产率，这是由于适度的环境分权给予地方政府立足于本地实际情况灵活选择环境规制方式的权利，形成地方政府竞争背景下区域间环境规制的“逐顶竞争”，该区域仅占样本总量的2/5，存在进一步拓展的空间；在中高度环境分权水平区域，即环境分权水平高于e -0.196(0.822)时，环境规制变量的估计系数均在1%的水平上显著为负，此时环境规制的碳生产效应的根本转变，是由于过高的环境分权会明显削弱中央对地方环境治理的约束，进而导致地方政府与高能耗企业合谋下环境规制的软化。因此，适度的环境规制水平有助于环境规制碳生产率的正面效应。
从表4、表5的控制变量来看，技术创新变量的估计系数均为正，表明本土技术创新能明显改善能源利用效率，进而促进碳生产率提升；人力资本变量的估计系数除模型BI外，其余模型均为正，说明对成熟技术的吸收、运用或扩散的能力也是制约能源利用效率改进的关键；产业结构变量的估计系数仅在表4的模型AII、表5的模型BII达到10%的显著性水平，其余模型均不显著，这一估计结果与产业结构调整的方向有关，2008年后产业结构调整方向转向淘汰高能耗、高排放的落后生产能力，因此环境规制约束下产业结构调整是促进碳生产率提升的驱动因素，正如模型BI中环境规制变量的估计结果，即产业结构调整的影响主要通过其中介效应机制得以体现；表4、表5中要素禀赋结构、对外开放变量的估计结果差异较大，结论不稳健。关于空间滞后项，除表4的AII模型、表5的BI模型外，其余模型碳生产率的空间滞后项均显著为正，表明地理邻近地区的碳生产率具有较强的局部聚集空间特征，此外，环境规制空间滞后项估计系数在表3、表4的AII、CII模型中显著为正，这意味着在高技术创新、中高环境分权水平区域，邻近地区环境规制强化表现出明显的正向溢出效应，即邻近地区碳生产率与环境规制均存在空间极化现象。
表5 空间面板杜宾模型估计结果(地理距离权重矩阵)
注：同表4。
五 、主要结论与政策启示
(一)主要结论
本文运用2008-2015年面板数据，以及非线性门槛面板模型、空间面板杜宾模型，结合地方政府竞争背景，从技术创新、产业结构及环境分权三个方面考察环境规制对区域碳生产率的中介效应，研究显示，技术创新、产业结构及环境分权的变化，会促使环境规制的碳生产率净效应发生根本转变。具体包括：第一，随着技术创新水平的提升，环境规制对碳生产率的影响呈现出“V型”的单一门槛特征，即每万人拥有的专利授权数大于11项时，环境规制的“波特效应”能抵消“成本遵循效应”的负面影响，由于该区域仅占样本总量的1/5，存在进一步提升的潜力；第二，随着产业结构的绿色化调整，环境规制对碳生产率的影响表现为“倒V型”的单一门槛特征，即产业结构偏向度低于36.31%时，环境规制内化的高成本倒逼企业从事生产末端环境治理活动，或实施生产过程的清洁技术改造，该区域占样本总量的1/5，可适度强化该区域环境规制力度；第三，随着环境分权水平的提高，环境规制的碳生产率效应呈现为“倒V型”的门槛特征，即环境分权水平低于0.822，更有助于促进环境规制的实施力度及其正面激励效应；第四，地理邻近地区环境规制、碳生产率存在局部聚集的空间特征，表现出明显的正向溢出效应。
(二)政策启示
地方政府竞争背景下环境规制的区域碳生产率净效应存在不确定性，我国应结合技术创新、产业结构、环境分权以及环境规制区域协调等多方面力促环境规制对区域碳生产率的正面效应，具体如下：
其实，罪魁祸首是沉迷，不是网游。沉迷便是兴趣过于浓厚，超出了自控力。未成年人自控力不强，容易陷于网游不能自拔，乃至走火入魔酿成悲剧。即使是成年人，如果陷于沉迷，后果同样不堪设想。
1.地方政府构建、落实创新孵化平台，拓宽风险资本来源以改善单一的投资结构，为企业节能减排技术创新提供充裕的风险资本支持，避免形成政府财力主导下孵化平台的行政化、非理性的现象。除技术创新支持外，还需要结合产学研平台联合培养相关专业的技术服务型人才，推动已有节能减排技术的扩散，促进技术创新到企业产品创新的转化，以更好地发挥技术创新或以技术创新为中介的环境规制碳生产率效应。
2.产业结构绿色化调整的中介效应影响有限，由于大部分地区已完成高能耗、高排放等产能过剩企业的转型，以服务业主导的产业结构限制了产业结构调整的空间，需要特别关注少数地区如河北、安徽、江西、河南、陕西、青海等地产业结构绿色化调整。
3.适度压缩地方政府环保事权、财权的自由裁量空间，规避地方政府与企业合谋导致的环境规制非完全执行，引导区域之间环境规制的逐顶竞争行为，或联合各地区共同减排，统一供给区域内环境公共物品的规模经济，以规模化生产降低环境公共物品的供给成本。
参考文献 ：
[1] 谌伟, 诸大建, 白竹岚.上海市工业碳排放总量与碳生产率关系[J].中国人口·资源与环境, 2010,20(9): 24-29.
[2] 何建坤,苏明山.应对全球气候变化下的碳生产率分析[J].中国软科学,2009(10): 42-47.
[3] 何小钢,张耀辉.行业特征,环境规制与工业 CO2 排放——基于中国工业 36 个行业的实证考察[J].经济管理,2011(11): 17-25.
[4] 李巍,郗永勤.创新驱动低碳发展了吗?——基础异质和环境规制双重视角下的实证研究[J].科学学与科学技术管理,2017,38(5): 14-26.
[5] 李小平,王树柏,郝路露.环境规制,创新驱动与中国省际碳生产率变动[J].中国地质大学学报: 社会科学版,2016,16(1): 44-54.
[6] 梁平汉,高楠.人事变更,法制环境和地方环境污染[J].管理世界,2014(6): 65-78.
[7] 刘传江,胡威,吴晗晗.环境规制,经济增长与地区碳生产率——基于中国省级数据的实证考察[J].财经问题研究,2015(10): 31-37.
[8] 聂辉华,李金波.政企合谋与经济发展[J].经济学(季刊),2006,6(1): 75-90.
[9] 齐晔.中国低碳发展报告(2013)[M].社会科学文献出版社,2013.
[10] 祁毓,卢洪友,徐彦坤.中国环境分权体制改革研究: 制度变迁,数量测算与效应评估[J].中国工业经济,2014(1): 31-43.
[11] 吴晓华,李磊.中国碳生产率与能源效率省际差异及提升潜力[J].经济地理,2014,34(5): 105-108.
[12] 徐盈之,杨英超,郭进.环境规制对碳减排的作用路径及效应——基于中国省级数据的实证分析[J].科学学与科学技术管理,2015(10): 135-146.
[13] 张成,蔡万焕,于同申.区域经济增长与碳生产率——基于收敛及脱钩指数的分析[J].中国工业经济,2013(5): 18-30.
[14] 张华,魏晓平.绿色悖论抑或倒逼减排——环境规制对碳排放影响的双重效应[J].中国人口资源与环境,2014,24(9): 21-29.
[15] 张军,吴桂英,张吉鹏.中国省际物质资本存量估算: 1952-2000[J].经济研究,2004(10): 35-44.
[16] 张先锋,韩雪,吴椒军.环境规制与碳排放:“倒逼效应” 还是 “倒退效应”——基于 2000～ 2010 年中国省际面板数据分析[J].软科学,2014,28(7): 136-139.
[17] Acemoglu D,Aghion P,Bursztyn L,et al.The environment and directed technical change[J].American economic review,2012,102(1): 131-66.
[18] Gerlagh R.Too much oil[J].CESifo Economic Studies,2011,57(1): 79-102.
[19] Gray W B,Shadbegian R J.‘Optimal’pollution abatement—whose benefits matter,and how much?[J].Journal of Environmental Economics and management,2004,47(3): 510-534.
[20] Hansen B E.Threshold effects in non-dynamic panels: Estimation,testing,and inference[J].Journal of econometrics,1999,93(2): 345-368.
[21] Hoel M.Is there a green paradox?[J].2010.
[22] Jaffe A B,Stavins R N.Dynamic incentives of environmental regulations: The effects of alternative policy instruments on technology diffusion[J].Journal of environmental economics and management,1995,29(3): S43-S63.
[23] Magnani E.The Environmental Kuznets Curve,environmental protection policy and income distribution[J].Ecological economics,2000,32(3): 431-443.
[24] Meng M,Niu D.Three-dimensional decomposition models for carbon productivity[J].Energy,2012,46(1): 179-187.
[25] Michielsen T O.Brown backstops versus the green paradox[J].Journal of Environmental Economics and Management,2014,68(1): 87-110.
[26] Millimet D L.Assessing the empirical impact of environmental federalism[J].Journal of Regional Science,2003,43(4): 711-733.
[27] Porter M E,Van der Linde C.Toward a new conception of the environment-competitiveness relationship[J].Journal of economic perspectives,1995,9(4): 97-118.
[28] Sinn H W.Public policies against global warming: a supply side approach[J].International Tax and Public Finance,2008,15(4): 360-394.
[29] Van der Ploeg F,Withagen C.Is there really a Green Paradox?[J].Journal of Environmental Economics and Management,2012,64(3): 342-363.
[30] Van der Werf E,Di Maria C.Imperfect environmental policy and polluting emissions: the green paradox and beyond[J].International Review of Environmental and Resource Economics,2012,6(2): 153-194.
The Nonlinear Effects of Environmental Regulation on Regional Carbon Productivity under Local Government Competition ：Based on the View of Threshold Characteristic and Spatial Spillover
LI Shan-shan, LUO Liang-wen
(School of Economics ,Zhongnan University of Economics and Law ,Wuhan 430074,China )
Abstract ：Based on the intermediate perspective, this paper studies the theoretical influencing mechanism from technological innovation, industrial structure and environmental decentralization.Then this paper estimates the threshold effects and spatial spillover effects of environmental regulation on regional carbon productivity based on the panel data of 30 provinces in China from 2008 to 2015.The empirical conclusion shows that, the regional carbon productivity effects of environmental regulation has a structural change with the changes of local technological innovation level, industrial structure deviation and environmental decentralization under the local governmental competition.In other words, under the different technical innovation levels,industrial structures and levels of environmental decentralization,the net effect of environmental regulation on regional carbon productivity shows a ′V′ type,inverted ′V′ type and inverted ′V′ type threshold characteristic respectively;both environmental regulation and regional carbon productivity in adjacent areas have spatial characteristics of local aggregation, which shows obvious positive spillover effect.In order to realize the target of emission reduction commitment, we should promote the practical effects of the environmental regulation on the regional carbon productivity from technological innovation, industrial structure adjustment, environmental decentralization and regional coordination, except for setting the suitable pattern and strength of environmental regulation.
Key words ：environmental regulation; carbon productivity; threshold characteristic; spatial Dubin model
中图分类号 ：F062.2
文献标识码： A
文章编号： 1001-148X(2019)01-0088-10
收稿日期： 2018-09-05
作者简介： 李珊珊(1980-)，女，湖北宜昌人，中南财经政法大学经济学院副教授，经济学博士，研究方向：环境经济学；罗良文(1965-)，男，河南商城人，中南财经政法大学经济学院教授，博士生导师，经济学博士，研究方向：国际经济与环境。
基金项目： 国家自然科学基金项目“环境规制对中国区域碳生产率影响的理论解释与非线性效应研究:基于地方政府竞争视角”，项目编号：71503272。
(责任编辑：周正)
标签：环境规制论文; 碳生产率论文; 门槛特征论文; 空间面板杜宾模型论文; 中南财经政法大学经济学院论文;