基于生物多样性保护的耕地入选基本农田布局

王建英，李江风，，邹利林

（1.中国地质大学（武汉）经济管理学院，武汉430074；2.中国地质大学（武汉）公共管理学院，武汉430074）

随着我国人口的持续增长和社会经济的快速发展，人类对耕地资源的压力和威胁与日剧增，耕地资源数量急剧下降和质量明显减退的现状引起了政府和学术界的广泛关注［1］，将优质耕地划为永久基本农田，实施特殊保护成为保障国家粮食安全和保护耕地的一种有效措施。我国耕地入选基本农田的相关研究已较为广泛，早期的研究多从耕地自然质量评价出发［2－3］，筛选出土壤肥沃、分布集中连片等高质量耕地划为基本农田保护区。然而这种基本农田划定思路缺少对耕地立地条件的考虑，造成基本农田被频繁调入和调出，不具有稳定性和长久性［2，4］。因此，考虑社会经济发展中建设占用耕地的需求，部分学者通过综合考虑耕地自然质量和社会经济影响因素，构建基本农田决策模型，并依据耕地综合质量排序结果划定基本农田［5－8］，以权衡社会经济发展需求与保障基本农田质量之间的矛盾。目前的相关研究中较少考虑基本农田，特别是永久性基本农田划定对区域生物多样性保护、生态流水平移动等景观生态安全格局的影响，特别是在生物多样性丰富地区，基本农田保护区内化肥农药的不当使用产生土壤污染、引起生物链断裂，对区域内生物多样性保护产生很大的负面作用，不利于协调区域人口、资源和环境之间的关系［9］。针对这一问题，本文以乐东黎族自治县为例，考虑研究区社会经济发展和生物多样性保护对基本农田保护的影响，构建耕地入选基本农田的决策模型，采用熵权理想点逼近法对耕地质量综合评价结果进行排序，筛选出符合入选条件的耕地划为基本农田，为新一轮土地利用总体规划编制中基本农田的科学划定提供思路借鉴。

1 研究区概况

乐东黎族自治县位于海南岛西南部，北纬18°24′—18°58′，东经108°39′—109°24′，属于典型的干湿季交替热带区域。全县东西长72km，南北宽58km，共有建制镇11个，2010年总人口为525 830人。据全国第二次土地调查结果显示：乐东黎族土地自治县总面积为2 763.22km2，耕地面积47 568.12 hm2，占土地总面积的17.20%。其中，水田面积为29 553.51hm2，水浇地面积为329.26hm2，旱地面积为17 685.34hm2，分别占耕地总面积的62.13%，0.69%，37.18%。人均耕地面积0.091hm2，略低于全国平均水平0.091 3hm2。

研究区森林覆盖率高达70%，区内尖峰岭国家自然保护区几乎涵盖了世界热带地区所有的植被类型，物种资源丰富，包括海南苏铁、水松、水杉、毛坡垒、望天树、金丝李、野生兰和野生稻等在内的2 835种国家级和省级重点保护植物和包括黑长臂猿、云豹、猕猴、黑熊、白鹭、黄鹂、画眉、八声杜鹃、原鸡、白鹇、孔雀雉、蛇雕等在内的2 300多种国家重点保护动物和珍稀动物。丰富多样的物种资源，不仅具有极高的听觉、视觉和观赏价值，还具有巨大的生态价值。

2 数据来源及研究思路

2.1 数据来源

研究使用的基础数据包括乐东黎族自治县耕地自然等级图，第二次农村土地利用现状调查图及数据（2009年），乐东黎族自治县土地利用总体规划（2006—2020年）文本和相关图件等，均来源于《海南省乐东黎族自治县土地利用总体规划（2006—2020）年》项目和《乐东黎族自治县基本农田划定（补划）》项目。

2.2 研究思路和数据处理

2.2.1 研究思路 耕地入选基本农田实质上是一个复杂的决策问题，需要综合考虑多方面因素的影响［5］，如耕地自身的条件（土壤肥力、区位条件、坡度等）、人为因素的影响（政策限制等）、社会经济发展需求（建设占用等），生态安全（生物多样性保护、景观格局等）。本文以海南省基本农田保护规划的相关理论和政策为指导，结合实地野外调研的结果，提出优先生物多样性保护的耕地入选基本农田的决策思路。首先将土地利用总体规划中建设占用和生物多样性极丰富区的耕地剔除；其次耕地综合质量评价体系是在总结前人研究成果的基础上［6，10－13］，构建耕地综合质量评价体系对耕地进行综合质量排序，从耕地自然质量、区位条件、水利设施条件、城镇村镇辐射影响、生态安全5个方面选取8个评价决策指标进行评价，运用层次分析法，对耕地入选基本农田评价指标进行权重确定（表1），以GIS软件为平台，采用熵权理想点逼近排序法确定耕地综合质量等级的高低，最终划定基本农田保护区的范围。具体技术路线如图1所示。

表1 耕地入选基本农田评价决策指标

图1 技术路线图

2.2.2 数据处理 数据处理主要采用ArcGIS软件，具体处理过程如下：（1）从乐东黎族自治县第二次农村土地利用现状调查图（2009年）中提取出耕地利用现状图（包括水田、水浇地和旱地），保存为shp格式，面积合计为47 568.12hm2；（2）将耕地自然质量等级图从北京54坐标校正为与耕地利用现状图一致的西安80坐标，并对两者进行空间分析，获得每一个耕地图斑对应的耕地自然质量等指数；（3）在3D Analyst模块中，利用等高线创建TIN文件，获得研究区坡度图，与耕地利用现状图进行相交分析得到每块耕地图斑对应的坡度；（4）将耕地以地理图斑为单元进行Buffer分析，设置缓冲半径参数为20m，得到耕地缓冲区图；将破碎的耕地缓冲区采用Dissolve工具进行融合；融合后的缓冲区与现状耕地进行相交分析，得到每个耕地图斑对应的缓冲区面积。耕地连片面积采用缓冲区面积表示；（5）利用Euclidean Distance工具对从乐东黎族自治县第二次农村土地利用现状调查图（2009年）中提取的水域、交通用地、城镇、农村居民点进行缓冲区分析，得到县水域缓冲区分析图；运用Zonal Statistics工具计算每个耕地图斑中心点到水域、交通用地、城镇、农村居民点中心点的最小距离；（6）将生物多样性丰富度分区图与耕地图斑进行相交分析，得到每个耕地图斑上的生物多样性丰富度。

3 研究过程

3.1 调出预留耕地

考虑区域经济发展中的建设用地、生物多样性保护及地方其他政策性因素占用基本农田，故在划定基本农田前将可能被占用的耕地调出。调出的耕地主要包括：（1）土地利用总体规划（2006—2020年）中新增建设用地和城乡建设用地扩展范围内的耕地；（2）国家和省重点建设项目规划范围内的耕地，已列入国家或者本省城镇建设规划近期建设控制区内的耕地，已列入国家或者本省规划的非农业开发区近期建设控制区内的耕地；（3）生物多样性丰富度极为丰富区内的耕地。第二次农村土地利用现状调查的耕地面积为47 568.12hm2，建设占用预留耕地544.89hm2，生物多样性保护占用耕地101.54hm2，剩余的耕地面积为46 921.69hm2，可以满足土地利用总体规划确定的基本农田保护区规模不得低于46 755hm2的要求。

3.2 构建基本农田决策模型

3.2.1 构造规范化的决策矩阵 乐东黎族自治县内14 857块耕地组成决策对象集A＝｛a1，a2，…，an，n＝14 857｝，根据实际情况以及耕地入选基本农田的决策指标体系，确定一组决策指标F＝｛f1，f2，…，fm，m＝8｝。决策对象ai（i＝1，2，…，m）在决策指标fi（i＝1，2，…，m）下的取值为xij，由xij组成决策矩阵：

由前面层次分析法确定的指标权重表示为：

W ＝ ［w1，w2，…，wm］T＝ ［0.173，0.1213，0.0825，0.2562］［0.0479，0.0731，0.1347，0.112］T

矩阵X乘以权重，决策矩阵X就转化为一个规范化的决策矩阵R。

式中ri＝（ri1，ri2，…，rim）（i＝1，2，…，n）可看作欧式空间里的一个“决策点”。

3.2.2 理想点法排序

（1）确定理想点和反理想点。一般来说，决策指标可分为收益型指标和损失型指标。收益型指标数值越大对于评估结果越有利，如耕地自然质量等指数；损失型指标数值越大对于评估结果越有害，如耕地到城镇用地的距离。用M1表示理想点，M2表示反理想点。当指标为收益型时，理想点为矩阵R的列向量的最大值，反理想点则为最小值；当指标为损失型时，理想点为矩阵R的列向量的最小值，反理想点为最 大 值。即：M1＝ （p1′，p2′，…，pm′），M2＝ （p1″，p2″，…，pm″）

（2）理想点评价函数。即计算评价单元到理想点M1和反理想点M2的“距离”，第i块耕地i（1，2，…，n）到M1和M2两点的“距离”，分别用S1i和S2i表示。

S1i值越小，说明评价单元距离理想点越近，耕地入选基本农田的优先程度越大；S2i值越小，说明评价单元距离反理想点越近，耕地入选基本农田的优先程度越小；因此，评价单元的S1i值越小，S2i值越大，即为划定基本农田的优选耕地。

（3）计算理想点的贴近度。

显然Ti∈［0，1］，用指标Ti进行评价时，Ti值最大的规划方案即为最佳规划方案。

（4）综合排序。计算14 857个耕地对应的贴近度并进行排序，Ti值大的耕地入选基本农田的优先度越高。若两块耕地的贴近度相同，则以距理想点的“距离”值较小者为优选耕地。

4 结果与分析

乐东黎族自治县耕地入选基本农田综合质量评价选择耕地自然质量等指数、耕地坡度、耕地连片程度、灌溉保证率、城镇辐射影响和村镇辐射影响、距道路距离、生物多样性丰富度8个评价指标，运用层次分析法确定权重，采用熵权理想点逼近法对14 857块耕地质量进行了综合评价，并根据贴近度的大小对耕地入选基本农田的优先程度排序。根据贴近度将耕地图斑分为5个区间，即：最优区间、较优区间、一般区间、较差区间和最差区 间，其贴近度的范围依次为（0.640 5，0.916 8］，（0.526 6，0.640 5］，（0.429 5，0.526 6］，（0.330 2，0.429 5］，［0.0000，0.330 2］，分别占耕地总面积的 12.60%，33.33%，29.90%，10.73%，13.45%。其中较优区间和一般区间的耕地面积最多，在各乡镇均有分布；最优区间的耕地主要分布在滨海佛罗镇和尖峰镇，最差区间的耕地主要分布在抱由镇、利国镇和千家镇（表2）。

按照土地利用总体规划（2006—2020年）中确定的基本农田保护区规模，依据各个耕地图斑的理想点贴近度进行排序，将排序结果中名次靠前且总面积之和大于等于46 755hm2的耕地划为基本农田，共有14 573块耕地划入基本农田，合计46 756.89hm2。

表2 评价结果hm2

耕地入选基本农田综合评价是采用熵权逼近理想点法计算贴近度的大小来衡量的。评价结果可以作为耕地入选基本农田的依据，但这必须建立在评价结果真实可靠的基础上，因此有必要对评价结果的准确合理性进行分析。

（1）贴近度最优区间（0.640 5，0.916 8］内无自然质量等级低的五等地和四等地，而与质量等别最高的一等耕地耦合度最高，约有70%左右的一等耕地划入最优区间，且集中分布在西南滨海地区和东部内陆山区，剩余的一等地耕地都划入了较优区间和一般区间；而在最差区间中，自然质量等级较差的五等耕地占该区间耕地总量的30%，四等耕地占该区间耕地总量的63%，主要分布在西部滨海的尖峰镇和佛罗镇，无自然质量高的一等地和二等地。基本农田综合排序结果与耕地自然质量等指标数值具有较强的一致性。

（2）研究区耕地多分布在坡度极为平缓的0～3°和3°～5°区，故贴近度分区中均以0～3°和3°～5°的平地和缓坡地为主，在分布于西南滨海乡镇的最优区间（0.640 5，0.916 8］中，88%左右的耕地坡度在0～3°之间，而最差区间中0～3°之间坡度比例仅占60.45%，说明排序结果与耕地坡度级别基本一致，与质量等别高的耕地分布也具有一致性的，排序结果十分合理。

（3）最优区间中耕地到河流水域的距离在1 000 m以内的耕地数量最多，约占耕地总量的67%，且该区间内没有距离超过4 000m的耕地；在较优区间中到河流水域距离在1 000m以内的耕地则占70%左右，超过4 000m的耕地仅占1.4%；在较差区间中耕地到河流水域的距离在6 000m以外的耕地数量较多，约占区间耕地总量的74.5%，而距离在1 000m以内的仅占11.2%；在最差区间中耕地到水域的距离在6 000m以外的占区间耕地最多，约占该区间耕地总量的99.85%，没有距离在1 000m以内的耕地。通过排序结果与水利基础设施水平指标的比较，说明排序结果与水利基础设施水平指标数值有很强的一致性。

（4）最优区间中耕地到道路用地的距离在1 000 m以内的耕地数量最多，占区间耕地总量的67%，主要分布在各建制镇的周边地区，距离在1 000～6 000 m以内的耕地占该区间耕地总量的33%，没有距离在6 000m以外的耕地；在较优区间中，距离在1 000 m以内的耕地占区间耕地总量的70.82%，距离在6 000m以外的耕地仅占1.4%；在最差区间中，距离在6 000m以外的耕地占区间耕地总量的99.85%，没有距离在1 000m以内的耕地；在较差区间中，距离在6 000m以外的耕地占区间耕地总量的74.49%，距离在1 000m以内的耕地仅占11.20%。通过比较分析，说明在交通区位条件指标的检验下，排序结果准确合理。

（5）生物多样性极丰富区内的耕地已经剔除，最优区间中基本无生物多样性丰富区，以生物多样性不丰富区为主，主要分布在滨海的尖峰镇、佛罗镇。生物多样性较丰富区与耕地分布的空间耦合性较低，约为172hm2，多分布在一般区间中。耕地入选基本农田的最差区间内，以生物多样性不丰富区为主，有少量的生物多样性一般区。通过比较分析，说明耕地综合质量排序受到生物多样性丰富度这一因子的影响，排序结果准确合理。

通过将耕地入选基本农田综合排序结果与耕地自然质量等级、耕地坡度、水利基础设施水平、交通区位条件、生物多样性丰富度5个评价指标进行比较分析，说明乐东县2009年土地利用现状耕地的综合质量评价结果准确可靠，可以作为科学合理划定基本农田的依据。

5 结论与讨论

（1）耕地入选基本农田的相关研究已较丰富，依据耕地质量评价结果划定基本农田已成为学者的共识［14］。在耕地质量评价过程中，本研究不仅考虑了建设预留耕地的需要，还加入了研究区生物多样性保护对基本农田布局的要求，避免基本农田划定及后续基本农田整治工程对区域物种生境保护和生物多样性维持产生负面影响，为生物多样性丰富区的基本农田划定提供了参考。

（2）应用GIS技术和熵权理想点排序法相结合的基本农田划定能够客观地反映每块耕地的综合质量，并依据上级下达的任务及时调整基本农田的数量和空间分布，为划定永久性基本农田提供依据，为区域耕地保护和粮食安全提供保障。

（3）乐东县现存耕地面积47 568.12hm2（调出不适宜的耕地之前），而海南省下达的基本农田保护区面积高达46 755hm2，占耕地总面积的98.29%，这一情况直接导致几乎所有的耕地都会调入基本农田。然而，综合考虑社会经济发展和生物多样性保护的要求的耕地适宜性评价结果表明，约有13.45%的耕地综合质量较差，不适宜划为基本农田，无法“保质”完成上级下达的指标，应重新协商避免劣质耕地划入基本农田。基于生物多样性保护的基本农田决策模型不仅为基本农田客观合理划定提供依据，还可以随时依据耕地的实际质量状况，对已划定为基本农田的耕地的空间位置、数量实施有效的动态监测。

［1］ 于亚男，郭晓辉，王秀茹.基本农田整理与我国粮食安全的浅析：以河北省香河县刘宋镇基本农田整理项目为例［J］.水土保持研究，2007，14（3）：318－320，323.

［2］ 杨树佳，郑新奇，王爱萍，等.耕地保护与基本农田布局方法研究：以济南市为例［J］.水土保持研究，2007，14（2）：4－7.

［3］ 王加恩，蔡子华，胡艳华，等.基于土壤生态环境的县域基本农田保护区划：以浙江省嘉善县为例［J］.中国农业资源与区划，2011，32（6）：85－89.

［4］ 袁枫朝，严金明，燕新程.GIS支持下的大都市郊区基本农田空间优化［J］.农业工程学报，2008，32（S1）：61－65，257.

［5］ 程锋，石英，朱德举.耕地入选基本农田决策模型研究［J］.地理与地理信息科学，2003，19（3）：50－53.

［6］ 吴飞，濮励杰，许艳，等.耕地入选基本农田评价与决策［J］.农业工程学报，2009，25（12）：270－277.

［7］ 郑新奇，杨树佳，象伟宁，等.基于农用地分等的基本农田保护空间规划方法研究［J］.农业工程学报，2007，23（1）：66－71，292.

［8］ 钱凤魁，王秋兵.基于农用地分等和LESA方法的基本农田划定［J］.水土保持研究，2011，18（2）：251－255.

［9］ 赵亚莉，吴群.基本农田保护研究综述［J］.国土资源科技管理，2007（6）：30－34.

［10］ 宇向东，郝晋珉，鲍文东.基于耕地分等的基本农田空间配置的方法［J］.农业工程学报，2008，24（S1）：185－189.

［11］ 张落成，王辰，李青.泰州市基本农田空间配置探讨［J］.长江流域资源与环境，2010，19（12）：1421－1426.

［12］ 陈茜，段建南，孔祥斌，等.北京市基本农田保护区内耕地数量提升潜力研究［J］.水土保持研究，2012，19（3）：200－203.

［13］ 黄丽梅，刘廷祥，鲍文东，等.基本农田空间配置的自动化研究［J］.水土保持研究，2008，15（4）：112－115.

［14］ 张蕾娜.农用地分等定级成果在基本农田保护中的应用［J］.经济地理，2006，26（S1）：128－130.