Lupei Zhu教授的FK工具包：水平分层模型格林函数计算与地震图合成教程

FK是陆培泽教授开发的一款开源代码包，专用于计算水平分层模型下的理论格林函数并合成理论地震图。它主要由以下几个关键部分组成： 1. **fk**：这是主程序，负责计算格林函数，是整个工具的核心功能。 2. **st_fk**：用于计算静态格林函数，通常在处理静校正或特定类型的地震分析时使用。 3. **trav**：辅助程序，用于计算P波和S波的初至时间，这对于确定地震波传播路径至关重要。 4. **sachd**：用于编辑SAC头段，可能与数据预处理或格式转换有关。 5. **fk.pl**：这个脚本是fk、st_fk、trav和sachd功能的整合，提供了一站式的使用体验，用户一般通过它来执行主要任务。 使用fk.pl时，需要准备一个一维水平分层速度模型作为输入，该模型格式包括以下内容： - 层的厚度（km）：第1列的数值。 - S波速度（km/s）：第2列，必要信息。 - 波速比或P波速度（可选）：第3列，如果使用第一种格式（-Mmodel/depth/k），则此列包含波速比，否则默认P波速度是S波速度的两倍。 - 密度（g/cm³）：第4列，可选，默认使用经验公式ρ=0.77+0.32*vp。 - S波和P波的Q值（可选）：第5和第6列，若未指定，分别取Qs=500和Qp=2*Qs。 在使用fk.pl时，模型文件名（-Mmodel/depth）的输入方式决定了速度模型的格式： - 使用-Mmodel/depth：模型以第二种格式输入，即第三列是波速比。 - 使用-Mmodel/depth/k：模型的第一种格式输入，即波速比单独列出。 - 使用-Mmodel/depth/f：模型需进行展平变换，适用于震中距较大的情况。 值得注意的是，fk.pl接受的速度模型格式是先Vs（S波速度）后Vp（P波速度），而fk内部处理的模型顺序则是相反的，因此在传递模型时需要确保正确对应。此外，如果第一层的厚度为零，那么这一层将被解析为上半空间的参数；若厚度不为零，则上半空间被视为真空。 FK是一个强大的地震学工具，能够处理复杂的地震模型计算，并提供了用户友好的脚本接口，使得地震学家和地球物理学家能够在处理地震数据时更加高效地利用该工具。