Modflow计算地下水流程式原理.doc

1.1 Modflow計算地下水流程式原理 　　運用PMWIN(processing modflow)程式中之Modflow來地下水模擬，模式採用有線差分法(FDM)來求得近似數值解。模擬之前需將模擬區域劃分成網格(Cell)，每一網格假設為均勻物質，以網格節點或網格中心之水頭代表網格水頭，含水層厚度在水平方向上可以變化，但在數值上會有較大的誤差。 　　在Modflow中使用的固定密度(Constant Density)地下水流偏微分方程式(McDonald and Harbaugh，1988，p.2-1)如下所示：　　　　　　式1-1 式中、、分別為沿著x、y、z座標軸的水力傳導係數。 　　h=壓力水頭(Potentiometric Head) 　　WW0.0表示流出地下水系統(Out Flow of Ground-Water System)，反之，W0.0表示流進地下水系統。 　　Ss=孔隙介質之儲水率。　　t = 時間(T)。 　　當式1-1結合邊界及初始條件時，描述在非均質性及非等向性介質(提供主軸方向的水力傳導係數)中的暫態三維地下水流。左側代表含水層在不同方向之入滲與水頭變化及地下水補給與抽水造成之三向度地下水流動模式；在均向性穩定流含水層系統中，因為流率為一常數，多孔介質之儲水量為定值，此式即為Laplace方程式，可由給定之含水層流況、初始水位狀況與邊界條件獲得