共查询到20条相似文献，搜索用时 31 毫秒
1.
A quasi-static scheme based on pore space spatial statistics is presented to simulate pore-scale two-phase capillary-dominant displacement processes. The algorithm is coupled with computational fluid dynamics in order to evaluate saturation functions. Wettability heterogeneity in partial and fractional/mixed-wet media is implemented using a contact angle map. The simulation process is pixel-wised and performed directly on binary images. Bypassing and snap-off are tackled as non-wetting phase trapping mechanisms. Post-processing results include residual saturations, effective permeability and capillary pressure curves for drainage and imbibition scenarios. The primary advantages of the proposed workflow are eliminating pore space skeletisation/ discretization, superior time efficiency and minimal numerical drawbacks when compared to other direct or network-based simulation techniques. 相似文献
2.
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. 相似文献
3.
We present a semi-analytical, combinatorial approach to compute three-phase capillary entry pressures for gas invasion into pore throats with constant cross-sections of arbitrary shapes that are occupied by oil and/or water. For a specific set of three-phase capillary pressures, geometrically allowed gas/oil, oil/water and gas/water arc menisci are determined by moving two circles in opposite directions along the pore/solid boundary for each fluid pair such that the contact angle is defined at the front circular arcs. Intersections of the two circles determine the geometrically allowed arc menisci for each fluid pair. The resulting interfaces are combined systematically to allow for all geometrically possible three-phase configuration changes. The three-phase extension of the Mayer and Stowe – Princen method is adopted to calculate capillary entry pressures for all determined configuration candidates, from which the most favorable gas invasion configuration is determined. The model is validated by comparing computed three-phase capillary entry pressures and corresponding fluid configurations with analytical solutions in idealized triangular star-shaped pores. It is demonstrated that the model accounts for all scenarios that have been analyzed previously in these shapes. Finally, three-phase capillary entry pressures and associated fluid configurations are computed in throat cross-sections extracted from segmented SEM images of Bentheim sandstone. The computed gas/oil capillary entry pressures account for the expected dependence of oil/water capillary pressure in spreading and non-spreading fluid systems at the considered wetting conditions. Because these geometries are irregular and include constrictions, we introduce three-phase displacements that have not been identified previously in pore-network models that are based on idealized pore shapes. However, in the limited number of pore geometries considered in this work, we find that the favorable displacements are not generically different from those already encountered in network models previously, except that the size and shape of oil layers that are surrounded by gas and water are described more realistically. The significance of the results for describing oil connectivity in porous media accurately can only be evaluated by including throats with more complex cross-sections in three-phase pore-network models. 相似文献
4.
《Advances in water resources》2001,24(8):887-898
Fluid–fluid interfacial areas play important roles in numerous subsurface processes such as dissolution, volatilization, and adsorption. Integral expressions have been derived to estimate both entrapped (discontinuous) and free (continuous) nonwetting fluid–wetting fluid specific interfacial areas in porous media. The expressions, compatible with widely used capillary head-saturation and entrapment models, require information on capillary head-saturation relation parameters, porosity, and fluid-pair interfacial tension. In addition, information on the maximum entrapped nonwetting fluid saturation as well as the main drainage branch reversal point for water and total liquid saturations is necessary to estimate entrapped fluid interfacial areas. Implementation of the interfacial area equations in continuum-based multifluid flow simulators is straightforward since no additional parameters are needed than those required by the simulators to complete the multifluid flow computations. A limited sensitivity analysis, based on experimentally obtained parameter values, showed that imposed variations resulted in logical and consistent changes in predicted specific interfacial areas for both entrapped and free nonwetting fluid–wetting fluid systems. A direct comparison with published experimental work to test the derived expressions was limited to free air–water systems and yielded reasonable results. Such comparisons are often not possible because of the lack of information given on retention parameters, and variables used to determine nonwetting fluid entrapment. This contribution is dedicated to John W. Cary. 相似文献
5.
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. 相似文献
6.
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. 相似文献
7.
Geological heterogeneities affect the dynamics of carbon dioxide (CO2) plumes in subsurface environments in important ways. Previously, we showed how the dynamics of CO2 plumes are influenced by the multiscaled sedimentary architecture in deep brine fluvial-type reservoirs. The results confirm that representing small-scale features and the corresponding heterogeneity in saturation functions, along with hysteresis in saturation functions, are all critical to understanding capillary trapping processes. Here, we show that when heterogeneity and hysteresis are represented, the two conventional approaches for defining saturation functions, Brooks-Corey and van Genuchten, represent fundamentally different physical systems. The Brooks-Corey approach represents heterogeneity in entry pressures, and leads to trapping by capillary pinning. The van Genuchten approach represents a network of pores transporting the nonwetting fluid, across rock types, with negligible capillary entry pressure, and leads to capillary retardation. These differences significantly affect the large-scale characteristics of CO2 plumes (i.e., their mass, shape, and position). 相似文献
8.
《Advances in water resources》1998,21(7):605-615
The interface of two immiscible fluids flowing in porous media may behave in an unstable fashion. This instability is governed by the pore distribution, differential viscosity and interface tension between the two immiscible fluids. This study investigates the factors that control the interface instability at the wetting front. The development of the flow equation is based on the mass balance principle, with boundary conditions such as the velocity continuity and capillary pressure balance at the interface. By assuming that the two-phase fluids in porous media are saturated, a covariance function of the wetting front position is derived by stochastic theory. According to those results, the unstable interface between two immiscible fluids is governed by the fluid velocity and properties such as viscosity and density. The fluid properties that affect the interface instability are expressed as dimensionless parameters, mobility ratio, capillary number and Bond number. If the fluid flow is driven by gravitational force, whether the interface undergoes upward displacement or downward displacement, the variance of the unstable interface decreases with an increasing mobility ratio, increases with increasing capillary number, and decreases with increasing Bond number. For a circumstance in which fluid flow is horizontal, our results demonstrate that the capillary number does not influence the generation of the unstable interface. 相似文献
9.
孔隙尺度的喷射流流动是引起地震波速度频散和衰减的重要机制之一.目前,大多数喷射流模型仅考虑硬孔隙与微裂隙之间的局部流动,而忽略了具有不同孔隙纵横比微裂隙间的喷射流作用.为了研究各种类型孔隙间的流体流动效应,本文对经典喷射流模型进行了扩展,通过结合等效介质理论和孔隙结构模型,根据从干燥岩石超声速度-压力曲线中提取的微裂隙孔隙纵横比分布,求取出岩石中各种微裂隙的体积压缩系数,并在此基础上,利用孔隙空间的压力松弛效应对微裂隙间的喷射流效应进行了模拟,并运用Biot理论描述了硬孔隙间的宏观流动效应.扩展后的理论模型不仅考虑了微裂隙与硬孔隙间的局部流动、硬孔隙与硬孔隙间的Biot宏观流,还加入了微裂隙与微裂隙间的喷射流作用,且模型的高、低频极限始终与Mavko-Jizba理论和Gassmann方程保持一致.模型应用分析发现,对于砂岩和大部分致密灰岩样品,扩展模型均能给出与超声实验测量数据吻合良好的估计结果.此外,扩展模型预测的速度频散及衰减表明,喷射流机制在地震和测井频段发挥着重要作用,其速度频散曲线由低频至高频呈逐渐增大趋势,不具有明显的快速变化特征,与经典频散曲线形态存在显著差异;在低有效压力下,频散和衰减程度较大,喷射流机制发挥主要作用,而随着有效压力的增加,Biot宏观流机制开始占主导,频散和衰减程度逐渐减小. 相似文献
10.
岩芯核磁共振（NMR）T2分布和毛管压力分析数据均在一定程度上反映了岩石的孔隙结构,理论分析表明,这两组数据具有相关性.应用NMR T2分布研究岩石孔径分布，关键是在分析两者的相关性的基础上，从T2分布构造出可靠的毛管压力曲线.但以前用饱和水T2分布构造的毛管压力曲线的方法与实际毛管压力曲线匹配性差.事实上，薄膜束缚水部分的存在引起T2分布反映的孔隙空间与毛管压力曲线反映的孔隙空间有差异.本文提出一种改进方法，在消除薄膜束缚水对T2分布的贡献后，用自由水T2分布构造毛管压力曲线.应用本文方法,对24块岩芯数据自由水T2分布构造的毛管力曲线及其孔喉半径分布与隔板毛管压力分析结果进行了对比.结果表明，改进方法对于毛管压力曲线的构造精度有明显改进，从而为NMR T2分布研究孔隙结构提供了可靠的理论和方法上的支持. 相似文献
11.
Although there is no assumption of pore geometry in derivation of Gassmann's equation, the pore geometry is in close relation with hygroscopic water content and pore fluid communication between the micropores and the macropores. The hygroscopic water content in common reservoir rocks is small, and its effect on elastic properties is ignored in the Gassmann theory. However, the volume of hygroscopic water can be significant in shaly rocks or rocks made of fine particles; therefore, its effect on the elastic properties may be important. If the pore fluids in microspores cannot reach pressure equilibrium with the macropore system, assumption of the Gassmann theory is violated. Therefore, due to pore structure complexity, there may be a significant part of the pore fluids that do not satisfy the assumption of the Gassmann theory. We recommend that this part of pore fluids be accounted for within the solid rock frame and effective porosity be used in Gassmann's equation for fluid substitution. Integrated study of ultrasonic laboratory measurement data, petrographic data, mercury injection capillary pressure data, and nuclear magnetic resonance T2 data confirms rationality of using effective porosity for Gassmann fluid substitution. The effective porosity for Gassmann's equation should be frequency dependent. Knowing the pore geometry, if an empirical correlation between frequency and the threshold pore‐throat radius or nuclear magnetic resonance T2 could be set up, Gassmann's equation can be applicable to data measured at different frequencies. Without information of the pore geometry, the irreducible water saturation can be used to estimate the effective porosity. 相似文献
12.
含裂缝多孔介质渗透率预测是非常规油气资源勘探开发的一个紧迫问题.现有多孔介质岩石物理模型通常利用圆形孔管模拟宏观岩石孔隙空间,难以定量描述软孔隙/裂缝在压力作用下的闭合情况,缺乏裂缝/孔隙间流量交换的连通机制.本文提出含三维裂缝/软孔隙网络多孔介质模型,将储层岩石裂缝/软孔隙表示为椭圆截面微管,建立了周期性压力作用下微观裂缝流量表达式,通过网络模型和流量守恒条件,得到含有三维裂缝/软孔隙网络的多孔介质渗透率计算方法.数值算例表明,预测结果与实验数据分布范围吻合很好,能够给出不同类型岩心对应孔隙纵横比的分布图.三维裂缝/软孔隙网络模型建立了宏观可观测量与裂缝参数之间关系,能够定量分析岩石渗透率随裂缝体密度、纵横比、孔隙流体类型和围压等因素的变化规律,为复杂条件下储层渗透率预测提供了一种有效方法. 相似文献
13.
介观尺度孔隙流体流动是地震频段岩石表现出较强速度频散与衰减的主要作用.利用周期性层状孔隙介质模型,基于准静态孔弹性理论给出了模型中孔隙压力、孔隙流体相对运动速度以及固体骨架位移等物理量的数学解析表达式,同时利用Biot理论将其扩展至全频段条件下,克服了传统White模型中介质分界面处流体压力不连续的假设. 在此基础上对准静态与全频段下模型介质中孔隙压力、孔隙流体相对运动速度变化形式及其对弹性波传播特征的影响进行了讨论,为更有效理解介观尺度下流体流动耗散和频散机制提供物理依据.研究结果表明,低频条件下快纵波孔压在介质层内近于定值,慢纵波通过流体扩散改变总孔隙压力, 随频率的增加慢波所形成的流体扩散作用逐渐减弱致使介质中总孔压逐渐接近于快纵波孔压,在较高频率下孔压与应力的二次耦合作用使总孔压超过快纵波孔压.介质中孔隙流体相对运动速度与慢纵波形成的流体相对运动速度变化形式一致;随频率的增加孔隙流体逐渐从排水的弛豫状态过渡到非弛豫状态,其纵波速度-含水饱和度变化形式也从符合孔隙流体均匀分布模式过渡到斑块分布模式,同时介质在不同含水饱和度下的衰减峰值与慢纵波所形成的孔隙流体相对流动速度具有明显的相关性. 相似文献
14.
15.
A novel semi-analytical model for computation of capillary entry pressures and associated fluid configurations in arbitrary, potentially non-convex, 2D pore space geometries at uniform wettability is developed. The model computes all possible centre positions of circular arcs, and physically sound criteria are implemented to determine the set of these arcs that correspond to geometrically allowed interfaces. Interfaces and pore boundary segments are connected to form closed boundaries of identified geometrical regions. These regions are classified as either oil regions, located in the wider parts of the pore space, or as water regions located in pore space constrictions. All possible region combinations are identified and evaluated for each radius value in an iterative procedure to determine the favourable entry radius and corresponding configuration based on minimisation of free energy. The model has been validated by comparison with known analytical solutions in idealised pore geometries. In cases where different analytical solutions are geometrically possible, the model generates several oil and water regions, and the valid solution is determined by the region combination that corresponds to the most favourable entry pressure, consistent with the analytical solution. Entry pressure radii and configurations are computed in strongly non-convex pore spaces extracted from an image of Bentheimer sandstone, which demonstrates that the model captures successfully well-known characteristics of capillary behaviour at different wetting conditions. The computations also demonstrate the importance of selecting the fluid configuration of minimum change in free energy. In some cases, a merged region formed by a combination of oil and water regions corresponds to the favourable entry configuration of oil, whereas in other cases, an individual oil region may correspond to the favourable oil entry configuration. It is also demonstrated that oil entry configurations may constitute merged regions for weakly water-wet conditions and individual oil regions for strongly water-wet conditions in the same pore space. The computations show that the ratio of pore area to perimeter multiplied by the cosine of the contact angle under-predicts the entry pressure radii in Bentheimer sandstone pore spaces. An alternative formula is proposed for prediction of entry radii for nonzero contact angles based on the entry radii obtained for zero contact angles. 相似文献
16.
低渗透砂岩油气藏已成为油气增储生产的重要勘探开发目标,但孔隙结构复杂使得储层及其有效性难以准确识别.笔者利用物性、压汞、核磁等资料,对东营凹陷南坡沙四段(Es4)低渗透砂岩孔隙结构进行分析,划分出了3种类型.核磁T2谱与毛管压力曲线都在一定程度上反映孔喉分布,但常规方法利用T2谱重构伪毛管压力曲线所得到的孔隙半径与压汞孔喉半径有较大误差,而岩石孔隙自由流体T2与压汞孔喉分布对应关系更好,以此建立了不同孔隙结构类型二者之间不同孔喉尺度对应的关系式(大尺度:线性;小尺度:分段幂函数),可在井筒剖面上通过识别孔隙结构类型,进而利用核磁共振测井(NML)定量反演孔径分布,省去了构建伪毛管曲线环节,为低渗透砂岩储层有效性评价提供了直接依据,也是测井用于定量反演储层微观孔隙结构信息的有益探索. 相似文献
17.
18.
Processes in porous media governed by capillary forces, such as drainage or imbibition of wetting phase, are of great importance in different branches of soil science, petroleum engineering and hydrology. In this work we describe a new way of modeling imbibition by considering the movement and merger of interfaces at the pore level. The model is based upon a physically consistent dynamic criterion for the imbibition of a single pore originally proposed by Melrose. 相似文献
19.
以鄂尔多斯盆地X地区长6储层为研究对象,利用多尺度CT成像技术、聚焦离子束扫描电镜技术,结合Avizo软件的强大数据处理和数值模拟功能,对研究区目的层岩石样品进行不同尺度孔喉分维数重构,建立三维超低渗透砂岩储层纳米级孔隙结构模型.研究表明,微米尺度下,超低渗透砂岩储层孔喉形态呈点状、球状和条带状及管状.储集空间类型以溶蚀微孔为主且多孤立分布,局部孔隙为片状,连通性较差.纳米尺度下,超低渗透砂岩储层孔喉系统整体较发育,孔喉形态呈球状、短管状.远离裂隙处多为孤立状,连通性较差,仅具有储集能力.微裂缝、粒间缝发育部位多为短管状,有一定连通性,相当于喉道.微观非均质性较强,岩样整体较致密,局部相互连通,溶蚀孔及裂隙对储集能力、渗滤能力具有控制作用.数值计算求得目的层A、Y、C三个样品的孔隙度分别为6.95%、5.55%、4.44%,渗透率分别为0.828×10^-3μm^2、0.115×10^-3μm^2、0.00065×10^-3μm^2.聚焦离子束扫描电镜与多尺度CT成像技术相结合能够定量表征超低渗透砂岩储层微、纳米级孔隙结构. 相似文献
20.
Yongyang Sun Boris Gurevich Maxim Lebedev Stanislav Glubokovskikh Vassili Mikhaltsevitch Junxin Guo 《Geophysical Prospecting》2019,67(4):888-899
Quantifying the effects of pore-filling materials on elastic properties of porous rocks is of considerable interest in geophysical practice. For rocks saturated with fluids, the Gassmann equation is proved effective in estimating the exact change in seismic velocity or rock moduli upon the changes in properties of pore infill. For solid substance or viscoelastic materials, however, the Gassmann theory is not applicable as the rigidity of the pore fill (either elastic or viscoelastic) prevents pressure communication in the pore space, which is a key assumption of the Gassmann equation. In this paper, we explored the elastic properties of a sandstone sample saturated with fluid and solid substance under different confining pressures. This sandstone sample is saturated with octadecane, which is a hydrocarbon with a melting point of 28°C, making it convenient to use in the lab in both solid and fluid forms. Ultrasonically measured velocities of the dry rock exhibit strong pressure dependency, which is largely reduced for the filling of solid octadecane. Predictions by the Gassmann theory for the elastic moduli of the sandstone saturated with liquid octadecane are consistent with ultrasonic measurements, but underestimate the elastic moduli of the sandstone saturated with solid octadecane. Our analysis shows that the difference between the elastic moduli of the dry and solid-octadecane-saturated sandstone is controlled by the squirt flow between stiff, compliant, and the so-called intermediate pores (with an aspect ratio larger than that of compliant pore but much less than that of stiff pores). Therefore, we developed a triple porosity model to quantify the combined squirt flow effects of compliant and intermediate pores saturated with solid or viscoelastic infill. Full saturation of remaining stiff pores with solid or viscoelastic materials is then considered by the lower embedded bound theory. The proposed model gave a reasonable fit to the ultrasonic measurements of the elastic moduli of the sandstone saturated with liquid or solid octadecane. Comparison of the predictions by the new model to other solid substitution schemes implied that accounting for the combined effects of compliant and intermediate pores is necessary to explain the solid squirt effects. 相似文献