Reflection imaging of complex geology in crystalline environment using virtual-source seismology: case study from the Kylylahti polymetallic mine, Finland

For the first time, we apply a full-scale 3D seismic virtual-source survey (VSS) for the purpose of near-mine mineral exploration. The data was acquired directly above the Kylylahti underground mine in Finland. Recorded ambient noise (AN) data is characterized using power-spectral density (PSD) and beamforming. Data has most energy at frequencies 25–90 Hz and arrivals with velocities higher than 4 km/s have wide range of azimuths. Based on the PSD and beamforming results, we created 10-days subset of AN recordings that were dominated by multi-azimuth high-velocity arrivals. We use illumination-diagnosis technique and location procedure to show that the AN recordings associated with high apparent velocities are related to body-wave events. Next, we produce 994 virtual-source gathers by applying seismic-interferometry processing by cross-correlating AN at all receivers resulting in full 3D VSS. We apply standard 3D time-domain reflection seismic data processing and imaging using both a selectively stacked subset and full passive data, and validate the results against a pre-existing detailed geological information and 3D active-source survey data processed in the same way as the passive data. The resulting post-stack migrated sections show agreement of reflections between the passive and active data and indicate that VSS provide images where the active-source data are not available due to terrain restrictions. We conclude that while the all-noise approach provides some higher quality reflections related to the inner geological contacts within the target formation and the general dipping trend of the formation, the selected subset is most efficient in resolving the base of formation.

Evaluation of Lake Sediment Thickness from Water-Borne Electrical Resistivity Tomography Data

Lakes are integrators of past climate and ecological change. This information is stored in the sediment record at the lake bottom, and to make it available for paleoclimate research, potential target sites with undisturbed and continuous sediment sequences need to be identified. Different geophysical methods are suitable to identify, explore, and characterize sediment layers prior to sediment core recovery. Due to the high resolution, reflection seismic methods have become standard for this purpose. However, seismic measurements cannot always provide a comprehensive image of lake-bottom sediments, e.g., due to lacking seismic contrasts between geological units or high attenuation of seismic waves. Here, we developed and tested a complementary method based on water-borne electrical-resistivity tomography (ERT) measurements. Our setup consisted of 13 floating electrodes (at 5 m spacing) used to collect ERT data with a dipole–dipole configuration. We used a 1D inversion to adjust a layered-earth model, which facilitates the implementation of constraints on water depth, water resistivity, and sediment resistivity as a priori information. The first two parameters were readily obtained from the echo-sounder and conductivity-probe measurements. The resistivity of sediment samples can also be determined in the laboratory. We applied this approach to process ERT data collected on a lake in southern Mexico. The direct comparison of ERT data with reflection seismic data collected with a sub-bottom profiler (SBP) showed that we can significantly improve the sediment-thickness estimates compared to unconstrained 2D inversions. Down to water depths of 20 m, our sediment thickness estimates were close to the sediment thickness derived from collocated SBP seismograms. Our approach represents an implementation of ERT measurements on lakes and complements the standard lake-bottom exploration by reflection seismic methods.

Lateral migration of large sedimentary bodies in a deep-marine system offshore of Argentina

Contourite features are increasingly identified in seismic data, but the mechanisms controlling their evolution remain poorly understood. Using 2D multichannel reflection seismic and well data, this study describes large Oligocene- to middle Miocene-aged sedimentary bodies that show prominent lateral migration along the base of the Argentine slope. These form part of a contourite depositional system with four morphological elements: a plastered drift, a contourite channel, an asymmetric mounded drift, and an erosive surface. The features appear within four seismic units (SU1–SU4) bounded by discontinuities. Their sedimentary stacking patterns indicate three evolutionary stages: an onset stage (I) (~ 34–25 Ma), a growth stage (II) (~ 25–14 Ma), and (III) a burial stage (< 14 Ma). The system reveals that lateral migration of large sedimentary bodies is not only confined to shallow or littoral marine environments and demonstrates how bottom currents and secondary oceanographic processes influence contourite morphologies. Two cores of a single water mass, in this case, the Antarctic Bottom Water and its upper interface, may drive upslope migration of asymmetric mounded drifts. Seismic images also show evidence of recirculating bottom currents which have modulated the system’s evolution. Elucidation of these novel processes will enhance basin analysis and palaeoceanographic reconstructions.

Three-dimensional shear velocity structure of the Mauleon and Arzacq basins (Western Pyrenees)

We present a 3-D shear wave velocity model of the Mauleon and Arzacq basins from the surface down to 10~km depth. This model is obtained by inverting phase velocity maps for periods from 2 to 9~s measured on coherent surface wavefronts extracted from ambient seismic noise by matched filtering. This new model, which is found in good agreement with local earthquake tomography, reveals the architecture of the Mauleon and Arzacq basins which were poorly imaged by conventional reflection seismic data. Combining these new tomographic images with surface and subsurface geological information allows us to trace major orogenic structures from the basement to the surface. In the basin, the models are successfully imaging first-order folds and thrusts at kilometric scale. The velocity structure within the basement and its geometrical relationship with the base of inverted rift basins supports a progressive northward exhumation of deep crustal and mantle rocks in the hanging wall of north-vergent Pyrenean thrusts. Our tomographic models image in 3-D orogen-perpendicular structures responsible for crustal segmentation as the Saison and Barlanes transfer zones. We propose that these steep structures consist in tear faults that accommodate the deepening of the Mauleon basin basement from west to east. To the west, this basement made of former hyper-extended rift domains (including mantle rocks) is anomalously sampled within the hanging-wall of north-directed orogenic thrusts, explaining its shallow attitude and its best preservation in comparison to the eastern segment of the study area. Eastward, the vertical shift of the basement makes that the former Mauleon basin hyper-extended rift basement remained in a footwall situation in respect of orogenic thrust and was underthrust. The comparison of the tomographic models obtained with surface wave tomography and local earthquake tomography shows that each approach has its own advantages and shortcomings but also that they are very complementary in nature, which would suggest to combine them in joint inversions to further improve passive imaging of the shallow crust and sedimentary basins.

Estimation of Shallow Sulphur Deposit Resources Based on Reflection Seismic Studies and Well Logging

In this article we present a novel method for the estimation of sulphur deposit resources based on high-resolution shallow reflection seismic survey and well logging. The study area was sited in the northern part of the Carpathian Foredeep (SE Poland), where sulphur ore occurs in carbonate rocks at a depth of about 120 m, with a thickness of approximately 25 m. The results of many years of seismic monitoring performed in the area of the sulphur deposit allowed us to determine the quantitative relationships between the amplitude of the seismic signal reflected from the top of the deposit and its petrophysical parameters such as porosity and sulphur content. The method of evaluating sulphur deposit is based on extensive statistics concerning the reservoir properties obtained from borehole data. We also discuss a methodology for conducting field acquisition and processing of seismic data in the aspect of mapping the actual amplitudes of the signal reflected from the top of a deposit. The results of estimating the abundance of carbonate sulphur deposits are presented based on the example of a seismic cross-section from the Osiek sulphur mine. Obtained results allow indicating the most prospective zones suitable for exploitation.