Document Type : Research Paper
Authors
 Syed Waqas Haider ^{} ^{1}
 Mustafa Yar ^{2}
 Raja Ahtisham Ghafoor ^{3}
 Tallat Majeed Khan ^{4}
^{1} Lecturer, Department of Geology, FATA University, Darra Adam Khel, F.R. Kohat, Pakistan
^{2} M.S. Student, Department of Earth and Environmental Sciences, Bahria University, Shangrilla Road E8, Islamabad, Pakistan
^{3} Graduates, Department of Earth and Environmental Sciences, Bahria University, Shangrilla Road E8, Islamabad, Pakistan
^{4} Department of Earth and Environmental Sciences, Bahria University, Shangrilla road E8, Islamabad, Pakistan
Abstract
The well Sarai Sidhu01 is located on Punjab Platform, Central Indus Basin, Pakistan. Punjab Platform is the eastern part of Central Indus Basin, and tectonically it is the stable portion of Indus Basin, which was least affected during Tertiary Himalayan orogeny. This study attempts to decipher reservoir potential for hydrocarbon exploration. It aims to delineate a subsurface hydrocarbon bearing zone and to estimate the reservoir properties. A complete suite of wireline logs containing Caliper log (CALI), gamma ray log (GR), spontaneous potential log (SP), neutron log (Ø_{N}), density log (Ø_{D}), and resistivity logs (MSFL, LLS, and LLD) with all drilling parameters and well tops were utilized. The methodology adopted to accomplish this task includes the calculation of volume of shale (V_{sh}) by using gamma ray log and effective porosity (Ø_{E}) by using density and neutron logs. Resistivity of water (R_{w}) was calculated by SPmethod, and the saturation of water (S_{w}) and the saturation of hydrocarbons (S_{h}) is calculated with the help of Archie’s equation. According to log signatures, Lumshiwal formation of early Cretaceous age encountered in well in the depth range of 5433 ft. to 5797 ft. was marked as a possible reservoir, and this zone was evaluated for its reservoir potential in detail using a set of equations. The average values calculated for different parameters are as follows: V_{sh}= 30%, Ø_{E}= 17%, S_{w}= 46%, and S_{h}= 54%. The analysis shows that S_{h} is low, so it is inferred that Lumshiwal formation has a low potential and is economically not feasible for hydrocarbons production.
Keywords
Main Subjects
Wireline logging is a practice to attain information about subsurface strata encountered in a borehole drilled for hydrocarbon exploration (Rider, 2002). Its key purpose is to evaluate the potential of a formation to act as a reservoir for accumulation and economical recovery of hydrocarbons.
The present study has been performed on well SaraiSidhu1 located in tehsil Kabirwala, District Khanewal, Punjab province, Pakistan (Figure 1); geologically, it is a part of Central Indus Basin. Sarai Sidhu1 was drilled by Amoco in 197374 on seismically delineated prospect up to TD 15303 ft., and initially lowquality gas was discovered in cretaceous sequence; now, it is an abandoned well (Quadri and Quadri, 1998).
The main objective of this study is the demarcation of hydrocarbon bearing zone and its evaluation for reservoir characteristics by interpreting wireline logs in well SaraiSidhu 01. Reservoir characterization includes the estimation of different physical parameters, i.e. the volume of shale (V_{sh}), porosity (∅), saturation of water (S_{w}), and hydrocarbons.

Figure 1
Tectonic map of the studied area (Kadri, 1995).
2. Regional geology and tectonic settings
The central Indus Basin from east to west comprises of three units (Figure 2); Punjab Platform, Sulaiman depression, and Sulaiman Fold belt (Raza et al., 1989). The well under study lies on Punjab platform (Figure 2). Punjab Platform marks the eastern segment of Central Indus Basin, and it is a broad monocline dipping gently towards Sulaiman depression (Kadri, 1995). In the north, Sargodha high separates Central Indus Basin from the Upper Indus Basin. Sukkur Rift comprises of Jacobabad and Mari KandhKot, and highs are present in the south separating it from the southern Indus Basin. It is locked by Indian shield in the east and by the marginal zone of the Indian plate in the west (Shah, 1977). The area was relatively less affected by Tertiary Himalayan orogenic event because of its greater distance from the collision zone (Raza et al., 1989). The Punjab Platform is covered by thick unconsolidated quaternary alluvial deposits of sand, silt, and clay with the maximum thickness of about 1600 ft. (Kadri, 1995; Raza et al., 1989). The oldest rocks encountered in Punjab Platform through drilling are InfraCambrian Salt range formation (Kadri, 1995; Humayun et al., 1991). Table 1 shows the stratigraphy of the study area encountered in the borehole of SaraiSidhu 01. The Punjab Platform further extends towards east into India, where it is called BikanerNagaur Basin; significant quantities of oil have been discovered from InfraCambrian evaporites therein (Aadil and Sohail, 2014). This discovery enhanced the exploration activities on Punjab Platform.

Figure 2
Tectonic map of Pakistan highlighting Central Indus Basin and its subdivisions (Modified after Kadri, 1995).
Table 1
Borehole stratigraphy of SaraiSidhu1.
Era 
Period / epoch 
Formation 
Lithology description 
Formation tops (feet) 
Hydrocarbon system 

Seal 
Source 
Reservoir 

CENOZIC 
PLIOCENE 
NAGRI 
Claystone, Siltstone 
00.0 



MIOCENE 
CHINJI 
Claystone, Siltstone 
2063.97 




EOCENE 
SAKESAR 
Limestone 
4373.68 




NAMMAL 
Shale, Limestone 
4582.67 




GHAZIJ 
Shale, Sandstone 
4681.75 




PALEOCENE 
DUNGHAN 
Limestone 
5178.80 




RANIKOT 
Limestone, Shale 
5291.99 




MESOZIC 
CRETACEOUS 
LUMSHIWAL 
Sandstone 
5433.07 



CHICHALI 
Sandstone 
5797.24 




JURASSIC 
SAMANA SUK 
Limestone, Shale 
5985.56 




SHINAWARI 
Limestone, Sandstone 
6434.71 




DATTA 
Shale 
6746.71 




TRIASSIC 
KINGRIALI 
Sandstone, Shale 
6807.74 




TREDIIAN 
Sandstone 
7185.03 




PALEOZOIC 
PERMIAN 
AMB 
Limestone 
7355.64 



SARDHAI 
Clays 
7664.69 




WARCHA 
Sandstone 
8202.1 



3. Material and methodology
Data acquired to accomplish this study consist of openhole wireline logs of well SaraiSidhu01 and well tops. Data are provided by Land Mark Resources (LMKR) with prior approval of Directorate General Petroleum Concessions (DGPC) Islamabad, Pakistan; they include Caliper log (CALI), gamma ray log (GR), spontaneous potential log (SP), resistivity logs (LLD, LLS, and MSFL), neutron log (Ø_{N}), density log (Ø_{D}), and well tops along with all drilling parameters. The workflow adopted for logs interpretation is illustrated in Table 2.
Table 2
Methodology adopted for research work.
4. Results and discussion
The evaluation of the reservoir properties of well SaraiSidhu01 was performed by the careful analysis of all the available wireline logs data. A zone is marked on the bases of GR log trend (Figure 3a) in a depth range of 5433.07 ft. to 5797.24 ft. as a possible reservoir. The total thickness of zone is 364 ft. Well tops of SaraiSadhu 01 (Table 1) confirmed this zone as Lumshiwal formation of a cretaceous age. Further analysis shows that the effective porosity calculated using neutron and density logs for the zone of interest is 17% (Figure 3c). The average saturation of water calculated using Archie’s equation is 46%, and the saturation of hydrocarbon is 54% (Figure 3c). After applying cut offs to minimum extractable hydrocarbons, net pay thickness is calculated 150 ft. for the whole zone in different depth intervals. Figure 3 shows the combined graphs of all the calculated parameters for the reservoir evaluation of Lumshiwal formation.
A zone of interest, having the fair potential to act as a reservoir for the accumulation of hydrocarbons, is marked with gamma ray log and correlated with the other logs in a depth range of 5433.07 ft. to 5797.24 ft. (Figure 3). This zone is tested for its hydrocarbon potential using a set of equations.

Figure 3
Composite graphs for the complete evaluation of our zone of interest; A) GR log curve in our zone of interest; B) the composite graph of V_{sh} and V_{sand}; and C) the combined graph of porosity, S_{w}, and S_{h}.
4.1. Calculation of volume of shale (V_{sh})
V_{sh} is also known as the dirtiness of reservoir, and GR log is used for the calculation of V_{sh}using the following equation (Rider, 1996).
(1) 
where,
 GR log = Gamma ray log reading.
 GR max = Maximum gamma ray deflection in the zone.
 GR min = Minimum gamma ray deflection in the zone.
4.2. Calculation of clean volume (V_{sand})
V_{sand} is the cleanliness of a reservoir, known as the net to gross ratio. V_{sand} is obtained by subtracting 1 from V_{sh}. Equation is given below.
(2) 
4.3. Porosity calculations
Porosity analysis is an important factor for reservoir evaluation. In the present study, neutron and density logs are used for porosity calculations. Neutron log (Ø_{N}) measures hydrogen ion present in formation, so it can be calculated directly from log curves; density porosity is calculated using the following equation (Rider, 1996).
(3) 
where,
 RHOM = Matrix density (constant values according to lithology)
 RHOF = Fluid density (from log header)
 RHOB = Bulk density (density log values from log curve)
Average porosity is calculated using combined neutrondensity logs with the following equation (Rider, 1996).
(4) 
where,
 Ø_{N} = Neutron porosity
 Ø_{D} = Density porosity
(5) 
where,
 Ø_{E} = Effective porosity
 Ø_{Avg} = Total porosity
 V_{clean} = Volume of sand
4.4. Resistivity of water (R_{w})
The calculation of R_{w} is very sensitive because the saturation of water depends on it. Spontaneous potential (SP) method is adopted for the determination of R_{w}; calculations involve the determination of SSP, which is calculated directly from the SP log curve. Formation temperature is also determined by using Equation 6:
(6) 
where,
 BHT = Bottom hole temperature
 Ts = Surface temperature
 TD = Total depth
 Df = Formation depth
Bottom hole temperature (BHT) and the resistivity of mud filtrate at the surface temperature (Rmf @ Ts) were taken from the header of wireline log. To calculate the resistivity of mud filtrate at formation temperature (Rmf @ Tf) the Schlumberger Gen9 chart was used. The Rmf @ Tf is converted to Rmf_{eqv}/Rw_{eqv}using Schlumberger’s SP1 chart (Schlumberger, 2009). Using these values, Rw is determined by Schlumberger’s SP2 chart (Schlumberger, 2009). Table 3 shows the systematic approach to found parameters for determining R_{w} by the SP procedure.
Table 3
Step by step parameters to be found for the determination of R_{w}.
4.5. Saturation of water
Archie’s equation is applied to the calculation of S_{w} (Rider, 1996).
(7) 
where,
 S_{w}= Saturation of water
 n = Saturation exponent
 a = Tortuosity factor
 R_{w} = Resistivity of formation water
 ∅ = Porosity
 m = Cementation exponent
 R_{t}= Resistivity of uninvaded zone
a is a tortuosity factor, which varies from 0.91 to 1.05, and it is normally taken as 1 (Elias and Steagall, 1996). m is cementation exponent, which ranges from 1.8 to 2 for sandstones (Salem and Chilingarian, 1999b). n value also ranges from 2.40 to 4.24 (Elias and Steagall, 1996).
4.6. Saturation of hydrocarbons
The determination of hydrocarbon saturation depicts the potential of reservoir to produce economically feasible hydrocarbons. It is calculated by Equation 8 (Rider, 1996)
(8) 
Saturations of water and hydrocarbons (Figure 3c) are opposite to each other; in other words, where the saturation of hydrocarbons increases, the saturation of water drops. The saturation of hydrocarbons calculated for our zone of interest is 30%, which is very low for an economic recovery of hydrocarbons.
5. Cut off
Applying cut off is a manual practice for petrophysical evaluation. Cut offs are employed to find an average value in the zone of interest, where we find the least extractable volume of hydrocarbons. Cut off grades used to determine economically feasible zone for hydrocarbons production is 30% for V_{sh}, 10% for porosity, and 50% for saturation of hydrocarbons.
6. Conclusions
The analysis shows although the porosity is high, the saturation of hydrocarbon is least economically feasible for recovery. Therefore, it is concluded that Lumshiwal formation of early Cretaceous age has a low potential for an economic recovery of hydrocarbons due to the low volume of hydrocarbons present in pore spaces.
Acknowledgements
We are extremely thankful to the Land Mark Resources (LMKR) and Director General Petroleum Concessions (DGPC) Pakistan for giving the data of the required wells. Furthermore, we are grateful to the Department of Earth and Environmental Sciences, Bahria University, Islamabad for providing us with the facilities to complete this project.
Nomenclature
CALI 
: Caliper log 
GR 
: Gamma ray log 
LLD 
: Laterolog deep 
LLS 
: Laterolog shallow 
MSFL 
: Microspherical focused log 
Rw 
: Resistivity of water 
SP 
: Spontaneous potential log 
S_{w} 
: Saturation of water 
V_{sh} 
: Volume of shale 
Ø_{D } 
: Density log 
Ø_{E} 
: Effective porosity 
Ø_{N } 
: Neutron log 
Aadil, N. and Sohail, G. M., 3D Geological Modeling of Punjab Platform, Middle Indus Basin, Pakistan, through Integration of Wireline Logs and Seismic Data, Journal Geological Society of Indian, Vol. 83, p. 211217, 2014.
Archie, G. E., The Electrical Resistivity Log as an Aid in Determining Some Reservoir Characteristics, Transaction of American Institute of Mining, Metallurgical and Petroleum Engineers (TAIME), Vol. 146, p. 5462, 1942.
Elias, V. L. G. and Steagall, D. E., The Impact of the Values of Cementation Factor and Saturation Exponent in the Calculation of Water Saturation for Macaé Formation, Campos Basin, SCA Conference Paper Number 9611, 1996.
Humayun, M., Lillie, R. J., and Lawrence, R. D., Structural Interpretation of Sulaiman Folds Belt and Foredeep, Pakistan, Tectonophysics Vol. 10, p. 299324, 1991.
Kadri, I. B., Petroleum Geology of Pakistan, P. 300, Pakistan Petroleum Limited, Karachi, Pakistan, 1995.
Raza, H. A., Ahmed, R., Alam, S., and Ali, S. M., Petroleum Zones of Pakistan, Pakistan Journal of Hydrocarbon Research, Vol. 1, p. 2156, 1989.
Quadri, V.N. and Quadri, S.M.G.J., Failuretosuccess Targets May Lie in Pakistan Basins, Oil and Gas Journal, Vol. 96, No. 12, p. 9297, 1998.
Rider, M. H., the Geological Interpretation of Well Logs, John Wiley and Sons, New York, 1996.
Rider, M. H., the Geological Interpretation of Well Logs, 2^{nd} Edition, Sutherland, RiderFrench Consulting, 2002.
Salem, H. S. and Chilingarian, G. V., The Cementation Factor of Archie’s Equation for Shaly Sandstone Reservoirs, Journal of Petroleum Science and Engineering, Vol. 23, P. 8393, 1999b.
Schlumberger, Log Interpretation Charts, Gen9, SP1, SP2 Charts, p.1920, 2009.
Shah, S. M. I., Precambrian and Paleozoic, Stratigraphy of Pakistan, Memoir, Vol. 12, p. 129, Geological Survey of Pakistan, Quetta, 1977.