宋延杰, 王鶴, 唐曉敏, 鄧鑫, 劉玥

(1.東北石油大學(xué)地球科學(xué)學(xué)院, 黑龍江 大慶 163318; 2.東北石油大學(xué)非常規(guī)油氣成藏與開(kāi)發(fā)省部共建國(guó)家重點(diǎn)實(shí)驗(yàn)室培育基地, 黑龍江 大慶 163318; 3.中石化勝利石油工程有限公司測(cè)井公司, 山東 東營(yíng) 257096)

0 引 言

產(chǎn)水率作為評(píng)價(jià)儲(chǔ)層產(chǎn)液性質(zhì)的重要參數(shù),能直觀反映儲(chǔ)層產(chǎn)油產(chǎn)水的相對(duì)比例,相對(duì)滲透率是產(chǎn)水率模型的重要組成部分。建立合理、實(shí)用的相對(duì)滲透率模型是確定儲(chǔ)層產(chǎn)液性質(zhì)的關(guān)鍵。巖石相對(duì)滲透率真值主要是通過(guò)實(shí)驗(yàn)方法獲得。實(shí)驗(yàn)室測(cè)定巖石相對(duì)滲透率的常用方法有穩(wěn)態(tài)法和非穩(wěn)態(tài)法。Leverett、Caudle等[1-2]在實(shí)驗(yàn)中采用穩(wěn)定流法,將油水按不同比例注入巖樣中,記錄壓差及油水流量,利用達(dá)西定律計(jì)算出各相的相對(duì)滲透率。后人[3-6]對(duì)穩(wěn)態(tài)法實(shí)驗(yàn)進(jìn)行改進(jìn),在實(shí)驗(yàn)中考慮了末端效應(yīng)和滯后效應(yīng),實(shí)驗(yàn)結(jié)果更準(zhǔn)確。然而,傳統(tǒng)穩(wěn)態(tài)法測(cè)試時(shí)間長(zhǎng)、效率低,測(cè)試過(guò)程易受環(huán)境溫度影響,測(cè)試結(jié)果存在相對(duì)較大的誤差。針對(duì)穩(wěn)態(tài)法的缺點(diǎn),Brace等[7]提出非穩(wěn)態(tài)法,它是通過(guò)在測(cè)試巖樣入口端施加一定壓力脈沖并記錄該壓力脈沖在巖樣中的衰減數(shù)據(jù),結(jié)合相應(yīng)的理論公式計(jì)算相對(duì)滲透率。后人[8-11]針對(duì)低滲透率油藏油水滲流特征,考慮毛細(xì)管力、重力及啟動(dòng)壓力梯度的影響,提出了低滲透率油藏非穩(wěn)態(tài)油水相對(duì)滲透率的數(shù)學(xué)模型。利用非穩(wěn)態(tài)實(shí)驗(yàn)方法確定油水相對(duì)滲透率,實(shí)驗(yàn)時(shí)間較短,獲取的數(shù)據(jù)較為可靠,但數(shù)學(xué)處理難度較大,如果對(duì)實(shí)驗(yàn)數(shù)據(jù)采用不恰當(dāng)?shù)臄M合方法其處理結(jié)果會(huì)出現(xiàn)較大偏差。人們又提出了利用壓汞法確定相對(duì)滲透率[12-14]。壓汞法測(cè)量速度快,資料豐富且測(cè)點(diǎn)多,但在壓汞結(jié)束后巖樣被汞污染不能重復(fù)使用。

針對(duì)實(shí)驗(yàn)方法的不足,一些學(xué)者提出利用實(shí)驗(yàn)測(cè)量的相對(duì)滲透率真值與含水飽和度、束縛水飽和度等參數(shù)建立經(jīng)驗(yàn)關(guān)系,如Purcell模型、Burdine模型、Corey模型等[15-17],從而達(dá)到利用測(cè)井資料計(jì)算的含水飽和度、束縛水飽和度等參數(shù)求取巖石相對(duì)滲透率的目的。經(jīng)驗(yàn)公式中回歸參數(shù)一般對(duì)于某一區(qū)塊取一個(gè)平均值很難描述巖石孔隙結(jié)構(gòu)變化影響?？衫秒娮杪屎暧^參數(shù)間接反映儲(chǔ)層孔隙的微觀變化。近年來(lái)分析家們提出了利用電阻率確定儲(chǔ)層的油水相對(duì)滲透率方法。Kewen Li[18]基于純巖石孔隙水流與電流流動(dòng)相似性原理,將含水和含油氣純巖石電導(dǎo)率與巖石絕對(duì)滲透率和水相有效滲透率等效類比,根據(jù)水相相對(duì)滲透率定義和物理界限,得出水相相對(duì)滲透率與含水飽和度及電阻率增大系數(shù)之間的關(guān)系;再根據(jù)Brooks-Corey模型,得出油相相對(duì)滲透率關(guān)系式。Attia A M[19]依據(jù)純巖石孔隙水流與電流流動(dòng)相似性原理,定義水力地層因素,將純巖石孔隙水分為可動(dòng)水與束縛水,并將可動(dòng)水孔隙等效為n根毛細(xì)管組成,依據(jù)可動(dòng)水流動(dòng)等效曲折度近似等于其導(dǎo)電等效曲折度,建立了水相相對(duì)滲透率與含水飽和度及電阻率增大系數(shù)之間的關(guān)系;結(jié)合比面積定義,建立了油相相對(duì)滲透率與含水飽和度及電阻率增大系數(shù)之間的關(guān)系。上述基于純巖石建立的相對(duì)滲透率與電阻率之間關(guān)系沒(méi)有考慮泥質(zhì)對(duì)于巖石導(dǎo)電性以及滲流的影響,因此,不能用于泥質(zhì)巖石。

大慶G地區(qū)葡萄花油層泥質(zhì)體積分?jǐn)?shù)范圍為10%～35%、孔隙度范圍為9%～19%、滲透率范圍為0.03～10 mD*非法定計(jì)量單位,1 mD=9.87×10-4 μm2,下同、微孔隙體積百分含量范圍為20%～50%,儲(chǔ)層具有含泥量高、泥質(zhì)附加導(dǎo)電性強(qiáng)、孔隙結(jié)構(gòu)復(fù)雜、中低孔隙度低滲透率的特征[20]。建立該區(qū)葡萄花油層油水相對(duì)滲透率與電阻率之間關(guān)系必須考慮泥質(zhì)對(duì)于巖石導(dǎo)電性以及滲流的影響。本文針對(duì)葡萄花油層的儲(chǔ)層特征,引用“三水”概念,將泥質(zhì)巖石總孔隙水分成可動(dòng)水、微孔隙水和黏土水,依據(jù)水流與電流流動(dòng)相似性原理,將可動(dòng)水孔隙等效為n根毛細(xì)管組成,結(jié)合泊肅葉方程和達(dá)西定律,推導(dǎo)水相相對(duì)滲透率與含水飽和度及可動(dòng)水流動(dòng)等效曲折度之間的關(guān)系,根據(jù)可動(dòng)水流動(dòng)等效曲折度近似等于其導(dǎo)電等效曲折度假設(shè)以及三孔隙導(dǎo)電模型[21-25],建立泥質(zhì)巖石水相相對(duì)滲透率與含水飽和度及電阻率增大系數(shù)之間的關(guān)系;結(jié)合比面積定義,建立泥質(zhì)巖石油相相對(duì)滲透率與含水飽和度及電阻率增大系數(shù)之間的關(guān)系。利用壓汞及巖電實(shí)驗(yàn)數(shù)據(jù)以及試油結(jié)果,評(píng)價(jià)公式的正確性。

1 水流與電流流動(dòng)相似性原理與水力地層因素定義

對(duì)于完全含水純巖石,假定地層水均為可動(dòng)水,且?guī)r石骨架不導(dǎo)電,則巖石導(dǎo)電來(lái)自巖石孔隙中水導(dǎo)電,其電流會(huì)沿著水流的流動(dòng)路徑傳導(dǎo),即可近似認(rèn)為巖石中水流流動(dòng)與電流導(dǎo)電具有相同的路徑。

對(duì)于只含有可動(dòng)水的圓柱形含水純巖石(見(jiàn)圖1),應(yīng)用達(dá)西定律以及歐姆定律和電阻定律可得巖石中水流量和電量表達(dá)式為

(1)

式中,Q為單位時(shí)間通過(guò)巖石的流量,cm3/s;μw為水黏度,×10-3Pa·s;L為巖石長(zhǎng)度,cm;Δp為巖石兩端的壓力差,×105Pa;R為巖樣半徑,cm;K為巖石滲透率,μm2;I0為通過(guò)巖石的電流強(qiáng)度,A或C/s;Δu為巖石兩端的電壓,V;R0為含水巖石電阻率,Ω·m。

對(duì)于半徑為R的單根毛細(xì)管(見(jiàn)圖2),應(yīng)用泊肅葉方程以及歐姆定律和電阻定律可得單根毛細(xì)管中水流量和電量表達(dá)式

(2)

式中,Q′為單位時(shí)間通過(guò)單根毛細(xì)管的流量,cm3/s;I′為通過(guò)單根毛細(xì)管的電流強(qiáng)度,A或C/s;Rw為水電阻率,Ω·m。

對(duì)于完全含水純巖石,地層電導(dǎo)率因素表示為

(3)

由水流與電流流動(dòng)相似性可知,Δp～Δu;對(duì)于同一圓柱形含水純巖石,I0～Q;對(duì)于同一單根毛細(xì)管,I′～Q′。與地層電導(dǎo)率因素類比,定義水力地層因素[19]為

(4)

圖1 飽含水純巖石模型

圖2 單根毛細(xì)管模型

2 含水泥質(zhì)巖石水力地層因素

對(duì)于飽含水泥質(zhì)巖石,結(jié)合“三水”概念[21-25],將孔隙中的水分為可動(dòng)水、微孔隙水、黏土水,而巖石骨架部分由巖石顆粒和干黏土顆粒組成,其等效體積模型見(jiàn)圖3。設(shè)可動(dòng)水孔隙由n根半徑為r的毛細(xì)管組成,引入可動(dòng)水孔隙毛細(xì)管曲折度τ,令毛細(xì)管長(zhǎng)度為l,有l(wèi)=Lτ,則有

n(πr2)l=πR2φt(1-Swit-Swb)L

(5)

式中,φt為泥質(zhì)巖石總孔隙度;Swit為泥質(zhì)巖石總微孔隙水飽和度;Swb為泥質(zhì)巖石總黏土水飽和度。

由泊肅葉方程和達(dá)西定律知,通過(guò)n根毛細(xì)管組成的可動(dòng)水孔隙的流量為

(6)

由式(5)和式(6)可得水力地層因素Fh的表達(dá)式為

(7)

圖3 含水泥質(zhì)巖石等效模型

圖4 含油氣泥質(zhì)巖石等效模型

3 含油氣泥質(zhì)巖石水相相對(duì)滲透率與電阻率增大系數(shù)關(guān)系建立

對(duì)于含油氣泥質(zhì)巖石,前述可動(dòng)水孔隙部分被油氣取代,其等效模型見(jiàn)圖4。類似式(7),可得含油氣泥質(zhì)巖石的水相水力地層因素Fhew表達(dá)式為

(8)

式中,Qw為單位時(shí)間通過(guò)含油氣泥質(zhì)巖石的水流量,cm3/s;new為含油氣泥質(zhì)巖石中可動(dòng)水孔隙的等效毛細(xì)管個(gè)數(shù);τew為含油氣泥質(zhì)巖石中可動(dòng)水孔隙的等效毛細(xì)管曲折度;Swt為含油氣泥質(zhì)巖石總含水飽和度。

對(duì)于含油氣泥質(zhì)巖石,應(yīng)用達(dá)西定律以及歐姆定律和電阻定律可得巖石中水流量和電量表達(dá)式為

(9)

(10)

由相對(duì)滲透率定義,將式(7)和式(8)代入式(10),可得含油氣泥質(zhì)巖石水相的相對(duì)滲透率表達(dá)式為

(11)

(12)

假設(shè)流體流動(dòng)與電流流動(dòng)相似,則式(11)中反映含水巖石和含油氣巖石中可動(dòng)水流動(dòng)的等效曲折度(τ、τew)可用其導(dǎo)電等效曲折度(τc和τcew)代替。

根據(jù)三孔隙導(dǎo)電理論,對(duì)于完全含水泥質(zhì)巖石和含油氣泥質(zhì)巖石(見(jiàn)圖3、圖4),在只有可動(dòng)水導(dǎo)電時(shí),微孔隙水和黏土水均可當(dāng)成不導(dǎo)電的骨架,其等效模型見(jiàn)圖5和圖6,按照并聯(lián)導(dǎo)電理論可得,在只有可動(dòng)水導(dǎo)電時(shí)的完全含水泥質(zhì)巖石和含油氣泥質(zhì)巖石地層因素為

(13)

(14)

(15)

(16)

由τ≈τc、τew≈τcew,將式(15)、式(16)以及Ir=Rft/Rf0代入式(11),則

(17)

式中,Ir為只有可動(dòng)水導(dǎo)電時(shí)的含油氣巖石電阻率增大系數(shù)。對(duì)式(17)進(jìn)行改進(jìn),則通用水相相對(duì)滲透率表達(dá)式為

(18)

(1) 當(dāng)Swt=1.0時(shí),有Swf=1,Ir=1,Krw=1,則c0=1,式(18)正確。

(2) 當(dāng)Swt=Swi+Swb時(shí),有Swf=0,Ir→∞,Krw=0,則式(18)正確。

將c0=1.0代入式(18),得

(19)

根據(jù)三孔隙導(dǎo)電理論[21-25],完全含水和含油氣泥質(zhì)巖石在只有可動(dòng)水導(dǎo)電時(shí),微孔隙水和黏土水均可當(dāng)成不導(dǎo)電的骨架;按照并聯(lián)導(dǎo)電理論與地層因素定義可得在只有可動(dòng)水導(dǎo)電時(shí)地層電阻率增大系數(shù)為

(20)

式中,φi、φc分別為泥質(zhì)巖石微孔隙水和黏土水孔隙度;mi、mc分別為泥質(zhì)巖石微孔隙水和黏土水孔隙指數(shù);ai、ac分別為泥質(zhì)巖石微孔隙水和黏土水孔隙系數(shù)。

圖5 只有可動(dòng)水導(dǎo)電的含水泥質(zhì)巖石等效模型

圖6 只有可動(dòng)水導(dǎo)電的含油氣泥巖石等效模型

4 含油氣泥質(zhì)巖石油相相對(duì)滲透率與電阻率增大系數(shù)關(guān)系建立

類似水相相對(duì)滲透率的推導(dǎo),可得出油相的相對(duì)滲透率公式為

(21)

式中,nnw為含油氣泥質(zhì)巖石中油氣孔隙的等效毛細(xì)管個(gè)數(shù);τnw為含油氣泥質(zhì)巖石中油氣孔隙的等效毛細(xì)管曲折度。

根據(jù)比表面定義,可知含油氣巖石中單位表面積對(duì)應(yīng)的可動(dòng)流體孔隙體積等于其比表面的倒數(shù);同理,含油氣巖石中單位表面積對(duì)應(yīng)的可動(dòng)水孔隙體積和油氣孔隙體積分別等于其比表面的倒數(shù)[19]。對(duì)于含油氣巖石,其可動(dòng)流體孔隙體積等于可動(dòng)水孔隙體積和油氣孔隙體積之和,則有

(22)

r=rew+rnw

(23)

由式(5)、式(6)可得可動(dòng)流體孔隙的毛細(xì)管半徑表達(dá)式為

(24)

同理,可得可動(dòng)水孔隙和油氣孔隙的毛細(xì)管半徑表達(dá)式為

(25)

式中,Kew為含油氣泥質(zhì)巖石水相有效滲透率,μm2。

(26)

式中,Knw為含油氣泥質(zhì)巖石油相有效滲透率,μm2。

將式(24)、式(25)和式(26)代入式(23),得

(27)

將式(15)、式(16)、式(17)和式(27)代入式(21),得

(28)

(29)

對(duì)式(29)改進(jìn),則通用油相相對(duì)滲透率表達(dá)式為

(30)

(1) 當(dāng)Swt=Swit+Swb時(shí),有Swf=0,Ir→∞,Krnw=1,則c=1.0,式(30)正確。

(2) 當(dāng)Swt=1.0時(shí),有Swf=1,Ir=1,Krnw=0,則式(30)正確。

將c=1.0和c0=1.0代入式(30),得

(31)

5 水相和油相相對(duì)滲透率與電阻率增大系數(shù)關(guān)系實(shí)驗(yàn)驗(yàn)證

5.1 巖石物理實(shí)驗(yàn)

為了定量研究水相和油相相對(duì)滲透率與電阻率增大系數(shù)關(guān)系,在G地區(qū)葡萄花油層取心井段內(nèi)選取6塊巖心樣本,且?guī)r樣兼顧了巖性、物性和電性的要求。為了確保用同一巖樣完成孔滲、巖電、粒度以及壓汞實(shí)驗(yàn)測(cè)量,設(shè)計(jì)了巖樣的實(shí)驗(yàn)流程,實(shí)驗(yàn)步驟如下。

(1) 對(duì)巖樣進(jìn)行洗油、洗鹽、烘干等預(yù)處理。在常溫條件下,應(yīng)用氣測(cè)法測(cè)量巖樣孔隙度和滲透率。

(2) 制備飽含水巖樣(水礦化度16 000 mg/L)。測(cè)試溫度為75 ℃,壓力為40 MPa,采用油驅(qū)水方式測(cè)量不同含水飽和度的巖樣電阻率。

(3) 對(duì)巖樣重新進(jìn)行預(yù)處理,切一塊巖樣,對(duì)巖樣進(jìn)行激光粒度分析;對(duì)剩余巖樣進(jìn)行壓汞實(shí)驗(yàn),測(cè)定巖樣的毛細(xì)管壓力曲線和孔喉大小分布。

5.2 實(shí)驗(yàn)驗(yàn)證

利用英×52井、英×6井、英×41井葡萄花儲(chǔ)層6塊巖樣孔滲、黏土體積分?jǐn)?shù)、壓汞及巖電實(shí)驗(yàn)數(shù)據(jù),得到φc=φclVcl、φt=φe+φc、Swb=φc/φt、Swit=Swiφe/φt、φf(shuō)=φe(1-Swit)、φi=φe-φf(shuō)。其中,φcl為黏土孔隙度,取值為0.2;φe為巖心分析有效孔隙度;Swi為有效微孔隙水飽和度,對(duì)應(yīng)壓汞實(shí)驗(yàn)中孔隙半徑小于0.1 μm的含水飽和度。令ai=ac=af=1,利用巖電實(shí)驗(yàn)數(shù)據(jù),對(duì)三孔隙電阻率模型采用最優(yōu)化技術(shù)求解Rt-Swf非相關(guān)函數(shù),可得到每塊巖樣mf、mi、mc、n值(見(jiàn)表1,Rtc為計(jì)算巖樣電阻率)。利用泥質(zhì)巖樣壓汞實(shí)驗(yàn)數(shù)據(jù)根據(jù)Burdine模型獲得水相和油相相對(duì)滲透率與含水飽和度實(shí)驗(yàn)關(guān)系曲線[19]。利用同一泥質(zhì)巖樣的相滲和巖電實(shí)驗(yàn)數(shù)據(jù),對(duì)式(19)、式(31)、式(20)采用最優(yōu)化技術(shù)求解Krw-Ir-Swf和Krnw-Ir-Swf的非相關(guān)函數(shù),可得到每塊巖樣n1、n2、n3、n4、n5、n6值(見(jiàn)表2,Krwc和Krnwc為計(jì)算巖樣水相和油相相對(duì)滲透率)。

表1 泥質(zhì)巖石三孔隙電阻率模型參數(shù)優(yōu)化值

表2 泥質(zhì)巖石水相和油相相對(duì)滲透率與含水飽和度及電阻率增大系數(shù)關(guān)系式參數(shù)優(yōu)化值

圖7 泥質(zhì)巖樣計(jì)算油和水相對(duì)滲透率與測(cè)量油和水相對(duì)滲透率對(duì)比圖(符號(hào)點(diǎn)代表巖心測(cè)量數(shù)據(jù),曲線代表巖心擬合結(jié)果)

圖7給出了6塊泥質(zhì)巖樣利用優(yōu)化參數(shù)和建立的相滲模型計(jì)算巖石油相和水相相對(duì)滲透率值與巖心測(cè)量值對(duì)比圖。圖7中,計(jì)算巖樣油相和水相相對(duì)滲透率與巖樣測(cè)量值吻合很好,平均相對(duì)誤差為4.9%和3.7%,誤差很低,說(shuō)明本文建立的相滲模型可以描述大慶G地區(qū)的葡萄花油層高含泥儲(chǔ)層的巖石相滲規(guī)律。

6 結(jié) 論

(1) 高含泥巖石建立其油水相對(duì)滲透率與電阻率之間關(guān)系必須考慮泥質(zhì)對(duì)巖石導(dǎo)電性以及滲流的影響。通過(guò)引入“三水”概念將泥質(zhì)巖石總孔隙水分成可動(dòng)水、微孔隙水和黏土水,依據(jù)可動(dòng)流體孔隙中水流流動(dòng)路徑與電流流動(dòng)路徑相同原理,利用泊肅葉方程和達(dá)西定律推導(dǎo)的水相相對(duì)滲透率關(guān)系式與利用三孔隙導(dǎo)電模型建立的電阻增大系數(shù)關(guān)系式結(jié)合,建立泥質(zhì)巖石水相相對(duì)滲透率與含水飽和度及電阻率增大系數(shù)之間關(guān)系式;結(jié)合比面積定義和可動(dòng)流體孔隙各組分體積等量關(guān)系,由水相相對(duì)滲透率關(guān)系式可推導(dǎo)出泥質(zhì)巖石油相相對(duì)滲透率與含水飽和度及電阻率增大系數(shù)之間的關(guān)系式。

(2) 在驗(yàn)證建立的泥質(zhì)巖石水相相對(duì)滲透率和油相相對(duì)滲透率與含水飽和度及電阻率增大系數(shù)之間關(guān)系式正確性中,必須采用同一巖樣而不采用并行樣進(jìn)行巖石物理實(shí)驗(yàn)測(cè)量,以確保儲(chǔ)層孔隙結(jié)構(gòu)和泥質(zhì)含量在巖電和壓汞實(shí)驗(yàn)測(cè)量中的一致性。其實(shí)驗(yàn)測(cè)量順序?yàn)榭诐B、巖電、激光粒度(取一薄片)、壓汞(剩余巖樣)。

(3) 利用泥質(zhì)巖樣的壓汞實(shí)驗(yàn)數(shù)據(jù)根據(jù)Burdine模型可以獲得水相和油相相對(duì)滲透率與含水飽和度實(shí)驗(yàn)關(guān)系曲線,利用同一泥質(zhì)巖樣的巖電實(shí)驗(yàn)數(shù)據(jù)根據(jù)三孔隙導(dǎo)電模型可以獲得對(duì)應(yīng)某一含水飽和度的假定只有可動(dòng)流體孔隙存在的巖石的電阻率增大系數(shù)值,將2種實(shí)驗(yàn)數(shù)據(jù)結(jié)合可驗(yàn)證建立的泥質(zhì)巖石水相相對(duì)滲透率和油相相對(duì)滲透率與含水飽和度及電阻率增大系數(shù)之間關(guān)系式正確性,實(shí)驗(yàn)證明基于水流與電流流動(dòng)相似性的泥質(zhì)巖石油水相對(duì)滲透率與電阻率關(guān)系模型能夠準(zhǔn)確求取儲(chǔ)層相對(duì)滲透率,可用于高含泥儲(chǔ)層產(chǎn)水率測(cè)井解釋。

參考文獻(xiàn):

[1] Leverett M C, Lewis W B. Steady Flow of Gas-oil-water Mixtures Through Unconsolidated Sands [C]∥Transactions of the AIME, 1941, 142(1): 107-116.

[2] Caudle B H. Further Developments in the Laboratory Determination of Relative Permeability [J]. Journal of Petroleum Technology, 1951, 3(5): 145-150.

[3] Saraf D N, Fatt I. Three-phase Relative Permeability Measurement Using Nuclear Magnetic Resonance Technique for Estimating Fluid Saturation [J]. SPE Journal, 1967, 7(3): 235-242.

[4] Schneider F N, Owens W W. Sandstone and Carbonate Two and Three-phase Relative Permeability Characteristics [J]. SPE Journal, 1970, 10(1): 75-84.

[5] Oak M J, Baker L E, Thomas D C. Three-phase Relative Permeability of Berea Sandstone [J]. Journal of Petroleum Technology, 1990, 42(8): 1054-1061.

[6] 周顯民, 袁慶峰, 劉桂芳, 等. 利用穩(wěn)態(tài)法測(cè)定油、氣、水三相相對(duì)滲透率實(shí)驗(yàn)研究 [J]. 石油學(xué)報(bào), 1991, 12(1): 69-78.

[7] Brace W F, Walsh J B, Frangos W T. Permeability of Granite Under High Pressure [J]. Journal of Geophysical Research, 1968, 73(6): 2225-2236.

[8] 吳克柳, 李相方, 樊兆琪, 等. 低滲-特低滲油藏非穩(wěn)態(tài)油水相對(duì)滲透率計(jì)算模型 [J]. 中國(guó)石油大學(xué)學(xué)報(bào), 2013, 37(6): 76-81.

[9] 李克文, 沈平平, 秦同洛. 考慮毛管壓力時(shí)計(jì)算油水相對(duì)滲透率新方法 [J]. 力學(xué)與實(shí)踐, 1994, 16(2): 46-48.

[10] 鄧英爾, 王允誠(chéng), 劉慈群, 等. 低滲非達(dá)西滲流相對(duì)滲透率計(jì)算方法及特征 [J]. 西南石油學(xué)院學(xué)報(bào), 2000, 22(3): 34-36.

[11] Civan F, Donaldson E C. Relative Permeability from Unsteady-state Displacements with Capillarypresure Included [J]. Journal of SPE Formation Evaluation, 1989, 4(2): 189-193.

[12] 王尤富. 用進(jìn)汞和退汞毛管壓力資料計(jì)算相對(duì)滲透率 [J]. 油氣采收率技術(shù), 1996, 3(2): 40-44.

[13] 王金勛, 楊普華, 劉慶杰, 等. 應(yīng)用恒速壓汞實(shí)驗(yàn)數(shù)據(jù)計(jì)算相隨滲透率曲線 [J]. 石油大學(xué)學(xué)報(bào), 2003, 27(4): 66-70.

[14] 孫煥泉, 洪世鐸, 孫士孝. 應(yīng)用壓汞資料計(jì)算油水相對(duì)滲透率的試驗(yàn)研究 [J]. 石油大學(xué)學(xué)報(bào), 1992, 16(1): 41-46.

[15] Purcell W R. Capillary Pressures—Their Measurement Using Mercury and the Calculation of Permeability [J]. Journal of petroleum technology, 1949, 1(2): 39 -48.

[16] Burdine N I. Relative Permeability Calculations from Pore Size Pistribution Data [J]. Journal of petroleum, 1953, 5(3): 71-78.

[17] Corey A T. The Interrelation Between Gas and Oil Relative Permeabilities [J]. Prod Mon, 1954, 19(1): 38-41.

[18] Kewen Li. A Semianalytical Method to Calculate Relative Permeability from Resistivity Well Logs [J]. SPE-95575-MS, 2005.

[19] Attia A M. Relative Permeability and Wettability of Rocks Obtained from Their Capillary Pressure and Electrical Resistivity Measurements [J]. SPE-106642-MS, 2007.

[20] 宋延杰, 姜艷嬌, 宋楊, 等. 古龍南地區(qū)低阻油層膠結(jié)指數(shù)和飽和度指數(shù)影響因素實(shí)驗(yàn) [J]. 吉林大學(xué)學(xué)報(bào), 2014, 44(2): 704-714.

[21] 莫修文, 賀鐸華, 李周波, 等. 三水導(dǎo)電模型及其在低阻儲(chǔ)層解釋中的應(yīng)用 [J]. 長(zhǎng)春科技大學(xué)學(xué)報(bào), 2001, 31(1): 92-95.

[22] 張麗華, 潘保芝, 李周波, 等. 新三水導(dǎo)電模型及其在低孔低滲儲(chǔ)層評(píng)價(jià)中的應(yīng)用 [J]. 石油地球物理勘探, 2010, 45(3): 431-435.

[23] 宋延杰, 李懷彪, 徐廣田, 等. 印度尼西亞方程的改進(jìn)及其在大慶長(zhǎng)垣低孔滲儲(chǔ)層的應(yīng)用 [J]. 大慶石油學(xué)院學(xué)報(bào), 2012, 36(2): 14-21.

[24] 宋延杰, 陳濟(jì)強(qiáng), 唐曉敏, 等. 低孔滲泥質(zhì)砂巖三水孔隙結(jié)合導(dǎo)電模型及其在海上油田的應(yīng)用 [J]. 大慶石油學(xué)院學(xué)報(bào), 2010, 34(6): 106-110.

[25] 張麗華, 潘保芝, 李寧, 等. 基于三水模型的儲(chǔ)層分類方法評(píng)價(jià)低孔隙度低滲透率儲(chǔ)層 [J]. 測(cè)井技術(shù), 2011, 3(1): 31-35.