宋加平，歐陽(yáng)高明，李 燦

（廣東中灝勘察設(shè)計(jì)咨詢有限公司,廣東 肇慶 526060）

1 引言

水動(dòng)力學(xué)特征反映了河道水流滲流現(xiàn)狀[1-2],與水利設(shè)施調(diào)控、運(yùn)營(yíng)密切相關(guān),研究河道水動(dòng)力學(xué)特征,有助于水工建筑設(shè)計(jì)與運(yùn)營(yíng)。水沙演變特征反映了河道泥沙遷移、懸浮特性,對(duì)河道清淤、排沙及沖砂設(shè)計(jì)具有參考價(jià)值,也有利于流域水土保持設(shè)計(jì)[3,4]。蔡玉鵬等[5]、王蕊[6]為研究河道水動(dòng)力學(xué)特征,采用模型試驗(yàn)方法,開展了河道原型比尺復(fù)制設(shè)計(jì),分析了模型試驗(yàn)過(guò)程中河床水位、流速及泥沙懸浮特征,豐富了工程設(shè)計(jì)參考成果。華中等[7]、程正飛等[8]、張瑜等[9]采用了流場(chǎng)離散元計(jì)算方法,基于Fluent 等平臺(tái)建立河道概化模型,探討河道在不同方案下流速、含沙量及水位變化特征,反映了河道運(yùn)營(yíng)現(xiàn)狀。針對(duì)水沙特征,馮淑琳等[10]、陳志軍[11]、劉富強(qiáng)等[12]采用MIKE 等水沙模擬平臺(tái)開展了泥沙淤積分析,獲得了河床泥沙淤積厚度演變及河道斷面含沙量變化特征,為河道清淤、沖排沙提供了參照。本文基于北江社崗流域河床特征,利用MIKE21 開展降雨、植被覆蓋下河道泥沙淤積與水沙演變特征分析。

2 設(shè)計(jì)模擬分析

2.1 工程概況

作為粵北地區(qū)重要地表水源,北江乃是廣東省境內(nèi)重要水系流域。北江流域內(nèi)包括有韶關(guān)、清遠(yuǎn)等地區(qū),全長(zhǎng)超過(guò)460 km,總流域面積達(dá)4.7×104km2,全流域內(nèi)建設(shè)水利設(shè)施包括有北江飛來(lái)峽水庫(kù)、北江清遠(yuǎn)水利樞紐工程及其他中小型水利工程,北江各水利工程承擔(dān)著流域內(nèi)水力發(fā)電、防洪排澇、通航、灌溉等水利功能,特別是流徑清遠(yuǎn)等地區(qū)的防洪聯(lián)圩工程,對(duì)北江中下游城市防洪帶來(lái)較大保障,堤頂最高處位于清遠(yuǎn)石角鎮(zhèn),距離河口面軸長(zhǎng)為350 m,全聯(lián)圩工程涉及流域面積超過(guò)3.7 萬(wàn)km2,聯(lián)圩所在區(qū)段內(nèi)降雨也是北江降水最多地區(qū),水位波幅段最大也是該區(qū)域,最高水位為38.5 m。與北江聯(lián)圩工程密切相關(guān)的水利設(shè)施為飛來(lái)峽水庫(kù),該水庫(kù)年降水量超過(guò)1800 mm,設(shè)計(jì)最大庫(kù)容量為3800 萬(wàn)m3,主壩最大高程達(dá)65.5 m,全壩身均設(shè)置有高度為2.5 m 的防浪墻,采用全階梯式溢流面設(shè)計(jì)形式的溢洪道,配置有可承擔(dān)最大泄流量300 m3/s 的泄洪閘,汛期可同時(shí)四孔泄流,閘室凈寬度達(dá)8.6 m,但由于受上游降水及區(qū)段內(nèi)泥沙懸浮影響,導(dǎo)致泄洪閘泥沙淤積過(guò)高,不利于通流。飛來(lái)峽水庫(kù)其一方面自身承擔(dān)著農(nóng)業(yè)灌溉、水資源調(diào)度等作用,同樣可對(duì)下游聯(lián)圩工程水位進(jìn)行控制,減弱水庫(kù)泄流對(duì)下游防洪設(shè)施影響。飛來(lái)峽水庫(kù)流域內(nèi)總?cè)丝诮咏?0 萬(wàn),轄區(qū)內(nèi)共有10 個(gè)鄉(xiāng)鎮(zhèn)。人口活動(dòng)對(duì)北江干流飛來(lái)峽水庫(kù)區(qū)段兩側(cè)生態(tài)植被影響帶來(lái)一定影響,特別是由于農(nóng)業(yè)灌渠的建設(shè),部分河道植被覆蓋密度受影響較大,導(dǎo)致河坡水土流失嚴(yán)重,根據(jù)下游聯(lián)圩工程監(jiān)測(cè)表明,飛來(lái)峽水庫(kù)社崗流域內(nèi)在汛期含沙量最大接近6.5 kg/m3,水質(zhì)狀態(tài)受影響較大,同時(shí)泥沙懸浮對(duì)下游聯(lián)圩工程的沖擊作用,導(dǎo)致堤身受沖刷面區(qū)域達(dá)到72%。綜合北江干流與飛來(lái)峽水庫(kù)現(xiàn)狀分析,北江流域內(nèi)水利設(shè)施一方面受降雨徑流活動(dòng)節(jié)制影響,同時(shí)也受到生態(tài)植被覆蓋現(xiàn)狀影響。為此,水文部門考慮以北江干流飛來(lái)峽水庫(kù)K3+005～K6～605所在區(qū)段的社崗流域?yàn)榉治鰧?duì)象,開展該區(qū)段內(nèi)降水、生態(tài)植被覆蓋對(duì)水沙特征影響分析。

2.2 模擬設(shè)計(jì)

本文采用MIKE21 對(duì)飛來(lái)峽水庫(kù)社崗流域水系河道進(jìn)行概化[11-13],概化圖中對(duì)部分水面寬度低于1m 的水渠進(jìn)行忽略,主要水系包括了飛來(lái)峽水庫(kù)上下游支流,而本文模擬的社崗流域區(qū)段也位于該區(qū)段內(nèi)。根據(jù)流域內(nèi)實(shí)際河道寬度、水位高程等信息輸入至河道斷面屬性特征中,區(qū)段河道起始點(diǎn)空間X～Z 數(shù)據(jù)均采用社崗流域內(nèi)河道實(shí)際數(shù)據(jù)。

降水尺度按照2010 年～2020 年10 年共185 場(chǎng)次降水強(qiáng)度數(shù)據(jù),按照尺度范圍內(nèi)降水?dāng)?shù)據(jù)分別降雨強(qiáng)度150 mm/d、300 mm/d、450 mm/d、600 mm/d、750 mm/d、900 mm/d,在MIKE21 中按照降雨與徑流的線性相關(guān)數(shù)據(jù),設(shè)定區(qū)段內(nèi)水系凈流參數(shù),降水與徑流關(guān)系見圖1,按照汛期與年平均降水量進(jìn)行設(shè)定,社崗流域內(nèi)最大徑流量為1575 萬(wàn)m3。研究區(qū)段內(nèi)生態(tài)植被覆蓋以環(huán)境科學(xué)中心監(jiān)測(cè)NDVI 數(shù)據(jù)進(jìn)行計(jì)算[14],并按照?qǐng)D2 中近十年NDVI 植被覆蓋數(shù)據(jù)分別設(shè)定有NVDI值0.42、0.48、0.6、0.66、0.72,MIKE21 中模擬計(jì)算時(shí)也相應(yīng)輸入NVDI 植被覆蓋特征。

圖1 降水與徑流關(guān)系圖

圖2 近十年植被覆蓋NDVI 值

基于MIKE 流域模型的建立,其邊界條件考慮水動(dòng)力演化與水沙擴(kuò)散特征,外部接觸面邊界直接采用飛來(lái)峽水庫(kù)泄流期最高水位55.5 m 作為計(jì)算工況,內(nèi)部水質(zhì)泥沙懸浮遷移等以能量等替交換形式,河段糙率設(shè)定為0.04,河道初始含沙量設(shè)定為飛來(lái)峽水庫(kù)入庫(kù)平均含沙量,為2.5 kg/m3。在流域內(nèi)泄流期工況,開展降水、植被覆蓋特征因素對(duì)研究區(qū)段內(nèi)水動(dòng)力學(xué)特征、水沙演化特征影響分析。

3 流域泥沙淤積特征影響

3.1 降水影響

根據(jù)不同降雨強(qiáng)度方案下流域水動(dòng)力學(xué)特征計(jì)算,獲得社崗流域區(qū)段內(nèi)泥沙淤積厚度變化特征,見圖3,為分析方便,以飛來(lái)峽水庫(kù)K3+005 樁號(hào)點(diǎn)為河道橫斷面起點(diǎn)（0 m）。分析泥沙淤積特征可知,泥沙淤積厚度與降雨強(qiáng)度為正相關(guān)關(guān)系,在運(yùn)營(yíng)期為100 d 下,降雨強(qiáng)度150 mm/d 下流域斷面680 m 處泥沙淤積厚度為0.35 m,而降雨強(qiáng)度為300 mm/d、600 mm/d、900 mm/d 下同斷面處泥沙淤積厚度較之分別增長(zhǎng)了26.4%、47.2%、71.7%；從社崗流域河道斷面整體泥沙淤積厚度來(lái)看,降雨強(qiáng)度150 mm/d 下泥沙淤積平均厚度為0.35 m,最大、最小厚度分別為0.42 m、0.3 m,分別位于斷面1800 m、2550 m,而在降水方案中降雨強(qiáng)度每增長(zhǎng)150 mm/d 下,河道斷面上泥沙淤積平均厚度可增長(zhǎng)11.3%,最大、最小厚度在降雨強(qiáng)度600 mm/d 前基本一致,均處于同一斷面,而降雨強(qiáng)度750 mm/d、900 mm/d 下最大、最小斷面均發(fā)生改變,分別位于1090 m、3250 m。筆者認(rèn)為,河道泥沙淤積厚度最大值斷面上水動(dòng)力學(xué)活動(dòng)較緩,不利于河道沖砂、排流活動(dòng),而淤積厚度最小的斷面上動(dòng)水沖擊作用較強(qiáng),泥沙受水流裹挾遷移力較大；當(dāng)降雨強(qiáng)度過(guò)大,河道內(nèi)水動(dòng)力學(xué)活動(dòng)的強(qiáng)弱斷面出現(xiàn)改變[15]。從降雨強(qiáng)度對(duì)泥沙淤積厚度影響變幅來(lái)看,降雨強(qiáng)度過(guò)大,泥沙淤積厚度的增幅實(shí)質(zhì)上在減弱,如運(yùn)營(yíng)期100 d 下,降雨強(qiáng)度600 mm/d 與150 mm/d、300 mm/d 的斷面平均淤積厚度分別具有增幅45.7%、15.9%,而降雨強(qiáng)度750 mm/d、900 mm/d 較之600 mm/d 下分別僅有增幅7.8%、14.5%。由此說(shuō)明,降雨強(qiáng)度增大,河道泥沙淤積趨勢(shì)在減弱,高強(qiáng)度降雨下的徑流活動(dòng)一方面帶來(lái)水土流失,另一方面河道徑流量也得到增長(zhǎng),河道泥沙遷移作用增強(qiáng),有助于控制泥沙的沉降。

圖3 降雨對(duì)河道斷面泥沙淤積厚度特征

當(dāng)運(yùn)營(yíng)期增長(zhǎng)10 倍,為1000 d 時(shí),泥沙淤積受時(shí)序效應(yīng)影響,斷面厚度均得到增大,在降雨強(qiáng)度300 mm/d 時(shí),其全斷面上泥沙淤積厚度較之運(yùn)營(yíng)100 d 時(shí),增幅分布為1.7～2.2倍,淤積厚度平均值的增幅為1.82 倍,淤積增長(zhǎng)最為顯著的屬斷面2550 m,該斷面也是水動(dòng)力學(xué)活動(dòng)處于較活躍。降雨強(qiáng)度增大,運(yùn)營(yíng)期1000 d 下泥沙淤積厚度增長(zhǎng)最顯著的為降雨強(qiáng)度600 mm/d 下,泥沙淤積平均厚度的增幅達(dá)2.05 倍,而各斷面上淤積厚度增幅分布為1.9～2.22 倍,而降雨強(qiáng)度150 mm/d及高降雨強(qiáng)度750 mm/d、900 mm/d 等方案中,其平均淤積厚度值的增幅分別為1.9 倍、1.85 倍、1.76 倍。分析表明,在過(guò)大或過(guò)低的降雨強(qiáng)度方案下,泥沙淤積厚度隨運(yùn)營(yíng)期變化均不是最敏感,而以中等降雨強(qiáng)度方案中泥沙淤積厚度隨運(yùn)營(yíng)期變化敏感度最大。

3.2 植被影響

根據(jù)對(duì)不同NDVI 植被覆蓋方案下泥沙淤積分析,獲得了植被覆蓋對(duì)流域泥沙淤積影響,見圖4。從圖中可知,植被覆蓋參數(shù)NDVI 值愈大,則泥沙淤積厚度愈小,即NDVI 參數(shù)值與泥沙淤積厚度具有負(fù)相關(guān)關(guān)系,運(yùn)營(yíng)期為100 d 時(shí),在NDVI 參數(shù)0.42 方案中,斷面380m 處泥沙淤積厚度為0.63 m,而植被覆蓋方案的NDVI 值每增大0.6,則該斷面處泥沙淤積平均可減少12.8%。NDVI 值0.42 方案中,河道全斷面上泥沙淤積厚度平均值為0.65 m,其隨每階次NDVI0.06 的變化,厚度值的平均降幅為12.5%,且各方案間平均厚度降幅較穩(wěn)定,基本均接近平均降幅,表明植被覆蓋不僅有利于固土護(hù)沙[16],且植被覆蓋密度的變化對(duì)水土流失效應(yīng)的影響變化較恒定。

圖4 植被對(duì)河道斷面泥沙淤積厚度特征

當(dāng)運(yùn)營(yíng)期增長(zhǎng)至1000 d 時(shí),各方案中泥沙淤積厚度值均有增長(zhǎng),在NDVI 參數(shù)0.48 方案中,其流域斷面上泥沙厚度的增幅分布為54.2%～57.9%,淤積厚度平均值的增幅為55.2%,河道斷面間泥沙淤積的變幅最大、最小均與運(yùn)營(yíng)期100 d 下一致,而其他NDVI 方案中,同樣斷面間淤積厚度的最大、最小變幅區(qū)間斷面均與運(yùn)營(yíng)期100 d 一致,未發(fā)生改變。分析泥沙淤積厚度的增幅斷面,各NDVI 參數(shù)方案中泥沙淤積厚度隨運(yùn)營(yíng)期的增長(zhǎng)最大、最低均一致,分別為1500 m～1800 m、1090 m～1500 m,即在該區(qū)間斷面上,水動(dòng)力學(xué)特征參數(shù)受運(yùn)營(yíng)時(shí)序效應(yīng)影響最為敏感,且植被覆蓋NDVI 參數(shù)的變化,不改變河道時(shí)序敏感區(qū)間所在斷面。

4 流域水沙演變影響特征

為研究流域水沙演變特征,對(duì)不同降雨強(qiáng)度、植被覆蓋方案下開展河道斷面含沙量參數(shù)計(jì)算,獲得運(yùn)營(yíng)期1000 d 時(shí)各方案中斷面含沙量特征,見圖5。

圖5 河道斷面含沙量影響變化

從圖中可知,降雨強(qiáng)度愈大,則含沙量愈高,但高強(qiáng)度強(qiáng)度方案下斷面含沙量較大的區(qū)間延伸至研究區(qū)段的中游,降雨強(qiáng)度150 mm/d 下斷面1090 m 處含沙量為1.84 kg/m3,而降雨強(qiáng)度每增長(zhǎng)150 mm/d,該斷面處含沙量平均增長(zhǎng)16.5%,降雨150 mm/d 方案下河道斷面平均含沙量為2.16 kg/m3,隨降雨強(qiáng)度增長(zhǎng),斷面平均含沙量的增幅13.2%。在降雨強(qiáng)度低于600 mm/d 時(shí),各斷面上含沙量從上游至下游均為遞增,但降雨強(qiáng)度在750 mm/d 時(shí)其斷面2550 m 后含沙量增幅較小,甚至在強(qiáng)度900 mm/d 方案下斷面2550 m 后含沙量出現(xiàn)稍有降低的現(xiàn)象。分析認(rèn)為,降雨強(qiáng)度與含沙量具有正相關(guān),但降雨強(qiáng)度過(guò)大,會(huì)影響河道含沙量在斷面上分布特征,改變上游至下游泥沙遷移、沉降的規(guī)律。

當(dāng)植被NDVI 參數(shù)值愈大,則河道斷面含沙量愈低,在NDVI 值0.42 時(shí),其斷面平均含沙量為4.17 kg/m3,而植被覆蓋參數(shù)NDVI 值每梯次增長(zhǎng)0.06,則其斷面平均含沙量可減少16.3%,表明植被覆蓋愈廣泛,有利于限制泥沙流失。在各NDVI 值參數(shù)方案中,同一斷面上含沙量的降幅在各方案中基本保持一致,如斷面1500 m 處含沙量在NDVI 值0.54～0.6 間、0.66～0.72 間的降幅均穩(wěn)定在25%,表明植被覆蓋參數(shù)對(duì)河道斷面含沙量影響具有恒定性。從水土保持及河道泥沙治理方面考量,在降雨強(qiáng)度處于中等時(shí),更應(yīng)重點(diǎn)考慮水土流失,而流域內(nèi)植被覆蓋應(yīng)盡量廣泛,抑制泥沙流失。

5 結(jié)論

本文主要得到以下三點(diǎn)結(jié)論：

（1）泥沙淤積厚度與降雨強(qiáng)度為正相關(guān)關(guān)系,降雨強(qiáng)度每增長(zhǎng)150 mm/d 下,運(yùn)營(yíng)期100 d 時(shí)河道斷面淤積平均厚度可增長(zhǎng)11.3%,降雨強(qiáng)度較大時(shí),會(huì)改變淤積厚度峰、谷值所處斷面,同時(shí)淤積厚度增幅也減?。荒嗌秤俜e具有運(yùn)營(yíng)時(shí)序增長(zhǎng)效應(yīng),且以中等強(qiáng)度降雨方案下泥沙淤積受此效應(yīng)更為敏感。

（2）植被NDVI 值與泥沙淤積厚度具有負(fù)相關(guān)關(guān)系,每階次NDVI0.06 的變化,運(yùn)營(yíng)期100 d 時(shí)淤積厚度平均降幅為12.5%,各NDVI 值方案間降幅穩(wěn)定；運(yùn)營(yíng)期100 d 與1000 d 下,淤積厚度變幅最大、最小區(qū)間斷面均保持一致,為1500 m～1800 m、1090 m～1500 m,該區(qū)間斷面水動(dòng)力特征受運(yùn)營(yíng)時(shí)序效應(yīng)最為敏感。

（3）降雨強(qiáng)度愈大,含沙量愈高,降雨150 mm/d 的階次增長(zhǎng),可引起斷面平均含沙量的梯次增幅為13.2%,降雨強(qiáng)度過(guò)大,上下游斷面泥沙遷移、沉降效應(yīng)受改變；植被NDVI值每梯次增長(zhǎng)0.06,則斷面平均含沙量可減少16.3%,且同一斷面含沙量的降幅保持恒定。