Published January 17, 2025 | Version 1.0.0

Dynamic Time Warping Proximity Analysis (DTW-PA)

  • 1. ROR icon Universidade de São Paulo
  • 2. Research Centre for Greenhouse Gas Innovation (RCGI)
  • 3. Universidade Federal de São Carlos (UFSCar)

Description

MANUSCRIPT SUPPLEMENTARY MATERIAL

This supplementary material is provided for the manuscript titled Full-waveform inversion cycle-skipping mitigation with dynamic time warping, Part 1: a method for a proximity analysis between the current and the observed data, submitted to the journal Pure and Applied Geophysics. For additional context and details, please refer to the manuscript.

Table of contents (English)

SUPPLEMENTARY RESULTS

The files provided include the results of the FWI numerical simulations presented in Table 2 of the manuscript. Due to Zenodo's size limitations, some files could not be uploaded, as:

  • Seismograms from the l_bfgs_b_... folders.
  • DTW best path results from the l_bfgs_b_..., initial_proximity_eval, and proximity_eval... folders.

The l_bfgs_b_... folders contain data for each multiscale step. The folder name indicates the low-pass frequency filter (FF) applied at each step and the intermediate step (IS) (e.g., l_bfgs_b_FF-10.0Hz_IS-1).

The initial_proximity_eval folders contain files from the Dynamic Time Warping Proximity Analysis (DTW-PA) conducted prior to the first multiscale step using a 5 Hz low-pass filter. These files assess the initial model's proximity to the target data.

The proximity_eval... folders contain DTW-PA results used to evaluate the low-pass frequency filters for subsequent multiscale steps. Each folder name specifies the filter tested (e.g., proximity_eval-10.0Hz).

CODE

Each simulation folder contains a code subfolder, which includes the scripts and code used to run the numerical simulations. The following files are included in each folder:

  • A .py script for running the FWI numerical simulation.
  • ..._pa_plt.py to generate DTW-PA figures.
  • ..._plt.py to produce additional figures.
  • A .npy velocity model, when needed.
  • The requirements.txt.

Additionally, .slurm files are provided, as the simulations were executed on a cluster using the SLURM job scheduler. These files are supplied AS-IS and must be adjusted to fit your specific computational environment.

Technical info (English)

SOFTWARE INSTALLATION

To run the provided codes, follow these steps to set up a virtual Python environment:
  1. Ensure Python 3.8 is installed and set as your default Python3 version.
  2. Install virtualenv if it is not already installed.
  3. Execute the following commands:
python3 -m venv my_venv
source my_venv/bin/activate
python -m pip install --upgrade pip
pip install -r requirements.txt
In each code folder is a requirements.txt file.

Notes (English)

HARDWARE REQUIREMENTS

These simulations require substantial computational resources. Running them on a personal computer is not recommended.
 
We used five Intel Xeon Gold 6148 nodes for the numerical simulations, each running at 2.4 GHz with 27 MB cache and 10.4 GT/s. Each node had 192 GB of DDR4 2,666 MHz ECC Reg memory. Due to the scheduler using one node, only four were used for FWI computations. Similarly, the scheduler limited the number of workers employed for shot parallelisation.

Files

bp_salt_model.zip

Files (12.2 GB)

Name Size
md5:b245b5abf5a36c560fa6fcd610ab36d4
5.4 GB Preview Download
md5:eb48c4a2897cd5698e4cb0f49abf6dee
1.5 GB Preview Download
md5:6aed5bac48e80e8439bffa386c6fabac
5.3 GB Preview Download

Additional details

Software

Repository URL
https://github.com/eikmeier-cn/dtw-pa-isa
Programming language
Python