US 4,864,129 AGrant
Logging Apparatus and Method
Issue Date:1989-09-05
•16 Claims
•5 Drawing Sheets
Abstract
A method of calibrating a radiation detector, positioned in apparatus including a source of gamma radiation such that the detector may sense radiation from the environment of the apparatus due to radiation emitted by the source and scattered by the environment, includes utilizing test environments including formations of known bulk densities. The apparatus may be positioned within the test environments at different distances from the surfaces of the formations. Count rates of the detector are compared to the formation bulk densities for the distances of the apparatus from the formation surfaces to calibrate the detector relative to the radiation source.
Metadata
Assignee
- Baroid Technology, Inc.
Inventors
- William C. Paske
- Richard F. Roesler
Application Information
Application Number:US 2439042
Filing Date:1988-09-13
Priority Date:1986-06-11
Art Unit:256
Classifications
IPC:
G01V 512
Field of Search:
250378252.1;269;262;264;265;266;25489
Patent Drawings (5 sheets)
Description
Background of the Invention
1. Field of the Invention
The present invention pertains to techniques for making measurements of subsurface conditions in and surrounding a borehole such as a well or the like. More particularly, the present invention relates to method and apparatus for logging a borehole traversing an earth formation to simultaneously determine two or more downhole parameters, such as a transverse dimension of the borehole, the density of the surrounding formation and the density of the borehole fluid, including drilling mud, and provides method and apparatus which can perform this logging generally while drilling the borehole through the formation.
2. Description of the Prior Art
When drilling a borehole in the earth in search of petroleum, it is necessary to obtain as much information as possible regarding the nature and structure of the formations through which the borehole is passing. This information is necessary to enable the drilling operator to determine the progress of the drilling operation and to control its direction so as to intercept the pay zone. In the past, most of the necessary measurements have been made by pre-boring geological surveying techniques and then by wireline logging of the borehole after it has been drilled. This approach has a number of obvious disadvantages including loss of drilling time, the expense and delay involved in tripping the drill string so as to enable the wireline to be lowered into the borehole and both the build up of a substantial mud cake lining the borehole and invasion of the formation by the drilling fluids during the time period between drilling and taking measurements. An improvement over these prior art techniques is the recently developing art of measuring-while-drilling in which many of the characteristics of the formation are determined substantially contemporaneously with the drilling of the borehole. Measuring-while-drilling logging either partly or totally eliminates the necessity of interrupting the drilling operation to remove the drill string from the hole in order to make the necessary measurements by wireline techniques.
In addition to the ability to log the characteristics of the formation through which the drill bit is passing, obtaining this information on a real time basis provides substantial safety advantages for the drilling operation. Change in the density of the drilling fluid or the density of the formation or diameter of the borehole, for instance, would indicate conditions which possibly would require immediate attention of the driller in order to prevent a possible blowout. The decrease in the density of the drilling fluid might indicate influx of gas into the borehole from the surrounding formation. It would therefore be necessary to take prompt corrective action in order to prevent a blowout, for example by changing the density of the drilling fluid. With the previous wireline techniques, tripping of the drill string under these conditions could greatly increase the chances of a blowout occurring.
As mentioned above, one of the more important pieces of information to ascertain downhole is the density of the formation through which the borehole is passing. A known technique incorporates the use of gamma ray density probes which are devices generally incorporating a gamma ray source and at least one gamma ray detector which is shielded from the source and which during operation of the probe counts the photons emanating from the source and interacting with the electrons of the material of the formation primarily by Compton scattering. The percentage of photons emitted from the source which eventually pass to the detector after having undergone Compton scattering through the formation depends upon the density of the formation. The photons reaching the detector are counted by standard associated signal processing and data counting equipment.
One of the major difficulties encountered by the previously known density measuring devices is the requirement that the logging device physically contact the formation at the borehole wall. This requirement was necessitated by the fact that all of the known devices were wireline devices and were employed hours or even days after drilling. During this delay, the drilling fluids can both invade the formation and build up a coating of substantial thickness on the borehole wall, either one of which would directly affect the accuracy of measurements taken. It would not be possible to trip a drill string and lower a wireline device fast enough to avoid invasion and/or mud cake problems. These problems are not as severe in a measuring-while-drilling situation since there most likely is no significant amount of mud cake deposited on the borehole walls in a drilling situation, and measurements would generally be taken within an hour of drilling through a location to be measured. Some of the attempts to compensate for the mud cake problem in wireline measurements have included the use of two different detectors spaced axially along a deployable pad which is pressed into and plows through part of the mud cake. The near detector receiving radiation which is scattered partially from the mud cake provides a mud cake correction to the further spaced detector receiving radiation which has scattered principally from the formation. Such multiple detectors are usually used in combination with a complex collimation scheme to narrowly define the emitted beam and direct it into a specific region of the formation and to receive, at a particular detector, only radiation coming from a particular region of the formation.
An example of a prior art wireline density probe which is claimed to function regardless of the thickness and the chemical composition of the materials that are located between the density probe and the samples is discussed in U.S. Pat. No. 3,846,631. The disclosed technique comprises passing two gamma ray beams from two intermittently operated sources into the formation, receiving the radiation back scattered from each of the two sources by two separate detectors, and building ratios of products of the four separate counting rates in such a manner that the numerical result is an indication of the density of the formation. The two-detector probes must be deployed against the borehole wall, as previously discussed, and the spacing between the detectors is a critical dimension. Any non-uniformity of the formation materials or the borehole wall between the detectors will cause an erroneous result.
There are currently no known patented gamma radiation density detecting devices which operate in a measuring-while-drilling condition.
Prior art calipers are usually mechanical devices which require physical contact with the borehole walls. This requirement is very difficult to meet while drilling without affecting the steering of the bottom hole assembly. Acoustical calipers would be equally difficult to use during drilling due to the noise caused by the drilling itself.
It would be very advantageous to overcome the limitations and inaccuracies of the prior art by having a system including a method and apparatus for measuring the density of formations while drilling a borehole through the formations without the necessity of defining narrow bands of the formation, performing collimation of the radiation, or physically contacting the borehole walls.
The instant invention will obviate the need for additional caliper measuring devices, or assumptions regarding the borehole dimensions, or assumptions regarding the composition of the formation being logged, or information on the formation derived from secondary sources, or assumptions concerning the fluids in the borehole at the time of the measurement.
Summary of the Invention
The present invention provides measuring apparatus for use in determining downhole parameters associated with a borehole, including a generally elongate housing carrying a source of gamma radiation and at least two pairs of radiation detectors. Each pair of detectors is positioned with the detectors diametrically opposite each other on opposite sides of the longitudinal axis of the housing. The geometric configuration of the combination of the first detector pair with the radiation source is different from the geometric configuration of the combination of the second pair and the radiation source. The radiation source and the detectors are positioned so that radiation from the source may be emitted into the borehole and surrounding formation, and scattered radiation may be detected by the detectors. The detectors produce output signals which may be processed to provide data signals indicative of the radiation detected by each of the detectors. Using such data signals, two or more downhole parameters indicative of the borehole and/or its surroundings may be determined. Additional pairs of detectors, diametrically opposed across the longitudinal axis as described above and in unique configurations with the radiation source, may be employed whereby additional downhole parameters may be simultaneously determined.
The housing carrying the radiation source and detectors may be positioned in the borehole by means of a wireline or other flexible line extending to the surface, or may be included in a pipe string used i drilling, for example. In the latter case, the apparatus may be operated in the measuring-while-drilling mode, that is, the detectors may sense the scattered radiation while the pipe string is stationary during an interruption in the drilling, or while the pipe string is being rotated, even while drilling is taking place. Similarly, the housing depending from the wireline may be intentionally rotated, or not, during data taking procedures.
The detectors may be calibrated on the tool utilizing formations of known bulk densities, even including liquids as calibration formations.
According to the present invention, a borehole may thus be logged to simultaneously obtain radiation count data for use in determining two or more parameters indicative of the downhole environment. Such parameters may include the density of the borehole fluid, the bulk density of the formation at the level of the housing where data is acquired, or a transverse dimension of the borehole at such level.
Brief Description of the Drawings
FIG. 1 is a generally schematic side elevation, partly in section, of the lower end of a drill string including a downhole sub incorporating one embodiment of the present invention;
FIG. 2 is an enlarged, generally schematic transverse section, taking along line 2--2 of FIG. 1 and partly broken away to show relative placement of radiation detectors and a radiation source;
FIG. 3 is a generally schematic transverse section similar to FIG. 2, but showing details of another embodiment of the present invention;
FIG. 4 is a partial perspective schematic illustration of a portion of the instrument package insert of another embodiment of the present invention, showing staggered axial placement of detectors;
FIG. 5 is a generally schematic transverse section similar to FIGS. 2 and 3, but showing the circumferential placement of the detectors of the embodiment shown in FIG. 4;
FIG. 6 is an end view of a detector, illustrating how such a detector may include a bank of individual sensing devices;
FIG. 7 is a generally schematic transverse section of a drill string, including the present invention, positioned in a borehole taken along line 7--7 of FIG. 1, showing standoff distances A and B, the tool diameter TD, and the borehole diameter BD;
FIG. 8 is a graphical representation of typical calibration curves for detectors of the present invention;
FIG. 9 is a plot of slopes of calibration curves of the type shown in FIG. 8, as a function of standoff measurements;
FIG. 10 is a graphical representation of the dependence of the apparent formation density on the geometric constant as a function of the borehole fluid density;
FIG. 11 is a graphical representation of the dependence of the apparent formation density on the geometric constant as a function of the borehole diameter; and
FIG. 12 is a graphical representation showing the correction between the apparent formation density observed and the true formation density, as a function of the borehole diameter.
Description of the Preferred Embodiments
Measuring apparatus according to the present invention is shown generally at 10 in FIG. 1 housed in a sub 12 included within a drill string 14 a short distance above the drill bit 16 in a borehole. The sub 12 may be provided by a modified drill collar having an internal throughbore 18 sufficiently broad to receive an instrument package insert 20 and still provide an elongate throughbore 22 which may align with the internal bore of other drill string components and effectively provide an uninterrupted flow passage for drilling fluid from the surface down to the drill bit 16. The sub 12 is further modified to include at least one gamma radiation source 24 comprising a threaded insert received within a threaded opening in the side wall of the sub. Such a source 24 of gamma radiation may be any conventional well logging source, such as cesium-137.
The measuring system 10 is constructed to generate a series of signals for telemetry to the surface, either by electrical communication or mud pulsing, or recording on a downhole recording system for retrieval upon withdrawal of the apparatus from the well bore. As discussed more fully hereinbelow, detectors are included whereby radiation may be sensed and the radiation counts provided in the form of electrical signals indicative of the environment surrounding the measuring apparatus 10. The measuring system 10 measures averages of quantities that a affect the scattering of gamma radiation, and thus measures parameters that affect the average density surrounding the measuring system. For example, the signals may include information concerning the formation density of the earth formation adjacent to the borehole, the average borehole diameter and the density of the fluid, generally drilling mud, located in the annular region within the borehole and surrounding the sub 12. The signal processing, telemetry and recording systems that may be included as part of or ancillary to the measuring apparatus 10 may be of conventional design, and are therefore not described in detail.
The measuring apparatus 10 includes within the instrument package insert 20 two pairs of detectors, although additional detector pairs may be included. In each instance, each set, or pair, of detectors features the two detectors positioned mutually diametrically opposed across the longitudinal axis of the sub 12 and the instrument package insert 20. For example, in the embodiment illustrated in FIGS. 1 and 2, two detectors 26 and 28 are disposed to the right and left of the central longitudinal axis of the sub 12 as viewed, and directly thereabove is a second pair of detectors 30 and 32 similarly situated to the right and left of the sub longitudinal axis, respectively, so that detector 30 is directly above detector 26 and detector 32 is directly above detector 28. It will be appreciated that the right side detectors 26 and 30 are centered circumferentially at an azimuthal angle of 90.degree. from the circumferential position of the gamma radiation source 24, and that the left side detectors 28 and 32 are each similarly positioned circumferentially at an azimuthal angle of 90.degree. from the circumferential location of the source, but in the opposite sense from the right side detectors. Consequently, the right side detectors 26 and 30 and the left side detectors 28 and 32 are symmetrically disposed relative to the gamma radiation source 24. However, the upper detectors 30 and 32 are axially displaced farther from the source 24 than the lower detectors 26 and 28, and therefore are, in general, farther from the source.
Another arrangement of two pairs of detectors according to the present invention is illustrated in FIG. 3 wherein a first pair of detectors 26' and 28' is provided with the detectors again positioned diametrically on opposite sides of the longitudinal axis of the sub 12, and a second pair of mutually diametrically opposed detectors 30' and 32' is provided with the detectors positioned directly above the lower detectors 26' and 28', respectively. In this case, however, while again the upper detectors 30' and 32' are axially displaced farther from the source 24 of gamma radiation than the lower detectors 26' and 28', neither of the pairs of detectors is positioned symmetrically relative to the source. One detector of each of the pairs of detectors, namely the detectors 26' and 30', is centered aligned circumferentially above the source, and the other detector of each of the pair of detectors, namely detectors 28' and 32', is centered circumferentially at an azimuthal angle of 180.degree. from the source. Consequently, each of the pairs of detectors is oriented asymmetrically relative to the radiation source 24.
Another arrangement of two pairs of diametrically opposed detectors is illustrated in FIGS. 4 and 5 wherein it may be seen that a first pair of detectors 26" and 28" is arranged in the instrument package insert 20" asymmetrically relative to the gamma radiation source 24 with the detector 26" centered circumferentially above the source and the detector 28" centered 180.degree. circumferentially relative to the source. The second pair of detectors is arranged symmetrically relative to the source 24, with one detector 30" centered 90.degree. circumferentially to the right of the source as illustrated, and the second detector 32" centered 90.degree. circumferentially to the left of the source. However, as indicated in FIG. 4, the symmetrically positioned detectors 30" and 32" are axially displaced farther from the axial level of the radiation source 24 than are the asymmetrically arranged detectors 26" and 28". The relative difference in axial displacement of the asymmetrically positioned detectors 26" and 28" compared to that of the symmetrically positioned detectors 30" and 32" need only be a few centimeters, but could be as large as tens of centimeters. Further, in yet another array of detectors, the difference in axial displacement of the four detectors 26"-32" may be zero.
It will be appreciated that still another formation of the two pairs of detectors may be achieved by rotating the position of the source 24 in the embodiment of FIGS. 4 and 5 90.degree. to the right or left so that the radiation source is circumferentially aligned with one or the other of the detectors 30" and 32" which are axially displaced farther from the radiation source than are the remaining two detectors 26" and 28".
In each case, appropriate material for shielding the detectors against direct radiation from the gamma ray source 24 may be utilized. For example, as shown in FIG. 1, a generally cylindrical shielding construction 34 may be positioned within the array of detectors, and an annular shielding ring may underlie the detectors with a further cylindrical shielding shaft 38 extending downwardly and even radially within the position of the source 24 to minimize the detection of gamma radiation from the source directly to the detectors by "streaming" of the radiation axially along and through the sub 12 and the bore 22. The material used to provide the shielding 34-38 may be any appropriate gamma radiation shielding material, such as tungsten.
Each of the detectors discussed herein is illustrated schematically, or in block form, being shown generally as an arcuate oblong object. It will be appreciated that any appropriate gamma radiation detector may be utilized, and in each case the detector may comprise any number of gamma radiation sensors appropriately connected to provide a single output signal indicative of the radiation detected by the total array of sensors. Typical conventional gamma radiation sensors which may be utilized include Geiger-Mueller tubes, scintillators such as sodium iodide detectors, or solid state sensors such as of the germanium family of detectors. FIG. 6 illustrates how the detector 26, as an example, may include a cluster of generally cylindrically shaped Geiger-Mueller tubes 40, positioned within the region designated by the schematic shape otherwise used to indicate the detector 26.
Each of the pairs of detectors in any of the detector arrangements discussed herein includes two detectors that are positioned symmetrically relative to each other about the longitudinal axis of the sub 12 at equal azimuthal angles of separation from one another, that is, 180.degree. on diametrically opposite sides of the sub. The output signals generated by the sensing devices of each of the two detectors in such a pair are electronically processed to produce data signals proportional to the count rates of the detector sensors and which carry information about the environment in which the sub is located. This geometry of detectors so disposed may provide accurate measurements of the sub surroundings even with the sub 12 eccentered within the borehole because the distances between the eccentered sub and the borehole wall may be compensated for mathematically, at least in some cases. Alternatively, the measuring apparatus may be rotated while radiation is detected so that each data signal is proportional to the average readings of radiation detected by the respective detector sweeping in a 360.degree. circle around the borehole, and therefore carries information of the average values of the borehole parameters for the environment surrounding the sub 12.
As an example, the formation density could be obtained utilizing a pair of opposed detectors, such as the lower detectors 26 and 28 in the embodiment of FIGS. 1 and 2. If two pairs of detectors are used, two parameters may be independently determined, such as formation density and the diameter of the borehole. If three pairs of detectors are utilized, three logging parameters may be determined, such as formation density, borehole diameter and drilling fluid density. For each additional pair of opposed detectors added, an additional parameter of the sub environment may be determined.
In practice, whether in a measurement-while-drilling operation or with the measuring apparatus in a sonde (not shown) suspended by a wireline or other flexible line for example, at any given depth in the borehole the electronics operating the detectors may be turned on and radiation counts obtained. Gamma radiation from the source 24 interacts with the drilling fluid in the borehole surrounding the sub 12, and with the formation in the vicinity of the borehole at the level of the sub 12. In the latter case, the radiation interactions occur not only with the solid structure of the formation but with the fluid present in the interstitial spaces within the formation. The interaction of interest is Compton scattering, and the detectors may be so adjusted to detect gamma radiation in a range from 100 kev to approximately 600 kev (in the case of a cesium-137 source) to ensure that nearly all of the detected radiation that reaches the detectors is due to the emitted radiation interacting with the environment of the sub by Compton scattering.
It will be appreciated that the radiation emission pattern from the source 24 will be affected by the presence of the sub 12 so that more radiation reaches the formation adjacent the source 24 than is available at the formation on the opposite side of the sub from the source. Consequently, with an asymmetric detector arrangement, such as illustrated in FIG. 3, and FIGS. 4 and 5 in the case of the detectors 26" and 28", the detectors on the same side of the sub 12 as the source 24 will detect more scattered gamma photons than those on the opposite side of the sub. Further, in the case of detectors positioned one above the other, as shown in FIGS. 1-3, a detector that is positioned axially displaced relatively far from the axial position of the source 24 will generally detect radiation scattered from a greater depth within the surrounding formation than will a detector that is positioned closer axially to the level of the source. This is a well known phenomenon, and is due to the fact that radiation penetrating to greater depths in the surrounding formation will generally scatter to greater axial distances from the source. In addition to the foregoing considerations of the geometry of the positions of the detectors relative to the source affecting the radiation sensed by the detectors, the radiation reaching the detectors may also be affected by the position of the sub 12 in the borehole.
Mathematical expressions can be obtained relating the density of the formation in the vicinity of the sub 12, the density of the fluid within the borehole surrounding the sub and the dimensions of the borehole, which expressions can be utilized with the radiation count data obtained by the detectors to provide values for these physical quantities. As an example, the case of right and left side detectors 26 and 28, respectively, of the array of FIGS. 1 and 2 may be considered in relation to FIG. 7 which shows the sub 12 at a level in the borehole wherein the left standoff distance from the sub's outer surface to the borehole wall adjacent the left detector 28 is A, and the right standoff distance from the sub's outer surface to the borehole wall adjacent the right side detector 26 is B. With the outer diameter of the tool, or sub 12, being TD, the width of the borehole at the right and left detectors 26 and 28, respectively, is then the sum of TD plus A and B, which is denoted BD. In a borehole of circular cross section, BD is the borehole diameter or a chord. If measurements are taken while the sub 12 is being rotated about its longitudinal axis, the average borehole transverse dimension may be obtained as BD for any shape borehole cross section and with the sub generally at any lateral configuration within the borehole.
Each of the detectors in general converts the detection of radiation events into an output voltage or current pulse signal. The pulse rate is proportional to the radiation flux at the detector, that is, the number of photons incident on the detector per unit time. The output signal from each detector may be appropriately processed using conventional methods to provide an analog voltage data signal proportional to the radiation count rate sensed by the detector, or a digital signal which gives the count rate. In either case, the data signal can be related to the density of the environment scattering radiation from the source 24 to the detector in question. Thus, for example, the density measured by the left side detector in FIG. 7, .rho..sub.1, is given by the expression:
where:
.rho..sub.m =formation matrix density;
.phi.=porosity of the formation;
.rho..sub.f =density of the fluid in the formation;
.rho..sub.b =density of the fluid (drilling mud) in the borehole; and
.alpha.=fraction of the gamma radiation which interacts in the borehole fluid.
The formation porosity .phi. may be defined in terms of the apparent formation density, or bulk density, of the formation, .rho.', as follows:
which can be written as:
Thus, Equation 1 may be written as:
The formation includes the solid matrix as well as the fluids interspersed in the interstitial space of the matrix. Thus, the density of the entire formation, including that of the solid matrix as well as the fluid, is the apparent formation density, or the bulk density. The density of only the solid portion of the formation is the formation matrix density, as opposed to the density of the fluid only in the formation.
Similarly, the density measured by the right side detector in FIG. 7 may be expressed as:
where .beta. is the fraction of the gamma radiation interacting in the drilling fluid on the right side of the sub.
Equations 1, 4 and 5 are equally valid for wireline applications as well as measurement-while-drilling applications As noted hereinabove, the geometry of the position of the two detectors carried by the tool and located in the borehole is constrained to first order by the aforementioned definition for BD. The probability of a photon traveling a distance A before colliding is represented by e.sup.-k 1.sup.A. The probability, then, of such a collision occurring at a distance less than A is defined as [1-e.sup.-k 1.sup.A ]. Similarly, the probability of a photon undergoing a collision within the distance B is given by [1-e.sup.-k 2.sup.B ]. The constants k.sub.1 and k.sub.2 are geometric constants related to the right and left side detectors, respectively, and are discussed further hereinbelow. For the specified geometry, the constants .alpha. and .beta. may be written as follows:
and
Thus, .alpha. and .beta. are probabilities that the photons will interact in the borehole fluid, and (1-.alpha.) and (1-.beta.) are the probabilities that the photons will interact in the formation and not within the borehole.
Combining Equations 4 and 6 yields:
The left side standoff A may be expressed as:
where k is the sum of the diagonally opposed standoffs A and B, or, with the outer diameter of the sub known, is a measure of the transverse borehole dimension along the line passing through the centers of the left side and right side detectors. Equation 8 can then be expressed in terms of the right side standoff, B, as:
Solving equations 5 and 7 for the exponential term in the standoff B yields:
or
Combining Equations 10 and 12 provides an expression in terms of the measured densities which does not depend on the standoff distances A and B:
Equations 13 can be solved for the apparent, or bulk, formation density:
Equations 13 and 14 are independent of the tool's position in the borehole.
If, for a given pair of detectors, the constants k.sub.1 and k.sub.2 are equal to a constant k.sub.12, Equation 14 reduces to:
This will be the case for two detectors symmetrically disposed relative to the radiation source, such as detector pairs 26 and 28, or 30 and 32, or 30" and 32", but not for detector pairs 26' and 28', or 30' and 32', or 26" and 28", which are asymmetric relative to the source.
Assuming that constant k.sub.12 may be known, and since .rho..sub.l and .rho..sub.r are measured using the two detectors in question, Equation 15 includes three quantities which may be unknown: k, which depends on the transverse dimension of the borehole, .rho..sub.b, which is the density of the fluid in the borehole, and .rho.', the apparent formation, or bulk, density. If the borehole diameter is known and the drilling fluid density is known, Equation 15 may be utilized to determine the bulk density of the formation.
If a second pair of opposed detectors, symmetric relative to the radiation source, is used, and the geometry of the second detectors relative to the radiation source is different so that the radiation measurements taken by the two sets of detectors are not redundant, two independent measurements may be used to determine two of the three unknown parameters present in Equation 15. By providing a third pair of mutually opposed detectors, symmetric relative to the source and with geometry different still from either of the two previous sets of detectors, all three unknown parameters of Equation 15 may be determined unambiguously. For example, Equation 15 may be written for any two detectors positioned mutually opposite each other across the longitudinal axis of the sub 12 and symmetrically with respect to the radiation source. Utilizing the two pairs of detectors arranged as shown in FIG. 2, for example, the following two equations result, wherein the subscripts 26, 28, 30 and 32 are used to identify the density measurements made by the respective detectors and the constant k.sub.12 is characteristic of each of the detectors 26 and 28, and similarly the constant k.sub.34 is characteristic of each of the detectors 30 and 32:
and
In general, the value of k will be the same for all detector pairs if the tool is rotating during measuring of the scattered radiation, or if the tool is centered in a borehole of circular cross section. In the special case of the pairs of detectors being stacked as in FIGS. 1 and 2, the value of k may also be the same for all detector pairs if the transverse dimension of the borehole through the detector pairs is constant. Otherwise, a different value for k would be considered for each equation of the form of Equation 14 or 15, for example.
The two Equations 16 and 17 may be simultaneously solved for two of the three unknown parameters contained therein, again assuming that the constant k.sub.34 may be determined. The addition of a third pair of diametrically opposed detectors, symmetric relative to the radiation source but with geometry relative to the radiation source 24 differing from the geometry of either of the other two pairs of detectors, provides yet a third equation similar to Equations 16 and 17 so that such three equations may be simultaneously solved for three unknown parameters characteristic of the environment in which the sub containing the measuring apparatus is located.
As noted hereinbefore, each detector 26-32 may comprise an array of individual sensors whose output signals may be combined to provide a single output signal from the respective detector taken as a whole. Similarly, the output signals from two detectors may be combined to provide a single output signal, treating the two detectors as one detector. For example, the two right side detectors 26 and 30 illustrated in FIG. 1 may have their output signals combined to provide a single output signal which, after appropriate signal processing, results in a data signal from which may be obtained the number of counts per unit time of photons sensed by the sensors of detectors 26 and 30, so that these two detectors function essentially as a single detector on the right side of the sub 12. In a similar fashion, the two left side detectors 28 and 32 may be operated as a single detector. In each case, the resulting combined detectors are effectively centered longitudinally at a point between the two detectors so combined, so that the center of the right side detector appears midway between detectors 26 and 30 and the center of the left side detector appears midway between detectors 28 and 32.
The resulting combined right side and left side detectors form a pair of detectors which are disposed symmetrically relative to the longitudinal axis of the sub 12 and in a unique geometric configuration with the radiation source, that is, they form a diametrically opposed pair of detectors with their own geometric constant, k.sub.56. With the combined right side and left side detectors centered at a different axial displacement relative to the radiation source 24 compared lo the axial displacement of the centers of the two original pairs of detectors 26, 28 and 30, 32, the pattern of scattered radiation sensed by the right side and left side detectors may be expected to be different from that sensed by either of the original two pairs of detectors, and the geometric constant k.sub.56 of each of the right side and left side detectors is different from k.sub.12 and from k.sub.34. A density measurement equation, similar to Equations 16 and 17, may be written for use with the output signals of the combined right side and left side detectors thus described, where the total density signal for the right side detector is .rho..sub.26+30 and that for the left side detector is .rho..sub.28+32 :
assuming the same value for k as in the case of Equations 16 and 17. Equations 16, 17 and 18 are mutually independent and may be simultaneously solved for the three unknown parameters of the environment surrounding the sub 12.
In cases of opposed detectors which are not positioned symmetrically relative to the radiation source, Equation 14 may be utilized. For example, in using the detector array shown in FIG. 3, Equation 14 may be applied to the lower front and rear detectors 26' and 28', respectively. The same form of equation may be separately applied to the upper front and rear detectors 30' and 32', respectively, and again to an opposed pair of detectors formed by combining the signals from the two front detectors 26' and 30' for one detector, and by combining the signals from the two rear detectors 28' and 32' for the second detector. The result equations, with subscripts relating to the specified detectors, are:
and
again assuming k is the same for all detector pairs.
It will be appreciated that each of the six detectors thus provided by the array of FIG. 3 is situated uniquely relative to the radiation source 24, and therefore has a unique geometric constant. Therefore, each of the three pairs of detectors thus formed features its own geometry relative to the radiation source 24, and, therefore, the three resulting Equations 19-21 are mutually independent and may be solved simultaneously for three unknown parameters reflecting the environment surrounding the sub 12.
Likewise, the output signals from the four detectors of the array illustrated in FIGS. 4 and 5 may be expressed in the form of Equation 14 to obtain the three unknown parameters reflecting the environment surrounding the sub 12. Thus, Equation 14 may be applied to the pair of lower front and rear detectors 26" and 28", and either Equation 14 or Equation 15 may be applied to the upper right and left detectors 30" and 32" (since these detectors are symmetric relative to the radiation source and have the same geometric constant), and Equation 14 may again be applied to the pair of detectors obtained by combining the output signals from detectors 26" and 30" to form one detector, and combining the output signals from the detectors 28" and 32" to form the second detector The resulting equations are
and
where k.sub.x,k.sub.y and k.sub.z are employed as the transverse borehole dimension across the three pairs of detectors 26", 28"; 30", 32"; and 26"+30", 28"+32", respectively.
The detector combination 26" and 30" is symmetric with the detector combination 28" and 32" relative to the longitudinal axis of the sub 12. The three detector pairs thus formed are axially centered at different displacements relative to the axial height of the radiation source 24, and each of the three pairs of detectors comprises two detectors diametrically opposed relative to the longitudinal axis of the sub 12. Further, each such pair of detectors displays a different geometry relative to the position of the radiation source 24, and therefore the three Equations 22-24 may be simultaneously solved for three unknown parameters of the environment surrounding the sub 12, it being understood that k.sub.x, k.sub.y and k.sub.z may constitute up to three unknown parameters.
It will be appreciated that other combinations of detectors may be formed utilizing two pairs of diametrically opposed single detectors. For example, the output signals of detectors 26" and 32" may be combined to effectively form one detector diametrically opposite a second detector formed by combining the output signals of the two detectors 28" and 30". Other detector arrays may be utilized. For example, detector pairs may be placed above and below the axial position of the radiation source. In various configurations of two pairs of diametrically opposed detectors, different combinations of detectors may be formed to provide data signals for use in three independent equations of the form of Equation 14 and/or 15 for simultaneous solution to obtain three unknown parameters of the environment of the measuring apparatus. Care must be taken in combining detectors to form a third pair, without the addition of more sensors, to avoid redundancies in the detector pairs. Such redundancies would be characterized by similar geometry relative to the radiation source for two or more selected detector pairs, and the corresponding two or more equations of the form of Equation 14 and/or 15 would not be mutually independent so that simultaneous solution of the equations would not provide unambiguous values for the parameters sought.
In practice, the density measurements can be made with the detectors while the measuring apparatus 10 is being rotated. Such rotation can be accomplished during normal drilling operations, with the sub 12 being rotated with the drill bit 16. Then, the count rate values obtained from processing the output signals produced by each of the detectors are average count rates of photons detected by detectors looking out from the tool over 360.degree.. With the count rate data obtained during such rotation, the parameters reflecting the environment surrounding the sub 12, obtained by simultaneous solution of Equations 16-18, 19-21 or 22-24, for example, are also average values for the environment surrounding the sub. For example, the value of k=BD-TD yields the average value of the borehole transverse dimension less the tool diameter rather than just the sum of the standoff distances in one direction.
The geometric constants for the individual detectors may be determined utilizing calibrator pit calibrations, but unlike the case for ordinary wireline tool calibrations, the measuring apparatus 10 of the instant invention need not be held against the borehole wall during calibration. To calibrate the measuring apparatus, the tool may be maintained centered in a calibration pit of circular cross section utilizing centralizers, or the tool may be rotated while measurements are being made so that average values are once again obtained regardless of the pit transverse cross section or of the lateral position of the tool in the pit. For both downhole logging as well as calibration purposes, the rotation rate should be at least five rotations per measurement by the detectors.
By knowing the borehole diameter, or average borehole transverse dimension, and the tool diameter, the average standoff of the measuring apparatus in the test pit is determined, and count rates are obtained for the particular detector being calibrated at several depths, with the formation bulk density being known at each depth. For a given such calibrator pit, the logarithm of the count rate obtained by the detector being calibrated is plotted as a function of the bulk density as shown in FIG. 8 (which illustrates representative curves). The detector may be thus calibrated in several pits of differing transverse borehole dimension, providing corresponding different actual non-zero standoffs. Thus, a family of curves as shown in FIG. 8 may be obtained for a given detector for different non-zero standoffs. The slope may be determined in normal fashion for each calibration curve in the family of curves of FIG. 8. The natural logarithms of the calibration curve slopes from FIG. 8 are plotted as functions of the standoff distances to obtain a straight line plot as shown in FIG. 9. From the discussion preceding Equations 6 and 7, it is clear that the slope of the straight line of FIG. 9 is equal to the geometric constant for the detector being calibrated. This geometric constant for the detector is fixed by the geometry of the measuring apparatus, that is, the configuration in the housing, or sub, of the combination of the detector and the radiation source, and needs to be determined only during the original tool design stages. If the type or placement of the detector is changed, the geometric constant will in general be different.
The ordinate intercept of the straight line of FIG. 9 is the slope of the calibration curve for the hypothetical case of zero standoff in FIG. 8. It will be appreciated that, while the zero standoff curve of FIG. 8 may be obtained as a by-product of the measurements used to obtain the geometric configuration constant, k, as described above, the curve may be obtained through direct measurements by immersing the tool in various fluids (of effective infinite extent) and known, differing densities since conformity of the fluid to the tool's surface is equivalent to zero standoff from a formation. Fluids with densities in the range of actual downhole formations may be utilized, in addition to water for example, to obtain the zero standoff curve directly. Then, the slope of the zero standoff curve of FIG. 8 may provide the ordinate intercept of the curve of FIG. 9, which could be constructed to find the geometric constant by providing one more calibration curve, for nonzero standoff, in FIG. 8. The "formations" used in obtaining the zero standoff calibration curve directly are liquids of different known densities. Appropriate conversion of density units, for example between electron density and limestone equivalent density, is well known.
This zero standoff curve, which is plotted in FIG. 8, provides, for a given measured count rate, the corresponding measured average density value of the environment of the detector, or .rho., for the detector in question, considered plotted along the abscissa of FIG. 8. Thus, count rate data obtained by processing the output signal from the detector in question can be translated into average environment density values (e.g. .rho..sub.l or .rho..sub.r in Equations 4 and 5) measured by the detector, by using the zero standoff calibration curve of FIG. 8.
It will be appreciated that the several curves of FIG. 8 intersect at the water point, which corresponds to a count rate obtained for a total environment of water which has a known electron density of 1.1 gm/cc. The tool may be place in a total water environment, and the detector count rate obtained and plotted on FIG. 8 at 1.1 gm/cc density.
The sensitivity of the measuring apparatus 10 may be selectively altered by changing the geometry of the detectors, that is, their placement relative to the radiation source 24, to change their geometric constants. For example, the variation of sensitivity of measuring apparatus depending on detector geometric constant is illustrated in the graph of FIG. 10 for a tool diameter of 7.0 inches in a 9.0 inch diameter borehole. Three different geometric k constants are considered, the geometric constants being 0.35, 0.65 and 1.35 corresponding to the detector being positioned progressively closer to the radiation source, respectively. A formation bulk density of 2.6 gm/cc surrounds the measuring apparatus, and mud of density 1.2 gm/cc surrounds the measuring apparatus in the borehole. The curves for the three different geometric k constants intersect at the point representing these two densities. Different values for the mud density were assumed, and the corresponding apparent density values calculated utilizing an equation of the form of Equation 14. The resulting apparent density values were plotted as a function of the assumed mud density values for each detector geometric constant.
The resulting graph of FIG. 10 shows the size in errors in apparent density values that would be obtained if incorrect mud density values were assumed, for example. Additionally, the graph demonstrates how the measuring apparatus may be designed to select the sensitivity of the tool. For example, detectors positioned relatively close to the radiation source, and therefore having higher geometric constants, will have a stronger dependence on the mud density in the borehole than detectors spaced farther from the geometric source, that is, particularly in axial displacement from the source, and therefore having lower geometric k constants. Thus, measuring apparatus may be constructed that is relatively insensitive to the density of the borehole fluid, or, conversely, measuring apparatus may be constructed which enhances the dependence of the output data on the mud density, whereby the measuring apparatus may be used to measure the borehole fluid density more accurately.
FIG. 11 provides a graph which is similar to that of FIG. 10, except that in FIG. 11 the sensitivity of the measuring apparatus in determining borehole diameters, as a function of the detector geometric constant, is demonstrated. Again, a 9.0 inch diameter borehole is indicated, the bulk density of the formation is given as 2.60 gm/cc, and the borehole fluid density is 1.20 gm/cc. The same detector geometric constants demonstrated in FIG. 10 are utilized for the presentation of FIG. 11. All three plotted curves intersected at the point representing the aforementioned bulk and mud density values. Different borehole diameters were then assumed, and the apparent matrix density was calculated using an equation of the form of Equation 14. The resulting curves once again demonstrate that a detector positioned closer to the radiation source, and therefore having a larger geometric constant, provides a greater sensitivity in the ability to determine the transverse dimension of the borehole.
FIG. 12 further demonstrates the dependence of the sensitivity of the measuring apparatus on the detector geometric constant for borehole diameter measurements. Data may be obtained from the graph of FIG. 11 for any given geometric constant, and additional graphs similar to that of FIG. 11 for the same geometric constant may be obtained (for the same or different borehole diameter) utilizing different values for the actual apparent formation density. For a given true apparent density (2.60 gm/cc in FIG. 11), different values for the borehole diameter are assumed, and the corresponding apparent density values obtainable from FIG. 11 or like graphs are noted, and plotted along the vertical axis of FIG. 12 as observed apparent density. For each selected borehole value, a straight line is obtained in FIG. 12 as illustrated. Since the correct borehole diameter is 9.0 inch, the straight line for the 9.0 inch borehole data has a slope of 1. It will be appreciated by reference to FIG. 12 that a very sensitive caliper may be obtained with the measuring apparatus of the present invention to determine borehole diameters, at least with an accuracy of 0.1 inch, for example.
The foregoing discussion illustrates that, in general, the sensitivity of the measuring apparatus may be enhanced by placing the detectors closer to the radiation source. However, the distance which radiation from the source travels into the formation and is scattered back to the detectors increases as the detectors are moved axially farther from the source. This depth of penetration for scattered and detected radiation ultimately peaks as the detector is moved farther away from the radiation source along the axis of the tool. Consequently, to obtain information concerning the formation at significant depths into the formation from the borehole, and to obtain desired sensitivity of the measuring apparatus, the axial position of the detectors relative to the radiation source may be selectively adjusted to optimize the two measuring conditions, namely sensitivity and formation penetration.
As discussed herein, the present measuring apparatus need not be held against the borehole wall to obtain any of the measurements noted. Thus, inaccuracies due to a tool following a major elliptical axis, for example, of a borehole cross section are avoided. Further, while the measuring apparatus may be utilized in a stationary mode, the benefits of data averaging may be appreciated by operating the detectors while the measuring apparatus is being rotated, for example at 100 revolutions per minute.
With the radiation source positioned on the side of the tool sub as illustrated, it will be appreciated that, if the sub is against the wall of the borehole at the level of the radiation source, misleading data counts may be obtained. For example, if the radiation source is facing away from the borehole wall with the sub against the borehole wall, an overabundance of the radiation scatter occurs in the borehole fluid, and the detector counts may be misleadingly high. If the radiation source is pressed against the borehole wall, virtually all of the radiation may enter the formation with little or no initial scattering in the borehole fluid, and the radiation will interact directly with the formation before having had any opportunity to travel axially in the borehole so that the data counts may be misleadingly low. Rotating the measuring apparatus while it is in operation to obtain averaging effects will alleviate such misleading count rate problems. A subgauge stabilizer may be used about the sub 12 to keep it from being positioned too close to the borehole wall to avoid such problems as well. The subgauge stabilizer, which may also be used even with the tool being rotated during measuring operations, would not interfere with steering of the drill bit, for example.
The present invention does not require the use of collimation to define radiation beams from the source or to the detectors, since there is no such dependence on direction built into Equation 14, for example. If a direction sensor is utilized in conjunction with the present measuring apparatus, and the measuring apparatus is not rotated during measurement operations, then details of the formation and borehole profile as functions of direction may be observed, even without collimation of the radiation.
The present invention thus provides a technique for simultaneous measurement of multiple downhole parameters by utilizing two or more pairs of diametrically opposed detectors, and may be employed in any appropriate fashion, such as in a measuring-while-drilling mode or suspended by a wireline, for example. The capability of making all measurements simultaneously to provide sufficient data to determine multiple downhole parameters avoids the inaccuracies that accompany measurements made after delays have allowed conditions in the borehole to change (mud caking, etc.), and where measurements are made during the drilling process or during wireline logging, avoids unduly axially extending the region in which measurements are being made, so that measurements to obtain all parameters are made at the same location.
The foregoing disclosure and description of the invention is illustrative and explanatory thereof, and various changes in the method steps as well as in the details of the illustrated apparatus may be made within the scope of the appended claims without departing from the spirit of the invention.
Claims
What is claimed is:
1. A method of calibrating a radiation detector positioned in apparatus including a source of gamma radiation such that the detector may sense radiation from the environment of the apparatus due to radiation emitted by the source and scattered by the environment, comprising the following steps: a. providing a test environment including a formation of known bulk density; b. positioning the apparatus in the test environment with the detector at a first distance from the surface of the formation; c. sensing, with the detector, radiation scattered from the environment due to radiation emitted by the source of gamma radiation; d. comparing the count rate of the detector to the bulk density of the formation to obtain a first data point: e. repeating steps a. through d., for a formation of different known bulk density and for the same distance between the detector and the formation, to obtain a second such data point; f. repeating the steps a. through d. utilizing a formation of known bulk density and a second distance between the detector and such formation to obtain a third such data point; g. repeating step f. for a formation of different known bulk density and the same second distance between the detector and the surface of the formation to obtain a fourth such data point; h. obtaining the slope of the functional relationship between the logarithm of the count rate of the detector as a function of the known densities of the formations utilizing the first and second data points, and likewise obtaining a second such slope utilizing the third and fourth data points; i. comparing the logarithm of the first and second slopes as a function of the corresponding distances between the detector and the surfaces of the formations to obtain a value for the slope for an effective zero distance between the detector and the surfaces of formations; and j. utilizing the slope for the effective zero distance between detector and formation surfaces as the functional relationship between the logarithm of the detector count rate and the density of the environment of the detector as a calibration function.
2. A method as defined in claim 1 wherein the slope of the relationship between the logarithm of the slopes of the relationships obtained by utilizing the data points, as a function of the distances between the detectors and the surfaces of the formations, provides a geometric constant indicative of the configuration of the detector relative to the source of gamma radiation.
3. A method as defined in claim 1 further comprising repeating steps a. through d. for additional formations of different known densities to obtain additional data points for any of the distances between the detector and the surfaces of the formations used.
4. A method as defined in claim 1 further comprising repeating steps a. through e. additional times for different known distances between the detector and the surfaces of the formations.
5. A method as defined in claim 1 wherein the environment is provided by means of a borehole passing through a formation of known bulk density.
6. A method as defined in claim 1 wherein the formation used to scatter radiation to obtain at least one data point is essentially water.
7. A method as defined in claim 6 wherein the distance between the water and the detector is effectively zero.
8. A method of calibrating a radiation detector positioned in apparatus including a source of gamma radiation such that the detector may sense radiation from the environment of the apparatus due to radiation emitted by the source and scattered by the environment, comprising the following steps: a. providing a test environment including a formation of known bulk density; b. positioning the apparatus in the test environment with the detector at a first distance from the surface of the formation; c. sensing, with the detector, radiation scattered from the environment due to radiation emitted by the source of gamma radiation; d. comparing the count rate of the detector to the bulk density of the formation to obtain a first data point; e. repeating steps a. through d., for a formation of different known bulk density and for the same distance between the detector and the formation, to obtain a second such data point; f. repeating the steps a. through d. utilizing a formation of known bulk density and a second distance between the detector and such formation to obtain a third such data point; g. repeating step f., if necessary, for a formation of different known bulk density and the same second distance between the detector and the surface of the formation to obtain a fourth such data point; h. obtaining the slope of the functional relationship between the logarithm of the count rate of the detector as a function of the known densities of the formations utilizing the first and second data points, and likewise obtaining a second such slope utilizing the third and fourth, if necessary, data points; and i. comparing the logarithm of the first and second slopes as a function of the corresponding distances between the detector and the surfaces of the formations to obtain a geometric constant indicative of the configuration of the detector relative to the source of gamma radiation.
9. A method as defined in claim 8 further comprising repeating steps a. through d. for additional formations of different known densities to obtain additional data points for any of the distances between the detector and the surfaces of the formations used.
10. A method as defined in claim 8 further comprising repeating steps a. through e. additional times for different known distances between the detector and the surfaces of the formations.
11. A method as defined in claim 8 wherein the environment is provided by means of a borehole passing through a formation of known bulk density.
12. A method as defined in claim 8 wherein the formation used to scatter radiation to obtain at least one data point is essentially water.
13. A method as defined in claim 12 wherein the distance between the water and the detector is effectively zero.
14. A method of calibrating a radiation detector positioned in apparatus including a source of gamma radiation such that the detector may sense radiation from the environment of the apparatus due to radiation emitted by the source and scattered by the environment, comprising the following steps: a. providing a test environment including a liquid of known density; b. positioning the apparatus in the test environment with the liquid effectively at the detector; c. sensing, with the detector, radiation scattered from the liquid due to radiation emitted by the source of gamma radiation; d. comparing the count rate of the detector to the density of the liquid to obtain a first data point; e. repeating steps a. through d., for a liquid of different known density, to obtain a second such data point; and f. utilizing the slope of the functional relationship between the logarithm of the count rate of the detector as a function of the known densities of the liquids as a calibration function.
15. A method as defined in claim 14 further comprising repeating steps a. through d. for additional liquids of different known densities to obtain additional data points.
16. A method as defined in claim 14 further comprising the additional steps of: g. carrying out the steps a. through e. of claim 14 for at least one formation of different known bulk density, spaced a known distance from the detector to obtain at least a third data point; h. obtaining the slope of the functional relationship between the logarithm of the count rate of the detector as a function of the known density of the formation; and i. utilizing the logarithms of the slope of the calibration function and of the slope of step h. compared with the distance between the formation and the detector to obtain a geometric constant indicative of the configuration of the detector relative to the source of gamma radiation on the apparatus.