關(guān)于水庫庫容和淤積測量技術(shù)研究與思考

摘要：作者應(yīng)用全球定位系統(tǒng)（GPS）和回聲測深技術(shù)，對水庫庫容和淤積測量進(jìn)行了研究，并提出了三角形構(gòu)網(wǎng)方法，經(jīng)實(shí)際運(yùn)用取得了滿意的效果。

關(guān)鍵詞：GPS；庫容；淤積；測量

水庫庫容和淤積量是水庫調(diào)度的重要參數(shù)，其精度直接影響到水庫的防洪安全與蓄水興利。隨著現(xiàn)代測控技術(shù)的迅速發(fā)展，將其應(yīng)用到水利電力領(lǐng)域是一項(xiàng)值得研究和推廣的課題。

1 常規(guī)庫容及淤積量的確定

常規(guī)的庫容計(jì)算方法多采用斷面法。其庫區(qū)容量的計(jì)算模型為

式中：Vi、Li為第i個斷面到第i+1個斷面間的庫容和距離；n為分段個數(shù)；Si、m、d、hi分別為第i個斷面的面積、測點(diǎn)個數(shù)、點(diǎn)間距和每個測點(diǎn)的深度測量值。采用斷面法雖然操作簡單，但受前提假設(shè)的制約，精度難以保證。淤積量是根據(jù)前后兩次的的庫容較差獲得，庫容不準(zhǔn)確，淤積量的計(jì)算精度就無從談起。

2 高精度水下地形測量技術(shù)[1，2]

2.1 水下地形測量 所謂水下地形測量，就是利用測量儀器來確定水底點(diǎn)的三維坐標(biāo)的過程。隨著GPS技術(shù)的迅速發(fā)展，水下地形測量方法取得了很大的進(jìn)展。目前，水下地形測量技術(shù)已定型于采用GPS獲取平面坐標(biāo)，測深儀獲取深度數(shù)據(jù)的基本模式。同時，為了獲得水下地物的海拔高程，以及消除潮汐、水位落差等諸因素的影響，進(jìn)行水位監(jiān)測也是一個重要環(huán)節(jié)。水下地形測量現(xiàn)狀示意圖如下。針對上述現(xiàn)有模式，文獻(xiàn)[1]對測量設(shè)備的選型、基準(zhǔn)點(diǎn)的布設(shè)、潮汐觀測的具體實(shí)施等諸多技術(shù)問題，進(jìn)行了深入探討。

2.2 GPS載波相位差分定位技術(shù)和回聲測深技術(shù) 隨著GPS技術(shù)的發(fā)展，GPS日益廣泛應(yīng)用于水利電力工程的各個方面。為了提高定位精度，一般均采用差分技術(shù)。在眾多的差分技術(shù)中，偽距差分和載波相位差分是最為常用的兩種測量模式，后者的定位精度較高（厘米級），通常用于高精度的測量工程和研究中。

載波相位差分測量的定位精度很大程度上依賴于整周模糊度能否在航精確確定。整周模糊度在航解算（OTF）是一種動態(tài)環(huán)境下的模糊度確定方法，它可省去在精密動態(tài)定位中的的靜態(tài)初始化過程。常規(guī)精密定位中復(fù)雜的整周跳變問題也因OTF的引入變得十分簡單。載波相位差分測量整周模糊度的確定模型為

式中：Xk=（dx dy dz x y z dn0 dn1…dnm）為狀態(tài)向量；Φk-1為狀態(tài)轉(zhuǎn)移矩陣；Hk，φ為載波相位的測量矩陣；Rk，φ為載波相位的方差陣；Qk為系數(shù)陣。

由式（3）計(jì)算得到整周模糊度N后，代入載波相位觀測方程，便可以獲得厘米級甚至毫米級的平面定位精度。

回聲測深儀是一種單波束測深設(shè)備，深度的測量是根據(jù)最小聲程決定。按照使用頻率個數(shù)的不同，又可分為單頻和雙頻。雙頻測深儀根據(jù)兩個頻率測量深度較差獲得淤積層厚度。

2.3 高精度庫容和淤積量測量方法 庫容和淤積量的精密測量采用現(xiàn)代水下地形測量方法，即利用GPS載波相位差分測量技術(shù)進(jìn)行平面定位，測深儀進(jìn)行深度測量，GPS和測深儀保證同步作業(yè)，獲取水底測點(diǎn)平面和深度信息的作業(yè)模式。

為了保證庫容和淤積量的計(jì)算精度，需要對庫區(qū)進(jìn)行測線設(shè)計(jì)，GPS和測深采樣也要按照水下地形測量規(guī)范等間隔或等時間采樣。設(shè)測量比例尺為1：Scale，測量船的平均速度為，則測線間距d和時間間隔Δt為

為了提高測量精度，在測線布設(shè)時，還應(yīng)該考慮水下地形的變化趨勢，若變化相對比較平坦，則測線間距可以適當(dāng)放寬，否則，需加密測線。這有利于使測點(diǎn)均勻分布于整個測區(qū)，同時在測區(qū)水下地形變化復(fù)雜的地區(qū)使測點(diǎn)深度或高程能更好地反映水下地形的真實(shí)面貌。

3 庫容和淤積量的計(jì)算方法

3.1 庫容計(jì)算方法 為了提高計(jì)算精度，充分利用水下地形測量數(shù)據(jù)，本文提出了一種三角柱計(jì)算庫容的方法。該法建立在實(shí)際測點(diǎn)的基礎(chǔ)上，根據(jù)圖2，相鄰三個測點(diǎn)可構(gòu)成的三角柱體積為

設(shè)n為整個庫區(qū)三角形個數(shù)，則整個庫區(qū)的庫容為

3.2 淤積量計(jì)算方法 當(dāng)庫底為基巖構(gòu)造時，采用雙頻（f1、f2）測深儀測深，淤積量的計(jì)算方法同庫容相似。設(shè)相鄰三個測點(diǎn)在淤積層表面利用f1測得深度分別為h1、h2和h3，利用f2，在對應(yīng)點(diǎn)基巖上測得的深度分別為h′1、h′2和h′3，若設(shè)淤積層表面面積為S2，基巖上的面積為S3，則淤積量為

式中淤積層上下面的面積S2和S3的計(jì)算方法同式（5）。

則庫區(qū)的淤積量為

對于上述情況，V′的計(jì)算還可采用模型

式中：Vf1、Vf2分別代表根據(jù)f1、f2測得的淤積表面、基巖表面上的深度計(jì)算得到的體積。

然而，對于淤積層下地質(zhì)是非基巖的情況，式（9）的庫區(qū)淤積量計(jì)算模型就不再適用，而需要根據(jù)建庫初期的原始床面（地形圖）計(jì)算空庫容，或者前期確定的庫容量，與根據(jù)本次利用f1頻率測量的水深（淤積層表面的水深）計(jì)算所得庫容Vf1較差得到實(shí)際的庫區(qū)淤積量。其計(jì)算模型為

V′=Vf1-V前期庫容（10）

4 問題討論

本文所提出的基于現(xiàn)代水下地形測量技術(shù)的水庫庫容和淤積量確定方法相對于傳統(tǒng)的斷面法具有許多優(yōu)點(diǎn)，諸如定位精度高、計(jì)算結(jié)果準(zhǔn)確、所得數(shù)據(jù)可用于水下地形圖的繪制及DTM的建立等。然而，相對于傳統(tǒng)庫容和淤積量的確定方法，由于采用了先進(jìn)的測控設(shè)備，無疑會增大測量和計(jì)算方法上的復(fù)雜度，但這些是可以通過計(jì)算機(jī)編程來自動化實(shí)現(xiàn)的。現(xiàn)將上述方法在實(shí)際數(shù)據(jù)處理中的幾個難點(diǎn)加以討論。

（1）對于比較大的庫區(qū)，如江河形成的自然庫區(qū)，數(shù)據(jù)量會隨水域面積的增加而急劇增大。在利用這些數(shù)據(jù)構(gòu)造庫區(qū)三角形時會因存儲量和搜索范圍過大，占用過多的計(jì)算機(jī)內(nèi)存，可能會導(dǎo)致計(jì)算速度過慢或者死機(jī)。為克服這一問題，在三角形構(gòu)造中可采用一種快速的三角形構(gòu)網(wǎng)方法，即局域搜索法。根據(jù)測區(qū)范圍和測點(diǎn)的數(shù)量，可事先對整個區(qū)域根據(jù)坐標(biāo)進(jìn)行劃分，然后在結(jié)合拓展三角形的范圍索引各個分割區(qū)，在小區(qū)域內(nèi)實(shí)現(xiàn)快速搜索。這樣可以大大的節(jié)約計(jì)算機(jī)內(nèi)存，提高三角形的構(gòu)網(wǎng)速度。

（2）通過水下地形測量可給出水面以下的深度，以及根據(jù)水面下的實(shí)測結(jié)果計(jì)算水底到水面高程變化的庫容曲線，而對于高于當(dāng)前水面的水位面庫容曲線無法進(jìn)行計(jì)算和繪制。為了得到一個全面反映庫區(qū)容量變化的庫容曲線，需要將庫區(qū)邊緣數(shù)字高程信息引入庫容計(jì)算中。庫區(qū)邊緣陸地的數(shù)字高程信息可通過兩種途徑獲得。一種是利用GPS載波相位差分技術(shù)進(jìn)行動態(tài)地形測量獲得；另一種方法是通過已有的地形圖或DTM獲得。若利用GPS載波相位差分測量技術(shù)獲得陸地?cái)?shù)字信息，則GPS天線相位中心的平面位置即為陸地測點(diǎn)的平面位置，相位中心的高程減去天線高便是陸地高程。

（3）在（2）中，已有地形圖與現(xiàn)有測量成果共同用于庫容曲線計(jì)算時，兩套資料的高程和坐標(biāo)基準(zhǔn)必須匹配。對于將水底點(diǎn)的深度轉(zhuǎn)換成高程問題，傳統(tǒng)的解決方法是，在進(jìn)行水下地形測量的同時，同步進(jìn)行水位觀測，以獲取水位面高程。當(dāng)測區(qū)的水位面隨時間（或距離）變化較大時，要定期（或定距離）的進(jìn)行水位觀測，并利用觀測所得時間（或距離）與潮位的對應(yīng)關(guān)系，內(nèi)插出每一時刻（或每一位置）的水位面高程；若水位變化微小或基本不發(fā)生變化，無須內(nèi)插，僅測量一個水位面高程即可。

參考文獻(xiàn)：

[1] 趙建虎，張紅梅.水下地形測量技術(shù)探討 [J].測繪信息與工程，1999，88（4）：22-26.

[2]梁開龍.水下地形測量 [M].北京：測繪出版社，1995.

[3]趙建虎，周豐年，張紅梅.水上GPS水位測量方法研究 [J].測繪通報(bào)，2001：1-3.

[4]趙建虎，劉經(jīng)南，等.GPS測定船體姿態(tài)方法研究 [J].武漢，武漢測繪科技大學(xué)學(xué)報(bào)，2000，25（4）：353-357.