Seismic Imaging in Crystalline Terrains of the Superior Province, Canada
Abstract
The imaging of structures in sedimentary basins has typically been carried out using conventional seismic reflection techniques. In complex geological architectures, such as those within the Superior Province in Canada, conventional seismic imaging yields suboptimal results due to factors stemming from complex energy scattering and crooked seismic geometries. Full-waveform inversion (FWI) is a nonlinear inverse technique capable of retrieving quantitative images of velocity structures. Judicious data preconditioning, and multiscale inversion strategies allow FWI to effectively estimate velocity variations in the first few kilometers of the subsurface. Non-conventional seismic reflection processing, based on azimuthal binning and enhanced migration velocity models, improves energy focusing and imaging quality of structures at depth.
I demonstrate that the serious nonlinearity of FWI in crystalline zones (in Larder Lake and in the northeastern portion of the Sudbury Structure), is alleviated by implementing a “multiscale layer-stripping” strategy. The strategy uses a i) combination of explosive and vibroseis sources to retrieve low-wavenumber features of the velocity background, ii) hierarchical minimization of logarithmic phase-only and conventional phase-amplitude residuals to mitigate large dynamic variations within the data, and iii) progressive inclusion of higher frequencies and late arrivals to obtain a natural transition between low- and high- wavenumber features. I illustrate that the implementation of optimum binning strategies effectively enhances signal alignment, generating in-phase stacked sections. The use of near-surface FWI P-wave velocity estimations in the construction of migration velocity models, improves the strength and continuity of reflections at depth.
With the application of these imaging techniques, I unravel the geometry of prominent structures in structurally-controlled mineralized zones. In Larder Lake, the Larder Lake-Cadillac Deformation Zone (LLCDZ), and the Lincoln-Nipissing Shear Zone (LNSZ), are imaged. The LNSZ is retrieved as a north-dipping fault, extending to depths not previously identified of 8 km. The Sudbury velocity model reveals the internal character of the structure, agreeing with known geology and borehole information. New estimates of thicknesses and dips are given for the northeastern portion of the Sudbury Structure.
Keywords: Full-waveform inversion, Finite-difference modeling, Seismic imaging, Azimuthal binning, Energy focusing, Reflection processing, Crystalline environments, Larder Lake, Sudbury.