李青峰，邢旭光，馬孝義

(1.西北農(nóng)林科技大學(xué)水利與建筑工程學(xué)院， 陜西 楊凌 712100； 2.西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)實(shí)驗(yàn)室， 陜西 楊凌 712100)

土壤持水性能評(píng)價(jià)對(duì)制定合理的農(nóng)田灌溉制度具有指導(dǎo)意義，土壤持水與保水能力可以通過土壤孔隙大小和分布特征進(jìn)行判定，然而土壤的持水性、土壤水分有效性以及孔隙分布特性均可以通過土壤水分特征曲線進(jìn)行評(píng)價(jià)[1-2]。許多學(xué)者對(duì)土壤水分特征曲線擬合[3-5]及其影響因素[6-8]進(jìn)行了研究，同時(shí)也有較多研究基于土壤水分特征曲線對(duì)土壤持水能力以及土壤水分有效性進(jìn)行評(píng)價(jià)[9-13]，縱觀有關(guān)土壤水分特征曲線的研究成果，較多針對(duì)添加物條件或是將非鹽堿土作為研究對(duì)象，而對(duì)鹽堿土壤持水性能的研究仍較少，關(guān)于不同鹽分離子對(duì)其影響效果的評(píng)價(jià)則更是鮮有報(bào)道。明晰不同鹽離子對(duì)土壤持水性能的影響，將有利于對(duì)不同類型鹽堿地制定合理洗鹽方式。

土壤水分特征曲線一般是通過離心機(jī)測(cè)得，然而隨著離心機(jī)預(yù)設(shè)吸力和平衡時(shí)間逐漸增加，土體發(fā)生失水，并伴隨沉降，容重逐漸增大，最終導(dǎo)致土體表面出現(xiàn)裂縫。土壤干縮開裂問題對(duì)于農(nóng)田和作物均具有較強(qiáng)的破壞作用，土壤裂縫在很大程度上可導(dǎo)致水分與養(yǎng)分流失，造成灌溉效率降低和水環(huán)境污染，同時(shí)又可能破壞作物根系，導(dǎo)致產(chǎn)量降低[14-15]。目前有關(guān)土壤裂縫的研究多是集中于土壤總鹽的影響[16-18]，或是單獨(dú)將土壤干縮開裂作為研究對(duì)象[14,19]，這在一定程度上忽略了不同鹽離子各自產(chǎn)生的效應(yīng)，且較少研究將土壤失水和土體收縮相結(jié)合，即土壤脫水過程中的土體收縮變化特征。土體收縮變化與吸力、含水率和容重變化具有同步性，故在土壤發(fā)生失水時(shí)，對(duì)其土壤水分特征曲線和土壤收縮特征同時(shí)展開研究具有實(shí)際意義；對(duì)不同鹽離子導(dǎo)致的土壤裂縫進(jìn)行評(píng)價(jià)，可為進(jìn)一步認(rèn)識(shí)鹽堿土的干縮開裂機(jī)理和農(nóng)田灌溉提供指導(dǎo)。

將土壤失水與土壤收縮相結(jié)合，在已有研究基礎(chǔ)上，探索4種陽性及陰性鹽離子對(duì)土壤持水及收縮特性的作用效果，以期更好的為不同鹽堿土壤的持水性評(píng)價(jià)奠定基礎(chǔ)，同時(shí)為土壤失水過程中的干縮開裂研究提供參考。

1 材料與方法

1.1 供試材料與試驗(yàn)設(shè)計(jì)

試驗(yàn)土壤取自灌溉試驗(yàn)站小麥—玉米輪作試驗(yàn)田，采集深度為30 cm農(nóng)田耕作層；供試土壤經(jīng)風(fēng)干且過2 mm篩后，采用激光粒度儀測(cè)定土壤顆粒組成，粒徑介于0～0.002 mm、0.002～0.02 mm和0.02～2 mm的顆粒含量分別為17.28%、44.32%和38.40%，由國(guó)際制質(zhì)地三角形可知土壤類型為粉黏壤土。

根據(jù)田間實(shí)測(cè)土壤容重(1.38～1.40 g·cm-3)將裝填干容重設(shè)置為1.40 g·cm-3，均勻裝入環(huán)刀(100 cm3)內(nèi)；用蒸餾水分別與氯化鉀、氯化鈉、氯化鈣、氯化鎂、碳酸氫鈉、碳酸鈉和硫酸鈉粉末/晶體配置成鹽溶液(濃度均約為1 g·L-1)。試驗(yàn)開始前，分別采用各鹽溶液對(duì)土樣進(jìn)行飽和處理48 h，并以蒸餾水浸泡為對(duì)照(CK)；試驗(yàn)結(jié)束后，采用干燥法測(cè)定土壤含水率；各處理均4次重復(fù)。

1.2 試驗(yàn)測(cè)定方法

1.2.1 土壤持水曲線測(cè)定 采用離心機(jī)法測(cè)定飽和環(huán)刀樣品的土壤水分特征曲線，高速恒溫冷凍離心機(jī)(CR21G II型)內(nèi)恒溫4℃；測(cè)定吸力為10、50、100、300、500、700、1 000、3 000、5 000 cm和7 000 cm，各吸力達(dá)到平衡時(shí)間后，用電子天平稱取土樣質(zhì)量進(jìn)而計(jì)算土壤含水率。采用van Genuchten模型(式1)對(duì)供試土樣進(jìn)行持水曲線擬合，進(jìn)而基于土壤水分特征曲線對(duì)比分析各鹽離子對(duì)土壤持水能力的影響。

(1)

式中，θ為體積含水率，cm3·cm-3；θs為飽和體積含水率，cm3·cm-3；θr為殘余體積含水率，cm3·cm-3；s為吸力，cm；α為進(jìn)氣吸力的倒數(shù)(即α=1/sa)；m和n為形狀系數(shù)。

1.2.2 土體收縮特性測(cè)定 離心過程中環(huán)刀內(nèi)土體產(chǎn)生沉降現(xiàn)象，并以垂向收縮為主；每達(dá)到平衡時(shí)間時(shí)，用游標(biāo)卡尺量取各轉(zhuǎn)速下環(huán)刀內(nèi)土體的軸向沉降高度進(jìn)而計(jì)算土體收縮量，采用線縮率(式2)描述土體在各吸力時(shí)的收縮度，采用軸向收縮應(yīng)變(式3)表征離心結(jié)束時(shí)土體的最終收縮形變量；同時(shí)測(cè)定各吸力下的土體容重，進(jìn)而對(duì)離心過程中的土壤容重進(jìn)行校正。

δsl=zi/H×100%

(2)

δs=Δh/H×100%

(3)

式中，δsl為土體線縮率，%；δs為土體軸向收縮應(yīng)變，%；zi為各吸力對(duì)應(yīng)的土壤收縮度，mm；H為土樣初始高度，mm；Δh為土壤脫水始末狀態(tài)高度差，mm。

1.2.3 土體開裂特征測(cè)定 在試驗(yàn)結(jié)束后用數(shù)碼相機(jī)(Canon IXUS 950IS，分辨率3264×2448)拍攝獲取各處理環(huán)刀樣品截面圖像，將相機(jī)置于固定支架以保證拍攝高度相同，且拍攝時(shí)遮擋外界光源、關(guān)閉閃光燈，僅用日光燈照明，從而使得各處理土壤裂縫圖像具有一致性；而后編寫Matlab程序進(jìn)一步對(duì)圖像中的裂隙參數(shù)進(jìn)行提取。

在原有的數(shù)字圖像中，裂隙區(qū)域圖像較無裂隙區(qū)域的色彩更黑且純度更高，鑒于飽和度分量圖可以更明顯地展現(xiàn)土壤裂隙，因此，將原有的RGB圖像進(jìn)行HSI轉(zhuǎn)換，提取飽和度分量圖并采用全局閾值法得到土壤裂縫的二值化圖。在此基礎(chǔ)上，為了盡量保留原始的短小裂縫，利用膨脹腐蝕方法對(duì)短小裂縫進(jìn)行適當(dāng)連接，并去除該圖像中的孤立單點(diǎn)或雜點(diǎn)。最后根據(jù)處理后的土壤裂縫二值化圖像提取裂縫總長(zhǎng)度和裂縫面積參數(shù)，并進(jìn)一步計(jì)算獲取裂縫長(zhǎng)度密度和裂縫面積密度參數(shù)。

2 結(jié)果與分析

2.1 鹽基離子對(duì)土壤水分特征曲線的影響

土壤含水率隨著吸力增加呈現(xiàn)逐漸降低趨勢(shì)，從圖1(a)可以看出，在脫水過程中陽離子處理土壤含水率在較低吸力段(h=0～1 000 cm)和較高吸力段(h=1 000～7 000 cm)分別平均減小了約39.25%和18.96%；陽離子處理的土壤水分特征曲線與CK處理呈現(xiàn)平行狀，可見各處理土壤比水容量差異較小，土壤孔隙結(jié)構(gòu)差異不大。圖1(a)顯示，在h=0～7 000 cm內(nèi)當(dāng)吸力相同時(shí)，陽性鹽離子處理的土壤含水率大小表現(xiàn)為K+>Na+>CK>Ca2+>Mg2+，即不同陽性鹽離子對(duì)土壤持水能力的影響程度不同，具體表現(xiàn)為K+>Na+>Ca2+>Mg2+，且低濃度(≈1 g·L-1)K+和Na+可增強(qiáng)土壤持水能力，而Ca2+和Mg2+使得土壤持水能力減弱。從圖1(b)可以看出，陰離子處理土壤含水率在較低吸力段和較高吸力段分別平均減小了約30.76%和15.83%；Cl-和SO42-處理土壤水分特征曲線基本重合，土壤含水率平均絕對(duì)誤差僅約±0.003 cm3·cm-3，且與CK處理呈現(xiàn)平行狀，表明此2種離子對(duì)土壤持水能力的影響效果相似，且對(duì)土壤孔隙結(jié)構(gòu)的改變不明顯；而CO32-和HCO3-處理土壤水分特征曲線明顯發(fā)生上移現(xiàn)象，在一定程度上改變了土壤孔隙結(jié)構(gòu)，進(jìn)而使得土壤持水能力能到顯著提升(圖1(b)，表1)。圖1(b)顯示，在h=0～7 000 cm內(nèi)當(dāng)吸力相同時(shí)，陰性鹽離子處理的土壤含水率大小表現(xiàn)為CO32->HCO3->SO42->Cl->CK，即不同陰性鹽離子對(duì)土壤持水能力的影響程度不同，具體表現(xiàn)為CO32->HCO3->SO42->Cl-，且4種離子均可增強(qiáng)土壤持水能力。與CK相比，各處理土壤含水率在不同吸力段的增減程度不同，如表1所示，除Ca2+和Mg2+之外，其余鹽離子在一定程度上均能增強(qiáng)土壤持水能力，且在h=0～1 000 cm范圍內(nèi)效果不明顯，而在h=1 000～7 000 cm范圍內(nèi)土壤持水能力顯著增強(qiáng)，表明吸力增加有利于該6種離子提高土壤持水能力，且陰陽離子中分別以CO32-和K+改善土壤持水性的效果最為顯著，HCO3-對(duì)于提高土壤持水能力也具有較為顯著的作用。

圖1 不同鹽離子處理土壤水分特征曲線

注：表中數(shù)據(jù)為鹽離子處理土壤含水率比對(duì)照處理高或低的百分比；“+”表示比對(duì)照處理高，“-”表示比對(duì)照處理低。

Note: The figures in Table 1 indicate the percentage of soil moisture for salt ion treatments higher or lower than for CK. “+” indicates that the soil moisture for salt ion treatments is higher than for CK, and “-” indicates that the soil moisture for salt ion treatments is lower than for CK.

2.2 鹽基離子對(duì)土壤收縮特征曲線的影響

土壤比容積和質(zhì)量含水率關(guān)系曲線即為土壤收縮特征曲線，可用來根據(jù)土壤干濕變化預(yù)測(cè)田間地表收縮與膨脹、計(jì)算土壤含水率剖面、制定灌溉計(jì)劃[20]。對(duì)比圖2(a)、2(b)可以發(fā)現(xiàn)，陰陽鹽離子對(duì)土壤收縮特征曲線具有不同程度的影響：(1) 對(duì)于陽性鹽離子，4種離子處理土壤收縮曲線形態(tài)類似，且數(shù)值差異較小，土壤比容隨含水率的變化趨勢(shì)一致。由此可知，K+、Na+、Ca2+和Mg2+均可維持土壤收縮曲線的3個(gè)收縮段無顯著變化且一致性較高；(2) 對(duì)于陰性鹽離子，SO42-和Cl-對(duì)土壤收縮曲線影響較小，且與CK處理基本重合，而CO32-和HCO3-導(dǎo)致土壤收縮曲線發(fā)生較大變化，即土壤含水率相同時(shí)與CK處理相比，CO32-和HCO3-導(dǎo)致土壤比容積出現(xiàn)不同程度降低，使其3個(gè)收縮段與其他處理不同。綜上，SO42-、Cl-以及4種陽性鹽離子可維持土壤收縮特征曲線基本穩(wěn)定，其中Na+、SO42-和Cl-對(duì)土壤收縮特征曲線的作用效果甚小，而CO32-和HCO3-對(duì)其的影響效果極為明顯。

圖2不同鹽離子處理土壤收縮特征曲線

Fig.2 Soil shrinkage characteristic curves for different salt ion treatments

2.3 鹽基離子對(duì)土壤收縮度和裂隙的影響

離心過程導(dǎo)致土體不斷失水，并出現(xiàn)軸向沉降現(xiàn)象、容重逐漸增加，不同離子導(dǎo)致土壤容重增加幅度介于0.39～0.47 g·cm-3之間。圖3顯示了測(cè)定過程中和測(cè)定結(jié)束時(shí)的土體線縮率變化和軸向收縮應(yīng)變(δs)，與CK處理(δs=22.50%)相比，K+、Na+、Ca2+和Mg2+4種陽鹽離子處理δs值分別為23.21%、24.18%、23.20%和21.93%；Cl-、HCO3-、CO32-和SO42-4種陰鹽離子處理δs值分別為24.18%、24.40%、20.57%和21.63%。由此可知，在本研究設(shè)定的鹽溶液濃度(≈1 g·L-1)條件下，Mg2+、CO32-和SO42-有利于緩解土壤發(fā)生沉降、減輕土壤垂向收縮程度；另一方面，在吸力為300 cm和3 000 cm時(shí)全部處理土壤線縮率發(fā)生較大幅度增長(zhǎng)，故可初步將失水過程中的土體軸向沉降分為3個(gè)階段，即h=10～100 cm、h=300～1 000 cm和h=3 000～7 000 cm(圖3)。

在離心作用下，土體不斷失水，逐漸產(chǎn)生裂隙(本文省略裂隙飽和度分量示意圖)。不同鹽離子處理土壤裂縫的長(zhǎng)度、寬度、密度、走向及分布特征均存在差異，采用裂隙總長(zhǎng)度、總面積、面積密度和長(zhǎng)度密度4個(gè)指標(biāo)對(duì)土壤裂隙特征進(jìn)行評(píng)價(jià)(見表2)，與CK處理相比：(1) 對(duì)于4種陽性鹽離子，Na+和Mg2+可同時(shí)減小土壤裂隙總長(zhǎng)度和總面積以及相應(yīng)的長(zhǎng)度密度和面積密度，其中Na+效果較好，而K+和Ca2+在一定程度上可加劇裂縫產(chǎn)生，導(dǎo)致裂隙總長(zhǎng)度和總面積增加；(2) 對(duì)于4種陰性鹽離子，均有利于減輕土壤開裂水平，使得土壤裂隙總長(zhǎng)度和總面積及其長(zhǎng)度密度和面積密度均有所減小，其中CO32-和HCO3-效果較好，以CO32-最優(yōu)。

圖3 鹽離子對(duì)土壤收縮和容重的影響

3 討 論

邢旭光等[21]研究表明，不同離子對(duì)土壤持水性能及其在失水過程中的軸向收縮和裂隙水平具有不同的作用效果，且與離子濃度存在相關(guān)性；在其研究中，離子濃度是基于咸水(TDS=3～10 g·L-1)、鹽水(TDS=10～50 g·L-1)和鹵水(TDS>50 g·L-1)標(biāo)準(zhǔn)而設(shè)定。另一方面，我國(guó)鹽堿地面積廣闊，且土壤中可溶性鹽基離子類型多，而當(dāng)前關(guān)于鹽堿土壤水分特征曲線和收縮特征曲線的研究多是針對(duì)富含Na+的土壤而展開[19,22-24]，存在一定局限性。為消除試驗(yàn)隨機(jī)誤差、提高研究可信度，同時(shí)增加研究結(jié)果適用范圍，本研究在此基礎(chǔ)上進(jìn)一步擴(kuò)展了對(duì)陰性鹽離子的研究，且基于微咸水(TDS=1～3 g·L-1)標(biāo)準(zhǔn)將離子濃度設(shè)置為1 g·L-1，即本文對(duì)比研究了濃度為1 g·L-1的陰陽性8種鹽離子對(duì)土壤持水能力的影響，同時(shí)對(duì)比評(píng)價(jià)了土壤失水過程中的收縮水平，使得研究結(jié)果較為豐富、完善。鹽離子對(duì)土壤持水能力的影響主要受土壤結(jié)構(gòu)、團(tuán)聚體、離子價(jià)態(tài)等因素影響[21]，不同鹽離子具有不同的化學(xué)性質(zhì)、與土壤的親和力存在差異、土壤孔隙分布特征不同[25-26]，進(jìn)而對(duì)土壤物理性質(zhì)的作用效果以及影響機(jī)理有所差異。有必要將土壤物理與土壤化學(xué)相結(jié)合進(jìn)一步探索其作用機(jī)理。

本研究設(shè)定鹽基離子濃度均約為1 g·L-1，將本研究與文獻(xiàn)[21]的研究結(jié)果進(jìn)行結(jié)合：(1) 除5 g·L-1Na+外，其余處理均在一定程度上降低土壤持水能力[21]；本研究則表明在濃度很低時(shí)，K+和Na+均具有提高土壤持水能力的作用，進(jìn)一步證實(shí)離子濃度與土壤持水能力之間的負(fù)相關(guān)關(guān)系。(2) 隨著Mg2+濃度降低，土壤軸向收縮程度呈現(xiàn)減小趨勢(shì)，具有利于減小土壤沉降潛力[21]；本研究進(jìn)一步證實(shí)，與CK相比，1 g·L-1Mg2+可減小土壤軸向沉降，印證了上述推斷。(3) 本研究指出，1 g·L-1Na+和Mg2+亦有利于降低土壤裂縫評(píng)價(jià)指標(biāo)值，這對(duì)已有研究成果具有補(bǔ)充作用。綜上，土壤持水能力和土壤裂隙評(píng)價(jià)指標(biāo)均與鹽離子濃度具有顯著相關(guān)關(guān)系?？梢姳狙芯繉?duì)已有研究成果具有補(bǔ)充作用，完善了微咸水、咸水、鹽水和鹵水條件下的不同鹽離子對(duì)土壤持水和收縮特征的作用效果；同時(shí)對(duì)陰性鹽離子作用的研究具有一定探索性。

4 結(jié) 論

該研究中離子濃度是基于微咸水標(biāo)準(zhǔn)而設(shè)定，對(duì)于咸水、鹽水和鹵水條件是一種補(bǔ)充與完善，在已獲得研究成果基礎(chǔ)上進(jìn)行擴(kuò)展，進(jìn)一步得到以下結(jié)論：

1) 不同鹽離子導(dǎo)致土壤持水能力存在差異，濃度為1 g·L-1的K+、Na+和4種陰性鹽離子均可提高土壤持水能力，陰陽鹽離子分別以CO32-和K+效果最為顯著，HCO3-亦可以有效改善土壤持水性能。

2) 土體在失水過程中發(fā)生收縮現(xiàn)象，濃度為1 g·L-1的Mg2+、SO42-和CO32-有利于減小土體軸向收縮度，不同離子收縮應(yīng)變表現(xiàn)為Na+>K+>Ca2+>CK>Mg2+和HCO3->Cl->CK>SO42->CO32-，其中CO32-和HCO3-可明顯改變土壤收縮特征曲線形態(tài)。

3) 濃度為1 g·L-1的Na+、Mg2+和4種陰性鹽離子均可同時(shí)減小土壤裂隙總長(zhǎng)度和總面積以及相應(yīng)的長(zhǎng)度密度和面積密度；陰陽鹽離子中分別以CO32-和Na+效果最為顯著。

參考文獻(xiàn)：

[1] 冉艷玲，王益權(quán)，張潤(rùn)霞，等.保水劑對(duì)土壤持水特性的作用機(jī)理研究[J].干旱地區(qū)農(nóng)業(yè)研究，2015，33(5)：101-107.

[2] 邢旭光，趙文剛，馬孝義，等.土壤水分特征曲線測(cè)定過程中土壤收縮特性研究[J].水利學(xué)報(bào)，2015，46(10)：1181-1188.

[3] Zhou W H, Yuen K V, Tan F. Estimation of soil-water characteristic curve and relative permeability for granular soils with different initial dry densities[J]. Engineering Geology, 2014,179:1-9.

[4] Mohammadi M H, Meskini-Vishkaee F. Predicting soil moisture characteristic curves from continuous particle-size distribution data[J]. Pedosphere, 2013,23(1):70-80.

[5] Shwetha P, Varija K. Soil water retention curve from saturated hydraulic conductivity for sandy loam and loamy sand textured soils[J]. Aquatic Procedia, 2015,4:1142-1149.

[6] Thyagaraj T, Rao S M. Influence of osmotic suction on the soil-water characteristic curves of compacted expansive clay[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2010,136(12):1695-1702.

[7] Gao H, Shao M. Effects of temperature changes on soil hydraulic properties[J]. Soil and Tillage Research, 2015,153:145-154.

[8] Ugarte-Nano C C, Nicolardot B, Quinche M, et al. Effects of integrated weed management based cropping system on the water retention of a silty clay loam soil[J]. Soil and Tillage Research, 2016,156:74-82.

[9] 王志超，李仙岳，史海濱，等.農(nóng)膜殘留對(duì)土壤水動(dòng)力參數(shù)及土壤結(jié)構(gòu)的影響[J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2015，46(5)：101-106，140.

[10] Leung A K, Garg A, Ng C W W. Effects of plant roots on soi-water retention and induced suction in vegetated soil[J]. Engineering Geology, 2015,193:183-197.

[11] 邢旭光，柳 燁，馬孝義.土壤添加物對(duì)土-水曲線和土體收縮的影響[J].水科學(xué)進(jìn)展，2016，27(1)：40-48.

[12] 高會(huì)議，郭勝利，劉文兆，等.不同施肥土壤水分特征曲線空間變異[J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2014，45(6)：161-165，176.

[13] Stoof C R, Wesseling J G, Ritsema C J. Effects of fire and ash on soil water retention[J]. Geoderma. 2010,159:276-285.

[14] Ren Jianhua, Li Xiaojie, Zhao Kai, et al. Study of an on-line measurement method for the salt parameters of soda-saline soils based on the texture features of cracks[J]. Geoderma, 2016,263:60-69.

[15] He B, Cai Y, Ran W, et al. Spatial and seasonal variations of soil salinity following vegetation restoration in coastal saline land in eastern China[J]. Catena, 2014,118:147-153.

[16] Pauchard L, Parisse F, Allain C. Influence of salt content on crack patterns formed through colloidal suspension desiccation[J]. Physical Review E, 1999,59:3737.

[17] Lima L A, Grismer M E. Soil crack morphology and soil salinity[J]. Soil Science, 1992,153(2):149-153.

[18] 張中彬，彭新華.土壤裂隙及其優(yōu)先流研究進(jìn)展[J].土壤學(xué)報(bào)，2015，52(3)：477-488.

[19] Fattah M Y, Yahya A Y, Al-Hadidi M Th, et al. Effect of salt content on total and matric suction of unsaturated soils[J]. European Scientific Journal, 2013,9:228-245.

[20] 楊紹鍔，黃元仿.關(guān)于土壤收縮特征曲線的探討[J].土壤通報(bào)，2007，38(4)：749-752.

[21] 邢旭光，馬孝義，康端剛.鹽陽離子類型及濃度對(duì)土壤持水及干縮開裂的作用效果[J].農(nóng)業(yè)工程學(xué)報(bào)，2016，32(9)：115-122.

[22] 栗現(xiàn)文，周金龍，靳孟貴，等.高礦化度土壤水分特征曲線及擬合模型適宜性[J].農(nóng)業(yè)工程學(xué)報(bào)，2012，28(13)：135-141.

[23] 孫凱強(qiáng)，唐朝生，王 鵬，等.鹽分對(duì)土體干縮開裂的影響[J].工程地質(zhì)學(xué)報(bào)，2015，23(增刊)：77-83.

[24] 栗現(xiàn)文，周金龍，趙玉潔，等.高礦化度對(duì)砂性土毛細(xì)水上升的影響[J].農(nóng)業(yè)工程學(xué)報(bào)，2011，27(8)：84-89.

[25] 李學(xué)垣.土壤化學(xué)[M].北京：高等教育出版社，2001.

[26] 李海云，王秀峰.不同陰離子化肥對(duì)設(shè)施土壤理化性狀的影響研究[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，2004，12(4)：126-128.