霍錦華，張舒航，秦智鵬，謝佳玲，李會芬，門佳坤，張世偉，張霄朋，黃為良

（廣西大學(xué) 機(jī)械工程學(xué)院，廣西 南寧 530000）

0 引言

隨著防砂完井技術(shù)、復(fù)合射孔優(yōu)化新技術(shù)、二次儲層完井和完井保護(hù)新技術(shù)以及智能完井技術(shù)的發(fā)展[1]，深?？睖y已經(jīng)成為人類找尋新化石能源資源的重要手段。從2012年-2018年，全球約50%的新發(fā)現(xiàn)的油田，約60%的新發(fā)現(xiàn)油氣儲量來自深海水域[2]。因此，進(jìn)一步完善改進(jìn)海洋油氣勘測技術(shù)對于緩解能源短缺至關(guān)重要。目前我國急需探索海底資源，但缺乏對于海床亞層油氣遷徙機(jī)理的了解，主要癥結(jié)在于不了解海床中碳?xì)浠狭黧w的遷徙機(jī)理。

由于海底亞層結(jié)構(gòu)的復(fù)雜性，并且油、氣、固多相流動的流動特性遠(yuǎn)比單相流動復(fù)雜，多相流相間存在著界面效應(yīng)和相對速度，所以導(dǎo)致多相流動過程是一個復(fù)雜多變量的非線性動態(tài)過程，其機(jī)理性研究目前還不夠充分?？v觀國內(nèi)外相關(guān)多相流多孔介質(zhì)方面的研究，主要的研究方向還以管道氣泡流和彈射流為主，無法模擬海底三維平面的復(fù)雜性。1979年，OLE HASSAGER研究了黏性牛頓液體中氣體上升體積和形狀的變化狀態(tài)[3]；1986年TABELING等[4]在前人基礎(chǔ)上研究了在極小間隙下氣體指進(jìn)狀態(tài)的變化，得出指狀的界面不穩(wěn)定與指進(jìn)速率有關(guān)；2015年，Germán Varas等[5]通過對Hele-Shaw氣體聯(lián)通后氣柱邊緣呈拋物線的流化床狀態(tài)，得出流化床的寬度隨著速率增加而變大，并且有分散逐漸集中。2019年周彥霞等[6]針對致密油藏用實驗的方法模擬了鹽水的滲析過程；在現(xiàn)有研究中無法直觀看出影響指進(jìn)現(xiàn)象的因素。

為此，在文獻(xiàn)[7-9]的研究基礎(chǔ)上，設(shè)計制作一套能清楚模擬展現(xiàn)海底油氣遷徙過程的實驗裝置。該裝置依靠Hele-Shaw流動模型建立，可調(diào)配不同成分的懸濁液以模擬不同海床的沉積巖成分[10]，通過注射泵注射氣體，用氣體的逸出過程模擬油氣的遷徙過程，并可以改變氣體流速以獲得不同情況下的實驗結(jié)果。實驗可以直觀展示油氣在不同海域海床亞層結(jié)構(gòu)中滲透及遷徙的過程。這對于我國海底油氣資源勘測及開采具有重大意義，同時也對于評估不同海域油氣開發(fā)的成本與風(fēng)險有借鑒意義，以及對國家海底資源“碳封存”有參考意義。

1 實驗裝置設(shè)計

1.1 實驗原理

本實驗設(shè)計的目的是通過實驗定量分析海底油氣遷徙的動力學(xué)過程，了解和預(yù)測滲透位移和行走趨勢，非固結(jié)介質(zhì)中氣泡的形成和上升的力學(xué)機(jī)制。即多相流中氣體在固、液混合懸浮液中的流動過程。為達(dá)到要求，實驗即實驗裝置需滿足3種需求：

（1）模擬海底油氣遷移的基本環(huán)境與條件。由于海床亞層表面是一塊平直，且厚度有幾百米之高的巨型板塊，所以海床板塊間氣體滲透過程，完全可以簡化成一種兩無限大平板間的多相流過程。

（2）模擬呈現(xiàn)油氣在遷徙過程中以孔隙延展分型的形式擴(kuò)大的過程。油氣在海底以油氣藏的形式儲存，隨著油氣藏的膨脹，頂部油氣壓力增大，最終突破土層的束縛在沉積巖中滲透，油氣在遷徙過程中以孔隙延展分型的形式擴(kuò)大，該實驗需要模擬此延展分型的過程。

（3）模擬兩相置換過程，即為氣相在平流流動的液相中不斷擴(kuò)張、置換和演替推進(jìn)的過程。

因為海床亞層的氣體滲透過程極其緩慢和復(fù)雜，無法使用近似模擬的方式與真實海床對比，流體力學(xué)中可以根據(jù)雷諾數(shù)相等的原理將無法模擬的海床亞層按比例縮小，即構(gòu)建縮小裝置來模擬。關(guān)于海床油氣遷徙的推演模擬實驗，目前在油氣遷徙過程模擬領(lǐng)域有兩種主流模型，一是Hele-Shaw cell模型，二是3D cell模型。3D cell是一種三維的模型，通常以樹脂玻璃圓柱體為容器，填充懸濁液后從底部進(jìn)氣，從而形成三維的孔隙結(jié)構(gòu)，但缺點在于難以觀測到其內(nèi)部的孔隙分布情況。因此基于Hele-Shaw流動模型設(shè)計實驗裝置以及實驗操作步驟。Hele-Shaw流動模型由英國機(jī)械工程師Henry Selby Hele-Shaw在1898年提出，是兩個緊密間隔的平行板之間，某相流體在其他流體間的緩慢流動模型，是黏性流動stokes流的一種特例[11-13]。其原理為：假設(shè)兩平板面積為無限大，且流體為流速極小的充分發(fā)展的定常層流、黏性不可壓縮流體，則在該系統(tǒng)中，流速僅為垂直于板方向位置的函數(shù)，壓力僅為平行平板方向位置的函數(shù)。研究多相流之間交融現(xiàn)象，由于流體間黏性不同，導(dǎo)致流體形成分叉多行的指狀形狀。此外，在Hele-Shaw模型中，進(jìn)入兩板間流體的速度滿足拋物線分布，流體的流動有旋[14]。但是由于流體垂直于平板方向的渦量分量為零，所以從平板上方觀察的時候流動是無旋的。粘度較低的相將比粘度較高的相移動或置換得更快，導(dǎo)致界面不穩(wěn)定，從而導(dǎo)致出現(xiàn)手指狀圖案[15]。

實驗擬通過兩塊透明板形成夾層，將海底土質(zhì)與油氣相互作用簡化為多孔介質(zhì)與空氣的相互作用，將復(fù)雜三維現(xiàn)象簡化成二維，從而使實驗現(xiàn)象更容易被觀察，而這種夾層以及流動的填充液，恰好能夠構(gòu)成Hele-Shaw流動狀態(tài)的形成條件。實驗所模擬的過程即為Hele-Shaw流動過程。在此基礎(chǔ)上，加入對于懸濁填充液以及氣體流動速率的控制方式，模擬不同環(huán)境和條件下的海床亞層氣體遷徙過程。可以實現(xiàn)Hele-Shaw流動過程的儀器被統(tǒng)稱為Hele-Shaw cell。

同時，實驗裝置還需能夠滿足注入填充液，保持氣密性，注射氣體。為此，本實驗裝置對主流的Hele-Shaw cell進(jìn)行了改良，以實現(xiàn)對海底油氣遷徙狀況的適配性調(diào)整：將空氣注入到二氧化硅懸濁液中，氣相流體進(jìn)入固液混合相，模擬二維平面流體在Hele-Shaw cell中的流動，觀察指進(jìn)現(xiàn)象，通過添加不同的懸濁液模擬不同海質(zhì)，以此來體現(xiàn)海底油氣過程的遷徙。

1.2 實驗裝置設(shè)計

完整的裝置由Hele-Shaw Cell模塊、電控模塊、監(jiān)測設(shè)備與圖形處理系統(tǒng)組成。

Hele-Shaw Cell模塊主要由兩塊平行的玻璃板組成（圖1）。將兩塊玻璃板（長w=350 mm，寬a=200mm，高h(yuǎn)=10 mm）通過玻璃膠粘接在一起。為了控制間隙，兩個厚度受控的硬質(zhì)塑料間隔件（厚度為1mm）被粘在兩個玻璃板之間。兩塊玻璃板中心形成流體通道，兩玻璃板一端封閉成V字形，另一端打開與大氣壓相通。同時在一層玻璃板的中心線上打通一個進(jìn)氣口。Hele-Shaw模型被一個鋁框架包圍，鋁框架固定在玻璃板上，硅密封件確保模型的密封性。使用流量計連接壓力緩沖器控制氣體流量，高像素攝像機(jī)記錄實驗現(xiàn)象。

圖1 Hele-Shaw Cell模塊裝配圖

在實驗時，將空氣通入壓力緩沖器中，用流量計連接壓力緩沖器與模型的進(jìn)氣口，在流量計與壓力緩沖器之間安裝電控模塊，該模塊是一個可以控制的閥門，用以調(diào)節(jié)氣體通斷，再將由流量計與數(shù)據(jù)采集器組成的檢測設(shè)備連接傳至電腦。圖像處理系統(tǒng)由相機(jī)及其組件和電腦組成，使用分辨率更高的相機(jī)，型號為MVXG2500GC-T，2500萬像素的攝像機(jī)安裝在距離機(jī)身30cm處的框架上，用以固定視野范圍，用LED光源板提供平行光源。攝像機(jī)與觸發(fā)信號發(fā)生器連接，再將觸發(fā)信號發(fā)生器與流量計連接，通過改變氣體流量的大小觸發(fā)攝像機(jī)，將Hele-Shaw cell模塊中二氧化硅顆粒通道的形成與變化記錄下來（實驗裝置如圖2所示）。

圖2 實驗裝置

使用該實驗裝置模擬不同情況下的海底油氣遷徙過程以及孔隙分布情況并直觀表現(xiàn)，嘗試復(fù)現(xiàn)海底油氣孔隙的形態(tài)并定性分析部分海床結(jié)構(gòu)特征。由于Hele-Shaw cell模型受到出油孔數(shù)目與孔距半徑的影響，氣體排出率隨出油孔數(shù)目的增大而增大，隨孔距半徑的增大而降低[16]，因此在此做出預(yù)測：流速越低，顆粒平均直徑越大，實驗指進(jìn)現(xiàn)象越明顯。

2 實驗設(shè)計與結(jié)果

2.1 實驗?zāi)康?/h3>

本實驗的目的是為了觀察氣體在多孔介質(zhì)下的遷徙狀態(tài)，以此來與數(shù)值模型進(jìn)行對比。海床是一個多孔結(jié)構(gòu)，油氣遷徙的過程，本質(zhì)上就是氣與液在多孔介質(zhì)中的移動過程。本項目組的實驗利用二氧化硅與鹽水溶液作為填充物模擬海床的多孔結(jié)構(gòu)，注入空氣模擬油氣。通過空氣的在多孔介質(zhì)的遷徙過程，使油氣與通道介質(zhì)的動力關(guān)系以及油氣在多孔介質(zhì)中的遷徙路徑可視化。

2.2 實驗設(shè)計

由于大部分海域天然氣水合物沉積物的覆蓋層沉積物主要由粘土（＜4 μm）和粉土（4～63 μm）和砂（＞63 μm）組成，其中粘土、砂占到了70%以上[17]，根據(jù)前人對于海洋土的研究可以得知海洋土的鈉、鎂離子較多并且普遍含有綠泥石、碳酸鹽砂層[18]，所以采用氯化鈉溶液與平均直徑分別為60 μm、100 μm二氧化硅顆粒。因為二氧化硅二氧化硅的的物理性質(zhì)、抗拉強(qiáng)度、含水率與沙土相同，其顆粒相較于其他顆粒更穩(wěn)定，更容易獲得，也更經(jīng)濟(jì)，且二氧化硅顆粒與飽和氯化鈉溶液混合時，對比與純凈水、酒精相混合，更容易在氯化鈉溶液中沉降，更容易出現(xiàn)顆粒與溶液分層，形成多孔介質(zhì)。因此選擇大小可控的二氧化硅顆粒與質(zhì)量分?jǐn)?shù)為5%氯化鈉溶液混合，來模擬海底巖土層與海水環(huán)境。

實驗選取了兩個控制變量，分別是二氧化硅的直徑（60 μm、100 μm）與注射速率（10 ml/min、5 ml/min）。二氧化硅直徑取值基于海底沉積顆粒的平均大小上下浮動，注射速率取值源于低雷諾數(shù)，即低流速、低密度、低特征長度和高粘度系數(shù)的邊界條件，此外后續(xù)調(diào)整量級變化時只需調(diào)整倍數(shù)，取值為5、10，以便能夠在制流量和圖形差異的圖像中獲得方形的直觀圖像。驅(qū)動氣體遷徙的動力學(xué)過程的主要因素是注射流速，理論上存在一個閾值流速，使得在流速閾值范圍內(nèi)氣體遷徙的指狀形狀和數(shù)量存在差異。為了方便后期計算，流速分別取0.5 ml/min、1 ml/min、5 ml/min、10 ml/min。

2.3 實驗步驟

（1）調(diào)配懸浮液。懸浮液主要成分是顆粒材料為平均直徑分別為60 μm、100 μm的二氧化硅顆粒與質(zhì)量分?jǐn)?shù)為3.5%NaCl溶液的混合體（二氧化硅與氯化鈉溶液的質(zhì)量比為2∶1），根據(jù)實驗需要調(diào)配體積分?jǐn)?shù)為48%的懸浮液。

（2）充分?jǐn)嚢杌旌弦?，將顆粒與氯化鈉混合物在磁力攪拌儀上攪拌3～5min至充分混合（懸濁液為均勻的乳白色狀）。

（3）添加懸浮液。從進(jìn)氣口通過注射器引入調(diào)配好的懸浮液，直到通道充滿顆粒懸浮液。

（4）注射氣體。注射泵連接注射器并從進(jìn)氣口注入氣體。注射泵可以調(diào)節(jié)注射的壓力大小和氣體流量。

（5）通過設(shè)定程序定時記錄實驗照片，使用電腦將照片比對觀察不同實驗現(xiàn)象，得出實驗結(jié)論。

2.4 實驗現(xiàn)象

第一組實驗選用平均直徑為60 μm的二氧化硅顆粒，調(diào)節(jié)注射泵注射速度為10 ml/min，實驗裝置在1 min20 s前無明顯變化，在1 min20 s時，注射口出現(xiàn)細(xì)小裂紋并不斷擴(kuò)張，但擴(kuò)張范圍十分有限，出現(xiàn)高度約4 cm的樹枝分支（圖3），該分支在短時間后與底部斷裂變成氣泡快速浮動上升。

圖3 第一組實驗結(jié)果

第二組實驗選用平均直徑為100 μm的二氧化硅顆粒，調(diào)節(jié)注射泵注射速度為10 ml/min，實驗裝置在2 min20 s前無明顯變化，2 min20 s時，氣體開始打通注射口處多孔介質(zhì)，出現(xiàn)指進(jìn)現(xiàn)象，并發(fā)展為類似樹枝結(jié)構(gòu)的遷徙通道，隨后該樹枝狀通道現(xiàn)象不斷成長至穩(wěn)定（圖4），分支多且細(xì)小。該結(jié)構(gòu)生長至距離玻璃板頂部約五分之一處停止生長。

圖4 第二組實驗結(jié)果

第三組實驗選用平均直徑為60 μm的二氧化硅顆粒，調(diào)節(jié)注射泵注射速度為5 ml/min，實驗裝置在2 min20 s前無明顯變化，在2 min20 s時，注射口出現(xiàn)細(xì)小裂紋并不斷擴(kuò)張，擴(kuò)張范圍相比2.2 ml/min情況下更明顯，但經(jīng)過短時間擴(kuò)張后，擴(kuò)張達(dá)到極限，分支停止衍生，經(jīng)過5 s后分支與裝置頂部空氣聯(lián)通。該組產(chǎn)生的樹枝分支面積較少，分支直徑較大（圖5）。

圖5 第三組實驗結(jié)果

第四組實驗選用平均直徑為100 μm的二氧化硅顆粒，調(diào)節(jié)注射泵注射速度為5 ml/min，實驗裝置在2 min20 s前無明顯變化，2 min20 s時，氣體開始打通注射口處多孔介質(zhì)，出現(xiàn)指進(jìn)現(xiàn)象，并發(fā)展為類似的樹枝結(jié)構(gòu)通道，隨后該樹枝狀通道現(xiàn)象不斷成長，形成兩個大小幾乎一樣的樹枝結(jié)構(gòu)，分支多且細(xì)小。左側(cè)結(jié)構(gòu)生長至距離玻璃板頂部約二分之一處停止生長，而右側(cè)結(jié)構(gòu)成長至聯(lián)通裝置頂部（圖6），整個聯(lián)通過程約18 s。

圖6 第四組實驗結(jié)果

2.5 實驗結(jié)果與分析

注射空氣，觀察空氣在填充物中的演變過程。注入的空氣以一種緩慢的、準(zhǔn)連續(xù)的指進(jìn)現(xiàn)象變化，在這個過程中，指進(jìn)通道側(cè)向分支和隨機(jī)生長方向產(chǎn)生了一個樹形結(jié)構(gòu)（圖7）。由前進(jìn)界面匯編而成的壓實顆粒前緣沿整個前緣可見為一條狹窄的暗帶。由于樹枝結(jié)構(gòu)內(nèi)部的空氣壓力與界面上的摩擦力之間的力平衡，會出現(xiàn)一個特征樹枝寬度。

圖7 樹狀分形形態(tài)

在不同條件下進(jìn)行了多次對比實驗，使用不同的填充液去模擬不同的海床亞層結(jié)構(gòu)，得到的油氣遷徙路徑都具有相似的樹狀結(jié)構(gòu)。以不同的注射速率進(jìn)行實驗，得到的樹狀通道雖然在寬度上有差異，但在形狀上都為樹狀。同時，不同實驗條件下樹狀結(jié)構(gòu)整體與局部上具有相似性。再者，為了避免由于操作人員的主觀性與事件偶然性對實驗結(jié)果的影響，獨(dú)立進(jìn)行了大量實驗，得到的結(jié)果都為相同。由此可以看出，油氣遷徙的樹狀結(jié)構(gòu)遷徙通道不受空間和時間的約束，符合分形理論[19]中的標(biāo)度不變性原則，證實了該實驗?zāi)Ｐ偷臏?zhǔn)確性和可行性。

實驗發(fā)現(xiàn)，樹枝結(jié)構(gòu)不斷演化的過程中隨著前端空氣被恒定的、緩慢的活塞驅(qū)動所壓縮，空氣壓力線性增加。隨著壓力的增加，前緣最薄弱的地方最終屈服，壓縮空氣在徑向擴(kuò)張之前通過狹窄的頸部迅速侵入。前端空氣減壓后，膨脹會減慢并停止。堆積在界面上的顆粒狀物質(zhì)沉淀下來，并逐漸壓縮成一個堵塞的前端，在這里，最薄弱的地方再次決定了下一個氣泡的位置。在不斷前進(jìn)、代替、壓縮、破裂過程中，空氣在填充物中不斷樹狀型前進(jìn)（如圖8所示）。這主要是因為壓力差為遷徙驅(qū)動力，壓力下降最快的方向為遷徙的優(yōu)先方向，豎向遷徙趨勢大于橫向遷徙趨勢[20]。

圖8 隨時間演化的粘性指狀結(jié)構(gòu)

在真實的油氣遷徙過程中往往以低速率進(jìn)行，低速率意味著靜態(tài)過程，樹枝周圍環(huán)境的影響可以忽略不計，從而有效地將樹枝與周圍環(huán)境隔開。因此，特征指標(biāo)與通道寬度無關(guān)，僅受到表面張力和材料屈服應(yīng)力的影響。屈服應(yīng)力和摩擦流體都分別受與屈服應(yīng)力和靜摩擦相關(guān)的閾限動力學(xué)控[21]。

2.5.1 注射速度對指進(jìn)現(xiàn)象的影響

在相同二氧化硅平均直徑的條件下，注射速度越慢，出現(xiàn)指進(jìn)現(xiàn)象所需要的注射量越少，在注射速度過快的情況下，甚至無法出現(xiàn)明顯指進(jìn)現(xiàn)象，僅僅是在注射口出現(xiàn)短暫細(xì)小的指進(jìn)現(xiàn)象，之后便成為一個個獨(dú)立的氣泡擴(kuò)散。將空氣注射速度增加到一個閾值以后，動態(tài)會發(fā)生顯著變化，出現(xiàn)不同的形態(tài)（圖9）。界面的發(fā)展不再緩慢和連續(xù)，相反，它的運(yùn)動變得高度間歇性，前端在很長一段時間內(nèi)保持靜止，隨后，粒狀流體混合物突然被氣泡形狀的空氣爆發(fā)所取代，位移繼續(xù)逐泡前進(jìn)。

圖9 蠕蟲形態(tài)

通過分析得知注射速度越大，氣體內(nèi)部的壓力越大。壓力增大，使氣體打破內(nèi)部壓力平衡，達(dá)到多孔介質(zhì)的屈服壓力，克服二氧化硅顆粒的摩擦力，克服毛細(xì)管力與邊緣粘滯力所需要的時間更短從而使樹狀結(jié)構(gòu)發(fā)生改變。

2.5.2 二氧化硅直徑對指進(jìn)現(xiàn)象的影響

通過對比多組實驗現(xiàn)象，最終發(fā)現(xiàn)平均直徑為100 μm的二氧化硅比平均直徑為60 μm的二氧化硅更容易產(chǎn)生更大的樹枝結(jié)構(gòu)。通過電鏡觀測多孔介質(zhì)結(jié)構(gòu)并分析多孔介質(zhì)顯微成像圖（圖10），得知增大二氧化硅顆粒大小，會導(dǎo)致多孔介質(zhì)中的空隙增大，較大的空隙更有利于空氣突破顆粒間的摩擦力，有更多的路徑能夠被輕易打破平衡，從而形成更多的樹枝狀分支。并且由于產(chǎn)生的樹枝分支數(shù)量多、面積大，氣體內(nèi)部壓力的增大不容易打破分支內(nèi)部壓力與外部張力、毛細(xì)管力之間的平衡，有利于分支結(jié)構(gòu)成長發(fā)展壯大，且不會使得通道突然斷裂形成氣泡。

該實驗對我國南海北部神狐海域天然氣儲藏海床結(jié)構(gòu)進(jìn)行了模擬并與該海域?qū)嶋H勘測的地震資料圖[22]進(jìn)行了對比。對比結(jié)果如圖11所示。

圖11 神狐海域結(jié)構(gòu)

從地震資料圖可以看出地震剖面上呈直立的、上小下大的煙囪狀通道，局部橫向擴(kuò)張呈囊狀、樹枝狀，底部結(jié)構(gòu)刺穿各沉積層最后聯(lián)通表面釋放油氣。通過實驗?zāi)M得到了與之相符合的結(jié)構(gòu)，實驗所得的樹枝狀結(jié)構(gòu)，從底部開始發(fā)育經(jīng)歷縱向、橫向發(fā)展的過程，最后與上部聯(lián)通形成上小下大的油氣遷徙通道。綜合對比，該模型成功模擬了此海域的油氣遷徙過程以及該海域的海床亞層結(jié)構(gòu)。

3 結(jié)語

本實驗裝置根據(jù)氣體和流體在兩平板間的流動特性，在Hele-Shaw模型和分形理論的基礎(chǔ)上設(shè)計制作，能模擬不同情況下的海底油氣遷徙過程以及孔隙分布情況并直觀表現(xiàn)。實驗表明：（1）該裝置能成功復(fù)現(xiàn)海底油氣遷徙孔隙的形態(tài)，如：樹狀形態(tài)、蠕蟲形態(tài)、粘性指形態(tài)等。（2）通過改變控制參數(shù)，該裝置能模擬不同海域不同海底亞層結(jié)構(gòu)狀態(tài)下油氣遷徙的路徑及發(fā)展，實驗結(jié)論對于我國海底油氣資源勘測及開采具有參考意義。（3）該實驗?zāi)転楹罄m(xù)對于油氣遷徙過程的分析研究提供理論研究依據(jù)和實驗設(shè)計思路，具有一定的實用意義和推廣價值。