Modeling of seismo-electromagnetic phenomena
N. Gershenzon and G. Bambakidis

Appendix A

Crustal deformation is accompanied by the process of pore water diffusion. The hydrodynamics of pore water was described in the classic papers of Frenkel [1944] and Biot [1956]. We can write the following equation describing the relation between the change in volume strain q=DV/V and the hydrostatic pressure change P in a porous, water-saturated medium [Frenkel, 1944]:

eqn085.gif(A1)

where b=(1-K/K0)/m; bprime=1+(b-1)K2/K0; K0, K and K2 are the bulk moduli of rock matrix, dry porous rock, and pore water, respectively; m is the porosity; mv, k0 and rw are the dynamic viscosity, permeability, and density of water, respectively. Note that the magnitude of q depends in general on P. However, for most rock parameter values, the changes in pore pressure have practically no effect on the volume strain. This is why we shall suppose that q depends only on time t and position r, but is independent of P. So q (t, r) acts simply as a source function in equation A1.

Let's compare the first two terms on the left in equation A1, assuming a monochromatic disturbance of frequency w and time dependence eiwt. Then these terms are -w2 P/K2 and iwmv P/K2k0rw. If wggw constequivmv/k0rw then the second term on the left can be ignored. By a similar argument the second term on the right can be ignored relative to the first term. This means that the pore pressure will be described by the inhomogeneous wave equation,

eqn086.gif(A2)

From this equation we can define the speed of propagation of compressional waves in pore water:

eqn087.gif(A3)

Compare this speed with the speed of compressional waves in the matrix [Landau and Lifshitz, 1986]:

eqn088.gif(A4)

where ms is shear modulus and r matrix is the density of the matrix. Using typical parameters (Table A1), we find V2 water/V2 matrixapprox1/10.

From equation A2, when a seismic wave of frequency wggw const and wavelength l=2pVmatrix/w propagates through a porous medium, it will drive the pressure fluctuations in the water at the same frequency and wavelength. The second term on the left in equation A2 is therefore essentially w2 P/V2water while the term nabla2 P is of order (2p/l)2 P or w2 P/V2matrix. This latter term will be much smaller than the former; ignoring it gives

eqn089.gif(A5)

Integration results in

eqn090.gif(A6)

(Under conditions of hydrostatic equilibrium, the constants of integration vanish.)

Now consider the situation when wllw0. In this case we can neglect the terms in equation A1 containing the second time derivatives, leading to a diffusion-type equation:

eqn091.gif(A7)

It is easy to show that, in this case also, the pressure fluctuations will track the volume fluctuations. To see this we again consider a seismic wave as the driving force. The first term on the left is mv iwP/k0rw or iw constwP, while the second term is essentially (K0/bprimerw)w2 P/V2 matrix or (V2 water/V2 matrix)w2 P. The latter term can be ignored compared to the former, so

eqn092.gif(A8)

Therefore

eqn093.gif(A9)

So in order to connect P and q for fast processes when wggw const we should use equation A6 and for slow processes (like water diffusion after sudden changes of volume strain) we should use equation A9.


Appendix B

We present here asymptotic formulas for the components of the above-ground electromagnetic field produced by horizontal and vertical electric and magnetic dipoles embedded in the conducting half-space z>0 (see Figure 2a in Section 4) with conductivity s, electric permeability e and magnetic susceptibility m [Banos, 1966]. The distance r from the coordinate origin is characterized as being in the near zone, intermediate zone or far zone (see Section 4 and Figure 2b). In order to save space, we present results only for the dominant components in the intermediate and far zones. For the near zone, all components are given. The quantities n, k1 and k2 are defined in Section 4.

1. Horizontal electric dipole

(a) Near zone

eqn094.gif(B1)

eqn095.gif(B2)

eqn096.gif(B3)

eqn097.gif(B4)

eqn098.gif(B5)

eqn099.gif(B6)

(b) Intermediate zone

eqn100.gif(B7)

eqn101.gif(B8)

(c) Far zone

eqn102.gif(B9)

eqn103.gif(B10)

2. Horizontal magnetic dipole

(a) Near zone

eqn104.gif(B11)

eqn105.gif(B12)

eqn106.gif(B13)

eqn107.gif(B14)

eqn108.gif(B15)

eqn109.gif(B16)

(b) Intermediate zone

eqn110.gif(B17)

eqn111.gif(B18)

(c) Far zone

eqn112.gif(B19)

eqn113.gif(B20)

3. Vertical electric dipole

(a) Near zone

eqn114.gif(B21)

eqn115.gif(B22)

eqn116.gif(B23)

(b) Intermediate zone

eqn117.gif(B24)

eqn118.gif(B25)

(c) Far zone

eqn119.gif(B26)

eqn120.gif(B27)

4. Vertical magnetic dipole

(a) Near zone

eqn121.gif(B28)

eqn122.gif(B29)

eqn123.gif(B30)

(b) Intermediate and far zones

eqn124.gif(B31)

eqn125.gif(B32)


Appendix C

In Section 4, consideration of the electrokinetic effect resulted in equation 21b, in which the second and third terms involve surface integrals instead of volume integrals. To see how this comes about, consider the volume integral (see equation 12),

eqn126.gif

where j0 is the electrokinetic current and the integral is over the crust. If the crust is homogeneous this integral equals zero, since P is a localized function, i.e. vanishes at large distances. Now suppose there is an inhomogeneity across the plane x =0, resulting in a discontinuity in the coefficient C: C = C1H(x) - C2H(-x). Then

eqn127.gif

eqn128.gif

eqn129.gif

If P(r, t)approx P(0)lcr exp[- rL - ( t - r/VDt)2]H(r - lc), I can be written in the form

eqn130.gif

eqn131.gif

eqn132.gif

which is the second term in equation 21b. The same reasoning can be applied to the third term.


This document was generated by TeXWeb (Win32, v.1.3) on October 28, 2000.