Diffusion Equation

Key Takeaways

• The diffusion equation predicts the radius of investigation during a pressure transient test — the distance to which a measurable pressure perturbation has propagated from the wellbore by a given time is approximately r = 0.029 × sqrt(k × t / (φ × μ × cₜ)), where r is in feet, k is in millidarcies, t is in hours, and the other terms are in consistent units; this radius of investigation formula tells the engineer how much of the reservoir is being "seen" by the pressure transient at any given time during the test, and establishes the minimum test duration required for the pressure transient to reach a specific reservoir feature of interest (a fault, a gas-water contact, a permeability barrier); in a high-permeability reservoir with k = 100 md, a 24-hour buildup test may investigate several thousand feet from the well; in a tight gas sand with k = 0.01 md, the same test may investigate only a few hundred feet, and months of transient testing would be required to see boundary effects from the drainage area that the well is actually exploiting.
• The logarithmic pressure derivative — the product of time and the pressure derivative with respect to the natural log of time — is the most powerful diagnostic tool derived from the diffusion equation solutions — on a log-log plot of pressure change and pressure derivative versus elapsed time, each flow regime prescribed by the diffusion equation leaves a characteristic signature: radial flow appears as a flat (horizontal) derivative, linear flow appears as a half-slope (slope 0.5), bilinear flow appears as a quarter-slope, spherical flow appears as a negative half-slope, and boundary effects appear as a unit slope (closed boundary) or a zero-slope (constant pressure boundary); the ability to identify these signatures on the derivative plot allows reservoir engineers to diagnose the flow geometry around the well (is it radial? linear, as in a fractured well? does it hit a fault?) and to extract reservoir properties from specific time windows corresponding to each regime; this derivative diagnostic methodology, developed in the 1980s primarily by Bourdet and collaborators, transformed pressure transient analysis from an art into a quantitative discipline.
• Superposition in time extends the diffusion equation to handle variable rate production histories — the diffusion equation is linear, meaning its solutions can be added (superposed) to handle wells that have produced at multiple different rates rather than a single constant rate; when a well is shut in for a buildup test after a period of production at a varying rate, the buildup pressure response is the superposition of the drawdown response from each rate period, with each period contributing its own diffusion equation solution weighted by the rate change at that time; this superposition principle is implemented in the Horner time ratio (for simple single-rate-then-shut-in cases) and in the equivalent time transformations (superposition time, multirate superposition) used for complex rate histories; accurate rate history is therefore as important as accurate pressure data for buildup test interpretation, because errors in the rate history distort the superposed pressure response and lead to incorrect reservoir property estimates.
• Tight reservoirs and shale formations require modified diffusion equation solutions to account for non-Darcy flow and dual-porosity media — in ultra-low-permeability reservoirs, gas flow near the wellbore occurs at high velocity relative to the Darcy-flow assumption, creating inertial pressure drop (Forchheimer flow) that the linear diffusion equation doesn't capture; additionally, naturally fractured tight reservoirs and shale gas formations exhibit dual-porosity behavior where gas stored in tight matrix blocks desorbs and diffuses into a connected natural fracture network on different timescales than flow through the fractures themselves; the Warren-Root dual-porosity model extends the diffusion equation to capture this behavior, introducing an interporosity flow term that creates characteristic "V-shaped" pressure derivative signatures during buildup tests; for shale gas wells with hydraulic fractures, complex fracture network flow models that combine diffusion through the matrix, linear flow through the hydraulic fractures, and radial flow in the reservoir beyond fracture tips require multi-region diffusion equation solutions that are implemented in commercial pressure transient analysis software.
• The diffusion equation connects pressure transient testing to reservoir simulation through the same governing physics — reservoir simulators discretize the diffusion equation (along with equations for multiphase flow, compositional effects, and gravity) into finite difference or finite element grids and solve numerically for pressure and saturation distributions throughout the reservoir over time; pressure transient analysis solves analytical solutions to the diffusion equation to extract average reservoir properties from wellbore pressure measurements; the two approaches are complementary: pressure transient analysis provides calibration data (permeability, skin, drainage area) for the simulator, while the simulator predicts future performance that can be validated against subsequent production history and additional transient tests; understanding that both approaches are implementations of the same fundamental physics helps engineers integrate them correctly rather than treating simulation and well testing as separate, disconnected disciplines.