3. Modélisation A#

3.1. Caractéristiques de la modélisation A#

Le repère global est choisi de manière à ce que la discontinuité soit horizontale dans ce repère (\(\alpha =0°\) donc un angle d’Euler de \(90°\) ).

La modélisation est réalisée en déformation plane avec 1092 éléments QU4 pour le massif et 91 éléments QU4 pour la discontinuité.

Discrétisation en temps :

  • 25 pas de temps pour les \(7,5\) premières années

  • 25 pas de temps pour les \(12,5\) années suivantes.

3.2. Grandeurs testées et résultats#

On compare les résultats obtenus avec Code_Aster avec ceux obtenus par le code LAGAMINE .

On présente le profil de pression le long de la discontinuité à \(7,5\) et \(20\) ans (figure ), les variations de l’ouverture le long de la discontinuité à \(7,5\) et \(20\) ans (figure ) et le flux sortant de la fissure au niveau du puits par jour (figure ).

_images/10000000000003980000028038EC468F3E07BBC9.jpg

Illustration 3.2.1: Pression le long de la discontinuité à 7,5 ans et 20 ans et comparaison avec les résultats de LAGAMINE

_images/10000000000003980000027BAC4044696631F394.jpg

Illustration 3.2.2: Ouverture de la discontinuité le long du réservoir après 7,5 et 20 ans et comparaison avec les résultats de LAGAMINE

Les résultats sont quasiment identiques à ceux obtenus avec LAGAMINE .

_images/100000000000037F0000027A70FF5A1E13094534.jpg

Illustration 3.2.3: Débit massique sortant de la fissure au niveau du puits par jour au cours du temps obtenu avec Code_Aster

Les valeurs de références sont celles obtenues avec LAGAMINE .

\(X(m)\)

\(Y(m)\)

Temps (années)

Référence \(\mathrm{PRE1}\) ( \(\mathrm{MPa}\) )

Tolérance ( \(\text{\%}\) )

256,97

0

7,5

39,95

1.0%

516,51

0

7,5

42,68

1.0%

256,97

0

20

38,45

1.0%

516,51

0

20

40,98

1.0%

\(X(m)\)

\(Y(m)\)

Temps (années)

Référence \(\mathrm{V1}\) ( \(m\) )

Tolérance ( \(\text{\%}\) )

516,51

0

7,5

2,20E-04

1.0%

516,51

0

20

2,04E-04

1.0%

On ajoute un cas de non régression :

\(X(m)\)

\(Y(m)\)

Temps (années)

Référence \(\mathrm{FH11X}\) \(({\mathrm{kg.m}}^{-1}.{s}^{-1})\)

Tolérance ( \(\text{\%}\) )

516,51

0

7,5

-7,0189E-03

0.10%

516,51

0

20

-5,3847E-03

0.10%