✍️ 朱邓达  •  📅 2025-09-22

static_syn

简介:

指定震源机制,根据静态格林函数合成三分量位移(及其空间导数)

语法

grt static syn -Gingrid -S[u]scale -Ooutgrid [ -Xx1/x2/dx ] [ -Yy1/y2/dy ] [ -Ffn/fe/fz ] [ -Mstrike/dip/rake ] [ -TMxx/Mxy/Mxz/Myy/Myz/Mzz ] [ -N ] [ -e ] [ -h ]

描述

关于输出分量、震源机制的部分与 syn 模块相同,但输入和输出均使用 NetCDF 网格。

static syn 模块计算时使用的 XY 网格默认继承 static_greenfn 模块设置的 -X-Y 。 也支持设置 -X-Y 来指定新的 XY 网格,此时每个节点将使用最近震中距的格林函数来进行合成, 因此这是一种近似方法,程序会打印震中距之差的统计信息,以方便评估近似效果。 但近似好处是一旦静态格林函数计算好,合成阶段使用新的网格也可以复用格林函数,无需再计算。

必选选项

-Gingrid

存有静态格林函数的网格文件。

-S[u]scale

放大系数 scale ,作用在合成结果上。

  • 对于爆炸源、剪切源和矩张量源, scale 即标量地震矩,单位为 dyne-cm 。

  • 对于单力源, scale 单位为 dyne 。

  • 由于标量地震矩公式中 \(M_0 = \mu\Sigma D\) 包含震源处的 \(\mu\) , 这里为了方便设置参数,用户可以使用 -Su 表示震源处的 \(\mu\) , 而 scale 就仅需设置 \(\Sigma D\) (面积 x 滑动距离,单位 cm^3)。

-Ooutgrid

输出网格文件,若同名文件已存在,则销毁后新建。

可选选项

-Xx1/x2/dx

指定北向等距网格点 (km) 。 x1,x2 分别为最小、最大值, dx 为间隔。

-Yy1/y2/dy

指定东向等距网格点 (km) 。 y1,y2 分别为最小、最大值, dy 为间隔。

-Ffn/fe/fz

设置单力源参数,每个方向的力将乘上 -S 中定义的放大系数。关于 z 方向,详见 Warning

  • fn - 北向力

  • fe - 东向力

  • fz - 垂直向下力

-Mstrike/dip[/rake]

设置剪切源/张裂源参数,角度单位为 °

  • strike - 断层走向, 0~360

  • dip - 断层倾角, 0~90

  • rake - 滑动角, -180~180 [不指定 rake 即设置张裂源]

-TMxx/Mxy/Mxz/Myy/Myz/Mzz

设置矩张量源的 6 个分量,其中下标代表 x (北向),y (东向),z (向下)。 关于 z 方向,详见 Warning

-N

输出的三分量为 ZNE [默认为 ZRT]。

-e

在合成位移的同时,也合成其空间导数 \(\dfrac{\partial (u_z, u_r, u_\theta)}{\partial (z, r, \theta)}\) 。 偏导对应在文件名/变量名开头添加了 z , rt 。合成的位移空间导数的单位均为 1 。

-h

打印帮助文档。

示例

详见教程: