You must log in to edit PetroWiki. Help with editing

Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. More information

# Fluid flow through permeable media

This article discusses the basic concepts of single-component or constant-composition, single phase fluid flow in homogeneous petroleum reservoirs, which include flow equations for unsteady-state, pseudosteady-state, and steady-state flow of fluids. Various flow geometries are treated, including radial, linear, and spherical flow.

## Ideal reservoir model

Virtually no important applications of fluid flow in permeable media involve single component, single phase 1D, radial or spherical flow in homogeneous systems (multiple phases are almost always involved, which also leads to multidimensional requirements). The applications given in this Chapter are based on a model that includes many simplifying assumptions about the well and reservoir, and are interesting mainly only from a historical perspective See "Reservoir Simulation" for proper treatment of multi-component, multiphase, multidimensional flow in heterogeneous porous media. The simplifying assumptions are introduced here as needed to combine the law of conservation of mass, Darcy’s law, and equations of state to obtain closed-form solutions for simple cases.

Consider radial flow toward a well in a circular reservoir. Combining the law of conservation of mass and Darcy’s law for the isothermal flow of fluids of small and constant compressibility yields the radial diffusivity equation, ^{[1]}

In the derivation of this equation, it is assumed that compressibility of the total system, *c*_{t}, is small and independent of pressure; permeability, *k* , is constant and isotropic; viscosity, *μ*, is independent of pressure; porosity, *ϕ*, is constant; and that certain terms in the basic differential equation (involving pressure gradients squared) are negligible. The grouping 0.0002637*k*/*ϕμc*_{t} is called the hydraulic diffusivity and is given the symbol *η*. None of these assumptions are valid in general.

## Line-source solution to the diffusivity equation

Assume that a well produces at constant reservoir rate, *qB*; the well has zero radius; the reservoir is at uniform pressure, *p*_{i}, before production begins; and the well drains an infinite area (i.e., that *p* → *p*_{i} as *r* → ∞). Under these conditions, the solution to **Eq. 1** is^{[1]}

where *p* is the pressure at distance *r* from the well at time *t*, and

the*Ei *function or exponential integral.

The *Ei*-function solution is an accurate approximation to more exact solutions to the diffusivity equation (solutions with finite wellbore radius and finite drainage radius) for 3.79 × 10^{5} *ϕμc*_{t}*r*_{w}^{2}/*k* < *t* < 948 *ϕμc*_{t}*r*_{e}^{2}/*k*. For smaller times, the assumption of zero well size (line source or sink) limits the accuracy of the equation; for larger times, the reservoir’s boundaries affect the pressure distribution in the reservoir, so that the reservoir is no longer infinite acting.

For the argument, *x*, of the*Ei *function less than 0.01, the*Ei *function can be approximated with negligible error by

For *x* > 10, the*Ei *function is zero for practical applications in flow through porous media. For 0.01 < *x* < 10,*Ei *functions are determined from tables or subroutines available in appropriate software. ^{[2]}

## Altered zone and skin factor

In practice, most wells have reduced permeability (damage) near the wellbore resulting from drilling or completion operations. Many other wells are stimulated by acidization or hydraulic fracturing. **Eq. 2** fails to model such wells properly. Its derivation includes the explicit assumption of uniform permeability throughout the drainage area of the well up to the wellbore. Hawkins^{[3]} pointed out that if the damaged or stimulated zone is considered equivalent to an altered zone of uniform permeability. *k*_{s}, and outer radius, *r*_{s}, the additional pressure drop, Δ*p*_{s}, across this zone can be modeled by the steady-state radial flow equation

**Eq. 5** states that the pressure drop in the altered zone is inversely proportional to *k*_{s} rather than to *k* and that a correction to the pressure drop in this region must be made. Combining **Eqs. 2** and **5**, we find that the total pressure drop at the wellbore is

For *r* = *r*_{w}, the argument of the*Ei *function is sufficiently small after a short time that we can use the logarithmic approximation; thus, the drawdown is

We can conveniently define a dimensionless skin factor, *s*, in terms of the properties of the equivalent altered zone:

Thus, the drawdown is

**Eq. 9** provides some insight into the physical significance of the algebraic sign of the skin factor. If a well is damaged (*k*_{s} < *k*), *s* will be positive, and the greater the contrast between *k*_{s} and *k* and the deeper into the formation the damage extends, the larger the numerical value of *s*, which has no upper limit. Some newly drilled wells will not flow before stimulation; for these wells, *k*_{s} = 0 and *s* → ∞. If a well is stimulated (*k*_{s} > *k*), *s* will be negative, and the deeper the stimulation, the greater the numerical value of *s*. Rarely does a stimulated well have a skin less than –7, and such skin factors arise only for wells with deeply penetrating, highly conductive hydraulic fractures. If a well is neither damaged nor stimulated (*k* = *k*_{s}), *s* = 0.

The altered zone near a well affects only the pressure near that well; that is, the pressure in the unaltered formation away from the well is not affected by the existence of the altered zone. Thus, use **Eq. 9** to calculate pressures at the sandface of a well with an altered zone, and **Eq. 2** to calculate pressures beyond the altered zone in the formation surrounding the well. See Fluid flow with formation damage for more information on damage.

## Inertial-turbulent flow and rate-dependent skin

The diffusivity equation, **Eq. 1**, assumes that Darcy’s law represents the relationship between flow velocity and pressure gradients in the reservoir, an assumption that is adequate for low-velocity or laminar flow. However, at higher flow velocities, deviations from Darcy’s law are observed as a result of inertial effects or even turbulent flow effects. In 1D radial flow, these inertial/turbulent effects (often called non-Darcy flow effects) are confined to the region near the wellbore in which flow velocities are largest. This results in an additional pressure drop similar to that caused by skin, but the additional pressure drop is proportional to flow rate. The apparent skin, *s*′, for a well with non-Darcy flow near the wellbore is given by^{[4]}

where *D* is the non-Darcy flow factor for the system. *D* can be regarded as constant, although, in theory, it depends slightly on near-well pressure. In practice, non-Darcy flow is ordinarily important only for gas wells, which have high-flow velocities near the wellbore, but it can be important for oil wells with high-velocity flow in some situations.

## Radius of investigation and stabilization time

Radius of investigation is the distance a pressure transient has moved into a formation following a rate change in a well. This distance is related to formation rock and fluid properties and time elapsed since a rate change in the well. Consider this concept by visualizing the pressure distributions at increasing times as **Fig. 1** shows for a well producing at constant rate from a reservoir initially at uniform pressure. (These pressure distributions were calculated using the *Ei*-function solution to the diffusivity equation.)

Important observations about this figure include the following:

- The pressure in the wellbore, at
*r*=*r*_{w}, decreases steadily as flow time increases; likewise, pressures at other fixed values of r also decrease as flow time increases. - The pressure drawdown (or pressure transient) caused by producing the well moves further into the reservoir as flow time increases. For the range of flow times shown, there is always a point beyond which the drawdown in pressure from the original value is negligible. This time-dependent point of "negligible drawdown" can be considered to be a radius of investigation.

Analysis shows that the time, *t*, at which a pressure disturbance reaches a distance, *r*_{i}, which is called the radius of investigation, is given by the equation^{[2]}

Investigators differ on the numerical constant in **Eq. 11**, but this difference is of little practical importance if the radius of investigation is used as a semiquantitative indicator of the distance into the reservoir to which formation properties have influenced the response of a well in a pressure-transient test.

The radius of investigation has several applications in pressure-transient test analysis and design. A qualitative use is to help explain the shape of a pressure buildup or drawdown curve. For example, a buildup test plot may have a complex shape at early times when the radius of investigation is in the altered zone near the wellbore, where the permeability is different from formation permeability. Or a buildup test plot may change shape at long times when the radius of investigation reaches the general vicinity of a reservoir boundary.

The radius-of-investigation concept provides a guide for well-test design. For example, you may want to sample reservoir properties at least 1,000 ft from a test well. The radius of investigation concept allows you to estimate the time required to achieve the desired depth of investigation.

**Eq. 11** also provides a means to estimate the time required to achieve "stabilized" flow; that is, the time required for a pressure transient to reach the boundaries of a tested reservoir. For example, if a well is centered in a cylindrical drainage area of radius *r*_{e}, then the time required for stabilization, *t*_{s}, is

For other drainage shapes, the time to stabilization can be quite different, as discussed later.

## Pseudosteady-state flow

The *Ei*-function solution to the radial diffusivity equation is valid only while a reservoir is infinite-acting; that is, until boundaries begin to affect the pressure drawdown at the well. For the constant rate flow of a well centered in its drainage area of radius, *r*_{e}, and modeled by the *Ei*-function solution, these effects begin at *t* = 948 *ϕμc*_{t}*r*_{e}^{2}/*k*. Before these boundary effects, the regime is called unsteady-state flow. After boundary effects are felt fully, the solution to the radial diffusivity equation for a well centered in a cylindrical drainage area and producing at constant rate is^{[2]}

This equation for calculating pressure in the wellbore becomes valid for *t* > 948 *ϕμc*_{t}*r*_{e}^{2}/*k* at the same time at which the *Ei*-function solution becomes invalid.

Another form of **Eq. 13** is useful for some applications. It involves replacing original reservoir pressure, *p*_{i}, with average pressure, , within the drainage volume of the well. The volumetric average pressure within the drainage volume of the well can be found from material balance. The pressure decrease resulting from removal of *qB* RB/D of fluid for *t* hours (a total volume removed of 5.615*qBt*/24 ft^{3}) is

Substituting in **Eq. 13**, the time-dependent terms cancel, and the result is

**Eqs. 13** and **15** are more useful in practice if they include skin factors to account for damage or stimulation. In **Eq. 15**,

## Productivity index

The productivity index, *J*, of an oil well is the ratio of the stabilized rate, *q*, to the pressure drawdown, , required to sustain that rate. For flow from a well centered in a circular drainage area, **Eq. 17** allows us to relate productivity index to formation and fluid properties:

Thus, if a well is tested at several different stabilized rates and the stabilized flowing bottomhole pressure (BHP), *p*_{wf}, is measured at each rate (that is, if pseudosteady-state is attained at each rate), **Eq. 19** implies that a plot of test data should produce a straight line with slope *J* and intercepts *q* = 0 when *p*_{wf} = and when *p*_{wf} = 0. (See **Fig. 2**.) In practice, actual field data will fall below the theoretical straight line for pressures below the bubblepoint pressure of the oil because of increasing gas saturations and oil viscosities that increase the resistance to flow.

## Generalized drainage area shapes

**Eq. 17** is limited to a well centered in a circular drainage area. A similar equation models pseudosteady-state flow in more general reservoir shapes^{[2]}:

where *A* is the drainage area in square feet, and *C*_{A} is the dimensionless shape factor for a specific drainage-area shape and configuration. **Table 1** gives values of *C*_{A}.

The productivity index, *J*, can be expressed for general drainage-area geometry as

Other numerical constants tabulated in **Table 1** allow us to calculate the maximum elapsed time during which a reservoir is infinite-acting (so that the *Ei*-function solution can be used), the time required for the for the pseudosteady-state solution to predict pressure drawdown within 1% accuracy, and time required for the pseudosteady-state solution to be exact. For a given reservoir geometry, the maximum time a reservoir is infinite acting can be determined using the entry in the column "Use Infinite System Solution With Less Than 1% Error for *t*_{DA} <." This *t*_{DA} is defined as 0.0002637*kt*/*ϕμc*_{t}*A*, so this means that the time in hours is calculated from

Time required for the pseudosteady-state equation to be accurate within 1% can be found from the entry in the column titled "Less Than 1% Error for t DA., " Finally, the time required for the pseudosteady-state equation to be exact is found in the entry in the column "Exact for *t*_{DA} >."

**Figs. 3 and 4** show the flow regimes that occur at various times. These figures show *p*_{wf} in a well flowing at constant rate, plotted as a function of time on both logarithmic and linear scales. In the transient region, the reservoir is infinite acting and is modeled by **Eq. 9**, which implies that *p*_{wf} is a linear function of log *t*. In the pseudosteady-state region, the reservoir is modeled by **Eq. 20** in the general case or **Eqs. 15** or **13** for the special case of a well centered in a cylindrical drainage area. **Eq. 13** shows a linear relationship between *p*_{wf} and *t* during pseudosteady-state flow. This linear relationship also exists in generalized reservoir geometries.

At times between the end of the transient region and the beginning of the pseudosteady-state region, there is a transition region, sometimes called the late-transient region. This region is, for practical purposes, nonexistent for wells centered in circular, square, or hexagonal drainage areas, as **Table 1** indicates. However, for a well off-center in its drainage area, the late-transient region can span a significant time region, as **Table 1** also indicates.

## Steady-state flow

Pseudosteady-state flow describes production from a closed drainage area (one with no-flow outer boundaries, either permanent and caused by zero-permeability rock or "temporary" and caused by production from offset wells). In pseudosteady-state, reservoir pressure drops at the same rate with time at all points in the reservoir, including at the reservoir boundaries. Ideally, true steady-state flow can occur in the drainage area of a well, but only if pressure at the drainage boundaries of the well can be maintained constant while the well is producing at constant rate. While unlikely, steady-state flow is conceivable for wells with edgewater drive or in repeated flood patterns in a reservoir. The solution to the radial diffusivity equation is based on a constant-pressure outer boundary condition, instead of a no-flow outer boundary condition. The steady-state solution, applicable after boundary effects have been felt, is

## Constant pressure in the well

Both the steady-state solution (constant pressure at the outer boundaries) and the pseudosteady-state solution (no-flow at the outer boundaries) assume constant rate production in the well. A well is actually more likely to be produced at something close to constant flowing BHP than constant rate. When pressure transients reach no-flow drainage area boundaries, the flow regime is not pseudosteady state; instead, it is more correctly called boundary-dominated flow. If the drainage boundaries are maintained at constant pressure, however, steady-state flow is achieved when the pressure transient reaches the reservoir boundaries.

These different flow regimes are clarified with figures showing pressure distributions in the drainage area of wells with constant flow rate and constant-pressure outer boundaries (**Fig. 5**); constant BHP and constant-pressure outer boundaries (also **Fig. 5**); constant flow rate and no-flow outer boundaries (**Fig. 6**); and constant BHP and no-flow outer boundaries (also **Fig. 6**).

## Wellbore storage

The*Ei-*function solution to the diffusivity equation assumes constant flow rate in the reservoir, starting at time zero. In practice, only the rate at the surface can be controlled. Under ideal conditions, a constant surface rate can be maintained, but the first fluid produced will be fluid that was stored in the wellbore, and, at first, the flow rate from the reservoir into the wellbore will be zero. As the wellbore is unloaded, the reservoir rate approaches the surface rate (**Fig. 7**). Only as the reservoir and surface rates become approximately equal does the*Ei-*function solution become valid. This wellbore unloading during flow tests is a special case of a general phenomenon called wellbore storage.

For a pressure buildup test, the surface rate is zero starting at the instant of shut-in. However, fluid continues to flow into the wellbore from the reservoir because of existing pressure gradients. Idealized models of pressure buildup tests assume a reservoir rate of zero starting at the time of shut-in for the test. This assumption is obviously violated because of the afterflow into the wellbore. As the afterflow rate diminishes, the downhole rate approaches the surface rate (zero), and only as the afterflow rate approaches zero closely can the idealized models closely approximate actual well behavior (**Fig. 8**). Afterflow during buildup tests is another special case of wellbore storage.

The relationship between changes in bottomhole pressure and wellbore unloading or afterflow rates can be modeled with mass balances on the wellbore. There are two special cases of interest: a wellbore completely filled with a single-phase fluid (**Fig. 9**, usually gas in practice) and a wellbore with a rising or falling liquid/gas interface in the well (**Fig. 10**).

For the wellbore filled with a single-phase fluid, ^{[2]}

For a well with a rising or falling liquid/gas interface, ^{[2]}

In most applications, *p*_{t} is assumed to be constant, a convenient but frequently inaccurate simplification. Both equations can be written in the general form

where, for a fluid-filled wellbore,

and, for a moving liquid/gas interface with unchanging surface pressure,

*C* is called the wellbore storage coefficient.

For special cases in which, at earliest times for a flowing well, all the production is coming from fluid stored in the wellbore and none is entering the wellbore from the formation (or, for a shut-in well, the rate of afterflow is equal to the rate before shut in), the integration of **Eq. 26** yields

where Δ*p* is the pressure change in the time because either the start of flow or shut in and Δ*t* is the elapsed time. On a log-log plot of Δ*p* vs. Δ*t* during these early times, a straight line with a slope of unity will result. For any point on this unit slope line, the wellbore storage coefficient, *C*, can be found from any point on the line (Δ*t*, Δ*p*) and **Eq. 29** (**Fig. 11**). Alternatively, the slope (*qB*/24*C*) of a plot of Δ*p* vs. Δ*t* on Cartesian coordinates also leads to an estimate of the wellbore-storage coefficient.

## Linear flow

Linear flow occurs in some reservoirs with long, highly conductive vertical fractures; in relatively long, relatively narrow reservoirs (channels, such as ancient stream beds); and near horizontal wells during certain times. For unsteady-state linear flow in an unbounded (infinite-acting) reservoir, ^{[2]}

## Spherical flow

Spherical flow occurs in wells with limited perforated intervals and into wireline formation test tools. The solution to the spherical/cylindrical, 1D form of the diffusivity equation, subject to the initial condition that pressure is uniform before production and the boundary conditions of constant flow rate and an infinitely large drainage area, is^{[5]}

where ....................(32)

and *r*_{sp} = the radius of the sphere into which flow converges.

## Superposition

The principle of superposition indicates that the total pressure at any point in a reservoir is the sum of the pressure drops at that point caused by flow in each of the wells in the reservoir. A simple illustration of this principle is the case of three wells in an infinite reservoir. Consider wells A, B, and C, that start to produce at times *t*_{A}, *t*_{B}, and *t*_{C} in an infinite-acting reservoir (**Fig. 12**). Application of the principle of superposition shows that^{[2]}

For an infinite-acting reservoir, use the*Ei-*function solutions, including the logarithmic approximation at Well A:

where *t*_{A}, *t*_{B}, and *t*_{C} are times at which wells A, B, and C will begin to produce. The skin factor for Well A is included in **Eq. 29**. The skin factors for other wells are not, because skin factors for individual wells affect only pressures measured inside altered zones for those wells.

Next, consider the use of superposition to model the effects of boundaries in bounded reservoirs. Consider the well in **Fig. 13**, a distance *L* from a single no-flow boundary (such as a sealing fault). Mathematically, this problem is identical to the problem of a well at distance 2*L* from an "image" well; that is, a well that has the same production history as the actual well. The reason that the two-well system simulates the behavior of a well near a boundary is that a line equidistant between the two wells can be shown to be a no-flow boundary. That is, along this line the pressure gradient is zero, which means that there can be no flow. Thus, this problem is a simple problem of two wells in an infinite reservoir:

The drawdown term of the image well does not include a skin factor.

As examples, extend the imaging technique to model wells between boundaries intersecting at 90° (**Fig. 14**); wells between two parallel boundaries (**Fig. 15**); wells near single constant-pressure boundaries (**Fig. 16**); and wells at various locations in closed reservoirs (**Fig. 17**).

One of the most frequently used applications of superposition is to model variable-rate production. Consider **Fig. 18**, in which a well in an infinite-acting reservoir produces at rate *q*_{1} from time 0 to time *t*_{1}; *q*_{2} from *t*_{1} to *t*_{2}, and *q*_{3} for times greater than *t*_{2}. To model the total drawdown for *t* > *t*_{2}, add three drawdowns: the drawdown because of a well producing at rate *q*_{1} starting at time zero and continuing to produce to time *t*; the drawdown because of a well producing at rate (*q*_{2} – *q*_{1}), starting at time *t*_{1} and continuing to time *t*; and the drawdown because of a well producing at rate (*q*_{3} – *q*_{2}) starting at time *t*_{2} and continuing to time *t*. The total drawdown is thus

Horner^{[6]} proposed a convenient alternative to superposition to model the many changes in rate in the history of a typical well. With this approximation, the sequence of*Ei *functions reflecting rate changes can be replaced with a single*Ei *function that contains a single producing time and a single producing rate. The single rate is the most recent nonzero rate at which the well has produced, *q*_{n}. The single producing time, called *t*_{p}, is the ratio of cumulative production, *N*_{p}, to *q*_{n}.

This approximation preserves the material balance in the drainage area of the well and properly gives greatest weight to most recent rate (as opposed to average rate), which dominates the pressure distribution near a well out to the radius of investigation achieved while the well was produced at rate *q*_{n}. The approximation is particularly useful for hand calculations. Given the widespread availability of computer software for analyzing flow and buildup tests on well, the use of more rigorous superposition to model variable-rate production histories is generally more appropriate.

## Semilog methods for flow tests

The logarithmic approximation to the*Ei-*function solution can be used as a basis for analysis of an ideal constant-rate flow test in a well. Written in terms of log_{10}, this equation, which models the BHP for a well in a homogeneous-acting formation with an infinite-acting drainage area and, in absence of wellbore unloading, becomes

This expression has the same form as the equation of a straight line, *y* = *mx* + *b*, with the analogies

These analogies suggest a graphical method of analysis. **Eq. 38** indicates that a plot of *p*_{wf} vs. log_{10}(*t*) should be a straight line with slope m that will allow an estimate of effective permeability to the single liquid phase flowing. (See **Fig. 19**.)

From the intercept, *b*, at *t* = 1 hr [log_{10}(1) = 0], *p*_{1hr}, calculate the skin factor.

In these equations, the slope, *m*, is given by

**Eq. 45** indicates that *m* is most easily determined by choosing values of times *t*_{1} and *t*_{2} that differ by powers of 10 and is especially easy if *t*_{1} and *t*_{2} differ by one log cycle. The intercept, *p*_{1hr}, is the pressure at a time of 1 hour on the best straight line through the data. It may be necessary to extrapolate the straight line to a time of one hour to read the intercept.

## Semilog methods for pressure buildup test

Consider the rate history for an idealized pressure test shown in **Fig. 20**. A well is produced at constant rate *q* for a time *t*_{p}, and then the well shut in (*q* = 0) for a pressure buildup test. The rate history is modeled as the sum of two constant flow rate periods, one at rate *q*, beginning at *t* = 0, and the other at rate –*q*, beginning at *t* = *t*_{p}, at which the time elapsed since shut-in, Δ*t*, is zero. Use the log approximation to the*Ei-*function solution to model the drawdown, and sum them as **Fig. 21** shows. Represented mathematically, the superposition process is

This can be simplified to

Like the drawdown equation, **Eq. 43** can be interpreted as the equation of a straight line. The analogies are

The group [(*t*_{p} + Δ*t*)/Δ*t*] is called the Horner time ratio (HTR) or sometimes simply the Horner time. Our simple model, which describes a buildup test in a homogeneous, infinite-acting reservoir, a well with one constant rate before shut in and without afterflow (wellbore storage), indicates that a graph of *p*_{ws} vs. the HTR should fall on a straight line. From the slope, *m*, of this line, the permeability to the single-phase liquid flowing into the wellbore can be estimated. The intercept, *b*, at log 10 [(*t*_{p} + Δ*t*)/Δ*t*] = 0 or [(*t*_{p} + Δ*t*)/Δ*t*] = 1 provides an estimate of original drainage area pressure, *p*_{i}.

Obtain the slope, *m*, from

where {[ (*t*_{p} +Δ*t*)/Δ*t*]_{1}, *p*_{ws1}} and {[ (*t*_{p} +Δ*t*)/Δ*t*]_{2}, *p*_{ws2}} are any two points on the straight-line (**Fig. 22**). Normally, choose [*t*_{p} + Δ*t*)/Δ*t*]_{1} and [(*t*_{p} + Δ*t*)/Δ*t*]_{2} to be powers of 10.

In **Fig. 22**, which is called a Horner plot, the HTR on the horizontal axis decreases from left to right, so that shut-in time increases from left to right. In some Horner plots, the HTR increases from left to right; in that case shut-in time increases from right to left.

Skin factor can be estimated from a pressure buildup test, even though the skin factor does not appear in the buildup equation, **Eq. 47**. Simultaneously solve the equation modeling the drawdown at the instant of shut in (at time *t*_{p}) with **Eq. 47**, discard terms that are ordinarily negligible, and arrive at the result

The radius-of-investigation concept is also useful for pressure buildup tests, as **Fig. 23** illustrates. The approximate position of the point at which the pressure has built up to a uniform level intersects the region in which the pressure is little affected by the shut-in is given by **Eq. 11**, with elapsed time, *t*, interpreted as shut-in time, Δt.

## Nomenclature

## References

- ↑
^{1.0}^{1.1}Matthews, C.S. and Russell, D.G. 1967. Pressure Buildup and Flow Tests in Wells, Vol. 1. Richardson, Texas: Monograph Series, SPE. - ↑
^{2.0}^{2.1}^{2.2}^{2.3}^{2.4}^{2.5}^{2.6}^{2.7}Lee, W.J. 1982. Well Testing. Dallas, Texas: Textbook Series, SPE. - ↑ Hawkins, M.F.J. 1956. A Note on the Skin Effect. J Pet Technol 8 (12): 65–66. SPE-732-G. http://dx.doi.org/10.2118/732-G
- ↑ Wattenbarger, R.A. and Ramey Jr., H.J. 1968. Gas Well Testing With Turbulence, Damage and Wellbore Storage. J Pet Technol 20 (8): 877-887. http://dx.doi.org/10.2118/1835-PA
- ↑ Joseph, J.A. and Koederitz, L.F. 1985. Unsteady-State Spherical Flow with Storage and Skin. SPE J. 25 (6): 804–822. SPE-12950-PA. http://dx.doi.org/10.2118/12950-PA
- ↑ Horner, D.R. 1967. Pressure Buildup in Wells. Proc., Third World Pet. Cong., The Hague (1951) Sec. II, 503–523; also Pressure Analysis Methods, 9, 25–43. Richardson, Texas: Reprint Series, SPE.

## Noteworthy papers in OnePetro

Use this section to list papers in OnePetro that a reader who wants to learn more should definitely read

## External links

Use this section to provide links to relevant material on websites other than PetroWiki and OnePetro

## See also

Flow equations for gas and multiphase flow

Fluid flow with formation damage

Fluid flow in horizontal wells

Fluid flow in hydraulically fractured wells