Raytracing algorithm for one-way travel-time calculation in passive seismic or microseismic applications. Written in Python3.
- shooting( ) - the core function of raytracing algorithm, actually a 2d algorithm
- psraytrace( ) - the core function to call shooting function and adapt it to 3D case. Unlike calculating two-way travel-time in reflection seismic raytracing program, only one-way travel-time for transmimitted P- or S-waves is calculated. You may modify this code to calculate the PS wave travel-time.
psraytrace can be used as a module by simply importing the core function contained in folder [psmodules].
You may refer to the two examples with demo data, i.e. example2d.py and example3d.py included in folder [examples] for how to use the core function and may adapt to your case by modifying these examples.
To run this program from source codes, you need to import the following dependencies:
- Matplotlib
- Pandas
By runing the main.py script, you would have the following images displayed.
2D raytracing example
3D raytracing example
Below is the psraytrace based raytracing animation generated by RayTracer, a module from the in-house DeepListen software.
MIT License.
psraytrace was written in 2016. The sole contributor is Zhengguang Zhao, who now works for DeepListen on part-time basis.