Application of Topsoil Stripping and Utilization Balance Analysis in Production and Construction Projects Based on ArcGIS
Topsoil stripping and utilization balance analysis is one of the important contents of soil and water conservation design in production and construction projects, and the accuracy requirements of its design results are constantly improving in the intra-industry. Taking the topsoil design of a large-scale irrigation area as an example, this paper applies the spatial data processing ability of ArcGIS to the project topsoil stripping and utilization design work under the specified slope requirements, which greatly improves the work efficiency and the accuracy of the results. Based on this, it can be applied to all production and construction projects.
Topsoil Stripping and Utilization
土壤作为陆地生物养分之源极其宝贵,但其形成极为缓慢。表土是指最表层的土壤,作为土壤中养分含量最高的部分,能有利于快速恢复地力和植物生长,珍贵程度更甚
本文基于ArcGIS对某大型灌区生产建设项目的表土剥离与利用进行平衡分析。由于该项目区地形高低起伏较大,坡度大于25˚的林地、草地、耕地表土层较薄且含有大量石块,无法使用机械进行剥离,采用人工剥离成本较高,因此该区域表土不进行剥离。利用ArcGIS的空间数据处理能力分析得出项目区表土量、可剥离表土量,并进行后期利用平衡分析。
利用ArcGIS的数据管理、数据编辑与处理、坐标系转换、矢量数据空间分析,可得出项目扰动范围内表土分布情况、分布数量。通过对项目区地形图分析得出DEM文件,经栅格数据处理与分析,将项目扰动区域生成以坡度25˚为界限的适量数据,在此基础上与三调图相交,即可得出项目区可剥离表土范围。将以上处理成果的属性表输出,统计项目区表土量、可剥离表土量和表土需求量,进行后续表土平衡分析。在ArcGIS中亦可制作输出表土分布范围图、可剥离表土分布图。
1) 将项目扰动范围划分为永久占地和临时占地两个图层
由于项目施工永久占地范围与临时占地范围有重合部分,为防止重复计算,需将两者进行擦除分析,得出“施工永久占地”和“擦除后的施工临时占地”两个图层。操作流程如下:将两者cad格式文件导入ArcGIS中,用分析工具/叠加分析/擦除功能得出“擦除后的施工临时占地”。注意赋予两图层与项目区三调图相同的坐标系,否则无法进行后续相交分析。具体方法为:数据管理工具/投影和变换/栅格/定义投影。
2) 计算项目扰动范围内永久占地和临时占地范围内所含表土量
在ArcGIS中导入项目区三调图,分别于“擦除后的施工临时占地”和“施工永久占地”进行相交得出项目扰动范围各地类面积。操作流程如下:分析工具/叠加分析/相交;计算各图斑面积后导出属性表;操作步骤如下:右键打开属性表/添加字段“area”(浮点型)/在“area”右键计算几何;随后导出属性表。同理得出永久占地地类表。
根据现场调查,表土主要分布在耕地、园地、林地、草地等占地类型区。表土厚度依据地类不同差别较大,耕地表土厚度为30 cm,园地表土厚度约35 cm,乔木林地表土厚度约25 cm,灌木林地、竹林地表土厚度约20 cm,其他林地表土厚度约15 cm,草地表土厚度约20 cm。结合现场调查结果可分析出工程区所含表土量(见
利用项目区地形图生成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/得出两种坡度面积。
1) 将生成坡度后的文件格栅转面:转换工具/由栅格转出/栅格转面。
2) 计算项目区可剥离表土量:添加“永久占地地类”和“临时占地地类”图层,将两者分别与转换为面的坡度图层分别相交,随即得出永久占地和临时占地范围内不同坡度不同地类的图层,分别在各图层的属性表中增加字段“面积”,计算几何面积,方法同上。
3) 将得到的最终图层属性表导出。利用Excel数据透视功能统计出项目区可剥离表土面积,结合现场调查表土厚度情况,计算出项目区可剥离表土量(见
根据《水土保持工程设计规范》(GB 51018-2014)要求,项目区临时占地范围内林地和草地回覆表土20 cm厚,耕地回覆表土50 cm厚;永久占地范围内绿化区域回覆表土50 cm厚,项目区共计需表土量89.12万m3(见
1) 由于地形条件限制,某大型灌区生产建设项目可剥离表土范围被限制在25˚以下。本文利用ArcGIS的数据管理、数据编辑与处理、坐标系转换、矢量数据空间分析、DEM分析及栅格数据处理与分析功能,分析得出该项目扰动区域内所含表土量、可剥离表土量及分布范围。在此基础上对项目表土剥离与利用进行平衡分析,做好各防治分区之间的表土相互调运设计。
2) 目前尚未将ArcGIS应用于生产建设项目的表土平衡分析中,传统的设计思路为选取具有代表性的区域现场调查表土厚度,再根据调查结果与工程扰动面积相结合得出项目区所含表土量,进而指导后续表土平衡分析工作,未考虑或无法引入坡度因素。本方法可推广应用于所有生产建设项目表土平衡分析中,较传统的实地调查方法,可引入坡度因素,同时极大地提高了工作效率和结果的准确性。
3) 本方法的应用将便于设计人员根据现场情况指导施工,使得设计成果更具实用性;实现设计成果与地图动态互联。
4) 本方法对项目区地形图完整度要求较高,若扰动范围周边无较完整的等高线,则会影响到坡度的计算,从而影响结果的准确性。
附表1. 项目区表土分布情况表
占地类型 |
防治分区 |
含表土区域(m2) |
含表土 量(m3) |
|||||||||
耕地 |
园地 |
林地 |
草地 |
|||||||||
水田 |
旱地 |
果园 |
其他园地 |
乔木林地 |
竹林地 |
灌木林地 |
其他林地 |
其他草地 |
||||
临时占地 |
主体工程区 |
管线工程区 |
17,014 |
73,526 |
187,587 |
136,427 |
1,824,522 |
11,300 |
318,872 |
8462 |
86,061 |
681,213 |
交通道路区 |
20,911 |
46,528 |
131,369 |
5203 |
45,986 |
12,187 |
14,042 |
1468 |
541 |
85,102 |
||
施工生产生活区 |
0 |
11,654 |
14,162 |
18,961 |
22,595 |
335 |
14,908 |
1865 |
1632 |
24,393 |
||
弃渣场区 |
1377 |
7822 |
18,675 |
12,997 |
264,982 |
55,817 |
9288 |
93,111 |
||||
小计 |
39,302 |
139,529 |
351,792 |
173,588 |
2,158,085 |
23,822 |
403,639 |
11,795 |
97,522 |
883,820 |
||
永久占地 |
主体工程区 |
泵站 |
3178 |
401 |
826 |
1378 |
||||||
管线工程 |
295 |
2210 |
3855 |
23,248 |
88,518 |
646 |
29,038 |
73 |
4483 |
39,211 |
||
小计 |
295 |
2210 |
7032 |
23,248 |
88,919 |
646 |
29,864 |
73 |
4483 |
40,589 |
||
合计 |
39,596 |
141,739 |
358,825 |
196,836 |
2,247,004 |
24,469 |
433,503 |
11,869 |
102,005 |
924,409 |
附表2. 项目区可剥离表土情况表
占地类型 |
防治分区 |
可剥离表土区域(m2) |
可剥离 表土量(m3) |
|||||||||
耕地 |
园地 |
林地 |
草地 |
|||||||||
水田 |
旱地 |
果园 |
其他园地 |
乔木林地 |
竹林地 |
灌木林地 |
其他林地 |
其他草地 |
||||
临时占地 |
主体工程区 |
管线工程 |
16,647 |
51,792 |
183,741 |
129,503 |
1,705,435 |
11,300 |
295,201 |
8462 |
86,061 |
83,373 |
交通道路区 |
20,911 |
46,001 |
128,731 |
5203 |
44,998 |
12,187 |
12,035 |
1468 |
541 |
23,061 |
||
施工生产生活区 |
11,146 |
14,162 |
16,195 |
22,071 |
335 |
14,508 |
1865 |
1632 |
93,111 |
|||
弃渣场区 |
1377 |
7822 |
18,675 |
12,997 |
264,982 |
55,817 |
9288 |
636,308 |
||||
小计 |
38,935 |
116,762 |
345,308 |
163,898 |
2,037,487 |
23,822 |
377,561 |
11,795 |
97,522 |
835,853 |
||
永久占地 |
主体工程区 |
泵站 |
0 |
0 |
3178 |
0 |
401 |
0 |
826 |
0 |
0 |
1378 |
管线工程 |
49 |
1121 |
3855 |
20,598 |
75,845 |
646 |
22,215 |
73 |
4118 |
33,277 |
||
小计 |
49 |
1121 |
7032 |
20,598 |
76,246 |
646 |
23,041 |
73 |
4118 |
34,655 |
||
合计 |
38,984 |
117,883 |
352,341 |
184,496 |
2,113,733 |
24,469 |
400,602 |
11,869 |
101,640 |
870,509 |
附表3. 项目表土需求情况表
占地类型 |
防治分区 |
表土需求数量 |
|||||||
复耕 |
植被恢复 |
小计(m3) |
|||||||
覆土面积(hm2) |
厚度(cm) |
数量(m3) |
覆土面积(hm2) |
厚度(cm) |
数量(m3) |
||||
临时占地 |
主体工程区 |
管线工程 |
41.46 |
50 |
207,277 |
224.92 |
20 |
449,843 |
657,120 |
交通道路区 |
20.40 |
50 |
102,006 |
7.42 |
20 |
14,845 |
116,850 |
||
施工生产生活区 |
4.48 |
50 |
22,389 |
4.13 |
20 |
8267 |
30,656 |
||
弃渣场区 |
4.09 |
50 |
20,436 |
33.01 |
20 |
66,017 |
86,453 |
||
小计 |
70.42 |
352,107 |
269.49 |
538,973 |
891,079 |
||||
永久占地 |
主体工程区 |
泵站 |
0.02 |
50 |
120 |
120 |
|||
合计 |
70.42 |
352,107 |
269.51 |
539,093 |
891,199 |
附表4. 项目表土回覆情况表
占地类型 |
防治分区 |
表土回覆数量 |
|||||||
复耕 |
植被恢复 |
小计(m3) |
|||||||
覆土面积(hm2) |
厚度(cm) |
数量(m3) |
覆土面积(hm2) |
厚度(cm) |
数量(m3) |
||||
临时占地 |
主体工程区 |
管线工程 |
41.46 |
50 |
207,277 |
224.92 |
20 |
449,843 |
657,120 |
交通道路区 |
20.40 |
40 |
81,316 |
7.42 |
20 |
14,845 |
96,161 |
||
施工生产生活区 |
4.48 |
50 |
22,389 |
4.13 |
20 |
8267 |
30,656 |
||
弃渣场区 |
4.09 |
50 |
20,436 |
33.01 |
20 |
66,017 |
86,453 |
||
小计 |
70.42 |
331,417 |
269.49 |
538,973 |
870,389 |
||||
永久占地 |
主体工程区 |
泵站 |
0.02 |
50 |
120 |
120 |
|||
合计 |
70.42 |
331,417 |
269.51 |
539,093 |
870,509 |
*第一作者。