基于多源遙感技術(shù)的水庫特征參數(shù)反演研究

李闖　李兵　王寧

摘要：以丹江口水庫為例，采用Sentinel-1A合成孔徑雷達(dá)（SAR）和Sentinel-2A/B多光譜成像儀（MSI）聯(lián)合提取水域面積，通過建立多源傳感器（SAR和MSI）所提取水域面積差異的定量函數(shù)關(guān)系，修正二者水域面積之間的差異，并將修正后的水域面積與Sentinel-3A雷達(dá)高度計（RA）所提取的水位相結(jié)合，開展水庫蓄水量反演研究，最后以“全國水雨情網(wǎng)”公布的實(shí)測數(shù)據(jù)等為參考，評估水庫特征參數(shù)的反演精度。結(jié)果表明：由RA數(shù)據(jù)提取的遙測水位與實(shí)測水位相比，R2達(dá)到0.99以上，RMSE達(dá)到0.093 3，遙測水位與實(shí)測水位的平均誤差僅有9.33 cm;相對于單一MSI影像或SAR影像所提取的水域面積，二者水域面積聯(lián)合參與建立的“遙測水位-多源遙感水域面積”模型的R2和RMSE均有所優(yōu)化，其中，R2分別提高了0.26%和0.60%，RMSE分別降低了2.63和1.29;由遙測水位和多源遙感水域面積結(jié)合水庫蓄水量模型反演的蓄水量與標(biāo)準(zhǔn)蓄水量相比，R2達(dá)到0.99，RMSE達(dá)到1.177，相對誤差平均為1.8%。研究結(jié)果表明：利用多源遙感技術(shù)開展水庫特征參數(shù)的反演研究，能夠更加準(zhǔn)確地反映水庫的運(yùn)行狀態(tài)。

關(guān)鍵詞：水庫特征參數(shù); 多源遙感; 雷達(dá)高度計; 多光譜成像儀; 合成孔徑雷達(dá)

中圖法分類號： TN958

文獻(xiàn)標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2021.11.015

0引 言

水庫是陸地水資源的聚集地，水庫特征參數(shù)主要包括水庫的水位、水域面積和蓄水量等。這些參數(shù)直接反映著水庫健康運(yùn)行的狀態(tài)，對區(qū)域的氣候、濕度、降水等有著重要的影響。研究水庫特征參數(shù)有助于掌握水庫的運(yùn)行狀態(tài)，對攔截洪水、保障供水、調(diào)配水資源等具有積極作用，因此，開展水庫特征參數(shù)的反演研究具有十分重要的現(xiàn)實(shí)意義[1]。

傳統(tǒng)方法主要依靠人力到現(xiàn)場布設(shè)各種監(jiān)測儀器并繪制庫區(qū)地圖等來反演水庫特征參數(shù);單一遙感方法主要利用一種遙感技術(shù)對水庫的一種參數(shù)進(jìn)行提取，其他參數(shù)仍采用現(xiàn)場調(diào)查的方式。李書杰[2]利用人力進(jìn)行外業(yè)測量獲得水庫的水位并繪制了庫區(qū)地圖，通過對內(nèi)業(yè)資料的整理，利用等高線、方格網(wǎng)等求得水庫的水域面積，受制圖過程中人為操作產(chǎn)生的精度損失影響，所反演的蓄水量只是粗略的結(jié)果。蘇業(yè)助等[3]為了科學(xué)化管理水庫，減少因外業(yè)測量和內(nèi)業(yè)資料整理而造成的數(shù)據(jù)偏差，根據(jù)蓄水量平衡原理修正了蓄水量反演曲線，提高了蓄水量反演精度。孫建蕓等[4]基于我國GF-1衛(wèi)星數(shù)據(jù)和實(shí)測水文數(shù)據(jù)，利用歸一化水體指數(shù)（NDWI）提取了丹江口水庫的水域面積，并根據(jù)水位、水域面積和蓄水量三者之間的相關(guān)性，分析了丹江口水庫的供水情況，是遙感數(shù)據(jù)與實(shí)測數(shù)據(jù)相結(jié)合的應(yīng)用。王慶健[5]利用Sentinel-1A提取的水域面積、“全國水雨情網(wǎng)”公布的實(shí)測水位和標(biāo)準(zhǔn)蓄水量反演了水庫參數(shù)，并通過蓄水量平衡公式計算水庫的滲透量等，最后，對影響蓄水量的各種因素進(jìn)行了分析。這些研究采集信息的時間長、工作量大、耗資多，所布設(shè)的儀器也常常會因?yàn)閻毫犹鞖夂腿粘Ｏ亩鴵p壞，從而導(dǎo)致數(shù)據(jù)的缺失;一段時期內(nèi)，獲取的可用遙感數(shù)據(jù)少，并且在相關(guān)模型建立之后仍需輸入實(shí)測數(shù)據(jù)才能獲得參數(shù)結(jié)果，對實(shí)測數(shù)據(jù)的依賴較大，精度相對較低。

相對于傳統(tǒng)方法和單一遙感方法，多源遙感方法能夠使各遙感技術(shù)優(yōu)勢互補(bǔ)，數(shù)據(jù)源豐富，水庫特征參數(shù)反演精度高。本文利用星載的雷達(dá)高度計（Rader Altimetry，RA）提取水庫的水位，以及星載的多光譜成像儀（Multi-Spectral Instrument，MSI）和合成孔徑雷達(dá)（Synthetic Aperture Radar，SAR）聯(lián)合提取水庫的水域面積。提取水域面積利用的是衛(wèi)星有效載荷生產(chǎn)的多源遙感數(shù)據(jù)，會存在一定的誤差。誤差來源主要有3個方面：① 不同衛(wèi)星所搭載的傳感器本身會有一定誤差;② 不同衛(wèi)星過境同一水庫的時間不同，會增加以小時、天等為時間單位的計算誤差;③ 不同衛(wèi)星傳感器的成像機(jī)理不同，所生產(chǎn)的數(shù)據(jù)需要采用特定的處理方法，這些處理方法精度不一，會造成一定的誤差[6-7]。雖然可以采用軌跡間和衛(wèi)星間等偏差的先驗(yàn)調(diào)整，但不同衛(wèi)星的修正方法使得水庫特征參數(shù)的計算變得繁瑣與復(fù)雜，因此，綜合考慮多源遙感技術(shù)在水庫特征參數(shù)上的研究并不多[8]。本文不去深究修正上述3種誤差的具體方法，而是在驗(yàn)證水域面積精度的基礎(chǔ)上，通過建立差異函數(shù)直接對提取的水域面積結(jié)果進(jìn)行修正，減少了中間過程中的復(fù)雜計算，節(jié)約了大量資源，為利用多源遙感技術(shù)開展水庫特征參數(shù)的反演研究提供了有益參考。

1研究區(qū)域與數(shù)據(jù)資料

1.1研究區(qū)域

丹江口水庫（32°36′～33°48′N、110°59′～111°49′E）位于漢江中上游，如圖1所示，是國家一級水源保護(hù)區(qū)和國家級生態(tài)文明示范區(qū)，作為我國南水北調(diào)中線工程的水源地，保障了河南、河北、北京、天津四?。ㄊ校┑纳a(chǎn)和生活用水[9]。

丹江口水庫從1958年9月開工，到1973年結(jié)束初期工程的建設(shè)，自建庫以來發(fā)揮了巨大的經(jīng)濟(jì)、社會和生態(tài)效益。2013年8月完成對二期工程的驗(yàn)收，設(shè)計壩頂高程為176.6 m，正常水位為170 m，水域面積可達(dá)1 050 km2，蓄水量可達(dá)290.5億m3，2014年12月正式向南水北調(diào)中線一期工程輸水[10-11]。

1.2數(shù)據(jù)資料

采用歐洲航天局的Sentinel系列衛(wèi)星作為遙感數(shù)據(jù)源，獲取地址為：https：∥scihub.copernicus.eu/dhus/，采用“全國水雨情網(wǎng)”上公布的實(shí)測數(shù)據(jù)作為地面驗(yàn)證的數(shù)據(jù)源，獲取地址為：http：∥xxfb.mwr.cn/，所有數(shù)據(jù)的獲取日期為2019年1月1日至2019年12月31日，其中，Sentinel-1A SAR影像29景、Sentinel-2A/B MSI影像14景和Sentinel-3A RA數(shù)據(jù)13份等。

重訪周期為12 d的Sentinel-1A衛(wèi)星和重訪周期為27 d的Sentinel-3A衛(wèi)星是發(fā)射電磁波的主動工作模式，所搭載的合成孔徑雷達(dá)，不受云霧和天氣的干擾，具有全天時、全天候的特點(diǎn)。重訪周期為5 d的Sentinel-2A/B衛(wèi)星是反射光波的被動工作模式，所搭載的多光譜成像儀，可提供媲美人眼的光學(xué)數(shù)據(jù)，相對于灰度值的成像模式，其對地物目標(biāo)的光譜差異性具有更加細(xì)微的捕捉能力[12-13]。

2研究方法

2.1水位反演

利用Sentinel-3A雷達(dá)高度計所生產(chǎn)的RA數(shù)據(jù)反演水庫水位，理想情況下，RA測量結(jié)果為衛(wèi)星到水面的瞬時距離[14]，如式（1）所示。

h=cτ2（1）

式中：h為RA到水面的高度;c為光速;τ為從RA發(fā)射信號到接收信號的往返時間。

RA在測量過程中會受許多因素的影響，例如大氣、電離層、波浪等，實(shí)際計算時必須消除這些影響，主要有干對流層Rangedt、濕對流層Rangewt、電離層Rangei、固體潮Rangeset、極潮Rangept以及水面橢球高GeoEGM96的修正，如式（2）所示。受誤差影響和參考水準(zhǔn)面的選取，該文通過Ice-1波形重跟蹤算法處理Sentinel-3A RA數(shù)據(jù)獲得以中國黃海水準(zhǔn)面為基準(zhǔn)的遙測水位，如式（3）所示，Range Cor為總誤差修正量，Height為水位，Altitude為衛(wèi)星到地球橢球面距離，Range為衛(wèi)星到水面距離。

2.2水域面積反演

2.2.1MSI影像反演水域面積

由MSI影像反演水域面積的處理流程如圖2所示，包括影像預(yù)處理、水體信息增強(qiáng)、水陸分割、后處理，以及水域面積計算等步驟。

其中，水體信息增強(qiáng)是為了加強(qiáng)水體與非水體之間的對比度，以便提高后續(xù)處理的精度;水陸分割用于水體提取，將影像特征相似的非水體區(qū)分開來;后處理主要用于對分割結(jié)果中可能存在的干擾信息進(jìn)行去除，例如橋梁、浮島、池塘等地物的影響，并獲得完整水體輪廓。

圖3展示了Sentinel-2A/B衛(wèi)星所生產(chǎn)的MSI影像反演水域面積過程中關(guān)鍵環(huán)節(jié)的處理結(jié)果。采用NDWI指數(shù)對圖3（a）所示的MSI影像圖以地物綠波段和近紅外波段反射率的差異性增強(qiáng)水體信息[15]，處理結(jié)果如圖3（b）所示。水體信息增強(qiáng)后的二值圖中，水體呈現(xiàn)較低灰度值的黑色，非水體呈現(xiàn)較高灰度值的灰色，具有良好的區(qū)分度。水陸分割利用決策樹將具有相似性的水體與地形陰影進(jìn)行區(qū)分，結(jié)果如圖3（c）所示。圖3（c）中可以明顯看出存在非水體的干擾信息，可以通過設(shè)置連通域閾值方式對小連通域干擾物進(jìn)行去除，從而獲得水體輪廓的提取結(jié)果，如圖3（d）所示。

基于提取到的水體，采用如式（4）所示的像元統(tǒng)計方法進(jìn)行計算，可獲得水域面積Smw。

式中：M為影像中像元總數(shù)目;N為提取到的水體像元數(shù)目;Sim為影像所對應(yīng)的地表區(qū)域面積，km2;Smw為最終反演的水域面積，km2。

2.2.2SAR影像反演水域面積

SAR影像反演水域面積的處理流程如圖4所示。首先要進(jìn)行輻射定標(biāo)、多視處理、斑點(diǎn)濾波、幾何校正等處理，以提高影像質(zhì)量，為獲取更加精準(zhǔn)水體信息奠定基礎(chǔ)[16-17]。

預(yù)處理后的Sentinel-1A SAR影像如圖5（a）所示。由細(xì)節(jié)放大圖可知，該影像具有較高分辨率。水陸分割是一個區(qū)分水體與非水體的二分類過程，模糊C均值算法（Fuzzy C-Means，F(xiàn)CM）是一種常用分類方法，該方法以隨機(jī)聚類中心開始，經(jīng)過不斷的迭代優(yōu)化完成地物類型歸屬劃分[18]。經(jīng)水陸分割和小連通域去除的水體輪廓提取結(jié)果如圖5（b）所示。

采用FCM算法進(jìn)行水陸分割的關(guān)鍵是確定影像中水體與非水體的臨界灰度值信息，從而獲得更加準(zhǔn)確的聚類中心。因此，該文對SAR影像像元的灰度值進(jìn)行了統(tǒng)計，其直方圖如圖5（c）所示。

可以明顯地看到直方圖中出現(xiàn)了2個波峰，第一個波峰主要為SAR影像中的水體，其像元灰度值分布在0～8內(nèi)，第二個波峰以及剩余灰度值主要代表SAR影像中的非水體?；叶戎?為水體與非水體的臨界值，為了避免漏提淺水水體，結(jié)合灰度統(tǒng)計直方圖和實(shí)踐研究，應(yīng)將FCM算法的一個關(guān)鍵聚類中心設(shè)置為8。

2.3多源遙感水域面積修正研究

MSI影像的觀測周期較短，影像信息豐富，但受天氣和光照的影響較大，故MSI影像實(shí)際可用的數(shù)據(jù)較少。SAR影像可全天時、全天候工作，但其重返周期較長，加上受限于軍事保密、技術(shù)瓶頸和成本控制，獲取到合適的SAR影像數(shù)據(jù)也相對較少。針對上述情況，將MSI影像和SAR影像聯(lián)合應(yīng)用于研究，可以豐富數(shù)據(jù)源、提高實(shí)驗(yàn)精度。

MSI影像和SAR影像屬于異源傳感器影像，前者是被動成像模式，利用光波的反射成像，后者是主動成像模式，利用電磁波的散射成像，兩者具有不同的成像機(jī)理，出現(xiàn)了“同物異譜”現(xiàn)象。為了獲得后續(xù)更精確的蓄水量反演結(jié)果，必須構(gòu)建兩者水域面積差異的定量函數(shù)關(guān)系，修正上述誤差。

在構(gòu)建定量函數(shù)關(guān)系之前，需要先驗(yàn)證水域面積的提取精度。水體輪廓對水域面積的提取精度有著至關(guān)重要的影響，只要驗(yàn)證水體輪廓的提取精度即可說明水域面積的提取精度。利用遙感數(shù)據(jù)提取水體輪廓，手動提取法精度最高，但其效率低下。以手動提取的水體輪廓為驗(yàn)證數(shù)據(jù)，通過漏警率、虛警率、準(zhǔn)確率和Kappa系數(shù)等指標(biāo)，對提取的水體輪廓進(jìn)行精度驗(yàn)證，其定量分析結(jié)果如表1所列。可知，2種影像水體輪廓的漏警率和虛警率都較低，準(zhǔn)確率優(yōu)于0.97，Kappa系數(shù)高于0.91，表明提取的水體輪廓足夠用于后續(xù)研究。

對SAR水域面積和MSI水域面積進(jìn)行指數(shù)、冪次、多項(xiàng)式等不同函數(shù)的擬合，并選取最優(yōu)函數(shù)來建立多源遙感水域面積差異的定量函數(shù)關(guān)系，所得定量分析的結(jié)果如表2所列。

根據(jù)R2越大越好、RMSE越小越好和多項(xiàng)式次數(shù)越低計算越簡單的定量分析原則，四次多項(xiàng)式函數(shù)具有最小的RMSE、較大的R2和較低的多項(xiàng)式次數(shù)。因此，采用四次多項(xiàng)式函數(shù)表示MSI影像和SAR影像所提取水域面積之間的差異。

四次多項(xiàng)式函數(shù)的表達(dá)式如式（5）所示。

F（x）=ax4+bx3+cx2+dx+e （5）

式中：x為水位，m;F（x）為同一水位下MSI影像和SAR影像提取水域面積的差異值，m2;a、b、c、d、e為多項(xiàng)式系數(shù)，擬合后的值分別為-0.000 476，0.302 6，-72.1，7 632，-0.000 030 28。

丹江口水庫位于北半球和南半球冷暖氣流的交匯處，易受副熱帶高壓北進(jìn)擺動和穩(wěn)定南退的影響，天氣變化劇烈，水庫上空云霧天氣的概率較大。對于選取的部分MSI影像，不可避免地存在少量的云霧，其遮擋效應(yīng)在一定程度上影響了水域面積提取的精度。相比之下，SAR影像具有全天時、全天候的特點(diǎn)，不存在上述誤差。另外，從水體輪廓提取結(jié)果來說，SAR影像所提取水體輪廓精度要優(yōu)于MSI影像。因此，考慮到上述情況，該文以SAR影像為基準(zhǔn)，修正水域面積差異，即將MSI影像提取的水域面積值P（x），加上面積誤差修正值F（x），得到修正后的MSI水域面積S（x），如式（6）所示。修正后的MSI水域面積與SAR水域面積聯(lián)合在一起稱為多源遙感水域面積。

S（x）=P（x）+F（x）（6）

2.4水庫蓄水量反演

水庫蓄水量是無法直接觀測的量，一般通過水位和水域面積進(jìn)行間接反演。該文把水庫理想化為棱臺，獲得鄰水位或者相鄰橫斷面之間的一小部分蓄水量，并把這些小蓄水量累加起來，再加上死水位對應(yīng)的蓄水量就能求得整個水庫的總蓄水量[19]，如式（7）～（9）所示。

3結(jié)果與分析

3.1水位結(jié)果及分析

根據(jù)Sentinel-3A RA數(shù)據(jù)提取遙測水位的獲取日期查找到對應(yīng)的實(shí)測水位，然后將遙測水位和實(shí)測水位進(jìn)行對比，結(jié)果如圖6所示，日期范圍是2019年1月1日至12月31日。

由圖6可知，2019年丹江口水庫的水位變化范圍為150～167m。隨著旱季和汛季的轉(zhuǎn)換，水位呈現(xiàn)規(guī)律性波動，在7～10月的汛期，水位上漲迅速，11月之后將逐漸進(jìn)入枯水期，水位開始慢慢下降，并且由Sentinel-3A RA提取的遙測水位與實(shí)測水位的波動趨勢相吻合。

采用R2和RMSE對遙測水位的精度進(jìn)行評價，得到R2=0.99，RMSE=0.0933，表明提取的遙測水位具有很高的精度，遙測水位與實(shí)測水位的平均誤差僅有9.33 cm。

3.2水域面積結(jié)果及分析

利用星載MSI影像和星載SAR影像分別提取的MSI水域面積和SAR水域面積，結(jié)果如圖7～8所示。

由圖7～8可知：2019年丹江口水庫的水域面積在5月初達(dá)到最小值，7月開始明顯增大，在11月中旬達(dá)到最大值，然后再次回落，最大水域面積與最小水域面積相差約250 km2。隨著時間推移，整體上MSI水域面積和SAR水域面積變化趨勢較為吻合，MSI水域面積和SAR水域面積的準(zhǔn)確率分別達(dá)到了97.95%和98.95%。

3.3建立“水位-水域面積”模型

由2.4小節(jié)可知：只要獲得水位和水域面積就能在不依靠地面實(shí)測數(shù)據(jù)的情況下獲得蓄水量。對于水位和水域面積，已知其中一個參數(shù)，想要快速獲得另外一個參數(shù)，需要建立“水位-水域面積”關(guān)系曲線。該關(guān)系曲線建立之后，再通過多源遙感技術(shù)獲得水位或水域面積任何一個參數(shù)，就能獲得包括蓄水量在內(nèi)的另外兩個參數(shù)。

以多源遙感技術(shù)獲取的遙測水位、MSI水域面積、SAR水域面積以及修正后的多源遙感水域面積分別建立“遙測水位-MSI水域面積”關(guān)系曲線、“遙測水位-SAR水域面積”

關(guān)系曲線和“遙測水位-多源遙感水域面積”關(guān)系曲線。以“遙測水位-多源遙感水域面積”關(guān)系曲線為例，建立過程如下：

首先構(gòu)建“遙測水位-多源遙感水域面積”關(guān)系曲線的指數(shù)函數(shù)、冪函數(shù)和多項(xiàng)式函數(shù)的數(shù)學(xué)模型，由于四次及四次以上多項(xiàng)式函數(shù)出現(xiàn)過擬合現(xiàn)象，選擇指數(shù)函數(shù)、冪函數(shù)和四次以下多項(xiàng)式函數(shù)對擬合效果作進(jìn)一步驗(yàn)證，如表3所列。

由表3可知：指數(shù)函數(shù)的R2最小，RMSE最大;冪函數(shù)與指數(shù)函數(shù)的R2和RMSE最接近;三次多項(xiàng)式函數(shù)的R2最大，但RMSE不是最小，其R2只比二次多項(xiàng)式函數(shù)的R2多0.0003，二次多項(xiàng)式函數(shù)的RMSE最小。根據(jù)R2越大越好，RMSE越小越好的定量分析原則，二次多項(xiàng)式函數(shù)最適合描述遙測水位與多源遙感水域面積之間的關(guān)系。

得到“遙測水位-多源遙感水域面積”模型如式（10）所示，關(guān)系曲線如圖9（a）所示;同理可得“遙測水位-MSI水域面積”模型如式（11）所示，關(guān)系曲線如圖9（b）所示;“遙測水位-SAR水域面積”模型如式（12）所示，關(guān)系曲線如圖9（c）所示。

“水位-水域面積”關(guān)系曲線建立之后，各模型最優(yōu)函數(shù)的R2和RMSE如表4所列。

由表4可知：各模型本身也已經(jīng)具有較高的R2和較低的RMSE，相比單獨(dú)使用MSI影像或SAR影像，聯(lián)合使用兩種遙感影像后，“遙測水位-多源遙感水域面積”模型的R2分別提高了0.26%和0.60%，RMSE分別降低了2.63和1.29，即R2和RMSE均有所優(yōu)化。

3.4蓄水量結(jié)果及分析

利用多源遙感技術(shù)反演得到水庫的遙測水位和多源遙感水域面積，進(jìn)而得到水庫的蓄水量，將反演的蓄水量與“全國水雨情網(wǎng)”上公布的標(biāo)準(zhǔn)蓄水量進(jìn)行對比，對比結(jié)果如圖10所示，日期范圍是2019年1月1日至12月31日。

由圖10可知，反演蓄水量和標(biāo)準(zhǔn)蓄水量的波動趨勢具有較高吻合度。在7～10月的汛期，相關(guān)地區(qū)降雨較多，丹江口水庫的蓄水量增長迅速，在11月17日左右達(dá)到最高峰，然后開始慢慢減少，逐漸進(jìn)入枯水期，在5月初達(dá)到最低值。整體上，隨著旱季和汛季的轉(zhuǎn)換，其蓄水量也呈現(xiàn)規(guī)律性的變化。

采用R2和RMSE進(jìn)行定量分析。分析結(jié)果表明，R2=0.99，說明反演蓄水量與標(biāo)準(zhǔn)蓄水量具有極高的吻合度，能夠真實(shí)地反映水庫蓄水量的變化情況;另外，RMSE=1.177，表明反演的水庫蓄水量標(biāo)準(zhǔn)差在1.177億m3以內(nèi)，相對誤差平均為1.8%。

4結(jié) 語

本文以丹江口水庫為例，利用多源遙感技術(shù)反演了水庫參數(shù)，其中，利用RA數(shù)據(jù)提取了水庫的水位，利用多源傳感器生產(chǎn)的MSI影像和SAR影像提取水庫的水域面積。通過構(gòu)建多源傳感器水域面積之間的定量函數(shù)關(guān)系，修正了二者水域面積之間的誤差，最后，基于修正后多源遙感水域面積，結(jié)合蓄水量模型反演水庫的蓄水量。

實(shí)驗(yàn)結(jié)果表明，利用多源遙感技術(shù)反演的水庫特征參數(shù)能夠真實(shí)地反映水庫健康運(yùn)行的狀態(tài)。在水庫特征參數(shù)反演過程中，除了利用地面實(shí)測數(shù)據(jù)對反演結(jié)果進(jìn)行驗(yàn)證之外，不再依賴實(shí)測數(shù)據(jù)的支持就能達(dá)到對水庫特征參數(shù)的反演，降低了成本，節(jié)約了資源，有效補(bǔ)充了實(shí)測數(shù)據(jù)缺失情況下的研究工作。同時，利用多源遙感技術(shù)反演水庫特征參數(shù)能夠達(dá)到對水庫異常預(yù)警和日常管理的最大效益化，特別是對堰塞湖或出現(xiàn)大面積違法的圍湖造田等的預(yù)警。該文關(guān)于水庫特征參數(shù)的反演方法不僅能夠克服當(dāng)前傳統(tǒng)方法和單一遙感方法的不足，也可以進(jìn)一步推廣到其它大型水庫特征參數(shù)的反演工作當(dāng)中。未來還可以考慮采用更多類型的遙感數(shù)據(jù)和更加優(yōu)良的水庫蓄水量模型等參與反演。

參考文獻(xiàn)：

[1]周建中，頓曉晗，張勇傳.基于庫容風(fēng)險頻率曲線的水庫群聯(lián)合防洪調(diào)度研究[J].水利學(xué)報，2019，50（11）：1318-1325.

[2]李書杰.水庫庫容的量算精度[J].東北水利水電，1997（7）：34-36.

[3]蘇業(yè)助，洪為善.根據(jù)水量平衡原理修正水庫庫容曲線方法[J].人民長江，1999，30（12）：39-41.

[4]孫建蕓，袁琳，王新生，等.基于GF_1衛(wèi)星的丹江口水庫水面面積-蓄水量-水位相關(guān)性研究[J].南水北調(diào)與水利科技，2017，15（5）：89-96.

[5]王慶健.基于SAR監(jiān)測的水庫主要參數(shù)計算系統(tǒng)的設(shè)計與實(shí)現(xiàn)[D].開封：河南大學(xué)，2018.

[6]BONNEFOND P，EXERTIER P，LAURAIN O，et al.Absolute calibration of Jason-1 and Jason-2 altimeters in Corsica during the for- mation flight phase[J].Marine Geodesy，2010，33（S1）：80-90.

[7]MERTIKAS S P，IOANNIDES R T，TZIAVOS I N，et al.Statistical models and latest results in the determination of the absolute bias for the radar altimeters of Jason satellites using the gavdos facility[J].Marine Geodesy，2010，33（S1） ：114-149.

[8]王文種，黃對，劉九夫，等.基于Landsat與Sentinel-3A衛(wèi)星數(shù)據(jù)的當(dāng)惹雍錯1988～2018年湖泊水位-水量變化及歸因[J].湖泊科學(xué)，2020，32（5）：1552-1563.

[9]王元超，王旭，雷曉輝，等.丹江口水庫入庫徑流特征及其演變規(guī)律[J].南水北調(diào)與水利科技，2015，13（1）：15-19.

[10]吳川，張玉龍，許秀貞，等.基于Landset TM/ETM和HJ-1A/B影像的丹江口水庫水域變化監(jiān)測研究[J].長江流域資源與環(huán)境，2013，22（9）：108-114.

[11]尹志杰，王容，趙蘭蘭，等.丹江口水庫以上流域夏秋汛期來水及相關(guān)性分析[J].水文，2019，39（4）：74-77.

[12]楊魁，楊建兵，江冰茹.Sentinel-1衛(wèi)星綜述[J].城市勘測，2015（2）：26-29.

[13]田穎，陳卓奇，惠鳳鳴，等.歐空局哨兵衛(wèi)星Sentinel-2A/B數(shù)據(jù)特征及應(yīng)用前景分析[J].北京師范大學(xué)學(xué)報（自然科學(xué)版），2019，55（1）：57-65.

[14]郭華東.雷達(dá)對地觀測理論與應(yīng)用[M].北京：科學(xué)出版社，2000.

[15]GAO B C.NDWI—A normalized difference water index for remote sensing of vegetation liquid water from space[J].Remote Sensing of Environment，1995，58（3）：257-266.

[16]潘習(xí)哲.星載SAR圖像處理[M].北京：科學(xué)出版社，1996.

[17]王超.SAR圖像處理新進(jìn)展[J].中國圖象圖形學(xué)報，2009，14（1）：1-2.

[18]BEZDEK J C，EHRLICH R，F(xiàn)ULL W.FCM：The fuzzy c-means clustering algorithm[J].Computers & Geoences，1984，10（2-3）：191-203.

[19]張國慶.青藏高原湖泊變化遙感監(jiān)測及其對氣候變化的響應(yīng)研究進(jìn)展[J].地理科學(xué)進(jìn)展，2018，（2）：214-223.

（編輯：黃文晉）

Abstract：Take Danjiangkou Reservoir as an example，the water area was extracted by using the Sentinel-1A Synthetic Aperture Radar （SAR） and Sentinel-2A/B Multi-Spectral Imager （MSI），the difference of water areas by SAR and by MSI was corrected by establishing a quantitative function relation of these areas，and the revised water area was combined with the inverted water level of Sentinel-3A Radar Altimeter （RA） to calculate the reservoir volume.Finally，the retrieval precision of the reservoir parameters was evaluated according to the measured data published by the National Water and rain regime network.The results showed that R2 was above 0.99 and RMSE was 0.0933，which meant that the average error between inverted water level and measured water level was only 9.33 cm.Compared with the water area extracted from single MSI image or SAR image，the R2 and RMSE of the "remote sensing water level-Multi-Source Remote Sensing water area" model were optimized，especially R2 increased by 0.26% and 0.60% ，while RMSE decreased by 2.63 and 1.29，respectively.Compared with the standard water storage，for the novel model，R2 reached 0.99 and RMSE reached 1.177，respectively，and the relative error was 1.8% on average.The result shows that the reservoir characteristic parameters inversion using multi-source remote sensing technology can more accurately reflect the operation of the reservoir.

Key words：reservoir characteristic parameters;multi-source remote sensing;Rader Altimetry;Multi-Spectral Instrument;Synthetic Aperture Radar