苟凯歌, 蒋 辉, 刘兆阳

(吉首大学 乡村振兴战略研究中心 湖南 吉首 416 000)

水资源贫困与经济贫困存在天然的内在联系，中国农村水资源贫困与经济贫困耦合度很高，提高水资源的利用效率，加强水资源的可持续管理，提升水资源的经济效益对解决贫困问题，促进经济社会的发展有重要的意义。中国的人均水资源拥有量只有世界平均水平的25%，是全球人均水资源最贫乏的国家之一。水资源的贫困在很大程度上影响了经济社会发展，在一些水资源极度贫困的地区尤其如此。虽然中国农村的反贫困工作上取得了历史性突破，但绝对贫困与相对贫困依然是今后需要面临的长期性问题。水资源的匮乏，农村水资源贫困是制约农村发展以及反贫困工作的巨大障碍。同时，由于城市化快速推进，农村面源污染严重，加上农村用水技术落后，用水效率不高，使得原本的农村水资源贫困问题更是“雪上加霜”，农村经济发展面临的挑战依然严峻。因此，正确认识农村水资源贫困与经济贫困的相互作用关系，可以有效实现水资源的可持续管理，提高水资源的利用效率和经济效益，从而找到一条协同减缓水资源贫困与经济贫困的良性发展路径。

近年来，水贫困评价指数(water poverty index，WPI)用于衡量水资源贫困程度得到普遍认可[1-5]。中国学者的研究集中于水资源贫困的现状、地区差异以及发展趋势等方面[2-3,5-8]。并将水贫困理论应用于农村生产活动中，不仅丰富和发展了水贫困理论，更为解决农村水资源可持续性问题提供了新的思路。学者们或构建水资源支撑指数(WSPI)来描述水资源短缺程度，探讨水资源条件对区域社会经济发展的支撑能力以及中国水资源支撑能力空间格局，或对中国31个省(市、区)农村地区的水贫困与经济贫困进行评价，验证了31个省(市、区)农村水贫困与经济贫困之间存在共生关系[5]，或建立WPI-EPI模型，从时序、空间分异两个角度对水贫困与经济贫困进行测度和研究[2]。《国家农业节水纲要(2012—2020年)》提出农村水资源的供需矛盾依旧是制约农村经济快速发展的瓶颈，随着中国城镇化进程加快、产业结构升级，农村耕地面积开始减少，农村用水量受到显著影响，农民为追求利润最大化，优化农作物种植结构或变动农产品价格，使得农村地区水资源与经济协调发展的压力升级；同时，国家加大农林水支出力度，为推动农村经济快速发展、提高农村用水效率起到积极作用。

本研究以中国31个省(市、区)为研究对象，在借鉴前人有关水资源贫困和经济贫困耦合协调研究成果的基础上，构建适用于中国国情的农村水资源贫困和经济贫困的评价指标体系，计算农村水资源贫困指数与经济贫困指数，运用系统耦合协调度模型，测算出农村水资源贫困指数和经济贫困指数的耦合度和耦合协调度，并运用空间杜宾模型找出影响耦合协调的重要因素，旨在为中国减轻农村水贫困与经济贫困提供理论依据和政策启示。

1 数据来源、研究方法和指标体系

1.1 数据来源

本文数据主要构建包括人均水资源量、节水灌溉类机械拥有量、化肥施用量等16个指标在内的农村水资源贫困评价指标体系所需要的数据，和构建包括农村居民家庭人均纯收入、农村居民家庭恩格尔系数、每千农村人口乡镇卫生院人员数等14个指标在内的农村经济贫困评价指标体系所需要的数据，相关数据为中国31省(市、区)2000—2017年共18 a的数据。数据主要通过官方公布的《中国统计年鉴》《中国环境统计年鉴》《中国农村统计年鉴》《中国水利统计年鉴》《中国经济社会发展年鉴数据》《中国财政年鉴》《中国卫生和计划生育统计年鉴》和《中国人口和就业统计年鉴》等收集整理。部分缺失数据采用周围地区近似替代法、临近年份插值法，Spssau回归估计法得到，以求数据的完整性和可靠性。

1.2 研究方法

1.2.1 水贫困测度模型 水贫困指数(WPI)是定量评价一个国家或地区相对缺水程度的一组综合性指标，可以反映一个地区水资源实际情况，以及工程、管理、经济、人类福利与环境状况。相关测度模型为：

WPI=wrR+waA+wcC+wuU+weE

(1)

在实际应用过程中，由于原始数据存在数量级间的差距，因此在计算过程中，对原始数据取其对数，然后进行无量纲化处理，本文无量纲化处理方法采用极值法，其中，对数值越低水贫困(或经济贫困)越深的指标带入正向指标计算公式，反之带入负向指标计算公式。

正向指标计算公式：

(2)

负向指标计算公式：

(3)

1.2.2 经济贫困测度模型 经济贫困是一个与发展环境、权利、机会、自然资源禀赋、经济水平紧密相关的综合现象，采用EPI方法更能客观全面反映出经济贫困的真实内涵，其测度模型为：

EPI=wiI+weE+wdD+whH+wvV+wsS

(4)

1.2.3 主客观综合赋权法 在本文中，主观赋权法选用层次分析法，客观赋权法选用熵值法。层次分析法主要是根据专家知识和已有经验来确定指标权重，具有比较大的主观性和随意性；基于调查数据和一定的数学模型，通过计算确定指标的重要程度，熵值法有很强的数理依据，却忽视了评价人员的经验。本文将两种赋权方法结合起来，克服单一赋权法的不足，使指标赋权更加客观、准确。借鉴Komnenic等[6]研究，其计算公式为：

(5)

参考张海涛等[7]关于主客观权重的研究,本文将采用α=0.3，β=0.7来计算综合权重。

1.2.4 水资源贫困与农村经济贫困的耦合协调度模型 考虑到水资源贫困与经济贫困的开放性、非平衡性和非线性特征，对二者的耦合程度测量主要借鉴现有成熟的方法，具体计算公式如下：

(6)

T=aWPI+bEPI

(7)

(8)

式中：C为水资源贫困与经济贫困耦合度,取值[0,1]。C越大，表明水资源贫困与经济贫困耦合度越高；WPI是水资源贫困指数；EPI是经济贫困指数；a，b为待定系数，参考孙才志等[8]关于耦合度模型的研究，本文认为解决水资源贫困问题和经济贫困问题同等重要，因此令a=b=0.5；k为调节系数，为增加区分度，令k=2；D为农村水资源贫困与经济贫困的耦合协调度，D取值[0,1]，D越大，表明水资源贫困与经济贫困耦合协调度越高。

1.2.5 农村水资源贫困经济贫困的耦合类型和耦合协调类型划分 本文参照Shi Tao等[9]关于耦合度与耦合协调度的研究成果，将耦合度分为4类；本文根据Sullivan等[10]关于耦合协调度的研究将耦合协调度分为4类(表1—2)。

表1 耦合类型划分及判断标准

表2 耦合协调类型划分及判断标准

1.2.6 耦合协调度的空间相关性检验 利用莫兰指数(Moran’sI)可以表明水资源贫困与经济贫困的耦合协调度的空间相关性，测度模型如下：

(9)

1.2.7 空间杜宾面板数据模型 构建空间面板数据模型可以分析影响农村水资源贫困与经济贫困的耦合协调程度的重要因素，以揭示耦合协调程度时空变化差异的成因。而由Elhorst首次提出的空间杜宾面板数据模型(SDM)，在空间滞后面板数据模型和空间误差面板数据模型的基础上引入被解释变量和解释变量的空间滞后项，可以更好的估计基于面板数据测得的空间效应。

Y=ρWY+β1X+β2WX+ε

(10)

式中：Y为省(市、区)农村水资源贫困与经济贫困的耦合协调度;WY是Y的空间滞后项;ρ用来衡量Y的溢出效应;X为一系列影响耦合协调度的因素，当β1反映X对Y的影响程度;WX是解释变量X的空间滞后项，β2用来衡量邻近省级区域的解释变量X对本省级区域Y的空间影响程度;ε是随机误差项向量;λ是空间误差相关系数。

1.3 指标体系构建

1.3.1 水资源贫困评价指标体系构建及权重确定 本文主要通过构建水资源贫困评价指标体系，来探讨农村水资源使用与可持续发展的协调关系。已有学者针对WPI模型构建了一套完整的评价指标体系，水资源贫困评价体系包括了由资源、设施、能力、使用和环境5个子系统的系列变量组成的指标[8]。水资源状况是指可以被利用的地表及地下水资源量及其可靠性或可变性；设施状况指自来水灌溉的普及率，节水、用水的设施水平，反映社会大众接近清洁水源的程度以及用水的安全性；利用能力考虑农村居民经济能力对水行业的影响；使用效率综合反映生活和农业各部门的用水效率；环境状况反映与水资源管理相关的环境状况。结合中国农村生产实际情况，在现有的农村水贫困指标体系上增添或删减子系统指标。运用层次分析法(根据专家知识和已有经验来确定指标权重)和熵值法(基于调查数据和一定的数学模型，通过计算确定指标权重)求得各评价指标的主观权重和客观权重，采用主客观综合赋权法求得综合权重(表3)。

1.3.2 经济贫困评价指标体系的构建及权重确定 经济贫困指标体系(EPI)主要由收入水平、支出水平、教育水平、医疗卫生、发展环境、家庭及就业6个子系统系列变量指标组成。在具体指标确定上，运用历史文献统计法，得到使用最频繁且相对有效的指标，结合中国31个省(市、区)农村实际情况，构建适合精准测度31个省(市、区)农村地区经济情况的EPI体系，以满足长期分析需要[11]。

在经济贫困指标体系中，收入水平由农村居民家庭人均纯收入来表征；支出水平由农村家庭人均消费性支出与农村居民家庭恩格尔系数来表征；教育水平主要用农村普通高中毕业生数、农村人口文盲率与农村人均教育经费来表征；医疗卫生的指标为每千农村人口乡镇卫生院人员数、乡镇医院床位数和乡镇卫生院数；发展环境的指标为农作物社会救济费、农作物受灾面积、农作物自然灾害救济费、通有线电视村比重、通自来水村比重和通公共交通农村比重。家庭及就业由第一产业从业人员占全社会从业人员比重与农村平均家庭户规模来表征。运用层次分析法(根据专家知识和已有经验来确定指标权重)和熵值法(基于调查数据和一定的数学模型，通过计算确定指标权重)求得各评价指标的主观权重和客观权重，采用主客观综合赋权法求得综合权重(表4)。

表3 中国省域农村水贫困评价指标体系及指标权重

表4 中国省域经济贫困评价指标体系及指标权重

1.3.3 解释变量说明 农村水资源贫困与经济贫困的相互作用的复杂性决定了二者的耦合协调受多种因素影响。从外部因素来看，国家对于农林水事务的财政支出是二者协调发展的调控推手，对于提高农村用水效率、改善农村基础设施等具有正向作用；城镇化进程加快对解决农村剩余劳动力问题，进而提高农民经济收益，拉动农村现代化发展，繁荣农村经济具有重要意义。从内部因素来看，农作物是农村生产与农民增收的重要经济纽带，农产品价格的大幅度变动，或不考虑当地实际情况，过度增加农作物种植面积，会损坏粮食安全，进而负向影响农民收入；产业结构优化对于加快当地经济快速发展，以及实现水资源优化配置都有重要意义。各解释变量如下：

(1) 城镇化率X1。Sun[12]提出城镇化率不断增加，城镇人口增加，城市用地规模扩大，导致耕地面积减少，进而影响农村用水。

(2) 农作物种植结构X2。该指标是指实际播种或移植有农作物的面积。王洁萍等[13]认为农作物在种植全过程中要消耗水资源，农作物播种面积的增加，会负向影响农村用水效率。借鉴王洁萍的做法，采用粮食作物播种面积比重来表示。

(3) 农产品价格X3。Ritu等[14]认为农民在追求利润最大化的驱动下，农产品价格是影响农民生产活动成本和预期收益，从而影响农民对水资源及其他要素的使用及种植方式的确定。

(4) 产业结构X4。借鉴栗清亚等[15]的研究，产业结构调整会影响农村用水量和用水结构，有利于节约农村水资源。Peneder M[16]认为提高产业结构合理化，促进农村人口向第二、三产业转移就业，促进转移人口的收入增长，有助于提高农村经济发展。具体计算指标为(二、三产业增加值)/GDP。

(5) 中国财政农林水事务支出X5。Ramakumar R[17]提出农林水事务支出是政府增加农业投入，保护农村发展的有效手段，推动农村经济快速发展。

综上所述，借鉴以上学者研究经验，考虑到全国性数据的可获得性，最终决定采用城镇化率、农作物种植结构、农产品价格、产业结构以及中国财政农林水事务支出作为农村水资源贫困与经济贫困耦合协调发展的影响因素。

2 模型结果与分析

为避免周期性波动，本文将研究的时间单元及相应的计算结果划分为2000—2004，2005—2008，2009—2013，2014—2017年4个时间段。将各子系统权重置入WPI，EPI模型，得到中国农村水资源贫困指数WPI和经济贫困指数EPI,利用双系统耦合模型得到各个时间段的耦合度和耦合协调度的均值。根据本研究对2000—2017年的4个研究单元，结合耦合度和耦合协调度计算结果及其分类标准，得到中国农村水贫困和经济贫困耦合特征及其时空演化趋势。

2.1 农村水资源贫困与经济贫困耦合度较高，且呈现上升趋势

整体来看，2000—2017年，31个省(市、区)水资源贫困与经济贫困耦合度呈上升趋势，虽然在2009—2013年有小幅下降，但在考察的18个年份间均达到中度耦合，且后3个时期均高于基期(2000—2004年)水平(表5)。

2000—2004年，有21个省(市、区)的农村水资源贫困与经济贫困耦合度超过平均水平，主要分布于北部和西部地区，有27个省(市、区)呈高度耦合状态；2005—2008年，共有28个省(市、区)的耦合度高于平均水平；2009—2013年，共有28个省(市、区)农村水贫困与经济贫困耦合度超过平均水平，呈东多西少分布；2014—2017年，共有28个省(市、区)的耦合度高于平均水平，与2000—2004年相比，河北、辽宁、浙江、福建、山西、海南、湖南、湖北、宁夏等18个省(市、区)农村水资源贫困与经济贫困的耦合度仅出现0～0.020范围内的微小变动。至2017年，除北京、天津、上海3个省(市、区)，其余省(市、区)耦合度均在0.800以上。表明水资源贫困与经济贫困二者间存在较强的相互关联性，改善农村水资源贫困能够有效促进农村贫困的缓解，且作用强度较大。

表5 中国31个省(市、区)农村31个省份农村水资源贫困与农村经济贫困耦合度

2.2 中国的农村水资源贫困与经济贫困的耦合协调度均为中度协调水平

由表6可以看出，2000—2017年，中国31个省(市、区)农村水资源贫困与经济贫困的耦合协调度均处于0.500～0.790之间，表明中国各省(市、区)耦合协调度均为中度耦合协调型。从时间跨度上来看，整体耦合协调度没有明显变化，除内蒙古、青海和陕西3省，北方其余各省(市、区)耦合协调度仅在0.001～0.023范围内微增，提升度普遍低于南方地区。空间差异上，中国农村水资源贫困与经济贫困的耦合协调度相对较高的地区主要集中在东北和西南地区，东南地区协调度相对较低。

耦合度反应了农村水资源贫困与经济贫困之间相互影响的强弱程度，耦合协调度则较全面的表达了二者相互作用过程中良性耦合程度的大小，体现了协调状况的好坏，可以反映出二者是在高水平上相互促进还是在低水平上相互制约；耦合协调度是由耦合度与耦合协调发展水平共同决定的，因耦合度变动幅度不明显，所以农村水资源贫困与经济贫困的协调发展水平成为影响耦合协调度的最终决定性因素，决定了耦合协调度的高低。研究表明，2000—2017年，可能主要因为近年来南水北调工程、西部大开发战略、乡村振兴战略的实施，一定程度上缓解了部分西部地区农村水资源匮乏的情况，同时也推动了西部地区经济社会快速发展，国家水利发展迅速，中国省级层面农村水资源配置水平与经济发展水平已经形成较强的协调态势，二者存在明显的共生关系，缓解水资源匮乏可以有效减缓经济贫困的程度(表6)。

表6 中国31个省(市、区)农村水资源贫困与农村经济贫困耦合协调度

2.3 中国农村水资源贫困与农村经济贫困的耦合协调度在空间上正向相关

由表7可知，所有年份的Moran’sI值均大于零且Z检验显著,在90%以上的水平，这表明中国农村水资源贫困与经济贫困的耦合协调度在空间上具有明显的正自相关关系，在空间上呈集聚分布态势，农村水资源贫困与经济贫困的协调发展水平不但与自身因素有关，同时亦与周围地区存在空间依赖性。由表7第二列可以看出，Moran’sI值从2000年的0.245降至2008年的0.187，在2017年达到最高值(0.456)，随着时间的推移，耦合协调度的空间分布的相关性先减弱后升高。

表7 各省级区域农村水资源贫困与农村经济贫困耦合协调度的Moran’s I指数及其显著性

2.4 中国农村水资源贫困与经济贫困耦合协调度的影响因素分析

本文选取了城镇化率、粮食作物播种面积比重、农产品价格、产业结构以及中国财政农林水事务支出5个要素入手，来探究其对农村水资源贫困与农村经济贫困的耦合协调发展的作用程度。空间杜宾面板数据模型检验包含Wald检验和LR检验，借鉴鲍超等[18]的研究，本文采用LR检验来判断SDM是否可以简化成SLM或SEM。LR检验结果分别为3.097，110.065，分别在10%和1%的水平上显著，即SDM不适合简化为SLM与SEM。

(1) 空间效应分析。从整体来看，空间杜宾模型的 达到了0.935，显示模型的整体拟合程度较高，解释变量与被解释变量之间有较好的拟合关系，能够较全面的表征二者之间的关系。空间自回归系数Rho为-0.025 3，在1%显著水平上显著，表明耦合协调程度存在空间溢出效应，一个省级区域农村水资源贫困与农村经济贫困耦合协调度的改变将负向影响相邻省(市、区)，使相邻省级区域的耦合协调度相应降低。

在模型结果中，回归系数的显著性检验表明，各省(市、区)农村水资源与农村经济的耦合协调度不同程度的受到当地农产品价格X3，产业结构X4，财政农林水事务支出X5这5个指标的影响，而城镇失业率X1，粮食作物播种面积比重X2对耦合协调度影响不明显。从影响显著的3个指标的回归系数值来看，农产品价格X3，产业结构X4的回归系数分别为-0.011 5，-0.083 3，即负向影响；而财政农林水事务支出X5的回归系数为0.028 8，对耦合协调度产生正向影响。城镇化率X1，产业结构X4与财政农林水事务支出X5这3个解释变量的空间滞后项系数为-0.019 9，0.037 8，-0.004 1，均在1%水平上显著，反映出相邻省(市、区)城镇失业率、产业结构与国家对当地农林水的财政支出的改变会影响到本省级区域农村水资源贫困与农村经济贫困的耦合协调度(表8)。

表8 空间杜宾面板数据模型回归估计结果

(2) 影响因素分析。为了能够对空间杜宾面板数据模型回归系数进行合理的解释，勒萨热和佩斯提出了空间回归模型偏微分方法，由直接效应、间接效应来反映解释变量X与被解释变量Y之间的相互关系(如表9所示)。

财政农林水事务支出X5的直接效应与间接效应分别在10%与1%的水平上显著，系数分别为0.032 0，-0.003 3，表明该指标对当地省级区域以及邻近省级区域的农村水资源贫困与农村经济贫困耦合协调度均产生影响，各省级区域国家财政农林水事务支出提高1%该省(市、区)的农村水资源贫困与农村经济贫困耦合协调度提高0.032%，同时邻近省(市、区)降低0.003 3%。

表9 各解释变量对协调度的效应

农产品价格X3和产业结构X4这两个指标的直接效应显著而间接效应不显著。农产品价格X3，产业结构X4的直接效应系数分别为-0.128 1，-0.926 0，分别在1%和5%的水平上显著，而间接效应不显著，表明农产品价格的提升和二三产业占比的提升会负向影响该省(市、区)的农村水资源贫困与农村经济贫困耦合协调度，但对邻近省(市、区)无影响。

总体上看，财政农林水事务支出X5具有较为显著的直接效应与间接效应，是影响农村水资源贫困与农村经济贫困耦合协调度的重要因素。

3 结论与建议

本研究运用耦合协调度模型探讨了2000—2017年中国内陆31个省(市、区)农村水资源贫困与经济贫困的时空分异演变趋势，得到以下主要结论：

(1) 相互影响明显。中国水资源贫困与经济贫困的耦合度基本趋于稳定，均处于高度耦合状态。这说明，目前中国农村水资源可用程度与农村经济发展水平之间存在明显的相互依存、彼此影响的关系。即农村水资源匮乏的地区极易产生经济贫困，反之，经济贫困地区也更容易出现农村水资源贫困的现象。

(2) 呈现协调互动态势。协调发展程度表现为中度协调。这意味着一个地区的农村水资源贫困得到改善时，其经济贫困问题也会得到相应改善，同样，当一个地区经济得到发展的同时，其水资源贫困现象也会得到一定改善，二者呈现出共同发展的良性态势。其中，东部沿海地区的江苏、浙江、福建和广东，中部地区的湖南省等几个省(市、区)的耦合协调度最高，这可能主要因为上述省(市、区)先天水资源条件较好，且水利事业发展迅速。相比之下，广大西部地区的耦合协调度相对较低，但随着近年来南水北调工程、西部大开发战略、乡村振兴战略的实施在一定程度上缓解了部分地区农村水资源匮乏的情况，同时也推动了西部地区经济社会快速发展，从而缩小了这些地区与沿海发达地区在农村用水和经济发展方面的差距。

(3) 耦合协调发展过程中的Moran’sI值均在90%显著水平下通过检验，在空间格局上具有正向空间自相关性。在分析影响因素时，空间计量模型与基准面板数据回归模型相比，考虑到空间滞后与空间误差项的综合作用，可以更准确的提取出耦合协调度影响因素，且空间杜宾模型对其模拟最优。

(4) 协调发展水平存在显著的负向溢出效应。相邻省(市、区)协调发展程度的提升会扼制本省(市、区)协调协调水平发展。财政农林水事务支出对各省级区域农村水资源与农村经济协调发展水平具有正向促进作用，农产品价格、二三产业比重过大则不利于农村水资源与农村经济协调发展。合理调整农产品价格、发挥产业集聚效应、增加国家财政对于农林水事务支出对于实现农村水资源与农村经济协调发展具有重要意义。

基于以上结论，为推动中国农村水资源贫困与经济贫困的耦合协调发展，应从市场角度出发，合理调整农产品价格。国家农业相关部门应制定农产品相关的政策与措施，避免农产品垄断现象，加大技术投入，使农业进入技术化与规模化。坚持市场化改革，基于科技进步的结构优化，加强地区间分工与合作，可使农村富余劳动力转移力度不断加大，促进产业升级。国家应加大或维持对农村原有的财政支持力度，同时要不断优化和完善经济政策工具，充分发挥政府和市场在资源配置中的调节作用，确保农村水资源的合理使用。