Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet

Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet

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

small eigenvalues. The parameter vector of a high eigenvalue shows the parameter combination which can be resolved well, the parameter eigenvector of a small eigenvalue gives the combination of para- meters for which only a poor resolution can be obtained. The generalized inverse is used both to invert a given data set and to estimate the information contents of this data set when the final solution is reached. hring the inversion procedure one has two too1.s to stabilize the notable unstable process of linearization: a) application of only a fraction of the con~puted parameter correc-- tion, leading to a trade-off between convergence rate and sta- bility; b) decrease of number of eigenvalues taken into account. In the final estimation of the data contents one might prescribe for exh parameter a maximum variance. Then one has to determine from (6.35) the number of eigenvalues leading to a value nearest to the prescribed one. Finally the row of the resolution matrix for this particular parameter is calculated. 6.3. 'mri.vation of the kernels for the linearized inverse problem of electromagnetic induction 6.3.1. The one-dimensional case Both for the Backus-Gilbert procedure and for the generalized linear inversion a knowledge of the change of the data due to a small change in the conductivity structure is required. In the one-dimensional case the pertinent differential equation and datum are , flr(z,w) = I K ~ + iwp0o(z)}f (zrw) Consider two conductivity profiles o and o with correspondicg 1 2 fields fl and f2. Multiplying the equation for fl with f2 and the equation for f2 with flr subtracting the resnltinq equations and - integrating the difference over z from zero to infinity, we obtain .I (f;'f2 - f:fl)dz = iwp o 1 (ol-cr2)flf2 dZ 0 0

Now m Division by f !, (0) :f; (0) yields If the difference 60 = a - o is small, f2 in the integral may 2 I be replaced by fl, since the difference f2 - fl is of the order of 60 = a2 - o leading to a second order term in 6C = C2 1' - C1. Hence to a first order in 6s m 6C(w) = - iwp l6a(z){ f (z,w) 12 dz O 0 f' (0,~) Therefore in the Backus-Gilbert procedure the FrGchet derivative of C is -iwpb{f(~)/f'(o)}~. In the generalized inversion the deri- vatives of C with respect to layer conductivities and layer thick- nesses is required, if a structure with uniform layers is assumed. Let there be L layers with conductivity am and thickness d in m the m-th layer, h < z 2 hm+l (hL+l = -). m - Then (6.38) yields since in the last case all layers below the m-th layer are dis- placed too.' ' '6':3 :2 : Pa,rt?al - derivatives in two- and 'three-dimensions For two-dimensions only the E-polarization case is considered. The pertinent equation is (cf . (3.5a) ) AE = iwp a E . 0 Considef again two conductivity structures a and a2. 1

Now<br />

m<br />

Division by f !, (0) :f; (0) yields<br />

If the difference 60 = a - o is small, f2 in the integral may<br />

2 I<br />

be replaced by fl, since the difference f2 - fl is of the order<br />

of 60 = a2 - o leading to a second order term in 6C = C2<br />

1'<br />

- C1.<br />

Hence to a first order in 6s<br />

m<br />

6C(w) = - iwp l6a(z){ f (z,w) 12 dz<br />

O 0 f' (0,~)<br />

Therefore in the Backus-Gilbert procedure the FrGchet derivative<br />

of C is -iwpb{f(~)/f'(o)}~. In the generalized inversion the deri-<br />

vatives of C with respect to layer conductivities and layer thick-<br />

nesses is required, if a structure with uniform layers is assumed.<br />

Let there be L layers with conductivity am and thickness d in<br />

m<br />

the m-th layer, h < z 2 hm+l (hL+l = -).<br />

m - Then (6.38) yields<br />

since in the last case all layers below the m-th layer are dis-<br />

placed too.'<br />

' '6':3 :2 : Pa,rt?al - derivatives in two- and 'three-dimensions<br />

For two-dimensions only the E-polarization case is considered. The<br />

pertinent equation is (cf . (3.5a) )<br />

AE = iwp a E .<br />

0<br />

Considef again two conductivity structures a and a2.<br />

1

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

Saved successfully!

Ooh no, something went wrong!