國立中興大學水土保持學系碩專班 專題討論(四)(9375)

Slides:



Advertisements
Similar presentations
側向延散對雙井循環流場追蹤劑試驗溶質傳輸的影響
Advertisements

地圖符號:圖例 地圖符號:圖例.
集水區水文特性變異分析 指導老師:陳文福 學生:李懿書.
因應氣候變遷中央管河川防洪標準之探討 授課老師:鄭皆達老師 指導教授︰陳文福老師 中興大學 水土保持學系碩專班 專題討論(二)
太陽_介面設計.
行銷研究 單元三 次級資料的蒐集.
探析台灣地區醫師繼續教育提供者之意見與態度
台灣國家公園巡禮 製作:翁顗棋.
課程大綱 衛星通信與導航 96年度第一學期(二技).
邊坡穩定數值模擬以 苗栗縣食水坑地滑地為例 專題報告 報告人:陳昇德 學號: 中華民國九十八年五月十六日.
題目:台中市大坑地滑區 可用地下水量之推估 學生:温偉建 學號: 日期:2010年1月9日
實驗計畫資料分析作業解答 何正斌 國立屏東科技大學工業管理系.
課程大綱 衛星通信與導航 96年度第一學期(四技).
第5章習作 水 文.
A CLOSED-FORM EQUATION FOR PREDICTING THE HYDRAULIC CONDUCTIVITY OF UNSATURATED SOILS 預測不飽和土壤水力傳導度之方程式 M.TH van GENUCHTEN. 指導教授: 林俐玲 教授 授課教授: 謝平城.
固體力學案例分析 2018/11/23.
國立台灣大學生物產業機電工程研究所 簡君良
埤塘之水資源與防災功能定位與應用 財團法人農業工程研究中心 簡傳彬 博士 民國102年6月5日,土木館五樓 埤塘的功能定位 水資源調配
Android App簡介及 App Inventor 2體驗 靜宜大學資管系 楊子青
SuperGIS 2.0 基本架構介紹.
網路安全技術 OSI七層 學生:A 郭瀝婷 指導教授:梁明章.
標籤、按鈕、工具列、狀態列 (Labels, Buttons, Tool Strips, and Status Strips)
建立一 function s (type) 可以用來繪製cyclic-harmonic curves
龍井區步道行 竹坑-南寮步道自然探究.
授課教授:石文傑 博士 報 告 人:施金橐 學 號: 日 期:98年03月06日
台灣的人口 教學者:李蕙娟 年級:二年級.
資源未來 佛光人文社會學院 第07章 台灣降雨長期變為化對水資源之衝擊 楊棋明 博士 中央研究院生物多樣性研究中心
1999 Chi-Chi Earthquake (ML=7.3 MW=7.6)
Definition of Trace Function
有關於股票報酬及匯率變化對台灣醫療產業市場收益的分析
信度分析 (11/7~11/13) 1.何謂『信度』 2.信度分析步驟.
註:應在1280x 色之顯示器,PowerPoint 50%之顯示比例,
趙義隆 教授 台大國際企業學研究所 趙義隆 教授 台大國際企業學研究所 2003 MBA How-to #2.
緩衝區溢位攻擊 學生:A 羅以豪 教授:梁明章
國立中興大學水土保持學系 專題討論 專題討論(二) 坡地農業沉砂滯洪設施及植生緩衝帶 配置區位優選之研究
電子量角器 撰寫人:董瑩蟬.
水與綠建設計畫 滾動式檢討報告 行政院經濟建設委員會 93年7月1日.
行的色彩 交通標誌(4).
國立成功大學(農業) 報告人 協同主持人 林翰佑
8公分 中文論文名稱 (大小80~136) English Thesis Title (size80~136) 中文名字 English Name (大小50~60) 學校學系/任職單位職稱 (大小50~60) 8公分 8公分 內容 (大小40~50) 內容 (大小40~50) 文字或圖示 文字或圖示.
人類的石油年代.
有限差分法在大地工程之應用 – 以壓密沉陷問題為例
PSoC(Programmable SoC)的架構與優勢
具繪畫藝術風格之非真實照片式繪圖 Non-Photorealistic Illustrations with Artistic Style
國際金融期末報告 授課老師:楊奕農副教授 學生:金恩惠 學號:
電腦概論考題分析 佛學資訊組 碩一 張榮顯.
(中文標楷、英文Times New Roman,85點,粗體,置中)
物理化學輔助學習工具 2018/12/04.
國立臺北科技大學機械系106學年度實務專題海報
台北盆地,充滿了活潑生命力的土地,且讓我們一同來探尋它的傳奇...
幼葉林地滑地調查研究 授課老師:林信輝老師 指導老師:林俐玲老師 學生: 鄭裕適 年12月11日.
聖地的位置.
受壓含水層中水平井於河邊抽水之解析解 An Analytical Solution for Horizontal Well Pumping by the Riverside in Confined aquifers 鄒佩蓉.
例題 1. 多項式的排列 1-2 多項式及其加減法 將多項式 按下列方式排列: (1) 降冪排列:______________________ (2) 升冪排列:______________________ 排列 降冪:次數由高至低 升冪;次數由低至高.
10328: Coin Toss ★★★☆☆ 題組:Problem Set Archive with Online Judge
4-1 變數與函數 第4章 一次函數及其圖形.
-創新學習 從教師開始- 主講人 國立臺灣師範大學 李忠謀 教授
CINII操作手冊 圖21.會員可檢索資料之範圍及入境 圖1.CINII基本檢索畫面 圖22.單位別逐一檢索的畫面及操作方式
一 可靠度問題.
第四組 停車場搜尋系統 第四組 溫允中 陳欣暉 蕭積遠 李雅俐.
自來水系統 生科、航太、機械、通訊及研發 自來水 下水道系統 循環水 總用水量 光電及半導體產業 污水處理廠 排放140640
NDL委託代工流程-PECVD 葉佳翰 2019/8/6.
單元三:敘述統計 內容: * 統計量的計算 * 直方圖的繪製.
后里鄉的地理位置.
認識比例尺 廖建榮.
主題研究架構.
2-3 地 殼 變 動.
第6課 資源問題與環境保育對策 一、水資源 (一)整體問題:不足、時空分布不均 (二)分區問題: 1.季風氣候區: (1)問題: a.水汙染
單一河道濃度傳輸驗證 Conversion+ Longitudinal Dispersion.
物理化學輔助學習工具 2018/12/04.
Presentation transcript:

國立中興大學水土保持學系碩專班 專題討論(四)(9375) 題目:以模組化三維地下水流模式(MODFLOW) 分析台中盆地地下水流之研究 授課教授:陳文福 教授 指導教授:謝平城 教授 學生:温偉建 學號:59742207 日期:2010年5月15日

簡報大綱 一、前言:緣起、目的 一、前言:緣起、目的 一、前言:緣起、目的 一、前言:緣起、目的 一、前言:緣起、目的 一、前言:緣起、目的 二、文獻回顧 二、文獻回顧 二、文獻回顧 二、文獻回顧 二、文獻回顧 二、文獻回顧 三、研究材料與方法 三、研究材料與方法 三、研究材料與方法 三、研究材料與方法 三、研究材料與方法 四、研究結果與討論 四、研究結果與討論 五、參考文獻 五、參考文獻

前言:緣起 人類因過度開發逐年超量排放二氧化碳導致全球暖化氣候異常極端,台灣本島乃屬於氣候異常的受害區域之一。近年來極端的降雨情形時有所聞,而超量的降雨形成過度之地表逕流所引發的山崩、地滑、土石流又造成大地及國人的二度傷害;同時,降雨季節的極端化更加劇人民陷入缺水的困境,故地下水的含量與開發就越顯得珍貴及迫切。

前言:目的 本研究便針對台中盆地的地下水層進行數值模擬,為了能評估台中盆地之地下水利用的可行性,必須充分了解其水文地質形成的機制與其邊界條件;本研究便以地質觀點切入,對台中地區水文地質形成機制進行探討,提出地質概念模型,做為建構水文地質模型之依據,並配合現有之水文資料了解區域內的地下水概況建立地下水流之數值模型進行模擬,以期可供後續研究對此區域之地下含水量及安全出水量進行推求而使相關單位可妥善利用地下水之珍貴水資源。

文獻回顧 經濟部水利署水利規劃試驗所(2007): 地形 台中盆地東緣的麓山帶丘陵區是南北走向車籠埔斷層逆衝形成,西緣的大肚台地及八卦台地則是彰化斷層(包含清水斷層)向西逆衝抬升所造成,地勢由南北兩端向中央傾斜,盆地內沖積扇有北部之豐原沖積扇(大甲溪古沖積扇),中部之太平聯合沖積扇,與南部之烏溪沖積扇。周圍地形區則有西北端之大甲溪沖積平原,北邊之后里台地,東側之東勢丘陵地、豐原丘陵地及南投丘陵地,以及西側之八卦台地和大肚台地。

台中盆地周圍地形分區圖

文獻回顧 地質與地層 經濟部水利署水利規劃試驗所(2007): 1.上新世至更新世岩盤(P) 2.頭嵙山層(PQ) 主要由砂岩、砂質泥岩、礫石質砂岩及厚層狀礫岩所組成,下段為香山段(PQs)、上段為火炎山段(PQc)。

3.紅土台地堆積層(Q3) 由礫石、砂、黏土組成,礫石以砂岩及石英砂岩為主,為古河相沉積物。本層以不整合覆於頭嵙山層之上,主要分佈於大肚台地及八卦台地東斜面及西斜面上。 4.山麓斜面堆積層(Q4) 本層分佈於八卦台地及大肚台地兩翼坡度較緩處,大部份為植被所掩覆。本層岩性以礫石、泥為主,是台地物質(主要為頭嵙山層岩層)受侵蝕再堆積而成,推估其年代為晚更新世至全新世。 5.沖積層(Q6) 為近期河流堆積之產物,岩性以礫石、砂及泥,主要分佈在現有河流的河床、台中盆地、大肚台地及八卦台地西側的沖積平原。

台中盆地地質圖(中油公司台中圖幅1/100,000)

文獻回顧 地層鑽探 台灣第區地下水觀測網第三期94年度計畫-台灣糖業公司,中央地質調查所委辦(2005): 1.大肚台地:依蒐集資料分別繪製南北向剖面及東西向剖面,根據剖面資料,背斜軸東翼並無阻水層存在,只有一自由含水層,而西翼並在此研究範圍內故不予討論。 2.台中盆地:大致由北邊台中縣豐原往南至南投名間,地形上以中間烏日-霧峰連線為低點,故暫依地形上之不同,將台中盆地在水文地質上分為南、北兩盆地。

北盆地:根據剖面資料得知水文地質上只有一自由含水層,往南至烏日在深度46~58公尺處才有一較厚之泥層出現成為阻水層,產生一自由含水層及一拘限含水層。 南盆地:因有泥層出現,故拘限含水層便明顯出現於南盆地中。南側營盤口及軍功寮兩井未有自由含水層出現,而泥層以下出現一拘限含水層;北側霧峰坑口井顯示淺層有一自由含水層存在,深層因數層泥層出現,形成拘限含水層,但由於阻水層厚度較薄(不到10公尺),拘限含水層之分層及分布仍待進一步調查;草屯以北至烏溪沿岸一帶地層資料較為匱乏,推測應與霧峰一帶相似(因同屬烏溪沖積扇)。

蘇清林(2007),於該論文中提及: 位於烏溪流域內之經濟部水利署觀測站網(2006)鑽探井分別為大肚站、烏日站與大里站,鑽探深度約為地表下200公尺深,鑽探資料顯示大里站鑽探資料大都為礫石層,烏日站除40公尺與100公尺處有部份黏土與粉土層外,其於大多為砂礫石層。大肚站則因接近烏溪出海口,其岩層受海進海退影響,屬砂與黏土互層現象。故鑽探深度60公尺位置以上之砂礫石層為主要含水層,亦即第一含水層,於該研究數值模式中設定該含水層為自由含水層形式,而鑽探深度60公尺位置以下至100公尺間為第二含水層,於數值模式中設定該含水層為拘限含水層形式。

台中地區觀測井及水文地質剖面位置圖(台糖-2005)

文獻回顧 地面水文、流量 1.經濟部水利署水利規劃試驗所(2007): 台中盆地內持續觀測水位流量之水文站有四處,以烏溪最下游大肚橋水文站為例,歷年統計平均流量約為116cms,平均年逕流量約為36.6億M3,豐水期為5至10月,流量約佔全年的70%,以6月份最多。 2.陳尉平(2006): 蒐集中央氣象局自1990~2003年台中測站之蒸發量數值,其年平均為1083.4公厘。

台中盆地各水文站歷年月流量統季表(水規所,2003)

文獻回顧 斷層概況 徐鐵良(1992)指出: 雙冬斷層位在本地質區之東界,位於火炎山及頭嵙山之東,為一逆斷層,東側為昇側,約呈南北至北30 度東走向,斷層東側出露地層為中新世晚期地層,西側地層為頭嵙山層,地層層位落差3000m 。 車籠埔斷層位在本區西緣鄰台中盆地地帶,自南投民間至台中豐原,約呈南北走向;斷層東側為昇側 。出露地層為上新世之錦水頁岩及卓蘭層,西側則出露更新世頭嵙山層及現代沖積層。

三、研究方法與材料 首先先針對台中盆地之水文地質資料進行蒐集並加以分析,建構該區水文地質概念模式,並研究該區內的目標水文地質架構,依據概念模式與地下水流控制方程式(Governing Equation)建立地下水流之模擬式。因地質作用與沉積作用的影響常使得含水層呈非均質性(Heterogeneity)與異向性(Anisotropy)及邊界條件的複雜性,故對其進行解析解有其困難性,一般仍以數值模擬為主。 再者將模擬所輸出的水位近似解依試誤法(Trial and Error)配合面積控制法求得最符合現地狀況的水文地質架構、地下水流模式及水力參數。

研究方法 數值模式 目前常用於地下水流模式的數值方法有:有限差分法(Finite Difference)、有限元素法(Finite Elements)、 積分有限法(Integrated Finite)、邊界積分方程法(Boundary Integral Equation Method)、解析元素法(Analytic Element)等五種方法,但較常使用的是有限 差分法與有限元素法;近年來,這些數值方法被應用及 撰寫成計算地下水流動及溶質傳輸現象的模擬套裝軟 體,而本研究則是採用PMWIN(Processing MODFLOW For Windows﹐Chiang&Kinzelbanch﹐1996)。 ﹐

研究流程

研究材料 初始地下水位 由經濟部水利署:水文水資源資料管理供應系統內挑選台中盆地數個涵蓋年份較長(2006~2008)之地下水位觀測井做為地下水位資料分析,並以各水位觀測井之2006年3月的地下水位做為初始水位進行分析。

研究材料 研究材料 研究材料 網格設定與邊界條件 網格設定與邊界條件 台中盆地依據其地形及地質特性,將其邊界以西北側大肚山台地,北側大甲溪,西南側八卦台地,東北側豐原丘陵,東南側南投丘陵為界,其邊界條件如下頁圖示,盆地內河川水系由北向南依序為筏子溪、大里溪及烏溪。網格所採元素為東西向35行,南北向94列,間距為500m ,總網格數為3290個,有效網格數為1763個。

台中盆地邊界範圍、邊界條件及網格劃分 大甲溪 特定流量邊界 豐原丘陵 特定流量邊界 大肚台地 特定流量邊界 大肚溪 River 八卦台地 南投丘陵 特定流量邊界 台中盆地邊界範圍、邊界條件及網格劃分

地質分層 本模擬區依蒐集之資料顯示全區可視為存在一自由含水層,此外在烏日地區及烏溪上游霧峰附近存在一拘限含水層,但因現階段對此拘限含水層之文獻及資料極為稀少,故本研究針對自由含水層進行地下水位之模擬;並將地層視為透水層依1974年水資會將台中盆地之滲透係數畫分為五區,如下圖表所示。

台中盆地透水層分區示意圖

研究材料 滲透係數及補注量 地區 K(m/s) 補注量 (1×108/year) 北 盆 地 西側 (大肚山台地) 0.00032 0.628~0.988 北側 (大甲溪) 0.00109 1.75~2.68 東側 (豐原丘陵) 0.00030 2.21~4.1 合計 4.6~7.77 南 西側 (八卦台地) 0.00275 3.9 東側 (南投丘陵) 0.00102 1.64~2.93 10.14~14.6 台中盆地地下水層滲透係數及補注量分區表

研究材料 抽水量 本研究區域範圍之抽水量參考水規所於2009年所做之調查,擷取2006年之各縣市抽水量再平均分配於網格中。 年份 縣市 總計 1×104m3 2006 台中縣 13,200 台中市 4,328 南投縣 5,896 彰化縣 681 24,231 2006年台中盆地地下水抽水量統計表

四、研究結果與討論 數值模式經過參數的敏感度分析及率定後,原靠近各邊界之格網模擬值與觀測值相差甚異,經側向補注量的調整之後,水位的模擬值可達到與觀測值接近的狀況,以下將列舉水位觀測站之模擬與觀測值折線圖說明其模擬成果與問題。

年/月 2008/1 2008/2 2008/3 數值(m) 測站 觀 測 值 模 擬 誤 差 誤差 新庄 潭子 大雅 141.57 141 0.57 140.59 140.5 0.09 140.33 -0.17 潭子 105.19 105 0.19 103.56 103 0.56 102.73 -0.27 大雅 121.59 121 0.59 120.06 120.5 -0.44 119.63 120 -0.37 觀測、模擬水位與誤差值一覽表

誤差可能原因:雨季或颱風強降雨所造成當月地下水位之增加明顯所致。 新庄站觀測與模擬水位折線圖

結論與建議 本研究之數值模式可提供後續評估該區域地下水量抽取量之依據。 透過MODFLOW可模擬地下水位,但後續若能考慮地下水之水質更能確定地下水之使用條件。 本研究因缺乏更深層之地質鑽探資料及相關參數,故僅針對自由含水層進行模擬,若能增加此方面之資料便可在進行拘限含水層之水位模擬。 本研究模擬時初次得到之地下水位於北側靠大甲溪之網格明顯低於觀測值,可能因溪水透過入滲補注至地下水層之水量高於輸入之補注量,故建議後續之研究人員應考慮大甲溪河畔之入滲量。

五、參考文獻 1.林時猷(2003),以地下水位之區位相關性輔助濁水溪 沖積扇地下水模擬之參數決定,逢甲大學土木及水利 工程研究所碩士論文。 2.蔡清研(2007),濁水溪沖積扇整合模式下之MODFLOW 地下水模擬研究,國立中正大學應用地球物理研究所 碩士論文。 3.潘文健(2002),屏東平原合適出水量分析之研究, 國立成功大學資源工程學系碩士論文。 4.黃清譽(2007),彰化地區地面水與地下水聯合供水管理 調配之研究,國立成功大學資源工程學系碩士論文。

5.蘇清林(2007),河畔取水所引致河川滲漏量之評估,國立 成功大學資源工程學系碩士論文。 6.丁澈士、蘇惠珍(2004),三維地下水模式:PMWIN, 五南圖書出版股份有限公司。 7.王如意、易任(1997),應用水文學,茂昌圖書有限公司。 8.經濟部水利署水利規劃試驗所(2007),台中盆地地下水資 源利用調查評估(1/3)。 9.經濟部水利署水利規劃試驗所(2007),台中盆地做為地下 水庫可能性之探討專提報告書。 10.經濟部水利署(2008),『97年度地下水觀測站井規劃、 調查試驗及監造計畫』成果報告。 11.經濟部水利署水工試驗所(2006),台灣地區地下水文圖 圖集繪製工作(4/4)。 12.David Keith Todd﹐Larry W.Mays(2004)﹐Groundwater Hydrology, 3rd Edition.

簡 報 完 畢 敬 請 斧 正 45