共查询到20条相似文献，搜索用时 31 毫秒
1.
An Approach for Modeling Rock Discontinuous Mechanical Behavior Under Multiphase Fluid Flow Conditions 总被引：1，自引：0，他引：1
Peng-Zhi Pan Jonny Rutqvist Xia-Ting Feng Fei Yan 《Rock Mechanics and Rock Engineering》2014,47(2):589-603
In this paper, the two computer codes TOUGH2 and RDCA (for “rock discontinuous cellular automaton”) are integrated for coupled hydromechanical analysis of multiphase fluid flow and discontinuous mechanical behavior in heterogeneous rock. TOUGH2 is a well-established code for geohydrological analysis involving multiphase, multicomponent fluid flow and heat transport; RDCA is a numerical model developed for simulating the nonlinear and discontinuous geomechanical behavior of rock. The RDCA incorporates the discontinuity of a fracture independently of the mesh, such that the fracture can be arbitrarily located within an element, while the fluid pressure calculated by TOUGH2 can be conveniently applied to fracture surfaces. We verify and demonstrate the coupled TOUGH–RDCA simulator by modeling a number of simulation examples related to coupled multiphase flow and geomechanical processes associated with the deep geological storage of carbon dioxide—including modeling of ground surface uplift, stress-dependent permeability, and the coupled multiphase flow and geomechanical behavior of fractures intersecting the caprock. 相似文献
2.
In the context of carbon capture and storage, deep underground injection of CO2 induces the geomechanical changes within and around the injection zone and their impact on CO2 storage security should be evaluated. In this study, we conduct coupled multiphase fluid flow and geomechanical modeling to investigate such geomechanical changes, focusing on probabilistic analysis of injection-induced fracture reactivation (such as shear slip) that could lead to enhanced permeability and CO2 migration across otherwise low-permeability caprock formations. Fracture reactivation in terms of shear slip was analyzed by implicitly considering the fracture orientations generated using the Latin hypercube sampling method, in one case using published fracture statistics from a CO2 storage site. The analysis was conducted by a coupled multiphase fluid flow and geomechanical simulation to first calculate the three-dimensional stress evolution during a hypothetical CO2 injection operation and then evaluate the probability of shear slip considering the statistical fracture distribution and a Coulomb failure analysis. We evaluate the probability of shear slip at different points within the injection zone and in the caprock just above the injection zone and relate this to the potential for opening of new flow paths through the caprock. Our analysis showed that a reverse faulting stress field would be most favorable for avoiding fracture shear reactivation, but site-specific analyses will be required because of strong dependency of the local stress field and fracture orientations. 相似文献
3.
We present a discussion of the state-of-the-art on the use of discrete fracture networks (DFNs) for modelling geometrical characteristics, geomechanical evolution and hydromechanical (HM) behaviour of natural fracture networks in rock. The DFN models considered include those based on geological mapping, stochastic generation and geomechanical simulation. Different types of continuum, discontinuum and hybrid geomechanical models that integrate DFN information are summarised. Numerical studies aiming at investigating geomechanical effects on fluid flow in DFNs are reviewed. The paper finally provides recommendations for advancing the modelling of coupled HM processes in fractured rocks through more physically-based DFN generation and geomechanical simulation. 相似文献
4.
Jonny Rutqvist Liange Zheng Fei Chen Hui-Hai Liu Jens Birkholzer 《Rock Mechanics and Rock Engineering》2014,47(1):167-186
This paper presents simulation results related to coupled thermal–hydraulic–mechanical (THM) processes in engineered barrier systems (EBS) and clay host rock, in one case considering a possible link to geochemistry. This study is part of the US DOE Office of Nuclear Energy’s used fuel disposition campaign, to investigate current modeling capabilities and to identify issues and knowledge gaps associated with coupled THMC processes and EBS–rock interactions associated with repositories hosted in clay rock. In this study, we simulated a generic repository case assuming an EBS design with waste emplacement in horizontal tunnels that are back-filled with bentonite-based swelling clay as a protective buffer and heat load, derived for one type of US reactor spent fuel. We adopted the Barcelona basic model (BBM) for modeling of the geomechanical behavior of the bentonite, using properties corresponding to the FEBEX bentonite, and we used clay host rock properties derived from the Opalinus clay at Mont Terri, Switzerland. We present results related to EBS host–rock interactions and geomechanical performance in general, as well as studies related to peak temperature, buffer resaturation and thermally induced pressurization of host rock pore water, and swelling pressure change owing to variation of chemical composition in the EBS. Our initial THM modeling results show strong THM-driven interactions between the bentonite buffer and the low-permeability host rock. The resaturation of the buffer is delayed as a result of the low rock permeability, and the fluid pressure in the host rock is strongly coupled with the temperature changes, which under certain circumstances could result in a significant increase in pore pressure. Moreover, using the BBM, the bentonite buffer was found to have a rather complex geomechanical behavior that eventually leads to a slightly nonuniform density distribution. Nevertheless, the simulation shows that the swelling of the buffer is functioning to provide an adequate increase in confining stress on the tunnel wall, leading to a stabilization of any failure that may occur during the tunnel excavation. Finally, we describe the application of a possible approach for linking THM processes with chemistry, focusing on the evolution of primary and secondary swelling, in which the secondary swelling is caused by changes in ionic concentration, which in turn is evaluated using a transport simulation model. 相似文献
5.
The sensitivity of coal permeability to the effective stress means that changes in stress as well as pore pressure within a coal seam lead to changes in permeability. In addition coal swells with gas adsorption and shrinks with desorption; these sorption strains impact on the coal stress state and thus the permeability. Therefore the consideration of gas migration in coal requires an appreciation of the coupled geomechanical behaviour. A number of approaches to representing coal permeability incorporate the geomechanical response and have found widespread use in reservoir simulation. However these approaches are based on two simplifying assumptions; uniaxial strain (i.e. zero strain in the horizontal plane) and constant vertical stress. This paper investigates the accuracy of these assumptions for reservoir simulation of enhanced coalbed methane through CO2 sequestration. A coupled simulation approach is used where the coalbed methane simulator SIMED II is coupled with the geomechanical model FLAC3D. This model is applied to three simulation case studies assembled from information presented in the literature. Two of these are for 100% CO2 injection, while the final example is where a flue gas (12.5% CO2 and 87.5% N2) is injected. It was found that the horizontal contrast in sorption strain within the coal seam caused by spatial differences in the total gas content leads to vertical stress variation. Thus the permeability calculated from the coupled simulation and that using an existing coal permeability model, the Shi–Durucan model, are significantly different; for the region in the vicinity of the production well the coupled permeability is greater than the Shi–Durucan model. In the vicinity of the injection well the permeability is less than that calculated using the Shi–Durucan model. This response is a function of the magnitude of the strain contrast within the seam and dissipates as these contrasts diminish. 相似文献
6.
Uncertainty quantification for geomechanical and reservoir predictions is in general a computationally intensive problem, especially if a direct Monte Carlo approach with large numbers of full-physics simulations is used. A common solution to this problem, well-known for the fluid flow simulations, is the adoption of surrogate modeling approximating the physical behavior with respect to variations in uncertain parameters. The objective of this work is the quantification of such uncertainty both within geomechanical predictions and fluid-flow predictions using a specific surrogate modeling technique, which is based on a functional approach. The methodology realizes an approximation of full-physics simulated outputs that are varying in time and space when uncertainty parameters are changed, particularly important for the prediction of uncertainty in vertical displacement resulting from geomechanical modeling. The developed methodology has been applied both to a subsidence uncertainty quantification example and to a real reservoir forecast risk assessment. The surrogate quality obtained with these applications confirms that the proposed method makes it possible to perform reliable time–space varying dependent risk assessment with a low computational cost, provided the uncertainty space is low-dimensional. 相似文献
7.
Aleksander Grm Tomaž Šuštar Tomaž Rodič Franci Gabrovšek 《Mathematical Geosciences》2017,49(5):657-675
Scallops and flutes are common dissolution rock forms encountered in karst caves and surface streams. Their evolution is only partially understood and no numerical model that simulates their formation has been presented. This work at least partially fills the gap by introducing a numerical approach to simulate the evolution of different initial forms of soluble surfaces embedded in a turbulent fluid. The aim is to analyze wall dissolution phenomena from basic principles and to identify stable profiles. The analysis is based on a finite volume moving boundary method. The underlying mathematical model is a \(k-\epsilon \) turbulent model for fluid flow coupled with turbulent scalar transport. The rock wall is treated as a moving boundary, where the normal wall retreat velocity is proportional to the under-saturation of the boundary fluid cells with respect to the mineral comprising the wall. As the flow time scale is several orders of magnitude smaller than the dissolution time scale, stationary flow field, concentration field and wall propagation velocity are calculated for each iteration. The boundary at all points is then moved by distracting minimal velocity along the entire boundary from the actual velocity at a certain location, and then normalized to the maximum allowed shift, which is equal to half the height of the boundary cell. In this way only deformation of the initial wall is calculated. The method was applied to several different initial profiles. During the evolution, the profiles progressively converged towards stable forms. In this work, a framework is proposed for a computation of the moving boundary problem related to slow dissolution of a soluble surface. 相似文献
8.
The pressure variations during the production of petroleum reservoir induce stress changes in and around the reservoir. Such changes of the stress state can induce marked deformation of geological structures for stress sensitive reservoirs as chalk or unconsolidated sand reservoirs. The compaction of those reservoirs during depletion affects the pressure field and so the reservoir productivity. Therefore, the evaluation of the geomechanical effects requires to solve in a coupling way the geomechanical problem and the reservoir multiphase fluid flow problem. In this paper, we formulate the coupled geomechanical‐reservoir problem as a non‐linear fixed point problem and improve the resolution of the coupling problem by comparing in terms of robustness and convergence different algorithms. We study two accelerated algorithms which are much more robust and faster than the conventional staggered algorithm and we conclude that they should be used for the iterative resolution of coupled reservoir‐geomechanical problem. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
9.
This paper presents a new, fully-coupled, hydro-mechanical (HM) formulation for a finite-discrete element method computer code. In the newly-developed, hydraulic solver, fluid flow is assumed to occur through the same triangular mesh used for the mechanical calculations. The flow of a viscous, compressible fluid is explicitly solved based on a cubic law approximation. The implementation is verified against closed-form solutions for several flow problems. The approach is then applied to a field-scale simulation of fluid injection in a jointed, porous rock mass. Results show that the proposed method can be used to obtain unique geomechanical insights into coupled HM phenomena. 相似文献
10.
Impact of single inclined faults on the fluid flow and heat transport: results from 3-D finite element simulations 总被引：1，自引：1，他引：0
Yvonne Cherubini Mauro Cacace Guido Blöcher Magdalena Scheck-Wenderoth 《Environmental Earth Sciences》2013,70(8):3603-3618
The impact of inclined faults on the hydrothermal field is assessed by adding simplified structural settings to synthetic models. This study is innovative in carrying out numerical simulations because it integrates the real 3-D nature of flow influenced by a fault in a porous medium, thereby providing a useful tool for complex geothermal modelling. The 3-D simulations for the coupled fluid flow and heat transport processes are based on the finite element method. In the model, one geological layer is dissected by a dipping fault. Sensitivity analyses are conducted to quantify the effects of the fault’s transmissivity on the fluid flow and thermal field. Different fault models are compared with a model where no fault is present to evaluate the effect of varying fault transmissivity. The results show that faults have a significant impact on the hydrothermal field. Varying either the fault zone width or the fault permeability will result in relevant differences in the pressure, velocity and temperature field. A linear relationship between fault zone width and fluid velocity is found, indicating that velocities increase with decreasing widths. The faults act as preferential pathways for advective heat transport in case of highly transmissive faults, whereas almost no fluid may be transported through poorly transmissive faults. 相似文献
11.
Large‐scale poroelastic fractured reservoirs modeling using the fast multipole displacement discontinuity method 
下载免费PDF全文
下载免费PDF全文 An effective approach to modeling the geomechanical behavior of the network and its permeability variation is to use a poroelastic displacement discontinuity method (DDM). However, the approach becomes rather computationally intensive for an extensive system of cracks, particularly when considering coupled diffusion/deformation processes. This is because of additional unknowns and the need for time‐marching schemes for the numerical integration. The Fast Multipole Method (FMM) is a technique that can accelerate the solution of large fracture problems with linear complexity with the number of unknowns both in memory and CPU time. Previous works combining DDM and FMM for large‐scale problems have accounted only for elastic rocks, neglecting the fluid leak‐off from the fractures into the matrix and its influence on pore pressure and stress field. In this work we develop an efficient geomechanical model for large‐scale natural fracture networks in poroelastic reservoirs with fracture flow in response to injection and production operations. Accuracy and computational performance of the proposed method with those of conventional poroelastic DDM are compared through several case studies involving up to several tens of thousands of boundary elements. The results show the effectiveness of the FMM approach to successfully evaluate field‐scale problems for the design of exploitation strategies in unconventional geothermal and petroleum reservoirs. An example considering faults reveals the impact of reservoir compartmentalization because of sealing faults for both geomechanical and flow variables under elastic and poroelastic rocks. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
12.
One of the major causes of instability in geotechnical structures such as dikes or earth dams is the phenomenon of suffusion including detachment, transport and filtration of fine particles by water flow. Current methods fail to capture all these aspects. This paper suggests a new modeling approach under the framework of the porous continuous medium theory. The detachment and transport of the fine particles are described by a mass exchange model between the solid and the fluid phases. The filtration is incorporated to simulate the filling of the inter-grain voids created by the migration of the fluidized fine particles with the seepage flow, and thus, the self-filtration is coupled with the erosion process. The model is solved numerically using a finite difference method restricted to one-dimensional (1-D) flows normal to the free surface. The applicability of the model to capture the main features of both erosion and filtration during the suffusion process has been validated by simulating 1-D internal erosion tests and by comparing the numerical with the experimental results. Furthermore, the influence of the coupling between erosion and filtration has been highlighted, including the development of material heterogeneity induced by the combination of erosion and filtration.
相似文献13.
A fully coupled formulation of a hydro‐thermo‐poro‐mechanical model for a three‐phase black oil reservoir model is presented. The model is based upon the approach proposed by one of the authors which fully couples geomechanical effects to multiphase flow. Their work is extended here to include non‐isothermal effects. The gas phase contribution to the energy equation has been neglected based on a set of assumptions. The coupled formulation given herein differs in several ways when compared to the earlier work and an attempt is made to link the flow based formulation and mixture theory. The Finite Element Method is employed for the numerical treatment and essential algorithmic implementation is discussed. Numerical examples are presented to provide further understanding of the current methodology. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
14.
Reservoir simulators typically use cell‐centered finite volume schemes and do not model directly the coupling of the flow processes with the geomechanics. Coupling of geomechanics with fluid flow can be important in many cases, but introducing fully coupled geomechanical effects in those simulators is not a trivial issue, because the geomechanics is better done by using the Galerkin vertex‐centered finite element methods by which the solid displacements are computed at the vertices of the cells. This creates difficulties in interfacing cell variables with nodal variables. Uncoupled or loosely coupled models are used by many researchers/practitioners by which a reservoir model is coupled to a geomechanical model by staggering in‐time flow and deformation via a sophisticated interface that repeatedly calls first flow and then mechanics. The method therefore requires projection of the reservoir cell variables onto the nodes of the geomechanics Galerkin finite element mesh. In this note, we attempt to quantify the errors associated with cell to node projection operations. For that purpose, we use a simple model of the pressure equation for a heterogeneous medium in one dimension. We are able to derive the exact analytical solution for this problem for both nodal and cell pressures. This allows us to compute the errors due to projection analytically, function of meshing refinement and permeability field variations. We compute upper and lower bounds for the errors, and analyze their magnitude for a variety of cases. We conclude that, in general, cell to node projection operations lead to substantial errors. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
15.
构造成矿非线性动力学:1.递增应力流变学模型 总被引：1，自引：1，他引：1
本文提出了构造成矿作用动力学研究的递增应力流变学方法。利用该方法将构造变形、应力、流体流动、地球化学反应及成岩成矿作用等多过程耦合起来，可以从多种地质过程的耦合与反馈作用对构造成矿的动力学演化过程进行１～３维数值模拟。模拟的主要内容是在各种过程耦合作用下，以下描述构造成矿体系的主要变量的时空演化：（１）与成矿流体的形成和性质有关的变量，如地层中矿物（包括成矿物质）的溶解速率、流体中各组分的浓度与饱和度、流体温度、压力、离子强度等；（２）与构造变形和流体运移有关的各变量，如应力与变形速率、岩石孔隙度、构造（断裂）渗透率等；（３）与沉淀成矿有关的变量，如矿物（金瞩矿物和脉石矿物）的成核速率、各矿物的沉淀量等；（４）上述各有关变量间的时空耦合关系，如断裂渗透率时空演化与流体流动、汇聚和成矿的耦合关系等。 相似文献
16.
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. 相似文献
17.
18.
19.
Multiphysics hillslope processes triggering landslides 总被引：4，自引：3，他引：1
In 1996, a portion of a highly instrumented experimental catchment in the Oregon coast range failed as a large debris flow from heavy rain. For the first time, we quantify the 3-D multiphysical aspects that triggered this event, including the coupled sediment deformation-fluid flow processes responsible for mobilizing the slope failure. Our analysis is based on a hydromechanical continuum model that accounts for the loss of sediment strength due to increased saturation as well as the frictional drag exerted by the moving fluid. Our studies highlight the dominant role that bedrock topography and rainfall history played in defining the failure mechanism, as indicated by the location of the scarp zone that was accurately predicted by our 3-D continuum model. 相似文献
20.
This paper proposes a three-dimensional coupled hydrothermal model for fractured rock based on the finite-discrete element method to simulate fluid flow and heat transport. The 3D coupled hydrothermal model is composed of three main parts: a heat conduction model for the rock matrix, a heat transfer model for the fluid in the fractures (including heat conduction and heat convection), and a heat exchange model between the rock matrix and the fluid in the fractures. Four examples with analytical solutions are provided to verify the model. A heat exchange experiment of circulating water in a cylindrical granite sample with one fracture is simulated. The simulation results agree well with the experimental results. The effects of the fracture aperture, fluid viscosity, and pressure difference on the heat exchange between the fluid and rock are studied. Finally, an application concerned with heat transport and fluid flow in fractured rock is presented. The simulation results indicate that the 3D fully coupled hydrothermal model can capture the fluid flow and temperature evolution of rocks and fluids. 相似文献