✍️ 陈箫翰 , 何星辰  •  📅 2026-01-03

grdvs30

官方文档:

grdvs30

简介:

计算 Vs30 速度模型

Vs30,即上部30米内的时间平均剪切波速度,是地震工程领域用来评估地震场地条件的关键指标。 grdvs30 模块通过地形数据计算某个范围的 Vs30 速度模型。

输入一个地理坐标系的地形网格文件,以及一个常量克拉通值或一个克拉通文件作为权重, 范围从0(对应稳定的地盾区,即克拉通)到1(对应活动构造区)。

一个可选参数 water_vel 表示水覆盖区域的速度值,默认值为600。

语法

gmt grdvs30 ingrid -Goutgrid -Cval|fname[+g] [ -Rregion ] [ -V[level] ] [ -Wwater_vel ] [ -iflags ] [ -rreg ] [ -:[i|o] ]

输入数据

ingrid[=ID|?varname][+bband][+ddivisor][+ninvalid][+ooffset][+sscale]

输入网格名。通过追加 =ID 可指定 网格格式 [默认为 =nf]。 追加 ?varname 可指定 NetCDF 变量 [默认为 GMT 找到的第一个 2-D 网格]。 参数详细介绍请参考 读 netCDF 文件

必须选项

-C
-Cval|fname[+g]

参数可以是以下三种之一:

  • 一个介于 0 和 1 之间的值 val ,其中 0 表示稳定的克拉通区,1 表示活动构造区。

  • 一个克拉通区域多边形的多段文件的名称 fname

  • 一个包含克拉通/活动构造区域的网格文件的名称 fname 。 在这种情况下,必须使用 +g 后缀来表明正在读取的是网格文件。

-G
-Goutgrid

输出的 Vs30 速度模型网格文件名

可选选项

-R
-Rxmin/xmax/ymin/ymax[+r][+uunit]

指定数据范围。 (参数详细介绍)

-W
-Wwater_vel

设置水域区域的 Vs30 值 [默认值=600]

-V
-V[level]

设置 verbose 等级 [w]。 (参数详细介绍)

-i
-icols[+l][+sscale][+ooffset][,...][,t[word]]

对输入的数据进行列选择以及简单的代数运算。 (参数详细介绍)

-r
-r[g|p]

设置网格配置方式 [默认为网格线配准]。 (参数详细介绍)

-:
-:[i|o]

交换输入或输出数据的前两列。 (参数详细介绍)

-^-

显示简短的帮助信息,包括模块简介和基本语法信息(Windows下只能使用 -

-++

显示帮助信息,包括模块简介、基本语法以及模块特有选项的说明

-? 或无参数

显示完整的帮助信息,包括模块简介、基本语法以及所有选项的说明

--PAR=value

临时修改GMT参数的值,可重复多次使用。参数列表见 配置参数

示例

计算 topo.grd 网格的 Vs30 估计值,并将克拉通值设为0:

gmt grdvs30 topo.grd -C0 -Gvs30.grd

参考资料

相关模块

grdshake