Pore scale study of permeability and tortuosity for flow through particulate media using Lattice Boltzmann method     

Ali Ghassemi  Ali Pak 《国际地质力学数值与分析法杂志》2011,35(8):886-901

In this paper, Lattice Boltzmann method (LBM) has been used to study the effects of permeability and tortuosity on flow through saturated particulate media and identify the relationships between permeability and tortuosity with other parameters such as particles diameter, grain specific surface, and porosity. LBM is a simple kinematic model that can incorporate the essential physics of microscopic and mesoscopic processes involved in flow through granular soils. The obtained results indicate that the 2D LB model, due to its inherent theoretical advantages, is capable of demonstrating that the porosity and specific surface are the most influential parameters in determining the intrinsic permeability of granular media. The obtained results show that particles' size diameter has a two‐fold effect on the coefficient of permeability: one is through specific surface and the other is by tortuosity factor. Numerical study also reveals that tortuosity of granular soils decreases almost linearly with increasing the porosity. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

致密砂岩水压瞬态致裂液化格子波尔兹曼及有限元数值研究     

朱伯靖  刘旭耀  程惠红  柳畅  刘善琪  任天翔  李永兵  石耀霖 《地震》2013,33(4):71-83

水压瞬态致裂液化在地应力测量、 地震破坏评估和机理研究、 油气(天然气、 页岩气)及地热资源开发等地学领域都具有重要的理论和应用价值。 自20世纪60年代至今在理论和实际应用方面取得了一些奠基性和开创性成果, 但鉴于问题复杂性, 瞬态致裂液化机理至今尚不清楚。 本研究应用格子波尔兹曼及有限元多孔介质流固耦合物理模型, 对地震波载荷作用下致密砂岩水压瞬态致裂液化过程进行数值模拟研究。 首先, 以鄂尔多斯盆地某油田延长组致密砂岩为例, 利用X射线CT断层成像技术, 应用基于量子力学第一性原理格子波尔兹曼方法, 建立致密砂岩数字岩芯模型。 进而, 推导格子波尔兹曼及有限元多孔介质流固耦合数值模型公式, 建立致密砂岩水压致裂液化物理模型。 最后, 应用地震波载荷边界和初始条件, 模拟瞬态致裂液化流固耦合过程, 讨论了地震波载荷幅值、 频率及作用时间对致密砂岩孔隙结构(孔隙度大小及连通性)、 致密砂岩破裂最大主应力之间关系, 得到了地震波作用下致密砂岩致裂液化准则。  相似文献   

Numerical experiments for the conductive properties of saturated rock     

《中国科学D辑(英文版)》2008,(Z2)

The reservoir evaluation as a key technology in oil exploration and production is based on the electrical transport property (ETP) of saturated rock that is described in a mathematical form with Arhcie’s equa-tion. But there have been increasing cases observed in many researches indicating that the ETP is non-Archie especially for the complex reservoir with low porosity and permeability. In this paper,the numerical experiments based on the Lattice Boltzmann method (LBM) have been employed to study the effect of porous structure and fluids on the ETP for revealing the nature of non-Archie phenomenon in micro-scale. The results of numerical experiments have proved that the saturation exponent n is a function of water saturation and porosity instead of being a constant in Archie’s equation. And then,a new formula has been developed for the EPT through combining the result of numerical simulation with that of laboratory measurements. The calculations from the new formula show very good agreement with laboratory measurements to demonstrate the efficiency of the new formula over the conventional methods in non-Archie rock.  相似文献   

Effect of macropore tortuosity and morphology on preferential flow through saturated soil: A Lattice Boltzmann study     

《Computers and Geotechnics》2014

Unlike micropores where water moves upward or downward based on hydraulic gradient, in macropores, water flows predominantly downward due to the gravity. Therefore, models based on capillary flow are not capable of simulating macropore flow. There are attempts to model the macropore flow using two domains, one for capillary flow and another one for macropores. These models use Richard’s equation for capillary flow and Poiseuille’s law for macropores in which the macropore is approximated to be cylindrical or planar. This study quantifies the magnitudes of the errors induced by this assumption. Influence of macropore shapes and tortuosity was quantified by using a 3D Lattice Boltzmann model, which is capable of simulating fluid flow in micropores as well as macropores of cracked clays. Artificial macropores of constant sectional area and volume, but different shapes were generated in 3D and the influence of macropore shapes, shape related parameters, and tortuosity were systematically investigated. Macropore flow rate decreases with different shapes compared to cylindrical macropores and increase in aspect ratio of sectional shape leads to decrease in macropore flow rate. The maximum effect of bends/turnings along the tortuous macropore was about 25% on overall decrease of flow rate due to tortuosity. However, more detailed study is required on the influence of bends on macropore flow rate. The macropore flow rate reduces by about 70% for tortuosity of 1.41. A prediction equation is verified to predict the flow rate of different shapes and tortuous macropores based on straight cylindrical macropore using aspect ratio and tortuosity factor.  相似文献   

Boundary effects on the Soil Water Characteristic Curves obtained from lattice Boltzmann simulations     

《Computers and Geotechnics》2016

  相似文献   

扩展SSH　Hubbard模型的晶格振动矩阵     

杨进 《成都信息工程学院学报》1996,(Z1)

在计及最近邻电子相互作用的扩展ＳＳＨＨｕｂｂａｒｄ模型时，得到了晶格振动矩阵，它对电子的正负自旋态对称和非对称的各类非线性元激发情形均适用。  相似文献   

孔隙尺度多孔介质流体流动与溶质运移高性能模拟          下载免费PDF全文

周鸿翔  郑延丰  吴劳生  陈铖  曾令藻 《水科学进展》2020,31(3):422-432

深入探究孔隙尺度下的流体流动特性和溶质运移规律对石油开采、农田养分管理、地下水污染修复有着重要意义。以人工构建的多孔介质结构和同步辐射X射线显微CT扫描的土壤团聚体(分辨率3.7μm)为研究对象,在空间节点数多达64 000 000的情况下,基于格子Boltzmann模型和GPU并行技术计算得到多孔介质流体运动和溶质运移过程的关键参数,并据此探究多孔介质空间异质性对水力学特性的影响。通过对3组不同结构的多孔介质比较发现,结构复杂程度最高的土壤样品和不规则堆叠的圆球结构的渗透率在100 mD(即10^-13m^2)量级,远低于规则堆叠的圆球结构(>20 000 mD);土壤的迂曲度为1.40~1.60,明显高于规则堆叠的圆球结构。研究结果表明,渗透率大的样品具有较小的迂曲度,这与结构的空间异质性有较强的关系;土壤的渗透率和迂曲度呈现各向异性;在水力梯度一定的前提下,渗透率较大的样品,纵向弥散系数也较大;同时,结构的异质性也会影响溶质的穿透曲线。本研究提出的模拟方法可在土壤结构中进行高效的水流运动和溶质运移模拟,可用于土壤多孔介质在孔隙尺度下的水力学特性研究。  相似文献   

Directional dependence of P- and S-wave propagation and polarization in foliated rocks from the Kola superdeep well: Evidence from laboratory measurements and calculations based on TOF neutron diffraction   总被引：1，自引：0，他引：1  

T.I. Ivankina  H.M. Kern  A.N. Nikitin 《Tectonophysics》2005,407(1-2):25-42

We have measured P- and S-wave velocities on two amphibolite and two gneiss samples from the Kola superdeep borehole as a function of pressure (up to 600 MPa) and temperature (up to 600 °C). The velocity measurements include compressional (Vp) and shear wave velocities (Vs1, Vs2) propagating in three orthogonal directions which were in general not parallel to inherent rock symmetry axes or planes. The measurements are accompanied by 3D-velocities calculations based on lattice preferred orientation (LPO) obtained by TOF (Time Of Flight) neutron diffraction analysis which allows the investigation of bulk volumes up to several cubic centimetres due to the high penetration depth of neutrons. The LPO-based numerical velocity calculations give important information on the different contribution of the various rock-forming minerals to bulk elastic anisotropy and on the relations of seismic anisotropy, shear wave splitting, and shear wave polarization to the structural reference frame (foliation and lineation). Comparison with measured velocities obtained for the three propagation directions that were not in accordance with the structural frame of the rocks (foliation and lineation) demonstrate that for shear waves propagating through anisotropic rocks the vibration directions are as important as the propagation directions. The study demonstrates that proper measurement of shear wave splitting by means of two orthogonal polarized sending and receiving shear wave transducers is only possible when their propagation and polarization directions are parallel and normal to foliation and lineation, respectively.  相似文献   

A scanning electron microscope study of the effects of dynamic recrystallization on lattice preferred orientation in olivine     

Kyung-Ho Lee  Zhenting Jiang  Shun-ichiro Karato   《Tectonophysics》2002,351(4)

Effects of dynamic recrystallization on lattice preferred orientation (LPO) in olivine were investigated through the combination of two SEM-based techniques, electron backscattered diffraction (EBSD) technique for crystallographic orientation measurement and backscattered electron imaging (BEI) for dislocation observation. Samples are experimentally deformed olivine aggregates in simple shear geometry. In the sample deformed at T=1473 K and high stresses (480 MPa), only incipient dynamic recrystallization is observed along grain-boundaries. Orientations of these small recrystallized grains are more random than that of relict grains, suggesting an important role of grain-boundary sliding at this stage of recrystallization. In the sample deformed at T=1573 K and low stress (160 MPa), dynamic recrystallization is nearly complete and the LPO is characterized by two [100] peaks. One peak is located at the orientation subparallel to the shear direction and is dominated by grains with high Schmid factor. The other occurs at high angles to the shear direction and is due to the contribution from grains with low Schmid factor. Grains with high Schmid factor tend to have higher dislocation densities than those with low Schmid factor. Based on these observations, we identify two mechanisms by which dynamic recrystallization affects LPO: (1) enhancement of grain-boundary sliding due to grain-size reduction, leading to the modification of LPO caused by the relaxation of constraint for deformation; (2) grain-boundary migration by which grains with lower dislocation densities grow at the expense of grains with higher dislocation densities. Based on the deformation mechanism maps and stress versus recrystallized grain-size relation, we suggest that the first mechanism always plays an important role whereas the second mechanism has an important effect only under limited conditions.  相似文献   

有限元离散模型中的出平面运动形式     

关慧敏  李小军 《华南地震》1993,13(2):1-6

波动在有限元离散模型中传播时,会出现许多不同于连续介质中波动的现象.本文引用经典物理学中晶格动力学的分析方法,由有限元离散形式的出平面波动方程得到频散关系,在此基础上,讨论了集中质量有限元模型中各种可能的出平面运动形式,并分析了考虑时域离散化影响下各种形式可能存在的频率范围.  相似文献