Title of Invention

A METHOD FOR EVALUATING AN INFLOW PERFORMANCE FOR COMPLETED INTERVALS IN A WELL USING A DOWNHOLE TOOL

Abstract A method for evaluating an inflow performance for completed intervals in a well includes pulsing neutrons into a fluid flow in the wellbore, detecting gamma rays from decays of nitrogen-16 in the well fluid flow with a first gamma ray detector, determining an in situ water velocity from the measurement of the gamma ray decays, and estimating the inflow performance of one or more competed intervals in the well from data that includes the in situ water velocity. H-Final Draft (3-31)
Full Text PATENT
201. One or more packers 203 may be used to isolate the interval containing the production zones 111, 112,113 inthe multi-layer reservoir and to stabilize the production tubing 201 in the well 101. In some cases, this is done because other production zones exist above or below the production zones 111,112,113 shown in FIG. 2. The packer 203 seals the annulus between the tubing and the casing so that only the formation fluids from the desired production zones 111,112,113 may flow through the tubing 201. Additional strings of tubing may be installed for producing from other intervals. [0006] After the well has been completed, production logging may be used to obtain in situ measurements of the nature and behavior of the fluids in or around the wellbore during production. The purpose of production logging is to analyze the flowing performance of a well, including the productivity of different production zones in the well. This enables the well operator to produce and obtain maximum production from the well in the most economical way possible. [0007] FIG. 3 shows one example of how production well logging may be performed. The tubing string 201 in positioned such that the lower end of the tubing string is located above the uppermost production zone 111 in the reservoir. A production logging tool 301 is lowered into the well 101 and data are taken as the tool 301 is moved through a vertical region where the measurements are to be taken. In FIG. 3, the production logging tool 301 is positioned between the two uppermost production zones 111, 112, where it may be used to measure flow rates, pressure, temperature, fluid mixture density, and fractions of liquid and gas, among other parameters. By positioning the production logging tool 301 in the position shown in FIG. 3, it will measure the cumulative flow that results from the inflow from production zones 112 and 113. The production logging tool 301 is shown schematically, and the tool 301 in FIG. 3 does not represent the shape or relative size of an actual production logging tool. There are multiple types of production logging tools known in the art, such as spinners, pitot tubes, and turbine meters, that are typically used in combination with other measurement probes, including optical scanning devices to measure the fractions of gas and liquid in well stream.
[0008] Using a production logging tool 301, the flow characteristics of the well fluid stream may be measured and recorded (or transmitted to the surface) at multiple points. For example, the production logging tool 301 may subsequently be moved to a position between the second and third production zones 112,113. In such a position, the production logging tool 301 would measure and record (or transmit to the surface) data based on the inflow from the lowest production zone 113. In

PATENT
this manner, the inflow performance may be determined for each of the production zones 111,112, 113 by measuring the composite flow characteristics at multiple locations. [0009] It is noted that a multi-layer interval may contain more than two or three production zones. A production logging tool may be relocated between each zone so that data corresponding to the properties of the inflow from each zone may be determined.
Summary of the Invention
[0010] In one aspect, the invention relates to a method for evaluating an inflow performance for completed intervals in a well that includes pulsing neutrons into a fluid flow in the wellbore, detecting gamma rays from decays of nitrogen-16 in the well fluid flow with a first gamma ray detector, determining an in situ water velocity from the measurement of the gamma ray decays, and estimating the inflow performance of one or more competed intervals in the well from data that includes the in situ water velocity. In some embodiments, determining the in situ water velocity may include determining a time for activated oxygen to flow from a pulsed neutron generator to the first gamma ray sensor and dividing a distance between the pulsed neutron generator and the gamma ray sensor by the time for the activated oxygen to flow from the pulsed neutron generator to the first gamma ray sensor.
[0011 ] In another aspect, the invention relates to a method for profiling multiphase fluid flow in a well that includes pulsing a first pulse of neutrons into the fluid flow, detecting gamma rays from inelastic collisions between the first pulsed neutrons and a nuclei in the fluid flow, determining an inelastic count rate from the detected gamma rays from the inelastic collisions, pulsing a second pulse of neutrons into the flow, detecting gamma rays from decays of nitrogen-16 in the flow, and determining a water velocity in the production fluid.
[0012] In another aspect, the invention relates to a downhole tool that includes a pulsed neutron generator and at least one gamma ray detector, wherein the downhole tool is configured to measure a velocity of water in a well fluid.
[0013] In another aspect, the invention relates to a method for profiling multiphase fluid flow in a well that includes pulsing a first pulse of neutrons into the fluid flow at a plurality of locations, detecting gamma rays from inelastic collisions between the first pulsed neutrons and a nuclei in the fluid flow at the plurality of locations, determining an inelastic count rate of the production from the

PATENT
detected gamma rays from the inelastic collisions at the plurality of locations, determining a liquid holdup at the plurality of locations based on the inelastic count rate at the plurality of locations, pulsing a second pulse of neutrons into the flow at the plurality of locations, detecting gamma rays from decays of nitrogen-16 in the flow at the plurality of locations, determining a water velocity in the production fluid at the plurality of locations, and determining the inflow performance of each of a plurality of production in the well based on the liquid holdup and the water velocity at the plurality of locations.
[0014] In another aspect, the invention relates to a method for evaluating an inflow performance for completed intervals in a well that includes pulsing neutrons into a fluid flow in the wellbore, detecting gamma rays from inelastic collisions between the pulsed neutrons and a nuclei in the fluid flow, determining an inelastic count rate of the fluid from the detected gamma rays from the inelastic collisions, determining a liquid holdup of the fluid flow from the inelastic count rate, and determining the inflow performance of one or more competed intervals in the well from data that includes the liquid holdup.
[0015] In another aspect, the invention relates to a method for correcting an in situ gas phase velocity approximation that includes estimating a critical gas phase velocity. When the critical gas phase velocity is greater than the in situ average gas phase velocity, the method may also include determining a flow regime, estimating a bubble rise velocity, and computing a corrected in situ gas phase velocity.
Brief Description of the Drawings
[0016] FIG. 1 shows a cross section of a typical well with a multilayer reservoir having multiple
production zones.
[0017] FIG. 2 shows a cross section of a well with a tubing string positioned in the well for
transporting well fluids to the surface.
[0018] FIG. 3 shows a prior art production logging tool for characterizing the flow properties of
a multilayer reservoir.
[0019] FIG. 4 shows a cross section of a well with a velocity string for unloading the well.
[0020] FIG. 5 A shows a cross section of a well with a pulsed neutron logging tool disposed
within the velocity string.

PATENT
[0021] FIG. 5B shows a cross section of a well with a pulsed neutron logging tool disposed
within the velocity string.
[0022] FIG. 6 shows a cross section of an example of a downhole tool that includes a pulsed
neutron logging tool.
[0023] FIG. 7 A shows a graph of the predicted count rate measurement from activated oxygen
versus an actual count rate for flowing activated oxygen.
[0024] FIG. 7B shows a graph of an actual count rate for activated oxygen.
[0025] FIG. 8 shows a graph comparing results from a pulsed neutron log analysis and from
conventional production log analysis.
[0026] FIG. 9 shows one example of a method for evaluating inflow performance for completed
intervals in a well.
[0027] FIG. 10 shows one example of a method for profiling multiphase production flow in a
well.
[0028] FIG. 11 shows one example of a method for evaluating inflow performance for completed
intervals in a well.
[0029] FIG. 12 show one example of a graph of the coefficient of friction on a spherical liquid
particle based on the Reynolds number.
[0030] FIG. 13 shows one example of a method for correcting an in situ gas phase velocity.
Detailed Description
[0031] The present invention is generally related to using a pulsed neutron logging tool to determine the water velocity in a well fluid flow, as well as the inelastic count rate, which may also be measured. This information may then be used in a multiphase fluid flow production or injection systems analysis to determine the inflow performance of the completed intervals in a multilayer reservoir.
[0032] FIG. 4 shows a cross section of a production tubing or velocity string 404 positioned in well 401 having multiple production zones 411,412, 413. A packer 403 may be used to seal the production zones 411,412, 413 from the wellbore above and to force the produced fluids to flow into the velocity string 404. Typically, a velocity string 404 is used to maintain the wellbore fluid outflow velocity at a sufficient level so that formation produced or wellbore condensed liquids do not

PATENT
accumulate or increase in quantity in the wellbore 401 such that the reservoir energy would not be sufficient to produce those fluids to the surface. The velocity string 404 has a smaller diameter than the casing 402 or other production tubing (not shown) that may be used, and the smaller diameter causes the fluid at a given flow rate to flow at a higher velocity. The higher velocity may more efficiently lift the liquids in the wellbore 401 to the surface to minimize or prevent their accumulation in the well 401 which may result in an increase in hydrostatic pressure acting against the productive formations or completed intervals 411,412,413 in the well, resulting in "loading up" of the well and its ceasing to flow naturally. Using a velocity string 404 in this manner is often called "unloading" a well because the velocity string 404 enables fluids to flow from a well 401 that may otherwise be "loaded up" with fluids. A velocity string is generally described, but should not be considered limiting. Other types of tubulars may be used, and in some examples, the described methods and tools may be used within a casing, outside of any tubulars.
[0033] In some cases, a velocity string 404 must be located below the top of the upper most production zone 411 or completed interval. This is so that the fluids flowing into the well 401 must flow at a higher velocity starting at a lower position. Thus, as shown in FIG. 4, in some cases the velocity string 404 may even be positioned such that the bottom of the velocity string 404 is below the lowest production zone 413 in a multilayer reservoir.
[0034] Use of conventional production logging tools in the well 401 shown in FIG. 4 is not effective because the measurement devices, when positioned in the velocity string 404, would measure only the composite flow from the inflow from all of the production zones 411, 412,413, and would be unable to measure fluid properties outside of the tubing. Moreover, it may not be possible to move the velocity string 404 vertically and then take measurements at selected locations. In many cases, this would cause the well 401 to become loaded up with fluids. In other cases, while there may still be a production flow, the inflow performance at each production zone 411,412,413 will vary based on the position of the velocity string 404. Thus, such a technique will not produce accurate results of inflow rates from the completed intervals 411,412,413 that are representative of their inflow contributions to the composite well production when the velocity string 404 is in a different vertical position.
[0035] FIG. 5 A shows a cross section of a well 401 with a pulsed neutron logging tool 501 is suspended from a wireline 502 and positioned within the velocity string 404. The inflow from each

PATENT
of thd production zones 411,412,413 flows downwardly to the end of the velocity string 404, then the composite flow travels up the velocity string 404. A packer 403 forces the production fluids from the production zones 411,412,413 to flow into the velocity string 404. It is noted that a packer may not be necessary if the hydrostatic pressure in the annulus between the velocity string 404 and the casing 402 is high enough to prevent the production fluids from flowing upwardly in the annulus or if upward migration of the produced fluids in the annulus is prevented in some other manner, such as no flow permitted from the annulus at the surface by closing the annular valves. [0036] The pulsed neutron logging tool 501, as will be described, may be used to measure the inelastic count rate of the fluids in the annulus between the velocity string 404 and the casing 402, as well as the water velocity in the fluid flow. Based on this information, collected at a plurality of positions, the inflow performance of each production zone 411,412, 413 may be estimated. [0037] Use of a pulsed neutron logging tool requires that there is at least some water in the fluid flow to be analyzed. This requirement is easily met. Even if there is not any water in the produced fluid stream, a low-rate pump may be used to introduce a small amount of water into the fluid flow, as is known in the art.
[0038] FIG. 5A, the pulsed neutron logging tool 501 is shown positioned between production zones 411 and 412. In this position, the inelastic count rate and water velocity measurements will represent the flow in the annulus between these production zones 411,412. Because these are the two upper-most zones, the measurements made at the position shown in FIG. 5 A will represent the inflow performance of the first production zone 411.
[0039] FIG. 5B shows the same well 401 as in FIG. 5A, but with the neutron logging tool 501 positioned between the second and third production zones 412, 413. In this position, the measurements will represent a composite flow of the inflow from both the first and second production zones 411, 412. The inflow performance of the second production zone 412 may be determined by subtracting the data collected relating to only the first production zone 411. Similar measurements and calculations may be performed over the vertical depth of the well 401 and between each of the production zones. Wells with more than two or three production zones may also be characterized in this manner, possibly requiring more measurements to characterize the additional zones. In some cases, the data collected from a pulsed neutron production log are used in a complete

PATENT
production systems analysis computational model to determine the fluid flow characteristics at a plurality of points in the well. The inflow performance may be determined from the model. [0040] Reference will be made to FIG. 6 to describe the neutron interactions and detector arrangements. FIG. 6 shows a schematic of a downhole tool 610 that is configured to operate as a pulsed neutron logging tool. The downhole tool 610 is positioned in velocity tubing 604 in a well 601. Production fluids, which may be a gas, a liquid, or both, flow downwardly through the annulus between the velocity tubing 604 and the casing 602, shown by the downward flow arrows 612. Inside the velocity tubing 604, the production fluids flow in an upward direction, shown by upward flow arrows 613. As noted above, the invention is not limited to velocity tubing. This is provided only as an example.
[0041] The flow 613 on the inside of the velocity tubing 604 may be a composite flow from all of the production zones in a multilayer reservoir. The flow 612 in the annulus of the well 601 may be comprised of the inflow from one or more of the production zones in a multilayer reservoir, depending on the position of the downhole tool 610 relative to the production zones, as discussed above with reference to FIGS. 5A and 5B. In addition, it is noted that even though this discussion references upward and downward flow, a downhole tool may be located in a horizontal or deviated well. The same principles will apply to a horizontally positioned tool, and the invention is not limited by the orientation of the tool and well. Up and down are used in this disclosure for convenience. Also, a pulsed neutron logging tool may be used to measure the properties of an upward flow. For example, production fluids generally flow upward in a tubing. In an injection well, the injection flow may flow downward through a tubing, but upward in the annulus. Other flow scenarios will be explained later.
[0042] The downhole tool 610, configured to operate as a pulsed neutron logging tool, is shown as a rough schematic diagram. Pulsed neutron logging tools are known in the art. Generally, a downhole tool is a modular tool that may include any number of sensors and detectors for a particular downhole application. In order to perform pulsed neutron logging functions, a downhole tool may include a pulsed neutron generator 621 along with appropriate gamma ray sensors 624,625, 628. Such a tool 610 may also be configured to perform in a number of different modes to gather different types of information. The gamma ray sensors 624,625,628 may also be used for different functions within the downhole tool, such as gamma ray spectroscopy. A pulsed neutron logging tool

PATENT
may be a downhole tool operating in a pulsed neutron logging mode. It is noted that the downhole tool 610 in FIG. 6 is large relative to the components that serve the pulsed neutron logging functions. This is because, as stated above, a typical down hole tool is modular and will include other components for other logging functions.
[0043] The downhole tool 610 shown in FIG. 6 includes a pulsed neutron generator 621 and several gamma ray detectors 624, 625,628. The first two detectors 624,625 typically are called the near detector 624 and the far detector 625. In a typical tool, the near detector 624 is located less than 1 foot from the pulsed neutron generator 621, and the far detector is located on the order of about 2 feet away from the pulsed neutron generator 621. The gamma ray detector 628 is located farther away from the pulsed neutron source 621. Depending on the application, it may be located 20 feet from the pulsed neutron generator, or more. The third detector 628 may be a gamma ray detector that is typically used for other tool functions, but is able to operate in conjunction with the pulsed neutron logging tool functions for purposes of determining the inflow performance of the production zones in a multilayer reservoir. In other cases, the third gamma ray sensor 628 may be positioned in the tool specifically for purposes of determining the inflow performance of the production zones in a multilayer reservoir.
[0044] One possible type of pulsed neutron generator is a particle accelerator that may be energized when a pulse of neutrons is desired. Such a neutron generator is generally not directional, meaning that the neutrons are emitted spherically in all directions. The gamma ray sensors are also non-directional sensors positioned to detect gamma rays coming from all directions. Other types of pulsed neutron generators and gamma ray sensors may be used.
[0045] In one example, a pulsed neutron logging tool may be used to determine the inelastic count rate, which provides an indication of the gas and liquid fractions in the flow stream. The inelastic count rate may be determined by pulsing the pulsed neutron generator 610 and sensing the gamma rays that result from inelastic collisions between the neutrons and nuclei of atoms. Generally, the gamma rays from inelastic collisions are sensed in the first several milliseconds after the pulse. In that time period, the neutrons are generally still fast neutrons because they have not had enough time to slow down to the point where neutron capture is significant. In one example, the inelastic count rate is determined from measurements by the far gamma ray detector 625.

PATENT
IUU46J The inelastic count rate and the liquid holdup may be determined from the gamma ray data using any means known in the art. For example, the liquid holdup in the flow stream may be directly computed from the inelastic count rate (Irat) values recorded in a pulsed neutron log. The minimum and maximum inelastic count rate values represent the case where the fluid is 100% water (Iratw) and 100% gas (IratgX respectively. These set points may be established by an interpretation engineer. The minimum set point for the analysis {Iratw) can often be readily obtained if there is a standing water column in the well, such as at the bottom of the well below the deepest completed interval. The maximum set point {Iratg) is generally more difficult to estimate, except in cases where the well production rate is sufficiently low enough that the well stream liquid falls back and essentially single phase gas is produced from the well at the surface. In such cases, the 100% gas point in the well flow stream can be observed in the log response. In most all other cases, the maximum set point (Iratg) must be reasonably estimated by an interpretation engineer. [0047] With the minimum (Iratw) and maximum (Iratg) set points established for the inelastic count rate analysis, the liquid holdup in the well (Yw) may be computed directly from the recorded inelastic count rate values at the logged depths {Irat(z))in the well using the relationship given in Equation 1:

[0048] The relationship in Equation 1 is only one example of a relationship that may be used to determine the liquid holdup based on the inelastic count rate. Other relationships may be used. [0049] The "liquid holdup" is the fraction of liquid that is present in an interval of pipe or casing. Because the liquid and gas in a multiphase flow will have different densities and viscosities, the different fluids will flow at different velocities, with the heavier or more viscous liquid phase moving slower, or being held up more.
[0050] In one example, a pulsed neutron logging tool may be used to determine the velocity of water flowing in the production fluid. This process will be explained using the downhole tool 610 in FIG. 6 as a reference. The pulsed neutron generator 621 is pulsed, and neutrons are incident on the fluid near the pulsed neutron generator 621. The neutrons will collide with nuclei in the fluid in an

PATENT
inelastic collision. Some of the energy from the neutron will be imparted to the nuclei. Collisions of neutrons with oxygen nuclei will cause a transmutation reaction. The neutron and an oxygen-16 nucleus collide to form a nitrogen-16 nucleus and a proton. This is called oxygen activation, because the stable oxygen nucleus is transmuted in to a different, radioactive nucleus. The general reaction equation is shown in Equation 2: [0051]
[0052] Nitrogen-16 has a half-life of 7.13 seconds, and when it decays, nitrogen-16 will most often decay by beta emission, whereby the nitrogen-16 nucleus emits a beta particle to become an oxygen-16 nucleus. In addition, the beta decay of nitrogen-16 also yields a very high-energy gamma ray, about 6.129 MeV. The gamma ray may be detected by gamma ray detectors. [0053] FIGS. 7A and 7B include graphs that illustrate how the oxygen activation may be used to determine the velocity of water in the production fluid flow. The graph in FIG. 7A includes a line 701 representing one example of the normalized predicted number of gamma rays that will be detected at the far sensor (625 in FIG. 6) as a result of oxygen activation versus time, if the fluid is not moving. The predicted number of counts decays exponentially as the nitrogen-16 decays away. The second line 702 represents one example of the actual counts that are detected by the far sensor. The second line 702 rises above the predicted line 701 because the activated water is moving towards the sensor, and the counts will increase according to the inverse square law. At the point 703, the detected number of counts reaches a local maximum. Following that, the second line 702 falls quickly, eventually falling below the predicted line 701.
[0054] The local maximum count rate 703 observed in the second line 702 represents the point where the activated oxygen (i.e., the nitrogen-16) is closest to the far detector. In FIG. 7A, that occurs at approximately 8 seconds after the pulse, which is taken to be at time 0. In order to determine the velocity, the distance between the far detector 625 and the pulsed neutron source 621 is divided by the time to the maximum — the time it took the activated oxygen to travel the distance between the pulsed neutron source 621 and the far detector 625.
[0055] FIG. 7B shows another example of an actual count rate 712 based on the third gamma ray detector 628 in FIG. 6. The local peak 713 in the actual count rate 712 occurs at about 17 seconds



PATENT
[0060] The water velocity required for the analysis of the inflow from each of the completed intervals in the well begins by the computation of the water velocity moving past the pulsed neutron logging tool, as described above. The computed water velocity must account for the effect of the logging tool in the well. In the case where the water velocity being measured is in the annulus between the tubing and casing, no correction is required. Where the water velocity being considered is the result of tubular flow (i.e., flow and tool are in the tubing or both are in the casing below the end of the tubing), however, the velocity reported by the typical water flow log analysis must be corrected to reflect the well bore water velocities in the absence of the logging tool. This may be accomplished using Equation 3:

where: vwc = corrected in situ water velocity
vw/ = log-reported water velocity
Apt - hydraulic area around tool
Ap = cross sectional area of tubular or casing
[0061] The cross sectional area for flow that is used in most production systems analysis computations are based on the hydraulic or effective diameter concept for fluid flow in an annulus. The area for tubular flow in casing below the end of a tubing string, or in a tubing string (Ap) is shown in Equation 4, where the area is in units of square feet and the inside diameter of the tubular is
in units of inches:

[0062] The hydraulic area to flow in the annulus between the casing or tubing and the pulsed neutron logging tool for which the log analysis reported water velocities are given is presented in Equation 5, where d is equal to {D-d)9 and d represents the outside diameter of the pulsed neutron logging tool and D represents the inside diameter of the casing or larger tubing, where the area is in units of square feet and the inside diameter of the tubular is in units of inches:



PATENT

[0067] When the equivalent pipe diameter relationship in Equation 9 is used for annular flow computations using tubular flow relationships, the equivalent diameter (de) maybe used to compute a fictitious velocity (v*) that may be used in production systems analysis pressure traverse computations. One example of a relationship that defines a fictitious velocity is provided in Equation 10, where the average velocity is in feet per second, the in situ flow rate is in cubic feet per second, and the equivalent diameter is in inches:

[0068] In order to get a more complete analysis of the multiphase flow system, it may be necessary to evaluate the interrelationships between liquid holdup, slip velocity, and the individual in situ fluid phase flow rates. The slip velocity is the difference between the average gas velocity and the average liquid velocity. The pressure traverse and multiphase flow correlations of Hagedorn and Brown, Duns and Ros, Orkiszewski, Beggs and Brill, Cullender and Smith, or various other multiphase flow pressure traverse correlations of interest maybe used to analyze the multiphase flow system. For example, the Duns and Ros technique established a fundamental relationship for multiphase flow between the slip velocity, the in situ gas velocity, the in situ liquid velocity, and the liquid holdup. That relationship is shown in Equation lla:


PATENT

where: vc = slip velocity
vsg = in situ superficial gas velocity
vsi = in situ superficial liquid velocity
vg = in situ average gas velocity
vi = in situ average liquid velocity
Yw = liquid holdup
[0069] As shown in Equations l1b and lie, the in situ average liquid velocity may be calculated by dividing the in situ superficial liquid velocity by the liquid holdup. Likewise, the in situ average gas phase velocity may be calculated by dividing the in situ superficial gas phase velocity by one minus the liquid holdup. In this disclosure, the term "/« situ gas phase velocity" is used generically to mean either the in situ average gas phase velocity or the in situ superficial gas phase velocity because one may be calculated if the other is known. Similarly, the term "m situ liquid phase velocity" is used generically to mean either the in situ average liquid phase velocity or the in situ superficial liquid phase velocity.
[0070] The above discussed examples include determining the liquid holdup (Yw) and the superficial velocity of water in the fluid flow stream (vs/). Even with these data, Equation lla still includes two unknowns, the in situ superficial gas velocity and the slip velocity. Both of these quantities are linked to the dimensional oil, gas, and water flow rates in the multiphase flow regime. A similar slip velocity relationship is also used between the oil and water phases in a three phase (oil, gas, and water) system analysis. In both slip velocity relationships (gas-liquid and oil-water), the slip velocity denoted by the particular relationship is the difference between the average velocity of the less dense phase and the more dense fluid phase. There are at least two options for solving for these variables.
[0071] One option that may be used is the implementation of a root-solving solution procedure to solve for the unknown in situ gas velocity that would satisfy the multiphase fluid flow relationships and result in a liquid holdup value that is comparable to that determined from the inelastic count rate

PATENT
of the Pulsed Neutron log measurements, for an assumed slip velocity. This technique has been found to work reasonably well for only the Duns and Ros multiphase flow correlation in practice since there are specific liquid holdup and slip velocity relationships for each flow regime. [0072] Another technique that may be used is an inverse interpolation. A series of forward solutions are computed with the previously determined multiphase pressure traverse, liquid holdup, and slip velocity correlations for a range of assumed gas flow rates and previously determined liquid flow rates and liquid holdup. A comparison is made not only of the computed liquid holdup values obtained, but also of the identification of the corresponding flow regime and slip velocity. The implementation of this type of in situ gas velocity and slip velocity evaluation procedure may be complex, but does result in a stable, reliable, and accurate solution procedure for simultaneously evaluating the in situ gas velocity and corresponding slip velocity.
[0073] A pulsed neutron log analysis, as described above, coupled with measured well pressures and temperatures, may be used to produce a reasonably direct production log evaluation of the multiphase flow in the well and inflow contributions of the individual completed intervals in the well. Most modern downhole tool assemblies include temperature and pressure sensors, enabling an independent verification of the conventional production log response. The verification may also be used to aid in the interpretation of conventional production log analyses where there are complexities, such as converging flow at the bottom of the tubing string or slug flow conditions. [0074] FIG. 8 shows a graph of data collected for a field well that shows the applicability and accuracy of pulsed neutron log inflow analysis. In this well, the end of the tubing is at a point above all of the completed intervals, thus permitting the use of a conventional production log as well as the pulsed neutron logging technique. A conventional production log with spinner, pressure and temperature probes was run, as well as an optical scanner to record the bubble counts in the flow stream. Multiple passes were made with the conventional production logging system, both with and against the direction of flow (uphill in this case).
[0075] Only a short distance existed between the end of the tubing string and the top of the shallowest completed interval in the well, and some convergence effects are seen in the conventional production log response of the well, increasing the difficulty in its interpretation, as well as some slugging effects in the production log response in the well. These complicating effects in the conventional production log analysis also help to demonstrate the advantage of using the pulsed



PATENT
[0078] The left-hand side of the graph in FIG. 8 shows the liquid flow rate, in stock tank barrels per day, as obtained with the pulsed neutron log analysis 801 and by conventional production log analysis 802. As shown in the graph, the two lines 801, 802 are in very close agreement. The right-hand side of the graph in FIG. 8 shows the gas flow rate, in thousands of standard cubic feet per day, as obtained with the pulsed neutron log analysis 803 and by the means of a conventional production log analysis 804. As shown in the graph, the two lines 803, 804 are in very close agreement. [0079] In fact, the pulsed neutron production log analysis results in better agreement with the measured gas production rate at the surface (1.7 MMscf/D) during the logging operations than was obtained with the conventional production log analysis (1.8 MMscf/D). Besides the problem previously stated regarding the flow convergence at the bottom of the tubing string, there was also a significant amount of slug flow evidenced in the production log behavior. Both the conventional measurement and pulsed neutron log analyses identified the flow regime over this interval of the well as being in the slug flow regime.
[0080] FIG. 9 shows one example of a method for determining an inflow performance for completed intervals in a well. The method first includes pulsing neutrons into a fluid flow in the wellbore, at step 901. In at least one example, the well fluid is production fluid flowing downwardly in an annulus. The well fluid flow may be any type of well fluid flow, including the examples described above. Next, the method may include detecting gamma rays from the decay of nitrogen-16 in the well fluid, at step 902. The neutron pulse will activate oxygen in the well fluid to form an amount of nitrogen-16. The decay of nitrogen-16 produces a gamma ray that may be detected by gamma ray detectors.
[0081] The method may include determining the velocity of the water in the well fluid. This is shown generally at steps 903 and 904. This may be performed in a single step by processor. In another example, determining the velocity of water in the fluid include measuring the time that it takes for the activated oxygen to flow from the pulsed neutron generator to the gamma ray detector, at step 903. In one example, the time is determined from a local maximum in the count rate at the detector. Determining the fluid velocity may then be accomplished by dividing the distance between the pulsed neutron generator and the gamma ray detector, at step 904.
[0082] Finally, the method may include estimating the inflow performance of one or more completed intervals in the well from data that includes the in situ water velocity, at step 905. This

PATENT
step may include using one of the various multiphase flow correlations known in the art, along with other data collected by sensors in the wellbore. This may include evaluating other flow properties, such as the in situ gas velocity, the slip velocity, the wellbore pressure and fluid mixture densities. These properties maybe used in a comprehensive production systems analysis computational model to determine the quantities at a plurality of locations in the wellbore, tubing string, and annulus between the casing and the tubing string. The wellbore fluid velocities, fluid mixture densities, and pressures at a plurality of locations in the wellbore, annulus, and tubing string are used to determine the inflow rates of fluid from or into each of the completed intervals for production zones encountered by the wellbore.
[0083] It is noted that the method may be used using more than one gamma ray detector. The data from the detectors may be combined in any manner known in the art to determine the velocity of the water. For example, the data may be averaged.
[0084] It is also noted that "inflow" is used to mean flow in any direction from a completed interval. In cases of injection or for a thief zone, the inflow may be flow from the well into the completed interval. The term inflow does not exclude such situations.
[0085] FIG. 10 shows a method for profiling multiphase fluid flow in a well. The method may include pulsing a first pulse of neutrons into the fluid flow, at step 1001. Next, the method may include detecting gamma rays produced from inelastic collisions between the neutrons and nuclei in the production flow, at step 1002. The gamma rays from inelastic collisions may be detected within the first few milliseconds following the first neutron pulse.
[0086] The method may then include determining the inelastic count rate of the fluid from the detected gamma rays, at step 1003. The flux and energy of the gamma rays is related to the inelastic count rate. The method may then include determining the liquid holdup of the fluid, at step 1004. In one example, estimating the liquid holdup is performed using Equation 1.
[0087] The method may include pulsing a second pulse of neutrons into the well fluid flow, at step 1005. In at least one example, the well fluid is production fluid flowing downwardly in an annulus. The well fluid flow may be any type of well fluid flow, including the examples described above. Next, the method may include detecting gamma rays from the decay of nitrogen-16 in the well fluid, at step 1006. The neutron pulse will activate oxygen in the well fluid to form an amount

PATENT
of nitrogen-16. The decay of nitrogen-16 produces a gamma ray that may be detected by gamma ray detectors.
[0088] The method may include estimating the velocity of the water in the well fluid. This is shown generally at step 1007. Specific methods for estimating the water velocity were described above with respect to steps 903 and 904 in FIG. 9. The method may include estimating the in situ gas velocity and the slip velocity, at step 1008. These two quantities may be separately determined, or they may be computed simultaneously using a multiphase fluid flow model. For example, the Duns and Ros correlations may be solved using an inverse interpolation or a root-solving method to determine the in situ gas velocity and the slip velocity.
[0089] It is noted that, while FIG. 10 includes a first neutron pulse and a second neutron pulse, the order of the pulses is not limited. For example, the second pulse and corresponding velocity measurements may be performed before the first pulse and the corresponding inelastic count rate determination. In addition, FIG. 10 shows a method for one particular location in a well. The method may be performed at a plurality of locations in the well so that the data may be used to resolve the inflow performance of each of a plurality of production zones in a reservoir. [0090] FIG. 11 shows one example of a method for determining an inflow performance for completed intervals in a well. The method first includes pulsing neutrons into a fluid flow in the wellbore, at step 1101. In at least one example, the well fluid is production fluid flowing downwardly in an annulus. The fluid flow may be any type of fluid flow, including the examples described above. Next, the method may include detecting gamma rays resulting from inelastic collisions between the neutrons and nuclei in the fluid flow, at step 1102.
[0091] The method may then include determining the inelastic count rate of the fluid from the detected gamma rays, at step 1103. The flux and energy of the gamma rays is related to the inelastic count rate. The method may then include estimating the liquid holdup of the fluid, at step 1104. [0092] Finally, the method may include estimating the inflow performance of one or more competed intervals in the well from data that includes the liquid holdup, at step 1105. This may include evaluating other flow properties, such as the in situ gas velocity, the slip velocity, the wellbore pressure and fluid mixture densities. These properties may be used in a comprehensive production systems analysis computational model to determine the quantities at a plurality of locations in the wellbore, tubing string, and annulus between the casing and the tubing string.

PATENT
[0093] There are a large number of well flow configurations for which a pulsed neutron log may be used for multiphase flow analysis. Among these configurations are the situations where: (1) there is production from the tubing, but the annulus is static at the surface; (2) there is production from the annulus, but the tubing is static at the surface; (3) there is production from both the tubing and annulus at the surface; (4) there is injection into the tubing and production from the annulus at the surface; (5) there is injection into the annulus and production from the tubing at the surface; (6) there is injection into the tubing, but the annulus is static at the surface; (7) there is injection into the annulus and the tubing is static at the surface; and (8) there is injection into both the tubing and annulus at the surface. In addition, pulsed neutron analysis maybe useful in a flow situation where there is production from several completed intervals in a multilayer reservoir, but where at least one interval lacks the pressure to produce. The production fluids from other intervals may flow into such an interval. A "thief zone," such as this, maybe identified using neutron log analysis. Neutron log analysis may present advantages in each of these different flow scenarios.
[0094] It is noted that injection relates to pumping fluids and materials from the surface into the well, and ultimately into the formation. As is known in the art, this is done for several purposes, for example, waste disposal or to charge a well. In such a case, the flow will be from the well and into the formation. Nonetheless, the term inflow is used generically to describe these situations, even though the fluid is flowing out of the well.
[0095] There are also multiple well tubular and completion scenarios for which neutron log analysis may be applicable. These include the cases of: (1) there is no tubing in the well and flow is only in the casing; (2) there is tubing in the well that is set above all of the completed intervals; (3) there is tubing in the well that is set below the top of the shallowest completed interval and above the bottom of the deepest completed interval; and (4) the tubing is set below all of the completed intervals in the well. The specification of a sealing packer can also be made in the analysis in which the packer can be above all of the completed intervals in the well or below one or more of the completed intervals. A pulsed neutron log production inflow analysis, as described, may be used with other completion scenarios as well.
[0096] The specification of a mixed direction (uphill and downhill flow) for the measured water velocity is permitted. The specification of the start and stop depths of the timing of the water (liquid) velocity measurements directly permits the proper assignment of the appropriate signs for the water

PATENT
velocity values specified in the analysis. The specified measured water velocity table is merged in the analysis with the specified completed intervals in the well that have been assigned to properly provide a realistic water velocity profile for interpolation in the analysis. This is particularly important in cases where there is a sparse set of water velocity measurements made in the well. [0097] It is also possible to take water flow velocity measurements above and below each of the completed intervals in the well, except the deepest completed interval in the well, which may only require the water velocity measurement above the completed interval. The surface flow rates of each of the fluid phases may also be measured at the surface, in both the tubing and annulus (if applicable). In addition, the temperatures and pressures in the tubing and annulus pressures and temperatures may also be measured during the pulsed neutron log measurements. [0098] In some cases, the velocity of the gas maybe insufficient to entrain the liquid particles and continuously remove them from the well. In such a case, a portion of the liquid may fall back down the well and accumulate at the bottom of the well. When this occurs, the measurement of the liquid holdup may be erroneously high due to the presence of additional liquid in the well, creating an error in the estimation of the in situ gas phase velocity and the slip velocity. An appropriate correction may be required to determine the true in situ gas phase velocity. [0099] The "critical velocity" of the gas as used in the context of this discussion is the minimum average gas velocity where the upward drag force on the largest droplets of entrained liquid is equal to the downward gravitational forces on the droplet. If the average gas velocity is super-critical, the entrained liquid will flow upwardly with the gas. If the average gas velocity is sub-critical, the liquid droplets will fall back down into the well.
[00100] The critical velocity for continuous transport of entrained liquid droplets (yg
where a = Interfacial (surface) tension between the liquid and gas;
Pi = Liquid density;
pg = Gas density; and
Cd = Coefficient of drag.



PATENT
Eq. 14
where a = Interfacial (surface) tension between the liquid and gas; fig = Gas viscosity; and vg = In situ average gas phase velocity.
[00104] Equations 12-14 include three equations with three unknowns (vg-crit Cd, and NRE), thus, a computational solution must be used to solve all three equations simultaneously, using Equation 12 as the basis function. The values of the unknowns may be evaluated until they converge to a solution. Once a solution is obtained, the critical gas phase velocity for continuous transport of entrained liquid droplets in Equation 12 may be compared with the in situ average gas phase velocity (vg). If the in situ average gas phase velocity (vg) is greater than the critical gas phase velocity (vg, ait), Ae flow of gas is sufficient to unload the liquids from the well, and no correction is required. If, on the other hand, the in situ average gas phase velocity (vg) is less than the critical gas phase velocity (vg.cri{) for continuous liquid transport from the well, then the calculation of the in situ gas phase velocity may have been affected by fall back of liquid in the well because of the possibility that the in situ average gas phase velocity (vg) is insufficient to unload liquids from the well. [00105] Before making a correction to either the in situ average gas phase velocity (vgg) or the in situ superficial gas phase velocity (vsg), an evaluation of the flow type must be made. In the bubble flow regime, the gas phase is distributed as bubbles throughout the liquid phase. In the slug flow regime, a series of liquid slugs are separated by relatively large pockets of gas. The slugs occupy most of the cross-sectional area of the tubing. As will be explained, a correction for the bubble flow regime may use the relationships in Equations 15-17, below, and a correction for the slug flow regime may use the relationships in Equations 18-21, below. Other flow regimes are not treated herein because other flow regimes will generally be above the critical gas phase velocity. [00106] In the bubble flow regime, the bubble rise velocity (v) may be computed using Eq. 15:
Eq. 15





PATENT
where g = acceleration of gravity;
dtt = Outside diameter of the tubing;
dc = Inside diameter of the casing;
pt = Liquid density; and
pg = Gas density.
[00115] Using the Taylor bubble rise velocity (v^) for annular slug flow, as determined in Equation 20, the liquid holdup (Yww) may be estimated by the relationship in Equation 21:

where vsg = in situ superficial gas phase velocity;
vssi = in situ superficial liquid phase velocity;
dt = outside diameter of the tubing; and
dc = inside diameter of the casing.
[00116] The right side of Equation 21 may then be equated to the right side of Equation 1, and the corrected in situ superficial gas phase superficial velocity (vsgg) may be determined by solving for vsg, using the liquid holdup values derived from Equation 1, the Taylor bubble rise velocity (vrt) computed using Equation 20, and the other known variables in Equation 21. [00117] FIG. 13 shows one example of a method for correcting an in situ gas phase velocity measurement. The method includes estimating the critical gas phase velocity, at step 1301. This is the minimum gas velocity at which the entrained liquid will be removed from the well. In one example, the critical gas phase velocity (vg-crit) may be estimated using Equation 12, above. The method next includes comparing the in situ average gas phase velocity with the estimated critical gas phase velocity, at step 1302. In some examples, the in situ gas phase velocity is estimated using one or more of the methods described above and show, for example, in Figures 9-11. [00118] If the in situ average gas phase velocity is greater than the estimated critical gas phase velocity, then the in situ average gas phase velocity is great enough to continuously unload the liquid from the well, and no correction for liquid fallback is required. Thus, if the in situ average gas phase velocity is greater than the estimated critical gas phase velocity, the decision may include ending the method, at step 1303.

PATENT
00119] If, however, the in situ average gas phase velocity is less than the critical gas phase velocity, then the in situ average gas phase velocity is not sufficient to continuously transport the entrained liquid particles and unload the well. The liquid particles may "fall back" into the well and collect at the bottom. In this case, the measurements and calculations for determining the inflow profile of the well, including the in situ gas phase velocity, may require correction for the liquid fall back.
[00120] In the case where the in situ average gas phase velocity is less than the estimated critical gas phase velocity, the method may next include determining the flow regime, at step 1304. In one example, the flow regime may be selected as either the bubble flow regime or the slug flow regime. [00121] In the case where the flow regime is the bubble flow regime, the method may next include estimating the bubble rise velocity, at step 1305. In one example, the bubble rise velocity is estimated using Equation 15, above. The method may also include computing a corrected in situ gas phase velocity, at step 1306. In one example for tubular flow, the corrected in situ gas phase velocity is determined by setting the right side of Equation 1, above, equal to the right side of Equation 16, above. In another example for annular flow, the corrected in situ gas phase velocity is determined by setting the right side of Equation 1, above, equal to the right side of Equation 17, above. [00122] In the case where the flow is in the slug flow regime, the method may include estimating the Taylor bubble rise velocity, at step 1307, and then computing the corrected in situ gas phase velocity, as step 1308. In one example for tubular flow, the Taylor bubble rise velocity is computed using Equation 18, above. In a further example, the corrected in situ gas phase velocity is determined by setting the right side of Equation 1, above, equal to the right side of Equation 19, above. In another example for annular flow, the Taylor bubble rise velocity is computed using Equation 20, above. In a further example, the corrected in situ gas phase velocity is determined by setting the right side of Equation 1, above, equal to the right side of Equation 21, above. [00123] Although this detailed description has shown and described illustrative embodiments of the invention, this description contemplates a wide range of modifications, changes, and substitutions. In some instances, some features of the present invention may be employed without a corresponding use of other features. Accordingly, it is appropriate that readers should construe the appended claims broadly, and in a manner consistent with the scope of the invention.

Documents:

688-che-2006 abstract.pdf

688-che-2006 claims.pdf

688-che-2006 correspondence others.pdf

688-che-2006 description(complete).pdf

688-che-2006 form-1.pdf

688-CHE-2006 FORM-13 09-02-2009.pdf

688-che-2006 form-3 09-02-2009.pdf

688-che-2006 other patent document 09-02-2009.pdf

688-CHE-2006 OTHER PATENT DOCUMENT 04-09-2009.pdf

688-che-2006-abstract.pdf

688-che-2006-claims.pdf

688-che-2006-correspondence-others.pdf

688-che-2006-description-complete.pdf

688-che-2006-drawings.pdf

688-che-2006-form 1.pdf

688-che-2006-form 3.pdf

688-che-2006-form 5.pdf


Patent Number 245229
Indian Patent Application Number 688/CHE/2006
PG Journal Number 2/2011
Publication Date 14-Jan-2011
Grant Date 10-Jan-2011
Date of Filing 13-Apr-2006
Name of Patentee PRAD RESEARCH AND DEVELOPMENT LIMITED
Applicant Address PO BOX 71, CRAIGMUIR CHAMBERS, ROAD TOWN, TORTOLA, BRITISH VIRGIN ISLANDS
Inventors:
# Inventor's Name Inventor's Address
1 POE, BOBBY, D., 13502 WESTPORT LANE, HOUSTON, TEXAS 77079
2 BUTSCH, ROBERT, J., 15406 GETTYSBURG, TOMBALL, TEXAS 77377
PCT International Classification Number G01V5/00
PCT International Application Number N/A
PCT International Filing date
PCT Conventions:
# PCT Application Number Date of Convention Priority Country
1 60/671,944 2005-04-15 U.S.A.
2 11/278, 444 2006-04-03 U.S.A.