达西速度 (空间分析)

摘要

计算蓄水层中稳流的地下水渗流速度矢量(方向和模)。

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

用法

语法

DarcyVelocity (in_head_raster, in_porosity_raster, in_thickness_raster, in_transmissivity_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_magnitude_raster

输出流向栅格。

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

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

Raster Dataset

返回值

名称说明数据类型
out_direction_raster

输出流向栅格。

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

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

Raster

代码示例

DarcyVelocity 示例 1(Python 窗口)

计算蓄水层中稳流的地下水渗流速度(方向和模)。

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

计算蓄水层中稳流的地下水渗流速度(方向和模)。

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

环境

相关主题

许可信息

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

7/10/2012