S.E. Brame, J.W. Castle, SPE, O.K. Fawumi*, SPE and R.W. Falta, Clemson University

1. History Matching of Heavy Oil Production for Comparing New Approaches to Generating Reservoir Property Distributions, West Coalinga Field, California S.E. Brame, J.W. Castle, SPE, O.K. Fawumi*, SPE and R.W. Falta, Clemson University * Now with Mobil Producing Nigeria Unltd. SPE 93469

2. This research was funded by the U.S. Department of Energy, Fossil Energy Oil Technology Program, through the National Petroleum Technology Office under contract number DE-AC26-98BC15119.

3. Objectives • Assess the suitability of different permeability distributions generated from geological and fractal modeling. • Examine the feasibility of creating realistic permeability distributions from fractal theory for reservoir simulation. • Use steam flood simulations of a portion of the West Coalinga oil field in California for model assessment.

4. Location of the West Coalinga Field

5. Heavy Oil Sands of Coalinga • The West Coalinga oil field in California produces from heavy oil sands of the Miocene Temblor Formation. • The oil in this field has low API gravity (12o to 15o API) and is highly viscous (900 cp) at the 40oC initial reservoir temperature. • This makes the field an ideal candidate for enhanced oil recovery through steam injection. Steam injection well at Coalinga

6. Methods • Characterize geology in outcrop and core • Identify lithofacies • Group lithofacies into Facies Groups • Identify Facies Tracts • Characterize Fractal Facies • Construct geologic models • Assign permeability distributions • Simulate steam flooding of different models and assess results

7. Geologic Characterization • Coalinga offers a unique opportunity to observe and characterize the reservoir rocks on the surface, immediately adjacent to the oil production area. • Thus it was possible to relate cores of the producing formation to nearby outcrops. Clemson Students examining outcrop of Temblor Fm. in hills north of Coalinga Field

8. Lithofacies Characterization • A total of fifteen lithofacies were identified from outcrops and cores based on observed lithological differences. • Statistical methods were used to consolidate the 15 different lithofacies groups into five facies groups.

9. Facies Tract Characterization Five facies tracts were interpreted based on detailed sedimentological analysis:

10. Facies CorrelationsBetween Wells Well 118A BS = bounding surface (depths are in feet) Correlation of well 239 with Type Well 118A to determine the location of bounding surfaces, facies tracts, and facies groups.

11. Model AreaThis map shows the 3 adjacent five-spot configurations that were modeled and the wells that were used. P r o d u c t i o n W e l l I n j e c t i o n W e l l 3 0 0 8 0 0 8 - 1 3 0 0 6 0 0 2 3 9 2 2 9 2 3 9 W 2 2 9 W 3 0 0 4 0 0 2 2 8 - 2 B 2 3 8 3 0 0 2 0 0 2 3 8 W 2 2 8 W 2 2 8 3 0 0 0 0 0 2 3 8 A 8 - 2 2 9 9 8 0 0 1 2 8 B 1 2 8 1 1 8 B 2 9 9 6 0 0 Northing 2 3 7 2 2 7 2 9 9 4 0 0 8 - 3 2 3 7 W 2 9 9 2 0 0 1 2 7 1 2 7 B 2 9 9 0 0 0 2 3 6 1 1 8 A 8 - 4 2 9 8 8 0 0 2 3 6 W 0 0 2 9 8 6 1 5 8 7 5 0 0 1 5 8 7 7 0 0 1 5 8 7 9 0 0 1 5 8 8 1 0 0 1 5 8 8 3 0 0 1 5 8 8 5 0 0 Section 36D E a s t i n g

12. Geologic Model Construction • 3-D geologic models were constructed using GOCAD. • Inputs included the bounding surface horizons, geophysical logs, and facies group and facies tract data. • The grid consisted of 9600 cells: - 32 vertical layers, - 10 cells in the x direction - 30 cells in the y direction

13. Geologic Models in GOCAD Facies Group 1 Facies Group 2 Facies Group 3 Facies Group 4 Facies Group 5 200 ft Facies Tract Facies Group

14. Fractal Group Model Development • The cores of five West Coalinga wells were analyzed foot by foot. • Analysis of the individual facies data revealed that a unique Gaussian fractal structure was present in each one. • These results and others led to the development of a new model for representing natural heterogeneity called the fractal/facies concept.

15. Reservoir Simulation • The geologic model grids were used as the framework for the flow simulation mesh. • Petrophysical properties were assigned to all cells of the mesh. • Permeability distributions of the facies tract, facies group, and fractal group models were assigned. • Numerical simulations of steam injection were used to assess the different models.

16. T2VOC Flow Simulator • Uses a general finite difference formulation. • Can solve multi-phase, multi-component mass and energy balance equations. • Has been used to simulate a variety of subsurface processes such as: • nonaqueous phase liquid (NAPL) migration, • soil vapor extraction, • air-sparging, • steam injection, and • direct pumping of water and NAPL.

17. Permeability Distribution of the Facies Tract Model

18. Permeability Distribution of the Facies Group Model

19. Fractal Permeability Distribution • Fractal permeabilities were stochastically generated for the 5 different facies groups. • The values were assigned to the flow grid but using a finer grid (~3,000,000 cells). • The fractal permeabilities were upscaled to the standard flow grid cell size (9600 cells). • An arithmetic mean was used for the horizontal permeability. • A harmonic mean was used for the vertical permeability.

20. Permeability Distribution of the Fractal Model

21. Simulation Period • A five year period (Oct. 1995 to Oct. 2000) was used. • The injected steam volume changed monthly. • Some injection wells were not online until 1997. • All changes in production and injection were honored. Volume of Steam (bbls of water) Time (years)

22. Parameters Adjusted for Better Fit • The oil-water relative permeability curves provided for the field were based on a data fit from a core. • The initial oil saturation was interpolated from well logs. • Problem: these values resulted in simulations where the water to oil ratio was off by a factor of 10 or more compared to field values!

23. Solution: The water relative permeability endpoint was reduced from .56 to .15, and the oil saturations were increased everywhere by 20% (upper limit of 70%). Normalized Water-oil Relative Permeabilities ■ Kro Field Data ● Krw Field Data Kro used in Model Krw used in Model Kro Krw Water Saturation (%)

24. Reservoir Temperatures at 5 Years Fractal Group Facies Tract Facies Group

25. Oil Saturations at 5 Years Fractal Group Facies Tract Facies Group

26. Simulated versus Field Production of Oil (Combined Production)

27. Simulated versus Field Production of Water (Combined Production)

28. Summary • Three models were assessed using simulations of steam injection and heavy oil production. • Lowering the water relative permeability curve endpoint decreased water production. • Increasing the initial oil saturation resulted in a better prediction of oil production. • The fractal realization matched oil production up to Year 3. • Additional realizations could be generated to improve the history match.

29. Conclusions • Permeability distributions from the geologic models provided reasonable matches of oil production. • The fractal/facies approach is a viable method of generating geologically realistic permeability distributions. • Steamflood simulations demonstrate that the fractal/facies approach is feasible for modeling heavy oil reservoirs.

30. Acknowledgments • Chevron contributed core data, production data, and geophysical logs. • Many thanks to Venton Shoemaker, George Anderson, Louis Klonsky, Paul Henshaw, and Mike Clark of Chevron. • Fred Molz and Silong Lu of Clemson University performed the fractal analysis • Ray Christopher of Clemson University assisted with the statistical treatment of the lithofacies.