数值模拟结合同位素技术划分地下水水源地保护区-河海期刊网-河海大学.PDF

第 卷增刊第 期 水 资 源 保 护 ! # 234 5! +6774898:; 13 5 # 年 月 !$$% #$ ’()* *)+,-*.)+ /*,().(0,1 ,; 5 !$$% 数值模拟结合同位素技术划分地下水水源地保护区 # # ! # B # .h/ % +0 %912 %345 %c67 %,89 （ 中国科学院地质与地球物理研究所，北京 ； 北京市水务局，北京 ； 煤炭科学研究总院西安研究院，陕 #5 #$$$!% ! 5 #$$$BA B 5 西 西安 ?#$$= ） 摘要：利用数值方法模拟地下水流入水源井的时间，并根据时间定量划分地下水水源地保护区，是目前国内 外可信度较高的保护区划分方法。计算过程中，采用同位素定年技术辅助判断水文地质参数的取值范围，以 提高模拟精度。以山西省灵石县静升盆地内的供水站为例，采用N,@OP, 和 N,@/’(Q 软件，结合 .O.J 同位素定年数据，对其进行了保护区划分，较好地将数值模拟方法和同位素技术联合应用到水源地保护区划 分中。 关键词：数值模拟；同位素；粒子示踪；地下水；水源地保护区 中图分类号： 文献标识码： 文章编号： （ ） /=# ’ #$$= %BB !$$% +# $$? $= ! ! ! 近几十年来，国内外地下水水源地保护区划分 ! 5 ! 5 ! 渗透系数取值范围的确定 方法迅速发展。在现有的多种方法中，数值模拟法 根据同位素定年数据，计算出沿流线方向上两 能较为准确地刻画复杂水文地质条件，其结果最为 取样点的地下水年龄差，估算两点间的水力梯度和 可靠。但由于水文地质参数的不确定性，往往会在 有效孔隙度，测量两取样点间距，根据达西定律推算 ［］ # 出该区间含水层的渗透系数。参考不同范围内渗透 建模过程中出现较大误差 ，而同位素方法可以有 效地缩小参数取值范围以减小误差。本文以山西省 系数的计算值，进行参数分区。 灵石县静升盆地为例，采用.O.J 定年技术控制关键 将同位素定年技术应用到水文地质参数的估算 参数———渗透系数的变化范围和分区，在用数值法 中，与传统的抽水试验等方法相比，在很大程度上提 模拟出地下水流场的基础上，对盆地内的静升供水 高了参数取值精度，对地下水流场模拟的准确性有 站定量划定了保护区。 极其重要的意义。 ! 5 ! 5 地下水流场的模拟 ! 方法与原理 充分了解研究区地质概况，分析区域地质地貌、 ! 5 ! 划分方法 气象、含水介质、地下水位、边界条件等与地下水流 用数值模拟法划分地下水水源地保护区，并结 动系统相关的因素，掌握包括水源地在内的较大范 合同位素定年数据确定渗透系数的取值范围，大致 围的地下水分布规律，并结合上一步重要参数的分 可分为以下几个步骤： 根据地下水定年数据及其 区结果，建立地下水流动的概念模型、数学模型及数