黃土丘陵區小流域壩地沉積泥沙顆粒分布特征

申震洲, 亢晨波, 張寧寧, 劉普靈

(1.黃河水利科學研究院 水利部 黃土高原水土流失過程與控制重點實驗室, 鄭州 450003;2.榆林學院 陜西省陜北礦區生態修復重點實驗室, 陜西 榆林 719000;3.中國科學院 水利部 水土保持研究所 黃土高原土壤侵蝕與旱地農業國家重點實驗室, 陜西 楊凌 712100)

土壤侵蝕是全世界所面臨的頭號環境問題[1]。黃土丘陵區由于溝壑縱橫、溝深坡陡的復雜地形地貌條件，導致該地區水土流失范圍廣、面積大、土壤侵蝕嚴重，侵蝕的泥沙成為黃河泥沙的主要來源[2]，進而導致生態失衡，旱澇災禍頻發。為保障黃河安全，改善生態環境，以小流域為侵蝕單元進行的小流域綜合治理措施成為黃土高原土壤侵蝕治理的重要手段[3]，而淤地壩是其中一項成效顯著的土壤侵蝕治理工程措施，它在攔泥、滯洪、淤沙造田等方面起著重要作用[4-5]。壩地系統所攔蓄泥沙顆粒賦存了大量土壤侵蝕產沙過程及侵蝕環境演替的相關信息，是研究小流域土壤侵蝕變化特征及反演侵蝕環境變化的良好載體，為研究小流域土壤侵蝕提供了新的途徑[6-7]。

淤地壩內泥沙淤積特征受區域降雨特性、水動力條件和研究小流域地形地貌特征等多種因素影響，因此壩前沉積物厚度和顆粒組成等在空間上存在很大的變異性，進而影響典型剖面的選擇，并進一步影響沉積量的估算和泥沙來源的追蹤和反演[8-9]。目前許多學者通過小流域淤地壩內的沉積物顆粒分析開展了小流域土壤侵蝕過程、產沙來源等方面的研究，并取得了一些研究成果[10-13]。目前研究主要基于137Cs進行時標，然后根據旋回層厚度，建立泥沙沉積過程的時空序列，并結合歷次暴雨、洪水資料及流域內的沉積泥沙量來分析流域內的泥沙來源和產沙強度[14-16]。Zhang等[17]根據淤地壩沉積泥沙所形成的沉積旋回層厚度，并結合137Cs和孢粉等相關技術分析了小流域在不同時間內的侵蝕產沙速率。然而由于建壩年代限制，并非所有壩地都可用137Cs示蹤斷代技術，但每一個沉積旋回都對應著一次明顯的侵蝕性降雨事件，也相應的會帶來一次泥沙的運移和分選過程，因此研究分析壩地的沉積泥沙顆粒的機械組成顯得十分必要。彌智娟等[18]通過對壩地沉積泥沙的分析，發現壩地沉積泥沙粒級組成中粉粒占主導地位。一些學者通過對壩地內不同分層沉積物的研究分析表明沉積泥沙具有分選性，壩前泥沙顆粒比壩尾的細[19-20]。然而有的則認為與此相反，即壩前的顆粒比壩尾的粗，或無明顯差異[21]。

鑒于此，本文選用延安湫溝的一座2003年興建“悶葫蘆壩”為研究對象，分析研究小流域控制壩沉積泥沙垂直剖面泥沙粒徑組成變化及分形特征，探討新的侵蝕背景下小流域侵蝕泥沙再分布特征，以期能夠為開展和研究小流域土壤侵蝕及產沙過程提供理論依據。

1 材料與方法

1.1 研究區概況

研究小流域位于陜西省延安市湫溝小流域，屬于燕溝流域的一條支溝，地理位置為36°29′96″N，109°31′11″E，海拔1 131 m；淤地壩控制面積為5.23 km2。該淤地壩始建于2003年，由于淤積時間較短，該淤地壩仍然處于有水狀態，且并未出現放水現象。因此，所有小流域侵蝕泥沙全部都被攔蓄在該淤地壩內。該地區的氣候特征為半干旱大陸性季風氣候，多年平均降水量為423 mm，主要集中在7—9月份，且多為暴雨。研究區的土壤侵蝕類型主要是水力侵蝕和重力侵蝕。流域內土壤類型以黃綿土為主，主溝道最底層為裸露的基巖。

1.2 土壤樣品采集和處理

1.2.1 樣點布設與土壤樣品采集 為了研究淤地壩降雨沉積過程，選擇延安市的湫溝小流域的湫溝大壩，由于淤積時間較短，該淤地壩仍然處于有水狀態。經過前期調查，于2012年選擇淤地壩淤積泥沙最厚的點進行淤積泥沙樣品的采集，采樣使用內徑為42 mm的食品級不銹鋼空心管進行水上采樣。在采樣的過程中，將空心管垂直打入，直到原始河床(即基巖位置)，隨后垂直取出，同時根據采樣點所處位置的水面深度，當場對鋼管進行焊接和截取，以確保采樣鋼管內土柱樣品的完整性，密封帶回實驗室。在采樣的過程中進行GPS定位，記錄采樣點的位置信息。

1.2.2 土壤樣品處理 將帶回的空心管土柱樣品取出上部密封物體后，垂直放置，自然風干，待空心管內樣品成半干燥狀態后沿著管道中心線垂直剖開(該過程要特別小心，盡量減少對空心管內土柱的影響)。將剖開的土柱樣品進行0.5 cm/樣的切割取樣。隨后將土樣處理后裝入密封袋內進行保存，用于土壤顆粒組成的分析，并注明采樣時間以及采集的土樣所處的剖面深度。

1.3 土壤樣品分析

將土柱分割所得的土壤樣品平鋪在通風口自然風干、除去草石雜物、研磨、過篩(2 mm的篩子)。在100 ml的燒杯中加入2～4 g土壤樣品，并加入過氧化氫(H2O2)溶液和10%的鹽酸(HCL)溶液以消除對土樣中的有機質和碳酸鹽。最后給每個土壤樣品加入六偏磷酸鈉(NaPO3)6溶液，進行分散處理。對處理的土樣使用Mastersizer 2000激光粒度儀進行粒度分析，每個土樣重復測量兩次，結果取平均值。土壤粒徑分類采用1987年的中國制標準，分別為黏粒(<0.002 mm)、細粉粒(0.002～0.01 mm)、粗粉粒(0.01～0.05 mm)、細沙粒(0.05～0.25 mm)、粗沙粒(0.25～1 mm)。本文采用土壤分形維數和土壤質地粗化度綜合描述顆粒分布特征。

土壤分形維數可表征土壤粒徑的大小組成，還能反映質地的均一程度，分形維數隨著土壤黏粒和粉粒含量的增加，其值越來越大；隨著土壤粒徑組成中沙粒含量的增大而逐漸減小。計算公式如下：

式中:D表示土壤粒徑分形維數;V(r

土壤質地粗化度指的是壩地剖面沉積旋回層中土壤的沙粒含量(>0.05 mm)與黏粒含量和粉粒含量之和(<0.05 mm)的比值。粗化度越大，則粗顆粒泥沙在土壤粒級組成中的成分越多；粗化度越小，則土壤粒級組成中的粗顆粒含量越少，土壤受侵蝕程度越大。

1.4 數據處理

本試驗原始數據的處理和制圖使用Excel 2010軟件完成，差異顯著性檢驗使用SPSS 19.0軟件完成。

2 結果與分析

2.1 淤地壩沉積泥沙再分布特征

流域壩地侵蝕類型主要以水力侵蝕與重力侵蝕為主，每一場侵蝕性降雨都會形成一次泥沙搬運沉積現象，同時形成一個泥沙沉積旋回層。從圖1A壩前沉積泥沙剖面中沙粒含量百分比隨深度變化的曲線來看，沙粒在壩前沉積泥沙垂直剖面上呈現高值低值交替出現的周期性變化規律，具體表現為：粗細相間分布，由上到下依次呈現細、粗、細、粗的分布規律，并且粗細分布的厚度差異較大。這主要是由于歷次侵蝕性降雨的水動力條件分選沉積作用而造成的，大水大沙下，淤積層厚度大，差異也大，水沙較小時，淤積厚度和差異減小。且由于降雨強度不同，流域內的侵蝕產沙量也不盡相同，進而導致剖面各沉積旋回層厚度的有所差異。

圖1 淤地壩壩前沉積泥沙沙粒含量變化特征(A)及沉積旋回厚度變化特征(B)

結合研究小流域降雨特性，根據壩前泥沙粗顆粒(沙粒)的周期變化特征[22]，從沙粒百分含量最小開始到沙粒百分含量逐漸增大到一定程度便開始突然減少為止即為一個完整的沉積旋回。因此，在本研究小流域淤地壩淤積期間，壩內沉積泥沙共存在19個沉積旋回層(圖1B)。旋回層的沉積厚度表示單次洪水事件過后，沉積在淤地壩前的泥沙總厚度。在該研究中，不同旋回層的沉積厚度差異顯著，平均厚度為12.7 cm，其中，第3旋回層厚度最小(1.5 cm)，第17旋回層厚度最大值(37 cm)。不同旋回層的沉積厚度標準差為10.48，變異系數為0.82，屬于中等強度變異。從沉積厚度變化趨勢分析發現，泥沙沉積較厚的層次均分布在建壩初期，隨著沉積時間，除去第1、第4和第5旋回層，其他旋回層厚度均低于10 cm，即壩前沉積泥沙層厚度隨著沉積時間顯著降低。

2.2 小流域侵蝕沉積泥沙粒徑特征分析

2.2.1 不同旋回土壤粒徑的變化特征 由表1可知，壩地剖面沉積旋回層中各粒徑組成以粗粉粒含量最多，所占比例達到51.73%；其次為黏粒，所占比例為19.88%；粗沙粒含量最少，所占比例僅僅為0.25%。研究小流域壩地剖面各粒徑含量由大到小分別是粗粉粒(51.73%)、黏粒(19.88%)、細粉粒(17.30%)、細沙粒(10.85%)、粗沙粒(0.25%)；各粒徑變異系數由大到小依次為粗沙粒>細沙粒>細粉粒>黏粒>粗粉粒，其中粗砂粒變異系數可達387.04%，而黏粒、細粉粒和粗粉粒僅僅分別為28.93%，36.16%和13.36%，說明粉粒與黏粒在研究小流域的輸移、沉積過程中分選性較差。

表1 研究區壩地剖面沉積旋回層各粒徑的統計分析

進一步從研究區壩地剖面不同沉積旋回層各粒徑的統計分析(圖2)可以看出，整個剖面共形成19個沉積旋回層，且每個旋回層中各級沉積泥沙在壩內的空間分布變化較大，變化差異顯著(p<0.05)。各沉積層粒徑組成均以粗粉粒為主，所占比例在42.85%～62.60%；不同旋回層中，粗沙粒僅在3個旋回層中出現(第1，9,12旋回層)；細沙粒(3.86%～27.02%)、細粉粒(8.20%～24.22%)和黏粒(12.05%～27.01%)的體積分數接近。而旋回層中粉粒與沙粒極值的空間分布則大致相反，即粉粒含量出現最大值時沙粒含量通常出現最小值，反之亦然。分析表明，退耕還林具有降低土壤細顆粒物質流失的作用，且年際變化顯著。

圖2 不同沉積旋回層粒徑的組成特征

2.2.2 不同旋回土壤粒徑分形維數的變化特征 由表2可以看出，該流域壩地剖面各沉積旋回層的土壤分形維數在2.60～2.78，平均值為2.70。其中第8旋回層分形維數最低(2.60)；第14旋回層最高(2.78)；分形維數的線性回歸的決定系數在0.91～0.96，相關性極好，說明所有的分形參數在計算上是可信的。但是不同旋回層次的分形維數隨著沉積深度有逐漸增大的趨勢，即隨著沉積時間不同旋回層次沉積泥沙顆粒的分形維數呈現增大趨勢，表明隨著退耕時間的增加，小流域侵蝕泥沙的細顆粒物質含量增加。但整個剖面上泥沙粒徑的分形維數隨沉積深度并無顯著性差異，這表明該流域土壤質量良好。

表2 小流域壩地剖面泥沙粒徑的分形維數的統計分析

為研究分形維數與土壤各粒級含量之間的關系，本文采用皮爾遜相關分析進行分析(表3)，分析結果表明，黏粒與分形維數的相關系數最大為0.99；粗砂粒的相關系數最小為-0.12。分形維數與黏粒、細粉粒、粗粉粒和細砂粒含量在99%檢驗水平上均呈現出極顯著相關關系，而與粗砂粒不存在顯著的相關關系；且黏粒和細粉粒與分形維數呈正相關關系，而粗粉粒和細砂粒與分形維數呈負相關關系。因此，具有較高的黏粒和細粉粒含量(即細物質顆粒含量)，較低的砂粒含量的旋回層具有較高的分形維數值。

表3 小流域壩地剖面泥沙粒徑的分形維數與各級粒徑相關分析

2.2.3 壩地土壤顆粒剖面粗化度的變化 壩地土壤顆粒組成主要以粗粉粒和細沙粒為主，說明徑流對小流域侵蝕泥沙顆粒具有一定的分選性。由于土壤侵蝕會引起粒徑的分選，因此采用粗化度可以用來說明沉積泥沙顆粒組成的變化。土壤質地粗化度指的是壩地剖面沉積旋回層中土壤的沙粒含量(>0.05 mm)與黏粒含量和粉粒含量之和(<0.05 mm)的比值。粗化度越大，則粗顆粒泥沙在土壤粒級組成中的成分越多；粗化度越小，則土壤粒級組成中的粗顆粒含量越少，土壤受侵蝕程度越大。

研究結果表明(圖3)，該研究區小流域壩地沉積泥沙質地粗化度分布范圍在0.04～0.37，平均值為0.16。且該壩地接近地表的剖面的質地粗化度均較大，而靠近最底下基巖位置處的質地粗化度都較小，即其泥沙顆粒粒徑組成呈現先細后變粗隨后又再次細化的趨勢。

圖3 壩地沉積泥沙質地粗化度的剖面變化

3 討 論

在黃土丘陵區，每場降雨或洪水都會將壩控小流域溝坡、溝谷及溝間地的土壤沖刷后，經過搬運、沉積，最后淤積在壩地。經過長期的研究與分析，進入流域出口的泥沙在淤積的過程中有一個不變的定律：即粗顆粒泥沙較黏粒和粉粒先行沉積，淤積在旋回層最下層，其次是粉粒和黏粒，同時，在淤地壩垂直剖面形成一個個具有時標功能的沉積旋回層[4]。而降雨量和降雨強度都會引起徑流量和水流流向的變化[23]。當流域內的水動力較強時，其沖刷和搬運能力也強，則更多的泥沙顆粒通過徑流被輸移至淤地壩進而發生沉積，沉積的泥沙粒徑較粗，淤積層厚度大；水動力條件較弱時，則小流域壩地沉積的泥沙顆粒粒徑較細，淤積層厚度相對較小，因此，每個沉積旋回層有明顯的粗細分層結構。且由于降雨強度不同，流域內的侵蝕產沙量也不盡相同，進而導致剖面各沉積旋回層厚度的有所差異且各沉積旋回層厚度的差異顯著。但隨著小流域水土流失治理林草措施(退耕還林草)的實施，小流域內土地利用和植被覆蓋逐年發生變化，研究表明退耕還林還草工程實施對土壤侵蝕具有明顯的減緩作用[24]，且在相同的植被蓋度等級下，耕地—林地轉化區和耕地—草地轉化區的土壤侵蝕模數呈減少趨勢[25]。本研究結果發現壩前沉積泥沙層厚度隨著沉積時間顯著降低，一方面是由于隨著植被的恢復，其對泥沙的攔截作用增強，導致進入溝道的泥沙量減少。另一方面是隨著泥沙的不斷淤積，原始溝道被泥沙所掩埋，在一定程度上減緩了溝道侵蝕。

淤積在壩前的泥沙經過了分選，相應的泥沙顆粒的組成也發生了變化，因此小流域壩地沉積物中泥沙顆粒粒徑組成對土壤侵蝕過程的演變具有一定的指示作用，本文的研究結果表明，該地區小流域淤地壩垂直剖面層中沉積泥沙不同粒徑的變異程度互不相同，變異性由大到小分別為粗沙粒(387.04%)>細沙粒(66.95%)>細粉粒(36.16%)>黏粒(28.93%)>粗粉粒(13.36%)。分析認為，沙粒和黏粒的變異性大小可能與其含量的大小有重要的關系。大于0.05 mm的粗顆粒泥沙的變異性與降雨、洪水特性、泥沙的來源，以及壩地所處的流域空間位置等因素密切相關，而由于黏粒的含量最大，因而變異性最小。當某次降雨強度較大時，徑流對表層土壤的沖刷能力增強，從溝邊與溝底進入流域內的泥沙顆粒特別是粗顆粒泥沙含量便會增多。造成壩地該次沉積泥沙粒徑組成中粗顆粒泥沙含量顯著增多，這一點與李勉等[8]的研究結果相符。

從粗粉粒含量隨剖面變化的曲線可以看出，深度越大，粗顆粒泥沙粒徑的變異程度越大，究其原因是由于淤地壩在建成之時，由于流域內河道比降較大，水動力條件較為復雜，進而導致沉積物泥沙顆粒粒徑的變異程度較大。隨著時間的推移，由于該淤地壩淤積的泥沙量和淤積面積的逐漸增大，致使流域內的水動力條件逐漸變得平穩，泥沙粒徑的變異程度呈顯著減小。且該流域壩地剖面各沉積旋回層的土壤分形維數在2.60～2.78，整個剖面上泥沙粒徑的分形維數隨有逐漸增大的趨勢，即隨著沉積時間不同旋回層次沉積泥沙顆粒的分形維數呈現增大趨勢，表明隨著退耕時間的增加，小流域侵蝕泥沙的細顆粒物質含量和分形維數均增大，這與董莉麗等[26]的研究結論一致。

進一步分析其質地粗化度，結果發現隨著剖面深度增大，其值變小，泥沙粒徑組成逐步細化，可能與周邊人為營造的樹木有關，因為植被不僅可以有效的防止水土流失，還可以攔截地表徑流，減小徑流流速，削減洪峰，致使沖刷過程中產生的泥沙減少，植被根系過濾了較粗的顆粒，致使進入流域內的粗顆粒泥沙量減少。

4 結 論

(1) 小流域壩地沉積泥沙剖面粒徑呈現高值低值交替出現的周期性變化規律：即粗細相間分布，由上到下依次呈現細、粗、細、粗的分布規律，且粗細分布的厚度差異較大，總體上沉積泥沙層厚度隨著沉積時間顯著降低。

(2) 淤地壩沉積泥沙粒徑主要以粗粉粒(0.01～0.05 mm)為主，平均百分含量為51.73%，且壩地剖面各分層的分形維數分布范圍在2.60～2.78，其中黏粒和細粉粒與分形維數呈極顯著正相關關系，而粗粉粒和細砂粒與分形維數呈極顯著負相關關系。隨著退耕時間的增加，小流域侵蝕泥沙中的細顆粒物質含量增加。

(3) 小流域壩地沉積泥沙質地粗化度在0.04～0.37，平均值為0.16。在接近地表的剖面的質地粗化度都較大，而靠近最底下基巖位置處的質地粗化度都較小。