The developed algorithm is applied to generate the dimensionless pressure and the dimensionless pressure derivative for a horizontal well in eight different systems whose analytical solutions are not available. These cases include (1) a horizontal well in a two-layer reservoir, (2) a horizontal well in a two-layer reservoir with gas-cap drive, (3) a horizontal well in a two-layer reservoir with bottom-water drive, (4) a horizontal well in a two-layer reservoir with gas-cap and bottom-water drives, (5) a horizontal well intersecting a two-layer reservoir, (6) a snake-shape horizontal well in a two-layer reservoir, (7) a horizontal well intersecting a three-layer reservoir, and (8) a vertical well intersecting a two-layer reservoir without cross flow. Their transient behaviors are studied in this work.

This dissertation presents a boundary element method for solving pressure transient behavior of multi-layer reservoirs. Two types of boundary conditions, Dirichlet and Neumann, can be handled by the proposed technique. The boundary element method yields several advantages over the conventional finite-different and finite-element methods; including the reduction of the dimension of the problem by one, yielding more accurate results by using the analytical solutions of the governing equations as the weighting functions, no gird orientation and numerical dispersion effects, and flexibility of handling complex geometries and boundary conditions.

The step-by-step procedures for calculating reservoir parameters are developed in this work. The direct synthesis is applied to interpret pressure transient behavior of a horizontal well in complicated systems without type-curve matching. The conventional log-log, semi-log, and the Cartesian plots are used in this work. The procedure is illustrated by numerical examples.

The developed algorithm is successfully applied to some simple cases with known analytical solutions including a finite linear aquifer, a horizontal well in closed reservoir, and a horizontal well in a multi-layer reservoir.

Pressure transient analysis for horizontal wells has gained a lot of attention recently due to the increasing number of horizontal wells. Interpretation of horizontal well test data is much more difficult than for a vertical well. Due to heterogeneity, application of single-layer models fails to interpret the reservoir properties. Numerical methods are essential for a horizontal well in multi-layer reservoir with complicated wellbore configuration and boundary conditions.