共查询到20条相似文献，搜索用时 531 毫秒
1.
Pore-scale forces have a significant effect on the macroscopic behaviour of multiphase flow through porous media. This paper studies the effect of these forces using a new volume-of-fluid based finite volume method developed for simulating two-phase flow directly on micro-CT images of porous media. An analytical analysis of the relationship between the pore-scale forces and the Darcy-scale pressure drops is presented. We use this analysis to propose unambiguous definitions of Darcy-scale viscous pressure drops as the rate of energy dissipation per unit flow rate of each phase, and then use them to obtain the relative permeability curves. We show that this definition is consistent with conventional laboratory/field measurements by comparing our predictions with experimental relative permeability. We present single and two-phase flow simulations for primary oil injection followed by water injection on a sandpack and a Berea sandstone. The two-phase flow simulations are presented at different capillary numbers which cover the transition from capillary fingering at low capillary numbers to a more viscous fingering displacement pattern at higher capillary numbers, and the effect of capillary number on the relative permeability curves is investigated. Overall, this paper presents a new finite volume-based methodology for the detailed analysis of two-phase flow directly on micro-CT images of porous media and upscaling of the results to the Darcy scale. 相似文献
2.
When nonwetting fluid displaces wetting fluid in a porous rock many rapid pore-scale displacement events occur. These events are often referred to as Haines jumps and any drainage process in porous media is an ensemble of such events. However, the relevance of Haines jumps for larger scale models is often questioned. A common counter argument is that the high fluid velocities caused by a Haines jump would average-out when a bulk representative volume is considered. In this work, we examine this counter argument in detail and investigate the transient dynamics that occur during a Haines jump. In order to obtain fluid–fluid displacement data in a porous geometry, we use a micromodel system equipped with a high speed camera and couple the results to a pore-scale modeling tool called the Direct HydroDynamic (DHD) simulator. We measure the duration of a Haines jump and the distance over which fluid velocities are influenced because this sets characteristic time and length scales for fluid–fluid displacement. The simulation results are validated against experimental data and then used to explore the influence of interfacial tension and nonwetting phase viscosity on the speed of a Haines jump. We find that the speed decreases with increasing nonwetting phase viscosity or decreasing interfacial tension; however, for the same capillary number the reduction in speed can differ by an order of magnitude or more depending on whether viscosity is increased or interfacial tension is reduced. Therefore, the results suggest that capillary number alone cannot explain pore-scale displacement. One reason for this is that the interfacial and viscous forces associated with fluid–fluid displacement act over different length scales, which are not accounted for in the pore-scale definition of capillary number. We also find by analyzing different pore morphologies that the characteristic time scale of a Haines jump is dependent on the spatial configuration of fluid prior to an event. Simulation results are then used to measure the velocity field surrounding a Haines jump and thus, measure the zone of influence, which extends over a distance greater than a single pore. Overall, we find that the time and length scales of a Haines jump are inversely proportional, which is important to consider when calculating the spatial and temporal averages of pore-scale parameters during fluid–fluid displacement. 相似文献
3.
Interplay between capillary, gravity and viscous forces in unsaturated porous media gives rise to a range of complex flow phenomena affecting morphology, stability and dynamics of wetting and drainage fronts. Similar average phase contents may result in significantly different fluid distribution and patterns affecting macroscopic transport properties of the unsaturated medium. The formulation of general force balance within simplified pore spaces yields scaling relationships for motion of liquid elements in which gravitational force in excess of capillary pinning force scales linearly with viscous force. Displacement fluid front morphology is described using dimensionless force ratios expressed as Bond and Capillary numbers. The concise representations of a wide range of flow regimes with scaling relations, and predictive capabilities of front morphology based on dimensionless numbers lend support to certain generalizations. Considering available experimental data, we are able to define conditions for onset of unstable and intermittent flows leading to enhanced liquid and gas entrapment. These results provide a basis for delineation of a tentative value of Bo ∼ 0.05 as an upper limit of applicability of the Richards equation (at pore to sample scales) and related continuum-based flow models. 相似文献
4.
In porous media, the dynamics of the invading front between two immiscible fluids is often characterized by abrupt reconfigurations caused by local instabilities of the interface. As a prototype of these phenomena we consider the dynamics of a meniscus in a corner as it can be encountered in angular pores. We investigate this process in detail by means of direct numerical simulations that solve the Navier–Stokes equations in the pore space and employ the Volume of Fluid method (VOF) to track the evolution of the interface. We show that for a quasi-static displacement, the numerically calculated surface energy agrees well with the analytical solutions that we have derived for pores with circular and square cross sections. However, the spontaneous reconfigurations are irreversible and cannot be controlled by the injection rate: they are characterized by the amount of surface energy that is spontaneously released and transformed into kinetic energy. The resulting local velocities can be orders of magnitude larger than the injection velocity and they induce damped oscillations of the interface that possess their own time scales and depend only on fluid properties and pore geometry. In complex media (we consider a network of cubic pores) reconfigurations are so frequent and oscillations last long enough that increasing inertial effects leads to a different fluid distribution by influencing the selection of the next pore to be invaded. This calls into question simple pore-filling rules based only on capillary forces. Also, we demonstrate that inertial effects during irreversible reconfigurations can influence the work done by the external forces that is related to the pressure drop in Darcy’s law. This suggests that these phenomena have to be considered when upscaling multiphase flow because local oscillations of the menisci affect macroscopic quantities and modify the constitutive relationships to be used in macro-scale models. These results can be extrapolated to other interface instabilities that are at the origin of fast pore-scale events, such as Haines jumps, snap-off and coalescence. 相似文献
5.
A significant body of current research is aimed at developing methods for numerical simulation of flow and transport in porous media that explicitly resolve complex pore and solid geometries, and at utilizing such models to study the relationships between fundamental pore-scale processes and macroscopic manifestations at larger (i.e., Darcy) scales. A number of different numerical methods for pore-scale simulation have been developed, and have been extensively tested and validated for simplified geometries. However, validation of pore-scale simulations of fluid velocity for complex, three-dimensional (3D) pore geometries that are representative of natural porous media is challenging due to our limited ability to measure pore-scale velocity in such systems. Recent advances in magnetic resonance imaging (MRI) offer the opportunity to measure not only the pore geometry, but also local fluid velocities under steady-state flow conditions in 3D and with high spatial resolution. In this paper, we present a 3D velocity field measured at sub-pore resolution (tens of micrometers) over a centimeter-scale 3D domain using MRI methods. We have utilized the measured pore geometry to perform 3D simulations of Navier–Stokes flow over the same domain using direct numerical simulation techniques. We present a comparison of the numerical simulation results with the measured velocity field. It is shown that the numerical results match the observed velocity patterns well overall except for a variance and small systematic scaling which can be attributed to the known experimental uncertainty in the MRI measurements. The comparisons presented here provide strong validation of the pore-scale simulation methods and new insights for interpretation of uncertainty in MRI measurements of pore-scale velocity. This study also provides a potential benchmark for future comparison of other pore-scale simulation methods. © 2012 Elsevier Science. All rights reserved. 相似文献
6.
7.
Upscaling pore-scale processes into macroscopic quantities such as hydrodynamic dispersion is still not a straightforward matter for porous media with complex pore space geometries. Recently it has become possible to obtain very realistic 3D geometries for the pore system of real rocks using either numerical reconstruction or micro-CT measurements. In this work, we present a finite element–finite volume simulation method for modeling single-phase fluid flow and solute transport in experimentally obtained 3D pore geometries. Algebraic multigrid techniques and parallelization allow us to solve the Stokes and advection–diffusion equations on large meshes with several millions of elements. We apply this method in a proof-of-concept study of a digitized Fontainebleau sandstone sample. We use the calculated velocity to simulate pore-scale solute transport and diffusion. From this, we are able to calculate the a priori emergent macroscopic hydrodynamic dispersion coefficient of the porous medium for a given molecular diffusion Dm of the solute species. By performing this calculation at a range of flow rates, we can correctly predict all of the observed flow regimes from diffusion dominated to convection dominated. 相似文献
8.
Chaotic advection at the pore scale: Mechanisms,upscaling and implications for macroscopic transport
The macroscopic spreading and mixing of solute plumes in saturated porous media is ultimately controlled by processes operating at the pore scale. Whilst the conventional picture of pore-scale mechanical dispersion and molecular diffusion leading to persistent hydrodynamic dispersion is well accepted, this paradigm is inherently two-dimensional (2D) in nature and neglects important three-dimensional (3D) phenomena. We discuss how the kinematics of steady 3D flow at the pore scale generate chaotic advection—involving exponential stretching and folding of fluid elements—the mechanisms by which it arises and implications of microscopic chaos for macroscopic dispersion and mixing. Prohibited in steady 2D flow due to topological constraints, these phenomena are ubiquitous due to the topological complexity inherent to all 3D porous media. Consequently 3D porous media flows generate profoundly different fluid deformation and mixing processes to those of 2D flow. The interplay of chaotic advection and broad transit time distributions can be incorporated into a continuous-time random walk (CTRW) framework to predict macroscopic solute mixing and spreading. We show how these results may be generalised to real porous architectures via a CTRW model of fluid deformation, leading to stochastic models of macroscopic dispersion and mixing which both honour the pore-scale kinematics and are directly conditioned on the pore-scale architecture. 相似文献
9.
10.
Multiphase flow in porous media is described by coupled nonlinear mass conservation laws. For immiscible Darcy flow of multiple fluid phases, whereby capillary effects are negligible, the transport equations in the presence of viscous and buoyancy forces are highly nonlinear and hyperbolic. Numerical simulation of multiphase flow processes in heterogeneous formations requires the development of discretization and solution schemes that are able to handle the complex nonlinear dynamics, especially of the saturation evolution, in a reliable and computationally efficient manner. In reservoir simulation practice, single-point upwinding of the flux across an interface between two control volumes (cells) is performed for each fluid phase, whereby the upstream direction is based on the gradient of the phase-potential (pressure plus gravity head). This upwinding scheme, which we refer to as Phase-Potential Upwinding (PPU), is combined with implicit (backward-Euler) time discretization to obtain a Fully Implicit Method (FIM). Even though FIM suffers from numerical dispersion effects, it is widely used in practice. This is because of its unconditional stability and because it yields conservative, monotone numerical solutions. However, FIM is not unconditionally convergent. The convergence difficulties are particularly pronounced when the different immiscible fluid phases switch between co-current and counter-current states as a function of time, or (Newton) iteration. Whether the multiphase flow across an interface (between two control-volumes) is co-current, or counter-current, depends on the local balance between the viscous and buoyancy forces, and how the balance evolves in time. The sensitivity of PPU to small changes in the (local) pressure distribution exacerbates the problem. The common strategy to deal with these difficulties is to cut the timestep and try again. Here, we propose a Hybrid-Upwinding (HU) scheme for the phase fluxes, then HU is combined with implicit time discretization to yield a fully implicit method. In the HU scheme, the phase flux is divided into two parts based on the driving force. The viscous-driven and buoyancy-driven phase fluxes are upwinded differently. Specifically, the viscous flux, which is always co-current, is upwinded based on the direction of the total-velocity. The buoyancy-driven flux across an interface is always counter-current and is upwinded such that the heavier fluid goes downward and the lighter fluid goes upward. We analyze the properties of the Implicit Hybrid Upwinding (IHU) scheme. It is shown that IHU is locally conservative and produces monotone, physically-consistent numerical solutions. The IHU solutions show numerical diffusion levels that are slightly higher than those for standard FIM (i.e., implicit PPU). The primary advantage of the IHU scheme is that the numerical overall-flux of a fluid phase remains continuous and differentiable as the flow regime changes between co-current and counter-current conditions. This is in contrast to the standard phase-potential upwinding scheme, in which the overall fractional-flow (flux) function is non-differentiable across the boundary between co-current and counter-current flows. 相似文献
11.
Alexander Y. Rozhko 《Geophysical Prospecting》2019,67(5):1404-1430
A key task of exploration geophysics is to find relationships between seismic attributes (velocities and attenuation) and fluid properties (saturation and pore pressure). Experimental data suggest that at least three different factors affect these relationships, which are not well explained by classical Gassmann, Biot, squirt-flow, mesoscopic-flow and gas dissolution/exsolution models. Some of these additional factors include (i) effect of wettability and surface tension between immiscible fluids, (ii) saturation history effects (drainage versus imbibition) and (iii) effects of wave amplitude and effective stress. We apply a new rock physics model to explain the role of all these additional factors on seismic properties of a partially saturated rock. The model is based on a well-known effect in surface chemistry: hysteresis of liquid bridges. This effect is taking place in cracks, which are partially saturated with two immiscible fluids. Using our model, we investigated (i) physical factors affecting empirical Brie correlation for effective bulk modulus of fluid, (ii) the role of liquids on seismic attenuation in the low frequency (static) limit, (iii) water-weakening effects and (iv) saturation history effects. Our model is applicable in the low frequency limit (seismic frequencies) when capillary forces dominate over viscous forces during wave-induced two-phase fluid flow. The model is relevant for the seismic characterization of immiscible fluids with high contrast in compressibilities, that is, for shallow gas exploration and CO2 monitoring. 相似文献
12.
For a high-velocity stable flow through a periodic corrugated channel representing an element of porous medium, we suggest splitting the overall nonlinear macroscopic effects into two kinds of different physical origin: a pure inertia effect produced by the convective term of Navier–Stokes equations and an inertia–viscous cross effect representing a variation of the viscous dissipation due to a streamline deformation by inertia forces. We will show that the inertia–viscous cross effects may be revealed by simulating a periodic flow, whilst the pure inertia effects are produced by the microscale flow nonperiodicity. We will reveal the individual flow law for each nonlinear component and analyze the relative role of both components numerically by using the finite element method applied to the Navier–Stokes equations. Both the pure inertia and the inertia–viscous cross effects are revealed to be exponential prior to quadratic or cubic ones. The influence of the dead volume is analyzed. The inertia–viscous cross phenomena are shown to be negligible when the flow structure is clearly nonperiodic. 相似文献
13.
We present the results of a pore-scale experimental study of residual trapping in consolidated sandstone and carbonate rock samples under confining stress. We investigate how the changes in wetting phase flow rate impacts pore-scale distribution of fluids during imbibition in natural, water-wet porous media. We systematically study pore-scale trapping of the nonwetting phase as well as size and distribution of its disconnected globules. Seven sets of drainage-imbibition experiments were performed with brine and oil as the wetting and nonwetting phases, respectively. We utilized a two-phase miniature core-flooding apparatus integrated with an X-ray microtomography system to examine pore-scale fluid distributions in small Bentheimer sandstone (D = 4.9 mm and L = 13 mm) and Gambier limestone (D = 4.4 mm and L = 75 mm) core samples. The results show that with increase in capillary number, the residual oil saturation at the end of the imbibition reduces from 0.46 to 0.20 in Bemtheimer sandstone and from 0.46 to 0.28 in Gambier limestone. We use pore-scale displacement mechanisms, in-situ wettability characteristics, and pore size distribution information to explain the observed capillary desaturation trends. The reduction was believed to be caused by alteration of the order in which pore-scale displacements took place during imbibition. Furthermore, increase in capillary number produced significantly different pore-scale fluid distributions during imbibition. We explored the pore fluid occupancies and studied size and distribution of the trapped oil clusters during different imbibition experiments. The results clearly show that as the capillary number increases, imbibition produces smaller trapped oil globules. In other words, the volume of individual trapped oil globules decreased at higher brine flow rates. Finally, we observed that the pore space in the limestone sample was considerably altered through matrix dissolution at extremely high brine flow rates. This increased the sample porosity from 44% to 62% and permeability from 7.3 D to 80 D. Imbibition in the altered pore space produced lower residual oil saturation (from 0.28 to 0.22) and significantly different distribution of trapped oil globules. 相似文献
14.
15.
The macroscopic modelling of two-phase flow processes in subsurface hydrosystems or industrial applications on the Darcy scale usually requires a constitutive relationship between capillary pressure and saturation, the Pc(Sw) relationship. Traditionally, it is assumed that a unique relation between Pc and Sw exists independently of the flow conditions as long as hysteretic effects can be neglected. Recently, this assumption has been questioned and alternative formulations have been suggested. For example, the extended Pc(Sw) relationship by Hassanizadeh and Gray [Hassanizadeh SM, Gray WG. Mechanics and thermodynamics of multiphase flow in porous media including interphase boundaries. Adv Water Resources 1990;13(4):169–86] proposes that the difference between the phase pressures to the equilibrium capillary pressure is a linear function of the rate of change of saturation, thereby introducing a constant of proportionality, the coefficient τ. It is desirable to identify cases where the extended relationship needs to be considered. Consequently, a dimensional analysis is performed on the basis of the two-phase balance equations. In addition to the well-known capillary and gravitational number, the dimensional analysis yields a new dimensionless number. The dynamic number Dy quantifies the ratio of dynamic capillary to viscous forces. Relating the dynamic to the capillary as well as the gravitational number gives the new numbers DyC and DyG, respectively. For given sets of fluid and porous medium parameters, the dimensionless numbers Dy and DyC are interpreted as functions of the characteristic length and flow velocity. The simulation of an imbibition process provides insight into the interpretation of the characteristic length scale. The most promising choice for this length scale seems to be the front width. We conclude that consideration of the extended Pc(Sw) relationship may be important for porous media with high permeability, small entry pressure and high coefficient τ when systems with a small characteristic length (e.g. steep front) and small characteristic time scale are under investigation. 相似文献
16.
Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based on the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This study provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods. 相似文献
17.
Sizeable amounts of connected microporosity with various origins can have a profound effect on important petrophysical properties of a porous medium such as (absolute/relative) permeability and capillary pressure relationships. We construct pore-throat networks that incorporate both intergranular porosity and microporosity. The latter originates from two separate mechanisms: partial dissolution of grains and pore fillings (e.g. clay). We then use the reconstructed network models to estimate the medium flow properties. In this work, we develop unique network construction algorithms and simulate capillary pressure–saturation and relative permeability–saturation curves for cases with inhomogeneous distributions of pores and micropores. Furthermore, we provide a modeling framework for variable amounts of cement and connectivity of the intergranular porosity and quantifying the conditions under which microporosity dominates transport properties. In the extreme case of a disconnected inter-granular network due to cementation a range of saturations within which neither fluid phase is capable of flowing emerges. To our knowledge, this is the first flexible pore scale model, from first principles, to successfully approach this behavior observed in tight reservoirs. 相似文献
18.
《Advances in water resources》2002,25(1):33-49
Experiments designed to elucidate the pore-scale mechanisms of the dissolution of a residual non-aqueous phase liquid (NAPL), trapped in the form of ganglia within a porous medium, are discussed. These experiments were conducted using transparent glass micromodels with controlled pore geometry, so that the evolution of the size and shape of individual NAPL ganglia and, hence, the pore-scale mass transfer rates and mass transfer coefficients could be determined by image analysis. The micromodel design permitted reasonably accurate control of the pore water velocity, so that the mass transfer coefficients could be correlated in terms of a local (pore-scale) Peclet number. A simple mathematical model, incorporating convection and diffusion in a slit geometry was developed and used successfully to predict the observed mass transfer rates. For the case of non-wetting NAPL ganglia, water flow through the corners in the pore walls was seen to control the rate of NAPL dissolution, as recently postulated by Dillard and Blunt [Water Resour. Res. 36 (2000) 439–454]. Break-up of doublet non-wetting phase ganglia into singlet ganglia by snap-off in pore throats was also observed, confirming the interplay between capillarity and mass transfer. Additionally, the effect of wettability on dissolution mass transfer was demonstrated. Under conditions of preferential NAPL wettability, mass transfer from NAPL films covering the solid surfaces was seen to control the dissolution process. Supply of NAPL from the trapped ganglia to these films by capillary flow along pore corners was observed to result in a sequence of pore drainage events that increase the interfacial area for mass transfer. These observations provide new experimental evidence for the role of capillarity, wettability and corner flow on NAPL ganglia dissolution. 相似文献
19.
Two-phase imbibition behavior of immiscible fluids was studied in dry and prewetted porous media using a laser-induced fluorescence technique. Imbibition was first investigated in two-dimensional (2-D) systems under conditions comparable to those for a study of drainage [Ovdat H, Berkowitz B. Pore-scale study of drainage displacement under combined capillary and gravity effects in index-matched porous media. Water Resources Research 2006;42:W06411. doi:10.1029/2005WR004553] in the capillary-dominated regime. The effect of initial wetting saturation (IWS) was then explored in 2-D and 3-D porous media under the combined effect of gravity, capillary and viscous forces, within and outside the capillary-dominated regime. Parameters that describe maximum vertical advance, volumetric fraction, total surface area and specific surface area of the invading fluid were used to quantify the behavior. Comparison of 2-D drainage and imbibition patterns demonstrates significant qualitative differences under analogous viscosity ratio, buoyancy number, and capillary number values. However, quantitative analyses show strong pore-scale similarities between these patterns. Invasion structures in 3-D, prewetted (IWS ≈ 8% of the pore volume) porous media are ramified, with lateral branching and regions containing trapped residual fluid. These structures are qualitatively and quantitatively different from the compact, branchless structures that develop in dry (IWS = 0) porous media. 相似文献
20.
Motivated by a wide range of applications from enhanced oil recovery to carbon dioxide sequestration, we have developed a two-dimensional, pore-level model of immiscible drainage, incorporating viscous, capillary, and gravitational effects. This model has been validated quantitatively, in the very different limits of zero viscosity ratio and zero capillary number; flow patterns from modeling agree well with experiment. For a range of stable viscosity ratios (μinjected/μdisplaced ? 1), we have increased the capillary number, Nc, and studied the way in which the flows deviate from capillary fingering (the fractal flow of invasion percolation) and become compact for realistic capillary numbers. Results exhibiting this crossover from capillary fingering to compact invasion are presented for the average position of the injected fluid, the fluid–fluid interface, the saturation and fractional flow profiles, and the relative permeabilities. The agreement between our results and earlier theoretical predictions [Blunt M, King MJ, Scher H. Simulation and theory of two-phase flow in porous media. Phys Rev A 1992;46:7680–99; Lenormand R. Flow through porous media: limits of fractal patterns. Proc Roy Soc A 1989;423:159–68; Wilkinson D. Percolation effects in immiscible displacement. Phys Rev A 1986;34:1380–90; Xu B, Yortsos YC, Salin D. Invasion Percolation with viscous forces. Phys Rev E 1998;57:739–51] supports the validity of these general theoretical arguments, which were independent of the details of the porous media in both two and three dimensions. 相似文献