缝洞型碳酸盐岩气藏产量递减特征曲线研究.docxVIP

缝洞型碳酸盐岩气藏产量递减特征曲线研究 双气体层广泛存在于缝煤岩和玄石盐酸盐岩的气藏中。程时清等人沿用沃伦-鲁特等人提出的双重介质油藏的模型, 建立了三重介质油藏三孔单渗、三孔双渗数学模型, 对其压力动态特征进行了分析。对均质、双孔介质、三重介质油气藏在定井底流压生产条件下的单井产量递减规律, 已经进行了很多常规方法研究。但这些研究, 或要求关井测试压力恢复, 这将影响油田的正常生产;或要求气井以定井底流压生产, 这在气井生产过程中难以实现。Blasingame产量递减曲线图版是现代产量递减分析中最具代表性的图版, 它克服了常规试井解释方法的缺点, 应用该方法研究三重介质气藏的生产动态数据以确定地层参数, 对缝洞型碳酸盐岩油气藏的勘探开发无疑具有重要指导意义。笔者在前人研究[1,2,3,4,5,6,7,8,9,10]的基础上, 建立考虑表皮效应的三重介质气藏产量递减模型, 绘制Blasingame产量递减曲线, 并对产量递减曲线影响因素进行分析。 1 流量模型的构建和求解 1.1 储层精细度根据 建立三重介质气藏渗流物理模型所做假设如下:①顶底不渗透、各向同性、水平等厚、径向圆形封闭气藏中心一口井, 以定井底流压生产;②气层全部打开, 径向流入井底;③单相微可压缩流体在地层中渗流服从达西定律;④地层中存在天然裂缝、基岩岩块和溶洞系统3种介质, 天然裂缝系统是主要的渗流通道, 基岩和溶洞是流体的主要储集空间, 不断向天然裂缝补给, 其补给方式为拟稳态窜流;⑤气层在打开之前, 储层各处的压力均相等, 且各种介质的孔隙度和另一介质的压力变化相对独立;⑥考虑表皮效应的影响, 忽略重力和毛管力的影响。图1为所建立的三重介质气藏渗流物理模型示意图。 1.2 拟压力函数及ct值 由以上模型假设, 定义如下无因次变量: ψ(p)=μiZipi∫pp02pμZdppjD=ψ(pi)?ψ(pj)ψ(pi)?ψ(pwf)j=m,f,vqD=1.274×10?2qscTkfh[ψ(pi)?ψ(pwf)]tD=3.6kftμir2w(?Ct)f+m+vωj=(?Ct)j(?Ct)m+v+fj=m,v,fλjf=αjkjr2wkfj=m,vrD=rrwreD=rerwψ(p)=μiΖipi∫p0p2pμΖdppjD=ψ(pi)-ψ(pj)ψ(pi)-ψ(pwf)j=m,f,vqD=1.274×10-2qscΤkfh[ψ(pi)-ψ(pwf)]tD=3.6kftμirw2(?Ct)f+m+vωj=(?Ct)j(?Ct)m+v+fj=m,v,fλjf=αjkjrw2kfj=m,vrD=rrwreD=rerw 式中,p为压力, MPa;pi为原始地层压力, MPa;μ为气体粘度, mPa·s;Z为气体偏差系数, 无量纲;ψ(p) 为标准拟压力函数;下标i表示原始;pjD为无因次压力, 无量纲;下标m表示基质;下标f表示天然裂缝;下标v表示溶洞;pwf为井底压力, MPa;下标wf表示井底;qD为无因次气产量, 无量纲;qsc为地面日产气量, 104m3/d;下标SC表示地面;T为气层温度, K;kf为裂缝渗透率, μm2;h为气层厚度, m;tD无因次时间, 无量纲;t为生产时间, h;Ct为综合压缩系数, 1/MPa; ?为孔隙度, 无量纲;rw为井半径, m;ωj为介质j的弹性储能比, 无量纲;λjf为介质j系统中的流体向裂缝系统中窜流能力大小的系数, 无量纲;αm为基质块的形状因子, 1/m2;αv为溶洞形状因子, 1/m2;r为径向距离, m;re为外边界半径, m;reD为外边界re对应的无因次距离, 无量纲;下标e表示外边界。 从物理模型出发, 根据渗流力学理论, 可建立如下无因次数学模型: 渗流微分方程: ?2pfD?r2D+1rD?pfD?rD+λmf(pmD?pfD)+λvf(pvD?pfD)=(1?ωm?ωv)?pfD?tD(1)ωm?pmD?tD=λmf(pfD?pmD)(2)ωv?pvD?tD=λvf(pfD?pvD)(3)?2pfD?rD2+1rD?pfD?rD+λmf(pmD-pfD)+λvf(pvD-pfD)=(1-ωm-ωv)?pfD?tD(1)ωm?pmD?tD=λmf(pfD-pmD)(2)ωv?pvD?tD=λvf(pfD-pvD)(3) 初始条件: pjD(rD, 0) =0j=m, v, f (4) 内边界条件: (pfD?SrD?pfD?rD)∣∣rD=1=1(5)(pfD-SrD?pfD?rD)|rD=1=1(5) 式中,S为表皮系数。 外边界条件: ?pfD?rD∣∣rD=reD=0(6)?pfD?rD|rD=reD=0(6) 无因次产量: qD=?rD?pfD?rD∣∣rD=1(7)qD=-rD?p