江磊， 秦富倉， 李龍， 王露， 趙金洋

（內(nèi)蒙古農(nóng)業(yè)大學(xué)沙漠治理學(xué)院，呼和浩特 100018）

土壤是由不同顆粒組成、具有不規(guī)則形狀自相似結(jié)構(gòu)的多孔介質(zhì)，并且具有一定的分形特征[1]。Tyler等[2]最早將分形理論引入土壤學(xué)研究；楊培嶺等[3]、王國梁等[4]分別用土壤顆粒重量分布和體積分布（土壤體積分形維數(shù)）來表征土壤顆粒的分形維數(shù)，提出土壤分形維數(shù)可反映土壤粒徑大小和土壤質(zhì)地的均一程度。土壤分形維數(shù)越高則土壤越緊實，反之則松散、通透性越好。土壤的分形特征不僅可以反映土壤質(zhì)地和均一程度，還越來越多地被用于定量描述土壤肥力特征、土壤侵蝕和入滲性能。肖東東等[5]通過分析不同林地土壤分形特征及其與土壤理化性質(zhì)的相關(guān)性，以分形維數(shù)值來表征土壤質(zhì)地和養(yǎng)分狀況，進而提出改良和優(yōu)化植被生長的建議。茹豪等[6]通過分析黃土丘陵溝壑區(qū)不同土地利用類型對土壤粒徑分形參數(shù)的影響，預(yù)測不同土地利用類型的侵蝕程度。徐永福等[7]通過膨脹土的土壤粒徑分布確定土壤表面分形維數(shù)，估算了土壤的水分特征曲線。戰(zhàn)海霞等[8]通過分析沂蒙山不同植物群落土壤顆粒的分形特征，提出植被是通過增加土壤粉粒、黏粒來改良土壤分形結(jié)構(gòu)，進而改善土壤團粒結(jié)構(gòu)和孔隙結(jié)構(gòu)，提高入滲性能，并得出土壤分形維數(shù)與入滲率呈顯著正相關(guān)的結(jié)論。

(1)礦級響應(yīng)處置內(nèi)容，即作業(yè)條件變化大、施工難度大、技術(shù)性強、或涉及面廣、對全礦的安全生產(chǎn)影響大的一級變化內(nèi)容。

石峽口水庫位于黃土丘陵區(qū)，粉砂土占比較大，且土壤侵蝕嚴重，是否能通過人工林增加土壤黏粉粒、改變其分形特征、提高入滲能力等均有待探討。土壤入滲過程指降水或灌溉后地面上的水滲入土壤形成土壤水的過程，土壤的入滲特性反映其對植被水分的補給能力，對雨水、徑流的抗蝕能力[9]。自然狀況下，不同植被的土壤入滲特性受到土壤容重、孔隙度、質(zhì)地、初始含水率、地表結(jié)皮、根系分布、動物巢穴等因素的影響，其中重要因素是土壤粒徑分布（particle size distribution，PSD)。土壤體積分形維數(shù)根據(jù)土壤PSD計算而來，反映了土壤PSD的特征[4]。根據(jù)土壤體積分形維數(shù)與土壤理化性質(zhì)的關(guān)系、入滲率與土壤體積分形維數(shù)及土壤理化性質(zhì)的關(guān)系可以定量描述土壤體積分形維數(shù)與土壤水分入滲的關(guān)系，研究結(jié)果可以直接用土壤體積分形維數(shù)來比較研究區(qū)不同林地表層土壤入滲規(guī)律。

石峽口水庫位于清水河縣韭菜莊鄉(xiāng)，設(shè)計效益以蓄水、防洪、灌溉、城鎮(zhèn)供水為主，是黃土丘陵溝壑區(qū)就地攔蓄降水的實例。其所在的黃土丘陵溝壑區(qū)是全世界土壤侵蝕和水土流失最嚴重的地區(qū)，水庫周邊的人工林發(fā)揮著保水固土、涵養(yǎng)水源的功能。入滲能力較強的土壤通過對降水的吸收、滲透及調(diào)節(jié)作用，將其變?yōu)槿乐辛骱偷叵铝?，對該地區(qū)雨水侵蝕和徑流的減弱效果良好[10]。由于長期的粗放管理，水庫周邊的人工林生長情況及林下土壤養(yǎng)分、涵養(yǎng)水源能力及土壤侵蝕狀況差異顯著，水庫長期出現(xiàn)蓄水不足等問題。因此，本研究采用環(huán)刀進行多點采樣對該地區(qū)3種典型人工林表層土壤（0—20 cm）的理化性質(zhì)進行測定，在中國制粒徑分級標準[11]下采用粒徑體積分布計算土壤體積分形維數(shù)，用透明亞克力管填土對不同人工林表層土壤進行入滲分析，以期了解研究區(qū)不同人工林的土壤結(jié)構(gòu)和養(yǎng)分狀況，探索不同人工林對土壤分形特征及入滲性能的影響，揭示土壤體積分形維數(shù)與土壤理化性質(zhì)及入滲的關(guān)系，為黃土丘陵區(qū)植被恢復(fù)和水源涵養(yǎng)林建設(shè)提供參考。

1 研究方法

1.1 研究區(qū)概況

石峽口水庫位于內(nèi)蒙古自治區(qū)清水河縣韭菜莊鄉(xiāng)（39°54′—39°55′N，111°54′—111°57′E），海拔1 400 m左右，屬于內(nèi)蒙古高原與黃土高原交接處，是典型的黃土丘陵區(qū)（圖1）。年平均降水量413.8 mm。黃土覆蓋較厚，土壤類型主要為栗褐土。水庫于1977年建成，總庫容1.38×108m3，水庫狹長約4 km，最寬處250 m，高66 m。由于水庫周邊植被稀疏，近年來以油松（Pinus tabuliformisCarr.）、落葉松〔Larix gmelinii（Rupr.）Kuzen.〕、檸條（Caragana korshinskiiKom.）、沙棘（Hippophae rhamnoidesLinn.）為主的樹種進行了密植。

圖1 研究區(qū)地理位置Fig.1 Geographical location of the study area

1.2 樣地選取

分別選取水庫周圍山頂?shù)挠退?、落葉松、檸條×沙棘混交林及草地。3種人工林立地類型一致、海拔高度均在1 400 m以上，林齡為13～14 a，草地在3種林地附近。在每種植被中選取3片長勢良好、近似正方形的樣地，共計12片樣地（表1）。

從土壤粒徑頻度分布（圖2）可以看出，細黏粒(0～1 μm)、粗黏粒(1～2 μm) 在3種林地及對照草地表層土壤的累計含量均極少；細粉粒（2～5 μm）、中粉粒（5～10 μm）在3種林地及對照草地表層土壤的累計含量均在4%～7%左右。3種典型人工林及草地土壤的黏粒、細粉粒、中粉粒差異不顯著。從油松到草地，土壤中粗粉粒（10～50 μm）、細砂粒（50～200 μm）的占比減少，而粗砂粒（200～1 000 μm）及石礫（1 000～3 000 μm）的占比增加。3種人工林及草地土壤均為粗粉粒（10～50 μm）累計含量最多，且不同類型間差異顯著，表現(xiàn)為油松（47%）＞落葉松（41%）＞檸條×沙棘（22.83%），草地土壤中石礫含量最多。總體上看，研究區(qū)土壤中粗粉粒（10～50 μm）累計含量最多，其次為細砂粒（50～200 μm），分別占土壤顆粒體積的38%和27%。

表1 樣地概況Table 1 Sample plot overview

1.3 采樣及測定方法

采樣時間為2021年7月。由于人工林種植較均勻，采用5點取樣法對油松、落葉松純林、檸條×沙棘混交林和草地分別進行取樣，樣方大小分別為20 m×20 m、5 m×5 m、1 m×1 m。將表層雜草、枯落物除去后，挖土壤剖面(0.5 m深×0.5 m寬)。采用環(huán)刀、鋁盒、塑封袋每5 cm（用于平均表層土壤理化性質(zhì)及入滲土柱的分層填土）采集一層土樣用于土壤理化性質(zhì)測定。環(huán)刀鋁盒法計算土壤容重、孔隙度、含水率。將塑封袋內(nèi)土壤晾干過不同網(wǎng)篩后，采用馬爾文3 000激光粒度分析儀測定土壤粒徑分布，并計算土壤體積分形維數(shù)(Dv)[4]。采用堿解擴散法測定有效氮含量；采用碳酸氫鈉侵提-鉬銻比色法測定速效磷含量；采用乙酸氨浸提-火焰光度法測定速效鉀含量；采用重鉻酸鉀氧化-油浴加熱法測定測有機質(zhì)含量[12]。

式中：k=p2μ0τl/g是與極對數(shù)、氣隙磁導(dǎo)率及電機幾何尺寸有關(guān)的常數(shù)。式(4)表明，電磁轉(zhuǎn)矩與定、轉(zhuǎn)子磁動勢的幅值Fs和Fr之乘積成正比，且是定、轉(zhuǎn)子磁動勢矢量之間交角λ的正弦函數(shù)。

式中，V為每次裝入土體體積(cm3)；ρb為土壤容重(g·cm-3)；θv為室內(nèi)土壤含水率(%)。

傳統(tǒng)觀點認為，金融中介和資本市場這兩個金融系統(tǒng)的不同組成部分對實體經(jīng)濟的影響是相互獨立的。而它們在影響技術(shù)進步時，由于自身的特點導(dǎo)致其在面對創(chuàng)新時的反應(yīng)大相徑庭。

1.4 數(shù)據(jù)處理

采用Excel 2010和SPSS 23.0對數(shù)據(jù)進行處理、方差分析、相關(guān)性分析及回歸分析，并計算下列參數(shù)。

干堆法中填土量(W) 的計算公式如下。

5.選人用人的導(dǎo)向性。競爭的實質(zhì)是參與干部德能勤績綜合素質(zhì)的比較。在競爭性選拔過程中，干部的升降去留完全以本人的德才、政績和民意為依據(jù)。這就教育、引導(dǎo)和激勵干部努力學(xué)習(xí)、勤奮工作、廉潔自律，以自己的真才實學(xué)取勝，從而形成正確的用人導(dǎo)向，營造良好的競爭環(huán)境。

式中，Ri表示相鄰2個粒級的平均粒徑（μm）；Rmax為最大粒級的平均粒徑（μm）；V(r＜Ri)為粒徑小于Ri的累積土壤顆粒體積百分含量（%）；VT為土壤各粒級總體積百分含量（%）。

土壤體積分形維數(shù)(Dv)的計算公式如下。

為研究土壤PSD對入滲的影響，要防止根系、動物巢穴和石塊的干擾，由此采用單環(huán)土柱法進行填土入滲。土柱為圓柱形透明的亞克力管，直徑15 cm，高40 cm，中間填裝表層0—20 cm土壤，上下10 cm分別為入水層和砂石層。填土?xí)r將下層填10 cm的砂石并放置鐵絲網(wǎng)，按照不同深度的土壤容重、室內(nèi)土壤含水率、填土高度，采用干堆法計算土柱每層應(yīng)填的土壤重量[13]。在之前挖土壤剖面的原地點，除掉表層植被后，一層一層的篩土稱重，每層5 cm，并夯實，層間打毛，保證最后填入的土壤剛好在5、10、15 cm刻度線，依次類推，盡可能保證所填土狀態(tài)接近原土狀態(tài)。填完后，在土柱入水層放入鐵絲網(wǎng)和孔板，防止水的擊濺和土壤運動；并將土柱邊緣涂上凡士林以防止形成大孔隙流。馬氏瓶直徑15 cm，高40 cm，保證連續(xù)供水，中途不加入水。

采用Kostiakov公式[14]對入滲過程進行計算及模型。

式中，f為入滲率(mm·min-1)；a為常數(shù)；t為時間（min）。

2 結(jié)果與分析

2.1 表層土壤的理化性質(zhì)比較分析

土壤理化性質(zhì)是研究土壤分形特征和水分入滲的前提，3種人工林及對照草地的表層土壤理化性質(zhì)如表2所示。3種人工林中，檸條×沙棘混交林土壤容重最大，而土壤孔隙度（指土壤總孔隙度）最??；油松和落葉松林間的表層土壤容重和土壤孔隙度差異不顯著。3種典型人工林的表層土壤自然含水率均在5%左右，且差異不顯著，草地土壤的自然含水率最低。落葉松、檸條×沙棘混交林土壤的有效氮含量顯著高于油松、草地，在土壤養(yǎng)分分級[15]中處于極高水平，而油松林、草地土壤的有效氮含量屬于中下水平。土壤速效磷含量表現(xiàn)為油松＞落葉松＞檸條×沙棘混交林、草地，但是在土壤養(yǎng)分分級中均處于低和較低水平。速效鉀含量在3種典型人工林間差異不顯著，草地土壤的速效鉀含量最低。落葉松林下土壤的有機質(zhì)含量最高，油松和檸條×沙棘混交林間差異不顯著，草地土壤的有機質(zhì)含量最低。

上海市崇明區(qū)2016至2017年肺結(jié)核發(fā)病率從50歲開始出現(xiàn)高峰，早于2008至2009年(60歲)[7]，這可能與上海中老年人群重視體檢并及時診療有關(guān)。2010年，全國第5次結(jié)核病流行病學(xué)抽樣調(diào)查報告[8]發(fā)現(xiàn)，老年肺結(jié)核患者中近40%無癥狀，有癥狀者未就診的比例也超過1/2。老年人群因年齡大、生理功能減退、抗病能力低下、愛集聚娛樂等原因，極易感染結(jié)核桿菌；老年肺結(jié)核發(fā)病率高也可能與其?；加刑悄虿?、心血管疾病及肺部其他疾病，導(dǎo)致機體抵抗力下降有關(guān)[9]。

圖2 土壤粒徑頻率分布Fig.2 Frequency distribution of soil particle size

2.2 土壤體積分形維數(shù)與土壤理化性質(zhì)的關(guān)系

土壤體積分形維數(shù)是反映土壤顆粒分布的重要指標，指土壤顆粒的排列和組合形式[15]。3種林地表層土壤體積分形維數(shù)均值如圖3所示，表現(xiàn)為草地（2.30）、檸條×沙棘混交林（2.28）＞落葉松（2.24）＞油松（2.18）。對土壤體積分形維數(shù)與粒徑分布及其他土壤理化性質(zhì)進行相關(guān)分析，結(jié)果（表3）表明，研究區(qū)土壤體積分形維數(shù)與細黏粒、粗黏粒相關(guān)不顯著，與細粉粒、中粉粒呈極顯著正相關(guān)，與粗粉粒、細砂粒呈顯著負相關(guān)；土壤體積分形維數(shù)與土壤容重呈正相關(guān)，與土壤孔隙度呈顯著負相關(guān)，與自然含水率相關(guān)不顯著。由此表明，研究區(qū)土壤中細粉粒、中粉粒、粗砂粒含量越多，則土壤體積分形維數(shù)越大；土壤容重越大，孔隙度越小，土壤體積分形維數(shù)越大。從土壤體積分形維數(shù)與養(yǎng)分關(guān)系來看，土壤體積分形維數(shù)與土壤有效氮含量呈正相關(guān)，與土壤速效磷含量呈顯著負相關(guān)，與有機質(zhì)、速效鉀含量相關(guān)不顯著。根據(jù)粒徑分布及差異來看，粗粉粒和細砂粒含量對研究區(qū)土壤體積分形維數(shù)有重要影響。

國內(nèi)外處理鉆井廢棄泥漿方法主要有：固液分離法、環(huán)空或地層注入法、土地分散耕施法、干燥焚燒法、化學(xué)固化法、生物降解法、轉(zhuǎn)化為其它用途如制磚等［3］，但是這些方法都存在成本高、應(yīng)用性差或處理量小的特點，不適合低成本大批量處理。

表3 土壤體積分形維數(shù)與土壤理化性質(zhì)的關(guān)系Table 3 Relationship between soil volume fractal dimension and soil physical and chemical properties

圖3 土壤體積分形維數(shù)Fig.3 Fractal dimension of soil volume

2.3 土壤的入滲特征分析

通過Kostiakov公式進行擬合(圖4)，3種人工林入滲擬合方程分別為：落葉松y=11.382X-0.58611（R2=0.97），油松y2=13.623X2-0.56（6R2=0.98），檸條×沙棘y=7.846X-0.54（1R2=0.96），草地y=6.643X-0.6943344（R2=0.95），相關(guān)系數(shù)R2均大于0.9，說明方程擬合效果極好。由于試驗誤差及水頭控制等因素，將2 min時的入滲率記為初滲速率，則初滲速率在不同類型間表現(xiàn)為油松（11.0 mm·min-1）＞落葉松（8.5 mm·min-1）＞檸條×沙棘混交林（6.2 mm·min-1）＞草地（4.7 mm·min-1）；穩(wěn)滲速率表現(xiàn)為油松（1.2 mm·min-1）＞落葉松（0.9 mm·min-1）＞檸條×沙棘混交林（0.8 mm·min-1）＞草地（0.45 mm·min-1）；1 h的累計入滲量表現(xiàn)為油松（116 mm）＞落葉松（80 mm）＞檸條×沙棘混交林（62 mm）＞草地（51 mm）。由此可見，3種入滲因子都表現(xiàn)為油松最大，表明油松表層土壤的入滲性能最強。

圖4 土壤水分入滲特征Fig.4 Soil water infiltration

2.4 入滲因子與土壤體積分形維數(shù)、物理性質(zhì)的關(guān)系

對土壤分形維數(shù)（變量X）與土壤入滲因子（應(yīng)變量y）進行回歸分析，結(jié)果（圖5）表明，土壤體積分形維數(shù)與初滲速率、穩(wěn)滲速率和1 h累計入滲量均呈顯著負相關(guān)。初滲速率與土壤體積分形維數(shù)的回歸方程為：y1=-35.392X1+87.407（R2=0.659 5，P=0.002）；穩(wěn)滲速率與土壤體積分形維數(shù)關(guān)系：y2=-3.915 6X2+9.638 6（R2=0.655 1，P=0.002）；1 h累計入滲量與土壤體積分形維數(shù)關(guān)系：y3=-360.32X3＋892.71（R2=0.726 3，P=0.001）。由此表明，研究區(qū)土壤體積分形維數(shù)越大，土壤入滲性能越強。將土壤入滲因子與土壤PSD、容重、孔隙度做相關(guān)分析，結(jié)果（表4）表明，初滲速率、穩(wěn)滲速率、1 h累計入滲量均與細粉粒、粗粉粒、細砂粒、粗砂粒呈一定相關(guān)性，且與土壤容重呈顯著負相關(guān)，與土壤孔隙呈顯著正相關(guān)。

表4 入滲特征與土壤物理性質(zhì)的關(guān)系Table 4 Relationship between infiltration characteristics and soil physical properties

圖5 土壤入滲特征與土壤體積分形維數(shù)的關(guān)系Fig.5 Relationship between infiltration characteristics and soil volume fractal dimension

3 討論

研究區(qū)3種人工林中，油松和落葉松林地的土壤容重差異不顯著，因而兩者對比入滲率時可能受土壤PSD影響較大。落葉松林地土壤的有機質(zhì)、速效鉀、有效氮含量明顯較高，說明其養(yǎng)分含量較油松、檸條×沙棘混交林高；檸條×沙棘混交林林地表層土壤的容重、孔隙度最小，說明檸條×沙棘混交林林地土壤緊實，通氣性差；除其他養(yǎng)分外，檸條×沙棘混交林林地土壤的有效氮含量較高，可能是由于豆科植物較強的固氮能力。但研究區(qū)3種典型人工林土壤的速效磷含量都屬于極低水平，不利于人工林生長，迫切需要進行人工干預(yù)來補充磷元素。對比草地土壤的理化性質(zhì)，說明植樹造林明顯改善了土壤結(jié)構(gòu)和養(yǎng)分狀況。土壤PSD是土壤物理性質(zhì)之一，與土壤通氣、保肥、保水關(guān)系密切[16]。對不同粒度進行分析表明，研究區(qū)的土壤屬于粉砂質(zhì)土，3種人工林及草地土壤中的黏粒含量都較少，且差異不顯著，其中油松、落葉松林地的粗粉粒含量最多，而檸條×沙棘混交林的粗砂粒含量最多，草地的石礫含量最多。因此，通過植被增加粗粉粒、細砂粒含量可能不科學(xué)，只能通過植物提高對土壤細顆粒的保護能力，一方面是由于油松、落葉松的枯落物較多，大部分降雨被枯落物儲存；另一方面是因為喬木比灌木林有更多的樹干截留和樹干徑流，減緩了地面徑流速度，因此枯落物下方較細的土壤顆粒被保存起來，而檸條×沙棘混交林表層土壤中的細顆粒易被沖刷，留下較粗的石礫；油松與落葉松相比，油松土壤細顆粒更多，可能是由于油松為常綠樹種，在冬春季節(jié)的北方山區(qū)，抗風(fēng)蝕能力和對表層土壤保水能力優(yōu)于落葉松；而草地的表層土壤大部分直接裸露在地表，因此細顆粒含量明顯減少，留下裸露的石礫。

研究區(qū)3種典型人工林及草地的土壤體積分形維數(shù)表現(xiàn)為：草地（2.30）、檸條×沙棘混交林（2.28）＞落葉松（2.24）＞油松（2.18）。引起不同人工林土壤體積分形維數(shù)差異的原因主要是植被對土壤粒徑分布的改變。本研究表明，土壤體積分形維數(shù)與黏粒（中國粒徑分級）無顯著相關(guān)性，這可能與研究區(qū)黏粒含量較少及3種人工林土壤黏粒差異不顯著有關(guān)。土壤體積分形維數(shù)與粉粒的關(guān)系存在爭議[17-19]，本研究表明土壤體積分形維數(shù)與細粉粒、中粉粒呈顯著正相關(guān)，而與粗粉粒呈顯著負相關(guān)；土壤體積分形維數(shù)與細砂粒呈顯著負相關(guān)，而與粗砂粒呈顯著正相關(guān)。曾憲勤等[20]研究表明，當(dāng)黏粒含量較少或一定時，最大粒級越大則分形維數(shù)越大。而檸條×沙棘混交林、草地含有石礫較多，因此其分形維數(shù)較大。研究區(qū)土壤體積分形維數(shù)與土壤容重也呈顯著正相關(guān)，即土壤體積分形維數(shù)越大，土壤越緊實，孔隙度變小。土壤體積分形維數(shù)與初始含水率無顯著相關(guān)性，可能是因為3種林地的初始含水率差異不顯著。土壤體積分形維數(shù)與土壤養(yǎng)分存在一定相關(guān)性，與有效氮含量呈顯著正相關(guān)，與速效磷含量呈顯著負相關(guān)，而與其他養(yǎng)分含量相關(guān)不顯著，表明在研究區(qū)內(nèi)土壤體積分形維數(shù)只能表征土壤有效氮和速效磷含量的多少，難以表征土壤整體的養(yǎng)分狀況。

3種典型人工林及草地土壤的入滲能力表現(xiàn)為油松＞落葉松＞檸條×沙棘＞草地，說明人工林明顯改善了土壤入滲能力，且油松林地土壤的入滲能力最強，能夠?qū)⒏嗟挠晁D(zhuǎn)化為壤中流和地下徑流，有效減少水庫邊坡徑流和侵蝕。土壤體積分形維數(shù)與初滲速率、穩(wěn)滲速率、1 h累計入滲量均呈顯著負相關(guān)，與戰(zhàn)海霞等[8]結(jié)果存在差異，這可能是由于不同研究區(qū)的土壤類型存在差異造成。由入滲因子與土壤體積分形維數(shù)的關(guān)系可知，研究區(qū)土壤體積分形維數(shù)越小入滲率越強，累計入滲量越多，主要受粗粉粒和細砂粒影響，且容重是影響土壤入滲的主要因素。油松與落葉松林地表層土壤的容重差異不顯著，但兩者粒徑分布及土壤體積分形維數(shù)差異顯著，說明兩者入滲率受粗粉粒和細砂粒影響較大，并與粗粉粒、細砂粒呈顯著正相關(guān)。徐萍等[21]研究表明。在黏粒較少的情況下，土壤細顆粒越多，土壤顆粒及孔隙分布的均勻度越高，有利于增加土壤入滲速率和持水能力。在本研究區(qū)范圍內(nèi)，土壤體積分形維數(shù)可以表征土壤水分入滲性能，土壤分形維數(shù)越小入滲越強。未來在分析不同植被的土壤分形維數(shù)對入滲能力影響時，應(yīng)選取更多樹種，且在相似的土壤容重和孔隙度條件下，這樣可以排除其他因子的干擾。