建議

1. 本研究方法目前僅應用於拘限含水層之模擬，但搭配徐昇氏網格 之模擬方式並非僅可應用於此，建議後續研究可以擴充至非拘限 含水層、未飽和問題等模擬。

2. 本數值模式內已經建立 Rana A. Fine(1973)提出之水流密度計算公 式，可以同時針對壓力與溫度計算對應的水流密度，然目前僅針 對水流議題進行探討，故僅固定為攝氏 20 ，建議未來可以搭配 熱流模式，使得水流密度亦可隨模擬溫度而變化。

oC

3. 土壤壓密係數 是用以呈現土體隨壓力變化之壓縮性參數，本研 究目前均設定為 0，建議未來可透過實驗或其他方法，進一步檢驗 各類土壤之係數

α/

α/

參考文獻

1. McDonald, M.G., and Harbaugh, A.W., (1988) A modular three-dimensional finite-difference ground-water flow model: U.S.

Geological Survey Techniques of Water-Resources Investigations, book 6, chap. A1, 586 p.

2. Trescott, P.C., Pinder, G.F., and Larson, S.P., (1976) Finite-difference model for aquifer simulation in two dimensions with results of numerical experiments: Techniques of Water-Resources Investigations of the United States Geological Survey, book 7, chap.

C1, 116 p.

3. Larson, S.P., (1978) Direct solution algorithm for the two-dimensional ground-water flow model: U.S. Geological Survey Open-File Report 79-202, 22 p.

4. Manteuffel, T.A., Grove, D.B., and Konikow, L.F., (1983) Application of the conjugate-gradient method to ground-water models: U.S. Geological Survey Water-Resources Investigations Report 83-4009, 24 p.

5. Ozbilgin, M.M., and Dickerman, D.C., (1984) A modification of the finite-difference model for simulation of two-dimensional ground-water flow to include surface-ground water relationships: U.S.

Geological Survey Water-Resources Investigations Report 83-4251, 98 p.

6. Weeks, J.B., (1978) Digital model of ground-water flow in the Piceance Basin, Rio Blanco and Garfield Counties, Colorado: U.S.

Geological Survey Water-Resources Investigations 78-46, 108 p.

7. Trescott, P.C., (1975) Documentation of finite-difference model for

simulation of three-dimensional ground-water flow: U.S. Geological Survey Open-File Report 75-438, 32 p.

8. Trescott, P. C., and Larson, S. P., (1976) supplement to Open-File Report 75-438--Documentation of finite-difference model for simulation of three-dimensional ground-water flow: U.S. Geological Survey Open-File Report 76-59l, 21 p.

9. Kuiper, L.K., (1987) Computer program for solving ground-water flow equations by the preconditioned conjugate gradient method: U.S.

Geological Survey Water-Resources Investigations Report 87-4091, 24 p.

10. Leake, S.A., and Prudic, D.E., (1991) Documentation of a computer program to simulate aquifer-system compaction using the modular finite-difference ground-water flow model: Techniques of Water-Resources Investigations of the United States Geological Survey, book 6, chap. A2, 68 p.

11. Grove, D.B., and Stollenwerk, K.G., (1984) Computer model of one-dimensional equilibrium-controlled sorption processes: U.S.

Geological Survey Water-Resources Investigations Report 84-4059, 58 p.

12. Konikow, L.F., and Bredehoeft, J.D., (1978) Computer model of two-dimensional solute transport and dispersion in ground water:

Techniques of Water-Resources Investigations of the United States Geological Survey, book 7, chap. C2, 90 p.

13. Goode, D.J., and Konikow, L.F., (1989) Modification of a method-of-characteristics solute-transport model to incorporate decay and equilibrium-controlled sorption or ion exchange: U.S. Geological

Survey Water-Resources Investigations Report 89-4030, 65 p.

14. Dunlap, L.E., Lindgren, R.J., and Carr, J.E., (1984) Projected effects of ground-water withdrawals in the Arkansas River Valley, 1980-99, Hamilton and Kearny Counties, southwestern Kansas: U.S.

Geological Survey Water-Resources Investigations Report 84-4082, 168 p.

15. Lewis, F.M., Voss, C.I., and Rubin, Jacob, (1986) Numerical simulation of advective-dispersive multisolute transport with sorption, ion exchange and equilibrium chemistry: U.S. Geological Survey Water-Resources Investigations Report 86-4022, 165 p.

16. Glover, K. C., (1987) A dual-porosity model for simulating solute transport in oil shale: U.S. Geological Survey Water Resources Investigations Report 86-4047, 88 p.

17. Cooley, R.L., (1992) Derivation of finite-element equations and comparisons with analytical solutions; Part 2 of A MODular Finite-Element model (MODFE) for areal and axisymmetric ground-water-flow problems: U.S. Geological Survey Techniques of Water-Resources Investigations, book 6, chap. A4, 108 p.

18. Torak, L.J., (1993) Design philosophy and programming details; Part 3 of A MODular Finite-Element model (MODFE) for areal and axisymmetric ground-water-flow problems: U.S. Geological Survey Techniques of Water-Resources Investigations, book 6, chap. A5, 243 p.

19. Lin, H. C., D. R. Richards, G. T. Yeh, J. P. Cheng and N. L. Jones.

(1996). "FEMWATER: a three-dimensional finite element computer model for simulating density dependent flow and transport.", US

Army Engineer Waterways Experiment Station, Pennsylvania State University, Brighan Young University.

20. J. N. Portela, and M. S. Alencar, (2004) “Outage contour using. a Voronoi diagram”, IEEE, Piscataway NJ, ETATS-UNIS.

21. A.H. Thiessen., (1911) Precipitation averages for large areas.

Monthly Weather Review, 39(7): 1082-1084

22. N. Sukumar, (2003) Voronoi cell finite difference method for the diffusion operator on arbitrary unstructured grids. International Journal for Numerical Methods in Engineering, 57(1):1-34.

23. Tucker, G.E., Lancaster, S.T., Gasparini, N.M., Bras, R.L., and Rybarczyk, S.M, (2001) An Object-Oriented Framework for Hydrologic and Geomorphic Modeling Using Triangulated Irregular Networks, Computers and Geosciences, 27(8), pp. 959-973.

24. A. Iske and M. Käser, (2004) Conservative semi-Lagrangian advection on adaptive unstructured meshes. Numerical Methods for Partial Differential Equations, 20:388–411.

25. J. Fuhrmann and H. Langmach, (1998) Stability and existence of solutions of time-implicit finite volume schemes for viscous nonlinear conservation laws. Preprint 437, Weierstraß-Institut, Berlin.

26. Georgy Voronoi, (1907) Nouvelles applications des paramètres continus à la théorie des formes quadratiques. Journal für die Reine und Angewandte Mathematik, 133:97-178.

27. Andreas Fabri, (2007) Voronoi Diagrams in CGAL, the Computational Geometry Algorithms Library. Voronoi Diagrams in Science and Engineering, 2007. ISVD '07. 4th International Symposium on 9-11 July 2007 Page(s):8 - 14

28. Steven Fortune, (1986) A sweepline algorithm for Voronoi diagrams, Algorithmica, 2:153-174, 1987.

29. van Genuchten, M.Th, (1980) A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Am. J.

44:892-898.

30. R.A. Fine and F.J. Millero, (1973) The compressibility of water as a function of temperature and pressure. Journal Chemical Physics, 59: 5529-5536.

31. 黃宏明(2005) Voronoi Diagram 互動教學程式之設計與研究，私立 中華大學碩士論文。

附錄 土壤壓密係數與傳統方法之關聯推導

式 A-6 為本模式描述控制體積內蓄水質量之方程式，其中 代表

層因應壓力水頭變化，土體與水體因壓縮或膨脹所釋放出的水量，式