DEM模型法與傳統(tǒng)等高線積分法計算庫容的比較

吳凱 何先寧 李勝才

摘? 要：水庫庫容是認(rèn)識水庫的一項十分重要的數(shù)據(jù)，水下的地形復(fù)雜多變，在野外采點計算的時候難以保證計算精準(zhǔn)度。該文探討了四川省某水電站水庫庫容的計算方法利用DEM模型法實現(xiàn)對庫容的精準(zhǔn)度計算，并與傳統(tǒng)等高線法計算庫容的方法進(jìn)行比較，以供同行參考。

關(guān)鍵詞：庫容;等高線法;數(shù)字高程模型

中圖分類號：P271? ? ? ? 文獻(xiàn)標(biāo)志碼：A

水庫庫容是認(rèn)識水庫的重要參數(shù)，對水庫的精準(zhǔn)度計算直接關(guān)系到整個水庫的防洪安全和發(fā)電灌溉效益。在計算水庫庫容的時候要對各項信息進(jìn)行綜合分析，既要保證數(shù)據(jù)信息的來源真實可靠，又要保證數(shù)據(jù)信息計算的科學(xué)準(zhǔn)確。該文對四川省某水電站水庫庫容進(jìn)行了探究，比較了在庫容計算中DEM模型法與傳統(tǒng)等高線法的區(qū)別。

1 工程項目概況

該項目名稱為四川省某水電站水庫庫容曲線復(fù)測。某水電站主要任務(wù)為發(fā)電，無供水、灌溉、防洪等，電站總裝機(jī)容量240 MW。水庫正常蓄水位2 930.00 m，設(shè)計洪水位2 930.09 m，校核洪水位2 930.93 m，死水位2 886 m。工程屬二等大（2）型規(guī)模，首部樞紐主要建筑物為2級，輸水工程為3級。

2 數(shù)字高程模型（DEM）法計算水庫庫容

2.1 數(shù)字高程模型（DEM）法介紹

數(shù)字高程模型（DEM）是用來描述地表起伏形態(tài)的空間數(shù)據(jù)模型，由地面規(guī)則格網(wǎng)點的高程值構(gòu)成矩陣，形成柵格結(jié)構(gòu)數(shù)據(jù)集，用來描述水庫底部信息空間分布的多種有序數(shù)值序列。高程數(shù)學(xué)模型主要有基于等高線、基于三角網(wǎng)和規(guī)則格網(wǎng)3種表現(xiàn)形式，特別是基于規(guī)則格網(wǎng)的DEM法，在實際測量工作中最為常用，因為其結(jié)構(gòu)簡單，使用起來十分方便。

2.2 DEM生產(chǎn)

計算水庫庫容是掌握水庫容量，計算水庫水量平衡的必要手段，計算結(jié)果是優(yōu)化水庫調(diào)度的數(shù)據(jù)基礎(chǔ)，庫容數(shù)據(jù)的準(zhǔn)確性直接影響到水庫的發(fā)電能力、防洪能力和灌溉水平等多種經(jīng)濟(jì)效益。過去，只能測量少量的斷面來進(jìn)行庫容計算，現(xiàn)在儀器設(shè)備、計算軟件的改進(jìn)則可較為精確的計算出庫容。

該項目水下測量采用GNSS-RTK配合多波束測深儀掃描的方法進(jìn)行。通過將所采集的多波束數(shù)據(jù)文件導(dǎo)入CARIS HIPS and SIPS軟件進(jìn)行內(nèi)業(yè)數(shù)據(jù)處理，輸出點云數(shù)據(jù)文件。對激光點云進(jìn)行去噪點、裁切，構(gòu)建地面三角網(wǎng)TIN模型，然后由地面三角網(wǎng)TIN模型融合計算出1 m格網(wǎng)間距的數(shù)字高程模型DEM。

2.3 采用DEM模型計算庫容

應(yīng)用DEM計算庫容時，其基本思路就是通過不規(guī)則的三角形面（TIN）盡可能地與水庫底下的地形曲面進(jìn)行精準(zhǔn)的結(jié)合，建立一套水下地形曲面模型，使建設(shè)的模型與被測水庫中的某一個高程平面進(jìn)行交互，進(jìn)而構(gòu)建成一套專屬于水庫庫容計算的實體模型。計算時由三角形的各個頂點向高程面引垂線，將實體模型分成好幾個三棱錐和三棱柱，之后再對分割而成的椎體進(jìn)行體積求和，所得到的體積便就是水庫庫容量。

將計算機(jī)構(gòu)建的DEM數(shù)據(jù)導(dǎo)入到ARCGIS軟件中，按照給定的要求計算出每個蓄水高程對應(yīng)的水庫面積及對應(yīng)的庫容量。該項目中計算蓄水高程范圍為2 859 m～2 930 m，按每1 m高程的步長計算庫容量。庫容計算高程起始值由測深儀所測水庫最低高程點高程值而定，截止高程值由甲方提供。

3 采用等高線積分法計算庫容

利用等高線積分法來計算庫容，需要以現(xiàn)有水下地形圖為依據(jù)，計算時主要是通過水庫等高線所圍成的面積和不同等高線之間的高度差異進(jìn)行，所采取的方法是通過數(shù)學(xué)微積分中的旋轉(zhuǎn)體積計算原理和計算方法來實現(xiàn)的。使用地形圖等高線法計算水庫庫容時，需要準(zhǔn)備的資料數(shù)據(jù)包括壩址、等高線、等高距和水位。利用等高線計算是一種測試方法簡單，經(jīng)濟(jì)投入較低的方法，計算精度比斷面法要高。

4 2種方法的比較

DEM模型計算法原理淺顯易懂，最大的優(yōu)點在于能充分利用全部測繪數(shù)據(jù)，精度較常規(guī)的等高線積分法和斷面法有很大提升，并且能夠計算出非整數(shù)高程的庫容量，計算精準(zhǔn)度也極高，而且工作效率也明顯提高。只需將整理好的點云數(shù)據(jù)輸入軟件，就可以自動進(jìn)行計算，這種計算方法適應(yīng)于所有庫容計算場景。另外庫容計算的精準(zhǔn)度和結(jié)果可靠性與DEM的精度有直接聯(lián)系。DEM的精度由采集原始數(shù)據(jù)時的誤差和高程內(nèi)插誤差決定，原始數(shù)據(jù)采集時會因為設(shè)備精度、人為操作和轉(zhuǎn)換坐標(biāo)等導(dǎo)致誤差產(chǎn)生，這種誤差可以采取一些方法減小，將誤差控制在適度范圍內(nèi);高程內(nèi)插誤差則是與采點方式和數(shù)學(xué)方法有關(guān)?；贒EM的庫容計算可大幅度降低工作人員的計算負(fù)擔(dān)，讓繁重的工作變簡單，提高了計算精準(zhǔn)度，計算工作自動化則減小了出錯的概率。

地形圖上的等高線主要是連接地圖等高程的點而形成的，2條等高線所圍成的體積根據(jù)已知的高差進(jìn)行求解，計算時要求等高線呈現(xiàn)閉合狀態(tài)。等高線積分法雖然使用較為簡單，但使用時需假設(shè)2點之間的坡度是均勻的，所以明顯是為了計算而假設(shè)的不存在的地形條件，因此只能用于地面起伏均勻的地形。以下為2種方法的計算成果比較。

從表1可知，采用DEM模型解析計算和等高線積分法計算2種結(jié)果的差值很小，其中死水位2 886 m高程相差比例為0.109%，正常蓄水位2 930 m高程相差比例為0.040%，說明2套數(shù)據(jù)可靠，考慮到等高線積分法計算時默認(rèn)為相鄰等高線間幾何形狀為線性，適合水下地形較為規(guī)整的水庫庫容計算，而DEM模型格網(wǎng)間距很小，數(shù)據(jù)量大，模型精度較高，適用于各種地形的庫容計算，所以本次計算成果采用DEM模型解析計算結(jié)果。

5 結(jié)語

水庫庫容事關(guān)附近百姓生計和水庫經(jīng)濟(jì)效益與社會效益，所以在計算過程中提升精確度能夠保證庫容結(jié)果的可靠性。該文介紹了在計算水庫庫容中的DEM模型法和等高線法，并對這2種方法進(jìn)行了比較。采用DEM模型法比傳統(tǒng)等高線法計算更具有優(yōu)勢，能更為準(zhǔn)確可靠地獲得庫容計算結(jié)果，為水庫運行提供真實有效的數(shù)據(jù)參考。

參考文獻(xiàn)

[1]楊玉棟.大型水利工程DEM模型構(gòu)建的方法[J].河南水利與南水北調(diào)，2018，47（8）：86-88.

[2]王耀革，朱長青，王志偉.數(shù)字高程模型（DEM）的整體誤差分析[J].武漢大學(xué)學(xué)報（信息科學(xué)版），2009，34（12）：1467-1470.