鄒宏榮,王　燁,陸國(guó)賓,李瓊芳,張　弘,虞美秀,陳蕓蕓,閆方秀

(1. 河海大學(xué)水文水資源學(xué)院, 江蘇 南京　210098;2. 河海大學(xué)國(guó)際河流研究所,江蘇 南京　210098;3. 河海大學(xué)水利水電學(xué)院, 江蘇 南京　210098)

?

淮河上游產(chǎn)沙模擬研究

鄒宏榮1,2,王燁2,3,陸國(guó)賓1,2,李瓊芳1,2,張弘1,2,虞美秀1,2,陳蕓蕓1,2,閆方秀1,2

(1. 河海大學(xué)水文水資源學(xué)院, 江蘇 南京210098;2. 河海大學(xué)國(guó)際河流研究所,江蘇 南京210098;3. 河海大學(xué)水利水電學(xué)院, 江蘇 南京210098)

摘要：以淮河上游息縣水文站以上流域?yàn)檠芯繀^(qū)域,將基于網(wǎng)格(1 km×1 km)的新安江模型與土壤侵蝕概念模型集成,構(gòu)建基于網(wǎng)格的淮河上游泥沙負(fù)荷模擬模型,以日為時(shí)間尺度,利用息縣水文站2000—2008年實(shí)測(cè)降雨、徑流和含沙量等資料,模擬計(jì)算得到息縣水文站以上流域的泥沙負(fù)荷。結(jié)果表明,建立的泥沙負(fù)荷模擬模型能有效模擬淮河流域產(chǎn)沙過程,且模擬精度較高。研究成果可為淮河上游水土流失治理提供依據(jù),也為其他流域的水土流失模擬提供借鑒。

關(guān)鍵詞：淮河上游;土壤侵蝕；含沙量；新安江模型;泥沙負(fù)荷

隨著人口的急劇增加和科學(xué)技術(shù)水平的提高,人類正以空前的速度和規(guī)模改變著賴以生存的地球環(huán)境,地球上水的循環(huán)正在發(fā)生變化,并引發(fā)了一系列環(huán)境和生態(tài)劣變,如水土流失加劇、水環(huán)境惡化、土地荒漠化等,對(duì)人類的生存和社會(huì)的發(fā)展構(gòu)成了極大的威脅。以我國(guó)水土流失為例,2000年水土流失給國(guó)家?guī)淼慕?jīng)濟(jì)損失至少在2 000億元以上,相當(dāng)于當(dāng)年全國(guó)GDP的2.25%。因此,水土流失的治理已成為世界各國(guó)的重要課題,很多專家學(xué)者開展水土流失治理相關(guān)方面的研究,包括水土流失模擬。國(guó)外已有很多模型對(duì)流域水土流失進(jìn)行預(yù)測(cè),具有代表性的有美國(guó)農(nóng)業(yè)部提出的SWAT[1-2]模型 (該模型主要應(yīng)用于水土流失及面源污染模擬、預(yù)報(bào))。近年來,美國(guó)農(nóng)業(yè)部又提出新一代水蝕預(yù)測(cè)模型 (WEPP)[3],這是一個(gè)土壤侵蝕的物理模型。我國(guó)許多學(xué)者在流域水土流失方面也開展了大量的研究工作,具有代表性的模型有包為民[4]提出的具有物理概念的土壤侵蝕模型。但是，將廣泛應(yīng)用于我國(guó)降雨徑流模擬的新安江模型與泥沙侵蝕概念模型進(jìn)行集成,來模擬流域水土流失過程,這類研究還有待深入。

據(jù)《淮河上游水土保持規(guī)劃》,淮河流域多年的平均土壤侵蝕量為1.58億t,水土流失面積為20.1萬km2,其中上游地區(qū)的多年平均侵蝕量為1.08億t,占全流域土壤侵蝕量的68.4%。因此,開展淮河上游產(chǎn)沙模擬研究十分必要。筆者選擇淮河上游息縣水文站以上流域?yàn)檠芯繀^(qū),將新安江模型與土壤侵蝕概念模型進(jìn)行集成融合,構(gòu)建基于網(wǎng)格的分布式泥沙負(fù)荷模擬模型,根據(jù)2000—2008年實(shí)測(cè)日降雨、日徑流和日含沙量等資料,對(duì)流域徑流量、泥沙負(fù)荷時(shí)空變化全過程進(jìn)行模擬,旨在為淮河上游水土流失治理提供依據(jù),為其他流域的水土流失模擬提供參考。

1　泥沙負(fù)荷模擬模型構(gòu)建

1.1產(chǎn)水模擬

本文采用基于網(wǎng)格(1 km×1 km)的新安江模型[5],以3層蒸散發(fā)模型計(jì)算蒸散發(fā)量,以網(wǎng)格作為計(jì)算單元[4]進(jìn)行產(chǎn)流計(jì)算,以“先演后合”的匯流計(jì)算方法進(jìn)行匯流演算。坡面匯流采用線性水庫(kù)法,在計(jì)算河道匯流時(shí),將流域劃分為25個(gè)子流域,將各柵格單元在子流域出水口的總流量按照基于匯流系數(shù)的洪水演算法演算至大流域出口,然后在大流域出口斷面進(jìn)行產(chǎn)流和匯流的疊加,最終得到出口斷面的流量過程[6]。

1.2泥沙負(fù)荷模擬模型的構(gòu)建

1.2.1產(chǎn)沙模型[7]

根據(jù)產(chǎn)沙機(jī)制的差異,將流域產(chǎn)沙概化成坡面產(chǎn)沙與溝道產(chǎn)沙,則流域產(chǎn)沙量S為坡面產(chǎn)沙量和溝道產(chǎn)沙量之和:

(1)

式中:SS為坡面產(chǎn)沙量；SG為溝道產(chǎn)沙量；SSi為i時(shí)段坡面產(chǎn)沙量;SGi為i時(shí)段溝道產(chǎn)沙量;N為時(shí)段數(shù);ΔT為時(shí)段長(zhǎng)。

坡面產(chǎn)沙方程為

(2)

式中:SC為坡面水流挾沙能力;REM為流域平均的抗侵蝕能力;REMM為流域最大的抗侵蝕能力;BS為抗侵蝕能力分布曲線指數(shù);CM為全流域坡面水流的最大可能含沙量,可取常數(shù),也可取變數(shù),主要取決于坡面水流因素變化幅度的大小;R為全流域坡面平均水流深度;A為全流域坡面面積;REC為土壤的抗侵蝕能力;α為抗侵蝕能力小于REC的面積比值;α0為抗侵蝕能力為0的面積比。

溝道產(chǎn)沙量可表示為

(3)

根據(jù)拜格諾的河道水流懸移質(zhì)泥沙公式[8]來推求CG的概念性表達(dá)式從而得到SG的表達(dá)式:

(4)

式中:CG為溝道水流含沙量;Q為溝道水流流量;CGM為平均溝道產(chǎn)沙濃度;LQ為 ln(Q+1)的時(shí)間平均值;BV為常參數(shù)。

1.2.2匯沙模型

泥沙的運(yùn)動(dòng)是建立在水流運(yùn)動(dòng)之上的,而水流中的泥沙含量又能夠影響水流的速度,進(jìn)而影響水流的運(yùn)動(dòng),二者是相輔相成的。因此,參考流域水流的匯流模型,建立流域匯沙模型。匯沙包括坡面匯沙與溝道匯沙兩個(gè)部分。坡面匯沙和溝道匯沙與水流相似,也存在著泥沙平衡方程和泥沙蓄泄關(guān)系。

泥沙平衡方程為

(5)

泥沙蓄泄方程為

(6)

將式(1)進(jìn)行差分并代入式(2),可得到坡面匯沙方程:

(7)

溝道匯沙方程考慮到河道內(nèi)泥沙的沖刷和淤積,因此仿照水流分段馬斯京根法可得:

(8)

其中

式中:WS為泥沙蓄量;IS為坡面時(shí)段產(chǎn)沙量;IS1、IS2分別為坡面時(shí)段1和時(shí)段2的產(chǎn)沙量;So為坡面時(shí)段出沙量;So1、So2分別為坡面時(shí)段1和時(shí)段2的出沙量;KS為坡面的匯沙時(shí)間;CSS為匯沙系數(shù);Sg1、Sg2分別為溝道時(shí)段1和時(shí)段2的出沙量；Ig為溝道時(shí)段產(chǎn)沙量；Ig1、Ig2分別為溝道時(shí)段1和時(shí)段2的產(chǎn)水得；b0、b1和b2均為過程系數(shù);KES為河段內(nèi)的泥沙傳播時(shí)間;XES為泥沙比重系數(shù)。

1.3產(chǎn)水與產(chǎn)沙模型的集成融合

河道泥沙演算方程中需要輸入河道流入斷面和流出斷面的流量,筆者將把基于網(wǎng)格(1 km×1 km)的新安江模型模擬計(jì)算的流量作為河道泥沙演算方程中的河道流入斷面和流出斷面的流量,從而實(shí)現(xiàn)泥沙負(fù)荷模擬模型與新安江模型的集成與融合。

2　泥沙負(fù)荷模擬模型在淮河上游的應(yīng)用

2.1研究區(qū)域概況

淮河流域地處我國(guó)南北氣候過渡帶,特殊的過渡性氣候和大規(guī)模人類活動(dòng)對(duì)地表的強(qiáng)烈擾動(dòng),造成了嚴(yán)重的水土流失。研究區(qū)息縣(圖1)以上流域位于東經(jīng)113°15′～114°46′,北緯31°31′～32°43′,流域面積為10 190 km2,多年平均水面蒸發(fā)量800～1 000 mm,多年平均降水量800～1 400 mm,多年平均徑流深約371 mm。流域內(nèi)主要土地利用類型是旱地、林地、水田和灌木叢。對(duì)淮河息縣水文站1964—2008年的泥沙資料進(jìn)行分析,發(fā)現(xiàn)年水土流失量主要取決于汛期月份降雨量和暴雨出現(xiàn)次數(shù),泥沙負(fù)荷汛期(6—9月)占全年的63.4%。以水土流失為載體的面源污染已成為流域上游主要水源區(qū)水質(zhì)污染的重要來源,危及流域水安全與流域生態(tài)系統(tǒng)健康。

圖1　研究區(qū)域

2.2資料處理

研究所需數(shù)據(jù)資料包括地形、土地利用、水文、氣象等方面的資料,數(shù)據(jù)來源與處理方法是:①采用地圖投影對(duì)模型中輸入的所有空間數(shù)據(jù)的坐標(biāo)系統(tǒng)進(jìn)行統(tǒng)一,考慮到研究區(qū)域?qū)儆谥懈呔暥鹊貐^(qū),采用在該地區(qū)變形很小的等角橫軸割圓柱投影[9],即UTM,Universal Transverse Mercator投影。②地形資料選用美國(guó)國(guó)家地球物理中心提供的、全球分辨率為1 km×1 km的數(shù)字高程數(shù)據(jù),運(yùn)用ArcGIS軟件中的掩膜工具,結(jié)合流域形狀圖,切割得到研究流域的數(shù)字高程圖。③土地利用類型空間分布及面積比例的確定,采用中國(guó)科學(xué)院提供的2 000 s分辨率為1 km×1 km的土地利用圖,結(jié)合流域邊界,切割得到研究流域2 000 s的土地利用圖。

2.3 模型率定和檢驗(yàn)

選定2000—2005年為率定期、2006—2008年為檢驗(yàn)期,以日為時(shí)間尺度對(duì)新安江模型和泥沙負(fù)荷模擬模型分別進(jìn)行率定和檢驗(yàn)。參數(shù)率定結(jié)果和模型模擬結(jié)果分別見表1～2和表3～4。

表1　息縣流域新安江模型參數(shù)率定結(jié)果

表2　息縣流域泥沙模擬參數(shù)率定結(jié)果

2.4結(jié)果分析

表2中各參數(shù)的敏感性分析采用比較函數(shù)對(duì)各參數(shù)的敏感度[10]進(jìn)行分析的方法。本文主要列出CM和CGM的敏感性分析圖(其余參數(shù)經(jīng)分析后不敏感,故此處不列出),見圖2。由表2可知,流域坡面水流最大可能含沙量CM=0.17 kg/m3,平均溝道產(chǎn)沙濃度CGM=0.6 kg/m3。由圖2可知這兩個(gè)參數(shù)的敏感性很強(qiáng),其微小的改變都能夠?qū)Y(jié)果產(chǎn)生較大的影響,且敏感性隨著參數(shù)的增大而增大,因此,調(diào)整參數(shù)時(shí)要特別注意這兩個(gè)參數(shù)。

表3　息縣流域日徑流模擬結(jié)果

表4　息縣流域泥沙日模型模擬結(jié)果

圖2　CM和CGM敏感性分析

由表4可知,在率定期與驗(yàn)證期模型模擬的年輸沙量的相對(duì)誤差基本都在20%以內(nèi)(除2001年以外),確定性系數(shù)基本都在0.5以上(除2001年以外)。分析其原因,首先是2001年的降水量?jī)H568.9 mm,小于研究區(qū)域多年平均降水量800～1 400 mm,屬于來水量較枯年份,在無人類活動(dòng)對(duì)流域進(jìn)行破壞干擾的前提下,人類對(duì)流域每年的需水量相對(duì)穩(wěn)定,但2001年屬于來水量較枯年份,在扣除人類的需水量后,來水量更少了,而流域的泥沙運(yùn)動(dòng)是建立在水流運(yùn)動(dòng)之上的,所以實(shí)測(cè)輸沙量也相應(yīng)較小,僅0.9萬t,遠(yuǎn)遠(yuǎn)小于其他年份(流域多年平均輸沙量為150萬～200萬t),因此其相對(duì)誤差較大,確定性系數(shù)較小,屬于模擬效果較差的年份,其他年份的模擬結(jié)果都能夠滿足精度要求(相對(duì)誤差在20%以內(nèi),確定性系數(shù)大于0.5)。

由圖3可知,模型模擬的相關(guān)系數(shù)為0.98,二者呈強(qiáng)正相關(guān)性,表明數(shù)據(jù)的規(guī)律比較可靠,故其日尺度的泥沙負(fù)荷模擬可滿足泥沙模型的精度要求。

圖3　泥沙模擬相關(guān)性分析

圖4　息縣以上流域2000、2006、2008年模擬與實(shí)測(cè)的日輸沙量過程線

在本模型中設(shè)為常數(shù),但在實(shí)際情況下,降雨強(qiáng)度能影響CM值的大小,降雨強(qiáng)度較大的時(shí)候,相應(yīng)的CM值也相應(yīng)較大,降雨強(qiáng)度與CM之間呈非線性關(guān)系。此外,比較表3和表4可知,泥沙模擬結(jié)果差于徑流模擬結(jié)果,這是因?yàn)槟嗌衬M是建立在徑流模擬的基礎(chǔ)之上的,徑流模擬的誤差在泥沙模擬時(shí)得到了放大,因此泥沙模擬相對(duì)誤差會(huì)比徑流模擬相對(duì)誤差大。

3　結(jié)　語

將新安江模型與土壤侵蝕概念模型進(jìn)行集成與融合,構(gòu)建了基于網(wǎng)格的淮河上游泥沙負(fù)荷模擬模型，并模擬了淮河上游息縣水文站以上流域的泥沙負(fù)荷變化過程。研究結(jié)果表明,模型參數(shù)除CM和CGM比較敏感外,其余參數(shù)皆不敏感;在率定期與驗(yàn)證期模型模擬的年輸沙量的相對(duì)誤差基本都在20%以內(nèi)(除2001年以外),確定性系數(shù)基本都在0.5以上(除2001年以外),模型模擬的輸沙量與實(shí)測(cè)輸沙量的相關(guān)系數(shù)為0.98,二者呈強(qiáng)正相關(guān)性?？梢?日尺度的泥沙負(fù)荷過程模擬效果較好,可用于該地區(qū)的產(chǎn)沙過程模擬,也可為其他地區(qū)的產(chǎn)沙負(fù)荷模擬提供參考,為流域水土流失治理和科學(xué)資源開發(fā)利用提供一定的依據(jù)。

參考文獻(xiàn)：

[1] NEITSCH S L, ARNOLD J G, KINIRY J R, et al. Soil and water assessment tool user’s manual version 2000[J]. GSWRL Report, 2002, 202：22-26.

[2]王中根,劉昌明,黃友波.SWAT模型的原理、結(jié)構(gòu)及應(yīng)用研究[J].地理科學(xué)進(jìn)展,2003,22(1):79-86.(WANG Zhonggen,LIU Changming,HUANG Youbo.Principle,structure and application of SWAT model[J]. Journal of Geographical Science, 2003, 22 (1): 79-86. (in Chinese))

[3]張玉斌,鄭粉莉,賈媛媛.WEPP模型概述[J].水土保持研究,2004,11(4):146-149.(ZHANG Yubin,ZHENG Fenli,JIA Yuanyuan.Summary of WEPP model[J]. Research of Soil and Water Conservation, 2004, 11 (4): 146-149. (in Chinese))

[4]包為民.水文預(yù)報(bào)[M].3版. 北京:中國(guó)水利水電出版社,2006.

[5]姚成.基于柵格的分布式新安江模型構(gòu)建與分析[D].南京:河海大學(xué),2007.

[6]李瓊芳,謝偉,薛運(yùn)宏,等.新安江模型在土壤侵蝕模擬中的應(yīng)用[J].水電能源科學(xué),2010, 28(3): 11-13.(LI Qiongfang,XIE Wei,XUE Yunhong,et al.Xin’anjiang model in the application of soil erosion simulation[J]. Water Research and Power, 2010,28 (3): 11-13. (in Chinese))

[7]包為民.黃土地區(qū)流域水沙模擬概念模型與應(yīng)用[M].南京:河海大學(xué)出版社,1995.

[8]BRUNASH J C,FERRAL R L,MCGUIRE A,et al.A generalized streamflow simulation system:conceptual modeling for digital computers[R]. Sacramento,CA:Joliet Federal State River Forecasts Center,1973.

[9]莊衛(wèi)東,汪春.精準(zhǔn)農(nóng)業(yè)中 UTM 投影及反算應(yīng)用研究[J].黑龍江八一農(nóng)墾大學(xué)學(xué)報(bào), 2005, 17(3): 47-50.(ZHUANG Weidong,WANG Chun.Studies on the UTM projection and it’s inverse calculation in the precision agriculture[J].Journal of Heilongjiang Bayi Agricultural University, 2015,17(3):47-50.(in Chinese))

[10]李毅,邵明安,王文焰,等.土壤非飽和導(dǎo)水率模型中參數(shù)的敏感性分析[J].水科學(xué)進(jìn)展,2003,14(5):593-597.(LI Yi,SHAO Ming’an,WANG Wenyan,et al.Unsaturated soil hydraulic conductivity the sensitivity analysis of parameters in the model[J]. Advance in Water Science, 2003, 14 (5): 593-597.(in Chinese))

由于篇幅有限,圖4只列出2000年、2006年和2008年的模擬與實(shí)測(cè)的日輸沙量過程線。由圖4可以看出,模型模擬的計(jì)算輸沙量存在一個(gè)普遍現(xiàn)象,即主峰值上不去。分析其原因,可能是因?yàn)閰?shù)CM(坡面水流最大含沙量,隨著坡面水流流速的變化而變化)

DOI：10.3880/j.issn.1004-6933.2016.04.013

基金項(xiàng)目:國(guó)家自然科學(xué)基金(41171220);長(zhǎng)江學(xué)者和創(chuàng)新團(tuán)隊(duì)發(fā)展計(jì)劃資助(IRT13062);水利部公益行業(yè)經(jīng)費(fèi)項(xiàng)目(201001069,201101052)

作者簡(jiǎn)介:鄒宏榮(1990—),男,碩士研究生,研究方向?yàn)樯鷳B(tài)水文及水環(huán)境保護(hù)。E-mail: 453551965@qq.com 通信作者:李瓊芳,教授,博士生導(dǎo)師。E-mail:qfli@hhu.edu.cn

中圖分類號(hào)：TV143+.4

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1004-6933(2016)04-0080-04

(收稿日期：2015-07-09編輯：彭桃英)

Simulation of sediment yield in upper reaches of Huaihe River

ZOU Hongrong1,2, WANG Ye2,3, LU Guobin1,2, LI Qiongfang1,2, ZHANG Hong1,2,YU Meixiu1,2, CHEN Yunyun1,2, YAN Fangxiu1,2

(1. College of Hydrology and Water Resources, Hohai University, Nanjing 210098, China;2.InstituteofInternationalRiversResearchAcademy,HohaiUniversity,Nanjing210098,China;3.CollegeofWaterConservancyandHydropowerEngineering,HohaiUniversity,Nanjing210098,China)

Abstract：This study was carried out in the basin located above the Xixian Hydrological Station in the upper reaches of the Huaihe River. A grid-based (1 km×1 km) Xin’anjiang model was coupled with the conceptual model of soil erosion to build the sediment load simulation model for the upper reaches of the Huaihe River. Using the daily measured data of rainfall, runoff, and sediment concentration at the Xixian Hydrological Station from 2000 to 2008, the sediment load in the basin above the Xixian Hydrological Station was calculated. The results show that the built model can be effectively used to simulate the sediment yielding process in the Huaihe Basin, and the model has high precision. This study provides a scientific basis for the control of water loss and soil erosion in the upper reaches of the Huaihe River and provides references for the simulation of water loss and soil erosion in other basins.

Key words：upper reaches of Huaihe River; soil erosion; sediment concentration; Xin’anjiang model; sediment load