弱透水层释水过程中水力参数响应规律.doc

弱透水层释水过程中水力参数响应规律 李兆峰1,2，周志芳3，李明远3，周翠英*1,2 （1.中山大学工学院，广东 广州510275；2.中山大学岩土工程和信息技术研究中心，广东 广州510275；3.河海大学地球科学与工程学院，江苏 南京211100） 摘要：为了提高地下水运移模拟、地下水资源评价和地面沉降预测的精度，利用室内试验弱透水层释水固结过程中水力参数的规律。基于相邻含水层降深恒定且初始时刻水流稳定条件下弱透水层释水量的解析解，提出了参数求解的配线法，利用自发研制的试验装置进行试验研究。试验结果表明：弱透水层释水过程中渗透系数和贮水率逐渐减小，固结系数变化不大。试验中土层的渗透系数和贮水率分别减小了52％和59％。弱透水层水力参数恒定不变的假设会对其释水量的计算造成较大误差，取弱透水层固结变形初始阶段的贮水率，计算结果，取固结变形结束阶段的贮水率。所得成果研究成果对地面沉降预测和地下水资源量计算具有一定的应用价值。 地下水资源是人类淡水资源的重要组成部分，主要储存在含水层和弱透水层组成的含水层系统中[1]。弱透水层作为含水层系统的重要组成部分，在地下水资源管理和预测中经常被忽视[2,3]。弱透水层在冲积平原和沉积盆地中分布广泛，且主要由黏土、亚黏土、粉质黏土等细粒沉积物组成，具有低渗透性和高储水性的重要性质[4-7]。弱透水层的渗透系数一般小于10–8 m/s，比含水层小几个数量级[8-10]，对地下水环境保护具有重要作用[11]。当抽水含水层水位下降时，弱透水层内部水力减小，有效应力增加，弱透水层固结变形并释水到抽水含水层[12,13]。由于弱透水层的贮水率比含水层大1—2个数量级，弱透水层释水量非常大，弱透水层释水量的计算对于地下水资源评价非常重要[1,2]。弱透水层参数是释水量计算的重要影响因素，因此弱透水层水力参数的确定是含水层系统水流运移模拟、地下水资源评价和地面沉降预测的关键。 弱透水层的水力参数包括渗透系数和贮水率，Field[14]指出固结系数的不确定性是传统固结理论计算变形速率的极限性。很多水文地质学者对弱透水层水文地质参数进行了大量研究[15,16]。余闯等[17]推导出正常固结和超固结状态下固结系数与有效应力之间的表达式；张明等[18]利用统计分析得到固结系数与有效应力的拟合关系；上述研究主要对软土的固结系数进行了分析。叶淑君等[19]利用图解法对上海含水层系统中的弱透水层参数进行了研究。Zhou等[20]人提出了配线法求解弱透水层的渗透系数和贮水率，并通过室内试验验证了方法的可靠性。弱透水层释水过程中水力参数变化规律的研究尚未发现，特别是弱透水层的贮水率。本文在固结理论的基础上推导了弱透水层定降深条件下释水量的解析解，利用该解析解提出了相应条件下求解弱透水层参数的配线法，并通过室内试验对弱透水层的释水过程中水力参数的响应规律进行研究。1原理 Terzaghi一维固结理论在解决软土地基变形控制和预测中发挥着重要作用，至今仍被广泛应用于计算各种荷载条件下土体的固结问题2016YFC0402803） 第一作者：李兆峰（1987-），男，山东安丘人，助理研究员，博士，主要从事地质资源与地质工程。E-mail: Lizhfzx@ 通讯作者：周(a)弱透水层是均质的，且渗透系数和贮水率不随时间变化；(b)弱透水层水平侧向无限延伸；(c)弱透水层始终是饱和状态；(d)弱透水层中的水流为垂向一维流，且服从达西定律。通常情况下弱透水层的渗透系数比含水层小两个数量级以上，因此，弱透水层中的水流可近似为一维流动。坐标轴原点位于弱透水层的上表面，厚度为，坐标向下为正。 图1含水层系统概念模型 Fig.1 Conceptual model of aquifer system 弱透水层的水压力用降深(s)代替，Terzaghi一维固结理论[21]建立数学模型 (1) 其中，是土样的固结系数或者水力扩散系数，是弱透水层的渗透系数是时间。初始时刻弱透水层内部水头线性分布，即内部水流是稳定流。 (2) 其中，是抽水含水层初始时刻的降深。假设弱透水层上部含水层水位不变，即上边界降深为0，下部含水层降深，则 (3) (4) 式(1)—(4)所示的偏微分方程，利用分离变量法可以求得其解析解[22]。 (5) 根据Darcy定律，求得弱透水层中的水流速度,弱透水层底面的水流速度为 (6) 上式无量纲化