共查询到20条相似文献，搜索用时 171 毫秒
1.
One‐dimensional analytical solution for mesoscopic flow induced damping in a double‐porosity dual‐permeability material 
下载免费PDF全文
下载免费PDF全文 Heterogeneities, such as fractures and cracks, are ubiquitous in porous rocks. Mesoscopic heterogeneities, that is, heterogeneities on length scales much larger than typical pore size but much smaller than the wavelength, are increasingly believed to be responsible for significant wave energy loss in the seismic frequency band. When a compressional wave stresses a material containing mesoscopic heterogeneities, the more compliant parts of the material (e.g., fractures and cracks) respond with a greater fluid pressure than the stiffer portions (e.g., matrix pores). The induced fluid flow, resulting from the pressure gradients developed on such scale, is called mesoscopic flow. In the present study, the double‐porosity dual‐permeability model is adopted to incorporate mesoscopic heterogeneities into rock models to account for the attenuation of wave energy. Based on the model, the damping effect due to mesoscopic flow in a one‐dimensional porous structure is investigated. Analytical solutions for several boundary‐value problems are obtained in the frequency domain. The dynamic responses of infinite and finite porous layer are examined. Numerical calculations show that the damping effect of mesoscopic flow is significant on the pore pressure response and the resulting effective stress. For the displacement, the effect is seen only at the very low frequency range or near the resonance frequencies. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
2.
Statistical characterization and stochastic modeling of pore networks in relation to fluid flow 总被引：18，自引：0，他引：18
R. D. Hazlett 《Mathematical Geology》1997,29(6):801-822
3.
《Physics and Chemistry of the Earth, Part A: Solid Earth and Geodesy》1999,24(7):593-599
We present spatial distributions for pore path length and coordination number, pore throat size and nodal pore volume obtained for a 1.53 mm3 volume of 12.1% porosity Fontainebleau sandstone. The sandstone was imaged using Synchrotron X-Ray computed microtomography at 6 micron resolution. The spatial distributions were computed based upon three dimensional medial axis analysis of the void space in the image. We also present vesicle size and vesicle-vesicle contact surface area distributions for a 1.36 mm length of a 6.36 mm diameter core of basalt from a vesiculated lava flow imaged at 20 micron resolution. 相似文献
4.
华南成矿省福建魁歧晶洞花岗岩样品孔隙
结构的工业X CT三维可视化研究 总被引：1，自引：0，他引：1
地球科学的发展与各种分析测试技术密不可分。尤其是,随着现代科技的长足进步,产生了更为先进的现代分析测试手段使其在地球科学中的作用日益为世人瞩目,工业X-CT(X-ray Computed Tomography)测量技术即是其中之一。它具有形象、直观显示地质岩矿样品内部的、无损的、原位的二维或三维(立体)信息等特点,可实现对地质岩矿样品CT二维断层图像进行全数字化的三维可视化处理与分析。本研究针对华南成矿省福建魁歧晶洞花岗岩KQ-6样品的孔隙结构,依据工业X-CT技术的原理及方法,获得了样品的CT二维序列图像,借助于数字图像处理分析工具,得到了样品的孔隙度和孔径分布统计结果,孔隙度测量值为8.14%,孔径尺寸0.4~2mm;基于此,开展了该地质样品的三维可视化孔隙结构研究,获得了样品的孔隙连通的拓扑图形,直观展现了样品内部孔隙连通性,这对精细刻画地质样品内部结构特征及成矿流体的研究具有重要意义。 相似文献
5.
弹塑性条件下岩土孔隙介质有效应力系数理论模型 总被引：3，自引：1，他引：2
分析了有效应力系数的物理机理及其主要影响因素,在此基础上,提出了等效孔隙连通率的概念,以此表征岩土类孔隙介质的结构和孔隙之间的连通性,建立了有效应力系数张量演化的普适性理论模型。基于大理岩峰后和砂岩不同塑性变形阶段的有效应力测试试验,分析了塑性条件下影响有效应力系数的主要因素,结果表明塑性条件下,影响有效应力系数的主要因素为等效孔隙连通率。利用试验数据,通过拟合得到了大理岩和砂岩的等效孔隙连通率随应变的演化规律,从而得到了有效应力系数与变形的关系。研究成果为弹塑性条件下的流固耦合研究提供了基础性的理论和方法支持。 相似文献
6.
Theoretical analyses of chemical dissolution‐front instability in fluid‐saturated porous media under non‐isothermal conditions 下载免费PDF全文
下载免费PDF全文 This paper mainly deals with the theoretical aspects of chemical dissolution‐front instability problems in two‐dimensional fluid‐saturated porous media under non‐isothermal conditions. In the case of the mineral dissolution ratio (that is defined as the ratio of the dissolved‐mineral equilibrium concentration in the pore fluid to the molar concentration of the dissolvable mineral in the solid matrix of the fluid‐saturated porous medium) approaching zero, the corresponding critical condition has been mathematically derived when temperature variation effects are considered. As a complementary tool, the computational simulation method is used to simulate the morphological evolution of chemical dissolution fronts in two‐dimensional fluid‐saturated porous media under non‐isothermal conditions. The related theoretical and numerical results have demonstrated that: (i) a temperature increase in a non‐isothermal chemical dissolution system can have some influence on the propagation speed of the planar chemical dissolution front in the system. Generally, the chemical dissolution front in the non‐isothermal chemical dissolution system propagates slower than that in the counterpart isothermal chemical dissolution system when the temperature of the non‐isothermal chemical dissolution system is higher than that of the counterpart isothermal chemical dissolution system; (ii) a temperature increase in the non‐isothermal chemical dissolution system can stabilize the chemical dissolution front propagating in the system, because it can cause a decrease in the Zhao number of the system but does not affect the critical Zhao number of the system; and (iii) the temperature gradient in the upstream direction of a chemical dissolution front is smaller than that in the downstream direction of the chemical dissolution front when the non‐isothermal chemical dissolution system is supercritical. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
7.
LI Peng JIA Chengzao JIN Zhijun LIU Quanyou BI He ZHENG Min WU Songtao HUANG Zhenkai 《《地质学报》英文版》2020,94(2):219-232
Pore distribution and micro pore-throat structure characteristics are significant for tight oil reservoir evaluation, but their relationship remains unclear. This paper selects the tight sandstone reservoir of the Chang 7 member of the Xin’anbian Block in the Ordos Basin as the research object and analyzes the pore size distribution and micro pore-throat structure using field emission scanning electron microscopy(FE-SEM), high-pressure mercury injection(HPMI), highpressure mercury injection, and nuclear magnetic resonance(NMR) analyses. The study finds that:(1) Based on the pore size distribution, the tight sandstone reservoir is characterized by three main patterns with different peak amplitudes. The former peak corresponds to the nanopore scale, and the latter peak corresponds to the micropore scale. Then, the tight sandstone reservoir is categorized into three types: type 1 reservoir contains more nanopores with a nanopore-to-micropore volume ratio of 82:18;type 2 reservoir has a nanopore-to-micropore volume ratio of 47:53;and type 3 reservoir contains more micropores with a nanopore-to-micropore volume ratio of 35:65.(2) Affected by the pore size distribution, the throat radius distributions of different reservoir types are notably offset. The type 1 reservoir throat radius distribution curve is weakly unimodal, with a relatively dispersed distribution and peak ranging from 0.01 μm to 0.025 μm. The type 2 reservoir’s throat radius distribution curve is single-peaked with a wide distribution range and peak from 0.1 μm to 0.25 μm. The type 3 reservoir’s throat radius distribution curve is single-peaked with a relatively narrow distribution and peak from 0.1 μm to 0.25 μm. With increasing micropore volume, pore-throat structure characteristics gradually improve.(3) The correlation between micropore permeability and porosity exceeds that of nanopores, indicating that the development of micropores notably influences the seepage capacity. In the type 1 reservoir, only the mean radius and effective porosity have suitable correlations with the nanopore and micropore porosities. The pore-throat structure parameters of the type 2 and 3 reservoirs have reasonable correlations with the nanopore and micropore porosities, indicating that the development of these types of reservoirs is affected by the pore size distribution. This study is of great significance for evaluating lacustrine tight sandstone reservoirs in China. The research results can provide guidance for evaluating tight sandstone reservoirs in other regions based on pore size distribution. 相似文献
8.
An analytical solution is proposed for transient flow and deformation coupling of a fluid‐saturated poroelastic medium within a finite two‐dimensional (2‐D) rectangular domain. In this study, the porous medium is assumed to be isotropic, homogeneous, and compressible. In addition, the point sink can be located at an arbitrary position in the porous medium. The fluid–solid interaction in porous media is governed by the general Biot's consolidation theory. The method of integral transforms is applied in the analytical formulation of closed‐form solutions. The proposed analytical solution is then verified against both exact and numerical results. The analytical solution is first simplified and validated by comparison with an existing exact solution for the uncoupled problem. Then, a case study for pumping from a confined aquifer is performed. The consistency between the numerical solution and the analytical solution confirms the accuracy and reliability of the analytical solution presented in this paper. The proposed analytical solution can help us to obtain in‐depth insights into time‐dependent mechanical behavior due to fluid withdrawal within finite 2‐D porous media. Moreover, it can also be of great significance to calibrate numerical solutions in plane strain poroelasticity and to formulate relevant industry norms and standards. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
9.
渗透率是评价多孔介质导水性能的主要参数,是导水介质微观孔隙结构在宏观上的表现特征。通过扫描电镜实验和图像处理技术,统计出砂岩微观孔隙几何参数,并结合分形理论探究了砂岩微观孔隙结构对渗透率的影响。同时,以实测渗透率为约束条件,确定了对渗流起作用的最小孔隙半径rmin区间,该区间可作为渗透率预测时最小孔隙半径rmin的最优取值区间,并获取了与之对应的孔隙半径累积百分数区间。结果表明:砂岩最小孔隙半径rmin的取值对渗透率影响较大,随着砂岩最小孔隙半径rmin取值的增大,其渗透率先急剧下降,后渐渐地趋于稳定,并伴有上升趋势。在致密砂岩中,理论渗透率与实测渗透率相等时,对应的最小孔隙半径rmin落在孔隙半径累积百分数为15%~25%内,该孔隙半径累积百分数区间对应的孔隙半径区间可作为计算砂岩理论渗透率时最小孔隙半径rmin的最优取值区间,并且随着砂岩渗透率的增大,最小孔隙半径rmin对应下的孔隙半径累积百分数逐渐增加。 相似文献
10.
The aim of this paper is to present a three‐dimensional (3D) finite element modeling of heat and mass transfer phenomena in partially saturated open porous media with random fields of material properties. Randomness leads to transfer processes within the porous medium that naturally need a full 3D modeling for any quantitative assessment of these processes. Nevertheless, the counterpart of 3D modeling is a significant increase in computations cost. Therefore, a staggered solution strategy is adopted which permits to solve the equations sequentially. This appropriate partitioning reduces the size of the discretized problem to be solved at each time step. It is based on a specific iterative algorithm to account for the interaction between all the transfer processes. Accordingly, a suitable linearization of mass convective boundary conditions, consistent with the staggered algorithm, is also derived. After some validation tests, the 3D numerical model is used for studying the drying process of a cementitious material with regard to its intrinsic permeability randomness. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
11.
孔隙是多孔介质内渗流的发生场所,与介质渗透性存在必然的联系。珊瑚砂因其特殊的物质来源和形成过程,较陆源砂具有截然不同的孔隙特性。通过一系列微观研究手段,从本质上揭示了珊瑚砂特殊孔隙性质的原因。研究发现,从孔隙形状、孔喉尺寸和整体连通性3个角度描述孔隙性质较为合理。其中,孔隙形状用形状因子度量,孔喉尺寸包括孔隙半径和喉道半径,珊瑚砂多孔介质整体连通性利用配位数进行描述。而影响孔隙形状、孔喉尺寸和整体连通性的主导因素包括颗粒形状和颗粒表面粗糙度两方面。其中颗粒形状主要影响孔隙形状、喉道尺寸、孔喉尺寸离散性和介质内部连通性的均匀分布情况。颗粒表面粗糙度主要影响孔隙形状、孔隙形状离散性、孔隙尺寸和介质整体连通性。 相似文献
12.
This paper deals with the theoretical aspects of nonaqueous phase liquid (NAPL)‐dissolution‐induced instability in two‐dimensional fluid‐saturated porous media including solute dispersion effects.After some weaknesses associated with the previous work are analyzed and overcome, a comprehensive dimensionless number, known as the Zhao number, is proposed to represent the main driving force and three controlling mechanisms of an NAPL‐dissolution system that has a finite domain. The linear stability analysis is carried out to derive the critical value of the comprehensive dimensionless number of the NAPL‐dissolution system in a limit case as the ratio of the equilibrium concentration to the density of the NAPL approaches zero. As a result, a theoretical criterion that can be used to assess the instability of planar NAPL‐dissolution fronts in two‐dimensional fluid‐saturated porous media of finite domains has been established. Not only can the present theoretical results be used for the theoretical understanding of the effect of solute dispersion on the instability of an NAPL‐dissolution front in the fluid‐saturated porous medium of either a finite domain or an infinite domain, but also they can be used as benchmark solutions for verifying numerical methods employed to simulate detailed morphological evolution processes of NAPL‐dissolution fronts in two‐dimensional fluid‐saturated porous media. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
13.
Determination of transport properties of geomaterials is an important issue in many fields of engineering analysis and design. For example, in petroleum engineering, in situ permeability of an oil reservoir may be crucial in establishing its viability for exploitation, whilst prevention of leakage from underground storage facilities for oil and gas, nuclear waste as well as viability of CO2 sequestration projects crucially depends on its long‐term values. Permeability is indirectly related to the porosity, pore‐size distribution and pore architecture of the porous media. These parameters evolve when a strain field is imposed. Physical measurement of permeability under a strain field in laboratory conditions is difficult, expensive and prone to a number of uncertainties. In the past, pore network models have been used to compute permeability of materials under stress/strain‐free conditions. In this paper, we propose an enhanced pore network model to compute permeability of rocks and rock‐like porous media under a stress/strain field. Data of pore‐size distribution obtained from mercury intrusion porosimetry are used to compute permeability of rock samples from various unspecified oilfields in the world. It is shown that the two permeabilities can be predicted from the model with sufficient accuracy. A hypothesis for change in porosity, pore‐size distribution and pore architecture as a result of imposed mechanical strains is then proposed. Based on this, permeability is computed again for one of the rock samples under uniaxial and triaxial compressive and tensile strain fields. It is shown that depending on the state of strain field imposed, permeability evolves in an anisotropic manner. Permeability under tensile strain field increases dramatically compared with the reduction that takes place under compressive strain field of the same magnitude. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
14.
15.
Visualizing complex pore structure and fluid flow in porous media using 3D printing technology and LBM simulation 下载免费PDF全文
下载免费PDF全文 Pore structures of porous media and properties of fluid flow are key factors for the study of non-Darcy groundwater flow. However, it is difficult to directly observe pore structures and flow properties, resulting in a “black box” problem of porous media. This problem has hindered the in-depth study of the groundwater flow mechanism at the pore scale. In recent years, 3D rapid prototyping technology has seen tremendous development. 3D printing provides digital models and printing models of porous media with clear internal structure. Thus, Lattice Boltzmann Method can be used to simulate the flow processes at the pore scale based on real pore structures. In this study, 3D printing cores and Lattice Boltzmann Method were coupled to conduct both laboratory and numerical experiments in spherical porous media with different sphere diameters and periodic arrays. The LBM simulation results show a good agreement with laboratory experimental results. With the advantages of LBM and 3D printing, this approach provides a visualization of the complex pore structure and fluid flow in pores, which is a promising method for studies of non-Darcy groundwater flow at the pore scale. 相似文献
16.
为揭示低渗孔隙型碳酸盐岩油藏储层微观孔隙结构特征和水驱剩余油分布特征,以中东某孔隙型碳酸盐岩油藏为例,基于油藏物性测试、恒速压汞孔喉参数统计、CT扫描孔隙参数分析和水驱油流体微观分布特征等,得出孔隙型碳酸盐岩储层的微观孔隙、喉道发育特征及其影响下的剩余油分布规律;在此基础上,将孔隙型碳酸盐岩与砂岩储层研究结果进行对比,进一步明确了两类储层孔隙结构及剩余油分布特征差异。结果表明：与低渗砂岩储层相比,低渗碳酸盐岩储层呈现出"孔隙大而数量少""喉道数量多且类型全"的特征,并建立孔喉发育特征模型;由于孔隙型微观孔隙结构的影响,低渗碳酸盐岩储层大孔隙中的剩余油仍然是需要研究的对象,同时明确了剩余油挖潜开发方式建议。 相似文献
17.
A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
18.
A multiphase coupled elasto‐viscoplastic finite element analysis formulation, based on the theory of porous media, is used to describe the rainfall infiltration process into a one‐dimensional soil column. Using this framework, we have numerically analyzed the generation of pore water pressure and deformations when rainfall is applied to the soil. A parametric study, including rainfall intensity, soil–water characteristic curves, and permeability, is carried out to observe their influence on the changes in pore water pressure and volumetric strain. From the numerical results, it is shown that the generation of pore water pressure and volumetric strain is mainly controlled by material parameters α and n′ that describe the soil–water characteristic curve. A comparison with the laboratory results shows that the proposed method can describe very well the characteristics observed during the experiments of one‐dimensional water infiltration into a layered unsaturated soil column. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
19.
A modified three‐dimensional discontinuous deformation analysis (3D‐DDA) method is derived using four‐noded tetrahedral elements to improve the accuracy of current 3D‐DDA algorithm in practical applications. The analysis program for the modified 3D‐DDA method is developed in a C++ environment and its accuracy is illustrated through comparisons with several analytical solutions that are available for selected problems. The predicted solutions for these problems using the modified 3D‐DDA approach all show satisfactory agreement with the corresponding analytical results. Results presented in this paper demonstrate that the modified 3D‐DDA method with discontinuous modeling capabilities offers a useful computational tool to determine stresses and deformations in practical problems involving fissured elastic media with reasonable accuracy. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
20.
In natural rock masses, the shapes of three‐dimensional (3‐D) blocks cut by arbitrary fracture networks may be very complex. Owing to the geometric complexity and difficulty of mesh discretization of 3‐D blocks and fracture facets, explicit consideration of fracture networks in flow analysis of fractured porous medium (FPM) is very challenging. Using the numerical manifold method based on independent covers (NMMIC), an FPM flow model was proposed that can deal with very complex 3‐D fracture networks. In this paper, the convergence of NMMIC was first demonstrated. The theoretical basis of the arbitrary refinement of computational meshes was proven. Moreover, three peculiarities of NMMIC meshes, that is, arbitrary shape, arbitrary connection, and arbitrary refinement of independent covers, were concluded. Finally, some two‐dimensional (2‐D) tunnel flow examples were analyzed and the numerical results were compared with the analytical results. 3‐D examples with complex fracture distributions were also analyzed. In addition, the computational scale of the developed program was tested by increasing the number of computational elements. The results show that our model can accurately analyze the groundwater flow of rocks surrounding tunnels with complex fracture distributions. 相似文献