均匀降雨时坡地之水流动力分析课件.ppt

1、LOGO國立中興大學水土保持學系碩士班國立中興大學水土保持學系碩士班均勻降雨時坡地之水理分析指導教授：謝平城指導教授：謝平城 教授教授學學 生：楊佳錫生：楊佳錫日日 期：期：20100507201005072簡報大綱前言前人研究材料與方法結果與討論結論3 前言 前人研究 材料與方法 結果與討論 結論與建議 n台灣地區可利用之土地有限，隨著人口逐漸增台灣地區可利用之土地有限，隨著人口逐漸增加往坡地尋求土地資源的利用實屬無可避免的加往坡地尋求土地資源的利用實屬無可避免的現象。因為地理位置的關係現象。因為地理位置的關係，常有豪雨發生，常有豪雨發生，大量降雨易造成地表破壞，或帶動不穩定之土大量降雨易造

2、成地表破壞，或帶動不穩定之土體大規模向下移動，不僅造成水土流失，更可體大規模向下移動，不僅造成水土流失，更可能危害人民生命財產的安全，因此對於坡地開能危害人民生命財產的安全，因此對於坡地開發應更謹慎的面對。如何降低逕流水深、減低發應更謹慎的面對。如何降低逕流水深、減低逕流流速，增加地表入滲以成為必須探討的課逕流流速，增加地表入滲以成為必須探討的課題。題。4 前言 前人研究 材料與方法 結果與討論 結論與建議 n開發往往會伴隨著地表的破壞，而本研究主要開發往往會伴隨著地表的破壞，而本研究主要將重點放在坡面排水上，因此將針對地表逕流將重點放在坡面排水上，因此將針對地表逕流加以深入探討。台灣年雨量豐

3、沛，如果無法有加以深入探討。台灣年雨量豐沛，如果無法有效排除地面積水，將會形成大量地表逕流，造效排除地面積水，將會形成大量地表逕流，造成土壤表面的沖蝕或帶動鬆動之土壤岩塊向下成土壤表面的沖蝕或帶動鬆動之土壤岩塊向下崩落，進而造成嚴重的土石災害。對於山坡地崩落，進而造成嚴重的土石災害。對於山坡地排水問題，必須先找出影響形成漫地流的因子排水問題，必須先找出影響形成漫地流的因子，推測水深的變化情形，針對因子，做出相對，推測水深的變化情形，針對因子，做出相對應的防護措施。應的防護措施。5 前言 前人研究 材料與方法 結果與討論 結論與建議 n一般研究漫地流之流速方面，常以經驗公式、一般研究漫地流之流速

4、方面，常以經驗公式、現地調查資料或實驗來進行回歸分析，對於影現地調查資料或實驗來進行回歸分析，對於影響水流流況的因子，缺乏考慮。因此將以水力響水流流況的因子，缺乏考慮。因此將以水力學的觀點，對土層之水流流動情形進行討論。學的觀點，對土層之水流流動情形進行討論。在明渠水力學中，地表通常被視為不透水層，在明渠水力學中，地表通常被視為不透水層，然而地表下水流亦為必須探討的一環，在此使然而地表下水流亦為必須探討的一環，在此使用孔隙介質流理論，解決土層水流之流況，同用孔隙介質流理論，解決土層水流之流況，同時結合地表流況來進行流速、流況之剖析，以時結合地表流況來進行流速、流況之剖析，以期待能在坡地道路的開

5、發中，提供參考依據。期待能在坡地道路的開發中，提供參考依據。6黏性流體黏性流體 材料與方法 結果與討論 結論與建議 前言 前人研究前人研究 7 材料與方法 結果與討論 結論與建議 前言 前人研究前人研究 Staffman(1971)Staffman(1971)以數值模擬及實驗證實了以數值模擬及實驗證實了Beavers and JosephBeavers and Joseph理論及其對邊界理論及其對邊界之處理之適用性。之處理之適用性。Beavers and Joseph(1967)Beavers and Joseph(1967)對二為普西流對二為普西流(PoiseuillePoiseuille

6、flow)flow)流經飽和、具滲透性的多孔介質流經飽和、具滲透性的多孔介質之問題進行分析。之問題進行分析。8 材料與方法 結果與討論 結論與建議 前言 前人研究前人研究 Biot(1956a,1956b)Biot(1956a,1956b)基於下列兩種假設：基於下列兩種假設：（1 1）固體結構應力應變關係滿足廣義虎克定律）固體結構應力應變關係滿足廣義虎克定律（2 2）介質孔隙中流體滿足）介質孔隙中流體滿足DarcyDarcy定律定律推導出多孔彈性介質勢流推導出多孔彈性介質勢流理論理論。宋長虹（宋長虹（19931993）根據根據BiotBiot之多孔彈性介質是流理論建立了多孔介質是流模式、之多孔

7、彈性介質是流理論建立了多孔介質是流模式、層流模式並對其驗證。層流模式並對其驗證。9黏性流體黏性流體 材料與方法 結果與討論 結論與建議 前言 前人研究前人研究 孫永達（孫永達（19981998）以多孔介質流模式配合以多孔介質流模式配合NeviorNevior-Stoke-Stoke描述明渠中之流體運動，描述明渠中之流體運動，求解水流經滲透性孔隙介質時之流速分佈。求解水流經滲透性孔隙介質時之流速分佈。許煜聖（許煜聖（20032003）將草視為孔隙介質，推求各種流況下水流流經植被地表時之速將草視為孔隙介質，推求各種流況下水流流經植被地表時之速度解。度解。10黏性流體黏性流體 材料與方法 結果與討論

8、 結論與建議 前言 前人研究前人研究 HuppertHuppert（19821982）探討黏性流體流經平板之問題，並因水平方向長度遠大於厚度探討黏性流體流經平板之問題，並因水平方向長度遠大於厚度，故忽略水平方向導數，求得水深與流速之關係，以數值方法，故忽略水平方向導數，求得水深與流速之關係，以數值方法求解液面線。並設計實驗驗證結果。根據求解液面線。並設計實驗驗證結果。根據BiotBiot之多孔彈性介質之多孔彈性介質是流理論建立了多孔介質是流模式、層流模式並對其驗證。是流理論建立了多孔介質是流模式、層流模式並對其驗證。Govindaraju et al.Govindaraju et al.（19

9、881988）探討降雨時漫地流於陡坡上之流動情形，假設初始水深為零且探討降雨時漫地流於陡坡上之流動情形，假設初始水深為零且無入流，分別以運動波模式、擴散波模式、數值解畫出流量與無入流，分別以運動波模式、擴散波模式、數值解畫出流量與時間之關係圖，並與以比較三者之差異。時間之關係圖，並與以比較三者之差異。11黏性流體黏性流體 材料與方法 結果與討論 結論與建議 前言 前人研究前人研究 Liu and Mei(1989)Liu and Mei(1989)研究賓漢流體再斜板上漫延的流速及液面，以淺水波理論求解研究賓漢流體再斜板上漫延的流速及液面，以淺水波理論求解水深與流速及流量之關係水深與流速及流量之

10、關係。Pascal(1991)Pascal(1991)結合結合HuppertHuppert及及Liu and MeiLiu and Mei之研究，對非牛頓性流體之重力流之研究，對非牛頓性流體之重力流流經斜面之現象進行研究，分析不同流體特性及流速對其所求流經斜面之現象進行研究，分析不同流體特性及流速對其所求得液面線之影響。得液面線之影響。12黏性流體黏性流體 材料與方法 結果與討論 結論與建議 前言 前人研究前人研究 Govindaraju and Govindaraju and KavvasKavvas（19901990）以淺水波方程式化簡漫地流流況之解析解，考量降雨情況下，以淺水波方程式化簡

11、漫地流流況之解析解，考量降雨情況下，分別給予擴散波與運動波不同的起始條件，求解漫地流之水面分別給予擴散波與運動波不同的起始條件，求解漫地流之水面線。線。薄膜流技術薄膜流技術(2005)(2005)對無窮寬度傾斜面非定常薄膜流之解析，其因無限寬而視為二對無窮寬度傾斜面非定常薄膜流之解析，其因無限寬而視為二維流場，以那維爾維流場，以那維爾-史托克斯方程式及連續方程式對流況作解析史托克斯方程式及連續方程式對流況作解析求解液面之高度變化的情形。求解液面之高度變化的情形。13 前人研究 結論與建議 前言 材料與方法材料與方法 結果與討論 0h-u1 水水 層層u2 土土 層層 xy0h-H u1 水水

12、層層u2 土土 層層 xyii14材料與方法材料與方法0h-u1 水水 層層u2 土土 層層 xy0h-H u1 水水 層層u2 土土 層層 xyii15材料與方法材料與方法水層0singxpyuw12120cosgypw1l ：重力加速度l ：流體動力黏滯係數l ：水層中方向之流體速度l ：水層中流體的壓力l ：水體密度l ：坡地傾角:方向x:y方向g gg1u1pw16材料與方法材料與方法0singnxpnuknyunw22p222220cosgnypnw土層:方向x:y方向l ：土層中方向之流體速度l ：土層中流體的壓力l ：土層中之孔隙率l ：土層中之比滲透係數2u2u2pnpk17材

13、料與方法材料與方法l自由水面處 l切線(x)方向流體應力連續l法線(y)方向流體應力連續l水層與土層之交接面l切線(x)方向流體速度連續l切線(x)方向流體應力連續l法線(y)方向流體應力連續hy 0p10dydu10y 221unu yuyu2121pp 18l土層無限厚之不透水底床土層無限厚之不透水底床l切線(x)方向流體速度為零l法線(y)方向流體速度為零l土層有限厚之不透水底床土層有限厚之不透水底床l切線(x)方向流體速度為零l法線(y)方向流體速度為零材料與方法材料與方法Hy0u2y0u20v 0v 19材料與方法材料與方法l源頭源頭l沒有入流l無窮遠處無窮遠處l水深會接近正常水深

14、，或者 。nhh 0dxdh0 x x0)t,0(u120材料與方法材料與方法)knkhhyy21)(tanxh(cosgupp2w1 在得知控制方程式、邊界條件後，便可對流速進行聯立求解。得到結果如下：一般土壤道路坡面：水層流速：土層流速：)enkhnk)(tanxh(cosguyknppw2p21材料與方法材料與方法)e1(ke1e1)nkhek(hyy21)tanxh(cosguHknpH2knH2knpHknp2w1pppp 在得知控制方程式、邊界條件後，便可對流速進行聯立求解。得到結果如下：土層有限厚之坡面：水層流速：透水層流速：nke1ee hnkee nk)tanxh(cosgu

15、peykn)yH2(knp)yH(kn)yH(knpw2H2pknpppp22材料與方法材料與方法 將水層的流速對水深做積分，土層的流速對土體的厚度做積分，可得到單位寬流量。得到結果如下：無限厚之坡地：水層單寬流量：土層單寬流量：總流量：)hknkhh31()costan-xh(gdyuqpp23wh011)costan-xh(ghkqwp2)kh2nkhh31()costan-xh(gqpp23w23材料與方法材料與方法 將水層的流速對水深做積分，透水層的流速對鋪面的厚度做積分，可得到單位寬流量。得到結果如下：土層有限厚之坡地：水層單寬流量：透水層單寬流量：總流量：)hknkhh31()co

16、stan-xh(gdyuqpp23wh011Hknkkhkcos)tan-xh(gqppppw2)Hknkkkh2nkhh31()costan-xh(gqppppp23w24連續方程式：總流量：土層無限厚之水深對距離的微分關係式：材料與方法材料與方法xqthi)kh2nkhh31()costan-xh(gqpp23w 將降雨i強度視為一側入流，在流況為穩態，初始位置x=0時，u1和u2為零的情況下，可將總流量式代入促入流之連續方程式化簡整理成水深對距離的微分關係式如下:tancosg)hk2nkhh31(x ixddhwpp2325材料與方法材料與方法 由於無法以解析解直接求得水深，因此以阮奇

17、-庫特法求得水深之解：Runge-Kutta Method:各位置之水深關係可表示如下：)h,x(Zxdhd)aa2a2a(61hh4321i1i)2ah,2lx(Zla1ii2)h,x(Zlaii1)ah,lx(Zla3ii4)2ah,2lx(Zla2ii326結果與討論結果與討論 27結果與討論結果與討論 28結果與討論結果與討論用Pi定理分析無因次因子的形式，將參數合併整理，找出能代表影響水深之無因次因子，依照不同的設計可改變的參數有坡度因子 、降雨強度因子1=土壤之滲透性因子2=。antS antS LiLnkp不同坡度之水面線：29結果與討論結果與討論用Pi定理分析無因次因子的形式，

18、將參數合併整理，找出能代表影響水深之無因次因子，依照不同的設計可改變的參數有坡度因子 、降雨強度因子1=土壤之滲透性因子2=。antS antS LiLnkp不同降雨強度之水面線：30結果與討論結果與討論用Pi定理分析無因次因子的形式，將參數合併整理，找出能代表影響水深之無因次因子，依照不同的設計可改變的參數有坡度因子 、降雨強度因子1=土壤之滲透性因子2=。antS antS LiLnkp基本參數之選用：31結果與討論結果與討論用Pi定理分析無因次因子的形式，將參數合併整理，找出能代表影響水深之無因次因子，依照不同的設計可改變的參數有坡度因子 、降雨強度因子1=土壤之滲透性因子2=。antS

19、 antS LiLnkp不同土壤滲透性之水面線：32結果與討論結果與討論用Pi定理分析無因次因子的形式，將參數合併整理，找出能代表影響水深之無因次因子，依照不同的設計可改變的參數有坡度因子 、降雨強度因子1=土壤之滲透性因子2=。antS antS LiLnkp不同土壤滲透性之水面線：33結果與討論結果與討論用Pi定理分析無因次因子的形式，將參數合併整理，找出能代表影響水深之無因次因子，依照不同的設計可改變的參數有坡度因子 、降雨強度因子1=土壤之滲透性因子2=。antS antS LiLnkp不同土壤滲透性之水面線：34結果與討論結果與討論用Pi定理分析無因次因子的形式，將參數合併整理，找出

20、能代表影響水深之無因次因子，依照不同的設計可改變的參數有坡度因子 、降雨強度因子1=土壤之滲透性因子2=。antS antS LiLnkp地表上之流速分佈：35結果與討論結果與討論用Pi定理分析無因次因子的形式，將參數合併整理，找出能代表影響水深之無因次因子，依照不同的設計可改變的參數有坡度因子 、降雨強度因子1=土壤之滲透性因子2=。antS antS LiLnkp地表下之流速分佈：36結果與討論結果與討論 antS 不同坡度：37結果與討論結果與討論 antS 不同降雨強度：38結果與討論結果與討論 antS 不同種類土壤：39l 不透水鋪面水深相當接近。結論結論 本研究主要在探討以水力學

21、的觀點，考慮土層的滲透性，運用半解析解方式來分析水流於坡面上的流動情形，並且觀察改變參數後，對水流流速和水面線的影響。l 固定降雨和土壤種類時，坡度越陡，水面線的改變量越小。固定土壤種類，改變降雨和坡度時，坡度越陡者，降雨強度改變，水面線彼此差異越小，坡度越緩者，降雨強度改變，水面線彼此差異越大。l 坡度因子越大，在不同滲透性因子的土壤，尤其土壤為粗砂時，其水深變化有較明顯的差異，但土壤為砏土、黏土與不透水鋪面水深相當接近。40l 不透水鋪面水深相當接近。結論結論 l 前人之研究一般考慮地表為不透水，因此地表之流速為零，而本研究除了可以描述地表流速與水深外，同時也考慮地下土層之滲透性質，詮釋透水土壤之水流流速，其值雖小，但接近地表時流速不為零。l 坡度越陡，土層單寬流量所佔的比率越大，當坡度接近水平時，其土層流量遠小水層之流量。同一坡度時，降雨強度因子越大，土層流量比率也越大，降雨強度越小，則土層流量比率越小。當不同種類土壤時，在固定的坡度和降雨強度下，粗砂之透水性較好，因此土層之所佔之流量比率也較大，黏土和不透坡面則因透水性低，水層之流量遠大於土層之流量。LOGO