A reconstruction method for three‐dimensional pore space using multiple‐point geology statistic based on statistical pattern recognition and microstructure characterization

To predict the macroscopic properties (e.g., transport, electromagnetic, and mechanical properties) of porous media, it is necessary to have a three‐dimensional (3D) representation of porous media. We reconstruct the geologically realistic 3D structure of Fontainebleau sandstone based on the two‐dimensional (2D) thin sections by using the multiple‐point statistics method. For this method, the size of template is an important parameter that reflects the perceived scale of spatial structure of porous media. In this paper, we take advantage of entropy method to obtain the appropriate size of the template, which is proven to be correct and feasible. The reconstruction method proposed by us combines successive 2D MPS simulations as well as 3D MPS simulation, which takes account into the pore structure information (e.g., heterogeneity and connectivity) both intralayer and interlayer. This reconstruction method is tested on Fontainebleau sandstone for which 3D images from micro‐CT scanning are available. Applying local percolation theory analysis, this new approach can depict the expected patterns of geological heterogeneities. In addition, it also can well reproduce a high degree of connectivity of the pore space better than other reconstruction methods based on lower‐order statistics. Copyright © 2011 John Wiley & Sons, Ltd.     

High order discontinuous Galerkin method for simulating miscible flooding in porous media

We present a high-order method for miscible displacement simulation in porous media. The method is based on discontinuous Galerkin discretization with weighted average stabilization technique and flux reconstruction post processing. The mathematical model is decoupled and solved sequentially. We apply domain decomposition and algebraic multigrid preconditioner for the linear system resulting from the high-order discretization. The accuracy and robustness of the method are demonstrated in the convergence study with analytical solutions and heterogeneous porous media, respectively. We also investigate the effect of grid orientation and anisotropic permeability using high-order discontinuous Galerkin method in contrast with cell-centered finite volume method. The study of the parallel implementation shows the scalability and efficiency of the method on parallel architecture. We also verify the simulation result on highly heterogeneous permeability field from the SPE10 model.     

饱和多孔介质材料的应变局部化萌生条件

在单相介质和非渗流饱和多孔介质应变局部化萌生条件的基础上,应用饱和多孔介质控制方程和Liapunov稳定理论,导出了渗流条件下的固相应力-应变描述和有效应力-应变描述的多孔介质固相部分的应变局部化的萌生条件。不同应力描述下的萌生条件的形式有一定变化。应用简单算例,讨论了Terzaghi有效应力描述的应变局部化萌生条件中两种固、液相对运动特例下的饱和多孔介质应变局部化破坏的形式。     

考虑孔隙变形的孔隙介质本构关系初探

孔隙存在是孔隙介质材料结构的本质特征,它决定了孔隙介质的本构关系呈现更为复杂的性质。首先给出了孔隙介质的视应力、视应变、骨架实际应力和骨架实际应变的定义,并在理想孔隙介质的物质模型假定下讨论了理想孔隙介质的本构关系,表明即使孔隙介质的骨架结构满足虎克弹性体本构关系的假定,即假定骨架的实际应力与实际应变之间呈现线性关系,但由于孔隙的存在,孔隙介质的视应力与视应变之间还是呈现非线性,说明应力与变形之间的非线性是孔隙介质的固有特性。     

考虑固结效应的溶质传输研究

传统的孔隙介质水动力学采用对流-扩散方程,研究溶质在流体中的迁移。在这个过程中,孔隙介质被认为是不变形的,因而是一个稳态问题。针对二维情况下孔隙介质变形对溶质传输的影响,给出了考虑孔隙介质固结效应的溶质传输方程,并且探讨了该类问题的求解方法。     

断裂多孔岩石中的复杂反应流体动力学

:断裂多孔介质中的流体流动和组分输运过程对于理解金属成矿过程、油气运移、污染物质的输运和分布十分重要 ,近 2 0年来得到了广泛的研究。文中综述了这一领域的进展 ,涉及断裂多孔介质中的流体流动和组分输运的连续介质模型、多孔介质和断裂多孔介质的渗透率关系 ,特别是重点介绍了渗滤理论及其在断裂多孔介质中的流体流动和组分输运中的应用     

流体饱和两相多孔介质动力反应计算分析

基于流体饱和两相多孔介质的弹性波动方程组,运用显式逐步积分格式与局部透射人工边界相结合的时域显式有限元方法对该波动方程组进行求解,对两相介质在输入地震波作用下的弹性动力反应进行计算和分析;对在是否考虑孔隙流体渗流的两种情况下计算得到的两相介质弹性动力反应结果的差异进行对比研究,从而揭示孔隙流体渗流对两相介质动力反应性质的影响。计算结果表明：两相介质弹性动力反应时程的波形与入射地震波的波形相同,且弹性动力反应的峰值出现的时刻对应于入射地震波的峰值出现的时刻;孔隙流体的渗流将对两相介质的弹性动力反应性质产生显著的影响。数值计算同时表明,时域显式有限元方法是进行流体饱和两相多孔介质弹性动力反应计算分析的一种有效的方法。     

孔隙介质润湿性的核磁共振刻度特征的测量

核磁共振测量技术已能够快速而无损地获得孔隙介质的物性信息, 但对于颗粒表面润湿性的测量还处于定性到半定量的水平.采用毛细管法和玻璃板法对经不同浓度二甲基二氯硅烷溶液处理的玻璃表面润湿性测量, 然后对核磁共振方法测量的孔隙介质润湿性结果进行标定, 进而得到孔隙介质润湿性系列的刻度特征.实验表明, 对于玻璃颗粒孔隙介质, 核磁共振测量的结果在作为表面处理剂的二甲基二氯硅烷溶液浓度0 %~ 0.7%变化范围内才与孔隙介质润湿性的变化有着明确的对应关系.在此范围之外, 对应的驰豫时间的变化虽然较大, 但其对于润湿性的刻度已没有明确的作用.      

A developed technique for measuring water content in oil-contaminated porous media

Water content is an important physical parameter for soil, vadose zone, and porous aquifer. Accurate measurement of water content in oil-contaminated porous media is critical for the research on oil pollution process and remediation in soil and groundwater systems. In this study, an improved water content calculation formula for oil-contaminated porous media was proposed based on the theory of oven-drying method, and laboratory experiments were conducted to test the applicability and accuracy of the formula for several types of manually prepared oil-contaminated porous media with different water contents. Furthermore, the measuring method and calculation formula, which can be used to determine the water content of porous media sampled from the oil-contaminated sites, were proposed for the first time in this study based on the improved formula. The experimental results showed that the improved formula was very accurate when used to calculate the water contents of diesel-contaminated sand, gasoline-contaminated mild clay, and engine oil-contaminated sand, indicating that it was widely applicable to oils with different volatile ability as well as porous media with different texture. This study meets the urgent need for accurate determination of water content in oil-contaminated porous media, and it solves the technical problem that the existing water content measuring methods cannot be applied directly in the field study.     

地球化学作用下饱和-非饱和介质水力-传质-传热耦合模型

地球化学作用是影响岩土类多孔介质中流体流动、传质、传热的重要因素之一。基于Fredlund所提出的非饱和土四相理论,将地球化学作用所产生的溶解/沉淀视为一个独立相,在一定的假设基础上,利用溶质浓度对介质孔隙度、含水率进行修正,从而建立起传质过程与介质水力性质之间的关系,改进了传统模型中传质过程与流体密度、黏度之间单一的耦合关系;并以此为基础建立了考虑地球化学作用下饱和-非饱和介质中多组分、多相流体渗流场-浓度场-温度场耦合的数学模型。通过在实验室内对反应性煤矸石进行渗透实验的结果对笔者所建模型进行了验证,结果表明：综合考虑溶解作用时所测点处渗流速度与未考虑耦合作用的结果相比增大0.58倍;单独考虑溶解/沉淀相时,所测点处渗流速度与未考虑密度耦合作用时的结果相比增大0.44倍。水流的实际渗透速度与在使用考虑溶解/沉淀相计算模型时的结果吻合较好。     

An analytical solution for the transient response of saturated linear elastic porous media

The motions of fluid and solid phases in saturated porous media are coupled by inertial, viscous and mechanical interactions as described by Biot's equations. A one-dimensional exact analytical solution of the Biot's equations for the completely general solution of the transient problem in saturated, linear, elastic, porous media is presented. The problem is solved by using the Fourier series. The transient response of porous media is shown for typical material properties of a natural granular deposit and for different degrees of viscous coupling. The analytical results show the mechanics of dispersive wave propagation in saturated porous media and they should provide a useful comparison term for the existing numerical solution methods.     

冰点下多孔介质中甲烷水合物的生成特性

冰点下水合物在多孔介质的生成是一个复杂的多相转化过程,为了研究冰点下多孔介质中水合物生成过程的水相转化率、气体消耗量与稳定压力等生成特性,在定容条件下,进行不同孔径与粒径的多孔介质中甲烷水合物在冰点下的生成实验。所使用的多孔介质平均孔径为1295 nm、1796 nm和3320 nm。研究结果表明：水合物生成结束时水的转化率随着初始生成压力的增大而增大,随着温度的升高而降低,随着孔径的增大而增大;在相同的孔径下,多孔介质粒径的增大降低了水合物的生成速率但对最终气体消耗量没有影响;在相同的温度下,随着初始生成压力的增大,实验最终压力、气体消耗量与最终水的转化率均随之升高;温度越高,不同的生成初始生成压力下体系的最终稳定压力与水的转化率相差越大;在多孔介质的毛细管作用力与结合水的共同作用下,冰点下水合物生成的水的转化率会大大地降低。在本实验条件下,水相转化为水合物的比例最高为32.39%。     

On multicomponent gas diffusion and coupling concepts for porous media and free flow: a benchmark study

Computational Geosciences - Multicomponent gas transport in porous media and at the interface between porous media and free flow occurs in a wide range of technical and environmental systems....     

孔隙尺度多孔介质流体流动与溶质运移高性能模拟

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

Porous Structure Reconstruction Using Convolutional Neural Networks

The three-dimensional high-resolution imaging of rock samples is the basis for pore-scale characterization of reservoirs. Micro X-ray computed tomography (µ-CT) is considered the most direct means of obtaining the three-dimensional inner structure of porous media without deconstruction. The micrometer resolution of µ-CT, however, limits its application in the detection of small structures such as nanochannels, which are critical for fluid transportation. An effective strategy for solving this problem is applying numerical reconstruction methods to improve the resolution of the µ-CT images. In this paper, a convolutional neural network reconstruction method is introduced to reconstruct high-resolution porous structures based on low-resolution µ-CT images and high-resolution scanning electron microscope (SEM) images. The proposed method involves four steps. First, a three-dimensional low-resolution tomographic image of a rock sample is obtained by µ-CT scanning. Next, one or more sections in the rock sample are selected for scanning by SEM to obtain high-resolution two-dimensional images. The high-resolution segmented SEM images and their corresponding low-resolution µ-CT slices are then applied to train a convolutional neural network (CNN) model. Finally, the trained CNN model is used to reconstruct the entire low-resolution three-dimensional µ-CT image. Because the SEM images are segmented and have a higher resolution than the µ-CT image, this algorithm integrates the super-resolution and segmentation processes. The input data are low-resolution µ-CT images, and the output data are high-resolution segmented porous structures. The experimental results show that the proposed method can achieve state-of-the-art performance.     

孔隙尺度多孔介质流体流动与溶质运移高性能模拟

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

A finite element method for modeling coupled flow and deformation in porous fractured media

Modeling the flow in highly fractured porous media by finite element method (FEM) has met two difficulties: mesh generation for fractured domains and a rigorous formulation of the flow problem accounting for fracture/matrix, fracture/fracture, and fracture/boundary fluid mass exchanges. Based on the recent theoretical progress for mass balance conditions in multifractured porous bodies, the governing equations for coupled flow and deformation in these bodies are first established in this paper. A weak formulation for this problem is then established allowing to build a FEM. Taking benefit from recent development of mesh‐generating tools for fractured media, this weak formulation has been implemented in a numerical code and applied to some typical problems of hydromechanical coupling in fractured porous media. It is shown that in this way, the FEM that has proved its efficiency to model hydromechanical phenomena in porous media is extended with all its performances (calculation time, couplings, and nonlinearities) to fractured porous media. Copyright © 2015 John Wiley & Sons, Ltd.     

多孔介质油水两相k~s~p关系数学模型的实验研究

数学模型是研究相对渗透率与饱和度关系曲线的重要方法。采用自行开发设计的人工平面多孔介质模型,测定了相对渗透率与饱和度的关系曲线。多孔介质选择粒径为0.5~1mm、1~2mm的标准砂,纯净的水为湿润相,用3号苏丹红染色的93#汽油为非湿润相,组成多孔介质油水两相流动系统。采用Van Genuchten and Mualeum(VGM)和Brooks-Corey-Burdine(BCB)两种数学模型计算相对渗透率与饱和度的关系曲线,通过比较两种数学模型计算结果之间和模型计算结果与实测结果的差异以及模型的应用、多相渗流系统自身特征,得出VGM、BCB两种数学模型计算结果符合实际情况,VGM模型应用过程更为简便,但VGM模型具有一定适用条件;在砂性多孔介质中,BCB模型计算相对渗透率与饱和度关系曲线更准确。     

A 2D mixed fracture–pore seepage model and hydromechanical coupling for fractured porous media

A novel two-dimensional mixed fracture–pore seepage model for fluid flow in fractured porous media is presented based on the computational framework of finite-discrete element method (FDEM). The model consists of a porous seepage model in triangular elements bonded by unbroken joint elements, as well as a fracture seepage model in broken joint elements. The principle for determining the fluid exchange coefficient of the unbroken joint element is provided to ensure numerical accuracy and efficiency. The mixed fracture–pore seepage model provides a simple but effective tool for solving fluid flow in fractured porous media. In this paper, examples of 1D and 2D seepage flow in porous media and porous media with a single fracture or multiple fractures are studied. The simulation results of the model match well with theoretical solutions or results obtained by commercial software, which verifies the correctness of the mixed fracture–pore seepage model. Furthermore, combining FDEM mechanical calculation and the mixed fracture–pore seepage model, a coupled hydromechanical model is built to simulate fluid-driven dynamic propagation of cracks in the porous media, as well as its influence on pore seepage and fracture seepage.

     

考虑扩散路径的宾汉姆流体渗透注浆机制

浆液在多孔介质中的扩散路径对渗透扩散范围和注浆效果具有非常重要的影响。采用理论分析,以分形特征与宾汉姆流体在多孔介质中的渗流运动方程为基础,揭示了考虑扩散路径的宾汉姆流体渗透注浆机制,并利用团队前期开展的渗透注浆试验对其进行了验证。分析了多孔介质孔隙率、宾汉姆水泥浆液水灰比、多孔介质渗透系数、注浆压力、地下水压力等对扩散半径的影响变化规律。同时,基于Comsol Multiphysics平台,采用计算机编程技术二次开发得到了考虑扩散路径的宾汉姆流体渗透注浆机制的渗透注浆三维数值模拟程序,并以此开展了宾汉姆水泥浆液在多孔介质中渗透扩散形态效果的数值模拟。研究结果表明：与不考虑扩散路径的宾汉姆流体渗透注浆球形扩散公式获得的扩散半径理论计算值相比,采用考虑扩散路径的宾汉姆流体渗透注浆机制得到的扩散半径理论计算值更接近试验值。该研究成果可为实践注浆工程提供一定的理论支撑。