从这里开始,你就正式学习 gCAP 了。 gCAP 反演震源机制需要三个步骤:
计算格林函数
格林函数表示 时刻作用在 处的单位力矢量在 时刻 处产生的位移矢量。
地震的波形记录是震源项、格林函数以及仪器响应三者的卷积。 波形反演震源机制解首先是将假定的震源项(震源机制和震源时间函数)与计算出的格林函数进行卷积, 然后将卷积得到的结果与去仪器响应之后的观测波形进行对比。 因而正确计算格林函数是关键的一步。
数据预处理
数据预处理是对观测到的地震波形数据进行一些处理以符合 gCAP 对数据的要求。
gCAP 要求地震波形数据满足如下要求:
- 必须是 SAC 格式
- 三分量波形数据文件名格式为
station.[rtz]
且必须放在一个事件目录内 - 必须去除仪器响应得到位移或速度记录,单位为 cm 或 cm/s
- SAC头段中
o
、az
、dist
必须定义
反演
gCAP 会截取波形中的 Pnl 波和面波的波形,并通过网格搜索, 遍历所有可能的震源参数以得到与实际波形最接近的理论地震图。