Semianalytical Stream-Function Solutions on Unstructured Grids for Flow in Heterogeneous Media
- Randy Doyle Hazlett (Potential Research Solutions) | D. Krishna Babu (Potential Research Solutions) | Larry Wayne Lake (University of Texas at Austin)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- June 2007
- Document Type
- Journal Paper
- 179 - 187
- 2007. Society of Petroleum Engineers
- 5.6.5 Tracers, 5.7.2 Recovery Factors, 5.1.1 Exploration, Development, Structural Geology, 5.5 Reservoir Simulation, 5.5.2 Core Analysis, 5.5.7 Streamline Simulation, 5.1 Reservoir Characterisation, 4.1.2 Separation and Treating, 4.1.5 Processing Equipment, 4.3.4 Scale, 5.6.4 Drillstem/Well Testing, 5.3 Reservoir Fluid Dynamics, 5.3.2 Multiphase Flow
- 3 in the last 30 days
- 526 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 10.00|
|SPE Non-Member Price:||USD 30.00|
This paper outlines a Boundary Element Method (BEM) for a piece-wise analytic solution of the Laplace (Poisson) equation for pseudosteady-state, single-phase flow on unstructured, rectangular grids. The method models flow through a reservoir that has been segmented into interacting homogeneous rectangular regions; no further discretization of the solution space analogous to grid refinement in numerical schemes is required for improved accuracy. Rather, boundary discretization allows for continuation of pressure and flux. Previous work on pressure distribution modeling is extended to analytically capture the stream function. Stream-function solutions can then form the basis for other performance measures, such as improved oil recovery efficiency estimation or tracer flow analysis. Moving beyond structured grids into unstructured grid geometry allows for advanced flexibility in problem development and improved efficiency in solution construction. The analytic approach avoids the need for numerical differentiation of the pressure field and particle tracking methods to recover streamlines. Capturing flow in highly heterogeneous media, without local grid refinement, is demonstrated to showcase the robustness of the technique in handling complex reservoir architecture, of particular interest in optimal well positioning and optimal well-pattern development.
The solution to fluid flow problems is typically a map of the driving force (i.e., potential or pressure). A more intuitive result is a map that shows actual trajectories of fluid elements, the stream function (Muskat 1937). Potential, F, and stream-function, ?, are related in 2D by their spatial gradients.
Curves of constant stream-function value are the so-called streamlines. Stream-function indexing is associated with integration, as the difference in stream-function indices represents the amount of fluid flowing between streamlines of those fixed values. By definition, no fluid convection occurs across a streamline. Unfortunately, the concept of stream function is restricted to two dimensions. While an orthogonal mesh to the pressure field can certainly be constructed in 3D, an equivalent of stream function is found in three dimensions only for flows exhibiting symmetry properties that effectively reduce the dimensionality of the flow problem.
|File Size||1 MB||Number of Pages||9|
Archer, R.A. and Horne, R.N. 2000. The Green Element Method forNumerical Well Test Analysis. Paper SPE 62916 presented at the SPE AnnualTechnical Conference and Exhibition, Dallas, 1-4 October. DOI:10.2118/62916-MS.
Babu, D.K. and Odeh, A.S. 1989. Productivity of a HorizontalWell. SPERE 4 (4): 417-421. SPE-18298-PA. DOI:10.2118/18298-PA.
Batycky, R.P., Blunt, M.J., and Thiele, M.R. 1997. A 3D Field-Scale Streamline-BasedReservoir Simulator. SPERE 12 (4): 255-263. SPE-36726-PA.DOI: 10.2118/36726-PA.
Bratvedt, F. et al. 1993. FRONTLINE and FRONTSIM. Two Full-Scale, Two-Phase,Black Oil Reservoir Simulators Based on Front Tracking. Surv. Math. Ind.3: 185-215.
Brown, P.N. and Hindmarsh, A.C. 1989. Reduced Storage Matrix Methods inStiff ODE Systems. J. Appl. Math. & Comp. 31: 40-91.
Deeds, N.E., McKinney, D.C., and Pope, G.A. 2000. LaboratoryCharacterization of Non-Aqueous Phase Liquid/Tracer Interaction in Support of aVadose Zone Partitioning Interwell Tracer Test. J. of ContaminantHydrology 41: 193-204.
Ding, Y. 1999. Using BoundaryIntegral Methods to Couple a Semi-Analytical Reservoir Flow Model and aWellbore Flow Model. Paper SPE 51898 presented at the SPE ReservoirSimulation Symposium, Houston, 14-17 February. DOI: 10.2118/51898-MS.
Duff, G.D.F. and Naylor, D. 1966. Differential Equations of AppliedMathematics. New York City: John Wiley & Sons.
Edwards, M.G., Agut, R., and Aziz, K. 1998. Quasi K-Orthogonal Streamline Grids:Gridding and Discretization. Paper SPE 49072 prepared for presentation atthe SPE Annual Technical Conference and Exhibition, New Orleans, 27-30September. DOI: 10.2118/49072-MS.
Emanuel, A.S. and Milliken, W.J. 1997. Application of Streamtube Techniquesto Full-Field Waterflood Simulation. SPEREE 12 (3): 211-217.SPE-30758-PA. DOI: 10.2118/30758-PA.
Hazlett, R.D. and Babu, D.K. 2005. Optimal Well Placement inHeterogeneous Reservoirs Through Semi-Analytic Modeling. SPEJ10 (3): 286-296. SPE-84281-PA. DOI: 10.2118/84281-PA.
Higgins, R.V. and Leighton, A.J. 1962. A Computer Method To CalculateTwo-Phase Flow in Any Irregularly Bounded Porous Medium. JPT14 (6): 679-683; Trans., AIME, 225. SPE-243-PA. DOI:10.2118/243-PA.
Jain, V., Bryant. S.L., and Sharma, M.M. 2003. Influence of Wettability andSaturation on Liquid-Liquid Inerfacial Area in Porous Media. Environ. Sci.Technol. 37 (3): 584-591.
Jasti, J.K., Penmatcha, V.R., and Babu, D.K. 1999. Use of Analytical Solutions toImprove Simulator Accuracy. SPEJ 4 (1): 47-56. SPE-55225-PA.DOI: 10.2118/55225-PA.
Jongkittinarukorn, K. and Tiab, D. 1998. Development of the Boundary ElementMethod for a Horizontal Well in Multilayer Reservoir. Paper SPE 39939presented at the SPE Rocky Mountain Regional/Low Permeability ReservoirsSymposium, Denver, 5-8 April. DOI: 10.2118/39939-MS.
King, M.J. and Datta-Gupta, A. 1998. Streamline Simulation: A CurrentPerspective. In Situ 22: 91-140.
Kikani, J. and Horne, R.N. 1993. Modeling Pressure-Transient Behaviorof Sectionally Homogeneous Reservoirs by the Boundary-Element Method.SPEFE 8 (2): 145-152. SPE-19778-PA. DOI: 10.2118/19778-PA.
Kirkup, S. 2001. The BEM for Helmholtz Problems.http://www.boundary-element-method.com/bemhelm/manual.htm.
Koh, L.S. and Tiab, D. 1993. ABoundary-Element Algorithm for Modeling 3D Horizontal Well Problems Using 2DGrids. Paper SPE 26228 presented at the SPE Petroleum Computer Conference,New Orleans, 11-14 July. DOI: 10.2118/26228-MS.
Linton, C.M. 1999. Rapidly Convergent Representations for Green's Functionsfor Laplace's Equation. Proc. R. Soc. Lond. A 455:1767-1797.
McCann, R C., Hazlett, R.D., and Babu, D.K. 2001. Highly AccurateApproximations of Green's and Neumann Functions. Proc. R. Soc. A457: 767-772.
Mahadevan, J., Lake, L.W., and Johns, R.T. 2003. Estimation of True Dispersivity inField-Scale Permeable Media. SPEJ 8 (3): 272-279.SPE-86303-PA. DOI: 10.2118/86303-PA.
Marshall, S.L. 1999. A Rapidly Convergent Modified Green's Function forLaplace's Equation in a Rectangular Region. Proc. R. Soc. Lond. A455: 1739-1766.
Martin, J.C. and Wegner, R.E. 1979. Numerical Solution of Multiphase, 2DIncompressible Flow Using Streamtube Relationship. SPEJ 19(5): 313-323; Trans., AIME, 267. SPE-7140-PA. DOI:10.2118/7140-PA.
Masukawa, J. and Horne, R.N. 1988. Application of the Boundary IntegralMethod to Immiscible Displacement Problems. SPERE 3 (3):1069-1077. SPE-15136-PA. DOI: 10.2118/15136-PA.
Muskat, M. 1937. The Flow of Homogeneous Fluids Through Porous Media.New York City: McGraw-Hill.
Numbere, D.T. and Tiab, D. 1988. An Improved Streamline-GeneratingTechnique That Uses the Boundary (Integral) Element Method. SPERE3 (3): 1061-1068. SPE-15135-PA. DOI: 10.2118/15135-PA.
Odeh, A.S. and Babu, D.K. 1990. Transient Flow Behavior of HorizontalWells: Pressure Drawdown and Buildup Analysis. SPEFE 5 (1):7-15; Trans., AIME, 289. SPE-18802-PA.DOI:10.2118/18802-PA.
Portella, R.C.M. and Hewett, T.A. 1998. Upscaling, Gridding, and SimulationUsing Streamtubes. Paper SPE 49071 prepared for presentation at the SPEAnnual Technical Conference and Exhibition, New Orleans, 27-30 September. DOI:10.2118/49071-MS.
Sato, K. and Horne, R.N. 1993. Perturbation Boundary Element Methodfor Heterogeneous Reservoirs: Part 1—Steady-State Flow Problems.SPEFE 8 (4): 306-314. SPE-25299-PA. DOI: 10.2118/25299-PA.
Seager, M.K. 1989. A SLAP for the Masses. In Parallel Supercomputing:Methods, Algorithms and Applications, ed. G.F. Carey. New York City: Wiley,135-155.
Stroud, A.H. 1971. Approximate Calculation of Multiple Integrals. EnglewoodCliffs, New Jersey: Prentice-Hall.
Wen, X-H., Deutsch, C.V., and Cullick, A.S. 2003. Inversion of DynamicProduction Data for Permeability: Fast Streamline-Based Computation ofSensitivity Coefficients of Fractional Flow Rate. J. Hydrology281: 296-312.
Zemel, B. 1995. Tracers in the Oil Field. New York City:Elsevier.