The manner in which the objectives of the invention and other desirable characteristics may be obtained is explained in the following description and attached drawings in which:
It is to be noted, however, that the appended drawings are not to scale and illustrate only typical embodiments of this invention, and are therefore not to be considered limiting of its scope, for the invention may admit to other equally effective embodiments.
In the following description, numerous details are set forth to provide an understanding of the present invention. However, it may be understood by those skilled in the art that the present invention may be practiced without these details and that numerous variations or modifications from the described embodiments may be possible.
All phrases, derivations, collocations and multiword expressions used herein, in particular in the claims that follow, are expressly not limited to nouns and verbs. It is apparent that meanings are not just expressed by nouns and verbs or single words. Languages use a variety of ways to express content. The existence of inventive concepts and the ways in which these are expressed varies in language-cultures. For example, many lexicalized compounds in Germanic languages are often expressed as adjective-noun combinations, noun-preposition-noun combinations or derivations in Romanic languages. The possibility to include phrases, derivations and collocations in the claims is essential for high-quality patents, making it possible to reduce expressions to their conceptual content, and all possible conceptual combinations of words that are compatible with such content (either within a language or across languages) are intended to be included in the used phrases.
The invention describes methods for obtaining a wellbore schematic, defined as the relationship between true vertical distance (TVD) and measured distance (MD) of a tubular in a wellbore. Currently, in wellbore cleanout procedures and other procedures where liquids are pumped into the wellbore via tubing and out through the annulus, if hydrostatic head pressure may be removed, one has an accurate estimate of the wellbore pressure at the bottom of (entrance to) the annulus. However, the only way to remove the hydrostatic component from downhole data is to have a copy of the wellbore schematic in advance of the job. This schematic could have been obtained while drilling the well via measurement-while-drilling data, or after drilling by lowering a wireline inclinometer tool such as a gyroscope. However, no tool that is currently used for stimulating reservoirs is known to have an internal inclinometry platform, nor is there known any previously existing method to determine TVD strictly from pressure data and flow rate information. Another challenge is in so-called wellbore cleanout operations, wherein various fill materials are carried by a fluid injected down the wellbore, typically through coiled tubing or other tubulars, and flowed out through the annulus. The cleanout fluid carrying solid particles along the annulus is a suspension whose density correlates with the concentration of solid particles. For an effective cleanout the suspended particles must be transported all the way out of the well. The hydrodynamic pressure in the annulus is directly proportional to the suspension density. It would be an advance in the art if methods could be devised that provide information about the relationship between TVD vs. MD, in other words the wellbore schematic, while flowing fluids into the wellbore. It would further be an advance in the art to use the obtained wellbore schematic to monitor and/or control wellbore operations, such as wellbore cleanout procedures, via information about the annulus. There is a continuing need for systems and methods that address one or more of these challenges.
As used herein “wellbore schematic” means the relationship between true vertical depth and measured depth, where measured depth is the depth measured at the wellhead of coiled tubing that has entered the wellbore. As used herein “annulus fluid” and “annular fluid” may be used interchangeably and refer to the fluid traversing past a coiled tubing back to the surface. As used herein “wellbore servicing” means any operation designed to increase hydrocarbon recovery from a reservoir, reduce non-hydrocarbon recovery (when non-hydrocarbons are present), or combinations thereof, involving the step of pumping a fluid into a wellbore, or into coiled tubing that is or will be placed into the wellbore. This includes pumping fluid into a reeled or spooled coil of coiled tubing. The fluid pumped may be a composition to increase the production of a hydrocarbon-bearing zone, a composition pumped into other zones to block their permeability or porosity, a composition designed to flush or cleanout a wellbore or portion thereof, and the like. Methods of the invention may include pumping fluids to stabilize sections of the wellbore to stop sand production, for example, or pumping a cementatious fluid down a wellbore, in which case the fluid being pumped may penetrate into the completion (e.g. down the innermost tubular and then up the exterior of the tubular in the annulus between that tubular and the rock) and provide mechanical integrity to the wellbore. As used here in the phrases “treatment” and “servicing” are thus broader than “stimulation”. In many applications, when the rock is largely composed of carbonates, one of the fluids may include an acid and the hydrocarbon increase comes from directly increasing the porosity and permeability of the rock matrix. In other applications, often sandstones, the stages may include proppant or additional materials added to the fluid, so that the pressure of the fluid fractures the rock hydraulically and the proppant is carried behind so as to keep the fractures from resealing. The details are covered in most standard well service texts and are known to those skilled in the well service art so are omitted here.
Methods within this aspect of the invention include sending real-time pressure data to the surface during wellbore servicing using one or methods selected from wireless methods (such as mud-pulse electromagnetic telemetry), wire methods via a data-carrying wire (such as an eline cable), and fiber-optic lines. The wireless methods may be used particularly when running in joints of tubing. In other embodiments the tubing is brought to the well spooled onto a reel with a telemetry cable already inserted into the spool, but the invention is not so limited. The wireline may be inserted into the tubing at the well site. An advantage of fiber-optic telemetry is that the bottom-hole pressure may be measured without the need for downhole electronics. Indeed, if one has downhole electronics, then an inclinometer may be added to the electronic package for minimal additional cost, so one of the prime advantages of this invention is for bottom-hole assemblies without an electronics package. Fiber-optic techniques to measure pressure are well-known in the industry. One common device relies on interferometry to identify the size of a cavity, that cavity itself changing size based on the external pressure applied to the cavity. Such devices are made, for example, by FISO Technologies in Montreal, Canada and have been implemented in the bottom-hole assemblies.
Exemplary methods of the invention rely on running tubing into the bottom of a wellbore while pumping a fluid therethrough at varying rates while running in. The fluid may be one in which the friction drop down the coiled tubing behaves according to a power-law relationship:
friction pressure=A*(flow rate)n,
where n is an exponent (typically between 1 and 2) and A depends on viscosity of the fluid, local friction effects and tubular internal diameter. The pressure measured at the bottom of the tubing will be given by the circulating pressure (measured at the surface) less the friction pressure through the tubing plus the hydrostatic pressure. The friction pressure in the coiled tubing may be best modeled as two components:
friction pressure=A1*(flow rate)n1+A2*(flow rate)n2
where the first term to the right of the equal sign represents the pressure drop along that part of the coil wound around a spool, and the second term to the right to the equal sign represents the pressure drop along the unspooled coil. This latter component may be taken to be proportional to the length of coil run into the wellbore, so surface measurement of this length will be needed. Apparatus for such measurements are commercially available and well-known in the industry. For example, small wheels may be pushed against the coil and the rotation of those wheels will give the length of the coil run in. One embodiment is that known under the trade designation UTLM, from Schlumberger. The first component of the friction pressure may be modeled either as a formula which takes into account the changing diameter of the spooled coil, or more simply may be taken as proportional to the length of coil wound around the spool. Thus if there is a total of LT feet brought to the rig and MD(t) has been run into the ground at time t, then the friction pressure may take the form:
friction pressure=a1*(LT−MD(t))*(flow rate(t))n1+a2*MD(t)*(flow rate(t))n2.
In order to determine the unknown coefficients a1 and a2, and the exponents n1 and n2, the flow rate and MD as the coil is run in may be varied with time. The hydrostatic pressure will be proportional to the density of the fluid times its TVD. In many embodiments the density of the pumped fluid varies with depth and flow rate; however, in some embodiments the density may be assumed to be fixed, so the hydrostatic term becomes:
hydrostatic pressure=TVD(t)*density*gravity.
One method of the invention is thus to find a best fit of the parameters which matches up the sum of the theoretical friction pressure and hydrostatic pressure against the difference of the measured circulating and bottomhole pressures. This best fit may be done with a number of techniques for non-linear optimization. Such programs are readily available in software packages, such as Matlab. The result is then a cross-plot of TVD(t) versus MD(t) at each time. This is precisely the wellbore schematic. The terminology Y(t) may be used to denote the difference between theoretical pressure drop in the coil against the measured pressure drop.
In the unlikely event that the density of the fluid is not known at the beginning of the job, it may be estimated if the wellbore schematic is at least known at the top of the wellbore, e.g., if the top of the wellbore is vertical. This estimate could then be used for the rest of the inversion.
The nature of nonlinear parameter estimation means that the plot of TVD(t) versus MD(t) will be quite noisy. This estimation may be made more robust by adding additional information such as the maximum dogleg angle of the wellbore. A second piece of information is that the borehole inclination may only change quite slowly with depth. A standard practice in the industry is to assume that the borehole schematic follows a so-called minimum radius of curvature. While drilling the well, periodic measurements of inclination are passed to the surface. The inclination between two such measurements is determined by fitting an arc of a circle of fixed radius such that the inclinations at the ends of the arc match the measured inclinations. In effect, the wellbore schematic is that combination of arcs that has a fixed radius between each measurement of inclination. We may use this methodology in the derivation of the wellbore schematic from pressure.
The unknown parameters become a1, a2, n1 and n2 and a series of inclinations, θ(MD), where θ is the inclination angle and MD is the length of coil run into the well. The nonlinear estimation will then minimize the sum of y(t)2+Z(t)2 where Z(t) is a weighting term constraining the rate of change of θ. There are well-known techniques to constrain rate of change. One standard formula is the sum of the absolute value:
rate of change=|θ(MD(j+1))−θ(MD(j))|
for a predetermined selection of depths MD(1), MD(2), . . . . A typical selection of depths would be fixed interval of 10 m or 30 ft along the length of the wellbore. The result of this optimization is not just the wellbore schematic. The parametric values in the friction expression are in themselves useful because they may give indications of viscosity and the nature of the flow—for example, the exponent of the flow is indicative of the flow profile, whether it is laminar or turbulent. See for example, Bird, et al., “Transport Phenomena”, Chapter 6, pp. 180-190, John Wiley & Sons (1960).
Once the wellbore schematic is estimated then the same information on the wellbore schematic and pressure of fluids may be used to analyze the annulus fluid around the coil. The pressure drop between the bottomhole and the wellhead is the sum of the hydrostatic and friction pressures in the annulus, plus the effect of the reservoir (e.g. whether it is causing a net increase in pressure in the annulus or a decrease). Also the hydrostatic pressure at a given depth may be subtracted from the annular bottomhole pressure to get directly the effect of the formation pressure (and the changes in that formation pressure vs. time). For example, if the tool is stationary then the hydrostatic pressure may be subtracted from pressure measurements during a fall-off and formation parameters may be estimated using standard well-testing techniques. If the tool is not stationary, then to be able to use such techniques requires subtracting of the varying hydrostatic pressure versus depth. Interestingly, if there is a small error in the input fluid density then there will be a corresponding error in estimated TVD, but this would not then translate into an error in the estimated hydrostatic versus depth.
It is important that the flow-rate be varied during the run in the well. If a fixed flow-rate is used then deriving the parameters a1, a2, n1 and n2 will be very unstable.
Note that there is a significant advantage in transmitting the bottom-hole pressure in real-time because then the wellbore schematic may be determined without having to extract the coiled tubing.
Further, note that in a typical coiled tubing operation, there will be repeated passes through the wellbore, so that during the course of the operation, the uncertainty in the wellbore schematic will be removed. The surface operator (or his computer) will need to monitor which fluids are being pumped, which in turn would allow parameters a1, a2, n1, n2 and density to vary from one fluid to the next.
Referring now to the drawing figures,
Measurements (see
Wellhead pressure: WHP, measured at surface at the flow exit.
Circulation pressure: Pcirc, measured at surface, inside the CT at the ‘in’ extremity.
Annulus bottom hole pressure: Pan, measured in the wellbore, at the end of the CT string.
CT bottom bole pressure: PCT, measured inside the CT, at the bottom end.
Parameters:
Total CT string length: LT
CT length in hole: MD
Wellbore, CT radii (resp. diameters): rw, rCT (resp. dw, dCT).
Friction coefficient: f
annulus fluid velocity: υan.
The four measured pressures are linked by the following relationships:
P
an
=WHP+F
an
+H
an (1)
P
CT
=P
circ
−F
CT
+H
CT (2)
P
CT
=P
an
+DP
nozzle (3)
Notations: F for friction pressure, H for hydrostatic pressure. The subscripts ‘an’ and ‘CT’ stand respectively for ‘in the annulus or wellbore’ and ‘inside the coiled tubing’.
With the annulus friction pressure—in theory calculable—and the measured annulus and wellhead pressures, the quantity of interest, the annulus hydrostatic pressure, is inferred from (1):
H
an
=P
an
−WHP−F
an (1a)
The hydrostatic pressure is also:
H
an
=ρ
an
·g·TVD (4)
wherein TVD is the vertical depth, equal to MD as defined above in a vertical well. We therefore obtain the average annulus fluid density ρan:
Obtaining the annulus friction pressure:
The friction in the wellbore is given by (with the usual assumptions):
Note: The density may vary along the annulus, i.e., ρ=ρ(MD). ρan in (5) is the average annulus fluid density given by:
and the friction pressure is:
The annulus friction pressure is a function of the annulus fluid density, i.e., the friction term in (1a) cannot be accessed without knowing the density. An estimate of the density can be used in (5) to get the friction loss, and then re-adjusted at each computation cycle after the set of equations (1a, 4a) has been solved. The friction term could be very inaccurate, one of the reasons being that it requires the friction coefficient f, which has large uncertainties.
The following scheme gives the variations of the annulus fluid density without having to calculate the friction pressure.
Case 1: Vertical well.
Equation (5) may be re-written:
F
an
=MD ρ
an
fk
geoυ2an (8)
where kgeo is a constant that depends on the geometry of the system. Note that equation (8) is not specific to the vertical case.
From equations (1, 4, 8) one obtains equations (9) and (9a):
The difference between the downhole annulus pressure and the wellhead pressure is proportional to the hydrodynamic pressure and density for any given flow rate. It follows that:
Even without knowing the friction in the annulus, the measured quantity (Pan−WHP)/MD gives the variations of the density in the annulus.
With f*kgeo known the method is quantitative (both f and kgeo are accessible, an experimental method for estimating the product f*kgeo is described further).
Case 2: Deviated well.
In a deviated well we lose the proportionality between Han and MD. Assuming a constant deviation, if m is the cosine (deviation angle), and reviewing
H
an
ρ
an
·g·[MD
0
+m·(MD−MD0)] (10)
and from equations (1, 8, and 10), equations 11 and 11a may be obtained:
Equation (11) may be solved for ρan, given the well configuration. Another option is a chart of (Pan−WHP) vs. MD with a set of pre-defined constant-density lines.
Friction test, an experimental method for estimating the product f*kgeo:
While in hole, a friction test could be performed based on equations (9) or (11).
Before starting cleaning, i.e., no particles in suspension, equations (9 or 11) may be solved for the quantity f*kgeo which characterizes the friction. This must be done before reaching the treatment zone so as to have a density well defined (density of the injected fluid).
Interrupted Flow Test
While flowing, short pump interruptions (van=0) will allow solving equations (9, 11) for “ρan” without the need of taking into account friction. Such pumping interruptions may only be possible when the particles' settling times are long enough, i.e., with gel fluids and the like.
Derivation of well deviation and TVD (cf
The well has a vertical section of length MD0, the wellbore deviation is a function of the measured depth MD. The friction pressure is still given by (7):
The hydrostatic pressure is:
H
an=∫ρ(MD)·g·cos [θ(MD)]·dMD (12)
P
an
−WHP=∫ρ(MD)·gcos [θ(MD)]·dMD+∫ρ(MD)·f·kgeoνan 2·dMD (13)
The left hand side of (14) is measured. Equation (14) may be solved for ρ(MD) given the well trajectory (i.e. cos [θ(MD)] vs. MD) or for cos [θ(MD) given the density (i.e. ρ(MD) vs. MD). After equation (14) is solved the TVD may be obtained through:
TVD=∫cos [θ(MD)]·dMD (15)
a, b, c) illustrate an example of application of the method. In this example, given the diameter of the wellbore and fluid flow rate, the friction term is negligible compared to the hydrostatic term. The density, which is constant, is determined while in the vertical portion of the well.
In many cases, it is advantageous to drain a reservoir with a multiplicity of wellbore branches connected together downhole to a main trunk wellbore, in the similar way that the roots of a plant retrieve water from the soil. Such wellbores are referred to as multilaterals, with each branch being referred to as a lateral. In such circumstances, it is important to know which branch of the reservoir has been penetrated by the coiled tubing. Using one or more embodiments of the invention described herein, a wellbore schematic can be determined from parameters measured on the coiled tubing. The derived wellbore schematic can be compared to a schematic of the multilateral well, and thereby identify which of the laterals has been penetrated. Note that only an approximate schematic is needed of the overall multilateral reservoir. As the coiled tubing penetrates a particular lateral, then a more accurate description of the multilateral reservoir can be obtained. With the help of an entry sub at the end of the coiled tubing, it is possible to enter many, or all, the laterals and so obtain a complete multilateral schematic.
Knowing which lateral has been penetrated is also important to optimize the reservoir stimulation. For example, if a water is being produced out of one lateral and hydrocarbon out of a second, then the operator will desire to pump a stimulating fluid, such as acid, into the hydrocarbon-containing lateral, and the operator will desire to pump a non-stimulating or viscous fluid, such as a gel, into the water-containing lateral. If these fluids were to be pumped into the wrong laterals, then overall hydrocarbon recovery would be ruined. Similarly, if many laterals are penetrating hydrocarbon, then it will be efficient to add stimulating fluids to each lateral. If the coiled tubing should accidentally re-enter an already stimulated lateral, then it is disadvantageous to pump more stimulating fluid into that lateral. In this way, it can be seen that increasing knowledge of the wellbore schematic penetrated by the coiled tubing is a means to increase overall hydrocarbon productivity. The ability to selectively choose fluids is only one such example of how to use wellbore information to increase overall hydrocarbon productivity and other applications will be immediately apparent to those skilled in the art.
In certain embodiments of the invention communication from the communication line to a surface data acquisition system may comprise wireless telemetry. The surface data acquisition system need not be at the well site, for example it may be a networked system including a computer at the well site and a second system at some remote location. The data transmitted may optionally be used to control the operation, whereby the pump rate or the composition of a treatment fluid is adjusted based purely upon the downhole data collected and transmitted by the communication line, or from a combination of downhole data and surface measurements.
As used herein, “pumping” means using a “pumping system”, which in turn means a surface apparatus of pumps, which may include an electrical or hydraulic power unit, commonly known as a powerpack. In the case of a multiplicity of pumps, the pumps may be fluidly connected together in series or parallel, and the energy conveying the pumped fluid may come from one pump or a multiplicity. The pumping system may also include mixing devices to combine different fluids or blend solids into the fluid, and the invention contemplates using downhole and surface data to change the parameters of the fluid being pumped, as well as controlling on-the-fly mixing.
By the phrase “surface acquisition system” is meant one or more computers at the well site, but also allows for the possibility of a networked series of computers, and a networked series of surface sensors. The computers and sensors may exchange information via a wireless network. Some of the computers do not need to be at the well site but may be communicating via a communication system such as that known under the trade designation InterACT™ or equivalent communication system. In certain embodiments a communication line may terminate at the wellhead at a wireless transmitter, and the downhole data may be transmitted wirelessly. The surface acquisition system may have a mechanism to merge the downhole data with the surface data and then display them on a user's console. The surface acquisition system may also include apparatus allowing communication to the downhole sensors.
Data transmitted from the communication line may be used to monitor subsequent stages of reservoir or wellbore treatment. The data transmitted may optionally be used to control some or all of the treatment operation, whereby for example a pump rate or composition of a fluid being injected is adjusted based purely on the downhole data obtained by the communication line, or from a combination of downhole data and surface measurements. The downhole data transmitted may be that from one or more sensors attached to the end of one or more communication lines, and may supplement or be supplemented by a variety of other measurements. The data may be from a distributed section of a communication line such as distributed temperature along an optical fiber. The data collected may be stored on the acquisition system and the information used to optimize and/or model subsequent stimulation runs.
Although only a few exemplary embodiments of this invention have been described in detail above, those skilled in the art may readily appreciate that many modifications are possible in the exemplary embodiments without materially departing from the novel teachings and advantages of this invention. Accordingly, all such modifications are intended to be included within the scope of this invention as defined in the following claims. In the claims, no clauses are intended to be in the means-plus-function format allowed by 35 U.S.C. §112, paragraph 6 unless “means for” is explicitly recited together with an associated function. “Means for” clauses are intended to cover the structures described herein as performing the recited function and not only structural equivalents, but also equivalent structures.