SPE Reservoir Evaluation & Engineering
Volume 10, Number 2, April 2007, pp. 122-131

SPE-92041-PA

Grid-System Requirements in Numerical Modeling of Pressure-Transient Tests in Horizontal Wells

View full textPDF ( 1,719 KB )

DOI  More information 10.2118/92041-PA http://dx.doi.org/10.2118/92041-PA

Citation

  • Al-Mohannadi, N., Ozkan, E., and Kazemi, H. 2007. Grid-System Requirements in Numerical Modeling of Pressure-Transient Tests in Horizontal Wells. SPE Res Eval & Eng  10 (2): 122-131. SPE-92041-PA.

Discipline Categories

  • 0
  • 6.6.4 Drillstem/Well Testing
  • 6.3.1 Flow in Porous Media

Summary

The common approach of using refined grids and small timesteps usually does not provide accurate and efficient numerical models to simulate horizontal-well pressure-transient tests. In most cases, this approach requires impracticably long computational times or creates artifacts on pressure and derivative responses that may be confused with the characteristics of certain flow regimes. This paper presents an investigation of the grid and timestep requirements for the conventional well index (WI) to provide accurate simulations of horizontal-well pressure-transient responses. The results of this study indicate that the conventional WI should be used with log-distributed grid and timesteps to obtain a closer match to the analytical solution at early times. There is a limit to grid refinement beyond which the results do not improve. For these cases, the transient well index (TWI) introduced in this paper allows larger grids while improving the accuracy of simulations, especially at early times.

Introduction

The need for numerical simulation of the pressure-transient behavior of horizontal wells often arises in the interpretation of complex structures and the calibration of static data used in reservoir models. To compute the bottomhole pressures, reservoir simulators require a relationship between the wellblock and wellbore pressures. This relationship is obtained by considering an appropriate well-deliverability equation, which expresses the flow rate as a function of the difference between the wellblock and wellbore pressures. The pressure calculated in numerical simulation is uniform within the wellblock, whereas deliverability equations consider the change of pressure as a function of distance from the wellbore. However, it is possible to find an equivalent distance in the wellblock at which the pressure corresponds to the simulated wellblock pressure. Therefore, the choice of the appropriate deliverability relationship and equivalent wellblock radius is essential for accurate numerical calculation of wellbore pressures. The selected deliverability relationship is invoked into the numerical algorithm in the form of WI.

Van Poollen et al. (1968) and Peaceman (1978, 1983) presented the first comprehensive discussions of the WI to simulate wellbore pressures. Peaceman (1978, 1983) introduced the probe radius formula to represent the equivalent wellblock radius in the WI. Peaceman’s approach has been commonly accepted and widely used in most numerical simulators.

Steady-state and radial-flow assumptions are inherent in Peaceman’s numerical WI. When this conventional WI concept is applied to wells with nonradial-flow regimes (e.g., wells off-centered in the grid and inclined or horizontal wells at late times) and used to calculate pressure-transient responses, the accuracy of wellbore pressures may not be satisfactory. The common approach for the solution of this problem is to reduce the grid size (and the timesteps) so that the flow convergence within smaller wellblocks becomes closer to the radial and steady-state assumption is more reasonable because of the smaller storativity of the wellblocks.

Peaceman’s probe radius formula (Peaceman 1978, 1983) has been commonly used for horizontal wells without much critical assessment of success for early times. As an alternative, Babu et al. (1991a, 1991b) presented another WI and equivalent wellblock radius formula for horizontal wells. Their formula extended Peaceman’s conventional WI to wellblocks with high aspect ratio and high anisotropy.

Blanc et al. (1999) proposed techniques including the use of a TWI to improve the grid-system requirements in modeling pressure responses of vertical wells, especially at early times. They implied the extension of their technique to horizontal wells but did not provide conclusive information.

In this paper, our objective is to provide a means of improving pressure-transient modeling in standard simulators. We intend to respond to two different categories of the modeler’s needs. If the modeler has the option of modifying WI in the simulator, then we recommend the use of the TWI as introduced for horizontal wells in this paper. If an existing (commercial) simulator is to be used without revising WI, then we provide recommendations for the selection of the grid system and timesteps.

View full textPDF ( 1,719 KB )

History

  • Original manuscript received: 13 September 2004
  • Meeting paper published: 8 November 2004
  • Revised manuscript received: 10 February 2007
  • Manuscript approved: 24 February 2007
  • Version of record: 20 April 2007