岷江上游农民收入时空异质性的地理因素及其演化格局

　　异质性( Heterogeneity) 是指生态学过程的不均匀性及其复杂性,广泛应用于生物和医学领域,就社会科学领域来说,异质性普遍存在于各种社会现象中。其中,人均收入的时空差异则是区域经济学领域异质性的一种特殊表现形式。区域收入异质性是当前我国区域发展面临的一个主要问题,国内外学者就不同空间尺度的区域经济异质性问题做了大量研究。从全国尺度看主要有魏后凯、杨开忠、胡鞍钢等以人均收入为指标研究区域间收入差异;Tsui KaiYuen 采用人均国民收入指标研究 1952—1985 年间中国区域发展差异情况; Lyons 利用人均净产值分析 1952—1987 年我国区域差异问题;李小建等以县域为基本单元,在全国尺度下对比县际经济差异。Ravi Kanbur 研究 1983—1995 年间中国农村与城市之间的经济差异演化。从区域尺度看,吴殿廷,李二玲和覃成林等对南北经济发展差异做了大量研究; 徐建华,鲁凤,陈国阶、魏后凯等人用 GDP 作为测度指标研究我国三大地带的差异演化问题并分析其成因; 此外还有其他学者从省域尺度研究经济差异问题。

　　综述当前研究,发现大部分学者集中在较大尺度或经济发达地区经济异质性问题。对于我国农村地区、西部山区以及少数民族集中区的人均收入异质性研究仍有待深入。回溯 20 世纪以来,中国农村、农业和农民( 三农) 问题在区域发展中一直占据十分重要的地位,学术界对相关领域的研究主要体现在农民收入差异测度演化、原因分析及其宏观调控,农村经济差异与发展对策以及农村公共服务水平区域差异,等等。研究区域主要集中在全国或典型农业省市,针对微观地域单元或江河上游地区农村收入时空异质性问题仍缺乏足量实证支撑。岷江上游地区位于我国西部省际边缘区,是当前我国贫困人口比较聚居的地区,生计状况极为脆弱,农户收入来源单一并不稳定,各乡镇农民收入在空间上存在明显的梯度效应。

　　因此,在此背景下,论文以岷江上游地区为例,以乡镇为单元,分析区域内农民收入时空异质性的地理因素及其演化格局,旨在为岷江上游地区经济发展决策提供参考。

　　1、 研究方法

　　1. 1 研究区概况

　　岷江是长江重要支流之一,发源于青藏高原东端岷山南坡。岷江上游指岷江流域都江堰以上区域( 31°26' ~ 33°16'N,102°59' ~ 104°14'E ) ,其范围与四川省阿坝藏族羌族自治州的汶川县、理县、茂县、黑水县和松潘5 县的行政辖区基本重合,全流域面积约为 23 078 km2。岷江上游地区属高山峡谷区( 图 1) ,地质构造复杂,断裂发育,支流众多,径流以降雨补给为主,冬干夏湿,冬冷夏凉,干湿季节分明,具有高原型季风气候特征; 流域内森林资源丰富,生物种类繁多,是一个具有高度生物多样性的生态敏感地区。岷江上游聚居了以羌族、藏族等为主的大量少数民族,是经济相对滞后的少数民族贫困山区,以旱作农业、林业和牧业等第一产业为主导产业; 2010 年,汶川县、理县、茂县、黑水县和松潘 5县共计地区生产总值 74. 69 亿元,三次产业增加值结构为 11. 55∶ 61. 54∶ 26. 91,农民人均纯收入约3657 元 / 人,共计人口约 39. 15 万,其中农业人口29. 13 万人。

　　

　　1. 2 研究方法

　　有关评估区域异质性的统计方法很多,为避免单一测定方法出现误差以及考虑到人口因素的影响,采用加权变异系数和泰尔指数来评价“十一五”期间岷江上游农民收入差异的变化情况。

　　1. 2. 1 加权变异系数

　　变异系数指样本标准差除以平均数,是一个反映总体差异的相对指标。考虑到了人口规模的影响,采用加权变异系数来衡量差异的变动趋势。加权变异系数计算公式为

　　式中 N 为乡镇数,Xi和 X 分别表示 i 乡镇和全部乡镇的平均收入,Pi和 P 表示 i 乡镇和总人口数,Pi/ P 用以反映 i 镇人口规模的影响,CV 越大表明某属性值的差异越大。

　　1. 2. 2 泰尔指数

　　泰尔( Theil) 指数又称泰尔熵,采用农村经济纯收入代替传统研究中所用的 GDP 指标,即可进行农村经济差异研究。其中泰尔 T 的计算公式为

　　式中 N 为乡镇数,Yi为 i 乡镇农村经济纯收入占岷江上游总体的份额,Pi为 i 乡镇的人口数占岷江上游总体的份额。T 值越大,表示各乡镇间农村经济收入发展水平差异越大; 反之越小。

　　通过泰尔指数的一阶分解,可以将总体差异分解为组间和组内差异,计算公式

　　式中 Tn,T*,T 分别表示总体差异、组间差异、组内差异,i 表示县数( i = 1,2,…n) ,j 为组内子区域( 乡、镇) ,Yi表示 i 县的收入占总体的比重,Pi表示i 县人口占总人口比重,Yij表示 j 乡镇在 i 县收入比重,Pij表示 j 乡镇人口在 i 县人口比重。

　　1. 2. 3 探索式空间数据分析方法( ESDA)

　　引入 Moran’sⅠ和 Local Moran’s Index 来测度全局和局部的空间关联特征。全局空间自相关分析采用 Moran’s I 来分析空间关联度和空间差异程度:

　　

　　n 表示乡镇数目,xi乡镇 i 农民人均纯收入,x为平均值,Wij为空间权重矩阵。局部空间自相关分析用来识别不同空间位置上可能存在的不同空间相关模式,一般采用统计量Local Moran’s I:

　　

　　Moran’s I 统计量取值介于 - 1 和 1 之间,其绝对值越趋近于 1,表示单元的空间相关性越显著。

　　Moran 散点图描述变量与其空间滞后量之间的相关关系,划分的四个象限,分别对应四种空间差异类型: 第一象限( H - H) Ii> 0,说明相邻空间的值都高于总体平均水平; 第二象限( L - H) Ii< 0,i 区域值低于总体平均水平,而相邻区域却高于平均水平; 第三象限( L - L) Ii> 0,相邻空间的值都低于平均水平; 第四象限( H - L) Ii< 0,i 区域值高于总体平均水平,但相邻区域值却低于总体平均水平。

　　1. 3 数据来源和时间尺度

　　选择岷江上游 5 县( 汶川、茂县、理县、黑水、松潘) 属于岷江流域的乡镇作为研究对象,采用农村经济纯收入、乡村人口及农民人均纯收入,研究2005—2010 年 6 a 的农民收入差异演化及其空间分布格局。数据来源于汶川县、茂县、理县、黑水县、松潘县统计年鉴。

　　2、 研究结果

　　2. 1 农民收入异质性演变特征

　　表 1 和图 2 分别给出了岷江上游总体,以及县域间和县域内部的农民人均纯收入 Theil 系数和变异系数,反映了以上三者之间差异的构成和变化趋势。

　　2. 1. 1 农民收入差异总体缩小

　　分析发现,两种方法所得到的研究结果相似: 从2005 年和 2010 年两个时间截面数据看,除松潘县农民收入差异呈小幅增加外,其余 4 县以及岷江上游总体的农民收入差异呈现缩小的趋势; 从连续的时间截面看( 2005—2010 年) ,各年份间的差距表现出增减相间的波动状态; 松潘县内各乡镇农民收入差距总体增大,但在 2006—2007 年间曾出现了收敛,此后又发散; 黑水县在 2005—2010 年间农民收入差异呈现稳定收敛趋势; 理县在 2008 年各乡镇收入差距达到最大,此后又逐渐缩小,呈现出一个倒“V”形的变化趋势; 茂县内各乡镇农民收入差异总体呈现收敛。汶川县在 2008—2009 年间出现反弹,农民收入差距急剧扩大,2009 年后又急剧缩小。

　　

　　2. 1. 2 县域内部差异是构成总体差异的主要因素

　　县域内部差异对农民收入总体差异的贡献程度稍大,是构成岷江上游总体差异的主要因素。其中,汶川和理县的差异贡献率始终处于高值区而茂县和黑水的贡献率始终处于低值区; 即汶川和理县内部各乡镇的农民收入差距在岷江上游的总体差异中有重要影响。松潘农民收入差异对总体贡献率稳步上升,贡献率从 2005 年的 6. 03% 上升到 2010 年的13. 02% ,对总体差距的贡献率仅次于理县( 17. 5%) ,而茂县贡献率呈现稳步下降,贡献率从2005 年的 7. 79% 下降到 2010 年的 1. 73% 。农民收入差异在“十一五”期间存在明显异常点,县域之间和汶川县内部差异在 2009 年呈现明显突变,县域间差异贡献率从 52. 2%降到 28. 9%,汶川县内部差异贡献率从 15%上升到 41. 3%,导致其成为构成岷江收入差异的最主要因素。

　　2. 2 农民收入异质性空间格局

　　通过对岷江上游各乡镇农民人均纯收入进行全局空间自相关分析,2005—2010 年农民人均纯收入的 Moran’sⅠ指数分别为: 0. 637 2、0. 608 9、0. 6497、0. 687 6、0. 387 0、0. 588 8。表明农民人均纯收入在不同年份呈现显著的空间正相关特性,收入水平相似的乡镇在空间上呈明显集聚趋势。通过 Mo-ran’sⅠ散点图和 LISA 集聚图分析发现,2005 年落入第三象限的乡镇数最多( 低 - 低集聚) ,为 41 个乡镇,占全部乡镇总数的 46%,落在第一象限的乡镇次之,有 35 个乡镇,占总数的 42%。到 2010 年,落在第一象限的乡镇变成了 47 个( 高 - 高集聚) ,增加 12 个乡镇,占全部总数的一半以上,达到52% ,第三象限的乡镇数减少了 7 个,为 28 个乡镇,占总数的比例减少到 31%。说明 2005—2010 年间,农村人均收入高值区的空间相关度有所提高。

　　从空间分布看,在 1%的显著性水平下,农民人均收入的高 - 高和低 - 低集聚空间特征明显,低值区主要集中在岷江中部的黑水和茂县,而高值集中在松潘东北部,汶川南部和理县西部,但随时间演变其集聚趋势又明显变化,2009 年岷江南部的高值集聚区消失,说明该年农村返贫现象明显( 图 3) 。

　　3、 农民收入差异的地理因素

　　3. 1 海拔决定收入格局

　　自然地理因素是中国农村贫困发生的第一因子。已有研究表明中国农村贫困是由特殊的区域环境、生产条件以及人口素质等一系列带有明显地域性特征的约束条件导致的“点上”贫困。一般来说,随着海拔的升高,有效积温降低,农民收入会随耕地单产水平降低而降低。

　

　　分析发现,岷江上游农民收入与海拔拟合的线性关系不明显。但如果以收入均值为基准,发现随着海拔的升高,高收入乡镇的收入水平有逐渐降低的趋势; 当海拔大于 2 500 m 以后,大部分乡镇处于低水平区域; 当海拔大于 3 500 m 后,大部分乡镇收入水平又整体位于高值区域( 图4) ,2 500 ~3 500 m的海拔带构成了岷江上游农民收入低值区。因此,岷江上游贫困空间集聚的原因与其所在的海拔有明显关系,低海拔乡镇的农村收入以种植业为主,加上其便利的交通条件,使该区域成为成都市的蔬菜供应基地。高海拔区域乡镇大多位于青藏高原面上,坡度平缓,具有天然的牧草优势,农村收入以畜牧业为主,加上独特的草原风光和旅游资源,带动当地旅游业发展,对农民增收具有明显促进作用。处于的农牧交错地带的乡镇,所在区域坡度较陡,且常伴山地灾害发生,生产力普遍偏低加上收入结构单一,形成中海拔区域的低收入集聚区。

　　3. 2 交通条件反馈资源优势体现程度

　　岷江上游一直就是我国北方游牧民族向南扩展的通道,其沿江地带成为北方和南方农牧业发展交汇地带,也成为商贸业发展的走廊。随着社会经济的发展,工业化和城镇化的加快,岷江上游的对外交通虽然有改善,但内部各乡镇的可达性差异较大,交通不便,接受外来投资困难的状况突出,本地区的特色资源难以与区外进行交换,导致各乡镇农民收入呈现差异。

　　对比 2010 年岷江上游农民收入和交通分布发现,横穿岷江上游的两条国道与农民收入分布有较好的空间耦合关系( 图 5) 。收入较高的乡镇主要沿213 和 318 国道分布。其中汶川、茂县具有与成都较近的区位优势,加上较为便利的交通条件,两县成为成都市的蔬菜供应基地。位于上游松潘县的牟尼乡和川主寺镇等乡镇,虽然海拔较高,但由于黄龙机场的建设减少了成都到松潘的通达时间,使岷江上游高海拔地区的旅游资源得到充分利用,促使该区域的农民收入结构多元,成为收入水平的高值集聚区。而黑水县大部分乡镇居民多处在高半山的农地边缘,境内主要交通公路沿河谷走向,除沿线分布的聚落能够享受便捷的交通和良好的基础设施外,大部分乡镇远离公路,可进入性差,也因此构成了农民低收入乡镇的一个集聚区域。

　　3. 3 旅游资源优势促进农民增收

　　自然资源是人口增长和经济发展的物质基础,它的数量及其分布对人口与经济作用不可低估。岷江上游生物资源丰富,珍稀动植物品种多。林地面积占全区土地总面积的 49. 85%,有卧龙自然保护区、宝顶沟自然保护区 2 处国家级自然保护区等,有大熊猫栖息地 1 处世界自然遗产地,以及四川省草坡自然保护区和阿坝州七盘沟保护区( 表 2) 。

　　对比岷江上游各乡镇农民收入的空间分布与风景区分布图( 图 5) ,发现在全国和四川省的主体功能分区中被明确列为禁止开发区的乡镇,如汶川卧龙镇、银杏乡、三江乡以及米亚罗镇等,农民人均收入都位于岷江上游的高值区,该区域独特的自然风光吸引外来观光者,带动当地经济发展,以至于农民收入甚至超过了汶川威州镇,映秀镇,理县杂谷脑镇等一些工业重镇。这说明充分发挥地区不可替代性的旅游资源优势加上便利的交通条件,可以促进资源优势转化为经济优势,成为农民增收的有效途径。

　　3. 4 汶川地震对收入结构影响显著

　　岷江上游位于青藏高原东缘,川西高原东北边缘与龙门山过渡地带,地质构造复杂,隆升运动强烈,历史上发生过多次地震,加上其高山峡谷地貌,泥石流、滑坡等山地灾害发生频繁。自然灾害的发生,在不同程度上造成农作物的减产和绝收,这势必影响到农业生产的经济效益,进一步影响区域农户再生产投入。岷江上游乡镇农民主要收入来源是家庭经营性收入,其他收入比重不高,尤其是在个别乡镇,家庭经营性收入仅能解决温饱问题,一旦遭受自然灾害,农户收入必然受到极大的影响。2005—2010 年的农民收入差异 Theil 指数和变异系数分析发现: 2007—2009 年间岷江上游县域之间的差异出现异常,尤其在 2008 和 2009 两个时间截面上表现得尤为突出,2008 年各县的农民收入差距总体减小,分析各县 2005—2010 年农民收入结构,可以发现汶川大地震导致岷江上游农民普遍返贫; 农民经营性收入大幅降低甚至出现负值,收入主要来源主要依靠转移性收入( 图 6) 。但由于灾后重建的投入力度差异决定了农户转移性收入的差异,最后直接导致 2009 年的汶川县各乡镇农民纯收入差异急剧增大。通过 Local Moran’I 散点图与 LISA 集聚图可以发现,2009 年岷江上游的空间自相关表现不显著,( L - H) 的数量增加,高值集聚区( H - H) 减少( 见图 3) ,这与 2008 年的汶川大地震在时间上有明显耦合关系,自然灾害造成的农村返贫现象明显。

　　4、 结论与讨论

　　岷江上游资源环境承载力较弱,属于全国生态脆弱地区,也是我国贫困人口集中区域,面临生态环境保护和经济发展的尖锐矛盾。论文利用 Theil 指数、变异系数以及探索式空间数据分析方法,对岷江上游五县农民收入差异进行研究,得到以下结论:

　　1.“十一五”年间,地区农民人均收入差距呈现总体缩小趋势; 从县域内部看,除松潘县农民收入差距呈现小幅度增大外,其余 4 县的农民收入差距在波动中呈现逐渐缩小的趋势; 从构成差距的贡献值看,县域内部差异对总体差异的贡献较大,说明农民收入差异主要原因是由于县域内部乡镇间差异过大所致。

　　2. 岷江上游农民收入差异变化趋势在 2009 年出异常,通过 Theil 指数分析,汶川县的差异贡献率从 2008 年 15%突变到 2009 年的 41. 3%,是构成农民收入差异的主要因素; 而农民收入的 Moran’I 指数在该年份也出现异常突变,在 LISA 集聚图上也呈现出空间格局的显著变化,这与 2008 年汶川地震有一定的时空耦合关系,反映自然灾害对农民收入及其空间格局有着重大影响。

　　3. 从农民人均纯收入的空间分布看,呈现显著的空间正相关,收入水平相似性的乡镇在空间上呈现高 - 高和低 - 低集聚状态; 高值集中在南北两侧,低值主要集聚在中部偏西,且不同的时间截面又表现出明显的差异; 低值集聚区大部分乡镇位于高半山的农牧接壤区,远离国道线,自然环境条件恶劣,高值区大部分集中在交通条件较好的低海拔乡镇或高原乡镇。

　　4. 卧龙镇、耿达乡、三江乡以及米亚罗镇等乡镇位于岷江上游的生态敏感区,这些乡镇位于国家自然保护区,世界自然遗产地区域,属于禁止大规模高强度工业化和城镇化开发区域,但都成为农民收入高水平乡镇。便利的交通条件和丰富的自然资源优势,是该类乡镇农民具有较高收入的主要原因之一。

　　说明生态保护地区争取国家财政转移支付的同时,有效保护兼合理利用区内资源优势,是该类地区提高农民收入的有效途径之一,为其他类似生态区农民增收提供借鉴作用。

　　论文以岷江上游为例,以乡镇数据为基本单元,探讨县域尺度背景下农民收入差异时空演变及其成因,限于数据的获取问题,未能就乡镇尺度下农民收入结构及其差异演化进行深入分析,而且缺乏长时间尺度的收入差异空间分析,对单一指标的研究仍显不足。现阶段我国农村问题的研究,已经从单纯的农民收入问题逐渐深入到农民民生质量评估与调控,就农民民生质量的差异问题来看,需要从农村经济收入、农村医疗和社会保障、农村公共服务设施拥有量、农村人均营养水平等方面差异进行综合深入分析。因此,下一步研究中,拟基于多指标综合的微观尺度,进行更长时间尺度的农民民生质量差异研究,探讨政策变化对农民民生质量差异的形成和影响机制,为区域差异评估调控和空间管制措施提供决策参考。