吳 勇，鐘永紅，杜 森，姚 麗，岳煥芳，孟范玉

（1.全國農(nóng)業(yè)技術(shù)推廣服務(wù)中心，北京100125；2.農(nóng)芯科技（北京）有限責(zé)任公司，北京100097；3.北京農(nóng)業(yè)技術(shù)推廣站，北京100125；4.農(nóng)業(yè)部農(nóng)業(yè)信息軟硬件產(chǎn)品質(zhì)量檢測重點(diǎn)實(shí)驗(yàn)室，北京100125）

0 引 言

我國是世界上13個(gè)貧水國之一，人均水資源占有量僅為世界人均水平的1/4。與此同時(shí)，我國是一個(gè)農(nóng)業(yè)大國，水資源稀缺很大程度上阻礙了農(nóng)業(yè)的健康發(fā)展，要實(shí)現(xiàn)農(nóng)業(yè)的可持續(xù)發(fā)展，應(yīng)該合理利用水資源[1]。根據(jù)不同作物對水分的需求，選擇不同的灌溉方式，可以實(shí)現(xiàn)對水分的高效利用[2，3]。因此，準(zhǔn)確了解土壤含水量，能夠?yàn)樽魑镞x擇恰當(dāng)?shù)墓喔确绞?，有效?jié)約水資源。

在土壤含水量的測量方法中，烘干法是最簡單、應(yīng)用最廣泛的方法。由于該方法是一種直接測量方法，因此可作為其他方法的參考標(biāo)準(zhǔn)[4]。采用烘干法時(shí)，通常采用環(huán)刀等工具在野外采集土樣，然后在實(shí)驗(yàn)室進(jìn)行烘干、稱重和計(jì)算。這種方法不方便取樣和儲存，且浪費(fèi)時(shí)間和精力，效率較低[5]。另外，采樣對土壤具有破壞性，長期測量需要不斷改變土壤采樣點(diǎn)，不利于長期連續(xù)定點(diǎn)監(jiān)測。土壤本身的變異性也會導(dǎo)致測量結(jié)果的差異。目前，隨著科學(xué)技術(shù)的發(fā)展，傳感器技術(shù)應(yīng)運(yùn)而生，使得土壤傳感器的應(yīng)用越來越廣泛。與烘干法相比，土壤水分傳感器具有在同一地點(diǎn)長期連續(xù)觀測含水量、安裝方便、對土壤結(jié)構(gòu)無損傷等優(yōu)點(diǎn)[6]，因此土壤水分傳感器代替烘干法測定土壤水分含水量成為了可能[7]。

市面上存在多種類型的土壤水分傳感器，由于測試原理、設(shè)計(jì)結(jié)構(gòu)的不同造成了不同類型土壤水分傳感器性能的差異，國內(nèi)外專家學(xué)者在不同類型土壤水分傳感器性能測試方面開展了許多研究。舒素芳等對FDR 土壤水分傳感器與人工土壤相對濕度觀測資料進(jìn)行了質(zhì)量對比分析，發(fā)現(xiàn)FDR 傳感器對土壤水分變化的敏感程度較低并且傳感器測定數(shù)據(jù)與人工測定數(shù)據(jù)相比普遍偏小[8]。吳日峰等對國產(chǎn)FDR 傳感器和美國的FDR 傳感器進(jìn)行室內(nèi)外同步觀測試驗(yàn)，結(jié)果表明在土壤水分含量較低時(shí)，國產(chǎn)FDR 傳感器的準(zhǔn)確性優(yōu)于美國FDR 傳感器[9]。葉志杰等對EC?5 和5TE 兩種FDR 傳感器進(jìn)行了性能測試與校正，認(rèn)為EC?5和5TE的電學(xué)性質(zhì)對溫度、電導(dǎo)率和pH值變化敏感，并給出了溫度效應(yīng)、電導(dǎo)率效應(yīng)和pH 效應(yīng)的補(bǔ)償模型[10]。Heidi Mittelbach 等對4種不同原理的土壤水分傳感器在瑞士野外進(jìn)行性能測試實(shí)驗(yàn)，發(fā)現(xiàn)在所研究的條件下，4種傳感器均不符合各自的制造商規(guī)范，需要針對特定場地重新校準(zhǔn)[11]。隨著物聯(lián)網(wǎng)技術(shù)的快速發(fā)展，目前國內(nèi)大力推廣傳感器測墑技術(shù)，國內(nèi)外學(xué)者研究主要以進(jìn)口設(shè)備為主，結(jié)構(gòu)均為探針式傳感器，探管式傳感器具有安裝方便等優(yōu)點(diǎn)，成為國內(nèi)應(yīng)用主流。因此針對國內(nèi)大力推廣測墑技術(shù)，有必要開展FDR探管式傳感器的研究?！?020年數(shù)字鄉(xiāng)村發(fā)展工作要點(diǎn)》中提出強(qiáng)化土壤墑情監(jiān)測工作，國內(nèi)FDR 傳感器應(yīng)用最為廣泛，但設(shè)備品牌眾多，測試精度性能不一，為了解不同F(xiàn)DR 傳感器的性能，選取市面上最常見的3種FDR 傳感器，對田間觀測試驗(yàn)觀測結(jié)果進(jìn)行分析，為傳感器監(jiān)測土壤水分技術(shù)應(yīng)用提供理論基礎(chǔ)。

1 材料方法

試驗(yàn)在北京市昌平區(qū)國家精準(zhǔn)農(nóng)業(yè)研究示范基地（116°46'E，40°18'N）進(jìn)行，試驗(yàn)所用的土壤水分傳感器是基于FDR 技術(shù)的3種國產(chǎn)傳感器，傳感器參數(shù)情況見表1。為保證安裝灌漿后，土壤水分充分運(yùn)移沉降，傳感器提前于2018年10月安裝，正式試驗(yàn)從2019年2月1日開始，10月20日完成數(shù)據(jù)收集，考慮覆蓋通常的3種土壤水分傳感器統(tǒng)一安裝在同一考核場地內(nèi)，相鄰傳感器的間隔距離為1 m，各傳感器每隔1 h 采集1次數(shù)據(jù)，并自動(dòng)上傳到全國土壤墑情監(jiān)測系統(tǒng)網(wǎng)站（www.soilmoisture.org.cn）。資料均以天為單位進(jìn)行統(tǒng)計(jì)整理，共收集了261 d 的土壤水分觀測數(shù)據(jù)、傳感器運(yùn)行總時(shí)長和傳感器正常運(yùn)行時(shí)長。

表1 傳感器基本參數(shù)Tab.1 Basic parameters of the sensor

另外，人工烘干法采集數(shù)據(jù)也在田間同期進(jìn)行，在傳感器監(jiān)測點(diǎn)周圍使用土鉆取不同深度（0～20、20～40、40～60 和60～80 cm）的土壤，每隔1 周取土1次。取土后利用烘干法測定不同土壤深度的土壤含水量，如遇有降雨量超過5 mm 以上情況，降雨停止后2 h 內(nèi)加測一次（如夜間降雨停止，則次日上午加測）。

1.1 傳感器差異性評價(jià)指標(biāo)

計(jì)算3種傳感器的各月的土壤含水量平均值及變異系數(shù)，通過比較不同月份的土壤含水量平均值以及變異系數(shù)的差值大小，來評價(jià)傳感器之間的差異性。其中變異系數(shù)可由下式計(jì)算得出式中：C.V為變異系數(shù)；SD為標(biāo)準(zhǔn)偏差；MN為平均值。

1.2 傳感器數(shù)據(jù)傳輸穩(wěn)定性評價(jià)指標(biāo)

正常率是評價(jià)傳感器數(shù)據(jù)傳輸穩(wěn)定性的一個(gè)指標(biāo)，由傳感器的運(yùn)行總時(shí)長和正常傳輸時(shí)長計(jì)算得出，計(jì)算方法如下式。

式中：T為總時(shí)長；N為正常傳輸總時(shí)長。

1.3 傳感器的二次標(biāo)定過程

人工烘干法一直被認(rèn)為是最精確的測量土壤含水量的方法，因此將傳感器所測數(shù)據(jù)與人工烘干法所測數(shù)據(jù)用多元回歸進(jìn)行線性擬合，得到3種傳感器的標(biāo)定方程。

2 結(jié)果分析

2.1 FDR土壤水分傳感器田間應(yīng)用評估

對3種土壤水分傳感器田間試驗(yàn)數(shù)據(jù)整理后，繪制0～20、20～40、40～60 和60～80 cm 4個(gè)土層的土壤含水量序列圖，如圖1所示。可以看出，3種土壤水分傳感器的觀測值隨時(shí)間呈波動(dòng)性變化，不同類型土壤水分傳感器變化趨勢一致。在時(shí)間變化方向，3種土壤水分傳感器的觀測值在2?6月曲線波動(dòng)較小，在7?10月曲線波動(dòng)較大，以Q1土壤水分傳感器在土層深度為0～20 cm 的曲線變化較為明顯，Q1 土壤水分傳感器在2?6月出現(xiàn)了5次明顯的峰值，而在7?10月出現(xiàn)了8次峰值。這是由于發(fā)生降雨事件時(shí)，土壤含水量增加到峰值，降雨量在7?8月較為集中，日降雨量分布在20～80 mm 范圍內(nèi)。以7月29日為例，該日降雨量為44.9 mm，傳感器Q1、Q2和Q3在0～20 cm 的土壤含水量增幅分別為81.58%、5.64%、16.41%；在20～40 cm 的土壤含水量增幅分別為14.49%、2.83%、14.81%；在40～60 cm 的土壤含水量增幅分別為31.95%、0.23%、2.89%；在60～80 cm 的土壤含水量增幅分別為41.18%、0.30%、10.41%。在空間深度方向，0～20 cm 土層土壤水分傳感器曲線波動(dòng)劇烈，出現(xiàn)了多次波峰、波谷特征，隨著土層深度的增加，土壤水分傳感器曲線波動(dòng)變緩，尤其傳感器Q1 變化較為明顯，在0～20 cm 出現(xiàn)13次明顯的峰值，在20～40 cm出現(xiàn)峰值的次數(shù)減少為5次。

綜上可以看出，3種土壤水分傳感器在時(shí)間和深度方向變化趨勢一致，對降雨量具有明顯的響應(yīng)過程，但不同傳感器之間的絕對數(shù)值有所差異。

圖1 3種FDR土壤水分傳感器土壤含水量隨時(shí)間變化過程圖Fig.1 Three FDR soil moisture sensors soil moisture change process with time

2.2 FDR土壤水分傳感器的差異性

2?10月3種傳感器各月土壤含水量平均值情況如圖2所示。土層深度為0～20 cm時(shí)，傳感器Q2 的整體數(shù)值最高，其次是傳感器Q3 和Q1。2?10月傳感器Q1、Q2 和Q3 的土壤含水 量 平 均值 在6.38%～28.84%、14.56%～37.92% 和11.45%～28.51%范圍內(nèi)變化，其中3種傳感器最高值出現(xiàn)在8月，Q1傳感器的土壤水分最低值在6月，Q2和Q3最低值出現(xiàn)在2月。土層深度為20～40 cm時(shí)，2?7月仍然是傳感器Q2的數(shù)值偏高，到8?10月傳感器Q1 的監(jiān)測數(shù)值升高，傳感器Q3 的數(shù)值整體上均偏低。2?10月傳感器Q1、Q2 和Q3 的土壤含水量平均值分別在3.28%～40.55%、23.38%～37.54%和18.84%～29.29%范圍內(nèi)變化，其中最高值均出現(xiàn)在8月，最低值均出現(xiàn)在2月。土層深度為40～60 cm時(shí)，3種土壤水分傳感器的差異性有所減小，傳感器Q1、Q2 和Q3 土壤含水量平均值的變化范圍分別為19.88%～42.12%、31.92%～36.71%和33.59%～39.74%。土層深度為60～80 cm時(shí)，3種土壤水分傳感器的差異性又有所增加，除8月份以外3種傳感器的數(shù)值由大到小排序依次為：Q2＞Q1＞Q3。

圖2 3種FDR土壤水分傳感器月均土壤含水量變化Fig.2 Monthly average soil moisture change of three FDR soil moisture sensors

2?10月3種土壤水分傳感器的變異系數(shù)如表2所示。從整體來看，3種土壤水分傳感器在0～20 cm 土壤深度的差異性較大，平均變異系數(shù)較高為0.36，而在其他3個(gè)土層的差異性相對較小，其中40～60 cm 土壤深度的平均變異系數(shù)最小為0.11。從不同月份來看，表層土壤（0～20 cm）傳感器的變異系數(shù)在6?7月明顯增加，增幅16.0%，其他土層變異系數(shù)隨時(shí)間呈微小波動(dòng)。

表2 傳感器差異性對比結(jié)果Tab.2 Results of sensor difference comparison

2.3 土壤水分傳感器的數(shù)據(jù)傳輸穩(wěn)定性

3種土壤水分傳感器的數(shù)據(jù)傳輸穩(wěn)定性分析結(jié)果如表3所示。由表3可知，傳感器Q1、Q2 和Q3 均運(yùn)行261 d 共計(jì)6 264 h，但數(shù)據(jù)正常傳輸時(shí)間有所差別。傳感器Q1與Q2的正常率較高，均在90%以上，傳感器Q3 的正常率較低為81%，造成數(shù)據(jù)異常傳輸?shù)闹饕蛴性O(shè)備缺電、物聯(lián)網(wǎng)欠費(fèi)等。

表3 傳感器數(shù)據(jù)傳輸穩(wěn)定性分析結(jié)果Tab.3 Results of the stability analysis of sensor data transmission

2.4 FDR土壤水分傳感器的原位二次標(biāo)定模型構(gòu)建

將一個(gè)土壤水分傳感器所有土層的數(shù)據(jù)與人工測得的數(shù)據(jù)進(jìn)行擬合得到一個(gè)標(biāo)定方程，具體方法為：以傳感器觀測的土壤含水量值作為x軸，人工烘干法測得的土壤含水量值作為y軸，分別進(jìn)行一次線性擬合，得到傳感器Q1、Q2和Q3的擬合效果如圖3所示，可以看出Q2 的標(biāo)定結(jié)果較其他兩個(gè)傳感器標(biāo)好，其相關(guān)系數(shù)為0.571 7。

圖3 FDR土壤水分傳感器的標(biāo)定方程Fig.3 Calibration equation of FDR soil moisture sensor

將傳感器Q1、Q2 和Q3 觀測數(shù)據(jù)分別帶入二次標(biāo)定方程得到標(biāo)定后土壤水分含水量，圖4為標(biāo)定后的2019年2月1日至10月20日3種土壤水分傳感器的土壤含水量序列和日降雨量序列。對比圖1和圖4可以看出，3種土壤水分傳感器的曲線變化趨勢比標(biāo)定前更加一致，波動(dòng)曲線幾乎重合，尤其是0～20 cm，標(biāo)定前3種土壤水分傳感器的絕對數(shù)值差異性較大，最大差值達(dá)到27.71%，而標(biāo)定后波動(dòng)曲線較為接近，最大差值不超過10.81%。此外，原位二次標(biāo)定校正了土壤水分傳感器的異常數(shù)據(jù)，傳感器Q1 在60～80 cm 土層，8?9月降雨量較大時(shí)，標(biāo)定前觀測值遠(yuǎn)大于其他傳感器觀測值[圖1（d）]，傳感器Q1與其他兩個(gè)傳感器之間的最大差值為44.03%，標(biāo)定后觀測值過高的情況改善，傳感器Q1 與其他兩個(gè)傳感器之間的最大差值減小為9.95%[圖4（d）]。

圖5為標(biāo)定后的傳感器各月土壤含水量平均值情況。對比圖2和圖5，標(biāo)定前3種土壤水分傳感器的土壤含水量平均值差異性較大，標(biāo)定后差異性明顯減小，尤其是在土層0～20 cm，標(biāo)定前傳感器Q1、Q2、Q3 的土壤含水量平均值變化范圍相差較大，分別為6.38%～28.84%、 14.56%～37.92%、11.45%～28.51%，而標(biāo)定后土壤含水量平均值變化范圍（除Q1 在2月的數(shù)值外）非常相近，分別變?yōu)?6.52%～21.25%、15.92%～21.53%、17.24%～20.60%。

3 討 論

通過分析傳感器田間評估，發(fā)現(xiàn)3種土壤水分傳感器的觀測值隨時(shí)間變化的趨勢相似，但絕對數(shù)值差異過大。雖然3種土壤水分傳感器采用同一測量原理—FDR技術(shù)[12]，但由于不同廠家傳感器的設(shè)計(jì)結(jié)構(gòu)以及制造工藝的差別，以及出廠標(biāo)定時(shí)不同廠家對影響傳感器性能的因素如溫度、容重等補(bǔ)償程度不同，標(biāo)定方程形式不同，導(dǎo)致不同傳感器在測試同一土壤含水量時(shí)讀數(shù)有所差異。蘇志誠等對2種FDR 傳感器的性能進(jìn)行分析，也發(fā)現(xiàn)2種傳感器存在一定的差異性，與烘干法比較的均方根偏差分別為2.8 和3.8，相對誤差分別為9.9%和13.8%[13]。

圖4 FDR土壤水分傳感器二次標(biāo)定后土壤水分隨時(shí)間變化過程Fig.4 Change process of soil moisture after the second calibration of FDR soil moisture sensors with time

圖5 3種FDR土壤水分傳感器標(biāo)定后月均土壤含水量變化Fig.5 Changes of monthly average soil moisture content after calibration of three FDR soil moisture sensors

通過對比變異系數(shù)，可以看出3種土壤水分傳感器在土壤表層和6?8月的變異系數(shù)較大。研究表明，土壤溫度在夏季較高，并隨土層深度的增加而減小[14]，表層土壤的溫度和水分變化較為劇烈且土壤變異系數(shù)較大[15，16]。溫度會影響土壤的表觀介電常數(shù)，而FDR 技術(shù)通過測量土壤的表觀介電常數(shù)進(jìn)而得到土壤含水量，當(dāng)溫度變化時(shí)，土壤的表觀介電常數(shù)的測量值也發(fā)生變化，導(dǎo)致FDR 傳感器的測量值也隨之發(fā)生變化[17?19]，具體表現(xiàn)為溫度越高，其測試的土壤含水量輸出值越大[20，21]。另外，由于不同土壤水分傳感器對溫度的敏感程度有所差別，最終造成了在高溫條件下設(shè)備差異性增加。為解決這個(gè)問題，許多研究人員通過采集樣本數(shù)據(jù)進(jìn)行分析處理，建立FDR 土壤水分傳感器溫度補(bǔ)償模型，對FDR 傳感器進(jìn)行溫度補(bǔ)償后，減小溫度對土壤水分傳感器的影響，補(bǔ)償結(jié)果與實(shí)際土壤含水量的差距減小[22，23]。

有學(xué)者研究表明，傳感器受土壤結(jié)構(gòu)變化、蒸發(fā)、土壤類型等因素的影響，觀測值與真實(shí)值存在一定的誤差[24，25]。陳海波等人對GStar?I 型FDR 土壤水分傳感器的標(biāo)定結(jié)果進(jìn)行分析，檢驗(yàn)結(jié)果表明標(biāo)定后的FDR 土壤水分傳感器觀測數(shù)據(jù)與人工觀測數(shù)據(jù)具有較好的一致性[26]。本文利用田間原位標(biāo)定方法對3種FDR 傳感器進(jìn)行標(biāo)定有效減小了FDR 傳感器之間的差異性，標(biāo)定后不同土層的傳感器變異系數(shù)如表4所示，在土層0～20、20～40、40～60 和60～80 cm 分別降低了64%、62%、53%和66%。因此，在FDR 傳感器的應(yīng)用中建議在原位進(jìn)行二次標(biāo)定，來消除誤差值。

表4 3種傳感器標(biāo)定前后變異系數(shù)對比Tab.4 Comparison of coefficient of variation of three sensors before and after calibration

4 結(jié) 論

3種土壤水分傳感器的觀測值的曲線變化規(guī)律一致，在空間上表現(xiàn)為隨著土層深度的增加，土壤含水量曲線波動(dòng)幅度減小，在時(shí)間上由于受降雨因素的影響，土壤含水量曲線在夏季（6?8月）波動(dòng)幅度較大。3種土壤水分傳感器的絕對數(shù)值差異較大，尤其是土壤表層（0～20 cm），變異系數(shù)比較大為0.42，而在20～40、40～60 和60～80 cm 的變異系數(shù)分別為0.29、0.15 和0.28。受土壤溫度的影響，土壤水分傳感器之間的差異性在夏季和土壤表層較大，夏季的變異系數(shù)相對于其他季節(jié)較大，增幅為16%，土層0～20 cm 的變異系數(shù)相對于其他土層較大，增幅為38%以上，建議在對土壤水分傳感器標(biāo)定時(shí)應(yīng)實(shí)時(shí)考慮溫度因素進(jìn)行補(bǔ)充校定。通過數(shù)據(jù)擬合，分別建立了3種土壤水分傳感器土壤水分的原位二次標(biāo)定方程，標(biāo)定后的3種土壤水分傳感器之間的絕對數(shù)值差異性明顯減小，其變異系數(shù)在0～20、20～40、40～60、60～80 cm 分別減小了64%、62%、53%和66%。在農(nóng)田應(yīng)用中，應(yīng)考慮土壤溫度、土壤類型因素對FDR 傳感器原位二次標(biāo)定的影響，以消除不同傳感器工廠標(biāo)定的差異性。