三维两相流数值试井模型

建立了基于混合PEBI网格的三维两相流数值试井模型，内容包括：物理模型的建立，网格剖分方法，控制方程的离散格式以及内、外边界的处理方法。用单相流双重介质和双渗介质解析解对模型的计算结果进行了验证。对3层油藏两相流试井曲线特征分析发现，多层两相流动试井曲线特征与单相流动双重介质的曲线特征相似。因此。有可能寻求一个参数关系将双重渗透介质试井问题与三维两相流试井问题联系起来，从而减小多相流试井分析的难度和多解性。径向PEBI混合网格与数值试井模型的结合应用，给出了三维两相流、有井筒存储和表皮效应的不稳定渗流产量隐式表达式。     

压裂水平井多裂缝系统不稳定压力动态特征

针对考虑无限导流裂缝的情况下,由于问题的内边界复杂化,水平井多段压裂形成的多裂缝系统,难以直接获得解析解,采用解析解与数值解相结合的半解析法构造计算模型,既保证有解析法的高精度,又具有数值方法的灵活性。对油气藏水平井压裂后的压力动态特征进行计算和不稳态压力特征进行分析和讨论,其不稳态压力特征的认识对于分析该类油气井的动态以及试井解释具有重要指导作用。     

三维两相流数值试井模型

建立了基于混合PEBI网格的三维两相流数值试井模型 ,内容包括 :物理模型的建立 ,网格剖分方法 ,控制方程的离散格式以及内、外边界的处理方法。用单相流双重介质和双渗介质解析解对模型的计算结果进行了验证。对 3层油藏两相流试井曲线特征分析发现 ,多层两相流动试井曲线特征与单相流动双重介质的曲线特征相似。因此 ,有可能寻求一个参数关系将双重渗透介质试井问题与三维两相流试井问题联系起来 ,从而减小多相流试井分析的难度和多解性。径向PEBI混合网格与数值试井模型的结合应用 ,给出了三维两相流、有井筒存储和表皮效应的不稳定渗流产量隐式表达式     

基于混合PEBI网格的复杂井数值模拟应用研究

针对笛卡儿正交网格在描述复杂井型以及复杂油藏边界时存在的局限性,提出了综合运用径向网格、笛卡儿网格和六边形PEBI网格的混合PEBI网格计算方法.该方法中,在直井区域采用径向网格,可实现网格体积由小到大的变化;在复杂井区域采用笛卡儿网格,可充分把握井筒附近的流动状态;在油藏区域采用六边形PEBI网格,可更好地逼近复杂油藏边界和断层的形态.应用结果表明,混合PEBI网格具有高度的灵活性,与有限差分五点格式相比,有效降低了网格取向效应,能够满足精细油藏数值模拟的要求.     

油水两相渗流压降数值试井模型的建立

数值试井问题不同于常规油藏数值模拟问题。介绍了数值试井模型中采用的非结构网格技术 ,给出了在非结构网格下油水两相流动控制方程的离散形式 ,而后耦合了井筒存储和表皮系数公式 ,得到了压力降落数值试井问题的产量模型。将压力法与数值试井模型进行了对比校验 ,发现随油井的生产 ,压力导数曲线呈上翘趋势。通过计算结果对比 ,对原压力法提出的无量纲方法进行了修正。修正后的双对数曲线均可与同参数的单相流双对数曲线重合 ,边界反应也可以从两相流的压降试井曲线中获得。     

基于OpenFOAM的周期性管内流动边界处理方法研究

为探究开源计算流体力学软件OpenFOAM中不同周期性边界条件处理方法的适用性,针对圆管内泊肃叶流动,分别在mapped,cyclicAMI,cyclic fan边界条件下开展数值模拟,对模拟得到的速度分布和压力分布进行分析.结果表明:3种边界条件设置下,数值模拟得到的速度分布均与解析解吻合,其中mapped和cyclicAMI边界条件下的数值解非常接近,且与解析解吻合更好;采用mapped边界能够得到与解析解符合的压力场,采用cyclic fan边界时在模拟区域进、出口存在很大的压力梯度,采用cyclicAMI边界得到的压力则恒为0.基于此,探讨了mapped边界条件下,数值模拟结果对计算网格的敏感性.研究发现:随着径向网格数目的增加,速度的计算精度先迅速提高而后趋于稳定;边界层层数的增加能够有效提高近壁区速度的计算精度;数值模拟时径向网格数目采用16个以上为宜,边界层层数不宜少于3层.     

基于移动粒子半隐式方法波传播模型的声传播数值求解方法研究

为扩展无网格移动粒子半隐式法(MPS)的应用范围,将其扩展到声学计算应用领域,在流声分离假设下,基于拉格朗日描述建立了声场控制方程,即MPS方法下的声波传播模型(MPS-WP)。实现了声学硬边界和吸收边界的离散粒子模型,研究了CFL数及粒子间距等关键参数对计算精度的影响,发现在一定的CFL数和粒子间距内,该模型可以保持较小的误差水平。采用高斯脉冲传播算例验证了MPS-WP在硬边界与吸收边界下的二维声场计算,声压的计算与解析解吻合较为良好,展示了该求解方法的有效性和准确性。在此基础上,对静止流场多个脉冲叠加及均匀来流下声波的传播特性进行了模拟和分析,计算结果与解析解吻合良好,为复杂流动条件下的声学计算提供了一种新的研究方法。     

底水稠油油藏水平井含水规律预测方法

针对底水稠油油藏水平井出水位置难以确定的问题,通过建立底水稠油油藏物理及数学模型,利用Green函数、Newman乘积方法和叠加原理推导底水稠油油藏水平井非均匀产液的压力响应解析解,采用Stehfest数值反演算法得到考虑井筒储集效应和表皮效应的底水稠油油藏水平井分段产液试井井底压力解.通过绘制水平井分段试井典型图版,进行流动阶段划分及分析,进而得到不同生产段长度、生产段数目、生产段位置分布、流量分布等影响因素影响规律,并绘制相关解释图版.矿场应用实例验证了研究结果的有效性及广阔的应用前景.     

基于多重网格法的MT正演模型边界截取

为了验证基于多重网格法边界截取的效果,通过数值模拟,分别对大地电磁法问题的均匀半空间模型和三层模型进行基于多重网格法的边界截取研究和分析.采用加权平均法产生粗网格上系数矩阵,并与传统的方法对比.研究结果表明:粗网格加权平均法的渐进收敛性明显比Galerkin法(采用简单的插值矩阵)的好,但比粗网格近似的几何法的渐进收敛性略差;收敛性的广义Fourier谱分析在大地电磁法中对介质的非连续性不敏感,不能区分模型参数变化所引起的收敛性变化;随着边界截取的增加,满足特定收敛标准的数值解与理论解偏差增大;当截取到的深度约为趋附深度的2倍时,截取后的计算结果接近理论值,说明基于多重网格法的边界截取有效.     

基于2.5维有限元算法的边界处理及网格划分

由于2.5维周期化技术能有效地处理列车振动波传播模拟过程中计算尺度与计算时间的矛盾,因而国内外学者在此理论基础上将FE(有限元)和粘滞边界以及IFE(无限元)等边界处理方法进行耦合.为讨论上述算法中有限元网格划分区域及单元大小对计算精度影响,通过运用软件Matlab编制程序读取Ansys中的模型信息,并将2.5维有限元和无限元边界耦合等算法在其中得以实现.在此基础上发现在相同的网格划分条件下,采用2.5维无限元边界较粘滞阻尼边界具有更高的精度,同时通过选取不同尺寸的有限元模拟区域及最大网格大小,将计算结果与经典解析解进行比较,并对上述网格划分参数的合理选择提出了建议.