The development of stylolites, from small-scale heterogeneities to multi-scale roughness

Stylolites are rough paired surfaces, indicative of localized stress-induced dissolution under a non-hydrostatic state of stress, separated by a clay parting which is believed to be the residuum of the dissolved rock. These structures are the most frequent deformation pattern in monomineralic rocks and thus provide important information about low temperature deformation and mass transfer. The intriguing roughness of stylolites can be used to assess amount of volume loss and paleo-stress directions, and to infer the destabilizing processes during pressure solution. But there is little agreement on how stylolites form and why these localized pressure solution patterns develop their characteristic roughness.rnNatural bedding parallel and vertical stylolites were studied in this work to obtain a quantitative description of the stylolite roughness and understand the governing processes during their formation. Adapting scaling approaches based on fractal principles it is demonstrated that stylolites show two self affine scaling regimes with roughness exponents of 1.1 and 0.5 for small and large length scales separated by a crossover length at the millimeter scale. Analysis of stylolites from various depths proved that this crossover length is a function of the stress field during formation, as analytically predicted. For bedding parallel stylolites the crossover length is a function of the normal stress on the interface, but vertical stylolites show a clear in-plane anisotropy of the crossover length owing to the fact that the in-plane stresses (σ2 and σ3) are dissimilar. Therefore stylolite roughness contains a signature of the stress field during formation.rnTo address the origin of stylolite roughness a combined microstructural (SEM/EBSD) and numerical approach is employed. Microstructural investigations of natural stylolites in limestones reveal that heterogeneities initially present in the host rock (clay particles, quartz grains) are responsible for the formation of the distinctive stylolite roughness. A two-dimensional numerical model, i.e. a discrete linear elastic lattice spring model, is used to investigate the roughness evolving from an initially flat fluid filled interface induced by heterogeneities in the matrix. This model generates rough interfaces with the same scaling properties as natural stylolites. Furthermore two coinciding crossover phenomena in space and in time exist that separate length and timescales for which the roughening is either balanced by surface or elastic energies. The roughness and growth exponents are independent of the size, amount and the dissolution rate of the heterogeneities. This allows to conclude that the location of asperities is determined by a polimict multi-scale quenched noise, while the roughening process is governed by inherent processes i.e. the transition from a surface to an elastic energy dominated regime.rn

塔里木盆地寒武－奥陶系白云岩类型、成因机理及其孔隙性研究

Cambrian-Ordovician dolostones in Tarim Basin are hydrocarbon reservoir rocks of vital importance. Under the guidance of the theories of sedimentology and the sedimentology of carbonate reservoir, based on the first-hand qualitative and quantitative data especially, combined with micro-study, geochemical and reservoir capacity analysis, and precursor research, the origin and reservoir characteristics of the dolostones were discussed. Based on detailed petrographic investigations, four types of dolostone have been recognized, which are, respecitively, mud-silt-sized dolostones, algal laminated dolostones (ALD), prophyritic dolostone, and neomorphic dolostone. Mud-silt-sized dolostones always presents as laminas together with evaporated signatures, its REE patterns and ΣREE are all close to that of the finely crystalline limestone. This kind of dolomite probably experienced relatively low fluid-rock ratio during diagenesis was formed in hypersaline and oxidizing environment and involved fast dolomitization process. It was dolomitized by evaporated seawater in sabkha environment.The main primary fabrics of algal lamination in algal laminated dolomite (ALD) can still be identified and its ΣREE (21.37) is very close to that of algae. This reveals that ALD was dolomitized during early diagenesis and algae possibly played an important role. The ALD was formed under mediation of organic matter and Mg2+ were supplied by magnesium concentrated algal laminites and sea water. Prophyritic dolostones presents mainly as patchy occurrence and yield the lowest δ13C and Z value. Its ΣREE is much less than that of the finely crystalline limestone. These characteristics reveal that the cloudy cores were dolomitized in shallow early diagenetic environments by pore fluids riched in Mg2+. Whereas the clear rims were likely formed in subsequent burial into deeper subsurface environments, and the Mg2+ needed for further dolomitization possibly was supplied by the transformations of clay minerals. Neomorphic dolostones consist of coarse, turbid crystals and exhibits sucrosic and mosaic textures. It has highest Fe2+ contents and average homogeneous temperature (110.2℃). Collectively, these characteristics demonstrate that the neomorphic dolostones was likely formed by recrystallization of pre-existing dolomites during deep burial. The ΣREE of the four types of dolostone distinctly differentiates from each other. However, their REE patterns are all enriched in LREE, depleted in HREE and have Eu negative anomalies. Its ΣREE 13.64 ppm, less than 1/4 of finely crystalline limestone, and ranks the lowest in the 4 types.These characteristics are comparable to those of finely crystalline limestone, and are mainly infuenced by the sea water. These four types of dolostone show similar REE mobility behaviour and no significant fractionation, althouth they have been subjected to evidently different diageneses. Seven main pore types are identified in the dolostones , which are fenestral, moldic, intercrystal, dissolved，breccia, dissolved breccia and stylolite pores. Fenestral pores are primary and the others are secondary. The dissolved pores and intercrystal pores are the most important reservoir spaces and followed by breccias and dissolved breccia pores, and the moldic and fenestral pores are less important. Stylolites can enhance permeability of reservoir rocks in one hand, for the other hand, the capacity of reservoir and permeability are enhanced and then better reservoir rocks can be formed when they are combined with patchy dolostones. The relationship between porosity and the type of dolostones is that the dissolved neomorphic dolostones have the highest porosity of 3.65%, than followed by dissolved Mud-silt-sized dolostones of 3.35%. The mud-silt-sized dolostones without dissolution have the lowest porosity of 0.90%. Moreover, the porosity of prophyritic dolostones and the neomorphic dolostones without dissolution are lower, respectively 1.675% and 1.41%. Although algal laminated dolostones consist of euhedral crystals and riched in intercrystal pores, its porosity just yields 1.20%. The relationship between permeability and the type of dolostones is that that algal laminated dolostones have the highest permeability of 0.462mD and followed by 0.065mD of prophyritic dolostones. Dissolution have no significant influence on the permeability of neomorphic dolostones and this presented by the permeability of dissolved and non-dissolved are very close, respectively 0.043mD and 0.062mD. No matter dissolved or not, mud-silt-sized dolostones are much less permeable. The permeability of non-dissolved and dissolved are 0.051mD and 0.016mD. Collectively, in the study area, neomorphic dolostones can form high quality reservoir.

The Influence of Fractures on Vulnerability to Pollution of a Carbonate Aquifer: Monte Conero (Italy)

The carbonate outcrops of the anticline of Monte Conero (Italy) were studied in order to characterize the geometry of the fractures and to establish their influence on the petrophysical properties (hydraulic conductivity) and on the vulnerability to pollution. The outcrops form an analog for a fractured aquifer and belong to the Maiolica Fm. and the Scaglia Rossa Fm. The geometrical properties of fractures such as orientation, length, spacing and aperture were collected and statistically analyzed. Five types of mechanical fractures were observed: veins, joints, stylolites, breccias and faults. The types of fractures are arranged in different sets and geometric assemblages which form fracture networks. In addition, the fractures were analyzed at the microscale using thin sections. The fracture age-relationships resulted similar to those observed at the outcrop scale, indicating that at least three geological episodes have occurred in Monte Conero. A conceptual model for fault development was based on the observations of veins and stylolites. The fracture sets were modelled by the code FracSim3D to generate fracture network models. The permeability of a breccia zone was estimated at microscale by and point counting and binary image methods, whereas at the outcrop scale with Oda’s method. Microstructure analysis revealed that only faults and breccias are potential pathways for fluid flow since all veins observed are filled with calcite. According this, three scenarios were designed to asses the vulnerability to pollution of the analogue aquifer: the first scenario considers the Monte Conero without fractures, second scenario with all observed systematic fractures and the third scenario with open veins, joints and faults/breccias. The fractures influence the carbonate aquifer by increasing its porosity and hydraulic conductivity. The vulnerability to pollution depends also on the presence of karst zones, detric zones and the material of the vadose zone.

Ultrasonic and X-ray computed tomography characterization of progressive fracture damage in low-porous carbonate rocks

This paper studies the fracturing process in low-porous rocks during uniaxial compressive tests considering the original defects and the new mechanical cracks in the material. For this purpose, five different kinds of rocks have been chosen with carbonate mineralogy and low porosity (lower than 2%). The characterization of the fracture damage is carried out using three different techniques: ultrasounds, mercury porosimetry and X-ray computed tomography. The proposed methodology allows quantifying the evolution of the porous system as well as studying the location of new cracks in the rock samples. Intercrystalline porosity (the smallest pores with pore radius < 1 μm) shows a limited development during loading, disappearing rapidly from the porosimetry curves and it is directly related to the initial plastic behaviour in the stress–strain patterns. However, the biggest pores (corresponding to the cracks) suffer a continuous enlargement until the unstable propagation of fractures. The measured crack initiation stress varies between 0.25 σp and 0.50 σp for marbles and between 0.50 σp and 0.85 σp for micrite limestone. The unstable propagation of cracks is assumed to occur very close to the peak strength. Crack propagation through the sample is completely independent of pre-existing defects (porous bands, stylolites, fractures and veins). The ultrasonic response in the time-domain is less sensitive to the fracture damage than the frequency-domain. P-wave velocity increases during loading test until the beginning of the unstable crack propagation. This increase is higher for marbles (between 15% and 30% from initial vp values) and lower for micrite limestones (between 5% and 10%). When the mechanical cracks propagate unstably, the velocity stops to increase and decreases only when rock damage is very high. Frequency analysis of the ultrasonic signals shows clear changes during the loading process. The spectrum of treated waveforms shows two main frequency peaks centred at low (~ 20 kHz) and high (~ 35 kHz) values. When new fractures appear and grow the amplitude of the high-frequency peak decreases, while that of the low-frequency peak increases. Besides, a slight frequency shift is observed towards higher frequencies.

Sedimentology, palaeontology and diagenesis of the Much Wenlock limestone formation

Lithofacies distribution indicates that the Much Wenlock Limestone Formation of England and South Wales was desposited on a shelf which was flat and gently subsiding in the north, but topographically variable in the south. Limestone deposition in the north began with 12m of alga-rich limestone, which formed an upward shoaling sequence. Deepening then led to deposition of calcareous silty mudstones on the northern shelf. The remainder of the formation in this area formed during a shelf-wide regression, culminating in the production of an E to W younging sandbody. Lithofacies distribution on the southern shelf was primarily controlled by local subsidence. Six bedded lithofacies are recognised which contain 14 brachiopod/bryozoan dominated assemblages, of which 11 are in situ and three consist of reworked fossils. Microfacies analysis is necessary to distinguish assemblages which reflect original communities from those which reflect sedimentary processes. Turbulence, substrate-type, ease of feeding and other organisms in the environment controlled faunal distribution. Reefs were built dominantly by corals, stromatoporoids, algae and crinoids. Coral/stromatoporoid (Type A) reefs are common, particularly on the northern shelf, where they formed in response to shallowing, ultimately growing in front of the advancing carbonate sandbody. Algae dominate Type B and Type C reefs, reflecting growth in areas of poor water circulation. Lithification of the formation began in the marine-phreatic environment with precipitation of aragonite and high Mg calcite, which was subsequently altered to turbid low Mg calcite. Younger clear spars post-date secondary void formation. The pre-compactional clear spars have features which resemble the products of meteoric water diagenesis, but freshwater did not enter the formation at this time. The pre-compactional spars were precipitated by waters forced from the surrounding silty mudstones at shallow burial depths. Late diagenetic products are stylolites, compaction fractures and burial cements.

Criação de zonas de alta permeabilidade em fraturas: modelagem numérica do Lajedo de Soledade, Formação Jandaíra, Bacia Potiguar

The fracturing in carbonate rocks has been attracting increasingly attention due to new oil discoveries in carbonate reservoirs. This study investigates how the fractures (faults and joints) behave when subjected to different stress fields and how their behavior may be associated with the generation of karst and consequently to increased secondary porosity in these rocks. In this study I used satellite imagery and unmanned aerial vehicle UAV images and field data to identify and map faults and joints in a carbonate outcrop, which I consider a good analogue of carbonate reservoir. The outcrop comprises rocks of the Jandaíra Formation, Potiguar Basin. Field data were modeled using the TECTOS software, which uses finite element analysis for 2D fracture modeling. I identified three sets of fractures were identified: NS, EW and NW-SE. They correspond to faults that reactivate joint sets. The Ratio of Failure by Stress (RFS) represents stress concentration and how close the rock is to failure and reach the Mohr-Coulomb envelopment. The results indicate that the tectonic stresses are concentrated in preferred structural zones, which are ideal places for carbonate dissolution. Dissolution was observed along sedimentary bedding and fractures throughout the outcrop. However, I observed that the highest values of RFS occur in fracture intersections and terminations. These are site of karst concentration. I finally suggest that there is a relationship between stress concentration and location of karst dissolution in carbonate rocks.

Criação de zonas de alta permeabilidade em fraturas: modelagem numérica do Lajedo de Soledade, Formação Jandaíra, Bacia Potiguar

The fracturing in carbonate rocks has been attracting increasingly attention due to new oil discoveries in carbonate reservoirs. This study investigates how the fractures (faults and joints) behave when subjected to different stress fields and how their behavior may be associated with the generation of karst and consequently to increased secondary porosity in these rocks. In this study I used satellite imagery and unmanned aerial vehicle UAV images and field data to identify and map faults and joints in a carbonate outcrop, which I consider a good analogue of carbonate reservoir. The outcrop comprises rocks of the Jandaíra Formation, Potiguar Basin. Field data were modeled using the TECTOS software, which uses finite element analysis for 2D fracture modeling. I identified three sets of fractures were identified: NS, EW and NW-SE. They correspond to faults that reactivate joint sets. The Ratio of Failure by Stress (RFS) represents stress concentration and how close the rock is to failure and reach the Mohr-Coulomb envelopment. The results indicate that the tectonic stresses are concentrated in preferred structural zones, which are ideal places for carbonate dissolution. Dissolution was observed along sedimentary bedding and fractures throughout the outcrop. However, I observed that the highest values of RFS occur in fracture intersections and terminations. These are site of karst concentration. I finally suggest that there is a relationship between stress concentration and location of karst dissolution in carbonate rocks.