基于氚梯度的地下含水层水文循环信息的识别方法技术

本发明专利技术提供的基于氚梯度的地下含水层水文循环信息的识别方法，步骤有将已知氚活度的示踪剂投入到特定地下含水层中，对特定含水层测定氚活度并绘制氚活度序列曲线，划分氚活度时间变化特征，理论氚梯度随衰变时间变化的曲线簇，计算地下含水层的参考氚梯度并标入曲线簇得到对应的梯度曲线，再将计算的实际氚梯度在梯度曲线上标出并获得该含水层的平均滞留时间。该识别方法有完善的物理背景支撑，操作路线清晰明确，能够定量识别地下含水层的多种水文循环滞时及速率信息，分析结果可信度较高，还可以预测未来地下含水层的氚活度变化情况。

【技术实现步骤摘要】
基于氚梯度的地下含水层水文循环信息的识别方法
本专利技术属于水文地质领域，具体地说是涉及基于氚梯度来确定地下含水层的水文循环特征及相应参数信息的技术方法。
技术介绍
30年来，反射性同位素地球化学方法在地下水水文学领域的研究中应用广泛。测定并分析地下含水层的氚活度能够有效地帮助理解地下水文循环。地下水中的氚(3H)成分主要来源于1950-1970之间的人工核试验，大气降水持续不断的补给地下水从而注入了各个年代的氚，因此，地下水氚活度可以被认为是地下水补给的天然示踪剂。尤其地，氚同位素化学性质非常的稳定，衰变周期为12.43年。目前，利用氚活度获取地下水的水文循环信息的方法出现了一些问题，一是地下水初始的氚活度难以确定，二是水源的混合影响着氚活度信号。现在地下水的定年计算方法主要有活塞模型(PFM)，全混模型(EM)及线性模型(LM)；在PFM中，基本假设是不同流线都有着同一滞留时间，忽视了分子的扩散运动和水动力分布；在EM中不同流线都有着同一滞留时间，不考虑流线的混合；在LM中仅将流线的滞留时间作线性增加考虑。当测试的含水层氚活度较低时，这些模型的计算结果都会存在极大的不确定性，并且计算得到的地下水年龄也是整个含水层的大致平均年龄，对于判断含水层的水文循环特征难以提供直观科学的信息，因此需要进一步研究。
技术实现思路
为了解决现有技术中存在的不足，本专利技术提供一种基于氚梯度的地下含水层水文循环信息的识别方法，将氚活度视为天然输入的同位素示踪剂，提出基于氚梯度的计算方法与氚活度曲线进行结合分析，从而可以实现准确识别所研究地下含水层水文变化的时间响应，为研究地下含水层水文循环提供新的技术方法。为解决上述问题，本专利技术具体采用以下技术方案：基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，包括以下步骤：步骤1：将一定数量的已知氚活度为a0T0的示踪剂投入到地下含水层中；步骤2：对地下含水层一定的时间间隔下，采集地下含水层样品，立即测定对应的氚活度；步骤3：绘制步骤2取得样品的氚活度序列曲线，划分氚活度时间变化特征，包括氚活度滞留开始时间SRT、氚活度滞留中期时间MRT、氚活度峰值出现时刻PT、氚活度最大滞留时间LRT以及氚活度峰值形成历时PFT，推断出地下含水层的多种水文循环滞时及速率信息；步骤4：理论氚梯度计算：选用不同的氚活度初始值aT0计算对应的理论氚梯度Tg，并且设置合适的衰变时间t步长，绘制理论氚梯度Tg随衰变时间t变化的曲线簇即Tg_t曲线簇；步骤5：利用步骤2中首次采样测定的氚活度aT1与最后一次采样测定的氚活度aT2分别计算地下含水层的参考氚梯度Tg_T和实际氚梯度Tg_O：步骤6：地下含水层氚梯度趋势曲线确定：将地下含水层的参考氚梯度Tg_T标入步骤4中的Tg_t曲线簇，查询其对应的Tg_t曲线；步骤7：地下含水层水文循环周期计算：将步骤5求得的实际氚梯度Tg_O代入在步骤6中所得的Tg_t曲线，求出对应的时间ta，ta即为该地下含水层的平均滞留时间。前述的基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，所述步骤1中的示踪剂的氚活度a0T0应当大于30TU，所述示踪剂的数量为1-5升。前述的基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，所述步骤2中的时间间隔为15-30天，持续时间为120-360天。该步骤2中氚活度测试选用低本底液闪仪及其前处理系统(型号：Tri-Carb3170TR/SL)，精度0.1TU；前述的基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，所述步骤3中的氚活度序列曲线是一条通过一系列样点的平滑曲线，样点的纵坐标为地下含水层样品的氚活度，横坐标为地下含水层样品对应的时间，其中时间为0时刻曲线对应的氚活度为(a0T0+aT1)/2。前述的基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，所述步骤3中的氚活度滞留开始时间SRT、氚活度滞留中期时间MRT、氚活度峰值出现时刻PT、氚活度最大滞留时间LRT以及氚活度峰值形成历时PFT通过图解法划分；其中，SRT为峰值开始形成的时间，MRT为SRT到峰值出现时间的中值，MRT＝(PT-SRT)/2，LRT为峰值开始形成到峰值结束的时间，PFT为峰值开始形成到峰值出现的时间，PFT＝LRT-SRT。其中，SRT指示出特定含水层对于外界补给的响应滞时，MRT指示出特定含水层对于外界补给的响应速度，LRT指示出水流运动经过特定含水层的最大历时，PFT指示出特定含水层循环的时间。前述的基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，所述步骤4中的理论氚梯度计算具体如下：选用不同的氚活度初始值aT0计算对应的理论氚梯度Tg：式中，Tg为理论氚梯度，atT为地下含水层的氚活度，t为衰变时间，atT＝aT0×e-t/17.93(2)。前述的基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，所述步骤4中选用不同的氚活度初始值aT0为30TU、40TU、50TU、60TU、70TU，衰变时间t步长为0.5-3年。前述的基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，所述步骤5中计算地下含水层的参考氚梯度Tg_T和实际氚梯度Tg_O的具体操作如下：将首次采样与最后一次采样测定的氚活度分别记为aT1、aT2，分别对应的采样时间为T1、T2，分别计算地下含水层的参考氚梯度Tg_T和实际氚梯度Tg_O：atT2＝aT1×e-t/17.93(5)t＝T2-T1(6)。前述的基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，所述步骤6中将Tg_T标入步骤4中的Tg_t曲线簇的操作是：令y＝Tg_T，x＝T2-T1，在Tg_t曲线簇中，最接近该点(x,y)的曲线为对应Tg_t曲线。通过该Tg_t曲线，可预测未来特定含水层的氚活度变化情况。本专利技术的有益效果：1、技术原理有完善的物理机制，通过人工投放氚活度示踪剂的方式，监测地下含水层氚活度的时程变化情况来反映地下含水层的水文循环过程；2、操作路线清晰明确，综合多源信息来识别氚活度变化，根据氚活度序列曲线及氚梯度趋势能够定量识别地下含水层的多种水文循环滞时及速率信息；3、分析手段全面，分析结果可信度较高，不仅可以准确得到目前地下含水层的水文循环特征，还可以预测未来地下含水层的氚活度变化情况，为未来的相关研究提供科学参考。附图说明图1是基于氚梯度的地下含水层水文循环信息的识别方法的流程示意图；图2是黄河源区的某一地下含水层的氚活度序列曲线图；图3是黄河源区的某一地下含水层理论氚梯度Tg随衰变时间t变化的曲线簇(Tg_t曲线簇)；附图标记如下：1—氚活度滞留开始时间(SRT)，2—氚活度滞留中期时间(MRT)，3—氚活度峰值出现时刻(PT)，4—氚活度最大滞留时间(LRT)，5—理论氚梯度Tg随衰变时间t变化的曲线簇(Tg_t曲线簇)，6—参考氚梯度Tg_T，7—参考氚梯度Tg_T对应的Tg_t曲线，8—实际氚梯度Tg_O。具体实施方式下面结合附图和具体实施例对本专利技术作进一步描述。本实施例是以黄河源区的某一地下含水层为研究对象，人工投放氚活度示踪剂，定时测定地下含水层的氚活度，综合识别地下含水层的水文循环信息。基于氚梯度的地下含水层水本文档来自技高网...

【技术保护点】
1.基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，包括以下步骤：步骤1：将一定数量的已知氚活度为a0T0的示踪剂投入到地下含水层中；步骤2：对地下含水层一定的时间间隔下，采集地下含水层样品，立即测定对应的氚活度；步骤3：绘制步骤2取得样品的氚活度序列曲线，划分氚活度时间变化特征，包括氚活度滞留开始时间SRT、氚活度滞留中期时间MRT、氚活度峰值出现时刻PT、氚活度最大滞留时间LRT以及氚活度峰值形成历时PFT，推断出地下含水层的多种水文循环滞时及速率信息；步骤4：理论氚梯度计算：选用不同的氚活度初始值aT0计算对应的理论氚梯度Tg，并且设置合适的衰变时间t步长，绘制理论氚梯度Tg随衰变时间t变化的曲线簇即Tg_t曲线簇；步骤5：利用步骤2中首次采样测定的氚活度aT1与最后一次采样测定的氚活度aT2分别计算地下含水层的参考氚梯度Tg_T和实际氚梯度Tg_O：步骤6：地下含水层氚梯度趋势曲线确定：将地下含水层的参考氚梯度Tg_T标入步骤4中的Tg_t曲线簇，查询其对应的Tg_t曲线；步骤7：地下含水层水文循环周期计算：将步骤5求得的实际氚梯度Tg_O代入在步骤6中所得的Tg_t曲线，求出对应的时间ta，ta即为该地下含水层的平均滞留时间。...

【技术特征摘要】
1.基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，包括以下步骤：步骤1：将一定数量的已知氚活度为a0T0的示踪剂投入到地下含水层中；步骤2：对地下含水层一定的时间间隔下，采集地下含水层样品，立即测定对应的氚活度；步骤3：绘制步骤2取得样品的氚活度序列曲线，划分氚活度时间变化特征，包括氚活度滞留开始时间SRT、氚活度滞留中期时间MRT、氚活度峰值出现时刻PT、氚活度最大滞留时间LRT以及氚活度峰值形成历时PFT，推断出地下含水层的多种水文循环滞时及速率信息；步骤4：理论氚梯度计算：选用不同的氚活度初始值aT0计算对应的理论氚梯度Tg，并且设置合适的衰变时间t步长，绘制理论氚梯度Tg随衰变时间t变化的曲线簇即Tg_t曲线簇；步骤5：利用步骤2中首次采样测定的氚活度aT1与最后一次采样测定的氚活度aT2分别计算地下含水层的参考氚梯度Tg_T和实际氚梯度Tg_O：步骤6：地下含水层氚梯度趋势曲线确定：将地下含水层的参考氚梯度Tg_T标入步骤4中的Tg_t曲线簇，查询其对应的Tg_t曲线；步骤7：地下含水层水文循环周期计算：将步骤5求得的实际氚梯度Tg_O代入在步骤6中所得的Tg_t曲线，求出对应的时间ta，ta即为该地下含水层的平均滞留时间。2.根据权利要求1所述的基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，所述步骤1中的示踪剂的氚活度a0T0大于30TU，所述示踪剂的数量为1-5升。3.根据权利要求2所述的基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，所述步骤2中的时间间隔为15-30天，持续时间为120-360天。4.根据权利要求3所述的基于氚梯度的地下含水层水文循环信息的识别方法，其特征在于，所述步骤3中的氚活度序列曲线是一条通过一系列样点的平滑曲线，样点的纵坐标为地下含水层样品的氚活度，横坐标为地下含水层样品对应的时间，其中时间为0时刻...

【专利技术属性】
技术研发人员：陈静，钱睿智，杭庆丰，万程炜，衣鹏，龚梦，熊岭，
申请(专利权)人：江苏省水文水资源勘测局，
类型：发明
国别省市：江苏,32