低滲透致密油藏可動(dòng)流體飽和度計(jì)算方法
——以板橋地區(qū)長6油層組致密油儲(chǔ)層為例

曹雷，孫衛(wèi)，盛軍 霍磊，陳強(qiáng)

大陸動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室(西北大學(xué))

西北大學(xué)地質(zhì)學(xué)系，陜西 西安710069

謝彩花

(中石油青海油田分公司采油三廠，青海 海西州 816400)

低滲透致密油藏可動(dòng)流體飽和度計(jì)算方法

——以板橋地區(qū)長6油層組致密油儲(chǔ)層為例

曹雷，孫衛(wèi)，盛軍 霍磊，陳強(qiáng)

大陸動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室(西北大學(xué))

西北大學(xué)地質(zhì)學(xué)系，陜西 西安710069

謝彩花

(中石油青海油田分公司采油三廠，青海 海西州 816400)

[摘要]針對(duì)核磁共振試驗(yàn)獲得低滲透致密油儲(chǔ)層可動(dòng)流體飽和度成本高、周期長、應(yīng)用受限等弊端，深入分析了低滲透致密油儲(chǔ)層可動(dòng)流體飽和度的影響因素。結(jié)合恒速壓汞、高壓壓汞等試驗(yàn)提出了計(jì)算致密油可動(dòng)流體飽和度的方法。結(jié)果表明：低滲透致密油儲(chǔ)層可動(dòng)流體飽和度與恒速壓汞總進(jìn)汞飽和度、黏土礦物絕對(duì)含量相關(guān)性較好，將計(jì)算值與實(shí)測(cè)值比較，2種方法獲得的可動(dòng)流體飽和度差值小于5%，計(jì)算誤差平均為3.76%。計(jì)算結(jié)果與實(shí)測(cè)結(jié)果較為符合，表明該方法適用可靠。計(jì)算研究區(qū)長6油層組致密油儲(chǔ)層可動(dòng)流體飽和度平均為37.6%，按照儲(chǔ)層可動(dòng)流體飽和度劃分標(biāo)準(zhǔn)，研究區(qū)以Ⅲ類儲(chǔ)層為主。

[關(guān)鍵詞]致密油；恒速壓汞；高壓壓汞；可動(dòng)流體飽和度；板橋地區(qū)

低滲透致密油藏通常是指常規(guī)空氣法測(cè)滲透率小于2mD或覆壓基質(zhì)滲透率小于0.2mD的砂巖、碳酸鹽巖和頁巖儲(chǔ)層，具有巖性致密，非均質(zhì)性強(qiáng)，油水關(guān)系復(fù)雜，地層壓力異常等特點(diǎn)[1]。

目前，對(duì)鄂爾多斯盆地中生界延長組低滲透致密砂巖儲(chǔ)層來說，除了常規(guī)的孔隙度、滲透率、微觀孔隙結(jié)構(gòu)參數(shù)研究之外，缺少相應(yīng)的可動(dòng)流體參數(shù)研究，且通過現(xiàn)有的孔隙度、滲透率資料無法有效地預(yù)測(cè)低滲透致密油藏可動(dòng)流體飽和度。筆者在前人研究的基礎(chǔ)之上，綜合應(yīng)用核磁共振、高壓壓汞以及恒速壓汞技術(shù)，在分析研究可動(dòng)流體飽和度差異影響因素基礎(chǔ)之上，最終建立了研究區(qū)致密砂巖儲(chǔ)層可動(dòng)流體飽和度計(jì)算方法。

1核磁共振試驗(yàn)及可動(dòng)流體飽和度測(cè)定

隨著國內(nèi)外非常規(guī)油氣資源的開發(fā)利用，核磁共振以其測(cè)量精度高、數(shù)據(jù)可靠，及對(duì)低滲透致密砂巖儲(chǔ)層微觀孔隙結(jié)構(gòu)及孔隙內(nèi)流體流動(dòng)性質(zhì)的準(zhǔn)確反映，已經(jīng)成為研究低滲透致密油藏特征的一項(xiàng)關(guān)鍵技術(shù)[2]。在對(duì)低滲透致密砂巖儲(chǔ)層可動(dòng)流體飽和度的研究中，通過對(duì)飽和水或飽和油的巖樣進(jìn)行核磁共振分析，得到的τ2(橫向弛豫時(shí)間)的長短主要取決于流體分子受到孔隙固體表面作用力的強(qiáng)弱。τ2長短主要受巖石孔隙(孔隙大小、孔隙形態(tài))、礦物(礦物成分、礦物表面性質(zhì))和孔隙內(nèi)流體性質(zhì)(流體類型、流體黏度)3種因素控制。當(dāng)儲(chǔ)層孔隙半徑足夠小時(shí)，巖石固體表面的毛細(xì)管力或黏滯力就會(huì)將孔隙內(nèi)流體束縛，使之處于不可動(dòng)狀態(tài)，相對(duì)應(yīng)孔隙流體的τ2就存在一個(gè)界限，當(dāng)孔隙流體的τ2小于τ2截止值時(shí)，流體稱之為束縛流體或不可動(dòng)流體，反之則為自由流體或可動(dòng)流體[3，4]。在確定τ2截止值的基礎(chǔ)上，通過對(duì)比離心前后樣品τ2譜分布，計(jì)算低滲透致密砂巖儲(chǔ)層可動(dòng)流體飽和度。

圖1　板橋地區(qū)長6油層組致密油儲(chǔ)層核磁共振τ2譜

該次核磁共振τ2的測(cè)量使用的是中國石油勘探開發(fā)研究院廊坊分院引進(jìn)的Magnet2000型核磁共振巖心分析儀，試驗(yàn)溫度恒溫20℃，對(duì)研究區(qū)長6油層組致密砂巖儲(chǔ)層的10塊巖樣進(jìn)行了核磁共振試驗(yàn)，對(duì)自旋回波衰減曲線進(jìn)行了τ2譜解譜分析，確定該次試驗(yàn)τ2截止值為13.895ms，并計(jì)算束縛水飽和度等關(guān)鍵參數(shù)，最終獲取可動(dòng)流體飽和度測(cè)試結(jié)果(見圖1和表1)。

表1　核磁共振試驗(yàn)可動(dòng)流體飽和度測(cè)試結(jié)果(τ2截止值13.895ms)

2影響可動(dòng)流體飽和度差異因素

可動(dòng)流體飽和度差異影響因素研究是建立研究區(qū)低滲透致密砂巖儲(chǔ)層計(jì)算方法的前提。同一口井不同深度、不同巖性可動(dòng)流體參數(shù)變化復(fù)雜。前人研究表明，儲(chǔ)層物性、黏土礦物含量及充填形式、次生孔隙發(fā)育程度、微裂縫發(fā)育程度等主要因素影響著可動(dòng)流體飽和度的參數(shù)變化[5]。

2.1儲(chǔ)層物性

通過樣品物性參數(shù)與可動(dòng)流體飽和度參數(shù)進(jìn)行相關(guān)性分析表明，研究區(qū)致密砂巖儲(chǔ)層孔隙度與可動(dòng)流體飽和度相關(guān)性較差，相關(guān)系數(shù)為0.0982(圖2(a))；而可動(dòng)流體飽和度與滲透率呈正相關(guān)，相關(guān)系數(shù)為0.4147(圖2(b))。表明儲(chǔ)層中不同大小孔喉分布比例對(duì)可動(dòng)流體飽和度有一定影響，但非主導(dǎo)因素。總體來說，儲(chǔ)層物性與可動(dòng)流體飽和度之間不存在必然因果聯(lián)系[6]。

圖2　可動(dòng)流體飽和度與孔隙度(a)、滲透率(b)關(guān)系

2.2黏土礦物含量及充填形式

板橋地區(qū)低滲透致密油藏長6油層組儲(chǔ)層黏土礦物類型多、分布廣，對(duì)可動(dòng)流體飽和度影響大(圖3(a)、(b))。主要包含伊利石、綠泥石、伊-蒙混層及少量高嶺石等(圖4)。其中伊利石平均體積分?jǐn)?shù)最高，達(dá)到71.22%，多呈彎片狀或絲縷狀垂直碎屑顆粒表面生長或呈搭橋狀充填于顆粒之間[7～11]。伊利石賦存于孔隙使孔隙空間切割分解，孔喉結(jié)構(gòu)發(fā)生變化，逐漸曲折迂回。綠泥石及伊-蒙混層平均體積分?jǐn)?shù)相對(duì)較低，分別為13.22%和10.35%。綠泥石呈典型孔隙襯墊式包繞在顆粒表面，使得可動(dòng)流體賦存空間顯著減少。伊-蒙混層含量對(duì)可動(dòng)流體飽和度的影響相對(duì)最弱。由于黏土礦物大量富集，一方面使得儲(chǔ)層比表面大大增加，流體與孔喉表面摩擦性增強(qiáng)，增加了孔隙與喉道間的水膜厚度，使得一部分可動(dòng)流體轉(zhuǎn)化為束縛流體，進(jìn)而進(jìn)一步增加了孔喉表面的束縛流體含量。另一方面，黏土礦物微孔的存在也使得賦存其中的流體由于毛細(xì)管力的束縛作用而變的不可動(dòng)從而成為束縛流體[12]。

圖3　黏土絕對(duì)體積分?jǐn)?shù)與總進(jìn)汞飽和度(a)、可動(dòng)流體飽和度(b)關(guān)系

圖4　黏土礦物掃描電鏡照片

2.3次生孔隙發(fā)育

長6油層組致密油儲(chǔ)層下有長7油層組富含有機(jī)質(zhì)的頁巖，有機(jī)質(zhì)生烴脫羧產(chǎn)生的有機(jī)酸和CO2使地層孔隙水酸性顯著增強(qiáng)，其中長石沿著解理粒緣及粒內(nèi)溶蝕形成的孔隙為研究區(qū)最發(fā)育的次生孔隙。研究區(qū)三角洲前緣水下分流河道微相中溶蝕作用發(fā)育明顯，長石、巖屑等礦物經(jīng)過溶蝕作用改造，形成了眾多的粒內(nèi)溶孔和鑄模孔(圖5)。半深-深湖濁積砂體中也發(fā)育著較多的巖屑、長石粒緣溶孔和粒內(nèi)溶孔。次生孔隙的發(fā)育不僅增加了孔隙之間連通程度，而且極大改善了低滲透致密砂巖儲(chǔ)層的滲流能力，相關(guān)可動(dòng)流體參數(shù)也隨之提高。

2.4微裂縫發(fā)育

微裂縫自身具有一定的可動(dòng)流體含量，低滲透致密砂巖儲(chǔ)層孔喉細(xì)小，微裂縫溝通了互不連通的孔隙，有效改善了儲(chǔ)層的滲流能力，孔隙內(nèi)可動(dòng)流體含量隨之增加[13]。

結(jié)合前人研究[14]，影響研究區(qū)低滲透致密砂巖儲(chǔ)層可動(dòng)流體飽和度的2個(gè)主要因素為儲(chǔ)層微觀孔隙結(jié)構(gòu)和黏土礦物含量及充填形式，但前人并未定量分析黏土礦物含量對(duì)可動(dòng)流體飽和度的影響。黏土礦物如伊利石的搭橋狀充填孔隙內(nèi)部從根本上影響了儲(chǔ)層的孔隙結(jié)構(gòu)，間接控制了可動(dòng)流體飽和度的高低。此外，黏土礦物比表面大以及其表面所帶電荷的強(qiáng)吸附性直接影響可動(dòng)流體飽和度。因此，筆者在前人研究基礎(chǔ)之上，重點(diǎn)定量分析了黏土礦物含量及充填形式這一影響因素。

3可動(dòng)流體飽和度計(jì)算方法及結(jié)果分析

通過高壓壓汞和恒速壓汞試驗(yàn)，不僅能準(zhǔn)確反映低滲透致密砂巖儲(chǔ)層的微觀孔隙結(jié)構(gòu)，而且利用相關(guān)τ2譜分布曲線也能建立壓汞曲線轉(zhuǎn)化模型[15]。因此，可動(dòng)流體飽和度與壓汞曲線之間存在著密切聯(lián)系。

鄂爾多斯盆地中生界延長組低滲透致密砂巖儲(chǔ)層研究實(shí)例表明，核磁共振可動(dòng)流體飽和度與主要壓汞參數(shù)之間存在較好的相關(guān)關(guān)系[16]。但極少有人進(jìn)行低滲透致密油藏可動(dòng)流體飽和度計(jì)算方法研究，筆者在對(duì)研究區(qū)10塊樣品高壓壓汞和恒速壓汞試驗(yàn)分析結(jié)果基礎(chǔ)之上，充分考慮黏土礦物含量及充填形式影響，建立了適合研究區(qū)低滲透致密砂巖儲(chǔ)層的可動(dòng)流體飽和度計(jì)算公式[17～21]。

試驗(yàn)結(jié)果表明：分選系數(shù)、中值半徑、有效喉道半徑、主流喉道半徑、最大連通喉道半徑、單位體積有效孔隙體積、總進(jìn)汞飽和度與可動(dòng)流體飽和度呈正相關(guān)關(guān)系，相關(guān)系數(shù)依次為0.6256、0.8202、0.6446、0.7624、0.7725、0.7325、0.8556(圖6(a)～(g))；孔喉半徑比、中值壓力、閾壓與可動(dòng)流體飽和度呈負(fù)相關(guān)關(guān)系，相關(guān)系數(shù)分別是0.7209、0.9319、0.9193(圖6(h)～(j))，由此可見，致密油藏可動(dòng)流體飽和度在一定程度上受儲(chǔ)層微觀孔隙結(jié)構(gòu)的控制作用[21～27]。

由于可動(dòng)流體飽和度與中值壓力、閾壓、黏土礦物絕對(duì)體積分?jǐn)?shù)以及總進(jìn)汞飽和度的相關(guān)性較好，因此可以根據(jù)上述參數(shù)的擬合公式來計(jì)算可動(dòng)流體飽和度。考慮到研究區(qū)致密油儲(chǔ)層中只有少量的樣品進(jìn)汞飽和度超過50%，故無法利用中值壓力來計(jì)算可動(dòng)流體飽和度；且由壓汞試驗(yàn)(恒速壓汞和高壓壓汞)獲得的閾壓數(shù)據(jù)易受試驗(yàn)樣品橫縱切面形態(tài)影響，無法準(zhǔn)確得出可動(dòng)流體飽和度數(shù)據(jù)[28]。研究區(qū)黏土礦物富集，其中黏土礦物絕對(duì)體積分?jǐn)?shù)與總進(jìn)汞飽和度和可動(dòng)流體飽和度的相關(guān)系數(shù)分別為0.7313和0.8977，在一定程度上可以反映可動(dòng)流體飽和度飽和度的高低，故在誤差平方和最小的前提下，采用二元線性回歸模型建立研究區(qū)致密油儲(chǔ)層可動(dòng)流體飽和度與黏土礦物絕對(duì)體積分?jǐn)?shù)、總進(jìn)汞飽和度之間的關(guān)系式：

圖6　壓汞參數(shù)與可動(dòng)流體飽和度關(guān)系圖

Sm=0.4402×SHg-5.605×φ(clay)+53.98

式中，Sm為可動(dòng)流體飽和度計(jì)算值，%；SHg為恒速壓汞試驗(yàn)總進(jìn)汞飽和度，%; φ(clay)為黏土礦物絕對(duì)體積分?jǐn)?shù)，%。

在實(shí)際計(jì)算過程中，鑒于恒速壓汞試驗(yàn)存在成本較高、周期長等因素，因此可以考慮使用方便易得的高壓壓汞參數(shù)來替代恒速壓汞試驗(yàn)數(shù)據(jù)。但是由于高壓壓汞試驗(yàn)的滯后效應(yīng)(進(jìn)汞壓力通常比恒速壓汞高0.5～1MPa)，所以在建立計(jì)算公式時(shí)應(yīng)充分考慮該重要因素[29]。該次恒速壓汞試驗(yàn)的最大進(jìn)汞壓力為6.2MPa,故恒速壓汞試驗(yàn)總進(jìn)汞飽和度可由高壓壓汞試驗(yàn)6.8MPa時(shí)的進(jìn)汞飽和度代替。

板橋地區(qū)長6油層組致密油儲(chǔ)層可動(dòng)流體飽和度計(jì)算公式為：

Sm=0.4402×SHg6.8-5.605×φ(clay)+53.98

式中， SHg6.8為高壓壓汞試驗(yàn)在進(jìn)汞壓力為6.8MPa時(shí)的進(jìn)汞飽和度，%。

將計(jì)算得到的可動(dòng)流體飽和度與實(shí)測(cè)可動(dòng)流體飽和度值進(jìn)行比較(見表2)，1號(hào)樣品出現(xiàn)較大誤差，誤差原因可能是由于試驗(yàn)誤差導(dǎo)致。其余9塊樣品計(jì)算值與實(shí)測(cè)值都較為接近，2種方法獲得的可動(dòng)流體飽和度差值小于5%，除去1號(hào)異常樣品其余樣品計(jì)算誤差平均為3.76%。計(jì)算結(jié)果與實(shí)測(cè)結(jié)果較為符合，表明該方法適用可靠。

利用上述建立的研究區(qū)長6油層組致密油儲(chǔ)層可動(dòng)流體飽和度計(jì)算公式，將板橋地區(qū)長6油層組儲(chǔ)層50塊樣品可動(dòng)流體飽和度進(jìn)行計(jì)算，結(jié)果表明:可動(dòng)流體飽和度最大為65.3%，最小為20.5%，平均為37.6%。

表2　可動(dòng)流體飽和度計(jì)算結(jié)果及試驗(yàn)誤差

通過對(duì)比分析計(jì)算結(jié)果與樣品實(shí)測(cè)可動(dòng)流體飽和度與孔隙度、滲透率、中值壓力、排驅(qū)壓力等參數(shù)的相關(guān)關(guān)系(圖7(a)～(d))，驗(yàn)證了計(jì)算方法的合理性及可靠性。根據(jù)核磁共振可動(dòng)流體飽和度劃分標(biāo)準(zhǔn)(表3)，可將板橋地區(qū)長6油層組低滲透致密砂巖儲(chǔ)層劃分為4類， 其中Ⅰ類儲(chǔ)層占2%，Ⅱ類儲(chǔ)層占21.6%，Ⅲ類儲(chǔ)層占50.2%,Ⅳ類儲(chǔ)層為26.2%(圖8)。研究區(qū)主要為Ⅲ類儲(chǔ)層和Ⅳ類儲(chǔ)層，Ⅱ類儲(chǔ)層以及Ⅰ類儲(chǔ)層相對(duì)較少。

4結(jié)論

1) 研究區(qū)低滲透致密油儲(chǔ)層可動(dòng)流體飽和度是進(jìn)行致密油儲(chǔ)層精細(xì)評(píng)價(jià)的重要參數(shù)之一，主要受黏土礦物含量及充填形式和儲(chǔ)層微觀孔隙結(jié)構(gòu)2種因素的影響，其中黏土礦物含量及充填形式的控制作用尤為顯著。

圖7　可動(dòng)流體飽和度計(jì)算值與各參數(shù)的關(guān)系

2) 建立研究區(qū)低滲透致密油儲(chǔ)層可動(dòng)流體飽和度計(jì)算模型，將計(jì)算值與實(shí)測(cè)值比較，2種方法獲得的可動(dòng)流體飽和度差值小于5%，計(jì)算誤差平均值為3.76%。計(jì)算結(jié)果與實(shí)測(cè)結(jié)果較為符合，表明該方法適用可靠。

圖8　　　　板橋地區(qū)長6油層組致密油儲(chǔ)層　　　計(jì)算可動(dòng)流體飽和度分布

可動(dòng)流體飽和度/%儲(chǔ)層分類>65Ⅰ50～65Ⅱ35～50Ⅲ20～35Ⅳ<20Ⅴ

3) 利用建立的研究區(qū)低滲透致密油儲(chǔ)層可動(dòng)流體飽和度計(jì)算模型，計(jì)算研究區(qū)長6油層組致密油儲(chǔ)層可動(dòng)流體飽和度最大為65.3%，最小為20.5%，平均為37.6%。其中儲(chǔ)層類型以Ⅲ類為主。

[參考文獻(xiàn)]

[1]魏海峰, 凡哲元, 袁向春. 致密油藏開發(fā)技術(shù)研究進(jìn)展[J]. 油氣地質(zhì)與采收率, 2013，20(2): 62～66.

[2]王為民, 孫佃慶. 核磁共振測(cè)井基礎(chǔ)實(shí)驗(yàn)研究[J]. 測(cè)井技術(shù), 1997, 21(6): 385～392.

[3]高輝, 孫衛(wèi), 田育紅, 等. 核磁共振技術(shù)在特低滲砂巖微觀孔隙結(jié)構(gòu)評(píng)價(jià)中的應(yīng)用[J]. 地球物理學(xué)進(jìn)展, 2011, 26(1): 294～299.

[4]陳冬霞, 龐雄奇, 姜振學(xué), 等. 利用核磁共振物理模擬實(shí)驗(yàn)研究巖性油氣藏成藏機(jī)理[J]. 地質(zhì)學(xué)報(bào), 2006, 80(3): 432～438.

[5]王瑞飛, 陳明強(qiáng). 特低滲透砂巖儲(chǔ)層可動(dòng)流體賦存特征及影響因素 [J]. 石油學(xué)報(bào), 2008, 29(4)：558～561.

[6]葉禮友, 高樹生, 熊偉, 等. 可動(dòng)水飽和度作為低滲砂巖氣藏儲(chǔ)層評(píng)價(jià)參數(shù)的論證[J]. 石油天然氣學(xué)報(bào)(江漢石油學(xué)院學(xué)報(bào)), 2011, 33(1): 57～59.

[7]Paxton S T, Szzbo J O. Ajdukiewicz J M, et al. Construction of an intergranular volume compaction curve for evaluating and predicting compaction and porosity loss in rigid-grain sandstone reservoirs[J]. AAPG Bulletin, 2002, 86(12):2047～2067.

[8]王金鵬, 彭仕宓, 趙艷杰,等. 鄂爾多斯盆地合水地區(qū)長6-8段儲(chǔ)層成巖作用及孔隙演化[J]. 石油天然氣學(xué)報(bào)(江漢石油學(xué)院學(xué)報(bào)), 2008, 30(2): 170～174.

[9]Ehrenberg S N. Assessing the relative importance of compaction processes and cementation to reduction of porosity in sandstones: discussion; compaction and porosity evolution of Pliocene sandstones, Ventura Basin, California: discussion[J]. AAPG Bulletin, 1989, 73(10):1274～1276.

[10] Paxton S T, Szzbo J O. Ajdukiewicz J M, et al. Construction of an intergranular volume compaction curve for evaluating and predicting compaction and porosity loss in rigid-grain sandstone reservoirs[J]. AAPG Bulletin, 2002, 86(12):2047～2067.

[11]Berger A, Gier S, Krois P. Porosity-preserving chlorite cements in shallow-marine volcaniclastic sandstones: Evidence from Cretaceous sandstones of the Sawan gas field, Pakistan[J]. AAPG Bulletin, 2009, 93(5):595～615.

[12]高輝, 孫衛(wèi). 特低滲透砂巖儲(chǔ)層可動(dòng)流體變化特征與差異性成因——以鄂爾多斯盆地延長組為例[J]. 地質(zhì)學(xué)報(bào), 2010, 84(8): 1223～1230.

[13]曾聯(lián)波. 低滲透砂巖油氣儲(chǔ)層裂縫及其滲流特征[J]. 地質(zhì)科學(xué), 2004, 39(1): 11～17.

[14]熊偉, 雷群, 劉先貴, 等. 低滲透油藏?cái)M啟動(dòng)壓力梯度[J]. 石油勘探與開發(fā), 2009, 36(2)：232～237.

[15]何雨丹, 毛志強(qiáng), 肖立志, 等. 利用核磁共振 T2分布構(gòu)造毛管壓力曲線的新方法[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2005, 35(2): 177～181.

[16]師調(diào)調(diào), 孫衛(wèi), 何生平. 低滲透儲(chǔ)層微觀孔隙結(jié)構(gòu)與可動(dòng)流體飽和度關(guān)系研究[J]. 地質(zhì)科技情報(bào), 2012, 31(4): 81～83.

[17]Luthi, Stefan M. Geological Well Logs[M].Berlin :Springer Berlin Heidelberg, 2001.

[18] Guan H, Brougham A D, Sorbieetal K S . Wettability effects in a sandstone reservoir and outcrop cores from NMR relaxation time distributions [J]. Journal of Petroleum Science and Engineering, 2002, 34:35～54.

[19]Chen J, Hirasaki G J, Flaum M. NMR wettability indices: Effect of OBM on wettability and NMR responses [J]. Journal of Petroleum Science and Engineering, 2006, 52:161～171.

[20] Al-Mahrooqi S H, Crattoni C A, Moss A K, et al. An investigation of the effect of wettability on NMR characteristics of sandstone rock and fiuld systems [J]. Journal of Petroleum Science and Engineering, 2003, 39:389～398.

[21] Mahrooqi S H, Crattoni C A , Muggeridge A H , et al. Pore-scale modeling of NMR relaxation for the characteristics of wettability [J]. Journal of Petroleum Science and Engineering, 2006, 52:172～186.

[22]Guan C J, Ying Y L, Zhang M. Evaluation of gas wettability and its effects on fluid distribution and fluid flow in porous media[J]. Petroleum Science，2013, 10(4): 512～527.

[23]Zhang V, Stallmatch C H, Skjetne T,et al. NMR response of non-reservoir fluids in sandstone and chalk [J]. Magnetic Resonance Imaging, 2001, 19:543～545.

[24]Jim D K. The influence factors of NMR logging porosity in complex fluid reservoir[J]. Science in China, 2008, 51(2 Supplement):212～217.

[25]Hirasaki G J, Lo S W, Zhang Y. NMR properties of petroleum reservoir fluids.[J]. Magnetic Resonance Imaging, 2003, 21(3-4):269～277.

[26]Davies S. Quantification of oil and water in preserved reservoir rock by NMR spectroscopy and imaging[J]. Magnetic Resonance Imaging, 1994, 12(2):349～353.

[27]Cheng Z, Luo S, Zhiwen D U, et al. The Method to Calculate Tight Sandstone Reservoir Permeability Using Pore Throat Characteristic Parameters[J]. Well Logging Technology, 2014,38(2):185～189.

[28]喻建, 楊孝, 李斌, 等. 致密油儲(chǔ)層可動(dòng)流體飽和度計(jì)算方法[J]. 石油實(shí)驗(yàn)地質(zhì), 2014, 36(6):767～779.

[29]闕洪培，雷卞軍. NMR 束縛流體體積確定方法及其應(yīng)用條件[J]. 西南石油大學(xué)學(xué)報(bào) (自然科學(xué)版), 2002, 24(2): 5～8.

[30]雷啟鴻, 孫華嶺, 李晶, 等. 長慶低滲透儲(chǔ)層可動(dòng)流體飽和度評(píng)價(jià)[J]. 低滲透油氣田, 2006, 11(1): 94～97.

[編輯]鄧?yán)?/p>

[收稿日期]2015-10-10

[作者簡介]曹雷(1991-)，男，碩士生，現(xiàn)主要從事低滲、特低滲儲(chǔ)層評(píng)價(jià)方面的學(xué)習(xí)與研究工作，287384095@qq.com。

[中圖分類號(hào)]TE122.2;TE311

[文獻(xiàn)標(biāo)志碼]A

[文章編號(hào)]1673-1409(2016)20-0001-08

[引著格式]曹雷，孫衛(wèi)，盛軍，等.低滲透致密油藏可動(dòng)流體飽和度計(jì)算方法[J].長江大學(xué)學(xué)報(bào)(自科版),2016,13(20):1～8.