刘佶鑫, 王可健, 刘 瑛, 刘德富, 冯晶红, 李 毅

(湖北工业大学河湖生态修复与藻类利用湖北省重点实验室， 湖北 武汉 430068)

河流中的生物体对热变化比较敏感，只能适应一定范围内的温度，因此温度是影响有机体和生态系统功能的最重要因素之一[1]。河流潜流带中，温度可以影响微生物活性、无脊椎动物的发育、部分鱼卵的存活，具有重要的生态意义[2]。而潜流交换控制着多种栖息地因子(如温度、砂砾等)，这些因子直接控制着生物对产卵栖息地的选择。如，一些蛙和鱼将卵产在下降潜流和上升潜流的交界位置，向上的水流会带来恒定的温度和丰富的营养，向下的水流则会带来充足的溶解氧[3]。上覆水和地下水之间的压力梯度和沉积物的渗透性能是影响潜流交换的两个主要因素，水压是局部潜流交换的驱动力[4]。Menichino等[5]研究发现跌水结构可以使潜流带形成局部水力梯度，在上游产生下行流，在下游产生上行流，进而产生潜流交换，对河流温度具有一定的调节作用。Sawyer等[6-7]采用实验室水槽试验和数值模拟方法对潜流和热交换进行量化，认为河道中大型枯木结构周围的河水流动会沿河床产生压力梯度，从而推动河水与地下水的混合以及潜流带内的热传输。Malard等[8]研究发现，在较大型河流系统中，由于局部地表-地下水的相互作用，河段内部存在明显的热非均质性现象。陈孝兵等[9-10]通过构建循环式水槽装置和数值模型，研究了不同河床地形、河床渗透特性的潜流交换与地表水动力、河床沉积物非均质性之间的关系。Moore等[11]研究表明，河道中原木引起的河床复杂地形会对地表水温度的昼夜变化产生显著影响。

为恢复河道的生态功能，近年来人们仿造河流自然形态提出并实施了生态跌水恢复工程，旨在营造河流的深潭-浅滩结构，增加河道空间异质性和形态的多样性，以此达到提高生物多样性的目的。原木跌水结构是生态跌水结构的一种，国外应用和研究的较多[12-13]，国内的研究主要集中在生境恢复效果方面[14]。

目前，模拟或者试验研究更多地集中于床面形态等因素对潜流交换的影响以及潜流交换过程对地表水温的调节作用[15-17]，很少有学者研究如何设置工程结构来增强潜流交换、增强潜流带的热缓冲能力。基于以上考虑，本研究应用有限元软件COMSOL Multiphysics构建地表水-地下水二维耦合模型，研究河道内原木跌水结构对河床潜流带温度场的影响，以原木高度、间距和数量为基础变量，探究其变化对地表水与潜流带温度传递过程的影响，以期能为河流生态修复提供参考。

1 地表-地下水耦合数值模型

1.1 模型建立

COMSOL Multiphysics可实现任意多物理场直接耦合分析，与其他有限元数值分析软件相比，具有十分强大的多物理场求解功能。该软件在化学反应、弥散、流体动力学、地球科学等领域得到了广泛应用。前人应用该软件对潜流交换和热交换等开展了较多的研究[18-19]，具有较强的实用性。因此，本文利用该软件建立地表水-地下水耦合模型(图1)。模型河床长×高为3 m×0.5 m，河床表面有效水位为0.4 m，原木高度为0.2 m，放置于河床1.5 m处。模型分为地表水部分和地下水部分，地表水部分共划分7 221个网格单元，地下水部分共划分12 426个网格单元。地表水部分利用基于雷诺平均方程(N-S方程)和k-ω湍流模型计算河床表面的壁边界压力分布，并将此压力作为地下水模型的上压力边界；地下水部分是通过Darcy方程及热传导方程建立温度和流动耦合。为探究原木跌水结构对潜流带温度场的影响，通过改变原木的高度、间距和数量来建立多种工况进行对比分析。

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

1.2 初始和边界条件

1.2.1初始条件

根据河流现场实测数据得到模型初始条件：河流初始流速为0.18 m/s，水深0.4 m，地表水温度为23 ℃，河床的初始温度为21 ℃(实测平均值)，河床底质的渗透率为3.4×10-9m2。

1.2.2边界条件

如图1所示地表水部分的水力边界为：AD为进口速度边界，BC为出口边界，CD为对称边界，AB为不透水边界，E为原木结构。地下水模型水力边界为：ab为压力边界，bc、cd、ad都为不透水边界。地下水热力学边界为：ab为温度边界，温度为河床表面温度，其他边均为绝热边界。

1.3 模型验证

本文采用河流现场监测点实测数据进行模型率定。河道位于武汉市黄陂区夏家寺河(E114°28′，N31°02′)，平均河宽约2 m，水深0.2～0.5 m。选取长约50 m的顺直河段，每隔15 m放置一根横跨河道的原木并固定，共设置三级原木跌水结构(图2)。参考前人试验[20]，试验中记录温度的传感器型号为DS1922L，采样频率为10 min/次。监测河床0 m、0.1 m、0.2 m、0.3 m深度的温度数据，监测周期为6 d。温度传感器的放置位置与上文模型的1、2、3点设置位置相同(图1)。

图2 在现场设置的单级原木结构

实测河床温度与模型模拟温度的对比结果如图3所示。通过均方根误差(RMSE)和决定系数(R2)对模型模拟结果进行验证，以此判断模型的精度，R2>0.6可认为模型模拟结果较好。监测点1的RMSE值为0.139，R2值为0.930；监测点2的RMSE值为0.142，R2值为0.810；监测点3的RMSE值为0.146，R2值为0.693。各监测点均方根误差均不大于0.146，其决定系数R2均大于0.6，说明模型可靠、模拟效果较好，用来研究潜流带温度变化规律是可行的。

图3 河床温度实测值与模拟值对比

1.4 工况设置

本文以原木数量、高度和间距为变量来探讨不同结构影响下潜流带温度场的变化，模拟时间为24 h，具体工况参数见表1。

表1 不同工况参数

2 结果与分析

2.1 原木跌水结构高度对潜流带温度场的影响

2.1.1压力和流场分析

在无原木结构及工况1、2、3情况下，得到的不同原木高度影响下潜流带压力和达西流线分布(图4)。原木跌水结构设置对潜流带压力场和流场有明显影响，在原木下方潜流带局部区域，压力沿x轴方向呈现分层分布，原木上游压力表现为正值，下游压力表现为负值。由于原木附近上下游的压力差，原木上游部分水流表现为下降流，下游部分为上升流。随着原木高度的增加，河床表面压力分布也发生变化，原木上游和下游间的压力差逐渐增大，压力差极大值从31.8 Pa增加至76.3 Pa。

图4 不同原木高度影响下的潜流带压力和流场分布

2.1.2温度场分析

模拟24 h后，不同原木高度影响下的潜流带温度等值线变化如图5所示。受原木的影响，等值线在其下方形成了一个半椭圆的较高温区域，该区域温度更接近地表水温度。并随着原木高度的增加，该区域的面积在逐渐增大，潜流带内部的温度范围也发生变化。结合潜流带压力和流场变化分析，设置原木跌水结构会增大局部地下压力差，使得该区域达西流速增加，因此原木下方更易受地表水温的影响。此外在原木下游的上升流部分存在一梯形区域，该区域温差范围小，温度相对比较稳定，在纵向上无明显温度分层，因此在极端温度下可以为一些无脊椎生物提供暂时的热避难所[21]。

图5 不同原木高度下影响下的潜流带温度等值线 (单位：℃)

2.2 原木跌水结构间距对潜流带温度的影响

2.2.1压力和流场分析

在无原木跌水结构及工况4、5、6情况下，不同原木间距影响下的压力分布和地下达西流线情况如图6所示。两组原木下方的地下部分的压力沿x轴方向均有分层现象，随着原木间距的增大，第二组原木下方的压力分层现象越明显。第一组原木上游压力表现为正值，第二组原木下游压力表现为负值。随着原木间距的增大，两组原木之间地下部分的压力值在增大，且上下游的压力差最大值也在增大，其最大压力差由53.47 Pa增大至68.2 Pa。且随着原木间距的增大，两组原木之间的潜流带压力分层现象越明显，出现了多次上升流和下降流。

图6 不同间距的原木影响下的潜流带压力和流场分布

2.2.2温度场分析

在经过24 h模拟后，不同原木间距影响下的潜流带温度等值线变化如图7所示。受原木影响，原木下方会形成一个半椭圆形的高温区域，在每组原木的下游部分存在一定大小且温度相对稳定的梯形区域。随着间距的增大，两组原木之间出现了多次上升流和下降流，且对应的潜流带区域温度分层现象越明显，并引起了局部的热异质性[7]。随着原木间距的增大，第一级原木下游的缓冲区域在慢慢变大，温度缓冲效果更好。研究结果表明增加两级原木间距，可以降低地表水水温对潜流带局部的影响程度，增强了该区域在温度传递过程中的热缓冲效果。

图7 不同间距的原木影响下的潜流带温度等值线 (单位：℃)

2.3 原木跌水结构数量对潜流带温度的影响

2.3.1压力和流场分析

在无原木结构及工况7、8、9情况下，得到的不同数量的原木影响下的压力分布和地下达西流线情况如图8所示。只有最上游一组原木下方的地下部分的压力沿x方向均有明显分层现象，第一组原木上游地下部分呈现正值，最后一组原木下游压力呈现负值。当原木数量增多时，地下部分压力差极大值也在增大，其最大压力差由53.47 Pa增大至68.2 Pa。其地下流场也随着潜流带压力的分层而变化，主要的压力差分布都集中于第一组原木下方处，在其上游地下部分表现为下降流，在其下游地下部分基本表现为上升流。

图8 不同数量原木影响下的潜流带压力和流场分布

2.3.2温度场分析

在经过24 h模拟后，不同数量原木影响下的潜流带温度等值线变化如图9所示。每组原木下方均形成一个接近地表水温的半椭圆形的较高温区域。在每组原木的下游部分仍然存在一定大小的温度稳定区域。由于第一组原木下方局部压力差最大，其下方较高温区域面积最大。随着河道内原木数量的增加，该区域温度极小值减小，最后一组原木下游潜流带局部区域受地表水水温的影响程度在降低，热缓冲作用在增强。

图9 不同数量原木影响下的潜流带温度等值线 (单位：℃)

3 讨 论

3.1 原木跌水结构对潜流带温度场的影响

原木高度、间距和数量增加时，河床的压力最大值增大，原木上下游总压力差增大，且越靠近原木，压力值越大。河道原木跌水结构添加将改变河床的潜流交换，从而在河床内形成空间上的热异质性。在自然河流中，地下水的平均温度通常与地表水不同，地下水温度相对较为恒定，而地表水温度受气温影响波动较为明显。在潜流交换过程中，当水流向下流动时，孔隙水温度会继承河水的波动特征；而当向上流动时，孔隙水温度会继承深层地下水的稳定特征[7]。因此在潜流交换的上升流区域可以为一些对热变化敏感的生物提供庇护场所[1]。研究结果表明，原木下方都可形成一个接近地表水温的半椭圆形的较高温区域，而在原木下游形成类似梯形的热缓冲区，且不同工况下热缓冲区有所不同。由于热缓冲区的温度相对稳定，这会增加底栖生物和水生生境的多样性以及热避难所的可用性[21]。本文模拟了夏季较高温的地表水与河床的热交换过程，原木结构会影响河床表面压力并促进潜流带局部的下降流和上升流的形成，原木跌水结构参数的变化也会引起潜流带局部在温度传递过程中的热异质性。随着原木高度增加，原木下方区域潜流带受地表水水温的影响范围扩大；随着原木间距的增加和原木数量的增加，潜流带热缓冲作用增强。

3.2 原木跌水结构在河流生态修复中的应用

大型原木是河流生态修复中的常用材料，其目的是增强河流的自然特征，如深潭和浅滩，来增加河流生态系统中生物生存环境的异质性，减少河岸侵蚀并提供鱼类栖息地[22]。单一原木结构对河流潜流带水交换的影响可能十分有限，本研究试图用多级原木结构，并通过控制其尺寸、间距及数量来研究对它们共同作用下潜流带温度场的变化规律。发现原木跌水结构参数的变化确实能引起潜流带在热传递过程中的热异质性。因此，为了取得更有效的热缓冲能力，建议在河流生态恢复中采取增加原木的高度和数量等措施来达到热缓冲效果。通过研究，可以推断在河道内设置的原木结构对河流的影响达到一定规模，其在地表水流和潜流带中引起的热异质性足够大，甚至可以对生态环境产生直接的影响。

如何量化结构参数与热异质性的关系以及产生的热异质性对河流生态系统具体产生怎样的影响则需要进一步的研究。此外，原位河流中潜流带受影响的因素众多，原木结构在河流生态修复过程中起到的效果在很大程度上取决于周围的地下水排放或补给率，河床渗透率和原木的阻塞率等[17]。河道中原木结构的添加是恢复河流水文连通性的很有前途的工具，因此在河流生态修复过程中也应更多地考虑在河道中如何设置原木结构来增强地表水与潜流带的相互作用来增加河流生态系统的多样性[6]。

4 结 语

河道中原木跌水结构的添加会改变潜流带局部的压力分布，增加原木高度、间距和数量时，都会使河床的压力最大值增大，使得河床上下游总压力差增大，促进了潜流带局部的下降流和上升流的形成。

在24 h内的温度传递过程中，原木跌水结构下游的上升流区域会存在一梯形的温度稳定区域，该区域的温差范围较小，不易受地表水温的影响，可能会为一些底栖生物提供热避难所。

原木跌水结构加快了地表水与河床的热交换，使原木下方区域形成接近地表水水温的半椭圆形高温区域。原木高度、间距及数量的不同均会影响潜流带在热传递过程中的温度场的变化，引起河床局部热异质性和热缓冲区的变化。随着原木高度增加，原木下方区域潜流带受地表水水温的影响范围扩大；随着原木间距的增加，原木间的热缓冲作用增强；随着原木数量的增加，最后一组原木下游区域热缓冲作用增强。

本文所构建模型为二维概化模型，对河床在空间上的异质性的影响没有进行分析。构建三维模型，考虑河床空间异质性及原木参数的优化需要进一步研究。