关于渗流中流线不封闭的特性和条件

本文对于流体在多孔介质中流动的特性进行理论研究和数值计算，提出两个关于渗流中流线不封闭的特性和条件，得到了在一般工程实际情况中的多孔介质区域内部不存在封闭流线的结论。本文以突变截面圆管中不可压缩渗流为算例，利用半人工瞬变方法进行数值计算，得到流体在充满多孔介质的突扩截面圆管和突缩截面圆管中流动时关于速度分布和压力分布的结果。由此表明，在突变截面附近的渗流区域中不存在回流和分离流，也不存在封闭的流线。渗流的这些流动特性不同于在无多孔介质的空间区域中的流动特性。     

裂缝性低渗透油藏流-固耦合理论与数值模拟

根据裂缝性低渗油藏的储层特征,建立适合裂缝性砂岩油藏渗流的等效连续介质模型。将渗流力学与弹塑性力学相结合,建立裂缝性低渗透油藏的流-固耦合渗流数学模型,并给出其数值解.通过数值模拟对一实际井网开发过程中孔隙度、渗透率的变化以及开发指标进行计算,并和刚性模型以及双重介质模型的计算结果进行了分析比较.     

可变形多孔介质渗透系数的测定方法

在Ｂｉｏｔ理论基础上给出可变形多孔介质耦合渗流基本方程；求出小试件一维定常耦合渗流问题的解答；表明在一维流固耦合情况下试件内部压力梯度有明显的非均匀性。因此通过实验确定可变形多孔介质渗透系数在数学上可归结为微分方程的反问题，传统的测试渗透系数的方法需要改进。介绍了可变形多孔介质渗透系数的测试原理和测试方法。对粒状多孔材料实验的结果表明，传统实验方法得到的渗透系数误差较大。     

基于物理界面的双重介质有限裂隙渗流模型

基于连续介质或者离散裂隙假设,含裂隙的多孔介质渗流问题有多种数学力学模型。受物理界面的启发,提出一种新的有限裂隙连续介质力学模型,可以为宏观裂隙-多孔介质内的流体输运问题等提供近似计算方案。该模型属于一类双重介质模型,将曲面上低维度的流场转化为三维空间的流场,并且与连续的多孔介质的流场耦合,在数学上表示为统一的输运控制方程和初始边界条件。这个近似模型为不方便实施高维度-低维度耦合求解的数值计算方法提供新的模拟思路,如光滑粒子流体动力学等无网格粒子类方法。     

考虑渗透系数变化的地下结构温度-渗流耦合分析

地下岩体的温度随着深度的增加呈逐渐上升的趋势，在地下水的作用下，岩体温度-渗流耦合作用明显。分析了地下岩体温度场与渗流场间的相互作用效应，考虑流体渗流和岩体骨架本身的热传导作用以及温度梯度所产生的渗流作用，并计入温度对水体运动粘滞系数及岩体渗透率的影响，给出了完整的基于连续介质模型的地下岩体温度-渗流耦合数学模型。借助...     

可压缩流体流动多尺度混合有限元数值方法研究

可压缩流体是天然油藏中广泛存在的一种流体,研究其在多孔介质中的渗流规律对于油藏开发具有重要意义。本文采用多尺度混合有限元方法,对可压缩流体渗流问题进行了研究。考虑流体的可压缩性以及介质形变,推导得到了可压缩流体渗流问题的多尺度计算格式。数值计算结果表明,多尺度混合有限元适于求解非均质性和可压缩流问题,具有节省计算量、计算精度高等优势,对于实际大规模油藏模拟具有重要意义。     

谈谈渗流

渗流是流体在多孔介质中的流动,渗流现象广泛地存在于自然界、工程材料、动物、植物中。多孔介质种类繁多,包括岩石(含各类矿藏)、土壤、生物材料和人工多孔介质材料等。渗流理论已经成为人类开发地下水、地热、石油、天然气、煤炭与煤层气等诸多地下资源的重要理论基础。本文从渗流的基本概念、渗流的分类、渗流的影响因素、渗流的特征以及渗流的研究意义等方面进行了阐述。     

论“岩体结构控制论”

经过长期实践和研究,作者于1984年提出岩体结构控制论是岩体力学的基础理论,并全面、系统地以岩体结构控制论为指导研究了岩体变形、岩体破坏及岩体力学性质的基本规律;提出岩体变形系山岩体材料变彤和岩体结构变形共同贡献的,岩体破坏系受岩体材料破坏和岩体结构破坏控制的;岩体力学性质不仅决定于岩体材料力学的性质,而且受控于岩体结构力学效应及环境因素力学效应。在此基础上,作者提出了岩体可以划分为连续介质、碎裂介质,块裂介质及板裂介质四种岩体力学介质,从而建立了完整的岩体结构力学理论体系。     

渗流力学的近况和展望

渗流力学研究流体在多孔介质的运动规律.渗流是指流体在多孔介质内的流动;孔隙介质、裂缝-孔隙介质以及各种类型的毛细管体系等均属多孔介质.渗流力学是流体力学与多孔介质理论和表面物理化学等学科交叉渗透产生的一个独立的学科领域,是多种工程技术的理论基础.     

裂隙岩体渗流模型研究现状与展望

裂隙岩体渗流对于边坡、地下工程及基础岩土体的承载能力有显著的制约作用。本文简要地介绍了多种裂隙岩体渗流模型研究现状,评述了几类比较有代表性的渗流模型特点以及存在的不足,为选取合理的数学模型用于求解具体的裂隙岩体渗流问题提供了参考依据,并在上述基础上提出了一些需要进一步研究的问题。指出就目前应用最为广泛的等效连续介质模型而言,裂隙岩体几何参数、有效孔隙度以及等效渗透张量的确定仍有待于更深入的研究,而从细观力学结构入手研究渗流耦合模型将具有重大的理论意义和实用价值。     

变开度岩体裂隙多相渗流实验与有效渗透率模型

岩体裂隙的有效渗透率是描述岩体非饱和或多相渗流的关键参数,而裂隙开度是影响有效渗透率的重要因素.通过自主研发的粗糙裂隙多相渗流可视化实验平台,针对天然岩体裂隙复制而成的裂隙模型开展变开度条件下的多相渗流可视化实验,研究开度变化对多相渗流流动结构以及有效渗透率的影响.研究表明:非湿润相流体运动通道,在低流量比条件下呈现出气泡流流动结构,而在高流量比条件下呈现较为稳定的通道流流动结构.随着开度的增加,非湿润相流动通道的分支变少、等效宽度增加,两相流体的有效渗透率均增大,流动结构趋于稳定.可视化结果还阐明了柱塞流流动结构下,两相流体交替占据裂隙空间的竞争机制:当非湿润相流体通道由连续转变为不连续时,裂隙进出口压差显著增加;反之,当该通道由不连续转变为连续时,压差显著减小.最后,基于分形理论以及渗透率统计建模方法,建立了考虑开度效应的岩体裂隙多相渗流有效渗透率理论模型,并通过实验测定的有效渗透率数据验证了该模型的正确性与有效性.     

裂隙岩体渗透系数确定方法研究

裂隙岩体渗透系数以及渗透主方向的确定对研究岩体渗透性大小及各向异性具有重要意义。高放废物地质处置库介质岩体的渗透性能将直接影响其使用安全性。本文运用离散裂隙网络模拟的方法对我国高放废物处置库甘肃北山预选区3#钻孔附近裂隙岩体进行了渗透性质分析。通过对3#钻孔1715～1780m段压水试验数据的反演,标定了离散裂隙网络渗流模型中的裂隙渗透参数(导水系数T)。利用标定的离散裂隙网络模型对场区裂隙岩体进行了渗流模拟,确定了该区域裂隙岩体的渗流表征单元体(REV)的尺寸大小以及渗透主值和主渗透方向。运用离散裂隙网络模型计算得出的渗透主值的几何均值与现场压水试验计算结果较接近,证明了计算结果的有效性。     

Geometric and Hydrodynamic Characteristics of Three-dimensional Saturated Prefractal Porous Media Determined with Lattice Boltzmann Modeling

Fractal and prefractal geometric models have substantial potential for contributing to the analysis of flow and transport in porous media such as soils and reservoir rocks. In this study, geometric and hydrodynamic parameters of saturated 3D mass and pore–solid prefractal porous media were characterized using the lattice Boltzmann model (LBM). The percolation thresholds of the 3D prefractal porous media were inversely correlated with the fraction of micro-pore clusters and estimated as 0.36 and 0.30 for mass and pore–solid prefractal porous media, respectively. The intrinsic permeability and the dispersivity of the 3D pore–solid prefractals were larger than those of the 3D mass prefractals, presumably because of the occurrence of larger solid and pore cluster sizes in the former. The intrinsic permeability and dispersivity of both types of structure increased with increasing porosity, indicating a positive relationship between permeability and dispersivity, which is at odds with laboratory data and current theory. This discrepancy may be related to limitations of the convection dispersion equation at the relatively high porosity values employed in the present study.     

Pore-Scale Simulation of Shear Thinning Fluid Flow Using Lattice Boltzmann Method

The present work attempts to identify the roles of flow and geometric variables on the scaling factor which is a necessary parameter for modeling the apparent viscosity of non-Newtonian fluid in porous media. While idealizing the porous media microstructure as arrays of circular and square cylinders, the present study uses multi-relaxation time lattice Boltzmann method to conduct pore-scale simulation of shear thinning non-Newtonian fluid flow. Variation in the size and inclusion ratio of the solid cylinders generates wide range of porous media with varying porosity and permeability. The present study also used stochastic reconstruction technique to generate realistic, random porous microstructures. For each case, pore-scale fluid flow simulation enables the calculation of equivalent viscosity based on the computed shear rate within the pores. It is observed that the scaling factor has strong dependence on porosity, permeability, tortuosity and the percolation threshold, while approaching the maximum value at the percolation threshold porosity. The present investigation quantifies and proposes meaningful correlations between the scaling factor and the macroscopic properties of the porous media.     

岩土介质的分形孔隙和分形粒子

岩土介质是晶粒状材料,存在大量的孔隙。这些孔隙的存在严重地影响岩土介质的力学、物理和化学性能。大量的研究表明,岩土介质的孔隙几何从原子尺度到晶粒尺寸范围内均表现出分形特征。目前分形几何已被广泛应用来研究岩土的孔隙率,输运特性和渗透性等等。本文从5个方面来讨论岩土介质的分形孔隙和分形粒子:①分形几何简介;②孔隙介质的分形模型;③岩土介质的分形孔隙特性和它们的分形量测方法;④岩土的分形粒子;⑤水土保持估计中的分形毛细管模型。      

膨胀性非牛顿流体多孔介质流动的模拟分析

在表征体元尺度采用格子Boltzmann方法分析膨胀性非牛顿流体在多孔介质中的流动,基于二阶矩模型在演化方程中引入表征介质阻力的作用力项,求解描述渗流模型的广义Navier-Stokes方程．采用局部法计算形变速率张量,通过循环迭代得到非牛顿粘度和松弛时间．对多孔介质的Poiseuille流动进行分析,通过比较发现结果与孔隙尺度的解析解十分吻合,并且收敛较快,表明方法合理有效．分析了渗透率和幂律指数对速度和压力降的影响,研究结果表明,膨胀性流体的多孔介质流动不符合达西规律,压力降的增加幅度小于渗透率的减小幅度．当无量纲渗透率Da小于10-5时,流道中的速度呈现均匀分布,并且速度分布随着幂律指数的减小趋于平滑．压力降随着幂律指数的增加而增加,Da越大幂律指数对压力降的影响越明显．     

Stress-induced fluid flow anisotropy in fractured rock

Anisotropic stress states are common in the upper crust and result in fracture apertures being dependent on fracture orientation. Fractured rocks should therefore display an anisotropic permeability determined by the aperture, length, and orientation of those fractures remaining open. In this paper, a numerical study of this effect is made for a rock containing two orthogonal fracture sets subject to a uniaxial compressive stress applied perpendicular to one of the sets. With increasing compressive stress, the decreasing aperture of fractures orientated perpendicular to the stress axis leads to a decrease in permeability both parallel and perpendicular to the stress. For flow parallel to the stress direction, this is a consequence of the finite length of the fractures, flow in fractures perpendicular to the stress being required to connect fractures orientated parallel to the stress direction. As the number of fractures is decreased towards the percolation threshold, the average permeability tensor is found to become increasingly isotropic. This behaviour results from the highly tortuous nature of the flow paths just at the percolation threshold.     

Low Pressure Gas Percolation Characteristic in Ultra-low Permeability Porous Media

Low pressure gas percolation characteristic in ultra-low permeability porous media is investigated in this article through core flow experiments. The results show that the wall-slip layer covers more than 10% of the average porous channel radius on account of minimum pore size when the permeability is below 0.1 × 10^?3μ m ² order, and seepage behavior is contrasted to that in mid-high permeability pore media. When the gas pressure is not high enough, the flow regime turns into transitional flow instead of slip flow, and nonlinear relationship between the measured gas permeability and the reciprocal of average pressure exists. The gas measuring permeability experiment would be influenced by the non-linear relationship. If Klinkenberg-corrected method is applied to speculate the equivalent liquid permeability, the extrapolated value will become less or minus. Simultaneously, actual gas flow velocity at the outlet is beyond the calculated value with Klinkenberg formula. A new gas seepage model based on the general slip boundary condition is derived from the homogenization technique in this article. At last the flow model is examined to be suitable for representing the gas flow behavior in ultra-low permeability media and estimating the absolute permeability from single-point, steady-states measurements.     

Fractal porous media I: Longitudinal stokes flow in random carpets

Two-dimensional porous media whose random cross-sections are derived from site percolation are constructed. The longitudinal flow of a Newtonian fluid in the Stokes approximation is then computed and the longitudinal permeability is obtained. Two methods are used and yield the same result when porosity is low. The Carman equation is shown to apply within ±7% when porosity is within the range from 0 to 0.75. Finally, random structures derived from stick percolation are investigated; results are qualitatively the same, but the Carman equation yields a poorer approximation.     

Investigation of Generalized Relative Permeability Coefficients for Electrically Assisted Oil Recovery in Oil Formations

Evaluation of relative permeability coefficients is one of the key steps in reliable simulation of two-phase flow in porous media. An extensive body of work exists on evaluation of these coefficients for two-phase flow under pressure gradient. Oil transport under an applied electrical gradient in porous media is also governed by the principles of two-phase flow, but is less understood. In this paper, relative permeability coefficients under applied electric field are evaluated for a specific case of two- phase fluid flow in water-wet porous media, where the second fluid phase is oil. It is postulated that the viscous drag on the oil phase, exerted by the electro-osmotic flow of the water phase, is responsible for the transport of oil in the absence of a pressure gradient. Reliable prediction of the flow patterns necessitates accurate representation and determination of the relative permeability coefficients under the electrical gradient. The contribution of each phase to the flow is represented mathematically, and the relative permeability coefficients are evaluated through electro-osmotic flow measurements conducted on oil bearing rock cores.