基于资源环境承载力和国土空间开发适宜性双评价的农业土地资源评价操作步骤

本

系列基于自然资源部2019.3月

和2020年1月印发的两版《资源环境承载能力和国土空间开发适宜性评价指南（试行）》（简称“双评价”），因为计算公式与方法在2020.1月版的少之又少，可以说不限制模型方法，只有规定大体的方向，但是都可以基于2019.6月的指标、模型、方法进行计算，因此一些定义、目标、流程、制图规范都基于新版的，但是方法都基于旧版的。

首先说明基于2020.1月版的“双评价”的定义、目标：

定义：（1）“资源环境承载力评价”：在一定发展阶段，经济技术水平和生产生活方式，一定地域范围内资源环境要素能够支撑的农业生产、城镇建设等人类活动的最大规模。

（2）“国土空间开发适宜性评价”：在维系生态系统健康的前提下，综合考虑资源环境要素和区位条件以及特定国土空间，进行农业生产城镇建设等人类活动的适宜程度。

目标：分析区域资源禀赋与环境条件，研判国土空间开发利用问题和风险，识别生态保护极重要区（含生态系统服务功能极重要区和生态极脆弱区），明确农业生产、城镇建设的最大合理规模和适宜空间，为编制国土空间规划，优化国土空间开发保护格局，完善区域主题功能定位，划定三条控制线，实施国土空间生态修复和国土综合整治重大工程提供基础依据，促成形成以生态优先、绿色发展为导向的高质量发展新路子。

然后说明于2020.1月版的“双评价”的工作流程以及要求：

工作流程：

通过工作流程可看出，本系列将分为五大步：

第一步，首先要进行本底评价，也就是以往版本双评价中的1土地资源评价、2水资源评价、3气候评价、4生态评价、5环境评价、6灾害评价六个单项评价。需要提出的是，虽然新版的双评价似乎可以直接开始生态保护重要性评价，但是本系列依然先从6个本底评价开始，各位可以当作一个旧版的研究评价。

第二步：生态保护重要性评价。生态保护省级评价从区域生态安全底线出发，在陆海全域，评价水源涵养、水土保持、生物多样性维护、防风固沙、海岸防护等生态系统服务功能重要性，以及水土流失、石漠化、土地沙化、海岸侵蚀及沙源流失等生态敏感性，综合形成极重要区和重要区。

第三步：农业生产适宜性评价。省级在生态保护极重要区以外的区域，开展种植业、畜牧业、渔业等农业生产适宜性评价，识别农业生产适宜区和不适宜区。

城镇建设适宜性评价。省级在生态保护极重要区以外的区域，优先考虑环境安全、粮食安全和地质安全等底线要求，识别城镇建设不适宜区。沿海地区针对海洋开发利用活动开展评价。

第四步：承载规模评价。基于现有经济技术水平和生产生活方式，以水资源、空间约束等为主要约束，缺水地区重点考虑水平衡，分别评价各评价单元可承载农业生产、城镇建设的最大合理规模。各地可结合环境质量目标、污染物排放标准和总量控制等因素，评价环境容量对农业生产、城镇建设约束要求。按照短板原理，取各约束条件下的最小值作为可承载的最大合理规模。

第五步：综合分析，如资源禀赋分析、问题和风险识别、潜力分析和情景分析等。

最后说明于2020.1月版的“双评价”的制图规范：

要求：评价统一使用2000国家大地坐标系（CGCS2000），高斯-克吕格投影，陆域部分采用1985国家高程基准，海域部分采用理论深度基准面高程基准。制图规范、精度等参考省级国土空间规划要求。

制图规范：

那么再参考《省级国土空间规划编制指南》，明确制图规范与精度的要求：

字体：汉字宋体，西字Times New Roma。

比例尺：1:50万或1:100万。

指北针：。

图例：图幅内左下角或右下角。

因为本系列只是起研究模拟作用，因此以广东省（省级）为例，地理坐标系使用WGS1984（许多网上得到的数据都是WGS84，所以不使用CGCS2000），投影使用WGS_1984_UTM_Zone_49N，空间分辨率设为90m。

首先，对广东省水资源、土地资源、气候、生态、环境、灾害进行评价。

那么我们先从土地资源评价开始，这样第一步就开始了！

1-1-1 广东省农业土地资源评价

公式：农业指向的土地资源评价：[农业耕作条件]＝f ([坡度], [土壤质地])

农业指向的步骤要求（坐标系与精度与双评价要求不同，主要操作步骤相同）：

坡度：利用DEM，计算地形坡度，按≤2°、2~6°、6~15°、15~25°、>25°划分为平地、平坡地、缓坡地、缓陡坡地、陡坡地5 个等级，生成坡度分级图。

土壤质地：将土壤的粉砂含量≥80%区域，农业土地资源直接取最低等；60%≤粉砂土含量＜80%的区域，将坡度分级降 1 级作为农业土地资源等级。

首先我们需要的数据：DEM数据、土壤质地数据。

广东省DEM数据如下，空间分辨率为90m：

首先坡度分析，得到slope：

重分类，（至于为什么取50,40,30,20,10，在土壤质地分类后会讲到）得到nyslope，那么50,40,30,20,10依次代表平地、平坡地、缓坡地、缓陡坡地、陡坡地：

接下来需要根据土壤质地进行分类，那么在2020.1月版的双评价中，没有指明划分标准，因此参照2019.6月版的：

可以看出主要是从粉砂的含量进行划分，那么需要使用到HWSD世界土壤数据库，然后使用到其中的clay字段，clay指黏土，那么粗略的可以认为，土壤中若粉砂≥80%，则黏土≤20%，若粉砂≥60%，则黏土≤40%，因此，将土壤质地按照黏土含量进行分级，≤20%，取最低，20%

首先按掩膜提取广东省域的，然后需要进行投影以及重采样至90m分辨率，再次不重复类似操作：

然后进行连接，连接HWSD_DA TA的MU_GLOBAL：

HWSD_DATA中的CLAY就是我们所需的黏土含量，然后我们新建字段，赋值后移除连接：

然后我们就可以根据得到的CLAY字段进行重分类了：