• 沒有找到結果。

數值試驗流程與案例設計

第四章 數值試驗

4.1 數值試驗流程與案例設計

準備觀測資料與 設定參數目標值

建立待檢定案例

分析測試結果 測試檢定

案例

建立一地下水 模式A (真值系統)

建立與模式A完 全相同之模式B

修改模式B中待 檢定參數值

將模式A之模擬 水位設定為模 式B之觀測水位

使用本參數檢定系統 對模式B進行參數檢定

將模式B檢定後之參 數與模式A之參數

(目標值)比較

輸出模式A模 擬水位

圖4.1-1 地下水模式參數檢定之正確性驗證流程

數值試驗的流程圖如 4.1-1 所示,無論是 DDS 與 MDDS 版本的 地下水模式參數檢定系統,均須透過數值試驗檢驗系統之正確性,待 確認兩系統均可完成檢定後,接著再進一步較兩者之檢定成果,以驗 證MDDS 確實可改善 DDS 檢定地下水模式之效率。數值試驗之流程 分為四個步驟:第一步驟為建立「真值系統」,此真值系統概念上代表 真實世界之地下水系統,因此真值系統之模擬水位將作為現地之觀測 水位。第二步驟則是建立「待檢定系統」,由於本研究著重於穩態模式 之淨補注量檢定,故待檢定系統中除了淨補注量以隨機方式給定外,

其他參數均直接給予真值系統之參數值。第三步驟即使用本研究建立 之參數檢定系統,對待檢定系統進行參數檢定,其檢定依據為真值系 統之模擬水位。第四步驟則是檢驗參數檢定成果,比較待檢定系統與 真值系統之淨補注量差異。

在地下水流模式的設計上,本研究設定為穩態地下水模擬,模擬 區域為 11(km)×11(km),以 1(km)×1(km)見方切割整個模擬區域,共 121 個網格,如圖 4.1-2 所示。此外,含水層共分為三層,如圖 4.1-3 所示,第一層及第三層為含水層,第二層為阻水層,但阻水層在空間 上並未涵蓋整個模擬區域,有部分區域之透水性近似含水層。故在模 擬上,第一層之含水層型態為非受壓水層(unconfined aquifer),第二、

393 個(11×11×3)之正方形網格。

不 透 水 邊 界

不透水邊界 不透水邊界

定 水 頭 邊 界

圖4.1-2 設計案例之模擬網格

定水頭

邊界 第一分層

第二分層 第三分層 含水層1

含水層2 阻水層1

圖 4.1-3 設計案例之概念分層

在邊界條件方面,第一分層的左邊界設定為定水頭邊界且水位高 在 0 公尺,用以模擬大海,右方、上方及下方為不透水邊界;第二、

三分層四個方向均為不透水邊界;在分區方面,本案例在第一層及第 三層均分為 9 個分區,共有 18 個分區,每個分區內皆為 9 個網格,

如圖4.1-4、4.1-5 與 4.1-6 所示。在地表高程方面,地勢由右向左遞減,

由120 公尺逐漸降至 0 公尺高。含水層厚度為 25 至 60 公尺不等,阻

水層厚度約為5 公尺。在初始條件方面,整體模擬區域初始水位高皆

K:0.004 K:0.004

K:0.004

K:0.004

K:0.004 K:120

K:0.004 K:0.004 K:0.004

圖4.1-5 第二分層邊界及 K 值(公尺/天)配置圖