2017 | 70/2 | 87–92 | 10 Figs. | 3 Tabs. | www.geologia-croatica.hr Journal of the Croatian Geological Survey and the Croatian Geological Society 1. INTRODUCTION AND GEOLOGICAL SETTING The investigated well is located in the southern part of the Da­ nube-Tisza Interflow, Hungary. Samples for analysis were taken from the lower Pannonian age basal marl succession, namely from the Endrőd Marl Formation (JuHáSz, 1992; 1998). This basal marl sequence is an initial member of the Pannonian-Plio­ cene age prograding – dominantly siliciclastic – lacustrine self- system (Magyar et al., 2013). In the early stage of deposition due to the rapid thermal subsidence of the Pannonian basin, a so called „starved” sedimentation phase formed in the basin. Dur­ ing this stage, the siliciclastic sediment input was limited to the edges of the depocentre and a carbonate rich sequence was de­ posited almost throughout the Pannonian basin area. This se­ quence generally began with a high carbonate content calcareous marl, to marl sequence (Tótkomlós Calcareous Marl Member) overlain by a deep water, pelagic clay-marl sequence (Nagykörű Claymarl Member), developed with a continuous transition from the underlying carbonate rich rocks (JuHáSz, 1992; 1998). The carbonate content decreases upwards, caused by the increasing terrigenous material input (due to the prograding shelf-system), providing a characteristic and easy to recognize wireline pattern (Figure 1). In the study area, the calcareous marl-marl sequence deve- loped almost everywhere except for the south-western, shallow- most part of the basin (Figure 2). The calcareous marl has several lithofacies, but in general dark grey, black, laminated-microlaminated facies were depo- sited in the deep basin areas, while under the shallow water con­ ditions light grey – yellowish or brownish grey, whitish coloured, massive or laminated facies accumulated (JuHáSz, 1992; 1998; LEMbErkovICS, 2016). The investigated samples represent Pore characterization, acoustic and permeability measurements on core plug triplets, from the Miocene Tótkomlós (Calcareous Marl) Formation of the Pannonian-basin, Hungary Péter Koroncz1, Péter Ács1, Viktor Lemberkovics2 and Ferenc Fedor1 1 GEOCHEM Ltd., Viola str. 55/1., H-7761 Kozármisleny, Hungary; (geochem@geochem-ltd.eu) 2 RAG Hungary Ltd., Bocskai rd. 134-146., H-1113 Budapest, Hungary; (Viktor.Lemberkovics@rag-hungary.hu) doi: 10.4154/gc.2017.09 Abstract The Tótkomlós Calcareous Marl (TCM), a Member of the Late Miocene Endrőd Formation was investigated as a potential cap and source rock and also as a local tight reservoir in the Panno- nian Basin. Only a limited dataset is available for petrophysical characterisation of this formation. The study reports on a complex measurement campaign performed on three core triplet samples of the Tótkomlós Calcareous Marl, including pore structure and petrophyiscal analysis at various pressure conditions. Direct laboratory measurements of compressional (P) and shear (S) wave velocities on oriented rock samples provide information about the anisotropic behaviour of the studied samples. A quantitative description of seismic anisotropy can improve the quality of seis- mic data processing. Pore structure investigations indicated plate-like materials with a typical pore throat size between 75-110 nm. Gas permeability measurements showed very low permeability values in the order of magnitude between 10-16-10-18 m2. Weak acoustic anisotropy is observed both for P- and S- waves. Decreasing anisotropy with increasing confining pressure indicated that the samples be- come less anisotropic with increasing effective stress. the shallow water facies (Figure 2. represents approximately the palaeo-geomorphology on which deposition of the calcareous marl occurred). Their colour is generally yellowish grey to grey, mainly massive, but sometimes laminated, occasionally organic rich bands, layers of few – ten centimetres scale carbonate rich sandstone layers are interbedded with the monotonous calcare­ ous marl in the core samples (Figure 3). Article history: Manuscript received January 31, 2017 Revised manuscript accepted May 26, 2017 Available online June 28, 2017 Keywords: calcareous marl, acoustic anisotropy, pore structure Figure 1. A typical log pattern of the Tótkomlós Calcareous Marl Member as re- corded in a Battonya borehole. 1. calcareous marl, 2. limestone, 3. marl, 4. meta- morphic basement, 5. clay-marl, 6. core sample with the measured carbonate content (modified after MAGYAR et al., 2001). G eo lo gi a C ro at ic a Geologia Croatica 70/288 2. THEORETICAL BACKGROUND OF THOMSEN’S ANISOTROPY PARAMETERS In the case of elastic materials, the relationship between stress and strain can be described by Hooke’s law (auLD, 1973): s eij ijkl klC= (1) where: s – the second-rank stress tensor e – the second-rank strain tensor i, j and k, l = 1, 2, 3 indicating one of three orthogonal axes C – the fourth-rank stiffness tensor with 81 tensor compo­ nents Symmetry of the stress and strain tensors, as well as internal energy arguments lead to reduction of the number of indepen dent elements to 21 (HEMSINg, 2007). Elastic properties of shales and shale-type rocks (as well as calcareous marls) are known to be anisotropic, typically consi- dered to have a Transversely Isotropic (or hexagonal) symmetry with a symmetry axis perpendicular to the bedding plane. The main sources of anisotropy are layering (TI symmetry), the pre­ ferred orientation of clay minerals, microcracks and pores. In the case of vertically transversely isotropic (vTI) mater- ials with hexagonal symmetry, the stiffness tensor could be de­ scribed by 5 independent tensor components. C C C C C C C C C C C C CIJ = − − 11 11 66 13 11 66 11 13 13 13 33 44 2 0 0 0 2 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 44 66 C C                     (2) Elastic constants (C11, C33, C44, C66, C13) that describe a vTI medium as a function of phase velocities are given by THoM­ SEN (1986): C VP11 90 2= r , (3) C VP33 0 2= r , (4) C VS44 0 2= r , (5) C VSH66 90 2= r , (6) C C V C C C C CP 13 44 45 2 11 33 44 11 33 1 2 2 4 = − + − − − − −        (4 ) ( )2 2r , / (7) Figure 2. Basal Pannonian (Base Calcareous Marl Member) TWT seismic depth map of the study area. The investigated well is highlighted with the red circle. Calcareous marl is non-deposited or eroded at the reddish coloured acreages (LEMBERKOVICS, 2016). Figure 3. Core photo and its macroscopic description of the investigated sediments. Various rock types or subfacies were identified in the calcareous marl sequence (LEMBERKOVICS, 2016). The colour code of the rectangles highlights the same rock facies. where vP,90, vP,0, vS,0 and vSH,90 are the compressional and shear wave velocities measured in the different orientations according to Figure 5. and r is the bulk density of the sample. Thomsen’s anisotropy parameters (e, g, d) which are usually used in depth migration in seismics are dimensionless combina­ tions of elastic moduli that characterize transversely isotropic materials. They are defined in the following way: e = −C C C 11 33 332 (8) g = −C C C 66 44 442 (9) d = + − − − ( ) ( ) ( ) 2 2C C C C C C 13 44 33 44 33 33 442C (10) where e is the P-wave anisotropy factor, g is the SH-wave aniso- tropy factor and δ is a more complicated expression that may be interpreted as a measure of the anellipticity of the P wave curve (CHoLaCH & SCHMITT, 2006). as THoMSEN (1986) de­ G eologia C roatica Koroncz et al.: Pore characterization, acoustic and permeability measurements on core plug triplets, from the Miocene Tótkomlós (Calcareous Marl) ... 89 scribed, δ is the difference between the smallest offset normal moveout (NMo) velocity and vertical velocity to interpret the small avo response (TSvaNkIN, 1997, oSTaDHaSSaN et al., 2012). Depth migration algorithms in seismic processing require Thomsen’s e and d values, and these parameters can reliably be calculated from direct laboratory measurements. In the case of sedimentary rocks, values of e, g and d are generally in the same order of magnitude, usually smaller than 0.2. Moreover, the va- lues are typically higher than 0 in most rocks but negative values are also possible (THoMSEN, 1986). 3. SAMPLES Three plug-triplet samples were taken from a homogenous part of the calcareous marl core sections and are macroscopically characterised in Table 1. Table 1. Macroscopic description of investigated core sections. S-01 Light brown calcareous marl bands (thickness 1-3 cm). Transitions from the lighter bands to the darker ones are continuous. S-02 Light yellowish-greyish brown calcareous marl with disperse organic matter. S-03 Greyish brown, homogenous calcareous marl. 4. TEST METHODS an ultrasonic pulse-transmission technique was used to deter­ mine the P- and S-wave velocities. In this case two transducers are placed on the end faces of the samples. The frequency of the transducers used for measurements was centred around 1 MHz, both for compressional and shear waves. Travel times for velo city data were determined with a „first-break” record (Figure 4). For an accurate and exact determination of flying times, a new algo­ rithm was developed based on the akaike Information Criteria (akaIkE, 1973, koroNCz & FEDor, 2016). Previous stud­ ies (e.g. DELLINgEr & vErNIk, 1994; HorNby, 1998) con­ firmed that ultrasonic laboratory measurements represent phase velocities. For the analysis of acoustic anisotropy, a multicore test method was used. In this case, multiple 1” diameter plugs were drilled parallel (90°), perpendicular (0°) and at 45° to bedding planes from a larger sample (Figure 5). Multi-core methods have been used by many authors to quantify anisotropy (e.g. JoHNSToN & CHrISTENSEN, 1995; vErNIk & LIu, 1997; HEMSINg, 2007) however there are dis­ advantages of the multi-core method. For example, coring may sample different heterogeneities of the material (HEMSINg, 2007). Figure 4. A typical record of compressional (above) and shear (below) waveforms with first-break detection using a modified Akaike Information Criteria (SAIC) algorithm. Figure 5. Schematic figure of multicore sampling with the polarization direc- tions of shear wave transducers used for velocity measurements (MARTÍNEZ & SCHMITT, 2011). For the investigation of the effect of effective stress, mea- surements were carried out in a high pressure triaxial cell in the range of 0-50 MPa effective stress using 5 MPa pressure incre­ ments. Plugs were dried at 60°C to preserve the chemically bound water in the lattice of clay minerals, and were stored in a desic­ cator before measurement. In order to prevent contamination of samples, coupling media were not used. Due to the low coupling, measurement error is slightly higher at 0 Mpa effective stress. at higher confining pressures, the rMS (root-mean-square) error of velocities was less than 1%. For investigating the stress sensitive behaviour of porosity, „boyle’s law Single Cell Method“ (aPI, 1998) was used to deter­ mine porosity using 10 MPa pressure steps. Helium gas perme­ abilities were also measured with the „Pressure-Falloff, axial gas Flow measurement technique” (aPI, 1998) at a constant pres­ sure of 15 MPa. G eo lo gi a C ro at ic a Geologia Croatica 70/290 For pore structure measurements, plugs with different dia- meters (1”, 9 mm) were studied. Every sample was dried at 60 °C in a drying oven for a week prior to measurement. The skeletal volume of plugs with 1” diameters were mea- sured at 25 °C with a Pentapyc 5200e (PPy-30T) instrument. knowing the geometric volume and the skeletal volume, the he­ lium porosity (φHe) of the plugs was calculated. Skeletal density (ρskeletal) and bulk density (ρbulk) were also calculated from the measured data. The N2 and Co2 adsorption-desorption isotherms of small plugs (9 mm diameter) were measured with a Quantachrome au­ tosorb aS-1 instrument. after the drying process the samples were outgassed for 24 h at 60 °C before measurement. Physisorp­ tion tests were carried out with N2 (at 77 k) and Co2 (at 273 k) gases, in the case of N2 gas bET (brunauer-Emmett-Teller method) surface area (SbET), total pore volume (vtot) parameters were calculated. The Co2 gas is very useful for exploring the mi­ cropore structure of the samples (CazorLa-aMoroS et al., 1996), the micropore surface area (Smic) and micropore volume (vmic) were determined with the Dubinin-radushkevich method. a Quantachrome Poremaster-60 gT instrument was used during the mercury intrusion porosimetry (MIP) tests. The maxi- Figure 6. Acoustic wave velocities of cores with different orientations at Sample-01 and -02. Figure 7. Calculated elastic constants and anisotropic parameters as a function of effective stress. G eologia C roatica Koroncz et al.: Pore characterization, acoustic and permeability measurements on core plug triplets, from the Miocene Tótkomlós (Calcareous Marl) ... 91 mum (420 MPa) and minimum (1.5 kPa) applied pressure corre­ sponds to cylindrical pore sizes of 3 nm and 800 μm, respectively. Contact angles were determined with a mercury contact anglo- meter (Quantachrome) as an average value of three independent measurements. MIP was used to investigate the pore throat size distribution in the macro- and mesopore region as well as the mer­ cury porosity (φHg) and surface area (SHg) of the samples. 5. RESULTS AND INTERPRETATION Figure 6. shows acoustic velocity measurements as a function of confining pressure. all of the cores show an increase in velocity with increasing pressure. This is the result of self-sealing, as pores and microcracks begin to close when a sample is pressu- rized. Plugs drilled parallel to bedding show higher velocities than cores drilled perpendicularly to bedding. The velocities of 45° cores were between these values. velocities perpendicular to bedding increase more with pres­ sure, which could suggest that pores in the samples are aligned horizontally (MarTÍNEz et al., 2012). Figure 7. presents the calculated elastic constants and aniso- tropy parameters as a function of effective stress in the case of Sample-01 and Sample-02. Weak anisotropy is observed in both samples, Sample-01 showing slightly higher anisotropy values compared to Sam­ ple-02. For Sample-01 compressional wave anisotropy (e) mea- sured at 5 MPa is up to 11.5%, while shear wave anisotropy (g) of 6.5% was observed. both P- and S-wave anisotropy (e, g) de­ creased with increasing confining pressure, so indicating that the samples become less anisotropic with increasing effective stress. In the case of the third plug-triplet (Sample-03) higher ani­ sotropy could be observed (Figure 8). although macroscopically no difference could be observed between the samples, the much lower porosity (5.07 %) in the case of sample S-03-90° (drilled parallel to bedding planes) indicated that plug-triplets may not represent the same rock type. The explanation of this phenome­ non requires a more detailed petrological characterization (in­ cluding thin section analysis, XrD, DTa) in the future. gas permeability measurements indicate very low permea­ bility values in the order of magnitude between 10-16-10-18 m2. Since pressure fall-off permeability data are less reliable in such Figure 8. Acoustic velocity measurements of Sample-03. Figure 9. Porosity reduction due to increasing effective stress in the case of samples perpendicular to bedding. Figure 10. Typical CO2 (left) and N2 (right) adsorption-desorption isotherms of plugs. G eo lo gi a C ro at ic a Geologia Croatica 70/292 a low permeability range, these measurements will be measured in a future work by a pressure pulse decay (PPD) technique that provides more accurate results in the case of permeabilities lower than 10-15 m2 (FEDor et al., 2008). The decrease of porosity is clearly observed during measure­ ments as Figure 9. shows. Initial porosity values of tested samples vary between 10-15 % (Table 2). Porosity reduction due to in­ creasing effective stress can reach 10% in the 50 MPa range. Table 2. Test results of He-pycnometry measurements. sample ID ρskeletal [g/cm3] ρbulk [g/cm3] φHe [%] S-01 2.6591±0.0072 2.259±0.013 15.06±0.71 S-02 2.6524±0.0043 2.296±0.012 13.43±0.45 S-03* 2.6573* 2.433* 10.16* *remark: only two samples were measured The adsorption isotherms of the samples can be classified in the Type II IuPaC (International union of Pure and applied Chemistry) class, which is related to macroporous materials (Figu re 10). Hysteresis was observed on all isotherms. The shape of the hysteresis loop is similar with H3 type according to the IuPaC classification. The H3 loop is often associated with Type II isotherms and is related to non-rigid aggregates of plate-like materials (THoMMES, 2015). In addition, macropores that are not completely filled with pore condensate also caused such a hysteresis loop. The main pore throat size is between 75-110 nm according to the mercury porosimetry measurements. Calculated Hg-porosity (φHg) values are lower than the He- porosity values (φHe) indicating micropores (< 2 nm pore diameter) in the pore structure. besides, the SbET values are also higher than SHg which also suggests the presence of micropores (Table 3). Table 3. Test results of pore structure investigations. sample ID Physisorption Hg-porosimetry N2 CO2 SBET [m2/g] Vtot [cm3/g] Vmic [cm3/g] Smic [m2/g] SHg [m2/g] PSDmax [μm] S-01 7.74 4.26E-02 2.82E-03 8.10 4.66 9.33E-02 S-02 5.60 3.38E-02 2.04E-03 5.84 4.24 7.79E-02 S-03 7.61 3.99E-02 2.93E-03 8.40 4.76 1.09E-01 6. CONCLUSIONS a detailed petrophysical investigation was performed on 3 plug- triplets drilled from the Tótkomlós Calcareous Marl section, in­ cluding acoustic anisotropy and complex pore structure chara- cterization. using the multicore test method, acoustic wave velocities were measured at different confining pressures. based on accurate velocity measurements on plugs with different orientations, ani­ sotropy could be described in a quantitative way. These parameters could be applied in seismic depth migration. Porosity measure­ ments showed porosity reduction as a function of increasing effec­ tive stress. Pore structure investigations indicated plate-like mate­ rials with a typical pore throat size between 75-110 nm. Micro pores (< 2 nm pore diameter) were also observed in the pore space. In conclusion, self-sealing in marl could be observed during the measurements. PPD permeability measurements under rese- rvoir conditions parallel with ultrasonic velocity measurements would help to track and understand the process and hence petro­ physically characterize marls of sedimentary basins. REFERENCES akaIkE, H. (1973): Information theory and an extension of the maximum likelihood principle.– In: PETrov, b.N. & CSakI, F. (eds.): Second international sympo­ sium on information theory. academiai kiado, budapest, 267–281. aMErICaN PETroLEuM INSTITuTE (1998): recommended Practices for Core analysis – aPI Publishing Services, Washington. auLD, b.a. (1973): acoustic fields and waves in solids.– John Wiley & Sons Inc, New york, 423 p. CazorLa-aMoroS, D., aLCaNIz-MoNJE, J. & LINarES-SoLaNo, a. (1996): Characterization of activated Carbon Fibers by Co2 adsorption.– Langmuir, 12. 2820–2824. doi: 10.1021/la960022s CHoLaCH, P.y. & SCHMITT, D.r. (2006): Intrinsic elasticity of a textured transverse­ ly isotropic muscovite aggregate: Comparisons to the seismic anisotropy of schists and shales.– Journal of geophysical. research, 111, 410–427. doi: 10.1029/2005Jb004158 DELLINgEr, J. & vErNIk, L. (1994): Do traveltimes in pulse-transmission experi­ ments yield anisotropic group or phase velocities?– geophysics, 59, 1774–1779. doi: 10.1190/1.1443564 FEDor, F., HáMoS, g., JobbIk, a., MáTHé, z., SoMoDI, g. & SzűCS, I. (2008): Laboratory pressure pulse decay permea Laboratory Pressure Pulse Decay Perme­ ability Measurement of boda Claystone, Mecsek Mts, SW Hungary.– Physics and Chemistry of the Earth, 33, S45–S53. HEMSINg, D.b. (2007): Laboratory determination of seismic anisotropy in sedimen­ tary rock from the Western Canadian Sedimentary basin: M.Sc. thesis, university of alberta. HorNby, b.E. (1998): Experimental laboratory determination of the dynamic elastic properties of wet, drained shales.– Journal of geophysical. research, 103, 29945–29964. doi: 10.1029/97Jb02380 JoHNSToN, J.E. & CHrISTENSEN, N.I. (1995): Seismic anisotropy of shales.– Jour­ nal of geophysical. research, 100, 5991–6003. doi: 10.1029/95Jb00031 JuHáSz, gy. (1992): a pannóniai (s.l.) formációk térképezése az alföldön: elterjedés, fácies és üledékes környezet [Pannonian (s.l.) lithostratigraphic units in the Great Hungarian Plain: distribution, facies and sedimentary environment – in Hungar­ ian].– Földtani közlöny, 122/2–4, 133–165. JuHáSz, gy. (1998): a magyarországi neogén mélymedencék pannóniai képződményeinek litosztratigráfiája [Litostratigraphy of Pannonian Formations of Hungarian Neogene basins – in Hungarian].– In: JáMbor et al.: Magyarország geológiai képződményeinek rétegtana, budapest, 469–483. koroNCz, P. & FEDor, F. (2016): Experimental investigation of stress-dependent petrophysical behaviour of reservoir rocks.– In: CvETkovIĆ et al. (eds.): 8th Croatian-Hungarian and 19th Hungarian geomathematical congress, Trakošćan, 83–87. LEMbErkovICS, v. (2016): Distribution and heterogeneity of Tótkomlós Calcareous Marl Member of Endrőd Marl Formation, kiskunhalas & kelebia Exploration are- as – Mészmárga ankét, Szolnok – verbal presentation, 11 slides. Magyar, I., JuHáSz, gy., SzuroMINé korECz, a. & SÜTŐNé SzENTaI, M. (2004): a pannóniai Tótkomlósi Mészmárga Tagozat kifejlődése és kora a batton­ ya-pusztaföldvári-hátság környezetében [The Tótkomlós Calcareous Marl Member of the Lake Pannon sedimentary sequence in the battonya-Pusztaföldvár region, SE Hungary – in Hungarian].– Földtani közlöny, 133/4, 521–540. Magyar, I., raDIvoJEvIĆ, D., SzTaNÓ, o., SyNak, r., uJSzáSzI, L. & PÓC­ SIk, M. (2013): Progradation of the paleo-Danube shelf margin across the Pan­ nonian basin during the Late Miocene and Early Pliocene.– global and Planetary Change, 103, 168–173. doi: 10.1016/j.gloplacha.2012.06.007 MarTÍNEz, J.M. & SCHMITT, D. (2011): Investigating anisotropy in rocks by using pulse transmission method.– Canadian Society of Exploration geophysicists re­ corder, 36/10, 38–43. MarTÍNEz, J.M., SCHMITT, D.r. & koFMaN, r. (2012): anisotropy measurements in a multi-faced core sample by using pulse transmission method, geoConvention, Calgary, alberta. oSTaDHaSSaN, M., zENg, z., JabbarI, H. (2012): anisotropy analysis in Shale using advanced Sonic Data – Extended abstract, bakken Case Study, aaPg an­ nual Convention and Exhibition, uSa. THoMMES, M., kaNEko, k., NEIMark, a.v., oLIvIEr, J.P., roDrIguEz-rEI­ NoSo, F., rouQuEroL, J. & SINg, k.S.W. (2015): Physisorption of gases, with special reference to the evaluation of surface area and pore size distribution (IuPaC Technical report).– Pure and applied Chemistry, 87, 9–10, 1051–1069. doi: 10.1515/pac-2014-1117 THoMSEN, L. (1986): Weak elastic anisotropy.– geophysics, 51, 1954–1966. doi:10.1190/1.1442051 TSvaNkIN, I. (1997): reflection move-out and parameter estimation for horizontal transverse isotropy.– geophysics, 62, 614–629. doi: 10.1190/1.1444170 vErNIk, L. & LIu, X. (1997): velocity anisotropy in shales: a petrophysical study.– geophysics, 62, 521–532. doi: 10.1190/1.1444162