摘要：探究石河子墾區(qū)土壤鹽分空間分布特征及其環(huán)境影響因素，可為該地區(qū)鹽堿地改良和土地可持續(xù)利用提供科學(xué)依據(jù)。以石河子墾區(qū)為基礎(chǔ)，通過均勻布點(diǎn)及樣本采集，利用方差分析和灰色關(guān)聯(lián)度分析等方法評估石河子墾區(qū)土壤鹽分的空間分布，并分析海拔、坡度、巖性、土地利用方式和土壤類型等環(huán)境因子對其影響。結(jié)果表明：（1）石河子墾區(qū)0～20、20～40、40～60 cm土層含鹽量均值分別為2.14、3.89、5.00 g·kg^-1，其中0～20 cm土層主要為非鹽化土，20～60 cm土層呈現(xiàn)輕度鹽漬化，各土層鹽分含量在空間分布上呈現(xiàn)中到強(qiáng)的變異性，并與結(jié)構(gòu)性和隨機(jī)性的因素都有關(guān)；三層土壤都為中等強(qiáng)度變異，在0～20 cm土層土壤具有中等程度相關(guān)性，由隨機(jī)性因素和結(jié)構(gòu)性因素共同作用，20～60 cm土層土壤具有強(qiáng)空間相關(guān)性，由結(jié)構(gòu)性因素主導(dǎo)。（2）0～20 cm土層鹽分高值區(qū)主要集中在墾區(qū)北部，20～40、40～60 cm土層鹽分高值區(qū)主要集中在墾區(qū)北部和西部；從空間分布上看，鹽漬化土壤主要在20～40、40～60 cm土層，且以輕度鹽化土為主，分別占總面積的69.14%、58.34%，中度和重度鹽化土則以條帶狀和斑狀鑲嵌于輕度鹽化土內(nèi)。（3）不同環(huán)境因子與各土層鹽分的關(guān)聯(lián)度關(guān)系為土壤類型gt;土地利用方式gt;海拔gt;巖性gt;坡度，研究區(qū)西北方向的土壤鹽漬化程度較東南方向的嚴(yán)重；土壤類型、土地利用方式和海拔是影響土壤鹽分的主控因子，裸地、海拔較低區(qū)域以及土壤類型為草甸土、潮土、鹽土和沼澤土的地區(qū)是未來土壤鹽漬化預(yù)防重點(diǎn)區(qū)域。

關(guān)鍵詞：石河子墾區(qū)；土壤鹽分；空間變異性；影響因素；灰色關(guān)聯(lián)度分析

中圖分類號(hào)：S151.9文獻(xiàn)標(biāo)志碼：文獻(xiàn)標(biāo)識(shí)碼A

DOI：10.13880/j.cnki.65-1174/n.2024.21.018

文章編號(hào)：1007-7383（2025）01-0092-10

Spatial variation of soil salinity and analysis of influencing factors in the reclamation area of Shihezi，Xinjiang

YUAN" Yushan^{1，2，3，4}，L Desheng^{1，2，3，4 *}，WANG" Zhenhua^{1，2，3，4} ，LI" Haiqiang^{1，2，3，4}，MA" Zhanli^{1，2，3，4}

（1 College of Water Conservancy amp; Architectural Engineering，Shihezi University，Shihezi，Xinjiang 832000，China;

2 Key Laboratory of Modern Water-saving Irrigation of Xinjiang Production amp; Construction Corps，Shihezi，Xinjiang 832000，China;

3 Technology Innovation Center for Agricultural Water and Fertilizer Efficiency Equipment of Xinjiang Production amp; Construction Corps，Shihezi，Xinjiang 832000，China;

4 Key Laboratory of Northwest Oasis Water-saving Agriculture， Ministry of Agriculture and Rural Affairs，Shihezi，Xinjiang 832000，China）

Abstract： To explore the characteristics of spatial distribution of soil salinity in Shihezi Reclamation Area and its influencing factors can provide scientific basis for saline and alkaline land improvement and land utilization in this area.Based on Shihezi Reclamation Area，the spatial distribution of soil salinity in Shihezi Reclamation Area was evaluated through uniform distribution and sample collection，using variance and gray correlation，etc，and the influence of factors such as elevation，slope，lithology，land use and soil type was also evaluated.The results show that：（1）the mean values of salinity in the 0～20，20～40，40～60 cm soil layers in Shihezi Reclamation Area were 2.14，3.89 and 5.00 g-kg^-1，respectively，of which the 0～20 cm soil layer was mainly non-salinized，while the 20～60 cm soil layer was mildly salinized，and the spatial distribution of salinity content of each soil layer showed moderate to strong variability，which was related to both structural and stochastic factors.All three soil layers showed medium-strength variability，with moderate correlation in the 0～20 cm soil layer，which was jointly contributed by stochastic and structural factors，and strong spatial correlation in the 20～60 cm soil layer，which was dominated by structural factors.（2）The salinity high value area of 0～20 cm soil layer is mainly concentrated in the northern part of the reclamation area，and the salinity high value areas of 20～40，40～60 cm soil layers are mainly concentrated in the northern and western parts of the reclamation area.In terms of spatial distribution，salinized soils are mainly found in the 20～40，40～60 cm soil layers，and are dominated by lightly salinized soils，accounting for 69.14%，58.34% of the total area respectively，while moderately and heavily salinized soils are embedded in lightly salinized soils in the form of stripes and patches.（3）The correlation relationship between different factors and salinity of each soil layer was soil type gt; land use gt; elevation gt; lithology gt; slope.Soil salinization was more serious in the northwest than in the southeast of the study area.Soil type，land use and elevation were the main controlling factors affecting soil salinity.Bare land，lower elevation areas，and areas with soil types of meadow soil，tidal soil，salt soil，and marsh soil are the key areas for soil salinization prevention in the future.

Key words： Key words;Shihezi reclamation area;Soil salinity;Spatial variability;Influencing factors;Grey correlation degree

新疆是我國最大的優(yōu)質(zhì)棉生產(chǎn)基地，中國約有31.1%的耕地受到鹽漬化威脅^［1］，其中，新疆是我國受鹽堿威脅最嚴(yán)重的省份之一；膜下滴灌的長期不當(dāng)實(shí)施將導(dǎo)致土壤次生鹽漬化^［2］，嚴(yán)重威脅農(nóng)業(yè)可持續(xù)發(fā)展。在采用灌溉的地區(qū)減輕和控制農(nóng)業(yè)用地的土壤鹽分至關(guān)重要^［3］，了解土壤鹽分積累的過程和空間變化規(guī)律，將有助于預(yù)防土壤退化，改善土壤質(zhì)量，促進(jìn)經(jīng)濟(jì)發(fā)展^［4］，這就要求在大尺度上準(zhǔn)確估計(jì)土壤鹽漬化的空間變化^［5］。

因此，研究石河子墾區(qū)土壤鹽分的空間分布特征對當(dāng)?shù)赝寥利}漬化防治、改良及農(nóng)業(yè)生產(chǎn)具有重要意義^［6］。

近年來，國內(nèi)外眾多學(xué)者對大尺度鹽漬化問題開展了大量研究，Wang等^［7］對1982—2015年中國西北干旱區(qū)三公河流域典型內(nèi)陸小流域0～20 cm表層土壤鹽分遷移的研究發(fā)現(xiàn)土壤鹽分有從上游向下游遷移的趨勢；Benyamini等^［8］通過在以色列耶斯列河谷灌溉區(qū)的研究發(fā)現(xiàn)，安裝排水系統(tǒng)2～3年后，土壤鹽度明顯降低，且在春季將地下水位降至1 m以下能很好的防治土壤鹽漬化問題；Huang等^［9］利用地理信息系統(tǒng)（GIS）對東營市土壤鹽分和地貌空間變化之間關(guān)系的研究結(jié)果表明，當(dāng)?shù)赝寥利}分含量從孤立洼地的區(qū)域性高值下降到有堤防的前背沼澤的相對中等值，在廢棄河道中出現(xiàn)最低值。

石河子墾區(qū)作為國內(nèi)最早使用膜下滴灌技術(shù)之一的地區(qū)，當(dāng)?shù)氐耐寥利}漬化問題受到眾多學(xué)者的關(guān)注。

李玉義^［10］等通過對石河子墾區(qū)不同地貌條件下鹽分分布特征的研究表明，土壤含鹽量為沖積扇邊緣gt;沖積扇平原gt;干三角洲；宗含^［11］等通過對石河子墾區(qū)土壤鹽分與地下水埋深關(guān)系的研究表明，地下水埋深與土壤鹽分含量成反比，且對土壤含鹽量影響較為明顯；鄭埼^［12］等研究表明，隨著滴灌年限的增加，土壤含鹽量呈現(xiàn)降低的趨勢；谷海斌^［13］等通過對土壤鹽分和河流流向關(guān)系的研究表明，河流下游的土壤含鹽量相對較高；顏安^［14］等對石河子墾區(qū)不同種植作物土壤鹽分含量的研究結(jié)果表明，不同種植作物下土壤鹽分狀況表現(xiàn)為棉花地gt;葡萄地gt;玉米地。

綜合分析以上研究結(jié)果后可知，大尺度土壤鹽漬化的空間分布特征受到地形地勢、地貌、種植制度、土地利用和土壤質(zhì)地及地下水等多種因素的影響。近年來針對石河子墾區(qū)土壤鹽分大尺度空間分布的研究較為局限，同時(shí)以往由膜下滴灌應(yīng)用初期所得到的結(jié)果可能與現(xiàn)在的實(shí)際情況有所差異，并且前人對石河子墾區(qū)土壤鹽分分布的研究大多針對地貌、地下水和種植制度展開，當(dāng)?shù)赝寥篮}量與其它因素之間的關(guān)系尚不明確，因此，本文綜合海拔、坡度、巖性、土壤類型和土地利用類型這5種環(huán)境因子對石河子墾區(qū)土壤鹽分剖面特征和空間分布特征的影響進(jìn)行對比與分析，通過田間取樣、化學(xué)分析和地統(tǒng)計(jì)學(xué)方法研究了石河子墾區(qū)土壤鹽分的空間分布特征，并利用方差分析與灰色關(guān)聯(lián)度分析研究區(qū)土壤含鹽量與環(huán)境因子的關(guān)系，旨在深入了解研究區(qū)土壤鹽漬化的影響因素，得出其潛在鹽漬化區(qū)域，從而為石河子墾區(qū)的鹽堿地治理與土地可持續(xù)利用提供一定的科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

石河子墾區(qū)（圖1）地處天山北麓中段（84°43′～86°35′E，43°21′～45°20′N），古爾班通古特沙漠南緣，全墾區(qū)面積為7.53×103 km2；該墾區(qū)下轄17個(gè)農(nóng)場、1個(gè)牧場，種植農(nóng)作物以棉花為主，采用膜下滴灌水肥一體化技術(shù)；該墾區(qū)5月（棉花出苗期）滴灌次數(shù)為1次，8月（鈴期）滴灌次數(shù)4～6次，單次灌水量270～330 m3·km^-2，10月（收獲期）停止灌溉，年施肥量1020～1125 kg·hm^-2［10］。該墾區(qū)地勢整體上由南向北微微傾斜，坡度在0°～15°之間，平均海拔450.8 m；干旱少雨，屬于典型的溫帶大陸性氣候，年均氣溫為4.7～5.7 ℃，年均蒸發(fā)量在1 500～2100 mm，年均降水量主要集中于100～200 mm，各地蒸發(fā)量與降雨量差異不大；土地利用類型多樣，土壤類型以草甸土、灰漠土、風(fēng)沙土、潮土、沼澤土和鹽土為主，巖性以砂泥、砂、砂石和礫石為主。

1.2 樣品采集與樣品處理

本研究于2023年10月24—31日使用網(wǎng)格法在研究區(qū)域內(nèi)按間距為8 km均勻布設(shè)取樣點(diǎn)，并根據(jù)實(shí)際情況進(jìn)行適當(dāng)調(diào)整，兼顧取樣點(diǎn)的代表性和均勻性，每個(gè)取樣地塊用“S”形取樣法確定4個(gè)取樣點(diǎn)，每個(gè)取樣點(diǎn)選取0～20、20～40、40～60 cm土層利用取土鉆進(jìn)行分層取樣，對4個(gè)取樣點(diǎn)土樣分層混勻，共采集432份土壤樣品，將其放入自封袋中密封保存（取樣點(diǎn)分布見圖2），同時(shí)使用GPS工具箱記錄樣點(diǎn)的經(jīng)緯度坐標(biāo)、海拔和坡度，同時(shí)調(diào)查記錄樣點(diǎn)巖性以及土壤類型和土地利用方式等信息。對所采集的土壤樣品風(fēng)干、研磨后過2 mm篩，按1∶5土水質(zhì)量比浸提，測定其電導(dǎo)率EC1∶5（DDS-11A型，上海大普儀器有限公司），選擇具有代表性的浸提液采用烘干法測定土壤中總鹽含量，并對其電導(dǎo)率值進(jìn)行擬合，得到含鹽量與電導(dǎo)率的回歸關(guān)系，擬合的結(jié)果如下：

SCC=0.004 6EC1∶5＋0.526 4R2=0.954，（1）

其中，SCC為含鹽量，g·kg^-1；EC1∶5為土壤浸提液的電導(dǎo)率，μS·cm^-1。

1.3 研究方法

本文采用多種統(tǒng)計(jì)分析方法對土壤鹽分?jǐn)?shù)據(jù)進(jìn)行探討。先利用Kolmogorov-Smirnov檢驗(yàn)進(jìn)行正態(tài)性檢驗(yàn)，在5%的置信度下判斷數(shù)據(jù)是否服從正態(tài)分布^［15］，再使用熱點(diǎn)分析確定數(shù)據(jù)中的顯著熱點(diǎn)區(qū)域^［16］；然后通過半變異函數(shù)擬合土壤鹽分樣本數(shù)據(jù)，并選取塊金值（Nugget，C0）、基臺(tái)值（Sill，C0+C）和變程（Range，A0）作為主要參數(shù)，在決定系數(shù)R2最大、殘差RSS最小的條件下得到最佳模型。

本研究采用普通Kriging方法^［17］對土壤鹽分進(jìn)行插值。為獲得最佳插值結(jié)果，在插值過程中采用多種方法且對結(jié)果進(jìn)行比較，以找出最適合特定應(yīng)用的數(shù)據(jù)插值方法；再將半變異函數(shù)的擬合參數(shù)導(dǎo)入ArcGIS進(jìn)行交叉驗(yàn)證，以檢驗(yàn)插值結(jié)果的精度；最后運(yùn)用方差分析評估環(huán)境因子對土壤含鹽量的影響，并采用灰色關(guān)聯(lián)度分析各因子與含鹽量之間的關(guān)聯(lián)程度。

1.4 數(shù)據(jù)分析

使用統(tǒng)計(jì)軟件SPSS24.0對土壤樣本數(shù)據(jù)統(tǒng)計(jì)分析和正態(tài)檢驗(yàn)，以及對異常值的剔除和方差分析，隨后通過GS+9.0地統(tǒng)計(jì)工具評估土壤鹽分的空間變異并確定模型，最后在ArcGIS10.8.1平臺(tái)上進(jìn)行熱點(diǎn)分析和應(yīng)用Kriging方法插值，以構(gòu)建土壤含鹽量高低值區(qū)和土壤鹽分空間分布圖，常規(guī)制圖在OriginPro 2021中完成。

2 結(jié)果與分析

2.1 土壤鹽分統(tǒng)計(jì)特征分析及正態(tài)分布檢驗(yàn)

通過統(tǒng)計(jì)軟件SPSS24.0對土壤樣本數(shù)據(jù)進(jìn)行基本統(tǒng)計(jì)分析，剔除0～20 cm土層含鹽量的離群值6.83、7.48、7.61、7.89、9.37 g·kg^-1和20～40 cm土層含鹽量的離群值10.50、11.17 g·kg^-1以及40～60 cm土層含鹽量的離群值15.55、16.32 g·kg^-1；同時(shí)根據(jù)《新疆土壤》分級標(biāo)準(zhǔn)對鹽漬土等級進(jìn)行劃分，含鹽量lt;3 g·kg^-1為非鹽化土，3～6 g·kg^-1為輕度鹽化土，6～10 g·kg^-1為中度鹽化土，10～20 g·kg^-1為重度鹽化土，gt;20 g·kg^-1為鹽土。按上述方法及過程得到研究區(qū)不同土層含鹽量統(tǒng)計(jì)分析結(jié)果見表1，可見：

（1）土壤含鹽量在空間上存在顯著的差異性，0～20、20～40 、40～60 cm土層土壤含鹽量的范圍分別在0.17～5.97、0.71～9.48、0.12～13.48 g·kg^-1；在0～60 cm土層，土壤含鹽量的均值和最大值隨土層深度的增加而增大，而土壤含鹽量最小值呈現(xiàn)先增大后減小的趨勢；研究區(qū)0～20 cm土層以非鹽化土壤（含鹽量lt;3 g·kg^-1）為主，而20～40、40～60 cm土層不同程度的鹽化土壤（含鹽量gt;3 g·kg^-1）占比大于50%。

（2）研究區(qū)各土層含鹽量都屬于中等強(qiáng)度變異（變異系數(shù)Cv反映土壤參數(shù)的變異程度，Cv≤10%時(shí)呈弱變異性，10%＜Cv≤100%時(shí)呈中等變異性，Cv＞100%時(shí)呈強(qiáng)變異性），0～20、20～40、40～60 cm的土層含鹽量變異系數(shù)分別為86.07%、60.57%、68.05%。

各土層土壤含鹽量經(jīng)對數(shù)轉(zhuǎn)換后服從正態(tài)分布，滿足地統(tǒng)計(jì)學(xué)分析基本要求，所以可采用ArcGIS進(jìn)行普通Kriging插值，并對研究區(qū)鹽分的分布進(jìn)行預(yù)測，然后采用皮爾遜相關(guān)性分析，結(jié)果（表2）表明：研究區(qū)0～60 cm不同土層間含鹽量均存在極顯著相關(guān)性（Plt;0.01），20～40 cm土層含鹽量與40～60 cm土層相關(guān)性最顯著（P=0.756）。

2.2 熱點(diǎn)分析

對研究區(qū)土壤含鹽量進(jìn)行熱點(diǎn)分析，結(jié)果（圖3）顯示：0～20 cm土層含鹽量高值（熱點(diǎn)）主要焦聚區(qū)分布在研究區(qū)的北部，低值（冷點(diǎn)）焦聚區(qū)分布在研究區(qū)的西南和北部，而西部和中部則表現(xiàn)為不顯著。20～40、40～60 cm土層?xùn)|南部為低值（冷點(diǎn)）焦聚區(qū)，且存在小面積的不顯著區(qū)域，西部和北部地區(qū)均為高值（熱點(diǎn)）焦聚區(qū)。

2.3 土壤含鹽量的空間變異函數(shù)

利用GS+9.0進(jìn)行半方差函數(shù)計(jì)算及模型擬合，得到最優(yōu)模型和塊金值C0、基臺(tái)值C0+C、塊金系數(shù)C0/（C0+C）、變程A0、殘差RSS、決定系數(shù)R2等模型參數(shù)（表3）。其中，塊金值C0反映研究區(qū)內(nèi)變量隨機(jī)性變異程度；基臺(tái)值C0+C反映研究區(qū)內(nèi)變量的變異程度；塊金系數(shù)C0/（C0+C）表示研究區(qū)內(nèi)變量由隨機(jī)性因素部分引起的空間變異程度，若＜25%時(shí)表明具有強(qiáng)烈的空間相關(guān)性，且由結(jié)構(gòu)因素（地形地貌、土壤類型、土地利用等）引起的變異程度大，若＞75%時(shí)表明空間相關(guān)性較弱，由隨機(jī)性因素（灌水、施肥、耕地等）引起的變異程度大，在25%和75%之間時(shí)，表明有中等空間相關(guān)性，由結(jié)構(gòu)性因素和隨機(jī)性因素共同主導(dǎo)；變程A0表示變量具有空間相關(guān)性的閾值大小，超出范圍則相關(guān)性減弱^［16］；模型評價(jià)以殘差RSS最小，決定系數(shù)R2最大為最優(yōu)。

由表3可知：根據(jù)RSS和R2得到0～20、40～60cm土層含鹽量最優(yōu)半方差理論模型是球形，20～40 cm土層最優(yōu)半方差理論模型符合高斯模型。0～20 cm土壤的塊金系數(shù)為49.95%，具有中等空間相關(guān)性，而20～40、40～60 cm土層塊金系數(shù)分別為0.23%、0.14%，具有強(qiáng)烈的空間相關(guān)性，由結(jié)構(gòu)性因素引起的變異占主體。各層次土壤含鹽量半方差模型的變程均大于本試驗(yàn)所采用的取樣間距8 km，表明該取樣尺度已滿足研究的要求。

2.4 土壤含鹽量的空間分布特征

為進(jìn)一步準(zhǔn)確且直觀地表現(xiàn)研究區(qū)土壤鹽分空間分布特征，利用擬合的變異函數(shù)模型進(jìn)行普通Kriging插值，得到土壤含鹽量分布圖（圖4），并對插值結(jié)果進(jìn)行精度驗(yàn)證，結(jié)果（表4）顯示：最優(yōu)擬合""模型的標(biāo)準(zhǔn)為平均誤差接近0，標(biāo)準(zhǔn)均方誤差趨近于1，均方根誤差達(dá)到最小，平均標(biāo)準(zhǔn)誤差和均方根誤差大致相等。上述結(jié)果表明模型擬合情況較好。

計(jì)算鹽漬土所占面積及比例，結(jié)果見表5。

由表5和圖4可知：研究區(qū)0～20 cm土層以非鹽化土（lt;3 g·kg^-1）為主，占研究區(qū)總面積的93.03%，僅有占比6.97%的輕度鹽化土（3～6 g·kg^-1）以塊狀分布在研究區(qū)東南部和西北部。在20～40 cm土層，非鹽化土占比為16.08%，重度鹽化土（10～20 g·kg^-1）占比極少，僅為0.05%，以輕度鹽化土居多，占總面積的69.14%，且分布無明顯規(guī)律，中度鹽化土（6～10 g·kg^-1）占比為14.74%，以條塊狀和點(diǎn)狀鑲嵌于輕度鹽化土中。在40～60 cm，相較于20～40 cm土層，非鹽化土和輕度鹽化土面積占比減少，分別為7.32%、58.34%，中度鹽化土和重度鹽化土占比增大，分別為32.67%、1.37%，各類鹽化土的分布格局與20～40 cm土層較為相似。

2.5 土壤含鹽量影響因素的分析

2.5.1 海拔對不同土層含鹽量的影響

各層土壤含鹽量隨海拔的升高均表現(xiàn)為先增大后減小的趨勢，其中0～20 cm土層含鹽量在各個(gè)海拔梯度無顯著性差異，而20～40、40～60 cm土層含鹽量均表現(xiàn)為海拔lt;600 m的區(qū)域含鹽量顯著高于gt;600 m的區(qū)域（Plt;0.05）。

2.5.2 坡度對不同土層含鹽量的影響

方差分析結(jié)果（圖5）顯示：0～20 cm土層含鹽量隨坡度的增加表現(xiàn)為先增大后減小的趨勢，而20～40、40～60 cm土層含鹽量隨坡度的增加呈現(xiàn)出逐漸減小的趨勢，其中0～20、40～60 cm土層含鹽量均表現(xiàn)為坡度lt;9°的區(qū)域顯著高于坡度gt;9°的區(qū)域（Plt;0.05），而20～40 cm土層無顯著性差異。

2.5.3 土地利用方式對不同土層含鹽量的影響

方差分析結(jié)果（圖5）顯示：不同土地利用方式中，0～20 cm土層含鹽量由大到小依次為裸地gt;林地gt;草地gt;耕地gt;城鄉(xiāng)，其中裸地含鹽量顯著大于其它土地利用方式（Plt;0.05），20～40、40～60 cm土層含鹽量由大到小依次為裸地gt;林地gt;耕地gt;草地gt;城鄉(xiāng)，裸地含鹽量顯著高于其它土地利用方式（Plt;0.05），耕地和林地含鹽量顯著高于草地和城鄉(xiāng)（Plt;0.05），草地含鹽量顯著高于城鄉(xiāng)（Plt;0.05）。

2.5.4 土壤類型對不同土層含鹽量的影響

方差分析結(jié)果（圖5）顯示：不同土壤類型中，0～20 cm土壤含鹽量由大到小依次為：草甸土gt;沼澤土gt;潮土gt;鹽土gt;風(fēng)沙土gt;灰漠土，灰漠土含鹽量顯著低于其他土壤類型（Plt;0.05），20～40、40～60 cm土層含鹽量由大到小依次為草甸土gt;沼澤土gt;鹽土gt;潮土gt;風(fēng)沙土gt;灰漠土，灰漠土含鹽量顯著低于風(fēng)沙土（Plt;0.05），且兩者含鹽量均顯著低于其他類型土壤含鹽量（Plt;0.05）。

2.5.5 巖性對不同土層含鹽量的影響

方差分析結(jié)果（圖5）顯示：在不同粒徑巖性顆粒組成的土壤中，由含粒徑lt;0.0625 mm（砂泥）巖石碎屑的土壤含鹽量顯著高于由粒徑gt;0.0625 mm（砂、砂石、礫石和卵石）巖石碎屑組成的土壤（Plt;0.05），而沖洪積卵礫石層含鹽量顯著高于后者（Plt;0.05），與前者無顯著差異。

2.6 灰色關(guān)聯(lián)度分析

結(jié)果（表6）顯示：在0～20、20～40、40～60 cm土層中，海拔、坡度、巖性、土地利用方式和土壤類型對不同土層含鹽量的關(guān)聯(lián)度表現(xiàn)一致。這表明這些因素在不同深度的土層中對土壤含鹽量的影響程度相似，并且它們在整個(gè)土壤剖面中具有較為穩(wěn)定的關(guān)聯(lián)特征，其中土壤類型與土壤含鹽量關(guān)聯(lián)度最高，相關(guān)系數(shù)為0.944，其次為土地利用方式、海拔和巖性，相關(guān)系數(shù)分別為0.916、0.884、0.678，坡度與土壤含鹽量關(guān)聯(lián)度最低，為0.461。

3 討論

3.1 土壤鹽分的空間分布特征

（1）從垂直方向看，研究區(qū)0～60 cm土層的非鹽化土（lt;3 g·kg^-1）和的分布狀態(tài)呈表聚形，中度鹽化土（6～10 g·kg^-1）和重度鹽化土（10～20 g·kg^-1）的分布狀態(tài)呈底聚形，而輕度鹽化土（3～6 g·kg^-1）的分布趨勢隨著深度增加呈現(xiàn)先增加后減小的趨勢（表1）。這與喬江飛等^［17］的研究結(jié)果略有不同。另外，吳英^［18］研究表明土壤剖面的含鹽量分布特征呈現(xiàn)季節(jié)性變化，表現(xiàn)為春季返鹽-夏季脫鹽-秋季積鹽-冬季穩(wěn)定的變化趨勢；陳霖明等^［19］對瑪納斯中游滴灌棉田土壤鹽分含量的研究發(fā)現(xiàn)，土壤鹽分峰值出現(xiàn)在6～8月，土壤表層含鹽量明顯下降，而中層和深層含鹽量上升?？梢姡狙芯拷Y(jié)果與前人的有差異，這可能是取樣時(shí)間不同造成的結(jié)果，本研究結(jié)果僅能反映秋季積鹽時(shí)期的鹽分特征，若要了解研究區(qū)鹽分隨季節(jié)變化的趨勢，則要進(jìn)一步分析春季、夏季和冬季的土壤鹽分特征。

變異系數(shù)直接反映土壤鹽分空間變異程度。本研究表明土壤鹽分變異系數(shù)表層大于底層（表1），這是因?yàn)楦乇韺又苯颖┞队诃h(huán)境之中，易受到自然和人為因素的影響，因此變異系數(shù)最大^［20］。0～20、20～40、40～60 cm土壤含鹽量的塊金系數(shù)分別為49.95%、0.23%、0.14%，其中表層土壤含鹽量塊金系數(shù)明顯大于宗含等^［11］研究的結(jié)果，這可能是因?yàn)楸狙芯咳佑谧魑锸斋@期，耕地表層受到農(nóng)機(jī)碾壓與翻地的影響，因此增大了隨機(jī)性因素引起的空間變異性。

（2）從水平方向看，石河子墾區(qū)土壤鹽分高值區(qū)整體上位于研究區(qū)北部和西部（圖3）。相關(guān)研究的結(jié)果表明：研究區(qū)西北部成土母質(zhì)含可溶性鹽分高^［21］；東北部與沙漠接壤，排水不暢且蒸發(fā)較為強(qiáng)烈，鹽分隨水上移到耕層使得含鹽量高^［22］；西南部由于地下水埋深淺，在持續(xù)的蒸發(fā)作用下深層土壤及地下水中的鹽分借助毛細(xì)管作用上升，造成土壤鹽分含量較高^［23］。本文研究所得的鹽分水平空間分布特征與前人研究結(jié)果基本一致。

3.2 環(huán)境因子對土壤鹽分分布的影響

（1）研究區(qū)土壤剖面塊金系數(shù)都lt;50%，表明研究區(qū)鹽分的空間變異性主要受結(jié)構(gòu)性因素的影響；研究區(qū)鹽分空間分布特征與海拔、土壤類型、土地利用方式和巖性都有關(guān)，而受坡度的影響較小。含鹽量隨海拔的增高表現(xiàn)為先增高（lt;400 m）再降低（gt;400 m）的趨勢，含鹽量最高值在海拔300～400 m，這與鄭琦等^［12］研究的結(jié)果基本一致，然而顏安等^［14］研究的結(jié)果是鹽分最高值出現(xiàn)在在400～500 m范圍內(nèi)，這一差異可能是因?yàn)楦鞣N研究區(qū)的樣本點(diǎn)分布與數(shù)量的不同。由于研究區(qū)多數(shù)位置海拔位于300～400 m，而本研究中采用網(wǎng)格取樣法取樣，300～400 m范圍內(nèi)的土壤樣本數(shù)量相對較多，這種取樣策略導(dǎo)致了對該海拔范圍內(nèi)土壤特征的覆蓋更為充分，涵蓋了更多不同的土壤類型和土地利用方式等特征，這些因素共同導(dǎo)致了本研究結(jié)果與前人研究，結(jié)果的差異。

（2）不同的時(shí)段也可能對土壤鹽分含量造成影響。在不同土地利用類型下，裸地土壤含鹽量顯著大于耕地、林地、草地和城鄉(xiāng)，說明對土地進(jìn)行開發(fā)利用可以減輕土壤鹽漬化，地表覆蓋度會(huì)對土壤含鹽量造成不同程度的影響^［24］。在不同土壤類型中，草甸土、沼澤土、鹽土和潮土含鹽量明顯大于風(fēng)沙土和灰漠土，這與王雪梅等^［25］研究結(jié)果相似；巖石碎屑顆粒大小對土壤含鹽量有一定的影響，較大的顆粒會(huì)使土壤的孔隙較大，導(dǎo)致鹽分在土壤中更容易擴(kuò)散和稀釋，從而土壤含鹽量較低，沖洪積卵礫石層雖以粒徑較大的卵石和礫石為主，但沖洪積層常處于地下水位較高的位置；另外，竇旭等^［26］、王國帥等^［27］研究表明，地下水埋深越淺，土壤積鹽越嚴(yán)重，因此沖洪積層含鹽量較高。

本研究表明5種環(huán)境因子對不同深度土層的土壤含鹽量的影響程度相似，關(guān)聯(lián)度大小表現(xiàn)為土壤類型（0.944）gt;土地利用方式（0.916）gt;海拔（0.884）gt;巖性（0.678）gt;坡度（0.461），說明5種環(huán)境因子中，土壤類型、土地利用方式和海拔是土壤含鹽量的主控因子，而巖性和坡度對土壤含鹽量的影響較小。

3.3 鹽漬土改良與展望

（1）本文以石河子墾區(qū)為研究區(qū)域，通過在該區(qū)域內(nèi)均勻布點(diǎn)，在作物收獲后采集土樣，利用普通Kriging插值法繪制土壤鹽分分布圖。這一方法能夠客觀地反映研究區(qū)土壤鹽漬化特征的整體情況，為在該區(qū)域內(nèi)實(shí)施針對性的土壤改良措施提供科學(xué)依據(jù)。

（2）本研究發(fā)現(xiàn)，土壤類型為草甸土、潮土、鹽土和沼澤土的土壤含鹽量顯著高于其他類型的土壤，研究區(qū)西南部草甸土和鹽土占比極高（圖1），而這兩種土壤類型的形成與當(dāng)?shù)氐叵滤裆钶^淺^［27］的水文條件有關(guān)。因此，通過建立排水系統(tǒng)，降低地下水位，可有效的改良當(dāng)?shù)氐耐寥拉h(huán)境。研究區(qū)北部與古爾班通古特沙漠接壤，接壤地帶未利用的裸地較多，而裸地含鹽量顯著高于其他土地利用方式（圖5），可在裸地選擇合適的樹種實(shí)施植樹造林，以增加土地覆蓋度降低土壤含鹽量，樹木通過改善土壤水分條件和土壤結(jié)構(gòu)，也可幫助降低附近耕地的土壤含鹽量^［28］。研究區(qū)內(nèi)海拔較低區(qū)域土壤含鹽量顯著高于高海拔地區(qū)，因此可對低海拔地區(qū)的耕地進(jìn)行深耕、秸稈還田和冬灌等措施^［29-30］。

（3）本研究由于單次取樣，不能明確研究區(qū)棉田土壤含鹽量的動(dòng)態(tài)變化，因此無法對灌溉管理的優(yōu)化提出針對性的建議。為了獲得更準(zhǔn)確、全面的土壤鹽漬化研究結(jié)果，建議采取多點(diǎn)多次取樣的方式，并結(jié)合長期監(jiān)測數(shù)據(jù)，以便于分析土壤鹽分變化的動(dòng)態(tài)特征和影響因素。這種方法不僅有助于更好地理解鹽漬化的時(shí)空變化特征，而且也能為后續(xù)的鹽漬化管理與治理提供更可靠的數(shù)據(jù)基礎(chǔ)。

4 結(jié)論

（1）研究區(qū)0～20、20～40、40～60 cm土層含鹽量均值分別為2.14、3.89、5.00 g·kg^-1，0～20 cm土層屬于非鹽化土，20～60 cm屬于輕度鹽化土，三層土壤都為中等強(qiáng)度變異，在0～20 cm土層土壤具有中等程度相關(guān)性，由隨機(jī)性因素和結(jié)構(gòu)性因素共同作用，20～60 cm土層土壤具有強(qiáng)空間相關(guān)性，由結(jié)構(gòu)性因素主導(dǎo)。

（2）研究區(qū)0～20 cm土層含鹽量熱點(diǎn)（高值區(qū)）主要集中在北部，20～40、40～60 cm熱點(diǎn)（高值區(qū)）分布相似，主要集中在北部和西部。從空間分布上看，鹽漬化土壤主要在20～40、40～60 cm土層，且以輕度鹽化土為主，分別占總面積的69.14%、58.34%，中度和重度鹽化土以條帶狀和斑狀鑲嵌于輕度鹽化土內(nèi)。

（3）研究區(qū)內(nèi)不同環(huán)境因子與各土層土壤含鹽量的關(guān)聯(lián)度依次為土壤類型gt;土地利用方式gt;海拔gt;巖性gt;坡度，土壤類型、土地利用方式和海拔是土壤含鹽量的主控因子，通過方差分析和灰色關(guān)聯(lián)度分析可推斷出裸地、海拔較低區(qū)域以及土壤類型為草甸土、潮土、鹽土和沼澤土的地區(qū)是今后土壤鹽漬化預(yù)防的重點(diǎn)區(qū)域。

參考文獻(xiàn)（References）

［1］ 李亞莉，王海江，喬江飛.索灣灌區(qū)土壤剖面鹽分空間分布特征［J］.湖北農(nóng)業(yè)科學(xué)，2016，55（10）：2488-2491.

LI Y L，WANG H J，QIAO J F.Characteristics of spatial distribution of soil profile salinity in Suowan irrigation area［J］.Hubei Agricultural Science，2016，55（10）：2488-2491.

［2］ 徐園.新疆棉花產(chǎn)業(yè)助力鄉(xiāng)村振興優(yōu)化路徑研究［J］.廣東蠶業(yè)，2023，57（11）：109-111.

XU Y.Research on the optimization path of Xinjiang cotton industry to help rural revitalization［J］.Guangdong Sericulture，2023，57（11）：109-111.

［2］ 楊未靜，虎膽·吐馬爾白，米力夏提·米那多拉.石河子墾區(qū)連作膜下滴灌棉田鹽漬化土壤鹽分空間變異性研究［J］.節(jié)水灌溉，2019，44（7）：35-40.

YANG W J，HUDAN T，MILIXIATI M.Study on spatial variability of salinity in salinized soils of cotton fields under continuous membrane drip irrigation in Shihezi reclamation area［J］.Water Saving Irrigation，2019，44（7）：35-40.

［3］ 王巍琦，楊海昌，王衛(wèi)超，等.新疆綠洲鹽漬化區(qū)域土壤質(zhì)量評價(jià)［J］.黑龍江農(nóng)業(yè)科學(xué)，2018，41（5）：26-30.

WANG W Q，YANG H C，WANG W C，et al.Soil evaluation of salinized regional soils in oasis of Xinjiang［J］.Heilongjiang Agricultural Science，2018，41（5）：26-30.

［4］ JALILA J，BALGHOUTHI M，EZZAOUIA H.Salt removal from soil using the argil porous ceramic［J］.Desalination，2016，379：53-67.

［5］ LUCA S，CARLOTTA F.The local-scale impact of soil salinization on the socioeconomic context;an exploratory analysis in italy［J］.Catena，2015，127：312-322.

［6］ DALIAKOPOULOS I N，TSANIS I K，KOUTROULIS A，et al.The threat of soil salinity;a european scale review［J］.Science of The Total Environment，2016，573：727-739.

［7］ WANG Y G，DENG C，LIU Y.Identifying change in spatial accumulation of soil salinity in an inland river watershed，China［J］.Science of The Total Environment，2018，621：177-185.

［8］ BENYAMINI Y，MIRLAS V，MARISH S，et al.A survey of soil salinity and groundwater level control systems in irrigated fields in the jezre’el valley，israel［J］.Agricultural Water Management，2005，76（3）：181-194.

［9］ HUANG F，LIU G，KEARNEY M，et al.Georelational analysis of soil type，soil salt content，landform，and land use in the yellow river delta，china［J］.Environmental Management，2005，35（1）：72-83.

［10］ 李玉義，張鳳華，潘旭東，等.新疆瑪納斯河流域不同地貌類型土壤鹽分累積變化［J］.農(nóng)業(yè)工程學(xué)報(bào)，2007，23（2）：60-64.

［12］ LI Y Y，ZHANG F H，PAN X D，et al.Changes in soil salinity accumulation in different landform types in the Manas River basin，Xinjiang［J］.Journal of Agricultural Engineering，2007，23（2）：60-64.

［11］ 宗含，高龍，王雅琴，等.膜下滴灌條件下鹽荒地土壤鹽分變化規(guī)律研究［J］.干旱地區(qū)農(nóng)業(yè)研究，2018，36（6）：7-12.

ZONG H，GAO L，WANG Y Q，et al.Research on the change rule of soil salinity in salt wasteland under the condition of drip irrigation under membrane［J］.Arid Region Agricultural Research，2018，36（6）：7-12.

［12］ 鄭琦，王海江，李萬濤，等.瑪納斯河流域土壤鹽漬化影響因素研究［J］.農(nóng)業(yè)資源與環(huán)境學(xué)報(bào)，2016，33（3）：214-220.

ZHRNG Q，WANG H J，LI W T，et al.Research on factors affecting soil salinization in Manas River Basin［J］.Journal of Agricultural Resources and Sciences，2016，33（3）：214-220.

［13］ 谷海斌，盛建東，武紅旗，等.灌區(qū)尺度土壤鹽漬化調(diào)查與評價(jià)——以石河子灌區(qū)和瑪納斯灌區(qū)為例［J］.新疆農(nóng)業(yè)大學(xué)學(xué)報(bào)，2010，33（2）：95-100.

GU H B，SHENG J D，WU H Q，et al.Survey and evaluation of soil salinization at the irrigation district scale-taking Shihezi and Manas irrigation districts as examples［J］.Journal of Xinjiang Agricultural University，2010，33（2）：95-100.

［14］ 顏安，王澤，盛建東.不同因子對瑪納斯河流域農(nóng)田土壤鹽分的影響［J］.新疆農(nóng)業(yè)大學(xué)學(xué)報(bào)，2013，36（4）：329-333.

YAN A，WANG Z，SHRNG J D.Effects of different factors on soil salinity in farmland in Manas River Basin［J］.Journal of Xinjiang Agricultural University，2013，36（4）：329-333.

［15］ 陳霖明，李艷紅，李發(fā)東，等.瑪納斯河流域出苗期棉田土壤膜下滴灌前后水分-鹽分-養(yǎng)分運(yùn)移分析［J］.甘肅農(nóng)業(yè)大學(xué)學(xué)報(bào)，2021，56（5）：110-119.

CHEN L M，LI Y H， LI F D，et al.Water-salt-nutrient transport in cotton field soils before and after membrane drip irrigation during emergence in Manas River basin［J］.Journal of Gansu Agricultural University，2021，56（5）：110-119.

［16］ 劉新，郝媛媛，花立民.石羊河下游民勤盆地土壤鹽分空間分異特征［J］.干旱區(qū)研究，2023，40（10）：1615-1624.

LIU X，HAO Y Y， HUA L M.Characteristics of spatial differentiation of soil salinity in Minqin Basin downstream of Shiyang River［J］.Arid Zone Research，2023，40（10）：1615-1624.

［19］ 付慧，劉艷軍，孫宏日，等.京津冀地區(qū)耕地利用轉(zhuǎn)型時(shí)空 分異及驅(qū)動(dòng)機(jī)制［J］． 地理科學(xué)進(jìn)展，2020，39（12）：1985-1998.

FU H，LIU Y J，SUN H R，et al.Spatio-temporal variation and driving mechanism of arable land utilization transition in Beijing-Tianjin-Hebei region［J］.Progress in Geoscience，2020，39（12）：1985-1998.

［21］ 張同娟，楊勁松，劉廣明，等.基于灰色關(guān)聯(lián)度法河口地區(qū)土壤鹽分影響因子分析［J］.土壤通報(bào)，2010，41（4）：793-796.

ZHANG T J，YANG J S，LIU G M，et al.Soil salinity influencing factors in estuaries based on the gray correlation method［J］.Soil Bulletin，2010，41（4）：793-796.

［22］ 孫波，趙其國，閭國年.低丘紅壤肥力的時(shí)空變異［J］.土壤學(xué)報(bào)，2002（2）：190-198.

SUN B，ZHAO Q G，LU G N.Spatial and temporal variability of fertility of red soil in low hills［J］.Journal of Soil Science，2002（2）：190-198.

［17］ 喬江飛，王海江，李亞莉，等.莫索灣灌區(qū)土壤剖面鹽分的季節(jié)性變化特征［J］.中國土壤與肥料，2016，53（3）：13-18.

QIAO J F，WANG H J，LI Y L，et al.Seasonal change characteristics of soil profile salinity in Moso Bay irrigation area［J］.China Soil and Fertilizer，2016，53（3）：13-18.

［18］ 吳英.松嫩平原低平易澇區(qū)土壤鹽分的季節(jié)性變化［J］.土壤，1997，40（2）：92-95.

WU Y.Seasonal changes of soil salinity in the low-flat，flood-prone area of the Songnen Plain［J］.Soil，1997，40（2）：92-95.

［19］ 陳霖明，李艷紅，李發(fā)東，等.瑪納斯河流域棉田膜下滴灌前后土壤水分-鹽分-養(yǎng)分運(yùn)移分析［J］.中國土壤與肥料，2022，59（8）：1-113.

CHEN L M，LI Y H，LI F D，et al.Soil water-salt-nutrient transport before and after membrane drip irrigation in cotton fields in the Manas River basin［J］.China Soil and Fertilizer，2022，59（8）：1-113.

［20］ 譚澤瓊，陳坤，傅軼，等.基于ArcGIS的帶狀數(shù)據(jù)插值模型對比研究［J］.海洋技術(shù)學(xué)報(bào)，2024，43（3）：46-52.

TAN Z Q，CHRN K，F(xiàn)U Y，et al.Comparative study of banded data interpolation models based on ArcGIS［J］.Journal of Marine Technology，2024，43（3）：46-52.

［21］ 張帆，錢剛.下野地灌區(qū)鹽漬化原因分析及防治對策［J］.山西水土保持科技，2005，32（2）：29-30.

ZHANG F，QIAN G.Causes of salinization and countermeasures for prevention and control of salinization in Shimanodi irrigation area［J］.Shanxi water and soil science and technology，2005，32（2）：29-30.

［22］ 管孝艷，王少麗，高占義，等.鹽漬化灌區(qū)土壤鹽分的時(shí)空變異特征及其與地下水埋深的關(guān)系［J］.生態(tài)學(xué)報(bào)，2012，32（4）：198-206.

GUAN X Y，WANG S L，GAO Z Y，et al.Characteristics of spatial and temporal variability of soil salinity in salinized irrigation areas and its relationship with groundwater burial depth［J］.Journal of Ecology，2012，32（4）：198-206.

［23］ 喬學(xué)瑾，石建初，譚麗麗，等.新疆安集海灌區(qū)膜下滴灌棉田根系層土壤鹽分多尺度空間分布特征［J］.土壤學(xué)報(bào)，2021，58（5）：1202-1213.

QIAO X J，SHI J C，TAN L L，et al.Characteristics of multi-scale spatial distribution of soil salinity in the rhizosphere of rhizosphere of drip-irrigated cotton fields under membrane in Andijhai Irrigation District，Xinjiang［J］.Journal of Soil Science， 2021，58（5）：1202-1213.

［24］ 翟江蕊，白云崗，加孜拉，等.塔河下游典型綠洲灌區(qū)不同土地利用類型土壤的鹽漬化特征［J］.水土保持通報(bào)，2023，43（3）：69-79，85.

ZHAI J R，BAI Y G， GAZIRA，et al.Salinization characteristics of soils of different land use types in a typical oasis irrigation area in the lower reaches of the Tar River［J］.Soil and Water Bulletin，2023，43（3）：69-79，85.

［25］ 王雪梅，劉瑞文.不同地形和土壤類型下精河縣耕層土壤鹽分特征分析［J］.陜西農(nóng)業(yè)科學(xué)，2014，60（3）：56-58.

WANG X M， LIU R W.Salinity characteristics of tillage soil in Jinghe County under different topography and soil types［J］.Shaanxi Agricultural Science，2014，60（3）：56-58.

［26］ 竇旭，史海濱，苗慶豐，等.鹽漬化灌區(qū)土壤水鹽時(shí)空變異特征分析及地下水埋深對鹽分的影響［J］.水土保持學(xué)報(bào)，2019，33（3）：246-253.

DOU X，SHI H B，MIAO Q F，et al.Characteristics of spatial and temporal variability of soil water salinity in salinized irrigation areas and the effect of groundwater burial depth on salinity［J］.Journal of Soil and Water，2019，33（3）：246-253.

［27］ 王國帥，史海濱，李仙岳，等.河套灌區(qū)不同地類鹽分遷移估算及與地下水埋深的關(guān)系［J］.農(nóng)業(yè)機(jī)械學(xué)報(bào)，2020，51（8）：255-269.

WANG G S，SHI H B，LI X Y，et al.Estimation of salinity migration of different land classes in the river-loop irrigation area and the relationship with groundwater burial depth［J］.Journal of Agricultural Machinery，2020，51（8）：255-269.

［28］ JENA A K.Tree plantation for implementing the land utilization policies incorporate with sustainable development towards environmental concerns in Assam，India［J］.Land Use Policy，2024，141：107085.

［29］ 鄭琦，王海江，李萬濤，等.瑪納斯河流域土壤鹽漬化影響因素研究［J］.農(nóng)業(yè)資源與環(huán)境學(xué)報(bào)，2016，33（3）：214-220.

ZHENG Q，WANG H J， Li W T，et al.Research on factors affecting soil salinization in Manas River Basin［J］.Journal of Agricultural Resources and Sciences，2016， 33（3）：214-220.

［29］ 趙波，王振華，李文昊.滴灌方式及定額對北疆冬灌棉田土壤水鹽分布及次年棉花生長的影響［J］.農(nóng)業(yè)工程學(xué)報(bào)，2016，32（6）：139-148.

ZHAO B，WANG Z H， LI W H.Effects of drip irrigation and quota on soil water and salt distribution in winter-irrigated cotton fields in Northern Xinjiang and cotton growth in the following year［J］.Journal of Agricultural Engineering，2016，32（6）：139-148.

［30］ 勞秀榮，吳子一，高燕春.長期秸稈還田改土培肥效應(yīng)的研究［J］.農(nóng)業(yè)工程學(xué)報(bào)，2002（2）：49-52.

LAO X R，WU Z Y，GAO Y C.Research on the effect of straw returned to the field for soil reformation and fertilization［J］.Journal of Agricultural Engineering，2002（2）：49-52.

（責(zé)任編輯：編輯張忠）

基金項(xiàng)目：第三次新疆綜合科學(xué)考察課題（2021xjkk0901）

作者簡介：袁裕山（1999—），男，碩士研究生，專業(yè)研究方向?yàn)楣?jié)水灌溉理論與技術(shù)，e-mail：528423221@qq.com。

*通信作者：呂德生（1963—），男，教授，主要從事建筑材料、節(jié)水灌溉等方面的研究，e-mail：13899536225@163.com