Le Nobel
Connexion
  • > Accueil >
  • InSAR MOI >
  • Thèse
Thèse
  • Sommaire
  • Contexte de l'étude
  • Chapitre 1
  • Chapitre 2
  • Chapitre 3
  • Chapitre 4
  • Chapitre 5
  • Chapitre 6
  • Chapitre 7
  • Chapitre 8
  • Chapitre 9
  • Chapitre 10
Toolbox MOI Matlab
Toolbox MOI Python
next up previous contents
Next: Complexification de la géométrie Up: Modélisation des déplacements post-éruptifs Previous: Modélisation des déplacements post-éruptifs   Contents

Déplacements dus à une source rectangulaire

Nous cherchons à vérifier dans un premier temps si une faille simple peut expliquer les déplacements. La faille a une forme rectangulaire, un pendage et un azimut. Elle subit exclusivement une contrainte cisaillante. La contrainte appliquée est donc coplanaire au plan de la faille mais peut avoir un azimut différent de celui de la faille (Fig. 9.3). Dans un second temps, nous chercherons à vérifier si les effets d'une ouverture sans cisaillement peuvent expliquer les déplacements observés. Nous imposerons alors au vecteur contrainte d'être perpendiculaire au plan de faille. Enfin, nous permettrons à la faille de subir un cisaillement et une ouverture.

Figure 9.3: Géométrie d'un rectangle à 9 ou 10 paramètres. On définit la position $xc$, $yc$ et $zc$ du rectangle, la longueur $L$ du petit côté, et le ratio entre le petit et le grand côté. Le rectangle appartient à un plan d'azimut et de pendage fixé ($az$ et $pd$ plan), et subit une rotation d'un angle $ang$ sur ce plan. Deux types de contraintes sont appliqués au rectangle. Pour modéliser une ouverture du rectangle, le vecteur contrainte qui lui est appliqué est orienté perpendiculairement au plan. Pour modéliser un cisaillement, le vecteur contrainte est coplanaire au rectangle et est orienté ($az$ contrainte). Si le vecteur contrainte est défini par un azimut et un pendage, alors le rectangle subira un cisaillement couplé à une ouverture.

Le tableau 9.3 donne les paramètres des meilleurs modèles trouvés pour chacune des configurations. La figure 9.4 montre les résultats pour le modèle prenant en compte un cisaillement et une ouverture, ou seulement un cisaillement, Les figures pages [*] à [*] donnent les PPD-1D et 2D des l'inversions. Dans le cas d'un vecteur contrainte coplanaire comme pour un vecteur contrainte perpendiculaire au rectangle, la source trouvée est profonde et pend vers l'est. Les déplacements modélisés ne permettent pas de reproduire convenablement les observations.

Figure 9.4: Résultats d'inversion pour un rectangle subissant soit un cisaillement (vecteur contrainte coplanaire), soit une ouverture, soit un cisaillement et une ouverture (vecteur contrainte 3D).

Le modèle prenant en compte un cisaillement et une ouverture permet de mieux reproduire les déplacements. Le vecteur contrainte pend de plus de 70 degrés vers l'est alors que la source ne pend que de 14 degrés. La composante cisaillante de la contrainte est donc plus faible que sa composante en ouverture. L'ouverture moyenne est négative, le plan de glissement perd donc du volume. Enfin, la source pend principalement vers le nord, la partie nord du rectangle est donc plus loin de la surface que la partie sud. Cette géométrie permet de reproduire au premier ordre la forme asymétrique des déplacements.


Table 9.3: Résultat d'inversion pour un rectangle. xc, yc et zc sont les coordonnées du centre, L la longueur du petit côté, rL le ratio entre les longueurs des côtés, ang est l'angle de rotation du rectangle sur le plan, az p, pd p, az t et pd t sont respectivement les azimuts et pendages du rectangle et du vecteur contrainte. Le nord correspond à un azimut de 0°, l'est de -90°.
modèle xc (m) yc (m) zc (m) L (m) rL ang(°) az p(°)
trac. coplan. $ 371109 $ $ 7649766 $ $ -2487 $ $ 2321 $ $ 1.18 $ $ 3.2 $ $ -96 $
trac. perp. $ 370932 $ $ 7649191 $ $ -1266 $ $ 2704 $ $ 1.24 $ $ 39 $ $ -120 $
trac. 3D $ 371095 $ $ 7649197 $ $ 96 $ $ 2997 $ $ 1.7 $ $ 4.0$ $ -10$
 
modèle pd p(°) az t(°) pd t(°) $\Delta V (m^3)$ ouverture moyenne (m) cisaillement moyen (m) misfit %
trac. coplan. $ 40 $ $ -46 $ - -2.52 $\times10^6$ -0.36 5.12 21.8
trac. perp. $ 10 $ - - -8.70 $\times10^6$ -0.97 0.006 24.2
trac. 3D $ 14.2 $ $ -90 $ $ 69.10 $ -7.30 $\times10^6$ -0.48 0.25 16.0



next up previous contents
Next: Complexification de la géométrie Up: Modélisation des déplacements post-éruptifs Previous: Modélisation des déplacements post-éruptifs   Contents
kunos 2014-07-01