scholarly journals On one method of numerical modeling of a two-phase fluid system in a fractured-porous reservoir

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.

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 ◽  
Vol 2131 (2) ◽  
pp. 022001
Author(s):  
Yu O Bobreneva ◽  
P I Rahimly ◽  
Yu A Poveshchenko ◽  
V O Podryga ◽  
L V Enikeeva

Abstract A method of numerical modeling based on splitting by physical processes of two-phase fluid transfer in a formation with fractured-porous reservoirs is described. Reservoirs of this type have a natural fracture system and are described by the dual porosity model. A four-block mathematical model of the fluid redistribution between a pore-type matrix and a natural fracturing pattern is proposed and studied. The resulting system is complex and entails a number of difficulties associated with a large number of variables and the absence of important properties of a linearized system of equations, such as self-adjointness and symmetry, which are present in the description of piezoconductive processes. The complete splitting by physical processes is carried out to solve this problem. The resulting split model is differentially equivalent to the discrete initial balance equations of the system (conservation of the mass components of the fluids and the total energy of the system), written in divergent form. This approach is associated with a nonlinear approximation of the grid functions in time, which depends on the fraction of the volume occupied by the fluids in the pores, and is easy to implement.


2011 ◽  
Vol 46 (1) ◽  
pp. 90-96
Author(s):  
G. R. Idrisova ◽  
L. A. Kovaleva ◽  
M. V. Mavletov ◽  
A. A. Musin

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.


Sign in / Sign up

Export Citation Format

Share Document