04.08.2013 Views

P. Schmoldt, PhD - MTNet - DIAS

P. Schmoldt, PhD - MTNet - DIAS

P. Schmoldt, PhD - MTNet - DIAS

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

6. Using magnetotellurics to gain information about the Earth<br />

mean cross-spectral density (S xy), which in turn is the product of the respective FT’s X( f )<br />

and Y( f ), i.e.<br />

C 2 xy( f ) =<br />

<br />

<br />

〈S xy( f )〉 2<br />

〈S xx( f )〉〈S yy( f )〉 =<br />

|〈X ∗ ( f )Y( f )〉| 2<br />

〈X ∗ ( f )X( f )〉〈Y ∗ ( f )Y( f )〉<br />

(6.3)<br />

where the asterisk denotes the complex conjugate of the function. Coherence weighting<br />

is implicitly applied during the impedance estimation using the remote reference method,<br />

in which modern processing schemes also consider partial and multiple coherence as indicators<br />

for the quality of an impedance estimate.<br />

To illustrate the application of this procedure to an MT dataset consider the explicit<br />

form of Equation 3.34 for the electric field component in x-direction (Ex), i.e.<br />

Ex = ZxxHx + ZxyHy<br />

(6.4)<br />

where both, the electric and magnetic channels Ex and Hy are assumed here to contain<br />

noise. Multiplication of this Equation with the complex conjugate of either the electric or<br />

magnetic field component in x-direction (R ∗ x) or in y-direction (R ∗ y) yields<br />

and<br />

〈ExR ∗ x〉 = Zxx〈HxR ∗ x〉 + Zxy〈HyR ∗ x〉 (6.5)<br />

〈ExR ∗ y〉 = Zxx〈HxR ∗ y〉 + Zxy〈HyR ∗ y〉, (6.6)<br />

respectively. Combining these two equations results in an estimate of the Zxy component<br />

Zxy = 〈ExR∗ y〉〈HxR∗ x〉 − 〈ExR∗ x〉〈HxR∗ y〉<br />

〈HxR∗ x〉〈HyR∗ y〉 − 〈HxR∗ y〉〈HyR∗ . (6.7)<br />

x〉<br />

Expressions for all components of the impedance tensor can be derived in a similar manner,<br />

therefore the general form can be stated as<br />

Zi j = 〈EiR∗ ∗<br />

∗<br />

j 〉〈HkRk 〉 − 〈EiRk 〉〈HkR∗ j 〉<br />

, (6.8)<br />

DET<br />

with DET = 〈HxR ∗ x〉〈HyR ∗ y〉 − 〈HxR ∗ y〉〈HyR ∗ x〉, and i, j, k ∈ [x, y], with k j. Whether processing<br />

with the electric or the magnetic component as reference provides superior results<br />

depends on the type of noise and its effect on each channel. The magnetic component is<br />

commonly assumed to be less contaminated by noise and usually chosen as remote, however<br />

the choice is often a rather subjective one. To use a combination of remote references<br />

from different components or stations can often be useful.<br />

112<br />

For the case of a perfectly 2D subsurface and according rotation of the magnetic field

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

Saved successfully!

Ooh no, something went wrong!