Line Source Solution
The line source solution is the analytical solution to the radial diffusivity equation that describes pressure transient behavior in a porous medium around a vertical wellbore that is modeled as an infinitely thin line source (or line sink) of fluid production or injection — derived by treating the wellbore as a mathematical line of infinitesimal radius rather than a cylinder of finite dimension, and applying the Laplace transform and inverse transform methods to the governing partial differential equation of pressure diffusion in a homogeneous, isotropic, single-phase reservoir, the line source solution provides the pressure distribution as a function of radial distance from the well and time in terms of the exponential integral function Ei(-x), which evaluates to a simple logarithmic approximation at small values of x (corresponding to the condition of radial flow established more than a few hours to days after well opening); in well test analysis, the line source solution is the mathematical foundation for the conventional semilog straight-line method (the Horner buildup plot, the Miller-Dyes-Hutchinson method, and related techniques) in which the logarithmic approximation to the Ei function produces a linear relationship between shut-in pressure and the logarithm of the Horner time ratio, with the slope of that line yielding the formation permeability and the extrapolated pressure providing an estimate of average reservoir pressure; the assumptions underlying the line source solution (infinite homogeneous reservoir, single-phase, slightly compressible flow, constant-rate production, no wellbore storage or skin) define the regime in which the semilog straight-line analysis is valid and the conditions under which more complex analytical or numerical models are required.
Key Takeaways
- The exponential integral Ei(-x) function that appears in the line source solution simplifies to a logarithmic approximation when the argument x is less than approximately 0.01, which corresponds to the condition that sufficient time has elapsed after well opening for the pressure transient to have propagated far enough from the wellbore to be in the radial flow regime — mathematically, this logarithmic approximation is: p(r,t) = pi - (qBmu/4*pi*kh) * [ln(4kt/phi*mu*ct*r^2) - 0.80907], where the term in brackets becomes the familiar semilog straight line when r is fixed at the wellbore radius and time t is varied; this logarithmic approximation is extremely accurate for well test analysis purposes once the wellbore storage effect and skin-affected near-wellbore period have ended and the reservoir is in radial flow, typically after one to two log cycles of time on the log-log pressure derivative plot; the practical implication is that when you see a flat derivative plateau on the log-log diagnostic plot (indicating radial flow), you are in the regime where the line source solution's logarithmic approximation is valid and the semilog straight-line method will give accurate permeability and skin estimates.
- Interference testing between two wells uses the full Ei function form of the line source solution rather than the logarithmic approximation, because the observation well at which the pressure response is measured may be at a radial distance from the active well where the argument x is not small enough for the logarithmic approximation to apply — in a well pair with 1,000 feet of spacing, the early time pressure response at the observation well is governed by the Ei function behavior where x is of order 1 or larger, and the logarithmic approximation would be inaccurate; the interference test analysis must use the complete Ei function matched to the observed pressure response at the observation well, and the match parameters (permeability, storativity) characterize the reservoir between the two wells rather than just the near-wellbore properties characterized by a single-well buildup test; the line source solution for interference testing was first formalized by Theis (1935) in groundwater hydrology and adapted for petroleum engineering by Ramey and colleagues in the 1950s and 1960s, demonstrating the deep connection between the mathematical tools used in water well testing and petroleum well test analysis.
- The infinite-acting radial flow period governed by the line source solution is recognizable on the pressure derivative plot as the characteristic flat (zero slope) derivative plateau — when the pressure derivative dp/d(lnt) is constant with time, the well is in radial flow and the reservoir between the wellbore and the outer boundary of the pressure transient is behaving as predicted by the line source solution; this diagnostic identification of radial flow is the first step in any well test analysis workflow and is prerequisite for applying the semilog straight-line method; any deviation from the flat derivative before radial flow (rising derivative indicating wellbore storage, falling derivative indicating a high-conductivity natural fracture or a skin transition) or after radial flow (falling derivative indicating sealing boundaries or interference from adjacent wells, rising derivative indicating a closed reservoir or a permeability decrease) must be identified and accounted for before the semilog method is applied; applying the semilog straight-line to a time period that is not in radial flow gives a permeability estimate that may be off by a factor of two to ten compared to the true formation permeability.
- The skin factor, which quantifies near-wellbore damage or stimulation, is extracted from the line source solution by comparing the measured wellbore pressure at a reference time with the pressure predicted by the undamaged (skin = 0) line source solution for the same formation permeability — the difference between the measured pressure and the undamaged prediction represents an additional pressure drop concentrated at the wellbore face (a delta-p caused by skin), and the skin factor S is defined as the ratio of this additional pressure drop to the pressure drop per unit of logarithmic time interval during radial flow; a positive skin indicates that the near-wellbore permeability is lower than the bulk formation permeability (damage from drilling, cementing, or scale deposition), while a negative skin indicates that the near-wellbore connection is better than the undamaged matrix (natural fractures, stimulated perforations, hydraulic fractures intersecting the wellbore); the quantification of skin from the line source solution provides the engineer with a numerical estimate of how much production is being lost to wellbore damage and what production improvement would result from a successful stimulation treatment that eliminates the skin.
- The line source solution fails to adequately describe wellbore behavior under several conditions that are common in real wells, and recognizing these deviations is as important as knowing when the solution applies — the line source solution assumes radial symmetry (violated by horizontal wells, hydraulically fractured wells, and wells near boundaries), infinite reservoir extent (violated when the pressure transient reaches a boundary before the test ends), homogeneous permeability (violated by naturally fractured reservoirs and layered systems), and constant rate production (violated during rate changes, pump-off events, and during wellbore storage); when any of these assumptions is violated, the semilog straight-line method based on the line source solution gives incorrect permeability estimates, and more sophisticated analytical models (dual-porosity models for fractured reservoirs, linear flow solutions for horizontal wells, composite models for layered systems) or numerical simulation must be used to match the observed pressure behavior and extract reliable formation parameters; the history of well test analysis since the 1960s is largely the history of recognizing when the line source solution is insufficient and developing the appropriate alternative models to handle the conditions it cannot describe.
Fast Facts
The mathematical foundation of the line source solution in petroleum engineering traces directly to work done by the physicist George Green in 1828 and to Charles-Eugène Delaunay's 19th-century heat conduction solutions. The petroleum-specific form was independently developed for groundwater by Charles Vernon Theis in 1935 (the "Theis equation" in hydrogeology is mathematically identical to the line source solution) and adapted for petroleum reservoirs by A.F. van Everdingen and W. Hurst in their landmark 1949 SPE paper "The Application of the Laplace Transform to Flow Problems in Reservoirs." That paper established the mathematical framework for pressure transient analysis that underlies every buildup test, drawdown test, and interference test performed globally for the subsequent 75 years.
What Is the Line Source Solution?
The line source solution is the mathematical engine behind most of what engineers read from a pressure buildup test. It describes how pressure spreads through a reservoir from a producing well, assuming the well is a geometric line (infinitely thin) rather than a cylinder with physical dimensions. That simplification — treating the wellbore as a line rather than a hole — turns out to work remarkably well at the radial distances and time scales relevant to well testing, because the wellbore radius is tiny compared to the radius of investigation of a pressure transient. The solution says that the pressure at any point in the reservoir depends on the logarithm of time when the transient is in radial flow, and that logarithmic relationship is exactly the straight line that Horner plots exploit to extract permeability. It is a century-old piece of mathematical physics, derived from heat conduction theory, applied to the economics of oil and gas reservoirs. The fact that it still works is testimony to how well the simplifying assumptions of radial flow in a homogeneous reservoir represent real geological behavior across a remarkable range of formation types.
Synonyms and Related Terminology
The line source solution is also called the Ei solution (referring to the exponential integral function it contains), the Theis equation (in groundwater engineering), or the Ei-function solution. Related terms include radial flow (the flow geometry that the line source solution describes, in which fluid converges toward the wellbore from all directions equally), pressure transient well tests (the application domain in which the line source solution is the primary analytical tool), Horner plot (the semilog analysis method based on the logarithmic approximation to the line source Ei function), skin (the near-wellbore damage or stimulation parameter derived by comparing measured pressure with the line source solution prediction), interference testing (the multi-well application of the full Ei function form of the line source solution), and radial flow regime (the specific time period during a well test when the line source logarithmic approximation is valid).
Why the Line Source Solution Is Still the Starting Point After 75 Years
More sophisticated reservoir models exist. Dual-porosity solutions handle naturally fractured reservoirs. Horizontal well solutions account for linear and bilinear flow periods. Boundary-dominated solutions describe depleting compartments. All of them build on the line source solution as their foundation, either by modifying its assumptions or by using it as the far-field boundary condition for a more complex near-wellbore model. When you look at a pressure derivative plot and see that flat plateau indicating radial flow, you are watching the line source solution manifest itself in real data. When you draw the Horner straight line and read off the slope, you are using the mathematical result that van Everdingen and Hurst published in 1949. The durability of the solution comes from its accuracy in the regime where its assumptions hold, and in the middle time range of most well tests, that regime is where the data live. It does not cover everything, and knowing its limits is as important as knowing its mathematics. But as a starting point for understanding what a reservoir is telling you through a pressure transient test, the line source solution remains unequaled.