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

It remains to show a way to minimize Q, subject to (6.15). The way, however, is well-known. One simply intrsduces as (N+1 )-st unknown a Lagrangian parameter X and minimizes the quantity a Q' = Q + h (Cai~li - 1 ) , ui = / Gi (r)dr (6.16) 0 The differentiation of (6.16) with respect to the (N4-1) unknowns yields the (N+1) linear equations where Qik = W Sik + (l-W) c Eik. This system of equations is easily solved. The meaning of X is revealed by multiplying (6.17a) by ak, adding and using (6.17b) and (6.15) . It results In the case W = 1, we have X = -2s. With a knowledge of a (roll i is obtained from (6.8) with gi = y. and E' from (6.9). 3- When is inserted in (6.1) instead of m(r) , it w ill in general not exactly reproduce the data. The minimization (6.15) N subject to the constraint C aiui = 1 admits for two data (N:=2) i=1 a simple geometrical interpretation: For constant s and E' these positive definite quantities are represented by ellipses, the con- straint is a line in the (al,a2)-plane. For uncorrelated errors, the principal axes of E~ are the al and a2 axis. hT1a W varies from 0 to 1 the combinatiomof (al ,a2) on the fat line are obtained. s and. E' are determined from the ellipses throughthese points,Sj-nce all s-(~~)-elli~ses are similar, s and^ 2

are proportional to the long axes of these ellipses. 6.1.3. The nonlinear inverse problem The Backus-Gilbert procedure applies only to linear inverse problems, where according to a the gross earth functionals gi have the property that This means for instance that the data are bui1.t up in an aclditive way from different parts of the model, i.e., that there is no coupling between these parts. This certainly does not hold for electromagnetic inverse problem, where each part of the conductor is coupled with all other parts. In nonlinear problems the data kernel Gi(r) will depend on m. Here it is in general possible to replace (6.18) by . a gi.(mr)-gi~m)=f(m'(r)-m(r))~i(x,m)dr + ~(m'-m)~, (6.18a) 0 where m and m' are two earth models. The data kernel Gi(r,m) is called the Fr6chet derivative or functional derivative at model m.

It remains to show a way to minimize Q, subject to (6.15). The way,<br />

however, is well-known. One simply intrsduces as (N+1 )-st unknown<br />

a Lagrangian parameter X and minimizes the quantity<br />

a<br />

Q' = Q + h (Cai~li - 1 ) , ui = / Gi (r)dr (6.16)<br />

0<br />

The differentiation of (6.16) with respect to the (N4-1) unknowns<br />

yields the (N+1) linear equations<br />

where Qik = W Sik + (l-W) c Eik.<br />

This system of equations is easily solved. The meaning of X is<br />

revealed by multiplying (6.17a) by ak, adding and using (6.17b)<br />

and (6.15) . It results<br />

In the case W = 1, we have X = -2s. With a knowledge of a (roll<br />

i<br />

is obtained from (6.8) with gi = y. and E' from (6.9).<br />

3-<br />

When is inserted in (6.1) instead of m(r) , it w ill in<br />

general not exactly reproduce the data.<br />

The minimization (6.15)<br />

N<br />

subject to the constraint C aiui = 1 admits for two data (N:=2)<br />

i=1<br />

a simple geometrical interpretation: For constant s and E' these<br />

positive definite quantities are represented by ellipses, the con-<br />

straint is a line in the (al,a2)-plane. For uncorrelated errors,<br />

the principal axes of E~ are the al and a2 axis.<br />

hT1a W varies from 0 to 1 the combinatiomof (al ,a2) on the fat<br />

line are obtained. s and. E' are determined from the ellipses<br />

throughthese points,Sj-nce all s-(~~)-elli~ses are similar, s and^<br />

2

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

Saved successfully!

Ooh no, something went wrong!