共查询到20条相似文献，搜索用时 31 毫秒
1.
目前国内关于CO2咸水层封存尚处于先导性和试验性研究阶段,对超临界CO2注入过程中岩层力学响应和流体运移的理论与技术方面的认识还不完善。为研究CO2注入下岩层变形和流体运移,基于两相流动数学模型,给出了超临界CO2和咸水质量守恒方程;采用毛细压力和有效饱和度的关系式,将质量守恒方程变换成以毛细压力为变量的表达式,以便于考虑流体压力对岩层的影响。提出了无流体压力影响下的岩层力学本构模型,该模型能够同时考虑岩层的塑性变形和损伤。分析了两相流体-岩层相互作用机制:一方面,采用有效应力原理,考虑流体压力对岩层的力学影响;另一方面,通过岩层固有渗透率变化考虑岩层变形对流体运移的影响。 相似文献
2.
An iterative method is presented for solving a fully coupled and implicit formulation of fluid flow in a porous medium. The mathematical model describes a set of fully coupled three-phase flow of compressible and immiscible fluids in a saturated oil reservoir. The finite element method is applied to obtain the simultaneous solution (SS) for the resulting highly non-linear partial differential equations where fluid pressures are the primary unknowns. The final discretized equations are solved iteratively by using a fully implicit numerical scheme. Several examples, illustrating the use of the present model, are described. The increased stability achieved with this scheme has permitted the use of larger time steps with smaller material balance errors. 相似文献
3.
A numerical simulation is presented for three-dimensional three-phase fluid flow in a deforming saturated oil reservoir. The mathematical formulation describes a fully coupled governing equation systen which consists of the equilibrium and continuity equations for three immiscible fluids flowing in a porous medium. An elastoplastic soil model, based on a Mohr–Coulomb yield surface, is used. The finite element method is applied to obtain simultaneous solutions to the governing equations where displacement and fluid pressures are the primary unknowns. The final discretized equations are solved by a direct solver using fully implicit procedures. The developed model is used to investigate the problems of three-phase fluid flow in a deforming saturated oil reservoir. 相似文献
4.
H. F. Wang 《Contributions to Mineralogy and Petrology》1993,114(3):357-364
The concept of a double porosity medium to model fluid flow in fractured rock has been applied to model diffusion in rock containing a small amount of a continuous fluid phase that surrounds small volume elements of the solid matrix. The model quantifies the relative role of diffusion in the fluid and solid phases of the rock. The fluid is the fast diffusion path, but the solid contains the volumetrically significant amount of the diffusing species. The double medium model consists of two coupled differential equations. One equation is the diffusion equation for the fluid concentration; it contains a source term for change in the average concentration of the diffusing species in the solid matrix. The second equation represents the assumption that the change in average concentration in a solid element is proportional to the difference between the average concentration in the solid and the concentration in the fluid times the solid-fluid partition coefficient. The double medium model is shown to apply to laboratory data on iron diffusion in fluid-bearing dunite and to measured oxygen isotope ratios at marble-metagranite contacts. In both examples, concentration profiles are calculated for diffusion taking place at constant temperature, where a boundary value changes suddenly and is subsequently held constant. Knowledge of solid diffusivities can set a lower bound to the length of time over which diffusion occurs, but only the product of effective fluid diffusivity and time is constrained for times longer than the characteristic solid diffusion time. The double medium results approach a local, grain-scale equilibrium model for times that are large relative to the time constant for solid diffusion. 相似文献
5.
6.
用流固耦合方法研究油藏压裂后应力应变和孔渗特性变化 总被引：10，自引：1，他引：10
油藏压裂后将引起地应力场发生变化，使岩石变形，导致孔隙度和渗透率变化，进而影响产量，为研究这一问题，作者建立了油藏压裂后流－固耦合渗流模型，考虑了以下因素：油藏岩石变形，地应力，孔隙度和渗透率变化，人工裂缝，流体渗流与岩石应变耦合，储藏渗流与裂缝渗流耦合，非达西效应等。较详细地给出了耦合方程及推导过程，控制方程包括的未知变量有压力，饱和度及位移，11个变量，和11个方程，用有限差分方法将流体渗流和岩石应变方程离散成主对角占优的七对角矩阵，可在修改已有三维二相渗流和三维固体力学程序的基础上，采用隐式迭代方法求解，示例分析表明，用此模型可以研究储层应力变变，孔隙度和渗透率随时间和空间变化规律，为开发方案制定，整体压裂设计，压后生产管理等方面提供定量分析技术。 相似文献
7.
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. 相似文献
8.
The coupled heat-fluid-stress problem of circular wellbore or spherical cavity subjected to a constant temperature change and a constant fluid flow rate is considered. Transient analytical solutions for temperature, pore pressure and stress are developed by coupling conductive heat transfer with Darcy fluid flow in a poroelastic medium. They are applicable to low permeability porous media suitable for liquid-waste disposal and also simulating reservoir for enhanced oil recovery, where conduction dominates the heat transfer process. A full range of solutions is presented showing separately the effects of temperature and fluid flow on pore pressure and stress development. It is shown that injection of warm fluid can be used to restrict fracture development around wellbores and cavities and generally to optimize a fluid injection operation. Both the limitations of the solutions and the convective flow effect are addressed. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
9.
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. 相似文献
10.
This paper presents a model which can be used for fast landslides where coupling between solid and pore fluid plays a fundamental role. The proposed model is able to describe debris flows where the difference of velocities between solid grains and fluid is important. The approach is based on the mathematical model proposed by Zienkiewicz and Shiomi, which is similar to those of Pitman and Le and Pudasaini. The novelty of the present work is the numerical technique used, the smoothed particle hydrodynamics (SPH). We propose to use a double set of nodes for soil and water phases, the interaction between them being described by a suitable drag law. The paper presents both mathematical and numerical models, describing the main assumptions and their limitations. Then, the model is applied to (1) a simple case where shocks and expansion waves appear, (2) a dam break problem on a horizontal plane with a frictional soil phase, and (3) a debris flow which happened in Hong Kong. The main conclusions that can be drawn from the applications are:
- Debris flows having 2 phases with important relative mobility present a rich structure of shocks and rarefaction waves, which has to be properly modeled. Otherwise, the model will have numerical damping or dispersion.
- Dambreak exercises provide interesting information in simple and controlled situations. We can see how both phases move relative to each other.
- Real debris flows can be simulated with the proposed model, obtaining reasonable results.
11.
The surface subsidence above a compacting saturated oil reservoir is the main topic of this paper. From a literature review, it is obvious that extensive efforts have been conducted for investigating this phenomenon in various situations. Herein, a numerical model, based on the finite element method, was used for simulating three-dimensional three-phase fluid flow in a deforming saturated oil reservoir. The mathematical formulation describes a fully coupled governing equation system which consists of the equilibrium and continuity equations for three immiscible fluids flowing in a porous media. An elastoplastic soil model, based on a Mohr Coulomb yield surface, was utilized. The finite element method was applied to obtain simultaneous solutions to the governing equations where the displacements and the fluid pressures are the primary unknowns. The final discretized equations are solved by a direct solver using fully implicit procedures. A linear analysis was used to study the stability conditions of the present model. Finally, a series of simulations were conducted to indicate the validity and the utility of the developed model. 相似文献
12.
Dynamic data integration for structural modeling: model screening approach using a distance-based model parameterization 总被引：2，自引：0，他引：2
This paper proposes a novel history-matching method where reservoir structure is inverted from dynamic fluid flow response.
The proposed workflow consists of searching for models that match production history from a large set of prior structural
model realizations. This prior set represents the reservoir structural uncertainty because of interpretation uncertainty on
seismic sections. To make such a search effective, we introduce a parameter space defined with a “similarity distance” for
accommodating this large set of realizations. The inverse solutions are found using a stochastic search method. Realistic
reservoir examples are presented to prove the applicability of the proposed method. 相似文献
13.
M. D. Sharma 《Journal of Earth System Science》2008,117(6):951-958
Biot’s theory for wave propagation in saturated porous solid is modified to study the propagation of thermoelastic waves in
poroelastic medium. Propagation of plane harmonic waves is considered in isotropic poroelastic medium. Relations are derived
among the wave-induced temperature in the medium and the displacements of fluid and solid particles. Christoffel equations
obtained are modified with the thermal as well as thermoelastic coupling parameters. These equations explain the existence
and propagation of four waves in the medium. Three of the waves are attenuating longitudinal waves and one is a non-attenuating
transverse wave. Thermal properties of the medium have no effect on the transverse wave. The velocities and attenuation of
the longitudinal waves are computed for a numerical model of liquid-saturated sandstone. Their variations with thermal as
well as poroelastic parameters are exhibited through numerical examples. 相似文献
14.
Hydraulic fracturing in unconventional gas reservoirs: risks in the geological system, part 2 总被引：2，自引：2，他引：0
Alexander Kissinger Rainer Helmig Anozie Ebigbo Holger Class Torsten Lange Martin Sauter Michael Heitfeld Johannes Klünker Wiebke Jahnke 《Environmental Earth Sciences》2013,70(8):3855-3873
Hydraulic fracturing is a method used for the production of unconventional gas resources. Huge amounts of so-called fracturing fluid (10,000–20,000 m3) are injected into a gas reservoir to create fractures in solid rock formations, upon which mobilised methane fills the pore space and the fracturing fluid is withdrawn. Hydraulic fracturing may pose a threat to groundwater resources if fracturing fluid or brine can migrate through fault zones into shallow aquifers. Diffuse methane emissions from the gas reservoir may not only contaminate shallow groundwater aquifers, but also escape into the atmosphere where methane acts as a greenhouse gas. The working group “Risks in the Geological System” as part of ExxonMobil’s hydrofracking dialogue and information dissemination processes was tasked with the assessment of possible hazards posed by migrating fluids as a result of hydraulic fracturing activities. In this work, several flow paths for fracturing fluid, brine and methane are identified and scenarios are set up to qualitatively estimate under what circumstances these fluids would leak into shallower layers. The parametrisation for potential hydraulic fracturing sites in North Rhine-Westphalia and Lower Saxony (both in Germany) is derived from literature using upper and lower bounds of hydraulic parameters. The results show that a significant fluid migration is only possible if a combination of several conservative assumptions is met by a scenario. 相似文献
15.
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. 相似文献
16.
We present a fully implicit formulation of coupled flow and geomechanics for fractured three-dimensional subsurface formations. The Reservoir Characterization Model (RCM) consists of a computational grid, in which the fractures are represented explicitly. The Discrete Fracture Model (DFM) has been widely used to model the flow and transport in natural geological porous formations. Here, we extend the DFM approach to model deformation. The flow equations are discretized using a finite-volume method, and the poroelasticity equations are discretized using a Galerkin finite-element approximation. The two discretizations—flow and mechanics—share the same three-dimensional unstructured grid. The mechanical behavior of the fractures is modeled as a contact problem between two computational planes. The set of fully coupled nonlinear equations is solved implicitly. The implementation is validated for two problems with analytical solutions. The methodology is then applied to a shale-gas production scenario where a synthetic reservoir with 100 natural fractures is produced using a hydraulically fractured horizontal well. 相似文献
17.
流固耦合问题及研究进展 总被引：16，自引：0，他引：16
传统的渗流理论一般假设流体流动的多孔介质骨架是完全刚性的,即在孔隙流体压力变化过程中,固体骨架不产生任何弹性或塑性变形,这时可将渗流作为非耦合问题来研究。这种简化虽然可以得到问题的近似解,但存在许多缺陷,而且也不切合生产实际。比如:在油田开采过程中,孔隙流体压力会逐渐降低,将导致储层内有效应力的变化,使储层产生变形。近年来,流固耦合问题越来越受到人们的重视,这方面的研究涉及许多领域。该文介绍了有关工程涉及到的流固耦合问题,重点针对油、气开采问题,介绍了储层流固耦合渗流的特点及研究方法和理论进展,包括单相、多相流体渗流的流固耦合数学模型及有限元数值模型。 相似文献
18.
A micro‐hydromechanical model for granular materials is presented. It combines the discrete element method for the modeling of the solid phase and a pore‐scale finite volume formulation for the flow of an incompressible pore fluid. The coupling equations are derived and contrasted against the equations of conventional poroelasticity. An analogy is found between the discrete element method pore‐scale finite volume coupling and Biot's theory in the limit case of incompressible phases. The simulation of an oedometer test validates the coupling scheme and demonstrates the ability of the model to capture strong poromechanical effects. A detailed analysis of microscale strain and stress confirms the analogy with poroelasticity. An immersed deposition problem is finally simulated and shows the potential of the method to handle phase transitions. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
19.
This paper presents a numerical model for simulating free surface flow in porous media with spatially varying porosity. The governing equations are based on the mixture theory. The resistance forces between solid and fluid is assumed to be nonlinear. A multiphase SPH approach is presented to solve the governing equations. In the multiphase SPH, water is modeled as a weakly compressible fluid, and solid phase is discretized by fixed solid particles carrying information of porosity. The model is validated by several numerical examples including seepage through specimen, fast flow through rockfill dam and wave interaction with porous structure. Good agreements between numerical results and experimental data are obtained in terms of flow rate and evolution of free surface. Parameter study shows that (1) the nonlinear resistance law provides more accurate results; (2) particle size and porosity have significant influence on the porous flow. 相似文献
20.
构造-流体-成矿体系的反应-输运-力学耦合模型和动力学模拟 总被引：12，自引：0，他引：12
建立了一个综合的构造流体成矿体系的反应输运力学耦合动力学模型。利用有限元方法求解岩石变形、断裂作用和断裂网络统计动力学、流体流动、有机和无机地球化学反应及成岩成矿作用、压力溶液和其它压实力学、热迁移的方程组 ,可以对构造流体成矿体系的动力学演化过程进行 1～ 3维数值模拟。模拟的主要内容是在各种过程耦合作用下描述构造流体成矿体系的主要变量的时空演化 :( 1)与成矿流体的形成和性质有关的变量 ,如地层中矿物 (包括成矿物质 )的溶解速率、流体中各组分的浓度与饱和度、流体温度、压力、离子强度等 ;( 2 )与构造变形和流体运移有关的各变量 ,如应力与变形速率、岩石孔隙度、构造 (断裂 )渗透率等 ;( 3 )与沉淀成矿有关的变量 ,如矿物 (金属矿物和脉石矿物 )的成核速率、各矿物的沉淀量等 ;( 4 )上述各有关变量间的时空耦合关系 ,如断裂渗透率时空演化与流体流动、汇聚和成矿的耦合关系等。以湖南沃溪金锑钨矿床为例 ,应用该模型和方法对成矿动力学过程和动力学机制进行了初步的模拟与分析。 相似文献