格尔木河河床沉积物渗透系数变异性研究-水文学及水资源专业论文.docx

万方数据 万方数据 本论文的研究得到了长安大学承担的以下三个项目的资助： （1）国家自然科学基金项目：《旱区地表—地下水系统界面动力 学与水循环研究》（项目编号； （2）中国地质调查局地调项目：西北大型盆地水文地质调查综合 研究工作项目专题《西北大型盆地关键水文地质参数试验研究》 （项目编号：12120113104200）； （3）中国地质调查局地调项目：《柴达木盆地循环经济试验区重 点地区（格尔木、乌兰）地下水勘查》（项目编号：1212011220975） 子课题《格尔木河流域与乌兰盆地地下水循环和数值模拟》. I 摘 要 河床沉积物的渗透系数 K 是关键的水文地质参数之一，它直接影响着河水与周围地 下水之间的转化能力，河流对含水层的入渗补给量，含水层向河流的排泄量以及傍河水 源地的开采能力。渗透系数代表孔隙介质传送水量的能力，准确的估算河床的渗透系数 对于河流一地下水相互关系及其流域的生态系统分析极其重要，同时有助于了解地下水 与河流之间的水量交换以及水质运移的程度,并能解决一系列相关的水文地质问题。本 文采用水头下降竖管试验法对格尔木河流域河床沉积物的渗透系数值进行了测定分析， 采用传统统计学方法及地统计学方法对实测数据的正态性分布和变异性特征做了深入 的分析。研究主要取得了以下几方面的研究成果： 一、河床沉积物渗透系数随深度的空间变异特征 60河床沉积物渗透系数 K 值随深度递减。表层 0~20cm 沉积物 K 值最大达 52.25m/d， 30cm 处最小，仅 0.69m/d，40~50cm 较稳定。河床沉积物渗透系数 K 值随沉积物厚度的 增加而逐渐减小。10~20cm 厚的柱样 K 值最大，30~50cm 厚的柱样 K 值变幅较小，较 稳定。河床 K 值随深度呈现对数正态分布，变异系数值为 1.0391，为强变异性；河床 沉积物渗透系数 K 值随沉积物厚度呈现对数正态分布，变异系数值为 1.0131，为强变 异性。河床 K 值随深度的变程值为 27cm，空间相关性为 9.94%，表现出强烈的空间相 关性。河床沉积物渗透系数 K 值随深度的变化主要受到沉积物岩性的影响。沉积物粒径 随深度递减，K 值与 d60 相关性最强，拟合度达 0.86，K=10.47d 2。 二、典型剖面河床沉积物渗透系数的空间变异特性 60 河床沉积物渗透系数值由河床中心向两岸递减，横剖面的变异程度大于纵剖面。 1. 上游冲洪积扇区，河床沉积物渗透系数值沿河床横纵剖面均发生明显变异。纵剖 面 K 的平均值为 68.62m/d，横剖面均值为 18.78m/d，总体在 0.75~233m/d 之间变化。受 岩性条件控制，上游冲洪积扇区颗粒粒径大小悬殊，分选性及均一性很差，河床沉积物 渗透系数值沿横纵剖面均服从对数正态分布，变异系数值均1，强变异性。变异函数计 算结果表明该区河床纵剖面 K 值的变程为 4m，空间相关性为 0.8%，强烈相关；横剖面 为 15m，空间相关性为 29%，中等相关。 2. 下游细土平原区，河床沉积物渗透系数值沿河床横纵剖面变化较明显，而沿纵剖 面几乎不变，纵剖面 K 均值 19.12m/d，横剖面 K 均值 13.93m/d，总体在 0.65~47.37m/d 之间变化。该区沉积物颗粒均一性较好，K 值差距不是很大，服从正态分布，变异系数 I 值在 0.1~1 之间，中等变异性。变异函数计算结果表明该区河床纵剖面 K 值的变程为 22.25m，空间相关性为 50%，中等相关；横剖面为 27m，空间相关性为 5.7%，强烈相 关。 3. 典型剖面河床沉积物 K 值主要受到河床形态、流速，潜流带的水文过程及沉积 物岩性的影响。研究结果表明冲洪积扇区由于沉积物颗粒分选差，沉积物 K 值与特征粒 1.1196 径 d30 之间相关性最强，K 的关系式为：K=4.1572d30 。细土平原区沉积物均一性较 0好，K 值与特性粒径 d102 相关性最强，K 的关系式为：K=2021.8*d1 2。当均匀系数大于 5 时，K 随 Cu 的增大而减小，小于 5 时，二者无关。当曲率系数大于 2 时，K 随 Cc 的 增大而减小，小于 2 时，二者无关。 三、流域尺度河床沉积物渗透系数的空间变异特征 0 流域尺度下河床沉积物渗透系数值沿流程总体呈现下降趋势，从出山口至入湖口河 段 K 的平均值为 42.587m/d，在 8.390m/d 至 98.570m/d 之间变化，服从对数正态分布， 变异系数值为 0.743，中等变异。流域尺度下河床沉积物渗透系数值的空间最大相关长 度为 15km 左右，空间相关性为 8.18%25%，空间相关性非常强。最优拟合模型为 LnK 下的球状模型块金值 C0 大于 0，说明该区