二维正交裂缝网络非线性渗流有限元分析方法.pdfVIP

《囊族蕊面；i磁牲渗渝葫百匾夏≯磊万暮王 二维正交裂缝网络非线性渗流的 有限元分析方法 王殿生1’2+卢占国”封b军2 (1．中国石油大学(华东)物理科学与技术学院，山东东营，257061； z，中国石油大学(华东)提高油气采收率研究中心，山东东营，257061) 摘要：当流体在裂缝中流动时，渗流速度与压力梯度之间常常呈现非缦性关系．这给裂缝渗流的数值模 拟带来了很大的困难。为了解凌裂缝非线性渗流的羲值模拟问题．在Forchbeimer公式的基础上提出了一 种简便处理非线性渗流问踅的有限元分析方法，推导了非线性渗流方程的等效积分弱形式，并对一十二维 裂缝网络模型曲非线性渗流闭建进行了嫩健模拟，计算结果与实验结果具有较好曲一致性。二者之闻酌 相对误差小于026％，表明这种处理非线性渗流问题的简便方浩是可行的和有兢的． 关键词：渗流；非线性l有限元分析法；裂缝；数值模拟 0 BI言 流体在介质中渗流时，粘滞阻力与惯性阻力比值的大小，即渗流阻力是以牯滞阻力为主 还是以惯性阻力为主，决定了渗流速度与压力梯度之间是否服从线性关系。当渗流速度较 低时，粘滞阻力占主导地位，而惯性阻力可以忽略不计，渗流速度与压力梯度之间呈线性关 系。称为线性渗流，渗流规律服从达西定律；当渗流速度增大时，惯性阻力不断增加，逐渐占 主导地位，渗流速度与压力梯度之间不呈线性关系，称为非线性渗流。达西定律不再适用。 介质的结构类型不同，所产生的粘滞阻力不一样。多孔介质的粒问孔踪产生粘滞阻力大，而 裂缝介质的裂缝产生粘滞阻力小。因此，在相同的压力梯度下，一般裂缝介质中的渗流速度 要比多孔介质中的快得多，在裂缝介质中更容易呈现非线性渗流。 FiroozabadI等“3分析了高速气 国内外学者对非线性渗流问题进行了一些研究。Abbas 体通过多孔介质的流动特性，提出了一个描述非线性渗流的三次多项式；J．s．Andrade 等”1在二维不规则多孔介质内进行了Navi”stokes方程的数值计算，研究了渗流偏离经典 diBelh 达西定律的原因；Ha Brinkman牿度项的非线性描述方程；于龙“3通过对岩体裂隙水流运动规律的研究，导出了 多研究均表明，在描述流体特别是液体的非线性渗流规律的各种公式中，Forchheimer公式 不仅具有较好的理论基础和明确的物理意义，而且具有更广泛的适用性。 本文采用Forchh￡unef公式描述非线性渗流规律，建立了单相不可压缩液体的稳定非 ·作者简介：王殿生．男。196s年生，黑龙江饭兰县入。教授。1988年毕韭于石油大学(华东)应用物理学 专业，现从事教学和缝洞型油藏渗流机理研究工作。 f四l————_＼查鱼±塑竺苎篓查叁塾墨璺堡垒堡垒叁墨i 线性渗流方程，提出了一种简便的处理非线性渗流问题的有限元分析方法，应用虚位移原理 推导渗流方程的等效积分弱形式，并对一个二维正交裂缝网络物理实验模j塞!进行了数值模 拟计算，验证了这种简便处理方法的可行性和有效性。 l非线性渗流方程 的关系： 一譬=∞．+幻f (1)“’ 一茁2“。十钿f 式中，p为压力；警表示渗流速度”-方向上的压力梯度；n，6表示与介质和流体性质相关的 系数泪的一次项代表了粘滞阻力的影响，二次项代表了惯性阻力的影响。 考虑二维情况，单相不可压缩液体在剐性介质中稳定渗流的连续性方程可表示为： 。 (2)。 等+篝=oaz。av 根据式(1)和式(2)可推导出： 。f一。+√蹈]。f一。+√霹] f 26 f f 26 f (3)