流體在巖石孔隙中的運(yùn)移規(guī)律及其流固耦合效應(yīng)是地下油氣儲(chǔ)備與開發(fā)的核心科學(xué)問題，也是導(dǎo)致不同工程災(zāi)害或工程難題的重要因素。精確表征巖石微觀孔隙結(jié)構(gòu)，揭示微觀孔隙結(jié)構(gòu)與流體輸運(yùn)特性的內(nèi)在關(guān)聯(lián)，是開展深部巖體相關(guān)工程研究的基礎(chǔ)。
近期，中國(guó)科學(xué)院武漢巖土力學(xué)研究所的宋睿副研究員、劉建軍研究員、楊春和研究員聯(lián)合西南科技大學(xué)的汪堯博士等人提出了一種利用3d打印和微ct成像技術(shù)實(shí)現(xiàn)致密砂巖復(fù)雜孔隙結(jié)構(gòu)定量表征和多相流體輸運(yùn)特性的可視化研究方法。研究團(tuán)隊(duì)利用新型的面投影微立體光刻技術(shù)（pμsl，nanoarch s130，摩方精密）實(shí)現(xiàn)了致密砂巖孔隙模型的原位尺度打印（~2μm光學(xué)分辨率），再現(xiàn)了致密砂巖復(fù)雜孔隙系統(tǒng)的三維拓?fù)浣Y(jié)構(gòu)特征與空間連通性。研究人員對(duì)比分析了3dp巖心與數(shù)字巖心（drp）模擬得到的孔徑分布（psd）、孔隙度和絕對(duì)滲透率的差異；同時(shí)結(jié)合原位ct成像技術(shù)開展了3dp巖心可視化co2驅(qū)油實(shí)驗(yàn)，并與實(shí)驗(yàn)基準(zhǔn)數(shù)據(jù)進(jìn)行了比較。研究成果為定量表征巖石復(fù)雜孔隙結(jié)構(gòu)特征及其中多相流體輸運(yùn)機(jī)制提供了新的工具，具有廣闊的應(yīng)用前景。論文研究工作得到國(guó)家自然科學(xué)基金，武漢市知識(shí)創(chuàng)新專項(xiàng)（基礎(chǔ)研究）和四川省自然科學(xué)基金等項(xiàng)目的支持。相關(guān)研究成果以“3d printing of natural sandstone at pore scale and comparative analysis on micro-structure and single/two-phase flow properties”為題發(fā)表在《energy》期刊上。
圖1.基于ct圖像與面投影微立體光刻技術(shù)的致密砂巖微觀孔隙結(jié)構(gòu)提取與3d打印制備流程(a)天然致密砂巖的微ct掃描；(b)數(shù)字圖像處理與巖心重建；(c)面投影微立體光刻3d打印成型
該研究中所采用的天然巖心樣本為海相致密砂巖。通過從原始巖心中鉆取直徑約為5mm的小巖心柱塞樣本，利用蔡司xradia microxct-400三維成像系統(tǒng)進(jìn)行微ct掃描成像，獲取天然巖心孔隙結(jié)構(gòu)的微ct圖像（如圖1a所示），并將其用于孔隙空間提取、數(shù)字巖心重建與模擬（如圖1b）；然后，基于數(shù)字圖像處理轉(zhuǎn)化為3d打印通用的.stl文件，利用bmf公司的面投影微立體光刻成型技術(shù)完成孔隙模型的3d打印（如圖1c所示）。
圖2. 3d打印巖心與天然巖心微觀孔隙結(jié)構(gòu)的對(duì)比分析(a)基于偏光顯微鏡和ct成像得3dp巖心孔隙結(jié)構(gòu)表征；(b)基于圖像校準(zhǔn)的3dp巖心與原始巖心孔隙結(jié)構(gòu)拓?fù)湫螒B(tài)特征的對(duì)比分析；(c)孔隙結(jié)構(gòu)特征參數(shù)的計(jì)算與分析
為表征3d打印巖心在復(fù)刻天然巖心孔隙結(jié)構(gòu)特征方面的準(zhǔn)確性，該團(tuán)隊(duì)分別采用偏光顯微鏡和微ct成像對(duì)3dp巖心的2d/3d微觀孔隙結(jié)構(gòu)特征進(jìn)行了定量表征（如圖2a所示）?；趫F(tuán)隊(duì)自行開發(fā)的數(shù)字圖像處理與模型重建技術(shù)，分別研究了3dp巖心孔隙分布特征，并與天然樣品的實(shí)驗(yàn)室測(cè)試結(jié)果進(jìn)行了對(duì)比分析，結(jié)果表明3dp巖心和原始樣品的psd分布總體上一致（如圖2c所示）。在對(duì)3dp巖心和原始巖心ct圖像手動(dòng)校準(zhǔn)的基礎(chǔ)上，團(tuán)隊(duì)采用開源圖像處理軟件（fijiyama）中的塊匹配算法（block-matching algorithm）實(shí)現(xiàn)了3dp巖心ct圖像與原始樣品ct圖像的自動(dòng)配準(zhǔn)，并作為后續(xù)分析的基準(zhǔn)數(shù)據(jù)（如圖2b所示）。結(jié)果表明，3dp巖心與原始巖心孔隙特征吻合較好，驗(yàn)證了3dp巖心在微米尺度下再現(xiàn)巖石微觀結(jié)構(gòu)的可行性和適用性。
在此基礎(chǔ)上，團(tuán)隊(duì)以分割的微ct圖像為數(shù)據(jù)藍(lán)本，引入峰值信噪比（peak signal-to-noise ratio, psnr）和結(jié)構(gòu)相似性指數(shù)度量（structural similarity index measure, ssim）兩個(gè)關(guān)鍵參數(shù)對(duì)3dp巖心孔隙結(jié)構(gòu)特征進(jìn)行表征，以量化3dp巖心與原始巖心孔隙結(jié)構(gòu)的保真度（如圖2c所示）。psnr用于衡量相同空間位置上孔隙特征參數(shù)（大小和坐標(biāo)位置）的絕對(duì)誤差。ssim用于測(cè)量?jī)蓚€(gè)圖像之間的相似性，用于評(píng)估相應(yīng)位置上的孔隙是否由3d打印機(jī)識(shí)別。計(jì)算結(jié)果表明：本文中3dp巖心的psnr值介于[9.010，14.983]之間，其ssim值介于[0.870，0.925]之間。大多數(shù)孔隙特征被打印識(shí)別，但一些孔隙并不在原始尺寸或位置上。由于后處理過程中，樣品近端部的液體樹脂更容易被去除，因此頂/底部結(jié)構(gòu)的打印精度優(yōu)于其他部分，顯示出更高的ssim值。
圖3. 基于原位ct成像的微觀可視化多相滲流試驗(yàn)(a)團(tuán)隊(duì)自行設(shè)計(jì)的用于原位ct成像的微觀可視化滲流試驗(yàn)系統(tǒng)；(b)3dp巖心飽和油狀態(tài)（上部）和co2驅(qū)油后（下部）3dp巖心中油相分布的微ct圖像；(c)co2驅(qū)油后3dp巖心中co2分布及對(duì)應(yīng)的孔隙網(wǎng)絡(luò)模型，以及3dp巖心和原始巖心中殘余油相原位潤(rùn)濕角計(jì)算結(jié)果的對(duì)比
在3dp巖心與原始巖心孔隙結(jié)構(gòu)特征對(duì)比分析的基礎(chǔ)上，團(tuán)隊(duì)針對(duì)3dp巖心的流體輸運(yùn)特性開展了進(jìn)一步的研究。利用自行設(shè)計(jì)的基于原位微ct成像的可視化滲流試驗(yàn)系統(tǒng)分別進(jìn)行了3dp巖心的飽和油和co2驅(qū)油試驗(yàn)（如圖3a所示）。分別采集了飽和油狀態(tài)與驅(qū)替完成時(shí)3dp巖心的微ct圖像（如圖3b所示）。為了消除不同掃描階段樣品放置的人為誤差，研究人員對(duì)獲取的ct圖像也進(jìn)行了手動(dòng)校準(zhǔn)和圖像配準(zhǔn)操作。分析結(jié)果表明：注入co2氣體主要沿孔隙中部流動(dòng)，導(dǎo)致顆粒表面出現(xiàn)大規(guī)模殘余油。考慮到制備3dp巖心使用的htl樹脂是強(qiáng)油濕性，殘余油相優(yōu)先附著到固體表面。當(dāng)注入流體發(fā)生突破時(shí)，樣品中會(huì)留下很大部分以油膜形式分布的殘余油。在油濕性巖心中，毛細(xì)管壓力是注入co2的阻力，導(dǎo)致大量殘留油塊被毛管力卡斷在小孔中。此外，研究團(tuán)隊(duì)對(duì)3dp巖心和原始巖心的原位接觸角進(jìn)行了計(jì)算與對(duì)比分析，討論了微觀潤(rùn)濕性在殘余流體捕獲機(jī)制中的影響（如圖3c所示），并進(jìn)一步提取了co2驅(qū)替后3dp巖心的孔隙網(wǎng)絡(luò)模型，對(duì)驅(qū)替過程中co2氣體的主要滲流通道以及微觀賦存狀態(tài)進(jìn)行了討論與分析。結(jié)果表明，注入氣體主要沿3dp巖心的左側(cè)分布，注入co2沿優(yōu)先通道突破，與剩余油分布一致?？紤]到注入co2的操作壓力低于最小混相壓力，驅(qū)替過程為不混相氣-液流，界面張力和注入流體粘度的降低有助于提高波及效率和采收率。（如圖3c所示）。