Superposition in Space: Image Wells, No-Flow Boundaries, and Multi-Well Pressure Interference
Superposition in space is the mathematical technique reservoir engineers use to add together the pressure responses of several wells, or of one real well and one or more imaginary image wells, to compute how production from multiple sources affects pressure at any chosen point and time in a reservoir. It works because the diffusivity equation that governs single-phase, slightly compressible flow through a porous medium is linear, meaning the individual pressure-drop solutions for each source can simply be summed at a given location without coupling terms. The most familiar application in well testing is representing a sealing fault or other impermeable barrier that sits near a producing well. A real no-flow boundary is reproduced exactly by placing an image well, an imaginary producer flowing at the same rate, mirror-positioned across the barrier at an equal perpendicular distance on the far side. The pressure interference from that image well forces the net flux across the barrier plane to zero, which is the mathematical definition of a no-flow surface. A constant-pressure boundary, such as a strong edge aquifer or a gas cap, is built the same way but with an image injector of equal and opposite rate, so the pressure on the boundary plane stays pinned at its initial value. By combining several image wells, complex shapes such as a well in a rectangular fault block, a channel sand, or the corner between two intersecting faults can be assembled from simple line-source solutions, a process that may require an infinite array of images for a fully bounded reservoir. Superposition in space is the spatial companion to superposition in time, which handles variable flow rate by stacking time-shifted rate steps. Together they let analysts interpret real pressure data through the lens of an idealized model, then read boundary distances directly from the timing of the derivative response. In the Western Canadian Sedimentary Basin, where multi-well pads on the Montney and Duvernay place six to sixteen horizontal wells within a single drilling spacing unit, spatial superposition is also the engine behind interference and frac-hit analysis: each offset producer or injector contributes a measurable pressure signal at the observation well, and summing those contributions reproduces the observed transient. Without superposition, every bounded or multi-well problem would require a fresh analytical or numerical solution; with it, the entire catalogue of classical well-test models is built from one fundamental building block, the radial line-source solution evaluated at a calculated distance.
Key Takeaways
- Linearity Is The Prerequisite: Spatial superposition is only valid because the diffusivity equation is linear for slightly compressible single-phase flow. That lets per-well pressure drops add directly. For strongly pressure-dependent gas viscosity-compressibility, engineers first linearize using pseudo-pressure and pseudo-time before superposing, otherwise the additive assumption breaks down in tight Montney and Duvernay gas.
- No-Flow Boundary Uses A Like-Sign Image: A sealing fault is modeled by an image producer of identical rate placed symmetrically across the fault. The doubling of pressure drop after the transient reaches the image produces the classic factor-of-two slope change on a semilog plot, from which distance to the fault d is read using the radius-of-investigation timing.
- Constant-Pressure Boundary Uses An Opposite-Sign Image: An aquifer or gas cap is represented by an image injector of equal and opposite rate, holding the boundary plane at initial pressure. The pressure derivative falls toward zero and stabilization appears, the opposite signature to a sealing fault, which doubles the derivative.
- Multi-Well Interference Sums Real Sources: On a crowded WCSB pad, the pressure at an observation well equals the initial pressure minus the summed drawdown contributions of every active offset well, each evaluated at its own distance and elapsed time. This is the basis of interference tests, pulse tests, and quantitative frac-hit diagnosis under AER Directive 083 spacing reviews.
- Image Arrays Build Complex Geometry: A well between two parallel faults, inside a fault-bounded block, or in a channel sand requires multiple or infinite image wells reflected repeatedly across each boundary. Modern PTA software assembles these arrays automatically, but the underlying idea remains stacking simple line sources at computed mirror distances.
Reading Fault Distance From The Derivative Doubling
When a horizontal Montney producer drilled near a sealing fault is shut in for a buildup, the pressure derivative first stabilizes on the radial flow plateau, then climbs to roughly twice that level once the transient reaches the image well. The time of that doubling fixes the distance d to the fault through the radius-of-investigation relationship, d being proportional to the square root of the product of permeability, time, and the reciprocal of porosity-viscosity-compressibility. A typical result might place a fault 180 m (about 590 ft) from the lateral, a number an operator like ARC Resources or Tourmaline would cross-check against 3D seismic before committing the next infill well, since drilling into compartment-limiting structure changes recoverable reserves and well economics directly.
Frac Hits And Parent-Child Interference
In stacked Duvernay and Montney development, a child well completed beside a depleted parent often registers an immediate pressure response, a frac hit, when its fracture network intersects the parent's drainage. Spatial superposition quantifies this: the observed pressure change at the parent equals the summed transient from each child stage treated as a source at its mapped location. Engineers use the magnitude and timing to estimate effective fracture half-length and spacing efficiency, feeding AER Directive 083 well-spacing and well-density arguments. A pad showing strong superposed interference at 250 m (820 ft) spacing may be downspaced no further, protecting per-well recovery and avoiding sterilized reserves.
Fast Facts
The image-well concept predates digital reservoir simulation by decades. It was borrowed directly from the method of images in classical electrostatics and heat conduction, where a charge or heat source near a boundary is mirrored to satisfy the boundary condition. Lord Kelvin formalized the electrostatic version in 1848, and petroleum engineers adapted the identical mathematics in the mid-twentieth century because the steady and transient flow equations share the same Laplacian form as electrostatic potential and Fourier heat flow, a striking case of one trick solving three unrelated physical problems.
Related Terms
Spatial superposition pairs with Pressure Transient Analysis, the broader discipline of reading reservoir properties from buildup and drawdown data, because every bounded-reservoir interpretation rests on summed image-well responses. It complements Radius of Investigation, which converts the timing of a boundary signature into a physical distance. The Diffusivity Equation supplies the linear governing physics that makes the additive sum legitimate, and Skin Factor often appears alongside in the same buildup analysis as a near-wellbore correction layered onto the superposed model.
Real-World WCSB Scenario: A Montney Buildup Near A Sealing Fault
A Tourmaline Montney horizontal in the Alberta Deep Basin is shut in for a 96-hour buildup after a year on production. The pressure derivative stabilizes near 1,400 kPa (203 psi) on the radial plateau, then doubles to roughly 2,800 kPa (406 psi) after about 30 hours. Plugging permeability of 0.08 mD, porosity of 6 percent, and the elapsed time into the radius-of-investigation equation places a no-flow boundary about 165 m (541 ft) from the lateral, matching a fault mapped on reprocessed 3D seismic. Engineering cost of the test, including deferred production, runs near CAD 240,000.
The confirmed compartment boundary led the operator to shift the planned infill 200 m updip into the open block rather than drilling toward the fault, preserving an estimated 0.9 Bcf (about 25 e3m3 times 10^3) of incremental recoverable gas. The superposition-based interpretation, costing a fraction of the well, directly steered a multi-million-dollar capital decision.