共查询到20条相似文献，搜索用时 46 毫秒
1.
Fluid injection–induced tensile opening is modeled using an extended finite volume method (XFVM). An embedded fracture strategy is used for the flow problem, that is, the fractures are discretized using finite volume segments without resolving the grid around them. Further, the discontinuities across fractures are modeled using special basis functions. The fracture openings due to enhanced fluid pressure and the associated shear slip due to traction free boundary condition on the fracture segments are both modeled using these special discontinuity basis functions. Mass transfer between fractures and matrix is modeled using the pressure difference. The enhancement of fracture storativity due to tensile opening leads to stronger coupling between flow and mechanics. An iterative scheme relying on the fixed-stress approach for fractures, which conserves the stress dependent terms over each iteration of the flow problem, has been introduced. Tensile opening has been simulated for single fractures embedded in two- and three-dimensional matrices. The convergence criterion for sequentially implicit fixed-stress scheme for fractures embedded in elastic media is established and has been validated numerically. Further, for 2D simulations, the effect of the matrix permeability for fracture propagation due to tensile opening has been studied. 相似文献
2.
In this paper, a numerical model to predict flow‐induced shear failure along pre‐existing fractures is presented. The framework is based on a discrete fracture representation embedded in a continuum describing the damaged matrix. A finite volume method is used to compute both flow and mechanical equilibrium, whereas specifically tailored basis functions are used to account for the physics at discontinuities. The failure criterion is based on a maximum shear strength limit, which changes with varying compressive stress on the fracture manifold. The displacements along fracture manifolds are obtained such that force balance is achieved under conditions, where shear stress of the failing fracture segment is constrained to the maximum shear strength at the segment. Simultaneously, the fluid pressure is computed independently of the shear slip. A relaxation model approach is used to obtain the maximum shear limit on the fracture manifold, which leads to grid convergence. 相似文献
3.
Fushen Liu 《国际地质力学数值与分析法杂志》2020,44(12):1634-1655
The paper presents an embedded strong discontinuity approach to simulate single hydraulic fracture propagation in the poroelastic medium under plane-strain conditions. The method enriches the strain field with the discontinuous deformation mode and allows the fracture to be modeled inside elements. The Mode-I fracture initiation and propagation are described by the trilinear cohesive law, which is implemented by the penalty method. The enhanced permeability inside the fractured elements is dependent on the fracture aperture. Hydraulic fracture propagation is driven by the high pressure gradient near the fracture. Fluid transfer between the fracture and bulk rock is automatically captured within the poroelastic framework. The numerical framework is verified by the comparisons with the asymptotic analytical solutions for single hydraulic fracture propagation. 相似文献
4.
The fluid flow in rock fractures during shear processes has been an important issue in rock mechanics and is investigated in this paper using finite element method (FEM), considering evolutions of aperture and transmissivity with shear displacement histories under different normal stress and normal stiffness conditions as measured during laboratory coupled shear-flow tests. The distributions of fracture aperture and its evolution during shearing were calculated from the initial aperture, based on the laser-scanned sample surface roughness results, and shear dilations measured in the laboratory tests. Three normal loading conditions were adopted in the tests: simple normal stress and mixed normal stress and normal stiffness to reflect more realistic in situ conditions. A special algorithm for treatment of the contact areas as zero-aperture elements was used to produce more accurate flow field simulations, which is important for continued simulations of particle transport but often not properly treated in literature. The simulation results agree well with the measured hydraulic apertures and flow rate data obtained from the laboratory tests, showing that complex histories of fracture aperture and tortuous flow fields with changing normal loading conditions and increasing shear displacements. With the new algorithm for contact areas, the tortuous flow fields and channeling effects under normal stress/stiffness conditions during shearing were more realistically captured, which is not possible if traditional techniques by assuming very small aperture values for the contact areas were used. These findings have an important impact on the interpretation of the results of coupled hydro-mechanical experiments of rock fractures, and on more realistic simulations of particle transport processes in fractured rocks. 相似文献
5.
Benjamin Ganis Mark E. Mear A. Sakhaee-Pour Mary F. Wheeler Thomas Wick 《Computational Geosciences》2014,18(5):613-624
We describe an algorithm for modeling saturated fractures in a poroelastic domain in which the reservoir simulator is coupled with a boundary element method. A fixed stress splitting is used on the underlying fractured Biot system to iteratively couple fluid and solid mechanics systems. The fluid system consists of Darcy’s law in the reservoir and is computed with a multipoint flux mixed finite element method, and a Reynolds’ lubrication equation in the fracture solved with a mimetic finite difference method. The mechanics system consists of linear elasticity in the reservoir and is computed with a continuous Galerkin method, and linear elasticity in the fracture is solved with a weakly singular symmetric Galerkin boundary element method. This algorithm is able to compute both unknown fracture width and unknown fluid leakage rate. An interesting numerical example is presented with an injection well inside of a circular fracture. 相似文献
6.
Discrete fracture model for simulating waterflooding processes under fracturing conditions 
下载免费PDF全文
下载免费PDF全文 In our study, we develop a model for simulating fracturing processes in a poroelastic medium. The proposed approach combines the discrete fracture model enriched with contact plane mechanics. The model captures mechanical interactions of fractures and a deformable medium, fluid, and heat transfer in fractures and in a porous medium. Both effects of poroelasticity and thermoelasticity are accounted in our model. Mass and heat conservation equations are approximated by the finite volume method, and mechanical equilibrium equations are discretized by means of the Galerkin finite element approach. Two‐dimensional grid facets between 3‐dimensional finite elements are considered as possible fracture surfaces. Most of these facets are inactive from the beginning and are activated throughout the simulation. A fracture propagation criterion, based on Irwin's approach, is verified on each nonlinear iteration. When the criterion is satisfied, additional contact elements are added into finite element and discrete fracture model formulations respectively. The proposed approach allows modeling of existing natural and artificially created fractures within one framework. The model is tested on single‐ and multiple‐phase fluid flow examples for both isothermal and thermal conditions and verified against existing semianalytical solutions. The applicability of the approach is demonstrated on an example of practical interests where a sector model of an oil reservoir is simulated with different injection and production regimes. 相似文献
7.
Poro‐mechanical and thermo‐mechanical processes change the fracture aperture and thus affect the water flow pattern in the fracture during the cold water injection into enhanced geothermal systems (EGS). In addition, the stresses generated by these processes contribute to the phenomenon of reservoir seismicity. In this paper, we present a three‐dimensional (3D) partially coupled poro‐thermoelastic model to investigate the poroelastic and thermoelastic effects of cold water injection in EGS. In the model, the lubrication fluid flow and the convective heat transfer in the fracture are modeled by the finite element method, while the pore fluid diffusion and heat conductive transfer in the reservoir matrix are assumed to be 3D and modeled by the boundary integral equation method without the need to discretize the reservoir. The stresses at the fracture surface and in the reservoir matrix are obtained from the numerical model and can be used to assess the variation of in situ stress and induced seismicty with injection/extraction. Application of the model shows that rock cooling induces large tensile stresses and increases fracture conductivity, whereas the rock dilation caused by fluid leakoff decreases fracture aperture and increases compressive total stresses around the injection zone. However, increases in pore pressure reduce the effective stresses and can contribute to rock failure, fracture slip, and microseismic activity. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
8.
The propagation of hydraulic fracture in elastic rocks has widely been investigated. In the paper, we shall focus on numerical modeling of hydraulic fracturing in a class of porous rocks exhibiting plastic deformation. The plastic strain of porous rocks is described by a non-associated plastic model based on Drucker–Prager criterion. The plastic deformation is coupled with fluid pressure evolution described by the lubrication theory. An extended finite element method is used for modeling the propagation of fracture. The fracture propagation criterion is based on the J-integral. The proposed numerical model is validated by comparisons with numerical and analytical results. The influence of plastic deformation on fracture propagation process is investigated.
相似文献9.
Magnus Wangen 《Computational Geosciences》2013,17(4):647-659
A procedure based on the finite element method is suggested for modeling of 3D hydraulic fracturing in the subsurface. The proposed formulation partitions the stress field into the initial stress state and an additional stress state caused by pressure buildup. The additional stress is obtained as a solution of the Biot equations for coupled fluid flow and deformations in the rock. The fluid flow in the fracture is represented on a regular finite element grid by means of “fracture” porosity, which is the volume fraction of the fracture. The use of the fracture porosity allows for a uniform finite element formulation for the fracture and the rock, both with respect to fluid pressure and displacement. It is demonstrated how the fracture aperture is obtained from the displacement field. The model has a fracture criterion by means of a strain limit in each element. It is shown how this criterion scales with the element size. Fracturing becomes an intermittent process, and each event is followed by a pressure drop. A procedure is suggested for the computation of the pressure drop. Two examples of hydraulic fracturing are given, when the pressure buildup is from fluid injection by a well. One case is of a homogeneous rock, and the other case is an inhomogeneous rock. The fracture geometry, well pressure, new fracture area, and elastic energy released in each event are computed. The fracture geometry is three orthogonal fracture planes in the homogeneous case, and it is a branched fracture in the inhomogeneous case. 相似文献
10.
Fractured-rock aquifers display spatially and temporally variable hydraulic conductivity generally attributed to variable fracture intensity and connectivity. Empirical evidence suggests fracture aperture and hydraulic conductivity are sensitive to in situ stress. This study investigates the sensitivity of fractured-rock hydraulic conductivity, groundwater flow paths, and advection-dominated transport to variable shear and normal fracture stiffness magnitudes for a range of deviatoric stress states. Fracture aperture and hydraulic conductivity are solved for analytically using empirical hydromechanical coupling equations; groundwater flow paths and ages are then solved for numerically using groundwater flow and advection-dispersion equations in a traditional Toth basin. Results suggest hydraulic conductivity alteration is dominated by fracture normal closure, resulting in decreasing hydraulic conductivity and increasing groundwater age with depth, and decreased depth of long flow paths with decreasing normal stiffness. Shear dilation has minimal effect on hydraulic conductivity alteration for stress states investigated here. Results are interpreted to suggest that fracture normal stiffness influences hydraulic conductivity of hydraulically active fractures and, thus, affects flow and transport in shallow (<1 km) fractured-rock aquifers. It is suggested that observed depth-dependent hydraulic conductivity trends in fractured-rock aquifers throughout the world may be partly a manifestation of hydromechanical phenomena. 相似文献
11.
Numerical modeling of hydraulic fracture propagation,closure and reopening using XFEM with application to in‐situ stress estimation 下载免费PDF全文
下载免费PDF全文 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. 相似文献
12.
岩溶地区地下发育着大量的溶洞和地下河管道,地下水流状态既有层流也有紊流,而紊流是溶洞管道形成的重要条件。紊流的形成受到岩石初始裂隙的影响,初始裂隙的张开度、分布、走向、迹长、密度等因素都影响着裂隙发育过程中水流状态的变化。通过对不同统计特征的初始裂隙网络进行水流和溶蚀的数值模拟发现,以张开度标准差反映的裂隙网络非均匀性越强,模拟紊流出现的时间就越早;主要裂隙的存在使裂隙网络的非均性增强,主要裂隙与水力梯度总方向的角度越小,紊流出现的时间就越早;当裂隙平均迹长过小时会导致裂隙连通性较差,影响裂隙水流和溶蚀作用;裂隙密度,尤其是主要裂隙密度,对岩溶发育的影响较大。相对于次要裂隙,如果主要裂隙密度偏小,紊流形成时间会大大增加,甚至很难形成紊流。当初始裂隙张开度小于0.001 cm,增大水力梯度仍没有紊流发生,岩溶几乎不发育。 相似文献
13.
危岩是三峡库区典型的地质灾害类型之一,而主控结构面受荷断裂扩展是危岩发育成灾的关键核心。将危岩主控结构面类比为宏观裂纹,利用扩展有限元法在模拟裂纹扩展方面的优势,基于考虑裂纹面水压力作用的虚功原理推导出了采用扩展有限元法分析水力劈裂问题的控制方程,给出了危岩主控结构面水力劈裂问题的扩展有限元实现方法,对重庆万州太白岩危岩主控结构面的水力劈裂进行了数值模拟分析。计算结果表明：暴雨是威胁危岩稳定性的最敏感因素,随着裂隙水压力上升,裂端拉应力会急剧升高,危岩的稳定性降低;I型裂纹扩展是危岩主要的结构面扩展形式,结构面一旦发生开裂,将处于非稳定扩展状态。 相似文献
14.
Liao Jianxing Gou Yang Feng Wentao Mehmood Faisal Xie Yachen Hou Zhengmeng 《Acta Geotechnica》2020,15(2):279-295
Although hydraulic fracturing has been massively studied and applied as a key technique to enhance the gas production from tight formations, some problems and uncertainties exist to accurately predict and analyze the fracture behavior in complex reservoirs, especially in the naturally fractured reservoirs like shale reservoirs. This paper presents a full 3D numerical model (FLAC3D) to study hydraulic fracturing behavior under the impact of preexisting orthogonal natural fractures. In this numerical model, the hydraulic fracture propagation direction is assumed perpendicular to the minimum principal stress and activated only by tensile failure, whereas the preexisting natural fractures can be activated by tensile or shear failure or a combination of them, and only tensile failure can open the natural fracture as well. The newly developed model was used to study the impact of preexisting orthogonal natural fractures on hydraulic fracturing behavior, based on a multistage hydraulic fracturing operation in a naturally fractured reservoir from the Barnett Shale formation, northwest of Texas in USA. In this multistage operation, two more representative stages, i.e., stage 1 with a relatively large horizontal stress anisotropy of 3.3 MPa and stage 4 with a comparatively small one of 1.3 MPa, were selected to conduct the simulation. Based on the numerical results, one can observe that the interaction between hydraulic and natural fracture is driven mainly by induced stress around fracture tip. Besides, the horizontal stress anisotropy plays a key role in opening the natural fracture. Thus, no significant opened fracture is activated on natural fracture in stage 1, while in stage 4 an opened fracture invades to about 90 m into the first natural fracture. Conversely, the hydraulic fracture length in stage 1 is much longer than in stage 4, as some fluid volume is stored in the opened natural fracture in stage 4. In this work, the shear failure on natural fractures is treated as the main factor for inducing the seismic events. And the simulated seismic events, i.e., shear failure on natural fractures, are very comparable with the measured seismic events.
相似文献15.
针对现有页岩气体积压裂层理裂缝扩展模型及其计算方法研究的不足,基于三点弯曲试验,结合数字图像法获取了页岩层理关键力学参数;利用弹性力学理论和线弹簧模型建立了页岩气体积压裂层理裂缝扩展拟三维数学模型,并得到了试验验证;开发了裂缝几何参数计算程序,计算分析了层理参数与压裂施工参数对水力裂缝扩展的影响规律。结果表明：当层理刚度小于10 GPa/m及大于30 GPa/m时,剪切滑移量达到极大值及极小值,且基本保持不变,当层理刚度在10～30 GPa/m之间时,层理刚度与剪切滑移量呈线性负相关关系;当层理密度为5～7条时,主裂缝会沟通更多的层理;当层理强度为5～8 MPa时,水力裂缝易穿层扩展,且能使层理产生剪切滑移,从而生成复杂裂缝;当压裂液排量和压裂液黏度分别在9～12 m3/min和2.5～5m Pa·s范围内时,水力裂缝易穿层扩展,最终形成十字型裂缝,有助于复杂裂缝的形成。该研究对认识页岩层理力学性能及其对层理裂缝扩展规律的影响有一定的指导意义。 相似文献
16.
Hydraulic Characteristics of Rough Fractures in Linear Flow under Normal and Shear Load 总被引：19，自引：3，他引：19
Summary
A hydro-mechanical testing system, which is capable of measuring both the flow rates and the normal and shear displacement
of a rock fracture, was built to investigate the hydraulic behaviour of rough tension fractures. Laboratory hydraulic tests
in linear flow were conducted on rough rock fractures, artificially created using a splitter under various normal and shear
loading. Prior to the tests, aperture distributions were determined by measuring the topography of upper and lower fracture
surfaces using a laser profilometer. Experimental variograms of the initial aperture distributions were classified into four
groups of geostatistical model, though the overall experimental variograms could be well fitted to the exponential model.
The permeability of the rough rock fractures decayed exponentially with respect to the normal stress increase up to 5 MPa.
Hydraulic behaviours during monotonic shear loading were significantly affected by the dilation occurring until the shear
stress reached the peak strength. With the further dilation, the permeability of the rough fracture specimens increased more.
However, beyond shear displacement of about 7 to 8 mm, permeability gradually reached a maximum threshold value. The combined
effects of both asperity degradation and gouge production, which prohibited the subsequent enlargement of mean fracture aperture,
mainly caused this phenomenon. Permeability changes during cyclic shear loading showed somewhat irregular variations, especially
after the first shear loading cycle, due to the complex interaction from asperity degradations and production of gouge materials.
The relation between hydraulic and mechanical apertures was analyzed to investigate the valid range of mechanical apertures
to be applied to the cubic law.
Received June 12, 2001; accepted February 26, 2002 Published online September 2, 2002 相似文献
17.
A procedure for numerical approximation to two-dimensional, hydraulically-driven fracture propagation in a poroelastic material is described. The method uses a partitioned solution procedurè to solve a finite element approximation to problems described by the theory of poroelasticity, in conjunction with a finite difference approximation for modelling fluid flow along the fracture. An equilibrium fracture model based on a generalized, Dugdale–Barenblatt concept is used to determine the fracture dimensions. An important feature is that the fracture length is a natural product of the solution algorithm. Two example problems verify the accuracy of the numerical procedure and a third example illustrates a fully-coupled simulation of fracture propagation. Photographs taken from a high-performance engineering workstation provide insight into the nature of the coupling among the physical phenomena. 相似文献
18.
An embedded fracture modeling framework for simulation of hydraulic fracturing and shear stimulation 总被引：1，自引：0，他引：1
Jack H. Norbeck Mark W. McClure Jonathan W. Lo Roland N. Horne 《Computational Geosciences》2016,20(1):1-18
A numerical modeling framework is described that is able to calculate the coupled processes of fluid flow, geomechanics, and rock failure for application to general engineering problems related to reservoir stimulation, including hydraulic fracturing and shear stimulation. The numerical formulation employs the use of an embedded fracture modeling approach, which provides several advantages over more traditional methods in terms of computational complexity and efficiency. Specifically, the embedded fracture modeling strategy avoids the usual requirement that the discretization of the fracture domain conforms to the discretization of the rock volume surrounding the fractures. As fluid is exchanged between the two domains, conservation of mass is guaranteed through a coupling term that appears as a simple source term in the governing mass balance equations. In this manner, as new tensile fractures nucleate and propagate subject to mechanical effects, numerical complexities associated with the introduction of new fracture control volumes are largely negated. In addition, the ability to discretize the fractures and surrounding rock volume independently provides the freedom to choose an acceptable level of discretization for each domain separately. Three numerical examples were performed to demonstrate the utility of the embedded fracture model for application to problems involving fluid flow, mechanical deformation, and rock failure. The results of the numerical examples confirm that the embedded fracture model was able to capture accurately the complex and nonlinear evolution of reservoir permeability as new fractures propagate through the reservoir and as fractures fail in shear. 相似文献
19.
Wang Mingyu The University of Arizona Tucson Arizona USA Department of Water Resources Environmental Engineering China University of Geosciences Beijing Chen Jinsong Wan Li Department of Water Resources Environmental Engineering 《中国地质大学学报(英文版)》2001,12(3)
ＩＮＴＲＯＤＵＣＴＩＯＮＧｒｏｕｎｄｗａｔｅｒｏｒｆｌｕｉｄｆｌｏｗｍｏｄｅｌｉｎｇｉｎｆｒａｃｔｕｒｅｄｒｏｃｋｓｉｓａｃｏｍｐｌｉｃａｔｅｄｔｈｅｏｒｅｔｉｃａｌａｎｄａｐｐｌｉｅｄｔｏｐｉｃ.Ｂｏｔｈｔｈｅｏｒｅｔｉｃａｌｌｙａｎｄｏｐｅｒａｔｉｏｎａｌｌｙ ,ｉｔｉｓｉｍｐｏｒｔａｎｔｉｎｍａｎｙｆｉｅｌｄｓｓｕｃｈａｓｇｅｏｌｏｇｉｃａｌａｎｄｈｙｄｒｏｇｅｏｌｏｇｉｃａｌｅｎｇｉｎｅｅｒｉｎｇ ,ｅｎｖｉｒｏｎｍｅｎｔａｌｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｅｔｒｏｌｅｕｍｅｎｇｉｎｅｅｒｉｎ… 相似文献
20.
We present a stabilized extended finite element formulation to simulate the hydraulic fracturing process in an elasto‐plastic medium. The fracture propagation process is governed by a cohesive fracture model, where a trilinear traction‐separation law is used to describe normal contact, cohesion and strength softening on the fracture face. Fluid flow inside the fracture channel is governed by the lubrication equation, and the flow rate is related to the fluid pressure gradient by the ‘cubic’ law. Fluid leak off happens only in the normal direction and is assumed to be governed by the Carter's leak‐off model. We propose a ‘local’ U‐P (displacement‐pressure) formulation to discretize the fluid‐solid coupled system, where volume shape functions are used to interpolate the fluid pressure field on the fracture face. The ‘local’ U‐P approach is compatible with the extended finite element framework, and a separate mesh is not required to describe the fluid flow. The coupled system of equations is solved iteratively by the standard Newton‐Raphson method. We identify instability issues associated with the fluid flow inside the fracture channel, and use the polynomial pressure projection method to reduce the pressure oscillations resulting from the instability. Numerical examples demonstrate that the proposed framework is effective in modeling 3D hydraulic fracture propagation. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献