反演参数
反演各个参数的意义在朱老师的脚本 src/gcap/cap.pl
中有官方的说明。
另外,请读者一定明确设置参数究竟是在做什么。
src/gcap/cap.pl
就是 gCAP 的反演脚本。在这个脚本中,各个反演参数已经有一个默认值。
设置反演参数就是去改这些默认值。学习参数设置时首先要注意了解什么情况需要设置参数,修改默认值。
这里对我有充分把握的参数进行说明。自己没有用过的参数,我虽然能或多或少揣摩意思,但也只是英译汉的层次,为了保证正确性,暂时不写,只用原英文说明代替。
-B
output misfit errors of all solutions for bootstrapping late ($bootstrap).
-Cf1_pnl/f2_pnl/f1_sw/f2_sw 滤波拐角频率
在反演前,gCAP 会先对观测数据做滤波。
-Dw1/p1/p2 Pnl 和面波的权重
w1 是 Pnl 的权重。p1 和 p2 是目标函数中 Pnl 和面波的归一化距离的指数。 If p1 or p2 is negative, all traces will be normalized.
-Fthr
include first-motion data in the search. thr is the threshold ($fm_thr). The first motion data are specified in $weight. The polaritiescan be specified using +-1 for P, +-2 for SV, and +-3 for SH afte the station name, e.g. LHSA/+1/-3 means that P is up and SH is CCW. The Green functions need to have take-off angles stored in the SAC header.
-Ggreen 格林函数库位置
green 是格林函数库所在的路径。
除非你没有像例子中那样把格林函数就放在 Glib/
下,否则 不要 设置这个参数。
-Hdt 采样周期
dt 是采样周期必须和数据预处理的 RESAMPLE 参数和格林函数的 DT 二者相同。一定要设置,除非采样周期恰好是 0.1。
-Idd[/dm] 搜索步长
dd 是断层面三个参数(strike、dip 和 rake)的搜索步长,dm 是震级的搜索步长。默认值是 10/0.1,如果想用其他步长则需要设置。
If dm<0, the gain of each station will be determined by inversion.
我没有用过负数的情况。
-Jiso[/diso[/clvd[/dclvd]]] 非双力偶成分设置
只有需要搜索非双力偶成分,才需要设置这个参数。iso 和 clvd 是 ISO 和 CLVD 的反演初始值。diso 和 dclvd 是 ISO 和 CLVD 搜索的步长。 如果需要搜索非双力偶成分,通常设置为 0/0.1/0/0.1。
-Kvpvs[/mu]
use the vpvs ratio and mu at the source to compute potency tensor parameters ISO and P0. (0/0, off).
-Ldura 震源持续时间
地震波形是格林函数和震源时间函数(不是震源持续时间)的卷积。 在 gCAP 中,认为震源时间函数是一个等腰梯形。等腰梯形由底边长和梯形上升阶段占据的底边长的比例决定。 震源持续时间就是这里的底边边长。dura 为震源持续时间的长度(单位:秒)。也可以用 sac 文件来指定,我没用过这种方式。
如果不设置这个参数,gCAP会使用 -M 参数中的震级,用经验公式计算。
-Mmodel_depth/mag 模型、深度、震级
gCAP 反演程序唯一的必选选项,必须设置
model 是模型名称。depth 是深度。gCAP 会在格林函数库的路径下去寻找相应模型、相应深度的格林函数文件。 mag 是震级,用预估的矩震级值。
-Nn
repeat the inversion n times and discard bad traces ($repeat).
-O 显示设置
没有参数值,只要有了 -O,就不会进行反演,而是把传给二进制程序 cap
的各种设置输出到屏幕。
-P 绘图设置
如果画图没有问题,就不需要设置。Yscale: 第一个 trace 的振幅值(英寸)。Xscale: 每英寸的时间长度。
append k if one wants to keep those waveforms.
-Qnof
number of freedom per sample ($nof)
-R 网格搜索的范围
默认值是 0/360/0/90/-90/90。strike2 建议改为 180,减少计算量。
-S 最大移动时间
Pnl 和面波的移动最大时间(s1 和 s2),tie=0:SV 和 SH 相互独立移动,tie=0.5:SV 和 SH 的移动时间保持一致
-T 时窗长度
Pnl 和面波的时窗长度
-Udirct
directivity, specify rupture direction on the fault plane (off).
-V 视速度
Pnl、Love 和 Rayleigh 的视速度
-Wi 反演的物理量
通常 i = 1,默认值是 1。如果是 1,则观测波形的振幅是速度,2 是位移。
-Xn
output other local minimums whose misfit-min<n*sigma ($mltp).
-Zstring 权重文件
string 是权重文件的文件名,如果不设置,默认值是 weight.dat。