syn
- 简介:
指定震源机制，根据动态格林函数合成三分量位移（及其空间导数）
语法
grt syn -Ggrndir -Aazimuth -S[u]scale -Ooutdir [ -Ffn/fe/fz ] [ -Mstrike/dip/rake ] [ -TMxx/Mxy/Mxz/Myy/Myz/Mzz ] [ -Dtftype/tfparams ] [ -Iodr ] [ -Jodr ] [ -N ] [ -e ] [ -s ] [ -h ]
描述
syn 模块默认合成的是脉冲型 （impulse-like）位移，单位 cm。三分量位移分别为：
Z - 垂直向上
R - 径向向外
T - 沿 R 方向顺时针旋转 90 °
保存 SAC 文件路径为 *{outdir}/{ch}.sac*，其中 ch 为分量名。
若设置 -N 则输出为 ZNE 分量，其中 Z 不变，NE分别为北向和东向。
目前合成位移支持的震源类型包括爆炸源(-S)，单力源(-F)，剪切源(-M)和矩张量源(-T)， 其中矩张量源是根据爆炸源和剪切源合成的。一次只能选择一种震源。
必选选项
- -Ggrndir
保存格林函数的目录，例如 greenfn 模块计算的单个震中距的格林函数保存在
{outdir}/{model}_{depsrc}_{deprcv}_{r}。
- -Aazimuth
方位角，单位°，北向为 0 °
- -S[u]scale
放大系数 scale ，作用在合成结果上。
对于爆炸源、剪切源和矩张量源， scale 即标量地震矩，单位为 dyne-cm 。
对于单力源， scale 单位为 dyne 。
由于标量地震矩公式中 \(M_0 = \mu\Sigma D\) 包含震源处的 \(\mu\) ， 这里为了方便设置参数，用户可以使用 -Su 表示震源处的 \(\mu\) ， 而 scale 就仅需设置 \(\Sigma D\) （面积 x 滑动距离，单位 cm^3）。
- -Ooutdir
输出目录名，不存在会自动新建。若存在同名保存路径则直接覆盖已有结果。
可选选项
- -Mstrike/dip/rake
设置剪切源参数，角度单位为 °
strike - 断层走向， 0~360
dip - 断层倾角， 0~90
rake - 滑动角， -180~180
- -TMxx/Mxy/Mxz/Myy/Myz/Mzz
设置矩张量源的 6 个分量，其中下标代表 x （北向），y （东向），z （向下）。 关于 z 方向，详见 Warning 。
- -Dtftype/tfparams
将合成结果卷积最大振幅为 1 的时间函数，支持以下选项：
-Dp/t0 - 抛物波，t0 (secs) 为持续时间，例如 -Dp/1.3
-Dt/t1/t2/t3 - 梯形波， 三个时刻 (secs) 分别表示爬升期、平台期和下降期的截止时刻， 例如 -Dt/0.1/0.2/0.4, -Dt/0.4/0.4/0.6 （退化为三角波）
-Dr/f0 - 雷克子波， f0 为主频 (Hz) 。
-D0/file - （实验性质）自定义时间函数，文件 file 仅包含一列波形幅值。
在卷积时为了匹配格林函数的时间间隔，会适当调整传入的时间函数参数。 生成的时间函数也会以 SAC 格式保存在 outdir 目录下。
使用 -D 不同参数生成的时间函数
- -Iodr
对波形做 odr 次时间积分，例如 -I1 得到阶跃型（step-like）位移。
- -Jodr
对波形做 odr 次时间微分，例如 -J1 得到速度记录。
- -N
输出的三分量为 ZNE [默认为 ZRT]。
- -e
在合成位移的同时，也合成其空间导数 \(\dfrac{\partial (u_z, u_r, u_\theta)}{\partial (z, r, \theta)}\) 。 偏导对应在文件名/变量名开头添加了
z,r和t。合成的位移空间导数的单位均为 1 。
- -s
静默输出，不在终端打印参数和进度条。
- -h
打印帮助文档。
示例
详见教程：