P. Schmoldt, PhD - MTNet - DIAS
P. Schmoldt, PhD - MTNet - DIAS P. Schmoldt, PhD - MTNet - DIAS
3. Mathematical description of electromagnetic relations Apparent resistivity and impedance phase are therefore related to electric impedance via ρai j 1 = µω Zi j 2 where Zi j represents the integrated impedance of the (heterogeneous) subsurface, and respectively, with i, j ∈ [x, y]. φi j = arctan Im(Zi j) Re(Zi j) 3.2.2. Relationships for simple subsurface cases 1D case (3.35) (3.36) For situations where the subsurface exhibits solely vertical conductivity changes, naturally no alignment of any coordinate system axis with lateral interfaces can be made, making the choice of the coordinate system direction completely arbitrary. In such cases, diagonal elements of the impedance tensor (Eq. 3.34) are zero since electric fields are independent of parallel magnetic fields: Ex Ey = 0 Zxy −Zxy 0 Hx Hy . (3.37) The off-diagonal elements are of the same amplitude, but with an inverted sign due to the change from a right-handed into a left-handed coordinate system when dealing with the yx-component 2D case When the subsurface possesses lateral conductivity interfaces in only one direction (with an arbitrary number of vertical conductivity changes), the original coordinate system used for the recording can be rotated to the interface directions using a transformation matrix cos(Θ) sin(Θ) R = . (3.38) − sin(Θ) cos(Θ) For the case of perfect adjustment the impedance matrix reduces to Ex 0 = Zyx Zxy Hx . 0 (3.39) 36 Ey Hy
3.3. Magnetotelluric induction area Induction arrows Magnitude Wiese Direction Parkinson Real arrow ||Re(Tx) + Re(Ty)|| arctan Re(Ty) Re(Tx) arctan − Re(Ty) Imaginary arrow ||Im(Tx) + Im(Ty)|| arctan Re(Tx) Im(Ty) Im(Tx) arctan − Im(Ty) Im(Tx) Tab. 3.1.: Induction arrows derived from the magnetic transfer function (T = (Tx, Ty)) with two different conventions regarding the direction, pointing towards the resistor (Wiese [Wiese, 1962]) and the conductor (Parkinson [Parkinson, 1959]), respectively. 3.2.3. Vertical magnetic transfer function The vertical magnetic transfer function (T = (Tx, Ty)), also referred to as Tipper, relates the vertical magnetic field (Hz) to the horizontal magnetic fields ( Hh = (Hx, Hy)) Hz( f ) = T( f ) Hh( f ). (3.40) For the case of a 2D subsurface and adequate rotation of the coordinate system Equation 3.40 reduces to Hz( f ) = Ty( f )Hy( f ) (3.41) with the y-axis parallel to the conductivity interface [e.g. Vozoff , 1987]. The magnetic transfer function can be displayed using the induction arrows defined in Table 3.1, with arrows pointing towards either resistive regions (Wiese convention [Wiese, 1962]) or conductive regions (Parkinson convention [Parkinson, 1959]). For the undisturbed 2D case the phases of Tx and Ty are equal, thus the magnitude of the imaginary induction arrow equals zero, and the direction of the real induction arrow is orthogonal to the geoelectric strike direction. 3.3. Magnetotelluric induction area In MT, the depth of investigation is directly dependent on period range and subsurface conductivity, as the intensity of a penetrating magnetic wave is proportional to each of these. As a measure of the MT method’s sensitivity, the distance to the Earth’s surface at which the amplitude of the penetrating wave is reduced by 1/e is commonly used, i.e. F(δs) = 1 e F0 = F0e −1 = F0e −Re(k)δs . (3.42) Therein, δs is referred to as skin depth (or induction depth). Using the expression for the magnetic field in Equation 3.22, the wave number as defined in Equation 3.21, and Re(k) = ωµσ/2, (3.43) 37
- Page 22 and 23: List of Acronyms FE finite element
- Page 25 and 26: List of Symbols Below is a list of
- Page 27 and 28: Symbol SI unit Denotation φ · pha
- Page 29: Abstract The Tajo Basin and Betic C
- Page 32 and 33: Publications Poster presentations x
- Page 34 and 35: Acknowledgements Team, namely Colin
- Page 37 and 38: Introduction 1 The Iberian Peninsul
- Page 39 and 40: ections from enhanced one-dimension
- Page 41: Part I Theoretical background of ma
- Page 44 and 45: 2. Sources for magnetotelluric reco
- Page 46 and 47: 2. Sources for magnetotelluric reco
- Page 48 and 49: 2. Sources for magnetotelluric reco
- Page 50 and 51: 2. Sources for magnetotelluric reco
- Page 52 and 53: 2. Sources for magnetotelluric reco
- Page 54 and 55: 2. Sources for magnetotelluric reco
- Page 56 and 57: 2. Sources for magnetotelluric reco
- Page 58 and 59: 2. Sources for magnetotelluric reco
- Page 60 and 61: 2. Sources for magnetotelluric reco
- Page 62 and 63: 2. Sources for magnetotelluric reco
- Page 64 and 65: 2. Sources for magnetotelluric reco
- Page 67 and 68: Mathematical description of electro
- Page 69 and 70: yields 3.2. Deriving magnetotelluri
- Page 71: 3.2. Deriving magnetotelluric param
- Page 75 and 76: Depth d s d 1 d 2 d n-2 d n-1 t 1 t
- Page 77 and 78: 3.4. Boundary conditions materials
- Page 79 and 80: 3.5. The influence of electric perm
- Page 81 and 82: 3.5. The influence of electric perm
- Page 83 and 84: 3.5. The influence of electric perm
- Page 85 and 86: Distortion of magnetotelluric data
- Page 87 and 88: 4.1. Types of distortion Fig. 4.1.:
- Page 89 and 90: 4.1. Types of distortion Fig. 4.3.:
- Page 91 and 92: J s 0 s 0 4.1. Types of distortion
- Page 93 and 94: 4.1. Types of distortion Fig. 4.7.:
- Page 95 and 96: Scale Type Terminology Example Atom
- Page 97 and 98: 4.1. Types of distortion the use of
- Page 99 and 100: 4.2. Dimensionality Fig. 4.12.: The
- Page 101 and 102: 1D 2D local 3D/1D 3D/2D regional 4.
- Page 103 and 104: 4.3. General mathematical represent
- Page 105 and 106: 4.4. Removal of distortion effects
- Page 107 and 108: Parameter Geoelectrical application
- Page 109 and 110: 4.4. Removal of distortion effects
- Page 111 and 112: 4.4.5. Caldwell-Bibby-Brown phase t
- Page 113 and 114: 4.4. Removal of distortion effects
- Page 115: Method Applicability Swift angle 2D
- Page 118 and 119: 5. Earth’s properties observable
- Page 120 and 121: 5. Earth’s properties observable
3. Mathematical description of electromagnetic relations<br />
Apparent resistivity and impedance phase are therefore related to electric impedance via<br />
ρai j<br />
1 <br />
=<br />
<br />
µω<br />
<br />
<br />
Zi j<br />
2<br />
where Zi j represents the integrated impedance of the (heterogeneous) subsurface, and<br />
respectively, with i, j ∈ [x, y].<br />
φi j = arctan<br />
<br />
Im(Zi j)<br />
Re(Zi j)<br />
3.2.2. Relationships for simple subsurface cases<br />
1D case<br />
(3.35)<br />
(3.36)<br />
For situations where the subsurface exhibits solely vertical conductivity changes, naturally<br />
no alignment of any coordinate system axis with lateral interfaces can be made,<br />
making the choice of the coordinate system direction completely arbitrary. In such cases,<br />
diagonal elements of the impedance tensor (Eq. 3.34) are zero since electric fields are<br />
independent of parallel magnetic fields:<br />
Ex<br />
Ey<br />
<br />
=<br />
0 Zxy<br />
−Zxy 0<br />
Hx<br />
Hy<br />
<br />
. (3.37)<br />
The off-diagonal elements are of the same amplitude, but with an inverted sign due to the<br />
change from a right-handed into a left-handed coordinate system when dealing with the<br />
yx-component<br />
2D case<br />
When the subsurface possesses lateral conductivity interfaces in only one direction (with<br />
an arbitrary number of vertical conductivity changes), the original coordinate system used<br />
for the recording can be rotated to the interface directions using a transformation matrix<br />
<br />
cos(Θ) sin(Θ)<br />
R =<br />
. (3.38)<br />
− sin(Θ) cos(Θ)<br />
For the case of perfect adjustment the impedance matrix reduces to<br />
<br />
Ex 0<br />
=<br />
Zyx<br />
<br />
Zxy Hx<br />
.<br />
0<br />
(3.39)<br />
36<br />
Ey<br />
Hy