梁晨璟+李春光+趙文娟

摘要：土壤水分特征曲線是土壤水吸力與含水率之間的關系曲線。此次試驗利用壓力膜儀測定銀川北部鹽漬土的土壤水吸力和含水率。對土壤水分特征曲線的Van-Genuchten模型采用Matlab軟件和Microsoft Excel軟件進行擬合。經(jīng)過比較和分析后發(fā)現(xiàn)，用Microsoft Excel軟件與Matlab軟件擬合土壤水分特征曲線精度近似相同，但是Microsoft Excel軟件對于不懂任何編程語言的人來說操作更為簡便。

關鍵詞：土壤水分特征曲線；Van-Genuchten模型；擬合方法；比較

中圖分類號：S152.7+1 文獻標識碼：A 文章編號：0439-8114（2014）01-0056-03

Comparative Analysis of Two Methods Fitting Soil Moisture Characteristic Curve

LIANG Chen-jing1，LI Chun-guang2，ZHAO Wen-juan1

（1. Civil Engineering and Water Conservancy， Ningxia University， Yinchuan 750021，China； 2. Institute of Numerical Computation and Engineering Applications， Beifang University of Nationalities，Yinchuan 750021，China）

Abstract： Soil moisture characteristic curve describes the relationship between water suction and moisture content in soil. This paper is aimed to measure the suction and moisture content by using pressure membrane analyzer in the northern area of Yin Chuan saline soil. Van - Genuchten model is used to fit soil moisture characteristic curve use by Matlab and Microsoft Excel. By comparative analysis， it was found that it had approximately the same accuracy of the soil moisture characteristic curve fitted by using Microsoft Excel and Matlab software. But it is more simple and convenient for people not knowing any programming language by using Microsoft Excel software.

Key words： soil moisture characteristic curve； Van-Genuchten model；fitting method； comparison

收稿日期：2013-03-29

基金項目：國家自然科學基金重大研究計劃培養(yǎng)項目（91230111）

作者簡介：梁晨璟（1988-），女，河北邢臺人，在讀碩士研究生，研究方向為旱區(qū)河流泥沙動力學理論與數(shù)值模擬，（電話）15909510629（電子信箱）

liangmeng88-12@163.com；通訊作者，李春光（1964-），男，教授，博士，主要從事計算機數(shù)學和流體力學研究，（電子信箱）

cglizd@hotmail.com。

土壤水分特征曲線是土壤水的基質勢或土壤水吸力隨土壤含水率變化的關系曲線，是研究土壤水分的保持和運動所用到的反映土壤水分基本特征的曲線。土壤水分特征曲線表示了土壤的一個基本特性，有重要的實用價值：①可以利用它進行土壤水吸力和含水率之間的換算；②它可以間接地反映出土壤中空隙大小的分布；③它可以用來分析不同質地土壤的持水性和土壤水分的有效性；④應用數(shù)學物理方法對土壤中的水分運動進行定量分析[1]。因而，對土壤水分特征曲線的研究是土壤物理學工作者關注的重點問題。

1 材料與方法

此次試驗是針對寧夏銀川北部地區(qū)的鹽漬土進行研究，因此選取銀川市西夏區(qū)雙渠口村、銀川市永寧縣望遠鄉(xiāng)望遠三隊、石嘴山市大武口區(qū)和石嘴山市平羅縣4個地方進行采樣并依次編號為1、2、3、4號。經(jīng)測定它們的土壤物理參數(shù)如表1。

測定土壤水分特征曲線的方法有很多種[2]。此次試驗采用的是直接測定法中的壓力膜儀法。壓力膜儀法與其他方法相比，不僅可完整地測定出整個土壤水分特征曲線，而且可測定原裝土或是擾動土。

1.1 研究選取的模型

由于土壤水分特征曲線受土壤質地、土壤結構、土壤容重、溫度及土壤的膨脹和收縮等影響，所以很難定性擬合土壤水分特征曲線。為了計算和分析方便，常把水分特征曲線近似處理為經(jīng)驗公式和理論模型，如目前比較常用的經(jīng)驗公式有Brooks-Corey模型[3]、Gardner模型[4，5]、Van-Genuchten模型[6]和Gardner-Russo模型[7]。

此次的曲線擬合采用Van-Genuchten模型。Van-Genuchten模型既連續(xù)又有連續(xù)斜率，得到的曲線光滑，對絕大多數(shù)土壤在相當寬的水勢或含水量范圍內具有普遍適用性，并可得到相對導水率的解析解，因而在土壤物理領域得到了最為廣泛的應用[8]，其公式表示為：

■=■■ （1）

式中，θ為體積含水率（cm3/cm3），θr為殘留含水率（cm3/cm3），θs為飽和含水率（cm3/cm3），h為土壤水吸力（cm），m，n，α為擬合參數(shù)，其中m=1-（1/n）。

1.2 擬合參數(shù)初值的確定

在擬合時，Van-Genuchten模型中含水率θ和土壤水吸力h為試驗實測數(shù)值，而殘留含水率θr，飽和含水率θs，α，m，n在擬合過程中都視為非線性曲線擬合參數(shù)。無論用什么方法都需要對這些擬合參數(shù)設置初值。

參數(shù)θr為殘留含水率，它的物理意義為土壤水分特征曲線導數(shù)為零時的土壤含水率。在實際的應用中，一般認為凋萎點處的含水率為殘留含水率，即在土壤水吸力為1 500 cm時的含水率可視為凋萎含水率。參數(shù)θs為飽和含水率，也就是說在土壤水吸力等于零時的含水率。θr和θs可以通過試驗測得，一般土壤的取值范圍為0<θr<θs<0.5 cm3/cm3。參數(shù)n沒有明顯的物理意義但是決定了土壤水分特征曲線的坡度。當n大時曲線緩，當n小時曲線陡。參數(shù)α一般認為是進氣值。吸力小時土壤無水排出，土壤含水率維持飽和；當吸力超過某一臨界值，土壤孔隙中水分開始外排，這一臨界負壓值稱為進氣值。重質黏性土壤進氣值較大，輕質土或結構良好土壤進氣值較小[9]。

1.3 利用Matlab軟件擬合土壤水分曲線

Matlab軟件是用于算法開發(fā)、數(shù)據(jù)可視化、數(shù)據(jù)分析以及數(shù)值計算的高級技術計算語言。在Matlab中內置許多命令可以直接調用來進行數(shù)據(jù)處理，其中命令lsqcurvefit可用于非線性曲線擬合。命令lsqcurvefit的原理是利用最小二乘法，即∑■■[F（x，xdatai）-ydatai]2最小，式中，xdata和ydata為向量，F(xiàn)（x，xdata）為向量函數(shù)值，x為擬合參數(shù)向量。根據(jù)實測得到的數(shù)據(jù)xdata和ydata，找到與函數(shù)F（x，xdata）最佳擬合參數(shù)x[10]。

想要使用命令lsqcurvefit，事先應定義M-函數(shù)文件，而后在腳本函數(shù)中調用命令lsqcurvefit。命令lsqcurvefit用法為a=lsqcurvefit[@fun，a0，xdata，ydata]，式中，a為非線性曲線擬合參數(shù)值，fun為定義的文件名稱，a0為非線性曲線擬合參數(shù)初值，xdata，ydata為試驗實測數(shù)值。

1.4 利用Microsoft Excel 2007擬合土壤水分特征曲線

在Microsoft Excel 2007中內置有“規(guī)劃求解”這一功能，可利用這一方法進行土壤水分特征曲線的擬合。Microsoft Excel 2007的“規(guī)劃求解”為非線性最優(yōu)化代碼。Microsoft Excel 2007 Solver程序代碼是以宏的方式提供調用的。使用時不需要關心其具體的實現(xiàn)方法，只需要和它的對話框進行交互即可[11]。

使用Microsoft Excel 2007具體方法如下：

1）在Microsoft Excel 2007中的一列輸入測得的土壤水吸力值，在另一列輸入含水率，同時把殘留含水率θr，飽和含水率θs，α，m，n在程序中都視為參數(shù)，單寫一行并假定初值。

2）利用Microsoft Excel 2007的計算功能編寫Van-Genuchten模型關系式，根據(jù)現(xiàn)有的數(shù)據(jù)計算出含水率，并求得殘差平方和。

3）利用Microsoft Excel 2007的“規(guī)劃求解”功能計算模型中的參數(shù)。“設置目標單元格”為所算出的殘差平方和并使它等于最小值。“可變單元格”選擇殘留含水率θr，飽和含水率θs，α，m，n所設的初值，“約束”為模型中的擬合參數(shù)并使其設定為非負。

4）在右側的“選項（Q）”可設置最長運算時間、迭代次數(shù)、精度等有關具體運算的相關信息，待選擇好之后確定。選擇“求解（s）”，便可求出擬合曲線的擬合參數(shù)。

5）根據(jù)上述計算結果，利用Microsoft Excel 2007的畫圖功能便可畫出Van-Genuchten模型的土壤水分特征曲線。

2 Matlab 與Microsoft Excel 2007擬合結果比較

圖1至圖4、表2為1～4號土壤的水分特征曲線擬合的結果。

從上述實例可以發(fā)現(xiàn)，利用不同的軟件擬合的結果是不同的。這主要歸結于它們各自算法的不同。用Microsoft Excel 2007軟件進行曲線擬合解決了大多數(shù)人對程序語言不了解的問題。操作人員僅需要了解Microsoft Excel 2007表格公式的計算方法和對“規(guī)劃求解”對話框內容的填充。待所有計算完畢，計算結果又會表現(xiàn)在Microsoft Excel 2007表格內。這樣使得操作人員可以清晰地看到擬合點的具體數(shù)值。但是在計算過程中會發(fā)現(xiàn)，用Microsoft Excel 2007擬合土壤水分特征曲線時，有時會將殘留含水率θr計算為零，這顯然是不對的。出現(xiàn)這種情況的原因是在設置約束條件時，發(fā)現(xiàn)約束條件只能設置5種情況，即大于等于、小于等于、等于、取整和二進制。但是在Van-Genuchten模型中要求擬合參數(shù)都大于零，這樣Microsoft Excel 2007的約束條件便無法完全滿足計算需要。因而用Microsoft Excel 2007只能適用大多數(shù)情況。除去上述缺點，Microsoft Excel 2007在土壤水的低吸力和高吸力下都能很好地擬合出曲線。

利用Matlab軟件擬合土壤特征曲線時，需要操作人對編程語言了解且能熟練應用。但是經(jīng)上述實例發(fā)現(xiàn)，Matlab軟件只是在土壤水的低吸力狀態(tài)下擬合很好。通過與Microsoft Excel 2007軟件比較各自所算得的殘差平方和和相關系數(shù)，可以發(fā)現(xiàn)用Microsoft Excel 2007軟件擬合曲線與Matlab軟件擬合曲線精度相近。

3 結論

無論用哪種軟件擬合曲線，都要注意初值的設定。Matlab與Microsoft Excel 2007相比，兩種軟件各有利弊?？偟膩碚f雖然用Microsoft Excel 2007軟件擬合曲線的精度近似于用Matlab軟件擬合曲線的精度，且Microsoft Excel 2007的操作不需要有任何語言基礎，但是Microsoft Excel 2007軟件在擬合上操作較繁瑣，在設置約束條件上有局限性，有可能使計算結果為零。而Matlab對懂編程語言的人來說更為簡單，僅需調用非線性曲線擬合命令lsqcurvefit便可精確擬合出曲線。

參考文獻：

[1] 雷志棟，楊師秀，謝森傳. 土壤水動力學[M].北京：清華大學出版社，1988.

[2] 徐紹輝，劉建立.土壤水力性質確定方法研究進展[J].水科學進展，2003，14（4）：494-501.

[3] MILLY P C D. Estimation of the Brooks-Corey parameters from water retention data[J].Water Resour，1987，23：1085-1089.

[4] GARDNER W R， HILLEL D， BENYAMIN Y. Post irrigation of soil water I. Redistribution[J]. Water Resour，1970，6：851-861.

[5] GARDNER W R，HILLEL D， BENYAMIN Y. Post irrigation of soil water II. Simultaneous redistribution and evaporation[J]. Water Resour，1970，6：1148-1153.

[6] VAN GENUCHTEN M TH. A closed-from equation for predicting the hydraulic conductivity of unsaturated soils[J]. Soil Sci Soc Am，1980，44：892-898.

[7] RUSSO D. Determining soil hydraulic properties by parameter estimation： on the selection of a model for the hydraulic properties[J].Water Resour，1988，24：453-459.

[8] 朱蔚利，肖自幸，牛健植，等. 兩種模型對土壤水分特征曲線擬合的比較分析[J].湖南農業(yè)科學，2011（17）：47-51.

[9] 李法虎，傅建平，孫雪峰. 土壤水分運動參數(shù)的確定及其靈敏性能分析[J].灌溉排水，1993，12（2）：6-14.

[10] 彭建平，邵愛軍.基于Matlab方法確定VG模型參數(shù)[J].水文地質工程地質，2006，33（6）：25-28.

[11] 余 亮. 利用Microsoft Excel軟件進行非線性擬合的非編程方法[J].微型機與應用，2000（5）：16-17.

（責任編輯 鄭 威）

3 結論

無論用哪種軟件擬合曲線，都要注意初值的設定。Matlab與Microsoft Excel 2007相比，兩種軟件各有利弊?？偟膩碚f雖然用Microsoft Excel 2007軟件擬合曲線的精度近似于用Matlab軟件擬合曲線的精度，且Microsoft Excel 2007的操作不需要有任何語言基礎，但是Microsoft Excel 2007軟件在擬合上操作較繁瑣，在設置約束條件上有局限性，有可能使計算結果為零。而Matlab對懂編程語言的人來說更為簡單，僅需調用非線性曲線擬合命令lsqcurvefit便可精確擬合出曲線。

參考文獻：

[1] 雷志棟，楊師秀，謝森傳. 土壤水動力學[M].北京：清華大學出版社，1988.

[2] 徐紹輝，劉建立.土壤水力性質確定方法研究進展[J].水科學進展，2003，14（4）：494-501.

[3] MILLY P C D. Estimation of the Brooks-Corey parameters from water retention data[J].Water Resour，1987，23：1085-1089.

[4] GARDNER W R， HILLEL D， BENYAMIN Y. Post irrigation of soil water I. Redistribution[J]. Water Resour，1970，6：851-861.

[5] GARDNER W R，HILLEL D， BENYAMIN Y. Post irrigation of soil water II. Simultaneous redistribution and evaporation[J]. Water Resour，1970，6：1148-1153.

[6] VAN GENUCHTEN M TH. A closed-from equation for predicting the hydraulic conductivity of unsaturated soils[J]. Soil Sci Soc Am，1980，44：892-898.

[7] RUSSO D. Determining soil hydraulic properties by parameter estimation： on the selection of a model for the hydraulic properties[J].Water Resour，1988，24：453-459.

[8] 朱蔚利，肖自幸，牛健植，等. 兩種模型對土壤水分特征曲線擬合的比較分析[J].湖南農業(yè)科學，2011（17）：47-51.

[9] 李法虎，傅建平，孫雪峰. 土壤水分運動參數(shù)的確定及其靈敏性能分析[J].灌溉排水，1993，12（2）：6-14.

[10] 彭建平，邵愛軍.基于Matlab方法確定VG模型參數(shù)[J].水文地質工程地質，2006，33（6）：25-28.

[11] 余 亮. 利用Microsoft Excel軟件進行非線性擬合的非編程方法[J].微型機與應用，2000（5）：16-17.

（責任編輯 鄭 威）

3 結論

無論用哪種軟件擬合曲線，都要注意初值的設定。Matlab與Microsoft Excel 2007相比，兩種軟件各有利弊?？偟膩碚f雖然用Microsoft Excel 2007軟件擬合曲線的精度近似于用Matlab軟件擬合曲線的精度，且Microsoft Excel 2007的操作不需要有任何語言基礎，但是Microsoft Excel 2007軟件在擬合上操作較繁瑣，在設置約束條件上有局限性，有可能使計算結果為零。而Matlab對懂編程語言的人來說更為簡單，僅需調用非線性曲線擬合命令lsqcurvefit便可精確擬合出曲線。

參考文獻：

[1] 雷志棟，楊師秀，謝森傳. 土壤水動力學[M].北京：清華大學出版社，1988.

[2] 徐紹輝，劉建立.土壤水力性質確定方法研究進展[J].水科學進展，2003，14（4）：494-501.

[3] MILLY P C D. Estimation of the Brooks-Corey parameters from water retention data[J].Water Resour，1987，23：1085-1089.

[4] GARDNER W R， HILLEL D， BENYAMIN Y. Post irrigation of soil water I. Redistribution[J]. Water Resour，1970，6：851-861.

[5] GARDNER W R，HILLEL D， BENYAMIN Y. Post irrigation of soil water II. Simultaneous redistribution and evaporation[J]. Water Resour，1970，6：1148-1153.

[6] VAN GENUCHTEN M TH. A closed-from equation for predicting the hydraulic conductivity of unsaturated soils[J]. Soil Sci Soc Am，1980，44：892-898.

[7] RUSSO D. Determining soil hydraulic properties by parameter estimation： on the selection of a model for the hydraulic properties[J].Water Resour，1988，24：453-459.

[8] 朱蔚利，肖自幸，牛健植，等. 兩種模型對土壤水分特征曲線擬合的比較分析[J].湖南農業(yè)科學，2011（17）：47-51.

[9] 李法虎，傅建平，孫雪峰. 土壤水分運動參數(shù)的確定及其靈敏性能分析[J].灌溉排水，1993，12（2）：6-14.

[10] 彭建平，邵愛軍.基于Matlab方法確定VG模型參數(shù)[J].水文地質工程地質，2006，33（6）：25-28.

[11] 余 亮. 利用Microsoft Excel軟件進行非線性擬合的非編程方法[J].微型機與應用，2000（5）：16-17.

（責任編輯 鄭 威）