地下水超抽引致地層下陷之三維解析研究

台灣地層下陷問題嚴重,致相關研究發展,已成為防治地層下陷之重要決策方向之一,惟台灣地區地下水含水層地質錯縱複雜,致數值模擬耗時而無法即時提供決策者所需之答案。故本研究方向乃發展一個簡捷之解析模式,俾供使用者能快速求解地下水超抽引致之地層下陷量,以預為地下水資源管理及水土資源保育分析之重要工具。根據相關文獻分析,了解目前對於水流耦合土體之模式尚未臻完整外,且多為數值解法,對於現實地層狀況諸多參數難以獲取及校定情況下,除模擬複雜困難外,模擬結果亦難以符合真實狀況。本研究目的依據台灣地層下陷防治之現況需求為導向,除藉由多孔彈性理論推導出土體變形與水體流動之三維耦合地層下陷方程式外,並尋求合理簡化分析方法,再推導出土體變形量與水壓非耦合方程式,以迅速反應計算地層下陷量。並藉由抽水行為探討與案例說明,使模式更臻完善。有關本研究成果可歸納臚列如下:一、成功推導出三維水流與土體耦合方程式,並再接續完成含孔隙水壓力與土體位移項之非耦合控制方程式,便於解決複雜問題及做為大型計算程式撰寫之支援。二、另研究成果於尋求一個可即時反應計算地層下陷量之簡化分析方法時,獲得一重大發現,即所得控制方程式型式與Stokes Flow相同,對於求解大區域三維之地層下陷問題係一突破。且相較於數值模擬可快速求得結果。另本研究除建立控制方程式外,對於含水層交界面之各類邊界條件,亦探討完成,極利於後續引申加值應用。三、含水層內之點抽水源,由本研究證明可藉由流體力學中無雷諾數黏流之Stokeslet現象來描述。並利用邊界積分法再次佐證為Radial-Stokeslet。四、地層內之水力傳導係數經由達西定律與本研究結果比對,顯現其非為定值,為距離之函數,其值將因土體位移而隨之改變,而改變值隨遠距離至抽水點之距離成反比現象。另因水力傳導係數值極小,乘上距離後,其改變值微小,惟已反應抽水後實際顆粒組成改變之物理意義。

作者:曾鈞敏