Journal of Geology & Geophysics

Journal of Geology & Geophysics
Open Access

ISSN: 2381-8719

+44 1478 350008

Research Article - (2019) Volume 8, Issue 1

Integrated Seismic Interpretation, 2D Modeling along with Petrophysical and Seismic Attribute Analysis to Decipher the Hydrocarbon Potential of Missakeswal Area, Pakistan

Ayesha Abbas1*, Hongtao Zhu1, Aqsa Anees1, Umar Ashraf2 and Nosheen Akhtar3
1Faculty of Earth Resources, China University of Geosciences, Wuhan 430074, China
2Institute of Geodynamics and Geophysics, China University of Geosciences, Wuhan, China
3Department of Earthscience, Quaid-i-Azam University, Islamabad, Pakistan
*Corresponding Author: Ayesha Abbas, Faculty of Earth Resources, China University of Geosciences, Wuhan 430074, China, Tel: +86-18502714190 Email:

Abstract

Missakeswal is an important hydrocarbon field, lying on active foreland fold and thrust belt of Himalayan orogeny in Potwar plateau. Integrated study of 2D seismic data (SEG-Y, Navigation and seismic velocities) and well logs helps us to delineate the potential reservoir rock of the area. Seismic interpretation based on stratigraphic studies and well tops, aids to mark four reflectors; Chorgali, Sakesar, Lockhart and Basement. Time sections are converted to depth section using velocity analysis system to delineate subsurface structure. Besides this, fault-bounded anticlines and crustal shortening analysis of the depth sections, revealed that folding in the sedimentary successions pre-date reverse faulting and regime of the Potwar basin, is suitable for hydrocarbon accumulation. 2D modeling of the interpreted seismic sections confirms reverse faulting in the sedimentary successions and normal faulting in the basement. Moreover, Seismic Attributes Analysis has carried out which helps in understanding the lateral continuity, bedding sequences and thickness of desired beds highlighted the petroleum system and affirmed the interpretation. The identified structural variations would help in the understanding of the regional tectonic settings, besides this, reservoir character in terms of lateral thickness variation, fault offsets and lithological dissimilarities are achieved. It also reveals that carbonate successions of the Sakesar and Chorgali formations acted as potential reservoirs in Missakeswal area.

Keywords: Interpretation; Velocity; Seismic attributes; Modeling; Reservoir

Introduction

Missakeswal oil field is located 60 km SE of Islamabad, lies on the Eastern extremity of the Potwar sub-basin, Punjab Province, Pakistan (Figure 1). The estimated terrain elevation of area is 500 m above mean sea level [1]. Potwar basin is an important host to hydrocarbon reserves, lies at the northern side of the Upper Indus basin [2]. Geographically study area bordered by Main Boundary thrust in North, Jhelum strike slip fault in East, Kala bagh strike slip fault lies towards West and Salt range in south. Pakistan has high potential of hydrocarbons in its northern (Potwar, Kohat) and southern (Badin, Mari etc.) parts. The topography Potwar plateau is undulating comprised of E-W trending low hills, parallel ridges and dissected valleys. Chorgali and Sakesar Formations predominantly comprised of limestone of Eocene age are the well-known oil producing reservoirs [1]. Due to presence of detachments at different levels, complex structures, anomalous evaporates thickness and technical drilling problems; so far most of the wells were unsuccessful in this area [3-5]. Two-dimensional seismic reflection data reveals the presence of normal faults in basement, and reverse faults in the cover sequence, which can be seen on the interpreted sections [5-7].

geology-geosciences-Missakesal-area

Figure 1: Satellite image of Pakistan showing Missakesal area (Khan et al., 2008).

This study is intended to evaluate the hydrocarbon potential by using seismic reflection, well log and velocity data, will be helpful in future for prospects generation and hydrocarbon exploration in Potwar Basin. The average, interval and RMS (Root Mean Square) velocities will helps in the study of subsurface lateral variations, time to depth conversion and their correlation with synthetic seismogram on depth section affirms interpretation. Exploration companies integrate geophysical, geological and well data to carry out reliable subsurface structural interpretation. Precise description of the subsurface geological horizons (Chorgali, Nammal and Sakesar), structures and seismic attributes provide an important insight to hydrocarbon exploration. The basement with normal faults and cover sequence with fault bounded anticlines displaying Jurassic rifting, Pangea splitting and associated diapirism of Himalayan orogeny in the Potwar plateau. The main discoveries of the Misskeswal area can be related with salt cored anticline, which are formed because of crustal shortening in the cover sequences. The calculated average energy and instantons phase attributes helps in study of reflection continuity, faults, lithological dissimilarities and data type. Similar approach used in the proposed study will also aid in delineation and exploration of regional subsurface structures in the Potwar basin and nearby areas.

Geological settings

The Himalayan collision system represents an active collision Orogen between Indian and Eurasian subcontinents. The collision is active since about 55 Ma [8] involves continuous uplifting erosion and deposition of the sediments. The Indus basin, including the Kohat- Potwar (study area) depression, belongs to the category of extracontinental downward basins, which account for 48% of the world petroleum resources [9]. The tectonic depression of Potwar sub-basin formed as result of continent-to-continent collision at NW margin of the Indian plate. The foreland zone; comprised of Salt range, Kohat and Potwar plateau and Hazara Ranges, the area is bounded by Panjal- Khairabad fault in North Salt Range thrust in South, whereas eastern end up with NS trending left lateral Jhelum fault [10]. Potwar sub-basin is predominated by structural traps and mostly seismic data used for the delineation of these structures. Missakeswal area is well known for its hydrocarbon potential, some of the exploration trends and oil field locations of the Potwar plateau is shown in Figure 2.

geology-geosciences-Geological-map

Figure 2: Geological map of Potwar plateau, oil fields and wells are higlighted (William et al., 1998).

Potwar plateau comprised of 5000-1000 m thick sediments of Precambrian to recent age, these sediments range from fluvial to marine environment of deposition [11,12]. The basement of Precambrian is mainly dominated of slates and phyllites. The Salt Range formation unconfirmably overlying the basement, composed of metamorphic and Volcanic rocks of Indian Shield [13], and overlaind unconfirmably with Khewra Sandstone of Early Cambrian. Contrary to the well data obtained from Adhi and Missakeswal Kussak (sandstone and shale), Jutana (dolomite) of early Cambrian and middle Cambrian Baghanwala (shale, sandstone, siltstone) lies conformably over Khewra sandstone [5]. The Nilawan group comprising of Dandot (sandstone) and early Permian Tobra (conglomerates), whereas Warcha and Sardhai formations are eroded in study area.

The unconfirmable boundary exposed the area from upper Permian carbonates to lower Paleocene [14]. Hangu (sandstone), Lockhart (limestone) and Patala (shales) are well developed Paleocene sequence and Nammal, Sakesar and Chorgali Formations of Eocene age conformably overlies Paleocene strata. The Mesozoic sediments are eroded and Oligocene sediments have not been found in the region [2]. Stratigraphy of the Potwar plateau is well developed along the Salt range; a generalized stratigraphic chart is given in Table 1.

Group Formation Lithology and Thickness (m) Age
CHARRET Chor Gali Limesrone (30-140) m EOCENE
Sakesar Limesrone (70-300) m EOCENE
Nammal Shales, marl lirnesionc (35-130) m EOCENE
MAKARWAL Patla Shales, marl, sandstone (27-200) m PALEOCENE
Lockhart Limestone (30-240) m PALEOCENE
Hangu Sandstone,shales (15-150) m PALEOCENE
  Samanna Suk Limestone, marl, shale (66-366) m PALEOCENE
Datta Sandstone (30-400) m J URASSIC
M USAKHAL Kingriali Limestone, dolomite (100) m TRIASSIC
ZALUCH Chhidru Sandstone, limestone (64) m PERMIAN
Wargal Limestone, dolomite (183) m PERMIAN
Amb Sandstone, limestone, shales (80) m PERMIAN
NILAWAN Sardhai Clay, sand, shales (42) m PERMIAN
Wareha Sandstone, shales (26-180) m PERMIAN
Dandot Sandstone, shales (50) m PERMIAN
Tobra Siltstone, shales (JJ) m PERMIAN
JHELUM Baghanwala Shales. Sandstone (40 )m CAM BRIAN
Jutana Dolomite (80 )m CAM BRIAN
Kussak Sandstone, siltstone,dolomite (6-53) m CAMBRIAN
Khewra Sandstone,shale ( 150) m CAMBRIAN
  Salt Range Marl, salt, clay, shale, dolomite (830) m PRECAMBRIAN

Table 1: Generalized stratigraphic Chart of the study area.

Potwar marine facies has great potential of hydrocarbon generation. Recent discoveries in Potwar result in delineation of deep subsurface crest [4]. The gray shales of the Datta (Jurassaic age) and Patala Formations (Paleocene age) are potential source rocks in Salt Range Potwar- Foreland Basin (SRPFB) [1]. The Cambrian, Permian, Jurassic, Paleocene and Eocene reservoirs are producing oil in Salt Range Potwar-Foreland Basin (SRPFB). Thick layers of evaporite and shale have good sealing potential for Infra-Cambrian reservoir. Most of the fields discovered in Kohat-Potwar geological province are overturned faulted anticlines, popup structures or fault-block traps. The latest trapforming thrust event began at approximately 5 to 2 Ma [15].

Presently two-fold division is envisaged for the Potwar sub-basin. Based on the deformation style; (1) the Northern Potwar Deformation Zone (NPDZ) and (2) the southern Potwar platform zone [1,4,8].

The Northern part of Potwar is more intensely deformed, known as the North Potwar Deformed Zone (NPDZ), characterized by East West, tight and complex folds, overturned to the south and sheared by steep-angle faults [10]. The central Potwar comprised of fault bounded structures mainly thrust and back thrusts, while asymmetric anticlines with single faults have been observed in different locations. Main structure in the Potwar area based on seismic interpretation is Salt cored anticlines bounded with faults on limbs and Triangle zones [16].

Work Flow and Methodology

Seismic, well log, Velocity and Navigation data of fifteen seismic lines were acquired from DGPC (Directorate General of Petroleum concession), Government of Pakistan. Missa-keswal area is projected in Universal Transverse Mercator (UTM, Zone 43) projection system, Base map generated using K-tron Precision Matrix is shown in Figure 3.

geology-geosciences-Base-map

Figure 3: Base map of the study area displaying wells and seismic lines (Khan, 2000)

Different modules of K-tron software used to accomplish various processes; workflow to carry out this research work is shown in Figure 4. Interpretation has been carried out by using X-works software, provides an interface for marking horizons and faults, velocity and depth data for contouring and for further analysis such as crustal shortening and 2D seismic modelling. VAS (Velocity Analysis System) was used to calibrate velocity data which is used in time to depth conversion [17]. Crustal shortening and 2D modeling performed to confirm compressional regime and structural interpretation. The seismic velocities used in different applications like time to depth conversion, contouring and Modeling. Using depth section and velocities, we can also depict structures, lithological boundaries and the lateral variation in the subsurface area. Seismic attributes analysis affirms interpretation along with that also provide information about bedding continuity, lithology and helps us to predict the probable reservoir of the area.

geology-geosciences-Work-Flow

Figure 4: Integrated Work Flow used in study.

Results and Discussion

Seismic interpretation

Initially different reflector was marked based on prominent reflections; well data used as secondary information from subsurface depths utilized to pick prominent reflections. Each reflector is marked with different color; so, that they can easily distinguished, main reflectors of on the sections are Chorgali, Lockhart and Sakesar Formations (Figure 5). Processing of Seismic velocities (conversion, interpolation and smoothing) carried out by using X-works software, which then used in different applications like time to depth conversion, contouring and Modelling [18]. Depth sections give real picture of the subsurface and the horizontal variation in the horizon, which occur due to the lateral variation of the velocity, because only time section is insufficient to delineate the results (Figure 6). Average velocities are multiplied with time of each reflector and the divided by two to acquire depth section (Figure 7). Horizon interpolated average velocity contour map of Chorgali Formation generated to understand the spatial velocity variations provide true depth imaging of the formation (Figure 8).

geology-geosciences-seismic-line

Figure 5: Interpreted time section of seismic line GNA-19.

geology-geosciences-Seismic-velocities

Figure 6: Seismic velocities of GNA-19 combinedly displayed at different CDP’s.

geology-geosciences-Depth-Section

Figure 7: Depth Section of seismic line GNA-19.

geology-geosciences-contour-map

Figure 8: Chorgali Formation Horizon Velocity contour map.

Crustal shortening

When two continental crusts collide and neither sub ducts, the material pushed towards the Earth’s surface resulting in the formation of mountains like Mount Everest. This causes the crust to reduce in size [19]. The crustal shortening of all the three formations of seismic line GNA-14, is determined by loading the depth section in X-Works that gives the length of each segment of horizon (Figure 9). The horizon divided into different segments because of reverse faulting. The crustal shortening calculated by using the formula given below;

geology-geosciences-cross-Section

Figure 9: Geological cross-Section for crustal shortening with length of each faulted part.

CS = L – Lo (1)

L is the total length of the horizon obtained by adding all the segments of the horizon. Lo the length of the seismic section calculated by multiplying the CDP (Common Depth Point) number with the CDP interval. Crustal shortening performed to confirm compressional regime and structural interpretation. The horizon length and their corresponding crustal extension are given in Table 2.

Horizon Name Horizon Length L(m) Section length Lo (m) Crustal Shortening L-Lo (m)
Chorgali 22359 21400 959
Sakesar 22343 21400 943
Lockhart 22336 21400 936

Table 2: The horizon lengths and their corresponding crustal shortening.

Volume of shale

Missa-01 well illustrates reservoir and source zones, the prominent shale signature is corresponding to the location of the Patala shales at depth 1750 and 2250 M as given by the information from the well tops, while the zone which can be seen above shale is reservoir at depth of about 1800-2020 m, which corresponds approximately to the depth of Chorgali and Sakesar Formation. The response of Gamma ray log along with clear description of shale and clean line is shown in Figure 10.

geology-geosciences-Gamma-Ray

Figure 10: Gamma Ray log (red) with calculated Volume of Shale curve (blue). [Cut-off line (sky blue), shale line (green) and sand line (yellow)].

Seismic modeling

1D- forward modeling: An artificial seismic reflection record, prepared from velocity log data by convolving the product of acoustic and density logs (reflectivity function) with seismic data derived wavelet is termed as Synthetic Seismogram [20]. Synthetic seismogram of Missa-02 well generated by using Wavelets software [21].

In this procedure, the petrophysical logs; Sonic (DT) and Bulk Density (RHOB), which provide the inverse velocity and density information of subsurface layers respectively, this velocity and density data used to compute a series of reflection coefficients called reflectivity series. Source Ricker wavelet with a dominant frequency of 35 Hz is generated. The reflectivity series then convolved with the source wavelet shown in Figure 11. Synthetic seismogram matched with the interpreted seismic depth section at the well point to correlate the succession of reflectors as shown in Figure 12. It also used to calibrate our seismic velocities. 1D modeling is important to study the area at different points and it does not give information of lateral variation that is computed by using 2D modeling.

geology-geosciences-Seismic-Modeling

Figure 11: A Ricker Wavelet Generated with frequency 35 Hz used as a source for Seismic Modeling.

geology-geosciences-Geological-Depth

Figure 12: Geological Depth Section of GNA-19 Correlated with Synthetic Seismogram generated from sonic log.

Depth section based seismic model (2D Seismic Modeling)

The reverse process of seismic interpretation is now as seismic modeling, interpretation is carried out by marking the horizons based on continuity of reflector on the Seismic section whereas, in 2D seismic modelling the synthetic seismic section, is generated by the convolution of Ricker wavelet of 35 Hz with the Geological section of seismic line, here we are using GNA-19 at 25 trace intervals as shown in Figure 11.

The time sections are converted to depth sections using velocity data and reflection coefficient were assigned to each horizon as an attribute. The geological section is convolved Ricker wavelet at every 10th CDP interval to generate a 2D synthetic seismic section as in Figure 13. This modeling is based on structural data, so the generated synthetic section corresponds to the subsurface structures. 2D Modeling is important in the areas or points where they arise lateral variation that may give rise to focusing problem, shadow zone, diffractions, and beneficial for adjustment of acquisition parameters foe neo surveys and seismic interpretation affirmation etc. [22].

geology-geosciences-cross-section

Figure 13: Geological cross-section with generated Synthetic Seismic Model.

Seismic attribute analysis

Seismic attributes are widely used for lithological and petro physical prediction of reservoirs. Various methodologies have developed for their application to broader hydrocarbon exploration and development decision-making [23]. There are no specific rules for computing attributes, a quantity that is directly calculated from seismic data is called seismic attribute [24]. Seismic attributes affirm interpretation along with that also provide information about bedding continuity, porosity, phase, and lithology and helps us to predict the probable reservoir location [25].

Envelope of trace (reflection strength) and average energy/ apparent polarity

Envelop as a measure of local energy level of seismic traces has related to the interval thickness of the corresponding geological unit.

Trace envelop helps in horizon and fault interpretation particularly smaller faults, to recognize the major changes in depositional environment and provide aid in Spatial correlation to porosity and other lithological variations [26].

Envelop trace is computed for seismic line GNA-19 to see the major changes in lithologies. The maximum reflection strength indicated by the thick yellow packages, and it also showing spatial patterns representing changes in the limestone thickness and breakage due to the faults (Figure 14). Weak reflections are due to the interbedded shales [27,28]. Even negative reflection coefficients such as limestone formation overlaid on clayey formation would generate a positive response in this attribute.

geology-geosciences-Envelop-attribute

Figure 14: An Envelop attribute map of GNA-19 seismic line [Chorgali (Red), Sakesar (Green), Lockhart (Blue)].

Average energy is a post-stack wavelet attribute, which is computed at the peak of envelop attribute within a specified window, in a defined zone the attribute of wavelets are represented by trace envelope minima. Thus, wavelet attributes are computed over each cycle of the seismic trace [29]. Apparent polarity shows the blocky nature as it represents the average value over each cycle in the seismic trace, which can be observed in Figure 15. Each peak in the reflection indicates the cycle. The Polarity base attribute can have used for better understanding in interpretation particularly for fault interpretation [30].

geology-geosciences-attribute-map

Figure 15: Average energy attribute map of GNA-19 [Chorgali (Red), Sakesar (Green), Lockhart (Blue)].

Instantaneous phase

The information of phase is related with the propagation of seismic wave front, and independent of trace amplitude. It computed from real and imaginary traces as given below;

Ө (t) = tan-1[h(t)/f(t)] (2)

Instantaneous phase attribute used to highlight the deeper horizons that are not visible in the normal amplitude section as shown in the Figure 16. The interpreted horizons lie over the minimum phase regions indicated by orange color, instantaneous phase attributes which changes from -1800 to +1800 [31]. This attribute further confirms the interpretation as the input data is minimum phase. The negative value in color bar is indicating negative phase. It can be observed in comparison to amplitude-based sections that the instantaneous phase shows much deeper horizons. Phase attribute is sensitive to lithological boundaries and good indicator of lateral continuity, so this can be used for marking boundaries in interpretation.

geology-geosciences-Average-energy

Figure 16: Average energy attribute map of GNA-19 [Chorgali (Red), Sakesar (Green), Lockhart (Blue)].

Conclusions and Perspectives

• Interpreted seismic reflection data of Missakeswal area indicates the presence of normal faults in the basement and salt cored anticlines bounded by reverse faults on both limbs in the cover sequence.

• The study area had experienced crustal shortening, as result of continent-continent collision at NW margin of Indian plate and the tectonic forces acting on North Potwar Deformed Zone (NPDZ) affirmed the compressional tectonic regime in the study area.

• ISO (International Organization for Standardization) Velocity contour maps indicate the velocities vary laterally and horizontally due to compaction, lithological contrast, faulting and different burial depth.

• Synthetic Seismogram matches with the marked reflector, which confirms stratigraphic interpretation, and interpreted 2D model of geological cross-section confirms the type and nature of sub-surface structures.

• The structural and lithological interpretation is further affirmed through seismic attributes maps i.e. envelope, instantaneous phase and average energy. Instantaneous phase attribute provided good results and the precise positions of faults. The overall results indicate the economic viability of Chorgali limestone as a reservoir.

Acknowledgements

I am thankful to the Director General, Petroleum Concessions (DGPC), Pakistan, for the provision of data for this research work and to Sir Khalid Amin Khan (OIST) for his thoughtful guidance and technical manner, I found it quite helpful for me.

References

  1. Khan AM, Ahmed R, Raza HA, Kemal A (1986) Geology of petroleum in Kohat- Potwar depression, Pakistan. Am Ass Petro Geol Bull 70: 396-414.
  2. Amir M, Siddiqui MM (2006) Interpretation and visualization of thrust sheets in triangular zone in eastern Potwar, Pakistan. The Leading Edge 25: 24-37.
  3. Jaume SC, Lillie RJ (1998) Mechanics of the salt range-Potwar Plateau, Pakistan: A fold-and-thrust belt underlain by evaporites. Tectonics 7: 57-71.
  4. Moghal M, Anwar (2007) Subsurface geometry of Potwar sub basin in relation to structuration and entrapment. Pak J Hydroc Res 17: 61-72.
  5. Lillie RJ, Johnson GH, Yousaf M, Zamin ASH, Yeats RS (1987) Structral development within the Himalayan fireland and fold thrust belt of Pakistan. Can Soc Petro Geol Mem 12: 379-392.
  6. Pennock ES, Lillie RJ, Zaman ASH, Yousaf M (1989) Structral interpetation of seismic reflection data from eastern salt range and Potwar Plateau, Pakistan. Am Ass Petro Geol Bull 73: 841-857.
  7. Jadoon IAK, Frisch W, Jaswal TM, Kemal A (1999) Triangle zone in the Himalayan foreland, north Pakistan. Geol Soc Am Spec Pap 328: 277-286.
  8. Hasnay ST, Saleem U (2001) An integrated subsurface, geological & engineering study of Meyal Field, Potwar Plateau, Pakistan. AAPG Search and Discovery Article, PAPG/SPE Ann Tech Conf.
  9. Kazmi AH, Jan MQ (1997) Geology and tectonics of Pakistan, Graphic publishers, Karachi, Pakistan.
  10. Gee ER (1980) Pakistan geological salt range series. Directorate of overseas surveys, United Kingdom, for Gov Pak Geo Surv Pak 6 sheets, scale 1,50,000.
  11. Kazmi AH, Abbasi AI (2008) Stratigraphy and historical geology of Pakistan. Dept. & Nat Cen Excel Geol, Peshawar Publ, 181-185.
  12. Yeats RS, Lawrence RD (1984) Tectonics of the Himalayan thrust belt in the northern Pakistan. US-Pakistan Workshop on Marine Sciences in Pakistan 11: 39.
  13. Shah SMI (2009) Stratigraphy of Pakistan. Government of  Pakistan Ministry of Petroleum and Natural Resources. Nat Resour Geol sur Pak 245-273.
  14. Jaswal TM, Robert JL, Lawernce RD (1997) Structure and evolution of the northern Potwar deformed zone, Pakistan. Am Ass Petro Geol Bul 81: 308-328.
  15. Mughal MA, Hameed A, Saqi MI, Bugti MN (2003) Subsurface geometry of Potwar sub-basin in relation to structuration and entrapment. ATC Conference & Oil Show Islamabad.
  16. Khan KA (2000) Integrated Geo-Systems- A computational environment for integrated management, analysis and presentation of petroleum industry data. Am Ass Petro Geol 215-226 pp.
  17. Khan KA, Akhter G (2011) Work flow shown to develop useful seismic velocity models. Oil & Gas J 109: 52-61.
  18. Van Hise CR (1898) Estimates and causes of crustal shortening. J Geol 6: 29-31.
  19. Peterson RA, Fillipone WR, Coker FB (1955) The synthesis of seismograms from well log data. Geophys 20: 516-538.
  20. Khan KA, Akhter G, Ahmed Z, Khan MA, Naveed A (2006) Wavelets - A computer based training tool for seismic signal processing, Pakistan. J Hydrocarbon Res 16: 37-43.
  21. Gadallah J, Fisher I (2009) A comprehensive introduction to Exploration Geophysics. Exploration geophysics, Springer.
  22. Chopra S, Marfurt K (2006) Seismic attributes-A promising aid for geologic prediction. Can Soc Exp Geophys 31.
  23. Coffeen JA (1986) Seismic exploration fundamentals. Pennwell Publishing Company, Tulsa, Oklahoma.
  24. Gee ER (1989) Overview of the geology and structures of the Salt Range, with observations on related areas of northern Pakistan. Geol Soc Am Spec Pap 232: 95-111.
  25. Khan KA, Akhter G, Ahmad Z, Rashid M (2008) Development of a projection independent multi-resolution imagery tiles architecture for compiling an image database of Pakistan. Proc 2nd Int Conf Adv Space Technol Islamabad, Pakistan, 164-170.
  26. Khan KA, Akhter G, Ahmad Z (2009) The Real meaning of geoscience data and process integration. Proc IAMG Int Conf.
  27. Khan KA (2010) Seismic the next step series: Seismic attributes. OIST, Islamabad.
  28. Taner MT (2001) Seismic attributes. Rock solid images. Can Soc Exp Geophys 26.
  29. William JS (1998) Wrench faulting in the northern Pakistan foreland. Am Ass Petro Geol 82: 2003-2030.
Citation: Abbas A, Zhu H, Anees A, Ashraf U, Akhtar N (2019) Integrated Seismic Interpretation, 2D Modeling along with Petrophysical and Seismic Attribute Analysis to Decipher the Hydrocarbon Potential of Missakeswal Area, Pakistan. J Geol Geophys 7: 455.

Copyright: © 2019 Abbas A, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Top