第9章 同位素水文地球化学方法.doc

第九章 同位素水文地球化学方法 同位素是指核内质子数相同（即核电荷相同、原子序数相同），而中子数不同（即原子量不同）的同一种元素。它们在元素周期表中居于同一位置。 　　同位素水文地球化学可以用来解决水文学、水文地质学、资源和环境学中的许多问题。在水文学中可用同位素来研究冰川的形成和年龄、河水流量、地表水的弥散作用、湖泊动力学、湖泊和水库的渗漏、悬浮泥沙、底砂搬运和沉积物；在土壤学中解决土壤含水量、土壤密度、非饱和带水的运移；水文地质学中的含水层特征、地下水起源、地下水年龄（滞留时间）、地表水与地下水的联系、含水层间越流、地下水的盐化、地热水的成因、地温；在资源开发中可用同位素来研究矿产的成因；在环境学中用同位素来研究污染源和污染途径；在工程中常用同位素来示踪等。 在同位素水文地球化学中区分有稳定同位素和放射性同位素；环境（天然）同位素和人工示踪同位素等。以下将对放射性同位素、环境稳定同位素和人工示踪同位素分别进行介绍。 第一节　同位素测定地下水年龄法 一、测定地下水年龄的同位素数学物理模型 1.　基本概念 　　根据放射性衰变定律，同位素测定年龄的基本方程为： t = (T1/2/ln2) ln(A0/At) (9.1.1) 　　式中，半衰期为T1/2 的放射性原子的衰变使它的放射性活度在时间段t内从A0减少为At。如果已知t=0时刻的初始放射性活度A0 和t时刻样品的放射性活度At ，则我们就可以测定年龄。但事情并非如此简单，一是五十年代热核试验以来，大气降水中的放射性同位素浓度（A0）已不是常量；二是地下水在渗流过程中会发生不同程度的弥散混合。因此，计算地下水年龄时，首先要建立水混合的物理模型，然后根据同位素物质平衡原理建立数学计算模型。用这种方法建立的地下水测年模型称之为测定地下水年龄的同位素数学物理模型。 　　为了建立数学模型，人们引入了地下水年龄分配函数g(t)的概念（Zuber, 1986, 1994）。g(t)又称系统的反应函数或权重函数，其表达式为： g(t) = CI（t）Q / M （9.1.2） 　　式中，CI（t）为因系统入口处t=0时刻的瞬时注入所导致的在实测点所观测到的示踪剂的浓度；Q为系统的流量；M为注入示踪剂的活度或质量。它表示任一给定年龄（t）的水在水体中所占的份额。这样，对于稳定流的地下水流系统，任何示踪剂的输出浓度Cout（t）与输入浓度Cin的关系可用褶积公式表示： Cout（t）= Cin（t - ?）g(t) exp(- ? ?) d? (9.1.3) 　　这里，t 为同位素输出时间系列（年代）；?为同位素传输时间（年龄）；（t - ?）为同位素输入时间系列（年代）； ?为放射性同位素的衰变常数。该式为同位素数学物理模型的基本形式，可用两种方法求解：（1）假定系统的年龄分配函数h（r），用已知的输入和观测到的输出，求出系统年龄分配函数的参量，这种方法又称混合参数法；（2）利用已知输入和系列观测的输出，计算出年龄分配函数。具体解法这里不予详述，读者可参阅有关文献（Zuber，1986；王恒纯，1991；IAEA，1994）。 2.　常见的计算地下水年龄的地下水数学物理模型 2.1.　活塞流模型（The piston flow model） 　　活塞流模型(PFM)是最简单的一种水流动模型，它假定水流在含水层中运动时，象在活塞筒中被活塞推动一样，完全不发生混合（图9.1.1）。这种情况下，系统输出的放射性同位素浓度是输入放射性同位素浓度与系统周转时间(系统周转时间T指的是在排泄量Q的作用下，系统中储存的流动水体积Vm全部被排出所需的时间，即T = Vm / Q)的函数，其年龄分配函数在数学上可用狄拉克?函数描述（图9.1.2）： g(?) = ?(t-?) (9.1.4) 　　将上式代入（9.1.3）式，积分后得： Cout( t )= Cin( t - ? ) exp(- ? ?) (9.1.5) 图9.1.1活塞流模型示意图 图9.1.2活塞流年龄分配函数图解 　　活塞流模型的地下水年龄等于系统周转时间，近似于地下水真实年龄。严格地说，地下水系统并不存在典型的活塞流，但在解决某些实际问题时，可以近似地应用。 例如，研究包气带均质土层中入渗水的垂向渗透和均质承压含水层中地下水的运动等。 2.2. 指数模型（The exponential model） 1958年Eriksson提出了指数模型(EM)，他的假定条件是基于含水层的渗透性随深度增大而减弱，地下水传输时间随深度呈指数分布。后来，Bredenkamp和Vogel（1970）也提出了类似的观点。Vogel（1970）通过数学推导证明，平面上无限延伸，厚度为H，具均匀线状补给的二维稳态流，