变形双重介质分形油藏不稳定渗流数学模型有限元法求解研究.pdfVIP

长江大学学报(自科版) 2013年8月号理工上旬刊第lo卷第22期 of Scl No．22 JournalYangtzeUniversity(NatEdit)Aug．2013，V01．10 变形双重介质分形油藏不稳定渗流 数学模型有限元法求解研究 张 勇 (泸州医学院生物医学工程系，四川泸州646000) 731) 何国良 (电子科技大学数理学院，四川成都611 [摘要]变形双重介质分形油藏的不稳定渗流数学模型，当边暴条件为第一类边界条件时。用预估一校正 法可以很好地解决；当边界条件含有第二类边界条件时，证明差分解的存在性和收敛性时却遇到巨大困 难。为了较为简单地解决上述问题，采用有限元方法求其数值解，并证明了有限元离散解的存在性和收 敛性，为其在石油工程中的应用提供理论依据。 [关键词]不稳定渗流；收敛；有限元法；收敛性 J 22—0013—04 [中圈分类号]029 [文献标志码]A [文章■号]1673—1409(2013 在石油勘探和开采中，利用计算机技术进行油藏数值模拟是了解地下储油变化的一项简便方法·为 此，研究者提出相关模型和算法。Warren-Root模型建立在均质油藏欧几里得几何基础上，不适用于具 数口，、融，考虑压力对具有分形特征的渗透率和孔隙度的影响，建立变形双重介质分形油藏不稳定渗流数 学模型如下： 静+(卢一1)细3u=e“COeⅦ一唧冲一韶+(1一出‘升2-…。茅] (1) 譬一一Ae_一。f-1n(1一口。7)+e] (2) 相应的定解初始条件为叩I。。。=亭I。：。=o，井底定压为17 I。；。=与}}，有界定压外边界为 y(1nRtr)2 e(1nRc’r)=o’定产量为韶I。：。=一1’封闭外边界为袈l剥氓．=00 上述模型可甩来刻画具有轴称特征的地下油藏，然而由于该模型的复杂性，如何高效求解成为一个 亟待解决的问题。 1含有第二类边界条件渗流数学模型的有限元离散 当边界条件为第一类边界条件(井底定压和有界定压外边界)时，用预估一校正法[23可以很好地解 决该类定解问题；当边界条件含有第二类边界条件(定产量和封闭外边界)时，对差分解的存在性以及 收敛性的证明却十分困难。因此，为了较为简单地解决上述问题，可以采用有限元方法来求其数值解。 在不引起混淆的情况下，令t=r，r—At，则将式(2)代入式(1)并记为如下形式： o“lnR。 (3) q-fs(u,叩，￡) o≤t≤T ，1(“，7)丑3t一号≥=^(“，叩)丑au 其中： f(u，r1)一(卢一1) ^(“，节)一∞e阱2h(1一口D可)场～ ’ ，3(M，7，e)=一(1一ou)Ae阱卜扣山I土ln(1_口D17)+车l ’Laoj [收祷日期]2013一05—14 [作者俺介】张勇(1976一)男，硕士。讲师，现主要从事数据挖掘方面的教学与研究工作。 万方数据 理工上旬刊*数理科学与应用 2013年8月 定解条件上述各条件令l—lnR。，记： o：to≤t1≤…≤，一Tr：△￡一￡n—f，r1 口=L壁h：吉 矿=j7(r)矿÷=吾(矿，+矿)盔j7一=掣 矿号=丢(矿1一j7=)Et矿=