A numerical model is presented that solves the partial differential equations describing the motion of salt water and fresh water separated by a sharp interface. The areal equations are based on the Dupuit approximation and are obtained from partial integration over the vertical dimension. Finite-difference techniques are applied and the utility of several solution schemes is tested. The most efficient and accurate solution scheme uses block line-successive over-relaxation. Examples are given to: (1) test the model, (2) evaluate the Dupuit approximation, and (3) demonstrate the application to a field situation. The results show that the model is in good agreement with an analytical solution, but under severe conditions the Dupuit approximation may be inappropriate. The model is applied to a field area near Kahului, Maui, Hawaii and results extend the analysis of the problem beyond previous efforts.
[1]
Gedeon Dagan,et al.
Moving Interface in Coastal Aquifers
,
1964
.
[2]
K. H. Coats,et al.
Three-Dimensional Simulation of Steamflooding
,
1974
.
[3]
K. H. Coats,et al.
Direct methods in reservoir simulation
,
1974
.
[4]
J. Gillis,et al.
Matrix Iterative Analysis
,
1961
.
[5]
Charles R. Faust,et al.
Finite-difference model to simulate the areal flow of saltwater and fresh water separated by an interface
,
1980
.
[6]
Uri Shamir,et al.
Motion of the Seawater Interface in Coastal Aquifers: A Numerical Solution
,
1971
.
[7]
G. Forsythe,et al.
Computer solution of linear algebraic systems
,
1969
.