农业水土资源即维持农业生产所必需的水资源和耕地资源，是人类赖以生存的基础资源和社会经济发展的物质载体。随着我国工业化和城镇化发展持续推进，农业水土资源所表现出的挤占、短缺和区域性结构破坏等资源非持续利用问题日益突出，对国家粮食生产安全形成重大威胁^[1-2]。因此，开展农业水土资源承载力评价和水土平衡研究具有重要的科学意义。农业水土资源承载力研究一直是土地科学、水资源科学等学科关注的热点^[3-5]，但前期多是将农业水资源和耕地资源作为独立系统，分别进行水资源承载力研究和土地资源承载力研究，这种思路忽视了农业水土资源复杂的综合作用与耦合关系，不利于水土资源系统整体效能的充分发挥^[6-7]。近年来，农业水土资源承载力综合研究得到越来越多学者的重视，产生了一些有价值的研究成果。施开放等^[8]利用熵权可拓决策模型对重庆三峡库区水土资源承载力进行了定量评价；陈红梅等^[9]构建时序全局因子分析模型研究了黑龙江省水土资源承载力空间分异特征；李慧等^[10]结合粒子群优化算法和投影寻踪模型对延安市各区县农业水土资源承载力进行了评价；南彩艳等^[11]针对水土资源承载系统的不确定性，采用改进集对分析对关中地区水土资源承载等级进行了研究。

农业水土资源系统是一个非线性的复杂巨系统，其影响因素具有复杂性和动态性特征，所以目前尚无统一的评价指标和评价模型。本文在全面分析湖北省水土资源禀赋和社会经济发展现状的基础上，将农业水土资源视为一个统一系统，采用DPSIR概念模型确定湖北省农业水土资源承载力评价指标因子，利用熵权TOPSIS评价法对湖北省17个地区(二级行政区)农业水土资源承载力进行评价研究，以期为合理开发利用水土资源、促进水土资源优化配置和农业经济发展提供科学依据。

1 研究区概况与数据来源 1.1 研究区概况

湖北省位于我国中部地区、长江中游，地理位置108°21′ ~ 116°07′E、29°01′ ~ 33°6′N，国土面积18.59万km²，下辖12个省辖市、3个省直管市、1个自治州和1个林区，常住人口为5 851.5万(2015年底)。地势呈东西北三面环山、中间低平、向南敞开的不完整盆地，平原湖区、丘陵岗地和山地各占总面积的31.78%、59.92%、8.30%。截至2015年底全省共有耕地面积343.62×10⁴ hm²，主要分布在江汉平原和鄂北岗地，其中水田面积201.04×10⁴ hm²，旱地面积142.58×10⁴ hm²，人均耕地面积0.08 hm²。湖北省地处北亚热带季风性湿润、半湿润气候区，光热充足，降水丰富，多年平均气温15 ~ 17 ℃，多年平均降水量1 177 mm。2015年全省水资源总量为1 015.63× 10⁸ m³，其中地表水资源量986.35×10⁸ m³，地下水资源量279.64×10⁸ m³。尽管湖北省水资源总量比较丰富，但存在降水变率大、空间分布不均的特点，导致水资源供给极不稳定，旱涝灾害频繁发生^[12]。

1.2 数据来源

本研究采用的农业水资源、耕地资源以及社会经济等方面的基础数据主要来自《湖北省统计年鉴》、《湖北省水资源公报》、《中国水资源公报》等资料，部分指标是利用原始数据计算获得，选取的资料主要包括4个时间断面，分别为2000年、2005年、2010年和2015年。

2 农业水土资源承载力评价方法 2.1 指标体系建立

指标体系的建立是农业水土资源承载力研究的重要依据和基础环节，其指标选择直接关系到评价结果的科学性和合理性。本文采用DPSIR评价指标体系概念模型^{[1, 13-14]}建立湖北省城市尺度的农业水土资源承载力评价指标体系。DPSIR模型将表征一个自然系统的评价指标分成驱动力(driving force)、压力(pressure)、状态(state)、影响(impact)、响应(response)5大要素，其中“驱动力”是决定农业水土资源配置和利用方式的原始自然和社会经济因素；“压力”是在驱动力作用下，反映水土资源系统结构和功能的关键信息；“状态”表征农业水土资源的使用现状和发展水平；“影响”主要是指在驱动力和压力作用下水土资源系统数量和质量发生的变化；“响应”则是人类为提高农业水土资源承载力水平采取的措施。DPSIR模型各要素之间具有明显的连续反馈机制，因此被认为是寻找人类活动及其环境影响之间因果联系的有效途径。在充分理解农业水土资源承载力内涵和水土资源耦合高效利用的基础上，基于数据的可获得性、指标的可操作性、主导性和科学性原则^[15]，选择了17个指标因子建立研究区农业水土资源承载力综合评价指标体系(表 1)。

2.2 熵权TOPSIS评价模型

TOPSIS即“逼近理想解排序方法”(technique for order preference by similarity to ideal solution)，多用来解决有限方案多目标决策问题，是一种在系统分析中常用的运用距离作为评价标准的综合评价方法^[16-18]。其基本思想是通过确定评价单元的正、负理想解(分别用最优向量和最劣向量表示)，分别计算各评价对象与正理想解和负理想解的欧氏距离，获得各评价对象与最优方案的相对接近程度，以此作为评价优劣的依据^[19]。使用TOPSIS模型的中心环节是权重的确定。而熵权法是根据各评价对象提供的原始信息确定其效用价值，是一种客观赋权方法。本文综合运用熵权TOPSIS模型评价湖北省农业水土资源承载力状况。

2.2.1 构建标准化评价矩阵

湖北省农业水土资源承载力问题的原始评价指标矩阵为：

$X = \left[ \begin{array}{l} {x_{11}}\;\;\;\;{x_{12}}\;\;\;\; \cdots \;\;\;\;{x_{1n}}\\ {x_{21}}\;\;\;\;{x_{22}}\;\;\;\; \cdots \;\;\;\;{x_{2n}}\\ \;\; \vdots \;\;\;\;\;\;\; \vdots \;\;\;\;\;\;\;\;\;\;\;\;\;\;\; \vdots \\ {x_{m1}}\;\;\;\;{x_{m2}}\;\;\;\; \cdots \;\;\;{x_{mn}} \end{array} \right]$

(1)

对于高优指标(越大越好)采用公式(2)进行标准化处理，对于低优指标(越小越好)采用公式(3)进行标准化处理，得到标准化评价矩阵(矩阵(4))。

${r_{ij}} = {x_{ij}}/{\rm{max}}\left( {{x_j}} \right)$

(2)

${r_{ij}} = {\rm{min}}\left( {{x_j}} \right)/{x_{ij}}$

(3)

$ R = \left[ \begin{array}{l} {r_{11}}\;\;\;\;{r_{12}}\;\;\;\; \cdots \;\;\;\;{r_{1n}}\\ {r_{21}}\;\;\;\;{r_{22}}\;\;\;\; \cdots \;\;\;\;{r_{2n}}\\ \;\; \vdots \;\;\;\;\;\;\; \vdots \;\;\;\;\;\;\;\;\;\;\;\;\; \vdots \\ {r_{m1}}\;\;\;\;{r_{m2}}\;\;\; \cdots \;\;\;{r_{mn}} \end{array} \right] $

(4)

式中：X为原始评价矩阵；x_ij为第i个评价对象第j个指标的初始值(i = 1, 2, …, m；j = 1, 2, …, n)；max(x_j)和min(x_j)分别代表第j个指标的最大值和最小值；R为标准化后的评价矩阵；r_ij为第i个评价对象第j个指标的标准化值。

2.2.2 确定指标权重

熵权计算公式为：

$ {W_j} = \frac{{1 - {H_j}}}{{m - \sum\limits_j^n {{H_j}} }} $

(5)

$ {H_j} = \frac{1}{{\ln m}}\sum\limits_i^m {{f_{ij}}\ln {f_{ij}}} $

(6)

$ {f_{ij}} = \frac{{{r_{ij}}}}{{\sum\limits_i^m {{r_{ij}}} }} $

(7)

式中：w_j为第j个指标的权重；H_j称为信息熵；f_ij称为指标的特征比重。确定好各指标的熵权之后，即可计算规范化加权评价矩阵。

$ Y = {({y_{ij}})_{m \times n}} = \left[ \begin{array}{l} {r_{11}} \times {w_1}\;\;\;\;\;{r_{12}} \times {w_1}\;\;\;\;\; \cdots \;\;\;\;\;{r_{1n}} \times {w_1}\\ {r_{21}} \times {w_2}\;\;\;\;\;{r_{22}} \times {w_2}\;\;\;\;\; \cdots \;\;\;\;\;{r_{2n}} \times {w_2}\\ \;\;\;\;\; \vdots \;\;\;\;\;\;\;\;\;\;\;\;\;\; \vdots \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; \vdots \\ {r_{m1}} \times {w_n}\;\;\;\;\;{r_{m2}} \times {w_n}\;\;\;\;\; \cdots \;\;\;\;\;{r_{mn}} \times {w_n} \end{array} \right] $

(8)

2.2.3 确定正负理想解

规范化加权评价矩阵各列的最大、最小值构成最优、最劣向量，分别记为：

${Y^ + } = ({y_{{\rm{max}}1}},{y_{{\rm{max}}2}}, \ldots ,{y_{{\rm{max}}n}})$

(9)

${Y^ - } = ({y_{{\rm{min}}1}},{y_{{\rm{min}}2}}, \ldots ,{y_{{\rm{min}}n}})$

(10)

式中：y_maxn，y_minn分别为规范化加权评价矩阵第n列的最大和最小值。

2.2.4 距离计算

采用欧氏距离公式计算第i个评价对象与正、负理想解的距离，计算方法为：

$ D_i^ + = \sqrt {\sum\limits_j^n {{{({y_{\max j}} - {y_{ij}})}^2}} } $

(11)

$ D_i^ - = \sqrt {\sum\limits_j^n {{{({y_{\min j}} - {y_{ij}})}^2}} } $

(12)

式中：y_ij为第i个评价对象第j个指标加权后的标准化值。

2.2.5 计算综合评价指数

研究区第i个评价对象农业水土资源承载力与最优承载力的贴进度即表示实际承载力的大小(C_i)，其值介于[0, 1]，C_i越大表明农业水土资源承载力水平越高，其计算公式为：

$ {C_i} = \frac{{D_i^ - }}{{D_i^ + + D_i^ - }} $

(13)

3 结果与分析 3.1 农业水土资源空间分布及匹配

区域水资源的丰缺可用单位土地面积拥有的水资源量表示，土地资源的开垦程度可用耕地面积占国土总面积的比例表示，而水土资源匹配程度用单位面积耕地可拥有的水资源量(水土资源匹配系数)来反映。以湖北省17个城市作为基本评价单元，以多年平均单位面积水资源量和耕地资源比例作为度量指标，分析湖北省农业水土资源空间分布特征及其丰缺状况(图 1)。在此基础上，参考刘彦随等^[20]提出的农业水土资源匹配系数测算方法，计算了各地区水土资源匹配系数，根据水土资源匹配程度利用ArcGIS自然断裂法将其划分为4个等级并实现空间可视化(图 2)。

湖北省多年平均单位面积水资源量为51.28× 10⁴ m³/km²，高于全国均值水平(30.23×10⁴ m³/km²)，但其空间分布不均衡，总体上南部地区单位面积水资源量多于北部地区，单位面积水资源量最多的5个地区依次是咸宁市、黄石市、鄂州市、恩施自治州和黄冈市，集中分布于湖北省东南部或西南部的山地、丘陵地区。而单位面积水资源量较少的荆门市、随州市、襄阳市、十堰市和孝感市集中分布于中部或北部，这与湖北省多年平均年降水量的分布特征具有较强的一致性。耕地资源的分布明显受到地形的影响，全省耕地资源比例超过30%的城市有4个，由高到低分别是潜江市、天门市、仙桃市和荆州市，均分布于江汉平原腹地，而西部地区的神农架林区、十堰市、宜昌市和恩施自治州以山地为主，地形起伏大，耕地资源比例较小。

农业水资源和耕地资源空间分布不均导致研究区水土资源匹配程度表现出较大的异质性，总体为“东西优，中部差”的匹配格局。水土资源匹配系数低于全省平均水平(1.45×10⁴ m³/hm²)的城市有10个，全部位于湖北省中部，水土资源匹配相对较差的原因是该地区地形以平原和丘陵岗地为主，土地开垦程度较高，土壤肥沃，优质耕地资源比较丰富，但就全省而言，其水资源相对短缺，水资源状况无法满足农业生产的需要。以2015年为例，中部10市耕地面积占全省的61.67%，而其水资源总量仅占全省的38.06%，因此单位面积耕地可拥有的水资源量较少。相对中部地区而言，东部的咸宁、黄石、黄冈和西部的十堰、宜昌、恩施和神农架林区大部分地区为降水高值区，水资源比较充足，但土地垦殖率较低，因此水土资源匹配程度比较均衡，水资源基础条件完全可以满足农业生产的需要。

3.2 农业水土资源承载力测算及时空分布格局 3.2.1 水土资源承载力测算及等级划分

根据上文构建的湖北省农业水土资源承载力评价指标体系，利用熵权法确定各评价指标的权重(表 2)，结合TOPSIS方法求解农业水土资源承载力综合评价指数，根据数值大小进行降序排序，其结果如表 3所示。利用熵值法确定权重时，某项指标的正、负理想解之间的差距越大，其熵值越小，说明该指标能提供更多的有用信息，因此被赋予的权重越大。由表 2可以看出，在17项评价指标中，人均水资源量(X₉)和水土资源匹配系数(X₁₁)平均权重都大于0.19，说明该两项指标解释了不同地区农业水土资源承载力差异的较多信息，是影响农业水土资源承载力大小的最主要因素；其次为农业水资源利用效率(X₁₄)，其权重大于0.1，说明农业水土资源承载力大小不仅取决于资源禀赋条件，还受到资源利用技术水平的较大影响。

从湖北省农业水土资源承载力综合评价指数(C_i)的计算结果可以看出，综合评价指数最高的是神农架林区，其平均值为0.786，远远高于其他地区的承载水平，并且其排名一直高居第1位；综合评价指数最低的为十堰市，其多年平均值仅为0.270，与神农架林区相差0.516。利用ArcGIS的探索数据功能计算C_i的变异系数，结果显示其变异系数为23.65% ~ 35.75%，属于中等空间变异性，表明不同地区水土资源承载力水平存在较大差异。从排名动态变化来看，水土资源承载力水平进步最大的是荆州市，其排名从2000年的第9名上升到2015年的第3名；与此同时，排名下降最多的是孝感市，由2000年的第7名下降到2015年的第13名。其他地区排名相对位置变化不大，一般在3个次序之内，说明湖北省农业水土资源承载力的相对水平具有一定的稳定性。

由于农业水土资源承载力影响因子相互作用的复杂性及其他不确定性因素的存在，综合评价指数的绝对数值并非与真实承载能力呈线性对应关系，因此有必要根据其集聚和离散的分异特征进行等级划分。根据湖北省农业水土资源禀赋条件和社会经济发展的阶段特点，并考虑科学技术进步和资源配置等因素对水土资源承载力可拓性的影响，将湖北省农业水土资源承载力定性地划分为4个等级：盈余承载(Ⅰ级)、安全承载(Ⅱ级)、临界承载(Ⅲ级)和超载(Ⅳ级)。

3.2.2 水土资源承载力时空分布格局

在GIS技术的支持下实现湖北省农业水土资源承载力等级时空分布格局的可视化表达(图 3)。由图 3可以看出，湖北省农业水土资源承载力水平存在较大的空间异质性。整体上，农业水土资源承载力最高的是神农架林区和恩施自治州，其承载力水平大部分时段处于盈余状态，表明两地区农业水土资源完全可以满足发展需要，未来进一步开发利用的潜力很大。神农架林区和恩施自治州位于湖北省西南部，地形以中高海拔的山地为主，森林覆盖率很高，年均降水1 100 mm以上，水资源十分丰富，虽耕地资源有限，但由于人口密度全省最低，社会经济活动对生态环境破坏较小，因此水土资源利用与社会经济发展非常协调。其次是荆门、仙桃、潜江、天门、武汉、荆州5个城市，水土资源承载力基本保持在Ⅱ级以上水平(安全承载为主)，水土资源利用和社会经济发展比较均衡，水土资源未来尚有较大的开发潜力。不难发现，这些地区刚好处于江汉平原的核心地带，自然条件优越，耕地资源充足，是湖北省最重要的粮食产区，只要对水土资源科学开发、有效利用，水土资源承载力将会继续维持在较高水平。处于第三梯队的是随州、襄阳、孝感、宜昌、黄冈等城市，水土资源承载力以临界超载为主，表明其社会经济快速发展与水土资源供给不足的矛盾开始凸显，水土资源面临的形势不容乐观。虽然襄阳、随州、孝感3市耕地资源比较丰富，但水资源相对短缺，其多年平均水土资源匹配系数分别为0.74×10⁴、1.05×10⁴、0.73×10⁴ m³/hm²，远远低于湖北省水土资源匹配系数的均值水平(1.45×10⁴ m³/hm²)；相反，宜昌和黄冈两地水资源较为丰富，但耕地面积相对不足。农业水土资源承载力最低的是十堰、鄂州、咸宁、黄石4个城市，2000—2015年一直呈超载状态，说明其社会经济发展与水土资源供给不平衡，现状条件下水土资源进一步开发利用的潜力较小。原因主要是这些地区以山地丘陵为主，耕地资源短缺，水土流失严重；工业、高耗水产业比重较大，对生态环境造成很大破坏。此外，非农产业挤占农业用地、水资源利用效率较低也加剧了上述地区水土资源紧张态势。

从农业水土资源承载力水平的时间变化趋势来看，虽然大部分地区综合承载力指数呈现波动，但其承载力等级较为稳定，仅荆州市实现了承载力状态的较大跨越。荆州市2000年和2005年的水土资源承载力为临界超载，2010年变为安全承载，2015年为盈余承载。原因可能是该市大力加强农田水利建设，保护耕地，发展节水节能技术，提高农业现代化水平，从而促进了农业水土资源承载力的较大提升。

4 讨论

各地区应根据农业水土资源承载力的现状采取不同的应对措施。以盈余承载为主的地区在满足经济社会发展对水土资源需求的前提下，应更加注重生态环境的保护，适度开发当地特色资源，促进居民增产增收，建设生态宜居城市。以安全承载为主的地区，应在保护农业水资源和耕地资源的前提下，增加农业资金和技术投入，提升农业现代化水平，挖掘粮食增产潜力；加强农田基本建设，大力推广节水技术，提高水资源利用效率。以临界超载和超载为主的区域，必须严格控制非农产业对耕地的占用，保护现有耕地数量和生产能力，预防潜在的水土流失；限制高耗水、高污染产业的发展，加强区域间水资源配置工程建设；特别注意优化产业结构，因地制宜地发展当地特色产业，促进水土资源可持续利用和社会经济可持续发展。

熵权TOPSIS将TOPSIS模型的简单实用和熵值法的客观赋权相结合，是进行农业水土资源承载力综合研究的一种可行方法，采用该方法评价湖北省农业水土资源承载力时空动态，其结论与实际情况相符，可为研究区农业水土资源可持续高效利用提供决策参考。但利用该模型只能测度区域各评价单元水土资源承载力的相对等级，无法反映真实的承载能力大小。因此，今后工作重点应是构建更为科学合理的、具有广泛适用性的水土资源承载力评价指标体系，并制定切实可行的评价指标等级划分标准，以加强不同区域农业水土资源承载力的对比分析和综合研究。

5 结论

湖北省农业水资源和耕地资源空间分布不均衡，多年平均水土资源匹配系数为1.45×10⁴ m³/hm²，总体呈现“东西优，中部差”的匹配格局。不同地区水土资源承载力水平存在较大差异，其中，神农架林区和恩施自治州最高，承载力水平基本处于盈余状态；其次是荆门、仙桃、潜江、天门、武汉、荆州5个城市，水土资源承载力基本以安全承载为主；再次是随州、襄阳、孝感、宜昌、黄冈等城市，水土资源承载力以临界超载为主；水土资源承载力最低的是十堰、鄂州、咸宁、黄石4个城市，2000—2015年一直呈超载状态。从农业水土资源承载力水平时间变化趋势来看，大部分地区综合承载力指数出现波动，但其承载力等级较为稳定，只有荆州市水土资源承载力水平出现显著提升。