A prototype electromagnetic vibrator, referred to here as E‐Vib, was upgraded and developed for broadband hardrock and mineral exploration seismic surveys. We selected the iron oxide mine in Blötberget, central Sweden, for a test site in 2019 for the newly developed E‐Vib because of the availability of earlier seismic datasets (from 2015 to 2016) for verification of its performance for hardrock imaging purposes. The two‐dimensional data acquisition consisted of a fixed geometry with 550 receiver locations spaced at every 5 m, employing both cabled and wireless seismic recorders, along an approximately 2.7 km long profile. The E‐Vib operated at every second receiver station (i.e. 10 m spacing) with a linear sweep of 2–180 Hz and with a peak force of 7 kN. The processing workflow took advantage of the broadband signal generated by the E‐Vib in this challenging hardrock environment with varying ground conditions. The processed seismic section shows a set of reflections associated with the known iron oxide mineralization and a major crosscutting reflection interpreted to be from a fault system likely to be crosscutting the mineralization. The broadband source data acquisition and subsequent processing helped to improve signal quality and resolution in comparison with the earlier workflows and data where a drophammer seismic source was used as the seismic source. These results suggest new possibilities for the E‐Vib source for improved targeting in hardrock geological settings.

Abstract. To discover or delineate mineral deposits and other geological features such as faults and lithological boundaries in their host rocks, seismic methods are preferred for imaging the targets at great depth. One major goal for seismic methods is to produce a reliable image of the reflectors underground given the typical discontinuous geology in crystalline environments with low signal-to-noise ratios. In this study, we investigate the usefulness of the reverse time migration (RTM) imaging algorithm in hardrock environments by applying it to a 2D dataset, which was acquired in the Ludvika mining area of central Sweden. We provide a how-to solution for applications of RTM in future and similar datasets. When using the RTM imaging technique properly, it is possible to obtain high-fidelity seismic images of the subsurface. Due to good amplitude preservation in the RTM image, the imaged reflectors provide indications to infer their geological origin. In order to obtain a reliable RTM image, we performed a detailed data pre-processing flow to deal with random noise, near-surface effects, and irregular receiver and source spacing, which can downgrade the final image if ignored. Exemplified with the Ludvika data, the resultant RTM image not only delineates the iron oxide deposits down to 1200 m depth as shown from previous studies, but also provides a better inferred ending of sheet-like mineralization. Additionally, the RTM image provides much-improved reflection of the dike and crosscutting features relative to the mineralized sheets when compared to the images produced by Kirchhoff migration in the previous studies. Two of the imaged crosscutting features are considered to be crucial when interpreting large-scale geological structures at the site and the likely disappearance of mineralization at depth. Using a field dataset acquired in hardrock environment, we demonstrate the usefulness of RTM imaging workflows for deep targeting mineral deposits.

No abstract

We develop a new data-driven algorithm that uses directional elastic wavepackets as seismic sources to image subsurface voids (i.e., cavities). Compared with a point source, the advantage of the new approach is that the wavepacket illuminates only a small volume of the medium around the ray path to significantly reduce multiple scattering effects in imaging. We take the difference of traces at identical source-receiver offsets from each of two neighboring source packets. The difference contains mainly the void scattering events but not the direct waves, the layer reflections, refractions, nor layer-related multiples. We use both P-to-P and P-to-S scattered waves to locate the voids and the results using scattered P and S waves can cross-validate each other to reduce the possibility of false detections. The directional wavepacket can be numerically synthesized using existing shot gathers so no special physical source is required. We demonstrate our method using data calculated using a boundary element method to model the seismic wavefield in an irregularly layered medium containing several empty voids. We show the robustness of our method using the same data but with 15% RMS random noise added. Furthermore, we compare our method with the RTM imaging method using the same data and find that our method provides superior results that are not dependent on the construction of a velocity model.