趙新華，王 沛，余華芬，陳連芳，姜安春

（1.國(guó)網(wǎng)新源水電有限公司新安江水力發(fā)電廠，浙江杭州，311608；2.浙江省第二測(cè)繪院，浙江杭州，310012）

基于航空攝影測(cè)量的新安江水庫(kù)水域面積及庫(kù)容變化分析

趙新華1，王 沛2，余華芬2，陳連芳1，姜安春2

（1.國(guó)網(wǎng)新源水電有限公司新安江水力發(fā)電廠，浙江杭州，311608；2.浙江省第二測(cè)繪院，浙江杭州，310012）

采用航空攝影測(cè)量技術(shù)獲得2001年和2011年新安江水庫(kù)攝影水位以上的地形數(shù)據(jù)，通過(guò)構(gòu)建DEM計(jì)算新安江水庫(kù)水域面積和庫(kù)容。從平面和立體兩個(gè)不同的角度，分析了新安江水庫(kù)2001年和2011年水域面積的變化與相應(yīng)區(qū)域的庫(kù)容變化。從平面上分析了水域面積變化，從立體上分析了庫(kù)容變化，對(duì)水域和庫(kù)容的變化都進(jìn)行了定量分析與定位分析，校核了新安江水庫(kù)現(xiàn)有的水位-庫(kù)容曲線模型，為新安江水庫(kù)調(diào)度提供了參考。

新安江水庫(kù)；航空攝影測(cè)量；庫(kù)容；水域面積；DEM

0 引言

水庫(kù)是蓄水發(fā)電、灌溉和防洪調(diào)度等的重要設(shè)施，水庫(kù)庫(kù)容作為一個(gè)非常重要的參數(shù)，其計(jì)算結(jié)果的精度高低直接影響水庫(kù)的防洪安全及發(fā)電、灌溉等效益。水庫(kù)建成后多年運(yùn)行造成庫(kù)區(qū)的淤積、庫(kù)岸人類活動(dòng)以及水流的沖刷，改變了原有的水位－庫(kù)容關(guān)系，岸線及蓄水位也產(chǎn)生了一定的變化。因此，有必要對(duì)水庫(kù)的水域面積、水位－庫(kù)容關(guān)系進(jìn)行復(fù)核，掌握水庫(kù)準(zhǔn)確的庫(kù)容變化，在水庫(kù)、大壩安全的前提下，最大限度地發(fā)揮水庫(kù)的經(jīng)濟(jì)、社會(huì)與生態(tài)效益[1]。

已有研究對(duì)水庫(kù)水域面積變化的分析較多[2-7]，但數(shù)據(jù)源多為L(zhǎng)andsat或者SPOT影像，空間分辨率不高，精度不能保證，且針對(duì)新安江水庫(kù)水域面積變化的研究目前較少。對(duì)于庫(kù)容計(jì)算，多為量算某一期的庫(kù)容曲線，而沒有對(duì)多年份的庫(kù)容變化進(jìn)行比較。

水庫(kù)庫(kù)容計(jì)算常用方法有地形圖量算法、水下地形測(cè)量法、數(shù)字高程模型（DEM）法等[8-13]。地形圖量算方法受地形圖比例尺、量測(cè)儀器及計(jì)算方法的影響，精度較低；水下地形測(cè)量法工作量大、成本投入高，數(shù)據(jù)量會(huì)隨水域面積的增加而急劇增大，并且對(duì)高于當(dāng)前水面的庫(kù)容曲線無(wú)法進(jìn)行計(jì)算和繪制[14]。數(shù)字高程模型（DEM）法計(jì)算庫(kù)容具有精度高、運(yùn)算速度快等特點(diǎn)，在庫(kù)容計(jì)算中更具優(yōu)勢(shì)。

本研究采用航空攝影測(cè)量技術(shù)獲得2001年和2011年新安江水庫(kù)攝影水位以上的地形數(shù)據(jù)，從平面和立體兩個(gè)不同的角度，分析了新安江水庫(kù)2001～2011年水域面積的變化與相應(yīng)區(qū)域的庫(kù)容變化。從平面上分析了水域面積變化，從立體上分析了庫(kù)容變化，對(duì)水域和庫(kù)容的變化都進(jìn)行了定量分析與定位分析，校核了新安江水庫(kù)現(xiàn)有的水位-庫(kù)容曲線模型，提高新安江水庫(kù)調(diào)度的效率。

1 研究區(qū)概況及數(shù)據(jù)源

新安江水庫(kù)位于錢塘江上游主干流新安江的下游，水庫(kù)總庫(kù)容216.26億m3，為多年調(diào)節(jié)水庫(kù)，設(shè)計(jì)、校核洪水標(biāo)準(zhǔn)分別為重現(xiàn)期1 000年和重現(xiàn)期10 000年加安全保證值。水庫(kù)正常蓄水位108.0 m，設(shè)計(jì)洪水位111.0m，校核洪水位114.0m，汛期防洪限制水位為106.5 m，死水位86.0 m。其中，安徽省境內(nèi)的新安江水庫(kù)108 m水位的水域面積為18.65 km2，占新安江水庫(kù)總體水域面積的3.29%，見圖1。

本研究所用的數(shù)據(jù)來(lái)源和資料包括：新安江水庫(kù)地區(qū)2001年和2011年的航空攝影資料（共涉及1∶10 000圖77幅）、相應(yīng)的1∶10 000數(shù)字線劃圖（DLG）產(chǎn)品、數(shù)字正射影像圖（DOM）產(chǎn)品、第一次國(guó)家地理國(guó)情普查地表覆蓋數(shù)據(jù)等。

圖1 新安江水庫(kù)區(qū)位圖Fig.1 Location of Xin'anjiang reservoir

2 研究方法

目前數(shù)字高程模型（DEM）的兩種主要結(jié)構(gòu)有基于規(guī)則格網(wǎng)的DEM和基于不規(guī)則三角網(wǎng)的DEM，為便于存儲(chǔ)、處理和建模，文中DEM構(gòu)建采用基于矢量數(shù)據(jù)構(gòu)建規(guī)則格網(wǎng)DEM的方法。

方格網(wǎng)法是利用已建立的庫(kù)區(qū)數(shù)字高程模型（DEM），將水庫(kù)水體微分成若干個(gè)正方體，通過(guò)對(duì)每個(gè)正方體的體積空間積分求和，即可求得整個(gè)水庫(kù)的庫(kù)容。其數(shù)學(xué)模型為：

式中：V為庫(kù)容；Ps為單個(gè)格網(wǎng)的面積；H為指定水位的高程值；hi為高程小于指定水位的格網(wǎng)高程值；n為高程小于H的DEM格網(wǎng)的個(gè)數(shù)。

采用數(shù)字?jǐn)z影測(cè)量技術(shù)獲取2001年、2011年攝影水位以上的地形數(shù)據(jù)，通過(guò)構(gòu)建DEM，利用方格網(wǎng)法對(duì)100～115 m高程范圍之間的新安江水庫(kù)水域面積和庫(kù)容進(jìn)行計(jì)算，并分析2001～2011年的變化情況。

3 結(jié)果與分析

3.1 水域變化

根據(jù)新安江水庫(kù)的設(shè)計(jì)值，正常高水位是108 m，繪制該水位線對(duì)應(yīng)的水域分布圖，如圖2所示。

新安江水庫(kù)變化區(qū)域主要分布在新登-淳安的一級(jí)公路沿線以及千島湖鎮(zhèn)的環(huán)湖北路、睦州大道、珍珠大道沿線等地。

2001年和2011年新安江水庫(kù)水位-水域面積曲線如圖3所示，并與新安江水庫(kù)的設(shè)計(jì)值進(jìn)行對(duì)比。

圖2 新安江水庫(kù)108 m水位對(duì)應(yīng)的水域變化范圍Fig.2 Corresponding change scope of water area as the water level is 108 m

圖3 新安江水庫(kù)水位-水域面積曲線Fig.3 Water level-area curve of Xin'anjiang reservoir

由圖3可知，2001～2011年，新安江水庫(kù)的水域面積呈現(xiàn)減小的趨勢(shì)，均小于水庫(kù)設(shè)計(jì)值且不同水位對(duì)應(yīng)的水域面積減小幅度不同。以108 m正常高水位線的水域面積為例，2001年水域面積為564.85 km2，2011年下降至555.51 km2，108 m水位的水域面積值在10年間減小了約9.34 km2。

對(duì)浙江省境內(nèi)各鄉(xiāng)級(jí)行政區(qū)范圍內(nèi)的水域變化進(jìn)行分析，2001～2011年，新安江水庫(kù)水域面積變化較大的鄉(xiāng)級(jí)行政區(qū)有千島湖鎮(zhèn)、文昌鎮(zhèn)、汾口鎮(zhèn)等，其中千島湖鎮(zhèn)變化最大，在108 m水位處水域面積減小了394.58 ha。水域面積基本沒有發(fā)生變化的鄉(xiāng)級(jí)行政區(qū)有大墅鎮(zhèn)和宋村鄉(xiāng)。水域面積變化較小的鄉(xiāng)級(jí)行政區(qū)有新安江街道、安陽(yáng)鄉(xiāng)、浪川鄉(xiāng)等，其中新安江街道在108 m水位處的水域面積減小了0.86 ha。

參照第一次國(guó)家地理國(guó)情普查地表覆蓋數(shù)據(jù)，對(duì)2001～2011年變化區(qū)域進(jìn)行地類分析，見表1。

表1 變化區(qū)域現(xiàn)狀地類分析Table 1 Current land types in the change area

由表1可知，變化區(qū)域中現(xiàn)狀地類為建設(shè)用地447.68 ha，占變化區(qū)域總面積的52.02%，是變化區(qū)域中最多的現(xiàn)狀地類，主要分布在淳安縣千島湖鎮(zhèn)；其次是耕地，面積151.98 ha，占變化區(qū)域總面積的17.66%，主要分布在淳安縣汾口鎮(zhèn)；裸地最少，面積為3.27 ha，占變化總面積的0.38%，主要分布在千島湖鎮(zhèn)。

3.2 庫(kù)容變化

本研究中水庫(kù)庫(kù)容以100 m等高線開始計(jì)算，即水庫(kù)水位在100 m時(shí)庫(kù)容為0，計(jì)算得到的水庫(kù)庫(kù)容不是庫(kù)容的絕對(duì)值，而是相對(duì)庫(kù)容。圖4為特征水位對(duì)應(yīng)的相對(duì)庫(kù)容。

由圖4可知，2001～2011年，新安江水庫(kù)的庫(kù)容總體是減小的，與水域面積變化趨勢(shì)一致。以108 m水位為例，2001年的相對(duì)庫(kù)容為41.49億m3，到2011年，相對(duì)庫(kù)容減小為41.06億m3，10年間減小了0.43億m3。

圖4 新安江水庫(kù)2001年、2011年特征水位對(duì)應(yīng)的相對(duì)庫(kù)容Fig.4 Corresponding relative storage capacity of the character?istic water level in 2001 and 2011 of Xin'anjiang reservoir

根據(jù)庫(kù)容計(jì)算數(shù)據(jù)分別構(gòu)建新安江水庫(kù)設(shè)計(jì)值、2001年和2011年的水位-相對(duì)庫(kù)容曲線，如圖5所示。

圖5 設(shè)計(jì)值、2001年、2011年的水位-相對(duì)庫(kù)容曲線Fig.5 Curve of water level and relative storage capacity in de?sign and in the year 2001 and 2011

依據(jù)最新的水位-庫(kù)容曲線，可以有效提高安排灌溉用水、發(fā)電用水、防洪調(diào)度等日常水庫(kù)管理工作的精度。

3.3 變化原因分析

水庫(kù)水域面積和庫(kù)容變化的原因可以是多方面的，通?？梢詺w結(jié)為兩大類：自然因素和人為因素。自然因素，通常是水庫(kù)沿岸崩塌或泥沙自然淤積造成水域面積和庫(kù)容減小，而本研究以相對(duì)庫(kù)容來(lái)對(duì)新安江水庫(kù)的庫(kù)容進(jìn)行分析，泥沙淤積對(duì)相對(duì)庫(kù)容的計(jì)算影響不大，故忽略此因素。對(duì)于人為因素，本研究將分析得到的新安江水庫(kù)水域變化區(qū)域與獲得的數(shù)字正射影像（DOM）數(shù)據(jù)進(jìn)行對(duì)比，從而分析新安江水庫(kù)水域及庫(kù)容變化的具體原因。

圖6是新安江水庫(kù)區(qū)域中的淳安縣千島湖鎮(zhèn)珍珠半島2001～2011年的變化區(qū)域，(a)、(b)分別為2001年、2011年的數(shù)字正射影像。對(duì)比可知，變化區(qū)域由2001年的水體變?yōu)?011的陸地。相應(yīng)區(qū)域2001年和2011年的DEM如圖7所示。

圖6 淳安縣千島湖鎮(zhèn)珍珠半島2001年和2011年水域變化對(duì)比Fig.6 Comparison of water area at Pearl peninsula in the year 2001 and 2011

圖7 千島湖鎮(zhèn)珍珠半島2001～2011年水庫(kù)變化區(qū)域DEM對(duì)比Fig.7 Comparison of DEM in the reservoir change area at Pearl peninsula in the year 2001 and 2011

人工填湖造陸是新安江水庫(kù)水域面積和庫(kù)容減小的主要原因。2001～2011年填湖造陸引起的108 m水位的水域面積減少9.34 km2，庫(kù)容減少量為0.43億m3。

4 結(jié)語(yǔ)

新安江水庫(kù)2001～2011年正常高水位108 m的水域面積減小了9.34km2，對(duì)應(yīng)庫(kù)容減小了0.43億m3，100～115m水位范圍的庫(kù)容量減小了1.02億m3。

綜合上述分析，本研究有以下結(jié)論：

（1）本研究探索性地將航空攝影測(cè)量技術(shù)應(yīng)用于特大型水庫(kù)的庫(kù)容測(cè)量與分析，精確地測(cè)量出“攝影水位”以上的庫(kù)區(qū)地形，填補(bǔ)了現(xiàn)有地形圖中“水庫(kù)高水位”與“攝影水位”之間的地形。與利用現(xiàn)有地形圖、遙感等方法計(jì)算庫(kù)容相比，本方法具有更高的精度。

（2）2001～2011年，新安江水庫(kù)的水域面積和庫(kù)容均呈現(xiàn)減小的趨勢(shì)，且基本上小于水庫(kù)原設(shè)計(jì)值。新安江水庫(kù)的水域面積和庫(kù)容不斷減小，主要是修建高級(jí)度假酒店和別墅項(xiàng)目、各種經(jīng)濟(jì)開發(fā)區(qū)、修建公路等人工填湖造陸活動(dòng)引起的。 ■

[1]丁志雄,顏廷松,屈吉鴻.多源遙感影像在水庫(kù)水位-庫(kù)容曲線復(fù)核中的應(yīng)用[J].華北水利水電學(xué)院學(xué)報(bào),2012,4 (33):32-35.

[2]曾忠平,盧新海.城市湖泊時(shí)空演變的遙感分析——以武漢市為例[J].湖泊科學(xué),2008,20(5):648-654.

[3]李景剛,李紀(jì)人,黃詩(shī)峰,等.Terra/MODIS時(shí)間序列數(shù)據(jù)在湖泊水域面積動(dòng)態(tài)監(jiān)測(cè)中的應(yīng)用研究[J].自然資源學(xué)報(bào), 2009,5(24):923-933.

[4]黃浦江,劉艷芳,劉暢,等.基于RS與GIS的武漢城市湖泊鹽化研究[J].生態(tài)環(huán)境學(xué)報(bào),2012,21(9):1588-1593.

[5]吳川,張玉龍,許秀貞,等.基于Landsat TM/ETM和HJ-1 A/B影像的丹江口水庫(kù)水域變化監(jiān)測(cè)研究[J].長(zhǎng)江流域資源與環(huán)境,2013,9(22):1207-1213.

[6]徐娜,賈建華,羅菊花,等.江蘇省湖泊遙感監(jiān)測(cè)及10年動(dòng)態(tài)變化分析[J].長(zhǎng)江流域資源與環(huán)境,2014,4(23):468-474.

[7]張毅,陳成忠,吳桂平,等.遙感影像空間分辨率變化對(duì)湖泊水體提取精度的影響[J].湖泊科學(xué),2015,27(2):335-342.

[8]楊中華,陳琳.基于GIS的水庫(kù)庫(kù)容測(cè)量方法的研究與實(shí)踐[J].測(cè)繪通報(bào),2002(11):28-30.

[9]馬文祥.基于GIS的中小水庫(kù)地形和庫(kù)容測(cè)量實(shí)踐[J].測(cè)繪通報(bào),2006(8):66-68.

[10]張建新,程琳,王光謙,等.堰塞湖最大庫(kù)容及庫(kù)容曲線分析計(jì)算[J].水文,2009,5(29):63-66.

[11]盧云峰,譚德寶,楊中華.基于空間信息技術(shù)的大型水庫(kù)庫(kù)容計(jì)算方法[J].長(zhǎng)江科學(xué)院院報(bào),2010,1(27):9-12.

[12]吳恒友.貴陽(yáng)市阿哈水庫(kù)庫(kù)容計(jì)算方法的探討[J].測(cè)繪工程,2010,2(19):78-80.

[13]楊伯鋼,張勁松,霍滿軍.城市湖泊庫(kù)容勘測(cè)方法對(duì)比研究[J].測(cè)繪通報(bào),2014(11):50-52.

[14]張紅梅,趙建虎.水庫(kù)庫(kù)容和淤積測(cè)量技術(shù)研究[J].水利學(xué)報(bào),2012(12):33-37.

作者郵箱：zhaoxinhua1@163.com

Analysis on change of water area and storage capacity of Xin′anjiang reservoir based on aerial photogrammetry

by ZHAO Xin-hua,WANG Pei,YU Hua-fen,CHEN Lian-fang and JIANG An-chun
Xin′anjiang Hydropower Plant

In this paper,water area and storage capacity of Xin′anjiang reservoir were calculated using digital elevation model which was constructed based on terrain data of Xin′anjiang reservoir(2001 and 2011).Then,the changes of water area and corresponding storage capacity of Xin′anjiang reservoir were analyzed quantitatively in the year of 2001 and 2011.Finally,water level-capacity curve model of Xin′anjiang reservoir was checked,which could provide a new reference for the regulation of Xin′anji?ang reservoir.

Xin′anjiang reservoir;aerial photogrammetry;storage capacity;water area;DEM

P23

B

1671-1092（2016）05-0046-04

2015-10-10；

2015-11-16

趙新華（1968-），男，浙江省文成縣人，工程碩士，高級(jí)工程師，主要從事水庫(kù)調(diào)度和大壩安全管理工作。