低速非达西渗流动边界问题的积的分解

研究了低渗透油藏低速非达西径向流的动边界问题,给出了高精度的积分解,分析了启动压力梯度对压力分布的影响,发现启动压力梯度越大,井底附近压力下降越快,外边界传播越慢．     

基于分形理论的煤层气非达西流动分析

在考虑了煤层的分形特征和启动压力梯度影响的基础上,建立了无限大煤层中气体低速非达西流动的数学模型,并求得了量纲为一的井底压力的Laplace空间解析解,并根据数值求解结果绘制了典型的井底压力动态曲线。研究结果表明:在定产量生产的情况下,分形维数和量纲为一的启动压力梯度对早期井底压力动态无显著影响;在生产的中后期,由于受二者的影响,压力导数曲线上的径向流水平直线段消失;量纲为一的井储系数的影响主要表现在续流阶段,而吸附系数则主要影响煤基质向裂隙扩散的过渡阶段。     

考虑二次梯度项及动边界的双重介质低渗透油藏流动分析

在传统试井模型的非线性偏微分方程中根据弱可压缩流体的假设,忽略了二次梯度项,对于低渗透油藏这种方法是有疑问的.低渗透问题一个显著的特点就是流体的流动边界随着时间不断向外扩展.为了更好地研究双重介质低渗透油藏中流体的流动问题,考虑了二次梯度项及活动边界的影响,同时考虑了低渗透油藏的非达西渗流特征,建立了双重介质低渗透油藏流动模型.采用Douglas-Jones预估-校正差分方法获得了无限大地层定产量生产时模型的数值解,分别讨论了不同参数变化时压力的变化规律及活动边界随时间的传播规律,还分析了考虑和忽略二次梯度项影响时模型数值解之间的差异随时间的变化规律,做出了典型压力曲线图版,这些结果可用于实际试井分析.     

各向异性双重介质垂直裂缝井两相流体渗流

建立了各向异性双重介质中垂直裂缝井两相渗流数学模型,用有限差分法求得了其解,进行了算例分析。得到了含水饱和度沿径向、最大及最小渗透率方向的分布和水驱油前缘等饱和度面随时间的变化,分别讨论了各向异性、吸渗作用对这种介质中两相流体渗流的影响。结果表明：渗透率较在的方向,前缘推进速度较快,油井见水较早;吸渗使前缘推进速度减慢,使油井见水晚。注水开发这种油藏时,见水前注入率不能太大以充分发挥吸渗作用,存在最佳注入率。这对于水力压裂垂直裂缝井注水开发碳酸盐岩裂缝性油藏有指导意义。     

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

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

多孔介质中的非达西自然对流的分岔研究

利用分岔理论研究了多孔介质底部加热所引起的非达西自然对流。用有限差分方法计算了对流的分岔；确定了Beta数与临界瑞利数的关系。结果表明：随着Be从0增大到1，出现分岔的单胞对流的临界瑞利数Rac从39．35单调地增大到41．15。双胞对流亦有类似的趋势。这说明惯性－湍流效应有使对流稳定性增强的趋势。     

非达西渗流拟启动压力梯度推算

岩体启动压力梯度的大小是石油开采,工程防渗处理等方面需要考虑的影响因素。非线性渗流现象已被广泛所认同,低渗透性岩渗流存在启动压力梯度,它基本代表了流体产生渗流时的压力梯度大小。试验对试样进水端水压进行长期稳定控制,出水端采用精确测定流出水体积变化量的方法进行流量测定; 当渗出端水体积的变化量与时间呈直线关系时,即认为该时段渗透特征符合达西定律。通过对具有不同渗透性的软弱岩系统的试验测试,根据实测渗透系数与压力梯度的正相关关系,提出了拟启动压力梯度推算的"三次平均法",即:首先统一调整实测渗透系数值,然后确定绝对渗透系数值,再根据该值确定V0值, 3次取平均推算拟启动压力梯度的方法。分析表明:拟启动压力梯度随围压升高而增大,随岩石的绝对渗透系数降低而增大; 在渗透压差增大过程测得的拟启动压力梯度大于降低过程的拟启动压力梯度,这些都与理论分析结果相吻合。实例说明拟启动压力梯度推算"三次平均法"较为合理,易于推广使用。这为通过室内试验确定工程岩体启动渗流压力梯度的大小,提出了一种新的思路。     

岩体系统渗流场与应力场耦合的广义双重介质模型的应用研究

本文通过实验资料分析, 提出了裂隙岩体应力与渗流之间的分形几何关系式。以研究区岩体系统结构为基础, 提出了岩体系统应力场与渗流场耦合的广义双重介质模型, 并将该模型应用于某水利水电工程的渗流与控制中, 对坝区岩体进行了应力场和渗流场的耦合分析, 取得了比较满意的结果。     

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

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

逆压梯度边界层未分离情况下的湍流模式适用性研究Turbulence model applicability for boundary layer flow with adverse pressure gradient in attached case

边界层逆压梯度作用下的流动是许多工程中的一个基础问题,由于逆压梯度作用,流动形态复杂,使得数值模拟有很大的难度。基于雷诺平均纳维‐斯托克斯RANS（Reynolds Averaged Navier‐Stokes）方程对二维平板逆压梯度边界层作数值计算研究,选取6种代表性的湍流模式,得到局部摩擦系数的数值解,与实验值比较,发现k‐ω模式具有很好的精度。基于该湍流模式,给出了湍动能分布,该结果有助于认识逆压梯度边界层流动的复杂特征。     

Effect of the boundary layer on the base pressure in a two-dimensional flow with a Mach number M = 5

New data on the base pressure in a two-dimensional ow with a Mach number M = 5 are obtained for a wide range of variation of the normalized boundary-layer thickness in the flow-separation cross section. The test results are compared with Tanner’s theory, and a conclusion is made that this numerical model has to be corrected. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 23–28, May–June, 2005.     

Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability简

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.     

Turbulence structure in non-equilibrium boundary layers with favorable and adverse pressure gradients

An experimental study was conducted to document the turbulence in boundary layers on smooth walls subject to a favorable pressure gradient followed by a zero pressure gradient recovery and an adverse pressure gradient. Two component velocity profiles were acquired along the spanwise centerline of the test section, and velocity fields were obtained at the same locations in streamwise wall-normal and streamwise–spanwise planes using PIV. The FPG was shown to reduce the turbulence in the outer part of the boundary layer, reducing the transport of this turbulence and the effect of sweeps toward the wall. This reduced the inclination angle of the large structures and increased their length scale, particularly in the streamwise and spanwise directions. Recovery from the FPG to a ZPG was rapid. The APG reduced the near wall shear, resulting in a reduced effect of ejections relative to sweeps. The APG had an opposite but smaller effect on the shape and size of structures compared to the FPG.     

Numerical investigation of Dufour and Soret effects on unsteady MHD natural convection flow past vertical plate embedded in non-Darcy porous medium

The Dufour and Soret effects on the unsteady two-dimensional magnetohydro-dynamics(MHD) double-diffusive free convective flow of an electrically conducting fluidpast a vertical plate embedded in a non-Darcy porous medium are investigated numeri-cally.The governing non-linear dimensionless equations are solved by an implicit finitedifference scheme of the Crank-Nicolson type with a tridiagonal matrix manipulation.The effects of various parameters entering into the problem on the unsteady dimension-less velocity,temperature,and concentration profiles are studied in detail.Furthermore,the time variation of the skin friction coefficient,the Nusselt number,and the Sherwoodnumber is presented and analyzed.The results show that the unsteady velocity,tem-perature,and concentration profiles are substantially influenced by the Dufour and Soreteffects.When the Dufour number increases or the Soret number decreases,both the skinfriction and the Sherwood number decrease,while the Nusselt number increases.It isfound that,when the magnetic parameter increases,the velocity and the temperaturedecrease in the boundary layer.     

Low- and high-speed structures in the outer region of an adverse-pressure-gradient turbulent boundary layer

Large- and very large-scale structures in the form of elongated regions of low and high streamwise momentum have been studied in the outer region of a turbulent boundary layer subjected to a strong adverse pressure gradient. Large sets of streamwise–spanwise instantaneous velocity fields are acquired by particle image velocimetry at three wall-normal positions (0.2δ, 0.5δ, 0.8δ) at three different streamwise locations and at 0.1δ at the last streamwise location which allows us to study the wall-normal and streamwise variations of the structures. Subsequently, a pattern-recognition method and a classification scheme are employed in order to detect, classify and characterize the structures in an efficient and rigorous manner. Like in the case of zero-pressure-gradient turbulent boundary layers, long meandering streaky regions of low and high momentum are observed in the outer region of the present flow but they appear less frequently; especially in the lower part (at 0.1δ and 0.2δ) of the large-velocity-defect zone, i.e. near detachment. The dimensions of these large structures scale on boundary-layer thickness (δ) and are generally comparable to those previously reported for such structures in the overlap region of zero-pressure-gradient turbulent boundary layers. Interestingly, the adverse pressure gradient does not significantly affect the dimensions and arrangement of the large-scale structures in the upper part (at 0.5δ and 0.8δ) a segment of the outer region where the scaled Reynolds stresses also remain fairly self-similar.     

Some new results of the observation on the coherent structures of turbulent boundary layer in flow with adverse pressure gradient

By a suitable manipulation of hydrogen bubble generation, some new results were obtained: (1) The long-streaks are generated along the interfaces between low and high-speed streaks. The long-streaks are generally stretching and are moving faster than its neighboring high-speed streaks. The hydrogen bubbles in long-streaks have longer life. (2) The stream-wise vortices are also generated along the interfaces. The project is supported by the National Natural Science Foundation of China.     

An experimental study of the dynamic characteristics in boundary layer in the flow past a two-dimensional circular cylinder at subcritical reynolds numbers

An experimental study of the dynamic characteristics of flow past a two-dimensional circular cylinder is described. The fluctuationsoof wall shear stress, surface-pressure and velocity of the flow are measured with hot-film, hot-wire and pressure transducer. The frequency feature of fluctuations of wall shear stress is given. The cross-correlation functions of these fluctuations at any two points are calculated. The experimental results reveal that there is an overall syncronous fluctuation, at the shedding frequency, in boundary layer in the flow past a two-dimensional circular cylinder at subcritical Reynolds number.