一种寒区化冻期土壤内源污染释放汇流过程监测方法
【技术领域】
[0001] 本发明属于寒区水资源污染监测与治理技术领域，特别涉及一种寒区化冻期土壤 内源污染释放汇流过程监测方法。
【背景技术】
[0002] 水资源由于其形成特点，通常其污染物来源于一个区域范围，而对于区域的污染 物的数量和过程的测定，则远远不像工业污染物排放那样，在排污口监测流量和浓度过程 即可。在区域上，污染物在产生机理、径流条件以及进入水资源承载体过程中都存在着很大 的变异性和不确定性。而在寒区，由于寒区特殊的环境，冻结过程中，污染物会在土壤中聚 集，形成内源，而如何在区域上实现内源释放过程监测，则涉及到三个方面的问题：⑴水 体的流动是内源污染释放的动力学条件，在区域范围内，不同下垫面影响水流运动的关键 因素一径流条件表现出显著的差异，对于内源释放过程需要首先考虑这种差异；(2)由于 化冻期土壤中的内源污染是低浓度的融化水体和高浓度的未冻结水体之间的部分掺混，传 统的对流弥散关系式不适用污染物对于水文过程的驱动响应，当流域存在多种土壤的情况 下，则由于不同土壤融化曲线的差异，内源污染释放过程会表现出显著的差异，而如何实现 内源污染过程的有效监测则面临着重大挑战，此外（3)在多孔介质中，特别是冰水共同存 在土壤中，如何直接进行污染物迀移的实时监测，则更是远未解决的问题。
[0003] 寒区冻土融化期内源污染释放过程监测与实时监控体系的建立离不开以下几个 方面技术的进步与支撑：
[0004] (1) 土壤、水文与污染物自动监测技术的集成，单一的水文自动观测和污染物自动 站均无法获取较准确可信的内源示范数据，而且在化冻期，由于区域包括不同的下垫面类 型和土壤类型，需要将土壤、水文和污染指数等数据统一实现数据集成；
[0005] (2)海量数据的存储和连续监测数据的实时传输技术也对寒区冻土融化期的水文 及面源污染物通量实时监测起到约束作用；
[0006] (3)GPS、DGPS与GIS技术的成熟和发展为建立多维的寒区化冻期水资源污染物的 土壤内源释放监测体系提供了可能；
[0007] (4)流域产汇污机理研究和模型开发是实现全域河流污染通量监测的重要依据， 也为建立河流污染通量监测体系大大节约了人力、物力和财力成本。
[0008] 由于以上四个方面的技术需求，在寒区冻土融化期内源污染释放过程监测中需要 将污染物数据通过算法模型、地理模型的集成，实现监测数据的精确计算和数据解析。目前 数据方法侧重于"汇"的监测，也就是污染物进行水体后对水体产生影响的监测，并没有针 对寒区内源污染释放，即"源"过程的系统监测方法，以及系统的监测数据解析方法。
[0009] 如上所述，一套完整的寒区冻土化冻期水资源污染物土壤内源释放监测体系还没 有建立起来，无法为解决日趋严重的寒区冻土融化期水资源污染问题提供科学支持。
【发明内容】
[0010] 为了克服现有技术的不足，本发明提供了一种寒区化冻期土壤内源污染释放汇流 过程监测方法，以获得冻土化冻期土壤内源释放规律，分析化冻过程中土壤的中的内源释 放过程和均衡过程及其变化规律，为实现寒区水资源污染物"源释放"的远程数字化监测提 供支持。
[0011] 本发明所提供的技术方案具体如下：
[0012] 一种寒区化冻期土壤内源污染释放汇流过程监测方法，包括以下步骤：
[0013] (1)利用流域30m精度的数字高程地图，进行均衡区划分，根据土地利用情况，进 行下垫面类型识别；
[0014] (2)分析化冻期区域尺度上的水文特性和污染物析出规律，确定区域尺度的土壤 污染物释放源强；
[0015] (3)根据质量平衡法分析各均衡区化冻期水量及污染物的均衡过程；
[0016] (4)根据数字高程地图获取、调整均衡区，并回转执行步骤（1)，直至所需均衡区 的水量和污染物质量均衡分析均执行完毕；
[0017] (5)寒区化冻期水资源污染物释放源强数据的计算、显示、数据输出及预警；所述 的污染物至少包括NH 4+、NO3。
[0018] 所述的确定区域尺度的土壤污染物释放源强的方法为：
[0019] 方法（一）：
[0020] ①通过式（1)和式⑵计算冻土中水流运动的通量qtot:
[0022] 其中，qtot为冻土中水流运动的通量，K(h)为水力传导度，h为土壤基质势，T为温 度，z为土壤深度，σ为水动力粘滞系数；
[0023] 冻土中土壤基质势与温度的关系表述为克劳修斯-克拉佩龙方程：
[0025] 式⑵中，Lf为潜热，T为温度，T η为土壤水的结冰温度，g为重力加速度；
[0026] ②现场测定污染物释放动力学系数k，通过式（3)计算化冻期土壤污染物释放浓 度c :
[0027] c = c〇e kt (3)
[0028] 式（3)中，c为化冻期土壤污染物释源强浓度，c。为化冻期土壤污染物本底浓度， t为污染释放时间，k为污染物释放动力学系数；
[0029] ③计算污染物的释放量
[0030] 污染物的释放量=qtotc (4)
[0031] ④计算源强
[0032] 源强=污染的释放量/释放时间 （5);
[0033] 或方法（二）：
[0034] 根据质量平衡法，在均衡区内，不考虑子流域河道水体中储量变化导致的通量变 化，均衡区出口位置的污染物通量为：
[0035] M01=Mp^Mnpi-AM1 (6)
[0036] 其中，Mcil为均衡区出口位置控制断面污染物通量，Mnpi为点源污染入河通量，M pi为 冻土源强污染入河通量，AM1S河道中污染物的衰减通量。
[0037] 所述的方法（一）采用多参数水质分析仪现场测定污染物浓度。
[0038] 所述的方法（二）采用多普勒流速仪测定均衡区控制断面流速。直接得到了寒区 化冻期水资源污染指数。
[0039] 所述的寒区化冻期水资源污染物释放源强数据的显示采用折线图方式，选取观测 均衡区的时间作为横坐标，污染物通量为纵坐标，绘制出指定河流指定区域冻结过程中的 污染指数的变化趋势，选取一年的12月份至下一年4月份作为横坐标，区域的污染物浓度 为纵坐标，直观反映区域的污染情况。
[0040] 所述的寒区化冻期水资源污染物释放源强数据的显示，其用户视图与标准地理经 炜度坐标相对应。河流坐标定位精确可靠。
[0041] 与现有技术相比，本发明具有以下优点和有益效果：
[0042] (1)采用30m精度的DEM (Digital Elevation Model)数字高程地图数据，以对小 流域尺度中子流域进行精确划分；
[0043] (2)采用多普勒流速仪型流速仪进行均衡区控制断面流速，采用多参数水质分析 仪现场测定污染物浓度，采用解析反演法根据同期进行的惰性离子水团迀移过程监测数据 确定污染物纵向弥散系数，实现了现场全部参数的实时测定，根据数学模型计算上、下游断 面污染物浓度，河道污染物降解系数；
[0044] (3)绘制指定河流指定区域冻结过程中的污染指数的变化趋势，采用的是折线图 方式，选取一年的12月份至下一年4月份作为横坐标，区域的污染物浓度为纵坐标，直观反 映区域的污染情况；
[0045] (4)寒区化冻期各均衡区的土壤污染物释放源强的数据显示采用的是折线图方 式，选取观测汇流区的时间作为横坐标，污染物释放源强为纵坐标，直观反映土壤污染物释 放的源强的情况；
[0046] (5)用户视图与标准地理经炜度坐标对应，河流坐标定位精确可靠。
【具体实施方式】
[0047] 为了便于本领域普通技术人员理解和实施本发明，下面结合实施例对本发明作进 一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本发明，并不用于限 定本发明。
[0048] 一种寒区化冻期土壤内源污染释放汇流过程监测方法，包括以下步骤：
[0049] (1)利用流域30m精度的数字高程地图，进行均衡区划分，根据土地利用情况，进 行下垫面类型识别；
[0050] (2)分析化冻期区域尺度上的水文特性和污染物析出规律，确定区域尺度的土壤 污染物释放源强；
[0051] (3)根据质量平衡法分析各均衡区化冻期水量及污染物的均衡过程；
[0052] (4)根据数字高程地图获取、调整均衡区