Software developed by Geotem Ingeniería in Python 3.7.x to compute the spectral ratio (H/V) using seismic noise data and also to compute the inversion of the H/V using the Diffuse Field Assumption.
The processing of ambient noise vibrations is made with the stacking and the mean of the selected windows considering the diffuse field assumption.
The inversion of the H/V runs an external C++ executable file. The modeler of the H/V phenomena used in the inversion algorithm is the one described in Sánchez-Sesma, et al. 2011.
The code is developed and tested on Windows and Ubuntu using Python 3.7.x. You have two options of installation, through the creation of an Anaconda enviroment or using the pip method.First we need to create an enviroment in the Anaconda prompt with the following specifications:
conda create -n terrawarehv python=3.7 numpy matplotlib=3.2.2 scipy=1.5.0 pandas=1.0.5 pyside2
Then we need to activate the enviroment
conda activate terrawarehv
In addition to Python and setting it as a Path variable, you need to install the following packages using pip.
pip install numpy==1.19.0
pip install matplotlib==3.2.2
pip install scipy==1.5.0
pip install pandas==1.0.5
pip install pyside2==5.15.0
Download or clone this repository, then open main.py file. For the anaconda option we need to activate the "terrawarehv" enviroment and run the file in it.
python main.py
Once the main window is opened, load the files in the folder "Ejemplos" with "Load V.A. file" option and double click a sounding to work on it. Then change the parameters and process it.
The processing parameters are described in the next table:
| Processing Parameters | Description |
|---|---|
| Normalización | Type of normalization in frequencies |
| Window | Type of tappered window |
| Traslape | Overlap percentage |
| Factor de Taper | Taper factor if the window is Tukey type |
| OneBit | Applies a OneBit normalization in time |
| Quitar tendencia | Remove mean trend of the signal |
| Muestreo tiempo | Time sampling rate. Takes the signal default sampling rate |
| Muestreo frecuencia | Frecuency sampling rate for the spectral ratio. |
| Segundos de ventana | Window length in seconds |
| Konno-Ohmachi | Konno-Ohmachi filter factor that goes from 0 to 90, being 90 the raw ratio and 0 the most smoothed ratio. |
| Frecuencias para H/V | Range of frequencies to visualize the spectral ratio |
| Filtro de Butterworth | Cut-off frequencies to apply a Butterworth filter if needed |
| H/V direccionales | Spectral ratio of the N-S and E-W components |
| Espectros | Amplitude spectrum of the raw signal |
The initial model parameters depend on the results you are expecting:
Modify the parameters until you get the results you expect.
Please read CONTRIBUTING.md for details on our code of conduct, and the process for submitting pull requests to us.
- Michelle Sainos Vizuett - Developer
- Juan Carlos Colchado Casas - Developer
- Francisco Huante Arana - Consultant
- Alfredo Sánchez Galindo - Chief Software Manager
- Francisco José Sánchez Sesma - External Consultant