Fracture Analysis in Marun Oilfield
Fracture Analysis in Marun Oilfield
net/publication/292129932
CITATIONS READS
10 6,902
3 authors, including:
All content following this page was uploaded by Ghasem Aghli on 27 January 2016.
Abstract
Fractures are considered as one of the important structures in fractured reservoirs due to their effect on fluid currents and reservoir
parameters such as porosity and permeability. Fracture parameters can only be directly calculated with core and image logs. Cores have
serious limitations, so image logs are the best method. The aim of this study is the systematic fractures analysis of the Asmari
Formation in the Marun field as one of the giant oilfields in world. The main objectives of image logs were evaluating structural dip,
characterizing natural fractures and field structure heterogeneity, and finally correlating the results with complimentary methods such
as Velocity Deviation Log (VDL), Repeat Formation Test (RFT), mud lost data, and isodip map in the carbonate Asmari Formation.
Generally, electric and ultrasonic imaging tools record vast amounts of high-resolution data. This enables geoscientists to describe in
detail the structural fracture networks. The results indicate that the highest fracture density is in the zones 1, 20, and 30 of the Asmari
reservoir that show high correlation with VDL and mud lost data. Image logs also show a range of bedding dips from 20˚ in the
northern limb to 30˚ in the southern limb with strikes ranging from 10˚ to 270˚N. Regarding the general pattern of fractures, it is
evident that they are related to the folding and are classified mainly as longitudinal, transverse, and oblique. The longitudinal pattern is
dominant and often forms open fractures. They are characterized by N50W-S50E and mainly observed in the upper Asmari zones.
Moreover, to find the vertical relation of the layer and fractures, RFT data were used. The findings revealed the presence of a vertical
relation in the upper horizons of the reservoir, especially in the eastern section due to the high fracture density.
Keywords: Asmari Reservoir, Complimentary Methods, Fracture Analysis, Image Logs, Marun Oilfield.
in the oil-based mud where oil base micro imager foreland basin (Blanc et al., 2003; Wennberg et al.,
(OBMI) and ultrasonic borehole imager (UBI) are 2007). Asmari deposition was followed by the
applied (Serra, 1989). Fractures diagnosis is not deposition of Lower Fars group evaporates known
valuable by itself because fractures study is useful as Gachsaran Formation (Aqrawi, 1993).The
when their role in the system of porosity and Asmari Formation comprises of fractured dolomites
permeability are determined (Martinez et al., 2002; and limestones and, locally, the Ahwaz sandstone
Brie, 1985). The present work is an attempt in order Member.
to determine the role of fractures and their effect on The source rock intervals were attributed to the
the reservoir properties as a vital subject in the Kazhdumi Formation, and the Gachsaran
carbonate reservoir. For this aim, VDL and RTF Formation is the cap rock. The Asmari Formation
data coupled with image log interpretations were in this field is generally subdivided into three parts.
used in our study. The lower part consists of open-marine facies
deposited mainly on a platform margin and its
Geological setting of the studied area seaward slope and basinal equivalents (50-150 m).
The giant Marun oilfield is located in the Dezful The middle and upper Asmari consists of relatively
Embayment, parallel to the general trend of the shallow-water facies (maximum 30-50 m) deposited
Zagros fold and trust belt (Fig. 1). It is one of the in the platform-top settings. The Asmari Formation
well-known oil provinces in the Middle East. The was divided into five main layers (1, 20, 30, 40, and
structural and stratigraphic features of the Dezful 50) and 5 sub-layers (11, 28, 36.30, 40.80, 50 and
Embayment have been investigated since the 60) based on alternating siliciclastic and carbonate
middle of the last century (e.g., Haynes and lithology (Gholipour, 1998). This classification is
McQuillan, 1974; Alavi, 2007). The Asmari still used for the Asmari Formation in the Iranian
Formation is the most important and proliferous National Oil Company (NIOC) and also in this
reservoir in the SW of Iran (Fig. 2). In the Dezful study.
Embayment, the formation was deposited from the
Early Oligocene (Rupelian) to the Early Miocene Methods and Materials
(Burdigalian) in a NW-SE trend basin (Stocklin, In this study, image logs as the main data and mud
1968). Continued subsidence during the deposition lost, RFT, and VDL as complimentary data are
of the Asmari Formation may be related to the available in 11 studied wells (Fig. 4). Also, an
latest stages of the closure of the Neo-Tethys, and isodip map was drawn for the Asmari Formation.
the Asmari basin may represent the early stages of a
Figure 1. Location map of the Marun Field in southwestern Iran (Red circle) and 3D seimic section for marun anticlin and other
adjacent fields provided by NISOC.
Systematic fractures analysis using image logs and complementary methods… 141
Figure 2. Cenozoic stratigraphic correlation chart of the Iranian Sector of the Zagros Basin, adopted from James and Wynd (1965).
The images are carefully described to the structural dip (Fig. 3). There are some intervals
characterize fractures and other geological features within the Asmari sequence where the image log
and finally the image results were correlated with reveals relatively thinly bedded formations.
other data. In this study, image logs were processed
and interpreted by the Geoframe 4.5 software. Complimentary methods
In this study, in order to achieve more accuracy, all
Image logs image results were correlated with complimentary
There is a general consensus that image logs are the methods including VDL and RFT. The Velocity-
best method for fractures and structural study. Deviation Log is a log calculated by combining the
Fractures are planar features with no apparent sonic log with the neutron-porosity or density log,
displacement of blocks along their planes (Fossen, which provides a tool to obtain downhole
2010). Generally, they have a steep dip in tensional information on the predominant pore type in
and wrench regimes. Their aperture may be open, carbonates. This log can also be used to trace the
tight (closed), or filled with some minerals downhole distribution of diagenetic processes and
(Schlumberger, 2005). On the image logs, fractures to estimate trends in permeability. The main
tend to occur as linear features that generally have a purpose of calculating the VDL log is creating an
dip steeper than the structural dip (Fig. 3). Open artificial log using the other petrophysical logs that
fractures, in a clay free formation, have a can detect porous types and fractured zones (Figs. 6
conductive appearance on the images due to the , 7), because thay cannot be determined by Sonic,
invasion of their aperture with the conductive Neutron, and density logs in a carbonate reservoir
drilling mud. However, the mineralized or sealed (Guadagno & Nunziata, 1993). This log is
fractures appear resistive if the filling material of calculated by first converting porosity-log data to a
their apertures is dense like calcite or anhydrite. In synthetic velocity log using a time-average
some cases, the fullest logs (PEF, RHOZ, Neutron equation. The difference between the real sonic log
and etc.) can also be very helpful for this kind of (DT) and the synthetic sonic log (DTsyn) can then
differentiation (Tokhmchi et al., 2010; Ja’fari et al., be plotted as VDL. This log shows the difference
2011). Between fracture parameters, the aperture between the real density wave values (Vpreal) and
has the greatest effect on reservoir properties artificial density wave (Vpsyn) (Anselmetti and
especially permeability (Aghli et al., 2014). The Eberli, 1999). A Repeat Formation Tester (RFT) is
main types of geological structures are bedding and a wireline tool used to sample fluid formation and
other sedimentary features. Generally, there are formation pressure. The tool is several tens of feet
strong resistivity contrasts and the sedimentary dip long with several chambers that hold several
information is of high quality where the formation gallons of fluid. The tool is lowered into the well to
is stratified. However, diagenetic processes have a producing zone. A shoe is extended from the tool
caused irregular bedding planes at few intervals and forcing it against the opposite wall of the well.
consequently hinder the precise determination of
142 Zahmatkesh et al. Geopersia, 5(2), 2015
By opening a value, fluid is allowed to flow into high confidence resistive fractures show a
the tool. As the fluid is flowing, the formation fluid relatively large scattering in their dip azimuth and
pressure is also measured. RFT data is used for strike. The most fractures were in zones 1, 20, and
revealing the presence of the vertical relation. The 30 in the Asmari reservoir (Figs. 6 and 7). In some
pressure is used to calculate formation wells, zone 1 has 300 open fractures and zone 2 has
permeability, as well. The repeat formation tester is 600 open fractures. In this field, considering the
run on an open hole. In addition, mud lost data is general pattern of fractures, it is evident that they
used as an auxiliary method because it is usable for are related to the folding and are mainly classified
all fracture studies as available data in drilled wells. as longitudinal, transverse, and oblique. The main
fracture set has the same orientation as the bedding
Results and Discussion which classifies it as a longitudinal set (Fig. 7).
After the analysis of electric and sonic image logs Strike in this set is almost N130E. The other
in the studied wells (Fig. 4), several features such fracture set is obliquely oriented relative to the
as fractures, beddings, induced fractures, borehole bedding and is considered as an oblique set. Strikes
breakouts, and stylolites were identified. Because in this set are N30E to N60E which can mostly be
of the importance of open fractures and bedding in seen in the eastern section (Fig. 4). It is notable that
this study (based on this study purpose), they were the high density fractures in the central and eastern
determined on image logs and fracture density was sections (well A and B, Fig. 4) are related to the
calculated for all wells. The results of each step are axis re-circulation by the basement fault which is
presented below: probably in the south section.
images and correspond to the surface or boundaries Confidence Bedding Dips (HC), and those from the
separating two beds of different lithology (Serra, uneven and vague boundaries are classified as Low
1989). So the dips computed from the first type of Confidence Bedding Dips (LC).
layer/bed boundaries are classified into High
Figure 4. Represent of fracture and bedding strike in studied field using image logs results. As is evident, there are two fractures set in
the Eastern section.
Figure 5. Represent of stereonet for fractures (blue points) and bedding (green pints) on Iso-mud loss map in the studied field. (Mud
loss data represents center and eastern as fractured section in the studied field which show high correlation with images data)
144 Zahmatkesh et al. Geopersia, 5(2), 2015
Figure 6. Results for Image log in the well A and its comparison with VDL log (Track 4)
Both classes of bedding dips compute nearly the Complimentary Methods Results
same average and dominant dip magnitude and dip Velocity Deviation Log
azimuth. Therefore, based on both types of bedding In this study, for evaluating the dominant pores
dips, a structural dip of 6 degrees N50W-S50E pattern and fractures detection, VDL log was drawn
strike can be computed for the studied wells (Fig. using Sonic and Neutron digital data for the Asmari
5). In this field, axis re-circulation has caused the reservoir in well A and B (Fig 4). Generally, there
bedding shift in the total area especially in the is an inverse relationship between porosity and
eastern section which in all reservoir parts fractures velocity, so velocity decreases with increasing
intercept the bedding. porosity (Wyllie et al., 1956). Three zones are
Systematic fractures analysis using image logs and complementary methods… 145
Figure 7. Results for Image log in the well B and its comparison with VDL log (Track 4)
146 Zahmatkesh et al. Geopersia, 5(2), 2015
Figure 8. Repeat formation test (RFT) result for well A in studied field. Result show vertical relation in the upper horizons of the
reservoir. Red circles indicate no conductivity between fractures due the presence of shale.
that have high correlation with RFT and image isodip map represents the changes of bedding dip
results (Fig. 5). and bending layers which can identify areas with
high potential for fractures (Roberts, 2001). The
Isodip Map isodip map for the Asmari Formation in the studied
In this study, an isodip map is generated from the well represents the southern edge as an area with
second derivative underground contour (UGC) to high slope and curvature (Fig. 10) which has a high
provide an overview on the fields’ structure and the potential for fractures; this confirms the results of
changes of axial bending. The first derivative of the this study.
Figure 9. Display of RFT data for a well in western section. Results indicate low vertical relation due the low fracture density in this
section.
Based on Figure 12, the stress direction in the Zagros general trend. This result confirms the
eastern section of this oil field is therefore different hypothesis of the presence of a basement fault in
from the western section and does not follow the the southern edge.
Figure 10. Isodip map for Asmari formation in the studied well
Figure11. Maximum and minimum main stresses direction in the studied field around wells (blue arrows indicate max stress and green
arrows indicate min stress) using the borehole breakout. Stress directions different from eastern section. Therefore, eastern section does
not follow Zagros trend.
Systematic fractures analysis using image logs and complementary methods… 149
The second question is about the effect of S60E and are mainly observed in the upper Asmari
fractures on the reservoir properties. What is the zones. Images indicate that the highest fracture
role of fractures in the reservoir? For resolving this density is in the zones 1, 20, and 30 of the Asmari
problem, VDL, RTF data, and other complimentary reservoir that show high correlation with other data.
methods such as mud lost data were used. Based on RFT data reveal that reservoir zones have low
VDL, we found that fractures are the most vertical relation in the western section of the field
important features in the reservoir especially in the due to the low fracture density. However, there is a
upper Asmari zones. On the other hand, production relatively high vertical relation in the central and
in this field is a combination of fractures and rock eastern sections owing to the high fracture density.
matrix. Also, RTF data demonstrate that in the VDL as a method for the detection of fracture zones
western section of the studied field, reservoir zones shows high correlation with fracture density log
have low vertical relation due to the low fracture (FVDC) and indicates that the production in the
density but in the central and eastern sections, Asmari reservoir of this field can be the
zones have high vertical relation because there is combination of fractures and rock matrix. Image
high fracture density. logs also show a range of bedding dips from 20˚ in
the northern limb to 30˚ in the southern limb with
Conclusion strikes ranging from 10˚ to 270˚N. Overall, fracture
The main aim of this study was the systematic sets are ambiguous in this study because fracture
fractures analysis and the determination of their sets in the eastern section were different from the
role in the reservoir properties of a carbonate western ones. In situ stress direction for the studied
Asmari Formation in Marun oilfield. Then, a well indicates that the stress direction in the eastern
comparison between image logs results and section of this oilfield is different from the western
complimentary methods such as RFT, VDL, mud section and does not follow the Zagros general
lost data, and isodip map has been drawn for more trend. This point can confirm the hypothesis of the
accuracy. Image logs interpretation shows that presence of a basement fault in the southern edge
fracturing is present at a number of places. There which could be created by other dominating
are more than 600 open fractures in some zones fracture set. In this field, other auxiliary methods
which represent the Asmari Formation as a such as mud lost data and isodip map, in
fractured reservoir. Between fracture sets, accordance with the other methods, demonstrate
longitudinal patterns are dominant and often form that the central and eastern parts are fractured.
open fractures. They are characterized by N60W-
References
Aghli, G., 2013. Fracture analysis of the Asmari reservoir in Balarud oilfield using the image logs: M.Sc dissertation,
Shahid Chamran University, Ahvaz, Iran.
Aghli, G., H. Fardin, R. Mohamadian, G. Saedi, 2014. Structural and Fracture analysis using EMI and FMI image Log in
the Carbonate Asmari Reservoir ( Oligo-Miocene ), SW Iran. 4( 2): 43–58.
Alavi, M., 2004, Regional stratigraphy of the Zagros fold-thrust belt of Iran and its proforeland evolution: American
Journal of Science. 304 (1): 1–20
Alavi, M., 2007, Structures of the Zagros fold-thrust belt in Iran: American Journal of science, 307 (9):1064–1095.
Anselmetti, F. S., and G. P. Eberli, 1999, The velocity-deviation log: a tool to predict pore type and permeability trends in
carbonate drill holes from sonic and porosity or density logs: AAPG bulletin. 83(3): 450- 466.
Aqrawi, A. A. M., 1993, Miocene evaporitic sequence of the southern Mesopotamian Basin: Marine and petroleum
geology. 10 (2): 172–179.
Blanc, E.-P., M. B. Allen, S. Inger, and H. Hassani, 2003, Structural styles in the Zagros simple folded zone, Iran: Journal
of the Geological Society.160 ( 3)401–412.
Brie, A., Johnson, D. L., Nurmi, R. D., 1985. Effect of spherical pores on sonic and resistivity measurements, in SPWLA
26th Annual Logging Symposium: Society of Petrophysicists and Well-Log Analysts.
Darling, T., 2005. Well Logging and Formation Evaluation: Elsevier Science, Gulf drilling guides.
Fossen, H., 2010. Structural Geology: Cambridge, Cambridge University Press.
Gardner, G. H. F., Gardner, L. W., Gregory, A. R., 1974. Formation velocity and density-the diagnostic basics for
stratigraphic traps: Geophysics. 39 (6): 770–780.
150 Zahmatkesh et al. Geopersia, 5(2), 2015
Gholipour, A. M., 1998, Patterns and structural positions of productive fractures in the Asmari Reservoirs, Southwest Iran:
Journal of Canadian Petroleum Technology. 37( 1): 44-50.
Guadagno, F. M., Nunziata, C., 1993. Seismic velocities of fractured carbonate rocks (southern Apennines, Italy):
Geophysical Journal International.113 (3): 739- 746.
Halliburton, 1996. Electrical Micro Imaging Service (Sales Kit): Halliburton.
Haynes, S. J., McQuillan, H., 1974. Evolution of the Zagros suture zone, southern Iran: Geological Society of America
Bulletin. 85 (5): 739–744.
Ja’fari, A., A. Kadkhodaie-Ilkhchi, Y. Sharghi, and K. Ghanavati, 2011, Fracture density estimation from petrophysical log
data using the adaptive neuro-fuzzy inference system: Journal of Geophysics and Engineering. 9 (1): 105–114.
Khoshbakht, F., Azizzadeh, M., Memarian, H., Nourozi, G. H., Moallemi, S. a., 2012: Comparison of electrical image log
with core in a fractured carbonate reservoir: Journal of Petroleum Science and Engineering. 86-87: 289–296.
Khoshbakht, F., Memarian, H., Mohammadnia, M., 2009. Comparison of Asmari, Pabdeh and Gurpi formations' fractures,
derived from image log: Journal of Petroleum Science and Engineering. 67 (1-2): 65–74,
Martinez, L. P., Hughes, R. G., Wiggins, M. L., 2002. Identification and Characterization of Naturally Fractured
Reservoirs Using Conventional Well Logs.
Mohebbi, A., M. Haghighi, and M. Sahimi, 2007, Conventional Logs for Fracture Detection & Characterization in One of
the Iranian Fields, in International Petroleum Technology Conference: International Petroleum Technology Conference.
Nelson, R., 2001, Geologic analysis of naturally fractured reservoirs: Gulf Professional Publishing.
Nie, X., Zou, C., Pan, L., Huang, Z., Liu, D., 2013. Fracture analysis and determination of in-situ stress direction from
resistivity and acoustic image logs and core data in the Wenchuan Earthquake Fault Scientific Drilling Borehole-2 (50–
1370m): Tectonophysics. 593:161–171
Rajabi, M., Sherkati, S., Bohloli, B., Tingay, M., 2010. Subsurface fracture analysis and determination of in-situ stress
direction using FMI logs: An example from the Santonian carbonates (Ilam Formation) in the Abadan Plain, Iran:
Tectonophysics. 492 (1-4): 192–200.
Rezaee, M. R., Chehrazi, A., 2005. Fundamentals of Well Log Interpretation: Tehran, University of Tehran (in Persian).
Roberts, A., 2001, Curvature attributes and their application to 3D interpreted horizons. 19: 85-100.
Saedi, G., 2010, Fracture analysis of Asmari reservoir in Lali oilfield using the FMI image log: Shahid Chamran
University, Ahvaz, Iran.
Schlumberger, 1994, FMI Fullbore Formation MicroImager: Houston, Schlumberger Educational Services.
Schlumberger, 2005. GeoFrame 4.2, BorView User’s Guide: Schlumberger Ltd.
Schlumberger, 2003. Using borehole imagery to reveal key reservoir features, in Reservoir Optimization Conference.
Serra, O., 1989. Formation MicroScanner image interpretation: Schlumberger Educational Services.
Serra, O., Serra, L., 2004.Well logging: data acquisition and applications.
Stocklin, J., 1968. Structural history and tectonics of Iran: a review: AAPG Bulletin. 52 (7): 1229–1258.
Thompson, L. B., 2000. Fractured reservoirs: Integration is the key to optimization: Journal of petroleum technology. 52
(2): 52–54.
Tingay, M., Reinecker, J. , Müller, B., 2008. Borehole breakout and drilling-induced fracture analysis from image logs:
World Stress Map Project. p. 1–8.
Tokhmchi, B., Memarian, H., Rezaee, M. R., 2010. Estimation of the fracture density in fractured zones using
petrophysical logs: Journal of Petroleum Science and Engineering. 72 (1): 206- 213.
Wennberg, O. P., Azizzadeh, M., Aqrawi, A. A. M., Blanc, E., Brockbank, P., Lyslo, K. B., Pickard, N., Salem, L. D.,
Svana, T., 2007. The Khaviz Anticline: an outcrop analogue to giant fractured Asmari Formation reservoirs in SW Iran.
Wyllie, M. R. J., Gregory, A. R., Gardner L. W., 1956. Elastic wave velocities in heterogeneous and porous media:
Geophysics. 21 (1): 41–70.