基于几何重建的气液两相流孔隙尺度模拟

本文基于二维micro CT扫描图像进行数值重构得到三维多孔介质几何模型,采用Shan-Chen多相格子Boltzmann方法(LBM)建立适用不规则几何边界区域内气液两相流动数值模型,在此基础上对复杂多孔介质内的两相流动过程进行孔隙尺度模拟。结果表明,气-液两相流动受限于多孔介质的复杂孔隙结构,气液表面张力、壁面润湿特性以及进出口压差对两相流动特性将产生明显影响。     

基于显微CT图像的岩芯孔隙分形特征

采用国内自主开发的高分辨率显微CT设备,对岩芯样本进行图像采集,运用数字图像分析方法获得图像中的岩芯孔隙目标,并以此为基础重建三维模型。然后基于分形理论,分别从二维图像和重建的三维模型角度计算孔隙的分形维数。最后,再采用基于多孔介质的分形模型对岩芯孔隙分形特征进行验证,得出岩芯的二维分形维数和三维分形维数之间的差值均值为1.000 3,最大偏差为0.004,很好地满足多孔介质分形模型中关于二维和三维分形维数的关系。     

水分子在微孔隙介质中的受限扩散模拟

在饱和多相流体孔隙介质中,利用核磁共振测得的扩散系数可以区分油气水,并用于饱和度等岩石物理参数的求取,但在微小孔隙(20μm)中,流体的表观扩散系数会偏离其自由扩散系数.通过有限差分和随机游走模拟方法对水分子在单孔隙及多孔隙介质中的扩散进行数值模拟,考察水分子扩散受限的程度.结合核磁共振在测井及实验室岩石物理应用的技术,进一步分析了受限扩散对求取油水饱和度等岩石物理参数造成的影响.同时,分析总结了两种模拟方法的使用条件及优缺点.     

基于数字图像相关的三维刚体位移测量方法

基于三维位移测量在工程技术领域的必要性和重要性,开展了单摄像机和数字图像相关相结合的三维刚体位移测试方法的研究.基于图像位移场矢量中心和斜率与面内和离面位移的分别对应关系,采用最小二乘拟合法分离图像位移场的常量项与一次项,据此,可实现物体三维位移分量的有效分离.以针孔摄像机成像模型为基础,开展了数值模拟及硅片平移实验,发展了与三维刚体位移对应的散斑图模拟方法,验证了基于仿射变换的相关迭代法的精度和适用性.数值模拟与实验结果验证了数字图像相关方法用于实现物体三维刚体化移重构的可行性和优越性,最大测量误差为5%.     

水分子在微孔隙介质中的受限扩散模拟

在饱和多相流体孔隙介质中,利用核磁共振测得的扩散系数可以区分油气水,并用于饱和度等岩石物理参数的求取,但在微小孔隙（< 20 μm）中,流体的表观扩散系数会偏离其自由扩散系数．通过有限差分和随机游走模拟方法对水分子在单孔隙及多孔隙介质中的扩散进行数值模拟,考察水分子扩散受限的程度．结合核磁共振在测井及实验室岩石物理应用的技术,进一步分析了受限扩散对求取油水饱和度等岩石物理参数造成的影响．同时,分析总结了两种模拟方法的使用条件及优缺点．     

基于LBM的多孔陶瓷有效热导率和相变传热分析

多孔氮化铝陶瓷具有高导热、耐熔盐腐蚀等优点在相变材料封装方面具有较大应用前景。本文采用微米计算机断层扫描(CT)得多孔氮化铝的灰度图像,并精准重构三维孔隙结构。并采用四参数随机生长、堆积颗粒、维诺结构等数值方法建立结构模型。基于介观尺度的格子玻尔兹曼方法 (LBM),计算了多孔陶瓷的有效热导率,结果表明维诺结构与CT真实重构的热导率接近。采用三维焓法LBM模拟了多孔骨架内的相变传热过程,获得了孔隙尺度对流相变的传热机理,并分析了不同孔隙率时的储热密度和复合材料的有效热导率。基于CT扫描和LBM的相变传热模拟有助于快速分析骨架的热特性,为多孔骨架的设计提供理论基础。     

碎片超高速撞击防护结构粒子场三维重构

为模拟空间碎片超高速撞击航天器防护结构表面材料喷射/溅射粒子场演化过程,并获取粒子场相关物理信息,基于粒子场同轴激光全息图像开展了碎片撞击过程的三维重构技术研究.首先对全息图像进行边缘剪切和缩放,将其划分为分辨率300×300左右的子图像以便于进行网格剖分;对于粒子堆叠区域子图像,采用基于三角化的网格剖分算法;剖分后形成的单一粒子采用Sobel算子提取其二维轮廓,然后将其投影到特定的三维空间形成三维形体;基于MAXScript语言实现了粒子场演化过程模拟.重构结果表明,无论是粒子场静态三维重构结果还是其演化过程均与撞击试验全息图像吻合较好,从而验证了该重构技术的有效性,为研究空间碎片对航天器防护结构的损伤效应提供了一种新的思路.     

碎片超高速撞击防护结构粒子场三维重构

为模拟空间碎片超高速撞击航天器防护结构表面材料喷射/溅射粒子场演化过程,并获取粒子场相关物理信息,基于粒子场同轴激光全息图像开展了碎片撞击过程的三维重构技术研究.首先对全息图像进行边缘剪切和缩放,将其划分为分辨率300×300左右的子图像以便于进行网格剖分;对于粒子堆叠区域子图像,采用基于三角化的网格剖分算法;剖分后形成的单一粒子采用Sobel算子提取其二维轮廓,然后将其投影到特定的三维空间形成三维形体;基于MAXScript语言实现了粒子场演化过程模拟.重构结果表明,无论是粒子场静态三维重构结果还是其演化过程均与撞击试验全息图像吻合较好,从而验证了该重构技术的有效性,为研究空间碎片对航天器防护结构的损伤效应提供了一种新的思路.     

数字重建孔隙结构内流体流动特性的计算研究

数字岩芯技术可重建储层介质的复杂孔隙结构,是从微观尺度上研究多孔介质内渗流过程及渗流机理的重要手段。本文基于真实岩芯的CT扫描图像,通过降噪、阈值分割、二值化处理获得岩芯CT图像的三维数字矩阵,并根据矩阵信息对真实岩芯孔隙结构进行重建,并输出为后续处理软件易读取的通用文件格式。利用Fluent软件对重建岩芯结构中的单相流动现象进行了数值模拟,对流体在复杂通道结构内的流动特性进行了详细分析,并对岩芯试样的宏观渗流率特性进行了预测。计算结果表明,本文采用的方法可高效重建孔隙结构,并对真实岩芯的宏观渗透率特性做出合理预测。     

基于ADI-FDTD模拟的岩石核磁共振横向弛豫特征分析

将复杂的骨架-孔隙系统抽象成等效双孔介质,根据Bloch方程构建数学模型,用交替隐式时域有限差分(ADI-FDTD)和联合反演迭代法(SIRT)进行横向宏观磁化矢量的数值模拟与核磁共振T₂谱的反演,定量研究扩散系数、弛豫速率、孔隙组分比和孔隙宽度对核磁响应的影响.结果表明:横向宏观磁化矢量衰减速率与扩散系数和微孔隙分量成正比,与孔隙宽度成反比,与表面弛豫速率基本无关.当扩散系数较大、孔隙宽度较小时,核磁共振T₂谱难以直观反映孔隙组分及孔隙结构.应用核磁共振评价孔隙结构时需特别注意扩散系数和孔隙尺寸的影响.     

基于图象处理和傅里叶变换的三维多孔介质重构方法

本文基于傅里叶变换的截断高斯场方法的理论,通过从实际多孔介质二维切片的二元化图像中提取孔隙率和自相关函数来构造整个三维多孔结构.模拟结果显示通过傅里叶变换技术得到的随机过程确实满足高斯分布,且构造出的三维多孔结构的统计特性基本与二维切片类似,说明了该技术能够构造出符合要求的三维多孔介质,从而为进一步从孔隙尺度研究其内部传热传质过程奠定了基础.     

Research on the reconstruction method of porous media using multiple-point geostatistics

The pore structural characteristics have been the key to the studies on the mechanisms of fluids flow in porous media. With the development of experimental technology, the modern high-resolution equipments are capable of capturing pore structure images with a resolution of microns. But so far only 3D volume data of millimeter-scale rock samples can be obtained losslessly. It is necessary to explore the way of virtually reconstructing larger volume digital samples of porous media with the representative stru...     

Direct pore-level modeling of incompressible fluid flow in porous media

We present a dynamic particle-based model for direct pore-level modeling of incompressible viscous fluid flow in disordered porous media. The model is capable of simulating flow directly in three-dimensional high-resolution micro-CT images of rock samples. It is based on moving particle semi-implicit (MPS) method. We modify this technique in order to improve its stability for flow in porous media problems. Using the micro-CT image of a rock sample, the entire medium, i.e., solid and fluid, is discretized into particles. The incompressible Navier–Stokes equations are then solved for each particle using the MPS summations. The model handles highly irregular fluid–solid boundaries effectively. An algorithm to split and merge fluid particles is also introduced. To handle the computational load, we present a parallel version of the model that runs on distributed memory computer clusters. The accuracy of the model is validated against the analytical, numerical, and experimental data available in the literature. The validated model is then used to simulate both unsteady- and steady-state flow of an incompressible fluid directly in a representative elementary volume (REV) size micro-CT image of a naturally-occurring sandstone with 3.398 μm resolution. We analyze the quality and consistency of the predicted flow behavior and calculate absolute permeability using the steady-state flow rate.     

Comparison of NMR simulations of porous media derived from analytical and voxelized representations

We develop and compare two formulations of the random-walk method, grain-based and voxel-based, to simulate the nuclear-magnetic-resonance (NMR) response of fluids contained in various models of porous media. The grain-based approach uses a spherical grain pack as input, where the solid surface is analytically defined without an approximation. In the voxel-based approach, the input is a computer-tomography or computer-generated image of reconstructed porous media. Implementation of the two approaches is largely the same, except for the representation of porous media. For comparison, both approaches are applied to various analytical and digitized models of porous media: isolated spherical pore, simple cubic packing of spheres, and random packings of monodisperse and polydisperse spheres. We find that spin magnetization decays much faster in the digitized models than in their analytical counterparts. The difference in decay rate relates to the overestimation of surface area due to the discretization of the sample; it cannot be eliminated even if the voxel size decreases. However, once considering the effect of surface-area increase in the simulation of surface relaxation, good quantitative agreement is found between the two approaches. Different grain or pore shapes entail different rates of increase of surface area, whereupon we emphasize that the value of the “surface-area-corrected” coefficient may not be universal. Using an example of X-ray-CT image of Fontainebleau rock sample, we show that voxel size has a significant effect on the calculated surface area and, therefore, on the numerically simulated magnetization response.     

The internal magnetic field distribution, and single exponential magnetic resonance free induction decay, in rocks

When fluid saturated porous media are subjected to an applied uniform magnetic field, an internal magnetic field, inside the pore space, is induced due to magnetic susceptibility differences between the pore-filling fluid and the solid matrix. The microscopic distribution of the internal magnetic field, and its gradients, was simulated based on the thin-section pore structure of a sedimentary rock. The simulation results were verified experimentally. We show that the 'decay due to diffusion in internal field' magnetic resonance technique may be applied to measure the pore size distribution in partially saturated porous media. For the first time, we have observed that the internal magnetic field and its gradients in porous rocks have a Lorentzian distribution, with an average gradient value of zero. The Lorentzian distribution of internal magnetic field arises from the large susceptibility contrast and an intrinsic disordered pore structure in these porous media. We confirm that the single exponential magnetic resonance free induction decay commonly observed in fluid saturated porous media arises from a Lorentzian internal field distribution. A linear relationship between the magnetic resonance linewidth, and the product of the susceptibility difference in the porous media and the applied magnetic field, is observed through simulation and experiment.     

Experiment and simulation study on construction of a three-dimensional network model

The construction of a network model is one of the key techniques in organic combination of microscopic flow experiment and simulation. The construction method of a three-dimensional network model is presented on the basis of CT scanning images in this paper. A series of CT slice images describing microscopic pore structure and fluid distribution of actual rock is obtained with the help of the industrial microfocus CT system. Based on the extraction of pore space skeleton, pore and throat information, the corresponding network model is established, and the conversion from three-dimensional CT image information to pore-throat size distribution and topological information is also achieved. The feature of this method lies in the fact that complicated pore space of rock may be characterized by pores and throats with a simple shape while keeping the geometry and flow characteristics. It is indicated that the calculated results of porosity, permeability, relative permeability curve and microscopic remaining oil distribution match very well the experimental results of water flooding and polymer flooding. This network model may fairly well characterize the rock microscopic pore-throat size and topological characteristics. Supported by the National Natural Science Foundation of China (Grant Nos. 10302021 and 10772200)     

Identification of radiative properties of reticulated ceramic porous inert media using ray tracing technique

The radiative properties of reticulated porous inert media are computationally identified using the real three-dimensional structural data of porous media. The computational grids data are reconstructed from three-dimensional computer tomography scans and magnetic resonance image scans of different reticulated porous media. A ray tracing algorithm is used to track the rays inside the grid structure. Statistically large numbers of rays are traced for their path length and incident angle, which are used to find the probability based equivalent extinction coefficient and scattering phase function. The equivalent extinction coefficients are found for porous media with different porosities and pore densities. The dependency of specular and diffuse scattering phase functions on the porous structure and surface reflectance are also studied.     

Understanding the micro structure of Berea Sandstone by the simultaneous use of micro-computed tomography (micro-CT) and focused ion beam-scanning electron microscopy (FIB-SEM)

Berea sandstone is the building block for reservoirs containing precious hydrocarbon fuel. In this study, we comprehensively reveal the microstructure of Berea sandstone, which is often treated as a porous material with interconnected micro-pores of 2-5 μm. This has been possible due to the combined application of micro-computed tomography (CT) and focused ion beam (FIB)-scanning electron microscopy (SEM) on a Berea sample. While the use of micro-CT images are common for geological materials, the clubbing and comparison of tomography on Berea with state-of-the-art microstructure imaging techniques like FIB-SEM reveals some unforeseen features of Berea microstructure. In particular, for the first time FIB-SEM has been used to understand the micro-structure of reservoir rock material like Berea sandstone. By using these characterization tools, we are able to show that the micro-pores (less than 30 μm) are absent below the solid material matrix, and that it has small interconnected pores (30-40 μm) and large crater-like voids (100-250 μm) throughout the bulk material. Three-dimensional pore space reconstructions have been prepared from the CT images. Accordingly, characterization of Berea sandstone specimen is performed by calculation of pore-structure volumes and determination of porosity values.