基于NMR技術研究邊界條件對致密砂巖油水、氣水系統滲吸的影響

程志林 ，王慶 *，寧正福 ，李銘琪 ，齊榮榮 ，余雄飛

1 中國石油大學(北京)油氣資源與探測國家重點實驗室，北京 102249

2 中國石油大學(北京)石油工程學院，北京 102249

0 引言

伴隨著日益嚴峻的能源危機以及非常規油氣資源勘探的不斷深入，致密油氣藏已成為重要的替代能源[1-2]。水平井壓裂技術已經廣泛地應用于非常規油氣資源的開發當中[3]，大量的壓裂液被注入地層中，研究發現壓裂液的返排率低(少于30%)[4]，壓裂液的滯留造成了儲層損害，流體依靠自發滲吸進入儲層是造成這一現象的主要原因[5]。另一方面，當壓裂液進入儲層后，巖石孔隙中的油氣能夠被置換出來，因此自發滲吸對致密儲層而言，也是一種重要的提高油氣采收率方法[6-7]。

自發滲吸是指在多孔介質中潤濕相流體依靠毛管力作用置換出非潤濕相流體的過程[8-9]。根據滲吸過程中濕相流體和非濕相流體流動方向，可以將滲吸過程分為正向滲吸和逆向滲吸。由于儲層經壓裂改造后巖石被裂縫中的壓裂液包圍，這種情況下流體一般發生逆向滲吸，而正向滲吸發生幾率較小，因此過去的研究重點主要集中在多孔介質逆向滲吸[10-11]。

自發滲吸取決于多孔介質的性質、流體類型以及流固之間相互作用。影響因素包括巖石孔隙度、滲透率、孔隙結構、巖樣的形狀和大小、邊界條件、流體黏度、初始含水飽和度、潤濕性條件和相滲等，滲吸的表征是一個復雜的過程[12]。Zhou等[13]研究了在不同潤濕狀態和不同初始含水飽和度條件下，Berea砂巖自發滲吸和水驅實驗中采收率的變化特征。發現隨著老化時間的增加，由于殘余油飽和度的降低，水驅采收率得到了提高。而滲吸采油大幅度的下降，最終采收率總體上呈現下降的趨勢。Li和Horne等[14]基于滲吸前緣活塞式推進假設，同時考慮了重力作用和流體相滲提出了一個標度模型，并用實驗驗證了模型的合理性。Standnes[15]通過實驗系統地研究了正向和逆向滲吸條件下，邊界條件和樣品形狀對滲吸的影響。邊界條件包括全部開啟(AFO)、兩端開啟(TEO)、兩端封閉(TEC)和一端開啟(OEO)邊界條件等。研究結果表明：標度模型無法應用于不規則形狀的巖石，巖石與濕相流體接觸面積對滲吸速度影響較大而對采收率影響不大。由于只有在樣品較高的情況下，重力對滲吸的影響才能得以體現，室內滲吸實驗實施較為困難，Mirzaei-Paiaman[16]通過數值模擬研究了重力對裂縫性儲層巖石逆向滲吸的影響，并基于檢視性分析提出了一個新的標度模型。Meng等[17-18]利用填充玻璃珠和石英砂結構分別研究了濕相和非濕相黏度對滲吸采收率的影響，發現在相同邊界條件下，濕相黏度對于滲吸采收率影響較小；而隨著非濕相流體黏度的增大，油相流體在均質性較差的石英砂結構滯留嚴重，滲吸采收率逐漸降低；而均質性較好的玻璃珠結構，非濕相黏度對采收率影響不大。

過去的研究對于常規儲層巖石強水濕滲吸特征認識的較為全面，提出的標度模型得到了廣泛的應用。而對巖石在混合潤濕及油濕條件下流體滲吸特征目前的認識還遠遠不夠，非常規油氣儲層以頁巖為例，由于其骨架既有機質也有無機質礦物，有機質部分表現出疏水性而無機質為親水性，因而頁巖往往呈現出混合潤濕的特點[19]，傳統強水濕滲吸理論將不再適用于這類儲層。Dutta等[20]基于CT技術對比了致密砂巖和中高滲砂巖氣水滲吸特點，發現致密砂巖極低的滲透率會對滲吸產生抑制作用，其孔隙結構的非均質性特征造成了不均勻氣水驅替前緣，滲吸進入巖樣的水包括黏土礦物中的束縛水和孔隙空間的自由水。Lai等[21]研究了不同種類巖石在不同滲吸流體介質條件下的滲吸規律，并分析了巖石孔隙結構、潤濕性及溫度對滲吸采收率的影響。Lin等[22]通過NMR(核磁共振)實驗研究了頁巖氣滲吸過程中，不同位置含水飽和度隨時間的變化關系，定性地分析了滲吸前緣長度與黏土含量和礦化度的關系。Dehghanpour等[23]通過一系列實驗詳細地研究了頁巖氣水和油氣滲吸特點，發現黏土的吸水膨脹以及產生的誘導裂縫造成了過量吸水現象。Yang等[24]基于實驗方法研究了三種致密儲層巖石滲吸特征，發現巖石吸水量與黏土含量成正相關。此外，通過巖石吸水量與時間的關系，評價了巖石孔隙結構的連通性。Sun等[25]研究了動態滲吸條件下，致密油砂巖孔隙內油水運移特征。

需要說明的是，近些年關于非常規儲層巖石滲吸特征的認識，主要集中在氣水滲吸系統，關于致密儲層巖石油水滲吸系統的研究仍不充分。本文基于NMR技術系統地研究了在不同影響因素下致密砂巖油水滲吸特征，并用氣水滲吸實驗與之進行對比分析。文章首先介紹了實驗所用的樣品信息及實驗方法，之后分析了滲吸過程中，油水和氣水運移特征及不同邊界條件下致密砂巖滲吸特點。此外，也探討了標度模型在致密砂巖流體滲吸條件下的適用性。

1 實驗樣品及方法

1.1 實驗材料

本文實驗所用巖樣來自鄂爾多斯盆地長8組儲層的露頭巖心，經鉆取并打磨加工成圓柱體，直徑在2.5 cm左右，長度約為5.0 cm。巖心孔隙度的測量利用氣測法，氣體介質選用氦氣。滲透率通過氣測并經克氏修正得到，基本物性數據見表1。

此外，XRD(X射線衍射技術)用來分析巖石礦物成分及其相對含量，見表2。可以發現致密砂巖礦物成分以石英、鉀長石、斜長石及黏土礦物為主，總量近70%。黏土礦物(占19.5%)主要以高嶺石為主，含有少量的伊利石和伊蒙混層。Dehghanpour等[23,26]發現頁巖具有強烈的吸水性，誘發了微裂縫的擴展，巖石吸水量與黏土含量成正比關系。Yang等[24]通過對比致密砂巖和泥頁巖的滲吸實驗，進一步發現造成黏土水化膨脹的主要礦物為蒙脫石和伊蒙混層。本文采用的致密砂巖黏土礦物以高嶺石為主，屬于非膨脹型礦物且不易水化[27]，因此黏土礦物對滲吸的影響較小。

對于油水系統，NMR實驗只需要獲取油中氫離子的信息，通過配置質量分數為40%的MnCl2溶液作為濕相流體(DIWm)[28]可以將水中的氫離子信號屏蔽掉，測得其密度為1.32 g/cm3，品氏黏度計測得其黏度為2.47 mPa·s；5#標準白油作為非濕相流體，密度為0.85 g/cm3，黏度為7.97 mPa·s；氣水系統中非濕相為空氣，濕相流體為去離子水(DIW)，其密度為1.0 g/cm3，黏度為1.1 mPa·s，以上流體黏度和密度的測量均在常壓室溫25℃下完成。此外，采用聚四氟乙烯分散液封堵巖石不同部位以構造不同滲吸邊界條件。

表1 實驗所用巖石樣品和流體信息Table 1 Properties of rock samples and fluids for experiments

表2 XRD礦物分析結果Table 2 Results of XRD mineralogy analysis

1.2 實驗原理及設備

低場NMR技術目前已廣泛應用于石油工程領域，如核磁測井[29-30]、巖石孔徑分布[31]、巖石表面潤濕性測定[32-33]及非常規儲層巖石油水、氣水滲吸特征研究[21-22,25,34]等。當核磁共振T2譜以表面弛豫為主時，橫向弛豫時間的大小可以反映巖石孔隙大小，弛豫時間越長，孔徑越大；信號強度反映了H質子的數量，與孔隙內流體體積成正比關系，詳細原理介紹見文獻[28, 30]。

本文所用實驗儀器為MesoMR23-60H-I中尺寸核磁共振分析儀，該設備兼具核磁共振波譜與成像技術，共振頻率23.408 MHz，磁體強度0.55 T，線圈直徑為32 mm，磁體溫度為32 ℃。回波時間TE設為0.2 ms，這樣較小的孔隙的信號將不會丟失[35]。

由于MnCl2溶液消除了滲吸過程中來自水的H質子信號，T2譜只體現了孔隙中白油的H質子信號，與橫軸所圍面積反映了巖心內部剩余油量，由此可以計算不同時間段滲吸采收率，計算公式為：

上式中，E為采收率，S0和S1分別表示初始條件下和不同滲吸時刻T2譜與橫軸所圍面積。

氣水系統實驗時，如果仍采用T2譜面積計算采收率，首先需要采取加壓飽和的方式將水注入巖樣后，測得其完全飽和水時的T2譜，之后再將巖樣烘干消除水分；對于致密巖樣而言，部分束縛水可能難以消除，對氣水滲吸造成一定影響，因此氣水實驗采用稱重法計量不同時間段進入巖樣的水體積，除以孔隙體積即可求得采收率。

1.3 實驗步驟

1.3.1 油水系統滲吸實驗

(1) 從露頭砂巖上鉆取直徑約為2.5 cm、長約5.0 cm的巖心，放置在恒溫箱里105℃下烘干48 h，測量孔隙度和滲透率。

(2) 在恒溫條件下分別對巖心飽和5#白油，飽和壓力25 MPa，飽和時間為48 h。根據巖樣飽和前后質量，計算得到含油飽和度，結果見表1。將聚四氟乙烯分散液均勻涂抹于巖樣特定滲吸面并靜置48 h，構造AFO、LEC、TEO和TEC邊界條件。

(3) 配制質量分數為40%的MnCl2滲吸液，分別倒入標有巖心編號的容器中備用。

(4) 將各巖樣放置在試管中，并放入核磁共振分析儀內，記錄初始核磁共振T2譜。

(5) 將巖心分別用繩懸掛浸入對應編號的容器中，保證巖心所有表面與溶液充分接觸，見圖1。分別記錄巖樣開始滲吸的時間，在特定的時間點取出巖樣，測得其核磁共振T2譜，之后再放入對應容器內，繼續滲吸實驗直至完成全部測試點。本文選取的時間間隔分別為0.5 h、3 h、6 h、19 h、43 h、73 h和120 h時，120 h時各巖樣T2譜基本不再變化，認為滲吸已結束。

1.3.2 氣水系統滲吸實驗

氣水實驗操作步驟與油水系統類似，需要在不同時刻將樣品取出記錄其T2譜并進行稱重，然后迅速放入容器中，設定的時間點分別為0 h、1 h、3 h、6 h、20 h、26 h、47 h、72 h、84 h、120 h，此時各樣品質量幾乎不再變化。

2 實驗結果分析與討論

本部分基于NMR T2譜實驗結果，分析了致密砂巖油水和氣水滲吸系統中流體運移特征。此外也討論了邊界條件對致密砂巖滲吸的影響和標度模型的應用。

2.1 致密砂巖自發滲吸流體運移特征

2.1.1 油水系統

由于MnCl2溶液消除了水中H質子信號，測得的T2譜反映了不同孔喉內剩余油的分布情況。在T2譜以表面弛豫為主時，弛豫時間與巖石孔徑大小存在對應關系[30,36]，根據Gao等[37]和Lai等[21]的研究結果，將巖樣T2譜弛豫時間分為四個區間，分別對應不同孔隙類型，見表3。

圖1 致密砂巖滲吸實驗Fig. 1 Imbibition experiment of tight sandstone

表3 T2弛豫時間與不同孔隙類型對應關系Table 3 Relationship between relaxation time and different types of pores

致密砂巖微觀孔隙結構具有多尺度分布特征，主要以納微米孔隙為主[38-39]，圖2為4個樣品滲吸過程中不同時間段T2數據半對數分布圖，弛豫時間從0.01 ms到1000 ms均有分布，符合致密砂巖孔隙結構特征，各個樣品分別對應不同的滲吸邊界條件，見表1。從圖中可以看出，初始時刻時巖心飽和白油，各樣品T2譜呈現出一致的分布特征，大體上存在三個波峰(見L4樣標注)。左峰處于0.01～8 ms之間，分布面積占總體分布的66.7%；中峰處于10～200 ms之間，約占總體分布的33%；200 ms之后存在第三個峰，占總體分布的0.3%左右，幾乎可以忽略不計，根據表3可知，L樣中的大孔隙數目極少，主要以微孔和小孔為主(67%左右)。

隨著滲吸的進行，巖石孔隙中的白油不斷被MnCl2溶液置換出去，T2譜信號幅度不斷下降，其與橫軸包圍面積不斷減小。從圖2可以看出，除L5外，其余巖樣的T2譜均由原來的多峰分布逐漸演變為單峰分布。由于L5樣滲吸邊界條件的設置，造成其采收率較低，T2曲線在滲吸過程中變化不明顯。以L2樣不同時間T2譜為例，左峰下降較快，而中峰下降幅度較小，左峰主要代表微小孔隙中含油量的變化，說明滲吸采油主要源自微小孔隙的貢獻。毛管力大小不僅與流體類型有關，也取決于巖石孔隙結構[40-41]。孔徑越小，毛管力越大，因此水會優先被吸入更小的孔隙中。Yang等[24]認為小的孔隙雖然有著較大的毛管力，但孔隙壁面的粗糙度、水膜厚度、納米孔隙邊界層厚度等可能不利于流體運移[42]，即滲吸速度與有效毛管力有關。由圖2中(a)可知，微孔和小孔中滲吸阻力與毛管力相比相對較小，依然是毛管力主導滲吸。中峰主要體現了中孔的滲吸貢獻，由于其孔徑較大，毛管力較小，因此滲吸速度相對較慢，波峰不斷下降直至消失，最終T2譜呈現單峰分布，從圖中可以看出非濕相主要滯留在小孔中。此外，直觀上看中孔(10～100 ms)滲吸采油速度甚至大于小孔(1～10 ms)，這是由于滲吸開始后，水會優先進入微孔隙，微孔與小孔和中孔隙之間通過細小的喉道進行連接，考慮到致密砂巖表面親水，較大孔隙表面部分被水相占據，而水相中H質子被Mn2+屏蔽，間接影響了油相流體的表面弛豫，造成中孔被識別小孔，因此滲吸結束后T2譜是以小孔為主的單峰分布。

圖2 不同樣品油水滲吸過程中T2譜的分布變化情況Fig. 2 Variation of T2 spectrum of different samples in different stages of imbibiton experiments, (a), (b), (c) and (d) represent that for sample L2, L4, L5 and L11 respectively

Liang等[30]以鹵水、煤油和絕緣油為實驗流體，發現當巖心內部油水共存時，致密和常規巖樣中流體依然以表面弛豫為主；而中低滲巖樣中，流體表面弛豫時間小于體積弛豫，此時油相流體的體積弛豫不可忽略，T2譜呈現出明顯的雙峰特征。圖3為本文滲吸實驗所用白油和去離子水的體積弛豫曲線，其中去離子水作為參照，可以看出白油的體積弛豫是在200 ～300 ms；而滲吸末期T2譜峰值對應的弛豫時間是在1～5 ms，顯然滲吸過程中核磁T2譜主要以白油的表面弛豫為主。

圖3 白油和去離子水的體積弛豫Fig. 3 The bulk relaxation T2 distributions of white oil and deionized water

2.1.2 氣水系統

氣水系統不同時間T2譜分布見圖4。氣水滲吸實驗類似于油水滲吸前巖心飽和油的過程，隨著水不斷滲吸進入巖樣，水中H質子信號幅度不斷增大，T2譜與橫軸所圍面積不斷增大。四個樣品除B8外，T2譜分布特點相同；以B6為例，初始時刻巖心內部只有微弱的H質子信號，主要來自黏土水或含結晶水的礦物，T2譜呈雙峰分布，右鋒位于100～1000 ms，與油水滲吸實驗類似，可以忽略不計；左峰位于微孔區間，微孔隙有著更強的毛管力，因此水更容易進入。滲吸早期左峰增長速度較快，之后速度逐漸變慢，主要因為水不斷進入巖心，毛管力作用逐漸被削弱，而右峰幾乎沒有變化。B8樣T2譜信號幅度增長速度明顯慢于其他巖樣，是由于巖樣的主要滲吸通道被封閉，滲吸速度變慢所致。

與油水系統最大不同在于直至滲吸結束，氣水系統T2譜始終是以微孔為主的單峰分布，而油水系統出現了以中孔為主的中峰，原因在于油水滲吸初始狀態是通過加壓飽和方式實現的，白油靠外力被驅替進入孔隙，飽和度均在95%左右，見表1。而水是依靠毛管力自發滲吸進入巖樣中，巖石內部實際上是由復雜的孔隙網絡組成，濕相與非濕相流體在不同孔喉中流動，如圖5所示。NWP和WP分別代表非濕相和濕相，由于孔徑不同小孔隙中毛管力較大，WP在較小孔隙中流速較快，見圖(a)所示，小孔隙中的濕相流體會比較大孔隙中流體提前到達相連的喉道處，進入較大孔隙后非濕相卡斷在孔隙中[17,43]，見圖5(b)。對于氣水滲吸而言，氣泡會滯留在較大孔隙中，因此這部分孔隙的信號在T2譜上未形成波峰。

圖4 不同樣品氣水滲吸過程中T2譜的分布變化情況Fig. 4 Variation of T2 spectrum of different samples under imbibition conditions

2.2 邊界條件對滲吸采收率的影響

由以上分析可知，邊界條件對滲吸采油的速度和最終采收率都有重要的影響。本節通過NMR技術研究邊界條件對致密砂巖滲吸的影響，各樣品采收率與滲吸時間關系見圖6。

圖5 較大孔隙中非濕相卡斷示意圖Fig. 5 Snap-off of NWP in larger pores

從圖6可以看出，對于同一滲吸系統不同邊界條件下滲吸初期采油速度較快，隨著滲吸的不斷進行，流體間的黏滯力作用愈加顯著，毛管力逐漸減小，滲吸速度逐漸減慢。AFO條件下，巖石有著最快的滲吸速度和采收率，其次是LEC和TEO，TEC最慢，采收率最低。圖2中，不同時間L5樣T2譜沒有明顯的變化，這是由于巖石整個側面都被封堵，而巖石側面是滲吸采油的主要作用通道，端面滲吸貢獻較小。然而，傳統觀點認為如果實驗采用的多孔介質和流體一致，那么邊界條件只影響滲吸的速度，而對最終采收率幾乎沒有影響[15,44-45]。近年來，一些學者[21,46-47]發現邊界條件不僅影響滲吸采油的速度，對采收率也有重要的影響，但是卻沒有給出合理的解釋。本文得到了與之相似的結論，從圖6可以明顯看出，L5樣滲吸速度和最終采收率均低于其他樣品。造成這一異常的原因，主要是由于常規中高滲巖樣均質性較好，其滲透率和相對滲透率沒有明顯的方向性，因此在巖樣與濕相流體接觸面積較小時(TEO邊界條件)，其滲吸速度顯著降低，但是采收率與其他邊界條件差別不大。而致密巖樣一般孔滲遠低于常規儲層巖樣，孔隙結構復雜且具有各向異性的特點，也就是說其滲透率、相滲和毛管力曲線等常常具有方向性。由于L5樣為TEO邊界條件，滲吸模式為一維線性流，而其他邊界條件主要以徑向流為主，L5采收率低的原因主要是軸向上巖石物性差于徑向方向，造成了TEO邊界下滲吸采收率遠低于其他條件。

圖6 致密砂巖滲吸采收率與時間關系Fig. 6 Recovery versus imbibition time under different boundary conditions

B8樣的邊界條件與L5均為側面封閉，從圖4可以看出其滲吸速度最慢。氣水系統跟油水系統區別在于，不同邊界條件下最終采收率大致相當，最大采收率(B6)與最小采收率(C5)相差僅9.3%，遠低于油水系統(L2與L5)的78.7%。造成這一問題的主要原因是：多孔介質滲吸過程由多種作用力共同主導，不考慮流體與礦物化學反應，毛管力為滲吸過程的主要動力；隨著濕相流體逐漸進入孔隙中，毛管力作用逐漸減弱，濕相與非濕相流體之間黏滯力作用顯著增大，削弱了毛管滲吸作用，因此滲吸速度減慢。與此同時，濕相流體飽和度只有達到一定程度時，才具有可觀的流動能力。而非濕相流體隨著滲吸的進行流動能力不斷下降，不利于滲吸采油。滲吸過程中各種因素相互競爭、相互制約，對氣水滲吸而言，流體之間的黏滯力、流度的變化及巖石各向異性造成的滲流阻力不足以對毛管力造成太大影響，因此不同邊界條件下最終采收率變化不大。

2.3 滲吸實驗數據的標度

滲吸標度模型用來分析滲吸數據，對于預測油藏滲吸產油量有著非常重要的參考價值。如果標度模型能夠將更多的因素考慮進去，那么將有更高的預測精度以及更廣的適用范圍。過去幾十年中，許多學者基于毛管擴散方程和達西定律來建立無因次標度模型[10-11,16,48-52]。其中最早的是由Mattax和Kyte[51]提出的標度模型，見式(2)，該模型有著嚴格的使用條件，適用性較差。Ma等[48]對該模型進行了修正，考慮了不同邊界條件下的特征長度，同時也考慮了非濕相流體的黏度，見式(3)。當兩種流體黏度比較大時，黏度項通過幾何平均處理將帶來較大誤差，Mason等[52]對此進行了修正并用大量的實驗數據驗證了該模型的合理性，見式(4)。值得一提的是，雖然Schmid和Geiger[11]提出了一個廣適性的解析解標度模型用來研究水濕和混合潤濕性油藏滲吸特征，但模型相對復雜需要迭代求解，此外模型還需要相滲曲線和毛管力曲線，這些參數對于致密儲層巖石來說都難以通過實驗獲取。本文采用式(4)對滲吸數據進行標度分析。

由于各式包含物理量基本一致，現統一對上面各式進行說明，σ為油水界面張力，N/m；K為巖石滲透率，m2；φ為孔隙度；L為樣品高度，m；μw和μnw分別為濕相和非濕相流體黏度，Pa·s；t為時間，s；Lc為特征長度，m；Vb為巖樣體積，m3；Ai為第i個滲吸面的面積，m2；li為第i個滲吸面到不流動邊界的距離，m。

Schechter[53]定義了無量綱數N，即毛管力與重力作用的比值來衡量重力對于滲吸的影響，該值大于5時，重力作用可以忽略。從表1計算結果來看，各樣品N均大于5，因此本文采用的標度模型忽略重力作用是合理的。

式中，C一般取為0.4；Δρ為濕相與非濕相密度差，kg/m3；g為重力加速度，m/s2，其余參數與上述定義一致。

為了更好地說明標度滲吸數據的必要性，作出未標度前各巖樣滲吸數據隨時間變化的曲線，見圖7。標度滲吸數據時，不同邊界條件下特征長度采用式(5)進行計算，然后將巖石與流體基本參數(見表1)分別代入式(4)，可以得到不同樣品采收率隨著無因次時間tD的關系曲線，最終結果見圖8。

圖7 滲吸采收率與時間關系Fig. 7 Recovery from imbibition versus time

圖8 滲吸采收率與無因次時間關系Fig. 8 Imbibition data plotted with dimensionless time

從圖7可以直觀看出油水和氣水系統不同邊界條件下，巖石滲吸速度及滲吸采收率的變化情況，但整體數據非常分散。圖8為基于式(4)進行標度之后的結果，可以看出數據的離散性得到了顯著的改善，各樣品滲吸數據點趨向于一條曲線。考慮到巖樣黏土水化較弱且各樣品潤濕性一致，采用Mason等[52]提出的模型來標度滲吸數據是合理的，側面反映出基于達西定律來預測致密儲層滲吸采收率在一定程度上是可以接受的。L5的偏離是由于在TEO邊界條件下，特征長度是基于巖心兩端均勻采油的假設來計算的，而實際實驗中由于巖樣兩端面孔隙大小的不同會造成毛管回壓不同，兩端面采油并不是對稱的[54]，因此滲吸數據的標度出現了較大的偏差。此外，致密巖樣各向異性

導致TEO邊界條件下的巖石滲透率與其他邊界條件下的滲透率差異也會導致標度的偏差。

3 結論

本文基于NMR技術研究了致密砂巖自發滲吸過程中，孔隙內流體運移特征及不同邊界條件對于滲吸的影響，主要取得以下認識：

(1) NMR T2譜能夠反映飽和流體巖樣內部孔隙大小及分布，致密砂巖微觀孔隙結構具有多尺度分布特征，主要以微孔和小孔為主，對應的T2弛豫時間小于10 ms。

(2) 滲吸開始后，微孔中的白油會優先被水置換出來，通過較小的孔隙流出。滲吸結束時，巖樣T2譜演變為以小孔為主的單一波峰，說明非濕相主要滯留在小孔中。這是因為巖石表面的親水性，造成中孔表面部分被水占據，影響了油相的表面弛豫，被NMR檢測為小孔所致，致密砂巖滲吸采油主要來自于微孔的貢獻。而氣水滲吸系統，由于卡斷效應，水相難以進入較大的孔隙中。

(3) 邊界條件對致密砂巖油水系統滲吸速度和最終采收率影響顯著，巖石表面與水接觸面積越大，滲吸速度越快，采收率也最高。對氣水系統而言，邊界條件的影響體現在滲吸速度上，對最終的滲吸采收率影響不大。

(4) 由于本文致密砂巖黏土水化作用微弱且各巖樣潤濕性相同，傳統的無因次標度模型能夠較好的標度油水和氣水滲吸數據，標度模型預測此類儲層滲吸采收率是合理的。