AI-Aided GPR Data Multipath Summation Using x-t Stacking Weights

The Ground Penetrating Radar (GPR) method can image dielectric discontinuities in subsurface structures, which cause the reflection of electromagnetic (EM) waves. These discontinuities are imaged as reflectors in GPR sections, often distorted by diffracted energy. To focus the diffracted energy within the GPR sections, migration is commonly used. The migration velocity of GPR data is a low-wavenumber attribute crucial for effective migration. Obtaining a migration velocity model, typically close to a Root Mean Square (RMS) model, from zero-offset (ZO) data requires analysis of the available diffractions, whose density and (x, t) coverage are random. Thus, the accuracy and efficiency of such a velocity model, whether for migration or interval velocity model estimation, are not guaranteed. An alternative is the multipath summation method, which involves the weighted stacking of constant velocity migrated sections. Each stacked section contributes to the final stack, weighted by a scalar value dependent on the constant velocity value used and its relation to its estimated mean velocity of the section. This method effectively focuses the GPR diffractions in the presence of low heterogeneity. However, when the EM velocity varies dramatically, 2D weights are needed. In this study, with the aid of an Artificial Intelligence (AI) algorithm that detects diffractions and uses their kinematic information, we generate a diffraction velocity model. This model is then used to assign 2D weights for the weighted multipath summation, aiming to focus the scattered energy within the GPR section. We describe this methodology and demonstrate its application in enhancing the lateral continuity of reflections. We compare it with the 1D multipath summation using simulated data and present its application on marble assessment GPR data for imaging cracks and discontinuities in the subsurface structure.