首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
In this paper, 3D steady‐state fluid flow in a porous medium with a large number of intersecting fractures is derived numerically by using collocation method. Fluid flow in the matrix and fractures is described by Darcy's law and Poiseuille's law, respectively. The recent theoretical development presented a general potential solution to model the steady‐state flow in fractured porous media under a far‐field condition. This solution is a hypersingular integral equation with pressure field in the fracture surfaces as the main unknown. The numerical procedure can resolve the problem for any form of fractures and also takes into account the interactions and the intersection between fractures. Once the pressure field and then the flux field in fractures have been determined, the pressure field in the porous matrix is computed completely. The basic problem of a single fracture is investigated, and a semi‐analytical solution is presented. Using the solution obtained for a single fracture, Mori‐Tanaka and self‐consistent schemes are employed for upscaling the effective permeability of 3D fractured porous media. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

2.
A mixed finite element–boundary element solution for the analysis of two-dimensional flow in porous media composed of rock blocks and discrete fractures is described. The rock blocks are modelled implicitly by using boundary elements whereas finite elements are adopted to model the discrete fractures. The computational procedure has been implemented in a hybrid code which has been validated first by comparing the numerical results with the closed-form solution for flow in a porous aquifer intercepted by a vertical fracture only. Then, a more complex problem has been solved where a pervious, homogeneous and isotropic matrix containing a net of fractures is considered. The results obtained are shown to describe satisfactorily the main features of the flow problem under study. © 1997 by John Wiley & Sons, Ltd.  相似文献   

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.
A new artificial boundary approach for transient seepage problems in unbounded domain is presented. The artificial boundary condition at the truncated boundary is derived from the analytical solutions for transient seepage problems in one dimension, including solutions, respectively, for flow in one‐dimensional infinite space and for radial flow in an infinite layer, and then it is tentatively applied for some two dimensional problems in addition to the one‐dimensional problems mentioned above. The boundary conditions derived relate the time‐dependent boundary flux with the time derivative of the hydraulic head at the truncated boundary, which makes the implementation much easier compared with the infinite element method. The accuracy and efficiency of the artificial boundary are validated by several numerical examples, which shows that the proposed boundary can give very good results for one‐dimensional transient seepage problems, as expected, whereas reasonable results can be also obtained for two‐dimensional problems, such as two‐dimensional axisymmetric flow and flow in an infinite plane. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

5.
This paper presents an analytical layer element solution to axisymmetric thermal consolidation of multilayered porous thermoelastic media containing a deep buried heat source. By applying the Laplace–Hankel transform to the state variables involved in the basic governing equations of porous thermoelasticity, the analytical layer elements that describe the relationship between the transformed generalized stresses and displacements of a finite layer and a half‐space are derived. The global stiffness matrix equation is obtained by assembling the interrelated layer elements, and the real solutions in the physical domain are achieved by numerical inversion of the Laplace–Hankel transform after obtaining the solutions in the transformed domain. Finally, numerical calculations are performed to demonstrate the accuracy of this method and to investigate the influence of heat source's types, layering, and the porous thermoelastic material parameters on thermal consolidation behavior. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

6.
A finite element procedure is developed to accurately locate the free surface of unconfined seepage flow through porous media. The free surface is taken as the boundary between wet and dry soils, with flow in the saturated region characterized by Darcy's law. The method involves equations and meshing which are fully consistent with a general formulation for geotechnical engineering problems involving simultaneous solution of pore fluid pressures and soil skeleton displacements. Accuracy and versatility of the proposed procedure are demonstrated by solving various unconfined seepage flow problems through earth structures. Free surfaces and flownets are presented for the calculated flow fields.  相似文献   

7.
This paper presents a novel analytical solution to the transient, z‐dependent, and asymmetric problem of an infinite wellbore drilled into a fluid‐saturated porous medium. The formulations are based on Biot's linear theory of poroelasticity, in which the dependency of poroelastic field variables to spatial coordinates as well as time domain is considered in the most general form. This gives flexibility to the solution in cases that cannot be analyzed using the conventional plane strain or symmetric models. One such case is when calculating the stress variations around an inclined wellbore where the far‐field stresses are acting over a finite vertical section. The results of our solution to this case with a three‐dimensional state of far‐field stress are used to analyze the stability of inclined wellbores passing through abnormally stressed formations. The presented solution is capable of finding expressions for fundamental solutions with stress or flow boundary conditions at the wellbore. These solutions are here adopted to analyze the pressure disturbances generated by multiprobe formation tester, a standard wireline device that is designed for downhole fluid sampling as well as estimating the directional permeabilities of subsurface earth formations. A comparison with the conventional solution for the relevant pressure diffusion equation indicates that the poroelastic effect is relatively significant in relation to the transient response of the pore pressure. Further, it is shown that the finite dimensions of sink probe would, to a great extent, contribute to the formation's pore pressure variations at its immediate proximity. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

8.
The theoretical aspects of fully coupled thermohydromechanical behaviour of saturated porous media are presented. The non-linear behaviour of soil skeleton is assumed. A new concept called ‘thermal void ratio state surface’ is introduced to include thermal effects, and the stress state level influence on volume changes. The fluid phase flows according to Darcy's law and energy transport is assumed to follow Fourier's law classically. Variation of water permeability, water and solid unit weight due to thermal effects and pore pressure changes are included. A finite element package is developed based on final matrix form obtained from discretization of integral form of field equations by finite element method and integration in time. A very good agreement between the theoretical predictions and the experimental results was obtained for the several simple problems proposed by other authors. © 1997 by John Wiley & Sons, Ltd.  相似文献   

9.
This paper presents a model for the analysis of plane waves diffraction at a cavity in an infinite homogeneous poroelastic saturated medium, lined by a lining composed of four equal segments. An elastic boundary layer is placed between the cavity lining and the infinite porous medium. The boundary layer is simulated by ‘elastic boundary conditions’ in which the bulk matrix stress is proportional to the relative displacement between the lining and the surrounding medium matrix boundary. In addition, fluid impermeability through the intermediate layer is assumed. For the frequencies, that differ from the pseudoresonanse frequencies, the problem was reduced to the problem of an ideal elastic medium. A closed‐form analytical solution of the problem was obtained using Fourier–Bessel series, the convergence of which was proven. It was shown that the number of series terms required to obtain a desired level of accuracy can be determined in advance. The influence of the medium porosity on the medium dynamic stress concentration was studied. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

10.
Mathematical modelling of the ascent of free fluid through relatively strong rock, deep in the Earth's mantle, presents a challenge in geomechanics. Here the medium is considered as fluid-saturated, porous, elastic and bounded, and the fluid enters at a point source. An explicit finite difference method is developed for the numerical solution to the problem of the dilatation of a fluid-saturated porous elastic sphere due to a point fluid source of constant strength at the centre of the sphere. A cubic spline interpolant is used to evaluate a definite integral which occurs in the boundary condition for the pore fluid pressure at the surface of the sphere. The numerical solutions for the dilatation and pore fluid pressure are compared with analytical solutions and the absolute and relative errors of the numerical solutions are calculated. When the fluid source is switched on, the pore fluid pressure starts to decrease, reaches a minimum value and then steadily increases. The initial time rate of decrease of the pore fluid pressure is independent of the radial distance from the source. It decreases as the radius of the sphere increases and vanishes for a point fluid source in an infinite porous elastic medium.  相似文献   

11.
We derive a macroscopic model for single-phase, incompressible, viscous fluid flow in a porous medium with small cavities called vugs. We model the vuggy medium on the microscopic scale using Stokes equations within the vugular inclusions, Darcy's law within the porous rock, and a Beavers–Joseph–Saffman boundary condition on the interface between the two regions. We assume periodicity of the medium and obtain uniform energy estimates independent of the period. Through a two-scale homogenization limit as the period tends to zero, we obtain a macroscopic Darcy's law governing the medium on larger scales. We also develop some needed generalizations of the two-scale convergence theory needed for our bimodal medium, including a two-scale convergence result on the Darcy–Stokes interface. The macroscopic Darcy permeability is computable from the solution of a cell problem. An analytic solution to this problem in a simple geometry suggests that: (1) flow along vug channels is primarily Poiseuille with a small perturbation related to the Beavers–Joseph slip, and (2) flow that alternates from vug to matrix behaves as if the vugs have infinite permeability.  相似文献   

12.
This paper presents a fracture mapping (FM) approach combined with the extended finite element method (XFEM) to simulate coupled deformation and fluid flow in fractured porous media. Specifically, the method accurately represents the impact of discrete fractures on flow and deformation, although the individual fractures are not part of the finite element mesh. A key feature of FM‐XFEM is its ability to model discontinuities in the domain independently of the computational mesh. The proposed FM approach is a continuum‐based approach that is used to model the flow interaction between the porous matrix and existing fractures via a transfer function. Fracture geometry is defined using the level set method. Therefore, in contrast to the discrete fracture flow model, the fracture representation is not meshed along with the computational domain. Consequently, the method is able to determine the influence of fractures on fluid flow within a fractured domain without the complexity of meshing the fractures within the domain. The XFEM component of the scheme addresses the discontinuous displacement field within elements that are intersected by existing fractures. In XFEM, enrichment functions are added to the standard finite element approximation to adequately resolve discontinuous fields within the simulation domain. Numerical tests illustrate the ability of the method to adequately describe the displacement and fluid pressure fields within a fractured domain at significantly less computational expense than explicitly resolving the fracture within the finite element mesh. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

13.
A high‐frequency open boundary has been developed for the transient seepage analyses of semi‐infinite layers with a constant depth. The scaled boundary finite element equation of pore water pressure is formulated first in the frequency domain. With the eigenvalue problem, the equation can be decoupled into modal equations whose modal dynamic permeability equation can be determined. The continued fraction technique is adopted to formulate the continued fraction solution in the frequency domain. All constants in the solution are determined recursively at the high‐frequency limit. By introducing auxiliary variables and the continued fraction solution to the relationship between the prescribed seepage flow and the pore water pressure in the frequency domain, the open boundary condition is obtained. After transformed to the time domain, the open boundary condition is expressed as a system of fractional differential equations. No convolution integral is required. The accuracy of the analysis results increases with the increasing orders of continued fraction. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

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

15.
Numerical simulation of fluid flow in a hydrocarbon reservoir has to account for the presence of wells. The pressure of a grid cell containing a well is different from the average pressure in that cell and different from the bottom-hole pressure for the well [17]. This paper presents a study of grid pressures obtained from the simulation of single phase flow through an isotropic porous medium using different numerical methods. Well equations are proposed for Darcy flow with Galerkin finite elements and mixed finite elements. Furthermore, high velocity (non-Darcy) flow well equations are developed for cell-centered finite difference, Galerkin finite element and mixed finite element techniques.  相似文献   

16.
A hybrid discrete‐continuum numerical scheme is developed to study the behavior of a hydraulic fracture crossing natural fractures. The fully coupled hybrid scheme utilizes a discrete element model for an inner domain, within which the hydraulic fracture propagates and interacts with natural fractures. The inner domain is embedded in an outer continuum domain that is implemented to extend the length of the hydraulic fracture and to better approximate the boundary effects. The fracture is identified to propagate initially in the viscosity‐dominated regime, and the numerical scheme is calibrated by using the theoretical plane strain hydraulic fracture solution. The simulation results for orthogonal crossing indicate three fundamental crossing scenarios, which occur for various stress ratios and friction coefficients of the natural fracture: (i) no crossing, that is, the hydraulic fracture is arrested by the natural fracture and makes a T‐shape intersection; (ii) offset crossing, that is, the hydraulic fracture crosses the natural fracture with an offset; and (iii) direct crossing, that is, the hydraulic fracture directly crosses the natural fracture without diversion. Each crossing scenario is associated with a distinct net pressure history. Additionally, the effects of strength contrast and stiffness contrast of rock materials and intersection angle between the hydraulic fracture and the natural fracture are also investigated. The simulations also illustrate that the level of fracturing complexity increases as the number and extent of the natural fractures increase. As a result, we can conclude that complex hydraulic fracture propagation patterns occur because of complicated crossing behavior during the stimulation of naturally fractured reservoirs. Copyright © 2017 John Wiley & Sons, Ltd.  相似文献   

17.
A superposition scheme is proposed to obtain a fundamental solution for boundary elements in multi‐layered elastic media. A three‐layered elastic region is obtained by superposing two sets of bonded half‐planes and subtracting one infinite plane. Therefore, the solution for an element in a layered media can be expressed in terms of bonded half‐plane solutions and an infinite‐plane solution. The major advantages of this superposition scheme are: (1) it is unnecessary to introduce elements at the interface, (2) it can be extended to higher‐order element, and (3) it may be applicable to three dimensions easily. The accuracy and performance of the developed model is illustrated by two examples. For the problem of a pressurized two‐dimensional crack within a three‐layered system, the comparison with other numerical results shows the model is quite accurate and efficient. The model is also used for a study of a practical two‐dimensional mining problem in South Africa, i.e. stoping through a dyke with material properties different from the host rock. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

18.
侯晓萍  陈胜宏 《岩土力学》2020,41(4):1437-1446
采用复合单元法建立了模拟裂隙多孔介质变饱和流动的数值模型。该模型具有以下特点：裂隙不需要离散成特定单元，而是根据几何位置插入到孔隙基质单元中形成复合单元；在复合单元中，分别建立裂隙流和孔隙基质流的计算方程，二者通过裂隙?基质界面产生联系并整合成复合单元方程；复合单元方程具有和常规有限单元方程相同的格式，因此，可以使用常规有限单元方程的求解技术。采用欠松弛迭代、集中质量矩阵以及自适应时步调节等技术，开发了裂隙多孔介质变饱和流动计算程序。通过模拟一维干土入渗和复杂裂隙含水层内的流动问题，验证了该模型的合理性和适用性。模拟结果为进一步认识非饱和裂隙含水层地下水流动特性提供了理论依据。  相似文献   

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

20.
By using a complete set of poroelastodynamic spherical wave potentials (SWPs) representing a fast compressional wave PI, a slow compressional wave PII, and a shear wave S with 3 vectorial potentials (not all are independent), a solution scheme based on the method of fundamental solution (MFS) is devised to solve 3‐D wave scattering and dynamic stress concentration problems due to inhomogeneous inclusions and cavities embedded in an infinite poroelastic domain. The method is verified by comparing the result with the elastic analytical solution, which is a degenerated case, as well as with poroelastic solution obtained using other numerical methods. The accuracy and stability of the SWP‐MFS are also demonstrated. The displacement, hoop stress, and fluid pore pressure around spherical cavity and poroelastic inclusion with permeable and impermeable boundary are investigated for incident plane PI and SV waves. The scattering characteristics are examined for a range of material properties, such as porosity and shear modulus contrast, over a range of frequency. Compared with other boundary‐based numerical strategy, such as the boundary element method and the indirect boundary integral equation method, the current SWP‐MFS is a meshless method that does not need elements to approximate the geometry and is free from the treatment of singularities. The SWP‐MFS is a highly accurate and efficient solution methodology for wave scattering problems of arbitrary geometry, particularly when a part of the domain extends to infinity.  相似文献   

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

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