1. 引言
党的二十大报告指出,促进区域协调发展,深入实施区域协调发展战略、区域重大战略、主体功能区战略、新型城镇化战略,优化重大生产力布局,构建优势互补、高质量发展的区域经济布局和国土空间体系。而我国幅员辽阔,区域经济与生态环境差异巨大,发展不平衡不充分仍是我国的基本国情,并且随着工业化与城市化的快速推进,生态环境的需求与供给矛盾日益尖锐,生态环境的约束逐渐成为阻碍经济增长的主要因素之一。为了推进区域高质量发展和协同交融,需量化区域发展不平衡不充分的程度,并将经济、资源和环境的空间异质性统筹考虑,以促进各类要素合理流动和高效集聚,推动形成主体功能明显、优势互补、高质量发展的国土空间开发保护新格局。
长江经济带横跨中国东中西三大区域,是中央重点实施的“三大战略”之一,而长三角城市群、长江中游城市群和成渝城市群是支撑和引领长江经济带高质量和一体化发展的重要引擎,也是我国区域发展战略格局中的重要功能区,然而由于长江经济带各省市之间的经济发展基础不同,资源与环境承载力也存在较大差异,城市群发展的异质性问题较为突出。因此,本文将以长江经济带三大城市群为研究对象,对环境基尼系数公式加以创新,衡量经济产出与生态环境的区域协调性。首先,基于区域内污染排放量与空间单元数的讨论,分别构建连续分布函数下的环境基尼系数。其次,基于copula函数建立空间单元污染排放量的总体分布函数,得到空间异质性条件下二元群组混合环境基尼系数公式。最后,运用改进后的二元群组混合环境基尼系数公式对长江经济带的区域协调发展进行评价,并结合空间自相关分析结果,与改进前的环境基尼系数测算结果进行比较,保证新公式评价结果的可靠性。
2. 文献综述
基尼系数是衡量一个国家或地区居民收入差距的常用指标之一。现也常用于分析资源分配的不均衡程度。田凤平 [1] 度量了中国三大城市群经济发展的区域差异,戴平生 [2] 提出一种新的基尼系数计算方法分析了我国教育经费配置的公平性,Papathanasopoulou [3] 运用基尼系数研究化石资源的均衡性问题,张龙耀 [4] 揭示了农村数字普惠金融发展水平的空间差异,于伟 [5] 测算了八大城市群内部生态效率的区域差异,王争磊 [6] 测度生态安全水平的区域差异与来源。基尼系数的应用领域越来越广,但基尼系数的计算方法不同,计算结果差异也较大,目前大多研究是基于非参数方法计算的,但Andrea [7] 等人认为具有厚尾特征的数据在使用非参数方法计算基尼系数时会出现向下偏差。
长期以来,我国的区域发展一直以过度的空间集聚为特征,不同区域的人口数量与土地面积不同,经济发展基础、资源与环境承载力也存在较大差异,造成我国经济与生态环境数据严重厚尾,考虑到区域发展的异质性问题较为突出,需要计算总体基尼系数。胡志军 [8] 指出连续分布函数方法比离散方法更适合于计算分组数据集的基尼系数。程永宏 [9] 创造性地利用城乡收入分布函数,基于洛伦兹曲线的参数方程,推导出新的基尼系数计算公式,后又推广至多群组混合基尼系数 [10] 。但人口数与区域空间数的量级存在较大差异,不同数据的分布函数也存在差异,不能直接将收入分配基尼系数应用至环境基尼系数领域中。艾小青 [11] 也研究了多群组的混合基尼系数,还实现了有限总体 [12] 与信息不完全下 [13] 的财富基尼系数的估算。多群组混合基尼系数的研究多是基于各子群体独立的前提下计算的总体收入分布函数,忽略了各子群体的分布函数存在不独立关系的可能,不同的总体分布函数会对基尼系数的计算结果产生影响。
为适应区域经济和生态环境的厚尾数据特征,在参数式基尼框架下构建了多群组混合环境基尼系数,并且考虑到各子群体分布函数的不独立性,建立了总体污染排放分布函数新的公式,构造了空间异质性条件下新的混合环境基尼系数公式,对长江经济带三大城市群产业经济和生态环境的区域协调性进行评价分析。
3. 研究设计
3.1. 基于连续分布函数的环境基尼系数的构建
3.1.1. 基于连续分布的环境基尼系数
环境基尼系数主要以资源消耗量与污染排放量为分配对象(下文主要介绍污染排放量)。记空间单元的污染物排放随机变量为S,
,其排放分布记为
,x则可表示为污染排放不超过s的空间单元数(设为N)占整个区域空间单元总数的比重,即
。y表示为这N个空间单元的累积排放占整个区域污染排放总量的比重,可推导出用参数方程表示的洛伦兹曲线:
(1)
其中,N0为整个区域内空间单元总数,S0为整个区域污染排放总量。
可得环境基尼系数的基本公式:
(2)
3.1.2. 不同样本间隔与样本大小下的环境基尼系数
由前文推导的环境基尼系数的基本公式可以看出,环境基尼系数与分配对象的分布以及空间单元的数量有关。若要计算环境基尼系数,还需得到污染排放分布函数的具体形式、区域中所有空间单元的最低污染排放量a与最高污染排放量b、区域污染排放总量S0,其中分布函数的具体形式可由数据拟合得到,为方便计算,其余问题将进行如下讨论:
1) 空间单元的污染排放量上下限
第一,若各空间单元污染排放量均存在差异,所有空间单元的污染排放S在
内连续。污染排放为a的空间单元数为1,则相应地,其概率为
;排放小于a的空间单元数为0,相应的概率为0;排放大于a的空间单元数为
,相应的概率为
。因此,污染排放分布函数
在
处有一个跳跃间断点,同理,在
处也存在不连续,即:
(3)
(4)
可由拟合的污染排放分布函数
推导污染物排放上下限的估计值,即:
(5)
第二,若各空间单元污染排放量不存在差异,则达到最低与最高污染排放量的空间单元数不唯一,记为
,其余空间单元的污染排放仍在
内连续。因此,污染排放为a的概率为
,排放小于a的概率为0,排放大于a的概率为
。同理,污染排放分布函数
在
与
处均存在不连续,即:
(6)
(7)
则污染物排放上限和下限的估计值为:
(8)
在计算收入分配基尼系数时,样本量
,分布函数在不连续点的左右极限非常接近。然而在计算环境基尼系数时,一个区域内的空间单元数是有限的,区域经济产出和生态效益分布函数在不连续点的左右极限不相等,采用不同的近似方法会对环境基尼系数的计算结果产生不同的影响。
2) 区域污染排放总量
第一,若区域内空间单元数量足够大,即
和
,那么污染排放分布函数与收入分布函数的差异则会非常小,主要差异体现在分布函数的区间上。此时,区域整体污染排放S0:
(9)
将式(3)、(4)、(9)代入式(2),则可得到环境基尼系数:
(10)
第二,若区域内空间单元数量有限,不连续点的左右界限对计算环境基尼系数的影响不可忽略。此时,区域整体污染排放S0:
(11)
或者
(12)
因此,环境基尼系数:
(13)
(14)
3.2. 空间异质性条件下二元群组混合环境基尼系数的构建
3.2.1. 空间单元污染排放量的总体分布函数
在不能判断变量间存在何种相关关系的情况下,Copula函数可将多元随机变量的联合分布函数与边缘分布函数进行连接,通过变量间相关关系的具体结构,提供了更多相关性变化的细节,能更加充分地解释变量间的相关性,因此使用Copula函数构造总体污染排放分布函数。
根据Sklar定理,假设有随机变量
,其边缘分布函数为
,他们的联合分布为
,则对任意
,一定存在一个Copula函数C,使得:
(15)
若
与
连续,则Copula函数C唯一存在 [14] [15] 。
研究中常见的Copula函数主要为椭圆Copula与Archimedean Copula函数族,常用Copula函数类型、分布函数公式及其参数范围具体见表1 [16] 。其次需对Copula函数进行参数估计和模型效果评价,根据评价指标选择最优Copula函数描述变量间的相关结构,本文选择极大似然估计法(ML),并以与经验Copula函数的平方欧式距离最小作为检验标准。
假设研究区域空间范围内存在M个不同的群组,代表不同经济水平与不同生态承载力的空间单元,随机变量S表示区域整体污染排放,其分布函数为
,其中
,
,
。第m个空间排放记为
,
,所以得到区域总体污染排放分布函数为:
(16)
使用Copula函数建立总体的污染排放分布函数不依赖于“各空间单元分布函数相互独立”的假定,这也与根据程永宏城乡混合居民收入分布函数推演的,以空间单元数占该区域空间单元总数比重为权重构造的总体污染排放分布函数表达式不同。
3.2.2. 二元群组混合环境基尼系数
式(16)是在空间异质条件下的区域污染排放的多组总体分布函数,代入式(10)、(13)与(14),可以得到在无限样本容量下,具有空间异质性的多群组环境基尼系数:
![](Images/Table_Tmp.jpg)
Table 1. Basic characteristics of Copula functions
表1. Copula函数基本特征
(17)
以及,在有限样本下的多群组环境基尼系数:
(18)
或者
(19)
二元群组环境基尼系数公式中三个重要元素的数值处理3.1.2节相同。首先,将不同功能区的空间单元分组,然后根据相应的分布函数拟合每组的排放量数据,再使用Copula函数建立总体的污染排放分布函数。其次,基于分布函数计算各空间单元的污染排放量上下限以及相应的区域总污染排放量。此外,在计算总体环境基尼系数的过程中,需注意各组污染排放上下限的交叉点以及总体污染排放的上下限。
考虑到分布函数的独立性以及区域空间单元的有限性,式(17)~(19)不同于传统的环境基尼系数。传统环境基尼系数是基于参数式基尼框架下构建的多元群组混合环境基尼系数,其中总体污染排放分布函数是以空间单元数占该区域空间单元总数比重为权重构造的,具体的环境基尼系数计算公式为:
(20)
其中,
表示第m个空间单元数占该区域空间单元总数的比重,
表示第t年第m个区域污染排放分布函数。
当各区域排放分布函数不独立时,仅以区域空间单元数占总单元数的比例为权重构造的总体分布函数不足以解释各区域间的相关性。在小样本情况下,计算多群组环境基尼系数不能忽略上下限数据分布的不连续性,各指标的最值对数据分布特征的影响也会更为显著,多群组环境基尼系数则可以更好地衡量区域发展的平衡程度。
4. 长江经济带经济产出与生态环境区域协调性发展评价
应用上述方法,对我国长江经济带的经济产出和生态环境区域协调性进行评价分析。将长江经济带以城市群进行分组,包括长三角、长江中游、成渝三个城市群,其中长三角共计26个城市,长江中游共计28个城市,成渝共计16个城市。本节整理了2005~2020年各城市的GDP (亿元)、社会用电量(亿千瓦时)、工业废水排放量(万吨)和工业SO2排放量(吨)的数据,数据均来源于2005~2020年的《中国统计年鉴》、《中国城市统计年鉴》以及各省市统计年鉴。以区域GDP作为经济产出指标,并且以2005年为基期,对各城市其余年份的GDP进行平减处理;以市社会用电量作为能源消耗的指标,其中2005~2016年使用市辖区用电量,2017~2020年使用全市用电量;以SO2排放量和工业废水排放量作为环境污染的指标,其中SO2排放代表大气污染,工业废水排放代表水资源污染。
4.1. 构建区域经济产出、能源消耗和污染排放的分布函数
4.1.1. 分组数据拟合区域工业产出、能源消耗和污染排放的分布
使用(对数)正态分布族、Logistic分布、Gamma分布以及Weibull分布等14种统计分布,分别对2005~2020年长江经济带三个城市群四类经济和生态效益指标数据拟合分布函数,并通过K-S检验,选取P值最大的模型作为最佳拟合分布。从拟合结果看,各类指标数据均呈现明显的“厚尾”特征,并且不同城市群间存在显著的时间和区域差异,经济和生态效益指标间也存在一定的差异。以对数正态分布拟合结果为例,如图1所示,GDP数据呈现明显的“厚尾”特征,能源消耗与污染排放数据呈现“峰值”与“厚尾”特征,随着时间的推移,GDP与能源消耗数据的“峰值”呈下降趋势,污染排放数据的“峰值”呈上升趋势;GDP数据的“厚尾”特征依然存在,能源消耗数据的“厚尾”特征逐渐增加,污染排放数据的“厚尾”特征逐渐减少。
![](//html.hanspub.org/file/18-2623541x69_hanspub.png?20231016170312373)
Figure 1. Histogram of grouped data of various indicators under lognormal distribution
图1. 对数正态分布下的各类指标分组数据直方图
4.1.2. 构建区域工业产出、能源消耗和污染排放的联合分布
本节首先使用Spearman和Kendall秩相关系数来判断长江经济带城市群间区域工业产出、能源消耗和污染排放指标的相关性。从图2的结果看,两种相关系数均显著大于0.8,说明不同城市群间存在高度相关,这也说明两者并不独立,不能使用传统的参数式环境基尼系数计算,应使用Copula函数解释变量间的非线性关系。
本节选择GumbelCopula、ClaytonCopula与FrankCopula,拟从这三个常用的Copula函数中选择最佳Copula函数 [17] ,建立长江经济带区域工业产出、能源消耗和污染排放的总体分布函数。本节采用极大似然估计方法进行参数估计,并计算平方欧氏距离,选取距离最小值作为最佳Copula函数,部分结果见表2,总体的污染排放分布函数表达式为:
(21)
(22)
(23)
其中,
为最优Copula函数的参数,数值见表3说明;u、v表示两个子群体的污染排放分布函数。
![](//html.hanspub.org/file/18-2623541x74_hanspub.png?20231016170312373)
Figure 2. Correlations among variables of urban agglomerations in the Yangtze River Economic Belt
图2. 长江经济带城市群变量间相关性
4.2. 产业经济、生态环境区域协调性的混合环境基尼系数评价
将式(21)~(23)代入式(18),以最优函数为GumbelCopula为例,得到二元群组混合环境基尼系数的计算公式:
![](Images/Table_Tmp.jpg)
Table 2. Fitting results of the pollution emission distribution function for the midstream and downstream of the Yangtze River Economic Belt urban cluster
表2. 长江经济带中下游城市群污染排放分布函数拟合结果
![](Images/Table_Tmp.jpg)
Table 3. Optimal copula functions and their parameters for various indicators of urban agglomerations in the Yangtze River Economic Belt
表3. 长江经济带城市群各类指标的最佳copula函数及其参数
*表示最优函数为Gumbel Copula时的参数,**表示最优函数为Clayton Copula时的参数,***表示最优函数为Frank Copula时的参数。
(24)
其中,
为第t年的混合环境基尼系数,
为第t年最优Copula函数时的参数,
、
分别表示第t年子群体的污染排放分布函数。
根据式(24)分析了空间异质性对长江经济带区域协调发展的影响及其动态演变趋势,计算的环境基尼系数结果如图3。一般认为,环境基尼系数小于0.2时,说明区域发展高度平均,0.2~0.3时较为平均,0.3~0.4时相对合理,0.4~0.5时差异较大,大于0.5时差异悬殊。
从区域发展看,空间异质性对于长江经济带区域平衡发展影响较大。长三角城市群与长江中游、成渝城市群的GDP基尼系数均超过了0.4,而长江中游城市群与成渝城市群的基尼系数在0.3~0.4之间;2005~2016年,长江经济带市辖区的能源消耗基尼系数已接近0.5,2017~2020年长江经济带全市能源消耗的基尼系数也接近0.4;环境方面,长江中上游的基尼系数均在0.35~0.45之间,中下游与上下游的工业废水基尼系数均超过了0.4,上下游的SO2基尼系数已超过0.45,中下游的SO2基尼系数在0.3~0.4之间。长江经济带区域经济发展与生态空间协调程度弱,主要是长三角城市群与其他城市群发展差异较大,长江中游与成渝城市群之间略有差异。
从演变趋势看,长江经济带区域发展差异逐渐缩小,区域经济发展与生态空间协调性增强,但这四类指标基尼系数仍大于0.3,表明区域发展不平衡问题依然存在。长江中下游与上下游的GDP基尼系数有小幅度的下降趋势,中上游的基尼系数略有上升;能源消耗的基尼系数整体呈下降趋势;环境方面,长江中上游的工业废水基尼系数变化不大,其他虽有波动但整体有下降趋势。
4.3. 区域协调发展时空联系
本节将计算全局莫兰指数与传统环境基尼系数。首先,计算全局莫兰指数以研究长江经济带四类指标的空间自相关性。其次,计算传统的参数式多群组环境基尼系数,最后,结合莫兰指数,将传统的与改进后的环境基尼系数进行比较分析。
4.3.1
. 区域协调发展空间自相关性
莫兰指数是度量空间相关性的一个重要指标,在一定程度上可以解释空间异质性对区域协调发展的影响 [18] [19] 。本节计算全局莫兰指数以表征空间聚集效应 [20] ,其取值范围
为,若数值小于0,说明区域协调发展能力存在空间负相关,即相邻区域之间协调发展能力有相异属性,数值越小说明空间离散性越大;若数值大于0,说明区域协调发展能力存在正相关,数值越大说明空间聚集性越强;若数值等于0,说明属性随机分布,不存在空间自相关性。首先构建长江经济带三个城市群的空间权重矩阵,再计算全局莫兰指数,结果见表4。
![](//html.hanspub.org/file/18-2623541x81_hanspub.png?20231016170312373)
Figure 3. Time series of BVE-Gini of GDP, Energy, industrial wastewater, and SO2 for urban clusters in the Yangtze River Economic Belt
图3. 长江经济带城市群GDP、Energy、工业废水和SO2的环境基尼系数时序图
长江经济带的产业经济与生态空间相关性显著。长江中下游城市群四类指标的全局莫兰指数大多通过了显著性检验,且数值均在0.1~0.6之间,说明长江中下游出现了正相关的聚集效应,即具有相同属性的相邻空间单元更容易聚集,表现为高–高聚集或低–低聚集,这更易导致区域发展不平衡。长江中上游的全局莫兰指数虽然小于0,但大多没有通过显著性P值检验,说明长江中上游虽然存在负向的空间相关性,但不够显著;长江上下游GDP与工业废水的全局莫兰指数存在显著正向的空间相关性,说明长江上下游的GDP与工业废水数据更易导致区域不平衡发展。
莫兰指数显著不等于0时,易有空间聚集或离散现象,更易导致不平衡,随着数值的减小,区域发展协调性增强,但仍存在发展不平衡的问题。随着时间的推移,长江中下游与上下游GDP的空间聚集程度呈现先加强后减弱的趋势,但变化非常微弱;工业废水的莫兰指数先后经过两次上升再下降,但整体呈上升趋势,说明空间聚集性增强,发展不平衡程度加剧,区域发展协调性减弱;SO2虽有波动,但最终呈下降趋势,说明空间聚集性减弱,区域发展协调性增强。
![](Images/Table_Tmp.jpg)
Table 4. Spatial autocorrelation of urban agglomerations in the Yangtze River Economic Belt
表4. 长江经济带城市群空间自相关性
*表示在0.05水平显著,**表示0.01水平显著,***表示0.001水平显著。
4.3.2
. 与传统环境基尼系数比较
本节将由程永宏的收入分配基尼系数推导所得的环境基尼系数作为传统的基尼系数,将不依赖于“各子群体分布函数相互独立”假定的环境基尼系数作为改进后的基尼系数。
将两者进行比较,结果如图4所示。结果表明,传统环境基尼系数的数值普遍低于改进后的环境基尼系数,采用传统环境基尼系数方法会低估长江经济带区域经济与生态的失衡。与传统环境基尼系数相比,本文提出的环境基尼系数方法对不独立数据的关系更敏感,能更好地拟合长江经济带区域经济与生态环境数据的特征,并更好地反映分组数据下的区域不平衡问题。
![](//html.hanspub.org/file/18-2623541x82_hanspub.png?20231016170312373)
Figure 4. Time series of improved and traditional environmental Gini coefficient of GDP, Energy, industrial wastewater, and SO2 in the Yangtze River Economic Zone (2005~2020)
图4. 不同计算方法下长江经济带GDP、Energy、工业废水和SO2的环境基尼系数时序图(200~2020年)
![](//html.hanspub.org/file/18-2623541x83_hanspub.png?20231016170312373)
Figure 5. Relationship between the Gini coefficients of mixed environments in the Yangtze River Economic Belt
图5. 长江经济带混合环境基尼系数的关系
结合全局莫兰指数与两种方法下计算的基尼系数结果分析,面对不独立数据时,传统的环境基尼系数衡量数据间不平衡的准确度降低。以长江中下游SO2为例,2016~2019年的长江中下游的SO2数据,莫兰指数有非常明显的下降,从3.66下降到0.183,说明空间聚集程度有明显的减弱趋势;改进后的基尼系数同样有下降的趋势,从0.388下降到0.296,说明区域不平衡现象在逐渐改善;但传统环境基尼系数下降趋势不明显,甚至有上升趋势,从0.309下降到0.289再上升至0.299,说明传统的环境基尼系数有向下的偏差,对于存在相关关系的数据间不平衡的估计不够准确。同样,2013~2019年,长江经济带中下游工业废水的全局莫兰指数从2013年的0.412上升至2017年的0.521,再下降至2019年的0.508,改进后的基尼系数也是同样的趋势,说明区域不平衡存在先增强后减弱的现象,但传统的环境基尼系数却是从2013年的0.341上升至2015年的0.365,再下降至2018年的0.350,后又上升至0.392,趋势明显不一致。其他城市群的传统环境基尼系数也存在类似低估的问题,便不再赘述。
4.4. 区域经济和生态平衡发展评价
由前面的分析可以看出长江经济带区域经济与生态环境发展存在明显的不平衡问题,区域发展协调性弱。为了研究分组数据下四类指标的基尼系数之间的演化关系,本节将以能耗基尼系数为横轴,经济产出基尼系数和污染排放基尼系数为纵轴,表示“清洁生产效率”的演变趋势如图5。
区域间清洁生产效率的差异具有阶段性特征。第一阶段,2005~2010年,长江经济带区域经济产出和污染排放的不平衡现象显著,说明各城市群之间清洁生产效率呈现相对较高的不平衡。第二阶段,2011~2015年,长江经济带中下游与上下游经济产出的不平衡性仍然较高,而中上游的经济产出较平衡;污染排放不平衡呈现持续增长的趋势。第三阶段,2016~2020年,长江经济带区域经济产出的不平衡仍在上升,区域间清洁生产效率平衡与上一阶段不同,污染排放不平衡呈下降趋势。
5. 结论与建议
本研究以长江经济带城市群为研究对象,为适应区域经济和生态环境数据明显的厚尾特征,基于参数式基尼框架下构建了多群组混合环境基尼系数,并在此基础上,考虑到各子群体的污染排放分布函数不独立,基于Copula函数建立了总体污染排放函数新的公式,构造了空间异质性条件下新的混合环境基尼系数公式,对2005~2020年长江经济带三大城市群的产业经济和生态环境区域协调性进行评价分析。与传统环境参数式基尼系数相比,改进后的环境基尼系数能更好地解释子群体间的相关性,计算结果也能看出:对于存在多个经济基础差异和多个生态承载力水平的群组,传统的环境基尼系数在评估区域发展的经济和生态空间平衡时会存在向下的偏差,改进后的环境基尼系数有更高的准确性。将2005~2020年的结果按照国民经济的五年规划期进行分组,结果表明:从“十一五”至“十三五”期间,长江经济带的产业经济与生态环境区域协调性整体呈上升趋势;产业经济、能源消耗与环境水污染方面,长江中上游与下游发展差异较大,长江中上游之间发展略有差异,大气污染方面,长江中下游与上游发展差异较大,长江中下游发展差异较小。“十一五”期间,长江经济带产业经济与生态承载力的区域协调性持续增强,而中上游经济产出与能源消耗的区域协调性下降;“十二五”期间,长江经济带水污染的区域协调性持续减弱;“十三五”期间,环境污染区域协调性仍在减弱。
基于以上的结论,据此建议,在“十四五”时期,重点发展长江中上游,充分发挥南京、武汉、重庆、成都等城市的核心载体作用,促使生产要素向周边城市流动,缩小长江经济带发展差异。对长江经济带的生态环境采取差异化管理,针对水污染严重的长三角城市群与大气污染严重的成渝城市群,采取更为精细化、科学化的管理措施。充分发挥互联网、大数据、人工智能等技术,实现各环节的节能减排,全面推动长江经济带工业产业绿色低碳转型,以形成上中下游优势互补、协作的互动格局,促进长江经济带成为引领我国经济高质量发展的生力军。
NOTES
*通讯作者。