P. Schmoldt, PhD - MTNet - DIAS

P. Schmoldt, PhD - MTNet - DIAS P. Schmoldt, PhD - MTNet - DIAS

mtnet.dias.ie
from mtnet.dias.ie More from this publisher
04.08.2013 Views

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

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!