无网格方法在油藏工程中的应用初探.pdfVIP

维普资讯 2003年 1O月 油 气 井 测 试 第 12卷 第5期 无网格方法在油藏工程中的应用初探。 曾亿山 卢德唐 晏忠 良 曾清红 (中国科学技术大学力学与机械工程系) 摘要 选用有解析解的圆形油藏，采用无网格方法之一的无单元 Galerkin(Element—FreeGalerkin， 简称EFG)法求解其压力分布。并将所得结果与其解析解及有限元方法的解进行比较。阐明了使用 EFG法的求解过程，同时强调，有技巧地布点，可以提高计算速度，降低对硬件的要求。对处理裂缝 等问题，EFG法的优势更为明显。 主题词 圆形油藏 EFG法 压力分布 件如下 ： 引 言 无网格方法 (MeshlessMethod)是近年国外发展 起来的新方法。目前在国内关于这方面的研究成果 ， I I ＼ 乓 还比较少。该方法建立在移动最小二乘法 (MLS)、 式中：r——半径 ，m；p——地层压力，MPa；——油 、J = = 核函数法和单位分解法基础之上…1，不需要划分网 井半径，m；R。——油藏半径，m。 格，而是在求解区域内布置大量离散分布的节点，这 O ， ．， 些节点可以是均匀分布，也可以是随机分布，然后通 移动最小二乘法插值 ．量 过节点形式建立近似模型。这样可以消除划分网格 所带来的困难，能解决许多其它方法不能解决的问 在求解区域 Q上，函数 u()近似地表示为 题 J。特别对诸如油藏中的压裂井、水平井等复杂 ()=∑ () ()暑pT()口() (2) 边界油井数值计算 ，若用有限元法、有限体积法或有 』 = 1 式中： ()——基函数 ；口()——待求系数。 限差分法，其网格划分往往很困难，计算量大，而结 考虑下面加权的误差平方和 果却不尽人意。本文将用无网格方法之一的无单元 Galerkin(Flement-FreeGaledi【n——EF1G)方法求解油 ．，：∑w(x一研)(Uh(X)一 (研)) ， 藏中的渗流力学方程。为了便于比较，选用有解析 解的简单几何形——圆形油藏的压力 ，使用线性基 ： (一，)【p(，)口()一 (3) I i 函数和指数权函数建立形函数，运用修正变分原理 式中：n—— 的影响区域 内结点数 目； 得到离散化方程 ]。再将其结果与解析解、有限元 (一X1)——权函数 [(一X1]≠O)； 方法的解进行比较。 — — 在结点 = ，处的 值 ，，= (X1)。 (2)式对 口()求导，可得到 口()的表达式 ，代 圆形油藏模型 入(2)式得到形函数 ，即 考虑图形油藏，并假设流动是稳定的。 }()=∑ ()(A一()曰())