Hydrocarbon prospectivity of the miocene-pliocene clastic reservoirs, Northern Taranaki basin, New Zealand: integration of petrographic and geophysical studies

In this study, it is aimed to characterize the Early pliocene sandstone (EP-SD) and the Late Miocene-Early Pliocene Mangaa sandstone reservoirs and the efficiency of their sealing cap rocks using the petrographical and petrophysical data of these sandstone zones in northern Taranaki basin, New Zealand. The prospective potential reservoirs were studied using impregnated thin sections, XRD data analysis, and well log data (self-potential, gamma-ray, sonic, density, neutron, shallow\deep resistivity and PEF) to characterize the reservoir zones, in addition to Mercury intrusion capillary pressure data (MICP) to check the efficiency of some potential seals. The EP-SD and the Mangaa sandstone units are typically poorly consolidated very fine sandstone to siltstone, with porosities averaging 25%. The sands are composed of quartz (38.3–57.4%), with common feldspars (9.9–15.2% plagioclase, and 2.7–6.3% K-feldspars) and up to 31.8% mica. In Albacore-1 well to the north of the Taranaki Basin, the Mangaa formation includes three separate for each of the EP-SD zones (EP-SD1, EP-SD2, and EP-SD3), and the Mangaa sequence (Mangaa-0, Mangaa-1, and Mangaa-2). The thin section studies indicate that, the studied samples are grouped into greywackes, arenites and siltstone microfacies with much lithic fragments and feldspars, sometimes with glauconite pellets. From the XRD data, it is achieved that the mineral composition is dominated by quartz, mica/illite, feldspars, and chlorite. The petrophysical investigation revealed absence of pay zones in the EP-SD zones, and presence of thin pay zone with net thickness 5.79 m and hydrocarbon saturation of about 25.6%. The effective porosities vary between 23.6 and 27.7%, while the shale volume lies between 12.3 and 16.9%. Although the shale content is relatively low, the relatively high API (50–112 API of average 75 API) is contributed by the relatively high K-feldspar content and intercalations with thin siltstone and muddy siltstone beds. Sealing units include the intra-formational seals within the Mangaa sequence, mudstones and fine grained units overlying the Mangaa and further intra-formational mudstones, within the shallower EP-SD units. The efficiency of these seals indicates the capability to trap 16.4–40.6 m gas or 17.4–43.0 m oil which is relatively low in correlation with their efficiency in the central parts of the Taranaki Basin Overlying the primary seals, mudstones of the Giant Foresets Formation provide additional regional seal.


Introduction
The Great Taranaki Basin represents one of the largest sedimentary basins along the western coast of New Zealand. The basin is considered the lone hydrocarbon-producing basin in New Zealand; it has an areal extent of about 330,000 km 2 . It is bounded to the south by the West Coast basins and to the north by Reinga Basin (King and Thrasher, 1996;Funnell et al., 2001;Uruski et al., 2003;Radwan et al., 2021aRadwan et al., , 2022b (Fig. 1). Up to date, more than 400 offshore and onshore exploration and production wells were drilled within the Taranaki Basin. However, the Taranaki Basin is still potential for more exploration and discoveries.
The seismic survey indicated the presence of thick sedimentary sequences filling a depocentre for at least 150 km along the axis of the basin (Sutherland et al, 2001;Uruski et al., 2003). Among these, the Upper Miocene-Lower Pliocene reservoirs in the Taranaki Basin were rarely studied. Thereby, this study concerns with the potentiality and the reservoir quality of the Miocene-Pliocene Mangaa sandstone reservoir which is considered a main gas-bearing target in the basin. The Upper Miocene-Lower Pliocene deep-water Mangaa sandstone is located in the North of Taranaki Basin and attains an average thickness of about 340 m; it is built of very fine sandstones to siltstone of seafloor turbidity fan (King and Thrasher, 1996;Radwan et al., 2022b).
The Mangaa sandstones sequence of the offshore Albacore-1 well in North Taranaki is chosen as a case study for the reservoir characterization of the Mangaa Formation in this basin. These sandstones are fine to medium-grained, poorly consolidated, and are highly porous of porosity values varying from place to place of average 27.5% (Radwan et al., 2022b).
The studied Albacore-1 well has been drilled through a Miocene-Pliocene structural trend associated with some volcanics and volcaniclastics. The Mangaa sandstones in this well can be separated into three separate seismic units; these are named Mangaa-0, Mangaa-1 and Mangaa-2.

Geologic setting
The Taranaki basin is one of the most important hydrocarbon-prospective sedimentary basins in the offshore of New Zealand. It is named after the Pleistocene Mount Taranaki volcano that forms the main parts of the Taranaki Peninsula which is surrounded by most of the oil and gas fields in New Zealand (King and Thrasher, 1996;Uruski and Wood, 1991;Uruski, 2000). Up to 90% of its extension is offshore, whereas only the Taranaki Peninsula and the Northwest Nelson area are located onshore (Isaac et al., 1994;Uruski et al., 2003). The onshore Taranaki Basin is primarily overlain by the Pleistocene Mount Taranaki andesite, the Kaitake, and the Pouakai ranges. It is bounded to the east by the Taranaki Fault and the Wanganui Basin (Radwan et al., 2022b) (Fig. 1). The main Taranaki Fault is considered the largest structural compressional feature within this area. It has shown remarkable activities during various episodes. Its most recent reactivation has been caused by the convergent movements that have been driven by the modern plate boundary development. The Neogene deposits extend to the east crossing the main Taranaki fault and merging with the Wanganui and King basins. Besides, the western boundary of the Taranaki Basin is generally considered as a shelf break merging into the New Caledonia Basin.
The Taranaki Basin is divided into two parts; an Eastern Mobile Belt and a Western Stable Platform. The Eastern Mobile Belt is a part of the Australian plate that has been deformed by the colliding between the Pale and the Pacific plate. Its complex morphology is due to its modification by many tectonic episodes including the Neogene marginrelated convergent events and the early rift-drift phase (King and Thrasher, 1996;Radwan et al., 2022b).
The Late Pliocene event reflects the present-day crustal geometry, where the Late Miocene phase indicates the start-up and geometry the present day tectonic setting of the Taranaki Basin (Knox, 1982).
Throughout the depositional history of Taranaki Basin, its predominant paleo-environments in the eastern and the southern portions were affected a complex eustatic and tectonic interplay, while the western and the northern areas retained connected to the sea (Palmer, 1985;King, 2000;Palmer and Geoff, 1991;Radwan et al., 2022b).
The Taranaki basin has been filled by a Cretaceous to Neogene sedimentary cover. In its onshore part, these sediments are overlain by Quaternary alluvial, aeolian, marine deposits, and volcaniclastics (Pilaar and Wakefield, 1978;Knox, 1982;Armstrong et al., 1996). The paleoenvironmental data indicates that these sands were deposited in bathyal environments as turbidites near the base of the Pliocene Foresets slope sequence (Radwan et al., 2022b).
The Miocene sediments in Taranaki Basin represent a set of regressive sequences accumulated in episodes of increasing rates of sediment influx that have been triggered by the tectonic events, where clastics reworked from the uplifted areas in the vicinity with large volumes and accumulated in the centre of the basin depocentre. Upwards, the Pliocene-Pleistocene sedimentary deposits in the Taranaki Basin comprising the Giant Foresets, Matemateaonga, Tangahoe, and Mangaa formations (King, 1988). In general, the sedimentary sequence of the Taranaki Basin is primarily regressive and controlled by the detrital supply volumes and the interplay between the tectonic subsidence and the sediment influx. The Mangaa sandstone and the Giant Foresets sequences are composed of deep water fan sandstone alternated with a slope channel sandstone turbidites and basin floor shale the basin floor (Radwan et al., 2022b).

Data and methodology
The present study aims at characterizing the Mangaa sandstone in Taranaki Basin by utilizing different data sets including petrographic, and petrophysical data. Thereby, the Albacore-1 well which has complete data sets was chosen as a case study for characterizing the Mangaa sandstone.

Fig. 2
A generalized stratigraphic column of the Neogene sediments in the Taranaki Basin (King and Thrasher, 1996). The Mangaa sandstone, the core of the present study, is outlined by a red box

Petrographic studies
For the petrographical studies, ten samples were selected, and slabbed, representatively for the Mangaa sandstone in the Albacore-1 well. To describe the different pore types and to assess the visual porosity rough estimation, the prepared slabs were impregnated by a blue-dyed epoxy. These slabs were then prepared as thin sections and then the visual porosity and the main mineral components were estimated using a polarized microscope supplied with a PC monitor enabling a detailed estimation for the different mineral constituents.
Also, the diagenetic features and the grain-to-grain relationships were defined. The longest dimension of 100 grains/ thin section was measured using a grain size comparator and then the grain size and sorting degree were accurately estimated.
Besides, a semi-quantitative estimate of the overall mineral composition of the studied samples (grinded into very fine size) was determined precisely for the whole rock and clay fraction (< 2 μm) using the X-ray diffractometer (XRD) with a scanning speed of 0.01•/sec in the range of 2θ = 5-85° at 40 kV.
The mineral composition of eight samples of the probable reservoir zones (EP-SD-1, EP-SD2, EP-SD3, Mangaa-0, Mangaa-1, and Mangaa-2) were examined using the X-ray diffraction (XRD) analysis. In combination, these two analyses give a quantitative estimate of the overall mineralogy of the samples. The mineralogy, estimated from the XRD data, is for each sample as a whole, so the identified mineral composition is for the main rock sample constituents and also for the multi-mineral lithic fragments, not necessarily for the rock sample itself.

Well log analysis
To examine the Early Pliocene and the Mangaa sandstone reservoir properties and to define the prospective flow units, net-pay thickness, lithology, porosity, shale volume, and the overall fluid types and saturation, a log-based petrophysical analysis has been carried out on the available well logging data of the Abacore-1 well.
For the present study, a complete set of conventional well log data was available including Gamma-ray (GR), Selfpotential (SP), neutron (NPHI), density (RHOB), photoelectric factor (PEF), sonic (DT), and resistivity logs (shallow and deep resistivity, LLS and LLD, respectively). To carry out the quantitative and qualitative reservoir evaluation, the las files of these logs were imported into the Interactive Petrophysics software (IP 4.5) of Senergy Co. Then, the sand zones were defined and separated from the shaly zones. The petrophysical analysis was performed using the relevant equations and standard procedures that have been published by many authors (e.g., Archie, 1942;Poupon et al., 1970;Asquith, 1985).
Lithology and the mineral composition of the Mangaa sandstone were estimated using the M-N cross-plot, where the M and N values can be calculated using the following equations.
where Δt F and Δt M = the sonic log readings for the fluid and matrix, respectively, ρ F and ρ M = the fluid and matrix densities, respectively; and.
∅ NF and ∅ NM = the neutron porosity for both the fluid and matrix, respectively.
The density-neutron combination was used for calculating the total porosity (∅ T ) values, whereas the effective porosity (∅ e ) was calculated by eliminating the shale volume share from the total porosity.
where where ∅ N = the apparent neutron porosity, ∅ D = the density porosity, ρ ma = the matrix density, ρ b = the bulk density, and ρ f = the fluid density.
Besides, the effective porosity was estimated as the total porosity corrected for the shale volume (Schlumberger, 1989) as follows.
The shale volume was estimated using a single-log method (gamma-ray), where the shale parameters that were applied to the Early Pliocene and the Mangaa sandstone reservoir are the sonic transit time (Δt sh = 100 μs/ft), electric resistivity (R sh = 2.0 Ω.m), neutron porosity (∅ Nsh = 39%), and the bulk density of shale (ρ sh = 2.39 g/cm 3 ). Then, the corrected shale volume (V sh ) was estimated using Larionov's Eq. (1969).
where I sh = gamma-ray index which is estimated as using the Asquith and Gibson formula (1982) as follows. (1) where GR log = the gamma-ray log reading at a given depth, as the GR max and GR min are the maximum and the minimum gamma-ray readings, respectively, through the logged zones.
The water saturation (Sw) was estimated using the Indonesian equation of Poupon and Leveaux (1971).
where R t = the formation true resistivity, V sh = is the shale volume, R sh = shale resistivity, ∅ e = effective porosity, R w = formation water resistivity, m = the porosity exponent, and n = the saturation exponent.
The formation water resistivity (Rw in ohm.m) and the porosity exponent (m) were estimated using the Pickett cross-plot (Pickett, 1972) by plotting the porosity values versus their corresponding deep resistivity (R t ) at the measuring levels. Based on this plot, the Rw was calculated as 0.098 Ω.m at 45 °C formation temperature.
To determine the net-pay thickness, the reservoir and non-reservoir intervals were discriminated from each other considering the porosity, shale volume and water saturation cut-off values, where 15% effective porosity was applied as a cut-off value, whereas 35% and 75% were applied as cut-off values for the shale volume and the water saturation.

Seal capacity
For measuring the seal capacity of the cap rocks overlying the Early Pliocene and the Mangaa sandstone, the pore throat distribution of six samples was measured using the Mercury injection capillary pressure technique (MICP). These samples were selected representatively for the different seals above the probable reservoir horizons (EP-SD1, EP-SD2, EP-SD3, Mangaa 0, Mangaa 1 and Mangaa 2 at depths 1260, 1510, 1588, 1777, 1918, 2035 m). In addition, two samples were selected representatively for the fine grained muddy siltstone lithofacies within the reservoir units (EP-SD3 and Mangaa 0 at depths 1655 and 1837, respectively).
To measure the seal capacity, eight samples were measured using the Mercury intrusion capillary pressure (MICP) technique. Then, the air-mercury capillary pressure data that has been measured at the lab conditions (Pc L ) was converted to the oil-water and gas-water dual systems at the reservoir conditions (Pc R ). Then, the equivalent gas and the oil columns were estimated at four points of mercury saturation (0, 5.0, 7.5, and 10% mercury saturation). For estimating the hydrocarbon columns, the gas density = 0.19 g/cm 3 , oil density = 0.75 g/c m 3 , and water density = 1.05 g/c m 3 , gas/water surface tension = 50 dyn/cm, and the oil-water interfacial tension = 30 dyn/cm were considered. The seal capacity or the hydrocarbon column thickness or height above the free water or oil level was calculated at the aforementioned four mercury saturation points, each with its corresponding Pc R values using the seal capacity equations of Hannon (1987). It is also required to know a number of other variables, such as the water and hydrocarbon pressure gradients in the reservoir (Gw and GH, respectively), and the interfacial tension and contact angles at the reservoir and the laboratory conditions (Vavra et al., 1992). As these variables are not known with any confidence in this study, the input values for these variables were taken from another publication on the Mangaa sandstone sequence (Awatea-1 well; Murray and de Bock, 1996).
where Pc R is the capillary pressure taken at reservoir conditions, Hc is the hydrocarbon height/column in the reservoir below the seal rock (in m).
PcL is the threshold capillary pressure measured at the laboratory conditions, (σ.Cosθ) R is the interfacial tension multiplied by the contact angle between the hydrocarbons and the water at the reservoir conditions (gas-water = 50; oil-water = 26), (σ.Cosθ) L is interfacial tension multiplied by the contact angle between the air and mercury at the laboratory conditions (air-mercury = 367), G w is the water pressure gradient (0.469 psi/ft), and G H is the hydrocarbon pressure gradient (0.282 psi/m for oil, and 0.088 psi/m for gas. The MICP data were then presented on mercury saturation-pore radius size plot and its rank was estimated following the classification ranks of Nabawy et al. (2018a) (Nabawy and Wassif, 2017;El Sawy et al., 2020;Abuamarah and Nabawy, 2021;Radwan et al., 2022b) and the threshold pressure as the first inflection point on the capillary pressure curve that is followed by a large gradient increase on the incremental pore volume curve following Dewhurst et al. (2002).

Reservoir characterization using well log interpretation
Due to the gas shows during the well testing processes during the exploration activities in north of the Taranaki Basin, the reservoir characterization of the Early Pliocene and the Mangaa sandstone has been achieved using the well log analysis and interpretation. It aims at getting an insight into the hydrocarbon potentiality of the Mangaa reservoir and the probable Early Pliocene reservoirs by evaluating its petrophysical characteristics. As explained in the methodology section, to achieve a detailed reservoir characterization the shale parameters of the Mangaa reservoir should be estimated using the cross-plots of the gamma-ray versus the bulk density (ρ sh ), neutron porosity (∅ Nsh ), acoustic transit time (Δt sh ), and deep resistivity values (R sh ) ( Table 1). From the present study, the well log data of the Early Pliocene and the Mangaa sandstones were plotted as a function of depth (Figs. 3,4). Based on these plots, it is indicated that though the relatively high porosity values, the Early Pliocene sandstone sequences are composed primarily of siltstone and muddy siltstone with much clay content and no pay zones were assigned due to the high gamma-ray and shale content, and the low resistivity values (Fig. 3). It seems to be just gas pockets which may be evolved from the underneath Mangaa sandstone zones. So, the EP-SD sand zones were excluded from further petrophysical processing and interpretation.
Besides, the average petrophysical parameters of the Mangaa Formation, considering the reservoir intervals of the Albacore-1 well, indicate the presence of about 5.79 m net-pay thickness of 25.26% effective porosity, while the shale volume is 14.3% (Table 1). Although the relatively low shale volume, the API values are primarily higher than 75° API for the different reservoir zones, which is attributed to a relatively high K-feldspar content and intercalations with siltstone (Table 1) (Molenaar et al., 2015). It is indicated that the Mangaa sandstone attains an average gross thickness of about 342 m with average net reservoir equals 224.2 m, representing the potential thickness of the producing hydrocarbon, and net-pay thickness equals 5.79 m (Table 1).
Plotting the different petrophysical parameters of Mangaa Formation in the Albacore-1 well indicates that it can be subdivided into three reservoir units; Mangaa-0, Mangaa-1 and Mangaa-2 units. Based on this plot, the effective porosity values of the Mangaa Formation vary from 23.6% in Mangaa-2 sand unit to 27.7% in the Mangaa-0 sand unit. Besides, the shale volume values range from 12.3% in the Mangaa-2 sand unit to 16.9% in the Mangaa-1 sand unit The study also shows that the Mangaa-0 and Mangaa-2 units are not prospective, while the Mangaa-1 unit in the middle parts of the Mangaa reservoir sequence is prospective due to the presence of net-pay zone of 5.79 m thickness (at 2034.5-2044 m depth interval) with hydrocarbon saturation of about 25.6% ( Table 2). The detailed quantitative  Table 2.
To check the lithological composition and porosity of the studied reservoir sequence, the neutron-density cross-plot was applied. However, shifting the measuring points into higher porosity values on the dolomite line and the relatively high gamma-ray values (50-150 API, Fig. 5a) is primarily attributed to the relatively high content of feldspars and lithic fragments, the very fine sand composition, and the intercalations with many siltstone and mudstone interbeds. Also, following Teama and Nabawy (2016) this shift may be contributed to the effect of increasing the calcareous shale and silt content/interbeds.
The M-N cross-plot of the Mangaa reservoir shows that the predominant mineral composition of the Mangaa reservoir is siliciclastic that is represented by very fine argillaceous sand with some siltstone (Fig. 5b).
This achievement is stated by integration with the lithology description, petrography and the XRD data. This integrated data set declares that the measured log responses indicate a predominant very fine siliceous sandstone rich in feldspar and lithic content.
Plotting the porosity versus the water saturation values (Buckles technique, Buckles, 1965;Dewan, 1983 Table 2 Quantitative results of the petrophysical analysis of all the identified sand zones, within the Mangaa Formation and through the studied well NTS is the net sand thickness, NRT is the net reservoir thickness, NPT is the net-pay thickness, N/G is the net to gross, ∅e is the effective porosity, Vsh is the volume of shale, Sw is the water saturation and S hc is the hydrocarbon saturation If reservoir is at irreducible water saturation (Sw irr ) and exhibits similar bulk volume of water values BVW (BVW = ∅ x Sw) through its different levels, accordingly, the reservoir data points would either cross or be parallel to one of the hyperbolic lines of very low BVW value. If the data points failed to match these criteria, thus the reservoir is not at Sw irr (Asquith, 1985;Holmes et al.. 2009;Nabawy et al., 2018b).
It indicates the irreducible water saturation of the given reservoir and the prospectivity of the given reservoir intervals, where higher water saturation and porosity values indicate a lower grain size and reservoir quality, while lower water saturation and high porosity values indicate higher grain size and reservoir quality (Nabawy et al., 2018b;Abuhagaza et al., 2021;Safa et al., 2021).
Shifting from high porosity and low water saturation values (high permeability) to the low porosity-high water saturation values (low permeability), indicates shifting from high to low reservoir quality, from higher to lower permeability. Accordingly, the lowest BVW values indicate the high ability of the reservoir to produce free-water hydrocarbons, whereas the higher BVW values indicate fine sediments and low capability to produce free-water hydrocarbons (Fig. 6).

Microfacies analysis and mineral composition
The petrological and microfacies analysis of the Mangaa sand reservoir intervals and their discrimination into some microfacies were applied to reveal the mineral composition, the diagenetic features and the reservoir quality to support the well log data analysis. Each microfacies or rock type has its diagnostic features including the pore types, mineral composition, grain-to-grain fabric, and the diagenetic features that can be identified petrographically and considering the XRD data.

Feldspathic greywacke microfacies
This microfacies represents the unit EP-SD1 at depth 1436.5 m; it is composed of very fine, well sorted and angular to subangular monocrystalline quartz grains with some potash feldspars and a few lithic fragments, mica flakes, and fine to very fine glauconite pellets that are embedded in clay matrix (Fig. 7a). The lithic fragments are primarily composed of fine grained mudstone and siltstone fragments that often compacted and deformed (Fig. 7b). These muddy clastics seem to be the source of the argillaceous matrix due to their deformation and breakdown. Also, a few foraminifera tests are also observed, while the sparry calcite cement is rarely observed. Mica is represented by some deformed biotite and muscovite flakes (Fig. 7c). Porosity of this microfacies is fair to good (11.6-19.3%) and can be described as (1) intergranular pores (dyed blue, Fig. 7a), (2) micro vugs, and (3) matrix porosity. Porosity is mostly occluded by the clay content (Fig. 7d) and rarely by the micro sparry calcite.

Micaceous lithic arenite microfacies
This microfacies represents the unit EP-SD2 at depth 1545.3 m; it is composed of very fine, well sorted, and angular to subangular monocrystalline quartz grains with common lithic fragments and mica flakes. Besides, a few potash feldspars, plagioclase, and fine to very fine glauconite pellets are embedded in clay matrix (Fig. 7d). The plagioclase is fresh and sometimes is highly altered to sericite. The mica flakes are primarily bended, deformed and altered (Fig. 7e, f). The lithic fragments are frequently of sedimentary origin and composed of fine grained mudstone and siltstone (Fig. 7e). Besides, a few foraminifera tests are recorded but frequently crushed. These phenoclasts are embedded in clay matrix which is mostly related to the crushing and breakdown of the mudstone lithic fragments. It is rich in organic material. Porosity is primarily connected and ranked as good porosity (19.0%) and described as (1) intergranular pores, (2) micro to meso vugs (dyed blue, Fig. 7e, f), and (3) matrix porosity. Porosity is mostly deteriorated due to compaction and richness in the mica flakes and clay matrix that occluded the pore volume and bended with the predominated compaction forces (Fig. 7f).

Micaceous lithic greywacke/siltstone microfacies
This microfacies represents the unit EP-SD3, Mangaa-0, and the Mangaa-1 units at depths 1600.0 m, 1826.7 m, and 1996.5 m, respectively; it is composed of silt-sized to very fine, well sorted, and angular monocrystalline quartz grains with common mica flakes and lithic fragments (Fig. 7g). The lithic fragments are fine grained, and frequently altered mudstone fragments, whereas the mica is represented by highly compacted muscovite and biotite flakes that are frequently altered to sericite. Besides, plagioclase, and foraminifera bioclasts are less commonly with traces of altered glauconite pellets, recorded compacted and scattered in a clayey matrix (Fig. 7g). It is rich in organic material. This microfacies has good to very good effective porosity values (19.6-25.6%) and is described as (1) micro intergranular pores (dyed blue, Fig. 7g, h), and (2) matrix porosity slightly reduced by clays as patches and filling materials. Fig. 7 Photomicrographs of the Mangaa sandstone showing a intergranular connected pore spaces occluded by the clay matrix, PPL, feldspathic greywacke microfacies, 1436.5 m depth, b some altered lithic fragments of cloudy surfaces with very good effective porosity (dyed blue), PPL, feldspathic greywacke microfacies, 1517.0 m depth, c some fine to very fine altered lithic fragments and deformed mica flakes compacted together and embedded in matrix-rich groundmass, PPL, feldspathic greywacke microfacies, 1517.0 m depth, d very fine grained sandstone embedded in matrix-rich groundmass, XPL, feldspathic greywacke microfacies, 1436.5 m depth, e some deformed and altered lithics fragments, and mica flakes with a very fine foraminifera test, PPL, micaceous lithic arenite microfacies, 1545.3 m depth, f fine to coarse angular grained lithic and wood fragments broken down, bended and embedded in pore-occluding matrix, PPL, micaceous lithic arenite microfacies, 1545.3 m depth, g very fine to silt-sized quartz grains well-preserved pores partially filled with some mica flakes, organic materials and silt-sized foraminifera tests embedded in clay matrix, PPL, micaceous lithic greywacke/ siltstone microfacies, 1826.7 m depth, h silt-sized quartz grains and a foraminifera test in the center of the image with many connected intergranular pores partially occluded by clay matrix, PPL, micaceous lithic greywacke/siltstone microfacies, 1826.7 m depth, i fine to very fine highly altered muddy lithic fragments embedded in silica cement (micro quartz) and clay matrix, XPL, siliceous greywacke microfacies, 1925.2 m depth, j highly altered medium-grained lithic clast, altered and broken down to clay matrix, PPL, feldspathic and lithic greywacke microfacies, 2045.0 m depth, k fine to very fine volcanic fragments with much connected pore spaces, PPL, feldspathic and lithic greywacke microfacies, 2048.0 m depth, and l fine to very fine quartz sandstone rich in fresh glauconite pellets with much pore intergranular spaces, PPL, glauconitic lithic arenite microfacies, 2045.0 m depth ◂

Siliceous greywacke microfacies
This microfacies represents the Mangaa-1 unit at depth 1925.2 m. It consists of well sorted, very fine to fine, and sub-angular quartz grains, mostly monocrystalline, with some plagioclase, K-feldspars, mica flakes, and lithic fragments (Fig. 7i). The lithic fragments of this microfacies are different from that of the other microfacies and levels, where some fine fragments of metamorphic and igneous origin are recorded in addition to the muddy fragments.
These clasts are commonly altered to clays presented as patches and filling materials. The alteration of this microfacies is sever and extend even to the quartz grains that have pitted surfaces and corroded borders, sometimes presents as grain ghosts (Fig. 7i).
Porosity of this siliceous greywacke microfacies is very good and effective (21.6%), it is primarily represented (1) micro to meso intergranular pores slightly reduced by clays due to the alteration of the main siliciclasts (Fig. 7), and (2) matrix porosity in between the scattered clay matrix and silica cement.

Feldspathic and lithic greywacke microfacies
Downward, the Mangaa sandstone mineral composition became more variable in the Mangaa-2 unit with much more diagenesis, lithic and feldspathic content changes. The Mangaa-2 unit at depths 2045.0 m, and 2048 m is mostly represented by this microfacies which is built of well sorted, very fine to fine, and sub-angular quartz grains, predominantly monocrystalline. It is rich in plagioclase, K-feldspars, and lithic fragments (Fig. 7j). Phenoclasts are frequently altered with pitted surfaces and corroded borders. The lithic fragments are moderately altered and primarily of metamorphic and volcanic in origin ( Fig. 7k) with some sedimentary fragments. Alteration increases the amount of the clayey matrix which has been accumulated as patches and shale streaks due to the overload compaction. Mica flakes of this microfacies are less common than the other samples and are represented mostly by very fine muscovite and biotite flakes. Richness in clay content slightly reduced porosity of this microfacies to 17.6%. Pore types are described as (1) micro intergranular pores occluded by the clay matrix (Fig. 7j), and (2) matrix porosity.

Glauconitic lithic arenite microfacies
It selected from the Mangaa-2 unit at depth 2048.0 m, and composed of very fine to fine grained quartz grains and glauconite pellets with some lithic fragments (Fig. 7l). The Table 3 The bulk sample-XRD results for the Mangaa sandstone zones, Albacore-1 well quartz grains are well sorted and angular to sub-angular with clean surfaces. Also, the glauconite pellets are primarily fresh and scattered in between the quarts grains. Most components of this microfacies are fresh; sometimes with a slight alteration giving rise to scattered clay patches that masked few parts of the groundmass and filled the pore spaces. The lithic fragments of this sample are frequently muddy and sub-rounded to rounded due to the relatively distant transportation. Compaction mostly affected the clastic components of this sample, where the grains are in point and side contact situation (Fig. 7l). These are commonly altered to clays presented as patches and filling materials. The alteration of this microfacies is sever and extend even to the quartz grains that have pitted surfaces and corroded borders, sometimes presents as grain ghosts (Fig. 7i). It is characterized by very good porosity (23.3%) that is represented primarily by (1) well-preserved micro intergranular pores (Fig. 7l), and less frequently (2) matrix porosity.

XRD studies
The XRD analysis has been applied to confirm the mineral composition results that have been obtained from the petrographical study which indicated that the studied samples are primarily described as very fine sandstone and siltstone that consist of quartz and feldspar grains, in addition to the lithic fragments. This composition is confirmed by the XRD analysis data (Table 3), which shows a predominant quartz composition (38.3% at EP-SD2 unit up to 57.7% at Mangaa-2 sand unit), and some plagioclase and K-feldspar content (9.9-15.2% plagioclase and 2.7-6.3% K-feldspar, Table 3). Based on the XRD data, it is indicated that the studied sand samples are mostly rich mica/illite (16.0-31.8%) with some chlorite content (4.1-8.4%) and few kaolinite (0.4-0.9%) and illite/smectite (0-1.7%). The very low content of the calcite/ dolomite (< 3.1% for both, Table 3) is primarily due to the presence of few disintegrated foraminifera tests as indicated from the petrographical studies.
It is indicated that the quartz content increases in two cycles with the highest content encountered for the EP-SD3 unit (57.4%) and the Mangaa-2 unit (57.7%). This cyclic increase is associated with a cyclic decrease in the mica content as show in Table 3, where the lowest mica content is assigned at the EP-SD3 (17.7%) and Mangaa-1 and Mangaa-2 (16.0%).
Based on the XRD data, it is indicated that the petrographically described mudstone lithics are probably composed of silt-sized sandstone, and plagioclase beside to a chlorite and mica, where the clay minerals are not common (Table3). This is a main problem to be solved in the Pliocene sandstones throughout the New Zealand gas reservoirs, which make identifying the sandstones composition is difficult, where in spite of the clay content is relatively low, the GR API values are relatively moderate to high (50-150, Fig. 4). This may be explained by the attributing the API values to the relatively high feldspar content and the presence of intercalation with some thin shale streaks due to the overload stress compaction. However, it was difficult to distinguish authigenic clays from those either within the lithic clasts or in matrix related to the breakdown of the detrital lithic clasts.

Implications of diagenesis on reservoir quality
The diagenetic processes have either reservoir qualitydeteriorating or enhancing implementation on the reservoir rocks. The reservoir quality-deteriorating features can be discriminated into physical and chemical compaction, cementation, and predominance of authigenic minerals, while the quality enhancing factors can differentiated into dissolution and fracturing as stated in many literature (Nabawy and El Sharawy, 2015;El Sharawy and Nabawy, 2016;El Sawy et al. 2020;Nabawy et al. 2020b;Abuhagaza et al. 2021).

Reservoir quality-deteriorating features
In this concern, the Pliocene Mangaa sandstone is slightly affected by cementation due to the scarcity of the calcite/ dolomite content and the clay content as indicated from the petrographical study and the XRD data.
Illite and illite/smectite (rich in smectite) are the dominant clay minerals with significant chlorite and only minor kaolinite and traces of sepiolite. Illite, illite/smectite and chlorite are consistent with the phases observed in some of the fine grained lithic clasts (mudstones, siltstones) and also in the degraded plagioclase ghosts. Some minor chlorite may be authigenic and related to the breakdown of biotite mica. There may also be some authigenic component to all of the clay phases, but none is volumetrically significant. It would appear that the compaction of labile clasts has redistributed the clay content into the matrix and intercalated streaks, but there is a little evidence for significant authigenic reprecipitation of clays.
Minor authigenic carbonate is seen in most samples, but also it is volumetrically insignificant and never has low ability to occlude pore spaces. Much calcite is probably in the form of deteriorated foraminifera. No clear diagenetic siderite has been observed, and it seems likely that siderite is present in the altered mudstone clasts. At least some are stained for ferroan calcite, although slightly ferroan dolomite was identified from the XRD data.
Minor authigenic pyrite content is present in most samples but with negligible amounts as confirmed from the XRD data, and is likely related to the biotite breakdown and other labile phases. No evidence was seen for quartz overgrowth.
Since the present burial depth is 1.4-2.0 km, so it was thought that the compaction to has slightly affected the present samples. However, the Mangaa sandstones have undergone a moderate amount of mechanical compaction as indicated from bending and fracturing of the labile mica flakes and the lithic clasts, which have broken down forming the matrix as indicated in the feldspathic greywacke microfacies (Fig. 7c), micaceous lithic arenite microfacies (Fig. 7e, f), and the micaceous lithic greywacke/siltstone microfacies (Fig. 7g). Compaction is also indicated from the wavy extinction of the quart grains as revealed in the micaceous lithic arenite microfacies (Fig. 7f), and the siliceous greywacke microfacies (Fig. 7i). The point contact and long straight grain-to-grain contact confirms a slight mechanical compaction of the feldspathic greywacke microfacies (Fig. 7d), and the glauconitic lithic arenite microfacies (Fig. 7l). Besides, it is worthy mentioned that the Mangaa sequence is dissected by many gravity-driven listric growth faults.

Reservoir quality-enhancing features
Fracturing and dissolution are the common porosity and reservoir quality-enhancing factors in different reservoir rocks (Nabawy and El Aal, 2019;Nabawy et al. 2020a;2020b). For the present study, dissolution is the most effective diagenetic feature. It seems to compensate the deteriorating implication of compaction on quartz grains, where in spite of the compaction effect, the quartz grains are still float within the matrix and just in grain-to-grain point and straight contact ( Fig. 7a:l). Dissolution is indicated by the relatively high visual porosity values (intergranular and vuggy pore spaces, 11.6-25.6%) and scarcity of cement in the fledspathic greywacke, micaceous lithic arenite, micaceous lithic greywacke/siltstone, siliceous greywacke, feldspathic and lithic greywacke microfacies (Figs. 7a,b,e,f,h,j,k). The corroded borders and pitted surface of the quartz grains are additional confirmation for the dissolution effect that developed on the grain borders and surfaces (Fig. 7i). The scarcity of the clay content as indicated from the XRD (Table 3) though the richness in deteriorated feldspars may be additional indicator that confirms leaching out of the authigenic clays that produced from this feldspars alteration.
Fracturing, on the other hand, is not clearly evident from the petrography due to the absence of cement and dominance of ductile flakes and muddy lithics. This may be attributed to the impact of dissolution and leaching out of the broken constituents. However, implementation of fracturing may be indicated by the presence of some broken wood fragments as in the micaceous lithic arenite microfacies (Fig. 7f), fragmented lithics as in the feldspathic and lithic greywacke microfacies (Fig. 7j), and deteriorated foraminifera tests in the micaceous lithic greywacke/siltstone microfacies (Fig. 7g, h).

Seal capacity results
The seal capacity for the present study has been achieved via measuring the pore throat distribution and its threshold pressure using the MICP technique. Samples were primarily selected representatively for the gas-bearing zones and from seals above the gas-bearing zones at the depths indicated in Table 4 and then their threshold pressures were estimated and their hydrocarbon-sealing efficiency were estimated (Table 4). These tested seals are mostly composed of siltstones and muddy siltstone. The grain Table 4 Summary of the MICP results of the studied seal samples and the potential hydrocarbon columns/thickness which the seals can retain before leakage PcL is the threshold pressure measured in the lab, R 35 is the pore radius measured at 5% mercury saturation, k is the permeability and R mode is the more dominant pore radius  (Table 3). In the EP-SD2, these constituents reach up to 31.8%, 6.9%, and 1.3% for the mica, chlorite and siderite, respectively (Table 3). The mica content which has a density equals to 2.88 g/cm 3 , the chlorite content of density may reach up to 3.3 g/cm 3 , and the siderite content which has a density equals to 3.96 g/cm 3 due to its iron content. Based on the mercury saturation-pore diameter plot (Fig. 8a), it is indicated that less than 15% of the pore throats are in the range of meso and macro pore sizes (1.0-100 μm), while the predominant pore spaces are in the range of micro and nano pore sizes (D < 1.0 μm). Besides, presenting the pore   Fig. 8b). From this plot, it is also indicated that the R 35 (mercury saturation at 35%) and the R mode values are in the range of the micro pore spaces with less values for the EP-SD3 mudstone, seal above Mangaa-0 and the seal above Mangaa-1 (at depths 1655, 1777, and 1918, respectively), i.e., higher efficiency for trapping oil and gas-bearing zones with good thickness (Table 4). This high efficiency is primarily attributed to the siltstone and muddy siltstone composition.
Following the reservoir quality classification of Nabawy et al. (2018b), these samples are tight and may act as seal beds. However, this plot indicates the presence of 5-15% pore spaces in the size of macro and meso pores which reduces the efficiency of these seal samples. The R 35 and the R mode values of the seal above the Mangaa-0 sand in the EP-SD3 and the mudstone sample of the EP-SD3 to be the best seals in this sequence (Table 4).
The estimated threshold pressure (PcL) for each sample can be used by importing into Eqs. 9 and 10 to estimate the efficiency of these samples to act as seal rocks and to estimate thickness of the gas and oil that can trap. The results of the calculations suggest that the potential hydrocarbon columns, that the sampled seal intervals in Albacore-1 well could hold 17.4 m up to 43.0 m hydrocarbon columns (Table 4, Fig. 9).
Correlating the estimated seal efficiencies (hydrocarbon thicknesses) with that published by Radwan et al. (2022b) indicated that the efficiency of seals in the Mangaa sandstone sequence to the north of Taranki Basin (at the Albacore-1 well) is less effective than that in the central parts of the basin at Karewa Field (reached up to 150 m and 237 m of gas and oil thickness at 10% Hg saturation, respectively).
To examine the ability of the studied seals to leak the fluids, their permeability values were estimated based on the R 35 using the Winland Eq. (1972) and a verified mathematical model has been estimated to calculate permeability of these seals as a function of their porosity as follows. k seal = 0.042 (∅ seal ) 0.69 (R 2 = 0.614) (11). Following the rock classification ranks of Abuamarah et al. (2019), the obtained permeability values (0.327-0.412 md, Table 4, Fig. 10) indicate poor permeable rocks (0.1 md < poor permeable rock < 1.0 md). Following this classification good seals are tight permeable rocks (0.01 md < tight permeable rock < 0.1 md), whereas the best seals are in the micro Darcy range (very tight permeable rock < 0.01 md). Therefore, based on the estimated permeability values of the EP-SD and Mangaa seals, the studied seals are in the range of poor permeability and considered of low potentiality.

Conclusions
The well log analysis of the studied Miocene-Early Pliocene sequence in north the Taranaki Basin indicates that the Early pliocene sandstone (EP-SD) units are not prospective, while a total of 5.79 m net-pay thickness is assigned in the Upper Miocene-Lower Pliocene Mangaa-1 sandstone unit. This prospective unit has an effective porosity equals 19.7%, shale volume equals 19.8% and hydrocarbon saturation about 25.6%. Besides, the Buckles plot indicates that, in spite of the relatively low reservoir quality in north Taranaki Basin, the BVW values > 0.18 for about 75% of the total estimated points. This indicates that the Mangaa sandstone units will primarily produce water with low hydrocarbons production.
The petrographical studies indicate that the studied Late Miocene-Early Pliocene zones are composed of very fine greywackes and lithic arenites/siltstones rich in feldspars, chlorite and mica. The mechanical compaction is the predominant reservoir quality-deteriorating factor. However, dissolution and leaching of the dissolved and altered components reduced the compaction implementation on the reservoir quality of the Mangaa sandstone units.
The (MICP) results of the seal beds above the different sand units indicate the presence of relatively low threshold pressures that can potentially hold only moderate hydrocarbon columns (not more than 43.0 m), which would imply poor to moderate seals. These low values are to be expected given that most of the potential seal units are siltstones and muddy siltstones but neither mudstones nor evaporites.
Therefore, the present study reveals that only a Mangaa-1 sandstone unit is considered prospective unit (5.79 m) with micaceous and lithic arenite and greywacke composition with some chlorite and feldspars content. Its hydrocarbon charge is sealed by the muddy siltstone of the overlying Mangaa-0.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.