Diverse Scale Data for Shale Gas Formation Description—Why is Digital Shale Rock Model Construction Difficult? The Polish Silurian and Ordovician Rocks Case Study

: The aim of the study was to show that the petrophysical parameters, characterizing the shale gas formation, obtained from the various scale well logging and laboratory methods, correlated among themselves. Relationships determined on the basis of mesoscale (logs) and microscale outcomes (laboratory experiments on plugs or crumbs) were also recognizable in nanoscale in the computed tomography results. Selected logs (spectral gamma ray, resistivity, density, neutron, geochemical and acoustic logs) and laboratory methods turned out to be effective in rock typing and description of petrophysical parameters. Nanoscale results processing and interpretation was supported by the sophisticated special software poROSE (version 3.18, AGH Kraków, Poland) to determine special parameters, which correlated with the standard laboratory outcomes. Results of the mercury injection porosimetry, together with adsorption/desorption of nitrogen at 77 K and pressure decay permeability, were used as the basal parameters for building a digital model of shale rock and a detailed description of the Silurian and Ordovician shale formation, treated as the hydrocarbon prospective unconventional reservoirs. Including the computed X-ray tomography results in the correlation analyses, gave the platform to extend the standard 2D approach in building the rock model to novel, 3D and more detailed presentations of rock characteristics.


Introduction
The broad discussion on shale gas as a potential source of hydrocarbon production took place in Poland in the first two first decades of the twenty-first century. Many industrial and scientific units were engaged in prospection, completion and exploitation works in order to provide information on the industrial production from the Polish shale gas formations [1,2]. Shale gas deposits belong to unconventional hydrocarbon resources. The most important issue regarding these type of deposits, is porosity volume and pore space structure determination, which allow the understanding of the quantity and location of shale gas resources. Nowadays, unconventional resources in Poland, mostly tight gas and shale gas, are under careful and detailed considerations regarding cognitive works. There is worldwide interest in prospection of unconventional deposits, because of the necessity to increase energy production and consciousness of geological limitations regarding conventional hydrocarbon deposits.
The multidimensional analysis of the available data of diverse scale is presented in the paper. The results from the well logging and laboratory experiments were included, especially designed for shale gas formation recognition. The research purpose is to show how the various scale measurements provide indispensable information and also to show how the data mutually complete the knowledge on this very complicated geological formation. Results of the standard applied geophysical methods (seismic and well logging) were used in introduction of the specific applications. They represent macroscale and mesoscale research and data, which were considered at the very beginning of the analysis to define the potential shale gas deposits. Laboratory results were the source of the microscale data. New technologies, such as computed X-ray tomography, represent additional, very useful information on the microstructure of shale gas rocks in microscale or nanoscale. Correlations of various laboratory outcomes with well logging data, i.e., logs and results of the quantitative interpretation, were considered as an effective tool for a shale gas model construction, among other solutions.
The complex relationship between porosity and permeability made estimations of flow behavior complicated, due to heterogeneous pore networks in shale gas. Even though permeability slightly increased with porosity, permeability of samples with similar porosity differed by up to three orders of magnitude, starting with very low values. Thus, the absolute permeability prediction yielded large uncertainty. The complicated structure of the pore space in shales and the gas presence, not only in pores but also in matrix, as well as the existence of different mechanisms of gas flow, gave the assumption to use many available methods and combine the results. The relationships between connected versus unconnected pores were especially of prime importance for quantitative reservoir models.
Examples discussed in this paper originated from known and broadly tested potential Silurian and Ordovician shale gas formations from the Baltic Basin, because a large data set was gathered and accessible during the numerous industry prospection works and scientific projects [3][4][5][6][7]. Attention was focused on the combination between laboratory results, obtained on standard size samples and using the information from computed X-ray tomography (CT) as a tool for the description of shale gas model in more detailed way.

Geological Setting and Materials
Petrophysical analyses were based on the results of laboratory measurements and well logging data. Samples from the selected depth intervals of the lithostratigraphic profiles and well logging data from the same sections in three wells (Well-1, Well-2 and Well-3) were considered. Wells are sited in the Peri-Baltic Syneclise (Figure 1), north of Poland. Detailed geological description, combined with well logging information and geochemical TOC data, were considered in primary analyses. Jantar Member (Ja Mb) of the Silurian Pasłęk Formation (Pa Fm) and Ordovician Sasino Formation (Sa Fm), as potential resources of unconventional hydrocarbons, were selected for the analyses (Table 1). Silurian and Ordovician shales in Poland spread along the western margin of the East European Platform in Lublin, Podlasie and Baltic basins, extending about 700 km in length [8]. On Łeba Elevation, where the wells are located, the sedimentation of the organic-rich black shales with graptolites started in the Late Llanvirnian, reaching Wenlock [9]. Silurian Pa Fm was firstly described by [10]. An average thickness of the formation is 20-40 m, increasing from east to west and not exceeding 70 m [11]. Ja Mb is built of black bituminous claystones and belongs to the lower part of the Pa Fm. It is characterized by high oil-prone organic matter content, but the thickness of this bed does not exceed 12 m. Ja Mb commonly contains pyrite. Ordovician Sa Fm was fully described by Modliński and Szymański in the works published in 1997 and 2008 [12,13]. Lithology of this formation is mainly comprised of black, dark grey and grey-greenish bituminous shales. The thickness of the Sa Fm increases from east to west and north-west, from 3.5 to 37 m in the land part of the Baltic basin.

Methods
Well logging is a basic method conducted in the well, in order to measure the physical and chemical parameters of the rocks in situ. The measuring set can consist of various logs. Thus, the information from the well logs can cover different parameters and be the input in calculation of reservoir parameters, such as porosity and permeability. The following well logs, which are based on the different physical or chemical laws, were taken into consideration in the analysis: spectral gamma ray (uranium, thorium and potassium), resistivity logs from dual laterolog and from Xtended Micro Imager (XRMI), bulk density, volumetric photoelectric absorption index, P wave slowness, neutron porosity, geochemical log (GEM: dry elemental weight percent of selected elements and the volumetric percent of major and minor minerals). The results of the well logs interpretation were also considered: volume of fluids, i.e., free water, gas, clay water and kerogen, saturation (total water saturation, effective water saturation) and porosity (total and effective).
The following laboratory measurement results were considered in the analysis: mercury injection porosimetry (MIP), nuclear magnetic resonance spectroscopy (NMR), pulse/pressure decay permeability (PPP), absorption/desorption of nitrogen at 77 K and X-ray diffraction (XRD).
The computed X-ray microtomography (CT) results are presented in the last step.
Mercury injection porosimetry (MIP) is a technique which determines petrophysical parameters, such as porosity and pore-throat size distribution, total pore area and bulk density through measuring the intruded mercury volume under increasing pressures [14,15]. Mercury intrusion into large macropores is realized with only slight pressure, whereas greater pressure is required to force mercury into small pores. The advantage of MIP is based on the large resolution, 400 µm, down to 3 nm. Nevertheless, this method cannot provide any information about isolated pores. MIP experiments were done at the AGH University of Science and Technology (AGH UST), Faculty of Geology, Geophysics and Environmental Protection (FGGEP) using AutoPore IV 9500 (Micromeritics Instrument Corporation). The next described method, a gas sorption, is a complementary technique to the mercury injection porosimetry.
The physical adsorption of nitrogen at 77 K method can reduce the smallest size measurement to about 0.00035 mm in diameter by investigating the intraparticle microstructure. Nitrogen adsorption/desorption isotherms, obtained with an automatic gas adsorption equipment (Quantachrome Autosorb-1C automatic gas adsorption apparatus), were used to characterize the porous texture, i.e., the pore size distribution and specific surface area of extremely low pore diameters. The Brunauer-Emmet-Teller (BET), Barrett-Joyner-Halenda (BJH) and quenched solid density functional theory (QSDFT) approaches are usually used for the interpretation of nitrogen adsorption data [16][17][18][19]. For shales, with the observed mesoporosity, the nitrogen adsorption isotherms are Type IV (isotherms with hysteresis loop), according to the International Union of Pure and Applied Chemistry (IUPAC) classification [20][21][22]. The samples were degassed overnight, under vacuum at 473 K, before the physical adsorption of nitrogen measurements and other experiments. Measurements were made at the AGH UST, Faculty of Energy and Fuels (FEF).
NMR relaxometry for hydrogen nuclei was realized also at the AGH UST, FGGEP using Maran Ultra 23 MHz spectrometer [23,24]. NMR signals, from water saturated plugs, carried information about water content and pore space volume, especially total and effective (dynamic) porosity and irreducible water saturation parameter.
Pressure decay permeability was measured on crushed grains of shale. Crushing of shales allowed for access of working fluid to the pore space. Bulk density was determined on the same plugs. Measurements were done by Terra Tek-Schlumberger Reservoir Laboratory [3].
Pulse decay permeability method was also applied for selected plugs. It was extremely difficult to plug shales. Samples, in the plug form, contained no natural microfractures from the coring and processing stage of the core analysis, which enhance or artificially increase permeability values and sometimes porosity values [25]. Permeability from pulse decay method confirmed high values of this parameter in several samples. Experiment was done at the AGH UST, Faculty of Drilling Oil and Gas. The minimum and maximum values of the parameter were interpreted considering two permeability values from pressure decay and pulse decay methods, respectively.
The computed X-ray microtomography images were acquired at the AGH UST, Faculty of Physics and Applied Computer Sciences by Nanotom S (General Electric) [26]. Resolution of the images was on the level of 0.6 × 0.6 × 0.6 µm 3 (volume of 1 voxel). Segmentation was done by threshold method, using a single threshold value determined by the choice of the interpreter, comparing the results of total porosity from CT with the results from the other laboratory techniques, such as FIB-SEM (Focused Ion Beam-Scanning Electron Microscope), helium pycnometry (HeP). Interpretation of the CT images was done using poROSE (poROus materials examination SoftwarE) computer system, developed at the AGH UST, FGGEP, Department of Geophysics [27][28][29].
A prototype set for water/light kerosene immersion porosimetry (dual liquid porosimetry) [30,31] was included. TOC and other parameters from the Rock-Eval pyrolysis were also considered.
XRD analyses were done in the Institute of Geological Sciences of Polish Academy of Sciences (ING PAS), Research Centre in Krakow. Measurements were carried out according to the methodology described in detail by Środoń et al [32], while interpretation used Q-Min software (Institute of Geological Sciences of Polish Academy of Sciences, Kraków, Poland) by Marek Szczerba (IGS PAS). Measured diffraction patterns were compared to the sum of pure standards, recorded for the identified minerals. Measured and synthetic diffraction patterns were compared until the difference was minimal [33]. Mineral content (normalized to 100%) included organic matter (OM), calculated on the basis of total organic carbon content (TOC) from the Eltra CS-500 IR analyzer, as a difference between total carbon content (TC) and total inorganic carbon content (TIC) [3].
The assortment of laboratory measurements was very broad. All mentioned experiments were not able to be done on all samples, because of limited available rock material. The number of samples was limited, as regards the pressure decay permeability determination (together with bulk density) and microtomography analyses. Mercury injection porosimetry (MIP) was the basal method and its outcomes were available for all rock samples in the study and form the main group of the analyzed quantities. MIP outcomes were completed with results of the physical adsorption/desorption of nitrogen at 77 K measured in an automatic gas adsorption apparatus to extend the range of the examined pore radii. Four samples were used for the computed X-ray microtomography experiments. CT images as well as various TOC and mineral composition measures represented both Ja Mb and Sa Fm from two wells: Well-2 (three samples) and Well-1 (one sample).
The laboratory measurements listed above were done on the samples in the form of standard laboratory cylinders: 1 inch in diameter and more than 1 inch in height (e.g., pulse decay permeability), on the pieces (e.g., MIP, NMR, pressure decay permeability) or powders (e.g., XRD, HeP, absorption/desorption). Microtomographic investigations were done on very small parts of rock samples, approximately 0.5 cm in length. The idea was to carry out the measurements on the same part of material and on the homogeneous powder (pieces). MIP porosimetry is the destructive measurement, so it was done at the end of the measurement trial. Integrated well logging and laboratory analyses workflow is presented in the Figure 2, used to improve the rock typing and pore network characterization.

Results
Well logs provided continuous data along well axis, but the standard logs available in well logging data set (Quad Combo) [34] turned out to be deficient in the vertical resolution ( Figure 3) in the shale formation analysis. Nevertheless, the known and effective relationship between total organic carbon (TOC) from the laboratory measurements (for instance Rock-Eval analyses [35]) and uranium content (from the spectral gamma ray log) were successfully included, giving the basic information on parts of shales in the total geological profile that were rich in organic matter (and potentially in natural gas) ( Figure 4). Electric borehole images, provided by the modern tools, such as FMI and XRMI [36,37], are a source of detailed lithology plots, illustrating changes of thickness of claystone and mudstone laminas in gas shales. These images, compared with photos of cores, were included to present the structure of the formation (Figure 3). Moreover, the results of the geochemical log (GEM or Litho Scanner) enriched well logging data with precise information of elemental structure of rocks, helping in differentiation between various elements, especially organic and nonorganic carbon ( Figure 5). Figure 5 presents the selected standard logs in the typical way. The following logs are plotted versus depth: in the first track there is depth; in the second-spectral gamma ray SGR curves, i.e., uranium URAN (ppm), thorium THOR (ppm) and potassium POTA (%); in the third, electrical resistivity logs from dual laterolog (LLD and LLS (ohm)); in the fourth track-bulk density RHOB (g/cc), volumetric photoelectric absorption index PE (barn/electron) and density of matrix RHOMA (g/cc), P wave slowness DT (µs/ft) and neutron porosity NPHI (ls, limestone porosity units). Results of geochemical log GEM are visualized as dry elemental weight percent of Mg, Al, Si, S, K, Ca (wt%) in the fifth track. Results of the volumetric percent of major and minor minerals: illite, Fe chlorite, plagioclase, quartz, dolomite, calcite and kerogen are included in the sixth track. Next, volume of fluids, i.e., free water, gas, clay water and kerogen (%) are shown in seventh track. Saturation: total water saturation (SWT), effective water saturation (SWE (dec)) and porosity (Por) tot (total) and Por eff (effective) (dec) are presented in the eighth and ninth tracks, according to Halliburton GEM Mineralogical & Petrophysical Well Log Analysis & Interpretation (2013) [38]. It can be observed that the mostly dolomitic limestone is in the lower part of the Ordovician formation in Figure 5. The upper part is mostly interbedded siltstone, claystone, shale and marl. The core XRD data showed that the primary minerals building up the Ordovician and Silurian strata were illite, quartz, calcite, dolomite and plagioclase. Minor amounts of Fe chlorite and pyrite were also identified. Some minerals identified in other depth sections of the geological profile were not taken in consideration in this analysis. The core geochemical TOC data from the offset well was used to improve the kerogen volume estimation in the GEM interpretation.
Standard log curves are compared with the borehole wall resistivity image (sixth track) and selected microresistivity curves EDD from the XRMI and GR (seventh track), as well as the photos of the cores, in Figure 3. Typical, very thin-bedded shale formation is distinctly visible. Additionally, several to over a dozen centimeters thick laminas are clearly visible. EDD curves in comparison with GR and other standard logs show that the vertical resolution of the XRMI electric image is much better. In the first two parts of the Figure 3 it can be observed that standard logs, although providing continuous curves, are not precise enough. The depth plots of uranium content in comparison with point TOC laboratory data (Figure 4) showed the usefulness of spectral gamma log SGR for quick identification of shale formations sections rich in organic matter. In the next steps, the analysis was focused only to the potentially productive beds. Layers rich in the organic matter (kerogen) were also indicated in the comprehensive interpretation of well logs (Figure 3).   Two basic datasets were considered for the analysis of the laboratory measurement results (Tables 2 and 3).  Two plots illustrate relationships between bulk density from MIP and TOC (from Rock-Eval) ( Figure 6a) and also organic matter from XRD/TOC analyses vs. total pore volume V_tot, determined from the nitrogen adsorption/desorption at 77 K experiment (Figure 6b). Data from Figure 6 are based on the more numerous dataset (Table 2), originating from the pressure decay permeability measurements, and comprising data from Well-1 and Well-3. The next two plots present the relationships between pressure decay permeability PD_Perm (in microdarcy) and total pore volume determined from the nitrogen adsorption/desorption at 77 K experiment (Figure 6c) and total pore area from the MIP (Figure 6d). In all plots data scattering is observed. The main reasons for the data scattering are different physical phenomena, being the basis of the measurement techniques, resulting in different resolution. Heterogeneity of shale formations is the next factor complicating the relationships [39,40]. Nevertheless, correlation of data is observed, which means that the recorded quantities related to shale formation are mutually dependent and may be used to characterize the rock. The outlier observed in the Figure 6a (marked in red) (Well-1, Sa Fm) is related to sample of the highest percent of pyrite, which resulted in the increase of bulk density in the presented level of TOC. One of the outliers in the Figure 6c (marked in red) belongs to Pr Fm and consists of the extreme value of Fe-chlorite, which resulted in grain density increase ( Table 2). PD_Perm coordinate shows that chlorite provided different influence on pore space in comparison to illite/montmorillonite, of which clay minerals dominated in the considered formations. The second outlier (in orange) represents Sa Fm and is rich in quartz and organic matter. Increased volume of quartz resulted in increased permeability. Visible relations between parameters, measured in various experiments, allow expecting that similar trends will be extended to extremely detailed parameters obtained from microtomography. The plots in Figure 7 illustrate MIP characteristics for the selected samples. Quantities, i.e., cumulative pore volume and incremental pore volume, are presented in percent, in relation to the maximal value of the corresponding parameter. There is observed broad pore-throat distribution, with the peaks at about 400 nm, which are in accordance with the comparably low permeability ( Table 2). Low cumulative pore volume is consistent with the determined porosity ( Table 2). Modal shape of the incremental pore volume plots in Figure 7 reflected the pore space construction, different in both samples. One spike, located in section of small pore diameter in Figure 7b, means that nanopores influenced the porosity in this case. Several spikes in Figure 7a mean that porosity is built by pores from broad range of size, up to almost 3 µm. Different courses of the cumulative pore volume curves (blue lines in Figure 7a,b) mean different activity of pores of related radii in the samples. Both samples, Well-1(26) and Well-2(7), represented Sa Fm, which confirms its heterogeneity regarding fluid flow ability.
Graphs in Figure 8 show how adsorption/desorption extended the porosity range, explored in laboratory methods. Samples presented in the Figure 8 are the same as those presented in Figure 7. In Figure 8a,c, a nitrogen volume at standard temperature and pressure (STP) is presented as a function of the relative pressure. The volume increment per increment of pore diameter, as a function of pore diameter, is presented in Figure 8b,d. It can be observed (Figure 8a,c) that a significant amount of gas, adsorbed at a low relative pressure, is indicative of microporosity. At the relative pressure of about 0.85, the volume of gas adsorbed increased significantly as a result of the monolayer adsorption shifting to multilayer adsorption. Capillary condensation is observed at a relative pressure above 0.85. Due to the capillary condensation, the adsorption and desorption branches of the isotherm do not coincide. Thus, a hysteresis loop is formed with its shape, dependent on the pore structure. The complex shape of the hysteresis loop was probably due to the combination of several pore types. The height of shift isotherms depends on porosity. As an example, a higher porosity sample (Well-2(7), Table 3) is characterized by the bigger shift. Differences in porosity values and activity of pores of various diameters are also visible in MIP plots (blue lines). Large anomalies in Figure 8b,d indicate that large groups of small pores that were included in the gas movement network. It is important to observe, that only the density functional theory model (DFT) provides additional information about the microporosity. Higher maximum of the DFT curve in the Figure 8d plot means higher porosity of the Well-2(7) sample. The next series of plots are related to data obtained from the microtomography measurements for four samples. One of them, Well-1 (26), is a junction between two sets (Tables 2 and 3). Results of the object (pores, microfractures) size classifications for the Well-1(26) and Well-2(7) samples are presented in the Table 4. Colored illustrations of the object (pores, microfractures) size classification are included in the Figure 9a,b. The Well-1(26) sample is characterized by the higher number of pores and microfractures in comparison to the Well-2(7) sample (as evidenced by the numbers in the SUM column). The most numerous volume class is class 100-999 voxels (orange, % column) for both samples, meaning that pores gathered within the volume range of 21.6-215 µm 3 . Table 4. Object size classification. Symbols: #-number of objects in volume class; %-percentage of volume class in the pore space.

CT Analysis
Volume  Relationships between parameters obtained from microtomography [41] and other laboratory methods are presented in Figure 10. At the beginning, analysis started with the mutual relation between parameters originating only from the microtomography (Figure 10a): average coordination number (calculated on the basis of the number of branches coming out from the one node) is on the horizontal axis, while average number of the junctions in the sample is on the vertical axis. The goal was to show that the correlation between parameters is visible and only one outlier is not from the Well-1, from which there was only one sample, but from the Well-2(1), from which there were three samples (Table 3). Plots in Figure 10a,b are similar in shape. In both cases, the outlier is the same (sample of CN Aver = 3.0820) and represents the group of three samples from the Well-2, belonging to the Silurian Ja Mb. The outlier in Figure 10c belongs also to the group of samples from Well-2. It is the Well-2(4) sample. It is observed that this was the highest value of the irreducible water, calculated from NMR experiment, in comparison to the others for this sample (Table 3). Highest Sw irr, together with the lowest permeability, resulted from the highest clay and organic matter content and lowest quartz content. This was detected on the basis of XRD/TOC analyses. Feret Coefficient is the next introduced parameter from the 3D analysis of the pore space and represents ratio of the maximum horizontal diameter to maximum vertical diameter of the pore (here, average value for the sample for all detected pores). High Feret Coefficient means meaningful difference between pore radii, measured in two perpendicular directions, which is related with the lowest elongation (ratio of the shortest chord to the longest chord) in the object. Comparing plots in Figure  10c,d, where Feret Coefficient in on the vertical axis, the same outlier can be observed-Well-2(4) sample. Mineral composition of the samples from Well-1(26) and Well-2(4 and 7) ( Figure 11) could be compared on the basis of the XRD/TOC results. They differed regarding organic matter OM and the sum of clay minerals. High value of OM corresponds to low bulk density from MIP and PD_Perm, and high TOC (Rock-Eval). Samples Well-2(7) and (4) also differ among themselves as regards TOC and permeability PD_Perm. Sample Well-2(7) is characterized by the highest permeability. Main mineral components in all samples comprised quartz, orthoclase and plagioclases, calcite and dolomite, iron (pyrite/marcasite and hematite) and titanium (anatase) containing minerals, clay minerals such as an illite/montmorillonite group and chlorite. Diverse content of minerals, especially quartz (the highest value for which was found in the Well-2(7) sample) and organic matter, among others, caused the outlier position of the samples. The highest quartz amount found in the Well-2(7) sample may explain the highest permeability PD_Perm (Figure 10c), which continued the trend line, visible in the plot (decreasing of the Feret Coefficient with increased permeability). The outlier in Figure 10d, with high value of Sw irr, is related to relatively low sum of clay minerals and relatively high value of organic matter and permeability. In the presented results, there is an observed correlation with Feret Coefficient, not with Feret diameter. This means that diameters measured in different directions better describe the pore shape than diameter measured in one direction does, which means that pores have various shapes. High Sw irr, together with high PD_Perm, can be explained by water covering the surface of quartz grains (high quartz content) and identification of hydrogen in organic matter, not only in water by NMR experiment. Here, in the discussion, only pressure decay permeability was considered, but several pulse decay permeability results, especially for the Well-1(26) sample ( Table 2), revealed that not only matrix micro/nanopores should be accounted for to illustrate the complex structure of shales. High values of the pulse and pressure decay permeability showed the presence of microfractures (partially originating naturally, partially being artifacts), increasing the porosity and permeability.
Comparison of pressure and pulse decay permeability presented the range of permeability variability, from minimal value (pressure decay) to maximal value (pulse decay).
Comprehensive analysis of the groups of plots in the Figures 6 and 10 show that correlation between parameters is not only observed in the groups of the results obtained from the same scale measurements.

Conclusions
Well logging and laboratory methods, together with the comprehensive workflow, were established, utilizing a combined analysis of petrophysical results to enable the rock characterization. In rock typing SGR, XRMI and GEM logs turned out to be the most effective in quick selection of formations of the highest hydrocarbon reservoir potential in the geological profile. All laboratory experiments provided important data. Some of the methods resulted in data redundancy, so, having the limited budget for laboratory experiments, some of the methods can be rejected. In the presented shale characterization MIP and adsorption/desorption of nitrogen at 77 K showed themselves indispensable, providing information on effective porosity (total pore volume), bulk density, total pore area and pore distribution. Combined MIP and adsorption/desorption of nitrogen at 77 K covered a large pore size range, from the fracture of millimeter to nanometer scale. This combination of laboratory methods provides a novel approach to porosity/permeability calculation; including the smallest pores changed the image of the rock. Combining the outcomes of MIP and nitrogen adsorption/desorption at 77 K brought qualitative data, for better understanding of the ability of media flow in shale formation. This is a great challenge, because the estimation of flow behavior in reservoirs often suffers from uncertainties associated with pore types and pore distributions. Pressure decay permeability was the important source of true matrix permeability. Microtomography analysis allowed a representative 3D assessment of total porosity, pore connectivity and pore shape. Correlations between measured petrophysical properties and parameters from the microtomography established the platform to combine the information from the diverse scale measurements. Trends observed in the relationships between parameters, determined from the outcomes from laboratory measurements and well logging, are promising results for the newest approach in well logging and standard laboratory experiments for digital rock model construction.
The main, detailed conclusions are: • Microresistivity curves form the modern tools, such as FMI and XRMI, together with a geochemical log (GEM or LithoScanner), are a source of detailed lithology information, reflecting changes of thickness of claystone and mudstone laminas in gas shales. • Spectral gamma ray log SGR is a key log for the quick identification of shale formations sections rich in organic matter. • Correlation of the organic matter OM vs. total pore volume V_tot (nitrogen adsorption/desorption at 77 K) was observed with the correlation coefficient R 2 = 0.8, which demonstrated that the recorded parameters in the shale formation are mutually dependent and may be used to characterize the rock. • MIP quantities illustrated the nanopores' influence (up to 3 µm) on the porosity of the shale formations and broad pore-throat distribution. • Adsorption/desorption of nitrogen at 77 K gave the information that a significant amount of gas, adsorbed at a low relative pressure, is indicative of microporosity. Moreover, only the density functional theory model (DFT) delivered the additional information about the microporosity in analyzed shale formations.
• Computed X-ray tomography showed that the samples differed in the number of pores and microfractures (by almost 6-fold), but almost 50% of the pores were concentrated in the volume range of 100-999 voxels. Feret coefficient correlated with the pressure decay permeability for the analyzed samples and can be associated with the quartz content in the sample. • Diverse content of minerals (from the XRD analysis) caused the outlier position of the samples for the different correlations of laboratory measurement results.