This study pertains to the evaluation of shale gas and rock physics properties of this area with respect to its total organic content of Sember Formation, Khiproarea, Pakistan. We use well logs data for this study. The Khipro area is prominent in the Lower Indus Basin for its hydrocarbon (oil and gas) structural traps. In shale gas evaluation, TOC of Sember Formation is estimated. The analysis has been done with the help of the wire line data of the well Bilal North-01. The presence of shale gas in the study area is analyzed with the help of different techniques. Rock physics and petrophysical analysis have been done in order to get the properties of the area related to the shale gas evaluation.
Muhammad Sohail Khan1*, Zainab Bibi2*
1School of Earth and Space Sciences, University of Science and Technology of China, Hefei, China. 2Department of Earth Sciences, Quaid-i-Azam University, Islamabad, Pakistan
Received 27 June 2016; accepted 14 August 2016; published 17 August 2016
Copyright © 2016 by authors
The study area, Khipro Block is tectonically situated in Southern Indus Basin and geographically in the Sanghar district Sindh province of Pakistan. Southern Indus Basin is very important for hydrocarbon exploration because it is the main producing oil and gas basin of Pakistan. Khipro is located in the Thar platform area of Southern Indus Basin. The Khipro area is characterized by a series of horst and graben structures present almost below the base Paleocene unconformity within the Cretaceous. The geographical co-ordinates of the area are listed:
- Latitude 25˚20’N to 26˚30’N.
- Longitude 68˚20’E to 69˚14’E.
Rifting and breaks up Gondwana land in Jurassic period is responsible for the formation of Khipro block baisan. The East Gondwana plate (India, Antarctica and Australia) separated from the West Gondwana plate (Africa and South America) in the Cretaceous period. In Aptian time (120 Ma), the Indian plate separated from east Gondwana  defines in article “A Speculative Tectonichistory of Pakistan and Surroundings” that at the end of Cretaceous, Seychelles and Madagascar separated from India with associated faulting resulting in basaltic flows (Deccan Volcanism) in the southern part of Lower Indus basin. After Paleocene there was continuous oblique convergence of Asian plate and Indian plate throughout Tertiary time and the collision results in tilting of the entire region. Deposition during the rifting is shown by the presence of Jurassic rocks in the area.
Due to rifting normal faulting and horst and graben structures are formed. The famous among these structures include “Sukkur Rift”. However, this localized rifting phase was unable to continue after the Paleocene-Eocene time . Below is the Geographic location map of study area  (Figure 1).
Figure 1. Geographic location of study area (Banks and Warburton, 1986).
2. Objectives and Stratigraphy of the Area
The studies carry out to identify the zone of organic rich shale from non-source shale and to estimate TOC of Sember Formation in area under study, using methodology based on sonic and density log with the combination of resistivity Passive method . Bilal North 01 well with complete well logging data is utilized for this purpose. Then to observe the effect of TOC on rock physics parameters to determine the shale characterization of Sember Formation of Lower Indus Basin Pakistan.
It is very critical to have knowledge about the stratigraphy of the area for the hydrocarbons prospecting, by this knowledge it is determined that what are the source, reservoirs and seals rocks of the area. The lithological setting and stratigraphic sequence of the study area which is our interest zone is as follows: Goru Formation, Upper Goru, Lower Goru, Sember Formation. Generalized stratigraphic column of the study area  (Figure 2).
Figure 2. Generalized stratigraphic column of the study area (Quadri, 1986).
3. Data and Methodology
For this research was acquired as part of a project to reassess the oil and gas exploration potential through a regional potential field data survey of the area and a seismic data survey of Khipro area Pakistan. To achieve the desired objective we used following methodologies.
- Reservoir characterization.
- Shale gas reservoir characterization.
- Well log response in shale gas.
- Methods to compute TOC.
- Rock physics.
- Model for shale gas.
4. Shale Gas Reservoir Characterization
In conventional petroleum system, carbonates and sandstones form the reservoir rocks and shales act as a source and seal rock. Since the last decades, shales have also been emerged as a good reservoir because of high porosity and very low permeability. Thus shales are considered as source rock as well as reservoir rock and seal too in terms of unconventional shale gas reservoirs. These shales have variability in term of mineralogical, elastic and physical properties from basin to basin and even within the basin . The characterization of conventional reservoirs are well established from the seismic and wireline logging techniques, but in case of shale gas reservoirs, thickness, Saturation, porosity, fracture, brittleness, etc. are the characteristics in the development of shale gas play.
The gas may be trapped as free gas in natural fractures and inter granular porosity, as gas absorbed into kerogen and clay or gas dissolved in kerogene and bitumen where gas may be biogenic and thermo genic and may be mixture of both. The source rock is mainly depending on the estimation of total organic content and level of maturity of source rock and done through the geochemical sedimentological methods. The organic content in the source rock has specific velocity, density and resistivity characteristics which affect log and seismic responses . Thus source-shales and hence shale-gas reservoirs can be characterized through logs and seismic. Kerogen is the insoluble fraction of organic matter which usually constitutes 80% to 90% organic content. Vitrinite reflectance is a key parameter that indicates the maturity of organic material to calibrate the thermal and burial history models (Figure 3) .
Figure 3. LOM vs vitrinite reflectance (Passey et al., 2010).
Factors affecting accurate Ro measurement:
- Rough textured vitrinite.
- Partially dissolved (pitted).
- Oxidized vitrinite.
- Other macerals.
- Too few readings (<20).
Gas producing shales are typically the source rocks that also function as reservoir rock. Apart from low permeability of shales, the key properties for a gas shale reservoirs far a successful shale gas play, the following characteristics need to be considered:
- Total organic carbon (TOC).
- Reservoir thickness.
- Thermal maturity (R%).
- Reservoir characteristics (brittleness etc.).
Shale which have sufficient TOC may generate hydrocarbon at favorable temperature and pressure conditions. Combination of these factors lead to favorable productivity and the depth of the shale gas formation should also be considered as it will have a bearing on the economics of the gas recovery.
5. Total Organic Carbon and Thermal Maturity
TOC is basically the amount of carbon bound in organic compounds of the rock. TOC has primarily been used to evaluate the quality of source rock although it is now also used in evaluation of unconventional reservoirs.
Table 1. Classification of source potential on TOC content.
TOC and kerogene are closely associated with shales and silt rich claystone however they can also be present in many carbonates. TOC is expressed as a percentage by weight present in rock . Economic shale gas development candidates require sample organic matter to generate sufficient volumes of hydrocarbons i.e., relatively high TOC. TOC values ranges from 2% to 10% have good exploration targets . Classify the shale rocks to be act as potential source of shale gas based on TOC contents. This classification is given in Table 1.
Thermal maturity is one of the most important parameters used in the evaluation of gas-shale and shale-oil plays. Vitrinite reflectance (Ro) is a commonly used thermal maturity indicator. Thermal maturity of organic matter represents the critical element of thermogenic Shale Gas. Vitrinite reflectance is the key factor for estimating the level of thermal maturity of organic matter (LOM). The measurement of Ro was originally developed to rank the coal maturity .
Above 0.6% Ro value is the indicator of immature kerogen not having sufficient thermal conditions for conversion of organic material to hydrocarbon shown in Table 2. The thermal maturity determines whether a source rock can produce oil, gas or condensate .
Table 2. Interpreted maturity stages (Alexander et al., 2011).
6. Well Log Response in Shale Gas
Organic-rich rocks are usually composed of three components i.e., rock matrix, solid organic matter and the fluid(s) filling the pore space, while non-source rocks are composed primarily of only two components i.e., the rock matrix and the fluid filling the pore spaces. As the source rock matures, a portion of the solid organic matter is transformed to liquid (or gaseous) hydrocarbons which move into the pore space, displacing the formation water. These physical transformations affect the porosity, density and resistivity of the rocks and hence reflected in their electric logs. To identify the organic rich rock, gamma ray log gives high values at the interval. This tool distinguishes the clay rich strata and sandy layers. The relationship between Natural gamma ray log and organic matter has been described by .
The quantitative interpretation of natural gamma ray log at organic rich interval suggests that such intervals do not have a high gamma ray value. So from analysis it is resulted that uranium con tent cause the high gamma ray value absorbed by the organic matter in marine sediments. Recently spectral gamma ray tool has been used to identify the organic matter because of empirical relationship between uranium and organic matter. Resistivity of rocks is depending on the fluid present in pore spaces of rock. In conventional reservoirs formation water is electrically conductive at least when he fluid is saline that allows the ionic conductivity. In these formations the resistivity value is very low due to the large volume of the formation water.
When the formation water is displaced by the significant amount of hydrocarbon is non-conductors in the shale gas reservoirs, the response of resistivity log is dramatically increased in organic rich source rock that indicates the generation of hydrocarbon. The resistivity log depends on the free hydrocarbons present in the pores and not on a solid organic matter present in it. The response of resistivity log is also associated with the degree of maturity that implies the amount of hydrocarbon in pores. The presence of low density grain (kerogen) organic matter, the density log response at organic rich source interval gives the low value.
The normal average value of density in mixture of clay minerals is about 2.7 g/cm3. Thus the overall organic rich shale bulk density marked the effect on density log as low value. Combined with the resistivity log, density log give the extremely valuable results in organic rich source interval. Due to the presence of organic matter such as hydrocarbons in shale, the sonic log response give the low value of P-wave velocity due to the high transit time value similar to the density log response. Shear wave velocity response is also affected by the TOC.
The use of sonic log can be enhanced to determine the TOC when combined with the resistivity log. Neutron log response gives the high value of porosity along the organic rich interval. But the log is not good indicator of organic matter as single tool because it is not affected by only the presence of hydrogen in organic matter, but also by the hydroxyl ion (OH−) in clay minerals, hydrogen present in the formation water and even any form of hydrocarbons present in pores. For shale gas formations, the overlay of density/neutron is useful for identification of the higher gas volume but this technique is only applicable when the formation is not clay rich.
7. Methods to Compute the TOC
Shale gas reservoirs have the capability of self-generated and self-stored hydrocarbons. Due to the unique physical properties of organic matter shale gas reservoirs, the well logging response will be definitely different relative to the non-hydrocarbons shale. The well logging response in these reservoirs shows the high acoustic transit time, high resistivity, high Gamma-ray and low density. As a result, integrated with certain techniques, using the different well logging response between non-hydrocarbon source rocks and hydrocarbon source rocks to identify and calculate TOC is available .
The TOC is the fundamental indicator for the evaluation of shale gas reservoirs. Usually the TOC is the measured by the cores, cutting or sidewall cores in the laboratory with source rock evaluation instruments. Due to the limited number of core samples, the experimental results are not continuous and it is impossible to show the whole face of a source rock bed. So the continuous logging tool information contributes to show whole face of source rock bed shown in (Figure 4).
Figure 4. Logging results of Bilal North-01, track description is as follows: Track 1: Depth; Track 2: Gamma Ray; Track 3: Resistivity Log; Track 4: Neutron Log; Track 5: Sonic Log; Track 6: Density Log.
From well logs, several methods have been published to calculate the TOC. Four common methods for TOC calculation such as 1) directly by regression of core TOC with core bulk density; 2) ΔlogR technique; 3) based on the petro-physics response model; 4) using an artificial neutral network). The bulk density regression method gives the satisfactory results when the formation is fairly homogeneous, but in the case of heterogeneity, this method will show result in poor corelations. In petro-physical model, to determine the amount of TOC/kerogen, this model requires accurate porosity values. Unfortunately the porosity and kerogen content cannot separate from neutron and density log, and the Gamma-ray log readings are affected by the variable kerogen content of shale gas formation.
To determine the effective porosity and shale content from shaleg as interval logs such as NMR and spectral-lite logs are required. The artificial neural network can use one or multi-log curves to estimate other log properties , but the technique needs extensive trial and error to select partial samples to overcome the over-learning problems. In a highly heterogeneous case, the method is difficult to predict the TOC because some meaningful core samples might not be available for the neural network training. The most common method which is widely used for TOC determination in shale gas plays is ΔlogR technique to calculate the ΔlogR based on the separation between the deep resistivity curve and the three porosity logs such as density, neutron and sonic logs can be used, which can be converted into TOC through the parameter of level of organic maturity (LOM).
LOM is obtained from various methods of samples analysis (e.g., vitrinite reflectance, thermal alteration index and Tmax), or estimated using burial and thermal history of the basin. LOM determine the source rock quality either it can produce oil, gas or condensates. The technique of Rock-Eval pyrolysis as a standard method is used in the petroleum exploration . This method computes the thermal maturation of the source rock. Thermal maturity parameters based on the maximum temperature (Tmax) to determine the dimensions of oil windows. A single scale which synthesizes several existing indices of organic maturity based on coal rank. The scale relates LOM to vitrinite reflectance among other indicators. LOM is unit less. We have analyzed our well: Bilal-01 lower Indus basin Pakistan.