吳　飛, 范宜仁, 王　帥, 鄧少貴, 邢東輝, 巫振觀, 楊培強(qiáng)

(1.中國(guó)石油大學(xué)地球科學(xué)與技術(shù)學(xué)院,山東青島 266580; 2.中國(guó)石油大學(xué)CNPC測(cè)井重點(diǎn)實(shí)驗(yàn)室,山東青島 266580;3.大港油田第五采油廠地質(zhì)研究所,天津 300280; 4.上海紐邁電子科技有限公司, 上海 200333)

一D-T2二維核磁共振脈沖序列改進(jìn)設(shè)計(jì)及性能對(duì)比

吳飛1,2, 范宜仁1,2, 王帥3, 鄧少貴1,2, 邢東輝1,2, 巫振觀1,2, 楊培強(qiáng)4

(1.中國(guó)石油大學(xué)地球科學(xué)與技術(shù)學(xué)院,山東青島 266580; 2.中國(guó)石油大學(xué)CNPC測(cè)井重點(diǎn)實(shí)驗(yàn)室,山東青島 266580;3.大港油田第五采油廠地質(zhì)研究所,天津 300280; 4.上海紐邁電子科技有限公司, 上海 200333)

針對(duì)現(xiàn)有D-T2二維核磁共振脈沖序列無法兼顧擴(kuò)散系數(shù)測(cè)量范圍和橫向弛豫分辨率的難題,在對(duì)比PFG(脈沖梯度場(chǎng))、STE-PFG(stimulated echo PFG)、BP-PFG(Bi-polar PFG)、改良式CPMG(Car-Purcel-Meiboom-Gill)、“擴(kuò)散編程”、多回波間隔CPMG脈沖序列技術(shù)特點(diǎn)的基礎(chǔ)上,提出一種基于脈沖梯度場(chǎng)的兩窗口、雙變量D-T2改進(jìn)脈沖序列。PFG、改良式CPMG、D-T2改進(jìn)脈沖序列在氣水、油水、稠油模型不同信噪比條件下的正反演數(shù)值模擬結(jié)果表明,D-T2改進(jìn)脈沖序列達(dá)到了平衡擴(kuò)散系數(shù)測(cè)量范圍和橫向弛豫分辨率的設(shè)計(jì)要求,在非常規(guī)油氣儲(chǔ)層的探測(cè)和測(cè)井評(píng)價(jià)中具有廣泛應(yīng)用前景,可為推動(dòng)D-T2二維核磁共振技術(shù)的發(fā)展提供有利條件。

二維核磁共振; 擴(kuò)散系數(shù); 橫向弛豫分辨率; 脈沖梯度場(chǎng); 雙變量; 改進(jìn)脈沖序列; 數(shù)值模擬

2002年,Sun等[1-4]為彌補(bǔ)一維T2譜油氣水信號(hào)重疊的短板,將核磁共振波譜學(xué)中的二維核磁共振概念應(yīng)用到石油測(cè)井領(lǐng)域,開創(chuàng)了以D-T2、T1-T2為典型代表的二維核磁共振測(cè)井技術(shù)[5]。然而隨著勘探對(duì)象的日益復(fù)雜,現(xiàn)有的D-T2二維核磁共振技術(shù)無法兼顧擴(kuò)散系數(shù)測(cè)量范圍和橫向弛豫分辨率[6-10]。在對(duì)比六種常用D-T2脈沖序列的基礎(chǔ)上,筆者提出一種基于脈沖梯度場(chǎng)的兩窗口、雙變量D-T2改進(jìn)脈沖序列,并構(gòu)造不同儲(chǔ)層流體類型的D-T2分布模型,結(jié)合正反演數(shù)值模擬,對(duì)比分析PFG、改良式CPMG、D-T2改進(jìn)脈沖序列對(duì)天然氣、中等黏度原油、稠油的識(shí)別能力,考察D-T2改進(jìn)脈沖序列的綜合性能。

1　D-T2基本原理

由于油、氣、水?dāng)U散系數(shù)存在明顯差異(一般Dg>Dw>Do),D-T2二維分布能夠有效區(qū)分油、氣、水,如圖1所示。當(dāng)核磁測(cè)量的等待時(shí)間TW足夠長(zhǎng)時(shí),梯度場(chǎng)下自旋回波的幅度可表示為

(1)

式中,b為回波幅度;T2為孔隙流體的固有弛豫(自由弛豫+表面弛豫),本文中T2如無特殊說明均指固有弛豫;f(D,T2)為氫核在(D,T2)二維空間的分布;kT為與橫向弛豫時(shí)間相關(guān)的核函數(shù);t為回波時(shí)間;kD為與擴(kuò)散系數(shù)相關(guān)的核函數(shù);X為加載擴(kuò)散弛豫的變量。

圖1　油氣水D-T2分布示意圖Fig.1　Petroleum, gas and water D-T2 distribution diagram

2　D-T2脈沖序列改進(jìn)設(shè)計(jì)

2.1D-T2脈沖序列發(fā)展現(xiàn)狀

2.1.1PFG脈沖序列

射頻場(chǎng)采集時(shí)序的第一個(gè)窗口持續(xù)時(shí)間為t0,在該窗口內(nèi)180°脈沖兩側(cè)施加一組對(duì)稱的脈沖梯度,用于加載流體擴(kuò)散弛豫信息;第二個(gè)窗口用最短回波間隔采集CPMG回波串,采集孔隙流體橫向弛豫信息,如圖2(a)所示[11-12]。PFG脈沖序列的回波幅度為

(2)

式中,bik表示脈沖梯度為Gk時(shí)第i個(gè)回波的幅度;f(Dp,T2j)為擴(kuò)散系數(shù)Dp、橫向弛豫時(shí)間T2j對(duì)應(yīng)的孔隙度分量;γ為氫核的旋磁比;δ為梯度脈沖持續(xù)時(shí)間;Δ為兩個(gè)梯度脈沖間的間隔;TE為后續(xù)CPMG序列的回波間隔。改變脈沖輸出梯度大小,采集多條CPMG回波串,解譜時(shí)將采集的回波串簇利用式(2)聯(lián)合反演[13-14],即可得到儲(chǔ)層巖石孔隙流體的D-T2分布,國(guó)產(chǎn)紐邁核磁共振巖心分析儀MR-DF的D-T2采集就是采用該序列。

2.1.2STE-PFG脈沖序列

將PFG脈沖序列的第一個(gè)180°脈沖替換為兩個(gè)90°脈沖,對(duì)稱脈沖梯度施加在第一個(gè)和第二個(gè)90°脈沖、第三個(gè)90°脈沖和第一個(gè)180°脈沖之間,如圖2(b)所示[15]。由于第二個(gè)和第三個(gè)90°脈沖之間的核磁信號(hào)衰減為縱向弛豫,則STE-PFG脈沖序列的回波幅度為

(3)

式中,f(T1q,Dp,T2j)為縱向弛豫時(shí)間T1q、擴(kuò)散系數(shù)Dp、橫向弛豫時(shí)間T2j對(duì)應(yīng)的孔隙度分量;t1為第二個(gè)90°和第三個(gè)90°脈沖之間的間隔,即縱向弛豫的時(shí)間;t0為第一個(gè)窗口中橫向弛豫的時(shí)間。

當(dāng)測(cè)量對(duì)象滿足T1?T2時(shí),式(3)中的縱向弛豫因子exp(-t1/T1)≈1,此時(shí)STE-PFG脈沖序列的回波信號(hào)衰減可表示為式(2);改變脈沖輸出梯度大小,采集多條CPMG回波串,解譜時(shí)將采集的回波串簇利用式(2)聯(lián)合反演,即可得到儲(chǔ)層巖石孔隙流體的D-T2分布,牛津巖心分析儀GeoSpec2的D-T2采集就是采用該序列。

當(dāng)測(cè)量對(duì)象中含有短弛豫組分(T2

圖2　六種常見D-T2脈沖序列示意圖Fig.2　Schematic diagrams of six D-T2 pulse sequences

2.1.3BP-PFG脈沖序列

在STE-PFG脈沖序列基礎(chǔ)上,將第一個(gè)窗口中橫向弛豫的時(shí)間進(jìn)一步劃分,脈沖梯度也增加到兩對(duì),并且兩對(duì)梯度脈沖的方向相反,其余的與圖2(b)的STE-PFG序列相同,如圖2(c)所示[16-18]。曲巖濤等[19]在自主研制的低場(chǎng)核磁巖心分析儀上利用BP-PFG脈沖序列研究了水驅(qū)油過程中的油水分布規(guī)律。

BP-PFG脈沖序列的雙梯度設(shè)計(jì)主要用于克服儲(chǔ)層巖石的高內(nèi)部磁場(chǎng)梯度影響,考慮內(nèi)部磁場(chǎng)梯度g0時(shí),BP-PFG脈沖序列的回波幅度為

i=1,2,3,….

(4)

式中,δ1為脈沖梯度距離左側(cè)緊鄰扳轉(zhuǎn)脈沖的時(shí)間間隔;δ2為脈沖梯度距離右側(cè)緊鄰扳轉(zhuǎn)脈沖的時(shí)間間隔;g0為巖石內(nèi)部磁場(chǎng)梯度。

考慮內(nèi)部磁場(chǎng)梯度g0時(shí),STE-PFG脈沖序列的回波幅度為

(5)

對(duì)比式(3)～(5)可知,由于內(nèi)部磁場(chǎng)梯度的作用,回波串幅度衰減增加了兩個(gè)與g0有關(guān)的增強(qiáng)弛豫因子E(Gg0)、E(g02)。在t0、t1、G、δ相同時(shí),BP-PFG脈沖序列通過設(shè)置δ1=δ2,使E(Gg0)=1,與之相對(duì)的是,STE-PFG脈沖序列的E(Gg0)<1,并且STE-PFG的E(g02)小于BP-PFG的E(g02),因此BP-PFG的雙脈沖設(shè)計(jì)有效降低了巖石內(nèi)部磁場(chǎng)梯度對(duì)D-T2數(shù)據(jù)采集的影響。

由于BP-PFG脈沖序列的獨(dú)特設(shè)計(jì),反演時(shí)可忽略內(nèi)部磁場(chǎng)梯度的影響,其回波幅度為

(6)

改變脈沖輸出梯度,采集多條CPMG回波串,解譜時(shí)將采集的回波串簇利用式(6)聯(lián)合反演,即可得到儲(chǔ)層巖石孔隙流體的D-T2分布。

2.1.4改良式CPMG脈沖序列

每條CPMG序列在時(shí)間軸上分為兩個(gè)窗口,第一個(gè)窗口長(zhǎng)度固定為t0,改變第一個(gè)窗口中的回波個(gè)數(shù)NE1使回波間隔從大變小,加載儲(chǔ)層孔隙流體擴(kuò)散弛豫信息;第二個(gè)窗口用儀器的最短回波間隔采集CPMG回波信號(hào),將擴(kuò)散弛豫影響降到最小,采集儲(chǔ)層孔隙流體橫向弛豫信息,如圖2(d)所示,該序列最早用于儲(chǔ)層巖石內(nèi)部磁場(chǎng)梯度的研究。恒定梯度場(chǎng)下,改良式CPMG脈沖序列的回波幅度為

(7)

式中,bik代表第一個(gè)窗口的回波個(gè)數(shù)為NE1k時(shí)第i個(gè)回波的回波幅度。改變第一個(gè)窗口的回波個(gè)數(shù),采集多條CPMG回波串,解譜時(shí)將采集的回波串簇利用式(7)聯(lián)合反演,即可得到儲(chǔ)層巖石孔隙流體的D-T2分布。

2.1.5“擴(kuò)散編程”脈沖序列

(8)

改變第一個(gè)窗口的回波間隔,采集多條CPMG回波串,解譜時(shí)將采集的回波串簇利用式(8)聯(lián)合反演,即可得到儲(chǔ)層巖石孔隙流體的D-T2分布。MR Scanner的D-T2二維核磁測(cè)井采用“擴(kuò)散編程”脈沖序列,并取得了不錯(cuò)的應(yīng)用效果。

2.1.6多回波間隔CPMG脈沖序列

利用常規(guī)CPMG脈沖序列,無需重新設(shè)計(jì)脈沖序列,只需改變回波間隔采集多條自旋回波串,如圖2(f)所示。恒定梯度場(chǎng)下,多回波間隔CPMG脈沖序列的回波幅度為

(9)

變化回波間隔,采集多條CPMG回波串,實(shí)現(xiàn)D-T2二維數(shù)據(jù)采集,解譜時(shí)將采集的回波串簇利用式(9)聯(lián)立反演[20-23],即可得到儲(chǔ)層巖石孔隙流體的D-T2分布。MREX的PP OIL、PP HEAVY OIL觀測(cè)模式、MRIL新增的四TE觀測(cè)模式均采用多回波間隔CPMG脈沖序列,沿用了成熟的CPMG序列,降低了D-T2核磁測(cè)井儀器的研發(fā)成本。

2.2D-T2現(xiàn)有脈沖序列性能對(duì)比

上述六種常見D-T2脈沖序列的基本信息及性能對(duì)比見表1。脈沖梯度場(chǎng)得益于高磁場(chǎng)梯度值,其D-T2脈沖序列具有較大的擴(kuò)散系數(shù)測(cè)量范圍,但是這類脈沖序列在脈沖梯度作用窗口內(nèi)沒有自旋回波信號(hào),導(dǎo)致其對(duì)短弛豫組分的分辨率較低。與之相對(duì)的是,恒定梯度場(chǎng)梯度值一般較小,其脈沖序列的擴(kuò)散弛豫加載能力不如脈沖梯度場(chǎng),這也是目前核磁共振測(cè)井儀器測(cè)遇稠油層時(shí)D-T2譜中D軸產(chǎn)生“拖尾”現(xiàn)象的根源,但是這類脈沖序列在擴(kuò)散弛豫加載窗口內(nèi)采集自旋回波信號(hào),或?qū)U(kuò)散弛豫信息加載到整條CPMG回波串上,因而這類脈沖序列的橫向弛豫分辨能力更強(qiáng)。

表1　六種常見D-T2脈沖序列的基本信息及性能對(duì)比

2.3D-T2脈沖序列改進(jìn)設(shè)計(jì)

本文中提出一種基于脈沖梯度場(chǎng)的兩窗口、雙變量D-T2改進(jìn)脈沖序列,如圖3所示。

圖3　脈沖梯度場(chǎng)D-T2改進(jìn)脈沖序列示意圖Fig.3　Schematic diagram of modified D-T2pulse sequence on pulsed field gradient

圖3中D-T2改進(jìn)脈沖序列的射頻采集時(shí)序與改良式CPMG序列相同,時(shí)間軸也分為兩個(gè)窗口,第一個(gè)窗口長(zhǎng)度固定為t0,并在該窗口內(nèi)采集自旋回波信號(hào),為D-T2反演譜中短弛豫組分(T2

由多孔介質(zhì)核磁共振弛豫理論,可得出圖3所示D-T2改進(jìn)脈沖序列的回波幅度為

(10)

式中,bikq代表第一個(gè)窗口的回波個(gè)數(shù)為NE1k、脈沖梯度為Gq時(shí)第i個(gè)回波的回波幅度;t0為第一個(gè)窗口的持續(xù)時(shí)間;TE為第二個(gè)窗口的回波幅度。改變第一個(gè)窗口內(nèi)的脈沖梯度值和回波間隔,采集多條CPMG回波串,解譜時(shí)將采集的回波串簇利用式(10)聯(lián)合反演,即可得到儲(chǔ)層孔隙流體的D-T2二維分布。

3　D-T2改進(jìn)脈沖序列性能對(duì)比

利用D-T2二維核磁共振正反演數(shù)值模擬綜合評(píng)估本文中D-T2改進(jìn)脈沖序列的性能,主要步驟:構(gòu)建不同儲(chǔ)層流體D-T2分布模型,設(shè)置采集參數(shù)正演D-T2回波數(shù)據(jù),回波數(shù)據(jù)加噪聲后使用截?cái)嗥娈愔捣纸夥ń庾V,對(duì)比構(gòu)造譜與反演譜并計(jì)算相對(duì)誤差。構(gòu)造氣水、油水、稠油D-T2模型,如圖4所示。對(duì)PFG、改良式CPMG、D-T2改進(jìn)脈沖序列進(jìn)行正反演數(shù)值模擬。PFG脈沖序列基本參數(shù):Gmax=0.8T/m,Gkmin=0%,Gkmax=100%,脈沖梯度最小間隔1%,TEmin=0.1ms。改良式CPMG脈沖序列基本參數(shù):G=0.4T/m,TEmin=0.1ms。D-T2改進(jìn)脈沖序列基本參數(shù):Gmax=0.8T/m,Gkmin=0%,Gkmax=100%,脈沖梯度最小間隔1%,TEmin=0.1ms。圖4所示模型的D-T2回波串見圖5,回波數(shù)據(jù)加噪聲后的反演譜見圖6～8(圖中T2譜、D譜的黑線是構(gòu)造譜,紅線是反演譜),不同信噪比條件下D-T2反演譜的相對(duì)誤差及含水飽和度的相對(duì)誤差見圖9和圖10。

圖4　不同儲(chǔ)層流體D-T2構(gòu)造模型Fig.4　D-T2structural models of different formation fluid

圖6　PFG脈沖序列D-T2反演譜Fig.6　D-T2inversion spectrums of PFG pulse sequence

圖7　改良式CPMG脈沖序列D-T2反演譜Fig.7　D-T2 inversion spectrums of modified CPMG pulse sequence

圖8　D-T2改進(jìn)脈沖序列D-T2反演譜Fig.8　D-T2inversion spectrums of modified D-T2pulse sequence

圖9　三種脈沖序列D-T2反演譜的相對(duì)誤差Fig.9　Relative error of inversion spectrums in three pulse sequences

由上述數(shù)值模擬結(jié)果可知,與改良式CPMG脈沖序列相比,D-T2改進(jìn)脈沖序列保持了脈沖梯度場(chǎng)高梯度值的優(yōu)點(diǎn),即使流體(天然氣、中等黏度原油、稠油)的擴(kuò)散系數(shù)跨越4個(gè)數(shù)量級(jí),D-T2改進(jìn)脈沖序列仍然能夠有效識(shí)別流體類型,對(duì)稠油的識(shí)別能力明顯優(yōu)于改良式CPMG脈沖序列;與PFG脈沖序列相比,D-T2改進(jìn)脈沖序列保持了恒定梯度場(chǎng)回波采集不間斷的優(yōu)點(diǎn),其橫向弛豫分辨率優(yōu)于PFG脈沖序列;同時(shí),D-T2改進(jìn)脈沖序列的雙變量(第一個(gè)窗口的脈沖梯度、回波間隔)設(shè)計(jì)能更好地適應(yīng)儲(chǔ)層巖石孔隙流體性質(zhì)(擴(kuò)散系數(shù)、橫向弛豫時(shí)間)的變化。

圖10　三種脈沖序列D-T2反演譜含水飽和度相對(duì)誤差Fig.10　Relative error of calculated saturation from inversion spectrums in three pulse sequences

4　結(jié)束語

提出一種基于脈沖梯度場(chǎng)的兩窗口、雙變量D-T2改進(jìn)脈沖序列。氣水模型、油水模型、稠油模型的數(shù)值模擬結(jié)果初步驗(yàn)證了D-T2改進(jìn)脈沖序列的可行性；與PFG、改良式CPMG脈沖序列的數(shù)值模擬對(duì)比結(jié)果表明，D-T2改進(jìn)脈沖序列達(dá)到了平衡擴(kuò)散系數(shù)測(cè)量范圍和橫向弛豫分辨率的設(shè)計(jì)目標(biāo)；D-T2改進(jìn)脈沖序列的雙變量設(shè)計(jì)使采集參數(shù)的設(shè)置更靈活，能更好地適應(yīng)儲(chǔ)層孔隙流體核磁共振性質(zhì)的變化，但是雙變量設(shè)計(jì)也增加了D-T2測(cè)前設(shè)計(jì)的工作量及難度，下一步研究重點(diǎn)是借助數(shù)值模擬摸索D-T2改進(jìn)脈沖序列的采集參數(shù)調(diào)節(jié)規(guī)律。

[1]SUN B Q, DUNN K J. Probing the internal field gradients of porous media [J]. Phys Rev E, 2002,65(5):051309.

[2]HURLIMANN M D, VENKATARAMANAN L. Quantitative measurement of two dimensional distribution functions of diffusion and relaxation in grossly inhomogeneous fields [J]. Journal of Magnetic Resonance, 2002,157(1):31-42.

[3]HURLIMANN M D, VENKATARAMANAN L, FLAUM C. The diffusion-spin relaxation time distribution function as an experimental probe to characterize fluid mixtures in porous media [J]. Journal of Chemical Physics, 2002,117(22):10223-10232.

[4]SONG Y Q, VENKATARAMANAN L, HURLIMANN M D, et al.T1-T2correlation spectra obtained using a fast two-dimensional Laplace inversion[J]. Journal of Magnetic Resonance, 2002,154(2):261-268.

[5]鄧克俊. 核磁共振測(cè)井理論及應(yīng)用[M].東營(yíng): 中國(guó)石油大學(xué)出版社, 2010:9-148.

[6]HURLIMANN M D, FLAUM M, VENKATARAMANAN L, et al. Diffusion-relaxation distribution functions of sedimentary rocks in different saturation states [J]. Magnetic Resonance Imaging, 2003,21(3):305-310.

[7]SUN B Q, DUNN K J. Two-dimensional nuclear magnetic resonance petrophysics [J]. Magnetic Resonance Imaging, 2005,23(2):259-262.

[8]謝然紅,肖立志. (T2,D)二維核磁共振測(cè)井識(shí)別儲(chǔ)層流體的方法[J]. 地球物理學(xué)報(bào), 2009, 52(9):2410-2418. XIE Ranhong, XIAO Lizhi. The (T2,D) NMR logging method for fluids characterization [J].Chinese Journal of Geophysics, 2009,52(9):2410-2418.

[9]譚茂金, 鄒友龍. (T2,D)二維核磁共振測(cè)井混合反演方法與參數(shù)影響分析[J]. 地球物理學(xué)報(bào), 2012,55(2):683-692. TAN Maojin, ZOU Youlong. A hybrid inversion method of (T2,D) 2D NMR logging and observation parameters effects [J]. Chinese Journal Geophysics, 2012,55(2):683-629.

[10]胡法龍,周燦燦,李潮流, 等. 基于弛豫-擴(kuò)散的二維核磁共振流體識(shí)別方法[J]. 石油勘探與開發(fā), 2012,39(5):552-558. HU Falong, ZHOU C C, LI Chaoliu, et al. Fluid identification method based on 2D diffusion-relaxation NMR [J].Petroleum Exploration & Development, 2012,39(5):552-558.

[11]STEJSKAL E O, TANNER J E. Spin diffusion measurement: spin echoes in the presence of a time-dependent field gradient [J]. The Journal of Chemical Physics, 1965,42(1):288-292.

[12]STEJSKAL E O. Use of spin echoes in a pulsed magnetic-field gradient to study anisotropic, restricted diffusion and flow [J]. The Journal of Chemical Physics, 1965,43(10):3597-3603.

[13]王忠東,肖立志,劉堂宴. 核磁共振弛豫信號(hào)多指數(shù)反演新方法及其應(yīng)用[J]. 中國(guó)科學(xué):G輯, 2003,33(4):323-332. WANG Zhongdong, XIAO Lizhi, LIU Tangyan. A new multi-exponential inversion method of NMR relaxation signal and its application [J].Science in China (Series G), 2003,33(4):323-332.

[14]VENKATARAMANAN L, SONG Y Q, HURLIMANN M D. Solving Fredholm integrals of the first kind with tensor product structure in 2 and 2.5 dimensions [J]. IEEE Transactions on Signal Processing, 2002,50(5):1017-1026.

[15]TANNER J E. Use of the stimulated echo in NMR diffusion studies [J]. The Journal of Chemical Physics, 1970,52(5):2523-2526.

[16]KARLICEK J R F, LOWE I J. A modified pulsed gradient technique for measuring diffusion in the presence of large background gradient[J]. Journal of Magnetic Resonance, 1980,37(1):75-91.

[17]COTTS R M, HOCH M J R, SUN T, MARKERT J T. Pulsed field gradient stimulated echo methods for improved NMR diffusion measurements in heterogeneous systems[J].Journal of Magnetic Resonance, 1989,83:252-266.

[18]VANDUSSCHOTEN D.Flexible PFG NMR desensitized for susceptibility artifacts,using the PFG multiple-spin-echo sequence[J]. Journal of Magnetic Resonance, Series A, 1995,112(2):237-240.

[19]曲巖濤,姜志敏,史京生, 等. 水驅(qū)油過程的核磁共振二維譜研究[J]. 波譜學(xué)雜志, 2012,29(1):51-59. QU Yantao, JIANG Zhimin, SHI Jingsheng, et al. Water flooding studied by measuring NMR two-dimensional distribution function [J]. Chinese Journal of Magnetic Resonance, 2012,29(1):51-59.

[20]SUN B Q, DUNN K J.A global inversion method of multi-dimensional NMR logging[J].Journal of Magnetic Resonance, 2005,172(1):152-160.[21]謝然紅,肖立志,劉家軍, 等. 核磁共振多回波串聯(lián)合反演方法[J]. 地球物理學(xué)報(bào), 2009,52(11):2913-2919. XIE Ranhong, XIAO Lizhi, LIU Jiajun, et al. A method for multiple echo trains jointing inversion of NMR relaxation measurements[J].Chinese Journal of Geophysics, 2009,52(11):2913-2919.

[22]TAN M J, ZOU Y L, ZHOU C C. A new inversion method for (T2,D) 2D NMR logging and fluid typing[J].Computers & Geosciences, 2013,51(2):366-380.

[23]高陽,肖立志,謝慶明. 基于BG理論確定Phillips方法的正則化參數(shù)并評(píng)價(jià)反演結(jié)果[J]. 中國(guó)石油大學(xué)學(xué)報(bào):自然科學(xué)版,2014,38(3):50-56. GAO Yang, XIAO Lizhi, XIE Qingming. Evaluating regular parameter of Phillips method and assessing inversions based on BG theory[J]. Journal of China University of Petroleum (Edition of Natural Science), 2014,38(3):50-56.

(編輯修榮榮)

An improved design ofD-T22D NMR pulse sequence and performance comparison

WU Fei1,2, FAN Yiren1,2, WANG Shuai3, DENG Shaogui1,2,ING Donghui1,2, WU Zhenguan1,2, YANG Peiqiang4

(1.SchoolofGeosciencesinChinaUniversityofPetroleum,Qingdao266580,China；2.CNPCKeyLaboratoryforWellLogginginChinaUniversityofPetroleum,Qingdao266580,China；3.GeologyInstituteoftheFifthOilProductionPlantinDagangOilfield,Tianjin300280,China；4.ShanghaiNiumaiElectronicTechnologyCompanyLimited,Shanghai200333,China)

In order to solve the difficulties inD-T22D NMR technology that cannot balance the dynamic range of the diffusion coefficient measurements and the transverse relaxation resolution in practical applications, an improved design ofD-T22D NMR pulse sequence, based on two "windows" and bivariate on pulsed filed gradient, was proposed, after summarizing the technical characteristics of PFG (pulsed filed gradient), STE-PFG (stimulated echo pulsed filed gradient), BP-PFG(Bi-polar pulsed filed gradient), modified CPMG (Car-Purcel-Meiboom-Gill), "diffusion editing", and multi echo intervals CPMG pulse sequences. The numerical simulation results of different structural models (gas-water, oil-water and heavy oil-water) under different signal-to-noise ratios show that the modifiedD-T2pulse sequence balances the diffusion coefficient measurement range and transverse relaxation resolution successfully. The new pulse sequence has wide applications in the exploration and formation evaluation for unconventional reservoir, and will benefit future development ofD-T22D NMR technology.

2D NMR; diffusion coefficient; transverse relaxation resolution; pulsed filed gradient; bivariate; modified pulse sequence; numerical simulation

2014-06-12

“十二五”國(guó)家油氣重大專項(xiàng)(2011ZX05020008-002)；國(guó)家重大科學(xué)儀器設(shè)備開發(fā)專項(xiàng)(2013YQ170463)；國(guó)家自然科學(xué)基金項(xiàng)目(41174009,41474100)

吳飛(1986-),男,博士研究生,研究方向?yàn)閮?chǔ)層巖石物理實(shí)驗(yàn)方法與應(yīng)用。E-mail:feizai123@126.com。

1673-5005(2015)01-0050-10

10.3969/j.issn.1673-5005.2015.01.007

P 631.813

A

吳飛,范宜仁,王帥,等.D-T2二維核磁共振脈沖序列改進(jìn)設(shè)計(jì)及性能對(duì)比[J].中國(guó)石油大學(xué)學(xué)報(bào):自然科學(xué)版,2015,39(1):50-59.

WU Fei, FAN Yiren, WANG Shuai, et al. An improved design ofD-T22D NMR pulse sequence and performance comparison [J]. Journal of China University of Petroleum (Edition of Natural Science), 2015,39(1):50-59.