張帆　王劍秦　鄭立華　盧海達(dá)

摘 要：作物所需的無(wú)機(jī)營(yíng)養(yǎng)絕大多數(shù)存在于土壤水分中，精確描述土壤水分在作物根系層的分布情況有利于監(jiān)測(cè)作物的生長(zhǎng)壞境。隨著農(nóng)業(yè)物聯(lián)網(wǎng)技術(shù)在土壤水分監(jiān)測(cè)方面應(yīng)用日趨成熟，土壤含水量在監(jiān)測(cè)點(diǎn)處數(shù)據(jù)采集較為方便，但整個(gè)土壤體中水分分布狀況無(wú)法直接獲取。本文利用克里金插值算法，結(jié)合監(jiān)測(cè)點(diǎn)所采集的土壤水分?jǐn)?shù)據(jù)，建立土壤水分分布的三維可視化模型，并通過(guò)采集的土壤含數(shù)量數(shù)據(jù)對(duì)該模型進(jìn)行了驗(yàn)證。實(shí)驗(yàn)結(jié)果表明，該模型可有效表征土壤體中水分的空間分布狀況。

關(guān)鍵詞：克里金插值算法 土壤水分建模 土壤水分分布

中圖分類號(hào)：TP311 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1672-3791（2013）03（c）-0138-02

土壤水分不僅影響土壤的物理性質(zhì)，關(guān)系著土壤中無(wú)機(jī)養(yǎng)分的溶解、轉(zhuǎn)移和微生物的活動(dòng)，而且也是作物對(duì)無(wú)機(jī)營(yíng)養(yǎng)的獲取的重要渠道，是作物生長(zhǎng)發(fā)育的基本條件。不同含水量的土壤層中，無(wú)機(jī)營(yíng)養(yǎng)（如NH4-N、NO3-N、O lesn-P等）的分布及含量也不盡相同。目前，不同學(xué)者在研究區(qū)域內(nèi)的土壤含水量分布預(yù)測(cè)方面做了大量的研究工作，克里金法在土壤含水量中的預(yù)測(cè)亦受到廣泛運(yùn)用，其合理性也在多篇文獻(xiàn)中獲得證實(shí)。本文運(yùn)用角點(diǎn)網(wǎng)格技術(shù)建立出土壤體的三維空間模型，并以此為基礎(chǔ)，利用克里金算法對(duì)各個(gè)土壤層進(jìn)行插值，得到連續(xù)的土壤水分?jǐn)?shù)據(jù)，從而建立土壤水分分布的三維模型。

1 理論與數(shù)據(jù)

1.1 克里金算法原理

克里金法是當(dāng)前空間插值方法中較為常用的方法之一，是由南非探礦學(xué)家D.Krige于1956年提出的一種經(jīng)驗(yàn)的方差，并經(jīng)法國(guó)科學(xué)家Matheson在理論上加以發(fā)展而來(lái)的。20世紀(jì)80年代，科學(xué)家Burgess將克里金法運(yùn)用到土壤科學(xué)中，后經(jīng)過(guò)多位科學(xué)家的理論論證及大量的實(shí)驗(yàn)論證，使其的可靠性也大幅增加。

假設(shè)xp（非實(shí)測(cè)點(diǎn)）為待估區(qū)域內(nèi)的一點(diǎn)，xi（i=1，2，…，n）為待估區(qū)域內(nèi)的n個(gè)實(shí)測(cè)點(diǎn)，其測(cè)量值表示為Z*（xi）（i=1，2，…，n）。則Z*（xp）。在xp處對(duì)Z（xp）（xp處的真實(shí)值）預(yù)測(cè)值可以通過(guò)n個(gè)有效實(shí)測(cè)點(diǎn)的線性組合表示，即

式中，λi表示Z*（xi）的權(quán)重因素。

克里金插值算法在理論上要求λi符合兩個(gè)條件：（1）計(jì)算結(jié)果必須滿足無(wú)偏估計(jì)；（2）估計(jì)方差必須最小。通過(guò)這兩點(diǎn)作為其理論約束條件，便能獲得克里金方程組。表述為矩陣形式即為：[K][λ]=[M]。式中，矩陣K和矩陣M與克里金法中重要的分析工具（變差函數(shù)）有關(guān)。

變差函數(shù)γ（h）本質(zhì)上僅僅是關(guān)于任意兩點(diǎn)之間距離h的函數(shù)，一般是通過(guò)γ（h）和h描繪成圖（即實(shí)驗(yàn)變差函數(shù)圖），再用數(shù)學(xué)公式擬合，就能夠得到理想的變差函數(shù)模型。

1.2 實(shí)驗(yàn)設(shè)計(jì)及數(shù)據(jù)

數(shù)據(jù)來(lái)源于于華北平原（2012年5月采集）玉米實(shí)驗(yàn)田，實(shí)驗(yàn)田長(zhǎng)約40，寬約40 m，總占地面積約1600 m2，其土壤顆粒主要由砂粒、粉粒和粘粒組成，透水性較好。主要包含兩個(gè)玉米品種（鄭單958和先玉335），玉米株距30 cm，行距70 cm，生育期大約在110天左右。土壤水分的數(shù)據(jù)是通過(guò)若干組監(jiān)測(cè)點(diǎn)獲得的。每組監(jiān)測(cè)點(diǎn)按測(cè)點(diǎn)深度排列共分四處，每處位于不同的土壤層中。第一處位于0～30 cm處（表層）；第二個(gè)處位于30～60 cm處（上層）；第三處位于60～90 cm處（下層）；最后一處位于90～120 cm處（底層）。檢測(cè)點(diǎn)各個(gè)的深度設(shè)計(jì)主要是依據(jù)該玉米品種的根莖整體深度，一般微微大于玉米根莖的整體深度。

通過(guò)獲得25組監(jiān)測(cè)點(diǎn)的數(shù)據(jù)（每組數(shù)據(jù)均包含表層、上層、下層和地層4組）分析，土壤含水量整體走勢(shì)平緩，服從正態(tài)分布。表層的土壤含水量數(shù)值上較小，而上層與下層數(shù)據(jù)整體較為近似，底層數(shù)據(jù)略大，基本符合當(dāng)日的農(nóng)田土壤水分特征。本文利用監(jiān)測(cè)點(diǎn)的采樣數(shù)據(jù)，進(jìn)行變差函數(shù)實(shí)驗(yàn)圖描繪，得出采樣數(shù)據(jù)的變差函數(shù)模型近似球狀模型（變程值：43.2，臺(tái)基值：2.1，塊金值：0.2）。因此本文選用該理論模型作為本次實(shí)驗(yàn)的克里金方法的分析工具。

2 三維土壤水分建模

當(dāng)前較為成熟的構(gòu)模方法約有20余種，大致可分為面模型、體模型和混合模型，本文依據(jù)土壤層的特點(diǎn)，分析比較網(wǎng)格各自的優(yōu)缺點(diǎn)，選用角點(diǎn)網(wǎng)格模型。

2.1 土壤空間結(jié)構(gòu)建模

在圈定的水平方向上建模范圍的基礎(chǔ)上，利用三維建模面模型（角點(diǎn)網(wǎng)格）技術(shù)，結(jié)合研究區(qū)域土壤水分建模精度要求，形成具有合理步長(zhǎng)（包括X軸方向和Y軸方向）的角點(diǎn)網(wǎng)格層，并通過(guò)土壤層中的一些離散點(diǎn)的位置信息，對(duì)角點(diǎn)網(wǎng)格中的每個(gè)網(wǎng)格點(diǎn)進(jìn)行合理的調(diào)整，使其能夠較為準(zhǔn)確地描述和逼近土壤層面的空間展布。此外，在同一次的土壤空間結(jié)構(gòu)建模中，由于水平方向上建模區(qū)域和土壤水分建模精度要求是一樣的，因此在建立出的多個(gè)離散的角點(diǎn)網(wǎng)格層面的空間拓?fù)潢P(guān)系是相同的。由于每層空間拓?fù)潢P(guān)系是一致的，所以上、下兩層中的網(wǎng)格節(jié)點(diǎn)關(guān)系是一一對(duì)應(yīng)的，利用這種關(guān)系，建立至上而下的空間結(jié)構(gòu)，形成土壤空間結(jié)構(gòu)模型。如圖1所示。

2.2 土壤水分分布建模

在土壤空間結(jié)構(gòu)模型中，上層面與下層面的節(jié)點(diǎn)建立了一一對(duì)應(yīng)關(guān)系。因此，在兩個(gè)層面之間的土壤層可以劃分出了若干組六面體（Grid）。土壤水分分布建模旨在通過(guò)克里金算法，計(jì)算出每個(gè)Grid含水量的預(yù)測(cè)值。由于克里金法是計(jì)算具有具體坐標(biāo)的某個(gè)點(diǎn)的值，因此在插值的前期，要對(duì)數(shù)據(jù)進(jìn)行一系列的預(yù)處理，包括土壤水分?jǐn)?shù)據(jù)的粗化處理，Grid內(nèi)部的加權(quán)平均值的計(jì)算等等。圖2為以土壤含水量為插值屬性所形成的分布示意圖。

2.3 三維土壤水分建模流程

三維土壤水分建模流程圖如圖3所示。

（1）建模區(qū)域信息指用戶所指定的土壤水分監(jiān)測(cè)范圍內(nèi)的土壤層信息。將離散的土壤層相對(duì)高程信息組合成角點(diǎn)網(wǎng)格，作為土壤層的架構(gòu)信息，為②中的擬合算法提供土壤的空間結(jié)構(gòu)；（2）監(jiān)測(cè)點(diǎn)的位置信息包含監(jiān)測(cè)點(diǎn)方位數(shù)據(jù)，即監(jiān)測(cè)點(diǎn)的坐標(biāo)（X、Y、Z）物理位置。在層次處理上，將每組監(jiān)測(cè)點(diǎn)的四處物理位置作為土壤層的分界線。通過(guò)對(duì)每層的角點(diǎn)網(wǎng)格擬合，并對(duì)上下相鄰的角點(diǎn)網(wǎng)格的網(wǎng)格點(diǎn)建立對(duì)應(yīng)關(guān)系，形成土壤信息的三維空間結(jié)構(gòu)，即土壤空間結(jié)構(gòu)模型；（3）土壤空間結(jié)構(gòu)模型主要為土壤水分模型提供空間架構(gòu)模型支持，包括土壤空間的拓?fù)潢P(guān)系、土壤層與土壤層之間的分界面的信息、土壤層內(nèi)部Grid分布情況等等；（4）對(duì)試驗(yàn)田中所采集的數(shù)據(jù)，在土壤空間結(jié)構(gòu)模型的基礎(chǔ)上，實(shí)現(xiàn)數(shù)據(jù)區(qū)域上的粗化，即將物理位置處在模型的某個(gè)Grid中的所有采樣點(diǎn)數(shù)據(jù)作加權(quán)平均處理，并將結(jié)果存放與該Grid的數(shù)據(jù)結(jié)構(gòu)中，作為整個(gè)Grid的屬性表現(xiàn)；（5）利用克里金算法，分別對(duì)每個(gè)土壤層中所有Grid的土壤水分進(jìn)行插值預(yù)測(cè)，預(yù)測(cè)結(jié)果存放在Grid的數(shù)據(jù)結(jié)構(gòu)中，作為該Grid的土壤含水量數(shù)值，以此形成以土壤含水量為建模要求的屬性模型；（6）顯示土壤水分分布的三維模型。

2.4 結(jié)果展示

依據(jù)圖4、圖5設(shè)計(jì)，建立以土壤含水量為插值屬性的三維土壤水分分布模型。

實(shí)驗(yàn)表明，該方案能達(dá)到監(jiān)測(cè)實(shí)驗(yàn)田各個(gè)層次的土壤含水量分布情況的要求。同時(shí)土壤水分分布模型查詢具體層面，具體位置的土壤含水量信息，也可以再次利用監(jiān)測(cè)點(diǎn)獲得新的采樣數(shù)據(jù)進(jìn)行實(shí)時(shí)更新，從而進(jìn)一步保證了模型的實(shí)時(shí)性和可持續(xù)性。

3 可靠性分析

為了比較測(cè)試結(jié)果的可靠性，我們?cè)谠囼?yàn)田中隨機(jī)采樣了7組數(shù)據(jù)，作為實(shí)驗(yàn)驗(yàn)證數(shù)據(jù)，并從已經(jīng)建立好的土壤水分分布模型中獲取與這7組數(shù)據(jù)坐標(biāo)相近的7組對(duì)應(yīng)數(shù)據(jù)，所有的數(shù)據(jù)數(shù)值如表1所示。

為了驗(yàn)證模型數(shù)據(jù)的可靠性，本文對(duì)測(cè)量值與預(yù)測(cè)值之間相對(duì)誤差進(jìn)行統(tǒng)計(jì)分析，其中相對(duì)誤差公式為：

分別對(duì)上述7組數(shù)據(jù)進(jìn)行誤差統(tǒng)計(jì)，結(jié)果如圖6所示。

如統(tǒng)計(jì)圖6所示，絕大部分的相對(duì)誤差基本控制在15%以下，能夠符合農(nóng)業(yè)生產(chǎn)的需求，達(dá)到了本次研究的目的。

4 結(jié)論

本文通過(guò)克里金內(nèi)插算法，結(jié)合土壤水分分布空間統(tǒng)計(jì)特征建立三維土壤水分分布模型，相對(duì)于過(guò)去傳統(tǒng)的預(yù)測(cè)手段，實(shí)現(xiàn)了土壤水分縱向分布的模型預(yù)測(cè)，從而進(jìn)一步細(xì)化描述了土壤水分的分布情況，具有更強(qiáng)的實(shí)用性，符合現(xiàn)代精細(xì)化節(jié)水灌溉的科學(xué)理念。通過(guò)華北平原測(cè)試數(shù)據(jù)驗(yàn)證，模型達(dá)到了農(nóng)業(yè)生產(chǎn)與灌溉的要求。土壤水分預(yù)測(cè)與三維可視化建模技術(shù)結(jié)合，對(duì)于直觀的研究和實(shí)時(shí)觀測(cè)土壤水分空間分布規(guī)律具有一定的實(shí)用價(jià)值。

參考文獻(xiàn)

[1]劉晶淼，安順清，廖榮偉.玉米根系在土壤剖面中的分布研究[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，2009，17（3）：517-521.