首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 781 毫秒
1.
Hydraulic fracturing (HF) of underground formations has widely been used in different fields of engineering. Despite the technological advances in techniques of in situ HF, the industry uses semi‐analytical tools to design HF treatment. This is due to the complex interaction among various mechanisms involved in this process, so that for thorough simulations of HF operations a fully coupled numerical model is required. In this study, using element‐free Galerkin (EFG) mesh‐less method, a new formulation for numerical modeling of hydraulic fracture propagation in porous media is developed. This numerical approach, which is based on the simultaneous solution of equilibrium and continuity equations, considers the hydro‐mechanical coupling between the crack and its surrounding porous medium. Therefore, the developed EFG model is capable of simulating fluid leak‐off and fluid lag phenomena. To create the discrete equation system, the Galerkin technique is applied, and the essential boundary conditions are imposed via penalty method. Then, the resultant constrained integral equations are discretized in space using EFG shape functions. For temporal discretization, a fully implicit scheme is employed. The final set of algebraic equations that forms a non‐linear equation system is solved using the direct iterative procedure. Modeling of cracks is performed on the basis of linear elastic fracture mechanics, and for this purpose, the so‐called diffraction method is employed. For verification of the model, a number of problems are solved. According to the obtained results, the developed EFG computer program can successfully be applied for simulating the complex process of hydraulic fracture propagation in porous media. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

2.
The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr–Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils.  相似文献   

3.
Geological environments, such as petroleum reservoirs, normally exhibit physical discontinuities, for example, fractures and faults. Because of the reduced thickness of these discontinuities, finite element formulations with strong discontinuity have been applied to the numerical modelling of geological environments. Until now, two relevant characteristics of petroleum reservoirs have not been addressed by these formulations. The first is the pore pressure jump in the direction normal to a discontinuity in a fluid‐mechanical coupling condition, which is present primarily in sealing faults owing to the contrast of permeability with the porous medium. The absence of this jump can affect the prediction of the deformability of a physical discontinuity. Furthermore, reservoir models frequently use coarse meshes. Thus, the method used to evaluate the pore pressure in the discontinuity may exhibit a strong dependence relative to the mesh refinement. Based on these characteristics, in this study, a formulation of an enriched finite element for application to coupled fluid‐mechanical problems with pre‐existing physical discontinuities saturated by a single fluid is presented. The formulation employs discontinuous interpolation functions and enables the reproduction of jumps of displacement and pore pressure associated with a discontinuity inside the element without the need to discretise it. An approximation to estimate the pore pressure in the discontinuity was developed, one which seeks to minimise the influence of refinement. The element's response is verified by comparison with a one‐dimensional analytical solution and simple examples that are simulated using commercial software. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

4.
To improve the stability and efficiency of explicit technique, one proposed method is to use an unconditionally stable alternating direction explicit (ADE) scheme. However, the standard ADE scheme is only moderately accurate and restricted to uniform grids. This paper derives a novel high‐order ADE scheme capable of solving the fluid diffusion equation in non‐uniform grids. The new scheme is derived by performing a fourth‐order finite difference approximation to the spatial derivatives of the diffusion equation in non‐uniform grid. The implicit Crank‐Nicolson technique is then applied to the resulting approximation, and the subsequent equation is split into two alternating direction sweeps, giving rise to a new high‐order ADE scheme. Because the new scheme can be potentially applied in coupled hydro‐mechanical (H‐M) simulation, the pore pressure solutions from the new scheme are then sequentially coupled with an existing geomechanical simulator in the computer program Fast Lagrangian Analysis of Continua. This coupling procedure is called the sequentially explicit coupling technique based on the fourth‐order ADE scheme (SEA‐4). Verifications of well‐known consolidation problems showed that the new ADE scheme and SEA‐4 can reduce computer runtime by 46% to 75% to that of Fast Lagrangian Analysis of Continua's basic scheme. At the same time, the techniques still maintained average percentage error of 1.6% to 3.5% for pore pressure and 0.2% to 1.5% for displacement solutions and were still accurate under typical grid non‐uniformities. This result suggests that the new high‐order ADE scheme can provide an efficient explicit technique for solving the flow equation of a coupled H‐M problem, which will be beneficial for large‐scale and long‐term H‐M problems in geoengineering.  相似文献   

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

6.
In a companion Part I of this paper (Int. J. Numer. Anal. Meth. Geomech. 2008; DOI: 10.1002/nag.735 ), a coupled hydro‐mechanical (HM) formulation for geomaterials with discontinuities based on the finite element method (FEM) with double‐node, zero‐thickness interface elements was developed and presented. This Part II paper includes the numerical solution of basic practical problems using both the staggered and the fully coupled approaches. A first group of simulations, based on the classical consolidation problem with an added vertical discontinuity, is used to compare both the approaches in terms of accuracy and convergence. The monolithic or fully coupled scheme is also used in an application example studying the influence of a horizontal joint in the performance of a reservoir subject to fluid extraction. Results include a comparison with other numerical solutions from the literature and a sensitivity analysis of the mechanical parameters of the discontinuity. Some simulations are also run using both a full non‐symmetric and a simplified symmetric Jacobian matrix. On top of verifying the model developed and its capability to reflect the conductivity changes of the interface with aperture changes, the results presented also lead to interesting observations of the numerical performance of the methods implemented. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

7.
In a previous publication we developed a new mixed least‐squares method for poro‐elasticity. The approximate solution was obtained via a minimization of a least‐squares functional, based upon the equations of equilibrium, the equations of continuity and weak forms of the constitutive relationships for elasticity and Darcy flow. The formulation involved four independent types of variables: displacements, stresses, pore pressures and velocities. All of them were approximated by linear continuous triangles. Encouraged by the computational results, obtained from the two‐dimensional implementation of the method, we extended our formulation to three dimensions. In this paper we present numerical examples for the performance of continuous linear tetrahedra within the context of the mixed least‐squares method. The initial results suggest that the method works well in the nearly and entirely incompressible limits for elasticity. For poro‐elasticity, the obtained pore pressures are stable without exhibiting the oscillations, which are observed when the standard Galerkin formulation is used. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

8.
This paper presents a numerical model for the elasto‐plastic electro‐osmosis consolidation of unsaturated clays experiencing large strains, by considering electro‐osmosis and hydro‐mechanical flows in a deformable multiphase porous medium. The coupled governing equations involving the pore water flow, pore gas flow, electric flow and mechanical deformation in unsaturated clays are derived within the framework of averaging theory and solved numerically using finite elements. The displacements of the solid phase, the pressure of the water phase, the pressure of the gas phase and the electric potential are taken as the primary unknowns in the proposed model. The nonlinear variation of transport parameters during electro‐osmosis consolidation are incorporated into the model using empirical expressions that strongly depend on the degree of water saturation, whereas the Barcelona Basic Model is employed to simulate the elasto‐plastic mechanical behaviour of unsaturated clays. The accuracy of the proposed model is evaluated by validating it against two well‐known numerical examples, involving electro‐osmosis and unsaturated soil behaviour respectively. Two further examples are then investigated to study the capability of the computational algorithm in modelling multiphase flow in electro‐osmosis consolidation. Finally, the effects of gas generation at the anode, the deformation characteristics, the degree of saturation and the time dependent evolution of the excess pore pressure are discussed. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

9.
A new mixed displacement‐pressure element for solving solid–pore fluid interaction problems is presented. In the resulting coupled system of equations, the balance of momentum equation remains unaltered, while the mass balance equation for the pore fluid is stabilized with the inclusion of higher‐order terms multiplied by arbitrary dimensions in space, following the finite calculus (FIC) procedure. The stabilized FIC‐FEM formulation can be applied to any kind of interpolation for the displacements and the pressure, but in this work, we have used linear elements of equal order interpolation for both set of unknowns. Examples in 2D and 3D are presented to illustrate the accuracy of the stabilized formulation for solid–pore fluid interaction problems. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

10.
A hierarchical mathematical model for analyses of coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes at high temperature is presented. The concretes are modelled as unsaturated deforming reactive porous media filled with two immiscible pore fluids, i.e. the gas mixture and the liquid mixture, in immiscible–miscible levels. The thermo‐induced desalination process is particularly integrated into the model. The chemical effects of both the desalination and the dehydration processes on the material damage and the degradation of the material strength are taken into account. The mathematical model consists of a set of coupled, partial differential equations governing the mass balance of the dry air, the mass balance of the water species, the mass balance of the matrix components dissolved in the liquid phases, the enthalpy (energy) balance and momentum balance of the whole medium mixture. The governing equations, the state equations for the model and the constitutive laws used in the model are given. A mixed weak form for the finite element solution procedure is formulated for the numerical simulation of chemo‐thermo‐hygro‐mechanical behaviours. Special considerations are given to spatial discretization of hyperbolic equation with non‐self‐adjoint operator nature. Numerical results demonstrate the performance and the effectiveness of the proposed model and its numerical procedure in reproducing coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes subjected to fire and thermal radiation. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

11.
This paper presents a coupled hydro‐mechanical formulation for the simulation of non‐planar three‐dimensional hydraulic fractures. Deformation in the rock is modeled using linear elasticity, and the lubrication theory is adopted for the fluid flow in the fracture. The governing equations of the fluid flow and elasticity and the subsequent discretization are fully coupled. A Generalized/eXtended Finite Element Method (G/XFEM) is adopted for the discretization of the coupled system of equations. A Newton–Raphson method is used to solve the resulting system of nonlinear equations. A discretization strategy for the fluid flow problem on non‐planar three‐dimensional surfaces and a computationally efficient strategy for handling time integration combined with mesh adaptivity are also presented. Several three‐dimensional numerical verification examples are solved. The examples illustrate the generality and accuracy of the proposed coupled formulation and discretization strategies. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

12.
In this paper, a fully coupled model is developed for numerical modeling of hydraulic fracturing in partially saturated weak porous formations using the extended finite element method, which provides an effective means to simulate the coupled hydro‐mechanical processes occurring during hydraulic fracturing. The developed model is for short fractures where plane strain assumptions are valid. The propagation of the hydraulic fracture is governed by the cohesive crack model, which accounts for crack closure and reopening. The developed model allows for fluid flow within the open part of the crack and crack face contact resulting from fracture closure. To prevent the unphysical crack face interpenetration during the closing mode, the crack face contact or self‐contact condition is enforced using the penalty method. Along the open part of the crack, the leakage flux through the crack faces is obtained directly as a part of the solution without introducing any simplifying assumption. If the crack undergoes the closing mode, zero leakage flux condition is imposed along the contact zone. An application of the developed model is shown in numerical modeling of pump‐in/shut‐in test. It is illustrated that the developed model is able to capture the salient features bottomhole pressure/time records exhibit and can extract the confining stress perpendicular to the direction of the hydraulic fracture propagation from the fracture closure pressure. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

13.
Artificial ground freezing (AGF) is a commonly used technique in geotechnical engineering for ground improvement such as ground water control and temporary excavation support during tunnel construction in soft soils. The main potential problem connected with this technique is that it may produce heave and settlement at the ground surface, which may cause damage to the surface infrastructure. Additionally, the freezing process and the energy needed to obtain a stable frozen ground may be significantly influenced by seepage flow. Evidently, safe design and execution of AGF require a reliable prediction of the coupled thermo‐hydro‐mechanical behavior of freezing soils. With the theory of poromechanics, a three‐phase finite element soil model is proposed, considering solid particles, liquid water, and crystal ice as separate phases and mixture temperature, liquid pressure, and solid displacement as the primary field variables. In addition to the volume expansion of water transforming into ice, the contribution of the micro‐cryo‐suction mechanism to the frost heave phenomenon is described in the model using the theory of premelting dynamics. Through fundamental physical laws and corresponding state relations, the model captures various couplings among the phase transition, the liquid transport within the pore space, and the accompanying mechanical deformation. The verification and validation of the model are accomplished by means of selected analyses. An application example is related to AGF during tunnel excavation, investigating the influence of seepage flow on the freezing process and the time required to establish a closed supporting frozen arch. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

14.
Thermo‐hydro‐mechanical responses around a cylindrical cavity drilled or excavated in a low‐permeability formation are studied when the cavity is subjected to a time‐dependent thermal loading. The cavity is considered backfilled after it is supported by casing or lining. Solutions of temperature, pore water pressure, stress, and displacement responses are analytically formulated based on Biot's consolidation theory with the assumption that the backfilling material, supporting material, and surrounding low‐permeability formation are poroelastic media. The solution is expressed in Laplace space, and numerical inversion techniques are used to find field variables in the real‐time domain. After the solution is verified with the numerical results, it is applied in a large‐scale in situ heating test – PRACLAY heating test – for a predictive reference calculation and an extensive parametric study. Another medium‐scale in situ heating test – ATLAS III heating test – is also analyzed using the solution, which provides reasonable agreement with measurements. The new analytical solution proves to be a convenient tool for a good understanding of the resulting coupled thermo‐hydro‐mechanical behavior and is therefore valuable for the interpretation of measured data in engineering practices and for a rational design of potential radioactive waste repositories. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

15.
In this paper, a series of multimaterial benchmark problems in saturated and partially saturated two‐phase and three‐phase deforming porous media are addressed. To solve the process of fluid flow in partially saturated porous media, a fully coupled three‐phase formulation is developed on the basis of available experimental relations for updating saturation and permeabilities during the analysis. The well‐known element free Galerkin mesh‐free method is adopted. The partition of unity property of MLS shape functions allows for the field variables to be extrinsically enriched by appropriate functions that introduce existing discontinuities in the solution field. Enrichment of the main unknowns including solid displacement, water phase pressure, and gas phase pressure are accounted for, and a suitable enrichment strategy for different discontinuity types are discussed. In the case of weak discontinuity, the enrichment technique previously used by Krongauz and Belytschko [Int. J. Numer. Meth. Engng., 1998; 41:1215–1233] is selected. As these functions possess discontinuity in their first derivatives, they can be used for modeling material interfaces, generating only minor oscillations in derivative fields (strain and pressure gradients for multiphase porous media), as opposed to unenriched and constrained mesh‐free methods. Different problems of multimaterial poro‐elasticity including fully saturated, partially saturated one, and two‐phase flows under the assumption of fully coupled extended formulation of Biot are examined. As a further development, problems involved with both material interface and impermeable discontinuities, where no fluid exchange is permitted across the discontinuity, are considered and numerically discussed. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

16.
In recent years, the authors have proposed a new double‐node zero‐thickness interface element for diffusion analysis via the finite element method (FEM) (Int. J. Numer. Anal. Meth. Geomech. 2004; 28 (9): 947–962). In the present paper, that formulation is combined with an existing mechanical formulation in order to obtain a fully coupled hydro‐mechanical (or HM) model applicable to fractured/fracturing geomaterials. Each element (continuum or interface) is formulated in terms of the displacements (u) and the fluid pressure (p) at the nodes. After assembly, a particular expression of the traditional ‘u–p’ system of coupled equations is obtained, which is highly non‐linear due to the strong dependence between the permeability and the aperture of discontinuities. The formulation is valid for both pre‐existing and developing discontinuities by using the appropriate constitutive model that relates effective stresses to relative displacements in the interface. The system of coupled equations is solved following two different numerical approaches: staggered and fully coupled. In the latter, the Newton–Raphson method is used, and it is shown that the Jacobian matrix becomes non‐symmetric due to the dependence of the discontinuity permeability on the aperture. In the part II companion paper (Int. J. Numer. Anal. Meth. Geomech. 2008; DOI: 10.1002/nag.730 ), the formulation proposed is verified and illustrated with some application examples. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

17.
The development of stability problems related to classical mixed methods has recently been observed. In this study, a soil-water coupled boundary-value problem, one type of stability problem, is presented using the element-free Galerkin method (EFG method). In this soil-water coupled problem, anomalous behavior appears in the pressure field unless a stabilization technique is used. The remedy to such numerical instability has generally been to adopt a higher interpolation order for the displacements than for the pore pressure. As an alternative, however, an added stabilization term is incorporated into the equilibrium equation. The advantages of this stabilization procedure are as follows: (1) The interpolation order for the pore pressure is the same as that for the displacements. Therefore, the interpolation functions in the pore pressure field do not reduce the accuracy of the numerical results. (2) The stabilization term consists of first derivatives. The first derivatives of the interpolation functions for the EFG Method are smooth, and therefore, the solutions for pore pressure are accurate. In order to validate the above stabilization technique, some numerical results are given. It can be seen from the results that a good convergence is obtained with this stabilization term.  相似文献   

18.
This paper presents an algorithm and a fully coupled hydromechanical‐fracture formulation for the simulation of three‐dimensional nonplanar hydraulic fracture propagation. The propagation algorithm automatically estimates the magnitude of time steps such that a regularized form of Irwin's criterion is satisfied along the predicted 3‐D fracture front at every fracture propagation step. A generalized finite element method is used for the discretization of elasticity equations governing the deformation of the rock, and a finite element method is adopted for the solution of the fluid flow equation on the basis of Poiseuille's cubic law. Adaptive mesh refinement is used for discretization error control, leading to significantly fewer degrees of freedom than available nonadaptive methods. An efficient computational scheme to handle nonlinear time‐dependent problems with adaptive mesh refinement is presented. Explicit fracture surface representations are used to avoid mapping of 3‐D solutions between generalized finite element method meshes. Examples demonstrating the accuracy, robustness, and computational efficiency of the proposed formulation, regularized Irwin's criterion, and propagation algorithm are presented.  相似文献   

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

20.
This paper is interested in the hydro‐mechanical behaviour of an underground cavity abandoned at the end of its service life. It is an extension of a previous study that accounted for a poro‐elastic behaviour of the rock mass (Int. J. Comput. Geomech. 2007; DOI: 10.1016/j.compgeo.2007.11.003 ). Deterioration of the lining support with time leads to the transfer of the loading from the exterior massif to the interior backfill. The in situ material has a poro‐visco‐elastic constitutive behaviour while the backfill is poro‐elastic, both saturated with water. This loading transfer is accompanied by an inward cavity convergence, thereby compressing the backfill, and induces an outward water flow. This leads to a complex space–time evolution of pore pressures, displacements and stresses, which is not always intuitive. In its general setting, a semi‐explicit solution to this problem is developed, using Laplace transform, the inversion being performed numerically. Analytical inversion leading to a quasi‐explicit solution in the time domain is possible by identifying the characteristic creep and relaxation times of volumetric strains with those of the deviatoric strains, on the basis of a parametric study. A few numerical examples are given to illustrate the hydro‐mechanical behaviour of the cavity and highlight the influence of key parameters (e.g. stiffness of backfill, lining deterioration rate, etc.). Further studies accounting for more general material behaviours for the backfill and external ground are ongoing. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

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

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