Étude et simulation d'un modèle stratigraphique advecto-diffusif non-linéaire avec frontières mobiles

Abstract

Retracer l’histoire d’un bassin est un préalable essentiel à toute recherche d’hydrocarbures. Pour cela, on a recours à un modèle stratigraphique, qui simule l'évolution des bassins sédimentaires sur de grandes échelles de temps (millions d'années) et d'espace (centaines de kilomètres). Le logiciel Dionisos, développé à IFPEN depuis 1992 et très apprécié par les compagnies pétrolières, permet d’effectuer ce type de calculs en prenant en compte deux grands processus physiques : (1) le transport gravitaire des sédiments dû à l’inclinaison du sol ; (2) l’écoulement de l’eau provenant des fleuves et des précipitations. Le transport gravitaire est décrit par une équation de diffusion dans laquelle le flux de sédiments dépend de la pente du sol. Initialement, cette dépendance est linéaire. Pour mieux s’approcher des observations réelles, on souhaite la rendre non-linéaire par l’intermédiaire d’un p-Laplacien. Ce changement nécessite la conception d’une nouvelle méthode de résolution numérique, qui doit offrir non seulement une grande rapidité d’exécution, mais aussi des garanties de robustesse et de précision des résultats. De plus, elle doit être compatible avec une contrainte sur le taux d’érosion présente dans le modèle. L’ajout de l’écoulement de l’eau est aussi une sophistication récente du modèle physique de Dionisos. Il se traduit par l’introduction d’une nouvelle équation aux dérivées partielles, couplée à celle du transport. Là encore, il est important d’élaborer une stratégie de résolution numérique innovante, en ce sens qu’elle doit être à la fois performante et bien adaptée au fort couplage de ces deux phénomènes. L'objectif de cette thèse est de moderniser le cœur numérique de Dionisos afin de traiter plus adéquatement les processus physiques ci-dessus. On cherche notamment à élaborer un schéma implicite par rapport à toutes les inconnues qui étend et améliore le schéma actuel. Les méthodologies retenues serviront de base à la prochaine génération du calculateur.