共查询到20条相似文献，搜索用时 15 毫秒
1.
The pore throats in a porous medium control permeability, drainage, and straining through their pore scale geometry and through the way they are connected via pore bodies on the macroscale. Likewise, imbibition is controlled through the geometry of the pore bodies (pore scale) and through the way the pore bodies are connected via pore throats on the macroscale. In an effort to account for both scales at the same time we recently introduced an image-based model for pore spaces that consists of two parts related by duality: (1) a decomposition of a polyhedral pore space into polyhedral pore bodies separated by polygonal pore throats and (2) a polygonal pore network that is homotopy equivalent to the pore space. In this paper we stick to the dual concept while amending the definition of the pore throats and, as a consequence, the other elements of the dual model. Formerly, the pore throats consisted of single two-dimensional Delaunay cells, while they now usually consist of more than one two-dimensional Delaunay cell and extend all the way into the narrowing ends of the pore channel cross sections. This is the first reason for naming the amended dual model “tight”. The second reason is that the formation of the pore throats is now guided by an objective function that always attains its global optimum (tight optimization). At the end of the paper we report on simulations of drainage performed on tight dual models derived from simulated sphere packings and 3D gray-level images. The C-code for the generation of the tight dual model and the simulation of drainage is publicly available at https://jshare.johnshopkins.edu/mhilper1/public_html/tdm.html. 相似文献
2.
Saeid Jamshidi Ramin Bozorgmehry Boozarjomehry Mahmoud Reza Pishvaie 《Advances in water resources》2009
In pore network modeling, the void space of a rock sample is represented at the microscopic scale by a network of pores connected by throats. Construction of a reasonable representation of the geometry and topology of the pore space will lead to a reliable prediction of the properties of porous media. Recently, the theory of multi-cellular growth (or L-systems) has been used as a flexible tool for generation of pore network models which do not require any special information such as 2D SEM or 3D pore space images. In general, the networks generated by this method are irregular pore network models which are inherently closer to the complicated nature of the porous media rather than regular lattice networks. In this approach, the construction process is controlled only by the production rules that govern the development process of the network. In this study, genetic algorithm has been used to obtain the optimum values of the uncertain parameters of these production rules to build an appropriate irregular lattice network capable of the prediction of both static and hydraulic information of the target porous medium. 相似文献
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.
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. 相似文献
5.
The purpose of this study is to quantify the dispersivity in the longitudinal direction by upscaling pore scale mixing over a network domain and to verify the dispersivity with that obtained through the more rigorous upscaling technique, the Brownian particle tracking model (BPTM). We model a porous medium with a network of pore-units that are comprised of pore bodies and bonds of finite volume. Such a pore-unit is assumed to be a mixing cell with the steady state flow condition for a single fluid. Dispersivity can be obtained by solving the mixing cell model (MCM) for the concentration in each pore-unit and by averaging the concentrations for a large number of pore units (as a function of time and space). A minimal size of network that ascertains an asymptotic value of dispersivity was determined and verified with large size pore networks. This numerically computed dispersivity is compared with the results from the BPTM for the same porous medium and flow conditions. We show that the dispersivity obtained from the MCM is equally reliable for the heterogeneous pore-networks and can be estimated as a function of pore size heterogeneity. For homogeneous networks with the MCM, the iteration time step plays an important role. On the other hand, for networks with the BPTM, the assumption of intra-bond velocity profile affects the results. 相似文献
6.
Improved network flow models require the incorporation of increasingly accurate geometrical characterization of the microscale pore structure as well as greater information on fluid–fluid interaction (interfaces) at pore scales. We report on three dimensional (3D) pore scale medium characterization, absolute permeability computations for throat structures, and pore scale residual fluid distribution in a Berea core. X-ray computed microtomography combined with X-ray attenuating dopants is used to obtain 3D images of the pore network and to resolve phase distributions in the pore space. 相似文献
7.
A.R. Mazaheri B. Zerai G. Ahmadi J.R. Kadambi B.Z. Saylor M. Oliver G.S. Bromhal D.H. Smith 《Advances in water resources》2005,28(12):1267-1279
For single-phase flow through a network model of a porous medium, we report (1) solutions of the Navier–Stokes equation for the flow, (2) micro-particle imaging velocimetry (PIV) measurements of local flow velocity vectors in the “pores throats” and “pore bodies,” and (3) comparisons of the computed and measured velocity vectors. A “two-dimensional” network of cylindrical pores and parallelepiped connecting throats was constructed and used for the measurements. All pore bodies had the same dimensions, but three-different (square cross-section) pore-throat sizes were randomly distributed throughout the network. An unstructured computational grid for flow through an identical network was developed and used to compute the local pressure gradients and flow vectors for several different (macroscopic) flow rates. Numerical solution results were compared with the experimental data, and good agreement was found. Cross-over from Darcy flow to inertial flow was observed in the computational results, and the permeability and inertia coefficients of the network were estimated. The development of inertial flow was seen as a “two-step” process: (1) recirculation zones appeared in more and more pore bodies as the flow rate was increased, and (2) the strengths of individual recirculation zones increased with flow rate. Because each pore-throat and pore-body dimension is known, in this approach an experimental (and/or computed) local Reynolds number is known for every location in the porous medium at which the velocity has been measured (and/or computed). 相似文献
8.
岩石物理岩石孔喉结构特征对核磁T2谱影响的数值模拟 总被引：2，自引：2，他引：0
We built a three-dimensional irregular network model which can adequately describe reservoir rock pore-throat structures. We carried out numerical simulations to study the NMR T2 distribution of water-saturated rocks. The results indicate that there is a good correlation between T2 distribution and the pore radius frequency histogram. The total T2 distribution can be partitioned into pore body and pore throat parts. The effect of parameters including throat radius, pore-throat ratio, and coordination number of the micro- pore structure on the T2 distribution can be evaluated individually. The result indicates that： 1 ） with the increase of the pore throat radius, the T2 distribution moves toward longer relaxation times and its peak intensity increases; 2） with the increase of the pore-throat ratio, the T2 distribution moves towards longer T2 with the peak intensity increasing and the overlap between pore body T2 and pore throat T2 decreasing; 3） With the increase of connectivity, the short T2 component increases and peak signal intensity decreases slightly. 相似文献
9.
岩石孔喉结构特征对核磁T_2谱影响的数值模拟(英文) 总被引：1，自引：0，他引：1
《应用地球物理》2008,(2)
本文建立了能够考虑孔喉特征的储层岩石三维非规整网络模型,用数值模拟方法进一步研究饱含水岩样的核磁T_2谱特征。通过模拟储层岩石孔喉核磁共振T_2谱的特征,表明核磁共振T_2谱同岩心孔喉半径分布具有很好的相似性,对T_2谱中孔隙体、喉道信号的分解可以评价储层岩石孔隙体、喉道的半径大小及分布情况。分别模拟喉道半径、孔喉比以及孔隙连通性等孔隙结构特征对T_2谱的影响,发现:1.随着喉道半径的增大,T_2谱峰值向右移动,而且峰值处的驰豫信号强度增大;2.随着孔喉比的增大,T_2谱峰值向右移动,峰值处的驰豫信号强度增大,而且孔隙体T_2谱和喉道T_2谱逐渐分开;3.随着连通性的增大,曲线左端上升,短驰豫时间信号所占比列增大,右边的峰值有所降低。 相似文献
10.
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 mm3. 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. 相似文献
11.
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. 相似文献
12.
岩芯核磁共振（NMR）T2分布和毛管压力分析数据均在一定程度上反映了岩石的孔隙结构,理论分析表明,这两组数据具有相关性.应用NMR T2分布研究岩石孔径分布，关键是在分析两者的相关性的基础上，从T2分布构造出可靠的毛管压力曲线.但以前用饱和水T2分布构造的毛管压力曲线的方法与实际毛管压力曲线匹配性差.事实上，薄膜束缚水部分的存在引起T2分布反映的孔隙空间与毛管压力曲线反映的孔隙空间有差异.本文提出一种改进方法，在消除薄膜束缚水对T2分布的贡献后，用自由水T2分布构造毛管压力曲线.应用本文方法,对24块岩芯数据自由水T2分布构造的毛管力曲线及其孔喉半径分布与隔板毛管压力分析结果进行了对比.结果表明，改进方法对于毛管压力曲线的构造精度有明显改进，从而为NMR T2分布研究孔隙结构提供了可靠的理论和方法上的支持. 相似文献
13.
《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. 相似文献
14.
A fundamental understanding of flow in porous media at the pore-scale is necessary to be able to upscale average displacement processes from core to reservoir scale. The study of fluid flow in porous media at the pore-scale consists of two key procedures: Imaging - reconstruction of three-dimensional (3D) pore space images; and modelling such as with single and two-phase flow simulations with Lattice-Boltzmann (LB) or Pore-Network (PN) Modelling. Here we analyse pore-scale results to predict petrophysical properties such as porosity, single-phase permeability and multi-phase properties at different length scales. The fundamental issue is to understand the image resolution dependency of transport properties, in order to up-scale the flow physics from pore to core scale. In this work, we use a high resolution micro-computed tomography (micro-CT) scanner to image and reconstruct three dimensional pore-scale images of five sandstones (Bentheimer, Berea, Clashach, Doddington and Stainton) and five complex carbonates (Ketton, Estaillades, Middle Eastern sample 3, Middle Eastern sample 5 and Indiana Limestone 1) at four different voxel resolutions (4.4 µm, 6.2 µm, 8.3 µm and 10.2 µm), scanning the same physical field of view. Implementing three phase segmentation (macro-pore phase, intermediate phase and grain phase) on pore-scale images helps to understand the importance of connected macro-porosity in the fluid flow for the samples studied. We then compute the petrophysical properties for all the samples using PN and LB simulations in order to study the influence of voxel resolution on petrophysical properties. We then introduce a numerical coarsening scheme which is used to coarsen a high voxel resolution image (4.4 µm) to lower resolutions (6.2 µm, 8.3 µm and 10.2 µm) and study the impact of coarsening data on macroscopic and multi-phase properties. Numerical coarsening of high resolution data is found to be superior to using a lower resolution scan because it avoids the problem of partial volume effects and reduces the scaling effect by preserving the pore-space properties influencing the transport properties. This is evidently compared in this study by predicting several pore network properties such as number of pores and throats, average pore and throat radius and coordination number for both scan based analysis and numerical coarsened data. 相似文献
15.
碳酸盐岩储层孔隙类型多样,各种孔隙的尺寸变化范围可以跨越几个数量级,孔隙结构非常复杂,这种复杂孔隙结构和不均匀分布的多元孔隙空间使得储层电性呈现明显非阿尔奇特性.为了了解影响电阻率变化的控制因素,本次研究选取中三叠世雷口坡组的8块全直径碳酸盐岩岩样,开展了核磁共振、岩电实验、孔渗实验、压汞实验及薄片等实验,并利用数字图像分析法定量分析了孔隙结构特征.研究结果表明：①孔隙度是影响电阻率高低的重要因素,但并非唯一因素,除孔隙度以外,孔隙尺寸和数量、孔隙网络复杂程度远比吼道大小对电阻率的影响大;②在孔隙度一定的条件下,胶结指数m随储层中孤立大孔隙占比的增多而增大,当孔隙度增大到一定程度后,胶结指数m又随大孔隙占比的增多而减小,微裂缝起重要沟通作用;③在给定孔隙度时,以简单大孔隙为主的岩样表现为胶结指数m值较大,而以复杂孔隙网络、细小孔隙为主的岩样表现为胶结指数m值较小,具分散、孤立大孔隙的岩样,胶结指数m值最高;④依据孔隙几何参数与电阻率和胶结指数之间的关系,可以利用测井资料间接判别储层类型,从而提高储层有效性和含水饱和度评价精度. 相似文献
16.
低渗透砂岩油气藏已成为油气增储生产的重要勘探开发目标,但孔隙结构复杂使得储层及其有效性难以准确识别.笔者利用物性、压汞、核磁等资料,对东营凹陷南坡沙四段(Es4)低渗透砂岩孔隙结构进行分析,划分出了3种类型.核磁T2谱与毛管压力曲线都在一定程度上反映孔喉分布,但常规方法利用T2谱重构伪毛管压力曲线所得到的孔隙半径与压汞孔喉半径有较大误差,而岩石孔隙自由流体T2与压汞孔喉分布对应关系更好,以此建立了不同孔隙结构类型二者之间不同孔喉尺度对应的关系式(大尺度:线性;小尺度:分段幂函数),可在井筒剖面上通过识别孔隙结构类型,进而利用核磁共振测井(NML)定量反演孔径分布,省去了构建伪毛管曲线环节,为低渗透砂岩储层有效性评价提供了直接依据,也是测井用于定量反演储层微观孔隙结构信息的有益探索. 相似文献
17.
18.
Todd H. Skaggs 《Advances in water resources》2011,34(10):1335-1342
Critical path analysis (CPA) is a method for estimating macroscopic transport coefficients of heterogeneous materials that are highly disordered at the micro-scale. Developed originally to model conduction in semiconductors, numerous researchers have noted that CPA might also have relevance to flow and transport processes in porous media. However, the results of several numerical investigations of critical path analysis on pore network models raise questions about the applicability of CPA to porous media. Among other things, these studies found that (i) in well-connected 3D networks, CPA predictions were inaccurate and became worse when heterogeneity was increased; and (ii) CPA could not fully explain the transport properties of 2D networks. To better understand the applicability of CPA to porous media, we made numerical computations of permeability and electrical conductivity on 2D and 3D networks with differing pore-size distributions and geometries. A new CPA model for the relationship between the permeability and electrical conductivity was found to be in good agreement with numerical data, and to be a significant improvement over a classical CPA model. In sufficiently disordered 3D networks, the new CPA prediction was within ±20% of the true value, and was nearly optimal in terms of minimizing the squared prediction errors across differing network configurations. The agreement of CPA predictions with 2D network computations was similarly good, although 2D networks are in general not well-suited for evaluating CPA. Numerical transport coefficients derived for regular 3D networks of slit-shaped pores were found to be in better agreement with experimental data from rock samples than were coefficients derived for networks of cylindrical pores. 相似文献
19.
In the water flooding of mixed-wet porous media, oil may drain down to relatively low residual oil saturations (Sor). Various studies have indicated that such low saturations can only be reached when oil layers in pore corners are included in the pore-scale modelling. These processes within a macroscopic porous medium can be modelled at the pore-scale by incorporating the fundamental physics of capillary dominated displacement within idealised pore network models. Recently, the authors have developed thermodynamic criteria for oil layer existence in pores with non-uniform wettability which takes as input geometrically and topologically representative networks, to calculate realistic Sor values for mixed-wet and oil-wet sandstones [16, 21]. This previous work is developed in this paper to include (i) the visualisation of the 3D structure of this residual oil, and (ii) a statistical analysis of this “residual/remaining” oil. Both the visualisation and the statistical analysis are done under a wide range of wettability conditions, which is reported for the first time in this paper.The structure of residual oil for strongly water wet systems is well known (where residual = remaining oil) and our model agrees with this but this structure changes radically for mixed wet systems (where residual ≠ remaining) and this has not yet been visualised experimentally. We find that for more water-wet systems high final residual oil saturations are reached at relatively small amounts of water injected and this oil is present in the pores as bulk oil. On the other hand, for more oil-wet systems we find a slow decrease of the amount of remaining oil with increasing amounts of injected water. During the process, the remaining connectivity of the oil phase is increasingly provided by oil layers only, hence the slow drainage. The final residual oil saturation, only reached in the theoretical limit of an infinite amount of injected water, is almost entirely contained in large number of (relatively low volume) oil layers, which are present in pores of most radius sizes. 相似文献