Horizontal redistribution of fluids in a porous medium: The role of interfacial area in modeling hysteresis

Recent advances in multi-phase flow theory have shown that the flow of several phases in a porous medium is highly influenced by the interfaces separating these phases. First modeling studies based on this new theory have been performed on a pore scale, as well as on a volume-averaged macro scale using balance equations and constitutive relations that take the role and presence of interfaces into account. However, neither experimental data nor analytical solutions are available on the macro scale so far, although their knowledge is essential for the verification of the new models.     

Thermodynamically constrained averaging theory approach for modeling flow and transport phenomena in porous medium systems: 3. Single-fluid-phase flow

This work is the third in a series of papers on the thermodynamically constrained averaging theory (TCAT) approach to modeling flow and transport phenomena in multiscale porous medium systems. Building upon the general TCAT framework and the mathematical foundation presented in previous works in this series, we demonstrate the TCAT approach for the case of single-fluid-phase flow. The formulated model is based upon conservation equations for mass, momentum, and energy and a general entropy inequality constraint, which is developed to guide model closure. A specific example of a closed model is derived under limiting assumptions using a linearization approach and these results are compared and contrasted with the traditional single-phase-flow model. Potential extensions to this work are discussed. Specific advancements in this work beyond previous averaging theory approaches to single-phase flow include use of macroscale thermodynamics that is averaged from the microscale, the use of derived equilibrium conditions to guide a flux–force pair approach to simplification, use of a general Lagrange multiplier approach to connect conservation equation constraints to the entropy inequality, and a focus on producing complete, closed models that are solvable.     

Pore-scale characteristics of multiphase flow in porous media: A comparison of air–water and oil–water experiments

Studies of NAPL dissolution in porous media have demonstrated that measurement of saturation alone is insufficient to describe the rate of dissolution. Quantification of the NAPL–water interfacial area provides a measure of the expected area available for mass transfer and will likely be a primary determinant of NAPL removal efficiency. To measure the interfacial area, we have used a synchrotron-based CMT technique to obtain high-resolution 3D images of flow in a Soltrol–water–glass bead system. The interfacial area is found to increase as the wetting phase saturation decreases, reach a maximum, and then decrease as the wetting phase saturation goes to zero. These results are compared to previous findings for an air–water–glass bead study; The Soltrol–water interfacial areas were found to peak at similar saturations as those measured for the air–water system (20–35% saturation range), however, the peak values were in some cases almost twice as high for the oil-water system. We believe that the observed differences between the air–water and oil–water systems to a large degree can be explained by the differences in interfacial tensions for the two systems.     

One-dimensional dispersion in unsteady flow in an adsorbing porous medium: An analytical solution

Analytical solutions describing the concentration distribution along one-dimensional unsteady seepage flow through adsorbing saturated finite porous medium have been obtained. An exponential function concentration is enforced at the source of the dispersion, while the change in the concentration is zero at the other boundary. A new time variable has been introduced to solve the unsteady flow problem and the solution is illustrated.     

Lattice-Boltzmann simulations of the capillary pressure–saturation–interfacial area relationship for porous media

Hysteresis in the relationship between capillary pressure (P_c), wetting phase saturation (S_w) and nonwetting–wetting interfacial area per volume (a_nw) is investigated using multiphase lattice-Boltzmann simulations of drainage and imbibition in a glass bead porous system. In order to validate the simulations, the P_c–S_w and a_nw–S_w main hysteresis loops were compared to experimental data reported by Culligan et al. [Culligan KA, Wildenschild D, Christensen BS, Gray WG, Rivers ML, Tompson AB. Interfacial area measurements for unsaturated flow through porous media. Water Resour Res 2004;40:W12413]. In general, the comparison shows that the simulations are reliable and capture the important physical processes in the experimental system. P_c–S_w curves, a_nw–S_w curves and phase distributions (within the pores) show good agreement during drainage, but less satisfactory agreement during imbibition. Drainage and imbibition scanning curves were simulated in order to construct P_c–S_w–a_nw surfaces. The root mean squared error (RMSE) and mean absolute error (MAE) between drainage and imbibition surfaces was 0.10 mm⁻¹ and 0.03 mm⁻¹, respectively. This small difference indicates that hysteresis is virtually nonexistent in the P_c–S_w–a_nw relationship for the multiphase system studied here. Additionally, a surface was fit to the main loop (excluding scanning curves) of the drainage and imbibition P_c–S_w–a_nw data and compared to the surface fit to all of the data. The differences between these two surfaces were small (RMSE = 0.05 mm⁻¹ and MAE = 0.01 mm⁻¹) indicating that the P_c–S_w–a_nw surface is adequately represented without the need for the scanning curve data, which greatly reduces the amount of data required to construct the non-hysteretic P_c–S_w–a_nw surface for this data.     

随机-均匀介质散射波特征相干分析研究

基于统计学方法建立了随机-均匀介质模型,利用有限差分法进行正演模拟,根据地震剖面中散射波同相轴随着模型横向自相关长度的增大其连续性逐渐变好的特点,通过相干技术研究了散射波同相轴的连续性.研究结果表明,当模型横向自相关长度与地震波长比值较小（a/λ<1）时,散射波同相轴连续性较差,其相干值变化不明显;当模型横向自相关长度与地震波长比值较大（a/λ>1）时,散射波同相轴连续性明显变好,其相干值明显增大.     

饱和多孔微极介质的波动方程及其势函数方程

土是由一定尺寸大小颗粒所构成的多孔介质,具有明显的颗粒特性,当土颗粒间的孔隙被流体（如水或油）充满时则成为饱和土.利用微极理论和Biot波动理论的研究成果,把饱和土中多孔固体骨架部分近似地视为微极介质,孔隙中的流体部分视为质点介质,获得饱和多孔微极介质的弹性波动方程.借鉴Greetsma理论,建立了饱和多孔微极介质弹性本构方程力学参数与相应单相介质弹性参数的相互关系,使饱和多孔微极介质弹性波动方程中的物理参数具有明确的物理意义,易于在试验中确定.运用场论理论把饱和多孔微极介质的波动方程简化为势函数方程,建立了饱和多孔微极介质中五种弹性波的弥散方程,数值分析了五种简谐体波在无限饱和多孔微极介质中的传播特性. 结果表明,P1波、P2波和剪切S1波的波速弥散曲线与经典饱和多孔介质基本相同,当频率小于临界频率ω0时旋转纵波θ波和横波S2波不存在,当频率大于临界频率ω0时,θ波和S2波的传播速度随频率增加而减小.     

Non-Fickian mass transport in fractured porous media

The paper provides an introduction to fundamental concepts of mathematical modeling of mass transport in fractured porous heterogeneous rocks. Keeping aside many important factors that can affect mass transport in subsurface, our main concern is the multi-scale character of the rock formation, which is constituted by porous domains dissected by the network of fractures. Taking into account the well-documented fact that porous rocks can be considered as a fractal medium and assuming that sizes of pores vary significantly (i.e. have different characteristic scales), the fractional-order differential equations that model the anomalous diffusive mass transport in such type of domains are derived and justified analytically. Analytical solutions of some particular problems of anomalous diffusion in the fractal media of various geometries are obtained. Extending this approach to more complex situation when diffusion is accompanied by advection, solute transport in a fractured porous medium is modeled by the advection-dispersion equation with fractional time derivative. In the case of confined fractured porous aquifer, accounting for anomalous non-Fickian diffusion in the surrounding rock mass, the adopted approach leads to introduction of an additional fractional time derivative in the equation for solute transport. The closed-form solutions for concentrations in the aquifer and surrounding rocks are obtained for the arbitrary time-dependent source of contamination located in the inlet of the aquifer. Based on these solutions, different regimes of contamination of the aquifers with different physical properties can be readily modeled and analyzed.     

基于时空守恒元和解元(CE/SE)方法的孔隙介质多相流动计算

将时空守恒元/解元(CE/SE)方法推广到二维孔隙介质多相流问题的数值计算中,采用人工压缩法耦合速度和压力,同时结合杂交粒子水平集方法捕捉物质界面.提出一套完整的二维欧拉型孔隙介质非稳态多相不可压缩黏性流动计算方案.通过对溃坝和液滴在重力作用下的运动和变形问题的数值模拟,验证了方法的精度和有效性.在此基础上,提出了一个新的孔隙介质两相流物理模型——双层流体顶盖驱动方腔流.     

Simulation of density-driven flow in fractured porous media

We study density-driven flow in a fractured porous medium in which the fractures are represented as manifolds of reduced dimensionality. Fractures are assumed to be thin regions of space filled with a porous material whose properties differ from those of the porous medium enclosing them. The interfaces separating the fractures from the embedding medium are assumed to be ideal. We consider two approaches: (i) the fractures have the same dimension, d, as the embedding medium and are said to be d-dimensional; (ii) the fractures are considered as (d − 1)-dimensional manifolds, and the equations of density-driven flow are found by averaging the d-dimensional laws over the fracture width. We show that the second approach is a valid alternative to the first one. For this purpose, we perform numerical experiments using finite-volume discretization for both approaches. The results obtained by the two methods are in good agreement with each other.     

含混合裂隙、孔隙介质的纵波衰减规律研究

地下多孔介质中的孔隙类型复杂多样,既有硬孔又有扁平的软孔.针对复杂孔隙介质,假设多孔介质中同时含有球型硬孔和两种不同产状的裂隙(硬币型、尖灭型裂隙),当孔隙介质承载载荷时,考虑两种不同类型的裂隙对于孔隙流体压力的影响,建立起Biot理论框架下饱和流体情况含混合裂隙、孔隙介质的弹性波动方程,并进一步求取了饱和流体情况下仅由裂隙引起流体流动时的含混合裂隙、孔隙介质的体积模量和剪切模量,随后,在此基础上讨论了含混合裂隙、孔隙介质在封闭条件下地震波衰减和频散的高低频极限表达式.最后计算了给定模型的地震波衰减和频散,发现地震波衰减曲线呈现"多峰"现象,速度曲线为"多频段"频散.针对该模型分析讨论了渗透率参数、裂隙纵横比参数以及流体黏滞性参数对于地震波衰减和频散的影响,表明三个参数均为频率控制参数.     

含两项不混合/单项黏性流体孔隙介质分界面上反透射系数特征研究

弹性孔隙介质分界面上的反透射系数特征,在岩性划分、流体识别、储层边界判识等方面有重要的应用.本文研究上层为含两项不混合黏性流体孔隙介质、下层为含单项黏性流体孔隙介质分界面上的反透射理论.首先根据两种孔隙介质分界面上的能量守恒得到边界条件,再将波函数、位移、应力与应变关系代入边界条件,推导出完全连通孔隙情况下,第一类纵波入射到孔隙介质分界面上的反透射系数方程.通过建立砂岩孔隙介质模型,分别分析不同孔隙流体类型、不同含油饱和度及不同入射角情况下,各类波的反透射系数特征.研究表明,第二、三类纵波反透射系数数值比第一类纵波小多个数量级,且两者对入射角的变化不敏感,但对孔隙流体性质、含油饱和度的变化较敏感,而横波反透射系数特征恰好与此相反;第一类纵波反透射系数特征比较复杂,入射角、孔隙流体的性质及含油饱和度的变化都对其产生影响.不同孔隙流体弹性物性的差异、孔隙介质中含油饱和度的变化及不同入射角引起垂向和切向应力分量的变化都会影响各类波的反透射系数特征,分析这些特征可以为研究储层含油气性提供理论基础.

     

Comment on the effect of anisotropy on the onset of convection in a porous medium

The effect of anisotropy on the onset of convection in a saturated porous medium is discussed. In particular, the case of time-dependent density-driven convection is examined. The applicability of the value of an equivalent Rayleigh number as the criterion for the onset of convection is discussed.     

Non-equilibrium interphase heat and mass transfer during two-phase flow in porous media—Theoretical considerations and modeling

Mass and heat transfer occurring across phase-interfaces in multi-phase flow in porous media are mostly approximated using equilibrium relationships or empirical kinetic models. However, when the characteristic time of flow is smaller than that of mass or heat transfer, non-equilibrium situations may arise. Commonly, empirical approaches are used in such cases. There are only few works in the literature that use physically-based models for these transfer terms. In fact, one would expect physical approaches to modeling kinetic interphase mass and heat transfer to contain the interfacial area between the phases as a variable. Recently, a two-phase flow and solute transport model was developed that included interfacial area as a state variable [36]. In that model, interphase mass transfer was modeled as a kinetic process.     

含黏滞流体各向异性孔隙介质中弹性波的频散和衰减

基于Biot理论,考虑液相的黏弹性变形和固液相接触面上的相对扭转,提出了含黏滞流体VTI孔隙介质模型.从理论上推导出,在该模型中除存在快P波、慢P波、SV波、SH波以外,还将存在两种新横波－慢SV波和慢SH波.数值模拟分析了6种弹性波的相速度、衰减、液固相振幅比随孔隙度、频率的变化规律以及快P波、快SV波的衰减随流体性质、渗透率、入射角的变化规律.结果表明慢SV波和慢SH波主要在液相中传播,高频高孔隙度时,速度较高;大角度入射时,快P波衰减表现出明显的各向异性,而快SV波的衰减则基本不变;储层纵向和横向渗透率存在差异时,快SV波衰减大的方向渗透率高.     

Experimental investigation of suspended particles transport in porous medium under variable temperatures

This paper presents a study on suspended particle transport in porous medium with the aid of a sand layer transportation–deposition testing system to determine the kinetic characteristics of particles in porous medium under variable temperatures. Quartz sand and quartz powder were chosen as the porous medium and particle in the tests, respectively. Four size compositions and two operational modes, that is, temperature reduction mode (changing from 18°C to 5°C) and temperature increment mode (changing from 18°C to 35°C), were adopted. The turbidity and concentration of quartz powder were measured under various conditions. We observed a high temperature‐independent correlation between them. Breakthrough curves under different conditions were analysed using this testing system. The results showed that changes in temperature affected the particle transport process to some extent, and the degree of influence was closely related to the time moment of the temperature change onset. Moreover, we found a hysteresis phenomenon in the breakthrough curve under both temperature reduction and increment conditions. The results also indicated that the temperature effect was particularly significant for smaller particles. The typical curves to represent particle transport process under variable temperatures were put forward according to the results. To explain the test results, four factors, that is, water viscosity, adsorption effect, double layer force, and particle kinetic energy, were considered and categorized as promotion or constraining factors.     

孔、裂隙并存地层中的声波测井理论及多极子声场特征

在油、气储层的勘探和开发中观察到的一个现象是储层岩石中普遍存在孔隙和裂隙.随着近年来孔、裂隙介质弹性波动理论的进展,我们可以将此理论应用于测井技术,以此来指导从声波测井中测量孔、裂隙地层的声学参数.本文计算了孔、裂隙地层里充流体井眼中的多极子声场,分析了声场随裂隙介质的两个主要参数(即裂隙密度和裂隙纵横比)的变化特征.井孔声场的数值计算表明裂隙密度可以大幅度地降低井中声波纵、横波的波速和振幅.随着裂隙密度的增加,在测井频段内也可以看到纵、横波速的频散现象(这种频散在孔隙地层中一般是观察不到的).本文还研究了多极子模式波 (即单极的Stoneley波、伪瑞利波以及偶极的弯曲波)随裂隙参数的变化特征.结果表明,这些模式波的振幅激发和速度频散都受裂隙密度的影响.裂隙密度越高影响越大.此外,裂隙还对模式波的传播造成较大的衰减.相对裂隙密度而言,裂隙纵横比是一个频率控制参数,它控制裂隙对声场影响的频率区间.本文的分析结果对裂缝、孔隙型地层的声波测井具有指导意义.     

Extraction of physically realistic pore network properties from three-dimensional synchrotron X-ray microtomography images of unconsolidated porous media systems

This paper presents application of a series of algorithms used to extract pore network structure from high-resolution three-dimensional synchrotron microtomography images of unconsolidated porous media systems. These algorithms are based on the three-dimensional skeletonization that simplifies the pore space to networks in the form of nodes connected to paths. Dilation algorithms were developed to generate inscribed spheres on the nodes and paths of the medial axis to represent pore-bodies and pore-throats of the network, respectively. The end result is a physically representative pore network structure, i.e. three-dimensional spatial distribution (i.e. x-, y-, and z-coordinates) of pore-bodies and pore-throats, pore-body size distribution, pore-throat size distribution, and the connectivity. Systems analyzed in this study include different glass bead systems and natural marine sand. The media ranged in size from 0.123 to 1.0 mm, while the image volumes ranged between 7.7 and 108.9 mm³. In addition to extracting the pore network structure, the porosity, specific surface area, and representative elementary volume analysis on the porosity were calculated. Spatial correlation between pore-body sizes in the network was investigated using semivariograms and integral scale concepts. The impact of resolution on the calculated property was also investigated.

In this work, we show that microtomography is an effective tool to non-destructively extract the structure of many systems. The quality of the datasets depends on photon energy, photon flux, size of the sample, type of the sample, and size of the sample ‘features’. Results show that the developed method of extracting pore network structure is applicable to ideal and natural porous media systems. The impact of resolution on the quantification of the network structure properties varies in its significance based on feature size of the system and the properties being calculated. Therefore, a thorough resolution sensitivity analysis should be carried out to determine the degree of error associated with a system imaged at a given resolution.     

The effect of microporosity on transport properties in porous media

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.     

Numerical simulation of three-dimensional saltwater–freshwater fingering instabilities observed in a porous medium

We consider saltwater–freshwater fingering instabilities in a saturated porous medium. In the first part, we present three-dimensional results obtained from a laboratory experiment using non-invasive imaging. In the second part, we define a set of model problems in which the performed laboratory experiments can be ranged in. Due to its highly non-linear behavior and inevitable modeling errors, a detailed numerical reproduction of the physical concentration measurements cannot be expected. Nevertheless, four criteria have been identified, two quantitative and two qualitative, which facilitate a substantiated comparison of the physical experiment and the numerical simulation. With respect to these criteria a high degree of similarity could be observed. The use of these features allows a deeper understanding of the physical processes and the influence of the initial conditions.