Mineral precipitation associated with vertical fault zones: the interaction of solute advection,diffusion and chemical kinetics

This article is concerned with chemical reactions that occur between two interacting parallel fluid flows using mixing in vertical faults as an example. Mineral precipitation associated with fluid flow in permeable fault zones results in mineralization and chemical reaction (alteration) patterns, which in turn are strongly dependent on interactions between solute advection (controlled by fluid flow rates), solute diffusion/dispersion and chemical kinetics. These interactions can be understood by simultaneously considering two dimensionless numbers, the Damköhler number and the Z‐number. The Damköhler number expresses the interaction between solute advection (flow rate) and chemical kinetics, while the Z‐number expresses the interaction between solute diffusion/dispersion and chemical kinetics. Based on the Damköhler and Z‐numbers, two chemical equilibrium length‐scales are defined, dominated by either solute advection or by solute diffusion/dispersion. For a permeable vertical fault zone and for a given solute diffusion/dispersion coefficient, there exist three possible types of chemical reaction patterns, depending on both the flow rate and the chemical reaction rate. These three types are: (i) those dominated by solute diffusion and dispersion resulting in precipitation at the lower tip of a vertical fault and as a thin sliver within the fault, (ii) those dominated by solute advection resulting in precipitation at or above the upper tip of the fault, and (iii) those in which advection and diffusion/dispersion play similar roles resulting in wide mineralization within the fault. Theoretical analysis indicates that there exists both an optimal flow rate and an optimal chemical reaction rate, such that chemical equilibrium following focusing and mixing of two fluids may be attained within the fault zone (i.e. type 3). However, for rapid and parallel flows, such as those resulting from a lithostatic pressure gradient, it is difficult for a chemical reaction to reach equilibrium within the fault zone, if the two fluids are not well mixed before entering the fault zone. Numerical examples are given to illustrate the three possible types of chemical reaction patterns.     

Following the Money: Corruption,Conflict, and the Winners and Losers of Suburban Land Acquisition in China

Land acquisition in suburban China is not merely the administrative means through which land is transferred from collectively owned rural land to state‐owned urban land; it is also a process of unlawful profit grabbing by rent‐seeking stakeholders. This licentiousness has been the source of intense controversy in recent years, becoming a major source of social unrest in contemporary China. Although standards of compensation are legally governed by strict policies, these standards are rarely followed as prescribed. As a result, we know very little about how revenue from land acquisition is shared in practice. This is because, as a contentious political issue rife with corrupt practices, accurate data are extremely difficult to obtain. This paper aims to address this vast empirical gap by tracing the flow of profit proliferation and distribution in a land acquisition case, using a suburban village in Shandong province, China, as an example. To expedite the drawn‐out acquisition process, the local government permits the corrupt behaviour of the village party secretary, who exploits his position to amass a great personal fortune. In contrast to previous research, which emphasises low compensation standards or excessive profit grabbing by local government, we find that villagers' protests are triggered by the corrupt practices of the village party secretary.     

Modelling the Lost City hydrothermal field: influence of topography and permeability structure

The Lost City hydrothermal field (LCHF) is hosted in serpentinite at the crest of the Atlantis Massif, an oceanic core complex close to the mid‐Atlantic Ridge. It is remarkable for its longevity and for venting low‐temperature (40–91°C) alkaline fluids rich in hydrogen and methane. IODP Hole U1309D, 5 km north of the LCHF, penetrated 1415 m of gabbroic rocks and contains a near‐conductive thermal gradient close to 100°C km^?1. This is remarkable so close to an active hydrothermal field. We present hydrothermal modelling using a topographic profile through the vent field and IODP site U1309. Long‐lived circulation with vent temperatures similar to the LCHF can be sustained at moderate permeabilities of 10^?14 to 10^?15 m² with a basal heatflow of 0.22 W m^?2. Seafloor topography is an important control, with vents tending to form and remain in higher topography. Models with a uniform permeability throughout the Massif cannot simultaneously maintain circulation at the LCHF and the near‐conductive gradient in the borehole, where permeabilities <10^?16 m² are required. A steeply dipping permeability discontinuity between the LCHF and the drill hole is required to stabilize venting at the summit of the massif by creating a lateral conductive boundary layer. The discontinuity needs to be close to the vent site, supporting previous inferences that high permeability is most likely produced by faulting related to the transform fault. Rapid increases in modelled fluid temperatures with depth beneath the vent agree with previous estimates of reaction temperature based on geochemical modelling.     

Direct inversion of Young's and Poisson impedances for fluid discrimination

Seismic inversion with prestack seismic data such as amplitude variation with offsets (AVO) inversion is an important tool in the estimation of elastic parameters for predicting lithology and discriminating fluid in conventional or unconventional hydrocarbon reservoirs. The product of Young's modulus and density (Young's impedance, YI) and the product of Poisson ratio and density (Poisson ratio impedance, PI) show great potential in lithology prediction and fluid discrimination of unconventional resources such as shale gas or oil. The high quality requirements for prestack data in density inversion render the estimation of YI and PI arduous and inaccurate with a conventional prestack inversion approach. In this study, a direct AVO inversion approach is proposed to estimate YI, PI, and density directly from P‐wave seismic data. The linearized P‐wave reflectivity approximate equation in terms of YI, PI, and density is initially derived. Five models, including four typical AVO classes, are utilized to verify the accuracy of the derived linearized P‐wave reflectivity equation in comparison with the exact P‐wave reflectivity equation and the frequently used linearized reflectivity approximate equation involving P‐ and S‐wave velocities and density. Parameter sensitivity analysis illustrates that YI and PI can reasonably be estimated from P‐wave reflectivity if a decorrelation scheme is utilized in the inversion algorithm. In addition, a pragmatic AVO inversion using a Bayesian scheme is suggested for the direct inversion of YI and PI from prestack seismic data. Synthetic and field data examples demonstrate the feasibility of the proposed inversion approach in the estimation of YI and PI and show the potential of this approach in fluid discrimination.     

Stress‐induced temperature variations in groundwater of the Monferrato area (north‐western Italy)

Although characterized by low seismicity, the Monferrato area of north‐western Italy was affected by earthquakes, of magnitude M5.1 and M4.8, in 2000 and 2001. At the same time, marked changes were recorded in water temperature and chemistry in several wells within the epicentral area. In May 2004, an automatic network for the continuous monitoring of groundwater was installed in selected wells to study the phenomenon. Here, we report on data collected during a 3‐year period of groundwater monitoring. During the first year, episodes of water heating (by up to 20°C) were observed in one monitored well. The temporal analysis of the seismic activity recorded in the area revealed as almost all seismic events occurred during the period of elevated water temperatures. The similar timing of earthquakes and groundwater‐temperature anomalies suggests that both may be triggered by the same processes acting in the crust.     

Mineralogy and fluid inclusion gas chemistry of production well mineral scale deposits at the Dixie Valley geothermal field,USA

At the Dixie Valley geothermal field, Nevada, USA, fluid boiling triggered the precipitation of carbonate scale minerals in concentric bands around tubing inserted into production well 28–33. When the tubing was removed, this mineral scale was sampled at 44 depth intervals between the wellhead and 1227 m depth. These samples provide a unique opportunity to evaluate the effects of fluid boiling on the scale mineralogy and geochemistry of the vapor and liquid phase. In this study, the mineralogy of the scale deposits and the composition of the fluid inclusion gases trapped in the mineral scales were analyzed. The scale consists mainly of calcite from 670–1112 m depth and aragonite from 1125 to 1227 m depth, with traces of quartz and Mg‐smectite. Mineral textures, including hopper growth, twinning, and fibrous growth in the aragonite and banded deposits of fine grained calcite crystals, are the result of progressive boiling. The fluid inclusion noncondensable gas was dominated by CO₂. However, significant variations in He relative to N₂ and Ar provide evidence that the geothermal reservoir consists of mixed source deeply circulating reservoir water and shallow, air saturated meteoric water. Gas analyses for many inclusions also showed higher CH₄ and H₂ relative to CO₂ than measured in gas sampled from this well, other production wells, and fumaroles. These inclusions are interpreted to have trapped CH₄‐ and H₂‐enriched gas resulting from early stages of boiling.     

Experimental study of the effect of variation in in‐situ stress on capillary residual trapping during CO2 geo‐sequestration in sandstone reservoirs

During a geo‐sequestration process, CO₂ injection causes an increase in reservoir pore pressure, which in turn decreases the reservoir net effective stress. Changes in effective stress can change all the reservoir and cap‐rock properties including residual saturations. This article presents the results of an experimental work carried out to understand the potential change in the volumes of residually trapped CO₂, while the porous medium tested underwent change in the net effective stress under in‐situ reservoir conditions of pore pressure and temperature. The experimental results obtained show that an initial 1725 psi (11.9 MPa) decrease in the net effective pressure caused 1.4% reduction in the volumes of residually trapped CO₂, while another 1500 psi (10.3 MPa) reduction caused a further 3.2% drop in the residual saturation of CO₂.     

Devolatilization‐induced pressure build‐up: Implications for reaction front movement and breccia pipe formation

Generation of fluids during metamorphism can significantly influence the fluid overpressure, and thus the fluid flow in metamorphic terrains. There is currently a large focus on developing numerical reactive transport models, and with it follows the need for analytical solutions to ensure correct numerical implementation. In this study, we derive both analytical and numerical solutions to reaction‐induced fluid overpressure, coupled to temperature and fluid flow out of the reacting front. All equations are derived from basic principles of conservation of mass, energy and momentum. We focus on contact metamorphism, where devolatilization reactions are particularly important owing to high thermal fluxes allowing large volumes of fluids to be rapidly generated. The analytical solutions reveal three key factors involved in the pressure build‐up: (i) The efficiency of the devolatilizing reaction front (pressure build‐up) relative to fluid flow (pressure relaxation), (ii) the reaction temperature relative to the available heat in the system and (iii) the feedback of overpressure on the reaction temperature as a function of the Clapeyron slope. Finally, we apply the model to two geological case scenarios. In the first case, we investigate the influence of fluid overpressure on the movement of the reaction front and show that it can slow down significantly and may even be terminated owing to increased effective reaction temperature. In the second case, the model is applied to constrain the conditions for fracturing and inferred breccia pipe formation in organic‐rich shales owing to methane generation in the contact aureole.     

Geodynamic processes in the NW Bohemian swarm earthquake region,Czech Republic,identified by continuous gas monitoring

The time series of two continuously operating gas monitoring stations at Oldřišská and Nový Kostel located along seismoactive faults in the epicentral area of the NW Bohemian swarm earthquakes (Czech Republic) are compared with water level fluctuations in two boreholes positioned along these faults and with gas flux variations of a mofette at the Soos mofette field at 9 km distance. The seasonal trend of the monitored CO₂ concentration with a maximum in November and a minimum in March/April is governed by groundwater temperatures, superimposed in spring by soil temperatures. CO₂ concentration variations identified at Oldřišská are also reflected in gas flux variations in the Soos mofette and/or water level fluctuations of two boreholes. Variations in the gas monitoring recordings of station at Nový Kostel are also linked with variations at Oldřišská. In all data sets, diurnal variations generated by earth tides occur, reflecting a daily stress – fault permeability cycle. Additional stress interferes with this cycle. Significant, abrupt changes are attributed to geodynamic processes linked with seismic events, as revealed by local seismicity or by the transient of waves of a strong remote earthquake. Simultaneous variations of the gas concentrations in the Nový Kostel area and in the gas flux in the Soos point to an interconnected hydraulic conductive fault systems present in the northern part of the Cheb Basin. Sharp falls in gas concentration, during or subsequent to, earthquake swarms may reflect fault compression associated with impeded gas migration. However, gas variations also occur in periods without seismic activity, indicating changes in fault permeability were caused by local aseismic fault movements, as revealed by events with opposite trends in the gas recordings at Oldřišská, Nový Kostel and the Soos. Therefore, a mathematical approach to establish a correlation between seismicity and gas geochemical variations is not possible.