王 珊，胡振華，張寶忠，陳 鶴

(1.山西農(nóng)業(yè)大學林學院，山西 晉中 030801；2.中國水利水電科學研究院 流域水循環(huán)模擬與調(diào)控國家重點實驗室，北京 100038；3.國家節(jié)水灌溉北京工程技術(shù)研究中心，北京 100048)

0 引 言

土壤含水量的準確獲取是制定灌溉制度，實現(xiàn)對作物適時適量灌溉的基礎。然而，由于農(nóng)田土壤特性的空間異質(zhì)性，導致土壤含水量監(jiān)測點的布設數(shù)量與位置對監(jiān)測精度存在很大影響，因此亟待開展考慮空間變異特點的土壤水分優(yōu)化布點方法研究。

近年來，國內(nèi)外許多學者在土壤水分空間變異定量表征[1-5]、土壤含水量監(jiān)測點合理數(shù)目及優(yōu)化布設位置[6-11]等方面做了大量研究。在空間變異表征方面，申祥民[12]等在新疆以膜下滴灌棉田為研究對象，利用地統(tǒng)計學理論研究土壤含水量變異規(guī)律，結(jié)果表明土壤含水量在20～40和40～60 cm的空間變異為弱變異，60～80和80～100 cm層的空間變異強度屬弱變異及中等變異。M.Herbst[13]等在德國以0.29 km2田塊為研究對象，對0～20 cm土層土壤水分的空間變異性進行了研究，在1 hm2內(nèi)采用10 m×20 m網(wǎng)格取樣，結(jié)果表明其變異函數(shù)符合球狀模型，相關(guān)距約53.8 m。Wilson D J[14]等在新西蘭某流域60 hm2內(nèi)采用40 m×40 m網(wǎng)格取樣，研究了0～30 cm層土壤水分的空間分布，結(jié)果表明其變異系數(shù)為19%，變異函數(shù)符合球狀模型，空間相關(guān)域為150 m。N.V.Prasolova[15]等利用經(jīng)典統(tǒng)計學研究了澳大利亞種植夏松地塊的土壤水分及氮等元素的時空變異，結(jié)果表明采樣策略與研究的對象、采樣手段等有關(guān)。

在布設數(shù)量優(yōu)化方面，李彥[16]等在新疆石河子以滴灌棉田為對象研究發(fā)現(xiàn)應用等值線圖及地統(tǒng)計學方法，土壤含水量的采樣效率比傳統(tǒng)統(tǒng)計學方法提高約3倍。張繼光、楊貴羽、吳宏平[17-19]等研究發(fā)現(xiàn)在確定合理取樣數(shù)時應同時考慮統(tǒng)計特征和土壤的空間結(jié)構(gòu)性。孫凱[20]等在北京通州區(qū)采用經(jīng)典統(tǒng)計學和地統(tǒng)計學方法研究區(qū)域含水量的合理采樣數(shù)目。陳坤[21]等研究發(fā)現(xiàn)在120 m×90 m區(qū)域內(nèi)的膜下滴灌番茄生育期，其平均地塊全方向的變程值為23.81 m，在置信水平95%，采樣誤差5%μ的條件下，合理的采樣數(shù)目為18～20 個。

在布設位置方面，楊風亮、李芳松[22,23]等研究發(fā)現(xiàn)在棉花膜下滴灌田含水量監(jiān)測點布置在距滴灌帶47 cm范圍內(nèi)為宜。Phene[24]等研究發(fā)現(xiàn)在滴灌條件下，傳感器應該安放在靠近滴頭的區(qū)域。Stegman[25]等研究表明傳感器可以埋設在不同根系分布層，并且較深處傳感器埋設深度是較淺處傳感器的2倍。Haman[26]等研究發(fā)現(xiàn)一般將滴灌傳感器放置在作物根系活躍層，且距滴頭一定距離時，測墑精度較高。申孝軍[27]等利用R型譜系聚類法對滴灌棉田的土壤水分測點布設方案進行研究，發(fā)現(xiàn)利用徑向上距滴灌帶0、32.5和50 cm處觀測點以及垂向0～10、20～30、40～50和60～80 cm處觀測點的體積含水量可以較好地反映整個土壤剖面的含水量情況。

然而，有關(guān)土壤含水量合理采樣數(shù)目、布設位置的研究多以土壤含水量實測數(shù)據(jù)為基礎確定，每次的實測數(shù)據(jù)僅能代表這段時期的土壤含水量，多次之間可能存在較大差異，從而導致合理采樣點數(shù)量和布設不斷變化，缺乏通用性。因此，本文提出以土壤有效最大含水量作為關(guān)鍵參數(shù)，結(jié)合經(jīng)典統(tǒng)計方法和地統(tǒng)計方法，構(gòu)建土壤含水量監(jiān)測優(yōu)化布設方法。

1 材料與方法

1.1 研究區(qū)概況

試驗在國家節(jié)水灌溉北京工程技術(shù)研究中心大興試驗基地(N39°37.25′，E116°25.51′)周邊3.645 km2試驗區(qū)開展。該試驗區(qū)位于北京市大興區(qū)，屬于半干旱大陸性暖濕季風氣候，這里四季分明，春季少雨多風，夏季炎熱多雨，秋季天高氣爽，冬季寒冷少雨，多年平均氣溫為12.1 ℃，極端溫度最高為39.5 ℃，最低為-25 ℃。區(qū)內(nèi)年際間降雨量變化較大，降雨量最大為971 mm，最小為206 mm，年內(nèi)降雨量分布不均，多年平均降雨為540 mm，降雨多集中在汛期的6-9月份，其降雨量占全年降雨量的80%以上[28]。該區(qū)光熱自然條件豐富，適合多種作物生長，包括小麥、玉米等。

1.2 測定指標與方法

將3.645 km2試驗區(qū)劃分為162個小區(qū)，每個小區(qū)范圍為150 m×150 m。根據(jù)種植結(jié)構(gòu)，本文選擇玉米種植的129個小區(qū)開展研究，每個小區(qū)的采樣點按照深度0～20、20～40、40～60、60～80 cm共分4層取樣。取0～20和20～40 cm的平均值作為耕作層(0～40 cm)數(shù)據(jù)，取0～20、20～40、40～60和60～80 cm的平均值作為作物根區(qū)(0～80 cm)數(shù)據(jù)。

1.2.1 土壤有效最大含水量

采用馬爾文激光粒度儀(MS2000)對129個小區(qū)的分層土樣進行顆粒分級分析，根據(jù)土樣中各粒徑顆粒的百分含量后，利用美國土壤質(zhì)地分類三角表進行土壤質(zhì)地類型的確定，并采用SPAW[]軟件，獲得土壤有效最大含水量。

wa=wf-wp

(1)

式中：wa為土壤有效最大含水量，%；wf為田間持水量，%；wp為凋萎系數(shù)，%。

1.2.2 土壤含水量

采用烘干法測定129個小區(qū)樣點的分層土壤重量含水量，25 d左右測定一次，共測5次。

1.3 土壤含水量監(jiān)測優(yōu)化布設方法

以土壤有效最大含水量為關(guān)鍵參數(shù)，考慮其統(tǒng)計特征和空間結(jié)構(gòu)性，采用經(jīng)典統(tǒng)計學方法確定某置信度的合理布設數(shù)目，根據(jù)地統(tǒng)計學方法的變程分析確定取樣點布置。

1.3.1 合理布設數(shù)目

根據(jù)經(jīng)典統(tǒng)計學方法，在樣本數(shù)量較多的情況下，可以用樣本的數(shù)學期望和方差代替總體的數(shù)學期望和方差，計算公式如下[30]：

N=(λa·S)2/Δ2

(2)

式中：λα為正態(tài)分布特征值；S為樣本標準方差；△為采樣誤差，通常取Δ=kμ，k可取5%、10%、15%、20%等。

如果計算所得出的合理采樣數(shù)大于樣本容量的10%，則應將計算公式轉(zhuǎn)換為不等重復抽樣公式再計算抽樣數(shù)N′：

(3)

式中：N總為樣本總數(shù)；N為等重復抽樣時計算的抽樣數(shù)。

1.3.2 合理布設位置

根據(jù)地統(tǒng)計分析方法，采用半方差函數(shù)(空間變異函數(shù))描述土壤特性空間變異結(jié)構(gòu)，計算公式如下：

(4)

式中：γ(h)為區(qū)域化變量Z(X)的樣本變異函數(shù)；h為兩樣本點空間分隔距離；Z(Xi)和Z(Xi+h)分別是區(qū)域化變量Z(X)在空間位置Xi和(Xi+h)的觀測值；N(h)為間距為的樣本對數(shù)。本文采用交叉驗證法，根據(jù)誤差最小原則確定最優(yōu)半方差函數(shù)模型。

應用surfer12.0軟件對土壤有效最大含水量進行空間插值，獲得土壤有效最大含水量等值線圖與平均值，計算每一個插值點處的含水率與平均含水率的相對誤差，選擇平均相對誤差最小的插值點即為最優(yōu)測墑點，根據(jù)半方差函數(shù)模型參數(shù)變程為基準做出扇形區(qū)域與均值線相交選擇樣點，確定采樣點的具體位置。

2 結(jié)果與討論

2.1 土壤有效最大含水量的空間變異特點

2.1.1 正態(tài)分布檢驗

圖1是0～40與0～80 cm土壤有效最大含水量的正態(tài)分布檢驗圖，可以看出，0～40、0～80 cm深度有效最大含水量數(shù)據(jù)絕大部分都在45°直線附近，說明它們均服從正態(tài)分布。因此，可以應用經(jīng)典統(tǒng)計學方法分析空間變異和確定合理布設數(shù)目。這與之前一些研究有些相似，陳坤[31]和張繼光[17]都在研究土壤含水量變異特征過程中發(fā)現(xiàn)土壤含水量需滿足正態(tài)分布。

圖1 土壤有效最大含水量正態(tài)分布檢驗圖(單位：%)Fig.1 Normal test of soil available maximum water content

2.1.2 土壤有效最大含水量空間變異

表1為不同深度田間土壤有效最大含水量統(tǒng)計特征值，表明0～40和0～80 cm深度的土壤有效最大含水量的變異系數(shù)分別是14.96%和13.56%，屬于中等程度變異。

圖2為0～40、0～80 cm土壤有效最大含水量的空間分布圖，從圖2中可以看出，0～40、0～80 cm的空間分布情況基本一致，在西南-西北方向的土壤有效最大含水量較大，其他地區(qū)較??；0～40 cm的土壤有效最大含水量在10%～23%之間變化，0～80 cm深度的土壤有效最大含水量在12%～23%之間變化。

表1 不同深度田間土壤有效最大含水量統(tǒng)計特征值Tab.1 Statistic characters of soil available maximum water content of the different depth

圖2 土壤有效最大含水量等值線圖(單位：%)Fig.2 Contour map of soil available maximum water content

2.2 土壤含水量監(jiān)測優(yōu)化布設數(shù)量與位置

2.2.1 半方差函數(shù)理論模型

采用交叉驗證法檢驗選擇適合的半方差函數(shù)理論模型，具體結(jié)果如表2和圖3所示。由表2可以看出，0～40與0～80 cm深度的半方差函數(shù)理論模型都為指數(shù)形式，且0～40與0～80 cm深度的塊金基臺比分別是25.56%、24.51%，0～40 cm表現(xiàn)為中等程度的空間相關(guān)性，受隨機性因素與結(jié)構(gòu)性因素的共同作用影響，0～80 cm表現(xiàn)為較強的空間相關(guān)性，主要受結(jié)構(gòu)性因素的影響。0～40 cm深度變程是780 m，0～80 cm深度變程是1 104 m，在變程范圍內(nèi)的監(jiān)測點之間具有相關(guān)性，在變程范圍外的監(jiān)測點是相互獨立。由圖3可看出0～40 cm的基臺值大于0～80 cm的基臺值，0～40與0～80 cm的變異函數(shù)值均隨著滯后距離的增加先增大后趨于平穩(wěn)狀態(tài)。

表2 土壤有效最大含水量半方差函數(shù)模型參數(shù)Tab.2 Semivariogram model parameter of soil available maximum water content

圖3 土壤有效最大含水量半方差函數(shù)模型Fig.3 Semivariogram model of soil available maximum water content

2.2.2 合理采樣數(shù)目與位置

根據(jù)經(jīng)典統(tǒng)計方法，取置信區(qū)間95%、90%，采樣誤差取5%、10%、15%、20%，根據(jù)式(2)和式(3)確定0～40和0～80 cm深度的合理采樣數(shù)目見表3所示。在同一置信水平情況下，隨著采樣誤差的增大而采樣數(shù)目減少；在一定的采樣誤差條件下，合理采樣數(shù)目也隨著置信水平的降低而逐漸減少。楊詩秀以133 hm2農(nóng)田為例，研究也表明在置信度p=80%，誤差Δ=0.1的條件下，取樣數(shù)目為10，將采樣誤差降為0.2時，取樣數(shù)目為4[32]。

表3 不同置信水平、不同采樣誤差條件下的合理采樣數(shù)目Tab.3 The reasonable number of samples under the different confidence level and different sampling error

通過129個樣點的土壤有效最大含水量進行空間插值，獲得0～40、0～80 cm的空間分布圖，根據(jù)空間插值獲得每個插值點的土壤有效最大含水量值與區(qū)域平均值，計算它們的相對誤差，選擇最小相對誤差為最優(yōu)測墑點，并以半方差函數(shù)模型參數(shù)變程為基準，在均值線上確定其他采樣點的具體位置。試驗區(qū)的土壤有效最大含水量在90%置信區(qū)間，10%的采樣誤差水平下，0～40 cm深度合理采樣點數(shù)目為6個，0～80 cm深度合理采樣點數(shù)目為5個，具體坐標如表4和圖4所示。

表4 不同深度合理采樣點坐標統(tǒng)計Tab.4 Statistics of reasonable sampling point coordinates at different depths

圖4 土壤有效最大含水量的空間分布(單位：%)Fig.4 Spatial distribution of soil available maximum water content

2.3 誤差分析

以129個樣點插值獲得的區(qū)域含水量作為實測值，優(yōu)化布設的采樣點插值獲得的區(qū)域含水量作為估算值，將其對應的土壤含水量平均值進行誤差分析。具體分析如表5、表6所示。

表5是0～40 cm土壤平均含水量估算值與實測值的誤差分析結(jié)果。在0～40 cm深度土壤含水量的誤差值在1.56%～8.95%范圍內(nèi)變化，均在10%的采樣誤差范圍內(nèi)，說明0～40 cm深度的6個固定監(jiān)測點滿足采樣精度，能夠代表區(qū)域的土壤含水量。

表5 0～40 cm土壤含水量誤差分析Tab.5 Analysis of soil water content error at 0～40 cm depth

表6 0～80 cm土壤含水量誤差分析Tab.6 Analysis of soil water content error at 0～80 cm depth

表6是0～80 cm土壤平均含水量估算值與實測值的誤差分析結(jié)果。在0～80 cm深度土壤含水量的誤差值在0.56%～5.43%范圍內(nèi)變化，均在10%的采樣誤差范圍內(nèi)，說明0～80 cm深度的5個固定監(jiān)測點滿足采樣精度，能夠代表區(qū)域的土壤含水量。

3 結(jié) 語

(1)試驗區(qū)土壤有效最大含水量在0～40、 0～80 cm深度的變異程度均表現(xiàn)為中等程度，0～40、0～80 cm的土壤有效最大含水量的空間分布情況基本一致，在西南-西北方向的土壤有效最大含水量較大，其他地區(qū)較??；且0～40 cm的變化較0～80 cm大。

(2)0～40 cm的土壤有效最大含水量具有中等程度的空間相關(guān)性，0～80 cm的土壤有效最大含水量具有較強的空間相關(guān)性，相關(guān)域分別是780、1 104 m，選擇的理論模型均為指數(shù)模型。在置信區(qū)間90%，采樣誤差10%的條件下，0～40與0～80 cm的合理采樣數(shù)分別是6個、5個，在土壤有效最大含水量等值線圖中確定具體坐標。

(3)根據(jù)試驗區(qū)129個樣點估算的區(qū)域土壤含水量與優(yōu)化布設的采樣點估算的區(qū)域含水量進行誤差分析，結(jié)果表明均滿足精度要求，說明土壤有效最大含水量可以代替實測含水量值進行測墑點的優(yōu)化布設。

□

參考文獻：

[1] 申祥民,雷曉云.棉花膜下滴灌條件下農(nóng)田土壤水分的空間變異性研究[D]. 烏魯木齊：新疆農(nóng)業(yè)大學,2010.

[2] Brocca L, Morbidelli R, Melone F, et al. Soil Moisture Spatial Variability in Experimental Areas of Central Italy[J].Journal of Hydrology(0022-1694). 2007,333:356-373.

[3] 周世軍,雷曉云,李芳松,等. 棉花膜下滴灌田間土壤水分的時空變異規(guī)律研究[J].中國農(nóng)村水利水電,2010,(10):22-24.

[4] Herbst M, Diekkruger B. Modeling the spatial variability of soil moisture in a micro-scale catchment and comparison with field data using geostatistics[J]. Phys. Chem. Earth 2003,28(6-7):239-245.

[5] 李芳松,雷曉云,周世軍,等.南疆地區(qū)棉花整個生育期內(nèi)土壤水分空間變異性分析[J].水資源與水工程學報,2010,21(2):63-67.

[6] 龔紹琦, 黃家柱, 李云梅, 等. 太湖梅梁灣水質(zhì)參數(shù)空間變異及合理取樣數(shù)目研究[J]. 地理與地理信息科學, 2006,22(2): 50-54.

[7] 謝恒星,張振華,劉繼龍,等.土壤含水量合理取樣數(shù)目影響因素的試驗研究----以煙臺蘋果園為例[J].干旱地區(qū)農(nóng)業(yè)研究, 2007,25(4):114-118.

[8] 信秀麗,張佳寶,朱安寧.區(qū)域土壤水分監(jiān)測點布設方式探討[J].灌溉排水學報,2008,27(3):16-19.

[9] 汪羽寧,樊 軍,李世清.小麥實時控制灌溉的土壤水分含量探頭合理埋設深度研究[J].灌溉排水學報,2009,28(5):10-12.

[10] Junior V V, Carvalho. Spatial variability of soil water content and mechanical resistance of Bra Zilian ferralsol[J].Soil and Tillage Research,2006,85 (1-2):166-177．

[11] Ferreyra R A,Apeztegufa H P,Sereno R. Reduction of soil water spatial sampling density using scaled semivariograms and simulated annealing[J].Geodema,2002,110(3-4),265-289.

[12] 申祥民,雷曉云,陳大春,等.不同布點方式的膜下滴灌棉田土壤水分的空間變異研究[J].新疆農(nóng)業(yè)大學學報,2010,33(4):363-368.

[13] Lydia Mumbi CHABALA，Augustine MULOLWA，Obed LUNGU. Application of Ordinary Kriging in Mapping Soil Organic Carbon in Zambia[J].Pedosphere,2017,27(2):338-343.

[14] Wilson D J,Western A W, Rodger B Grayson, et al. Spatial distribution of soil moisture over 6 and 30 cm depth, Mahurangi river catchment, New Zealand[J]. Journal of Hydrology, 2003,276(1-4):254-274.

[15] N.V.Prasolova. Spatial-temporal variability of soil moisture,nitrogen availability indices and other chemical properties in hoop pine (Araucaria cunninghamii) plantations of subtropical Australia[J].Forest Ecology and Management, 2000,136(1-3):1-10.

[16] 李 彥,鄭國玉,雷曉云.滴灌棉田灌溉單元內(nèi)土壤墑情監(jiān)測點的采樣效率研究[J].灌溉排水學報.2014,33(2):90-93.

[17] 張繼光,陳洪松,蘇以榮,等. 喀斯特地區(qū)典型峰叢洼地表層土壤水分空間變異及合理取樣數(shù)研究 [J]. 水土保持學報, 2006,20(2):114-117.

[18] 楊貴羽,陳亞新. 土壤水分鹽分空間變異性與合理采樣數(shù)研究[J].干旱地區(qū)農(nóng)業(yè)研究, 2002,20(4):64-66.

[19] 吳宏平,王樹仿,繳錫云. 膜下滴灌墑情監(jiān)測點布設方案研究[J].水利與建筑工程學報,2015,13(2):144-147.

[20] 孫 凱,王一鳴,楊紹輝.墑情監(jiān)測取樣方法的研究[J].農(nóng)業(yè)工程學報,2004,20(4):74-78.

[21] 陳 坤,雷曉云,李 彥,等.加工番茄膜下滴灌墑情監(jiān)測點的數(shù)目研究[J].節(jié)水灌溉，2012,(1)：4-7.

[22] 楊風亮,繳錫云,劉一休,等.棉花膜下滴灌墑情監(jiān)測點的定位[J].灌溉排水學報,2010,29(3):29-31.

[23] 李芳松,雷曉云,周世軍,等.膜下滴灌棉田墑情監(jiān)測點的橫向定位研究[J].節(jié)水灌溉,2011,(3):4-6.

[24] Phene C J, Howell T A. Soil sensor control of high-frequency irrigation systems[J]. Transaction of the ASAE, 1984,27(2):102-104.

[25] Stegman E C. Irrigation Water Management[M].Chicago, Illinois, USA :American Society of Agricultural Engineers, 1982:763-816.

[26] Haman D Z,Ross Phil, Tacker D S. Critical management issues for SDI systems in North Carolina[J].North Carolina Cooperative Extension Service,2008,(6):45-56.

[27] 申孝軍,孫景生,張寄陽,等.滴灌棉田土壤水分測點最優(yōu)布設研究[J].干旱地區(qū)農(nóng)業(yè)研究,2012,30(3):90-95.

[28] 劉國水,劉 鈺,許 迪,等.農(nóng)田不同尺度蒸散量的尺度效應與氣象因子的關(guān)系[J].水利學報.2011,42(3):284-289.

[29] 關(guān)作正,石 嵐,徐麗娜.黃河流域托萬區(qū)間SWAT模型土壤屬性數(shù)據(jù)庫的建立[J].內(nèi)蒙古水利.2009,(4):5-7.

[30] Cochran W G. Sampling Techniques[M]. Third Edition. New York: John Wiley and Sons, Inc., 1977.

[31] 陳 坤,雷曉云.自動化膜下滴灌加工番茄墑情特征分析及墑情預報模型研究[D]. 烏魯木齊：新疆農(nóng)業(yè)大學,2012.

[32] 楊詩秀,雷志棟,吳婉如,等.農(nóng)田尺度土壤水分的監(jiān)測[J].水科學進展.1996,7(1):14-19.