周　維,郜會(huì)彩,吳　騰

(1.貴州省水利水電勘測(cè)設(shè)計(jì)研究院 夾巖設(shè)計(jì)處，貴陽(yáng)　550002 ；2.紹興文理學(xué)院 土木工程學(xué)院, 浙江 紹興　310020；3.河海大學(xué) 港口海岸與近海工程學(xué)院，南京　210098)

?

山區(qū)河流二維推移質(zhì)數(shù)學(xué)模型及應(yīng)用研究
——以望謨河整治工程為例

周維1,郜會(huì)彩2,吳騰3

(1.貴州省水利水電勘測(cè)設(shè)計(jì)研究院 夾巖設(shè)計(jì)處，貴陽(yáng)550002 ；2.紹興文理學(xué)院 土木工程學(xué)院, 浙江 紹興310020；3.河海大學(xué) 港口海岸與近海工程學(xué)院，南京210098)

為了充分利用我國(guó)西南山區(qū)豐富的水能資源，同時(shí)減少望謨河因山洪引起的高含沙水流洪澇和泥石流災(zāi)害，應(yīng)用二維水流泥沙數(shù)學(xué)模型常用的控制方程式，采用有限體積法對(duì)方程進(jìn)行離散，并使用“凍結(jié)法”對(duì)模型中陸地邊界進(jìn)行處理，建立了具有一定通用性的山區(qū)河流的二維推移質(zhì)數(shù)學(xué)模型。利用建立的模型對(duì)貴州省望謨河新屯段進(jìn)行了流場(chǎng)和推移質(zhì)沖淤模擬計(jì)算，并進(jìn)行了望謨河整治工程影響分析。結(jié)果表明：對(duì)于望謨河水位和河床沖淤變形驗(yàn)算，數(shù)值模型計(jì)算值與經(jīng)驗(yàn)公式推算變化趨勢(shì)一致，數(shù)值也較為接近；整治工程實(shí)施后望謨河沿程的水位絕大部分發(fā)生不同程度的降低，僅局部存在微小的抬高；望謨河整治后，同流量下，水位降低，流速歸槽泥沙淤積減小，有利于河道防洪。

山區(qū)河流；二維推移質(zhì)；數(shù)學(xué)模型；沖淤變化;河道整治工程

1　研究背景

推移質(zhì)輸移是河流動(dòng)力學(xué)研究的一個(gè)關(guān)鍵問題，也是研究推移質(zhì)運(yùn)動(dòng)和河床沖淤規(guī)律的理論基礎(chǔ)，同時(shí)又是解決河床演變計(jì)算、物理模型試驗(yàn)以及數(shù)學(xué)模擬計(jì)算等問題不可缺少的工具。國(guó)內(nèi)外許多學(xué)者對(duì)推移質(zhì)進(jìn)行了研究，據(jù)統(tǒng)計(jì)，目前各種推移質(zhì)輸沙率公式已超過(guò)50個(gè)。梅葉-彼得公式是以拖曳力為主要參數(shù)的典型公式[1-2]，該公式吸收了大量的試驗(yàn)資料；R.A.Bagnold[3]、張瑞瑾等[4]、張紅武等[5]從能量平衡的觀點(diǎn)出發(fā)，根據(jù)推移質(zhì)輸沙特性所得公式，適用于小至細(xì)沙,大至卵石的低、高輸沙強(qiáng)度的推移質(zhì)輸沙率計(jì)算；以岡卡洛夫?yàn)榇淼那疤K聯(lián)學(xué)者建立的一系列公式[6]和我國(guó)成都勘測(cè)設(shè)計(jì)院的推移質(zhì)輸沙率公式[7]，以流速為主要參數(shù)，此類公式中輸沙率與流速4次方成比例，對(duì)流速變化十分敏感。過(guò)去幾十年來(lái)，國(guó)內(nèi)外許多學(xué)者對(duì)已有的輸沙率公式進(jìn)行了驗(yàn)證和比較，但關(guān)于推移質(zhì)輸沙率公式的比較尚無(wú)定論。一方面，對(duì)推移質(zhì)運(yùn)動(dòng)過(guò)程的認(rèn)識(shí)并不清晰；另一方面，局限于當(dāng)前推移質(zhì)的測(cè)量手段，尚無(wú)較為準(zhǔn)確的推移質(zhì)實(shí)測(cè)資料，導(dǎo)致推移質(zhì)輸沙率公式中的參數(shù)并不適用于天然河道，因此，給實(shí)際應(yīng)用帶來(lái)了巨大的困難。

望謨河位于貴州省望謨縣境內(nèi)，歷來(lái)洪澇災(zāi)害頻繁。2011年望謨縣 “6·06”特大洪水過(guò)后，河道及兩岸淤積十分嚴(yán)重，河道部分區(qū)域淤積厚度達(dá)2～3 m，降低了河道的過(guò)洪能力；洪水期水位抬高，水流漫出河道，導(dǎo)致兩岸農(nóng)田及建筑物損失嚴(yán)重。為保障人民生命財(cái)產(chǎn)安全和地方經(jīng)濟(jì)穩(wěn)定發(fā)展，近年實(shí)施了望謨河河道治理工程。整治工程修建后在一定程度上改變河道原有的水沙運(yùn)動(dòng)及河床變形過(guò)程。由于望謨河為山區(qū)河流，河床以沙卵石為主，為推移質(zhì)造床，因此，有必要加深對(duì)推移質(zhì)輸移及河道演變規(guī)律的理解和認(rèn)識(shí)，為防災(zāi)減災(zāi)決策制定提供可靠支撐。

本文根據(jù)山區(qū)河流的特點(diǎn)，建立了具有一定通用性的山區(qū)河流二維推移質(zhì)數(shù)學(xué)模型，利用建立的模型對(duì)貴州省望謨河新屯段進(jìn)行了流場(chǎng)和推移質(zhì)沖淤模擬計(jì)算。

2　模型的建立

2.1控制方程

二維水流泥沙數(shù)學(xué)模型可用于描述水沙各運(yùn)動(dòng)要素沿水深的平均值、河床高程隨流程沿河寬2個(gè)方向的變化情況，基于控制體概念并遵循質(zhì)量守恒定律和牛頓第二定律，可導(dǎo)出二維水流泥沙數(shù)學(xué)模型經(jīng)常采用的控制方程式，即

(1)水流連續(xù)方程：

(1)

式中：Z為水位；t為時(shí)刻；h為水深，h=Z-Zb，Zb為河底高程；u，v分別為垂線平均流速沿x，y方向的分量。

(2)水流運(yùn)動(dòng)方程：

(2)

(3)

式中：f為柯氏力；g為重力加速度；νt表示紊動(dòng)黏性系數(shù)；n為糙率系數(shù)。

(3)懸移質(zhì)不平衡輸沙基本方程：

(4)

(4)推移質(zhì)不平衡輸沙基本方程：

本次計(jì)算采用竇國(guó)仁導(dǎo)出的推移質(zhì)不平衡輸沙方程，即

(5)

(5)河床變形方程(不平衡輸沙法)：

(6)

2.2方程的離散

本模型采用有限體積法對(duì)方程進(jìn)行離散，其基本思想是：將計(jì)算區(qū)域劃分為若干規(guī)則或不規(guī)則相互連接但不相互重疊的單元，將變量定義在控制體形心，對(duì)基本方程在每個(gè)控制體內(nèi)積分。得出一組以計(jì)算節(jié)點(diǎn)上物理量為未知數(shù)的代數(shù)方程組來(lái)求解。

2.3邊界處理

陸地邊界的處理方法是：令閉邊界的法向流速為0，而沿切線方向的流速為非0值。隨著水位的變化，陸地邊界位置不斷變化。本模型采用“凍結(jié)法”處理這種問題，根據(jù)計(jì)算水位和湖底高程判斷網(wǎng)格單元是否露出水面，即定義臨界水深h=0.005～0.01 m，當(dāng)水深h>h時(shí)，糙率取正常值，反之糙率取一大值(1010量級(jí))。

3　模型的驗(yàn)證

3.1典型條件下數(shù)學(xué)模型的驗(yàn)證

實(shí)際工程中潰壩大部分出現(xiàn)局部、部分的潰決。部分潰壩有著間斷流，急流，非恒定流等特性，是非常經(jīng)典的例子。本文研究的模型尺寸如圖1所示，壩的中間有75 m的潰決口，初始水位為10 m。圖2為部分潰壩后7 s時(shí)的水面圖。由圖2可以看出，上游兩側(cè)水流均向中間流動(dòng)，在水面中間處形成低谷。隨著波不斷向上傳播，低谷也不斷往上游移動(dòng)。此外，在潰口兩側(cè)拐點(diǎn)處產(chǎn)生明顯回流。

圖1　潰壩模型尺寸Fig.1　Size of dam-break model

圖2　潰壩后t=7 s水面Fig.2　Water surface at 7 s after dam break

為了進(jìn)一步驗(yàn)證算法，把本文計(jì)算結(jié)果與S.Fagherazzi等[8]計(jì)算的結(jié)果進(jìn)行比較，如圖3所示?？傮w而言，兩者計(jì)算值較為接近，可以模擬潰壩的流態(tài)。

圖3　水面等勢(shì)線計(jì)算結(jié)果對(duì)比Fig.3　Comparison of calculated results of water surface contours

3.2望謨河的驗(yàn)證

3.2.1水位驗(yàn)證

本次計(jì)算選取的范圍為望謨河新屯至東巖段設(shè)計(jì)防洪治理河道長(zhǎng)度為1.647 km的整治河段，網(wǎng)格數(shù)為289×95，網(wǎng)格長(zhǎng)度為7 m，計(jì)算過(guò)程中時(shí)間步長(zhǎng)取10 s。由于望謨河整治河段內(nèi)無(wú)沿程的實(shí)測(cè)水位、流速等數(shù)據(jù)，因此，本次驗(yàn)證采用經(jīng)驗(yàn)公式推算的水位與本模型計(jì)算的水位進(jìn)行類比(該成果已通過(guò)專家驗(yàn)收)。圖4(a)、圖4(b)分別為5%頻率洪水驗(yàn)證水位和2011年洪水沿程水位，可以看出，數(shù)值模型計(jì)算值與經(jīng)驗(yàn)公式推算值變化趨勢(shì)一致，數(shù)值也較為接近，初步說(shuō)明本文建立的模型能應(yīng)用于望謨河的研究。

圖4　 經(jīng)驗(yàn)公式和數(shù)值模型計(jì)算水位值對(duì)比Fig.4　Comparison of water levels obtained from empirical formula and simulation

3.2.2河床沖淤變形驗(yàn)證

根據(jù)搜集到的實(shí)測(cè)資料，望謨河2011年洪水過(guò)后河道內(nèi)泥沙淤積厚度為2～3 m，限于實(shí)測(cè)資料，驗(yàn)證計(jì)算模擬該場(chǎng)次洪水后河道的沖淤變化。圖5為數(shù)值模型計(jì)算的2011年洪水后望謨河河道沖淤變化分布，由淤積分布可知，該洪水條件下河段的絕大部分發(fā)生淤積，淤積厚度為2～3 m，局部發(fā)生沖刷，沖刷范圍和幅度均較小。計(jì)算結(jié)果表明，該模型也可以初步運(yùn)用于模擬望謨河的河床沖淤變形。

圖5　2011年洪水后望謨河河道沖淤變化Fig.5　Scouring and deposition of Wangmo river after the flood period in 2011

4　望謨河整治工程影響分析

4.1數(shù)值模型計(jì)算條件

為了更好分析望謨河整治后的效果，選取10 a一遇洪水，計(jì)算分析整治工程對(duì)水位、流速的影響。數(shù)值模型計(jì)算的條件為：上游來(lái)流量400 m3/s，下游水位為631.58 m。圖6為數(shù)值模型水位測(cè)點(diǎn)及測(cè)量斷面分布圖。

圖6　工程區(qū)測(cè)點(diǎn)及測(cè)量斷面分布Fig.6　Layout of observation points and observation sections

4.2望謨河整治對(duì)水位的影響

表1為望謨河整治河段工程實(shí)施前后測(cè)點(diǎn)水位變化。由表1可知，10 a一遇洪水條件下，整治工程后沿程絕大部分河段的水位發(fā)生不同程度降低，僅在D2斷面和D3斷面水位存在微小的抬高，分別抬高0.06 m和0.01 m。具體而言，進(jìn)口斷面至D3斷面水位差有逐漸減小趨勢(shì)，從D3斷面至D5斷面水位差有逐漸增大趨勢(shì)，D1斷面的最大水位差為0.72 m，D2斷面最大水位差0.21 m，同時(shí)該斷面部分測(cè)點(diǎn)水位差出現(xiàn)負(fù)值，D3斷面整治前后水位變化不大，D4斷面及其下游水位差增大迅速，水位差最大值出現(xiàn)在D4-1斷面，為2.44 m。該計(jì)算結(jié)果表明，望謨河整治后，由于河槽變深，同水位下過(guò)水面積增大，水位降低，望謨河整治有利于河道防洪。

4.3望謨河整治對(duì)流態(tài)的影響

圖7為10 a一遇洪水條件下整治前后流速比較。由圖7可以看出，D1斷面上游，河道未進(jìn)行疏浚，整治前后流速和流向變化不大，但整治后有向主河槽集中的趨勢(shì)；D1至D1-2河段，主河槽流速增大，灘地流速減?。籇1-2至D1-3河段，整治后流速較整治前流速減小，主要由于河道疏浚后，河槽過(guò)水面積增大導(dǎo)致流速降低，但水流均有歸槽趨勢(shì)；D1-3至D3-1河道內(nèi)，灘面流速歸槽明顯，尤其是在D2河段附近，未疏浚時(shí)有明顯水流，但河道疏浚后，水流均在河槽中運(yùn)動(dòng)，水流未漫灘，有利于汛期的防洪；D3-1至D5-1河道內(nèi)流速變化不大，流向偏向疏浚的河槽；D5-1斷面下游，水流流向變化不明顯，流速稍有減小。工程區(qū)域內(nèi)，整治前最大流速出現(xiàn)在D2斷面，流速大小為4.33 m/s，整治后最大流速出現(xiàn)在D5斷面，最大流速為4.08 m/s；整治后與整治前流速差異最大值出現(xiàn)在D2斷面，流速差為2.86 m/s，主要由于主流歸槽的影響。

表1　整治前、后測(cè)點(diǎn)水位變化

注：水位差為整治前水位減整治后水位，“—”為不過(guò)水。

注：實(shí)心箭頭為整治前，空心箭頭為整治后圖7　10 a一遇洪水條件下整治前、后流速比較Fig.7　Comparison of flood velocity before and after river project under flood probability of 10%

4.4望謨河整治后洪水條件下河道演變分析

根據(jù)望謨河的實(shí)測(cè)洪水資料分析，本次計(jì)算模擬洪水的歷時(shí)為5 h，河床的中值粒徑為5 cm。圖8為修建工程后遇2011年洪水河道的沖淤分布。

圖8　河道整治后遇2011年洪水河道沖淤分布Fig.8　Scouring and deposition of Wangmo river encountered with flood in 2011 after river regulation

由圖8可知，該條件下望謨河道主要以淤積為主，局部區(qū)域存在沖刷。淤積厚度的變化規(guī)律為：洪水頻率越大，望謨河泥沙淤積的厚度也越大，尤其是D1斷面至D2斷面間以及D4斷面至D5-1斷面間河道淤積較大。該變化規(guī)律的主要原因是該河段的淤積主要由于上游來(lái)的推移質(zhì)泥沙，當(dāng)洪水越大時(shí)上游的來(lái)沙量越大，當(dāng)流速稍有減小則可能發(fā)生淤積，故在該區(qū)域的淤積也會(huì)越大。此外，根據(jù)計(jì)算結(jié)果也表明，當(dāng)發(fā)生5 a一遇洪水時(shí)，該計(jì)算區(qū)域的泥沙淤積量并不大，僅在河道的D1斷面、D4斷面和D5-1斷面產(chǎn)生較大的淤積，其余區(qū)域則淤積較小。河道整治后，河道淤積量較整治前(圖5)明顯減小。

5　結(jié)　語(yǔ)

基于山區(qū)河流的水沙運(yùn)動(dòng)特性，建立了適用于望謨河水沙數(shù)學(xué)模型，并采用該模型對(duì)望謨河整治工程的效果進(jìn)行了計(jì)算分析。計(jì)算結(jié)果表明：

(1)整治工程實(shí)施后望謨河沿程的水位絕大部分發(fā)生不同程度的降低，僅在D2斷面至D3斷面水位存在微小的抬高。水位差的變化規(guī)律為：D1斷面至D3斷面水位差有逐漸減小趨勢(shì)，從D3-1斷面至D5-1斷面水位差有逐漸增大趨勢(shì)。水位差最大值均出現(xiàn)在D4-1斷面，10 a一遇洪水條件下整治后水位降低最大為2.44 m。

(2)10 a一遇洪水條件下，D1斷面上游，整治前后流速大小差異不大；D1至D1-3河段，整治后流速較整治前流速稍?。籇1-3至D3-1河道內(nèi)，灘面流速歸槽明顯，尤其是在D2河段附近，整治前右側(cè)過(guò)流，整治后右側(cè)不過(guò)流，該條件下整治后洪水的淹沒范圍較整治前減小較大；D3-1至D5-1河段內(nèi)流速變化不大，水流流向疏浚的河槽；D5-1斷面下游，水流流向變化不明顯，流速稍有減小，同時(shí)灘地上的過(guò)流量減小。

(3)望謨河整治后，同流量下，水位降低、水流歸槽泥沙淤積減小，有利于河道防洪。

[1]MEYER-PETER E ,MULLER R.Formulas for Bed Load Transport[C]∥International Association of Hydraulic Research .Proceedings of 2nd Meeting of the International Association for Hydraulic Structures Research.Delft,Netherlands,June 7,1948:39-64.

[2]MEYER-PETER E,FAVRE H,EINSTEIN H A.Neuere Versuchsresutate Uber den Geschiebetrieb[J].Schweizerische.Bauzeitung,1934,103(12):147-150.

[3]BAGNOLD R A.An Approach to the Sediment Transport Problem from General Physics[R].Washington,USA:United States Government Printing Office ,1996.

[4]張瑞瑾，謝鑒衡．河流泥沙動(dòng)力學(xué)[M]．北京：水利電力出版社，1998.

[5]張紅武,張俊華,卜海磊,等.試論推移質(zhì)輸沙率公式[J].南水北調(diào)與水利科技,2011，9(6)：140-145.

[6]謝鑒衡．河流泥沙工程學(xué)(上冊(cè))[M]．北京：水利電力出版社,1982.

[7]李昌華，金德春．河工模型試驗(yàn)[M]．北京：人民交通出版社，1981.

[8]FAGHERAZZI S,RASETARINERA P,HUSSAINI M Y,et al.Numerical Solution of the Dam-Break Problem with a Discontinuous Galerkin Method[J].Journal of Hydraulic Engineering,2004,130(6):532-539.

(編輯：占學(xué)軍)

Application of 2-D Numerical Model of Bed Load in Riversof Mountain Area：Regulation Project in Wangmo River as An Example

ZHOU Wei1,GAO Hui-cai2,WU Teng3

(1.Design Department of Jiayan Project ,Guizhou Investigation and Design Institute of Water Conservancy and Hydropower,Guiyang550002,China; 2.College of Civil Engineering,Shaoxing University,Shaoxing 310020,China; 3.College of Harbour,Coastal and Offshore Engineering,Hohai University, Nanjing210098,China)

In order to utilize rich hydropower resources in southwest mountain area of China and decrease flow with high sand proportion and debris flow disasters caused by flash flood,we should establish suitable two-dimensional bed load mathematical models for rivers in mountain area.On the basis of finite volumetric method,we dispersed common controlling equations in sediment numerical model of two-dimensional flow,dealt with dry land boundary by using frozen method,and established a general two-dimensional bed load mathematical model.Xintun section in Wangmo river located in Guizhou province was taken as an example.We carried out numerical calculation of flow field,exploring scouring and deposition of bed load.Test results showed that,as for water level calculation and verification of scouring and deposition in the river bed,simulated data was in consistent with data from empirical formula; after regulation project of the river,most of water levels along the flow direction decreased,except for slight increase in few sections; meanwhile,in the condition of given discharge，water level after the regulation was smaller than that before the regulation,and sediment deposition into the water flume decreased,which is in favor of river flood control.

river in mountain area; two-dimensional bed load; numerical model; scouring and deposition; river regulation project

2015-06-04;

2015-08-19

國(guó)家自然科學(xué)基金項(xiàng)目(51309084);貴州省水利科技項(xiàng)目(KT201303);浙江省公益技術(shù)應(yīng)用研究計(jì)劃項(xiàng)目(2015C31134)

周維(1982-)，男，貴州仁懷人，高級(jí)工程師，主要從事水工建筑物設(shè)計(jì)和河流泥沙研究，(電話)0851-85584296(電子信箱)287043164@qq.com。

郜會(huì)彩(1979-)，女，河南鄢陵人，高級(jí)工程師，碩士，主要從事河流泥沙和泥石流研究，(電話)15215998120(電子信箱)gaohuicai@sina.com。

10.11988/ckyyb.201504732016,33(08):6-10

TV92

A

1001-5485(2016)08-0006-05