冻土水盐时空变异性条件模拟研究.pdfVIP

农1V31U3土环境保护 冻土水盐时空变异性条件模拟的研究+ 徐英1 陈亚新1 史海滨1 魏占民1 李延林2 张艺强2 (1|内蒙古农业大学水利与土木建筑工程学院，呼和浩特010018： 2．沙壕渠试验站，内蒙古临河015000) 摘要： 本文用地质统计学的变差函数讨论了试验区冻融土中不同时期非饱和带土壤垂 向水分和盐分的时空变异性；运用条件模拟理论用较少的采样点再现其水盐的波动过程、 并同普通Kriging最优插值法进行了对比分析和评价：同时研究了影响条件模拟Monte— Carlo实验计算精度和效率的两个重要因素一随机数种子和随机数序列长度的取值规律 得到了比较平稳的取值范围。 关t词： 条件模拟冻土土壤水一盐时空变异性 中国北方农业冻土属季节性冻土，有“秋汇”、“冬灌”储水灌溉经验，其土壤冻融 过程明显地影响着水分和盐分的再分配，而冻土中水盐的特殊运动规律和再分配特性是土 壤盐渍化发生、发展和演变的重要因素，也是干旱地区土壤水资源转换、储水保墒的主要 机制，研究冻土中水盐迁移规律具有重要意义。对于冻融期非饱和带土壤水和地下水运动 属于垂向蒸发一人渗补给型的灌区，如河套大型灌区，把握其垂向土壤冻融过程又是研究 该地区冻土水盐迁移规律的关键。 地质统计学条件模拟理论由于不仅能保持空间变量的白相关性，而且能使实测点处 的模拟值等于实测值；可用较少的采样点模拟许多未知信息及时空分布，从而大幅度减轻 野外工作量，回避常规数值方法的求参问题。本文试探采用这种方法研究冻土中垂向水盐 的时空变异性，以指导节水灌溉田间采样与水盐监测工作。 1．模拟的基本原理及应用方法 假设已有ZG)在n个点的实测数据z。O。)，。=1．2，3，…n，条件模拟要求将模拟条件 化到实测数据，即实测点处模拟值等于实测值。 Kriging(简称OK) 对于非实测值的任意一点石处的真实值z0)可表示为Ordinary 估计值z：G)与其误差RG)之和，即： (1) z0)=Z：O)+kO)一Z：O)J=z：O)+RG) 如果用与RG)同构且独立的非条件模拟OK误差k，0)一z：O)』来代替Ro)，即可得 条件模拟的理论计算方程： ’本文获国家自然科学基金项目资助 农业高效用水与水土环境保护 Z。O)=z：G。)+【z，O)一z：0)J (2) 其中z：o)：至^。zo。)，z：G)；￡a。z，厶。)，这样(2)式又可记作 4‘I 4。】 (3) o。)一Z，o。)】 z，，G)；z，G)+￡土。【z d。I (3)式是常用的条件模拟计算方程，用它可得任意多个条件模拟现实。 由式(3)知条件模拟关键是求非条件模拟项，变差函数模型为球状模型的一维非条 件模拟可以作为具有权函数，矗)的一种滑动平均来实现。实践中将具有Dirac协方差测 度的随机测度r(出、用沿着一条直线位于等间距上的一系列随机变量进行离散逼近的。这 些随机变量具有独立同分布，且E仃．)：0，E0：)：，：。将正的现实l按等间隔b分 i+k处用随机 配到直线的各离散点上。令r。，…，r。，…，t。表示在直线上点i-k，…，i…t 数种子Rd产生的诸现实．然后在每个点陂b求出指定权函数，函)加权的滑动平均数z。。 实践中只须求有限多项离散和，球状模型取奇数(即2R+1)个基本值t。(k取整数， 从．R到+R，R为自然数)，z；可表示为： 出 zi=]t…毋出)