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

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

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

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

分形油藏非牛顿幂律流体低速非达西不稳定渗流的组合数学模型

结合分形理论与渗流理论,对分形油藏非牛顿幂律流体低速非达西不稳定渗流的试井分析问题的数学模型进行了推导.该分形油藏模型由内域为非牛顿幂律流体低速非达西渗流,外域为非牛顿幂律流体达西渗流的同心圆域组成.在考虑井筒储存、表皮效应影响下,建立了该油藏的不稳定渗流有效井径组合数学模型,在3种外边界条件下求出了两个区域内压力在Laplace空间的解析解,应用Stehfest数值反演方法求得井底的无因次压力,分析了井底压力动态特征和参数影响.非牛顿幂律流体的幂律指数、分形参数均对典型曲线产生较大的影响,呈现出与牛顿流体和均质油藏明显不同的特征.这对非均质油藏非牛顿流体的不稳定试井分析及研究其非线性渗流特征均十分重要.     

致密油藏微纳米喉道中的边界层特征

通过引入吸引力修正耗散粒子动力学（DPD）方法,实现流体和固体的相互吸引作用,模拟纳米喉道中的微尺度流动,探讨边界层的产生机理,结合微圆管实验,定量表征微纳米喉道中边界层的特征,明确微纳米喉道中边界层的影响因素.研究发现：分子尺度,热运动对速度影响很大;超过分子尺度,压差占主导作用.热运动使粒子在原位置振动,不改变粒子的整体移动方向.随着喉道半径的增大,泊肃叶流动的抛物线特征越来越明显.边界层厚度受压力梯度、喉道半径和流体粘度的影响.当压力梯度增大或流体粘度减小时,边界层厚度增大;当喉道半径减小时,边界层厚度先增大后减小.边界层厚度是导致非线性渗流特征的根本原因.随着边界层厚度增大,非线性渗流特征越来越明显.     

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

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

稠油油藏试井理论研究

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

聚合物驱两相流试井参数敏感性

考虑油水两相、生产历史、油藏平面非均质性、井筒储存和表皮效应等因素,建立了生产历史阶段聚合物驱数学模型和不稳定试井阶段的流线模型,用流管法对解释模型进行了数值求解.研究表明:随着油水粘度比的增大,压力及压力导数曲线向上平移,随着生产时间的增加,储层的有效渗透率降低,当高渗透条带沿主流线方向分布时,注水井压降导数曲线反映不出油水前缘的影响,而随着聚合物注入浓度的增大,压力导数曲线下凹出现的越来越早.     

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

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

各向异性油藏水平井多角度人工裂缝线性单元计算方法

考虑各向异性油藏渗透率张量的表征,利用Green函数和拉普拉斯变换建立裂缝流动的一维单元,裂缝的流量分布采用节点线性插值,裂缝内的流动处理为线性积分,耦合地层与人工裂缝的流动,建立有限导流裂缝井底压力的求解方法.结果表明:多裂缝压裂水平井存在压裂裂缝线性流、地层线性流、系统径向流3种流动形态,压裂裂缝条数越多,相同的生产时间,无因次井底压降越小;裂缝条数对流动影响明显.随着裂缝条数的增加,压降变化减小;裂缝长度和导流能力有相似的变化.人工裂缝与井筒角度越大,产能越大,当裂缝垂直于井筒时,产量最大;地层最大渗透率方向垂直于人工裂缝时产量最大,平行裂缝时产量最小.当人工裂缝垂直于井筒,并同时垂直于地层最大渗透率方向时,达到最大产量值.     

用拉普拉斯变换差分法求解试井分析中的一维渗流问题

结合拉普拉斯变换和有限差分法给出求解试井分析中一维渗流问题的拉普拉斯变换差分法:首先对渗流方程采用拉普拉斯变换消去时间变量得到拉普拉斯空间数学模型,采用有限差分法求解拉普拉斯空间数学模型,最后通过拉普拉斯反演算法得到井底压力或产量.通过与有限差分法结果和解析解对比,拉普拉斯变换差分法比有限差分法计算误差小.虽然单步计算耗时长,但计算任意时刻结果时对空间网格的适应性和不依赖其它时刻计算结果的特性使得拉普拉斯变换差分法在试井分析中有非常好的应用前景.     

Numerical investigation of a coupled moving boundary model of radial flow in low-permeable stress-sensitive reservoir with threshold pressure gradient

The threshold pressure gradient and formation stress-sensitive effect as the two prominent physical phenomena in the development of a low-permeable reservoir are both considered here for building a new coupled moving boundary model of radial flow in porous medium. Moreover, the wellbore storage and skin effect are both incorporated into the inner boundary conditions in the model. It is known that the new coupled moving boundary model has strong nonlinearity. A coordinate transformation based fully implicit finite difference method is adopted to obtain its numerical solutions. The involved coordinate transformation can equivalently transform the dynamic flow region for the moving boundary model into a fixed region as a unit circle, which is very convenient for the model computation by the finite difference method on fixed spatial grids. By comparing the numerical solution obtained from other different numerical method in the existing literature, its validity can be verified. Eventually, the effects of permeability modulus, threshold pressure gradient, wellbore storage coefficient, and skin factor on the transient wellbore pressure, the derivative, and the formation pressure distribution are analyzed respectively.     

聚合物驱试井分析叠加原理适用性的探索

建立考虑聚合物剪切变稀特性的聚驱试井数学模型,基于矩形网格和井筒周围径向加密的复合网格。采用有限体积方法求得数值解,将数值解与叠加原理得到的注聚井停注压降进行对比。结果表明：叠加原理计算得到的注聚井关井井底压力值远低于数值解,说明叠加原理不能用于聚合物驱试井解释。在实测注聚井压降资料解释中,叠加原理得到的渗透率明显偏低,进一步证明了叠加原理不适合非线性聚驱试井模型,数值方法更加适用于聚驱试井解释。研究结论适用于类似的非线性试井问题。     

高应变率拉伸下纯铝中氦泡长大的动力学研究

 分析了高应变率加载下纯铝中氦泡长大的动力学过程，给出了含内压氦泡长大的动力学方程，并且分别研究了氦泡内压、基体材料惯性、粘性、表面张力以及基体环境温度对初始半径为1 nm氦泡长大的影响。研究结果表明：（1）初始内压可以促使氦泡快速长大，当氦泡直径超过1 μm时，内压对氦泡长大的影响可以忽略不计。（2）表面张力在氦泡整个长大过程中的影响都很小。（3）材料惯性对氦泡长大起抑制作用，并且随着氦泡半径的增长，抑制效应越来越明显。（4）在所有因素中，温度对氦泡长大的影响最为明显，温度升高，材料的粘性降低，氦泡的内压增加，促使氦泡加速长大。     

Oscillation and Migration of Bubbles within Ultrasonic Field

The oscillation and migration of bubbles within an intensive ultrasonic field are important issues concerning acoustic cavitation in liquids.We establish a selection map of bubble oscillation mode related to initial bubble radius and driving sound pressure under 20 kHz ultrasound and analyze the individual-bubble migration induced by the combined effects of pressure gradient and acoustic streaming.Our results indicate that the pressure threshold of stable and transient cavitation decreases with the increasing initial bubble radius.At the pressure antinode,the Bjerknes force dominates the bubble migration, resulting in the large bubbles gathering toward antinode center,whereas small bubbles escape from antinode.By contrast,at the pressure node,the bubble migration is primarily controlled by acoustic streaming,which effectively weakens the bubble adhesion on the container walls,thereby enhancing the cavitation effect in the whole liquid.     

RHIC能区Au+Au碰撞中椭圆流的横半径依赖

用RQMD模型模拟了质心系束能量为s_NN=200GeV的Au+Au非对心碰撞, 研究了椭圆流对末态粒子冻出位置的关系. 研究表明, 随着冻出半径的增加, 椭圆流逐渐增强, 在源的初始表面处达到最大值, 随着冻出半径的继续增加, 椭圆流逐渐减小. 椭圆流的这种变化趋势反映了随着半径的增加, 压力梯度在各个方向上不同的变化. 本文提出了利用平均径向速度分析压力梯度的方法.