1. 引言

土壤作为陆地生物养分之源极其宝贵，但其形成极为缓慢。表土是指最表层的土壤，作为土壤中养分含量最高的部分，能有利于快速恢复地力和植物生长，珍贵程度更甚[1]-[3]。表土作为资源已成为共识，生产建设项目在进行主体工程设计及水土保持方案报告书编制时，表土分布范围、剥离及利用等相关内容设计已成为水土保持设计的重中之重[4]。然而，目前生产建设项目表土可剥离量主要通过实地勘察，并结合三调结果匡算而来，准确性低。部分项目由于自身情况的特殊性，要求仅对某一坡度以下扰动区域进行表土剥离，若采用传统方法，需根据地形图逐一测量估算，工作量巨大，且准确性无法保证，此时ArcGIS的空间数据处理与分析能力能够完美解决此类问题。

本文基于ArcGIS对某大型灌区生产建设项目的表土剥离与利用进行平衡分析。由于该项目区地形高低起伏较大，坡度大于25˚的林地、草地、耕地表土层较薄且含有大量石块，无法使用机械进行剥离，采用人工剥离成本较高，因此该区域表土不进行剥离。利用ArcGIS的空间数据处理能力分析得出项目区表土量、可剥离表土量，并进行后期利用平衡分析。

2. 应用思路

利用ArcGIS的数据管理、数据编辑与处理、坐标系转换、矢量数据空间分析，可得出项目扰动范围内表土分布情况、分布数量。通过对项目区地形图分析得出DEM文件，经栅格数据处理与分析，将项目扰动区域生成以坡度25˚为界限的适量数据，在此基础上与三调图相交，即可得出项目区可剥离表土范围。将以上处理成果的属性表输出，统计项目区表土量、可剥离表土量和表土需求量，进行后续表土平衡分析。在ArcGIS中亦可制作输出表土分布范围图、可剥离表土分布图。

3. 应用过程

3.1. 项目区所含表土量分析

1) 将项目扰动范围划分为永久占地和临时占地两个图层

由于项目施工永久占地范围与临时占地范围有重合部分，为防止重复计算，需将两者进行擦除分析，得出“施工永久占地”和“擦除后的施工临时占地”两个图层。操作流程如下：将两者cad格式文件导入ArcGIS中，用分析工具/叠加分析/擦除功能得出“擦除后的施工临时占地”。注意赋予两图层与项目区三调图相同的坐标系，否则无法进行后续相交分析。具体方法为：数据管理工具/投影和变换/栅格/定义投影。

2) 计算项目扰动范围内永久占地和临时占地范围内所含表土量

在ArcGIS中导入项目区三调图，分别于“擦除后的施工临时占地”和“施工永久占地”进行相交得出项目扰动范围各地类面积。操作流程如下：分析工具/叠加分析/相交；计算各图斑面积后导出属性表；操作步骤如下：右键打开属性表/添加字段“area”(浮点型)/在“area”右键计算几何；随后导出属性表。同理得出永久占地地类表。

根据现场调查，表土主要分布在耕地、园地、林地、草地等占地类型区。表土厚度依据地类不同差别较大，耕地表土厚度为30 cm，园地表土厚度约35 cm，乔木林地表土厚度约25 cm，灌木林地、竹林地表土厚度约20 cm，其他林地表土厚度约15 cm，草地表土厚度约20 cm。结合现场调查结果可分析出工程区所含表土量(见附表1)。

3.2. 项目区可剥离表土量分析

3.2.1. 项目区坡度分析

利用项目区地形图生成DEM文件，分别得出项目区>25˚和≤25˚的区域。具体操作步骤如下：

1) CAD图层分离：将项目区地形图中DGX和GCD图层复制到新的CAD文件中，将dwg文件转化为可编辑的要素类文件(shapefiles)：在Polyline文件右键数据/导出数据/导出成shp文件。注意定义坐标系。

2) 高程数据修正：在该shp文件右键打开属性表，双击高程(Elevation)字段排序发现其中有部分等高线数据错误，利用编辑器删除错误选项。操作流程如下：点击编辑器/开始编辑/选中高程(Elevation)字段结果为0的等高线右键/删除所选项。

3) 生成TIN文件：3D Analyst工具/数据管理/TIN/创建TIN(注意此处增加坐标系)。

4) 生成DEM文件：3D Analyst工具/转换/由TIN转出/TIN转栅格。

5) 生成坡度：3D Analyst工具/栅格表面/坡度。

6) 坡度分类：3D Analyst工具/栅格重分类/重分类/分类，将坡度分类成为“≤25˚”和“>25˚”。

7) 计算两种坡度面积：打开坡度分级结果属性表/添加字段“area”/设置为浮点型/在“area”字段右键选择“字段计算器”/双击“count”字段后输入*30*30/得出两种坡度面积。

3.2.2. 计算项目区可剥离表土量

1) 将生成坡度后的文件格栅转面：转换工具/由栅格转出/栅格转面。

2) 计算项目区可剥离表土量：添加“永久占地地类”和“临时占地地类”图层，将两者分别与转换为面的坡度图层分别相交，随即得出永久占地和临时占地范围内不同坡度不同地类的图层，分别在各图层的属性表中增加字段“面积”，计算几何面积，方法同上。

3) 将得到的最终图层属性表导出。利用Excel数据透视功能统计出项目区可剥离表土面积，结合现场调查表土厚度情况，计算出项目区可剥离表土量(见附表2)，绘制项目区可剥离表土分布图(见图1)。注意：“gridcode”为坡度分类后的值，1代表坡度 ≤ 25˚，2代表坡度 > 25˚。

3.3. 项目区表土剥离与利用平衡分析

根据《水土保持工程设计规范》(GB 51018-2014)要求，项目区临时占地范围内林地和草地回覆表土20 cm厚，耕地回覆表土50 cm厚；永久占地范围内绿化区域回覆表土50 cm厚，项目区共计需表土量89.12万m³(见附表3)。本项目实际可剥离表土量为87.05万m³，无法满足要求。项目区周边合理运距范围内无其他生产建设项目存在多余表土，只能通过调整回覆厚度，达到表土平衡。本项目将交通道路区复耕回覆表土厚度从50 cm调整至40 cm，剥离的表土量即可满足项目回覆需求(见附表4)，达到挖填平衡。

注：1为表土分布图，2为可剥离表土分布图。

Figure 1.Map of topsoil distribution in localized areas of the project and map of stripping topsoil distribution

图1.项目局部区域表土分布图及可剥离表土分布图

4. 结语

1) 由于地形条件限制，某大型灌区生产建设项目可剥离表土范围被限制在25˚以下。本文利用ArcGIS的数据管理、数据编辑与处理、坐标系转换、矢量数据空间分析、DEM分析及栅格数据处理与分析功能，分析得出该项目扰动区域内所含表土量、可剥离表土量及分布范围。在此基础上对项目表土剥离与利用进行平衡分析，做好各防治分区之间的表土相互调运设计。

2) 目前尚未将ArcGIS应用于生产建设项目的表土平衡分析中，传统的设计思路为选取具有代表性的区域现场调查表土厚度，再根据调查结果与工程扰动面积相结合得出项目区所含表土量，进而指导后续表土平衡分析工作，未考虑或无法引入坡度因素。本方法可推广应用于所有生产建设项目表土平衡分析中，较传统的实地调查方法，可引入坡度因素，同时极大地提高了工作效率和结果的准确性。

3) 本方法的应用将便于设计人员根据现场情况指导施工，使得设计成果更具实用性；实现设计成果与地图动态互联。

4) 本方法对项目区地形图完整度要求较高，若扰动范围周边无较完整的等高线，则会影响到坡度的计算，从而影响结果的准确性。

附 录

TableS1.Table of topsoil distribution in the project area

附表1.项目区表土分布情况表

TableS2.Table of stripping topsoil in the project area

附表2.项目区可剥离表土情况表

TableS3.Table of topsoil requirements for the project

附表3.项目表土需求情况表

TableS4.Table of topsoil replenishments for the project

附表4.项目表土回覆情况表

NOTES

^*第一作者。

参考文献