微波加熱對致密砂巖孔隙水的影響

姚俊輝，陶明，郭陳響

(1.中南大學(xué)資源與安全工程學(xué)院，湖南長沙，410083；2.新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院，新疆烏魯木齊，830001)

致密砂巖氣已經(jīng)成為我國非常規(guī)天然氣的重要組成部分[1-3]。由于致密砂巖儲層發(fā)育微米或納米級孔喉和富含高黏土礦物，含水飽和度較低，毛管力強，因此，侵入的液相極易在孔隙結(jié)構(gòu)中吸附滯留造成水鎖損害[4-6]。致密砂巖氣藏在開發(fā)過程中的水鎖損害率一般為70%～90%[7-9]，液相侵入會導(dǎo)致儲層滲透率大幅度降低，如當(dāng)儲層含水飽和度為80%時，氣相相對滲透率幾乎降為0[10]。水鎖損害會導(dǎo)致氣體產(chǎn)能大幅度下降，甚至發(fā)生不產(chǎn)氣的極端情況。

目前，解除儲層水鎖損害的常見方法是通過化學(xué)手段對吸附滯留水相進行返排，但致密砂巖儲層的物性特征導(dǎo)致化學(xué)溶劑難以注入儲層，因此，化學(xué)手段對解除致密砂巖儲層水鎖的效果并不理想[11]。為了尋找有效解除水鎖的方法，研究者嘗試通過預(yù)熱地層產(chǎn)生的蒸汽壓力來驅(qū)替儲層中吸附滯留的水相。但傳統(tǒng)的地層加熱是通過熱傳導(dǎo)的方式將熱量層層傳遞至水鎖區(qū)域，能量利用率較低，并且當(dāng)其遇到斷層等地質(zhì)阻礙時，加熱范圍非常有限。微波加熱作為一種綠色、高效的熱處理技術(shù)具有體積加熱的特點，因此，可以忽略不良地層地質(zhì)因素對熱量傳導(dǎo)過程帶來的影響。此外，微波會優(yōu)先對多孔巖石介質(zhì)中的水分進行加熱，這會進一步提高微波的能量利用效率。SEEHRA 等[12]發(fā)現(xiàn)微波加熱的干燥效率比傳統(tǒng)的傳熱加熱效率高1 個數(shù)量級，而FU 等[13]則發(fā)現(xiàn)在相同加熱效果下，微波加熱的單位能耗遠比傳熱加熱的能耗低。此外，許多研究已經(jīng)證實微波輻射能夠顯著改變巖石的孔隙結(jié)構(gòu)，增大巖石的孔隙體積，這會進一步為水分的轉(zhuǎn)移和滲透提供更大的活動空間[14-18]。微波加熱具有高效、節(jié)能、易控、安全等優(yōu)點，通過微波預(yù)熱地層解除水鎖損害是一種理想的熱處理方法。目前，關(guān)于微波干燥特性的研究主要采用多模諧振腔進行微波實驗，放置于腔內(nèi)的巖石樣品由于吸收不同方向的反射微波導(dǎo)致巖石樣品被過度加熱，對應(yīng)的實驗結(jié)果難以反映微波注熱技術(shù)的現(xiàn)場效果。在微波預(yù)熱地層的工程現(xiàn)場，微波源前方是無限的巖石介質(zhì)，微波在巖體內(nèi)單向傳播的同時將巖石內(nèi)的水分轉(zhuǎn)變?yōu)樗魵?，?dǎo)致微波在很短的傳播距離內(nèi)被吸收，但隨著儲層含水飽和度降低，微波的傳播范圍又會進一步擴大[19]。因此，研究滿足工程應(yīng)用場景的微波干燥特性對解除儲層的水鎖損害具有重要的實際意義。本研究為模擬微波預(yù)熱地層的工程場景，實驗前對致密砂巖樣品進行預(yù)處理。在微波加熱過程中，利用核磁共振技術(shù)測試巖石樣品的含水量，以揭示微波作用下致密砂巖孔隙水的遷移規(guī)律，研究微波對飽水致密砂巖的干燥效果及微波的干燥機制。

1 核磁共振的原理

核磁共振測試(NMR)利用低頻磁場和短脈沖間隔計算高頻電磁波作用下氫質(zhì)子的弛豫時間，是一種研究材料中水分狀態(tài)和含量的定量測試方法[20]。與常規(guī)測量材料中水分含量的質(zhì)量法相比，核磁共振測試可以消除加熱過程中由于個別礦物質(zhì)揮發(fā)而帶來的實驗誤差，同時探測巖石內(nèi)部的水分分布。目前，核磁共振技術(shù)已經(jīng)在巖土工程領(lǐng)域得到廣泛應(yīng)用，但主要研究集中在巖石孔隙結(jié)構(gòu)的表征[16]、儲層滲透率[15]和巖石損傷[21]的分析方面，而利用核磁共振技術(shù)研究微波作用下巖石中水分遷移的研究較少，有限的研究也只局限于調(diào)查微波輻射時間、微波功率和樣品質(zhì)量對微波干燥能力的影響[22-23]，缺少對微波作用下巖石內(nèi)部水分分布隨時間和空間變化的深入研究。在核磁共振測試過程中，外部磁場的施加使氫磁矩定向排列并產(chǎn)生偶極矩，偶極矩的振幅與氫原子數(shù)成正比，可以用橫向弛豫時間圖譜(T2圖譜)進行表征[17]。根據(jù)核磁共振基本理論[24]，弛豫時間與擴散弛豫時間、體積弛豫時間和表面弛豫時間的關(guān)系為

式中：T2為總弛豫時間，ms；T2B為體積弛豫時間，ms；T2S為表面弛豫時間，ms；T2D為擴散弛豫時間，ms。

在低均勻磁場和短脈沖環(huán)境下，核磁共振的總弛豫時間T2主要為在水巖界面處的表面弛豫時間，而表面弛豫時間為與孔隙表面體積比有關(guān)的函數(shù)。在這種情況下，式(1)可以簡化為

式中：ρ為橫向弛豫強度的常數(shù)；S為孔隙表面積，μm2；V為孔隙體積，μm3。

2 實驗材料和方法

2.1 巖樣和設(shè)備

實驗中采用的紅砂巖樣品來自于四川盆地川西地區(qū)的致密砂巖氣藏，核磁孔隙率為6.99%～7.90%。微波實驗采用由4個功率為1.5 kW的磁控管和1 個封閉金屬腔組成的多源多模微波加熱系統(tǒng)，該系統(tǒng)的微波功率在1～6 kW 范圍內(nèi)連續(xù)可調(diào)。

2.2 樣品預(yù)處理

在微波預(yù)熱地層的工程現(xiàn)場，微波在地層內(nèi)的傳播符合微波在自由空間傳播的特點，即：微波在傳播過程中不斷被巖石中的水分和礦物質(zhì)所吸收，導(dǎo)致微波能量不斷發(fā)生衰減；由于金屬腔對微波的反射作用，放置于多模諧振腔內(nèi)的巖石樣品吸收來自不同方向的反射微波，導(dǎo)致巖石樣品被均勻加熱。因此，研究符合現(xiàn)場實際的微波除水效果，重點是讓微波在巖石樣品內(nèi)部進行單向傳播。

為控制微波在巖石樣品內(nèi)的傳播路徑，實驗中對致密砂巖樣品進行以下預(yù)處理：

1)從完整性良好的砂巖塊中提取巖芯，并將其切割打磨成9個直徑為50 mm、高度為20 mm的圓柱狀樣品，為便于區(qū)分，簡稱此圓柱狀樣品為巖芯塞。巖芯塞無可見裂紋，加工精度為±1 mm。

2)將巖芯塞按順序堆疊為圓柱狀砂巖樣品，并分別用厚度為10 mm 和0.1 mm 的陶瓷纖維棉和銅箔進行包裹。

陶瓷纖維棉的作用是在測溫環(huán)節(jié)中對樣品進行保溫以保證數(shù)據(jù)的準(zhǔn)確性，其本身不具備吸水性，因此，不會在后續(xù)稱重環(huán)節(jié)產(chǎn)生實驗溫差。處理后的砂巖樣品在微波腔體內(nèi)僅受到樣品軸向方向(Z方向)的微波穿透，原本作用于樣品側(cè)面(X和Y方向)的微波被銅箔反射，如圖1所示。此外，由于堆疊后的巖芯塞間距遠小于微波波長(12.2 cm)，因此，微波在樣品中傳播時，在巖芯塞之間的能量損失可以忽略不計。

圖1 實驗中的微波傳播路徑Fig.1 Microwave propagation path during microwave heating

2.3 實驗流程

已有研究表明，高功率的微波加熱更有利于微波干燥，因此，本實驗采取功率為6 kW的微波加熱。在微波加熱過程中，利用蘇州紐邁公司生產(chǎn)的核磁共振分析與成像系統(tǒng)(AniMR-150)探測樣品中的水分。核磁實驗在共振頻率為21 MHz和溫度為32 ℃的環(huán)境下進行。具體的實驗設(shè)備和主要流程如圖2所示。具體流程為：

圖2 實驗設(shè)備與流程Fig.2 Experimental equipment and procedure

1)對紅砂巖樣品在50 ℃的溫度下干燥24 h，隨后立即進行核磁測試和稱質(zhì)量；

2)對干燥紅砂巖樣品真空飽水6 h，隨后立即進行核磁測試和稱質(zhì)量；

3)對飽水砂巖樣品分別加熱30，64，98，132和156 s，隨后立即進行核磁測試和稱質(zhì)量。

需要說明的是，64 s為累計加熱時長，其含義是在上一次加熱30 s后取出樣品完成核磁測試，然后對樣品加熱34 s再次進行核磁測試。后續(xù)加熱至98，132和156 s分別進行核磁測試。

3 結(jié)果分析

3.1 樣品預(yù)處理法的實驗驗證

微波作為電磁波譜的一部分，在通過自由空間中的介電材料時發(fā)生衰減。微波從材料表面衰減到初始能量1/e的深度被定義為穿透深度，它與電磁波的頻率和材料的介電常數(shù)有關(guān)。對于巖石介質(zhì)，損耗因子小于介電常數(shù)，微波對巖石的穿透深度可以用下式計算[25]：

式中：z為穿透深度(m)；λ0為微波波長(m)；ε′為材料的介電常數(shù)；ε″為材料的損耗因子。

巖石材料吸收微波后導(dǎo)致溫度升高，根據(jù)熱力學(xué)定律，溫度升高到一定程度所需的能量由下式計算[26]：

式中：P為所需能量，J；ρ為材料的密度，kg/m3；cp為材料的比熱容，J/(kg·K)；ΔT為材料的溫度差，K；Δt為加熱的時間差，s。

當(dāng)巖石暴露于開放波導(dǎo)發(fā)射的微波時，微波穿透深度內(nèi)巖石材料中的能量密度呈指數(shù)下降[27]：

式中：P(z)為深度z處的能量密度，w；P0為入射能量密度，w；z0為微波的穿透深度，m。

微波的穿透深度與功率無關(guān)(式(3))，微波場中的巖石材料由于吸收微波能而溫度升高，并隨微波穿透深度增加呈指數(shù)衰減(式(4)和(5))。為驗證樣品預(yù)處理方法的有效性，對干燥狀態(tài)下20 個直徑為50 mm、高度為20 mm 的巖芯塞按照2.2 節(jié)中步驟進行預(yù)處理，隨后分別在1～6 kW的微波功率下對處理后的巖石樣品加熱120 s，實驗結(jié)束后，立即利用手持式測溫槍(Hima AT800)在90 s內(nèi)測量巖芯塞表面溫度。測溫實驗重復(fù)3次，最終溫度取平均值，得到的測溫曲線如圖3所示。

圖3 樣品內(nèi)部的溫度梯度Fig.3 Temperature gradient inside sample

從圖3可見：巖芯塞的溫度隨微波在巖石樣品中的傳播呈指數(shù)下降；在不同微波功率作用下，實測的微波穿透深度保持恒定，大約為16 cm(圖3中的黑色實線)。實驗結(jié)果和理論分析結(jié)果一致，表明實驗中對樣品預(yù)處理法獲得的溫度分布符合微波在自由空間傳播的特點。值得注意的是，微波在樣品底部的反射產(chǎn)生了駐波效應(yīng)，提高了加熱效果(圖3中的黑色虛線區(qū)域)，這與文獻[28]中的數(shù)值模擬結(jié)果一致。

3.2 水分遷移規(guī)律和機制

對飽水巖芯塞微波加熱的實驗結(jié)果表明，位于砂巖樣品底部的4個巖芯塞溫度偏高，因此，本文只對位于樣品上端的5 個巖芯塞(2，4，6，8 和10 cm處)進行研究(見圖4)。

圖4 樣品在腔體中的位置Fig.4 Location of sample in cavity

圖5所示為樣品中不同位置巖芯塞在不同干燥時間下的T2圖譜。為了研究方便，本次研究定義弛豫時間小于10 ms的孔隙為較小孔隙，弛豫時間大于10 ms 的孔隙為較大孔隙。從圖5可見：實驗采用的5 塊巖芯塞的T2圖譜呈現(xiàn)典型的雙峰型分布，說明砂巖樣品具有相似的孔隙分布特征；圖譜中右側(cè)波峰高于左側(cè)波峰，說明實驗樣品的孔隙系統(tǒng)以較大孔隙為主；波峰之間的波谷稍低于左側(cè)波峰，表明實驗樣品的孔隙系統(tǒng)具有良好的連通性；2 cm 處的巖芯塞在微波干燥30 s 后，較大孔隙對應(yīng)的孔隙信號明顯衰減，且T2曲線的右側(cè)分布范圍收縮，說明較大孔隙中的水分在微波作用下由液相轉(zhuǎn)變?yōu)闅庀?，并在蒸汽壓力作用下首先從砂巖中擴散，失去水分的孔隙又為后續(xù)水蒸氣的遷移和滲透提供了活動空間；微波干燥64 s后，較大孔隙中的水分含量不斷減少，較小孔隙中的水分也開始蒸發(fā)；隨著微波干燥繼續(xù)進行，砂巖樣品內(nèi)賦存的水分不斷被去除，但水分含量的變化程度不斷減?。划?dāng)微波干燥進行至132 s時，砂巖內(nèi)的水分含量基本保持恒定，說明6 kW微波所攜帶的能量無法將砂巖樣品內(nèi)剩余的水分完全去除，此時，微波的干燥能力已達到極限。值得注意的是，T2圖譜的曲線形狀在微波干燥過程中基本保持不變，基本呈現(xiàn)雙峰型的分布特征，即使在完全干燥的狀態(tài)下依然存在大孔隙中水分對應(yīng)的孔隙信號，導(dǎo)致這種現(xiàn)象的原因可能是：致密砂巖孔隙系統(tǒng)較復(fù)雜，其中較大的盲孔或者閉孔鎖住了部分水分，現(xiàn)有6 kW微波的能量水平較低，無法使被困的水分突破孔隙壁的封鎖而滯留于孔隙內(nèi)部；此外，在核磁共振測試過程中，大孔隙粗糙表面包含的微小孔隙體積也常被作為大孔隙體積的一部分，導(dǎo)致實驗結(jié)果存在誤差[29]。

與水相比，砂巖樣品中的石英和長石是較弱的微波吸收介質(zhì)，微波實驗中的微波能主要用于砂巖內(nèi)水相的變化。在微波作用下，較大孔隙中的水分率先突破孔隙結(jié)構(gòu)的阻礙而散失，而較小孔隙中的水分依然滯留于砂巖內(nèi)部。致密砂巖中的較大孔隙主要為分布于長石和石英顆粒周圍的粒間孔隙[30-32]以及由長石顆粒或者巖石碎片內(nèi)部溶蝕產(chǎn)生的粒內(nèi)孔隙[33-34]；較小孔隙則賦存于黏土內(nèi)部[35-36]。由于黏土的親水性和小孔隙吼道較強的毛管力，導(dǎo)致砂巖微孔里的水分在微波的刺激下較難脫離孔隙結(jié)構(gòu)的束縛而繼續(xù)滯留于砂巖內(nèi)部，因此，微波對大孔隙中的水分去除效果更好。當(dāng)微波實驗進行到98 s時，砂巖樣品中的較大孔隙基本處于干燥狀態(tài)。隨著砂巖中的孔隙水轉(zhuǎn)變?yōu)樗魵鈴纳皫r樣品中不斷逸出，微波的干燥范圍進一步擴大。圖5表明：隨著微波干燥時間增加，巖芯塞中水分含量不斷減少，但減少的幅度隨微波作用范圍擴大而不斷降低。

圖5 樣品不同位置的T2圖譜Fig.5 T2 spectrum of different locations in sample

利用核磁共振技術(shù)測定巖石中水分的含量，必須揭示T2圖譜積分面積與含水量之間的關(guān)系[37]。實驗中，通過高精度電子秤(精度為0.01)獲得加熱后巖芯塞中的含水量(見表1)。本文以2 cm處的巖芯塞為例，求解獲得T2圖譜積分面積與含水量之間的關(guān)系：

表1 2 cm處巖芯塞不同加熱時間后的含水量Table 1 Water content of core plug at 2 cm after different heating time

式中：S為圖譜積分面積；A為Ti處的孔隙信號振幅；Ti為水溶液的橫向弛豫時間，ms。

圖6所示為樣品2 cm 處的巖芯塞在不同加熱時長下的T2圖譜以及圖譜面積和含水量之間的擬合結(jié)果。從圖6可見含水量與圖譜面積呈線性關(guān)系，相關(guān)系數(shù)為0.997 05，表明擬合效果很好，線性擬合方程可以表示為

圖6 2 cm處巖芯塞的T2圖譜與含水量的對應(yīng)關(guān)系Fig.6 Correspondence between T2 spectra and water content of core plug located at 2 cm of sample

式中：M為巖石樣品中的含水量，g；S為T2圖譜的積分面積。

3.3 微波的干燥機制和效果

微波是頻率為300 MHz 至300 GHz 的電磁波，其加熱效果與目標(biāo)材料的介電性能密切相關(guān)，微波加熱相較于常規(guī)加熱具有選擇性加熱的特點。因此，對于巖石這種混合材料，微波會優(yōu)先被巖石內(nèi)具有高介電損耗的物質(zhì)所吸收，同時釋放出熱量。實驗中，砂巖樣品的主要礦物質(zhì)成分為石英、長石和黏土類礦物，均屬于較弱的微波吸收介質(zhì)，微波通過它們時發(fā)生的損耗可以忽略不計?？紫吨械乃肿鳛橐环N高損耗介質(zhì)，吸收微波后水分子瞬間發(fā)生極化，然后以每秒數(shù)億次的速度作來回往復(fù)運動，分子間相互碰撞和摩擦產(chǎn)生的熱能使孔隙水轉(zhuǎn)變?yōu)樗魵狻？紫端鄳B(tài)改變導(dǎo)致巖石內(nèi)外的氣壓不再平衡，巖石內(nèi)部的高壓水蒸氣通過孔隙之間的吼道向較低壓力的梯度方向(巖石外)運移，直至巖石內(nèi)外壓力一致。因此，在微波干燥初期，微波會在很短的傳播距離內(nèi)被巖石中的水分所吸收，但隨著孔隙水變成水蒸氣后不斷外排，微波的干燥范圍進一步擴大。因此，微波的干燥能力可以由微波的傳播距離和干燥時間進行表征。

為研究微波對致密砂巖的干燥效果，利用式(7)和(8)計算微波干燥后巖石樣品的含水率。

式中：ω為巖石含水率，%；Mi為微波干燥后的樣品中的含水量，g；M0為飽水狀態(tài)下的樣品中的含水量，g。

圖7(a)和圖7(b)所示為微波干燥時間對致密砂巖含水率和干燥速率(vt)的影響，而圖7(c)和圖7(d)所示為微波傳播距離對致密砂巖含水率和干燥速率(vs)的影響，其中，干燥速率vt和vs計算式為：

式中：vt為考慮時間因素的干燥速率，單位為10-2%/s；vs為考慮距離因素的干燥速率，單位為10-2%/cm；t為微波的干燥時間，s；s為微波在致密砂巖中的傳播距離，cm。

從圖7(a)和圖7(b)可以看出：微波的干燥速率在加熱前60 s內(nèi)急劇降低，巖石樣品的含水率在微波的傳播范圍內(nèi)也由飽水狀態(tài)的100%降低至50%以下；隨著干燥過程繼續(xù)，干燥速率曲線變緩并逐漸趨于平穩(wěn)(接近0)，因此，在隨后100 s內(nèi)樣品的含水率最多只降低了20%，剩下10%～30%的孔隙水則滯留于致密砂巖的孔隙內(nèi)部，無法被功率為6 kW 的微波去除。圖7(c)和圖7(d)表明加熱初期的100 s 致密砂巖的含水率在10 cm 的微波傳播距離內(nèi)顯著降低，加熱100 s后巖石樣品的底部(6～10 cm)的干燥速率基本為0；繼續(xù)加熱時，含水率也基本保持不變(25%左右)。

圖7 微波對砂巖樣品的干燥效果Fig.7 Drying effect of microwave on sandstone sample

以上結(jié)果表明：巖石中的孔隙水對時長的敏感程度更高，在功率為6 kW的微波作用下，可以在100 s 內(nèi)基本去除致密砂巖內(nèi)部10 cm 之內(nèi)的水分；若需要繼續(xù)去除剩余10%～30%的水分，繼續(xù)增加微波的干燥時間不可行，需要往目標(biāo)砂巖內(nèi)注入更高能量密度的微波。

上述研究結(jié)果表明：在一定能量密度微波作用下，微波的干燥效果主要由微波的干燥時間和傳播距離所決定，換句話說，致密砂巖含水率是干燥時間和傳播距離的函數(shù)。為了確定這種關(guān)系，采用poly 函數(shù)對致密砂巖含水率在時間和空間中的分布進行曲面擬合，結(jié)果如圖8所示。擬合優(yōu)度為0.999 81，表明變量之間相關(guān)性良好，致密砂巖含水率的函數(shù)關(guān)系式為

圖8 曲面擬合結(jié)果Fig.8 Surface fitting results

需要注意的是，本次實驗的研究結(jié)果只針對6 kW微波功率作用下雙峰型孔隙結(jié)構(gòu)的砂巖樣品，而微波功率、巖石內(nèi)部的孔隙結(jié)構(gòu)都會影響微波的干燥效果。因此，全面評價微波技術(shù)對含水致密砂巖儲層的除水效果還需針對特定的目標(biāo)巖體進行具體分析，從而為微波除水技術(shù)的工程應(yīng)用提供指導(dǎo)。

4 結(jié)論

1)由于毛管力的作用和黏土的親水性導(dǎo)致砂巖中的水分難以脫離較小孔隙的束縛。對一定能量密度的微波而言，其干燥能力存在1 個時間閾值，超過時間閾值的微波加熱無法繼續(xù)提高微波的干燥效果，應(yīng)注入更高能量密度的微波以提高微波干燥效果。

2)建立了致密砂巖T2圖譜面積與含水量之間的線性關(guān)系，發(fā)現(xiàn)巖石中的含水率隨干燥時間和傳播距離的增加而減少，而且?guī)r石中的水分對時間的敏感程度更高。功率為6 kW的微波加熱100 s后可以基本去除10 cm 致密砂巖內(nèi)70%以上的水分。通過對干燥時間、傳播距離和含水率的三維空間分布進行擬合，確定了以poly 函數(shù)為基礎(chǔ)的致密砂巖含水率與微波干燥時間和傳播距離的關(guān)系。