High‐resolution 2.5D multifocusing imaging of a crooked seismic profile in a crystalline rock environment: Results from the Larder Lake area, Ontario, Canada

Abstract

ABSTRACTA high‐resolution seismic reflection transect was acquired over a hard‐rock geological setting along an existing roadway in the Larder Lake area of the Superior Craton of Canada for the Metal Earth project in 2017. This profile, as well as other Metal earth transects, primarily aims to enhance the knowledge and to better understand the subsurface geology of the Abitibi Greenstone Belt within the Canadian Shield. The complex geological settings of the study area as well as the tribulations caused by the survey geometry have made the imaging and velocity field estimation more challenging. A recently introduced 2.5D multifocusing stacking method is one potential solution for processing crooked‐line seismic data with a poor signal‐to‐noise ratio. The 2.5D multifocusing approach offers more realistic modelling of the zero‐offset wavefield by explicitly accounting for the midpoint dispersion and cross‐dip effects. The main practical problem of the 2.5D multifocusing implementation is the simultaneous determination of the optimal wavefield parameters for each image point and time location. We address this optimization problem using a multidimensional constrained differential evolution global optimization algorithm, as this improves the efficiency and accuracy of the estimation. We have also designed an efficient processing sequence for multifocusing seismic imaging. The performance of the 2.5D multifocusing procedure has been examined on a synthetic model, generated using the same real acquisition geometry. Numerical tests demonstrate that the 2.5D multifocusing technique can produce a more focused stack with the primary reflections appearing at their poststack correct locations, and the procedure can also provide reliable estimates of interface dips. Due to the importance and difficulty of imaging the data, several conventional and advanced processing strategies have been attempted on the transect, specifically: 2D phase‐shift time migration of a dip moveout corrected stack; 2D prestack Kirchhoff time migration; swath 3D poststack migration; and our 2.5D multifocusing imaging algorithm. We found that applying the 2.5D multifocusing stacking algorithm followed by a poststack time migration approach improved the resolution of the image significantly compared to all the conventional and advanced methods and identified new reflections. The 2.5D multifocusing method also focused the steeply dipping reflections more coherently, which resolved ambiguities in geological architecture by understanding the location and continuity of structures. The method also accurately extracts 3D structural information and results in an improved signal‐to‐noise ratio.