Magnetotelluric (MT) Method

Size: px
Start display at page:

Download "Magnetotelluric (MT) Method"

Transcription

1 Magnetotelluric (MT) Method Dr. Hendra Grandis Graduate Program in Applied Geophysics Faculty of Mining and Petroleum Engineering ITB Geophysical Methods Techniques applying physical laws (or theory) to the study of the solid earth Estimation of subsurface physical property distribution by measuring relevant parameters : density gravity susceptibility magnetic field resistivity electric field, electromagnetic field seismic wave travel time velocity 1

2 Geo-electromagnetic electromagnetic Methods Dependence of electric and magnetic phenomena on the conductivity of the medium can be exploited to study the solid Earth Methods to estimate subsurface electrical property (resistivity)) distribution by measuring EM fields: Magnetotellurics (MT), Controlled-Source Audio-freq. MT (CSAMT), Transient EM, Very Low Freq. EM (VLF-EM), Ground Probing Radar (GPR) Electrical resistivities of rocks 2

3 What is magnetotellurics (MT) )? MT is a geophysical method to estimate subsurface electrical property (resistivity( or conductivity) distribution by measuring natural EM fields Source of MT signals comes from interaction of Earth s permanent magnetic field with particles from the solar wind and with atmospheric lightning no need for transmitter, simplifies the logistics random signals, low S/N (dead band ~ 1 Hz) Electromagnetic Induction transmitter generates time varying EM fields induce Eddy currents in the conductor (Earth) generate secondary magnetic field electric and magnetic fields sensed at the receiver 3

4 natural electromagnetic field f > 1 Hz f < 1 Hz Characteristics of MT method Infinite distance of source sounding site plane wave assumption, time invariance of the source simplifies analysis of the governing equations Frequency domain and wide frequency bands intermediate to deep investigation depth Wide range of applications regional scale geological studies/tectonics mineral, geothermal and oil exploration 4

5 EM induction ~ wave diffusion incident waves surface transmitted waves reflected waves spectral content, not propagation parameters MT field set-up 5

6 MT field set-up MT time series 6

7 MT time series MT time series 7

8 Electric (E) and magnetic (H) fields relationship For a homogeneous or layered (1-D) medium E x = Z H y Z = scalar impedance For a medium with 2-D 2 D symmetry E x = Z xy H y E y = Z yx H x Z xy Z yx For a general 3-D 3 D medium E x = Z xx H x + Z xy H y E = Z H E y = Z yx H x + Z yy H y Z = tensor impedance Data processing To extract impedance tensor Z from observed EM fields (time series of E and H) Spectral analysis and transfer function estimation Analysis of sub- surface properties contained in Z 8

9 MT framework EM induction theory relationship between electric (E) and magnetic (H) phenomena dependence of EM fields on electrical property (resistivity or conductivity) of the medium Observation of EM fields correlation between horizontal components of EM fields expressed as tensor impedance (Z) Extraction of subsurface parameters from Z EM theory (Maxwell s s equations) electric - magnetic phenomena represented by the Maxwell s s equations B(t) t E(t) E(t) = B(t) t 9

10 electric - magnetic phenomena represented by the Maxwell s s equations J(t) + D(t) t H(t) H(t) = J(t) + D(t) t Resolving the Maxwell s s equations Constitutive equations: B = µh, J = σe, D = εe Time dependency: exp(+iωt ), ω = 2π f = 2π/ T Neglecting electric displacement term (D) and variation of µ = µ 0 = 4π 10-7 H/m, ε = ε 0 Consider 1-D 1 D medium resistivity varies only with depth: ρ(z) no vertical components of EM fields horizontal components (E x, E y ) and (H x, H y ) vary with z (but do not vary with x and y) 10

11 component of ( E(t) = ) The x component of 2 E x z 2 = iωµ 0 σ E x The y component of E x z (1) component of ( ) = iωµ 0 H y E(t) = (2) B(t) Elementary solution to diffusion equation (1) t B(t) E x = A exp( kz) + B exp(+kz) k = (iωµ 0 /ρ) 1/2 t Elementary solution to Maxwell s s equations in 1-D1 E x = A exp ( k( z) + B exp (+k( z) k H y = (A exp ( k( z) B exp (+k( z)) i ωµ 0 EM fields as function of depth (z( > 0 downwards) terms with A represent attenuation of EM fields with increasing depth terms with B represent attenuation of EM fields with decreasing depth ("reflected" waves) A and B are constants to be determined from boundary conditions 11

12 Amplitude attenuation of EM fields attenuation with increasing depth (z( > 0) attenuation with decreasing depth (z( < 0) Homogeneous half-space medium (with external sources, i.e. MT) no resistivity interface at depth, i.e. no "reflection EM fields tend to zero at great depth no terms with B E x = A exp ( k( z) k H y = A exp ( k( z) i ωµ 0 attenuation with increasing depth (z( > 0) 12

13 Homogeneous half-space medium k = (iωµ 0 /ρ) 1/2 = α + i α α = ( 0.5 ωµ 0 /ρ) 1/2 E x = A exp ( α( z) ) exp ( iα( z) amplitude decay A(z) sinusoidal term amplitude of EM fields decays exponentially with depth and becomes negligible at certain depth Skin effect and penetration depth Skin effect = exponential EM wave attenuation with depth Skin depth (δ)( ) = depth in a homogeneous medium at which the amplitude becomes 1/e of that at the surface A(δ) = A exp ( α( δ) = A exp ( 1)( δ = (2ρ/ωµ 0 ) 1/2 500 (ρt ) 1/2 δ in meter, ρ in Ohm.m, Τ in seconds Skin depth is associated to penetration depth of EM 13

14 Skin effect and penetration depth δ 500 (ρt ) 1/2 Lower frequency (or higher period) and higher resistivity ~ slower attenuation ~ deeper penetration Principles of MT sounding, i.e. wide frequency band measurement probes different parts (depths) of the subsurface z Homogeneous half-space medium impedance, proportionality between E and H E x = A exp ( k( z) k H y = A exp ( k( z) i ωµ 0 Z xy E x H y i ωµ xy = = 0 = k (i ωµ 0 ρ) 1/2 Z 0 = intrinsic impedance 14

15 Homogeneous half-space medium intrinsic impedance Z 0 = (i ωµ 0 ρ) 1/2 Z 0 = ( 0.5 ωµ + i ( 0.5 ωµ 0 ρ) 1/2 0 ρ) 1/2 resistivity and impedance phase 1 ρ = Z 0 2 ωµ 0 φ = tan Im (Z 0 ) tan [ ] 1 Re (Z 0 ) Homogeneous half-space medium intrinsic impedance Z 0 = (i ωµ 0 ρ) 1/2 impedance (theoretical) homogeneous half-space resistivity and phase 1 ρ = Z 0 2 ωµ 0 φ = tan Im (Z 0 ) tan [ ] 1 Re (Z 0 ) resistivity (inferred) impedance (measured) 15

16 Measurement and data processing Measurement of orthogonal EM fields (time series) E x, E y, H x, H y Data processing to extract impedance tensor E x = Z xx E y = Z yx H x + Z xy H y H x + Z yy H y E = Z H Apparent resistivity and phase ρ a(ij ij) = 1 ωµ 0 Z ij 2 φ (ij ij) tan [ ] 1 = tan Im (Z ij ij ) Re (Z ij ij ) Apparent resistivity and phase sounding curves plot of log 10 ρ a vs. T and φ vs. T qualitatively represent resistivity as function of depth 16

17 Apparent resistivity and phase sounding curves APP. RESISTIVITY (Ohm.m) obs. data calc. data DEPTH (m) RESISTIVITY (Ohm.m) ? 1000 PHASE (deg.) PERIOD (sec.) Summary Electromagnetic induction concepts MT framework: Measured parameters model parameters relationship Maxwell s s equations: solution for simplest medium skin effect and penetration depth MT sounding principles MT raw data MT sounding curve 17

A Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods

A Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods A Brief Introduction to Magnetotellurics and Controlled Source Electromagnetic Methods Frank Morrison U.C. Berkeley With the help of: David Alumbaugh Erika Gasperikova Mike Hoversten Andrea Zirilli A few

More information

MT Prospecting. Map Resistivity. Determine Formations. Determine Structure. Targeted Drilling

MT Prospecting. Map Resistivity. Determine Formations. Determine Structure. Targeted Drilling MT Prospecting Map Resistivity Determine Formations Determine Structure Targeted Drilling Cross-sectional interpretation before and after an MT survey of a mineral exploration prospect containing volcanic

More information

Geophysics for Geothermal Exploration

Geophysics for Geothermal Exploration Geophysics for Geothermal Exploration Dr. Hendra Grandis Geofisika - ITB Agenda Basic concept of Geophysics: : Review Geophysical signatures of a geothermal system Geophysical methods for geotermal exploration

More information

Hendra Grandis 1 and Prihadi Sumintadireja 2

Hendra Grandis 1 and Prihadi Sumintadireja 2 A BRIEF REVIEW FOR THE PROPER APPLICATION OF MAGNETOTELLURIC (MT) AND CONTROLLED-SOURCE AUDIO-FREQUENCY MAGNETOTELLURIC (CSAMT) IN GEOTHERMAL EXPLORATION Hendra Grandis 1 and Prihadi Sumintadireja 2 ABSTRACT

More information

GLE 594: An introduction to applied geophysics

GLE 594: An introduction to applied geophysics GL 594: An introduction to applied geophysics Ground Penetrating Radar Fall 005 Ground Penetrating Radar Reading Today: 309-316 Next class: 316-39 Introduction to GPR Using the reflection (and sometimes

More information

Geophysical Applications GPR Ground Penetrating Radar

Geophysical Applications GPR Ground Penetrating Radar Overview: Basics of GPR Radar-wave velocity, attenuation and skin depth Modes of acquisition The Radar-range equation Dielectric properties of materials and relation to porosity Case studies [Archeology,

More information

THEORETICAL ANALYSIS AND NUMERICAL CALCULATION OF LOOP-SOURCE TRANSIENT ELECTROMAGNETIC IMAGING

THEORETICAL ANALYSIS AND NUMERICAL CALCULATION OF LOOP-SOURCE TRANSIENT ELECTROMAGNETIC IMAGING CHINESE JOURNAL OF GEOPHYSICS Vol.47, No.2, 24, pp: 379 386 THEORETICAL ANALYSIS AND NUMERICAL CALCULATION OF LOOP-SOURCE TRANSIENT ELECTROMAGNETIC IMAGING XUE Guo-Qiang LI Xiu SONG Jian-Ping GUO Wen-Bo

More information

Electromagnetic Waves For fast-varying phenomena, the displacement current cannot be neglected, and the full set of Maxwell s equations must be used

Electromagnetic Waves For fast-varying phenomena, the displacement current cannot be neglected, and the full set of Maxwell s equations must be used Electromagnetic Waves For fast-varying phenomena, the displacement current cannot be neglected, and the full set of Maxwell s equations must be used B( t) E = dt D t H = J+ t D =ρ B = 0 D=εE B=µ H () F

More information

5 Electromagnetic Waves

5 Electromagnetic Waves 5 Electromagnetic Waves 5.1 General Form for Electromagnetic Waves. In free space, Maxwell s equations are: E ρ ɛ 0 (5.1.1) E + B 0 (5.1.) B 0 (5.1.3) B µ 0 ɛ 0 E µ 0 J (5.1.4) In section 4.3 we derived

More information

General review: - a) Dot Product

General review: - a) Dot Product General review: - a) Dot Product If θ is the angle between the vectors a and b, then a b = a b cos θ NOTE: Two vectors a and b are orthogonal, if and only if a b = 0. Properties of the Dot Product If a,

More information

Introducing the Stratagem EH4. Electrical Conductivity Imaging System Hybrid-Source Magnetotellurics

Introducing the Stratagem EH4. Electrical Conductivity Imaging System Hybrid-Source Magnetotellurics Introducing the Stratagem EH4 Electrical Conductivity Imaging System Hybrid-Source Magnetotellurics Stratagem EH4 Hybrid-Source Magnetotellurics Frequency range of 10 Hz to 90k Hz Approx. depth of investigation

More information

Peter R. Milligan. Stephan Thiel. Geoscience Australia. The University of Adelaide. MT acquisition & processing

Peter R. Milligan. Stephan Thiel. Geoscience Australia. The University of Adelaide. MT acquisition & processing Magnetotelluric acquisition & processing, with examples from the Gawler Craton, Curnamona Province and Curnamona- Gawler Link transects in South Australia Peter R. Milligan Geoscience Australia Stephan

More information

Uniform Plane Waves Page 1. Uniform Plane Waves. 1 The Helmholtz Wave Equation

Uniform Plane Waves Page 1. Uniform Plane Waves. 1 The Helmholtz Wave Equation Uniform Plane Waves Page 1 Uniform Plane Waves 1 The Helmholtz Wave Equation Let s rewrite Maxwell s equations in terms of E and H exclusively. Let s assume the medium is lossless (σ = 0). Let s also assume

More information

Uniform Plane Waves. Ranga Rodrigo. University of Moratuwa. November 7, 2008

Uniform Plane Waves. Ranga Rodrigo. University of Moratuwa. November 7, 2008 Uniform Plane Waves Ranga Rodrigo University of Moratuwa November 7, 2008 Ranga Rodrigo (University of Moratuwa) Uniform Plane Waves November 7, 2008 1 / 51 Summary of Last Week s Lecture Basic Relations

More information

Plane Waves Part II. 1. For an electromagnetic wave incident from one medium to a second medium, total reflection takes place when

Plane Waves Part II. 1. For an electromagnetic wave incident from one medium to a second medium, total reflection takes place when Plane Waves Part II. For an electromagnetic wave incident from one medium to a second medium, total reflection takes place when (a) The angle of incidence is equal to the Brewster angle with E field perpendicular

More information

Towards a Dimensionally Adaptive Inversion of the Magnetotelluric Problem

Towards a Dimensionally Adaptive Inversion of the Magnetotelluric Problem 1 Towards a Dimensionally Adaptive Inversion of the Magnetotelluric Problem J. Alvarez-Aramberri 1 D. Pardo 2 H. Barucq 3 1 University of the Basque Country (UPV/EHU), Bilbao, Spain, Inria team-project

More information

1 Fundamentals of laser energy absorption

1 Fundamentals of laser energy absorption 1 Fundamentals of laser energy absorption 1.1 Classical electromagnetic-theory concepts 1.1.1 Electric and magnetic properties of materials Electric and magnetic fields can exert forces directly on atoms

More information

Module 5 : Plane Waves at Media Interface. Lecture 36 : Reflection & Refraction from Dielectric Interface (Contd.) Objectives

Module 5 : Plane Waves at Media Interface. Lecture 36 : Reflection & Refraction from Dielectric Interface (Contd.) Objectives Objectives In this course you will learn the following Reflection and Refraction with Parallel Polarization. Reflection and Refraction for Normal Incidence. Lossy Media Interface. Reflection and Refraction

More information

M.Sc. Thesis. Author: Merijn de Block

M.Sc. Thesis. Author: Merijn de Block M.Sc. Thesis An alternative method for estimating the porosity and the permeability of (potential) geothermal reservoirs using magnetotelluric data; A case study of the Sherwood Sandstone Group in the

More information

ECE 107: Electromagnetism

ECE 107: Electromagnetism ECE 107: Electromagnetism Notes Set 1 Instructor: Prof. Vitaliy Lomakin Department of Electrical and Computer Engineering University of California, San Diego, CA 92093 1 Introduction (1) atom Electromagnetism

More information

Electromagnetic Waves

Electromagnetic Waves Electromagnetic Waves Maxwell s equations predict the propagation of electromagnetic energy away from time-varying sources (current and charge) in the form of waves. Consider a linear, homogeneous, isotropic

More information

Lesson 1 Introduction to geophysical methods Emanuele Forte

Lesson 1 Introduction to geophysical methods Emanuele Forte Lesson 1 Introduction to geophysical methods Emanuele Forte 3 rd September 2016 Naxos, Greece 1 Outline - What is Geophysics - Base concepts: (Geo)Physical parameters Sensitivity Resolution - Active and

More information

EECS 117 Lecture 20: Plane Waves

EECS 117 Lecture 20: Plane Waves University of California, Berkeley EECS 117 Lecture 20 p. 1/2 EECS 117 Lecture 20: Plane Waves Prof. Niknejad University of California, Berkeley University of California, Berkeley EECS 117 Lecture 20 p.

More information

3. Magnetic Methods / 62

3. Magnetic Methods / 62 Contents Preface to the Second Edition / xv Excerpts from Preface to the FirstEdition / xvii Mathematical Conventions / xix 1. Introduction / 1 Reference / 5 2. Gravity Methods / 6 2. I. Introduction /

More information

Kimmo Silvonen, Transmission lines, ver

Kimmo Silvonen, Transmission lines, ver Kimmo Silvonen, Transmission lines, ver. 13.10.2008 1 1 Basic Theory The increasing operating and clock frequencies require transmission line theory to be considered more and more often! 1.1 Some practical

More information

E : Ground-penetrating radar (GPR)

E : Ground-penetrating radar (GPR) Geophysics 3 March 009 E : Ground-penetrating radar (GPR) The EM methods in section D use low frequency signals that trael in the Earth by diffusion. These methods can image resistiity of the Earth on

More information

3 Constitutive Relations: Macroscopic Properties of Matter

3 Constitutive Relations: Macroscopic Properties of Matter EECS 53 Lecture 3 c Kamal Sarabandi Fall 21 All rights reserved 3 Constitutive Relations: Macroscopic Properties of Matter As shown previously, out of the four Maxwell s equations only the Faraday s and

More information

Geophysics for Environmental and Geotechnical Applications

Geophysics for Environmental and Geotechnical Applications Geophysics for Environmental and Geotechnical Applications Dr. Katherine Grote University of Wisconsin Eau Claire Why Use Geophysics? Improve the quality of site characterization (higher resolution and

More information

INTRODUCTION TO APPLIED GEOPHYSICS

INTRODUCTION TO APPLIED GEOPHYSICS INTRODUCTION TO APPLIED GEOPHYSICS EXPLORING THE SHALL0W SUBSURFACE H. Robert Burger Anne F. Sheehan Craig H.Jones VERSITY OF COLORADO VERSITY OF COLORADO W. W. NORTON & COMPANY NEW YORK LONDON Contents

More information

DHANALAKSHMI SRINIVASAN INSTITUTE OF RESEARCH AND TECHNOLOGY

DHANALAKSHMI SRINIVASAN INSTITUTE OF RESEARCH AND TECHNOLOGY DHANALAKSHMI SRINIVASAN INSTITUTE OF RESEARCH AND TECHNOLOGY SIRUVACHUR-621113 ELECTRICAL AND ELECTRONICS DEPARTMENT 2 MARK QUESTIONS AND ANSWERS SUBJECT CODE: EE 6302 SUBJECT NAME: ELECTROMAGNETIC THEORY

More information

E E D E=0 2 E 2 E (3.1)

E E D E=0 2 E 2 E (3.1) Chapter 3 Constitutive Relations Maxwell s equations define the fields that are generated by currents and charges. However, they do not describe how these currents and charges are generated. Thus, to find

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION Supplementary online material for Bai et al., (2). EHS3D MT data collection Broadband magnetotelluric (MT) data were recorded on profiles P, P2 and P4 in the frequency band -.5

More information

Numerical simulation of magnetotelluric fields at Stromboli

Numerical simulation of magnetotelluric fields at Stromboli Numerical simulation of magnetotelluric fields at Stromboli A. Franke, S. Kütter, R.-U. Börner and K. Spitzer TU Bergakademie Freiberg, Germany 1 Summary Stromboli is a small volcanic island in the Mediterranean

More information

PART A: Short-answer questions (50%; each worth 2%)

PART A: Short-answer questions (50%; each worth 2%) PART A: Short-answer questions (50%; each worth 2%) Your answers should be brief (just a few words) and may be written on these pages if you wish. Remember to hand these pages in with your other exam pages!

More information

Simple medium: D = ɛe Dispersive medium: D = ɛ(ω)e Anisotropic medium: Permittivity as a tensor

Simple medium: D = ɛe Dispersive medium: D = ɛ(ω)e Anisotropic medium: Permittivity as a tensor Plane Waves 1 Review dielectrics 2 Plane waves in the time domain 3 Plane waves in the frequency domain 4 Plane waves in lossy and dispersive media 5 Phase and group velocity 6 Wave polarization Levis,

More information

RESISTIVITY SURVEYING AND ELECTROMAGNETIC METHODS

RESISTIVITY SURVEYING AND ELECTROMAGNETIC METHODS Presented at Short Course VII on Surface Exploration for Geothermal Resources, organized by UNU-GTP and LaGeo, in Santa Tecla and Ahuachapán, El Salvador, March 14-22, 2015. GEOTHERMAL TRAINING PROGRAMME

More information

ELE3310: Basic ElectroMagnetic Theory

ELE3310: Basic ElectroMagnetic Theory A summary for the final examination EE Department The Chinese University of Hong Kong November 2008 Outline Mathematics 1 Mathematics Vectors and products Differential operators Integrals 2 Integral expressions

More information

Chapter 9. Electromagnetic waves

Chapter 9. Electromagnetic waves Chapter 9. lectromagnetic waves 9.1.1 The (classical or Mechanical) waves equation Given the initial shape of the string, what is the subsequent form, The displacement at point z, at the later time t,

More information

Electromagnetically Induced Flows in Water

Electromagnetically Induced Flows in Water Electromagnetically Induced Flows in Water Michiel de Reus 8 maart 213 () Electromagnetically Induced Flows 1 / 56 Outline 1 Introduction 2 Maxwell equations Complex Maxwell equations 3 Gaussian sources

More information

a. The distance from a shot point in which no refracted energy is recorded is termed the.

a. The distance from a shot point in which no refracted energy is recorded is termed the. GE 3040 Midterm #2 4/08/2009 Name 1. Fill in the blank: (20 pts) a. The distance from a shot point in which no refracted energy is recorded is termed the. b. In a two-layer case, the upward moving wave

More information

Principles of Applied Geophysics

Principles of Applied Geophysics Principles of Applied Geophysics Fifth edition D.S. Parasnis Professor emeritus of Applied Geophysics Department of Applied Geophysics, University ofluled, Lulea, Sweden Fellow of the Royal Swedish Academy

More information

Chapter 7. Time-Varying Fields and Maxwell s Equations

Chapter 7. Time-Varying Fields and Maxwell s Equations Chapter 7. Time-arying Fields and Maxwell s Equations Electrostatic & Time arying Fields Electrostatic fields E, D B, H =J D H 1 E B In the electrostatic model, electric field and magnetic fields are not

More information

Electromagnetic Waves Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space

Electromagnetic Waves Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space Electromagnetic Waves 1 1. Retarded potentials 2. Energy and the Poynting vector 3. Wave equations for E and B 4. Plane EM waves in free space 1 Retarded Potentials For volume charge & current = 1 4πε

More information

KINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING QUESTION BANK

KINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING QUESTION BANK KINGS COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING QUESTION BANK SUB.NAME : ELECTROMAGNETIC FIELDS SUBJECT CODE : EC 2253 YEAR / SEMESTER : II / IV UNIT- I - STATIC ELECTRIC

More information

1 Macroscopic Maxwell s equations

1 Macroscopic Maxwell s equations This lecture purports to the macroscopic Maxwell s equations in the differential forms and their revealation about the propagation of light in vacuum and in matter. Microscopic Maxwell s equations and

More information

CHAPTER 9 ELECTROMAGNETIC WAVES

CHAPTER 9 ELECTROMAGNETIC WAVES CHAPTER 9 ELECTROMAGNETIC WAVES Outlines 1. Waves in one dimension 2. Electromagnetic Waves in Vacuum 3. Electromagnetic waves in Matter 4. Absorption and Dispersion 5. Guided Waves 2 Skip 9.1.1 and 9.1.2

More information

ELE 3310 Tutorial 10. Maxwell s Equations & Plane Waves

ELE 3310 Tutorial 10. Maxwell s Equations & Plane Waves ELE 3310 Tutorial 10 Mawell s Equations & Plane Waves Mawell s Equations Differential Form Integral Form Faraday s law Ampere s law Gauss s law No isolated magnetic charge E H D B B D J + ρ 0 C C E r dl

More information

POTASH DRAGON CHILE GEOPHYSICAL SURVEY TRANSIENT ELECTROMAGNETIC (TEM) METHOD. LLAMARA and SOLIDA PROJECTS SALAR DE LLAMARA, IQUIQUE, REGION I, CHILE

POTASH DRAGON CHILE GEOPHYSICAL SURVEY TRANSIENT ELECTROMAGNETIC (TEM) METHOD. LLAMARA and SOLIDA PROJECTS SALAR DE LLAMARA, IQUIQUE, REGION I, CHILE POTASH DRAGON CHILE GEOPHYSICAL SURVEY TRANSIENT ELECTROMAGNETIC (TEM) METHOD LLAMARA and SOLIDA PROJECTS SALAR DE LLAMARA, IQUIQUE, REGION I, CHILE OCTOBER 2012 CONTENT Page I INTRODUCTION 1 II FIELD

More information

Chapter 7. Time-Varying Fields and Maxwell s Equation

Chapter 7. Time-Varying Fields and Maxwell s Equation Chapter 7. Time-Varying Fields and Maxwell s Equation Electrostatic & Time Varying Fields Electrostatic fields E, D B, H =J D H 1 E B In the electrostatic model, electric field and magnetic fields are

More information

Class 15 : Electromagnetic Waves

Class 15 : Electromagnetic Waves Class 15 : Electromagnetic Waves Wave equations Why do electromagnetic waves arise? What are their properties? How do they transport energy from place to place? Recap (1) In a region of space containing

More information

Integrated Geophysical Model for Suswa Geothermal Prospect using Resistivity, Seismics and Gravity Survey Data in Kenya

Integrated Geophysical Model for Suswa Geothermal Prospect using Resistivity, Seismics and Gravity Survey Data in Kenya Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 Integrated Geophysical Model for Suswa Geothermal Prospect using Resistivity, Seismics and Gravity Survey Data in Kenya

More information

Inductive source induced polarization David Marchant, Eldad Haber and Douglas W. Oldenburg, University of British Columbia

Inductive source induced polarization David Marchant, Eldad Haber and Douglas W. Oldenburg, University of British Columbia David Marchant, Eldad Haber and Douglas W. Oldenburg, University of British Columbia Downloaded /9/3 to 37.8.7.. Redistribution subject to SEG license or copyright; see Terms of Use at http://library.seg.org/

More information

Theoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 9

Theoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 9 WiSe 202 20.2.202 Prof. Dr. A-S. Smith Dipl.-Phys. Ellen Fischermeier Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg

More information

Attenuation of Multiples in Marine Seismic Data. Geology /23/06

Attenuation of Multiples in Marine Seismic Data. Geology /23/06 Attenuation of Multiples in Marine Seismic Data Geology 377 10/23/06 Marine Seismic Surveying Marine seismic surveying is performed to support marine environmental and civil engineering projects in the

More information

ambiguity in earth sciences IESO Geophysics Section Eddy hartantyo, Lab Geofisika FMIPA UGM

ambiguity in earth sciences IESO Geophysics Section Eddy hartantyo, Lab Geofisika FMIPA UGM ambiguity in earth sciences IESO Geophysics Section Eddy hartantyo, Lab Geofisika FMIPA UGM Pelatihan Tahap II IESO Teknik Geologi UGM Februari 2009 1 Introduction Photos from http://www.eegs.org/whatis/

More information

Geophysics Course Introduction to DC Resistivity

Geophysics Course Introduction to DC Resistivity NORAD supported project in MRRD covering Capacity Building and Institutional Cooperation in the field of Hydrogeology for Faryab Province Afghanistan Geophysics Course Introduction to DC Resistivity By

More information

Mapping of Mesozoic Sediments below Deccan trap using Magnetotelluric Method.

Mapping of Mesozoic Sediments below Deccan trap using Magnetotelluric Method. P - 364 Mapping of Mesozoic s below Deccan trap using Magnetotelluric Method. Surajit Gorain*(sgorain@gmail.com), Puja Prakash, Dr. Prabhakar Thakur Hindustan Petroleum Corporation Limited Summary Magnetotelluric

More information

Physics 106a/196a Problem Set 7 Due Dec 2, 2005

Physics 106a/196a Problem Set 7 Due Dec 2, 2005 Physics 06a/96a Problem Set 7 Due Dec, 005 Version 3, Nov 7, 005 In this set we finish up the SHO and study coupled oscillations/normal modes and waves. Problems,, and 3 are for 06a students only, 4, 5,

More information

Pseudo-seismic wavelet transformation of transient electromagnetic response in engineering geology exploration

Pseudo-seismic wavelet transformation of transient electromagnetic response in engineering geology exploration GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L645, doi:.29/27gl36, 27 Pseudo-seismic wavelet transformation of transient electromagnetic response in engineering geology exploration G. Q. Xue, Y. J. Yan, 2 and

More information

Wave Phenomena Physics 15c. Lecture 15 Reflection and Refraction

Wave Phenomena Physics 15c. Lecture 15 Reflection and Refraction Wave Phenomena Physics 15c Lecture 15 Reflection and Refraction What We (OK, Brian) Did Last Time Discussed EM waves in vacuum and in matter Maxwell s equations Wave equation Plane waves E t = c E B t

More information

1. Resistivity of rocks

1. Resistivity of rocks RESISTIVITY 1) Resistivity of rocks 2) General principles of resistivity surveying 3) Field procedures, interpretation and examples 4) Summary and conclusions INDUCED POLARIZATION 1) General principles

More information

Some Remarks on Shielding. Herbert Kapitza (FLA) (using slides from a talk by Mike Thuot) DESY,

Some Remarks on Shielding. Herbert Kapitza (FLA) (using slides from a talk by Mike Thuot) DESY, Some Remarks on Shielding Herbert Kapitza (FLA) (using slides from a talk by Mike Thuot) DESY, 09.10.2006 A shield may be used to confine the radiated field from a noise source. Shields are metallic partitions

More information

ECE357H1F ELECTROMAGNETIC FIELDS FINAL EXAM. 28 April Examiner: Prof. Sean V. Hum. Duration: hours

ECE357H1F ELECTROMAGNETIC FIELDS FINAL EXAM. 28 April Examiner: Prof. Sean V. Hum. Duration: hours UNIVERSITY OF TORONTO FACULTY OF APPLIED SCIENCE AND ENGINEERING The Edward S. Rogers Sr. Department of Electrical and Computer Engineering ECE357H1F ELECTROMAGNETIC FIELDS FINAL EXAM 28 April 15 Examiner:

More information

toroidal iron core compass switch battery secondary coil primary coil

toroidal iron core compass switch battery secondary coil primary coil Fundamental Laws of Electrostatics Integral form Differential form d l C S E 0 E 0 D d s V q ev dv D ε E D qev 1 Fundamental Laws of Magnetostatics Integral form Differential form C S dl S J d s B d s

More information

An Experimental Magnetotelluric Short Sounding During High Sun Spot Activity Cycle in 2013

An Experimental Magnetotelluric Short Sounding During High Sun Spot Activity Cycle in 2013 Proceedings World Geothermal Congress 2015 Melbourne, Australia, 19-25 April 2015 An Experimental Magnetotelluric Short Sounding During High Sun Spot Activity Cycle in 2013 Lendriadi AGUNG and Imam B.

More information

Site Characterization & Hydrogeophysics

Site Characterization & Hydrogeophysics Site Characterization & Hydrogeophysics (Source: Matthew Becker, California State University) Site Characterization Definition: quantitative description of the hydraulic, geologic, and chemical properties

More information

Estimation of Cole-Cole parameters from time-domain electromagnetic data Laurens Beran and Douglas Oldenburg, University of British Columbia.

Estimation of Cole-Cole parameters from time-domain electromagnetic data Laurens Beran and Douglas Oldenburg, University of British Columbia. Estimation of Cole-Cole parameters from time-domain electromagnetic data Laurens Beran and Douglas Oldenburg, University of British Columbia. SUMMARY We present algorithms for the inversion of time-domain

More information

Geophysical mapping and imaging of soil structures: basic overview

Geophysical mapping and imaging of soil structures: basic overview Geophysical mapping and imaging of soil structures: basic overview Stéphane Garambois Landslide team, LGIT, Université Joseph Fourier 3D imaging of the water table (seismic) Geophysical methods 1 Plan

More information

Electromagnetic Waves

Electromagnetic Waves Physics 8 Electromagnetic Waves Overview. The most remarkable conclusion of Maxwell s work on electromagnetism in the 860 s was that waves could exist in the fields themselves, traveling with the speed

More information

TECHNO INDIA BATANAGAR

TECHNO INDIA BATANAGAR TECHNO INDIA BATANAGAR ( DEPARTMENT OF ELECTRONICS & COMMUNICATION ENGINEERING) QUESTION BANK- 2018 1.Vector Calculus Assistant Professor 9432183958.mukherjee@tib.edu.in 1. When the operator operates on

More information

1 Electromagnetic concepts useful for radar applications

1 Electromagnetic concepts useful for radar applications Electromagnetic concepts useful for radar applications The scattering of electromagnetic waves by precipitation particles and their propagation through precipitation media are of fundamental importance

More information

Overview in Images. S. Lin et al, Nature, vol. 394, p , (1998) T.Thio et al., Optics Letters 26, (2001).

Overview in Images. S. Lin et al, Nature, vol. 394, p , (1998) T.Thio et al., Optics Letters 26, (2001). Overview in Images 5 nm K.S. Min et al. PhD Thesis K.V. Vahala et al, Phys. Rev. Lett, 85, p.74 (000) J. D. Joannopoulos, et al, Nature, vol.386, p.143-9 (1997) T.Thio et al., Optics Letters 6, 197-1974

More information

Basic principles of the seismic method

Basic principles of the seismic method Chapter 2 Basic principles of the seismic method In this chapter we introduce the basic notion of seismic waves. In the earth, seismic waves can propagate as longitudinal (P) or as shear (S) waves. For

More information

INSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad Electronics and Communicaton Engineering

INSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad Electronics and Communicaton Engineering INSTITUTE OF AERONAUTICAL ENGINEERING Dundigal, Hyderabad - 00 04 Electronics and Communicaton Engineering Question Bank Course Name : Electromagnetic Theory and Transmission Lines (EMTL) Course Code :

More information

Magnetic and Gravity Methods for Geothermal Exploration

Magnetic and Gravity Methods for Geothermal Exploration Magnetic and Gravity Methods for Geothermal Exploration Dr. Hendra Grandis Geophysics - ITB method and survey procedure Aero- or ground magnetic (covers a large area) Schlumberger resistivity mapping and

More information

TIME DOMAIN ELECTROMAGNETICS IN MARS ANALOG ENVIRONMENTS: COMPARING TWO FIELD STUDIES. Joern A. Jernsletten 1

TIME DOMAIN ELECTROMAGNETICS IN MARS ANALOG ENVIRONMENTS: COMPARING TWO FIELD STUDIES. Joern A. Jernsletten 1 AAS 06-256 TIME DOMAIN ELECTROMAGNETICS IN MARS ANALOG ENVIRONMENTS: COMPARING TWO FIELD STUDIES Joern A. Jernsletten 1 The purpose of this study is to evaluate the use of (diffusive) Time Domain Electromagnetics

More information

Early Analogue Modeling Experiments and Related Studies to Today s Problems of Geo-electromagnetic Exploration

Early Analogue Modeling Experiments and Related Studies to Today s Problems of Geo-electromagnetic Exploration Journal of Earth Science, Vol. 20, No. 3, p. 618 625, June 2009 ISSN 1674-487X Printed in China DOI: 10.1007/s12583-009-0051-y Early Analogue Modeling Experiments and Related Studies to Today s Problems

More information

Geophysics of Exploration for Water. Geoelectric methods III. edited by Peter Vass

Geophysics of Exploration for Water. Geoelectric methods III. edited by Peter Vass Geophysics of Exploration for Water Geoelectric methods III edited by Peter Vass The phenomenon of induced polarization was already noticed at the early stage of applying direct current electrical methods

More information

Research Article Error Analysis in Measured Conductivity under Low Induction Number Approximation for Electromagnetic Methods

Research Article Error Analysis in Measured Conductivity under Low Induction Number Approximation for Electromagnetic Methods ISRN Geophysics Volume 2013, Article ID 720839, 4 pages http://dx.doi.org/10.1155/2013/720839 Research Article Error Analysis in Measured Conductivity under Low Induction Number Approximation for Electromagnetic

More information

Plane Waves GATE Problems (Part I)

Plane Waves GATE Problems (Part I) Plane Waves GATE Problems (Part I). A plane electromagnetic wave traveling along the + z direction, has its electric field given by E x = cos(ωt) and E y = cos(ω + 90 0 ) the wave is (a) linearly polarized

More information

Eddy Current Testing using the Bode 100

Eddy Current Testing using the Bode 100 Page 1 of 12 using the Bode 100 Lukas Heinzle Abstract: (ET) is a commonly used technique surface inspections of conducting materials. An eddy current sensor, namely a probe coil that produces an alternating

More information

Waves in Linear Optical Media

Waves in Linear Optical Media 1/53 Waves in Linear Optical Media Sergey A. Ponomarenko Dalhousie University c 2009 S. A. Ponomarenko Outline Plane waves in free space. Polarization. Plane waves in linear lossy media. Dispersion relations

More information

Lecture 10. Surface Energy Balance (Garratt )

Lecture 10. Surface Energy Balance (Garratt ) Lecture 10. Surface Energy Balance (Garratt 5.1-5.2) The balance of energy at the earth s surface is inextricably linked to the overlying atmospheric boundary layer. In this lecture, we consider the energy

More information

Physics 54 Lecture March 1, Micro-quiz problems (magnetic fields and forces) Magnetic dipoles and their interaction with magnetic fields

Physics 54 Lecture March 1, Micro-quiz problems (magnetic fields and forces) Magnetic dipoles and their interaction with magnetic fields Physics 54 Lecture March 1, 2012 OUTLINE Micro-quiz problems (magnetic fields and forces) Magnetic dipoles and their interaction with magnetic fields Electromagnetic induction Introduction to electromagnetic

More information

UNIT I ELECTROSTATIC FIELDS

UNIT I ELECTROSTATIC FIELDS UNIT I ELECTROSTATIC FIELDS 1) Define electric potential and potential difference. 2) Name few applications of gauss law in electrostatics. 3) State point form of Ohm s Law. 4) State Divergence Theorem.

More information

I) Resistivity. Literatur: Telford, Geldart und Sheriff (1990): Applied Geophysics, Second Edition, Cambridge University Press, NY.

I) Resistivity. Literatur: Telford, Geldart und Sheriff (1990): Applied Geophysics, Second Edition, Cambridge University Press, NY. I) Resistivity Literatur: Telford, Geldart und Sheriff (1990): Applied Geophysics, Second Edition, Cambridge University Press, NY. Bender, F. (Hg.) (1985): Angewandte Geowissenschaften, Bd. II, Enke Verlag,

More information

Hydrocarbon Fluid Wave and it's resonance G.K. Batta*, D.S. Manral, G.V.J. Rao, and S.K. Basha, Oil India limited, India

Hydrocarbon Fluid Wave and it's resonance G.K. Batta*, D.S. Manral, G.V.J. Rao, and S.K. Basha, Oil India limited, India Hydrocarbon Fluid Wave and it's resonance G.K. Batta*, D.S. Manral, G.V.J. Rao, and S.K. Basha, Oil India limited, India gnangp@gmail.com Keywords Hydrocarbon fluid wave, seismic acoustic impedance, law

More information

2 Decay of electromagnetic fields in the ocean

2 Decay of electromagnetic fields in the ocean Some notes on bathymetric effects in marine magnetotellurics, motivated by an amphibious experiment at the South Chilean margin Gerhard Kapinos, Heinrich Brasse 1 Freie Universität Berlin, Fachrichtung

More information

ECE357H1S ELECTROMAGNETIC FIELDS TERM TEST March 2016, 18:00 19:00. Examiner: Prof. Sean V. Hum

ECE357H1S ELECTROMAGNETIC FIELDS TERM TEST March 2016, 18:00 19:00. Examiner: Prof. Sean V. Hum UNIVERSITY OF TORONTO FACULTY OF APPLIED SCIENCE AND ENGINEERING The Edward S. Rogers Sr. Department of Electrical and Computer Engineering ECE357H1S ELECTROMAGNETIC FIELDS TERM TEST 2 21 March 2016, 18:00

More information

Worked Examples Set 2

Worked Examples Set 2 Worked Examples Set 2 Q.1. Application of Maxwell s eqns. [Griffiths Problem 7.42] In a perfect conductor the conductivity σ is infinite, so from Ohm s law J = σe, E = 0. Any net charge must be on the

More information

Sensibility Analysis of Inductance Involving an E-core Magnetic Circuit for Non Homogeneous Material

Sensibility Analysis of Inductance Involving an E-core Magnetic Circuit for Non Homogeneous Material Sensibility Analysis of Inductance Involving an E-core Magnetic Circuit for Non Homogeneous Material K. Z. Gomes *1, T. A. G. Tolosa 1, E. V. S. Pouzada 1 1 Mauá Institute of Technology, São Caetano do

More information

Antennas and Propagation

Antennas and Propagation Antennas and Propagation Ranga Rodrigo University of Moratuwa October 20, 2008 Compiled based on Lectures of Prof. (Mrs.) Indra Dayawansa. Ranga Rodrigo (University of Moratuwa) Antennas and Propagation

More information

Acoustic Velocity, Impedance, Reflection, Transmission, Attenuation, and Acoustic Etalons

Acoustic Velocity, Impedance, Reflection, Transmission, Attenuation, and Acoustic Etalons Acoustic Velocity, Impedance, Reflection, Transmission, Attenuation, and Acoustic Etalons Acoustic Velocity The equation of motion in a solid is (1) T = ρ 2 u t 2 (1) where T is the stress tensor, ρ is

More information

Problem set 3. Electromagnetic waves

Problem set 3. Electromagnetic waves Second Year Electromagnetism Michaelmas Term 2017 Caroline Terquem Problem set 3 Electromagnetic waves Problem 1: Poynting vector and resistance heating This problem is not about waves but is useful to

More information

The Interaction of Light and Matter: α and n

The Interaction of Light and Matter: α and n The Interaction of Light and Matter: α and n The interaction of light and matter is what makes life interesting. Everything we see is the result of this interaction. Why is light absorbed or transmitted

More information

Chapter 2 Basics of Electricity and Magnetism

Chapter 2 Basics of Electricity and Magnetism Chapter 2 Basics of Electricity and Magnetism My direct path to the special theory of relativity was mainly determined by the conviction that the electromotive force induced in a conductor moving in a

More information

Application of sensitivity analysis in DC resistivity monitoring of SAGD steam chambers

Application of sensitivity analysis in DC resistivity monitoring of SAGD steam chambers The University of British Columbia Geophysical Inversion Facility Application of sensitivity analysis in DC resistivity monitoring of SAGD steam chambers S. G. R. Devriese and D. W. Oldenburg gif.eos.ubc.ca

More information

Time-harmonic form Phasor form. =0 (1.11d)

Time-harmonic form Phasor form. =0 (1.11d) Chapter 2 Wave in an Unbounded Medium Maxwell s Equations Time-harmonic form Phasor form (Real quantity) (complex quantity) B E = Eˆ = jωbˆ (1.11 a) t D= ρ Dˆ = ρ (1.11 b) D H = J + Hˆ = Jˆ+ jωdˆ ( 1.11

More information

An Introduction to Geophysical Exploration

An Introduction to Geophysical Exploration An Introduction to Geophysical Exploration Philip Kearey Department of Earth Sciences University of Bristol Michael Brooks Ty Newydd, City Near Cowbridge Vale of Glamorgan Ian Hill Department of Geology

More information

APPENDIX A: Magnetotelluric Data in Relation to San Pedro Mesa Structural. The San Pedro Mesa structural high (discussed in main text of paper) was

APPENDIX A: Magnetotelluric Data in Relation to San Pedro Mesa Structural. The San Pedro Mesa structural high (discussed in main text of paper) was Page of DR for GSA Special Paper 9, Chapter, Geophysical constraints APPENDIX A: Magnetotelluric Data in Relation to San Pedro Mesa Structural High The San Pedro Mesa structural high (discussed in main

More information