Synthetic Offset DAS VSP in a Two‑Layer Earth: Data and Code

This repository contains the code used to generate the two‑layer (layer‑over‑half‑space) offset DAS (Distributed Acoustic Sensing) VSP (Vertical Seismic Profile) synthetic example in our manuscript. The workflow is Devito‑based for forward modeling and includes post‑processing to synthesize the DAS response from particle‑velocity outputs. Contents job.py — Main simulation driver (Devito). Sets up the two‑layer model and offset VSP geometry, runs the elastic wave simulation, and writes particle‑velocity components rec_x.pkl and rec_z.pkl to disk. utils_mod.py and model_mod.py — Adapted from Devito utilities; modified to expose and save the desired particle‑velocity outputs (vx, vz) needed for DAS synthesis. plot_data.py — Post‑processing script that reads vx/vz, synthesizes the DAS strain/strain‑rate response for a user‑defined gauge length, and produces the figures used in the paper. How to use Run job.py to generate the particle‑velocity outputs. Run plot_data.py to compute the DAS response and reproduce the figures.