用边界元分析任意形状裂缝性油藏不稳定渗流.pdfVIP

2007 2007年5月 西安石油大学学报(自然科学版) May ofXi’all Science Vd．22No．sl 第22卷第s1期 Jomud ShiyouUniversity(NaturalEdition) 文章编号：1673—064X(2007)sl-0043-03 用边界元分析任意形状裂缝性油藏不稳定渗流 何应付，刘先贵，杨正明 (中国科学院渗流流体力学研究所，河北廊坊065007) 摘要：在考虑裂缝系统各向异性的基础上，建立了任意开头双孔隙度双重介质油藏不稳定渗流的数 学模型，利用坐标转换、Laplace变换和边界元方法获得油藏内任意点的压力，进而由杜哈美原理得 到考虑井筒储存和表皮效应的井底压力．绘制了考虑各向异性、复杂边界以及油藏内存在不渗透区 域等因素影响的井底不稳定压力典型曲线，并分析了曲线特征．与解析方法的结果对比表明，本文 的边界元方法计算精度很高，适用于任意形状双重介质油藏的不稳定压力动态分析． 关键词：边界元方法；双重介质油藏；各向异性；典型曲线； 中图分类号．-TE312文献标识码：A 常见的边界元方法主要用于解决均质油藏的渗 流问题【1-4]，虽然文献[1]采用边界元方法研究了 双孔隙度双重介质油藏井底压力的计算方法，但其 警=等_0，(刎)童蹦=2’3，．．，，m； 并没有考虑裂缝系统的各向异性以及油藏内不渗透 (5) 率区域的影响．为此，本文采用边界元方法分析了裂 缝系统各向异性和不渗透区域双重介质油藏的压力 P。=力=Pi，t=0 (6) 动态特征的影响． 其中：11l是油藏外边界，且r1—1U11卜2=11，／-1—1 n 1数学模型 界，且nnrs=0，(i≠_『)；m是油藏边界的个数； 艿是DiracDelta函数，用来表明源(汇)的位置；z，y 假设渗流场内有源(汇)分布，地层水平、等厚， 表示直角坐标；K为渗透率，矗；／z为流体黏度，Pa·s； 流体单相微可压缩，拟稳态窜流，则考虑裂缝系统各 q1为源(汇)z的强度，l／s；声为孔隙度；c为裂缝压缩 向异性含有不渗流区域的双重介质油藏二维不稳定 系数，1／Pa；t表示时间，s．下标，表示裂缝，m表示 渗流数学模型为 基质，i表示初始状态． ‰磐+％奢+％(锄一力)+ 标变化，并无因次化和Laplace变换，则有 2 (1) J¨7巧D g秘(z—xt)3(y一抑)=概警； ．0a2硼p29+案卅s)珈 aX毛。a踞 (2) flK。m(夕。一办)|+惦警=o； (7) ：导宝口D庐(XD—y肼)艿(yD一‰)． 3 (3) f=I P。=力=P