張圣山，衛(wèi)云波，顧 熙

1. 湖州中交投資發(fā)展有限公司，湖州 31300；2. 河海大學(xué) 地球科學(xué)與工程學(xué)院，南京 210013；3. 南京大學(xué) 地球科學(xué)與工程學(xué)院，南京 210023

1 引言

排水和吸水過(guò)程是非飽和帶中兩個(gè)相反的水氣驅(qū)替過(guò)程。排水過(guò)程指的是多孔介質(zhì)中含水量減少（即氣驅(qū)水）的過(guò)程，而吸水過(guò)程指的是含水量增加（即水驅(qū)氣）的過(guò)程（Lenormand and Zarcone,1984）。而根據(jù)吸水過(guò)程中水相壓強(qiáng)（pw）與氣相壓強(qiáng)（pa）之間的大小關(guān)系，吸水過(guò)程又可以被細(xì)分為三類(lèi)：自吸（即 free imbibition，此時(shí)pw=pa）、受迫吸水（即forced imbibition，此時(shí)pw>pa）以及受限吸水（即constrained imbibition，此時(shí)（pw

因?yàn)椤捌款i效應(yīng)”的存在，孔隙尺度上的排水和吸水過(guò)程存在很大差別，如排水過(guò)程中存在不穩(wěn)定的“海恩斯跳躍”現(xiàn)象（Haines, 1930）、臨近孔隙干擾排水（Smiles et al., 1981）現(xiàn)象，以及吸水過(guò)程中存在的相鄰孔隙的“協(xié)同充水”（cooperative pore filling）現(xiàn)象（Hughes and Blunt, 2000）、毛細(xì)阻斷（snap off）現(xiàn)象（Berg et al., 2013）等等。這些孔隙尺度的排水/吸水現(xiàn)象上的差異是否會(huì)對(duì)達(dá)西尺度下的水分分布產(chǎn)生影響，文獻(xiàn)中尚存在爭(zhēng)議（Smiles et al., 1981; Ferrari and Lunati, 2014）。

本研究使用光透射法監(jiān)測(cè)了二維砂箱中含水量的分布，并定性和定量地分析了排水 吸水過(guò)程中含水量分布的特征差異，目的是通過(guò)室內(nèi)實(shí)驗(yàn)探討孔隙尺度下的水/氣驅(qū)替過(guò)程對(duì)于達(dá)西尺度下水分運(yùn)移過(guò)程的可能影響，并加深人們對(duì)于非飽和帶水氣驅(qū)替過(guò)程的理解。

2 實(shí)驗(yàn)設(shè)計(jì)

2.1 實(shí)驗(yàn)材料與裝置

實(shí)驗(yàn)中使用的介質(zhì)，為美國(guó)尤尼明公司（Unimin Corporation, Le Sueur, MN）生產(chǎn)的高純石英砂，商品名稱(chēng)Accusand。這種石英砂的SiO2含量大于92% 且有機(jī)質(zhì)含量很低，因此透光性能良好，適合用作光透法實(shí)驗(yàn)的介質(zhì)。實(shí)驗(yàn)所用的40/50目石英砂，顆粒密度2.66 g/cm3，粒徑介于40~50目（0.300~0.425 mm）之間，平均粒徑（d50）為0.359 mm，均一系數(shù)1.21（鄭菲，2015）。實(shí)驗(yàn)中使用的二維砂箱，如圖1a所示，內(nèi)部尺寸為55 cm×45 cm×2.5 cm。砂箱結(jié)構(gòu)由兩塊鋼化玻璃板夾一個(gè)中心鋁框構(gòu)成。玻璃板與中心鋁框之間使用硅酮密封膠黏合，以防止流體泄漏。中心鋁框的左右邊界完全封閉，上邊界設(shè)置了六個(gè)進(jìn)氣口與大氣相連通，下邊界設(shè)置了六個(gè)進(jìn)/出水口與蠕動(dòng)泵相連，以控制砂箱中的水位（圖1c）。為了盡可能避免出現(xiàn)水/氣流動(dòng)的優(yōu)勢(shì)通道，需采用濕法裝填砂箱：將石英砂逐層裝填進(jìn)砂箱，使用小錘輕擊砂面壓實(shí)，同時(shí)調(diào)整水位，以保證水面始終高于砂面1 cm左右（若水位過(guò)高，則石英砂在沉降過(guò)程中有可能發(fā)生分選，出現(xiàn)分層的現(xiàn)象）。在裝填過(guò)程中記錄石英砂的總質(zhì)量，結(jié)合石英砂顆粒密度及砂箱總體積可以計(jì)算得到實(shí)驗(yàn)所裝填砂箱的平均孔隙度。

圖1 實(shí)驗(yàn)裝置示意圖Fig. 1 The experimental setup (a) and side/front sketch view of the experimental setup (b, c)

2.2 光透法監(jiān)測(cè)系統(tǒng)

實(shí)驗(yàn)中使用光透法監(jiān)測(cè)系統(tǒng)（圖1b）記錄砂箱中的含水量變化。光透法監(jiān)測(cè)系統(tǒng)由燈箱、砂箱和CCD（Charge Coupled Device）相機(jī)組成。燈箱作為光源放在砂箱一側(cè)，其發(fā)出的光透射過(guò)砂箱以后，被放置在另一側(cè)的CCD相機(jī)所接收。為了避免外界光源對(duì)于實(shí)驗(yàn)結(jié)果的干擾，整個(gè)實(shí)驗(yàn)需要放置在暗箱中進(jìn)行。實(shí)驗(yàn)所用的暗箱尺寸約為2.5 m（長(zhǎng)）×1.0 m（寬）1.0 m（高）。燈箱作為暗箱中的唯一光源，放置在暗箱前端，距砂箱約20 cm。CCD相機(jī)（AP2E, ApogeeInstruments, Auburn,CA)放置在暗箱后端，距砂箱約1.8 m。相機(jī)通過(guò)數(shù)據(jù)線(xiàn)與計(jì)算機(jī)相連，使用相機(jī)控制軟件MaxIm DL（Diffraction limited, Ottwa）實(shí)時(shí)傳輸并記錄透過(guò)砂箱的光強(qiáng)變化情況。CCD相機(jī)所記錄的砂箱圖片分辨率為1200像素（長(zhǎng)）*980像素（寬），即實(shí)際砂箱中每1 cm長(zhǎng)度在監(jiān)測(cè)照片中約對(duì)應(yīng)22個(gè)像素。

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

2.3.1 連續(xù)排水 吸水循環(huán)

實(shí)驗(yàn)開(kāi)始前，將砂箱處于完全飽水狀態(tài)。因此在實(shí)驗(yàn)開(kāi)始前，CCD相機(jī)所記錄的光強(qiáng)值，可以視作為完全飽水下的光強(qiáng)值Iws。完全飽水時(shí)，水所占的體積與孔隙體積相等，因此相應(yīng)的飽和體積含水量θws就等于孔隙體積與總體積之比。此處體積含水量θw定義為：

式中，θw為體積含水量；V水為水的體積；V孔隙為孔隙的體積；V顆粒為顆粒的體積。飽水光強(qiáng)記錄完成后，打開(kāi)蠕動(dòng)泵，控制砂箱經(jīng)歷三次排水 吸水的循環(huán)。首先是初次排水，以10 mL/min的流量從砂箱底部向外抽水，歷時(shí)150分鐘。初次排水結(jié)束以后靜置50分鐘，避免由于邊界條件突然逆轉(zhuǎn)造成的砂箱內(nèi)含水量波動(dòng)。靜置結(jié)束后進(jìn)行第一次吸水，使用蠕動(dòng)泵以10 mL/min的速度從砂箱底部向內(nèi)注水，歷時(shí)100分鐘。吸水結(jié)束以后，靜置50分鐘以完成第一次排水 吸水循環(huán)。隨后，以相同的步驟又進(jìn)行了兩次排水 吸水的循環(huán)，流量與歷時(shí)設(shè)置見(jiàn)表1。實(shí)驗(yàn)中，CCD相機(jī)的拍照間隔設(shè)置為1分鐘，每次拍照的曝光時(shí)間設(shè)置為6 s。從完全飽水/干燥的條件下開(kāi)始排水/吸水，稱(chēng)為初次排水/吸水。

表1 排水 吸水循環(huán)實(shí)驗(yàn)條件Table 1 Experimental conditions of drainage and imbibition

2.3.2 不同流量下的吸水 排水實(shí)驗(yàn)

為了研究流速對(duì)于含水量分布特征的影響，參照前人室內(nèi)及野外實(shí)驗(yàn)的流速設(shè)置（Lehmann et al., 1998; Cartwright, 2014; Shoushtari et al., 2018），在同一個(gè)砂箱中使用三種其他流量（5 mL/min，100 mL/min，400 mL/min）又做了三次吸水 排水實(shí)驗(yàn)（表2）。在每一次吸水 排水實(shí)驗(yàn)開(kāi)始之前，首先將砂箱中的水盡可能排干并靜置 50分鐘。接著使用設(shè)定的流量從砂箱底部向內(nèi)注水至飽和，隨后靜置 50分鐘并使用同樣的流量從砂箱底部向外排水，直到將砂箱排干。對(duì)于5 mL/min 流量下的實(shí)驗(yàn)，CCD 相機(jī)的拍照間隔設(shè)置為1分鐘，每張照片的曝光時(shí)間設(shè)置為 6 s。而對(duì)于100 mL/min和400 mL/min流量下的實(shí)驗(yàn)，CCD 相機(jī)的拍照間隔設(shè)置為3 s，每張照片的曝光時(shí)間設(shè)置為1 s。

表2 不同流量下吸水 排水實(shí)驗(yàn)條件Table 2 Experimental conditions of drainage and imbibition at different flow rates

2.4 分析方法及精度評(píng)價(jià)

CCD相機(jī)拍攝到的光強(qiáng)數(shù)據(jù)可以使用Niemet和Selker（2001）提出的方法轉(zhuǎn)換成飽和度數(shù)據(jù)，其基本原理在文獻(xiàn)中已有闡述（鄭菲，2015；章艷紅等，2014），本文在此就不做介紹。為了評(píng)估該光透法監(jiān)測(cè)系統(tǒng)的測(cè)量精度，可以對(duì)比實(shí)際流入/流出砂箱的水分體積與光透法監(jiān)測(cè)到的砂箱內(nèi)水分體積的變化，如圖2。這里使用均方根誤差（RMSE）來(lái)表征實(shí)驗(yàn)過(guò)程中測(cè)量體積與實(shí)際體積的平均誤差，即

圖2 排水吸水循環(huán)實(shí)驗(yàn)中，實(shí)際流入/流出砂箱的水體積（actual volume）與光透法監(jiān)測(cè)到的水體積（measured volume）變化對(duì)比Fig. 2 Comparison of the actual volume with measured volume of inflow/outflow water

式中，N為測(cè)量次數(shù)；Vactu為實(shí)際流入/流出體積；Vmeas為光透法測(cè)得的砂箱內(nèi)水分體積變化。計(jì)算得到的RMSE值見(jiàn)表3，可見(jiàn)光透法監(jiān)測(cè)數(shù)據(jù)的整體誤差約在5%附近，最大不超過(guò)7%。這部分誤差為光透法監(jiān)測(cè)系統(tǒng)的系統(tǒng)誤差，來(lái)源可能有很多方面，比如由于不透明邊界的阻擋，貼近邊界部分的介質(zhì)飽和度無(wú)法被準(zhǔn)確測(cè)量。另外，由于石英砂是自下而上逐層填入砂箱，這就導(dǎo)致砂箱下部的石英砂比砂箱上部更加密實(shí)一點(diǎn)。砂箱上下密實(shí)度的差異也會(huì)給光透法的測(cè)量引入系統(tǒng)誤差(在計(jì)算飽和度及含水量時(shí)我們假設(shè)整個(gè)砂箱是均質(zhì)的)。除了系統(tǒng)誤差，光透法的測(cè)量結(jié)果也會(huì)受到隨機(jī)誤差的影響。光透法的隨機(jī)誤差主要由電噪聲、光子散粒噪聲、砂顆粒顏色差異等因素造成（Azooz et al., 1996）。與系統(tǒng)誤差不同，由于隨機(jī)誤差會(huì)在求和時(shí)互相抵消，因此不會(huì)對(duì)整體的水量平衡造成影響。隨機(jī)誤差會(huì)造成單點(diǎn)測(cè)量結(jié)果的波動(dòng)，因此影響到光透法測(cè)量結(jié)果的單點(diǎn)精度。這里使用殘余水區(qū)域以及完全飽水區(qū)域的含水量波動(dòng)標(biāo)準(zhǔn)差 (±2%S.D.)來(lái)評(píng)估光透法測(cè)量結(jié)果的隨機(jī)誤差。在殘余水區(qū)域以及完全飽水區(qū)域，含水量一般被認(rèn)為是常數(shù)，分別等于θwr以及θws。因此這兩個(gè)區(qū)域中含水量的波動(dòng)可以被歸因于光透法測(cè)量的隨機(jī)誤差。實(shí)驗(yàn)過(guò)程中，光透法測(cè)得殘余水以及完全飽水區(qū)域的含水量波動(dòng)標(biāo)準(zhǔn)差見(jiàn)表3?？梢钥吹?，在實(shí)驗(yàn)過(guò)程中隨機(jī)誤差造成的測(cè)量結(jié)果波動(dòng)約在2%附近，最大不超過(guò)3%，因此其對(duì)測(cè)量結(jié)果的影響基本可以被忽略。

表3 光透法監(jiān)測(cè)結(jié)果的均方根誤差及隨機(jī)誤差Table 3 Comparison of the RMSEs and random errors during the experiments

3 實(shí)驗(yàn)結(jié)果與分析討論

3.1 排水 吸水過(guò)程含水量分布對(duì)比

圖3顯示了排水 吸水循環(huán)實(shí)驗(yàn)中六個(gè)階段的含水量分布曲面。從圖3a（初次排水）中可以看到，在砂箱的不同高度處，含水量的分布呈現(xiàn)出不同的規(guī)律。根據(jù)含水量的分布規(guī)律，可以將曲面分為三個(gè)部分。在圖中的I號(hào)區(qū)域（即濕潤(rùn)鋒以上的區(qū)域），含水量值小于 0.10，區(qū)域內(nèi)的含水量基本沒(méi)有橫向波動(dòng)，曲面相對(duì)平滑。而當(dāng)含水量介于 0.10~0.25之間時(shí)（區(qū)域II，即濕潤(rùn)鋒所在的區(qū)域），含水量出現(xiàn)了比較明顯的橫向波動(dòng)，曲面變得不再平滑，并能觀察到“指狀”的起伏。當(dāng)含水量增加到0.25以上時(shí)（區(qū)域III，即濕潤(rùn)鋒以下的區(qū)域），含水量的橫向波動(dòng)又變得不明顯起來(lái)。在圖3b、c（第二和第三次排水）的濕潤(rùn)鋒上也可以觀察到“指狀”起伏；但是對(duì)于吸水過(guò)程，在圖3b、d、f 中，類(lèi)似的“指狀”起伏卻不是特別明顯。排水過(guò)程和吸水過(guò)程中濕潤(rùn)鋒上的這種含水量分布上的差異可以被歸結(jié)于以下原因：排水過(guò)程是一個(gè)相對(duì)不穩(wěn)定的過(guò)程，由于孔隙尺度下海恩斯跳躍以及臨近孔隙干擾排水這些瞬變現(xiàn)象的存在，排水過(guò)程介質(zhì)中含水量的波動(dòng)會(huì)大于相對(duì)穩(wěn)定的吸水過(guò)程。

圖3 含水量分布曲面Fig. 3 The moisture content distribution for each period

非穩(wěn)定的排水過(guò)程，即孔隙尺度上的海恩斯跳躍現(xiàn)象，很早就有人研究（Klimentos and Mccann,1990）。海恩斯跳躍一般發(fā)生在氣體驅(qū)替水分從孔喉進(jìn)入孔腹之時(shí)。在排水過(guò)程中，由于毛細(xì)力的作用，水氣交界的彎液面往往不會(huì)均勻地流過(guò)孔隙，而是直接從所處的孔喉迅速地躍變到下一個(gè)孔喉處。這種躍變往往也伴隨著毛細(xì)壓力的瞬間下降（Reynolds et al., 2009），造成相鄰孔隙里水分分布的不穩(wěn)定波動(dòng)。Rücker等（2009）指出，在達(dá)西尺度下的排水過(guò)程，可以視作是一系列孔隙尺度下“雪崩式”不停發(fā)生的海恩斯跳躍事件的集合。從圖3中排水與吸水過(guò)程濕潤(rùn)鋒的對(duì)比可以推測(cè)，孔隙尺度下的不穩(wěn)定排水過(guò)程影響到了濕潤(rùn)鋒上的含水量分布，并且這類(lèi)影響能夠在達(dá)西尺度下被觀察到。

為了更進(jìn)一步研究排水與吸水過(guò)程中含水量波動(dòng)的差異，將濕潤(rùn)鋒上含水量的橫向分布提取出來(lái)加以對(duì)比。圖4顯示的是第二次排水與吸水過(guò)程中，濕潤(rùn)鋒上含水量橫向分布的對(duì)比。從圖4可以看出來(lái)，不管是排水過(guò)程還是吸水過(guò)程，其含水量的橫向分布都存在一定的波動(dòng)。在排水過(guò)程中含水量的波動(dòng)范圍更大（0.145~0.175），而吸水過(guò)程中含水量的波動(dòng)范圍（0.150~0.165）僅為排水過(guò)程的 1/2。排水過(guò)程中濕潤(rùn)鋒上含水量的大范圍波動(dòng)，可以視作是導(dǎo)致圖3中“指狀”起伏出現(xiàn)的直接因素。

圖4 第二次排水與吸水過(guò)程中，濕潤(rùn)鋒上含水量橫向分布的對(duì)比Fig. 4 Comparison of the lateral moisture distribution during main drainage and main imbibition. The insets show the locations in the sand chamber where the moisture content data is extracted

3.2 含水量空間分布相關(guān)性分析

圖5中顯示的是第二、第三次排水/吸水過(guò)程中濕潤(rùn)鋒上含水量橫向分布的對(duì)比。由圖5a可以發(fā)現(xiàn)，在第二次排水與第三次排水過(guò)程中，含水量波動(dòng)的趨勢(shì)非常相似。同樣的，在圖5b中，第二次吸水與第三次吸水的含水量波動(dòng)趨勢(shì)也比較一致。根據(jù)圖5中的數(shù)據(jù)，計(jì)算了這四個(gè)階段中濕潤(rùn)鋒上含水量橫向分布的斯皮爾曼相關(guān)系數(shù)r （表4）。從表4可知，第二與第三次排水過(guò)程的含水量分布顯示出極顯著的強(qiáng)相關(guān)性（r=0.71，p<0.01），同樣的，第二與第三次吸水過(guò)程的含水量分布也顯示極顯著的強(qiáng)相關(guān)性（r=0.79，p<0.01）。

表4 不同排水/吸水過(guò)程，濕潤(rùn)鋒上含水量分布的斯皮爾曼相關(guān)系數(shù)Table 4 Spearman correlation coefficients (r) for the moisture distribution during drainage and imbibition.

圖5 (a) 第二、第三次排水過(guò)程及 (b) 第二、第三次吸水過(guò)程中，濕潤(rùn)鋒上含水量分布對(duì)比Fig. 5 Comparison of the moisture distribution on the wetting front at the 50th min of (a) main drainage and secondary drainage,(b) main imbibition and secondary imbibition

圖5和表4中不同排水(或吸水)過(guò)程含水量分布的強(qiáng)相關(guān)性可以被歸因于石英砂介質(zhì)本身的空間分布特征。Suekane等（2009）使用X射線(xiàn)層析成像法研究了砂巖巖芯以及玻璃珠介質(zhì)中壓實(shí)程度與含水量分布的關(guān)系，并發(fā)現(xiàn)在平衡狀態(tài)下壓實(shí)度越大的地方含水量越大。Zhou等（2010）在隨后的研究中指出，平衡狀態(tài)下的含水量分布與介質(zhì)孔隙度分布呈現(xiàn)較強(qiáng)的相關(guān)性，并且這種相關(guān)性不隨實(shí)驗(yàn)流速的改變而改變。雖然本次實(shí)驗(yàn)砂箱內(nèi)的局部孔隙度分布難以測(cè)量，但是從不同排水（或吸水）過(guò)程含水量分布的強(qiáng)相關(guān)性可以推測(cè)，實(shí)驗(yàn)中觀察到的含水量波動(dòng)并不是隨機(jī)出現(xiàn)的，而是受到了介質(zhì)本身的空間分布特征（局部孔隙度分布/局部壓實(shí)度差異）的控制。此外，從表4中還可以看到，排水過(guò)程含水量分布與吸水過(guò)程含水量分布之間的相關(guān)性較弱，說(shuō)明含水量的分布不僅與介質(zhì)本身的空間分布特征有關(guān)，還會(huì)受到水分運(yùn)動(dòng)方向的影響。

3.3 含水量分布的Hurst指數(shù)

為了檢驗(yàn)砂箱中含水量的橫向分布是否符合分?jǐn)?shù)布朗運(yùn)動(dòng)（fractional Brownian motion，fBm）模型分布，在砂箱的不同高度處提取出含水量的橫向分布，并分別使用差值法和V/S分析法作雙對(duì)數(shù)散點(diǎn)圖。圖6顯示了三個(gè)不同位置的分析示例。從圖6中可以發(fā)現(xiàn)，不管是差值法還是V/S分析法得到的散點(diǎn)圖中，窗口長(zhǎng)度（?x）的對(duì)數(shù)與該窗口上統(tǒng)計(jì)量（|?θw|或V/S2）的對(duì)數(shù)均呈較好的線(xiàn)性關(guān)系，并且兩種方法估計(jì)得到的H值非常接近。由此可知，fBm可以很好地描述砂箱中含水量的橫向波動(dòng)。

圖6 (a) 差值法及 (b) V/S分析法估計(jì)的Hurst指數(shù)值的對(duì)比Fig. 6 Comparison of Hurst coefficients estimated by (a) Intensity Difference method and (b) V/S analyses

圖7顯示了每個(gè)排水/吸水過(guò)程，砂箱中不同高度處含水量橫向分布的Hurst指數(shù)值。由圖7可見(jiàn)，不管是排水還是吸水過(guò)程，Hurst指數(shù)的峰值都出現(xiàn)在濕潤(rùn)鋒上，排水過(guò)程的峰值H在0.5附近，而吸水過(guò)程的峰值H只有0.2左右。H值越小意味著含水量更傾向于局限在均值附近波動(dòng)，因此波動(dòng)范圍越??；相對(duì)的，更大的H值意味著含水量的波動(dòng)范圍也更大。H值在濕潤(rùn)鋒上達(dá)到峰值，意味著濕潤(rùn)鋒上的含水量分布波動(dòng)范圍最大，這一點(diǎn)與我們?cè)趫D3中觀察到的現(xiàn)象一致。另外，吸水過(guò)程中0.2的H峰值顯示了在吸水過(guò)程中含水量的波動(dòng)受到了較強(qiáng)的抑制，其分布有強(qiáng)烈的回歸均值的趨勢(shì)；而排水過(guò)程中0.5的H峰值則說(shuō)明濕潤(rùn)鋒上含水量的波動(dòng)更接近隨機(jī)波動(dòng)，相比吸水過(guò)程受到的抑制較少，更接近一個(gè)不穩(wěn)定的過(guò)程。Moebius和Or（2012）在他們的孔隙尺度排水實(shí)驗(yàn)中指出，達(dá)西尺度下看似連續(xù)運(yùn)動(dòng)的濕潤(rùn)鋒實(shí)際上在孔隙尺度上看來(lái)是非常不規(guī)則和不穩(wěn)定的。從圖7中可以看到，使用fBm模型及Hurst指數(shù)可以較好刻畫(huà)這種不穩(wěn)定排水現(xiàn)象造成濕潤(rùn)鋒上含水量波動(dòng)的精細(xì)結(jié)構(gòu)。

圖7 排水 吸水過(guò)程中砂箱不同高度處含水量分布的Hurst指數(shù)值Fig. 7 Comparison of Hurst coefficients in the sand chamber during each period

3.4 不同流量條件下含水量分布的Hurst指數(shù)

前文提到，較大的Hurst指數(shù)與孔隙尺度不穩(wěn)定排水過(guò)程有關(guān)。因此，如果增大排水流量，加快排水速度，應(yīng)該也能觀察到 Hurst 指數(shù)的升高。圖8中顯示了不同流量條件下的Hurst指數(shù)對(duì)比。從圖8中可以看到，排水過(guò)程中隨著流量從5 mL/min升高到400 mL/min，濕潤(rùn)鋒區(qū)域（區(qū)域II）的H值變化很小，而濕潤(rùn)鋒以上區(qū)域（區(qū)域I）的H值變化較為明顯，尤其是當(dāng)流量增加到400 mL/min時(shí)，區(qū)域I的H值從0升高到了0.3附近。區(qū)域II的H值變化不明顯是因?yàn)樵搮^(qū)域H峰值本來(lái)就已經(jīng)很高，在0.5左右的H值意味著含水量的波動(dòng)已經(jīng)接近隨機(jī)波動(dòng)，因此排水流量的進(jìn)一步增加無(wú)法造成區(qū)域II中H值的增加。而區(qū)域I中H值的明顯升高可以被歸因于快速排水造成的不穩(wěn)定現(xiàn)象。隨著排水流量的增加，越來(lái)越多的水分被滯留于濕潤(rùn)鋒上方來(lái)不及及時(shí)排出。由于這部分滯留的水分無(wú)法在濕潤(rùn)鋒上方穩(wěn)定存在，會(huì)隨著排水過(guò)程的進(jìn)行逐步排出。正是這些仍然在緩慢排出的滯留水分造成了區(qū)域I中Hurst指數(shù)的升高。另一方面，觀察吸水過(guò)程的Hurst指數(shù)可以發(fā)現(xiàn)，隨著吸水流量增大，吸水速度加快，計(jì)算得到的Hurst指數(shù)反而出現(xiàn)了些許下降。從圖8中可以看到，隨著吸水流量從5 mL/min上升到400 mL/min，砂箱中Hurst指數(shù)的峰值從0.2左右下降到了0.1左右。H峰值的下降意味著隨著吸水流量升高，濕潤(rùn)鋒上含水量的波動(dòng)有更強(qiáng)烈的均值回歸的傾向，其波動(dòng)范圍受到了抑制。關(guān)于吸水速度對(duì)于濕潤(rùn)鋒形態(tài)的影響，前人也做了不少研究（Berg et al., 2013, Nguyen et al.,2006）。Hughes和Blunt指出，在吸水速度較快的時(shí)候，粘性力占主導(dǎo)而毛細(xì)力受到抑制，因此水分更傾向于優(yōu)先充填臨近的孔隙，而不是在毛細(xì)力的作用下上升形成指狀毛細(xì)水（Berg et al., 2013）。此時(shí)相鄰孔隙的“協(xié)同充水”現(xiàn)象發(fā)生得很頻繁，在孔隙尺度下來(lái)看，水氣交界面往往是連續(xù)且平滑的。而當(dāng)吸水速度較低時(shí)，毛細(xì)力占據(jù)了主導(dǎo)地位，在毛細(xì)力的作用下水分沿著孔隙直徑較小的路徑上升，形成了一處處不規(guī)則的指狀毛細(xì)水。根據(jù)圖8中的結(jié)果可以推測(cè)，這種吸水速度對(duì)于水氣交界面的影響同樣能在濕潤(rùn)鋒上體現(xiàn)出來(lái)。由于介質(zhì)本身存在一定的非均質(zhì)性，低流量下濕潤(rùn)鋒上的水分在毛細(xì)力的作用下不規(guī)則上升，體現(xiàn)為含水量在橫向上的波動(dòng)。然而隨著吸水速度加快，毛細(xì)力的作用受到抑制，濕潤(rùn)鋒上的水分分布變得更加規(guī)則，所以此時(shí)含水量的波動(dòng)也受到了抑制，H值也隨之下降。

圖8 不同流量條件下吸水/排水過(guò)程含水量分布Hurst指數(shù)對(duì)比Fig. 8 Comparison of Hurst coefficients under different flow rates

4 結(jié)論

通過(guò)二維砂箱實(shí)驗(yàn)研究了排水及吸水過(guò)程中含水量分布的差異，并使用分?jǐn)?shù)布朗運(yùn)動(dòng)模型定量分析和討論了含水量的波動(dòng)特征。在排水過(guò)程中，由于“海恩斯跳躍”等孔隙尺度不穩(wěn)定排水現(xiàn)象的存在，導(dǎo)致含水量在橫向上出現(xiàn)較大范圍的波動(dòng)，并且在含水量曲面上能觀察到含水量波動(dòng)造成的指狀起伏。而在吸水過(guò)程中，含水量的波動(dòng)范圍更小，含水量曲面也更為平整。含水量分布在排水過(guò)程和吸水過(guò)程之間的差異說(shuō)明了孔隙尺度的一些常見(jiàn)過(guò)程能夠?qū)_(dá)西尺度下的水分運(yùn)動(dòng)過(guò)程造成影響，而不是在提升尺度的過(guò)程中被簡(jiǎn)單抵消掉。