曹雪梅 詹永康 袁磊 舒中潘 張軍

摘要：核磁共振技術(shù)（NMR）憑借其便捷性、低成本、無損傷的優(yōu)勢被廣泛應(yīng)用于黃土孔隙空間研究領(lǐng)域中。但目前對于重塑黃土核磁共振T2譜的分析應(yīng)用成果過于單一，從而造成了一定程度上的數(shù)據(jù)浪費(fèi)。有鑒于此，基于核磁共振技術(shù)，利用數(shù)學(xué)分析手段對不同干密度重塑黃土的T2譜開展了定量分析;將分析得到的數(shù)學(xué)參數(shù)應(yīng)用于石油領(lǐng)域中的SDR模型，并構(gòu)建了重塑黃土的滲透率計(jì)算模型。結(jié)果表明：① 不同干密度重塑黃土的核磁共振T2譜在數(shù)學(xué)分析參數(shù)（均值、幾何均值、標(biāo)準(zhǔn)差以及峰度）中存在不同程度的差異，反映到物理意義上，則是重塑黃土干密度會(huì)對土體內(nèi)部孔隙大小、分布及分選特征產(chǎn)生影響。相比傳統(tǒng)的方法僅獲取土體孔隙體積，數(shù)學(xué)統(tǒng)計(jì)方法獲得的數(shù)學(xué)參數(shù)對孔隙的描述更全面也更精確。② 借助于構(gòu)建的重塑黃土的SDR模型，求解了試驗(yàn)黃土的滲透率，計(jì)算的最小滲透率約為0.3 mD，計(jì)算結(jié)果與前人總結(jié)的規(guī)律接近。③ 干密度會(huì)影響多個(gè)土體微觀空間特征，從而引起土體滲透特性的改變。④ 從相關(guān)關(guān)系方面的分析來看，T2譜的均值和幾何均值對土體滲透率影響較大。研究成果在一定程度上可為核磁共振技術(shù)在黃土領(lǐng)域中的應(yīng)用提供新的思路，可為黃土滲透性的預(yù)測提供新的方法，也可為了解土體滲透機(jī)制提供新的思考。

關(guān) 鍵 詞：重塑黃土; 孔隙空間; 核磁共振技術(shù)（NMR）; 微觀特征; 定量參數(shù)

中圖法分類號(hào)： P66

文獻(xiàn)標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2022.05.033

0 引 言

在中國西北地區(qū)，黃土極為重要，分布面積達(dá)63萬km2，主要分布在山西、陜西、甘肅和寧夏地區(qū)[1]。典型黃土是一種具有開放和亞穩(wěn)定結(jié)構(gòu)的風(fēng)沙沉積[2]。一些黃土表現(xiàn)出水敏性，在濕潤[3]時(shí)，由于內(nèi)部結(jié)構(gòu)塌陷趨于不可控沉降，黃土地基塌陷等巖土工程問題與壓實(shí)黃土的水分滲流特性密切相關(guān)，而作為流體在壓實(shí)黃土中的流動(dòng)通道，孔隙空間起著至關(guān)重要的作用[4-6]。因此，研究重塑黃土微觀孔隙特征，是了解工程中黃土宏觀滲透表現(xiàn)的基礎(chǔ)。

核磁共振技術(shù)（NMR）憑借其便捷性、低成本、無損傷的優(yōu)勢，被廣泛應(yīng)用于研究分析巖土體孔隙空間結(jié)構(gòu)的巖土工程領(lǐng)域，其精確性與穩(wěn)定性均得到了研究人員的驗(yàn)證。田慧會(huì)等[7]利用自由水與結(jié)合水冰點(diǎn)不同的特性獲取了T2截止值，通過T2截止值得到了T2譜代表的結(jié)合水含量。Yuan等[8]通過核磁共振T2譜，研究了不同干密度和干濕循環(huán)條件下壓實(shí)黃土的滲透特性及結(jié)構(gòu)演化。還有大量學(xué)者通過核磁共振技術(shù)研究了不同巖土體的孔隙特征，進(jìn)而分析了土體某些宏觀表現(xiàn)行為，獲得了許多科研成果，效果良好[9-11]。

T2譜作為核磁共振技術(shù)的重要成果，可以反映土體內(nèi)孔隙分布、組成等特點(diǎn)，通過對試樣T2譜的特征進(jìn)行定性與定量分析，可進(jìn)一步了解重塑黃土孔隙的空間結(jié)構(gòu)。有學(xué)者將T2譜與壓汞試驗(yàn)相結(jié)合，從而提出了非飽和土的滲透系數(shù)的預(yù)測方法[12];還有學(xué)者對T2譜進(jìn)行了數(shù)學(xué)分析，將數(shù)學(xué)參數(shù)應(yīng)用于石油開采率上，從新的角度了解了砂巖儲(chǔ)油層中孔隙空間特征[11]。

針對黃土的孔隙結(jié)構(gòu)空間，已有部分學(xué)者利用核磁共振技術(shù)對其開展了較多的研究[13-14]，但對于T2譜的應(yīng)用往往僅局限于獲取土體孔隙直徑大小的分布上，并未對T2譜本身所含有的其他空間信息進(jìn)行進(jìn)一步提取，更未進(jìn)行定量分析，成果應(yīng)用過于單一，這就使得核磁共振T2譜在黃土領(lǐng)域中的應(yīng)用不夠充分，而且在一定程度上也造成了數(shù)據(jù)的浪費(fèi)。

本文基于核磁共振技術(shù)獲取了不同干密度重塑黃土的T2譜，并利用數(shù)學(xué)分析手段對其開展了定量分析，并將定量分析得到的數(shù)學(xué)參數(shù)應(yīng)用到求取不同干密度重塑黃土的滲透率計(jì)算模型的構(gòu)建中。該方法在一定程度上可為核磁共振技術(shù)在黃土領(lǐng)域中的應(yīng)用提供一個(gè)新的思路，也可為黃土滲透性的計(jì)算提供新的方法;同時(shí)，也有助于了解土體的滲透機(jī)制。

1 NMR技術(shù)原理

核磁共振（NMR）技術(shù)是基于H原子的自旋磁性和外加磁場之間的相互作用對艙內(nèi)試樣中的水進(jìn)行檢測的。信號(hào)大小與氫原子的個(gè)數(shù)成正比，與橫向磁化矢量的衰減成正比。核磁共振弛豫時(shí)間與試樣的孔隙半徑成正比，這是因?yàn)榇罂字袃?chǔ)水量大，而微孔中儲(chǔ)水量小。核磁共振是指在外部磁場作用下能級(jí)之間的共振躍遷[15]。在這種磁場下，可以得到偏轉(zhuǎn)后的質(zhì)子數(shù)和質(zhì)子自轉(zhuǎn)軸恢復(fù)到初始平衡態(tài)所需要的時(shí)間（T2）[16]。

多孔介質(zhì)的孔隙內(nèi)部結(jié)構(gòu)、孔隙中的分子運(yùn)動(dòng)及其與黃土的反應(yīng)過程是多孔介質(zhì)研究領(lǐng)域的重要課題。對于土壤孔隙中的流體，存在3種類型的弛豫機(jī)制，即表面流體和自由流體的弛豫、分子的自擴(kuò)散。因此，孔隙流體流動(dòng)的弛豫時(shí)間T2可表示如下[17]：

1T2=1T2f+1T2s+1T2d（1）

式中：T2為核磁共振橫向弛豫時(shí)間，ms;T2f為大容器中孔隙流體的弛豫時(shí)間，ms;T2s和T2d分別為表面流體弛豫引起的孔隙流體弛豫時(shí)間和磁場梯度（ms）擴(kuò)散引起的孔隙流體弛豫時(shí)間。

1T2s=ρ2（SV）（2）

式中：ρ2是橫向弛豫系數(shù)，取決于土壤性質(zhì)[17];S/V是孔隙表面積與流體體積的比值。

當(dāng)只考慮飽和流體時(shí)，自由弛豫可以忽略不計(jì)，其值比表面弛豫小得多。同時(shí)，當(dāng)磁場均勻且回波時(shí)間間隔足夠短時(shí)，擴(kuò)散弛豫也可以忽略不計(jì)。因此，流體的松弛主要來源于土壤顆粒的表面松弛，式（2） 可近似為[18]

1T2=ρ2（SV）=Fsρ2R（3）

當(dāng)孔隙為球形時(shí)，T2與孔隙半徑有關(guān)，F(xiàn)s為孔隙形狀因子（當(dāng)孔隙為圓柱形時(shí)Fs=3）。因此，通過式（3）可知核磁共振T2譜線能夠反演出土體內(nèi)部的孔隙空間結(jié)構(gòu)與大小分布特征[13]。

2 試驗(yàn)過程

2.1 試樣準(zhǔn)備

本次試驗(yàn)所用黃土取自陜西省延安市某挖方邊坡，其基本物理性質(zhì)如表1所列。

依據(jù)GB/T 50123-2019《土工試驗(yàn)方法標(biāo)準(zhǔn)》，將土樣在108 ℃下烘干24 h，然后碾細(xì)后過2 mm篩備用。正如表1所列，試驗(yàn)所用黃土最大干密度為1.8 g/cm3，最優(yōu)含水量為16.3%，據(jù)此本文設(shè)置了5種不同干密度（1.80，1.74，1.68，1.62，1.56 g/cm3）的重塑黃土樣，為了減少試驗(yàn)偶然誤差，另設(shè)置1組平行試驗(yàn)。重塑黃土試樣含水量均為最優(yōu)含水量：16.3%。

試樣配制時(shí)，為了控制制樣均勻，首先，將定量散狀篩后黃土與相應(yīng)的蒸餾水混合攪拌均勻，控制含水量為16.3%，用保鮮膜密封靜置24 h，待土中水分分布均勻。然后，根據(jù)目標(biāo)干密度計(jì)算得到單個(gè)試樣所需的濕土質(zhì)量，先后分3層擊實(shí)進(jìn)入環(huán)刀內(nèi)，每層土間刮毛防止分層現(xiàn)象。由于核磁共振中磁場較強(qiáng)，傳統(tǒng)鋼制環(huán)刀對其影響較大，本試驗(yàn)中所用環(huán)刀為特制聚四氟乙烯材質(zhì)環(huán)刀，聚四氟乙烯可有效消除金屬對NMR信號(hào)的干擾。按照上述步驟，依次完成了2組5個(gè)干密度重塑黃土試樣。

2.2 試驗(yàn)過程

采用核磁共振掃描法測量重塑黃土的T2圖譜，所用核磁共振儀器型號(hào)為MacroMR12-150H-I。該核磁共振儀永久磁鐵磁場強(qiáng)度為（0.3±0.05）T，射頻發(fā)射功率的峰值輸出大于300 W，最大線性失真小于0.5%。本次試驗(yàn)所用到的線圈探頭的直徑為60 mm。

在核磁共振中，要得到能完整反映土體空間結(jié)構(gòu)的T2譜線，試樣內(nèi)孔隙必須要先飽水。所以，首先對聚四氟乙烯內(nèi)的土樣進(jìn)行飽和，將試樣兩端用透水石封住，防止土體孔隙在飽和過程中發(fā)生改變，也可使飽和過程更加均勻。然后將試樣置于真空飽和器中，抽取真空2 h后再進(jìn)行真空飽和12 h，飽和后立刻進(jìn)行NMR常溫掃描，以獲取2組不同干密度重塑黃土的T2圖譜。圖1為飽和中的試樣。

3 核磁共振T2譜圖形分布特征

在核磁共振對土體的測試中，首先對試樣施加磁場，讓土體孔隙中的水的H原子產(chǎn)生共振，然后撤銷磁場。由于H原子從共振的狀態(tài)恢復(fù)到激發(fā)前的狀態(tài)時(shí)會(huì)發(fā)出一定的電信號(hào)，孔隙越大，恢復(fù)到激發(fā)前狀態(tài)的橫向弛豫時(shí)間也就越久。同一個(gè)孔徑大小的孔隙越多，體積就越多，產(chǎn)生的電信號(hào)也就越強(qiáng)。利用傳感器收集這些電信號(hào)，然后經(jīng)過數(shù)學(xué)反演就得到了T2譜線，橫坐標(biāo)為弛豫時(shí)間，其物理意義是孔徑大小，弛豫時(shí)間越大，代表的孔隙越大?？v坐標(biāo)是信號(hào)強(qiáng)度，某個(gè)大小的孔隙越多，其信號(hào)強(qiáng)度越強(qiáng)，其物理意義為某個(gè)大小孔隙的體積。

接下來，再將縱坐標(biāo)的信號(hào)強(qiáng)度進(jìn)行歸一化，即用單個(gè)弛豫時(shí)間下的信號(hào)強(qiáng)度去除所有弛豫時(shí)間下的信號(hào)強(qiáng)度的總和，這樣就將信號(hào)強(qiáng)度轉(zhuǎn)化為了孔隙度分量。最大孔隙度分量約為3%。

在核磁共振試驗(yàn)中，不同干密度樣品的橫向弛豫時(shí)間T2分布曲線均為雙峰分布（見圖2），弛豫時(shí)間短則代表了小孔隙，弛豫時(shí)間長則代表了大孔隙。雙峰峰型分布表明重塑黃土體內(nèi)孔隙具有宏觀、中、微觀的特征。

根據(jù)T2譜的特征可分為圖形分布特征與數(shù)學(xué)分布特征，圖形分布特征可用：最大最小弛豫時(shí)間、譜峰弛豫時(shí)間、區(qū)間孔隙分量、最大孔隙度分量和歪度這幾個(gè)特征參數(shù)進(jìn)行描述，下面依次分析其圖形主要特征。

由于2組試樣為平行試驗(yàn)，以下參數(shù)分析均取這2組試驗(yàn)的平均值。

3.1 最大最小弛豫時(shí)間

T2譜中的弛豫時(shí)間越大，其代表的孔隙也就越大，因此，最大最小弛豫時(shí)間代表了重塑黃土內(nèi)部最大與最小的孔隙。由圖2可知：不同干密度的重塑黃土最小弛豫時(shí)間基本一致，為0.01 ms，這是由于儀器本身的限制，所以最小弛豫時(shí)間沒有對比意義。最大弛豫時(shí)間隨干密度的增大而減?。ㄒ妶D3），呈現(xiàn)出非線性的趨勢。

3.2 半弛豫時(shí)間

半弛豫時(shí)間指的是累計(jì)孔隙度達(dá)到50%時(shí)的橫向弛豫時(shí)間，該參數(shù)可在一定程度上反映土體中孔隙比的大小。如表2所列，試樣的半弛豫時(shí)間隨干密度的減小而增大，兩者呈負(fù)相關(guān)的關(guān)系。這也表明，隨著干密度的減小，土體內(nèi)總孔隙體積與孔隙比在不斷的增加。

3.3 譜峰弛豫時(shí)間

譜峰弛豫時(shí)間指的是在核磁共振T2譜中信號(hào)強(qiáng)度最大的弛豫時(shí)間，代表了孔隙大小最為集中的弛豫時(shí)間分量。不同干密度試樣譜峰弛豫時(shí)間基本相差不大，說明重塑黃土在弛豫時(shí)間在0.73 ms左右所代表的孔隙孔徑最為集中，與試樣的密度無明顯關(guān)系（見表3）。

3.4 區(qū)間孔隙分量

核磁共振T2譜線的橫軸弛豫時(shí)間代表了土體中不同大小的孔隙，孔隙越大，其弛豫時(shí)間越大，而譜線縱軸的信號(hào)強(qiáng)度則可以反映該孔隙的體積大小。將不同大小的孔隙分成大中小3類，3類孔隙的體積分別在全孔隙體積中所占的比例被稱為區(qū)間孔隙分量。

根據(jù)譜線的形態(tài)將孔隙分為3類（見圖4）：最小弛豫時(shí)間至1.80 g/cm3試樣半弛豫時(shí)間（0.79 ms），是小孔隙S1;1.80 g/cm3試樣半弛豫時(shí)間（0.79 ms）至次譜峰弛豫時(shí)間（11.1 ms），為中孔隙S2;次譜峰弛豫時(shí)（11.1 ms）至最大弛豫時(shí)間（1 000 ms）為大孔隙S3。

T2譜與橫軸圍成的面積（無量綱）可代表不同大小孔隙的相對體積，所以通過統(tǒng)計(jì)3個(gè)區(qū)間（S1、S2、S3）的譜峰面積，可以計(jì)算得到不同大小的孔隙的相對體積，再分別計(jì)算出3類孔隙體積所占比例，計(jì)算結(jié)果如圖5所示。

圖5中，S2所代表的中孔隙比例隨干密度的變化影響不大，而S1代表的小孔隙，S3代表的大孔隙隨干密度的增大分別呈正線性相關(guān)和負(fù)線性相關(guān)，這表明干密度對于土體內(nèi)部孔隙大小分布影響巨大，干密度越大的土體孔隙越少，孔隙也越集中。

4 核磁共振T2譜數(shù)學(xué)分布特征

在傳統(tǒng)地質(zhì)數(shù)學(xué)統(tǒng)計(jì)模型中，最常用的概率統(tǒng)計(jì)模型是正態(tài)分布統(tǒng)計(jì)。前人應(yīng)用正態(tài)概率分布里的數(shù)學(xué)特征參數(shù)比如均值、方差、中值、變異系數(shù)等，去描述各類巖石與土體的孔隙大小分布情況;而T2譜作為土體內(nèi)部孔隙大小及分布的直接物理體現(xiàn)。因此，T2曲線也能使用數(shù)學(xué)統(tǒng)計(jì)參數(shù)來對其進(jìn)行數(shù)學(xué)分布特征的定量分析。

4.1 均 值

均值是數(shù)學(xué)統(tǒng)計(jì)中非常常用的數(shù)學(xué)位置特征參數(shù)，可由核磁共振T2譜與孔隙度加權(quán)平均求得[11]：

T2=Ni=1T2ii（4）

式中：i為第i個(gè)弛豫時(shí)間對應(yīng)的核磁共振孔隙度分量;

T2i為第i個(gè)弛豫時(shí)間，ms。

由式（4）可求得各重塑黃土試樣的T2譜的均值，能反映全孔隙中的平均孔隙大小。由表4可以看出：試樣T2譜均值隨干密度的減小而增大，呈負(fù)線性相關(guān)關(guān)系。

4.2 幾何均值

幾何均值和均值都是求一組數(shù)據(jù)中平均值的方法[19]，而幾何均值則更多用于數(shù)據(jù)的平均增長率的研究中。幾何均值可由下式定義[20]：

T2gm=（∏Ni=1T2i）i（5）

式中：T2gm為T2的幾何均值;

i為第i個(gè)弛豫時(shí)間對應(yīng)的核磁共振孔隙度分量;

T2i為第i個(gè)弛豫時(shí)間，ms。

按式（5）求得每組試樣的幾何均值。各試樣T2譜幾何均值如表5所列。

幾何均值隨干密度的增大而減小，呈負(fù)相關(guān)關(guān)系，表明重塑黃土的干密度越大，T2曲線的平均增長率越小。

4.3 標(biāo)準(zhǔn)差

標(biāo)準(zhǔn)差是描述離均差平方的算術(shù)平均數(shù)。在孔隙系統(tǒng)中，可用于描述孔隙體積大小的分選性，標(biāo)準(zhǔn)差越小，說明孔隙大小的分選性越好，孔隙大小越集中。可由下式計(jì)算：

σ=（iT2i-T2）2（6）

式中：σ為標(biāo)準(zhǔn)差;

i為第i個(gè)弛豫時(shí)間對應(yīng)的核磁共振孔隙度分量;

T2i為第i個(gè)弛豫時(shí)間，ms;

T2為均值;

為均值對應(yīng)的孔隙度分量。

按式（6）求得每組試樣的標(biāo)準(zhǔn)差。各試樣T2譜標(biāo)準(zhǔn)差如表6所列。

標(biāo)準(zhǔn)差隨干密度的增大而減小，呈負(fù)相關(guān)關(guān)系，表明干密度越大的重塑黃土，其內(nèi)部孔隙分選性越好，孔隙大小越集中在某個(gè)孔徑區(qū)間。

4.4 變異系數(shù)

標(biāo)準(zhǔn)差與均值的比值被稱為變異系數(shù)，是對數(shù)據(jù)值相對變化的體現(xiàn)。變異系數(shù)越小，孔隙結(jié)構(gòu)、分選性也就越好。其計(jì)算式為

C=σ/T2（7）

按式（7）求得每組試樣的變異系數(shù)。各試樣T2譜變異系數(shù)如表7所列。

從表7可以得知：不同干密度試樣核磁共振T2譜的變異系數(shù)基本一致，說明干密度的大小并不會(huì)影響重塑黃土內(nèi)部變異系數(shù)。

4.5 峰 度

峰度在核磁共振中是用來衡量T2譜曲線的陡峭程度，可以反映重塑黃土內(nèi)部不同大小孔隙的集中程度，是用來評(píng)價(jià)孔隙結(jié)構(gòu)很好的一個(gè)定量參數(shù)，計(jì)算公式如下：

KG=（2iT2i-T2）4σ4（8）

按式（8）求得每組試樣的峰度。各試樣T2譜峰度如表8所列。

由表8可以看出：峰度大致隨干密度增大而增大，表明干密度越大，T2曲線越陡，孔隙越集中。

5 定量參數(shù)分析與應(yīng)用

5.1 定量參數(shù)分析

前文對不同干密度重塑黃土的核磁共振T2譜進(jìn)行了定性與定量研究，得到了各類豐富的定量參數(shù)與圖形信息，各類參數(shù)在T2譜中所表達(dá)的物理意義如圖6所示。這些特征參數(shù)對于描述孔隙空間也更為全面，充分挖掘了T2譜中蘊(yùn)含的物理及數(shù)學(xué)信息。通過對每組試樣的核磁共振T2譜的分解，研究了不同干密度所導(dǎo)致重塑黃土孔隙結(jié)構(gòu)空間的異同，更進(jìn)一步地了解了重塑黃土微觀組構(gòu)。

通過對T2譜的圖形分析與定量表征，核磁共振特征參數(shù)可根據(jù)其代表的物理意義的不同而劃分為以下3類。

（1） 均值、幾何均值、半弛豫時(shí)間、譜峰弛豫時(shí)間以及最大弛豫時(shí)間，這類參數(shù)可表征重塑黃土的孔隙體積的大小;

（2） 大中小孔隙的區(qū)間孔隙分量，可代表每種孔隙以及孔喉的分布范圍與占比，與流體運(yùn)動(dòng)直接相關(guān);

（3） 變異系數(shù)、標(biāo)準(zhǔn)差，則表達(dá)了重塑黃土孔隙分布的分選特征與分布規(guī)律，間接影響了水在重塑黃土中的滲流過程。

圖7展示了重塑黃土的數(shù)學(xué)特征參數(shù)與圖像特征參數(shù)隨干密度的變化而變化的分布特征，這些參數(shù)直觀地描述了試樣內(nèi)部孔隙空間特征的具體差異，并將這種差異大小進(jìn)行了定量，有利于對孔隙空間特征進(jìn)行下一步的分析與應(yīng)用。由圖7可以發(fā)現(xiàn)：在圖形特征參數(shù)中，譜峰弛豫時(shí)間沒有太大差異，而半弛豫時(shí)間卻隨干密度的增大而明顯減小，這表明干密度越大，孔隙更加集中，大孔隙也越少，這種差異在區(qū)間孔隙分量這一參數(shù)中得以體現(xiàn)。在數(shù)學(xué)統(tǒng)計(jì)分析參數(shù)中，均值、幾何均值、峰度以及標(biāo)準(zhǔn)差均有顯著變化，而這些差異可以體現(xiàn)在下文對重塑黃土滲透率的影響中。

5.2 定量參數(shù)在滲透率模型求取中的應(yīng)用

作為流體在土體中的流動(dòng)通道，孔隙的大小、形態(tài)、分布等空間特征對于土體滲透性的影響巨大。前文獲取到了可定量表征重塑黃土內(nèi)孔隙空間特征差異的數(shù)學(xué)參數(shù)，這節(jié)將這些參數(shù)應(yīng)用到土體滲透率計(jì)算模型的構(gòu)建中，可為理解土體滲透機(jī)制提供新的思路。

目前，在黃土領(lǐng)域利用核磁共振構(gòu)建滲透率模型的文獻(xiàn)偏少或幾乎沒有，而在石油領(lǐng)域，利用核磁共振構(gòu)建滲透率模型，有多種構(gòu)建方法，Coates經(jīng)典模型的建立需要求得土體中的束縛水體積，而另一種由Klenberg[21]提出的SDR模型則需要求取T2譜的幾何均值。綜合前文獲取的數(shù)學(xué)參數(shù)，本研究在SDR模型基礎(chǔ)上，引入在石油儲(chǔ)層滲透率模型中應(yīng)用較好的SDR模型[11]來對重塑黃土的滲透率計(jì)算模型進(jìn)行研究，該計(jì)算模型公式如下：

K=0.01T2gmσT2hi，maxKG（9）

式中：K為滲透率，mD;

為孔隙度，%;

σ為標(biāo)準(zhǔn)差;

i，max為最大孔隙度分量，%;

T2h為譜峰弛豫時(shí)間，ms;

T2gm為幾何均值，ms;

Kg為峰度。

式（9）中，每組試樣孔隙度可根據(jù)土的三相公式換算得到，其他參數(shù)在前文均已計(jì)算得到，對不同干密度重塑黃土的滲透率模型進(jìn)行構(gòu)建，根據(jù)計(jì)算公式得到每組試樣的滲透率，計(jì)算結(jié)果如表9所列。

重塑黃土的干密度與滲透率呈現(xiàn)良好的負(fù)指數(shù)相關(guān)關(guān)系（見圖6），擬合相關(guān)性好，與前人對重塑黃土滲透性的研究[22-23]基本一致，一般認(rèn)為干密度是對重塑黃土滲透性最大的因素，因?yàn)楦擅芏葧?huì)在很大程度上改變土體內(nèi)部孔隙空間結(jié)構(gòu)，導(dǎo)致流體在流動(dòng)通道內(nèi)的運(yùn)動(dòng)過程發(fā)生改變，從而引起滲透特性的變化。

6 結(jié) 論

基于核磁共振技術(shù)，獲取了不同干密度重塑黃土的T2譜;并利用數(shù)學(xué)分析手段對其進(jìn)行定量分析，將定量分析得到的數(shù)學(xué)參數(shù)應(yīng)用于構(gòu)建重塑黃土的滲透率計(jì)算模型中，探索了一套T2譜在黃土孔隙空間分析中應(yīng)用的新思路?？梢缘贸鲆韵轮饕Y(jié)論。

（1） 利用數(shù)學(xué)統(tǒng)計(jì)手段對核磁共振T2譜進(jìn)行定量分析，可以更直觀與更精確地分析不同試樣間的孔隙空間差異。不同干密度重塑黃土的核磁共振T2譜在數(shù)學(xué)統(tǒng)計(jì)參數(shù)中存在不同程度的差異，這些差異反映到物理意義上則是：重塑黃土干密度不同會(huì)對土體內(nèi)部孔隙大小、分布及分選特征產(chǎn)生影響，這些影響的大小則均可通過參數(shù)的數(shù)值差異進(jìn)行更為準(zhǔn)確的定量表達(dá)。相比傳統(tǒng)僅獲取土體孔隙體積，數(shù)學(xué)統(tǒng)計(jì)方法獲得的數(shù)學(xué)參數(shù)對孔隙的描述更全面與精確。

（2） 將石油領(lǐng)域中SDR模型引入到了重塑黃土的滲透率模型的構(gòu)建中，得到了重塑黃土的滲透率模型;并根據(jù)模型求取了不同干密度重塑黃土的滲透率。根據(jù)模型計(jì)算得到試樣黃土的滲透率為0.30～3.81 mD。其次，干密度會(huì)導(dǎo)致土體滲透特性的改變：干密度越大，滲透率越低，且呈負(fù)指數(shù)的相關(guān)關(guān)系。另外，模型還反映了不同微觀孔隙特征參數(shù)對土體滲透率的貢獻(xiàn)大小，進(jìn)一步理解了干密度變化導(dǎo)致土體滲透率變化的微觀因素。模型計(jì)算結(jié)果與前人獲得規(guī)律相似度高，這表明本文研究方法可信度較高，可行性較好。

本文研究成果在一定程度上可為核磁共振在黃土領(lǐng)域中的應(yīng)用提供新的思路，可為黃土滲透性的預(yù)測提供新的方法，更可為了解土體滲透機(jī)制提供新的思考，具有較高的理論創(chuàng)新性與應(yīng)用實(shí)用性。

參考文獻(xiàn)：

[1] 王數(shù).地質(zhì)學(xué)與地貌學(xué)[M].北京：中國農(nóng)業(yè)大學(xué)出版社，2013.

[2] ROLLINS K M，ROGERS G M.Mitigation measures for small structures on collapsible alluvial soils[J].Journal of Geotechnical Engineering，1994，120（9）：1533-1553.

[3] ASSALLAY A M，ROGERS C D F，SMALLEY I J.Formation and collapse of metastable particle packings and open structures in loess deposits[J].Engineering Geology，1997，48（1-2）：101-115.

[4] SIRVITHAYAPAKORN S，KELLER A.Transport of colloids in saturated porous media：A pore-scaleobservation of the size exclusion effect and colloid acceleration[J].Water Resources Resarch，2003，39（4）：1109-1120.

[5] 李鵬舉，谷宇峰.核磁共振T2譜轉(zhuǎn)換偽毛管壓力曲線的矩陣方法[J].天然氣地球科學(xué)，2015，26（4）：700-705.

[6] 劉勇健，李彰明，郭凌峰，等.基于核磁共振技術(shù)的軟土三軸剪切微觀孔隙特征研究[J].巖石力學(xué)與工程學(xué)報(bào)，2018，37（8）：1924-1932.

[7] 田慧會(huì)，韋昌富.基于核磁共振技術(shù)的土體吸附水含量測試與分析[J].中國科學(xué)：技術(shù)科學(xué)，2014，44（3）：295-305.

[8] YUAN K Z，LU X F，WANG X J.Permeability characteristics and structural evolution of compacted loess under different dry densities and wetting-drying cycles[J].PLoS ONE，2021，16（6）：e0253508.

[9] 運(yùn)華云，趙文杰，劉兵開，等.利用T2分布進(jìn)行巖石孔隙結(jié)構(gòu)研究[J].測井技術(shù)，2002，26（1）：18-21，89.

[10] 李海波，朱巨義，郭和坤.核磁共振T2譜換算孔隙半徑分布方法研究[J].波譜學(xué)雜志，2008，25（2）：273-280.

[11] 白松濤，程道解，萬金彬，等.砂巖巖石核磁共振 T2譜定量表征[J].石油學(xué)報(bào)，2016，37（3）：282-391.

[12] 陳銀.基于核磁共振技術(shù)的非飽和土滲透系數(shù)預(yù)測方法[D].武漢：湖北工業(yè)大學(xué)，2019.

[13] 何攀，許強(qiáng)，劉佳良，等.基于核磁共振技術(shù)的結(jié)合水含量對重塑黃土抗剪強(qiáng)度影響試驗(yàn)研究[J].山地學(xué)報(bào)，2020，38（4）：571-580.

[14] 楊更社尤，吳迪.凍融環(huán)境下原狀黃土孔徑分布與其力學(xué)特性關(guān)系的試驗(yàn)研究[J].煤炭工程，2019，51（3）：107-112.

[15] SCHAUMANN G E.Proton nuclear magnetic resonance（NMR）relaxometry in soil science[J].Encyclopedia of Agrophysics，2011：543-669.

[16] CONTE P，STEFANO C D，F(xiàn)ERRO V，et al.Assessing hydrological connectivityinside a soil by fast-fieldcycling nuclear magnetic resonance relaxometry and its link to sediment delivery processes[J].Environ.Earth Sci.，2017，76：526.

[17] XU J，LI Y，REN C，et al.Damage of saline intact loess after dry-wet and its interpretation based on SEM and NMR[J].Soils and Foundations，2020，60（4）：911-928.

[18] 何雨丹，毛志強(qiáng)，肖立志.核磁共振T2分布評(píng)價(jià)巖石孔徑分布的改進(jìn)方法[J].地球物理學(xué)報(bào)，2005，48（2）：373-378.

[19] 何更生，唐海.油層物理[M].2 版.北京：石油工業(yè)出版社，2011.

[20] 楊正明，李治碩，王學(xué)武，等.特低滲透油田相對滲透率曲線測試新方法[J].石油學(xué)報(bào)，2010，31（4）：629-632.

[21] KLEINBERG R L.Utility of NMR T2 distributions，connection with capillary pressure，clay effect，and determination of the surface relaxivity parameter ρ2[J].Magnetic Resonance Imaging，1996，14（7-8）：761-767.

[22] 高燕燕.重塑黃土的滲透特性研究[D].西安：長安大學(xué)，2016.

[23] 高文靜.細(xì)粒土滲透特性的影響因素試驗(yàn)研究[D].西安：西安理工大學(xué)，2018.

（編輯：趙秋云）

Quantitative analysis on NMR T2 spectrum in remolded loess and application

CAO Xuemei1，ZHAN Yongkang2，YUAN Lei3，SHU Zhongpan4，ZHANG Jun4

（1.Sichuan Vocational and Technical College of Communications，Chengdu 611130，China; 2.Sichuan Aba Traffic Construction Investment Co.，Ltd.，Aba 624000，China; 3.Transportation Bureau of Aba Tibetan and Qiang Autonomous Prefecture，Aba 624000 ，China; 4.Sichuan Institute of Geological Engineering Investigation Group Co.，Ltd.，Chengdu 610072 ，China）

Abstract：

Nuclear Magnetic Resonance （NMR） is widely used in the analysis of pore space characteristics of loess due to its convenience，low cost，and no damage.However，the previous analysis and application of remolded loess T2 spectrum were limited to obtaining the distribution of pore size，and other spatial information contained in the T2 spectrum was not further extracted，which made the application results of the NMR T2 spectrum too single，resulting in a certain degree of data waste.Based on this，we used mathematical analysis method to quantitatively analyze the T2 spectrum of remolded loess with different dry densities，and applied the mathematical parameters obtained by quantitative analysis to the SDR model，so as to construct a permeability calculation model for remodeled loess.The main conclusions were as follows：（1） There are various differences in the statistical analysis parameters （mean，geometric mean，standard deviation and kurtosis） of the NMR T2 spectrum of loess remolded with different dry densities，which has physical meanings：different dry densities of remolded loess will affect the size，distribution and sorting characteristics of pores in the soil.Compared with the traditional method to obtain the pore volume of soil，the mathematical parameters obtained by numerical statistics method can describe the empty pores more comprehensively and accurately.（2） The permeability of each remolded loess was calculated by constructing the SDR model of the remolded loess，the minimum permeability was about 0.3 mD，and the calculated results were close to the previous experimental results.（3） Dry density affected micro-spatial characteristics of soil，and causing the change of soil permeability characteristics.（4） From the correlation analysis，the mean and geometric mean of T2 spectrum had great influence on soil permeability.To a certain extent，the research results can provide new ideas for the application of NMR in the loess study field，new methods for predicting the permeability of loess，and new thinking for understanding the mechanism of soil permeability.

Key words：

remolded loess;pore space;Nuclear Magnetic Resonance （NMR）;microscopic characteristics;quantitative parameters