洪泛區(qū)河道行洪能力數(shù)值模擬研究

呂 敏（遼寧省丹東水文局，遼寧丹東118000）

?

洪泛區(qū)河道行洪能力數(shù)值模擬研究

呂 敏
（遼寧省丹東水文局，遼寧丹東118000）

由于我國的防洪措施不完善，人口密度相對較大，洪水發(fā)生時會造成嚴(yán)重的損失。如何合理的對河道進(jìn)行設(shè)計，計算河道的洪水水位，預(yù)測洪水的演進(jìn)路線和破壞力，對抗洪搶險有著重要意義。隨著數(shù)值模擬在水利工程中的應(yīng)用，對水力特性研究和洪水演進(jìn)提供了較好的模擬數(shù)據(jù)。以公路橋到水利樞紐間的44km河道為研究對象，首先建立了淺水運(yùn)動控制方程，并對該方程進(jìn)行離散化。隨后采用MIKE11HD和HEC-RAS兩種差分格式對洪泛區(qū)河道行洪能力進(jìn)行數(shù)值模擬，分析了兩種軟件計算水位情況。研究表明：由于兩種軟件的水力半徑不同，HEC-RAS軟件的計算水位較高。

洪泛區(qū)；河道；行洪能力；數(shù)值模擬

DOI：10.3969 /j.issn.1672-2469.2016.01.011

洪水是由于急劇降水、急驟冰雪融化、風(fēng)暴引起的江河湖海水位迅速提升的水流現(xiàn)象。研究表明：洪水具有流量大、流速高、來勢兇猛、破壞力強(qiáng)等特點(diǎn)。由于我國的防洪措施不完善，人口密度相對較大，洪水發(fā)生時會造成嚴(yán)重的損失。泄洪是水利建設(shè)中的重要技術(shù)問題，水利工程的修建勢必引起河道水流特性的改變，如果不采取有效措施，將會沖刷下游河床，危及水利樞紐的安全。如何合理的對河道進(jìn)行設(shè)計，計算河道的洪水水位，預(yù)測洪水的演進(jìn)路線和破壞力，對抗洪搶險有著重要意義。

由于數(shù)值仿真較傳統(tǒng)物理模型更為簡單，修改模型十分方便。隨著數(shù)值模擬在水利工程中的應(yīng)用，對水力特性研究和洪水演進(jìn)提供了較好的模擬數(shù)據(jù)。目前，國內(nèi)外都會對洪泛區(qū)進(jìn)行洪水風(fēng)險分析，并對一定洪水位下的淹沒面積進(jìn)行數(shù)值計算。以公路橋到水利樞紐間的44km河道為研究對象，首先建立了淺水運(yùn)動控制方程，并對該方程進(jìn)行離散化。隨后采用MIKE11HD和HEC-RAS兩種差分格式對洪泛區(qū)河道行洪能力進(jìn)行數(shù)值模擬，分析了兩種格式計算水位情況。希望為今后的洪水水位計算提供參考。

1 淺水運(yùn)動控制方程

在數(shù)值模擬普遍之前，計算河道行洪能力的方法主要有：水文學(xué)法、系統(tǒng)方法、水力學(xué)法。水文方法對初始條件的要求十分簡單，僅需要洪水的測量數(shù)據(jù)，不需要河道的結(jié)構(gòu)參數(shù)，應(yīng)用最為廣泛。系統(tǒng)學(xué)方法是將河道內(nèi)的水流視為一個完整的系統(tǒng)，將進(jìn)水量作為系統(tǒng)入口，出流量作為輸出結(jié)構(gòu)，河道則為系統(tǒng)的作用過程。水力學(xué)法主要依據(jù)圣維南定理得出，物理基礎(chǔ)很強(qiáng)。上世紀(jì)50年代，計算水力學(xué)誕生，其將水力學(xué)和數(shù)值計算相結(jié)合，對流體力學(xué)的求解起到了很大的促進(jìn)。

為了計算研究河道的行洪能力，首先建立了淺水運(yùn)動控制方程，并對該方程進(jìn)行離散化。淺水運(yùn)動控制方程主要針對不可壓、粘性流體，采用納維斯-托克斯方程進(jìn)行建模。由于淺水運(yùn)動十分復(fù)雜，在進(jìn)行建模時需滿足以下假設(shè)：

（1）洪泛區(qū)河道的水體水平速度要比垂直速度大的多，因此忽略垂直方向的加速度，認(rèn)為垂直方向的壓強(qiáng)分布為靜水壓強(qiáng)分布。

（2）由于淺水運(yùn)動中，水流的密度變化不大，較大的密度變化只在特殊情況出現(xiàn)，因此動量方程中認(rèn)為密度恒定。

（3）在正常溫度和壓力下，流體的彈性模量基本恒定，忽略其壓縮性，認(rèn)為其為不可壓流體。

由基本假設(shè)可知，流體速度沿垂向均勻分布，當(dāng)水流在橫向上分布均勻時，即可采用一維運(yùn)動方程進(jìn)行求解，將二維淺水運(yùn)動控制方程沿著y軸積分，并按照河道寬度取平均值，即可得到一維淺水運(yùn)動控制方程。

連續(xù)性方程：

動量方程：

其中：Q為流量值；A為河道截面積；x為橫坐標(biāo)；t為時間坐標(biāo)；h為洪水位；q為旁側(cè)入流流量；C、n為系數(shù)；R為水力半徑。

MIKE11HD差分格式采用Q和H交叉計算的方法對一維淺水運(yùn)動控制方程進(jìn)行離散，六點(diǎn)隱式方法進(jìn)行離散，該方法將H和Q進(jìn)行交替計算，計算示意圖見圖1。

圖1 MIKE11HD差分格式離散方法

HEC-RAS差分格式采用四點(diǎn)隱式方法進(jìn)行離散，這種方法主要針對變量f，在因變量一階偏微分方程的基礎(chǔ)上，采用相鄰點(diǎn)間加權(quán)平均離散方法，對相同時間層t的偏微分取j和j+1兩個插商值，對相鄰點(diǎn)x的偏微分取t和t+1插商進(jìn)行平均。HEC-RAS差分格式離散方法見圖2。

圖2 HEC-RAS差分格式離散方法

2 行洪能力數(shù)值模擬

2.1 模型建立

在模擬之前首先建立數(shù)學(xué)模型。在模型建立前，首先以河道參數(shù)為計算基礎(chǔ)，結(jié)合圣維南定理，確定相應(yīng)的邊界條件和初始條件。網(wǎng)格劃分在GAMBIT軟件中完成，并設(shè)置好邊界條件和入口參數(shù)，為了保證計算精度，對網(wǎng)格進(jìn)行加密處理，橫斷面間距為0.1km。河道的地形數(shù)據(jù)可以反映出沿程斷面的變化規(guī)律，據(jù)此得出計算區(qū)域水流量和水位的關(guān)系。計算采用分離式求解器，并進(jìn)行壓力場與速度場的耦合。

2.2 求解結(jié)果

采用兩種差分格式對計算河道內(nèi)的渡河大橋進(jìn)行行洪計算，兩種格式的計算方法都建立在能量方程的基礎(chǔ)上。采用半經(jīng)驗(yàn)公式計算渡河大橋橋墩處的水位壅水情況，并與兩種差分格式的計算結(jié)果進(jìn)行對比。HEC-RAS差分格式計算結(jié)果見圖3，MIKE11HD差分格式計算結(jié)果見圖4。

圖3 HEC-RAS差分格式計算結(jié)果

圖4 MIKE11HD差分格式計算結(jié)果

由數(shù)值計算結(jié)果可知，MIKE11HD格式計算的壅水高度為0.038m；HEC-RAS差分格式計算的壅水高度為0.041m；Henderson半經(jīng)驗(yàn)公式計算的壅水高度為0.027m；無坎寬頂堰半經(jīng)驗(yàn)公式計算的壅水高度為0.03m；規(guī)范公式計算的壅水高度為0.017m。兩種計算值與半經(jīng)驗(yàn)公式結(jié)果相近，分析產(chǎn)生誤差的原因，半經(jīng)驗(yàn)公式計

算時采用的比降為全局比降，而數(shù)值計算中采用的比降與每個斷面的水位高程、河底深度均相關(guān)，因此計算出的壅水高度較大。

為了便于對比兩種差分格式的計算誤差，圖5給出了44km河道的兩種模型計算結(jié)果對比。

圖5 兩種模型計算結(jié)果對比

由圖5可以看出，兩種差分格式得出的計算結(jié)果基本一致，只是在中部有點(diǎn)偏差。但是當(dāng)計算距離大于35km后，HEC-RAS差分格式的計算結(jié)果明顯大于MIKE11HD差分格式，最大偏差值為0.55m。分析產(chǎn)生偏差的原因：兩種格式的水力半徑不同，HEC-RAS軟件的計算水位較高。

3 結(jié)語

洪水具有流量大、流速高、來勢兇猛、破壞力強(qiáng)等特點(diǎn)。水利工程的修建勢必引起河道水流特性的改變，如果不采取有效措施，將會沖刷下游河床，危及水利樞紐的安全。通過建立淺水運(yùn)動控制方程，并對該方程進(jìn)行離散化，采用MIKE11HD和HEC-RAS兩種差分格式對洪泛區(qū)河道行洪能力進(jìn)行數(shù)值模擬。研究表明：

（1）兩種計算值與半經(jīng)驗(yàn)公式結(jié)果相近，但半經(jīng)驗(yàn)公式計算時采用的比降為全局比降，而數(shù)值計算中采用的比降與每個斷面的水位高程、河底深度均相關(guān)，因此計算出的壅水高度較大。

（2）兩種差分格式得出的計算結(jié)果基本一致，但是當(dāng)計算距離大于35km后，HEC-RAS格式結(jié)果明顯大于MIKE11HD格式，這是由于兩種格式的水力半徑不同。

［1］張驗(yàn)科，王麗萍，劉方，等.模擬最大熵法及其在河道行洪風(fēng)險計算中的應(yīng)用［J］.中國農(nóng)村水利水電，2011（02）：121-124.

［2］蘇崇峰.沂河華沂至鐵路橋段河道治理方案研究［J］.水利技術(shù)監(jiān)督，2011（04）：51-53 +65.

［3］姜樹海.基于隨機(jī)微分方程的河道行洪風(fēng)險分析［J］.水利水運(yùn)科學(xué)研究，1995（02）：127-137.

［4］郭夢京，李婧，周孝德，等.水動力模擬在渣場河道行洪分析中的應(yīng)用［J］.水土保持通報，2013（02）：156-159.

［5］楊婷.淺析城市水土流失與行洪排澇——以珠江三角洲城市為例［J］.水利技術(shù)監(jiān)督，2013（04）：31-32 +44.

［6］李計龍.瑪納斯河河道采砂規(guī)劃分析［J］.水利規(guī)劃與設(shè)計，2015（03）：11-12.

［7］余彥群，楊曉梅.淮河干流王家壩至臨淮崗段河道治理［J］.水利規(guī)劃與設(shè)計，2013（11）：62-65.

［8］賀芳丁，馮毅，張濤.尚莊水庫除險加固工程溢洪道設(shè)計方案分析［J］.水利技術(shù)監(jiān)督，2011（05）：45-47.

TV135

A

1672-2469（2016）01-0033-02

2015-06-19

呂 敏（1963年—），女，工程師。