Radial waves in a borehole and stoneley waves for measuring formation permeability and electroacoustic constant

Abstract

A Stoneley wave is generated in a borehole in a saturated porous earth formation. Measurements are made of the velocity of motion of the formation and the fluid in the formation. The difference in the velocities is indicative of formation permeability. An additional measurement of the electric field at the borehole wall enables determination of the electroacoustic constant.

Claims

1. A method of evaluating an earth formation intersected by a borehole, the method comprising: using an acoustic transmitter to generate acoustic signals at a frequency producing a resonance of fluid wherein said resonance generates a velocity of the fluid in the borehole and a velocity of the formation at an interface between the fluid and the formation such that a difference between the velocity of the fluid in the borehole and the velocity of the formation is dependent upon permeability of the formation; and making a first measurement of the fluid velocity; making a second measurement of the velocity of the formation; making a third measurement indicative of an electric charge on a wall of the borehole taken at the frequency; and using a processor to estimate an electroacoustic constant of the earth formation using the first measurement, the second measurement, and the third measurement.

2. The method of claim 1 wherein the acoustic transmitter comprises a swept frequency source, and wherein using the acoustic transmitter further comprises using the swept frequency source in the borehole and identifying the frequency by identifying from a plurality of frequencies in a frequency sweep the frequency at which a corresponding difference between a corresponding fluid velocity in the borehole and a corresponding velocity of the formation is at a maximum.

3. The method of claim 2 further comprising using a geophone for measuring the velocity of the formation and using a flow rate sensor to measure the fluid velocity.

4. The method of claim 2 further comprising: making an additional first measurement of the fluid velocity and an additional second measurement of the velocity of the formation at at least one additional frequency different from the frequency; and using the additional first measurement and the additional second measurement to estimate a permeability of the earth formation.

5. The method of claim 4 further comprising using a measured value of a first compressional wave in the formation, a second compressional waves in the formation, and a shear velocity of the formation to estimate the permeability.

6. The method of claim 2 wherein using the swept frequency source further comprises using a monopole source.

7. The method of claim 2 further comprising conveying the swept frequency source into the borehole on a conveyance device selected from: (i) a wireline, and (ii) a drilling tubular.

8. The method of claim 1 wherein making the third measurement indicative of the electric charge further comprises at least one of (i) using a grounded ammeter, and (ii) making a measurement in a pilot hole in the borehole wall.

9. An apparatus configured to evaluate an earth formation intersected by a borehole, the apparatus comprising: a device including an acoustic transmitter configured to generate acoustic signals at a frequency producing a resonance of fluid wherein said resonance generates a velocity of the fluid in the borehole and a velocity of the formation at an interface between the fluid and the formation such that a difference between the velocity of the fluid in the borehole and the velocity of the formation is dependent upon permeability of the formation; a first sensor configured to make a first measurement of the fluid velocity; a second sensor configured to make a second measurement of the velocity of the formation; a third sensor configured to make a third measurement indicative of an electric charge on a wall of the borehole taken at the frequency; and a processor configured to estimate an electroacoustic constant of the earth formation using the first measurement, the second measurement, and the third measurement.

10. The apparatus of claim 9 wherein the acoustic transmitter comprises a swept frequency source; and wherein the processor is further configured to identify the frequency by identifying from a plurality of frequencies in a frequency sweep the frequency at which a corresponding difference between a corresponding fluid velocity in the borehole and a corresponding velocity of the formation is at a maximum.

11. The apparatus of claim 10 further comprising a geophone configured to measure the velocity of the formation and a flow rate sensor configured to measure the fluid velocity.

12. The apparatus of claim 10 wherein using the swept frequency source further comprises using a monopole source.

13. The apparatus of claim 10 further comprising a conveyance device configured to convey the swept frequency source into the borehole, the conveyance device being selected from: (i) a wireline, and (ii) a drilling tubular.

14. The apparatus of claim 9 wherein the third sensor is selected from: (i) a grounded ammeter, and (ii) a sensor in a pilot hole in the borehole wall.

15. The apparatus of claim 9 wherein the processor is further configured to estimate a permeability of the earth formation using an additional first measurement of the fluid velocity and an additional second measurement of the velocity of the formation at at least one additional frequency of the fluid different from the frequency.

16. The apparatus of claim 15 wherein the processor is further configured to use a measured value of a first compressional wave in the formation, a second compressional wave in the formation, and a shear velocity of the formation to estimate the permeability.

Description

BRIEF DESCRIPTION OF THE FIGURES

(1) For detailed understanding of the present disclosure, references should be made to the following detailed description of exemplary embodiment(s), taken in conjunction with the accompanying drawings, in which like elements have been given like numerals, wherein:

(2) FIG. 1 shows an acoustic well logging instrument as it is used to acquire signals useful with the method of this invention;

(3) FIG. 2 is an illustration of an exemplary configuration of acoustic transmitters, a flow sensor, a geophone and an electrical charge sensor;

(4) FIG. 3 is a plot of the difference between the velocities of the saturating fluid and the porous matrix a function of frequency;

(5) FIG. 4 shows an exemplary dependence of the difference in velocities of the saturating fluid and the porous matrix as a function of permeability; and

(6) FIG. 5 shows the dependence of the ratio of the electric field to the velocity of deformation of the borehole upon permeability for different frequencies.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(7) FIG. 1 shows an acoustic well logging instrument as it is used to acquire signals suitable for processing according to the method of this disclosure. The instrument 10 is inserted into a wellbore 2 drilled through earth formations 20. The instrument 10 can be inserted into and withdrawn from the wellbore 2 by means of an armored electrical cable 14 spooled from a winch (not shown) or any similar conveyance known in the art.

(8) The wellbore 2 is typically filled with a liquid 4, which can be drilling mud or any similar fluid usually used for drilling or completion of wellbores. The instrument includes a plurality of acoustic transmitter 8. The transmitters 8 are periodically energized and emit acoustic energy that radiate from the tool 10. The instrument 10 typically includes a telemetry module, shown generally at 12, which converts the electrical signals from the sensors (discussed with reference to FIG. 2) into a form suitable for recording and transmission over the cable 14 to a surface processor 18. It should be understood that the number of transmitters 8 as shown in FIG. 1 is only an example and is not meant to limit the invention.

(9) Shown in FIG. 2 are an array of transmitters 241a . . . 241n. The arrays may include elements that are axially and/or circumferentially disposed. While the illustration shows them on a single housing, this is not to be construed as a limitation to the disclosure; a common configuration uses transmitters and receivers on more than one sub. Also shown in the figure is the borehole 26, and the logging tool 100. Also shown are a sensor 271 to measure the electric charge or field at the borehole wall, a sensor 273, such as a flow rate sensor, to measure the borehole fluid velocity and a motion sensor such as a geophone 275 to measure the velocity of the deformation of the borehole wall. The present disclosure includes two methods of measuring the surface charge (or the electric field in the porous medium). One is to connect an ammeter to the borehole wall as sensor 271 with the other end of the ammeter grounded: the measured current is indicative of the surface charge. A second method is drill pilot holes 226 in the formation and measure the potential difference between the borehole wall and the bottom of the pilot hole using a sensor 277 in the bottom of the pilot hole 226.

(10) It should be noted that the disclosure is with reference to a wireline conveyed logging tool. This is not to be construed as a limitation, and the transmitter array for generating Stoneley waves could be on a bottomhole assembly (BHA) and used for imaging ahead of the borehole. This is discussed in U.S. Pat. No. 8,055,448 of Mathiszik et al., having the same assignee as the present disclosure and the contents of which are incorporated herein by reference.

(11) With a monopole excitation of the elements of the transmitter array 241a . . . 241n, radial excitation of the fluid in the borehole results. The theory of the radial oscillations, as it is shown below, shows that the electric charge appearing on the surface of the interface is determined by the following ratio:

(12) $\begin{array}{cc}Q{|}_{\mathrm{boundary}}=\frac{E_r{|}_{\mathrm{boundary}}}{4}=.Math.\frac{{}_l\left(v_r-u_r\right){|}_{\mathrm{boundary}}}{4}.& \left(1\right)\end{array}$
In the above formula, is the electroacoustic constant of the porous medium (the earth formation), is electric conductivity, .sub.l is partial density of fluid in the saturated porous medium (earth formation), v.sub.ru.sub.r is the difference between velocities of the saturating fluid and the porous matrix in the radial oscillations mode outside of the resonance zones. Having measured the difference between velocities v.sub.ru.sub.r (to this end, we could perform independent measurements of the fluid velocity v.sub.r.sup.(0) and velocity u.sub.r.sup.(1) of porous matrix next to the interface between the media, and based on these data, calculate the difference between velocities u.sub.r.sup.(1)v.sub.r.sup.(1)=(u.sub.r.sup.(1)v.sub.r.sup.(0))/) and the surface electric charge using a special sensor, we can calculate the electroacoustic constant via the formula (1) above:

(13) $\begin{array}{cc}=\frac{4Q{|}_{\mathrm{boundary}}}{{}_l\left(v_r-u_r\right){|}_{\mathrm{boundary}}}.& \left(2\right)\end{array}$
In order to estimate the electroacoustic constant using this method, the method includes a measurement of formation porosity using, for example, a neutron porosity logging tool. The partial density of the liquid .sub.l is given by the product .sub.l.sup.ph where the physical fluid density .sub.l.sup.ph is measured by a density logging tool, such as a gamma ray density tool. The conductivity of the formation fluid may be measured by known methods.

(14) An important comment should be made here. In the first approximation, the difference v.sub.ru.sub.r does not depend on permeability of the porous space k. FIG. 3 shows the dependence of the difference between velocities v.sub.ru.sub.r 301 (ordinate) upon the frequency of the signal (abscissa) for an exemplary model. The model parameters are: borehole liquid density .sub.f=1.0 g/cm.sup.3; sound velocity in borehole liquid c.sub.p0=1.45.Math.10.sup.5 cm/s; porous matrix density .sub.s=2.2 g/cm.sup.3; saturated liquid density .sub.l=1.0 g/cm.sup.3; porosity =0.2; first longitudinal sound velocity c.sub.p1=2.2.Math.10.sup.5 cm/s; second longitudinal sound velocity c.sub.p2=1.2.Math.10.sup.5 cm/s; share sound velocity c.sub.t=1.6.Math.10.sup.5 cm/s; liquid viscosity =1.05.Math.10.sup.2 P; permeability k=50 mD; source radius r.sub.1=1.0 cm; borehole radius r.sub.2=10 cm.) The intersection of 301 with the line 303 gives an indication of the width of the resonance zones. The derivation of the curves follows from the theory presented below, discussed in Appendix I.

(15) For the case where the formation is impermeable, then the configuration shown in FIG. 2 may be recognized as a waveguide for which the resonance frequencies depend upon the radius of the tool, the radius of the borehole and the velocity of compressional waves in the borehole. When the formation has a non-zero permeability, the configuration in FIG. 2 acts as a leaky waveguide with the leakage characteristics being dependent upon the permeability. This is the basis for the determination of formation permeability in the present disclosure.

(16) FIG. 4 shows the dependence of the resonance on permeability. The abscissa is the permeability. The curve 401 shows the dependence of the difference in velocities on the permeability in a narrow resonance zone about 300 Hz wide at a frequency of 27.39 kHz. The curve 403 shows the dependence of the difference in velocities on permeability outside the resonance zone. Thus, outside of the resonance zone, the difference between velocities depends weakly on permeability, and we may assume, in the first approximation, that the surface electric charge formed on the interface does not depend on permeability as the surface electric charge is a function primarily of the electroacoustic constant. Similar results can be obtained for the other resonances in FIG. 3. We next discuss the electroacoustic ratio of the Stoneley wave.

(17) As derived in the theory presented below in Appendix II, for the Stoneley wave, the ratio of the electric field E.sub.z amplitude to the amplitude V.sub.z in the acoustic wave which determined this field has the following form:

(18) $\begin{array}{cc}\frac{E_z}{V_z}=.Math.k.Math.f\left(,k\right),& \left(3\right)\end{array}$
where (, k) is a function weakly dependent on its arguments. This dependence is illustrated in FIG. 5. Shown therein is a plot of

(19) $\frac{E_z}{V_z}$
as a function of permeability for frequencies of 30 kHz 501, 15 kHz 503 and 3 kHz 505.

(20) When permeability changes by three orders of magnitude, the function changes only by the factor of two. Having the theoretical value of f(, k), measured value of E.sub.z/V.sub.z, and computed value of (see eqn. (2)), we can find permeability:

(21) $\begin{array}{cc}k=\frac{E_z/V_z}{f\left(,k\right)}.& \left(4\right)\end{array}$
Thus, having performed two types of measurements: one with the radial waves and one with the electroacoustic ratio of the Stoneley wave, we first measure the electroacoustic constant and then, permeability. In order to improve accuracy of subsequent measurements, in one embodiment of the disclosure, an iterative procedure is used taking into account the weak functional dependence (,k) and (v.sub.ru.sub.r)|.sub.boundary(, k). In such an iterative process, a first estimate of permeability is made ignoring the functional dependence, and the first estimate is used to calculate the functional dependence of f(, k) and a revised value of the permeability is calculated.

(22) The description above has been in terms of a device conveyed on a BHA on a drilling tubular into a borehole in the earth formation. The method and apparatus described above could also be used in conjunction with a logging string conveyed on a wireline into the earth formation. For the purposes of the present disclosure, the BHA and the logging string may be referred to as a downhole assembly. It should further be noted that while the example shown depicted the transmitter assembly and the receiver assembly on a single tubular, this is not to be construed as a limitation of the disclosure. It is also possible to have a segmented acoustic logging tool to facilitate conveyance in the borehole. Once the formation permeability has been estimated, it can be used for further reservoir development operations.

(23) Implicit in the processing of the data is the use of a computer program implemented on a suitable machine readable medium that enables the processor to perform the control and processing. The machine readable medium may include ROMs, EPROMs, EAROMs, Flash Memories and Optical disks. The determined formation velocities and boundary locations may be recorded on a suitable medium and used for subsequent processing upon retrieval of the BHA. The determined formation permeability and electroacoustic constant may further be telemetered uphole for display and analysis.

APPENDIX I

(24) To measure permeability in the saturated porous medium, a theory that describes the radial waves based on the linearized version of the filtration theory is used. To describe the acoustic field in fluid, the following equation is used:
{umlaut over (v)}c.sub.p0.sup.2 div v=0.(I.1)
Here c.sub.p0 is the velocity of sound in borehole fluid. The electrical field in porous medium is defined by

(25) $\begin{array}{cc}E=-\frac{{}_l}{}\left(u-v\right),& \left(I\mathrm{.2}\right)\end{array}$
where is the electroacoustic constant of the porous medium, and is the electric conductivity of the porous medium. Velocities u, v of the matrix of the porous medium and the fluid contained therein satisfy equations

(26) $\begin{array}{cc}\stackrel{\mathrm{.Math.}}{u}-c_t^{2}u-a_1\mathrm{div}u+a_2\mathrm{div}v+\frac{{}_l}{{}_s}b\left(\stackrel{.}{u}-\stackrel{.}{v}\right)=0,& \left(I\mathrm{.3}\right)\\ \stackrel{\mathrm{.Math.}}{v}+a_3\mathrm{div}u-a_4\mathrm{div}v-b\left(\stackrel{.}{u}-\stackrel{.}{v}\right)=0.& \left(I\mathrm{.4}\right)\end{array}$
Here .sub.l, .sub.s are partial densities of the saturating fluid and the porous matrix, respectively, =.sub.l+.sub.s is density of the saturated porous medium. The dissipative coefficient is b=.sub.l=/(k), where is dynamic viscosity the saturating fluid, k is permeability.

(27) Equation coefficients .sub.i (j=1, . . . , 4) are determined by the elasticity moduli , , and :

(28) $\begin{array}{cc}{a}_1=\frac{1}{{}_s}\left(\frac{{}_s^2}{{}^2}+\frac{{}_l^2}{{}^2}K+\frac{1}{3}\right),a_2=\frac{{}_l}{{}_s}\left(\frac{{}_l}{{}^2}K-\frac{{}_s}{{}^2}\right),& \left(I\mathrm{.5}\right)\\ a_3=\frac{{}_l}{{}^2}K-\frac{{}_s}{{}^2},a_4=\frac{{}_l}{{}^2}K+\frac{{}_l}{{}^2},& \left(I\mathrm{.6}\right)\end{array}$
where K=+2/3. Three elasticity moduli , , and are determined univalently by three measurable velocities of sound in saturated porous medium (two compressional velocities c.sub.p1, c.sub.p2 and one shear velocity c.sub.t)

(29) $\begin{array}{cc}={}_s{c}_t^2,& \left(I\mathrm{.7}\right)\\ K=\frac{1}{2}\frac{{}_s}{{}_l}\left(c_{p1}^2+c_{p2}^2-\frac{8}{3}{}_l{c}_t^2-\sqrt{{\left(c_{p1}^2-c_{p2}^2\right)}^2-\frac{64}{9}{}_s{}_l{c}_t^4}\right),& \left(I\mathrm{.8}\right)\\ =\frac{1}{2}\left(c_{p1}^2+c_{p2}^2-\frac{8}{3}{}_s{c}_t^2+\sqrt{{\left(c_{p1}^2-c_{p2}^2\right)}^2-\frac{64}{9}{}_s{}_l{c}_t^4}\right).& \left(I\mathrm{.9}\right)\end{array}$

(30) The acoustic source in the system is a finite-size harmonic pressure generator p=P.sub.0e.sup.it located in the center of the wellbore. The acoustic response of the system under study is also harmonic (u, v).fwdarw.(u, v)e.sup.it. The latter fact enables us to keep the following coordinate dependence of equations (I.1, I.2, I.3) in Fourier amplitudes of the radial wave velocities u(u.sub.r, 0, 0), v(v.sub.r, 0, 0):
.sub.0 div v+v=0(I.10)
And
.sub.0u.sub.2v+(1+b)ubv=0, (r.sub.2<r<),(I.11)
.sub.3u+.sub.4vbu+(1+b)v=0, (r.sub.2<r<).(I.12)
Here we have dimensionless equations where velocities u, v are in units of c.sub.t, coordinates are in units of c.sub.t/; coefficients are transformed: .sub.0=c.sub.p0.sup.2/c.sub.t.sup.2, .sub.i=a.sub.i/c.sub.t.sup.2 (i=0, . . . , 4), b=ib/; =.sub.l/.sub.s, .sub.0=.sub.1+1.

(31) The electrical field in dimensionless form is defined by

(32) 0$\begin{array}{cc}E=-\frac{c_e}{c_t}\left(u-v\right),& \left(I\mathrm{.13}\right)\end{array}$
where c.sub.e is the velocity of light, E is in units of c.sub.e.sub.l/. Let us find common solutions to equations (I.10-I.12).

(33) To reduce the equations described the acoustic field in porous medium to the diagonal form, the set (I.11-I.12) could be conveniently expressed in the matrix form.
AU+WU=0(I.14)
Here

(34) $\begin{array}{cc}U=\left(\begin{array}{c}{u}_r\\ v_r\end{array}\right),A=\left(\begin{array}{cc}{\stackrel{\_}{a}}_0& -{\stackrel{\_}{a}}_2\\ -{\stackrel{\_}{a}}_3& {\stackrel{\_}{a}}_4\end{array}\right),W=\left(\begin{array}{cc}\left(1+\mathrm{.Math.}\stackrel{\_}{b}\right)& -\mathrm{.Math.}\stackrel{\_}{b}\\ -\stackrel{\_}{b}& \left(1+\stackrel{\_}{b}\right)\end{array}\right).& \left(I\mathrm{.15}\right)\end{array}$
Having multiplied set (I.12) by the inverse matrix

(35) $W^{-1}=\left(\begin{array}{cc}\frac{\left(1+\stackrel{\_}{b}\right)}{\det W}& \frac{\mathrm{.Math.}\stackrel{\_}{b}}{\det W}\\ \frac{\stackrel{\_}{b}}{\det W}& \frac{\left(1+\mathrm{.Math.}\stackrel{\_}{b}\right)}{\det W}\end{array}\right),$
where det W=1+(1+)b, we arrive at the set of equations:
HU+U=0(I.16)
with this matrix of coefficients: H=W.sup.1A
h.sub.11=((1+b).sub.0b.sub.3)/det W, h.sub.12=((1+b).sub.2+b.sub.4)/det W,
h.sub.21=(b.sub.0(1+b).sub.3)/det W, h.sub.22=(b.sub.2+(1+b).sub.4)/det W.(I.17)

(36) To reduce set (I.14) to the diagonal form, we multiply it by the matrix

(37) $R=\left(\begin{array}{cc}{R}_1& R_2\\ R_3& R_4\end{array}\right),$
as found from this problem solution for the eigenvalues of the matrix
H:HR=R(I.18)
as the matrix compiled of the right-hand side eigenvectors,

(38) $R:\left(\left(\begin{array}{c}{R}_1\\ R_3\end{array}\right)\left(\begin{array}{c}{R}_2\\ R_4\end{array}\right)\right).$
Components of the transform matrix R (or eigenvectors of the matrix H) have the following:
R.sub.1=1, R.sub.2=(.sub.2h.sub.22)/h.sub.21, R.sub.3=h.sub.21/(.sub.1h.sub.22), R.sub.4=1.(I.19)

(39) The matrix of eigenvectors

(40) $=\left(\begin{array}{cc}{}_1& 0\\ 0& {}_2\end{array}\right)$
is diagonal, so that R.sup.1H R=, and eigenvalues .sub.1, .sub.2 can be expressed as follows:
.sub.1=tr H+{square root over ((tr H).sup.2det H)}, .sub.2=tr H{square root over ((tr H).sup.2det H)},(I.20)
where tr H=h.sub.11+h.sub.22 is the trace of the matrix H, det H=h.sub.11h.sub.22h.sub.12h.sub.21 is the determinant of matrix H.
The vector-column of functions is transformed as follows:

(41) $\begin{array}{cc}{R}^{-1}u=V,\mathrm{where}V=\left(\begin{array}{c}{v}_1\\ v_2\end{array}\right).& \left(I\mathrm{.21}\right)\end{array}$
After transformation, set (I.16) becomes diagonal:
V+V=0,(I.22)
or, in coordinate notations (prime denotes the derivative with respect to the radius):

(42) $\begin{array}{cc}{v}_1^{}+\frac{1}{r}{v}_1^{}+\left(\frac{1}{{}_1}-\frac{1}{r^2}\right)v_1=0& \left(I\mathrm{.23}\right)\\ v_2^{}+\frac{1}{r}{v}_2^{}+\left(\frac{1}{{}_2}-\frac{1}{r^2}\right)v_2=0& \left(I\mathrm{.24}\right)\end{array}$

(43) The initial velocities can be expressed via functions v.sub.1, v.sub.2 (I.21):
u.sub.r=v.sub.1+R.sub.2v.sub.2, v.sub.r=R.sub.3v.sub.1+v.sub.2.(I.25)

(44) The Acoustic equation (I.10) for the borehole fluid can be expressed similarly:

(45) $\begin{array}{cc}{v}_r^{}+\frac{1}{r}{v}_r^{}+\left(\frac{1}{{}_0}-\frac{1}{r^2}\right)v_r=0.& \left(I\mathrm{.26}\right)\end{array}$
The set of equations (I.23, I.24, I.26) taking into account (I.25) describes the acoustic field of the radial waves in the wellbore and in the formation beyond the wellbore. The boundary conditions for this geometry (in the dimensionless case) are the following (the index in brackets denotes values pertaining to the fluid .sup.(0) and the saturated porous medium .sup.(1)): 1) On the surface of the source exciting the oscillations in the fluid, (r=r.sub.1, FIG. 2), a harmonic source is set:
p.sup.(0)=P.sub.0(I.27) 2) On the interface between the borehole and the saturated porous medium r=r.sub.2, FIG. 2), continuity of the medium is assumed, as well as continuity of the full stress tensors (.sub.rr.sup.(1), .sub.rr.sup.(0)) and the partial pressures:
(1)u.sub.r.sup.(1)+v.sub.r.sup.(1)=v.sub.r.sup.(0),(I.28)

(46) $\begin{array}{cc}\underset{\mathrm{rr}}{\overset{\left(1\right)}{.Math.}}=\stackrel{\left(0\right)}{\underset{\mathrm{rr}}{.Math.}},& \left(I\mathrm{.29}\right)\\ \frac{p^{\left(1\right)}}{{}^{\left(1\right)}}-\frac{{}^{\left(1\right)}}{{}_s^{\left(1\right)}}\stackrel{\_}{b}\left({}_u^{\left(1\right)}-{}_v^{\left(1\right)}\right)=\frac{p^{\left(0\right)}}{{}^{\left(0\right)}}.& \left(I\mathrm{.30}\right)\end{array}$
Here is porosity, .sub.u,v are potentials which determine the velocity of the radial waves in the medium u=.sub.u, v=.sub.v, is a conditional coefficient which enables us to change the friction force at the boundary, P.sub.0=P.sub.0/(.sub.sc.sub.t.sup.2).

(47) The final boundary condition follows from a more general continuity condition for the partial stress tensor in fluid.
.sub.rr.sup.l(1)=.sub.rr.sup.(0),(I.30)
where .sub.rr.sup.(1)=.sub.rr.sup.s(1)+.sub.rr.sup.l(1).
Components of stress tensors and pressure are determined via displacement velocities as follows:
{dot over (p)}.sup.(1)=.sub.1 div u.sup.(1)+.sub.2 div v.sup.(1),(I.31)
{dot over ()}.sub.rr.sup.(1)=.sub.3 div u.sup.(1)+.sub.4 div v.sup.(1)+2u.sub.r,(I.32)
{dot over ()}.sub.rr.sup.(0)={dot over (p)}.sup.(0)=.sub.0 div v.sup.(0),(I.33)
where the following is determined:

(48) 0$\begin{array}{cc}{}_1=\left(\frac{{}_s}{}-\frac{{}_l}{}K\right)/\left({}_s{c}_t^2\right),{}_2=\left(\frac{{}_s}{}-\frac{{}_l}{}K\right)/\left({}_s{c}_t^2\right),& \left(I\mathrm{.34}\right)\\ {}_3=\left(\frac{{}_s}{}-\frac{2}{3}\right)/\left({}_s{c}_t^2\right),{}_4=\left(\frac{{}_l}{}\right)/\left({}_s{c}_t^2\right),& \left(I\mathrm{.35}\right)\\ {}_0={}^{\left(0\right)}{c}_{p0}^2/\left({}_s{c}_t^2\right).& \left(I\mathrm{.36}\right)\end{array}$

(49) The set of equations (I.23, I.24, I.26) with boundary conditions (I.27-I.30) for geometry shown in FIG. 1 has a solution that is expressed via the Hankel functions H.sub..sup.(1)(z), H.sub..sup.(2)(z). In the saturated porous medium, this solution in a dimensionless form is

(50) for pressure:
p.sup.(0)iB.sub.1(.sub.1+.sub.2R.sup.3)H.sub.1.sup.(1)(r/{square root over (.sub.1)})/{square root over (.sub.1)}iB.sub.2(.sub.1R.sub.2+.sub.2)H.sub.1.sup.(1)(r/{square root over (.sub.2)})/{square root over (.sub.2)},(I.37) and
for the solid matrix velocities and velocities of the saturating fluid:
u.sub.r.sup.(1)=B.sub.1H.sub.1.sup.(1)(r/{square root over (.sub.1)})+B.sub.2R.sub.2H.sub.1.sup.(1)(r/{square root over (.sub.2)}),(I.38)
v.sub.r.sup.(1)=B.sub.1R.sub.3H.sub.1.sup.(1)(r/{square root over (.sub.1)})+B.sub.2H.sub.1.sup.(1)(r/{square root over (.sub.2)}),(I.39)
and, consequently, for the difference of the matrix and fluid velocities:
u.sub.r.sup.(1)v.sub.r.sup.(1)=B.sub.1(1R.sub.3)H.sub.1.sup.(1)(r/{square root over (.sub.1)})+B.sub.2(R.sub.21)H.sub.1.sup.(1)(r/{square root over (.sub.2)}),(I.40)
and electrical field

(51) $\begin{array}{cc}{E}_r^{\left(1\right)}=-\frac{c_e}{c_t}\left(u_r^{\left(1\right)}-v_r^{\left(1\right)}\right).& \left(I\mathrm{.41}\right)\end{array}$

(52) In the wellbore, this solution in a dimensionless form is:

(53) for pressure in fluid
p.sup.(0)=iB.sub.0.sub.0H.sub.1.sup.(1)(r/{square root over (.sub.0)})/{square root over (.sub.0)}iC.sub.0.sub.0H.sub.1.sup.(2)(r/{square root over (.sub.0)})/{square root over (.sub.0)},(I.42)
for the fluid velocity
v.sub.r.sup.(0)=B.sub.0H.sub.1.sup.(1)(r/{square root over (.sub.0)})+C.sub.0H.sub.1.sup.(2)(r/{square root over (.sub.0)}).(I.43)
Integration constants B.sub.1, B.sub.2, B.sub.0, C.sub.0 are determined by the condition that boundary problem solutions (I.20, A.21, A.23) are to satisfy boundary conditions (I.24-I.27). Substituting solutions (I.34-I.40) into boundary conditions (I.24-I.27), we arrive at this set:
G=P,(I.44)
where one can see the column vector of the integration constants , the source part P, and the coefficient matrix G:

(54) $\begin{array}{cc}=\left(\begin{array}{c}{B}_1\\ B_2\\ B_0\\ C_0\end{array}\right),P=\left(\begin{array}{c}{\stackrel{\_}{P}}_0\\ 0\\ 0\\ 0\end{array}\right),G=\left(\begin{array}{cccc}{g}_{11}& g_{12}& g_{13}& g_{14}\\ g_{21}& g_{22}& g_{23}& g_{24}\\ g_{31}& g_{32}& g_{33}& g_{34}\\ g_{41}& g_{42}& g_{43}& g_{44}\end{array}\right),& \left(I\mathrm{.45}\right)\end{array}$
Where
g.sub.11=0, g.sub.12=0,(I.46)
g.sub.13=.sub.0H.sub.0.sup.(1)(r.sub.1/{square root over (.sub.0)})/{square root over (.sub.0)}, g.sub.14=.sub.0H.sub.0.sup.(2)(r.sub.1/{square root over (.sub.0)})/{square root over (.sub.0)},(I.47)
g.sub.21=((1)+R.sub.3)H.sub.1.sup.(1)(r.sub.2/{square root over (.sub.1)}), g.sub.22=((1)R.sub.2+)H.sub.1.sup.(1)(r.sub.2/{square root over (.sub.2)}),(I.48)
g.sub.23=H.sub.1.sup.(1)(r.sub.2/{square root over (.sub.0)}), g.sub.24=H.sub.1.sup.(2)(r.sub.2/{square root over (.sub.0)}),(I.49)
g.sub.31=(.sub.3+.sub.4R.sub.3+2)H.sub.0.sup.(1)(r.sub.2/{square root over (.sub.1)})/{square root over (.sub.1)}2H.sub.1.sup.(1)(r.sub.2/{square root over (.sub.1)})/r.sub.2,(I.50)
g.sub.32=(.sub.3R.sub.2+.sub.4+2)H.sub.0.sup.(1)(r.sub.2/{square root over (.sub.2)})/{square root over (.sub.2)}2R.sub.2H.sub.1.sup.(1)(r.sub.2/{square root over (.sub.2)})/r.sub.2,(I.51)
g.sub.33=.sub.0H.sub.0.sup.(1)(r.sub.2/{square root over (.sub.0)})/{square root over (.sub.0)}, g.sub.34=.sub.0H.sub.0.sup.(2)(r.sub.2/{square root over (.sub.0)})/{square root over (.sub.0)},(I.52)
g.sub.41=((.sub.1+.sub.2R.sub.3).sub.1+ib(1R.sub.3))H.sub.0.sup.(1)(r.sub.2/{square root over (.sub.1)})/{square root over (.sub.1)},(I.53)
g.sub.42=((.sub.1R.sub.2+.sub.2).sub.2+ib(R.sub.21))H.sub.0.sup.(1)(r.sub.2/{square root over (.sub.2)})/{square root over (.sub.2)},(I.54)
g.sub.43=(.sup.(1)/.sup.(0)).sub.0H.sub.0.sup.(1)(r.sub.2/{square root over (.sub.0)})/{square root over (.sub.0)},(I.55)
g.sub.44=(.sup.(1)/.sup.(0)).sub.0H.sub.0.sup.(2)(r.sub.2/{square root over (.sub.0)})/{square root over (.sub.0)},(I.56)

(55) The solution to set (I.44) yields integration constants B.sub.1, B.sub.2, B.sub.0, C.sub.0:
=G.sup.1P.(I.57)
It should be noted that pressure is in the units of the elastic modulus =.sub.sc.sub.t.sup.2, and velocities u, v are in the units of the shear sound velocity c.sub.t.

APPENDIX II

(56) The theory describing propagation of the guided Stoneley waves at the interface of the fluid and the deformed saturated porous medium is, as in Appendix I, based on a linearized version of the filtration theory.

(57) The acoustic and electrical fields in fluid are described by the following linear equation:

(58) $\begin{array}{cc}\stackrel{\mathrm{.Math.}}{v}-c_{p0}^2\mathrm{div}v=0,& \left(\mathrm{II}\mathrm{.1}\right)\\ \stackrel{.}{B}-\frac{c_e^2}{4}B=0,& \left(\mathrm{II}\mathrm{.2}\right)\\ E=\frac{c_e}{4}\mathrm{rot}B.& \left(\mathrm{II}\mathrm{.3}\right)\end{array}$
The acoustic and electrical fields in the saturated porous medium are described by the following set of linear equations:

(59) $\begin{array}{cc}\stackrel{\mathrm{.Math.}}{u}-c_t^{2}u-a_1\mathrm{div}u+a_2\mathrm{div}v+\frac{{}_l}{{}_s}{b}_{}\left(\stackrel{.}{u}-\stackrel{.}{v}\right)+\frac{{}_l}{{}_s}{b}_e\mathrm{rot}B=0,& \left(\mathrm{II}\mathrm{.4}\right)\\ \stackrel{\mathrm{.Math.}}{v}+a_3\mathrm{div}u-a_4\mathrm{div}v-b_{}\left(\stackrel{.}{u}-\stackrel{.}{v}\right)-b_e\mathrm{rot}B=0,& \left(\mathrm{II}\mathrm{.5}\right)\\ \stackrel{.}{B}-\frac{c_e^2}{4}B-\frac{c_e{}_l}{}\mathrm{rot}\left(u-v\right)=0,& \left(\mathrm{II}\mathrm{.6}\right)\\ E=\frac{c_e}{4}\mathrm{rot}B-\frac{{}_l}{}\left(u-v\right).& \left(\mathrm{II}\mathrm{.7}\right)\end{array}$

(60) Here .sub.l, .sub.s, are partial densities of the saturating fluid and the porous matrix, respectively, =.sub.l+.sub.s is density of the saturated porous medium, u is the velocity of the porous matrix, v is the velocity of the saturating fluid or the borehole fluid, c.sub.p0 is the velocity of sound in the borehole fluid. The dissipative coefficients are b.sub.=.sub.l(1.sup.2/()), =/(k.sub.l), b.sub.e=c.sub.e/(4), where is dynamic viscosity the saturating fluid, k is permeability, is the electroacoustic constant of the porous medium, is electric conductivity, c.sub.e is the velocity of light. Coefficients .sub.j in the equations above are determined by three elastic moduli , , :

(61) $\begin{array}{cc}{a}_1=\frac{1}{{}_s}\left(\frac{{}_s^2}{{}^2}+\frac{{}_l^2}{{}^2}K+\frac{1}{3}\right),a_2=\frac{{}_l}{{}_s}\left(\frac{{}_l}{{}^2}K-\frac{{}_s}{{}^2}\right),& \left(\mathrm{II}\mathrm{.8}\right)\\ a_3=\frac{{}_l}{{}^2}K-\frac{{}_s}{{}^2},a_4=\frac{{}_l}{{}^2}K+\frac{{}_l}{{}^2}\mathrm{where}K=+\frac{2}{3}.& \left(\mathrm{II}\mathrm{.9}\right)\end{array}$

(62) Three elastic moduli , , and are found via the acoustic velocities of the saturated porous medium: the first compressional velocity c.sub.p1, the second compressional velocity c.sub.p2, shear velocity c.sub.t, measured at high frequency:
=.sub.sc.sub.t.sup.2,(II.10)

(63) $\begin{array}{cc}K=\frac{1}{2}\frac{{}_s}{{}_l}\left(c_{p1}^2+c_{p2}^2-\frac{8}{3}{}_l{c}_t^2-\sqrt{{\left(c_{p1}^2-{}_{p2}^2\right)}^2-\frac{64}{9}{}_s{}_l{c}_t^4}\right),& \left(\mathrm{II}\mathrm{.11}\right)\\ =\frac{1}{2}\left(c_{p1}^2+c_{p2}^2-\frac{8}{3}{}_s{c}_t^2+\sqrt{{\left(c_{p1}^2-{}_{p2}^2\right)}^2-\frac{64}{9}{}_s{}_l{c}_t^4}\right).& \left(\mathrm{II}\mathrm{.12}\right)\end{array}$
These are considered known quantities.

(64) Let us present equations (II.1-II.7) in the dimensionless form, like in Appendix I:

(65) in the borehole fluid:

(66) $\begin{array}{cc}{}_0\mathrm{div}v+v=0,& \left(\mathrm{II}\mathrm{.13}\right)\\ \mathrm{ig}B-B=0,& \left(\mathrm{II}\mathrm{.14}\right)\\ E=\frac{c_t}{c_e}g\mathrm{rot}B& \left(\mathrm{II}\mathrm{.15}\right)\end{array}$
in the saturated porous medium:

(67) $\begin{array}{cc}u+{\stackrel{\_}{a}}_1\mathrm{div}u-{\stackrel{\_}{a}}_2\mathrm{div}v+\left(1+\mathrm{.Math.}{\stackrel{\_}{b}}_{}\right)u-\mathrm{.Math.}{\stackrel{\_}{b}}_{}v+\mathrm{.Math.}{\stackrel{\_}{b}}_e\mathrm{rot}B=0,& \left(\mathrm{II}\mathrm{.16}\right)\\ -{\stackrel{\_}{a}}_3\mathrm{div}u+{\stackrel{\_}{a}}_4\mathrm{div}v-{\stackrel{\_}{b}}_{}u+\left(1+{\stackrel{\_}{b}}_{}\right)v-{\stackrel{\_}{b}}_e\mathrm{rot}B=0,& \left(\mathrm{II}\mathrm{.17}\right)\\ \mathrm{ig}B+i\mathrm{rot}\left(u-v\right)-B=0,& \left(\mathrm{II}\mathrm{.18}\right)\\ E=\frac{c_t}{c_e}g\mathrm{rot}B-\frac{c_t}{c_e}\left(u-v\right).& \left(\mathrm{II}\mathrm{.19}\right)\end{array}$
Here, velocities u, v are in units of the shear velocity of sound c.sub.t; coordinatesin units of c.sub.t/; Bin units of c.sub.e.sub.l/. Dimensionless coefficients in these equations are: .sub.0=c.sub.p0.sup.2/c.sub.t.sup.2, .sub.i=a.sub.i/c.sub.t.sup.2, b.sub.=ib.sub./, =.sub.l/.sub.s, g=c.sub.e.sup.2/(4 c.sub.t.sup.2), b.sub.e=i.sup.2.sub.lc.sub.e.sup.2/(4 .sup.2c.sub.t.sup.2).

(68) To study dispersive properties of the Stoneley waves, one needs the dispersion ratio. To obtain it, let us introduce potentials which will enable us to separate compressional and shear waves in the medium under consideration, with the help of the following equations:
u=.sub.urot .sub.u,(II.20)
v=.sub.vrot .sub.v.(II.21)
To study the Stoneley waves with the source in the center of the wellbore, it suffices to keep the cylindrical symmetry of the problem. The hydrodynamic velocities of the media should have components u=(u.sub.r, 0, u.sub.z) and v=(v.sub.r, 0, v.sub.z), while vector potentials have only a single component .sub.u=(0, .sub.u, 0) and .sub.v=(0, .sub.v, 0). The magnetic field should have component B=(0, B.sub., 0).

(69) In cylindrical coordinates, equations (II.20, II.21) have the following form:

(70) $\begin{array}{cc}{u}_r={}_r{}_u+{}_z{}_u,u_z={}_z{}_u-{}_r{}_u-\frac{1}{r}{}_u,& \left(\mathrm{II}\mathrm{.22}\right)\\ v_r={}_r{}_v+{}_z{}_v,v_z={}_z{}_v-{}_r{}_v-\frac{1}{r}{}_v.& \left(\mathrm{II}\mathrm{.23}\right)\end{array}$
Two independent sub-systems describing compressional and shear waves in the saturated porous medium are obtained via substituting (II.22, II.23) in (II.16-II.18):
(.sub.1+1).sub.u.sub.2.sub.v+(1+b.sub.).sub.ub.sub..sub.v=0,(II.24)
.sub.3.sub.u+.sub.4.sub.vb.sub..sub.u+(1+b.sub.).sub.v=0,(II.25)
.sub..sub.u+(1+b.sub.).sub.ub.sub..sub.vb.sub.eB.sub.=0,(II.26)
b.sub..sub.u+(1+b.sub.).sub.v+b.sub.eB.sub.=0,(II.27)
ig.sub.B.sub.+i.sub..sub.ui.sub..sub.vB.sub.=0.(II.28)

(71) Here we define the operators

(72) 0$={}_r^2+\frac{1}{r}{}_r+{}_z^2\mathrm{and}$${}_{}{}_i={}_r^2{}_i+\frac{1}{r}{}_r{}_i-\frac{1}{r^2}{}_i+{}_z^2{}_i.$
Along the z axis, for the symmetry assumed, the Stoneley wave is a plane wave:
(.sub.u,.sub.v,.sub.u,.sub.v,B).fwdarw.(.sub.u,.sub.v,.sub.u,.sub.v,B).Math.exp(ikzit)(II.29)
Let us find common solutions to equations (II.1, II.2, II.4-II.6).

(73) Let us consider set of equations (II.24, II.25) determining the evolution of the compressional waves in the saturated porous medium:
(.sub.1+1).sub.u.sub.2.sub.v+(1+b.sub.).sub.ub.sub..sub.v=0,(II.30)
.sub.3.sub.u+.sub.4.sub.vb.sub..sub.u+(1+b.sub.).sub.v=0,(II.31)
Taking into account the form of the solution, set of equations (II.30-II.31) determining the evolution of the compressional waves may be written out in the matrix form:
A+W=0.(II.32)
Here

(74) $A=\left(\begin{array}{cc}{\stackrel{\_}{a}}_0& -{\stackrel{\_}{a}}_2\\ -{\stackrel{\_}{a}}_3& {\stackrel{\_}{a}}_4\end{array}\right)$
is the matrix of coefficients,

(75) ${\stackrel{\_}{a}}_0={\stackrel{\_}{a}}_1+1,W=\left(\begin{array}{cc}\left(1+\mathrm{.Math.}{\stackrel{\_}{b}}_a\right)& -\mathrm{.Math.}{\stackrel{\_}{b}}_a\\ -{\stackrel{\_}{b}}_a& \left(1+{\stackrel{\_}{b}}_a\right)\end{array}\right),$
whereas the other coefficients are defined above; =(.sub.u .sub.v).sup.T is the vector-column of functions.

(76) We transform set of equations (II.32) via multiplying it by the matrix W.sup.1:

(77) $W^{-1}=\left(\begin{array}{cc}\frac{1+{\stackrel{\_}{b}}_a}{\det W}& \frac{\mathrm{.Math.}{\stackrel{\_}{b}}_a}{\det W}\\ \frac{{\stackrel{\_}{b}}_a}{\det W}& \frac{1+\mathrm{.Math.}{\stackrel{\_}{b}}_a}{\det W}\end{array}\right),$
where det W=1+(1+)b.sub., and we arrive at the set of equations with the diagonal matrix in front of the free term:
H+=0, where H=W.sup.1A.(II.33)
The components of the matrix H are:
h.sub.11=((1+b).sub.0b.sub.3)/det W, h.sub.12=((1+b).sub.2+b.sub.4)/det W,(II.34)
h.sub.21=(b.sub.0(1+b).sub.3)/det W, h.sub.22=(b.sub.2+(1+b).sub.4)/det W.(II.35)
Let us bring set of equations (II.33) obtained above to the diagonal form:

(78) $\begin{array}{cc}+=0,\mathrm{where}& \left(\mathrm{II}\mathrm{.36}\right)\\ =\left(\begin{array}{cc}{}_1& 0\\ 0& {}_2\end{array}\right),={\left({}_1{}_2\right)}^T& \left(\mathrm{II}\mathrm{.37}\right)\end{array}$
using the transform matrix R

(79) $\begin{array}{cc}R=\left(\left(\begin{array}{c}{R}_1\\ R_3\end{array}\right)\left(\begin{array}{c}{R}_2\\ R_4\end{array}\right)\right),& \left(\mathrm{II}\mathrm{.38}\right)\end{array}$
which is made of the right-hand side eigenvectors of the matrix A: R.sup.1AR=. Here
R.sub.1=1, R.sub.2=(.sub.2h.sub.22)/h.sub.21, R.sub.3=h.sub.21/(.sub.1h.sub.22), R.sub.4=1.(II.39)

(80) The column vector of these functions is transformed as follows:

(81) $\begin{array}{cc}{R}^{-1}=,\mathrm{where}=\left(\begin{array}{c}{}_1\\ {}_2\end{array}\right)& \left(\mathrm{II}\mathrm{.40}\right)\end{array}$
The diagonal matrix =R.sup.1AR has eigenvalues .sub.1 and .sub.2 of the matrix A as its diagonal elements, which are determined as follows:
.sub.1=tr H{square root over ((tr H).sup.2det H)}, .sub.2=tr H+{square root over ((tr H).sup.2det H)},(II.41)
where tr H=h.sub.11+h.sub.22 is the trace of the matrix H, det H=h.sub.11h.sub.22h.sub.12h.sub.21 is the determinant of the matrix H.

(82) Thus, set of equations (II.30, II.31) in cylindrical coordinates is expressed in the diagonal form, (II.29)

(83) $\begin{array}{cc}{}_r^2{}_1+\frac{1}{r}{}_r{}_1+\frac{1}{{}_1}{}_1+{}_z^2{}_1=0,& \left(\mathrm{II}\mathrm{.42}\right)\\ {}_r^2{}_2+\frac{1}{r}{}_r{}_2+\frac{1}{{}_2}{}_2+{}_z^2{}_2=0,& \left(\mathrm{II}\mathrm{.43}\right)\end{array}$
or, taking into account the form of solution (II.29),

(84) $\begin{array}{cc}{}_1^{}+\frac{1}{r}{}_1^{}-l_1^2{}_1=0,& \left(\mathrm{II}\mathrm{.44}\right)\\ {}_2^{}+\frac{1}{r}{}_2^{}-l_2^2{}_2=0,& \left(\mathrm{II}\mathrm{.45}\right)\end{array}$
where it is determined that

(85) $l_1^2=k^2-\frac{1}{{}_1},l_2^2=k^2-\frac{1}{{}_2}.$
The prime sign marks the derivative with respect to radius.

(86) Solutions to equations (II.44, II.45) are expressed via the modified Bessel functions K.sub.0(l.sub.ir), I.sub.0(l.sub.ir) and, taking into account the boundedness of this solution in the infinity, they have the following form:
.sub.1=C.sub.1K.sub.0(l.sub.1r), .sub.2=C.sub.2K.sub.0(l.sub.2r)(II.46)
Here C.sub.1, C.sub.2 are integration constants.
The solution to set (II.44, II.45) and the solution to the original set (II.30, II.31) are related as in (II.40):
=R(II.47)
Or
.sub.u=.sub.1+R.sub.2.sub.2,(II.48)
.sub.v=R.sub.3.sub.1+.sub.2.(II.49)
Thus, the solution for the potentials .sub.u, .sub.v can be written out as follows:
.sub.u=C.sub.1K.sub.0(l.sub.1r)+C.sub.2R.sub.2K.sub.0(l.sub.2r), .sub.v=C.sub.1R.sub.3K.sub.0(l.sub.1r)+C.sub.2K.sub.0(l.sub.2r). (II.50)

(87) Let us next consider the set of equations (II.26-II.28) determining the evolution of the shear waves:
.sub..sub.u+(1+b.sub.).sub.ub.sub..sub.vb.sub.eB.sub.=0,(II.51)
b.sub..sub.u+(1+b.sub.).sub.v+b.sub.eB.sub.=0,(II.52)
ig.sub.B.sub.+i.sub..sub.ui.sub..sub.vB.sub.=0.(II.53)
The second equation yields the following relationship between the potentials:

(88) 0$\begin{array}{cc}{}_v={\stackrel{\_}{b}}_a{}_u-{\stackrel{\_}{b}}_e{B}_{},\mathrm{where}=\frac{1}{1+{\stackrel{\_}{b}}_a},& \left(\mathrm{II}\mathrm{.54}\right)\end{array}$
which allows us to exclude .sub.v from set (II.52):
.sub..sub.u+b.sub..sub.ub.sub.eB.sub.=0,(II.55)
i.sub..sub.u+i(g+b.sub.e).sub.B.sub.B.sub.=0.(II.56)

(89) Taking into account the form of the solution, set of equations (II.55, II.56) determining the evolution of the compressional waves may be written out in the matrix form:
D.sub.+W=0.(II.57)
Here

(90) $D=\left(\begin{array}{cc}1& 0\\ i& i\left(g+{\stackrel{\_}{b}}_e\right)\end{array}\right)$
is the matrix of coefficients,

(91) $W=\left(\begin{array}{cc}\mathrm{.Math.}{\stackrel{\_}{b}}_a& -\mathrm{.Math.}{\stackrel{\_}{b}}_e\\ 0& -1\end{array}\right),$
while the other coefficients are defined above; =(.sub.u B.sub.).sup.T is the vector-column of functions.
We transform set of equations (II.57) via multiplying it by the matrix W.sup.1:

(92) $W^{-1}=\left(\begin{array}{cc}\frac{-1}{\det W}& \frac{\mathrm{.Math.}{\stackrel{\_}{b}}_e}{\det W}\\ 0& \frac{\mathrm{.Math.}{\stackrel{\_}{b}}_{}}{\det W}\end{array}\right),$
where det W=b.sub.e, and we arrive at the set of equations with the diagonal matrix in front of the free term:
N.sub.+=0, where N=W.sup.1D.(II.58)
The components of the matrix N are:
n.sub.11=(1+b.sub.e.sup.2)/det W, n.sub.12=ib.sub.e(g+b.sub.e)/det W,(II.59)
n.sub.21=ib.sub..sup.2/det W, n.sub.22=ib.sub.(g+b.sub.e)/det W.(II.60)

(93) Let us bring set of equations (II.58) obtained above to the diagonal form:

(94) $\begin{array}{cc}{}_{}+=0,\mathrm{where}& \left(\mathrm{II}\mathrm{.61}\right)\\ =\left(\begin{array}{cc}{}_3& 0\\ 0& {}_4\end{array}\right),={\left({}_1{}_2\right)}^T& \left(\mathrm{II}\mathrm{.62}\right)\end{array}$
using the transform matrix P

(95) $\begin{array}{cc}P=\left(\begin{array}{cc}{P}_1& P_2\\ P_3& P_4\end{array}\right),& \left(\mathrm{II}\mathrm{.63}\right)\end{array}$
which is made of the right-hand side eigenvectors of the matrix D: P.sup.1DP=. Here
P.sub.1=1, P.sub.2=(.sub.4n.sub.22)/n.sub.21, P.sub.3=n.sub.21/(.sub.3n.sub.22), P.sub.4=1.(II.64)
The column vector of these functions is transformed as follows:

(96) $\begin{array}{cc}{P}^{-1}=,\mathrm{where}=\left(\begin{array}{c}{}_1\\ {}_2\end{array}\right).& \left(\mathrm{II}\mathrm{.65}\right)\end{array}$

(97) The diagonal matrix =P.sup.1DP has eigenvalues .sub.3 and .sub.4 of the matrix A as its diagonal elements, which are determined as follows:
.sub.3=tr N{square root over ((tr N).sup.2det N)}, .sub.4=tr N+{square root over ((tr N).sup.2detN)},(II.66)
where tr N=n.sub.11+n.sub.22, det N=n.sub.11n.sub.22n.sub.12n.sub.21.
Thus, set of equations (II.55, II.56) in cylindrical coordinates is expressed in the diagonal form, (II.61)

(98) $\begin{array}{cc}{}_r^2{}_1+\frac{1}{r}{}_r{}_1-\frac{1}{r^2}{}_1+\frac{1}{{}_3}{}_1+{}_z^2{}_1=0,& \left(\mathrm{II}\mathrm{.67}\right)\\ {}_r^2{}_2+\frac{1}{r}{}_r{}_2-\frac{1}{r^2}{}_2+\frac{1}{{}_4}{}_2+{}_z^2{}_2=0,& \left(\mathrm{II}\mathrm{.68}\right)\end{array}$
or, taking into account the form of its solution (II.29), set of equations (II.67, II.68) determining the evolution of the shear waves can be reduced to the following set of equations:

(99) $\begin{array}{cc}{}_1^{}+\frac{1}{r}{}_1^{}-\left(l_3^2+\frac{1}{r^2}\right){}_1=0,& \left(\mathrm{II}\mathrm{.69}\right)\\ {}_2^{}+\frac{1}{r}{}_2^{}-\left(l_4^2+\frac{1}{r^2}\right){}_2=0,& \left(\mathrm{II}\mathrm{.70}\right)\end{array}$
where it is determined that

(100) $l_3^2=k^2-\frac{1}{{}_3},l_4^2=k^2-\frac{1}{{}_4}.$
The prime sign marks the derivative with respect to radius.
Solutions to equations (II.67, II.68) are expressed via the modified Bessel functions K.sub.1(l.sub.ir), I.sub.1(l.sub.ir) and, taking into account the boundedness of this solution in the infinity, they have the following form:
.sub.1=C.sub.3K.sub.1(l.sub.3r), .sub.2=C.sub.4K.sub.1(l.sub.4r).((II.71)
Here C.sub.3, C.sub.4 are integration constants.
The solution to set (II.67, II.68) and the solution to the original set (II.55, II.56) are related as in (II.65):
=P(II.72)
Or
.sub.u=.sub.1+P.sub.2.sub.2,(II.73)
B.sub.=P.sub.3.sub.1+.sub.2.(II.74)

(101) The solution for the potentials .sub.u, .sub.v and magnetic field B.sub. are expressed via the solution for the potentials .sub.1, .sub.2:
.sub.u=C.sub.3K.sub.1(l.sub.3r)+C.sub.4P.sub.2K.sub.1(l.sub.4r),(II.75)
B.sub.=C.sub.3P.sub.3K.sub.1(l.sub.3r)+C.sub.4K.sub.1(l.sub.4r),(II.76)
.sub.v=C.sub.3.sub.3K.sub.0(l.sub.3r)+C.sub.4.sub.4K.sub.0(l.sub.4r),(II.77)
where .sub.3=(b.sub.b.sub.eP.sub.3), .sub.4=(b.sub.P.sub.2b.sub.e).

(102) Solutions (II.75, I.76, I.77) obtained above lead to expressions for velocities and time derivatives of the stress tensors. Due to the choice of geometry in this problem, velocity components can be expressed as follows:

(103) 0$\begin{array}{cc}\begin{array}{cc}{u}_r={}_r{}_u+{}_z{}_u,& u_z={}_z{}_u-{}_r{}_u-\frac{1}{r}{}_u,\end{array}& \left(\mathrm{II}\mathrm{.78}\right)\\ \begin{array}{cc}{v}_r={}_r{}_v+{}_z{}_v,& v_z={}_z{}_v-{}_r{}_v-\frac{1}{r}{}_v.\end{array}& \left(\mathrm{II}\mathrm{.79}\right)\end{array}$
Because, from (II.48, II.49, II.73, II.74):
.sub.u=.sub.1+R.sub.2.sub.2, .sub.v=R.sub.3.sub.1+.sub.2,(II.80)
.sub.u=.sub.1+P.sub.2.sub.2, B.sub.=P.sub.3.sub.1+.sub.2, .sub.v=.sub.3.sub.1+.sub.4.sub.2,(II.81)
we have

(104) $\begin{array}{cc}{u}_r={}_r{}_1+R_2{}_r{}_2+{}_z{}_1+P_2{}_z{}_2,& \left(\mathrm{II}\mathrm{.82}\right)\\ u_z={}_z{}_1+R_2{}_z{}_2-{}_r{}_1-P_2{}_r{}_2-\frac{1}{r}{}_1-\frac{1}{r}{P}_2{}_2,& \left(\mathrm{II}\mathrm{.83}\right)\\ v_r=R_3{}_r{}_1+{}_r{}_2+{}_3{}_z{}_u+{}_4{}_z{}_{},& \left(\mathrm{II}\mathrm{.84}\right)\\ v_z=R_3{}_z{}_1+{}_z{}_2-{}_3{}_r{}_v-{}_4{}_r{}_{}-\frac{1}{r}{}_3{}_u-\frac{1}{r_e}{}_4{}_{}.& \left(\mathrm{II}\mathrm{.85}\right)\end{array}$

(105) In accordance to (II.50, II.75-I.77) and taking into account the recurrent relationships between the Bessel functions, we arrive at the final form of the solution to set of equations (II.16-II.18):
u.sub.r=C.sub.1l.sub.1K.sub.1(l.sub.1r)C.sub.2R.sub.2l.sub.2K.sub.1(l.sub.2r)+ikC.sub.3K.sub.1(l.sub.3r)+ikC.sub.4P.sub.2K.sub.1(l.sub.4r),(II.86)
u.sub.z=ikC.sub.1K.sub.0(l.sub.1r)+ikC.sub.2R.sub.2K.sub.0(l.sub.2r)+C.sub.3l.sub.3K.sub.0(l.sub.3r)+C.sub.4P.sub.2l.sub.4K.sub.0(l.sub.4r),(II.87)
v.sub.r=C.sub.1R.sub.3l.sub.1K.sub.1(l.sub.1r)C.sub.2l.sub.2K.sub.1(l.sub.2r)+ikC.sub.3.sub.3K.sub.1(l.sub.3r)+ikC.sub.4.sub.4K.sub.1(l.sub.4r),(II.88)
v.sub.z=ikC.sub.1R.sub.3K.sub.0(l.sub.1r)+ikC.sub.2K.sub.0(l.sub.2r)+C.sub.3.sub.3l.sub.4K.sub.0(l.sub.4r),(II.89)
B.sub.=C.sub.3P.sub.3K.sub.1(l.sub.3r)+C.sub.4K.sub.1(l.sub.4r)(II.90)

(106) Let us next consider the equation determining the evolution of the borehole fluid. We can introduce potentials for the fluid velocity just as in the discussion above, via the following expression:
v=.sub.frot .sub.f,(II.91)
which has the following components in the given geometry:
v.sub.r=.sub.r.sub.f,(II.92)
v.sub.z=.sub.z.sub.f.(II.93)
Then, in accordance to equation (II.13), we arrive at:
.sub.0.sub.f+.sub.f=0.(II.94)

(107) The solution to equation (II.94) for the borehole fluid is expressed via the modified Bessel functions I.sub.(z):

(108) $\begin{array}{cc}{}_v=C_5{I}_0\left(l_{0}r\right),l_0^2=k^2-\frac{1}{{}_0}.& \left(\mathrm{II}\mathrm{.95}\right)\end{array}$
Here we take into account the fact that no singularity is found in the center of the borehole.

(109) Let us next consider the equation determining the evolution of magnetic field in the borehole (II.14):
ig.sub.B.sub.B.sub.=0,(II.96)
The solution to equation (II.96) for the borehole fluid is expressed via the modified Bessel functions I.sub.(z):

(110) $\begin{array}{cc}{B}_{}=C_6{I}_1\left(l_{6}r\right),l_6^2-\frac{i}{g}.& \left(\mathrm{II}\mathrm{.97}\right)\end{array}$
Here we take into account the fact that no singularity is found in the center of the borehole. The solutions to equations (II.13, II.14) for the borehole are:
v.sub.r=C.sub.5l.sub.0I.sub.1(l.sub.0r),(II.98)
v.sub.z=ikC.sub.5I.sub.0(l.sub.0r),(II.99)
B.sub.=C.sub.6I.sub.1(l.sub.6r).(II.100)
We arrive at the dispersive relationship for the wave vector k in accordance to the condition of solutions (II.86-II.90, II.98-II.100) satisfying the boundary conditions at the borehole boundary:

(111) The continuity condition for the component of the magnetic field and z component of the electrical field:
B.sub..sup.(1)=B.sub..sup.(0)(II.101)
E.sub.z.sup.(1)=E.sub.z.sup.(0)(II.102)
The continuity condition for the normal mass flow:
(1)u.sub.r.sup.(1)+v.sub.r.sup.(1)=v.sub.r.sup.(0)(II.103)
The continuity condition for the components of the normal projection of the stress tensor .sub.ik
.sub.rr.sup.(1)=.sub.rr.sup.(0),(II.104)
.sub.rz.sup.(1)=0,(II.105)

(112) The equality condition for the components of the normal projection of the partial stress tensor of the saturating fluid .sub.rr.sup.l(1) and the borehole fluid .sub.rr.sup.(0)
.sub.rr.sup.l(1)=.sub.rr.sup.(0).(II.106)
The latter condition in terms of pressure takes the following form:

(113) $\begin{array}{cc}\frac{p^{\left(1\right)}}{{}^{\left(1\right)}}-\frac{{}^{\left(1\right)}}{{}_s^{\left(1\right)}}\stackrel{\_}{b}\left({}_u^{\left(1\right)}-{}_v^{\left(1\right)}\right)=\frac{p^{\left(0\right)}}{{}^{\left(0\right)}}.& \left(\mathrm{II}\mathrm{.107}\right)\end{array}$
Here the coefficient is introduced because there is no reason to believe that the friction coefficient at the boundary has to coincide with the friction coefficient in the porous medium. Consequently, the continuity condition is satisfied for the time derivatives of the components of the normal projection of the elastic tensor:
{dot over ()}.sub.rz.sup.(1)=0, {dot over ()}.sub.rr.sup.(1)={dot over ()}.sub.rr.sup.(0), {dot over ()}.sub.rr.sup.(1)={dot over ()}.sub.rr.sup.(0).(II.108)
Here
{dot over (p)}.sup.(1)=.sub.1 div u.sup.(1)+.sub.2 div v.sup.(1),(II.109)
{dot over ()}.sub.rr.sup.(1)=.sub.3 div u.sup.(1)+.sub.4 div v.sup.(1)+2.sub.ru.sub.r,(II.110)
{dot over ()}.sub.rz.sup.(1)=.sub.ru.sub.z+.sub.zu.sub.r,(II.111)

(114) $\begin{array}{cc}{\stackrel{.}{}}_{\mathrm{rr}}^{\left(0\right)}=-{\stackrel{.}{p}}^{\left(0\right)}={}_0{\mathrm{divv}}^{\left(0\right)},& \left(\mathrm{II}\mathrm{.112}\right)\\ {}_1=\left(\frac{{}_s}{}-\frac{{}_l}{}K\right)/\left({}_s{c}_t^2\right),{}_2\left(\frac{{}_l}{}+\frac{{}_l}{}K\right)/\left({}_s{c}_t^2\right),& \left(\mathrm{II}\mathrm{.113}\right)\\ {}_3=\left(\frac{{}_s}{}-\frac{2}{3}\right)/\left({}_s{c}_t^2\right),{}_4\left(\frac{{}_l}{}\right)/\left({}_s{c}_t^2\right),& \left(\mathrm{II}\mathrm{.114}\right)\\ {}_0={}^{\left(0\right)}{c}_{p0}^2/\left({}_s{c}_t^2\right).& \left(\mathrm{II}\mathrm{.115}\right)\end{array}$

(115) Substituting equations (II.93-II.100) into the boundary conditions (II.101-II.106), we can then present the boundary conditions in the matrix form:
G=0.(II.116)
where the vector-column of coefficients =(C.sub.1 C.sub.2 . . . C.sub.5 C.sub.6).sup.T and the matrix G with the components

(116) $\begin{array}{cc}{g}_{11}=\left(\left(1-\right)+R_3\right)l_1{K}_1\left(l_1{r}_2\right),g_{12}=\left(\left(1-\right)R_2+\right)l_2{K}_1\left(l_2{r}_2\right),& \left(\mathrm{II}\mathrm{.117}\right)\\ g_{13}=-k\left(\left(1-\right)+{}_3\right)K_1\left(l_3{r}_2\right),g_{14}=-k\left({}_4-\left(1-\right)P_2\right)K_1\left(l_3{r}_2\right),& \left(\mathrm{II}\mathrm{.118}\right)\\ g_{15}=-l_0{l}_1\left(l_0,r_2\right),g_{16}=0,& \left(\mathrm{II}\mathrm{.119}\right)\\ g_{21}=\left(\left({}_3+{}_4{R}_3\right)\left(l_1^2-k^2\right)+2l_1^2\right)K_0\left(l_1{r}_2\right)+2l_1{K}_1\left(l_1{r}_2\right)/r_2,& \left(\mathrm{II}\mathrm{.120}\right)\\ g_{22}=\left(\left({}_3{R}_2+{}_4\right)\left(l_2^2-k^2\right)+2R_2{l}_2^2\right)K_0\left(l_2{r}_2\right)+2R_2{l}_2{K}_2\left(l_2{r}_2\right)/r_2,& \left(\mathrm{II}\mathrm{.121}\right)\\ g_{23}=-2k\left(l_3{K}_0\left(l_3{r}_2\right)+K_1\left(l_3{r}_2\right)/r_2\right),& \left(\mathrm{II}\mathrm{.122}\right)\\ g_{24}=-2{\mathrm{kP}}_2\left(l_4{K}_0\left(l_4{r}_2\right)+K_1\left(l_4{r}_2\right)/r_2\right),& \left(\mathrm{II}\mathrm{.123}\right)\\ g_{25}=-{}_0\left(l_0^2-k^2\right)I_0\left(l_0{r}_2\right),g_{26}=0,& \left(\mathrm{II}\mathrm{.124}\right)\\ g_{31}=2{\mathrm{kl}}_1{K}_1\left(l_1{r}_2\right),g_{32}=2{\mathrm{kR}}_2{l}_2{K}_1\left(l_2{r}_2\right),& \left(\mathrm{II}\mathrm{.125}\right)\\ g_{33}=\left(l_3^2+k^2\right)K_1\left(l_3{r}_2\right),g_{34}=P_2\left(l_3^2+k^2\right)K_1\left(l_3{r}_2\right),& \left(\mathrm{II}\mathrm{.126}\right)\\ g_{35}=0,g_{36}=0,& \left(\mathrm{II}\mathrm{.127}\right)\\ g_{41}=\left(\left({}_1+R_3{}_2\right)\left(l_1^2-k^2\right)-{\stackrel{\_}{b}}_{}\frac{{}^{\left(1\right)}}{{}_s^{\left(1\right)}}\left(1-R_3\right)\right)K_0\left(l_1{r}_2\right),& \left(\mathrm{II}\mathrm{.128}\right)\\ g_{42}=\left(\left(R_2{}_1+{}_2\right)\left(l_2^2-k^2\right)-{\stackrel{\_}{b}}_{}\frac{{}^{\left(1\right)}}{{}_s^{\left(1\right)}}\left(R_2-1\right)\right)K_0\left(l_2{r}_2\right),& \left(\mathrm{II}\mathrm{.129}\right)\\ g_{43}=0,g_{44}=0,g_{45}=-\frac{{}^{\left(1\right)}}{{}^{\left(0\right)}}{}_0\left(l_0^2-k^2\right)I_0\left(l_0{r}_2\right),g_{46}=0,& \left(\mathrm{II}\mathrm{.130}\right)\\ g_{51}=k\left(1-R_3\right)K_0\left(l_1{r}_2\right),g_{52}=k\left(R_2-1\right)K_0\left(l_2{r}_2\right),& \left(\mathrm{II}\mathrm{.131}\right)\\ g_{53}=\left(1+{\mathrm{gP}}_3-{}_3\right)l_3{K}_0\left(l_3{r}_2\right),g_{54}=\left(P_2+g-{}_4\right)l_4{K}_0\left(l_4{r}_2\right),& \left(\mathrm{II}\mathrm{.132}\right)\\ g_{55}=0,g_{56}={\mathrm{gl}}_6{I}_0\left(l_6{r}_2\right),& \left(\mathrm{II}\mathrm{.133}\right)\\ g_{61}=0,g_{62}=0,g_{63}=P_3{K}_1\left(l_3{r}_2\right),g_{64}=K_1\left(l_4{r}_2\right),& \left(\mathrm{II}\mathrm{.134}\right)\\ g_{65}=0,g_{66}=-I_1\left(l_6{r}_2\right).& \left(\mathrm{II}\mathrm{.135}\right)\end{array}$