An example of Receiver function inversion of Hunan province
- Receiver_function_HunanProvince
- Table of Contents
- Background
- Install
- Manual
- Maintainers
- Contributing
- Contributors
- License
This is my first project about Geophysics (Receiver function) which is managed in Github. During the research, I found that receiver function is much popular in the inversion of nature seismic, but there is not a complete process of receiver function inversion. So I decided to create a example of receiver function inversion. I hope that every small process during the inversion could be shown in this project, and it may help for young scientists who work on receiver function.
This project is base on the data of Hunan Province, China. According to the restrictions of data, only limited data is provided. Besides, this package contain serveral preseved code of specific research group, and I will provide some substitute code on Github. All the inversion tools are public, you can find these code or packages below.
This package consists of several software, make sure that you have install these software before:
I strongly recommend Obspy for data preprocessing. The scripts written by python are read much easier than those written by sh/bash. If you install Obspy by anaconda, please make sure that the enviorment setting is correct.
In order to run some Linux script, this serveral programs was demand
git clone https://github.com/SuwenJunliu/Receiver_function_HunanProvince.git
Accoring to the purpose of this project, I will introduce the package structure, including what every folder is, and how to use them. Next, the whole inversion workflow will be present. You could follow the instruction of "QuickStart" and realize them easily. For more details, you could find the "readme" of corresponding folder.
The folder "./data" consists the data of Hunan Province in SAC format.
cd /prepareEventFile/
python
Then, you will find the event index in this folder.
python
After that, the data of nature seismic event was created.
The first thing for RF inversion is preparing SAC file of nature seismic event. Usually, I use Obspy to download the event index. The code could be found in
./prepareEventFile/
Note that if your study area is in China, and the data was gather recently, the event may not be included by ISC. In this situation, you can download the event index from CENC.
Once you prepare the event index, you could cut the event data form original record. Usually, the original record last one hour or much longer. In this project, I use python script to cut the event. The code could be found in
./prepareEventFile/
Feel free to dive in! Open an issue or submit PRs.
GPL3.0 © Junliu Suwen