ACTAUNIVERSITATIS

UPSALIENSISUPPSALA

2015

Digital Comprehensive Summaries of Uppsala Dissertationsfrom the Faculty of Science and Technology 1269

Application of the SeismicReflection Method in MineralExploration and Crustal Imaging

Contributions to Hardrock Seismic Imaging

OMID AHMADI

ISSN 1651-6214ISBN 978-91-554-9290-8urn:nbn:se:uu:diva-259396

Dissertation presented at Uppsala University to be publicly examined in Hambergsalen,Geocentrum, Villavägen 16, Uppsala, Friday, 25 September 2015 at 10:00 for the degree ofDoctor of Philosophy. The examination will be conducted in English. Faculty examiner: Dr.Gilles Bellefleur (Geological Survey of Canada).

AbstractAhmadi, O. 2015. Application of the Seismic Reflection Method in Mineral Exploration andCrustal Imaging. Contributions to Hardrock Seismic Imaging. Digital ComprehensiveSummaries of Uppsala Dissertations from the Faculty of Science and Technology 1269. 76 pp.Uppsala: Acta Universitatis Upsaliensis. ISBN 978-91-554-9290-8.

The seismic reflection method has been used extensively in mineral exploration and for imagingcrustal structures within hardrock environments. In this research the seismic reflection methodhas been used and studied to address problems associated with hardrock settings. Papers I andII, address delineating and imaging a sulfide ore body and its surrounding rocks and structuresin Garpenberg, central Sweden, at an active mine. 3D ray-tracing and finite-difference modelingwere performed and the results suggest that although the detection of the ore body by the seismicreflection method is possible in the area, the presence of backfilled stopes in the mine makesseismic imaging of it difficult. In paper III the deeper structures of the Pärvie fault system innorthern Sweden were revealed down to about 8 km through 2D seismic reflection profiling.The resulting images were interpreted using microearthquake data as a constraint. Based on theinterpretation, some locations were suggested for future scientific deep drilling into the faultsystem. In paper IV, the seismic signature of complex geological structures of the Cue-WeldRange area in Western Australia was studied using a portion of a deep 2D seismic reflectionprofile. The pronounced reflections on the seismic images were correlated to their correspondingrock units on an available surface geological map of the study area. 3D constant velocity ray-tracing was performed to constrain the interpretation. Furthermore, the proposed structuralmodel was tested using a 2D acoustic finite-difference seismic modeling method. Based on thisstudy, a new 3D structural model was proposed for the subsurface of the area. These studieshave investigated the capability of the seismic reflection method for imaging crustal structureswithin challenging hardrock and complex geological settings and show some its potential, butalso its limitations.

Keywords: Seismic reflection, Hardrock, Mineral exploration, Crustal imaging, Interpretation,Modeling

Omid Ahmadi, Department of Earth Sciences, Geophysics, Villav. 16, Uppsala University,SE-75236 Uppsala, Sweden.

© Omid Ahmadi 2015

ISSN 1651-6214ISBN 978-91-554-9290-8urn:nbn:se:uu:diva-259396 (http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-259396)

Dedication

To my beloved wife, Maryeh.

List of Papers

This thesis is based on the following papers, which are referred to in the text by their Roman numerals.

I Ahmadi O, Juhlin C, Malehmir A, Munck M (2013) High-resolution

2D seismic imaging and forward modeling of a polymetallic sulfide deposit at Garpenberg, central Sweden. Geophysics 78:B339–B350. doi: 10.1190/geo2013-0098.1

II Ahmadi O, Juhlin C (2015) The effect of the backfilled stopes on seismic imaging of a sulfide deposit in Garpenberg, central Sweden. Manuscript.

III Ahmadi O, Juhlin C, Ask M, Lund B (2015) Revealing the deeper structure of the end-glacial Pärvie fault system in northern Sweden by seismic reflection profiling. Solid Earth 6:621–632. doi: 10.5194/se-6-621-2015

IV Ahmadi O, Koyi H, Juhlin C, Gessner K (2015) Seismic signatures

of complex geological structures in the Cue-Weld Range area, Mur-chison domain, Yilgarn Craton, Western Australia. Submitted to Tectonophysics.

The published papers were reprinted with the permission from the corre-sponding publishers. In addition, below is a list of other publications during my PhD studies, which are not included in the thesis:

Ahmadi O, Juhlin C, Gessner K (2014) New Insights from Seismic Imaging

Over the Youanmi Terrane, Yilgarn Craton,Western Australia. Energy Procedia 59:113–119. doi: 10.1016/j.egypro.2014.10.356

Ahmadi O, Hedin P, Malehmir A, Juhlin C (2013) 3D seismic interpretation and forward modeling – an approach to providing reliable results from 2D seismic data. In: 12th Biennial SGA Meeting. Society for Geology Applied to Mineral Deposits (SGA), Uppsala, Sweden.

Ahmadi O, Juhlin C (2013) Seismic Forward Modeling of a Poly-metallic Massive sulfide Deposit at Garpenberg, central Sweden. In: 75th EAGE Conference & Exhibition incorporating SPE EUROPEC 2013. The Eu-ropean Association of Geoscientists and Engineers (EAGE), London, United Kingdom.

Contents

1.  Introduction ......................................................................................... 11 

2.  Hardrock environments ....................................................................... 13 2.1.  Lappberget sulfide ore body ....................................................... 14 2.2.  The end-glacial Pärvie fault system ........................................... 14 2.3.  Cue-Weld Range greenstone unit ............................................... 15 

3.  Methods ............................................................................................... 19 3.1.  Seismic reflection method .......................................................... 19 

3.1.1.  Seismic data acquisition ......................................................... 20 3.1.2.  Seismic data processing ......................................................... 21 

3.2.  Interpretation and forward modeling .......................................... 25 3.2.1.  Ray-tracing ............................................................................ 26 3.2.2.  Finite-difference method ........................................................ 28 

3.3.  Measurements of acoustic properties ......................................... 31 

4.  Summary of the papers ........................................................................ 34 4.1.  Paper I: Seismic Imaging at Garpenberg .................................... 35 

4.1.1.  Summary ................................................................................ 35 4.1.2.  Conclusions ............................................................................ 39 

4.2.  Paper II: Seismic responses from the Lappberget ore body and backfilled stopes ....................................................................................... 41 

4.2.1.  Summary ................................................................................ 41 4.2.2.  Conclusions ............................................................................ 44 

4.3.  Paper III: Seismic imaging of the Pärvie fault system ............... 46 4.3.1.  Summary ................................................................................ 46 4.3.2.  Conclusions ............................................................................ 52 

4.4.  Paper IV: Seismic signature of a complicated folding system in Western Australia ................................................................................. 53 

4.4.1.  Summary ................................................................................ 53 4.4.2.  Conclusions ............................................................................ 62 

5.  Conclusions ......................................................................................... 63 

6.  Summary in Swedish ........................................................................... 65 

Acknowledgments......................................................................................... 68 

References ..................................................................................................... 70 

Abbreviations

BIF Banded iron-formation CDP Common depth point CMP Common midpoint DMO Dip moveout FD Finite difference f-k Frequency-wavenumber domain Ga Billion years ago Hz Hertz kg Kilogram km Kilometer NMO Normal moveout m Meter Ma Million years ago MHz Mega Hertz ms Millisecond m/s Meter per second VMS Volcanogenic massive sulfide VSP Vertical seismic profiling 2D Two-dimensional 3D Three-dimensional s Second µs Microsecond

11

1. Introduction

Exploration seismology has been established since 1923 and it has been used extensively in natural resources exploration (Sheriff and Geldart 1995; Dentith and Mudge 2014). The basic idea behind it is that seismic waves travel through the earth and may be refracted, reflected or diffracted once they reach an inter-face with different seismic impedance properties compared to their surrounding rocks. The reflected or refracted seismic waves bring us very useful information about the rocks and geological structures in the subsurface. There are generally different methods in exploration seismology which can be implemented in order to address specific problems. One of the most commonly used method is the seismic reflection method (Sheriff and Geldart 1995; Yilmaz 2001) in which reflected or diffracted waves are studied (Telford et al. 1990; Sheriff and Geldart 1995). The geometry of the subsurface structures and their accurate delineation are of great interest to exploration geoscientists. Knowledge about geological structures may have different applications in various projects such as geotechni-cal invistigations, mine-planning, reservoir estimation, oil, gas and mineral ex-ploration (Malehmir et al. 2012a; Telford et al. 1990; Dentith and Mudge 2014). In mineral exploration, the seismic reflection method has had proven capability for almost six decades (e.g. Schmidt 1959; Wu et al. 1995; Eaton et al. 1996; Milkereit et al. 2000; Adam et al. 2003; Bohlen et al. 2003; Malehmir et al. 2012a). The strength of the method lies in its good signal penetration, which makes the method capable of imaging deep targets and structures (e.g. Malehmir et al. 2012a and references therein). This benefit led the mining industry to take advantage of the method and seek solutions to their prospecting problems with the seismic reflection method (Malehmir et al. 2012a). Apart from its applica-tion in natural resource exploration, the seismic reflection method may be used to image crustal structures (e.g. Lie and Andersson 1998; Juhojuntti et al. 2001; Hobbs et al. 2006; Martínez Poyatos et al. 2012). Understanding the geometry of crustal structures is important and may have applications in risk assessment of natural hazards such as earthquakes, studies of geodynamical processes and detecting mineral deposits associated with geological structures (e.g. Williams et al. 1995; Hurich 1996; Larkin et al. 1996; Haugen and Schoenberg 2000; Juhlin et al. 2000; Juhojuntti et al. 2001; Hobbs et al. 2006; Goleby et al. 2006; Ji and Long 2006; Juhlin et al. 2010; Campbell et al. 2010; Juhlin and Lund 2011; Cheraghi et al. 2011; Lundberg et al. 2012; Ehsan et al. 2014; Jackson et al. 2015). In this thesis, three main problems have been addressed using the seismic reflection method. The first problem is delineating a sulfide deposit at an active

12

mine in central Sweden. The second is revealing the deeper structure of an end-glacial fault system in northern Sweden, and the third is an interpretation of a 2D seismic reflection profile through a complex geological setting in Western Aus-tralia. All these case studies are located within crystalline or hardrock environ-ments. Seismic imaging in hardrock settings is challenging and has been ad-dressed by many researchers around the world (e.g. Nelson 1984; Harman et al. 1987; Read 1989; Milkereit et al. 2000; Adam et al. 2003; Bohlen et al. 2003; Nedimovic and West 2003; Salisbury et al. 2003; Li and Eaton 2005; Dehghannejad et al. 2010; Bongajum et al. 2012; Malehmir et al. 2012a; Greenhalgh and Manukyan 2013; Bellefleur et al. 2015; Buske et al. 2015; Koivisto et al. 2015; Manzi et al. 2015). This dissertation is a comprehensive summary of four papers and consists of four chapters. In the first chapter, an introduction and background of the study is given. The second chapter deals with a definition of the hardrock environment. I have tried to summerize all the methods that I have used in chapter 3, including numerical and experimental methods. Finally, in chapter 4, all the papers upon which the thesis is based on, are summarized.

In paper I, the seismic reflection method was used in attempt to delineate and enhance the image of a sulfide ore body and its surrounding rocks and structures in Garpenberg, central Sweden. A high resolution 2D seismic reflection profile was acquired at an active mine site, a challenging and noisy environment. The seismic data were processed along a straight CDP line and an interpretation was made based on the migrated section. 3D ray-tracing modeling was performed to test the interpretation. A few rock samples from the ore body were provided by Boliden Mines and these were used to estimate the P-wave velocity of the sul-fide ore.

In paper II, I expanded the research by including the backfilled stopes in the modeling that had been mined out and refilled during the mining process. Addi-tional rock samples were collected from the mine and their physical properties measured using an ultrasonic instrument. In order to obtain more reliable results and a better approximation to the wave equation, a 3D acoustic finite-difference method was used for modeling.

In paper III, the deep structures of the Pärvie fault system were revealed down to about 8 km via a new 2D seismic profile. Furthermore, an interpretation was made using microearthquake data as a constraint. Based on the interpreta-tion, some locations were suggested for future scientific deep drilling into the fault system.

In paper IV, the seismic signature of the complex geology of the Cue-Weld Range area in Western Australia was studied. Enhanced structures on a 2D seismic section were correlated to their corresponding rock units on a surface geological map. 3D constant velocity ray-tracing was performed to constrain the interpretation. Furthermore, the interpretation was tested using a 2D acoustic FD method.

13

2. Hardrock environments

Rocks can be divided into three main groups: i) sedimentary rocks, ii) igne-ous rocks and iii) metamorphic rocks. Hardrock is a common term that is used for both igneous and metamorphic rocks while for sedimentary rocks, softrock is frequently used (Hall 1988). In general, hardrocks (e.g. granite, gneiss, quartzite) may be harder than sedimentary rocks (e.g. limestone, shale). Although this is true in the broad sense, there are some cases, which do not follow this rule. For instance, marble and limestone are classified as a hardrock and softrock respectively. But the calcite in limestone is as hard as that in marble. Moreover, the quartz grains in sandstone, which is a sedima-trary rock, are as hard as those in quartzite, a metamorphic rock. In addition, the quartz in sandstone which is classified as a softrock is harder than the calcite mineral in marble or even the feldspar or mica grains in granitic rocks (Hall 1988).

Hardrocks and softrocks may be found together in one place, but their origins will be quite different. The formation of hardrocks is controlled by magma flow, deformations, metamorphic and chemical processes in their localities while the deposition of sedimentary rocks is governed by surface or near-surface processes (such as precipitation, morphological processes, weathering etc.). Hardrock environments are often an assembly of complex structures which are formed during ductile and/or brittle deformations and magma emplacement (e.g. Van Kranendonk and Ivanic 2009). Hardrocks often have irregular spatial extent and appear in rough shapes, unlike the sedimentary rocks which commonly have regular dip and bedding. Addi-tionaly, hardrock terrains are also characterized with more discontinuities and greater heterogeneities (e.g. L’Heureux et al. 2009).

A variety of mineral deposits are hosted within hardrock settings. For ex-ample, volcanogenic massive sulfide deposits (VMS) are associated with extrusive igneous processes within intra- and inter-continental formation. Both supergene and hypogene processes may form iron ore bodies associ-ated with banded iron-formations (BIF) and gabbroic rocks (Robb 2005; Allen et al. 2008; Duuring and Hagemann 2013a). There are many other deposits which can be associated with hardrocks (e.g. Au, Cu-Ni, PGE, U and etc.).

In this thesis, three hardrock terrains will be discussed; i) the Lappberget sulfide ore body in central Sweden, ii) the end-glacial Pärvie fault system

14

within Paleoproterozoic greenstones in northern Sweden, and iii) the Cue-Weld Range, a neoarchean greenstone unit in Western Australia.

2.1. Lappberget sulfide ore body The Lappberget sulfide ore is an ore body within the Garpenberg deposit which is located in one of the most important mining regions of Sweden, known as Bergslagen (Stephens et al. 2009) (see Fig. 2.1). The region is well-known for hosting over 6000 mineral occurrences and has a long histo-ry of mining. Bergslagen formed during the Svecokarelian orogeny and is a part of the Fennoscandian Shield (Jansson 2011). Several stages of ductile-brittle deformation have affected the rocks in the region. The majority of the mineral deposits, including the Garpenberg deposit, are found within supracrustal volcanic inliers (Allen et al. 2008). The sulfide mineralogy of the Lappberget ore body consists of sphalerite (ZnS), galena (PbS), pyrite (FeS2) and lesser pyrrhotite (FeS), tetrahedrite ((Cu,Fe)12Sb4S13), argentite (Ag2S) and alabandite (MnS). The ore body is located inside a domal struc-ture in the central parts of the Garpenberg Syncline within an altered dacitic-rhyolitic rock unit (volcanic ash-siltstone). Figure 2.1 shows a geological map of the Garpenberg area.

2.2. The end-glacial Pärvie fault system The Pärvie fault is located in the intraplate Baltic Shield, mostly composed of Paleoproterozoic rocks on the boundary to the Archean terrain of north-ernmost Sweden. The general strike of these rocks is parallel with, and only some 10-50 km away from the Ordovician Caledonian mountain front. Rock types range from acidic volcanic to ultrabasic intrusive rocks. Granitic rocks are extensively exposed at the surface (Fig. 2.2). In some places, sedimen-tary rocks such as sandstone and limestone crop out. The establishment of permanent seismic stations in northern Sweden in 2003-2004, as part of the Swedish National Seismic Network (Bödvarsson and Lund 2003), showed that the Pärvie fault, and the other end-glacial faults in northernmost Swe-den, are still seismically active. Lindblom et al. (2015) analyzed the earth-quakes recorded along the Pärvie fault between 2003 and 2013. Their study showed that the earthquakes occur all along the mapped fault scarp, mostly on the eastern side of the fault in accordance with an east-dipping reverse fault (Lindblom et al. 2015). At four locations along the main Pärvie fault recent investigations on outcrop scale provide the kinematics of brittle de-formation structures and identification of four paleostress fields. Three of the paleostress fields are also found in other major shear zones in the area (Edfelt et al. 2006), whereas the fourth, and most recent field, seems to be

15

unique for the Pärvie fault. This structure indicates a reverse stress field with a strike-slip component. Figure 2.2 shows a geological map of the Pärvie area.

2.3. Cue-Weld Range greenstone unit The Weld Range area is located about 60 km NNW of the town of Cue in the northwestern part of Murchison domain which forms a granite-greenstone unit of the Youanmi terrane within the neoarchean Yilgarn carton in Western Australia (Fig. 2.3). The greenstones in the area have been divided into three major groups (Van Kranendonk and Ivanic 2009): i) mafic volcanic rocks, felsic volcaniclastic sandstones and BIF of the Norie group, 2825-2800 Ma ii) the Polelle group of mafic volcanic rocks, felsic volcanic and volcaniclastic sedimentary rocks and banded iron-formation (BIF) with an age of 2785-2800 Ma and iii) the Glen group of coarse clastic sedimentary rocks and komatiitic basalt, 2724-2700 Ma (Figure 2.3). Voluminous granit-ic rocks were emplaced extensively from 2716 to 2592 Ma. Supercrustal rocks in the area have been metamorphosed to lower greenschist facies (Duuring et al. 2012). The two deformation phases that control geological relationships in the area are: 1) the D1 deformation, which produced the E-trending folds and occurred at the same time as granitic rock emplacement at about 2676 Ma, 2) the regional D2 shortening in the east-west direction which formed the N- to NNE-trending folds from 2660 to 2637 Ma. The second deformation (D2) produced foliation and lineation in the supercustal rocks as well as NNW-striking sinistral shear zones and NE-striking dextral shear zones in the area (Duuring and Hagemann 2013a). The main iron de-posits are hosted within the gabbroic complex of the Weld Range area and there is a high demand to mine them out.

F

igur

e 2.

1. B

edro

ck g

eolo

gica

l map

of

the

Gar

penb

erg

area

(m

odif

ied

afte

r A

llen

et a

l. 20

03).

The

inse

t map

sho

ws

the

loca

tion

of th

e ar

ea in

S

wed

en. T

he d

ashe

d li

ne s

how

s th

e se

ism

ic p

rofi

le o

ver

the

Lap

pber

get o

re b

ody.

F

igur

e 2.

2. B

edro

ck g

eolo

gy o

f th

e P

ärvi

e ar

ea in

nor

ther

n S

wed

en (

© L

antm

äter

iet I

2014

/006

01).

The

inse

t map

sho

ws

the

loca

tion

of th

e ar

ea

in S

wed

en. T

he b

lack

line

sho

ws

the

loca

tion

of th

e se

ism

ic p

rofi

le.

F

igur

e 2.

3. G

eolo

gica

l map

of

the

Cue

-Wel

d R

ange

are

a (I

vani

c et

al.

2013

b). T

he in

set s

how

s th

e lo

catio

n of

the

area

in W

este

rn A

ustr

alia

.

19

3. Methods

In order to address the problems in all the case studies I have been involved in I used the 2D seismic reflection method. The success of the method in hardrock enviroments has been proven by many reported studies around the world (e.g. Wu et al. 1995; Wu et al. 1995; Eaton et al. 1996; Milkereit et al. 2000; Juhlin et al. 2000; Adam et al. 2003; Bohlen et al. 2003; Nedimovic and West 2003; Dehghannejad et al. 2010; Manzi et al. 2012; Harrison and Urosevic 2012; Malinowski et al. 2012; White et al. 2012; Bellefleur et al. 2012; Malehmir et al. 2013; Bellefleur et al. 2015; Cheraghi et al. 2015). I processed three seismic data sets and imaged subsurface structures of each area along corresponding 2D seismic profiles. In the Lappberget case study, some samples from the target ore and surrounding host rock were collected. This allowed me to measure and obtain physical properties of the samples using an ultrasonic system. In the first paper, in order to provide more con-straints for the interpretation of the seismic data, 3D ray-tracing was per-formed. Later, to improve the modeling, 3D finite difference modeling was carried out in the second paper. In the third paper, a seismic reflection profile was processed to delineate the geometry of the Pärvie fault system. In the fourth paper, I performed 2D acoustic finite-difference modeling and a 3D constant velocity ray-tracing. In the following sections, I briefly explain each method that I used throughout my research.

3.1. Seismic reflection method The seismic reflection method is sensitive to the seismic impedance contrast of the rocks through which the waves propagate. In general, a reflection form a structure may be observed if there is a significant impedance contrast between the structure and the surrounding media. In addition, the signal-to-noise ratio of the seismic data and the fold are other important factors in seismic imaging. Normally, the seismic reflection method is carried out in two phases; seismic data acquisition, and then seismic data processing. These are discussed in the following sections.

20

3.1.1. Seismic data acquisition Enhancing a clear reflection originating from a geological structure, requires sorting source gathers into a set of gathers in which their seismic traces have a common midpoint (CMP). These gathers are named CMP or CDP (com-mon depth point) gathers and the number of traces within each CDP gather equals the fold of that gather. In order to achieve enough fold-coverage for imaging, each reflecting point must be recorded more than once (Sheriff and Geldart 1995). In CMP or CDP recording a seismic spread should move continuously down the seismic line to record and obtain the desired fold-coverage. There are a number of seismic spread types that can be used (Sher-iff and Geldart 1995). Out of these, three of them were employed to acquire the data used in the papers: i) fixed-spread ii) split-spread and iii) end-on spread.

In the fixed-spread the receiver line is fixed and the source moves along the seismic profile. This type of spread is suitable when the length of the profile is relatively short. Moving the source results in a good fold coverage. In paper I, an example of a fixed-spread is demonstrated.

Split-spread is a common survey type where the source is located in the center of the receiver line (Fig. 3.1a). Receivers are spaced regularly depend-ing upon the desired lateral resolution. For example, in a deep seismic sur-vey, 40 m group geophone spacing can be used (see paper IV). Another common seismic survey type is the end-on spread. In the end-on configura-tion, the source is located at one end of the receiver line (Sheriff and Geldart 1995) and receivers are spaced regularly (Fig. 3.1b).

Figure 3.1. Two types of seismic spreads: a) split-spread and b)

end-on spread.

Seismic data acquisition may be acquired along a straight or crooked line and, depending on which, the seismic data may be processed in different

21

ways. The apparent dip of the structures may change if the direction of the seismic line changes. Therefore to avoid ambiguity during interpretation, there should be an effort to design straight seismic lines as much as possible (Sheriff and Geldart 1995). However, in practice this is not always possible because of access and structural complications, especially in hardrock envi-ronments which are usually ruggy and difficult terrains. Therefore, most seismic data are acquired along crooked lines.

3.1.2. Seismic data processing In the seismic reflection method, we aim to image subsurface structures or singularities. In the conventional seismic data processing, the imaged struc-tures in the final seismic section should be migrated to their true locations at depth. Most migration algorithms are based on the single-diffracting as-sumption, while in fact seismic waves are multiple diffracted in the presence of random noise in the medium when propagating through the heterogeneous subsurface (Fleury 2012; Broggini and Snieder 2012). This implies that seismic data should be processed to provide the data ideally free of noise, multiples and unwanted wave modes (e.g. S-waves, ground roll and surface waves) before migration. In general, conventional seismic data proceesing can be divided into three main steps: 1) deconvolution 2) CDP stacking and 3) migration (Yilmaz 2001). In addition, there are some other common proc-essing operations that should be applied before the main procedures. Pre-processing, noise attenuation and filtering, static corrections and cross-dip correction are some of these. Before explaining the main processing steps, the preprocessing operations are briefly summarized.

Pre-processing is required to convert the format of the seismic data. Usu-ally, the data are acquired in SEGD format and in order to process it, the data need to be converted to SEGY format, which is a standard format for seismic data processing (Sheriff and Geldart 1995; Dentith and Mudge 2014). The geometry of the seismic line should be set up. Accurate coordinates of source and receivers must be defined since these play an important role in static corrections and stacking of the data correctly.

Filtering is an essential step in seismic data processing. In order to re-move or attenuate unwanted waveforms such as ground roll and airwaves, an appropriate filter can be designed and applied to the data. For example a high-pass filter may be applied to remove low-frequency surface waves (Dentith and Mudge 2014). Dip-filtering (f-k filtering) can also be used if carefully designed. Sometimes, because of the complexity of the subsurface, it is not straightforward to obtain a clean shot gather with only one filter. Therefore, we may need to design and apply a time-varing filter and enhance reflections appropriately at different arrival times.

22

Static corrections have a great effect on boosting refletions. These correc-tions are divided into two groups: i) refraction static corrections ii) reflection (residual) corrections (Yilmaz 2001). The refraction static corrections try to correct for near surface low velocities. These low velocities correspond to weathered or the loose sediments near the surface and may cause local varia-tions in arrival times (Dentith and Mudge 2014). Therefore, we need to ap-ply a time shift to restore the seismic traces as if they were recorded in an ideal terrain with constant topographic elevation (Dentith and Mudge 2014). Furthermore, we may apply reflection (residual) static corrections where a time shift based on the lateral coherceny of the reflections is calculated and applied to increase coherency of the reflections in the seismic data (Dentith and Mudge 2014).

The cross-dip correction is necessary to correct for out-of-plane compo-nents or reflections in the seismic data (Sheriff and Geldart 1995). This cor-rection is more relevant when the seismic line is crooked (see Figure 3.2) (Larner et al. 1979; Nedimovic and West 2003; Rodriguez-Tablante et al. 2007). The procedure consists of applying calculated cross-dip corrections on normal moveout (NMO)-corrected CDP-binned gathers and then stacking the CDP gathers. Equation (3.1) expresses the calculation of the cross-dip correction where dx is the perpendicular distance of a midpoint from the stacking line:

. (3.1)

23

Figure 3.2. A crooked seismic line (black line) and a layer with a cross-dip compo-nent. Cross-dip corrections can be calculated by equation (3.1) for this layer. The blue line depicts the CDP stacking line.

3.2.1.1. Deconvolution Ideally, we would like to extract the earth’s reflectivity from seismograms. In order to do this, deconvolution can be applied. Deconvolution tries to sharpen the wavelet using an inverse filter and thus increase temporal resolu-tion of the data (Yilmaz 2001). In theory, convolving the deconvolution op-erator with the seismic trace results in a series of spikes which represent the subsurface reflectivity. In this concept, it is assumed that the seismic wavelet is a minimum-phase wavelet and therefore the source wavelet is already known. This type of deconvolution is called spiking deconvoltion (Yilmaz 2001). But in practice, because of the complexity of the surface, the source wavelet may vary during the seismic recording. In this case, predictive de-convoltion can be used. In this type of deconvolution, an effort is needed to to design a reasonable inverse operator (Yilmaz 2001) by changing different parameters (such as gap, length and random noise level).

3.2.1.2. CDP stacking Stacking along a CDP line is another complimentary procedure which will normally attenuate noise and multiples (Yilmaz 2001). In order to design a CDP line, coordinates of receivers and sources must be defined accurately. A desired CDP line needs to be defined depending upon the geometry of the seismic profile (Fig. 3.3). If the line is crooked, then it is best to split the CDP line into a number of straight segments which run within mid-points

24

(Wu 1996). CDP bins are usually designed perpendicular to the CDP line, but can also be at some angle to the stacking line where the CDP bins are parallel to strike of the geological structures (Kashubin and Juhlin 2010; Lundberg and Juhlin 2011). The seismic data are then sorted in CDP gathers.

Figure 3.3. Straight CDP line (red line) along the 2D cooked seismic profile (green line) in Garpenberg, Sweden (Paper I). The black cloud depicts mid-points of the seismic data.

Since conventional seismic imaging and migration is based on zero-offset traces, it is necessary to correct the seismic traces to as if they had been ac-quired with the sources and receives at the same location (zero-offset). The normal moveout correction (NMO) is the simplest correction that can be made to obtain the zero-offset traces (Yilmaz 2001). A drawback of the NMO correction is that it does not take into account confilicting dips in the subsurface and therefore a bias arises in achieving good zero-offset gathers (Yilmaz 2001). Additionaly, the moveout of steeply-dipping layers is in con-flict with the moveout associated with gently-dipping structures (Yilmaz 2001). In order to enhance possible dipping layers, dip moveout (DMO) is another pre-stack process which may be applied (Hale 1984; Deregowski 1991). After applying DMO, the NMO correction should follow and the

25

resulting zero-offset CDP gathers be stacked to provide a final stacked sec-tion in the time domain.

3.2.1.3. Migration Migration is the final process applied to re-position reflections to their true geological positions at depth (Yilmaz 2001). Migration also collapses dif-fractions and improves the quality of the final image (Yilmaz 2001). In con-ventional seismic data processing, migration is applied after stacking and is called post-stack migration, while it may be applied on pre-stack data as well (pre-stack migration). The accuracy of the final migrated section highly de-pends on the input velocity model (Jones 2010). Hence, a great effort should be made to obtain a good velocity model, which is a major challenge in complex and hardrock environments. Based on the type of algorithm, migra-tion methods can be divided into two major groups: 1) ray-based migration methods and 2) differential equation-based migration methods (Jones 2010). Kirchhoff migration is a common ray-based (integral) method which can be applied even when steeply-diping layers are present (Yilmaz 2001). A draw-back of this method is the difficulty to handle lateral velocity variations (Yilmaz 2001). Differential equation-based methods can be divided into two subgroups: 1) finite-difference (FD) based methods (e.g. downward con-tinuation or FD method) and 2) frequency-wave number based methods (e.g. Stolt method) (Yilmaz 2001). The FD method handles well any type of sub-surface velocity model, but one needs to be very careful on choosing the right dip approximation for the method (Yilmaz 2001). In contrast, the fre-quency-wave number methods are limited in handling lateral velocity varia-tions (Yilmaz 2001). Therefore, each migration method has its own advan-tages and drawbacks and the selection of the method to be used is case de-pendent.

3.2. Interpretation and forward modeling Final migrated seismic sections should be interpreted and all imaged struc-tures should, when possible, be correlated to mapped rock units on surface geological maps or stratigaphical columns extracted from well logs. Occa-sionally, there is a lack of information and it is not possible to interpret re-flections with confidence. However, even with few constraints, it is still fea-sible to provide an explanation for prominent reflections. Forward modeling has an important role to evaluate and test interpretations. A synthetic data can be produced with the same survey configuration as the real data. The real and synthetic data can then be compared (e.g. Blake et al. 1999; Bohlen et al. 2003; L’Heureux et al. 2009; Cheraghi et al. 2011; Malinowski et al. 2012; Bellefleur et al. 2012; Squelch et al. 2013). In practice, interpreted structures and their corresponding physical properties can be discretized on

26

1D, 2D or 3D grids. Furthermore, a numerical experiment can be performed throughout the grid and simulate wave propagation with the same survey as the real layout. There are many methods to model seismic wave propagation. Ray-tracing and finite-difference methods are the most common techniques that I have also used in the course of my research.

3.2.1. Ray-tracing Asymptotic ray-tracing is a high frequency based technique to simulate seismic rays and calculate arrival times for the modeled reflections. In order to avoid the scattering problem with long wavelength seismic waves, we need to use the high frequency approximation to the wave equation (Gjøystdal et al. 2002; Chapman 2004). By employing the ray-tracing tech-nique, we aim to calculate traveltimes to reflectors in a model. The ampli-tude of the reflections can also be calculated (Cerveny 2005). The essential governing equations of the method for acoustic media are briefly presented here. Acoustic wave propagation through a smoothed inhomogeneous media can be represented with the following equation (Cerveny 2005):

, (3.2)

where ϕ is the pressure wavefield, ϕtt is the second derivative of the pres-

sure field with respect to time and c(xi) is variable velocity function with constant density assumed throughout the medium. The source term is not shown in the equation. The following time-harmonic high frequency solution satisfies the above equation:

, . (3.3) Here is the high frequency term, while P(xi) and T(xi) are assumed to be

smooth scalar functions of coordinates. Equation (3.3) forms a plane wave solution and plugging it into the equation (3.2) yields:

1

2 .

0. (3.4) In order for equation (3.4) to be valid, we arrive at the Eikonal equation:

, (3.5)

27

and the transport equation:

2 . 0. (3.6) Equations (3.5) and (3.6) are very important and should be solved in the

ray-tracing (Cerveny 2005). The last term∇ 2P in equation (3.4) also contains the amplitude of the waves and needs to be considered in some cases (Scales 1997). This amplitude term, may be written as a function in a series of in-verse powers of frequency and then can be solved for (Cerveny 2005). Fig-ure 3.4 shows the Marmousi model in which a source at 4 km distance is excited. The corresponding seismic rays are traced and overlaid on the model. Ray-trace modeling has been widely used around the world (e.g. Julian and Gubbins 1977; Cerveny 1989; Farra 1993; Blake et al. 1999; Gjøystdal et al. 2002). The governing equations can be expanded for more complex and heterogeneous media. In addition, elastic waves can also be considered (Cerveny 2005).

Figure 3.4. Ray-tracing through the Marmousi model. A seismic source is located at 4 km distance and at the top of the model (sea surface). The ray-tracing was per-formed using the CREWES toolbox in Matlab (Margrave 2000).

A simplified type of ray-tracing is 3D constant velocity ray-tracing, which was used in paper IV. This type of ray-tracing is suitable to model reflec-tions in pre-stack data and then track them on an unmigrated stacked section (Ayarza et al. 2000; Lundberg et al. 2012). The method is based on Fermat’s principle for traveltime calculations and the amplitude of the reflections can be calculated using the formulas given by Aki and Richards (2009). Using constant velocity ray-tracing, traveltime curves of planar reflections are cal-culated (Fig. 3.5). Dip and azimuth of a planar reflector are changed with an iterative approach in small increments until the best fit is achieved. Use of this method should result in the calculation of reasonable traveltimes for the reflections when velocity variations are small (Lundberg et al. 2012).

28

Figure 3.5. Calculated traveltimes using 3D constant velocity ray-tracing on a shot gather of the 10GA-YU1 data set in the Cue-Weld Range area, Western Australia: a) two reflections are prominent in the shot gather where reflection A1 appears to come to the surface while reflection A2 does not; b) travel times to a plane with a strike of N45°E and a dip of 30° NW match well with both reflections A1 and A2 and c) travel times of the modeled reflections on the unmigrated stacked section.

3.2.2. Finite-difference method The finite-difference method is the most frequently used method for model-ing of seismic wave propagation (Fichtner 2011). This method provides a good approximation to the solution of the wave equation and it has been used for seismic modeling through various hardrock environments (e.g.

29

Yerneni et al. 2002; Bohlen et al. 2003; Saenger and Bohlen 2004; L’Heureux et al. 2009; Cheraghi et al. 2011; Bellefleur et al. 2012; Squelch et al. 2013). The wave equation has two general modes; the acoustic wave and elastic wave equations. The earth is an imperfect elastic media and the acoustic wave equation is only truly valid in a fluid media where the shear moduls equals zero (Aki and Richards 2009; Fichtner 2011). Therefore, in order to simulate wave propagation through the subsurface, it is preferable to use visco-elastic or elastic wave equations (Juhlin 1995; Saenger and Bohlen 2004). However, the acoustic wave equation is computationally less expen-sive than the elastic wave equation (Fichtner 2011). Using the acoustic mode of the wave equation is strictly only justifiable when the data processing is limited to the first arrivals of P-waves and also when an active seismic source (such as explosive) is used during recording that generates less S-wave energy (Fichtner 2011). However, the acoustic expression of the wave equation has been used extensively in active source seismic modeling (Alkhalifah 2000; Wang and Liu 2007; Sava and Poliannikov 2008; Bean and Martini 2010; Sava and Vlad 2011; Squelch et al. 2013). In this thesis, 2D and 3D acoustic FD methods were used. To understand the methods bet-ter, the mathematics of the 2D FD method is demonstrated below.

Equation (3.2) for a 2D acoustic problem may be expressed as follows:

, ,,

, , . (3.7)

Now, the right side of the above equation can be written using the second

order finite-difference approximation in the following form:

, ,,

, , Δ2 , , , , Δ

,(3.8)

and then we have:

, , Δ = 2 Δ , , , , , Δ . (3.9)

Based on equation (3.9), we are able to calculate the seismic wavefield at

t+∆t using two earlier wavefields at t and t-∆t (Margrave 2003) (see Fig. 3.6).

The numerical experiment has to be stable and accurate as much as possi-ble. Spatial sampling (grid spacing) must be chosen adequately to avoid grid dispersion. The grid spacing (h) needs to satisfy the following general condi-tion, in order to have an accurate experiment and be as free of grid disper-sion as possible:

30

, (3.10)

where is the minimum wavelength in the model and n depends on

the order and length of the FD operator (Bohlen et al. 2015). For the 4th-order in space and 2nd-order in time approximation we have (Bohlen et al. 2003):

. (3.11)

If the amplitude of the solution grows without bound during the simula-

tion, then the experiment is not stable. In order to avoid this, time-stepping ∆ should also fulfil the following condition (Fichtner 2011):

Δ , (3.12)

in which, and denote the minimum grid spacing and maximum

velocity in the model. The constant coefficient depends on the method used for space and time discretisation (Fichtner 2011). For example, for a finite-difference algorithm based on the 4th-order and 2nd-order approximation in space and time, respectively, and in a staggered grid, the expression (3.11) can be represented as (Blanch 1996; Bohlen et al. 2003):

Δ√

, (3.13)

with d as the dimension of the model.

Figure 3.6. A snapshot of P-wave propagation through the Marmousi model at 1 s. The source is located at 4 km distance and at the sea surface. The FD modeling was performed using Madagascar software (Fomel et al. 2013).

31

The grid model used to discretize the subsurface physical properties is limited by boundaries. Therefore, during seismic modeling, reflections from the boundaries interfere with responses from structures within the model. This effect needs to be removed or attenuated. One way to do this is by add-ing a number of zero nodes to the boundaries to absorb the seismic energy (Fichtner 2011). Although the grid padding technique is not a perfect solu-tion for the boundary problem, it has been used widely (e.g. Juhlin 1995; Bohlen 2002; Margrave 2003; Fichtner 2011). However, this technique will not attenuate all seismic waves with arbitrary incident angles (Fichtner 2011). Furthermore, using the grid padding, makes the model bigger and therefore computationally more expensive. Fortunately, seismic modeling can be coded in parallel and be performed on cluster or multi-node comput-ers within a reasonable computational time (Bohlen 2002; Weiss and Shragge 2013; Fomel et al. 2013).

3.3. Measurements of acoustic properties In order to carry out seismic modeling using either ray-tracing or FD meth-ods, we need to define the physical properties (velocities and densities) of the rocks as input. One way to do this is to collect some rock samples from the area of interest. The samples should represent the lithology range within the area as much as possible. Seismic velocities of the samples can be meas-ured using ultrasonic instruments (Vanorio et al. 2002; Salisbury et al. 2003; Stanchits et al. 2006; Malehmir et al. 2013). An ultrasonic device consists of two transducers (one receiver and the other source) connected to a pulse generator. The transducers induce a high frequency signal (e.g. 1 MHz) through the rock sample. A signal is passed through an oscilloscope and finally a computer to record the waveform (Fig. 3.7). The procedure can be done both at room condition and/or at elevated confining pressures (Salis-bury et al. 2003; Malehmir et al. 2013). The rock samples should be cut in a way to make flat surfaces on each side of the samples. Sometimes samples are taken from core where it is easy to cut the top and the bottom of the sam-ple. If the sample has an irregular shape, an effort should be made to form a cubic shape out of the sample (Malehmir et al. 2013). After cutting, the sam-ples are ready for the ultrasonic measurements. I performed the measure-ments at atmospheric condition due to the instrumental limitations.

The density of the rock samples should also be measured. Although the rock samples are cut for the ultrasonic measurements, they still appear with ir-regular shapes and the calculation of their densities is not straightforward. There are different methods to obtain the density values such as suspension, level and overflow methods (Hughes 2005).

32

Among these, the suspension method is the most accurate and precise method (Hughes 2005). This method is based on Archimedes’ principle. This principle states that when an object is suspended in a liquid, it displaces the liquid and the weight of the displaced liquid equals the force which buoyes up the object (Halliday et al. 2000). The suspension device consists of an electronic scale and a container with enough water to suspend the rock samples entirely in it (Fig. 3.8). Based on Archimedes’ principle, the volume of the immersed sample equals the difference in weight divided by the den-sity of the fluid (Hughes 2005):

V , (3.14)

Figure 3.7. Ultrasonic measurement: a) an illustration of the ultrasonic apparatus at room conditions; b) a recorded waveform through a sulfide ore sample from the Lappberget ore body using the ultrasonic instrument at room conditions.

33

where V is the volume of the sample, ∆w the difference in weight and ρf is the density of the fluid. If we use pure water as the fluid in the experiment with a density of about 1 (kg/m3), we can write:

V

Δ . (3.15)

Figure 3.8. Density measurement based on Archimedes' principle using the suspen-sion device.

The mass of the sample can be measured using the electronic scale. We

are then able to calculate the density of the sample as:

. (3.16)

34

4. Summary of the papers

In this section, each paper is summarized to provide a brief description about the main goals, methods and results of the studies. My contributions in these studies were as follows:

Paper I: I participated in the data acquisition in 2011. I processed the seis-mic data with advice from my supervisor. Boliden Mines provided the rock samples from the mine and I conducted the lab measurements using the ul-trasonic instrument. I performed the 3D ray-tracing using the Norsar-3D program (www.norsar.no/seismod/Products/NORSAR-3D/). The paper was written initially by me and later improved by my co-authors.

Paper II: I performed the numerical modeling using Madagasacr software (www.ahay.org/wiki/Main_Page). The synthetic data were processed by me. Boliden Mines provided us with additional rock samples and I measured their elastic properties. I wrote the main text of the paper. My supervisor helped with the numerical modeling and improving the paper.

Paper III: I participated in the field work and seismic data acquisition. I processed the data. My supervisor has given me practical advice throughout the data processing. I wrote the processing and interpretation parts of the paper and also synthesized my co-authors’ contributions in the conclusion of the paper. The co-authors had a significant role in improving the initial drafts of the paper.

Paper IV: I reprocessed the real seismic data and processed the synthetic data. I owe the core of the interpretation to my co-authors. I built the 3D geological model based on the interpretation and performed the numerical modeling. I wrote the main parts of the paper.

35

4.1. Paper I: Seismic Imaging at Garpenberg

4.1.1. Summary Garpenberg is a well-known mining site in central Sweden because of its high potential for sulfide ores. The Lappberget ore body is the largest known ore reserve of the Garpenberg area (see Fig. 2.1). In May 2011, a 2D seismic profile was acquired over Lappberget in an attempt to delineate some of the lithological boundaries and, perhaps, image the ore body (Table 4.1). The seismic data were processed along a straight CDP line using a conventional processing method (Table 4.2). Since the seismic profile was designed to cross over the Lappberget ore body, the data suffered from noise due to the mining activity and traffic. Therefore, an attempt was made to attenuate the noise and enhance reflections. Following a refraction static correction, an f-k filter boosted a reflection in the data (Fig. 4.1). Furthermore, a cross-dip correction with a dip of -15° resulted in some pronounced reflections appear-ing on the stacked section (Fig. 4.2). The DMO correction also improved the image significantly (Fig. 4.3).

Table 4.1. Main seismic data acquisition parameters used for the Garpenberg seismic test profile.

Date May, 2011

Geometry 2D crooked line

Sampling rate 1 ms

Recording length 30 seconds – 4 sweeps

Geophones 28 Hz

Receiver spacing 5 m

Shot point spacing 5-10 m

Source VIBSIST

Number of active channels 394

36

Figure 4.1. A raw shot gather from the Garpenberg data and the corresponding im-provement after processing: (a) the raw shot gather; the ellipse encompasses the strong noise generated from infrastructures at the mining site; (b) the shot gather after static corrections and spherical divergence correction; the arrow marked by A shows P-wave first arrivals alignment before and after the correction, shear waves are also marked by B; (c) the shot gather after applying dip filtering and a median filter; the noise and shear waves are attenuated; (d) the shot gather after applying spectral whitening and deconvolution; airwaves are removed, shear waves are fur-ther attenuated, and the arrows show a reflection that is enhanced.

The reflections were interpreted as corresponding to geological rock units in the area that partly interfere with the potential ore body signal. To further investigate the seismic response of the ore body, 3D ray-tracing was applied using the ore body geometry. The ore body model was based on in-mine and drill core geological observations provided by Boliden Mines. A few rock samples from the ore body and its host rock were also collected. An ultra-sonic instrument was used to measure the physical properties of the rocks at

37

room conditions. The measurements showed a considerable contrast between the ore body and host rock. We assumed two units to be present in the mod-el, the ore body (with its internal structure) and the host rock.

Table 4.2. Data processing steps applied to the Garpenberg data.

1. Decoding VIBSIST data

2. Preprocessing

3. Spherical divergence correction + scaling 1 dB/s

4. Refraction static correction

5. F-K filtering (Dip ranges 1.5-3.5 ms/trace)

6. Spectral equalization (30-50-200-250 Hz)

7. Deconvolution (length 170 ms – gap 5 ms – added noise 0.1%)

8. Band-pass filtering (30-50-200-250 Hz)

9. Cross-dip correction ( Dip -15 degrees)

10. Residual ctatics correction

11. Velocity analysis

12. NMO and DMO corrections

13. Velocity analysis

14. NMO correction

15. Muting

16. Stacking

17. Band-pass filtering (50-60-110-120 Hz)

18. FX-decon

19. Migration

By comparing the modeled and observed data, we found that the high ampli-tude signal in the real seismic section most likely emanates from near the top of one concentrated ore, which lies inside the larger mapped ore body (see Figs. 4.3 and 4.4). The base of the ore body was only observed in the syn-thetic data (Fig. 4.3c). In order to investigate the effectiveness of the seismic signal, a signal penetration analysis was applied. The analysis suggested that the seismic signal penetrated efficiently along the entire survey line, a proxy to determine the penetration depth. The presence of disseminated ore and lower fold towards the northern end of the profile (see Fig. 3.3) could be combined reasons that made imaging the base of the ore body difficult.

38

Figure 4.2. Results of the cross-dip analysis for a few dip angles. Best results were obtained by using a cross-dip of -15°. The arrow shows the location of the ore body on the profile. Note that the cross-dip 0° section corresponds to the stacked section without applying a cross-dip correction.

Figure 4.3. Stacked sections. (a) NMO stacked section; (b) the real DMO stacked section (below) with the CDP fold coverage plotted over it; (c) the synthetic DMO stacked section. The four main reflection groups which are enhanced are marked on the real DMO stacked section.

39

Figure 4.4. Migration results. (a) The real migrated section; (b) The synthetic mi-grated section. The projected border of the ore body is shown on the sections (view from the east).

4.1.2. Conclusions A 2D high-resolution seismic profile was acquired in the Garpenberg area and it was processed with mainly standard methods after applying a cross-dip correction of -15°. The processing strategies were successful in eliminat-ing noise and boosting signal in the data. A bright-spot anomaly is interpret-ed to originate from near the top of the ore body with the location of this anomaly consistent with the working geological model of the Lappberget ore body. Furthermore, DMO processing improved the image quality of the sub-surface. In order to constrain the seismic interpretations, 3D ray-based mod-

40

eling was carried out using the true acquisition geometry of the seismic pro-file and the geological model of the ore body based on the mine data and drilling. Given the geometry of the ore body, seismic modeling results sug-gested that imaging the top and bottom of the ore body is possible. Migration of the synthetic data and comparison with the real data indicate that the bright reflection does not come from the top of the modeled ore body, but rather from within the ore volume. This study demonstrates that imaging of the sulfide deposits at the site is possible, but their delineation is challenged with only one 2D survey and the absence of additional geological con-straints.

41

4.2. Paper II: Seismic responses from the Lappberget ore body and backfilled stopes

4.2.1. Summary The Lappberget sulfide ore body has been mined out using the transverse sublevel stoping method (Hamrin 1980; Eklind et al. 2007). The seismic response of the ore body was previously studied with 2D seismic data and 3D ray tracing. Ahmadi et al. (2013) enhanced a bright diffraction, or possi-bly a reflection, which was interpreted as originating within the ore volume near to the top of the ore body. In this research, we study the effect of the backfiled stopes on seismic imaging of the Lappberget ore body. The back-filled stopes were backfilled with a mixture of waste rocks and concrete. We perform 3D acoustic finite-difference (FD) modeling for two scenarios; in one case the ore body with the stopes are included and in the other only the ore model is considered. Figure 4.5 shows a cross section of the ore body with the mined-out stopes at the time of the seismic survey in early 2011.

Figure 4.5. A cross section of the Lappberget ore body (view from the east). Beige colored stopes are mined out and included in our modeling. The red line depicts the boundary of the ore body and numbers correspond to the level numbers in the mine (based on Boliden Mines information).

42

Based on the average velocity of the first arrivals of the seismic data and their ultrasonic measurements, Ahmadi et al. (2013) assumed a velocity of 5500 m/s and 4606 m/s for the host rock and the ore, respectively, in the modeling. Here, we have also used these parameters for the modeling. In addition, a P-wave velocity of 3000 m/s was assumed for the paste-fill mate-rial within the stopes (Saleh and Milkereit 2014). The model has a total size of 1330×2240 ×2008 m3 with a cell size of 5×5×6.5 m3. A Ricker wavelet with a frequency of 100 Hz was used in the modeling. In order to absorb and attenuate reflected waves from the boundaries, the sides of the model were padded with 60 nodes, except for the top. Here, the free surface boundary condition was invoked. Source gathers were only generated for every other source point, requiring about two days to produce 200 source gathers on a standard desktop computer (Fig. 4.6). The synthetic data were processed in the same way as the real data (see Table 4.2). In this study, we have re-migrated the stacked sections using the Stolt migration method. This im-proved the final seismic image compared to the one that Ahmadi et al. (2013) provided (see Fig. 4.4a). Comparison of the generated synthetic source gathers with the real ones (Fig. 4.6) show that the top and bottom of the ore body are clearly observed in the model without the stopes.

Figure 4.6. Comparison of real and synthetic shot gathers. a) the real raw shot gather which suffers from infrastructure related noise; b) the processed shot gather, the ellipse shows an enhanced reflection or diffraction; c) the synthetic gather without the stopes in the model; d) the synthetic gather with the stopes included.

43

In contrast, the top of the ore body is not very visible in the gather of the model with the stopes included (Fig. 4.6d). The bottom of the ore body is entirely obscured (Fig. 4.6d). The response from the stopes dominates the source gather and the ore body response is masked and not readily detecta-ble. In the synthetic migrated section with the stopes included the top of the ore body can be barely identified while, in contrast, the response of the stopes is very strong (Fig. 4.7b). The migrated diffractions in the model without the stopes collapse to a scatterer which corresponds to the top of the ore body (Fig. 4.7c). Figure 4.8 shows both the synthetic migrated sections with the stopes and the real migrated section, with a vertical slice through the model overlaid on them. The location of the bright anomaly (reflections A) is close to the top of the ore body in the real data, but the modeling re-sults suggest that it is likely that these responses originate from the mine related structures.

Figure 4.7. Migrated and depth converted sections using the Stolt migration method a) the real migrated section; the arrow shows the bright anomaly (reflections A) which was enhanced; b) the synthetic migrated section with the stopes in the model; the arrow shows the recorded response off the top of the ore body which is disguised by the response of the stopes in the model; c) the synthetic migrated section without the stopes.

44

Figure 4.8. Migrated sections and a vertical slice of the velocity model with the stopes overlain: a) the real migrated section; b) the synthetic migrated section. The red arrow shows the weak response off the top of the ore body model.

The diffractions from a depth of about 600 m down to 900 m correspond to the paste-filled stopes in the model (Fig. 4.8b). In the real migrated section, some enhanced scattered responses match nicely with the stopes in the same depth range (Fig. 4.8a). Reflections A appear at a slightly shallower depth of about 580 m on the real migrated section. The response of the top of the ore body is relatively weak compared with the mined-out stopes on the synthetic migrated section (Fig. 4.8b). The modeling shows that seismic characteriza-tion of the ore body is probably difficult due to the strong seismic response of the backfilled stopes in the mine.

4.2.2. Conclusions 3D acoustic modeling was performed to understand the effect of the back-filled stopes on the seismic response from the ore body and the synthetic data were compared with the 2D real data acquired over the ore body in 2011. The modeling results suggest that the ore body response is disguised

45

by the seismic response of the mined-out stopes. The paste-filled stopes gen-erate significant scattering within the ore body and, therefore, the bright diffractions which were enhanced on the real data are likely to correspond to infrastructures in the mine. However, further studies are needed to better understand the real data acquired.

46

4.3. Paper III: Seismic imaging of the Pärvie fault system

4.3.1. Summary The largest fault scarp in Scandinavia belongs to the Pärvie fault system in northern Sweden, both in terms of its length and calculated magnitude of the earthquake that generated it (Lindblom et al. 2015).

Microearthquake activity is currently ongoing in the vicinity of the fault scarps in the area (Fig. 4.9) The relation between the earthquakes and the structures mapped at the surface has been unclear. A new 2D seismic reflec-tion profile for imaging the deeper levels of the Pärvie fault system was ac-quired in June 2014. Important questions concerning the Pärvie end-glacial fault system are i) what is the geometry of the faults at depth and their ex-tent, ii) what was the stress field that caused these ruptures during the last deglaciation, and iii) what is the stress field at depth that is generating the current earthquake activity along the faults.

A seismic profile across the Pärvie Fault system acquired in 2007, with a mechanical hammer as a source, showed a good correlation between the surface mapped faults and moderate to steeply dipping reflections (Juhlin et al. 2010). In the new seismic survey an explosive source with a maximum charge size of 8.34 kg in 20 m deep shot holes was used (Table 4.3) for deeper penetration.

Table 4.3. Seismic acquisition parameters used during the seismic recording.

Profile length 22 km

Seismic source Explosive: 8.34 kg/shot hole Receiver spacing 25 m Nominal shot spacing 400 m Sampling rate 1 ms Recording length 25 s Active channels 300-390

47

Figure 4.9. Location of the seismic profile over the Pärvie fault system in northern Sweden. Earthquakes are shown for three different time periods (after Lindblom et al., 2015). Inset shows the location of the Pärvie area in northern Sweden.

The new seismic profile was laid almost in the same location as the old one. The data were processed using a conventional processing method (Ta-ble 4.4) along a straight CDP line (Figure 4.10). An FX-decon filter with a window of 7 traces resulted in an enhanced reflection which corresponds to the main fault in the area (Figure 4.10 and Figure 4.11).

48

Table 4.4. Seismic data processing steps.

1 Reading SEGD files and sorting the data 2 Geometry setting 3 Picking first arrivals 4 Inner trace muting 5 Refraction static correction 6 Spectral balancing 10-20-100-150 Hz 7 Spherical divergence correction 8 Band pass filtering 10-20-100-150 Hz 9 Reflection residual correction 10 FXDECON 7 traces window 11 Time variant filtering 0-650 ms: High pass 3-70 Hz

655-5000 ms: Band pass 15-25-90-120 12 Velocity analysis 13 NMO correction Constant velocity 13000 m/s 14 Velocity analysis 15 First arrival and ground roll muting 16 Stacking 17 FXDECON 7 traces 18 Band pass filtering 15-25-90-120 Hz 19 Migration

Figure 4.10. The seismic profile location across the Pärvie fault system. The mid-point distribution shows the irregularity of the data. The location where a weak re-flection projects to the surface that was recorded in a few shot gathers (see Fig. 4.11)marked on the map by the blue circle. This location corresponds to the main fault with acquisition line crossing at receiver 105.

49

Figure 4.11. Application of FX-Decon in pre-stack processing; a) before; b) after FX-decon. The arrows show an enhanced reflection which projects to receiver peg of 105 at the surface and corresponds to the main east-dipping Pärvie fault (see the location of the receiver in Fig. 4.10)

Figure 4.12. The NMO stacked section. Reflections R1 and R3 (shown by the ar-rows) correspond to the subsidiary west-dipping fault, and the main Pärvie fault, respectively. Reflection R4 corresponds most likely to another west-dipping fault in the subsurface.

50

Furthermore, a time-varying filter improved the pre-stack data. An exten-sive velocity analysis was carried out and a constant velocity of 13000 m/s was chosen to be used in the NMO correction (Fig. 4.12). The enhanced reflections were migrated using the FD migration method and depth-converted using a variable velocity model obtained from 3D local earth-quake tomography (Lindblom et al. 2015). Figure 4.13 shows the migrated and depth-converted section. The main Pärvie fault (reflection R3) and the west-dipping reflection R1 in Figures. 4.12 and 4.13 have different reflec-tivity characteristics. The main fault is a weakly reflective structure, while the subsidiary west-dipping fault indicates the presence of a strong imped-ance contrast.

Figure 4.13. Cross-section of the microearthquakes (red dots) overlaid on the inter-preted migrated section. Red dashed lines depict the interpreted faults. Fault R2 was not imaged in the new seismic survey. All events within a maximum offset of 10 km perpendicular to the seismic profile to the north and south are projected onto the CDP plane. A dense cluster of seismic activity is located at a depth of about 11.5 km. The listric appearance of the R1 and R3 reflections at depth is most likely due to the migration process and not a feature of the faults. We have interpreted the faults to be planar based on the unmigrated sections.

The difference between the thicknesses of the fracture zones along the faults is a possible explanation for the difference between the reflectivity of the faults. However, if secondary mineralization has occurred along the faults, then the fault R1 has more capacity to host the secondary minerals than the apparently thinner R3 fault. The three steeply dipping faults are

51

imaged down to about 7.5 km in the new dataset while Juhlin et al. (2010) showed four faults in their seismic images down to about 2-3 km. The R2 fault from Juhlin et al. (2010) has not been imaged in our seismic processing. Well-constrained microearthquakes within an area of 400 km2, and with a maximum offset of 10 km at right angles to the seismic line, were projected onto the migrated section (Fig. 4.13). A denser cluster of seismic events is located below the imaged faults, suggesting that the intersection at 11.5 km depth of the main Pärvie fault and the west-dipping R1 fault is the source of the seismic activity in the area. The geometry of the R1 and R3 faults is un-certain below this depth. The main fault (R3) may continue at a similar dip to the east, as indicated more clearly by seismicity on the northern and southern segments of the fault (Lindblom et al. 2015), or it may become listric at deeper levels (Fig. 4.14).

Figure 4.14. A simplified model of the Pärvie fault system. A larger area compared with Figure 4.13 with projected earthquakes (red dots) and the location of some possible boreholes is shown. Below 8 km depth, the main fault (R3) may continue at a similar dip (shown by the red dashed line), become listric (marked by the dotted line) or form the lower part of a flower structure (marked by the black dashed lines). During the first phase of the drilling project (DAFNE), boreholes DA-1-1 and DA-1-2 will penetrate the R1 and R3 faults at about 1 km depth respectively. Based on the result of the first phase, one deep borehole (either DA-2-1 or DA-2-2) will be decid-ed to drill into R1 or R3 at a depth of about 2.5 km.

As indicated in Figure 4.14, the earthquake activity below 13-14 km depth is sparse in the vicinity of the reflection profile and the continuation of the fault is poorly constrained by the earthquakes. Another scenario is the

52

existence of a flower structure (Juhlin et al. 2010) in which the R1 and R3 faults converge to a master fault at depths of about 11-12 km. As seen in Figure 4.14, there are deeper seismic events east of the study area that are in a better agreement with a listric or generally east-dipping fault structure. In addition, the seismicity in the northern and southern sections of the fault does not support a fault wide flower structure, but rather an east-dipping fault zone (Lindblom et al., 2015).

Based on the enhanced seismic images, locations of some boreholes were suggested for future scientific drilling (Fig. 4.14). The drilling will be per-formed in two phases; first one borehole will be drilled into each of the faults R1 and R3 at a depth of about 1 km and based on the results of the first phase, a borehole will be drilled into one of the faults to intersect it at a depth of about 2 km (Fig. 4.14).

4.3.2. Conclusions Three reflections, corresponding to three steeply dipping faults and some sub-horizontal reflections, were imaged using data from the new seismic reflection profile. The reflections were migrated to their approximate true locations using an FD migration method with a velocity model extracted from the 3D local seismic network tomography and were traced down to a depth of about 7.5 km. Although weak, the main Pärvie fault is imaged dip-ping at approximately 65° to the east. A subsidiary, 60° west-dipping fault is significantly more reflective. We interpret the subsidiary west-dipping fault as consisting of a thicker fracture zone than the main Pärvie fault, explaining the difference in the reflectivity of the faults. Alternatively, the faults can have similar thickness, but the impedance contrast to the host rock may be greater for R1 than R3. Recent earthquakes cluster at a depth of about 11.5 km. The location and depth of the earthquake epicenters are consistent with the intersection of the main fault and the subsidiary fault at depth, assuming the faults are planar. Below 8 km depth the main fault may continue at a similar dip, become listric or form the lower part of a flower structure. Cur-rent microearthquakes support the two former alternatives.

53

4.4. Paper IV: Seismic signature of a complicated folding system in Western Australia

4.4.1. Summary The Cue-Weld Range area is located in the Murchison domain which forms the north-west part of the Youanmi terrane (Fig. 4.15), a tectonic unit within the Neoarchean Yilgarn Craton in Western Australia. The complex magmat-ic and deformation history of the Cue-Weld Range area (see Fig. 2.3) result-ed in a transposed array of greenstone belts that host significant iron, gold and base metal deposits. Crustal structures in the area and their relation to ore deposits are of major interest (Ivanic et al. 2013a; Romano et al. 2013; Zibra, et al. 2013). Banded iron formation (BIF) deposits are of particular interest to mine (Duuring et al. 2012; Duuring and Hagemann 2013b). In order to study the crustal structures and explore potential ore deposits at depth, the Geological survey of Western Australia (GSWA) in collaboration with Geoscience Australia acquired the deep seismic profile 10GA-YU1 with a length of about 300 km over the Youanmi terrane in 2010 (Ivanic et al. 2013a; Romano et al. 2013). An initial interpretation of the seismic data was carried out by the GSWA research team (Fig. 4.16) and they correlated pronounced reflections on the seismic image to regional structures (Romano et al. 2013). The previous interpretation of the data focused on providing a regional interpretation of subsurface structures down to Moho or deeper. In this study, we now investigate in detail a portion of the 10GA-YU1 profile where the geology is particularly complex and rocks outcrop. The studied segment on the 10GA-YU1 profile lies between the town of Cue and the Weld Range, with a length of about 22 km (Figs. 4.15 and 4.16). Outcrops along this portion of the line allow correlation of the rock units at the surface to their corresponding seismic signatures on the migrated seismic section.

The rocks in the area have been divided into three major groups (Van Kranendonk and Ivanic 2009): 1) mafic volcanic rocks, felsic volcaniclastic sandstones and banded iron-formations (BIF) of the Norie group, 2825-2800 Ma, 2) the Polelle group of mafic volcanic rocks, felsic volcanic and volcaniclastic sedimentary rocks and (BIF) with an age of 2785-2800 Ma, and 3) the Glen group of coarse clastic sedimentary rocks and komatiitic basalt, 2724-2700 Ma. Voluminous granitic rocks were emplaced extensive-ly from 2716 to 2592 Ma (Van Kranendonk and Ivanic, 2009) (Fig. 4.17). The two folding phases control the geological relationships in the area, spe-cifically: 1) the D1 deformation that produced the E-trending folds and oc-curred at the same time as granitic rock emplacement at about 2676 Ma and 2) the regional D2 shortening in the east-west direction which formed the N- to NNE-trending folds from 2660 to 2637 Ma.

54

BIF-hosted iron deposits in the Weld Range (Duuring and Hagemann 2013b) occur within a gabbroic complex (Fig. 3a). Duuring and Hagemann (2013a) suggested that the BIF deposits containing silica- and Fe oxide-rich bands formed in the BIF layers and later were transposed by the F1 and F2 folds and faults.

Figure 4.15. A general geologic map of the Youanmi region (after Wyche et al. 2013), the rectangle depicts the Cue-Weld Range area. The yellow stars show loca-tion of some valuable mineralizations in the region. The inset shows the location of the study area in Western Australia.

Romano et al. (2013) interpreted the north-western part of the 10GA-YU1 seismic profile (Fig. 4.16b) and focused on the deep reflections to conclude that Moho is located at the base of moderate reflective structures at a depth of about 33 km. Although there are some pronounced reflections near the surface, the lower crust is generally more reflective than the upper crust (Romano et al. 2013). Minor and major faults and shear zones were correlat-ed to the strong upper crustal reflections which may be traced down to depths greater than 15 km (Fig. 4.16b) (Romano et al., 2013). In the previous interpretation based on the available geological map (Fig. 4.17), shallow reflections were linked to mafic volcanic rocks of the Glen group, basaltic

55

rocks of the Polelle group and dolerites that had intruded into granitic rocks in the CDP range of 8000-9000 (Fig. 4.16b). Some reflective structures were not included in the interpretation provided by Romano et al. (2013) and marked as Form lines on the section (Fig. 4.16b).

Figure 4.16. A section of the previous seismic profile 10GA-YU1 across the north-western part of the Youanmi terrane: a) the processed section by Geoscience Aus-tralia. Location of the portion with 22 km length used in this study is shown by ar-rows; b) the interpreted section provided by GSWA (after Romano et al. 2013).

F

igur

e 4.

17. G

eolo

gica

l m

ap o

f th

e ar

ea (

mod

ifie

d af

ter

Ivan

ic a

t al

., 20

13b)

: a)

ove

rvie

w o

f th

e ar

ea. R

epet

ition

of

the

laye

rs i

n th

e so

uthe

rn s

truc

ture

s is

evi

denc

e of

the

fir

st f

old

syst

em w

hich

tre

nds

east

-wes

t (l

ocal

ly p

rese

rved

) in

the

area

. Not

e th

at th

ere

are

two

arro

w p

atte

rns

in th

e no

rthe

rn a

nd s

outh

ern

part

s of

the

area

. Bot

h th

ese

patte

rns

are

indi

ca-

tive

of

a N

to N

E-t

rend

ing

seco

nd f

old

syst

em. T

he r

ecta

ngle

sho

ws

the

loca

tion

of

the

inse

t map

and

Fig

ure

4.17

b. T

hein

set

show

s a

stac

k of

thr

ust

imbr

icat

es;

b ) t

he g

eolo

gy o

f th

e st

udy

area

is

show

n in

mor

e de

tail.

The

roc

ks g

ener

ally

dip

stee

ply.

The

fir

st (

F1)

and

seco

nd (

F2)

gen

erat

ion

of f

olds

are

sho

wn

on th

e m

ap.

57

The relationship between the interpreted shear zones and their surrounding lithologies in a 3D framework was not provided. We use a portion of the same line, but focus even more on linking the shallow reflections to the sur-face geology and taking into account the 3D structure that we can determine from the available geological map. In this study, we link these available field data and the most recent regional geological map of the area to the repro-cessed seismic line recently published by Ahmadi et al. (2014) (Fig. 4.18). The seismic profile (Fig. 4.18) shows key characteristic features of a super-imposed fold pattern which has resulted from the different folding phases. We suggest that these folded reflections represent a complex fold/thrust pat-tern. Two phases of superimposed folding are indicated on the geological map, as well as an imbricate stack of thrusts (Fig. 4.17). A conceptual model in Figure 4.19 shows how shortening may have occurred during the defor-mation where an imbricate stack of three layers/thrusts is folded in two di-rections. Based on the geologic map of the Cue-Weld Range (Fig. 4.17), the rock units of the Glen and Polelle groups appear to have been repeated in the NE and NW directions. In addition, some sub-horizontal reflections are ob-served along the seismic profile (marked by arrows in Fig. 4.18a). These reflections may represent segments of layers which appear sub-horizontal in the troughs of the superimposed folding phases.

Figure 4.20 shows a 3D geological model of three thrusts and three folded horizons which have been identified on the geological map and the seismic section. The thrusts with markers T1 and T2 are most likely associated with the folded reflections on the seismic section and therefore are constrained by those reflections (Fig. 4.18) and their traces at the surface (Fig. 4.17). These pronounced reflections can be linked to lithologies mapped at the surface using 3D constant velocity ray-tracing (see Fig. 3.5). Horizon A2 almost certainly corresponds to the foliated mafic volcanic rock of the Glen Group (Fig. 4.17), while horizon A1 and A3 likely correspond to the folded basaltic or metabasaltic rocks of the Polelle Group which are covered by granitic rocks and not exposed along the trace of the seismic profile (Figs. 4.20c and 4.20d).

58

Figure 4.18. The migrated and depth converted section: a) the migrated section with location of shotid 1908 marked (see Fig. 4.22), subhorizontal reflections are shown with arrows; b) the interpreted structures overlain on the migrated section. The no-menclatures correspond to the interpreted reflections in Figures 4.20 and 4.22.

Figure 4.19. A conceptual structural model of an imbricate stack of three lay-ers/thrusts which are folded in two directions: a) the thrusts are imbricated as a result of a SW-NE trending deformation (D1); b) the thrusts are refolded in the NW direc-tion as a result of the second deformation phase (D2).

To test the validity of the proposed model, seismic acoustic finite-difference forward modeling was carried out based on the interpreted profile (Fig. 4.18b). We have used a 2D FD code since cross-dip analysis of the data indicated that reflectivity coming from out-of-the-plane of the profile is mi-

59

nor (Ahmadi et al., 2014). Even though the working model is 3D (Fig. 4.20), the cross-dip component of the strongly reflecting horizons is quite gentle (+5°). The interpreted structures were used to build the model, which has a total size of 22× 6 km2 with a grid size of 27×27 m2. A Ricker wavelet with a dominant frequency of 60 Hz was used to perform the modeling. Figure 4.21 shows the density and velocity models used in the modeling process. Com-parison and analysis between the synthetic data generated using the 2D acoustic FD modeling method and the real data further supports the suggest-ed interpretation (Fig. 4.22). A number of hyperbolic reflections are obvious in the modeled gather, corresponding to the folded structures in the model (Fig. 4.22). Reflections A1 dip to the northwest and appear at 1 s down to 1.5 s and are consistent with similar reflections observed in the real data (Fig. 4.22b).

Figure 4.20. The 3D geological model of three folded thrusts (red) and reflections (green horizons): a) an bird’s eye view from the northwest, the folded reflections A1, A2 and A3 are picked and correlated to corresponding rock units mapped at the surface; b) same figure as Fig. 4.20a with transparent geological map overlaid; c) view from the top. Because of the complexity and lack of constraint laterally, we do not extend the layers too much; d) same figure as Fig. 4.20c with the transparent geological map overlaid. layer A2 is interpreted to be as the mafic volcanic rocks of Glen Group and layer A1 and A3 are interpreted to be as the basaltic and metabasaltic rocks of Polelle Group.

In the synthetic data two linear reflections can be observed (reflections C1 and C2) which project to the surface at about channels 25 and 50, corre-sponding to about CDP numbers 8530 and 8554, respectively (Figs. 4.22a and 4.17b). These linear reflections are from the near-vertical faults in the model at a distance of about 4 and 5 km, respectively, along the processed profile (Fig. 4.18). Reflection set B corresponds to the southeast-dipping structures in the model. Although not obvious, reflections corresponding to

60

group A2 in the synthetic data can be discerned in the real data (Fig. 4.22). The migrated sections show the structures in approximately their true spatial location in the subsurface (Fig. 4.23). The sections are depth-converted, al-lowing for a more realistic comparison between the modeled and the real data. The folded structures in the model occur at a depth of about 1 km down to 4 km and show good correspondence with the real data (Fig. 4.23).

Figure 4.21. The rock property models used in the forward modeling: a) the velocity model; b) the density model. The folded reflections are shown in red color and the faults are shown in blue in the models along with their corresponding physical prop-erties.

The synthetic reflections appear at the same depth as the reflections in the real data. The layered pattern and curvature of the reflections in the synthetic data follow the same patterns as in the real data (Fig. 4.23). The southeast-dipping structures are more pronounced in the synthetic data section at about CDP 8500 and a depth of about 1.5 km (Fig. 14.23a) than in the real data (Fig. 4.23b). However, indications of the reflecting structures are seen in the real section at the same location. The reflections from the near-vertical faults in the synthetic stacked section migrate to points near the surface and are not imaged at depth in the synthetic migrated section (Fig. 4.23a).

61

Figure 4.22. Shot gather comparison: a) the synthetic shot gather; b) the real shot gather. The recorded reflections group A1 are in a good agreement with their corre-sponding reflections in the real data (also marked as A1).

Figure 4.23. Comparison of migrated sections: a) the synthetic migrated section; b) the real migrated section. The sections are depth converted. The folded reflections in the synthetic data are consistent with the real data.

62

Comparison between the real and synthetic migrated sections suggests that the geometry of the seismic model may be representative of the subsurface in the Cue-Weld Range area. Duuring and Hagemann (2013a) classified the BIF deposits as having formed during the main deformation phases in the Cue-Weld Range district. These ore deposits include: 1) Quartz-magnetite meta-BIF, 2) Quartz-magnetite BIF and 3) Chert-Hematite BIF. The first group, Quartz-magnetite meta BIF, was interpreted to have formed during the first two deformation phases (Duuring and Hagemann 2013a), one short-ening in the N-S direction which resulted in an imbricate stack and the se-cond, shortening and folding in the E-W direction (F2 folds). This group is associated with greenschist facies metamorphism in the area and also hypogene hydrothermal processes (Duuring and Hagemann, 2013a). The other ore deposit groups are mainly associated with later brittle deformation phases and weathering (Duuring and Hagemann, 2013a). This evolution of the ore deposits is consistent with our interpretation and implies that the Quartz-magnetite meta BIF ore types might be found at different levels in the sub-surface since they are folded at least twice in the area. Therefore, the seismic reflection survey can potentially be used to identify fold structures that are associated with BIF ore bodies in the area, an indirect product of this study for mineral exploration.

4.4.2. Conclusions Different deformation phases were involved in forming the complex geology in the Cue-Weld Range area. We propose at least two folding deformation phases to be responsible for the observed structural pattern. The lithologies are folded in two near orthogonal directions (NE and NW directions) form-ing a class-2 fold interference pattern which is visible on the geologic map. Based on 3D constant velocity ray-tracing, some shallow reflections are linked to rocks outcropping at the surface. Forward seismic modeling was carried out to test the interpretation using a 2D acoustic finite difference code. The synthetic data generated correspond well with the acquired data, suggesting that the new 2D interpretation is consistent with the subsurface model proposed for the area. However, given that the outcrop is limited, as is the mapping, and that the seismic profile is only 2D other interpretations and models are possible. This study shows how the seismic reflection method can be used to better understand the structures at depth, both folds and de-formation zones, that can potentially be associated with ore bodies.

63

5. Conclusions

The seismic reflection method was used to image upper crustal structures in four papers. In each study, the capability of the method was studied, and corresponding problems associated with hardrock environments were ad-dressed.

In paper I, a 2D high resolution seismic profile was acquired over the Lappberget sulfide ore body in Garpenberg, central Sweden. The seismic data appeared noisy because of the mine-related activities at the time of the survey. In order to process the data, a conventional processing method was employed. Applying a cross-dip correction of -15° significantly enhanced a bright anomaly. Furthermore, the DMO correction was successful in han-dling conflicting dips in the seismic image. 3D ray-tracing was performed to provide an aid for the data interpretation. The bright reflections on the stacked section were interpreted to originate within and near the top of the ore body. The study confirmed the possibility of detecting sulfide deposits in the area using the seismic reflection method. However, interpretation of a single 2D seismic profile is challenging. In order to understand the 3D framework of the subsurface, a minimum of one 2D seismic profile perpen-dicular to the profile used in this work is necessary to provide some more constraints.

In paper II, the previous work presented in paper I was extended. The ef-fect of the backfilled stopes on seismic imaging in the Lappberget mine was studied. In order to simulate reflections/diffractions with a better approxima-tion, 3D acoustic FD modeling was performed. The modeling was carried out for both the ore body with the stopes included and without the stopes. Comparison between the synthetic data for both cases and the real data sug-gested it was likely that the response of the ore body was obscured by the response from the mined-out stopes. In addition, the true origin of the bright reflections/diffractions on the migrated section is not entirely obvious and further modeling with a detailed map of the mine infrastructures is needed to address this problem.

In paper III, the deep geometry of the Pärvie fault system in northern Sweden was revealed using a 2D seismic reflection profile acquired in 2014. The seismic data were processed using a conventional processing method. The NMO correction with a high constant velocity resulted in the best im-age. The main east-dipping fault and a subsidiary west-dipping fault were traced down to a depth of about 8 km. These two faults appeared with differ-

64

ent reflectivity characteristics. The main fault showed weak reflectivity while the west-dipping subsidiary fault presented strong reflectivity. The difference in the reflectivity of the faults may be because of the difference in the thickness of the faults and their corresponding fracture zones. However, the impedance contrast to the surrounding rocks may be greater in the west-dipping fault than the main Pärvie fault. An extrapolation of the two faults to depth shows that they intersect at about 11.5 km, coincident with a cluster of microearthquakes in the subsurface. The occurrence of the seismic events at these depths may be associated with the intersection of the faults. Based on the study, possible locations of scientific boreholes were suggested. Further investigations using several sensors to measure petrophysical and chemical properties of rocks and fluids along the faults will take place in the bore-holes. In order to image the faults and their host rocks with more resolution, the VSP method can be used.

In paper IV, a portion of a deep 2D seismic section was studied and de-tailed 3D and 2D interpretations were made for the subsurface of the Cue-Weld Range area in Western Australia. The seismic section displays a key characteristic of a superimposed folding system. Interpretation of the geolog-ical map and repetition of greenstone units in roughly two perpendicular directions indicates the occurrence of two generations of deformation in the area. Furthermore, 3D constant velocity ray-tracing allowed the correlation of shallow reflections visible in pre-stack data with mafic volcanic rocks of the Glen Group and basaltic rocks of the Polelle Group in the area. In order to test the model, a 2D FD method was performed and comparison between the synthetic and real data confirmed the proposed structural model. The study suggested using the seismic reflection method to help locate mineral deposits associated with the folding in the area is plausible. Acquiring more 2D data or 3D seismic data in the area may reveal more information about the complicated structures and deformation in the Cue-Weld Range area.

65

6. Summary in Swedish

Prospekteringsseismologi etablerades som geofysisk metod redan 1923 och har använts i stor utsträckning för prospektering av naturresurser (Sheriff and Geldart 1995; Dentith and Mudge 2014). Grundprincipen är att seismiska vågor färdas genom jordskorpan och kan refrakteras, reflekteras eller diffrakteras om de stöter på gränsskikt där den seismiska impedansen ändras. Reflekterade och refrakterade seismiska vågor kan ge viktig information om berggrunden och geologiska strukturer. Det finns flera olika metoder inom prospekteringsseismologi som kan användas beroende på syftet med undersökningen. En av de vanligaste metoderna är reflektionsseismik, vid vilken reflekterade eller diffrakterade seismiska vågor analyseras (Telford et al 1990; Sheriff and Geldart 1995).

Kunskap om geologiska strukturer och dess orientering är av stor vikt vid olika typer av projekt, till exempel för geotekniska undersökningar, mineralprospektering, gruvplanering eller för volymberäkning av olje- och gas-reservoarer (Malehmir et al. 2012a; Telford et al. 1990; Dentith and Mudge 2014). Även om reflektionsseismik ursprungligen användes huvudsakligen i områden med sedimentära bergarter för olje- och gas-prospektering så har tekniken visat sig användbar också för mineralprospektering (till exempel Malehmir et al. 2012a; Schmidt 1959; Wu et al. 1995; Eaton et al. 1996; Milkereit et al. 2000; Adam et al. 2003; Bohlen et al. 2003). En av styrkorna med metoden är bra signalpenetration vilket gör att djupt liggande kroppar och strukturer kan avbildas (Malehmir et al. 2012a). Reflektionsseismik kan också användas för att avbilda jordskorpan (till exempel Lie and Andersson 1998; Juhojuntti et al. 2001; Hobbs et al. 2006; Martínez Poyatos et al. 2012). Att förstå geometrin av jordskorpans strukturer är viktigt och kan användas i riskbedömning för jordbävningar, studier av geodynamiska processer och för prospektering efter mineraliseringar relaterade till geologiska strukturer (till exempel Williams et al. 1995; Hurich 1996; Larkin et al. 1996; Haugen and Schoenberg 2000; Juhlin et al. 2000; Juhojuntti et al. 2001; Hobbs et al. 2006; Goleby et al. 2006; Ji and Long 2006; Juhlin et al. 2010; Campbell et al. 2010; Juhlin and Lund 2011; Cheraghi et al. 2011; Lundberg et al. 2012; Ehsan et al. 2014; Jackson et al. 2015).

I denna avhandling har tre olika typer av problemställningar genomlysts med hjälp av reflektionsseismik. Det första problemet var att avbilda en sulfidmineralisering i en aktiv gruva i centrala Sverige. Fall nummer två

66

fokuserar på djupa förkastningsstrukturer relaterade till tiden vid slutet av glaciationen i norra Sverige. Det tredje problemet som avhandlas är nytolkning av en 2D reflektionsseismisk profil över komplex geologi i västra Australien.

Gemensamt för alla tre projekt är att de seismiska data har samlats in över kristallint berg. Seismisk avbildning i en kristallin miljö är utmanande och har avhandlats av många forskare jorden runt (till exempel: Malehmir et al. 2012a; Adam et al. 2003; Bohlen et al. 2003; Salisbury et al. 2003; Dehg-hannejad et al. 2010; Bongajum et al. 2012; Greenhalgh and Manukyan 2013).

Denna avhandling är en sammanfattning av fyra papers och består av fyra kapitel. Kapitel ett ger en introduktion och bakgrund till studien. Kapitel två presenterar geologisk bakgrund till de tre studieområdena. I kapitel tre sammanfattar jag alla metoder som jag har använt, inklusive numeriska och experimentella metoder. Slutligen, i kapitel fyra, sammanfattas alla papers som avhandlingen är baserad på.

I Paper I användes reflektionsseismik för att avbilda en sulfidmalmkropp och omgivande bergarter och strukturer i Garpenberg i södra Dalarna. Profiler med högupplöst 2D reflektionsseismik mättes, vilket var en utmaning eftersom mätningarna utfördes i en bullrig miljö med en aktiv gruva. De seismiska data processerades längs en rak CDP linje och tolkningen baserades på en migrerad profil. En 3D-modellering av de seismiska vågornas gångväg utfördes för att testa tolkningen. Ett antal stuffer tillhandahölls från gruvbolaget Boliden, dessa användes för att med ett ultraljudinstrument mäta typiska P-våghastigheter på sulfidmalmen. Studien visade att det är möjligt att detektera malmkroppen, men att det är en utmaning att göra det med 2D-seismiska profiler.

Paper II är en uppföljning till Paper I, här inkluderades de utbrutna och igenfyllda hålrummen i gruvan. Hålrummen i gruvan fylls igen med en pasta bestående av gråberg och cement. Ytterligare stuffer samlades in från gruvan för mätning av P-vågshastigheten. För att erhålla en bättre approximation till vågekvationen användes den finita differensmetoden. Modelleringen visar att de återfyllda hålrummen i gruvan gör det svårt att seismiskt avbilda malmkroppen Lappberget.

I Paper III studerades de djupa strukturerna i Pärvie-förkastningen. De avbildades ner till ungefär 8 km med en 2D-reflektionsseismisk profil. En samtolkning med data från små jordbävningar gjordes. Lämpliga platser för framtida vetenskaplig borrning igenom förkastningssystemet föreslogs baserat på tolkningen.

I Paper IV studerades den seismiska signaturen av den komplexa geologin i bergskedjan Cue-Weld Range i västra Australien. De tydliga strukturerna i den tvådimensionella seismisk profilen stämmer överens med motsvarande berggrundsenheter på den geologiska kartan. En tredimensionell modellering av de seismiska vågornas gångväg (under antagandet konstant hastighet)

67

utfördes för att avgränsa tolkningsmöjligheterna. Tolkningen testades också i 2D med den finita differensmetoden.

68

Acknowledgments

I am very grateful to God for the good health and prosperity that was neces-sary to complete my PhD project.

I would like to express deep gratitude to my supervisor Christopher Juhlin for his patient guidance, kind support and encouragement. I also appreciate him for allowing me to participate in workshops and conferences. His care-ful supervision and useful critiques have illuminated my scientific journey.

I am grateful to Alireza Malehmir for all the advice he has given me. I have always found it constructive to discuss scientific problems with him and his comments and suggestions have immensely promoted my research work.

I would like to offer my special thanks to Hemin Koyi, who taught me how to look at geological maps carefully. I was very lucky to have the chance to learn more about structural geology from the discussions we had.

I wish to acknowledge my co-authors for their collaboration on all of the papers which this thesis is based on: Alireza Malehmir, Björn Lund, Chris-topher Juhlin, Hemin Koyi, Klaus Gessner (Geological Survey of Western Australia), Maria Ask (Luleå University of Technology) and Mie Munk (Bo-liden). Likewise, the anonymous researchers from the Geological Survey of Western Australia should be acknowledged for their constructive comments on paper IV. I wish to thank Ernst Schetselaar (Geological Survey of Cana-da) for his advice on geological modeling and Lars Dynesius for his help with the ultrasonic measurements in this thesis. I am particularly grateful to Hans Palm for his support and advice during the seismic reflection surveys with which I was involved. I would also like to acknowledge Boliden Mines for providing information and rock samples regarding the Garpenberg pro-ject and papers. In addition, Geoscience Australia is acknowledged for providing the seismic data across the Younmi terrane in Western Australia.

It is a pleasure to thank Alireza Malehmir, Bojan Brodic, Christopher Juhlin, Monika Ivandic and Suman Mehta for providing critical suggestions which improved the dissertation. My great thanks go to Magnus Andersson for his kind help with the Swedish summary. I would like to acknowledge Karin Högdahl for her advice on some geological terms in this thesis.

69

The Society of Exploration Geophysicists (SEG) and the Solid Earth (SE) journal are acknowledged for their permission to reprint published papers. I would further like to acknowledge the University of Kurdistan and the Min-istry of Science, Research and Technology of Iran for their support and fund-ing, which allowed me to study abroad. I would also like to thank Christo-pher Juhlin, Roland Robert and the board members of the department of Earth Sciences at Uppsala University for allowing me to extend my study.

I am indebted to Mahdieh (Azita) Dehghannejad, Saeid Cheraghi and Taher Mazloomian who advised and helped me during the first days of settling in Uppsala. I wish to thank the following colleagues and friends who have in-spired and helped me throughout my PhD course: Ali (Ashkan) Dorostkar, Amir Abdi, Ari Tryggvason, Bjarne Almqvist, Bojan Brodic, Daniel Sopher, Dibakar Ghosal, Emil Lundberg, Ershad Gholamrezaie, Faramarz Nilfouroushan, Fatemeh (Sahar) Sharifi Brojerdi, Fei Huang, Fengjiao Zhang, Georgiana Maries, Hamzeh Sadeghisorkhani, Harbe Anwar Muhamad, Ida Lydersen, Joachim Place, Jochen Kamm, Magnus Andersson, Mahdieh (Azita) Dehghannejad, María de los Angeles García Juanatey, Monika Ivandic, Pouya Ahmadi, Remi Vachon, Peter Hedin, Peter Schmidt, Ping Yan, Saeid Cheraghi, Samar Amini, Sayyed Mohammad Abtahi, Silvia Salas Romero, Shunguo Wang, Suman Mehta, Théo Berthet and Zeinab Jeddi.

I owe a lot to my family, who has supported me generously. I believe this success would not have been possible without their help and encouragement.

Last, but not least, my deepest gratitude goes to my wife, Maryeh, who has stood by me throughout my scrambles and absences. She has given me in-valuable support and help, inspiring me in my career.

70

References

Adam E, Perron G, Grant A, Matthews L, Milkereit B (2003) 3D Seismic Imaging for VMS Deposit Exploration, Matagami, Quebec. In: Hardrock Seismic Explo-ration. Society of Exploration Geophysicists, pp 229–246

Ahmadi O, Juhlin C, Gessner K (2014) New Insights from Seismic Imaging Over the Youanmi Terrane, Yilgarn Craton,Western Australia. Energy Procedia 59:113–119. doi: 10.1016/j.egypro.2014.10.356

Aki K, Richards PG (2009) Quantitative seismology. University Science Books, Sausalito, CA

Alkhalifah T (2000) An acoustic wave equation for anisotropic media. Geophysics 65:1239–1250. doi: 10.1190/1.1444815

Allen R, Ripa M, Jansson N (2008) Palaeoproterozoic volcanic- and limestone-hosted Zn-Pb-Ag-(Cu-Au) massive sulphide deposits and Fe oxide deposits in Bergslagen, Sweden. IGCP Project 502

Ayarza P, Juhlin C, Brown D, Beckholmen M, Kimbell G, Pechnig R, Pevzner L, Pevzner R, Ayala C, Bliznetsov M, Glushkov A, Rybalka A (2000) Integrated geological and geophysical studies in the SG4 borehole area, Tagil Volcanic Arc, Middle Urals: Location of seismic reflectors and source of the reflectivity. J Geophys Res Solid Earth 105:21333–21352. doi: 10.1029/2000JB900137

Bancroft JC (2007) A Practical Understanding of Pre- and Poststack Migrations: Prestack. SEG Books

Bean CJ, Martini F (2010) Sub-basalt seismic imaging using optical-to-acoustic model building and wave equation datuming processing. Mar Pet Geol 27:555–562. doi: 10.1016/j.marpetgeo.2009.09.007

Bellefleur G, Malehmir A, Müller C (2012) Elastic finite-difference modeling of volcanic-hosted massive sulfide deposits: A case study from Half Mile Lake, New Brunswick, Canada. Geophysics 77:WC25–WC36. doi: 10.1190/geo2011-0445.1

Bellefleur G, Schetselaar E, White D, Miah K, Dueck P (2015) 3D seismic imaging of the Lalor volcanogenic massive sulphide deposit, Manitoba, Canada. Geophys Prospect 63:813–832. doi: 10.1111/1365-2478.12236

Blake B, Figueroa D, Manceda R, Oller J, Hofland G (1999) 3D forward ray trace seismic modeling of strike lines in complex geology. In: SEG Technical Pro-gram Expanded Abstracts 1999. Society of Exploration Geophysicists, pp 1871–1874

Blanch JO (1996) A study of viscous effects in seismic modeling, imaging, and inversion: Methodology, computational aspects, and sensitivity. PhD thesis, Rice University

Bödvarsson R, Lund B (2003) The SIL Seismological data Acquisition System -as Operated in Iceland and in Sweden. In: Methods and Applications of Signal Processing in Seismic Network Operations. Springer, Berlin, pp 131–148

Bohlen T (2002) Parallel 3-D viscoelastic finite difference seismic modelling. Comput Geosci 28:887–899. doi: 10.1016/S0098-3004(02)00006-7

71

Bohlen T, Müller CF, Milkereit B (2003) Elastic wave scattering from massive sul-fide orebodies: on the role of composition and shape. In: Hardrock Seismic Ex-ploration, Eaton, D., Milkereit, B. and Salisbury, M. Edited by Eaton, D., Milkereit, B. and Salisbury, M., SEG Developments in Geophysics series, 10, pp 70–89

Bohlen T, Nil DD, Köhen D, Jetschny S (2015) SOFI3D - Seismic Modeling with Finite differences, 3D-acoustic and viscoelastic version- User Guide.

Bongajum E, Milkereit B, Adam E, Meng Y (2012) Seismic imaging in hardrock environments: The role of heterogeneity? Tectonophysics 572–573:7–15. doi: 10.1016/j.tecto.2012.03.003

Broggini F, Snieder R (2012) Connection of scattering principles: a visual and mathematical tour. Eur J Phys 33:593. doi: 10.1088/0143-0807/33/3/593

Buske S, Bellefleur G, Malehmir A (2015) Introduction to special issue on “hard rock seismic imaging.” Geophys Prospect 63:751–753. doi: 10.1111/1365-2478.12257

Campbell FM, Ghisetti F, Kaiser AE, Green AG, Horstmeyer H, Gorman AR (2010) Structure and evolution of the seismically active Ostler Fault Zone (New Zea-land) based on interpretations of multiple high resolution seismic reflection pro-files. Tectonophysics 495:195–212. doi: 10.1016/j.tecto.2010.09.016

Cerveny V (1989) Ray tracing in factorized anisotropic inhomogeneous media. Geophys J Int 99:91–100.

Cerveny V (2005) Seismic Ray Theory. Cambridge University Press Cheraghi S, Craven JA, Bellefleur G (2015) Feasibility of virtual source reflection

seismology using interferometry for mineral exploration: A test study in the Lalor Lake volcanogenic massive sulphide mining area, Manitoba, Canada. Geophys Prospect 63:833–848. doi: 10.1111/1365-2478.12244

Cheraghi S, Malehmir A, Bellefleur G (2011) Crustal-scale reflection seismic inves-tigations in the Bathurst Mining Camp, New Brunswick, Canada. Tectonophysics 506:55–72. doi: 10.1016/j.tecto.2011.04.011

Dehghannejad M, Juhlin C, Malehmir A, Skytta P, Weihed P (2010) Reflection seismic imaging of the upper crust in the Kristineberg mining area, northern Sweden. J Appl Geophys 71:125–136.

Dentith PM, Mudge ST (2014) Geophysics for the Mineral Exploration Geoscientist, 1 edition. Cambridge University Press, Cambridge, United Kingdom

Deregowski S (1991) What Is DMO? In: Dip Moveout Processing. Society of Ex-ploration Geophysicists, p A–1–A–18

Duuring P, Hagemann S (2013a) Genesis of superimposed hypogene and supergene Fe orebodies in BIF at the Madoonga deposit, Yilgarn Craton, Western Austral-ia. Miner Deposita 48:371–395. doi: 10.1007/s00126-012-0429-0

Duuring P, Hagemann S (2013b) Leaching of silica bands and concentration of magnetite in Archean BIF by hypogene fluids: Beebyn Fe ore deposit, Yilgarn Craton, Western Australia. Miner Deposita 48:341–370. doi: 10.1007/s00126-012-0428-1

Duuring P, Hagemann SG, Novikova Y, Cudahy T, Laukamp C (2012) Targeting Iron Ore in Banded Iron Formations Using ASTER Data: Weld Range Green-stone Belt, Yilgarn Craton, Western Australia. Econ Geol 107:585–597. doi: 10.2113/econgeo.107.4.585

Eaton D, Guest S, Milkereit B, Bleeker W, Crick D, Schmitt D, Salisbury MH (1996) Seismic imaging of massive sulfide deposits; Part III, Borehole seismic imaging of near-vertical structures. Econ Geol 91:835–840. doi: 10.2113/gsecongeo.91.5.835

72

Edfelt Å, Sandrin A, Evins P, Jeffries T, Storey C, Elming S-Å, Martinsson O (2006) Stratigraphy and tectonic setting of the host rocks to the Tjårrojåkka Fe-oxide Cu-Au deposits, Kiruna area, northern Sweden. GFF 128:221–232. doi: 10.1080/11035890601283221

Ehsan SA, Carbonell R, Ayarza P, Martí D, Pérez-Estaún A, Martínez-Poyatos DJ, Simancas JF, Azor A, Mansilla L (2014) Crustal deformation styles along the reprocessed deep seismic reflection transect of the Central Iberian Zone (Iberian Peninsula). Tectonophysics 621:159–174. doi: 10.1016/j.tecto.2014.02.014

Eklind M, Ericsson P, Jönsson J, Lewen M, Nord G, Samuelsson B, Potts A, Casteel K, Ericsson M (2007) Mining Methods in Underground Mining. edited by Smith E., Atlas Copco

Farra V (1993) Ray tracing in complex media. J Appl Geophys 30:55–73. doi: 10.1016/0926-9851(93)90018-T

Fichtner A (2011) Full Seismic Waveform Modelling and Inversion. Springer Berlin Heidelberg, Berlin, Heidelberg

Fleury C (2012) Wave propagation in complex media, scattering theory, and appli-cation to seismic imaging. PhD thesis, Colorado School of Mines

Fomel S, Sava P, Vlad I, Liu Y, Bashkardin V (2013) Madagascar: open-source software project for multidimensional data analysis and reproducible computa-tional experiments. J Open Res Softw. doi: 10.5334/jors.ag

Gjøystdal H, Iversen E, Laurain R, Lecomte I, Vinje V, Åstebøl K (2002) Review of ray theory applications in modelling and imaging of seismic data. Stud Geophys Geod 46:113–164.

Goleby BR, Blewett RS, Fomin T, Fishwick S, Reading AM, Henson PA, Kennett BLN, Champion DC, Jones L, Drummond BJ, Nicoll M (2006) An integrated multi-scale 3D seismic model of the Archaean Yilgarn Craton, Australia. Tectonophysics 420:75–90. doi: 10.1016/j.tecto.2006.01.028

Greenhalgh S, Manukyan E (2013) Seismic Reflection for Hardrock Mineral Explo-ration: Lessons from Numerical Modeling. J Environ Eng Geophys 18:281–296. doi: 10.2113/JEEG18.4.281

Hale D (1984) Dip‐moveout by Fourier transform. Geophysics 49:741–757. doi: 10.1190/1.1441702

Halliday D, Resnick R, Walker J (2000) Fundamentals of Physics, 6th Edition, 6th edition. Wiley, New York

Hall RB (1988) Hardrock versus softrock geology. In: General Geology. Springer US, pp 331–331

Hamrin H (1980) Guide to underground mining methods and applications. Atlas Copco

Harman P, Leung L, Downey M (1987) Cross-hole seismic survey for mineral ex-ploration in the west Kimberley area, Western Australia. Explor Geophys 18:80–83.

Harrison CB, Urosevic M (2012) Seismic processing, inversion, and AVO for gold exploration — Case study from Western Australia. Geophysics 77:WC235–WC243. doi: 10.1190/geo2011-0506.1

Haugen GU, Schoenberg MA (2000) The echo of a fault or fracture. Geophysics 65:176–189. doi: 10.1190/1.1444708

Hobbs RW, Drummond BJ, Goleby BR (2006) The effects of three-dimensional structure on two-dimensional images of crustal seismic sections and on the in-terpretation of shear zone morphology. Geophys J Int 164:490–500. doi: 10.1111/j.1365-246X.2006.02814.x

73

Hughes SW (2005) Archimedes revisited: a faster, better, cheaper method of accu-rately measuring the volume of small objects. Phys Educ 40:468. doi: 10.1088/0031-9120/40/5/008

Hurich CA (1996) Kinematic evolution of the lower plate during intracontinental subduction: An example from the Scandinavian Caledonides. Tectonics 15:1248–1263. doi: 10.1029/96TC00828

Ivanic TJ, Korsch RJ, Wyche S, Jones LEA, Zibra I, Blewett RS, Jones T, Milligan P, Costelloe RD, Van Kranendonk MJ, Doublier MP, Hall CE, Romano SS, Pawley MJ, Gessner K, Patison N, Kennett BLN, Chen SF (2013a) Preliminary interpretation of the 2010 Youanmi deep seismic reflection lines and magnetotelluric data for the Windimurra Igneous Complex. In: Youanmi and Southern Carnarvon seismic and magnetotelluric (MT) workshop. compiled by S Wyche, TJ Ivanic and I Zibra: Geological Survey of Western Australia, Rec-ord 2013/6, pp 93–108

Ivanic TJ, Zibra I, Doublier MP, Wyche S (2013b) Geological Interpretation of the northwest Yilgarn Craton. In: Youanmi and Southern Carnarvon seismic and magnetotelluric (MT) workshop, Record 2013/6, Plate 1A. compiled by S Wyche, TJ Ivanic and I Zibra: Geological Survey of Western Australia, Record 2013/6, Plate 1A.

Jackson CA-L, Jackson MPA, Hudec MR, Rodriguez CR (2015) Enigmatic struc-tures within salt walls of the Santos Basin—Part 1: Geometry and kinematics from 3D seismic reflection and well data. J Struct Geol 75:135–162. doi: 10.1016/j.jsg.2015.01.010

Jansson NF (2011) The Origin of Iron Ores in Bergslagen, Sweden,and their Rela-tionships with Polymetallic Sulphide Ores. PhD thesis, Luleå University of Technology

Ji S, Long C (2006) Seismic reflection response of folded structures and implica-tions for the interpretation of deep seismic reflection profiles. J Struct Geol 28:1380–1387. doi: 10.1016/j.jsg.2006.05.003

Jones IF (2010) An Introduction to: Velocity Model Building. EAGE Publications, Houten

Juhlin C (1995) Finite-difference elastic wave propagation in 2D heterogeneous transversely isotropic media. Geophys Prospect 43:843–858. doi: 10.1111/j.1365-2478.1995.tb00284.x

Juhlin C, Dehghannejad M, Lund B, Malehmir A, Pratt G (2010) Reflection seismic imaging of the end-glacial Pärvie Fault system, northern Sweden. J Appl Geophys 70:307–316.

Juhlin C, Lund B (2011) Reflection seismic studies over the end-glacial Burträsk fault, Skellefteå, Sweden. Solid Earth 2:9–16. doi: 10.5194/se-2-9-2011

Juhlin C, Wahlgren C-H, Stephens MB (2000) Seismic imaging in the frontal part of the Sveconorwegian orogen, south-western Sweden. Precambrian Res 102:135–154. doi: 10.1016/S0301-9268(00)00064-4

Juhojuntti N, Juhlin C, Dyrelius D (2001) Crustal reflectivity underneath the Central Scandinavian Caledonides. Tectonophysics 334:191–210.

Julian BR, Gubbins D (1977) Three Dimensional Seismic Ray Tracing. J Geophys 43:95–113.

Kashubin AS, Juhlin C (2010) Mapping of crustal scale tectonic boundaries in the Ossa-Morena Zone using reprocessed IBERSEIS reflection seismic data. Tectonophysics 489:139–158. doi: 10.1016/j.tecto.2010.04.010

Koivisto E, Malehmir A, Hellqvist N, Voipio T, Wijns C (2015) Building a 3D model of lithological contacts and near-mine structures in the Kevitsa mining

74

and exploration site, Northern Finland: constraints from 2D and 3D reflection seismic data. Geophys Prospect 63:754–773. doi: 10.1111/1365-2478.12252

Larkin SP, Levander A, Okaya D, Goff JA (1996) A deterministic and stochastic velocity model for the Salton Trough/Basin and Range transition zone and con-straints on magmatism during rifting. J Geophys Res Solid Earth 101:27883–27897. doi: 10.1029/96JB02535

Larner KL, Gibson B, Chambers R, Wiggins RA (1979) Simultaneous estimation of residual static and crossdip corrections. Geophysics 44:1175–1192.

L’Heureux E, Milkereit B, Vasudevan K (2009) Heterogeneity and seismic scatter-ing in exploration environments. Tectonophysics 472:264–272. doi: 10.1016/j.tecto.2008.04.001

Lie JE, Andersson M (1998) The deep-seismic image of the crustal structure of the Tornquist Zone beneath the Skagerrak Sea, northwestern Europe. Tectonophysics 287:139–155. doi: 10.1016/S0040-1951(98)80065-6

Lindblom E, Lund B, Tryggvason A, Uski M, Bödvarsson R, Juhlin C, Roberts R (2015) Microearthquakes illuminate the deep structure of the endglacial Pärvie fault, northern Sweden. Geophys J Int 201:1704–1716. doi: 10.1093/gji/ggv112

Li T, Eaton D (2005) Delineating the Tuwu porphyry copper deposit at Xinjiang, China, with seismic-reflection profiling. Geophysics 70:B53–B60. doi: 10.1190/1.2122409

Lundberg E, Juhlin C (2011) High resolution reflection seismic imaging of the Ullared Deformation Zone, southern Sweden. Precambrian Res 190:25–34. doi: 10.1016/j.precamres.2011.07.012

Lundberg E, Juhlin C, Nasuti A (2012) High resolution reflection seismic profiling over the Tjellefonna fault in the Møre-Trøndelag Fault Complex, Norway. Solid Earth 3:175–188. doi: 10.5194/se-3-175-2012

Malehmir A, Andersson M, Lebedev M, Urosevic M, Mikhaltsevitch V (2013) Ex-perimental estimation of velocities and anisotropy of a series of Swedish crystal-line rocks and ores. Geophys Prospect 61:153–167. doi: 10.1111/j.1365-2478.2012.01063.x

Malehmir A, Durrheim R, Bellefleur M, Urosevic M, Juhlin C, White D, Milkerit B, Campbell G (2012a) Seismic methods in mineral exploration and mine plan-ning: A general overview of past and present case histories and a look into the future. Geophysics 77(5):WC173–WC190.

Malinowski M, Schetselaar E, White DJ (2012) 3D seismic imaging of volcanogenic massive sulfide deposits in the Flin Flon mining camp, Canada: Part 2 — For-ward modeling. Geophysics 77:WC81–WC93. doi: 10.1190/geo2011-0474.1

Manzi M, Cooper G, Malehmir A, Durrheim R, Nkosi Z (2015) Integrated interpre-tation of 3D seismic data to enhance the detection of the gold-bearing reef: Mponeng Gold mine, Witwatersrand Basin (South Africa). Geophys Prospect 63:881–902. doi: 10.1111/1365-2478.12273

Manzi MSD, Gibson MAS, Hein KAA, King N, Durrheim RJ (2012) Application of 3D seismic techniques to evaluate ore resources in the West Wits Line goldfield and portion of the West Rand goldfield, South Africa. Geophysics 77:WC163–WC171.

Margrave GF (2000) New seismic modelling facilities in Matlab. CREWES Re-search Report 12, University of Calgary

Margrave GF (2003) Numerical Methods of Exploration Seismology with algo-rithms in MATLAB. CREWES, Department of Geology and Geophysics, Uni-versity of Calgary

Martínez Poyatos D, Carbonell R, Palomeras I, Simancas JF, Ayarza P, Martí D, Azor A, Jabaloy A, González Cuadra P, Tejero R, Martín Parra LM, Matas J,

75

González Lodeiro F, Pérez-Estaún A, García Lobón JL, Mansilla L (2012) Im-aging the crustal structure of the Central Iberian Zone (Variscan Belt): The ALCUDIA deep seismic reflection transect. Tectonics 31. doi: 10.1029/2011TC002995

Milkereit B, Berrer EK, King AR, Watts AH, Roberts B, Adam E, Eaton DW, Wu J, Salisbury MH (2000) Development of 3-D seismic exploration technology for deep nickel-copper deposits—A case history from the Sudbury basin, Canada. Geophysics 65:1890–1899.

Nedimovic MR, West GF (2003) Crooked-line 2D seismic reflection imaging in crystalline terrains: Part 1, data processing. Geophysics 68:274–285.

Nelson RG (1984) Seismic reflection and mineral prospecting. Explor Geophys 15:229–250.

Read JJ (1989) Seismic reflection investigations of the Hellyer orebody and Que-Hellyer volcanics, North West Tasmania. Explor Geophys 20:159–162.

Robb LJ (2005) Introduction to ore-forming processes. Blackwell Pub, Malden, MA Rodriguez-Tablante J, Tryggvason A, Malehmir A, Juhlin C, Palm H (2007) Cross-

profile acquisition and cross-dip analysis for extracting 3D information from 2D surveys, a case study from the western Skellefte District, northern Sweden. J Appl Geophys 63:1–12. doi: 10.1016/j.jappgeo.2007.03.001

Romano SS, Ivanic TJ, Korsch RJ, Wyche S, Van Kranendonk MJ, Jones LEA, Zibra I, Blewett RS, Jones T, Milligan P, Costelloe RD, Doublier MP, Pawley MJ, Gessner K, Hall CE, Patison N, Kennett BLN, Chen SF (2013) Preliminary interpretation of the northern section of deep seismic line 10GA-YU1: Narryer Terrane to Murchison Domain of the Youanmi terrane. In: Youanmi and South-ern Carnarvon seismic and magnetotelluric (MT) workshop. compiled by S Wyche, TJ Ivanic and I Zibra: Geological Survey of Western Australia, Record 2013/6, pp 119–124

Saenger EH, Bohlen T (2004) Finite-difference modeling of viscoelastic and aniso-tropic wave propagation using the rotated staggered grid. Geophysics 69:583–591. doi: 10.1190/1.1707078

Saleh R, Milkereit B (2014) A Modelling Study Of Broadband Seismic Wave Prop-agation In A Deep Mine. In: 76th EAGE Conference and Exhibition 2014. Am-sterdam, Netherland,

Salisbury M, Harvey CW, Matthews L (2003) The acoustic properties of ores and host rocks in hardrock terranes. In: Hardrock Seismic Exploration, Eaton, D., Milkereit, B. and Salisbury, M. Edited by Eaton, D., Milkereit, B. and Salisbury, M., SEG Developments in Geophysics series, 10, pp 9–19

Sava P, Poliannikov O (2008) Interferometric imaging condition for wave-equation migration. Geophysics 73:S47–S61. doi: 10.1190/1.2838043

Sava P, Vlad I (2011) Wide-azimuth angle gathers for wave-equation migration. Geophysics 76:S131–S141. doi: 10.1190/1.3560519

Scales JA (1997) Theory of Seismic Imaging. Smizdat Press Schmidt G (1959) Results of Underground-Seismic Reflection Investigations in the

Siderite District of the Siegerland. Geophys Prospect 7:287–290. doi: 10.1111/j.1365-2478.1959.tb01470.x

Sheriff RE, Geldart LP (1995) Exploration seismology. Cambridge University Press, Cambridge

Squelch A, Madadi M, Urosevic M (2013) A reproducible framework for 3D acous-tic forward modelling of hard rock geological models with Madagascar. In: ASEG Extended Abstracts 2013: 23rd Geophysical Conference. Australian So-ciety of Exploration Geophysicists (ASEG), pp 1–4

76

Stanchits S, Vinciguerra S, Dresen G (2006) Ultrasonic Velocities, Acoustic Emis-sion Characteristics and Crack Damage of Basalt and Granite. In: Dresen G, Zang A, Stephansson O (eds) Rock Damage and Fluid Transport, Part I. Birkhäuser Basel, pp 975–994

Stephens MB, Ripa M, Lundström I, Persson L, Bergman T, Ahl M, Wahlgren C-H, Persson P-H, Wickström L (2009) Synthesis of the bedrock geology in the Bergslagen region, Fennoscandian Shield, south-central Sweden. Geological survey of Sweden

Telford WM, Telford WM, Sheriff RE (1990) Applied Geophysics. Cambridge University Press

Van Kranendonk MJ, Ivanic TJ (2009) A new lithostratigraphic scheme for the northeastern Murchison Domain, Yilgarn Craton. Geological Survey of Western Australia, Geological Survey of Western Australia, Perth

Vanorio T, Prasad M, Patella D, Nur A (2002) Ultrasonic velocity measurements in volcanic rocks: correlation with microtexture. Geophys J Int 149:22–36. doi: 10.1046/j.0956-540x.2001.01580.x

Wang X, Liu X (2007) 3-D acoustic wave equation forward modeling with topogra-phy. Appl Geophys 4:8–15. doi: 10.1007/s11770-007-0002-2

Weiss R, Shragge J (2013) Solving 3D anisotropic elastic wave equations on parallel GPU devices. GEOPHYSICS 78:F7–F15. doi: 10.1190/geo2012-0063.1

White DJ, Secord D, Malinowski M (2012) 3D seismic imaging of volcanogenic massive sulfide deposits in the Flin Flon mining camp, Canada: Part 1 — Seis-mic results. Geophysics 77:WC47–WC58. doi: 10.1190/geo2011-0487.1

Williams PF, Goodwin LB, Lafrance B (1995) Brittle faulting in the Canadian Ap-palachians and the interpretation of reflection seismic data. J Struct Geol 17:215–232. doi: 10.1016/0191-8141(94)E0046-2

Wu J (1996) Potential pitfalls of crooked‐line seismic reflection surveys. Geophysics 61:277–281. doi: 10.1190/1.1443949

Wu J, Milkereit B, Boerner DE (1995) Seismic imaging of the enigmatic Sudbury Structure. J Geophys Res 100:4117–4130.

Wyche S, Ivanic T, Zibra I, Gessner K, Doublier M, Korsch RJ, Blewett RS (2013) Youanmi and Southern Carnarvon traverses: implications for mineral systems. In: Youanmi and Southern Carnarvon seismic and magnetotelluric (MT) work-shop. compiled by S Wyche, TJ Ivanic and I Zibra: Geological Survey of West-ern Australia, Record 2013/6, pp 159–171

Yerneni S, Bhardwaj D, Chakraborty S, Phadke S (2002) Finite difference forward modeling for complex geological models. In: SEG Technical Program Expand-ed Abstracts 2002. Society of Exploration Geophysicists, pp 1987–1990

Yilmaz Ö (2001) Seismic Data Analysis: Processing, Inversion, and Interpretation of Seismic Data. Society of Exploration Geophysicists

Zibra, I, Gessner, K, Korsch RJ, Blewett RS, Jones T, Ayarza P, Jones LEA, Costelloe RD, Wyche, S, Doublier MP, Hall CE, Chen SF, Romano SS, Ivanic TJ, Pawley MJ, Patison N, Kennett BLN, Van Kranendonk MJ (2013) Main crustal-scale features of the Youanmi seismic transect. In: Youanmi and South-ern Carnarvon seismic and magnetotelluric (MT) workshop. compiled by S Wyche, TJ Ivanic and I Zibra: Geological Survey of Western Australia, Record 2013/6, pp 75–76

Acta Universitatis UpsaliensisDigital Comprehensive Summaries of Uppsala Dissertationsfrom the Faculty of Science and Technology 1269

Editor: The Dean of the Faculty of Science and Technology

A doctoral dissertation from the Faculty of Science andTechnology, Uppsala University, is usually a summary of anumber of papers. A few copies of the complete dissertationare kept at major Swedish research libraries, while thesummary alone is distributed internationally throughthe series Digital Comprehensive Summaries of UppsalaDissertations from the Faculty of Science and Technology.(Prior to January, 2005, the series was published under thetitle “Comprehensive Summaries of Uppsala Dissertationsfrom the Faculty of Science and Technology”.)

Distribution: publications.uu.seurn:nbn:se:uu:diva-259396

ACTAUNIVERSITATIS

UPSALIENSISUPPSALA

2015