Abstract:Seismic diffraction events, often overlooked as noise in conventional seismic imaging methods, contain valuable high-resolution information about subsurface structures and heterogeneities. Furthermore, the diffractive component of the wavefield is masked by the stronger reflection events. Although the literature presents innovative methods for separating diffracted events from the total wavefield, there are currently limited options of open-source alternatives in seismic diffraction processing and imaging. In this study, we introduce DiffraPy, an open-source Python package specifically designed for performing both conventional and diffraction-oriented migrations. The main workflow involves implementing an anti-stationary phase filter in the Kirchhoff migration for enhancing energy outside the Fresnel zone. Additionally, we propose using the semblance matrix obtained from the dip field calculation as an additional weight in conventional migration. The code provides functionality for forward modeling seismic data, allowing for new testing and experimentation. We demonstrate the general workflow of the code using a toy model and evaluate its capabilities on a synthetic salt model. In both cases, the program successfully suppressed the majority of reflected energy in the diffraction-oriented image. The final image highlights the responses of small lenses, discontinuities, and faults, providing valuable insights into the subsurface for interpretation. Despite the method sensitivity to variations in the correct velocity model, we demonstrate that our code is capable of imaging diffractors even when utilizing incorrect migration velocities of up to 5%. Our wavenumber analysis reveals the diffraction image preservation of smaller wavelengths, indicating higher resolution compared to conventional products. Future implementations may include extension to 3D datasets and improve the computational efficiency of our code. We expect our code to offer the geosciences community the first freely accessible and well-documented Python alternative for testing, reproducing, and exploring diffraction imaging examples.

Imaging 3-D faults using diffractions with modified dip-angle gathers

Wave-equation diffraction imaging using pseudo dip-angle gather

Diffraction Imaging Method Using the Pattern Operator in Migrated Dip-angle Gathers

Reflected wave least squares reverse time migration with angle illumination compensation

Separating and imaging diffractions of seismic waves in the full-azimuth dip-angle domain

3D Seismic Attributes for Structural Analysis in Compressional Context: A Case Study from Western Sichuan Basin

Diffraction Separation and Least-Squares Imaging Based on Multiscale and Multidirectional Wavefield and Image Decomposition

3-D Ultrasonic Image Reconstruction in Frequency Domain Using a Virtual Transducer Model.

A computational method for wide-azimuth 3D dip-angle gathers using Gaussian beam migration

3-D Seismic Diffraction Separation and Imaging Using the Local Rank-Reduction Method

Diffraction Separation and Imaging Using an Improved Singular Value Decomposition Method

Azimuthal imaging of rock fractures by incorporating single borehole radar and optical data

Diffraction separation and imaging based on double sparse transforms

Acoustic and Elastic Reverse-Time Migration With an Angle-Related Imaging Condition for Imaging Steeply-Dipping Structures

Three-dimensional passive seismic waveform imaging around the SAFOD site, California, using the generalized Radon transform

Wavefield decomposition in arbitrary direction and an imaging condition based on stratigraphic dip

DiffraPy: An open-source Python software for seismic diffraction imaging

Frechet Kernels for Imaging Regional Earth Structure Based on Three-Dimensional Reference Models

A Distortion Matrix Framework for High-Resolution Passive Seismic 3-D Imaging: Application to the San Jacinto Fault Zone, California

Method to quantify the dynamic near-fault full-field stress evolution associated with rough fault shear deformation based on 3D printed models and photoelastic measurements

Three‐Dimensional Teleseismic Elastic Reverse‐Time Migration With Deconvolution Imaging Condition and Its Application to Southwest Japan