气体动理学格式模拟多孔介质流动的可行性

将气体动理学格式（GKS）拓展到模拟多孔介质内的低速渗流,并检验在孔隙尺度上模拟不可压缩低速流动的可行性与有效性.结果表明：GKS具有二阶空间精度,能够较精确地计算多孔介质的渗透率;相比于单松弛格子玻尔兹曼方法,GKS能够精确实现壁面无滑移边界条件,从而正确反映渗透率与黏性无关的特性;对于Berea砂岩切片结构中的复杂流动,模拟结果与实验吻合较好,能较精确地计算渗透率.给出GKS模拟达西渗流的马赫数选取准则,为研究多孔介质流动提供新的工具.     

多孔介质中流体输运的孔尺度SPH模拟

本文采用光滑粒子流体动力学(SPH)方法[1]对流体流过多孔介质的过程进行数值模拟,方法与程序用经典的Poiseuille流与多孔介质流动的Kozeny解进行了验证.PH进一步用于数值模拟跨膜流动,获得了跨膜流动中速度场随时间的演化过程.多孔膜孔径在1 靘到200 靘范围内,流体跨膜渗透速度的SPH结果与K-K方程的计算值吻合较好.文中还讨论了粒子数目与计算精度的关系.本文的计算结果表明SPH方法具有模拟多孔介质微流动的能力.     

分形多孔介质中流体流动的LBM模拟

基于分形理论构造了规则Sierpinski地毯结构来模拟多孔介质孔隙结构,采用Lattice Boltzmann Method(LBM)分别对不同孔隙率与进出口压差下的分形多孔介质中的流体流动进行了数值模拟.结果表明,分形多孔介质的某些流场结构也具有分形特征;体积流量和流场进出口压差呈线性关系,符合Darcy定律;流体穿过分形多孔介质的体积流量与孔隙率成指数函数关系.     

梯度穿缝型双孔隙率多孔材料的吸声性能

在传统单一孔隙率多孔材料中引入宏观尺度的周期性梯度穿缝结构设计,构造出梯度穿缝型双孔隙率多孔材料,其包含多孔材料基体微孔尺度与穿缝尺度两个尺度。采用分层等效的理论建模方法,将复杂梯度渐变问题变为多层均匀等效层叠加问题。针对不同特征尺寸的多孔材料薄层,分别采用低、高两种渗透率对比度双孔隙率理论,给出了其等效密度和动态压缩系数,再应用传递矩阵方法得到了相邻薄层之间的声压和质点速度传递关系并求得其表面声阻抗,从而建立了梯度穿缝型双孔隙率多孔材料的吸声理论模型。发展了多尺度材料声学有限元数值模型,在所考虑的100~3000 Hz频段范围内数值模拟结果完全吻合理论模型结果。理论与模拟分析了多尺度结构参数对双孔隙率多孔材料吸声性能的影响,结果表明引入多尺度梯度结构设计能够显著提高单一孔隙率多孔材料的吸声性能,且穿缝尺度比穿缝梯度影响更为显著;精细数值模拟获得的声压和能量密度分布云图揭示了多尺度结构设计的吸声增强机制。该工作可用于指导双孔隙率多孔材料的多尺度结构设计,从而提高多孔材料的中低频吸声性能。      

裂缝性介质多尺度深度学习模型

结合人工神经网络建立裂缝介质多尺度深度学习流动模型.基于一套粗网格和一套细网格,通过在粗网格上训练数据,多尺度神经网络能够以较少的自由度训练出准确的神经网络.并在粗网格上通过求解局部流动问题获得多尺度基函数,结合神经网络进一步得到精细网格的解.基于离散裂缝的流动方程可视为多层网络,网络层数依赖于求解时间步数.阐述裂缝介质多尺度机器学习数值计算格式的建立,介绍如何使用多尺度算法构建离散裂缝模型的多尺度基函数,并采用超样本技术进一步提高计算准确性.数值结果表明,多尺度有限元算法与机器学习结合是一种有效的流体流动模拟算法.     

声波远探测测井与地震缝洞精细评价技术研究*

地震勘探探测范围广、空间尺度大,数据采集、处理、解释体系完善,资料运用率高,但纵向分辨率低。声波远探测可实现井周数十米范围内的异常体探测,克服了常规测井方法探测深度浅的局限,具有很高的纵向、周向、径向分辨率,但偶极横波远探测具有180°方位不确定性,径向探测范围不及地震勘探。本研究将两者优势结合,通过对缝洞型地质异常体远探测和地震响应特征进行综合分析,约束远探测方位解释结果,精细刻画地质异常体形态及发育情况,完善远探测测井的横向延伸信息,有效提高远探测解释的精准度,最终形成声波远探测测井与地震一体化缝洞综合评价技术,提升缝洞油气藏钻遇率和开发效率。在塔河等缝洞发育地区进行应用,实现多尺度、全方位和高精度的地质异常体描述与解释,为解决油气勘探过程中由于地震分辨率不足造成的脱靶等现象提供了新的思路。     

多孔介质中含溶解反应的互溶驱替过程格子Boltzmann研究

采用格子Boltzmann方法研究了孔隙尺度下多孔介质内含流固溶解反应的互溶驱替过程,重点研究了被驱替流体与驱替流体黏性差异较大的情况下,溶解反应引起的多孔介质内部结构变化对驱替过程的影响;定量分析了不同达姆科勒数及佩克莱数下多孔介质孔隙率和驱替过程驱替效率随时间的演变.研究结果表明:达姆科勒数较大时,溶解反应的发生会在多孔介质内部生成虫洞,导致一部分被驱替流体不能被波及,驱替流体沿虫洞离开多孔介质,造成驱替效率的减少.在此基础上,随着达姆科勒数的增大,孔隙率变化越大,生成的虫洞越宽,最终驱替效率变大,但仍小于无溶解反应时的驱替效率;随着佩克莱数的增大,指进增长速度越快,孔隙率变化越小,驱替效率越小.     

随机多孔介质的蒙特卡罗重构与渗流特性模拟

根据多孔介质微观结构的分形尺度标度特征,采用蒙特卡罗方法分别重构随机多孔介质的微观颗粒和孔隙结构,并基于分形毛管束模型研究多尺度多孔介质的气体渗流特性,建立多孔介质微观结构和宏观渗流特性的定量关系。结果表明：分形蒙特卡罗重构的多孔介质微细结构接近真实介质结构,气体渗流特性的计算结果与格子玻尔兹曼模拟数据较为吻合; 多孔介质气体渗透率随着克努森数的增加而增大,孔隙分形维数对于气体渗流的微尺度效应具有显著影响,而迂曲度分形维数对于表观渗透率和固有渗透率的比值影响可以忽略。提出的分形蒙特卡罗方法具有收敛速度快且计算误差与维数无关的优点,有利于深入理解多尺度多孔介质的渗流机理。     

多孔介质中流体流动及扩散的耦合格子Boltzmann模型

多孔介质中高Péclet数和大黏性比下混溶流体的流动和扩散广泛存在于二氧化碳驱油、化工生产等工业过程中.用数值方法对该问题进行研究时,关键在于如何正确描述高Péclet数和大黏性比下多孔介质内流体的行为.为此,提出了一种基于多松弛模型和格子动理模型的耦合格子Boltzmann模型.通过Chapman-Enskog分析,证明该模型能有效求解不可压Navier-Stokes方程和对流扩散方程.数值结果表明,该模型不仅具有二阶精度和良好的稳健性,而且对于高Péclet数和大黏性比的问题具有良好的数值稳定性,为模拟此类问题提供了有效工具.     

多孔介质结构对燃烧室内流场和喷雾特性影响的数值研究

以多孔介质发动机为背景,用数值模拟方法考察气缸内加入多孔介质蓄热体后对燃烧室内湍流流场及混合气形成的影响.计算基于Antohe和Lage的适用于多孔介质的κ-ε模型,其中引入了Darcy项和Forchheimer项.建立了一种可描述多孔介质随机结构的简化模型.针对燃油喷雾在不同结构多孔介质内的流场进行了数值计算.计算表明,多孔介质结构对燃烧室内燃油液滴的蒸发和混合气的形成过程具有重要的影响.     

XFEM modeling of hydraulic fracture in porous rocks with natural fractures

Hydraulic fracture (HF) in porous rocks is a complex multi-physics coupling process which involves fluid flow, diffusion and solid deformation. In this paper, the extended finite element method (XFEM) coupling with Biot theory is developed to study the HF in permeable rocks with natural fractures (NFs). In the recent XFEM based computational HF models, the fluid flow in fractures and interstitials of the porous media are mostly solved separately, which brings difficulties in dealing with complex fracture morphology. In our new model the fluid flow is solved in a unified framework by considering the fractures as a kind of special porous media and introducing Poiseuille-type flow inside them instead of Darcy-type flow. The most advantage is that it is very convenient to deal with fluid flow inside the complex fracture network, which is important in shale gas extraction. The weak formulation for the new coupled model is derived based on virtual work principle, which includes the XFEM formulation for multiple fractures and fractures intersection in porous media and finite element formulation for the unified fluid flow. Then the plane strain Kristianovic-Geertsma-de Klerk (KGD) model and the fluid flow inside the fracture network are simulated to validate the accuracy and applicability of this method. The numerical results show that large injection rate, low rock permeability and isotropic in-situ stresses tend to lead to a more uniform and productive fracture network.     

压敏三重介质油藏压力响应特征

建立由基岩系统、裂缝系统和溶洞系统组成,并考虑溶洞渗透率随压降的增加呈指数减小的压敏三重介质油藏试井解释模型,采用隐式差分格式对考虑井筒储存和污染效应的情况进行了求解.研究表明:无因次渗透率模数导致压力及压力导数明显增大,窜流系数影响窜流阶段出现的早晚,弹性储容比影响压力导数曲线"凹陷"的宽度和深度,外边界条件的影响则和普通三重介质明显不同,表皮系数主要影响全部压力和压力导数早期的"凸起",而无因次渗透率模数则主要影响压力和压力导数的后期.     

随机四参数法生成多孔介质及渗流模拟

本文基于多孔介质内流动的广泛应用背景,模拟了多孔介质渗流。采用随机四参数生长法生成了多孔介质模型,并对模型中的微小固相进行了过滤处理。将随机四参数生长法生成的多孔介质结合格子玻尔兹曼方法D2Q9模型模拟了压力驱动的多孔介质中的渗流。通过将模拟结果与著名的CK理论公式的计算结果进行比较,发现在孔隙率高时误差较小(φ=0.9时,误差约为6%)。研究表明在高孔隙率时采用文中的方法得到的多孔介质渗透率及孔隙尺度的流动结果较为准确。     

Theoretical analysis of porous radiant burners under 2-D radiation field using discrete ordinates method

This paper describes a theoretical study to investigate the heat transfer characteristics of porous radiant burners (PRBs). In the present work, a 2-D rectangular model is used to solve the governing equations for porous medium and gas flow before the premixed flame to the exhaust gas. The gas and the solid phases are considered in non-local thermal equilibrium and combustion in the porous medium is modeled by considering a non-uniform heat generation zone. The homogeneous porous media, in addition to its convective heat exchange with the gas, may absorb, emit and scatter thermal radiation. The radiation effect in the gas flow is neglected but the conductive heat transfer is taken into account. In order to analyze the thermal characteristics of porous burners, the coupled energy equations for the gas and porous medium in steady condition are solved numerically and the discrete ordinates method (DOM) is used to obtain the distribution of radiative heat flux in the porous media. Finally, the effects of various parameters on the performance of porous radiant burners are examined. The present results are compared with some reported theoretical and experimental results by other investigators and good agreement is found.     

A hierarchical fracture model for the iterative multiscale finite volume method

An iterative multiscale finite volume (i-MSFV) method is devised for the simulation of multiphase flow in fractured porous media in the context of a hierarchical fracture modeling framework. Motivated by the small pressure change inside highly conductive fractures, the fully coupled system is split into smaller systems, which are then sequentially solved. This splitting technique results in only one additional degree of freedom for each connected fracture network appearing in the matrix system. It can be interpreted as an agglomeration of highly connected cells; similar as in algebraic multigrid methods. For the solution of the resulting algebraic system, an i-MSFV method is introduced. In addition to the local basis and correction functions, which were previously developed in this framework, local fracture functions are introduced to accurately capture the fractures at the coarse scale. In this multiscale approach there exists one fracture function per network and local domain, and in the coarse scale problem there appears only one additional degree of freedom per connected fracture network. Numerical results are presented for validation and verification of this new iterative multiscale approach for fractured porous media, and to investigate its computational efficiency. Finally, it is demonstrated that the new method is an effective multiscale approach for simulations of realistic multiphase flows in fractured heterogeneous porous media.     

On the structure of the steady-state flow velocity field near the interface between a homogeneous liquid and a Brinkman porous medium

The structure of a homogeneous liquid flow at the interface with a porous medium saturated with the same liquid has been studied in the boundary layer approximation. The porous medium is described by the Brinkman model. Self-similar equations of motion in the form of Blasius equation have been found, and their numerical solutions have been presented. An expression for the force exerted by the flow on the porous medium is derived.     

Direct measurement of porous media local hydrodynamical permeability using gas MRI.

The concept of hydraulic permeability is at the core of modeling single phase or multi-phase flow in heterogeneous porous media, as it is the spatial distribution of the permeability that primarily governs the behavior of fluid flow in the medium. To date, the modeling of fluid flow in porous media has been hampered by poor estimates of local permeability. Magnetic Resonance Imaging is well known for its ability to measure non-invasively the local density and flow rate of different fluids saturating porous media [1,2]. In this paper we demonstrate the first non-invasive method for the direct measurement of a single projection of the local permeability tensor of a porous medium using gas-phase MRI. The potential for three-dimensional imaging of the medium permeability is also discussed. The limitations of the method are listed and results are presented in a model porous medium as well as in a real oil reservoir rock.     

Influence of various shapes of nanoparticles on unsteady stagnation-point flow of Cu-H2O nanofluid on a flat surface in a porous medium: A stability analysis

The nanofluid and porous medium together are able to fulfill the requirement of high cooling rate in many engineering problems. So, here the impact of various shapes of nanoparticles on unsteady stagnation-point flow of Cu-H₂O nanofluid on a flat surface in a porous medium is examined. Moreover, the thermal radiation and viscous dissipation effects are considered. The problem governing partial differential equations are converted into self-similar coupled ordinary differential equations and those are numerically solved by the shooting method. The computed results can reveal many vital findings of practical importance. Firstly, dual solutions exist for decelerating unsteady flow and for accelerating unsteady and steady flows, the solution is unique. The presence of nanoparticles affects the existence of dual solution in decelerating unsteady flow only when the medium of the flow is a porous medium. But different shapes of nanoparticles are not disturbing the dual solution existence range, though it has a considerable impact on thermal conductivity of the mixture. Different shapes of nanoparticles act differently to enhance the heat transfer characteristics of the base fluid, i.e., the water here. On the other hand, the existence range of dual solutions becomes wider for a larger permeability parameter related to the porous medium. Regarding the cooling rate of the heated surface, it rises with the permeability parameter, shape factor (related to various shapes of Cu-nanoparticles), and radiation parameter. The surface drag force becomes stronger with the permeability parameter. Also, with growing values of nanoparticle volume fraction, the boundary layer thickness (BLT) increases and the thermal BLT becomes thicker with larger values of shape factor. For decelerating unsteady flow, the nanofluid velocity rises with permeability parameter in the case of upper branch solution and an opposite trend for the lower branch is witnessed. The thermal BLT is thicker with radiation parameter. Due to the existence of dual solutions, a linear stability analysis is made and it is concluded that the upper branch and unique solutions are stable solutions.     

A model for strong attenuation and dispersion of seismic P-waves in a partially saturated fractured reservoir

In this work we interpret the data showing unusually strong velocity dispersion of P-waves (up to 30%) and attenuation in a relatively narrow frequency range. The cross-hole and VSP data were measured in a reservoir, which is in the porous zone of the Silurian Kankakee Limestone Formation formed by vertical fractures within a porous matrix saturated by oil, and gas patches. Such a medium exhibits significant attenuation due to wave-induced fluid flow across the interfaces between different types of inclusions (fractures, fluid patches) and background. Other models of intrinsic attenuation (in particular squirt flow models) cannot explain the amount of observed dispersion when using realistic rock properties. In order to interpret data in a satisfactory way we develop a superposition model for fractured porous rocks accounting also for the patchy saturation effect.     

含低渗透夹层的非均质双重介质垂向干扰的精确解

本文建立了单一介质、双重介质中由两个渗透层被一个致密低渗透层隔开的多层油藏渗流的数学模型,并求得了无穷大地层的精确解和长时渐近解。利用这个解可以在双重介质层状油藏的单井、多井试井中解释压力恢复曲线、垂向干扰试井和垂向脉冲试井。