Procédé pour quantifier les déplacements dans un modèle analogique représentant l’évolution de la déformation d’une formation souterraine au cours des temps géologiques, à partir d’informations géologiques et mécaniques relatives à chacune des couches géologiques de ladite formation desquelles on déduit au moins un ordre de dépôt desdites couches géologiques, des épaisseurs desdites couches géologiques, un comportement mécanique desdites couches géologiques, et des conditions aux limites subies par ladite formation, au moyen d’une boite de déformation pour construire un modèle analogique de ladite formation, ladite boite de déformation pouvant ou non comporter au moins une paroi mobile, une étape de mise à l’échelle étant préalablement appliquée dans laquelle on définit des épaisseurs des couches dudit modèle analogique représentatives desdites épaisseurs desdites couches géologiques, des conditions aux limites à l’échelle dudit modèle analogique représentatives desdites conditions aux limites subies par ladite formation, et des matériaux des couches dudit modèle analogique représentatifs dudit comportement mécanique desdites couches géologiques, ledit procédé étant caractérisé en ce qu’on applique au moins les étapes suivantes : A) on construit ledit modèle analogique à un temps initial de la manière suivante : on superpose, dans ladite boite de déformation, selon ledit ordre de dépôt et selon lesdites épaisseurs desdites couches dudit modèle analogique, une succession de couches dudit modèle analogique formées dans lesdits matériaux desdites couches dudit modèle analogique, au moins une desdites couches dudit modèle analogique comprenant une pluralité de marqueurs détectables par des rayons X ; B) au moyen de ladite boite de déformation, on applique lesdites conditions aux limites à l’échelle dudit modèle analogique audit modèle analogique audit temps initial ainsi construit, et on acquière simultanément une succession d’images tridimensionnelles entre ledit temps initial et un temps final pour une succession de pas de temps, au moyen d’une imagerie à rayons X, de préférence au moyen d’un scanner à rayons X couplé à une tomographie ; C) Pour chaque couple d’images tridimensionnelles consécutives dans ladite succession, on applique une méthode de flux optique, et on détermine un champ de déplacement pour chacun desdits couples desdites images tridimensionnelles. Procédé selon la revendication 1, dans lequel, pour chacune desdites couches dudit modèle analogique, si ladite couche géologique correspondante a un comportement mécanique ductile, on mélange au préalable ledit matériau de ladite couche dudit modèle analogique à ladite pluralité de marqueurs, et on dépose ledit matériau dans ladite boite de déformation. Procédé selon l’une des revendications précédentes, dans lequel, pour chacune desdites couches dudit modèle analogique, si ladite couche géologique correspondante a un comportement mécanique fragile, on verse ledit matériau de ladite couche dudit modèle analogique dans ladite boite de déformation, et on insère, au fur et à mesure dudit versement dudit matériau, ladite pluralité de marqueurs. Procédé selon l’une des revendications précédentes, dans lequel, pour chacune desdites couches dudit modèle analogique, si ladite couche géologique correspondante a un comportement mécanique fragile, on utilise du sable, de la poudre de verre ou des microbilles de verre comme matériau de ladite couche dudit modèle analogique. Procédé selon la revendication 4, si on utilise du sable comme matériau de ladite couche dudit modèle analogique, on utilise des marqueurs formés par des particules de verre expansé, de préférence de taille comprise entre 500 µm et 1000 µm. Procédé selon l’une des revendications précédentes, dans lequel, pour chacune desdites couches dudit modèle analogique, si la ladite couche géologique correspondante a un comportement mécanique ductile, on utilise de la silicone comme matériau de ladite couche dudit modèle analogique. Procédé selon la revendication 6, on utilise des marqueurs formés par des particules de verre expansé enrobées d’iodure de potassium. Procédé selon l’une des revendications précédentes, dans lequel ladite méthode de flux optique est choisie parmi la famille des méthodes différentielles d'estimation du flux optique, ladite méthode de flux optique de la famille des méthodes différentielles d'estimation du flux optique pouvant être une méthode globale ou une méthode locale. Procédé selon la revendication 8, dans lequel ladite méthode de flux optique de la famille des méthodes différentielles d'estimation du flux optique est une méthode locale, et dans lequel, on utilise une méthode de détection de coins pour sélectionner des points desdites images tridimensionnelles sur lesquels ladite méthode de flux optique de la famille des méthodes différentielles locale est appliquée, on obtient des vecteurs de déplacement auxdits points sélectionnés pour une desdites images tridimensionnelles de chacun desdits couples desdites images tridimensionnelles, et on interpole lesdits vecteurs de déplacement obtenus auxdits points sélectionnés pour une desdites images tridimensionnelles de chacun desdits couples desdites images tridimensionnelles afin de déterminer ledit champ de déplacement pour chacun desdits couples desdites images tridimensionnelles. Procédé selon l’une des revendications précédentes, dans lequel, à partir au moins dudit champ de déplacement déterminé pour chacun desdits couples desdites images tridimensionnelles, on construit un modèle géologique numérique de ladite formation, on réalise au moins une simulation géomécanique et/ou une simulation d’écoulement dans ladite formation au moyen dudit modèle géologique numérique, et à partir au moins de résultats de ladite au moins une simulation, on fore au moins un puits dans ladite formation pour exploiter une ressource de ladite formation ou pour stocker un fluide ou un solide dans ladite formation.