Cartes de Variogrammes
Généralités
Un aspect important dans la réalisation d’évaluations géostatistiques consiste à comprendre
l’anisotropie des données, ou quelle direction a la plus importante continuité et comment les
valeurs des données changent en fonction de cette direction, mais aussi en fonction des deux
directions mutuellement perpendiculaires.
Une carte de variogramme est un outil de Surpac qui vous permet de visualiser l’anisotropie
dans un plan et de calculer les paramètres anisotropes ellipsoïdes afin de les utiliser dans des
estimations.
Vous découvrirez:
• Le variogramme primaire
• Le variogramme secondaire
• Le calcul des paramètres anisotropes ellipsoïdes.
Conditions
Afin de comprendre ces informations vous devez:
• être familiarisé avec les fichiers string de Surpac
• Savoir comment calculer et modéliser un variogramme avec Surpac
• Comprendre le concept de l’ellipsoïde anisotropique
• Comprendre quels sont les paramètres qui définissent une ellipsoïde anisotropique
Le variogramme primaire
Tâche : Créer une carte de variogramme horizontale
1. Choisissez « Geostatistics > variogram modelling » afin d’ouvrir la fenêtre de
modélisation du variogramme.
2. Choisissez « variogram map > New variogram map ».
3. Entrez les informations comme ci-dessous :
Le nombre de variogrammes sélectionnés déterminera l’angle d’incrémentation. Dans notre
exemple, 16 variogrammes engendrent une incrémentation angulaire de 22,5 degrés
(360/16=22.5). Si le nombre de variogramme tend vers 36, Vous obtiendrez un incrément de
10 degrés (360/36=10).
Les paramètres de dispersion et les paramètres limites de dispersion sont les mêmes que ceux
dans un model de variogramme normal.
Vous devez considérer la relation entre l’angle d’incrémentation et l’angle de dispersion.
Pour certains, il peut paraitre déraisonnable de définir une tolérance de dispersion supérieure
à la moitié de la valeur de l’angle d’incrémentation. Pour ces données, à cause du faible
nombre de pairs, Si un angle de dispersion de 7,5 degrés est utilisé (la moitié des 15 degrés de
l’angle d’incrémentation entre les variogrammes adjacents), le nombre de données paires sera
si faible qu’aucun variogramme viable ne pourra être utilisé. Une dispersion de 30 degrés est
utilisée afin d’assurer un échantillonnage suffisant.
Considérant une dispersion de 30 degrés, nous pourrons nous demander si le nombre de
variogrammes ne pourrait pas être réduit afin de minimiser les chevauchement des cônes pour
les variogrammes adjacents. Bien que cet argument soit plausible en théorie, dans la pratique,
utiliser une carte de variogrammes afin de déterminer visuellement l’anisotropie dans un plan
exclus la réduction de l’angle d’incrémentation.
C’est vous, ainsi que les données que vous utilisez, qui déterminerez comment les deux
tâches pourront être complétées simultanément. Voici maintenant un autre exemple qui
prouve que la géostatistique est une science inexacte. Un essai avec des données vous
permettra de savoir quels paramètres sont acceptables.
Au bas de la fenêtre, le décalage, la distance maximum et les paramètres du variogramme
sont spécifiés, exactement de la même manière que dans la modélisation du variogramme.
Détail important : la distance maximum sera le rayon de la carte de variogramme. Vous vous
apercevrez que vous devrez faire varier cette valeur pour trouver celle qui convient le mieux.
Cliquez sur l’onglet « Advanced », entrez les informations comme ci-dessous, puis cliquez
sur « Apply » :
4. Cliquez droit sur le variogramme, et sélectionnez « Tile Windows ».
5. Choisissez « Display > Set axis limits ».
6. Entrez les informations comme ci-dessous, puis cliquez sur « Apply » :
Tous les variogrammes s’affichent :
Tâche : Déterminer l’anisotropie en 2 dimensions
Dans le cas de la 2D, les deux axes, principal et semi-principal, seront représentés dans le
plan de la carte du variogramme primaire. L’orientation de l’axe principal est déterminée par
celle du variogramme qui a la plus longue série pour une valeur de seuil donnée.
1. Utilisez le curseur de décalage afin de vous rendre compte des zones de fortes et de
faibles variances de la carte. En d’autres termes, bougez le curseur d’avant en arrière
et regarder les variations de couleurs du variogramme. Vous pourrez, par ce biais, être
plus à même de visualiser l’éventail de valeurs de décalage, il apparaitra des zones de
la carte élevées de façon consistante et d’autres moins. Dans l’exemple précédent,
l’orientation de 45 degrés fait apparaître de façon consistante des couleurs
représentant de faibles variances.
2. Dès que vous savez ce qui semblera être l’orientation avec la plus longue série, utiliser
le curseur de décalage pour améliorer la qualité du variogramme expérimental pour
cette direction. Dans l’image précédente, un décalage de distance de 4,3 permet
d’obtenir un variogramme satisfaisant avec un angle de 45 degrés.
3. Quand vous obtenez en variogramme expérimental convenable, créez un modèle de
variogramme pour cette orientation.
4. Choisissez « Variogram > Model ».
5. Cliquez-glissez le modèle afin de faire correspondre le variogramme pour cette
direction. Comme ci-dessous, faites correspondre le modèle à une orientation de 45
degrés :
6. Vérifiez ce modèle de variogramme pour toutes les autres orientations. L’axe
principal doit être celui du variogramme qui possède la plus faible variance pour la
plus grande distance. Dans notre cas, c’est le variogramme d’orientation 22,5 degrés
qui possède la gamme la plus importante.
7. Si une autre orientation semble avoir une plus grande gamme et une variance plus
faible, modifier le modèle pour ajuster le variogramme expérimental. Modifier l’écart
de distance pour cette direction peut aider à obtenir un meilleur ajustement. Modifier
le variogramme pour une orientation de 22,5 degrés comme ci-dessous :
8. Répétez les deux étapes précédentes jusqu’à ce que vous soyez satisfait de
l’orientation de l’axe principal que vous obtenez.
9. Après avoir défini l’axe majeur, demandez-vous ainsi qu’à d’autres personnes
familières avec la géologie, si l’orientation semble correcte. Dans notre cas, une
valeur de 22,5 degrés correspond bien à l’orientation du gisement présenté ci-
dessous :
Pour la zone de gisement #1, l’inclinaison de l’axe principal est de 22,5 degrés. Sauvegardez
le variogramme.
10. Choisissez « File > Save > Experimental variogram and model »
11. Entrez les informations comme ci-dessous, puis cliquez sur « Apply » :
Comme vous le constatez, non seulement la modélisation de variogramme peut être
considérée comme non-scientifique mais l’orientation de l’axe principal peut aussi être sujette
à interprétations et débats.
L’axe semi principal est dans le même plan et perpendiculaire à l’axe principal. Par
conséquent, l’inclinaison de l’axe semi principal pour nos données est de 22,5 + 90 = 112,5.
Le ratio d’anisotropie est la gamme de l’axe principal / la gamme de l’axe semi-principal pour
une valeur de seuil constante.
La gamme de l’axe principal est affichée, ici 15,1, dans le coin supérieur droit de la fenêtre de
modélisation du variogramme.
Déterminer la série de l’axe semi-principal pour un variogramme avec une structure unique et
simple, il suffit de changer la gamme afin ajuster le variogramme tout en conservant la même
valeur de seuil pour l’axe semi-principal.
12. Choisissez « Variogram > Model »
13. Cliquez-glissez le variogramme afin d’ajuster le variogramme expérimental pour une
inclinaison de 112,5 degrés.
Remarque : Le seuil du variogramme est exactement le même mais la gamme est
désormais de 6,1.
Par conséquent, le ratio d’anisotropie pour la zone 1 est de : 15.1 / 6.1 = 2.48.
Il faut admettre que ce n’est pas la manière d’utilisation habituelle du logiciel mais c’est la
méthode la plus rapide afin de déterminer l’anisotropie en 2D.