Multiphase flow simulation through porous media with explicitly resolved fractures

Accurate simulation of multiphase flow in fractured porous media remains a challenge. An important problem is the representation of the discontinuous or near discontinuous behaviour of saturation in real geological formations. In the classical continuum approach, a refined mesh is required at the interface between fracture and porous media to capture the steep gradients in saturation and saturation‐dependent transport properties. This dramatically increases the computational load when large numbers of fractures are present in the numerical model. A discontinuous finite element method is reported here to model flow in fractured porous media. The governing multiphase porous media flow equations are solved in the adaptive mesh computational fluid dynamics code IC‐FERST on unstructured meshes. The method is based on a mixed control volume – discontinuous finite element formulation. This is combined with the P_N+1DG‐P_NDG element pair, which has discontinuous (order N+1) representation for velocity and discontinuous (order N) representation for pressure. A number of test cases are used to evaluate the method's ability to model fracture flow. The first is used to verify the performance of the element pair on structured and unstructured meshes of different resolution. Multiphase flow is then modelled in a range of idealised and simple fracture patterns. Solutions with sharp saturation fronts and computational economy in terms of mesh size are illustrated.     

Fluid mixing induced by hydrothermal activity in the ordovician carbonates in Tarim Basin,China

Permian hydrothermal activity in the Tarim Basin may have been responsible for the invasion of hot brines into Ordovician carbonate reservoirs. Studies have been undertaken to explain the origin and geochemical characteristics of the diagenetic fluid present during this hydrothermal event although there is no consensus on it. We present a genetic model resulting from the study of δ¹³C, δ¹⁸O, δ³⁴S, and ⁸⁷Sr/⁸⁶Sr isotope values and fluid inclusions (FIs) from fracture‐ and vug‐filling calcite, saddle dolomite, fluorite, barite, quartz, and anhydrite from Ordovician outcrops in northwest (NW) Tarim Basin and subsurface cores in Central Tarim Basin. The presence of hydrothermal fluid was confirmed by minerals with fluid inclusion homogenization temperatures being >10°C higher than the paleo‐formation burial temperatures both in the NW Tarim and in the Central Tarim areas. The mixing of hot (>200°C), high‐salinity (>24 wt% NaCl), ⁸⁷Sr‐rich (up to 0.7104) hydrothermal fluid with cool (60–100°C), low‐salinity (0 to 3.5 wt% NaCl), also ⁸⁷Sr‐rich (up to 0.7010) meteoric water in the Ordovician unit was supported by the salinity of fluid inclusions, and δ¹³C, δ¹⁸O, and ⁸⁷Sr/⁸⁶Sr isotopic values of the diagenetic minerals. Up‐migrated hydrothermal fluids from the deeper Cambrian strata may have contributed to the hot brine with high sulfate concentrations which promoted thermochemical sulfate reduction (TSR) in the Ordovician, resulting in the formation of ¹²C‐rich (δ¹³C as low as ?13.8‰) calcite and ³⁴S‐rich (δ³⁴S values from 21.4‰ to 29.7‰) H₂S, pyrite, and elemental sulfur. Hydrothermal fluid mixing with fresh water in Ordovician strata in Tarim Basin was facilitated by deep‐seated faults and up‐reaching faults due to the pervasive Permian magmatic activity. Collectively, fluid mixing, hydrothermal dolomitization, TSR, and faulting may have locally dissolved the host carbonates and increased the reservoir porosity and permeability, which has significant implications for hydrocarbon exploration.     

Polyphasic hydrothermal and meteoric fluid regimes during the growth of a segmented fault involving crystalline and carbonate rocks (Barcelona Plain,NE Spain)

A polyphasic tectonic‐fluid system of a fault that involves crystalline and carbonate rocks (Hospital fault, Barcelona Plain) has been inferred from regional to thin section scale observations combined with geochemical analyses. Cathodoluminescence, microprobe analyses and stable isotopy in fracture‐related cements record the circulation of successive alternations of hydrothermal and low‐temperature meteoric fluids linked with three main regional tectonic events. The first event corresponds to the Mesozoic extension, which had two rifting stages, and it is characterized by the independent tectonic activity of two fault segments, namely southern and northern Hospital fault segments. During the Late Permian‐Middle Jurassic rifting, these segments controlled the thickness and distribution of the Triassic sediments. Also, dolomitization was produced in an early stage by Triassic seawater at shallow conditions. During increasing burial, formation of fractures and their dolomite‐related cements took place. Fault activity during the Middle Jurassic–Late Cretaceous rifting was localized in the southern segment, and it was characterized by hydrothermal brines, with temperatures over 180°C, which ascended through this fault segment precipitating quartz, chlorite, and calcite. The second event corresponds to the Paleogene compression (Chattian), which produced exhumation, folding and erosion, favouring the percolation of low‐temperature meteoric fluids which produced the calcitization of the dolostones and of the dolomite cements. The third event is linked with the Neogene extension, where three stages have been identified. During the syn‐rift stage, the southern segment of the Hospital fault grew by tip propagation. In the relay zone, hydrothermal brines with temperature around 140°C upflowed. During the late postrift, the Hospital fault acted as a unique segment and deformation occurred at shallow conditions and under a low‐temperature meteoric regime. Finally, and possibly during the Messinian compression, NW‐SE strike‐slip faults offset the Hospital fault to its current configuration.     

内乡县衙建置沿革与现存遗迹考

内乡县衙所存遗迹是清末章炳焘时代留下的。清末县衙的东西两路与明代及清光绪重建前大不相同。今天按照明清两代地方志复原的建筑群,不仅与章炳焘所留的清末旧衙署基本没有什么关系,与历史记载中的明内乡县衙也相去甚远。鉴于该衙正成为史学界研究明清官僚政体的新史料,弄清楚它的沿革、遗迹和现状间的关系是必要的。中国衙署建筑所存遗迹极稀贵,如内乡县衙可反映行政制度变迁者更少,故应特别慎重对待。     

Heat transport at boiling,near‐critical conditions

Two‐phase flow and near‐critical phenomena are likely to enhance energy transport in high‐temperature hydrothermal systems. We present a series of two‐dimensional simulations of two‐phase flow of pure water at near‐critical conditions. The results show that at near‐critical conditions, two‐phase convection can be more efficient in transporting energy than single‐phase convection. The highest heat fluxes are attained when two‐phase heat‐pipes form near the bottom boundary, recharging the root of the upflow zone and thereby enabling the formation of broad upflow regions. When the system becomes more vapor‐dominated, it loses this ability, upflow zones become narrower and the energy efficiency drops to more moderate values.     

Long‐term chemical evolution and modification of continental basement brines – a field study from the Schwarzwald,SW Germany

Highly saline, deep‐seated basement brines are of major importance for ore‐forming processes, but their genesis is controversial. Based on studies of fluid inclusions from hydrothermal veins of various ages, we reconstruct the temporal evolution of continental basement fluids from the Variscan Schwarzwald (Germany). During the Carboniferous (vein type i), quartz–tourmaline veins precipitated from low‐salinity (<4.5wt% NaCl + CaCl₂), high‐temperature (≤390°C) H₂O‐NaCl‐(CO₂‐CH₄) fluids with Cl/Br mass ratios = 50–146. In the Permian (vein type ii), cooling of H₂O‐NaCl‐(KCl‐CaCl₂) metamorphic fluids (T ≤ 310°C, 2–4.5wt% NaCl + CaCl₂, Cl/Br mass ratios = 90) leads to the precipitation of quartz‐Sb‐Au veins. Around the Triassic–Jurassic boundary (vein type iii), quartz–haematite veins formed from two distinct fluids: a low‐salinity fluid (similar to (ii)) and a high‐salinity fluid (T = 100–320°C, >20wt% NaCl + CaCl₂, Cl/Br mass ratios = 60–110). Both fluids types were present during vein formation but did not mix with each other (because of hydrogeological reasons). Jurassic–Cretaceous veins (vein type iv) record fluid mixing between an older bittern brine (Cl/Br mass ratios ~80) and a younger halite dissolution brine (Cl/Br mass ratios >1000) of similar salinity, resulting in a mixed H₂O‐NaCl‐CaCl₂ brine (50–140°C, 23–26wt% NaCl + CaCl₂, Cl/Br mass ratios = 80–520). During post‐Cretaceous times (vein type v), the opening of the Upper Rhine Graben and the concomitant juxtaposition of various aquifers, which enabled mixing of high‐ and low‐salinity fluids and resulted in vein formation (multicomponent fluid H₂O‐NaCl‐CaCl₂‐(SO₄‐HCO₃), 70–190°C, 5–25wt% NaCl‐CaCl₂ and Cl/Br mass ratios = 2–140). The first occurrence of highly saline brines is recorded in veins that formed shortly after deposition of halite in the Muschelkalk Ocean above the basement, suggesting an external source of the brine's salinity. Hence, today's brines in the European basement probably developed from inherited evaporitic bittern brines. These were afterwards extensively modified by fluid–rock interaction on their migration paths through the crystalline basement and later by mixing with younger meteoric fluids and halite dissolution brines.     

Hydrothermal activity associated with the Ries impact event, Germany

Combined field studies, optical and scanning electron microscopy, and electron microprobe studies of impactites from the Ries impact structure, Germany, have allowed a clearer picture of the hydrothermal system associated with the Ries impact event to be made. Hydrothermal alteration is concentrated within impact‐generated suevites in the interior of the crater (crater suevites) and around the periphery (surficial suevites), with minor alteration in the overlying sedimentary crater‐fill deposits. The major heat source for the Ries hydrothermal system was the suevite units themselves. Hydrothermal alteration of crater‐fill suevites is pervasive in nature and comprises several distinct alteration phases that vary with depth. An early phase of K‐metasomatism accompanied by minor albitization of crystalline basement clasts and minor chloritization, was followed by pervasive intermediate argillic alteration (predominantly montmorillonite, saponite, and illite) and zeolitization (predominantly analcite, erionite, and clinoptilolite). Hydrothermal fluids were typically weakly alkaline during the main stage of alteration. In contrast to the crater‐fill suevites, alteration within surficial suevites was typically restricted to montmorillonite and phillipsite deposition within cavities and fractures. The pervasive nature of the alteration within the crater‐fill suevites was likely due to the presence of an overlying crater lake; whereas alteration within surficial suevites typically occurred under undersaturated conditions with the main source of water being from precipitation. There are exceptional outcrops of more pervasively altered surficial suevites, which can be explained as locations where water pooled for longer periods of time. Hydrothermal fluids were likely a combination of meteoric waters that percolated down from the overlying crater lake and groundwaters that flowed in from the surrounding country rocks.     

建国初期中国共产党对私营金融业的社会主义改造

建国初期,我党参照、借鉴前苏联的经验,根据本国实际情况,通过加强对私营金融业的管理和整顿、积极引导私营金融业的资金投放和业务经营、组织私营行庄联合经营和联合管理、进行全行业社会主义改造四个阶段,最终完成了对私营金融业的社会主义改造任务,建立了社会主义的金融体系.     

The potential for abiotic organic synthesis and biosynthesis at seafloor hydrothermal systems

Calculations are presented of the extent to which chemical disequilibria are generated when submarine hydrothermal fluids mix with sea water. These calculations involve quantifying the chemical affinity for individual reactions by comparing equilibrium compositions with the compositions of mixtures in which oxidation–reduction reactions are inhibited. The oxidation–reduction reactions that depart from equilibrium in these systems provide energy for chemotrophic microbial metabolism. Methanogenesis is an example of this phenomenon, in which the combination of carbon dioxide, hydrogen and methane induced by fluid mixing is far from equilibrium, which can be approached if more methane is generated. Similarly, the production of other organic compounds is also favorable under the same conditions that permit methanogenesis. Alkanes, alkenes, alcohols, aldehydes, carboxylic acids and amino acids are among the compounds that, if formed, would lower the energetic state of the chemical composition generated in mixed fluids. The resulting positive values of chemical affinity correspond to the thermodynamic drive required for abiotic organic synthesis. It is also possible that energy release accompanies biosynthesis by chemotrophic organisms. In this way, hydrothermal ecosystems differ radically from familiar ecosystems at Earth's surface. If captured, the energy released may be sufficient to drive biosynthesis of carbohydrates, purines, pyrimidines and other compounds which require energy inputs.     

Chemical modification of pyroclastic rock by hot water: an experimental investigation of mass transport at the fluid–solid interface

Hydrothermal water–(pyroclastic) rock interactions were examined using flow-through experiments to deduce the effect of mass transport phenomena on the reaction process. A series of experiments were conducted over the temperature range 75–250°C, with a constant temperature for each experiment, and at saturated vapour pressure, to estimate the apparent rate constants as a function of temperature.
Based on the chemistry of analysed solutions, the water–rock interaction in the experiments was controlled by diffusion from the reaction surface and by the existence of a surface layer at the rock–fluid interface, which regulated the chemical reaction rate. The reaction progress depended to a high degree on flow velocity and temperature conditions, with element abundances in the fluid significantly affected by these factors. Mass transport coefficients for diffusion from the rock surface to the bulk solution have been estimated. Ca is selectively depleted under lower temperature conditions ( T  < 150°C), whereas Na is greatly depleted under higher temperature conditions ( T  > 150°C), and K reaction rates are increased when flow velocity increases. Using these conditions, specific alkali and alkali earth cations were selectively leached from mineral surfaces. The 'surface layer' comprised a 0.5–1.8 mm boundary film on the solution side (the thickness of this layer has no dependence on chemical character) and a reaction layer. The reaction layer was composed of a Si, Al-rich cation-leached layer, whose thickness was dependent on temperature, flow velocity and reaction length. The reaction layer varied in thickness from about 10⁻⁴ to 10⁻⁷ mm under high temperature/low fluid velocity and low temperature/high fluid velocity conditions, respectively.