不依赖径向流直线段的试井分析方法

专利名称：不依赖径向流直线段的试井分析方法
技术领域：
本发明涉及油田试井领域所用的方法，属于不依赖径向流直线段的试井分析方法。
可采用的技术方案是该试井分析方法可根据不同的压力降落、压力恢复(稳定后的关井恢复、未稳定后的关井恢复、拟稳定后的关井恢复)测试情况，采用k^n=70.6qμh[∂Δpw(tn)/∂lnt]]]>计算出实测压力资料的“瞬时渗透率” ；然后利用rn=2.0C1k^ntnφμct]]>计算出探测半径rn；再利用Δ(1kn)=1k^n-1kn-1+Σi=2n-1[Δ(1ki)∫z^0Dz^i-1DΩ(z^D)dz^D]1-∫z^0Dz^n-1DΩ(z^D)dz^D]]]>计算出压降(或恢复)时的Δ(1/kn)；最后利用kn=11kn-1+Δ(1kn)]]>计算出渗透率kn，便可得到地层的渗透率分布k(r) kn，rn-1≤r＜rn。
该方法的有益效果是使用该方法编制出的试井分析解释软件，可为非均质油藏测压资料的解释提供新的解释手段，打破了现有手段利用试井曲线直线段特征进行分析的解释方法，为现场大量不出半对数直线的试井资料的解释提供了手段，不仅能从实际测压资料解释出地层的渗透率分布，也能给出地层的压力分布、平均地层压力等，既适用于不出现径向流试井资料的解释，同时也能适合于有径向流直线段的试井资料解释，对油田开发中后期大量不出半对数直线的试井资料的解释及如何评价增产措施效果提供了一种有效的分析方法。缩短工作时间，大大提高了工作效率，适合推广应用。
利用式3计算Kn的方法如下定义r0＝rwr0D＝1 (6)Δ(1ki)=1ki-1ki-1---(7a)]]>Δ(1k1)=1k1---(7b)]]> 利用式6至式8，当t＝tn时式3变成1k^n=∫1∞2K(r′D,t^D)dr′D[1k1+Δ(1k2)u(r′D-r1D)+Δ(1kD)u(r′D-r2D)]]>+...+Δ(1ki)u(r′D-ri-1D)+...+Δ(1kn)u(r′D-rn-1D)]dr′D---(9)]]>整理后得1k^n=1k1∫1∞2K(r′D,t^D)dr′D+Δ(1k2)∫r1D∞2K(r′D,t^D)dr′D]]>+Δ(1k3)∫r2D∞2K(r′D,t^D)dr′D+...+Δ(1ki)∫ri-1D∞2K(r′D,t^D)dr′D]]>+...+Δ(1kn)∫rn-1D∞2K(r′D,t^D)dr′D---(10)]]>注意到∫1∞2K(r′D,t^D)dr′D---(11)]]>∫rjD∞2K(r′D,t^D)dr′D=∫1∞2K(r′D,t^D)dr′D-∫1rjD2K(r′D,t^D)dr′D---(12)]]>=1-∫1rjD2K(r′D,t^D)dr′D]]>把式10写为1k^n=1k1+Δ(1k2)[1-∫1r1D2K(r′D,t^D)dr′D]+Δ(1k3)[1-∫1r2D2K(r′D,t^D)dr′D]+....]]>+Δ(1ki)[1-∫1ri-1D2K(r′D,t^D)dr′D]+…+Δ(1kn)[1-∫1rn-1D2K(r′D,t^D)dr′D]---(13)]]>经进一步简化，可得1k^n=Σi=1n-1Δ(1ki)-Σi=2n-1Δ(1ki)∫1ri-1D2K(r′D,t^D)dr′D---(14)]]>+Δ(1kn)[1-∫1rn-1D2K(r′D,t^D)dr′D]]]>由式7a，有Σi=1n-1(1ki)=1kn-1---(15)]]>代入到式14，得1k^n=1kn-1-Σi=2n-1Δ(1ki)∫1ri-1D2K(r′D,t^D)dr′D---(16)]]>+Δ(1kn)[1-∫1rn-1D2K(r′D,t^D)dr′D]]]>利用上式就可以得到{1/Kn}，首先求出1k1=1k^(t1)---(17)]]>
其中t1为存在可靠的压力导数值的第一时间点，经推导式16有Δ(1kn)=1k^n-1kn-1+Σi=2n-1[Δ(1ki)∫1ri-1D2K(r′D,t^D)dr′D]1-∫1rn-1D2K(r′D,t^D)dr′D---(18)]]>如果已知 ，利用式18便可求出Δ(1/Kn)，而式18中的积分项在考虑到式29时，可以写成∫1rjD2K(r′D,t^D)dr′D=∫1rjDπr′Dt^Dexp(-rD′22t^D)W12,12(rD′2t^D)dr′D---(19)]]>进行如下变量替换z^D=rDt^D⇒drD=t^Ddz^D---(20)]]>rDt^Dt^Ddz^D=rDt^Ddz^D=z^Ddz^D---(21)]]>rD→1⇒z^D→1t^D≡z^0D---(22)]]>rD→rjD⇒z^D→rjDt^D≡z^jD---(23)]]>把式20～式23代入到式19，得∫1rjD2K(r′D,t^D)dr′D=∫z^0Dz^jDπz^Dexp(-z^D22)W12,12(z^D2)dz^D---(24a)]]>∫1rjD2K(r′D,t^D)dr′D=∫z^0Dz^jDΩ(z^D)dz^D---(24b)]]>其中z^0D=1t^D---(25)]]>z^jD=rjDt^D---(26)]]>且Ω(z^D)=πz^Dexp(-z^D22)W12,12(z^D2)---(27)]]>Oliver的研究成果表明当时间不在区间 时，K(rD，tD)≈0，因此z^0D=min{1t^D,0.12}---(28)]]>将式24b代入到18中，得Δ(1kn)=1k^n-1kn-1+Σi=2n-1[Δ(1ki)∫z^0Dz^i-1DΩ(z^D)dz^D]1-∫z^0Dz^n-1DΩ(z^D)dz^D]---(29)]]>而，kn=11kn-1+Δ(1kn)---(30)]]>即根据不同的压力降落、稳定后的关井恢复、未稳定后的关井恢复、拟稳定后的关井恢复测试情况，采用k^n=70.6qμh[∂Δpw(tn)/∂lnt]---(5)]]>式5计算出实测压力资料的“瞬时渗透率”
然后利用rn=2.0C1k^ntnφμct---(4)]]>式4计算探测半径rn；再利用Δ(1kn)=1k^n-1kn-1+Σi=2n-1[Δ(1ki)∫z^0Dz^i-1DΩ(z^D)dz^D]1-∫z^0Dz^n-1DΩ(z^D)dz^D]---(29)]]>式29计算压降(或恢复)时的Δ(1/kn)；最后，利用kn=11kn-1+Δ(1kn)---(30)]]>式30计算渗透率kn，便可得到地层的渗透率分布k(r) kn，rn-1≤r＜rn。
实施例一、对大庆油田芳190-124井的解释成果图，这是打破目前依赖径向流直线段解释的试井分析方法，也就是对不出径向流直线段的试井资料进行分析，能从实际测压资料解释出地层的渗透率分布，也能给出地层的压力分布、流度分布，这是常规方法不能作到的，如

图1、2、3所示。
实施例二对大庆油田芳226-100井的解释成果图，如图4、5、6所示，说明使用该方法分析时能给出随井径变化的渗透率分布、流度分布及压力分布，而常规方法只能给出平均的渗透率、平均地层压力等。
实施例三、对大庆油田北2-d3-470井的解释成果图，如图7、8、9所示，说明使用该方法分析时能给出随井径变化的渗透率分布、流度分布及压力分布，而常规方法只能给出平均的渗透率、平均地层压力等。
可见，使用该方法可为非均质油藏测压资料的解释提供新的解释手段，打破了现有手段利用试井曲线直线段特征进行分析的解释方法，为现场大量不出半对数直线的试井资料的解释提供了手段。
权利要求
1.一种涉及油田试井领域所用的不依赖径向流直线段的试井分析方法，其特征在于a.根据不同的压力降落、压力恢复(稳定后的关井恢复、未稳定后的关井恢复、拟稳定后的关井恢复)测试情况，采用k^n=70.6qμh[∂Δpw(tn)/∂lnt]---(5)]]>式5计算出实测压力资料的“瞬时渗透率” b.利用rn=2.0C1k^ntnφμct---(4)]]>式4计算探测半径rn；c.利用Δ(1kn)=1k^n-1kn-1+Σi=2n-1[Δ(1ki)∫z^0Dz^i-1DΩ(z^D)dz^D]1-∫z^0Dz^n-1DΩ(z^D)dz^D]---(29)]]>式29计算压降(或恢复)时的Δ(1/kn)；d.最后，利用kn=11kn-1+Δ(1kn)---(30)]]>式30计算渗透率kn，便可得到地层的渗透率分布k(r) kn，rn-1≤r＜rn，其中式中的Ct——地层综合压缩系数，MPa-1；h——油层厚度，ft； ——瞬时渗透率，md； q——井底流量，RB/D；r——径向半径，ft； t——时间，days；μ——粘度，cp； φ——孔隙度，小数。
全文摘要
一种涉及油田试井领域所用的不依赖径向流直线段的试井分析方法。解决了现有试井解释方法无法对不出径向流直线段的试井资料进行分析的问题。其特征在于根据不同的实测压力降落、压力恢复资料，计算出实测压力资料的“瞬时渗透率”k
文档编号E21B47/00GK1438404SQ0310266
公开日2003年8月27日 申请日期2003年2月17日 优先权日2003年2月17日
发明者张同义, 杨永智, 张奇斌 申请人:大庆油田有限责任公司