西柳溝丘陵區土壤組成及對流域產沙的影響

王普慶 侯素珍

摘?要：為進一步了解和認識不同地貌單元組成土壤對流域產沙的影響，選擇西柳溝為典型流域，針對地表物質分布開展現場調研取樣，分析西柳溝丘陵區坡面土壤組成以及不同土壤產沙對流域輸沙的影響。運用泥沙動力學理論，利用河道斷面資料和泥沙參數，結合流域的洪水輸沙特性，求出相應條件下西柳溝流域出口段河流懸沙上限粒徑值為0.12 mm。在一般洪水條件下，西柳溝丘陵地區粗骨土、栗鈣土進入黃河的比例分別為18%、20%。

關鍵詞：西柳溝;土壤組成;顆粒分析;產沙

中圖分類號：P333;TV882.1?文獻標志碼：A

doi：10.3969/j.issn.1000-1379.2020.02.001

Abstract： In order to further understand the influence of different geomorphic units to the sediment yield of watersheds， Xiliugou was selected as a typical watershed for field survey on the distribution of superficial substances to analyze the slope soil grain sizes in the gullied Rolling loess area of Xiliugou and the influence of different soils to the sediment yield of watersheds. Based on the dynamic theory of sediment， the river section and sediment data， combined with the sediment transport characteristics in floodwater of the watershed， the maximum grain size of river suspended sediment is 0.12 mm under corresponding conditions through calculation. Under the normal conditions of flood flow of 3 300 m3/s with a sediment concentration of 600 kg/m3， about 20% of skeletal soil，18% Chestnut soil and 6.8% sand of the Xiliugou tributary are carried into the Yellow River by flood.

Key words： Xiliugou; soil grain size; grain analysis; sediment yield

西柳溝是鄂爾多斯高原上的“十大孔兌（孔兌為蒙語，即山洪溝）”之一，位于內蒙古河套平原的中部，發源于鄂爾多斯市東勝區，是黃河的一級支流（見圖1），流域總面積1 356.43 km2。地勢南高北低，相對高差為500 m，上陡下緩，河道平均比降為0.36%。河道地形地貌從南向北呈現不同特征，南部（上游）為丘陵溝壑區，面積876.30 km2，水力侵蝕比較嚴重;中游為庫布齊沙漠風沙區，面積280.70 km2;下游為沿黃平原區，面積為199.43 km2。這樣的地貌格局和地表物質分布特征形成了特有的泥沙侵蝕和輸移模式：上游為水力作用下的高含沙水流主導區;中游沙漠區河段為風力作用主導區，區域內泥沙通過風力作用被搬運至床面;下游為泥沙輸移與沉積區。

大量研究成果[1-2]表明，孔兌洪水泥沙對干流影響深遠。同時，孔兌洪水驟漲驟落、含沙量高度集中輸移的特點，極易在干支流交匯處大量淤積形成沙壩阻塞黃河，帶來嚴重災害。因此，正確認知孔兌洪水泥沙的來源，科學合理地調控洪水和泥沙，實現洪水和泥沙資源化，對于流域綜合治理和區域經濟發展具有重要的意義。

坡面和溝道來沙是河流泥沙的主要來源，經過水土保持生態治理措施的實施，坡面土壤侵蝕已得到有效控制。依據部分小流域的土壤侵蝕資料分析[3]，從20世紀90年代到2010年，坡面和溝道侵蝕產沙比分別為70%、80%，溝道侵蝕作為流域產沙的主要來源呈增大趨勢，流域不同位置對流域產沙的貢獻也不同。另外，人類活動相對頻繁的道路、排水溝等也是河流泥沙的重要來源和泥沙輸送的重要通道。盡管溝道侵蝕量在不斷增大，但坡面侵蝕依然是溝道侵蝕的初始來源地，研究坡面土壤組成及其對河流泥沙的貢獻一直是人們關注的主要問題，這對于調整土地利用方式及流域規劃是十分有意義的。

1?丘陵區不同土壤分布及粒配

西柳溝流域位于十大孔兌中部，水文氣象、洪水泥沙和地質地貌資料比較齊全，流域地形地貌特征、表層土壤組成與其他孔兌也相似，在水力侵蝕和風力侵蝕等方面具有一定的代表性，故把西柳溝流域作為典型流域進行研究。

1.1?地表土壤分布

西柳溝流域地貌類型復雜多樣，依據地貌發育的內外營力及地表形態特征，可以把西柳溝丘陵區劃分為蓋土砒砂巖區、蓋沙砒砂巖區和裸露砒砂巖區3個地貌類型區。土壤主要是由砒砂巖衍生而形成的栗鈣土、粗骨土、風沙土和其他類型土[4-5]。

圖2為依據遙感圖像解析的西柳溝表層土壤分布情況，自下而上分布有草原風沙土、新積土、栗鈣土、淡栗鈣土、粗骨土、潮土6類土壤（見表1）。草原風沙土位于西柳溝流域的中游和上游，占西柳溝流域面積的16%，其中上游部分為流域發源地，形似啞鈴，橫向長5.5～6.0 km;栗鈣土分布于流域干流兩側，占流域面積的53.7%，屬于水力侵蝕的主要區域;粗骨土分布于流域中上游兩側，呈南北長、東西窄的帶狀形態，覆蓋了西柳溝大部分支流的上游，再加上該區域植被覆蓋率較低，荒坡面積大，屬于水力侵蝕強度最大的地區。

1.2?不同土壤粒配

1.2.1?取樣位置和方式

大量研究成果[6-7]表明，侵蝕泥沙主要來源于溝谷地坡面、溝道。對于坡面沙樣選取，可以參考土壤分布圖預定的位置，并結合現場情況，根據洪水痕跡確認溝坡和溝床坡面，由于土壤本身存在著空間分布的不均一性，因此為使沙樣更好地反映區域土壤性狀，采樣點選擇應遵從影響土壤侵蝕發生、發展的主要條件的相似性原則、地貌部位相似性原則、侵蝕方式一致性原則，同時注意采樣點避開田邊、溝邊、路邊等。

取樣線路采用蛇形分布，布點做到盡量均勻和隨機。采用混合樣品采集布點方法，以塊地為單位多點取樣，然后再混合成一個混合樣品，取樣深度為0～20 cm。

1.2.2?土壤粒配

考慮不同地表土壤組成，按照取樣位置和方式的要求，分別在西柳溝東岸的支溝艾來五庫溝、昌漢溝及西岸的馬利昌漢溝等處取樣，編號依次為D1～D7（見圖2）。圖3為各沙樣粒配組成曲線，其中風沙土組成最粗，中值粒徑為1.10 mm（見表2），其次為編號D4的栗鈣土，中值粒徑為0.42 mm，最細的為編號D7的粗骨土，中值粒徑為0.09 mm，土壤粒徑組成并不完全取決于其分類組成，而是與其所處的地理位置和土層深度有關。將丘陵區同一類型土壤的沙樣合并后得到不同土壤的分組粒徑含量（見表3），總體看，風沙土粒配組成最粗，栗鈣土次之，粗骨土相對較細。丘陵區內栗鈣土、粗骨土、草原風沙土三者合計占95.3%，是主要的水力侵蝕區，也是對流域產沙影響最大的區域，是本文研究的重點。

2?孔兌洪水特性

“十大孔兌”雖然為季節性河流，但是其地貌格局和地表物質分布組成對高含沙水流的形成和風水兩相的侵蝕—泥沙搬運作用影響很大，因此形成了特有的泥沙侵蝕和輸移模式。

在洪水泥沙方面，文獻[8]統計了西柳溝龍頭拐站1960—1990年洪水流量和沙量的多年變化，在這31 a中，最大洪峰流量為6 940 m3/s，最大含沙量為1 550 kg/m3，其中有12 a的最大含沙量超過1 000 kg/m3，場次洪水平均含沙量為580 kg/m3。文獻[8]對于西柳溝部分場次洪水特征值進行了統計，而且多數場次洪水過程歷時較短（不超過24 h）。

為分析洪水泥沙組成，借助發源于砒砂巖地區的皇甫川洪水泥沙測驗資料，點繪了皇甫川洪水泥沙懸沙中值粒徑與含沙量的關系（見圖4），可看出洪水挾帶的泥沙粒徑隨含沙量的增大而增大[9]，中值粒徑可以達到0.3 mm。

在孔兌河床淤積物組成方面，根據西柳溝河床淤積物沙樣粒配分析，結果顯示床沙組成沿程細化[1]。上游河段床沙組成偏粗，隨著河床的不斷淤積抬升，到下游平原河段則細顆粒泥沙偏多，泥沙中值粒徑從上游的2 mm減小到下游的0.157 mm，這顯然是高含沙洪水沿程淤積分選的結果。綜合以上分析，孔兌洪水具有洪峰流量大、含沙量高、歷時短的基本特征。在懸沙組成方面，粒徑隨含沙量的增大而增大，在河床淤積物組成方面表現為沿程細化。

3?流域泥沙來源分析

3.1?泥沙來源研究現狀

砒砂巖不同侵蝕類型區小流域泥沙來源的研究主要體現在兩方面：一為物源說，即泥沙的來源問題;二是貢獻率問題，即水力和風力侵蝕沙量占比及對黃河的影響。王曉[10]對砒砂巖不同侵蝕類型區小流域泥沙來源的分析表明，泥沙主要來源于溝谷地;楊明義等[11]應用137Cs研究了安塞紙坊溝流域的泥沙來源，表明小流域泥沙主要來源于溝谷;張勝利[12]認為只有風蝕、水蝕交錯地區才是強烈侵蝕產生粗泥沙的中心;史學建等[13]對比分析了黃河中游皇甫川等流域河流粗泥沙含量和平均粒徑的年內變化規律，認為流域內基巖產沙對河流泥沙貢獻不大，懸移質中的粗泥沙主要來源于風沙區。

關于孔兌泥沙入黃貢獻率的研究成果相對有限，文獻[1，13]從不同的角度較全面地分析闡述了小流域泥沙來源。文獻[1]利用床沙礦物成分分析、斷面分析等方法來推算西柳溝與孔兌總輸沙量，從而得到孔兌輸沙量中沙漠沙的占比，整體上十大孔兌沙漠沙的貢獻率接近20%，其他80%則來自丘陵地區的溝道、坡地等。文獻[14]利用137Cs估算出流域內各類用地單元的平均侵蝕模數，再與相應的面積相乘，得到整個小流域的多年平均侵蝕量，從而獲得流域泥沙的主要來源，結果表明對流域產沙貢獻最大的是荒坡地和沙地，兩者占比達76%，與文獻[1]得出的結論接近。綜上所述，西柳溝流域泥沙來源約20%來自中游的庫布齊沙漠，其余來沙則來源于荒坡、溝道等，然而針對不同土壤的坡面產沙對黃河干流的影響目前尚無具體的結論性成果。

3.2?泥沙輸移理論分析

考慮到孔兌下游屬于平原地區寬淺河道，符合多沙河流輸水輸沙的基本規律，依據泥沙動力學相關理論，應用多沙河流輸沙平衡計算方法[15]，可推求流域產沙后進入黃河干流的沙量。式（1）為維持河流懸沙不淤的最小流速公式，該式經大量水槽試驗及渠道、河道實測資料驗證，表明該式用于計算明渠高含沙水流挾沙力有一定普遍意義，由于驗證資料包含侵蝕強烈的黃土丘陵地區，因此可用于孔兌地區河道高含沙水流。應用此式，結合河道斷面水力要素和泥沙參數，可求出相應條件下懸沙上限粒徑，進一步估算孔兌進入黃河干流的泥沙量。

當河流處于輸沙平衡狀態時，水流含沙量與水流挾沙力相同。依據西柳溝平原區河段的過流斷面情況，應用式（1）～式（5）可得到不同含沙量條件時洪水流量與懸沙粒徑的關系（見圖5）。參照平原河段設防標準為20 a一遇的洪水流量（約3 300 m3/s）及多年平均含沙量（600 kg/m3），與上述水沙條件對應的懸沙上限粒徑值為0.12 mm（見圖5），即大于該粒徑的懸沙將會在河床上沉積[16]。按照該粒徑計算，結合表3粒配資料，得到西柳溝孔兌約有15%的粗骨土、17%的栗鈣土、6.8%的風沙土進入黃河。

結合西柳溝歷史洪水情況，將流量2 000 m3/s和7 000 m3/s分別作為中小洪水和極端特大洪水情況考慮，對應洪水含沙量變化范圍通常為300～1 000 kg/m3。由圖5可以查取孔兌流域出口斷面的粒徑，其中，中小洪水時孔兌出口斷面的粒徑為0.113～0.115 mm，與該粒徑對應的粗骨土、栗鈣土、風沙土含量分別為13.5%、16.0%、7.5%;特大洪水時，孔兌出口斷面的粒徑為0.128～0.14 mm，對應的粗骨土、栗鈣土、風沙土含量分別為22%、25%、8.2%。將兩類洪水產沙的情況平均后可得到粗骨土、栗鈣土、風沙土含量分別約為18%、20%、7.9%。因此，一般情況下，西柳溝丘陵地區在水力侵蝕條件下，占比份額較大的粗骨土和栗鈣土進入黃河的沙量大約不超過侵蝕總量的20%。

4?結?語

泥沙從侵蝕到沉積是一個復雜的過程，通過對西柳溝流域實地查勘、取樣、測試、分析計算，針對不同地貌單元的土壤組成，分析計算了西柳溝流域不同地貌單元的粒配組成。結合流域的洪水輸沙特性及泥沙動力學理論，應用河道斷面資料和泥沙參數，求出相應條件下河流懸沙的上限粒徑，在洪水流量為3 300 m3/s、含沙量為600 kg/m3時，對應懸沙上限粒徑值為0.12 mm。西柳溝約有15%的粗骨土、17%的栗鈣土和6.8%的風沙由洪水挾帶進入黃河。在其他洪水條件下，西柳溝丘陵地區粗骨土、栗鈣土進入黃河的含量分別約為18%、20%。

西柳溝流域是典型的風水兩相作用地區，由于流域內地形地貌復雜，支流眾多，再加上缺乏系統的洪水水文泥沙資料，因此目前條件下精確定量計算流域產沙情況比較困難。本文旨在通過不同的思路探討孔兌泥沙來源，并為其他類似條件下的泥沙輸移計算提供理論參考和技術途徑。

參考文獻：

[1]?侯素珍，王平，林秀芝，等.黃河內蒙古段孔兌高濃度挾沙洪水調控技術研究[R].鄭州：黃河水利科學研究院，2015：218-257.

[2]?梁其春，駱鴻固，王英順，等.十大孔兌區水土流失現狀與防治對策[J].中國水土保持，1996（3）：3-6.

[3]?馬瞳宇，李譚寶，張平倉.近20年水蝕風蝕交錯區小流域侵蝕及坡溝產沙演變的粒徑對比分析[J].水土保持學報，2013，27（4）：83-87.

[4]?范昊明，蔡強國.流域侵蝕產沙平衡研究進展[J].泥沙研究.2004（4）：76-81.

[5]?許炯心.黃河中游多沙粗沙區高含沙水流的粒度組成及其地貌學意義[J].泥沙研究，1999（5）：13-17.

[6]?侯素珍，王平，楚衛斌.黃河上游水沙變化及成因分析[J].泥沙研究，2012（4）：47-53.

[7]?王平，田勇，侯素珍，等.黃河內蒙古段孔兌水沙特點分析[J].人民黃河，2012，34（11）：39-42.

[8]?王普慶.內蒙古西柳溝入河口流路變遷分析[J].中國防汛抗旱，2017，27（6）：57-62.

[9]?王平，侯素珍，張原鋒，等.黃河上游孔兌高含沙洪水沖淤特性[J].泥沙研究，2013（1）：67-73.

[10]?王曉.“粒度分析法”在小流域泥沙來源研究中的應用[J].水土保持研究，2002，9（3）：42-43.

[11]?楊明義，田均良，劉普靈，等.應用137Cs研究小流域泥沙來源[J].土壤侵蝕與水土保持學報，1999，5（3）： 49-53.

[12]?張勝利.黃河中游多沙粗沙區風蝕產沙對黃河粗泥沙的影響分析[J].中國水土保持，1997（9）：13-18.

[13]?史學建，劉宇梁，黃靜，等.再談黃河中游粗泥沙的來源[J].人民黃河，2007，29（2）：62-63.

[14]?姚文藝.黃河沙漠寬谷河道水沙變化及驅動機理[M].北京：科學出版社，2018：329-335.

[15]?費祥俊，邵學軍.泥沙源區溝道輸沙能力的計算方法[J].泥沙研究，2004（1）：1-8.

[16]?侯素珍，李婷，郭彥，等.西柳溝河床泥沙組成分析[J].人民黃河，2018，40（7）：7-11.

【責任編輯?翟戌亮】