Permeability evolution across carbonate hosted normal fault zones

Acknowledgements: The authors would like to thank Total E&P and BG Group for project funding and support, and the Industry Technology Facilitator for facilitating the collaborative development (grant number 3322PSD). The authors would also like to express their gratitude to the Aberdeen Formation Evaluation Society and the College of Physical Sciences at the University of Aberdeen for partial financial support. Raymi Castilla (Total E&P), Fabrizio Agosta and Cathy Hollis are also thanked for their constructive comments and suggestions to improve the standard of this manuscript as are John Still and Colin Taylor (University of Aberdeen) for technical assistance in the laboratory. Piero Gianolla is thanked for his editorial handling of the manuscript.

Mechanical behavior and localized failure modes in a porous basalt from the Azores

Basaltic rocks are the main component of the oceanic upper crust, thus of potential interest for water and geothermal resources, storage of CO2 and volcanic edifice stability. In this work, we investigated experimentally the mechanical behavior and the failure modes of a porous basalt, with an initial connected porosity of 18%. Results were acquired under triaxial compression experiments at confining pressure in the range of 25-200 MPa on water saturated samples. In addition, a purely hydrostatic test was also performed to reach the pore collapse critical pressure P*. During hydrostatic loading, our results show that the permeability is highly pressure dependent, which suggests that the permeability is mainly controlled by pre-existing cracks. When the sample is deformed at pressure higher than the pore collapse pressure P*, some very small dilatancy develops due to microcracking, and an increase in permeability is observed. Under triaxial loading, two modes of deformation can be highlighted. At low confining pressure (Pc < 50 MPa), the samples are brittle and shear localization occurs. For confining pressure > 50 MPa, the stress-strain curves are characterized by strain hardening and volumetric compaction. Stress drops are also observed, suggesting that compaction may be localized. The presence of compaction bands is confirmed by our microstructure analysis. In addition, the mechanical data allows us to plot the full yield surface for this porous basalt, which follows an elliptic cap as previously observed in high porosity sandstones and limestones.

Large-scale dynamic hydrofracturing, healing and fracture network characterization

Permeability of a rock is a dynamic property that varies spatially and temporally. Fractures provide the most efficient channels for fluid flow and thus directly contribute to the permeability of the system. Fractures usually form as a result of a combination of tectonic stresses, gravity (i.e. lithostatic pressure) and fluid pressures. High pressure gradients alone can cause fracturing, the process which is termed as hydrofracturing that can determine caprock (seal) stability or reservoir integrity. Fluids also transport mass and heat, and are responsible for the formation of veins by precipitating minerals within open fractures. Veining (healing) thus directly influences the rock’s permeability. Upon deformation these closed factures (veins) can refracture and the cycle starts again. This fracturing-healing-refacturing cycle is a fundamental part in studying the deformation dynamics and permeability evolution of rock systems. This is generally accompanied by fracture network characterization focusing on network topology that determines network connectivity. Fracture characterization allows to acquire quantitative and qualitative data on fractures and forms an important part of reservoir modeling. This thesis highlights the importance of fracture-healing and veins’ mechanical properties on the deformation dynamics. It shows that permeability varies spatially and temporally, and that healed systems (veined rocks) should not be treated as fractured systems (rocks without veins). Field observations also demonstrate the influence of contrasting mechanical properties, in addition to the complexities of vein microstructures that can form in low-porosity and permeability layered sequences. The thesis also presents graph theory as a characterization method to obtain statistical measures on evolving network connectivity. It also proposes what measures a good reservoir should have to exhibit potentially large permeability and robustness against healing. The results presented in the thesis can have applications for hydrocarbon and geothermal reservoir exploration, mining industry, underground waste disposal, CO2 injection or groundwater modeling.

Early profiles of clinical evolution after intravenous thrombolysis in an unselected stroke population.

BACKGROUND: Intravenous recombinant tissular plasminogen activator (rt-PA) is the only approved pharmacological treatment for acute ischaemic stroke. The authors aimed to analyse potential causes of the variable effect on early course and late outcome. METHODS AND RESULTS: 136 patients (42% women, 58% men) treated with intravenous rt-PA within 3 h of stroke onset in an acute stroke unit over a 3-year period, were included. Early clinical profiles of evolution at 48 h were divided into clinical improvement (CI) (decrease >4 points in the National Institute of Health Stroke Scale (NIHSS)); clinical worsening (CW) (increase >4 points NIHSS); clinical worsening after initial improvement (CWFI) (variations of >4 points in the NIHSS). Patients with clinical stability (no NIHSS modification or <4 points) were excluded. The patients showed in 66.9% CI, 13.2% CW 8.1 % CWFI and 11.8% remained stable. Female sex, no hyperlipaemia and peripheral arterial disease were associated with CW. Male sex and smoking were associated with CI. Absence of arterial occlusion on admission (28.4%) and arterial recanalisation at 24 h were associated with CI. Main causes of clinical deterioration included symptomatic intracranial haemorrhage (sICH), persistent occlusion and cerebral oedema. 23.5% developed ICH, 6.6% of which had sICH. At 3 months, 15.5% had died. Mortality was increased in CW, mainly related to sICH and cerebral oedema. The outcome of CWFI was intermediate between CW and CI. CONCLUSIONS: Early clinical profiles of evolution in thrombolysed patients vary considerably. Even with CI, it is critical to maintain vessel permeability to avoid subsequent CW.

A model of the three-dimensional evolution of Arctic melt ponds on first-year and multiyear sea ice

During winter the ocean surface in polar regions freezes over to form sea ice. In the summer the upper layers of sea ice and snow melts producing meltwater that accumulates in Arctic melt ponds on the surface of sea ice. An accurate estimate of the fraction of the sea ice surface covered in melt ponds is essential for a realistic estimate of the albedo for global climate models. We present a melt-pond–sea-ice model that simulates the three-dimensional evolution of melt ponds on an Arctic sea ice surface. The advancements of this model compared to previous models are the inclusion of snow topography; meltwater transport rates are calculated from hydraulic gradients and ice permeability; and the incorporation of a detailed one-dimensional, thermodynamic radiative balance. Results of model runs simulating first-year and multiyear sea ice are presented. Model results show good agreement with observations, with duration of pond coverage, pond area, and ice ablation comparing well for both the first-year ice and multiyear ice cases. We investigate the sensitivity of the melt pond cover to changes in ice topography, snow topography, and vertical ice permeability. Snow was found to have an important impact mainly at the start of the melt season, whereas initial ice topography strongly controlled pond size and pond fraction throughout the melt season. A reduction in ice permeability allowed surface flooding of relatively flat, first-year ice but had little impact on the pond coverage of rougher, multiyear ice. We discuss our results, including model shortcomings and areas of experimental uncertainty.

Modeling the summertime evolution of sea-ice melt ponds

1] We present a mathematical model describing the summer melting of sea ice. We simulate the evolution of melt ponds and determine area coverage and total surface ablation. The model predictions are tested for sensitivity to the melt rate of unponded ice, enhanced melt rate beneath the melt ponds, vertical seepage, and horizontal permeability. The model is initialized with surface topographies derived from laser altimetry corresponding to first-year sea ice and multiyear sea ice. We predict that there are large differences in the depth of melt ponds and the area of coverage between the two types of ice. We also find that the vertical seepage rate and the melt rate of unponded ice are important in determining the total surface ablation and area covered by melt ponds.

Prediction of water–rock interaction and porosity evolution in a granitoid-hosted enhanced geothermal system, using constraints from the 5 km Basel-1 well

Numerical simulations based on plans for a deep geothermal system in Basel, Switzerland are used here to understand chemical processes that occur in an initially dry granitoid reservoir during hydraulic stimulation and long-term water circulation to extract heat. An important question regarding the sustainability of such enhanced geothermal systems (EGS), is whether water–rock reactions will eventually lead to clogging of flow paths in the reservoir and thereby reduce or even completely block fluid throughput. A reactive transport model allows the main chemical reactions to be predicted and the resulting evolution of porosity to be tracked over the expected 30-year operational lifetime of the system. The simulations show that injection of surface water to stimulate fracture permeability in the monzogranite reservoir at 190 °C and 5000 m depth induces redox reactions between the oxidised surface water and the reduced wall rock. Although new calcite, chlorite, hematite and other minerals precipitate near the injection well, their volumes are low and more than compensated by those of the dissolving wall-rock minerals. Thus, during stimulation, reduction of injectivity by mineral precipitation is unlikely. During the simulated long-term operation of the system, the main mineral reactions are the hydration and albitization of plagioclase, the alteration of hornblende to an assemblage of smectites and chlorites and of primary K-feldspar to muscovite and microcline. Within a closed-system doublet, the composition of the circulated fluid changes only slightly during its repeated passage through the reservoir, as the wall rock essentially undergoes isochemical recrystallization. Even after 30 years of circulation, the calculations show that porosity is reduced by only ∼0.2%, well below the expected fracture porosity induced by stimulation. This result suggests that permeability reduction owing to water–rock interaction is unlikely to jeopardize the long-term operation of deep, granitoid-hosted EGS systems. A peculiarity at Basel is the presence of anhydrite as fracture coatings at ∼5000 m depth. Simulated exposure of the circulating fluid to anhydrite induces a stronger redox disequilibrium in the reservoir, driving dissolution of ferrous minerals and precipitation of ferric smectites, hematite and pyrite. However, even in this scenario the porosity reduction is at most 0.5%, a value which is unproblematic for sustainable fluid circulation through the reservoir.

Assessing Initial Conditions for Chloride Transport Across Low- permeability Argillaceous Rocks, Wellenberg, Switzerland

Information about fluid evolution and solute transport in a low-permeability metamorphic rock sequence has been obtained by comparing chloride concentrations and chlorine isotope ratios of pore water, groundwater, and fluid inclusions. The similarity of d37Cl values in fluid inclusions and groundwater suggests a closed-system evolution during the metamorphic overprint, and signatures established at this time appear to form the initial conditions for chloride transport after exhumation of the rock sequence.

Implementation and evaluation of permeability-porosity and tortuosity-porosity relationships linked to mineral dissolution-precipitation

Changes of porosity, permeability, and tortuosity due to physical and geochemical processes are of vital importance for a variety of hydrogeological systems, including passive treatment facilities for contaminated groundwater, engineered barrier systems (EBS), and host rocks for high-level nuclear waste (HLW) repositories. Due to the nonlinear nature and chemical complexity of the problem, in most cases, it is impossible to verify reactive transport codes analytically, and code intercomparisons are the most suitable method to assess code capabilities and model performance. This paper summarizes model intercomparisons for six hypothetical scenarios with generally increasing geochemical or physical complexity using the reactive transport codes CrunchFlow, HP1, MIN3P, PFlotran, and TOUGHREACT. Benchmark problems include the enhancement of porosity and permeability through mineral dissolution, as well as near complete clogging due to localized mineral precipitation, leading to reduction of permeability and tortuosity. Processes considered in the benchmark simulations are advective-dispersive transport in saturated media, kinetically controlled mineral dissolution-precipitation, and aqueous complexation. Porosity changes are induced by mineral dissolution-precipitation reactions, and the Carman-Kozeny relationship is used to describe changes in permeability as a function of porosity. Archie’s law is used to update the tortuosity and the pore diffusion coefficient as a function of porosity. Results demonstrate that, generally, good agreement is reached amongst the computer models despite significant differences in model formulations. Some differences are observed, in particular for the more complex scenarios involving clogging; however, these differences do not affect the interpretation of system behavior and evolution.

(Table T1) Porosity and permeability of samples from ODP Leg 180 sites

Woodlark Basin, an area of continental extension, is an ideal location to study the evolution of permeability and the development of overpressures within an active rift basin. In this investigation, we measured sediment permeabilities of cores from Woodlark Basin and used numerical modeling to determine if pore fluid overpressures are likely at the base of the rift basin. Constant-rate flow tests were conducted on cores from Site 1108, located in the rift basin, and Sites 1115 and 1118, located on the northern margin of the basin. Results of the laboratory tests indicated permeabilities that range from 1.5 x 10**-18 to 1 x 10**-16 m**2. Results of the numerical modeling of Site 1108 suggest that overpressures due to sedimentation are unlikely.

(Table T1) Density, porosity, void ratio, diffusivity and permeability of ODP Site 185-1149 sediments

The subduction of oceanic plates regulates crustal growth, influences arc volcanism, and refertilizes the mantle. Continental growth occurs by subduction of crustal material (seawater components, marine sediments, and basaltic crust). The geochemical and physical evolution of the Earth's crust depends, in large part, on the fate of subducted material at convergent margins (Armstrong, 1968, doi:10.1029/RG006i002p00175; Karig and Kay, 1981, 10.1098/rsta.1981.0108). The crustal material on the downgoing plate is recycled to various levels in the subduction zone. The recycling process that takes place in the "Subduction Factory" is difficult to observe directly but is clearly illuminated using chemical tracers. Von Huene and Scholl (1991, doi:10.1029/91RG00969) and Plank and Langmuir (1993, doi:10.1038/362739a0) preliminarily calculated a large flux of subducted materials. By mass balancing the chemical tracers and measuring the fractionations that occur between them, the Subduction Factory work and the effect on the Earth's evolution can be estimated. In order to elucidate this mass balance, Ocean Drilling Program Leg 185 drilled two deepwater shales into the oceanic crust situated in the Mariana-Izu Trenches and recovered core samples of incoming oceanic crust. The calculations of mass circulation in the subduction zone, however, did not take into account the mass transfer properties within subducted oceanic crust, although the dewatering fluid and diffused ions may play an important role in various activities such as seismogeneity, serpentine diapiring, and arc volcanism. Thus, this paper focuses on the quantitative measurements of the physical and mass transfer properties of subducted oceanic crust.

Evolution of strength and physical properties of ultramafic and carbonate rocks under hydrothermal conditions

Interaction of rocks with fluids can significantly change mineral assemblage and structure. This so-called hydrothermal alteration is ubiquitous in the Earth’s crust. Though the behavior of hydrothermally altered rocks can have planet-scale consequences, such as facilitating oceanic spreading along slow ridge segments and recycling volatiles into the mantle at subduction zones, the mechanisms involved in the hydrothermal alteration are often microscopic. Fluid-rock interactions take place where the fluid and rock meet. Fluid distribution, flux rate and reactive surface area control the efficiency and extent of hydrothermal alteration. Fluid-rock interactions, such as dissolution, precipitation and fluid mediated fracture and frictional sliding lead to changes in porosity and pore structure that feed back into the hydraulic and mechanical behavior of the bulk rock. Examining the nature of this highly coupled system involves coordinating observations of the mineralogy and structure of naturally altered rocks and laboratory investigation of the fine scale mechanisms of transformation under controlled conditions. In this study, I focus on fluid-rock interactions involving two common lithologies, carbonates and ultramafics, in order to elucidate the coupling between mechanical, hydraulic and chemical processes in these rocks. I perform constant strain-rate triaxial deformation and constant-stress creep tests on several suites of samples while monitoring the evolution of sample strain, permeability and physical properties. Subsequent microstructures are analyzed using optical and scanning electron microscopy. This work yields laboratory-based constraints on the extent and mechanisms of water weakening in carbonates and carbonation reactions in ultramafic rocks. I find that inundation with pore fluid thereby reducing permeability. This effect is sensitive to pore fluid saturation with respect to calcium carbonate. Fluid inundation weakens dunites as well. The addition of carbon dioxide to pore fluid enhances compaction and partial recovery of strength compared to pure water samples. Enhanced compaction in CO2-rich fluid samples is not accompanied by enhanced permeability reduction. Analysis of sample microstructures indicates that precipitation of carbonates along fracture surfaces is responsible for the partial restrengthening and channelized dissolution of olivine is responsible for permeability maintenance.

The Cratylia Mollis Seed Lectin Induces Membrane Permeability Transition In Isolated Rat Liver Mitochondria And A Cyclosporine A-insensitive Permeability Transition In Trypanosoma Cruzi Mitochondria.

Previous results provided evidence that Cratylia mollis seed lectin (Cramoll 1,4) promotes Trypanosoma cruzi epimastigotes death by necrosis via a mechanism involving plasma membrane permeabilization to Ca(2+) and mitochondrial dysfunction due to matrix Ca(2+) overload. In order to investigate the mechanism of Ca(2+) -induced mitochondrial impairment, experiments were performed analyzing the effects of this lectin on T. cruzi mitochondrial fraction and in isolated rat liver mitochondria (RLM), as a control. Confocal microscopy of T. cruzi whole cell revealed that Cramoll 1,4 binding to the plasma membrane glycoconjugates is followed by its internalization and binding to the mitochondrion. Electrical membrane potential (∆Ψm ) of T. cruzi mitochondrial fraction suspended in a reaction medium containing 10 μM Ca(2+) was significantly decreased by 50 μg/ml Cramoll 1,4 via a mechanism insensitive to cyclosporine A (CsA, membrane permeability transition (MPT) inhibitor), but sensitive to catalase or 125 mM glucose. In RLM suspended in a medium containing 10 μM Ca(2+) this lectin, at 50 μg/ml, induced increase in the rate of hydrogen peroxide release, mitochondrial swelling, and ∆Ψm disruption. All these mitochondrial alterations were sensitive to CsA, catalase, and EGTA. These results indicate that Cramoll 1, 4 leads to inner mitochondrial membrane permeabilization through Ca(2+) dependent mechanisms in both mitochondria. The sensitivity to CsA in RLM characterizes this lectin as a MPT inducer and the lack of CsA effect identifies a CsA-insensitive MPT in T. cruzi mitochondria.

Prognostic Value Of Dna And Mrna E6/e7 Of Human Papillomavirus In The Evolution Of Cervical Intraepithelial Neoplasia Grade 2.

This study aimed at evaluating whether human papillomavirus (HPV) groups and E6/E7 mRNA of HPV 16, 18, 31, 33, and 45 are prognostic of cervical intraepithelial neoplasia (CIN) 2 outcome in women with a cervical smear showing a low-grade squamous intraepithelial lesion (LSIL). This cohort study included women with biopsy-confirmed CIN 2 who were followed up for 12 months, with cervical smear and colposcopy performed every three months. Women with a negative or low-risk HPV status showed 100% CIN 2 regression. The CIN 2 regression rates at the 12-month follow-up were 69.4% for women with alpha-9 HPV versus 91.7% for other HPV species or HPV-negative status (P < 0.05). For women with HPV 16, the CIN 2 regression rate at the 12-month follow-up was 61.4% versus 89.5% for other HPV types or HPV-negative status (P < 0.05). The CIN 2 regression rate was 68.3% for women who tested positive for HPV E6/E7 mRNA versus 82.0% for the negative results, but this difference was not statistically significant. The expectant management for women with biopsy-confirmed CIN 2 and previous cytological tests showing LSIL exhibited a very high rate of spontaneous regression. HPV 16 is associated with a higher CIN 2 progression rate than other HPV infections. HPV E6/E7 mRNA is not a prognostic marker of the CIN 2 clinical outcome, although this analysis cannot be considered conclusive. Given the small sample size, this study could be considered a pilot for future larger studies on the role of predictive markers of CIN 2 evolution.

Effects Of Therapeutic Approach On The Neonatal Evolution Of Very Low Birth Weight Infants With Patent Ductus Arteriosus.

To analyze the effects of treatment approach on the outcomes of newborns (birth weight [BW] < 1,000 g) with patent ductus arteriosus (PDA), from the Brazilian Neonatal Research Network (BNRN) on: death, bronchopulmonary dysplasia (BPD), severe intraventricular hemorrhage (IVH III/IV), retinopathy of prematurity requiring surgical (ROPsur), necrotizing enterocolitis requiring surgery (NECsur), and death/BPD. This was a multicentric, cohort study, retrospective data collection, including newborns (BW < 1000 g) with gestational age (GA) < 33 weeks and echocardiographic diagnosis of PDA, from 16 neonatal units of the BNRN from January 1, 2010 to Dec 31, 2011. Newborns who died or were transferred until the third day of life, and those with presence of congenital malformation or infection were excluded. Groups: G1 - conservative approach (without treatment), G2 - pharmacologic (indomethacin or ibuprofen), G3 - surgical ligation (independent of previous treatment). Factors analyzed: antenatal corticosteroid, cesarean section, BW, GA, 5 min. Apgar score < 4, male gender, Score for Neonatal Acute Physiology Perinatal Extension (SNAPPE II), respiratory distress syndrome (RDS), late sepsis (LS), mechanical ventilation (MV), surfactant (< 2 h of life), and time of MV. death, O2 dependence at 36 weeks (BPD36wks), IVH III/IV, ROPsur, NECsur, and death/BPD36wks. Student's t-test, chi-squared test, or Fisher's exact test; Odds ratio (95% CI); logistic binary regression and backward stepwise multiple regression. Software: MedCalc (Medical Calculator) software, version 12.1.4.0. p-values < 0.05 were considered statistically significant. 1,097 newborns were selected and 494 newborns were included: G1 - 187 (37.8%), G2 - 205 (41.5%), and G3 - 102 (20.6%). The highest mortality was observed in G1 (51.3%) and the lowest in G3 (14.7%). The highest frequencies of BPD36wks (70.6%) and ROPsur were observed in G3 (23.5%). The lowest occurrence of death/BPD36wks occurred in G2 (58.0%). Pharmacological (OR 0.29; 95% CI: 0.14-0.62) and conservative (OR 0.34; 95% CI: 0.14-0.79) treatments were protective for the outcome death/BPD36wks. The conservative approach of PDA was associated to high mortality, the surgical approach to the occurrence of BPD36wks and ROPsur, and the pharmacological treatment was protective for the outcome death/BPD36wks.