ダルシー流速(Darcy Velocity) (Spatial Analyst)

サマリ

帯水層の安定した流れにおける地下水流の浸透速度のベクトル(方向と大きさ)を計算します。

[ダルシー フロー(Darcy Flow)] ツールと [ダルシー流速(Darcy Velocity)] ツールの仕組みの詳細

使用法

構文

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 を掛け合わせたもので、単位は長さの二乗を時間で割ったものと定義されます。この特質は一般に揚水試験のような現地実験データから推測されます。「Darcy Flow および Darcy Velocity の仕組み」の表 1 と表 2 に地質素材の透水係数の範囲が示されています。

Raster Layer
out_magnitude_raster

出力方向ラスタ。

各セル値は、セルの中心部の浸透速度ベクトル(平均線形速度)の方向を表します。セルの 4 辺を通過する浸透速度の平均値として計算されます。

これはフロー ベクトルを記述する出力マグニチュード ラスタと共に使用されます。

Raster Dataset

リターン

名前説明データ タイプ
out_direction_raster

出力方向ラスタ。

各セル値は、セルの中心部の浸透速度ベクトル(平均線形速度)の方向を表します。セルの 4 辺を通過する浸透速度の平均値として計算されます。

これはフロー ベクトルを記述する出力マグニチュード ラスタと共に使用されます。

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