基于微焦点CT与数据约束模型的古代家猪牙齿显微结构研究

本文对一例考古发掘出土的古代家猪下颌骨及牙齿进行了微焦点CT成像,通过加置滤波片有效减弱了线束硬化对CT图像的影响。基于数据约束模型（DCM）对样品三维结构进行了定量计算,获取了包含小于CT图像体元尺寸在内的样品跨尺度结构信息。对样品中孔隙及高吸收组分的联通特性进行了跨尺度计算,得到这两种组分的三维团簇分布。与常见的图像阈值分割法相比,DCM模型有效减弱了部分体积效应对图像分析的影响。本文研究结果为考古样品的内部形态测量、结构分析提供了数据基础,为考古样品的数字建模存档提供一种新的方法。      

3D image-based characterization of fluid displacement in a Berea core

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.     

Accuracy of formation factors for three-dimensional pore-scale images of geo-materials estimated by renormalization technique

The formation factor, the dimensionless electric resistivity of porous rock/sediment saturated with conductive fluid, is an important quantity in geophysical exploration for petroleum reservoirs and groundwater aquifers. The renormalization technique is a promising approximation method for the quick estimation of the formation factors from large three-dimensional images of porous geo-materials obtained by X-ray microtomography. In the present study, we applied the renormalization technique to various pore-scale image sets of real geo-materials (sandstones, pumice, lava, and sandy sediments). The purpose is to explore the factors controlling the estimation accuracy of the formation factor. The results revealed that the accuracy increases with increasing (i) porosity, (ii) degree of the pore elongation along the direction of the applied electric field, and (iii) size of the initial subsystem with which the renormalization step starts and with decreasing (iv) pore or grain size. Most importantly, a high degree of elongation of the pore structure along the applied field ensures good accuracy even if the porosity is low, the initial subsystem is small, and the pore or grain size is large.     

Numerical estimation of carbonate rock properties using multiscale images

Characterizing the pore space of rock samples using three‐dimensional (3D) X‐ray computed tomography images is a crucial step in digital rock physics. Indeed, the quality of the pore network extracted has a high impact on the prediction of rock properties such as porosity, permeability and elastic moduli. In carbonate rocks, it is usually very difficult to find a single image resolution which fully captures the sample pore network because of the heterogeneities existing at different scales. Hence, to overcome this limitation a multiscale analysis of the pore space may be needed. In this paper, we present a method to estimate porosity and elastic properties of clean carbonate (without clay content) samples from 3D X‐ray microtomography images at multiple resolutions. We perform a three‐phase segmentation to separate grains, pores and unresolved porous phase using 19 μm resolution images of each core plug. Then, we use images with higher resolution (between 0.3 and 2 μm) of microplugs extracted from the core plug samples. These subsets of images are assumed to be representative of the unresolved phase. We estimate the porosity and elastic properties of each sample by extrapolating the microplug properties to the whole unresolved phase. In addition, we compute the absolute permeability using the lattice Boltzmann method on the microplug images due to the low resolution of the core plug images. In order to validate the results of the numerical simulations, we compare our results with available laboratory measurements at the core plug scale. Porosity average simulations for the eight samples agree within 13%. Permeability numerical predictions provide realistic values in the range of experimental data but with a higher relative error. Finally, elastic moduli show the highest disagreements, with simulation error values exceeding 150% for three samples.     

裂缝-孔隙介质储层渗透率表征及其影响因素分析

含裂缝多孔介质渗透率预测是非常规油气资源勘探开发的一个紧迫问题.现有多孔介质岩石物理模型通常利用圆形孔管模拟宏观岩石孔隙空间,难以定量描述软孔隙/裂缝在压力作用下的闭合情况,缺乏裂缝/孔隙间流量交换的连通机制.本文提出含三维裂缝/软孔隙网络多孔介质模型,将储层岩石裂缝/软孔隙表示为椭圆截面微管,建立了周期性压力作用下微...     

Texture and anisotropy analysis of Qusaiba shales

Scanning and transmission electron microscopy, synchrotron X‐ray diffraction, microtomography and ultrasonic velocity measurements were used to characterize microstructures and anisotropy of three deeply buried Qusaiba shales from the Rub’al‐Khali basin, Saudi Arabia. Kaolinite, illite‐smectite, illite‐mica and chlorite show strong preferred orientation with (001) pole figure maxima perpendicular to the bedding plane ranging from 2.4–6.8 multiples of a random distribution (m.r.d.). Quartz, feldspars and pyrite crystals have a random orientation distribution. Elastic properties of the polyphase aggregate are calculated by averaging the single crystal elastic properties over the orientation distribution, assuming a nonporous material. The average calculated bulk P‐wave velocities are 6.2 km/s (maximum) and 5.5 km/s (minimum), resulting in a P‐wave anisotropy of 12%. The calculated velocities are compared with those determined from ultrasonic velocity measurements on a similar sample. In the ultrasonic experiment, which measures the effects of the shale matrix as well as the effects of porosity, velocities are smaller (P‐wave maximum 5.3 km/s and minimum 4.1 km/s). The difference between calculated and measured velocities is attributed to the effects of anisotropic pore structure and to microfractures present in the sample, which have not been taken into account in the matrix averaging.     

Micro-computed tomography pore-scale study of flow in porous media: Effect of voxel resolution

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.     

Automated contact angle estimation for three-dimensional X-ray microtomography data

Multiphase flow in capillary regimes is a fundamental process in a number of geoscience applications. The ability to accurately define wetting characteristics of porous media can have a large impact on numerical models. In this paper, a newly developed automated three-dimensional contact angle algorithm is described and applied to high-resolution X-ray microtomography data from multiphase bead pack experiments with varying wettability characteristics. The algorithm calculates the contact angle by finding the angle between planes fit to each solid/fluid and fluid/fluid interface in the region surrounding each solid/fluid/fluid contact point. Results show that the algorithm is able to reliably compute contact angles using the experimental data. The in situ contact angles are typically larger than flat surface laboratory measurements using the same material. Wetting characteristics in mixed-wet systems also change significantly after displacement cycles.     

Using synchrotron X-ray microtomography to characterize the pore network of reservoir rocks: A case study on carbonates

In this work we propose a new methodology to calculate pore connectivity in granular rocks. This method is useful to characterize the pore networks of natural and laboratory compaction bands (CBs), and compare them with the host rock pore network.Data were collected using the synchrotron X-ray microtomography technique and quantitative analyses were carried out using the Pore3D software library. The porosity was calculated from segmented tridimensional images of deformed and pristine rocks. A process of skeletonization of the pore space was used to obtain the number of connected pores within the rock volume. By analyzing the skeletons the differences between natural and laboratory CBs were highlighted. The natural CB has a lower porosity than to the laboratory one. In natural CBs, the grain contacts appear welded, whereas laboratory CBs show irregular pore shape. Moreover, we assessed for the first time how pore connectivity evolves as a function of deformation, documenting the mechanism responsible for pore connectivity drop within the CBs.     

Micro-scale experimental investigation of the effect of flow rate on trapping in sandstone and carbonate rock samples

We present the results of a pore-scale experimental study of residual trapping in consolidated sandstone and carbonate rock samples under confining stress. We investigate how the changes in wetting phase flow rate impacts pore-scale distribution of fluids during imbibition in natural, water-wet porous media. We systematically study pore-scale trapping of the nonwetting phase as well as size and distribution of its disconnected globules. Seven sets of drainage-imbibition experiments were performed with brine and oil as the wetting and nonwetting phases, respectively. We utilized a two-phase miniature core-flooding apparatus integrated with an X-ray microtomography system to examine pore-scale fluid distributions in small Bentheimer sandstone (D = 4.9 mm and L = 13 mm) and Gambier limestone (D = 4.4 mm and L = 75 mm) core samples. The results show that with increase in capillary number, the residual oil saturation at the end of the imbibition reduces from 0.46 to 0.20 in Bemtheimer sandstone and from 0.46 to 0.28 in Gambier limestone. We use pore-scale displacement mechanisms, in-situ wettability characteristics, and pore size distribution information to explain the observed capillary desaturation trends. The reduction was believed to be caused by alteration of the order in which pore-scale displacements took place during imbibition. Furthermore, increase in capillary number produced significantly different pore-scale fluid distributions during imbibition. We explored the pore fluid occupancies and studied size and distribution of the trapped oil clusters during different imbibition experiments. The results clearly show that as the capillary number increases, imbibition produces smaller trapped oil globules. In other words, the volume of individual trapped oil globules decreased at higher brine flow rates. Finally, we observed that the pore space in the limestone sample was considerably altered through matrix dissolution at extremely high brine flow rates. This increased the sample porosity from 44% to 62% and permeability from 7.3 D to 80 D. Imbibition in the altered pore space produced lower residual oil saturation (from 0.28 to 0.22) and significantly different distribution of trapped oil globules.     

Application of GA in optimization of pore network models generated by multi-cellular growth algorithms

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.     

Computer simulation of flow through a lattice flow-cell model

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).     

基于微CT技术的砂岩数字岩石物理实验

数字岩石物理技术可弥补传统岩石物理实验的诸多不足,为岩石物理学研究提供一个新平台.本文以常规砂岩为研究对象,利用微CT扫描结合先进的图像处理技术建立了具有真实孔隙结构特征的三维数字岩芯模型;应用Avizo软件内含的多种形态学算法进行数字岩芯孔隙结构量化及表征研究,统计获取了孔隙度、孔隙体积分布及孔径分布特征,建立了等价孔隙网络模型;将Avizo与多场耦合有限元软件Comsol完美对接,实现了孔隙尺度的渗流模拟并计算获得绝对渗透率,对于考虑固相充填孔隙的情况,模拟计算了岩石有效弹性参数,并与近似Gassmann方程良好验证.本文所提出的将Avizo与Comsol结合使用的方法丰富了现有的数字岩石物理研究手段,为其大规模发展提供了一条新途径.     

Direct numerical simulation of pore-scale flow in a bead pack: Comparison with magnetic resonance imaging observations

A significant body of current research is aimed at developing methods for numerical simulation of flow and transport in porous media that explicitly resolve complex pore and solid geometries, and at utilizing such models to study the relationships between fundamental pore-scale processes and macroscopic manifestations at larger (i.e., Darcy) scales. A number of different numerical methods for pore-scale simulation have been developed, and have been extensively tested and validated for simplified geometries. However, validation of pore-scale simulations of fluid velocity for complex, three-dimensional (3D) pore geometries that are representative of natural porous media is challenging due to our limited ability to measure pore-scale velocity in such systems. Recent advances in magnetic resonance imaging (MRI) offer the opportunity to measure not only the pore geometry, but also local fluid velocities under steady-state flow conditions in 3D and with high spatial resolution. In this paper, we present a 3D velocity field measured at sub-pore resolution (tens of micrometers) over a centimeter-scale 3D domain using MRI methods. We have utilized the measured pore geometry to perform 3D simulations of Navier–Stokes flow over the same domain using direct numerical simulation techniques. We present a comparison of the numerical simulation results with the measured velocity field. It is shown that the numerical results match the observed velocity patterns well overall except for a variance and small systematic scaling which can be attributed to the known experimental uncertainty in the MRI measurements. The comparisons presented here provide strong validation of the pore-scale simulation methods and new insights for interpretation of uncertainty in MRI measurements of pore-scale velocity. This study also provides a potential benchmark for future comparison of other pore-scale simulation methods. © 2012 Elsevier Science. All rights reserved.     

Fractal generation of surface area of porous media

Many natural porous geological rock formations, as well as engineered porous structures, have fractal properties, i.e., they are self-similar over several length scales. While there have been many experimental and theoretical studies on how to quantify a fractal porous medium and on how to determine its fractal dimension, the numerical generation of a fractal pore structure with predefined statistical and scaling properties is somewhat scarcer. In the present paper a new numerical method for generating a three-dimensional porous medium with any desired probability density function (PDF) and autocorrelation function (ACF) is presented. The well-known Turning Bands Method (TBM) is modified to generate three-dimensional synthetic isotropic and anisotropic porous media with a Gaussian PDF and exponential-decay ACF. Porous media with other PDF's and ACF's are constructed with a nonlinear, iterative PDF and ACF transformation, whereby the arbitrary PDF is converted to an equivalent Gaussian PDF which is then simulated with the classical TBM. Employing a new method for the estimation of the surface area for a given porosity, the fractal dimensions of the surface area of the synthetic porous media generated in this way are then measured by classical fractal perimeter/area relationships. Different 3D porous media are simulated by varying the porosity and the correlation structure of the random field. The performance of the simulations is evaluated by checking the ensemble statistics, the mean, variance and ACF of the simulated random field. For a porous medium with Gaussian PDF, an average fractal dimension of approximately 2.76 is obtained which is in the range of values of actually measured fractal dimensions of molecular surfaces. For a porous medium with a non-Gaussian quadratic PDF the calculated fractal dimension appears to be consistently higher and averages 2.82. The results also show that the fractal dimension is neither strongly dependent of the porosity nor of the degree of anisotropy assumed.     

Fractal generation of surface area of porous media

Many natural porous geological rock formations, as well as engineered porous structures, have fractal properties, i.e., they are self-similar over several length scales. While there have been many experimental and theoretical studies on how to quantify a fractal porous medium and on how to determine its fractal dimension, the numerical generation of a fractal pore structure with predefined statistical and scaling properties is somewhat scarcer. In the present paper a new numerical method for generating a three-dimensional porous medium with any desired probability density function (PDF) and autocorrelation function (ACF) is presented. The well-known Turning Bands Method (TBM) is modified to generate three-dimensional synthetic isotropic and anisotropic porous media with a Gaussian PDF and exponential-decay ACF. Porous media with other PDF's and ACF's are constructed with a nonlinear, iterative PDF and ACF transformation, whereby the arbitrary PDF is converted to an equivalent Gaussian PDF which is then simulated with the classical TBM. Employing a new method for the estimation of the surface area for a given porosity, the fractal dimensions of the surface area of the synthetic porous media generated in this way are then measured by classical fractal perimeter/area relationships. Different 3D porous media are simulated by varying the porosity and the correlation structure of the random field. The performance of the simulations is evaluated by checking the ensemble statistics, the mean, variance and ACF of the simulated random field. For a porous medium with Gaussian PDF, an average fractal dimension of approximately 2.76 is obtained which is in the range of values of actually measured fractal dimensions of molecular surfaces. For a porous medium with a non-Gaussian quadratic PDF the calculated fractal dimension appears to be consistently higher and averages 2.82. The results also show that the fractal dimension is neither strongly dependent of the porosity nor of the degree of anisotropy assumed.     

碳酸盐岩储层孔隙结构对电阻率的影响研究

碳酸盐岩储层孔隙类型多样,各种孔隙的尺寸变化范围可以跨越几个数量级,孔隙结构非常复杂,这种复杂孔隙结构和不均匀分布的多元孔隙空间使得储层电性呈现明显非阿尔奇特性.为了了解影响电阻率变化的控制因素,本次研究选取中三叠世雷口坡组的8块全直径碳酸盐岩岩样,开展了核磁共振、岩电实验、孔渗实验、压汞实验及薄片等实验,并利用数字图像分析法定量分析了孔隙结构特征.研究结果表明：①孔隙度是影响电阻率高低的重要因素,但并非唯一因素,除孔隙度以外,孔隙尺寸和数量、孔隙网络复杂程度远比吼道大小对电阻率的影响大;②在孔隙度一定的条件下,胶结指数m随储层中孤立大孔隙占比的增多而增大,当孔隙度增大到一定程度后,胶结指数m又随大孔隙占比的增多而减小,微裂缝起重要沟通作用;③在给定孔隙度时,以简单大孔隙为主的岩样表现为胶结指数m值较大,而以复杂孔隙网络、细小孔隙为主的岩样表现为胶结指数m值较小,具分散、孤立大孔隙的岩样,胶结指数m值最高;④依据孔隙几何参数与电阻率和胶结指数之间的关系,可以利用测井资料间接判别储层类型,从而提高储层有效性和含水饱和度评价精度.     

Online service for monitoring the ionosphere based on data from the global navigation satellite system

A service is described that makes possible the effective construction of a three-dimensional ionospheric model based on the data of ground receivers of signals from global navigation satellite positioning systems (GNSS). The obtained image has a high resolution, mainly because data from the IPG GNSS network of the Federal Service for Hydrometeorology and Environmental Monitoring (Rosgidromet) are used. A specially developed format and its implementation in the form of SQL structures are used to collect, transmit, and store data. The method of high-altitude radio tomography is used to construct the three-dimensional model. The operation of all system components (from registration point organization to the procedure for constructing the electron density three-dimensional distribution and publication of the total electron content map on the Internet) has been described in detail. The three-dimensional image of the ionosphere, obtained automatically, is compared with the ionosonde measurements, calculated using the two-dimensional low-altitude tomography method and averaged by the ionospheric model.     

砾岩储层地震波传播方程:三重孔隙结构模型

针对砾岩储层的砂、砾、泥三重孔隙结构特征,本文分析砾岩孔隙区域、砂岩孔隙区域以及泥岩孔隙区域相互之间的孔隙流体流动机制,将静态的砾岩骨架本构方程与动态的孔隙流体运动方程联立,提出了复杂砾岩储层的弹性波传播理论方程.采用实测砾岩储层参数,在算例中与双重孔隙介质理论进行对比分析,验证了本文理论方程的合理性;基于三重孔隙介质模型,分析不同储层环境下纵波的传播特征,结果显示:随流体黏滞系数增大,在衰减-频率轴坐标系中,砾与砂、砂与泥孔隙区域间局域流导致的两个衰减峰向低频端移动,而Biot全局流导致的衰减峰向高频端移动;嵌入体尺寸及背景相介质渗透率的变化,主要影响纵波速度频散曲线沿频率轴左、右平移,不影响波速低频、高频极限幅值;嵌入体含量及孔隙度的变化改变了岩石干骨架的弹性、密度参数,不仅影响速度频散曲线沿频率轴平移,而且影响其上、下限幅值;砾包砂包泥三重孔隙介质模型所预测的衰减曲线中,低频段"第一个衰减峰"主要由砾岩孔隙区域与砂岩孔隙区域之间的局域流导致,中间频段"第二个衰减峰"主要由砂岩孔隙区域与泥岩孔隙区域之间的局域流导致,超声频段"第三个衰减峰"由Biot全局流导致.对慢纵波传播特征的分析显示,砂岩骨架(局部孔隙度较大)内部的宏观孔隙流体流动造成的耗散明显强于砾岩与泥岩骨架.     

Dual models of pore spaces

We present a model for pore spaces that consists of two parts related by duality: (1) a decomposition of an open polyhedral pore space into open contractible pore bodies separated by relatively open interfaces and (2) a pore network that is homotopy equivalent to the pore space. The dual model is unique and free of parameters, but it relies on regularity conditions for the pore space. We show how to approximate any pore space by the interior of a polyhedral complex such that the regularity conditions are fulfilled. Thus, we are able to calculate the dual model from synthetic porous media and images of real porous media. The pore bodies are unions of relatively open Delaunay cells with respect to the corners of the pore boundary, and the pore network consists of certain at most two-dimensional (2D) Voronoi cells with respect to the corners of the pore boundary. The pore network describes the neighborhood relations between the pore bodies. In particular, any relatively open 2D Delaunay face f separating two pore bodies has a unique (relatively open) dual network edge. In our model, f is a pore throat only if it is hit by its dual network edge. Thus, as opposed to widespread intuition, any pore throat is convex, and adjacent pore bodies are not necessarily separated by pore throats. Due to the duality between the pore network and the decomposition of the pore space into pore bodies it is straightforward to store the geometrical properties of the pore bodies [pore throats] as attributes of the dual network vertices [edges]. Such an attributed network is used to perform 2D drainage simulations. The results agree very well with those from a pore-morphology based modeling approach performed directly on the digital image of a porous medium. Contractibility of the pore bodies and homotopy equivalence of the pore space and the pore network is proven using discrete Morse theory and the nerve theorem from combinatorial topology.