# Patent

## Title
Neutron through-pipe measurement, device, system and use thereof

## Abstract

A method for evaluating a characteristic of a formation adjacent a wellbore comprises emitting a neutron field into the formation from a neutron source, obtaining data comprising a long neutron neutron (LNN), short neutron neutron (SNN), long neutron gamma (LNG) and short neutron gamma data (SNG) from the neutron field emitted into the formation and combining the long data with the short data to reduce error associated with the characteristic to be evaluated. A detector for evaluating a characteristic of a formation includes a neutron source and detectors for detecting LNN, SNN, LNG, SNG.

## Background
FIELD OF THE INVENTION
The invention relates to systems, devices and methods for taking through-pipe measurements of a formation, and to characterization of the formation based on those measurements.
BACKGROUND
Borehole measurement devices may be used to determine formation characteristics surrounding the borehole and are typically used in wellbores drilled for the purpose of extracting natural resources, such as hydrocarbon, from the formations surrounding the borehole. Borehole measurement devices may use different types of measurements, for example, a borehole measurement device may use gamma measurements, thermal neutron measurements, resistivity measurements or other types of measurements.
At present there are a number of physics being employed to perform thru pipe formation evaluation and other wellbore measurements. For neutron based measurements, there are pulsed neutron thermal neutron, pulsed neutron gamma, gamma gamma, neutron thermal neutron, neutron epi-thermal neutron etc. It is noted that naming convention for down hole geophysical devices are based on source—detection physics. For example neutron thermal neutron indicates a neutron source and thermal neutron detection. Most of these systems utilize one source per physics, this will be referred to as single physics measurement. Systems that utilize one source for two physics will be referred to as a dual physics measurement. Of the currently available dual physics systems, there is one that utilizes a combination neutron thermal neutron, neutron epi-thermal neutron. In addition there is a neutron thermal neutron, neutron gamma dual physics measurement device built in Azerbaijan.
Some shortcomings or disadvantages of using single physics measurements is the lack of corrections available for factors such as borehole rugosity, annular fluid changes, mineralogy, tubulars etc. To compensate for these short comings, dual detectors are used. These devices are commonly referred to as compensated devices. Another solution is to combine multiple single physics measurement devices, including compensated devices, during analysis. An example is neutron thermal neutron and gamma gamma physics, commonly referred to as neutron density measurements. Measuring thru pipe also limits the effectiveness of some of the physics. For example, gamma gamma measurements are limited because the pipe itself shields gamma and therefore there are losses as the gamma photons travel from the source thru the pipe and then again as the photon returns back to the detector. This results in low count rates and increases error in the measurements.
Neutrons can easily penetrate pipe and therefore is a logical choice for thru pipe measurements. Of the neutron physics based measurements; pulsed neutron devices do not lend themselves to dual physics measurements. The reason is that the length of the pulsed neutron source does not allow for effective measurement spacing for the required detectors. Chemical neutron sources are much smaller and therefore can be used effectively for dual physics based measurements. The neutron epi-thermal measurement is highly sensitive to borehole rugosity and therefore not an ideal choice to determine formation parameters.
A measuring tool or logging device together with a method of carrying out formation evaluation using same is therefore needed.
SUMMARY OF THE INVENTION
Methods, systems and logging tools are provided for evaluating characteristics of a formation or formations adjacent a wellbore, down which a logging tool is positioned. The characteristics are evaluated based on the detection of a neutron thermal neutron field and a neutron gamma field generated by a neutron source emitting into the formation and returning to detectors, generally four detectors, positioned in an orientation relative the neutron source to reduce error associated to the detected fields. This may be accomplished through a combining formula as it was determined that the errors associated to the neutron gamma field are in the opposite direction to the errors associated to the neutron thermal neutron field.
The logging tool may be used in wellbores where a casing has been put in place as measurements may be taken through the casing and the various characteristics of the formation may be evaluated thereby allowing for evaluation after installation of the casing and as a result, monitoring of the formation at various stages throughout production.
In one embodiment, the present invention provides for a method of evaluating a characteristic of a formation adjacent a wellbore comprising:

    
    
        emitting a neutron field into the formation from a neutron source;
        obtaining data comprising a long neutron neutron (LNN), short neutron neutron (SNN), long neutron gamma (LNG) and short neutron gamma data (SNG) from the neutron field emitted into the formation;
        combining the long data with the short data to reduce error associated with the characteristic to be evaluated.
    
    


In a further embodiment of the method above, the method further comprises the step of:

    
    
        optimizing a spacing of a long neutron detector for detecting LNN and a short neutron detector for detecting SNN relative the neutron source to minimize the error associated to the characteristic to be evaluated.
    
    


In a further embodiment of the method above, the method further comprises the step of:

    
    
        optimizing a spacing of a long gamma detector for detecting LNG and a short gamma detector for detecting SNG relative the neutron source to minimize the error associated to the characteristic to be evaluated.
    
    


In a further embodiment of the method above the data is obtained through a wellbore casing.
In a further embodiment of the method above the combining step comprises the cancelation of the error associated to the long data with the error associated to the short data.
In a further embodiment of the method above the characteristic is neutron porosity (QTP) and the formula for evaluating the characteristic QTP is:




  
    Q
    ⁢
    
      
    
    ⁢
    T
    ⁢
    
      
    
    ⁢
    P
  
  =
  
    
      A
      *
      
        ln
        ⁡
        
          (
          
            
              SNN
              *
              LNG
            
            
              SNG
              *
              LNN
            
          
          )
        
      
    
    +
    B
  




wherein A and B are determined empirically and are influenced by bit size and the amount of high dense material in a volume of investigation.

In a further embodiment of the method above the characteristic is neutron clay (QNC) and the formula for evaluating the characteristic QNC is:




  QNC
  =
  
    
      A
      *
      
        (
        
          
            C
            
              (
              
                SNG
                LNG
              
              )
            
          
          +
          
            D
            *
            
              (
              
                SNN
                LNN
              
              )
            
          
        
        )
      
    
    +
    B
  




wherein A, B, C and D are empirically determined coefficients.

In a further embodiment of the method above A, B, C and D are 0.004, 0, 1.9 and −1.5, respectively.
In a further embodiment of the method above the characteristic is neutron liquid (QNL) and the formula for evaluating the characteristic QNL is:




  QNL
  =
  
    
      A
      *
      
        ln
        ⁡
        
          (
          
            
              LNG
              *
              SNG
            
            
              LNN
              *
              SNN
            
          
          )
        
      
    
    +
    B
  




wherein A and B are chosen to create the best overlay of QNL to QTP.

In a further embodiment of the method above the characteristic is water saturation (Sw) and the formula for evaluating the characteristic Sw is:




  Sw
  =
  
    1
    -
    
      kf
      ⁢
      
        
          (
          
            QEP
            -
            QEL
          
          )
        
        QEP
      
    
  




wherein

    
    
        k=fluid factor;
        f=formation factor;
        QEP=Quad Effective Porosity (Clay Free); and
        QEL=Quad Effective Liquid porosity (Clay Free); and
        wherein
        QEP and QEL are determined as follows:

QEP=QTP−QC and QEL=QL−QC
        
            wherein
            QTP=Quad Total Porosity;
            QC=Quad Clay Porosity; and
            QL=Quad Liquid Porosity;

the fluid factor k=1/(MaxPor*QLgain)

            wherein
            MaxPor=WaterPor−OilPor;
            WaterPor is determined as follows:

WaterPor=(12.73+√{square root over (12.732−4*0.0966*(403.04−Waterkppm)))}/(2*0.0966)

        
        
        OilPor is determined as follows:

OilPor=0.1333*oilAPI+71 and

        QLgain is computed as

QLgain=QNL A value/32.5.

    
    


In another embodiment, the present invention provides for a method of determining the relative bulk density of a formation adjacent a wellbore through the wellbore casing comprising:

    
    
        emitting a neutron field into the formation from a neutron source;
        obtaining long neutron neutron (LNN) field readings and long neutron gamma (LNG) field readings;
        converting the LNN and LNG for porosity to provide LNNpor and LNG por respectively; and
        determining a relative bulk density by subtracting LNGpor from LNNpor,
        wherein LNNpor and LNGpor represent the porosity converted counts.
    
    


In another embodiment, the present invention provides for a method of determining the nuclear caliper of a formation adjacent a wellbore through the wellbore casing comprising:

    
    
        emitting a neutron field into the formation from a neutron source;
        obtaining long neutron neutron (LNN) field readings and short neutron neutron (SNN) field readings;
        converting the LNN and SNN for porosity to provide LNNpor and SNN por respectively; and
        determining a nuclear caliper by subtracting LNNpor from SNNpor,
        wherein SNNpor and LNNpor represent the porosity converted counts.
    
    


In another embodiment, the present invention provides for a method of determining a chemical effect of a formation adjacent a wellbore through the wellbore casing comprising:

    
    
        emitting a neutron field into the formation from a neutron source;
        obtaining long neutron gamma (LNG) field readings and short neutron neutron (SNN) field readings;
    
    


converting the LNG and SNN for porosity to provide LNGpor and SNN por respectively; and
determining the chemical effect by subtracting the LNGpor from the SNNpor.


    
    
        wherein SNNpor and LNGpor represent the porosity converted single detector counts.
    
    


In another embodiment, the present invention provides for a system for evaluating a characteristic of a formation adjacent a wellbore comprising:

    
    
        a logging tool for downhole use collecting data for evaluating a characteristic of a formation adjacent a wellbore, the logging tool comprising:
        
            a neutron source for emitting a neutron field into the formation;
            a short detector for detecting a neutron gamma field and a neutron thermal neutron field; and
            
                a long detector for detecting a neutron gamma field and a neutron thermal neutron field; and
            
            
        
        
        a computer for combining the long data with the short data to reduce error associated to the characteristic to be evaluated
    
    


In another embodiment, the present invention provides for a logging tool for downhole use collecting data for evaluating a characteristic of a formation adjacent a wellbore, the logging tool comprising:

    
    
        a neutron source for emitting a neutron field into the formation;
        a short detector for detecting a neutron gamma field and a neutron thermal neutron field; and
        a long detector for detecting a neutron gamma field and a neutron thermal neutron field.
    
    


In a further embodiment of the logging tool above the short detector and the long detector are situated at an optimal distance from the neutron source so that the error associated to the short detector substantially cancels the error associated to the long detector.
In a further embodiment of the logging tool above, the logging tool further comprises a housing comprising beryllium copper.



BRIEF DESCRIPTION OF THE DRAWINGS
Embodiments are described herein with references to the appended drawings, in which:
FIG. 1 is a diagram of a single source neutron logging tool for carrying out neutron detection in a wellbore, optionally thru a casing, of formation characteristics adjacent or nearly adjacent the wellbore;
FIG. 2 is a diagram of a drilling system for taking neutron measurements using a single neutron source;
FIG. 3 is a diagram of an example of a single neutron source logging tool having both long and short detectors capable of taking measurements thru a wellbore casing;
FIG. 4 depicts a further example of a single source logging tool having both long and short detectors capable of taking measurements thru a wellbore casing; and
FIG. 5 is a graph showing test data illustrating the slopes of the natural log of the fields of both neutron thermal neutron and neutron gamma;
FIG. 6 is a graph illustrating the effect of casing on both neutron thermal neutron and neutron gamma fields;
FIG. 7 is a graph illustrating the relationship on QTP gain (A) vs. various casing weights in characterizing quad neutron porosity;
FIG. 8 is a graph showing test data illustrating the effect of varying water salinity to 32 API oil in characterizing quad neutron liquid;
FIG. 9 is a plot showing test data illustrating a comparison of quad neutron clay free porosity curve plotted against core porosity data along with open hole neutron density porosity on both sandstone and limestone matrix done through-casing;
FIG. 10 is a plot of data collected from the same well as FIG. 5, of the core clay volume plotted against the quad neutron clay curve done through-casing;
FIG. 11 is a plot of test data showing the quad Relative Bulk Density curve (labelled as CE) compared to the open hole bulk density curve (RHOB) and the bulk density obtained from core analysis;
FIG. 12 is a plot of test data showing current in situ saturation done through-casing;
FIG. 13 is a data set collected from a horizontal well initially drilled and left standing because of poor to no oil shows during drilling. Tubing was run in the well with a perforated pup joint on the bottom. A memory quad neutron was pumped down on slickline to the end of the tubing and then logged out with slickline;
FIG. 14 is a plot of the well from FIG. 13 showing the production results showing that the quad neutron correctly identified the region of highest oil saturation wherein the measurements were taken through-casing; and
FIG. 15 is a plot showing the natural log of the detector counts against detector spacing.



DETAILED DESCRIPTION
The use of a borehole logging tool or measurement device in determining formation characteristics of an existing vertical borehole is shown for illustrative purpose in FIG. 1. As depicted, a borehole 102 is drilled through a plurality of formations 104a, 104b, 104c, 104d, 104e. A logging tool 106 or measurement device may be lowered down the borehole to determine one or more characteristics within the formation(s). The logging tool 106 may include, for example, a neutron source 108 and a plurality of neutron detectors including short neutron detectors 110a and long neutron detectors 110b, such as gamma detectors, thermal neutron detectors, or epithermal neutron detectors. The positioning of the neutron detectors will be discussed in further detail below and with reference to FIGS. 3 and 4. The measurement characteristics may be associated with the height of the measurement within the borehole. Measurements in vertical, horizontal or deviated boreholes may be made after the borehole is drilled or as it is being drilled and may also be made after installation of the casing.
The neutron analysis technology using four different types of data set collection including two data sets for neutron thermal neutron field detectors (one closer to the neutron source referred to as the short neutron thermal neutron field or data set and one further from the neutron source referred to as the long neutron thermal neutron field or data set) and two data sets from neutron gamma field detectors (one closer to the neutron source referred to as the short neutron gamma field or data set and one further from the neutron source referred to as the long neutron gamma field or data set) is referred to herein is as quad neutron technology as the neutron analysis technology includes four types of data set and the logging tool which both emits and detects the neutron field includes, in this case a single neutron source, and two sets of two detectors, namely short thermal neutron and gamma detectors and long thermal neutron and gamma detectors. An illustrative embodiment of a quad single neutron source logging tool and drilling system will be discussed in more detail below with reference to FIGS. 2, 3 and 4.
A single neutron source with multiple detectors can provide useful information on formation characteristics. As described further herein, it is possible to obtain formation characteristics from neutron field measurements acquired by neutron detectors including thermal neutron detectors and gamma detectors positioned in the logging tool at short and long distances from the neutron source.
FIG. 2 is a diagram of a drilling system for taking formation measurements using a single neutron source logging tool. As depicted, a borehole 202 is drilled through the ground, which may comprise numerous different formations 204a, 204b, 204c (referred to collectively as formations 204). As will be appreciated, different formations will have different characteristics such as porosity, density, clay levels, oil saturation, water saturation and other characteristics.
The borehole 202 is formed by a drilling system 206.
The drilling system 206 comprises a number of components for drilling the borehole 202 as well as measuring and logging data of the borehole 202 and surrounding formations 204. As will be appreciated, the drilling system 206 is representative of a typical drilling system and numerous specific drilling systems are possible, with their particular suitability to particular applications known to one skilled in the art. The drilling system 206 comprises a drilling string in the borehole. The drilling string has at its head a drill bit 208 and a mud motor 210. Drilling mud is supplied to the mud motor 210 from the surface under pressure. The drilling mud causes the mud motor 210 to rotate the drill bit 208. It is noted that when using a mud motor as depicted, up-string components, that is components of the drilling string connected between the mud motor 210 and the surface do not rotate with the drill bit 208.
One or more pipes 212 may be included in the drill string. The pipe 212 is depicted as having a number of sensor devices 216, 218, 226 within its interior 214. As will be appreciated, numerous different sensor devices can be incorporated into the drill string at various different locations. The sensor devices 216, 218, 226 are depicted as being housed within the same pipe 212 for convenience. The different sensor devices may provide various functionality. For example, sensor device 216 may provide telemetry measurements, inclination measurements, magnetic field strength measurements, pressure measurements, vibration measurements, or other types of measurements as desired or appropriate as would be appreciated.
The sensor device 218 also provides measurements, however as described further below, it provides neutron-based measurements. These measurements can provide an indication of formation characteristics surrounding the borehole. The formation characteristics that can be determined from the measurements depend upon the analysis carried out following collection of the field information by the short and long sensors incorporated into the device, and may include measurements such as density, porosity, saturation, borehole caliper, etc. The sensor device 218 is depicted as having a neutron source 222 and a long neutron detector 220 and a short neutron detector 224. The short and long neutron detectors may comprise a gamma detector and a thermal neutron detector.
In one embodiment, the quad logging tool comprises both thermal neutron detectors on one side of the neutron source and the neutron gamma detectors on the other side. The layout shown in the Figures shows the short neutron on the opposite side of the source as compared to long neutron which is not intended to be limiting as, outlined above, the thermal neutron detectors may be on one side of the neutron source and the neutron gamma detectors on the other side.
The sensor devices 216, 218 may store the measurement information internally in memory for later retrieval and processing at the surface. Additionally or alternatively, the measurement information, or a portion of the measurement information, may be transmitted in real-time or near real-time to the surface using various communication techniques. For example, an additional device 226 may receive the measurement information from a plurality of sensor devices 216, 218 and provide storage of the measurements and/or communicate the measurements to the surface. For example, the additional device 226 can provide a wire line communication interface or other type of down-hole communication interface such as mud-pulse communication. Although the storage and/or communication of measurements by sensor devices 216, 218 have been described as being provided by an additional device 226, it is contemplated that the storage and/or communication interface may be provided in each individual device.
FIG. 3 is a diagram of a logging tool capable of taking neutron field measurements using a single neutron source. The logging tool 900 may be incorporated into a drill string in order to provide measurements while drilling or alternatively may simply be put down the wellbore following drilling to provide a characterization of the formations adjacent the wellbore at various stages of operation or production. The logging tool 900 comprises a housing 902 or casing. The housing may be made from beryllium copper which has been determined to have very good gamma transmissivity properties while still retaining the required mechanical properties for down hole tool housings, such as for example resistance to pressure, temperature and corrosion.
The logging tool 900 further includes a neutron source 908 for emitting a neutron field into the formation(s) adjacent the wellbore and neutron-based detectors 910,912, 914, 916, such as gamma detectors, thermal neutron detectors and/or epi-thermal neutron detectors. The neutron detectors are situated from the neutron source 908 so that certain detectors are closer to the neutron source, referred to as short detectors such as detectors 912 and 914 and certain detectors are further way from the neutron source, referred to as long detectors such as detectors 910 and 916. As depicted, the neutron source 908 and neutron-based detectors 910, 912, 914, 916 are arranged off-centre from the logging tools centre axis, however, any suitable orientation may be used so long as a short and long measurement may be made. Further discussion relating to the orientation of the short and long detectors will be outlined in more detail below.
FIG. 4 depicts a further logging tool capable of taking measurements for characterising a formation adjacent a wellbore using a neutron source. The logging tool 1000 may be received within a drill pipe or a casing. The logging tool 1000 may be rotated with the drill pipe, which in turn is rotated by the drilling rig. Alternatively, an electric motor may be used to rotate the logging tool if power can be supplied to the logging tool from the surface. It will be appreciated that other techniques for rotating the logging tool may be employed.
The logging tool 1000 may comprise a cylindrical housing 1002 for housing the various detectors and components of the tool. Although the particular dimensions may vary, the logging tool 1000 is depicted as having a length of 179 inches. The main cylindrical housing may have an outer diameter of approximately 1.7 inches and may be made from beryllium copper which has been determined to have very good gamma transmissivity properties while still retaining the required mechanical properties for down hole tool housings, such as for example resistance to pressure, temperature and corrosion. As depicted, the logging tool 1000 includes various sensors and detectors located along the length of the tool. An up-hole temperature sensor 1004, which may also measure fluid resistivity, may be located 6.25 inches from an up-hole end of the tool. A high gamma ray detector 1006 and a gamma ray detector 1008 may be located 17.5 inches and 38.1 inches from the up-hole end, respectively. A casing collar locator 1010 may be located at 71.5 inches from the up-hole end. Long and short neutron gamma detectors 1012, 1014 may be located 101 inches and 110 inches from the up hole end, respectively. A neutron source 1016, such as an Americium Beryllium source, may be located at 121.7 inches from the up-hole end. Short and long neutron neutron detectors 1020, 1022 may be located at 136.2 inches and 144.4 inches from the up-hole end, respectively. A bottom hole temperature and fluid resistivity sensor 1024 may be located 172.8 inches from the up hole end of the logging tool. It should be appreciated that the logging tool depicted in FIG. 4 is intended to provide one illustrative arrangement of components. It is contemplated that various components can be omitted, their positions within the tool changed including the lengths from the up-hole end and other components included.
Due to the statistical nature of neutron measurements, a dual measurement physics system has been developed to reduce statistical differences due to variations in neutron output from separate neutron sources.
To create a balanced measurement of the fields, both neutron thermal neutron and neutron gamma fields were studied. It was determined that with proper detector spacing between the short detectors and the source and the long detectors and the source the slopes of the natural log of the measurements could be made to match in different fluid types. Test data illustrating this is shown FIG. 5.
Research done on the effects of taking measurements through the wellbore casing were then studied. The observation was that the neutron thermal neutron slopes remained unchanged with different casings and that the neutron gamma slopes started to show deviations in slope. The other observation was that the effect of the count change on the neutron gamma detectors was minimal if the correct detector spacing was chosen.
The ideal detector spacing varies based on the measurement environment. Since the measurement environment is not consistent, the detector spacing should be chosen based on the best average fit.
The effect of casing on both neutron thermal neutron and neutron gamma is shown in FIG. 6.
Since the effect of shielding was noticeable on the neutron gamma measurements, research was performed on the different logging tool housing materials in an effort to develop a housing that reduced or minimized the effect of the housing on the neutron gamma measurements. It was determined that beryllium copper had very good gamma transmissivity properties while still retaining the required mechanical properties for down hole tool housings, such as for example resistance to pressure, temperature and corrosion.
It was found that certain materials used in the photomultipliers on the neutron gamma sensors may cause artificial counts when subjected to a neutron flux. As a result, Hamamatsu R1288A photomultipliers may be used. This particular detector responds predictably and un-dramatically to exposure to near neutron source flux.
Quad Neutron Physics Model
Upon completion of the original single source quad neutron logging tool, testing was performed in a known environment to confirm the theoretical results. Results were not as expected. Testing revealed that the quad neutron responses did not meet expectations based on publicly available theoretical material related to the measurement physics. The entire design was reviewed and errors were found in the response function of the semiconductors used in the neutron logging tool. The logging tools were then redesigned with military grade semiconductor components. For a second time the quad neutron responses did not meet theoretical expectations during experimental testing. A second design review was conducted and no errors were found with the measurement device. The conclusion was that the common theoretical models used were incorrect for the neutron logging tool measurement device. As a result a bespoke measurement model for the quad neutron characterization of formations using data collected by a quad neutron logging tool was defined.
Research was conducted to determine the properties for the neutron measurement model of the quad neutron logging tool. Observations during this research led to the conclusion that the nuclear properties for the quad neutron logging tool may be characterized with the physics of particles and fields. This realization led to the development of the Quad NN (Neutron Neutron) String model.
Without wishing to be bound by theory or limited in any way, the non-limiting quad NN string model may be explained as follows. In any matrix, a neutron has a maximum travel distance from the neutron source before it is captured. The travel length from the source to the furthest capture point is considered a straight line which can be replaced by a string with one point permanently fixed at the source and the other at the capture point. The string then represents the maximum distance a neutron can travel from the source within the matrix. If a neutron can travel the entire length of the straight string, in should then be able to travel half the distance and then back again to the source. This can be represented by fixing both ends of the string at the source. Since detectable neutrons are of primary concern, one end of the string is left at the neutron source and the other end of the string is placed at the detector. If all the possible paths that the fixed string will allow are traced, this represents all the paths that detectable neutrons will travel and a volume of investigation is defined. The volume of investigation is the volume from which the quad neutron logging tool gathers data. As such, all materials lying within that volume contribute to the recorded response.
The shape of the volume of investigation is ellipsoidal with the source and detector as the foci of the major axis. The minor axis represents the depth of investigation and is dominated by the thermalization and neutron capture properties within the volume of investigation. The length of the string is directly related to the thermalization and thermal neutron cross section capture area properties of the matrix. The thermalization properties of a matrix will be referred to as the Thermalization index and the thermal neutron cross section capture area properties as the Capture index.
From a theoretical perspective, the NN string model was defined but the quad neutron logging tool also incorporated a neutron gamma (NG) component. Research was conducted to determine the properties for the neutron gamma measurement model of the quad neutron logging tool. Originally the conclusion was that the NG model should be spherical in nature as the source of the gamma radiation originated as a result of the capture event of the neutron. Testing revealed that the gamma field is similar in shape to the neutron field that generated it and comparisons between neutron neutron and neutron gamma counts confirm this relationship. Since the measurement is truly a gamma measurement it was also known that the gamma photons can be stopped by increasing the density of the matrix it travels through. This relationship then defines an elliptical volume of investigation with the foci lengths directly proportionate to the density of the matrix.
Based on the models above the following measurements were derived.
Since both fields of investigation appeared to have the same geometrical shape, the combined measurements may be used to determine formation properties. To determine porosity, the following observation method was used.
A plot of the natural log of the detector count rates is shown in FIG. 5. Borehole and liquid changes cause the NN and NG lines to shift. Porosity changes cause the lines to change slope.
The amount of slope change is related to the porosity. Slope change can be computed by Y−X or

Y−X=ln LNG−ln LNN−(ln SNG−ln SNN)

This lead to the derivation of quad porosity (QTP).
Quad Neutron Porosity
One of important formation characteristic to be measured is porosity and the quad neutron logging tool allows for porosity to be characterized based on the field information collected by the quad neutron logging tool. The effect of porosity on neutron gamma and neutron thermal neutron was studied. The results showed that the porosity error for neutron gamma logging data were in the opposite direction from the porosity error for neutron thermal neutron logging data. A combining formula was created to minimize porosity error. The formula for Quad porosity (QTP) is as follows:




  
    Q
    ⁢
    
      
    
    ⁢
    T
    ⁢
    
      
    
    ⁢
    P
  
  =
  
    
      A
      *
      
        ln
        ⁡
        
          (
          
            
              SNN
              *
              LNG
            
            
              SNG
              *
              LNN
            
          
          )
        
      
    
    +
    B
  




where A and B are determined empirically and are influenced by the bit size and the amount of iron or high density material in the volume of investigation. The graph of FIG. 7 shows the relationship on QTP gain (A) vs. various casing weights.

Changing the bit size and therefore the borehole create shifts in computed porosity. The following values outlined in Table 1 were determined empirically for A values at 65 and are commonly used as a starting point for porosity. Offset porosity data should be used to properly set the shift if the current casing/bit size configuration has not been previously encountered.







TABLE 1







Empirical values for A and B based on bit size and casing diameter.












Csg
Bit
A
B
















41⁄2′′
61⁄4′′
65
−6



51⁄2′′
77⁄8′′
65
−9



7′′
83⁄4′′
65
−12










Quad Neutron Clay
The presence of clay can cause errors in porosity measurement. Research in diamond drilled mining holes led to the derivation of the Quad neutron clay calculation. The derivation is based on the fact that most clays contain aluminum and that the aluminum can become radioactive when subjected to a neutron flux. The increased gamma activity can therefore be related to clay. The equation used to compute Quad Neutron Clay (QNC) in porosity units is as follows:




  QNC
  =
  
    
      A
      *
      
        (
        
          
            C
            
              (
              
                SNG
                LNG
              
              )
            
          
          +
          
            D
            *
            
              (
              
                SNN
                LNN
              
              )
            
          
        
        )
      
    
    +
    B
  




where A, B, C and D are empirically determined coefficients. Typical values for A, B, C and D are 0.004, 0, 1.9 & −1.5 respectively. Determination of absolute value can be determined from offset core data or local geological experience.

By absolute it is meant best fit relative to the core clay volume and core porosity. Clay volume may be measured in a laboratory to determine the physical amount of clay. The values from the core data can be plotted against the QNC curve. The parameters A, B, C and D may be adjusted until the QNC curve has a best fit relation with the core data. To do this process, B is left at 0, and then C and D are adjusted until the trend is matched and then A is adjusted such that the Core Porosity equals QTP−QNC.
Comparisons to core data confirm the accuracy of the measurement.
Quad Neutron Liquid
One method of determining fluid types and determining saturation, relies on the determination of porosity error. In June of 1983, Roke presented a paper at the SPWLA 24th Annual Logging Symposium titled “The use of Multiple Through Casing Porosity Logs to Quantify Water Saturation”. The paper is based on the concept of determining saturation based on porosity error due to liquid changes. The Quad Neutron Liquid, QNL, calculation combines the neutron thermal neutron and neutron gamma physics to maximize the porosity error. The formula is based on the observation that liquid changes cause the NN and NG lines in FIG. 15 to shift. This led to the X+Y calculation which led to the QNL calculation as follows:




  
    Q
    ⁢
    
      
    
    ⁢
    N
    ⁢
    
      
    
    ⁢
    L
  
  =
  
    
      A
      *
      
        ln
        ⁡
        
          (
          
            
              LNG
              *
              SNG
            
            
              LNN
              *
              SNN
            
          
          )
        
      
    
    +
    B
  




where A and B are chosen to create the best overlay of QNL to QTP. The effect of varying water salinity to 32 API oil are shown in the plot of FIG. 8 and are based on laboratory test data. The equation displayed in the plot of FIG. 8 can be used to determine the expected porosity value x based on a given kppm salinity y.

The Quad Neutron Liquid curve is sensitive to density changes and should be corrected before use in saturation calculations.
Dual Detector Neutron
The Dual Detector Neutron, DDN, is used to determine relative hydrogen index when overlaid with the QTP curve outlined above. The calculation is similar to the basic neutron porosity calculations used in the industry. It is presented here for completeness of the entire process. The Dual Detector Neutron is determined as follows:




  DDN
  =
  
    
      A
      *
      
        ln
        ⁡
        
          (
          
            SNN
            LNN
          
          )
        
      
    
    +
    B
  




where A and B are chosen to create the best overlay of DDN to QTP.

Saturation
One of the objectives of formation and reservoir evaluation is to determine the Oil Saturation. Typically, the water saturation, Sw, is determined and oil is assumed as 1−Sw. Laboratory testing has shown that the saturation relationship for the quad neutron data can be approximated with a linear relationship. The Sw calculation for the Quad Neutron is determined as follows:




  Sw
  =
  
    1
    -
    
      kf
      ⁢
      
        
          (
          
            
              Q
              ⁢
              
                
              
              ⁢
              E
              ⁢
              
                
              
              ⁢
              P
            
            -
            
              Q
              ⁢
              
                
              
              ⁢
              E
              ⁢
              
                
              
              ⁢
              L
            
          
          )
        
        
          Q
          ⁢
          
            
          
          ⁢
          E
          ⁢
          
            
          
          ⁢
          P
        
      
    
  




where, k represents fluid factor, f represents formation factor, QEP represents Quad Effective Porosity (Clay Free) and QEL represents Quad Effective Liquid porosity (Clay Free). QEP and QEL are determined as follows:

QEP=QTP−QC and QEL=QL−QC

where, QTP represents Quad Total Porosity, QC represents Quad Clay Porosity as referred to herein as QNC and QL represents Quad Liquid Porosity also referred to herein as QNL.

The fluid factor is determined as k=1/(MaxPor*QLgain)

where, MaxPor=WaterPor−OilPor and WaterPor is given as

WaterPor=(12.73+√{square root over (12.732−4*0.0966*(403.04−Waterkppm)))}/(2*0.0966)

OilPor is given as

OilPor=0.1333*oilAPI+71

QLgain is determined as

QLgain=QNL A value/32.5

The formation factor is determined through local experience in the same manner as the formation factor in Archie's equation, the industry standard Sw calculation. The Quad Liquid curve used in the above formulas should be corrected for density changes. The Quad liquid curve will measure too low porosity in high dense materials and too high porosity in low dense materials. The correction is typically done utilizing the chemical effect curve or a hybrid between the relative bulk density curve and the chemical curve.
Ancillary Measurements
Given that Quad Neutron Porosity is computed and accepted as porosity, the individual sensor readings are converted into porosity by creating a best fit relationship between the detector counts and porosity. Once this has been established, calculations can be made in porosity units between the single detectors. Of these, Relative Bulk Density, Nuclear Caliper and Chemical Effect are determined and elaborated upon.
Relative Bulk Density
One method of determining bulk density according to the prior art utilizes gamma gamma measurements that are influenced by the casing. The Quad relative bulk density measurement is based on relationships between neutron neutron and neutron gamma fields. Using observation and empirical methods, a relationship was discovered between the LNN field readings and the LNG field reading. Once both detectors have been converted into porosity, the relative bulk density is determined as follows:

Relative Bulk Density=LNNpor−LNGpor

where, LNNpor and LNGpor represent the porosity converted counts.

An absolute bulk density can be determined if offset or bulk density data is available to calibrate with.
Nuclear Caliper
Another formation evaluation made capable by the Quad Neutron system and method is the ability to measure the borehole size behind pipe. The method was derived again by utilizing the Quad Neutron nuclear model, observations and empirical methods. The observation was that the LNN detector was very sensitive to near wellbore effects and therefore changes in the borehole size will influence the count rates detected. The Nuclear Caliper, NC, is determined as follows:

NC=SNNpor−LNNpor

where SNNpor and LNNpor represent the porosity converted counts. The raw SNN and LNN counts are converted to calibrated filtered values as will be appreciated by one skilled in the art. The calibrated filtered values are then converted to porosity by adjusting gain and shift values for each of the single detector counts until the curve has a best fit relationship with QTP over a chosen interval. The interval is chosen such that the computation of NC and CE create a may west effect in low porosity calcite intervals.

Chemical Effect
The chemical effect curve was derived to measure the chemical changes in the formation. Through observation and empirical methods, the curve was seen to be sensitive to hydrocarbon and was subsequently used to provide a correction to the Quad Neutron Liquid curve.
The Chemical Effect Curve, CE, is determined as follows:

CE=SNNpor−LNGpor

where SNNpor and LNGpor represent the porosity converted single detector counts.

Flowability (Permeability)
Another measurement made possible by the quad neutron system and method is the dual gamma ray method of determining permeability. The method is based on work done on water well flowability, where it was observed that water well flow rates could be related to gamma ray levels. The basic theory is that fine highly transportable radioactive material, fines, becomes stuck in the muds of low energy depositional environments. The mud reduces permeability in the formation and therefore being able to measure the fines will allow the ability to infer permeability. The gamma energy level of these fines were thought to exist between 80 kev and 180 kev. To isolate these energies, two gamma detectors are used. A mechanical filter was designed to minimize this energy range. The filter was applied to one detector and the other detector was left unshielded. The measurements are compared to each other and as the two measurements approach, the amount of fines detected decreased which in turn infers permeability. Comparisons to open hole SP measurements has shown accurate results. Where the measurement is limited is in calcified sands as the calcification does not create a change in fines. In these instances however the quad porosity can be seen to decrease and the relative bulk density increase.
Attempts to replicate the shielding with an energy cut-off filter did not work. The reason is that the transmissivity of the mechanical filter is not a simple linear energy filter and therefore cannot be replicated with as simple energy cut-off. It has been theorized that an electrical or software filter could be made to replicate the shielding used on the filtered gamma ray but at time of writing had not been implemented. The notebooks on the design of the filters is included as part of the patent application.
Some advantages associated with the methods, systems and logging tools disclosed herein is the accurate formation evaluation carried out through-pipe. The methods, systems and logging tools of the present invention have the ability to recognize oil and gas where other technologies cannot. The other advantage is that it is safer to acquire data in horizontal wells because measurements through-pipe can be carried out which reduces the risk of losing a radioactive source in the well as a result of caving and hole sloughing. Should the hole slough or cave and the drill pipe becomes stuck, the Quad Neutron tool may be retrieved out of the pipe. Another advantage is in re-evaluating existing cased wells and identifying hydrocarbon and missed pay opportunities.
Examples
The following examples are a collection of measurements done using Quad Neutron technology and quad neutron logging tools and systems.
Case History 1. Quad Neutron identifies by passed oil in low resistivity formation. Increased oil production from 8 cubic meters per day to 20 cubic meters per day with no increase in water. Core comparisons. The example below shows the Quad Neutron Clay free porosity curve plotted against core porosity data along with open hole neutron density porosity on both sandstone and limestone matrix and is illustrated in FIG. 9. The quad neutron data was collected through the casing.
From the same well, the core clay volume is shown plotted against the Quad Neutron Clay curve which is illustrated in FIG. 10. Again, there is excellent correlation to clay volumes and the Quad Neutron measurement was done through the casing. Investigating the difference at 2388 m, it was learned that the core sample was taken in a thin tight carbonate layer that was approximately 5 cm thick which was much smaller than the 1 m vertical bed resolution of the Quad Clay measurement and hence the measurement was not fully resolved during logging.
In another core example, the Quad Relative Bulk Density curve (labelled as CE) is shown compared to the open hole bulk density curve (RHOB) and the bulk density obtained from core analysis in FIG. 11. There is good correlation between core and the Quad measurement and again the Quad measurement was made through the pipe.
In the example below, the Quad Neutron was logged to determine the current in situ saturation. Initial production was reported as 277 t/day of water and 7.4 t/day of oil for a 97.4% Sw. The Quad log showed that the upper interval had higher oil saturations than the lower. After performing a remediation on the well, only the upper 7 m remained open. The well was put back on production and was producing 25 t/day of water and 11 t/day of oil for a 69.4% Sw. Water was reduced by 90% and oil was increased by 48%. Again, the Quad measurement was done through the casing. The log is illustrated in FIG. 12.
The last example data set is shown below as FIG. 13. This horizontal well was initially drilled and left standing because of poor to no oil shows during drilling. Tubing was run in the well with a perforated pup joint on the bottom. A memory quad neutron was pumped down on slickline to the end of the tubing and then logged out with slickline. The resulting log is shown in FIG. 13.
Analysis of the log showed oil in the heel of the well and then large amounts of water after 1950 m. To test the results, a pump was installed and the well was put on production. Initial oil production climbed to over 100 bbls of oil/day at a water cut of 68% and then settled to 50 bbls/day oil at a water cut of 80% at the end of 45 days. Production was halted and an inflatable packer was set at 1950 m. The well was put back on production and after two weeks was producing 50 bbls/day of oil at a water cut of 50%. A decision was made to increase the pumping rate and production went to over 100 bbls/day oil while maintaining the 50% water cut. The high pumping rate continued until the packer failed dropping the oil production to less than 10 bbls/day with a water cut in excess of 95%.
As can be seen from the production results, outlined in FIG. 14, the Quad Neutron system correctly identified the region of highest oil saturation and allowed the producer to increase oil production while reducing their water production. And again the measurement was done through the pipe.
It will be appreciated that the theories outlined herein are non-limiting and it is not the inventors intent to be bound by theory. The analysis including the comparative analysis demonstrate the accuracy and utility of the method, system and/or logging device without the need to be limited by the theory.

## Claims

1. A method of evaluating a characteristic of a formation adjacent a wellbore comprising:
emitting a neutron field into the formation from a neutron source;
obtaining data comprising a long neutron neutron (LNN), short neutron neutron (SNN), long neutron gamma (LNG) and short neutron gamma data (SNG) from the neutron field emitted into the formation:
combining the long data with the short data to reduce error associated with the characteristic to be evaluated, wherein the characteristic is neutron porosity (QTP) and the formula for evaluating the characteristic QTP is:





  QTP
  =
  
    
      A
      *
      
        ln
        ⁡
        
          (
          
            
              SNN
              *
              LNG
            
            
              SNG
              *
              LNN
            
          
          )
        
      
    
    +
    B
  




wherein A and B are determined empirically and are influenced by bit size and the amount of high dense material in a volume of investigation.

2. The method of claim 1, further comprising the step of:
optimizing a spacing of a long neutron detector for detecting LNN and a short neutron detector for detecting SNN relative the neutron source to minimize the error associated to the characteristic to be evaluated.

3. The method of claim 1, further comprising the step of:
optimizing a spacing of a long gamma detector for detecting LNG and a short gamma detector for detecting SNG relative the neutron source to minimize the error associated to the characteristic to be evaluated.

4. The method of claim 1, wherein the data is obtained through a wellbore casing.

5. The method of claim 1, wherein the combining step comprises the cancelation of the error associated to the long data with the error associated to the short data.

6. A method of evaluating a characteristic of a formation adjacent a wellbore comprising:
emitting a neutron field into the formation from a neutron source;
obtaining data comprising a long neutron neutron (LNN), short neutron neutron (SNN), lone neutron gamma (LNG) and short neutron gamma data (SNG) from the neutron field emitted into the formation;
combining the long data with the short data to reduce error associated with the characteristic to be evaluated, wherein the characteristic is neutron clay (QNC) and the formula for evaluating the characteristic QNC is:





  QNC
  =
  
    
      A
      *
      
        (
        
          
            C
            
              (
              
                SNG
                LNG
              
              )
            
          
          +
          
            D
            *
            
              (
              
                SNN
                LNN
              
              )
            
          
        
        )
      
    
    +
    B
  




wherein A, B, C and D are empirically determined coefficients.

7. The method of claim 6, wherein A, B, C and D are 0.004, 0, 1.9 and −1.5, respectively.

8. The method of claim 1, wherein the characteristic is neutron liquid (QNL) and the formula for evaluating the characteristic QNL is:




  
    Q
    ⁢
    
      
    
    ⁢
    N
    ⁢
    
      
    
    ⁢
    L
  
  =
  
    
      A
      *
      
        ln
        ⁡
        
          (
          
            
              LNG
              *
              SNG
            
            
              LNN
              *
              SNN
            
          
          )
        
      
    
    +
    B
  




wherein A and B are chosen to create the best overlay of QNL to QTP.

9. The method of claim 1, wherein the characteristic is water saturation (Sw) and the formula for evaluating the characteristic Sw is:




  Sw
  =
  
    1
    -
    
      Kf
      ⁢
      
        
          (
          
            
              Q
              ⁢
              
                
              
              ⁢
              E
              ⁢
              
                
              
              ⁢
              P
            
            -
            
              Q
              ⁢
              
                
              
              ⁢
              E
              ⁢
              
                
              
              ⁢
              L
            
          
          )
        
        
          Q
          ⁢
          
            
          
          ⁢
          E
          ⁢
          
            
          
          ⁢
          P
        
      
    
  




wherein
k=fluid factor;
f=formation factor;
QEP=Quad Effective Porosity (Clay Free); and
QEL=Quad Effective Liquid porosity (Clay Free); and
wherein
QEP and QEL are determined as follows:


QEP=QTP−QC and QEL=QL−QC

wherein
QTP=Quad Total Porosity;
QC=Quad Clay Porosity; and
QL=Quad Liquid Porosity;


the fluid factor k=1/(MaxPor*QLgain)


wherein
MaxPor=WaterPor−OilPor;
WaterPor is determined as follows:


WaterPor=(12.73+√{square root over (12.732−4*0.0966*(403.04−Waterkppm)))}/(2*0.0966)


OilPor is determined as follows:


OilPor=0.1333*oilAPI+71 and


QLgain is computed as


QLgain=QNL A value/32.5.

10. The method of claim 6, further comprising the step of:
optimizing a spacing of a long neutron detector for detecting LNN and a short neutron detector for detecting SNN relative the neutron source to minimize the error associated to the characteristic to be evaluated.

11. The method of claim 6, further comprising the step of:
optimizing a spacing of a long gamma detector for detecting LNG and a short gamma detector for detecting SNG relative the neutron source to minimize the error associated to the characteristic to be evaluated.

12. The method of claim 6, wherein the data is obtained through a wellbore casing.

13. The method of claim 6, wherein the combining step comprises the cancelation of the error associated to the long data with the error associated to the short data.