ArcMap

Vitesse de Darcy

Disponible avec une licence Spatial Analyst.

Résumé

Calcule le vecteur de vitesse de la nappe phréatique d'infiltration (direction et magnitude) pour un flux constant dans un aquifère.

Pour en savoir plus sur les outils Flux de Darcy et Vitesse de Darcy

Utilisation

  • Différences entre les outils Flux de Darcy et Vitesse de Darcy :

    • L'outil Flux de Darcy produit un raster de volume en sortie, alors que l'outil Vitesse de Darcy ne le fait pas.
    • L'outil Vitesse de Darcy génère des rasters de direction et de grandeur comme sortie requise uniquement, alors que l'outil Flux de Darcy génère ces sorties de manière optionnelle.

  • Tous les rasters en entrée doivent avoir les mêmes étendue et taille de cellule.

  • Tous les rasters en entrée doivent être à virgule flottante.

  • La direction du vecteur de vitesse est enregistrée dans les coordonnées du compas (en degrés dans le sens horaire à partir du Nord) et la grandeur dans les unités de longueur dans le temps.

  • Aucun système d'unités particulier n'est spécifié par cet outil. Les données doivent être cohérentes et toutes les données doivent utiliser la même unité de temps (secondes, jours, années) et de longueur (pieds, mètres).

  • Le raster d'altitude peut provenir de diverses sources. Il peut être interpolé à partir des données d'un puits d'observation à l'aide d'un des outils d'interpolation de surface, tel que Krigeage ou Spline. Les valeurs de hauteur peuvent également être obtenues à partir des résultats d'un programme de modélisation distinct.

    Quelle que soit la manière dont le raster d'altitude est obtenu, la hauteur doit être cohérente avec le raster de transmissivité. Autrement dit, la hauteur doit refléter le flux à travers le champ de transmissivité. Il ne suffit pas d'utiliser des valeurs obtenues par un mesurage et des tests sur le terrain. La cohérence des valeurs rastérisées doit être analysée à l'aide d'un programme de flux en milieu poreux approprié. La cohérence implique que les hauteurs soient produites réellement par le champ de transmissivité modélisé. Comme les champs de transmissivité véritable et modélisé diffèrent souvent en pratique, les champs de hauteur véritable et modélisé diffèrent également. Vérifiez la cohérence des hauteurs en examinant le raster résiduel généré par l'outil Flux de Darcy. Le raster résiduel reflète la cohérence du jeu de données. Toute analyse utilisant l'outil Vitesse de Darcy sur des jeux de données incohérents produira des résultats dépourvus de signification.

  • Le champ de porosité effective, une propriété physique de l'aquifère, est estimé en général à partir de données géologiques. Il est défini comme le volume d'espace vide qui contribue au flux de liquide divisé par le volume entier. La porosité est exprimée sous la forme d'un nombre compris entre 0 et 1, avec des valeurs standard avoisinant 0,35, et est une valeur sans dimension. Une valeur de porosité effective de 0,35 signifie que 35% du volume du milieu poreux contribuent au flux de liquide. Les 65% restants, composés d'une matrice solide et de pores non connectés, ne contribue pas au flux de liquide.

  • L'épaisseur saturée, mesurée en unités de longueur, est interprétée à partir d'informations géologiques. Pour un aquifère confiné, cette mesure correspond à l'épaisseur de la formation entre les couches de confinement supérieure et inférieure. Pour un aquifère non confiné, l'épaisseur saturée correspond à la distance entre la nappe phréatique et la couche de confinement inférieure.

  • Les rasters en sortie sont de type virgule flottante.

  • Pour plus d’informations sur les environnements de géotraitement qui s’appliquent à cet outil, reportez-vous à la rubrique Environnements d’analyse et Spatial Analyst.

Syntaxe

DarcyVelocity(in_head_raster, in_porosity_raster, in_thickness_raster, in_transmissivity_raster, out_magnitude_raster)
ParamètreExplicationType de données
in_head_raster

Raster en entrée dans lequel chaque valeur de cellule représente l'altitude de nappe à cet emplacement.

La valeur représente généralement une altitude au-dessus d'un autre datum tel que le niveau moyen de la mer.

Raster Layer
in_porosity_raster

Raster en entrée dans lequel chaque valeur de cellule représente la porosité de formation effective à cet emplacement.

Raster Layer
in_thickness_raster

Raster en entrée dans lequel chaque valeur de cellule représente l'épaisseur saturée à cet emplacement.

La valeur de l'épaisseur est interprétée à partir des propriétés géologiques de l'aquifère.

Raster Layer
in_transmissivity_raster

Raster en entrée dans lequel chaque valeur de cellule représente la transmissivité de formation à cet emplacement.

La transmissivité d'un aquifère est définie comme la conductivité hydraulique K fois l'épaisseur saturée de l'aquifère b, en unités de longueur au carré dans le temps. Cette propriété est généralement évaluée à partir de données expérimentales comme des tests de pompage. Les tableaux 1 et 2 de la section sur le fonctionnement du flux et de la vitesse de Darcy répertorient les plages de conductivité hydraulique de certains matériaux géologiques généralisés.

Raster Layer
out_magnitude_raster

Raster de direction de flux en sortie.

Chaque valeur de cellule représente la direction du vecteur de vitesse d'infiltration (vitesse linéaire moyenne) au centre de la cellule, calculée en tant que pourcentage de la vitesse d'infiltration à travers les quatre faces de la cellule.

Il est utilisé avec le raster de magnitude en sortie pour décrire le vecteur de flux.

Raster Dataset

Valeur renvoyée

NomExplicationType de données
out_direction_raster

Raster de direction de flux en sortie.

Chaque valeur de cellule représente la direction du vecteur de vitesse d'infiltration (vitesse linéaire moyenne) au centre de la cellule, calculée en tant que pourcentage de la vitesse d'infiltration à travers les quatre faces de la cellule.

Il est utilisé avec le raster de magnitude en sortie pour décrire le vecteur de flux.

Raster

Exemple de code

Premier exemple d'utilisation de l'outil DarcyVelocity (fenêtre Python)

Calcule la vitesse de la nappe phréatique d'infiltration (direction et magnitude) pour un flux constant dans un aquifère.

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outDarcyVelocity = DarcyVelocity("gwhead", "gwporo", "gwthick", "gwtrans", 
                            "C:/sapyexamples/output/outdarcymag")
outDarcyVelocity.save("c:/sapyexamples/output/outdarcyvel")
Deuxième exemple d'utilisation de l'outil DarcyVelocity (script autonome)

Calcule la vitesse de la nappe phréatique d'infiltration (direction et magnitude) pour un flux constant dans un aquifère.

# Name: DarcyVelocity_Ex_02.py
# Description: Calculates the groundwater seepage velocity 
#              vector (direction and magnitude) for steady 
#              flow in an aquifer.
# Requirements: Spatial Analyst Extension

# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *

# Set environment settings
env.workspace = "C:/sapyexamples/data"

# Set local variables
inHeadRaster = "gwhead"
inPorosityRaster = "gwporo"
inThicknessRaster = "gwthick"
inTransmissivityRaster = "gwtrans"
outMagnitudeRaster = "C:/sapyexamples/output/outdarcymag"

# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")

# Execute DarcyVelocity
outDarcyVelocity = DarcyVelocity(inHeadRaster, inPorosityRaster, inThicknessRaster,
                            inTransmissivityRaster, outMagnitudeRaster)

# Save the output 
outDarcyVelocity.save("C:/sapyexamples/output/outdarcyvel")

Environnements

Informations de licence

  • Basic: Requiert Spatial Analyst
  • Standard: Requiert Spatial Analyst
  • Advanced: Requiert Spatial Analyst

Rubriques connexes