Sub-Surface Master Development Plan-Phase 12
Sub-Surface Master Development Plan-Phase 12
Oil&GasEngineeringManagement
PHASE 12
SUB-SURFACE MASTER DEVELOPMENT PLAN
JUNE2008
TABLE OF CONTENTS
1
1.1
1.2
Geology, Geophysics and Reservoir................................................................................ 1
1.2.1 Structure ............................................................................................................................. 2
1.2.2 Reservoir Characteristics ............................................................................................... 2
1.2.3 Hydrocarbons in Place .................................................................................................... 3
1.2.4 Reservoir Fluids ................................................................................................................ 3
1.3
Reservoir Development Plan .............................................................................................. 4
1.3.1 Reservoir Simulation ....................................................................................................... 4
1.3.2 Plateau Length .................................................................................................................. 4
1.4
1.5
2.1
Geological History................................................................................................................. 5
2.2
2.3
2.4
Top-Seal ................................................................................................................................... 9
2.5
2.6
2.7
2.8
Scope of Work for Geophysical Interpretation for Phase 12 of South Pars ........ 12
2.8.1 3D Seismic Interpretation ............................................................................................. 12
2.8.2 3D Seismic Inversion ..................................................................................................... 13
2.8.3 Quantitative Reservoir Characterization Based on Seismic Attribute.............. 13
2.8.4 Interim Report .................................................................................................................. 14
2.8.5 Deliverables...................................................................................................................... 14
2.9
Litho-stratigraphy ............................................................................................................... 14
2.10
2.11
Sedimentary Environment................................................................................................. 17
2.12
Zonation ................................................................................................................................. 19
3.1
Methodology and Procedures of Interpretation .......................................................... 25
3.1.1 Multimin Model ................................................................................................................ 25
3.1.2 Lithological Parameter Picking ................................................................................... 26
3.1.3 Petrophysical Parameters ............................................................................................ 26
3.1.4 Selection of Water Saturation Model ......................................................................... 28
3.1.5 Clay Volume Calculations ............................................................................................ 30
3.1.6 Gas Correction ................................................................................................................ 31
3.2
Petrophysical Interpretation ............................................................................................. 31
3.2.1 Depth Matching ............................................................................................................... 31
3.2.2 Precalc ............................................................................................................................... 32
3.2.3 Environmental Correction ............................................................................................ 32
3.2.4 Preliminary Multimin Modeling ................................................................................... 33
3.2.5 Evaluation of Vertical Wells ......................................................................................... 33
3.2.6 Evaluation of Deviated Wells ....................................................................................... 33
3.3
3.4
Scope of Work for Petrophysical Logging for Phase 12 of South Pars ................ 35
4.1
4.2
4.3
4.4
Volumetric ............................................................................................................................. 40
4.5
4.6
5
5.1
5.2
Data Availability ................................................................................................................... 42
5.2.1 Core Data .......................................................................................................................... 43
5.2.2 Well Test Data .................................................................................................................. 43
5.2.3 Reservoir Fluid Data ...................................................................................................... 43
5.3
Reservoir Dynamic Modeling ........................................................................................... 44
5.3.1 Dynamic Model Grid....................................................................................................... 44
5.3.2 Full Field Fluid Model .................................................................................................... 44
5.4
5.5
5.6
5.7
5.8
Well Characteristics ............................................................................................................ 48
5.8.1 Well Perforation .............................................................................................................. 49
5.8.2 Well Simulation ............................................................................................................... 49
5.8.3 Well Constraints.............................................................................................................. 49
5.8.4 Reservoir Simulation Results ...................................................................................... 50
6.1
Introduction........................................................................................................................... 50
6.2
6.3
6.4
6.5
Technical Well Specification ............................................................................................ 53
6.5.1 Well Trajectories ............................................................................................................. 53
6.5.2 Casing Design ................................................................................................................. 54
6.5.2.1 32" Drilling Phase 26" Conductor Pipe ........................................................ 54
6.5.2.1 24" Drilling Phase 18" Surface Casing ..................................................... 54
6.5.2.2 16" Drilling Phase 13" Intermediate Casing ............................................ 55
6.5.2.3 12" Drilling Phase 9" Production Casing .............................................. 55
6.5.2.4 8" Drilling Phase 7" Production Liner ....................................................... 55
6.5.3 Cementing Design .......................................................................................................... 55
6.5.3.1 The 26" CP .............................................................................................................. 56
6.5.3.2 The 18" Casing ................................................................................................... 56
6.5.3.3 The 13" Intermediate Casing........................................................................... 56
6.5.3.4 The 9"-10" Intermediate Casing .................................................................. 56
6.5.3.5 The 7" Production Liner ...................................................................................... 56
6.5.4 Completion Design ......................................................................................................... 57
6.5.4.1 Completion Accessories ..................................................................................... 57
6.5.4.2 Wellhead and Xmas Tree..................................................................................... 57
6.5.5 Drilling and Completion Fluid...................................................................................... 58
6.5.5.1 32" Hole Section .................................................................................................... 58
6.5.5.2 24" Hole Section .................................................................................................... 58
6.5.5.3 16" Hole Section .................................................................................................... 58
6.5.5.4 12" Hole Section................................................................................................. 59
6.5.5.5 8" Hole Section ................................................................................................... 59
6.5.5.6 Completion Fluid ................................................................................................... 59
6.5.6 Solid Control Management ........................................................................................... 59
6.5.7 Perforation and Stimulation Strategy ........................................................................ 60
6.5.8 Well Clean-up Strategy .................................................................................................. 60
Appendixes
Appendix A
Appendix B
Appendix C
Appendix D
Tables
Figures
Well Trajectories
Well Programs
1 Executive Summary
1.1
Introduction
The South Pars Gas field is located in the Persian Gulf, at some 100 KM from
shore. Four huge condensate rich gas bearing reservoirs have been identified
(K1, K2, K3 and K4) in this field.
NIOC has envisaged developing this field by phases of some 1000 MMSCFD
wet gas production each. TOTAL / GasProm / Petronas have already
developed two such phases (referred to as "Phase 2" and "Phase3"). NIOC
developed South Pars development "Phase 1".
ENI / PETROPARS consortium developed Phases 4 & 5 and PETROPARS /
STATOIL joint ventures developed phases 6 / 7 /8. Phases 9 & 10 are
developing by NIOC / LG In each contract, respecting phases are integrated
in order to share common facilities and thus reduce the overall cost of the
development.
NIOC is pursuing further development under Phase 11 and "Phases 13 &14"
and "Phases 15 and 16", and "Phases 17 and 18", "Phases 19-21" and
"Phases 22-24" developments which are at tendering stage. It is believed that
developing the whole field will lead to 24 development phases.
This Master Development Plan provides the basic support data involving
geophysical, geological, petrophysical, reservoir and drilling / completion
information for the development of Phase 12 of South Pars development
project.
1.2
The current state of knowledge on the South Pars field relies on seismic
surveys shot between 1988 and 2003 and 19 exploration / appraisal wells
drilled from 1991 until 2007 together with the information gained through
drilling of development wells in Phases 1 (12 wells), Phase 2 & 3 (20 wells),
Page1
Phase 4 & 5 (24 wells), Phase 6 & 7 & 8 (30 wells) and Phases 9 & 10 (24
wells which are drilling). With the data acquired, the Kangan-Upper Dalan (K1
to K4) reservoirs can be described as structurally simple reservoirs on which
some geological and reservoir unknowns still exist such as the extension of
the reservoir, hydrocarbon source, gas-water level, H2S content in reservoir
fluid at different locations at different reservoir depth levels, etc.
1.2.1
Structure
1.2.2
Reservoir Characteristics
Page2
1.2.3
Hydrocarbons in Place
Zone
K1
K2
K3
K4
TOTAL
Current GIIP(TCF)
Current
OOIP (MMSTB)
109.82
74.15
117.96
175.83
477.76
3151
3137
4991
7740
19019
1.2.4
Reservoir Fluids
The initial condensate yield for K2, K3 and K4 from what has been produced
over the past couple of years after production start-up of Phases 2 & 3 and 4
& 5 is more than 42SBBL condensate per MMSCF of gas (wellhead
equivalent). It is expected that by introduction of leaner K1 (CGR for this layer
Page3
1.3
The Phase 12 development of the South Pars field includes all four reservoir
units K1 to K4. The development is based on an average plateau gas
production rate of 3000MMSCFD.
1.3.1
Reservoir Simulation
1.3.2
Plateau Length
1.4
Reservoir Uncertainties
Page4
A number of uncertainties are recognized which have an impact on gas-inplace estimate. These uncertainties are as follows, in decreasing order of
importance:
structural definition on the flanks of the field, position and configuration of the
gas-water-contact, actual distribution of the reservoir properties over the field
(rock and fluid), H2S content of reservoir fluid at different locations and depth
levels and its origin, faulting and possible effect on reservoir segmentation.
1.5
The completion studies have shown the interest reducing the pressure losses
in the tubing. A 7" monobore completion is recommended. The following
casing program is only tentative and could be modified depending upon the
well deviation: 26" CP x 18 5.8" x 13 3.8:" x 9.5/8" -10 3/4" x 7" monobore
completion. The casing program and setting depths for reaching the K4 or the
Kangan-Dalan K2 / K3 seem to be identical. Corrosion resistant Alloy (CRA)
should be chosen for the liner lap and the tubing in order to protect against
corrosion during production as H2S and CO2 is expected in the well stream.
On the completion side, the liners will be cemented and perforated and
cleaned. The clean up of the wells will be further completed through matrix
acidizing right after wells are completed at surface. It is recommended to run
matrix acid stimulations using 28% HCl acid and diverting agents through bull
heading.
Geological History
(120KM * 80KM). South Pars-North Dome Gas field is part of the huge NNESSW trending Qatar-Fars Arch which is one of the major structural elements
in central Persian Gulf area (figure 1).
The main structure of the South Pars / North Field at the Kangan and Dalan
level is an anticline fold. Its axis is elongated and curved from the South (In
Qatar) to the North East and dips 0.5 to the North East. The flanks are
approximately symmetrical, dipping 0.5 to the South East, and 0.6 to the
North West. The minimum Depth lies within the Qatar Portion of the field,
some 25KM from the Iran / Qatar border.
The Qatari part was discovered by Shell in 1971 by the North West Dome-1
(ND-1) well and the South Pars Field was subsequently discovered by NIOC
in 1991 by drilling South Pars-1 well (SP-01) that encountered gas reservoir in
the Permo-Triassic Kangan-Dalan carbonate formations. In 1991-92 the field
was confirmed by drilling three additional wells SP-02, SP-03 and SP-04 in
the original phases (The development of South Pars is being implemented by
subdividing the field into several development areas or phases); further fifteen
appraisal / delineation wells (SP-5 to SP-15, SPD2-04, SPD4-01, SPD10-08
and SPD11-08) which some of them are now as producing wells were drilled
in 1998-2007 on the structure. Also seventy eight development wells of which
eleven wells in phase 1, nine wells in phase 2, nine wells in phase 3, twenty
two wells in phases 4 & 5 and twenty seven wells in phases 6 & 7 & 8 have
been drilled in South Pars area. It must be mentioned that the drilling of
development wells of phases 9 & 10 has been started since the beginning of
2007.
The reservoir consists of about 400 meters of limestone and dolomites with
subordinate anhydrite, divided into the Dalan and Kangan formations
(together comprise the Khuff Formation of Arabian terminology). The latter
formations are separated by Permian-Triassic boundary and are divided into
five main units. The Kangan Formation is divided into K1 and K2 units and
Dalan Formation into K3, K4 and K5 units. K3 and K4 units are called as
Upper Dalan and K5 unit is called as Lower Dalan which is separated from K4
Page6
by a massive anhydrite body known as Nar member. This unit exhibit poor
reservoir characteristic (figure 2). The Dalan and Kangan succession has
been subdivided into significant large scale depositional sequences
representing major 3rd order cycles consisting KS1 to KS4 which are bounded
by distinct stratigraphic surfaces. In general sequence stratigraphic framework
coincides with main reservoir units which are lithological based.
2.2
Regional Geology
South Pars field is located on the Qatar-Fars Arch, one of the major structural
elements of the Central Persian Gulf Area. Persian Gulf is part of what is
referred to in the Plate Tectonic literatures as the Arabian Plate and / or the
Middle East Sedimentary Basin which is approximately 3000KM in length
and 2000KM in width. It is bounded from the north by Turkey Bitlis Suture,
from the west by the Red Sea, from the east by Zagros Mountain (Zagros
Thrust), and from the south by the Gulf of Oman and Owen Fracture Zone in
the Arabian Sea (figure 1).
The major structural features of the area are the results of two major tectonic
phases:
The first tectonic phase is the Amar Collision which was occurred 620-
640 million years ago along the north-south trending Amar weak zone in the
Arabia Shield.
The second tectonic phase corresponds to the Najd Rift System with
Page7
years ago. These structures are formed due to the collision along the northerly
trending Amar Suture in the central part of the Arabian Shield. The major
north-south trending features are Qatar / Fars Arch, En Nala Anticline,
Khurias-Burgan Anticline, and Maaqala Anticline. These features have shaped
most of the major oil and gas field structures in Oman, Saudi Arabia, Kuwait,
southern Iraq, and southwestern Iran and also the Persian Gulf area. The
super giant structures of South Pars / North Dome, Ghawar, Safaniya-Khafji,
and Burgan owe their formation to these north-south structural features.
This region holds over two third of the oil and one third of the gas reserves
discovered in the world. This is due to a combination of some favorable
factors in this region. The most important factor is the prolonged and almost
uninterrupted
sedimentation
history
since
the
Precambrian
time.
and
southern
basins
with
somewhat
different
sedimentary
2.3
Among the exploration / delineation and development wells drilled in this area
to date, only in one well, SP-01, a short penetration of 19 meters was made
into the Faraghan Formation. All the other wells drilled in the South Pars have
been terminated in the upper parts of the Nar Member, as the objective of the
drilling was gas production, reservoir definition and / or fluid contact
identification in the Kangan-Upper Dalan reservoirs. This provided reasonable
stratigraphic knowledge down to the top of the Nar Member of the Dalan
Formation. But, stratigraphic information of the deeper horizons has to come
from the nearby areas of the sedimentary basin. The expected regional
stratigraphic sequence is shown in figure 4.
2.4
Top-Seal
The top seal of the Kangan and Dalan gas reservoirs is provided by the
Triassic Dashtak Formation and especially its Aghar Shale Member (Sudair
Shale) which directly overlies the Kangan formation. This shale is reddishbrown shale and believed to have a terrestrial origin. Additional sealing
capacity is provided by numerous anhydrites and tight carbonates in Dashtak
Formation.
2.5
Source Rocks
The possible source rocks for the South Pars-North Field accumulation are:
Based on the recent studies it is accepted here that the most likely source are
Silurian shales of Sarchahn (Qusaiba of Saudi Arabia) Formation.
Page9
2.6
The South Pars Gas Field reservoir is shared with Qatar in median line and
the portion, falling within Iranian water, was detected through the seismic
survey shot in 1989 and 1992. Since then, eighteen exploration and appraisal
/ delineation wells (SP-1 to SP-15, SPD10-08, SPD11-08, SPD2-04 and
SPD4-01) have been drilled. The top of reservoir roughly sets around 2700M
from the sea level and the wet gas (gas and condensate) is in a gross rock
thickness of about 400M immediately below the reservoir top.
K1 mainly consisted of dolomite, Packstone-Wackstone limestone and streaks
of shale. There are also three rather thin anhydrite sub layers at the middle
and in the lower half of the layer. K2 mainly consisted of Wackstone,
Packstone and Grainstone limestone interbedded with thick dolomite layers.
K3 mostly consisted of Mudstone, Grainstone limestone interbedded with thick
dolomite layers significant anhydrite interval occurs at the bottom part of K3.
K4 mostly consisted of Packstone, Oolitic Grainstone limestone, dolomite and
rare thin layers of anhydrite.
K1 thickness is about 110M. The rock quality is poor and its average porosity
and water saturation is about 12 and 18 percent. K2 thickness is about 43M.
and its average porosity and water saturation is about 11 and 7 percent. K3
thickness is about 120M. and its average porosity and water saturation is
about 10 and 8 percent. K4 thickness is about 145M. and its average porosity
and water saturation is about 16 and 7 percent respectively.
The difference in GWC at different well locations indicates a tilted gas-watercontact. Based on the available information, there is a tilt in GWC from westsouth-west to east-north-east.
2.7
National Iranian Oil Company in 1992-1993. This survey covers northern part
of South Pars Gas Field.
9
National Iranian Oil Company in 1994. This data acquisition covers an area of
85 square kilometers including the area between wells SP-01 and SP-02. This
survey was carried out to improve the structural definition of the field structure
for Phase 1 development area.
9
ASA Company (GGS) for AGIP over an area of 1015 square kilometers,
covers phases 4, 5, 6, 7 and 8 development areas.
9
The third and the most recent 3D seismic survey is carried out by GGS
for POGC over an area of 2665.90 square kilometers that is under processing
by SPECTRUM.
Main structural results found from seismic interpretation are:
dome which has a more gentle dip on the southeast flank (0.5) as compared
with the northwest flank (1).
Page11
Kangan / U.Dalan reservoir has shown a lot of variation from east to west due
to the geological lateral variations.
There are some small graben-like features to the east of Phase 2 & 3,
direction along the edge of the crestal area towards the 3H-1 structure.
Some faults with ENE-WSW orientations are observed to the north of
Phase 2 & 3, at the edge of the crestal area and at the east of the Phases 12
and 8. They are linked to a series of downward faulting to the north that are
observed above the Cretaceous level.
2.8
There is only one well (SP-13) in the study area with full set logs and VSP
data, but consultant should use well information of adjacent wells SP-05 and
SP-06 in the study. General workflow of the 3D interpretation and inversion for
the Phase 12 of South Pars gas field is as follows.
Page12
Fault interpretation
Depth conversion
Page13
2.8.5 Deliverables
At the end of the geophysical study, consultant should deliver all generated
data and maps. Deliverables should be including the following items:
2.9
Digital database
Litho-stratigraphy
Page14
The average thickness of the Kangan Formation is about 193M in the South
Pars. It consists of two major reservoir layers, KG-1 and KG-2 (K1 & K2):
f The average thickness of KG-2 (K2) reservoir layer is about 43M. It
includes two litho-units that are from bottom to top as follows: porous to
tight limestone to dolomitic limestone unit (KG_2-2) and porous to tight
dolomitic limestone grading upward to dolomite unit (KG_2-1).
f The average thickness of K1 reservoir layer is 111M in the South Pars. It
consists of three litho-units that are from bottom to top as follows: tight
anhydritic dolomite (KG_1-3), porous to tight dolomitic limestone (KG_12) and tight anhydritic dolomite (KG_1-1).
The Nar-Upper Dalan boundary should be picked where the dominant
lithology changes from dolomite dominated to anhydrite dominated. Therefore,
particular care must be exercised for the chronological correlations.
According to the bio-stratigraphic and stable isotope studies performed to
date, the Dalan top (the onset of Triassic age and Kangan Formation) is
characterized with the appearance of a stromatolite / thrombolite unit over a
breccia interval.
Thickness variations of the Upper Dalan and Kangan are minor over the field,
both in the layer-scale and the reservoir-scale as indicated for some of the
vertical wells in table 2.
Review of all available core descriptions indicates that lateral variations of
litho-facies are minor over the field. This is well shown and confirmed by the
correlation of electro-facies (rock types) and petrophysical multi-mineral
evaluations for the vertical wells of the field.
2.10
Petro-Facies Classification
2.11
Sedimentary Environment
Page17
Page18
2.12
Zonation
A new zonation was first assessed on wells SP-05 and SP-06 and at any step,
correlations were made to all the other vertical wells to ensure having a full
view over the field. In order to respect the sequence stratigraphy hierarchy,
the new geological-reservoir zonation includes four major zones, which are in
fact the same as those mentioned above (K1 to K4). These four major
reservoir zones are named KG1, KG2, UD3 and UD4 whose prefixes refer to
the Kangan and Upper Dalan Iranian nomenclature. Top and bottom of all
these zones are limited to sequence boundaries, except the boundary
between KG2 and UD3 (K2-K3) that is a maximum flooding surface. This
issue originated from the sequence stratigraphy study of TOTAL, in which no
evidence of sequence boundary is recognized within the Permian-Triassic
boundary interval. Hence the nearest stratigarphic top to the PTB interval has
been used as the startigraphic marker to differentiate the two major units.
Page19
In the new zonation, each of four major zones are subdivided onto minor
zones, so that K1 is divided into 3, K2 into 2, K3 into 5, and K4 into 6 minor
zones. Table 3 presents the hierarchy of the zonation along with the type of
boundaries and dominant lithology of each zone. The nomenclature used for
these 16 minor zones are as follows:
KG1
KG2
(KG2-1, KG2-2)
Page20
2.13
Litho-typing provides different rock types with specific reservoir and geological
properties necessary for Geo-cellular Modeling. The data required for lithotyping include core data and well logs. To compare the two sources of data,
Petro-facies determination was a prerequisite. An unsupervised litho-typing
applied in this study using a new method called MRGC (Multi-Resolution
Graph-based Clustering).
Permeability prediction is a necessary step of each reservoir characterization
and modeling and MRGC method proved to be the best technique to
permeability prediction in this study.
2.13.1 Litho-Typing
Litho-typing is a method to determine different rock types present in the
reservoir using log and core data in reference wells. The rock types can then
be predicted in all other wells having the same set of log data. This is done
through the data analysis of log and core data and matching the two in the
reference wells. The model built for the reference wells is then used to
analyze the log data for other wells. For each depth the model allocates the
rock type analogous to the log values at the same depth.
Based on the abundance and quality of the data available, Unsupervised
Approach method was applied in litho-typing. In this method the data structure
is analyzed to define the number and the limits of the statistically
homogeneous clusters (electro-facies). Then the defined electro-facies are
having close reservoir and / or lithological characteristics or properties are
merged and finally, the remaining electro-facies are introduced as main types
of the reservoir. Carbonate rocks with different textures, i.e. mudstones to
grainstones, can not be easily differentiated by the responses of the
conventional logs. They mainly differ in porosity, but this is not characteristic
enough for the Supervised litho-typing. The best evidence of the complexity of
Page21
the South Pars Gas Reservoir rocks is the complex and varying relationship
between porosity and permeability.
Among the variety of the methods for electro-facies determination such as
Fuzzy Logic in Fuzzy Logic Module and, Ascendant Hierarchical Clustering
(AHC), Dynamic Clustering (Dyn Clustering) and Self-Organization Map
Clustering (SOM), in Facimage Module in Geolog 6.5, the new method
called Multi Resolution Graph-based Clustering (MRGC), based on the
following advantages, proved to be the most efficient:
MRGC begins by making small clusters and after that it merges them to
form larger clusters. This technique produces a number of models for
the data set, each representing a different resolution. Thereby the user
can choose between the models the resolution required for the project.
key wells in this study but just some wells could be used to Learning due to
software / hardware limitations.
As twelve core-facies had been already distinguished during the study of core
thin-sections, and considering the wide rang of porosity values, the suitable
number of electro-facies was expected to be about 12 or few more. Different
ranges of maximum and minimum number of clusters were tested and
compared. It was concluded that as MRGC clustering does not take into
account the geological considerations by itself, to improve the conceptual
quality of the litho-typing, some geological concepts had to be applied through
user. To do so, it was decided to create models with more clusters (higher
resolutions) than predicted / required and then merge those that geologically
and at the same time statistically alike.
Among five different clustering models (different resolutions) that MRGC
algorithm made available (31, 35, 40, 43 and 49 clusters), the one with 35
clusters was found the most suitable. Figure 11 illustrates 35-cluster model for
wells SP-01, SP-06 and SP-08. As it is seen, it has reached to enough
resolution in order to separate limestones, dolomites and anhydrites for all the
porosity ranges. Figure 12 illustrates 35-cluster model that is applied to all
vertical wells.
35-cluster model bears a resolution higher than what is required for this study.
Hence, as it is done in most of litho-typing studies, in the next step, clusters
with close properties were merged together. This merging step was performed
with this consideration that the lithological (core-facies) character, while the
character should worth defining a new cluster (electro-facies). This was done
referring to the core data and petro-facies.
After the required suitable merges were applied, the optimal number of
clusters was determined as 16 clusters (electro-facies).
Page23
Dolomite
49.38%
Limestone
31.69%
Anhydritic Dolomite
13.76%
applied to all the vertical wells as well as development wells. For Phase 1
wells, due to the lack of RXO data, PHIE & EFAC model was applied. This
model (PHIE, EFAC) was chosen for Phase 1 wells due to its higher CC of
prediction and also due to the fact that both PHIE and EFAC data have been
corrected and quality checked with other wells while other log data in Phase 1
wells show some deviations compared to other wells data. The permeability
log that was predicted this way is called predicted permeability.
3 Petrophysical Evaluations
3.1
3.1.1
Preliminary modeling with single MULTIMIN model for entire reservoir section
showed that the general model would not provide best fit for the
characteristics of KG1-1, KG1-3, UD3-1, and UD4-5. These reservoir layers
are mainly composed of anhydrite and anhydritic dolomite. In order to
Page25
3.1.2
3.1.3
Petrophysical Parameters
Petrophysical parameters were measured and identified for wells SP-05 and
SP-06 in special core analysis reports (TOTAL SOUTH PARS, Kangan /
Dalan formation, special core analysis, wells SP-05 and SP-06, TOTALFINA,
May 2000). In these reports, the cementation factor m is identified as m=2 for
dolomite. It is stated in the report, that for limestone, 'm' parameter may be as
large as 2.4. Another finding of the report is that, the cementation factor varies
with isostatic stress. To investigate the effect of variable m versus constant
Page26
Statistically, near 3000 water saturation and porosity values are calculated in
this evaluation for SP-06. The average water saturation for the entire well
interval by constant and variable m is calculated and shown in the following
table.
For all vertical wells and Phase 2 and 3 wells if the volume of calcite
was over 50%, then the porosity calculated using variable m formula.
Otherwise (if calcite volume was less that 50%), cementation factor m
assumed to be equal to 2.
3.1.4
Core studies show the presence of some shaly and argillaceous intervals in
the South Pars gas reservoir layers. Therefore in petrophysical interpretation,
presence of clay minerals should be taken into account in the model, as a
minor constituent. To model this phenomenon, the Dual Water Model was
chosen for the interpretation. Archie model on the other hand is suitable for
water saturation calculation in clean formations. The advantage of using Dual
Water model versus the other model, Archie, is that, if the formation is clean
(i.e. no clay is present), then the Dual Water formula reduces to the Archie
equation .Another advantage of using Dual Water model is that, this model is
Page28
suitable for brine waters having salinity greater than 20,500 PPM, as it is the
case in the South Pars wells. Archie and Dual Water equations are as follows:
Archie:
Ct = Cw
m
a
S wn
Where:
Cw = Brine Conductivity (S/m)
Ct = Total Rock Conductivity (S/m)
t = Core Analysis Total Porosity (v/v)
Sw = Initial Water Saturation (v/v)
m = Cementation Exponent
n = Saturation Exponent
a = Intercept
Dual Water:
3.1.5
Page30
3.1.6
Gas Correction
3.2
Petrophysical Interpretation
In petrophysical interpretation of South Pars wells, the main steps will be;
Depth
Matching,
PRECALC,
Environmental
Corrections,
Preliminary
3.2.1
Depth Matching
In wire line logging, logs are recorded in multiple runs, and then different sets
of logs can be recorded off depth from one another for a variety of reasons.
So wire line logs were depth matched to RHOB and GR as base logs. Core
depths were also depth matched to NPHI log as the base (tools sticking, poor
motion compensation, engineer error, etc.). In PEX and LWD logging, all logs
are recorded during the same run so, in these types of loggings there will be
Page31
no need for depth matching. In addition to depth matching for logs, core
depths should also be matched with a base log.
3.2.2
Precalc
Down-hole
mud
properties
(Rmf,
Rm,
Rmc)
from
sample
measurements
3.2.3
Environmental Correction
3.2.4
3.2.5
The final MULTIMIN model was applied to all vertical wells. Some changes
and adoptions needed to fit the model with special characteristics of each well
to achieve the best interpretation result. Each well was evaluated using the
final model. As mentioned before, the model for each well was adjusted to
achieve the best quality matching with core-measured lithology and porosity
To cross check the validity of single well interpretations, the calculated
lithology of the wells were correlated across the field by creating crosssections. The cross-sections were selected such that all vertical wells to
correlate with each other from different directions.
In figure 17, some of the vertical wells are correlated. The petrophysical
evaluations of the vertical wells were completed by making complementary
evaluations including lithotype determination and permeability prediction.
3.2.6
The final model for each vertical well was used as the preliminary model for
petrophysical evaluation of deviated wells clustered around the subject
Page33
vertical well. Deviated wells, which were logged with LWD, show a lower
resistivity when compared to wire line logs. These wells needed adjustment to
match the resistivity of the vertical wells after preliminary evaluation.
The same model as mentioned in previous section was used for the wells that
were logged with PEX tools. No sonic logs were available for wells logged
with LWD also the photoelectric logs ran for these wells found to have poor
quality. Due to lack of sonic log and poor quality of photoelectric logs in wells
logged with LWD, the final model could not directly be run in Geolog6 to
compute all minerals and fluid saturations. To resolve this problem, sonic log
was computed (reconstructed) for each well from RHOB by the GardnerGregory equations. To enhance the accuracy of petrophysical evaluation of
deviated wells, it was tried to adjust evaluation of the deviated wells with the
nearby vertical well(s), drilled from the same platform. This is demonstrated in
different cross-sections made between the vertical well(s) and the nearby
deviated wells.
From the cross-sections, calculated lithologies and porosities show a good
correlation between wells for which PEX has been run and the vertical wells,
but to allow acceptable correlation some adjustments were needed for the
wells having LWD logs. These adjustments were necessary to improve the
calculated lithology and porosity.
3.3
Net-Pay Summation
Page34
3.4
For each deviated well MWD / GR is suggested to run from KOP to T.D
and CDR / GR is to run from reservoir entrance point and KOP to
surface.
Full suite logs, including Full-Wave Sonic, TLD, HRLA, HGNS, and
NGS are suggested to acquire for the formations of Kangan and Upper
Dalan for each vertical well under wireline operation.
Full suit logs are suggested to acquire for deviated wells under LWD
operation.
FMI logging for the objective formations including Kangan and upper
Dalan is suggested to perform in all vertical wells.
Page35
Structural modeling
3D Grid construction
3D Facies modeling
Page36
(figures 19 & 20). Upon completing structural modeling and quality control
(QC) of the results, 3D geological grid was built with the following
specifications:
Development Wells
For identifying the cells intersected by the well tracks and giving each cell the
average value of various log properties located in 3D grid cell, Blocking Well
was done. Zone Log (discrete log, in which each sub-zone, the isochors
between specified horizons, is represented by a unique discrete code)
treatment in up-scaling process is shift and scale the logs to match sub
grids. Facies discrete log was up-scaled using weight function to select the
dominant value code to represent each cell. Continuous logs including PHIE,
Perm_Log, SWE, VOL_ANH, VOL_CAL, VOL_DOL and VOL_ILL were upscaled Using arithmetic average and biased to facies log (EFAC_Disc). To
understand the quality of matching between the raw and blocked well data,
Page37
they were statistically compared for main properties such as Facies, PHIE,
SWE, etc.
The priority in a facies model is the treatment of the larger scale architecture.
Studies and measurements show that variability in petrophysical properties is
often greater among bodies of different facies types than among facies bodies
of the same type. The way in which the petrophysical properties vary, is often
the function of the facies type. A proper facies model is therefore required in
order to obtain a sound and possibly detailed petrophysical model for
heterogeneous reservoirs. In the current study, a pixel-based simulation
technique is selected using Facies Indicator module from software for facies
modeling. Then Statistical analyses of facies data in blocked wells were done
considering:
Comparing distribution of facies raw and blocked well data to check the
correct up-scaling process
Distribution of each facies type in any / all sub grids to see the
frequency and presence of facies in the sub grids
The results showed that there was good match between raw and blocked well
data but, no specific trend was found in any directions; therefore vertical
proportion curve for South Pars blocked well data was created as shown in
figure 21. The VPC was used as facies trend and volume fractions for the
facies modeling in each sub grid. 16 variorums created for each facies type
and facies stochastic simulation was done using facies indicator module
subsequently (figure 22).
The sequence of stochastic property modeling consists of the following steps:
Page38
The simulation algorithm used for creating the realizations was based on the
assumption that the variable to be modeled has a Gaussian distribution with a
mean value of zero and a standard deviation of one. It was therefore
necessary to make appropriate transformation to the data to remove any
trends and matching distribution requirements. In next step, transformation of
porosity and permeability was applied for each type of facies separately using
the following sequences:
Page39
4.4 Volumetric
One of the main objectives of preparing the 3D geological model is to have
the best estimate of Initial Hydrocarbon Inplace. Different values for the Initial
Gas Inplace (IGIP) were made available using different statistical realizations
results for reservoir properties.
Page40
5 Reservoir Engineering
5.1
Introduction
The South Pars-North Field straddles the Qatari-Iranian maritime border in the
central part of the Persian Gulf. It is believed to be the worlds largest gas
field, with relatively dry, sour gas in the Permo-Triassic Kangan and Dalan
carbonates. The Qatari part of the field was discovered by Shell in 1971 by
the North West Dome-1 well and the South Pars field was subsequently
discovered by N.I.O.C. in 1991 (South Pars#1 so called SP1).
The field structure is a relatively simple, NE-SW elongated flat dome,
approximately 120KM x 80KM of which 3700KM2 is located in Iranian side of
the field. The reservoir consists of about 400 meters of limestones and
dolomites with subordinate anhydrite, divided into the Dalan Formation
(Permian) and the Kangan Formation (Triassic), which together comprise the
Khuff Formation of Arabian terminology. These formations are traditionally
further classified into 5 main deposition sequences named K1 to K5.
Additionally a tilted gas-water contact towards the east is observed, which is
believed to have been caused by hydrodynamic processes.
The reservoir quarelity is reasonable for gas, with an average porosity of
around 8% and an average permeability of around 30MD. Regionally, lateral
correlation of the depositional sequences is good. However, even for
carbonates, the reservoir properties are highly heterogeneous, both vertically
and laterally, as a consequence of early and late diagenetic processes.
Page41
The development of Phase 12 of the South Pars field includes all fourreservoir units K1 to K4 wherever they are available. The development is
based on an average plateau production rate of 3000MMSCF wellhead gas
(1000MMSCF per each platform).
These data have been worked through in-house dynamic reservoir model that
was constructed to establish the likely range in production performance for the
proposed development scheme. This dynamic model was constructed from
fully integrated geological model that incorporate all data acquired over the life
of the field till a year ago. The potential impact of adjacent developments such
as development Phases 13 and 11 has also been taken into account.
Technical evaluation indicates that the required production rates can be
achieved with the proposed scheme consisting of 36 production wells, drilled
from three single wellhead platforms. Subsequent satellite platform and
compression may be required to maintain plateau production that has to be
considered inside the scope of this MDP.
Plateau length is predicted to be 7, 11 and 12 years for Platforms SPD16,
SPD15 and SPD14 respectively without any capacity maintenance operations
(either satellite platform or compressor).
SP-13 has been drilled at the center of the phase 12 area and alleviated some
vague points in this area. On the modeling it has been assumed that there will
be a well delivery strategy that shall aim at drilling and completing of reliable,
fit-for-purpose, high-potential wells, at optimum life-cycle value, in a timely and
HSE-conscious manner and applying field proven technology. For this 7"
nominal monobore completion design has been taken for VLQ calculations
planned with commingled production from the K1, K2, K3 and K4 reservoir
units. On the areas the following nominal casing configuration is considered
as typical: 26" CP x 18 "x 13 "x 9 " 10 " x 7" monobore completions.
5.2
Data Availability
Page42
Along with drilling 19 appraisal wells and more than 80 development wells in
different areas of South Pars, a lot of basic information has been gathered
that would be very helpful in better knowing of reservoir. These data are
classified as:
5.2.1
Core Data
5.2.2
Well tests and pressure data provides adequate information for reservoir
characterization and simulation. The available DSTs gathered from appraisal
wells are presented in table 12.
5.2.3
More than 60 samples from individual Kangan / Dalan layers in SP-1 to SP15, SPD2-04, SPD4-01 and SPD10-08 have been collected. The summary of
the available results is as come in tables 13 & 14.
Page43
5.3
Average Net to Gross (NTG) ratio: 0.89 (for phases 12: 0.81)
Average Water Saturation for pay zone: 0.15 (for phase 12: 0.18)
C5+ are split into the four heaviest components and represent the
condensate
The 7 pseudo components and reservoir compositions are shown in table 15.
Also the Peng-Robinson EOS with Peneloux Volume Correction parameters is
shown in table 16.
There is a significant spread in the observed GWC levels throughout the field,
with the GWC being some 2000FT deeper in the East-Northeast direction of
the field. The GWC was determined from log saturation profiles, RFT / MDT
data and well test data. The uncertainty range on the gas water contact
determination is + 20 meters.
The gas-water contact was firstly determined by using the test data to define
the "gas down to" and 'water up to" limits and then within these limits,
examining the log water saturation profiles. The most westerly well, SP-14
appears to be water-bearing in the K4 and K3 with proven GDT from MDT at
2853MSS in the border of K2 / K3.
Both SP-1 and SP-2 have tested gas in the K4 with, at slightly deeper levels,
reasonably sharp increases in water saturation which was interpreted as
being the presumed gas down to level. The proven GDT in SP-3 is within the
Kangan-Dalan 3 (K3) while the suggested contact from the water saturation
profiles is just above the K4 / K3 boundary.
The level of GWC in SP-14 and SP-15 which were drilled recently showed
that the closure of reservoir is going to be closed in north and west part of field
by SP-15 and SP-14 respectively.
The differing gas-water levels in the South Pars wells do show, however, an
apparent deepening trend eastwards. This could be due to an easterly flow of
aquifer waters which has resulted in the hydrodynamic tilting of the gas-water
contact.
To model the tilted gas water contact the entire pane were segmented into 45
regions to cover the entire reservoir closure. For each segment the
appropriate capillary pressure curve introduced to model the OGWC in correct
position.
The estimated reservoir pressure is based on the RFT / MDT survey carried
out in appraisal wells SP-1 to SP-15, which gives a pressure gradient of
0.114PSIA / FT (figure 84).
Page46
5.4
The calculation of the Gas Initially in Place (GIIP) volumes for Phase 12 and
South Pars as a whole assumes that the K1 to K4 layers are a single pressure
system and are in pressure equilibrium in the gas phase, and that the tilted
GWC is common to all layers. The distribution of GIIP between the reservoir
layers in the said Phase is shown in following table:
In Place Volume
Reservoir
Kangan
Dalan
Layer
K1
K2
K3
K4
Total
Phase 12
GIIP
Condensate
(TSCF)
(MMSBBL)
5.6
162
4.6
195
6.63
280.5
17.5
773
34.33
1410.5
5.5
Model Initialization
The model has been initialized in equilibrium conditions with the following
parameters: At GWC depth the capillary pressure is equal to zero and the
static pressure at this level is provided as data (i.e. @ 9514M the initial
pressure is 5300PSIA). The pressure and fluid saturation above the GWC is
Page47
populated following the capillary curves corresponding to each cell, its height
from the free water level and the fluid density. The field also is divided into
four vertical equilibrium regions for different gas compositions in KG1, KG2,
UD3, and UD4 reservoir layers.
5.6
5.7
Production Scheme
Vertical lift curves were implemented for Phase12 with a minimum THP of
900PSIA for onshore compression and a minimum THP of 1800PSIA for no
compression. Formation water production is limited to 2000BBL / DAY per
platform. Gas production is limited to 1000MMSCFD per platform according to
the platform design specifications. Based on individual well gas production
potential, Phase 12 wells are set up in order to achieve contractual area
production target without violating platform constraints.
5.8
Well Characteristics
The wells for the proposed development scheme have the following
characteristics:
Page48
All other wells deviated and placed on a specific pattern (figure 85) with
a maximum horizontal departure from vertical to TD of 3500M and a
maximum deviation of 64DEG
Acid stimulation
5.8.1
Well Perforation
To maximize the well productivity and condensate recovery, the wells are
perforated through K1 to K4 (Up to GWC) for commingled production from all
reservoir layers.
5.8.2
Well Simulation
VFP tables have been calculated by VFPi and are used in ECLIPSE to
monitor the wellhead pressures during the runs and decide the
compression date and end of plateau.
5.8.3
Well Constraints
Well step out from the main platform should not exceed 3500M
Introduction
The base plan is to drill a total of 36 wells from 3 platforms in Phase 12.
Within mentioned phase drilling operations are planed to be carried out from
the jackets with a temporary drilling deck, using cantilevered jack-up drilling
rigs. At least two rigs will be required, though they could start at different
times, depending on jacket installation timing and rig availability.
6.2
Appraisal Well
The vertical well which in turn is the first well drilled on each platform, shall be
treated as learning well. The following data gathering program is planned for
closer acquaintance with development area in learning wells:
Page50
Wire line logging: The program should include the basic logging while
drilling tool (GR and Deep Resistivity). Across the reservoir section
following tools should be run:
Well test program: Individual production tests for the K1, K2, K3 and K4
reservoirs are planned for these appraisal wells.
Page51
VSP walkaway program: A well seismic survey will add value by using
new velocity data (the former velocity model was based on regional
geology considerations) to calibrate phase 12, 2D seismic lines and the
South Pars 3D seismic cube. The calibration is required to improve the
success of 3 dimensional volume visualization techniques for reservoir
quality prediction.
6.3
Development Wells
A limited data gathering program is proposed for these vertical and the
subsequent deviated development wells. Only data required for adequate well
reservoir evaluation is included in the program as following:
6.4
The following well construction risks and uncertainties have been identified
from the available data:
Page52
Sulphurous water flows from the Jahrum can be source of well control,
health hazard and corrosion attack.
Losses during cementing of the 18", 13" & 9" 10" casing
strings affect the ability to achieve to full zonal isolation. This may lead
to annular pressure build up and migration to surface of sour gases.
Shale problems especially from the Aghar shales, could lead to hole
problems and extended well durations.
6.5
6.5.1
Well Trajectories
Nudge off bellow the conductor to stay outside collision limits and then
build up to 10-15 degrees above the loss zone in the Jahrum
Page53
6.5.2
Casing Design
The actual setting depth of the casing may be different for each of the
individual well in the platform. However, the setting depth will be decided upon
the formation encountered.
Despite learning wells, all other wells will be drilled directionally in order to
reach the proposed displacements at reservoir depth. The following is the
proposed well architecture:
6.5.2.1
The 32" section will be drilled vertically with a standard rotary assembly using
a 17" bit followed by a 32" hole opener and the 26" conductor pipe shall be
set and cemented at approximately 70M below sea bed. Conductor hanging
and offline cementing shall be applied.
The size, wall thickness and grade of the conductor is specified as 26" OD,
0.625" wall thickness, RL-4S or equivalent connector and X52 grade material.
6.5.2.1
The 24" hole will be drilled directionally with two BHAs, with the initial kick-off
in the first run which will drill down to top of Jahrum formation (The Kick off
point will be below the 26" CP). The second BHA, a rotary BHA including
MWD will be used to drill down to the Illam formation, total depth (TD) for the
section.
The 18" casing shall cover this most troublesome hole section where
circulation losses, formation instability and sulphurous water shows are
common problems. The recommended 18" casing is 114LBS/FT and X55
grade with a modified Buttress or equivalent connector.
Page54
6.5.2.2
The 16" hole shall be drilled directionally down to the Hith anhydrite. The 13"
casing shall cover the sometimes thief Cretaceous formations. The casing
shoe shall be set in the competent anhydrite to obtain a good fracture margin
for the next hole phase.
The recommended 13" casing is 68LBS/FT and L80 grade material with a
modified Buttress or equivalent connector.
6.5.2.3
The 12" hole will be drilled down to the top of the Kangan reservoir. The 9"
production casing will be set right at top of the Kangan reservoir, to case off
the Aghar shale in the Lower Dashtak and to allow for drilling into the
reservoir. It will be crossed over to 10" to accommodate the sub surface
safety valve.
The recommended 9" casing is 47LBS/FT and L80 grade material with Vam
Top or equivalent connection. The recommended 10" casing is 55.5LBS/FT
and L80 material with Vam Top or equivalent connections.
6.5.2.4
The 8" hole will be drilled through the reservoir levels (K1-K4) down to the
planned TD. The 7" production liner will be set in the K4 level, 3M above
planned TD. The recommended 7" liner is 29LBS/FT, 28% Cr with Vam Top
or equivalent connections down to the top of K1 level and 7" liner is
26LBS/FT, 13% Cr with Vam Top or equivalent connections to cover K1, K2,
K3 and K4 levels. Top of liner shall feature an arrangement in order to allow
for mono bore 7" production string.
6.5.3
Cementing Design
Page55
The final slurry recipe, slurry volume and operation details will be provided
after field and lab tests results, prior to commence the cementing job, and will
be dictated by the actual operation requirements. The planned TOC for
various casings are listed below:
Page56
The 7" production liner will be cemented up to the liner hanger using neat
slurry. The Shoe Track must be contain of 7" float shoe, 1 joints of 7" liner, 7"
float collar, 3 joint of 7" liner and landing collar.
6.5.4
Completion Design
A 7" mono bore completion will be used. The liner will be tied back to surface
by means of 7" 23 26LBS/FT, CRA tubing. A sub-surface safety valve will
be installed near the surface.
Liner hanger
Top isolation tie-back packer with a tie back polished bore receptacle
7" production tubing will be 23 26LBS/FT CRA, 28% Cr with tight gas
connections, type Vam Top or equivalent
(or split multi bowl), which offers the advantages of faster make-up because of
fewer flanged connections thus exposing fewer leak paths and reducing
height and weight. The Xmas-tree would be 7 1/16", single block 90 degree
bend type, rated for sour service PSL 3G, 6500PSI WP. The lower master
valve would be manually operated, while the upper master valve would be
fitted with a hydraulic actuator with wire cutting capability.
6.5.5
The following mud systems and completion fluid (table 18) are selected for
drilling and completion the wells in the platforms SPD14 to SPD16. However,
the recommended mud types here described are not mandatory, but a
guideline.
The section will be drilled with a KCL / Polymer / Glycol mud system with
additives. The hole section traverses through limestone, clays, chert, marl,
shale, sandstone, siltstone, dolomite and anhydrite. Main hazards are oil
water influx, shale swelling, bit balling, partial to total losses from Kazhdumi
and Fahliyan. H2S may occur in Dariyan and Hith.
6.5.6
Page59
Solids in the mud system shall be kept under control though the most
appropriate use of the solid removal equipment. The optimum, suggested
system solid control system should be composed of:
Mud cleaner
2 centrifuges
6.5.7
6.5.8
Well cleanup will be conducted with the drilling rig on location. The objective
of the well cleanup operation is the removal of completion fluid and spent acid
from the well bore, whilst avoiding contamination of the surface facilities and
the sea. Major hazards associated with well cleanup operations are releases
of sour gas, exposure to heat radiation and spillages.
6.5.9
Page60
Taking the offset performance Timing information and experience based flat
timing gives the following time estimate for a max 3500M displacement
development well (figure 91).
2.
3.
4.
Page61
6.6.1
6.6.2
The full ESHIA to be completed prior to the start of operations in the field
Page62
will address both drilling and construction related issues. Main hazards,
from a well construction perspective, are related to loss of containment,
from the well, the rig or from manse operations. Hazards will be studied to
evaluate impact and to assess mitigations and recovery plans. An
emergency response exercise will be carried out as part of the spud
preparations.
In addition, procedures / programs for basic drilling operations developed and
maintained by the chosen Contractor shall be referenced as the industry
standards.
Page63