Evaluation of Neoproterozoic Source Rock Potential in SE Pakistan and Adjacent Bikaner- Nagaur Basin India Based on Integrated Geochemical, Geological, and Geophysical Data


 Discoveries of heavy crude oil in the Neoproterozoic rocks (Infracambrian rock sequence) from the Bikaner-Nagaur Basin of India emphasizes the significance to study and explore the Neoproterozoic source rocks potential in the southeastern part of Pakistan. This study evaluates the potential of the source rock in the Infracambrian rock sequence (Salt Range Formation) based on surface geochemical surveys, Rock-Eval pyrolysis, source biomarkers, geophysical characterization, and seismic inversion using machine learning for maturity index estimation. Core samples of Infracambrian rock were extracted for Rock-Eval pyrolysis and biomarker characterization. Also, 81 geo-microbial soil and gas samples were collected from the surface to explore the petroleum system and potential source rocks in the subsurface. We followed the standard laboratory procedures to investigate the origin and concentration of hydrocarbons gases at the surface, thermal maturity, the source facies, and the environment of deposition of organic matter. The results show that the investigated samples are characterized by restricted marine clay devoid of carbonate source facies with thermal maturity in the early-stage of the oil generation window. Surface geochemical samples also confirm higher concentrations of thermogenic C2-C4 hydrocarbons over the vicinity of anticlinal structures proving the existence of an effective migration path along deep-seated faults to the surface. The inverted maturity index profile demonstrates a reasonable correlation of thermal maturity with the surface geochemical survey, source biomarkers, and Rock-Eval pyrolysis. It validates the reliability of multilayer linear calculator and particle swarm optimization algorithms for inverting seismic reflection data into a maturity index profile. The obtained results indicate a higher probability of heavy and light oil along the eastern flank of Pakistan, where Infracambrian rocks are thicker and more thermally mature, and deep-seated pledged structural closures occur, in comparison to the Bikaner-Nagaur Basin, India.


Introduction
Neoproterozoic rocks are among the oldest rock units in the world, and several successful discoveries have been reported from these sequences in Salt Basins of Oman and Siberia [1][2][3][4][5] . In the Bikaner-Nagaur (BK-N) Basin India, 935 million barrels of heavy oil reserves have been discovered in the Infracambrian Jodhpur sandstone, the Bilara dolostone, and the thin siltstone layers in the Hanseran evaporite which the porosity ranges between 16-25%, 7-15%, and 3-12%, respectively 5,6 . The heavy oil discovered is nonbiodegraded, and biomarker ratios support crude oil from the source rock in the early-stage of oil generation 2,7,8 . In the BK-N Basin, the biomarker study from the Baghewala-1 well reveals that mature hydrocarbons have been generated from a deeper buried source rock 2 . Several authors reported that the same quantities of petroleum could be generated from the Infracambrian source rock of Salt Range (SR) Formation in southeastern (SE) Pakistan, where Infracambrian source rock beds are much deeper and thicker than that in the BK-N Basin 2,6,9 . In recent years, the Infracambrian play in SE Pakistan has been explored to assess the hydrocarbon potential by geological modelling or seismic data characterization 6,9−12 . However, the previous attempts The Infracambrian SR Formation is exposed in the northeastern part of Pakistan (Potwar Plateau) 9,10 .
Exposures of Precambrian rocks occur as many hillocks (volcanic and volcaniclastics) located at Kirana Hills (Hachi volcanic) in the Sargodha-Shahkot-Sangla-Chiniot region of SE Pakistan 15,16 . These Hachi volcanic are analogous to Infracambrian rocks in the Nagar Parkar (Tosham-Malani volcanic India), as these crystalline rocks, with or without the cover of Phanerozoic, lie at the bottom of the Indus Basin 9 .

Lithostratigraphy
The lithostratigraphy of the study area was investigated using outcrops and the well logs data from BK-N Basin India and SE Pakistan. The lithostratigraphy of the BK-N Basin is dominantly composed of nearly 1,000 m of shales, sandstones, carbonates, and evaporates. This lithology is quite similar to the Potwar sub-basin (upper Indus Basin) 6,10 . Several regional unconformities, erosion, and hiatuses were observed in the study area through the sedimentary succession. The Infracambrian Salt Range (SR) Formation unconformably overlies a Precambrian basement 18 . Another erosional unconformity surface is observed between the Early-Middle Cambrian SR Formation and the Middle/Upper Permian Tobra Formation (Fig. 2). This widespread Early-Middle Cambrian and the Middle/Upper Permian unconformity occurred due to localized uplift and a global eustatic fall in sea level. The upper part of the SR Formation is sometimes absent due to erosion during the Cambrian and Permian unconformity. The maximum thickness of the SR Formation is 906 m.
The Permian rocks are covered by about 400 m of thick marine Cenozoic sediments, gradually thickening westward. The Paleocene-Eocene succession comprises foraminiferal sandstone, limestone, and shale. It is overlain by Neogene uvial clastic consisting of thick alluvium silt, sand, and clay layers 13 . Fig. 2 shows a generalized stratigraphic column based on data from the Karampur-1 well.
Four stratigraphic units of the SR Formation are identi ed in the study area that are equivalent to the Infracambrian rocks in the BK-N Basin. The lithostratigraphy of the SR Formation is based on the depositional environment and lithological log information from Suji-1, Marot-1, Karampur-1, and Bahawalpur East-1 wells 10,15,19 . The SR Formation is subdivided, from oldest to youngest, into four formations: Sonia, Jodhpur, Bilara, and Hanseran (Fig. 3). The Hachi volcanic (equivalent to Tosham Malani igneous suite) forms the basement. Both Hachi and Tosham Malani volcanic unconformably underlie the Basal conglomerates, composed of dolomites and volcano-clastic. The clastic sediments of the Jodhpur Formation form the primary potential reservoir in the basin and can be classi ed into Sonia shale and sandstone members. The Bilara-Hanseran sequence lies conformably over the Jodhpur Formation. The Bilara-Hanseran acts as a source rock with TOC ranges from 5 to 6% 20 . The calcareous Bilara Group comprises dolomite, limestone, stromatolitic limestone, and occasional clay beds. The Hanseran Evaporite Group comprises mainly of halite, anhydrite, dolomite, and reddish sandstone with thin shale beds. The Bilara-Hanseran sequence has total estimated thickness of about 200 m 20 . The comparison of Infracambrian and Cambrian succession between the BK-N Basin and SE Pakistan is described in Table 1.  Figure 5 shows the Infracambrian play concept in Kirana-Malani-Tsham Basin 6 . From rifting to collision and underthrusting of the Indian Plate during Infracambrian to Mesozoic, the basin is associated with deep-seated normal faults, salt-induced structures, horst and graben structural highs and lows. The previous drilling campaign in the BK-N Basin revealed that wells located on paleo-highs did not encounter Infracambrian strata than those drilled on paleo-lows 9 .

Data description
Logs and core data, including well-cuttings from nine wells, as shown in Table 2, were used for estimating the total organic carbon (TOC), rock-eval pyrolysis, biomarkers, mineralogy, saturation, porosity, and permeability in the Infracambrian to Permian rocks. Also, surface geo-microbial soil and gas samples were investigated to evaluate the petroleum potential of source rocks in the subsurface. Moreover, 1000 km-long 2D seismic data were interpreted to mark the stratigraphic horizons and structural style, focusing mainly on the top and bottom re ectors of the Infracambrian-Cambrian rocks.

Methods
The following sequence was adopted to achieve the study's objectives.

Surface geochemical survey
Surface geo-microbial soil and soil gases were studied to explore the prospective areas of hydrocarbon microseepage and its origin. For this purpose, we collected 81 soil and gas samples along seismic lines over the study area, as shown in Fig. 6. The samples were extracted from varying depths, 0.8 to 1.5 m using a soil sampler, mechanical auger, sterilized soil sampling bags, and vacuum gas sampling bags. To suck up the gas, a syringe was injected into the septum after sealing the annulus between the sampling tube and the borehole, as described by Tedesco 22 . While the area of interest is situated among arid climate zone, all necessary measures were taken into consideration during sampling, e.g., avoid sampling the soils from water-saturated as well as excavated areas, land contaminated by hydrocarbons, animal feces, the ow of chemicals, marshes, and the areas under water table 23 .
We incubate soil samples with light hydrocarbon gases in the mineral salt medium for geo-microbial analysis and count the hydrocarbon oxidizing bacteria that have developed. For gas analysis, hydrocarbon gases desorbed from soil samples through acid-base extraction were analyzed by gas chromatography. The results were shown using concentration distribution maps to identify areas with anomalously high concentrations. Anomalies were marked through a statistical analysis approach, where the mean plus half standard deviation was assumed as the background value.

Rock-Eval pyrolysis
Rock-Eval pyrolysis was used to identify the maturity of OM and generative potential of the Infracambrian rock. Four basic parameters, such as oxygen index (OI), production index (PI), hydrogen index (HI), and maximum temperature (Tmax), were obtained by pyrolysis. The chosen samples were measured using Vinci Rock-Eval 6 apparatus by following the previously reported procedures [24][25][26] . Approximately ten well cuttings samples from the Bijnot-1 were investigated to measure the OI, HI, PI, Tmax, and TOC. LECO's CS230-series carbon and sulfur determinator was used to measure the TOC content of well cuttings samples by following the formal standards process (CNS GB/T19145-2003).

Biomarkers analysis
The ten well-cuttings samples were also extracted for sedimentary OM characterization. Extracts were separated by liquid chromatography into aliphatic and aromatic fractions of hydrocarbons and further analyzed by gas chromatography-mass spectrometry (GC-MS). Both aliphatic and aromatic biomarkers were analyzed to evaluate the organic matter's source facies, thermal maturity, and deposition environment. The methods, techniques, and guidelines are similar to those reported in the literature [27][28][29] .

Maturity index (MI) calculation
Zhao et al. 30 developed a relationship between MI and Tmax to determine the kerogen type. They established two baselines, i.e., MI < 5 referring to the oil with some dissolved gases and MI > 7 indicating dry gases without any condensate. The values of MI between 5 and 6 indicate oil (or condensate) and wet gas.
Bulk density, neutron porosity, shear and compressional sonic, and deep resistivity logs, which are sensitive to maturity variations, were utilized in Bijnot-1 well to compute the MI in the Infracambrian rock by employing Equ. (1). In this equation, Zhao et al. 30 evaluated the total and effective porosities of the formation from the core samples and correlated it with porosity calculated from the density log in the same interval in order to determine the threshold values of MI.
In this study, data of samples with log density porosity greater than 9% with water saturation less than 75% were included in the calculations of MI by the following equation: where 'N' denotes the number of those samples with log density porosity greater than 9% and saturation of water less than 75%, φ n9i denotes the neutron porosity of those samples having log density porosity greater than 9%, and S w75i is the saturation of water of such samples.

Calibration of the seismic data
Both conventional logs (sonic and bulk density) and check shot data (time-depth relation) of Bijnot-1 well along with seismic line were loaded into the computer program to calculate synthetic seismogram. The acoustic impedance (AI) was obtained as the product of sonic velocity (V P ) and bulk density log, as given in Equ. (2). Using AI, the re ection coe cient (re ectivity) at each interface was then calculated by Equ.
(3). In the process, we extract the wavelet from the seismic dataset, speci cally seismic traces from seismic data. The wavelet extraction derived from digitized AI and bulk density logs was convolved with the series of re ectivity to compute a synthetic seismogram. The re ectivity was converted from the depth domain to the time domain using available check shot data from wells. Finally, the synthetic seismogram was then correlated with the seismic traces to build an impedance model. The AI curve closely matched with the lithofacies (e. g., high AI is related to dolomite and shale). Particularly for this dataset, we selected Ricker wavelet with a 2-millisecond sample rate, 25 Hz frequency, and 128-millisecond sample lengths to generate the synthetic trace. Substitution like sonic data calibration, datum, drift corrections, and seismic traces correlation was also performed during the construction of synthetic seismogram. Fig.  7 shows the synthetic seismogram, AI, re ectivity series, and the computed extracted wavelet for seismic line FABS-11.

Seismic data inversion
We applied a joint inversion strategy in this work that incorporates MLC and PSO algorithms to convert seismic re ection data to spatial variability in the AI and MI pro les in the Infracambrian SR formation over the studied region 31,32 . An open-source rock star seismic inversion software was used to predict the MI pro le. Seismic line and MI data (input) were entered into the inversion software, and the parameters according to the provisions of input data were adjusted to estimate MI (output).

Surface geochemical and geo-microbial survey
The abundance of hydrocarbon oxidizing bacteria ranges from 0 to 525×103 (cfu/g of soil sample) (Fig. 8a). The hydrocarbon gases measured by gas chromatography show concentrations of C2 through C4 ranging from 0 to 900 ppb (Fig. 8b).
Distribution maps of geo-microbial and soil gas concentrations show that the anomalies match and are well inversely correlated with each other (Fig. 8a, b). In fact, hydrocarbon micro-seeps provide suitable conditions for developing highly specialized microbial populations that consume light hydrocarbon gas as their only food source. As a result, these microbes multiply in anomalous concentrations in the nearsurface soils along with the geological structures containing hydrocarbons. Here we note that the obtained result represents the normal distribution between gas and microbial concentrations. This reverse correlation results from the consumption of the soil gas by the bacteria that oxidize hydrocarbons. The light hydrocarbons micro-seepage is consumed, particularly over the highest seepage area, allowing a high bacterial activity rate and depletion of hydrocarbon gas concentrations. The bacteria found just above the ow chimney would transform the petroleum gases, while the edges will show high gas concentration as non-utilized by bacteria where all microbial activity is prolonged.
The soil probe gas analyses show variable concentrations where C1 is absent, and the samples contain C2-C4 hydrocarbons. However, the amounts of migrated gases usually decrease in the order given as methane > ethane > propane > butane. We can attribute the de ciency of methane to its total consumption by the methanotrophic bacteria in the soil. Because it is the simplest hydrocarbon molecule, methane would be the rst and easiest component to decompose by bacteria causing its complete depletion in the soil gas samples in proportion to the other components.
The presence of heavier components of hydrocarbon gases, mainly propane and butane, favours a thermogenic origin of the detected gases. Nevertheless, another source of hydrocarbon gases could be detected on the surface, such as the biogenic origin.
To further investigate the origin of these surface gases, we included other studies such as geophysical characterization, including gravity modelling and seismic data interpretation, to explore if possible deep pathways exist, in addition to maturity estimation and biomarker analysis.
The soil probe gas concentrations are plotted and overlapped on the gravity anomaly map (also ties with seismic data interpretation and geological cross-section). The composite map of C2-C4 hydrocarbons and gravity anomaly shows a good correlation between high concentrations of C2-C4 hydrocarbons and the zone of gravity high, as well as major faults (Fig. 9). The C2-C4 hydrocarbons show high concentrations along the faults separating gravity highs and gravity lows.
The seismic lines (Fig. 6) along the surface geochemical study were also interpreted to show the structural features and possible migration pathways. The stratigraphic horizons were marked and tied using check shot velocities. The cross-sections along seismic line GG(n or th → south) and II from (east to west) were prepared and converted to depth. We note that salt-induced structural leads identi ed at the Infracambrian level show a continuous rise from the basement overlying stratigraphic horizons from west to east.
A geo-seismic depth section along II' lines reveals a basement high located at a distance of 40 km (Fig. 10a). The overlying Mesozoic-Tertiary sequence follows the rise in the basement. A four-way structural closure of 80 msec (~160 m) (vertical closure) is mapped in the eastern part of the study area.
A normal fault can also be observed in the seismic section. The Infracambrian top in this seismic section is identi ed at 1600 to 2200 m depth from east to the west.
A geo-seismic cross-section along with GG` is oriented in the NS direction (strike line). The cross-section shows that the basement starts rising from a distance of 15 km to the highest crustal part located at 25 km (Fig. 10b). The Infracambrian to Tertiary strata overlying the bulge in the basement shows prominent anticlinal four-way structural closure (as shown by the dip line from east to west). The thickness of the stratigraphic horizons remains constant. The Infracambrian top in this seismic section is identi ed at 1400 to 1600 m depth from north to the south.
We conclude that deep-seated faults overlying the basement high as well as near-surface diffusion provided pathways for C2-C4 hydrocarbon migration to the surface, supporting the deep thermogenic origin of these gases. Meanwhile, this indicates the presence of an active petroleum system beneath the earth's surface.

vitrinite re ectance
When modelling the burial history, we considered the regional erosional unconformity between the Cretaceous and Paleocene, the most critical erosional surfaces in the study area. In models development, we considered total erosions between 40 and 180 m. However, after conducting a comprehensive sensitivity analysis and trying out several scenarios, we found a negligible impact of erosions on the present-day maturity and temperature trends.
In Bijnot-1 well, a constant heat ow of approximately 72 mW/m 2 results in the best calibration between the measured and calculated vitrinite re ectance ( Fig. 11a and 12). Burial history and thermal maturity modelling of Infracambrian rock demonstrated that Sonia, Jodhpur, Bilara, and Hanseran reached the early stages of the oil generation window, with a maximum burial temperature of approximately 105 o C and calculated vitrinite re ectance of 0.60% R o (Fig. 11a).
For the maturity pro le, the Infracambrian rock (Sonia, Jodhpur, Bilara, and Hanseran Evaporite) fall within the early-stage of the oil generation window (Fig. 11b). The present-day calculated TTI is 35. Vitrinite re ectance (R o ) in the Infracambrian source rocks encountered in the studied wells of SE Pakistan suggests immature to the early stage of oil generation window (Fig. 12).

Rock-Eval pyrolysis
The results of Rock-Eval pyrolysis of Bijnot-1 well cuttings (  (Fig. 13). The HI versus Tmax plot reveals that more than half of the samples have reached an early mature stage of the oil window generation. A mixture of types II and III is present in SR Formation source rock with predominantly type III kerogen (Fig. 13a). This mixture with the predominance of type III explains the enrichment of the original OM in oxygenated components and, consequently, OI's high values. In the PI versus Tmax plot, the PI values fall between 0.08 to 0.46 and Tmax from 420 to 435°C, further con rming that Infracambrian SR Formation is thermally immature to the early-mature (stage of oil generation window) (Fig. 13b).  Fig. 14a. The n-alkanes distribution pro le shows the higher abundance of even-carbon number n-alkanes than odd n-alkanes. This suggest a reducing environment of organic matter deposition where even-carbon n-alkanes preserved preferred to odd n-alkanes. An exclusive abundance of lower carbon range alkanes ~C 16

Possible source facies
The distribution of tricyclic terpanes, steranes, and hopanes in cuttings extracts (1742 m) was applied to determine the variation of source facies in the SR Formation (Fig. 14b, c). Hopanes are present in higher relative abundance than tricyclic terpanes. The C 23 compound is the predominant peak among the tricyclic terpanes and is particularly abundant comparing to C 24 tetracyclic. The C 23 /C 24 is about 1.5.
This ratio re ects the importance of tricyclic terpanes concerning tetracyclic terpanes and is signi cant for the abundance of precursors of these compounds (tricyclohexaprenol) synthesized by bacteria as a stabilizer of microorganism membranes 36 . This indicates that the OM of SR Formation source rock could originate from a marine source facies enriched in microorganisms. Even so, the presence of C 24 in small but signi cant concentrations re ects a possible minor contribution of terrestrial OM sources. Peters and Moldowan 37 studied that the C 24 tetracyclic terpane is related to source rocks containing continental organic matter in respectable amounts. Also, C 20 and C 21 are relatively abundant, and C20 tricyclic terpane abundance is related to terrestrial plant sources . The abundance of C 25 and C 26 hopanes are almost similar, with a C 26 /C 25 tricyclic terpane ratios less than 0.9 36 , indicating a prevalent marine depositional environment. Similarly, the C 24 tetracyclic terpane is more abundant than the C 26 tricyclic terpanes at almost all depths, revealing the marine carbonates to marl facies.
The stable 17α(H), 18α(H), 21β(H)-28, 30-bisnorhopane were detected with low concentrations comparing to hopane ratios. This con rms that the depositional environments, in SE Pakistan, during the Infracambrian are generally suboxic. In contrast, the high C 28 bisnorhopane is related to reducing physicochemical conditions in organic matter deposition environments 38 . In the pentacyclic terpanes,

Thermal maturity
Aliphatic and aromatic hydrocarbon parameters were applied to evaluate the thermal maturity of SR Formation sediments. Isomerization at C-22 positions in C 31 to C 35 17 α(H)-hopanes 42 occurs prior to most biomarker reactions and is used for thermal maturity assessment as isomerization at C 20 position in regular steranes. In response to temperature augmentation, the biologically existing 22R con guration transforms to 22R and 22S mixture. As a result, the 22S/(22R + 22S) homohopane ratio could generally be used in determining maturity. However, homohopane isomerization attains equilibrium at the early maturity stage; therefore, 22S/(22R + 22S) homohopane ratio cannot be used to determine highermaturity stages 29  On the other hand, the Ts/(Ts+Tm) ratio shows values around 0.35, indicating that the OM has not reached the oil generation window yet, given that Ts is more stable to thermal maturity than Tm and the ratio Ts/(Ts + Tm) increases with increasing maturity. These immature Precambrian rocks are consistent with previous studies from the basin 35 . However, if we examine the rest of the thermal maturity biomarkers parameters of the samples from the SR Formation, especially the deeper ones, we see a gradual increase in thermal maturity. The pregnane (C 21 -sterane) and methyl-pregane (C 22 -sterane) are also used as thermal maturity biomarkers since they are characteristic of light oils 43,44 . In the laboratory, these steranes are produced by thermal cracking at 300°C for a long period 43 . In the OM of the SR Formation, C 21 and C 22 steranes were detected in high concentrations (Fig. 14c).  AI inversion transforms seismic data into pseudo acoustic impedance logs at every trace. Fig. 15a, b shows the interpreted seismic line of FABS-11 and inverted AI using the MLC and PSO inversion strategy.
It has been demonstrated that the AI in the zone of interest, i.e., Cambrian and Infracambrian sequence (900 to 1460 ms), varies from 9000 to 13500 (m/s) × (g/cc). These variations of AI are associated with signi cant lithological changes to sandstone, shale, limestone, dolomite, and sandstone with subordinate shale lithofacies. It is to be noted that the low AI (9000-11500 (m/s) × (g/cc)) between time interval 1250 to 1450 ms indicates perhaps a potential hydrocarbon-bearing zone at this particular interval (arrow).
However, the overlying red-yellow horizons just above the blue horizons (arrow) reveal higher values of AI (12500-13500 (m/s) × (g/cc)) and is, hence, most likely an impermeable seal rock for a potential reservoir. The harder rocks, i.e., compact limestone, have higher AI than sandstone and clay 31 . Note that irregular patterns of low and high AI re ecting interbedded limestone and dolomite layers are wellcaptured by MLC and PSO inversion strategy (indicate with arrow) (Fig. 15b). Since a nonlinear relationship exists between the interbedded thin layers and seismic waveforms, PSO could provide a nonlinear projection correlation to determine the interbedded layers from the seismic waveform 45 . Good reservoir quality areas and thin layer beds show up clearly in the Cambrian and Infracambrian interval (arrow).

MI prediction
MI is an essential parameter for determining the thermal maturity levels of hydrocarbon source rock and correlating it with thermal maturity evaluated from geochemical methods. Fig. 16 shows a graphical comparison between the calculated MI using Equ. (1) and measured Tmax in the Infracambrian SR Formation in Bijnot-1. It is shown that the measured Tmax (red dots) falls within the early-stage of the oil generation window in the range of 420-433°C with MI (maroon line) 2-3%.
In Fig. 17, the inverted AI surface is converted into an inversion solution for a MI map using the MLC and PSO inversion strategy 30,46 . It is shown that the deeper horizons reveal fair to good MI values, i.e., 3 to 5% with low AI, and correspond to good reservoir quality area (arrow). It should be noted that blue and sky blue colour layers in the middle and upper part of pro le with MI < 2 correspond to the immature oil generation stage, whereas yellow and green colour layers in the lower part of pro le with MI > 4 correspond to the early oil generation stage. The results of MI concur with the earlier research of Abdizadeh et al. 46 and correlate well with the thermal maturity measured from rock-eval pyrolysis and biomarker in the studied interval. Generally, thermal maturity and MI increase as a function of depth; however, structural features such as structural highs and lows, salt-induced anticlines, and paleo highs could produce local thermal anomalies.
The comparison of original logs with inverted AI and MI from MLC and PSO for FABS-11 seismic sections is shown in Fig. 18. The plot indicates good concordance between the inverted (red line) and the computed (blue line) AI and MI. The overall correlation coe cients for both inverted and computed values are equal to 0.96 and 0.89 for AI and MI, respectively. This high correlation coe cient validates the calibration between inverted AI and MI from seismic and well log data and con rms the robustness and accuracy of the MLC and PSO inversion strategy.

Research Synthesis
The Infracambrian play in SE Pakistan has been previously explored to study the distribution and modelling of the source and reservoir rocks. Many deep wells targeting to penetrate Infracambrian rock were drilled on paleo-highs but did not encounter Infracambrian strata 6,9 . All wells were drilled to the uppermost column of the oil window and allowed minor in-situ heavy oil shows within the source rock beds. The timing of hydrocarbon generation versus trap formation is another reason for the failure of deep wells. The potential charging episode and maximum burial depth are likely pre-Permian, whereas the trap formation in the drilled structures is post-Permian.
The literature study from various authors revealed that the heavy oil shows from the Bijnot-1, Fort Abbas- steranes, monoaromatic, and diasteranes steroids, negative stable carbon isotope ratio, low pristane/phytane ratio, low diasteranes, low API gravity, and high sulfur) of Karampur-1, Baghewala-1, and southern Oman suggest similar environmental conditions with the same age, as shown in Fig. 19. sometimes absent due to erosion during the Cambrian and Permian unconformity, which could explain why the upper part is more oxygenated. It is uplifted and, as a result, more exposed. This exposure makes the OM in the upper part subject to alteration and explains the immature biomarkers features we detected in some samples and the possible minor terrestrial provision recorded.

Conclusions
We document the evaluation of source parameters from Neoproterozoic rocks (Infracambrian rock sequence) in SE Pakistan and adjacent BK-N Basin India, using geophysical data, cores, rock-eval pyrolysis, biomarker, and surface geochemical samples. Based on available data and techniques, we can draw the following conclusions: 1. The seismic and well log data interpretation identi ed the Infracambrian rock sequence at 1600 to       Conceptual play type of Infracambrian rock sequence in BK-N Basin India 6,9.

Figure 6
Page 27/31 The location of sampling sites along with the seismic and gravity pro les (blue lines show seismic crosssection) over the study area. Synthetic seismogram (red traces) for Bijnot-1 well, showing from left to right: true vertical depth in meter, sonic velocity (µs/ft) and bulk density (g/cm3) logs, re ectivity (re ection coe cient), acoustic impedance (AI), lithology, synthetic seismogram, traces from a part of seismic line FABS-11, and two-way time in ms.   The overlapping of C2-C4 hydrocarbons anomalies on the interpreted gravity anomaly map. The location of gravity pro les is shown in gure 6.  Burial history and thermal modelling of Infracambrian rock, (a) vitrinite re ectance, (b) maturity pro le.

Figure 20
Subsurface correlation of Neoproterozoic sequence (Nagaur, Hanseran, Bilara, and Jodhpur Formations) in the SE Pakistan and BK-N Basin India.

Figure 21
Distribution of various source and reservoir rock properties along with certain biomarkers parameters in Infracambrian rock sequence from SE Pakistan, India, and southern Oman.