国立中兴大学水土保持学系专题讨论四授课老师陈文

土壤水力性質轉換函數之研究

國立中興大學水土保持學系專題討論(四)授課老師：陳文福教授指導老師：林俐玲教授報告學生：林明毅學號：7098042009報告日期：2011

/

05

/

271研究目的文獻回顧研究流程研究方法研究結果2研究目的隨著各地氣候環境與不同性質的土壤條件下，土壤水分特性曲線會呈現不同特性，一般利用實驗方式求得水分特性曲線，但其過程繁瑣費時。本研究利用近年提出之土壤轉換函數(pedotransferfunction)，運用簡單之土壤物理性質推求水力函數，省去繁雜的實驗步驟與時間的花費。3文獻回顧土壤水份特性關係vanGenuchtenModelPTFs模式之發展迴歸分析統整4土壤水分特性關係Corey(1986)提出之毛細理論，在土壤排水過程中，壓力略低於大氣壓力時，仍不致水份流出，直至吸力超過某一界限值大孔隙才開始排水，且土壤之間孔隙較大，毛細壓力較小。5土壤水分特性關係基質勢能與土壤濕度之關係，依獲得之方法不同分成排水及吸水曲線兩種。為避免複雜計算忽略遲滯現象，一般實驗多採用排水過程取得之水份特性曲線。6vanGenuchtenModelSe：相對飽和度θ：土壤水分含量θr：土壤殘餘體積含水比θs：土壤飽和體積含水比h：基質勢能頭n、α：曲線代定函數m7土壤轉換函數由於土壤中張力與孔隙幾何效應複雜，不易以簡單之公式完全表示，已被提出的土壤水分特性曲線模式多是半經驗半物理模式(physico-empiricalmodel)，因此便有學者嘗試從土壤性質，使得土壤參數具有物理上的意義，更能正確且有效的推估土壤水分特性。8土壤轉換函數Bouma(1989)提出土壤轉換函數(pedotransferfunction,PTFs)，以擁有的資料轉換為計算所需的資料，使用可靠且具根據性的土壤性質作為預測曲線之工具。Minasny(1999)指出延伸非線性迴歸在參數估計與預測水分含量較類神經網路優良。9土壤轉換函數Cornelis等(2001)依需求的條件及方式不同分為以下三種：使用複線性迴歸分析或類神經網路估計特定壓力頭之土壤水分含量。同樣使用複線性迴歸或類神經網路預測方程式之參數。研究具有物理概念模式之水分保持現象或碎形數學及尺度相似性之使用。10土壤轉換函數近年PTFs的發展專注於不同地質區域或不同土壤類型之土壤水力性質，使其有較好的推估效果。國內也開始拓展PTFs，洪靖惠(2008)預測相異質地之水份特性曲線；林佳燕(2008)以AryaandParisModel

為主體來推估土壤水分特性。11迴歸分析統整觀察各文獻採用之基本物理參數，發現大多數研究均會選用總體密度、顆粒粒徑、有機質等與土壤水力性質相關且可經由簡單實驗即可取得之基本性質。其餘尚有某些特殊性質，如：topsoil、砂粒顆粒大小中間值(M50)、pH值與採樣點高程、鋁鐵離子含量等。12迴歸分析統整選定所需之變數後，因無法使用原始數據直接推導出高準確率之迴歸式，必須將變數作數值轉換，常見的轉換為取自然對數、倒數、平方、相互相乘、相加等方法。各研究根據當地土壤性質與資料作為變數取捨，故推估之迴歸式均為獨一無二的。13迴歸分析統整作者名稱年代應變數因變數Teepeetal.2003θs,lnα,ln(n-1),PWPD,C,S,lnSa,lnSa2,S2,Sa/D,S/D,S0.5,Sa0.5,C0.5Kalmanetal.2004θs,lnα,lnnD,OM,C,S,FC,lnC,lnS,Sa2,C2,D2,D*Sa,Sa/S,D2*C2,D*S,D2*C,D2*S2Matulaetal.2007θs,lnα,ln(n-1)D,OM,C,S,lnD,lnOM,lnS,Sa2,S2,OM2,D2,C2,C-1,S-1,OM-1,D-1,OM*C,D*C,D*OM,topsoil,topsoil*C,topsoil*SMajouetal.2007θs,α,nD,OM,C,S,lnD,lnOM,lnS,Sa2,S2,OM2,D2,C2,C-1,S-1,OM-1,D-1,OM*C,D*C,D*OM,D*SSantraetal.2008lnKs,lnα,nD,OM,C,S,pH,elevation14研究流程15研究方法室內試驗：基本物理性質水分特性曲線實驗vanGenuchten模式計算參數α、nPTFs參數彙整選定參數迴歸分析得形狀參數α、n與土壤基本性質迴歸式16研究方法基本物理性質土壤水分係數(MF)因烘乾後土壤理化性質將改變，故先將重量水分含量測出，作為往後數據轉換之依據。MF=1+(AD+OD)/OD土壤總體密度(BulkDensity,ρb)使用金屬管法，由烘乾土重和金屬管體積及可求出乾總體密度。17研究方法土壤顆粒密度(ParticleDensity,ρp)採用比重瓶量測。粒徑分布(ParticleSizeDistribution)採用比重計法，推求出砂粒、坋粒、粘粒之重量百分比組成。對照質地三角圖即可得知。有機質含量(OrganicMatterContentofSoils)運用氧化還原滴定法測定。18研究方法水份特性曲線描述土壤保水特性需要兩個基本參數：土壤保持水分含量與水分基質勢能，水分基質勢能通常以水力頭、張力或壓力表示。透過增加於壓力鍋的空氣壓力，測定不同勢能下不同質地土壤所能保存的水分含量。19研究方法各壓力下體積水份含量公式如下：θi：壓力時之體積水分含量(％)Mi：壓力時土壤壓克力環重(cm3/hr)M：空壓克力環重(gm)mi：對照組之不織布吸水重(gm)Ms：烘乾土重(gm)V：壓克力環體積(cm3)ρw：水的密度(gm/cm3)ρw20研究方法vGModel水分特性曲線擬合藉由vGModel，以最小平方法結合實驗結果，取得研究所需之水分特性參數α與n。PTFs參數彙整由前人文獻中選取重複性大的土壤參數，如：總體密度D、砂粒含量Sa、黏粒含量C、田間容水量FC及有機質含量OM。21研究方法迴歸分析迴歸分析為研究一個應變數與一個或多個因變數之間是否存在某種關係的統計方法。上述vGModel擬合後可得輸出變數(土壤水力參數)，接著將輸入變數(基本土壤性質)個別做數值轉換作為因變數。經由SPSS程式選取符合標準之因變數後，根據實用性及方便性，選取最佳迴歸式為代表。22研究結果壓力鍋實驗中，土樣平衡的時間依壓力高低而有所不同，提高bar數所需時間相對減少，甚至到達8bar與12bar時，加壓時間縮短至一天即可。而後將測得之基質勢能與土壤水分含量相對關係進行vGModel擬合，求得水份特性曲線形狀參數α、n。23研究結果在土樣預濕過程中常達到過飽和而使數據高估，故飽和水分含量利用理論由總體密度及顆粒密度推求；殘餘水分含量則採用風乾土水份體積含量求得。計算與擬合後土壤水分特性曲線參數，包含θs、θr、α和n。24張力水分含量336.6cmH2O1020cmH2O2040cmH2O4080cmH2O8160cmH2O12240cmH2O新竹峨眉A13.28611.1379.6348.5957.3786.519新竹峨眉C24.93721.90518.87317.06415.09613.128苗栗大湖A17.76914.77811.9029.6018.6247.991苗栗大湖B21.99619.70617.33616.38815.24313.979南投中寮B30.88227.33323.94622.33420.45218.731南投中寮C22.49520.27817.30314.97813.35611.625台南南化B34.37330.34626.31824.56922.66122.078台南南化C38.00232.20527.18123.70321.93620.777台東卑南A40.78732.81229.70726.97825.69324.783台東卑南B49.7844.05341.26936.50535.00732.812花蓮壽豐A42.98238.86035.64931.2629.76128.744花蓮壽豐C39.17934.39234.41129.50527.72526.64825土樣名稱θsθrαn新竹峨眉A50.861.90.11.275新竹峨眉C49.3420.11.195苗栗大湖A42.252.20.02071.2711苗栗大湖B43.0030.03451.1415南投中寮B40.041.70.01531.1515南投中寮C41.311.50.06631.1929台南南化B53.781.90.01731.1944台南南化C52.975.50.01691.1793台東卑南A55.045.10.00421.1501台東卑南B54.544.30.05641.1256花蓮壽豐A50.363.40.05721.1555花蓮壽豐C49.192.00.00191.237226利用求得之θs、θr、α和n，代回vanGenuchten公式得知推估水分含量與土壤水分張力關係。27研究結果根據影響土壤保水能力的可能性與參考前人文獻，選定五個土物性質作為回歸分析中的因變數。選擇變數確定後，可將vGModel之α、n參數以函數式表示：α=function(D,Sa,C,FC,OM)n=function(D,Sa,C,FC,OM)28土樣名稱總體密度D(g/cm3)有機質OM(%)田間容水量FC(%)砂粒含量Sa(%)粘粒含量C(%)新竹峨眉A1.2561.318.6674.5768.255新竹峨眉C1.3071.424.9466.659.607苗栗大湖A1.5212.224.9566.0134.447苗栗大湖B1.4623.030.8969.5387.509南投中寮B1.5781.230.8853.83235.497南投中寮C1.5361.822.4960.58232.235台南南化B1.3712.44623.67656.162台南南化C1.3952.444.4822.52656.383台東卑南A1.3332.036.7628.65947.475台東卑南B1.2981.436.3