一种用于快速确定非饱和土渗透性函数的方法及装置的制造方法
【技术领域】
[0001] 本发明涉及非饱和土渗透系数实验测定领域，特别是一种用于快速确定非饱和土 渗透性函数的方法及装置。
【背景技术】
[0002] 非饱和渗透系数与土体含水量之间的函数关系称为渗透性函数。地球表面覆盖着 的土体大部分处于非饱和状态，随着饱和度的减小，土的非饱和渗透系数减小的幅度可以 产生几个甚至更多数量级的变化。因此，如何确定土的非饱和渗透系数，则成为常见的边坡 工程、基础工程、大坝等岩土工程以及污染物迀移等非饱和渗流分析中的一个关键问题。
[0003] 目前，土的非饱和渗透性函数的获取主要通过经验公式（间接方法）和渗透试验 (直接方法）确定。由于非饱和土渗透性试验对于试验仪器、试验过程及试验周期都要求很 高，实际应用中几乎都采用非饱和土的土水特征曲线进行预测，这种半经验的确定方法并 不具有普适性，常常出现这种半经验的非饱和土渗透性函数被大量滥用，而导致计算分析 结果与实际工程相差甚远。因此，非饱和土的渗透性函数最好通过试验测出。
[0004] 常用的试验方法分为稳态试验方法（即流速不随时间变化）和非稳态试验方法 (即流速随时间变化）。由于非饱和土在高基质吸力作用下的低渗透性，使得稳态试验方法 极其耗时，因此非稳态试验方法更受青睐，其中最常用的为瞬态剖面法，在试验中量测非稳 态水流中土样的张力及体积含水量分布及其随时间的变化，试验过程中水力梯度和流速都 随时间变化，然后利用体积含水量计算流速，由流速和水力梯度之比给出渗透系数。瞬态剖 面法试验过程十分复杂，要控制流量边界条件，要防止试样的任何部位有明显的浸润面或 饱和，时间周期也较长。相对于流量边界条件，水头边界条件要容易控制和稳定得多，因此 研发一种通过控制水头边界条件、能快速进行非饱和土渗透性试验并能获得非饱和土渗透 性函数的实验装置及方法，对于将非饱和土理论应用于实际工程意义重大。
[0005] 通过对土壤水分运动基本控制方程的偏微分形式进行离散，获得非饱和土渗透系 数的数值求解式，根据数值求解式，结合简单定水头一维垂直非饱和入渗试验来快速精确 地确定非饱和土渗透性函数，是本发明的基本思想。
【发明内容】
[0006] 本发明所要解决的技术问题是提供一种用于快速确定非饱和土渗透性函数的方 法及装置，利用前向欧拉法和后向欧拉法通过对土壤水分运动基本方程Richards偏微分 方程进行离散化和反演求解，得到两种非饱和土渗透系数表达式并研发定水头一维垂直非 饱和入渗实验装置；通过定水头一维垂直非饱和入渗试验测得试样中A z处的体积含水量 Θ t及张力ht随时间变化的一系列数据，将实验数据代入所得的两种渗透系数表达式，联合 二种表达式的计算结果可快速精确地获得非饱和土的渗透性函数。
[0007] 为解决上述技术问题，本发明所采用的技术方案如下：
[0008] -种用于快速确定非饱和土渗透性函数的方法，其特征是包括以下步骤：
[0009] 1)利用前向欧拉法和后向欧拉法通过对土壤水分运动基本方程Richards偏微分 方程进行离散化和反演求解，获得了两种非饱和渗透系数表达式：
[0012] 式中：
[0013] Kt表示体积含水量Θ t对应的渗透系数；
[0014] Ks表示饱和渗透系数；
[0015] = 表示基质势h对于体积含水量Θ求导；
[0016] Δ z表示测定体积含水量Θ t随时间t变化测点处距离水土接触面的距离；
[0017] Θ 3表不饱和体积含水量；
[0018] 式（I)含有体积含水量Θ和空间距离Δ z两个变量，式（II)含有体积含水量Θ、 时间A t、及空间距离Δ z三个变量；
[0019] 2)利用非饱和土渗透实验装置，通过一维垂直非饱和入渗试验测得土样测试点 A z处的体积含水量Θ t及张力ht随时间变化的一系列数据，并将所得到的数据代入式（I) 及式（II)，快速确定非饱和土的渗透性函数。
[0020] -种用于快速确定非饱和土渗透性函数的实验装置，包括试样筒体、量测系统、集 水系统和供水系统，其特征是：
[0021] 量测系统包括用于量测体积含水量Θ t的土壤水分传感器和用于测量张力ht的 土壤张力传感器，以及用于数据采集的数据采集装置，土壤水分传感器和土壤张力传感器 测试端埋设于试样筒内部Az处，另一端与数据采集装置连接；
[0022] 试样筒体内设有上过滤层和下过滤层，上过滤层由上至下分别为上过滤板、上钢 丝网、上滤纸，下过滤层由上至下分别为下滤纸、下钢丝网、下过滤板，上过滤层和下过滤层 之间的空间内用于填充土样，在试样筒体上反滤层下方一段距离上开有两个开口，用于设 置土壤水分传感器和土壤张力传感器；
[0023] 供水系统包括供水桶、潜水栗、进水管和溢流管，潜水栗放置在供水桶内，潜水栗 通过进水管与试样筒体连接，进水管一端与潜水栗连接，另一端位于试样筒体内壁的上过 滤层上方，溢流管一端与试样筒体连接，连接位置位于进水管与试样筒体连接位置的下方， 溢流管另一端位于供水桶内；
[0024] 集水系统包括排水管、蓄水瓶和天平，天平放置在集水系统底座上，蓄水瓶放置在 天平上，蓄水瓶用于收集试验流出的积水，排水管一端置于试样筒体内下过滤层的下方，另 一端伸入蓄水瓶内。
[0025] 优选的方案中，所述的溢流管与试样筒的连接位置位于上过滤层上方5cm处，用 于提供5cm的定水头边界条件。
[0026] 优选的方案中，所述的试样筒体下方设有试验柜台，试样筒体放置在试验柜台上 方，集水系统放置在试验柜台下面。
[0027] 优选的方案中，所述的蓄水瓶上设有蓄水瓶盖，排水管穿过蓄水瓶盖伸入蓄水瓶 中，蓄水瓶盖与蓄水瓶、排水管紧密贴合保证蓄水瓶的密封。
[0028] 优选的方案中，所述的试样筒体上方设有顶盖，顶盖能够开启或关闭。
[0029] 优选的方案中，所述的进水管位于试样筒体内上过滤层的上方的端头上设有水喷 头。
[0030] 本发明提供的一种用于快速确定非饱和土渗透系数的方法及装置，通过采用上述 方法及装置结构，具有以下有益效果：
[0031] (1)利用前向欧拉法和后向欧拉法对土壤水分运动基本方程进行离散和反演求 解，获得了两种非饱和土渗透性函数表达式，直接由得出的表达式结合简单一维入渗试验 便可计算得出非饱和土渗透系数，与瞬态剖面法通过实时计算变动的流速及水力梯度获得 渗透系数相比，方法更加简单，函数表达式意义更加清楚；
[0032] (2)只需要知道固定位置Δ z处体积含水量Θ和张力ht随时间t的关系，代入上 述非饱和土渗透性函数表达式便可以求得土的非饱和渗透性函数，将两个计算表达式在较 低体积含水量与较高体积含水量下的计算结果相结合便可以获得较准确的土的非饱和渗 透性函数；
[0033] (3)将供水系统设计成一种恒定水头的方法，利用潜水栗将不锈钢圆桶中的无气 水通过进水管末端的喷头，均匀喷到过滤板上面，通过过滤板将水均匀分布渗入圆形土体 中，随着水位的升高，当水位达到溢水口高度时，多余的水将从出水口溢出，从而到达控制 恒定水头压力的目的，为试验提供稳定的边界条件。
【附图说明】
[0034] 下面结合附图和实施例对本发明作进一步说明：
[0035] 图1为本发明的装置整体结构示意图。
[0036] 图2本发明的简化入渗模型示意图。
[0037] 图中：顶盖1，水喷头2,进水管3,试样筒体4,水5,上过滤板6,上钢丝网7,上滤 纸8, 土壤水分传感器9, 土壤张力传感器10, 土样11，下滤纸12,下钢丝网13,下过滤板14， 开口 15,排水管16,蓄水瓶盖17,蓄水瓶18,积水19,天平20,试验台21，集水系统底座22， 溢流管23,供水桶24,潜水栗25,水土接触面26,数据采集装置27。
【具体实施方式】
[0038] 关于非饱和土渗透系数表达式的推导：
[0039] 1. 土壤水分运动基本控制方程
[0040] 土壤水分运动基本方程的一维垂直入渗形式如下：
[0042] 式中：
[0043] Θ表不体积含水量；
[0044] t表示时间；
[0045] Κ( Θ )表示对应于Θ的渗透系数；
[0046] h表示基质势。
则式（III)可以变化为如下形式：
[0049] 2. 土壤水分运动基本控制方程的离散和反演求解
[0050] 在一维非饱和入渗模型中，设i代表空间位置，j代表时间位置，时间间隔为At， 空间间隔为A z。可分别