How to calculate Bed Load

Slides:



Advertisements
Similar presentations
楊學成 老師 Chapter 1 First-order Differential Equation.
Advertisements

工職數學 第四冊 第一章 導 數 1 - 1 函數的極限與連續 1 - 2 導數及其基本性質 1 - 3 微分公式 1 - 4 高階導函數.
不定積分 不定積分的概念 不定積分的定義 16 不定積分的概念 16.1 不定積分的概念 以下是一些常用的積分公式。
1 專利師考照班 師資介紹 高苑科技大學 研究發展處 專線服務電話 陳小姐 校內分機 1701.
單元九:單因子變異數分析.
地理信息系统的空间特性 空间实体及其描述 空间问题论述 空间处理方法 北京大学遥感与GIS研究所 程承旗.
Excel - 九十七年度教職員工資訊教育訓練 董建弘.
第六章 专题地图 本 章 要 点 1、掌握专题地图、地图集、电子地图集的定义、分类及其基本特征。
二岸工作經驗分享座談 蔡翔宇.
第六期 IPMA D 級專案副理 種子講師口試 心得報告(題號:46)
不動產投資分析-教學大綱 郭燿禎 博士.
現階段治水對策與新思維 王瑞德 副署長 1 102年12月13日 1.
教學優良教師分享 資訊管理系 陳昌助.
聯合大學 環境與安全衛生工程學系 教 師 簡 介.
實驗計畫資料分析作業解答 何正斌 國立屏東科技大學工業管理系.
第 二 章 水文循環 本章重點 瞭解水文循環的定義 瞭解水文環境因子 瞭解地面水與地下水的特性.
程式設計概論 1.1 程式設計概論 程式語言的演進 物件導向程式 程式開發流程 1.2 C++開發工具
Supplement Data Mining 工具介紹 楊立偉教授 台灣大學工管系 2014 Fall 1.
Java簡介.
Chapter 1 Introduction.
固體力學案例分析 2018/11/23.
(Circular Linked Lists)
机器人学基础 第四章 机器人动力学 Fundamentals of Robotics Ch.4 Manipulator Dynamics
Visual Basic 物件導向程式設計簡介.
OpenID與WordPress使用說明
文學院、管理學院、海科院 暨各系所 及通識教育中心
程式設計專題.
授課教授:林俐玲 鄭皆達 陳鴻烈 林德貴 報告學生:安軒霈
指導老師:黃貞芬 老師 專題組員:B 黃育宇 B 魏志軒 B 平震宇
第十八單元 平面上之參數方程式.
DHCP for W2K.
網路安全技術 OSI七層 學生:A 郭瀝婷 指導教授:梁明章.
虎克定律與簡諧運動 教師:鄒春旺 日期:2007/10/8
指導老師: 王思文 行銷二A 第二組 黃新強 黃秀菁 楊麗馨
A high payload data hiding scheme based on modified AMBTC technique
資料結構 優點 缺點 1 陣列 (Array) 沒有額外變量 (例head, next,...) 運作/操作較簡單 更新資料時,若要保持順序,需要移動較大量資料 靜態結構Static (宣告時已決定了陣列元素多少,不能在程式執行期間增減元素) 2 隊列Queue (FIFO) 容易更新 加入enqueue:
----直線運動 應用力學by志伯 ----直線運動
Ch2多項式函數 2-2 多項式的運算與應用 影音錄製:陳清海老師 資料提供:龍騰文化事業股份有限公司.
PowerPoint 2019/4/9.
授課老師:楊維邦教授 組長:劉秋良 成員:李政均、郭瀚文、鄒震耀
Working Model 2D 朝陽科技大學 工業設計系 邱相文.
網頁程式概論 建國科技大學資管系 饒瑞佶 2015/9 V1 2016/4 V2 2016/9 V3.
2011清大電資院學士班 「頂尖企業暑期實習」 經驗分享心得報告 實習企業:工研院 實習學生:電資院學士班 楊博旭.
第一次Labview就上手 參考書籍: LabVIEW for Everyone (Jeffrey Travis/Jim Kring)
期末考.
緩衝區溢位攻擊 學生:A 羅以豪 教授:梁明章
資訊網路專題 Special Topics on Information Networks
位移與向量(Displacement and Vector)
產品設計與流程選擇-服務業 等候線補充資料 20 Oct 2005 作業管理 第六章(等候線補充資料)
電腦概論考題分析 佛學資訊組 碩一 張榮顯.
流程控制:Switch-Case 94學年度第一學期‧資訊教育 東海大學物理系.
AI-ECON Research Center Department of Economics National Chengchi University.
SQL語法教學 2015/10/15 John.
原子/分子系統架構 Quantum Chemistry Dynamics Monte Carlo.
課程時間:星期二下午2:20-5:20 -> 1:20-4:10 ? 授課教師 逄愛君, 辦公室: 資訊系館 417室 先修課程
道路交通標誌、標線、號誌 與駕駛行為的互動關係 屏東縣私立96駕訓班製作.
例題 1. 多項式的排列 1-2 多項式及其加減法 將多項式 按下列方式排列: (1) 降冪排列:______________________ (2) 升冪排列:______________________ 排列 降冪:次數由高至低 升冪;次數由低至高.
國立台灣大學 關懷弱勢族群電腦課程 By 資訊工程 黃振修
全組專題討論紀要 題目: 講者: 時地: 題 地形性逕流模式之發展與應用 李光敦 教授•國立臺灣海洋大學河海工程學系
非負矩陣分解法介紹 報告者:李建德.
決策支援系統 實例簡介.
補充 數值方法 數值方法.
廖志宏 助理教授 學歷 經歷 公部門與民間產學研究計劃案 證照 De La Salle University, Philippines
多站台網路預約系統之 AJAX即時資料更新機制
一 可靠度問題.
第四組 停車場搜尋系統 第四組 溫允中 陳欣暉 蕭積遠 李雅俐.
10303: How Many Trees? ★★☆☆☆ 題組:Contest Archive with Online Judge
11621 : Small Factors ★★☆☆☆ 題組:Problem Set Archive with Online Judge
主 講 人: 任 立 中 博士 國立臺灣大學管理學院國際企業學系暨研究所副教授
單一河道濃度傳輸驗證 Conversion+ Longitudinal Dispersion.
Presentation transcript:

How to calculate Bed Load 自我介紹 NETSTARS背景與理論 Bed Load或土砂產量研究狀況 模式近五年來應用及發展概況 動畫Show 報告人:謝慧民 致遠管理學院資訊管理學系助理教授

自我介紹─學歷 1992-1996 國立台灣大學土木工程學研究所工學博士 1988-1992 國立台灣大學農業工程學研究所農業碩士 1992-1996 國立台灣大學土木工程學研究所工學博士 1988-1992 國立台灣大學農業工程學研究所農業碩士 1984-1988 國立台灣大學農業工程學系農學士

自我介紹─經歷 服務機關 服務部門/系所 職稱 起訖年月 現職:致遠管理學院 資訊管理學系 助理教授 2002 / 2  2002 / 2   經歷:財團法人水利研究發展中心 副研究員 1999/10 至 2002/3 國立成功大學防災研究中心 工程組 組長 1997/7 至 1999/9 台灣省水利局第十工程處 工務課 副工程師 1996/7 至 1997/7 台灣省水利局第二工程處 工程員 1992/9 至 1993/5

自我介紹─國家考試 (1).81年公務人員特種考試技術人員考試乙等考試,水利工程科中等及格 (2).81年專門職業及技術人員高等考試及格,水利工程技師 (3).84年全國性公務人員高等考試一級考試農業工程科中等及格(榜首) (4).85年專門職業及技術人員高等考試及格,水土保持技師

自我介紹─專長(一) 河川水理輸砂數值模擬(NETSTARS模式)、都市及河川沿岸淹水模擬(二維綜合淹水模式) 灌溉排水工程(含水車理論)、水利工程、水資源工程、水土保持工程、天然災害防災(含土石流)、水質分析、水質模擬、築壩相關工程與管理、水文水理模式應用等 結構化及物件導向程式設計( FORTRAN, Visual Basic ,Visual C++, Delphi, Java) 地理資訊系統GIS[ Arcview、Arc/Info、Mapinfo、Mapguide、IMS等]

自我介紹─專長(二) 防災資訊管理系統規劃設計(曾參與「水利設施防救災管理系統」、「山坡地防災地理資訊系統」及「流域防洪預警系統」等規劃設計工作) 網頁規劃製作管理[Frontpage、Dreamweaver、Flash、動態網頁Dhtml設計、ASP、JSP、XML等] 電子商務應用研究(B-to-B,B-to-C) 資料庫伺服器架設及資料庫規劃設計[SQL Server、Oracle、Dbase 及Access等]

NETSTARS背景與理論1 本模式為一個擬似二維沖淤模式,凡主支流、複雜河系、陡坡、緩坡、水躍、定量流及變量流之水理及相對應底床沖淤特性均可模擬,同時藉由流管之重新分配進行輸砂演算,可以反應河床橫斷面之不規則變化。 一般河川輸砂模式大多以總輸砂量公式去計算輸砂量,所以在不平衡輸砂流況下會有較大的誤差,本模式採用推移載(Bed Load)及懸浮載(Suspended Load)分開計算方式可以反映懸浮質控制之不平衡輸砂狀況。同時因具有處理節點之水理及輸砂分配能力,所以適用於一般網路型河川及水工模型試驗沖淤問題的模擬。

NETSTARS背景與理論2 本模式水理模擬分定量流、迴水演算與變量流三種演算法,定量流模式(Steady-Energy)係根據一維能量方程式配合節點連續關係差分求解,求得水位及流量值,主要取自BRALLUVIAL模式有關水理部份理論;變量流模式(Unsteady-Momentum)係根據de Saint Venant之一維渠道緩變量流連續及動量方程式再配合節點連續關係差分求解,求解水位及流量值,主要採用CHARIMA模式的水理理論及解法;迴水演算方面,係採用GSTARS模式水理部份加入網路節點處理修改而成,具有處理陡坡及水躍的功能,其適用範圍較廣。

NETSTARS背景與理論3 輸砂模擬以水理模擬所得河川網路的水位、流量配合流管理論分成多個流管,假設所分的每個流管在已知斷面輸水能力均相同,並分別計算輸砂量,再代入輸砂連續方程式中求解每個流管中平均底床沖淤情形。因每個時間段流量並非固定,所以每次流管分配的邊界也會有所不同,如此即可模擬橫向的河床沖淤變化。其中輸砂量計算又可分由輸砂公式直接計算總輸砂量及底床載、懸浮載分開計算後再合併為總輸砂量兩種方法,尤其在高懸浮載的河川或不平衡輸砂現象明顯的區域即需使用後者來計算方能符合實際物理現象。

發展NETSTARS模式的目的 建立一個沖積河川模式可預測或估計複雜河系河道底床之沖淤變化 提供多種流況(變量流、定量流;超臨界、亞臨界流、水躍、跌水)的水理計算方法 分流管及分組粒徑計算輸砂量以反應擬似二維之橫斷面底床變化 提供適合河段屬於懸浮載控制的輸砂演算方法 合理解決網路型河川會分流處之流況、輸砂量及懸浮質濃度的分配問題 本研究所發展之模式俾能有助於台灣河川或各種水工模型沖淤試驗規畫設計之參考

文獻回顧(相關模式 ) 一維擬似定量流模式 CHAR-2(Cunge and Perdreau,1973) HEC-6(Tomas and Prashun,HEC,1977,主支流) KUWASER(Simons and Li Assoc,1979,擬似二維沖淤) HEC2SR(Simons and Li Assoc,1980) IALLUVIAL(Karim and Kennedy,1982) BRALLUVIAL(Holly,1985,網路型) GSTARS(Molinas and Yang,1986,流管,河寬變化) SEDICOUP(Holly & Rahuel,1990)

文獻回顧(相關模式 ) 一維變量流模式 CHAR-3(Cunge and Simons,1975) SEDIMENT-4H(Ranjan Ariathurai,1977) CHAR-4(Chollet,1977) FLUVIAL-12(Chang & Hill,1976,河寬變化) CHAR-1(Cunge,1980) CHARIMA(Yang,1986,網路型) UUWSR(Simons and Chen,1979,網路型) 二維模式(變量流) TABS-2(Tomas,1984) MOBED2(Spasojevic,1988)

NETSTARS模式假設及限制 對河川蜿蜓橫斷面上因二次流作用造成之泥砂偏向分佈的情況不適用。 對河川垂直向之二次流無法模擬。 河床床形(bed form)變化也不做預測。 在水理計算中,斷面浸水區均假設為動床部份,所以在每個時距每個流管的浸水區內底床做均勻的沖淤變形,並藉由每個演算時距內重新計算流管邊界來反映不均勻的底床橫向變化。

NETSTARS模式假設及限制 由於總輸砂量公式只適用於輸砂平衡之狀況下,為反映懸浮質運動的機制,本模式將懸浮質運移行為從總輸砂量公式中分離出來,以較合乎懸浮載運動現象之對流擴散方程式來模擬,可適用於不平衡河段之沖淤模擬。 糙度係數的表示可由資料檔輸入設定,也可用河川流量的函數計算代替,或由半經驗式計算求得,在缺乏資料時可先由公式計算其概略值再行檢定或修正。 若要求精度高而調小t 時,所有演算所使用的邊界條件,程式會在歷線資料內自動做線性內插計算。

NETSTARS模式 複雜河系示意圖

The definition sketch of the channel network

NETSTARS計算流程

NETSTARS Computational Process

水理演算 定量流:根據一維能量方程式配合節點連續關係差分求解,求得水位及流量值。 變量流:根據de Saint Venant之一維渠道緩變量流連續及動量方程式再配合節點連續關係差分求解,求解水位及流量值。 迴水演算:也以一維能量方程式求解,適度加入比力方程式求解能量損失的水躍問題,節點以遞迴解決求解平橫水位值。

水理演算---變量流方程式 連續方程式: 動量方程式: A =河道通水橫斷面積=f(y) (); Q =流量(cms) ; t =時間(sec); x =沿水流方向之水平座標 (m); q =單位河段長度之側入流量(cms/m) ; α=動量修正係數; g =重力加速度(); y=水位(m) ; Sf= 摩擦坡度; K=輸水容量(cms); n =漫寧係數 R = 水力半徑

普立斯蒙四點差分法

Preissmann four point finite difference scheme 變量流求解差分式 Preissmann four point finite difference scheme 連續方程式: 動量方程式:

Cunge 及 Liggett 利用泰勒展開式 雙掃法求解

Stream Tube Computation 等輸水能力

Armouring scheme the parameter Nt is related to the duration of simulation time step

輸砂演算---連續方程式 總輸砂量法 連續方程式: 總輸砂量 公式:4個 差分求解:

總輸砂量(Total Load)公式 Yang‘s(1973,1984):則可用於礫石(Gravel),最大粒徑可為10mm; 計算若超出範圍粒徑時,超過150mm及小於0.06mm輸砂量以起動條件控制,而在10mm-150mm之間則有適當參數判斷是否應計算輸砂量。 Ackers & White(1973):粒徑範圍為0.0625mm-2mm,粒徑小於0.0005ft時,所計算之輸砂量誤差也較大,超過2mm也有參數判斷是否應計算輸砂量。 Engelund-Hansen (1967):粒徑範圍為0.0625mm-2mm,粒徑超過150mm及粒徑小於0.06mm的輸砂量均可由起動條件控制。 Van Rijn(1984):適用粒徑範圍為0.0625mm-10mm,粒徑超過150mm及粒徑小於0.06mm的輸砂量均可由起動條件控制。

輸砂演算---連續方程式 分開計算 連續方程式: 推移載 求法:3個 差分求解:

推移載(Bed Load)公式 Meyer-peter & Muller(1948) : Schoklitsch (1935): Van Rijn (1984):

輸砂演算---懸浮載濃度方程式 Suspended Load 對流延散方程式(Convection-Dispersion Equation): Van Rijn (1984) and Holly and Rahuel (1990) Resuspension: Deposition :

Split Operator Approach 懸浮載方程式分離及解法 Split Operator Approach 對流(Convection-step) : Holly-Preissmann two-point four-order scheme Courant Number U/Δx/Δt is less than 1 ( Cr<=1)

Split Operator Approach 懸浮載方程式分離及解法 Split Operator Approach Courant Number U/Δx/Δt is greater than 1( Cr>1)

NETSTARS模式對流計算方法修改 by李森淵博士 懸浮載計算之流移及延散下游邊界條件採控制方程式離散化自動求得。 流移網路演算法當逆流通過節點(變量流漲潮迴水)及Cr>1時,採用順向及逆向掃瞄演算解決此問題,順向掃瞄演算分為1.已知上游邊界掃瞄2.河段順向掃瞄3.節點掃瞄演算等三個步驟重複計算直到全部河段及節點都完成應有的演算為止。逆向掃瞄演算分為1.已知下游邊界掃瞄2.河段逆向掃瞄3.節點掃瞄演算等三個步驟重複計算直到全部河段及節點都完成應有的演算為止。

天然河道中仍有無法計算的案例如下:

解決之道為算不到的部分採分割△t使Cr<1,則問題迎刃而解:

Split Operator Approach 懸浮載方程式分離及解法 Split Operator Approach 縱向延散(Longitudinal dispersion step ) : Tee Scheme finite difference method, using Gaussian elimination method to solve the tri-diagonal matrix

Split Operator Approach 懸浮載方程式分離及解法 Split Operator Approach 橫向延散(Transverse dispersion step ) : Tee Scheme finite difference method, using Gaussian elimination method to solve the tri-diagonal matrix

Split Operator Approach 懸浮載方程式分離及解法 Split Operator Approach 反應(Reaction step ) : Analytical solution

節點(Node)水位流量處理 節點流量和=0,水位相等: Taylor's series expansion link forward sweep node matrix loading link backward sweep

節點(Node)Bed Load分配 依輸砂能力分配進入的砂量

節點(Node)Suspended Load分配 濃度在節點處充分混合

節點(Node)Suspended Load分配 Convection Longitudinal dispersion

Prepare Data 河系斷面位置分佈圖 河系每一橫斷面測量點位置對應高程(X,Y)資料 斷面累距 各斷面糙度係數(如曼寧n值) 河系水力邊界條件(上游流量、下游水位及側流) 各斷面粒徑組成資料 河系輸砂量邊界條件(上游實測懸浮載或率定曲線)

Analysis Data 檢定、驗證(有二次以上之大斷面測量資料) 決定流管數 決定可沖刷層厚度(倍數ALT*最大粒徑) 水理(是否感潮)及輸砂計算方法 輸砂公式選擇 ΔX已知,決定ΔT 決定糙度係數(檢定或用規畫值)