分形复合油藏非牛顿幂律流体不稳定渗流的数学模型

对不稳定渗流的数学模型进行了推导并建立了分形复合油藏不稳定渗流模型.在无限大地层、有界定压、有界封闭三种外边界条件下分别求出了它们在Laplace空间的解析解.对于两区的特殊情况,分析了井底压力动态特征和参数影响,制作了典型曲线.非牛顿幂律流体的幂律指数,分形参数均对典型曲线产生较大的影响,呈现出与牛顿流体和均质油藏明显不同的特征.这对于非均质油藏非牛顿流体的不稳定试井分析和研究非线性渗流特征都是十分重要的.     

分形低渗透缝洞型油藏非线性渗流规律

引入分形理论,建立考虑低速非达西效应的分形三重介质缝洞型油藏数学模型,通过Laplace变换及Stehfest数值反演方法求出井底压力,借助Matlab编程绘制压力动态曲线,划分渗流阶段,分析渗流规律,进行非线性参数敏感性分析.最后结合实际算例,验证模型的正确性.结果表明：分形三重介质油藏渗流过程分为早期纯井储,过渡流,缝洞窜流,拟径向流,基质与溶洞、裂缝窜流及总体径向流6个渗流阶段;低速非达西效应对渗流的影响随时间的推移逐渐增大;启动压力梯度越大,总径向流阶段压力动态曲线上翘幅度越大;分形系数影响整个渗流过程,随着分形系数的增大,裂缝迂曲程度随之增大,致使渗流阻力增加,引起压力动态曲线整体上移幅度增大.     

考虑渗透率张量的非均质油藏有限元数值模拟方法

针对具有混合边界的非均质油藏,考虑全张量形式的渗透率,建立弹性微可压缩单相流体不稳定渗流问题的数学模型.根据变分原理,将压力微分方程的边值问题转化为泛函的极值问题,建立渗流模型的有限元方程.针对典型的均质和非均质渗流问题进行模拟计算,得到油藏内压力动态分布曲线,并分析曲线特征.研究表明,有限元法计算精度很高,适用于求解利用渗透率张量表征的非均质油藏渗流问题.为非均质油藏的开发和精细油藏数值模拟提供了理论依据.     

考虑动半径和粘度变化的稠油非牛顿不稳定渗流数学模型

针对稠油非牛顿特征,在Bingham流体渗流方程基础上,通过对动半径和粘度进行表征,建立同时考虑启动压力梯度、动半径变化和粘度变化的非牛顿稠油不稳定渗流数学模型,完善Bingham型稠油渗流数学模型.通过空间、时间离散差分及Matlab数值计算,得到非牛顿稠油非稳态渗流地层压力分布.结果表明,相同产量下,随启动压力梯度增大,动半径向井方向移动;启动压力梯度越大,压降曲线越陡,相应近井压降越大;相同启动压力梯度下,产量越大,不同吞吐周期压力差距越大.将半径和粘度动态变化相结合,弥补了现行非牛顿稠油渗流数学模型的-个缺陷.     

使用混合网格计算非达西渗流

针对垂直裂缝井的特殊流动模式,从非达西定律出发,建立二维平面的非达西渗流方程.通过建立一组无量纲量,最终得到无量纲的渗流方程及其定解条件.假定外边界为圆形,用PEBI网格及混合网格对求解区域进行网格划分,用有限差分法对无量纲的方程进行离散,最终得到垂直裂缝井的井底压力数值解.根据此数值解并考虑井筒存储和表皮因子的影响,得到真实垂直裂缝井的井底压力.对计算结果的分析表明,使用混合网格求解非达西渗流井底压力相当准确,该方法也适用于水平井等更复杂井型及复杂边界的问题求解.     

侧钻井压力分析的边界元法及储层伤害的评价

随着侧钻井的增加,非均质油藏侧钻井的压力动态特征和对储层伤害进行评价已变得越来越重要了。为了研究在非均质油藏中侧钻井的压力动态,建立了三维非均质油藏侧钻井的不稳定渗流数学模型。用边界元方法求得模型的解并分析了压力动态特征。通过拟表皮系数、无因次污染半径、无因次井筒储积系数和污染区渗透率等参数的合理组合,绘制了能用于侧钻井储层伤害评价的典型曲线。并对辽河油田某侧钻井的压力恢复试井数据进行了拟合。     

考虑相重新分布的双重介质分形油藏不稳定渗流压力动态特征

基于分形油藏渗流力学,针对存在井筒储集、表皮效应和井筒相重新分布影响下的不稳定渗流问题,建立双重介质分形油藏有效井径数学模型,并采用Douglas-Jones预估校正法求得无限大地层定产量生产条件下的非线性数值解.由数值解可知,该系统完整的压降曲线由四个流动区组成.最后分析各个分形参数、相重新分布参数和双重介质参数变化时压力的变化规律,做出典型压力曲线图版,为研究复杂渗流系统中相重新分布影响下的不稳定渗流规律提供参考.     

特低渗透油藏面积井网见水时间计算

特低渗透油藏在水驱开发过程中常表现出渗透率各向异性和流体非达西渗流特征,运用非达西渗流公式及流线积分法,得到特低渗透油藏五点井网、反九点井网和菱形反九点井网在油水两相非活塞式驱替条件下的油井见水时间.计算长庆鄂尔多斯盆地某特低渗透油藏井网见水时间,结果与实际动态符合较好.分析油水黏度比、渗透率各向异性及启动压力梯度对油井见水时间的影响,结果表明：油水黏度比越大,油井的见水时间越早;当菱形反九点井网长轴方向井距与短轴方向井距之比与渗透率各向异性强度大致相等时,长轴方向与短轴方向上的角井能实现均衡驱替;启动压力梯度延缓了油井的见水时间,且生产井距越大,启动压力梯度的影响越显著.     

低渗透多孔介质非达西渗流动边界界面追踪

基于低渗透多孔介质非达西不稳定渗流的动边界数学模型,推导动边界移动速度的微分表达式,揭示动边界移动速度与动边界上地层压力关于径向距离的二次导数成正比;由此利用拉格朗日三点插值公式求得动边界附近控制方程的有限差分格式,并对下一时刻动边界的精确位置进行追踪.有限差分方法的数值结果表明界面追踪法可较好地反映低渗透多孔介质非达西不稳定渗流动边界的移动规律.     

部分射开压裂直井渗流压力动态分析

建立符合非均质油藏部分射开地层实际的物理模型和三维各向异性矩形油藏的不稳定渗流数学模型,考虑不渗透顶、底边界和定压顶、底边界等边界条件的组合,通过无因次量纲变换、拉普拉斯变换、傅里叶余弦变换和分离变量等方法,得到拉普拉斯域的解析解,利用斯蒂芬森数值反演方法,得出实数域的压力数值解.绘制压力动态曲线,并进行敏感性分析.计算结果与数值模拟基本吻合,证实方法的可靠性.敏感性分析表明:压力动态曲线可分为早期线性流、中期径向流、晚期球形流、边界控制流四个流动期.裂缝长度主要影响早期线性流,渗透率各向异性主要影响中期径向流,储层射开程度和裂缝方位主要影响晚期球形流,边界条件和油藏宽度主要影响边界控制流.该方法可以确定最优射开程度、垂向渗透率等参数,为油藏工程分析和压裂工艺设计提供指导.     

Analysis of the flow of non-Newtonian visco-elastic fluids in fractal reservoir with the fractional derivative

In both the oil reservoir engineering and seepage flow mechanics, heavy oil with relaxation property shows non-Newtonian rheological characteristics. The relationship between shear rate g& and shear stress t is nonlinear. Because of the relaxation phenomena of heavy oil flow in porous media, the equation of motion can be written as[1] 2,rrvpqkppqtrrtll秏骣+=-+琪抖桫 (1) where lv and lp are velocity relaxation and pressure retardation times. For most porous media, the above motion equation (1)…     

Tracer dispersion in power law fluids flow through porous media: Evidence of a cross-over from a logarithmic to a power law behaviour

An analytical model is presented to describe the dispersion of tracers in a power-law fluid flowing through a statistically homogeneous and isotropic porous medium. The model is an extension of Saffman's approach to the case of non-Newtonian fluids. It is shown that an effective viscosity depending on the pressure gradient and on the characteristics of the fluid, must be introduced to satisfy Darcy's law. An analytical expression of the longitudinal dispersivity is given as a function of the Peclet number Pe and of the power-law index n that characterizes the dependence of the viscosity on the shear rate . As the flow velocity increases the dispersivity obeys an asymptotic power law: . This asymptotic regime is achieved at moderate Peclet numbers with strongly non-Newtonian fluids and on the contrary at very large values when n goes to 1 ( for n=0.9). This reflects the cross-over from a scaling behaviour for towards a logarithmic behaviour predicted for Newtonian fluids (n=1). Received: 22 July 1997 / Revised and Accepted: 2 July 1998     

Fractal Analysis of Power-Law Fluid in a Single Capillary

The fractai expressions for flow rate and hydraulic conductivity for power-law fluids in a single capillary are derived based on the fractai nature of tortuous capillaries. Every parameter in the proposed expressions has clear physical meaning. The flow rate and hydraulic conductivity for power-law fluids are found to be related to the tortuosity fractal dimension and the power-law index. The flow rate for power-law fluids increases with the increasing power-law index but decreases with the increasing tortuosity fractal dimension. Good agreement between the model predictions for flow in a fractai capillary and in a converging-diverging duct is obtained. The results suggest that the fractal capillary model can be used to model the power-law fluids with different rheologicai properties.     

稠油油藏试井理论研究

分析研究国内外稠油热采模型文献资料,总结归纳了9种稠油热采模型,主要包括非牛顿流体均质地层模型、二区复合模型、三区复合模型、多区复合模型、考虑重力超覆的数学模型以及干扰分析等,并对部分压力变化曲线及压力导数变化曲线做了简要的分析.比较各模型的优缺点,并指出目前的模型存在的问题以及稠油试井的发展趋势.     

Marangoni flow and mass transfer of power-law non-Newtonian fluids over a disk with suction and injection

We scrutinize the approximate analytical solutions by the optimal homotopy analysis method (OHAM) for the flow and mass transfer within the Marangoni boundary layer of power-law fluids over a disk with suction and injection in the present paper. Concentration distribution on the surface of a disk varies in a power-law form. The non-Newtonian fluid flow is due to the surface concentration gradient without considering gravity and buoyancy. According to the conservation of mass, momentum and concentration, the governing partial differential equations are established, and the appropriate generalized Kármán transformation is found to reduce them to the nonlinear ordinary differential equations. OHAM is used to access the approximate analytical solution. The influences of Marangoni the number, suction/injection parameters and power-law exponent on the flow and mass transfer are examined.     

Mechanical Properties and Flow Behavior of Polymers for Enhanced Oil Recovery

The mechanical properties and flow behavior in porous media of three different polymer systems including a hydrophobically modified acrylamide-based copolymer (HMSPAM), a partially hydrolyzed polyacrylamide (HPAM), and a polysaccharide (xanthan gum) were evaluated to establish their functional differentiation as mobility control agents in enhanced oil recovery (EOR). The rheological properties of the polymers were described by the power-law model to investigate their non-Newtonian behavior. The first normal stress difference (N₁) and Weissenberg number (W_e) were also used to compare their elastic properties. The experimental results showed that, at comparable shear viscosity, HMSPAM exhibited significant elasticity compared to HPAM and xanthan gum. Shear resistance tests indicated that all of the polymers experienced an extra stress when converging into a capillary tube due to the “entrance effect.” Xanthan gum was the most mechanically stable polymer. Moreover, HMSPAM showed the superior reformability which was quantified by the regained viscosity with relaxation time. This could be explained by the rapid re-association of the hydrophobic interactions. Sandpack flood tests indicated that HMSPAM rendered extremely high mobility control ability during polymer flooding suggesting its potential in EOR. However, this polymer also experienced significant retention within the porous media (potential injectivity and plugging problems), which may be attributed to the formation of bulky associative polymer networks. In this work, UV spectrometry was employed to monitor the concentration of the produced polymer solutions and quantify the polymer retention within porous media. This analytical approach offers great reliability and simplicity. It was concluded that the use of a particular polymer system depends on the oil reservoir conditions and the target EOR application.     

The inertial dynamics of thin film flow of non-Newtonian fluids

Consider the flow of a thin layer of non-Newtonian fluid over a solid surface. I model the case where the viscosity depends nonlinearly on the shear-rate; power law fluids are an important example, but the analysis here is for general nonlinear dependence. The modelling allows for large changes in film thickness provided the changes occur over a relatively large enough lateral length scale. Modifying the surface boundary condition for tangential stress forms an accessible foundation for the analysis where flow with constant shear is a neutral critical mode, in addition to a mode representing conservation of fluid. Perturbatively removing the modification then constructs a model for the coupled dynamics of the fluid depth and the lateral momentum. For example, the results model the dynamics of gravity currents of non-Newtonian fluids when the flow is not creeping.