Pulsed neutron well logging

Abstract: A method of measuring high vertical resolution earth formation thermal neutron capture cross sections of earth formations in the vicinity of a well borehole, comprising the steps of: (a) repetitively emitting in a well borehole relatively short duration bursts of fast neutrons which are rapidly moderated by interaction with nuclei of materials in the borehole and surrounding earth formations and slowed to thermal energy, creating repetitive time dependent distributions of the thermal neutron population in the borehole and surrounding earth formations; (b) detecting, as a function of depth, in the borehole during time intervals between said repetitive bursts of fast neutrons, radiations representative of the thermal neutron population in the borehole and surrounding earth formations in at least two time subintervals subsequent to each burst of fast neutrons and generating at least two count rate signals representative of said thermal neutron populations during said at least two time subintervals; (c) filtering said count rate signals to reduce statistical fluctuations in subsequent computations; (d) combining said at least two filtered count rate signals according to a predetermined relationship to obtain a neutron die-away based measurement of the formation thermal neutron capture cross section Σ F ; (e) selecting at least one of said at least two filtered gate count rate signals, or a sum of selected filtered count rate signals, to obtain a combined gate count rate signal G and using a predetermined functional relationship f(G) between G and the previously determined filtered formation capture cross section Σ F to obtain as a function of depth a normalization parameter C o ; (f) obtaining a combined count rate signal G′ from said at least one of said at least two count rate signals in step (b), prior to filtering in step (c), in a manner analogous to obtaining filtered count rate signal G; (g) determining from f(G′) and the normalization parameter C o in step (e) a high resolution thermal neutron capture cross section Σ HR with the approximate vertical resolution of the count rate signals in step (b); and (h) recording Σ HR as a function of borehole depth to produce an enhanced depth resolution measurement of thermal neutron capture cross section, relative to the Σ F measurement in step (d).

Description: This invention relates to in situ measurements of earth formations traversed by a well borehole and, more particularly, to the measurement of thermal neutron decay times or neutron lifetimes in earth formations in and around a well borehole, using pulsed neutron emissions. A well logging system which measures simultaneously the thermal neutron decay components of the earth formations in the vicinity of a well borehole and the borehole component itself, is known. U.S. patent specification No. 4,409,481 describes one successful commercial system for making such simultaneous measurements. Briefly in the invention described in the U.S. patent specification, a well logging tool is moved through a borehole. The tool includes a pulse source of fast neutrons and two radiation detectors. The pulsed neutron source generates a pulse of neutrons of approximately 14 Mev energy for a duration of between 10 and 150 microseconds at a pulse rate of approximately 1,200 pulses per second. The neutrons are introduced into the media comprising the well borehole and the surrounding earth formations and result in a thermal neutron population being generated from the slowing down of the, fast neutrons in the earth formation media and in the borehole media. After a very short pause to allow moderation of the fast neutrons following the neutron pulse, the detectors are gated on and capture gamma radiation resulting from the capture of thermal neutrons in the borehole, and the earth formations in the vicinity of the borehole, are measured essentially continuously until the next neutron burst is about to begin. Multiple time gates which occur during this essentially continuous interval are used for this purpose. In U.S. patent specification No. 4,409,481, the use of six such time gates is described. The gates, being essentially contiguous in time and of variable length, have the shorter gates closer to the neutron pulse and the longer duration gates further removed in time from the neutron pulse. The contiguous nature of the gates tends to reduce statistical errors, but neither the number of nor the continuous nature of the gates is essential to this invention, so long as the number and position of gates are adequate to provide data to solve for the desired components. The multiple time gate measurements of the counting rates of the gamma rays measured in each gate are supplied to a thermal neutron lifetime computer which computes formation and borehole neutron lifetime components by means of an iterative least squares fitting technique of the count rate data taken during four or more of the time gates following each neutron burst. The thermal neutron lifetime computer is enabled to calculate both the borehole thermal neutron lifetime component and the earth formation thermal neutron lifetime component simultaneously, and also simultaneously can compute the magnitude of the initial borehole and initial formation neutron components of thermal neutron population. Approximately once per second and for approximately 5% of the one second operating cycle, the neutron source pulsing sequence is turned off completely and the detectors are used to establish any relatively long lived background counting rate due to source neutron induced gamma ray activity within the gamma ray detector, the formation, the borehole, the logging sonde or any natural gamma radiation in the vicinity of the borehole. This background gamma radiation information is then properly normalized and subtracted from the count rates made in each of the time gates following the neutron bursts. The percentage of the one second cycle used for background can be varied but the approximate 5% amount mentioned has been found to be suitable for this purpose. A problem which has arisen due to the statistical nature of the measurement of the gamma rays generated by captured thermal neutrons has been that certain mathematical filtering techniques must be applied both to the raw measurement data prior to the data processing and then to the computed thermal neutron decay time or lifetime parameters for both the borehole and the formation which are produced by the logging system described previously. Moreover, the detectors in the system described for thermal multigate decay time logging are spaced as close as conveniently possible to the neutron source in the downhole logging instrument but, because of the relatively small diameter of the logging instrument, the near detector is spaced a distance of approximately 12 inches (30 cm) from the neutron source and the far detector is approximately 12 inches (30 cm) further from the neutron source than the near detector. A rough rule of thumb of vertical resolution of a logging instrument is that resolution is proportional to the distance between the source and the detector simply because of the physics of source propagation (of most types) and the received stimulus observed by the detector. This is not strictly true in the instance of pulsed neutron logging tools such as the thermal multigate decay log in that the vertical resolution of the long spaced detector is not quite twice that of the short spaced detector, but both have resolutions in the range of 18-30 inches (46 to 76 cm). It is thus seen that, because of the statistical count data filtering and the spacing of source to detectors in the thermal multigate decay time logging system, vertical resolution on the order of 4 to 5 feet (1.2 to 1.5 m) is to be expected with the conventionally processed data from this type of logging device. Because of the commercial realities of life, however, it is becoming more important to be able to evaluate thin bed stringers inside of shaley sand formations for potential hydrocarbon production. Higher vertical resolution than the 4 to 5 feet (1.2 to 1.5 m) resolution provided by conventionally processed data from the logging instrument in a thermal multigate decay system is quite desirable. In the present invention, such higher resolution is provided by unique processing techniques which are based upon the physics of the measurement, but which employ mathematical processing techniques in their implementation to improve the vertical resolution offered by the measuring instrument and processing software. Common practice to gain improved resolution for cased hole pulsed neutron logs has been to reduce or to eliminate filtering applied to the measured data from the detectors. While this does improve vertical resolution, it increases the statistical uncertainty to a level which is potentially undesirable in the acquired log. In the present invention, a technique is developed that makes use of the strong correlation between the formation cross section Σ F derived from the thermal multigate decay log processing and the overall system count rate. The count rate curve from, for example, gates 3 through 6 is much more repeatable than has been observed from the curve for Σ F . Hence, a Σ derived from the count rate curve will have less statistical variation than a Σ F curve produced from the raw unfiltered detector counts. It has been discovered that an unfiltered Σ F based on the gate 3 through 6 count rate of the thermal multigate decay log system has adequate repeatability for use directly in thin bed analysis. The gate 3 through 6 count rate, however, also varies with changes in borehole conditions as is noted in the aforementioned U.S. patent covering the thermal multigate decay log. It also varies with neutron generator output. Therefore, in the present invention a normalization procedure is used to correct for these non-formation cross section induced characteristics which may be termed a normalization procedure. The mathematical technique uses a continuous normalization of the Σ F derived from the unfiltered gate 3 to 6 count rate to the conventionally derived and filtered Σ F from standard thermal multigate decay logging processes. This procedure can yield a high resolution formation capture cross section Σ HR curve in which the statistical uncertainties are reduced to levels comparable to the standard filtered Σ F log but having a better thin bed resolution characteristic. In US-A-4267447, thermal neutron decay characteristics are measured by detecting indications of the thermal neutron concentration in the formation during a selected set of two different measurement intervals between neutron bursts and forming a ratio of such measurements, as accumulated for an appropriate length of time, from which the measured value of the decay characteristic is derived. The particular set of measurement intervals used is selected from among a number of possible sets as a function of a previously measured value of the decay characteristic. Each measurement interval set is used over only a specific range of decay characteristic values for which it has been determined, in accordance with a previously established relationship between the decay characteristic value and the ratio of the thermal neutron concentration measurements for the set, to afford enhanced statistical accuracy in the measured value of the decay characteristic. The measured value of the decay characteristic is then determined from the actual measured value of the ratio for the selected measurement interval set and the previously established relationship for that set. Each measurement interval preferably consists of a plurality of contiguous, discrete time gates whose durations and times of occurrence are adjustable in common by a finite number, e.g. four, of discrete scale factor values F. Each F value corresponds to a specific range of decay characteristic values, with the F value used for any given measurement being selected on the basis of a previously measured value of the decay characteristic. Regions of overlap are provided between the decay characteristic value ranges for adjacent F values so that F need not be changed merely as a result of statistical variations in the decay characteristic. GB-A-2173591 discloses a method of measuring earth formation characteristics, comprising in combination: repetitively irradiating the earth formations with pulses of neutrons; detecting gamma rays emanating from the irradiated formations; cumulatively counting for a preselected counting period the total gamma rays detected in a plurality of separate time bins, the cumulative counts for each pulse defining points on a cumulative count versus time curve which has a mathematical expression; and computing from the cumulative counts of the time bins variables in the mathematical expression which correspond to characteristics of the earth formations. According to the present invention, there is provided a method of measuring high vertical resolution earth formation thermal neutron capture cross sections of earth formations in the vicinity of a well borehole, comprising the steps of: (a) repetitively emitting in a well borehole relatively short duration bursts of fast neutrons which are rapidly moderated by interaction with nuclei of materials in the borehole and surrounding earth formations and slowed to thermal energy, creating repetitive time dependent distributions of the thermal neutron population in the borehole and surrounding earth formations; and (b) detecting, as a function of depth, in the borehole during time intervals between said repetitive bursts of fast neutrons, radiations representative of the thermal neutron population in the borehole and surrounding earth formations; characterized in that said radiations are detected in at least two time subintervals subsequent to each burst of fast neutrons and at least two count rate signals representative of said thermal neutron populations are generated during said at least two time subintervals; and wherein the method comprises the further steps of: (c) filtering said count rate signals to reduce statistical fluctuations in subsequent computations; (d) combining said at least two filtered count rate signals according to a predetermined relationship to obtain a neutron die-away based measurement of the formation thermal neutron capture cross section Σ F ; (e) selecting at least one of said at least two filtered gate count rate signals, or a sum of selected filtered count rate signals, to obtain a combined gate count rate signal G and using a predetermined functional relationship f(G) between G and the previously determined filtered formation capture cross section Σ F to obtain as a function of depth a normalization parameter C o ; (f) obtaining a combined count rate signal G′ from said at least one of said at least two count rate signals in step (b), prior to filtering in step (c), in a manner analogous to obtaining filtered count rate signal G; (g) determining from f(G′) and the normalization parameter C o in step (e) a high resolution thermal neutron capture cross section Σ HR with the approximate vertical resolution of the count rate signals in step (b); and (h) recording Σ HR as a function of borehole depth to produce an enhanced depth resolution measurement of thermal neutron capture cross section, relative to the Σ F measurement in step (d). The invention also provides a method of measuring high vertical resolution earth formation thermal neutron capture cross sections of earth formations in the vicinity of a well borehole, comprising the steps of: (a) repetitively emitting in a well borehole relatively short duration bursts of fast neutrons which are rapidly moderated by interaction with nuclei of materials in the borehole and surrounding earth formations and slowed to thermal energy creating repetitive bursts of thermal neutron population in the borehole and surrounding earth formation; (b) detecting, in the borehole during a time interval between said repetitive burst of fast neutrons, radiations representative of the thermal neutron population in the borehole and surrounding earth formations in at least six time subintervals subsequent to a burst of fast neutrons, and generating at least six count rate signals, G1, G2, G3, G4, G5, and G6, representative of said thermal neutron populations during said at least six time subintervals; (c) combining said at least six count rate signals according to an iterative fitting technique in which exponential relationships are assumed to exist for a borehole and a formation component of thermal neutron capture cross section to simultaneously obtain measurement signals representative of Σ F the formation component and Σ B the borehole component of thermal neutron capture cross section; (d) combining at least two of said gate count rate signals to obtain a combined gate count rate parameter (G3-6) and using a predetermined functional relationship between (G3-6) defined as f(G3-6) and the previously determined formation component of capture cross section Σ F to obtain a normalization parameter C o ; (e) determining as a function of the combined gate count rate parameter (G3-6) and the normalization parameter C o a high resolution thermal neutron capture cross section parameter Σ HR ; and (f) recording Σ HR as a function of borehole depth to produce an enhanced depth resolution measurement of thermal neutron capture cross section. In order that the present invention can be understood in more detail, more particular description of the invention, briefly summarized above, will be made with reference to the embodiments thereof which are illustrated in the appended drawings. The drawings illustrate only illustrative embodiments of this invention and are therefore not to be considered limiting of its scope. In the drawings: Figure 1 is a schematic representation of one embodiment of well logging system in accordance with concepts of the present invention; Figure 2 is a schematic block diagram depicting one embodiment of electronics system for a well logging system of the present invention; Figure 3 is a graphical relationship illustrating a composite thermal neutron population decay curve and time gates according to an embodiment of the present invention; Figure 4 is a graphical representation illustrating a telemetry sequence occurring as a function of time in connection with an embodiment of the present invention; Figure 5 is a graphical representation illustrating a cross-plot of short spaced gate 3 through 6 count rates versus Σ F in an interval having constant borehole conditions; Figure 6 is a cross plot showing short space gate 3 through 6 count rates versus Σ F in an interval having variable borehole conditions; and Figure 7 is a simplified flow diagram illustrating computational sequences for improving vertical resolution according to the concepts of the present invention. It will be recalled that the thermal multigate decay logging system of U.S. patent 4,409,481 assumes that the thermal neutron population in the vicinity of a well borehole following a burst of high energy neutrons may be described as the sum of a formation component, a borehole component and a background component. This may be expressed mathematically as set out in Equation 1: (1) C(t) = A B e -t/TB + A F e -t/TF + C B where C(t) is the count rate at any time t measured from a reference time; and A B and A F are constants which may be interpreted as A F being the initial amplitude of the formation component, A B representing the initial amplitude of the borehole component and wherein C B represents the background component. In equation (1) T B and T F represent the borehole and formation components of thermal neutron decay time, respectively. The formation cross section Σ F and Σ B are reciprocally related to the formation and borehole decay times as given in Equation 2: (2) T = 1/VΣ where v represents the velocity of thermal neutrons, and T and Σ are defined previously. It will further be recalled that the thermal multigate decay system of U.S. patent 4,409,481 measures the count rates of gamma rays resulting from the decay of nuclei formed by the capture of thermalized neutrons in the borehole and in the earth formation surrounding the borehole in six essentially contiguous time gates following each pulse of fast neutrons generated by the neutron generator. These count rates are then combined according to a predetermined mathematical relationship which involves the iterative solution of the exponential decay model given in Equation 1 by iterative techniques to obtain a borehole Σ B and a formation Σ F at each of the near and far detectors of the well logging tool. The time gating arrangement of the thermal multigate decay logging system is illustrated schematically in Figure 3. A fast neutron burst is first emitted and six essentially contiguous gates, labeled 1 through 6 in Figure 3, follow each neutron burst up to the beginning of the next neutron burst. The gates closer in time to the neutron burst are narrower in time, and the gates farther away are wider in time. This is to provide as nearly as possible by first approximation equal counting statistics for data from each of the six gates with the gates positioned as shown in Figure 3. It will be observed that the portion of the composite decay curve which may be approximated by the formation component exponential expression in Equation 1 is nearly a straight line on the logarithmic vertical axis of Figure 3. That is to say, if a sufficient time has elapsed after the neutron burst prior to the opening of gate 3, the gates 3, 4, 5 and 6 will contain primarily formation and background components, with only a modest contribution from the borehole component of the thermal neutron decay. It is also possible that the background component can be removed by using a subtraction technique well known in the art. The present invention exploits these features in making an approximation to be described. Referring now to Figure 1, a well logging system for performing thermal multigate decay time logging in accordance with the concepts of the present invention is illustrated schematically. A well borehole 10 is filled with a borehole fluid 11 and penetrates earth formations 20 to be investigated. A downhole well logging sonde 12 is suspended in the borehole 10 from a conventional armored logging cable 13 in a manner known in the art. The sonde 12 may be raised and lowered through the borehole as desired. The well logging cable 13 passes over a sheave wheel 14 at the surface. The sheave wheel 14 is electrically or mechanically coupled as indicated by a dotted line 15 to a well logging recorder 18 which may comprise a strip chart recorder or a magnetic tape recorder or both as known in the art. The flow of measurements made by the downhole sonde 12 may thus be recorded as a function of the depth in the borehole of the well logging sonde 12. In the downhole sonde 12, a neutron generator 21 is supplied with high voltage (approximately 100 kilovolts) for its operation by a high voltage power supply 22. Control and telemetry electronics 25 are utilized to supply control signals to the high voltage power supply and the neutron generator 21 and to telemeter information measured by the downhole instrument to the surface on conductors of the logging cable 13. In the sonde 12, and spaced from the neutron generator 21, are two radiation detectors 23 and 24 which may comprise, for example, thallium activated sodium iodide crystals which are optically coupled to photomultiplier tubes. The detectors 23 and 24 serve to detect gamma radiation produced in the surrounding formations 20 resulting from the action of the neutron generator 21 in emitting neutrons. A neutron shielding material 28 having a high density matter content or large scattering cross section is interposed between the neutron generator 21 and the dual spaced detectors 23 and 24 in order to attenuate direct irradiation of the detectors by neutrons emitted from the generator 21. Shielding 29 may also be interposed between the detectors 23 and 24 if desired. Upon activation of the neutron generator 21, a relatively short duration burst or pulse of neutrons of about 20 to 125 microseconds duration is initiated and is emitted into the well borehole 10 to pass through the borehole fluid 11, the steel casing 26 and cement layer 27 (assuming a cased well) surrounding the steel casing into the earth formations 20 being investigated. The neutron burst is rapidly moderated or slowed down by scattering interactions such that statistically essentially all the neutrons are at thermal energy in a relatively short time. The thermalized or thermal neutrons may then be captured by nuclei of constituents of the borehole 10 and the earth formations 20 and including the pore spaces contained in the earth formations. The capture of neutrons by nuclei of elements comprising the borehole constituents 11, 26, 12 and 27 and the earth formations 20 and their pore spaces result in capture gamma rays which are emitted and impinge upon detectors 23 and 24. A voltage pulse is produced from the photomultipliers of detectors 23 and 24 for each gamma ray so detected. These voltage pulses are supplied to the electronics section 25 and are counted in a digital counter and then telemetered to the surface via conductors 16 of the well logging cable 13. At the surface, a surface electronics package 17 detects the telemetered information from the downhole sonde 12 and performs the iterative least squares fitting technique to determine the parameters Σ F , Σ B , A F , and A B , as previously described with respect to Equation 1. The surface electronics package then supplies signals representative of the measured quantities to the recorder 18 where they are recorded as a function of borehole depth. Referring now to Figure 2, a schematic block diagram illustrates the electronics portion of the subsurface and surface electronic systems illustrated in more detail but still schematically. Power for the operation of the subsurface electronics is supplied by conductors of the well logging cable 32 to a conventional low voltage power supply 31 and a high voltage power supply 34. The high voltage power supply 34 may be of the Cockcroft Walton multistage type and it supplies approximately 100 kilovolts for the operation of the neutron generator tube 33. The neutron generator tube 33 is of the deuterium tritium accelerator type. An ion source 36, which is maintained at a potential near ground, is used to generate deuterium ions from a mixture of deuterium and tritium gas filling the envelope of tube 33. A replenisher heater 37 is impregnated with additional deuterium and maintains a pressure level of gas inside the tube 33 envelope sufficient to supply the ion source 36 with deuterium gas for ionization. A target 35 is impregnated with tritium and is maintained at a relatively high negative potential of 100 kilovolts. The ion source is controlled by an ion source pulser 41. When supplied with a relatively low level voltage pulse (2 kilovolts, for example), the ion source causes gas in the tube 33 envelope to become ionized and accelerated toward the target material 35. Upon impinging on the target material of target 35, the deuterium ions interact thermonuclearly with the tritium ions in the target to produce neutrons. The neutrons are emitted in a generally spherically, symmetrical fashion from the neutron generator tube 33 into the borehole 10 and surrounding earth formations 20 of Figure 1. A replenisher control circuit 39 is supplied with samples of the neutron generator target current by sampling circuit 38 and utilizes this to compare with a reference signal to control replenisher current and thereby the gas pressure in the envelope of the neutron generator tube 33. Timing circuits 42, which comprises a master timing oscillator operating at a relatively high frequency and an appropriate divider chain supply 1.25 kilohertz pulses to the ion source pulser 41 and also supply one second clock pulses to the neutron generator startup control circuit 40. Moreover, timing circuit 42 supplies 2.0 megahertz clock pulses to a micro processor and data storage array 44 and supplies timing pulses to the background circuit 45 and counters 52 and 53. Similarly timing signals are supplied to a pair of gain control circuits 48 and 49. During the background portion of the detection cycle, background circuit 45 is supplied with counts from the counters 52 and 53. This circuit also provides a disable pulse to the ion source 41 to prevent pulsing of the neutron generator during the background counting portion of the cycle. The background correction circuit 45 supplies background count information to the microprocessor and data storage 44. The background may be stored and averaged for longer periods. Since there are low discriminator thresholds, most gamma ray background is from neutron activation of the iodine in the detector crystal which has a half life of 27 minutes, better statistics in the subtracted signal result. It will be noted that adaptive or other filtering can be applied to background as well as capture gates prior to executing the iterative program. It should also be noted that if background count rate is not significant during the count rate gates used in the computation of the high resolution Σ F , then background correction may not be necessary to the implementation of this invention. The digital count information from the counters 52 and 53 in the background correction circuit 45 are supplied to the microprocessor and data storage circuit 44. These circuits format the data and present it in a serial manner to the telemetry circuit 43 which is used to telemeter the digital information from the counters and background correction circuit to the surface by conductors along the logging cable 32. At the surface, a telemetry interface unit 54 detects the analog voltage signals from the logging cable 32 conductors and supplies them to a telemetry processing unit 55 which formats the digital count rate information representing the counting rate from the counters 52 and 53 in the subsurface equipment in terms of the time gating scheme previously discussed with respect to Figure 3. Referring now to Figure 4, the telemetry stream just described along with synchronization or sync pulses is illustrated schematically. This would represent the voltage levels as seen on the conductors of the logging cable 32 of Figure 2. A sync pulse preceding each neutron burst is sent and then a signal indicating the neutron burst is sent. Count rate signals representative of the signals in gates 1 through 6 illustrated in Figure 4 as G₁ through G₆ then succeed as a function of time. These are digital numbers representative of the count rate summed over a plurality of counting intervals in each of the gates typically over a one second interval. The digital numbers representative of the count rates are supplied to a computer 56 which is programmed to interpret the six time gates and the background counting rate information in terms of thermal neutron decay times, or lifetimes, of the borehole and formation components from each detector. Output signals representing the parameters of interest are supplied from the computer 56 to a film recorder 57 and a magnetic tape recorder 58 for recording as a function of borehole depth. The computer 56 is programmed in accordance with the techniques shown in U.S. patent 4,409,481 and to be subsequently described herein to perform the high resolution determination of the Σ F . It should be noted that at this point the Σ F can be corrected for diffusion effects if desired. Further high resolution computations can be made with or without applying diffusion corrections, as desired. HIGH RESOLUTION PROCESSING Experience in logging many wells with the thermal multigate decay logging system has shown that a relationship exists between the detector count rates of the near and far detectors and the formation capture cross section, at least in an approximate sense. This relationship, in spite of its inherent borehole sensitivity, can be exploited to improve the repeatability of high resolution capture cross section processing. Even though Σ F is obtained from the same count rates as used in gates 3 through 6, it has much more statistical variation since it is derived from a least squares fit to the slope of a curve passing through each of the individual count gates. It is primarily statistically sensitive, not to the total of the counts in the gates but rather to the potentially low count rates in gates 5 and 6 which are far from the neutron burst. The gates 3 through 6 count rate curve is a much more repeatable curve than the processed Σ F since it is a direct measure of thousands of counts per second summed for all the gates used to compute Σ F . A Σ measurement based on a statistically superior overall count rate measurement, even when unfiltered, should have adequate repeatability for direct use in thin bed analysis. Several wells have been analyzed over intervals where the borehole conditions did not change significantly and a correlation between the gate 3 through 6 count rate and Σ F has been observed. For example, in Figure 5, a cross-plot of the short spaced detector gates 3 through 6 count rate versus Σ F in an interval with constant borehole conditions is illustrated. A general function of the form is given in Equation 3: (3) Σ F =f(G3-6) + C O where Σ F is the conventionally derived, decay based formation capture cross section; C O is a normalization parameter; and f(G3-6) represents a generalized function of the gate 3 through 6 count rates. Unfortunately, the gate 3 through 6 count rates also vary with changes in borehole conditions and the neutron generator output. A cross-plot is shown in Figure 6 which is an expansion of the depth interval plotted in Figure 5. In Figure 6, the borehole fluid changed from brine in the bottom (the portion labeled trend A) to a mixture of brine and oil, labeled trend B. The well logging instrument started in 17.78 cm (7 inch) casing and came into a 7.30 cm (2 7/8 inch) tubing with brine in the casing tubing anulus, see trend C. These changes in borehole conditions are evidenced by the shift of the trends along the x-axis of Figure 6. A normalization scheme to compensate for the non-Σ induced count rate changes has been developed and uses the normally filtered Σ F and gate 3 through 6 count rate data to solve for the slowly varying parameter C O as follows: (4) C O =Σ FM - f(G3-6) This normalization parameter is then used in Equation 3 with the unfiltered (quarter foot) gate 3 through 6 count rate data to perform the computation of a high resolution formation cross section. The output from this process is labeled Σ HR which has roughly the precision of the standard filtered Σ F log and which has a response to thin bedding similar to the Σ F derived from the unfiltered data which would otherwise have a high statistical variation without this unique normalization procedure. Formation diffusion effects can be compensated for during the normalization process since the standard diffusion corrected Σ F may be used as the filtered base for data normalization. Alternatively, diffusion effects can be corrected after the normalization process if the Σ F used in the normalization procedure is taken directly from the short spaced and/or long spaced formation decay component(s). Effects on Σ HR from direct borehole counts are reduced both by considering data only in gates 3 through 6 where the borehole counts are greatly reduced and also because the Σ F from the two exponential model of Equation 1 has only minimal borehole sensitivity. The actual selection of the specific gates used, however, is not critical so long as there is a good correlation between the selected count rate and the measured Σ FM . In addition, the Σ FM used in the normalization procedure does not neccessarily have to be derived using the technique of USP 4,409,481. Any conventionally produced decay based Σ may be employed. It is recognized that Σ HR changes caused by sudden anomalous temporary changes in the unfiltered gate 3 through 6 count rate due to packers located in the borehole, casing collars, etc. may not be adequately normalized using this procedure. This could create unusual Σ HR logs across depths containing such anomalies. To reduce these effects in the high resolution Σ HR presentation, we have chosen to place checks and bounds on Σ HR , based on the computed statistical uncertainty in the unfiltered decay based Σ F , and in the rate of change of Σ HR with depth. The computed Σ HR is compared to the Σ F computed conventionally from the unfiltered count rate decay data. These two Σ's will have the same vertical resolution, but the decay based Σ F will be more statistical, as discussed earlier. Therefore Σ HR should agree with Σ F to within the statistics of Σ F (unfiltered). We have defined an operator selected criterion for their comparison. If Σ HR falls within two standard deviations (selectable) of the unfiltered Σ F , the Σ HR is output. If Σ HR falls outside the statistical bound, then Σ F (or a minimally filtered Σ F ) is substituted for Σ HR over the interval where the two disagree, the assumption being that Σ HR is responding incorrectly to a borehole anomaly or generator output fluctuation, neither of which effects Σ F (unfiltered). Realistic bounds have also been placed on the maximum rate of change of the computed Σ HR as a function of depth (3.0cu/0.25ft.) (0.25 ft. = 0.076 m), to limit anomalies not trapped by this statistical comparison process. It may be desirable to indicate on the log as a quality control feature the intervals where these checks and bounds have been exceeded. Referring now to Figure 7, the high resolution Σ processing technique is illustrated in a flow diagram. This flow diagram is in the nature of a macro diagram illustrating the broad principles of the technique. Upon entry, the entire desired depth interval log is recomputed at block 71 conventionally with both normal filtering and no filtering using the raw telemetry data as input. Of course, if this data were recorded on the real time log and stored on tape, this step would not be necessary. The output from the process includes both filtered and unfiltered gates 3-6 count rates from each detector and Σ F curves. Control is then passed to block 72 in which the high resolution algorithm is calibrated over a selected depth interval. To perform this calibration, a log analyst chooses from either the short or long spaced gates 3-6 count rate and selects a Σ F computed directly from either detector or from the diffusion corrected Σ F to make the calculations. Borehole conditions should be consistent across the calibration well section, and the Σ F and gate 3 through 6 count rate curve should show good dynamic range for that section. The filtered gate 3 through 6 count rate data and the filtered Σ F data computed in the normal thermal multigate decay log manner are then used for the calibration interval to determine the best fit calibration coefficients in f(G3-6), Equation 3, using a standard error minimization routine. For quality control, these coefficients and a calculated fit error can be displayed for the log analyst to verify that the calibration is valid before continuing. Control is then passed to the block 73 in which a standard deviation is computed for the unfiltered Σ F . Σ HR is also computed over the entire depth interval in the well using the unfiltered gate 3-6 count rates in an algorithm according to Equation 3, as calibrated in block 72, and a normalization factor C O computed for each depth increment from the filtered gate 3-6 count rate and Σ F using Equation 4. These Σ HR computations are made for each successive depth increment using the techniques previously discussed. A rate of change limit is applied to Σ HR to reduce any unwanted anomaly spikes not handled by the standard deviation check, which is also conducted as discussed earlier using Σ HR , Σ F (unfiltered) and the standard deviation in Σ F (unfiltered) as inputs. Finally, control is transferred to block 74 in which environmental corrections (and possibly also formation diffusion corrections) to the high resolution Σ HR and unfiltered porosity ratios can be applied. If borehole parameters are specified then the same procedures may be used in correcting the Σ HR and unfiltered R N/F as are used in correcting the conventionally derived and filtered Σ F and R N/F curves of the thermal multigate decay process, where the R N/F is a porosity ratio of count rates in gates 3 through 6 taken from the near and far detectors over the same depth interval. This environmental procedure is outlined in Paper 55, 1988 SPWLA Symposium Transactions, entitled "Obtaining Intrinsic Formation Capture Cross Section With Pulsed Neutron Capture Tools". A limitation of the foregoing technique for determining high resolution Σ is that artificial thin beds can be created by count rate anomalies. Casing collars, packers, neutron generator output fluctuations, etc. can produce brief count rate changes in the gate 3 through 6 count rates which are not related to formation cross section effects. These effects are not completely normalizable using the procedures and the bounding limits and checks described earlier. However, these anomalies will generally cause Σ HR deflection toward higher neutron cross sections. Experienced log analysts will recognize such anomalous behavior when it occurs, especially if other logs such as a casing collar log are available, and to take it into account in interpreting the logs over the depth intervals in question.

Claims: A method of measuring high vertical resolution earth formation thermal neutron capture cross sections of earth formations in the vicinity of a well borehole, comprising the steps of: (a) repetitively emitting in a well borehole relatively short duration bursts of fast neutrons which are rapidly moderated by interaction with nuclei of materials in the borehole and surrounding earth formations and slowed to thermal energy, creating repetitive time dependent distributions of the thermal neutron population in the borehole and surrounding earth formations; and (b) detecting, as a function of depth, in the borehole during time intervals between said repetitive bursts of fast neutrons, radiations representative of the thermal neutron population in the borehole and surrounding earth formations; characterized in that said radiations are detected in at least two time subintervals subsequent to each burst of fast neutrons and at least two count rate signals representative of said thermal neutron populations are generated during said at least two time subintervals; and wherein the method comprises the further steps of: (c) filtering said count rate signals to reduce statistical fluctuations in subsequent computations; (d) combining said at least two filtered count rate signals according to a predetermined relationship to obtain a neutron die-away based measurement of the formation thermal neutron capture cross section Σ F ; (e) selecting at least one of said at least two filtered gate count rate signals, or a sum of selected filtered count rate signals, to obtain a combined gate count rate signal G and using a predetermined functional relationship f(G) between G and the previously determined filtered formation capture cross section Σ F to obtain as a function of depth a normalization parameter C o ; (f) obtaining a combined count rate signal G′ from said at least one of said at least two count rate signals in step (b), prior to filtering in step (c), in a manner analogous to obtaining filtered count rate signal G; (g) determining from f(G′) and the normalization parameter C o in step (e) a high resolution thermal neutron capture cross section Σ HR with the approximate vertical resolution of the count rate signals in step (b); and (h) recording Σ HR as a function of borehole depth to produce an enhanced depth resolution measurement of thermal neutron capture cross section, relative to the Σ F measurement in step (d). A method according to Claim 1, wherein the cross section Σ HR is determined according to a functional relationship of the general form Σ HR = f(G′) + C O where f(G′) is an empirically derived predetermined functional relationship and C O is said normalization parameter. A method according to Claim 1 or 2 wherein all steps are performed for count rate signals taken at each of two detectors, a near detector and a far detector, spaced longitudinally from the source of said repetitive neutron bursts, and wherein the two resulting Σ HR measurements are combined with information about formation porosity to obtain a formation diffusion corrected high resolution formation cross section. A method according to Claim 1, 2 or 3 wherein said predetermined functional relationship in step (e) is parameterized from filtered ΣF and G values observed over an interval of the well with reasonably constant borehole conditions. A method according to claim 1, 2 or 3, wherein said predetermined functional relationship is step (e) is parameterized from unfiltered ΣF and G′ values observed over an interval in the well with reasonably constant borehole conditions. A method according to claim 1, wherein radiations representative of the thermal neutron population in the borehole and surrounding earth formations are detected in at least six time subintervals subsequent to a burst of fast neutrons, and at least six count rate signals, G1, G2, G3, G4, G5 and G6, representative of said thermal neutron populations during said at least six time subintervals are generated; said at least six count rate signals are combined according to an iterative fitting technique in which exponential relationships are assumed to exist for a borehole and a formation component of thermal neutron capture cross section to simultaneously obtain measurement signals representative of Σ F the formation component and Σ B the borehole component of thermal neutron capture cross section; at least two of said gate count rate signals are combined to obtain a combined gate count rate parameter (G3-6) and using a predetermined functional relationship between (G3-6) defined as f(G3-6) and the previously determined formation component of capture cross section Σ F to obtain a normalization parameter C o ; a high resolution thermal neutron capture cross section parameter Σ HR is determined as a function of the combined gate count rate parameter (G3-6) and the normalization parameter C o ; and Σ HR is recorded as a function of borehole depth to produce an enhanced depth resolution measurement of thermal neutron capture cross section. A method according to Claim 6, wherein the parameter Σ HR is determined according to a functional relationship of the general form Σ HR = f(G3-6) + C O where f(G3-6) is a predetermined functional relationship empirically derived and C O is said normalization parameter. A method according to Claim 7, wherein C O the normalization parameter is determined according to a functional relationship of the general form of C O = Σ f - f(G3-6) where f(G3-6) is as previously defined and Σ F is the formation component at the corresponding depth as determined by combining said at least two count rate signals. A method according to Claim 8, wherein steps (d) and (e) are performed using the Σ F for the far detector to determine C O which is used with the combined gate count parameter (G3-6) from the near detector to define said normalization parameter. A method according to Claim 6, 7, 8 or 9, wherein steps (a) to (c) are performed for count rate signals taken at each of two detectors, a near detector and a far detector, spaced longitudinally from the source of said repetitive neutron bursts.

Inventors: SMITH HARRY D JR, SMITH MICHAEL P, WYATT DENNIS F JR, SMITH, HARRY D., JR., SMITH, MICHAEL P., WYATT, DENNIS F., JR.
Applicants: HALLIBURTON CO, HALLIBURTON COMPANY
Application year: 1990
Main classification: G01V5
IPC codes: G01V5