共查询到20条相似文献，搜索用时 62 毫秒
1.
While the tortuosity coefficient is commonly estimated using an expression based on total porosity, this relationship is demonstrated to not be applicable (and thus is often misapplied) over a broad range of soil textures. The fundamental basis for a correlation between the apparent diffusion tortuosity coefficient and hydraulic conductivity is demonstrated, although such a relationship is not typically considered. An empirical regression for estimating the tortuosity coefficient based on hydraulic conductivity for saturated, unconsolidated soil is derived based on results from 14 previously reported diffusion experiments performed with a broad range of soil textures. Analyses of these experimental results confirm that total porosity is a poor predictor for the tortuosity coefficient over a large range of soil textures. The apparent diffusion tortuosity coefficient is more reliably estimated based on hydraulic conductivity. 相似文献
2.
The Maqarin site in Jordan has been investigated for three decades as a natural analogue for the long term changes of materials in contact with hyper-alkaline solutions. Similar processes are expected in radioactive waste disposal sites, where cement based materials are in contact with natural rocks or other e.g. clay based materials. In this context, a numerical reactive transport model was used to study local geochemical alterations and induced porosity changes for the Maqarin marl rock in contact with the hyper-alkaline solution. The geochemical setup for the rock mineralogy and the pore water was calibrated to match measurements from the Maqarin site. The setup includes several clay and zeolite minerals, considers cation exchange processes, and a state-of-the-art model for cement phases. Similar to earlier calculations by Steefel and Lichtner (1998) who used a much simpler geochemical model, the pore clogging occurred after several hundred years at a distance of 5–10 mm from the contact to the hyper-alkaline solution. In our calculations, this was caused by a massive precipitation of ettringite and C–S–H minerals. We performed a sensitivity study by varying the intrinsic diffusion coefficient, the Archie’s law exponential factor, and the mineral surface area available for dissolution and precipitation. We found that the dissolution of clay minerals controls the availability of Al, which is needed for ettringite and C–S–H phase precipitation. Thus, the amount and kinetically controlled dissolution of clay minerals controls the spatial and temporal evolution of porosity changes. The simulations reveal that neither cation exchange processes nor the formation of zeolite minerals strongly influence the geochemical evolution of the system. 相似文献
3.
利用特征函数展开方法，给出了基于扩散效应的核磁共振Bloch控制方程的弛豫模式解表达式, 在此基础上分别给出了水饱和与油水双相饱和孔隙模型的核磁共振特性模拟结果. 结果表明，孔隙流体的核磁弛豫与孔隙大小、孔隙表面弛豫率、孔隙流体的扩散能力以及含油饱和度等有密切的关系. 对于水饱和孔隙，弛豫时间主要由孔隙大小控制. 当孔隙较大时，即使均匀大小孔隙，孔隙流体的弛豫也会表现为多指数弛豫. 而且最小模式弛豫时间与孔隙大小为非线性关系. 对于油水双相饱和孔隙，在孔隙较小时，含油饱和度对弛豫的影响主要表现在弛豫时间随含油饱和度的增加而线性减小；但在孔隙较大时，含油饱和度的增加对弛豫影响表现在两个方面，其一，孔隙水弛豫由多个弛豫模式控制逐渐转变为由最小弛豫模式控制；其二，孔隙水弛豫时间与含油饱和度表现为非线性关系. 对由实际岩芯抽象出的孔隙模型，采用本文获得的理论公式，在油水双相饱和时进行了正反演模拟. 计算结果与已有的实验结果较为一致. 相似文献
4.
The semi-empirical Kozeny–Carman (KC) equation is the most famous permeability–porosity relation, which is widely used in the field of flow in porous media and is the starting point for many other permeability models. However, this relation has many limitations from its inception, and the KC constant is an empirical parameter which was proved to be not a constant. In this paper, we briefly reviewed the KC equation, its modifications and various models for the KC constant. We then derived an analytical expression for the permeability in homogeneous porous media based on the fractal characters of porous media and capillary model. The proposed model is expressed as a function of fractal dimensions, porosity and maximum pore size. The analytical KC constant with no empirical constant is obtained from the assumption of square geometrical model. Furthermore, a distinct linear scaling law between the dimensionless permeability and porosity is found. It is also shown that our analytical permeability is more closely related to the microstructures (fractal dimensions, porosity and maximum pore size), compared to those obtained from conventional methods and models. 相似文献
5.
SAMI SOLIMAN MOHAMED 《Geophysical Prospecting》1970,18(Z1):654-665
During the last decade many hypotheses were suggested to explain the phenomenon of induced electrical polarization in ionic conductive media. The most reliable of these is Fredricksberg's. Fredricksberg (1962) supposed that the pore spaces of a rock is composed of successively narrow (active zones) and wide (inactive zones). He simulated these pore spaces by a synthetic material that has an extremely high resistance. The pore spaces were generally in tube forms which exhibited some constrictions. He saturated these tubes with an electrolyte of a given concentration. An electric current was passed through this model. He observed an induced polarization voltage after current interruption. He attributed the formation of this voltage to a concentration gradient which took place due to the presence of excess charges in the active zones. Fredricksberg introduced a parameter (9) which described the relation among the lengths and cross sectional areas of the wide zones, the number of ions within each zone after current interruption with the recorded polarizability. The aim of this work is to correlate Fredricksberg's parameter with a parameter determined for natural rocks and to show experimentally the validity of this hypothesis when applying for some varieties of sandstones and volcanic rocks. The new parameter will help to evaluate a relationship between the polarizability and the water-collecting properties of rocks. Herein, we used the tortuosity T of sandstone samples instead of the parameter φ which was used by Fredricksberg to represent the pore geometry within his model (tortuosity of the passes within the model). It was shown that both φ and T have the same relationship with the polarizability ν of the rock samples and if φ or T have very low or very high values the polarizability ν tends to its minimum value, i.e. the curve representing the relation between ν and T has a maximum point corresponding to an intermediate value of T. This result supports Fredricksberg's hypothesis and confirm his results on synthetic models. For volcanic rocks the formation factor F was used since it was difficult to determine the porosity of the samples and consequently to calculate the tortuosity T as for sandstone samples. Experimental results confirm those obtained from sandstone. The grain constituents of sandstone samples were represented on equilateral triangle and the magnitude of induced polarization ν of each sample was deduced and represented on this triangle. Equipolarizability values ν drawn on this triangle showed that TJ will increase as the silty fractions of the rock increase, where the center of this triangle (represents minimum porosity) has polarizability less than 0.25%. An attempt was made to determine the coefficient of anisotropy of volcanic rock samples using the induced polarization method. For this reason the polarizability was deduced by measuring the induced polarization voltage for two perpendicular directions in a fractured cubes of andesitic basalt samples the coefficient of anisotropy was found to be equal 1.18. 相似文献
6.
The selective radius shift model was used to relate changes in mineral volume due to precipitation/dissolution reactions to changes in hydraulic properties affecting flow in porous media. The model accounts for (i) precipitation/dissolution taking place only in the water-filled part of the pore space and further that (ii) the amount of mineral precipitation/dissolution within a pore depends on the local pore volume. The pore bundle concept was used to connect pore-scale changes to macroscopic soil hydraulic properties. Precipitation/dissolution induces changes in the pore radii of water-filled pores and, consequently, in the effective porosity. In a time step of the numerical model, mineral reactions lead to a discontinuous pore-size distribution because only the water-filled pores are affected. The pore-size distribution is converted back to a soil moisture characteristic function to which a new water retention curve is fitted under physically plausible constraints. The model equations were derived for the commonly used van Genuchten/Mualem hydraulic properties. Together with a mixed-form solution of Richards’ equation for aqueous phase flow, the model was implemented into the geochemical modelling framework PHREEQC, thereby making available PHREEQC’s comprehensive geochemical reactions. Example applications include kinetic halite dissolution and calcite precipitation as a consequence of cation exchange. These applications showed marked changes in the soil’s hydraulic properties due to mineral precipitation/dissolution and the dependency of these changes on water contents. The simulations also revealed the strong influence of the degree of saturation on the development of the saturated hydraulic conductivity through its quadratic dependency on the van Genuchten parameter α. Furthermore, it was shown that the unsaturated hydraulic conductivity at fixed reduced water content can even increase during precipitation due to changes in the pore-size distribution. 相似文献
7.
Analytical approximations of bulk and shear moduli for dry rock based on the differential effective medium theory 总被引：2，自引：0，他引：2
Differential effective medium theory has been applied to determine the elastic properties of porous media. The ordinary differential equations for bulk and shear moduli are coupled and it is more difficult to obtain accurate analytical formulae about the moduli of dry porous rock. In this paper, in order to decouple these equations we first substitute an analytical approximation for the dry‐rock modulus ratio into the differential equation and derive analytical solutions of the bulk and shear moduli for dry rock with three specific pore shapes: spherical pores, needle‐shaped pores and penny‐shaped cracks. Then, the validity of the analytical approximations is tested by integrating the full differential effective medium equation numerically. The analytical formulae give good estimates of the numerical results over the whole porosity range for the cases of the three given pore shapes. These analytical formulae can be further simplified under the assumption of small porosity. The simplified formulae for spherical pores are the same as Mackenzie's equations. The analytical formulae are relatively easy to analyse the relationship between the elastic moduli and porosity or pore shapes and can be used to invert some rock parameters such as porosity or pore aspect ratio. The predictions of the analytical formulae for experimental data show that the formulae for penny‐shaped cracks are suitable to estimate the elastic properties of micro‐crack rock such as granite, they can be used to estimate the crack aspect ratio while the crack porosity is known and also to estimate the crack porosity evolution with pressure if the crack aspect ratio is given. 相似文献
8.
本文运用随机游走方法模拟了砂岩储层中流体的核磁共振(NMR)响应及其受限扩散现象.通过改变数字岩心的分辨率模拟生成不同孔隙尺寸的砂岩,研究了不同孔隙尺寸砂岩饱含水时流体扩散系数随扩散时间的变化关系,同时模拟了砂岩饱和单相流体和两相流体的NMR响应;研究了流体的受限扩散系数与横向弛豫时间T2的关系,分析了表面弛豫率和胶结指数对润湿相流体受限扩散系数线位置的影响,并将其用于解释砂岩储层的D-T2分布.结果表明:孔隙流体的扩散系数会随扩散时间的增加而逐渐减小并趋于定值.随着岩石孔隙尺寸的减小,受限扩散现象越明显,受限扩散对岩石NMR响应的影响也越大.润湿相流体受限扩散系数线的位置受岩石胶结指数和表面弛豫率的影响较大.由于润湿相流体扩散系数减小,导致D-T2分布中润湿相流体信号偏离其自由扩散系数线,需要利用流体的受限扩散系数线准确识别D-T2分布中的润湿相流体. 相似文献
9.
In this paper,the mechanical response of saturated geological rock under tidal force is explored by poroelastic theory.First,we use the free energy formula of saturated rock under a tidal force to study the relationships of pore pressure with stress,and stress with strain.Then we analyze the relationship between rock strain and tidal potential by the equilibrium differential equations of saturated rock under tidal force.Finally,we derive the physical relationship between the two parameters (pore pressure and tidal mean stress) of saturated rock and tidal potential.The relationship shows that:pore pressure is directly proportional with tidal potential,but tidal mean stress of saturated rock is inversely proportional with tidal potential.The ratio coefficient is related not only to the Lame coefficients of rock skeletons,but also to the Biot modulus.By using this model to analyze observational well water level of C-18 well which locates in Huili,Sichuan Province,the well level response coefficient (D) was estimated.This way,we derive the Skempton coefficient (B),the coefficient A and C which refer to the response coefficients of pore pressure and tidal stress to tidal potential respectively.Then we compare the differences among each coefficient in coupling and uncoupling conditions.It shows that for saturated rocks,the response of stress and pore pressure to earth tides is a product of coupling,and it is necessary to take into account the coupling effect when we study the mechanical response.The model will provide the basis not only for the study of mechanics and hydrodynamics of well-confined aquifer systems,and the mechanics of faulting under tidal force,but also for quantitative research of the triggering mechanism of tidal forces. 相似文献
10.
本文应用孔隙弹性理论,探讨了引潮力作用下饱和地质岩体的力学响应.首先通过引潮力作用下饱和岩体的自由能表达式,得到岩体孔压与应力、应力与应变之间的关系;然后从引潮力作用下饱和岩体的平衡微分方程出发,结合流-固耦合理论,分析了饱水岩体应变与引潮位之间的关系;最后推导出饱和岩体的两大力学物理量--孔压和潮汐应力(平均应力)与引潮位之间的物理关系.模型表明:饱和岩体孔压与引潮位成反比,平均潮汐应力与引潮位成正比;比例系数不仅与岩体骨架的Lame系数有关,而且与Biot模量有关.将模型应用于会理川-18井水位变化分析,估计出水位响应系数D,并以此为基础,求得岩体孔压、潮汐应力与引潮位的相关系数(A和C)及Skempton系数B.最后对比分析了耦合条件下与不考虑耦合时得到的各参数之间的差异,分析表明:对饱和地质岩体而言,应力、孔压对引潮力的响应是流-固耦合作用的产物;研究其力学响应时必须充分考虑耦合效应.模型的建立,为研究引潮力作用下井-承压含水层系统力学、水动力学、与地震有关的断层力学以及引潮力触发机制的定量研究提供了基础. 相似文献
11.
《Journal of Applied Geophysics》2006,58(3):232-252
Several parameters of pore geometry are needed for estimating permeability which is a key parameter for the characterization of reservoir sandstones. Powerful techniques for probing the pore space are the self-diffusion and the relaxation time NMR methods. However, the quality of results depends on the petrophysical model which underlies the interpretation of measurements. We applied the pulsed field gradient nuclear magnetic resonance (PFG-NMR) technique and measured time-dependent self-diffusion coefficients, D(Δ), of water in anhydrite cemented sandstones with low porosity and high tortuosity. The conventional method of fitting data with a function (Padé approximation) of the surface-to-volume ratio and the tortuosity yielded uncertain results. As part of a novel approach, we developed a numerical simulation code based on physical principles and a fractal pore space model. We compared our method with the Padé approximation and tested with data from the literature. For porous media with low tortuosity values and simple geometry such as randomly packed glass beads, both methods are in good agreement and give similar results. In sedimentary rocks, however, the new method of calculation is able to determine more accurate details of the pore geometry. However, its main advantage occurs in porous media with fractal geometry or with constricted pores, where the D(Δ)/D0 curves deviate from the characteristic shape that is found for randomly packed glass beads. 相似文献
12.
A pore‐scale model based on measured particle size distributions has been used to quantify the changes in pore space geometry of packed soil columns resulting from a dilution in electrolyte concentration from 500 to 1 mmol l?1 NaCl during leaching. This was applied to examine the effects of particle release and re‐deposition on pore structure and hydraulic properties. Two different soils, an agricultural soil and a mining residue, were investigated with respect to the change in hydraulic properties. The mining residue was much more affected by this process with the water saturated hydraulic conductivity decreasing to 0·4% of the initial value and the air‐entry value changing from 20 to 50 cm. For agricultural soil, there was little detectable shift in the water retention curve but the saturated hydraulic conductivity decreased to 8·5% of the initial value. This was attributed to localized pore clogging (similar to a surface seal) affecting hydraulic conductivity, but not the microscopically measured pore‐size distribution or water retention. We modelled the soil structure at the pore scale to explain the different responses of the two soils to the experimental conditions. The size of the pores was determined as a function of deposited clay particles. The modal pore size of the agricultural soil as indicated by the constant water retention curve was 45 µm and was not affected by the leaching process. In the case of the mining residue, the mode changed from 75 to 45 µm. This reduction of pore size corresponds to an increase of capillary forces that is related to the measured shift of the water retention curve. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
13.
Differential effective medium (DEM) theory is applied to determine the elastic properties of dry rock with spheroidal pores. These pores are assumed to be randomly oriented. The ordinary differential equations for bulk and shear moduli are coupled and it is more difficult to obtain accurate analytical formulae about the moduli of dry porous rock. In this paper, we derive analytical solutions of the bulk and shear moduli for dry rock from the differential equations by applying an analytical approximation for dry-rock modulus ratio, in order to decouple these equations. Then, the validity of this analytical approximation is tested by integrating the full DEM equation numerically. The analytical formulae give good estimates of the numerical results over the whole porosity range. These analytical formulae can be further simplified under the assumption of small porosity. The simplified formulae for spherical pores (i.e., the pore aspect ratio is equal to 1) are the same as Mackenzie's equations. The analytical formulae are relatively easy to analyze the relationship between the elastic moduli and porosity or pore shapes, and can be used to invert some rock parameters such as porosity or pore aspect ratio. The predictions of the analytical formula for the sandstone experimental data show that the analytical formulae can accurately predict the variations of elastic moduli with porosity for dry sandstones. The effective elastic moduli of these sandstones can be reasonably well characterized by spheroidal pores, whose pore aspect ratio was determined by minimizing the error between theoretical predictions and experimental measurements. For sandstones the pore aspect ratio increases as porosity increases if the porosity is less than 0.15, whereas the pore aspect ratio remains relatively stable (about 0.14) if the porosity is more than 0.15. 相似文献
14.
15.
We used the 3D continuum-scale reactive transport models to simulate eight core flood experiments for two different carbonate rocks. In these experiments the core samples were reacted with brines equilibrated with pCO2 = 3, 2, 1, 0.5 MPa (Smith et al., 2013 [27]). The carbonate rocks were from specific Marly dolostone and Vuggy limestone flow units at the IEAGHG Weyburn-Midale CO2 Monitoring and Storage Project in south-eastern Saskatchewan, Canada. Initial model porosity, permeability, mineral, and surface area distributions were constructed from micro tomography and microscopy characterization data. We constrained model reaction kinetics and porosity–permeability equations with the experimental data. The experimental data included time-dependent solution chemistry and differential pressure measured across the core, and the initial and final pore space and mineral distribution. Calibration of the model with the experimental data allowed investigation of effects of carbonate reactivity, flow velocity, effective permeability, and time on the development and consequences of stable and unstable dissolution fronts.The continuum scale model captured the evolution of distinct dissolution fronts that developed as a consequence of carbonate mineral dissolution and pore scale transport properties. The results show that initial heterogeneity and porosity contrast control the development of the dissolution fronts in these highly reactive systems. This finding is consistent with linear stability analysis and the known positive feedback between mineral dissolution and fluid flow in carbonate formations. Differences in the carbonate kinetic drivers resulting from the range of pCO2 used in the experiments and the different proportions of more reactive calcite and less reactive dolomite contributed to the development of new pore space, but not to the type of dissolution fronts observed for the two different rock types. The development of the dissolution front was much more dependent on the physical heterogeneity of the carbonate rock. The observed stable dissolution fronts with small but visible dissolution fingers were a consequence of the clustering of a small percentage of larger pores in an otherwise homogeneous Marly dolostone. The observed wormholes in the heterogeneous Vuggy limestone initiated and developed in areas of greater porosity and permeability contrast, following pre-existing preferential flow paths.Model calibration of core flood experiments is one way to specifically constrain parameter input used for specific sites for larger scale simulations. Calibration of the governing rate equations and constants for Vuggy limestones showed that dissolution rate constants reasonably agree with published values. However the calcite dissolution rate constants fitted to the Marly dolostone experiments are much lower than those suggested by literature. The differences in fitted calcite rate constants between the two rock types reflect uncertainty associated with measured reactive surface area and appropriately scaling heterogeneous distribution of less abundant reactive minerals. Calibration of the power-law based porosity–permeability equations was sensitive to the overall heterogeneity of the cores. Stable dissolution fronts of the more homogeneous Marly dolostone could be fit with the exponent n = 3 consistent with the traditional Kozeny–Carman equation developed for porous sandstones. More impermeable and heterogeneous cores required larger n values (n = 6–8). 相似文献
16.
经典的微分等效介质(DEM)理论可用于确定多孔介质的弹性性质,但由于缺乏多重孔DEM方程,其估计的多重孔岩石的等效弹性模量依赖于包裹体(即不同孔隙纵横比的孔或缝)的添加顺序.本文首先从Kuster-Toksöz理论出发建立了Zimmermann和Norris两种形式的多重孔DEM方程.Norris形式的多重孔DEM方程预测的等效弹性模量总是位于Hashin-Shtrikman上下限内,而Zimmermann形式的多重孔DEM方程有时会越界.然后,通过使用干燥岩石模量比的解析近似式,对两个相互耦合的Norris形式DEM方程进行解耦得到干燥多重孔岩石的体积和剪切模量解析式.用全DEM方程的数值解对解析近似式的有效性进行了测试,解析公式的计算结果在整个孔隙度分布区间与数值解吻合良好.对实验室测量数据在假设岩石含有双重孔隙的情形下用双重孔DEM解析公式对岩石的弹性模量进行了预测,结果表明,解析式准确地预测了弹性模量随孔隙度的变化.双重孔(即软、硬孔)DEM解析模型可用来反演各孔隙类型的孔隙体积比,它可以通过实验室测量与理论预测之间的平方误差最小反演得到.砂岩样品的反演结果揭示,软孔的孔隙体积百分比与粘土含量没有明显的相关性. 相似文献
17.
通过融合Gassmann方程和由微分等效介质理论建立的干岩石骨架模型--DEM解析模型,本文提出根据纵波(和横波)速度反演岩石等效孔隙纵横比进行储层孔隙结构评价和横波速度预测的方法.首先,利用Gassmann方程和DEM解析模型建立岩石的纵、横波速度与密度、孔隙度、饱和度和矿物组分等各参数之间的关系;其次,将岩石孔隙等效为具有单一纵横比的理想椭球孔,应用非线性全局寻优算法来寻找最佳的等效孔隙纵横比使得理论预测与实际测量的弹性模量之间的误差最小;最后,将反演得到的等效孔隙纵横比代入到Gassmann方程和DEM解析模型中构建横波速度.实验室和井孔测量数据应用表明,反演得到的等效孔隙纵横比可准确反映储层的孔隙结构,对于裂缝型储层如花岗岩,其孔隙纵横比通常小于0.025,而对于孔隙型储层如砂岩,其孔隙纵横比通常大于0.08.只利用纵波与同时利用纵、横波反演得到的孔隙纵横比结果几乎完全一致,而且由纵波构建的横波与实测横波吻合良好,说明本文提出的等效孔隙纵横比反演及其横波速度预测方法是有效的. 相似文献
18.
Luong Duy Thanh Phan Van Do Nguyen Van Nghia Nguyen Xuan Ca 《Geophysical Prospecting》2018,66(4):753-766
Streaming potential is the result of coupling between a fluid flow and an electric current in porous rocks. The modified Helmholtz–Smoluchowski equation derived for capillary tubes is mostly used to determine the streaming potential coefficient of porous media. However, to the best of our knowledge, the fractal geometry theory is not yet applied to analyse the streaming potential in porous media. In this article, a fractal model for the streaming potential coefficient in porous media is developed based on the fractal theory of porous media and on the streaming potential in a capillary. The proposed model is expressed in terms of the zeta potential at the solid?liquid interface, the minimum and maximum pore/capillary radii, the fractal dimension, and the porosity of porous media. The model is also examined by using another capillary size distribution available in published articles. The results obtained from the model using two different capillary size distributions are in good agreement with each other. The model predictions are then compared with experimental data in the literature and those based on the modified Helmholtz–Smoluchowski equation. It is shown that the predictions from the proposed fractal model are in good agreement with experimental data. In addition, the proposed model is able to reproduce the same result as the Helmholtz–Smoluchowski equation, particularly for high fluid conductivity or large grain diameters. Other factors influencing the streaming potential coefficient in porous media are also analysed. 相似文献
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.
Prediction of diffusion coefficients in porous media using tortuosity factors based on interfacial areas 总被引：1，自引：0，他引：1
Determination of aqueous phase diffusion coefficients of solutes through porous media is essential for understanding and modeling contaminant transport. Prediction of diffusion coefficients in both saturated and unsaturated zones requires knowledge of tortuosity and constrictivity factors. No methods are available for the direct measurement of these factors, which are empirical in their definition. In this paper, a new definition for the tortuosity factor is proposed, as the real to ideal interfacial area ratio. We define the tortuosity factor for saturated porous media (tau5) as the ratio S/S(o) (specific surface of real porous medium to that of an idealized capillary bundle). For unsaturated media, tortuosity factor (tau(a)) is defined as a(aw)/a(aw),o (ratio of the specific air-water interfacial area of real and the corresponding idealized porous medium). This tortuosity factor is suitably measured using sorptive tracers (e.g., nitrogen adsorption method) for saturated media and interfacial tracers for unsaturated media. A model based on this new definition of tortuosity factors, termed the interfacial area ratio (IAR) model, is presented for the prediction of diffusion coefficients as a function of the degree of water saturation. Diffusion coefficients and diffusive resistances measured in a number of saturated and unsaturated granular porous media, for solutes in dilute aqueous solutions, agree well with the predictions of the IAR model. A comparison of permeability of saturated sands estimated based on tau(s) and the same based on the Kozeny-Carman equation confirm the usefulness of the tau(s) parameter as a measure of tortuosity. 相似文献