复杂微通道气体流动的格子Boltzmann模拟

构建了一个模拟复杂微通道内气体流动的多松弛格子Boltzmann模型。该模型采用动力学曲面滑移边界,考虑了微尺度效应和努森层影响。此外,为了更准确地描述微通道内气体的滑移速度,在模型中引入孔隙局部Kn数来代替平均Kn数。之后采用Poiseuille流对模型进行验证,模拟结果与用直接模拟蒙特卡洛方法和分子模拟结果吻合较好,证明了该模型模拟微通道内处于滑移区和过渡区气体流动的有效性。最后,采用该模型模拟多孔介质内气体渗流过程。结果表明,随着孔隙平均Kn数的增加,多孔介质内的高渗区域增加,且优先从小孔隙中开始增加,这是由于小孔隙中微尺度效应更加明显,相对大孔隙流动阻力更小所致。     

HPAM溶液在石英微圆管中的流动实验

针对HPAM溶液在油藏孔隙中的复杂渗流特性问题,在内径为10～350μm微圆管中进行了部分水解聚丙烯酰胺(HPAM)溶液的流动实验.实验结果表明:在本文实验条件下,HAPM溶液在管径尺度较小的微圆管内的流动规律明显偏离常规尺度下的非牛顿流体力学流动规律,其截面平均流速高于按照常规尺度流动推算出的结果,且管径越小,偏离程度越大,微尺度效应越强.     

低渗透多孔介质渗流动边界模型的解析与数值解

考虑启动压力梯度的低渗透多孔介质非达西渗流模型属于强非线性动边界问题, 分别利用相似变量变换方法和基于空间坐标变换的有限差分方法, 对内边界变压力情况下、考虑启动压力梯度的一维低渗透多孔介质非达西渗流动边界模型进行了精确解析与数值求解研究. 研究结果表明:该动边界模型存在唯一的精确解析解, 且所求得的精确解析解可严格验证数值解的正确性;且当启动压力梯度值趋于零时, 非达西渗流动边界模型的精确解析解将退化为达西渗流情况下的精确解析解. 由求解结果作出的非零无因次启动压力梯度下的地层压力分布曲线表现出紧支性特点, 其与达西渗流模型的有显著不同. 因此, 研究低渗透多孔介质中非稳态渗流问题时, 应该考虑动边界的影响. 研究内容完善了低渗透多孔介质的非达西渗流力学理论, 为低渗透油气藏开发的试井解释与油藏数值模拟技术提供了理论基础.      

非极性小分子有机液体在微管道中的流量特性

用去离子水及有机液体在内径约为25μm的石英圆管内进行了流量特性实验．液体分子量范围为18~160,动力黏性系数的范围为0.5~1mPa·s.实验雷诺数范围为Re＜8．所用有机液体为:四氯化碳、乙基苯及环己烷都是非极性液体,其分子结构尺度小于1nm．实验结果表明,在定常层流条件下,圆管内的液体流量与两端压力差成正比,其压力－流量关系仍符合经典的Hagen－Poiseuille流动．这说明非极性小分子有机液体在本实验所用微米尺度管道中其流动规律仍符合连续介质假设．鉴于微尺度流动实验的特殊性,文中还介绍了微流动实验装置,分析了微尺度流动测量误差来源及提高测量精度的措施．     

基于渗透率连续变化的低渗透多孔介质非线性渗流模型

基于低渗透多孔介质渗透率的渐变理论,确定了能精确描述低渗透多孔介质渗流特征的非线性运动方程,并通过实验数据拟合．验证了非线性运动方程的有效性。非线性渗流速度关于压力梯度具有连续-阶导数,方便于工程计算;由此建立了低渗透多孔介质的单相非线性径向渗流数学模型,并巧妙采用高效的Douglas-Jones预估一校正有限差分方法求得了其数值解。数值结果分析表明：非线性渗流模型为介于拟线性渗流模型和达西渗流模型之间的一种中间模型或理想模型,非线性渗流模型和拟线性渗流模型均存在动边界;拟线性渗流高估了启动压力梯度的影响,使得动边界的移动速度比实际情况慢得多;非线性越强,地层压力下降的范围越小,地层压力梯度越陡峭,影响地层压力的敏感性减弱,而影响地层压力梯度的敏感性增强。     

孔隙介质中稠油流体非线性渗流方程

为揭示稠油流体在油藏孔隙中渗流特性,基于力学平衡方程,建立了描述稠油流体渗流特征的非线性渗流方程,对油藏孔隙中稠油渗流过程及启动机理进行了深入分析,着重分析了边界层、流体屈服应力以及表面力对渗流过程的影响.结果表明,Hagen-Poiseuille定律需经修正方能描述稠油流体流动,边界层、流体屈服应力以及表面力对稠油渗流影响非常显著.孔隙中,稠油启动压力梯度来源于其屈服应力、表面力,边界层加剧了渗流非线性程度,实际稠油油藏开发中,要充分掌握稠油渗流非线性特征.      

水力压裂中裂纹扩展的数值模拟

水力压裂是在高压粘滞流体或清水作用下地层内裂缝起裂与扩展的过程。由于包含岩石断裂和流-固耦合等复杂问题,对该过程的数值模拟具有相当大的挑战性。本文建立基于有限元与离散元混合方法的裂纹模型,模拟岩石裂纹扩展,实现了连续向非连续的转化;建立双重介质流动模型,裂隙流作为孔隙渗流的压力边界,孔隙渗流反作用裂隙的压力求解,处理了流体在基岩与人工裂缝中的协调流动;将裂纹模型与流体流动模式进行结合,建立断裂-应力-渗流耦合形式的力学模型,进一步分析了水力压裂的基本过程,综合多种数值计算方法,编写程序,在验证岩体裂纹模型与双重介质流动模型有效性的基础上,对压裂过程进行复现,将模拟结果与文献结果进行了对比,并讨论了所构建模型的优缺点。     

破碎岩石气体渗透性的试验测定方法

为了测试采矿工程中松散破碎岩体的气体渗透特性,设计了与MTS815.02岩石力学伺服机配套的气体渗透仪和测试系统。依靠MTS伺服机改变压力和位移条件,渗透仪用以稳定破碎岩石并连接气体测试系统,接入气源、组成回路和测试气透特性等由测试系统完成。通过对几组破碎砂岩岩样氮气渗透性能的测试,结合考虑状态方程的渗透率换算公式,得出破碎砂岩气体渗透率随压力和粒径条件的变化规律。该测试方法和测试结果在研究采动破碎岩体瓦斯等气体流动规律方面具有重要意义。     

渗流力学的近况和展望

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

基于规则网络的碳酸盐岩多尺度网络模型构建方法研究

碳酸盐岩油藏非均质性强,孔隙大小变化可达好几个数量级,描述碳酸盐岩油藏多尺度孔隙特征具有重要意义.本文首先基于三维规则网络模型建立了不同物理尺寸的溶洞网络、大孔隙网络和微孔隙网络;然后提出一种耦合算法,以溶洞网络为基础,通过添加适当比例的大孔隙和微孔隙,构建出碳酸盐岩多尺度网络模型;最后对比分析了各网络模型的几何性质、拓扑性质和绝对渗透率.结果表明,碳酸盐岩多尺度网络模型能够同时描述不同尺度孔隙的几何和拓扑特征;且相比各单一尺度的孔隙网络模型,多尺度网络模型有着较高的绝对渗透率,这是由于各尺度孔隙之间的相互连通极大地提高了网络的整体连通性和流动能力,为碳酸盐岩油藏微观渗流模拟提供了重要的研究平台.     

An overview on nonlinear porous flow in low permeability porous media

This paper gives an overview on nonlinear porous flow in low permeability porous media, reveals the microscopic mechanisms of flows, and clarifies properties of porous flow fluids. It shows that, deviating from Darcy's linear law, the porous flow characteristics obey a nonlinear law in a low-permeability porous medium, and the viscosity of the porous flow fluid and the permeability values of water and oil are not constants. Based on these characters, a new porous flow model, which can better describe low permeability reservoir, is established. This model can describe various patterns of porous flow, as Darcy's linear law does. All the parameters involved in the model, having definite physical meanings, can be obtained directly from the experiments.     

气液混输悬链线立管系统两相流特性实验

严重段塞流是海洋工程气液混输管线--立管系统中常见的一种特殊有害流动现象, 采用水平--下倾--悬链线立管气液混输组合管道系统, 通过系列实验在悬链线立管中获得了严重段塞流、间歇流和震荡流等流型, 阐述了这些流动现象的形成机理, 提出了能够产生严重段塞流的判定准则. 结果表明, 悬链线立管严重段塞流具有明显周期性, 在一个周期内的流动特征可分为液塞形成、液体出流、液气喷发及液体回流等4个阶段, 进而给出了各阶段中相关流动参数的变化规律. 在实验中同时还对悬链线与垂直立管中严重段塞流形成机理进行了比较分析, 发现两者在液塞形成阶段有显著差别. 其中, 在悬链线立管中液塞形成之前首先需要经历一个气液混合液塞形成过程, 而垂直立管则没有这个过程.      

微尺度气体流动

了解微尺度气体流动特点是微机电系统设计和优化的基础.有关的研究可以上溯到20世纪初Knudsen的平面槽道流动质量流量的测量和Millikan的小球阻力系数的测量,实验结果揭示了稀薄气体效应即尺度效应对气体运动的重要影响.由于流动特征长度很小,微尺度气流经常处于滑流区甚至过渡领域,流动的相似参数为Knudsen数和Mach数.因此可以考虑利用相似准则,通过增大几何尺寸、减小压力的途径,解决微机电系统实验观测遇到的困难.为解决直接模拟Monte Carlo方法分析微机电系统中低速稀薄气流遇到的统计涨落困难,我们提出了信息保存法(IP),该方法能够有效克服统计散布,并已成功用于多种微尺度气流.      

气体运动论数值算法在微槽道流中的应用研究

简要介绍基于Boltzmann模型方程的气体运动论数值算法基本思想及其对二维微槽道流动问题数值计算的推广，并阐述适用于微尺度流动问题的气体运动论边界条件数值处理方法。通过对压力驱动的二维微槽道流动问题进行数值模拟，将不同Knudsen数下的微槽道流计算结果分别与有关DSMC模拟值和经滑移流理论修正的N—S方程解进行比较分析，表明基于Boltzmann模型方程的气体运动论数值算法对微槽道气体流动问题具有很好的模拟能力。     

Experimental validation of a two equation RANS transitional turbulence model for compressible microflows

Laminar-to-turbulent flow transition in microchannels can be useful to enhance mixing and heat transfer in microsystems. Typically, the small characteristic dimensions of these devices hinder in attaining higher Reynolds numbers to limit the total pressure drop. This is true especially in the presence of a liquid as a working medium. On the contrary, due to lower density, Reynolds number larger than 2000 can be easily reached for gas microflows with an acceptable pressure drop. Since microchannels are used as elementary building blocks of micro heat exchangers and micro heat-sinks, it is essential to predict under which conditions, the laminar-to-turbulent flow transition inside such geometries can be expected. In this paper, experimental validation of a two equations transitional turbulence model, capable of predicting the laminar-to-turbulent flow transition for internal flows as proposed by Abraham etal. (2008), is presented for the first time for microchannels. This is done by employing microchannels in which Nitrogen gas is used as a working fluid. Two different cross-sections namely circular and rectangular are utilized for numerical and experimental investigations. The inlet mass flow rate of the gas is varied to cover all the flow regimes from laminar to fully turbulent flow. Pressure loss experiments are performed for both cross-sectional geometries and friction factor results from experiments and numerical simulations are compared. From the analysis of the friction factor as a function of the Reynolds number, the critical value of the Reynolds number linked to the laminar-to-turbulent transition has been determined. The experimental and numerical critical Reynolds number for all the tested microchannels showed a maximum deviation of less than 12%. These results demonstrate that the transitional turbulence model proposed by Abraham etal. (2008) for internal flows can be extended to microchannels and proficiently employed for the design of micro heat exchangers in presence of gas flows.     

Planar dynamics of inclined curved flexible riser carrying slug liquid–gas flows

Flexible risers transporting hydrocarbon liquid–gas flows may be subject to internal dynamic fluctuations of multiphase densities, velocities and pressure changes. Previous studies have mostly focused on single-phase flows in oscillating pipes or multiphase flows in static pipes whereas understanding of multiphase flow effects on oscillating pipes with variable curvatures is still lacking. The present study aims to numerically investigate fundamental planar dynamics of a long flexible catenary riser carrying slug liquid–gas flows and to analyse the mechanical effects of slug flow characteristics including the slug unit length, translational velocity and fluctuation frequencies leading to resonances. A two-dimensional continuum model, describing the coupled horizontal and vertical motions of an inclined flexible/extensible curved riser subject to the space–time varying fluid weights, flow centrifugal momenta and Coriolis effects, is presented. Steady slug flows are considered and modelled by accounting for the mass–momentum balances of liquid–gas phases within an idealized slug unit cell comprising the slug liquid (containing small gas bubbles) and elongated gas bubble (interfacing with the liquid film) parts. A nonlinear hydrodynamic film profile is described, depending on the pipe diameter, inclination, liquid–gas phase properties, superficial velocities and empirical correlations. These enable the approximation of phase fractions, local velocities and pressure variations which are employed as the time-varying, distributed parameters leading to the slug flow-induced vibration (SIV) of catenary riser. Several key SIV features are numerically investigated, highlighting the slug flow-induced transient drifts due to the travelling masses, amplified mean displacements due to the combined slug weights and flow momenta, extensibility or tension changes due to a reconfiguration of pipe equilibrium, oscillation amplitudes and resonant frequencies. Single- and multi-modal patterns of riser dynamic profiles are determined, enabling the evaluation of associated bending/axial stresses. Parametric studies reveal the individual effect of the slug unit length and the translational velocity on SIV response regardless of the slug characteristic frequency being a function of these two parameters. This key observation is practically useful for the identification of critical maximum response.     

Parameter Determination Using 3D FIB-SEM Images for Development of Effective Model of Shale Gas Flow in Nanoscale Pore Clusters

A large amount of nano-pores exists in pore clusters in shale gas reservoirs. In addition to the multiple transport regimes that occur on the nanoscale, the pore space is another major factor that significantly affects the shale gas recoverability. An investigation of the pore-scale shale gas flow is therefore important, and the results can be used to develop an effective cluster-scale pore network model for the convenient examination of the process efficiency. Focused ion beam scanning electron microscope imaging, which enables the acquisition of nanometre-resolution images that facilitate nano-pore identification, was used in conjunction with a high-precision pore network extraction algorithm to generate the equivalent pore network for the simulation of Darcy and shale gas flows through the pores. The characteristic parameters of the pores and the gas transport features were determined and analysed to obtain a deeper understanding of shale gas flow through nanoscale pore clusters, such as the importance of the throat flux–radius distribution and the variation of the tortuosity with pressure. The best parameter scheme for the proposed effective model of shale gas flow was selected out of three derived schemes based on the pore-scale prediction results. The model is applicable to pore-scale to cluster-scale shale gas flows and can be used to avoid the multiple-solution problems in the study of gas flows. It affords a foundation for further study to develop models for shale gas flows on larger scales.     

Meso-level simulation of gas hydrate dissociation inlow-permeability sediments

This paper presents a meso-level simulation of gas hydrate dissociation in low-permeability marine sediments. Interstitial pores are defined to describe fluid flow and particle movement. The proposed model couples multiphase fluid flow with particle movement to simulate the thermodynamics of gas hydrate dissociation triggered by sharp temperature rises. Hydrates respond quickly to temperature rise in low-permeability sediments. Dissociation causes pore pressure to rise rapidly to equilibrium then steadily increase above equilibrium pressure. Lower permeability sediment builds up greater excess pore pressure as the dissipation of pore pressure is constrained.