油藏渗流问题的无网格法分析

无网格法作为一种新型的数值方法,因其近似函数不依赖于网格而受到广泛关注.基于无网格法对油藏渗流问题的求解进行了研究.对无网格法的基本原理进行了详细的阐述;并针对油藏单相渗流问题推导了无网格法计算格式,通过实例计算验证了该方法的有效性.     

三维蒸汽辅助重力驱油过程的自适应网格法

针对蒸汽辅助重力热采稠油过程中,蒸汽腔周围存在的物理量变化十分剧烈的锋面及计算计算网格必须足够的小的问题,为了有效提高对蒸汽辅助重力驱油问题的数值计算效率,研究了考虑井筒与油藏间相互耦合的三维均匀油藏的蒸汽辅助重力驱油问题的自适应网格方法,自适应网格法和均匀精细网格法之间的计算结果对比表明自适应网格法可以大幅度地提高对于蒸汽辅助重力驱油这一复杂流场的计算速度并具有很好的精度。     

利用油藏工程原理描述优势渗流通道的新方法

针对水驱砂岩油藏高含水期优势渗流通道普遍存在的问题,提出了应用油藏工程原理来识别优势渗流通道并计算其体积的方法。该方法是利用水驱特征曲线对优势渗流通道进行识别,并利用模糊聚类分析对优势通道窜流程度进行分类;然后利用剩余可采储量计算出优势渗流通道体积,进而计算出优势通道的渗透率和孔喉半径。利用VB语言编制了相应的软件,并用该软件对胜坨油田T143区块的优势渗流通道进行了识别和计算,通过其中三口井的示踪剂监测结果证实了该方法的准确性。     

径向基函数配点法求解边坡降雨瞬态渗流问题

根据土壤水分特征曲线Gardner模式推导线性化非饱和渗流Richard方程,应用无网格法中的多元二次径向基底函数配点法对空间域进行离散,利用龙格库塔法对时间域进行离散. 通过离散点满足控制方程式与边界条件,建立求解非饱和渗流问题的数值模型. 考虑土层组合以及降雨强度,计算分析非饱和无限边坡降雨入渗的瞬态渗流场的变化情况,得到边坡内不同时刻、不同深度的孔隙水压力. 其数值计算结果与解析解相符,较有效地解决了传统数值方法在模拟非饱和土渗流产生的数值病态问题.     

瞬态有自由面渗流分析的不动网格──高斯点有限元法

对于瞬态有自由面渗流问题，建立了一种简单、实用的不动网格─—高斯点有限 元法。该法与变网格法及已有的不动网格法相比，计算量小、格式简单。算例表明了 该法的有效性。     

低渗亲水油藏剖面数值模拟计算参数的研究

介绍了长庆油田低渗亲水砂岩油藏剖面数值模拟方法．建立油水两相剖面数值模拟的地质模型，然后根据地质模型中的渗透率级别，将渗流参数无因次化为流动水饱和度的有效模拟参数．经计算验证，这一方法的计算精度较高．对认识低渗亲水油藏的水驱油机理，具有现实意义．     

低渗透油藏弹性驱和溶解气驱采收率计算

针对低渗透油藏流体渗流存在启动压力梯度,基于达西定律的常规采收率计算方法难以适用的问题,结合非达西渗流理论,采用油藏工程方法,提出了考虑启动压力梯度的低渗透油藏弹性驱和溶解气驱采收率计算新方法.介绍了该方法的基本原理和计算过程,并进行了算例分析.研究结果表明:低渗透油藏启动压力梯度的存在降低了弹性驱和溶解气驱采收率,加密井网是提高弹性驱和溶解气驱低渗透油藏开发效果的有效途径.提出的计算方法拓展了非线性渗流理论的应用,可用于低渗透油田的开发设计.     

非线性渗流条件的低渗油藏产能计算方法

为了准确评价低渗透油藏的产能动态规律,针对低渗透油藏的渗流特点,建立了考虑非线性渗流特征的产能计算模型,并用Newton-Raphson方法对模型进行求解.结果表明:用启动压力梯度方法计算的产量将比用非线性渗流方法计算的产量低,特别是在井底流压较大时;而用达西方法计算的产量比用非线性渗流方法计算的产量高,特别是在井底流压较小时.非线性渗流方法比传统的启动压力梯度方法更能够反映低渗透油藏的渗流规律.     

油藏水平井数值模拟研究中一种新的求解方法

提出了底水油藏水平井油水两相渗流问题的物理模型和数学模型，并建立了差分方程，然后利用ＩＭＰＥＳ求解方法的思路，并对其不足之处进行了改进，从而形成了一种新的求解方法．这种方法是采用上游网格饱和度插值计算相邻网格间的水量交换，用本阶段时间中点的饱和度对压力方程进行再迭代．保持了ＩＭＰＥＳ方法简单方便、占计算机内存少、计算速度快的优点，而精度明显提高，结果更趋合理．实际计算表明，新方法效果良好，可广泛用于底水、边水或人工注水油藏油水两相渗流问题的数值模型研究，对水平井和垂直井均适用．     

径向基点插值无网格法与有限元耦合法

为了发挥无网格法和有限元法各自的优势,提出了径向基点插值无网格法与有限元直接耦合的计算方法.比较了点插值无网格方法与无网格Galerkin法(element-free Galerkin method, EFG)的优缺点:无网格法只需要结点信息,无需单元信息,克服了有限元计算中网格畸变和重新生成带来的困难,故其在分析裂纹扩展和局部大变形等问题方面具有优势.对悬臂梁受集中荷载问题进行了计算,两种计算结果一致,但是本算法效率较高.用本方法与FEM耦合方法对半无限平面受集中荷载问题进行了计算,结果表明:本文的耦合方法可行,计算效率较高.     

冲击动力学中无网格法应用

无网格法是一种新兴发展的数值计算方法,作为有限元法的一种重要补充在冲击动力学领域得到了广泛的应用.本文从求解过程上对无网格法和有限元法进行了比对分析,针对SPH,RKPM,EFG3种典型无网格法的基本原理及其在冲击动力学方面的应用进行了阐述,并提出无网格法的发展方向.     

Helmholtz问题具径向基函数的无网格法

构造了Helmholtz方程具径向基函数的无网格方法.通过引入多种径向基函数构造了Galerkin型的无网格方法.文末给出了数值算例,并与有限元方法进行了比较.讨论了无网格方法的数值精度以及径向基函数中参数对其数值解的影响.结果表明具径向基函数的Galerkin型无网格方法是求解Helmholtz方程的一种有效且精度高的方法.     

基于节点敏感性分析的无网格法节点布置研究

在应用无网格法进行数值计算时,由于节点的布置对于计算结果的精度有直接影响,因此节点布置方案以及节点性态特点是无网格法中的研究重点。把有限元和形状优化设计中的节点敏感性这一概念应用于无网格法,选择势能密度作为响应变量对无网格法中的节点进行敏感性分析,得到了节点布置方案与节点敏感性系数之间的关系。通过悬臂梁的数值算例,说明了该方法能够有效地指导并改进节点布置方案,使计算精度明显提高。     

冲击动力学中无网格法应用概述

无网格法是一种新兴发展的数值计算方法,作为有限元法的一种重要补充在冲击动力学领域得到了广泛的应用。本文从求解过程上对无网格法和有限元法进行了比对分析,针对SPH、RKPM、EFG三种典型无网格法的基本原理及其在冲击动力学方面的应用进行了阐述,并提出无网格法的发展方向。     

关于配点型无网格法边界条件处理技术

由于无网格数值方法具有传统的有限差分法和有限元法不可比拟的优点,着重介绍了配点型无网格法格式及其特点.在总结配点型无网格法处理导数边界条件的各种技术的基础上,提出了基于积分插值的新处理技术.通过对基于点插值的配点型无网格法解Helmholtz问题的研究,验证了该技术的优越性.     

浅水不规则波的无网格数值模拟

采用基于配点和径向基函数的无网格数值模拟方法,建立数值海洋工程水池,将源点布置在模拟波浪场之外,沿边界布置配点而不是划分网格;对深水不规则波经斜坡过渡后在均匀水深浅水域中传播的过程进行无网格数值模拟和波谱的迭代修正.不规则波的数值模拟结果显示,该方法形式简单、计算速度快、稳定性好.经物理水池的造波实验验证,波浪谱、波浪参数、波浪时历都与实验结果吻合良好,并具有浅水波特征,最终模拟得到的浅水不规则波令人满意.     

无网格法中一种新型权函数研究和应用

文章基于概率论和误差理论,提出了一种新的权函数——正态权函数(Normal weight function),从理论和实践上证明了它的实用、可行性。通过一维杆和二维梁实例,把正态权函数与现有流行的权函数进行比较,说明它是一种受影响域半径变化的影响较小、高效实用的权函数。最后给出了正态权函数影响域半径的确定规则。     

次弹性材料的无网格分析方法研究

次弹性材料在实际工程中是常见的,传统计算中大多数采用有限元方法。利用无网格伽辽金法对次弹性材料进行数值计算,并通过罚参数来实现本质边界条件,推导了求解此类问题的无网格伽辽金法离散格式。算例结果表明,无网格伽辽金法处理次弹性材料时,具有较高的计算精度,是一种有效的数值计算方法。     

具有动态接触边界自动识别功能的刚(粘)塑性无网格RKPM法

将刚(粘)塑性流动理论与再生核质点法(RKPM)相结合,提出了基于刚(粘)塑性不可压缩材料的无网格RKPM法.分别采用边界奇异权法和修正的罚函数法处理本质边界条件和体积不可压缩条件,根据刚(粘)塑性材料的不完全广义变分原理,推导了金属塑性成形过程无网格RKPM法数值模拟的刚度方程.将动态接触边界自动识别技术引入无网格法.最后,对金属塑性成形过程进行了数值模拟,并将模拟结果与刚塑性有限元商品软件Deform2D及实验结果作了比较.无网格RKPM法计算结果与有限元结果和实验结果均吻合良好,表明了该方法的正确性.     

弹性薄板弯曲问题的无网格解法

基于弹性薄板的线性弯曲理论和重调和方程的一般解理论构造的基本解,建立无网格法求解薄板弯曲问题的数值计算格式.采用径向基函数近似表示横向分布荷载,获得问题的特解,而齐次解答由基本解的线性组合得到.将边界条件用于确定未知系数,获得可以数值求解的线性方程组.在均布荷载情况下,计算并给出四边简支矩形薄板弯曲解,和解析解进行比较,证明无网格方法的收敛性和计算精度.