基于糙率參數不確定性的河道防洪能力風險分析*

殷 丹 石鳳君 王 凱

(遼寧省水利水電科學研究院， 遼寧 沈陽 110003)

基于糙率參數不確定性的河道防洪能力風險分析*

殷 丹 石鳳君 王 凱

(遼寧省水利水電科學研究院， 遼寧 沈陽 110003)

為減少糙率參數不確定性給河道水面線推求帶來的影響，提高河道防洪能力分析成果可靠性，提出基于參數不確定性的防洪能力風險分析技術。以柴河水庫下游段為例，采用改進的摩爾斯分類篩選法進行糙率靈敏度分析，得出不同河段的糙率靈敏度，為模型計算中的糙率值選取提供依據。通過對柴河現狀河道防洪能力以及堤段設計洪水標準下的河道行洪風險進行計算分析，得出行洪風險高的堤段，為防洪應急預案編制提供參考。

糙率； 不確定性； 防洪能力； 風險分析

河道防洪能力分析是根據不同洪水流量條件下的洪水淹沒范圍和淹沒水深，得出各河段某一保證水位下可以安全下泄的洪水流量。河道防洪能力分析結果決定著堤防高程的確定，同時也關系著堤防保護區的安全及其防洪標準的確定，是影響流域防洪規劃的重大問題。隨著橋梁、閘壩等涉河工程的建設以及采砂等人類活動，河道形態不斷發生變化，從而使河道的防洪能力受到不同程度的影響。在河道防洪能力分析計算中，控制斷面的初始水位、河道糙率、河道幾何形狀等輸入條件存在一定的不確定性，傳統的河道水面線推求未考慮各種隨機因素，而以定值方法進行計算，使河道水面線的不確定性難以得到定量的描述。河道糙率n是洪水計算分析中的一個重要參數，其取值是河道一維數值模擬的關鍵，糙率n取值準確與否與直接影響著水動力模型的計算精度，糙率率定是模型計算的重要步驟。本文以柴河為例，將糙率參數的不確定性影響引入水面線的推求過程，進行局部參數靈敏度分析，以降低河道行洪風險，提高河道防洪能力計算結果可靠度。

1 研究區概況

柴河為遼河干流的左側支流，發源于遼寧省清原縣枸乃甸鄉北樂山天橋嶺，于鐵嶺水文站上游處匯入遼河。流域面積1441km2，河流長度133km，全河比降為1.91‰。柴河上建有柴河水庫，水庫控制面積1355km2，占流域面積的94%。

柴河水庫下游屬低丘平原河流，地勢較平坦，主槽和灘地不明顯。柴河水庫下游是鐵嶺城市防洪的重點，河段兩岸村屯密布，人口聚集，特別是柴河入遼河河口段為鐵嶺市郊，該次分析范圍選取柴河水庫溢洪道入柴河處至長大鐵路橋，河道長度11.4km，具體見下圖。柴河水庫下游段河道寬約60～270m，治理型式主要為堤防和護岸，其中柴河大橋下游左岸為鐵嶺市城區，屬遼河回水段，堤防設計標準100年一遇，右岸為50年一遇；柴河大橋至柴河橡膠壩屬于城郊段，防洪標準50年一遇；柴河橡膠壩以上為農村段，基本無堤，防洪標準為20年一遇。隨著河道采沙以及城市的建設發展，不同程度的影響河道過流能力。

2 河道糙率參數靈敏度分析

2.1 靈敏度分析方法

靈敏度分析是模型參數識別過程的重要步驟之一，對于驗證及改進數值模型有重要的作用。在數值模型計算中，靈敏度分析可識別參數靈敏度的高低，定性或定量地評價參數不確定性對模擬結果的影響程度，能有效降低高靈敏度參數選擇不當產生的誤差，提高計算精度。參數靈敏度分析方法包括全局靈敏度分析和局部靈敏度分析，其中局部靈敏度分析主要用于計算單個參數的變化對模擬結果的影響，計算時僅改變一個固定參數的值，其他參數保持不變[1]。局部靈敏度分析只針對某一個參數，操作簡單、易于實施。本研究重點介紹糙率參數在模型計算中的靈敏度，為局部靈敏度分析，即檢驗糙率在最佳估計值附近的微變化對模擬結果的影響。

摩爾斯分類篩選法(morris screening method)在局部靈敏度分析中應用較為廣泛，方法簡單有效[2]。修正的摩爾斯篩選法，自變量以固定步長變化，靈敏度判別因子取多個平均值[3]，計算公式如式(1)所示。

(1)

式中S——靈敏度判別因子；

n——模型運行次數；

yi——控制斷面在模型第i次運行時的水位值；

y0——水位初值；

yi+1——模型第i+1次運行的水位值；

Pi——第i次模型運算參數值變化百分率；

Pi+1——第i+1次模型運算的參數值的變化百分率[4]。

S值越大，則該參數對計算結果越靈敏，影響越大。參照靈敏度分級標準，分為四級：0≤|S|<0.05為不靈敏參數；0.05≤|S|<0.2為中等靈敏參數；0.2≤|S|<1為靈敏參數；|S|>1為高靈敏參數[5]。

2.2 柴河糙率參數靈敏度分析

選用丹麥DHI水力模擬系列軟件的MIKE11模塊建立河道一維水動力模型，計算河道水面線，MIKE軟件在水力計算方面應用較為廣泛，其計算原理及過程在此不再贅述。該次計算橫斷面資料主要采用最新實測地形，柴河水庫下游河段計算長度11.4km，計算斷面20個。柴河水庫下游至河口段無較大支流匯入，有部分平原區及城區，洪水主要來自柴河水庫的放流，故該區段設計洪水直接采用柴河水庫設計情況下的泄流量。設計洪水采用《柴河水庫除險加固工程初步設計報告》中成果，見表1。下邊界采用起始斷面的水位流量關系。

為分析柴河糙率離散性對水位的影響，對糙率靈敏度進行分析。根據河道實際情況分段計算，分為上游鄉村段，溢洪道—熊官村；中游農田段，熊官村—沈四高速橋；下游城鎮段，沈四高速橋—長大鐵路橋等三段。

采用局部靈敏度分析法檢驗糙率在最佳估計值附近的微變化對模擬結果的影響。當糙率均值取0.03時，保持模型邊界條件、初始條件等其他因素不變，僅改變上、中、下游三段的糙率值，進行水面線計算。分別將糙率值以增大或減小10%、20%、30%、40%變化，計算河段沿程水位變化，并采用摩爾斯分類篩選法進行糙率靈敏度計算，結果見表2。

從計算結果可以看出，河道不同位置糙率值改變對水位計算結果影響不同，各河段的靈敏度順序依次

為|S|上游>|S|下游>|S|中游。根據靈敏度分級情況，柴河水庫溢洪道—熊官村、沈四高速橋—長大鐵路橋段的糙率值為靈敏參數，而中游的熊官村—沈四高速橋段的糙率為中等靈敏參數。柴河水庫溢洪道—熊官村的糙率靈敏，主要是由于河道兩岸多坑塘，河道地形變化較大；而下游城市段的糙率參數靈敏是由于城市河道灘槽明顯，灘地上的草木植物較多。對于糙率參數靈敏河段，在計算過程中要更加注意參數值的選取，減少計算步長，提高計算精度。

3 防洪能力風險分析

柴河長大鐵路橋至柴河大橋有堤防，設計防洪標準為50年一遇；柴河大橋以上現狀無堤防，分別進行50年一遇、20年一遇設計洪水的水位計算，以及現狀斷面保證流量計算，結果見表3。有堤段防洪能力基本達到50年一遇設計要求，由于無堤段兩岸地勢較高，泄流能力較強，該次計算中將其斷面保證流量確定為不超過100年一遇設計洪水流量，即1651m3/s。

由于河道行洪能力計算過程中參數不確定，使計算結果存在隨機性，該次研究對有堤防段進行河道行洪風險率計算。計算洪水頻率為50年一遇，根據堤防設計情況，堤防安全超高為2m。該次研究考慮糙率參數不確定性，計算發生設計洪水時河道行洪風險情況，結果見表4。

行洪風險率的數值為0～1，為方便比較，根據經驗確定行洪風險級別，將行洪風險率劃分為5個等級，見表5。

根據表4及表5對照分析，發生50年一遇洪水時，柴河長大鐵路橋斷面右岸堤防行洪風險級別較高，按正常糙率參數取值時，已達到設計防洪標準，但當發生特殊情況時，如大洪水將上游的樹木、莊稼或建筑物沖毀，將引起下游河道糙率大幅增加，引起水位上漲，有洪水漫溢或潰堤風險。

4 結 論

采用摩爾斯分類篩選法進行柴河水庫下游段的河道糙率靈敏度分析，得出柴河水庫溢洪道—熊官村以及沈四高速橋—長大鐵路橋段的糙率值為靈敏參數，而中游的熊官村—沈四高速橋段的糙率為中等靈敏參數。河道防洪風險分析結果表明，發生50年一遇設計洪水時，長大鐵路橋處右岸堤防的行洪風險率高，說明該段堤防在大洪水時，有洪水漫溢或潰堤風險，在汛期應注重該處河段的防汛檢查與應急搶險準備工作。

[1] 劉玉珍，程世迎.靈敏度分析法確定水文地質參數的基本模型及其應用[J].水利學報,2006,37(7):846-850.

[2] 薄會娟，董曉華，鄧霞.新安江模型參數的局部靈敏度分析[J].人民長江,2010,41(1):25-28.

[3] 張利茹，管儀慶，葉彬,等.新安江模型參數敏感性分析的實證研究[J].水電能源科學，2008，26(5)：16-17.

[4] 許春東，殷丹.基于MIKE11的河道糙率靈敏度分析[J].水電能源科學，2014，32(11)：101-103.

[5] Lenhart T,Eckhardt K,Fohrer N,et al.Comparison of two different approaches of sensitivity analysis[J].Physics and Chemistry of the Earth,2002,27(9-10):645-654.

Risk analysis of river flood control capacity based on roughness parameter uncertainty

YIN Dan, SHI Fengjun, WANG Kai

(LiaoningWaterConservancyandHydropowerResearchInstitute,Shenyang110003,China)

The flood control capacity risk analysis technology based on parameter uncertainty is proposed in order to reduce the influence of roughness parameter uncertainty on the investigation of river water surface profiles, and improve the reliability of river flood control capacity analysis results. The downstream section of Chaihe River is adopted as an example. Improved Morris screening method is adopted for analyzing roughness sensitivity. The roughness sensitivity in different river sections is obtained, thereby providing basis for selecting roughness value in model calculation. The embankment sections with high flood draining risks are discovered through calculating and analyzing current river flood control ability in Chaihe River and river flood draining risks under embankment design flood standard. Reference is provided for compiling flood control emergency preplans.

roughness; uncertainty; flood control capacity; risk analysis

10.16616/j.cnki.11- 4446/TV.2017.04.012

遼寧省科學技術計劃項目(2014212003)

TV82

A

1005-4774(2017)04- 0045- 04