許 傲，何宗斌，葉才駿，葉宇晗

（1.長(zhǎng)江大學(xué) 地球物理與石油資源學(xué)院，湖北 武漢 430100；2.長(zhǎng)江大學(xué) 油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430100）

0 引言

核磁共振技術(shù)自20 世紀(jì)80 年代問(wèn)世以來(lái)，被廣泛應(yīng)用于測(cè)井領(lǐng)域[1?3]，其將儀器測(cè)量到的信號(hào)擬合成回波衰減曲線，繼而反演得到橫向弛豫（T2）譜[4?5]，該T2譜可用來(lái)評(píng)價(jià)目的層的巖石物理參數(shù)，識(shí)別地層孔隙內(nèi)的流體類型[6?9]。

對(duì)于核磁共振T2譜的評(píng)價(jià)和預(yù)測(cè)，前人已做了一些研究。如白松濤等對(duì)T2譜進(jìn)行了定量表征的研究，對(duì)于儲(chǔ)層綜合分類等具有良好的評(píng)價(jià)效果[10]。鐘吉彬等人采用信號(hào)分析的方法對(duì)T2譜進(jìn)行了分解，對(duì)于復(fù)雜油水層的識(shí)別有較強(qiáng)的適應(yīng)性[11]。張家成等人利用連續(xù)小波變換和非對(duì)稱的高斯函數(shù)將T2譜分解成代表不同流體組分的分量譜，能夠準(zhǔn)確地識(shí)別孔隙內(nèi)不同賦存狀態(tài)的流體[12]。張哲等人通過(guò)長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)對(duì)核磁共振T2譜進(jìn)行了預(yù)測(cè)[13]。但是對(duì)于測(cè)井前的T2譜預(yù)測(cè)卻少有學(xué)者進(jìn)行研究，因此本文從測(cè)井前的T2譜預(yù)測(cè)著手，研究思路是在獲取地層的物性參數(shù)后，輸入仿真軟件，得到束縛水和可動(dòng)流體的T2信號(hào)；然后設(shè)置回波個(gè)數(shù)等參數(shù)生成回波，通過(guò)加噪模擬測(cè)量回波時(shí)產(chǎn)生的噪聲，反演得到T2譜。同時(shí)根據(jù)地層特性參數(shù)（如孔隙大小、地層溫度、壓力、油氣水黏度等）估算油氣水的T1值、T2值、地層最短完全極化時(shí)間、測(cè)量回波次數(shù)以及測(cè)量等待時(shí)間等，可以作為測(cè)量模式的參數(shù)參考。

以往的學(xué)者不僅沒(méi)有在測(cè)井前預(yù)測(cè)過(guò)T2譜的形態(tài)以及對(duì)應(yīng)的油氣水信號(hào)的大致位置，而且在核磁測(cè)井儀器測(cè)量時(shí)，并未考慮測(cè)前參數(shù)的合理設(shè)計(jì)，導(dǎo)致測(cè)量等待時(shí)間過(guò)長(zhǎng)或過(guò)短，影響了測(cè)井的開(kāi)發(fā)速率。選擇合適的測(cè)量參數(shù)能夠在不失精度的條件下，加速井下作業(yè)。在實(shí)際測(cè)井過(guò)程中，核磁資料的人力物力成本更高，但該軟件的優(yōu)勢(shì)在于通過(guò)其他測(cè)井資料得到的地層信息也可以作為輸入，能夠有效解決測(cè)前參數(shù)設(shè)置問(wèn)題以及預(yù)測(cè)地層的T2譜形態(tài)，從而對(duì)鄰近地層測(cè)井起到引導(dǎo)作用，這對(duì)于當(dāng)前地層的二次開(kāi)發(fā)也有著重要的應(yīng)用研究意義。

1 原 理

1.1 構(gòu)造T2譜信號(hào)

一維核磁T2譜由自旋回波脈沖序列反演得到，對(duì)于評(píng)價(jià)目的地層的巖石物理參數(shù)、識(shí)別地層孔隙流體起著至關(guān)重要的作用。構(gòu)造T2弛豫譜時(shí)，采用無(wú)序孔隙介質(zhì)的T2弛豫速率公式[14?15]進(jìn)行擬合，可寫(xiě)成如下形式：

式中：T2為橫向弛豫時(shí)間；T2S為表面弛豫時(shí)間，表征由自旋氫核與巖石孔壁碰撞引起能量衰減的過(guò)程；T2B為自由弛豫時(shí)間，其值取決于孔隙介質(zhì)本身；γ為氫核旋磁比；G為磁場(chǎng)梯度；TE為CPMG 脈沖序列回波間隔；D為流體的自擴(kuò)散系數(shù)；p為橫向表面弛豫率，一般砂巖的表面弛豫強(qiáng)度在1～10 μm/s；S為巖心表面積；V為孔隙體積；R為孔隙半徑；C為常數(shù)，C=1，2，3，分別用于平板模型、毛細(xì)管束模型和球狀模型；PPHI 為孔隙度；PERM 為滲透率。

由式（1）可知，巖石孔隙中流體的T2弛豫同時(shí)受三種弛豫過(guò)程的影響，相較可動(dòng)流體，束縛流體受孔隙表面的影響更大，因而該部分流體所受表面弛豫更強(qiáng)，弛豫時(shí)間更短。

自由弛豫T2B與流體本身及其流體所處的環(huán)境有關(guān)，為流體固有的弛豫特性。在實(shí)際地層中，由于測(cè)量時(shí)磁場(chǎng)的不均勻所帶來(lái)的影響不可忽略[16]，因此測(cè)前設(shè)計(jì)軟件考慮了非均勻磁場(chǎng)下引起的擴(kuò)散弛豫，引入了擴(kuò)散系數(shù)因子Dfact，可根據(jù)地層實(shí)際情況填寫(xiě)。各流體在巖石孔隙中的擴(kuò)散系數(shù)公式如下：

式中：Dfact為擴(kuò)散系數(shù)因子；Tk為溫度；press 為地層壓力；ηoil為油的黏度；pgas為氣體密度。

1.2 孔隙流體成分構(gòu)造

結(jié)合式（1）、式（2），可得到孔隙空間內(nèi)束縛流體和可動(dòng)流體的T2弛豫時(shí)間。由Timur/Coates 模型公式[17]可知，根據(jù)所給地層的滲透率和地層的沉積系數(shù)，可獲取到可動(dòng)流體與束縛流體的比值，公式如下：

式中：K為滲透率；A為滲透率乘積因子，與地層沉積等因素有關(guān)；?為地層孔隙度；m為孔隙度指數(shù)；FFI 為可動(dòng)流體體積；BVI 為束縛流體體積；n為可動(dòng)流體束縛流體比指數(shù)。

滲透率因子可反映不同地層的沉積情況，從而模擬不同的地層。已知某一地層的孔隙度后，結(jié)合式（3）可得各組成成分的孔隙大小，再結(jié)合式（1）可得各孔隙的橫向弛豫信號(hào)，如圖1 所示。

圖1 各流體組分信號(hào)大小

1.3 回波生成

為預(yù)測(cè)地層的T2譜，需由上述所得的流體信號(hào)正演得到衰減回波，再由回波反演得到T2譜?；夭ㄋp過(guò)程[18]公式為：

式中：M(t)為時(shí)刻t、第i個(gè)孔隙的信號(hào)；pi為孔隙；t=i·TE，TE為回波間隔；T2i為第i個(gè)孔隙的橫向弛豫時(shí)間。

回波串蘊(yùn)含了地層孔隙空間下各類流體的信息，通常大孔隙衰減速率更慢，小孔隙衰減速率更快，如圖2所示。

圖2 各孔隙流體的衰減速率

1.4 回波加噪

實(shí)際井下測(cè)量時(shí)，無(wú)法避免噪聲的影響，因此采用Box?Muller transformation 算法模擬井下接收到共振信號(hào)時(shí)產(chǎn)生的噪聲，將均勻分布的隨機(jī)數(shù)轉(zhuǎn)換成具有特定信噪比（SNR）的高斯分布噪聲，公式如下：

式中：rand 為符合高斯分布的隨機(jī)數(shù)；U1、U2為均勻分布的隨機(jī)數(shù)。

為使均勻分布的噪聲符合特定的信噪比，需對(duì)噪聲的分布進(jìn)行改進(jìn)。采用信噪比公式構(gòu)造，公式為：

式中：SNR 為信噪比；Ps為原始信號(hào)功率；Pn為噪聲信號(hào)功率；Si為第i個(gè)信號(hào)；ni為第i個(gè)噪聲；N為有效信號(hào)或噪聲的個(gè)數(shù)。

1.5 回波擬合方法

回波擬合方法將得到的孔隙組分以及T2弛豫信號(hào)反演成T2譜，實(shí)現(xiàn)了SVD（奇異值分解）反演算法，方便與井下所得到的T2譜進(jìn)行比較；同時(shí)預(yù)留了反演算法的擴(kuò)展接口，方便開(kāi)發(fā)人員擴(kuò)展新的反演算法。

2 軟件設(shè)計(jì)與實(shí)現(xiàn)

本文模擬軟件基于.NET4.8 框架，使用C#語(yǔ)言在Visual Studio 2022 環(huán)境下開(kāi)發(fā)。

2.1 結(jié)構(gòu)設(shè)計(jì)

軟件設(shè)計(jì)結(jié)構(gòu)整體分為算法、應(yīng)用模塊、數(shù)據(jù)及圖形顯示三大部分，如圖3 所示。

圖3 軟件結(jié)構(gòu)設(shè)計(jì)

2.2 功能實(shí)現(xiàn)

2.2.1T2譜信號(hào)構(gòu)造模塊

由式（1）、式（2）計(jì)算可得黏土束縛水、毛管束縛水、可動(dòng)水以及烴的弛豫時(shí)間。由式（3）計(jì)算可得各組分流體的孔隙大小。主要預(yù)留的參數(shù)包括流體類型選擇、地層孔隙度、滲透率、含油飽和度、氣體含氫指數(shù)、地層壓力和溫度、泥漿濾液侵入程度、磁場(chǎng)梯度、原油黏度、氣體密度以及受限擴(kuò)散系數(shù)。T2譜信號(hào)構(gòu)造模塊參數(shù)設(shè)置如圖4 所示，允許多種流體，包括泥漿濾液的綜合弛豫過(guò)程。

圖4 T2譜構(gòu)造參數(shù)對(duì)話框

用戶可以根據(jù)實(shí)際地層測(cè)量到的信息，或在缺少核磁資料的情況下，選用其他測(cè)井資料得到的地層孔隙和滲透率等物性參數(shù)作為測(cè)前設(shè)計(jì)軟件的輸入，內(nèi)部算法會(huì)實(shí)時(shí)根據(jù)用戶所提供的地層信息，計(jì)算生成目的層各流體的T1、T2弛豫時(shí)間，以及目的層各組分流體的孔隙大小，從而得到T2譜各組分信號(hào)。數(shù)值模擬結(jié)果如圖5所示。

圖5 各流體組分信號(hào)大小

2.2.2 回波生成及加噪模塊

回波生成模塊主要根據(jù)式（4）進(jìn)行編寫(xiě)，主要預(yù)留的參數(shù)為：回波間隔TE、回波個(gè)數(shù)NE、T2布點(diǎn)個(gè)數(shù)及其范圍，一般選用0.3～3 000 ms。參數(shù)設(shè)置如圖6 所示。NE的大小可由儀器參數(shù)建議模塊給出的NE或?qū)嶋H地層測(cè)量時(shí)儀器的NE大小來(lái)設(shè)置。

圖6 回波參數(shù)設(shè)置

T2譜信號(hào)構(gòu)造完后，點(diǎn)擊Echo 按鈕，生成的回波如圖7 所示。

圖7 回波數(shù)據(jù)

通過(guò)修改信噪比（SNR）對(duì)原始回波數(shù)據(jù)進(jìn)行加噪，模擬井下測(cè)量回波信號(hào)時(shí)伴隨產(chǎn)生的噪聲。加噪采用均勻分布的特定信噪比的高斯白噪聲，由式（5）和式（6）構(gòu)造得出：用戶可隨時(shí)更改信噪比，多次模擬，加噪后的回波信號(hào)如圖8 所示。

圖8 加噪后的回波信號(hào)

2.2.3 反演模塊

采用SVD 反演算法，將特定信噪比的回波數(shù)據(jù)反演成T2譜，預(yù)留的參數(shù)設(shè)置為T(mén)2布點(diǎn)個(gè)數(shù)、T2布點(diǎn)范圍以及平滑因子。平滑因子可使T2譜曲線更光滑，其數(shù)值越大，譜越光滑。T2譜構(gòu)造模塊示例如圖9 所示。

圖9 T2譜構(gòu)造模塊示例

2.2.4 儀器參數(shù)建議模塊

儀器參數(shù)建議模塊如圖10 所示。圖中：Result 欄為參數(shù)計(jì)算結(jié)果，展示了各組分的弛豫時(shí)間；Suggestion欄根據(jù)儀器選擇的測(cè)井模式和回波衰減的速率，結(jié)合各組分的弛豫時(shí)間得到TE和Tw以及回波個(gè)數(shù)的設(shè)置建議，對(duì)同一地層或相鄰地層再次測(cè)量時(shí)，能夠指導(dǎo)井下工作。

圖10 儀器參數(shù)建議模塊

3 實(shí)驗(yàn)比對(duì)

3.1 實(shí)驗(yàn)步驟

本文實(shí)驗(yàn)對(duì)比流程如圖11 所示，具體步驟如下：

圖11 實(shí)驗(yàn)對(duì)比流程

1）利用巖石物理實(shí)驗(yàn)，得到巖芯的核磁孔隙度以及滲透率等物性參數(shù)。

2）將步驟1）所得參數(shù)輸入測(cè)前設(shè)計(jì)軟件，得到回波與T2譜，記作A 組。

3）實(shí)際測(cè)量該地區(qū)的井下數(shù)據(jù)，計(jì)算得到回波與T2譜，記作B 組。

4）比對(duì)A、B 兩組，輸入測(cè)前設(shè)計(jì)軟件的其他物性參數(shù)，使得A 組的回波衰減與反演后的T2譜接近于B 組。

5）將A組所得物性參數(shù)視作鄰近地層的預(yù)測(cè)物性參數(shù)，可驗(yàn)證鄰近地層的油水信號(hào)及其分布的弛豫時(shí)間。

3.2 實(shí)驗(yàn)結(jié)果分析

將孔隙度（12.99）、滲透率（6.97 mD）、回波間隔（0.4 ms）、回波個(gè)數(shù)（2 048）、磁場(chǎng)梯度（13.5 gs/cm）等參數(shù)輸入至該軟件，經(jīng)過(guò)計(jì)算生成的弛豫信號(hào)的位置以及各孔隙組成成分如表1 所示。

表1 弛豫信號(hào)位置及各孔隙組成成分

圖12 所示為該深度段利用測(cè)前設(shè)計(jì)軟件得到的T2譜數(shù)據(jù)（實(shí)線）和實(shí)際井下測(cè)量得到的T2譜數(shù)據(jù)（虛線）比對(duì)?？梢钥闯鰷y(cè)前設(shè)計(jì)軟件預(yù)測(cè)得到的T2譜與井下T2譜的峰值位置、包絡(luò)線面積及形態(tài)都非常接近，進(jìn)一步驗(yàn)證了本文研制仿真軟件的可靠性。圖13 為深度在3 762.5～3 772.5 段的測(cè)井解釋成果圖，其中測(cè)前設(shè)計(jì)軟件預(yù)測(cè)的T2譜與井下T2譜比對(duì)，可以看到與實(shí)驗(yàn)結(jié)果十分接近，再次驗(yàn)證了測(cè)前設(shè)計(jì)軟件的可靠性。

圖12 T2譜對(duì)比結(jié)果

圖13 測(cè)前設(shè)計(jì)軟件預(yù)測(cè)的T2譜與井下T2譜比對(duì)

4 結(jié)論

本文研發(fā)的測(cè)前設(shè)計(jì)正演軟件能夠分析地層的油氣水信號(hào)及其橫縱向弛豫時(shí)間，得到一組測(cè)井模式的儀器參數(shù)建議，提高井下測(cè)量的速度，還可以針對(duì)儲(chǔ)層預(yù)測(cè)其T2譜，便于與實(shí)際井下測(cè)量所得的T2譜做參照，尋找誤差。通過(guò)改變回波間隔、等待時(shí)間以及回波個(gè)數(shù)等參數(shù)，所設(shè)計(jì)軟件能夠模擬多種地層下的T2譜，對(duì)核磁共振測(cè)井的正反演算法的改進(jìn)有重要意義。