SPE 59464 Full Field Model Study of A Dense Highly Fractured Carbonate Reservoir
SPE 59464 Full Field Model Study of A Dense Highly Fractured Carbonate Reservoir
SPE 59464 Full Field Model Study of A Dense Highly Fractured Carbonate Reservoir
Abstract
The purpose of this paper is to present the results of a
comprehensive subsurface study of a field in South of Iran.
The reservoir is typical of several carbonate fields in south
west Iran, comprising low porosity carbonate , but highly
fractured. After several years of natural depletion a gas
injection was started at a time that reservoir became saturated
at the crest to prevent further shrinkage of the oil and keep the
reservoir oil pressure constant. The source of gas was a deeper
reservoir where the gas was produced and injected by its own
pressure through one producer and one injector.
A three dimensional model was set up using a 3D static
model as the basis for geological input. All available data
including, fluid losses during drilling, well test and flowmeter
surveys were used to define the fracture properties. A good
history match was obtained for pressure and fluid levels with a
commercial black oil simulator using the model in dual
porosity mode. The calibrated model was then used in
prediction mode to forecast the reservoir behavior under gas
injection pressure maintenance and also repressuring close to
original crestal pressure. The simulation results show large
increase in ultimate recovery due to gas injection as compared
with natural depletion.
Introduction
This field is an elongated antyclinal structure about 20 km.
long and 5 to 7 km. wide. The oil production horizon is
Asmari limestone belonging to Oligocene-Miocene age having
a low matrix porosity with net average of about 9% and less
than 1 md average permeability.However extensive fracturing
particularly on the crestal areas of the field has increased the
effective rock permeability to several darcy so that
productivity index of more than thousand bbl/d/psi has been
Fracture Modeling
High well productivities and fluid losses during drilling
suggested existence of a well connected network of open
fractures on major part of the structure. Fracture occurrence is
controlled by matrix properties and by areal variations in
intensity of the fracturing process. As an example of degree of
fracturing well AS-08 has to be completed after 100 ft of
drilling into Asmari formation due to complete losses. The
section drilled is tight carbonate with low porosity. In this
study there has been limited core data on three wells and no
image log. However fluid loss data, flowmeter surveys and
build up tests have been used to define the fracture intensity
and its areal distribution.
The analysis of different data has shown that the fracturing
in this field is mainly fold related and wells with very high
productivities are located in the hinge area (highest curvature)
of the structure. Therefore for a qualitative estimate of
fracturing curvature analysis of the top structure map was
considered to define fractured area in different zones. Other
information like fluid loss, well productivity and flowmeter
surveys were incorporated for further refinement of the
fractured area as well as estimating fracture properties such as
aperture, size (block height), porosity and permeability.
Dynamic Modeling
The dynamic reservoir simulation was initialized by exporting
the 3D static geological model into a commercial black oil
reservoir simulator. The fluid properties of the Asmari which
was constant with depth except for the effect of temperature
was defined at the crest of the structure as the reservoir was
initially undersaturated by more than 550 psi. The oil which is
light (35 API) has had a bubble point pressure of 2260 psia,
oil formation factor of 1.35 bbl/stb and gas oil ratio of 700
scf/stb at initial condition.
The rock properties were taken from the available core
data in this field. However where the data was scarce
information from similar Asmari reservoirs has to be
substituted to complete all the necessary data needed for 3D
simulation. The relative permeability and capillary pressure
data were taken from other reservoirs with adjustment of the
end points to define 8 different rock types having either a
porosity range and/or water saturation range. The character of
these curves leans towards mixed to oil wet conditions. The
porosity-permeability relation developed for matrix rock
showed a permeability of 1 md at 11% porosity.
The early 3D simulation that contained the geological and
petrophysical details of all layers, including structure,
distribution of reservoir properties as well as production, PVT
and relative permeability data was initialized with a Cartesian
grid system. However, to reduce the number of grid blocks a
non-orthogonal grid system was developed with the
dimensions of 221415 or more than 9200 cells in dual
porosity mode. The fracture porosities were based on the best
estimate from the geological study and also experience from
other fields in the area. The fracture porosity was on average
about 0.1 percent of the bulk volume which in this low
SPE 59464
SPE 59464
criteria has to be specified. The constraints for producing gasoil ratio and water cut and also well tubing pressure were
selected based on field operation facilities limitations. As there
is no desalination facility in the field the constraint on watercut was very tight.
In order to select best areas in the field for infill well
drilling the net oil column map at the end of history was
generated. The infill wells were defined in areas of highest
remaining oil column. A constant rate of 25000 stb/d which
was close to the final rate of production during history was
specified and the model was run in prediction mode. Rate cut
back, shutting off the most offending layer producing at high
gas-oil ratio or water cut and opening up the new infill wells
to meet the production level as long as possible and then
declining the production to the predefined abandonment rate
of 3000 stb/d were all performed automatically by the model.
The gas injection to keep the oil pressure constant was
defined by injecting gas 10 percent more than the voidage
replacement. With this volume of gas the oil pressure
increased by about 140 psi during 16 years of forecast period
and the gas pressure at the crest was very close to the original
crestal pressure.
The cumulative oil recovery from the field down to 2016
reached 341 million stock tank barrel which is equivalent to
24 percent recovery factor. This is about 5 percent more than
recovery factor obtained under natural depletion. The gas
injection requirement for the above gas injection repressuring
will be about 380 billion ft 3 which can be supplemented by
deeper Bangestan reservoir and reinjection of the produced
gas. Production, injection and pressure behavior of the
reservoir during prediction period have been shown in Figs. 10
to 12.
Conclusions
The integrated geological-engineering study has resulted in
better understanding of in place volumes, aquifer strength,
rock and fluid characteristic and behavior of This field.
Through this study the superiority of gas injection
repressuring close to original crestal pressure was verified and
due to the availability of gas can be implemented without any
further investment as the gas injection by natural flow and
reinjection of produced gas are already in place and running.
The future infill well locations have been optimized based on
maximum remaining oil column throughout the field that can
be used for future development the field that has 40 percent of
its remaining oil reserves to be produced from now onward.
Acknowledgments
We would like to thank National Iranian Oil Company-South
Oil Fields management for their permission of publishing this
paper.
200000
2600
180000
Gas pressure
160000
2500
140000
2400
120000
100000
2300
80000
2200
60000
40000
2100
20000
0
2000
0
10
15
20
25
30
Year
10
SPE 59464
3100
25000
2900
Model
Actual
Pressure
20000
2700
15000
2500
10000
2300
Rate
5000
2100
1900
0
0
10
15
20
25
30
Year
3100
2600
100000
2400
50000
2200
2000
0
10
15
20
25
30
Year
Fig.3
History of oil production and individual wells oil
pressures as compared with the model average oil pressure
potential
2800
150000
2900
oil pressure psi
200000
40000
Model
Actual
Pressure
35000
30000
2700
25000
2500
20000
15000
2300
Rate
10000
2100
5000
1900
0
0
10
15
20
25
Year
30
3200
250000
3100
35000
30000
2700
25000
2500
20000
Rate
15000
2300
10000
2100
Pressure
350
Oil rate
2900
Oil pressure psi
30000
40000
Model
Actual
25000
325
20000
15000
275
10000
250
5000
225
Water rate
0
0
10
15
20
25
300
Cumulative oil
5000
1900
SPE 59464
200
25
30
30
35
40
45
Year
Year
Model
Actual
4500
Oil pressure potential psi
5000
5500
Water oil level
6000
1400
2950
1200
2900
1000
GOR
2850
800
2800
600
2750
Oil pressure
2700
400
200
2650
2600
6500
0
10
15
20
25
0
25
30
30
1000
30000
GOR
Productio rate stb/d
800
700
20000
600
Model
500
Actual
400
10000
300
Rate
5000
900
25000
200
100
0
5
10
15
45
20
25
30
Year
Fig.8 History of gas oil and water oil levels and comparison
with model calculated values
40
Year
Year
15000
35
during
400000
3400
350000
3200
300000
3000
250000
Gas pressure
2800
200000
2600
150000
pressure psi
Levels ft.s.s.
3000
4000
2400
100000
Gas injection rate
2200
50000
0
2000
25
30
35
40
45
Year