1998—2018年四川省县域茶叶产量时空演变及影响因素

周佳松，马黛玉，夏建国*，欧定华

(1.四川省国土整治中心，四川成都 610045；2.四川农业大学资源学院，四川成都 611130)

茶叶是集经济效益、社会效益和生态效益于一体的全球性商品，茶叶产业是植茶地区重要的经济型产业[1]。四川省位于中国长江上游，承担着生态保护的重要功能，既是世界茶树的发源地，也是中国重要的茶叶主产区之一，茶叶产业已经成为四川的一张名片。随着东部产业带人力和土地成本的不断上涨，“东茶西移”成为中国茶叶生产的必然趋势[2]，加上新一轮西部大开发战略的实施，川茶产业已迎来新的发展机遇。近年来，川茶产业被四川省委省政府作为第一大优势特色产业培育，在四川乃至全国茶产业发展中具有重要地位。从1998年开始，四川省茶叶产量逐年递增，产茶实力不断增强。到2018年，四川茶叶种植面积逾 37万hm2，茶叶产量破 30万t，实现毛茶产值 242 亿元，综合产值达到 726 亿元,产业综合实力位居全国第二[3]。四川省75%以上的茶园分布在盆周山区和丘陵地区，川茶产业的发展有力促进了地方经济发展和茶农持续增收，已成为乡村振兴、精准扶贫和建设美丽乡村的重要抓手，在促进农村全面建设小康社会中发挥着重要作用[4]。除经济效益外，茶园还具有涵养水源、保持土壤和固碳等生态效益。发展川茶产业不仅有利于发挥四川省在长江经济带发展上的引领作用，还有利于筑牢长江上游生态屏障。因此，分析茶叶产量时空格局演变及影响因素可为政府相关部门制定相关产业政策、发展规划和环境保护提供理论依据。

近年来，国内学者关于茶叶生产的研究主要集中在分析产业发展现状及对策[5-7]、产业扶贫[8-10]和成分研究[11]等方面。如张玉芬等[7]结合温泉镇茶叶生产的现状，从提质和增效两方面探讨了茶产业发展的对策和措施；崔理博[8]基于中国茶产业精准扶贫的现状和“互联网+”的技术背景，提出了精准扶贫的策略。当前，从时空角度对茶叶生产展开研究的较少且主要集中在全球[12]、国家[13]和经济带[12]等大尺度上，难以揭示茶叶生产的空间异质性。从中小尺度上分析茶叶生产的时空演变及影响因素的研究较为薄弱[14-15]，多借助探索性空间数据分析和重心转移模型等方法对地理事物的空间集聚-离散程度及分布格局进行描述，缺乏对影响因素的定量分析。因此，笔者以1998—2018年四川茶叶生产数据为基础，结合空间自相关分析、重心模型和地理探测器模型等方法，对四川省茶叶产量时空格局演变及影响因素进行分析，以期为优化四川茶叶生产布局提供科学依据。

1 材料与方法

1.1 数据来源研究数据包括属性数据和空间数据。其中，属性数据来源于《四川农业统计年鉴》(1998—2018年)，并从中选取1998、2003、2008、2013和2018年作为研究断面，以县域茶叶产量为指标来深入分析四川省茶叶生产的空间演变。空间数据来源于国家基础地理数据库，包括四川省各县级行政区划图和各县级政府的地理坐标，由于研究期内四川省行政区划多次调整，为消除研究期间县域行政区划变动对数据分析的影响，选用2018年行政区划数据。研究所需数据来源统一，可消除统计口径不同带来的误差。

1.2 研究方法

1.2.1茶叶产量空间集聚特征分析。全局空间自相关分析用于衡量某种单元属性在研究区内是否存在集聚特征，常用Moran’sI系数来衡量，计算公式如下[16]：

(1)

局部空间自相关可用来评价地理事物空间自相关的局部异质性，Getis-OrdGi*主要用于识别邻近的区域要素是否具有相关性，可以检验高值区和低值区在空间上的集聚，即热点区和冷点区的空间分布格局，计算公式如下[17]：

(2)

式中，Gi*表示局部空间关联性指标；其余参数与公式(1)相同。

1.2.2茶叶产量重心演变分析。重心模型是研究区域发展过程中空间要素变动的重要分析工具，茶叶产量重心的移动可以用于分析茶叶生产的空间动态发展轨迹。近年来，学者们多运用重心统计模型对经济重心[18]、人口重心[19]和粮食生产重心[20]进行分析。该研究依据重心理论模型，构建区域茶叶产量重心模型。

假设研究区由m个次级单元构成，第j个次级单元包含某种属性G的重心坐标为(Xj,Yj),表示该次级单元某属性重心的地理坐标，Gaj为第j个次级单元第a年某种属性的值,则该研究区某属性的重心坐标(Xa,Ya)和重心移动距离计算公式如下[19]：

(3)

(4)

式中,Da-b表示第a、b年份间某种属性重心移动的距离，(Xa,Ya)、(Xb,Yb)分别表示第a年和b年的区域重心所在空间的地理坐标；此处C取值111.11，是将地理坐标转为平面距离的常数。

1.2.3影响因子分析。当变量具有空间相关性时，普通最小二乘法回归模型不再适用于此类数据的分析。地理探测器是探测空间分异性及其驱动力的工具，该方法的一大优点是对变量共线性免疫,这主要是因为如果不对q值进行显著性检验，q值仍具有明确的物理含义，没有线性假设[21]。地理探测器作为有力的驱动力分析工具被广泛应用于公共健康[22]、粮食生产[23-24]和生态环境[25-26]等领域。该研究选用因子探测器识别各影响因子(X)对茶叶产量(Y)空间分异的解释力大小，计算公式如下：

(5)

式中，h=1,…,L为茶叶产量Y或影响因子X的分层；Nh和N分别为层h和全区的单元数；σh2和σ2分别是层h和全区的Y值的方差。q为影响因子对茶叶产量的影响程度，值域为[0,1]，q值越大表示影响因子X对茶叶产量Y的解释力越强，反之则越弱。

茶叶生产受多因素共同影响，考虑到数据的可获得性和科学性，该研究共选取农村人口(X1)、农林牧渔总产值(X2)、人均国民生产总值(X3)、茶园面积(X4)、茶叶单产(X5)和路网密度(X6)等6个指标，通过自然断点法将所选因子分为9个等级，使用公式(5)分析不同影响因子(X)对茶叶产量(Y)空间分布差异的影响强度。

2 结果与分析

2.1 茶叶生产时空演变分析由表1和图1可知：1998—2018年，产茶县域数量变化不大，但县域间产量变化较大。茶叶产量低于1 000 t的县域数量由80个减少至57个，这类县域的茶叶产量占全省茶叶总产量的比例(以下简称茶叶产量占比)大幅下降，低产县逐年减少，主要分布在川西高原地区和部分川中丘陵区。1998—2018年，茶叶产量在≥1 000～3 000 t和≥3 000～5 000 t的县域数量变化不大，茶叶产量占比则明显降低，这2类茶叶生产水平的县域空间分布较为均匀。茶叶产量在≥5 000～10 000 t和≥10 000～20 000 t的县域数量从0分别增长到12和6个，茶叶产量占比均呈上升趋势，集中分布在川西和川南地区。研究期内，茶叶产量达2万t及以上的县域极少，1998—2018年仅有2个县域达到此生产水平，茶叶产量占比呈小幅增长。

总体而言，1998—2018年，茶叶产量≥1 000 t的县域数量由17个增加至38个，其茶叶产量占全省总产量的比例由66.74% 增加至95.98%，表明全省茶叶产量显著增加；茶叶产量≥5 000 t的高产县域数量由0增加至20个，茶叶产量占全省总产量的比例由0增加到82.07%，主要集中分布在位于四川盆地西南缘和南缘的川西产业带和川南产业带，成为四川省茶叶生产的优势区域；川东北产业带高产县域分布分散，成为茶叶生产相对劣势区。

2.2 空间集聚特征分析利用公式(1)计算出1998—2018年四川省县域茶叶产量5个研究断面Global Moran’sI值依次为0.152 9、0.182 8、0.100 2、0.093 7和0.343 5，全部为正值且通过0.01水平显著性检验，这表明茶叶生产在空间分布上并非随机，而是呈现出集聚特征，即茶叶产量较高(低)的县域其邻近县域茶叶产量也比较高(低)，空间上存在高高集聚、低低集聚的分布特征。1998—2013年四川省县域茶叶生产空间关联度逐渐降低，2013—2018年又显著上升，远高于1998年，表明近5年茶叶产量相似的县域空间关联度显著增强。

表1 1998—2018年四川省县域茶叶生产分布Table 1 Distribution of tea output in counties of Sichuan Province from 1998 to 2018

图1 1998—2018年四川省县域茶叶产量时空格局Fig.1 Spatiotemporal pattern of county tea output in Sichuan Province from 1998 to 2018

全局自相关指数解释了茶叶产量在空间上具有集聚特性，但不能解释其空间分布规律。因此，根据公式(2)计算四川省县域茶叶产量Getis-Ord ，将产茶县域划分为3级热点、3级冷点和不显著7个类型，其中“一级、二级、三级”分别表示通过置信度为99%、95%、90%的统计显著性检验。由图2可知:茶叶产量冷点分布区域变化不大，但年际差异大；热点区域分布变化较大且年际变化显著，1998—2003年，茶叶产量冷点主要分布在四川盆地东北部且不断缩减，2008—2018年冷点分布急剧减少。1998—2003年，茶叶产量热点在成都市、雅安市、眉山市、乐山市和宜宾市均有分布，热点聚类显著性增强，且呈现出向西集中的态势；盆地南部热点聚类显著性减弱，热点范围缩小。2008—2013年，热点区域聚类显著性进一步增强，热点范围扩大但地理位置更为集中，呈球状分布在四川盆地西南部，主要包括川西产业带中的雅安市名山区、雨城区、天全县、芦山县、荥经县和宝兴县，成都市邛崃市和蒲江县，眉山市东坡区、丹棱县、洪雅县和彭山区，乐山市西北部夹江县、沙湾区和峨眉山市等15个县域。2018年，盆地西南部热点区域向西退缩，主要集中于雅安市东北部芦山县、中部雨城区和名山区，成都市大邑县、邛崃市和蒲江县，眉山市东坡区、丹棱县和洪雅县，以及乐山市和峨眉山市等10个县域。

图2 1998—2018年四川省县域茶叶产量空间集聚热点演化Fig.2 Hotspots in spatial concentration of tea output at county level in Sichuan Province from 1998 to 2018

2.3 茶叶产量重心变化分析由表2和图3可知，茶叶产量重心位于四川盆地西南缘，在川西和川南产业带范围内，地理位置为103°54′16″～104°18′56″E、29°43′43″～29°49′14″N。茶叶产量重心的经向变化更为明显，主要表现为由西北向西向东南再向东迁移，纬向变化较小，说明四川省茶叶产量在东西方向上变化较大。1998—2008年，茶叶产量重心由东南向西南摆动，平均移动速度为4.64 km/a，表明在该时段内，西南部茶叶产量增加明显，茶叶生产水平较高；2008—2018年，茶叶产量重心主要自西南向东摆动，平均移动速度为2.77 km/a，年际移动速度放缓，表明该时期内四川盆地西南部偏东位置茶叶产量稳步增加，川西和川南产业带之间的产量差异逐渐缩小。总的来说，1998—2018年，四川省茶叶产量重心主要呈现出由快到慢的移动特征，移动方向为“东—西—东”，空间变动范围不大，平均移动速度为3.70 km/a。

表2 1998—2018年四川省茶叶产量重心变化Table 2 Changes in the gravity center of tea output in Sichuan Province from 1998 to 2018

2.4 茶叶生产空间格局变化影响因素分析由表3可知，不同时期茶叶产量空间分异影响因素的解释力差异较大。研究时段内农村人口和茶园面积作用强度随时间变化先增后减，总体解释力较大，变化明显；农林牧渔总产值、茶叶单产和路网密度作用强度先减后增，解释力呈波动变化；表征人民生活水平的人均国民生总值解释力强度持续减弱。1998—2003年，表征茶叶生产规模的茶园面积和表征茶叶生产水平的茶叶单产成为影响茶叶产量空间分异的主要因素，其余因素按作用强度分别为农村人口、路网密度、人均国民生产总值和农林牧渔总产值；2008—2013年，表征劳动力的农村人口因子解释强度逐渐增加，成为仅次于茶园面积的第二大影响茶叶产量空间分异的主要因素，茶园面积在2008年对茶叶产量空间分异的解释强度达90.3%，远超其他因子；2018年，茶叶产量空间分异的驱动因素再次变化，表征地区交通发展水平的路网密度取代农村人口成为仅次于茶园面积的茶叶产量空间分异影响因子。

总体而言，近20年来茶园面积是决定茶叶产量空间分异的主要因素，但与1998—2013年相比，2018年茶园面积对茶叶产量的影响明显降低，即当茶园面积达到较高水平后，单纯依靠茶园面积扩张来提高茶叶产量的效率降低。茶叶单产、农村人口及路网密度各自在不同的时间断面成为茶叶产量空间分异的第二大主要影响因素，但这些因素解释力波动变化较大，表明未来茶叶产量稳产应注重增加科技投入，吸纳农村人才回流并加强专业技能人员培训投入，完善乡村道路布设，以保证茶叶单产的稳步提高、劳动力数量增加和素质提高以及交通水平的提高。

3 讨论

四川省茶叶产量空间集聚显著，形成明显的川西产业带、川南产业带和川东北产业带三大茶叶生产空间格局，与林正雨等[15]研究结果一致，且与当前四川省茶产业发展状况基本一致，表明四川省茶叶生产空间格局发展趋于稳定，当前产业布局较为合理。具体来看，三大茶叶生产空间的产量差异显著，其中川西产业带为热点集聚核心区，拥有全国最大的无性系繁育基地。川南产业带为次热点区，近年来与川西产业带的茶叶产量差距逐渐缩小，二者均为茶叶生产优势区。川东北产业带茶叶产量较低且为冷点集聚区，是四川省茶叶生产的相对弱势区，形成了“西南高东北低”的生产格局。茶园面积是造成茶叶生产“西南高东北低”的决定因素，这是因为茶园面积决定茶叶生产规模[13]，进而影响茶叶产量。茶园面积的作用机制与四川省生态建设工程密切相关。

图3 1998—2018年四川省茶叶产量重心演变轨迹Fig.3 The evolution of the gravity center of tea output in Sichuan Province from 1998 to 2018

表3 不同年份各影响因素解释力变化

四川省自20世纪30年代以来，高强度垦殖坡地，导致自身水土涵养能力下滑，其生态变动对下游影响巨大。1998年长江洪水汹涌而来，四川成了长江干流的主要泥沙输入地。1999年，四川省开始启动退耕还林还草试点工作，成效卓然。20年来，通过茶叶种植及其他经济作物种植逐渐走出了一条“山上长叶子，农户得票子”的农村发展新路径[27]。川西和川南产业地处盆周山地，气温适宜、云雾多湿度大，漫射光多，昼夜温差较大，有利于茶叶含氮化合物的合成和累积，自然资源禀赋优于川东北地区，是茶叶生产的适宜气候区[28]，不同区域的茶园面积扩张差异较大导致区域间茶叶产量差距也不断增大。茶园面积是影响茶叶产量的主导因素，这与旷成华等[12]的研究结果相似，但与肖智等[2,13]的研究结果不同,其研究显示茶园标准化管理水平、茶园机械化水平、茶数品类及新技术推广等因素对茶叶产量有重要影响。同时，有研究表明，一味依靠茶园扩张换取产量增长容易导致生态环境恶化[29]，因此，未来应考虑将影响茶叶产量的第一驱动力转变为茶园机械化水平和茶园标准化建设水平等，增强茶叶产业的可持续发展。当前，四川茶园机械化水平低于10%，仍属于劳动密集型产业[30]，在茶叶生产过程中主要靠大量的劳动力投入提高茶叶产量，因此茶叶产量空间分布与农村人口空间分布高度契合，而近年来城镇化发展使得越来越多农村人口转向城市，茶叶生产的劳动力投入减少，使得人口因素影响减弱。

该研究以县域为单元进行分析，采用空间分析和地理探测器等方法，揭示了1998—2018年四川省茶叶产量的时空演变特征及驱动机制，弥补了前人定性分析的部分缺陷。但受数据可获得性的限制，该研究在指标选取上未考虑农业机械总动力、气温和降水等因素对茶叶产量空间分布的影响，未来应深入研究上述因素对茶叶生产的影响，并进一步探讨各影响因素在不同研究区域内的作用强度。

4 结论

(1)1998—2018年，四川省茶叶产量呈增长态势，空间分异明显，存在显著的空间集聚特征。川西和川南产业带茶叶产量热点广布，是茶叶生产的优势区域；川东北产业带为茶叶产量冷点集聚区，成为茶叶生产的相对劣势区。

(2)1998—2018年四川省茶叶生产重心呈由快到慢的移动特征，总体移动范围变化较小，以3.70 km/年的速度在四川盆地西南缘内移动并表现出东移的趋势，茶叶生产“西南高东北低”的态势更加显化。

(3)县域茶叶产量空间分异的主要影响因素包括农村人口、茶园面积、茶叶单产和路网密度，各因素在不同时间段内对茶叶产量的空间分异作用强度和作用机制不同，茶园面积是茶叶产量空间分异的决定因素，考虑到茶园扩张的危害，未来应转化茶叶产量空间分异的第一驱动力，促进茶叶产业的可持续发展。