Bacterial methane is a dominant hydrocarbon component in the Northern Croatia’s Adriatic offshore proven hydrocarbon reservoirs. As onshore reserves are declining the potential of shallow gas accumulation, analogue to the Adriatic ones, are being tested. A part of the Lonja Formation (Pliocene Pleistocene and Holocene), the Ravneš Member (Early Pliocene age) is analysed for its maturation and expulsion regarding bacterial methane as potential source and reservoir rocks, especially as it is thermally immature. Two approaches were used for the initial lithology modelling processes - the convergent interpolation and sequential Gaussian simulation. Possibility for trapping and forming of accumulations was also modelled. Results show that selected member has a large source rock potential for bacterial methane with a total volume of 411.53 × 109 m3 for the Convergent interpolation model and 520.56 × 109 m3 for the sequential Gaussian simulation model of expulsed bacterial methane
The Miocene sediments within the area of research are relatively well explored and with proven hydrocarbon systems of Miocene age. There are a total of 17 oil and gas fields, majority of which are still in the production but with declining reserves . Largest amounts hydrocarbon reserves belong to the oil discovered in the most active exploration period within the Sava Depression (1950 - 1990). The natural gas was considered as a far less worthy resource when compared to oil, but today, such accumulations are explored and considered as a valuable resource . Analogy for the research of Pliocene and Pleistocene sediments, regarding hydrocarbon potential, can be also found in Adriatic gas accumulations. They are of biogene origin, i.e. result of bacterial degradation of organic (dominantly terrestrial origin) matter in a thermally immature stage of hydrocarbon maturation . The type of organic matter in the Croatian part of the Pannonian Basin System Pliocene sediments is similar in composition to the Adriatic ones regarding kerogene (type III and IV) and maturity [4–7]. Lithological composition of Ravneš member and its age can also be correlated to the Northern and Middle Adriatic Pliocene-Pleistocene where poorly consolidated sandstones, sands, marls and clays of thicknesses up to 2000m can be observed [5, 7, 8].
Gas accumulations in the Adriatic hold large reserves . For example, the largest Ivana Gas Field started production with 8 × 109 m3 of recoverable reserves of natural gas [9, 10] and the probability of finding a similar amount of these volumes in the Pannonian Basin System Sava Depression in resembling sediments is a good enough reason for detailed explorations. A preliminary assessment of the generation of bacterial methane from the Ravneš Member was already published  but it regarded only the generative potential with no migration and accumulation.
2 Geographical and geological settings
Investigation area is located in the Croatian part of the Pannonian Basin System(CPBS), more accurately Sava Depression. It covers roughly around 340 square kilometres and is located between Ivanic Grad on the West, Novska on the East, Moslavacka gora Mt. on the Northeast and Petrinja on the South (Figure 1).
Investigated stratigraphic section belongs to the youngest formation of the Neogene-Quaternary infill of the CPBS - namely Ravneš Member within Lonja Formation. The whole Neogene sequence can be up to 3500mthick in the explored area . It is composed of conglomerate and breccia (pre-Badenian and Lower Badenian), limestone and calcarenite (Badenian, Sarmatian), mostly sandstone with marls (Pannonian and Pontian) and poorly consolidated sandstones, sands and clays (Pliocene and Quaternary) as described in [6, 12–14]. Ravneš Member is part of the Lonja Formation along with five other members (Figure 2). It was selected for this preliminary analysis of hydrocarbon potential because of its lithological composition as it is in some parts almost exclusively composed of coal and in other parts holds large percentage of coal in its total thickness.
For this kind of Basin analysis in a preliminary research, there are several steps that need to be made or to be obtained from previous research. Those were: (1) structural modelling, (2) lithology determination (3) organic geochemistry and kinetics definition, (4) calibration. Lithology definition was a result of lithology analysis within 31 well. Part of the process was made using the Schlum-berger Petrel software (structural modelling, lithofacies mapping) while the basin modelling was carried out using the Schlumberger PetroMod software.
3.1 Structural modelling
A subsurface model of Pliocene, Pleistocene and Holocene sediments divided by five e-log markers and borders and into six members was built based on . These e-log markers and borders were defined by standard deviation curves of well log data (spontaneous potential and resistivity curves) and represent a general lithological changes controlled by tectonics and environment influence . The targeted Ravneš Member is located between e-log markers K’ and I’ or border J’ (Figure 2), depending on the location within the depression. Members’ thickness varies from 0 m on the marginal area of the depression to 167 m in its deepest parts. Initial modelling was made in Petrel  but for the basin modelling surfaces of each well log markers and the topography were transferred to PetroMod software in which the cell based model was re-built (Figure 3).
3.2 Lithology determination
Initial lithology analysis defined the Ravneš Member as impermeable, shale like dominant facies with substantial amounts of coal (lignite) and organic rich clay . This kind of approach regarded the two types of lithology with transitional mix-type lithology in between. The lithology was controlled  by the percentage of coal within the Ravneš Member for wells within the researched area. Also, the Ravneš member was as single vertical cell of variable thickness so the average cell height within the model was from 0 to 167 m.
What was initially disregarded was the fact that there are areas within the member which have substantial amount of sand, especially when the member was layered in three sublayers in equal proportions - Ravneš I - III. Layering was performed with a maximum cell height of 56 m, obtaining more detailed results than in . The lithology of a cell was controlled by either the ratio (>0.25) of permeable and impermeable component (e.g., ss/sh as used in ), or by the coal content percentage (%Co) if the lithology was dominantly impermeable (Table 1). If ss/sh was greater than 0.25 than the lithology was defined exclusively by that ss/sh value regardless of the %Co value.
Two mapping techniques were used. First was a relatively classic mapping approach using the Convergent interpolation (CI) algorithm . It obtains relatively smooth transition and is a generally good algorithm for mapping structural surfaces in the subsurface . The maps show relatively smooth transitions (Figure 4a, Figure 5a) depicting a less realistic image when compared to real transition of facies, especially regarding period of the Ravneš Member deposition . Second approach was appliance of the sequential Gaussian simulations (SGS) [18, 19] by which a more realistic map was made with abrupt facies transitions (Figure 4b, Figure 5b)which are characteristic for the Pliocene environments in the Sava Depression . The results were two sets of maps, each containing three maps - one per layer. All are used as facies maps in the basing modelling.Atotal of ten lithology classes were defined (Table 1).
3.3 Organic geochemistry properties and kinetics definition
Organic geochemistry data for coals and organic rich clays of the Ravneš Member were obtained from [4, 6]. The coals and clays show dominantly type III gas prone kerogene in a thermally immature stage. At the moment, biogenic reactions should be at the maximum stage regarding expulsion of bacterial methane in the deepest parts of the depression (Figure 7a, b) defined by the transformation ratio (TR) gained from maturity modelling which correspond to the expulsion of bacterial methane in the initial stages of burial . The TR value is presented only for the source rock facies defined in Table 3. Organic geochemistry properties (kerogene type, maturity) of the source rocks are similar to the ones in the Northern Adriatic [5, 7].
Moreover, properties of selected coals from the Ravneš Member are shown in Table 2 along with the several samples obtained from the member in the top and base of the Ravneš Member together with organic rich clays. These properties were then assigned to the Ravneš Member facies defined in PetroMod (Table 3). Algorithm for transformation of the organic matter to hydrocarbons (kinetics) was set for biogenic reactions.
|TOC (%)||HI (mg HC/g Corg)||OI (mg CO2/g TOC)|
|Lithology||PetroMod lithology||TOC||Kinetics||HI value||PSE|
|Ravnes_clay||Shale (typical)||1||Biogenic reaction||70||Source Rock|
|Ravnes_Sh80Co20||Sh80Co20||5||Biogenic reaction||80||Source Rock|
|Ravnes_Sh60Co40||Sh60Co40||10||Biogenic reaction||85||Source Rock|
|Ravnes_Sh40Co60||Sh40Co60||15||Biogenic reaction||90||Source Rock|
|Ravnes_Sh20Co80||Sh20Co80||20||Biogenic reaction||100||Source Rock|
|Ravnes_coal||Coal||35||Biogenic reaction||120||Source Rock|
3.4 Organic geochemistry properties and kinetics definition
The model output is generated based on all of the data provided in the input. Input data for the model also include “boundary conditions” and paleo geometry. Boundary conditions define the basic energetic conditions for temperature and burial history of the source rock and, consequently, for the maturation of organic matter through time1. These are surface water interface (SWIT), heat flow (HF) and the paleo water depth (PWD). The data can be altered during each modelling step so the modelled temperature and vitrinite reflectance (VR) correspond to the actual ones in the wells. Heat flow was set at a 70 mW/m2 as shown in . SWIT was determined from regional values from  and PWD was set according to the reconstructions in  in the range from 0 m (recent) to 50 m regarding environments that existed for deposition of members of the Lonja Formation.
Calibration was made only with the fitting of the temperature values (Figure 8), because VR is a poor indicator when dealing with immature source rocks. The model was made for temperature, generation, migration and accumulation for both CI and SGS models.
An estimate of possible generated bacterial methane was obtained for the CI and the SGS maps. Both showed substantial amount of bacterial methane generated - 411.53× 109 m3 for the CI model and 520.56 × 109 m3 for the SGS. The majority of the expulsed volume occurred between 3.5 and 2.5 Ma (Figure 9) after which only small ammounts were expulsed. This was the result of small rate of burial during Quaternary as sedimentation of more than a 100 m thick sediments rarely occurred outside the depocenters, which were at peak expulsion before 2.5 Ma.
The TR values are at about 50% for the amount that can be generated in the biogenic reactions domain and the remaining potential is related to the source rocks in the shallow areas (Figure 7). In contrast to large generated volumes of biomethane, the accumulated volumes are rather poor - 0.20 × 109 m3 for the SGS and 0.06 × 109 m3 for CI model. Small accumulation volumes wereexpected within the Ravneš member as it was regarded as dominantly a source rock facies with sporadic occurrences of reservoir facies (Ravnes_St20Ss80 and Ravnes_sand lithology defined in Table 1). There are a large number of accumulations within the SGS model with very small amounts and only a few with volumes of over 1 × 106 m3 (Figure 10). Most of the accumulations, especially of the smaller volumes are a result of stratigraphic trapping due to the lateral lithology change (Figure 10, enlarged area). The volume of these accumulations can be underestimated as the cell grid of the model was relatively coarse (100 × 100 m) in regard to the area of the accumulation (largest one 0.50 km2) according to  There are only a few accumulations in the CI model which are a result of a gradual lithology change in the up dip direction and a small number of possible structural traps.
A substantial volume of generated and expulsed bacterial methane is determined by the basin modelling process from the Ravneš Member. Both the CI and SGS lithology approaches gave similar results when volumes of generated methane was regarded - 411.53 × 109 m3 for the CI model and 520.56 × 109 m3 for the SGS. On contrary, volumes of the possible trapped bacterial methane in Ravneš Member are very small and significantly differ - 0.20 × 109 m3 for the SGS and 0.06 × 109 m3 for CI model. The differences in generated and trapped volumes resulted from different lithology distribution between the CI and the SGS derived lithofacies maps.
Although the convergent interpolation provides a safer result (probabilistic approach), the SGS (stochastic approach) is more corresponding to the dynamics and size of the depositional environment during the Lower Pliocene in the Sava Depression. Presented analysis supports the conclusions  that Ravneš Member is a good source rock interval for bacterial methane generation, but commercial size accumulations within the member itself are not to be expected except if regarded as secondary targets when drilling for deeper accumulations. Technical problems with inflow of sand into well which can arise from producing in poorly consolidated sands and sandstones, similar to the ones in North Adriatic, should also be taken into consideration when evaluating the possibility of economical accumulations within the Ravneš Member as well as all other accumulations in Lonja Formation.
Author would like to thank the Schlumberger Company for providing the Academic licenses of Petrel and PetroMod Software to the Faculty without which this kind of analysis, in its current extent, would not be possible. Authors would also thank the University support through funding of the research “Geomath-ematical research and mapping of selected Croatian depositional environments from Holocene to Lower Miocene” in 2015.
 Velić J., Malvić T., Cvetković M., Vrbanac B., Reservoir geology, hydrocarbon reserves and production in the Croatian part of the Pannonian Basin System. Geol. Croat., 2012, 65(1), 91-101.Search in Google Scholar
 Križan J., Proizvodnja prirodnog plina iz hrvatskog dijela pod-morja Jadrana [Exploitation and production of natural gas from the Croatian part of the Adriatic offshore - in Croatian]. Nafta i plin, 2014, 34(140), 24-33.Search in Google Scholar
 Whiticar M.J., Correlation of Natural Gases with Their Sources. In: Magoon L.B., Dow W.G. (Eds.), The Petroleum System-From Source to Trap, AAPG Memoir, 2009, 60, 261-284.Search in Google Scholar
 Barić G., Coal-bed hydrocarbon characteristics of the Vrbak - Benkovica exploration locations. Nafta, 1993, 44(10), 543-554.Search in Google Scholar
 Vulama I., Pliocene Source Rocks, Miocene Reservoir Rocks and Origin of the Gas Accumulation of the Irma Field (Northern Adriatic, Croatia) Based on Wireline-Logging. Geol. Croat., 1997, 50(1), 89-103.Search in Google Scholar
 Cvetković M., Lithostratigraphic units of the third neogene-quarternary megacycle in the Sava Depression and their petroleum potential. PhD Thesis, University of Zagreb, Faculty of Mining, Geology and PetroleumEngineering, Zagreb, 2013 (in Croatian).Search in Google Scholar
 Marić Đureković Ž., Lithofacies and stratigraphy of Pleistocene deposits in North Adriatic Offshore by using high-resolution well logs. PhD Thesis, University of Zagreb, Faculty of Mining, Geology and Petroleum Engineering, Zagreb, 2011.Search in Google Scholar
 Velić J., Malvić T., Depositional conditions during Pliocene and Pleistocene in Northern Adriatic and possible lithostratigraphic division of these rocks. Nafta, 62(1-2), 25-32.Search in Google Scholar
 Zelic M., Mlinaric Ž., Jelic Balta J., Croatian Northern Adriatic Ivana gas field ready for development (Reservoir characteristics and gas inflow conditions into the well). Nafta, 1999, 50(1), 1937.Search in Google Scholar
 Velić J., Malvić T., Cvetković M., Velić I., Stratigraphy and petroleum geology of the Croatian part of the Adriatic Basin. J. Pet. Geol., 2015, 38(3), 281-300.Search in Google Scholar
 Cvetković M., Velić J., Biogenic reactions and methane expulsion modelling from source rocks of Ravneš Member, Sava Depression. In: Horvath J., Cvetković M., Gabor Hatvani I. (Eds.), The Geomathematical Models: The Mirrors of geological reality or science fictions, Szeged, Hungary. Hungarian Geology Society - University of Szeged, Geology and Paleontology, 2015, 5057.Search in Google Scholar
 Saftic B., Velić J., Szanto O., Gyorgyi J., Ivkovic Ž., Tertiary subsurface facies, source rocks and hydrocarbon reservoirs in the SW part of the Pannonian Basin (northern Croatia and southwestern Hungary). Geol. Croat., 2003, 56(1), 101-122.Search in Google Scholar
 Malvic T., Sucic A., Cvetković M., Resanovic F., Velić J., Low permeability Neogene lithofacies in Northern Croatia as potential unconventional hydrocarbon reservoirs. Cent. Eur. J. Geosci., 6, 2, 2014, 182-194.Search in Google Scholar
 Troskot-Corbic T., Velić J., Malvić T., Comparison between the Middle Miocene and the Upper Miocene source rock formations in the Sava Depression (Pannonian Basin, Croatia). Geol. Croat., 2009, 62(2), 123-133.Search in Google Scholar
 Cvetković M., Malvić T., Defining electro-log markers in poorly consolidated, heterogeneous clastic sediments using standard deviation data trends - an example from the Sava Depression, Pannonian Basin System. In: Horváth J., Wágenhoffer A., Geiger J., Cvetković M., Malvić T. (Eds.), 16th Hungarian & 5th Croatian and Hungarian congress abstract and program book, Szeged, Hungary. Hungarian Geology Society, INA Oil Industry Plc., 2013, 1-4.Search in Google Scholar
 Forgotson J.M., Review and classification of quantitative mapping techniques. AAPG Bull., 1960, 44(1), 83-100.Search in Google Scholar
 Mastroianni G., Milanovic G., Interpolation Processes. Springer-Verlag, Berlin, 2008.Search in Google Scholar
 Novak Zelenika K., Malvić T., Stochastic simulations of dependent geological variables in sandstone reservoirs of Neogene age: A case study of Kloštar Field, Sava Depression. Geol. Croat., 2011, 64(2), 173-183.Search in Google Scholar
 Novak Zelenika K., Cvetković M., Malvić T., Velić J., Sremac J., Sequential Indicator Simulations maps of porosity, depth and thickness of Miocene clastic sediments in the Kloštar Field, Northern Croatia. J. Maps, 2013, 9(4), 550-557.Search in Google Scholar
 Lenkey L., Dövényi P., Horváth F., Cloetingh S.A.P.L., Geothermics of the Pannonian basin and its bearing on the neotectonics. EGU Stephan Mueller Special Publication Series, 2002, 3, 29-40.Search in Google Scholar
 Wygrala B., Integrated Study of an Oil field in the Southern Po Basin, Northern Italy, Berichte der Kernforschungsanlage Jülich 2313, 1989.Search in Google Scholar
 Slavinic P., Cvetković M., Volume calculation of subsurface structures and traps in hydrocarbon exploration - a comparison between numerical integration and cell based models. Open Geosci., 2016, in press, 10.1515/geo-2016-0003.Search in Google Scholar
©M. Cvetković, published by De Gruyter Open.
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.