Proposals to use induced gamma ray spectroscopy for detecting carbon/oxygen ratio in a well borehole dates back several decades. Commercial operations in this field have dated from the mid-1970's. The primary applications today remain the cased hole determination of oil saturation in reservoirs of low or unknown water salinity. Recent efforts have expanded the use of the technique to geochemical analysis in open hole. Unfortunately, the intrinsic dynamic range of the carbon/oxygen ratio measurement is small. This requires very high statistical precision for reasonable oil saturation uncertainty. Until recently, the technique has been very slow and often required stationary measurements. This was in part because of the small intrinsic dynamic range of the carbon oxygen ratio, but, also because the detectors used in well logging instruments had significant limitations.
Over the past decade or so, several new high density scintillators have become available which can provide improved gamma ray detection capability. Newer carbon/oxygen and geochemical logging tools can use these new high density scintillators. One of these materials is bismuth germanate (BGO). The merit of the high density scintillator relative to conventional sodium iodide detectors (previously used in induced gamma ray spectrometers) is that it provides substantially better detection efficiency. That is, a larger fraction of the total gamma ray flux is in the full energy or first escape peak. Less gamma ray energy is disposed in the structureless Compton scattered low energy tail, in the more dense bismuth germanate detectors relative to sodium iodide. Additionally, in a well logging system according to concepts of the present invention, a different and unique gain stabilization system is employed which further insures more accurate energy representation of the detected gamma rays by more precisely controlling the system gain.
Brief Description of the Invention
In accordance with one aspect of the present invention, there is provided a system for use in cased well boreholes for measuring gamma ray spectra induced by the bombardment of earth formations in the vicinity of the well borehole with fast neutrons. The system comprises a neutron source, a scintillation detector of bismuth germanate (BGO) optically coupled to a photomultiplier, a pulse height analyzer, and timing circuitry.
The neutron source comprises an electronic source of monoenergetic fast neutrons which are emitted in repetitive pulses. The fast neutrons penetrate the borehole fluid, the well casing, and its cement sheath and enter the earth formations in the vicinity of the borehole. The fast neutrons are moderated by elastic and inelastic scattering to thermal energy and are captured by the nuclei of elements in and about the borehole.
The scintillation detector, which is optically coupled to the photomultiplier, detects gamma radiation induced by the neutrons emitted by the neutron source and in turn produces output electrical pulses having a magnitude representative of the energy of the gamma rays impinging upon the scintillation detector.
The pulse height analyzer, which is connected to the output of the scintillation detector, separates the output electrical pulses into a spectrum of components representative of the gamma ray energy spectrum of gamma rays impinging upon the scintillation detector. The pulse height spectra are accumulated in a downhole memory, which is part of the pulse height analyzer, and are transmitted to the surface via a digital telemetry system at regular intervals. A surface computer is programmed to analyze the spectra and to record various quantities of interest, such as the carbon/oxygen ratio, the silicon to calcium ratio, and other parameters of interest in neutron induced gamma ray spectroscopy logging.
The timing circuitry generates time gate signals to time separate the output electrical pulses of the scintillation detector into pulses representative of fast neutron induced gamma rays, thermal neutron capture gamma rays, and background gamma rays. The time gate signals are generated prior to the presentation of the output electrical pulses to the pulse height analyzer in a timing cycle comprising a fast subcycle portion and a slow subcycle portion. During the fast subcycle portion, output pulses representative of substantially only fast neutron induced gamma rays and thermal neutron capture gamma rays are presented to the pulse height analyzer. During the slow subcycle portion, output pulses representative of substantially only thermal neutron capture gamma rays and background gamma rays are presented to the pulse height analyzer.
The fast subcycle portion of the timing cycle includes at least one capture gamma ray time gate passing output pulses to the pulse height analyzer essentially only during quiescent intervals between repetitive pulses of fast neutrons from the electronic neutron source.
The slow subcycle includes plural capture gamma ray time gates and a background gamma ray time gate. The plural capture gamma ray time gates and the background gamma ray time gate pass the output pulses to the pulse height analyzer only during time intervals when the electronic neutron source is not producing fast neutrons.
The above referred to and other cited features of the present invention are better understood by reference to the detailed description to follow when taken in conjunction with the accompanying drawings in which:
Brief Description of the Drawings
FIG. 1 is a schematic diagram illustrating a well logging system in accordance with concepts of the present invention disposed in a well borehole;
FIG. 2A is a timing diagram showing the timing of a short spaced burst subsequence for both the neutron source and detector when used in accordance with concepts of the present invention;
FIG. 2B is a timing diagram showing the timing subsequence for a slow or long spaced neutron burst subsequence for both the neutron source and the detector when used in accordance with concepts of the present invention;
FIG. 3 is a diagram illustrating an inelastic gamma ray spectrum taken with an instrument according to the concepts of the present invention; and
FIG. 4 is a schematic diagram illustrating a capture gamma ray spectrum taken according to the concepts of the present invention and showing gain stabilization of energy windows.
Detailed Description of the Invention
Referring initially to FIG. 1, a well logging system in accordance with concepts of the present invention is illustrated in a schematic manner in a cased well borehole. A well borehole 10 lined with a steel casing 30 and filled with a borehole fluid 11 penetrates earth formations 17, 18, 19 and 20. A well logging cable 12, which is of the usual armored multiconductor cable type, supports in the borehole 10, a well logging sonde 21 for making neutron induced gamma ray spectral measurements in accordance with the concepts of the present invention. The well logging cable 12 passes over a sheave wheel 13 which maybe electrically or mechanically coupled, as indicated by a dotted line 14, to a spectral analysis computer 15 and a data recorder 16. The spectral analysis computer 15 may comprise, for example, any suitably programmed general purpose computer capable of receiving signals from the downhole equipment. The recorder 16 may comprise a strip chart recorder, as indicated having a record medium 17, or may comprise a digital tape or a disk file memory as desired.
The downhole sonde 21 comprises a power supply portion 22, a timing circuit portion 23, a telemetry portion 24, a gain stabilizer portion 25, a pulse height analyzer portion 26, a detector portion 27, a source portion 28, and a high voltage power supply portion 29. In operation, the power supply portion 22 takes electrical power from one or more of the conductors of logging cable 12 and converts it to appropriate voltage and current levels for supplying internal circuits of the downhole sonde 21 with electrical power as required. The timing circuits, as will be discussed in more detail subsequently in conjunction, with FIGS. 2A and 2B, provide timing pulses for opening and closing various data gates and thereby operating various detection schemes according to the concepts of the invention. The telemetry circuitry 24 is a conventional digital telemetry system capable of communicating to and from the surface supplied equipment at a rate of 50 kilobits/sec.
The gain stabilizer circuit 25, which will be described in more detail subsequently, operates in conjunction with the pulse height analyzer 26. The analyzer 26 receives signals from the detector portion of the tool 27. The detector signals are representative of gamma rays impinging upon the detector and which have been induced by a neutron bombardment from a neutron source 28 which is supplied for its operation by approximately 100 kilovolt high voltage power supply 29 or from a background gamma ray source. The neutron source 28 includes a neutron source controller 32 of conventional design which receives timing pulses from the timing circuit portion 23 to thereby control the emission of neutron pulses.
Referring now to FIGS. 2A and 2B, these Figures illustrate the timing applied to both the neutron source 28 and the detector circuits 26 and 27. FIG. 2A illustrates a fast timing subcycle, while FIG. 2B illustrates a slower timing subcycle. Referring initially to FIG. 2b, a timing cycle of 5 milliseconds in duration and repeated 50 times/sec. is illustrated. The neutron source 28 is turned on for a very short 14 MeV neutron pulse having a time duration of approximately 15 microseconds at the beginning of this 5 millisecond duration sub-cycle.
During the 1985 microseconds after the neutron source is turned off, a plurality of capture gamma ray time gates are opened to measure the thermal neutron population die-away, also commonly known as thermal decay, as a function of time. Up to six, eight or even more gates may be used as desired according to known techniques such as that of U.S. Pat. No. 4,424,444 (these gates are indicated as TD1-TD6 in the timing diagram). Using this technique the borehole and formation thermal neutron capture crosssections are measured. Thus the die-away of thermal neutrons produced by the moderation of the fast neutrons emitted by the neutron source 28 are detected for the first two milliseconds of this 5 millisecond subcycle. For the remaining 3 milliseconds of the subcycle portion of FIG. 2b, then, the detector 27 is used to detect background gamma radiation due to naturally occurring gamma rays in the earth formations surrounding the well bore. A background gate (indicated as BKGD in the timing diagram) provides the measurement of the background gamma radiation. All measurements may thus be corrected for background based on this background measurement by subtraction of the appropriately scaled count rate.
During the fast pulse portion of the timing diagram illustrated in FIG. 2a, it may be observed that the neutron source 28 is pulsed on for fifteen microsecond duration pulses (indicated as NB in the timing diagram). In FIG. 2a, while the neutron bursts starting at time T=zero and T=100 microseconds are taking place an inelastic scattering gamma ray time gate (indicated as IG in the timing diagram) is opened allowing pulses from the detector 27 occurring during the neutron burst to be gated to the pulse height analyzer circuitry 26. At the close of the neutron burst, the detector 27 is gated off for approximately 5 microseconds. This allows thermalization of the fast neutrons. Then a first capture gamma ray time gate (indicating as CG1 in the timing diagram) having a duration of thirty-five 4 microseconds (from 20 to 55 microseconds) is opened. A second capture gamma ray time gate (indicated as CG2 in the timing diagram) having a duration of thirty-five microseconds (from 55 to 90 microseconds.) is opened thereafter. The entire short or fast pulse cycle is then repeated beginning at the time labelled 100 microseconds of FIG. 2a and the neutron source 28 and detector timing are repeated as previously described for the next 100 microseconds. The entire sequence which is illustrated in FIG. 2a is repeated for a period of 1500 microseconds. The short or fast pulse subintervals of 100 microseconds duration are time interlocked with the longer 5 millisecond intervals illustrated in FIG. 2B. Timing pulses suitable for this purpose are produced in Timing Circuit section 23 of FIG. 1. In order to correct all measurements for background radiation the background pulses output to the pulse height analyzer are appropriately analyzed and scaled for different gate widths and then subtracted from the spectral data gathered during the fast subcycle time gates and the total count rates gathered in the slow subcycle time gates prior to the further processing of the spectral data and count rates.
Referring now to FIG. 3, an inelastic gamma ray energy spectrum which comprises gamma ray pulses detected during the gates labelled IG (inelastic gamma rays) in FIG. 2a and supplied to the pulse height analyzer 26 during this time when the neutron bursts are on is shown schematically for three different formations (for ease of illustration the spectra are shown both at a normal scale and at 5X scale to enable a clear illustration of all of the peaks). It will be observed that in the inelastic gamma ray spectrum that inelastic scattering peaks attributable to carbon, oxygen, silicon and calcium appear in the inelastic gamma ray spectrum. The inelastic spectrum of course is caused by a fast neutron scattering from the nucleus of an element in the vicinity of the borehole and exciting these nuclei to higher energy levels thereby losing some of the energy of the scattered neutrons (ie. scattering with energy loss). The excited nuclei returning to the ground state emit gamma rays characteristic of each of the nuclear species. Energy windows labelled C, O, Ca, Si are shown in the spectrum of FIG. 3 for the carbon, oxygen, silicon and calcium inelastic scattering energy peaks.
Referring now to FIG. 4, a capture gamma ray spectrum is shown schematically for three different formations (for ease of illustration the spectra are shown both at a normal scale and 5X scale to enable a clear illustration of all of the peaks). The capture gamma ray spectrum is taken in gates labelled CG1 and CG2 in FIG. 2a. This spectrum is attributable to the emission of gamma rays by the decay to the ground state of nuclei which have been excited by the capture of thermalized neutrons. As it takes perhaps several microseconds for all of the neutrons to become thermalized following the emission of a fast neutron burst from the source 28 of the instrument 21 shown in FIG. 1, the waiting period of 5 microseconds between the end of the neutron burst to the beginning of the capture gamma ray gates illustrated in FIG. 2a allows for thermalization of the fast 14 MeV neutrons emitted by the source 28 to achieve a thermal neutron population.
Referring again to FIG. 4, it will be seen that the iron peak (labelled Fe in FIG. 4) occurs at an energy of 7.65 MeV. Energy windows labelled W1 and W2 are positioned to straddle the Fe peak and are used for system gain control purposes. The iron peak will of course always be present in a cased well borehole in which the pulsed spectral gamma ray tool according to the concepts of the present invention is deployed.
The iron edge is always a prominent feature in cased well borehole spectrum, because the casing and the tool housing both contain iron. The iron peak is located at 7.65 MeV and is discernable, but not very prominent, in the capture gamma ray spectrum of FIG. 4. This peak becomes less prominent if chlorine infested or salt water is present in the earth formations surrounding the borehole. However, all of the capture gamma ray spectra show a precipitous drop at energies just above the iron 7.65 MeV line. Thus, two equal width adjacent windows labelled W1 and W2 spanning this edge are indicated in FIG. 4. Equality of width of these energy windows is not a critical requirement but is only made for convenience.
A spectral procedure for gain stabilization is based upon the count rates in the energy windows labelled W1 and W2. This procedure which is implemented in a microprocessor program in gain stabilizer 25 adjusts the photomultiplier high voltage (contained in detector 27) to maintain the ratio of counts in the energy window labelled W1 to that in energy window labelled W2 equal to the integer two. Again, this value of 2.0 is not particularly necessary as any value in the range from two to three would probably work just as well. It is merely convenient to use. Thus, the gain of the photomultiplier portion of detector 27 is increased whenever the ratio becomes greater than 2.0 and the gain of the photomultiplier is decreased whenever the ratio becomes less than 2.0. This places the juncture of the two windows near the upper count rate portion of the iron (Fe) edge. Because there are two gain settings where the window ratio will be 2, a start up procedure is required to assume selection of the correct gain value. This is accomplished by raising the photomultiplier gain until the count rate in the upper window W2 exceeds a predetermined minimum value prior to beginning the gain stabilization process as previously described.
Laboratory and field experience have demonstrated the utility and reliability of this technique of gain stabilization. Important and innovative features of this gain stabilization are that no isotopic source is required which would contribute a background to the primary measurement. Also, the stabilization is made on a feature of a capture gamma ray spectrum which is always present and which is unperturbed by formation matrix composition or borehole or formation fluid composition. Finally, stabilization of the pulse height analyzer spectral gain on a high energy feature of the spectrum renders the spectrum much less sensitive to zero offset.
It should also be noted that while the preferred embodiment previously described varies the high voltage on the photomultiplier tube itself inorder to effect gain stabilization, that this stabilization technique could equally well be applied to a variable gain amplifier which would come subsequent to the photomultiplier tube in the electronic circuitry comprising the detector 27 portion of the sonde 21 of FIG. 1.
The foregoing descriptions may make other alternative embodiments of the invention apparent to those having skill in the art. It is the aim of the appended claims to cover all such changes and modifications as fall within the true spirit and scope of the invention.