sIMULACIÓN YACIMIENTOS
sIMULACIÓN YACIMIENTOS
sIMULACIÓN YACIMIENTOS
Simulation Studies
D. G. Harris, SPE-AIME, Exxon Production Research Co.
Introduction
This paper describes the geological activities required for Reservoir Simulation Models
constructing realistic mathematical reservoir models. The Over the past few years, improved technology4,S has made
methods used by a geologist are pointed out and' the it possible to develop detailed and complex computer
typical data provided by these methods are presented. programs for simulating fluid flow in new and in de-
Also, the geological controls on reservoir properties and veloped reservoirs. These dynamic models are based on
the significance of these controls in reservoir performance the Darcy equation for fluid flow and on material-balance
are summarized. The information presented should be calculations. The simulation may be in one, two, or three
helpful to both the engineer and the geologist involved dimensions and may involve one, two, or three fluid
with reserves estimation and performance prediction. phases. Engineering literature contains further informa-
Reservoir simulation models are used increasingly to tion on the variety and historical development of simula-
evaluate both new and mature fields,2,4 and to determine tion models. 4,5,11
the most efficient management scheme. Part of the infor- A common simulation model is the two-dimensional,
mation contained in the model is provided by the multicell cross-section shown in Fig. 1. Such a model is
geologist, based on his studies of the physical framework capable of handling vertical variations in porosity, per-
of the reservoir. To be useful, such studies mustdevelop meability, and capillary properties for a limited (or typi-
quantitative data. The geologist must also understand the cal) area. The variation is represented by assigning values
type of engineering data needed and how the engineer can to the essential parameters (Table 1) characterizing each
aid the quantification of the framework data. cell or element. These and other data are entered into a
In the following discussion, reservoir simulation mod- computer program, and changes in fluid volumes, fluid
els are described briefly. This provides a base for discus- state, etc., are calculated for each cell on a predetermined
sion of the geologist's activities in simulation work and production (or injection) schedule.
the assistance that he requires from the engineer. This The distributions of the reservoir- and nonreservoir-
discussion is restricted to the reservoir-description phase rock types and of the reservoir fluids determine the
of simulation studies. geometry of the model and influence the type of model to
The main topic is demonstration of the geological ac- be used. For example, the number and scale of the shale
tivities used in a study of a pilot-test site. The illustrations (or dense carbonate) breaks in the physical framework
indicate the information that is often required for simula- determine the continuity of the reservoir facies and influ-
tion modeling and demonstrate how geologic data can be ence the vertical and horizontal dimensions of each cell.
synthesized to provide quantified information. The gen- Real variations in reservoir parameters may require sev-
eral application of the methodology is also described. The eral cross-sections or a three-dimensional model. Other
technology required for future studies is discussed. influences on cell dimension include computing cost, well
Reservoir simulation models are becoming increasingly sophisticated, resulting in partfrom the
development ofboth geological and engineering reservoir-description methods. The geologist's'
research usually follows a four-step procedure that is applicable to most exploration and
production operations.
DEPOSITIONAL S.P.
--.... LL-8
ENVIRONMENTS DEPTH (FT) ~
AND UNITS
z
!!to CORE LOG CAL. ~
_r= TEX.
Z.
LITH. SED. CEM.
<
.... ::I: ..0
STRUC (%)
CL. '" !::z
«
~
40
o
«
....
0 100
.... ~ :::>
LLI
Q
'"
Z a.n
« ....
....
~
z:::>
>< ~ f - -
.....
LLI
~
CL. «
~ I:Q
0 :I:
U ....
:::> !::
~
60~
.... 0 Z
Z
0 ~ :::>
~ a::J.
~
« ....~
........ c.n
1===
f-- I---
~118070 ~
W
Q vi
Z M
<t ....
....
~ z:::>
..i
- N I----
<t ....
.... «
w
* ~ 80~
-
I-
w ::I:
.... \CLAY
Q *DELTA FRONT COARSE
0 '" Z
~ c :::>
SILTSTONE AND
CL.
0..: CLAYSTONE
Fig. 4-Example core-description graph showing results of rock studies.
LOWER TRANSITION
Fig. 5-Typical photographs of the six I.ithologic units (the photo for Unit 3 is about 1 ft long),
APPROXIMATE
PILOT LOCATION
UNITS· 3-5 DISTRIBUTARY CHANNEL
CARBONACEOUS SHALE
O SILTY SHALE
F3 LAMINATED. CARBONACEOUS
Ed SANDSTONE
it.
.:..c"1~ INTERBEDDED SANDSTONE
~~.=.~:. AND SILTY CLAYSTONE
.-
10
related deposits form thinner, discontinuous sheet de-
-
·posits. The discontinuous nature of the delta-front de-
posits is best seen at the 10-acre spacing of the field.
Because the Loudon field is largely north of a sand 1.0
o
thick, only the extreme southern end of the field received
much aquifer support. Most of the field, situated in the • •
SILTSTONE AND SI TSTONE WISH
SANDSTONE, V. F. GRAINED
30
begun to increase oil recovery .
POROSITY - % BULK VOLUME
The Upper Mississippian strata in the Illinois basin Fig. 7-Porosity vs permeabil ity cross-plot for various rock types
consist of several deltaic cycles, each similar to the aids identification of reservoir rocks.
Cypress. Shallow- water carbonate deposits are inter-
spersed with the deltaic clastics. The wide areal extent
of these carbonates render them excellent units for cor-
relation at both the regional and the pilot scale. Three
deltaic cycles (Weiler, Paint Creek, and Bethel) have
been productive at Loudon.
To ascertain the three-dimensional reservoir continuity
in a sedimentary basin, field, or pilot, well logs (and
seismic sections, where appropriate) are systematically
analyzed, using a depositional model as a guide. The
objective of correlation studies is to determine the vertical
and areal limits of fundamental reservoir bodies and the
degree of contact (stacking) of the bodies in time and
space. When cores are lacking, the depositional units
often can be estimated by comparing well-log cUrves in
cored and uncored sections or wells.
o 18
The method used in the case study is shown in Fig. 9. I MILE'S'
For easier representation, the environmental data (shown
in Fig. 4) have been removed. Wells in the pilot and in the • 200-250 FT. @ 50-100 FT.
160-acre area surrounding the pilot were correlated over
the interval extending from the limestone datum to the
e 150-200 FT. 00-50 FT.
limestone marker below the Cypress formation. The sys- o 100-150 FT. ~NO SAND
The pattern of lines in the interval above and below the WELL NO.1
Weiler, together with the depositional units (see Fig. 4), 13~ 13S0
aided in determining the vertical and lateral limits of
S.P.
reservoir units in various wells. The general parallelism of DATUM:
depositional interpretation.
Shale breaks of two scales were identified in the pilot:
(1) those with radii of at least 200 ft but less than 1,000 ft -~
(see Marker 8, Fig. 9) and (2) those with radii less than.50 ~~:::;<=================~===:;::-~-=:-
~
_ _ === __ ~4~50'0'
ft (see Unit 4, Fig. 4). Tight streaks of carbonate cement
(see Well 2, Fig. 9) occurred in five wells, but were not
correlative between wells. 'WEILER BASE
T.D. 1481 ( LOGGER ) VERTI~:E\ SCALE [~O
On the 160-acre scale, sand developments similar to
NO HORIZONTAL SCALE
Unit 4 in the pilot were found to owe their general con-
T.D. 1498 ( LOGGER)
tinuity to the stacking and overlapping of individual sand Fig. 9-Example of well-log correlation in the pilot indicating detail
"blobs. " This suggested that several distributaries were required for sand continuity studies.
DEPOSITIONAL
ENVIRONMENTS KEROSENE tn DEPTH (FT) BRINE
AND UNITS POROSITY - % B. V.~O CORE LOG CAL. PERMEABILITY - MD
z 20 10~;= TEX. LITH. SED. CEM. 100 10 1.0
<C :E: ..0 1
21 .1' 14'1.7 101'.4 i ---40 STRUC (%0) i II I i
-oJ
Q,. '~
" I- 94.1......-AV. UNIT 4
et: « Z /' ~/ 0 0100 I
I-
-oJ
:E :;:) AV. AV. AV. - I
LLI
Q UNIT 4 UNIT 3 UNIT 5 \ ~_~ PERMEABILlTY I
1\ --~ CUT-OFF~
...- ·1·: - -. --
-
I_T
-;+ --,
II. t
~
- -
:\.
!.\ 60~ I -
:1,
I. •
t
,T
l
_-l -
--
.-
--~-
I -
_: __I .
et:
to-
~~- .-
• -I .. -T-.. -
1. .- - •
.-- .•
Q ~ M
« !::
ar: Z
.... ;:) CORE
/ ·5 · -I
~
~ ---
c(
....i
*~
u.I
ANALYSIS DATA ·1
~
---
80 - -
~ .-
....
-oJ
-oJ
«
w
c :r.: I- *DELTA FRONT
0 '" Z SILTSTONE AN D
ar: ci :l
Q,.
~
CLA YSTONE
Fig. ll-Example reservoir-quality profile showing porosity and permeability relationships to lithology.