A lattice Boltzmann investigation of steady-state fluid distribution, capillary pressure and relative permeability of a porous medium: Effects of fluid and geometrical properties

2018 ◽  
Vol 116 ◽  
pp. 153-166 ◽  
Author(s):  
Zi Li ◽  
Sergio Galindo-Torres ◽  
Guanxi Yan ◽  
Alexander Scheuermann ◽  
Ling Li
SPE Journal ◽  
2010 ◽  
Vol 15 (04) ◽  
pp. 917-927 ◽  
Author(s):  
Thomas Ramstad ◽  
Pål-Eric Øren ◽  
Stig Bakke

Summary We present results from simulations of two-phase flow directly on digitized rock-microstructure images of porous media using a lattice Boltzmann (LB) method. The implemented method is performed on a D3Q19 lattice with fluid/fluid and fluid/solid interaction rules to handle interfacial tension and wetting properties. We demonstrate that the model accurately reproduces capillary and wetting effects in pores with a noncircular shape. The model is applied to study viscous coupling effects for two-phase concurrent annular flow in circular tubes. Simulated relative permeabilities for this case agree with analytical predictions and show that the nonwetting-phase relative permeability might greatly exceed unity when the wetting phase is less viscous than the nonwetting phase. Two-phase LB simulations are performed on microstructure images derived from X-ray microtomography and process-based reconstructions of Bentheimer sandstone. By imposing a flow regulator to control the capillary number of the flow, the LB model can closely mimic typical experimental setups, such as centrifuge capillary pressure and unsteady- and steady-state relative permeability measurements. Computed drainage capillary pressure curves are found to be in excellent agreement with experimental data. Simulated steady-state relative permeabilities at typical capillary numbers in the vicinity of 10−5 are in fair agreement with measured data. The simulations accurately reproduce the wetting-phase relative permeability but tend to underpredict the nonwetting-phase relative permeability at high wetting-phase saturations. We explain this by pointing to percolation threshold effects of the samples. For higher capillary numbers, we correctly observe increased relative permeability for the nonwetting phase caused by mobilization and flow of trapped fluid. It is concluded that the LB model is a powerful and promising tool for deriving physically meaningful constitutive relations directly from rock-microstructure images.


SPE Journal ◽  
2017 ◽  
Vol 22 (03) ◽  
pp. 940-949 ◽  
Author(s):  
Edo S. Boek ◽  
Ioannis Zacharoudiou ◽  
Farrel Gray ◽  
Saurabh M. Shah ◽  
John P. Crawshaw ◽  
...  

Summary We describe the recent development of lattice Boltzmann (LB) and particle-tracing computer simulations to study flow and reactive transport in porous media. First, we measure both flow and solute transport directly on pore-space images obtained from micro-computed-tomography (CT) scanning. We consider rocks with increasing degree of heterogeneity: a bead pack, Bentheimer sandstone, and Portland carbonate. We predict probability distributions for molecular displacements and find excellent agreement with pulsed-field-gradient (PFG) -nuclear-magnetic-resonance (NMR) experiments. Second, we validate our LB model for multiphase flow by calculating capillary filling and capillary pressure in model porous media. Then, we extend our models to realistic 3D pore-space images and observe the calculated capillary pressure curve in Bentheimer sandstone to be in agreement with the experiment. A process-based algorithm is introduced to determine the distribution of wetting and nonwetting phases in the pore space, as a starting point for relative permeability calculations. The Bentheimer relative permeability curves for both drainage and imbibition are found to be in good agreement with experimental data. Third, we show the speedup of a graphics-processing-unit (GPU) algorithm for large-scale LB calculations, offering greatly enhanced computing performance in comparison with central-processing-unit (CPU) calculations. Finally, we propose a hybrid method to calculate reactive transport on pore-space images by use of the GPU code. We calculate the dissolution of a porous medium and observe agreement with the experiment. The LB method is a powerful tool for calculating flow and reactive transport directly on pore-space images of rock.


Author(s):  
Pål Ø. Andersen

Steady state relative permeability experiments are performed by co-injection of two fluids through core plug samples. Effective relative permeabilities can be calculated from the stabilized pressure drop using Darcy’s law and linked to the corresponding average saturation of the core. These estimated relative permeability points will be accurate only if capillary end effects and transient effects are negligible. This work presents general analytical solutions for calculation of spatial saturation and pressure gradient profiles, average saturation, pressure drop and relative permeabilities for a core at steady state when capillary end effects are significant. We derive an intuitive and general “intercept” method for correcting steady state relative permeability measurements for capillary end effects: plotting average saturation and inverse effective relative permeability (of each phase) against inverse total rate will give linear trends at high total rates and result in corrected relative permeability points when extrapolated to zero inverse total rate (infinite rate). We derive a formal proof and generalization of the method proposed by Gupta and Maloney (2016) [SPE Reserv. Eval. Eng. 19, 02, 316–330], also extending the information obtained from the analysis, especially allowing to calculate capillary pressure. It is shown how the slopes of the lines are related to the saturation functions allowing to scale all test data for all conditions to the same straight lines. Two dimensionless numbers are obtained that directly express how much the average saturation is changed and the effective relative permeabilities are reduced compared to values unaffected by end effects. The numbers thus quantitatively and intuitively express the influence of end effects. A third dimensionless number is derived providing a universal criterion for when the intercept method is valid, directly stating that the end effect profile has reached the inlet. All the dimensionless numbers contain a part depending only on saturation functions, injected flow fraction and viscosity ratio and a second part containing constant known fluid, rock and system parameters such as core length, porosity, interfacial tension, total rate, etc. The former parameters determine the saturation range and shape of the saturation profile, while the latter number determines how much the profile is compressed towards the outlet. End effects cause the saturation profile and average saturation to shift towards the saturation where capillary pressure is zero and the effective relative permeabilities to be reduced compared to the true relative permeabilities. This shift is greater at low total rate and gives a false impression of rate-dependent relative permeabilities. The method is demonstrated with multiple examples. Methodologies for deriving relative permeability and capillary pressure systematically and consistently, even based on combining data from tests with different fluid and core properties, are presented and demonstrated on two datasets from the literature. The findings of this work are relevant to accurately estimate relative permeabilities in steady state experiments, relative permeability end points and critical saturations during flooding or the impact of injection chemicals on mobilizing residual phase.


1969 ◽  
Vol 9 (02) ◽  
pp. 221-231 ◽  
Author(s):  
R. Ehrlich ◽  
F.E. Crane

Abstract A consolidated porous medium is mathematically modeled by networks of irregularly shaped interconnected pore channels. Mechanisms are described that form residual saturations during immiscible displacement both by entire pore channels being bypassed and by fluids being isolated by the movement of an interface within individual pore channels. This latter mechanism is shown to depend strongly on pore channel irregularity. Together, these mechanisms provide an explanation for the drainage-imbibition-hysteresis effect. The calculation of steady-state relative permeabilities, based on a pore-size distribution permeabilities, based on a pore-size distribution obtained from a Berea sandstone, is described. These relative permeability curves agree qualitatively with curves that are generally accepted to be typical for highly consolidated materials. In situations where interfacial effects predominate over viscous and gravitational effects, the following conclusions are reached.Relative permeability at a given saturation is everywhere independent of flow rate.Relative permeability is independent of viscosity ratio everywhere except at very low values of wetting phase relative permeability.Irreducible wetting-phase saturation following steady-state drainage decreases with increasing ratio of nonwetting- to wetting-phase viscosity.Irreducible wetting-phase saturation following unsteady-state drainage is lower than for steady-state drainage.Irreducible nonwetting-phase saturation following imbibition is independent of viscosity ratio, whether or not the imbibition is carried out under steady- or unsteady-state conditions. Experiments qualitatively verify the conclusions regarding unsteady-state residual wetting-phase saturation. Implications of these conclusions are discussed. Introduction Natural and artificial porous materials are generally composed of matrix substance brought together in a more or less random manner. This leads to the creation of a network of interconnected pore spaces of highly irregular shape. Since the pore spaces of highly irregular shape. Since the geometry of such a network is impossible to describe, we can never obtain a complete description of its flow behavior. We can, however, abstract those properties of the porous medium pertinent to the type of flow under consideration, and thus obtain an adequate description of that flow. Thus, the Kozeny-Carmen equation, by considering a porous medium as a bundle of noninterconnecting capillary tubes, provides an adequate description of single-phase provides an adequate description of single-phase flow. With the addition of a saturation-dependent tortuosity parameter in two-phase flow to account for flow path elongation, the Kozeny-Carmen equation has been used to predict relative permeabilities for the displacement of a wetting permeabilities for the displacement of a wetting liquid by a gas. It has long been recognized that relative permeability depends not only on saturation but permeability depends not only on saturation but also on saturation history as well. Naar and Henderson described a mathematical model in which differences between drainage and imbibition behavior are explained in terms of a bypassing mechanism by which oil is trapped during imbibition. Fatt proposed a model for a porous medium that consisted of regular networks of cylindrical tubes of randomly distributed radii. From this he calculated the drainage relative permeability curves. Moore and Slobod, Rose and Witherspoon, and Rose and Cleary each considered flow in a pore doublet (a parallel arrangement of a small and pore doublet (a parallel arrangement of a small and large diameter cylindrical capillary tube). They concluded that, because of the different rates of flow in each tube, trapping would occur in one of the tubes; the extent of which would depend upon viscosity ratio and capillary pressure. SPEJ p. 221


1964 ◽  
Vol 4 (03) ◽  
pp. 195-202 ◽  
Author(s):  
P.M. Blair

Abstract This paper presents numerical solutions of the equations describing the imbibition of water and the countercurrent flow of oil in porous rocks. The imbibition process is of practical importance in recovering oil from heterogeneous formations and has been studied principally by experimental means. Calculations were made for imbibition of water into both linear and radial systems. Imbibition in the linear systems was allowed to take place through one open, or permeable, face of the porous medium studied. In the radial system, water was imbibed inward from the outer radius. The effects on rate of imbibition of varying the capillary pressure and relative permeability curves, oil viscosity and the initial water saturation were computed. For each case studied, the rate of water imbibition and the saturation and pressure profiles were calculated as functions of time. The results of these calculations indicate that, for the porous medium studied, the time required to imbibe a fixed volume of water of a certain viscosity is approximately proportional to the square root of the viscosity of the reservoir oil whenever the oil viscosity is greater than the water viscosity. Results are also presented illustrating the effects on rate of imbibition of the other variables studied. Introduction The process of imbibition, or spontaneous flow of fluids in porous media under the influence of capillary pressure gradient s, occurs wherever there exist in permeable rock capillary pressure gradients which are not exactly balanced by opposing pressure gradients (such as those resulting from the influence of gravity). The importance of such capillary movement in the displacement of oil by water or gas was recognized in early investigations and described by Leverett, Lewis and True in 1942. Methods advanced by these authors for studying the process using dynamically scaled models were rendered more general and flexible by the research of later workers. The influence of capillary forces in laboratory water floods has also been discussed by several authors. While imbibition plays a very important role in the recovery of oil from normal reservoirs, Brownscombe and Dyes pointed out that imbibition might be the dominant displacement process in water flooding reservoirs characterized by drastic variations in permeability, such as in fractured- matrix reservoirs. In water-wet, fractured-matrix reservoirs, water will be imbibed from fractures into the matrix with a countercurrent expulsion of oil into the fractures. If the imbibition occurs at a sufficiently rapid rate, a very successful water flood can result; if the imbibition proceeds slowly the project might not be economically attractive. Scaled-model studies have demonstrated the vital importance of imbibition in secondary recovery in fractured reservoirs. It is therefore important in the evaluation of waterflooding prospects to develop a thorough understanding of the quantitative relationships of the factors which control the rapidity of capillary imbibition. The imbibition process serves reservoir engineers in still another important way by providing a technique for studying the wettability of reservoir core samples. Such experiments are usually conducted by observing the rate of expulsion of oil or water from core samples submerged in the appropriate fluid. Several papers have been published on the experimental techniques involved. Although Handy has recently published a method for calculating capillary pressures from experiments with gas-saturated cores, it has not yet been possible to deduce quantitative information regarding water-oil relative permeability and capillary pressure characteristics of the rock from the experimental results. Thus a technique is needed for studying the quantitative dependence of imbibition rate on oil and water viscosity, initial water saturation, relative permeability-saturation, and capillary pressure-saturation relations. The development of such information, including saturation and pressure profiles by laboratory experiments, would be very difficult. SPEJ P. 195ˆ


Sign in / Sign up

Export Citation Format

Share Document