CWEEE  Vol.11 No.1 , January 2022
Modeling of the Saltwater Intrusion Using the Level Set Method. Application to Henry’s Problem
Abstract: The salt intrusion phenomenon is caused by overexploitation of aquifers in coastal areas. This physical phenomenon has been the subject of numerous studies and numerous methods have been proposed, with the aim of protecting the quality of the water in these aquifers. This work proposes a two-dimensional saline intrusion model using the sharp interface approach and the level set method. It consists of a parabolic equation modeling the underground flow and a hyperbolic Equation (the level set equation) which makes it possible to track the evolution of the interface. High-order numerical schemes such as the space scheme WENO5 and the third-order time scheme TVD-RK were used for the numerical resolution of the hyperbolic equation. To limit the tightening of the contour curves of the level set function, the redistanciation or reinitialization algorithm proposed by Sussma et al. (1994) was used. To ensure the effectiveness and reliability of the proposed method, two tests relating to the standard Henry problem and the modified Henry problem were performed. Recall that Henry’s problem uses the variable density modeling approach in a confined and homogeneous aquifer. By comparing the results obtained by the level set method with reinitialization (LSMR) and those obtained by Henry (1964), and by Simpson and Clement (2004), we see in the two test cases that the level set method reproduces well the toe, the tip and the behaviour of the interface. These results correspond to the results obtained by Abarca for Henry’s problem with constant dispersion coefficients. The results obtained with LSMR, reproduced the interface with a slight spacing compared to those obtained by Henry. According to Abarca (2006), this spacing is due to the absence of the longitudinal and transversal dispersion coefficients in the model.
Cite this paper: Loua-Bouayi, J. , Tathy, C. and Manounou, A. (2022) Modeling of the Saltwater Intrusion Using the Level Set Method. Application to Henry’s Problem. Computational Water, Energy, and Environmental Engineering, 11, 11-33. doi: 10.4236/cweee.2022.111002.

[1]   Bear, J., Cheng, H., Sorek, S., Ouazar, D. and Herrera (1999) Seawater Intrusion in Coastal Aquifers Concepts, Methods and Practices. Kluwer Academic Publishers, Dordrecht.

[2]   Werner, A., Bakker, M., Post, V., Alexander, V., Lu, C., Behzad, A.A., Craig, S. and Barry, D. (2013) Seawater Intrusion Processes, Investigation and Management: Recent Advances and Future Challenges. Advances in Water Resources, 51, 3-26.

[3]   Mohamed, B., Margaret, S., Ramachandran, G. and Shelat, K. (2014) Vulnerability of Agriculture, Water and Fisheries to Climate Change.

[4]   Alix, L. (2006) Current Variations in Sea Level: Observations and Climatic Causes. Laboratory of Geophysics and Oceanographic Space, Toulouse Meteorology, 59.

[5]   Solomon, S., Qin, D., et al. (2007) The Scientific Elements. Contribution of Working Group I to Fourth Assessment Report of the Intergovernmental Panel on Climate Change. GIES, Climatic Changes.

[6]   Cheng, A. and Ouazar, D. (2004) Coastal Aquifer Management. Monitoring, Modeling, and Case Studies.

[7]   Bear, J. and Cheng, A. (2010) Modeling Groundwater Flow and Contaminant Transport. Part of the Theory and Applications of Transport in Porous Media Book Series (TATP, Volume 23).

[8]   Sherif, M. (1999) The Nile Delta Aquifer in Egypt, Chapter 17. In: Seawater Intrusion in Coastal Aquifers: Concepts, Methods and Practices, Book Series: Theory and Application of Transport in Porous Media, Kluwer Academic Publishers, Dordrecht, Vol. 14, 559-590.

[9]   Rifai, M., et al. (1956) Dispersion Phenomena in Laminar Flow through Porous Media. University of California Institute of Engineering, San Diego, Progress Report 2.

[10]   Essaid, H. (1986) A Comparison of the Coupled Fresh Water—Salt Water Flow and the Ghyben-Herzberg Sharp Interface Approaches to Modeling of Transient Behavior in Coastal Aquifer Systems. Journal of Hydrology, 86, 169-193.

[11]   Ogata, A. (1970) Theory of Dispersion in Granular Medium. U.S. Geological Survey Professional Paper 411-I.

[12]   Hubbert, M. (1940) The Theory of Ground Water Motion. Journal of Geology, 48, 785-944.

[13]   Sark, M. (1980) Seawater Intrusion and Potential Yield of Aquifers in the Soquel Aptos Area, Santa Cruz Country, California. Geological Survey, Water Resources Investigations, 80-84.

[14]   Badon-Ghyben (1888) Nota in verband met de voorgenomen putboring nabij Amsterdam (Notes on the Probable Results of Well Drilling near Amsterdam). Tijdschrift van het Koninklijk Institut van Ingenieurs, The Hague, 8-22.

[15]   Herzberg, A. (1901) Die Wasserversorgung einiger Nordseebader (The Water Supply of Parts of the North Sea Coast in Germany). Journal Gasbeleuchtung u. Wasserversorgung, 44, 815-819.

[16]   Henry, H. (1959) Salt Intrusion into Fresh-Water Aquifers. Journal of Geophysical Research, 64, 1911-1919.

[17]   Bear, J. and Verruijt, A. (1987) Modeling Groundwater Flow and Pollution. Theory and Applications of Transport in Porous Media.

[18]   Cooper, H. (1959) A Hypothesis Concerning the Dynamic Balance of Fresh Water and Salt Water in a Coastal Aquifer. Journal Geophysics Resources, 64, 461-467.

[19]   Fetter, C. (1972) Position of the Saline Water Interface beneath Oceanic Islands. Water Resources Research, 8, 1307-1315.

[20]   Mualem, Y. and Bear, J. (1974) The Shape of the Interface in Steady Flow in a Stratified Aquifer. Water Resources Research, 10, 1207-1215.

[21]   Amirat, Y., Hamdache, K. and Ziani, A. (1996) Mathematical Analysis for Compressible Miscible Displacement Models in Porous Media. Mathematical Models and Methods in Applied Science, 6, 729-747.

[22]   Choquet, C. (2010) Parabolic and Degenerate Parabolic Models for Pressure-Driven Transport Problems. Mathematical Models and Methods in Applied Science, 20, 543-566.

[23]   Mory, M. (2015) Mixed Net/Diffuse Interface Approach for Saline Intrusion Problems in Basement: Modeling, Mathematical Analysis and Digital Illustrations. PhD Dissertation, Rochelle University, Rochelle.

[24]   Abudawia, A. (2015) Numerical Analysis of a Finite Element Approximation for an Intrusion Model Saline in Coastal Aquifers. Ph.D. Thesis, Littoral Côte d’Opale University, Dunkerque.

[25]   Shamir, U. and Dagan, G. (1971) Motion If the Seawater Interface in Costal Aquifers: A Numerical Solution. Water Resources Research, 7, 644-657.

[26]   Anderson, M. (1976) Unsteady Groundwater Flow beneath Strip Oceanic Islands. Water Resources Research, 20, 640-644.

[27]   Sethian, J. (1996) Level Set Methods. Evolving Interface in Geometry, Fluid Mechanics, Computer Vision, and Materials Science. Cambridge Monographs on Applied and Computational Mathematics.

[28]   Osher, S. and Fedkiw, R. (2003) Level Set Methods and Dynamic Implicit Surfaces. Applied Mathematical Sciences. Springer, Berlin.

[29]   Zhu, J. and Sethian, J. (1992) Projection Methods Coupled to Level Set Interface Techniques. Journal of computational Physics, 102, 128-138.

[30]   Osher, S. and Sethian, J. (1988) Fronts Propagating with Curvature-Dependent Speed: Algorithms Based on Hamilton-Jacobi Formulations. Journal of Computational Physics, 79, 12-49.

[31]   Narkbuakaew, W., Nagahashi, H., Aoki, K. and Kubota, Y. (2014) An Efficient Liver-Segmentation System Based on a Level-Set Method and Consequent Processes. Journal of Biomedical Science and Engineering, 7, 994-1004.

[32]   Segol, G., Pinder, G. and William, G. (1975) A Galerkin-Finite Element Technique for Calculating the Transient Position of the Saltwater Front. Water Resources Research, 11, 343-347.

[33]   Hamidi, M., Reza, S. and Yazdi, S. (2006) Numerical Modeling of Seawater Intrusion in Coastal Aquifer Using Finite Volume Unstructured Mesh Method. Proceedings of the 9th WSEAS International Conference on Applied Mathematics, Istanbul, 27-29 May 2006, 572-579.

[34]   Sethian, J. and Smereka, P. (2003) Level Set Methods for Fluid Interfaces. Annual Review of Fluid. Mechanics, 35, 341-372.

[35]   Garzon, M., Sethian, J., et al. (2000) Numerical Simulation of Non-Viscous Liquid Pinch-Off Using a Coupled Level Set Boundary Integral Method. Journal of Computational Physics, 228, 6079-6106.

[36]   Thibaut, M. (2007) Development of a Level Set Method for Interface Monitoring, Application to Liquid Jet Rupture. Ph.D. Thesis, Rouen University, Rouen.

[37]   Sussman, M., Smereka, P. and Osher, S. (1994) A Level Set Approach for Computing Solutions to Incompressible Two-Phase Flows. Journal of Computational Physics, 114, 146-159.

[38]   Brackbill, J., Kothe, D. and Zemach, C. (1992) A Continuum Method for Modeling Surface Tension. Journal of Computational Physics, 100, 335-354.

[39]   Sefollahi, M. and Shirani, E. (2004) An Improved Method for Calculation of Interface Pressure Force in PLICVOF Methods. Abdus Salam International Centre for Theoretical Physics, Trieste.

[40]   Boliveau, A., Fortin, A. and Demay, Y. (1998) A Two-Dimensional Numerical Method for the Deformation of Drops with Surface Tension. International Journal of Computational Fluid Dynamics, 10, 225-240.

[41]   Popinet, S. (2018) Numerical Models of Surface Tension. Annual Review of Fluid Mechanics, 50, 49-75.

[42]   Chi-Wang (1997) Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws. Institute for Computer Application in Science and Engineering. NASA Langley Research Center, Hampton, Report, 97-65.

[43]   Sébastien, T. (2004) Development of an Interface Monitoring Method. Flow Applications Diphasic. Ph.D. Thesis, Rouen University, Rouen.

[44]   Harten, A., Engquist, B., Osher, S. and Chakravarthy, S. (1987) Uniformly High-Order Accurate Essentially Non-Oscillatory Schemes III. Hydrology. Journal of Computational Physics, 71, 231-303.

[45]   Biswarup, B. and Kumar, R. (2007) Accuracy Preserving ENO and WENO Schemes Using Novel Smoothness Measurement. Abdus Salam International Center for Theoretical Physics, Trieste.

[46]   Henry, H. (1964) Effects of Dispersion on Salt Encroachment in Coastal Aquifers, Sea Water in Coastal Aquifers. U.S. Geological Survey. Water Supply, 1613-C.

[47]   Frind, E. (1982) Simulation of Long-Term Transient Density-Dependent Transport in Groundwater. Advances in Water Resources, 5, 73-88.

[48]   Huyakorn, P., Andersen, P., Mercer, J. and White, H. (1987) Saltwater Intrusion in Aquifers: Development and Testing of a Three-Dimensional Finite Element Model. Water Resources Research, 23, 293-312.

[49]   Voss, C. and Souza, W. (1987) Modeling a Regional Aquifer Containing a Narrow Transition between Freshwater and Saltwater Using Solution Transport Simulation: Part I Theory and Methods. Geological Survey, U.S., 493514.

[50]   Croucher, A. and O’Sullivan, M. (1995) The Henry Problem for Saltwater Intrusion. Water Resources Research, 31, 1809-1814.

[51]   Thomas, A., Eldho, T. and Rastogi, A. (2016) Simulation of Seawater Intrusion in Coastal Confined Aquifer Using a Point Collocation Method Based Meshfree Model. Journal of Water Resource and Protection, 8, 534-549.

[52]   Lee, C. and Cheng, R. (1974) On Seawater Encroachment in Coastal Aquifers. Water Resources Research, 10, 1039-1043.

[53]   Simpson, M. and Clement, P. (2004) Improving the Worthiness of the Henry Problem as a Benchmark for Density-Dependent Groundwater Flow Models. Water Resources Research, 40, W01504.

[54]   Abarca, E., Carrera, J., Xavier, S.V. and Dentz, M. (2006) Anisotropic Dispersive Henry Problem. Advances in Water Resources, 30, 913-926.