床面形態(tài)滲透率對(duì)潛流交換特性的影響

趙 亮 王學(xué)勇 王岱峰 樊靖郁

(1.上海大學(xué)力學(xué)與工程科學(xué)學(xué)院, 上海 200444;2.上海大學(xué)上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所, 上海 200072)

潛流交換是河流水環(huán)境系統(tǒng)中, 各種物質(zhì)(水、溶質(zhì)、污染物、膠體和細(xì)顆粒等) 在與地表水相鄰的飽和沉積物層(潛流帶) 中發(fā)生的一種重要對(duì)流過程[1-3].潛流交換對(duì)潛流帶中溶解氧、氮/磷營養(yǎng)鹽、重金屬、有機(jī)污染物以及膠體、細(xì)顆粒物(如微塑料) 的遷移轉(zhuǎn)化和通量演變起控制和調(diào)節(jié)作用[4-5], 影響潛流帶的生物地球化學(xué)特征[6-8]對(duì)河流生態(tài)系統(tǒng)的健康維持和功能發(fā)揮起著非常重要的作用[4,9].

潛流交換可由不同尺度和類型的水靜力/水動(dòng)力因子驅(qū)動(dòng), 在河道局部尺度上各種類型的床面形態(tài)是潛流交換的重要驅(qū)動(dòng)因子[10-11].上覆水流與床面形態(tài)的相互作用會(huì)使床面壓力變化不均勻, 驅(qū)動(dòng)上覆水及其所攜帶的各種物質(zhì)產(chǎn)生進(jìn)出底床的潛流交換[1,10].床面形態(tài)驅(qū)動(dòng)的潛流交換涉及到眾多影響因素.已有研究表明, 水動(dòng)力作用(流量或流速/水深、非恒定性)、床沙組成和床面形態(tài)物理特性(滲透率、孔隙度和非均質(zhì)性) 對(duì)潛流交換的通量、潛流路徑和停留時(shí)間的分布具有重要影響[12-14].

非均質(zhì)底床結(jié)構(gòu)在天然河流中是普遍存在的.近年來河床的非均質(zhì)性在不同的時(shí)空尺度上對(duì)潛流交換的影響日益受到關(guān)注[15-17], 如不同粒徑床沙組成底床[18]、分層底床[19]以及低滲透率透鏡體[20]等.對(duì)于存在床面形態(tài)的潛流交換, 本工作側(cè)重于考慮河道局部尺度上的非均質(zhì)性、自然或人工形成的床面形態(tài)組成單元, 由于其物理特性與其下底床往往存在明顯差異,如砂質(zhì)河道中自然存在的倒伏樹木[21]、肋狀/簇狀礫石堆積體[11]等障礙物, 人工構(gòu)建的圓木壩(坎、堰)[22]、階梯-深潭[23]等潛流型河床結(jié)構(gòu), 因此這些床面形態(tài)組成單元的物理特性(如滲透率、孔隙度) 與底床泥沙相比呈現(xiàn)出河道局部尺度上的空間變化.即使對(duì)于存在床面形態(tài)(如沙波) 的均質(zhì)底床結(jié)構(gòu), 由于上覆水體中細(xì)顆粒泥沙的局部沉積作用[24], 或床面形態(tài)泥沙顆粒表面生物膜生長形成的生物堵塞效應(yīng)[25], 故床面形態(tài)組成單元的滲透率會(huì)明顯改變, 形成空間滲透率強(qiáng)烈不連續(xù)變化的非均質(zhì)底床結(jié)構(gòu).對(duì)于這種類型的非均質(zhì)底床結(jié)構(gòu), 床面形態(tài)滲透率的變化對(duì)近底區(qū)域上覆水/孔隙水耦合流動(dòng)結(jié)構(gòu)具有顯著的影響[26], 在與之相關(guān)的潛流交換特性方面, 目前較少有學(xué)者研究床面形態(tài)組成單元滲透率的相對(duì)變化對(duì)潛流交換通量所帶來的定量差異, 因此床面形態(tài)組成單元的物理特性(如滲透率、孔隙度) 對(duì)潛流交換特性的影響規(guī)律尚待深入研究[27-28].

本工作通過實(shí)驗(yàn)室環(huán)形水槽實(shí)驗(yàn), 采用概化的2 維離散床面形態(tài), 主要研究存在床面形態(tài)砂質(zhì)底床條件下, 床面形態(tài)組成單元與底床滲透率的相對(duì)變化與潛流交換通量之間的定量關(guān)聯(lián)和變化特征, 并采用參數(shù)化方法分析床面形態(tài)組成單元的物理特性(滲透率變化) 對(duì)有效擴(kuò)散系數(shù)的影響趨勢.

1 實(shí)驗(yàn)裝置和測量方法

1.1 實(shí)驗(yàn)裝置

本實(shí)驗(yàn)在上海大學(xué)力學(xué)所自行研發(fā)的環(huán)形水槽中進(jìn)行[28-29], 實(shí)驗(yàn)裝置如圖1 所示.環(huán)形水槽由有機(jī)玻璃制成的內(nèi)、外徑不同的環(huán)形槽道, 以及剪力環(huán)、操作臺(tái)、電機(jī)、減速機(jī)、底座等構(gòu)成, 槽道內(nèi)、外徑分別為0.6 和1.0 m, 槽道寬度為0.2 m, 高度為0.4 m, 有效體積約為0.7 m3.電機(jī)驅(qū)動(dòng)安裝于水槽上部的剪力環(huán)轉(zhuǎn)動(dòng), 帶動(dòng)水體表面產(chǎn)生單向流動(dòng).通過調(diào)節(jié)剪力環(huán)的轉(zhuǎn)速(R) 可以控制上覆水的平均流速(U).

圖1 環(huán)形水槽示意圖Fig.1 Schematic diagram of an annular flume

在實(shí)驗(yàn)過程中剪力環(huán)轉(zhuǎn)速控制在低速范圍(R≤20 r/min), 以避免底床泥沙的再懸浮, 同時(shí)可近似忽略水槽內(nèi)形成的2 次流.每個(gè)實(shí)驗(yàn)工況上覆水的水深(H= 15.0 cm) 和體積(Vw)均保持不變.與循環(huán)直水槽相比, 環(huán)形水槽的優(yōu)點(diǎn)在于把直槽長度轉(zhuǎn)化為時(shí)間尺度, 可模擬斷面水流狀態(tài)相同的無限長的明渠水流, 同時(shí)無需回水系統(tǒng)和消能措施, 結(jié)構(gòu)簡單且易于操作.此外, 環(huán)形水槽上覆水體積與底床泥沙體積(Vs) 之比可方便調(diào)節(jié), 在實(shí)驗(yàn)過程中上覆水溶質(zhì)的濃度變化范圍相對(duì)較大, 可有效提高上覆水溶質(zhì)濃度的測量精度.

1.2 測量方法

實(shí)驗(yàn)底床泥沙按粒徑(dg) 大小分為3 種: 細(xì)沙、中沙和粗沙.3 種泥沙樣品采用粒徑較為均勻(粒徑分布范圍較窄) 的模型沙(石英砂), 通過篩分法測量3 種泥沙的粒徑分別為細(xì)沙(dg= 0.425 mm)、中沙(dg= 0.950 mm) 和粗沙(dg= 2.000 mm).通過水蒸發(fā)法對(duì)泥沙樣品的孔隙度進(jìn)行測定, 結(jié)果分別為0.42(細(xì)沙)、0.44(中沙) 和0.45(粗沙), 滲透率(K) 可根據(jù)Kozeny-Carmen 公式估算[1,29], 相應(yīng)的滲透率值分別為2.21×10-6cm2(細(xì)沙)、1.36×10-5cm2(中沙) 和6.69×10-5cm2(粗沙)[29].在實(shí)驗(yàn)過程中, 將15 cm 厚的模型沙(分別采用中沙和粗沙) 均勻鋪設(shè)在水槽底部構(gòu)成平整底床, 底床泥沙厚度(db)、底床表面積(As) 和體積(Vs=Asdb) 均保持不變.床面形態(tài)采用概化的2 維離散床面形態(tài), 包括人工構(gòu)建的具有不同滲透率的方形2 維離散床面形態(tài)和人工鋪設(shè)的沙波床面形態(tài), 其幾何特征分別如圖2 和3 所示.由不同滲透率床面形態(tài)組成單元的高度(h= 2.0 cm) 保持相同, 通過其間距(w) 的變化將不同數(shù)量的沙袋或沙波均勻排列在水槽底床表面來實(shí)現(xiàn)(3 ≤w/h≤11).

圖2 方形2 維離散床面形態(tài)示意圖Fig.2 Schematic diagram of square two-dimensional discrete bedform

圖3 沙波床面形態(tài)示意圖Fig.3 Schematic diagram of sand wave bedform

采用非吸附性氯化鈉(NaCl) 作為代表性溶質(zhì)(上覆水中分子擴(kuò)散系數(shù)為Dm).實(shí)驗(yàn)開始時(shí), 待水槽啟動(dòng)穩(wěn)定運(yùn)行一段時(shí)間后, 將配制的高濃度NaCl 溶液添加于上覆水中(底床孔隙水初始濃度為0), 混合均勻后的上覆水溶質(zhì)濃度作為初始濃度(C0), 通過虹吸采樣和測量得到上覆水NaCl 濃度隨時(shí)間(t) 的變化(C(t)), 即可確定不同工況的潛流交換通量(J)[28].在實(shí)驗(yàn)過程中, 為避免上部剪力環(huán)轉(zhuǎn)動(dòng)對(duì)采樣的擾動(dòng), 采樣點(diǎn)設(shè)置于中間水深位置(見圖1), 測量結(jié)果表明此區(qū)域上覆水溶質(zhì)濃度的垂向分布較為均勻.為便于比較, 本實(shí)驗(yàn)不同工況的上覆水溶質(zhì)濃度C(t) 的實(shí)驗(yàn)數(shù)據(jù)由量綱一變量C?(C?=C/C0) 的時(shí)間序列給出.有效擴(kuò)散系數(shù)的分析采用Deff/D′m量綱一形式, 其中D′m(D′m=Dm/(1+3(1+θ))) 為孔隙水中NaCl 分子擴(kuò)散系數(shù).根據(jù)上覆水量綱一溶質(zhì)濃度C?～t1/2曲線的初始斜率S(S= dC?/dt1/2), 即可定量確定有效擴(kuò)散系數(shù)Deff= π(VwS/(2Asθ))2, 由實(shí)驗(yàn)數(shù)據(jù)擬合得到[29].每個(gè)實(shí)驗(yàn)工況上覆水初始濃度(C0=0.05 mol/L) 和采樣歷時(shí)(180 min) 均保持相同.

1.3 實(shí)驗(yàn)工況

本實(shí)驗(yàn)通過改變床面形態(tài)構(gòu)成單元以及底床泥沙組成來實(shí)現(xiàn)滲透率的相對(duì)變化.本實(shí)驗(yàn)共設(shè)置了30 組實(shí)驗(yàn)工況, 其中平整底床分別由中沙和粗沙構(gòu)成, 存在床面形態(tài)的實(shí)驗(yàn)工況則相應(yīng)設(shè)置了不可滲透(亞克力棒)、可滲透2 維離散床面形態(tài)(由紗布包裹不同粒徑泥沙制備的細(xì)長條沙袋) 和沙波(由粗沙人工鋪設(shè)), 以改變床面形態(tài)組成單元的滲透率.實(shí)驗(yàn)參數(shù)和工況如表1 所示.

表1 實(shí)驗(yàn)參數(shù)和工況Table 1 Experimental parameters and runs

2 結(jié)果和分析

2.1 底床滲透率變化對(duì)潛流交換特性的影響

床面形態(tài)驅(qū)動(dòng)的潛流交換不僅與水動(dòng)力作用(如流量或流速/水深、非恒定性和近底紊動(dòng)特性) 有關(guān), 而且還涉及到底床和床面形態(tài)的物理特性(如滲透率、孔隙度).以往研究結(jié)果表明, 對(duì)于可滲透的砂質(zhì)或礫石底床, 床沙組成對(duì)潛流交換通量影響較大, 在床面附近區(qū)域上覆水流與床面形態(tài)的相互作用將會(huì)驅(qū)動(dòng)各種物質(zhì)產(chǎn)生進(jìn)出底床的泵吸交換[1,28].與以往均質(zhì)底床條件下的相關(guān)水槽實(shí)驗(yàn)相比, 本工作側(cè)重于當(dāng)定量分析床面形態(tài)和底床的物理特性存在差異時(shí), 即非均質(zhì)底床結(jié)構(gòu)條件下的潛流交換特性, 因此為對(duì)比需要也設(shè)置了相應(yīng)的均質(zhì)底床結(jié)構(gòu)實(shí)驗(yàn)工況, 包括平整和存在床面形態(tài)這2 種情況.

圖4 給出了在上覆水平均流速U= 10.6 cm/s、中沙和粗沙平整均質(zhì)底床條件下(工況1和2), 上覆水量綱一溶質(zhì)濃度C?隨時(shí)間t的變化.

圖4 中沙和粗沙平整底床條件下上覆水濃度隨時(shí)間的變化Fig.4 Variation of overlying water concentration with time under the conditions of flat beds with medium and coarse sands

由圖4 可見, 對(duì)于平整均質(zhì)底床, 在初始交換階段粗沙平整底床條件下的上覆水溶質(zhì)濃度C?曲線下降速度明顯快于中沙平整底床情形, 即工況2 在t= 50 min 左右就已達(dá)到交換平衡, 而工況1 在采樣過程結(jié)束時(shí)仍未達(dá)到交換平衡.由圖4 中C?曲線初始斜率S確定的量綱一有效擴(kuò)散系數(shù)Deff/D′m, 工況2 所得結(jié)果為Deff/D′m=1.89×104, 與工況1 所得結(jié)果(Deff/D′m= 1.90×103) 相比, 其增幅為894.7%, 表明隨底床滲透率的增大, 潛流交換通量和有效擴(kuò)散系數(shù)有明顯增大.

圖5 給出了在上覆水平均流速U=10.6 cm/s, 存在可滲透2 維離散床面形態(tài)中沙和粗沙底床條件下(w/h=3, 6, 8, 11), 上覆水量綱一溶質(zhì)濃度C?隨時(shí)間t的變化.

圖5 中沙和粗沙底床條件下(存在可滲透2 維離散床面形態(tài)) 上覆水濃度隨時(shí)間的變化Fig.5 Variation of overlying water concentration with time under the conditions of medium and coarse sand beds in the presence of permeable 2D discrete bedform

由圖5 可見, 當(dāng)存在相同滲透率的2 維離散床面形態(tài)時(shí), 底床滲透率變化對(duì)潛流交換通量的影響趨勢總體上類似于平整床面情形, 對(duì)于較高滲透率的粗沙底床, 其潛流交換通量明顯大于較低滲透率的中沙底床情形, 這說明底床滲透率的增大對(duì)潛流交換(包括湍流滲透和泵吸交換) 起促進(jìn)作用.在非均質(zhì)粗沙底床條件下的潛流交換通量隨底床滲透率的變化, 與均質(zhì)中沙底床條件下的實(shí)驗(yàn)結(jié)果相比, 各工況的增幅并不完全相同.比較分析w/h變化范圍(3 ≤w/h≤11) 內(nèi), 各相應(yīng)工況的量綱一有效擴(kuò)散系數(shù)Deff/D′m定量數(shù)據(jù).w/h= 3 時(shí)工況10 所得結(jié)果(Deff/D′m= 1.41×104) 相比工況3 (Deff/D′m= 2.39×103) 增幅為490.0%;w/h=6 時(shí)工況13 所得結(jié)果(Deff/D′m=1.82×104) 相比工況6 (Deff/D′m=3.94×103) 增幅為361.9%;w/h= 8 時(shí)工況15 所得結(jié)果(Deff/D′m= 1.96×104) 相比工況8 (Deff/D′m=2.39×103) 增幅為720.1%;w/h=11 時(shí)工況16 所得結(jié)果(Deff/D′m=2.11×104) 相比工況9 (Deff/D′m=2.12×103) 增幅為895.3%.由此可見, 非均質(zhì)粗沙底床條件下的有效擴(kuò)散系數(shù)與均質(zhì)中沙底床條件下所得結(jié)果相比, 其增幅總體上要低于平整均質(zhì)底床情形.可見, 對(duì)非均質(zhì)底床結(jié)構(gòu)而言, 在w/h變化范圍(3 ≤w/h≤11) 內(nèi), 底床和床面形態(tài)組成單元的滲透率變化均會(huì)對(duì)潛流交換通量和有效擴(kuò)散系數(shù)產(chǎn)生不同程度的影響.隨著底床和床面形態(tài)組成單元滲透率的增大, 底床表面和床面形態(tài)組成的單元內(nèi)部的湍流滲透相應(yīng)增強(qiáng), 使得潛流交換作用得以增強(qiáng).同時(shí), 床面形態(tài)驅(qū)動(dòng)的潛流交換還受到床面形態(tài)幾何特征(w/h) 的影響, 非均質(zhì)粗沙底床條件下的有效擴(kuò)散系數(shù)隨w/h的增大并不完全呈單調(diào)變化.

2.2 床面形態(tài)滲透率變化對(duì)潛流交換特性的影響

為進(jìn)一步分析當(dāng)床面形態(tài)組成單元的滲透率與底床存在差異時(shí), 這種局部尺度上的非均質(zhì)性對(duì)潛流交換通量的定量影響, 本工作利用粗沙底床條件下不同滲透率床面形態(tài)組成單元的各實(shí)驗(yàn)工況所得數(shù)據(jù), 比較分析床面形態(tài)組成單元滲透率變化對(duì)潛流交換特性的影響.

圖6 給出了在粗沙底床條件下(w/h=3, 6, 8, 11), 存在不可滲透、可滲透2 維離散床面形態(tài)和沙波各實(shí)驗(yàn)工況, 上覆水量綱一溶質(zhì)濃度C?隨時(shí)間t的變化.

圖6 粗沙底床條件下(存在不同滲透率床面形態(tài)) 上覆水濃度隨時(shí)間的變化Fig.6 Variation of overlying water concentration with time under the conditions of coarse sand beds in the presence of the bedforms with different permeabilities

由圖6 可見, 對(duì)于相同滲透率的粗沙底床隨床面形態(tài)組成單元滲透率的變化, 潛流交換通量呈較為一致的變化趨勢, 即隨著床面形態(tài)組成單元滲透率的降低, 初始交換階段上覆水溶質(zhì)濃度C?隨時(shí)間的推移下降較緩, 達(dá)到最終交換平衡的時(shí)間相對(duì)變長, 相應(yīng)潛流交換通量逐漸減小.由分析w/h變化范圍(3 ≤w/h≤11) 內(nèi)各相應(yīng)工況的量綱一有效擴(kuò)散系數(shù)Deff/D′m定量數(shù)據(jù)可知, 雖然床面形態(tài)組成單元滲透率的變化對(duì)潛流交換通量的影響程度不如底床滲透率的變化, 但對(duì)有效擴(kuò)散系數(shù)仍存在不可忽視的影響.w/h= 3 時(shí)工況10 所得結(jié)果(Deff/D′m= 1.41×104) 相比工況17 (Deff/D′m= 2.05×104) 降低了31.2%,工況24 所得結(jié)果(Deff/D′m= 7.95×103) 則降低了61.2%;w/h= 6 時(shí)工況13 所得結(jié)果(Deff/D′m=1.82×104)相比工況20(Deff/D′m=2.26×104)降低了19.5%,工況27 所得結(jié)果(Deff/D′m=1.24×104)則降低了45.1%;w/h=8 時(shí)工況15 所得結(jié)果(Deff/D′m=1.96×104)相比工況22(Deff/D′m=2.36×104)降低了16.9%, 工況29 所得結(jié)果(Deff/D′m=1.42×104)則降低了39.8%;w/h= 11 時(shí)工況16 所得結(jié)果(Deff/D′m= 2.11×104) 相比工況23(Deff/D′m= 2.19×104) 降低了3.7%, 工況30 所得結(jié)果(Deff/D′m= 1.70×104) 則降低了22.4%.上述分析表明, 床面形態(tài)組成單元滲透率相對(duì)較低的非均質(zhì)底床結(jié)構(gòu)(存在不可滲透、可滲透2 維離散床面形態(tài)) 與均質(zhì)底床結(jié)構(gòu)(存在沙波床面形態(tài)) 相比, 潛流交換趨于減弱.

2.3 水動(dòng)力影響機(jī)理分析

一般而言, 床面形態(tài)驅(qū)動(dòng)的潛流交換受到多種物理驅(qū)動(dòng)機(jī)制的共同作用, 包括分子擴(kuò)散、剪切離散、湍流滲透和泵吸交換等[29-30].對(duì)于具有一定滲透率的砂質(zhì)底床, 泵吸交換和湍流滲透是影響潛流交換特性的主要驅(qū)動(dòng)機(jī)制, 而有效擴(kuò)散系數(shù)綜合反映了多種物理機(jī)制的共同作用.根據(jù)不同實(shí)驗(yàn)工況得到的量綱一有效擴(kuò)散系數(shù), 不僅可定量分析床面形態(tài)組成單元的滲透率變化對(duì)有效擴(kuò)散系數(shù)的影響, 同時(shí)還可對(duì)比分析均質(zhì)和非均質(zhì)底床結(jié)構(gòu)條件下的潛流交換水動(dòng)力影響機(jī)理.圖7 給出了不同實(shí)驗(yàn)工況下的量綱一有效擴(kuò)散系數(shù)Deff/D′m柱狀圖, 圖中M 表示為中沙底床, C 表示為粗沙底床, 2 條虛線分別表示中沙和粗沙在平整底床條件下的Deff/D′m值.

圖7 各工況有效擴(kuò)散系數(shù)柱狀圖Fig.7 Histogram of effective diffusion coefficients for all experimental runs

由圖7 可見, 對(duì)于存在床面形態(tài)的均質(zhì)底床結(jié)構(gòu)(包括中沙底床和沙波床面形態(tài)對(duì)應(yīng)工況), 在w/h變化范圍(3 ≤w/h≤11) 內(nèi)與平整底床相比, 其有效擴(kuò)散系數(shù)均有不同程度的增大, 說明對(duì)于存在床面形態(tài)的均質(zhì)底床結(jié)構(gòu)來說, 床面形態(tài)的存在引入了附加的泵吸交換, 潛流交換受到湍流滲透和泵吸交換的共同作用, 其有效擴(kuò)散系數(shù)相比平整底床情形總體上得以增大.從分析均質(zhì)底床結(jié)構(gòu)條件下各工況的Deff/D′m定量數(shù)據(jù)可知, 與相應(yīng)的平整底床相比,存在床面形態(tài)的中沙底床典型工況(w/h=6) 的Deff/D′m增幅為107.4%, 粗沙底床典型工況(w/h= 8) 的Deff/D′m增幅為24.9%, 這也說明了泵吸交換對(duì)潛流交換的相對(duì)貢獻(xiàn)隨床沙滲透率的增大趨于減弱, 湍流滲透的相對(duì)貢獻(xiàn)趨于增強(qiáng).

對(duì)于存在床面形態(tài)的非均質(zhì)底床結(jié)構(gòu)(粗沙底床對(duì)應(yīng)工況), 隨著床面形態(tài)組成單元滲透率的降低, 各工況有效擴(kuò)散系數(shù)總體上趨于減小(見圖7).對(duì)比分析存在不可滲透、可滲透床面形態(tài)組成單元各工況的Deff/D′m定量數(shù)據(jù), 與平整底床相比其有效擴(kuò)散系數(shù)呈現(xiàn)不一致的變化趨勢, 即不可滲透床面形態(tài)對(duì)應(yīng)工況在w/h變化范圍(3 ≤w/h≤11) 內(nèi), 其有效擴(kuò)散系數(shù)均小于平整底床所得結(jié)果.而可滲透床面形態(tài)對(duì)應(yīng)工況在床面形態(tài)組成單元相對(duì)稀疏時(shí)(w/h= 8, 11), 其有效擴(kuò)散系數(shù)大于平整底床的Deff/D′m值; 當(dāng)床面形態(tài)組成單元相對(duì)密集時(shí)(w/h= 3, 6), 其有效擴(kuò)散系數(shù)小于平整底床Deff/D′m值.上述情況說明, 在非均質(zhì)底床結(jié)構(gòu)條件下, 湍流滲透和泵吸交換對(duì)潛流交換通量的相對(duì)貢獻(xiàn)不僅與底床滲透率有關(guān), 同時(shí)還受到床面形態(tài)組成單元滲透率變化的影響.具體來看, 對(duì)于存在不可滲透床面形態(tài)典型工況(w/h=8), 其有效擴(kuò)散系數(shù)Deff/D′m(1.42×104) 相比平整底床(1.89×104) 降低了24.9%, 其原因在于①不可滲透床面形態(tài)組成單元遮擋了部分底床面積, 使得上覆水/孔隙水物質(zhì)交換表面積減小, 減弱了湍流滲透; ②存在不可滲透床面形態(tài)典型工況的有效擴(kuò)散系數(shù)Deff/D′m相比沙波工況(w/h=8) 降低了39.8%, 其降幅大于平整底床情形, 由此可知泵吸交換也有一定程度的減弱.

綜合來看, 對(duì)于床面形態(tài)與底床滲透率存在差異的非均質(zhì)底床結(jié)構(gòu), 床面形態(tài)組成單元滲透率的減小對(duì)潛流交換起抑制作用, 尤其對(duì)于存在不可滲透床面形態(tài)非均質(zhì)底床結(jié)構(gòu), 與均質(zhì)底床結(jié)構(gòu)(平整或存在床面形態(tài)) 相比, 湍流滲透和泵吸交換均趨于減弱, 相應(yīng)地有效擴(kuò)散系數(shù)總體上呈不同程度的減小趨勢.

3 結(jié)束語

通過實(shí)驗(yàn)室環(huán)形水槽實(shí)驗(yàn), 采用概化的2 維離散床面形態(tài), 測量得到床面形態(tài)組成單元滲透率與潛流交換通量之間的定量關(guān)聯(lián)和變化特征, 對(duì)比分析了床面形態(tài)組成單元的滲透率變化對(duì)有效擴(kuò)散系數(shù)的影響.實(shí)驗(yàn)結(jié)果表明, 對(duì)于床面形態(tài)與底床滲透率存在差異的非均質(zhì)底床結(jié)構(gòu)來說, 床面形態(tài)組成單元的滲透率變化對(duì)潛流交換起著不可忽視的作用.在本實(shí)驗(yàn)參數(shù)變化范圍內(nèi), 床面形態(tài)的存在引入了附加的泵吸交換, 泵吸交換對(duì)潛流交換通量的相對(duì)貢獻(xiàn)與床面形態(tài)組成單元的滲透率變化密切相關(guān).隨著床面形態(tài)組成單元滲透率的降低, 底床表面和床面形態(tài)組成單元內(nèi)部的湍流滲透作用趨于減弱, 同時(shí)床面形態(tài)誘導(dǎo)的泵吸交換受到一定程度的抑制, 與平整均質(zhì)底床相比有效擴(kuò)散系數(shù)總體上呈不同程度的減小趨勢.