共查询到20条相似文献，搜索用时 15 毫秒
1.
Production from geothermal reservoirs in volcanic-igneous hydrothermal systems may be disturbed owing to the formation of
authigenic minerals that reduce primary porosity and infill fissure systems. Crystallization may be induced by natural processes
or human activity related to the reservoir exploitation. In volcanic and volcaniclastic rocks, zeolites and related authigenic
silicate minerals commonly develop. Two selected study sites—Pauzhetka geothermal field in Russia and Smrekovec volcanic complex
in Slovenia are characterized by extensive development of progressive stage medium- and low-temperature propylitization and
zeolitization with similar mineral assemblages that reduced original reservoir rock porosity and permeability. Retrograde
and overprinting reactions of argillitization commonly enhance porosity and permeability of rocks, but some other reactions,
like from prehnite to laumontite, from laumontite to heulandite, and from laumontite to analcime, encountered in the Smrekovec
volcanic complex reduced it appreciably. Retrograde reactions recognized in the study sites were not induced by the temperature
drop only, but were accompanied by the change in chemical composition of reacting fluids. Chemical composition of interstitial
waters produced from exploitation wells in the Pauzhetka geothermal field indicate that propylitic and zeolite facies alteration
is related to Na–Cl, slightly alkaline waters, while argillitization involved acidic thermal waters with more complex ion
composition. 相似文献
2.
Nathan Benkemoun Hayder Al Khazraji Philippe Poullain Marta Choinska Abdelhafid Khelidj 《国际地质力学数值与分析法杂志》2018,42(3):449-468
In this paper, mesoscale hydromechanical simulations are performed to study (1) fracture features and (2) crack‐gas permeability coupling evolution in the context of the tensile splitting test. The mesostructure is based on a 2‐phase 3‐D representation of heterogeneous materials, such as concrete, where stiff aggregates are embedded into a mortar matrix. To take into account these heterogeneities without any mesh adaptation, a weak discontinuity is introduced into the strain field. In addition, a strong discontinuity is also added to take into account microcracking. This mechanical model is cast into the framework of the enhanced finite element method. Concerning the coupling with gas permeability, a double‐porosity method is used to simulate the flow through the cracks and the porosity. The apparent gas permeability is afterwards evaluated by a homogenization method. On the basis of finite element simulations, influence of aggregate size on ultimate crack opening, macroscopic ultimate tensile stress, total dissipated energy, and gas permeability evolution is numerically investigated. Furthermore, gas permeability evolution is also compared with experimental results from the literature. In addition, in the spirit of a sequential multiscale approach, macroscale gas permeability equations are identified from the hydromechanical results coming from the mesoscale computations. These equations lead to a relation between macroscale gas permeability evolution and crack opening. Besides, we show how the aggregate size influences the percolation threshold and that after this threshold, a cubic relation between macroscale gas permeability and crack opening is obtained. 相似文献
3.
构造-流体-成矿体系的反应-输运-力学耦合模型和动力学模拟 总被引：12，自引：0，他引：12
建立了一个综合的构造流体成矿体系的反应输运力学耦合动力学模型。利用有限元方法求解岩石变形、断裂作用和断裂网络统计动力学、流体流动、有机和无机地球化学反应及成岩成矿作用、压力溶液和其它压实力学、热迁移的方程组 ,可以对构造流体成矿体系的动力学演化过程进行 1～ 3维数值模拟。模拟的主要内容是在各种过程耦合作用下描述构造流体成矿体系的主要变量的时空演化 :( 1)与成矿流体的形成和性质有关的变量 ,如地层中矿物 (包括成矿物质 )的溶解速率、流体中各组分的浓度与饱和度、流体温度、压力、离子强度等 ;( 2 )与构造变形和流体运移有关的各变量 ,如应力与变形速率、岩石孔隙度、构造 (断裂 )渗透率等 ;( 3 )与沉淀成矿有关的变量 ,如矿物 (金属矿物和脉石矿物 )的成核速率、各矿物的沉淀量等 ;( 4 )上述各有关变量间的时空耦合关系 ,如断裂渗透率时空演化与流体流动、汇聚和成矿的耦合关系等。以湖南沃溪金锑钨矿床为例 ,应用该模型和方法对成矿动力学过程和动力学机制进行了初步的模拟与分析。 相似文献
4.
Gour-Tsyh Yeh Yilin Fang Fan Zhang Jiangtao Sun Yuan Li Ming-Hsu Li Malcolm D. Siegel 《Computational Geosciences》2010,14(1):149-170
Subsurface contamination problems of metals and radionuclides are ubiquitous. Metals and radionuclides may exist in the solute
phase or may be bound to soil particles and interstitial portions of the geologic matrix. Accurate tools to reliably predict
the migration and transformation of these metals and radionuclides in the subsurface environment enhance the ability of environmental
scientists, engineers, and decision makers to analyze their impact and to evaluate the efficacy of alternative remediation
techniques prior to incurring expense in the field. A mechanistic-based numerical model could provide such a tool. This paper
communicates the development and verification of a mechanistically coupled fluid-flow thermal-reactive biogeochemical-transport
model where both fast and slow reactions occur in porous and fractured media. Theoretical bases, numerical implementations,
and numerical experiments using the model are described. A definition of the “rates” of fast/equilibrium reactions is presented
to come up with a consistent set of governing equations. Two example problems are presented. The first one is a reactive transport
problem which elucidates the non-isothermal effects on heterogeneous reactions. It also demonstrates that the rates of fast/equilibrium
reactions are not necessarily greater than that of slow/kinetic reactions in the context of reactive transport. The second
example focuses on a complicated but realistic advective–dispersive–reactive transport problem. This example exemplifies the
need for innovative numerical algorithms to solve problems involving stiff geochemical reactions. It also demonstrates that
rates of all fast/equilibrium reactions are finite and definite. Furthermore, it is noted that a species-versus-time curve
cannot be used to characterize the rate of homogeneous fast/equilibrium reaction in a reactive transport system even if one
and only one such reaction is responsible for the production of this species. 相似文献
5.
孔隙结构是制约页岩气赋存与运移的重要因素。页岩致密、低渗,其孔隙度极低。酸液可溶蚀页岩中的部分矿物,使其孔隙增加,促进页岩内部的气体流动。论文利用稀盐酸、氢氟酸、混合酸分别对页岩试样进行了不同时间的酸化处理,应用扫描电镜(SEM)观测了酸化前后湖南牛蹄塘组页岩试样表面的孔隙。基于MATLAB对页岩酸化前后的电镜扫描结果进行了孔隙度分析,页岩试样孔隙度为10.68%,纯水、稀盐酸、氢氟酸、混合酸分别浸泡24 h后,试样孔隙度分别增加了1.21倍、1.3倍、2.61倍、2.99倍。利用扫描电镜能谱仪测试了酸化前后页岩的元素变化,分析了页岩试样与不同酸液反应的机理和反应程度。结果表明:1 h酸化后,氢氟酸作用后的页岩试样表面残有大量石英和黏土矿物;而混合酸作用后的页岩试样表面虽残有大量石英,但黏土矿物含量明显减小,混合酸液作用下黏土矿物更容易被溶蚀。混合酸溶液酸性最强,24 h酸化后,混合酸与石英和黏土矿物剧烈反应。混合酸作用后页岩的孔隙最为发育,促进了页岩渗透率的增加。 相似文献
6.
Hydrothermal alteration of a quartz‐K‐feldspar rock is simulated numerically by coupling fluid flow and chemical reactions. Introduction of CO2 gas generates an acidic fluid and produces secondary quartz, muscovite and/or pyrophyllite at constant temperature and pressure of 300°C and 200 MPa. The precipitation and/or dissolution of the secondary minerals is controlled by either mass‐action relations or rate laws. In our simulations the mass of the primary elements are conserved and the mass‐balance equations are solved sequentially using an implicit scheme in a finite‐element code. The pore‐fluid velocity is assumed to be constant. The change of rock volume due to the dissolution or precipitation of the minerals, which is directly related to their molar volume, is taken into account. Feedback into the rock porosity and the reaction rates is included in the model. The model produces zones of pyrophyllite quartz and muscovite due to the dissolution of K‐feldspar. Our model simulates, in a simplified way, the acid‐induced alteration assemblages observed in various guises in many significant mineral deposits. The particular aluminosilicate minerals produced in these experiments are associated with the gold deposits of the Witwatersrand Basin. 相似文献
7.
Alexis Navarre-Sitchler Carl I. Steefel Susan L. Brantley 《Geochimica et cosmochimica acta》2011,75(23):7644-7667
Saprolite formation rates influence many important geological and environmental issues ranging from agricultural productivity to landscape evolution. Here we investigate the chemical and physical transformations that occur during weathering by studying small-scale “saprolites” in the form of weathering rinds, which form on rock in soil or saprolite and grow in thickness without physical disturbance with time. We compare detailed observations of weathered basalt clasts from a chronosequence of alluvial terraces in Costa Rica to diffusion-reaction simulations of rind formation using the fully coupled reactive transport model CrunchFlow. The four characteristic features of the weathered basalts which were specifically used as criteria for model comparisons include (1) the mineralogy of weathering products, (2) weathering rind thickness, (3) the coincidence of plagioclase and augite reaction fronts, and (4) the thickness of the zones of mineral reaction, i.e. reaction fronts. Four model scenarios were completed with varying levels of complexity and degrees of success in matching the observations. To fit the model to all four criteria, however, it was necessary to (1) treat diffusivity using a threshold in which it increased once porosity exceeded a critical value of 9%, and (2) treat mineral surface area as a fitting factor. This latter approach was presumably necessary because the mineral-water surface area of the connected (accessible) porosity in the Costa Rica samples is much less than the total porosity (Navarre-Sitchler et al., 2009). The model-fit surface area, here termed reacting surface area, was much smaller than the BET-measured surface area determined for powdered basaltic material. In the parent basalt, reacting surface area and diffusivity are low due to low pore connectivity, and early weathering is therefore transport controlled. However, as pore connectivity increases as a result of weathering, the reacting surface area and diffusivity also increase and weathering becomes controlled by mineral reaction kinetics. The transition point between transport and kinetic control appears to be related to a critical porosity (9%) at which pore connectivity is high enough to allow rapid transport. Based on these simulations, we argue that the rate of weathering front advance is controlled by the rate at which porosity is created in the weathering interface, and that this porosity increases because of mineral dissolution following a rate that is largely surface-reaction controlled. 相似文献
8.
In CO2 geological storage (CGS) context, the evolution of the caprock sealing capacity has received increasing attention, particularly on a geological time span (thousands of years). At this time scale, geochemical reactions may enhance or weaken the caprock quality. It is widely recognized that, for the reservoir, geological heterogeneities affect the concentration and spatial distribution of CO2, and then affect the extent of gas–water–rock interactions, which in turn alters the hydrogeological properties of the reservoir. However, much less attention of these effects has been paid to the caprock. In this study, we presented and applied a novel approach to evaluate the effects of permeability and porosity heterogeneities on the alteration of minerals, the associated evolution of the caprock sealing efficiency and the containment of supercritical CO2 (scCO2) within the caprock. Even though this is a generic study, several conditions and parameters such as pressure, permeability, and mineral composition, were extracted from a caprock layer of the Shiqianfeng Formation in the Ordos Basin demonstration site in China. For the sake of simplification, a 2-dimensional model was designed to represent the caprock domain. We firstly generated an appropriate heterogeneous random field of permeability with the average permeability taken from the uppermost mudstone layer of the Shiqianfeng Formation, and then the heterogeneity in porosity was incorporated using a joint normal distribution method based on the available data. Homogeneous mineral compositions of the reservoir and caprock were used in all simulations. Simulations of three cases were performed, including a homogeneous case, a case with only permeability heterogeneity and a case with both permeability and porosity heterogeneities. The results demonstrate dramatic influences of permeability and porosity heterogeneities on the migration of scCO2 within the caprock, the alteration of minerals, and therefore the evolution of the caprock sealing quality. Specific to the data used in this study, hydrogeological heterogeneities facilitated the overall penetration of scCO2 within the caprock and promoted the alteration of minerals, thereby weakening the caprock sealing efficiency over the simulation time. 相似文献
9.
10.
结合化学反应方程式,并应用多孔多相介质溶混污染物输运过程的数值模型,对多孔多相介质中含均相/非均相化学反应传质过程进行了数值模拟。化学反应主要包含均相快速/慢速和非均相快速/慢速等5种化学反应过程,溶质输运行为的控制机制主要考虑对流、扩散及降解、吸附等。基于原有的隐式特征线Galerkin离散化的有限元方法,求解模型控制方程的边值初值问题,求解过程中把均相化学反应物质中按照反应物和生成物分开,非均相反应物质按照固相和液相分开,对均相反应物及非均相液相物质浓度耦合求解,而均相生成物和非均相固相物质独立求解。使方程组按照其不同类型进行分类,同时可减少未知数的个数。对于含有非线性内状态变量的右端项进行迭代求解。数值例题结果验证了所提出的数值方法的有效性、计算精度和稳定性。 相似文献
11.
12.
13.
《Applied Geochemistry》2003,18(10):1555-1571
In the context of the proposed low- and intermediate-level radioactive waste repository at Wellenberg (Switzerland), calculations simulating the interaction between hyperalkaline solutions and a fractured marl, at 25 °C, have been performed. The aim of these calculations is to evaluate the possible effects of mineral dissolution and precipitation on porosity and permeability changes in such a fractured marl, and their impact on repository performance. Solute transport and chemical reaction are considered in both a high-permeability zone (fracture), where advection is important, and the wall rock, where diffusion is the dominant transport mechanism. The mineral reactions are promoted by the interaction between hyperalkaline solutions derived from the degradation of cement (a major component of the engineered barrier system in the repository) and the host rock. Both diffusive/dispersive and advective solute transport are taken into account in the calculations. Mineral reactions are described by kinetic rate laws. The fluid flow system under consideration is a two-dimensional porous medium (marl, 1% porosity), with a high-permeability zone simulating a fracture (10% porosity) crossing the domain. The dimensions of the domain are 6 m per 1 m, and the fracture width is 10 cm. The fluid flow field is updated during the course of the simulations. Permeabilities are updated according to Kozeny's equation. The composition of the solutions entering the domain is derived from modeling studies of the degradation of cement under the conditions at the proposed underground repository at Wellenberg. Two different cases have been considered in the calculations. These 2 cases are representative of 2 different stages in the process of degradation of cement (pH 13.5 and pH 12.5). In both cases, the flow velocity in the fracture diminishes with time, due to a decrease in porosity. This decrease in porosity is caused by the precipitation of calcite (replacement of dolomite by calcite) and other secondary minerals (brucite, sepiolite, analcime, natrolite, tobermorite). However, the decrease in porosity and flow velocity is much more pronounced in the lower pH case. The extent of the zone of mineral alteration along the fracture is also much more limited in the lower pH case. The reduction of porosity in the fractures would be highly beneficial for repository performance, since it would mean that the solutions coming from the repository and potentially carrying radionuclides in solution would have to flow through low-conductive rock before they would be able to get to higher-conductive features. The biggest uncertainty in the reaction rates used in the calculations arises from the surface areas of the primary minerals. Additional calculations making use of smaller surface areas have also been performed. The results show that the smaller surface areas (and therefore smaller reaction rates) result in a smaller reactivity of the system and smaller porosity changes. 相似文献
14.
XUE Quan ZHANG Liwei XU Liang Matthew MYERS Cameron WHITE MEI Kaiyuan WANG Hanwen LI Qi LI Xiaochun 《《地质学报》英文版》2023,97(3):946-955
In order to improve CO2 capture, utilization and storage (CCUS) to solve carbon emission, sandstone from the Triassic Liujiagou Formation (LF) from the Ordos Basin in China was investigated using permeability tests and computed X-ray tomography (CT) scanning. The presence of reactive minerals within the geological CO2 sequestration target storage formation can allow reaction with injected CO2, which changes the porosity and permeability of the LF beds, affecting storage effectiveness. To investigate the effect of chemical reactions on the pore structure and permeability of sandstone cores representing the LF CO2 storage, tests were conducted to analyze the changes in porosity and permeability of sandstone cores induced by CO2-saturated brine at different reaction times (28-day maximum reaction period). Porosity and permeability of the sandstone increased after reaction with CO2-saturated brine due to mineral dissolution. The sandstone exhibited an increase in porosity and permeability after 15 days of reaction with CO2-saturated brine. Moreover, there was an increase in the volume of large pores in the sandstone after the 28-day period. The pore network of the sandstone was established through CT results, and the porosity calculated based on the obtained pore network was close to that measured in the test, demonstrating the feasibility to use CT to study the evolution of the microstructure of sandstone after long-time exposure to CO2-saturated brine. 相似文献
15.
A numerical model of coupled saturated/unsaturated water flow, heat transfer and multi-component reactive solute transport
is presented to evaluate the long-term geochemical evolution in bentonite, concrete and clay formation for a potential geological
radioactive waste repository. Changes in formation porosity caused by mineral dissolution/precipitation reactions are taken
into account. Simulations were carried out with a general-purpose multicomponent reactive transport code, CORE2D V4. Numerical results show that pH in the bentonite porewater can vary from neutral to up to 13 over a time scale of 1 Ma
although dissolution of silica minerals and precipitation of secondary calcium silicate hydrate (CSH) minerals in bentonite
buffer the effect of the hyperalkaline plume. Mineral precipitation reduces the volume of pore space in bentonite close to
the bentonite–concrete interface due to the precipitation of CSH minerals. Model results indicate that bentonite porosity
decreases less than 25%. The hyperalkaline plume from the concrete only extends to a distance of 0.7 m in the clay formation
over the time range of 1 Ma. 相似文献
16.
The transport of chemically reactive solutes (e.g. surfactants, CO2 or dissolved minerals) is of fundamental importance to a wide range of applications in oil and gas reservoirs such as enhanced
oil recovery and mineral scale formation. In this work, we investigate exponential time integrators, in conjunction with an
upwind weighted finite volume discretisation in space, for the efficient and accurate simulation of advection–dispersion processes
including non-linear chemical reactions in highly heterogeneous 3D oil reservoirs. We model sub-grid fluctuations in transport
velocities and uncertainty in the reaction term by writing the advection–dispersion–reaction equation as a stochastic partial
differential equation with multiplicative noise. The exponential integrators are based on the variation of constants solution
and solve the linear system exactly. While this is at the expense of computing the exponential of the stiff matrix representing
the finite volume discretisation, the use of real Léja point or the Krylov subspace technique to approximate the exponential
makes these methods competitive compared to standard finite difference-based time integrators. For the deterministic system,
we investigate two exponential time integrators, the second-order accurate exponential Euler midpoint (EEM) scheme and exponential
time differencing of order one (ETD1). All our numerical examples demonstrate that our methods can compete in terms of efficiency
and accuracy compared with standard first-order semi-implicit time integrators when solving (stochastic) partial differential
equations that model mixing and chemical reactions in 3D heterogeneous porous media. Our results suggest that exponential
time integrators such as the ETD1 and EEM schemes could be applied to typical 3D reservoir models comprising tens to hundreds
of thousands unknowns. 相似文献
17.
We use a reactive diffusion model to investigate what happens to CO2 injected into a subsurface sandstone reservoir capped by a chlorite- and illite-containing shale seal. The calculations simulate reaction and transport of supercritical (SC) CO2 at 348.15 K and 30 MPa up to 20,000 a. Given the low shale porosity (5%), chemical reactions mostly occurred in the sandstone for the first 2000 a with some precipitation at the ss/sh interface. From 2000 to 4000 a, ankerite, dolomite and illite began replacing Mg–Fe chlorite at the sandstone/shale interface. Transformation of chlorite to ankerite is the dominant reaction occluding the shale porosity in most simulations: from 4000 to 7500 a, this carbonation seals the reservoir and terminates reaction. Overall, the carbonates (calcite, ankerite, dolomite), chlorite and goethite all remain close to local chemical equilibrium with brine. Quartz is almost inert from the point of its dissolution/precipitation. However, the rate of quartz reaction controls the long-term decline in aqueous silica activity and its evolution toward equilibrium. The reactions of feldspars and clays depend strongly on their reaction rate constants (microcline is closer to local equilibrium than albite). The timing of porosity occlusion mostly therefore depends on the kinetic constants of kaolinite and illite. For example, an increase in the kaolinite kinetic constant by 0.25 logarithmic units hastened porosity closure by 4300 a. The earliest simulated closure of porosity occurred at approximately 108 a for simulations designed as sensitivity tests for the rate constants.These simulations also emphasize that the rate of CO2 immobilization as aqueous bicarbonate (solubility trapping) or as carbonate minerals (mineral trapping) in sandstone reservoirs depends upon reaction kinetics – but the relative fraction of each trapped CO2 species only depends upon the initial chemical composition of the host sandstone. For example, at the point of porosity occlusion the fraction of bicarbonate remaining in solution depends upon the initial Na and K content in the host rock but the fraction of carbonate mineralization depends only on the Ca, Mg, Fe content. Since ankerite is the dominant mineral that occludes porosity, the dissolved concentration of ferrous iron is also an important parameter. Future efforts should focus on cross-comparisons and ground-truthing of simulations made for standard case studies as well as laboratory measurements of the reactivities of clay minerals. 相似文献
18.
Semi-empirical equations for the systematic decrease in permeability with depth in porous and fractured media 总被引：1，自引：0，他引：1
Permeability loss with depth is a general trend in geological media and plays an essential role in subsurface fluid flow and solute transport. In the near surface zone where groundwater movement is active, the decrease in permeability with depth is dominated by the mechanical compaction of deformable media caused by the increase in lithostatic stress with depth. Instead of using empirical equations from statistical analysis, by considering the well-defined relationships among permeability, porosity, fracture aperture and effective stress under lithostatic conditions, new semi-empirical equations for the systematic depth-dependent permeability are derived, as well as the equations for the depth-dependent porosity in a porous medium and the depth-dependent fracture aperture in a fractured medium. The existing empirical equations can be included in the new equations as special cases under some simplification. These new semi-empirical equations perform better than previous equations to interpret the depth-dependent permeability of the Pierre Shale (with a maximum depth of approximately 4,500 m) and the granite at Stripa, Sweden (with a maximum depth of about 2,500 m). 相似文献
19.
20.
Surface reaction versus diffusion control of mineral dissolution and growth rates in geochemical processes 总被引：2，自引：0，他引：2
The kinetics of chemical reactions at mineral surfaces and the rates of diffusion of species in an aqueous phase are coupled in many geochemical systems. Analytical solutions to equations describing coupled mineral dissolution/growth and solute transport in both transient and steady-state systems are used to delimit regimes of pure reaction control, pure transport control and mixed kinetic control of mass-transfer rates. The relative significance of the two processes depends on the magnitudes of the diffusion coefficients and rate constants as functions of temperature, and the degree of disequilibrium in the system. In addition, the system geometry, the ratio of mineral surface area to diffusion cross-section, and the porosity and tortuosity of the medium through which aqueous species diffuse affect reaction vs. diffusion control. In general, diffusion control increases with increasing temperature and increasing distance over which diffusion occurs. Calculations for the mixed kinetic regime in transient systems demonstrate that the relative significance of diffusion and surface reaction varies with reaction progress, and approaches a limiting value as equilibrium is approached. This limiting value may be appropriate to natural water-rock interactions that occur at conditions that are close to equilibrium. This result permits extension of simple models for irreversible mass transfer in homogeneous systems to systems in which mass-transfer kinetics are controlled by coupled surface reactions and mass transport. Criteria are established for time and length scales and fluid velocity limits on the validity of the continuum hypothesis and the local equilibrium assumption in mass-transport modeling. 相似文献