不依賴徑向流直線段的試井分析方法
2023-05-05 01:05:36
專利名稱：不依賴徑向流直線段的試井分析方法
技術領域：
本發明涉及油田試井領域所用的方法，屬於不依賴徑向流直線段的試井分析方法。
可採用的技術方案是該試井分析方法可根據不同的壓力降落、壓力恢復(穩定後的關井恢復、未穩定後的關井恢復、擬穩定後的關井恢復)測試情況，採用k^n=70.6qh[pw(tn)/lnt]]]>計算出實測壓力資料的「瞬時滲透率」 ；然後利用rn=2.0C1k^ntnct]]>計算出探測半徑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＝rwr0D＝1 (6)(1ki)=1ki-1ki-1---(7a)]]>(1k1)=1k1---(7b)]]> 利用式6至式8，當t＝tn時式3變成1k^n=12K(rD,t^D)drD[1k1+(1k2)u(rD-r1D)+(1kD)u(rD-r2D)]]>+...+(1ki)u(rD-ri-1D)+...+(1kn)u(rD-rn-1D)]drD---(9)]]>整理後得1k^n=1k112K(rD,t^D)drD+(1k2)r1D2K(rD,t^D)drD]]>+(1k3)r2D2K(rD,t^D)drD+...+(1ki)ri-1D2K(rD,t^D)drD]]>+...+(1kn)rn-1D2K(rD,t^D)drD---(10)]]>注意到12K(rD,t^D)drD---(11)]]>rjD2K(rD,t^D)drD=12K(rD,t^D)drD-1rjD2K(rD,t^D)drD---(12)]]>=1-1rjD2K(rD,t^D)drD]]>把式10寫為1k^n=1k1+(1k2)[1-1r1D2K(rD,t^D)drD]+(1k3)[1-1r2D2K(rD,t^D)drD]+....]]>+(1ki)[1-1ri-1D2K(rD,t^D)drD]++(1kn)[1-1rn-1D2K(rD,t^D)drD]---(13)]]>經進一步簡化，可得1k^n=i=1n-1(1ki)-i=2n-1(1ki)1ri-1D2K(rD,t^D)drD---(14)]]>+(1kn)[1-1rn-1D2K(rD,t^D)drD]]]>由式7a，有i=1n-1(1ki)=1kn-1---(15)]]>代入到式14，得1k^n=1kn-1-i=2n-1(1ki)1ri-1D2K(rD,t^D)drD---(16)]]>+(1kn)[1-1rn-1D2K(rD,t^D)drD]]]>利用上式就可以得到{1/Kn}，首先求出1k1=1k^(t1)---(17)]]>
其中t1為存在可靠的壓力導數值的第一時間點，經推導式16有(1kn)=1k^n-1kn-1+i=2n-1[(1ki)1ri-1D2K(rD,t^D)drD]1-1rn-1D2K(rD,t^D)drD---(18)]]>如果已知 ，利用式18便可求出Δ(1/Kn)，而式18中的積分項在考慮到式29時，可以寫成1rjD2K(rD,t^D)drD=1rjDrDt^Dexp(-rD22t^D)W12,12(rD2t^D)drD---(19)]]>進行如下變量替換z^D=rDt^DdrD=t^Ddz^D---(20)]]>rDt^Dt^Ddz^D=rDt^Ddz^D=z^Ddz^D---(21)]]>rD1z^D1t^Dz^0D---(22)]]>rDrjDz^DrjDt^Dz^jD---(23)]]>把式20～式23代入到式19，得1rjD2K(rD,t^D)drD=z^0Dz^jDz^Dexp(-z^D22)W12,12(z^D2)dz^D---(24a)]]>1rjD2K(rD,t^D)drD=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.6qh[pw(tn)/lnt]---(5)]]>式5計算出實測壓力資料的「瞬時滲透率」
然後利用rn=2.0C1k^ntnct---(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.6qh[pw(tn)/lnt]---(5)]]>式5計算出實測壓力資料的「瞬時滲透率」 b.利用rn=2.0C1k^ntnct---(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日
發明者張同義, 楊永智, 張奇斌 申請人:大慶油田有限責任公司