水庫淤積測量技術分析論文

1、水庫淤積測量技術分析論文近年來，我國洪水災害頻繁發(fā)生，給國民經濟和人民生命財產帶來嚴重威脅和重大損失。水庫作為人類蓄水發(fā)電、灌溉和防洪調度等的重要設施，發(fā)揮著越來越大的作用，并取得了巨大的社會效益和經濟效益。水庫庫容和淤積量是水庫調度的重要參數，其精度直接影響到水庫的防洪安全與蓄水興利。隨著現(xiàn)代測控技術的迅速發(fā)展，將其應用到水利電力領域是一項值得研究和推廣的課題。本文在分析傳統(tǒng)水庫庫容及淤積量測量所存在問題的基礎上，依靠高精度GPS(GlobalPositioningSystem,簡稱GPS)定位和回聲測深技術，對湖南某庫區(qū)水下地形進行了測量，并提出了根據三角形構網方法，利用“三角柱”的水柱體

2、積和淤積體積來獲得庫容和淤積量的新見解，經實際運用，取得了滿意效果。 1常規(guī)庫容及淤積量的確定 常規(guī)的庫容計算方法多采用斷面法。其庫區(qū)容量的計算模型為： (1) 式中：Vi、Li為第i個斷面到第i+1個斷面間的庫容和距離；n為分段個數；Si、m、d、hi分別為第i個斷面的面積、測點個數、點間距和每個測點的深度測量值。 采用斷面法雖然操作簡單，但受前提假設的制約，精度難以保證。淤積量是根據前后兩次的的庫容較差獲得，庫容不準確，淤積量的計算精度就無從談起。 2高精度水下地形測量技術1,2 2.1水下地形測量所謂水下地形測量，就是利用測量儀器來確定水底點的三維坐標的過程。隨著GPS技術的迅速發(fā)展，水

3、下地形測量方法取得了很大的進展。目前，水下地形測量技術已定型于采用GPS獲取平面坐標，測深儀獲取深度數據的基本模式。同時，為了獲得水下地物的海拔高程，以及消除潮汐、水位落差等諸因素的影響，進行水位監(jiān)測也是一個重要環(huán)節(jié)。水下地形測量現(xiàn)狀示意圖如下。針對上述現(xiàn)有模式，文獻1對測量設備的選型、基準點的布設、潮汐觀測的具體實施等諸多技術問題，進行了深入探討。 2.2GPS載波相位差分定位技術和回聲測深技術隨著GPS技術的發(fā)展，GPS日益廣泛應用于水利電力工程的各個方面。為了提高定位精度，一般均采用差分技術。在眾多的差分技術中，偽距差分和載波相位差分是最為常用的兩種測量模式，后者的定位精度較高(厘米級)

4、，通常用于高精度的測量工程和研究中。 圖1水下地形測量現(xiàn)狀示意 載波相位差分測量的定位精度很大程度上依賴于整周模糊度能否在航精確確定。整周模糊度在航解算(OTF)是一種動態(tài)環(huán)境下的模糊度確定方法，它可省去在精密動態(tài)定位中的的靜態(tài)初始化過程。常規(guī)精密定位中復雜的整周跳變問題也因OTF的引入變得十分簡單。載波相位差分測量整周模糊度的確定模型為： Xk=k-1Xk-1+k-1Wk-1WkN(0,Qk) Zk,=Hk,Xk+Vk,Vk,N(0,Rk,) (2) 式中：Xk=(dxdydzxyzdn0dn1dnm)為狀態(tài)向量；k-1為狀態(tài)轉移矩陣；Hk,為載波相位的測量矩陣；Rk,為載波相位的方差陣；Q

5、k為系數陣。 =CCTQk=ffT=minf=CT(DN-D) (3) 由式(3)計算得到整周模糊度N后，代入載波相位觀測方程，便可以獲得厘米級甚至毫米級的平面定位精度。 回聲測深儀是一種單波束測深設備，深度的測量是根據最小聲程決定。按照使用頻率個數的不同，又可分為單頻和雙頻。雙頻測深儀根據兩個頻率測量深度較差獲得淤積層厚度。 2.3高精度庫容和淤積量測量方法庫容和淤積量的精密測量采用現(xiàn)代水下地形測量方法，即利用GPS載波相位差分測量技術進行平面定位，測深儀進行深度測量，GPS和測深儀保證同步作業(yè)，獲取水底測點平面和深度信息的作業(yè)模式。 為了保證庫容和淤積量的計算精度，需要對庫區(qū)進行測線設計，

6、GPS和測深采樣也要按照水下地形測量規(guī)范等間隔或等時間采樣。設測量比例尺為1:Scale，測量船的平均速度為，則測線間距d和時間間隔t為： d=Scale10-4 t=d/ (4) 為了提高測量精度，在測線布設時，還應該考慮水下地形的變化趨勢，若變化相對比較平坦，則測線間距可以適當放寬，否則，需加密測線。這有利于使測點均勻分布于整個測區(qū)，同時在測區(qū)水下地形變化復雜的地區(qū)使測點深度或高程能更好地反映水下地形的真實面貌。 3庫容和淤積量的計算方法 3.1庫容計算方法為了提高計算精度，充分利用水下地形測量數據，本文提出了一種三角柱計算庫容的方法。該法建立在實際測點的基礎上，根據圖2，相鄰三個測點可構

7、成的三角柱體積為： (5) 設n為整個庫區(qū)三角形個數，則整個庫區(qū)的庫容為： (6) 3.2淤積量計算方法當庫底為基巖構造時，采用雙頻(f1、f2)測深儀測深，淤積量的計算方法同庫容相似。設相鄰三個測點在淤積層表面利用f1測得深度分別為h1、h2和h3，利用f2，在對應點基巖上測得的深度分別為h1、h2和h3，若設淤積層表面面積為S2，基巖上的面積為S3，則淤積量為： (7) 式中淤積層上下面的面積S2和S3的計算方法同式(5)。 則庫區(qū)的淤積量為： (8) 對于上述情況，V的計算還可采用模型： V=Vf2-Vf1 (9) 式中：Vf1、Vf2分別代表根據f1、f2測得的淤積表面、基巖表面上的深

8、度計算得到的體積。 圖2相鄰三個測點構成的三角柱示意 圖3相鄰三個測點構成的淤積三角柱示意 然而，對于淤積層下地質是非基巖的情況，式(9)的庫區(qū)淤積量計算模型就不再適用，而需要根據建庫初期的原始床面(地形圖)計算空庫容，或者前期確定的庫容量，與根據本次利用f1頻率測量的水深(淤積層表面的水深)計算所得庫容Vf1較差得到實際的庫區(qū)淤積量。其計算模型為： V=Vf1-V前期庫容 (10) 4問題討論 本文所提出的基于現(xiàn)代水下地形測量技術的水庫庫容和淤積量確定方法相對于傳統(tǒng)的斷面法具有許多優(yōu)點，諸如定位精度高、計算結果準確、所得數據可用于水下地形圖的繪制及DTM的建立等。然而，相對于傳統(tǒng)庫容和淤積量

9、的確定方法，由于采用了先進的測控設備，無疑會增大測量和計算方法上的復雜度，但這些是可以通過計算機編程來自動化實現(xiàn)的?，F(xiàn)將上述方法在實際數據處理中的幾個難點加以討論。 (1)對于比較大的庫區(qū)，如江河形成的自然庫區(qū)，數據量會隨水域面積的增加而急劇增大。在利用這些數據構造庫區(qū)三角形時會因存儲量和搜索范圍過大，占用過多的計算機內存，可能會導致計算速度過慢或者死機。為克服這一問題，在三角形構造中可采用一種快速的三角形構網方法，即局域搜索法。根據測區(qū)范圍和測點的數量，可事先對整個區(qū)域根據坐標進行劃分，然后在結合拓展三角形的范圍索引各個分割區(qū)，在小區(qū)域內實現(xiàn)快速搜索。這樣可以大大的節(jié)約計算機內存，提高三角形

10、的構網速度。 (2)通過水下地形測量可給出水面以下的深度，以及根據水面下的實測結果計算水底到水面高程變化的庫容曲線，而對于高于當前水面的水位面庫容曲線無法進行計算和繪制。為了得到一個全面反映庫區(qū)容量變化的庫容曲線，需要將庫區(qū)邊緣數字高程信息引入庫容計算中。庫區(qū)邊緣陸地的數字高程信息可通過兩種途徑獲得。一種是利用GPS載波相位差分技術進行動態(tài)地形測量獲得；另一種方法是通過已有的地形圖或DTM獲得。若利用GPS載波相位差分測量技術獲得陸地數字信息，則GPS天線相位中心的平面位置即為陸地測點的平面位置，相位中心的高程減去天線高便是陸地高程。 (3)在(2)中，已有地形圖與現(xiàn)有測量成果共同用于庫容曲線

11、計算時，兩套資料的高程和坐標基準必須匹配。對于將水底點的深度轉換成高程問題，傳統(tǒng)的解決方法是，在進行水下地形測量的同時，同步進行水位觀測，以獲取水位面高程。當測區(qū)的水位面隨時間(或距離)變化較大時，要定期(或定距離)的進行水位觀測，并利用觀測所得時間(或距離)與潮位的對應關系，內插出每一時刻(或每一位置)的水位面高程；若水位變化微小或基本不發(fā)生變化，無須內插，僅測量一個水位面高程即可。根據文獻3和4，現(xiàn)代水下地形測量，省去了上述煩瑣的過程，直接利用GPSRTK技術獲得水底點高程。根據圖1和GPS載波相位測量技術，只要量取GPS天線相位中心到換能器之間的垂距hG-T得水位面的高程hsurface

12、，進而獲得水底點的高程hb。設h為測量的水深，GPS相位中心的高程為hGPS，則hb可表達為： hb=hsurface-h=(hGPS-hG-T)-h (11) 式(11)是在作業(yè)條件相對較好情況下計算水底點高程的模型，若由于波浪、船體的運動，上述條件很難滿足，為此，在實際測量中需要引入姿態(tài)測量的內容。姿態(tài)測量通常采用波浪補償儀或姿態(tài)儀，但由于儀器費用昂貴，這里引進GPS姿態(tài)測量技術。只需增加船載的2臺GPS接收機，使可獲得船體的姿態(tài)。儀器的架設如圖4。 圖4測姿GPS天線安放 根據文獻4，GPS測姿完全可以滿足IHO的精度要求。姿態(tài)參數(橫搖r、縱搖p、動態(tài)吃水ds)測定后，便可對式(11)

13、中的hG-T和h實施修正。設實際測量值分別為hG-T和h,則修正后的hG-t和h為： h=h-hr-hp-ds=h(cosp+cosr-1)-ds hG-T=hG-T-hG-T,r-hG-T,p=h-(cosp+cosr-1) (12) 這樣，利用式(11)和(12)便可獲得水底點的高程。這種方法無須進行水位改正，直接得到同陸地高程基準一致的高程。水下地形測量的平面坐標系統(tǒng)在測量時便可設置為同一系統(tǒng)；若不為同一系統(tǒng)，還要進行坐標轉換。 5實驗及結論 該方法在湖南某“水庫淤積測量及庫容曲線修正研究”課題中得到了應用和驗證。該水庫為山區(qū)的一個天然水庫(地質為巖石結構)，主要用于電廠的發(fā)電和蓄洪，水

14、庫面積約150km2。1998年，由武漢大學測繪學院(原武漢測繪科技大學地測學院)承擔了此項測量工作，總計測量34萬個測點。利用該方法計算僅花費不到2秒的時間，完成了庫容的計算，不同水位面庫容數據和庫容曲線計算結果如表1所示： 表1不同起算面的庫容量 高程起算面/m 庫容/m3 高程起算面/m 庫容/m3 84.00 86.00 88.00 90.00 92.00 94.00 96.00 98.00 100.00 102.00 0.00 236.83 2605.17 4989.83 13215.12 33995.46 116176.44 425711.85 1202286.67 2594917.95 104.00 106.00 108.00 110.00 112.00 114.00 116.00 118.00 120.00 4693450.36 7870823.68 12656827.05 18927380.47 26959658.83 37062902.96 48781192.88 62376448.42 77688490.85 將這種方法計算所得庫容曲線與結合已有淤積資料，并根據1963年所測得庫容來推算而得到的庫容進行比較，