劉佶鑫, 王可健, 劉 瑛, 劉德富, 馮晶紅, 李 毅

(湖北工業(yè)大學(xué)河湖生態(tài)修復(fù)與藻類利用湖北省重點(diǎn)實(shí)驗(yàn)室， 湖北 武漢 430068)

河流中的生物體對(duì)熱變化比較敏感，只能適應(yīng)一定范圍內(nèi)的溫度，因此溫度是影響有機(jī)體和生態(tài)系統(tǒng)功能的最重要因素之一[1]。河流潛流帶中，溫度可以影響微生物活性、無脊椎動(dòng)物的發(fā)育、部分魚卵的存活，具有重要的生態(tài)意義[2]。而潛流交換控制著多種棲息地因子(如溫度、砂礫等)，這些因子直接控制著生物對(duì)產(chǎn)卵棲息地的選擇。如，一些蛙和魚將卵產(chǎn)在下降潛流和上升潛流的交界位置，向上的水流會(huì)帶來恒定的溫度和豐富的營養(yǎng)，向下的水流則會(huì)帶來充足的溶解氧[3]。上覆水和地下水之間的壓力梯度和沉積物的滲透性能是影響潛流交換的兩個(gè)主要因素，水壓是局部潛流交換的驅(qū)動(dòng)力[4]。Menichino等[5]研究發(fā)現(xiàn)跌水結(jié)構(gòu)可以使?jié)摿鲙纬删植克μ荻?，在上游產(chǎn)生下行流，在下游產(chǎn)生上行流，進(jìn)而產(chǎn)生潛流交換，對(duì)河流溫度具有一定的調(diào)節(jié)作用。Sawyer等[6-7]采用實(shí)驗(yàn)室水槽試驗(yàn)和數(shù)值模擬方法對(duì)潛流和熱交換進(jìn)行量化，認(rèn)為河道中大型枯木結(jié)構(gòu)周圍的河水流動(dòng)會(huì)沿河床產(chǎn)生壓力梯度，從而推動(dòng)河水與地下水的混合以及潛流帶內(nèi)的熱傳輸。Malard等[8]研究發(fā)現(xiàn)，在較大型河流系統(tǒng)中，由于局部地表-地下水的相互作用，河段內(nèi)部存在明顯的熱非均質(zhì)性現(xiàn)象。陳孝兵等[9-10]通過構(gòu)建循環(huán)式水槽裝置和數(shù)值模型，研究了不同河床地形、河床滲透特性的潛流交換與地表水動(dòng)力、河床沉積物非均質(zhì)性之間的關(guān)系。Moore等[11]研究表明，河道中原木引起的河床復(fù)雜地形會(huì)對(duì)地表水溫度的晝夜變化產(chǎn)生顯著影響。

為恢復(fù)河道的生態(tài)功能，近年來人們仿造河流自然形態(tài)提出并實(shí)施了生態(tài)跌水恢復(fù)工程，旨在營造河流的深潭-淺灘結(jié)構(gòu)，增加河道空間異質(zhì)性和形態(tài)的多樣性，以此達(dá)到提高生物多樣性的目的。原木跌水結(jié)構(gòu)是生態(tài)跌水結(jié)構(gòu)的一種，國外應(yīng)用和研究的較多[12-13]，國內(nèi)的研究主要集中在生境恢復(fù)效果方面[14]。

目前，模擬或者試驗(yàn)研究更多地集中于床面形態(tài)等因素對(duì)潛流交換的影響以及潛流交換過程對(duì)地表水溫的調(diào)節(jié)作用[15-17]，很少有學(xué)者研究如何設(shè)置工程結(jié)構(gòu)來增強(qiáng)潛流交換、增強(qiáng)潛流帶的熱緩沖能力?；谝陨峡紤]，本研究應(yīng)用有限元軟件COMSOL Multiphysics構(gòu)建地表水-地下水二維耦合模型，研究河道內(nèi)原木跌水結(jié)構(gòu)對(duì)河床潛流帶溫度場的影響，以原木高度、間距和數(shù)量為基礎(chǔ)變量，探究其變化對(duì)地表水與潛流帶溫度傳遞過程的影響，以期能為河流生態(tài)修復(fù)提供參考。

1 地表-地下水耦合數(shù)值模型

1.1 模型建立

COMSOL Multiphysics可實(shí)現(xiàn)任意多物理場直接耦合分析，與其他有限元數(shù)值分析軟件相比，具有十分強(qiáng)大的多物理場求解功能。該軟件在化學(xué)反應(yīng)、彌散、流體動(dòng)力學(xué)、地球科學(xué)等領(lǐng)域得到了廣泛應(yīng)用。前人應(yīng)用該軟件對(duì)潛流交換和熱交換等開展了較多的研究[18-19]，具有較強(qiáng)的實(shí)用性。因此，本文利用該軟件建立地表水-地下水耦合模型(圖1)。模型河床長×高為3 m×0.5 m，河床表面有效水位為0.4 m，原木高度為0.2 m，放置于河床1.5 m處。模型分為地表水部分和地下水部分，地表水部分共劃分7 221個(gè)網(wǎng)格單元，地下水部分共劃分12 426個(gè)網(wǎng)格單元。地表水部分利用基于雷諾平均方程(N-S方程)和k-ω湍流模型計(jì)算河床表面的壁邊界壓力分布，并將此壓力作為地下水模型的上壓力邊界；地下水部分是通過Darcy方程及熱傳導(dǎo)方程建立溫度和流動(dòng)耦合。為探究原木跌水結(jié)構(gòu)對(duì)潛流帶溫度場的影響，通過改變?cè)镜母叨?、間距和數(shù)量來建立多種工況進(jìn)行對(duì)比分析。

圖1 地表-地下水耦合模型示意圖

1.2 初始和邊界條件

1.2.1初始條件

根據(jù)河流現(xiàn)場實(shí)測(cè)數(shù)據(jù)得到模型初始條件：河流初始流速為0.18 m/s，水深0.4 m，地表水溫度為23 ℃，河床的初始溫度為21 ℃(實(shí)測(cè)平均值)，河床底質(zhì)的滲透率為3.4×10-9m2。

1.2.2邊界條件

如圖1所示地表水部分的水力邊界為：AD為進(jìn)口速度邊界，BC為出口邊界，CD為對(duì)稱邊界，AB為不透水邊界，E為原木結(jié)構(gòu)。地下水模型水力邊界為：ab為壓力邊界，bc、cd、ad都為不透水邊界。地下水熱力學(xué)邊界為：ab為溫度邊界，溫度為河床表面溫度，其他邊均為絕熱邊界。

1.3 模型驗(yàn)證

本文采用河流現(xiàn)場監(jiān)測(cè)點(diǎn)實(shí)測(cè)數(shù)據(jù)進(jìn)行模型率定。河道位于武漢市黃陂區(qū)夏家寺河(E114°28′，N31°02′)，平均河寬約2 m，水深0.2～0.5 m。選取長約50 m的順直河段，每隔15 m放置一根橫跨河道的原木并固定，共設(shè)置三級(jí)原木跌水結(jié)構(gòu)(圖2)。參考前人試驗(yàn)[20]，試驗(yàn)中記錄溫度的傳感器型號(hào)為DS1922L，采樣頻率為10 min/次。監(jiān)測(cè)河床0 m、0.1 m、0.2 m、0.3 m深度的溫度數(shù)據(jù)，監(jiān)測(cè)周期為6 d。溫度傳感器的放置位置與上文模型的1、2、3點(diǎn)設(shè)置位置相同(圖1)。

圖2 在現(xiàn)場設(shè)置的單級(jí)原木結(jié)構(gòu)

實(shí)測(cè)河床溫度與模型模擬溫度的對(duì)比結(jié)果如圖3所示。通過均方根誤差(RMSE)和決定系數(shù)(R2)對(duì)模型模擬結(jié)果進(jìn)行驗(yàn)證，以此判斷模型的精度，R2>0.6可認(rèn)為模型模擬結(jié)果較好。監(jiān)測(cè)點(diǎn)1的RMSE值為0.139，R2值為0.930；監(jiān)測(cè)點(diǎn)2的RMSE值為0.142，R2值為0.810；監(jiān)測(cè)點(diǎn)3的RMSE值為0.146，R2值為0.693。各監(jiān)測(cè)點(diǎn)均方根誤差均不大于0.146，其決定系數(shù)R2均大于0.6，說明模型可靠、模擬效果較好，用來研究潛流帶溫度變化規(guī)律是可行的。

圖3 河床溫度實(shí)測(cè)值與模擬值對(duì)比

1.4 工況設(shè)置

本文以原木數(shù)量、高度和間距為變量來探討不同結(jié)構(gòu)影響下潛流帶溫度場的變化，模擬時(shí)間為24 h，具體工況參數(shù)見表1。

表1 不同工況參數(shù)

2 結(jié)果與分析

2.1 原木跌水結(jié)構(gòu)高度對(duì)潛流帶溫度場的影響

2.1.1壓力和流場分析

在無原木結(jié)構(gòu)及工況1、2、3情況下，得到的不同原木高度影響下潛流帶壓力和達(dá)西流線分布(圖4)。原木跌水結(jié)構(gòu)設(shè)置對(duì)潛流帶壓力場和流場有明顯影響，在原木下方潛流帶局部區(qū)域，壓力沿x軸方向呈現(xiàn)分層分布，原木上游壓力表現(xiàn)為正值，下游壓力表現(xiàn)為負(fù)值。由于原木附近上下游的壓力差，原木上游部分水流表現(xiàn)為下降流，下游部分為上升流。隨著原木高度的增加，河床表面壓力分布也發(fā)生變化，原木上游和下游間的壓力差逐漸增大，壓力差極大值從31.8 Pa增加至76.3 Pa。

圖4 不同原木高度影響下的潛流帶壓力和流場分布

2.1.2溫度場分析

模擬24 h后，不同原木高度影響下的潛流帶溫度等值線變化如圖5所示。受原木的影響，等值線在其下方形成了一個(gè)半橢圓的較高溫區(qū)域，該區(qū)域溫度更接近地表水溫度。并隨著原木高度的增加，該區(qū)域的面積在逐漸增大，潛流帶內(nèi)部的溫度范圍也發(fā)生變化。結(jié)合潛流帶壓力和流場變化分析，設(shè)置原木跌水結(jié)構(gòu)會(huì)增大局部地下壓力差，使得該區(qū)域達(dá)西流速增加，因此原木下方更易受地表水溫的影響。此外在原木下游的上升流部分存在一梯形區(qū)域，該區(qū)域溫差范圍小，溫度相對(duì)比較穩(wěn)定，在縱向上無明顯溫度分層，因此在極端溫度下可以為一些無脊椎生物提供暫時(shí)的熱避難所[21]。

圖5 不同原木高度下影響下的潛流帶溫度等值線 (單位：℃)

2.2 原木跌水結(jié)構(gòu)間距對(duì)潛流帶溫度的影響

2.2.1壓力和流場分析

在無原木跌水結(jié)構(gòu)及工況4、5、6情況下，不同原木間距影響下的壓力分布和地下達(dá)西流線情況如圖6所示。兩組原木下方的地下部分的壓力沿x軸方向均有分層現(xiàn)象，隨著原木間距的增大，第二組原木下方的壓力分層現(xiàn)象越明顯。第一組原木上游壓力表現(xiàn)為正值，第二組原木下游壓力表現(xiàn)為負(fù)值。隨著原木間距的增大，兩組原木之間地下部分的壓力值在增大，且上下游的壓力差最大值也在增大，其最大壓力差由53.47 Pa增大至68.2 Pa。且隨著原木間距的增大，兩組原木之間的潛流帶壓力分層現(xiàn)象越明顯，出現(xiàn)了多次上升流和下降流。

圖6 不同間距的原木影響下的潛流帶壓力和流場分布

2.2.2溫度場分析

在經(jīng)過24 h模擬后，不同原木間距影響下的潛流帶溫度等值線變化如圖7所示。受原木影響，原木下方會(huì)形成一個(gè)半橢圓形的高溫區(qū)域，在每組原木的下游部分存在一定大小且溫度相對(duì)穩(wěn)定的梯形區(qū)域。隨著間距的增大，兩組原木之間出現(xiàn)了多次上升流和下降流，且對(duì)應(yīng)的潛流帶區(qū)域溫度分層現(xiàn)象越明顯，并引起了局部的熱異質(zhì)性[7]。隨著原木間距的增大，第一級(jí)原木下游的緩沖區(qū)域在慢慢變大，溫度緩沖效果更好。研究結(jié)果表明增加兩級(jí)原木間距，可以降低地表水水溫對(duì)潛流帶局部的影響程度，增強(qiáng)了該區(qū)域在溫度傳遞過程中的熱緩沖效果。

圖7 不同間距的原木影響下的潛流帶溫度等值線 (單位：℃)

2.3 原木跌水結(jié)構(gòu)數(shù)量對(duì)潛流帶溫度的影響

2.3.1壓力和流場分析

在無原木結(jié)構(gòu)及工況7、8、9情況下，得到的不同數(shù)量的原木影響下的壓力分布和地下達(dá)西流線情況如圖8所示。只有最上游一組原木下方的地下部分的壓力沿x方向均有明顯分層現(xiàn)象，第一組原木上游地下部分呈現(xiàn)正值，最后一組原木下游壓力呈現(xiàn)負(fù)值。當(dāng)原木數(shù)量增多時(shí)，地下部分壓力差極大值也在增大，其最大壓力差由53.47 Pa增大至68.2 Pa。其地下流場也隨著潛流帶壓力的分層而變化，主要的壓力差分布都集中于第一組原木下方處，在其上游地下部分表現(xiàn)為下降流，在其下游地下部分基本表現(xiàn)為上升流。

圖8 不同數(shù)量原木影響下的潛流帶壓力和流場分布

2.3.2溫度場分析

在經(jīng)過24 h模擬后，不同數(shù)量原木影響下的潛流帶溫度等值線變化如圖9所示。每組原木下方均形成一個(gè)接近地表水溫的半橢圓形的較高溫區(qū)域。在每組原木的下游部分仍然存在一定大小的溫度穩(wěn)定區(qū)域。由于第一組原木下方局部壓力差最大，其下方較高溫區(qū)域面積最大。隨著河道內(nèi)原木數(shù)量的增加，該區(qū)域溫度極小值減小，最后一組原木下游潛流帶局部區(qū)域受地表水水溫的影響程度在降低，熱緩沖作用在增強(qiáng)。

圖9 不同數(shù)量原木影響下的潛流帶溫度等值線 (單位：℃)

3 討 論

3.1 原木跌水結(jié)構(gòu)對(duì)潛流帶溫度場的影響

原木高度、間距和數(shù)量增加時(shí)，河床的壓力最大值增大，原木上下游總壓力差增大，且越靠近原木，壓力值越大。河道原木跌水結(jié)構(gòu)添加將改變河床的潛流交換，從而在河床內(nèi)形成空間上的熱異質(zhì)性。在自然河流中，地下水的平均溫度通常與地表水不同，地下水溫度相對(duì)較為恒定，而地表水溫度受氣溫影響波動(dòng)較為明顯。在潛流交換過程中，當(dāng)水流向下流動(dòng)時(shí)，孔隙水溫度會(huì)繼承河水的波動(dòng)特征；而當(dāng)向上流動(dòng)時(shí)，孔隙水溫度會(huì)繼承深層地下水的穩(wěn)定特征[7]。因此在潛流交換的上升流區(qū)域可以為一些對(duì)熱變化敏感的生物提供庇護(hù)場所[1]。研究結(jié)果表明，原木下方都可形成一個(gè)接近地表水溫的半橢圓形的較高溫區(qū)域，而在原木下游形成類似梯形的熱緩沖區(qū)，且不同工況下熱緩沖區(qū)有所不同。由于熱緩沖區(qū)的溫度相對(duì)穩(wěn)定，這會(huì)增加底棲生物和水生生境的多樣性以及熱避難所的可用性[21]。本文模擬了夏季較高溫的地表水與河床的熱交換過程，原木結(jié)構(gòu)會(huì)影響河床表面壓力并促進(jìn)潛流帶局部的下降流和上升流的形成，原木跌水結(jié)構(gòu)參數(shù)的變化也會(huì)引起潛流帶局部在溫度傳遞過程中的熱異質(zhì)性。隨著原木高度增加，原木下方區(qū)域潛流帶受地表水水溫的影響范圍擴(kuò)大；隨著原木間距的增加和原木數(shù)量的增加，潛流帶熱緩沖作用增強(qiáng)。

3.2 原木跌水結(jié)構(gòu)在河流生態(tài)修復(fù)中的應(yīng)用

大型原木是河流生態(tài)修復(fù)中的常用材料，其目的是增強(qiáng)河流的自然特征，如深潭和淺灘，來增加河流生態(tài)系統(tǒng)中生物生存環(huán)境的異質(zhì)性，減少河岸侵蝕并提供魚類棲息地[22]。單一原木結(jié)構(gòu)對(duì)河流潛流帶水交換的影響可能十分有限，本研究試圖用多級(jí)原木結(jié)構(gòu)，并通過控制其尺寸、間距及數(shù)量來研究對(duì)它們共同作用下潛流帶溫度場的變化規(guī)律。發(fā)現(xiàn)原木跌水結(jié)構(gòu)參數(shù)的變化確實(shí)能引起潛流帶在熱傳遞過程中的熱異質(zhì)性。因此，為了取得更有效的熱緩沖能力，建議在河流生態(tài)恢復(fù)中采取增加原木的高度和數(shù)量等措施來達(dá)到熱緩沖效果。通過研究，可以推斷在河道內(nèi)設(shè)置的原木結(jié)構(gòu)對(duì)河流的影響達(dá)到一定規(guī)模，其在地表水流和潛流帶中引起的熱異質(zhì)性足夠大，甚至可以對(duì)生態(tài)環(huán)境產(chǎn)生直接的影響。

如何量化結(jié)構(gòu)參數(shù)與熱異質(zhì)性的關(guān)系以及產(chǎn)生的熱異質(zhì)性對(duì)河流生態(tài)系統(tǒng)具體產(chǎn)生怎樣的影響則需要進(jìn)一步的研究。此外，原位河流中潛流帶受影響的因素眾多，原木結(jié)構(gòu)在河流生態(tài)修復(fù)過程中起到的效果在很大程度上取決于周圍的地下水排放或補(bǔ)給率，河床滲透率和原木的阻塞率等[17]。河道中原木結(jié)構(gòu)的添加是恢復(fù)河流水文連通性的很有前途的工具，因此在河流生態(tài)修復(fù)過程中也應(yīng)更多地考慮在河道中如何設(shè)置原木結(jié)構(gòu)來增強(qiáng)地表水與潛流帶的相互作用來增加河流生態(tài)系統(tǒng)的多樣性[6]。

4 結(jié) 語

河道中原木跌水結(jié)構(gòu)的添加會(huì)改變潛流帶局部的壓力分布，增加原木高度、間距和數(shù)量時(shí)，都會(huì)使河床的壓力最大值增大，使得河床上下游總壓力差增大，促進(jìn)了潛流帶局部的下降流和上升流的形成。

在24 h內(nèi)的溫度傳遞過程中，原木跌水結(jié)構(gòu)下游的上升流區(qū)域會(huì)存在一梯形的溫度穩(wěn)定區(qū)域，該區(qū)域的溫差范圍較小，不易受地表水溫的影響，可能會(huì)為一些底棲生物提供熱避難所。

原木跌水結(jié)構(gòu)加快了地表水與河床的熱交換，使原木下方區(qū)域形成接近地表水水溫的半橢圓形高溫區(qū)域。原木高度、間距及數(shù)量的不同均會(huì)影響潛流帶在熱傳遞過程中的溫度場的變化，引起河床局部熱異質(zhì)性和熱緩沖區(qū)的變化。隨著原木高度增加，原木下方區(qū)域潛流帶受地表水水溫的影響范圍擴(kuò)大；隨著原木間距的增加，原木間的熱緩沖作用增強(qiáng)；隨著原木數(shù)量的增加，最后一組原木下游區(qū)域熱緩沖作用增強(qiáng)。

本文所構(gòu)建模型為二維概化模型，對(duì)河床在空間上的異質(zhì)性的影響沒有進(jìn)行分析。構(gòu)建三維模型，考慮河床空間異質(zhì)性及原木參數(shù)的優(yōu)化需要進(jìn)一步研究。