DOI：10.19981/j.CN23-1581/G3.2023.32.014

摘" 要：在水利工程前期規(guī)劃與設(shè)計(jì)中，土方量的計(jì)算方法不僅影響到方案布置，而且其計(jì)算結(jié)果還會(huì)直接影響到工程成本預(yù)算與項(xiàng)目組織問(wèn)題?；贕IS技術(shù)計(jì)算工程土方量的方法是利用DEM提供的高程信息，結(jié)合ArcGIS平臺(tái)的分析工具實(shí)現(xiàn)土方量快速、高效的計(jì)算，尤其是對(duì)于數(shù)據(jù)量大的計(jì)算尤為明顯。針對(duì)傳統(tǒng)DEM法計(jì)算工程土方量過(guò)程中存在重復(fù)作業(yè)、自動(dòng)化程度低等問(wèn)題，該文在DEM法基礎(chǔ)上，在ArcGIS平臺(tái)中基于Python語(yǔ)言結(jié)合特定的ArcPy函數(shù)編寫一種工作流模式的應(yīng)用工具，實(shí)現(xiàn)工程土方量計(jì)算的自動(dòng)化處理。通過(guò)優(yōu)化后的計(jì)算方法對(duì)洪澤湖（泗洪縣）退圩環(huán)湖工程中清退圩埂的土方量進(jìn)行計(jì)算和分析，并得到可靠的計(jì)算結(jié)果，為水利工程提供一種全新的土方算量解決方案。

關(guān)鍵詞：數(shù)字高程模型；土方量計(jì)算；ArcPy；自動(dòng)化；DEM法

中圖分類號(hào)：P208" " " " 文獻(xiàn)標(biāo)志碼：A" " " " " 文章編號(hào)：2095-2945（2023）32-0054-04

Abstract： In the early planning and design of water conservancy project， the calculation method of earthwork not only affects the scheme layout， but also directly affects the project cost budget and project organization. The method of calculating engineering earthwork based on GIS technology is to use the elevation information provided by DEM， combined with the analysis tools of ArcGIS platform to achieve fast and efficient earthwork calculation， especially for the calculation of large amount of data. In order to solve the problems of repeated work and low automation in the process of calculating engineering earthwork by traditional DEM method， on the basis of DEM method， this paper compiles an application tool of workflow pattern based on Python language and specific ArcPy function in ArcGIS platform to realize the automatic processing of engineering earthwork calculation. Through the optimized calculation method， the earthwork volume of the retreating polder ridge in the project around Hongze Lake （Sihong County） is calculated and analyzed， and reliable calculation results are obtained， which provides a new earthwork calculation solution for water conservancy projects.

Keywords： digital elevation model; earthwork calculation; ArcPy; automation; DEM method

數(shù)字高程模型（Digital Elevation Model，簡(jiǎn)稱DEM）是對(duì)地表形態(tài)的數(shù)字化表達(dá)，蘊(yùn)含了豐富的地學(xué)應(yīng)用分析所必需的地形地貌信息[1]。隨著GIS逐漸成為地理學(xué)研究的第三代語(yǔ)言，DEM也逐步替代等高線成為地形描述與分析的主要信息載體[2]。基于DEM數(shù)據(jù)與GIS技術(shù)相結(jié)合的理論、方法與應(yīng)用，是地理學(xué)、土木工程、地質(zhì)災(zāi)害分析等領(lǐng)域研究的熱點(diǎn)。土方量計(jì)算是工程勘測(cè)設(shè)計(jì)的重要內(nèi)容之一，準(zhǔn)確快速地計(jì)算土方量對(duì)詳細(xì)設(shè)計(jì)、工程量估算、資金分配等具有重要意義。近年來(lái)，DEM被越來(lái)越多地應(yīng)用在土方量計(jì)算中[3-4]，羅德仁等[5]、徐志敏等[6]的研究表明DEM法計(jì)算土方量精度較高、適用場(chǎng)景多，具有較高可行性。

洪澤湖是淮河流域最大的淺水型湖泊，地處蘇北平原中部偏西，位于淮河中下游結(jié)合部。近年來(lái)，大面積的圍墾使得湖泊的調(diào)蓄、行洪、生態(tài)和供水等功能受到了影響，制約了區(qū)域經(jīng)濟(jì)社會(huì)的可持續(xù)發(fā)展。為落實(shí)《江蘇省洪澤湖保護(hù)規(guī)劃》要求，提高洪澤湖防洪排澇能力，改善湖泊的水環(huán)境，解決歷史遺留問(wèn)題，推動(dòng)區(qū)域經(jīng)濟(jì)社會(huì)和諧發(fā)展，開展洪澤湖退圩還湖工程。本文擬綜合運(yùn)用DEM數(shù)據(jù)和GIS技術(shù)，及工程規(guī)劃和設(shè)計(jì)的理論，在ArcGIS軟件平臺(tái)上結(jié)合特定的ArcPy函數(shù)編寫的土方計(jì)算工具，最后以洪澤湖（泗洪縣）退圩還湖工程為例，進(jìn)行自動(dòng)化測(cè)算該工程中蓄水范圍線內(nèi)圩埂的土方量，為該工程規(guī)劃設(shè)計(jì)方案選優(yōu)、工程費(fèi)用概算與施工組織提供數(shù)據(jù)基礎(chǔ)。

1" DEM法

1.1" DEM法的原理

除傳統(tǒng)的方格網(wǎng)法、斷面法外，DEM法是當(dāng)前最常見(jiàn)的土方量計(jì)算方法。DEM的建立主要基于不規(guī)則分布采樣點(diǎn)、規(guī)則格網(wǎng)分布采樣點(diǎn)和等高線分布采樣點(diǎn)。本文DEM的建立是嘗試在不規(guī)則分布采樣點(diǎn)基礎(chǔ)上加入具有高程屬性的線狀要素作為不規(guī)則三角網(wǎng)（TIN）模型的特征約束條件，然后在TIN上進(jìn)行格網(wǎng)點(diǎn)的高程內(nèi)插形成柵格數(shù)據(jù)。基于DEM土方量計(jì)算的原理實(shí)質(zhì)上是計(jì)算現(xiàn)狀地表面與規(guī)劃設(shè)計(jì)地表面體積之差，負(fù)值在工程中為挖方量，反之正值稱為填方量。

1.2" DEM法的具體步驟

DEM法常利用研究區(qū)內(nèi)現(xiàn)有數(shù)字線畫圖（DLG）成果提取高程要素，包括高程點(diǎn)、等高線和具有高程屬性的其他線狀要素。在ArcGIS軟件平臺(tái)中對(duì)高程要素進(jìn)行矢量化處理，為確保數(shù)據(jù)質(zhì)量，套合DOM及外業(yè)實(shí)測(cè)數(shù)據(jù)等成果再對(duì)矢量數(shù)據(jù)進(jìn)行高程異常值的檢查。使用ArcGIS的3D Analyst對(duì)研究區(qū)矢量數(shù)據(jù)進(jìn)行擬合生成TIN模型，在此基礎(chǔ)上，通過(guò)線性內(nèi)插獲取DEM。對(duì)需要進(jìn)行范圍約束的柵格數(shù)據(jù)需要進(jìn)行裁剪。通過(guò)Spatial-Analyst模塊的疊加分析和統(tǒng)計(jì)分析等功能對(duì)現(xiàn)狀DEM和設(shè)計(jì)DEM進(jìn)行分析和統(tǒng)計(jì)得到該區(qū)域的填挖方成果，本文總結(jié)的DEM法操作流程如圖1所示。

2" 自動(dòng)化土方量計(jì)算的實(shí)現(xiàn)

2.1" 研究思路和研究方法

常規(guī)的基于DEM模型在ArcGIS平臺(tái)上進(jìn)行土方量計(jì)算有創(chuàng)建模型簡(jiǎn)單、可視化程度高等優(yōu)勢(shì)。但水利工程開展規(guī)劃和設(shè)計(jì)過(guò)程中有著方案反復(fù)調(diào)整的特點(diǎn)，其中作業(yè)范圍、設(shè)計(jì)高程的調(diào)整更是需要重新計(jì)算相應(yīng)的土方量，生產(chǎn)實(shí)際過(guò)程中存在重復(fù)作業(yè)的問(wèn)題。此外，當(dāng)作業(yè)范圍內(nèi)分區(qū)較多時(shí)，人工方式作業(yè)效率較低、自動(dòng)化程度不高。本文針對(duì)常規(guī)DEM法計(jì)算的這些不足，利用ArcGIS中的Arcpy站點(diǎn)包進(jìn)行土方量自動(dòng)計(jì)算工具編寫，較大程度上避免了重復(fù)作業(yè)，減輕工作強(qiáng)度的同時(shí)也提高了作業(yè)效率和成果質(zhì)量。

站點(diǎn)包ArcPy是ArcGIS10.0之后版本推出的一個(gè)針對(duì)所有許可均可使用的python腳本模塊，這個(gè)站點(diǎn)包把ArcGIS中所有功能都封裝成一個(gè)個(gè)類庫(kù)，只要在Python的編輯器中導(dǎo)入Arcpy站點(diǎn)包，使用者就可以調(diào)用里面的類庫(kù)從而實(shí)現(xiàn)其需要的功能，并根據(jù)自己需要開發(fā)相應(yīng)的工具[7]。文中所涉及工具大部分以ArcGIS中自帶分析工具為主，容易實(shí)現(xiàn)，可以為其他單位提供參考。

2.2" 土方量自動(dòng)計(jì)算的實(shí)現(xiàn)

本研究以ArcGIS為平臺(tái)，首先將計(jì)算區(qū)范圍、現(xiàn)狀DEM數(shù)據(jù)、約束范圍數(shù)據(jù)（可選）預(yù)處理完畢，其中計(jì)算區(qū)范圍矢量數(shù)據(jù)中需要新建屬性字段，將工程方案提供的設(shè)計(jì)高程值根據(jù)分區(qū)名稱批量賦值到該字段中，系統(tǒng)將根據(jù)設(shè)計(jì)高程值自動(dòng)創(chuàng)建設(shè)計(jì)DEM。在土方量自動(dòng)計(jì)算工具中只需要輸入分區(qū)范圍矢量數(shù)據(jù)、現(xiàn)狀DEM數(shù)據(jù)和約束范圍數(shù)據(jù)（可選），最后設(shè)置好成果輸出路徑并輸出處理成果。最終成果為Excel格式的工程各分區(qū)土方量統(tǒng)計(jì)表，過(guò)程數(shù)據(jù)統(tǒng)一存放在設(shè)定好的臨時(shí)文件夾中。土方自動(dòng)計(jì)算工具技術(shù)路線圖如圖2所示。土方量自動(dòng)計(jì)算工具在洪澤湖（淮陰區(qū)）退圩還湖工程先導(dǎo)段中的應(yīng)用驗(yàn)證了方法的可靠性，有較好的實(shí)際應(yīng)用價(jià)值。

3" 實(shí)例分析

洪澤湖（泗洪縣）退圩還湖工程規(guī)劃主要對(duì)洪澤湖蓄水范圍內(nèi)的違章堤圩和埂圍的圩埂土方進(jìn)行清退，本次擬清退分區(qū)54個(gè)，主要涉及泗洪縣臨淮鎮(zhèn)、半城鎮(zhèn)、龍集鎮(zhèn)和界集鎮(zhèn)，面積約為125.5 km2。筆者參與了該工程土方量計(jì)算工作，考慮到開挖范圍內(nèi)高程不一，靠近岸邊高程相對(duì)較高，遠(yuǎn)離岸邊相對(duì)較低，清退區(qū)底高程變化一般在11.5～12.5 m。根據(jù)泗洪縣洪澤湖地形、防洪、供水及水生態(tài)恢復(fù)需要，圩埂清退標(biāo)準(zhǔn)原則上清退至現(xiàn)狀湖底高程。

清退范圍內(nèi)養(yǎng)殖魚塘密布，圩埂縱橫交錯(cuò)，根據(jù)給定的設(shè)計(jì)高程值，采用了基于DEM的土方量自動(dòng)計(jì)算工具進(jìn)行土方量計(jì)算。

3.1" 數(shù)據(jù)預(yù)處理

建立DEM的第一步是獲取地形數(shù)據(jù)，本文采用的DEM地形數(shù)據(jù)主要來(lái)源于2019年研究區(qū)1∶1 000比例尺地形圖。對(duì)地形圖進(jìn)行高程要素提取，獲取特征點(diǎn)線并進(jìn)行矢量化。利用AcrGIS軟件的 3DAnalyst模塊根據(jù)現(xiàn)狀圖中矢量數(shù)據(jù)生成研究區(qū)域現(xiàn)狀TIN，線性內(nèi)插獲取現(xiàn)狀DEM。由于該工程規(guī)劃未涉及到養(yǎng)殖魚塘等水體部分的清退，因此需要對(duì)計(jì)算范圍進(jìn)行去除水體部分的處理，即創(chuàng)建圩埂矢量數(shù)據(jù)作為約束范圍用于DEM的裁切。

3.2" 工程填挖方計(jì)算

在土方量自動(dòng)計(jì)算工具中，首先輸入研究區(qū)各分區(qū)范圍（含有各分區(qū)設(shè)計(jì)高程屬性字段），再輸入前期構(gòu)建好的現(xiàn)狀DEM，由于本次清退主體是圩埂，因此在圩埂范圍中輸入繪制好的圩埂矢量數(shù)據(jù)（可選項(xiàng)）用于裁剪，最后設(shè)置好過(guò)程數(shù)據(jù)存放路徑和統(tǒng)計(jì)數(shù)據(jù)表存放路徑點(diǎn)擊運(yùn)行完成，計(jì)算工具操作界面如圖3所示。

3.3" 數(shù)據(jù)成果與分析

土方量自動(dòng)計(jì)算工具在內(nèi)部對(duì)疊加分析得到的填挖方柵格數(shù)據(jù)面積屬性字段數(shù)值進(jìn)行統(tǒng)計(jì)。由于項(xiàng)目要求只需要計(jì)算挖方量，因此只需要選擇凈損失區(qū)域的土方量即可，工具對(duì)自動(dòng)選擇土方量數(shù)值大于0的求和，最終成果見(jiàn)表1，表中包含序號(hào)、分區(qū)名稱、設(shè)計(jì)標(biāo)高以及土方量數(shù)據(jù)成果。

由表1可知，項(xiàng)目區(qū)總的挖方量為2 530.2萬(wàn)m3，其中分區(qū)LH2、BC1、BC2、BC4、CW1、LJ2、LJ4、379BC3、25CT7挖方量較大，超過(guò)100萬(wàn)m3，清退重點(diǎn)主要集中在臨淮鎮(zhèn)和半城鎮(zhèn)范圍內(nèi)。

4" 結(jié)束語(yǔ)

基于GIS技術(shù)的工程土方計(jì)算方法是以DEM法為基礎(chǔ)，較之傳統(tǒng)的土方量計(jì)算法有更真實(shí)反映現(xiàn)場(chǎng)地形地貌、計(jì)算簡(jiǎn)單快速、數(shù)據(jù)成果可視化的優(yōu)勢(shì)?；贕IS技術(shù)的工程土方計(jì)算方法優(yōu)化是在工程規(guī)劃設(shè)計(jì)方案反復(fù)調(diào)整需要重復(fù)計(jì)算、人工作業(yè)效率不高的背景下，較之常規(guī)的DEM法有自動(dòng)化程度高、操作簡(jiǎn)單等優(yōu)勢(shì)。土方量自動(dòng)計(jì)算工具以現(xiàn)狀DEM數(shù)據(jù)為基礎(chǔ)，通過(guò)一系列處理，再對(duì)每一個(gè)區(qū)域土方量進(jìn)行統(tǒng)計(jì)，最后統(tǒng)計(jì)分析出整個(gè)工程的土方量。在應(yīng)對(duì)工程規(guī)劃范圍調(diào)整、設(shè)計(jì)高程調(diào)整時(shí)更加靈活、高效，土方量成果為退圩還湖工程制定合法、合規(guī)、合理的補(bǔ)償、補(bǔ)助標(biāo)準(zhǔn)和安置方式提供依據(jù)。

參考文獻(xiàn)：

[1] 湯國(guó)安，李發(fā)源，劉學(xué)軍.數(shù)字高程模型教程[M].3版.北京：科學(xué)出版社，2016.

[2] 湯國(guó)安.我國(guó)數(shù)字高程模型與數(shù)字地形分析研究進(jìn)展[J].地理學(xué)報(bào)，2014（9）：1305-1325.

[3] 劉錦霞，張超，邵林，等.Civil 3D在退圩還湖設(shè)計(jì)項(xiàng)目中的應(yīng)用[J].江蘇水利，2019（5）：12-15.

[4] 廖偉華，徐彬.基于空間內(nèi)插的場(chǎng)地平整土方工程量計(jì)算方法研究[J].黑龍江工程學(xué)院學(xué)報(bào)（自然科學(xué)版），2007（4）：14-17.

[5] 羅德仁，鄒自力，湯江龍.工程土方量計(jì)算比較分析[J].東華理工學(xué)院學(xué)報(bào)，2005（1）：59-64.

[6] 徐志敏，林志勇，李雯靜，等.基于填挖方分析的DEM精度評(píng)價(jià)模型[J].武漢大學(xué)學(xué)報(bào)：信息科學(xué)版，2017（8）：1167-1171.

[7] 余詠勝，彭艷麗，尹言軍，等.基于Arcpy的影像地圖自動(dòng)處理技術(shù)研究[J].測(cè)繪通報(bào)，2015（3）：82-85.

第一作者簡(jiǎn)介：鄭亞運(yùn)（1990-），男，碩士，工程師。研究方向?yàn)樗こ虦y(cè)量、水利信息化。