專利名稱：不依賴徑向流直線段的試井分析方法
技術(shù)領(lǐng)域：
本發(fā)明涉及油田試井領(lǐng)域所用的方法，屬于不依賴徑向流直線段的試井分析方法。
可采用的技術(shù)方案是該試井分析方法可根據(jù)不同的壓力降落、壓力恢復(fù)(穩(wěn)定后的關(guān)井恢復(fù)、未穩(wěn)定后的關(guān)井恢復(fù)、擬穩(wěn)定后的關(guān)井恢復(fù))測(cè)試情況，采用k^n=70.6qμh[∂Δpw(tn)/∂lnt]]]>計(jì)算出實(shí)測(cè)壓力資料的“瞬時(shí)滲透率” ；然后利用rn=2.0C1k^ntnφμct]]>計(jì)算出探測(cè)半徑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]]]>計(jì)算出壓降(或恢復(fù))時(shí)的Δ(1/kn)；最后利用kn=11kn-1+Δ(1kn)]]>計(jì)算出滲透率kn，便可得到地層的滲透率分布k(r) kn，rn-1≤r＜rn。
該方法的有益效果是使用該方法編制出的試井分析解釋軟件，可為非均質(zhì)油藏測(cè)壓資料的解釋提供新的解釋手段，打破了現(xiàn)有手段利用試井曲線直線段特征進(jìn)行分析的解釋方法，為現(xiàn)場(chǎng)大量不出半對(duì)數(shù)直線的試井資料的解釋提供了手段，不僅能從實(shí)際測(cè)壓資料解釋出地層的滲透率分布，也能給出地層的壓力分布、平均地層壓力等，既適用于不出現(xiàn)徑向流試井資料的解釋，同時(shí)也能適合于有徑向流直線段的試井資料解釋，對(duì)油田開(kāi)發(fā)中后期大量不出半對(duì)數(shù)直線的試井資料的解釋及如何評(píng)價(jià)增產(chǎn)措施效果提供了一種有效的分析方法。縮短工作時(shí)間，大大提高了工作效率，適合推廣應(yīng)用。
利用式3計(jì)算Kn的方法如下定義r0＝rwr0D＝1 (6)Δ(1ki)=1ki-1ki-1---(7a)]]>Δ(1k1)=1k1---(7b)]]> 利用式6至式8，當(dāng)t＝tn時(shí)式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)]]>經(jīng)進(jìn)一步簡(jiǎn)化，可得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為存在可靠的壓力導(dǎo)數(shù)值的第一時(shí)間點(diǎn)，經(jīng)推導(dǎo)式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中的積分項(xiàng)在考慮到式29時(shí)，可以寫成∫1rjD2K(r′D,t^D)dr′D=∫1rjDπr′Dt^Dexp(-rD′22t^D)W12,12(rD′2t^D)dr′D---(19)]]>進(jìn)行如下變量替換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的研究成果表明當(dāng)時(shí)間不在區(qū)間 時(shí)，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)]]>即根據(jù)不同的壓力降落、穩(wěn)定后的關(guān)井恢復(fù)、未穩(wěn)定后的關(guān)井恢復(fù)、擬穩(wěn)定后的關(guān)井恢復(fù)測(cè)試情況，采用k^n=70.6qμh[∂Δpw(tn)/∂lnt]---(5)]]>式5計(jì)算出實(shí)測(cè)壓力資料的“瞬時(shí)滲透率”
然后利用rn=2.0C1k^ntnφμct---(4)]]>式4計(jì)算探測(cè)半徑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計(jì)算壓降(或恢復(fù))時(shí)的Δ(1/kn)；最后，利用kn=11kn-1+Δ(1kn)---(30)]]>式30計(jì)算滲透率kn，便可得到地層的滲透率分布k(r) kn，rn-1≤r＜rn。
實(shí)施例一、對(duì)大慶油田芳190-124井的解釋成果圖，這是打破目前依賴徑向流直線段解釋的試井分析方法，也就是對(duì)不出徑向流直線段的試井資料進(jìn)行分析，能從實(shí)際測(cè)壓資料解釋出地層的滲透率分布，也能給出地層的壓力分布、流度分布，這是常規(guī)方法不能作到的，如

圖1、2、3所示。
實(shí)施例二對(duì)大慶油田芳226-100井的解釋成果圖，如圖4、5、6所示，說(shuō)明使用該方法分析時(shí)能給出隨井徑變化的滲透率分布、流度分布及壓力分布，而常規(guī)方法只能給出平均的滲透率、平均地層壓力等。
實(shí)施例三、對(duì)大慶油田北2-d3-470井的解釋成果圖，如圖7、8、9所示，說(shuō)明使用該方法分析時(shí)能給出隨井徑變化的滲透率分布、流度分布及壓力分布，而常規(guī)方法只能給出平均的滲透率、平均地層壓力等。
可見(jiàn)，使用該方法可為非均質(zhì)油藏測(cè)壓資料的解釋提供新的解釋手段，打破了現(xiàn)有手段利用試井曲線直線段特征進(jìn)行分析的解釋方法，為現(xiàn)場(chǎng)大量不出半對(duì)數(shù)直線的試井資料的解釋提供了手段。
權(quán)利要求
1.一種涉及油田試井領(lǐng)域所用的不依賴徑向流直線段的試井分析方法，其特征在于a.根據(jù)不同的壓力降落、壓力恢復(fù)(穩(wěn)定后的關(guān)井恢復(fù)、未穩(wěn)定后的關(guān)井恢復(fù)、擬穩(wěn)定后的關(guān)井恢復(fù))測(cè)試情況，采用k^n=70.6qμh[∂Δpw(tn)/∂lnt]---(5)]]>式5計(jì)算出實(shí)測(cè)壓力資料的“瞬時(shí)滲透率” b.利用rn=2.0C1k^ntnφμct---(4)]]>式4計(jì)算探測(cè)半徑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計(jì)算壓降(或恢復(fù))時(shí)的Δ(1/kn)；d.最后，利用kn=11kn-1+Δ(1kn)---(30)]]>式30計(jì)算滲透率kn，便可得到地層的滲透率分布k(r) kn，rn-1≤r＜rn，其中式中的Ct——地層綜合壓縮系數(shù)，MPa-1；h——油層厚度，ft； ——瞬時(shí)滲透率，md； q——井底流量，RB/D；r——徑向半徑，ft； t——時(shí)間，days；μ——粘度，cp； φ——孔隙度，小數(shù)。
全文摘要
一種涉及油田試井領(lǐng)域所用的不依賴徑向流直線段的試井分析方法。解決了現(xiàn)有試井解釋方法無(wú)法對(duì)不出徑向流直線段的試井資料進(jìn)行分析的問(wèn)題。其特征在于根據(jù)不同的實(shí)測(cè)壓力降落、壓力恢復(fù)資料，計(jì)算出實(shí)測(cè)壓力資料的“瞬時(shí)滲透率”k
文檔編號(hào)E21B47/00GK1438404SQ0310266
公開(kāi)日2003年8月27日 申請(qǐng)日期2003年2月17日 優(yōu)先權(quán)日2003年2月17日
發(fā)明者張同義, 楊永智, 張奇斌 申請(qǐng)人:大慶油田有限責(zé)任公司