Unsteady-State Flow

Unsteady-state flow, also called transient flow, is a reservoir or wellbore flow condition in which pressure and flow rate are changing with time at every point in the reservoir because pressure disturbances have not yet reached the drainage boundaries, with the radius of investigation expanding outward from the wellbore as the square root of time and pressure behavior governed by the diffusivity equation, providing the theoretical basis for pressure transient analysis methods including drawdown and buildup tests used to determine reservoir permeability, skin factor, and drainage area.

Key Takeaways

During the infinite-acting radial flow (IARF) period, pressure at the wellbore declines as the logarithm of time, and the semi-log slope of the pressure-time relationship is inversely proportional to permeability-thickness (kh product), forming the basis for the conventional straight-line analysis used in pressure transient testing.
Wellbore storage, caused by fluid compressibility in the wellbore volume after shut-in, masks the early-time formation response in buildup tests and must be recognized and accounted for before the true reservoir response is interpretable; the end of wellbore storage distortion is identified on the log-log diagnostic plot as the departure from the unit-slope early-time trend.
The Horner plot (pressure versus the Horner time ratio log[(tp+delta-t)/delta-t]) linearizes the pressure buildup response during infinite-acting radial flow and extrapolates to the initial reservoir pressure (p*) when the time ratio approaches unity, giving both permeability from the slope and initial pressure without requiring long shut-in times.
Boundary effects transition the unsteady-state period to pseudo-steady-state (closed boundaries) or steady-state (constant-pressure boundaries) when the pressure disturbance reaches the drainage boundary; identification of boundary effects from the late-time behavior of a pressure transient test defines the drainage area and shape.
In tight gas and shale gas wells, unsteady-state production decline persists for years to decades because the extremely low matrix permeability means the pressure front advances so slowly that drainage boundaries are not reached within commercial production lifetimes, making transient decline analysis the relevant production forecasting tool.

Fast Facts

The radius of investigation (ri) during a pressure transient test is approximately ri = 0.032 x sqrt(k x t / (phi x mu x ct)), where k is permeability in millidarcies, t is time in hours, phi is porosity, mu is viscosity in centipoise, and ct is total compressibility in psi^-1. For a well in a 10 md conventional reservoir, the radius of investigation reaches 300 metres after approximately 24 hours. For a 0.001 md shale well, the same radius requires nearly 300,000 hours (34 years) of transient flow, illustrating why shale wells produce in transient flow throughout their economic life. Pressure buildup tests typically run for 1-10 times the production time to adequately sample the reservoir.

Tip: When interpreting a pressure buildup test, always construct the log-log diagnostic plot of pressure change and pressure derivative versus elapsed shut-in time before performing semi-log or Horner analysis. The derivative plot is the most powerful diagnostic tool in pressure transient analysis because it identifies flow regimes visually: a horizontal derivative indicates radial flow (use for kh), a half-slope indicates linear flow (hydraulic fracture or channel), and a unit slope indicates wellbore storage or a closed boundary. Attempting semi-log analysis without first confirming that the radial flow regime is present leads to incorrect permeability estimates.

What Is Unsteady-State Flow

In reservoir engineering, flow regimes are classified by whether the pressure distribution in the reservoir changes with time. Unsteady-state (transient) flow occurs when the pressure disturbance caused by production or injection is still propagating outward through the reservoir and has not encountered any boundary. During this period, every point in the reservoir is at a different pressure from its condition at the previous instant, and the entire pressure profile is evolving. This regime begins when a well starts producing and ends when the radius of investigation reaches the drainage boundary.

This contrasts with steady-state flow, where pressure at every point remains constant with time because a constant-pressure boundary (such as a strong aquifer or an injection well) supplies fluid at the same rate as production withdraws it. The intermediate case, pseudo-steady-state flow, occurs when a well has reached all closed boundaries of its drainage volume simultaneously, after which pressure declines uniformly throughout the drainage area at a constant rate proportional to the production rate divided by the drainage volume.

How Unsteady-State Flow Works

The mathematical description of unsteady-state flow in a porous medium is the diffusivity equation, which in radial coordinates for a slightly compressible fluid is: (1/r) d/dr [r dP/dr] = (phi mu ct / k) dP/dt. This equation states that the spatial pressure curvature drives the temporal pressure change, with the rate of diffusion governed by the hydraulic diffusivity kappa = k / (phi mu ct). A high diffusivity (high permeability, low compressibility) means pressure disturbances travel quickly through the reservoir, and transient behavior is short-lived. A low diffusivity (tight rock, high compressibility) means pressure disturbances move slowly, and transient behavior dominates production for a long time.

For a well producing at constant rate with no wellbore storage, the analytical solution to the diffusivity equation at early time (after the logarithmic approximation becomes valid, approximately when kt/phi mu ct rw^2 exceeds 100) gives the famous logarithmic pressure response: Pwf = Pi - (162.6 q mu B / kh) x [log(t) + log(k / phi mu ct rw^2) - 3.2275 + 0.86859 S], where S is the skin factor. This equation is the basis for the conventional straight-line analysis on a semi-log plot: the slope m = 162.6 q mu B / kh yields permeability-thickness, and the intercept at one hour gives skin factor.

Wellbore storage complicates early-time buildup analysis: after surface shut-in, fluid continues flowing from the reservoir into the wellbore, masking the formation response. On the log-log diagnostic plot, wellbore storage manifests as a unit slope in both pressure change and derivative. The end of wellbore storage is identified as the departure from unit slope, followed by the horizontal derivative plateau characteristic of infinite-acting radial flow.

Pressure derivative analysis, introduced by Bourdet and colleagues in 1983, revolutionized pressure transient interpretation by plotting the logarithmic derivative dP/d(ln delta-t) alongside the pressure change on the same log-log axes. Each flow regime has a characteristic derivative signature: radial flow shows a horizontal derivative at the level kh/141.2 q mu B; linear flow (from a hydraulic fracture or linear channel) shows a half-slope derivative; bilinear flow (from a finite-conductivity fracture) shows a quarter-slope; spherical flow shows a negative half-slope; wellbore storage shows a unit slope; and closed boundary pseudo-steady state shows a unit slope at late time. This diagnostic power allows simultaneous identification of flow regimes, reservoir heterogeneity, fracture characteristics, and boundaries from a single test.

Unsteady-State Flow Across International Jurisdictions

In Canada, the AER governs pressure transient testing under Directive 040 (Pressure and Deliverability Testing Oil and Gas Wells), requiring buildup tests on new wells in specified formations. The AER's reserves evaluation guidelines, published jointly with the Saskatchewan and BC regulators, accept rate-transient analysis (RTA) in unsteady-state flowing conditions as a valid reserves estimation method for tight Montney and Duvernay wells where conventional buildup tests require impractically long shut-in times.

In the United States, the SEC's 2009 Modernization of Oil and Gas Reporting Rules explicitly accepted rate-transient analysis for tight gas and shale wells, acknowledging that conventional well testing in sub-millidarcy formations is economically impractical. BSEE requires bottomhole pressure surveys and deliverability tests for OCS gas wells under 30 CFR Part 250 Subpart K. The SPE's Petroleum Reserves and Resources Definitions (2018) provide the global technical framework for applying transient flow concepts to reserves classification.

In Norway, Sodir requires well test data demonstrating reservoir productivity and connectivity in plan for development and operation (PDO) submissions under the Petroleum Resources Act. The complex geology of North Sea chalk and turbidite reservoirs, where natural fractures and compartmentalization strongly affect transient responses, has driven investment in pressure derivative interpretation workflows combining multi-well interference testing with single-well transient analysis.

In the Middle East, Saudi Aramco's reservoir surveillance for the Ghawar, Shaybah, and Khurais mega-fields includes systematic pressure buildup and falloff test programs to monitor connectivity and water front advancement. The Arab-D formation produces characteristic dual-porosity transient signatures from matrix-fracture systems that require specialized workflows. Kuwait Oil Company and ADNOC maintain multi-decade transient test archives used for history matching of full-field simulation models guiding billion-dollar production optimization.

Unsteady-state flow is synonymous with transient flow in most reservoir engineering usage. The specific period when no boundary has been reached is called infinite-acting radial flow (IARF). The transition regime after boundaries are partially sensed is sometimes called transient boundary-dominated flow. Related concepts include pressure buildup test, the standard well test method used to analyze unsteady-state behavior during shut-in; the Horner plot, the semi-log analysis tool for buildup data; skin factor, the dimensionless near-wellbore damage or stimulation parameter estimated from pressure transient analysis; and permeability, the primary formation property measured from the slope of the radial flow period. Pseudo-steady-state flow is the regime that follows unsteady-state once all drainage boundaries have been reached. Rate-transient analysis (RTA) applies unsteady-state principles to long-term production decline data for tight and shale wells.

FAQ

How long must a well be shut in to reach the end of infinite-acting radial flow in a tight formation?
The time required to reach a specific radius of investigation scales as ri^2 x phi x mu x ct / k. For a tight Montney well with 0.01 md permeability, 7 percent porosity, 0.02 cp gas viscosity, and 1x10^-4 psi^-1 total compressibility, reaching 300 metres radius of investigation requires approximately 3,000 hours (125 days) of shut-in. This is economically impractical for most operators, which is why rate-transient analysis using long production histories under flowing conditions has replaced conventional buildup tests in many tight formation development programs.

What is the difference between pseudo-steady-state and unsteady-state production decline?
Unsteady-state decline occurs while the pressure disturbance is still propagating outward and no boundary has been reached. Production rate versus time follows a pattern governed by the diffusivity equation, and the decline exponent b in Arps' decline curve analysis approaches 2 for linear flow (hydraulic fractures) or 0.5 for transitional radial flow periods. Pseudo-steady-state decline occurs after all boundaries are sensed simultaneously, and pressure declines uniformly across the drainage area; the Arps decline exponent b approaches 0 (exponential decline) for a single closed-boundary reservoir producing under depletion. Recognizing which regime applies is essential for selecting the correct decline curve analysis method for production forecasting.

Why Unsteady-State Flow Matters

Unsteady-state flow theory is the foundation of pressure transient analysis, one of the most powerful tools for characterizing reservoir properties without core measurement. Extracting permeability, skin factor, fracture length, drainage area, and reservoir pressure from a single buildup test is unique to transient analysis and irreplaceable for appraisal wells with sparse formation data. For unconventional wells producing in transient flow for their entire economic life, applying exponential decline will dramatically underestimate ultimate recovery. Well spacing, lateral length, and stimulation design decisions in major shale plays rely on accurate unsteady-state characterization using rate-transient analysis.