Superposition: Definition, Pressure Transient Analysis, and Well Testing
What Is Superposition in Well Testing?
Superposition is a mathematical principle applied in pressure transient analysis (PTA) that allows the pressure response of a reservoir to multiple wells or multiple rate changes to be computed by algebraically summing the individual pressure responses of each well or each rate period, treated as if they were independent. The principle holds because the diffusivity equation governing pressure propagation in a reservoir is linear — allowing individual solutions to be added together. In practice, superposition enables engineers to analyse pressure buildup tests conducted after variable-rate production (using the Horner method extended to multi-rate history), to account for interference between nearby producing wells, and to model the effect of boundaries (via image well superposition). Superposition is the foundation of virtually all analytical pressure transient analysis methods used in reservoir engineering.
Key Takeaways
- Superposition holds because the pressure diffusivity equation is linear — multiple pressure disturbances add independently without interfering with each other.
- Time superposition converts variable-rate production history into an equivalent constant-rate problem using the superposition time function — enabling correct Horner analysis after rate-varying production periods.
- Spatial superposition (method of images) accounts for no-flow boundaries by placing mirror-image wells across the boundary, so the combined solution automatically satisfies the zero-flux boundary condition.
- Multi-well interference is modelled by superimposing each well's pressure response at the observation point — enabling interference test analysis and estimation of inter-well transmissibility.
- Superposition breaks down when the diffusivity equation becomes non-linear — this occurs at high gas flow rates (non-Darcy flow), in reservoirs with pressure-dependent permeability, and in dual-porosity matrix-fracture systems at early time.
Time Superposition and the Horner Method
The Horner method for pressure buildup analysis assumes a constant-rate production period followed by shut-in. In practice, wells produce at variable rates between well tests and workovers. Time superposition corrects for this by converting the actual variable-rate history into an equivalent producing time tp using the superposition principle: tp = Np / qlast, where Np is the cumulative production and qlast is the final production rate before shut-in. This equivalent time is substituted into the Horner time ratio [(tp + Δt)/Δt] — the approximation is exact for a single previous shut-in but is a simplification for complex multi-rate histories. For complex rate histories, the multi-rate superposition time function (a sum of logarithmic terms for each rate step) replaces the Horner approximation and provides an exact analytical solution.
Spatial superposition (image well method) is used to model no-flow boundary effects in well test analysis. When a well produces near a sealing fault, the fault imposes a zero-flux boundary condition. The image well method places a virtual well of equal production rate on the opposite side of the fault, mirroring the real well. The combined pressure response of the real well and its image automatically satisfies the zero-flux condition at the fault. On a pressure buildup log-log derivative, the fault appears as a doubling of the derivative (the derivative doubles to 1.0 from its radial flow value of 0.5) — the distance to the fault can be calculated from the time at which the derivative begins to double.
- Principle basis: linearity of the pressure diffusivity equation — solutions add
- Time superposition use: correct Horner analysis for variable-rate production history
- Spatial superposition use: image wells to model faults and boundaries in PTA
- Multi-well superposition use: interference testing — pressure at observer from multiple producers
- Horner equivalent time: tp = Np / qlast (simplification for single-rate history)
- Fault distance formula: d = 0.0328 × √(k × t_boundary / φμct) in Darcy units
- Failure conditions: non-Darcy flow, pressure-dependent permeability, dual-porosity early time
- Software: Kappa Saphir, Ecrin, Fekete Harmony — all use superposition internally
Always account for recent production rate changes when analysing a pressure buildup — do not assume the Horner single-rate approximation is valid if the well has had rate fluctuations in the days before shut-in. A well that produced at 500 bbl/day for 60 days, then 800 bbl/day for 2 days before shut-in, has a very different equivalent time than the single-rate approximation suggests. The correct superposition time function using the full rate history will show radial flow at a different time and slope than the Horner approximation, and will give different permeability and skin. In Kappa Saphir or similar PTA software, input the full rate history (not just the final rate period) and let the software compute the superposition time function — this is the most common source of error in field BU analyses done with simplified Horner calculations.
Superposition Synonyms and Related Terminology
Superposition in well testing is also referred to as:
- Superposition principle — the underlying mathematical concept
- Method of images — the spatial superposition technique for boundary representation
- Multi-rate analysis — the application of time superposition to variable-rate production histories
- Superposition time function — the combined time variable that replaces clock time in multi-rate PTA plots
Related terms: Pressure Buildup, Flow Regime, Reservoir Pressure, Skin Factor
Frequently Asked Questions About Superposition
How does superposition account for a second well producing nearby?
When an observation well (shut in for a pressure interference test) is surrounded by producing wells, the pressure observed at the observation well is the sum of the pressure drawdown contributions from each surrounding producer — superimposed at the observation well location. The pressure drop contribution from well i at distance r from the observation well is: ΔP_i = (70.6 × q_i × B × μ / kh) × Ei(-948φμctr² / kt), where Ei is the exponential integral and the other terms are reservoir and fluid properties. By monitoring the pressure at the observation well while the surrounding producers follow a known rate schedule, inter-well transmissibility and reservoir continuity between the wells can be determined. Interference testing is the only field measurement that directly confirms hydraulic communication between wells — particularly valuable for identifying compartment boundaries and fault transmissibilities in complex reservoirs.
Why does superposition break down at high gas flow rates?
At high gas velocities in tight formations or near the wellbore, flow departs from Darcy's Law — inertial effects cause a pressure drop proportional to velocity squared (Forchheimer equation), not just to velocity linearly (Darcy). This velocity-squared term makes the governing equation non-linear, violating the linearity assumption that superposition requires. In practice, this non-Darcy flow component appears as an additional "rate-dependent skin" (D × q) in the wellbore equation: total apparent skin = S_mechanical + D × q. At the producing rate q, this rate-dependent skin adds to the mechanical skin from formation damage. A pressure buildup test eliminates rate-dependent skin (since rate is zero), but the drawdown test at producing rate includes it. Comparing buildup and drawdown skin estimates and dividing the difference by q gives the non-Darcy coefficient D — an important design parameter for high-rate gas wells.
What is the method of images and when is it used?
The method of images places virtual "image wells" at the mirror-image locations of real wells across reservoir boundaries to enforce the correct boundary conditions. A no-flow boundary (sealing fault) requires zero flux across it — achieved by placing an equal-rate image producer across the fault from the real well. The combined pressure field of real + image wells automatically has zero gradient at the fault plane. A constant-pressure boundary (aquifer or gas cap) requires zero pressure change — achieved by placing an equal-rate image injector across the boundary. Complex boundary geometries (intersecting faults forming a wedge or rectangle) require multiple image wells, each of which is itself an image of previous image wells — the series is terminated when subsequent image contributions become negligible. The method of images is exact for straight boundaries and provides accurate approximate solutions for irregular boundaries encountered in seismic-mapped fault systems.
Why Superposition Matters in Oil and Gas
Superposition is the mathematical machinery that makes analytical pressure transient analysis work in real wells with variable production rates, nearby wells, and bounded drainage areas. Without superposition, PTA would be restricted to idealised constant-rate single-well systems that do not exist in practice. Every Horner plot analysis, every interference test interpretation, every boundary distance calculation in well testing software implicitly uses superposition. Understanding superposition allows reservoir engineers to correctly apply PTA methods, to recognise when the linearity assumption breaks down (non-Darcy flow, dual-porosity early time), and to design well test programmes that yield interpretable data rather than rate artefacts that masquerade as reservoir responses.