Heterogeneous Rock Research Articles

Mechanical Characterization of Main Minerals in Carbonate Rock at the Micro Scale Based on Nanoindentation

The mechanical characterization of carbonate rock is crucial for the development of a hydrocarbon reservoir and underground gas storage. As a kind of natural composite material, the mechanical properties of carbonate rock exhibit multiscale characteristics. The macroscopic mechanical properties of carbonate rock are determined by the mineral composition and structure at the micro scale. To achieve a mechanical investigation at the micro scale, this study designed a scheme for micromechanical characterization of carbonate rock. First, scanning electron microscope observation and energy dispersive spectroscopy analysis were combined to select the appropriate micromechanical test areas and to identify the mineral types in each area. Second, the selected test area was positioned in the nanoindentation instrument through the comparison of different-type microscopic images. Finally, quasi-static nanoindentation was carried out on the surface of different minerals in the selected test area to obtain quantitative mechanical evaluation results. A typical carbonate rock sample from the Huangcaoxia gas storage was investigated in this study. The experimental results indicated apparent micromechanical heterogeneity in the carbonate rock. The Young’s modulus of pyrite was over 200 GPa, while that of clay minerals was only approximately 50 GPa. In addition, the proposed micromechanical characterization scheme was discussed based on experimental results. For minerals with an unknown Poisson’s ratio, the maximum error introduced by the 0.25 assumption was lower than 15%. To discuss the effectiveness of the nanoindentation results, the characterization abilities constituted by lateral spatial resolution and elastic response depth were analyzed. The analysis results revealed that the nanoindentation measurement of clay was more susceptible to influence by the surrounding environment as compared to other kinds of minerals with the experimental setup in this study. The micromechanical characterization scheme for clay minerals can be optimized in future research. The mechanical data obtained at the micro scale can be used for the interpretation of the macroscopic mechanical features of carbonate rock for the parameter input and validation of mineral-related simulation and for the construction of a mechanical upscaling model.

Numerical Investigation on Microcrack Propagation Behavior of Heterogeneous Rock under Ultrasonic Vibration Load

Numerical Investigation on Microcrack Propagation Behavior of Heterogeneous Rock under Ultrasonic Vibration Load

Microbial growth within porous gravity currents

The effect of microbial activity on buoyancy-driven flow within a porous layer is analysed. The input fluid provides an energy source for the growth of biofilms on the porous rock. At each location within the porous layer, the porosity and permeability begin to decrease once the input fluid has invaded. This leads to an evolving rock heterogeneity that depends on the passing time of the input fluid. Hence, the evolution of the flow is partly controlled by its own history. We present an axisymmetric gravity current model, accounting for this effect. In general, a reduction in permeability leads to the flow having a lesser extent in the radial direction and greater thickness (extent in the cross-flow direction), whilst a reduction in porosity has negligible effect on the thickness but leads to a much greater radial extent. The flow is fastest near the free surface where the permeability is greatest. In the case where the porosity and permeability reduce as power-law functions of fluid residence time, the evolution of the flow and the rock properties are self-similar. Consumption of the input fluid by the microbes is also incorporated in the model and it generally leads to flows with lesser radial extent but little change in the thickness. The three impacts of microbial growth (volume loss owing to consumption and the reduction in permeability and porosity) each influence the flow in substantially different ways and the interplay is analysed. A motivation of the study, the underground storage of hydrogen, is briefly discussed.

Adaptive Mesh Generation and Numerical Verification for Complex Rock Structures Based on Optimization and Iteration Algorithms

ABSTRACTThe modelling of rock structure is of great significance in characterizing rock characteristics and studying the failure laws of rock samples. In order to construct a high‐fidelity model of the rock structure efficiently, this paper proposes an adaptive mesh dissection algorithm based on the Voronoi structure. Image processing techniques, including greyscale, threshold segmentation and edge detection, are applied to simplify the original rock image into a feature edge image. Then, a probability density diagram of the feature image is generated, which provides a probabilistic basis for the subsequent spreading of mesh seed points. Moreover, the concept of polygonal representation rate and the mesh quality evaluation system of four‐dimensional metrics are established to suggest values for the seed point parameters of the initial mesh. The initial mesh is continuously optimized and iterated by barycentric iteration and gradient descent optimization methods to form mesh structural models with high representational performance efficiently. The model tests on particle, fracture and multi‐phase rock images show that the optimized mesh model is highly similar to the original image in terms of similarity and edge fit, and the algorithm significantly reduces the short‐edge rate and improves the shape regularity of the mesh structure. Finally, numerical tests of uniaxial compression are carried out based on the optimized mesh model. The results show that the model has computational potential in numerical calculations. This method builds a procedural structure from digital images to numerical models, which can provide a reliable model basis for simulating the physico‐mechanical behaviour of heterogeneous rocks.

Investigation on the transmission of Love waves due to an impulsive line source in a heterogeneous double porous rock structure

The region around the reservoir sometimes encounters the seismic response of ground. The phenomena of Love waves propagation also caused by seismic activity which can further significantly influenced by an impulsive line source present in the rock structures of ground. The current study elucidates the propagation behavior of Love waves induced by an impulsive line source within the anisotropic layered heterogeneous double porous rock structures commonly found in proximity to the reservoir areas. The expression of the dispersion relation for the Love waves propagation due to an impulsive line source has been derived by using the Fourier transformation and Green’s function technique. The above-mentioned rock composite structure is consisted of two different portions: upper transversely isotropic double porous (TIDP) layer medium and lower heterogeneous isotropic double porous half-space (IDP) medium. In the aforesaid half-space medium, both types of heterogeneities viz. linear and quadratic have been considered. The impacts of heterogeneity parameters associated with the lower heterogeneous IDP half-space, porosity parameter (corresponds to both upper TIDP layer and lower IDP half-space) and anisotropic parameter of TIDP layer on Love waves phase velocity have also been studied. Some special outlines have also been depicted graphically.

Study on the Mechanical Properties and Failure Characteristics of Heterogenous Shale Based on CT Scanning

Summary Shale reservoirs, as a significant type of unconventional reservoir, have always been a focal point in oil and gas exploration and development. The precise determination of shale mechanical properties is fundamental to the stimulation of shale oil and gas reservoirs. The heterogeneity of rock has a significant impact on its mechanical properties. Computed tomography (CT) scanning technology is an important method for observing the internal microstructure of rocks, and digital cores constructed based on CT scans can truly reflect the heterogeneity of shale. Numerical models of shale were established using image processing technology; the basic mechanical parameters of minerals were obtained through nanoindentation experiments; and mineral content was determined by scanning electron microscopy (SEM) and energy dispersive spectroscopy (EDS). Uniaxial and triaxial compression simulations were conducted to study the impact of mineral composition and porosity on shale mechanical performance. The results indicate that the mechanical properties of shale are the outcome of the combined effects of pore distribution, mineral arrangement, and porosity. Initial natural pores significantly influence the initiation and expansion of fractures during the loading process. For models with obvious through-going joints, fractures mainly expand along the joint planes. For models with uneven pore distribution, fractures start at the pores and expand along the loading direction, eventually connecting different pores, leading to failure. In cases where a certain type of mineral is abundant or concentrated in the mineral composition, its mechanical properties will be significantly influenced by that type of mineral. In this simulation model with a high quartz content, the direction of fracture propagation during fracture was altered by the quartz. Porosity also has a significant impact on mechanical properties. As porosity increases, the model’s compressive strength decreases. Under triaxial loading conditions, at lower confining pressures, the model primarily fails due to tensile stresses; as the confining pressure increases, the proportion of tensile failures decreases, while the proportion of compressive failures increases. Models with a high content of quartz maintain a relatively stable proportion of tensile failures under different confining pressures. Meanwhile, dolomite in the model, due to its strong deformation capability, is better able to withstand tensile stress initially, but as loading continues, the proportion of tensile failures gradually increases. The composition of shale plays a crucial role in determining its mechanical properties, serving as a key reference for analyzing the mechanical behavior of shale.

Deep learning-based acoustic emission source localization in heterogeneous rock media without prior wave velocity information

Abstract Acoustic emission (AE) source localization is crucial for monitoring but often relies on prior information, such as wave velocity and arrival time. This study introduces a novel method for locating AE sources in rocks without such information, addressing challenges posed by heterogeneous sensor arrays. Experiments involving pencil led break (PLB) tests on sandstone cubes collected AE waveforms and their coordinates. A ResNet-50 based deep learning model was developed to correlate the time-frequency spectra of AE with PLB locations, expressed as spatial Gaussian distributions. The model, achieved a 79% prediction accuracy for AE localization in complex environments. While there is room for improvement in training data quantity and diversity, the results validate the model’s effectiveness, particularly in coal mines and tunnel engineering.

Determination of the REV size for heterogeneous rocks with different grain sizes: Deep learning and numerical approaches

Accurate determination of the representative elementary volume (REV) size plays a pivotal role in analysing the mechanical properties and failure processes of heterogeneous rocks in complex engineering environments. In this study, a novel microstructure modelling strategy (NMMS) for determining the REV size is proposed by combining deep learning and an improved phase-field method (PFM). Micro- and macroscale experiments are systematically conducted to determine the real microstructural characteristics and mechanical properties of heterogeneous rocks with different grain sizes. On the basis of this experimental evidence, geometric models of different sizes were reconstructed through deep learning to avoid the limitations of human-based methods, and an improved PFM was used for numerical calculations. These models were then employed to perform numerical tests under uniaxial loading conditions, and the coefficient of variation was introduced to determine the REV size of heterogeneous rocks with different grain sizes. The research findings indicate that the final REV size is the maximum value of the REVs defined by the evaluation properties within an acceptable coefficient of variation. At a criterion of 5% for the coefficient of variation, the REV sizes are 60 mm×60 mm, 70 mm×70 mm, and 90 mm×90 mm for fine-medium-grained (FMG), medium-grained (MG), and coarse-grained (CG) rocks, respectively. Furthermore, the REV determined by the NMMS was applied to investigate the effects of microstructure on macromechanical properties and damage evolution under triaxial loading conditions. The numerical results show that the NMMS can accurately predict the macromechanical properties and microcracking patterns of heterogeneous rocks, especially the intracrystalline cracks in feldspar, the interfacial cracks in gravel, and the “voids” of cracks in biotite. This research can provide some basic references for the optimal choice of the REV size of heterogeneous rocks.

Hydrogeological and physico-mechanical responses of physically heterogeneous waste rock to storage time and vertical pressure

Hydrogeological and physico-mechanical responses of physically heterogeneous waste rock to storage time and vertical pressure

Geochemistry of shales of middle Permian Barren Measures Formation, West Bokaro Basin, India: implications on provenance, paleodepositional and paleoclimatic conditions

Geochemical analyses of fine-grained clastic sedimentary rocks, such as shales, are significant since they portray the nature of the control factors from source to sink in a sedimentary basin, and help characterizing the heterogeneity in potential source and/or reservoir rocks. Bulk rock geochemistry (including major element oxides, trace elements and rare earth elements) of 20 shale samples from the Permian Barren Measures Formation, West Bokaro Basin, Peninsular India, is presented to decipher the provenance, paleo-weathering, paleotectonic, paleo-depositional and paleoredox conditions of the basin within the Lower Gondwana paleogeographic framework. X-ray diffraction (XRD) peaks and scanning electron microscopy (SEM) reveals an abundance of quartz, muscovite and clay minerals, viz., illite, kaolinite, etc., with less abundant feldspar and glauconite. The concentrations of the major element oxides, the trace elements and rare earth elements (REEs) of the shale samples are compared with predefined standards such as PAAS and UCC. The CIA (chemical index of alteration) and the A–CN–K values, based on the ratios of the major element oxides, signify intense chemical weathering. The ratios of Sr/Cu versus Rb/Sr and K2O/Al2O3 versus Ga/Rb, and high ΣREE values indicate a temperate paleoclimate. Ratios of the major element oxides (e.g., K2O/Na2O versus SiO2) suggest a passive tectonic setting. The Eu/Eu∗ and (Gd/Yb)n values depict predominantly post-Archean granitic and gneissic source rocks. The ratios of trace elements, viz., Th/Sc, Th/Co, Th/Cr, La/Sc and Eu/Eu∗, and the distribution pattern of the REEs (higher LREE in comparison to HREE) attest to a felsic-intermediate source from a nearby provenance. The redox sensitive trace elements, such as V and Ni, indicate a dysoxic condition that prevailed during deposition. The ratios of Fe2O3 versus MgO and Log (MgO/Al2O3) versus Log (K2O/Al2O3), and the ternary plots of Fe2O3–MgO–SiO2/Al2O3 suggest a non-marine to deltaic transitional environment. Such evidence of marine influences during the sedimentation of the Permian Barren Measures Formation confirms the fluvio-marine transgressive paleogeography inferred from sedimentological and paleontological inputs from the Gondwana Basins in Peninsular India, and opens up a scope for regional correlations of paleogeographic changes during the Middle Permian (Guadalupian) across the Gondwanaland.

Molybdenum isotope behavior during subduction zone metamorphism

The molybdenum (Mo) isotope composition (defined as δ98Mo measured per mil relative to NIST-3134) of many modern arc systems and the upper continental crust is heavier than the mantle and most subducting slab lithologies. This observation has led to a model whereby fluids leaving the slab transfer isotopically heavy Mo preferentially to the mantle wedge, leaving the residual slab isotopically lighter. We explore this model via an Mo isotope study of the metasedimentary and mélange lithologies of the Catalina Schist in California. These rocks record subduction zone metamorphism over a wide range of high-pressure/low-to-medium temperature conditions.Mo isotope compositions of the metasedimentary rocks decrease with increasing metamorphic grade, from a δ98Momean of −0.04 ± 0.13 ‰ (1σ, n = 3) for the lawsonite albite facies to −0.38 ± 0.07 ‰ (1σ, n = 5) for the epidote–amphibolite facies. The highest grade [amphibolite facies] samples show a slight uptick in δ98Mo values, with a mean of −0.19± 0.14 ‰ (1σ, n = 4). When coupled with major and trace element variations, the changes in δ98Mo appear to reflect metamorphic effects rather than sedimentary source rock heterogeneity. The positive relationship between δ98Mo and [Mo] and negative relationship with Ce/Mo argue for progressive loss of Mo with isotope fractionation during increasing P-T conditions; the reversal of this trend, seen in the increase in δ98Mo between the epidote amphibolite and amphibolite facies may reflect a subsequent partial re-fertilization by fluids carrying isotopically heavy Mo. δ98Mo values in the mélange across all grades are highly variable, ranging from −1.75 to −0.19‰, consistent with large scale mobilization of Mo with isotope fractionation. A metasomatized amphibolite block and reaction rind show similar light isotope compositions and signs of Mo depletion and transport. Together these three whole-rock data sets demonstrate pervasive open system behavior and mobility of Mo in the Catalina Schist. LA-ICP-MS Mo measurement of minerals in the differing metasedimentary metamorphic grades indicates that Fe oxides and/or hydroxides, titanite, rutile, and epidote are important reservoirs of Mo but when comparing mineral Mo content and modal abundance to whole-rock Mo concentrations, several samples were found to have “missing Mo”, something that has been observed in other Mo inventory studies. We attribute this to sample heterogeneity between the thin section and whole-rock powder scales.

Visualized analysis of microscale rock mechanism research: A bibliometric data mining approach

Rock mechanics is an indispensable discipline in diverse sectors, from resource retrieval to disaster mitigation. Diving deeper into this field, particularly into microscale rock mechanics, offers strategic insights and potential advancements for rock engineering practice. The objective of this research is to map the scientific production tied to microscale rock mechanics to date. In doing so, we perform a bibliometric analysis to look over and discuss the performance of the related literature. The conceptual fabric of microscale rock mechanics research, constituted by four central themes, was revealed and visualized through a text mining analysis. These areas include (1) modeling and simulation of fluid flow and transport within porous media, (2) characterizing fracture and failure in rocks, (3) understanding the deformation mechanism in response to geological processes, and (4) studying the mechanical properties of rocks subjected to extreme conditions. Lastly, an upcoming research agenda for microscale rock mechanics was proposed, centered on addressing three identified research gaps: (I) the integration of geological processes to characterize mineral properties, (II) the augmentation of fracture predictions achieved through multiscale modeling of rock heterogeneity, and (III) the exploitation of Artificial Intelligence technologies for anticipating complex fracturing scenarios. This comprehensive approach promises to enrich our understanding of rock mechanics and its applications.

Қазба айналасындағы тау массивінің жағдайын ескере отырып, контурға жақын топырақ жыныстарын бекіту

The maintenance of excavation workings will also require significant repair costs both before and after they are put into operation, especially in the zone of influence of stoping operations. Maintaining and increasing the volume of underground coal mining is possible only if there is a highly efficient technology for conducting and maintaining development workings. The purpose of the research was to evaluate the stability control parameters of the contours of mine workings, fixed with anchor bolting, in order to create a technology for intensive and safe excavation of excavation workings based on the identified patterns of behavior of adjacent rock masses. The idea of the approach is to use the technogenic stress-strain state to develop an effective technology for fastening the near-contour rock mass. The mechanism of deformation, displacement and collapse of rocks in a structurally disturbed heterogeneous rock mass was studied to assess the state of the rock mass around mine workings. A technology for fastening near-contour soil rocks has been developed, taking into account the state of the rock mass around the working, and the parameters for the operation of rock bolts in mines for fixing rods in workings in order to ensure the safety of mining operations in the mines of the Karaganda coal basin have been determined.

A new combined finite-discrete element method for stability analysis of soil-rock mixture slopes

PurposeThe purpose of this paper is to propose a new combined finite-discrete element method (FDEM) to analyze the mechanical properties, failure behavior and slope stability of soil rock mixtures (SRM), in which the rocks within the SRM model have shape randomness, size randomness and spatial distribution randomness.Design/methodology/approachBased on the modeling method of heterogeneous rocks, the SRM numerical model can be built and by adjusting the boundary between soil and rock, an SRM numerical model with any rock content can be obtained. The reliability and robustness of the new modeling method can be verified by uniaxial compression simulation. In addition, this paper investigates the effects of rock topology, rock content, slope height and slope inclination on the stability of SRM slopes.FindingsInvestigations of the influences of rock content, slope height and slope inclination of SRM slopes showed that the slope height had little effect on the failure mode. The influences of rock content and slope inclination on the slope failure mode were significant. With increasing rock content and slope dip angle, SRM slopes gradually transitioned from a single shear failure mode to a multi-shear fracture failure mode, and shear fractures showed irregular and bifurcated characteristics in which the cut-off values of rock content and slope inclination were 20% and 80°, respectively.Originality/valueThis paper proposed a new modeling method for SRMs based on FDEM, with rocks having random shapes, sizes and spatial distributions.

Analysis of rock-breaking mechanisms of high-voltage pulsed electric electrode bits

High-voltage pulsed electric rock-breaking technology is an innovative, green, and efficient method with substantial potential in the field of rock fragmentation. The efficiency of this technology is primarily determined by the design of the electrode bit. To investigate the impact of electrode bit design on rock fragmentation, this study developed a three-dimensional electro-rock breaking model based on the coupling of multiple fields: current field, electrostatic field, breakdown field, heat transfer field, and solid mechanics field. Using this comprehensive three-dimensional model, we conducted dynamic electrical breakdown simulations of granite, incorporating five different electrode bit structures and six degrees of rock heterogeneity. The simulation results elucidate the effects of pulsed peak voltage, granite heterogeneity H , and electrode bit structure on the efficiency of high-voltage pulsed electric rock breaking. To validate the simulation results, laboratory experiments on electro-rock breaking were performed. The experimental findings indicate that the conical electrode bit exhibited the highest rock-breaking efficiency, while the pentagonal prism-shaped electrode bit showed the poorest performance. The tip of prismatic electrodes generates a tip discharge effect; for the triangular prism, this effect often results in irregular rock fragmentation, which is detrimental to drilling efficiency. These results highlight the significant influence of electrode shape on rocks’ electrical breakdown and fragmentation. This study provides valuable insights into the engineering application of high-voltage pulsed electric rock-breaking technology.

Geochemical Features of Organo-Accumulative Soils of Subtaiga and Subtaiga-Forest-Steppe Light Coniferous Forests of Northern Mongolia

Geochemical features of organo-accumulative (Eutric Regosols (Laomic, Ochric), Cambic Someric Phaeozems (Loamic)) soils widely distributed in the soil cover of the subtaiga and subtaiga-forest-steppe light coniferous forests forming the lower boundary of the forest belt in the mountain structures of Northern Mongolia are considered. Data on the microelement composition of soil-forming rocks are given. It was found that the paragenetic association of trace elements in them is represented by Pb, Cu, Zn, Co, V, Cr, Ni, Mn, Mo, Ba, Sr, Zr and B. It was found that, compared with the average content in the lithosphere within the subtaiga and subtaiga-forest-steppe forest-growing belt, the residual and re-deposited weathering crusts of igneous rocks are enriched with Zn, Cr, Mo, B, at the same time they contain less Pb, Co, Mn, Ba, Sr, Zr. The residual and re-deposited weathering crusts of carbonate rocks are enriched with Pb, Cu, Zn, V, Cr, Sr, B, they contain little Co, Ni, Mn, Mo, Ba, Zr. Data on the morphological structure of soils, their physico-chemical and chemical properties, as well as on the content of trace elements and their radial distribution in the soils under consideration are discussed. The data obtained indicate the accumulation of most trace elements in the surface organogenic and humus-accumulative horizons of soils, which is associated with both the heterogeneity of soil-forming rocks and the influence of soil processes that cause the accumulative redistribution of elements and their deposition on organo-sorption and carbonate geochemical barriers. It is shown that the studied soils differ not only in the absolute values of trace elements involved in the biological cycle, but also in the intensity of their involvement in biogenic migration.

Depressurization-induced production of shale gas in organic-inorganic shale nanopores: A kinetic Monte Carlo simulation

Methane gas production from unconventional reservoirs, such as shale formations, is in high demand due to the global energy needs. However, maximizing production remains challenging due to the ultra-tight porosity, heterogeneity, and complex nature of shale rocks under high pressure. To address this issue, we employ a novel kinetic Monte Carlo (kMC) simulation to investigate the molecular-level behavior of shale gas in both homogeneous and heterogeneous organic-inorganic shale nanopores. This approach not only provides accurate computations of shale gas under high pressure but also aims to uncover the mechanisms of shale gas storage at reservoir pressure and production during pressure drawdown. Our findings indicate that organic shale ultramicropores (pore width < 0.7 nm) contribute significantly to the highest storage capacity of shale gas but pose challenges for production capacity solely through depressurization. In contrast, the free gas zone is the primary source of shale gas production from mesoporous shales, which has a high recovery efficiency but a low production capacity due to less pronounced interaction effects from the shale surface. The heterogeneous nature of shale surfaces leads to asymmetric distributions of density and potential energy across pore widths, with methane molecules favoring locations near organic pore walls due to stronger attractive interactions, while inorganic pore walls facilitate shale gas migration. Interestingly, the optimal ultramicropore size yields the highest recovery efficiency of shale gas via depressurization, characterized by a transition from near-commensurate to incommensurate molecular packing between reservoir and post-drawdown pressures. Based on these detailed molecular simulations, further research is necessary to develop innovative techniques for enhancing shale gas recovery, especially in micropores.

Failure mechanism of a tunnel in a soil–rock mixture based on a new combined finite–discrete element method

The mechanical properties and failure characteristics of soil–rock mixtures (SRMs) directly affect the stability of tunnels constructed in SRMs. A new SRM modelling method based on the combined finite–discrete element method (FDEM) was proposed. Using this new SRM modelling method based on the FDEM, the mechanical characteristics and failure behaviour of SRM samples under uniaxial compression, as well as the failure mechanism of SRMs around a tunnel, were further investigated. The study results support the following findings: (1) The modelling of SRM samples can be achieved using a heterogeneous rock modelling method based on the Weibull distribution. By adjusting the relevant parameters, such as the soil–rock boundaries, element sizes and modelling control points, SRM models with different rock contents and morphologies can be obtained. (2) The simulation results of uniaxial compression tests of SRM samples with different element sizes and morphologies validate the reliability and robustness of the new modelling method. In addition, with increasing rock content, VBP (volumetric block proportion), the uniaxial compressive strength and Young’s modulus increase exponentially, but the samples all undergo single shear failure within the soil or along the soil–rock interfaces, and the shear failure angles are all close to the theoretical values. (3) Tunnels in SRMs with different rock contents all exhibit X-shaped conjugate shear failure, but the fracture network propagation depth, the maximum displacement around the tunnel, and the failure degree of the tunnel in the SRM roughly decreases via a power function as the rock content increases. In addition, as the rock content increases, such as when VBP = 40%, large rocks have a significant blocking effect on fracture propagation, resulting in an asymmetric fracture network around the tunnel. (4) The comparisons of uniaxial compression and tunnel excavation simulation results with previous theoretical results, laboratory test results, and numerical simulation results verify the correctness of the new modelling method proposed in this paper.

Potential of polymer’s viscosity and viscoelasticity for accessible oil recovery during low salinity polymer flooding in heterogeneous carbonates

Despite the usual improvement in sweep efficiency during polymer flooding, the polymer’s viscoelastic characteristics were reported to contribute to residual oil recovery in sandstone reservoirs at high flux conditions. In carbonate reservoirs, besides high-salinity, high-temperature, and non-water-wet nature, heterogeneity is a key variable that may impair the viscoelastic fluid flow in the porous media and recovery performance. However, no attempts were made to study the polymer’s viscoelastic influence on residual oil recovery in heterogeneous carbonate formations and the objective of this paper is to systematically examine the recovery performance of high molecular weight viscoelastic sulfonated polymer (AN-125-SH) in comparison to viscous xanthan gum polymer with the low salinity injection water of 5957 ppm total dissolved solids (TDS) in heterogeneous outcrop carbonate rocks at fluxes ranging from shear thinning to shear thickening using the combination of bulk shear rheometry, NMR, single-phase in-situ studies and two-phase core flooding experiments.Nuclear magnetic resonance (NMR) imaging ensured that the chosen cores are characterized with similar pore-size distribution. The steady-shear rheological experiments showed that 2500 ppm AN-125-SH polymer exhibits a lower viscosity at the lower shear rate and an early shear thickening (indicating higher non-linear viscoelasticity) than 2000 ppm xanthan gum. Then single-phase flow experiments are conducted using two polymer systems to determine the desired flux rates before the eventual multi-rate two-phase core flood experiments. Core flooding results revealed following an extensive secondary water flooding, at the shear thinning flux of 4 ft/day, 2000 ppm high viscous xanthan gum injection leads to an incremental recovery factor of 8.3 % with the generated pressure gradient of 20.48 psi/ft compared to 3.6 % recovery obtained using 2500 ppm AN 125-SH polymer with a pressure gradient of 16.2 psi/ft. However, at shear thickening fluxes of 20 ft/day and 30.5 ft/day, the higher elastic polymer enhanced the recovery by 3.3 % and 1.2 %, respectively. Conversely, the less elastic xanthan gum showed minimal to no incremental recovery at these fluxes.The novel findings of this work convey that the viscosity of xanthan gum provides a comparative advantage over the viscoelasticity of synthetic sulfonated polymers on oil recovery at relatively low fluxes in the studied heterogeneous formation at low salinity conditions.

Deriving Flow Velocity and Initial Concentration From Liesegang‐Like Patterns

AbstractZebra rocks, characterized by their striking reddish‐brown stripes, rods, and spots of hematite (Fe‐oxide), showcase complex self‐organized patterns formed under far‐from‐equilibrium conditions. Despite their ease of recognition, the underlying mechanisms of pattern‐forming processes remain elusive. We introduce a novel advection‐dominated phase‐field model that effectively replicates the Liesegang‐like patterns observed in Zebra rocks. This numerical model leverages the concept of phase separation, a well‐established principle governing Liesegang phenomena in a two‐dimensional setting. Our findings reveal that initial solute concentration and fluid flow velocity are critical determinants in pattern morphologies. We quantitatively explain the spacing and width of a specific Liesegang‐like pattern category. Furthermore, the model demonstrates that vanishingly low initial concentrations promote the formation of oblique patterns, with inclination angles influenced by rock heterogeneity. Additionally, we establish a quantitative relationship between band thickness and geological parameters for orthogonal bands. This enables the characterization of critical geological parameters based solely on static patterns observed in Zebra rocks, providing valuable insights into their formation environments. The diverse patterns in Zebra rocks share similarities with morphologies observed on early Earth and Mars, such as banded iron formations and hematite spherules. Our model, therefore, offers a plausible explanation for the formation mechanisms of these patterns and presents a powerful tool for deciphering the geochemical environments of their origin.