✍️ 朱邓达  •  🗓️ 2026-08-20 (创建时间)

static_sproj

简介:

根据接收断层形态将静态应力张量投影为法向应力变化和剪应力变化

语法

grt static sproj -Gingrid [ -Mgeometry | -Qfile ] [ -h ]

描述

static_sproj 模块读取 static_stress 模块生成的六个应力张量分量, 根据接收断层的走向、倾角和滑动角,计算每个接收点的法向应力变化 \(\Delta \sigma_\text{n}\) (sigma_n) 和沿 rake 方向的剪应力变化 \(\Delta \tau_\text{s}\) (tau_s), 并将这两个变量写回输入文件。如果变量已经存在,模块会覆盖其数据并给出警告;其他已有变量和属性保持不变。

sigma_n 采用张开为正的约定, tau_s 是牵引力在接收断层滑动方向上的有符号投影, 与 rake 方向同向为正。结果单位与输入应力相同,即 dyne/cm²(0.1 Pa)。

模块根据输入文件的 layout 属性区分 grid、points 和 faults 三种布局, 并使用对应的接收断层形态来源。

  • grid 布局没有接收断层形态,必须设置 -Mstrike/dip/rake。

  • 普通 points 布局可以在文件中包含 strike(point)、dip(point) 和 rake(point)。 如果这些变量完整且每个点的值都已定义,则不设置 -M 时使用文件中的形态。 若设置 -M,则忽略文件中的形态并给出警告。 如果文件没有完整的接收断层形态,则必须设置 -M。

  • 普通 points 布局还可以使用 -Qfile 为每个点传入新的形态。 -Q 文件每行必须包含六列 north east depth strike dip rake,且点数、坐标和顺序必须与输入文件一致。 -Q 与 -M 不能同时设置。

  • faults 布局表示有限接收断层,包含 point 和 nfault 维度。 此时输入文件中已经记录每个有限断层的 strike 和 dip,但 rake 可能未定义。 若任意断层的 rake 未定义,必须设置 -Mrake,此时未定义的 rake 使用手动值设置的值, 已定义的 rake 仍使用文件值。若所有 rake 都已定义,设置 -Mrake 会报错。 使用 -Mrake+f 可以强制所有点使用手动 rake。

必选选项

-Gingrid

输入静态合成文件。文件必须包含 static_stress 模块写入的六个 stress_?? 变量,模块会将投影结果写回该文件。

可选选项

-Mgeometry

手动指定接收断层形态。 对 grid 和普通 points 布局,格式为 strike/dip/rake,单位为度。 对有限接收断层 faults 布局,格式为 rake 或 rake+f,单位为度。 +f 表示强制所有点使用该 rake。

-Qfile

为普通 points 布局指定逐点接收断层形态。 每行格式为 north east depth strike dip rake,单位分别为 km 和度。 点数、坐标和顺序必须与输入文件一致,且必须恰好为六列。 在逐点投影前检查全部接收位置,任一位置不一致时直接中止,不改写输入文件。

-h

打印帮助文档。