节能减排约束下的工业全要素生产率.doc

1、节能减排约束下的工业全要素生产率摘要：本文运用基于熵的 USBM 模型建立了包含能源、环境因素的曼氏指数，对我国 30 个省际 1992-2011 年的工业节能减排生产率及其分解进行测度，结果发现：传统的忽略节能减排因素的生产率测度会高估我国工业的全要素生产率；中部地区在节能减排约束下的工业生产率增长最快，且全部来自于技术进步；生产前沿面几乎都来自东部地区，它们是节能减排技术的引领者；全国省际工业效率水平与前沿面的差异性正在扩大，未考虑节能减排约束的省际工业效率测度低估了这一差异的严重性。 关键词：节能减排；环境规制；熵；全要素生产率 中图分类号：F0622 文献标识码：A 一、引言 我国能源

2、资源有限，煤炭和水力资源人均拥有量相当于世界平均水平的 50%，石油、天然气拥有量仅为世界水平的 1/15 左右，但能源等资源的消耗量却十分巨大。英国石油（BP） 世界能源统计年鉴显示，我国已经超过美国，成为了世界第一能源消费者。作为能源消费大户的工业，虽然贡献了一半的 GDP，却消耗了 70%的能源并制造了 70%的污染。因此，合理测度节能、减排约束下的工业全要素生产率，对于推动我国产业转型升级具有十分重要的意义。 国内外对于环境规制及中国工业绩效评价的文献众多，大都涉及了能源、环境对我国工业发展的约束。Guang-Ming Shi 等1运用 2000-2006 年的数据，发现中国部分地区的

3、工业仍然依赖能源的大量投入。陈诗一2运用改革开放以来的数据，发现能源是驱动工业增长的源泉之一。陈军等3运用 30 个省的数据，发现能源效率的改进有利于我国工业化水平的提高。李世祥和成金华4运用 1990-2006 年的数据，发现我国 13个工业省区能源效率差异较大。唐玲和杨正林5运用 1998-2007 年的数据，发现中国工业能源效率的总体水平较低，具有较大的节能和提升空间。邵军和管驰名6运用 1999-2005 年的工业数据，发现我国工业能源效率在 50%-60%之间，效率不高。王少平和杨继生7运用 1985-2002 年的时序数据，发现中国工业大多数行业对能源具有高消费和强依赖的特点。可见

4、，能源作为工业的重要投入，在工业 TFP 测算中必须加以重视。与此同时，还有一些学者将研究重点放在了工业的产出上，即纳入了影响环境的污染排放因素。涂正革8的研究发现，区域间环境、工业协调性极不平衡。吴军9的研究发现各地区 TFP 增长均主要来源于技术进步，且东部地区是推动环境技术创新的主要地区。杨俊和邵汉华10运用1998-2009 年的地区工业数据，发现忽略环境因素会高估我国工业全要素生产率的增长。 综上所述，我国工业受能源和环境的双重约束。为此，本文将运用基于熵的非径向、非角度、无方向的数据包络分析（Data Envelopment Analysis， DEA）模型分别建立及包含能源、环境

5、的工业 TFP 生产率模型，并在此基础上对测度出的生产率进行分解，以便更为准确地测度我国工业的全要素生产率。 二、熵权背景下 SBM 节能减排技术效率测度模型的设定 （一）非参数 DEA 生产系统中投入产出要素熵权的引入 本部分在现有的非参数 DEA 中引入了包含能源、环境因素的节能减排技术模型。该模型可以测度在能源、环境的双重约束下我国地区工业的相对效率及生产率。 区别于传统径向的 DEA 方法，为了充分反映节能减排工业生产中各能源投入、污染排放的比例关系，本文引入了熵来为能源、非合意产出加权。熵加权是一种精度更高、客观性更强，能够反映各要素变异程度的赋权方法。现有文献11-12对于熵的运用

6、大多是代替主成分分析，以免出现负值的情况，本文则运用在对于能源、非合意产出的加权上。某一要素的变异程度越大，提供的信息就越多，在整个生产过程中起的作用也就越大，其权重应该越大；反之，某一要素的变异程度越小，其提供的信息量越少，在整个生产过程中的作用也就越小，其权重也应越小。（二）省际工业投入产出系统生产可能性集性质的确定 本文以省为决策单元，以 n 代表省的个数，其中有 m 个投入要素 xRm+，生产了 s1 个合意产出 ydRS1+和 s2 个非合意产出（污染排放）yuRS2+。关于规模报酬方面，我们认为，规模报酬可变更符合我国工业现状，故本文选择规模报酬可变，其生产可能性集可定义为： P=

7、?x（xe，xn，yd，yu）|xeX，xnX，yg （三）节能减排技术模型的构建 公式（1）技术条件下的相对效率可由包含非合意产出的 SBM13求得。定义松弛变量 se-Rm1、sn-Rm2、suRs2 分别为能源投入、非能源投入、非合意产出可能压缩的量。定义松弛变量 sdRs1 为合意产出可能扩张的量。这些松弛变量都大于或等于零。定义Rs1、Rs2 分别为能源投入和非合意产出的熵权重。因此包含能源、环境因素的节能减排相对效率测度模型可通过构建下面的优化问题来解决。 min =1-1m1+m2（m1i=1ise-ixei0+m2i=1sn-ixni0）1+1s1+s2（s1r=1sdrydr

8、0+s2r=1rsuryur0） s.t.C1：xe0=Xe+se- C2：xn0=Xn+sn- C3：yd0=Yd-sd C4：yu0=Yu+suC5：e=1C6：se-0，s-0，sd0， su0，0（2） 其中，C1、C2、C3 和 C4 描述了一个决策单元 （Decision Making Unit，DMU） （xe0，xn0，yd0，yu0） 。 定义一个 t 时期的 DMU 为 （xe0，xn0，yd0，yu0，t） ，在 s 技术条件下的相对相率为 s（xe0，xn0，yd0，yu0，t） 。运用公式（2）所计算出的相对效率，就可以求出从 t 期到 t+1 期的基于 Malmqu

9、ist 指数的生产率（以下简称 EEUSMt，t+1）： EEUSMt，t+1= t（xe0，xn0，yd0，yu0，t+1）t（xe0，xn0，yd0，yu0，t）t+1（xe0，xn0，yd0，yu0，t+1）t+1（xe0，xn0，yd0，yu0，t）（3） 其中，t（xe0，xn0，yd0，yu0，t） 和 t+1（xe0，xn0，yd0，yu0，t+1） 分别测度了 t 和 t+1 时期内的相对效率；t（xe0，xn0，yd0，yu0，t+1） 和 t+1（xe0，xn0，yd0，yu0，t） 则测度了跨期情况下的比较。t（xe0，xn0，yd0，yu0，t+1）t（xe0，xn0，

10、yd0，yu0，t） 和 t+1（xe0，xn0，yd0，yu0，t+1）t+1（xe0，xn0，yd0，yu0，t） 分别表示了 t 和 t+1 时期的效率变化，EEUSMt，t+1 反映了合意产出的扩张及能源、非能源投入、非合意产出的压缩，可以表述从 t 到 t+1 时期生产率的增长、停滞和衰退。如果 EEUSMt，t+11，代表了 DMU 从 t 到t+1 时期的生产率的增长；如果 EEUSMt，t+1=1 或者 EEUSMt，t+11 时，表示从 t 到 t+1 时期技术效率的增加；否则，当 EEUSM_Ct，t+1=1 或 EEUSM_Ct，t+11 时，表示从 t 到 t+1 时期

11、生产前沿面向有更多的合意产出、更少的能源、非能源投入和污染排放的方向移动，即技术的进步；否则，当 EEUSM_Ft，t+1=1 或 EEUSM_Ft，t+11 时，则分别代表了技术的停滞或生产前沿面向有更少的合意产出、更多的能源、非能源投入和污染排放的方向移动，即技术退步。 三、实证分析 （一）数据 文中的投入产出原始数据主要来源于中国统计年鉴 、 中国工业经济统计年鉴以及中国能源统计年鉴等。其中，能源投入为各地区能源总消耗量，非能源投入包括劳动投入和资本投入，劳动为从业人员数，资本为实际固定资产净值，产出的合意部分是实际工业增加值（IAV） ，产出的非合意部分为 SO2 和 COD。 表 2

12、 是本研究中相关数据的描述性统计。就传统投入、产出而言，劳动的平均增长率只有东部地区在增加（2239） ，中、西部都在降低；三大地区的资本平均增长率都超过了 10%；工业增加值的平均增长率则是东部最高（15522%） ，中部最低（13308%） 。就能源、污染排放而言：能源消耗的平均增长率西部最高（8730%） ，中部最低（6766%） ；SO2 的平均增长率中部最高（4348%） ，东部最低（1861%） ；COD 的平均增长率除西部以外均为负值。 （二）实证结果及其分析 本文运用 Matlab 软件编程计算出了节能减排模型及传统模型。 1.熵权的结果 各省际的能源、SO2、COD 熵权的结

13、果如表 3 所示。 2.效率 全国省际工业的效率如图 1 所示。总体来说，1992-2011 年这 20 年中，无论是否考虑能源、污染排放因素，全国工业平均效率值都在下降。而且，节能减排约束下的效率值下降更为明显。这说明，全国省际工业效率水平与前沿面的差异性正在扩大，且未考虑节能减排约束的省际工业效率测度低估了这一差异的严重性。从趋势看，新世纪以来两种效率的差距有扩大的趋势。直到“十一五”期间，国家把节能减排提升到国家层面上来，这种差距才重新缩小。 从省际工业效率均值来看（见图 2） ，上海、江苏、海南、陕西的年度均值为 1，说明这四个省际的工业一直处于节能减排效率的前沿面上。这其中只有陕西省

14、地处中部地区，其他三个均为东部省（市） 。效率排在前十名同时也是效率值高于 07 的省际（上海、江苏、海南、陕西、广东、福建、浙江、云南、北京、天津）基本都来自东部。可见，1992-2011这 20 年中，能源、环境约束下的工业效率受东部地区引领。 3.工业全要素生产率及其分解 根据公式（2） 、 （5） 、 （6） 、 （7）可以计算出节能减排技术模型的 M氏指数及其分解，以下将其结果与传统模型的结果进行对比分析。 传统模型的结果表明，1992-2011 年间全国工业年均增长率为1002%。三大地区中中部地区的工业全要素生产率最高（1253%） ，西部次之（985%） ，东部最低（828%）

15、 。传统模型的分解表明，1992-2011 年间全国工业技术效率的增速为-003%，技术进步为 1025%，因此技术进步是工业全要素生产率的主要来源。样本期内，西部地区的技术效率增速最快（034%） ，东部次之（019%） ，中部则为负值（-023%） 。技术进步方面，中部地区增速最快（1335%） ，西部次之（988%） ，东部最慢（827%） 。可见，中部地区工业全要素生产率最高，且全部来源于技术进步。 1992-2011 年间，全国工业节能减排 M 氏指数的年均变化为 1304%。这其中包括了年均增长 018%的技术效率，以及年均增长 1390%的技术进步。因此，能源、污染排放约束下的工

16、业全要素生产率也主要来源于技术进步。三大地区中，东部地区的工业节能减排技术生产率最高（3019%） ，中部次之（766%） ，西部为负值（-127%） 。样本期内，东部地区的技术效率增速最高（088%） ，西部次之（031%） ，中部为负值（-092%） 。技术进步方面，东部最高（3000%） （这主要得益于广东地区 34808%的超高增速） ，中、西部的技术进步分别为 890%和 042%。样本期内只有陕西的技术进步为负值（-5542%） 。省际层面上（参见图 3） ，广东省的节能减排 M 氏指数最高，为 34；湖南、安徽、内蒙古、四川、山东、福建、北京的 M 氏指数介于 11 到 12 之

17、间；黑龙江、江苏、甘肃、重庆、贵州、云南、山西、浙江、海南、广西、江西、河南、天津、辽宁、河北、吉林、湖北的 M氏指数介于 10 到 11 之间；陕西、上海、宁夏、青海、新疆的 M 氏指数小于 10。 从纵向来看（参见图 4） ，1992-1993 到 2010-2011 的这 19 个区间段全国工业节能减排 M 氏指数基本上总是低于未考虑能源、污染排放的传统 M 氏指数，只是在“十一五”中期的 2007-2008 这一区间段出现了明显地超越现象。但是，如果剔除广东省的数据，全国工业节能减排 M 氏指数基本上还是低于未考虑能源、污染排放的传统 M 氏指数。 30 个省际 19 个区间的节能减排

18、 M 氏指数和传统的 M 氏指数的箱图如图 5 所示。通过观察，我们可以发现大部分的区间内节能减排 M 氏指数的中位数、上四分位数、下四分位数都低于传统的未考虑能源、污染排放的 M 氏指数。我们可以发现，如果忽略能源消耗与污染排放对环境的污染，会高估我国工业的真实全要素生产率。 以上结果与 Kumar 和 Surender14的研究结论一致。这两位学者发现，如果考虑非合意产出，发展中国家的全要素生产率将会降低。但是这与另一些学者15-18的研究结论相矛盾。在这些研究中，考虑非合意产出计算出的全要素生产率高于传统的未考虑非合意产出的情况，即不考虑非合意产出会低估真实生产率。这种矛盾的产生是因为：

19、（1）在环境规制之下，本来用于生产合意产出的资源被用来减少非合意产出的量，而文献15-18中的传统模型忽略了投入对于减少非合意产出的正面作用，因此其结果低估了真实的生产率。然而，在中国很长一段时间以来，我们以发展为优先目标，对环境方面的规制相对于发达国家而言是较少的，执行的也不严格，自然很难产生波特效应，即未能实质影响企业行为。况且，波特效应的产生也是需要时间的。 （2）本文的节能减排模型还增加了非能源投入对能源投入的替代因素，这是上述文献所没有考虑到的。四、结论与启示 本文基于熵加权建立了考虑能源、污染排放的 M 氏指数，并将这一节能减排模型结果与传统模型的结果进行比较，结果发现忽略节能减排

20、因素会高估我国工业的全要素生产率。当然，这其中也有例外，2007-2008 年间广东的节能减排 M 氏指数达到了的 4 70696%，且这一增长全部来自于 技术进步。样本期内，中部地区在节能减排约束下的工业生产率增长最快，且全部来自于技术进步。但是生产前沿面几乎都来自东部，可见东部地区是节能减排技术的引领着。同时，全国省际工业效率水平与前沿面的差异性正在扩大，且未考虑节能减排约束的省际工业效率测度低估了这一差异的严重性。 注释： 2010 年，我国的能源消耗占世界的 203%，而同期美国的能源消耗只占世界的 19%。 参考文献： 1Shi， G.， Bi， J. and Wang， J.Chinese Regional Industrial Energy Efficiency Evaluation Based On a DEA Model of Fixing Non-Energy InputsJ.Energy Policy，2010，38（10）： 6172-6179. 2陈诗一. 能源消耗、二氧化碳排放与中国工业的可持续发展J.经济研究， 2009（4）：41-55. 3陈军，成金华，吴巧生.工业化水平区域差异与中国能源消费J.中