閔俊杰，蔣立新

(1.國電環(huán)境保護(hù)研究院有限公司，南京 210031；2.無錫市城市防洪工程管理處，江蘇 無錫 214000)

1 研究區(qū)域的概況和方法

1.1 研究區(qū)域的概況

本次用于研究的區(qū)域?qū)儆趤啛釒У募竟?jié)氣候，每年的四季比較分明，且氣候?qū)儆跍睾停磕甑钠骄鶜鉁鼐３衷?6℃左右，而每年的降水量則為1203.5mm，一年的降雨基本維持在6月至9月間，此時段的降雨量點(diǎn)全年的60%，年日照數(shù)則為2137.6h，本區(qū)域內(nèi)的光熱能量十分充足，對農(nóng)作物的生長非常有利。本地區(qū)的海拔則牌15-35m間，屬于低丘傾斜的平原地帶，兩個研究區(qū)域平均的坡度大概在12%左右，土地主要用來種茶和竹林，監(jiān)測點(diǎn)的土壤類型則屬于薄層的粗骨土。該土壤的土層比較淺，一般情況下厚度≤1m，并且土質(zhì)十分松散，其中含有30%的礫石[1]。

通過依據(jù)現(xiàn)場的實(shí)際情況以及土地的利用性質(zhì)決定按照8m的間距設(shè)立共計77個監(jiān)測點(diǎn)，在茶園區(qū)域內(nèi)設(shè)立39個，具體編號則從1-40號，并且將30號予以空缺，而在竹林區(qū)域中設(shè)立38，具體編號則從41-79號，并將69號予以空缺。本次監(jiān)測儀器運(yùn)用剖面土壤水分的傳感器，其主要以是TDR時域反射儀技術(shù)，其可以對土壤進(jìn)行測量介電常數(shù)，然后將其加以轉(zhuǎn)變得出土壤中的水份含量。利用剖面土壤水分傳感器對監(jiān)測點(diǎn)進(jìn)行定期的測量，主要測量深度為10cm、30cm深土壤內(nèi)的水分，并且每月一次和覆蓋干濕季節(jié)，對于降雨特別多的雨季更加強(qiáng)監(jiān)測的資料，比如在5月份的10號和13號以及15號開展監(jiān)測，未得到監(jiān)測13號、21號、44號、65號點(diǎn)的30cm位置水分的含量。而在2016年12月到2017年12月間總共測量的14次左右。通過利用在2017年3月到l0月間測量的10次詳細(xì)數(shù)據(jù)融入建立樣點(diǎn)測量土壤水分回歸的模型之中，然后運(yùn)用2016年12月和2017年1月、11月、12月的4次監(jiān)測數(shù)據(jù)對結(jié)果進(jìn)行驗(yàn)證。

1.2 具體研究的方法

研究區(qū)域土壤水分的時間穩(wěn)定性計算公式為：

(1)

(2)

(3)

如果δi=0則表示此監(jiān)測點(diǎn)代表研究區(qū)域內(nèi)的平均濕潤；δi為正值則表示監(jiān)測點(diǎn)的濕潤度要超出研究區(qū)域平均的濕潤度，δi為負(fù)值則與之相反。Sδi越大則表時監(jiān)測樣點(diǎn)的土壤水分時間方面非常的不穩(wěn)定[2]。

通過以時間穩(wěn)定性為重要的研究基礎(chǔ)，并加以利用統(tǒng)計學(xué)軟件對所有監(jiān)測點(diǎn)的土壤中水分的特征進(jìn)行分析。本次使用SPSS22.0統(tǒng)計軟件并采用降維因子進(jìn)行相應(yīng)的分析，主要對茶園區(qū)域與竹林區(qū)域的相關(guān)監(jiān)測點(diǎn)在3月至10月間對土壤進(jìn)行監(jiān)測的詳細(xì)數(shù)據(jù)加以分析研究，利用主成分相關(guān)性的分析方法進(jìn)而分析得出成分圖。在圖中一些相關(guān)性非常強(qiáng)的監(jiān)測樣點(diǎn)聚集在一起，然后與通過時間穩(wěn)定性分析得出典型樣點(diǎn)加以有效的結(jié)合進(jìn)而選擇合理的監(jiān)測點(diǎn)，最終實(shí)現(xiàn)選擇的監(jiān)測典型點(diǎn)可以體現(xiàn)出本次研究區(qū)域中全部樣點(diǎn)。通過有效的利用逐步多元線性回歸(P<0.1)，然后選擇自變量為典型監(jiān)測樣點(diǎn)的土壤水分，而因變量則是其它監(jiān)測樣點(diǎn)，由此最終建立用于研究區(qū)域典型樣點(diǎn)和樣點(diǎn)的土壤水分具有關(guān)系的相應(yīng)模型[3]。利用典型監(jiān)測樣點(diǎn)獲取的監(jiān)測數(shù)值，并經(jīng)過關(guān)系模型詳細(xì)的計算，最得出本次研究區(qū)域內(nèi)所有樣點(diǎn)的土壤水分含量，然后通過實(shí)際的測量數(shù)據(jù)加以校驗(yàn)，最終對預(yù)測結(jié)果的精度利用均方根誤差開展相對應(yīng)的綜合評定。

2 結(jié)果與分析

2.1 時間的穩(wěn)定性分析

通過對2017年3月到10月間對茶園區(qū)域和竹林區(qū)域的所有監(jiān)測點(diǎn)進(jìn)行監(jiān)測得來的詳細(xì)據(jù)開展時間穩(wěn)定性特征的相應(yīng)分析。研究區(qū)域所有樣點(diǎn)平均的土壤含水相對偏差具體如圖1顯示進(jìn)行分布，對所有樣點(diǎn)標(biāo)準(zhǔn)差的標(biāo)注采取垂直誤差線加以顯示。土壤內(nèi)的平均含水相對偏差，茶園區(qū)域基本分布-0.5-0.5之間，如圖1(A)所示，其中有15個監(jiān)測點(diǎn)的時間穩(wěn)定性比較好，其土壤含水相對偏差的標(biāo)準(zhǔn)差是Sδi<0.1；有20個監(jiān)測點(diǎn)的時間穩(wěn)定性一般，土壤含水相對偏差的標(biāo)準(zhǔn)差是Sδi在0.1-0.2；有4個監(jiān)測點(diǎn)的時間穩(wěn)定性較差，Sδi>0.2。而竹林區(qū)域內(nèi)則主要分布在-0.4-0.4之間，具體如圖1B所示，有22個監(jiān)測樣點(diǎn)的土壤含水量具有較好的時間穩(wěn)定性，13個監(jiān)測樣點(diǎn)時間穩(wěn)定性則為一般，還有3個監(jiān)測樣點(diǎn)的時間穩(wěn)定性表現(xiàn)較差。

圖1 茶園區(qū)域(A)和竹林區(qū)域(B)監(jiān)測點(diǎn)平均土壤的水分相對偏差

通過選擇茶園區(qū)域中的監(jiān)測點(diǎn)17號較為典型的樣點(diǎn)作為濕潤的代表，而40號監(jiān)測點(diǎn)的時間穩(wěn)定性特別好，其可作為接近土壤水分平均含量的樣點(diǎn)，而干旱樣點(diǎn)則由9號監(jiān)測點(diǎn)加以代表；時間穩(wěn)定性最差的監(jiān)測樣點(diǎn)則由l0號表示。與上述基本相同，在竹林區(qū)域內(nèi)由79號監(jiān)測點(diǎn)作為濕潤的代表，而接近土壤平均含量最好的則為61號監(jiān)測點(diǎn)，干旱監(jiān)測樣點(diǎn)為53號，時間穩(wěn)定性比較差的樣點(diǎn)為63號監(jiān)測點(diǎn)、通過以時間穩(wěn)定性為基礎(chǔ)進(jìn)行選取典型樣點(diǎn)4個點(diǎn)建立分析的線性回歸模型，其無法有效的對研究區(qū)域內(nèi)的全部監(jiān)測樣點(diǎn)的土壤水分進(jìn)行精確模擬。本次在茶園區(qū)域內(nèi)，在30cm深處共有22個監(jiān)測點(diǎn)取得了較好的模擬精度(如表1所示)，而在10cm的位置，有29個監(jiān)測取得較好的模擬精度，另外還有3個監(jiān)測樣點(diǎn)比較差。在竹林區(qū)域內(nèi)，，在30cm深處共有26個監(jiān)測點(diǎn)取得了較好的模擬精度，另外還有2個監(jiān)測樣點(diǎn)的精度比較差(如表2所示)；在10cm的位置，有35個監(jiān)測取得較好的模擬精度，另外還有2個監(jiān)測樣點(diǎn)比較差。由此表時，當(dāng)分析自變量選擇由時間穩(wěn)定性加以選擇的4個典型樣點(diǎn)時，無法對研究區(qū)域內(nèi)的所有監(jiān)測樣點(diǎn)的土壤水分進(jìn)行精確的模擬[4]。

表1 茶園區(qū)域典型監(jiān)測樣點(diǎn)回歸模型所制定系數(shù)R2的分布數(shù)據(jù)

表2 竹林區(qū)域典型監(jiān)測樣點(diǎn)回歸模型所制定系數(shù)R2的分布數(shù)據(jù)

2.2 降維因子分析

通過對2017年3月至10月間的10次土壤水分的數(shù)據(jù)開展降維因子分析而得出圖2(旋轉(zhuǎn)空間成分圖)。從該圖中可以看出，監(jiān)測樣點(diǎn)的位置和其之間的距離遠(yuǎn)近可以表明監(jiān)測點(diǎn)具有相關(guān)性的高低，經(jīng)過綜合性的考慮和分析所有的監(jiān)測樣點(diǎn)的相關(guān)性和典型樣點(diǎn)，竹林區(qū)域選擇的監(jiān)測樣占為55和45號，而茶園區(qū)域則選擇23和12以及3號。

圖2 茶園區(qū)域(A)和竹林區(qū)域(B)監(jiān)測點(diǎn)旋轉(zhuǎn)空間成分圖

通過在已經(jīng)選擇監(jiān)測典型的樣點(diǎn)為基礎(chǔ)之上加以融入利用降維因子進(jìn)行分析得出的監(jiān)測樣點(diǎn)當(dāng)作建立回歸模型時的自變量，觀察回歸模型的決定系數(shù)R2數(shù)據(jù)。通過對茶園區(qū)域融入添加典型監(jiān)測樣點(diǎn)之后，R2見表1。通過將3號、9號、10號、17號、23號、40號、12號等監(jiān)測樣占當(dāng)作回歸模型的自變量時，在30cm深度位置有29個較好的模擬精度，具體見表1所示，而在10cm深處則有35個較好的模擬精度。并且在總監(jiān)測中，10cm深處位置和30cm深處位置較好的模擬精度分別占總量的89.7%和78.4%，由此表明可以對本次研究區(qū)域內(nèi)的土壤水分加以有效的監(jiān)測。

而在竹林區(qū)域之中融入55號監(jiān)測樣后，見表2，其在30cm深度位置則有32個監(jiān)測樣點(diǎn)具有較好的模擬精度，在10cm深度位置則有36個監(jiān)測樣占具有較好的模擬精度，以及1個較差的模擬精度。在總監(jiān)測中，10cm深處位置和30cm深處位置較好的模擬精度分別占總量的94.7%和88.89%，由此表明竹林區(qū)域中只有5個監(jiān)測樣點(diǎn)便具有較好的模擬精度。

3 結(jié) 論

綜上所述，利用時間穩(wěn)定性與降維因子的結(jié)合而對土壤水分展開相應(yīng)的分析，并選擇科學(xué)合理的典型監(jiān)測樣點(diǎn)的相應(yīng)組合，最終可以有效的利用多元線性回歸模型對監(jiān)測點(diǎn)的土壤水分進(jìn)行準(zhǔn)確的模擬。當(dāng)只運(yùn)用時間穩(wěn)定性作為回歸模型的自變量時，模擬的精度并不能實(shí)現(xiàn)對研究區(qū)域內(nèi)的土壤水分進(jìn)行準(zhǔn)確的模擬，而當(dāng)時間穩(wěn)定性與因子進(jìn)行結(jié)合后，建立模型時所選擇的典型監(jiān)測樣點(diǎn)可以實(shí)現(xiàn)該模型具有非常準(zhǔn)確的模擬精度。