Abstract:
The invention considers the frequency distributions of singles, doubles and triple neutron emission events from a sample under assay. The count rates are equated to mathematical functions related to the spontaneous fission rate, self-induced fission rate, detection efficiency and α,n rate with probability distribution assigned to each of those factors, the value of the product of all the probability distributions being increased to give an optimized solution and so provide a value of the spontaneous fission rate which is linked to the mass of the neutron source. The technique aims to provide increased accuracy and certainty compared with neutron coincidence counting based techniques.

Description:
CROSS-REFERENCE TO RELATED APPLICATIONS  
       [0001]     This application is a continuation of U.S. patent application Ser. No. 10/930,445, filed Aug. 30, 2004, which is a divisional of U.S. patent application Ser. No. 10/075,965, filed Feb. 13, 2002, now abandoned, which is a continuation of U.S. patent application Ser. No. 09/554,142, filed Oct. 12, 2000, now abandoned, which claims priority to International Application No. PCT/GB98/03656, filed Dec. 11, 1998, which claims priority to United Kingdom Application Nos. 9810433.4, filed May 15, 1998 and 9726270.3, filed Dec. 12, 1997, which for purposes of disclosure are incorporated herein by specific reference. 
     
    
     BACKGROUND OF THE INVENTION  
       [0002]     1. The Field of the Invention  
         [0003]     This invention is concerned with improvements in and relating to monitoring, in particular, but not exclusively, monitoring of nuclear materials.  
         [0004]     2. The Relevant Technology  
         [0005]     Monitoring of nuclear waste, and plutonium waste in particular, is important both for accountancy purposes and in ensuring that criticality safety standards are maintained. It is important to maintain an inventory of the plutonium in a plant and during disposal, for instance, to know the exact amount of plutonium in a given volume of waste.  
       SUMMARY OF THE INVENTION  
       [0006]     Determination of the mass of plutonium in a waste sample presently conducted using neutron coincidence counting (NCC). On occasion it is necessary to supplement these investigations with high resolution gamma spectrometry measurement to check the validity of assumptions on the isotopic make up of the sample. For waste samples such as drums, boxes and other handleable packages a detection chamber is normally constructed into which the samples to be measured are introduced.  
         [0007]     The detection chamber principally consists of a series of neutron detecting tubes placed in a polyethylene moderating material. The signals produced by the detectors are analyzed electronically to determine those signals attributable to spontaneous fissions of the plutonium rather than background radiation and other sources. The spontaneous fissions of the plutonium produce 2, 3, 4 and higher numbers of neutrons for each event. The spontaneous fission rate is known to be proportional to the mass of certain isotopes and, therefore, to the neutron source material as a whole. Signals detected within a certain time window are taken as being indicative of a pair, rather than a single neutron, which can neither be attributed to the background or the plutonium with any certainty.  
         [0008]     As the component of the signal relating to the pairs is indicative of the plutonium mass of certain isotopes, principally  240 Pu then by making certain assumptions about the isotopic makeup of the sample the total mass can be suggested.  
         [0009]     As well as requiring such assumptions as to the material makeup which are inappropriate given the likely variation over time of the isotopes in the waste, NCC also suffers significant errors or inaccuracies due to the variable shielding effects of the different waste materials in which the plutonium of interest may be distributed as well as due to the different materials associated with the containers for the waste. For instance water, which will be present to varying degrees in the waste, and polyethylene or PVC, which may be associated with the container, are each strong neutron moderating materials. Variation in any of these can drastically affect the detection efficiency and hence the result. Problems with the effects of changes in the detection efficiency due to the materials actual position within the sample and as a consequence within the chamber are also encountered.  
         [0010]     As an improvement it has been proposed to undertake multiplicity counting (MC) in which the events producing single and triplet events are considered alongside the doublets. However, the count rates arising relate to four unknown parameters; the spontaneous fission rate (reflecting the isotopic make up of the material); the self-induced fission rate (known as multiplication); the detection efficiency; and the α,n reaction rate. Several other nuclear constants are also involved. The presence of three experimentally determinable factors, the single, double and triple count rates, still leaves the system unsolvable as there are four unknowns. Detection of the quadlet counting rate is impractical due to the very limited number of detectable occurrences for such events. In even the best prior art systems, therefore, it is necessary to set one the unknowns at a predetermined value in reaching a solution.  
         [0011]     Attempts to set any one of the self-induced fission rate; detection efficiency; or α,n reaction rate are not accurate in waste samples.  
         [0012]     The self-induced fission rate varies due to variations in materials and locations with the result that even a small error has a substantial effect on the end result due to the factors role in the equations. The detection efficiency differs from that obtained in calibration due to variations in the location of the plutonium in the waste sample and due to variations within the accompanying materials. The α,n rate changes with the isotopic, chemicals and materials in the waste matrix which will vary between samples. In each case errors in the end result arise resulting in inefficient and costly waste disposal schemes or an inaccurate inventory.  
         [0013]     According to a first aspect of the invention we provide a method of monitoring a sample containing a neutron source in which: 
        i) signals from a plurality of neutron detectors are analyzed and the count rates for single, double and triple incidence of neutrons on the detectors are determined;     ii) the single, double and triple count rates are equated to a mathematical function related to the spontaneous fission rate, and self-induced fission rate, detection efficiency and α,n rate;     iii) a probability distribution is assigned to each of the self-induced fission rate, detection efficiency, α,n reaction rate, and each of the counting rates to provide a probability distribution factor for any given value;     iv) and the value of the product of all the probability distribution factors is increased to give an optimized solution and so provide a value for the spontaneous fission rate which is linked to the mass of the neutron source.        
 
         [0018]     In this way an accurate value for the mass of a neutron source within a sample of material is provided without having to make hard assumptions about any of the variables within the system.  
         [0019]     Quadruplet or high counts could be substituted for one of the single, pair of triple counts, but the number of events unit time would be reduced accordingly.  
         [0020]     The sample may be in the form of a piece of equipment, such as a glove box, ventilation unit or the like; a component of a piece of equipment, such as a filter; or a container, such as a drum, box or other package. The container may be sealed. The sample may be of waste material, such as material destined for disposal, or it may be material with a future active life for which an inventory is required.  
         [0021]     The sample may include other materials besides the neutron source. For instance other materials such as non-neutron source elements, metals or compounds, water, plastics, glass and other sealing materials may be present.  
         [0022]     The neutron source may comprise one or more elements or compounds, one or more isotopes of an element or mixtures of both. The neutron source may be naturally occurring and/or arise from fission reaction products. Plutonium and  240 Pu in particular are neutron sources which may require such monitoring.  
         [0023]     The neutron detectors may be of the  3 He type. The detectors may be provided in polyethylene or other hydrogen providing material. In this way the neutron source to be monitored can be controlled and the neutrons moderated to detectable energy levels. Alternative or additional neutron absorbing materials, such as boron or cadmium, may be provided to shield the detectors, for instance by positioning around the detectors, against background events.  
         [0024]     The detectors may be provided around all sides and most preferably above and/or below the sample during monitoring. Between 20 and 125 detectors may be provided. Preferably between 30 and 50 detectors are provided around the sides of the sample. Preferably between 8 and 16 detectors are provided above and below the sample.  
         [0025]     Preferably each detector is provided with an amplifier (generally to increase the pulse to a signal level suitable for further processing) and a discriminator (generally to reject or eliminate noise events).  
         [0026]     The detector signals are preferably summed and fed to sequence analyzing means. Preferably each pulse causes a time period to be considered, with other pulses being received in that period being associated with the initial pulse. In this way sequences of single, double, triple and greater numbers of neutron detections are obtained. Neutrons detected within the time period are accepted as originating from the same spontaneous fission within the sample as the initial pulse. Preferably the detections are subjected to a correction factor. The correction factor may account for accidental coincidences, for instance detections from different source simultaneous emissions which give the impression of a pair or greater.  
         [0027]     Preferably the time period lasts between 10 and 500 micro seconds and most preferably between 50 and 100 micro seconds.  
         [0028]     Preferably the start of the time period occurs between 5 and 10 micro seconds after the initial pulse.  
         [0029]     Preferably the singlet count rate is related to the spontaneous fission rate, the self-multiplication factor, where  
       m   =       1   -   p         (     1   -   p     )     ⁢     u   I             
 
 and p=probability first neutron causes induced fission; the detection efficiency and the α, n reaction rate by the function 
 
 R   1   =ε·F   ε   ·M·V   ε1 ·(1+α) 
 
         [0030]     Preferably the doublet counting rate is related to the spontaneous fission rate, the self-multiplication factor, where  
       m   =       1   -   p         (     1   -   p     )     ⁢     u   I             
 
 and p=probability first neutron causes induced fission; the detection efficiency and the α,n reaction rate by the function  
         R   2     =       ɛ   2     ·     F   s     ·     M   2     ·     v     s   ⁢           ⁢   2       ·     (     1   +       (     M   -   1     )     ⁢     (     1   +   α     )     ⁢         v     s   ⁢           ⁢   1       ⁢     v     I   ⁢           ⁢   2             v     s   ⁢           ⁢   2       ⁡     (       v     I   ⁢           ⁢   1       -   1     )             )           
 
         [0031]     Preferably the triplet counting rate is related to the spontaneous fission rate, the self-multiplication factor, where  
       m   =       1   -   p         (     1   -   p     )     ⁢     u   I             
 
 and p=probability first neutron causes induced fission; the detection efficiency and the α, n reaction rate by the function  
         R   3     =       ɛ   3     ·     M   3     ·     v     s   ⁢           ⁢   3       ·     (     1   +     2   ⁢     (     M   -   1     )     ⁢         v     s   ⁢           ⁢   2       ⁢     v     I   ⁢           ⁢   2             V     s   ⁢           ⁢   3       ⁡     (       v     I   ⁢           ⁢   1       -   1     )           +       (     M   -   1     )     ⁢     (     1   +   α     )     ⁢         v     s   ⁢           ⁢   1       ⁢     v     I   ⁢           ⁢   3             v     s   ⁢           ⁢   3       ⁡     (       v     I   ⁢           ⁢   1       -   1     )         ⁢     (     1   +     2   ⁢     (     M   -   1     )     ⁢       v     I   ⁢           ⁢   2     2         v     I   ⁢           ⁢   3       ⁡     (       v     s   ⁢           ⁢   1       -   1     )             )         )           
 
         [0032]     Preferably the probability distribution assigned individual variables or counting rates is a normal distribution or a flat distribution or a triangular distribution. Normal distributions are preferably used for one or more, and most preferably all, the counting rates. Triangular distributions are preferably used for one or more, and most preferably all, the individual variables, such as detector efficiency, fission rate, multiplication distribution and alpha distribution. A flat distribution is preferably used for the fission rate.  
         [0033]     The probability distributions may be symmetrical. The probability distributions may be skewed.  
         [0034]     The distribution may be constrained within certain applied constraints/boundaries such that the probability distribution factor is zero beyond the constraints, particularly for triangular distributions.  
         [0035]     The distribution may be modified such that the probability distribution factor rapidly tends to zero beyond certain values, particularly for normal distributions.  
         [0036]     The constraints beyond which the distribution is zero or values beyond which a rapid transition to zero occurs, may be selected from one or more of the following: 
        for the detection efficiency 0≦E≦1;     for the self-induced fission rate M≧1     for the α, n reaction α≧0.        
 
         [0040]     The ranges of the distributions of one or more of the single count rate; the double count rate; and the triple count rate are preferably based on the measured count rates and their standard deviation.  
         [0041]     The lower and/or upper values of the distributions, may be provided according to one or more or all of the following ranges: 
        a) for the self-induced fission rate: according to the anticipated plutonium content; between 1.0 and 1.3, preferably between 1.0 and 1.2; particularly for a high plutonium content (i.e., &gt;100 g Pu) between 1.0 and 1.2; particularly for a low plutonium content (i.e., &lt;10 g Pu) between 1.00 and 1.01; and for intermediate plutonium contents, intermediate ranges;     b) for the detection efficiency: according to the anticipated moderator content; between 0 and 0.3; between 0 and 0.15 for anticipated maximum efficiency levels; between 0.05 and 0.15 for moderator containing materials;     c) for the α,n reaction rate: according to the form or forms of plutonium anticipated; between 0 to 10, where plutonium fluorides are anticipated as present; between 0 to 2, or 0 to 1 where plutonium metal is anticipated; greater than 0.0 where no plutonium metal is anticipated, for instance between 0.4 and 1.0 or 0.5 and 1.0.        
 
         [0045]     One or more of the constraints may be set according to information gathered from a preceding isotopic consideration or analysis of the sample. Gamma spectrometry may be employed. Information about the presence of plutonium metal and/or plutonium oxides and/or plutonium fluorides and/or moderators and/or other species or components may be obtained and applied.  
         [0046]     The permissible values for the spontaneous fission rate may be constrained. Preferably the spontaneous fission rate is constrained to F≧0.  
         [0047]     The distributions may be provided according to the general form:  
         pdf   ⁡     (   ɛ   )       =       1       2   ⁢   π         ⁢     exp   ⁡     [       -     1   2       ·         (     ɛ   -     μ   ɛ       )     2       σ   ɛ         ]             
 
         [0048]     The initial set of values assigned to the variables may be selected by the operator, predetermined or a function of prior samples which have been analyzed.  
         [0049]     The increasing, and preferably maximization, of the product of the probability distribution factors (pdf&#39;s) is preferably performed as an iterative process. The values of one, two, three and most preferably all four of the variables are varied in each iteration. The variation or correction applied may be of a fixed value between iterations but a variable correction is preferred.  
         [0050]     One or more of the probability distribution factors and/or the corresponding level of the variable from an optimized or maximized solution may be used to redefine the range of the applicable distribution used. One or more redefined distributions may be used in a subsequent optimization or maximization according to the method. Tighter errors apply as a result as tighter constraints would be applied. The redefining of the distribution(s) and their use in optimization may be performed repeatedly.  
         [0051]     The correction may be a correction vector. Most preferably the correction vector is derived from d(pdf)/d(x) where x is the variable, solved for zero.  
         [0052]     In one method, because there are four variables, the part function may be partially differentiated to yield a set of four equations:  
         [0053]     e.g. 
        μ ε =0.12 σ ε =0.03     μ m =1.02 σ m =0.02     μ α =0.4 σ α =0.1        
 
         [0057]     The product pdf may be partially differentiated to provide four simultaneous equations. The equations may be  
               ∂     (       ∂   pdf       ∂   ɛ       )         ∂   ɛ       ⁢   Δɛ     +         ∂     (       ∂   pdf       ∂   ɛ       )         ∂     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         ∂     (       ∂   pdf       ∂   ɛ       )         ∂   M       ⁢   Δ   ⁢           ⁢   M     +         ∂     (       ∂   pdf       ∂   ɛ       )         ∂   α       ⁢   Δα       =     -       ∂   pdf       ∂   ɛ             
               ∂     (       ∂   pdf       ∂     F   s         )         ∂   ɛ       ⁢   Δɛ     +         ∂     (       ∂   pdf       ∂     F   s         )         ∂     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         ∂     (       ∂   pdf       ∂     F   s         )         ∂   M       ⁢   Δ   ⁢           ⁢   M     +         ∂     (       ∂   pdf       ∂     F   s         )         ∂   α       ⁢   Δα       =     -       ∂   pdf       ∂     F   s               
               ∂     (       ∂   pdf       ∂   M       )         ∂   ɛ       ⁢   Δɛ     +         ∂     (       ∂   pdf       ∂   M       )         ∂     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         ∂     (       ∂   pdf       ∂   M       )         ∂   M       ⁢   Δ   ⁢           ⁢   M     +         ∂     (       ∂   pdf       ∂   M       )         ∂   α       ⁢   Δα       =     -       ∂   pdf       ∂   M             
               ∂     (       ∂   pdf       ∂   α       )         ∂   ɛ       ⁢   Δɛ     +         ∂     (       ∂   pdf       ∂   α       )         ∂     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         ∂     (       ∂   pdf       ∂   α       )         ∂   M       ⁢   Δ   ⁢           ⁢   M     +         ∂     (       ∂   pdf       ∂   α       )         ∂   α       ⁢   Δα       =     -       ∂   pdf       ∂   α             
 
         [0058]     Alternatively or additionally the effects of small variations can be determined based on linearized equations in four dimensions. The use of equations  
               δ   ⁡     (       δ   ⁢           ⁢   pdf     δɛ     )       δɛ     ⁢   Δɛ     +         δ   ⁡     (       δ   ⁢           ⁢   pdf     δɛ     )         δ   ⁢           ⁢     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         δ   ⁡     (       δ   ⁢           ⁢   pdf     δɛ     )         δ   ⁢           ⁢   M       ⁢   Δ   ⁢           ⁢   M     +         δ   ⁡     (       δ   ⁢           ⁢   pdf     δɛ     )       δα     ⁢   Δα       =     -       δ   ⁢           ⁢   pdf     δɛ           
               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s         )       δɛ     ⁢   Δɛ     +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s         )         δ   ⁢           ⁢     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s         )         δ   ⁢           ⁢   M       ⁢   Δ   ⁢           ⁢   M     +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s         )       δα     ⁢   Δα       =     -       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s               
               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M       )       δɛ     ⁢   Δɛ     +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M       )         δ   ⁢           ⁢     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M       )         δ   ⁢           ⁢   M       ⁢   Δ   ⁢           ⁢   M     +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M       )       δα     ⁢   Δα       =     -       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M             
               δ   ⁡     (       δ   ⁢           ⁢   pdf     δα     )       δɛ     ⁢   Δɛ     +         δ   ⁡     (       δ   ⁢           ⁢   pdf     δα     )         δ   ⁢           ⁢     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         δ   ⁡     (       δ   ⁢           ⁢   pdf     δα     )         δ   ⁢           ⁢   M       ⁢   Δ   ⁢           ⁢   M     +         δ   ⁡     (       δ   ⁢           ⁢   pdf     δα     )       δα     ⁢   Δα       =     -       δ   ⁢           ⁢   pdf     δα           
 
 where  
             ∂     (       ∂   pdf       ∂   ɛ       )         ∂     F   s         ≈       δ   ⁡     (       δ   ⁢           ⁢   pdf     δɛ     )         δ   ⁢           ⁢     F   s           =       (         δ   ⁢           ⁢     pdf   ⁡     (     ɛ   ,     F   s     ,   M   ,   α     )         δɛ     -       δ   ⁢           ⁢     pdf   ⁡     (     ɛ   ,       F   s     +     Δ   ⁢           ⁢     F   s         ,   M   ,   α     )         δɛ       )       δ   ⁢           ⁢     F   s             
 
 is the 2 nd  order derivative of  
           (       ∂   pdf       ∂   ɛ       )     ≈     (       δ   ⁢           ⁢   pdf     δɛ     )       =               pdf   ⁡     (       ɛ   +   δɛ     ,     F   s     ,   M   ,   α   ,     μ     R   1       ,     μ     R   2       ,     μ     R   3         )       -               pdf   ⁡     (       ɛ   +   δɛ     ,     F   s     ,   M   ,   α   ,     μ     R   1       ,     μ     R   2       ,     μ     R   3         )             δɛ         
 
 etc. may be employed. These functions may be solved to give a correction vector, for instance  
         [           Δ   ɛ               Δ     F   s                 Δ   m               Δ   α           ]             
 
         [0059]     The method may handle the solution in matrix form.  
         [0060]     Preferably the process is repeated until the correction vector is below a given threshold. The threshold may be predetermined.  
         [0061]     The initial pdf values for the allocated initial values are preferably evaluated as positive and negative correction vectors. Preferably each subsequent iterations pdf values are similarly evaluated. In this way the effects of poor starting criteria are alleviated.  
         [0062]     For a pdf the correction vector may be multiplied by a constant factor and divided by a constant factor and the corresponding new pdf compared with the preceding pdf, if the pdf value is greater than the preceding pdf the new correction vector is applied, if the pdf value is less than the preceding value the new correction factor is divided by the constant once more and the new pdf far this further correction factor compared with the new preceding pdf, with the stages further being repeated if necessary. Preferably the multiplication factor is 32 and/or the divisive factor is 0.5. Most preferably these steps are performed for both positive and negative correction vectors.  
         [0063]     If one or more of the initially allocated values corresponds to a pdf of zero a new value may be allocated. Alternatively the count rate standard deviations may be multiplied by a constant factor, greater than 1. In this way a pdf value can be provided for which derivatives can be calculated. Preferably once the correction vector is below a certain threshold the standard deviation can be divided by the initial constant or a lower figure and the solution maximized once more. This process is preferably repeated until the inflationary factor is 1.  
         [0064]     Preferably the method includes the provision of associated error estimates. Total uncertainty within the system may be determined.  
         [0065]     The associated error estimates may be determined by equations  
           (     σ     ɛ   ^       )     2     =         (       (       ∂     ɛ   ^         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂     ɛ   ^         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂     ɛ   ^         ∂     R   3         )     ⁢     σ     R   3         )     2           
           (     σ       F   ^     s       )     2     =         (       (       ∂       F   ^     s         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂       F   ^     s         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂       F   ^     s         ∂     R   3         )     ⁢     σ     R   3         )     2           
           (     σ     M   ^       )     2     =         (       (       ∂     M   ^         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂     M   ^         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂     M   ^         ∂     R   3         )     ⁢     σ     R   3         )     2           
           (     σ     α   ^       )     2     =         (       (       ∂     α   ^         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂     α   ^         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂     α   ^         ∂     R   3         )     ⁢     σ     R   3         )     2           
 
         [0066]     Alternatively the associate error estimates may be obtained from the equations  
           (     σ     ɛ   ^       )     2     =         (       (       ∂     ɛ   ^         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂     ɛ   ^         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂     ɛ   ^         ∂     R   3         )     ⁢     σ     R   3         )     2           
           (     σ       F   ^     s       )     2     =         (       (       ∂       F   ^     s         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂       F   ^     s         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂       F   ^     s         ∂     R   3         )     ⁢     σ     R   3         )     2           
           (     σ     M   ^       )     2     =         (       (       ∂     M   ^         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂     M   ^         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂     M   ^         ∂     R   3         )     ⁢     σ     R   3         )     2           
           (     σ     α   ^       )     2     =         (       (       ∂     α   ^         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂     α   ^         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂     α   ^         ∂     R   3         )     ⁢     σ     R   3         )     2           
 
 where  
           (       ∂     ɛ   ^         ∂     R   1         )     ≈     (       δ   ⁢     ɛ   ^         δ   ⁢           ⁢     R   1         )       =           ɛ   ^     ⁡     (       μ     R   1       ,     δ   ⁢           ⁢     μ     R   2         ,     μ     R   2       ,     μ     R   3         )       -       ɛ   ^     ⁡     (       μ     R   ⁢           ⁢   1       ,     μ     R   ⁢           ⁢   2       ,     μ     R   ⁢           ⁢   3         )           δ   ⁢           ⁢     μ     R   ⁢           ⁢   1               
 
 etc. where δμ R     2    a small change in μ R     2    etc. ε(μ R     1   , μ R     2   , μ R     3   )=final estimate from solution to (μ R     1   , μ R     2   , μR   3   ) count rates set. 
 
         [0067]     According to a second aspect of the invention we provide: 
        i) apparatus for monitoring a sample containing a neutron source comprising one or more neutron detectors;     ii) processing means to receive signals from the detectors; and wherein the signals are processed to provide a count rate value indicative or the number of single, double and triple incidences of neutrons on the detectors; 
            wherein the counting rates are mathematically related to the spontaneous fission rate, self-induced fission rate, detection efficiency and α, n reaction rate;     wherein a probability distribution is assigned by the processing means to each of the self induced fission rate, detection efficiency and α,n reaction rate and the three counting rates; and     wherein the product of the probability distribution factors for values of the variables is increased by the processing means, the solution providing a value proportional to the mass of the neutron source present.    
               
 
         [0073]     The sample may be in the form of a piece of equipment, such as a glove box, ventilation unit or the like; a component of a piece of equipment, such as a filter; or a container such as a drum, box or other package. The container may be sealed. The sample may be of waste material, such as material destined for disposal, or it may be material with a future active life for which an inventory is required.  
         [0074]     The sample may include other materials besides the neutron source. For instance other materials such as non-neutron source elements, metals or compounds, water, plastics, glass and other sealing materials may be present.  
         [0075]     The neutron source may comprise one or more elements or compounds, one or more isotopes of an element or mixtures of both. The neutron source may be naturally occurring and/or arise from fission reaction products. Plutonium and  240 Pu in particular are neutron sources which may require such monitoring.  
         [0076]     The neutron detectors may be of the  3 He type. The detectors may be provided in polyethylene or other hydrogen providing material. In this way the neutron source to be monitored can be controlled and the neutrons moderated to detectable energy levels. Alternative or additional neutron shielding materials such as boron or cadmium may be provided to shield the detectors, for instance by positioning around the detectors.  
         [0077]     The detectors may be provided around all sides and most preferably above and/or below the sample during monitoring. Between 20 and 125 detectors may be provided. Preferably between 30 and 50 detectors are provided around the sides of the sample. Preferably between 8 and 16 detectors are provided above and below the sample.  
         [0078]     The detector signals are preferably summed and fed to sequence analyzing means. Preferably each pulse causes a time period to be considered, with other pulses being received in that period being associated with the initial pulse. In this way sequences of single, double, triple and greater numbers of neutron detections are obtained. Neutrons detected within the time period are accepted as originating from the same spontaneous fission within the sample as the initial pulse.  
         [0079]     Other features and options available for the apparatus and its operating procedure are provide elsewhere within the application.  
         [0080]     According to a third aspect of the invention we provide a sample monitored using the method of the first aspect of the invention and/or using the apparatus of the second aspect of the invention.  
     
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0081]     Various embodiments of the invention will now be described by way of example only and with reference to the accompanying drawings in which:  
         [0082]      FIG. 1  is an illustration of a detection chamber layout;  
         [0083]      FIG. 2  illustrates a distribution assigned to the detection efficiency;  
         [0084]      FIG. 3   a  illustrates the probability distribution for a variable under consideration;  
         [0085]      FIG. 3   b  shows the partial derivative function and a correction vector for the probability distribution of  FIG. 4   a;    
         [0086]      FIG. 4   a  illustrates the probability distribution for an alternative determination;  
         [0087]      FIG. 4   b  shows the partial derivative function for the probability distribution of  FIG. 5   a  illustrating the problem with poorly selected initial values;  
         [0088]      FIG. 5  provides a comparison of NCC and multiplicity counting measurements of a  250 Cf source in various matrix filled drums;  
         [0089]      FIG. 6  illustrates a comparison of NCC and multiplicity counting measurements of a Pu source in various matrix filled drums;  
         [0090]      FIG. 7  provides a comparison of NCC and multiplicity counting limits of detection for a drum containing 50 kg of PVC; and  
         [0091]      FIG. 8  provides a comparison of NCC and multiplicity counting limits of detection for a drum containing 270 kg of steel hulls.  
     
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS  
       [0092]     The detection chamber illustrated in  FIG. 1  comprises a ten sided chamber  2  with openable wall portions  4  to provide access for a drum  6  or the like to be monitored. Four  3 He detectors are provided in each side with twelve similar detectors in both the ceiling and floor  10  of the chamber.  
         [0093]     Drums, 200 litre drums can be accommodated, may be lifted into the chamber, positioned there by forklift truck or conveyed on a bed of rollers. Suspending the drum is preferred to avoid having to strengthen the base of the chamber. Any such bearing surface could interfere with the detection efficiency of the base detectors. Once positioned the chamber is shut.  
         [0094]     The neutron detectors are of conventional type in tube configuration  1   m  long, 2.5 cm in diameter and operated at 4 atmosphere pressure and 960V. The outside of the chamber  2  is provided with an 24 cm thick layer of polyethylene to act as a neutron shield to any background radiation. The detectors are arranged vertically in bores in an inner layer of 8 cm thick polyethylene. A 1 mm layer of cadmium is provided on the inside and outside of the layer to prevent neutrons returning to the chamber and to prevent thermalize neutrons reaching the detectors.  
         [0095]     Each detector is provided with a preamplifier, amplifier and discriminator to minimize dead time effects on the signals. The signals are all summed and then fed to a frequency analyzer.  
         [0096]     The frequency analyzer works on each pulse in the following manner. A received pulse opens after a preset period; an observation interval with the number of pulses falling within this period being counted. A frequency distribution for the number of pulses detected within the observation window is generated as a result. Singlets, doublets, triplets, quadruplets and higher sets can be expected with decreasing occurrences.  
         [0097]     The singlet, doublet and triplet counting rates are determined according to the following procedure.  
         [0098]     The total neutron count rate observed during the measurement, N T , and is obtained simply by counting the number of signals. The first and second factorial moments of the frequency distribution, M m(1)  and N m(2) , are also calculated from the frequency table, using the following equations:  
               N     m   ⁡     (   1   )         =       ∑     x   =   1       x   =     N   max         ⁢     x   ⁢           ⁢     N   x                 (   1   )                 N     m   ⁡     (   2   )         =         ∑     x   =   2       x   =     N   max         ⁢       (         x           2         )     ·     N   x         =       ∑     x   =   2       x   =     M   max         ⁢         x   ⁡     (     x   -   1     )       2     ·     N   x                   (   2   )             
 
 where: 
        N x  is the normalized frequency to have x signals in the observation interval.     N max  is the maximum number of signals observed in any observation interval.        
 
         [0101]     The total count rate and the first and second factorial moments are then background corrected and used in the following equations to calculate the singlet count rate, and the rate of correlated doublets and triplets, R 1 , R 2  and R 3  respectively.  
               R   1     =     N   T             (   3   )                 R   2     =           N   T     ·     N     M   ⁡     (   1   )           -       R   1   2     ·   τ       f             (   4   )                 R   3     =           N   T     ·     N     M   ⁡     (   3   )           -       R   2     ·     R   1     ·   τ   ·     {     f   +     W   2       }       -       1   2     ·     R   1   3     ·     τ   2           f   2               (   5   )             
 
 where: 
        τ=length of the observation interval     l/λ=neutron die away time     T=pre-delay  
             f   =       ⅇ       -   λ     ⁢           ⁢   T       ⁡     (     1   -     ⅇ     1   ⁢           ⁢   λ   ⁢           ⁢   τ         )               (   6   )                 W   2     =     1   -       1     λ   ⁢           ⁢   τ       ⁢     (     1   -     ⅇ       -   λ     ⁢           ⁢   τ         )                 (   7   )             
       
 
         [0105]     The counting rates obtained can then be applied to the following equations which relate the counting rates to the spontaneous fission rate, the self-induced fission rate, the detection efficiency and the α,n reaction rate.  
         R   1     =     ɛ   ·     F   3     ·   M   ·     v     s   ⁢           ⁢   1       ·     (     1   +   α     )           
         R   2     =       ɛ   2     ·     F   s     ·     M   2     ·       v     s   ⁢           ⁢   2       ⁡     (     1   +       (     M   -   1     )     ⁢     (     1   +   α     )     ⁢         v     s   ⁢           ⁢   1       ⁢     v     1   ⁢           ⁢   s             v     s   ⁢           ⁢   2       ⁡     (       v     I   ⁢           ⁢   1       -   1     )             )             
         R   3     =       ɛ   3     ·     F   s     ·     M   3     ·       v     s   ⁢           ⁢   3       ⁡     (     1   +     2   ⁢     (     M   -   1     )     ⁢         v     s   ⁢           ⁢   1       ⁢     v     I   ⁢           ⁢   2             v     s   ⁢           ⁢   2       ⁡     (       v     I   ⁢           ⁢   1       -   1     )           +       (     M   -   1     )     ⁢     (     1   +   α     )     ⁢         v     s   ⁢           ⁢   1       ⁢     v     I   ⁢           ⁢   3             v     s   ⁢           ⁢   3       ⁡     (       v     I   ⁢           ⁢   1       -   1     )         ⁢     (     1   +     2   ⁢     (     M   -   1     )     ⁢       V     I   ⁢           ⁢   2     2         v     I   ⁢           ⁢   3       ⁡     (       v     s   ⁢           ⁢   1       -   1     )             )         )             
 
         [0106]     where 
        ε=detector efficiency     F 3 =the spontaneous fission rate of the sample     M=the leakage multiplication     α=(α,n) to spontaneous fission ratio     n sn =the nth spontaneous fission factorial moment     n ln =the nth induced fission factorial moment     V Sn =the nth spontaneous fission factorial moment (for plutonium)     V ln =the nth induced fission factorial moment (for plutonium)        
 
         [0115]     Rather than fixing one of the unknown values to obtain a solution to the equations the present invention seeks a best fit solution to the information obtained based on certain known absolute constraints and certain probability distributions assigned to each of the unknowns. Information about the mean values and standard deviations for the counting rate are determined experimentally. Information on the distributions for the other parameters can be obtained by calibration experiments as well as being based on experience and past tests.  
         [0116]     The distributions are provided in the following manner:  
         [0000]     Detection Efficiency Distribution 
 
ε˜ N (μ ε , σ ε ) 
 
0≦ε≦1 
 
         [0117]     where 
        μ ε =mean value of the efficiency distribution     σ ε =standard deviation of the efficiency distribution 
 
 Fission Rate Distribution 
       
 
         [0120]     Is constrained to be F ε ≧0 but is otherwise left free floating as the variable to be obtained.  
         [0000]     Multiplication Distribution 
 
 M˜N (μ M σ M ) 
 
 M≧ 1 
 
         [0121]     where 
        μ M =mean value of the multiplication distribution     σ M =standard deviation of the multiplication distribution 
 
 Alpha Distribution 
 
α˜ N (μ x , σ α ) 
 
α≧0 
       
 
         [0124]     where 
        μ α =mean value of the alpha distribution     σ α =standard deviation of the alpha distribution 
 
 Singles Count Rate Distribution 
 
 R   1   ˜N (μ r     1   , σ R     1   ) 
       
 
         [0127]     where 
        μ R     1   =mean value of the singles distribution obtained by measurement     σ R     1   =standard deviation of the singles distribution 
 
 Doubles Count Rate Distribution 
 
 R   2   ˜N (μ R     2   , μ R     2   ) 
       
 
         [0130]     where 
        μ R     2   =mean value of the doubles distribution obtained by measurement     μ R     2   =standard deviation of the doubles distribution 
 
 Triples Count Rate Distribution 
 
 R   3   ˜N (μ R     3   , σ R     3   ) 
       
 
         [0133]     where 
        μ R     3   =mean value of the triples distribution obtained by measurement     σ R     3   =standard deviation of the triples distribution        
 
         [0136]     The probability distribution function (pdf) is assigned a normal distribution with the pdf of any trial value being calculated by  
       z   =     (         y   1     -     μ   y         σ   y       )         
         pdf   ⁡     (     y   i     )       =       1       2   ⁢   π         ·     exp   ⁡     (       -   0.5     ⁢           ⁢     z   2       )             
 
         [0137]      FIG. 2  illustrates such a normal distribution for the detection efficiency. Strictly speaking the pdf  20  should change abruptly to zero when either of the constraints for the efficiency is reached, i.e., greater than or equal to zero and less than or equal to one, however, if such strict limits are applied an adverse effect on the solution process results. As a consequence a constrained pdf  22  is provided but in a form which crosses the constraints to an extent.  
         [0138]     The overall pdf for the variables can be determined by first calculating the count rates and then determining the product of all the individual pdf&#39;s. The solution is the set or values for the unknown parameters which give the maximum pdf product value, calculated according to 
 
 pdf ( e, F   s   , M, α, μ 1 R     1   , μ R     2   , μ R     3   )= pdf (ε)· pdf ( F )· pdf ( M )· pdf (α) ·pdf ( μ   R     1   )· pdf (μ R     2   )· pdf (μ R     3   ) 
 
         [0139]     The method for determining the maximum solution to such a function is to differentiate the function and solve for d(function)/d(parameter)=zero. As the pdf product has four input parameters it needs to be partially differentiated to yield four simultaneous equations which can then be solved. Linearization method or Taylor series method is used to determine corrections (Δε; ΔF s1 ; ΔM; Δα) to be made to the trial value set selected in order to reduce the observed partial derivatives to zero.  
         [0140]      FIG. 3   a  shows a probability distribution  40  for a variable with the partial derivative of this distribution being provided in  FIG. 3   b . This illustrates in a single dimension the correction applied to an initial estimate  42  to give the next estimate  44 . This in turn produces a further correction factor and a subsequent further estimate  46  and so on.  
         [0141]     In four dimensions the linearized equations are  
               ∂     (       ∂   pdf       ∂   ɛ       )         ∂   ɛ       ⁢   Δ   ⁢           ⁢   ɛ     +         ∂     (       ∂   pdf       ∂   ɛ       )         ∂     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         ∂     (       ∂   pdf       ∂   ɛ       )         ∂   M       ⁢   Δ   ⁢           ⁢   M     +         ∂     (       ∂   pdf       ∂   ɛ       )         ∂   α       ⁢   Δ   ⁢           ⁢   α       =     -       ∂   pfd       ∂   ɛ             
               ∂     (       ∂   pdf       ∂     F   s         )         ∂   ɛ       ⁢   Δ   ⁢           ⁢   ɛ     +         ∂     (       ∂   pdf       ∂     F   s         )         ∂     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         ∂     (       ∂   pdf       ∂     F   s         )         ∂   M       ⁢   Δ   ⁢           ⁢   M     +         ∂     (       ∂   pdf       ∂     F   s         )         ∂   α       ⁢   Δ   ⁢           ⁢   α       =     -       ∂   pfd       ∂     F   s               
               ∂     (       ∂   pdf       ∂   M       )         ∂   ɛ       ⁢   Δ   ⁢           ⁢   ɛ     +         ∂     (       ∂   pdf       ∂   M       )         ∂     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         ∂     (       ∂   pdf       ∂   M       )         ∂   M       ⁢   Δ   ⁢           ⁢   M     +         ∂     (       ∂   pdf       ∂   M       )         ∂   α       ⁢   Δ   ⁢           ⁢   α       =     -       ∂   pfd       ∂   M             
               ∂     (       ∂   pdf       ∂   α       )         ∂   ɛ       ⁢   Δ   ⁢           ⁢   ɛ     +         ∂     (       ∂   pdf       ∂   α       )         ∂     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         ∂     (       ∂   pdf       ∂   α       )         ∂   M       ⁢   Δ   ⁢           ⁢   M     +         ∂     (       ∂   pdf       ∂   α       )         ∂   α       ⁢   Δ   ⁢           ⁢   α       =     -       ∂   pfd       ∂   α             
 
         [0142]     In order to simplify the analytical evaluation of these equations, good approximations can be determined from the pdf product equation by observing the effects of small changes in the individual parameters. As an example involving a small change in the detection efficiency, ε, the 1st order derivative  
           (       ∂   pdf       ∂   ɛ       )     ≈     (       δ   ⁢           ⁢   pdf     δɛ     )       =               pdf   ⁡     (       ɛ   +     δ   ⁢           ⁢   ɛ       ,     F   s     ,   M   ,   α   ,     μ     R   1       ,     μ     R   2       ,     μ     R   3         )       -               pdf   ⁡     (     ɛ   ,     F   s     ,   M   ,   α   ,     μ     R   1       ,     μ     R   2       ,     μ     R   3         )               δ   ⁢           ⁢   ɛ           
 
 is arrived at. This expands to give the second order derivative approximation  
             ∂     (       ∂   pdf       ∂   ɛ       )         ∂     F   s         ≈       δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   ɛ       )         δ   ⁢           ⁢     F   s           =       (               δ   ⁢           ⁢     pdf   ⁡     (     ɛ   ,     F             ⁢   s       ,   M   ,   α     )           δ   ⁢           ⁢   ɛ       -                 δ   ⁢           ⁢     pdf   ⁡     (     ɛ   ,       F   s     +     Δ   ⁢           ⁢     F   s         ,   M   ,   α     )           δ   ⁢           ⁢   ɛ             )       δ   ⁢           ⁢     F   s             
 
 with the result that the linearization equations are approximated by  
               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   ɛ       )         δ   ⁢           ⁢   ɛ       ⁢   Δ   ⁢           ⁢   ɛ     +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   ɛ       )         δ   ⁢           ⁢     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   ɛ       )         δ   ⁢           ⁢   M       ⁢   Δ   ⁢           ⁢   M     +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   ɛ       )         δ   ⁢           ⁢   α       ⁢   Δ   ⁢           ⁢   α       =     -       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   ɛ             
               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s         )         δ   ⁢           ⁢   ɛ       ⁢   Δ   ⁢           ⁢   ɛ     +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s         )         δ   ⁢           ⁢     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s         )         δ   ⁢           ⁢   M       ⁢   Δ   ⁢           ⁢   M     +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s         )         δ   ⁢           ⁢   α       ⁢   Δ   ⁢           ⁢   α       =     -       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s               
               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M       )         δ   ⁢           ⁢   ɛ       ⁢   Δ   ⁢           ⁢   ɛ     +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M       )         δ   ⁢           ⁢     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M       )         δ   ⁢           ⁢   M       ⁢   Δ   ⁢           ⁢   M     +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M       )         δ   ⁢           ⁢   α       ⁢   Δ   ⁢           ⁢   α       =     -       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M             
               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   α       )         δ   ⁢           ⁢   ɛ       ⁢   Δ   ⁢           ⁢   ɛ     +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   α       )         δ   ⁢           ⁢     F   s         ⁢   Δ   ⁢           ⁢     F   s       +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   α       )         δ   ⁢           ⁢   M       ⁢   Δ   ⁢           ⁢   M     +         δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   α       )         δ   ⁢           ⁢   α       ⁢   Δ   ⁢           ⁢   α       =     -       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   α             
 
         [0143]     These equations are most readily handled in matrix notation; D2×ΔX=D1 
         Where   ⁢           ⁢   D   ⁢           ⁢   1     =     (           -       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   ɛ                   -       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s                     -       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M                   -       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   α               )         
         D   ⁢           ⁢   2     =     (             δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   ɛ       )         δ   ⁢           ⁢   ɛ               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   ɛ       )         δ   ⁢           ⁢     F   s                 δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   ɛ       )         δ   ⁢           ⁢   M               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   ɛ       )         δ   ⁢           ⁢   α                   δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s         )         δ   ⁢           ⁢   ɛ               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s         )         δ   ⁢           ⁢     F   s                 δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s         )         δ   ⁢           ⁢   M               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢     F   s         )         δ   ⁢           ⁢   α                   δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M       )         δ   ⁢           ⁢   ɛ               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M       )         δ   ⁢           ⁢     F   s                 δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M       )         δ   ⁢           ⁢   M               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   M       )         δ   ⁢           ⁢   α                   δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   α       )         δ   ⁢           ⁢   ɛ               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   α       )         δ   ⁢           ⁢     F   s                 δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   α       )         δ   ⁢           ⁢   M               δ   ⁡     (       δ   ⁢           ⁢   pdf       δ   ⁢           ⁢   α       )         δ   ⁢           ⁢   α             )         
         Δ   ⁢           ⁢   X     =     (           Δ   ⁢           ⁢   ɛ               Δ   ⁢           ⁢     F   s                 Δ   ⁢           ⁢   M               Δ   ⁢           ⁢   α           )         
 
         [0144]     Successive estimates are given by X n+1 =X n +ΔX n  which process is repeated until the correction vector is below a certain predetermined threshold which is accepted as negligible. A final solution:  
         (           ɛ   ^                 F   s     ^               M   ^               α   ^           )     =     X   final         
 
 results. This solution provides a very good fit to the experimental results obtained. From the spontaneous fission rate value so obtained the mass of the neutron emitting isotopes in the sample can be obtained. This in turn can be linked to the overall neutron source mass. 
 
         [0145]     In certain cases due to unusual sample conditions or due to poor selection of the starting criteria the correction vector provided may be such that a true solution cannot be obtained. As illustrated in  FIG. 4   a  with the pdf  50  shown, the differentiated pdf,  FIG. 5   b , may not change sign between the initial estimate  52  and the solution. If the initial estimate  54  is fortunately located then a true solution  56  will be tended towards, but if the initial estimate  52  is not so fortunate a false solutions  58  will be tended towards.  
         [0146]     To counter this trial pdf is evaluated for both positive and negative correction factors.  
         [0147]     The correction vector may substantially overestimate or underestimate the distance to a solution, because the function is non-linear. To counter this, the correction vector for a pdf value is multiplied by 32 and then halved and the pdf value recalculated. This pdf value is compared against the original. If the new pdf value is less than the previous the process is repeated and the new value compared with the new previous value. If the new value is greater than the previous then the correction vector is applied at that value. Otherwise the process is repeated until such a situation is reached.  
         [0148]     To avoid the situation where the initial pdf is zero and as a consequence derivatives cannot be calculated, due to poorly selected starting conditions for instance, the count rate standard deviation is inflated in such a case. Multiplication by a large constant is employed to this end. After the process has converged towards a solution the value of this factor is reduced in stages and the process repeated until a solution is obtained with the inflationary factor set at zero.  
         [0149]     In any experimentally derived result it is important to know the error possible in the result. This is particularly so for the monitoring situations with which this method principally concerned as the implementation of the waste disposal based upon it must always act on the worst possible case in meeting the critical safety factors.  
         [0150]     The associated error estimates are predetermined as follows. In each case the precision of the final solution is dependant on the precision of the count rates for the singlets, doublets and triplets. Assuming these count rates are independent the variance is given by  
           (     σ     ɛ   ^       )     2     =         (       (       ∂     ɛ   ^         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂     ɛ   ^         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂     ɛ   ^         ∂     R   3         )     ⁢     σ     R   3         )     2           
           (     σ       F   ^     s       )     2     =         (       (       ∂       F   ^     s         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂       F   ^     s         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂       F   ^     s         ∂     R   3         )     ⁢     σ     R   3         )     2           
           (     σ     M   ^       )     2     =         (       (       ∂     M   ^         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂     M   ^         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂     M   ^         ∂     R   3         )     ⁢     σ     R   3         )     2           
           (     σ     α   ^       )     2     =         (       (       ∂     α   ^         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂     α   ^         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂     α   ^         ∂     R   3         )     ⁢     σ     R   3         )     2           
 
         [0151]     The partial derivatives again need not be determined but can be approximated by  
           (       ∂     ɛ   ^         ∂     R   1         )     ≈     (       δ   ⁢           ⁢     ɛ   ^         δ   ⁢           ⁢     R   1         )       =           ɛ   ^     ⁡     (         μ     R   1       +     δ   ⁢           ⁢     μ     R   2           ,     μ     R   2       ,     μ     R   3         )       -       ɛ   ^     (           ⁢       μ     R   1       ,     μ     R   2       ,     μ     R   3         )         δ   ⁢           ⁢     μ     R   1               
 
         [0152]     etc. where δμ R     1   =a small change in μ R     1    etc.  
         [0153]     ε(μ R     1   , μ R     2   , μ R     3   )=final estimate from solution to (μ R     1   , μ R     2   , μ R     3   ) count rates set with the partial derivatives being determined for each of the count rates based on the observed rates and small deviations in each.  
         [0154]     The approximate error estimates are thus given by  
           (     σ     ɛ   ^       )     2     =         (       (       ∂     ɛ   ^         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂     ɛ   ^         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂     ɛ   ^         ∂     R   3         )     ⁢     σ     R   3         )     2           
           (     σ       F   ^     s       )     2     =         (       (       ∂       F   ^     s         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂       F   ^     s         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂       F   ^     s         ∂     R   3         )     ⁢     σ     R   3         )     2           
           (     σ     M   ^       )     2     =         (       (       ∂     M   ^         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂     M   ^         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂     M   ^         ∂     R   3         )     ⁢     σ     R   3         )     2           
           (     σ     α   ^       )     2     =         (       (       ∂     α   ^         ∂     R   1         )     ⁢     σ     R   1         )     2     +       (       (       ∂     α   ^         ∂     R   2         )     ⁢     σ     R   2         )     2     +       (       (       ∂     α   ^         ∂     R   3         )     ⁢     σ     R   3         )     2           
 
 which in matrix notation becomes  
       VX   =     DR   ×   VR         
         Where   ⁢           ⁢   VX     =     (             (     σ     ɛ   ^       )     2                 (     σ       F   ^     s       )     2                 (     σ     M   ^       )     2                 (     σ     α   ^       )     2           )         
       VR   =     (             (     σ     R   1       )     2                 (     σ     R   2       )     2                 (     σ     R   3       )     2           )         
       DR   =     (             (       δ   ⁢           ⁢     ɛ   ^         δ   ⁢           ⁢     R   1         )     2             (       δ   ⁢           ⁢     ɛ   ^         δ   ⁢           ⁢     R   2         )     2             (       δ   ⁢           ⁢     ɛ   ^         δ   ⁢           ⁢     R   3         )     2                 (       δ   ⁢           ⁢       F   ^     2         δ   ⁢           ⁢     R   1         )     2             (       δ   ⁢           ⁢       F   ^     2         δ   ⁢           ⁢     R   2         )     2             (       δ   ⁢           ⁢       F   ^     2         δ   ⁢           ⁢     R   3         )     2                 (       δ   ⁢           ⁢     M   ^         δ   ⁢           ⁢     R   1         )     2             (       δ   ⁢           ⁢     M   ^         δ   ⁢           ⁢     R   2         )     2             (       δ   ⁢           ⁢     M   ^         δ   ⁢           ⁢     R   3         )     2                 (       δ   ⁢           ⁢     α   ^         δ   ⁢           ⁢     R   1         )     2             (       δ   ⁢           ⁢     α   ^         δ   ⁢           ⁢     R   2         )     2             (       δ   ⁢           ⁢     α   ^         δ   ⁢           ⁢     R   3         )     2           )         
 
         [0155]     A series of experimental measurements were conducted to test the performance of the maximum-likelihood multiplicity analysis described above. These trials were carried out using a passive neutron counting test rig containing 48 3 He neutron detector tubes encased in polyethylene moderator. The chamber was sufficient to accommodate 200 litre drums and had a neutron detecting efficiency of 15% for an empty drum and a neutron die-away-time of 600 m/s. A variety of matrix filled drums were used in conjunction with Pu and  252 Cf standard sources to simulate waste measurements.  
         [0156]     In the first set of experimental trials, the reduction in systematic errors (largely due to positional and matrix effects) achieved by multiplicity counting compared with NCC using a  252 Cf standard source to simulate a large quantity of plutonium at various positions within waste matrix filled drums was employed. The specified probability distributions for this analysis were (ε, ε 0 , ε + )=(4%, 10%, 16%),(M, M 0 , M + )=(1.00, 1.00, 1.20) and (α, α 0 , α + )=(0.00, 0.00, 1.00). These distributions were designed to represent a measurement scenario where possibly large amounts of plutonium were present within a matrix with highly variable neutron moderating and/or absorption properties.  
         [0157]     The results of these measurements (3000 s count times) are presented in Table 1. The first column of Table 1 gives details of the matrix and the source position for each measurement. The second column shows the measured net (i.e., background corrected) neutron singles, doubles and triples count rates and their associated standard deviations. The third column gives the results of a maximum-likelihood multiplicity analysis of the measured data. The fourth column shows the results that have been obtained from a conventional NCC analysis in which the measured doubles rate is converted directly to a  240 Pu eff  mass by means of a calibration factor. In this case, the calibration factor assumes ε=10% (i.e., the most likely value of the probability distribution specified for the multiplicity analysis).  
                                                     TABLE 1                                   NCC       Measure   Net Count Rates/       Results       Number   ε − 1   Multiplicity Results     240 Pu eff /g                                #1   R 1  = 22773   +/− 3      240 Pu eff  = 132.7 g   100.8   g       PVC (85 kg)   R 2  = 3472   +/− 8    ε = 9.6%       Source @ B   R 3  = 269   +/− 12   M = 1.00                   α = 0.00       #2   R 1  = 20311   +/− 3      240 Pu eff  = 132.0 g   80.6   g       PVC (85 g)   R 2  = 2778   +/− 8    ε = 8.6%       Source @ A   R 3  = 185   +/− 13   M = 1.00                   α = 0.00       #3   R 1  = 31861   +/− 4      240 Pu eff  = 132.2 g   197.6   g       Steel (270 kg)   R 2  = 6806   +/− 12   ε = 13.5%       Source @ B   R 3  = 686   +/− 26   M = 1.00                   α = 0.00       #4   R 1  = 37947   +/− 4      240 Pu eff  = 132.1 g   280.8   g       Steel (270 kg)   R 2  = 9674   +/− 19   ε = 16.0%       Source @ A   R 3  = 1315   +/− 39   M = 1.00                   α = 0.00       #5   R 1  = 28169   +/− 3      240 Pu eff  = 131.9 g   155.1   g       PVC (50 kg)   R 2  = 5344   +/− 12   ε = 11.9%       Source @ B   R 3  = 484   +/− 24   M = 1.00                   α = 0.00       #6   R 1  = 33051   +/− 5      240 Pu eff  = 131.9 g   213.3   g       PVC (50 kg)   R 2  = 7348   +/− 13   ε = 14.0%       Source @ A   R 3  = 831   +/− 28   M = 1.00                   α = 0.00       #7   R 1  = 28105   +/− 5      240 Pu eff  = 132.3 g   153.8   g       Paper (20 kg)   R 2  = 5304   +/− 11   ε = 11.9%       Source @ B   R 3  = 511   +/− 24   M = 1.00                   α = 0.00       #8   R 1  = 35463   +/− 39     240 Pu eff  = 132.9 g   244.1   g       Paper (20 kg)   R 2  = 8410   +/− 22   ε = 14.9%       Source @ A   R 3  = 1091   +/− 52   M = 1.00                   α = 0.00       #9   R 1  = 20604   +/− 10     240 Pu eff  = 132.1 g   82.9   g       Water (105 kg)   R 2  = 2856   +/− 9    ε = 8.7%       Source @ B   R 3  = 181   +/− 13   M = 1.00                   α = 0.00       #10   R 1  = 13345   +/− 2      240 Pu eff  = 131.8 g   34.8   g       Water (105 kg)   R 2  = 1200   +/− 5    ε = 5.7%       Source @ A   R 3  = 51   +/− 6    M = 1.00                   α = 0.00       #11   R 1  = 29895   +/− 4      240 Pu eff  = 131.7 g   174.6   g       Polyethylene   R 2  = 6015   +/− 11   ε = 12.7%       (22 kg)       Source @ B   R 3  = 569   +/− 23   M = 1.00                   α = 0.00       #12   R 1  = 34554   +/− 21     240 Pu eff  = 132.5 g   232.4   g       Polyethylene   R 2  = 8004   +/− 15   ε = 14.6%       (22 kg)       Source @ A   R 3  = 992   +/− 34   M = 1.00                   α = 0.00                  
 
         [0158]     The  240 Pu etf  masses determined by the two techniques are compared graphically in  FIG. 5 . The error bars shown on the multiplicity results have been determined by the maximum-likelihood analysis and include both the systematic and random components. The error bars shown on the MC results also include both random and systematic components. The systematic term has been determined by assuming the same probability distribution for the detection efficiencies were specified for the multiplicity analysis and the random term is simply the standard deviation on the measured doubles rate.  
         [0159]     The consistently accurate results obtained from the maximum-likelihood MC analysis have been obtained by (correctly) varying ε (rather than M and α) to match the predicted and measured count rates. This is true even in extreme cases such as the measurement number 10 where the source was placed at the centre of a water-filled drum.  
         [0160]     In a second set of tests a small Pu source (nominal  240 Pu eff  mass=0.42 g) was measured under similar conditions to establish what benefits the maximum-likelihood MC analysis would provide for the measurement of smaller quantities of Pu. The specified probability distributions for the maximum likelihood MC analysis were (ε, ε 0 , ε + )=(10%, 13%, 16%),(M, M 2 , M + )=(1.00, 1.00, 1.01) and (α, α 0 , α + )=(1.40, 1.70, 2.00). The relatively tight distribution specified for the M reflects the fact that for small quantities of Pu, it is safe to assume that M will be close to unity. The relatively wide distribution specified for a reflects the fact that this parameter will be difficult to estimate for real wastage of variations in chemical composition and impurity content.  
         [0161]     The results of these measurements (3000 s count times) are presented in Table 2 and a comparison of the  240 Pu eff  mass is determined by the maximum-likelihood. MC and NCC techniques are illustrated in  FIG. 6 .  
                                                 TABLE 2                                   NCC       Measure   Net Count Rates/   Multiplicity   Results       Number   ε − 1   Results     240 Pu eff /g                                #1   R 1  = 124.34   +/− 0.31     240 Pu eff  = 0.40 g   0.26 g       PVC (85 kg)   R 2  = 3.88   +/− 0.08   ε = 13.3%       Source @ A   R 3  = 0.32   +/− 0.05   M = 1.00                   α = 1.93       #2   R 1  = 125.93   +/− 0.34     240 Pu eff  = 0.36 g   0.30 g       Wood (48 kg)   R 2  = 4.5   +/− 0.07   ε = 11.8%       Source @ B   R 3  = 0.29   +/− 0.03   M = 1.00                   α = 1.89       #3   R 1  = 156.94   +/− 0.33     240 Pu eff  = 0.41 g   0.47 g       Wood (48 kg)   R 2  = 7.10   +/− 0.10   ε = 13.8%       Source @ A   R 3  = 0.69   +/− 0.10   M = 1.00                   α = 1.70       #4   R 1  = 121.46   +/− 0.33     240 Pu eff  = 0.38 g   0.29 g       PVC (50 kg)   R 2  = 4.41   +/− 0.08   ε = 11.3%       Source @ B   R 3  = 0.48   +/− 0.14   M = 1.00                   α = 1.76       #5   R 1  = 159.82   +/− 0.34     240 Pu eff  = 0.38 g   0.47 g       PVC (50 kg)   R 2  = 7.09   +/− 0.09   ε = 14.2%       Source @ A   R 3  = 0.61   +/− 0.04   M = 1.01                   α = 1.92       #6   R 1  = 137.75   +/− 0.34     240 Pu eff  = 0.37 g   0.37 g       Paper (20 kg)   R 2  = 5.62   +/− 0.07   ε = 13.0%       Source @ B   R 3  = 0.47   +/− 0.05   M = 1.00                   α = 1.79       #7   R 1  = 173.61   +/− 0.35     240 Pu eff  = 0.39 g   0.58 g       Paper (20 kg)   R 2  = 8.77   +/− 0.10   ε = 15.5%       Source @ A   R 3  = 0.86   +/− 0.06   M = 1.00                   α = 1.83       #8   R 1  = 133.44   +/− 0.33     240 Pu eff  = 0.40 g   0.35 g       Polyethylene   R 2  = 5.22   +/− 0.09   ε = 12.0%       (22 kg)       Source @ B   R 3  = 0.66   +/− 0.16   M = 1.00                   α = 1.72       #9   R 1  = 172.47   +/− 0.36     240 Pu eff  = 0.43 g   0.54 g       Polyethylene   R 2  = 8.09   +/− 0.11   ε = 14.3%       (22 kg)       Source @ A   R 3  = 0.84   +/− 0.09   M = 1.00                   α = 1.73                  
 
         [0162]     As can be seen from  FIG. 6 , the maximum-likelihood MC analysis results are again consistently closer to the correct value than the conventional NCC results. This emphasizes the ability of the maximum-likelihood technique to find the correct solution to the multiplicity equations without requiring any of the unknown parameters to be specified exactly.  
         [0163]     The fact that the error bars shown for the MC analysis of the Pu source is significantly larger than those for the  252 Cf measurements reflects the fact that the relative standard deviations on the measured count rates for the far smaller Pu source are significantly greater than those of the  252 Cf source.  
         [0164]     A third set of measurements were also made on drums filled with inactive matrices to evaluate the improvement in the limit of detection that would result from the use of the maximum-likelihood analysis technique compared with NCC. Several sets of hour long measurements were performed on the matrix filled drums and similar measurements were made with an empty chamber to determine the ambient background. The net count rates were determined and then analyzed by both maximum-likelihood MC and conventional NCC techniques. The specified probability distributions for the maximum-likelihood MC analysis were (ε, ε 0 , ε + )=(10%, 13%, 16%), (M, M 0 , M + )=(1.00, 1.00, 1.01) and (α, α 0 , α + )=(1.40, 1.70, 2.00). The comparisons of the results obtained from the two techniques are shown graphically in  FIG. 7  for a drum containing 50 kg of scrap PVC and  FIG. 8  for a drum containing 270 kg of steel hulls. Also shown on the graphs is an indication of the TRU/LLW segregation boundary at 100 nCi/g (calculated by assuming 94% 235 Pu; 6% 240 Pu).  
         [0165]     The NCC results in  FIG. 7  are scattered to either side of zero because they depend solely on the variation in the measured doubles rate from the matrix drum to the empty chain of measurements. The multiplicity results, as discussed earlier, have been constrained to the physically meaningful situation of  240 Pu eff  greater or equal to zero. More specifically, the magnitude and the measurement uncertainty on the multiplicity results is consistently less than half that of the NCC results. This translates directly to a limited protection for the maximum-likelihood MC technique which is less than half that of a conventional NCC analysis performed under the same measurement conditions. This improvement over NCC detection limits is as expected given that more useful information is incorporated into the maximum-likelihood MC analysis.  
         [0166]     The results as shown in  FIG. 8  for the drum containing the steel again show reduced measurement uncertainties for the MC results compared with the NCC results, but now show significant positive bias in the assay results as evident despite the evidence of any spontaneous fissile material. Cosmic ray induced neutrons (which are generated primarily in high Z materials) are the likely cause of this effect. At its current level the bias still allows re-categorization against 100 nCi/g LLW limit. This is despite the lack of background shielding on experimental systems without the use of matrix specific background measurements.  
         [0167]     The improvement in limit of detection shown in  FIGS. 6 and 7  indicates that the maximum likelihood multiplicity analysis technique will provide a consistent and reliable method for the segregation of TRU/LLW at the 100 nCi/g boundary.  
         [0168]     The present invention therefore provides a technique whereby the unknown parameters for a waste sample can be determined with far greater accuracy whilst avoiding undue assumptions about the system. In waste investigations accurate Pu mass measurements can be made over the full range of wastes down to and below the 100 nCi/g LLW boundary. The result is a more accurate and cost effective monitoring procedure with advantages in inventory control and with significant cost savings in organizing waste disposal. The technique allows the low cost and high reliability associated with passive neutron counting systems to be retained. These advantages are obtained despite variations in the make-up of the waste and in the location of the neutron source within that waste. This improved monitoring is backed up by a firm indication of the errors encountered. The parameter envelopes applied to the various functions can be updated as the monitoring of samples progresses. Thus, the results from previous samples can improve the modeling of subsequent samples in the waste stream.