基于模拟有限差分的低渗透油藏非达西油水两相流动数值模拟

在模拟有限差分研究单相非线性流动基础上,进一步推广到低渗透油藏非达西油水两相流动数值模拟中,以拟线性流动模型对该方法的基本原理进行详细阐述,建立相应的数值计算格式,并采用IMPES方法对其两相流问题进行求解,通过数值算例验证方法的正确性和适用性.     

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

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

精细分层注水约束的油藏数值模拟

考虑分层注水的井筒约束方程作为内边界条件,基于黑油模型建立精细分层注水约束的油藏数值模型,形成带有加边七对角稀疏系数矩阵的封闭线性系统,采用分块矩阵乘法与预处理共轭梯度法相结合的方法进行求解.通过实例计算五点井网分层注水与笼统注水的开发效果,评价不同分段和配注模式对开发效果的影响.     

欧拉计算程序实现非均质炸药冲击起爆的数值模拟

爆轰的数值模拟研究多采用Lagrangian方法，这是因为Lagrangian网格随质点运动，介质之间的界面用滑移面来描述，故能精确地计算介质的界面。但在遇到网格扭曲和网格变形时，计算往往无法继续进行，这时就要用到网格重分等技术，这无疑增加了人为干预和计算误差。Eulerian方法能计算大变形问题，但不能精确地描述物质的界面，对于模拟平面撞击实验，爆轰波在单种介质中传播，对物质界面几乎没有要求。所以，Eulerian方法在模拟此类问题较Lagrangian方法具有明显的优势。     

基于聚类分析的降阶模型在低渗透油藏数值模拟中的应用

基于全隐式法求解考虑低渗透油藏流动过程中启动压力梯度影响的二维油水两相流动数学模型;利用本征正交分解法从获取的快照中提取基函数,将原模型投影到基函数展成的低维子空间中构造相应的降阶模型;考虑到依据等时间间隔获取的快照的缺陷,利用CVT-Lloyd算法对快照进行聚类分析.实例计算表明:通过本征正交分解法构造降阶模型能很好地对原模型进行近似;通过基于聚类分析处理快照后,一方面减少数据冗余,生成空间分布更均匀的快照集合从而提高降阶模型的精度,另一方面过少的聚类数可能导致快照信息的缺失,影响最后模型的精度.     

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

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

各向异性渗透率油藏数值模拟

建立三维三相各向异性油藏渗流数值模拟模型并应用于油藏开发.模型包含完全各向异性渗透率张量,采用十九点有限差分格式和全隐式迭代方法求解.可以描述不同区域不同主方向的各向异性渗透率分布,准确计算各向异性油藏内流体运动过程、剩余油分布和压力分布.     

基于模拟有限差分法的水驱油藏渗透率时变数值模拟

针对高含水油田注入水的长期冲刷导致储层物性参数逐渐变化从而影响油田开发及剩余油分布预测的问题,提出一种基于模拟有限差分法的水驱油藏渗透率时变模拟方法。首先拟合现场油井试井渗透率,建立通过冲刷时间表征渗透率的时变模型,然后对模型采用模拟有限差分法进行离散并采用IMPES方法对其求解,最后分析不同渗透率倍数对剩余油、压力、流线分布和开发效果的影响。结果表明：考虑渗透率时变后,水窜加剧、波及系数减小、开发效果变差。与传统的油藏描述方法相比,该方法符合水驱油藏的实际情况,能够为水驱油田后期深度挖潜提供帮助。     

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

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

基于离散裂缝模型的裂缝性油藏注水开发数值模拟

针对目前裂缝性油藏数值模拟方法存在的问题,基于单裂缝等效的概念,建立离散裂缝模型.对宏观大裂缝进行显式降维处理,在保证高计算效率的同时能够真实地反映裂缝对整个油藏流体流动的影响.详细阐述模型的基本原理,基于Galerkin加权残量法推导有限元数值计算格式,实现二维问题的数值模拟,通过算例验证模型和算法的正确性.以此为基础,分析裂缝对注水开发效果的影响.计算结果表明,对于裂缝发育程度不高尤其是当油藏中存在数条控制着流体流动方向的大裂缝时,离散裂缝模型具有很好的适用性.     

柔性管内油水环状流流动结构数值模拟

海洋中柔性管道受到外界猛烈撞击,撞击物的不同会导致管道内流道结构发生不同的改变.为探寻流道结构变化对环状流动结构的影响规律,采用N-S方程和湍流模型描述流体运动,用VoF模型与CSF模型追踪油水界面,建立数值方程与分析模型,分析流道不同形式的结构变形、不同程度变形和不同流速下油水环状流在管内的湍动能和漩涡结构演变,及环状流流型的演变和压力场的演变.分析结果表明:凹陷变形程度越大,漩涡的强度及面积越大,V形变形影响远大于○形变形影响,环状流的流型发生改变时仍然保持着水相包裹油相的稳定状态,最大压缩区域的油水界面压强也随之增大,且增速变大.所得结论对研究油水环状流在柔性管道运输中的稳定性具有一定指导意义.      

Numerical Simulation of Flow Patterns in Cyclones of Different Cone Dimensions

This paper presents the numerical simulation of the flow patterns within three cyclones of different cone dimensions. The simulation was performed with commercial CFD package FLUENT V5.5. The accuracy of the simulation was tested by comparison with the available experimental data. Simulation results show that the cone dimension significantly influences the flow pattern within cyclones. The comparison of the flow patterns within the cyclones of different cone dimensions explained well the different separation behavior of these cyclones. The results suggest that the CFD simulation is capable of predicting changes in flow patterns resulting from cyclone configuration modification.     

一维双流体模型的压力修正算法

针对一维双流体模型,通过推导双流体模型的压力修正方程,及相含率修正方程,将原本应用于单相不可压流动的压力修正系列算法推广至双流体模型求解,提出双流体模型的压力修正算法.在离散过程中运用高阶有界格式,在保证二阶以上精度的基础上克服了由相含率分布的阶跃所造成的数值结果非物理震荡.与公开发表算例进行对比,验证求解的可靠性.     

多孔介质两相流的局部守恒有限元分析

用局部守恒有限元法研究多孔介质两相渗流问题.详细阐述局部守恒有限元法的基本原理,推导两相渗流问题的局部守恒有限元计算格式并编制相应的计算程序.通过一维Buckley-Leverett两相渗流算例验证该方法的正确性.应用局部守恒有限元法和混合有限元法分别对2个模型进行分析对比.计算结果表明局部守恒有限元法具有良好的鲁棒性及适用性,相较于混合有限元法,处理过程简单,计算时间缩短,为标准有限元法应用于复杂渗流问题提供了一种途径.     

涡流二极管三维强旋湍流流动的数值模拟

为了揭示涡流二极管内三维强旋流动的特性与其阻抗性能的关系,进一步提高涡流二极管正反向流动阻抗比,本文采用Reynolds应力模型对涡流腔内非对称强旋转湍流流动进行了数值模拟,分析了腔内矢量场与标量场的分布特性,基于涡流腔内的涡量分析,得到了涡流腔内旋流的强制涡与自由涡结构及各自区域范围,并研究了涡流腔径宽比对整体性能的...     

Numerical Simulation of Swirl Flow Characteristics of CO2 Hydrate Slurry by Short Twisted Band

The development of oil and gas resources is gradually transferring to the deep sea, and the hydrate plugging of submarine pipelines at high pressures and low temperatures is becoming an important problem to ensure the safety of pipeline operations. The swirl flow is a new method to expand the boundary of hydrate safe flow. Numerical simulation of the hydrate slurry flow characteristics in a horizontal pipeline by twisted band has been carried out, and the flow of CO₂ hydrate slurry in low concentration has been simulated by the RSM and DPM models. The results show that the heat transfer efficiency is also related to Re and particle concentration. The velocity distribution has the form of symmetrical double peaks, and the peaks finally merge at the center of the pipeline. Vortexes firstly appear on both sides of the edge of the twisted band, and then move to the middle part of the twisted band. Finally, the vortex center almost coincides with the velocity center. The rotation direction of hydrate particles is the same as the twisted direction of the twisted band, twist rate (Y) is smaller, Re is larger, and the symmetric vortex lines merge farther away. The initial swirl number is mainly related to Y, but not Re. The swirl flow attenuates exponentially, and its attenuation rate is mainly related to Re, but not Y. Compared with ordinary pipelines, the swirl flow can obviously improve the transportation distance of hydrate slurry.     

往复搅拌流场中气体射流两相流数值研究

应用欧拉-欧拉双流体模型,与动网格技术相耦合数值分析往复搅拌流场中气体射流两相流的流动特性及气泡在往复流场内的分布规律。通过调整过程中吹气头的往复运动频率、气体射流速度来分析这些因素对于气泡分布和流场流动特性的影响。数值模拟结果表明,往复搅拌流场中气体射流两相流呈现与搅拌频率一致的周期性,气泡分布随吹气头位置的移动呈螺...     

二维非均匀多孔介质中不可压两相驱替的有限分析算法

为提高油藏数值模拟算法的计算效率,在求解单向稳态渗流的有限分析算法基础上,构建二维非均匀多孔介质中不可压两相渗流的有限分析算法.算法中,网格界面上的平均渗透率不是简单地取为相邻网格渗透率的调和平均值,而是通过奇点邻域解析解积分求得.相比于传统的数值算法,有限分析算法随着网格的加密,能够很快地收敛(仅需将原始网格细分至2×2或3×3),并且其计算精度和收敛性不依赖于介质的非均匀强度,从而计算效率得到提高.