謝　瑞，姬昌輝，王永平，葛　慧

(1:南京水利科學研究院，水文水資源與水利工程科學國家重點實驗室，南京 210024)(2:南京信息工程大學，水文氣象學院，南京 210044)

?

波浪作用下湖泊底泥的輸沙試驗研究*

謝瑞1，姬昌輝1，王永平1，葛慧2

(1:南京水利科學研究院，水文水資源與水利工程科學國家重點實驗室，南京 210024)(2:南京信息工程大學，水文氣象學院，南京 210044)

湖泊底泥的運動過程產(chǎn)生內(nèi)源污染，加劇湖泊生態(tài)環(huán)境的惡化. 研究湖泊底泥在波浪作用下的輸沙規(guī)律，可為研究湖泊的水質(zhì)變化成因及生態(tài)環(huán)境治理提供參考依據(jù).利用大型波浪水槽，在波浪作用條件下對太湖、龍感湖、巢湖的底泥進行了起動和輸沙試驗，對多組波浪水文條件下湖泊底泥的輸沙試驗結(jié)果進行分析，詳細地闡述了太湖、龍感湖、巢湖底泥在波浪作用下的輸沙變化規(guī)律. 分析整理試驗數(shù)據(jù)，得出輸沙率變化公式，為3個湖泊的水質(zhì)變化成因分析及生態(tài)環(huán)境治理提供參考依據(jù).

太湖；龍感湖；巢湖；底泥；波浪；輸沙率

湖泊底泥在波浪作用下的運動過程產(chǎn)生了內(nèi)源污染，是湖泊水體中營養(yǎng)鹽來源的重要方面，從而加劇了湖泊的富營養(yǎng)化. 前人對水流泥沙運動規(guī)律進行了大量的試驗和分析，有著豐富的研究成果. 1960s Jonsson就對波浪作用下的泥沙起動進行了大量的研究，建立了邊界層等理論[1]，竇國仁院士對于波浪作用下泥沙的起動規(guī)律進行了深入的研究，得出了泥沙起動公式[2]. 黃建維對新港浮泥作了專題研究[3]，嚴愷等在港口波浪潮流泥沙運動規(guī)律的研究中形成了一整套完整的理論[4]. 王軍等對淤積固結(jié)條件下粘性細泥沙起動沖刷研究進行了綜述[5]，在沖刷水槽試驗中環(huán)形水槽各種侵蝕沖刷水槽均被廣泛使用. 曹叔尤等利用環(huán)形水槽試驗建立了淤積物的沖刷率公式[6]等等. 研究表明，波浪作用下的細顆粒泥沙與粗顆粒泥沙的起動機理存在著較大的差異，粗顆粒泥沙的起動以形成沙波運動為主要特征；淤泥質(zhì)泥沙的起動以泥沙的卷起、懸揚為主要特征；細顆粒泥沙的起動則介于粗顆粒與淤泥質(zhì)泥沙之間，以底泥的卷起、懸揚為主，同時又有部分底質(zhì)的輸移，其運動特征與泥沙顆粒粒徑和顆粒級配有密切關(guān)系. 因此一些學者對細顆粒泥沙的卷起特性作了一些研究，由于細顆粒泥沙相對于粗顆粒泥沙的表面面積較大，顆粒間的相互摩擦、吸附將形成一個框架結(jié)構(gòu)，往往比粗顆粒泥沙有更多的空隙，其達到密實穩(wěn)定的歷時相對較長，重度變化也較大，故對起動流速的影響也頗大. 因此，在研究細顆粒泥沙的起動時，除應考慮波浪要素對起動流速的影響外，還應考慮細顆粒泥沙的特性及其在密實過程中重度對起動流速的影響[7]. 國內(nèi)外學者對于粗顆粒泥沙在波浪作用下的起動、輸移的研究較多，也較為成熟. 而對于淤泥質(zhì)泥沙的研究則相對薄弱.

同時，對于泥沙運動規(guī)律的研究主要集中在對于海岸[8]、河口[9]、河流[10]、閘下[11]泥沙起動及沉降方面，對于波浪作用下泥沙起動后的輸沙研究主要局限于利用數(shù)值模擬計算的方法，對于河床的輸沙研究也是局限于有建筑物作用下的河床、海床的沖刷變化，對于湖泊底泥輸沙研究則相對缺乏. 而在河口湖泊地區(qū)，當淤積粘性細泥沙吸附部分污染物時，這些污染物的遷徙擴散又與粘性細顆粒泥沙的起動輸沙過程緊密聯(lián)系，因此研究波浪作用下粘性細顆粒泥沙起動輸沙問題具有重要實際意義[12]. 太湖位于江蘇省南部、浙江省北部，是中國第二大淡水湖；龍感湖位于湖北省和安徽省交界處，也是一個淡水湖泊；巢湖位于江淮之間，地處合肥、巢湖二市和肥東、肥西、廬江三縣境內(nèi)，是中國著名的五大淡水湖之一. 雖然3個湖泊從地理位置來看相距較遠，但底泥的泥沙特性相似，均屬于細顆粒泥沙范疇[13]. 本文通過波浪水槽試驗，研究3個湖泊底泥起動以后的輸沙規(guī)律，可以預測一定風浪條件下湖泊底泥的輸沙量，使我們認識到在大風條件下湖泊水體含沙濃度的時空變化，可以進一步預測底泥營養(yǎng)鹽的釋放，為研究湖泊水質(zhì)波動成因分析以及風浪作用下水環(huán)境變化的預測分析提供參考依據(jù).

1　底泥特性

太湖的底泥樣品采自于太湖的梅梁灣，龍感湖的底泥樣品采自于龍感湖的湖心部位，巢湖的底泥樣品采自于兩個位置：巢湖(J)靠近湖岸；巢湖(Y)靠近湖中心位置. 樣品采集時間是2004年夏季，當時太湖梅梁灣水質(zhì)屬于重度富營養(yǎng)狀態(tài)[14]；巢湖西半湖水質(zhì)屬于重度富營養(yǎng)狀態(tài)，東半湖水質(zhì)相對好于西半湖[15]；龍感湖自1974年以來，流域大范圍化學肥料的使用增加了入湖可溶性營養(yǎng)鹽的排放量，從而導致水質(zhì)逐漸惡化，總體而言水質(zhì)好于太湖和巢湖[16]. 3個湖泊均屬于大型淺水湖泊，對于研究內(nèi)源污染問題具有共性. 樣品采集使用彼得森抓泥器，采集湖底表層20 cm以內(nèi)的底泥.

表1 淤泥粒徑

從淤泥外觀顏色來看，太湖底泥呈暗灰色，龍感湖和巢湖底泥呈土黃色. 3個湖泊內(nèi)底泥樣品的粒度分析結(jié)果見表1，根據(jù)中值粒徑大小，3個湖泊的底泥樣品均屬于粘性細顆粒泥沙范疇.

2　試驗和計算方法

試驗裝置如圖1，試驗水槽全長25 m，高0.8m，寬0.5m，水槽進口設有平推式造波設備，以產(chǎn)生規(guī)則波，水槽尾端設有消浪網(wǎng)，中間的鋪沙段裝有間距為5 m的兩個采集波高的波高儀. 整個系統(tǒng)采用計算機自動控制. 鋪沙段長度為6 m，鋪設泥層厚度為10 cm. 試驗前底泥在水槽中經(jīng)攪拌均勻后自然沉降，2天后測定底泥容重，然后開始試驗. 在天然湖泊中，底泥容重的垂直分布規(guī)律是隨著深度增加容重逐漸增大. 鋪沙段底泥樣品經(jīng)自然沉降2 d后表層底泥容重達到1.30 g/cm3左右，接近天然湖泊底泥的表層密實狀態(tài)，因此可以用來研究湖泊表層底泥的起動和沖刷. 試驗時波高先從小波開始，然后逐漸增大，直至底泥處于起動狀態(tài)下測定其波高. 在淤泥起動試驗中，觀察和測定了不同水深不同波周期下淤泥起動的波高Hc. 在此試驗的基礎上，重點研究3個湖泊4種泥樣波浪作用下的輸沙，在一種泥樣一種容重下，試驗水深為30 cm，波周期為1.0 s試驗條件下，共進行了12組試驗. 考慮實測太湖在大風天8 m/s 的風速持續(xù)1 h 以后的水體含沙濃度最大可以達到2.58 kg/m3[14]，在輸沙試驗開始時通過測定垂線含沙濃度，1 h后垂線平均含沙濃度與天然湖泊相當，由此確定輸沙試驗的最大的波高H5，在起動波高與最大波高之間劃分4個不同波高等級H1～H4，并使H1>H2>H3>H4>H5>Hc，從而研究輸沙率與波高的變化規(guī)律.

圖1 波浪水槽示意圖(單位：m)Fig.1　Layout of wave flume

通過測定水體的含沙濃度來計算不同波高下的輸沙率. 以太湖底泥為研究重點，研究5組不同波高下的輸沙變化過程，對于龍感湖和巢湖底泥分別實測2組波高下的輸沙變化過程，輸沙試驗的波浪水文條件見表2，試驗結(jié)果按如下公式計算：

(1)

(2)

式中，E為輸沙率，ρ為水體平均含沙濃度，V為水體總體積，s為淤泥床面面積，t為掀沙時間，ρ1、ρ2、ρ3為垂線含沙濃度分布，V1、V2、V3為分層水體體積.

表2　試驗波浪水文條件*

*試驗裝置水深30 cm,波周期為1.0 s,淤泥自然，沉降密實.

3　結(jié)果與分析

波浪水槽試驗先從小波開始，隨著波高的逐漸加大，在臨底水流切應力的作用下，床面上的一些輕質(zhì)絮狀物首先開始前后滾動，隨著波浪動力的加大，床面泥沙由少量動到普遍運動，但泥沙顆粒做往復振動，基本不懸浮，此時判定為泥沙起動[1]. 對粘性細顆粒泥沙，粘結(jié)力和靜水壓力遠超過重力，起動后再加大波浪動力，因突然失去粘結(jié)力和附加靜水壓力而懸浮，床面泥沙開始懸起，瞬時含沙量隨波浪底部水質(zhì)點速度的變化而變化，短時間內(nèi)底層含沙量明顯高于上層，底部形成高含沙層. 當波浪作用時間延長后，泥沙開始向上層水體擴散. 當波浪作用停止后，懸起的泥沙迅速沉降，高含沙層隨即消失. 上述現(xiàn)象表明，波浪作用下底泥的輸移不同于較粗泥沙的底部推移質(zhì)輸沙，也不同于水流作用下淤泥的全斷面懸浮輸移，其輸移主要以臨底高含沙層泥沙沿波浪底部水質(zhì)點運動方向往復運動. 試驗中發(fā)現(xiàn)，開始時，高含沙層水體的厚度和平均含沙量均隨波浪動力的增大而增大，但當厚度達到一定值后，即使再加大波浪動力，其厚度也基本不變，而此時，高含沙層平均含沙量迅速增大. 由于波浪非對稱性引起的傳質(zhì)速度較小，盡管波浪興起的含沙量較大，但在一個周期內(nèi)，其凈輸沙也非常小. 同時，在波高較小的條件下，泥沙的懸揚和擴散往往達不到水體上部，而只限于臨底一定高度內(nèi)，泥沙擴散高度隨波高的增大而增大，3個湖泊4種泥樣在波浪作用下的輸沙過程規(guī)律基本相似.

泥沙起動后，加大波浪動力，泥沙自床面懸浮，向上層水體擴散. 含沙量沿水深分布自泥沙床面向上逐漸減小. 通過實測垂線含沙量分布，根據(jù)公式計算波浪輸沙率. 太湖底泥的輸沙率變化過程為：當波高小于12 cm時，波浪的輸沙率比較?。划敳ǜ叽笥?2 cm時，波浪的輸沙率明顯增大(圖2a). 試驗中最大的波高為13.59 cm，當波浪達到最大波高時，起始輸沙率比較大，沖刷5 min后輸沙率下降明顯，10 min后輸沙率呈緩慢遞減態(tài)勢. 龍感湖底泥的輸沙率變化過程為：當波高小于10 cm時，波浪的輸沙率比較?。划敳ǜ叽笥?0 cm時，波浪的輸沙率明顯增大(圖2b). 試驗中最大的波高為11.90 cm時，當波浪達到最大波高時，起始輸沙率比較大，輸沙試驗在30 min以內(nèi)輸沙率下降明顯，在90 min的試驗過程中輸沙率呈緩慢遞減態(tài)勢. 巢湖底泥的輸沙率變化過程為：當波高小于10 cm時，波浪的輸沙率比較小；當波高大于10 cm時，波浪的輸沙率明顯增大(圖2c、2d). 試驗中最大的波高為13.11 cm，當波浪達到最大波高時，起始輸沙率比較大，輸沙試驗在20 min以內(nèi)輸沙率下降明顯，在60 min以后輸沙率基本趨于穩(wěn)定.

圖2 太湖(a)、龍感湖(b)、巢湖(J)(c)和巢湖(Y)(d)的底泥輸沙率變化過程Fig.2　Scour rates of sediments from Lakes Taihu (a), Longgan (b), Chaohu (J) (c) and Chaohu (Y) (d)

從圖2可以看出，波高越大輸沙率越大，同樣波高下，輸沙率隨著時間呈遞減的趨勢. 波浪作用下，泥沙自床面起始懸浮強度大，隨著波浪作用時間的延長，水體含沙量逐漸增大，水流夾帶泥沙的能力也有所降低，輸沙強度逐漸減小，波浪沖刷60 min以后，輸沙率基本近似穩(wěn)定.

4　公式的建立

前人波浪泥沙研究結(jié)果表明[17-18]，Shields參數(shù)是研究波浪泥沙輸移的良好指標，波浪的挾沙能力應與其底部波浪Shields參數(shù)呈正比. 按Madsen等的方法[19]，波浪Shields參數(shù)可寫為：

(3)

式中，ρs和ρ分別為泥沙顆粒和水的比重：g為重力加速度；d為泥沙顆粒粒徑；τm為床面剪切應力的最大瞬時值，定義為：

(4)

式中，ub為波浪水質(zhì)點近底水平速度的最大值，由微幅波理論確定；fw為Jonsson波浪摩擦系數(shù)，可按下式計算，即：

(5)

當日αb/γ≤1．59時，取fw=0.30. 其中，γ為床面糙率，對于平底床面通常取泥沙中值粒徑值，由于波浪作用下湖底泥面的沙紋尺度非常小，在此次計算中均取泥沙中值粒徑值D50；ab為近底水質(zhì)點運動的有效振幅. 用上述公式對本次試驗資料進行計算，3個湖泊4種底泥波浪作用下的起動切應力起動切應力在0.27～0.32 N/m2之間，平均值為0.30 N/m2，4種泥樣差別不大(表3)，可以看出湖泊底泥的特性具有相同的特征和起動的相似性. 不同波高下淤泥沖刷時切應力計算結(jié)果如表4所示，可以看出，在相同的水深和波周期條件下，當?shù)啄嗥饎雍?，隨著波高的增大，水體中含沙量亦逐漸增大，通過計算結(jié)果也可以看出水流對底泥沖刷的切應力亦逐漸增大.

表3　波浪作用下淤泥起動切應力

表4　波浪作用下淤泥沖刷時的切應力

(6)

其中指數(shù)相關(guān)求得隨沉積物結(jié)構(gòu)而變的系數(shù)E0=0.0107、α=2.0676. 通過天然的水深和波要素可以計算沖刷時切應力，進而計算天然情況下的底泥輸沙率(表5).

圖3 輸沙率曲線Fig.3　Curve of scour rate

表5 波浪作用下底泥的輸沙率

5　結(jié)論

1)3個湖泊4個泥樣的中值粒徑基本在0.01 mm以下，屬于黏性細顆粒泥沙范疇. 在泥沙特征以及起動、輸沙的過程方面具有相似性.

2)在30 cm水深、波周期為1.0 s條件下，3個湖泊底泥的起動波高介于5.5～7.2cm之間，平均值為6.5cm；起動切應力介于0.27～0.32 N/m2之間，平均值為0.30 N/m2.

3)泥沙起動后，波高越大，泥沙懸浮越強，輸沙率越大，水體含沙量沿水深分布自泥沙床面向上逐漸減小.

4)同樣波高作用下，輸沙率隨著波浪作用時間的延長呈遞減的趨勢. 隨著水體含沙量逐漸增大，水流夾帶泥沙的能力也有所降低，其強度逐漸減小，然后輸沙率達到近似穩(wěn)定. 在已知水深和波要素的條件下可以計算波浪作用下的沖刷切應力，進而根據(jù)公式計算輸沙率，從而得到大風天時波浪對湖泊床面的侵蝕程度.

[1]Jonsson IG, Wave boundary layer and friction factors. In: Proc. 10th Int Conf Coastal Eng. ASCE, Tokyo, 1966: 127-148.

[2]竇國仁. 竇國仁論文集. 北京: 中國水利水電出版社, 2003.

[3]黃建維. 連云港泥沙靜水沉降試驗及粘性泥沙沉降特性的探討. 南京水利科學研究院報告, 1980.

[4]嚴愷, 梁其荀. 海洋工程. 北京: 海洋出版社, 2002.

[5]王軍, 談光鳴, 舒彩文. 淤積固結(jié)條件下粘性細泥沙起動沖刷研究綜述. 泥沙研究, 2008, (3): 75-80.

[6]曹叔尤， 杜國翰. 粘性土沖淤的試驗研究. 泥沙研究, 1986, (4): 73-82.

[7]羅肇森, 羅勇. 浮泥挾沙力和輸沙規(guī)律的研究和應用. 南京: 南京水利科學研究院報告, 1998.

[8]張瑋, 濮勛, 廖迎娣. 淤泥質(zhì)海岸近海風電塔基局部沖刷計算研究. 海岸工程, 2010, 28(2): 105-109.

[9]孫志林, 張翀超, 黃賽花等. 黏性非均勻沙的沖刷. 泥沙研究, 2011, (3): 44-48.

[10]王萬杰, 韓龍喜, 孫玉龍等. 京杭運河底泥起動和沖刷率試驗研究. 中國農(nóng)村水利水電, 2007, (9): 30-35.

[11]田小平, 王斌, 陳劍. 烏牛新閘下游沖刷模型試驗及經(jīng)驗公式適用性. 水利水電科技進展, 2017, 34(4): 75-78.

[12]秦伯強. 太湖水環(huán)境面臨的主要問題、研究動態(tài)與初步進展. 湖泊科學, 1998, 10(4): 1-9. DOI 10. 18307/1998. 0401.

[13]錢寧, 萬兆惠. 泥沙運動力學. 北京: 科學出版杜, 1983.

[14]陳潤, 王躍奎, 高怡等. 2004～2008年太湖水質(zhì)變化原因及治理對策. 水電能源科學, 2010, 28(11): 35-37.

[15]劉偉. 巢湖藍藻重度發(fā)生時的水質(zhì)特征分析. 中國環(huán)境監(jiān)測, 2007, 23(4): 101-102.

[16]陳詩越, 吳愛琴, 于興修. 長江中下游湖泊富營養(yǎng)化過程探析——以龍感湖為例. 臨沂師范學院學報, 2004, 26(6): 58-61.

[17]Sato S, Horikawa K. Sand ripple geometry and sand transport mechanism due to irregular oscillatory flow. Proc. 21st Coastal Eng. Conf., ASCE, 1988: 1748-1762.

[18]Zhou YR, Chen YP, Ma QN. Threshold of sediment movement in different wave boundary layers.ChinaOceanEngineering, 2001, 15(4): 50-520.

[19]Madsen OS, Grant WD. Quantitative description of sediment transport by waves. Proc. 15th Coastal Eng. Conf., ASCE, 1976: 1093-1112.

Transporting experimental study of the bottom silt-clay particles under the lake wave movement

XIE Rui1, JI Changhui1, WANG Yongping1& GE Hui2

(1:StateKeyLaboratoryofHydrology-WaterResourcesandHydraulicEngineering,NanjingHydraulicResearchInstitute,Nanjing210024,P.R.China)(2:CollegeofHydrometeorology,NanjingUniversityofInformationScience&Technology,Nanjing210044,P.R.China)

Inter-sourced pollution coming from lake sediment movement has aggravated deterioration of the ecological environment. Sediment transport rate under the lake wave movement provides reference for investigating causes of lake water quality change and ecological environment governance. Sediment transporting experiments from Lakes Taihu, Longgan and Chaohu under wave movements were performed in laboratory wave tanks. The results of silt-clay variations under multi-group wave conditions in the lake were described in detail and then analyzed for the three lakes. The formula of the sediment transporting rate was obtained by analyzing experimental data. These would give reference for causes of water quality change and ecological environment governance.

Lake Taihu; Lake Longgan; Lake Chaohu; sediment; wave; sediment transport rate

J.LakeSci.(湖泊科學)， 2016， 28(3)： 669-675

10.18307/2016.0324

?2016 byJournalofLakeSciences

*國家自然科學青年基金項目(41301550)和水利部公益性行業(yè)科研專項經(jīng)費項目(201301068)聯(lián)合資助.2015-04-28收稿；2015-09-17收修改稿. 謝瑞(1964～)，男，高級工程師；E-mail：rxie@nhri.cn.