CA1096939A - Apparatus and method for determination of subsurface porosity - Google Patents

Abstract

60.453 ABSTRACT OF THE DISCLOSURE

The disclosure is directed to an improved apparatus and method for determining the water-filled porosity, designated ?.omega., of formations surrounding a borehole. A first quantity is obtained which is a measure of the travel time of microwave electromagnetic energy propagating over a predetermined distance in the formations, this first quantity being designated as tp1. A second quantity is also obtained which is a measure of the attenuation of the microwave electromagnetic energy propagating over the predetermined distance in the formations, the second quantity being the attenuation constant, .alpha.. A
third quantity, designated tpw1, is representative of the travel time of microwave electromagnetic energy propagating over the predetermined distance in water having a lossiness determined from ?. Finally, tp1 and tpw1 are combined to generate the water-filled porosity ?w.

Description

1~16~
BACKGROUWD OF THE INVENTION

This invention relates to determination of parameters of subsurface formations and, more particularly, to an apparatus and method for determining the water-filled porosity of formations surrounding a borehole.
A number of techniques presently exist for determining the porosity of subsurface formations.
Typically, electrical, sonic, and/or nuclear logging devices are lowered into a borehole and information from one or more of these devices is utilized to obtain a measure of the porosity of subsurface formations. Typically, readings or samples indicative of other subsurface parameters, for example lithology and saturation of fluids, is also determinea, and the determination of porosity is made using a combination of available information.
1~ There are a number of different porosity parameters, and it is generally useful to obtain information abou~ each one, if possible. Total porosity, designated ~t~ generally means the total fraction of formation unit volume which is occupied by pore spaces or interstices. Water-filled porosity, designated ~w' is the fraction of the formation unit volume occupied by water, and hydrocarbon porosity, designated ~h is the fraction of formation unit volume occupied by hydrocarbons. Generally, ~t equals the sum of ~w and ~hc ; ~t is evident that in situations when two of these three porosity parameters are known with reasonable reliability, valuable conc~usions can be drawn concerning the nature of the subsurface formations under investigation. It is also '~

:~Q~6~3~

known that comparison of porosity determinations from different logging devices, even ones which supposedly yield information concerning the same porosity parameter, can be highly useful in determining the nature of subsurface formations. Accordingly, it would be useful to have an additional reliable measure of a porosity parameter, such as ~w' which will typically be utilized in conjunction with other logging information in visual-izing the overall nature of subsurface formations traversed by a borehole.
In the U.S. Patent No. 3,944,910 of Rama Rau issued March 16, 1976 there is disclosed an apparatus for injecting microwave electromagnetic energy into subsurface formations to investigate properties of the formations. The Rau patent discloses, inter alia, a technique for determining the loss-corrected velocity (where velocity is inversely related to travel time) of microwave electromagnetic energy propagating through subsurface formations. For lossless media, the velocity of propagation of electromagnetic energy is inversely proportion-al to the square root of the dielectric constant. A velocity correction is generally made, however, when the propagating media exhibits lossiness.
It has been proposed that water-filled porosity of a subsurface formation, consisting of water in a matrix, could be determined by obtaining the loss-corrected velocity of the formation, and then utilizing a time-average relationship to determine the water-filled porosity. In particular, the following relationship has been set forth:

!
~cOrr = (1 ~ ~w) ~m + ~w ~ (1) &~39 where corr represents the loss-correc~ed dielectric constant of the formations, ~ is the dielectric constant of lossless water, and iS the dielectric constant of the formation matrix. ~ r can be obtained using the techni~ues set forth in the referenced Rau patent.
WO
is determinable, for a particular frequency and temperature, experimentally in the laboratory. If lithology is known, m of the matrix, generally assumed to be lossless, is also known. The square root of dielectric constant is inversely proportional to velocity, so the relationship (1) is in a time-average form. Relationship (1) indicates that the composite velocity is expressed as an average of the velocity in the water multiplied by the fractional volume of water plus the velocity in the matrix multiplied by the frac~ional volume of the matrix. A11 terms in the _elationship (1) consider the various propagation media (matrix, fluid, and composite) as being lossless, bearing in mind that the measured dielectric constant has been loss-corrected in this expression. Solving for ~ yields:
,_ ~w 2 ~ corr ~ m (2) ~ ~ ~ ' .
It is an object of the present invention to pro~ide an apparatus and method which yields a determination of water-fi~led porosity which is an improvement over prior art techniques, such as the one relating to expression (2).

. _ . _ . , .. .. . _ _~ . ... .. ....

6~39 SUMMARY OF THE INVENTION

Applicant has discovered that a possible disadvantage of the technique set forth in the Background section hereof is that the loss correction applied to obtain a measure of the loss-corrected velocity or transit time in the composite material may introduce undesirable inaccuracies. It is believed that a reason for the inaccuracies is that the correction does not take into account the fact that a component of the measured value (which is corrected to obtain COrr) is attributable to the matrix portion of the formations, the matrix generally being virtually lossless material. Accordingly, the correction to obtain COrr tends to introduce overall loss correction to a value when at least a portion of the ; value does not actually require loss correction.
Applicant's invented apparatus and method utilizes a different approach wherein it is not attempted to determine the loss-corrected velocity or propagation time through the composite formation.
One aspect of the invention is directed to a method for determining the water-filled porosity of formations surrounding a borehole; said formations including fluid in a matrix, comprising the steps of: deriving a first signal which is a measure of the travel time of micro-wave electromagnetic energy propagating over a pre-determined distance in the formations; deriving a secondsignal which is a measure of the attenuation of said microwave electromagnetic energy propagating over said predetexmined distanse in the formations;

~6~3~

generating a third signal as a function of the second signal which is indicative of losses of microwave electro-magnetic energy attributable to fluid in said formations and useful, when combined with said first signal, for providing an indication of the lossy-water-filled porosity.
Another aspect of the invention is directed to apparatus for determining the water-filled porosity of formations surrounding a borehole, said formations including fluid in a matrix, comprising: means for deriving a first signal which is a measure of the travel time of microwave electromagnetic energy propagating over a predetermined distance in the formations; means for deriving a second signal which is a measure of the attenuation of said microwave electromagnetic energy propagating over said predetermined distance in the formations; means responsive to said second signal for generating a third signal which is indicative of losses of microwave electromagnetic energy attribuatable to fluid in said formations and useful, when combined with said first signal, for providing an indication of the lossy-water-filled porosity.

`" ~Q~ 39 BRIEF DESCRIPTION OF T~E DRAWINGS

.

FIG. 1 is a schematic representatio~, partially in block form, of an apparatus incorporating an embodiment - of the invention.

FIG. 2 illustrates, in simplified form, the S nature of propagation of a microwave electromagnetic surface wave in the formations.

FIG. 3 is a block diagram of the amplitude comparator of FIG. 1.

FIG. 4 is a block diagram of the computing module -10 of FIG. 1.

l~lq~
DESCRIPTION OF T~E PREFERRED EMBOl~ qENT
:

Before proceeding to describe the apparatus and method of the present invention in detail, it is helpful in understanding the invention to review subject matter developed in the abovereferenced U. S. Patent No.
3,944,910 of Rau.
Consider a plane electromagnetic wave propagat-ing in a lossless dielectric medium. The wave propagates with a velocity where ~ is the magnetic permeability and the dielectric constant of the medium. If the type of materials under consideration are generally non-magnetic, ~ can be considered as a known constant and ~ can be determined from the relationship ~ 2 (4) llV
Next, consider two points in fixed space relation along the direction of propagation of the wave. For a given angular frequency, ~, the phase difference of t'ne wave as between the two points is ~ _ vL = ~L (5) where r is the distance separating the two points and ~ is the phase constant of the wave as defined by ~ Z U
Substituting from equation (5) it is seen that ~ can be obtained after determination of the phase constant from the relationship ~ = ~ ~ as = 2 (6 ) ~1 .
B

In terms of the phase measured between two points, the appropriate expression is

2 2 ~ L ~ ~
The preceding relationships are valid for a lossless material, but subsurface media of lnterest generally have appreciablé lossiness. The propagation constant, y, of a plane electromagnetic wave propagating in a lossy medium is a complex quantity of the form y _ ~ ~ ~ (8) where o is the conductivity of the medium. For the case where a is zero or ~ery small, the "loss tangent" term ~a can be neglected and we have y = ~ = ~ ~ which is consistent with equation (6) fox the lossless case. When o is significant, however, the loss tangent term can be kept relatively small by having ~ very large. In such case, the equation (6) is again approximately valid. For example, when ~ = 0.2, a calculation of in accordance with equation (4) yields only about a 0.5 error as compared to the case where ~ is zero.
As discussed in the referenced patent, there are practical limits on how high a frequency can be used in making meaningful measurements in a borehole. Thus, for formations~of appreciable conductivity, the use of the highest practical frequencies may still result in a significant loss tangent which, if ignored, wili lead to exror in measuring the apparent value of dielectric constant.
In the refexenced patent, measurements of dielectric constant are corrected for loss tangent, for example to obtain the .

3:~

quantity ~ referred to in the Background hereof.
To better understand the techniques utilized herein, it is helpful to initially represent the real and imaginary parts of the propagation constant, y, as ~ and -5 a, respectively, so we have Y = ~ + ~ t9) where ~ is associated with wave attenuation or loss.
~Note that the propagation constant is used in the well known wave equation in the form e , so the real part of the propagation constant bacomes the imaginary part of the exponent and vice versa.) Squaring equations (8) and (9) and equating the real parts of each gives - a2 = tlll12 (10) or ~ 2 + a2 (11) The equation ~11) will be referred to hereinbelow.
Applicant has discovered that an improved determination o the water-filled porosity of formations surrounding a borehole can be obtained by utilizing a time average relatio~ship wherein the measured travel tLme of microwave electromagnetic energy through the composite formation is not itself corrected. Instead, the term or signal representing t~e travel time of the energy in the water (only) is corrected. This is done by using a measured value of attenuation to obtain a quantity which represents the travel time of the energy in water having an appropriate degree of lossinesq. In particular, the following relationship is set forth:

1~6~3~

pl ~w~pWl + (1 - ~W) ~
whera ~ is the water-filled porosity of the formations, tpl is the measured transit time of the microwave electro-magnetic energy (uncorrected) ,~is the propagation time through the matrix (generally considered lossless), and t 1 is the propagation time through water having the lossiness of the formation water. In all cases, the propagation times are for a particular travel distance in the formation, determined from the spacing between a pair of receivers. Solving (12) for ~w yields = tpl ~ tPm (13) pwl pm In the present embodiment, the quantity t is pwl determined from the relationship ~ ~ (14) where tpwo is the transit time of the microwave electro-magnetic energy through lossless water, which is determinable, for a particular frequency and temperature, in a laboratory, is the measured attenuation constant~
and R is a constant. The relationship (14) is seen to be of the form of the equation (11~ and indicates that the propagation time through a lossy fluid is greater than the propagation time through an eguivalent lossless fluid by an amount which depends upon a. A relatively accurate determination of ~ is thus obtained by loss-correcting only the quantity which has lossiness fully attributable to it. It should be noted that the equation (12) considers = ~ in setting forth a time average relationship; i.e., w t / l) iQ~ 39 any hydrocarbon content is considered as part of the matrix from the standpoint of propagation tLme determina~ion.
It should be noted that the relationship (14) set forth abo~e, which applicant has found to be effecti~e in practice, deviates from expected theory in one respect, which can be described as follows: Assume that the measured attenuation of the bulk formation is designated as ac and can be expressed as aC ~waW + (1 ~W) am (15) where aw is the attenuation attributable to the formation watex and am is the attenuation attributable to the ormation matrix. Since am is close to zero, equation (15) reduces to w (16) a or a = ~ ~17) which merely express that the bulk formation losses are volumetrically "adjusted" (by a factor of ~w) to take account of the fact that loss ~s essentially only occurring in that fraction of the bulX formation occupied by the water. Using the general form of equation (11) one might expect a suitable ~.
expression for tpWl to be pwl ~wo w (18) where each of the terms is related to the water volume; i.e., t and t , as well as a , are for some given consistent pwl pwo w volume of water through which the energy has propagated.
Since the measured ac is associated with energy which has propagated through the bulX formation (matrix as well as fluid), .: .

/~

3~ 1!

~6~3~

it would stand to reason that equation (17~ could be utilized to obtain a volumetrically consistent ~w Substituting into equation (18) yields the expected relationship:
I Ka ~ 2 c (19) As stated, howe~er, tne relationship (14), which excludes the division of the second term under the radical by the factor ~w ; has been found to yield effecti~e and accurate results.
Referring to FIG. 1, there is shown a representative embodiment of an apparatus in accordance with the present in~ention for investigating subsurface formations 31 traversed by a borehole 32~ The borehole 32 is typically filled with a drilling fluid or mud which contains finely di~ided solids in suspension. The investiga~ing apparatus or logging device 30 is suspended in the borehole 32 on an armored cable 33, the length of which ~ubstantially determines the relative depth of the device 30. The cable length is controlled by suitable means at the surface, such as a drum and winch mechanism (not shown).
The logging device 30 includes an ~longated cylindrical support member 34, the interior portion of which has a fluid-tight housing containing the bulk of the downhole electronics. Mounted on support member 34 are a pair of bowed springs 35 and 36. The spring 35 has mounted thereon a pad member 37 which contains, ~3 .

3~

inter alia, a transmitting antenna T and vertically spaced receiving antennas Rl and R2. Mounted on the spring 36 is a secondary pad member 38 which may be an inactive pad that facilitates smooth vertical movement of the device 30 through the borehole. If desired, however, the pad 38 may contain electrodes or like additional means for investigating the surrounding formations. Electronic signals indicative of the information obtained by the logging device are trans-mitted through the cable 32 to a computing module 85and recorder 95 located at the surface of the earth.
The particular means shown in FIG. 1 for maintaining the antennas in engagement with the borehole wall is illustrative, and it will be appreciated that other suitable means for accomplishing this objective, such as hydraulic means, can be utilized.
FIG. 2 illustrates, in simplified form, the nature of propagation of the electromagnetic wave to be measured with the apparatus of FIG. 1. (For a more detailed description of the wa~e propagation path, reference can be made to the Rau U. S. Patent No.
3,944,~10.) In FIG. 2 the pad 37 is shown positioned ~LQ ~r 6~ ~ 39 against the side of the borehole 32 which, as above-stated, is filled with a drilling mud. Generally, the fluid pressure in the form~tions traversed by- a borehole is less than the hydrostatic pressure of the S col-~n of mud in the borehole, so that the mud and mud fil~rate flows somewhat into the formations. The formations tend to screen the small particles suspended in the mud so that a mudcake is ~ormed on the walls of the borehole. The thickness of the mudcake varies wi~h formation parameters such as permeability, but at least a very thin mudcake is usually present on the borehole wall. In FIG. 2, the pad 37 contacts a mudcake 40 which is shown as being of exaggerated scale thickness for illustrative clarity.
lS The transmitting antenna T emits microwave electromagnetic energy into the formation as represented by the arrow A. A resultant lateral wave propagating --~n the ~ormation is represented by the arrow B and it~ extension, arrow C. m e lateral wave continuously sheds energy back into the more lossy media (the mud-cake), and the portions of energy which are shed at the approximate locations of the receivers Rl and /`! . , ~ .
' /~ ' .

R2 are represented by the arrows D and E, respectively.
If the pathlengths represented by arrows D and E are assumed to be substantially e~ual, it is seen that the difference in pathlength between the energy received at Rl (via path A-B-D) and the energy received at R2 (via pa~h A-B-C-E) is the distance represented by arrow C; i.e., the distance between the receivers. Accordingly, a differential receiver arrange-ment allows investigation of the por~ion of the formation lying approximately opposite the separation between Rl and R2. Typically, but not necessarily, the investigated formation will be the "flushed" or "invaded" zone which surrounds the mudcake in the borehole and contains fluids from the mud which filter through the mudcake.
Referring again to FIG. 1, the downhole electronics contained within the member 34 are shown, for convenience of illustration, at the side of the borehole. A solid state oscillator 45 provides output energy in the microwave region of the spectrum. The microwave region is defined herein as including the range of frequencies between about 300 MHz and 300 G~z. The oscillator 45 may operate at the suitable frequency of 1.1 G~z; i.e., 1.1 x 10 cycles per second.
The output of oscillator 45 is coupled through an isolator 46 to the transmitting antenna T. Microwave energy is trans-mitted into the surrounding formations and propagates through the formation~ in the manner previously described. The energy which arrives at the receiving antennas Rl and R2 i5 respectively coupled to input terminals of the mixers 47 and 48. The signals which arrive from Rl and R2 are out of phase with each other by an amount which depends on the phase . ~

, ~ .

~ 39 constant 3 and have an amplitude ratio which depends upon the attenuation constant a. Secondary input terminals of the mixers are supplied with microwave energy at a frequency that is separated from the transmitter frequency by some relatively low frequency which is typically in the radio frequency range. In the embodiment shown, a solid state oscillator 49 supplies microwave energy to mixers 47 and 48 at a frequency of 1.1001 GHz, or 100 KHz ~ove the trans-mitter frequency. The output signals 47A and 48A of the mixers 47 and 48 therefore contain the difference frequency of 100 RC. In accordance with well ~nown principles, the signals 47A and 48A maintain the phase and amplitude relationships of the signals from Rl and R2, but the task of phase detection is greatly facilitated at the lower frequency of the mixed signals. To insure that the difference frequency between the outputs of the oscillators 45 and 49 remains at 100 KHz, the oscillator outputs are sampled and fed to a mixer 50. The output of the mixer is received by a frequency stabilization circuit 51 which detects drifts from the 100 RC standard and generates a correction signal 51A which controls oscillator 49 ln the manner of a conventional "phase-locked loop."
The signals 47A and 48 are applied to a_phase detector circuit 53 and to an amplitude comparator 54. The output of the phase detector 53 is a signal level which is proportional to the phase difference ~ betweèn the signals received at Rl and R~, and thus proportional to ~ in accordance with ~ = ~/L, where L is the distance separating the two receivers. For a particular frequency of operation, ~, the phase difference ~ is also proportional to the tra~el ^

l~q~6~39 t~me through the formations over the distance L in accordance with tpl = L~' where tpl is the travel time.
The output of a~plitude comparator 54 is a signal level which is proportional to the attenuation constant ~. A
convenient circuit 54 for obtaining an output signal pro-portional to a is shown in FIG. 3. The signals 47A and 48A are respectively applied to the logarithmic amplifiers 55 and 56 whose outputs are fed to the difference amplifier 57~ The ou~put of the difference amplifier 5? is a signal level propor~ional to ~. This can be visualized by representing the amplitude of the wave energy received at R as Ae az, where A is an amplitude constant and z is the distance separating T and Rl. It follows that the amplitude of the wave energy received at R2 i5 Ae ( ), where L is the distance separating the receivers Rl and R2. The ratio of the wave amplitudes at the two receivers is thereore - a(z+L) Ae _ -aL
~ -~az ~ e The log of the ratio of the wave amplitudes is therefore proportional to a. It will be appreciated that the circuit 54 of FIG. 3 accomplishes the same mathematical result by taking the difference of the logs of the wave amplitudes.
The output representative of tpl and a are transmitted to the surface over the conductor pair 53a and 54a which in actuality pass through the armored cable 33.
Typically, these signals may be DC levels which are stepped-up by amplification before transmission to the surface. At the surface of the earth the signals on the lines 53a and 54a are applied to a computing module 85 which computes the water-filled J~

l~.S~3~

porosity, ~w~ of the formations in accordance with the relationships (13) and (14). The computed porosity is recorded by a recorder 95 that is conventionally driven as a function of borehole depth by mechanical coupling to S a rotating wheel 96. The wheel 96 is coupled to the cable 33 and rotates in synchronism therewith so as to move as a function of borehole depth. ~hus, the water-filled porosity is recorded as a function of borehole depth by the recorder 95.
FIG. 4 is a block diagram of the computing module 85 which receives the signals on lines 53a and 54a that are indicative of the measured values of tpl and a, respectively. The signal representative of a is coupled to a squaring circuit 91 and the signal representative of t i5 coupled to the positive input terminal of a pl , difference amplifier 95. The output of squaring circuit 91, which is representative of a , is coupled to one input of a multiplier circuit 92, the other input to which is an adjustable signal representative of a constant, K. The output of circuit 92, which is a signal having a value representative of ~a2 is one input to a summiny circuit 93.
The other input to summing circuit 93 is a signal level representative of the value of tpwo. As previously stated, t is the travel time of microwave electromagnetic energy, pwo at the frequency ~, in substantially lossless water over the distance ~. This value is determinable in a laboratory as, for example, 28.3 nanoseconds per meter at 1.1 GHz a~
lQ0 F. The output of summing circuit 93 is a signal representative of tpwo + R~ , and this signal is coupled /~
_~ -3~

to a square root circuit 94, the output of which is a signal representative of tpWl, in accordance with the relationship (14).
The output of square root circuit 94 is coupled to the positive input terminal of a difference amplifier 96. The negative input terminal of each of the difference amplifiers 95 and 96 receives a signal representative of tpm; i.e., representative of the transit time of the microwave electromagnetic energy of frequency ~, over the distance L, for the particular matrix.
It is assumed that lithology is known from other logging or coring information. An example of a typical value for tpm for the frequency of interest is 8 nanoseconds per meter. The out-puts of difference amplifiers 95 and 96 are coupled to a ratio circuit 97 which produces an output representative of ~w in accordance with the relationship (13).
The invention has been described with reference to a particular embodiment, but variations within the spirit and scope of the invention will occur to those skilled in the art.
For example, while circuitry has been described for generating analog signals representative of the desired quantities, it will be understood that a general purpose digital computer could readily be programmed to implement the techniques as set forth herein. Also, it should be noted that the advantageous princi-ples of known borehole compensation techniques and/or of redun-dant processing channels, such as are disclosed in U.S. Patent 3,849,721 of Thomas Calvert issued November 19, 1974, can be utilized, if desired, in conjunction with the present invention.
Further, it will be understood that the measured values can, if desired, be corrected for mudcake effect, spreading, or tempera-ture variations, as is known in the art.

~ 39 Finally, although the illustrative embodiment shows - ~arious quantities as being derived directly from a logging device, these quantities may alternatively be derived from storage media or communicated from a logging location~

~d .
.~ .

Claims (14)

60.453

1. A method for determining the water-filled porosity of formations surrounding a borehole, said formations including fluid in a matrix, comprising the steps of:

deriving a first signal which is a measure of the travel time of microwave electromagnetic energy propagating over a predetermined distance in the formations;
deriving a second signal which is a measure of the attenuation of said microwave electromagnetic energy propagating over said predetermined distance in the formations;
and characterized by:
generating a third signal as a function of the second signal which is indicative of losses of microwave electromagnetic energy attributable to fluid in said formations and useful, when combined with said first signal, for providing an indication of the lossy-water-filled porosity.

2. The method as defined by Claim 1, characterized in that the third signal is representative of the travel time of microwave electromagnetic energy propagating over said predetermined distance in water having a lossiness determined from said second signal.

3. The method as defined by Claim 2, characterized by the step of generating a porosity-indicative signal as a function of said first and third signals.

4. The method as defined by Claim 3 further characterized by the step of deriving a fourth signal representative of the travel time of microwave energy propa-gating over said predetermined distance in said matrix, with said step of generating a porosity-indicative signal comprising generating said porosity-indicative signal as a function of said first, third and fourth signals.

5. The method as defined by Claim 4 characterized in that the step of generating a porosity-indicative signal comprises generating said porosity-indicative signal as a ratio of the difference between said first and fourth signals and the difference between said third and fourth signals.

6. The method as defined by Claim 1 characterized in that said third signal is also a function of a fifth signal, said fifth signal being representative of the travel time of microwave electromagnetic energy propagating over said predetermined distance in substantially lossless water.

7. The method as defined by Claim 6 characterized in that said third signal is generated as a function of the square root of a quantity which is the square of said fifth signal plus the square of said second signal.

8. Apparatus for determining the water-filled porosity of formations surrounding a borehole, said formations including fluid in a matrix, comprising:

means for deriving a first signal which is a measure of the travel time of microwave electromagnetic energy propagating over a predetermined distance in the formations;
means for deriving a second signal which is a measure of the attenuation of said microwave electromagnetic energy propagating over said predetermined distance in the formations;
and characterized by:

means responsive to said second signal for generating a third signal which is indicative of losses of microwave electro-magnetic energy attributable to fluid in said formations and useful, when combined with said first signal, for providing an indication of the lossy-water-filled porosity.

9. Apparatus as defined by Claim 8, characterized in that the third signal is representative of the travel time of microwave electromagnetic energy propagating over said predetermined distance in water having a lossiness determined from said second signal.

10. Apparatus as defined by Claim 9, characterized by means responsive to said first and third signals for generating a porosity-indicative signal as a function of said first and third signals.

11. Apparatus as defined by Claim 10 further characterized by means for deriving a fourth signal representa-tive of the travel time of microwave energy propagating over said predetermined distance in said matrix, and wherein said means for generating a porosity-indicative signal is also responsive to said fourth signal for generating said porosity-indicative signal as a function of said first, third and fourth signals.

12. Apparatus as defined by Claim 11 characterized in that means for generating a porosity-indicative signal generates said porosity-indicative signal as a ratio of the difference between said first and fourth signals and the difference between said third and fourth signals.

13. Apparatus as defined by Claim 8 characterized in that said means for generating said third signal is also responsive to a fifth signal, said fifth signal being representative of the travel time of microwave electromagnetic energy propagating over said predetermined distance in substantially lossless water.

14. Apparatus as defined by Claim 13 characterized in that said means for generating said third signal is opera-tive to generate said third signal as a function of the square root of a quanity which is the square of said fifth signal plus the square of said second signal.