水驱粘塑性非牛顿原油渗流特性的数值模拟

针对水驱粘塑性非牛顿原油在地层下呈强非线性渗流特征,建立了描述其流动过程的数学方程; 采用全隐式离散、牛顿－ 拉弗松方法迭代线化,稳定双共轭梯度方法求解,对特定初边条件下二维油藏油水两相流动问题作了数值模拟．结果表明,水驱非牛顿流体驱油效率明显低于水驱牛顿流体     

幂律流体不稳定径向流的数值解法

对微可压缩非牛顿幂律流体在均质多孔介质中的非稳定径向渗流方程进行了回顾和对比，用有限差分方法建立了微分方程的数值解法．研究了时间步长、空间步长对数值解精度和稳定性的影响，并结合算例绘出了不同流性指数的瞬时压力曲线．该研究结果可用于非牛顿流体渗流机理研究及试井分析．     

两相渗流驱动问题有限元—特征线修正差分格式及理论分析

研究多孔介质中可压缩可混溶两相渗流驱动问题的计算方法：压力方程用有限元方法求解，饱合度方程用特征线修正有限差分方法求解，构造了全离散数值计算格式，证明了最佳收敛阶     

两相渗流驱动问题有限元——特征线修正差分格式及理论分析

研究多孔介质中可压缩可混溶两相渗流驱动问题的计算方法，压力程用有限元方法求解，饱合度方程用特征线修正有限差分方法求解，构造了全离散数值计算格式，证明了最佳收敛阶。     

基于线法的功能梯度材料断裂分析

为了克服材料非均匀性引起的数值困难,一种半解析数值方法——线法,被引入功能梯度材料的断裂分析.通过有限差分将问题的控制方程半离散为定义在沿梯度方向离散节线上的常微分方程组,然后应用B样条高斯配点法求解该常微分方程组.为了演示线法功能梯度材料断裂分析的方法,给出了指数梯度含裂纹功能梯度材料板分别在恒定位移、弯曲载荷作用下...     

一种非均匀网格上的高精度紧致差分格式

提出了一种形式简单、网格剖分灵活、具有一定通用性的非均匀网格上的三点四阶紧致差分格式,对格式的截断误差进行了分析.采用文中提出的格式对Burgers方程和对流方程进行数值求解,并与均匀网格上的三点四阶紧致差分格式所得数值解对比,结果证明本文提出的格式对于大梯度问题的数值模拟有更高的精度.     

真实气体高速非达西流的单井数值模拟

在考虑天然气高速非达西渗流特征和真实气体PVT参数随压力变化基础上,建立了平面径向流下单井高速非达西不稳定渗流数学模型,模型中气体渗流符合Forchheimer二项式方程,并引入了天然气PVT参数的计算和拟压力的计算,采用数值差分方法对模型进行求解,并编制了计算软件。计算了非达西渗流以及真实气体PVT参数随压力变化对产能的影响,为气藏产能计算和动态分析提供依据。     

解渗流问题的一种新算法及分析

提出了解微可压缩两相渗流驱动问题的一种新算法,压力方程用混合有限元求解,饱合度方程用特征线修正有限差分方法求解,证明了最优误差收敛阶。     

无相间物质传递化学驱浓度方程算子分裂隐式解法

为了改进化学驱数学模型 U TCHEM显式求解浓度方程计算速度慢、计算结果精度低的缺点 ,研究了隐式求解组份浓度方程的方法。根据具有无相间物质传递关系的化学驱油藏流体渗流过程满足的相行为 ,推导出了化学驱数学模型 UTCHEM物质守恒方程的等价形式 :饱和度方程和组份浓度方程。利用算子分裂技术将组份浓度方程分裂为扩散方程和对流方程 ,隐式交替求解对流方程和扩散方程得到组份浓度方程的隐式解。扩散方程采用隐式局部一维格式差分离散 ,利用追赶法求解 ;对流方程选用了隐式迎风格式差分 ,并且结合油藏模拟问题的流场是有势场的特点 ,实现了对流方程隐式差分显式求解。所建立的隐式求解浓度方程的方法提高了计算精度 ,可以加大计算时间步长 ,加快计算速度。     

对流扩散方程的二次单调插值特征差分方法

特征差分方法适合于求解对流占优扩散问题,但也存在着精度低或非物理振荡等缺点.为克服其不足,建立和研究了二次单调插值特征差分方法.此方法具有较高精度并消除了非物理振荡.方法是无条件稳定的.用交替分组显式(AGE)方法求解了差分方程,方法便于并行计算.数值结果表明,方法可成功地求解对流占优扩散问题.     

裂缝多孔介质中达西流动的有限差分方法

用差分方法来模拟二维裂缝多孔介质中的单相达西流动问题.采用降维模型对二维区域内的裂缝进行建模,相比整个区域而言裂缝的宽度很小并且把裂缝视为一个一维的界面,流体通过裂缝会和周围多孔介质发生耦合现象.在裂缝和周围多孔介质中,流体流动均遵循达西定律和守恒定律.采用差分方法来求解降维模型中推导的流体流动方程.通过数值实验验证了该方法的有效性,并证明了裂缝是快速通道还是地质屏障取决于裂缝处渗透率张量的大小.     

可压多孔介质流的扩展混合元解法

讨论多孔介质中两种可压缩流体混溶驱动问题数值方法,假定介质是各向异性的,渗透率系数为张量形式。压力方程采用扩展混合元方法求解压力变量、梯度变量,以及速度变量;浓度方程采用标准有限元方法求解,这一方法对各向异性渗透率多孔介质流可以获得更可靠的数值解。构造了半离散数值格式,通过理论分析得到了压力、速度以及浓度等变量的最优L2模误差估计,对浓度变量获得了H1模最优误差估计。     

有启动压力梯度渗流的数学模型

围绕低渗透问题和Bingham流体渗流问题，从推广的达西定律入手，详细推导了其渗流控制方程，讨论了初始条件和各种边界条件．并且考虑到表皮效应和井筒储集常数，引入有效井筒半径，将控制方程和边界条件无量纲化，建立了有启动压力梯度渗流问题的数学模型，为有启动压力梯度油藏开采提供理论依据．并给出了无网格方法的数值模拟结果．     

非牛顿幂律流体球向不定常渗流问题的近似解

文献１用影响半径的概念和平均值方法求得了弱压缩非牛顿幂律流体球向不定常渗流问题的近似解．此解不满足连续性方程．本文在更加严密的理论基础上求得了近似满足所有基本方程和内、外边界条件的新的近似解．对于影响半径，近似解和文献１的近似解的相应公式的结构是相似的．当流态指数ｎ＝１时，无因次压力函数和文献１的相应公式基本一致．     

混合边界油藏压力分布的有限元法模拟

建立了稳定渗条件下的有限元方程，针对定边界和封闭边界合作用下的渗流问题进行了模拟计算，给出了圆形地层和矩形地层在以上条件下的压力分布，结果表明，混合油藏外边界将产生复杂的地层压力分布，压力分布随地层方位的不同而不同，同时说明，有限元法用于求解复杂边界条件的渗流问题是方便有效的。     

完全可压熔混流驱动问题的特征差分方法及其理论分析

研究多孔介质中完全可压缩渗流驱动问题的计算方法,对压力方程用有限差分离散,对浓度方程用基于局部二次插值的特征线差分方法离散,构造并分析了数值计算格式,证明了最优H_1模误差估计。     

双重介质分形油藏井筒耦合模型及数值模拟研究

针对双孔介质地质特征和油气的实际流动过程,同时考虑油藏的分形特征以及油藏和井筒耦合流动的影响,将双重介质分形油藏渗流与井筒管流作为整体进行研究。建立了双重介质分形油藏中非牛顿幂律流体径向渗流与井筒耦合流动模型。并采用Douglas-Jones预估-校正法求得模型的数值解,最后探讨压力动态变化规律,绘制压力图版。     

竖直圆柱空间多孔介质中非牛顿流体自然对流

建立了封闭圆柱空间多孔介质中非牛顿流体自然对流的数学模型，通过理分析和数值计算，获得了在不同达西－瑞利数Ｒａ和不同流变指数ｎ下，非牛顿流体的流变特性，Ｎｕ随Ｒａ和ｎ的变化规律，以及Ｒａ和ｎ对流场和温度场的影响。     

两点边值抛物型耦合方程组的数值方法研究

研究一类非线性两点边值问题,其方程具有抛物型特征,速度场与温度切需耦合求解．提出了该类方程组的两种数值求解方法．以幂次流体沿竖直平板层流自然对流为例,进行了数值模拟,并将求解方法与“盒式”方法进行了对比．结果表明,数值方法适合该类非线性两点边值抛物型耦合方程组的求解,并具有推导简单、数值求解易于收敛,且计算稳定的特点．     

渗流方程的渗透率自适应权重网格粗化算法

为了提高渗流方程的计算速度和精度,将渗透率自适应网格技术应用于三维非均匀非稳态渗流方程的网格粗化算法中。对于渗透率或孔隙度变化异常的区域,采用精细网格直接求解其压力分布;而在其他区域,采用不均匀网格粗化的方法计算其压力分布。用自适应权重网格粗化算法计算了三维非均匀非稳态渗流场的压力分布。结果表明,三维非均匀非稳态渗流方程的三维不均匀自适应网格粗化算法的解在渗透率或孔隙度异常区域的压力分布规律非常逼近精细网格算法的解,在其他区域的压力分布规律非常逼近粗化算法的解。与采用精细网格算法相比,其计算速度大大提高。