MARTHE (acronym for Modélisation d'Aquifères avec un maillage Rectangulaire, Transport et HydrodynamiquE) is designed for 2D or 3D modelling of flows and transfers in aquifer systems, including climatic, human influences and possible geochemical reactions.
MARTHE: Modelling software for groundwater flows
Model using the MARTHE code of flows in the northern Aquitaine aquifer system, comprising 15 aquifer layers:Saltel and Pédron (2012).
Model of the Poitou region produced with MARTHE.
Several views with the WinMarthe pre-processor. Top left: simulated piezometry and pathlines. Top right: 3D model of the Alsace plain using an irregular tartan grid containing a regular nested grid. Bottom right: vertical cross section showing bevelled layers3.
Complex hydrogeological configurations: a multilayer system extending across an entire region.
TCE saturation calculated to t = 180 (the colours represent computed saturation).
Extent of the TCE cloud observed (between glass plates) at t = 184.Kueper experiment.Infiltration of TCE (Tetra-Chloro-Ethylene), which is non miscible in initially water-saturated heterogeneous sands.
Diphase simulation of water injection at 4 points in an initially hydrocarbon-satured formation.The colours represent the percentage of water saturation calculated after 1209 days.
Methane saturation calculated after 8840 days.
Methane pressure (in m of water) calculated after 8840 days.Diphase simulation of methane storage in an aquifer.
pH profile according to distance from the injection point.
Variation in the volumes of certain minerals in the matrix. Injection of CO2-saturated water into the Dogger aquifer (Paris basin). Radial 2D grid. Geochemical REACT solver for the LBNL TOUGHREACT code (Berkeley) (MARTHE-TREACT coupling not available outside BRGM).
Description of the MARTHE computer code
MARTHE is an application for computing pollutant and energy flows and transport in porous environments (in aquifers) developed at the BRGM. This application, based on the finite volumes method, is a complete code incorporating:
- single-layer and multi-layer aquifers and 3D,
- hydrographic networks (rivers, other waterways, drainage),
- exchanges with the atmosphere (rain, snow, evapotranspiration).
The calculations may also include:
- mass transfer, for water quality protection,
- temperature effects,
- influence of salinity,
- degradation of pollutants,
- transfers through the vadose zone,
- geochemical interactions with Phreeqc or React modules.
MARTHE: areas of application
MARTHE is designed to model fluid flow and transfer issues in different contexts:
Management of aquifer resources
- Assessment of water budget components in an aquifer system: rainwater recharge, lateral input from catchment areas, underground circulation and associated flows, year-on-year fluctuations, seasonal storage and release,
- Hydrodynamic impacts of existing or planned structures: pumping, irrigation, drainage, gravel beds, retention basins,
- management and optimisation of catchments,
- studies of the influence of climate variability and predicted climate change impacts,
- diagnoses of groundwater pollution,
- Pollutant infiltration into the vadose zone, percolation to the aquifer followed by groundwater migration,
- simulation of a pollution plume released from a contaminated zone: pathlines, velocities, downstream concentrations.Modelling hydraulic scenarios for confinement or decontamination,
- impacts of domestic and industrial landfill on groundwaters,
- studies of underground storage confinement.
Civil engineering and mining
- dewatering of excavations
- hydraulic effects of impermeable walls,
- below-ground worksites (underground, car parks, tunnels),
- calculation of mine dewatering and corresponding drawdown.
- 1D, 2D or 3D ("full 3D") grids.
- Calculations by finite volumes on structured grids with possibilities for nested grids.
- Single or multilayer aquifers (layered aquifers, sometimes separated by a semi-permeable aquitard).
- Unconfined, confined or semi-confined aquifers in steady or non-steady state.
- Calculations may include discontinuities such as open water bodies (lakes, gravel pits), local aquifer depletion (and subsequent recharge), including in multilayer aquifers, groundwater flooding (rivers, springs, drains), impermeable walls (sheet pile walls, etc.).
- Automatic limitation of pumping rates as pump heads are uncovered.
- Calculation of well hydraulic head according to well diameter compared to the size of the grid cell containing it.
- Full coupling with hydrographic systems (rivers, drains).
- Coupled hydro-climatic calculations for rain, evapotranspiration, infiltration and runoff: the GARDENIA computing procedure.
- Horizontal and vertical anisotropism of permeable areas.
- Calculation of (direct and reverse) pathlines in steady or non-steady state.
- The field corresponding to each layer to be modelled is discretised into a structured rectangular tartan grid.
- The width of each line and column may differ with local heterogeneity, the density of available information and the level of accuracy required.
- The grid can be further refined locally with nested sub-grids.
- Grids can be in vertical cross-sections or radial (cylindrical coordinates).
- Several million cells can be handled, the maximum cell number depends purely on the physical memory installed.
- Maximum potential capacity:
- 3000 columns, 3000 rows, 999 layers, 99 nested layers.
- Unlimited calculation time steps.
- Hydrodispersive migration of an effuent in the aquifer and the vardose zone.
- Calculations use the TVD method ("Total Variation Diminishing"), the MOC method (Method of Characteristics) or the Finite Differences method.
- Chain-reaction degradation
- Retardation factor, kd partition coefficient (adsorption-desorption)
- Sorption by Langmuir or Freundlich isotherm.
- Possibilities for control of injection well concentration by production well concentration,
- Reactive transport with PHREEQC (USGS) geochemical modules.
Non Saturated Zone, Density effects, Temperature
- Continuous flow and transfer in the saturated zone and the vadose zone.
- Energy transfers, water temperature calculations.
- Density effects due to variations in salinity and/or water temperatures.
- Influence of temperature on viscosity and degradation coefficient taken into account.
- Possibility for control of injection well temperature by production well temperature (geothermal doublets).
Automatic calibration, sensitivity analyses
- Automatic calibration of model parameters, defined per homogeneous zone or per grid cell.
- Sensitivity analysis to the calibrated parameters.
- Vertical fracturing simulated by equivalent transmissivity.
- Semi-pervious vertical or horizontal fractures
- Gallery networks taken into account.
- Diphase flows with non-miscible phases: freshwater and salt water, water and air, water and "oil".
- Gas flows.
- Vegetation growth taken into account.
- Nitrate or nitrogen balance resulting from crops.
Visualisation, import and export
- Pre-processor and post-processor through graphic interface
- Detailed "help" file
- Input interface with 3D modellers:Multilayer ®, GDM ®, EarthVision ®, Éclipse ®, Surfer ®.
- Export to Mapinfo ®, Paraview, Tecplot ®, Winteracter-3Dview ®, Vrml QGis, ArcGis ® for 3D displays.
- Export to Mapinfo ®, QGis, ArcGis ®, Surfer ® for 2D displays.
- Export to Microsoft Excel ® to display series.
Calculations in non-steady state
Any grid-cell information can be modified at each calculation time step:
- hydraulic head, temperature, concentration, salinity,
- pumping and injection flow rates, rainfall, potential evapotranspiration, air temperature, recharge,
- permeability, porosity, riverbank clogging, extent of hydrographic network or drains,
- topographic elevation,
- setting or release of prescribed hydraulic heads,
Conditions / pricing (June 2021)
Including transport, temperature, density, non-saturated zone, PHREEQC geochemistry, etc.
- Single-user licence:7900 € (excl. VAT)
- Subsequent licences:4650 € (excl. VAT)
- Email assistance and maintenance (if required):15 % per year (*)
- Licences for universities and public research organisations: please contact us
(*) Email assistance: advice, and answers to specific questions, mainly for starting up the software: e.g. choosing modelling parameters, digital patterns and limit conditions; meanings of terms used in the files, etc.
Email assistance does not include questions concerning the site to be modelled, hydrogeological issues, etc.
Maintenance consists of supplying more recent upgraded versions (about one every 6 months).
Limited "free access" version
Calculation of hydraulic heads, flow, pollutant and heat transport and PHREEQC geochemistry only.
- Free licence
Email or telephone assistance is not provided for the free version.
Because MARTHE is a very comprehensive code, you will need a little time to master it. We therefore regularly organise training sessions at Orleans, for beginner and advanced levels.
- CEA - Cadarache
- CEA - IPSN
- EDF / SQR Aix en Provence
- LECES Consultancy
- ECOSYST Consultancy
- SOLVAY - Brussels
- URS France, Paris office
- SUEZ Environnement
- Société Anonyme des Eaux Minérales d'ÉVIAN. DANONE
- EGIS Structures & Environnement - Egis Géotechnique
- SÉVÊQUE ENVIRONNEMENT, Nord Pas de Calais agency
- Wroclaw Municipal Water and Sewage Company, Pologne
- CFG SERVICES
- Seine Normandie Water Agency - Direction Territoriale et Maritime des Rivières de Basse Normandie
- DIREN Nord-Pas-de-Calais
- Agence de Bassin Seine Normandie AESN Caen
- DREAL, Caen
- Lille Métropole Communauté urbaine (LMCU)
- Chambre d’Agriculture de l’Aisne
- Communauté urbaine du Grand Lyon
- Università degli Studi di Milano, Dipartimento di Scienze della Terra
- Chott-Mariem Higher School of Agronomy, Tunisia
- University of São Paulo, Brazil
- University of Pernambuco, Recife, Brazil
- University of Mons – Fundamental and Applied Geology, Belgium
- Institute of Geology, Estonia
- EGID Bordeaux 3
- ENSG - INPL
- UTC Compiègne
- IEEA Wroclaw University, Poland
- École des Mines de Paris / CIG
- École des Mines de Nancy / LAEGO
- University of Avignon and the Vaucluse
- University of Neuchâtel
- EPFL, Lausanne
- Andra Pradesh Groundwater Dept., India
- Ministry of Defence of Ukraine
- Institute of Geological Sciences of Ukraine
- Agence de Bassin Seine Normandie (AESN), Paris
- REGIDESO (Democratic Republic of Congo)
- SYSTRA - Mott MacDonald Paris
This list does not include organisations using the limited "free access" version of MARTHE.
Note: Dialog box, inputs and outputs of the software are in English; however the documentation and help files are in French, with the exception of the English version of the tutorial.