Linear Flow: Definition, Fractured Well Analysis, and Unconventional Reservoir Production

What Is Linear Flow?

Linear flow is a flow regime in which reservoir fluids move in parallel, straight-line paths toward the wellbore or fracture — as opposed to radial flow, where fluid converges from all directions around a cylindrical wellbore. In a hydraulically fractured vertical well, linear flow occurs when the pressure transient has propagated away from the fracture face but has not yet reached the fracture tip or the reservoir boundary — fluid flows perpendicularly from the low-permeability matrix into the high-permeability fracture plane, and from the fracture into the wellbore. In a horizontal well (with or without hydraulic fractures), linear flow describes fluid flowing perpendicularly from the formation into the horizontal drain along the full length of the lateral. Linear flow is the dominant production mechanism in hydraulically fractured tight gas and shale oil wells for the majority of their producing life — the low matrix permeability (<0.1 md) means radial flow to the fracture, which takes years to decades to develop in conventional wells, is replaced by linear flow from the matrix into the fracture system. Identifying and characterising linear flow on production data plots is the primary method for estimating stimulated reservoir volume (SRV) and original gas in place (OGIP) in unconventional wells.

Key Takeaways

  • Linear flow is identified by a half-slope (slope = 0.5) on a log-log plot of rate vs time, or equivalently a straight line on a square-root-of-time plot (√t vs normalised pressure) — the half-slope diagnostic is the primary tool for identifying linear flow on production data.
  • In a hydraulically fractured well, linear flow to the fracture (matrix linear flow) is characterised by the product √k × x_f × h (square root of permeability times fracture half-length times pay thickness) — estimating this product from the linear flow slope is how fracture half-length x_f is derived from production data analysis.
  • In unconventional shale and tight gas wells, linear flow often persists for years to decades — the low matrix permeability means the pressure transient propagates extremely slowly, and the well may never reach boundary-dominated flow (pseudo-steady state) within its economic life.
  • Multiple linear flow regimes occur in hydraulically fractured horizontal wells: early fracture linear flow (flow within the fracture from tip to wellbore), formation linear flow (flow from matrix into the fractures), and late-time compound linear flow (when all fractures drain simultaneously from a laterally extensive connected volume).
  • End of linear flow (ETLF) is the most important diagnostic event in unconventional well analysis — it marks when the pressure transient reaches the boundary of the stimulated volume, and allows estimation of the drained area (and therefore OGIP) from the time at which the linear flow ends.

Linear Flow Physics and Diagnostic Plots

In linear flow, the governing pressure diffusivity equation reduces to 1D — flow is perpendicular to the fracture plane, all streamlines parallel. Wellbore flowing pressure varies proportionally to the square root of time: ∆p = p_i − p_wf ∝ √(πk)/(x_f × h × √(φμc_t)) × q√t, where k is matrix permeability, x_f is fracture half-length, h is pay thickness, φ is porosity, μ is viscosity, and c_t is total compressibility. This square-root-of-time response means: (1) on a log-log diagnostic plot, linear flow appears as slope = 0.5 (half-slope); (2) on a Cartesian plot of normalised pressure (∆p/q) vs √t, linear flow is a straight line. In rate-transient analysis (RTA), the equivalent is rate normalised pressure (RNP = ∆p/q) vs material balance time (t_mb = Np/q) — linear flow is a half-slope on log-log RNP vs t_mb, or a straight line on √t_mb vs RNP.

The transition from linear flow to boundary-dominated flow (BDF) marks the end of linear flow — when the pressure transient has fully propagated through the stimulated reservoir volume. Cumulative production at ETLF is approximately equal to the initial gas or oil in the drained volume, allowing drainage volume estimation. Identifying end of linear flow on production data plots provides a direct estimate of how much reservoir each well has drained — and by subtraction, how much remains between wells accessible by infill drilling.

Fast Facts: Linear Flow
  • Diagnostic signature: half-slope (0.5 slope) on log-log plot of rate-normalised pressure vs material balance time; straight line on √time vs RNP Cartesian plot
  • Duration in unconventional wells: months to decades, depending on matrix permeability (lower k = longer linear flow period)
  • Dominant flow regime in: hydraulically fractured tight gas, shale oil (Permian Wolfcamp, Eagle Ford, Bakken), shale gas (Marcellus, Haynesville, Duvernay, Montney)
  • Formation linear flow slope intercept: proportional to √k × x_f × h — the "linear flow parameter" from which fracture size and permeability are estimated
  • End of linear flow time: t_ETLF ∝ (x_f²φμc_t)/k — longer in low-k formations; marks when pressure transient reaches the drainage boundary
  • Other linear flow regimes: bilinear flow (flow in both the fracture and formation simultaneously, slope = 0.25), compound linear flow (multiple fractures draining cooperatively at late time)
  • Rate-transient analysis (RTA) tools: IHS Harmony (now IHS Markit/Enverus), Kappa Topaz, PHDwin (TRC) — all include linear flow analysis modules specifically for unconventional wells
  • Key insight: if a well's entire production history plots as a half-slope line, the well has never reached its drainage boundary — the SRV is larger than the estimated drainage area, or the fractures are not as large as originally designed
Production Engineering Tip:

Plot all your unconventional well production data on log-log rate-normalised pressure vs material balance time plots as the primary diagnostic — not rate vs time or GOR vs cumulative. The half-slope linear flow signature is the most reliable indicator of hydraulic fracture performance and fracture half-length, and the end-of-linear-flow event is the only point at which you can directly estimate drainage area from production data alone. If a well that should reach the end of linear flow (based on fracture spacing, stage count, and estimated fracture half-length) never shows the transition from half-slope to unit-slope (boundary-dominated flow), one of two things is true: the fractures are smaller than designed (and the stage spacing is not fully draining the interwell region) or the matrix permeability is so low that linear flow will persist for decades and the well will decline on a gradual half-slope trend without ever reaching pseudo-steady state in its economic life. Distinguishing between these two scenarios requires comparing wells drilled at different spacings — if tighter spacing wells do show earlier end of linear flow, fracture half-length was limited, and the infill wells are accessing undrained rock. If tighter spacing does not accelerate end of linear flow, the fractures are large and the interwell region is already being drained.

Linear flow is also referred to as:

  • Formation linear flow — specifically the phase where fluid flows from matrix into fractures; distinguished from fracture linear flow (flow within the fracture plane itself) and bilinear flow (simultaneous linear flow in both matrix and fracture)
  • Half-slope — the diagnostic signature on the log-log diagnostic plot (slope = 0.5) that identifies linear flow; "the well is showing half-slope" is equivalent to "the well is in linear flow"
  • Square-root-of-time plot — the Cartesian plot of normalised pressure vs √t used to identify and quantify linear flow; the straight line through this plot characterises the linear flow parameter
  • Transient linear flow — distinguishes the transient (before boundary effects) linear flow regime from compound linear flow (after interference between fracture stages)

Related terms: Hydraulic Fracturing, Pressure Transient Analysis, Rate Transient Analysis, Tight Gas

Frequently Asked Questions About Linear Flow

How is linear flow used to estimate fracture half-length?

Fracture half-length (x_f) is estimated from the linear flow parameter F_CD derived from the slope of the normalised pressure vs √time Cartesian plot. If matrix permeability k is independently known (core measurements, PTA), x_f is the only unknown and is solved directly. In practice, matrix permeability in unconventional wells is often unknown (ranging 0.0001–0.01 md), so only the product √k × x_f can be determined from the slope alone — meaning fracture half-length estimates always carry significant uncertainty without an independent permeability measurement. Some operators use microseismic monitoring to estimate x_f independently, then use the linear flow slope to back-calculate k. For a multi-stage well with n fractures, the total linear flow parameter is n × x_f × h × √k — stage count and fracture height from design parameters help resolve the uncertainty.

What is bilinear flow and how does it differ from linear flow?

Bilinear flow occurs in finite-conductivity hydraulic fractures — it is characterised by simultaneous linear flow in two perpendicular directions: formation-to-fracture (perpendicular to the fracture plane) and fracture-to-wellbore (parallel to the fracture, tip toward wellbore). Both flows are active simultaneously because the fracture resistance to flow is significant, causing pressure to drop noticeably along the fracture length. The diagnostic signature is a quarter-slope (0.25) on the log-log diagnostic plot — distinctly shallower than the half-slope (0.5) of pure linear flow. Bilinear flow precedes linear flow at early times when fracture resistance dominates, transitioning to linear flow as the pressure transient grows deep enough that formation resistance dominates. Infinite-conductivity fractures (no fracture resistance) skip bilinear flow entirely. Most real shale and tight gas fractures have finite conductivity and show bilinear flow before the linear flow half-slope.

How does linear flow duration relate to well spacing in unconventional plays?

Linear flow duration is directly related to stimulated reservoir volume relative to well spacing — the key diagnostic for whether spacing is optimal, too tight (causing interference), or too wide (leaving undrained rock). Linear flow ends when the pressure transient reaches a no-flow boundary — the physical reservoir edge or the interference front from an adjacent well. Widely-spaced wells (660 ft+ in tight formations) may stay in linear flow for their entire economic life. Closely-spaced wells (250–330 ft) terminate linear flow early as interference from adjacent wells takes hold. End-of-linear-flow time scales as t_ETLF ∝ (x_f)² / k — in a 0.0001 md reservoir with 300 ft fracture half-length, ETLF may take 10–20 years; in a 0.01 md reservoir, 3–6 months. If all wells reach BDF within 2–3 years without overlapping drainage areas, spacing is too wide and infill wells can access undrained rock.

Why Linear Flow Matters in Oil and Gas

Linear flow is the most commercially significant production mechanism in unconventional oil and gas development — it describes the fundamental mode by which hydraulically fractured tight rock delivers hydrocarbons to the wellbore during the majority of a well's producing life. Understanding linear flow is essential for every aspect of unconventional well performance analysis: estimating EUR from production data, calibrating fracture designs against actual performance, optimising well spacing to balance early-time performance against long-term recovery, and deciding whether infill drilling will access undrained rock or simply cannibalise production from existing wells. The Permian Basin, Marcellus Shale, Montney Formation, Eagle Ford, and Bakken — collectively responsible for millions of barrels per day of production — are all developed based on the economic logic of sustainable linear flow from hydraulically created fracture systems. The ability to diagnose linear flow correctly, measure its duration, and extract fracture parameters from the diagnostic plots separates high-performance operators from average performers in unconventional development — a difference that compounds across thousands of wells into billions of dollars of EUR difference at the portfolio level.