Parameter Identification by High-Resolution Inverse Numerical Model Based on LBM/CMA-ES: Application to Chalk Aquifer (North of France)
The present paper proposes the numerical solution of an inverse problem in groundwater flow (Darcy’s equation). This solution was achieved by combining a high-resolution new code HYSFLO-LBM (Hydrodynamic of Subsurface Flow by Lattice Boltzmann Method), based on LBM, to solve the direct problem, and...
| Main Authors: | , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
MDPI AG
2021-06-01
|
| Series: | Water |
| Subjects: | |
| Online Access: | https://www.mdpi.com/2073-4441/13/11/1574 |
| _version_ | 1797531623921549312 |
|---|---|
| author | Lahcen Zouhri Sami Kaidi Hassan Smaoui |
| author_facet | Lahcen Zouhri Sami Kaidi Hassan Smaoui |
| author_sort | Lahcen Zouhri |
| collection | DOAJ |
| description | The present paper proposes the numerical solution of an inverse problem in groundwater flow (Darcy’s equation). This solution was achieved by combining a high-resolution new code HYSFLO-LBM (Hydrodynamic of Subsurface Flow by Lattice Boltzmann Method), based on LBM, to solve the direct problem, and the metaheuristic optimization algorithm CMA-ES ES (Covariance Matrix Adaptation-Evolution Strategy) to solve the optimization step. The integrated optimization algorithm which resulted from this combination, HYSFLO-LBM/CMA-ES, was applied to the hydrogeological experimental site of Beauvais (Northern France), instrumented by a set of sensors distributed over 20 hydrogeological wells. Hydrogeological parameters measured by the sensors are necessary to understand the aquifer functioning and to serve as input data for the identification of the transmissivity field by the HYSFLO-LBM/CMA-ES code. Results demonstrated an excellent concordance between the integrated optimization algorithm and hydrogeological applied methods (pumping test and magnetic resonance sounding). The spatial distribution of the transmissivity and hydraulic conductivity are related to the heterogeneous distribution of aquifer formations. The LBM and CMA-ES were chosen for their proven excellent performance and lesser cost, in terms of both money and time, unlike the geophysical survey and pumping test. The model can be used and developed as a decision support tool for integrated water resources management in the region. |
| first_indexed | 2024-03-10T10:47:19Z |
| format | Article |
| id | doaj.art-eb80e1a14ff14b498d3119db4c4b60c4 |
| institution | Directory Open Access Journal |
| issn | 2073-4441 |
| language | English |
| last_indexed | 2024-03-10T10:47:19Z |
| publishDate | 2021-06-01 |
| publisher | MDPI AG |
| record_format | Article |
| series | Water |
| spelling | doaj.art-eb80e1a14ff14b498d3119db4c4b60c42023-11-21T22:34:03ZengMDPI AGWater2073-44412021-06-011311157410.3390/w13111574Parameter Identification by High-Resolution Inverse Numerical Model Based on LBM/CMA-ES: Application to Chalk Aquifer (North of France)Lahcen Zouhri0Sami Kaidi1Hassan Smaoui2AGHYLE, Institut Polytechnique UniLaSalle Beauvais, SFR Condorcet FR CNRS 3417, 19 Rue Pierre Waguet, BP 30313, CEDEX, 60026 Beauvais, FranceCEREMA Risques Eaux et Mer (REM) EPR HA, 134 Rue de Beauvais, 60280 Margny-Les-Compiègne, FranceCEREMA Risques Eaux et Mer (REM) EPR HA, 134 Rue de Beauvais, 60280 Margny-Les-Compiègne, FranceThe present paper proposes the numerical solution of an inverse problem in groundwater flow (Darcy’s equation). This solution was achieved by combining a high-resolution new code HYSFLO-LBM (Hydrodynamic of Subsurface Flow by Lattice Boltzmann Method), based on LBM, to solve the direct problem, and the metaheuristic optimization algorithm CMA-ES ES (Covariance Matrix Adaptation-Evolution Strategy) to solve the optimization step. The integrated optimization algorithm which resulted from this combination, HYSFLO-LBM/CMA-ES, was applied to the hydrogeological experimental site of Beauvais (Northern France), instrumented by a set of sensors distributed over 20 hydrogeological wells. Hydrogeological parameters measured by the sensors are necessary to understand the aquifer functioning and to serve as input data for the identification of the transmissivity field by the HYSFLO-LBM/CMA-ES code. Results demonstrated an excellent concordance between the integrated optimization algorithm and hydrogeological applied methods (pumping test and magnetic resonance sounding). The spatial distribution of the transmissivity and hydraulic conductivity are related to the heterogeneous distribution of aquifer formations. The LBM and CMA-ES were chosen for their proven excellent performance and lesser cost, in terms of both money and time, unlike the geophysical survey and pumping test. The model can be used and developed as a decision support tool for integrated water resources management in the region.https://www.mdpi.com/2073-4441/13/11/1574groundwatermodelingLattice Boltzmann methodchalk aquifertransmissivityheterogeneity |
| spellingShingle | Lahcen Zouhri Sami Kaidi Hassan Smaoui Parameter Identification by High-Resolution Inverse Numerical Model Based on LBM/CMA-ES: Application to Chalk Aquifer (North of France) Water groundwater modeling Lattice Boltzmann method chalk aquifer transmissivity heterogeneity |
| title | Parameter Identification by High-Resolution Inverse Numerical Model Based on LBM/CMA-ES: Application to Chalk Aquifer (North of France) |
| title_full | Parameter Identification by High-Resolution Inverse Numerical Model Based on LBM/CMA-ES: Application to Chalk Aquifer (North of France) |
| title_fullStr | Parameter Identification by High-Resolution Inverse Numerical Model Based on LBM/CMA-ES: Application to Chalk Aquifer (North of France) |
| title_full_unstemmed | Parameter Identification by High-Resolution Inverse Numerical Model Based on LBM/CMA-ES: Application to Chalk Aquifer (North of France) |
| title_short | Parameter Identification by High-Resolution Inverse Numerical Model Based on LBM/CMA-ES: Application to Chalk Aquifer (North of France) |
| title_sort | parameter identification by high resolution inverse numerical model based on lbm cma es application to chalk aquifer north of france |
| topic | groundwater modeling Lattice Boltzmann method chalk aquifer transmissivity heterogeneity |
| url | https://www.mdpi.com/2073-4441/13/11/1574 |
| work_keys_str_mv | AT lahcenzouhri parameteridentificationbyhighresolutioninversenumericalmodelbasedonlbmcmaesapplicationtochalkaquifernorthoffrance AT samikaidi parameteridentificationbyhighresolutioninversenumericalmodelbasedonlbmcmaesapplicationtochalkaquifernorthoffrance AT hassansmaoui parameteridentificationbyhighresolutioninversenumericalmodelbasedonlbmcmaesapplicationtochalkaquifernorthoffrance |