黎 明 周福建 田亞凱 李明輝 汪 杰

1. 中國石油大學(北京)非常規(guī)油氣科學技術研究院, 北京 102200;2. 中國石油大學(北京)油氣資源與探測國家重點實驗室, 北京 102200

0 前言

致密儲層油氣資源豐富且開采潛力巨大,是未來全球油氣資源開發(fā)的主力軍。然而,由于致密儲層復雜的巖石礦物組成和孔隙結構,導致在開發(fā)油氣資源的過程中儲層滲透率應力敏感性(即滲透率隨孔隙流體壓力的變化規(guī)律)極其復雜。致密砂巖儲層應力敏感性評價一直是國內外學者研究的焦點,但至今仍未得到統(tǒng)一的結論。

國內外研究者[1-5]評價儲層應力敏感性主要是通過變圍壓測試,得到有效應力與滲透率關系數據,再根據應力敏感模型公式計算應力敏感系數和滲透率損失率。因此,確定有效應力和選取應力敏感模型及公式是評價儲層應力敏感的關鍵。

如何將實驗室結果和模型計算或生產參數設計相結合,是應力敏感研究的關鍵問題。多年來,研究者提出多種模型來關聯實驗結果與內在規(guī)律[6-8]，這些模型推導的前提為假設儲層的滲透率較高,并且一般都有1個或多個參數需要確定,所以在擬合低滲、特低滲儲層時與實驗數據有一定差距,存在分析誤差。為填補這一空白,本文以鄂爾多斯盆地致密砂巖油氣藏巖石實驗為基礎,從滲透率應力敏感的物理概念出發(fā),推導了一個修正的冪律滲透率模型,可以集成到儲層模擬器中,有助于提高致密砂巖儲層中的油氣產量預測精度。

1 致密砂巖儲層應力敏感數學模型

在儲層生產過程中,巖石所受有效應力增大,巖石骨架顆粒及裂縫端面發(fā)生變形,進而引發(fā)滲透率變化。其他因素包括流體的靜水壓力和動水壓力將導致裂縫面擴展,以及儲層流體的化學潛蝕和弱化作用同樣會對滲透率造成影響。但目前的研究無法將上述因素納入考慮范圍,主要還是通過滲流實驗來獲得巖石滲透率參數。在考慮外部圍壓和內部孔隙流體壓力相互作用的情況下,依據實驗測試數據,以有效應力為變量,建立有效應力與滲透率模型的數學關系。

綜合多位研究者的研究,通過擬合巖心實驗數據,得到現有的有效應力與滲透率模型的數學關系，主要分為5種，見表1。

表1中的模型都是在一定條件下針對不同巖樣推出的,故有其自身的適用范圍,在范圍之外擬合并不準確。對于低滲儲層來說,之前的學者在一定精度范圍內運用了冪律和指數模型擬合,而筆者推導的修正冪律模型對于低滲儲層有著更好的擬合效果。

致密砂巖滲透率敏感性的影響因素主要包括巖石的礦物含量、微裂縫和孔縫的發(fā)育連通等,在考慮低滲透致密地層的特征情況下可以推得改進模型,下文將介紹模型的推導過程。

按之前學者的研究,巖石所受有效應力與其滲透率之間關系可以用式(1)表示：

(1)

(2)

從圖1巖心電鏡掃描照片可以發(fā)現,致密砂巖巖心微裂縫發(fā)育和孔喉連通情況較差。對于一般儲層可以使用的冪律式(1)計算,但對于致密砂巖巖心,需要推導出更復雜的模型。

b)盒8組巖心,放大1 000倍b)He 8 cores,magnified 1 000 times圖1 巖心電鏡掃描照片Fig.1 SEM images of core

將式(1)移項取對數整理可得:

(3)

(4)

代入式(3)可得：

(5)

f(lnkσeff)=f(lnkσinit)+f′(lnkσinit)(lnσ-lnσ)+…+

(6)

(7)

整理得：

kσeff

(8)

(lnσ-lnσ)

(9)

整理得:

(10)

2 實驗方法

致密砂巖致密低滲,常規(guī)氣測穩(wěn)態(tài)法很難進行測定且結果不準確,本實驗使用壓力脈沖衰減法滲透率測試儀PDP-200,儀器結構見圖2[14]。

圖2 PDP-200儀器結構圖Fig.2 PDP-200 device structure diagram

在生產過程中,孔隙流體壓力將顯著降低,而上覆壓力將不會發(fā)生大變化。這將導致巖石承受比初始應力更高的有效應力,有效應力可以通過式(11)計算：

σ=p-αp

(11)

國內外研究者的大量研究結果表明[15-21],Biot系數α受巖石物性影響,而在致密低滲儲層中,α值與1相近。本文中的研究區(qū)塊滲透率極低,故近似取α值為1,即巖石所受有效應力等于圍壓減去孔隙流體壓力。將實驗孔隙流體壓力設為5 MPa,不斷增加圍壓,使有效應力依次為5 MPa、10 MPa、15 MPa、20 MPa、25 MPa、30 MPa、35 MPa、40 MPa。在每個壓力點應力老化20 min后進行滲透率測試。實驗所用致密砂巖巖心的基礎數據見表2。

表2 實驗參數及結果表

表2中實驗所用巖心主要取自鄂爾多斯長7組(以下簡稱長7組)和蘇格里盒8組(以下簡稱盒8組),其中ch、hh和cv、hv分別是取自水平和垂向層理方向,8塊巖樣均表現出低孔、低滲特征。其中致密砂巖油藏長7組孔隙度為4.11%～5.23%,滲透率為 0.032 9×10-3～0.094 5×10-3μm2。致密砂巖氣藏盒8組孔隙度為2.19%～2.81%,是長7組的一半,滲透率為0.000 7×10-3～0.001 8×10-3μm2,大致比長7組低一個數量級,滲透率在平行層理方向為垂直層理方向的2～3倍。

3 實驗結果分析與討論

式(10)能很好地對實驗數據進行回歸分析,運用最小二乘法擬合曲線,將改進模型與舊模型的實驗數據比較,結果見圖3。并用回歸系數R來評價擬合效果,R值0～1,越接近1擬合效果越好,Marrow冪律模型的回歸系數R值為0.964～0.985,平均值為0.977,其值低于改進模型(R值為0.990～0.995,平均值為0.993)。進一步分析:1)長7組水平巖心見圖3-a),改進模型回歸系數R=0.992,高于Marrow冪律模型(R=0.974);2)長7組垂向巖心見圖3-b),改進模型回歸系數R=0.990,高于Marrow冪律模型(R=0.964);3)盒8組水平巖心見圖3-c),改進模型回歸系數R=0.994,高于Marrow冪律模型(R=0.985);4)盒8組垂向巖心見圖3-d),改進模型回歸系數R=0.995,高于Marrow冪律模型(R=0.983)。通過對比明顯可以得出,巖心層理水平方向滲透率較高,為巖心垂向層理方向2～3倍本文提出的改進模型比Marrow冪律模型具有更高的回歸精度。

從圖3可看出,當有效應力為40 MPa時,長7組巖心滲透率降低至10-5μm2數量級,盒8組巖心滲透率降幅更顯著,低至10-7μm2數量級。在原始地層應力狀態(tài)下,水平和垂向層理面的巖心滲透率相差1～2倍,但隨著有效應力增大,滲透率差逐漸減小,當有效應力在 40 MPa 時,二者基本相同。這說明相比垂向層理面,水平層理面的致密砂巖滲透率應力敏感性更強?？梢缘贸?致密砂巖應力敏感性傷害更容易發(fā)生在水平滲透率方向。

為驗證改進冪律擬合效果,模型取雙對數擬合,見圖4。

a)ch巖心回歸擬合比較a)Comparison of ch core regression fitting

b)cv巖心回歸擬合比較b)Comparison of cv core regression fitting

c)hh巖心回歸擬合比較c)Comparison of hh core regression fitting

d)hv巖心回歸擬合比較d)Comparison of hv core regression fitting圖3 巖心滲透率回歸擬合比較圖Fig.3 Comparison of regression fitting of core permeability

a)ch巖心雙對數擬合比較a)Comparison of ch core double logarithmic fitting

b)cv巖心雙對數擬合比較b)Comparison of cv core double logarithmic fitting

c)hh巖心雙對數擬合比較c)Comparison of hh core double logarithmic fitting

d)hv巖心雙對數擬合比較d)Comparison of hv core double logarithmic fitting圖4 巖心滲透率雙對數擬合比較圖Fig.4 Comparison of double logarithmic fitting of core permeability

由圖4可看出，盒8組巖心擬合效果>長7組巖心擬合效果,需考慮到盒8組比長7組更致密(滲透率低一個數量級)。本文推導的改進模型對低滲透率的巖心擬合性能更優(yōu)。

儲層巖石骨架和填隙物具有不同硬度和力學性質，石英、長石等剛性礦物可壓縮性小,而云母、黏土等塑性礦物易受力變形破碎。巖石受力時,剛性礦物承壓,孔喉體積壓縮,塑性礦物破碎產生的細小顆粒運移進入巖石孔隙或堵塞孔喉,導致樣品滲透率和孔隙度降低。巖心滲透率和有效應力呈現出很好的冪律相關。

運用改進模型擬合表1中巖心所得模型參數見圖5。其中,回歸的ηeff為負,范圍為-0.84～-0.22。而回歸的為正值,范圍為0.01～0.08。絕對值比ηeff低一個數量級。圖5展示出了每個巖心的ηeff和結果。

圖5 ηeff和 η的回歸結果Fig.5 Regression results of ηeff and

圖6展示了運用改進模型擬合表1巖心的雙對數曲線,擬合直線的斜率由ηeff和共同作用,其絕對值越大表示應力敏感效應越強。斜率的絕對值大小排列依次為ch>cv>hh>hv,致密砂巖油藏應力敏感大于致密砂巖氣藏應力敏感,結果與實驗數據相同。

圖6 改進模型雙對數曲線Fig.6 Double logarithmic curve of the improved model

4 結論

本文應用PDP-200壓力脈沖衰減法滲透率測試儀研究了致密砂巖油氣藏應力敏感。致密砂巖油藏巖心取自鄂爾多斯盆地延長組長7組,致密砂巖氣藏巖心取自蘇格里氣田盒8組,取心包括水平和垂向層理兩個方向。推導出一個冪律改進模型,用于研究更加致密的儲層條件，得出以下結論。

1)致密砂巖油藏長7組巖心初始滲透率比致密氣藏盒8組巖心初始滲透率高一個數量級,并且在高有效應力下依然保持穩(wěn)定。表明長7組巖心驅替阻力低于盒8組巖心流動阻力。

3)通過擬合效果來看,對于滲透率更低的巖心,改進模型的擬合效果更好,在數值模擬中,可將改進模型整合到模擬器中,對于油氣生產開發(fā)和油井工作制度的確定有良好效果。