基于地下水性质空间变异特征的优化地统计模式研究-农业工程学报.DOCVIP

黄河三角洲地区浅层地下姚荣江杨劲松※（中国科学院南京土壤研究所，南京 210008） 摘 要：为研究黄河三角洲地区地下水作用条件下耕层土壤积盐规律，区土壤；浅层地下水；土壤积盐；S156.4 文献标识码：A 0 引言 盐渍化是黄河三角洲地区的主要障碍因子，区独特的气候、水文水资源以及土壤等决定了在空间上并不完全独立，而是存在空间上的有关这方面的研究已有很多详细的(1)五十年代，主要是应用数理统计法。通过大面积土壤和水文地质调查，统计地下水埋深(矿化度)与表土盐渍化现象的关系。该方法不涉及溶质迁移的动力机制，只能作自然现象的宏观描述。(2)六十年代我国用土壤强烈毛管水上升高度加主要根系活动层深度，计算安全深度，从构成安全深度的微观揭示和量化方面较传统理论有较大进步[9]。(3)七十至八十年代，随着地中渗透仪方法的广泛应用，开始将大气－土壤－地下水作为统一体，使不同地下水埋深、不同土质、有无作物覆盖的潜水蒸发影响的研究得以系列化、标准化，促进了水均衡研究，但大多数地中渗透仪内无专门的水、盐传感设备，不能获得其中含水量、盐分含量的时空分布。(4)八十至九十年代，在土壤水盐运动的能量转换理论和地下水数值模拟技术的推动下，在传统地中渗透仪或大型室内土柱内部增加水分或盐分传感器，并进行近代技术的更新改造，使数学模型的数值解与实验技术得以结合。进而用土壤水盐数据采集系统, 对大气－土壤－地下水统一体的水盐运移进行动态模拟。 研究表明，是研究土壤属性的空间变异或其它农田特征变异定量化的有效。提供一定的理论和依据区土壤1 研究区概况及研究方法 1.1 研究区概况 位于山东省垦利县永安镇东义和，地处北温带大陆性季风型气候带，东临渤海，属典型黄河下游三角洲地区。四季气候变化明显，降水主要集中在7~8月份，占全年降水量的70%，年均蒸降3.22。由于土壤直接发育于沉积物，再加上成陆过程中受海水浸渍侧渗，使得研究区土壤含有较多的可溶性盐，滨海潮盐土是该区最主要的盐土类型。平均矿化度g/L；由于当地缺乏蓄水工程，雨季水量丰富但不能调蓄，非雨季淡水资源紧缺，土壤含盐量呈现“积盐—淋溶”动态变化。 图1 研究区地理位置采样点分布Fig.1 The geographic position of study area and the distribution of sampling sites 研究方法 1.2.1 采样布置 选取一块面积约520 hm2的连续地块(东西长约4.1 km，南北宽约1.1~1.6 km)作为研究区域，研究区域地势相对平坦，其南部边界为一条自西向东的灌溉渠，东南部位毗邻大片的养殖鱼塘。采样地界于东经118°47’~118°50’，北纬37°33’~37°34’之间根据当地土质、植被类型等因素确定采样点的位置数量，共挖取个剖面由于区内地下水埋深较浅，每个都挖到地下水流出各采样点的相对坐标采用差分GPS定位技术确定2004年10月中下旬，棉铃吐絮期样品处理与分析地下水埋深：每个剖面都挖到地下水流出，待其水位稳定以后，用刻度尺测量其埋深。 地下水：取水样置于100 m塑料瓶中，在冰箱中保存。地下水含盐量由地下水中各离子组成的含量计算得到，离子组成含量的测定采用常规分析法测定方法参考《土壤农业化学分析方法》。采集的土样带回实验室内自然风干，磨碎、过2 mm筛后备用。采用常规分析法确定土壤各离子组成含量，计算出相应的土壤全盐含量测定方法2 结果与分析 表1 地下水埋深、矿化度及耕层土壤盐分的统计特征值 Table 1 Statistical feature values of groundwater depth, mineralization and soil salinity in plough layer 变量 分布 类型 样本数最小值最大值平均值标准差变异系数K-S p 合理采样数 埋深/ m 127 0.50 2.70 1.79 2.00 0.57 31.87 0.058* 42 矿化度/ g·-1 N 120 3.49 70.53 32.43 32.46 14.16 43.67 0.567 77 盐分/ g·kg 注：N表示正态分布；LN表示对数正态分布；*显著性水平在p0.05p0.10；合理采样数N＝λ2α,f(CV/k)2，其中λα, f为t分布特征值α为显著水平，f为自由度)，CV为变异系数，k为相对误差 地下水埋深、矿化度及耕层盐分统计特征值表1。从变幅来看：地下水矿化度与耕层土壤盐分的变幅较大，分别达67.04 g·L-1和34.15 g·kg-1，其最大值与最小值之比分别为20.21和88.56；地下水埋深的变化范围在0.50~2.70 m之间，变幅最小。从平均值来看，耕层土壤盐