首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 718 毫秒
1.
Rate‐dependent behaviour of chalk and other porous rocks has undergone widespread study in geomechanics due to its implications on the performance of engineering structures. We present a rate‐dependent constitutive model for chalk and other porous rocks with several new features. The model formulation is based on a viscoplastic rate‐lines approach in which the axial strain rate depends on the proximity of the stress point to an elliptical reference surface. A non‐associated viscoplastic potential surface and an axial scaling algorithm are used to determine the viscoplastic strain components. The model predicts that axial yields stress varies as a power function of applied axial strain rate, as shown by published laboratory data. Comparisons with published experimental data indicate that the model is capable of reproducing observed rate‐dependent behaviour of chalk under a variety of loading conditions. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

2.
Late Palaeozoic‐age strata from the Capitan Reef in west Texas show facies‐dependent heterogeneity in the sulphur isotopic composition of carbonate‐associated sulphate, which is trace sulphate incorporated into carbonate minerals that is often used to reconstruct the sulphur isotopic composition of ancient seawater. However, diagenetic pore fluid processes may influence the sulphur isotopic composition of carbonate‐associated sulphate. These processes variously modify the sulphur isotopic composition of incorporated sulphate from syndepositional seawater in shelf crest, outer shelf, shelf margin and slope depositional settings. This study used a new multicollector inductively‐coupled plasma mass spectrometry technique to determine the sulphur isotopic composition of samples of individual depositional and diagenetic textures. Carbonate rocks representing peritidal facies in the Yates and Tansill formations preserve the sulphur isotopic composition of Guadalupian seawater sulphate despite alteration of the carbon and oxygen isotopic compositions by meteoric and dolomitizing diagenetic processes. However, sulphur isotopic data indicate that limestones deposited in reef and slope facies in the Capitan and Bell Canyon formations largely incorporate sulphate from anoxic marine‐phreatic pore fluids isotopically modified from seawater by microbial sulphate reduction, despite generally preserving the carbon and oxygen isotopic compositions of Permian seawater. Some early and all late meteoric calcite cements have carbonate‐associated sulphate with a sulphur isotopic composition distinct from that of Permian seawater. Detailed petrographic and sedimentary context for carbonate‐associated sulphate analyses will allow for improved reconstructions of ancient seawater composition and diagenetic conditions in ancient carbonate platforms. The results of this study indicate that carbonate rocks that diagenetically stabilize in high‐energy environments without pore fluid sulphate gradients can provide a robust archive of ancient seawater's sulphur isotopic composition.  相似文献   

3.
Heterogeneities, such as fractures and cracks, are ubiquitous in porous rocks. Mesoscopic heterogeneities, that is, heterogeneities on length scales much larger than typical pore size but much smaller than the wavelength, are increasingly believed to be responsible for significant wave energy loss in the seismic frequency band. When a compressional wave stresses a material containing mesoscopic heterogeneities, the more compliant parts of the material (e.g., fractures and cracks) respond with a greater fluid pressure than the stiffer portions (e.g., matrix pores). The induced fluid flow, resulting from the pressure gradients developed on such scale, is called mesoscopic flow. In the present study, the double‐porosity dual‐permeability model is adopted to incorporate mesoscopic heterogeneities into rock models to account for the attenuation of wave energy. Based on the model, the damping effect due to mesoscopic flow in a one‐dimensional porous structure is investigated. Analytical solutions for several boundary‐value problems are obtained in the frequency domain. The dynamic responses of infinite and finite porous layer are examined. Numerical calculations show that the damping effect of mesoscopic flow is significant on the pore pressure response and the resulting effective stress. For the displacement, the effect is seen only at the very low frequency range or near the resonance frequencies. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

4.
Determination of transport properties of geomaterials is an important issue in many fields of engineering analysis and design. For example, in petroleum engineering, in situ permeability of an oil reservoir may be crucial in establishing its viability for exploitation, whilst prevention of leakage from underground storage facilities for oil and gas, nuclear waste as well as viability of CO2 sequestration projects crucially depends on its long‐term values. Permeability is indirectly related to the porosity, pore‐size distribution and pore architecture of the porous media. These parameters evolve when a strain field is imposed. Physical measurement of permeability under a strain field in laboratory conditions is difficult, expensive and prone to a number of uncertainties. In the past, pore network models have been used to compute permeability of materials under stress/strain‐free conditions. In this paper, we propose an enhanced pore network model to compute permeability of rocks and rock‐like porous media under a stress/strain field. Data of pore‐size distribution obtained from mercury intrusion porosimetry are used to compute permeability of rock samples from various unspecified oilfields in the world. It is shown that the two permeabilities can be predicted from the model with sufficient accuracy. A hypothesis for change in porosity, pore‐size distribution and pore architecture as a result of imposed mechanical strains is then proposed. Based on this, permeability is computed again for one of the rock samples under uniaxial and triaxial compressive and tensile strain fields. It is shown that depending on the state of strain field imposed, permeability evolves in an anisotropic manner. Permeability under tensile strain field increases dramatically compared with the reduction that takes place under compressive strain field of the same magnitude. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

5.
At sub‐arc depths, the release of carbon from subducting slab lithologies is mostly controlled by fluid released by devolatilization reactions such as dehydration of antigorite (Atg‐) serpentinite to prograde peridotite. Here we investigate carbonate–silicate rocks hosted in Atg‐serpentinite and prograde chlorite (Chl‐) harzburgite in the Milagrosa and Almirez ultramafic massifs of the palaeo‐subducted Nevado‐Filábride Complex (NFC, Betic Cordillera, S. Spain). These massifs provide a unique opportunity to study the stability of carbonate during subduction metamorphism at P–T conditions before and after the dehydration of Atg‐serpentinite in a warm subduction setting. In the Milagrosa massif, carbonate–silicate rocks occur as lenses of Ti‐clinohumite–diopside–calcite marbles, diopside–dolomite marbles and antigorite–diopside–dolomite rocks hosted in clinopyroxene‐bearing Atg‐serpentinite. In Almirez, carbonate–silicate rocks are hosted in Chl‐harzburgite and show a high‐grade assemblage composed of olivine, Ti‐clinohumite, diopside, chlorite, dolomite, calcite, Cr‐bearing magnetite, pentlandite and rare aragonite inclusions. These NFC carbonate–silicate rocks have variable CaO and CO2 contents at nearly constant Mg/Si ratio and high Ni and Cr contents, indicating that their protoliths were variable mixtures of serpentine and Ca‐carbonate (i.e., ophicarbonates). Thermodynamic modelling shows that the carbonate–silicate rocks attained peak metamorphic conditions similar to those of their host serpentinite (Milagrosa massif; 550–600°C and 1.0–1.4 GPa) and Chl‐harzburgite (Almirez massif; 1.7–1.9 GPa and 680°C). Microstructures, mineral chemistry and phase relations indicate that the hybrid carbonate–silicate bulk rock compositions formed before prograde metamorphism, likely during seawater hydrothermal alteration, and subsequently underwent subduction metamorphism. In the CaO–MgO–SiO2 ternary, these processes resulted in a compositional variability of NFC serpentinite‐hosted carbonate–silicate rocks along the serpentine‐calcite mixing trend, similar to that observed in serpentinite‐hosted carbonate‐rocks in other palaeo‐subducted metamorphic terranes. Thermodynamic modelling using classical models of binary H2O–CO2 fluids shows that the compositional variability along this binary determines the temperature of the main devolatilization reactions, the fluid composition and the mineral assemblages of reaction products during prograde subduction metamorphism. Thermodynamic modelling considering electrolytic fluids reveals that H2O and molecular CO2 are the main fluid species and charged carbon‐bearing species occur only in minor amounts in equilibrium with carbonate–silicate rocks in warm subduction settings. Consequently, accounting for electrolytic fluids at these conditions slightly increases the solubility of carbon in the fluids compared with predictions by classical binary H2O–CO2 fluids, but does not affect the topology of phase relations in serpentinite‐hosted carbonate‐rocks. Phase relations, mineral composition and assemblages of Milagrosa and Almirez (meta)‐serpentinite‐hosted carbonate–silicate rocks are consistent with local equilibrium between an infiltrating fluid and the bulk rock composition and indicate a limited role of infiltration‐driven decarbonation. Our study shows natural evidence for the preservation of carbonates in serpentinite‐hosted carbonate–silicate rocks beyond the Atg‐serpentinite breakdown at sub‐arc depths, demonstrating that carbon can be recycled into the deep mantle.  相似文献   

6.
An analytical solution is proposed for transient flow and deformation coupling of a fluid‐saturated poroelastic medium within a finite two‐dimensional (2‐D) rectangular domain. In this study, the porous medium is assumed to be isotropic, homogeneous, and compressible. In addition, the point sink can be located at an arbitrary position in the porous medium. The fluid–solid interaction in porous media is governed by the general Biot's consolidation theory. The method of integral transforms is applied in the analytical formulation of closed‐form solutions. The proposed analytical solution is then verified against both exact and numerical results. The analytical solution is first simplified and validated by comparison with an existing exact solution for the uncoupled problem. Then, a case study for pumping from a confined aquifer is performed. The consistency between the numerical solution and the analytical solution confirms the accuracy and reliability of the analytical solution presented in this paper. The proposed analytical solution can help us to obtain in‐depth insights into time‐dependent mechanical behavior due to fluid withdrawal within finite 2‐D porous media. Moreover, it can also be of great significance to calibrate numerical solutions in plane strain poroelasticity and to formulate relevant industry norms and standards. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

7.
In the present study, we have developed a numerical method which can simulate the dynamic behaviour of a seabed ground during gas production from methane hydrate‐bearing sediments. The proposed method can describe the chemo‐thermo‐mechanical‐seismic coupled behaviours, such as phase changes from hydrates to water and gas, temperature changes and ground deformation related to the flow of pore fluids during earthquakes. In the first part of the present study, the governing equations for the proposed method and its discretization are presented. Then, numerical analyses are performed for hydrate‐bearing sediments in order to investigate the dynamic behaviour during gas production. The geological conditions and the material parameters are determined using the data of the seabed ground at Daini‐Atsumi knoll, Eastern Nankai Trough, Japan, where the first offshore production test of methane hydrates was conducted. A predicted earthquake at the site is used in the analyses. Regarding the seismic response to the earthquake which occur during gas production process, the wave profiles of horizontal acceleration and horizontal velocity were not extensively affected by the gas production. Hydrate dissociation behaviour is sensitive to changes in the pore pressure during earthquakes. Methane hydrate dissociation temporarily became active in some areas because of the main motion of the earthquake, then methane hydrate dissociation brought about an increase in the average pressure of the fluids during the earthquake. And, it was this increase in average pore pressure that finally caused the methane hydrate dissociation to cease during the earthquake. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

8.
Theoretical analysis and computational simulations have been carried out to investigate how medium and pore‐fluid compressibility affects the chemical‐dissolution front propagation, which is associated with a fully‐coupled nonlinear problem between porosity, pore‐fluid pressure, pore‐fluid density and reactive chemical‐species transport within a deformable and fluid‐saturated porous medium. When the fully‐coupled nonlinear system is in a subcritical state, some analytical solutions have been derived for a special case, in which the ratio of the equilibrium concentration to the solid molar density of the chemical species is approaching zero. To investigate the effect of either medium compressibility or pore‐fluid compressibility on the evolutions of chemical dissolution fronts in supercritical chemical dissolution systems, numerical algorithms and procedures have been also proposed. The related theoretical and numerical results have demonstrated that: (i) not only can pore‐fluid compressibility affect the propagating speeds of chemical dissolution fronts in both subcritical and supercritical systems, but also it can affect the growth and amplitudes of irregular chemical dissolution fronts in supercritical systems; (ii) medium compressibility may have a little influence on the propagating speeds of chemical dissolution fronts, but it can have significant effects on the growth and amplitudes of irregular chemical dissolution fronts in supercritical systems; and (iii) both medium and pore‐fluid compressibility may stabilize irregular chemical‐dissolution‐fronts in supercritical chemical dissolution systems. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

9.
沉积盆地流体活动及其成岩响应   总被引：6，自引：0，他引：6  
盆地流体活动通过改变成岩场的温压和孔隙水介质条件,进而影响到盆地内成岩作用序列、水—岩相互作用以及储层孔隙变化。对于碎屑岩而言,流体活动仅能导致母岩中少部分矿物的溶蚀以及孔隙内溶解与沉淀作用发生,富含二氧化碳和有机酸的流体活动为次生孔隙发育提供条件,超压环境下流体幕式释放既可导致超压体内次生孔隙的发育,又可以导致热流体活动影响区域内浅部储层的成岩异常;对于碳酸盐岩而言,由于碳酸盐岩具有很强的溶解性,盆地流体活动导致碳酸盐岩储层组构及空间发生巨大的变化,其中天水的下渗以及与深部热液或热流体沿断裂或裂缝释放均可导致碳酸盐岩成分、组构、储集空间发生根本性变化,故与盆地流体活动密切相关的风化壳型岩溶储层和构造热液白云岩储层构成碳酸盐岩中最为重要的储层类型。然而,盆地流体活动及其成岩响应研究仍存在许多亟待解决的难点问题,急待深入的研究探索。  相似文献   

10.
Air sparging (AS) is an in situ soil/groundwater remediation technology, which involves the injection of pressurized air/oxygen through an air sparging well below the zone of contamination. Characterizing the mechanisms governing movement of air through saturated porous media is critical for the design of an effective cleanup treatment system. In this research, micromechanical investigation was performed to understand the physics of air migration and subsequent spatial distribution of air at pore scale during air sparging. The void space in the porous medium was first characterized by pore network consisting of connected pore bodies and bonds. The biconical abscissa asymmetric concentric bond was used to describe the connection between two adjacent pore bodies. Then a rule‐based dynamic two‐phase flow model was developed and applied to the pore network model. A forward integration of time was performed using the Euler scheme. For each time step, the effective viscosity of the fluid was calculated based on fractions of two phases in each bond, and capillary pressures across the menisci was considered to compute the pressure field. The developed dynamic model was used to study the rate‐dependent drainage during air sparging. The effect of the capillary number and geometrical properties of the network on the dynamic flow properties of two‐phase flow including residual saturation, spatial distribution of air and water, dynamic phase transitions, and relative permeability‐capillary pressure curves were systematically investigated. Results showed that all the above information for describing the air water two‐phase flow are not intrinsic properties of the porous medium but are affected by the two‐phase flow dynamics and spatial distribution of each phase, providing new insight to air sparging. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

11.
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.  相似文献   

12.
A three-phase hydro-mechanical model for hydraulic fracturing is proposed. Three phases include: porous solid, fracturing fluid and host fluid. Discontinuity is handled using extended finite element method (XFEM) while cohesive crack model is used as fracturing criterion. Flow through fracture is defined as one-dimensional laminar flow, and flow through porous medium (host reservoir) is defined as two-dimensional Darcy flow. Coupling between two fluids in each space, fracture and pore, is captured through capillary pressure–saturation relationship, while the identical fluids in fracture and pore are coupled through a so-called leak-off mass transfer term. Coupling between fluids and deformation is captured through compatibility of volumetric strain of fluids within fracture and pore, and volumetric strain of the matrix. Spatial and temporal discretisation is achieved using the standard Galerkin method and the finite difference technique, respectively. The model is verified against analytical solutions available from literature. The leaking of fracturing fluid into the medium and suction of porous fluid into the fracture around the tip, are investigated. Sensitivity analyses are carried out for cases with slow and fast injection rates. It is shown that the results by single-phase flow may underestimate the leak-off.  相似文献   

13.
Niobium and Ta concentrations in ultrahigh‐pressure (UHP) eclogites and rutile from these eclogites and associated high pressure (HP) veins were used to study the behaviour of Nb–Ta during dehydration and fluid‐rock interaction. Samples were collected through a ～2 km profile at the Bixiling complex in the Dabie orogenic belt, Central‐Eastern China. All but one eclogite away from veins (EAVs) display nearly constant Nb/Ta ratios ranging from 16.1 to 19.2, with an average of 16.9 ± 0.8 (2 SE), similar to that of their gabbroic protolith from the Yangtze Block. Nb/Ta ratios of rutile from the EAVs range from 12.7 to 25.3 among different individual grains, with the average values close to those of the corresponding bulk rocks. These observations show that Nb and Ta were not significantly fractionated by prograde metamorphism up to eclogite facies when no significant fluid‐rock interaction occurs. In contrast, Nb/Ta ratios of rutile from eclogites close to veins (ECVs) are highly variable from 17.8 to 49.8, which are systematically higher (by up to 17) than those of rutile from the veins. These observations demonstrate that Nb and Ta were mobilized and fractionated during localized fluid flow and intensive fluid‐rock interaction. This is strongly supported by Nb/Ta zoning patterns in single rutile grains revealed by in situ LA‐ICP‐MS analysis. Ratios of Nb/Ta in the ECV‐hosted rutile decrease gradually from cores towards rims, whereas those in the EAV‐hosted rutile are nearly invariable. Furthermore, the vein rutile shows Nb/Ta zoning patterns that are complementary to those in rutile from their immediate hosts (ECVs), suggesting an internal origin for the vein‐forming fluids. The Nb/Ta ratios of such fluids evolved from low values at the early stage of subduction to higher values at later supercritical conditions with increased temperature and pressure. Quantitative modelling was conducted to constrain the compositional evolution of metamorphic fluids during dehydration and fluid‐rock interaction focusing on Nb–Ta distribution. The modelling results based on our proposed multistage fluid phase evolution path can essentially reproduce the natural observations reported in the present study.  相似文献   

14.
The effective stress concept for solid‐fluid 2‐phase media was revisited in this work. In particular, the effects of the compressibility of both the pore fluid and the soil particles were studied under 3 different conditions, i.e., undrained, drained, and unjacketed conditions based on a Biot‐type theory for 2‐phase porous media. It was confirmed that Terzaghi effective stress holds at the moment when soil grains are assumed to be incompressible and when the compressibility of the pore fluid is small enough compared to that of the soil skeleton. Then, isotropic compression tests for dry sand under undrained conditions were conducted within the triaxial apparatus in which the changes in the pore air pressure could be measured. The ratio of the increment in the cell pressure to the increment in the pore air pressure, m, corresponds to the inverse of the B value by Bishop and was obtained during the step loading of the cell pressure. In addition, the m values were evaluated by comparing them with theoretically obtained values based on the solid‐fluid 2‐phase mixture theory. The experimental m values were close to the theoretical values, as they were in the range of approximately 40 to 185, depending on the cell pressure. Finally, it was found that the soil material with a highly compressible pore fluid, such as air, must be analyzed with the multi‐phase porous mixture theory. However, Terzaghi effective stress is practically applicable when the compressibilities of both the soil particles and the pore fluid are small enough compared to that of the soil skeleton.  相似文献   

15.
In this paper, a fully coupled numerical model is presented for the finite element analysis of the deforming porous medium interacting with the flow of two immiscible compressible wetting and non-wetting pore fluids. The governing equations involving coupled fluid flow and deformation processes in unsaturated soils are derived within the framework of the generalized Biot theory. The displacements of the solid phase, the pressure of the wetting phase and the capillary pressure are taken as the primary unknowns of the present formulation. The other variables are incorporated into the model using the experimentally determined functions that define the relationship between the hydraulic properties of the porous medium, i.e. saturation, relative permeability and capillary pressure. It is worth mentioning that the imposition of various boundary conditions is feasible notwithstanding the choice of the primary variables. The modified Pastor–Zienkiewicz generalized constitutive model is introduced into the mathematical formulation to simulate the mechanical behavior of the unsaturated soil. The accuracy of the proposed mathematical model for analyzing coupled fluid flows in porous media is verified by the resolution of several numerical examples for which previous solutions are known. Finally, the performance of the computational algorithm in modeling of large-scale porous media problems including the large elasto-plastic deformations is demonstrated through the fully coupled analysis of the failure of two earth and rockfill dams. Furthermore, the three-phase model is compared to its simplified one which simulates the unsaturated porous medium as a two-phase one with static air phase. The paper illustrates the shortcomings of the commonly used simplified approach in the context of seismic analysis of two earth and rockfill dams. It is shown that accounting the pore air as an independent phase significantly influences the unsaturated soil behavior.  相似文献   

16.
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.  相似文献   

17.
We have examined the data of 600 geothermal wells and re‐interpreted approximately 500 km seismic profiles through the field of Larderello, Italy. We conclude that the two main seismic reflectors present below the geothermal area host two different fluids: (1) superheated steam in the upper H‐horizon (reached by drillholes) and (2) supercritical fluid in the deeper K‐horizon (reached by few unproductive or damaged wells). The superheated steam has the physical and chemical connotation of the geothermal fluid exploited so far at Larderello, whereas the supercritical fluid represents a potential unconventional deep‐seated resource still to be assessed. The high temperatures existing in correspondence of the K‐horizon suggest that the silica‐rich rocks are close to a plastic state and the fluids should remain confined in a medium sealed to the confining rocks, unless occasional fluid overpressure and abrupt high strain rates occur.  相似文献   

18.
The chemical‐dissolution front propagation problem exists ubiquitously in many scientific and engineering fields. To solve this problem, it is necessary to deal with a coupled system between porosity, pore‐fluid pressure and reactive chemical‐species transport in fluid‐saturated porous media. Because there was confusion between the average linear velocity and the Darcy velocity in the previous study, the governing equations and related solutions of the problem are re‐derived to correct this confusion in this paper. Owing to the morphological instability of a chemical‐dissolution front, a numerical procedure, which is a combination of the finite element and finite difference methods, is also proposed to solve this problem. In order to verify the proposed numerical procedure, a set of analytical solutions has been derived for a benchmark problem under a special condition where the ratio of the equilibrium concentration to the solid molar density of the concerned chemical species is very small. Not only can the derived analytical solutions be used to verify any numerical method before it is used to solve this kind of chemical‐dissolution front propagation problem but they can also be used to understand the fundamental mechanisms behind the morphological instability of a chemical‐dissolution front during its propagation within fluid‐saturated porous media. The related numerical examples have demonstrated the usefulness and applicability of the proposed numerical procedure for dealing with the chemical‐dissolution front instability problem within a fluid‐saturated porous medium. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

19.
This paper mainly deals with the theoretical aspects of chemical dissolution‐front instability problems in two‐dimensional fluid‐saturated porous media under non‐isothermal conditions. In the case of the mineral dissolution ratio (that is defined as the ratio of the dissolved‐mineral equilibrium concentration in the pore fluid to the molar concentration of the dissolvable mineral in the solid matrix of the fluid‐saturated porous medium) approaching zero, the corresponding critical condition has been mathematically derived when temperature variation effects are considered. As a complementary tool, the computational simulation method is used to simulate the morphological evolution of chemical dissolution fronts in two‐dimensional fluid‐saturated porous media under non‐isothermal conditions. The related theoretical and numerical results have demonstrated that: (i) a temperature increase in a non‐isothermal chemical dissolution system can have some influence on the propagation speed of the planar chemical dissolution front in the system. Generally, the chemical dissolution front in the non‐isothermal chemical dissolution system propagates slower than that in the counterpart isothermal chemical dissolution system when the temperature of the non‐isothermal chemical dissolution system is higher than that of the counterpart isothermal chemical dissolution system; (ii) a temperature increase in the non‐isothermal chemical dissolution system can stabilize the chemical dissolution front propagating in the system, because it can cause a decrease in the Zhao number of the system but does not affect the critical Zhao number of the system; and (iii) the temperature gradient in the upstream direction of a chemical dissolution front is smaller than that in the downstream direction of the chemical dissolution front when the non‐isothermal chemical dissolution system is supercritical. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

20.
A finite element algorithm for frictionless contact problems in a two‐phase saturated porous medium, considering finite deformation and inertia effects, has been formulated and implemented in a finite element programme. The mechanical behaviour of the saturated porous medium is predicted using mixture theory, which models the dynamic advection of fluids through a fully saturated porous solid matrix. The resulting mixed formulation predicts all field variables including the solid displacement, pore fluid pressure and Darcy velocity of the pore fluid. The contact constraints arising from the requirement for continuity of the contact traction, as well as the fluid flow across the contact interface, are enforced using a penalty approach that is regularised with an augmented Lagrangian method. The contact formulation is based on a mortar segment‐to‐segment scheme that allows the interpolation functions of the contact elements to be of order N. The main thrust of this paper is therefore how to deal with contact interfaces in problems that involve both dynamics and consolidation and possibly large deformations of porous media. The numerical algorithm is first verified using several illustrative examples. This algorithm is then employed to solve a pipe‐seabed interaction problem, involving large deformations and dynamic effects, and the results of the analysis are also compared with those obtained using a node‐to‐segment contact algorithm. The results of this study indicate that the proposed method is able to solve the highly nonlinear problem of dynamic soil–structure interaction when coupled with pore water pressures and Darcy velocity. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

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

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