Seismic methods have been used in sedimentary environment for almost 80 years. During that time, exploration geophysicists have developed a number of techniques to handle specific aspects of working in sedimentary areas. This is not the case for studies in the hardrock environment, where significantly less time and money have been invested on seismic investigations. Therefore, there is still a need to develop the right techniques appropriate for working in hardrock environments. The research presented here, covers aspects of acquisition, processing and interpretation in hardrock environments. A cost-effective alternative for two-dimensional data acquisition is presented. Acquisition parameters are also discussed and recommendations for future work are given. The main effort of this thesis, however, was to find appropriate processing methods to address some of the different problems present in datasets acquired in the hardrock environment. Comparison of two computer programs for first arrival seismic tomography was performed in order to find the most suitable one for processing crooked line geometries. Three-dimensional pre-stack depth migration was also tested to find a detailed near-surface image. A processing method geared to enhance the signal-to-noise ratio was applied to the dataset with the lowest signal amplitudes to improve the quality of the stack. Finally, cross-dip analysis and corrections were performed on two of the three datasets included in this thesis. Cross-dip analysis was also applied as an interpretation tool to provide the information needed for estimation of the true dip of some of the reflectors related to geological structures. The results presented in this thesis indicate that cross-dip analysis and corrections are one of the most powerful tools for processing and interpretation in the presence of complex geology. Therefore, it is recommended to include this method as a standard step in the processing and interpretation sequence of data acquired in hardrock environments.

The study of deformation zones is of great geological interest since these zones can separate rocks with different characteristics. The geometry of these structures with depth is important for interpreting the geological history of an area. Paper I to III present 2D reflection seismic data over deformation zones targeting structures in the upper 3-4 km of the crust. These seismic profiles were acquired with a crooked-line recording geometry. 2D seismic processing assumes a straight recording geometry. Most seismic processing tools were developed for sub-horizontally layered structures. However, in the crystalline rocks in Scandinavia more complex structures with contrasting dip directions and folding are common. The crooked-line recording geometries have the benefit of sampling a 3D volume. This broader sampling can be used to gain knowledge about the true geometry of subsurface structures. Correlation with geological maps and other geophysical data along with seismic data modeling can be used to differentiate reflections from faults or fracture zones from other reflectivity, e.g. mafic bodies. Fault and fracture zones may have a large impedance contrast to surrounding rocks, while ductile shear zones usually do not. The ductile shear zones can instead be interpreted based on differing reflectivity patterns between domains and correlations with geology or magnetic maps. Paper IV presents 3D reflection seismic data from a quick-clay landslide site in southern Sweden. The area is located in a deformation zone and structures in unconsolidated sediments may have been influenced by faults in the bedrock. The main target layer is located at only 20 m depth, but good surface conditions during acquisition and careful processing enabled a clear seismic image of this shallow layer to be obtained.The research presented in this thesis provides increased knowledge about subsurface structures in four geologically important areas. The unconventional processing methods used are recommended to future researchers working with data from crooked-line recording geometries in crystalline environments. The imaging of shallow structures at the quick-clay landslide site shows that the 3D reflection seismic method can be used as a complement to other geophysical measurements for shallow landslide site investigations.

<p>Seismic methods have been used in sedimentary environment for almost 80 years. During that time, exploration geophysicists have developed a number of techniques to handle specific aspects of working in sedimentary areas. This is not the case for studies in the hardrock environment, where significantly less time and money have been invested on seismic investigations. Therefore, there is still a need to develop the right techniques appropriate for working in hardrock environments. The research presented here, covers aspects of acquisition, processing and interpretation in hardrock environments. A cost-effective alternative for two-dimensional data acquisition is presented. Acquisition parameters are also discussed and recommendations for future work are given. The main effort of this thesis, however, was to find appropriate processing methods to address some of the different problems present in datasets acquired in the hardrock environment. Comparison of two computer programs for first arrival seismic tomography was performed in order to find the most suitable one for processing crooked line geometries. Three-dimensional pre-stack depth migration was also tested to find a detailed near-surface image. A processing method geared to enhance the signal-to-noise ratio was applied to the dataset with the lowest signal amplitudes to improve the quality of the stack. Finally, cross-dip analysis and corrections were performed on two of the three datasets included in this thesis. Cross-dip analysis was also applied as an interpretation tool to provide the information needed for estimation of the true dip of some of the reflectors related to geological structures. The results presented in this thesis indicate that cross-dip analysis and corrections are one of the most powerful tools for processing and interpretation in the presence of complex geology. Therefore, it is recommended to include this method as a standard step in the processing and interpretation sequence of data acquired in hardrock environments. </p>

The South Portuguese Zone (SPZ), which host world-class massive sulphide deposits, forms the southern fold-and-thrust belt of the Iberian Variscan orogeny. This thesis focuses on seismic-reflection and seismic-refraction processing efforts on a subset of the IBERSEIS deep seismic-reflection data set aiming at resolving the SPZ upper crust in high resolution. A comparison of different crooked-line seismic-reflection imaging schemes showed that a processing sequence involving dip-moveout corrections, a common-midpoint projection, and poststack time migration of common-offset gathers provided the most coherent images considering the crooked acquisition geometry. Correlation with surface-geological data allows four units of different reflection character to be identified: the ~0–2 km deep Upper Carboniferous Flysch group, the highly reflective ~2–4 km thick and up to ~5 km deep Volcano-Sedimentary Complex (VSC) group, and two deep Paleozoic metasedimentary units, with the shallower Phyllite-Quartzite group exposed in an antiform. Prominent diffracted energy was enhanced using a modified Kirchhoff imaging routine. High reflectivity and distinct diffractions mark extensive dike bands at 6–12 km depth, possibly related to the intense hydrothermal activity that led to the formation of the ore-bearing VSC group. Source-generated noise obscures potential signals from depths shallower than ~500m depth on the seismic-reflection sections. P- and SV-wave first-arrival traveltimes were inverted for velocity models imaging the shallowest crust. Overall, the velocity models correlate well with surface-geological data marking high (&gt;5.25 km/s) and uniform P-velocities for the Flysch unit in the southern SPZ. A prominent P-wave low-velocity body (~4.5 km/s) is resolved where the Phyllite-Quartzite unit forms the core of an antiform. P-velocities fluctuate the most in the northern SPZ with Flysch group units exhibiting high velocities (&gt;5.25 km/s) and VSC group bodies showing intermediate velocities (~5 km/s). Low VP/VS-ratios (~1.8) computed for the southern profile part are interpreted as less deformed Flysch-group units, whereas high VP/VS-ratios (~1.9) indicate fractured units.