周 鴻
(油氣藏地質(zhì)及開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室·西南石油大學(xué),四川 成都 610500)
氣藏欠平衡鉆井過(guò)程產(chǎn)生的鉆井液自吸效應(yīng)會(huì)使近井地帶產(chǎn)生水鎖現(xiàn)象,儲(chǔ)層滲透率會(huì)急劇降低。目前,國(guó)外對(duì)自吸的研究多是考慮毛管力的作用[1-3],雖然得出了部分結(jié)論,但都是考慮在自然狀態(tài)下毛管垂直吸水的情況[4]。而國(guó)內(nèi)對(duì)欠平衡鉆井狀態(tài)下自吸的研究也沒(méi)有定量的分析[5],討論毛管力作用十分困難,因此筆者避開(kāi)此討論,采用氣體滲流理論,分析氣藏欠平衡鉆井自吸穩(wěn)態(tài)時(shí)壓力和滲流速度在巖心內(nèi)的分布及其對(duì)滲透率的影響。
此理論假設(shè)條件為:① 氣體為單相達(dá)西滲流;② 巖心為均勻介質(zhì);③ 滲流過(guò)程等溫;④ 滲透率為常數(shù);⑤ 平行流動(dòng)中重力忽略不計(jì)。
氣體單相穩(wěn)定滲流微分方程為:
由于實(shí)驗(yàn)室為低壓等溫狀態(tài),所以還可假設(shè)μg與Z均為常數(shù);自吸實(shí)驗(yàn)可以看做直線單向平行滲流,因此式(1)可化為:
在滿足邊界條件下式(2)解得:
在獲得實(shí)驗(yàn)室狀態(tài)下氣體單相平行流動(dòng)的壓力與流速分布函數(shù)后,考慮自吸過(guò)程。毛管自吸過(guò)程是1個(gè)類似活塞運(yùn)動(dòng)的過(guò)程[6],所以可以將自吸端與未自吸端考慮為兩種均勻介質(zhì)。隨著時(shí)間的延長(zhǎng),自吸過(guò)程最后將趨于動(dòng)態(tài)平衡[7]。
自吸穩(wěn)態(tài)時(shí)壓力分布方程為:
邊界條件為:
在等溫穩(wěn)定狀態(tài)下,由玻意耳定律推得:
將式(12)代入式(11)得出:
因此,只要獲得Kg1、Kg2和 p0中任意1個(gè)參數(shù)就可以通過(guò)式(7)和式(8)以及邊界條件下式(10)和式(13)得出自吸穩(wěn)態(tài)時(shí)壓力和流速在巖心中的分布及對(duì)應(yīng)的滲透率。
蘇里格氣田位于鄂爾多斯盆地依陜斜坡西北側(cè),其東區(qū)為一大面積分布的砂巖巖性氣藏。氣田地質(zhì)條件非常復(fù)雜,是一個(gè)低孔、低滲、低壓、低豐度的砂巖氣藏。選用4塊蘇里格氣藏東區(qū)低滲砂巖巖心進(jìn)行實(shí)驗(yàn),干燥巖心氣測(cè)滲透率為0.01~0.2 mD。實(shí)驗(yàn)步驟為:① 將巖心烘干,測(cè)量L和r后將其放入巖心夾持器中;② 打開(kāi)氮?dú)馄亢烷y門,讓氮?dú)庠谝欢▔毫ο峦ㄟ^(guò)巖心,同時(shí)保持圍壓穩(wěn)定;③在氣體流速穩(wěn)定后讀出壓差Δp和流量v,然后向模擬井筒水管中加入一定量水,水面高于巖心橫截面最高點(diǎn);④ 自吸開(kāi)始,定時(shí)測(cè)量氣體流量,直到流量穩(wěn)定不變時(shí),自吸達(dá)到動(dòng)態(tài)平衡狀態(tài);⑤ 測(cè)量出平衡后的p1、p2和v2,并剖開(kāi)巖心用圖像掃描儀測(cè)得m值;⑥ 烘干后再將巖心注入蒸餾水,然后在不同初始含水飽和度下重復(fù)步驟②~⑤,測(cè)量出不同飽和度下巖心自吸穩(wěn)態(tài)平衡后所對(duì)應(yīng)的p1、p2、v2和m,見(jiàn)圖1。
從實(shí)驗(yàn)可以看出,當(dāng)自吸開(kāi)始時(shí),滲流速度急劇下降,隨后下降變緩直至最后趨于平衡狀態(tài)。在不考慮滑脫效應(yīng)的情況下,通過(guò)Δp和v可得出Kg1,進(jìn)而通過(guò)式(7)、(8)、(10)和(13)得出自吸穩(wěn)態(tài)時(shí)壓力和流速在巖心的分布情況,并通過(guò)式(9)得出水浸部分的氣體滲透率。
在初始含水飽和度分別為20%、30%和50%的條件下,通過(guò)實(shí)驗(yàn)得出自吸穩(wěn)態(tài)下的壓力和流速的分布情況,即含水飽和度越高,壓力和滲流速度較自吸前的變化越小,見(jiàn)圖2和圖3。
圖1 巖心實(shí)驗(yàn)步驟圖
圖2 不同初始含水飽和度下自吸穩(wěn)態(tài)壓力分布圖
圖3 不同初始含水飽和度下自吸穩(wěn)態(tài)氣體流速分布圖
1)實(shí)驗(yàn)證明,自吸開(kāi)始時(shí)滲透率急劇下降,但在外部條件不變的情況下,只要自吸時(shí)間足夠長(zhǎng),最終一定能夠達(dá)到動(dòng)態(tài)平衡,形成穩(wěn)定狀態(tài)。
2)在穩(wěn)態(tài)情況下可以通過(guò)式(7)、式(8)及其邊界條件得出巖心中壓力和滲流速度的分布情況。
3)在不考慮滑脫效應(yīng)的情況下可以通過(guò)式(9)得出水浸部分的氣體滲透率。
4)實(shí)驗(yàn)證明了巖心初始含水飽和度越高,壓力和流速分布越接近于自吸前壓力和流速分布的情況。
符號(hào)說(shuō)明
Kg為氣測(cè)滲透率,D;μg為氣體黏度,mPa·s;Z為氣體偏差因子;p為氣體單相平行流動(dòng)壓力,MPa;x為距滲流出口的端面距離,cm;L為巖心長(zhǎng)度,cm;Δp為壓差,MPa;p1為進(jìn)口端壓力,MPa;p2為出口端壓力,MPa;v為滲流速度,cm/s;p0為自吸穩(wěn)態(tài)后干濕界面的對(duì)應(yīng)壓力,MPa;v1為進(jìn)口滲流速度,cm/s;v2為出口滲流速度,cm/s;v0為p0所對(duì)應(yīng)的滲流速度,cm/s;l為自吸穩(wěn)態(tài)后巖心干燥部分側(cè)面長(zhǎng)度,cm;m為自吸穩(wěn)態(tài)后巖心潤(rùn)濕部分側(cè)面長(zhǎng)度,cm;Kg1為l段氣體滲透率,D;Kg2為m段氣體滲透率,D;Q為氣體流量,cm3/s;pc為標(biāo)準(zhǔn)狀態(tài)下的大氣壓,MPa;Qc為標(biāo)準(zhǔn)狀態(tài)下的氣體流量,cm3/s;r為巖心截面半徑,cm。
[1]L.L.Handy.Determination of Effective Capillary Pres?sures for Porous Media from Imbibition Data[R].KansasCity:Joint SPE-AIChE Meeting,1959.
[2]C.C.Mattax,J.R Kyte.Imbibition Oil Recovery from Frac?tured,Water-Drive Reservoir[J].SPE Journal, 1962,2(2):177-184.
[3]S.Ma,N.R.Morrow.Generalized Scaling Of Spontaneous Imbibition Data For Strongly Water-Wet Systems[R].Re?gina:The Sixth Petroleum Conference Of The South Sas?katchewan Section,1995.
[4]D.Prey,E.Lefebvre.Gravity and Capillarity Effects on Im?bibition in Porous Media[J].SPE Journal,1978,18(3):195-206.
[5]游利軍,康毅力.油氣儲(chǔ)層巖石毛細(xì)管自吸研究進(jìn)展[J].西南石油大學(xué)學(xué)報(bào):自然科學(xué)版,2009,31(4):112-116.
[6]S.Akin,A.R.Kovscek.Imbibition Studies of Low-Permea?bility Porous Media[R].Anchorage:1999 Western Region Meeting,1999.
[7]KEVEN L,ROLAND N.Characterization of Spontaneous Water Imbibition into Gas-Saturated Rocks[R].Long Beach:the 2000 SPE/AAPG Western Regional Meet?ing,2000.