Procédé (400), mis en œuvre par ordinateur, de correction de l’invalidité d’un maillage non-structuré volumique dans le contexte de l’insertion d’un modèle géologique (100) dans un maillage non-structuré volumique valide à mailles en forme de tétraèdres, dit maillage tétraédrique, le modèle géologique (100) étant décrit par une pluralité de données géométriques qui incluent des données de points, chaque maille (300) du maillage tétraédrique présentant une région intérieure délimitée par des éléments de forme (310, 320, 330) qui incluent des sommets (310), des arêtes (320) et des faces (330), l’insertion de chaque point incluant : - l’identification d’une seule maille (300), dite maille d’accrochage, dans laquelle se trouve géographiquement le point, - le calcul d’une distance intérieure de la maille d’accrochage (300), entre le point et des éléments de forme (310, 320, 330) de la maille d’accrochage (300), - l’identification d’un élément de forme (310, 320, 330) dit élément de forme à découper, qui se trouve à une distance intérieure qui est inférieure à une valeur de tolérance prédéterminée, et - l’application d’un opérateur de modification topologique, sur l’élément de forme à découper, au niveau du point, de manière à subdiviser la maille d’accrochage (300) en une pluralité de sous-mailles, le procédé (400) comprenant : - une étape de détection (410), par l’arithmétique exacte, d’au moins une sous-maille dont le volume n’est pas strictement positif, dite sous-maille invalide (300a, 300b, 300c), - une étape de correction locale (420) d’au moins une sous-maille invalide (300a, 300b, 300c) qui comprend, - une étape de calcul (421) d’une distance intérieure de la sous-maille invalide (300a, 300b, 300c) entre des éléments de forme (310, 320, 330) disjoints de la sous-maille invalide (300a, 300b, 300c), pris deux à deux, - une étape d’identification (422) de deux éléments de forme (310, 320, 330) de la sous-maille invalide (300a, 300b, 300c) dits éléments de forme identifiés, qui sont associés à une distance intérieure, dite distance d’identification, qui est inférieure à la valeur de tolérance prédéterminée, - une étape de formation (423) d’une arête virtuelle qui relie les deux éléments de forme identifiés et dont la longueur vaut la distance d’identification, et - une étape d’application d’au moins un opérateur d’effondrement topologique par fusion d’arête (424), qui contracte une arête (320) qui coïncide avec l’arête virtuelle en un sommet d’effondrement, de manière à supprimer la sous-maille invalide (300a, 300b, 300c). Procédé (400) selon la revendication 1, lorsque les deux éléments de forme identifiés comprennent un sommet (310) et une face (330), - l’étape de formation (423) de l’arête virtuelle comprend, - une première étape d’établissement (4231) d’une première arête virtuelle entre le sommet 320 et un point fictif disposé sur la face (330), - une deuxième étape d’établissement (4232) de deuxièmes arêtes virtuelles entre le point fictif et les sommets (310) de la face (330), le procédé (400) comprenant en outre, - une étape d’application d’un opérateur de découpe topologique (425), sur la face (330), le long de la première arête virtuelle et le long des deuxièmes arêtes virtuelles, de manière à subdiviser la sous-maille invalide (300a, 300b, 300c) en trois sous-sous-mailles. Procédé (400) selon la revendication 1, lorsque les deux éléments de forme identifiés comprennent une première arête (320) et une deuxième arête (320), - l’étape de formation (423) de l’arête virtuelle comprend, - une première étape d’établissement (4233) d’une première arête virtuelle entre un premier point fictif de la première arête (320) et un deuxième point fictif de la deuxième arête (320), - une deuxième étape d’établissement (4234) de deuxièmes arêtes virtuelles entre le premier point fictif et les sommets (310) qui ne sont pas portés par la première arête (320), - une troisième étape d’établissement (4235) de troisièmes arêtes virtuelles entre le deuxième point fictif et les sommets (310) qui ne sont pas portés par la deuxième arête (320), et le procédé (400) comprenant en outre, - une première étape d’application d’un opérateur de découpe topologique (426), sur deux faces (330) adjacentes de la première arête (320), le long des deuxièmes arêtes virtuelles, de manière à subdiviser la sous-maille invalide (300a, 300b, 300c) en deux sous-sous-mailles, - une deuxième étape d’application d’un opérateur de découpe topologique (427), sur deux faces (330) adjacentes de la deuxième arête (320), le long des troisièmes arêtes virtuelles, de manière à subdiviser la sous-maille invalide (300a, 300b, 300c) en deux sous-sous-mailles. Procédé (400) selon l’une quelconque des revendications 1 à 3, dans lequel l’étape de détection (410) comprend l’utilisation d’au moins un prédicat géométrique qui implémente des algorithmes de manipulation d’objets géométriques en arithmétique exacte. Procédé (400) selon l’une quelconque des revendications 1 à 4, dans lequel l’étape de calcul (421) comprend le calcul d’une distance intérieure choisie parmi : une distance entre un sommet (310) et une face (330), une distance entre une première arête (320) et une deuxième arête (320) ou une distance entre un premier sommet (310) et un deuxième sommet (310). Procédé (400) selon l’une quelconque des revendications 1 à 5, dans lequel lorsqu’aucune distance intérieure n’est inférieure à la valeur de tolérance prédéterminée, l’étape d’identification (422) comprend en outre l’identification de deux éléments de forme (310, 320, 330), dits éléments de forme identifiés, qui sont associés à la distance intérieure la plus courte. Procédé (400) selon l’une quelconque des revendications 1 à 6, dans lequel la valeur de tolérance prédéterminée est inférieure à m. Procédé (400) selon l’une quelconque des revendications 1 à 6, dans lequel la valeur de tolérance prédéterminée est inférieure à m. Procédé (400) selon l’une quelconque des revendications 1 à 8, comprenant en outre une étape de disposition du sommet d’effondrement entre les extrémités de l’arête (320) qui coïncide avec l’arête virtuelle. Procédé selon l’une quelconque des revendications 1 à 8, comprenant en outre une étape de disposition du sommet d’effondrement au niveau d’une des extrémités de l’arête (320) qui coïncide avec l’arête virtuelle.