scholarly journals Numerical modeling of multiphase mass transfer processes in fractured-porous reservoirs

2021 ◽  
Vol 2131 (2) ◽  
pp. 022002
Author(s):  
Yu O Bobreneva ◽  
P I Rahimly ◽  
Yu A Poveshchenko ◽  
V O Podryga ◽  
L V Enikeeva

Abstract The paper presents an algorithm for solving the problem of the process of mass transfer of a two-phase fluid in a fractured-porous reservoir in a one-dimensional formulation. The presence of natural fractures in such reservoirs impedes various types of exploration and field development. Fractured-porous reservoirs are characterized by intense exchange fluid flow between fractures and porous blocks. Each system has its own individual set of filtration-capacity parameters, and this fact complicates the problem under consideration. To study the mass transfer of a two-phase fluid in a medium with double porosity, a four-block mathematical model with splitting by physical processes is proposed. The model is described by a system of partial differential equations. The splitting method forms two functional blocks on the water saturation and the piezoconductivity. For the numerical solution of this system, an absolutely stable implicit finite-difference scheme is constructed in the one-dimensional case. On the basis of the proposed difference scheme, pressures and saturations in the fracture system and matrix are obtained.

2021 ◽  
pp. 1-20
Author(s):  
Yulia Olegovna Bobreneva ◽  
Parvin Ilgar gizi Rahimly ◽  
Victoria Olegovna Podryga ◽  
Svetlana Sergeevna Bazhitova ◽  
Ahmed Elsaid Ezeldin Bakeer Ali Bakeer ◽  
...  

In this work, the authors propose an algorithm for solving the problem of the process of mass transfer of a two-phase fluid in a fractured-porous reservoir in a one-dimensional formulation. The presence of natural fractures in such reservoirs impedes various types of exploration during field development. Fractured porous reservoirs are characterized by intense exchange fluid flow between fractures and porous blocks. Each system under consideration has its own individual set of filtration-capacity parameters that complicates the problem. To study the mass transfer of a two-phase liquid in a medium with double porosity, a four-block mathematical model with splitting by physical processes is proposed. The model is described by a system of partial differential equations. The method of splitting by physical processes forms two functional blocks: by water saturation and piezoconductivity. For the numerical solution of this system, an absolutely stable implicit finite-difference scheme is made in the spatially one-dimensional case. On the basis of the proposed difference scheme, pressures and saturations in the matrix and fracture system are calculated.


1972 ◽  
Vol 12 (02) ◽  
pp. 171-184 ◽  
Author(s):  
N. Van-Quy ◽  
P. Simandoux ◽  
J. Corteville

Abstract This paper describes a general multicomponent two-phase flow model, taking into account convection, diffusion and thermodynamic exchange between phases. The main assumptions are: isothermal one-dimensional flow; two-phase flow (gas and liquid); each phase may be represented by a mixture of three components or groups of components. Actually, a great many recovery problems cannot be pictured by usual models because the oil and, in many cases, the injected fluid are not simple fluids and may bring about exchanges of components that considerably modify their characteristics. This is why efforts are now being made to develop "compositional" or "multicomponent" models capable of solving such situations. Generalization of the model to more complex systems can be considered. Cases treated may be any type of single- and two-phase flow, in particular any miscible process (e. g., high-pressure gas drive, condensing gas drive, slug displacement) and any diphasic processes with high mass exchange (e.g., displacement by carbon dioxide or flue gas). This model is working and has been successfully checked by experiments. Introduction Many investigations, broth experimental and theoretical, have been made on the recovery of oil from reservoirs. With regard to mathematical models, most of those conceived up to now have dealt with oil recovery by the injection of a fluid that is miscible or immiscible with the oil. For miscible drives, single-phase flow with a binary mixture and miscibility in all proportions is involved. In such an ideal situation oil recovery is theoretically total. For immiscible displacements flow is diphasic. Capillary pressure and relative permeability play a preponderant role. Since irreducible oil saturation preponderant role. Since irreducible oil saturation is inevitable, oil recovery can never be total. Actually, a great many recovery problems cannot be pictured by such models because the oil and, in many cases, the injected fluid are not simple fluids and may bring about exchanges of components that considerably modify their characteristics. This is why efforts are now being made to develop "compositional" or "multicomponent" models capable of solving such situations. Such a model is described here. It is designed to handle the most general case of the displacement of one fluid by another. This model offers the following possibilities.The fluids may be made up of more than two components.Flow may be entirely monophasic, entirely diphasic, or partially monophasic and diphasic.Miscibility may be partial or total.The material exchange between phases may take place under specific thermodynamic conditions. A model that is much closer to reality should provide more thorough knowledge of mass transfer provide more thorough knowledge of mass transfer mechanisms in a complex mixture as well as better oil recovery forecasting with the injection of a second fluid. DESCRIPTION OF THE MODEL In a porous formation, we will consider the displacement of a liquid hydrocarbon complex in place by another fluid that is injected into the place by another fluid that is injected into the formation. The injected fluid may be a gas or a liquid, containing or not containing hydrocarbons. We assume that the mass transfer in the transition zone between the displacing fluid and the displaced fluid occurs according to three mechanism: convection, diffusion and thermodynamic exchange between phases. We propose to study the flow thus described. The main assumptions are:flow is isothermal and one-dimensional;the porous medium is homogeneous and isotropic;there is no effect of gravity;there is a two-phase flow, i.e., oil and gaseach phase may be represented by a mixture of three components or three groups of components (e.g., C1, C2-6, C7+); SPEJ P. 171


2019 ◽  
pp. 39-46
Author(s):  
Oleg B. Bocharov ◽  
Igor G. Telegin

It has been repeatedly noted in the literature that the physic-chemical properties of oil extracted at different stages of exploitation can vary significantly. We have investigated a modification of the Buckley — Leverett model (BL model) in our previous works to study the influence of oil viscosity change on the solutions of twophase filtration problems. In this model the change in oil viscosity was modeled using the dependence of the viscosity ratio on dynamic water saturation. The solutions to the linear displacement problem at a given total filtration rate were analyzed. This article considers a modification of the two-phase isothermal Muskat — Leverett model (ML model), which takes into account capillary forces. The influence of oil viscosity change in the process of oil production on the structure of the solution is studied numerically a case study of the problem of radial displacement with a given pressure drop.


Energies ◽  
2019 ◽  
Vol 12 (10) ◽  
pp. 1854
Author(s):  
Guoshu Huang ◽  
Huolin Ma ◽  
Xiangyun Hu ◽  
Jianchao Cai ◽  
Jiabin Li ◽  
...  

Temperature is one of the most prominent factors affecting production operations, predicting the accurate wellbore-formation temperature in a water-flooding production well is of great importance for multiple applications. In this paper, an improved coupled model of oil–water two-phase fluid flow and heat transfer was developed to investigate the transient temperature behavior for a producing well with multiple pay zones. Firstly, a novel method was derived to simulate the water saturation and the water breakthrough time (WBT) for tubing, which are key monitoring parameters in the process of water flooding. Then, we incorporated water saturation and an equation set for immiscible displacement to calculate the seepage velocity and the pressure of the two-phase fluid in the pay zones. Next, the upward seepage velocity of the tubing fluid change with depth was focused on, and the proper coupled initial and boundary conditions are presented at the interfaces, therewith the implicit finite difference method was used to compute the transient temperature with the input of the seepage characteristics for the reservoirs. Meanwhile, the validity of the proposed model has been verified by the typical model. Finally, a sensitivity analysis delineated that the production rate and the production time had a significant impact on the tubing fluid temperature. The overburden was hotter with a lower volumetric heat capacity or a higher thermal conductivity. In addition, the sensitivity of the porosity and the irreducible water saturation to formation temperature was significantly different before and after the WBT. The coupled model presented herein helps to advance the transient seepage characteristics analysis of pay zones, the precise temperature prediction is very useful for reservoir characterization and production analysis purposes and provides insight for designing the exploitation scheme in deep reservoirs and geothermal resources.


2005 ◽  
Vol 128 (1) ◽  
pp. 79-87 ◽  
Author(s):  
Yaqing Gu ◽  
Dean S. Oliver

This paper reports the use of ensemble Kalman filter (EnKF) for automatic history matching. EnKF is a Monte Carlo method, in which an ensemble of reservoir state variables are generated and kept up-to-date as data are assimilated sequentially. The uncertainty of reservoir state variables is estimated from the ensemble at any time step. Two synthetic problems are selected to investigate two primary concerns with the application of the EnKF. The first concern is whether it is possible to use a Kalman filter to make corrections to state variables in a problem for which the covariance matrix almost certainly provides a poor representation of the distribution of variables. It is tested with a one-dimensional, two-phase waterflood problem. The water saturation takes large values behind the flood front, and small values ahead of the front. The saturation distribution is bimodal and is not well modeled by the mean and variance. The second concern is the representation of the covariance via a relatively small ensemble of state vectors may be inadequate. It is tested by a two-dimensional, two-phase problem. The number of ensemble members is kept the same as for the one-dimensional problem. Hence the number of ensemble members used to create the covariance matrix is far less than the number of state variables. We conclude that EnKF can provide satisfactory history matching results while requiring less computation work than traditional history matching methods.


Sign in / Sign up

Export Citation Format

Share Document