程啟亮 張 薔

（長治市水文水資源勘測站）

隨著水文信息自動化不斷發(fā)展，許多新的測量儀器配備應用于基本水文站，極大緩解人力成本。傳統(tǒng)的水庫水下地形測量方法是人工測量水深，推求河底高程，從而計算庫容曲線[1]。采用ADCP 采集庫底高程數(shù)據(jù)，經(jīng)過處理之后，使用ArcGIS 建立數(shù)字高程模型，編寫腳本文件計算固定水位分級相應庫容是一種新興計算方法。這一方法可以大量的節(jié)約計算時間，提高計算結(jié)果的精度。本文采用ADCP 對漳澤庫區(qū)進行水下地形測量，然后基于ArcGIS 軟件構(gòu)建了該庫數(shù)字高程模型。

1 水庫基本情況

漳澤水庫樞紐位于濁漳河南源干流上，控制流域面積3 176 km2，上游年均來水量為2.25 億m3，水庫總庫容為4.27 億m3[2]，平均水深為5.49 m，為淺水水庫[2]（見圖1）。水庫年內(nèi)降水不均且變化較大，多年平均降水量為573.3 mm，多年平均氣溫為8.6 ℃。漳澤水庫是長治市重要水源地，是一座以工業(yè)、農(nóng)業(yè)灌溉、防洪為主，兼顧養(yǎng)殖和旅游等綜合利用的大型水庫。近年來，隨著國家對環(huán)保、水污染治理的高度重視，漳澤水庫水質(zhì)污染情況得到有效緩解。依據(jù)山西省水環(huán)境監(jiān)測中心長治分中心2019年、2020年對漳澤水庫的長期監(jiān)測資料分析，水質(zhì)整體分別達到為Ⅲ、Ⅳ類水質(zhì)。

圖1 漳澤水庫位置及上游水利工程分布圖

2 水庫水下地形測量

2.1 數(shù)據(jù)采集

在施測水域測量時，可先通過預測或查閱相關(guān)資料得到最大水深，然后使用ADCP 測量模式。ADCP 在測量中采用的是自動采集模式，設(shè)置單元層厚度為0.10 m。工作前，利用帶有CORS 功能的GPS RTK 將測量中的某一垂線的左右控制點坐標輸入到流動站手薄中，以兩點建立一條直線，進行放樣操作，允許設(shè)置偏離直線左右距離0.3 m～0.8 m，并將實時動態(tài)顯示在手薄屏幕上，施測人員根據(jù)偏移直線距離隨時校正航向，實現(xiàn)精準導航，確保沿垂線航行。在測量啟動前，將每一條垂線用GPS RTK流動站精確定位起始水邊，記錄并在計算機終端輸入換能器與水邊之間距離，待到達對岸后，用同樣方法記錄并輸入相應距離。整個斷面的水深和垂線位置等數(shù)據(jù)，通過數(shù)傳電臺傳輸?shù)接嬎銠C，系統(tǒng)軟件連續(xù)實時采集、記錄、儲存。測量結(jié)束后，利用WinRiver2 軟件以ASCII 文本格式導出需要的水深和距離數(shù)據(jù)，然后對數(shù)據(jù)進行處理，最終形成各個垂線的Excel 表。

考慮到近岸邊的水淺且實際中水草、亂石等多因素的影響，ADCP 換能器測量精度很難有效保障，需采用激光測距儀，測量漳澤水庫水面邊界點的坐標，得到水面邊界線的輪廓。

2.2 坐標轉(zhuǎn)換

坐標系統(tǒng)轉(zhuǎn)換是將PTK 實際測量的WGS84 大地坐標轉(zhuǎn)換為西安80 大地坐標。轉(zhuǎn)換計算需三個步驟，其數(shù)學模型是：

WGS84 大地坐標轉(zhuǎn)換空間直角坐標：

式中：N—橢球面卯酉圈的曲率半徑；e—橢球的第一偏心率；a、b—橢球的長短半徑；B、L、H—WGS84大地坐標；X、Y、Z—WGS84 空間平面直角坐標。

WGS84 坐標系轉(zhuǎn)換西安80 坐標系：

其中：

式中：R—旋轉(zhuǎn)矩形；ω—旋轉(zhuǎn)參數(shù)；m—尺度參數(shù)；X、Y、Z—空間平面直角坐標。

西安80 空間直角坐標轉(zhuǎn)換大地坐標：

式中：N—橢球面卯酉圈的曲率半徑；e—為橢球的第一偏心率；B、L、H—西安80 大地坐標；X、Y、Z—西安80 空間平面直角坐標。

本次共29 個控制點，采用RTK 進行放樣測量，測量使用GPS 雙頻接收機，型號為萊卡GS14。然后利用MAPINFO 軟件在庫區(qū)內(nèi)布置測深垂線，求得左右控制點坐標。將地理坐標輸入到GPS RTK 中，使用GPS RTK 定位控制點，通過左右控制點，連成一條直線，確保ADCP 在該垂線上測量。每一垂線測量啟動前，用GPS RTK 流動站定位起始水邊，并記錄水邊的距離，到達對岸后，同樣用GPS RTK 流動站定位結(jié)束水邊，并記錄到水邊的距離，從而采集垂線斷面的水深和位置等數(shù)據(jù)[2]。由于ADCP 測量需要一定的入水水深，近岸邊的水深無法測量，需采用激光測距儀，測量漳澤水庫水面邊界點的坐標，得到水面邊界線的輪廓。

本次漳澤水庫水下地形測量共施測59 條垂線，以及水面邊界線，如圖2 所示。數(shù)據(jù)成果包括垂線上各點的坐標經(jīng)緯度以及高程。數(shù)據(jù)結(jié)果以Excel格式存儲。

圖2 漳澤水庫水下地形垂線分布圖

3 庫容分析

3.1 降雨分析

收集研究區(qū)1956-2018年19 個雨量站的降雨資料，并計算面平均雨量。1965-1985年系列降水均值為587.1 mm，最大降水為911.2 mm，最小降水為367.1 mm；1985-2004年系列降水均值為571.0 mm，最大降水為993.1 mm，最小降水為343.9 mm；2004-2016年降水量平均值571.3 mm，其中，最大值713.5 mm，最小值446.5 mm。

3.2 徑流分析

利用漳澤水庫東大關(guān)、高河、西蓮、湛上等進庫站徑流數(shù)據(jù)，分析水庫入庫徑流，實測數(shù)據(jù)均從1992年開始收集整理。1965—1992年入庫水量計算采用下式計算：

式中：W入庫—入庫水量；W出庫—出庫水量（包括灌溉、工業(yè)供水和棄水）；W蒸發(fā)—庫區(qū)蒸發(fā)損失；W滲漏—庫區(qū)滲漏損失；ΔV—當月庫蓄變量。

庫區(qū)的蒸發(fā)損失按蒸發(fā)深和水庫水面面積計算，庫區(qū)滲漏量按水庫月蓄水量的1%計，庫容曲線采用相應年份實測或內(nèi)插的庫容曲線。1993—2018年直接采用四個入庫站資料。1965年至1985年平均徑流量為18 405 萬m3，1985年至2004年平均徑流量10 578 萬m3，2004年至2018年平均徑流量12 154 萬m3。

3.3 庫容計算

采用以下計算庫容：

式中：d—DEM 格網(wǎng)間距；hi第i 個格網(wǎng)的高程；H—水位高程；n—水位高程H 以下的格網(wǎng)總數(shù)；Vi—第i個格網(wǎng)水位為H 的水面間四棱柱的體積；V—水位H 對應的庫容；S—水位H 對應的面積。

采用AICGIS-3D Analyst 工具中的功能性表面體積工具和生成的高程數(shù)字模型計算漳澤水庫庫容。本次計算水位分級為0.2 m，利用腳本文件計算可以得到每個水位分級對應下的庫容值，庫容曲線如圖3。

圖3 漳澤水庫庫容曲線

漳澤水庫自建庫以來，1985-2004年施測過水下地形，將2018年施測的水下地形數(shù)據(jù)與歷史數(shù)據(jù)進行對比（見圖4）。

圖4 漳澤水庫水位與庫容

1965-1985年，漳澤水庫總淤積量為1 693 萬m3，漳澤水庫平均淤積量為84.7 萬m3/年。1985-2004年，漳澤總淤積量為500 萬m3，漳澤水庫平均淤積量為26.3 萬m3/年。2004-2018年，漳澤水庫總淤積量為396 萬m3，漳澤水庫平均淤積量為28.3 萬m3/年。

從圖4 中可以看出，隨著水位的增加，淤積量增大，達到一定水位后淤積量減少，所以，2018年與歷年曲線間距為中間胖兩邊瘦。從水庫的淤積形態(tài)可以看出，漳澤水庫屬于三角洲形態(tài)[3]。

3.4 合理性分析

收集降水、徑流量資料，通過水沙相關(guān)關(guān)系，分析水庫淤積情況，驗證庫容曲線的合理性。水流中的含沙量主要是受暴雨影響，降水量越大，徑流量越大，含沙量則越大。1965-1985年間，降水量比較大，徑流量相應較大，徑流泥沙含量多，入庫的泥沙含量高。1986-2004年與2005-2018年兩系列比較，2005-2018年入庫徑流量均值較大。因此，1965-1985年入庫含沙量最大，2005-2018年入庫含沙量居中，1986-2004年入庫含沙量最小。經(jīng)計算，1965-1985年平均水庫淤積為84.7 萬m3，1986-2004年平均水庫淤積為26.3 萬m3，2005-2018年平均水庫淤積為28.3 萬m3，基本符合降水-徑流-產(chǎn)沙規(guī)律。

4 結(jié)論

本文采用ADCP 系統(tǒng)進行漳澤水庫水下地形測量，并將測量結(jié)果導入ARCGIS 軟件中進行處理，構(gòu)建水庫數(shù)值高程模型，計算出水位分級為0.2 m 的庫容曲線，得到1965-2018年期間漳澤水庫最大水位為900.86 m，相應庫容8 021 萬m3，2004年庫容總淤積量為396 萬m3，年平均水庫淤積量為28.3 萬m3。