首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 724 毫秒
1.
Safety assessment of geosequestration of CO2 into deep saline aquifers requires a precise understanding of the study of hydro‐chemo‐mechanical couplings occurring in the rocks and the cement well. To this aim, a coupled chemo‐poromechanical model has been developed and implemented into a research code well‐suited to the resolution of fully coupled problems. This code is based on the finite volume methods. In a 1D axisymmetrical configuration, this study aims to simulate the chemo‐poromechanical behaviour of a system composed by the cement well and the caprock during CO2 injection. Major chemical reactions of carbonation occurring into cement paste and rocks are considered in order to evaluate the consequences of the presence of CO2 on the amount of dissolved matrix and precipitated calcium carbonates. The dissolution of the solid matrix is taken into account through the use of a chemical porosity. Matrix leaching and carbonation lead, as expected, to important variations of porosity, permeability and to alterations of transport properties and mechanical stiffness. These results justify the importance of considering a coupled analysis accounting for the main chemical reactions. It is worth noting that the modelling framework proposed in the present study could be extended to model the chemo‐poromechanical behaviour of the reservoir rock and the caprock when subjected to the presence of an acidic pore fluid (CO2‐rich brine). Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

2.
The mechanical properties of calcarenites are known to be significantly affected by water saturation: both stiffness and strength decrease for wetting in the short term and for chemical dissolution in the long term. Both processes mainly affect bonds among grains: immediately after inundation depositional bonds fall in suspension, whereas diagenetic bonds dissolve more slowly. In this paper, the authors started from the micro‐structural analysis of the weathering processes to conceive a strain hardening hydro‐chemo‐mechanical coupled elastoplastic constitutive model. The concept of extended hardening rules is here enriched: weathering functions have been determined by employing a micro to macro simplified upscaling procedure. Chemical damage is incorporated into the formulation by means of a scalar damage function. Its evolution is also described by using a multiscale approach. A new term is added to the strain rate tensor in order to incorporate the dissolution induced chemical deformations developing once the soft rock is turned into a granular material. A calibration procedure for the constitutive parameters is suggested, and the model is validated by using both coupled and uncoupled chemo‐mechanical experimental test results. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

3.
In this paper, a numerical model is developed for the fully coupled hydro‐mechanical analysis of deformable, progressively fracturing porous media interacting with the flow of two immiscible, compressible wetting and non‐wetting pore fluids, in which the coupling between various processes is taken into account. The governing equations involving the coupled solid skeleton deformation and two‐phase fluid flow in partially saturated porous media including cohesive cracks are derived within the framework of the generalized Biot theory. The fluid flow within the crack is simulated using the Darcy law in which the permeability variation with porosity because of the cracking of the solid skeleton is accounted. The cohesive crack model is integrated into the numerical modeling by means of which the nonlinear fracture processes occurring along the fracture process zone are simulated. The solid phase displacement, the wetting phase pressure and the capillary pressure are taken as the primary variables of the three‐phase formulation. The other variables are incorporated into the model via the experimentally determined functions, which specify the relationship between the hydraulic properties of the fracturing porous medium, that is saturation, permeability and capillary pressure. The spatial discretization is implemented by employing the extended finite element method, and the time domain discretization is performed using the generalized Newmark scheme to derive the final system of fully coupled nonlinear equations of the hydro‐mechanical problem. It is illustrated that by allowing for the interaction between various processes, that is the solid skeleton deformation, the wetting and the non‐wetting pore fluid flow and the cohesive crack propagation, the effect of the presence of the geomechanical discontinuity can be completely captured. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

4.
5.
6.
This paper presents a biconcave bond model to investigate the effect of the cementation between grains on the mechanical behavior of rock. The proposed model considers the shape of the bonds among particles that have a biconcave cement form, based on observations of microscopic rock images. The general equations of the proposed model are based on Dvorkin theory. The accuracy and efficiency of the bond model is improved in three ways. After the biconcave bond model is implemented in the discrete element method software Particle Flow Code in 2 Dimensions, a series of numerical uniaxial compression tests were performed to investigate the relationships between the micro‐ to macro‐parameters. The simulations revealed that the biconcave bond model reflects the effect of micro‐parameters, such as the elastic modulus and Poisson's ratio of the cement, on the macroscopic deformation of cemented granular material. Variations in the bond geometry caused extremely diverse macro‐mechanical behaviors. Experimental results concerning rock demonstrate that the biconcave bond model accurately captures the mechanical behavior of intact rock and supports an innovative method for investigating the relationships between the micro‐ and macro‐parameters of cemented granular material. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

7.
A significant increase of the permeability of concrete upon micro‐cracking and a good correlation between the evolution of damage (material stiffness) and permeability are observed experimentally. The present contribution investigates this correlation theoretically, with the help of lattice analyses. Scaling analysis of lattices which contain elastic brittle bonds has shown that the material degradation should be described by the evolution of the material stiffness, or compliance, in a continuum setting (damage models). This result is reviewed and further documented in the first part of the paper. In the second part, hydro‐mechanical problems are considered with the construction of a hydraulic lattice, dual to the mechanical one. We observe that the average permeability upon micro‐cracking is the lattice scale‐independent controlling variable in the hydraulic problem. Additionally, results show that in a continuum poro‐mechanical approach, the evolution of the material permeability ought to be related to the elastic unloading stiffness, described e.g. with the help of continuum damage variables. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

8.
This study describes the microstructural and chemical development of symplectites, obtained in fluid‐mediated mineral replacement experiments. During the experiments polymineralic feldspar‐rich samples were exposed to aqueous Na–SiO2 solution at 600 °C and 2 kbar confining pressures for durations of 12 h to 20 days. The resulting reaction rims display high mineralogical and structural complexity and contain two varieties of symplectites, represented by nanometre‐scale intergrowths of gehlenite–zeolite and grossular–zeolite grains. The experimental fluid was enriched in 18O isotope in order to trace oxygen redistribution during the reaction. The elevated 18O concentration in the reaction products and the heterogeneity in its distribution suggest that symplectite formation was controlled by dissolution–precipitation mechanisms rather than volume‐diffusion processes. Microstructural and chemical observations suggest that symplectite formation occurred in multiple stages in response to spatially heterogeneous and temporarily evolving fluid composition at the reaction interfaces. Hence, our results shed light on the fundamental processes involved in symplectite formation improving our ability to interpret symplectite microstructures.  相似文献   

9.
Optical, cathodoluminescence and transmission electron microscope (TEM) analyses were conducted on four groups of calcite fault rocks, a cataclastic limestone, cataclastic coarse-grained marbles from two fault zones, and a fractured mylonite. These fault rocks show similar microstructural characteristics and give clues to similar processes of rock deformation. They are characterized by the structural contrast between macroscopic cataclastic (brittle) and microscopic mylonitic (ductile) microstructures. Intragranular deformation microstructures (i.e. deformation twins, kink bands and microfractures) are well preserved in the deformed grains in clasts or in primary rocks. The matrix materials are of extremely fine grains with diffusive features. Dislocation microstructures for co-existing brittle deformation and crystalline plasticity were revealed using TEM. Tangled dislocations are often preserved at the cores of highly deformed clasts, while dislocation walls form in the transitions to the fine-grained  相似文献   

10.
A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

11.
This paper presents an investigation of the reactive transport of multicomponent chemicals in clays under coupled thermal, hydraulic, chemical and mechanical framework, considering the diffusion processes in detail. More specifically, combined effects due to the electrochemical and the thermal diffusion potentials are investigated. A theoretical framework for coupling thermal diffusion, i.e. the Soret effect, with electrochemical diffusion in a multi-ionic system is provided. An explicit form of a definition for the thermal diffusion coefficient in a multicomponent chemical transport model is developed. Chemical transport is linked to an advanced geochemical model, PHREEQC (version 2), in order to include chemical reactions. The effects of the combined diffusion potentials on the reactive transport of multicomponent chemicals are investigated by a series of numerical simulations of coupled thermal, hydraulic and chemical behaviour.  相似文献   

12.
This paper presents an upscaling concept of swelling/shrinking processes of a compacted bentonite/sand mixture, which also applies to swelling of porous media in general. A constitutive approach for highly compacted bentonite/sand mixture is developed accordingly. The concept is based on the diffuse double layer theory and connects microstructural properties of the bentonite as well as chemical properties of the pore fluid with swelling potential. Main factors influencing the swelling potential of bentonite, i.e. variation of water content, dry density, chemical composition of pore fluid, as well as the microstructures and the amount of swelling minerals are taken into account. According to the proposed model, porosity is divided into interparticle and interlayer porosity. Swelling is the potential of interlayer porosity increase, which reveals itself as volume change in the case of free expansion, or turns to be swelling pressure in the case of constrained swelling. The constitutive equations for swelling/shrinking are implemented in the software GeoSys/RockFlow as a new chemo‐hydro‐mechanical model, which is able to simulate isothermal multiphase flow in bentonite. Details of the mathematical and numerical multiphase flow formulations, as well as the code implementation are described. The proposed model is verified using experimental data of tests on a highly compacted bentonite/sand mixture. Comparison of the 1D modelling results with the experimental data evidences the capability of the proposed model to satisfactorily predict free swelling of the material under investigation. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

13.
武文华  李锡夔 《岩土力学》2008,29(5):1152-1158
结合化学反应方程式,并应用多孔多相介质溶混污染物输运过程的数值模型,对多孔多相介质中含均相/非均相化学反应传质过程进行了数值模拟。化学反应主要包含均相快速/慢速和非均相快速/慢速等5种化学反应过程,溶质输运行为的控制机制主要考虑对流、扩散及降解、吸附等。基于原有的隐式特征线Galerkin离散化的有限元方法,求解模型控制方程的边值初值问题,求解过程中把均相化学反应物质中按照反应物和生成物分开,非均相反应物质按照固相和液相分开,对均相反应物及非均相液相物质浓度耦合求解,而均相生成物和非均相固相物质独立求解。使方程组按照其不同类型进行分类,同时可减少未知数的个数。对于含有非线性内状态变量的右端项进行迭代求解。数值例题结果验证了所提出的数值方法的有效性、计算精度和稳定性。  相似文献   

14.
This paper is the first part of a work that aims at developing a mechanical model for the behaviour of propellant‐like materials under high confining pressure and strain rate. The behaviour of a typical material is investigated experimentally. Several microstructural deformation processes are identified and correlated with loading conditions. The resulting behaviour is complex, non‐linear, and characterized by multiple couplings. The general structure of a relevant model is sought using a thermodynamic framework. A viscoelastic‐viscoplastic‐compaction model structure is derived under suitable simplifying assumptions, in the framework of finite, though moderate, strains. Model development, identification and numerical applications are given in the companion paper. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

15.
Compaction and associated fluid flow are fundamental processes in sedimentary basin deformation. Purely mechanical compaction originates mainly from pore fluid expulsion and rearrangement of solid particles during burial, while chemo‐mechanical compaction results from Intergranular Pressure‐Solution (IPS) and represents a major mechanism of deformation in sedimentary basins during diagenesis. The aim of the present contribution is to provide a comprehensive 3D framework for constitutive and numerical modeling of purely mechanical and chemo‐mechanical compaction in sedimentary basins. Extending the concepts that have been previously proposed for the modeling of purely mechanical compaction in finite poroplasticity, deformation by IPS is addressed herein by means of additional viscoplastic terms in the state equations of the porous material. The finite element model integrates the poroplastic and poroviscoplastic components of deformation at large strains. The corresponding implementation allows for numerical simulation of sediments accretion/erosion periods by progressive activation/deactivation of the gravity forces within a fictitious closed material system. Validation of the numerical approach is assessed by means of comparison with closed‐form solutions derived in the context of a simplified compaction model. The last part of the paper presents the results of numerical basin simulation performed in one dimensional setting, demonstrating the ability of the modeling to capture the main features in elastoplastic and viscoplastic compaction. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

16.
Based on relevant experimental data of a petroleum cement paste under mechanical loading and chemical leaching, an elastic‐plastic model is first proposed by taking into account plastic shearing and pore collapse. The degradation of mechanical properties induced by the chemical leaching is characterized by a chemical damage variable which is defined as the increase of porosity. Both elastic and plastic properties of the cement paste are affected by the chemical damage. The proposed model is calibrated from and applied to describe mechanical responses in triaxial compression tests respectively on sound and fully leached samples. In the second part, a phenomenological chemical model is defined to establish the relationship between porosity change and calcium dissolution process. The dissolution kinetics is governed by a diffusion law taking into account the variation of diffusion coefficient with calcium concentration. The chemical model is coupled with the mechanical model, and both are applied to describe mechanical response of cement paste samples subjected to progressive chemical leaching and compressive stresses. Comparisons between experimental data and numerical results are presented.  相似文献   

17.
低温下非静水应力作用固体颗粒间流体因其参与流体—矿物反应和流体传质而在地壳岩石变形和同构造地球化学过程中起着重要作用。粒间流体的性状、粒间流体—矿物反应和粒间流体传质都与固体颗粒所受应力密切相关。流体膜是非静水应力作用下固体颗粒间流体存在的主要形式,流体膜的厚度和成分决定着究竟是流体膜中水化力、双电层排斥力,还是因双电层而引起的流体渗透压在流体膜传递应力的机制中起主导作用。应力通过改变固—液界面上固体物质的自由能而影响固—液界面反应和流体的传质,并对这一应力感生过程的热力学和动力学方程进行了推演。  相似文献   

18.
The propagation characteristic of Rayleigh waves in a fluid‐saturated non‐homogeneous poroelastic half‐plane is addressed. Based on Biot's theory for fluid‐saturated media, which takes the inertia, fluid viscosity, mechanical coupling, compressibility of solid grains, and fluid into account, the dispersion equations of Rayleigh waves in fluid‐saturated non‐homogeneous soils/rocks are established. By considering the shear modulus of solid skeleton variation with depth exponentially, a small parameter, which reflects the relative change of shear modulus, is introduced. The asymptotic solution of the dispersion equation expressing the relationship between the phase velocity and wave number is obtained by using the perturbation method. In order to analyze the effects of non‐homogeneity on the propagation characteristic of Rayleigh waves, the variation of the phase velocity with the wave number is presented graphically and discussed through numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

19.
Metamorphic equilibration requires chemical communication between minerals and may be inhibited through sluggish volume diffusion and or slow rates of dissolution in a fluid phase. Relatively slow diffusion and the perceived robust nature of chemical growth zoning may preclude garnet porphyroblasts from readily participating in low‐temperature amphibolite facies metamorphic reactions. Garnet is widely assumed to be a reactant in staurolite‐isograd reactions, and the evidence for this has been assessed in the Late Proterozoic Dalradian pelitic schists of the Scottish Highlands. The 3D imaging of garnet porphyroblasts in staurolite‐bearing schists reveals a good crystal shape and little evidence of marginal dissolution; however, there is also lack of evidence for the involvement of either chlorite or chloritoid in the reaction. Staurolite forms directly adjacent to the garnet, and its nucleation is strongly associated with deformation of the muscovite‐rich fabrics around the porphyroblasts. “Cloudy” fluid inclusion‐rich garnet forms in both marginal and internal parts of the garnet porphyroblast and is linked both to the production of staurolite and to the introduction of abundant quartz inclusions within the garnet. Such cloudy garnet typically has a Mg‐rich, Mn‐poor composition and is interpreted to have formed during a coupled dissolution–reprecipitation process, triggered by a local influx of fluid. All garnet in the muscovite‐bearing schists present in this area is potentially reactive, irrespective of the garnet composition, but very few of the schists contain staurolite. The staurolite‐producing reaction appears to be substantially overstepped during the relatively high‐pressure Barrovian regional metamorphism reflecting the limited permeability of the schists in peak metamorphic conditions. Fluid influx and hence reaction progress appear to be strongly controlled by subtle differences in deformation history. The remaining garnet fails to achieve chemical equilibrium during the reaction creating distinctive patchy compositional zoning. Such zoning in metamorphic garnet created during coupled dissolution–reprecipitation reactions may be difficult to recognize in higher grade pelites due to subsequent diffusive re‐equilibration. Fundamental assumptions about metamorphic processes are questioned by the lack of chemical equilibrium during this reaction and the restricted permeability of the regional metamorphic pelitic schists. In addition, the partial loss of prograde chemical and textural information from the garnet porphyroblasts cautions against their routine use as a reliable monitor of metamorphic history. However, the partial re‐equilibration of the porphyroblasts during coupled dissolution–reprecipitation opens possibilities of mapping reaction progress in garnet as a means of assessing fluid access during peak metamorphic conditions.  相似文献   

20.
Reactive fluid-flow experiments in fractures subjected to normal stress suggest the potential for either increased or decreased permeability resulting from fracture-surface dissolution. We present a computational model that couples mechanical deformation and chemical alteration of fractures subjected to constant normal stress and reactive fluid flow. The model explicitly represents micro-scale roughness of the fracture surfaces and calculates elastic deformation of the rough surfaces using a semi-analytical approach that ensures the surfaces remain in static equilibrium. A depth-averaged reactive transport model calculates chemical alteration of the surfaces, which leads to alteration of the contacting fracture surfaces. The mechanical deformation and chemical alteration calculations are explicitly coupled, which is justified by the disparate timescales required for equilibration of mechanical stresses and reactive transport processes. An idealized analytical representation of dissolution from a single contacting asperity shows that under reaction-limited conditions, contacting asperities can dissolve faster than the open regions of the fracture. Computational simulations in fractures with hundreds of contacting asperities show that the transition from transport-limited conditions (low flow rates) to reaction-rate-limited conditions (high flow rates) causes a shift from monotonically increasing permeability to a more complicated process in which permeability initially decreases and then increases as contacting asperities begin to dissolve. These results are qualitatively consistent with a number of experimental observations reported in the literature and suggest the potential importance of the relative magnitude of mass transport and reaction kinetics on the evolution of fracture permeability in fractures subjected to combined normal stress and reactive fluid flow.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号