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

# grdvs30

官方文档:

  grdvs30

简介:

  计算 Vs30 速度模型

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

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

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

## 语法

gmt grdvs30 ingrid `-G` outgrid `-C` val|fname[+g] [ `-R` region ] [ `-V`[level] ] [ `-W` water_vel ] [ `-i` flags ] [ `-r` reg ] [ `-:`[i|o] ]

## 输入数据

ingrid[=ID|?varname][+b band][+d divisor][+n invalid][+o offset][+s scale]

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

## 必须选项

-C

-C val|fname[+g]

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

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

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

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

-G

-G outgrid

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

## 可选选项

-R

-R xmin/xmax/ymin/ymax[+r][+u unit]

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

-W

-W water_vel

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

-V

-V[level]

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

-i

-i cols[+l][+s scale][+o offset][,...][,t[word]]

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

-r

-r[g|p]

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

-:

-:[i|o]

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

-^ 或 -

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

-+ 或 +

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

-? 或无参数

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

--PAR=value

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

## 示例

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

    gmt grdvs30 topo.grd -C0 -Gvs30.grd

## 参考资料

- https://github.com/usgs/earthquake-global_vs30/blob/master/src/grad2vs30.c

## 相关模块

grdshake