Darcy Flow (Spatial Analyst)
Zusammenfassung
Berechnet den Restwert für den Grundwasservolumenausgleich und andere Ausgaben für gleichmäßige Fließgeschwindigkeit in einer grundwasserführenden Schicht.
Weitere Informationen zur Funktionsweise von Darcy Flow und Darcy Velocity
Verwendung
-
Darcy Flow und Darcy Velocity unterscheiden sich wie folgt:
- Darcy Flow erzeugt im Gegensatz zu Darcy Velocity ein Ausgabe-Raster für das Volumen.
- Darcy Velocity erzeugt Richtungs- und Betrag-Raster als erforderliche Ausgabe, während Darcy Flow diese Ausgaben optional erzeugt.
-
Alle Eingabe-Raster müssen die gleiche Ausdehnung und Zellengröße haben.
-
Alle Eingabe-Raster müssen Gleitkomma-Raster sein.
-
Die Richtung des Geschwindigkeitsvektors wird in Kompasskoordinaten aufgezeichnet, (Gradzahl im Uhrzeigersinn von Norden aus), und der Betrag wird in Längeneinheiten über die Zeit angegeben.
-
Durch dieses Werkzeug wird kein bestimmtes System von Einheiten angegeben. Die Daten müssen konsistent sein und dieselbe Einheit für Zeit (Sekunden, Tage, Jahre) und Länge (Fuß, Meter) verwenden.
-
Das Haupthöhen-Raster kann aus unterschiedlichen Quellen stammen. Es kann mithilfe der Werkzeuge für die Oberflächeninterpolation, z. B. Kriging oder Spline, aus Beobachtungsbrunnendaten interpoliert werden. Die Haupthöhenwerte können aber auch anhand der Ergebnisse eines separaten Modellierungsprogramms ermittelt werden.
Egal, wie das Haupthöhen-Raster erzeugt wird, die Haupthöhe muss dem Transmissivitäts-Raster entsprechen, d. h. sie muss den Fluss durch das Transmissivitätsfeld widerspiegeln. Es reicht nicht aus, durch Messungen und Feldversuche ermittelte Werte zu verwenden. Die gerasterten Werte müssen mithilfe eines geeigneten Programms für den Fluss poröser Medien analysiert und auf Konsistenz überprüft werden. Konsistenz bedeutet, dass die Haupthöhen durch das modellierte Transmissivitätsfeld erzeugt werden. Da das tatsächliche und das modellierte Transmissivitätsfeld in der Praxis häufig nicht identisch sind, unterscheiden sich auch das tatsächliche und das modellierte Haupthöhenfeld. Überprüfen Sie die Haupthöhen auf Konsistenz, indem Sie das von Darcy Flow erzeugte Fehler-Raster untersuchen. Der Fehler spiegelt die Konsistenz des Datasets wider. Jede Analyse inkonsistenter Daten mit Darcy Velocity führt zu irrelevanten Ergebnissen.
-
Das effektive Porositätsfeld, eine physische Eigenschaft der grundwasserführenden Schicht, wird im Allgemeinen aus geologischen Daten geschätzt. Es ist definiert als das Volumen des leeren Raumes, der zum Flüssigkeitsfluss beiträgt, dividiert durch das Gesamtvolumen. Die Porosität wird als Zahl zwischen 0 und 1 ausgedrückt, mit typischen Werten um 0,35, und sie ist dimensionslos. Ein Wert von 0,35 für die effektive Porosität bedeutet, dass 35 Prozent des Volumens des porösen Mediums zum Flüssigkeitsfluss beitragen. Die restlichen 65 Prozent, die aus fester Matrix und unverbundenen Poren bestehen, tragen nicht zum Flüssigkeitsfluss bei.
-
Die in Längeneinheiten gemessene gesättigte Stärke wird aus geologischen Informationen interpretiert. Für eine beschränkte grundwasserführende Schicht ist dieser Messwert die Stärke der Formation zwischen der oberen und unteren Begrenzungsschicht. Für eine nicht beschränkte grundwasserführende Schicht ist die gesättigte Stärke die Entfernung zwischen dem Wasserspiegel und der unteren Begrenzungsschicht.
Syntax
| Parameter | Erläuterung | Datentyp |
in_head_raster |
Das Eingabe-Raster, in dem jeder Zellenwert die Grundwasserhaupthöhe an dieser Position darstellt. Die Haupthöhe ist in der Regel eine Erhebung über einem bestimmten Datum, z. B. dem mittleren Meeresspiegel. | Raster Layer |
in_porosity_raster |
Das Eingabe-Raster, in dem jeder Zellenwert die effektive Porosität der Formation an dieser Position darstellt. | Raster Layer |
in_thickness_raster |
Das Eingabe-Raster, in dem jeder Zellenwert die gesättigte Stärke an dieser Position darstellt. Der Wert für die Stärke wird anhand geologischer Eigenschaften der grundwasserführenden Schicht interpretiert. | Raster Layer |
in_transmissivity_raster | Das Eingabe-Raster, in dem jeder Zellenwert die Transmissivität der Formation an dieser Position darstellt. Die Transmissivität einer grundwasserführenden Schicht ist als die hydraulische Leitfähigkeit K multipliziert mit der Stärke der gesättigten wasserführenden Schicht definiert. Die Einheit ist eine Längeneinheit zum Quadrat über die Zeit. Diese Eigenschaft wird im Allgemeinen anhand experimentelle Felddaten geschätzt, z. B. aus Pumptests. Die Tabellen 1 und 2 in Funktionsweise von "Darcy Flow" und "Darcy Velocity" führen Bereiche hydraulischer Leitfähigkeit für einige generalisierte geologische Materialien auf. | Raster Layer |
out_direction_raster (optional) | Das Ausgabe-Raster für die Fließrichtung. Jeder Zellenwert gibt die Richtung des Sickergeschwindigkeitsvektors (durchschnittliche lineare Geschwindigkeit) im Mittelpunkt der Zelle an, der als Durchschnitt der Sickergeschwindigkeit durch die vier Flächen der Zelle berechnet wird. Er wird zusammen mit dem Ausgabe-Raster für den Betrag verwendet, um den Flussvektor zu beschreiben. | Raster Dataset |
out_magnitude_raster (optional) | Ein optionales Ausgabe-Raster, in dem jeder Zellenwert den Betrag des Sickergeschwindigkeitsvektors (durchschnittliche lineare Geschwindigkeit) im Mittelpunkt der Zelle angibt, der als Durchschnitt der Sickergeschwindigkeit durch die vier Flächen der Zelle berechnet wird. Er wird zusammen mit dem Ausgabe-Raster für die Richtung verwendet, um den Flussvektor zu beschreiben. | Raster Dataset |
Rückgabewert
| Name | Erläuterung | Datentyp |
| out_volume_raster |
Das Ausgabe-Raster für den Volumenausgleich-Restwert. Jeder Zellenwert stellt den Restwert für den Grundwasservolumenausgleich für gleichmäßige Fließgeschwindigkeit in einer grundwasserführenden Schicht dar, wie durch das Darcy-Gesetz bestimmt. | Raster |
Codebeispiel
In diesem Beispiel werden das Grundwasservolumenausgleich-Raster sowie die Fließrichtung und die Sickergeschwindigkeit einer grundwasserführenden Schicht berechnet.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outDarcyFlow = DarcyFlow("gwhead", "gwporo", "gwthick","gwtrans",
"C:/sapyexamples/output/outdarcydir",
"C:/sapyexamples/output/outdarcymag")
outDarcyFlow.save("C:/sapyexamples/output/outdarcyflo")
In diesem Beispiel werden das Grundwasservolumenausgleich-Raster sowie die Fließrichtung und die Sickergeschwindigkeit einer grundwasserführenden Schicht berechnet.
# Name: DarcyFlow_Ex_02.py
# Description: Calculates the groundwater volume balance residual and other
# outputs 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"
outDirectionRaster = "C:/sapyexamples/output/outdarcydir"
outMagnitudeRaster = "C:/sapyexamples/output/outdarcymag"
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute DarcyFlow
outDarcyFlow = DarcyFlow(inHeadRaster, inPorosityRaster, inThicknessRaster,
inTransmissivityRaster, outDirectionRaster,
outMagnitudeRaster)
# Save the output
outDarcyFlow.save("C:/sapyexamples/output/outdarcyflow")