水平井及大斜度井砾石充填α波砂床平衡高度随机概率模型及实验研究

α波砂床平衡高度预测是水平井砾石充填过程数值模拟及施工参数优化设计的核心问题之一。根据水流紊动理论,将砂床表面流体瞬时流速视为正态分布的随机变量,将颗粒运动状态转换视为随机事件,分析状态转换事件概率及最终状态概率的数字特征。根据最终状态概率代表颗粒运动趋势的原理,砂床高度的变化是床层表面颗粒运动状态转换的结果。由平衡砂床高度和摩阻流速及平均流速的关系,建立了同时适用于水平井和大斜度井砾石充填的α波砂床平衡高度预测的随机概率模型。利用研制的大斜度井砾石充填全尺寸实验模拟装置,进行井筒倾角、筛管偏置度、流量、砂浓度对α波砂床平衡高度影响的实验研究,对随机概率模型的计算结果进行了对比验证。随机概率模型与实验结果比较吻合,平均误差10.5%,进行必要的修正可以用于大斜度井砾石充填α波砂床平衡高度预测。     

防砂井挡砂介质堵塞规律实验及堵塞程度定量预测模型

防砂井挡砂介质堵塞造成降产已成为困扰油气井防砂的关键问题。已开展的大量堵塞规律研究主要涉及机理和定性描述,尚未形成可用的堵塞程度定量预测方法。针对这一问题,本文在广泛应用的砾石充填防砂条件下,利用挡砂介质堵塞评价驱替实验装置,进行了砾石层堵塞机理和规律实验,着重考察了砾石层堵塞渗透率比随驱替时间、砾砂中值比、流体流速、流体粘度、粘土含量、细质含量等生产条件的定量变化规律。根据实验数据,首先拟合砾石充填层堵塞渗透率比随时间的定量关系,引入综合堵塞系数和单项堵塞系数表征各生产条件对堵塞规律的影响,然后拟合了堵塞系数与各生产参数的经验关系,整合得到一个根据防砂参数和生产条件预测砾石充填层堵塞程度随时间变化规律的定量模型。结果表明,不同生产条件砾石层渗透率在投产早期下降较快,然后下降速度减缓并趋于基本不再变化,达到堵塞平衡状态。具体变化规律由防砂参数和生产条件控制。使用实际防砂井资料进行了模型应用案例分析,结果和规律与现场实际基本相符。该模型进一步修正后可用于防砂井动态产能预测与评价。     

石油防砂割缝筛管的冲蚀磨损性能研究

以直径为0.1～0.6 mm的石英砂为磨料,黏度为10 mPa·s的原油为携砂液,在自制的冲蚀磨损试验装置上研究了含砂原油对石油防砂割缝筛管的冲蚀磨损性能,采用LW200-4CS数码倒置金相显微镜观测了筛管缝口表面的冲蚀磨损形貌;利用曲线拟合方法建立了割缝筛管冲蚀磨损率与原油流速、砂粒直径、砂粒浓度、冲蚀角之间关系的数学模型.结果表明:筛管的冲蚀磨损率随含砂原油流速的增加而成指数关系增加;开始时随着原油中所携带砂粒直径和浓度的增加,筛管的冲蚀磨损率增加,达到一定数值后砂粒直径和浓度再增加,割缝筛管的冲蚀磨损率反而降低;随着冲蚀角的增大,割缝筛管的冲蚀磨损率增加,在冲蚀角为90°时达到最大,这种现象不同于塑性材料的冲蚀磨损机理;所建立的数学模型与试验结果较吻合.     

注塑成型充填／后充填过程统一的可压缩流动分析

注塑成型是重要的塑料成型工艺,成型过程中熔体在模腔中的流动和传热对最终制品的性能和质量有重要的影响,因此,精确预测注塑过程的流动及传热历史,并进一步预测注塑制品的收缩、翘曲和机械性能等性能和质量指标具有重要意义。为了精确地描述成型过程中材料的流动及传热行为,本文针对注塑成型过程的工艺特点,将充填后充填过程作为一个统一的过程,考虑材料可压缩性及相变对充填和后充填过程的影响,建立了充填后充填过程的统一数学模型。采用有限元/有限差分/控制体积混合数值方法,实现了注塑成型充填后充填一体化模拟。数值模拟结果与实验结果的对比,验证了本文模型和算法。     

偏心环空层流螺旋流速度的近似解析解

本文分析计算了石油钻井工程偏心环空中内管转动的螺旋流动。将偏心环空的三维流动简化为环空薄层内的二维流动，从黎曼空间曲面张量表达的动量方程出发，求得了偏心环空层流螺旋流动速度的近似解析解，用此方程的同心环空表达式与精确解比较，精度较高。     

????????????????????????????о?

基于射孔完井水平裂缝中支撑剂均匀充填，考虑充填层强度、油层和流体 物性以及裂缝的渗流条件，研究了地应力及渗流引起的附加应力联合作用的支撑剂充填层渗 流模型和力学模型，并进行现场应用. 研究表明，所建立的模型对压裂油井支撑剂回采出砂 状况分析具有较高的符合率，油井总体符合率、不出砂井符合率和出砂井符合率分 别为84.66\%, 84.18\%和88.10\%, 为大庆油田压裂油井防砂措施实施方案制定提 供了决策支持，有利于提高措施有效率.     

黏性不可压缩流体流动前沿的数值模拟

提出了模拟注射成型中黏性、不可压缩流体流动前沿的新方法. 将Hele-Shaw流动应用于非 等温条件下的黏性、不可压缩流体，建立了流动分析模型，用充填因子的输运方程描述流动 前沿. 应用高阶Taylor展开式计算每一时间步长的充填因子，用Galerkin方法导出了计算 充填因子各阶导数的递推公式. 给出了时间增量的选取方法，证明了它的稳定性. 针对Han 设计的试验模具，用相同的材料及工艺条件模拟充填过程，比较了传统方法和该方法的模 拟结果与实验结果的差异. 算例分析表明，该方法可以有效地提高注射成型中流动前沿的 模拟精度和计算效率.     

环空后台阶管道内气固两相流动的数值模拟研究

环空管道后台阶突扩流动是空气正循环钻井过程中十分重要的关键部分,直接决定了钻探岩屑是否能够顺利上返地面.本文采用Euler-Lagrange两相流研究方法,气相湍流采用Realizable k-ε模型,固相采用离散相模型(DPM),固相的湍流耗散采用随机轨道模型,对环空后台阶突扩管道内气固两相流动进行了数值模拟研究.得出了气相场大涡演变规律,在此基础上研究了不同粒径时颗粒在流道中的浓度分布规律、运动轨迹,以及速度场分布规律.这为细观框架下研究气固两相相互作用规律提供了依据.     

用代数应力模型模化型腔充填过程中的紊流流动

针对型腔充填过程中的紊流流动，用代数应力模型研究紊流流动现象．成功地将代数应力模型引入ＰＨＯＥＮＩＣＳ软件中，完成了恒温流体充填过程三维时均速度和自由表面的计算模拟．与实验结果相比说明，本文基于ＰＨＯＥＮＩＣＳ软件开发的充填模拟具备了对复杂型腔充填过程的数值计算．     

注塑成型充模过程的温度场计算

对注塑成型过程中充填阶段的非等温效应的数值模拟进行了研究 ,模拟了非牛顿流体在任意形状型腔内的流动。这里 ,数学模型是基于非牛顿流体在非等温状态下的广义 Hele-Shaw流动 ,利用有限元 /有限差分混合数值方法求解流场中的压力和温度分布     

Experimental Conditions Favoring the Formation of Fluid Banks during Counter-Current Flow in Porous Media

The purpose of this study is to investigate factors that affect the formation of fluid banks during gravity-driven counter-current flow in porous media. To our knowledge, development of a fluid bank has been observed in only one previous counter-current flow experiment, although there are some hints of fluid banks in other experiments. We have undertaken experimental and simulation studies to confirm the presence of such banks and to delineate factors which enhance or inhibit their formation. Experiments were performed using glass bead packs and X-ray Computed Tomography to monitor saturation distribution as a function of time. The simulation approach considers saturation history at every point in the sample, defining conditions at each time point from hysteresis in capillary pressure and relative permeability. The model proved to reproduce experimental observations accurately. The experiments and associated model show that a minimal vertical sample height is needed for the development of a fluid bank. In addition, round sample boundaries and higher average nonwetting phase saturation tend to prevent the formation of a bank. The validated model can improve our ability to predict and optimize counter-current flow processes, both in the laboratory and in the field (e.g. exploration and hydrocarbon extraction).     

非牛顿流体在非均质油藏渗流压力场实验

在非均质油藏模型上进行非牛顿流体流动物理模拟实验，对比研究水驱、聚合物驱和交联聚合物对提高石油采收率的影响．通过布置高精度的压差传感器测量不同驱替过程模型中的渗流压力场的动态变化，成胶后的交联聚合物封堵了高渗条区，改变了油藏内流体流动方向，驱替出低渗区内油，提高了采收率．     

Studies on the law of laminar helical flow of power-law fluid in eccentric annuli

According to the principle of fluid mechanics,the law of laminar,helical flow of power-law fluid in eccentric annuli is studied extensively in this paper.The apparent viscosity,velocities distribution of laminar helical flow of power-law fluid are discussed and calculating methods of flow rate and pressure loss are presented.The factors influencing flow are also analysed.On the basis of theoretical studies some new results of the present paper are compared with the theories of the helical flow of the power-law fluid in concentric annuli.The test verified that theoretical formulas in this article are proper and general.They can be used for calculating hydraulic parameters in drilling engineering.     

An experimental study of the effect of wall temperature nonuniformity on natural convection in an enclosure heated from the side

An experimental study of natural convection in a parallelepipedal enclosure induced by a single vertical wall is described. The upper half of this wall was warm and the lower half cold. The other enclosure walls were insulated. The temperature and flow measurements were performed in the high Rayleigh number regime (10¹⁰<Ra<5×10¹⁰) by using water as the working fluid. The Rayleigh number was based on the enclosure height and the temperature difference between the warm and the cold part of the driving wall. The flow field featured two flat cells, one filled with warm fluid along the top horizontal wall, and the other filled with cold fluid along the bottom horizontal wall. Each of these cells was surrounded by an additional cell as tall as half the enclosure height. The above flow structure prohibited extensive thermal contact between warm and cold fluid, thus limiting the role of convection on the heat transfer process in the cavity. The findings of this study differ significantly from the findings of previous studies based on the ‘classical’ enclosure model possessing two isothermal vertical walls, the one warm and the other cold, and support the view that the use of ‘more realistic‘ temperature boundary conditions in enclosure natural convection needs careful examination.     

Dynamics of a long tubular cantilever conveying fluid downwards,which then flows upwards around the cantilever as a confined annular flow

A theoretical model is developed for the dynamics of a hanging tubular cantilever conveying fluid downwards; the fluid, after exiting from the free end, is pushed upwards in the outer annular region contained by the cantilever and a rigid cylindrical channel. This configuration thus resembles that of a drill-string with a floating fluid-powered drill-bit. The linear equation of motion is solved by means of a hybrid Galerkin–Fourier method, as well as by a conventional Galerkin method. Calculations are conducted for a very slender system with parameters appropriate for a drill-string, for different degrees of confinement of the outer annular channel; and also for another, bench-top-size experiment. For wide annuli, the dynamics is dominated by the internal flow and, for low flow velocities, the flow increases the damping associated with the presence of the annular fluid. For narrow annuli, however, the annular flow is dominant, tending to destabilize the system, giving rise to flutter at remarkably low flow velocities. The mechanisms underlying the dynamics are also considered, in terms of energy transfer from the fluid to the cantilever and vice versa, as are possible applications of this work.     

Experimental investigation on flow characteristics in a narrow annulus

Experiments are carried out to investigate the flow characteristics with/without heat exchange in a narrow annulus. In the experiments, directions of flow include horizontal, upstream and downstream flow. Experimental results show that the flow characteristics of water through the narrow annulus are different from those in normal tubes. Flow directions have little influences on the flow friction for the fluid flow in the narrow annulus with/without heat exchange. The flow characteristics in the narrow annulus have relations to the liquid temperature difference at the inlet and outlet of the annulus. Their influences on the flow characteristics are relatively obvious in the laminar flow area. When the Reynolds number is larger than 10⁴, there are little differences between the flow friction factors with/without heat exchange. It is also found that the asymmetrical flow can make the friction factor increase, whereas the symmetrical flow can reduce the flow friction. In the experiments, the transition from laminar to turbulent flow is carefully observed. In the narrow annulus, the flow transition is initiated earlier than that in normal pipes at a Reynolds number range from 1,100 to 1,500, which is different from the heat transfer transition. The results are gained to provide bases for the further investigations on the two-phase flow in narrow annuli.     

MHD liquid metal flow and heat transfer between vertical coaxial cylinders under horizontal magnetic field

Direct numerical simulations are presented of MHD liquid metal flow and heat transfer in vertical annuli. Three annular gaps and four ratios of annular height to annular gap are considered. The walls of the external and internal cylinders are isothermal with the temperature of the outer cylinder being higher and, thus, buoyancy is the driving force. The results show that the fluid motion increases as the aspect ratio and the annular gap become larger. The presence of the magnetic field results to fluid deceleration and, thus, to flow stabilization. Additionally, non symmetric flow patterns develop, due to the magnetic field, resulting in differently sized normal and parallel wall layers, namely the Hartmann and the Roberts layers, respectively. For all annular gaps considered, the highest spatially averaged heat transfer rates are obtained for aspect ratios equal to 1.     

Developing laminar forced convection in eccentric annuli

The paper presents a boundary-layer model describing the laminar forced convection heat transfer in the entry region of eccentric annuli. A finite-difference numerical algorithm is developed for solving this model. Numerical results are presented for the developing velocity profiles and the pressure drop in annuli of radius ratio 0.5 and 0.9 over a dimensionless eccentricity ranged from 0.1 to 0.8. Heat transfer parameters are presented for a fluid of Pr=0.7 under the conditions of an isothermally heated inner wall while the outer wall is kept at the inlet fluid temperature. Received on 18 March 1997     

Numerical investigation of fully developed laminar flow in irregular annuli: Square‐circular ducts

Numerical study is carried out on fully developed laminar axial flow of incompressible Newtonian fluid in irregular annuli of different cross‐sectional geometry. Accurate prediction of the friction factor of this flow was obtained using high‐order finite element method. Investigations are performed for the annular flow between rectangular and circular ducts under all possible arrangements. The governing equations are solved by an efficient finite element technique. It was found that a higher annular area‐ratio will lead to a monotonic increase in (f Re) value in the case of regular annuli, and will lead to an increase followed by a decrease in (f Re) value in the case of irregular annuli. Furthermore, this study illustrates that irregular annuli have lower (f Re) value than regular annuli, and that the square‐in‐circle case (SC‐case) has the lowest (f Re) value whereas the circle‐in‐circle case (CC‐case) has the highest (f Re) value. The finite element technique developed in this study was validated against the already reported numerical results in the literature, and the results presented in this work show excellent agreement. Copyright © 2009 John Wiley & Sons, Ltd.