Carboniferous Research Articles

Investigating the changes in the strength of carbonate rocks exposed to microwave energy

Investigating the changes in the strength of carbonate rocks exposed to microwave energy

Pore-scale fluid distribution and remaining oil during tertiary low-salinity waterflooding in a carbonate

The utilization of low-salinity waterflooding as a promising enhanced oil recovery method has exhibited exciting results in various experiments conducted at different scales. For carbonate rock, pore-scale understanding of the fluid distribution and remaining oil after low-salinity waterflooding is essential, especially the geometry and topology analysis of oil clusters. We performed the tertiary low-salinity waterflooding and employed X-ray micro-CT to probe the pore-scale displacement mechanism, fluid configuration, oil recovery, and remaining oil distribution. We found that the core becomes less oil-wet after low-salinity waterflooding. Furthermore, we analyzed the oil-rock and oil-brine interfacial areas to further support the wettability alteration. By comparing images after high-salinity waterflooding and low-salinity waterflooding, it is proven that wettability alteration has a significant impact on the behavior of the two-phase flow. Our research demonstrates that low-salinity waterflooding is an effective tertiary enhanced oil recovery technology in carbonate, which changes the wettability of rock and results in less film and singlet oil.

Experimental analysis of primary factors controlling carbonate rock dissolution capability and its impact on geothermal reservoir modification

Experimental analysis of primary factors controlling carbonate rock dissolution capability and its impact on geothermal reservoir modification

Generation and evolution of the Choiyoi granitic magmatism based on U-Pb zircon studies, Cordón del Portillo, Frontal Cordillera (Argentina)

Based on a relevant geochronological U-Pb zircon dataset (n = 47) from a sample (PBL-109) of the Cerro Punta Blanca pluton (CPB), which is part of a calc-alkaline suite, we corroborate the development of protracted magmatic activity with three major crystallization events for this Permian magmatism: 278 ± 1, 283 ± 2, and 289 ± 2 Ma, which outcrop in the Cordón del Portillo, Cordillera Frontal (CF) of Argentina. These ages can be assigned to the lower section of the magmatic record of the Choiyoi magmatism (ca., 290-265 Ma), while the age of 289 Ma represents the oldest known age for the Choiyoi, indicating the start of this magmatism during the Artinskiense. Considering these geochronological data, we postulate the presence of a deep mush reservoir where protracted magmatic activity permitted the prolonged crystallization of antecrysts (ca. 283–289 Ma). Migration of the parental magma from the mush reservoir zone occurred near the time of emplacement and culminated in the formation of an ephemeral magma chamber at shallow levels, where zircon autocrysts crystallized (ca. 278 Ma). Age spectra reported within individual samples support the idea of massive magma migration when conditions were favorable (e.g., thermally matured crust). In this view, the studied “older” Choiyoi magmatism represents a continuous magmatic event lasting 11 Ma and corresponds to a single magmatic episode rather than different periods of magmatic activity and subsequent emplacements. A later alkali-calcic magmatic event is recorded at 265 ± 4 Ma from a sample of the Cerro Bayo pluton (MH-0113), which could represent the end of the lower section of the Choiyoi magmatism. Whole-rock Sm-Nd, Rb-Sr, and Lu-Hf data in zircon, along with ages reported for the studied igneous and inherited zircon from the CPB; together with isotopic data and ages from the detrital zircon found in the Carboniferous accretionary complex of Chile, indicate that the source of the Permian parental magma in this region was a heterogeneous continental crust mainly formed by Devonian and Carboniferous rocks, related probably to a magmatic arc. However, some contribution from the Carboniferous accretionary complex of Chile to the parental magma should be consider.

The degree of soil arsenic background enrichment by carbonate weathering is mainly controlled by climate in large spatial scale

Spatial distribution of soil arsenic (As) is heterogeneous. Making clear the dominant factor(s) controlling its spatial variation contributes to the differentiation of its natural background from anthropogenic pollution. Recent studies have found that the “high background” of soil heavy metals may be induced by the process of carbonate weathering. However, the extent, process, and factors that controls the degree of this enrichment, is still unclear, especially in large spatial scales. In this study, this problem is to be revealed by a compilation and spatial data mining of soil geochemical data for Australia, North America, Europe, and China with a collection and chemical analysis of bedrock and soil profiles in a typical region. The results indicated that the process of carbonate weathering strengthened the As enrichment in the soil; therefore, the carbonate rock underlaid regions could be distinguished from those of other bedrocks in the soil As maps. In the southwest region of China, an area larger than 12,000 km2 exceeding the intervention value of 100 mg/kg could be produced according to the national standard. Iron minerals play an important role for this enrichment as they serve as carriers for As during the stages of both carbonate leaching and the decomposition of silicate minerals. Geodetector results indicated a significant dependence on climate for the degree of enrichment, and both higher temperature and precipitation were necessary, which caused a decreasing trend of soil As on carbonate region from the warm and humid south to the cold and dry north in China. The dependence of soil As concentration on the mean average precipitation (MAP) and mean average temperature (MAT) was formulated using the data of the wider geographical regions across the world, and a theoretically predicted map has been produced to show the extent caused by this causation.

Characteristics and Control Factors of a High-Quality Deeply Buried Calcareous Sandstone Reservoir, the Fourth Member of the Upper Xujiahe Formation in the Western Sichuan Basin, China

A special type of sandstone in which carbonate rock fragments (CRFs) dominate the composition developed in the Upper Triassic Xujiahe Formation’s fourth member (Xu4) in the western Sichuan Basin, known as calcareous sandstone. Calcareous sandstones are widely distributed in the western Sichuan and is the main production target of tight sandstone gas in the Sichuan Basin. In this study, thin sections, porosity–permeability testing, scanning electron microscopy, and X-ray diffraction are applied to examine the characteristics and control factors for high-quality reservoirs in the calcareous sandstones, with a view to providing guidance for natural gas exploration and development in calcareous sandstones. The results show that the calcareous sandstone belongs to litharenite, with an average framework grain composition of 30% quartz, 1% feldspar, and 69% rock fragments, while the Xu4 sandstone has a high quartz content (average content of 71%). Primary intergranular pores are the main storage space, and the reservoir quality is quite poor. Under the influence of different parent rock properties of sandstones, there are obvious differences in the composition of framework grains between the calcareous sandstone and the ordinary Xu4 sandstone, which in turn affects the reservoir storage space, diagenesis, and reservoir quality. High-energy depositional conditions, low content of late cements, and the development of fractures are the main controlling factors for the formation of high-quality reservoirs in Xu 4 calcareous sandstones.

Acoustic emission characteristics of coal and limestone failure based on MFCC

Fractures in the rock surrounding karst caves can result in substantial surface subsidence, especially in mining areas with abundant carbonate rocks. Acoustic emission (AE) is a crucial signal associated with rock fracture mechanics and lithology. This study aimed to distinguish the AE features between coal and limestone using the MFCC-based method and universal AE analysis method in conjunction with uniaxial compression tests. The results indicated that the MFCC-4–6 of AE could distinguish between coal and limestone prior to rock failure. The MFCC-1–6 of limestone exhibited a tail-like pattern change on the stress sawtooth and a step-like pattern change after reaching peak stress. The fluctuation degree of MFCC-2–5 exhibited an initial transition from small to large, followed by a subsequent transition from large to small, reaching its maximum value during the elastic deformation stage. Additionally, the power law distribution index of AE energy in the magnitude range of 102–104 for coal was larger than that for limestone. The AE events in coal were mainly characterized by shear fractures, while limestone exhibited more tension fractures than coal, as indicated by the AF/RA analysis. Furthermore, the lower the AE peak frequency, the larger the scale of fracturing in coal and limestone. These findings are valuable for ground control and early warning of rock failure.

CO2-foam flood with wettability alteration for oil-wet carbonate reservoirs

Carbon dioxide (CO2) flooding in carbonate reservoirs often suffers from poor sweep efficiency, attributed to unfavorable mobility ratio, gravity segregation, and reservoir heterogeneity. CO2-foam injection is a promising strategy to control CO2 mobility. The stability of CO2-foam is affected by reservoir wettability. This study aims to systematically investigate the effectiveness of CO2-foam flood coupled with wettability alteration in an initially oil-wet carbonate rock at a pressure below the minimum miscibility pressure. Three types of surfactants, including a nonionic surfactant Soloterra 843, a cationic surfactant CETAC-30, and a zwitterionic surfactant LMDO, were evaluated as potential foaming agents. Contact angle tests were conducted to assess the wettability alteration capability of these surfactants. Bulk foam stability tests and foam rheology tests were conducted under reservoir conditions. Finally, foam core flood experiments were performed using oil-wet carbonate core samples. Results from contact angle experiments indicated that only CETAC-30 was effective at altering rock wettability from oil-wet to water-wet by removing the negatively charged organic acid components from the rock surface. Foam stability and rheological tests revealed that LMDO exhibited highly stable CO2-foam in the absence of crude oil. However, the presence of crude oil adversely affected CO2-foam stability for all the surfactants. Core flooding experiments demonstrated that injecting wettability alteration surfactant, CETAC-30, in form of foam could not only alter the wettability of rocks, but also enhance the foam strength. Continuous gas injection (without foam) increased oil recovery by 14.4 % due to the near-miscibility of CO2 and oil. Both CETAC-30 and LMDO surfactants increased the oil recovery by about 35 % after waterflood, but the wettability altering surfactant showed the highest mobility reduction factor (MRF). These findings provide valuable insights into the complex interplay between near-miscibility, wettability alteration and foam strength in carbonate reservoirs, offering a basis for optimizing CO2-foam flooding strategies for enhanced oil recovery and carbon sequestration.

Two-episode mineralization in the Haerdaban Pb–Zn deposit, NW China: Insights from sulfide trace elements, in situ S–Pb isotopes, and Rb–Sr geochronology

The Haerdaban Pb–Zn deposit, located in the eastern West Tianshan Orogen, hosts stratiform and vein-type mineralization within Precambrian carbonate rocks. However, there has been limited research on the distinctions and relationships between these two mineralization styles. In this study, we compared trace element distributions, fluid conditions, material sources, and mineralization ages between stratiform and vein-type mineralization and reconstructed a detailed genetic model. Two episodes and four stages of mineralization were identified, including five generations of pyrite and two generations of sphalerite. The sedimentary exhalative episode represents stratiform mineralization including Stage I pyrite–pyrrhotite layers (Py-1, Py-2a, Py-2b, and Py-3) and Stage II sphalerite–galena layers (Sph-1). The magmatic–hydrothermal episode represents vein-type mineralization including Stage III pyrite–quartz–calcite veins (Py-4) and Stage IV sphalerite–galena–quartz–calcite veins (Sph-2). LA–ICP–MS analysis of Se–Co–Ni–As–Ag–Sb–Bi contents in pyrite suggests that Py-1 to Py-3 formed under relatively low temperatures (280 ± 8 °C) and high ƒO2 sedimentary conditions. Py-4 formed under relatively high temperatures (339 ± 18 °C) and low ƒO2 hydrothermal conditions. Variations in Fe–Mn–In–Ga–Ge–Cd–Cu contents in sphalerite indicate that Sph-1 formed under relatively low temperature (211 ± 7 °C), intermediate ƒS2 (logƒS2 = −11.1 ± 0.5), and moderate pH sedimentary exhalative conditions. Sph-2 formed under relatively high temperature (292 ± 5 °C), high ƒS2 (logƒS2 = −7.4 ± 0.2), and low pH hydrothermal conditions. In situ analysis of sulfide S–Pb suggests that ore-forming materials for stratiform mineralization were primarily derived from marine sediments, while those for vein-type mineralization primarily originated from magmatic sources. Sph-1 from stratiform mineralization yielded an Rb–Sr isochron age of 719 ± 16 Ma, while Sph-2 from vein-type mineralization exhibited an Rb–Sr isochron age of 380.3 ± 7.7 Ma. Random Forest classification of trace elements in pyrite and sphalerite predicted that stratiform mineralization is of sedimentary genesis, whereas vein-type mineralization is of magmatic–hydrothermal genesis. Based on our results, we identified a two-episode mineralization history for the Haerdaban Pb–Zn deposit: Neoproterozoic syngenetic sedimentary exhalative mineralization overprinted by Devonian epigenetic magmatic–hydrothermal mineralization. The results of this study highlight the importance of considering multiple geological events in ore-forming processes and will facilitate the exploration of similar deposits in the West Tianshan Orogen.

Use digital rock physics to characterize velocity and attenuation signatures of deep cavity-fracture carbonate reservoirs

Abstract Conventional rock physics experiments and well logging analyses face challenges in deep carbonate formation containing meter-scale cavity-fracture system, since it is unfeasible to extract representative samples from targeted formation. In this study, differing from the conventional digital rock models at the core scale, we construct a digital rock model from realistic outcrop images with meter-scale heterogeneity of cavities and fractures. Using a dynamic stress–strain simulation, we investigate how reservoir type, fracture density and spacing, cavity porosity and size, fluid saturation, and connectivity influence wave dispersion and attenuation. By selecting the simulation frequency of 0.3 MHz approximately corresponding to the field seismic measuring frequency of 30 Hz, the simulation results can reflect the realistic seismic velocity and attenuation characteristics of deep carbonate reservoir rocks. We find that the P-wave velocity decreases linearly with increasing fracture and cavity porosity. The P-wave attenuation is sensitive to variations in fracture density but shows limited sensitivity to changes in cavity porosity and size. Fractures with narrower spacing result in higher attenuation when fracture density remains constant. Deep carbonate rocks with cavity-fracture systems saturated with two immiscible fluids display enhanced attenuation compared to fully saturated conditions. The velocity demonstrates an almost linear increase with water saturation, while attenuation initially increases before decreasing at a peak water saturation of 70%–80%. The connectivity between fractures and cavities significantly influences the attenuation behavior. Our simulation results provide valuable guidance for seismic data processing and interpretation of deep carbonate reservoir rocks with cavity-fracture systems.

Soil organic carbon and nitrogen in a carboniferous spoil heap as a function of vegetation type and reclamation treatment

AbstractEvaluating the impact of vegetation types and reclamation methods on soil organic carbon and nitrogen in carboniferous spoil heaps is critical for selecting the best vegetation type and reclamation method to improve ecosystem services in a changing climate. This paper presents the relationship between vegetation types (woodland, forbland, and grassland) and reclamation techniques (barren rock, topsoil application, succession, and cultivation) on soil organic carbon (SOC) and total nitrogen (TN) in developing soils on carboniferous rocks in coal mine heaps. Soil samples were collected from the litter layer (Oi + Oe) and the A horizons (0–10 cm). The results revealed that vegetation types and reclamation methods significantly affected SOC and TN stocks. Woodland exhibited higher SOC and TN in the Oi + Oe horizons than other vegetation types. Topsoil application and cultivation resulted in the highest SOC and TN stocks in the A horizons (0–10 cm) under woodland and forbland compared to succession on bare carboniferous rock. In grassland, there was no significant difference in SOC stock under topsoil application and cultivation; however, significantly higher TN stock was observed in the 0–10 cm areas with topsoil application compared to succession on bare carboniferous rock. Based on the results, topsoil application is recommended to improve SOC if the mining site is restored using woodland. Conversely, grassland exhibits a similar amount of SOC stock with or without topsoil application. Considering the difficulty of obtaining topsoil, we suggest that grasses are optimal for SOC stock in the studied mining sites, followed by forbs.

Validation of Mg Isotopic Measurements for the Characterisation of Ten Carbonate Reference Materials with Ultra‐Low Mg Mass Fractions

Magnesium isotope ratios in carbonate rocks and minerals play an important role in tracing geological and biological processes. We report δ26MgDSM3 values for the first time in ten carbonate reference materials (GBW07865, GBW07114, GBW(E)070159, GBW07136, GBW07108, GBW03109a, GBW07120, GBW07214a, IAEA‐B‐7 and IAEA‐CO‐8) with calcium to magnesium mass ratios around 1300 g g−1 and ultra‐low Mg mass fractions of 295 μg g−1. A combination of AG MP‐50 (100–200 mesh) and AG 50W‐X12 (200–400 mesh) resins for the matrix extraction and multi‐collector inductively coupled plasma‐mass spectrometry (MC‐ICP‐MS) were used for the isotopic measurements. This measurement procedure is applicable to materials with high and low Mg mass fractions (e.g., carbonates, lunar highlands anorthosites, ice core). It was validated using various types of reference material. In‐house synthetic solutions with 1000 and 1300 g g−1 [Ca]/[Mg] mass ratios yielded δ26MgDSM3 values at 2s = ±0.05‰ (n = 62), ‐4.89‰ and ‐4.89‰, respectively, indistinguishable from those of the pure Mg solutions, at ‐4.91‰. δ26MgDSM3 values in well characterised carbonatite and carbonate reference materials (such as JDo‐1, COQ‐1, GBW07133, GBW07217a and GBW07129), with varying MgO and CaO mass fractions, were in agreement with literature values.

Scaling-up dynamic elastic logs to pseudo-static elastic moduli of rocks using a wellbore stability analysis approach in the Marun oilfield, SW Iran

Wellbore stability analysis is a critical component of petroleum engineering, evaluating the risks of sanding, reservoir compaction, and casing failures. Laboratory rock mechanical measurements must be scaled up to reservoir scales to achieve accurate results. One challenge lies in upscaling dynamic measurements from petrophysical logs to pseudo-static elastic properties, which has significant implications for oil and gas operations. We present a novel approach that combines laboratory rock mechanical measurements with well-log data to develop a mechanical earth model (MEM) for an Iranian oilfield with over 350 wells. We conducted static elastic property measurements on 40 core samples from various layers and depths of carbonate and sandstone rocks, demonstrating the practical application of our approach. By integrating these measurements with dynamic log data and static-dynamic correlations, we established a framework for evaluating the mechanical properties of different layers. Our findings indicate that the safe mud weight window ranges from 41.5 to 118.59 pcf, while the stable mud weight window ranges from 41.5 to 156 pcf. We demonstrate the importance of conducting parallel rock mechanical studies on cores and logs to reduce uncertainties, costs, and risks during oil and gas operations. We also propose a novel methodology combining lithological characteristics, abnormally high pressure, and borehole instability mechanisms to evaluate the stability of borehole walls. This framework provides a fresh perspective on wellbore stability analysis and offers practical solutions for the industry. Essential novel techniques include developing a geomechanical model that integrates laboratory rock mechanical measurements with well-log data to evaluate mechanical properties and calculate safe and stable mud-weight windows. Our study advances wellbore stability analysis by providing a new method for addressing this long-standing challenge. It offers valuable insights for petroleum engineers working in the oil and gas industry.

Fine-Scale Lithogeochemical Features Influence Plant Distribution Patterns in Alpine Grasslands in the Western Alps of Italy.

Bedrock geology is crucial in structuring alpine plant communities. Old studies mainly focused on the compositional differences between alpine plant communities on carbonate rocks and crystalline rocks, i.e., calcareous vs. siliceous vegetation. Increasing attention is being paid to bedrock types other than calcareous or siliceous ones, viz. those which have intermediate geochemical characteristics between pure calcareous and pure siliceous ones. Among these types of 'intermediate' bedrocks, calc-schists and serpentines are generally characterized by vegetation comprised of a mixture of basiphilous and acidophilous species. We selected several sites in alpine grasslands in the Western Italian Alps, on calc-schist and serpentine bedrocks, located at 2500 ± 100 m above sea level. X-ray fluorescence quantification of major and trace elements, combined with stereomicroscopic examination of bedrock samples with a petrographic approach, revealed a much broader range of bedrock types than recognized by inspection of geological maps. The vegetation investigated in our study was mostly composed of a set of species found more or less frequently in alpine silicicolous or calcicolous plant communities of the Alps and other European mountains. The carbonate content in the bedrock was one of the main drivers of variation in grassland vegetation, not necessarily related to soil pH. There were no distinctive species uniquely characterizing grassland vegetation on serpentines or calc-schists.

3D Microporosity Characterizations in the Heterogeneous Middle Jurassic Tuwaiq Mountain Formation: Insights into an Unconventional Sweet Spot.

Carbonate microporosity can vary significantly across depositional lithofacies and cycles, owing primarily to the high degree of heterogeneity in their pore sizes, pore throat radius, geometry, and connectivity. This is further compounded by the complex diagenetic alterations during various stages of burial. In addition, the presence of micropores, which are abundant in carbonate rocks, but not visible using conventional techniques, is challenging to characterize. To address this issue, our study focused on the Middle Jurassic Tuwaiq Mountain Formation (TMF) due to its importance as an analogue to subsurface conventional and unconventional reservoirs. Here, we utilized eight samples and performed pore network modeling to quantify microporosity distribution and connectivity from high-resolution microcomputed tomography images of different microfacies (MF) in the TMF from both mud- and grain-dominated facies. These results were then validated with petrographic, SEM images, and porosity-permeability measurements. Our study revealed that, in the high-energy, grain-supported microfacies of the shallow lagoon depositional cycle, micropores dominated by interparticle and microvug types were abundant and well-connected, with mean pore and throat sizes of 7 and 4 μm. Conversely, micropores within the low-energy mud-dominated microfacies of the deep lagoon depositional cycle dominated by intraparticle and intercrystalline types were isolated and rarely connected, even at the microscale (1-4 μm in diameter, with an average of 2 μm). This result suggests that pore connectivity at the microscale is not always related to matrix porosity, and the pore connectivity is present in submicron scale, which goes against common concepts in unconventional carbonate reservoirs. Furthermore, our observations indicate that primary depositional processes play a major role in controlling the distribution and connectivity of the Tuwaiq Mountain Formation's microporosity, while diagenetic processes only have minor controls. Our study emphasizes the importance of characterizing microporosity and its connectivity in heterogeneous carbonate rocks, which may reduce the uncertainty in exploring the properties of complex carbonate reservoirs worldwide.

Experimental Investigation of Slow Relaxation of Sound Velocity in Carbonate Rock

Experimental Investigation of Slow Relaxation of Sound Velocity in Carbonate Rock

Hierarchies of stratigraphic discontinuity surfaces in siliciclastic, carbonate and mixed siliciclastic‐bioclastic tidalites: Implications for fluid migration in reservoir quality assessment

AbstractThe hierarchies of the stratigraphic discontinuity surfaces observed in ancient tidalites are qualitatively assessed, aiming to evaluate their role as possible preferential conduits for fluid migration. Three outcrop examples are presented from microtidal settings of southern Italy: (i) siliciclastic tidalites consisting of quartz‐rich cross‐stratified sandstones generated by strong two‐directional tidal currents flowing along a tidal strait; (ii) carbonate tidalites, which accumulated in a Cretaceous lagoon and tidal flat where peritidal cycles formed vertically‐stacked sequences of biopeloidal and fenestral packstones, wackestones and bindstones during repeated phases of Milankovitch‐scale sea‐level changes; (iii) mixed, siliciclastic‐bioclastic tidalites, deposited in a bay and recording offshore‐transition, to shoreface wave‐dominated and tide‐influenced environments. Observations made during this study suggest that fluid movement can be controlled by the presence of main bounding surfaces that occur at different dimensions, from large (hectometre)‐scale, to medium (decametre)‐scale, to smaller (metre)‐scales. These surfaces produced either by depositional or erosional processes, are characterised by different features and geometries in siliciclastic, carbonate and mixed siliciclastic‐bioclastic tidalites arguably revealing complex internal pathways for fluid flows. These results suggest that fluids propagating along the main discontinuities follow a dominant sub‐horizontal direction of propagation, associated with minor sub‐vertical movements, due to local internal surface geometries and interconnections and a general lack of fractures. This surface‐based approach to the study of fluid‐flow transmission within stratified rocks represents a conceptual attempt to predict fluid mobility and reservoir potential in tidalite‐bearing siliciclastic, carbonate and mixed reservoir rocks.

The Phase Change Heat of Water in the Pore Space of Rocks Based on DSC Studies.

This research investigates the phase change behavior of water within the pore space of Devonian carbonate rock samples using Differential Scanning Calorimetry (DSC) across a temperature range of -80 to 0 °C. This study focuses on dolomite and limestone samples, all with porosities below 3%, an area not extensively covered in previous literature. Significant endothermic effects were observed at temperatures below -2 °C, challenging conventional understanding. The study reveals that the latent heat of phase change in these systems can exceed 334.2 J/g, the known value for bulk water, indicating unique thermodynamic properties of water in confined spaces. For the dolomite rock sample, observed endothermic heat effects below -2 °C were 23.5% and 26.7% of total phase change energy. The cumulative pore volume calculated using the thermoporometry method was found to be higher than expected from water occupancy alone, independent of assumptions about the thickness of the adsorbed unfreezable water layer or pore shape (spherical or cylindrical). This research provides novel insights into unfrozen water content calculations, significantly enhancing frost durability assessments and the geoengineering industry.

Machine Learning-Based Prediction of Pore Types in Carbonate Rocks Using Elastic Properties

Machine Learning-Based Prediction of Pore Types in Carbonate Rocks Using Elastic Properties

Multiscale discrete fracture network modelling of shallow-water carbonates: East Agri Valley Basin, Southern Italy

We investigate the multiscale geometrical and dimensional properties of fracture and fault networks crosscutting Lower Jurassic to Cretaceous shallow-water carbonate rocks. The carbonates are exposed along the flanks of the Viggiano Mt., on the eastern side of the High Agri Valley Basin of southern Italy. Furthermore, we employ the results of field and digital structural analyses to derive the input parameters for subsequent Discrete Fracture Network modelling (DFN) of geocellular volumes whose dimension and architecture resemble those of the investigated sites. DFN modelling is carried out for geocellular volumes representing the studied bed packages (5m-side volumes), outcrops (50m-side volumes), and reservoir-scale carbonate cliffs (500m-side volumes). Notably, modelling is obtained considering the minimum mechanical aperture value obtained in the field for porosity computations and theoretical values of fracture hydraulic aperture for permeability computations. This is because the mechanical aperture values and roughness profiles collected in the field along single fractures are unreliable due to weathering and localized karst-related dissolution. Our findings reveal that a scale-dependent geometry characterizes the Cretaceous carbonates over three orders of magnitude. These results support previously published data, and contrast with those obtained for the Lower Jurassic carbonates, which suffer of some bias due to the quality of the exposures. After DFN modelling, we show that the amount of computed fracture porosity is mainly due to the SB and NSB fractures, and that equivalent permeability is greater within faulted rock volumes. In terms of horizontal permeability, which is the most reliable result after DFN modelling, we document near isotropic horizontal permeability ellipses at all scales of observations. At individual outcrops, we note that the permeability ellipses are elongated parallel to the dominant fault sets that denotes how localized strain affects the directionality of fluid flow. Finally, we summarize the original data in a conceptual model illustrating the modalities of meteoric fluid infiltration through the vadose zone, and then subsequent horizontal flow in the phreatic zone present within the fracture carbonates at shallow depths.