Method of measuring horizontal fluid flow behind casing in subsurface formations with sequential logging for interfering isotope compensation and increased measurement accuracy
CROSS-REFERENCE TO RELATED APPLICATIONS 
     The present application is directed to radioactive logging of wells to measure horizontal fluid flow in formations, as are co-pending U.S. Patent Applications Ser. No. 698,394, filed June 21, 1976 now U.S. Pat. No. 4,051,368, and Ser. No. 808,422, filed June 20, 1977, each of which is assigned to the assignee of the present application. 
     BACKGROUND OF INVENTION 
     1. FIELD OF INVENTION: 
     The present invention relates to radioactive well logging techniques to measure the lateral flow of fluid in subsurface earth formations. 
     2. DESCRIPTION OF PRIOR ART: 
     In secondary and tertiary recovery of petroleum deposits many of the recovery techniques employ the injection of water or chemical solutions into the earth formations comprising the reservoir. In planning the recovery operation, the injection of water or chemical has in the past been limited by certain assumptions and/or approximations concerning the mobility of fluids in the formation comprising the reservoir. A crucial factor in such fluid injection programs is the vertical conformity of the producing formation as well as its horizontal permeability and uniformity. In some reservoirs, formation lensing or horizontal partitioning by permeability barriers such as faults can occur. In such instances, apparently correlevant intervals of permeability may be separated from one well to another in the field by such formation lensing or permeability barriers being interposed across the interval of formation between the wells. 
     It is therefore apparent that large amounts of costly chemicals or water could be injected before it is established that continuity between injection and producing wells is partially or totally absent. If, on the other hand, continuity does exist, it is known that the water in the vicinity of the producing well begins to flow laterally soon after injection is initiated and long before the injected fluid actually arrives at the producing well. The detection of lateral water flow is, therefore, indicative of formation continuity. The speed of the flow, when combined with injection rates, formation thickness, formation porosity, and well spacings, can be used to determine the degree of continuity. This early definition of formation continuity could prevent the expenditure of large sums of money, time and effort in a fruitless project to recover secondary or tertiary problem deposits. 
     A second application of the detection of lateral water-flow is the mapping of the total flow throughout a petroleum reservoir to help in the operational planning of injecting chemicals or water and to assist in determining optimum withdrawal rates. Moreover, a knowledge of the lateral water flow characteristics of a particular formation in a producing field can help greatly in general understanding of the reservoir dynamics of the particular reservoir being produced. 
     It is sometime desirable in a reservoir with multiple producing intervals for a reservoir engineer to be able to delineate those producing zones which provide the most water influx or water drive to the production of petroleum. The mapping of lateral water movement in all zones both above and below the expected water table in the producing formation should supply this information to the reservoir engineer. 
     In the past, reservoir engineers have been provided with relatively few and often inaccurate well logging instrumentation in order to determine the vertical conformance characteristics of the earth formations comprising a reservoir. This has led to resultant confusion as to the properties of the earth formations comprising a reservoir. Radioactive tracer studies of the movement of fluids in the vicinity of a well borehole can be misleading in this respect because of the lack of uniform absorption of the tracer element into the flowing stream of formation water. Also, it is difficult to provide tracer isotopes with sufficient half life to be injected at an injection well and observe their movement days or even weeks later at a monitoring or producing well, in order to obtain some idea of the lateral flow speed or velocity of fluids in the formation comprising the reservoir. 
     In co-pending application Ser. No. 698,394, now U.S. Pat. No. 4,051,368 referenced above, liquids in the formation adjacent the well were bombarded with high energy neutrons. Where the fluid was at least partially saline, as in salt water, radioactive Na 24  was produced by the thermal neutron capture Na 23  (n,γ)Na 24  reaction. By observing the decrease in gamma radiation from the Na 24  and 2.75 MeV or greater energy levels with time, a measure of the horizontal speed of the liquid was obtained. 
     Further, co-pending application Ser. No. 808,422 referenced above, discloses a method of compensating for the effect of interfering gamma radiation of the radioactive isotope manganese 56 resulting from neutron irradiation of the steel casing. 
     However, it has been found that gamma ray activities other than Na 24  from flowing formation water behind casing and Mn 56  from the casing are produced during neutron irradiation of the borehole and casing and formation adjacent thereto. Further, it has been found that these other gamma ray activities must be considered as interfering activities and compensation must be made for their presence in order to obtain more accurate, and thus more meaningful, horizontal water flow velocities in formations. 
     Examples of these other interfering gamma radiation activities typically present in irradiated formations include that of calcium 49 from the activation of calcium within the formation and additional sodium 24 resulting from the activation of interstitial saline formation water or from saline water which is often used to mix the borehole casing cement or saline borehole water. This additional sodium 24 activity is distinguished from the activity due to the horizontal flow of formation fluid in that the radioactive sodium casing the additional activity is not moving with the flowing formation water. It is referred to as &#34;fixed&#34; sodium 24 activity. 
     SUMMARY OF INVENTION 
     Briefly, with the present invention, a method is disclosed for determining the flow rate of earth formation fluid moving in a horizontal direction behind casing using a sequential logging technique to improve the accuracy of flow velocity measurement. The earth formations over an interval of depth in the vicinity of the well borehole at a particular depth are initially logged sequentially, that is a number of times, at a fixed logging speed to determine the number of counts for sub-intervals in the depth interval. From these initial sequential logs, a measure of the background gamma radiation of the interval is obtained. 
     The earth formations in the interval are then irradiated in a sequential number of passes by bombardment with neutrons to neutron activate elements in the formation, formation fluid, casing and borehole. Count rate signals representative of gamma radiation are then detected during a series of sequential counting passes to obtain a cumulative number of counts for each sub-interval in the interval under investigation. The time duration of the intervals during which the count rate signals are obtained is also measured. 
     Based on the detected count rate signals, the measured time intervals, the suspected interfering activities and corresponding half-lives, the flow velocity of the fluid is obtained, as well as a measure of the amount of gamma radiation attributable to the tracer element in the fluid, typically sodium isotope 24, along with a measure of the interfering gamma radiation attributable to elements in the formation and casing. 
     Where desired, the number of additional count rate measurements may be increased, together with the number of measured time intervals to obtain a more statistically precise measure. As another feature of the present invention, the obtaining of a measure of the amount of gamma radiation attributable to elements in the formation, casing and to the tracer element in the fluid lends itself to an iterative process wherein an initial or test flow speed is assigned and initial measures of the amount of gamma radiation attributable to elements in the formation, casing and to the tracer element in the fluid are obtained. The test flow speed is then adjusted based on the results of the initial measures obtained, and subsequent measures of the presence of the elements in the formation, casing and tracer elements in the fluid repeated until a statistically acceptable fluid flow velocity is obtained. 
     Accordingly, with the present invention, the lateral movement of fluids in a well borehole is accurately obtained and determined by neutron activation of the element sodium present in salt water as a portion of the fluids present in the formations adjacent the borehole. Furthermore, the effect of interfering activities of elements in the formation and the casing is taken into account and compensation for the otherwise interfering effect of neutron induced gamma radiation of these elements effected. Thus, according to the present invention, a more statistically accurate and precise measure of lateral movement and horizontal flow speed of formation fluid in the vicinity of a well borehole is obtained. 
     Other objects, features and advantages of the present invention will become apparent to those skilled in the art when considering the following detailed description of the invention. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     FIG. 1 is an illustration showing schematically a well logging sonde for horizontal water flow detection in accordance with the principles of the present invention; 
     FIG. 2 is a graphical representation illustrating measured gamma radiation counting rates in a borehole as a function of time; 
     FIG. 3 is a graphical representation illustrating measured gamma radiation counting rates in a borehole as a function of time and of various horizontal fluid flow velocities; 
     FIG. 4 is a logic flow diagram of process steps suitable for performance in a digital computer according to the present invention; and 
     FIG. 5 is a graphical representation of a measure, obtained in accordance with the present invention, of the amount of gamma radiation attributable to elements in a subsurface formation and to trace elements in moving formation fluid. 
    
    
     DESCRIPTION OF THE PREFERRED EMBODIMENT 
     I. APPARATUS 
     Referring now to FIG. 1, a horizontal flow measuring system in accordance with the present invention is shown schematically. A downhole sonde 10 is shown suspended by a well logging cable 12 in a well borehole 14 which is filled with borehole fluid 16 and surrounded by earth formations 18. As is typical, a steel alloy casing 20 and cement lining 22 are interposed between the formation 18 and the sonde 10. The casing 20 is usually alloy steel, containing manganese as one of the component elements. 
     The well logging cable 12 passes over a sheave wheel 24 which is mechanically or electrically coupled, as indicated by the dotted line 26, to a pulse-height analyzer/recorder 28 so that measurements from the downhole sonde 10 may be recorded as a function of depths in a well borehole 14. A digital computer 30, such as a PDP-11 computer, receives output signals from the pulse height analyzer 28, for reasons to be set forth. 
     Housed in the downhole sonde 10 is, at its lower end, a neutron source 32 which may be a continuous chemical neutron source such as Actinium Berrylium source, an Americium Beryllium source or a Californium 252 source as may be desired. For best results, the neutron source should have an intensity as close as possible to 10 8  neutrons per second. 
     Spaced about five feet from the neutron source is a single gamma ray scintillation detector 34. The detector 34 comprises a sodium iodide thalium activated crystal or a cesium iodide thallium activated crystal approximately 2 inches by 4 inches in extent and cylindrical in shape. The scintillation crystal of detector 34 is optically coupled through a photomultiplier tube (not shown) which functions to count scintillations or light flashes, occurring in the crystal from impingement thereon by high energy gamma rays from radioactive materials in the vicinity. 
     As is known in the art, the pulse height of voltage pulses produced by the photomultiplier of detector 34 are proportional to the energy of the gamma rays impinging upon the crystal of the detector 34. Thus, a succession of pulses from the detector whose pulse height is proportional to the energy of the impinging gamma rays is produced and is coupled to the surface pulse height analyzer 28 via a conductor of the well logging cable 12. Appropriate power sources (not shown) are supplied at the surface and connected to the downhole electronic equipment via conductors of cable 12 in order to supply operational power for the downhole detector 34 in a manner conventional in the art. 
     The space between the neutron source 32 and the detector 34 in the downhole sonde 10 is shielded by a shielding material 36 of suitable type to prevent direct irradiation of the detector crystal with neutrons from the neutron source 32. Shielding materials with high hydrogen content such as paraffin or other poly-molecular hydrocarbon structure may be utilized for this purpose. The high hydrogen content serves to slow down or rapidly attenuate the neutron population from the neutron source and prevent this thermalized neutron population from reaching the vicinity of the detector crystal. To this end, strong thermal neturon absorbers such as cadmium may be interposed in layers with the hydrogenate shielding material in order to make up the shield portion 36. 
     In logging operations, the sonde 10 is moved through a subsurface formation interval D of interest at a specified logging speed so that the interval D is partitioned into a number of sub-intervals, one of which is shown as S in FIG. 1. As will be set forth in more detail below, the sonde 10 first makes a plurality of sequential background logging passes through the interval with source 32 removed and the detector 34 active so that a measure of background gamma radiation from the interval D may be obtained. The source 32 is then inserted within the sonde and the detector 34 deenergized and the sonde 10 makes a plurality of sequential irradiation passes through the interval D during which the formation 18, casing 20 and the tracer element in the formation fluid are bombarded with fast neutrons. Thereafer, source 32 is removed and detector 34 again energized, with the sonde 10 then moved through the interval D for a plurality of sequential detecting passes to detect gamma radiation resulting from radioactive decay of those elements in the formation 18, casing 20 and formation fluid which have become neutron activated during the irradiating passes of the sonde 10. 
     Signals from the downhole detector 33 are transmitted to the surface via the logging cable 12 and are provided as input to the pulse height analyzer/recorder 28. A suitable energy window threshold is set, such as from 2.0 MeV to 3.8 MeV, so that the Na 24  peak at approximately 2.65 MeV is present in the pulse height analyzer/recorder 28, for reasons to be set forth. The computer 30 receives count rate signals from the pulse height analyzer 28 and processes such signals in a manner to be set forth, to determine the lateral horizontal flow velocity v and also the relative interfering activities of formation elements. 
     II. PHYSICAL PRINCIPLES INVOLVED 
     U.S. Patent Application Ser. No. 698,394, now U.S. Pat. No. 4,051,368, referred to hereinabove discloses a technique of measuring v, the linear speed of water moving past a cased well bore at an angle of approximately 90° with respect to the axis of the well bore. The knowledge of v so obtained is very useful in current secondary recovery, future tertiary recovery, and other field operations. 
     According to this technique, it is possible to &#34;manufacture&#34; radioactive isotopes within certain liquids by irradiating the moving liquid with neutrons. As an example, if the formation water is saline, radioactive Na 24  can be produced by the thermal neutron capture Na 23  (n,γ)Na 24  reaction. 
     Accordingly, a logging sonde containing a neutron source is positioned within the well borehole adjacent a formation containing horizontally moving water. The neutron source irradiates the water producing radioactive Na 24  which decays by the emission of gamma radiation. When a gamma ray detector is moved to a position previously activated by the neutron source, a decrease in intensity with time of the induced activity is observed. If the liquid is not moving and radioactive Na 24  is the only source of gamma radiation other than background or natural gamma radiation, the observed decrease in activity with time t, when corrected for background or natural gamma radiation, will follow the exponential decay e - λ Na t  where λ Na  is the decay constant of Na 24 . If, however, the liquid is moving in a horizontal direction, the observed decrease in activity will be due to the exponential decay e - λ Na t  plus an additional decrease caused by the induced activity being swept away from the vicinity of the detector by the moving liquid. The observed decrease in induced activity above the expected exponential decay e - λ Na t  is thus used to determine the horizontal linear speed of the moving liquid. 
     In the other co-pending U.S. Patent Application referenced hereinabove, Ser. No. 808,422, filed June 20, 1977, it was found that radioactive isotope Mn 56 , which occurs due to activation in the steel well casing containing manganese, limits the accuracy of measurement of the flow velocity v. 
     So far as is known, there is no way to adjust energy bias levels of a gamma ray detector and effectively differentiate between Na 24  and Mn 56  gamma radiation. As is evident from FIG. 3, the observed counting rates from Mn 56  are observed even with the energy bias of a gamma ray detector set at 2.65 MeV. Relative contributions from Mn 56  also increase at lower energy biases. One proposed method of eliminating this Mn 56  &#34;interference&#34; would be to delay the counting to allow the shorter lived Mn 56  to decay to a negligible level. In situations where either the linear flow velocity of the fluid is relatively high or the salinity of the water is relatively low, or both, long delays in counting results in loss of the Na 24  gamma radiation as well as that of Mn 56 , prevent meaningful data concerning fluid flow velocity from being obtained. 
     With the techniques of this other co-pending application, the counting rate C(t) in counts/minute recorded by the gamma ray detector at time t measured from the termination of irradiation is 
     
         C(t)=K.sub.Na e.sup.-λ.sbsp.Na.sup.t f(v[t+g(T)])+K.sub.Mn e.sup.-λ.sbsp.Mn.sup.t +B                          (1) 
    
     where 
     λ Na  =0.04651 hrs -1  =decay constant for Na 24   
     λ Mn  =0.2707 hrs -1  =decay constant for Mn 56   
     B=count rate from naturally occurring radioactive elements within the formation 
     K na  =observed gamma ray activity from Na 24  induced by neutron irradiation and measured at the end of irradiation (t=0) if there were no water movement 
     K mn  =observed gamma ray activity from Mn 56  induced by neutron irradiation and measured at the end of irradiation. 
     and further where: 
     
         f(v[t+g(T)])=1+av[t+g(T)]+b(v[t+g(T)]).sup.2+ c(v[t+g(T)]).sup.3 +d(v[t+g(T)]).sup.4                                       (2) 
    
     where for a 10 inch borehole with 7 inch casing in a 33% porosity sand formation, 
     a=-1.43·10 -2   
     b=-5.78·10 -3   
     c=2.25·10 -4   
     d=-4.00·10 -7   
     v=horizontal flow velocity in inches per hour; and 
     
         g(T)=p+qT+rT.sup.2                                         (3) 
    
     where 
     p=0 
     q=0.4553 
     r=0.01067 
     T=irradiation time (hours) 
     Physically, f(v[t+g(T)]) is a term which represents the decrease in observed Na 24  activity due to horizontal water movement after irradiation. G(T) represents a decrease in Na 24  build-up due to water movement during irradiation. As disclosed in this second application, the velocity v is measured while compensating for the presence of Mn 56  and its interfering effect with detection of Na 24  gamma radiation. 
     According to the present invention, the foregoing radioactivity well logging techniques are improved in a new and improved process for measuring v, the quantity of interest, from C(t), a measured nuclear counting rate. Since the nuclear decay process is statistical in nature, C(t) has an associated statistical uncertainty. The uncertainty in C(t) is eventually reflected as an uncertainty in v. It is advantageous, therefore, to induce as much Na 24  activity (K Na ) within the formation water as possible. For a single irradiation pass of the logging sonde 10, K Na  can be increased by increasing the neutron source strength. Considering transportation and handling problems, the maximum practical neutron source strength is approximately 10 8  neutrons/second. 
     On the other hand, K Na  can be increased according to the present invention appreciably by making multiple sequential passes of the logging sonde 10 over the interval D of interest. Defining L I  as the logging speed and H Na  as the half life of Na 24 , it is apparent that K Na  will build up considerably as a result of sequential irradiation passes if D/L I  &lt;H Na . Unfortunately, the interfering radiation from Mn 56 , K Mn , also builds up. However, the build up rate of K Mn  is less than that of K Na  since H Mn  &lt;H Na . 
     It can be shown, and verified, from test results obtained using test cell experimental data, that ##EQU1## where λ Na , λ Mn  =decay constants of Na 24  and Mn 56 , respectively (hrs -1 ); 
     L i  =logging speed during irradiation phase (feet/minute); 
     D=vertical extent of interval being investigated (feet); 
     Z mn , Z Na  =constants depending upon reaction cross sections, sonde design, formation porosity, water salinity, and borehole conditions which are obtained from and verified by using test well data; and 
     N i  =number of irradiation passes. 
     Equations (4) and (5) assume (a) the effective vertical extent of irradiation at a given depth, Q, is one foot, and (b) the irradiation passes over the vertical interval D are made sequentially in time (since Q equals one foot, it will be ignored in the remaining Equations). Substituting Equations (4) and (5) into Equation (1) yields ##EQU2## 
     At this point, the amount of Na 24  (and also the amount of Mn 56 ) has been increased in the vicinity of the borehole by utilizing the technique of sequential irradiation. In principle, three detecting or counting passes can now be made with the logging sonde 10 at three different times t 1 , t 2 , and t 3  ; record the counting rates C(t 1 ), C(t 2 ), and C(t 3 ), obtain three independent equations of the form of Equation (6); and solve, in the manner set forth in co-pending U.S. Patent Application Ser. No. 808,422, these three Equations. Of course, B must be measured previously, in a manner to be set forth. 
     Solution of these equations yields results for Z Na , Z Mn , and f(v(t+g[T]). From the resulting value of f(v(t+g(T))), Equation (2) may then be solved to obtain v, the quantity of interest. 
     However, with the present invention it has been found that the statistical accuracy of v can be improved further by making a series of sequential counting passes following a series of sequential irradiation passes. Once the entire interval D has been activated sequentially, the neutron source 32 is removed or de-activated, and the counting phase of the logging operation begins. Considering now the response of the gamma ray detector 34 within a subinterval of vertical extent S (FIG. 1), a number N c  of sequential counting passes during a time interval t a  to t b  are made. I(t a ,t b ), the cummulative number of counts recorded in the interval S during the time interval t a  to t b  (see FIG. 2), is ##EQU3## where u j  and r j  are the times, measured from the end of the irradiation sequence, at which the detector 34 enters and leaves zone S, respectively. The subscripts j indicate the counting pass number. Referring again to C(t) as defined in Equation (6), it can be seen that Equation (7) is a sum of time integrals including fourth-order polynomial functions of time, f(v[t+g(T)]) as defined in Equation (2), which is the function containing the flow velocity v, which is the quantity of interest, as a portion thereof. However, in order to enhance the statistical accuracy of the measurement of I(t a ,t b ), logging passes with the sonde 10 are made sufficiently often that 
     
         t.sup.a -t.sup.b &lt;&lt;0.693/λ.sub.Na                   (8) 
    
     and the time of logging delayed so that, 
     
         K.sub.Mn e.sup.-λ.sbsp.Mn.sup.t &lt;&lt;K.sub.Na e.sup.-λ.sbsp.Na.sup.t f(v[t+g(T)]), then          (9) ##EQU4## where ##EQU5## is the average count rate recorded during the interval t.sup.a to t.sup.b (FIG. 2). The uncertainty in Equation (11) is the percent standard deviation. The assumptions of Equations (8) and (9) are not necessary conditions for the sequential logging operation but will greatly simplify the discussion of the statistical accuracy of the measurements. Field observations, however, have shown that assumptions of Equations (8) and (9) are usually valid. 
    
     III. BACKGROUND GAMMA RADIATION MEASUREMENT 
     Before initiating the irradiation phase of the logging operation, the background counting rate B is determined. This is accomplished by logging sequentially the entire depth interval D with the neutron source 32 removed or inactive. N B  sequential passes are made at a logging speed L B . I B , the cumulative number of counts from a subinterval S, is recorded. I B  is related to B through the Equation 
     
         I.sub.B =BN.sub.B S/L.sub.B , or                           (12) 
    
     
         B=I.sub.B L.sub.B /N.sub.B S+/-(L.sub.B /BN.sub.B S).sup.1/2 (13) 
    
     where again the uncertainty is the percent standard deviation. Substituting Equations (13) and (11) into Equation (6) yields ##EQU6## The terms on the left hand side of Equation (14) are either measured (I(t a ,t b ), I B ) or are known (L B ,L C ,N B ,N C ,S) with a percent standard deviation of ##EQU7## Again, using techniques described in co-pending Patent Application Ser. No. 808,422, three groups of sequential counting passes are made during the time intervals t a   j  to t b   j  (j=1,2,3), the corresponding values I(t a   j ,t b   j ) are recorded, three equations of the form (14) are solved for f(v[t+g(T)]), Z Mn , and Z Na , and v is obtained from f(v[t+g(T)]). Examining (14-a), it can be seen that the percent standard deviation of the measured quantities is decreased by increasing N B ,N C ,S and/or decreasing L B  and L C . 
     IV. EXAMPLE LOGGING OPERATION 
     In logging an example interval D (FIG. 1), 50 feet thick, pertinent information concerning the logging sonde 10 is tabulated in FIG. 3. Ten sequential background logging passes are made at 3 feet/minute (N B  =10, L B  =3). After delays of 12, 14, and 24 hours, the interval D is logged with 12 sequential passes each at three feet per minute t a   1  =12, N C ,1 =12, L.sbsb.C 1  =3; t a   2  =14, N C ,2 =12, L c ,2 =3; t a   3  =24, N C ,3 =12, L C ,3 =3. The counting sequence is shown graphically in FIG. 3. The counts per minute curves are based upon actual field data normalized to the specified irradiation and count sequences. Typical values for the left hand side of Equation (14) along with standard deviations calculated with Equation (14-a) are shown in FIG. 3 for v=0. It is apparent from the error bars that horizontal water flows less than 1 foot/day can be detected with the sequential logging technique of this invention. This technique is thus effective in scanning relatively thick intervals D for horizontal water flows of this magnitude or greater. 
     V. COMPENSATION FOR FORMATION ELEMENT NEUTRON ACTIVATED GAMMA RADIATION 
     As disclosed hereinabove, the total gamma radiation intensity C&#34;(t i ) present within the borehole at a time t i  measured from the end of neutron irradiation can be expressed as ##EQU8## where K Na  =intensity of moving Na 24  present at the end of irradiation 
     λ Na  =decay constant of Na 24  ; 
     v=the horizontal flow velocity of the water; 
     T=irradiation time; 
     K j  =intensity of the jth interfering activity at the end of irradiation; 
     λ j  =the decay constant of the jth interfering activity; 
     n=the total number of interfering activities; and the functions f(v·[t i  +g(T)]) and g(T) are defined in Equations 2 and 3 above, respectively. 
     If (n+2) counting rates C(t i ) are measured at times t i , (n+2) equations of the form of Equation (15) are obtained. It is possible, therefore, to solve this system of (n+2) equations for the (n+2) unknown quantities of v, K Na , and K j  (j=1 to n). In practice, however, the direct solution of the system of equations is extremely complex due to the f function in Equation (2) which is a fourth order polynomial in v. For situations such as this, iteration techniques are preferably used to obtain the desired solution. 
     According to the present invention, therefore, the computer 30 operates according to a sequence of steps (FIG. 4) to obtain accurate flow velocities at the well site. For a given value of v, the unknown quantities K Na , K j  (j=1 to n) are determined by least squares fitting the n+1 Equations (15) to the measured counting rates C(t i ). For each iterated value of v, the quantity ##EQU9## is computed where m is the number of measured counting rates. χ 2  approaches a minimum as the iterated value of v approaches the true horizontal flow velocity. The minimization of χ 2  is, therefore, used as the criterion for determining the &#34;goodness of fit&#34; of the iterated velocity v. 
     The user first enters the irradiation time and background count rate during process step 50. Next, the measured counting rates C(t i ) and the corresponding times t i  are entered into the computer during step 52. The number of suspected interfering activities n and the corresponding half lives are then entered during step 54. Provisions for editing the count data and the interfering reaction parameters may be made, if desired. Finally, an initial assumed value of the velocity v and the velocity iteration parameters of the iteration increments and number of increments to be made are entered during step 56. For each velocity value, the terms K Na  e - λ.sbsp.Na t  f(v[t i  +g(T)]); K j  e - λ j t  i; C&#39;(t i )--computed count rate; and χ 2  are computed and printed during step 58. 
     A decision instruction 60 determines whether an optimum least squares fit is had between the computed count rate C&#39;(t i ), based on the iteration velocity v, and the measured count rate C(t i ), by determining whether the quantity χ 2  is a minimum for the present iteration velocity v. If not, the velocity v is incremented in the amount specified by the iteration parameters during process step 62 and process step 58 repeated and a new decision step 60 made. These steps are repeated for the iteration cycle. After the iteration cycle is completed, the best fit value of v along with K Na  and K j  (j=1 to n) are printed during process step 64. 
     Upon examination of the fit parameters, it is possible to again edit the raw data, vary the iteration value of v, vary the number and half lives of the interfering reactions, or subtract specified amounts of interfering activities in order to improve the fit. 
     FIG. 5 is a graphical representation of results obtained with the present invention from field data and processed using the techniques of the present invention. As can be seen, the interfering effect of Ca 49  and Na 24  (fixed) have been determined so that compensation can be made therefor. 
     The foregoing disclosure and description of the invention are illustrative and explanatory thereof, and various changes in the size, shape and materials as well as in the details of the illustrated construction can be made without departing from the spirit of the invention.