重庆地势由南北向长江河谷逐级降低，西北部和中部以丘陵、低山为主，东南部靠大巴山和武陵山2座大山脉，涉及秦巴山区和武陵山区2个集中连片特困地区，扶贫开发任务较重，重庆下辖的38个区县中涉及扶贫开发任务的有33个区县，2014年建档立卡贫困人口达165.9万人，贫困发生率为7.1%，截至2018年底151.99万人实现稳定脱贫，贫困发生率降至0.70%.

空间自相关分析是理论地理学的基本方法之一，包括全局空间自相关分析和局部空间自相关分析.使用全局空间自相关分析可以判断重庆市县域农村贫困化空间分布是聚类模式、均匀模式还是随机模式.根据区县位置及其贫困发生率属性使用Global Moran's I统计量测度重庆市农村贫困化空间自相关性. Moran's I值计算公式如下：

式中：n为样本量，即空间位置的个数；x_i，x_j是空间位置i和j的观察值，为所有空间单元的观测值的平均值；w_ij表示空间位置i和j的临近关系，当i和j为临近的空间位置时，w_ij=1，反之则为0.全局Moran's I系数的取值范围为(-1，1).

局部空间自相关指标LISA用于度量某空间单元与周围邻近单元的空间差异程度.用局部Moran's I系数来计算，其公式为：

式中：I_i为空间单元i的局部Moran's I值，其他变量的含义与式(1)相同.

传统的线性回归模型只是对参数进行“平均”或“全局”估计，如果自变量为空间数据，且自变量间存在空间自相关性，就无法满足传统回归模型(OLS模型)残差项独立的假设，那么用最小二乘法进行参数估计将不再适用.引入地理加权回归(GWR)模型对不同区域的影响进行估计，能够反映参数在不同空间的空间非平稳性，使变量间的关系可以随空间位置的变化而变化，其结果更符合客观实际，因此研究引入了GWR分析，在全局回归模型的基础上进行局部的参数估计，模型结构如下：

式中：y_i为样本i的因变量；β₀为截距；(u_i，v_i)是第i个样本空间单元的地理中心坐标；β_k(u_i，v_i)是连续函数β_k(u，v)在i样本空间单元的值.

本研究采用Arc GIS 10.2中的GWR工具来实现GWR模型的构建.相比OLS模型，GWR模型具有以下优点：①在处理空间数据时，模型的参数估计和统计检验较OLS模型更加显著，并且具有更小的残差；②每个样本空间单元对应一个系数值，模型结果更能反映局部情况，能够还原OLS模型所忽略的变量间关系的局部特性；③能够通过GIS对模型的参数估计进行空间表达，探索空间变异特征和空间规律.

研究数据包括社会经济数据、土地利用数据、农村贫困数据和数字高程模型，其中社会经济原始统计数据主要采集于《重庆市统计年鉴》(2014-2018年)、重庆市各区县同期统计年鉴以及EPS数据平台专业数据库，数据来源于重庆市统计局.土地利用数据数据来源于重庆市国土资源和房屋管理局.农村贫困数据主要包含贫困发生率、贫困人口信息、贫困村相关信息数据，数据来源于全国扶贫开发信息系统中重庆市相关数据.数字高程模型(DEM)，数据精度为25 m×25 m，数据来自NASA网站下载后运用ArcGIS裁切工具(Extact by Mask)按照研究区重庆市范围进行裁切，以此来做地形条件分析.

重庆市共有个14个国家级贫困区县，其中城口县、云阳县、巫山县、巫溪县、奉节县5个国家级贫困县分布于秦巴山区连片特困地区；丰都县、石柱土家族自治县、秀山土家族苗族自治县、酉阳土家族苗族自治县、彭水苗族土家族自治县、黔江区、武隆区7个国家级贫困区县分布于武陵山区连片特困地区.重庆市集中连片特困地区区县总面积达到40 717.95 km²，占所有贫困县面积的84.59%.

从扶贫开发主管部门提供的重庆市2014-2018年各县贫困人口的统计数据来看，全市农村贫困发生率由2014年的7.1%下降到2018年底的0.70%，降幅接近90%，另外，全市33个具有扶贫开发任务的区县自脱贫攻坚以来贫困发生率总体呈现下降趋势.具体来看，2014年全市贫困发生率为6.51%，其中贫困发生率高于10%的区县涉及城口、巫山、酉阳等10个区县，2018年底，重庆市各区县贫困发生率均低于10%，贫困发生率最高的是巫溪县，为3.25%.

为使区县之间农村贫困化具有可比性，本研究借鉴已有文献和现行政策，并结合数据样本特征，将区域农村贫困人口与区域当年农业人口的比值作为衡量区域农村贫困化水平的指标，即农村贫困发生率(PI).根据重庆市实际情况将农村贫困发生率(PI)分为3类：0＜PI≤3%为第一类，农村贫困化程度较小，属于低度贫困区；3%＜PI＜10%为第二类，农村贫困化程度较大，属于中度贫困区；PI≥10%，农村贫困化程度大，属于深度贫困地区.

受自然资源禀赋、社会经济发展水平以及其他因素的影响，重庆市下辖区县贫困发生率空间分布地域性差异显著.本研究在Arc GIS 10.2软件平台支撑下对重庆市2014年、2015年、2016年、2017年和2018年农村贫困发生率的差异进行了空间表征(图 1)，由图 1可知：①从总体上来看，研究时段内重庆市县域贫困发生率热点区域整体空间格局稳定，高值区主要分布于巫溪、城口、酉阳和彭水等渝东北和渝东南区域的相关区县，而主城区和渝西地区的区县贫困发生率相对较低(图 1)，呈现显著的空间差异性；②从时间序列上来看，各区县农村贫困发生率明显呈现下降态势，区县间的差距在缩小.

运用Geoda软件对重庆市县域农村贫困化进行空间自相关分析，得到基于2014年到2018年重庆市33个具有扶贫开发任务的区县农村贫困化的全局自相关莫兰指数(图 2)，图 2显示，2014-2018年重庆市县域农村贫困化莫兰指数均为正数，且通过了显著性检验，其中，2014年莫兰指数为0.839 788，2018年莫兰指数为0.563 036，可以看出，研究时段内重庆市县域农村贫困化在空间上趋于集聚分布，即农村贫困化高的区县显著集聚，也表明重庆市区县农村贫困在空间上的分布不是随机的，而是存在空间集聚效应，即农村贫困发生率高与低的区县都在地理空间上显著集聚；同时比较莫兰指数发现，研究时段内县域农村贫困化莫兰指数总体呈下降趋势，表明重庆市县域农村贫困化与空间分布的相关性在减弱.各县域农村贫困化在空间上并不是独立的，而是呈现一定的空间集聚特征，这为后面研究县域农村贫困化空间分异影响机制采用GWR模型奠定了基础，也为模型结构的有效性提供了必要的保障.

全局空间自相关揭示重庆市各区县农村贫困化在整体上呈现出显著的空间自相关性，为进一步揭示县域农村贫困化的高值和低值空间集聚状态，了解局部的空间差异性，需要进行局部空间的关联格局分析.根据局部Moran's I统计值，将重庆市33个研究单元划分为4种类型：①高高集聚区(H—H)，表明区县自身与相邻区县的农村贫困化水平高，两者之间的空间差异程度较小；②高低集聚区(H—L)，区县自身农村贫困化水平高，但相邻区县的农村贫困化水平叫低，呈现出负的相关性；③低高集聚区(L—H)，区县自身农村贫困化水平低，但相邻的区县农村贫困化水平高，也呈现出负的空间相关性；④低低集聚区(L—L)，区县自身的农村贫困化水平和相邻区县的农村贫困化水平都较低，两者空间差异较小，呈现出显著的正相关.

运用Geoda软件绘制重庆市2014-2018年县域农村贫困化的Moran散点图，图 3象限按其性质分为“H—H”关联(第一象限)、“L—H”关联(第二象限)、“L—L”关联(第三象限)、“H—L”关联(第四象限)，其中第一、三象限对应的空间关联是正的空间自相关，即有均质性，第二、四象限对应的空间关联是负的空间自相关，即有异质性. 2014年县域农村贫困化散点落入第一象限主要为国家级贫困县和市级贫困县，主要分布于渝东南和渝东北地区，涉及的区县有城口、丰都、奉节、开州、南川、彭水、黔江、石柱、万州、巫山、巫溪、武隆、秀山、酉阳、云阳、忠县，共16个区县，占比48.48%，落入第三象限的散点主要为非重点区县，分布于环主城周边区县，涉及的区县有巴南、北碚、璧山、大足、垫江、合川、梁平、綦江、荣昌、铜梁、潼南、永川、渝北、长寿、江津，共15个区县，占比45.45%，合计空间正相关的区县有31个，占93.94%；2018年落入第一、第三象限的散点有所降低，表明重庆市农村贫困接近的区县在空间上的集聚效应在不断减小.

由于全局Moran's I指数为总体自相关统计量，并不能表明具体地区的空间集聚特征强度，要研究重庆市农村贫困化是否存在局部集聚现象，则需用Moran散点图和局部Moran's I指数.在一定的显著性水平下，2014年高—高值集聚区主要分布在渝东北和渝东南地区，包括开州、云阳、巫山、巫溪、奉节、丰都、彭水、黔江、酉阳9个区县，均为贫困程度较深的国贫县，表明这些区县自身与相邻区县的农村贫困化水平高，两者之间的空间差异程度较小，联系紧密，辐射带动作用较强，会受到相邻区县贫困化的正向影响. 2014年至2018年，高—高值集聚区存在着明显的由渝东南和渝东北双中心逐渐向渝东北单中心转移的趋势，渝东北高—高值集聚区新增城口县，通过扶贫开发工作的推进，渝东南高—高值集聚区则逐年依次减少了丰都县、酉阳县.值得注意的是，秀山县于2016年新增为高—高值集聚区，并在2017年转变为低—高值集聚区，表明近5年来，秀山县扶贫压力较大. 2014年低—低值集聚区主要分布在主城区及其周边区县，包括合川、铜梁、北碚、渝北、璧山、沙坪坝、渝中、江北、九龙坡、南岸、大渡口、永川、江津、巴南14个区县，表明这些区县自身的农村贫困化水平和相邻区县的农村贫困化水平都较低，两者空间差异较小，集聚效益明显，呈现出显著的正相关. 2015年至2018年，低—低值集聚区区县数量逐渐减少，由14个区县减少为11个区县. 2015年新增大足区，2016年减少合川区、铜梁区和大足区，2017年减少永川区，新增合川区，2018年减少合川区和江津区.合川区、大足区变化较为频繁，在低—低值集聚区与无值区不断转换，表明合川区与大足区贫困化水平不稳定，值得关注.

衡量农村贫困化的指标有农民人均纯收入、农村贫困人口比重指标、贫困人口收入与贫困线差距的总和与达到贫困线收入的总和的比率、地区居民收入差距指标等^[24].其中农村贫困人口比重指标是目前比较常用的统计县域农村贫困化的方法，可以最直接、最客观反映县域农村贫困化水平^[9].因此本研究采用2014年重庆市县域农村贫困人口比重(农村贫困发生率)作为衡量县域农村贫困化的指标，即因变量.

本研究在参考前人研究成果的基础上^[24-26]，综合考虑指标的系统性、重点性和可获得性，结合重庆市社会经济发展和区域贫困特性，从自然地理环境、经济和社会3个维度构建县域农村贫困化空间分异的影响因素指标体系(表 1).

在进行回归分析之前，为探索影响因子对农村贫困化的解释能力，筛选关键指标，本研究运用SPSS 22.0软件，以2014年县域农村贫困发生率为因变量、13项影响因素指标为解释变量，计算Pearson相关系数，检验各影响因素与农村贫困发生率的线性关系，并根据自变量的方差膨胀因子(Variance Inflation Factor，VIF)诊断影响因素之间的多重共线性情况.

从表 2可以看出，13个解释变量中与县域农村贫困化存在显著相关性的有海拔高度(X₁)、地形起伏度(X₂)、喀斯特岩溶地貌面积比重(X₃)、人均耕地资源面积(X₄)、25°以上耕地面积比重(X₅)、区位水平(X₆)、未通客运班车村比重(X₇)、城镇化率(X₉)、文盲人口比重(X₁₀)、少数民族人口比重(X₁₁)和因残致贫人口比重(X₁₃)，以上自变量均在0.01的水平上显著，而外出务工人口比重(X₁₂)未能通过相关性检验，予以剔除；在ArcGIS的GWR模型中，当VIF(方差膨胀因子)值高于10时，多重共线性会使模型变得不稳定.因此以VIF值＞10为解释变量多重共线性的界定标准，因此剔除地形起伏度(X₂)、喀斯特岩溶地貌面积比重(X₃)、城镇化率(X₉)和少数民族人口比重(X₁₁).

运用Arc GIS 10.2软件的GWR模块，以各调查村的地理中心为中心坐标构建模型，对表 3中显著相关的8个因素进行估计，重点分析各影响因素对农村贫困化在空间区位上的表现差异.在GWR模型中，每一个单元都有特定的系数，表 3是对各研究单元回归系数的描述性统计，给出了各影响因素对农村贫困化影响程度的平均值、最小值、下四分位数、中位数、上四分位数和最大值.从回归系数的情况来看(表 4)，农村贫困化影响因素的回归系数有正有负，从最大值和最小值看，所有解释变量在不同的空间单元上对因变量的影响作用性质均一致.

1) 海拔高度. GWR模型结果显示(图5)，研究区内海拔高度与农村贫困发生率呈正相关，即平均海拔高度越大的区域，其贫困发生率越高.海拔高度的回归系数呈现由西南向东北逐渐增大的分布特征，从贫困类型分区上来看，深度贫困区受海拔高度因素影响更为显著.海拔高度作为整体地理条件，决定了县域气候、降水、地形地貌、生产资源条件等其他自然地理状况，进而影响区域社会、经济发展.高度贫困区主要位于渝东北地区和渝东南地区，两地地形尤为复杂，平均海拔较高，对社会经济发展和农户生产约束性更强.

2) 人均耕地资源面积.人均耕地资源面积对重庆市县域农村贫困化有正向影响，其回归系数在空间上呈现出由西北向东北和东南2个方向逐渐增大的趋势，贫困类型分区上，深度贫困地区回归系数较大，主要是由于该类地区喀斯特岩溶地貌广布，地形起伏度大，农业灌溉条件差，无法将人均耕地资源数量优势转化为质量优势所致；另外，人均耕地资源面积回归系数也呈现出区域间波动不大的特点，也说明当前重庆市贫困农户收入来源不仅是来源于种植业，人均耕地资源对县域农村贫困化空间分异影响有限.

3) 25°以上耕地面积比重和区位水平. 25°以上耕地面积比重和区位水平对重庆市县域农村贫困化呈正相关. 2个回归系数在空间上均呈现西南低、东北高的分布特征，且存在一定的波动，从贫困类型分区上来看，深度贫困地区受25°以上耕地面积比重和区位水平影响更为显著. 25°以上耕地面积比重过高，则农业资源利用的条件也较差.距离市中心的距离远近，反映出县域对外连接的便利程度，地处偏远的区县，由于交通时间成本更高，外界的资本、技术等资源不易进入，农户陷入贫困的几率也更大.

4) 未通客运班车村比重.从未通客运班车村比重回归系数来看，该因素与重庆市县域农村贫困化存在正向影响，但从回归系数上来看却呈现出由东北向西南递减的趋势，再看贫困类型分区，低度贫困区和中度贫困区回归系数较深度贫困区大，与实地调查结论相符，这主要是由于深度贫困地区百姓出行更多依赖于公共交通工具，因此，深度贫困地区未通客运班车村比重较高；相反，低度和中度贫困区对边远村的交通连接未能引起重视，改善程度不够，导致回归系数较高.因此，为防止农村贫困化“孤岛效应”的深化，扶贫过程中应注重改善公共交通，尤其加强低中度贫困区的公共交通建设，提高减贫效应，有效防止“灯下黑”.

5) 农村居民人均纯收入.农村居民人均可支配收入与农村贫困化呈负向关联作用，其回归系数的绝对值在空间上由西南向东南递增，但是波动较小，不同贫困类型分区受该因素影响的差异不大.

6) 文盲人口比.文盲人口比与农村贫困化呈负向关联作用，其回归系数在空间上的分布呈现由西南向东北逐渐增大的分布特征，从贫困类型分区上来看，深度贫困地区受文盲人口比影响更为显著，文盲人口比反映出区域人口受教育情况，实地调查发现贫困户未上过学的现象集中于深度贫困地区，且总的来看全市受访贫困人口教育文化水平以初中及小学文化水平为主.因此，深度贫困地区在推进扶贫工作过程中，应高度重视教育扶贫工作，通过教育来斩断贫困代际传递.

7) 因残致贫人口比重.因残致贫人口比重，与农村贫困化呈负向关联作用，其回归系数在空间上的分布呈现由东北向西南逐渐增大的分布特征，从贫困类型分区上来看，低度和中度贫困区回归系数较深度贫困地区要高，这主要是低、中度贫困地区贫困人口属“插花式”个体贫困，其中因残致贫的农户比例较高.因此，低、中度贫困地区在推进扶贫工作过程中，高度重视此类特殊困难群体的扶贫和社会保障工作，做到“应扶尽扶，应保尽保”.

综上，通过显著性检验的解释变量对县域农村贫困化起到了直接或间接的影响作用，且在不同类型的区域上表现出的影响力也不尽相同.各区域应高度关注显著影响区域农村贫困化的因素，在推进扶贫工作过程中，做到有的放矢，提高扶贫工作效率.

本研究开展县域尺度农村贫困化的时空分异和影响因素2个方面的内容研究，刻画了农村贫困化的空间关联格局，构建了影响县域农村贫困化空间分异的指标体系，并运用空间计量模型对县域农村贫困化的空间分异影响因素进行了回归分析和探讨.研究结果总体符合研究区实际情况，具有较强的理论和实践指导意义.

1) 在Arc GIS 10.2软件平台支撑下对重庆市2014-2018年县域农村贫困化进行空间表征，从总体上来看，研究时段内重庆市县域农村贫困化热点区域整体空间格局稳定，高值区主要分布于巫溪、城口、酉阳和彭水等渝东北和渝东南区域的相关区县，而主城区和渝西地区的区县贫困发生率相对较低，呈现显著的空间差异性；从时间序列上来看，研究时段内各区县农村贫困发生率明显呈现下降态势，区县间的差距在缩小.

2) 基于ESDA方法，重庆市2014-2018年全局Moran's I指数均达到显著的空间正相关，县域农村贫困化呈现出明显的空间集聚现象. 2014-2018年，全局Moran's I值总体呈下降趋势，表明重庆市县域农村贫困化与空间分布的相关性在减弱；高高集聚区分布在渝东北和渝东南地区，低低集聚区主要分布于主城和环主城周边区县；在时间序列上，空间聚集的高值中心逐渐由双中心向单中心过渡，低值中心空间分布变化相对较小，部分区县逐步退出低值中心，呈现不显著状态，表明随着扶贫开发工作的深入推进，部分区县农村贫困化水平有所下降.

3) 县域农村贫困化空间分异影响因素中，通过Peason相关性分析和多重共线性诊断，剔除了地形起伏度、喀斯特岩溶地貌面积比重、城镇化率、少数民族人口比重、外出务工人口比重等因素，通过GWR模型进行回归分析发现，8个显著的解释变量对县域农村贫困化起到了直接或间接的影响，其中海拔高度、人均耕地资源面积、25°以上耕地面积比重、区位水平和未通客运班车村比重与县域农村贫困化呈现正相关，农村居民人均纯收入、文盲人口比和因残致贫人口比重与县域农村贫困化呈负相关.

对重庆市县域农村贫困化空间分异与影响因素的分析可以看出，深度贫困地区受自然地理条件和资源禀赋等因素的制约，亟需优先补齐发展过程中的短板，如改善深度贫困地区农田水利条件，加强贫困地区与外界的沟通和联系，提升贫困群众文化素质，激发群众内生动力；另外在低度和中度贫困区域需进一步关注农村公共交通的改善和特殊困难人群的保障工作，确保各项措施能够精准到人到户，最终实现贫困人口持续稳定脱贫.

本研究在县域尺度上对农村贫困化的空间分异影响因素进行了探讨分析，对县域脱贫攻坚工作具有宏观指导意义，下一步的研究可进一步缩小研究尺度，在镇域尺度或村域尺度探讨农村贫困化空间分异的具体影响因素；另外，本研究探讨了单个因素对于农村贫困化的影响，下一步的研究可进一步探讨因素之间的相互作用对农村贫困化的影响.