A New Two Phase Extension of Modified Brinkman Formulation for Fluid Flow through Porous Media

2005 ◽  
Author(s):  
Hadi Belhaj ◽  
Shabbir Mustafiz ◽  
Fuxi Ma ◽  
M. R. Islam

In porous media research, Modified Brinkman’s equation is a very recent development. It is important as it incorporates the concept of viscous effect to inertial effect in a fluid flow system when Darcy’s, Forchheimer’s and Brinkman’s terms are brought all together. So far, researchers have developed the modified equation in its two-dimensional forms; however, limited to only one phase. In reality, petroleum reservoirs experience the multiphase conditions. Therefore, the simulation of a multidimensional, multiphase scenario is mostly desired, the highlight of this paper. The paper presents the formulation of two-dimensional, transient pressure and saturation equations for oil and water phases, one equation for each phase. The difference between phases is noticeable explicitly in their respective saturation, permeability, viscosity and velocity terms. The equations are then solved numerically to generate relative permeability curves. The simultaneous solution of pressure and saturation terms in the governing equations required additional relationships: the phase saturation constraint and capillary pressure as function of saturation. Finally, the numerical results are compared and validated with the experimental results. The implication of this study is manifold. The formulated equations including the solution part for the multiphase conditions are new. The new comprehensive model will describe fluid flow in reservoirs prone to high velocity or fractures more accurately than ever described by Darcy’s or other aforementioned equations.

1971 ◽  
Vol 11 (04) ◽  
pp. 419-425 ◽  
Author(s):  
Carlon S. Land

Abstract Two-phase imbibition relative permeability was measured in an attempt to validate a method of calculating imbibition relative permeability. The stationary-liquid-phase method was used to measure several hysteresis loops for alundum and Berea sandstone samples. The method of calculating imbibition relative permeability is described, and calculated relative permeability curves are compared with measured curves. The calculated relative Permeability is shown to be a reasonably good Permeability is shown to be a reasonably good approximation of measured values if an adjustment is made to some necessary data. Due to the compressibility of gas, which is used as the nonwetting phase, a correction to the measured trapped gas saturation is necessary to make it agree with the critical gas saturation of the imbibition relative permeability curve. Introduction The existence of hysteresis in the relationship of relative permeability to saturation has been recognized for many yews. Geden et al. and Osoba et al. called attention to the occurrence of hysteresis and the importance of the direction of saturation change on the relative permeability-saturation relations. It is generally believed that relative permeability is a function of saturation alone for a permeability is a function of saturation alone for a given direction of saturation change, but that there is a distinct difference in relative permeability curves for saturation changes in different directions. The reservoir engineer should be aware of this hysteresis, and he should select the relative permeability curve which is appropriate for the permeability curve which is appropriate for the recovery process of interest. The directions of saturation change have been designated "drainage" and "imbibition" in reference to changes in the wetting-phase saturation. In a two-phase system, an increase in the wetting-phase saturation is referred to as imbibition, while a decrease in wetting-phase saturation is called drainage. The solution-gas-drive recovery mechanism is controlled by relative permeability to oil and gas in which the saturation of oil, the wetting phase, is decreasing. In waterflooding a water-wet reservoir rock, the saturation of water, the wetting phase, is increasing. These two sets of relative permeability curves, gas-oil and oil-water, do not have the same relationship to the wetting-phase saturation. This difference is not due to the difference in fluid properties, but is a result of the difference in properties, but is a result of the difference in direction of saturation change. The flow properties of the drainage and imbibition systems differ because of the entrapment of the nonwetting phase during imbibition. As drainage occurs, the nonwetting phase occupies the most favorable flow channels. During imbibition, part of the nonwetting phase is bypassed by the increasing wetting phase, leaving a portion of the nonwetting phase in an immobile condition. This trapped part phase in an immobile condition. This trapped part of the nonwetting phase saturation does not contribute to the flow of that phase, and at a given saturation the relative permeability to the nonwetting phase is always less in the imbibition direction phase is always less in the imbibition direction than in the drainage direction. The concept that some of the nonwetting phase is mobile and some is immobile during a saturation change in the imbibition direction previously was used to develop equations for imbibition relative permeability. In this development, it was assumed permeability. In this development, it was assumed that the amount of entrapment at any saturation can be obtained from the relationship between initial nonwetting-phase saturations established in the drainage direction and residual saturations after complete imbibition. The equations for imbibition relative permeability were not verified by laboratory measurements. The purpose of this report is m give the results of a laboratory study of imbibition relative permeability and to present a comparison of calculated relative permeability with relative permeability from laboratory measurements. permeability from laboratory measurements. In two-phase systems, hysteresis is more prominent in the relative permeability to the nonwetting phase than in that to the wetting phase. The hysteresis in the wetting-phase relative permeability is believed to be very small, and thus difficult to distinguish tom normal experimental error. SPEJ P. 419


Poromechanics ◽  
2020 ◽  
pp. 333-338
Author(s):  
M. Bai ◽  
F. Meng ◽  
J.-C. Roegiers ◽  
Y. Abousleiman

2019 ◽  
Vol 142 (4) ◽  
Author(s):  
Hamed Movahedi ◽  
Mehrdad Vasheghani Farahani ◽  
Mohsen Masihi

Abstract In this paper, we present a computational fluid dynamics (CFD) model to perform single- and two-phase fluid flow simulation on two- and three-dimensional perforated porous media with different perforation geometries. The finite volume method (FVM) has been employed to solve the equations governing the fluid flow through the porous media and obtain the pressure and velocity profiles. The volume of fluid (VOF) method has also been utilized for accurate determination of the volume occupied by each phase. The validity of the model has been achieved via comparing the simulation results with the available experimental data in the literature. The model was used to analyze the effect of perforation geometrical parameters (length and diameter), degree of heterogeneity, and also crushed zone properties (permeability and thickness) on the pressure and velocity profiles. The two-phase fluid flow around the perforation tunnel under the transient flow regime was also investigated by considering a constant mass flow boundary condition at the inlet. The developed model successfully predicted the pressure drop and resultant temperature changes for the system of air–water along clean and gravel-filled perforations under the steady-state conditions. The presented model in this study can be used as an efficient tool to design the most appropriate perforation strategy with respect to the well characteristics and reservoir properties.


2013 ◽  
Vol 88 (5) ◽  
Author(s):  
Marion Erpelding ◽  
Santanu Sinha ◽  
Ken Tore Tallakstad ◽  
Alex Hansen ◽  
Eirik Grude Flekkøy ◽  
...  

1966 ◽  
Vol 6 (04) ◽  
pp. 350-362 ◽  
Author(s):  
K.H. Coats ◽  
M.H. Terhune

Abstract Analysis and example applications have been performed to compare the accuracy and computing speed of alternating-direction explicit and implicit procedures (ADEP and ADIP) in numerical solution of reservoir fluid flow problems. ADIP yields significantly greater accuracy and requires about 60 per cent more computing time than ADEP, not 300 or 500 per cent more as reported elsewhere. Introduction Several recent papers discuss an alternating-direction explicit difference approximation (ADEP) to the diffusion equation. Example applications of ADEP and ADIP were reported to support conclusions that ADEP is comparable in accuracy to ADIP and requires one-fifth to one-third the computing time of ADIP. Applications of ADEP in calculation of two-phase flow in reservoirs was also proposed. This study was performed to compare further the relative merits of ADEP and ADIP in simulation of two-dimensional flow of one and two fluid phases in reservoirs. Since two-phase flow equations are often essentially elliptic rather than parabolic, the efficiency of ADEP in solving the elliptic equation was also examined. ADIP AND ADEP DIFFERENCE EQUATIONS The diffusion equation: ...................(1) governs heat conduction, molecular diffusion and slightly compressible fluid flow through porous media for the case of homogeneous, isotropic media. The ADEP procedure involves replacement of Eq. 1 at odd time steps by: ,.................(2) and at even time steps by: ,.................(3) where Sweeping a two-dimensional grid from southwest to northeast using Eq. 2 and from northeast to southwest using Eq. 3 allows direct (explicit) calculation of u at the new time step at each grid point. SPEJ P. 350ˆ


Sign in / Sign up

Export Citation Format

Share Document