Characterization of Pliocene Biogenic Gas Reservoirs from the Western Black Sea Shelf (Romanian Offshore) by Integration of Well Logs and Core Data

: The successful interpretation of open-hole well logging data relies on jointly using all available petrophysical and geological information. This paper presents relevant case studies related to the integration of well logs with core measurements for exploration wells drilled in the Romanian continental shelf area of the Western Black Sea basin. The analyzed wells targeted gas-bearing sands and silts complexes of Early Pliocene (Dacian) age, developed in a deltaic to shallow marine sedimentary environment in two distinct ﬁelds. The wireline logging programs included conventional formation evaluation logs, pressure surveys, nuclear magnetic resonance, and borehole electrical imaging logs. The core dataset comprised routine and special measurements (porosity, grain density, permeability, water saturation, and Archie parameters) carried out at quasi-reservoir conﬁning pressure. The wireline logging suites were interpreted via a deterministic workﬂow, including core-derived interpretation parameters. Other core-derived parameters were used for constraining and validating the log interpretations. The results show that a problem related to the ambiguity of formation water resistivity can be overcome through resistivity–porosity dependencies constructed to include potential aquifer zones in the proximity of the Dacian gas-bearing reservoirs. This study also revealed and quantiﬁed uncertainties regarding the estimation of gas–water contacts from formation pressure surveys, which can be mitigated by the conﬁrmation or correction of pressure-derived ﬂuid contacts via the well log interpretation results. Lastly, we identiﬁed a probable resistivity logs suppression effect related both to high contents of capillary-bound water and also to the limited resolution of electrical logging tools in the presence of sand-shale thin bedding or laminations.


Introduction
Considered as one of the most important hydrocarbon-bearing areas in SE Europe, the Western Black Sea basin has already demonstrated its potential through oil and gas fields discoveries in the continental shelf areas of Romania, Bulgaria, and Ukraine.Detailed studies regarding the geology, tectonics, hydrocarbon systems, and hydrocarbon plays from the Western Black Sea basin, including the Romanian offshore, have been carried out by Robinson et al. [1,2], , Dinu et al. [6], Bega and Ionescu [7], Crânganu et al. [8], Georgiev [9], Tari et al. [10], Nikishin et al. [11,12], Oaie et al. [13], and Boote [14].
In the Romanian Black Sea shelf sector, the most important hydrocarbon fields are located in the central area, the Histria Basin/Depression (Figures 1 and 2).Ample geophysical investigations have been performed in this area starting in the 1970s, by means of seismic, gravity, and magnetometric surveys, for the identification of favorable geological structures and hydrocarbon accumulations.Nevertheless, the volume of geophysical exploration and drilling activities in shallow water and deepwater perimeters is still limited.To date, oil was encountered mainly in Cretaceous, Eocene, and Oligocene formations.
Commercial gas accumulations or gas shows were identified in Cretaceous and Eocene formations but also in Late Miocene-Pliocene deposits (Pontian and Dacian stages) in the majority of the wells drilled [4,5].(adapted from [4,5,13,15]).

2.
Figure 2. A-A' Geological cross-section with W-E direction through the Romanian Black Sea shelf area, based on seismic reflection and well data (redrawn and modified from [9]).
The Early Pliocene (Dacian stage) dry gas discoveries analyzed in this study, hereafter denoted "A" and "B" fields, are situated in the Histria Basin/Depression area of the Romanian offshore.The gas charge is exclusively biogenic, being reservoired in marginal marine (deltaic) to shallow marine sands and silts complexes.Both fields are four-way dip closures lying within a NE-SW trending fault terrace, bounded to the NW by an up-thrown fault and to the SE by a down-thrown fault, and are less than 20 km apart.
The petrophysical evaluation of "A" and "B" gas fields has presented several challenges since their discovery.The reservoir sands are unconsolidated, fine to very fine grained, becoming progressively more silty and muddy downwards.All available water samples collected in the exploration wells were contaminated by drilling mud filtrate (high contents of potassium detected in the samples) and, consequently, were not representative for quantitative log interpretation.Additionally, within the reservoir sections crossed by the wells, there are no distinct clean water-bearing sands from which to derive an unambiguous water resistivity determination.The main problems encountered in previous studies were the adequate estimation of water and gas saturations and the identification of gas-water contacts, which are vital inputs for a realistic gas reserves evaluation.
Few studies have been published on topics related to geophysical well logging and formation evaluation for exploration wells drilled in the Romanian offshore area of Western Black Sea [8, [15][16][17].This paper attempts to address this knowledge gap by reporting and discussing particular issues regarding the well log interpretation (including Nuclear Magnetic Resonance data) and its integration with routine and special core measurements, for several exploration wells that intercepted the Dacian gas-bearing reservoirs.The geological and tectonic background of the region is also presented, and the wireline logging programs are analyzed in correlation with the petrophysical interpretation results.The specific approaches used in this research to mitigate or overcome the uncertainties related to fluid saturations and fluid contacts estimation provide a suitable reservoir characterization methodology to be used for other hydrocarbon discoveries in the Black Sea shelf area.

Geological Setting and Hydrocarbon Systems
The Black Sea is generally considered to be a back-arc extensional basin with active rifting starting at the end of Early Cretaceous [1,11,12,18].The geological units from the Romanian offshore area are continuations towards east of the main continental structural units of the Dobrogea territory: the Moesian Platform (South and Central Dobrogea), the North Dobrogean Orogen, and the Scythian Platform [19][20][21].The present structure of the Western Black Sea basin, in the Romanian sector, is a result of tectonic movements along major crustal faults with NW-SE (or WNW-ESE) strike and which extend towards the East Carpathians Bend Zone: Sfântu Gheorghe fault, Peceneaga-Camena fault, Capidava-Ovidiu fault, and Intramoesian fault (Figures 1 and 2).
Besides the continuation of these major crustal faults into the shelf, secondary NW-SE faults from the same system have been delineated through seismic surveys.Another system consists of normal or strike-slip faults oriented approximately parallel to the Black Sea coast, such as Razelm fault, Lacul Roşu fault, and West Midia fault [13].
To date, the Histria Basin/Depression area of the Romanian Black Sea shelf is the most important from the standpoint of hydrocarbons potential and discoveries.This posttectonic depressionary area is superimposed over the North Dobrogea Orogen and it is limited by the Heracleea fault to the north and by the Peceneaga-Camena major crustal fault to the south.The border of the basin is marked by a structural feature called "Euxinic Threshold/Edge", which could correspond to the limit of shelf deposits in the Late Cretaceous.The syn-rift and post-rift sedimentary deposits of Histria Basin include Cretaceous, Paleogene, Neogene, and Quaternary formations (Figure 2).
The Neogene age package (up to 5 km thickness) is dominated by fluvio-deltaic and marginal marine clastics.The structural style of these offshore formations generally resulted from gravity-driven tectonics.Deposition was controlled both by sediment supply via paleo-rivers (e.g., proto-Danube) and by sea-level changes [7].Of particular interest in the context of this study is the Late Miocene-Pliocene deposits, which include claystones, mudstones/shales, siltstones, sandstones, and sands.
According to Moroşanu [5], four deep thermogenic hydrocarbon systems and one shallow biogenic gas system are present in the Romanian continental shelf area.For the biogenic system, the source rocks are considered the Middle-Late Miocene (Sarmatian-Pontian) pelitic deposits, whereas the seals are provided by the pelitic intervals of the Pliocene.The gas is reservoired in Late Miocene-Early Pliocene sandstones and sands.The traps are structural (drape anticlines, roll-over anticlines, beds that have undergone listric faulting) or stratigraphic (pinch-outs and submarine fans) [5].Besides the offshore shallow gas accumulations, gas seepage phenomena (methane chimneys) are a common occurrence within the Pliocene-Quaternary sediments.
The main gas-bearing formation for the "A" and "B" fields is an Early Pliocene (Dacian stage) shallow marine sequence, which lies above Late Miocene (Pontian stage) deposits.The reservoirs comprise predominantly fine to very fine-grained, occasionally thinly bedded, muddy-silty sands (frequently micaceous and with carbonaceous material), the top seal being provided by a mudstone unit, which may represent a flooding surface.
The likely depositional environment is represented by a shallow marine (lower shoreface-upper offshore) regime for field "A" to a marginal marine (deltaic) setting for field "B", with sediment influx provided from a delta system located towards northwest or north.In industry terminology, the reservoirs were conventionally divided into a top "Sand" unit (thicker and better-quality sands) and an underlying "Silt" unit (lower quality silt-dominated facies).This facies separation is clearer in field "A" than in field "B", where the sands are thinner and inter-bedded with silt units.The thickness of the "Sand" unit ranges from 26-27 m in field "A" to 15-28 m in field "B".The "Silt" unit has a thickness ranging from 38-41 m in field "A" to 14-47 m in field "B".

Well Logging Programs and Core Measurements Data
This study uses data from six exploration and appraisal wells (denoted A-1, A-2, B-1, B-2, B-3, and B-4), which crossed the Early Pliocene reservoirs of "A" and "B" fields, with the vertical depth of the wells varying from 1250 to 3000 m.The wells were drilled with water-based KCl mud and the borehole size was 8.5 inches in all the reservoir/final sections.
The openhole wireline investigation program run in the final sections of the newer wells included a logging suite of PEX-Platform Express type (Schlumberger Ltd., Houston, TX, USA), comprising HALS-High-Resolution Azimuthal Laterolog, TLD-Three-Detector Lithology-Density, MCFL-MicroCylindrically Focused Log, and HGNS-Highly Integrated Gamma Ray Neutron Sonde tools (Figures 3 and 4).The coring program undertaken in the wells from the "A" and "B" fields and the petrophysical measurements carried out are presented in Table 1.Core to log depth shifts were applied to optimize the match between the core lithology description and the log responses.

4.
The routine core analyses (RCAL) included water saturation (S w ) determination by the Dean-Stark extraction method [22], gas permeability (K g ) measurements with nitrogen permeameters, Klinkenberg-corrected gas permeability (K k ) [22,23], helium gas expansion porosimetry (φ), and matrix/grain density (ρ ma ) determinations, performed on plugs and reported at a net overburden/confining pressure of 1400 psig (pounds per square inch gauge) ≈ 98 bar, considered representative.This value was selected based on analyses of porosity and gas permeability variations as functions of overburden pressure, raised in steps from 200 to 2000 psig (15-139 bar).The largest decrease from ambient condition values for both porosity and permeability occurred between 200 and 1000 psig (15-70 bar), with permeability showing the greatest sensitivity.From 1000 to 2000 psig (70-139 bar) the rate of decline was low, with a discernible break of slope at 1400 psig, in most cases.For two wells (A-1 and B-3), a special core analysis study (SCAL) was performed on cut plugs.Besides K g , K k , φ, ρ ma determinations and X-ray diffraction (XRD) composition analyses, electrical measurements were conducted on brine-saturated plugs at overburden pressures increased incrementally from 400 to 1200 psig (29-84 bar).The saturated plugs underwent the determination of formation resistivity factor (F), cementation exponent m, resistivity index (I R ), and saturation exponent n [24,25], expressed by: where R o [Ω m] is the resistivity of a water-saturated core plug, R w [Ω m] is the formation water (brine) resistivity, φ [V/V] is the fractional porosity, a is the tortuosity factor (set to 1.0), m is the cementation exponent, R t [Ω m] is the resistivity of a core plug partially water saturated, S w [V/V] is the fractional water saturation, and n is the saturation exponent.The Archie coefficients m and n were calculated via linear least-squares regressions from F = f (φ) and I R = f (S w ) dependencies.Table 2 presents the main RCAL results obtained for the reservoir sections of three wells from "A" and "B" fields, at 1400 psig (≈ 98 bar) confining pressure.The relatively high density of the matrix, especially for wells B-2 and B-3, can be related to the calcareous and micaceous character of the reservoir sands.The higher permeabilities obtained for well B-3 might be caused by the poorly consolidated core plug material.Table 3 shows the SCAL results obtained for the "Sand" reservoir intervals of wells A-1 and B-3, at 1200 psig (84 bar) confining pressure.The mean values of Archie's cementation exponent and saturation exponent measured on core plugs were m = 1.71, n = 1.67 for well A-1 and m = 1.54, n = 1.32 (unusually low saturation exponent) for well B-3.
Table 3. Results of special core measurements performed at quasi-reservoir confining pressure for A-1 and B-3 wells-"Sand" intervals.

Well
Depth Wireline formation pressure measurements (p [psia]-pounds per square inch absolute) were acquired in all the wells over various depth intervals, which included the "Sand" and "Silt" units, but sometimes without or with few valid pressure readings within the reservoirs (Table 1).The formation pressure measurements taken within the reservoir sections were in the range of 1686-1720 psia (116-119 bar) for field "A" and 1634-1661 psia (113-115 bar) for field "B" and the recorded reservoir temperatures were in the 34-39 • C range.
The gas samples collected during the wireline testing or drill stem tests in wells from both fields show the same composition, with very high methane content (>99.7%),traces of N 2 and CO 2 , without H 2 S, and with a gas gravity of 0.557 relative to air.According to the ideal gas law, at the reservoirs' pressure and temperature conditions, the gas density should be ≈0.08 g/cm 3 .
The X-ray diffraction (XRD) analysis carried out on the 12 m core extracted from the "Sand" unit of well A-1 indicated a total clay-mineral content of 12.8-53.0%by weight (mean content: 28.7%).The mean relative abundance of the identified clay minerals was: 35.6% mixed-layer illite/smectite, 31.6%illite and mica, 29.8% chlorite, and 3.0% kaolinite.Such clay-mineral contents do not allow the reservoirs to be considered clean for petrophysical evaluation purposes.

Pre-Interpretation/Preliminary Processing of Well Logs
No representative formation water samples were obtained from the "Sand" and "Silt" units of the wells drilled in "A" and "B" fields, hence R w can only be inferred from indirect sources.Furthermore, within the reservoir units, there are no obvious clean and waterbearing sand beds that could provide a reliable R w reference.To overcome this problem, resistivity-porosity dependencies, such as the Hingle crossplot [26,27]-Equation (3) and the Pickett crossplot [27,28]-Equation (4), were employed over slightly larger Pliocene depth intervals including the reservoirs and post-reservoir sections in their proximity, where clean and potentially water-bearing sands are encountered: Applicable to clean (non-shaly) formations, these dependencies allow an estimation of R w at formation temperature, matrix parameters (ρ ma , φ Nma , and ∆t ma ), and m, especially when the reservoirs show large enough porosity and resistivity ranges.The crossplots are constructed with R t as a function of a computed φ curve or as a function of measured ρ b , φ N or ∆t c .Water-bearing formations are identified as distinctive linear trends of ) data and a linear least-squares regression (the S w = 1 "water line") through the clean data points provides the unknown parameters.The deep investigation Laterolog R A curves were used as a substitute for R t in Equations ( 3) and ( 4) because in high porosity reservoirs, the mud filtrate invasion is very shallow and R A invasion corrections are minor.For the R t ˆ−(1/m) = f (φ) or log(φ) = f (log(R t )) variants of the crossplots used in the pre-interpretation phase, φ was derived from density logs (φ D [V/V] is the density porosity) using the mean ρ ma values measured on core plugs and pore fluid (mud filtrate) densities ρ mf = 1.02-1.06g/cm 3 depending on filtrate's salinity: The NMR data were acquired only in A-2 well, over the 720-1590 m depth interval (pre-reservoir, "Sand" and "Silt" units, and post-reservoir deposits).The tool measures the buildup and decay of the polarization of hydrogen nuclei (protons) in the pore fluids, the measurements being lithology independent.The total (φ t , PHIT or TCMR) and effective (φ e or PHIE) porosities can be unequivocally distinguished, and it is possible to discriminate between free reservoir fluids (mobile water and hydrocarbons) and clay-bound or capillarybound water and, also, to assess the pore space distribution [29].The recorded T 2 [ms] (transverse relaxation time) NMR array data are a measure of the pore surface to fluid volume ratio (S/V); for small pores, S/V is large and corresponds to short T 2 values (rapid signal decay), whereas for large pores, S/V is small and corresponds to long T 2 values (slow decay rate of the signal).
The T 2 distribution in the 0.3-3000 ms range was calibrated using the input total porosity curve TCMR in terms of porosity fractions, for the determination of total bound fluid (water) fraction (BFT), clay-bound water fraction (CBW), capillary-bound water fraction (BVI-Bound Volume Irreducible), and producible porosity (FFI-Free Fluid Index): PH IT (TCMR) = FFI + BFT (7) The CBW and BFT fractions were defined using T 2 cut-offs of 3 ms and 33 ms (suitable for clastics), respectively, [29].The formation permeability K NMR was also estimated from NMR data, using the Coates model [29,30]: where K NMR is expressed in [mD]; PHIT, FFI, and BVI are expressed in [%]; and a, b, c are statistical parameters derived from permeability measurements on cores.When no such measurements are available to perform a calibration via polynomial regression, as in the case of the A-2 well, default empirical values can be assigned to the coefficients (a = 4, b = 2, c = 10) [30].

Quantitative Interpretation of Well Logs
A deterministic quantitative interpretation workflow was adopted, based on the combined responses of the density and neutron logs (Figure 5).For this ρ b = f (φ N ) log combination, the hydrocarbon (gas) effects can be accurately corrected and taken into account and φ e is obtained after an iterative and convergent hydrocarbon correction, at each depth level.Unlike sonic logs, neutron and density responses are not affected by the lack of compaction characteristic for young and shallow-depth sands, like the Pliocene reservoirs from "A" and "B" fields.
The "shaly-sands" model used (Figure 6) consists of clean matrix (V ma ) + wet clay (V clay ) + effective porosity (φ e ), being defined in terms of fractional volumes as: where V CBW is the volume of clay-bound water and φ clay is the wet clay porosity.A clear distinction is made between the terms "clay" (fraction composed of dry clay minerals and clay-bound water) and "shale" (in the sense of sedimentary rock consisting of a clay fraction and silt-sized particles).The ratio of clay and silt-sized particles in shales is highly variable and one can define a clay:shale ratio (CSR = V clay /V shale ) or, alternatively, a silt index (I SILT = V silt /V shale ).Usually, the clay content in shales may be 40-80%, the silt fraction (V silt ) being treated as part of the clean matrix and having similar properties [31].
6 The general responses of the density and neutron logs are expressed as: where ρ clay [g/cm 3 ] is the wet clay density, ρ hApp [g/cm 3 ] is the hydrocarbon (gas) apparent density, φ Nclay [V/V] is the wet clay neutron porosity, The S xo and S w water saturations were evaluated in the effective porosity system using the "Indonesia" model [34]: where R mf [Ω m] and R w [Ω m] are the mud filtrate resistivity and the formation water resistivity at reservoir temperature and R clay [Ω m] is the wet clay resistivity.These expressions, which account for the excess conductivity due to clays, are well suited for water saturation evaluation in formations with high clay content and low-salinity formation waters, as expected in the "A" and "B" fields.In Equations ( 14) and ( 15), R xo and R t were approximated by the MCFL or MLL R A curves and the deepest-reading Laterolog R A curves, respectively.
Because φ e and S xo cannot be determined independently in gas-bearing intervals, an iterative correction is performed using Equations ( 12)- (14) until φ e and S xo convergence, then S w is derived from Equation (15) and the converged φ e solution.
To provide an independent and continuous V clay estimation as input in Equations ( 12) and ( 13), GR logs were used as main clay indicators.Several studies [31, 35,36] have shown that the GR response in shaly formations should increase linearly with the clay volume fraction.The linear V clayGR [V/V] clay volume estimator is defined via a gamma-ray index 27,37]: where GR clean [gAPI] is the gamma-ray background radioactivity of clean reservoir rocks and GR clay [gAPI] is the gamma-ray radioactivity of clays.To verify that the relatively low-contrast GR logs recorded in the analyzed wells can be used as valid clay volume estimators, V clayGR were checked against independently derived neutron-density V clayND [V/V] estimations: where φ N is the matrix-corrected neutron apparent porosity (referenced to sandstone lithology), φ D is given by Equation ( 5   The wet clay parameters (ρ clay , φ Nclay, R clay ) were statistically derived from ρ b , φ N , and R t histograms and GR = f (ρ b ), GR = f (φ N ), and GR = f (R t ) frequency crossplots, as the most repeatable readings in the shaliest intervals delineated through a GR > GR cut-off criterion.The mean ρ ma, m, and n values measured on cores at quasi-reservoir overburden pressure were used for all the wells, along with the R w values estimated in the pre-interpretation phase and confirmed or slightly refined during the interpretation.Gas density (ρ h ) was estimated considering the true vertical depths (TVDs) of the reservoirs, the hydrostatic pressures, and the measured temperatures.Table 4 synthesizes the main log-derived and core-derived interpretation parameters.
Table 4. Main log-derived and core-derived parameters used for the effective porosity (φ e ) and water saturations (S xo , S w ) evaluation in the reservoir intervals of "A" and "B" fields.

Well
ρ ma [g/cm 3 ] R mf @ T mf -Mud filtrate resistivity at surface measurement temperature.

Wireline Formation Pressure Data Processing and Analysis
In favorable conditions, wireline formation pressure data allow the identification of the reservoir fluid contacts (gas-water contact-GWC, oil-water contact-OWC and gas-oil contact-GOC), the classification of fluid types (in situ densities), and, also, the assessment of fluids separation via permeability barriers.
The depths of the fluid contacts are critical for calculating the volume of hydrocarbons in the reservoirs and important for increasing the accuracy of well log interpretations, such as porosity evaluation in the case of reservoirs with multiple fluids and different physical properties.The available methods for locating the fluid contacts include fluid sampling (from drill stem tests or wireline formation testers), analysis of log-derived or core-derived fluid saturations variation with depth, and pressure profile surveys.
The normal (hydrostatic) pore pressure p measured at a particular depth z = TVDSS (true vertical depth subsea) is given by: where ρ f is the fluid density and g is the acceleration of gravity (≈9.81 m/s 2 ).If the pressures p 1 and p 2 are measured at different depths z 1 and z 2 , one can define a pressure-depth gradient (trend) ∆p/∆z, where ∆p = p 2 − p 1 and ∆z = z 2 − z 1 .
If a single fluid is present in the pore space of the formations, successive pressure measurements carried out at various depths will define a single pressure-depth gradient.
For multiple fluids with significantly different densities, a p = f (z) plot will show changes in the slope corresponding to multiple pressure-depth gradients, e.g., (∆p/∆z) 1 and (∆p/∆z) 2 .
The depths where the pressure trends are changing define the position of the fluid contacts (Figure 8) and the density of the reservoir fluids can be inferred from the gradients (e.g., ρ f1 = (∆p/∆z) 1 /g and ρ f2 = (∆p/∆z) 2 /g).Regression analysis was performed using the wireline pressure datasets of the wells from "A" and "B" fields, for deriving meaningful pressure-depth trends related to reservoir fluids type and their in situ densities.The trends were subsequently used to locate the likely GWCs, which were compared with the results of quantitative well log interpretation and of NMR data processing.

Permeability Modeling
The analysis of core data from "A" and "B" fields shows a strong correlation between K k and φ (Figure 9a), indicating that porosity is the main factor controlling the formations permeability.The inverse correlation between permeability and shaliness (i.e., permeability decreases with increasing clay content) is illustrated via well logging data by the K NMR = f (GR) dependency from Figure 9b, which also shows the inverse correlation between K NMR and S w .
For all the wells from "A" and "B" fields, continuous permeability curves were modeled via a multiple linear regression (MLR) technique, by fitting in the least-squares sense exponential functions of the type K = f (φ e , V clay ) or K = f (φ e , V clay , S w ).The predicted permeability curves that were tested had the form: where a, b, c, and d are computed polynomial regression coefficients that allow the best possible reconstruction of the input control (reference) permeability datasets.The control data for the best fit of permeability functions was provided by: (a) the sets of depth-shifted and log-matched permeability measurements on core plugs available for the A-1, B-1, B-2, and B-3 wells; (b) the NMR-derived formation permeability estimated for the A-2 well.The permeability modeling was carried out both for each well and also at the field level, by jointly fitting synthetic K curves to the multiple discrete-continuous (coreand NMR-derived) or solely discrete (core-derived) control datasets from fields "A" and "B", respectively.(a) Core permeability-core porosity relationship for the pre-reservoir, "Sand", "Silt", and post-reservoir deposits of "A" and "B" fields; (b) NMR permeability-gamma ray intensity dependence for the same sequence of deposits in the A-2 well (567 data levels).The color code corresponds to S w values resulting from the petrophysical interpretation of well logs.

Results and Discussion
The Early Pliocene gas-bearing reservoirs from the "A" and "B" fields (Figures 3 and 4) appear on well logging data with a characteristic signature, i.e., a coarsening upward and suppressed gamma-ray configuration (limited contrast between the reservoir and nonreservoir intervals), few intervals with neutron-density gas crossover as a consequence of clay content, and locally suppressed resistivity response but with a clear increased sonic ∆t c response (low P-wave velocity).There is a large discrepancy between the maximum recorded deep R A resistivities in the "Sand" main reservoir unit from the two fields: A-1 well-170 Ω m, A-2 well-371 Ω m, B-1 well-16 Ω m, B-2 well-17 Ω m, B-3 well-26 Ω m, B-4 well-82 Ω m.The likely underestimated R A readings may be explained by an averaging effect upon the Laterolog tools, when investigating thin beds with a thickness less than the vertical resolution (a core extracted from the "Sand" unit in well B-3 shows sand/silt laminations with thicknesses ranging from few mm to 10-15 cm, interbedded with muds).This limitation could be mitigated in future studies by investigating the ability of advanced resistivity logging tools, such as the Rt Scanner service (Schlumberger Ltd.), to resolve the thinly laminated reservoir sand beds and provide more realistic resistivity readings for accurate water and gas saturations determination.
The main results of the wireline logging data processing for the analyzed wells are presented and discussed in the following sections.

Formation Waters
Figure 10a,b show examples of Hingle-type crossplots corresponding to B-2 and B-4 wells, constructed in the pre-interpretation phase and including the "Sand", "Silt", and short post-reservoir Pliocene intervals comprising clean sand beds.Data points are color-coded according to the GR intensity.The input porosity φ = φ D was computed for a matrix density ρ ma = 2.71-2.72 g/cm 3 (mean of core-derived measurements) and a pore fluid density ρ mf ≈ 1.04 g/cm 3 ; the Archie parameters used were derived from SCAL measurements (m = 1.54, n = 1.32).Rt was approximated by the deep investigation Laterolog R A curves (LLD and RLA5).The datasets were filtered via gamma-ray and caliper cut-offs (GR < 70 gAPI, d < 9 in), to retain reasonably clean levels with good hole conditions.Forced least-squares regressions through the matrix point (φ = 0, R t → ∞) and the linear trend of NW located data points identified the water-bearing levels (S w = 1 "water line") and yielded from the maximum slope 1/R w ˆ(1/m): R w = 1.04 Ω m (equivalent salinity: 3932 ppm NaCl) for B-2 well and R w = 1.175Ω m (equivalent salinity: 3475 ppm NaCl) for B-4 well.The data points that define the linear water-bearing trends belong not only to clean sand beds in the post-reservoir sections, but also to few levels within the "Sand" and "Silt" units.This indicates that the main reservoirs and the post-reservoir Pliocene deposits in their proximity likely host similar or identical formation water in the pore space.Consequently, the R w values obtained from the resistivity-porosity analyses are suitable for the quantitative evaluation of the main reservoirs in terms of fluid saturations.11a presents an example of Pickett-type crossplot for B-3 well, including data from the "Sand", "Silt", and post-reservoir intervals (caliper cut-off applied to the data: d < 9 in).The final φ e resulting from log interpretation was used and Rt was approximated by the RLA5 deep Laterolog curve.Figure 11b shows the same crossplot with a GR < 70 gAPI cut-off applied, to retain only the relatively clean levels.A robust least-squares regression through the linear trend of the lowermost data clearly defined the S w = 1 line, confirmed the SCAL-derived cementation exponent from the slope (−1/m) of the best-fit line, and provided R w = 0.95 Ω m (equivalent salinity: 4320 ppm NaCl).The data points defining the linear aquifer trend correspond to clean sand beds from the post-reservoir section and also to some levels from the "Sand" and "Silt" reservoirs, suggesting again that the formation waters have a constant resistivity and salinity throughout the analyzed sequence.In Figure 11a, it can be noticed that the location of the S w = 1 line coincides with the onset of the log(φ) = f (log(R t )) change of trend, from a quasi-vertical data distribution in claystones and mudstones/shales towards a slight R t increase in the gas-bearing reservoir.In the specific petrophysical context of fields "A" and "B", this feature has a methodological importance if such crossplots are used over intervals with very limited porosity and resistivity variation range, when it is difficult to identify a clear linear trend of the waterbearing data.
A comparable result was obtained for B-1 well using the same approach: R w = 1.10 Ω m at formation temperature (equivalent salinity: 3719 ppm NaCl).One may observe that the formation water resistivities and salinities obtained for the wells of field "B" are almost identical, regardless of: (a) the different porosity and resistivity ranges that defined the water-bearing trend (Figure 10) and (b) the input porosity type or the type of resistivityporosity crossplot.
Figure 12a,b show Pickett-type crossplots for A-1 and A-2 wells, including data from the "Sand", "Silt", and post-reservoir intervals (caliper cut-off applied to the data: d < 9 in).The crossplots were constructed using the final φ e from log interpretation, Rt approximated by the RLA5 deep Laterolog R A curve, and Archie parameters obtained from the SCAL measurements.An S w = 1 linear regression with a fixed slope (m = 1.71) fitted through the start of the data distribution change of trend provided R w = 0.6 Ω m for both wells (equivalent salinity: 7000 ppm NaCl).The presented R w results, determined over slightly larger depth intervals than the main reservoir units of the "A" and "B" fields, may be considered as representative for the Pliocene deposits including the reservoirs.These R w values were used for S w evaluation for all the analyzed wells.

Figures 13-15 illustrate examples of deterministic quantitative interpretation results
for the wireline logging data recorded in the final sections of A-1, A-2, and B-2 wells.The main interpretation parameters were presented in Table 4, with ρ ma , m, and n obtained from measurements on core plugs at quasi-reservoir overburden pressure.
The final petrophysical solutions are presented in Figure 13-tracks 7-9, Figure 14tracks 6-8, and Figure 15-tracks 7-9: computed S w and S xo water saturations, bulk volumes of formation fluids (water in the uninvaded zone φ e S w , water in the flushed zone φ e S xo , movable hydrocarbons φ e (S xo − S w ), residual hydrocarbons φ e (1 − S xo ), and the lithological fractions (V clay , V silt , V ma , φ e ).Tracks 9-10 from Figure 14 show the raw NMR data (normalized T 2 transverse relaxation time distribution) and the NMR interpretation results in terms of clay-bound water (CBW), capillary-bound water (BVI), and free fluid (FFI) volumes, as well as the NMR-derived permeability (NMR_K).Two distinct gas-bearing reservoir intervals separated by permeability barriers were identified in the conventional "Sand" and "Silt" units from both fields (R1 and R2 in Figures 13-15), the reservoirs' separation being more evident in field "A".The upper reservoir located in the "Sand" unit is better developed, with a 29 m thickness in well A-1 and 24.5 m in well A-2, and characterized by higher overall porosities (well A-1: maximum φ e = 37%, mean φ e = 21.6%;well A-2: maximum φ e = 32%, mean φ e = 15.7%).The thickness of the secondary reservoir ranges from 17 m in well A-1 to 6 m in well A-2, the porosities being lower (well A-1: maximum φ e = 26%, mean φ e = 14.0%; well A-2: maximum φ e = 24%, mean φ e = 14.8%).In well B-2 (field "B"), the gas-bearing reservoirs delineated within the "Sand" and "Silt" conventional units are comparable, with a 12.5-13 m thickness, mean φ e of 14.4-15.2%,and maximum φ e of 27-36%.
The minimum water saturations (maximum gas saturations) computed in the reservoir intervals from the "Sand" units were: A-1 well-S w,min = 8.2% (S h,max = 91.8%),A-2 well-S w,min = 4.4% (S h,max = 95.6%),B-1 well-S w,min = 34.8%(S h,max = 65.2%),B-2 well-S w,min = 33.8%(S h,max = 66.2%),B-3 well-S w,min = 34.9%(S h,max = 65.1%),B-4 well-S w,min = 9.1% (S h,max = 90.9%).Rapid alternations of resistive (gas-bearing sands) and conductive (claystones or mudstones/shales) thin beds or laminations, as the ones observed in cores extracted from the field "B" wells, may not be correctly resolved by the Laterolog tools used.Due to this resistivity suppression (the maximum R A readings in the gasbearing intervals of field "B" field are one order of magnitude lower than those recorded in the wells from field "A"), the S w,min values computed for B-1, B-2, and B-3 wells may be overestimated and the S h saturations underestimated.
The NMR log interpretation results revealed additional sources of resistivity suppression and possible S h underestimation even in clean, thick reservoirs.In the wells A-2, B-2, B-3, and B-4, the deep R A maximum readings and the corresponding S w,min and S h,max computed saturations do not occur at the tops of the gas-bearing intervals but below them (e.g., the top of R1 reservoir-tracks 5 and 6 in Figure 14).This is caused by a large amount of irreducible capillary-bound water trapped in small pores and hosted at the top of the reservoir intervals (note the 13% capillary-bound water volume at 1125 m TVDSS, in Figure 14-track 10).The presence of silt layers on top of underlying sands, possibly due to underwater sediment gravity flows followed by graded bedding, could explain such abnormal low-resistivity zones, which, however, may produce water-free gas.
The validity of the quantitative log interpretations was evaluated by: (a) comparison between log-derived petrophysical parameters and the same parameters resulted from RCAL and SCAL core measurements at quasi-reservoir overburden pressure; (b) reconstruction of the theoretical response of water-saturated formations by means of a "wet resistivity" curve R 0 = F R w = a φ e −m R w and its comparison with the recorded deep investigation Laterolog R A curves.An example of the comparison and good agreement between log-derived S w and φ e and the core-derived equivalents (CORE_SW-water saturation measured on core plugs at 1200 psig pressure, CORE_PHI-porosity measured on core plugs at 1400 psig pressure) is shown in Figure 13-tracks 11 and 12.The resistivity comparison control criterion is illustrated in Figure 13-track 6, Figure 14-track 5, and Figure 15-track 6.The close match between the measured deep resistivities (RLA5, RLLD) and the theoretical R 0 curves, except in gas-bearing intervals, indicates that the log-derived (R w , ρ clay , φ Nclay ) and core-derived (ρ ma , m, n) interpretation parameters were adequate, and the saturation model used was realistic-Equations ( 14) and (15).
At the scale of the analyzed reservoirs, the MLR technique was found to be effective in predicting continuous permeability curves, with control provided by any available set of permeability measurements.Track 13 from Figure 13 illustrates the fit of a function K = f (φ e , V clay , S w ) simultaneously for A-1 and A-2 wells, by jointly using multiple control datasets: the Klinkenberg-corrected permeabilities (CORE_K) measured on core plugs from well A-1 at 1400 psig overburden pressure and the continuous K NMR curve from A-2 well (NMR_K in Figure 14-track 10).The modeled permeability curve MLR_K [mD] = 10ˆ(8.730999− 11.858403 φ e − 6.703942 V clay − 5.804627 S w ) is shown in track 12. Another example is presented in Figure 15-track 10, where a synthetic permeability curve MLR_K [mD] = 10ˆ(0.44883338+ 9.57922687 φ e + 0.18209513 V clay − 3.25371721 S w ) provided a close fit to the set of Klinkenberg-corrected permeabilities (CORE_K) measured on core plugs from well B-2 at quasi-reservoir pressure.
It is noticeable that slight background gas saturations resulted from the log interpretation in the pre-reservoir (Late Miocene/Pontian, generally argillaceous deposits) sections of the analyzed wells (e.g., Figures 14 and 15), confirming gas shows reported while drilling.Taking into account the widespread methane seepage phenomena in the younger sediments of the Black Sea offshore area, we consider these minor gas saturations as petrophysical evidences of an active/ongoing gas migration process from the deeper levels of the Miocene pelitic deposits.

Fluid Contacts
Figures 16 and 17 show the wireline formation pressure datasets from fields "A" and "B" (three outlying pressure readings from well B-1 at TVDSS 1784.4,1787.4,and 2461.9 m were removed).For both fields, the depth coverage with pressure readings (including pre-reservoir Pontian deposits and post-reservoir Dacian-Romanian deposits) allowed an adequate definition of the formation water trend.The contrasting gas trend was better outlined in field "A" than in field "B", due to the larger amount of data and higher gas saturations.At a field-level scale, the intersection of the hydrostatic pressure trends corresponding to the two reservoir fluids suggest GWC depths of 1157.5 m TVDSS for field "A" and 1123.1 m TVDSS for field "B" (for relatively permeable and gas-bearing reservoirs, the GWC depths normally correspond to the free water levels-FWLs).The pressuredepth gradients computed for field "A" were (∆p/∆z) 1 = 0.396 psia/ft (0.089 bar/m) and (∆p/∆z) 2 = 0.071 psia/ft (0.016 bar/m), corresponding to in situ fluid densities ρ f1 = 0.913 g/cm 3 (water) and ρ f2 = 0.163 g/cm 3 (gas).For field "B", the computed gradients were (∆p/∆z) 1 = 0.439 psia/ft (0.099 bar/m) and (∆p/∆z) 2 = 0.116 psia/ft (0.026 bar/m), corresponding to fluid densities ρ f1 = 1.013 g/cm 3 (water) and ρ f2 = 0.267 g/cm 3 (gas).Except the water density obtained for field "B" (aquifer trend defined by p [psia] = 1.441TVDSS [m] + 24.267), the other inferred fluid densities are either too low or too high with respect to the likely ρ w (≥ 1 g/cm 3 ) and to the ρ h predicted by the ideal gas law (≈ 0.08 g/cm 3 ).The inaccuracy of reservoir fluid densities estimation implies the uncertainty of GWC depths estimation from wireline pressure data.
The results of a detailed pressure data analysis performed for each well of the "A" and "B" fields are synthesized in Table 5.
For wells B-1 and B-3, the insufficient data coverage in the "Sand" and "Silt" reservoir units allowed the delineation of a single pressure trend related to formation water.For the wells of field "A", the estimated density of "fluid 1" (water) was lower than the anticipated ρ w and the estimated density of "fluid 2" (gas) was significantly higher than the likely gas density ρ h .For the wells of field "B", the estimated ρ w was veridical but ρ h was lower or higher than expected.This uncertainty can be partially explained by the small number of pressure measurements in the reservoir intervals of field "B".However, for field "A", the deviations of estimated ρ w and ρ h from the expected values cannot be explained by insufficient pressure data.The quantitative well log interpretation proved extremely useful for evaluating the validity of GWC estimations based on pressure data.The formation pressure measurements and the pressure-depth trends for the wells A-1, A-2, and B-2 are displayed in Figure 13-track 10, Figure 14-track 11, and Figure 15-track 11 (the horizontal dashed line marks the pressure trends intersection, i.e., the GWC estimated from pressure data).The petrophysical analysis results indicate that distinct fluid contacts should be considered within each of the reservoirs, instead of a single "averaging" GWC obtained exclusively from a pressure survey.The log-based fluid contact depths (GWC1, GWC2) were delineated at the base of the gas column in each reservoir, where a significant decrease in the computed S w values occurs.As observed, there are significant depth differences between the two types of estimations, ranging from 4.1 to 15.4 m.Besides the questionable fluid densities, an additional indication regarding the uncertainty of pressure-based fluid contact depths was provided by the NMR log available for well A-2 (Figure 14-tracks 9, 10).The GWC estimated from the formation water and gas pressure trends at 1158.4 m TVDSS is located in an impermeable mudstone/shale interval that separates the two gasbearing reservoirs.At that depth, the NMR results show no free fluids but only bound (immobile) water.
Especially for field "A", the ρ w underestimation and ρ h overestimation may be explained by considering that the wireline testers did not read pressures corresponding to true formation water and gas, but to mixtures of water (mud filtrate) and gas or gas and water, respectively.Additionally, the small number of pressure readings outside the "Sand" and "Silt" units, in the pre-reservoir and post-reservoir intervals of field "A", did not allow a better definition of the water trend.The more realistic ρ w values estimated for field "B" and the GWC depth obtained for well B-2 (1124.4m TVDSS), closer to the log-based fluid contacts, could be attributed to actually lower gas saturations in the "Sand" and "Silt" units of that field.

Conclusions
This study aimed to present and discuss the main issues related to the petrophysical evaluation of two biogenic dry gas fields of Early Pliocene (Dacian) age from the Romanian continental shelf-Western Black Sea basin (conventionally denoted "A" and "B").The sands and silts reservoirs have a better quality in field "A", developed in a shallow marine environment, than in the marginal marine (deltaic) field "B", which shows a higher silt content, fewer thick sand intervals, and frequent thinly laminated sand-shale sequences.A vertical, coarsening upward, variability of the reservoirs generally exists in both fields, with silty deposits in the lower part of the reservoir intervals overlain by sands.
The quantitative interpretation of the wireline well logs recorded in six exploration wells was performed using a deterministic workflow.Some of the interpretation parameters (ρ ma , m, and n) were derived from routine and special analyses carried out on core plugs at quasi-reservoir confining pressure, whereas the other parameters (ρ clay , φ Nclay , R clay , R w ) were estimated statistically or via crossplot techniques from the logs.The formation pressure data measured in all the wells were processed and interpreted in terms of probable gas-water contact (GWC) depths, which were evaluated by comparison with the results of well log interpretation.
A significant problem that affected the well log interpretation in the analyzed fields was the uncertainty related to R w values and, consequently, to S w and S h evaluation.Without representative water samples collected (uncontaminated by drilling mud filtrate) no accurate and direct R w information was available.On the other hand, no obvious water-bearing beds were discernible within the Dacian reservoir intervals.The approach used in this study for obtaining realistic R w values was the analysis of resistivity-porosity dependencies over slightly extended depth intervals, which included the Dacian reservoirs and short sections of post-reservoir Pliocene deposits where clean and probably waterbearing sand beds were present.
The main conclusions emerging from this study are: 1.
The integration of core measurements in the well log interpretation methodology had a major impact on the validity of the obtained results.Core-derived petrophysical measurements were used both as input computation parameters and also to check and validate the main reservoir parameters resulted from the interpretation (φ, S w , K).
Additionally, the core-derived ρ ma and m provided the necessary constraints for the realistic estimation of R w from resistivity-porosity dependencies; 2.
The approach used for R w estimation, i.e., the use of resistivity-porosity dependencies and the extension of the analysis interval to segments of Pliocene deposits above the gas-bearing reservoirs, proved to be effective.The capability of R t -φ dependencies to reveal linear data trends in clean water-bearing formations with constant R w showed that parts of the Dacian reservoirs (the very limited bottom water zones underneath the gas columns) and the adjacent post-reservoir Pliocene sections host similar formation waters.This allowed the determination of realistic R w values, which were used for S w evaluation in the analyzed wells; 3.
Particularly in the "B" field, the Laterolog (DLL and HALS) resistivity curves are likely suppressed to varying degrees in each well, leading to possible S w overestimation and S h underestimation and negatively impacting gas reserve evaluation.One cause of this problem may be represented by the alternance of thin (millimeter to decimeter thick) resistive and conductive layers of sand and mudstone/shale, which are averaged by the resistivity tools due to their limited vertical resolution; 4.
An additional source of resistivity logs suppression, both in "A" and in "B" fields, is represented by a high content of capillary-bound water, probably trapped in the small pores of silt intervals.The NMR logging performed in well A-2 from field "A" was essential for understanding the cause of these low-resistivity intervals, sometimes located at the top of resistive gas-bearing reservoirs sands; 5.
The estimation of GWC depths using formation pressure surveys (frequently considered the main and preferred source of data for defining the fluid contacts) should always be checked and validated using the well log interpretation results.There is a significant degree of uncertainty in using the hydrostatic pressure trends identified in the analyzed wells to estimate the fluid contact position, due to the possibility that the pressures read by the wireline testers are not representative of formation water and gas, but might reflect mixtures of water (or mud filtrate) and gas in various ratios.The NMR log recorded in well A-2 provided valuable insight into the intervals with free fluids and with bound (immobile) water and allowed an assessment of the pressure-derived GWC validity; 6.
The well log interpretation results indicate that the Dacian reservoirs from fields "A" and "B" cannot be treated as single units, because they include two distinct reservoir intervals separated by permeability barriers of various thicknesses.Consequently, separate fluid contacts should be considered for each reservoir interval, instead of a single GWC obtained from the pressure gradients analysis.

Figure 4 .
Figure 4. Example of wireline geophysical logs recorded in an exploration well (B-4) from the Early Pliocene gas field "B" (adapted from [15]).

For
the older wells, the openhole logging program in the final sections was carried out with Atlas Wireline Services equipment and included DLL-Dual Laterolog, MLL-MicroLaterolog, GR-Gamma Ray, ZDL-Compensated Z-Densilog, CN-Compensated Neutron, MAC-Multipole Array Acoustilog, and FMT-Formation Multi-Tester tools.

Figure 5 . 6 Figure 6 .
Figure 5. Neutron-density crossplot corresponding to the Pliocene ("Sand", "Silt", and post-reservoir) deposits of field "A" and including data from A-1 and A-2 wells (7180 data levels).A caliper cut-off (d < 9 in) was used to eliminate the bad hole levels.The approximate location of a clay point (φ Nclay , ρ clay ) representative for the entire Pliocene sequence is indicated.SS, LS, and DOL are theoretical ρ b = f (φ N ) response curves for the main reservoir lithologies (sandstone, limestone, dolomite) [32].
), and φ Dclay [V/V] is the wet clay density porosity.The neutron-density combination is one of the best clay indicators in good hole conditions and water-bearing intervals but is affected by washouts or rugose zones and computes negative V clayND values in gas-bearing intervals where distinctive crossovers (φ N < φ D ) occur.The clay indicators comparison allowed the tuning of the GR clean and GR clay endpoints, using as a reference (V clay = 1) the maximum and repeatable separations φ Nclay − φ Dclay , to obtain V clayGR ≈ V clayND over the largest part of the intervals.For the analyzed wells, the GR clean range was 45-60 gAPI, the GR clay range was 100-114 gAPI, and the reference φ Nclay − φ Dclay maximum separation range was 0.23-0.29.Figure7shows an example of the V clay indicators comparison for the A-2 well, over an interval that includes the post-reservoir, reservoir, and pre-reservoir sections.No nonlinearity was observed on V clayGR = f (V clayND ) crossplots for any of the wells, thus validating the adopted linear V clayGR estimator.

Figure 7 .
Figure 7. Clay volume derived from GR log (track 4) and from the neutron-density logs combination (track 6).V clayGR (track 5) is referenced against the φ Nclay − φ Dclay separation (tracks 6 and 7) to provide V clayGR ≈ V clayND in the shaly water-bearing intervals (track 8) and is used as the final clay indicator (track 9).TNPH_SS is the matrix-corrected (sandstone) neutron log.

4 7 8 8 .
Figure Estimation of fluid contact depths by means of formation pressure measurements, even if the actual contact is not intercepted by exploration wells.

Figure 9 .
Figure 9.(a) Core permeability-core porosity relationship for the pre-reservoir, "Sand", "Silt", and post-reservoir deposits of "A" and "B" fields; (b) NMR permeability-gamma ray intensity dependence for the same sequence of deposits in the A-2 well (567 data levels).The color code corresponds to S w values resulting from the petrophysical interpretation of well logs.

Figure
Figure 11apresents an example of Pickett-type crossplot for B-3 well, including data from the "Sand", "Silt", and post-reservoir intervals (caliper cut-off applied to the data: d < 9 in).The final φ e resulting from log interpretation was used and Rt was approximated by the RLA5 deep Laterolog curve.Figure11bshows the same crossplot with a GR < 70 gAPI cut-off applied, to retain only the relatively clean levels.A robust least-squares regression through the linear trend of the lowermost data clearly defined the S w = 1 line, confirmed the SCAL-derived cementation exponent from the slope (−1/m) of the best-fit line, and provided R w = 0.95 Ω m (equivalent salinity: 4320 ppm NaCl).The data points defining the linear aquifer trend correspond to clean sand beds from the post-reservoir section and also to some levels from the "Sand" and "Silt" reservoirs, suggesting again that the formation waters have a constant resistivity and salinity throughout the analyzed sequence.

Figure 16 . 17 Figure 17 .
Figure16.Estimation of gas-water contact (GWC) for field "A" by jointly using pressure measurements from available wells.

Table 1 .
Conventional coring program and wireline formation pressure measurements conducted in the analyzed gas exploration wells.

Table 2 .
Summary statistics of routine core measurements performed at quasi-reservoir confining pressure for A-1, B-2, and B-3 wells-"Sand" and "Silt" intervals.

Table 5 .
Results of wireline formation pressure data processing and interpretation: pressure-depth gradients, fluid densities, and GWC depths estimated from the gradients' slope and intersection.