本文采用2003至2017年江西省环境经济数据,选取其中具有代表性的环境污染与经济发展指标,构建了江西省人均废水排放量、人均工业废气排放量、人均工业固体废物排放量与人均GDP的计量模型,从而对江西省环境污染物排放的库兹涅茨曲线特征进行了实证分析。结果表明:江西省2003至2017年人均废水排放量、人均工业废气排放量、人均工业固体废物排放量与人均GDP拟合曲线呈现倒“U”型曲线关系,环境状况总体向好发展。 This paper adopts the environmental economic data of Jiangxi province from 2003 to 2017, and selects the representative indicators of environmental pollution and economic development. Then this paper constructs a measurement model of Jiangxi province's per capita waste water discharge, per capita industrial waste discharge, per capita industrial solid waste discharge and per capita GDP, so as to make an empirical analysis of the Kuznets Curve characteristics of Jiangxi province’s environmental pollutant discharge. The results show that Jiangxi province from 2003 to 2017 per capita waste water emissions, per capita industrial emissions, per capita industrial solid waste emissions and per capita GDP fitting curve shows an inverted “U” curve relationship, the overall environmental situation is good development.
廖恒坚,李昱明,冯亮梅,商荣廷,阎虎勤
厦门国家会计学院,福建 厦门
收稿日期:2020年4月7日;录用日期:2020年4月30日;发布日期:2020年5月7日
本文采用2003至2017年江西省环境经济数据,选取其中具有代表性的环境污染与经济发展指标,构建了江西省人均废水排放量、人均工业废气排放量、人均工业固体废物排放量与人均GDP的计量模型,从而对江西省环境污染物排放的库兹涅茨曲线特征进行了实证分析。结果表明:江西省2003至2017年人均废水排放量、人均工业废气排放量、人均工业固体废物排放量与人均GDP拟合曲线呈现倒“U”型曲线关系,环境状况总体向好发展。
关键词 :环境库兹涅茨曲线,江西省,环境污染,经济发展
Copyright © 2020 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
1995年诺贝尔奖获得者、经济学家库兹涅茨提出收入不均现象会随经济增长先升后降而呈倒U型曲线关系。基于此,1993年帕纳约托首次提出环境质量与人均收入之间的关系曲线——环境库兹涅茨曲线(EKC),该曲线得出环境质量刚开始随着人均GDP的上升而恶化,但人均GDP上升到一定程度后又随收入增加而改善,即环境质量与收入间呈倒U型关系。此后,EKC进入广大学者视野,实证研究层出不穷且研究结论不一,得出倒U型、U型、N型等形状曲线,这也正面印证了环境与经济发展间的关系是个复杂的问题,因为曲线形状可能因对象、环境污染指标、时期的不同而不同。根据相关定义及研究不难得出,EKC是一条经验曲线,描述的是一种普遍规律而不是必然规律,并不能说明人均GDP和广泛存在的各种环境之间存在必然联系。虽然环境具有一定的自动净化功能,但毕竟能力有限,优化环境更有赖于全社会环保意识的提高、限制污染的环境政策的实施,同时还需要技术进步的支持。所以本文将对江西省经济和环境指标间的关系进行拟合,并对拟合结果进行分析,以期对江西省的经济合理发展提供相关建议。
通过查阅相关文献得出,目前关于库兹涅茨曲线的研究多跟环境有关。颜蕾(2014)引用重庆市1995~2012年的经济增长与工业三废的数据构建模型后,得出重庆市的环境库兹涅茨曲线(EKC)呈倒“U”型和倒“N”型,并对该市提出调整能源结构的建议 [
根据环境库兹涅茨曲线的倒“U”型假说的理论、相关文献,以及江西省环境污染现状可以推断出,一个区域的环境污染与经济增长之间并非简单的线性关系,故本文选取EKC的一次曲线模型、二次曲线模型、三次曲线模型分别对自变量(经济指标:调整得到的江西省2003~2017年实际GDP)和因变量(环境指标:工业废气排放量、工业固体废物产生量、废水排放量)进行拟合,其计量模型为:
y = α 1 x + α 0 (1)
y = α 2 x 2 + α 1 x + α 0 (2)
y = α 3 x 3 + α 2 x 2 + α 1 x + α 0 (3)
在上述方程中,y代表江西环境三大指标值,x代表江西省实际GDP, α 0 为常数项, α 1 、 α 2 、 α 3 为待定系数,
待定系数 | EKC曲线形状 | 环境和GDP间的关系 | |
---|---|---|---|
1 | α1= α2= α3 | 水平线 | 不相关,无影响 |
2 | α1> 0, α2= α3= 0 | 线性正相关 | GDP增长,环境恶化 |
3 | α1< 0, α2= α3= 0 | 线性负相关 | GDP增长,环境改善 |
4 | α1> 0, α2< 0, α3= 0 | 倒“U”型 | 环境随GDP增加先恶化后改善 |
5 | α1< 0, α2> 0, α3= 0 | “U”型 | 环境随GDP增加先改善后恶化 |
6 | α1> 0, α2< 0, α3> 0 | “N”型 | 一定范围内,跟第4种情况类似,GDP超出范围后环境进一步恶化 |
7 | α1< 0, α2> 0, α3< 0 | 倒“N”型 | 一定范围内,跟第5种情况类似,GDP超出范围后环境进一步改善 |
表1. 环境库兹涅茨曲线的系数、曲线形状、环境和GDP之间的关系
在了解模型原理后,找出相关数据,然后进行拟合,通过看相关系数R_Square和预测结果来选定最优模型和方程。
本文选取2003~2017年江西省环境与经济数据为样本数据,采用环境库兹涅茨曲线来实证分析江西省的经济增长与其环境污染之间的关系,并预测其在2018~2023年的演变趋势,以寻求二者的平衡点,进而为政府实施宏观经济调控以及制定环境保护政策提供建议,助推江西省经济的高质量发展。
衡量一个地区经济发展状况的指标主要有宏观和微观两大类,国际上通常使用宏观经济指标中的人均GDP作为衡量标准 [
· 使用从世界银行网站1查询的中国2003~2017GDP平减指数,用于计算实际GDP:
实际GDP = 名义GDP/(1 + GDP平减指数/100) (4)
· 从《江西省统计年鉴》中查询到以1978年为基期的地区GDP指数(1978年的GDP指数为100),用地区GDP指数换算得到2003~2017年江西省实际GDP:
实际GDP = 基期江西省实际GDP * (当期地区GDP指数/100) (5)
本文借用Python软件,建立一般多项式,分别以上述两种方法计算得到的江西省2003~2017年实际GDP作为自变量,对相应的环境指标进行预测。表2的结果显示,使用两类数据得出的多项式方程,其拟合度(R2)都在0.9以上,但采用方法一的实际GDP数据时,拟合度更高一些;且考虑到以1978年为基期的地区GDP指数调整得到实际GDP金额与名义GDP金额相差较大,在后续查看预测结果时存在一定不便之处。
年份 | 名义GDP (亿元) | GDP平减指数 (中国) | 地区GDP指数 (江西) | 实际GDP (方法一,亿元) | 实际GDP (方法二,亿元) |
---|---|---|---|---|---|
2003 | 2,087.41 | 2.6054 | 981.9000 | 2736.12 | 854.25 |
2004 | 3,456.70 | 6.9543 | 1,111.5000 | 3231.94 | 967.01 |
2005 | 4,056.76 | 3.9026 | 1,254.9000 | 3904.38 | 1,091.76 |
2006 | 4,820.53 | 3.9279 | 1,409.2000 | 4638.34 | 1,226.00 |
2007 | 5,800.25 | 7.7492 | 1,595.3000 | 5383.10 | 1,387.91 |
2008 | 6,971.05 | 7.7918 | 1,807.4000 | 6467.14 | 1,572.44 |
2009 | 7,655.18 | −0.2105 | 2,046.0000 | 7671.33 | 1,780.02 |
2010 | 9,451.26 | 6.8811 | 2,332.5000 | 8842.77 | 2,029.28 |
2011 | 11,702.82 | 8.0756 | 2,621.7000 | 10828.36 | 2,280.88 |
2012 | 12,948.88 | 2.3351 | 2,910.1000 | 12653.40 | 2,531.79 |
2013 | 14,410.19 | 2.1610 | 3,204.0000 | 14105.37 | 2,787.48 |
2014 | 15,714.63 | 0.7912 | 3,514.8000 | 15591.27 | 3,057.88 |
2015 | 16,723.78 | 0.0627 | 3,834.6000 | 16713.30 | 3,336.10 |
2016 | 18,499.00 | 1.0728 | 4,179.7000 | 18302.65 | 3,636.34 |
2017 | 20,006.31 | 3.8842 | 4,547.6000 | 19258.28 | 3,956.41 |
表2. 2003~2017年江西省实际GDP测算
资料来源:江西省统计年鉴、世界银行。
所以,本文将采用方法一(中国2003~2017GDP平减指数)调整得到的江西省2003~2017年实际GDP作为自变量。
环境指标选取方面,本文选择“工业废气排放量、工业固体废物产生量、废水排放量”三个典型指标,相关数据如表3所示:
年份 | 实际GDP (亿元) | 工业废气排放量 (亿m3) | 工业固体废物产生量 (万吨) | 废水排放量(万吨) |
---|---|---|---|---|
2003 | 2736.12 | 3,202.00 | 6,181.82 | 118,739.00 |
2004 | 3231.94 | 3,972.00 | 6,524.26 | 120,092.00 |
2005 | 3904.38 | 4,378.00 | 7,006.71 | 123,320.00 |
2006 | 4638.34 | 5,096.00 | 7,392.64 | 130,183.00 |
2007 | 5383.10 | 6,103.00 | 7,777.28 | 141,266.69 |
2008 | 6467.14 | 7,456.00 | 8,190.35 | 138,909.00 |
2009 | 7671.33 | 8,286.00 | 8,898.18 | 147,080.54 |
2010 | 8842.77 | 9,812.00 | 9,407.30 | 160,661.00 |
2011 | 10828.36 | 16,102.00 | 11,372.43 | 194,431.55 |
2012 | 12653.40 | 14,814.00 | 11,133.60 | 201,189.71 |
2013 | 14105.37 | 15,573.00 | 11,518.19 | 207,137.56 |
2014 | 15591.27 | 15,613.00 | 10,821.21 | 208,289.04 |
2015 | 16713.30 | 17,054.00 | 10,776.68 | 223,232.28 |
2016 | 18302.65 | 15,162.00 | 12,665.33 | 221,091.88 |
2017 | 19258.28 | 15,065.00 | 12,340.92 | 189,362.32 |
表3. 2003~2017年江西省环境指标数据
资料来源:江西省统计年鉴表格2。
图1显示,2003~2017年江西省实际GDP总体呈直线上升趋势;废水排放量总体呈上升趋势,在2015年达到峰值223,232.28万吨,而后开始逐年下降。
图1. 2003~2017年江西省实际GDP与废水排放量
图2显示,2003~2017年,江西省工业废气排放量与工业固体废物产生量两个指标都总体呈上升趋势。工业废气排放量从2015年开始,逐年下降;工业固体废物产生量在2011~2015年缓慢下降后,又开始呈上升趋势。
图2. 2003~2017年江西省工业废气排放量与工业固体废物产生量
首先,我们用切比雪夫多项式对江西省人均GDP进行拟合,来预测其未来人均GDP的走势,根据表4可得,切比雪夫多项式对人均GDP拟合结果理想。由图3可以看出,江西省人均GDP随着时间的推移将会一直上涨,保持稳定的增长势头。
方程式 | R_Square | R_Sqrt | Adjusted R_Square | Spread of the Error | Root Mean Squared Error | 切比雪夫方程式 | 一般方程式 |
---|---|---|---|---|---|---|---|
切比雪夫 | 0.99 | 0.99 | 0.99 | 1139 | 1019 | GDP = 3352T0+ 31606T1+ 3630T2 | GDP = 29722 + 31606T + 7260T2 |
表4. 人均GDP拟合方程式
图3. 江西省人均GDP切比雪夫模型拟合图
由表5可以看出,在江西省人均GDP与人均废水排放量的3种拟合方程中,三次曲线模型中的R_Square值最大为0.91,说明三次曲线模型在此次拟合中拟合结果比较理想。但是,由于样本数据较小,利用三次曲线模型预测未来人均废水排放量时会出现负值,同时在二次曲线模型与三次曲线模型拟合程度相差不大的情况下,我们在这里选择二次曲线模型作为此次江西省人均GDP与人均废水排放量最终的回归方程:人均废水排放量 = −0.0000167X2+ 1.378X + 17470。
通过对图4、图5和图6的拟合曲线对比,我们选择图5作为最终拟合曲线。江西省人均GDP与人均废水排放量之间呈“倒U型”关系,基本和库兹涅兹曲线相吻合。意味着随着江西省人均GDP的增加,人均废水排放量历经先增加后减少的过程,根据二次曲线模型拟合得出的极值点可知,当江西省人均GDP达到41481元/每人时,人均废水排放量达到46229千克/每人的最大值,之后便随着人均GDP的增长而减少。时间大概对应的是2017年前后达到顶峰。
方程式 | R_Square | R_Sqrt | Adjusted R_Square | Spread of the Error | Root Mean Squared Error | 回归方程 | 极值点 |
---|---|---|---|---|---|---|---|
一次 | 0.86 | 0.93 | 0.85 | 3107.58 | 2893.00 | Y = 0.5963X + 244420 | / |
二次 | 0.91 | 0.95 | 0.89 | 2570.42 | 2299.06 | Y = −0.0000167X2+ 1.378X + 17470 | [41481, 46229] |
三次 | 0.96 | 0.98 | 0.95 | 1691.62 | 1448.61 | Y = −0.000000001785X3+ 0.0001111X2− 1.268X + 32410 | [34656, 47584] [6832, 28363] |
表5. 人均废水排放量拟合方程式
图4. 江西省人均废水排放量与人均GDP线性模型拟合图
图5. 江西省人均废水排放量与人均GDP二次曲线模型拟合图
图6. 江西省人均废水排放量与人均GDP三次曲线模型拟合图
三种拟合方程分别对江西省人均GDP与人均工业废气排放量进行了拟合分析,结果见表6。其中,三次曲线模型的R_Square达到0.97,拟合优度最佳。但是还是因为样本数据较小的原因,利用三次曲线模型预测未来人均工业废气排放量的时候会出现负值,因此,我们在这里仍然是选取拟合效果同样很好的二次曲线模型来预测江西省人均GDP和人均工业废气排放量的回归方程式:Y = −0.00003091X2+ 2.31X − 8341。
根据方程拟合结果以及图7、图8和图9对比可以看出,江西省人均GDP与人均工业废气排放量两者之间呈现倒“U”型曲线特征,基本与环境库兹涅茨曲线相吻合。意味着随着人均GDP的增长,江西省工业废气排放量会先增加、再减少。同时,根据二次曲线模型拟合得出的极值点可知,当江西省人均GDP达到37365元/每人时,人均工业废气排放量达到34814千立方米/每人的最大值,之后便随着人均GDP的增长而减少。时间大概对应的是2015~2016年之间达到顶峰。
方程式 | R_Square | R_Sqrt | Adjusted R_Square | Spread of the Error | Root Mean Squared Error | 回归方程 | 极值点 |
---|---|---|---|---|---|---|---|
一次 | 0.86 | 0.93 | 0.85 | 4420 | 4114 | Y = 0.8483X + 4515 | / |
二次 | 0.95 | 0.97 | 0.94 | 2825 | 2527 | Y = −0.00003091X2+ 2.31X − 8341 | [37365, 34814] |
三次 | 0.97 | 0.98 | 0.96 | 2294 | 1964 | Y = 0.000000001589X3+ 0.00008289X2− 0.05409X + 4966 | [329, 4956] [34438, 36491] |
表6. 人均工业废气排放量拟合方程式
图7. 江西省人均工业废气排放量与人均GDP线性模型拟合图
图8. 江西省人均工业废气排放量与人均GDP二次曲线模型拟合图
图9. 江西省人均工业废气排放量与人均GDP三次曲线模型拟合图
由表7可以看出,以下三种拟合方程中,二次曲线模型与三次曲线模型的拟合度优于线性方程,结合拟合结果倒“U”型图形来看,最终将二次曲线方程设定为人均GDP与人均工业固体废弃物产生量的回归方程,方程表达式如下:Y = 0.000000883X2+ 0.07566X + 997.5。
方程式 | R_Square | R_Sqrt | Adjusted R_Square | Spread of the Error | Root Mean Squared Error | 回归方程 | 极值点 |
---|---|---|---|---|---|---|---|
一次 | 0.89 | 0.94 | 0.88 | 150 | 140 | Y = 0.0339X + 1365 | / |
二次 | 0.94 | 0.97 | 0.93 | 117 | 105 | Y = 0.000000883X2+ 0.07566X + 997.5 | [42839, 2618] |
三次 | 0.94 | 0.97 | 0.92 | 121 | 103 | Y = 0.00000000001687X3− 0.000002091X2+ 0.1007X + 856.3 | / |
表7. 人均工业固体废物产生量拟合方程式
经图10、图11和图12对比可以看出,江西省人均GDP与人均工业固体废弃物之间均呈现环境库兹涅茨曲线关系,显示于倒“U”型曲线。由此可以看出随着人均GDP的增长,工业固体废弃物产生量经历了迅速上升随后下降的过程。根据二次曲线模型拟合得出的极值点可知,当江西省人均GDP达到42839元/每人时,人均工业固体废物产生量达到2618千克/每人的最大值,时间大概对应的是2016~2017年之间。
图10. 江西省人均工业固体废物产生量与人均GDP线性模型拟合图
图11. 江西省人均工业固体废物产生量与人均GDP二次曲线模型拟合图
图12. 江西省人均工业固体废物产生量与人均GDP三次曲线模型拟合图
江西省2003至2017年人均废水排放量、人均工业废气排放量、人均工业固体废物排放量三种指标与人均GDP之间有相关关系且均符合环境库兹涅茨曲线。三种类型的污染物排放量与人均GDP之间呈现倒“U”型曲线关系,表明安徽省环境污染问题正在得到改善。这与江西省公民环保意识的加强、政府部门的重视是密不可分的,针对江西省目前的环境问题提出以下建议:
第一,加强公民的环保意识。提高公民对环境保护的认识水平和自觉程度,法律、法规与媒体宣传是是公民提高环境意识的有力保障,政府应建立和完善相关的法律法规,并充分利用公共媒体加强宣传,增强公民的环保意识。
第二,调整优化产业与能源结构。重工业是江西省环境污染的最大来源,政府需要优化产业结构与能源结构,以改善环境污染。在产业结构方面,有关部门应当支持节能环保的新兴产业发展,淘汰重污染企业,对高耗能、高排放的项目严格把关。此外,优化产业结构不能仅限于淘汰重污染企业,强制性的限产停产,更重要的是发展高新技术,改造传统产业,倒逼传统企业转型升级。在能源结构方面,加快水电、核电、风电等新型能源替代传统能源的步伐,产业结构与能源结构优化双管齐下使环境质量改善。
第三,调整优化工业布局。目前江西省各市常见的环境问题是工业用地与生活用地缺乏合理划分。例如,南昌市经开区麦园垃圾处理场常在夜间焚烧垃圾,气味刺鼻呛人,严重影响周边住户的正常生活。江西其他城市也存在类似问题,可见调整优化工业布局的重要性。相关部门应在工业用地及生活用地之间加设环保隔离带,逐步有序地推进重污染企业的搬迁改造,进而优化工业布局。
第四,推广PPP模式治污。PPP (Public-Private-Partnership)是政府与民间资本对某一具体项目签订相关协议,并相互合作、共担风险、共负责任、共享收益的模式。PPP模式利用民间资本,能够有效地改善政府资金不足的问题,江西省财政收入较低,用于环境保护及环境治理的专项资金匮乏,采用PPP模式使社会资本介入,能够填补环境保护及环境治理的资金缺口 [
本论文得到了厦门国家会计学院2019年“云顶课题:YD20190101Python财务数据分析”项目的支持。
廖恒坚,李昱明,冯亮梅,商荣廷,阎虎勤. 基于库兹涅茨曲线的江西省经济发展与环境污染关系研究 An Empirical Study on the Relationship between Economic Growth and Environmental Pollution in Jiangxi Province Based on Kuznets Curve[J]. 低碳经济, 2020, 09(02): 47-58. https://doi.org/10.12677/JLCE.2020.92006
https://doi.org/10.1080/14786451.2019.1686380
https://doi.org/10.1016/j.scitotenv.2019.135768