基于三角網(wǎng)數(shù)據(jù)文件的水庫(kù)庫(kù)容計(jì)算

楊勝保，何嬋軍，周 盼

（1.湖北省水利水電規(guī)劃勘測(cè)設(shè)計(jì)院，湖北 武漢430064）

基于三角網(wǎng)數(shù)據(jù)文件的水庫(kù)庫(kù)容計(jì)算

楊勝保1，何嬋軍1，周 盼1

（1.湖北省水利水電規(guī)劃勘測(cè)設(shè)計(jì)院，湖北 武漢430064）

提出基于三角網(wǎng)數(shù)據(jù)文件的計(jì)算庫(kù)容方法，推導(dǎo)三角網(wǎng)庫(kù)容計(jì)算嚴(yán)密公式。通過(guò)與南方CASS軟件計(jì)算結(jié)果對(duì)比表明，不僅計(jì)算結(jié)果精度高、速度快，還能一次計(jì)算多個(gè)水位的庫(kù)容量等。

庫(kù)容；三角網(wǎng)；南方CASS軟件；精度

大型水庫(kù)庫(kù)容曲線是水庫(kù)運(yùn)行調(diào)度與管理的基本依據(jù)，關(guān)系到水庫(kù)長(zhǎng)期效益的發(fā)揮。水庫(kù)庫(kù)容是指水庫(kù)蓄水后的存水水體，是水庫(kù)設(shè)計(jì)中的一項(xiàng)重要技術(shù)指標(biāo)。比較常用的庫(kù)容計(jì)算方法有：等高線法、斷面法、規(guī)則格網(wǎng)DEM法及TIN三角網(wǎng)法等[1-7]。這幾種方法均存在一定局限性。本文從三角網(wǎng)庫(kù)容計(jì)算模型出發(fā)，推導(dǎo)出嚴(yán)密的庫(kù)容計(jì)算公式，最大程度地消除模型誤差；基于三角網(wǎng)數(shù)據(jù)文件進(jìn)行庫(kù)容計(jì)算，使大區(qū)域劃分為若干個(gè)小區(qū)域變得簡(jiǎn)單快捷，避免劃分誤差，解決了硬件設(shè)施跟不上計(jì)算需求的問(wèn)題。

1 庫(kù)容模型公式

本文基于三角網(wǎng)進(jìn)行庫(kù)容計(jì)算，即將整個(gè)庫(kù)容分解為若干個(gè)拓?fù)潢P(guān)系相連的頂面為空中三角形的五面體，當(dāng)水面高程為H時(shí)，假定每個(gè)空中三角形貢獻(xiàn)的水體積為Vi、水面積為Si，則整個(gè)庫(kù)容及水面積分別為：

假定空中三角形△ABC為三角網(wǎng)中任意三角形，首先根據(jù)頂點(diǎn)三維坐標(biāo)易求得△ABC三邊長(zhǎng)a，b，c及水平投影三角形△A1B1C1三邊長(zhǎng)a1，b1，c1；再根據(jù)余弦定理可求得三頂角∠A、∠B、∠C及∠A1、∠B1、∠C1；根據(jù)三角形面積公式進(jìn)一步求得三角形面積SABC

及SA1B1C1

假定△ABC三頂點(diǎn)高程HA≥HC≥HB，當(dāng)水位為H時(shí)，△ABC所貢獻(xiàn)的水面積Si及水體積Vi有以下4種情況。

1）當(dāng)H≤HB時(shí)，如圖1所示。

圖1 無(wú)頂點(diǎn)低于水面

圖2 一個(gè)頂點(diǎn)低于水面

根據(jù)三角形面積公式：

根據(jù)三角形邊長(zhǎng)等比定理可求得：

將式（6）、（7）代入式（5）可得：

根據(jù)三棱錐體積公式：

3）當(dāng)HC≤H≤HA時(shí)，如圖3所示。

圖3 兩個(gè)頂點(diǎn)低于水面

根據(jù)三角形邊長(zhǎng)等比定理可求得：

將式（11）、（12）代入式（10）的三角形面積公式中可得：

利用立體幾何割補(bǔ)法可得：

根據(jù)三棱錐體積公式：

ha1為三角形△A1B1C1邊長(zhǎng)a1的垂高，將式（16）、（17）代入式（15）得：

再根據(jù)立體幾何割補(bǔ)法可得：

根據(jù)三棱柱體積公式：

根據(jù)三棱錐體積公式：

將式（20）、（21）代入式（19）得：

將式（18）、（22）代入式（14）得：

（4）當(dāng)H≥HA時(shí)，示意圖如圖4所示。

圖4 3個(gè)頂點(diǎn)低于水面

根據(jù)立體幾何割補(bǔ)法可得：

類比式（18）可得：

根據(jù)三棱柱體積公式可得：

將式（26）、（27）代入式（25）中得：

2 三角網(wǎng)數(shù)據(jù)文件

南方CASS軟件可根據(jù)DTM法、斷面法、方格網(wǎng)法、等高線法計(jì)算庫(kù)容，唯有三角網(wǎng)數(shù)據(jù)可以自由讀出與寫(xiě)入，且數(shù)據(jù)格式簡(jiǎn)單。三角網(wǎng)構(gòu)網(wǎng)是否合理直接影響到庫(kù)容計(jì)算的精度，一般采用軟件自動(dòng)生成配合人工干預(yù)的方式生成三角網(wǎng)。編輯處理好三角網(wǎng)后，在“等高線”命令欄選擇-“三角網(wǎng)存取”-“寫(xiě)入文件”-輸入文件名即可獲得三角網(wǎng)數(shù)據(jù)文件，三角網(wǎng)示意圖如圖5所示。

圖5 三角網(wǎng)

圖6 三角網(wǎng)數(shù)據(jù)文件

三角網(wǎng)數(shù)據(jù)文件格式如圖6所示，每一個(gè)坐標(biāo)數(shù)據(jù)占用一行、每三行數(shù)據(jù)為一個(gè)點(diǎn)的坐標(biāo)、每九行數(shù)據(jù)為一個(gè)三角形的3個(gè)頂點(diǎn)數(shù)據(jù)，以“nil”標(biāo)志文件結(jié)束。圖6中紅色方框代表一個(gè)三角形單元的坐標(biāo)數(shù)據(jù)，如三角網(wǎng)數(shù)據(jù)量過(guò)大，可從第一行數(shù)據(jù)開(kāi)始，以9的整數(shù)行將大數(shù)據(jù)文件劃分為若干個(gè)小數(shù)據(jù)文件，既避免了在三角網(wǎng)圖形上劃分子區(qū)域時(shí)易重復(fù)或遺漏三角形的問(wèn)題，也解決了內(nèi)存占用大、計(jì)算時(shí)導(dǎo)致死機(jī)的問(wèn)題。

3 庫(kù)容計(jì)算程序設(shè)計(jì)

通過(guò)分析推算的計(jì)算模型及三角網(wǎng)數(shù)據(jù)文件，可通過(guò)編寫(xiě)程序計(jì)算庫(kù)容?？紤]到繪制庫(kù)容曲線需多個(gè)水位的庫(kù)容量，可設(shè)置起止水面高程值、步長(zhǎng)值等輸入框。

4 算例分析

某庫(kù)區(qū)采用“GPS RTK+測(cè)深儀”方式測(cè)量了部分區(qū)域水下地形，通過(guò)CASS軟件自動(dòng)構(gòu)網(wǎng)及人工編輯生成一個(gè)由175 004個(gè)三角形構(gòu)成的三角網(wǎng)，寫(xiě)入三角網(wǎng)數(shù)據(jù)文件后，采用自行設(shè)計(jì)程序ResCap計(jì)算不同水位的庫(kù)容值，與CASS軟件的計(jì)算結(jié)果對(duì)比如表1所示。

從表1中可以看到：ResCap與CASS兩者計(jì)算的庫(kù)容結(jié)果幾乎相同，存在的微小差值可能由于進(jìn)位誤差所致。此外，ResCap還能根據(jù)設(shè)置的水位步長(zhǎng)計(jì)算多個(gè)庫(kù)容量和水面積，進(jìn)一步擬合可得到各種關(guān)系曲線，如圖7、8、9所示。

表1 庫(kù)容計(jì)算結(jié)果

圖7 庫(kù)容量-水位關(guān)系曲線

圖8 水面積-水位關(guān)系曲線

圖9 庫(kù)容量-水面積關(guān)系曲線

5 結(jié) 語(yǔ)

提出一種新的水庫(kù)庫(kù)容計(jì)算方法，推導(dǎo)了庫(kù)容模型及計(jì)算公式，基于三角網(wǎng)數(shù)據(jù)文件給出開(kāi)發(fā)庫(kù)容計(jì)算程序的解決辦法，解決了計(jì)算機(jī)硬件設(shè)施跟不上計(jì)算需求的問(wèn)題，提高了計(jì)算速率。

[1] 杜玉柱．地形法計(jì)算庫(kù)容的公式分析[J]．水文,2008,28(4):54-56

[2] 譚德寶,申邵洪．基于規(guī)則格網(wǎng)DEM的庫(kù)容計(jì)算與精度分析[J]．長(zhǎng)江科學(xué)院學(xué)報(bào),2009,26(3):49-52

[3] 王媛媛,段建剛．基于規(guī)則格網(wǎng)DEM計(jì)算水庫(kù)庫(kù)容方法研究[J]．吉林水利,2013(11):9-10

[4] 劉迎新．淺談利用南方cass計(jì)算三門峽水庫(kù)庫(kù)容的方法[J]．陜西水利,2012(6):133-134

[5] 楊廷利．基于ArcGIS計(jì)算水庫(kù)庫(kù)容的方法研究[J]．現(xiàn)代測(cè)繪,2013,36(4):51-52

[6] 武長(zhǎng)松,常軍．三維地面模型在大藤峽靜態(tài)庫(kù)容計(jì)算中的應(yīng)用[J]．城市勘測(cè),2011,(5):110-112

[7] 程劍剛． AutoCAD Civil 3D在湖泊庫(kù)容計(jì)算中的應(yīng)用[J]．測(cè)繪與空間地理信息,2014,37(3):137-118

P208

B

1672-4623（2017）10-0082-03

10.3969/j.issn.1672-4623.2017.10.027

2017-08-18。

楊勝保，高級(jí)工程師，主要從事水利工程測(cè)量與行政管理工作。