Overpressure mechanism prediction based on well log and mineralogy analysis from drill cuttings of well NSE-001 in the North Sumatra Basin area, Indonesia

The well NSE-001 is an exploration well located in the North Sumatra Basin, Indonesia, which penetrates several formations with overpressure zones that have the potential as a drilling hazard in this area. The results of the pore pressure fracture gradient model show that the response of the wireline log data of well NSE-001 based on the normal compaction trend has two types of overpressure mechanisms, namely loading and unloading. Rapid sedimentation in the development of the Lower Baong Formation and Upper Baong Formation as thick and massive shale sequences in the North Sumatra Basin is the main factor in the generation of an overpressure zone with a loading mechanism. This correlates with the travel time in the sonic log which is relatively constant with response to the normal compaction trend at depth intervals of 620–1500 mMD (from the Middle Baong Sand Formation to the Top Gebang Sand Formation). The presence of sillimanite and kaliophilite minerals in bulk analysis from drill cuttings at depth intervals of 2200–3820 mMD contributed to the generation of hydrocarbon, while clay analysis showed that the mineral content of kaolinite + illite as alteration product (clay diagenesis) was higher than the mineral content of smectite. The results of the mineralogy analysis (bulk and clay) correlated with the response of the wireline log data of well NSE-001 to the normal compaction trend, which showed that at depth interval of 1500 m to 3200 mMD (from the Top Gebang Sand Formation to the Top Upper Belumai Formation) is an overpressure zone with an unloading mechanism, which was characterized by fluid expansion caused by clay diagenesis and possibly hydrocarbon generation.


Introduction
The North Sumatra Basin is a unique basin, which has high oil and gas potential, but is also an overpressure basin that affects drilling operations. According to Aziz and Bolt (1984) that the top of overpressure is mostly located in the massive shale sequence of the Post-Rift of the Lower Keutapang Formation. The diapir shale is evidence of the presence of overpressure, and indicates a stress regression in the Arun Limestone which lies beneath the Baong Formation under very high pressure. The characteristics of overpressure in the North Sumatra Basin are divided into two overpressure conditions in the Baong Formation, namely the Upper Baong Formation as a mild overpressure zone and the Lower Baong Formation as a high overpressure zone.
From the results of previous studies, an interesting phenomenon was found, namely in several wells in the North Sumatra Basin there were sonic log and resistivity logs reversals, each of which became high sonic and low resistivity (Hutasoit et al. 2013). It is proven that the well drilling in the North Sumatra Basin is in the overpressure zone, although the use of low mud density does not cause significant problems which indicates a mild overpressure zone. This phenomenon can be explained that the shale in the shallow section of hydrostatic pressure has experienced shallow compaction due to under-compacted carbonate cementation in the past, so that the overpressure in the zone is not too high, and is referred to as a mild overpressure zone (Syaiful et al 2020).
In identifying the overpressure mechanism in addition to using well log data analysis, it can also be done using mineralogy analysis, namely the hypothesis about the effect of temperature on clay diagenesis that causes the smectiteillite/kaolinite transformation (Storvoll and Brevik 2008;O'Conner et al. 2011).
The drilling target of well NSE-001 is to produce shale gas from the Lower Baong Formation and the Bampo Formation as source rock (Barber et al. 2005;Subroto et al. 2013;Yosandian et al. 2014). The presence of sillimanite mineral, which is brittle mineral due to alteration of clay at high temperatures (600-900 °C), which supports the process of hydrocarbon generation in shale (Buntoro et al. 2020). From the results of XRD analysis (bulk) shows the presence of sillimanite mineral from a depth of 2200 mMD in the Lower Baong Formation as source rock, which is also an overpressure zone.
The discussion in this paper begins by examining the research problem and the hypothesis to be proven. The methodology used in this study is then described, followed by wireline log response analysis, mineralogy analysis, and analysis of the pore pressure fracture gradient (PPFG) model.

Geology
The North Sumatra Basin contains several oil and gas fields that have drilling problems related to the overpressure zone (Aziz and Bolt 1984). Studies on regional overpressure are still very limited due to the limited data and the variation in the geological conditions of each area (sub-basin). The Baong Formation is a regional reservoir that has many overpressure problems. The sagging tectonic phase in the early Miocene made the sub-basin inundated and became a deep-sea environment. Deep-sea deposits such as shale and claystone are very well developed to form Lower Baong. In the Middle Miocene there was a regional uplift where sandstone deposits were formed which then rapidly subsided to produce fine-grained deposits such as shale and claystone (Upper Baong). The potential overpressure zone is in the Upper Baong Formation and the Lower Baong Formation because it has thick and massive deep-sea shale lithology.
The location of well NSE-001 is in the Tamiang Sub-Basin, North Sumatra Basin, Indonesia (Asset-1 PT. Pertamina EP Area). In the Tamiang Sub-basin the Baong Formation is very well formed. The development of the Lower Baong Formation and the Upper Baong Formation as thick and massive shale sequences is a major factor in the presence of an overpressure zone (Gurning and Ramdhan 2017).
There are two geological factors that may trigger the overpressure zone, namely the loading mechanism caused by rapid sedimentation and the unloading mechanism caused by changes in clay minerals and hydrocarbon generation. However, for hydrocarbon generation, geochemical analysis data of source rock is needed so that it can be more valid.
The North Sumatra Basin is one of the three back arc basins in Sumatra. Regionally, the three basins were formed at the beginning of the Tertiary period with Pre-Tertiary bedrock. The geometry of the North Sumatra Basin tends to be oval with a general direction of Southeast-Northwest and North-South extending eastward to off the coast of Sumatra. This basin is bounded by the Bukit Barisan to the west, the Asahan High to the south-southeast, the Mergui Arc to the north and the Malacca Shelf to the east.
The development and tectonic history of the North Sumatra Basin began in the late Paleozoic-early Mesozoic where the Indo-Australian Plate and the Eurasian Plate were already actively converging. The activity of magmatism and volcanism resulted in the intrusion of granite and granodiorite in Paleozoic and Mesozoic sedimentary rocks to form metamorphic rocks which will act as bedrock. The Indian and Eurasian collisions caused a change in the direction of subduction from the dominant direction north to northeast so that it was oblique to Sumatra.
According to Pulunggono et al. (1992), there are three tectonic regimes that involve the development of the basin. First, the extensional regime that dominates the basin configuration and early filling from the Late Cretaceous to the early Tertiary. Activities in this regime form a system of graben and horst, which traps and grows from reefs which will later become a hydrocarbon kitchen (kitchen area). Second, the intra Miocene regime that occurred during the Miocene. In this regime, uplift at the edges of the basin is followed by the deposition of clastic sediments. In the latter regime, namely the compressional regime that occurred during the Plio-Pleistocene, it caused an anticlinorium en-echelon pattern, northwest-northwest-east-southeast, and traps were formed. This ended Tertiary deposition in the North Sumatra Basin (Satyana 2009).
The regional stratigraphy of the North Sumatra Basin ( Barliana et al. 1999) from oldest to youngest can be seen in Fig. 1. The lithostratigraphic sequence of the North Sumatra Basin from oldest to youngest is as follows: The Tampur Formation consisting of limestone and massive dolomite in the Late Eocene to Early Oligocene. The Bruksah Formation is interspersed with coarse-grained carbonate sandstones and thin intercalations of coal. The Bampo Formation is dominated by black shale, generally associated with pyrite and calcareous nodules, which were deposited during the Early Miocene. The Peutu Formation consists of claystone and calcareous siltstone which contains many fossils and gluconite of Early Miocene-Middle Miocene age. The Baong Formation consists of shale with sandstone inserts of Middle Miocene age. The Keutapang Formation consists of finegrained sandstone and shale of Middle Miocene-Upper Miocene age. The Seurula Formation consists of alternating thin layers of sandstone and shale of Early Pliocene age. The Julu Rayeu Formation is an alternation of sandstone, shale, gray clay, tuff, plant remains, and lignite, as well as a small amount of conglomerate of Late Pliocene age.
Based on the explanation of the Regional Tectonic Order and Regional Stratigraphy of the North Sumatra Basin above, it is believed that well NSE-001 is in the Post-Rift to Syn-Rift phase, where in general the folding pattern is in the WNW-ESE direction, forming an en-echelon anticlinorium fold pattern and a horizontal fault (wrench fault). The direction of the horizontal fault movement (wrench fault) is NW-SE according to the basement fault pattern, this is related to the Paleogene graben system (Harding 1974). This condition will cover the blocks that are lifted due to the compression regime of the subduction system, so that the stress that works up to now is Pliocene-Pleistocene stress (Moulds 1989). Thus, it can be concluded that the North Sumatra Basin tectonic is a type of Simple Shear stress (Harding 1974), so that the position of well NSE-001 is believed to have the potential for overpressure, because it is in a simple shear stress system which is influenced by the NS structural pattern with the main stress direction NE-SW (Fig. 2).

Study area
The well NSE-001 is located in the North Sumatra Basin, Indonesia, where the well NSE-001 is an exploration well with a target shale gas reservoir in the Lower Baong Formation and Bampo Formation which are the main source rock in the North Sumatra Basin (Barber et al. 2005;Subroto et al. 2013).
This well was drilled on October 4, 2016, reaching a final depth of 3796.1 mTVD on February 22, 2017. The location of well NSE-001 can be seen in Fig. 3. The following is the formation penetrated in the drilling of well NSE-001 according to the stratigraphy of the North Sumatra Basin from the youngest to the oldest (Fig. 1) as follows: • Keutapang Formation (0-404 mTVD). The Keutapang Formation consists of fine-grained sandstone and shale.  Basin is a type of simple shear stress (Harding, 1974 The lithology penetrated in drilling of well NSE-001 is dominated by shale rock types, with some formations found in sandstone and carbonate rock lithology.

Literature review
Pore pressure fracture gradient (PPFG) model prediction method Many correlations and models have been developed for the prediction of abnormal pressure of various parameters (Hottman and Johnson 1965;Matthews, W.R. and Kelly 1967;Pennebaker 1968;McClendon 1971;Zamora 1974;Eaton 1975). The method has limitations; for example, some models can only be used in clean shales (Table 1).
Some pore pressure prediction methods only can be applied to the pressure generated by the under-compaction mechanism or some of them do not apply in unloading formations. Therefore, the researchers tried to use AI (Acoustic Impedance) from seismic data to predict pore pressure with good precision. To predict the pore pressure, four input parameters are used, namely: depth, Gamma Ray log, density log, and sonic log.
From some literature shows that the parameters used to predict pore pressure from well log data. Most methods require both normal and abnormal trends, to estimate pore pressures.

Overburden gradient prediction calculation
The formation overburden stress gradient (OBG) in onshore drilling can be estimated to vary from 1.0 to 1.1 psi/ft (0.0227-0.025 MPa/m). But in this study, OBG was calculated from the density log (RHOB) data. Because the density log data do not start from the surface/seabed, so in some equations to calculate the pseudo-RHOB shallow section by integrating the LWD/wireline RHOB log data.

Pore pressure (PP) prediction calculation
There are several methods for predicting pore pressure calculations, namely: Equivalent Depth, Eaton's, and Bower's methods.

Bower's method
The Bower's method uses sonic velocity empirically to calculate the value of the vertical effective stress, which is then subtracted from the overburden (total vertical stress) to determine the pore pressure.
This method can be applied to predict the pore pressure caused by compaction disequilibrium or due to some source mechanism.
"Loading" method In the "loading" method only two empirical parameters are required, if the overpressure is caused by compaction disequilibrium "loading".
The values of the two parameters can be empirically determined in the compaction trends analysis or selected based on experience in offset wells.
"Unloading" method Further information and additional empirical parameters are required when stress is caused by the source mechanism. (1)

3
The "Maximum velocity depth" parameter, d maxv controls whether "unloading" has occurred or not. If d maxv > depth, "unloading" has not occurred. If d maxv < depth then the assumed "unloading" behavior of pore pressure is calculated as follows: where PP = Pore pressure gradient (psi/ft, lb/gal), (kPa/m, g/cc). OBG = Overburden gradient (psi/ft, lb/gal), (kPa/m, g/cc) DT = Sonic travel time (microsec/ft, microsec/m) DT ml = Sonic travel time on V max . A, B, U = Empirical values. V max = The velocity at which unloading occurred for sediments buried at depths. d maxv = Depth at which unloading has occurred. Depth = TVD in appropriate units.
In this case, the vertical effective stress of the sediment will be less than in the past and is called the "unloaded state".
To perform this analysis, it is necessary to know the sediment's previous maximum effective stress σ max , and it is necessary to determine the effective stress behavior of the sediment "unloading" velocity, which is determined by the unloading parameter, U. The value of σ max is calculated from the normal compaction response and the user defined value of; and the value of U is determined empirically.

Fracture gradient (FG) prediction calculation
In predicting fracture gradient (FG) there are several methods that can be used, including: Eaton's Method, Daine's Method, and Matthew and Kelly Method.

Eaton's method
Fracture pressure is the pressure that forms the fracture, when the minimum compressive stress and tensile strength are exceeded by the pore fluid pressure. Practically speaking, the fracture pressure is the limit of pressure that a formation can withstand from a given column of mud.
Prediction of fracture pressure, LOT trend can be valid and representative using the Eaton calculation method.
Daine's method (1982) The Daine's method in predicting the fracture gradient by superposing the horizontal tectonic stress t to the Eaton equation. Expressing in terms of stress, as "the minimum stress in the borehole that must be resisted opening and extending the existing fault", which can be written in the following equation: where σ = Vertical effective stress. σ 1 = Superimposed tectonic stress. ν = Eaton's Poisson's ratio. P = Pore pressure.
In determining the value of the effective pressure ratio (K), Daine uses an additional parameter to describe the lithological condition due to the influence of tectonic activity. For example, there are folds and faults. Therefore, the fracture pressure value obtained from the Daine's method will depend on the lithology, depth, and tectonic activity.

Matthew and Kelly method
The Matthew and Kelly method is a fracture gradient predication method using a variable of "matrix stress coefficient (K i )", equivalent to the effective stress coefficient to calculate the fracture gradient of sedimentary formations: where FG = Fracture gradient (psi/ft or lb/gal), (kPa/m or g/cc). PP = Pore pressure (psi/ft or lb/gal), kPa/m or g/cc). OBG = Overburden gradient (psi/ft or lb/gal), kPa/m or g/ cc). K i = Matrix stress coefficient (dimensionless).

Determination of overpressure generating mechanism
The wireline log data response will show two basic criteria for the overpressure formation mechanism, namely "loading" and "unloading".

Characteristics of loading mechanism
The wireline log shows that the porosity data is constant and the pore pressure is parallel to the lithostatic stress under top overpressure, which is based on simple assumptions, including uniform lithology, no fluid outflow, no temperature change, and no diagenesis. A pore pressure profile parallel to the lithostatic pressure associated with constant porosity is usually estimated as the young mud-rich sequences undergoing rapid burial in sedimentary basins and is generally considered as disequilibrium compaction. The loading mechanism will show a constant trend of sonic log and density log, which is caused by the retention rate of the pore fluid compaction rate of the rock.
Accumulation of sediment causes an increase in vertical stress, or gravitational loading, which can potentially lead to the presence of overpressure due to unbalanced compaction (disequilibrium compaction). In the unbalanced compaction process, the pressure increase is the result of the competition between the fluid liberation rate and the rate of vertical compaction that occurs due to gravitational loading due to the ongoing sedimentation process.
The increase in lateral stress as a result of tectonic compression can also cause overpressure, as occurs in the disequilibrium compaction mechanism. This is also caused by gravitational loading. Sediment can be compacted horizontally in the same way as vertical compaction, but if the fluid that fills the rock pores cannot escape (migrate due to rock compaction), the pore pressure will increase and the sediment compaction rate will decrease. Figure 4 shows the pressure depth profile and wireline log response, where overpressure is caused by loading mechanism/disequilibrium compaction (Osborne and Swarbrick 1997).

Characteristics of unloading mechanism
The unloading mechanism can cause high levels of overpressure through processes that cause load transfer from the contact between rock grains to the fluid that fills the rock pores (Swarbrick et al. 2002). The influx of pore fluids or conversion of a solid matrix to a fluid will spontaneously cause an increase in pore pressure, if the fluid that fills the rock pores cannot escape/isolate. If some of the load on the rock grains turns into fluid, then some of the stress that was previously accommodated by the contact between the rock grains will also move to the pore fluid. In other words, the pore fluid volume increases due to load transfer. This causes an increase in pore pressure. The unloading mechanism causes a decrease in the value of the effective stress. In contrast to the loading mechanism which does not cause a decrease in the value of the effective stress, but only holds the rock to compact normally. The unloading mechanism will show a sonic log where the direction of the log response is opposite from the normal trend because the pore fluid is expanding faster than the compaction rate.
Geological processes that cause pore fluids to increase by load transfer and lead to unloading mechanism include clay diagenesis processes, such as the process of transforming smectite into illite (Lahann and Swarbrick 2011), hydrocarbon generation (Swarbrick et al 2002), and lateral transfer or vertical transfer (Bowers 1995). Figure 5 shows the pressure depth profile and wireline log response, where overpressure is caused by unloading mechanism/fluid expansion (Osborne and Swarbrick 1997).

Mineralogy analysis method
The results obtained from measurements using the X-ray diffraction (XRD) instrument are diffractogram graphs. The diffractogram is the output which is a graph between 2θ (diffraction angle) on the X-axis versus the intensity on the Y-axis. The intensity of the continuously diffracted X-rays is recorded as an example and the detector rotates through their respective angles (Fig. 6).
A peak in intensity occurs when the mineral contains a lattice with d-spacings corresponding to X-ray diffraction Fig. 4 The cartoon shows a depth profile and wireline log response, where overpressure is caused by loading mechanism/disequilibrium compaction: constant porosity values in the neutron porosity log, con-stant travel time values in the sonic log, and constant density values in the density log (Osborne and Swarbrick 1997) at a value of although each peak consists of two separate reflections (K α1 and K α2 ), at values smaller than 2θ locations to the location of the overlapping peak with K α2 appears as a mound on the K α1 side. Greater separation occurs at higher values of θ (theta angle).
The reading of mineral peak XRD (Fig. 7) is by reading the theta angle (deg) on the X-axis and intensity (cps) on the Y-axis at the highest peak/peak for which you want to know the name of the mineral. In determining the name of a mineral based on the measured peak (theta, deg), which is calculated by using the Bragg equation (1913), so that the d-spacing is obtained in Angstroms (Å). Furthermore, the d-spacing value can be used to determine the type of mineral from the Handbook of Mineralogy (Bladh et al. 2001).
After identifying the type of mineral, for example: Quartz, a semiquantitative calculation of the percent of the Quartz mineral can be carried out using a percentage calculation based on the count per second (cps) of the total identified minerals (Fig. 8).

Research method
The research method workflow is schematically shown in Fig. 9. The research workflow describes the research stages starting from data collection, literature review, determination of overpressure mechanism, prediction method selection, pore pressure validation, fracture pressure estimation, fracture pressure validation, and the final result is a pore pressure fracture gradient (PPFG) model.  (Nurcholis 2017) From the drill cuttings data were performed XRD analysis (bulk and parallel clay-oriented method) to determine bulk minerals and clay minerals and used to identify the phenomenon of clay diagenesis phenomenon (smectite to kaolinite/illite transformation) and hydrocarbon generation which has implications for the overpressure mechanism.

Wireline log response analysis
The response of wireline log data of well NSE-001 based on the normal compaction trend (NCT) shows that there are two types of overpressure mechanisms, namely loading and unloading (Fig. 10). The loading and unloading intervals are detailed in Table 2. The loading mechanism is shown by a constant trend of sonic log and density log, which is caused by rapid sedimentation (disequilibrium compaction).
Meanwhile, the unloading mechanism is shown by trend of sonic log and density log in the opposite direction from the normal trend response, which is thought to be due to the influence of fluid expansion, caused by the smectite to kaolinite/illite transformation (clay diagenenis) and hydrocarbon generation.
From the measurement data of well NSE-001 with Mean Sea Level 0 m, Drill Floor 34.35 m above Datum, the NSE-001 Max Recorded Bottom Hole Temperature (BHT) well is 126.67 °C at a depth of 2565.35 mMD, so it can be estimated that the temperature of 100 °C is at a depth of 2035.71 mMD (Fig. 10).
In the North Sumatra Basin, there is an overpressure phenomenon in the Middle Miocene Lower Baong and Keutapang formations at depth intervals of 4000-8000 ft, which is caused by rapid sedimentation/burial of sediments (Aziz and Bolt 1984). This can be seen very clearly in the wireline log data response of well NSE-001 based on the normal compaction trend (NCT), where at a depth interval  Fig. 9 The research workflow explains starting from data collection, determining the overpressure mechanism and the transformation of clay minerals which have implications for the overpressure mechanism of 620-1500 mTVD (from the Middle Baong Formation to the Top Gebang Sand Formation) indicating that the travel time in the sonic log is relatively constant, and referred to as the loading mechanism (Fig. 10).

Mineralogy analysis
The very thick rapid sedimentation of the Baong Formation is coupled with a high thermal gradient, so it is believed to have implications for the presence of overpressure in the North Sumatra Basin (Aziz and Bolt 1984). Based on the crossplot density (gr/cc) versus transit time (msec/m) it shows that the Upper Baong Formation is located in the smectite rich mud rock as a low overpressure zone, while the overpressure is present in the Lower Baong Formation which is located in the kaolinite + illiterich. Thus, it can be concluded that the overpressure in the Lower Baong Formation is caused by the unloading mechanism, due to the smectite to illite/kaolinite transformation (Hutasoit et al. 2013). Overpressure characteristics in the Tamiang Sub-basin, North Sumatra Basin indicates that there are two overpressure mechanisms, namely the loading mechanism caused by the rapid sedimentation rate in the basin and the unloading mechanism caused by clay diagenesis and hydrocarbon generation (Yosandian et al. 2014).
The results of XRD analysis of drill cuttings of well NSE-001 for bulk analysis (Table 3) show the significant presence of sillimanite and kaliophilite minerals at depth intervals of 2200-3820 mMD (Gebang Sand Formation, Lower Baong Formation, Belumai Formation, and Bampo Formation) which has implications for clay diagenesis (smectite to illite/ kaolinite transformation) and hydrocarbon generation (Buntoro et al 2020). From the results of clay analysis of drill cuttings at depths ranging from 2200 mMD downwards with temperature conditions reaching 100 °C (Table 4), it shows that the percentage of smectite mineral content is low (low smectite content), while the total content of kaolinite + illite is very significant (high kaolinite + illite content). The presence of significant kaolinite + illite minerals is an alteration product (Fig. 11) which has implications for the overpressure mechanism in the well NSE-001.
Based on mineralogy analysis (XRD) at depth intervals of 2200-3820 mMD (from the Gebang Sand Formation to the Bampo Formation), both bulk analysis and clay analysis showed a correlation with the response of the wireline log data of well NSE-001 to the normal compaction trend (NCT). The presence of sillimanite and kaliophilite minerals  Unloading in bulk analysis contributed to the generation of hydrocarbons, and clay analysis showed that the mineral content of kaolinite + illite as alteration products (clay diagenesis) was higher than the mineral content of smectite. It can be concluded that the depth interval of 1500 m to 3200 mMD is an overpressure zone with an unloading mechanism, which is characterized by fluid expansion caused by clay diagenesis and the possibility of hydrocarbon generation.

Summary and conclusions
• Based on the response of the wireline log data of well NSE-001 to the normal compaction trend (NCT), it shows that there are two types of overpressure mechanisms, namely loading and unloading. At a depth interval of 620-1500 mMD (from the Middle Baong Sand Formation to the Upper Gebang Sand Formation), where the travel time of the sonic log is relatively constant with response to the normal compaction trend (NCT), and is referred to as a loading mechanism caused by rapid sedimentation. Furthermore, at a depth interval of 1500-3200 mMD (from the Upper Gebang Sand Formation to the Upper Top Bawah Formation) the response of the wireline log data of well NSE-001 to normal compaction (NCT), where the travel time value is opposite to the sonic log, and the density value increases and decreases slightly at the bottom of the density log. Based on the response from the wireline log to the NCT, this overpressure depth interval is referred to as the unloading mechanism. Bampo Fm.
• The results of mineralogy analysis of drill cuttings at depth intervals of 2200-3820 mMD showed the presence of sillimanite and kaliophilite minerals in bulk analysis that contributed to the generation of hydrocarbon, while clay analysis showed that the mineral content of kaolinite + illite as alteration product (clay diagenesis) was higher than the mineral content of smectite. It can be concluded that the depth interval of 1500 m to 3200 mMD is an overpressure zone with an unloading mechanism, which is characterized by fluid expansion caused by clay diagenesis and the possibility of hydrocarbon generation.  • Rapid sedimentation in the development of the Lower Baong Formation and Upper Baong Formation as thick and massive shale sequences in the North Sumatra Basin is the main factor for the overpressure zone. Based on the Bottom Hole Temperature (BHT) measurement in the well NSE-001, it is estimated that a temperature of 100 °C at a depth of 2035.71 mMD will have implications for clay diagenesis and hydrocarbon generation which contributes to overpressure.