共查询到20条相似文献，搜索用时 734 毫秒
1.
We have conducted a series of high-resolution numerical experiments using the Pair-Wise Force Smoothed Particle Hydrodynamics (PF-SPH) multiphase flow model. First, we derived analytical expressions relating parameters in the PF-SPH model to the surface tension and static contact angle. Next, we used the model to study viscous fingering, capillary fingering, and stable displacement of immiscible fluids in porous media for a wide range of capillary numbers and viscosity ratios. We demonstrated that the steady state saturation profiles and the boundaries of viscous fingering, capillary fingering, and stable displacement regions compare favorably with micromodel laboratory experimental results. For a displacing fluid with low viscosity, we observed that the displacement pattern changes from viscous fingering to stable displacement with increasing injection rate. When a high viscosity fluid is injected, transition behavior from capillary fingering to stable displacement occurred as the flow rate was increased. These observations are also in agreement with the results of the micromodel laboratory experiments. 相似文献
2.
— We studied experimentally the stability of the front between two immiscible fluids in a Hele-Shaw cell, which consists of two parallel glass plates separated by a narrow spacing. The fluid flow in a Hele-Shaw cell is governed by the Laplace's equation (pressure field) and Darcy's law (fluid transport), which makes of this experimental setup a two-dimensional analog of a porous medium (FABER, 1997). SAFFMAN and TAYLOR (1958) demonstrated that when the injected fluid is less viscous than the displaced one, the front between the fluids can be unstable and fingers may develop along the interface. We investigated the possibility of influencing the development of the Saffman-Taylor instability by modulating the injection pressure with weak periodic oscillations of low frequency (5 to 15 Hz). A series of experiments were performed, in which the evolution of the interface between the injected air and the displaced silicon oil was recorded with the help of a numerical photo camera. We observed that the oscillations of the injection pressure inhibit the development of fingers. We quantified the observed stabilization effect by performing a harmonic analysis of the shape of the front, and suggested an interpretation in terms of “stiffening” of the interface. 相似文献
3.
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. 相似文献
4.
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. 相似文献
5.
Stochastic analysis is commonly used to address uncertainty in the modeling of flow and transport in porous media. In the stochastic approach, the properties of porous media are treated as random functions with statistics obtained from field measurements. Several studies indicate that hydrological properties depend on the scale of measurements or support scales, but most stochastic analysis does not address the effects of support scale on stochastic predictions of subsurface processes. In this work we propose a new approach to study the scale dependence of stochastic predictions. We present a stochastic analysis of immiscible fluid–fluid displacement in randomly heterogeneous porous media. While existing solutions are applicable only to systems in which the viscosity of one phase is negligible compare with the viscosity of the other (water–air systems for example), our solutions can be applied to the immiscible displacement of fluids having arbitrarily viscosities such as NAPL–water and water–oil. Treating intrinsic permeability as a random field with statistics dependant on the permeability support scale (scale of measurements) we obtained, for one-dimensional systems, analytical solutions for the first moments characterizing unbiased predictions (estimates) of system variables, such as the pressure and fluid–fluid interface position, and we also obtained second moments, which characterize the uncertainties associated with such predictions. Next we obtained empirically scale dependent exponential correlation function of the intrinsic permeability that allowed us to study solutions of stochastic equations as a function of the support scale. We found that the first and second moments converge to asymptotic values as the support scale decreases. In our examples, the statistical moments reached asymptotic values for support scale that were approximately 1/10000 of the flow domain size. We show that analytical moment solutions compare well with the results of Monte Carlo simulations for moderately heterogeneous porous media, and that they can be used to study the effects of heterogeneity on the dynamics and stability of immiscible flow. 相似文献
6.
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. 相似文献
7.
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. 相似文献
8.
Wavefields in porous media saturated by two immiscible fluids are simulated in this paper. Based on the sealed system theory,
the medium model considers both the relative motion between the fluids and the solid skeleton and the relaxation mechanisms
of porosity and saturation (capillary pressure). So it accurately simulates the numerical attenuation property of the wavefields
and is much closer to actual earth media in exploration than the equivalent liquid model and the unsaturated porous medium
model on the basis of open system theory. The velocity and attenuation for different wave modes in this medium have been discussed
in previous literature but studies of the complete wave-field have not been reported. In our work, wave equations with the
relaxation mechanisms of capillary pressure and the porosity are derived. Furthermore, the wavefield and its characteristics
are studied using the numerical finite element method. The results show that the slow P3-wave in the non-wetting phase can
be observed clearly in the seismic band. The relaxation of capillary pressure and the porosity greatly affect the displacement
of the non-wetting phase. More specifically, the displacement decreases with increasing relaxation coefficient. 相似文献
9.
《Advances in water resources》1998,21(6):451-461
In this paper we consider one-dimensional capillary redistribution of two immiscible and incompressible fluids in a porous medium with a single discontinuity. We study a special time-dependent solution, a similarity solution, which is found when the initial saturation is discontinuous at the same point as the permeability and porosity, and is constant elsewhere. The similarity solution can be used to validate numerical algorithms describing two-phase flow in porous media with discontinuous heterogeneities. We discuss the construction of the similarity solution, in which we pay special attention to the interface conditions at the discontinuity, both for media with positive and zero entry pressure. Moreover, we discuss some qualitative properties of the solution, and outline a numerical procedure to determine its graph. Examples are given for the Brooks-Corey and Van Genuchten model. We also consider similarity solutions for unsaturated water flow, which is a limit case of two-phase flow for negligible nonwetting phase viscosity. 相似文献
10.
地下岩石由岩石骨架和孔隙流体组成,通常流体含黏性.地震波在地下介质中传播时受岩石骨架和黏性流体的影响会呈现出复杂的变化.本文将流、固体位移和应力连续作为边界条件,推导出含黏性流体孔隙介质分界面上反透射系数方程;通过建立上层为饱油、下层为饱盐水的砂岩孔隙介质模型,开展反透射系数特征研究,分别分析不同频率、不同黏滞系数条件下,含黏性流体孔隙介质分界面上反透射系数随入射角的变化.研究表明,孔隙介质分界面上和等效介质分界面上的反透射系数分别随入射角的变化趋势基本一致,说明方程推导和数值计算的正确性;快纵波反透射系数受频率、流体黏性的影响较小,而快横波反透射系数在一定入射角范围内受频率、流体黏性的影响比较大;由于黏性孔隙流体的作用,慢纵波和慢横波的反透射系数受入射角、频率及流体黏性的影响都很大. 相似文献
11.
Multiphase flow in capillary regimes is a fundamental process in a number of geoscience applications. The ability to accurately define wetting characteristics of porous media can have a large impact on numerical models. In this paper, a newly developed automated three-dimensional contact angle algorithm is described and applied to high-resolution X-ray microtomography data from multiphase bead pack experiments with varying wettability characteristics. The algorithm calculates the contact angle by finding the angle between planes fit to each solid/fluid and fluid/fluid interface in the region surrounding each solid/fluid/fluid contact point. Results show that the algorithm is able to reliably compute contact angles using the experimental data. The in situ contact angles are typically larger than flat surface laboratory measurements using the same material. Wetting characteristics in mixed-wet systems also change significantly after displacement cycles. 相似文献
12.
We report here on recent developments and advances in pore-scale X-ray tomographic imaging of subsurface porous media. Our particular focus is on immiscible multi-phase fluid flow, i.e., the displacement of one immiscible fluid by another inside a porous material, which is of central importance to many natural and engineered processes. Multiphase flow and displacement can pose a rather difficult problem, both because the underlying physics is complex, and also because standard laboratory investigation reveals little about the mechanisms that control micro-scale processes. X-ray microtomographic imaging is a non-destructive technique for quantifying these processes in three dimensions within individual pores, and as we report here, with rapidly increasing spatial and temporal resolution. 相似文献
13.
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. 相似文献
14.
The flow of two immiscible fluids through a porous medium depends on the complex interplay between gravity, capillarity, and viscous forces. The interaction between these forces and the geometry of the medium gives rise to a variety of complex flow regimes that are difficult to describe using continuum models. Although a number of pore-scale models have been employed, a careful investigation of the macroscopic effects of pore-scale processes requires methods based on conservation principles in order to reduce the number of modeling assumptions. In this work we perform direct numerical simulations of drainage by solving Navier–Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid–fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and model the transition from stable flow to viscous fingering, we focus on the macroscopic capillary pressure and we compare different definitions of this quantity under quasi-static and dynamic conditions. We show that the difference between the intrinsic phase-average pressures, which is commonly used as definition of Darcy-scale capillary pressure, is subject to several limitations and it is not accurate in presence of viscous effects or trapping. In contrast, a definition based on the variation of the total surface energy provides an accurate estimate of the macroscopic capillary pressure. This definition, which links the capillary pressure to its physical origin, allows a better separation of viscous effects and does not depend on the presence of trapped fluid clusters. 相似文献
15.
A dominant mechanism for residual trapping of a nonwetting fluid in porous media during imbibition is snap-off or the disconnection of a continuous stream of the nonwetting fluid when it passes through pore constrictions and when a criterion based on capillary pressure imbalance is met. While quasi-static criteria for Roof snap-off have been defined for pores based on the imbalance between capillary pressure across the front/tail meniscus and local capillary pressure at the pore throat, and expressed in terms of pore body to pore throat ratio for simplification, we extended the previous quasi-static snap-off criterion by considering the local capillary pressure imbalance between the pore body and the pore throat for both circular and noncircular pores when the wetting film exists. We then used the criterion to analyze results from computational fluid dynamics (CFD) simulations of multi-phase flow with supercritical CO2 as the nonwetting fluid and water as the wetting fluid. The extended criterion successfully described most situations we modeled. Furthermore, we compared fluid interface shape for a noncircular 3D pore predicted by the minimum surface energy (MSE) theory against 3D CFD simulations. While the fluid interface shape at the pore throat for 3D simulation was consistent with the shape predicted by MSE theory, the shape could not be successfully predicted by the MSE theory at the upstream and downstream pore body. Moreover, film flow existed for the noncircular pore at the downstream pore body. 相似文献
16.
Contrast in capillary pressure of heterogeneous permeable media can have a significant effect on the flow path in two-phase immiscible flow. Very little work has appeared on the subject of capillary heterogeneity despite the fact that in certain cases it may be as important as permeability heterogeneity. The discontinuity in saturation as a result of capillary continuity, and in some cases capillary discontinuity may arise from contrast in capillary pressure functions in heterogeneous permeable media leading to complications in numerical modeling. There are also other challenges for accurate numerical modeling due to distorted unstructured grids because of the grid orientation and numerical dispersion effects. Limited attempts have been made in the literature to assess the accuracy of fluid flow modeling in heterogeneous permeable media with capillarity heterogeneity. The basic mixed finite element (MFE) framework is a superior method for accurate flux calculation in heterogeneous media in comparison to the conventional finite difference and finite volume approaches. However, a deficiency in the MFE from the direct use of fractional flow formulation has been recognized lately in application to flow in permeable media with capillary heterogeneity. In this work, we propose a new consistent formulation in 3D in which the total velocity is expressed in terms of the wetting-phase potential gradient and the capillary potential gradient. In our formulation, the coefficient of the wetting potential gradient is in terms of the total mobility which is smoother than the wetting mobility. We combine the MFE and discontinuous Galerkin (DG) methods to solve the pressure equation and the saturation equation, respectively. Our numerical model is verified with 1D analytical solutions in homogeneous and heterogeneous media. We also present 2D examples to demonstrate the significance of capillary heterogeneity in flow, and a 3D example to demonstrate the negligible effect of distorted meshes on the numerical solution in our proposed algorithm. 相似文献
17.
M. R. Deinert J. -Y. Parlange T. Steenhuis J. Throop K. Ünlü K. B. Cady 《Journal of Hydrology》2004,290(3-4):192-201
Few methods exist for measuring rapidly changing fluid contents at the pore scale that simultaneously allow whole flow field visualization. We present a method for using real-time neutron radiography to measure rapidly changing moisture profiles in porous media. The imaging technique monitors the attenuation of a thermal neutron beam as it traverses a flow field and provides measurements every 30 ms with an image area >410 cm2 and a spatial resolution 0.05 cm. The technique is illustrated by measuring the variation in moisture content across a wetting front moving at constant velocity through SiO2 sand. The relative contributions of the hydraulic conductivity and diffusivity terms in Richards' equation to the total fluid flux within the wetting front region were also measured. The diffusivity was found to rise from zero to a peak value within the wetting front region before falling off while the conductivity was found to rise monotonically. The reliability of the technique was checked via mass balance. 相似文献
18.
The phenomenon of reflection and transmission of plane harmonic waves at the plane interface between two dissimilar poroelastic solids saturated with two immiscible viscous fluids is investigated. Both porous media are considered dissipative due to the presence of viscosity in pore‐fluids. Four attenuated (three dilatational and one shear) waves propagate in such a dissipative porous medium. A finite non‐dimensional parameter is used to define the effective connections between the surface‐pores of two media at their common interface. Another finite parameter represents the gas‐share in the saturation of pores. An attenuated wave in a dissipative medium is described through the specification of directions of propagation and maximum attenuation. A general representation of an attenuated wave is defined through its inhomogeneous propagation, i.e., different directions for propagation and attenuation. Incidence of an inhomogeneous wave is considered at the interface between two dissipative porous solids. This results in four reflected and four transmitted inhomogeneous waves. Expressions are derived for the partition of incident energy among the reflected and transmitted waves. Numerical examples are studied to determine the effects of saturating pore fluid, frequency, surface‐pore connections and wave inhomogeneity on the strengths of reflected and transmitted waves. Interaction energy due to the interference of different (inhomogeneous) waves is calculated in both the dissipative porous media to verify the conservation of incident energy. 相似文献
19.
《Advances in water resources》1998,21(5):351-362
Flow and displacement of non-Newtonian fluids in porous media occurs in many subsurface systems, related to underground natural resource recovery and storage projects, as well as environmental remediation schemes. A thorough understanding of non-Newtonian fluid flow through porous media is of fundamental importance in these engineering applications. Considerable progress has been made in our understanding of single-phase porous flow behavior of non-Newtonian fluids through many quantitative and experimental studies over the past few decades. However, very little research can be found in the literature regarding multi-phase non-Newtonian fluid flow or numerical modeling approaches for such analyses.For non-Newtonian fluid flow through porous media, the governing equations become nonlinear, even under single-phase flow conditions, because effective viscosity for the non-Newtonian fluid is a highly nonlinear function of the shear rate, or the pore velocity. The solution for such problems can in general only be obtained by numerical methods.We have developed a three-dimensional, fully implicit, integral finite difference simulator for single- and multi-phase flow of non-Newtonian fluids in porous/fractured media. The methodology, architecture and numerical scheme of the model are based on a general multi-phase, multi-component fluid and heat flow simulator — TOUGH2. Several rheological models for power-law and Bingham non-Newtonian fluids have been incorporated into the model. In addition, the model predictions on single- and multi-phase flow of the power-law and Bingham fluids have been verified against the analytical solutions available for these problems, and in all the cases the numerical simulations are in good agreement with the analytical solutions. In this presentation, we will discuss the numerical scheme used in the treatment of non-Newtonian properties, and several benchmark problems for model verification.In an effort to demonstrate the three-dimensional modeling capability of the model, a three-dimensional, two-phase flow example is also presented to examine the model results using laboratory and simulation results existing for the three-dimensional problem with Newtonian fluid flow. 相似文献
20.
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. 相似文献