✍️ 朱邓达  •  📅 2025-12-19

(9) 计算海底模型中的动态位移和应力

下载示例： ex09.tar.gz

海底模型中存在固液界面，相关公式见 含液体层的反射透射系数矩阵 。

示例海底模型如下：

1.0   1.5   0.0   1.1   300   150
2.0   3.0   1.5   2.0   300   150
3.0   5.0   2.8   2.5   500   250
5.0   6.0   3.5   2.8   600   300
10.0   7.0   4.0   3.0   800   400
10.0   7.8   4.5   3.2   1000   500
10.0   8.2   4.7   3.3   1200   600
10.0   8.5   4.8   3.4   1500   700

Shell 脚本如下。其中有必要说明的是，固液边界处仅有垂向位移和垂向应力连续， 其余分量不连续，此时若将震源或台站正好放置在固液边界处则会产生歧义，程序会警告并退出。 因此对于此情况需适当在深度上设置一个极小量的偏离，以明确表明其所处层位。

dist=30
depsrc=3

nt=1500
dt=0.05

modname="mod_sea"
out="GRN_sea"

# =================================================
# compute displacement in solid
deprcv=1.001
grt greenfn -M${modname} -O${out} -N${nt}/${dt} -D${depsrc}/${deprcv} -R${dist}

# =================================================
# compute displacement and its derivatives in liquid
deprcv=0.999
grt greenfn -M${modname} -O${out} -N${nt}/${dt} -D${depsrc}/${deprcv} -R${dist} -e
# just give a test focal mechanism
grt syn -G${out}/${modname}_${depsrc}_${deprcv}_${dist} -A40 -S1e20 -M100/55/120 -Osyn -e
# compute stress
grt stress syn

  • 位移

../../_images/disp.svg

  • 应力张量（由于接收点位于液体层，仅有完全相同的 ZZ、RR、TT 分量，即液体压强）

../../_images/stress1.svg