Separation of intrinsic and scattering Q by envelope inversion
Qopen is a script in seismology that determines shear wave scattering and intrinsic attenuation parameters by inversion of seismogram envelopes.
The method will be described in an upcoming publication. Please read our DGG 2015 poster in the meantime.
Dependencies of Qopen are:
- ObsPy>=0.10
- future, NumPy>=1.7, SciPy>=0.11 and matplotlib>=1.3 (itself dependencies of ObsPy)
- statsmodels
- joblib (for multi-core support)
It's probably best to install ObsPy first and then let pip take care of the rest. Install Qopen and its dependencies with :
pip install qopen
Qopen provides the two scripts qopen and qopen-runtests. The installation can be tested with the second script:
qopen-runtests
The code is run by the installed command line script qopen. A tutorial can be created with the appropriate flag:
qopen --tutorial
This command copies an example configuration file in JSON format and the corresponding data files into the current directory. The configuration file is heavily commented and should be rather self-explanatory. Now you can perform the inversion by simply running :
qopen
which will calculate the results and create different plots.
To use the script with your own data you need 1. an inventory (StationXML or other ObsPy readable format) of your stations, 2. the earthquake catalog (QuakeML or other ObsPy readable format) preferable with P and S picks and 3. the waveforms. Waveforms may exist in data files of various formats or can be fetched from a webservice. A custom solution for waveform retrieval is also possible (e.g. mixing of data files and web requests). An example configuration file can be created with :
qopen --create-config
This file has to be adapted to your needs (time window selection, etc.). The inversion is started by simply running qopen again.