孔隙扩散 (空间分析)

摘要

计算与时间相关的二维浓度分布,形式为在某一离散点瞬时注入垂直混合蓄水层的单位体积溶质质量。

了解有关“孔隙扩散”工作原理的详细信息

用法

语法

PorousPuff (in_track_file, in_porosity_raster, in_thickness_raster, mass, {dispersion_time}, {longitudinal_dispersivity}, {dispersivity_ratio}, {retardation_factor}, {decay_coefficient})
参数说明数据类型
in_track_file

输入粒子追踪路径文件。

这是一个 ASCII 文本文件,其中包含沿路径移动的位置、局部速度矢量以及累积长度和时间的信息。

此文件将使用粒子追踪工具进行生成。

File
in_porosity_raster

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

Raster Layer
in_thickness_raster

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

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

Raster Layer
mass

源点处瞬间释放的质量数量值(以质量单位为单位)。

Double
dispersion_time
(可选)

表示溶质扩散的时间范围的值(以时间单位为单位)。

该时间必须小于或等于追踪文件中的最大时间。如果请求的时间超过追踪文件中的可用时间,该工具将中止。默认时间为追踪文件中的最晚时间(对应于终点)。

Double
longitudinal_dispersivity
(可选)

表示与流向平行的扩散性的值。

有关默认值的确定方法及其与研究范围的关系的详细信息,请参阅本文档的孔隙扩散的工作原理部分。

Double
dispersivity_ratio
(可选)

表示纵横扩散性比的值。

横向扩散性在同一水平面中垂直于流向。默认值为三。

Double
retardation_factor
(可选)

表示蓄水层中溶质延迟的无维度值。

延迟可以是一至无穷大之间的任意值,一对应于无延迟。默认值为一。

Double
decay_coefficient
(可选)

经历一阶指数衰减的溶质(如放射性核素)的衰减系数,以反时间单位为单位。

默认值为零,对应于无衰减。

Double

返回值

名称说明数据类型
out_raster

浓度分布的输出栅格。

每个像元值都表示该位置的浓度。

Raster

代码示例

PorousPuff 示例 1(Python 窗口)

此例对所需输入执行该工具,然后输出浓度分布栅格。

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outPorousPuff = PorousPuff("trackfile.txt", "gwporo", "gwthick", 50, 10000, "", 3,
                         "", "")
outPorousPuff.save("c:/sapyexamples/output/outporpuff")
PorousPuff 示例 2(独立脚本)

此例对所需输入执行该工具,然后输出浓度分布栅格。

# Name: PorousPuff_Ex_02.py
# Description: Calculates the time-dependent, two-dimensional 
#              concentration distribution in mass per volume of a 
#              solute introduced instantaneously and at a discrete 
#              point into a vertically mixed 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
inTrackFile = "trackfile.txt"
inPorosityRaster = "gwporo"
inThicknessRaster = "gwthick"
mass = 50
dispersionTime = 10000
longitudinalDispersivity = ""
dispersivityRatio = 3 
retardationFactor = "" 
decayCoefficient = 0


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

# Execute PorousPuff
outPorousPuff = PorousPuff(inTrackFile, inPorosityRaster, inThicknessRaster, 
                        mass, dispersionTime, longitudinalDispersivity,
                        dispersivityRatio, retardationFactor, 
                        decayCoefficient)

# Save the output 
outPorousPuff.save("c:/sapyexamples/output/outporpuff")

环境

相关主题

许可信息

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

7/10/2012