張 超，常占鑫，馬小杰

(四川農(nóng)業(yè)大學(xué)水利水電學(xué)院，四川 雅安 625014 )

沖積河流的河床演變會(huì)受到諸多因素的影響，河型的成因以及轉(zhuǎn)化機(jī)理一直是學(xué)界研究和關(guān)注的重點(diǎn)[1]。受到流域乃至全球?qū)用嫒祟惢顒?dòng)的影響，河流系統(tǒng)的邊界條件可能會(huì)發(fā)生較大的變化，特別是一些大型水利工程的建設(shè)，勢(shì)必會(huì)導(dǎo)致下游河道水沙條件的變化，從而誘發(fā)河流的穩(wěn)定性以及河型的轉(zhuǎn)化[2]。例如，河流上游修建水庫可以大幅提升下游的防洪標(biāo)準(zhǔn)，同時(shí)提升流域內(nèi)水資源的配置和利用效率[3]。但是，下游的水沙過程和水沙條件也會(huì)發(fā)生較大的變化，而水沙條件的變化也會(huì)對(duì)河床演變?cè)斐芍匾绊?，進(jìn)而發(fā)生河型轉(zhuǎn)化，對(duì)下游河道航運(yùn)以及兩岸居民的生活造成潛在的影響[4]。因此，水沙條件對(duì)河型轉(zhuǎn)化的影響研究具有重要的理論和工程意義。基于此，文章通過建立二維水沙模型，利用數(shù)值模擬的方法對(duì)不同水沙條件下的河床縱向地形變化規(guī)律進(jìn)行研究，以便為相關(guān)的水利工程設(shè)計(jì)和河道工程建設(shè)提供理論支持。

1 研究河段的選擇

本次研究選擇的是含沙量較大的四川某山區(qū)河道中比較蜿蜒的一段，試驗(yàn)?zāi)Ｐ偷拈L度為300m，寬為200m，以模型的長度和彎曲程度為劃分依據(jù)，將模型試驗(yàn)河段劃分為14個(gè)不同的斷面，本次試驗(yàn)選擇其中的D3、D5、D7和D9斷面進(jìn)行研究，所有斷面的橫坐標(biāo)見表1。

表1 試驗(yàn)河段斷面橫坐標(biāo) 單位：m

(續(xù)表)

2 二維水沙模型的構(gòu)建

2.1 模型的控制方程

二維水沙模型的控制方程主要包括兩部分，一是二維水流基本方程，二是二維泥沙基本方程[5]。其中，二維水流方程主要由水流連續(xù)方程和水流動(dòng)量方程組成；二維泥沙方程主要包括泥沙連續(xù)方程與河床變形方程，具體表達(dá)式如下所示：

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

(1)

式中，h—水位高度，m；p、q—x方向和y方向上的流速分量，m/s；ξ—壓強(qiáng)，Pa。

(2)水流動(dòng)量方程

(2)

(3)

式中，H—河道水深，m；g—重力加速度，取9.8N/s2；τxx、τxy、τyy分別為有效剪切力分量，MPa；C—謝才系數(shù)；f—風(fēng)阻力系數(shù)；fw—無因次風(fēng)應(yīng)力系數(shù)；WWx—表面風(fēng)引起的x向阻力，N；WWy—表面風(fēng)引起的y向阻力，N。

(3)泥沙連續(xù)方程

(4)

式中，u、ν—x方向和y方向上的垂線平均流速分量，m/s；S—水體的平均含沙量，kg/m3；c0—泥沙的擴(kuò)散系數(shù)；S*—平均挾沙力；ω—泥沙的沉速，cm/s；a—恢復(fù)飽和系數(shù)。

(4)河床變形方程

(5)

式中，γ′—泥沙干容重，kg/m3；Z0—河底高程，m。

2.2 模型的建立

考慮到二維水沙模型構(gòu)建的具體要求，文章以選擇河段的原始地形為基礎(chǔ)進(jìn)行幾何模型的構(gòu)建[6]。首先以研究河段的CAD資料為基礎(chǔ)，利用ArcGIS技術(shù)建立起研究河段的原始地形圖。對(duì)構(gòu)建的幾何模型進(jìn)行必要的概化處理，以便提升模型運(yùn)算的效率和穩(wěn)定性。在構(gòu)建的幾何模型中提取水深點(diǎn)數(shù)據(jù)，然后插值為地形并導(dǎo)入MIKE21 FM 軟件，利用MIKE21 FM軟件完成計(jì)算模型圖的生成和網(wǎng)格剖分[7]。

在模型網(wǎng)格剖分過程中，需要對(duì)模型計(jì)算的效率和精度進(jìn)行平衡。網(wǎng)格劃分采用的是非結(jié)構(gòu)三角形網(wǎng)格，能夠較好地模擬深槽、淺灘等比較復(fù)雜的河道地形，有利于提高模型的計(jì)算精度。模型網(wǎng)格的尺寸在3～5m之間，最終，模型劃分為1243個(gè)計(jì)算單元。模型的網(wǎng)格剖分示意圖如圖1所示。

圖1 研究河段網(wǎng)格剖分示意圖

2.3 計(jì)算方案設(shè)計(jì)

研究中首先利用水流的攜沙力公式計(jì)算出河道的最大挾沙力，公式如下：

(6)

式中，S—水流最大挾沙力，kg/m3；H—平均水深，m；w—水力粗度，cm/s；U—水流流速，m/s；K，m—常數(shù)系數(shù)和指數(shù)。

結(jié)合研究河段的資料，利用上述公式計(jì)算獲取河道水流挾沙能力的最大值為12.0kg/m3，然后結(jié)合實(shí)際模型的尺寸以及模型的抗沖能力等因素，對(duì)水槽進(jìn)水流量的大小進(jìn)行設(shè)計(jì)，最終獲得的計(jì)算工況見表2。

表2 計(jì)算工況設(shè)計(jì)

3 試驗(yàn)結(jié)果與分析

要想獲取水沙條件對(duì)河道演變特征和規(guī)律的影響，首先主要獲得河床的演變特征，也就是河床的縱向和橫向變化特征。本次研究中，利用上節(jié)構(gòu)建的數(shù)值計(jì)算模型，計(jì)算不同方案下的河床縱向變化數(shù)據(jù)，結(jié)果顯示，不同方案下的河槽形狀也明顯不同。

3.1 G1方案計(jì)算結(jié)果與分析

利用上節(jié)構(gòu)建的模型對(duì)G1方案下研究河段的河床縱向變化特征進(jìn)行計(jì)算，根據(jù)模型計(jì)算獲得的數(shù)據(jù)，獲得不同時(shí)段的河床高程線圖，其中12h河床高程線圖如圖2所示。由圖2可知，在G1方案下，研究段的河槽隨著試驗(yàn)時(shí)間的推移呈現(xiàn)出不斷變寬的特征，最后階段河槽的形態(tài)與原始河型相比已經(jīng)發(fā)生巨大差異。這說明上游水沙條件的變化會(huì)對(duì)河床變形造成十分顯著的影響[9- 10]。

圖2 G1工況12h河床高程線圖

3.2 G2方案計(jì)算結(jié)果與分析

利用上節(jié)構(gòu)建的模型對(duì)G2方案下研究河段的河床縱向變化特征進(jìn)行計(jì)算，根據(jù)模型計(jì)算獲得的數(shù)據(jù)，獲得不同時(shí)段的河床高程線圖，其中12h和24h的河床高程線圖如圖3和圖4所示。由此可知，在G2工況條件下，研究河段的流量有所增大，河槽進(jìn)口部位水流的總能量也明顯增大，而從其對(duì)應(yīng)的河床等高線圖可以看出，河槽隨著試驗(yàn)時(shí)間的推移也呈現(xiàn)出不斷變寬的特征，本來十分蜿蜒的河道幾乎變成了順直型的河道[11]。究其原因，主要是近水口水流能量增大以后，河槽會(huì)通過變寬的方式，來消耗河槽能水流的富余能量，進(jìn)而實(shí)現(xiàn)能量的平衡。由于研究河段河床的材料抗沖刷能力較差，產(chǎn)生了比較嚴(yán)重的河槽沖刷。另一方面，隨著河槽的變寬，河槽內(nèi)水流的流速會(huì)慢慢變小，水流的挾沙能力也會(huì)不斷降低，水流中挾帶的泥沙會(huì)沉積在河槽之內(nèi)而造成淤積。隨著試驗(yàn)時(shí)間的推移，河床也會(huì)逐漸抬高，同時(shí)河槽兩側(cè)淤積抬高的速率并不相同，右岸的河床抬高速率明顯偏大。此外，隨著時(shí)間的推移，河槽內(nèi)出現(xiàn)了泥沙淤積形成的灘地，但是這些灘地星星點(diǎn)點(diǎn)，出現(xiàn)之后沒過多久就被沖刷掉了[12]。

圖3 G2工況12h河床高程線圖

圖4 G2工況24h河床高程線圖

3.3 G3方案計(jì)算結(jié)果與分析

利用上節(jié)構(gòu)建的模型對(duì)G3方案下研究河段的河床縱向變化特征進(jìn)行計(jì)算，根據(jù)模型計(jì)算獲得的數(shù)據(jù)，獲得不同時(shí)段的河床高程線圖，其中6h、12h、18h和24h的河床高程線圖如圖5～圖8所示。由此可知，在G3方案下，河槽的變化呈現(xiàn)出更為特殊的特點(diǎn)。剛開始，河槽的地形變化特征與G2方案比較類似，隨著試驗(yàn)河槽進(jìn)水口流量的增加，總能量迅速增大，造成下游河槽沖刷變寬，同時(shí)河床逐漸淤積抬高。隨著時(shí)間的逐步推移，河槽兩岸產(chǎn)生了比較嚴(yán)重的泥沙淤積，進(jìn)而形成了新的河灘，受到河灘的約束，河道的河槽又開始變窄，進(jìn)而造成河槽內(nèi)水流能量的不平衡。之后，河槽內(nèi)由于泥沙的淤積作用形成若干面積大小不等的灘地，其中最大的一塊處于河道中央，將原來的河槽一分為二，形成了兩個(gè)繞流的岔道，其余的小塊灘地則將河槽分為幾條寬度不一的岔道。

圖5 G3工況6h河床高程線圖

圖7 G3工況18h河床高程線圖

圖8 G3工況24h河床高程線圖

河槽內(nèi)的地形特點(diǎn)之所以會(huì)發(fā)生上述變化，主要原因是河槽內(nèi)水流的能量過大，河槽要達(dá)到能量的平衡狀態(tài)，就需要通過某種方式降低河槽內(nèi)水流的能量。在開始時(shí)，主要是通過對(duì)兩側(cè)河岸的沖刷作用使河槽變寬，進(jìn)而增大水流的阻力，但是這也會(huì)造成水流速度降低，挾沙能力不足。因此，隨著過流時(shí)間的增長，河槽內(nèi)會(huì)產(chǎn)生大量的泥沙淤積，并在河槽中央形成大小不同的灘地，在灘地的兩側(cè)則形成了兩個(gè)或多個(gè)流速的集中區(qū)域，使水流的沖刷作用繼續(xù)增強(qiáng)。所以，隨著時(shí)間的推移，河槽內(nèi)形成的灘地亦被沖刷，同時(shí)在河槽內(nèi)沉積形成新的灘地，并把河槽進(jìn)一步劃分為幾個(gè)岔道，進(jìn)而造成了十分顯著的游蕩河型特點(diǎn)。

3.4 G4方案計(jì)算結(jié)果與分析

利用上節(jié)構(gòu)建的模型對(duì)G4方案下研究河段的河床縱向變化特征進(jìn)行計(jì)算，根據(jù)模型計(jì)算獲得的數(shù)據(jù)，獲得不同時(shí)段的河床高程線圖，其中12h河床高程線圖如圖9所示。由此可知，開始河槽的變化狀況與前幾個(gè)工況比較類似。隨著進(jìn)水口能量的不斷增大，河槽持續(xù)被沖刷、變寬，最后形成了比較典型的順直型河道。隨著試驗(yàn)時(shí)間的進(jìn)一步推移，兩側(cè)的河床受泥沙沉積的作用不斷抬高并出現(xiàn)河灘，最終造成試驗(yàn)河段的河槽逐漸變窄、變深。

圖9 G4工況12h河床高程線圖

4 結(jié)論

文章通過數(shù)值模擬的方法，選擇典型的山區(qū)高含沙河段，對(duì)不同水沙條件下的河道地形變化進(jìn)行了研究。結(jié)果顯示，蜿蜒河槽隨著流量的增加，斷面河寬也在逐漸變大，當(dāng)流量增加到一定程度時(shí)，斷面河寬變化得比較復(fù)雜，呈現(xiàn)出時(shí)而變寬、時(shí)而變窄的特征。在河槽的中心部位，河床也表現(xiàn)為時(shí)而沖刷、時(shí)而淤積的特征。如果流量進(jìn)一步增大，河槽持續(xù)被沖刷、變寬，最后形成了比較典型的順直型河道。文章的研究結(jié)論對(duì)河流上游水庫調(diào)整下泄水量，通過改變下游河道的水沙條件保證下游河道形態(tài)穩(wěn)定性具有重要意義。當(dāng)然，本次研究河段的河床存在比較顯著的易沖刷性，后續(xù)研究中需要針對(duì)各種不同性質(zhì)河道進(jìn)行研究，以便獲得更為廣泛的結(jié)論。