利用放射性氚估算太行山前地下水更新速率研究.pdfVIP

利用放射性氚估算太行山前地下水更新速率 陈京生n陈宗宇¨ 费宇红2’张兆吉2’张翠云2’ 1)中国地质科学院地下水科学与工程重点开放实验室石家庄050061 2)中国地质科学院水文地质环境地质研究所河北正定050803 摘 要：利用地下水的放射性同位素氚直接估算太行山前浅层地下水的平均更新速率，讨论地下水系统的更新速 率变化。目前条件下，冲积扇顶部浅层地下水的平均更新速率为5—15％，地下水降落漏斗区更新速率卜2％a～，研 究区相应的平均补给强度为224mm／a。这些信息对于确定地下水可持续开采量是十分重要的。 关键词：太行山前地下水系统氚更新速率 引 言 几十年来，由于对华北平原地下水的更新能力认识不足，对地下水的可开采资源量估计过高， 导致开采强度火火超过地下水的更新速度，地下水资源开发产生的一系列问题日益严重。从地下水 资源可持续开发利用的角度出发，人们的目光更加关注水循环交替积极、具有开放性和可恢复性的 浅层地下水系统，强调浅层地下水的改造与合理开发利用、保持浅层地下水的采补均衡【1’2】。这就需 要了解浅层地下水的更新能力，目前对浅层地下水系统的认识基本一致：浅层地下水系统具有可恢 复性【2’3】，一日降水向偏丰发展，尚有恢复的可能[41。但是，地下水系统的更新速率是多大?人们对 此问题尚无明确的结论，从少数研究者早期对华北平原区域地下水的年龄研究结果来看，含水层中 地下水的滞留时间为小于5年至112年筘~7l。但是，经过长期开发，石家庄山前地下水位埋深超 过15米，漏斗中心水位埋深超过40m，这些变化使含水层系统的更新性发生变化，因此这些年龄 数据不能完全代表目前地下水系统的实际情况。因此，必须回答地下水系统的更新速率是多少，并 据此确定可持续的开采量。一般来说，最常用的方法是包气带现场示踪试验、数值模拟等手段评 估地。F水的更新速率。包气带现场示踪试验通常获得的是点上的数据，缺乏区域代表性，受试 验时间的限制，难以了解较长期的平均更新速率；数值模拟在理论上可以计算较长时间的更新 速率，但由下受参数估计和监测资料的限制，往往存在不确定性。尤其在干旱区，由于补给发 生的时空变化人，补给量小，方法很难应用。环境同位素在过去40年里广泛应用于解决与地下水 补给和滞留时间相关的问题【8~11】，特别是放射性同位素(3H和14C)为直接估算干旱和半干早区浅 层地下水系统的更新速率【l2】提供了独特的方法。本文以滹沱河冲积扇为研究区，根据2001年浅层地 一F水放射性同位素氚的调查数据估算地下水系统的平均更新速率，并就其变化进行探讨，为浅层地 下水的可持续开发管理和环境保护提供有益信息。 1研究区概况 研究区位于太行山中段山前倾斜平原地带的滹沱河冲积扇(图1)，地面高程60—1lOm，地势 西高东低，向东倾斜，地形坡度1．6—2．5％，属温带大陆性季风气候区，多年平均气温13。C，降水 分布不均，多年平均降水量570mm，70％左右集中在7—9月份，区内蒸发作用较强，年均水面蒸发 量为900--1200mm。 第一作者简介：陈京生，1957年生，副研究员，从事水文地质研究，cjsl957917@msn．com 372 图l研究区及采样点位置 本区包气带岩性以砂性土为主，降水入渗系数O．3—0．4。含水层为第四系松散沉积物，按岩性禾I 沉积年代、隔水层分布状况、水动力条件等特征，划分为四个含水组，分别对应于全新统、上更新 统、中更新统、下更新统。从山前向东部，含水层由单层结构向双层及多层结构过渡，垂向连续性好。 泄方式主要是人jr开采和天然蒸发。研究区地下水位埋深一般在10--20m，在地下水位降落漏斗区 降水量大小而波动。天然条件下地下水补给主要为火气降水入渗，其次为河渠渗漏。 2方法原理与样品采集测试 地下水更新速率(Renewalrate)定义为地下水系统中补给水体积与含水层总水体积的比。 GalLaSalle 地下水中的放射性同位素氚通常作为降水补给和在含水层中运移的标记，最近，Le 等¨圳提出了利用氚来计算含水层的平均更新速率，其原理简要介绍如下： 假设含水层是稳定状态，即水的输入等于输出，地下水连续补给且完全混合。则地下水中 的放射性氚含量通过放射性衰变