沈五偉，陳一帆，申振東

（1.浙江省水利河口研究院，浙江 杭州 310020；2.浙江廣川工程咨詢(xún)有限公司，浙江 杭州 310020）

河道糙率反問(wèn)題研究回顧與展望

沈五偉1，2，陳一帆1，2，申振東1，2

（1.浙江省水利河口研究院，浙江 杭州 310020；2.浙江廣川工程咨詢(xún)有限公司，浙江 杭州 310020）

河道糙率是河流水力計(jì)算的關(guān)鍵靈敏參數(shù)，也是眾多因素影響的綜合性參數(shù)。工程上常用的查表法、公式法和試錯(cuò)法得出的糙率值往往難以滿(mǎn)足水力計(jì)算要求，因此河道糙率反演理論和方法研究引起了國(guó)內(nèi)外學(xué)者的廣泛重視。在回顧和總結(jié)河道糙率確定基本方法和優(yōu)化反演方法研究進(jìn)展的基礎(chǔ)上，指出了河道糙率反問(wèn)題研究面臨的困難和挑戰(zhàn)，展望了河道糙率反問(wèn)題研究的發(fā)展方向。

河道糙率；反演分析；糙率率定

1 問(wèn)題的提出

河道糙率，本質(zhì)上是反映床面岸灘粗糙程度、河道空間不規(guī)則性、水流水力狀態(tài)變化等眾多因素對(duì)水流運(yùn)動(dòng)阻力影響的重要水力參數(shù)，是河流水力計(jì)算的關(guān)鍵靈敏參數(shù)，其取值的合理性直接關(guān)系到河流水力計(jì)算結(jié)果的精確性和可靠性。雖然河道糙率時(shí)空變化已有一些經(jīng)驗(yàn)規(guī)律可循，但由于影響河道糙率的因素眾多、復(fù)雜多變，因此對(duì)其進(jìn)行量測(cè)計(jì)算相當(dāng)困難，常采用間接方法進(jìn)行確定，包括河道糙率確定基本方法和優(yōu)化反演方法2大類(lèi)。

針對(duì)河道糙率，本文在回顧和總結(jié)河道糙率計(jì)算研究成果的基礎(chǔ)上，提出了河道糙率優(yōu)化反演方法的評(píng)價(jià)標(biāo)準(zhǔn)，指出了河道糙率反問(wèn)題研究面臨的困難和挑戰(zhàn)，最后展望了河道糙率反問(wèn)題研究的發(fā)展方向。

2 河道糙率確定基本方法

由于天然河流的復(fù)雜性，僅依據(jù)糙率影響因素對(duì)糙率進(jìn)行定量分析尚不可行。通過(guò)長(zhǎng)期理論和實(shí)踐研究，河道糙率確定已有一些經(jīng)驗(yàn)公式、圖表對(duì)照或調(diào)試演算等基本方法供工程技術(shù)人員參考使用。目前，糙率確定的基本方法見(jiàn)表1。

表1 河道糙率確定基本方法表

上述糙率確定基本方法具有較大的盲目性和不確定性，具體表現(xiàn)為：①糙率確定受人為因素影響明顯，率定結(jié)果往往因人而異，具有一定的隨意性；②隨著河段數(shù)量的增多，糙率確定工作量和率定難度顯著增大，對(duì)工程技術(shù)人員的經(jīng)驗(yàn)要求隨之提高；③由于缺乏系統(tǒng)的理論指導(dǎo)，率定結(jié)果的可靠性難以有效評(píng)估，結(jié)果的準(zhǔn)確性得不到合理保障。因此河道糙率反演研究（或稱(chēng)為河道糙率辨識(shí)研究）引起了國(guó)內(nèi)外學(xué)者的廣泛重視。

3 河道糙率確定優(yōu)化反演方法

3.1 河道糙率優(yōu)化反演含義

河道糙率優(yōu)化反演是根據(jù)水文測(cè)站測(cè)得的水位、流量等水力要素觀測(cè)資料，采用系統(tǒng)優(yōu)化方法，將河道糙率的率定過(guò)程編制成優(yōu)化反演程序來(lái)實(shí)現(xiàn)糙率的自動(dòng)率定，這樣既保證了水力要素觀測(cè)值與模擬值的合理匹配，又可顯著降低人工試算的工作量。由于描述一維非恒定明渠水流運(yùn)動(dòng)的圣維南方程組是非線性雙曲型偏微分方程組，斷面糙率和水力要素之間呈現(xiàn)明顯的非線性關(guān)系，同時(shí)由于河道水流是復(fù)雜多變的動(dòng)態(tài)系統(tǒng)，因此河道糙率優(yōu)化反演問(wèn)題為動(dòng)態(tài)系統(tǒng)非線性反演問(wèn)題。

3.2 河道糙率優(yōu)化反演研究現(xiàn)狀

Becker等[1]采用結(jié)合敏感度分析的單純形法進(jìn)行了河道糙率的反演研究；WasanthaLal等[2]建立了基于奇異值分解的河道糙率反演模型，并提出了糙率分辨率的概念；Atanov等[3]采用拉格拉日算子變分法對(duì)規(guī)則棱柱形河道的糙率進(jìn)行了反演研究；Khatibi等[4]采用高斯 — 牛頓優(yōu)化方法進(jìn)行河道糙率的反演研究，分析了準(zhǔn)則函數(shù)和樣本數(shù)量對(duì)反演結(jié)果的影響；Ramesh R.等[5]建立了基于順序二次規(guī)劃算法的河道糙率反演模型，分析了不同噪聲水平下反演模型的率定效果；董文軍等[6]利用Frechet微分概念和構(gòu)造相應(yīng)的協(xié)態(tài)方程來(lái)確定目標(biāo)函數(shù)的下降方向，并用牛頓一辛普森迭代法求解糙率優(yōu)化反演模型；李光熾等[7]較先探索采用卡爾曼濾波進(jìn)行河道糙率反問(wèn)題的研究；程偉平等[8]基于廣義逆和Backus - Gilbert反演理論，構(gòu)造了用于河道糙率反演的牛頓 — 廣義逆方法和自然逆方法，并推導(dǎo)了2種算法反演解估計(jì)的分辨率矩陣和單位協(xié)方差矩陣；Ding，Y.等[9]建立了基于誤差最小化分析的有限記憶準(zhǔn)牛頓法，并在淺水河道中進(jìn)行糙率的反演分析；霍光等[10]首先采用多相模糊統(tǒng)計(jì)法確定河道糙率的一個(gè)小范圍鄰域，再采用決策方案優(yōu)選法以進(jìn)一步確定糙率優(yōu)選值；Bilgil，A.等[11]將人工神經(jīng)網(wǎng)絡(luò)法應(yīng)用于河道糙率的優(yōu)化率定，得到了較好試驗(yàn)結(jié)果；吳曉玲等[12]將卡爾曼濾波技術(shù)用于河道糙率的實(shí)時(shí)修正以提高了洪水預(yù)報(bào)精度；辛小康等[13]嘗試建立基于改進(jìn)遺傳算法的河網(wǎng)糙率反演模型；李麗等[14]測(cè)試了自適應(yīng)隨機(jī)搜索(ARS)算法在資料受限條件下河網(wǎng)模型糙率反演的表現(xiàn)，得到了較理想的反演效果；陳一帆等[15]從河道糙率空間分布特性和專(zhuān)家經(jīng)驗(yàn)知識(shí)2方面出發(fā)，分別提出了糙率空間分布緩變性和糙率經(jīng)驗(yàn)估值2種先驗(yàn)知識(shí)，用于天然河流的糙率反演；包紅軍等[16]以糙率修正量為狀態(tài)變量，觀測(cè)值與模擬值之差為量測(cè)變量，建立了基于卡爾曼濾波的河道糙率反演模型；陳素紅等[17]建立了基于多親遺傳算法的河道糙率率定模型，將其應(yīng)用南水北調(diào)中線工程總干渠的部分渠道斷面糙率的率定；陳一帆等[18]以糙率和水力狀態(tài)量作為河網(wǎng)非線性動(dòng)態(tài)系統(tǒng)狀態(tài)變量，采用擴(kuò)展卡爾曼濾波構(gòu)建了結(jié)合糙率動(dòng)態(tài)校正的河網(wǎng)水情數(shù)據(jù)同化模型，分析了水位動(dòng)態(tài)噪聲水平、糙率動(dòng)態(tài)噪聲水平、糙率初始值及測(cè)站個(gè)數(shù)對(duì)模型校正的影響。

4 河道糙率優(yōu)化反演方法評(píng)價(jià)標(biāo)準(zhǔn)

在現(xiàn)有河道糙率反演方法中，多數(shù)反演模型以最小化觀測(cè)值與模擬值之差作為模型的目標(biāo)函數(shù)，研究如何能夠搜索到目標(biāo)函數(shù)的最優(yōu)值，但對(duì)模型的數(shù)值穩(wěn)定性、結(jié)果的合理性等反演特性考慮并不多。評(píng)價(jià)河道糙率反演方法性能的主要標(biāo)準(zhǔn)包括以下3方面。

4.1 初始值影響評(píng)價(jià)

對(duì)于欠定或混定的反演模型，反演結(jié)果受初始條件的影響較大，經(jīng)常會(huì)因初值給得不恰當(dāng)而出現(xiàn)結(jié)果發(fā)散的現(xiàn)象或者得到不同的反演結(jié)果。一個(gè)穩(wěn)定的數(shù)學(xué)模型或數(shù)值求解模型應(yīng)該能在不同初始條件下給出相近甚至相同的結(jié)果，它包括2個(gè)方面的意義[15]：①?gòu)臄?shù)學(xué)模型的角度來(lái)講，不會(huì)因?yàn)槌踔颠x取的原因而導(dǎo)致不同解，即滿(mǎn)足解的唯一性；②從數(shù)值求解的角度來(lái)講，不會(huì)因?yàn)槌踔颠x取的原因而導(dǎo)致計(jì)算結(jié)果失效或發(fā)散。

4.2 監(jiān)測(cè)信息量影響評(píng)價(jià)

實(shí)際工程中，能夠用于糙率率定的水文監(jiān)測(cè)信息往往較少，導(dǎo)致反演模型常為欠定或者混定。對(duì)于一個(gè)穩(wěn)健的反演模型而言，它應(yīng)具有較好的魯棒性，即在可用監(jiān)測(cè)信息量較少時(shí)亦能給出較為合理的反演結(jié)果，并隨信息量的增加，反演結(jié)果逐漸逼近于真實(shí)解[15]。監(jiān)測(cè)信息的增加有2種方式：①單個(gè)觀測(cè)站點(diǎn)時(shí)間維信息量的角度，增加其歷時(shí)觀測(cè)數(shù)據(jù)；②從河網(wǎng)空間維的角度來(lái)講，增加觀測(cè)站的布置個(gè)數(shù)。

4.3 信息誤差影響評(píng)價(jià)

在真實(shí)世界中，觀測(cè)資料總是存有誤差。而誤差的存在可能會(huì)使反演結(jié)果不合理，難以反映河道真實(shí)糙率。對(duì)于一個(gè)穩(wěn)健的反演模型而言，它應(yīng)具有較好的誤差抗噪性，即在一定監(jiān)測(cè)誤差范圍內(nèi)能夠給出較為合理的反演結(jié)果，不因微小的擾動(dòng)導(dǎo)致巨大的數(shù)值振蕩[15]。

5 河道糙率優(yōu)化反演面臨的困難與挑戰(zhàn)

在回顧國(guó)內(nèi)外研究現(xiàn)狀的基礎(chǔ)上，結(jié)合客觀現(xiàn)實(shí)，認(rèn)為當(dāng)前河道糙率優(yōu)化反演面臨的困難和挑戰(zhàn)主要包括以下幾點(diǎn)：

（1）資料不完備。資料是反演分析的基礎(chǔ)，是制約反演分析的一大瓶頸，一般很難全面詳盡地獲取河流相應(yīng)的充足信息，且資料精度也難以得到保證，如斷面測(cè)量精度粗略、水力要素歷史觀測(cè)資料短缺、水工建筑物調(diào)度情況不明等等。

（2）水力要素觀測(cè)資料存在干擾和誤差。這可能使得優(yōu)化反演不易穩(wěn)定，即觀測(cè)數(shù)據(jù)中微小誤差可能導(dǎo)致反演過(guò)程的較大振蕩、反演結(jié)果的較大變動(dòng)，因此在進(jìn)行糙率優(yōu)化反演時(shí)，盡可能引入糙率先驗(yàn)知識(shí)（包括糙率時(shí)空分布特性、糙率經(jīng)驗(yàn)估值等信息），以減少反演過(guò)程的振蕩性，提高結(jié)果的合理性及可靠性。

（3）糙率反演問(wèn)題是多類(lèi)型參數(shù)共同率定的反問(wèn)題。在現(xiàn)有河網(wǎng)水力參數(shù)反演模型中，多數(shù)反演模型僅考慮河道糙率的率定，很少同步考慮水利工程過(guò)流特征參數(shù)的率定。然而，河網(wǎng)內(nèi)往往存在一定數(shù)量的水利調(diào)控工程，其是河網(wǎng)防洪除澇系統(tǒng)的重要組成部分，對(duì)水流運(yùn)動(dòng)影響較大。若忽視水利調(diào)控工程過(guò)流特征參數(shù)的率定，將難以保障河道糙率反演結(jié)果的合理性。

（4）糙率反演問(wèn)題是空間相關(guān)的反問(wèn)題。河道糙率沿著河流流向不斷變化，將其在空間方向上離散后，需要辨識(shí)的參數(shù)較多。然而，由于種種因素的限制，布置在河流上的水文測(cè)站通常比較稀疏，水文測(cè)站的個(gè)數(shù)往往少于河道糙率的個(gè)數(shù)，因而所建立的河道糙率反演模型通常是欠定的，加之缺乏糙率先驗(yàn)知識(shí)的應(yīng)用，極易因觀測(cè)誤差或數(shù)值擾動(dòng)導(dǎo)致糙率反演失敗。

（5）糙率反演問(wèn)題是時(shí)間相關(guān)的反問(wèn)題。河流水情狀態(tài)隨著時(shí)間不斷變化，而河道糙率及水利調(diào)控工程過(guò)流特征參數(shù)往往隨著水位的變化而變化，表現(xiàn)為一個(gè)過(guò)程，若按過(guò)程辨識(shí)的思維方式進(jìn)行水力參數(shù)的反分析，將不可避免地出現(xiàn)維數(shù)災(zāi)難。在現(xiàn)有研究成果中，多數(shù)反演模型將河道糙率或水利調(diào)控工程過(guò)流特征參數(shù)視為定值，進(jìn)行參數(shù)靜態(tài)反演，較少考慮參數(shù)的時(shí)間變化性，即參數(shù)的動(dòng)態(tài)反演問(wèn)題。

6 結(jié) 語(yǔ)

（1）由于水文測(cè)站數(shù)量有限或空間布置不合理，導(dǎo)致水情監(jiān)測(cè)信息量難以滿(mǎn)足河道糙率反演的需求，加之監(jiān)測(cè)信息往往存在干擾和誤差，這些都會(huì)使優(yōu)化反演不易穩(wěn)定、反演結(jié)果缺乏合理性，且監(jiān)測(cè)數(shù)據(jù)中微小誤差可能導(dǎo)致反演過(guò)程的較大振蕩。因此，探索研究河道糙率先驗(yàn)知識(shí)，建立具有一定物理意義的糙率修正機(jī)制，用于解決水情監(jiān)測(cè)數(shù)據(jù)不足條件下反演結(jié)果的不確定性，以期提高河道糙率反演過(guò)程的魯棒性和反演結(jié)果的合理性，將是河道糙率反問(wèn)題研究的發(fā)展方向之一。

（2）河道糙率往往隨時(shí)間不斷變化，如河道沖淤等下墊面變化引起的糙率變化、糙率隨水位漲落的變化等，然而當(dāng)前多數(shù)河道糙率反演模型將河道糙率視為定值，進(jìn)行糙率靜態(tài)反演，較少考慮糙率的時(shí)間變化性。因此，研究和解決河道糙率動(dòng)態(tài)反演問(wèn)題，體現(xiàn)糙率的時(shí)間可變性將是河道糙率反問(wèn)題研究的發(fā)展方向之一。

（3）河道糙率的率定，不僅要考慮糙率本身的率定，還要考慮水利工程過(guò)流特征參數(shù)的率定。然而，河道糙率和水利工程過(guò)流特征參數(shù)在進(jìn)行共同反演時(shí)，由于量綱、數(shù)量級(jí)、量變范圍的不同，極易引起反演過(guò)程出現(xiàn)劇烈振蕩，從而背離實(shí)際情況，使后續(xù)計(jì)算失效。因此，探索研究河道糙率和水利工程過(guò)流特征參數(shù)之間的反演機(jī)理和作用效應(yīng)，提出合理可行的共同反演方法，將是河道糙率反問(wèn)題研究的發(fā)展方向之一。

[1] Becker L，Yeh W W G．Identification of multiple reach channel parameters[J]．Water Resources Research，1973，9(2)：326 - 335．

[2] Wasantha Lal A M．Calibration of riverbed roughness[J]．Journal of Hydraulic Engineering，1995，121(9)：664 - 671．

[3] ATANOV G A，EVSEEVA E G，MESELHE E A．Estimation of roughness profile in trapezoidal open channel[J]．Journal of Hydraulic Engineering，ASCE，1999，125(3)：309 - 312．[4] Khatibi R H，Williams J J，Wormleaton P R．Identification problem of open - channel friction parameters[J]．Journal of Hydraulic Engineering，1997，123(12)：1078 - 1088．

[5] Ramesh R，Datta B，Bhallamudi S M，et al．Optimal estimation of roughness in open - channel flows[J]．Journal of Hydraulic Engineering， 2000，126(4)：299 - 303．

[6] 董文軍，姜亨余．一維水流方程中曼寧糙率的參數(shù)辨識(shí)[J]．天津大學(xué)學(xué)報(bào)(自然科學(xué)與工程技術(shù)版)，2001，34(2)：201 - 204．

[7] 李光熾，周晶晏，張貴壽．用卡爾曼濾波求解河道糙率參數(shù)反問(wèn)題[J]．河海大學(xué)學(xué)報(bào)(自然科學(xué)版)，2003，31(5)：490 - 493．

[8] 程偉平，劉國(guó)華．基于廣義逆理論的河網(wǎng)糙率反演研究[J]．浙江大學(xué)學(xué)報(bào)，2005，39(10)：1063 - 1068．

[9] Ding Y，Wang S S．Identification of Manning’s roughness coefficients in channel network using adjoint analysis[J]．International Journal of Computational Fluid Dynamics，2005，19(1)：3 - 13．

[10] 霍光，王義剛．基于多因素模糊綜合評(píng)判的河網(wǎng)糙率求解 [J]．河海大學(xué)學(xué)報(bào)(自然科學(xué)版)，2006，34(5)：518 - 521．

[11] Bilgil A，Altun H．Investigation of fl ow resistance in smooth open channels using artif i cial neural networks[J]．Flow Measurement and Instrumentation，2008，19(6)：404 - 408．

[12] 吳曉玲，王船海．基于水動(dòng)力學(xué)模型的實(shí)時(shí)糙率反推在洪水預(yù)報(bào)中的應(yīng)用[J]．水電能源科學(xué)， 2008，26(5)：43 - 45．

[13] 辛小康，劉剛，張向東，等．基于遺傳優(yōu)化的河網(wǎng)數(shù)學(xué)模型糙率參數(shù)反演[J]．水利水電科技進(jìn)展， 2009， 29(6)：21 - 24．

[14] 李麗，王加虎，王建群，等．自適應(yīng)隨機(jī)搜索算法在河網(wǎng)數(shù)學(xué)模型糙率反演中的應(yīng)用[J]．水利水電科技進(jìn)展，2011，31(5)：64 - 67．

[15] 陳一帆，程偉平，蔣建群．一種穩(wěn)健的河流糙率反演方法[J]．浙江大學(xué)學(xué)報(bào)(工學(xué)版)，2013，47(8)：1361 - 1365，1378．

[16] 包紅軍，趙琳娜．基于Kalman濾波糙率反演模型的河道洪水實(shí)時(shí)預(yù)報(bào)研究[J]．水力發(fā)電學(xué)報(bào)，2012，31(3)：59 - 64．

[17] 陳素紅，劉孟凱，邢領(lǐng)航．基于多親遺傳算法的渠道非恒定流糙率率定模型及其應(yīng)用[J]．水電能源科學(xué)，2013，31(5)：81 - 83．

[18] 陳一帆，程海洋，萬(wàn)曉麗，等．結(jié)合糙率校正的河網(wǎng)水情數(shù)據(jù)同化[J]．水科學(xué)進(jìn)展．2015，26(5) ：731 - 738．

（責(zé)任編輯 姚小槐）

Inverse Problems of River Roughness: Research Review and Prospect

SHEN Wu - wei1，2，CHEN Yi - fan1，2，SHEN Zhen - dong1，2
(1. Zhejiang Institute of Hydraulics & Estuary，Hangzhou 310020，Zhejiang，China；2. Zhejiang Guangchuan Engineering Consultation Co., Ltd., Hangzhou 310020，Zhejiang，China)

River roughness is a critical sensitive parameter of hydraulic calculation for rivers，which is a comprehensive one inf l uenced by many factors. Table - checking，formula，trial - and - error and other commonly used methods cannot meet the demands of hydraulic calculation. Therefore，the inversion theories and methods of river roughness have drawn extensive attention of scholars from home and abroad. After reviewing and summarizing the basic methods and research updates of inverse methods for river roughness，this paper brought forward the difficulties and challenges as well as prospective directions for future study in this fi eld.

river roughness；inversion analysis；roughness calibration

TV131.4

A

1008 - 701X(2017)04 - 0001 - 03

10.13641/j.cnki.33 - 1162/tv.2017.04.001

2016-02-21

國(guó)家自然科學(xué)基金資助項(xiàng)目(51609213)。

沈五偉(1983 - )，男，工程師，碩士，主要從事水文水動(dòng)力學(xué)研究、水利工程規(guī)劃設(shè)計(jì)工作。

E - mail：214884860@qq.com