350 likes | 547 Views
國立中興大學 水土保持學系 專題討論 ( 四 ) 授課老師:陳文福 教授 指導老師:林俐玲 教授 報告學生:林明毅 學號: 7098042009 報告日期: 2011 / 05 / 27. 土壤水力性質轉換函數之研究. 研究目的 文獻回顧 研究流程 研究方法 研究結果. 研究目的. 隨著各地氣候環境與不同性質的土壤條件下,土壤水分特性曲線會呈現不同特性,一般利用實驗方式求得水分特性曲線,但其過程 繁瑣費時 。
E N D
國立中興大學 水土保持學系 專題討論 (四) 授課老師:陳文福 教授 指導老師:林俐玲 教授 報告學生:林明毅 學號:7098042009 報告日期:2011/05/27 土壤水力性質轉換函數之研究
研究目的 • 文獻回顧 • 研究流程 • 研究方法 • 研究結果
研究目的 • 隨著各地氣候環境與不同性質的土壤條件下,土壤水分特性曲線會呈現不同特性,一般利用實驗方式求得水分特性曲線,但其過程繁瑣費時。 • 本研究利用近年提出之土壤轉換函數(pedotransfer function),運用簡單之土壤物理性質推求水力函數,省去繁雜的實驗步驟與時間的花費。
文獻回顧 • 土壤水份特性關係 • van Genuchten Model • PTFs模式之發展 • 迴歸分析統整
土壤水分特性關係 • Corey(1986)提出之毛細理論,在土壤排水過程中,壓力略低於大氣壓力時,仍不致水份流出,直至吸力超過某一界限值大孔隙才開始排水,且土壤之間孔隙較大,毛細壓力較小。
土壤水分特性關係 • 基質勢能與土壤濕度之關係,依獲得之方法不同分成排水及吸水曲線兩種。 • 為避免複雜計算忽略遲滯現象,一般實驗多採用排水過程取得之水份特性曲線。
van Genuchten Model • Se:相對飽和度 • θ:土壤水分含量 • θr:土壤殘餘體積含水比 • θs:土壤飽和體積含水比 • h:基質勢能頭 • n、 α:曲線代定函數 m
土壤轉換函數 • 由於土壤中張力與孔隙幾何效應複雜,不易以簡單之公式完全表示,已被提出的土壤水分特性曲線模式多是半經驗半物理模式(physico-empirical model),因此便有學者嘗試從土壤性質,使得土壤參數具有物理上的意義,更能正確且有效的推估土壤水分特性。
土壤轉換函數 • Bouma (1989)提出土壤轉換函數(pedotransfer function,PTFs),以擁有的資料轉換為計算所需的資料,使用可靠且具根據性的土壤性質作為預測曲線之工具。 • Minasny (1999) 指出延伸非線性迴歸在參數估計與預測水分含量較類神經網路優良。
土壤轉換函數 • Cornelis 等(2001)依需求的條件及方式不同分為以下三種: • 使用複線性迴歸分析或類神經網路估計特定壓力頭之土壤水分含量。 • 同樣使用複線性迴歸或類神經網路預測方程式之參數。 • 研究具有物理概念模式之水分保持現象或碎形數學及尺度相似性之使用。
土壤轉換函數 • 近年PTFs的發展專注於不同地質區域或不同土壤類型之土壤水力性質,使其有較好的推估效果。 • 國內也開始拓展PTFs,洪靖惠(2008)預測相異質地之水份特性曲線;林佳燕(2008)以Arya and Paris Model為主體來推估土壤水分特性。
迴歸分析統整 • 觀察各文獻採用之基本物理參數,發現大多數研究均會選用總體密度、顆粒粒徑、有機質等與土壤水力性質相關且可經由簡單實驗即可取得之基本性質。其餘尚有某些特殊性質,如:topsoil、砂粒顆粒大小中間值(M50)、pH值與採樣點高程、鋁鐵離子含量等。
迴歸分析統整 • 選定所需之變數後,因無法使用原始數據直接推導出高準確率之迴歸式,必須將變數作數值轉換,常見的轉換為取自然對數、倒數、平方、相互相乘、相加等方法。各研究根據當地土壤性質與資料作為變數取捨,故推估之迴歸式均為獨一無二的。
文獻收集 研究流程 土壤物理性質分析 1.水分係數 2.總體密度 3.顆粒密度 4.粒徑分布 5.有機質分析 複迴歸分析 求得實驗組形狀參數α、n與 土物基本性質迴歸方程式 驗證組土樣代入 迴歸式確定正確性 現地採樣 室內壓力鍋排水試驗 繪製土壤水分特性曲線 vG Model擬合形狀參數α、n 資料與數據統整 PTFs參數及 分析方法彙整 結論 PTFs參數及 分析方法選定
研究方法 • 室內試驗: • 基本物理性質 • 水分特性曲線實驗 • van Genuchten模式計算參數α、n • PTFs參數彙整 選定參數 • 迴歸分析 得形狀參數α、n與土壤基本性質迴歸式
研究方法 • 基本物理性質 土壤水分係數(MF) 因烘乾後土壤理化性質將改變,故先將重量水分含量測出,作為往後數據轉換之依據。 MF=1+(AD+OD)/OD 土壤總體密度(Bulk Density, ρb) 使用金屬管法,由烘乾土重和金屬管體積及可求出乾總體密度。
研究方法 土壤顆粒密度(Particle Density, ρp) 採用比重瓶量測。 粒徑分布(Particle Size Distribution) 採用比重計法,推求出砂粒、坋粒、粘粒之重量百分比組成。對照質地三角圖即可得知。 有機質含量(Organic Matter Content of Soils) 運用氧化還原滴定法測定。
研究方法 水份特性曲線 描述土壤保水特性需要兩個基本參數: 土壤保持水分含量與水分基質勢能,水分基質勢能通常以水力頭、張力或壓力表示。 透過增加於壓力鍋的空氣壓力,測定不同勢能下不同質地土壤所能保存的水分含量。
研究方法 • 各壓力下體積水份含量公式如下: • θi:壓力時之體積水分含量(%) • Mi:壓力時土壤壓克力環重(cm3/hr) • M:空壓克力環重(gm) • mi:對照組之不織布吸水重(gm) • Ms:烘乾土重(gm) • V:壓克力環體積(cm3) • ρ w:水的密度(gm/cm3) ρw
研究方法 vG Model 水分特性曲線擬合 藉由vG Model,以最小平方法結合實驗結果,取得研究所需之水分特性參數α與n。 PTFs參數彙整 由前人文獻中選取重複性大的土壤參數,如:總體密度D、砂粒含量Sa、黏粒含量C、田間容水量FC及有機質含量OM。
研究方法 迴歸分析 迴歸分析為研究一個應變數與一個或多個因變數之間是否存在某種關係的統計方法。 上述vG Model擬合後可得輸出變數(土壤水力參數),接著將輸入變數(基本土壤性質)個別做數值轉換作為因變數。 經由SPSS程式選取符合標準之因變數後,根據實用性及方便性,選取最佳迴歸式為代表。
研究結果 • 壓力鍋實驗中,土樣平衡的時間依壓力高低而有所不同,提高bar數所需時間相對減少,甚至到達8bar與12bar時,加壓時間縮短至一天即可。 • 而後將測得之基質勢能與土壤水分含量相對關係進行vG Model擬合,求得水份特性曲線形狀參數α、n。
研究結果 • 在土樣預濕過程中常達到過飽和而使數據高估,故飽和水分含量利用理論由總體密度及顆粒密度推求;殘餘水分含量則採用風乾土水份體積含量求得。 • 計算與擬合後土壤水分特性曲線參數,包含θs、θr、α和n。
利用求得之θs、θr、α和n,代回van Genuchten公式得知推估水分含量與土壤水分張力關係。
研究結果 • 根據影響土壤保水能力的可能性與參考前人文獻,選定五個土物性質作為回歸分析中的因變數。 • 選擇變數確定後,可將vG Model之α、n參數以函數式表示: α=function(D,Sa,C,FC,OM) n=function(D,Sa,C,FC,OM)
研究結果 • 為了表示土壤水力參數與基本土壤性質間之基礎關連性,運用SPSS相關性分析之功能來探討。 • 結果顯示,α、lnα與D、FC有較大相關性,n、lnn則與FC關係性較大。另外,水力參數經過自然對數轉換後,相關係數數值普遍提高。
尚缺 • 土壤轉換模式迴歸式 • 轉換模式驗證
簡報結束 謝謝聆聽