达西流 (空间分析)

摘要

计算蓄水层中稳流的地下水量平衡残差以及其他输出。

了解有关“达西流”和“达西速度”工作原理的详细信息

用法

语法

DarcyFlow (in_head_raster, in_porosity_raster, in_thickness_raster, in_transmissivity_raster, {out_direction_raster}, {out_magnitude_raster})
参数说明数据类型
in_head_raster

所包含的每一单元值都代表该处地下水位高程的输入栅格。

水位通常是高于某些基准面(如平均海平面)的高程。

Raster Layer
in_porosity_raster

所包含的每一单元值都代表该处有效地层孔隙度的输入栅格。

Raster Layer
in_thickness_raster

所包含的每一单元值都代表该处饱和厚度的输入栅格。

厚度值根据蓄水层的地质属性进行解释。

Raster Layer
in_transmissivity_raster

所包含的每一单元值都代表该地层导水系数的输入栅格。

蓄水层的导水系数定义为导水率 K 乘以饱和蓄水层厚度 b,长度单位随时间变化进行乘方。此属性通常从字段实验数据(例如抽水测试)中估计得出。达西流和达西速度的工作原理中的表 1 和 2 列出了一些常见地质材料的导水率范围。

Raster Layer
out_direction_raster
(可选)

输出流向栅格。

每一单元值都表示单元中心渗流速度矢量(平均线速度)的方向,以通过单元四个面渗流速度平均值的形式进行计算。

该栅格与输出量级栅格一起来描述流向矢量。

Raster Dataset
out_magnitude_raster
(可选)

可选输出栅格,其中的每一单元值都表示单元中心渗流速度矢量(平均线速度)的量级,以通过单元四个面渗流速度平均值的形式进行计算。

该栅格与输出方向栅格一起来描述流向矢量。

Raster Dataset

返回值

名称说明数据类型
out_volume_raster

输出水量平衡残差栅格数据。

其中各个像元均用于表示蓄水层中稳流的地下水量平衡残差,具体取决于达西定律。

Raster

代码示例

DarcyFlow 示例 1(Python 窗口)

本例将计算蓄水层的地下水量平衡栅格以及流向和渗流速度。

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")
DarcyFlow 示例 2(独立脚本)

本例将计算蓄水层的地下水量平衡栅格以及流向和渗流速度。

# 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")

环境

相关主题

许可信息

ArcView: 需要 Spatial Analyst
ArcEditor: 需要 Spatial Analyst
ArcInfo: 需要 Spatial Analyst

7/10/2012