IE48252B1 - Apparatus and method for determining velocity of acoustic waves in earth formations - Google Patents

Abstract

To determine the velocity of acoustic wave energy propagating in formations surrounding a borehole (14), a transmitter (16) periodically transmits acoustic wave energy into the surrounding formations. This energy is received at a plurality of spaced receivers (18.1-18.4). A first portion of the signal from the receiver 18.1 is then correlated with respective first portions of the signals from the other receivers (18.2-18.4) the time reference for the correlation of each of the first portions of the signals from the other receivers (18.2-18.4) being selected, relative to the time reference of the signal from the receiver 18.1, in accordance with the spacing of the receivers (18.1-18.4) and a first estimated velocity. The correlation is repeated for four more estimated velocities, and a velocity value which yields the best correlation is selected by interpolation as a first provisional velocity. The correlation for several assumed velocities and interpolation is then performed on second portions of the signals from the receivers (18.1-18.4) to obtain a second provisional velocity. An output velocity is then determined as a function of the first and second (and possibly other similarly-generated) provisional velocities. The velocity of the compressional component and/or the shear component of the acoustic wave energy can be obtained.

Description

This invention relates to determination oi properties oi earth formations surrounding a borehole and, more particularly, to an apparatus and method ior determining the velocity oi propagation oi acoustic wave energy propagating in such formations. Acoustic well logging is widely used to provide information concerning the characteristics of earth formations. Generally, measurements are made of the velocities of acoustic waves to reveal properties, such as porosity, of the formations surrounding a borehole. An acoustic well logging tool for measuring the velocity of acoustic waves typically employs a sonic pulse transmitter and a plurality of sonic receivers selectively spaced from the transmitter. The sonic receivers include transducers to convert the incident acoustic wave to an electrical waveform and suitable amplifiers to transmit the waveforms to surface located processing equipment.

It is well known that the compressional components of the acoustic wave energy generally travel at a higher velocity than the shear components of the acoustic wave energy It is desirable, at each depth level,to obtain a reading oi both the compressional velocity and shear velocity oi the acoustic wave energy propagating through the formations since both contribute useful information concerning the formations. When a plurality oi spaced receivers are employed, a correlation technique can be used to determine the desired velocity by correlating the signals obtained at the different receivers to obtain an output velocity value which optimally accounts for the difference in arrival times of the signals at the different receiver locations. Briefly, most correlation techniques are based on the assumption that substantially the same signal waveform will arrive at the different receiver locations at different times, the delay as between successively further receiver locations depending upon the distance between receivers and the velocity of propagation of the wave energy in the formations as between the receiver locations. Thus, a correlation technique can be employed to compare a delayed version of a signal arriving at a closer receiver location with a signal arriving at a more remote receiver location, the amount of delay (for a given known distance as between the receivers) depending on an assumed trial value of velocity of the acoustic wave energy as between the receiver locations in question.

In simplified terms, the correlation technique involves trying various trial velocity values and determining whioh one provides the best match.

To visualize a correlation procedure, assume that two signals to be correlated are plotted one-beneath-another on a common time axis, with the lower plotted signal being delayed (advanced in a direction cf increasing time on the time axis) with respect to the upper plotted signal.

Assume further that the time differential as between the signals (resulting from signal propagation time through formations of unknown nature over a known distance) is unknown. In performing the correlation, the signals are compared using a selected slope to account for the relative 8 2 5 2 time difference between the signals. The selected slope is representative of a trial velocity for the particular correlation since, as previously noted, the time delay as between the two signals (and therefore the slope as between corresponding points on the two signals plotted on the same time axis) is a function of the velocity of propagation of the signals. The comparison is generally done by multiplying the instantaneous values of the two signals at each corresponding point thereof and summing all the products to obtain a single correlation figure. Further correlation f-cures may then be determined for different trial velocities (or slopes), and the correlation figures may be plotted versus velocity. Ideally, this plot,, called a correlogram will have a single, well defined peak, the peak (i.e., the single highest correlation figure) indicating the true velocity of propagation of the wave energy. If the wave energy has components which travel at different velocities, multiple peaks may be encountered. Noise and other complications in hostile borehole environment will, of course, also introduce difficulties in interpreting correlograms.

In U.K. Patent Specification Mos. 1,534,004, 1,534,005 and 1,534,006 there is disclosed an acoustic logging technique wherein a multiple-fold correlation is utilized; that is, the instantaneous values of three or more signals are multiplied by each other in obtaining each correlation figure. This and other techniques set forth in that application have led to substantial advances in the art, but it is an object of the present invention to provide still further improvements in obtaining accurate values for the velocity cf the acoustic wave components of sonic waves propagating in formations surrounding a borehole.

The main aspects of the present invention are concerned with techniques for obtaining velocityrepresentative values from the plurality of received signals, and not with any particular technique for obtaining these signals in the first place (the method of obtaining these signals being known in the art and disclosed, for example, in the above-referenced UK Patent Specifications). Accordingly, it will be understood that the principles of the invention to be described are applicable to the processing of signals obtained directly from a borehole or to derived signals which were previously recorded or communicated from a remote location. Also, as used herein, the term signals is intended generically to include any waveforms that are in a form suitable for machine processing.

According to one aspect of this invention there is provided a method for determining the velocity of acoustic wave energy propagating in formations surrounding a borehole, comprising the steps of: placing in the borehole an acoustic logging tool having a source of acoustic wave energy and a plurality of acoustic wave energy receivers disposed at respective locations spaced from said source; causing said source to emit acoustic wave energy into the formations surrounding the borehole; deriving a plurality of signals respectively representative of acoustic wave energy received by said plurality of receivers at said locations; correlating a first portion in time of the signal from one of said locations with first portions in time of signals from the other locations, the first portions of the signals from the other locations being determined from an assumed velocity and the correlation being performed for a number of different assumed velocities to obtain a resultant first provisional velocity; correlating a second portion in time of the signal from said one of said locations with second portions in time of signals from the other locations, the second portions of the signals from the other locations being determined from an assumed velocity and the correlation of second portions being performed for a number of different assumed velocities to obtain a resultant second provisional velocity; and generating an output velocity as a function of said first and second provisional velocities.

According to another aspect of the invention there is provided an apparatus for determining the velocity of acoustic wave energy propagating in formations surrounding a borehole, comprising: means arranged to derive a plurality of signals respectively representative of acoustic wave energy received at a plurality of locations in the borehole spaced from a source of acoustic wave energy; means arranged to correlate a first portion in time of the signal from one of said locations with first portions in time of signals from the other locations, the first portions of the signals from the 25 other locations being determined from an assumed velocity and the correlation being performed for a number of different assumed velocities to obtain a resultant first provisional velocity; means arranged to correlate a second portion in time of the signal from said one of said locations with second portions in time of signals from the other .locations, the second portions of the signals from the other locations being determined from an assumed velocity and the correlation of second portions being performed for a number of different assumed velocities to obtain a resultant second provisional velocity; and means arranged to generate and output velocity as a function of said first and second provisional velocities.

According to a third aspect of the invention there is provided an apparatus for determining the velocities of different modes of propagation of acoustic wave energy propagating in formations surrounding a borehole, comprising: means arranged to derive a plurality of signals respectively representative of the wave energy received at a plurality of spaced locations in the borehole; means arranged to identify the location and duration of a first characteristic portion in time of the signal from one of said locations; means arranged to correlate said first portion of the signal from said one of said locations with first portions in time of signals from the other locations, the first portions of the signals from the other locations being determined from an assumed velocity and the correlation being performed for a number of different assumed velocities to obtain a resultant first provisional velocity; means arranged to identify subsequent characteristic portions in time of the signal from said one of said locations; means arranged to correlate said subsequent portions of the signal from said one of said locations with subsequent portions in time of signals from the other locations, the subsequent portions of the signals from the other locations being determined from an assumed velocity and the correlation of each of said subsequent portions being performed for a number of different assumed velocities to obtain resultant further provisional velocities; and means arranged to generate 48352 output velocities of said different modes of propagation as a function of said first provisional velocity and said subsequent provisional velocities.

The first and second portions of the signals may be compressional wave components thereof, and the output velocity will then be an output compressional wave velocity. Alternatively, the first and second portions of the signals may be shear wave components thereof, in which case the output velocity will be an output shear wave velocity. Of course, the technique may be utilized on both the compressional wave components and shear wave components of the signals to obtain both an output compressional wave velocity and an output shear wave velocity. The first and second portions preferably each have a duration of an integral number of half-cycles, not exceeding two full cycles, of the signal from said one location.

Generally, the output velocity (compressional and/or shear) is determined for each of a number of adjacent depth levels. In a disclosed embodiment, the generated output velocity, in addition to being a function of the previously mentioned provisional velocities, is selected taking into account the output velocity that had been generated at an adjacent depth level.

A preferred embodiment of the present invention, which involves determining provisional wave velocities by correlating, over a limited range of velocities (e.g. a range of about fifty microseconds per foot - 164 microseconds per metre), relatively small portions of the different receiver waveforms, tends to reduce prior art problems associated with multiple correlation peaks that could occur when longer portions of the signals are correlated. Also, individual portions of the signals containing spurious information have less tendency to upset the selection of accurate output values if, as in an embodiment to be described, selection is from a plurality of provisional values so that an individual spurious value is likely to be discarded without affecting the output values.

Methods and apparatus in accordance with this invention will now be described, by way of example, with reference to the accompanying drawings, in which:FIG. 1 is a block diagram of the apparatus, and includes a recording or log of acoustic velocity versus depth that can be set forth using the invention; FIG. 2 is a graph which illustrates the nature of the signals or waveforms received at the four receivers of the logging device in FIG. 1; FIG. 3 is another recording or log of acoustic velocity versus depth that can be set forth using the invention; FIG.s 4A, 4B, 4C and 4D, when placed one-belowanother, illustrate a flow chart suitable for programming the processor 40 of FIG. 1; FIG. 5 illustrates the manner in which an initial 8 2 5 2 estimate of slope is obtained from a pair of waveforms; FIG. 6 illustrates the manner in which the correlation window is initially positioned with respect to the waveform; FIG. 7 is a flow chart representative functions of the block 120 of FIG. 4 and sets technique of fourfold correlation; of the forth a FIG.s 8A and 8B illustrate the manner in which a fourfold correlation is performed, with the correlations of FIG.s 8A and 8B each utilizing different slopes; FIG. 9 illustrates a technique for obtaining the slope or velocity which yields a maximum correlation; and FIG. 10 is a plot of provisional compressional 15 velocity values and provisional shear velocity values at a particular depth level.

Referring to FIG. 1, a sonic borehole logging tool 10, with centering elements not shown, is illustrated as being suspended from a cable 12 in a borehole 14. The tool 10 includes a transmitter 16 located at the bottom of the tool and preferably four or more selectively spaced sonic receivers 18.1, 18,2, 18.3 and 18.4. The cable 12 is shown fed over a pulley 20 from a cable supply 22 and is electrically connected to a panel 24. The panel 24 includes suitable amplifiers, switching circuits and electrical supplies for tool 10 and the firing of sonic transmitter 16. The tool 10 includes suitable amplifiers and controls needed to enable receivers 18 to sequentially detect sonic waves and provide panel 24 with waveforms representative of the acoustic waves incident upon transducers in receivers 18.

In the operation of the sonic logging tool 10, the transmitter 16 is regularly energized (about ten times per second) to produce sonic pulses such as 26 in FIG. 2 while the tool is moved upwardly at a speed of the order of cne foot per second (.305 metres per second). Each pulse 26 is directed at the formation in which an acoustic wave such as 28 is launched.

The acoustic wave 28 has components which travel through the formation towards the receivers 18 and in turn produce wave components 29 which are incident upon the receivers 13 to produce a set of signals or waveforms such as 30.1, .2, 30.3 and 30.4 shown in FIG. 2. Waveforms 30 may also 48253 include those sonic waves which, after travel through the borehole, are incident upon receivers 18.

The transmitter is periodically activated, and the receivers 18 preferably are alternately enabled to generate waveforms in the sequence as shown by waveforms .1-30.4. In this manner, waveforms which are of several milliseconds duration can be analyzed to detect velocity in a manner to be described. In addition, the sequential enabling of receivers 18 allows their outputs to be IQ multiplexed onto a common line towards panel 24. This is an advantage when the analog form of the outputs of receivers 18 is to be conducted to panel 24 since the analog waveforms 30 are all exposed to the same electronic amplification and cable attenuation to enable amplitude comparisons. The advance of the tool 10 during the time needed to generate one complete set of waveforms 30 can be controlled and is not sufficient to appreciably affect velocity measurements.

The receivers 13 each typically but not necessarily bear a known spacing d from each other (i ,e. d^ = d^ = d^ = d in FIG. 1) while the distance ”S” between receiver 18.1 and transmitter 16 is known to further accurately locate all the receivers from transmitter 16. Thus, for example, if the arrival of the compressional wave could be precisely detected at each receiver, then the time' between respective receiver arrivals would provide an accurate determination of the velocity of the compressional wave. As described in U.K. Patent Specification Nos. 1,534,004, 1,534,005 and 1,534,006, employing a multiple number of receivers 18, such as four, allows a generally more accurate and reliable determination of wave velocities , although increased distance of the far receiver from the transmitter results in a substantially attenuated waveform at this receiver.

The signals 30 may arrive at the panel 24 either in analog form or sampled form or may be digitized downhole. For example, the tool 10 may be provided with an analog to digital (A/D) converter (not shown) to sample waveforms 30. The samples would then be transmitted along cable 12 to panel 24 for further processing. The waveforms 30 may be derived from a previously obtained record, communicated from another location, or derived directly from receivers 18. In the system depicted in FIG. 1, the waveforms 30 arrive in analog form at panel 24 where they may be coupled to a cathode ray tube 32 for real-time display and to a magnetic tape recorder 34 to provide an analog record. In order to utilize a sampled form of waveforms 30 they are shown coupled along a line 35 to an analog to digital converter 36 which generates samples of the waveforms at a high sampling rate on an output line 38.

Although not specifically illustrated in FIG. 1, panel 24 also provides depth signals from a depth sensor operatively coupled to the cable 12. The depth signals preferably are presented on line 37 in a format which is compatible with the format employed for the samples on line 483^3 38. The generation of depth signals is well known in the art of well logging and shall not be further described.

The output from waveform sampler 36 is applied to an apparatus 40 which processes the samples to produce velocity—representative values from waveforms 30. The apparatus 40 preferably is a general purpose digital computer, but may be any general purpose or special purpose digital or analog processor or special purpose circuitry. Also, with a recording capability present, processing need not necessarily be performed in real time.

T. -.e output of apparatus 40 may be, for example, a plot 42 of the velocity of the compressional and shear wave components of the sonic wave as a function of well depth in feet. Curve 43 shows the velocity of the compressional wave, designated ΔΤ , and curve 44 shows the velocity of c the shear wave, designated ATs. The plot 42 is made by moving suitable logging paper in the direction indicated by arrow 45 while the tool is raised as reflected by the indicated depth values. (It is conventional in acoustic logging to present velocity in terms of microseconds per foot, which is actually the inverse of true velocity. This convention vri.ll be followed herein. 1 ps per foot equals 3.28 ps per metre.) As described in the abovereferenced copending U. S. application Serial No. 581,381, and as illustrated in FIG. 2, the receivers 18 are sequentially enabled for successive firings of the transmitter 16. A portion of the electronics in device 10 provides waveform identification signals, which may be coded pulses such as a two bit code word, and which determine the particular receiver that is enabled. Such receiver waveform identification code, together with a pulse to identify the firing time of the transmitter 16 are delivered along cable 12 to panel 24.

A magnetic recorder 47 is provided to record the waveform samples produced by A/D converter 36 as well as the waveform identification code transmitted along cable 12.

As described in the referenced U.K. Patent Specifications, the panel 24 may be provided with a conversion control network 46 which selectively delays initiation of A/D converter 36, depending upon which receiver waveform is to be sampled. For example, if receiver waveforms 30 are each sampled by starting the same amount of time after the firing of transmitter 16 (or upon occurrence of its output pulse 26), the amount of useful information obtained varies as a result of receiver spacings d. In order to preserve greater portions of the waveforms for subsequent analysis, a conversion control network can be used, but will not be treated herein.

A/D converter 36 commences conversion at the end of a delay Δτ which is a function of ths distance between receiver 18 and transmitter 16 and the fastest expected velocity of the acoustic wave. The conversion process continues for a sufficient time period to provide samples of the waveforms of interest. A sampling rats of ten microseconds, for example, may be employed, and a total of 512 samples, or about five milliseconds, of sampled waveform duration will be obtained and be sufficient in most cases to determine the compressional and shear velocities of the sonic waves, as will be described hereinafter. However, it will be understood that the number of samples used in any one group may be varied with 512 being generally used herein as an illustrative example. A/D converter 36 is provided with a suitable counter (not shown) which terminates the conversion process when the desired number of samples 60 (see FIG. 2) have been generated. The A/D conversion of each waveform 30 results in the generation of a group identified at 62 in FIG. 2 of 512 samples.

Since each sample 60 occurs at a known sampling Xq rate, any one sample has an index value which can be directly related to the time interval measured from the time of occurrence of the sonic pulse 26 which resulted in the sampled waveform. Thus the first sample 60.1 in group 62.1 occurs at a time equal to ώτ + (NxSR) , where N is the index position value (N = 0 for the first sample) and SR is the sampling rate in microseconds. In a similar manner, each sample 60 in the other groups 62.2, 62.3 and 62.4 can be precisely related in time to their associated transmitted sonic pulses 26.2, 26.3 and 26.4. 2q An example of how these index values can be used is as follows: Assume a given reference (e.g. first motion) on the compressional wave is detected for the sample in group 62.1 at an index value NI, and the same given reference is detected in group 62.2 at index value N2. The time IT for that wave to travel the distance between receivers 13.1 and 18.2 would then be equal to (N2 - NI) SR.

At a ten microsecond sampling rate, the velocity of the compressional wave in microseconds per foot would be determined as (N2 - Nl) X 10/d.

The samples applied to processor 40 are assembled in groups of 512 samples with each group coded to identify a waveform from a particular receiver. The processor 40 is provided with a buffer to enable accumulation of a pair of sets of waveforms, wherein each set represents waveforms from all four receivers during a full operational cycle of the transmitter-receiver. Access to the buffer is under control of A/D converter 36 to enable transfer of t.-.e samples as they are produced.

The description set forth in conjunction with FIG.s 1-3, and the specific technique of sampling and using index numbers as described in U.S. Patent Application Serial No. 581,381 are intended to illustrate a manner in which a plurality of signals or waveforms can be obtained from a plurality of spaced acoustic receivers and loaded into a processor at known sampling times. It will become understood, however, that the novel aspects of this invention deal with apparatus and methods for obtaining useful information, such as velocityrepresentative signals or graphs, from the plurality of signals or waveforms 30 which are loaded into processor 40, and any suitable technique can be employed to obtain the signals or waveforms 30 (the terms signals and waveforms being utilized interchangeably herein in this context). For example, all four receiver signals may be digitized, after a single transmitter firing, by a downhole A/D converter.

A buffer memory and data transmission system may also be provided downhole. Also, the described indexing 8 2 5 2 system sets forth an example of how timing registration between the sampled waveforms is established; but it will be understood that alternate methods may be employed.

In this respect, and for clarity of explanation, the index number system set forth above will not be referred to each time timing registration considerations are .described below, although it will be understood that this index number technique can be utilized in each instance to keep track of the timing registration as between samples from the different receiver waveforms.

Referring to FIG. 4, there is shown a flow diagram suitable for programming the processor 40 to obtain output values of ΔΤ= and ATg at each depth level; e.g. the values shown in the graphs 43 and 44 of FIG. 1. The block 105 represents the initiation of processing for the next depth level. Block 110 is then entered, this block representing the determination of the first motion of each of the four receiver waveforms. As is known in the art, first motion can be determined by detecting when each waveform exceeds a predetermined threshold level.

Block 110 represents the determination of an initial estimate of compressional velocity or slope, designated ΔΤ , obtained ce by averaging the slopes as between the first motion of the close receiver (R^) waveform and the first motion of each of the other waveforms. (The term slope, as used in this context, represents a time difference as between waveform sample points received at different receivers that are a known fixed distance apart. Since this characteristic time difference is proportional to the inverse cf velocity, slope determinations are equivalent to determinations of velocity. Accordingly, the terms slope and velocity are, in this context, used interchangeably herein.) The signal received at the closest receiver, R^, is generally the strongest and cleanest signal and it is therefore used as a main reference for obtaining the slopes from which the initial estimate of compressional velocity is determined. In particular, ΔΤ is computed as: ce (fm, - fm.)/d + (fm, - fm.)/2d + (fia, - fm.)/3d at = _2_i_3_1_4_i_ (1) where fm, through fm are the times of the first motions 1 4 respectively detected in the waveforms of R1 through R^, respectively. The three numerator terms are seen to be the three slopes which are averaged to obtain iTcg. FIG. 5 illustrates, in simplified form, how one of the slopes, in particular (fm2 - fm^/d, is obtained.

Returning to FIG. 4, block 115 is next entered, this block representing the establishment of a correlation window whose end points are selected as respectively being at the zero crossing which precedes the firsn motion fm^ and at the second zero crossing thereafter; i.e. a window substantially corresponding to the first full cycle of the compressional wave received at R . This may be readily done by scanning the samples of the R^ waveform from point fm^ until the appropriate zero crossings are obtained.

FIG. S illustrates the correlation window on the R^ waveform.

Two indices, r and m {to be utilized later), are next set equal to unity (block 116). A fourfold correlation is then performed, as represented by the block 120, to determine a provisional oompressional velocity or slope, designated AT(The basic concept of a fourfold correlation is described in the abovereferenced 0. S. patent application Serial No. 581,381.) The present technique of fourfold correlation is set forth in FIG. 7 which is a flow chart representative of the functions of block 120 of FIG. 4, and can be referred to during the present description of the fourfold correlation. An initial assumed slope is first selected (block 701) as being the estimated compressional wave velocity or slope, AT^. The value of each sample point of the Rj_ signal within the correlation window is then multiplied by the value of a sample point on each of the other waveforms, the location of the sample point on the other waveforms being determined using the assumed slope (block 702 and loop 703). For example, the value at the leftmost point of the R^ waveform within the correlation window is multiplied by the value of a point on the R2 waveform that is displaced in time from the point on the Rt waveform by AT (d); i.e. the time that it would have taken ce the wave energy to travel the distance d between receivers R, and R_ (FIG. 1) at an (inverse) velocity AT (with άΣ = d-, = d3 = d for assumed equal spacings). The resultant partial product is then again multiplied by the value of a point on the Rj waveform that is displaced in time frcm rhe point on the R^ waveform by A+ce(2d) and then again multiplied by the value of a point on the R, waveform that is displaced in time from the point on the Rj waveform by λΤ (3d). The resultant product, designated PR, ce is multiplied by weighting function which is utilized to take account of the signs of the values multiplied together. A convenient weighting function is simply unity when the four multiplicands have the same sign and is otherwise zero. The use of the weighting function avoids making a positive contribution to the correlation when two of the signals are positive and two are negative; i.e., the worst possible mismatch. The determination of the product PR is represented by block 704 of FIG. 7. In the simplified illustration of FIG. 8A, this first product is PR^. The same procedure is then performed with respect to the next point of the R^ waveform within the correlation window (diamond 705 and block 706 of FIG. 7) to obtain a product designated as PR2< also shown in the simplified FIG. 8A. This procedure is then repeated for each point of the correlation window to obtain products PR,, PR,....PR , where n 5 3 4 n is the number of points in the correlation window. The summation, designated as Z, of all the products is then determined (block 707 of FIG. 7) as; 112 n where the subscript 1 designates that this summation is a correlation figure associated with the first estimated slope. FIG. 8A is again referred to for an understanding of how Σ is generated.

A different estimated slope is next assumed (diamond 708 and block 709 of FIG. 7), and the procedure for obtaining Σ (in this case designated ΐ ) repeated for Ene new assumed sloae. FIG. 83 illustrates the manner in which Σ * 2 is obtained using a slope which is different than that of FIG. 7A. The number of different slopes which are tried is a matter of choice and, in Ene present embodiment, five different sictes are used co obtain five values cf -, 252 e.g. £ through The five slopes utilized are: ΔΤ (the originally assumed slope, rounded off to the ce nearest multiple of the sampling rate); ΔΤ +10us; c© ATcg-l0ps; ATce+20ys; and ATce-20ps, the ten microsecond multiples being multiples of the sampling rate. The optimum slope can then be determined as a slope at which Σ is a maximum. In the present embodiment, an interpolation technique is preferably utilized to obtain the optimum slope, designated as ΔΤ^ (block 710 of FIG. 7).

This may be done by fitting a second order polynomial through the largest value of I and the two adjacent values of I. An example of the technique is illustrated in FIG. wherein the maximum of the polynomial curve, selected as ΔΤ^, is determined, in a given case, to be between ΔΤ_ -logs and ΔΤ , as shown by the dashed line. In the cs es event that the maximum value of I is found to be at one of the assumed slope extremes (viz., at either ΔΤ -20ps or cs ΔΤ +20vs) , diamond 711 of FIG. 7 causes block 712 to be C8 entered via diamond 712. The slope is then incremented toward the direction of the extreme and the same procedure is followed. Eventually, either an appropriate maximum is obrained (called ΔΤ^ in this case) or a no answer to decision diamond 712 indicates that the attempt at obtaining ΔΤ' has been unsuccessful, c Referring again to FIG. 4, decision diamond 130 is next entered and, if a value of ATj, has been successfully obtained (routine of block 120) , diamond 130 causes entry to another decision diamond 135. If, however, a value of ΔΤ1 was not successfully obtained the block 160 is entered, c In accordance with the criterion of diamond 135, ΔΤ^ is tested as against the value of ΔΤ^ obtained at the previous depth level (except, of course, where this is the first depth level being processed whereupon suitable provision is made for bypassing diamond 135), to determine / y if the value of ΔΤ^ at the present depth level is within a particular range, for example 10us, of the value of ΔΤ), which had been obtained at the previous depth level. If the test result is affirmative, the block 140 is entered, this block representing the acceptance of the value ΔΤ^, IP the accepted value being an output value designated as ΔΤ^ (e.g. curve 43 of FIG. 1). If, however, the answer to the inquiry of diamond 135 is negative, diamond 145 is entered and inquiry is made as to whether ΔΤ' has been rejected c (in accordance with the criterion of diamond 135) for five depth levels in a row. If not, block 160 is entered, this block representing the adoption of a ΔΤ^ for this depth level which is the same as the one from the previous depth level. The block 160, it is seen, is also directly entered from the negative output line of decision diamond 130; i.e.

Ene situation where the attemoted determination of ΔΤ' was c unsuccessful.

Having established an output value of ATC for this particular depth level, an initial estimate is next made of the shear wave velocity or slope, this estimated value being designated as ΔΤ » In the present embodiment, the initial se estimated slope or velocity of. the shear wave is computed as : ΔΤ = 1.35 ΔΤ +3 (sampling rate) se c as indicated by block 150 of FIG. 4. Thus, if, as above, the sampling rate is selected as 10 usee the estimated velocity or slope of the shear wave would be: ΔΤ = 1.35 ΔΤ +30 psec/ft. se c Block 156 is then entered and a lower limit of the^'shear wave velocity or slope, designated ΔΤ , is computed.

In the present embodiment, this lower limit, to be used in a manner described below in rejecting any values of shear wave velocity below the lower limit, is selected as beinc 1.5 ΔΤ , as indicated by the block 156. c The next portion of the flow chart is, in the present embodiment, concerned with movement along the time scale of the R, waveform until such time when useful information concerning the shear wave velocity or slope can be expected. The portion of the R^ waveform used in the previously described correlation is stripped away (block 165); i.e., the time reference is moved forward to the rightmost end point of the previously utilized correlation window. From this point in time, a first motion detection is again performed (as described above) , as represented by the block 170. The nearest zero crossing is then detected, and the beginning of the correlation window is moved to the time of zero crossing detection, as represented by the block 175. The correlation window size is then adjusted to encompass one full cycle of the waveform, e.g. by setting the far end thereof at the second zero crossing past the beginning thereof. The time of occurrence of the latest first morion detection is next compared (diamond 130) against a maximum allowable time, designated T , which reoresents the w approximate expected arrival time of the water or mud waves, the value being selected based on the transmitterreceiver spacing and an assumed water velocity of about 200 us per foot (656 με per metre). If T^ is exceeded ( of course, it is unlikely that x·'· it will be exceeded during an early pass through the loop being described), a meaningful value of shear velocity is unlikely to be obtained at this depth level, and diamond 240, described below, is entered. Decision diamond 190 is next entered and the index r (originally set to unity — block 116 above) is tested to see if it is greater than or equal to K, where K is the number of waveform cycles to be skipped, for example five cycles, from the original cycle which had been used to determine the compressional velocity or slope. This inquiry is represented by the diamond 190.

If the answer is in the negative, block 195 is entered, this block representing essentially the same function as the block 165 above; viz., of stripping away the R^ waveform within the present correlation window. 31ock 170 is then re-entered and, by action of the functions of blocks 170 and 175, the beginning of the correlation window is then moved to the zero crossing closest to the next first motion, and the size of the correlation window is readjusted, as previously described. This action continues, in the loop designated by reference numeral 199, until the requisite number of cycles (e.g. five cycles) have been skipped over.

The primary purpose of loop 199, in the present embodiment, is to save processing time since it is not generally expected that meaningful shear arrivals will occur within five cycles of the first strong compressional cycle. However, 8 2 5 2 it will be understood that, if desired, the technique to be subsequently described for determination of shear wave velocity or slope could be initiated directly after the output value of compressional velocity or slope is 'X determined.

When the requisite number of cycles have skipped, block 200 is entered and a fourfold correlation is performed with respect to Ene present correlation window. The fourfold correlation was described above in conjunction with block 120 of FIG. 4 and, in further derail, in conjunction wiEn the flow chart of FIG. 7. In the case of block 200, the initial estimated value of shear wave velocity or slope, designated ATse, is as derived above in conjunction with block 150. In other respects, the correlation may be the same as described with reference to FIG. 7, with Σ values being obtained at 10 ys spacings of slope and an interpolation being utilized to obtain a resultant provisional velocity or slope. In this case, the result is a provisional value for shear wave velocity cr slooe, designated as 4T1 . Since m had been initially set at unity (block 116), this is a first provisional value of shear velocity or slope, i.e. ΔΤ^.

Decision diamond 205 is next entered and the pro- - visional value of shear wave velocity or slope (the first provisional value for the first pass Enrough) is tested against Ene lower limit for shear wave velocity or slope (esrabiished above — block 156) as represented by the decision diamond 205. If ΔΤ1 is not below this lower limit, sm decision diamond 210 is entered and ΔΤ' sm is tested as againsr the expected velocity of the mud or water waves, ATW, which is typically about 200 useconds. If ATgm does not exceed ΔΤν, block 215 is entered, this block representing the acceptance of ΔΤ' as an acceptable provisional value of sm .- shear wave velocity or slope, designated ATsm< In other words, removal of the prime from ΔΤ' designates that it is sm within a range which makes it a feasible value for further consideration. The diamond 220 is then entered and Δτ is sm tested as against the value of shear wave velocity or slope which had been determined at the previous depth level. In rhe present embodiment, if ΔΤ is within 5 us of the sm previously stored value of shear wave velocity or slope, then the value ΔΤ is accepted (block 225) as the outnut sm shear wave velocity or slope at this depth level, designated as ΔΤ-. If, however, T is not within 3 μ s of the output » sm value at the previous depth level, diamond 230 is entered and inquiry is made as to whether a particular number cf values of T have now been obtained at the current death sm level. In the present embodiment, the particular number utilized is four. Thus, if four provisional values, Τ , sm have not yet been obtained, block 235 is entered, this block also being entered, as can be seen from the diagram, when the tests of diamonds 205 or 210 had indicated that Τ' was sm outside the acceptable range (as previously described).

In accordance with the function represented by the block 235, the beginning cycle of the waveform used in the most recent correlation is stripped away. 31ock 170 is then re-entered, via block 171, which represents the incrementing of index m, and operation of blocks 170 and 175 establish the new position of the correlation window with respect to the R waveform. This effectively means that both ends of the correlation window are advanced to the next zero crossing and the new correlation window overlaps the previous correlation window by one-half cycle. Testing^ against T is then performed (diamond 180), but the loom w 199 skipping routine is now ineffective since index r exceeds the maximum (diamond 190) . Block 200 is then entered with a resultant fourfold correlation which yields the next provisional value of shear wave velocity or slope, AT' > which is AT' in this example. The criteria of blocks sm s2 205 through 230 are next applied as previously set forth.

The described procedure continues until either an outout value of AT is obtained bv operation of the block s * 225, until four provisional values AT have been obtained sm (an affirmative answer to the inquiry of diamond 230) , or until the water arrival time T is exceeded (diamond ISO). w In each case, diamond 240 is entered and inquiry is made as to whether cr not there is an accepted value AT . If not, s block 245 is entered, this block representing the selection for AT of the value of AT which is nearest to the value of s sm AT that had been determined at the Drevious ceDth level, s Block 250 is then entered (block 250 also being entered from the affirmative answer branch of diamond 240) , this block representing the recording of the predetermined values of AT and AT for the present depth level. Block 110, at tne c s beginning, is then re-entered and the technique is repeated for the next depth level.

In the embodiment set forth in conjunction with FIG. 4, a plurality of provisional shear wave velocities are developed and an output shear wave velocity is determined as a function of the provisional shear wave · velocities. (An exception in this embodiment is when a particular orovisional shear wave velocity is sufficiently close to the value at the previous depth level whereupon it is immediately accepted.) Each of these orovisional shear wave velocities is determined by correlating relatively small portions of the waveforms over a limited range of slopes (of the order of 50 U s) , so prior art problems associated with multiple correlation peaks, that might occur using longer waveforms and wider slope ranges, are minimized. Also, individual portions of the waveform containing spurious information have less tendency to upset the selection of accurate output values since, in the described embodiment, selection is from among a plurality of provisional values and an individual spurious value is likely to be discarded without affecting the output values.

The comparison of provisional velocity values (both shear and compresssional) to the values obtained at the previous depth level provides a degree of protection against spurious perturbations in the output graphs, but genuine abrupt changes in formation characteristics (such as at a bed boundary) will be properly indicated within a few depth level increments since the technique allows for abrupt modifications which appear genuine (e.g. diamonds 145 and 220, 230) .

It will be understood that the technique of the invention of obtaining an output shear wave velocity from among individual provisional shear waVe velocities obtained by correlating individual relatively Short portions of the waveforms can be applied equally we determination of an output compressional wave velocity, ''-s. The invention has been described with reference to a particular embodiment, but variations within the scope of the claims will occur to those skilled in the art. For example, the provisional velocity values can be utilized, in the case of both compressional and shear wave velocity determination, without the need for utilizing information from the previously processed depth level, to arrive at a decision concerning velocity at the depth level being processed. In fact, the present invention has been found to operate advantageously in situations where prior art techniques exhibited inability to track shear or through casing signals arrivals when continuity. of these arrivals is absent. It will be understood that alternate techniques for rendering a decision, using the provisional velocity values obtained with the technique of the invention, can be set forth. For example, FIG. 10 illustrates the provisional compressional velocity values and provisional shear velocity values obtained at a particular depth level. The resultant clusters cf values can be used to render a decision concerning an output compressional and/or shear wave velocity graphically, or alternative automotive techniques car. be employed. In FIG. 3, there is shown a plot of shear velocity versus depth which is similar to Ene one of FIG. 1, but wherein the range of provisional values obtained at each depth level also appear in the plot in the form of a vertical line, at each depth level, which indicates the range of provisional values obtained. For example, the dashed block 101 in FIG. 4 and the loop 102 associated therewith could be used to repeat the compressional velocity-determining correlation of block 120 for different individual cycles (which may or may not overlap) near the beginning of the R^ waveform (instead of saving time by skipping through these cycles as described in conjunction with the loop 199). In such case, a plurality of provisional compressional wave velocities would be determined and an output compressional wave velocity may be obtained from these provisional values as is the case for the output shear wave velocity determination. Also, it will be understood that alternative relatively short portions of the waveforms, such as two cycle portions, could be used for the correlations. Alternate correlation techniques can also be employed, consistent with the principles of the invention.

Claims (32)

1. CLAIMS:1. A method for determining the velocity of acoustic wave energy propagating in formations surroundinga^borehole, comprising the steps of:5 placing in the borehole an acoustic logging tool having a source of acoustic wave energy and a plurality of acoustic wave energy receivers disposed at respective locations spaced from said source; causing said source to emit acoustic wave energy 10 into the formations surrounding the borehole; deriving a plurality of signals respectively representative of acoustic wave energy received by said plurality of receivers at said locations; correlating a first portion in time of the signal 15 from one of said locations with first portions in time of signals from the other locations the first portions of the signals from the other locations being determined from an assumed velocity and the correlation being performed for a number of different assumed velocities to obtain a 20 resultant first provisional velocity; correlating a second portion in time of the signal from said one of said locations with second portions in time of signals from the other locations, the second portions of the signals from the other locations being 25 determined from an assumed velocity and the correlation of second portions being performed for a number of different assumed velocities to obtain a resultant second provisional velocity; and generating an output velocity as a function of 30 said first and second provisional velocities.

2. The method as defined by claim 1 wherein said first and second portions of the signals are compressional wave components and the output velocity is an output compressional wave velocity, 5

3. The method as defined by claim 1 wherein said first and second portions of the signals are shear wave components and the output velocity is an output shear wave velocity.

4. The method as defined by claim 1 wherein said first 10 and second portions each have a duration of substantially an integral number of half-cycles, not exceeding two full cycles, of the signal from said one location.

5. The method as defined by claim 1 wherein said output velocity is determined for each of a number of adjacent 15 depth levels.

6. The method as defined by claim 5 wherein said generated output velocity takes into account the output velocity previously generated at an adjacent depth level.

7. The method as defined in claims 1, 4, 5 or 5 com20 prising correlating a third portion in time of the signal from one of said locations with third portions in time of signals from the other locations, the third portions of the signals from the other locations being determined from an assumed velocity and the correlation being performed for a 25 number of different assumed velocities to obtain a resultant output compressional wave velocity, wherein said first and 48253 second signal portions are later-occuring portions than said third signal portions, with said first'and second resultant provisional velocities being, respectively, first and second resultant shear wave velocities, and-χ said output velocity being an output shear wave velocity. The method as defined by claim 7 wherein said third portions of said signals are substantially equal in duration to the first and second portions. 9. The method as defined in claim 1, wherein each of said plurality of signals is comprised of components indicative, respectively, of a characteristic velocity of the formations, and said first and second signal portions in each signal are parts of the same component.

8. 10. The method as defined in claim 9, wherein one of said signal components is indicative of compressional wave velocity and another, later-occurring, of said signal components is indicative of shear wave velocity, and comprising: correlating a portion in time of said one signal component from the location closest to said source with a portion in time from said one signal component of each of the signals from the other locations, said portions of the signals from the other locations being determined from an assumed compressional wave velocity and the correlation being performed for a number of different assumed compressional wave velocities to obtain an output compressional wave velocity, said first and second signal portions being part of said other component of the signals, with the one of said locations being that closest to the source, the assumed velocity being the shear wave velocity, and the first and second provisional velocities and the output velocity being, respectively, first and second provisional shear wave velocities and an output shear wave velocity.

9. 11. The method as defined by claim 10 wherein said porticos of the cne signal component and the first and second portions of the other component each have a duration of substantially an integral number of half-cycles, not exceeding two full cycles of the signal from said closest location.

10. 12. The method as defined by claim 10 or 11 wherein said output velocities are determined for each of a number of adjacent depth levels.

11. 13. The method as defined by claim 12 wherein said generated output shear wave velocity takes into account the output shear wave velocity previously generated at an adjacent depth level.

12. 14. The method as defined by claim 12 or 13 wherein said generated output compressional wave velocity takes into account the compressional wave velocity previously generated at an adjacent depth level.

13. 15. Apparatus for determining the velocity of acoustic wave energy propagating in formations surrounding a borehole , comprising: 482 52 means arranged to derive a plurality of signals respectively representative of acoustic wave energy received at a plurality of locations in the borehole spaced from a source of acoustic wave energy; means arranged to correlate a first portion in time of the signal from one of said locations with first portions in time of signals from the other locations, the first portions of the signals from the other locations being determined from an assumed velocity and the correlation being performed for a number of different assumed velocities to obtain a resultant first provisional velocity; means arranged to correlate a second portion in time of the signal from said one of said locations with second portions in time of signals from the other locations, the second portions of the signals from the other locations being determined from an assumed velocity and the correlation of second portions being performed for a number of different assumed velocities to obtain a resultant second provisional velocity; and means arranged to generate an output velocity as a function of said first and second provisional velocities.

14. 16. Apparatus as defined by claim 15 wherein said first and second portions of the signals are compressional wave components and the output velocity is an output compressional wave velocity.

15. 17. Apparatus as defined by claim 15 wherein said first and second portions of the signals are shear wave components and the output velocity is an output shear wave velocity.

16. 18. Apparatus as defined by claim 15, 16 or 17 wherein said first and second portions each have a duration of substantially an integral number of halfcycles, not exceeding two full cycles, of the signal from said one location.

17. 19. Apparatus as defined by any one of claims 15 to 18 wherein said output velocity is determined for each of a number of adjacent depth levels.

18. 20. Apparatus as defined by claim 19 wherein said 10 generated output velocity takes into account the output velocity previously generated at an adjacent depth level.

19. 21. Apparatus as defined by claim 15 including means arranged to correlate a third portion in time of the signal from one of the locations with a third portion in time of 15 each of the signals from the other locations, said third portions of the signals from the other locations being determined from an assumed compressional wave velocity and the correlation being performed for a number of different assumed compressional wave velocities to obtain an output 20 compressional wave velocity; wherein said first and second signal portions are later-occurring portions than said third signal portions, with said first and second resultant provisional velocities being, respectively, first and second resultant shear wave velocities, and said output 25 velocity being an ouput shear wave velocity.

20. 22. Apparatus as defined by claim 21 wherein said first, second and third portions each have a dilation V of substantially an integral number of half-cycles/- not exceeding two full cycles, of the signal from said one x 5 location.

21. 23. Apparatus as defined by claim 21 or 22 wherein said output velocities are determined for each of a number of adjacent depth levels.

22. 24. Apparatus as defined by claim 23 wherein said 10 generated output shear wave velocity takes into account the output shear wave velocity previously generated at an adjacent depth level.

23. 25. Apparatus as defined by claim 23 or 24 wherein said generated output compressional wave velocity takes 15 into account the compressional wave velocity previously generated at an adjacent depth level.

24. 26. Apparatus for determining the velocities of different modes of propagation of acoustic wave energy propagating in formations surrounding a borehole, compris 20 ing: means arranged to derive a plurality of signals respectively representative of the wave energy received at a plurality of spaced locations in the borehole; means arranged to identify the location and duration of a first characteristic portion in time of the signal/from one of said locations; means arranged to correlate said first portion of the signal from said one of said locations with first portions in time of signals from the other locations, the first portions of the signals from the other locations being determined from an assumed velocity and the correlation being performed for a number of different assumed 10 velocities to obtain a resultant first provisional velocity; means arranged to identify subsequent characteristic portions in time of the signal from said one of said locations; means arranged to correlate said subsequent portions 15 of the signal from said one of said locations with subsequent portions in time of signals from the other locations, the subsequent portions of the signals from the other locations being determined from an assumed velocity and the correlation of each of said subsequent portions being performed for a 20 number of different assumed velocities to obtain resultant further provisional velocities; and means arranged to generate output velocities of said different modes of propagation as a function of said first provisional velocity and said subsequent provisional 25 velocities.

25. 27. Apparatus as defined by claim 26 wherein said first portions of the signals are compressional wave components, the subsequent portions of the signals are shear wave components, and said output velocities are compressional and shear velocities. X v

26. 28. Apparatus as defined by claim 26 or 27 wherein said first and subsequent portions each have a duration of substantially an integral number of halfcycles, not exceeding two full cycles, of said compressional and shear wave components, respectively.

27. 29. Apparatus as defined by claim 26, 27, or 10 28 wherein said output velocities are determined for each of a number of adjacent depth levels.

28. 30. Apparatus as defined by claim 29 wherein said generated output velocities take into account the output velocities previously generated at an adjacent 15 depth level.

29. 31. A method for determining the velocity of acoustic wave energy propagating in formations surrounding a borehole, substantially as hereinbefore described with reference to Figures 1 and 4A to 4D of the accompany20 ing drawings.

30. 32. A method for determining the velocity of acoustic wave energy propagating in formations surrounding a borehole, substantially as hereinbefore described with reference to Figures 1, 4A to 4D and 7 of the 25 accompanying drawings.

31. 33. Apparatus for determining the velocity of acoustic wave energy propagating in formations surrounding a borehole, arranged to operate substantially as hereinbefore ' described with reference to Figures 4A to 4D of the accompanying drawings.

32. 34. Apparatus for determining the velocity of acoustic wave energy propagating in formations surrounding a borehole, arranged to operate substantially as hereinbefore described with reference to Figures 4A to 4D and 7 of the 10 accompanying drawings.