Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet
Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet
4. Model calculations. for tl1.r,ee--d~nei~sioi~a1. structures - - 4.1. Introduction - In the three-dimensional case the TE- and TM-mode become mixed and cannot longer be treated separately. Now the differential equation for a vector field instead of a scalar field is to be solved. In numerical solutions questions of storage and computer time become important. Assume as example that in approach A a basic domain with 20 cells in each direction is chosen. In this case, only the storage of the electric field vector would require 48000 locations. For an iterative improvement of one field component at least 0.0005 sec are needed for each cell. This yields 12 sec for a complete itera- tion, and 20 min for 100 iterations. This appears to be the ].east time required for this model. Hence methods for a reduction of com- puter time and.storage are particularly appreciated in this case. The equation to be solved is cur1213 - (r) - -t k2 (r) E (r) = - i ~ l(g) i ~ ~ ~ --- where k2.(r) = i~p~c(g). - 1 (r) is the source current density. -e - After the splitting where E is that solution of -11 curlZE (r) + k2(r)F, (r) = -impo& -11 - n--n- (4.3) which vanishes at infinity, we obtain for the anomalous field the two alternative equations c ~ r l . + ~ k2E ~ = - k2E -a n-a a Eq. (4.4a) is the starting point for the volume integral or integral equation approach, Eq. (4.4b) is the point or-' starting for the sur- face integral approach.
4.2. In3ral - --- equation method Let -1 G.(r --a [r), -- i = 1,2,3 be a solution of cur12c. .-I (r --0 lr)+k2 - n (r)G. - --I (r --0 Ir) ,. A - = zi6(g-%), (4.5) vanishing at infinity. Here, the x. are unit vectors along the A h -l h A A A cartesian coordinate axes: x =x, x = -1 - -2 1, z3 = - 2. 14u'iul&31y (4.5) by E (r) and (4.4a) by G. -a - (r I?:) and integrate the difference with re- -1 -0 - spect to - r over the whole space. Green's vector theorem where d-r i.s a volume element, dA a surface element and - 1 the out- ward normal vector, yields since and G vanish at infinity. After combining all three compo-- -i nents and introducing - E instead of E , the vector integral equation -a 1 -gtg)d~ (4.7) -- is obtained. Here is Green's tensor being defined as wj(rJ % --O 3 3 A h A t(r Ir)= C x.G. (r [r)= Z G.. (r .Ir)x: x. - -1-1 -0 - ir j=1 13 -0 -- -1 -3 i= I (using dyadic notation). The tensor elements G.. admit a simple = 7 physical interpretation: Gij (r _o [r) - is the j-th electric field coin- ponent of an oscillating electric di-pole of unit moment pointlng in x.-direction, placed in the normal conductivity structure at &; 1 -- the point of observation is - r. Note that the first subscript and argument refer to the source, the second subscript and argument to the receiver. Because of the fundamental reci.p~:ocity in electro- mag~~etl.sm, source and observer parameters axe interchangeable, i.e. For a proof replace in (4.5) - r by - r', write an analogous equation for G.(rlrl), multiply cross-wise by G and Gi, integrate the -3 - - - j - difference wit11 respect to - r' over the whole space, and obtain (4.9) on using (4 .G) . Due to (4.9) , the equation (4.7) is alternatj.vely written
- Page 1 and 2: Electromagnetic Induction in the Ea
- Page 3 and 4: 6.2. Generalized matrix inversion 6
- Page 5 and 6: A1~Lernativel.y p = 30. m, where T
- Page 7 and 8: The e1ectrica:L effect - of the cha
- Page 9 and 10: The vari.ables x,y, and t which do
- Page 11 and 12: - d IufM12 > 0 . dz - - On the othe
- Page 13 and 14: a) TM-mode From (2.251, (2.26); (2.
- Page 15 and 16: with " the abbreviation + - 1 - ? =
- Page 17 and 18: 2n -i~cr cos (8-$1 -iKrcU J e ~B=J
- Page 19 and 20: In the 1irnl.t a -+ o, I +- m, M =
- Page 21 and 22: -- 2.5. Definition -- of the transf
- Page 23 and 24: we arrive at - v The same appl-ies
- Page 25 and 26: ) Computation of ---- C for a laxez
- Page 27 and 28: The approximate interpretation of C
- Page 29 and 30: I ~ispersi-on relations I Dispersio
- Page 31 and 32: where L is a positively oriented cl
- Page 33 and 34: 1 2 3 4 CPD 1 2 3 4 CPD - g) Depend
- Page 35 and 36: The TE-mode has no vertical electri
- Page 37 and 38: i I Earth Anomalous domain 3.2. Air
- Page 39 and 40: Hence, the conductivity is to be av
- Page 41 and 42: The RHS i.s a closed line integral
- Page 43 and 44: 4. Having determined B;, the coeffi
- Page 45 and 46: 3.4. Anomalous region as basic doma
- Page 47 and 48: - 6 and 6= can be so adjusted that
- Page 49 and 50: From the generalized Green's theore
- Page 51: and y can again be so adjusted that
- Page 55 and 56: The element GZx is needed for all z
- Page 57 and 58: With this knowledge of the behaviou
- Page 59 and 60: After having determined Qzr VJ,; @,
- Page 61 and 62: 4.3. The surface inteyral approach
- Page 63 and 64: F At the vertical boundaries the co
- Page 65 and 66: The four equations A A A A H = i sg
- Page 68 and 69: 6. Approaches to the inverse proble
- Page 70 and 71: to minimize the quantity a s = 12 /
- Page 72 and 73: It remains to show a way to minimiz
- Page 74 and 75: Agai-n, from a finite erroneous dat
- Page 76 and 77: Here lJ - is a N x P matrix contain
- Page 78 and 79: small eigenvalues. The parameter ve
- Page 80 and 81: Then - 77 - A(E2 - E ) = iwu U (E -
- Page 82 and 83: whence 2k d -2k d where a = CA:(A;)
- Page 84 and 85: . 7. Basic concepts of geomagnetic
- Page 86 and 87: orders of magnitude smaller' than t
- Page 88 and 89: Elimination of - E or .,. H yields
- Page 90 and 91: Observing that rot pot rot g = - ro
- Page 92 and 93: Two special types of such anomalies
- Page 94 and 95: Model : wo+ Solution for uniform ha
- Page 96 and 97: parameter u and that the pressure d
- Page 98 and 99: (=disturbed)-variations: After magn
- Page 100 and 101: with 4 as geographic latitude. From
4.2. In3ral - --- equation method<br />
Let -1 G.(r --a [r), -- i = 1,2,3 be a solution of<br />
cur12c. .-I (r --0 lr)+k2 - n (r)G. - --I (r --0 Ir)<br />
,.<br />
A<br />
- = zi6(g-%), (4.5)<br />
vanishing at infinity. Here, the x. are unit vectors along the<br />
A h -l h A A A<br />
cartesian coordinate axes: x =x, x =<br />
-1 - -2 1, z3 = - 2. 14u'iul&31y (4.5) by<br />
E (r) and (4.4a) by G.<br />
-a - (r I?:) and integrate the difference with re-<br />
-1 -0 -<br />
spect to - r over the whole space. Green's vector theorem<br />
where d-r i.s a volume element, dA a surface element and - 1 the out-<br />
ward normal vector, yields<br />
since and G vanish at infinity. After combining all three compo--<br />
-i<br />
nents and introducing - E instead of E , the vector integral equation<br />
-a<br />
1 -gtg)d~ (4.7)<br />
--<br />
is obtained. Here is Green's tensor being defined as<br />
wj(rJ<br />
% --O<br />
3 3 A h<br />
A<br />
t(r Ir)= C x.G. (r [r)= Z G.. (r .Ir)x: x.<br />
- -1-1 -0 - ir j=1<br />
13 -0 -- -1 -3<br />
i= I<br />
(using dyadic notation). The tensor elements G.. admit a simple<br />
= 7<br />
physical interpretation: Gij (r _o [r) - is the j-th electric field coin-<br />
ponent of an oscillating electric di-pole of unit moment pointlng<br />
in x.-direction, placed in the normal conductivity structure at &;<br />
1 --<br />
the point of observation is - r. Note that the first subscript and<br />
argument refer to the source, the second subscript and argument to<br />
the receiver. Because of the fundamental reci.p~:ocity in electro-<br />
mag~~etl.sm, source and observer parameters axe interchangeable, i.e.<br />
For a proof replace in (4.5) - r by - r', write an analogous equation<br />
for G.(rlrl), multiply cross-wise by G and Gi, integrate the<br />
-3 - - - j -<br />
difference wit11 respect to - r' over the whole space, and obtain (4.9)<br />
on using (4 .G) . Due to (4.9) , the equation (4.7) is alternatj.vely<br />
written