简介
裁剪-粘贴法(Cut And Paste,简称 CAP) 是一种将事件波形切割成多个波形片段来反演震源机制解的方法。 该方法主要由 朱露培 发展和实现。 一般而言,CAP 适用于矩震级 3 级以上,7 级以下, 在 500 公里内至少有 3 个宽频带波形数据(有仪器响应文件)的地震事件。
朱露培发布的 CAP 程序有三个不同的版本:CAP、gCAP 和 gCAP3D。
- CAP 是早期实现的代码,可反演双力偶解。该版本仅通过私下交流传播,并未在网络上开源;
- gCAP 在 CAP 的基础上增加了反演完整矩张量解的功能。gCAP 于 2013 年在网络上开源。
- gCAP3D 在 gCAP 的基础上可使用三维模型的格林函数进行反演。gCAP3D 于 2016 年在网络上开源。
另外,很多学者发布了自己的 CAP 版本,各有特色。
本项目基于 gCAP 建立。 之所以命名为 Oh My CAP,是模仿 Oh My Zsh 的命名方式,意思为 gCAP 的配置脚本。 本项目提供使用 gCAP 的教程和调用脚本。 下面不完全罗列本项目的优势:
- 对 gCAP 的核心代码仅微量修改。用 Oh My CAP,本质上就是用朱露培的官方原版 gCAP;
- 提供安装程序源文件、自动安装脚本和手把手的安装指南;
- 提供手把手的示例演示,和每一个步骤的详细讲解;
- 提供调用 gCAP 的现成 Perl 脚本,易用性上,远超官方版本;
- 可以依托 Github 的版本对照功能清楚查看对原 gCAP 代码做了哪些具体修改;
Oh My CAP 2.0 的更新
本项目的第一个版本 v1.0 发布在 2016 年 10 月 14 日。 目前,Oh My CAP 最新是 2.0 版本。 这是 2018 年 4 月的 v1.2 版之后近五年来的第一次版本更新。 Oh My CAP 不再依赖 gmt4 和 pssac 了,而只使用 gmt6 进行做图。 因为改变了软件依赖,版本号直接前进到 2.0。 这一次更新还纠正了 fk 的部分错误。 我相信,这一次的版本更新会为本项目重注一次生命力。
授权协议与论文引用
本项目中所使用的 FK 和 gCAP 的源码修改自朱露培的原始代码,按照 GPL 协议授权。 你可以免费获取、使用、修改和再发布代码。 如果你修改了 FK 和 gCAP 的,也必须按照 GPL 协议公开。 本项目中的其余源码以及网页内容采用更加宽松的 Apache 协议, 即在尊重本项目的署名权的前提下,可以选择不公开自己的修改。
如果本项目帮助了你,你应该引用我的这篇论文:
王亮, 薛霆虓, 季海磊. 2016. 集集强余震震源机制解分析[J]. 地球物理学进展. 31(5): 1998-2004
Wang L, Xue T X, Ji H L. 2016. Focal mechanisms of Taiwan ChiChi earthquake aftershocks[J]. Progress in Geophysics. 31(5): 1998-2004
我这篇论文不是英文的,也没有发表在顶级期刊。 三位作者都不是大牛。 因为这些,你就用我的程序,而不引用我的论文?
做人不好不讲良心的。
如果你公开引用了我的论文,无论是否是期刊发表,都请告诉我。这是对我的莫大鼓励。
你可以在文献介绍中了解文献阅读和文献引用的更多信息。
学习指南
读者应该具备的最基本的理论知识准备有:
- 弹性回跳理论
- 双力偶和沙滩球。
- 若要做矩张量解,还需要知道矩张量。
如果你完全不具备上述知识储备,你应该暂停本文档的学习,先补一补课。
通常,在学习地震学程序之前,应该先充分阅读文章,理解方法。 但是,阅读本文档基本可以代替阅读文献。 本文档会穿插应该在文献中获得的知识。 所以读者只需要在具备前述的弹性回跳等基本知识后, 按照顺序照本文档介绍的做就好了。
可以通过访问网站,在线阅读本文档。 也可以离线阅读。 文档在项目的根路径下的docs目录中。
文件树
本项目的整体文件树如下:
oh-my-cap
├── docs 文档
├── example 例子
├── src 源文件
├── Glib 格林函数库
├── LICENSE 授权协议
└── README.md 自述文件