CCS-monitoring_VSP-FWI is a Python-based source code package developed by Dr. Xu Shibo from ENEOS Xplora. It provides a comprehensive and modular framework for conducting forward seismic modeling and Full Waveform Inversion (FWI) specifically for Carbon Capture and Storage (CCS) monitoring applications. The package supports a range of FWI techniques, including acoustic FWI, elastic FWI, and DAS-FWI (Distributed Acoustic Sensing), making it suitable for evaluating different physical assumptions and sensing technologies in subsurface monitoring.
The toolkit is particularly optimized for Vertical Seismic Profiling (VSP) acquisition geometries, which provide high-resolution seismic data in the vicinity of injection wells. By relating synthetic modeling and inversion workflows, users can simulate realistic monitoring scenarios, assess inversion accuracy under different survey designs, and explore the detectability of time-lapse CO₂ plume migration. The code is organized for flexibility and extensibility, enabling users to perform sensitivity tests, geometry optimization, and hybrid modeling approaches, with a focus on realistic field conditions and numerical efficiency.
To install CCS-monitoring_VSP-FWI, first clone the repository:
git clone https://github.com/XU-SB/CCS-monitoring_VSP-FWI.gitThe package is designed to run on Linux systems.
💡 For full functionality, it is recommended to use JupyterLab within a virtual environment.
Core Dependencies:
pyseisnumpymatplotlibpandasscipy
Optional (but useful) Modules:
papermill(for Parallel Computation ★ Very Important)lasio(for LAS well log file reading)nbformat(for notebook handling)
Note: This package is designed for GPU-based computation.
At least one NVIDIA GPU with CUDA support is required to run the inversion code.
For optimal performance—especially for elastic Full Waveform Inversion (FWI)—multiple GPUs are recommended.
The package is organized into multiple submodules, each targeting a specific test or analysis scenario in CCS monitoring:
Perform differential waveform inversion using acoustic approximation across three representative subsurface models.
Explore the impact of various source and receiver geometries on Acoustic DD-FWI performance.
Assess the robustness of inversion under:
- Near-surface velocity variation
- Background noise
- Initial velocity model uncertainty
- Geometry optimization
Test the influence of full elastic modeling in time-lapse imaging and inversion.
Implement DAS-specific inversion by incorporating vertical particle velocity or stress component to approximate DAS responses.
Visualize inversion results and changes using high-resolution imaging maps.
For any questions, issues, or contributions, please contact the developer: 📧 [email protected] or 📧 [email protected]