巖心微流動的核磁共振可視化研究

狄勤豐，華 帥，顧春元，＊，葉 峰，龐東山，蔣 帆，楊培強

巖心微流動的核磁共振可視化研究

狄勤豐1，2，華 帥1，2，顧春元1，2，＊，葉 峰1，2，龐東山1，2，蔣 帆1，2，楊培強3

（1.上海大學(xué)上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所，上海 200072；2.上海大學(xué)上海市力學(xué)在能源工程中的應(yīng)用重點實驗室，上海 200072；3.上海紐邁電子科技有限公司，上海 200333）

巖心微流動可視化是研究化學(xué)驅(qū)油微觀機理的一項重要流體實驗新方法，重點介紹采用低場核磁共振成像技術(shù)研究天然巖心中流體分布可視化的新進展。提出和分析了國產(chǎn)低場核磁共振成像巖心驅(qū)替裝置面臨的圖像不清、材料干擾等問題，通過合理選材、優(yōu)化參數(shù)，從硬件和軟件2方面進行了改進與優(yōu)化，消除了干擾因素。開展了天然巖心的驅(qū)替實驗，采集了油水的實時NMR信號與MRI成像信號，以及不同驅(qū)替階段油水的NMR－T2譜，得到分辨率較高的油水分布圖像。結(jié)果顯示殘余油隨著驅(qū)替PV數(shù)的增加而減少，具有初期減少明顯而后趨緩的特點，并發(fā)現(xiàn)巖心中存在端部油殘滯現(xiàn)象，其范圍距端部4mm左右。研究了通過獲得的油水分布圖像計算含油飽和度的方法，其結(jié)果與傳統(tǒng)方法一致，誤差在10%以內(nèi)，這也說明了油水分布圖像的可靠性。這不僅為計算飽和度提供了一種新方法，而且該方法的另一個優(yōu)勢是可以分析任意局部位置的油水飽和度。研究結(jié)果表明，在研究巖心微流動過程中流體分布的可視化方面，核磁共振成像技術(shù)是值得深入研究的新方法。

微流動可視化；低場核磁共振；油水分布；端部殘滯；飽和度

0 引 言

探明化學(xué)驅(qū)油微觀作用機理是研究化學(xué)驅(qū)替提高原油采收率技術(shù)的關(guān)鍵，巖心可視化方法是研究微管道中驅(qū)油機理的重要手段之一，也是當(dāng)今石油工程、波譜學(xué)與流體力學(xué)領(lǐng)域的研究熱點方向之一［1－2］?？梢暬夹g(shù)與巖心流動實驗相結(jié)合已經(jīng)成為研究各種注劑驅(qū)替機理及實時效果可信的重要實驗流體力學(xué)新方法。傳統(tǒng)流動實驗中，不可視的巖心就如同一個黑箱，流體在“黑箱”內(nèi)流動的狀態(tài)以及分布狀況無法得知，只能通過檢測進、出口壓力、流量和巖心的電阻率等宏觀參數(shù)來實現(xiàn)巖心內(nèi)流體狀態(tài)的間接表征［3］，因此，驅(qū)替劑的作用機理比較含糊。另外，人們認為在管流及多孔介質(zhì)的流動過程中，存在端部效應(yīng)，而由于傳統(tǒng)檢測方法的局限性，其形態(tài)特征與形成機理尚不清楚，只能采取改造端部管路以及改變注入速度、反向驅(qū)替等手段來減輕阻塞程度。

目前，巖心內(nèi)部流體分布的可視化方法主要有平面微觀驅(qū)替模型、計算機層析掃描技術(shù)（CT）以及核磁共振技術(shù)（NMR）等［4－6］。

微觀驅(qū)替模型主要包括玻璃壓片模型和玻璃蝕刻模型。玻璃蝕刻模型保留了以巖心薄片的孔隙系統(tǒng)為基礎(chǔ)的孔隙網(wǎng)絡(luò)，具有與巖心相似的平面孔隙結(jié)構(gòu)形態(tài)，常用來研究各種驅(qū)替劑的作用機理［7］。但是其缺點也相當(dāng)明顯：（1）制作模型的過程復(fù)雜，需要將巖心二維特征刻蝕到特殊材質(zhì)的表面；（2）模型所用的壁面材料與真實情況相差較大，不能準確地反應(yīng)驅(qū)替劑與壁面的作用特征；（3）模型為二維流動模型，難以很好地表現(xiàn)空間結(jié)構(gòu)中的流動狀態(tài)［8－9］。

CT技術(shù)是一種分辨率較高的成像技術(shù)，其無損的特點使其能夠得到巖心內(nèi)部真實狀態(tài)的孔隙結(jié)構(gòu)，目前已經(jīng)成為研究巖心可視化的一項重要方法，在構(gòu)建數(shù)字巖心、獲取巖心骨架信息方面發(fā)揮了重要作用［10－11］。對于CT成像技術(shù)來說，分辨率和觀察視野是一對矛盾，要達到納米級的分辨率，其測試區(qū)域就大約在微米級別的范圍，對巖心來說這個測試范圍偏?。涣硗?，CT的效果還取決于操作者對錄取資料缺陷的認識與分析能力，以及對基于灰度閾值的體積分割是否合適［12］。

核磁共振成像技術(shù)（MRI）也是一種無損可視化方法，與CT一樣，在醫(yī)學(xué)領(lǐng)域具有重要的作用。近年來，核磁共振及成像技術(shù)在國內(nèi)的石油領(lǐng)域得到了長足的發(fā)展，為測錄井和儲層評價等提供了大量的實時檢測資料［13］。Feng Y X、Jiang L L等利用核磁共振技術(shù)進行了原油、巖心基本參數(shù)以及流體的分布研究，獲得了填砂管中CO2／水的分布圖像［14－15］。2011年，劉長福利用核磁共振的T2譜測試了含油飽和度和孔隙結(jié)構(gòu)的關(guān)系以及孔隙結(jié)構(gòu)的變化［16］。曲巖濤等利用高場核磁及離心驅(qū)替設(shè)備用核磁共振擴散－弛豫二維譜研究了水驅(qū)油的分布，獲得了油水在不同孔隙中的分布狀況以及潤濕性等信息［17］。劉瑜、宋永臣等以玻璃砂填砂管為介質(zhì)，用高場核磁共振設(shè)備得到了切片圖像，并以這些切片圖像為基礎(chǔ)，通過軟件重組了低分辨率的三維圖像［18］。2012年，郎東江、呂成遠等對巖心斷面采集信號，得到了巖心中裂紋的二維圖像［19］。

雖然多孔介質(zhì)中的核磁共振成像已經(jīng)有了一定的突破，但是實時驅(qū)替實驗中采集到的圖像分辨率較低，還不能很好地表達驅(qū)替過程中的流體分布情況，分析認為主要原因有以下幾個方面：

（1）目前的核磁共振成像驅(qū)替系統(tǒng)和配套分析軟件還不成熟。硬件上存在信號干擾因素，軟件上數(shù)據(jù)采集與分析功能存在缺陷，難以得到較好的實時檢測流體分布情況。

（2）天然巖心的微通道構(gòu)造與玻璃砂填砂管有很大的不同，玻璃填砂管中的流體分布一般比較均勻，而天然巖心非均質(zhì)性強，微細孔較多，流體分布不均。

（3）高場核磁在成像方面有著突出的優(yōu)勢，但是天然巖心的組分較為復(fù)雜，常含有少量的鐵磁性物質(zhì)，容易干擾高場核磁共振設(shè)備的磁場，從而難以得到真實巖心內(nèi)的流體分布圖像。而低場核磁共振儀在一定的范圍內(nèi)可以減弱這種影響。

3種巖心可視化方法各有特點，相比而言，核磁共振成像更適合測試巖心中的流體分布。

本文重點介紹國產(chǎn)低場核磁共振成像驅(qū)替裝置在使用過程中出現(xiàn)的問題及解決方法，利用核磁共振成像技術(shù)研究天然巖心內(nèi)流體的分布特征，以及利用核磁共振圖像計算飽和度的新方法。

1 低場核磁共振成像驅(qū)替系統(tǒng)及參數(shù)優(yōu)化

1.1核磁共振成像驅(qū)替系統(tǒng)存在的問題和解決措施

測試儀器為紐邁電子科技有限公司與上海大學(xué)聯(lián)合開發(fā)的低場核磁共振成像巖心驅(qū)替系統(tǒng)，為具有自主知識產(chǎn)權(quán)的國產(chǎn)設(shè)備。從硬件和軟件的完善程度看，該系統(tǒng)也是可用于進行實時檢測巖心流體分布的國產(chǎn)裝置，包括巖心驅(qū)替系統(tǒng)和小型核磁共振成像（MINI－MRI）系統(tǒng)2個部分。驅(qū)替系統(tǒng)由HA－IV型核磁驅(qū)替裝置、巖心夾持器以及驅(qū)替泵組成。整個系統(tǒng)除了采集核磁信號外，還可以采集進、出口壓力和流量，用來輔助評價液體的流動狀態(tài)。

利用裝置開展了系列實驗，發(fā)現(xiàn)硬件和配套軟件都存在較多的問題，并進行了改進。主要硬件問題有：

（1）首先是圍液的選擇，以確保圍液不會影響核磁共振采集信號。初期使用水作為圍液，干擾非常大，難以取得滿意的圖像，即使混入少量水，也會造成顯著的干擾，如圖1（a）所示；改用純氟油做圍液后，巖心及其中流體的信號能夠被很好地顯示出來，如圖1（b）所示。

圖1 混入純水和純氟油做環(huán)壓時信號對比Fig.1 Signals of different confining liquids

（2）巖心夾持器是重要的部件，既不能有干擾信號，又必須能承受一定的壓力。因此，采用何種生產(chǎn)材料非常重要。通過多次測試，選擇了一種非金屬聚合無磁材料。同時，與夾持器配套的快速接頭等連接部分也需要采用一種特殊轉(zhuǎn)換接頭進行信號屏蔽。

（3）第一代夾持器的結(jié)構(gòu)簡單，驅(qū)替水容易混入圍液中，而且不方便拆卸；常規(guī)巖心夾持器的結(jié)構(gòu)比較復(fù)雜，聚合無磁材料難以兼顧壓力與結(jié)構(gòu)，通過優(yōu)化設(shè)計、改進后的二代夾持器，解決了圍液與驅(qū)替液易混的問題，拆卸也十分方便。

（4）外部鐵磁性及信號等也會干擾內(nèi)部信號的采集，需要在夾持器兩端采用金屬網(wǎng)隔絕信號。

1.2核磁共振成像方法及參數(shù)優(yōu)化

采集核磁共振信號的主要步驟為：（1）在巖心驅(qū)

當(dāng)進行T1加權(quán)成像時，需要使用較短的Repetition Time（TR）值和較短的Echo Time（TE）值；當(dāng)進行T2加權(quán)成像時，需要使用較長的TR值和較長的TE值；當(dāng)進行一幅質(zhì)子密度成像時，需要使用較長的TR值和較短的TE值。

流體在多孔介質(zhì)中處于短弛豫的狀態(tài)，在這種情況下，較短的TE時間能夠使自旋回波的產(chǎn)生和檢測更早更快，而且數(shù)據(jù)點密度增加，有效信噪比也增加，提高了多孔介質(zhì)中流體分布的成像效果。

同時，在驅(qū)替過程中，利用MINI－MRI得到各個驅(qū)替階段的巖心成像。真實巖心圖像信噪比較差，巖心信號基本等于巖心夾持器的信號，亮暗無法很好的顯示其規(guī)律，偽圖也無法達到預(yù)期效果。針對這些問題，除了硬件上進行改進外，還通過軟件處理，增強其明暗差異性，獲得滿意偽彩圖。替過程中利用核磁共振軟件采集同一時間的T2譜與二維圖像，對驅(qū)替過程進行監(jiān)測；（2）對測試參數(shù)進了優(yōu)化，在成像脈沖序列上加90°射頻（Radio Frequency，RF）脈沖，經(jīng)過實時間τ再加180°脈沖，再過時間τ就出現(xiàn)一個自旋回波（Spin Echo，SE），采集信號數(shù)據(jù)；（3）對數(shù)據(jù)進行后處理，形成圖像。

本實驗核磁共振成像所用的脈沖序列為多層自旋回波（Multiple Spin Echo，MSE）。MSE是常用的一種核磁成像序列，不但可以實現(xiàn)SE的基本功能，而且可以一次進行多個層面的成像。

在一幅自旋回波圖像中，圖像上每點的信號強度可由式（1）來描述：

2 巖心驅(qū)替過程中流體的核磁共振成像實驗

2.1巖心驅(qū)替核磁共振成像實驗

利用低場核磁共振巖心流動成像裝置開展了水驅(qū)油流體分布成像實驗。實驗所用的多孔介質(zhì)為取自油田的天然巖心，其長度為5.752cm左右，直徑為2.5cm，孔隙度為17.70%，氯化錳溶液的濃度為2.5g／L，實驗溫度為24℃。

實驗流程如圖2所示，巖心夾持器水平放置，驅(qū)替方向從左到右。

實驗過程如下：

（1）將標樣水平放入巖心夾持器內(nèi)，調(diào)試成像系統(tǒng)參數(shù)（包括中心頻率、硬脈沖寬度、電子勻場和軟脈沖脈幅），以及標樣的MSE序列參數(shù)，得到合適的成像序列參數(shù)。

圖2 核磁共振驅(qū)替流程示意圖Fig.2 Schematic of set－up for examining core flow

（2）將標樣換成天然巖心，水平放入夾持器內(nèi)。將氯化錳溶液與柴油分別裝入2個儲液罐。開啟測試軟件，調(diào)節(jié)好流程，將柴油從左到右驅(qū)入巖心，關(guān)閉閥門使柴油在巖心內(nèi)老化24h。在保持環(huán)壓的情況下，完成飽和柴油巖心的T2譜以及MRI成像。

（3）用氯化錳溶液以0.5mL／min流速驅(qū)替柴油。驅(qū)替到0.1PV數(shù)時，關(guān)閉進出口閘閥，保持壓力，采集巖切片的圖像和T2弛豫，記錄出油、出水量。

（4）打開閘閥，繼續(xù)驅(qū)替，驅(qū)替每隔0.1PV，重復(fù)步驟（3）。測試直至T2譜與圖像基本不再變化，完成驅(qū)替，保存數(shù)據(jù)。

（5）在整個實驗過程中，采用環(huán)壓跟蹤技術(shù)，環(huán)壓始終大于驅(qū)替壓力2～5MPa。最后，利用專業(yè)軟件對信息數(shù)據(jù)進行處理，得到巖心不同驅(qū)替時刻的T2譜線和不同位置的圖像，利用圖像計算不同位置的殘余油飽和度。

2.2核磁共振成像檢測結(jié)果與分析

2.2.1水驅(qū)油過程中的T2譜

對這組數(shù)據(jù)進行反演，將回波串轉(zhuǎn)換成T2分布。巖心在不同驅(qū)替孔隙體積（PV）時的T2分布如圖3所示。

一般來說，氯化錳溶液的弛豫時間大致為0.1～10ms，而柴油的弛豫時間約為100～1000ms［19］。圖3給出了驅(qū)替量在0、1.5、2和12PV數(shù)時的T2譜，曲線的左峰代表氯化錳溶液，右峰代表柴油。根據(jù)圖中不同驅(qū)替PV數(shù)的T2譜，可以發(fā)現(xiàn)，隨著驅(qū)替PV數(shù)的增加，T2譜的左峰升高，右峰降低，說明巖心內(nèi)水的含量不斷增加，相應(yīng)的柴油含量不斷減小，表明巖

在實驗中，巖心內(nèi)同時存在著2種流體，這2種流體的弛豫信號不同，而不同弛豫時間的流體成分能夠通過反演的方法區(qū)分出來。CPMG脈沖序列測量磁化強度隨時間衰減由下式表示：心中的柴油不斷地被驅(qū)替出來。由圖4可見，在驅(qū)替2PV之前驅(qū)油效果明顯，2～3PV時，效果減緩，3PV之后，油水含量基本變化不大。

圖3 驅(qū)替過程中的T2弛豫譜Fig.3 T2spectrum during displacement

圖4 驅(qū)替過程中的油水T2譜峰面積變化Fig.4 The change of the peak area of the oil／water T2spectrum during displacement

圖5和6分別是水驅(qū)油在驅(qū)替不同PV數(shù)時巖心內(nèi)流體的MRI平視成像圖和俯視投影成像圖，左端是驅(qū)替液進口，右端為出口。

圖5 驅(qū)替過程中巖心內(nèi)流體的MRI平視投影成像Fig.5 MRI plan view images during forward displacement

圖6 驅(qū)替過程中巖心內(nèi)流體的MRI俯視投影成像Fig.6 MRI vertical views during forward displacement

由圖5可見，隨著驅(qū)替量的增加，巖心中流體的顏色由紅色（油）變成藍色（水），說明油逐漸被水所替代，含油量在減少。在4PV時，前2／3段基本為水，殘余油很少，這也說明該巖心中流體未出現(xiàn)明顯的重力分異作用。結(jié)合圖6看，驅(qū)替2PV時，2角部分含油明顯減少；驅(qū)替到4PV時，一側(cè)面附近已形成明顯的窄條水流通道，說明在驅(qū)替過程中巖心內(nèi)右側(cè)形成了優(yōu)勢通道，殘余油較多的分布在左側(cè)。整體上，在驅(qū)替4PV后，前段部分的含油量已經(jīng)較少，這與T2譜的結(jié)果基本一致。

同時，2圖都顯示，出口端的油信號直到8PV時依然比較強，并形成一個明顯的橫條帶狀分布，且主要處于中間，兩側(cè)有水流通道。研究發(fā)現(xiàn)端口原油殘留不是個別現(xiàn)象，在很多巖心成像中都存在。分析認為，這是由于端部效應(yīng)所引起的剩余油殘滯現(xiàn)象，當(dāng)水驅(qū)原油到接近末端時，在兩側(cè)已經(jīng)形成水流優(yōu)勢通道，由于末端阻力突然減小，形成水繞過阻力較大的原油的端部流動特征。該成像圖也是端部現(xiàn)象存在的直觀證據(jù)，并顯示了對驅(qū)替油產(chǎn)生的殘滯特點。通過測算，這種端部現(xiàn)象的范圍在4mm以內(nèi)。

2.2.3含油飽和度的核磁共振成像測試法

直觀的油水分布圖是研究驅(qū)油機理的重要依據(jù)，但是油水分布是否可靠，尚難以證實。含油飽和度是傳統(tǒng)流動實驗容易獲得的可靠參數(shù)，如果通過核磁成像測得的油水分布圖計算得到含油飽和度，那么不僅能獲得一種新的含油飽和度方法，而且可以驗證油水分布圖是否可靠。

核磁圖像的灰度是由信號值直接換算而來，每一個像素代表著這個像素范圍內(nèi)的流體信號幅度，也就包含了這個范圍內(nèi)的含油飽和度信息。飽和油時巖心內(nèi)的核磁信號累加就是所有油的核磁共振信號，而驅(qū)替后所有像素點的核磁信號的累加就是剩余油的信號。由此計算某時刻巖心內(nèi)不同位置的含油飽和度。

將飽和油的巖心設(shè)置為飽和狀態(tài)，并認為充分飽和油的巖心樣品，油飽和度為100%，而無油狀態(tài)的巖心，油飽和度設(shè)置為0。通過一定的序列采樣、灰度統(tǒng)一和數(shù)據(jù)處理，計算了天然巖心的油水分布圖。

巖心沿縱截面分3個切片成像，切片位置如圖7所示。圖8是其中巖心的水驅(qū)油過程中含油飽和度的變化曲線，包含了由3塊切片成像圖計算得到的3條飽和度曲線和1條實驗測試的飽和度曲線。

由圖8可見，用核磁共振圖像法計算得到的飽和度與傳統(tǒng)巖心含油飽和度的測試結(jié)果相比，變化趨勢一致，中間段吻合較好，前后端有一定的誤差，整體上誤差值在10%以內(nèi)。就各個切片的飽和度信息而言，切片1和切片3是位于巖心兩側(cè)對稱的切面，其飽和度的數(shù)值較為相似，而位于巖心中間的切片2在驅(qū)替進行到0.3PV時，飽和度下降得較快，這是由于在水驅(qū)油的過程中，中間形成了優(yōu)勢通道，位于通道中的柴油被快速驅(qū)替，而其余的柴油則殘留在巖心中。

圖8 水驅(qū)油過程中的含油飽和度Fig.8 Oil saturation during displacement

2種方法得到的飽和度的誤差主要來自3個方面：一是采用切片法，每塊切片的厚度為7mm，總厚21mm，而巖心直徑為25mm，有部分間隔未計算在內(nèi)；二是實驗技術(shù)的系統(tǒng)誤差，比如采集時間短，存在信號的漏捕，而且天然巖心的孔隙分布范圍廣，極小孔道的圖像信號很容易漏捕；三是實驗數(shù)據(jù)采用的計量法，存在計量誤差（如滯后、不穩(wěn)定等），尤其是前后出液不穩(wěn)定時，容易出現(xiàn)誤差。

整體上看，核磁共振可視化方法獲得的飽和度與測試計量獲得的飽和度相近，可以作為一種無損的巖心含油飽和度快速檢測方法，而且可以計算指定位置的局部含油飽和度。

綜上所述，通過對驅(qū)替實驗中同一時間采集的T2譜以及MRI成像圖的分析，可以從整體和局部的范圍共同觀測巖心中流體總量以及局部分布的變化。由不同方向獲取的圖像，可以對殘余油進行空間立體分析，并可以獲得空間立體的殘余油飽和度。

3 結(jié) 論

（1）通過對裝置軟、硬件的改造，解決了信號干擾等問題，使國產(chǎn)低場核磁共振成像巖心驅(qū)替系統(tǒng)達到了實時檢測驅(qū)替巖心中油水分布的要求。

（2）水驅(qū)替過程中含油量初期減少明顯，然后減幅變緩，并形成了水流優(yōu)勢通道。同時，該技術(shù)證實了端部現(xiàn)象的存在，顯示了端部效應(yīng)所引起的剩余油的殘滯，其影響范圍4mm左右，這對確定用于驅(qū)替實驗的巖心必須具備的長度具有一定的指導(dǎo)作用。

（3）核磁共振成像飽和度圖像法與傳統(tǒng)飽和度計量法的測試結(jié)果趨勢一致，相差在10%以內(nèi)，兩者基本吻合，這表明核磁共振可視化技術(shù)可以成為一種無損、快速的含油飽和度檢測法，而且該方法可以獲得局部指定位置的含油飽和度，這一點傳統(tǒng)的飽和度方法無法做到。

［1］張明安.二元復(fù)合體系微觀驅(qū)油機理可視化實驗［J］.油氣地質(zhì)與采收率，2013，20（3）：79－82.Zhang M A.Experiments on visualization of the microscopic mechanism of surfactant／polymer binary compound system［J］.PGRE，2013，20（3）：79－82.

［2］Song Y Q.Novel NMR techniques for porous media research［J］.Cement and Concrete Research，2007，37（3）：325－328.

［3］顧春元，狄勤豐，施利毅，等.納米顆粒在儲層微流道中的減阻機理實驗研究［J］.實驗流體力學(xué)，2010，24（6）：6－13.Gu C Y，Di Q F，Shi L Y，et al.Experiments on the drag－reduction mechanism of hydrophobic nanomaterial in reservoir microchannel［J］.Journal of Experiments in Fluid Mechanics，2010，24（6）：6－13.

［4］滕起，楊正明，劉學(xué)偉，等.特低滲透油藏平面滲流規(guī)律研究［J］.科學(xué)技術(shù)與工程，2012，12（36）：9820－9823.Teng Q，Yang Z M，Liu X W，et al.2Dseepage flow rule of the ultra－low permeability reservoir［J］.Science Technology and Engineering，2012，12（36）：9820－9823.

［5］Arns C H.A comparison of pore size distributions derived by NMR and X－ray－CT techniques［J］.Physica A：Statistical Mechanics and its Applications，2004，339（1－2）：159－165.

［6］徐祖新.基于CT掃描圖像的頁巖儲層非均質(zhì)性研究［J］.斷塊油氣藏，2014，26（6）：46－49.Xu Z X.Heterogeneity of shale reservoirs based on CT images［J］.Lithologic Reservoirs，2014，26（6）：46－49.

［7］陳濤平，宋文玲，高明，等.超低界面張力體系驅(qū)油方式與微觀驅(qū)油效果［J］.科學(xué)技術(shù)與工程，2010，（17）：4259－4262.Chen T P，Song W L，Gao M，et al.The research an displacement ways and microscopic displacement effects of ultralow interfacial tension system［J］.Science Technology and Engineering，2010，（17）：4259－4262.

［8］于明旭，朱維耀，宋洪慶.低滲透儲層可視化微觀滲流模型研制［J］.遼寧工程技術(shù)大學(xué)大學(xué)學(xué)報：自然科學(xué)版，2013，32（12）：1646－1650. Yu M X，Zhu W Y，Song H Q.Development of microscopic visualization flow model of low－permeability reservoir［J］.Journal of Liaoning Technical University（Natural Science），2013，32（12）：1646－1650.

［9］Sun W，Tang G Q.Visual study of water injection in low permeable sandstone［J］.Journal of Canadian Petroleum Technology，2006，45（11）：21－26.

［10］梁亞寧，張士誠，葉銀珠，等.CT掃描技術(shù)在低滲透砂巖巖心試驗中的應(yīng)用［J］.石油鉆采工藝，2011，33（2）：98－101.Liang Y N，Zhang S C，Ye Y Z，et al.Application of computed tomography scanning technology in low－permeability sandstone core test［J］.Oil Drilling and Production Technology，2011，33（2）：98－194101.

［11］Simjoo M，Dong Y，Andrianov A，et al.A CT scan study of immiscible foam flow in porous media for EOR［R］.SPE 155633，2012.

［12］Cnudde V，Boone M N.High－resolution X－ray computed tomography in geosciences：a review of the current technology and applications［J］.Earth Science Reviews，2013，23（4）：1－17.

［13］劉玉明，梁燦.復(fù)雜巖性儲層核磁共振測井巖心分析［J］.江漢石油學(xué)院學(xué)報，2014，36（6）：75－78.Liu Y M，Liang C.NMR core analysis in complex lithologic reservoir［J］.Journal of Oil and Gas Technology，2014，36（6）：75－78.

［14］Feng Y X，Chu X L，Xu Y P，et al.Corresponding factors influencing crude oils assay using low－field nuclear magnetic resonance［J］.China Petroleum Processing and Petrochemical Technology，2014，16（2）：34－39.

［15］Jiang L L，Song Y C，Liu Y，et al.Measurement of two phase flow in porous medium using high－resolution magnetic resonance imaging［J］.Chinese Journal of Chemical Engineering，2013，21（1）85－93.

［16］許長福，劉紅現(xiàn)，錢根寶，等.克拉瑪依礫巖儲集層微觀水驅(qū)油機理［J］.石油勘探與開發(fā)，2011，38（6）：725－732.Xu C F，Liu H X，Qian G B，et al.Microcosmic mechanisms of water－oil displacement in conglomerate reservoirs in Karamay Oilfield，NW China［J］.Petroleum Exploration and Development，2011，38（6）：725－732.

［17］曲巖濤，姜志敏，史京生，等.水驅(qū)油過程的核磁共振二維譜研究［J］.波譜學(xué)雜志，2012，29（1）：51－59.

［18］劉瑜，宋永臣，趙佳飛，等.多孔介質(zhì)斷面可視化與孔隙度測量的MRI方法［J］.測試技術(shù)學(xué)報，2011，25（03）：195－199.Liu Y，Song Y C，Zhao J F.Cross sectional visualization and porosity measurement of porous media by using MRI technique［J］.Journal of Test and Measurement Technology，2011，25（3）：195－199.

［19］郎東江，呂成遠，倫增珉，等.核磁共振成像巖心分析方法研究［J］.CT理論與應(yīng)用研究，2012，（01）：27－35.Lang D J，Lvy C Y，Lun Z M，et al.Study of method in core analysis by NMRI［J］.Computerized Tomography Theory and Applications，2012，21（1）：27－35.

Study of micro flow visualization with nuclear magnetic resonance in core

Di Qingfeng1，2，Hua Shuai1，2，Gu Chunyuan1，2，＊，Ye Feng1，2，Pang Dongshan1，2，Jiang Fan1，2，Yang Peiqiang3
（1.Shanghai Institute of Applied Mathematics and Mechanics，Shanghai University，Shanghai 200072，China；2.Shanghai Key Laboratory of Mechanics in Energy Engineering，Shanghai University，Shanghai 200072，China；3.Shanghai Niumag Corporation，Shanghai 200333，China）

Core micro flow visualization is an important method to study the microcosmic mechanism of the chemical flooding.The latest results of core micro flow visualization by MRI technology in our research group was introduced in detail.Some questions on low field MRI core displacement system（MRICDS）were put forward，and improvement and optimization were made from the two aspects of hardware and software to eliminate interference factors.The experiments of fluid image test during water－flooding process were successfully carried out using the improved MRICDS.The real－time signals of NMR and MRI of the oil－water distribution during the flow test were collected，and NMR－T2spectrums and oil or water distribution images were obtained in different stages.The images with high resolution can meet the request of visually analyzing the fluid distribution in the core.Experimental results show that residual oil decreases with PV increasing.The trendline of oil distribution drops rapidly first and then slowly.Images also show the phenomenon of residual oil gathering in the end area in core flow tests and its influence range is about 4mm from the end.A new method of oil saturation by the oil－water distribution map of MRI was introduced，and its results are consistent with that of the traditional method and the error is within 10%.This also shows that the reliability of the oil－water distribution map.This provides a new method for the study of oil saturation.Especially the other advantageof the method is that it can analyze the oil－water saturation at any local position.The results also indicate that the new method of MRI is worthy of further development for fluid distribution visualization during core micro flow.

micro flow visualization；low－field nuclear magnetic resonance；end effect of residual oil；oil or water distribution image；oil or water saturation

TE39

：A

（編輯：張巧蕓）

1672－9897（2016）03－0098－07

10.11729／syltlx20150107

2015－08－10；

2015－11－26

國家自然科學(xué)基金（51274136，50874071）；上海市科委重點項目（071605102）；上海高校創(chuàng)新團隊；上海市重點學(xué)科建設(shè)項目（530106）；上海市力學(xué)在能源工程中的重點實驗室和上海市教委重點項目（高峰學(xué)科建設(shè)項目）

＊通信作者E－mail：cygu＠shu.edu.cn

Di Q F，Hua S，Gu C Y，et al.Study of micro flow visualization with nuclear magnetic resonance in core.Journal of Experiments in Fluid Mechanics，2016，30（3）：98－103.狄勤豐，華 帥，顧春元，等.巖心微流動的核磁共振可視化研究.實驗流體力學(xué)，2016，30（3）：98－103.

狄勤豐（1963－），男，江蘇溧陽人，教授，博士生導(dǎo)師。研究方向：石油工程技術(shù)和力學(xué)問題研究。通信地址：上海市延長路149號189信箱（200072）.E－mail：qinfengd＠sina.com