Modèles mathématiques optimisent l’exploitation des réservoirs géothermiques

Contexte géothermique

Le gradient géothermique moyen est d’environ 2,5–3 °C par 100 m de profondeur, ce qui signifie que la température augmente de plusieurs centaines de degrés à quelques kilomètres sous la surface. Dans les zones de forte activité tectonique, comme le « Ring of Fire », ces gradients sont supérieurs, créant des réservoirs de chaleur exploitable. À 6 000 km de profondeur, le manteau atteint près de 5 000 °C, mais seule une fraction de cette énergie remonte à la surface sous forme de vapeur, alimentant geysers et sources chaudes. Ces caractéristiques physiques imposent des contraintes thermiques et mécaniques que les modèles doivent intégrer.

Modélisation mathématique des réservoirs

Les ingénieurs utilisent des modèles Lagrangien‑Eulerien pour suivre le déplacement du fluide à travers les milieux poreux du sous‑sol. Ce cadre combine la description des particules (Lagrangien) avec les champs de vitesse (Eulerien), permettant de simuler la convection et la diffusion de chaleur dans les roches déformables. En parallèle, des modèles stochastiques optimisent les stratégies d’injection et de production en évaluant les incertitudes liées aux propriétés de perméabilité. Les modèles géométriques quantifient la forme des fractures et des conduits, essentielles pour estimer la connectivité du réservoir. Un exemple concret est l’extension du modèle MODFLOW, publié en 1984, qui a été adapté pour intégrer les échanges thermiques et les contraintes mécaniques propres aux systèmes géothermiques.

Applications et limites techniques

Les pompes à chaleur géothermiques exploitent le gradient thermique en transférant la chaleur du sol vers les bâtiments en hiver, et inversement en été. Cette technologie, développée dès les années 1920 en Islande, repose sur des boucles enterrées au‑delà de la ligne de gel, où la température reste stable. Cependant, le pompage à grande échelle modifie la pression du sous‑sol, comme l’a montré le cas de Venise où le retrait d’eau souterraine a entraîné un affaissement du sol. De même, l’injection massive d’eau dans un champ géothermique peut provoquer des variations de porosité, affectant la productivité à long terme. Ainsi, la précision des modèles mathématiques est cruciale pour anticiper ces effets et limiter les risques de subsidence ou de perte de capacité de génération d’électricité.