岩溶水柜集流坡面的径流模式和调控建议

姜光辉,郭 芳

(中国地质科学院岩溶地质研究所自然资源部/广西岩溶动力学重点实验室,广西 桂林 541004)

南方岩溶大石山区采取修建水柜来调蓄山坡径流。大石山的山坡径流持续时间短,水柜普遍存在季节性储水量不足的问题[1],加上气候湿热,极易发生细菌和浮游生物污染[2]。目前水柜的研究局限在工程设计建造上[3],对于集流场关注不够。山坡的蓄水及产流机制不仅在水柜集流管理上发挥重要作用,同时也是揭示石漠化过程和构建流域水文模型的基础。国内多通过开展水土流失观测分析坡面产流影响因子[4],而对于包括表层岩溶带和包气带的山坡整体水文过程缺少探索[5]。国际上岩溶地区水文建模也由于缺少试验数据,而多沿用经典产流模式[6]。

桂林丫吉试验场以南方裸露型岩溶水文调蓄机制研究为目标,构建了模拟岩溶泉流量过程的丫吉模型[7]。本研究以丫吉试验场为例,获取一处表层岩溶泉的水文和水化学高密度监测数据,借此探讨峰丛洼地裸露型岩溶的山坡径流模式,结合水柜工程有针对性地提出增加集流效率的调控措施。

1 研究区概况

桂林丫吉试验场创建于1986年。场地包含1个完整的岩溶泉流域。位于补给区的1号洼地底部标高271 m,高出作为地下水排泄出口的S31岩溶泉120 m；其下部地下水水位176 m,低于洼地底部95 m (图1)。洼地高悬在地下水位之上,产流经由地表、表层岩溶带、包气带等多种路径汇集至落水洞,然后进入岩溶管道,最终排泄于岩溶泉。

山坡汇流过程中,水流从岩溶洼地的坡顶向坡脚运动,水量不断漏失,并反复发生流入和流出坡面的现象,形成串流和回归流。集中流出时构成了泉水,称为表层岩溶泉。1号洼地有多处水流复现点,其流量为0～40 L/s,维持时间有几天至几个月。在所有表层岩溶泉中S55维持时间最长。

图1 研究区示意

表层岩溶泉S55出露于1号洼地底部,泉水补给区位于泉水东边山坡。山坡垭口与泉口的高差达到80 m,山坡长为140 m,补给区坡度为25°～35°。坡面基岩裸露,土壤分布在溶蚀沟、槽、坑、洞等地表负形态之中。在洼地底部则有土层覆盖(图1(b))。土壤为棕色石灰土,容重0.9～1.4 g/cm3。洼地底部为乔木林,往上逐渐变为灌木林,分水岭地带分布稀疏低矮的灌木。袁道先等[7]构建丫吉模型时曾将S55泉作为观测点,用其水位变化校正模型中代表表层岩溶带水箱的参数(丫吉模型是线性水箱模型)。

2 研究方法

本研究选择在有工作基础的丫吉试验场地开展,水文地质条件明确,而且有其他的研究成果作为参考。S55作为表层岩溶泉观测站,泉口利用澳大利亚GREENSPAN公司生产的 CTDP300仪器自动记录水温、电导率和pH值、水压和降雨量,记录时间间隔为15 min,测量精度分别为0.01 ℃、0.01 μS/cm、0.01个pH单位、0.01 m和0.5 mm。自动记录时间为2004年3月至2005年5月,期间共有52场降雨。后改用HOBO水位计监测水位,获得了2017—2018年观测数据,作为对照。

在1996年丫吉模型的基础上[7],经过野外观察完善表层岩溶泉概念模型。研究还参考了山坡截流沟、洞穴滴水观测站等多种观测设施和试验手段取得的结果,综合不同层次水流的水文及水化学动态,概括成为山坡径流模式。采用历时曲线法进行水情分析。利用PEAKFIT软件的正态函数进行电导率分布曲线的分峰拟合,分析径流成分。根据水- 岩- 气(CO2)相互作用原理建立起来的水化学方法在径流分解中得到完善和应用。

S55表层岩溶泉与溶蚀裂隙连通，暴雨时岩溶管道的排泄能力小于来水量,径流溢流成泉。剧烈的水文动态和复杂的排泄方式使得泉水流量计算困难,而水化学数据完整且分辨率高。虽然在坑内修建了堰(图2),但在绝大部分的时间里,水流从堰的底部漏掉,或者淹没堰顶。本研究采用水位来表示径流的水文动态,主要运用水化学方法指示径流来源、组分和过程。

图2 S55表层岩溶泉的水文动态模式

3 结果与分析

S55表层岩溶泉作为岩溶洼地中一处山坡的集总水源,泉水的水文和水化学动态能够反映岩溶石质山坡产汇流过程。

3.1 山坡产水量动态

在1个水文年内S55表层岩溶泉的产水量变化可分为3个阶段:3—7月属于雨季,泉水波动频繁；8—10月进入旱季,泉水进入持续稳定的衰减阶段；11—3月为从旱季到雨季的过渡期,径流处于断断续续的状态(图3)。泉水有6个月时间水位低于全年整体水平(图4)。

图3 S55表层岩溶泉水文和水化学的动态

图4 S55泉水位历时曲线

泉管道水位变化范围为0至1.51 m,但管道处于不同水位的时间并不平均。将发生概率小于5%作为洪水水位,用W5表示,对于S55泉,W5=0.41 m。水位的中位数W50=0.11 m,要比0和1.51 m的中间值0.76 m低得多,同时也低于洪水位。在这里定义低于W50的水位为低水位,此时水流缓慢,水量较小。当水位接近零时(H≤0.02 m),水流为残留水,定义此时的状态为枯水位。在这里水位大于枯水位发生的概率是87%,W87=0.02 m,比通常采用的枯水概率标准95%要小。泉水连续处于枯水状态长达40多天。在水位超过0.45 m的时候,水位随流量增大呈现急剧上升,反映径流从明渠流到有压管道流状态的转变。在1个水文年中S55泉处于有压管道流的概率是1.5%,表示为W1.5=0.45 m。当泉水处于最高水位1.51 m时,泉水溢出地表,形成地面径流(图2)。

水位处于不同状态的概率凸显了集流场径流的水文模式(图4,2004—2005年数据,2004- 10曲线因断流未显示)。在1个月内S55泉水位可以经历从有压管道流水位(W1.5)至枯水位(W87)的变化。在1个水文年中水资源的分配高低悬殊。8月有64%的时间处于低水位之下。9月、10月和12月,这3个月泉水的水位完全处于低水位之下,其中10月的水位完全处于枯水位。泉水的水文动态说明山坡存在较长的干旱期,需要采取人工调控。

3.2 水化学对径流的指示作用

借助于水温、pH和电导率指标来判断径流组成和来源的变化。水温总体走势曲线的峰值出现在9月(图3),与气温之间存在1个月时间的延迟。峰值开始形成的时间正好是泉水水位持续降低的时间。降雨补给引起水温波动,显示坡面快速的产汇流过程。降雨补给的水化学效应表现为电导率、pH与水温的同步变化。在季节尺度上电导率和pH代表岩溶作用的强度,它们的变化受地表层生物活动强度和土壤CO2浓度的控制[7]。

水化学指标对坡面径流来源具有指示作用。由于不同来源的径流滞留时间以及水化学特征的差异,使得水化学曲线出现明显的波动信号,对波动信号的解译可以推断出其中的成分或来源。以电导率指标为例,通过对全年电导率分布曲线的分峰拟合[8],分离出包含的6类径流组分(图5)。电导率从低到高表示径流在山坡上停留的时间逐渐增加。处于中间位置的峰体为主峰,代表泉水处于低水位时期的水化学状态,所占的比例最大。而在其左边的2个峰体代表了大雨或暴雨的影响。处于峰体右边的3个峰体则是处于枯水位附近的水化学状态,主峰左右5个峰的占比较低。

图5 S55表层岩溶泉电导率分布曲线的分峰拟合结果

3.3 山坡的产流和汇流

本次试验观测到S55泉在10月断流,至11月复涌。经历近1个月的断流后集流场的水分亏缺接近其最大值。根据径流复出的降雨量条件判断集流场的水分亏缺量为38 mm,表示山坡覆被层都达到饱和或者田间持水量所需要的降雨量。这里存在一个假设的条件是整个集流场满足饱和产流的模式。集流场各种下垫面的产流有先后顺序,所以观测到的结果比理论上的最大蓄水容量偏小。该观测结果与中国南方其他峰丛洼地地区可对比[9]。

大石山地区集流面裸岩大片出露的地方蓄水容量最小,径流也最先出现,并且所产生的径流的水化学与其他径流类型差别显著,是所有径流中最活跃的组分。2004至2005年观测到的降雨有52场,时间跨度超过1个水文年(图3)。期间水化学的变化有季节和场雨2个周期,季节变化对水化学初始状态影响较大,这使得利用水化学指标判断径流来源不宜采用1个固定的阈值,可分为枯季和雨季分别讨论。

枯季山坡集流场土壤干燥,裸岩分布区相比其他下垫面类型更易产流。选取2004年12月的1次降雨事件为例说明枯季土壤初始蓄水量小的情况下集流场的产流模式(图6)。12月水位一直处于低水位,此次降雨量为12.5 mm,降雨历时17 h,属于降雨强度小的中雨。降雨引起水位有小幅度的升高,pH值下降,温度出现峰值。电导率变化范围430～480 μS/cm,根据电导率分峰拟合的结果(图5),判断此时水流滞留时间较长。电导率在开始有升高的现象,反映“滞留水”的干扰,随后电导率降低,并且因为“老水”被新水不断推出,电导率会不断下降。水位的变化与水化学的变化同步。在水分含量少且降雨量和强度都低的情况下,坡面流、壤中流很难发生,泉水的变化由石面流引起。

选取2004年4月的1次降雨事件说明雨季土壤水分含量高且降雨量大条件下的产流机制(图6)。此次降雨前4 d经历1次暴雨过程,集流场满足上述条件。此次降雨过程在2 h 45 min内降雨量达23.5 mm。S55泉的初始水位高于低水位,降雨后水位升高迅速,同时电导率、pH值和水温快速下降,表明石面流出现。电导率变化范围290～340 μS/cm,根据电导率分峰拟合的结果(图5),此时水流为滞留时间较短的新水。电导率和水温的变化基本一致,但是pH值在出现短暂的降低后出现一个峰值。原因是石面流中含有的少量的侵蚀性CO2很快被岩溶作用消耗掉,造成pH的升高；从水化学变化判断石面流影响持续了约20 h。此后电导率恢复的速度减慢,而温度曲线和水文曲线出现了一个缓慢的峰值,pH值保持稳定,表示壤中流逐渐占主导地位。在石面流和壤中流交替的时候,水位变化并没有出现显著的波动。主要原因是降雨过程中形成的所有径流都要经过岩溶管道排泄,在汇流的过程中管道形成了一定的缓冲作用。

图6 S55泉水化学曲线显示径流来源变化

综合分析坡地径流的水化学和水文曲线,可以推断降雨时石面流快速形成,且其水化学的特点是具有较低的电导率和pH值,水温与其他类型径流有较明显的差异。水文曲线的峰值与水化学曲线的变化在时间上对应,表示洪峰是由石面流组成的。在水文曲线的衰减阶段,径流主要是由壤中流构成。

山坡集流场在降雨后主要产生2种来源的径流。第1种是石面流,它能够顺利进入表层岩溶带,并引起水化学显著的变化,并且其产生条件容易满足；第2种为壤中流,多发生在土壤与岩石的界面附近,它产生的条件较高,但能维持较长的时间。2种径流产生后都进入表层岩溶带,汇集成表层岩溶带径流。

3.4 山坡集流场径流模式的形成和调控方法

图7 2004—2005时段与2017—2018时段水位曲线的对比

2004和2005年桂林市降雨量分别为1 843和1 877 mm,接近多年平均降雨量1 885 mm(1951—2014年)。2004年9月和10月的降雨量偏低,这直接导致同期泉水断流。2004年泉水断流并非偶然现象,与桂林发生秋旱有关。有关研究显示桂林20世纪90年代秋旱发生的年份占到了40%[15]。2004年华南多个地区出现了秋旱,本次研究获得的数据是受到了干旱的影响。2017—2018年S55泉处于枯水位的时长达到84 d,是2004—2005年的2倍(图7),这主要是因为该时段秋冬季降雨量偏低。2004—2005时段的观测结果受到干旱影响,但还不属于最严重的情况。

大石山坡地集流场的径流模式的形成与前期山坡水分状态、降雨类型以及下垫面物理结构有关。地质、地貌、植被等条件决定了集流场结构多样,可以采取物探的手段对此进行详细的勘探[10]。从产汇流功能上坡地集流场的结构概括为3种基本类型的组合。第1种类型称为裸岩区,岩石的溶蚀程度特别高,地表透水性非常好,主要分布在山脊和陡坡；第2种类型是岩石和土壤的组合,称之为半覆盖区,主要分布在山坡中部至坡脚；第3种类型土壤在上层,岩石在下层,称为覆盖区,分布在洼地(图1(b))。

裸岩区超渗产流的门槛容易满足,但径流形成后很快形成下降流转移至溶蚀裂缝中形成深层渗流,只有小部分能够顺着坡面往下运动。半覆盖区坡地的蓄水容量和入渗率增大,超渗产流减少。而土壤厚度和体积一般较小,在下伏基岩的顶托下易形成饱和径流，石面流汇集到土壤,进一步促使土壤达到饱和。在覆盖区土壤的厚度超过1 m,产流门槛进一步增加。顺坡而下的径流汇集在岩土界面附近,形成滞水或岩溶管道流[11]。

发挥集流面的水源潜力、增加水柜的储水量关键在于提高山坡产水效率并且延长径流在表层岩溶带中的停留时间。根据大石山坡地的径流模式,提出以下调控措施：① 裸岩区调控目标为减少径流深层渗漏,采取人工防渗和表面固化并举的措施。② 半覆盖区调控的目标为保持土壤水分,避免水分快速消耗。参考已有研究提出可用的措施为恢复植被、增加石面附着生物量和凋落物层,提高表面缓冲能力[12]；改变土壤的持水性,提高土壤的蓄水能力[13]。③ 覆盖区调控以截留更多径流为目标，通过植物根系的发育改善含水介质的通透性[14],延长含水体内水流的滞留时间。3种措施同时开展能够提高山坡产水率和水资源利用效率。

4 结 论

本研究对位于桂林峰丛洼地发育在纯灰岩上的表层岩溶泉进行了2个水文年2个时段的自动化观测，揭示了山坡水文动态和不同降雨条件下的坡地产汇流过程，为大石山缺水地区水利开发提供支撑。主要结论如下：

(1) 大石山坡地以表层岩溶泉的形式向水柜提供水源。水源动态极不稳定,存在累计长达半年的低水位期,其连续时长超过3个月,枯水期连续时长超过40 d,这可作为水柜容量设计的依据 。

(2) 大石山坡地能够提供的水流是不同界面和介质中由降雨转化形成的多种径流成分的组合。在岩溶高度发育的地质条件下,不论是在雨季还是枯季石面流成为最活跃的因素。大石山坡地蓄水容量小,径流容易得到补充同时也易干涸,未来有必要采取工程蓄水措施并对坡地集流面进行管理。

(3) 大石山坡地分为裸岩区、半覆盖区和覆盖区,对此应该分别设置不同的管理目标。裸岩区管控石面流的漏失,半覆盖区通过生态手段增强土壤保水能力岩溶介质通透性,覆盖区采用水文地质手段揭露地下径流,多种措施并举提高坡地径流的利用率。