Integrated Geophysical Study of the Lower Indus Platform Basin Area of Pakistan ()
1. Introduction
The Lower Indus Platform basin is one of the Pakistan’s tectonostratigraphic provinces, which is geographically located in the Sindh province of Pakistan between longitudes 23˚ and 23˚31 N and latitudes from 66˚E respectively. The location of the Lower Indus Platform is shown in Figure 1 [1]. Pakistan is divided tectonically into two major divisions, the Western Makran basin having an oceanic crust (Arabian Plate) and Eastern Indus basin with a continental crust (Indian Plate) with Murray Ridge/ Owen Fracture as the dividing line between these two [2]. The Eastern Indus Basin is subdivided into three geological Provinces [3]. The Lower Indus Platform basin belongs to the Indus platform and fore deep zone of tectonics framework. The vast sedimentary basins cover about 80% of the total country’s area with widespread oil and gas seepages. Despite the long history of oil and gas exploration and an attractive petroleum policy, the sedimentary basins still contain sizeable areas, classified as frontiers, which cover more than 60% of the total prospective area including the Pishin Basin, the Maluan-I Charan Basin, Fold and Thrust Belt parts of the Platform of the Indus Basin and the entire Offshore [4,5]. Presently, about 421 exploratory wells have been dug all across Pakistan. In the Southern Indus Basin 55 exploratory wells were drilled during the period 1955-1984, but 27 of these wells were based on the results of multifold seismic surveys [6]. In the past hundreds of years of petroleum research, compared to an average of one hundred wells in other similar basins, Pakistan remained one of the least explored basins of which the well density is less than five per 10,000 sq·Kms. The country desperately needs to exploit new oil and gas potentials across the new frontiers to increases self-reliance. For sustainable exploration, more efforts should be directed towards these frontier areas which after exploration may prove to be of greatest potential [7]. This paper briefly highlights the hydrocarbon potential and deep crustal structure of Lower Indus Platform basin deduced from the study of seismic and gravity data. The identification of existing petroleum systems and assessment of undiscovered oil and gas potential analogies with other basins suggest that there is significant potential for many hydrocarbon discoveries present in this basin [7]. But it requires a commitment from both the oil companies and the government to achieve success. New plays concepts need to be utilized with new technology along with advanced geological and geophysical studies that highlight the potential of these areas. This study analyzes this area in order to describe the potentialities of this basin which has emerged
Figure 1. Representing the location of the area (GSP Pakistan, 1982).
in the past few years as a major petroleum province due to many major oil and gas discoveries along with many other discoveries in Pakistan. A variety of structural styles and depositional regimes are associated with this area.
2. Geological Background and Stratigraphy of the Area
The Lower Indus Platform Basin is bounded to the North by the Central Indus Basin, to the Northwest by the Sulaimanfoldbelt Basin and the Kirthar Fold Belt Basin in the Southwest. The area lies in the tectonic zone of the “Indus plate form and fore deep” which comprises of many structural zones. The structural features which are present in the area are tilted fault blocks, thrust-faulted anticlines, and anticlines. The anticlines and thrusted anticlines occur in the foreland portions of the Indus basin as a result of compression related to collision of the Indian and Eurasian plates [8]. This tectonic province is underlain by infra-Cambrian to recent clastics and carbonates. It retained a passive margin until the Late Cretaceous, then became part of the complex suture between the Indian Plate and the Afghan Block. The tectonic activity affecting the structures and sedimentology of the Lower Indus Plate form basin including rifting of the Indian Plate from Gondwanaland (Jurassic or Early Cretaceous), which probably created NE-SW to N-S rift systems, and isostatic uplift or ridge-push at the margins [9]. At Cretaceous-Tertiary boundary during Cretaceous age, the thermal doming, strike-slip faulting, and hotspot activity occurred due to the separation of the Indian and Madagascan plates. This resulted in uplift, erosion, extrusion of the Deccan flood basalts and probably the NNWstriking normal faults [9]. There are Paleocene-Eocene emplacements of the Bela Ophiolites which caused gentle folding, the passive margin conditions of Eocene caused structural quiescence and carbonate deposition [10]. The Himalayan Oligocene to present-day collision caused sinistral transpression in the west of the Lower Indus Basin, due to which with fold-thrust structures overprinted by sinistral flower structures [9]. The stratigraphic sequence of the Middle and Lower Indus Basin with full details is also as shown in the Table 1 [11]. The stratigraphic succession changes from east to west. The Precambrian basement is westward. The important unconformities occur at the base of the Permian and the Tertiary. In the eastern part of the basin, the Tertiary has direct contact with the Jurassic sequence [12]. The Indus Lower Indus Platform basin reflects sedimentation associated with rifting.
Table 1. Representing the stratigraphic details of the area (Kadri, 1995).
3. Data and Research Methodology
Data 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 in the Lower Indus Plate form basin area of Pakistan. We have applied the potential field data to understanding the basement character as interpretation of the basement geology which is important for developing and evaluating oil and gas plays for the deeper parts of the basin as well as understanding the regional tectonic history. Integrating the potential field data has helped in refining the interpretation of the prospective part of the basin. The research methodology is carried out by using Interpretation of the potential field data and seismic data of the Lower Indus Platform basin area of Pakistan.
4. Interpretation of the Potential Field Data
The potential field data in the study area consists of regional bouger gravity, residual gravity, total magnetic intensity, second derivatives regional maps and a series of prospect-oriented gravity and magnetic profiles or maps. The regional maps represented in Figure 2 the gradient of bouger anomaly gravity, isostatic anomalies and magnetic depth and basement elements. Through the available data it was concluded that the Lower Indus Platform basin represents anomaly appears that have reflect actual density intensity data shown as distribution of rocks. It can be concluded that there is steep positive gravity gradient towards the southwest above thick sedimentary strata and the depression areas attributed to a rise in the mantle while the depth model indicates that the overall gravity effect is typical. The crustal thickness is more than 20 km thick towards the Northeast and is reduced to 7 Km on the mantle rise. The available intensity magnetic data shows magnetic field variations from over −100 is the Southwest to −50 gamma over much of the North Eastern area [13], which reflect seismic structure in south and southeastern areas. The crust is transitional. The presence of thick sedimentary rocks of various lithologies which commonly form source-reservoir seal trilogies and optimum geothermal gradients analogies shows that there are a number of structural traps whose age range from Cretaceous to Miocene and provide evidence of hydrocarbon. Due to thinning of the crust, the younger sediments may have attained maturity. Source rock material of organic clays is deposited in an anaerobic environment, and peaty formation forms residual coal deposits and source material for petroleum and natural
Figure 2. Representing gradient of bouger anomaly gravity, isostatic anomalies; magnetic depth and basement element map.
gas. The reservoir rocks, which are deposited are porous and permeable rocks and are carbonates and sandstones. Some carbonates reservoirs owe their porosity and permeability to digenetic and tectonic modification of the original sediments and are interceded with potential cap rocks. For the crustal modeling, the following parameters were used: the average density distribution of rocks in the Lower Indus Platform basin, which is about 2.8 g/cm and the crustal velocity, taken as 6.4 k/m [14]. The mean mantle velocity is about 8.1 k/m as by [15]. The Isostatic behavior of the area has been studied by calculating Isostatic Mohofor normal crust of 30 Km. [16]. The interpreting anomalies which indicate the presence of basic rocks and various rocks types having structural features such as faults and fold have been traced. Due to lacking of integration with corresponding seismic and well information for this study, qualitative and quantitative interpretation of the gravity and magnetic data was performed in conjunction with regional reflection profiles and exploratory drill hole information. The results of this synergetic approach show a provisional correlation of the potential field with reflection data exposing the tectonic-structural interpretation.
5. Seismic Data and Integrated Interpretation
Our analysis shows that most of the misinterpretations suggested above can be avoided by incorporating seismic information, particularly for isolating the real sources of local magnetic and gravity anomalies. The acquired seismic data base consists of regional and detail seismic of various vintages. Migration was also applied to the acquired data and then processed. These lines constitute the main grids from which the refined structure and paleotectonic evolution of the area were deciphered. The time-to-depth conversion of available seismic data using Layer Cake method after marking good reflectors assigns ages on the basis of well and long range correlation and velocity analysis which can be comparable to [13,17]. The [18] interval velocities were computed for each interval (layer) in the sedimentary wedge above the basement using Gardner relationship [19], and bulk densities were computed as given in Table 2. The sedimentary structure was interpreted from the seismic reflection data, the density contrasts used for each seismic sequence have been taken against the crust of mean density 2.8 gm/cm3. There is a complex structural deformation of the Jurassic starta and many different fault systems. Faults are represented on the seismic sections (Figures 3 and 4) as a discontinuous reflection along a preferred orientation of reflectors. The seismic sections that belong to the Lower Indus platform basin are represented in Figures 3 and 4. The interpretation of seismic lines demonstrated the presence of two main seismic events, viz. “syn-rift” and
Table 2. Average densities for different lithologies.
Figure 3. Representing the seismic data from Lower Indus Platform basin area.
Figure 4. Representing the seismic data from Lower Indus Platform basin area.
post-rift. The Late Miocene regression in the area caused large scale channeling followed by the growth fault related structures, evident on the seismic section [20]. The subsurface structures are related to normal faults in the east and listric normal faults and diapiric structures in the west [21]. The divergent pattern of seismic reflections from the strata within the basin suggests that the processes of continental rifting and deposition were contemporaneous. The younger sequence seems to have been affected by minor reactivation of an older fault. The fault bounding the basin is regarded as normal and listric in form indicating extension and accompanied with rotation of the basement block. The cap rocks are present for all potential reservoirs such as intra-formational shale for Lower Cretaceous reservoirs. The reservoirs are Cretaceous to Eocene clastics and carbonates, proven due to the presence of rich hydrocarbon discoveries in the Lower Indus platform basin. The Cretaceous rocks are potential source rocks for hydrocarbon in the area. The reservoirs which are older are present towards the east and northeast. The mature source rocks which are relatively younger are present towards the west. The structural evolution of a Lower Indus Basin is therefore controlled mainly by the amount of initial crustal thinning that occurs during the rifting phase as well as subsequent tectonic subsidence. The tectonic subsidence is a result of the cooling of thinned lithosphere (thermal contraction). The process of thinning produces a thermal anomaly due to passive upwelling of hot asthenosphere. The regional seismic lines clearly show very thick sediments in the area which delineate the structural configuration. The quality of existing seismic has allowed for conclusive interpretation. However, combined information from seismic, potential field, and well data enables identification of potential fields, while well data enables the identification of major structural and tectonic elements, estimation of their size and geometry, and understanding of their spatial and temporal relationship.
6. Conclusions
We have used the integrated geophysical approach in this study of the Lower Indus Platform basin which has proven to be efficient using reflection seismic profiles and potential field maps.
We have concluded that the trapping mechanis midentified from our interpretation indicates that Lower Indus Platform basin reservoir rocks have tremendous reserves of hydrocarbons and the tectonics modification of the basin has formed variety of trap structures. We have also concluded that the reservoir rocks which are deposited are porous and permeable rocks and are present in carbonates and sandstones, and the sources rocks are allshaley in nature. The reflected seismic profiles and potential field maps can be integrated into the practice for hydrocarbon exploration. We hope that our research study will be very useful for the exploration and production activities to explore and to enhance the energy deficiency and to fulfill needs requirements for the economy of the country and also can be used for anyone who interested in research geophysics. We suggest that there is need to do the research on the basis of integrating well logging data with seismic data by using neural network technology.
NOTES