1 引言
消费可拉动经济增长并将带动投资相应的增长，是社会再生产的终点和新的起点。它使经济增长的动力逐渐由供给转向需求，由生产转向消费，从而带来经济持久增长。消费结构通过其所创建出的需求结构影响产业结构，并在一定程度上决定着产业结构的变化。在当前我国经济增长方式转变和产业结构调整的背景下，促进消费需求的持续稳定增长是保持地区经济健康发展的重要因素，是推动社会经济结构调整的重要力量。
区域消费水平的研究内容广泛，涉及消费水平的测度和消费结构的优化调整、消费行为特征及关联影响因素、消费动因等相关方面。改革开放以来我国区域消费水平的研究主要包括省区间城镇居民消费结构的差异、变动趋势及消费水平的格局变化。但其多从时序特征展开，方法上多采用传统的统计学研究方法，缺乏对全国总体消费水平格局的空间格局与空间关联性研究。
伴随着社会主义市场经济的建立，区域之间的相互作用日益明显，这势必会带来相邻区域消费水平的相互影响，形成趋同的消费水平空间格局。本文采用马尔科夫链方法，通过构建改革开放以来各省人均消费水平的空间马尔科夫转移概率矩阵，对我国消费水平的时间与空间演变特征进行综合考虑，探讨区域消费水平的趋同过程与不同区域背景之间的内在联系，揭示我国消费水平演变的时空特征与规律。
2 数据与方法
2.1 数据源分析
社会消费品零售额是表现区域消费需求最直接的数据，反映了一定时期内人民物质文化生活水平，及社会商品购买力的实现程度与零售市场的规模状况，其中，人均消费额是整个社会消费水平的综合反映。本文选取全国 31 个省区市（不含港、澳、台）人均社会消费品零售额作为区域社会消费水平趋同分析的指标，主要数据利用《新中国五十年统计资料汇编》和历年《中国统计年鉴》计算所得。参照相关研究，将我国人均消费额划分为以下4 种类型：（1）低水平，人均社会消费品零售总额低于全国平均的75%；（2）中低水平，人均社会消费品零售总额介于全国的75%~100%之间；（3）中高水平，人均社会消费品零售总额介于全国的100%~125%之间；（4）高水平，人均社会消费品零售总额高于全国平均的125%。
2.2 研究方法
2.2.1 加权变异系数与泰尔指数
区域差异的测度指标较多，虽然各指标形式不一，但其原理基本一致。区域社会消费品零售额与地区人口规模密切相关，因此，本文选用加权变异系数与泰尔指数来计算我国地区间社会消费水平的差异，反映了人口数量对区域消费水平差异的影响作用，其计算公式为：【1-2】
式中，CV 为我国省区间消费水平的加权变异系数；T 为消费水平差异的泰尔指数；Xi为第 i 个省份的人均社会消费品零售额；Xˉ 为 Xi的均值；Pi第 i 个省份的人口数量；P 为全国人口总量。CV 值和T 值越大，则表明区域间社会消费水平差异越大。
2.2.2 Mann-Kendall 方法
Mann-Kendall 作为一种时间序列突变点划分的非参数统计检验方法，不需要样本遵从一定的分布，不受少数异常值的干扰，具有检测范围宽、人为影响小、定量化程度高的优点。其是根据时间序列{x1，x2，…，xn}构造标准正态分布的统计变化量 UFi和UBi，若UF或UB的值大于0，则表明时间序列呈上升趋势；小于0则表明呈下降趋势。在给定显着性水平临界直线下，若UF或UB的变化出现交点，且交点在临界直线内，其交点便是寻找的突变点。
2.2.3 马尔科夫链
马尔科夫链是一种时间和状态均为离散的随机过程，通过将连续的数据离散化为k种类型，分别计算相应类型的概率分布及其时间变化，以近似逼近事物演变的整个过程。通常，将 t 时刻的属性类型的概率分布表示为一个1×k的状态概率向量Pt，记为Pt=[P1, t，P2, t…，Pk, t]，而不同时刻事物状态转移过程可以用一个k×k的概率矩阵M表示。
概率矩阵M中的元素 mij表示 t 年份属于类型i 的区域在 t+1年份转移到 j 类型的一步转移概率值，其计算公式为：【3】
式（3）中，nij表示在整个研究期间内，由 t 年份属于 i 类型转变为 t+1年份 j 类型的空间单元数量总和；ni是研究时段内所有年份 i 类型的空间单元数量总和。如果某区域的属性类型在初始年份为i ，在下一年份保持不变，则区域类型转移状态平稳；如果属性类型提高，则区域向上转移；否则，区域向下转移。
2.2.4 空间马尔科夫链
空间马尔可夫转移概率矩阵以区域 i 在初始 t年份的空间滞后类型（k 个类型）为条件，将传统的马尔可夫矩阵分解为 k 个 k×k 条件转移概率矩阵。对第 k 个条件矩阵而言，元素 mij(k) 表示以区域在 t 年份的空间滞后类型 k 为条件，该年份属于类型 i 转变为下一年份类型 j 的一步空间转移概率。一个地区的空间滞后类型由属性值的空间滞后算子来确定。空间滞后算子是该区域邻接地区属性值的空间加权平均，通过区域属性值值向量（X）和空间权重矩阵（W）的乘积来计算，具体计算公式为：【4】
式（4）中，Lag表示空间滞后算子；Xi为区域单元的属性值；Wij为空间权重矩阵W的第i行第j的元素。本文采用邻接原则进行权重定义。空间马尔可夫链将传统的马尔可夫链与“空间滞后”相结合，不仅考虑区域发展的时间因素，而且也考虑其空间背景。通过比较马尔科夫矩阵元素与相应的空间马尔科夫矩阵对应元素大小，可了解一个区域消费水平向上或向下转移的可能性与邻域的关系，从而有效揭示区域消费水平的时空演变过程与区域背景间的内在联系。
3 1978-2011年间中国区域消费水平时空演变分析
3.1 1978 以来我国区域消费水平差异变化
利用加权变异系数和泰尔指数分别计算了改革开放以来我国省区间社会消费水平的差异变化（图1）。【图1】
由图1可看出，1978-2011 年间我国区域消费水平的加权变异系数和泰尔指数变化趋势大致相同。整体来看，区域消费水平差异呈现先增大后缩小变化的趋势，且具有明显的阶段性特征。具体来说，1978-1993 年间我国区域消费水平差距不断拉大，1993年差异值最大；1993-2008 年区域消费水平差异呈波动下降态势；2008年以后表现出迅速下降趋势，说明区域间消费水平差距不断缩小。
3.2 区域消费水平差异变动阶段
1978-2011 年间我国区域消费水平差异变动表现出从低到高，然后又到低的状态，表现为不同的阶段变化过程，呈现出明显的波动变化态势。为此，利用Mann-Kendall法判断我国区域消费水平的差异的时间突变点，划分其具体的阶段特征。以1978-2011 年间区域消费水平差异的加权变异系数值和泰尔指数值为时间序列，分别绘制出改革开放以来我国区域消费水平差异的加权变异系数和泰尔指数的Mann-Kendall统计量曲线（图2、3）。图2中UF和UB的交点位置大致位于1985年，图3中UF 和 UB 的交点位置大致位于 1986 年，因此，综合考虑选取突变点为1986年，并将改革开放以来我国区域消费水平差异变动划分为 1978-1986 年和1986-2011 年 2 个阶段。【图2-3】
3.3 1978 年以来我国区域消费水平的时间演变特征
根据我国区域消费水平差异变动的阶段划分，利用马尔科夫链分别计算 1978-1986 年和 1986-2011 年我国区域消费水平类型的马尔科夫转移概率矩阵（表1）。【表1】
由表1知，2个时段内主对角线上的元素值相对较高，表明区域单元在经济发展过程中，初期人均水平类型具有较高的稳定依赖性。（1）所有对角线元素数值均都大于非对角线元素值。其中，对角线元素最小值为0.8315，最大值为0.9933。这表明随着区域自身经济的发展，区域消费水平虽有所提高，但其在区域总体消费水平的分布状态与初始年份相比没有较大变化，其保持初始状态水平的可能性至少为85.71%。（2）区域消费水平类型间的转移概率非常小。非对角线元素最大值为0.0977，远小于对角线元素值。而且，大于0的元素值基本均位于对角线两侧；其中1978-1986 年间存在由低消费水平向高消费水平转移的概率为0.0093，表明区域经济发展过程中，区域消费水平实现跨越式发展的可能性是存在的。（3）1986-2011 年间的低消费水平及高消费水平地区保持稳定状态的概率值较1978-1986 年间的值更大，说明 1986-2011 年间区域消费水平的“俱乐部趋同”现象更显着，即区域消费水平具有较强的内生演化特性。
综上分析，无论是低消费水平地区还是高消费水平地区，均沿着最初类型稳步演进，但是，不同时段内区域消费水平类型的数量存在着增加与减少。
3.4 1978年以来我国区域消费水平的空间变化特征
长期以来，受各地区人均GDP、居民人均收入水平、交通通达水平和第三产业产值比例等因素的影响，我国消费水平存在着较大的区域差异性，表现为不同的空间分布特征。根据不同时段的地区消费水平类型，利用 ArcGIS 对 1978-1986 年和1986-2011 年 2 个时段内区域消费水平类型的空间变化进行可视化（图4（a）、（b））。【图略】
从图4可看出：（1）1978-1986 年间区域消费水平类型转移空间变化不显着，其中，江苏、浙江和广东3个省份的消费水平发生了向上转移，总体分布于我国东部沿海地区；区域消费水平向下转移的省份同样有3个，分别是内蒙古、新疆和青海，主要集中于我国西北部地区；而其他大部分区域消费水平类型未发生转移，这些地区分布集中且大部分位于我国的中部。该时期向上与向下转移的区域总体上呈现出集中分布的特征，分别集中于东部沿海和内陆西北，反映出1978-1986 年间我国区域消费水平存在着集聚与趋同现象。（2）1986-2011 年间区域消费水平类型转移具有较大的空间变化，其中，向上转移的区域数量增加至12个，东部沿海地区均呈上升态势，整体空间分布呈现出由沿海向内陆中部地区邻域扩展的趋势；向下转移的区域数量增加至7 个，空间上集中分布于西部的新疆、西藏及青海，零散分布于中部的黑龙江、山西、宁夏、海南；其他大部分未发生类型转移的地区分布集中分布于我国的中部。这一时期向上和向下转移的区域数量有所增多，空间分布更为集中分布，说明1986-2011年期间存在着更明显的趋同现象。（3）2个时段内我国东部沿海地区均不存在消费水平类型向下转移的现象，至少为平稳演进状态，这与地区的经济发展水平等因素密切相关。
3.5 空间作用背景下的我国区域消费水平演变特征分析
3.5.1 区域消费水平的时空转移特征
传统的马尔科夫链可计算区域消费水平的时间演化特征，但无法探测邻域区域消费水平对本地区的影响作用，随着区域间基础设施的完善，各种联系通道的建立使得商品流通，跨区消费、消费模仿等消费行为的区域间相互作用日益明显，区域消费水平的变化越来越受到周边区域的影响。为此，考虑区域背景作用的影响，以人均消费额的空间滞后为条件，分别构建1978-1986 年和 1986-2011 年的空间马尔科夫转移概率矩阵，探究区域间消费水平的相互影响与相互作用方向（表2）。【表2略】
由表2可知，邻域地区消费水平状态对区域自身消费水平的影响作用是存在的，但2个时段具有显着不同的影响作用，其中1986-2011 年间地区消费水平相互作用更明显。具体分析可知：不同的邻域消费水平对区域消费水平的转移变化所起的作用各不相同。通常而言，在高消费水平的地区邻域范围，区域消费水平往往向上转移的概率将会增加；相反，若以消费水平低的地区为邻，区域消费水平向下转移的概率将增加。例如，1978-1986 年间，在不同邻域消费水平的状态下，区域自身的类型数量及类型转移概率变化存在较大差距，其中，区域同等水平间转移的数量变化随着与其邻域水平高低基本呈负相关；一个低水平区域若以同等水平的地区为邻，其向上转移的概率为0.0373，而当其与中低水平的地区为邻时，其向上转移的概率增加至0.0616。同样，在 1986-2011 年间，一个高水平地区与不同水平的地区相邻，那么其同等水平间转移的数量变化趋势也显着不同，这种变化随着与其邻域水平高低基本呈正相关；一个高水平区域若以同等水平的地区为邻，其向上转移的概率为0.9933，而当其与高水平的区域为邻时，其向上转移的概率增加至1.0000。
综上可知，地区消费水平的演变受到来自邻域区域消费水平状态的影响。通常一个地区若以低消费水平的区域为邻，那么，其消费水平类型向上转移的概率将减小，不利于改变落后的消费水平面貌；若以高消费水平的区域为邻，那么，其消费水平类型向上转移的概率将增加，有助于实现消费水平跨越式地提高。
3.5.2 区域相互作用下的消费水平空间分布格局
区域消费水平类型转移在地理空间上并不是孤立的，往往受到邻域区域消费水平高低的影响。
图4揭示了我国区域消费水平类型转移的空间分布格局，但没有考虑到邻域的转移状况。为此，通过计算空间马尔科夫链，利用ArcGIS对地区及其邻域间的消费水平类型的空间变化进行可视化，为便于分析，本文对区域自身与邻域仅作3层分类，结果如图5所示。
由图5可看出：（1）区域自身消费水平与邻域消费水平同时向下转移的区域主要集中分布于我国西部地区，而区域自身和邻域地区消费水平同时向上变化的空间分布表现并不一致，表明其受区域背景的影响作用较小。（2）区域和相邻地区消费水平同时向上和向下转移的地区数量在1986-2011 年期间明显大于1978-1986 年期间。（3）区域消费水平的空间趋同变化过程分析要考虑来自周围地区的影响。【图略】
4 结论
区域消费水平不仅具有时间演变的历史性，且空间上呈现出显着的地区分异性。随着社会主义市场经济的建立，地理邻近使得区域间人流、物流、信息流、技术流、旅游流等要素的相互交流、相互作用日益频繁，为加快区域一体化进程，促进地区经济实力的提高，政府积极实施推动跨区消费，鼓励地区消费市场不断开放等一系列措施，使得相邻地区间的消费水平相互影响，地区间呈现出消费模仿与消费水平趋同的局面。故此，本文采用马尔科夫链及空间马尔科夫链方法，分别对我国区域1978-1986 年和 1986-2011 年 2 个阶段的人均消费水平的时空格局变化进行分析。结果表明：
（1）传统的马尔科夫链时间变化表明，自1978年以来，无论是低消费水平与高消费水平区域均沿着其最初类型稳步演进，区域消费水平存在着“俱乐部趋同”现象。但是，不同时段内区域消费水平类型数量的变化存在着增加与减少；空间马尔科夫转移概率表明，地区间的消费水平变化受到区域背景影响，其在空间上的趋同过程不独立。总的来说，空间单元若以低消费水平地区为邻，则其消费水平类型向上转移的概率将减小，若与高消费水平区域相邻，那么，其消费水平类型向上转移的概率将提高。
（2）我国区域消费水平呈现出明显的东高西低分异格局，这与区域自身与邻域地区消费水平同时向上、向下转移变动的空间分布格局基本一致。区域自身消费水平与邻域消费水平同时向下转移的地区主要集中分布于我国西部；区域自身和邻域消费水平同时向上转移的地区多分布于东部沿海；中部广大地区多为平稳分布，其相邻区域消费水平状态变化较大。
（3）我国区域消费水平的趋同主要表现为东部地区与西部地区内部空间的趋同，这与我国东西间经济发展水平、商品生产与流通、城乡居民人均收入、城乡居民人均储蓄和社会保障制度等方面存在较大差距有关。因此，今后在制定区域消费政策时不仅要加快经济发展、提高居民人均收入、改善区域消费环境、刺激消费欲望，还应充分考虑周围地区消费背景的影响作用，加强与高消费水平地区的联系，积极发挥其带动刺激效应来改善落后地区的消费水平。
我国区域消费水平在各个研究阶段的变化存在趋同现象，空间马尔科夫链方法可用来解释这种趋同现象。限于研究时间跨度较长，市域层面的数据获取较为困难，本文仅对省际层面的趋同现象进行了初步探究；此外，采用演化树进行我国区域消费水平的时空演变特征分析，将对发现区域消费水平格局新现象新规律具有重要意义。
参考文献：
[ 1 ] 张蕊,田澎.中国城镇居民 10 年消费结构变化实证研究[J].上海理工大学学报,2005,27(2):151-156.
[ 2 ] 蒋云飞,罗守贵.改革开放以来中国城镇居民消费结构变动及区域差异[J].经济地理,2008,28(3):415-418.
[ 3 ] 焦志伦.中国城市消费的空间分布于空间相关关系研究[J].经济地理,2013,33(7):41-46.
[ 4 ] 胡美娟,周年兴,李在军,等.江苏省县域消费水平的空间格局及影响因素分析[J].经济地理,2014,34(6):48-53,74.
[ 5 ] 文启湘,冉净斐.消费结构与产业结构的和谐：和谐性及其测度[J].中国工业经济,2005(8):14-19.
